Journal of Modern Power Systems and Clean Energy

ISSN 2196-5625 CN 32-1884/TK

网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

A Unified Control of Super-capacitor System Based on Bi-directional DC-DC Converter for Power Smoothing in DC Microgrid  PDF

  • Xialin Li
  • Pengfei Li
  • Leijiao Ge
  • Xunyang Wang
  • Zhiwang Li
  • Lin Zhu
  • Li Guo
  • Chengshan Wang
the Key Laboratory of Smart Grid of Ministry of Education, Tianjin University, Tianjin 300072, China; the School of Information and Electronic Engineering, Hebei University of Engineering, Handan 056038, China; the Binhai Power Supply Bureau of State Grid Tianjin Electric Power Company, Tianjin 300480, China

Updated:2023-05-22

DOI:10.35833/MPCE.2021.000549

  • Full Text
  • Figs & Tabs
  • References
  • Authors
  • About
CITE
OUTLINE

Abstract

To improve the equivalent inertia of DC microgrids (DCMGs), a unified control is proposed for the first time for a bi-directional DC-DC converter based super-capacitor (SC) system, whereby power smoothing and SC terminal voltage regulation can be achieved in a DCMG simultaneously. The proposed control displays good plug-and-play features using only local measurements. For quantitative analysis and effective design of the critical parameter of unified control, two indices, equivalent power supporting time and inertia contributed by the unified controlled SC system, are introduced firstly. Then, with a simple but effective reduced-order model of a DCMG, analytical solutions are obtained for the two indices. In addition, a systematic design method is presented for the proposed unified control. Finally, to verify the proposed unified control, a switching model is developed for a typical DCMG in PSCAD/EMTDC, and theoretical analyses are conducted for different operating conditions.

I. Introduction

IN remote areas and islands, AC/DC microgrids (MGs) integrated with conventional diesel generator set (DGS), energy storage systems (ESSs), and highly penetrated renewable power generations (RPGs) are emerging as promising solutions to power supply problems. However, when subjected to uncertain power disturbances, dynamic stability and power smoothing issues are the main challenges presented by the limited dynamic response and power creep rate of the DGS [

1]-[4]. In weak-grid-connected DCMG, no appropriate power smoothing methods are found to deal with the instability of DC-AC converters [5]. Therefore, it is critical to smooth the transient power and enhance the equivalent inertia of DCMGs.

In low-voltage DCMGs, system inertia can be improved by using the energy stored in DC bus capacitors via adaptive virtual inertia control with variable droop gain [

6] or DC virtual synchronous control [7], [8]. As the energy stored in DC bus capacitors is limited, the control methods in [6]-[8] cannot achieve the desired power smoothing performance when large power disturbances occur. Thus, power-type ESSs such as super-capacitors (SCs) can be used to improve the immunity to disturbances and to enhance the transient stability of DCMGs. To this end, this paper focuses on how to use SC systems to improve the equivalent inertia of DCMGs.

The main publications on control methods can be categorized into direct methods based on power perturbation by divider filtering [

9]-[11] and indirect methods based on the equivalent output impedance of hybrid ESSs with different dynamics [12]-[15]. In direct methods, the high-frequency and low-frequency components of power disturbances can be obtained by filtering the measured power disturbance [9], the output current reference from a centralized DC voltage controller [10], or the locally measured DC bus voltage [11]. These components can be set as the power or current references of the SC and DC voltage control units, respectively. In indirect methods, the equivalent output impedances of the SC and DC voltage control units are designed with different dynamics to realize multiscale decoupling of power smoothing. In [12]-[16], the methods such as virtual capacitors, integral droop control, and virtual resistive-capacitive (RC) droop control have been proposed to make the equivalent output impedance of the SC system capacitive, thereby realizing an effective power smoothing performance in case of a large power disturbance.

To the best of our knowledge, the methods in [

9]-[16] still have the following disadvantages.

1) Poor plug-and-play feature. Additional sensors, high bandwidth, reliable communication [

9], or a centralized controller [10] are necessary in direct methods, which share the risk of a single point of failure. Additional communication can also affect the plug-and-play features of an SC system. In [12]-[16], although both the DC voltage control unit and the SC system adopt droop control, the power smoothing performance is strongly dependent on the equivalent output impedance. Therefore, the plug-and-play capacity of an SC system is limited.

2) Ignoring the impact of SC terminal voltage regulation on the power smoothing performance. The power smoothing performance of the SC system is mainly related to the terminal voltage. Most existing studies suggest that the dynamics of terminal voltage regulation and power smoothing control of SC systems be decoupled. When the power smoothing transient is over, the steady terminal voltage of the SC system is designed to recover to its rated value [

14], [15] or to within a certain range such as 30%-70% of the rated value [16], or 50% of the rated value [17]. If restored to the rated value, the SC system loses the ability to balance disturbances in the case of a power increase of RPGs. Although the SC system can deal with uncertain power disturbances when its terminal voltage is recovered to 50% or to within a certain range, there is little work on determining whether these voltage settings are optimal for power smoothing control over a wide operating range of the DCMG.

To overcome these issues, the main contributions of this paper can be summarized as follows.

1) First, a unified control has been proposed for a bidirectional DC-DC converter based SC system, which can improve the equivalent inertia of the DCMG, smooth the active power of grid-forming units, and control the SC terminal voltage simultaneously. It should be noted that the proposed control has good plug-and-play features by only using local measurements.

2) Second, to represent and quantify the power smoothing capability more clearly, two indices, i.e., equivalent power supporting time and equivalent inertia contributed by the unified controlled SC system, are introduced. The analytical solutions for these two indices can be obtained using a simple but effective reduced-order model of the DCMG.

3) Third, based on the two analytical indices, a systematic design method is presented for the critical parameters of the proposed unified control, which is further verified by detailed simulations based on a typical DCMG switching model constructed in PSCAD/EMTDC.

