A Dynamic and Cooperative Control Strategy for Multi-Hybrid Energy Storage System of DC Microgrid Based on SOC

With the increasingly serious crisis of fossil energy and environmental pollution, clean renewable energy becomes the inevitable choice of energy structure adjustment. The instability of output power of distributed renewable energy system greatly affects the operation of DC microgrid. The hybrid energy storage system (HESS) composed of High-Energy Battery (HEB) and High-Power Battery (HPB) can solve the above problems. Thus, this paper proposes a dynamic and cooperative control strategy for multi-HESS based on state of charge (SOC). Based on the traditional LPF method and droop control, this paper proposes a control strategy that requires no communication among multiple hybrid energy storage (HES) modules. This method can realize the stable control of HEB current, reduce the change times of HEB charging-discharging mode, prevent HEB from overcharging and overdischarging, prolong the service life of HEB and balance different energy storage SOC, so as to improve the operation stability and economy of DC microgrid. In addition, the method has certain robustness against sudden failures. Simulation and experiment results show the effectiveness of the proposed method.


INTRODUCTION
The increasing penetration rate of renewable energy such as photovoltaic and wind power promotes the development of DC microgrid (Kathiresan et al., 2020;Zhou et al., 2020). Compared with AC microgrid, DC microgrid has no problems such as synchronization, reactive power transmission, harmonic current and converter loss, so it has attracted more and more attention (Rahimi and Ghadiriyan, 2019;Song et al., 2019).
However, due to the intermittent impact of renewable energy distributed generation and the short-term impact of some large loads, the stability of DC microgrid is severely challenged. Therefore, the corresponding energy storage system should be equipped to enhance the antiinterference ability of DC microgrid (Singh and Lather, 2021). On the one hand, combining High-Energy Battery (Such as lithium battery, lead-acid battery, sodium sulfur battery) with High-Power Battery (Such as supercapacitor, flywheel energy storage, super-magnetic energy storage) to form HESS is an important way of energy storage configuration based on the existing energy storage technology and satisfying the demand for DC microgrid (Kotra and Mishra, 2019;Mathews and Rajeev, 2020). On the other hand, in order to improve the access flexibility of energy storage system, and avoid the cost increase of using large-capacity converters in centralized energy storage system and the problems of operational reliability, distributed multi-HESS should be used (Su et al., 2018).
In terms of applying HESS to suppress grid-connected power fluctuations of distributed power generation, many scholars have studied various control strategies according to the characteristics of different types of energy storage, such as LPF method (Manandhar et al., 2015), fuzzy logic control (FLC) (Musilek et al., 2017;Mathews and Rajeev, 2020), wavelet decomposition method (Chiang et al., 2017), layered drooping control (Li et al., 2016;Zhang et al., 2018), virtual inertial control (Ming et al., 2017;Wang et al., 2017), machine learning (Chen et al., 2020) and so on. The above method mainly focuses on centralized HESS, when extended to distributed multi-HESS, they have some inapplicability. With the increase of using time, the differences among HES modules become larger and larger. If not restrained and controlled, the performance of HES modules will be seriously affected.
The concept of a virtual power rating was introduced in Literature (Hoang and Lee, 2019) to achieve accurate power sharing between batteries. However, this control method requires communication, and the scope of its application has certain limitations. In literature (Zhou et al., 2018), multi-HEB SOC is divided into five regions, and the whole system is divided into six working modes. However, this method may switch modes frequently under the circumstance of high-frequency fluctuation, leading to system instability. Most importantly, the above literatures mainly focus on the cooperative control of multiple energy storage units with the same medium and the same property.
In view of the collaborative control problem of distributed multi-HESS in DC microgrid, a power control strategy for multi-HESS was proposed based on the consensus protocol in literature . The communication was established among different energy storage units, and the SOC of different energy storage units was balanced by the consensus protocol, thus effectively solving the problem of balanced power distribution among multi-HESS. In literature (Wu et al., 2020), the control system automatically switches the operating mode based on the DC bus voltage, the supercapacitor voltage and the accumulator state. However, these methods still depend on external communication and are not suitable for situations requiring high reliability. To solve this problem, literature (Chen et al., 2016) proposed voltage-current droop control for HEB and voltage-change rate-current droop control for HPB. When the system only relies on the local information of the energy storage, it realizes the automatic power distribution between the components of different types of energy storage. However, this method is controlled by difference, and the deviation of bus voltage is large.
In order to solve the problem of high dependence on external communication in power distribution, this paper takes the island operation of DC microgrid as the research object, based on the traditional LPF method and droop control, proposes a SOC based power secondary distribution method without communication among HES modules. The control system automatically redistributes power according to SOC of different energy storage to achieve SOC balance. The smooth control of HEB current can be realized by using the characteristic that HPB voltage cannot change suddenly, so as to reduce the number of charging-discharging mode changes of HEB and extend the service life of HEB. The effectiveness of the proposed control strategy for distributed multi-hybrid energy storage module parallel system is verified by simulation and experiment.

