Abstract:This article proposes a distributed secondary control scheme that drives a dc microgrid to an equilibrium point where the generators share optimal currents, and their voltages have a weighted average of nominal value. The scheme does not rely on the electric system topology nor its specifications; it guarantees plug-and-play design and functionality of the generators. First, the incremental model of the microgrid system with constant impedance, current, and power devices is shown to admit a port-Hamiltonian (pH) representation, and its passive output is determined. The economic dispatch problem is then solved by the Lagrange multipliers method; the Karush-Kuhn-Tucker conditions and weighted average formation of voltages are then formulated as the control objectives. We propose a control scheme that is based on the Control by Interconnection design philosophy, where the consensus-based controller is viewed as a virtual pH system to be interconnected with the physical one. We prove the regional asymptotic stability of the closed-loop system using Lyapunov and LaSalle theorems. Equilibrium analysis is also conducted based on the concepts of graph theory and economic dispatch. Finally, the effectiveness of the presented scheme for different case studies is validated with a test microgrid system, simulated in both MATLAB/Simulink and OPAL-RT environments.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the uneven current distribution of distributed generators (DGs) and the deviation of voltage from the nominal value in DC Microgrids. Specifically, the paper proposes a distributed secondary control scheme, aiming to drive the DC microgrid to an equilibrium point, at which:
1. **Current sharing among generators**: Each generator can share the optimal current.
2. **The weighted average value of voltage is close to the nominal value**: The voltages of generators form a weighted average value, which is close to the nominal voltage value.
### Background and Motivation
With the development of green technologies, power systems are gradually integrating more renewable energy sources, energy storage systems, and electrical loads. These devices are interfaced with the power grid through power - electronic converters and are managed in the form of microgrids for ease of control and management. In particular, DC Microgrids have received extensive attention in recent years due to their compatibility with DC power sources and the simplification in control and management.
However, although the existing control method such as Droop Control is simple and robust, it cannot ensure the ideal current sharing among generators and voltage formation. Therefore, a new control strategy needs to be designed to solve these problems.
### Main Contributions
1. **Modeling and Representation**:
- The paper first shows that the incremental model of the DC microgrid system considering the dynamics of transmission lines and shunt capacitors can be represented as a Port - Hamiltonian (pH) system and defines its passive output.
- Then, the economic dispatch problem is solved by the Lagrange multiplier method, and the Karush - Kuhn - Tucker (KKT) conditions and the formation of the weighted average value of voltage are taken as control targets.
2. **Control Strategy**:
- A control scheme based on the Control by Interconnection (CbI) design philosophy is proposed, in which the consensus controller is regarded as a virtual pH system and is connected to the physical system.
- The regional asymptotic stability of the closed - loop system is proved, and the analysis is carried out using Lyapunov and LaSalle theorems.
3. **Equilibrium Analysis**:
- Based on the concepts of graph theory and economic dispatch, an equilibrium analysis is carried out to verify the effectiveness in different case studies.
### Mathematical Model
- **System Dynamics**:
\[
\begin{aligned}
L_{Ge_i} \dot{I}_{Ge_i} &= V_i - \sum_k b_{Ge_ki} V_{Ne_k} - R_{Ge_i} I_{Ge_i}, \\
L_{Ee_j} \dot{I}_{Ee_j} &= -\sum_k b_{kj} V_{Ne_k} - R_{Ee_j} I_{Ee_j}, \\
C_{Ne_k} \dot{V}_{Ne_k} &= \sum_j b_{kj} I_{Ee_j} + \sum_i b_{Ge_ki} I_{Ge_i} - I_L_k, \\
I_L_k &= G_{cte_k} V_{Ne_k} + I_{cte_k} + \frac{P_{cte_k}}{V_{Ne_k}}, \\
V_i &= V_{ref_i} = V_{nom} - R_{D_i} I_{Ge_i} + u_i.
\end{aligned}
\]
- **Hamiltonian Function**:
\[
H(x) = 0.5 x^\top Q x,
\]
where \( Q = \begin{bmatrix} L_{Ge}^{-1} & 0 & 0 \\ 0 & L_{Ee}^{-1} & 0 \\ 0 & 0 & C_{Ne}^{-1} \end{bmatrix} \), \( x = \begin{bmatrix} \phi_{Ge} \\ \phi_{Ee} \\ q_{Ne} \end{bmatrix} \)