The remainder of this paper is organized as follows. Section II presents a typical DCMG structure. The unified control for the SC system is proposed in Section III. An analytical solution of equivalent power support time Ts is presented in Section IV by a simple but effective reduced-order model of DCMG. The systematic analysis, design of the proposed unified control, and simulation verification are provided in Section V. Finally, Section VI concludes the paper.

II. Typical DCMG Structure

A typical DCMG structure is shown in Fig. 1, including grid-forming units, grid-following units, and grid-supporting units [

18]. For a grid-connected DCMG, the DC bus voltage is typically controlled by a bi-directional DC-AC converter. In the islanded mode, an ESS such as a battery system with a bi-directional DC-DC converter is generally used to control the DC bus voltage. These units can be considered as grid-forming units of the DCMG. Grid-following units can consist of RPGs such as wind power (WP) and photovoltaics (PV) in the maximum power point track (MPPT) mode, controllable DGs in power dispatching mode, and loads.

Fig. 1  Typical DCMG structure.

In remote areas or islands, DC-AC converters powered by a weak grid or DGS can be considered as the grid-forming units of the DCMG. Since these converters do not have the ability to respond quickly, a large power disturbance from the grid-following units may have a significant impact on the dynamic stability of the DC bus voltage. Thus, the SC system was adopted as a grid-supporting unit to improve the equivalent inertia of the DCMG, suppress power disturbance, and smooth the active power of the grid-forming unit. In this paper, a simple but effective control method is proposed for the SC system, referred to as unified control. udc is the DC bus voltage, and Ps, Pp, and Psc are the power outputs of grid-forming units, grid-following units, and grid-supporting units, respectively.

III. Unified Control for SC System

In this section, the requirements of the power smoothing control in the DCMG are initially introduced. Then, a simple unified control is proposed for the SC system, to improve the equivalent inertia of the DCMG, smooth the active power of the grid-forming unit, and control the SC terminal voltage simultaneously.

A. Requirements of Power Smoothing Control Performance in DCMG

In this paper, the expected power smoothing control performances of the DCMG are illustrated, as shown in Fig. 2(a) and (b). Considering a large power disturbance ΔPp (e.g., a power variation from RPGs or loads) at time t0, the SC system is expected to compensate for the unbalanced power immediately, with a dynamic response ΔPsc. To this end, the power of the grid-forming unit ΔPs, DC bus voltage udc, and terminal voltage of the SC system usc smoothly reach an acceptable equilibrium point with usc1 and usc2 denoting the terminal voltages of the SC system in the pre- and post-disturbance states, respectively. ΔEsmooth is the released energy, and Ts is the equivalent power supporting time.

Fig. 2  Ideal dynamic responses of power and voltage. (a) Ideal power dynamic response. (b) Ideal voltage dynamic response.

In some published studies, the DC bus voltage udc was regulated to a constant reference value. In contrast, the DC voltage droop control strategy of DC-AC converter [

19], [20] is adopted for the grid-forming unit of the DCMG in this paper, as shown in Fig. 3. All values of the variables are per-unit, i.e., normalized by the base values of power and DC voltage (denoted as PB and UdcB); udc,set and udc,pu are the reference value and measurement of the DC voltage, respectively; Ps,pu is the active power of DC-AC injected into DCMG; Rd is the droop gain; kpu,s and kiu,s are the proportional and integral gains of the DC voltage controller, respectively; idref,pu and iqref,pu are the reference values of the currents projected onto the d- and q-axis, respectively; and dabc is the the duty ratio.

Fig. 3  DC voltage droop control of DC-AC converter.

With droop control, the relationship between the power of grid-forming or grid-following units (referred to Ps,pu and Pp,pu, respectively) and DC bus voltage udc,pu in steady state can be expressed as:

udc,pu=udc,set-Ps,puRd=udc,set-Pp,puRd (1)

Hence, by measuring udc,pu, the SC system can obtain the power of grid-forming or grid-following units readily, without any additional sensors or communication, which is critical for the plug-and-play feature. In this paper, udc,set=1, and given that the variation range of the total power of the grid-following units in DCMG is -1Pp,pu1, the operation range of the DC bus voltage will be 1-Rdudc,pu1+Rd.

As depicted by the gray shaded area in Fig. 2(a), the released energy ΔEsmooth can be represented by the power smoothing ability of the SC system, whereby the detailed expression is given by:

ΔEsmooth=t0tΔPscdt=12Cscusc22-usc12 (2)

where Csc is the capacitance value of SC; and usc1 and usc2 are the terminal voltages of the SC system in the pre- and post-disturbance states, respectively.

Regarding the DGS in the isolated DCMGs, which is limited by the response or power climbing rate, a certain duration of time is required by the SC system to support the DCMG. In other emergency scenarios, a certain duration of time is also critical for the power-coordinated control response of the DCMG such as load shedding or power curtailment of RPGs. Thus, to evaluate the power smoothing ability of the SC system more intuitively, the concept of equivalent power supporting time Ts is proposed and defined as:

Ts=ΔEsmooth/ΔPp (3)

The main requirements of power smoothing control can be summarized as follows.

1) Control performance, as shown in Fig. 2(a), can be realized using only local measurements. Moreover, a given set value of Ts can be satisfied via the specific design of the critical control parameters.

2) The existing works in [

14]-[17] do not discuss how to set the steady terminal voltage value of the SC system from the perspective of being more effectively to cope with uncertain power disturbances more effectively. If the SC terminal voltage control can be coupled dynamically with power smoothing control on the same timescale, the design of the SC control system can be simplified.

Motivated by these considerations, a novel unified control of an SC system is proposed for power smoothing applications.

B. Proposed Unified Control of SC System

1) Unified Control Structure

The proposed unified control of the SC system unifies the SC terminal voltage control and power smoothing control, as shown in Fig. 4, in two parts consisting of DC voltage consensus control and base value regulation of SC terminal voltage. It should be noted that only local measurements of the DC voltages and currents are used here, indicating good plug-and-play feature.