SYSTEM MODEL
The schematic diagram of DC microgrid with multi-HESS is shown in Figure 1, which mainly includes renewable energy power generation unit, AC/DC load and energy storage unit. Each part is a distributed structure, and each unit is connected to the DC bus through the corresponding converter. The power relationship of each part is as follows: P Gen is the power generated by renewable energy, P HESS is the total power of multi-HESS, and P Load is the total load power, including AC load and DC load.
In system analysis and calculation, the load of renewable energy power generation and power electronic equipment access can be equivalent to a controllable current source, which can be positive or negative. When it is positive, the renewable energy power generation is greater than the load; when negative, the load is greater than the renewable energy power generation. HESS in Figure 1 consists of High-Energy Battery and High-Power Battery. High-Energy Battery has high energy density and long storage time, but low power density and low cycle times. High-Power Battery has high power density and high cycle times, but low energy density and short energy storage time.
The single HES module topology adopted in this paper is shown in Figure 2. U e and U p , R e and R p are the voltage and resistance of High-Energy Battery and High-Power Battery respectively. U dc is the DC bus voltage. The converter in the Figure 2 is a bidirectional Buck-Boost converter. C e and C p , L e and L p are respectively the filter capacitance and inductance of the converters. R le and R lp are the internal resistance of respective inductors. S1-S4 are power tubes. C dce and C dcp are filter capacitors of converters respectively. C o is the bus voltage regulator capacitor. (Subscript "e" means High-Energy Battery, and subscript "p" means High-Power Battery.)

COOPERATIVE CONTROL STRATEGY FOR MULTI-HESS
The overall architecture of multi-HESS collaborative control strategy is shown in Figure 3. The improved LPF control structure is adopted to regulate the power distribution within the single HES module. Based on the HPB SOC, the output power of HEB and HPB inside the module is adjusted to maintain the HPB SOC in a healthy range as far as possible. Among the multiple HES modules, the output of the whole module is regulated through the reference bus voltage adaptive regulator based on HEB SOC, so as to maintain the overall available capacity of each HES module relatively consistent. Since each HES module adopts V-I droop control, the coordinated control among HES modules does not require communication, which has high reliability and scalability, plug and play, and is convenient for installation and deployment.

Basic Control Structure
In microgrid, there are many control methods for converters, among which droop control is a widely used decentralized control method. It does not need communication, and only uses local information to realize the coordination control between distributed power supplies. It is suitable for DC microgrid with high reliability requirements, such as islands, reefs and ships, etc. Therefore, in order to maximize the autonomous control of each HES module, and enhance the system reliability and "plug and play" capability, the converter of HESS adopts droop control to stabilize the DC bus voltage.
V-I droop control (Li et al., 2016;Zhang et al., 2018) can be expressed as: U ref is the value of DC bus voltage under no load, U op is the specified value of converter output voltage after correction, and R V is the droop coefficient.
The value of R V is determined by Eq. 3: FIGURE 2 | Single HES module topology.
FIGURE 3 | The overall architecture of multi-HESS collaborative control strategy.
Frontiers in Energy Research | www.frontiersin.org January 2022 | Volume 9 | Article 795513 U omin is the allowable minimum steady state value of DC bus voltage, and I omax is the maximum current value that DC/DC converter can output/input.
Ignoring the line impedance, the output current relationship of several DC/DC converters with V-I droop control is shown as follows.
The control structure of High-Power Battery consists of a droop control outer loop and a voltage and current double closedloop, as shown in Figure 4. U ref , U bus and I ref_p respectively represent the reference voltage of the bus, the actual voltage of the bus and the reference current of the High-Power Battery.
HEB is mainly used to follow the current instructions given by the power distribution controller, as shown in Figure 5. I ref_e and i e respectively represent the reference current and the actual current of the High-Energy Battery. Figure 6 is the control structure diagram of LPF-based. Reasonable power distribution method is the key to the stability and reliability of HESS. The traditional LPF method has less computation and parameters, and the control is simple. It can realize the power distribution based on frequency without any communication. Thus, this paper adopts LPF method for power distribution. On the basis of the traditional LPF control, the adaptive power redistribution module is added to regulate the output of different energy storage units in real time and balance the remaining electric quantity. The adaptive power redistribution modules represented by the two blue squares are described in detail in the next section.