The easiest way to smooth the power disturbance and improve the inertia of the DCMG is to connect the SC directly to the DC bus as part of the DC bus capacitor. However, as the terminal voltage level of the SC is usually different from that of the DC bus, the SC is always connected to the DC bus via a bidirectional DC-DC converter. Thus, to make the SC system an equivalent DC bus capacitor, a simple unified control is proposed, as shown in Fig. 4, where ilsc is the inductor current; and ihsc and Chsc are the output current and DC capacitor on the DC bus side, respectively.

Fig. 4  Proposed unified control of SC system with plug-and-play feature.

With this simple consensus control, the per-unit value of the SC terminal voltage can be dynamically tracked to that of the DC bus voltage, such that:

uscUscBudc/UdcB=UscBudc,pu (4)

where usc and UscB are the terminal voltage and base value of the SC, respectively; and udc and UdcB are the actual and base values of the DC bus voltage, respectively. To improve the response of the consensus control in the current control loop, the current feedforward (F is the feedforward gain) is adopted.

It is known that when subjected to power disturbance, the SC system can improve the equivalent inertia of DCMG and smooth the active power of DC voltage control unit without any delay, by simply tracking the DC bus voltage.

The function of the control loop of the base value regulation of SC terminal voltage in Fig. 4 is to provide an adaptive base value UscB to realize SC terminal voltage control and satisfy a given set value of equivalent power supporting time Ts simultaneously, which is described in detail in the following subsection.

2) Base Value Regulation of SC Terminal Voltage

Adjusting UscB online is the key to realizing the dynamic coupling of the SC terminal voltage control in power smoothing applications, which is also critical for DCMG to suppress uncertain power disturbances. To this end, the following method is proposed to set UscB online:

UscB=Uratedusc,set-uwork2Rd(1+Rd-udc,pu) (5)

where Urated is the rated value of the SC terminal voltage; and usc,set and uwork (in per-unit) are the setting value and allowable working area of the SC terminal voltage, respectively. In this paper, the setting value usc,set is selected to be 1/(1+Rd) in order to prevent usc from exceeding the rated value Urated.

In (5), there is a droop feature between UscB and the per-unit value of the DC bus voltage udc,pu. Specifically, the higher the DC bus voltage udc,pu is, the greater UscB will be. In this scenario, the operation state is mainly related to the low power level of the grid-following units Pp,pu (e.g., the output power of the RPGs is greater than the load consumption). Then, the main cause of power perturbation is from a reduction in the output of RPGs or an increase in the load. Therefore, with a larger base value of UscB, the SC system has more stored energy to cope with such transient power disturbances. In contrast, a lower DC bus voltage value udc,pu indicates that the DCMG is in a state of lower RPG output or that the load is heavier. Obviously, a lower value of UscB can ensure that the SC system has sufficient charging capacity after being subjected to an RPG power rush or a sudden load decrease. By substituting (5) into (4), the final relationship between the SC terminal voltage and DC bus voltage can be derived as:

usc=Urateduworkudc,pu22Rd+usc,set-1+Rd2Rduworkudc,pu (6)

It should be noted that the working area uwork in Fig. 4 and (6) play a key role in the real power smoothing application of the SC system, especially to satisfy a given set value of Ts. To this end, the value of uwork should be designed systematically, which is introduced in the following section.

IV. Analytical Solution of Ts by a Simple but Effective Reduced-order Model of DCMG

In this section, a simple but effective reduced-order model of DCMG based on the DC bus voltage response is derived and verified. Then, the analytical solutions to the equivalent inertia and equivalent power supporting time Ts contributed by the unified controlled SC system are obtained, which provide a theoretical foundation for the systematic design of the unified control of the SC system in Section V.

A. Reduced-order Model of DCMG

The dynamics of the DC bus voltage can be expressed as:

CbusUdcB2PBdΔudc,pudt=ΔPs,pu+ΔPsc,pu-ΔPp,pu (7)

where Cbus is the lumped capacitance; and CbusUdcB2/PB is regarded as the original inertia of the DCMG denoted as Hdc0.

With the droop control in Fig. 3, and by ignoring the dynamics of the inner current control loop, the dynamic response of the DC-AC converter can be obtained as:

ΔPs,pu=-Gudc(s)1+RdGudc(s)Δudc,pu (8)

where Gudc(s) is the proportional-integral (PI) controller of the DC voltage control loop.

To obtain the response of Δudc,pu to ΔPs,pu, (9) is defined based on the conservation of energy.

12Csc,0usc22-usc12=12Csc,equdc22-udc12 (9)

where Csc,0 is the actual capacitance of SC; and Csc,eq is the equivalent capacitance of SC from the DC bus side. Assuming that the DC bus voltage changes from udc1 to udc2 after the disturbance, the SC terminal voltage is changed from usc1 to usc2 accordingly.

Based on (4), (6), and (9), the detailed expression of Csc,eq can be obtained as:

Csc,eq=AratioCsc,0 (10)

where Aratio is a transformation coefficient expressed as:

Aratio=Urated2UdcB2uwork2Rd2udc1,pu2+udc2,pu2+usc,set-1+Rd2Rduwork2+uworkRdusc,set-1+Rd2Rduworkudc1,pu2+udc1,puudc2,pu+udc2,pu2udc1,pu+udc2,pu (11)

Then, according to (9)-(11), the dynamic performance of the power smoothing control of the SC system can be described in a very concise form as:

ΔPsc,pu=-Csc,eqUdcB2PBdΔudc,pudt=-HscdΔudc,pudt (12)

where Csc,eqUdcB2/PB is defined as the equivalent inertia contributed by the unified controlled SC system denoted as Hsc.

Combining (7), (8), and (12), the equivalent DC voltage response model of DCMG subjected to a power disturbance can be expressed as:

Δudc,pu=-Gs+zs2+2ξdωns+ωn2ΔPp,pu (13)

The related parameters in (13) are expressed as:

G=Rdkpu,s+1HeqRdkpu,s+Heqz=Rdkiu,sRdkpu,s+1ξd=HeqRdkiu,s+kpu,sHeqRdkpu,s+Heq/2kiu,sHeqRdkpu,s+Heqωn=kiu,sHeqRdkpu,s+HeqHeq=Hdc0+Hsc (14)

Using model (13), the dynamic performance index of the power smoothing control of the SC system can be calculated analytically.

B. Verification of Second-order Model

To validate the second-order model of the DCMG in (13), a detailed switching model of the DCMG is developed in PSCAD/EMTDC, with the main circuit and control parameters listed in Table I, where Cvsc is the DC capacitor of grid-forming unit at DC bus side; and Chl and Chs are the DC capacitors of grid-following unit and grid-supporting unit at DC bus side, respectively. For the switching model in PSCAD/EMTDC, the control period of the converters and simulation time step are selected as 100 μs and 1 μs, respectively. The comparison results of reduced-order model (13) and the switching model under different conditions and parameters are shown in Fig. 5.

TABLE I  Main Circuit and Control Parameters of DCMG System
UnitItemValue
AC-side circuit AC RMS line voltage Vs 380 V
Transformer 380 V/220 V
Rated frequency 50 Hz
DC-AC (grid-forming unit) LCL filter 2 mH/10 μF, 0.5 Ω/0.12 mH
Rated power PB 10 kW
Rated RMS line voltage Vinv 220 V
udc,set/iqref,pu 1/0
kpu,s/kiu,s/kpi,s/kii,s 16/160/2/145
Phase locked loop kpp/kip 10/100
Cvsc 4000 μF
DC-side circuit Rated DC voltage UdcB 400 V
SC DC-DC kpu,sc/kiu,sc/kpi,sc/kii,sc 4/30/0.05/50
Lsc/Chs 2 mH/2000 μF
Grid-following unit Rated power 10 kW
kpu,l/kiu,l/kpi,l/kii,l 1/20/0.05/50
Uref 300 V
L1/Chl/Cll 2 mH/2000 μF/2500 μF

Fig. 5  Comparison results of reduced-order model (13) and switching model under different conditions and parameters. (a) With different values of ΔPp,pu. (b) With different values of Csc,0. (c) With different values of uwork. (d) With different values of Rd.

Given the basic parameters Urated=200 V, Csc,0=3 F, uwork=0.3, the simulation results under different power disturbances are presented in Fig. 5(b), (c), and (d) which compare the results for different values of Csc,0, uwork, and Rd, respectively. It is observed that for a wide range of variations in the main parameters, the step responses of model (13) are somewhat consistent with the simulation results based on the detailed switching model, which verifies the effectiveness of model (13).

C. Further Reduction of Second-order Model

Assuming that the dynamics of second-order model in (13) are designed to be overdamped, the characteristic equation of (13) will have two negative real poles (denoted by -1/T1 and -1/T2) and derived as:

1T1=q2+12q2-4ωn21T2=q2-12q2-4ωn2 (15)

where the variable q is expressed as:

q=HeqRdkiu,s+kpu,sHeqRdkpu,s+Heq (16)

The limits of q and ωn as Heq goes to infinity can be readily obtained as:

limHeq+q=Rdkiu,s+kpu,sHeqRdkpu,s+1=Rdkiu,sRdkpu,s+1=zlimHeq+ωn=kiu,sHeqRdkpu,s+Heq=0 (17)

With increasing Heq, pole -1/T1 gets closer to zero z while another pole -1/T2 moves toward the origin, which is verified by the changing trajectory of the zero and poles of model (13), as shown in Fig. 6(a).

Fig. 6  Influence of system parameters. (a) Changing trajectory of zeros and poles when Hsc changes. (b) DC voltage response with different values of Csc,0 when uwork=0.3 and Rd=0.05. (c) DC voltage response with different values of uwork when Csc,0=5 F and Rd=0.05. (b) DC voltage response with different values of Rd when Csc,0=5 F and uwork=0.2.

Since the pole -1/T1 and zero -z are very close to each other when DCMG is overdamped, model (13) can be further reduced to a simple first-order model:

Δudc,pu=-zGT1s+1/T2ΔPp,pu=-Rd1/T2s+1/T2ΔPp,pu (18)

Figure 6(b), (c), and (d) compares the results of the step responses of (13) and (18) for different values of Csc,0, uwork, and Rd, respectively. It is observed that for a wide range of variations in the main parameters, the step responses of (18) are very consistent with those of (13), which verifies the first-order model (18). To this end, the simple model (18) can be used for the systematic analysis and design of the proposed unified control of the SC system.

According to (18), the analytical expression for the dynamic response of the DC bus voltage can be derived as:

Δudc,pu=-Rd(1-e-t/T2)ΔPp,pu (19)

Substituting (19) into (12), the power response of the SC system is obtained as:

ΔPsc,pu=(HscRd/T2)e-t/T2ΔPp,pu (20)

The energy released or absorbed by the SC system to support DCMG can be calculated as:

ΔEsmooth,pu=0ΔPsc,pudt=HscRdΔPp,pu (21)

Then, using (2) and (21), the analytical solution of the equivalent supporting time Ts can be derived as:

Ts=HscRd (22)

Finally, with the simple models of DCMG derived in this section, the analytical solutions to equivalent inertia Hsc and equivalent supporting time Ts to the SC system are obtained.

V. Systematic Analysis, Design of Proposed Unified Control, and Simulation Verification

In this section, with a detailed DCMG system and associated parameters, a systematic analysis and design method for the proposed unified power smoothing control is presented and further verified by detailed simulations.

A. DCMG Configuration

To verify the proposed unified control and the related theoretical analysis, a detailed DCMG was built using PSCAD/EMTDC, as shown in Fig. 7.

Fig. 7  DCMG system for PSCAD/EMTDC simulation.