HES Module Control Based on LPF
iDC HESSi is the sum of the current at the bus end of the DC/DC converter in the ith HES module, as shown in Eq. 5. i ei and i pi are the current of High-Energy Battery and High-Power Battery in the ith HES module, while d ei and d pi are the control signals of DC/DC converter. The limiter prevents the reference current from exceeding the converter maximum current.
iDC HESSi iDC ei + iDC pi (5) iDC ei and iDC pi are the current at the bus end of the DC/DC converter connected with High-Energy Battery and the High-Power Battery in the ith HES module.

Overall Idea of Adaptive SOC Layered Control
During the use of the HESS, there are different initial SOC and rated capacity of each energy storage unit connected by the converter in a certain working condition. If SOC are not considered, some energy storage units may overcharge or overdischarge, affecting the stable operation of the whole system. Therefore, the balanced control of energy storage SOC is the key to the stable and economical operation of HESS.  In order to prevent energy storage units from overcharging/ overdischarging and frequent charging/discharging change, this paper adopts adaptive SOC layered control strategy. According to the SOC level of each energy storage unit, it is divided into different operating range as shown in Figure 7.
(1) Free charging/discharging zone. When SOC low ≤SOC≤SOC high , the energy storage unit can conduct normal charge/discharge in accordance with the droop relationship, without imposing other control operations based on SOC.
(2) Charging limit zone. When SOC high ≤SOC≤SOC max , the energy storage charging will be limited, and as SOC gets closer and closer to SOC max , the charging limit becomes more and more serious until no more charging is allowed. (3) Discharging limit one. When SOC min ≤SOC≤SOC low , the energy storage discharging will be limited, and as SOC gets closer and closer to SOC min , the discharging limit becomes more and more serious until no more discharging is allowed.

Adaptive SOC Control Method Inside the HES Module
In HESS, power fluctuation is usually divided into low frequency and high frequency, which are borne by High-Energy Battery and High-Power Battery respectively. The response speed of High-Power Battery is obviously faster than that of High-Energy Battery. DC microgrid is operating in island mode, the High-Power Battery will quickly compensate the power shortage in case of sudden change of load. When the High-Power Battery responds quickly, its SOC fluctuates. Then, the High-Energy Battery responds slowly to maintain the stability of the High-Power Battery SOC, achieving a dynamic adaptive balance. Through such a master-slave double structured with adaptive control, the system will gradually transfer the power deficiency firstly borne by the High-Power Battery to the High-Energy Battery, which can make the High-Energy Battery charge/discharge change relatively smooth, reduce the depth of the charge/ discharge, prolong the life of High-Energy Battery, and avoid overcharge/overdischarge of High-Power Battery. The ampere-hour integral method is used to calculate the SOC of High-Power Battery (Manandhar et al., 2015).
SOC p is the real-time SOC of High-Power Battery, SOC p0 is the initial SOC, T is the running time, μ p is the charging/discharging efficiency, i p is the charging/discharging current, C Np is the rated capacity of the HPB.
According to the layered control principle introduced in Figure 5, partial charging/discharging current of High-Power Battery is transferred to the High-Energy Battery. This paper designs I SOC_adj that is the current of the High-Power Battery transferred to the High-Energy Battery based on the concept of partition operation as shown in Figure 7. The specific expression is shown in Eq. 7: I ref_P refers to the initial reference current of the High-Power Battery, and V diff refers to the difference between the reference voltage and the actual voltage of the bus, which is used to characterize the charging and discharging state of energy storage unit at this time.
V diff > 0 means the energy storage unit is in charge, and V diff < 0 means it is in discharge.
I ref_p* and I ref_e* represent the regulated High-Power Battery and High-Energy Battery reference current. According to Eqs 7-10, when the High-Power Battery is in the charging limit zone, a certain proportion of the charging current is transferred to the High-Energy Battery according to the current SOC p . The higher the SOC p is, the larger the proportion of the transferred current will be. On the contrary, when the High-Power Battery is in the discharging limit zone, situation is the opposite. Finally, the High-Power Battery SOC is always in a reasonable operating range.

Among the HES Modules
High-Power Battery has low energy density. The adaptive regulating frequency of power among HES modules is much lower than that in a single HES module, so the High-Energy Battery SOC can be used as an indicator of the whole HES module available capacity.
The ampere-hour integral method is used to calculate the SOC of High-Energy Battery (Manandhar et al., 2015).
SOC e is the real-time SOC of High-Energy Battery, SOC e0 is the initial SOC, T is the running time, μ e is the charging/ discharging efficiency, i e is the charging/discharging current, C Ne is the rated capacity of the HEB.
The proposed adaptive regulation strategy of reference voltage based on SOC is used to regulate the reference voltage of droop control model.
V SOC_adj is the reference voltage regulation. In order to improve the regulation performance of the strategy and avoid system oscillation caused by re-regulation of critical SOC, the layered regulation framework shown in Figure 7 is adopted. The specific algorithm is as follows: k is the proportional amplification gain, which is determined by the maximum power of the converter and the requirement of the system on the regulation ability. According to Eqs 12, 13 when the SOC of an energy storage unit is in the charging limit zone, the reference voltage value of the DC/DC convert output should be increased to reduce the charging current of the energy storage unit. On the contrary, when the SOC is in the discharging limit zone, the situation is the opposite. Finally, SOC of different energy storage units tends to be relatively consistent. The specific process is shown in Figure 8, where the X-axis is the SOC of energy storage unit and the Y-axis is the regulator of the reference voltage value of the DC/DC convert output. The model is simple and easy to implement, and the switching between modes is relatively smooth.

Stability Analysis Model
Because the control strategy proposed in this paper presents a coupling state inside a single HES module, it is difficult to adopt impedance analysis method. Thus, this paper intends to adopt small-signal analysis method to establish a mathematical analysis model for a single HES module as a whole.
Since resistive load increases system damping, constant power load (CPL) decreases system damping. So in the worst case, assume that the load on the system contains only CPL.
Once the system is stable in the worst case, it is stable in all cases. Thus, in the study, the load will select the CPL, which can be expressed as: i Load is load current, P CPL is load power, u dc is bus voltage transient value. The small signal model can be obtained as follows: U dc is steady-state value of bus voltage at equilibrium point. The state equation of the main circuit is shown in Eq. 16: d e and d p are the duty cycle instantaneous value of HEB and HPB control signal. The small signal model can be obtained as follows: Eq. 17 is written in matrix form: The small-signal model of the HEB control loop in Figure 6: X v and X ie are the steady-state values of the input of the voltage loop and the current loop in the control loop of HEB. I DCref is the total current reference value, I ref_e is the current reference value of HEB, T f is the filtering time constant of the low-pass filter, k vp and k vi , k ep and k ei are the PI parameters of the voltage loop and the current loop in the control loop of HEB.
The deformation of Eq. 21 can be obtained as follows: The duty cycle of HEB converter can be expressed as: The control loop of HPB can be expressed as: X ip is the steady-state input value of the current loop in the HPB control loop, I ref_p is the current reference value of the HPB, and k pp and k pi are the PI parameters of the current loop in the HPB control loop.
In combination with Eqs 24, 27, the small-signal model can be written as: The state vector Δx 2 of the controller variable can be expressed as: According to Eqs 12-22 and 25, 26, we can get: The upper control input vector Δx 3 can be expressed as: The energy storage unit has an integral relation between SOC and time. Despite the existence of such an integral relationship, the present state value of the SOC is calculated and output by sampling, so the changes of the output current and voltage are difficult to reflect the changes of the SOC in very short time. Therefore, the SOC is only treated as an input signal in small signal modeling. At the same time, considering that the change of SOC is not as rapid Frontiers in Energy Research | www.frontiersin.org January 2022 | Volume 9 | Article 795513 as the change of voltage and current, delay sampling strategy is generally adopted, and LPF can be used to replace short delay.
As V SOC_adj and I SOC_adj are both piecewise functions in the optimization of upper power allocation considering SOC. In the stability analysis, only the case with the worst stability is considered, that is, they are all in the discharge limit region: Small signal transformation is performed for Eqs 31, 32: w is the cutoff frequency of LPF. In summary, the following formula can be obtained by combining Eqs 19, 20 and 25:    Figure 9 is the closed-loop pole distribution diagram of Eq. 35 when different PI parameters are selected (the black arrow indicates the increasing direction of the parameters taken). SOC limited boundaries of High-Energy Batteries and High-Power Batteries SOC min , SOC max , SOC low , SOC high respectively are 50, 70, 56, 63, and 40, 80, 55, 65%. Other specific simulation parameters are shown in Table 1. It can be seen from Figure 9 that when k vi and k vp change in the specified range, the system pole is always on the left of the Y-axis, so the system is always stable. As shown in Figure 9A, the poles close to the imaginary axis are conjugate poles when the voltage outer loop proportionality parameter k vp is less than 6, and it is an underdamped system. At the same time, with the increase of k vp , the damping ratio of the system gradually increases, and the overshoot decreases. Therefore, when k vp = 5, the overshoot of the system is small and the response speed is fast. As shown in Figure 9B, when k vi is greater than 30, the closed loop dominates the pole conjugate. At this time, the system has good dynamic response characteristics. As k vi continues to increase, the peak value of the system decreases, but the response time increases. Therefore, the trade-off is k vi = 35.