The DC bus voltage is controlled via the grid-connected DC-AC converter with the droop method described in Fig. 3, which can be regarded as the grid-forming unit of the DCMG. The grid-following unit is simulated using a DC-DC converter and an equivalent load resistance. The load terminal voltage is controlled to a constant value, and the power disturbance is simulated by changing the equivalent resistance. The SC system is connected to the DC bus by a bidirectional DC-DC converter with the proposed unified control shown in Fig. 4.

B. Systematic Analysis and Design of Proposed Control

This subsection initially presents the influence of physical and control parameters on the equivalent inertia Hsc. Then, the method of determining the feasible region of the equivalent supporting time Ts of the SC system is presented. Finally, according to specific power smoothing requirements within the feasible region of Ts, a simple method is developed to design an appropriate working area uwork for the SC system.

1) Equivalent Inertia Hsc

In (10)-(12), it is observed that the equivalent inertia Hsc contributed by the SC system is mainly related to: ① physical parameters, including the capacitance Csc,0 of SC and its rated voltage Urated, power of grid-following unit Pp0,pu in the pre-disturbance state, and the variation ΔPp,pu; ② control parameters, including droop gain Rd and the working area uwork of the SC system.

1) Influence of physical parameters on Hsc

Let UdcB=400 V, Rd=0.05, and uwork=0.3. Figure 8(a) shows that SC can contribute a greater equivalent inertia Hsc, with increasing Urated and Csc,0. In contrast, increasing PB leads to a decrease in Hsc. The power supporting ability of the SC system decreases with increasing capacity of the DCMG.

Let Rd=0.05, uwork=0.3, Csc,0=10 F, Urated=200 V, and PB=10 kW. The influence of the DC bus voltages udc1,pu and udc2,pu (specifically associated with Pp0,pu and ΔPp,pu, respectively) on Hsc is shown in Fig. 8(b). As the droop gain Rd is 0.05, and the per-unit value of Pp,pu is assumed to be -1Pp,pu1, the variation in the DC voltages udc1,pu and udc2,pu is between 0.95 and 1.05, as shown in Fig. 8(b). It is observed that a high level of udc1,pu and udc2,pu will increase the value of Hsc. This conclusion will help in the specific design of the equivalent supporting time Ts and working area uwork of the SC system.

Fig. 8  Influence of physical parameters on Hsc. (a) Relationship between Hsc and Urated, PB, and Ccs,0. (b) Relationship between Hsc and udc1,pu and udc2,pu.

2) Influence of Control Parameters on Hsc

Let Csc,0=10 F, Urated=200 V, UdcB=400 V, and PB=10 kW. The influence of droop gain Rd and working area uwork on Hsc is shown in Fig. 9. In Fig. 9(a) and (b), when udc1,pu and udc2,pu are both selected to be 1, the value of Hsc increases as uwork increases from 0 to 0.9.

However, within the same variation range of uwork, there is an optimized value of uwork to make the SC system provide the largest equivalent inertia Hsc when the DC voltage is at a low level (related to udc1,pu=udc2,pu=0.95 p.u.). As observed in Fig. 9(b), the main influence of Rd on Hsc is that Hsc decreases with a larger value of Rd. However, when udc1,pu=udc2,pu=0.95 p.u., and uwork is close to 0.9, a larger value of Rd may increase the value of equivalent inertia Hsc.

2) Feasible Region of Equivalent Supporting Time Ts

In Section III-A, the concept of the equivalent power supporting time Ts is provided to evaluate the power smoothing ability of the SC system, as shown in (3). Then, using the proposed unified control and reduced-order model, an analytical solution of Ts is derived, as depicted in (22). For a given DCMG with a configured SC system, assuming that the main parameters UdcB, PB, Rd, Csc,0, and Urated are well-designed, appropriate control parameters (specifically referred to as the working area uwork of the SC system) are chosen for the proposed unified control, so that a given set value of Ts is satisfied is an issue. In this paper, the feasible region of Ts is first determined. According to (22), the critical point for determining the feasible region of Ts is to find the available region of Hsc.

Fig. 9  Influence of control parameters on Hsc.

Based on (11) and (12), Hsc can be represented as a function constrained by:

Hscf(udc1,pu,udc2,pu,uwork) (23)

As shown in Fig. 8(b), Hsc increases with increasing udc1,pu and udc2,pu. To this end, if udc1,pu is selected as the minimum value (referred to as 1-Rd), Hsc satisfies:

Hscf(1-Rd,udc2,pu,uwork) (24)

If the equivalent smoothing time Ts0 is calculated using (24), the real power smoothing performance can ensure that TsTs0.

In addition, with udc1,pu=1-Rd, the variation range of |ΔPp,pu|, udc2,pu, and uwork can be obtained as:

0|ΔPp,pu|2 p.u. (25)
1-Rdudc2,pu=udc1,pu+|ΔPp,pu|Rd1+Rd (26)
0uwork1/(1+Rd) (27)

Finally, Algorithm 1 is provided to calculate the feasible region of Ts.

Algorithm 1  : calculation of feasible region of Ts

Step 1: give the known parameters UdcB, PB, Rd, Csc,0, and Urated.

Step 2: set udc1,pu=1-Rd, and get the variation range of |ΔPp,pu|, udc2,pu, and uwork using (25)-(27).

Step 3: calculate the feasible region of Ts using (11), (12), and (22) according to the different values of |ΔPp,pu| (or udc2,pu) and uwork within the applicable variation range.

For example, for a given DCMG with UdcB=400 V, PB=10 kW, Rd=0.05, Csc,0=33 F, and Urated=200 V, the feasible region of Ts is obtained using Algorithm 1, as shown in Fig. 10.

Fig. 10  Feasible region of Ts with different values of |ΔPp,pu| and uwork.

As shown in Fig. 10, for a given working area uwork of the SC system, the value of Ts increases as the power disturbance |ΔPp,pu| increases. When ΔPp,pu remains unchanged, there is an optimized value of uwork to ensure that the SC system has the largest value of Ts.