Stability Analysis and Control Parameter Selection
The same method is adopted for other control parameters and their stability analysis. The PI parameters of the proposed method are obtained and summarized as shown in Table 2.

SIMULATION AND EXPERIMENT
In this paper, PSCAD/EMTDC software will be used for simulation. In order to fully verify the performance of the proposed HESS control strategy, three groups of HESS will be set up. HESS1 is overdischarged, HESS2 is in the free charging/discharging zone, and HESS3 is overcharged. The structure of each group is shown in Figure 2. The overall simulation model is shown in Figure 10. In this simulation, lithium battery is used for HEB and supercapacitor is used for HPB. Simulation parameters are shown in Tables 1, 2.
Step Load Change In this case, the response of the system to consecutive step load changes is simulated to observe the net power decomposition between High-Energy Battery and High-Power Battery. The specific waveform is shown in Figure 11A. At 10 s, the load current increases from 50 to 125 A, then drops back to 50 A at 25 s, and finally jumps to 150 A at 50 s. The initial SOC of three High-Energy Batteries and High-Power Batteries are 55, 62, 67, 60, 70%. To verify the effectiveness of the proposed strategy, it will be compared with the LPF-based HESS control method without SOC regulation, that is, the control method shown in Figure 5 removes the SOC regulation module represented by the blue square. The results are shown in Figures 11B-K. It can be seen from Figures 11B,C that both methods can effectively cope with the step change of load and make the deviation of bus voltage stable within 2.5%. However, the voltage deviation rate of the proposed method is reduced by about 50%.
Figures 11D-G shows that SOC of each energy storage unit cannot be balanced without SOC regulation. Considering the practical differences of energy storage units, SOC difference will become larger and larger, eventually leading to the overcharge or overdischarge of some energy storage units. However, under the SOC coordinated control method proposed in this paper, SOC can tend to be consistent which can effectively prevent some energy storage units from overcharging/overdischarging, and improve the stability and reliability of the system.
According to the comparison among Figures 11H-K, it can be seen that under the SOC regulation, the energy storage unit with SOC higher than the limit increases its output power when discharging. The energy storage unit with SOC lower than the limit decreases output power when discharging. Under this regulation, SOC of each energy storage unit will gradually tend to be consistent. At the same time, High-Power Battery quickly responds to sudden load change, then gradually falls back, and gradually transfers the power shortage to the High-Energy Battery. In other words, frequency distribution is realized in the HESS, and the output/input current of the High-Energy Battery is smoothed, which effectively extends the service life of the High-Energy Battery and improves the economy and reliability of the energy storage system.

Random Fluctuations in Renewable Energy Generation
Considering the random fluctuation of renewable energy, the response characteristics of High-Energy Battery and High-Power Battery in the control strategy proposed in this paper are studied. The initial state of each energy storage unit is consistent with the previous section. Figure 12 shows the output curve of renewable energy adopted in simulation, the power response and SOC changing of High-Energy Batteries and High-Power Batteries.
From Figures 12B,C, it can be seen that multi-HESS can effectively make real-time power compensation in the case of random fluctuation of renewable energy. The partial power compensated by the High-Energy Batteries is more gradual than that of the High-Power Batteries and achieves the expected frequency distribution effect.
At the same time, from Figures 12D,E, the output of different High-Energy Battery and High-Power Battery varies due to SOC differences and the energy storage output power with high SOC is larger, so that SOC tends to be consistent. Moreover, due to the setting of free charging/discharging zone, the regulation of energy storage unit is smoother, reducing the number of charging/ discharging mode changes, and effectively extending the service life of energy storage unit. Simulation results show that the proposed control strategy can ensure the