In the real operation of a DCMG, the power disturbance |ΔPp,pu| is always uncertain, but can be known within a certain range. Thus, the working area uwork becomes a critical control parameter for the SC system to satisfy a given set value of Ts.

3) Design of Working Area uwork to Satisfy Power Smoothing Performance

Within the feasible region of Ts, the focus is on the following issue: when the value of an uncertain disturbance |ΔPp,pu| is greater than ΔPp,set (referred to a threshold value of power disturbance in the DCMG to be smoothed), how is an appropriate value of uwork designed to make the SC system to support the DCMG with the expected requirement denoted as TsTmin? Here, Tmin is the minimum equivalent supporting time requirement, which should be within the feasible region of Ts.

The design of an appropriate value of uwork is summarized in Algorithm 2. Here, two cases with different power smoothing performance requirements are considered for the application of Algorithms1 and 2, with the basic parameters listed in Table I.

Algorithm 2  : design of an appropriate value of uwork

Step 1: set |ΔPp,pu|=ΔPp,set, and obtain the relationship between uwork and Ts based on Algorithm 1.

Step 2: set the supporting time constraint Tmin, then obtain the feasible range of uwork (which can make Ts greater than the Tmin).

Step 3: select the largest value within the derived feasible range of uwork in Step 2 as the final design value.

1) Case 1: set Rd=0.05. And Ts satisfies TsTmin=15 s when the uncertain power disturbance |ΔPp,pu|ΔPp,set=0.25 p.u..

2) Case 2: set Rd=0.075. And Ts satisfies TsTmin=20 s when the uncertain power disturbance |ΔPp,pu|ΔPp,set=1.0 p.u..

From Step 1 in Algorithm 2, the feasible ranges of uwork satisfying the requirements of Case 1 and Case 2 are obtained as (0.2622, 0.6430) and (0.396, 0.696), respectively, as shown in Fig. 11(a). As long as the value of uwork is within the derived range, the preset power supporting time in the two cases is satisfied. Here, a method is proposed for selecting an optimal value.

Fig. 11  Influence of uwork. (a) Relationship between uwork and Ts. (b) Variation of A¯ratio associated with uwork.

In (10)-(12), Hsc is mainly associated with the transformation coefficient Aratio, which is related to the operation states udc1,pu and udc2,pu. In this paper, an average index of the transformation coefficient A¯ratio is defined as:

A¯ratio=1-Rd1+Rd1-Rd1+RdAratio4Rd2dudc1,pududc2,pu (28)

The parameter A¯ratio can be used to represent the average equivalent inertia and average power smoothing ability of the SC system under all possible operating conditions of the DCMG. Figure 11(b) shows the relationship between A¯ratio and uwork in Case 1. It is observed that the value of A¯ratio is positively proportional to that of uwork. Therefore, after determining the feasible range of uwork, the largest value of uwork within the feasible range can be selected for a more effective utilization of the SC system, which corresponds to Step 3 in Algorithm 2.

4) Actual Equivalent Power Supporting Time Ts with Designed Parameter of Working Area uwork

With the designed value of uwork, the actual equivalent power supporting time Ts contributed by the SC system can be determined according to (10)-(12) and (22) under all operating conditions of the DCMG. Taking Case 1 as an example, the feasible range of uwork is (0.2622, 0.6430). Then, uwork=0.6430 is considered according to Algorithm 2. To verify the effectiveness of this selection, the actual power smoothing performance represented by Ts is compared with uwork=0.6430 and uwork=0.2622, as shown in Fig. 12.

Fig. 12  Actual equivalent supporting time Ts with uwork=0.6430 and uwork=0.2622 in Case 1.

In Fig. 12, for most operating states of the DCMG, the value of the actual equivalent power supporting time Ts of the SC system for uwork=0.6430 is larger than that of uwork=0.2622, which means that the SC system can provide more equivalent inertia, verifying the effectiveness of Step 3 in Algorithm 2.

C. Simulation Verification

Based on the switching model in PSCAD/EMTDC, the proposed unified control of the SC system and previous theoretical analysis results are verified in this subsection.

1) Case 1

In this case, the expected power smoothing performance and equivalent supporting time Ts are the same as those in Section V-B, which is summarized as: setting Rd=0.05, Ts satisfies TsTmin=15 s when the uncertain power disturbance |ΔPp,pu|ΔPp,set=0.25 p.u..

Based on the previous analysis, the feasible range of uwork satisfying the requirement of Case 1 is (0.2622, 0.6430), as shown in Fig. 11(a). Figure 13 shows the simulation results of Case 1, for uwork=0.6430 and uwork=0.2622. Figure 13(a)-(e) shows the simulation results of Ps,pu (active power of DC-AC converter), udc,pu (DC bus voltage), Psc,pu (power of SC system), usc,pu (SC terminal voltage), and Ts, respectively.

When t<200 s, the DC bus voltage udc,pu is close to 1 p.u.. When uwork is selected to be 0.6430 and 0.2622, the SC terminal voltage value usc,pu is approximately 0.631 and 0.821 p.u. (the rated value of the SC terminal voltage Urated is 200 V), respectively, as shown in Fig. 13(d).

At t=200 s, the power of grid-following unit Pp,pu is changed to 0.5 p.u., with ΔPp,pu of 0.5 p.u.. As shown in Fig. 13(c), the SC system responds immediately to support the DCMG and balance the related power disturbance, which means that the expected power smoothing performance shown in Fig. 2(a) is realized with the proposed unified control in Fig. 4. Moreover, as shown in Fig. 13(a), (b), and (d), the power of the grid-forming units Ps,pu, DC bus voltage udc,pu, and the terminal voltage of the SC system usc,pu reach a new steady state smoothly, which indicates that with the proposed unified control, active power smoothing of the DC voltage control unit and SC terminal voltage control are realized simultaneously. In Fig. 13(c), when uwork1=0.6430 and uwork2=0.2622, the actual equivalent power supporting time during this transient is Ts1=24.88 s and Ts2=17.36 s, respectively. Obviously, Ts satisfies the given requirement TsTmin=15 s when the power disturbance |ΔPp,pu|ΔPp,set=0.25 p.u., verifying the effectiveness of Algorithm 2 in Section V. Moreover, from the result Ts1>Ts2, it is observed that with a larger value of working area uwork, the SC system can operate for a larger value of equivalent supporting time, which is consistent with the theoretical analysis shown in Fig. 12.