Failure of a High-Power Battery
It is assumed that a High-Power Battery breaks down and runs out, and all energy storage units are in free charging/discharging   Figure 13. Figure 13A shows that when it comes to 10 and 30 s, the cut and input of No. 1 High-Power Battery only causes a small impact on the bus voltage. Figure 13B shows that SOC of No. 1 High-Power Battery remains unchanged after failure, while the decline rate of other High-Power Battery increases. After the No. 1 High-Power Battery is repaired and put into operation in 30 s, the SOC of the three groups of High-Power Battery is still in the free charging/discharging zone, so it quickly keeps in sync. It can be seen from Figure 13C that the other two groups of High-Power Battery respond quickly at the moment of the cut and input of No. 1 High-Power Battery. Simulation results show that the remaining energy storage under the proposed control strategy can still work normally and respond to the expected power when some High-Power Batteries fail.

Failure of a High-Energy Battery
Under the same condition of section 4.2, replace the fault unit with No. 1 High-Energy Battery. The simulation results are shown in Figure 14.
The results in Figure 14 are similar to Figure 13, but more volatile. The main reason is that the High-Energy Battery has to bear more power in normal operation, and the impact on the system during cutting and input is also greater. But bus voltage deviation is still within a reasonable range, it's not more than 5% and soon returns to the normal reference voltage. The proposed method is robust to partial energy storage element faults.

Experimental Verification and Result Analysis
A microgrid model was established in the laboratory for verification. The analog diesel generator cabinet and analog wind generator cabinet are used as the power supply, the switching frequency of the energy storage bidirectional DC/DC converter is 10 kHz, the load is a programmable DC load, HEB is lithium iron phosphate battery, HPB is supercapacitor and other control parameters are the same as the simulation model. The experimental results are shown in Figures 15, 16.
It can be seen from Figures 15, 16A and Figures 15, 16B that after power fluctuation, the bus voltage can be quickly restored to stability with small voltage fluctuation. Figures 15,  16E and Figures 15, 16F show the output power changes of HEB and HPB. As can be seen from the figure, the change of HEB is relatively slow, while the response of HPB to power fluctuation is fast, which is in line with the design goal of the control strategy. Combined with the SOC change curves in Figures 15, 16C and Figures 15, 16D, it can be seen that the power curve of energy storage when SOC is higher than the limited range moves upward, that is, the discharge power is higher while the charge power is lower. When SOC is below the limit range, the power curve of energy storage moves downward, that is, the discharge power is lower while the charge power is higher. This shows that the proposed control strategy can effectively carry out SOC balancing.

CONCLUSION
Aiming at the problem that the traditional distributed multi-HESS cooperative control method of DC microgrid relies heavily on external communication, a communication free dynamic cooperative control strategy based on SOC is proposed. The simulation and experiment results show that the proposed control strategy has the following advantages: 1) being able to adjust the power distribution within HES module and among HES modules meet the needs of large-scale distributed new energy generating sets and loads; 2) using HPB to respond to the high frequency component of the bus voltage fluctuation, and using HEB to compensate the insufficient power of HPB to realize the optimization of HEB charge-discharge current, so as to effectively reduce the number of HEB charge-discharge mode change and prolong its service life; 3) reducing the dependency on communication among HES modules, thus improving the Frontiers in Energy Research | www.frontiersin.org January 2022 | Volume 9 | Article 795513 reliability of DC microgrid in island mode. However, this method does not take into account the impact of energy storage element aging, such as SOC estimation errors and reduced maximum output power due to declining energy storage life. How to consider these aspects under the condition of minimal communication requirements will be the focus of future research.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
HL completed the theoretical derivation, modeling and simulation, data analysis and full text drafting of the project. YZ and YX provided research direction, theoretical guidance and paper revision. LF provided theoretical guidance, financial and site support.

FUNDING
This work was supported by the National Natural Science Foundation of China under Grant 51877211.