Fig. 13  Simulation results of Case 1. (a) Active power of DC-AC converter. (b) DC bus voltage. (c) Power of SC system. (d) SC terminal voltage.

At t=350 s, Pp,pu is increased from 0.5 p.u. to 0.75 p.u., with ΔPp,pu=0.25 p.u.. From the results shown in Fig. 13(c), it can be observed that the actual Ts of the SC system is approximately Ts1=18.31 s and Ts2=15.49 s, respectively. The results also indicate that the designed value of uwork can achieve the expected power smoothing performance in this scenario.

Then, Pp,pu is increased from 0.75 p.u. to 1 p.u. at t=450 s and recovered back to 0.75 p.u. at t=550 s. In Fig. 13(c), it is observed that the actual Ts is approximately 14.3 s, which is slightly less than the set time of 15 s. During t=450-550 s, the DC bus voltage udc,pu reaches the lowest value of 0.95 p.u., and the actual power smoothing performance is affected, which is consistent with the theoretical analysis in Fig. 12.

2) Case 2

This case has the same conditions as that of Case 2 in Section V-B, which is re-emphasized as follows.

Setting Rd=0.075, Ts satisfies TsTmin=20 s when the uncertain power disturbance |ΔPp,pu|ΔPp,set=1 p.u..

From the theoretical analysis shown in Fig. 11(a), the feasible range of uwork that satisfies the requirement is between 0.396 and 0.696. In particular, the equivalent supporting time Ts reaches the maximum value of 20.96 s, when uwork is selected to be 0.552. Figure 14 shows the simulation results of Case 2, considering the three conditions of uwork1=0.696, uwork2=0.396, and uwork3=0.522. Figure 14(a)-(d) shows the simulation results of Ps,pu (active power of DC-AC), udc,pu (DC bus voltage), Psc,pu (power of SC system), usc,pu (SC terminal voltage), and Ts, respectively.

Fig. 14  Simulation results of Case 2. (a) Active power of DC-AC converter. (b) DC bus voltage. (c) Power of SC system. (d) SC terminal voltage.

When t<200 s, the DC bus voltage udc,pu is approximately 1 p.u.. When uwork is selected to be 0.696, 0.396, and 0.522, the SC terminal voltage values are usc1,pu=0.5822 p.u., usc2,pu=0.7322 p.u., and usc3,pu=0.6542 p.u., respectively (the rated value of the SC terminal voltage Urated is 200 V), as shown in Fig. 14(d).

At t=200 s, Pp,pu is increased to 1 p.u. with ΔPp,pu of 1 p.u.. As illustrated in Fig. 14(c), with the proposed unified control, the SC system discharges rapidly to support the DCMG at the moment of power disturbance, without any delay. In Fig. 14(a) and (b), the power output of the DC-AC converter (referred to as the grid-forming unit) and the DC bus voltage are smoothed as expected. According to the simulation results in Fig. 14(c), when uwork is selected to be 0.696, 0.396, and 0.522, the actual power supporting time Ts is Ts1=19.34 s, Ts2=19.27 s, and Ts3=20.2 s, respectively, which is very close to the designed value of Tmin=20 s. Thus, the simulation results are consistent with the theoretical analysis demonstrated in Fig. 11(a), which verifies the effectiveness of the systematic analysis and design method for the proposed unified control in Section V-B.

At t=350 s, Pp,pu is reduced to be 0.5 p.u. from 1 p.u., with a power disturbance |ΔPp,pu| of 0.5 p.u.. As the preset uncertain power disturbance is ΔPp,set=1.0 p.u., the actual Ts under a relatively small disturbance |ΔPp,pu| may not satisfy the requirement TsTmin=20 s. The simulation results in Fig. 14(c) show that the actual Ts associated with the control parameter uwork is about Ts1=14.18 s, Ts2=16.8 s, and Ts3=16.43 s, respectively, which is consistent with the previous theoretical analysis.

D. Discussion About Control

Two cases with different parameters and conditions are applied to verify the proposed unified control of the SC system and the theoretical analysis in this paper. From the simulation results, the following conclusions can be drawn.

1) The proposed unified control realizes power smoothing, and SC terminal voltage control, simultaneously on the same timescale. Compared with the methods in [

14]-[17], the unified control has a more simple but effective structure. Moreover, it should be noted that only local DC voltages and currents are used for power smoothing control, indicating good plug-and-play features.

2) The proposed simple reduced model of DCMG in Section IV and theoretical analysis in Section V are effective since the presented Algorithms1 and 2 derived from the analytical solutions of equivalent inertia and power supporting time have been verified by detailed simulation results.

In addition, it should be noted that the analytical solutions to equivalent inertia and power supporting time can be used not only for the design of the working area uwork of the SC system, but also for the optimal configuration of the SC system. For example, to satisfy the power smoothing performance, the questions such as how much the capacitance of the SC system is needed and how to determine the rated terminal voltage of the SC system can also be addressed with the proposed control, simple model, and theoretical results.

VI. Conclusion

In this paper, a novel unified control of the SC system is proposed for power smoothing in DCMG. Without any communication, it can simultaneously realize the smoothing of large power disturbance and terminal voltage control of the SC system, indicating satisfactory plug-and-play characteristics. Moreover, using the defined equivalent power supporting time and equivalent inertia contributed by the unified controlled SC system, a systematic design method for critical control parameters is presented and further verified by detailed simulation results. With the proposed unified control, reduced model, and design method, the expected power smoothing performance of the SC system represented by the equivalent power supporting time can be achieved in the feasible region, which is very useful for the practical operation of DCMGs.

References

1

H. Zhang, W. Xiang, W. Lin et al., “Grid forming converters in renewable energy sources dominated power grid: control strategy, stability, application, and challenges,” Journal of Modern Power Systems and Clean Energy, vol. 9, no. 6, pp. 1239-1256, Nov. 2021. [Baidu Scholar] 

2

K. Ahmed, M. Seyedmahmoudian, S. Mekhilef et al., “A review on primary and secondary controls of inverter-interfaced microgrid,” Journal of Modern Power Systems and Clean Energy, vol. 9, no. 5, pp. 969-985, Sept. 2021. [Baidu Scholar] 

3

C. Wang, Y. Mi, Y. Fu et al., “Frequency control of an isolated micro-grid using double sliding mode controllers and disturbance observer,” IEEE Transactions on Smart Grid, vol. 9, no. 2, pp. 923-930, Mar. 2018. [Baidu Scholar] 

4

Z. Zhao, P. Yang, J. M. Guerrero et al., “Multiple-time-scales hierarchical frequency stability control strategy of medium-voltage isolated microgrid,” IEEE Transactions on Power Electronics, vol. 31, no. 8, pp. 5974-5991, Aug. 2016. [Baidu Scholar] 

5

Y. Li, G. Fu, T. An et al, “Power compensation control for interconnection of weak power systems by VSC-HVDC,” IEEE Transactions on Power Delivery, vol. 3, no. 4, pp. 1964-1974, Aug. 2017. [Baidu Scholar] 

6

Y. Wang, C. Wang, L. Xu et al., “Adjustable inertial response from the converter with adaptive droop control in DC grids,” IEEE Transactions on Smart Grid, vol. 10, no. 3, pp. 3198-3209, May 2019. [Baidu Scholar] 

7

W. Wu, Y. Chen, A. Luo et al., “A virtual inertia control strategy for DC microgrids analogized with virtual synchronous machines,” IEEE Transactions on Industrial Electronics, vol. 64, no. 7, pp. 6005-6016, Jul. 2017. [Baidu Scholar] 

8

S. Samanta, J. P. Mishra, and B. K. Roy, “Virtual DC machine: an inertia emulation and control technique for a bi-directional DC-DC converter in a DC microgrid,” IET Electric Power Applications, vol. 12, no. 6, pp. 874-884, Jul. 2018. [Baidu Scholar] 

9

K. Bellache, M. B. Camara, and B. Dakyo, “Transient power control for diesel-generator assistance in electric boat applications using supercapacitors and batteries,” IEEE Journal of Emerging and Selected Topics in Power Electronics, vol. 6, no. 1, pp. 416-428, Mar. 2018. [Baidu Scholar] 

10

L. He, Y. Li, J. M. Guerrero et al., “A comprehensive inertial control strategy for hybrid AC/DC microgrid with distributed generations,” IEEE Transactions on Smart Grid, vol. 11, no. 2, pp. 1737-1747, Mar. 2020. [Baidu Scholar] 

11

Z. Jin, L. Meng, J. M. Guerrero et al., “Hierarchical control design for a shipboard power system with DC distribution and energy storage aboard future more-electric ships,” IEEE Transactions on Industrial Informatics, vol. 14, no. 2, pp. 703-714, Feb. 2018. [Baidu Scholar] 

12

Q. Xu, X. Hu, P. Wang et al., “A decentralized dynamic power sharing strategy for hybrid energy storage system in autonomous DC microgrid,” IEEE Transactions on Industrial Electronics, vol. 64, no. 7, pp. 5930-5941, Jul. 2017. [Baidu Scholar] 

13

P. Lin, P. Wang, J. Xiao et al., “An integral droop for transient power allocation and output impedance shaping of hybrid energy storage system in DC microgrid,” IEEE Transactions on Power Electronics, vol. 33, no. 7, pp. 6262-6277, Jul. 2018. [Baidu Scholar] 

14

Q. Xu, J. Xiao, P. Wang et al., “A decentralized control strategy for autonomous transient power sharing and state-of-charge recovery in hybrid energy storage systems,” IEEE Transactions on Sustainable Energy, vol. 8, no. 4, pp. 1443-1452, Oct. 2017. [Baidu Scholar] 

15

X. Chen, J. Zhou, M. Shi et al., “A novel virtual resistor and capacitor droop control for HESS in medium-voltage DC system,” IEEE Transactions on Power Systems, vol. 34, no. 4, pp. 2518-2527, Jul. 2019. [Baidu Scholar] 

16

J. Chen and Q. Song, “A decentralized dynamic load power allocation strategy for fuel cell/supercapacitor-based APU of large more electric vehicles,” IEEE Transactions on Industrial Electronics, vol. 66, no. 2, pp. 865-875, Feb. 2019. [Baidu Scholar] 

17

J. Pegueroles-Queralt, F. D. Bianchi, and O. Gomis-Bellmunt, “A power smoothing system based on supercapacitors for renewable distributed generation,” IEEE Transactions on Industrial Electronics, vol. 62, no. 1, pp. 343-350, Jan. 2015. [Baidu Scholar] 

18

J. Rocabert, A. Luna, F. Blaabjerg et al., “Control of power converters in AC microgrids,” IEEE Transactions on Power Electronics, vol. 27, no. 11, pp. 4734-4749, Nov. 2012. [Baidu Scholar] 

19

L. Che, M. Shahidehpour, A. Alabdulwahab et al., “Hierarchical coordination of a community microgrid with AC and DC microgrids,” IEEE Transactions on Smart Grid, vol. 6, no. 6, pp. 3042-3051, Nov. 2015. [Baidu Scholar] 

20

X. Li, L. Guo, Y. Li et al., “A unified control for the DC-AC interlinking converters in hybrid AC/DC microgrids,” IEEE Transactions on Smart Grid, vol. 19, no. 6, pp. 6540-6553, Nov. 2018. [Baidu Scholar]