A Multiport Electric Energy Routing Scheme Applied to Battery Energy Storage System

In this paper, the research status of topology and control strategy of energy storage grid-connected system is analyzed, and aiming at the working characteristics of the repurposed battery, a cascade power electronic transformer (CPET) with independent DC output is proposed. The working principle of current fed isolated bidirectional DC-DC converter (CF-IBDC) and cascaded H-bridge (CHB) is analyzed, and the decoupling control strategy is designed. In this paper, a hierarchical control strategy is designed for the repurposed battery energy storage (RBES) grid-connected system based on CPET, which consists of three layers: energy layer, power layer, and state of charge (SOC) layer. The energy layer responds to active and reactive power scheduling instructions, the power layer controls the grid-connected current and tracks the grid voltage, and the SOC layer equates the charged state of repurposed batteries. A 3 MVA/12 kV three-phase grid-connected simulation system was established, and a 1 kW single-phase system experiment platform was designed. The simulation and experimental results can verify the correctness of the theoretical analysis and the feasibility of the control strategy.


Introduction
In recent years, the proportion of renewable energy in the power system has gradually increased, but its output power is characterized by volatility and intermittency, which limits the capacity of renewable energy generation to connect to the grid on a large scale [1,2]. Supporting battery energy storage system can effectively improve the ability of power grid to accept renewable energy [3][4][5][6]. e cost factors of large-capacity converters and energy storage batteries limit the promotion and application of the battery energy storage system, while repurposed lithium batteries and link topology provide a new idea to solve the cost problem [7]. e echelon utilization of the retired power battery of electric vehicles is to design the repurposed battery cells into a battery pack after screening. Due to the limitation of the production process of the battery itself, the actual capacity, internal resistance, voltage, and other performance parameters of the battery are different to some extent, which is collectively referred to as the inconsistency of the battery [8][9][10]. After long-term use of repurposed battery packs, the battery performance degrades, and the single battery parameters will show obvious dispersion, which is manifested as obvious capacity difference, internal resistance difference, and voltage difference between batteries [11][12][13]. e existence of inconsistencies in battery cells leads to uneven charging and discharging within the battery pack. If not controlled, longterm operation will greatly reduce the reliability and security of the battery energy storage system [14]. In the design of large-capacity energy storage system, the battery pack needs to be continuously charged and discharged, and the inconsistency of discreteness will lead to a large difference in the decay speed of battery parameters, thus accelerating and amplifying the inconsistency of the battery pack [15,16]. In the echelon utilization of repurposed batteries, the inconsistency problem will be solved from two aspects. First, we select batteries with small differences to form the battery pack and build the battery management system to reduce the inconsistency of batteries within the same group. Second, we select the appropriate electrical topology structure and design the balanced control, which can slow down the rate of battery performance decline.
PET integrates voltage transformation, high-frequency electrical isolation, and flexible power flow control, making it a research direction of topological structure of battery energy storage system. For the PET used in distributed energy storage methods, the control targets are mostly DC link voltage equalizing control and isolation level parallel submodule power equalizing control. Reference [17] presented cell and phase balance for power and energy of a 10-cell multilevel inverter cascaded H-bridge. e CHB inverter control for each stage leg used the type of phase-shifted carrier. e control strategy based on the idea of a virtual synchronous generator is used on the way to balance the power and energy distributed among cells in the same phase and between phases. Reference [18] described a 14.14 kV, 2 MW, and 1000 Ah system with fortyfive Li-(lithium-) ion battery units. e system is based on a cascaded H-bridge multilevel PWM converter with star configuration focusing on the active power balancing of individual converter cells and a control system is proposed which consists of active power control and SOC balancing control (interphase SOC balancing control and phase-phase SOC balancing control). In [19], two SOC balancing techniques are proposed for an electrical vehicle charging station which is based on a grid-tied cascaded CHB multilevel converter. e first proposed technique uses the redundant states of the CHB converter to generate different AC voltages to balance the SOCs of the CHB cells. In the second proposed technique, the information of AC input current is employed to design the switching states at each quarter of the period. Hu et al. [20] proposed a novel hybrid power converter topology based on CHB and matrix converter. e CHB, as the energy control unit, can be used to balance the voltage/SOC of battery packages individually and significantly improve the modularity, stability, and safety of EV system. Kandasamy et al. [21] investigated a module level SOC balancing control in cascaded CHB based BESS using multidimensional modulation technique. SOC can be equalized between the battery modules in a phase leg of CHB with the proposed control strategy.
Aiming at the inconsistency problem of repurposed batteries, this paper presents a cascade power electronic transformer with independent DC outputs. e RBES combines the control target and topology structure to design the three-layer control architecture including energy layer, power layer, and SOC layer. e energy layer is used to respond to scheduling instructions and suppress power fluctuations, peak load clipping, and reactive power compensation. e power layer distributes active power based on SOC and divides reactive power equally to make full use of system capacity. SOC layer equalizes the charged state of decommissioned lithium battery to avoid short plate effect and stabilize the voltage of high voltage bus. e security and reliability of the RBES are enhanced.

System Composition.
e echelon RBES grid-connected system of the CPET proposed in this paper is shown in Figure 1. e system adopts a two-stage structure, with an independent DC output terminal to access the repurposed battery, and the topology is divided into the DC cell group consisting of cascade CF-IBDC module and the DC-AC consisting of CHB module. Based on a 3 MVA/12 kV medium voltage grid system, considering the trade-off among system cost, life, passive components, switching devices, frequency, and power quality during design progress, the number of each phase containing H-bridge cell is N � 4, and the number of CF-IBDC cell modules contained in the DC group corresponding to each H-bridge is M � 4. Within the DC-DC unit, the high-voltage side of the CF-IBDC adopts series structure to obtain the high-voltage bus voltage, V dcia � 3000 V, and the low-voltage side of the converter is connected to the repurposed battery module, respectively, to realize independent energy balance control. e H-bridge module constitutes a multilevel cascade system, realizing the grid-connected control.

Cascade H-Bridge Analysis.
e CHB topology is shown in Figure 1. All modules are in series structure; that is, the output current of the modules is the same, and the power distribution between the modules will be determined by their output voltage. In order to realize the independent control of active and reactive power, a discrete Fourier transform phase-locked loops (PLL) method is adopted. is method is only based on single-phase grid voltage signal and can extract basic phase, frequency, and amplitude information from any signal. e active power and reactive power can be decoupled by coordinate transformation, and then they can be controlled independently. e distribution of active power and reactive power of phase A among N H-bridge modules is illustrated in Figure 2. e same analysis is also carried out for phase B and phase C, which will not be repeated in this paper.
In the coordinate analysis, considering the relative stability of the grid voltage, the synchronous signal is V ga . e α-axis in the αβ coordinate system is in phase with the grid voltage, and the β-axis lags the α-axis by 90°. e d-axis in the dq coordinate system is aligned with the grid voltage through PLL control technology, and q-axis lags d axis by 90°, as shown in Figure 2(a). e grid voltage components in the αβ coordinate system and the dq coordinate system are given by the following equations, respectively: where ω is the grid voltage frequency, V ga is the amplitude of phase-A grid voltage, and v ga d � V ga ; V ga_q � 0.

Journal of Engineering
A vector diagram of the AC voltage is given in Figure 2(a) to illustrate the power distribution principle between the CHB modules of phase A. e same analytical method can be applied to phase B and phase C. e d′-q′ coordinate system is constructed in the vector diagram, where the d′ axis is in phase with the grid current, and the q′ axis lags the d′ axis by 90°. Obviously, the d′ axis component of the inverter output voltage V sa d′ determines the active power, while the q′ axis component V sa q′ determines the reactive power. In Figure 2(b), the power distribution of CHB topology in different grid-connected power situations is described. e grid-connected voltage V sa output by CHB is synthesized from the output voltage V ia (i    Figure 2: e relationship between the α-β coordinate system, the d-q coordinate system, and the d′-q′ coordinate system. (a) e relationship between phase A grid voltage V ga , grid current I ga , and inverter output voltage V sa . (b) H-bridge inverter output voltage distribution.
V ia d′ and V ia q′ can be decoupled to realize independent control of active power and reactive power of each H-bridge module. In the coordinate transformation of d-q coordinate system and d′-q′ coordinate system, the grid voltage and current reference angle θ ga are the key to the transformation. e grid current i ga can be measured and projected onto the coordinate system to obtain i ga_α and i ga_β . e grid current flowing through αβ-dq coordinate system is transformed into the dq coordinate system as follows: erefore, the reference angle of θ ga of grid voltage and current can be obtained: We substitute the result obtained from (4) into dq-d′q′: e voltage component based on d′q′ coordinate system can be obtained by (5), the voltage signals u sa d′ and u sa q′ of active power and reactive power output are distributed to each H-bridge module through the energy balance control algorithm, and then the d′ axis voltage reference value of each H-bridge module can be obtained, so as to control the active power component of corresponding module, and the obtained q′ axis voltage reference value controls the reactive power component. In order to convert the voltage component of the d′q′ coordinate system into the modulation ratio control quantity of each module of the CHB topology, we transform the output voltage reference signal of the active power and reactive power of the d′q′ coordinate system obtained by the equalization control to the αβ coordinate system through the d′q′-αβ coordinate system: As shown in Figure 2(a), the output reference signals of each module are u * ia � u * ia α , the modulation ratio of each H-bridge module can be calculated: m * ia � u * ia /V dcia . en, the decoupling control of active power and reactive power can be completed, and the power balance control of each H-bridge module can be realized easily. In the control, the modulation ratio should not exceed the modulation ratio control limit to avoid output instability. is paper does not analyze the overmodulation control.

Cascading CF-IBDC Analysis.
e structure of CF-IBDC is shown in Figure 3. e low-voltage side of the converter is a multiplexed half-bridge structure. e low-voltage halfbridge and the input inductor L b constitute the front Buck/ Boost circuit to obtain the capability of wide input voltage range.
e input inductance makes the converter have current source characteristics to reduce the input current ripple. At the same time, the half-bridges and capacitorbridges in both primary and secondary sides constitute the dual active bridge. e capacitor arm clamping the voltage of switching tube can solve the problem of high voltage spike. e center tap on the high voltage side of the transformer yields two high voltage half-bridge structures, and interbridge circulation can solve the soft switch problem under light load and improve the efficiency in the full power range.

Analysis of Hierarchical Control Strategy
e CPET based RBES grid-connected system designed in this paper mainly realizes the control target of smooth power fluctuation and peak cut of renewable energy. e system is divided into 2-layer structure, composed of DC-DC part and DC-AC part, which can realize the independent control of the two parts. Combined with the grid-connected RBES topology, the control requirements and objectives of the energy storage system are decomposed, and the three-layer control structure of the RBES system with CPET structure is designed to achieve their respective control objectives, namely, the energy layer, the power layer, and the SOC layer.
In the RBES grid-connected system, the DC-DC stage is divided into charging state and discharging state, and the DC-AC stage is divided into inverting state and rectifying state. When the RBES grid-connected system works in the inverter state, the repurposed battery discharges and the energy is transmitted forward. e whole system is powered by the repurposed battery to the grid. When the RBES gridconnected system works in the rectifier state, the repurposed battery is charged and the energy is transmitted in reverse, and the whole system absorbs the excess energy from the grid side to the repurposed battery.
Whether the energy is transmitted forward or backward, the control objectives of the DC-AC stage are grid-connected current tracking of the grid voltage and the completion of active and reactive power scheduling instructions, while the control objectives of the DC-DC stage are to stabilize the voltage of the HVDC bus and balance the charged state of repurposed batteries. e DC-DC stage and the DC-AC stage are connected by the HVDC bus. e voltage of the HVDC bus V dcik will be affected by the two stages of the circuit. In the case of rectifier, the DC-AC stage raises the voltage of the DC bus, while in the case of inverter, the DC-AC stage pulls down the voltage of the DC bus. In this paper, the HVDC bus voltage is controlled by the DC/DC stage. Based on the above analysis, it can be known that the automatic switching between forward flow and reverse flow of energy can also be realized in the stable voltage control implementation of CF-IBDC converter; two-stage circuit control can also be completely decoupled which only focus on their respective control objectives. Hierarchical control strategy diagram is shown in Figure 4.

Energy Layer Control Strategy.
According to the power output of renewable energy and the expected design value of system power at the current time of the power grid, the power dispatching center obtains the active power dispatching instruction P g and reactive power dispatching instruction Q g and sends the dispatching instructions to the RBES grid-connected system. e RBES grid-connected system achieves the control target of energy layer according to the control strategy shown in Figure 5.
In the energy layer control, the RBES system detects the reactive power QRBES and active power PRBES of PCC nodes and takes them as a feedback value. e reactive power dispatch instruction Q g and the active power dispatch instruction P g are, respectively, differentiated from the feedback value, getting the next cycle adjustment of power value ΔP. en, the power closed-loop control is completed to obtain the reference value of active power and the reference value of reactive power. By transferring the reference value to the power layer, the RBES grid-connected system based on CPET is controlled to realize the tracking control of scheduling instructions. e energy layer regulates the dynamic response speed of the RBES grid-connected system, controls the start-stop and power switching of the system, which realizes the nonimpact and flexible switching of operating conditions, and controls the active power and reactive power output of the RBES system to reach the grid dispatching value.

Power Layer Control Strategy.
e control objective of CHB topology is to track the active power of the energy layer by referring to the given value P * ref and the reactive power by referring to the given value Q * ref , which realizes the dispatching instruction of the grid dispatching center. e gridconnected current tracks the voltage on the side of the grid which is the standard sine wave. e control strategy of the power layer is shown in Figure 5. e above energy layer control strategy calculates the given values of active power and reactive power and transfers the reference values to the power layer. Firstly, according to (7), the power layer control obtains the given reference quantity of d axis current and the given reference quantity of q axis current. e active and reactive power value of the  Figure 3: CF-IBDC topology.  Journal of Engineering RBES grid-connected system is evenly divided into the three-phase system. According to single-phase PLL dq vector transformation technology, the voltage current detection signal of the system is obtained by abc/dq vector transformation in the dq static coordinate system, and the active power and reactive power are decoupled-controlled in the static coordinate system. e system error signal is obtained by the difference between the given reference value of dq axis current and the actual current flow, which is calculated by PI controller to obtain the voltage signals of d axis and q axis: e voltage signal of dq axis and the charged state SOC ik value of each DC unit group are input into the power distribution module, respectively. According to the active power and reactive power distribution (8), the voltage given signal of dq axis of each module can be obtained. Although the active power according to SOC ik value distribution can inhibit DC unit charged state significantly difference between groups and avoid the short board effect, it actually cannot completely eliminate the differences of DC units charged state between groups. Considering the capacity and control complexity of system, the small difference of charged state in engineering application is acceptable, so the control strategy of allocating active power according to the charged state of each DC unit group is feasible. e reactive power is an equilibrium control, and each H-bridge module has the same parameters and reactive power capacity, which makes the reactive power divided equally among all H-bridge units: e modulation ratio of each module is different, which makes the RMS value of output AC voltage inconsistent, while the link structure makes the RMS value of each module the same, so each module has different power output value, which indirectly realizes the SOC value equalization control, and controls the total power of the cascade module and the change of the scheduling value of the system.

SOC Layer Control Strategy.
In the SOC layer, the charging and discharging power of each module is evenly controlled according to the SOC status of the repurposed battery, which enables the balanced control of the SOC value of the repurposed battery to ensure the safe and stable operation of the battery pack. Combined with the working principle of the CF-IBDC converter, it can be seen that the power transmission direction is controlled by the phase shift angle between bridges, while the energy flow direction is automatically switched under the voltage stabilization control of the HVDC bus. Since the cascaded structure of each module makes the module have the same current, it is necessary to control different output voltage to regulate the output power of each module, so SOC equalization is to control the output voltage value of each module.
In Figure 6, voltage loop is used in CF-IBDC to realize the output voltage control, the control parameter is the phase shift angle ϕ between bridges, and voltage matching is achieved by controlling duty cycle D, which can reduce converter loss. SOC equalization control strategy requires that different voltage values of each module be given, respectively, and output power be adjusted to realize SOC equalization. e voltage error signal can be obtained by making the difference between the voltage given value and voltage sampling V H . e phase shift angle can be obtained by PI regulator. e PWM driving signal can be generated by the phase shift square wave controller, and then the voltage stabilizing control can be completed by adjusting the CF-IBDC. e voltage distribution controller is shown in Figure 7, the difference between the average SOC ik value of the DC unit group and the state of charge value of each module to obtain the error signal of the SOC value, which is corrected by the PI controller for the corresponding module voltage. In the SOC equalization control, a correction coefficient calculator is designed to distinguish the different equilibrium objectives of different SOC based on charge and discharge, which is actually a hysteresis control strategy. We multiply the obtained voltage correction value and correction coefficient to obtain the SOC equalization correction voltage increment and then add the given value of the mean voltage and the SOC equalization correction voltage increment to obtain the given value of the voltage of each module. e corresponding last module compensates the voltage regulator value of the whole DC unit group; that is, it satisfies equation (9) and controls the whole DC unit group by the voltage given value of the module-M, which realizes the voltage regulator control of the HVDC bus: e correction factor calculator is shown in Figure 8. In the process of voltage rise, when the voltage of HVDC bus reaches 3050 V, the system enters the charging control area and the correction coefficient gradually increases to 1. When it is higher than 3100 V, the correction coefficient remains at 1. In the process of voltage reduction, when the voltage is reduced to 2950 V, the correction factor gradually decreases to −1 and maintains −1 after 2900 V. e voltage of the HVDC bus ranges from 2950 V to 3050 V, which is the hysteresis area controlled by the system. In other words, the voltage of the HVDC bus will fluctuate in the voltage hysteresis area and stabilize at a certain voltage value.

Simulation Results
In order to verify the effectiveness of the CPET topology and its control strategy involved in this paper, a simulation model of the RBES system was built by Matlab software. e parameters used in the simulation are shown in Table 1.
In the power layer simulation, the variation of power scheduling value is shown in Figure 9. e whole schedule consists of six steady-state operating points and dynamic processes between them, which can reflect different operating conditions of the system, respectively. As can be seen in this figure, the grid voltage V ga of phase A is taken as the reference signal, the three-phase grid current can ensure tracking the grid  Figure 5: Power layer control strategy block diagram.
Phase shi square wave generator Figure 6: SOC layer control strategy block diagram.
Journal of Engineering voltage as a standard sine wave after the change of active and reactive power, and the RBES system can output the dispatching power value with no difference. In the start-up and power switching state of the RBES system, the impact of the power grid is small, and the scheduling instructions can be quickly tracked and adjusted in a short time to meet the engineering requirements of the large-capacity system. In reactive or active scheduling instruction switching, the two can reach a stable state without mutual influence. Figure 10 shows the modulation waveforms and the output voltage waveforms of each module when the DC unit group has different SOC values. In the simulation setting, the power of the energy storage system is 1.5 MW/1.5 MVar, and the SOC values of the DC unit group are 0.7, 0.7, 0.8, and 0.75 in turn.
It can be seen that each modulated wave has the same phase but different amplitude. In the detail figure, SOC values of H 1a module and H 2a module are the same, and the modulation wave is the same; the SOC value of H 3a module and H 4a module are both larger than that of H 1a module, so it can be seen that the amplitude of modulation is greater than that of H 1a module, and the difference of amplitude of modulation is equal. e corresponding output waveforms of each module are ±3 kV and 0 V level, while within a cycle, the phase shift between the drive signals of adjacent modules is 45°, and different modulation amplitude corresponds to the conduction time of different switching tube.
In the SOC layer simulation of the RBES system, the parameters of capacity and SOC value of repurposed lithium batteries within the DC unit group are shown in Table 2. e entire charging and discharging processes are simulated to verify the correctness of the SOC balancing strategy.
From the comparison and analysis of the SOC value in the charging and discharging simulation waveform diagrams in Figure 11, it can be seen that for repurposed lithium batteries with the same initial SOC value and actual capacity, the output voltage setting values of the corresponding modules are always the same (e.g., Ei k_1 and Ei k_2 ). For repurposed lithium batteries with equal initial SOC but different actual capacity, with the increase of discharge time, the difference of output voltage setting values of modules gradually increases, but the SOC values remain the same (e.g., Ei k_1 and Ei k_3 ). For repurposed lithium batteries with different initial SOC values and the same actual capacity, with the increase of discharge time, the output voltage setting values corresponding to the modules gradually tend to be consistent and finally equal, and the difference in SOC values also gradually decreases, tending to be consistent and finally equal (e.g., Ei k_1 and Ei k_4 ). Under different working conditions of the RBES system, the voltage of high-voltage bus can always maintain a steady state of 3 kV, which verifies the effectiveness of the control.

Experimental Results
In order to test the control algorithm, a single-phase RBES grid-connected experimental system with rated power of 1 kW was designed by referring to the parameter design and device selection method of the aforementioned 3 MVA/ 12 kV system, so as to verify the correctness of the conclusions obtained from the above analysis. e experimental platform is designed as a single-phase AC 220 V grid-connected system, in which the number of cascaded H-bridges is N � 2, and the number of cascaded CF-IBDC is set as M � 2.  Journal of Engineering e design of the experimental platform is shown in Figure 12.

Correction factor calculation
e rated power of the single CF-IBDC converter is 250 W, the rated input voltage is 50 V, the output voltage is controlled at 100 V, and the two modules cascade to obtain 200 V. e rated power of a single H-bridge module is 500 W, the rated input voltage is 200 V, and the detailed parameters of the CF-IBDC converter and CHB topology are shown in Table 3. Figure 13 shows the experimental platform diagram of the DC unit group. In the experiment, two DC stabilized power supplies are connected to the CF-IBDC converter    Figure 14 shows the soft switching waveforms of the CF-IBDC converter under the no-load condition. In the no-load state, the phase shift angle ϕ between the bridges is close to 0°, and the converter generates a circulating current through the phase shift angle ϕ s inside the bridge to achieve ZVS. Figure 15 shows the soft switching waveforms at half load. When the load reaches a certain condition, the output current can realize the soft switching of the high-voltage side switching tubes. In order to improve the overall efficiency of the converter, the SPS phase shift method is adopted. e phase angle ϕ s is always 0, and the converter can realize ZVS depending on the operating current. Figure 16(a) shows the steady-state output waveforms at half load, the output voltage of each module has slight fluctuations, and the peak-to-peak output voltage is less than 5 V, which meets the system requirements. e input voltage of the module D 1a_1 is 45 V and the output voltage obtained by the equalization calculation is 95 V. e input voltage of the module D 1a_2 is 50 V, and the output voltage obtained by the equalization calculation is 105 V. e measured output voltage can be seen that the steady-state output voltage of module D 1a_1 is 95.1 V, the steady-state output voltage of module D 1a_2 is 104 V, the cascaded output voltage is 200 V, and the output power is 250 W. Figure 16(b) shows the dynamic waveforms switching from half load to 3/4 load. e steady-state output voltage of module D 1a_1 is 95.1 V, the steady-state output voltage of module D 1a_2 is 104 V, and the steady-state cascaded output voltage is 200 V. When the load is switched, the output voltage of the module D 1a_2 fluctuates, which further affects the cascade output voltage, but the output gets the stability soon. It can be seen that the output voltage stabilizes before and after load switching of each module, which can realize different power output of each module under the premise of stable output voltage on the cascaded side. Figure 17 shows the experimental waveforms when the output power of the CHB is 1.1 kW, and the current waveform is in phase with the voltage waveform. e effective value of the voltage is 220 V, the effective value of the current is 5 A, the frequency of voltage and current is 50 Hz, and the voltage and current waveforms are present. It is a standard sine wave with no phase difference.    Figure 18 shows the waveforms of the power switching. e CHB output voltage v sa is not affected by power switching. e AC current i ga has almost no distortion, and the voltage and current waveforms have no phase difference. Figure 19 shows the waveforms of the off-grid operation with load and the grid-connected operation with resistive load. After grid-connected, the CHB converter delivers energy to both the off-grid load and the grid.     Figure 15: ZVS waveforms at half load.
By comparing the experimental waveforms and simulation waveforms, the correctness and effectiveness of the theory can be verified. e whole system can realize active and reactive power grid connection in the cascade state of modules.

Conclusions
Based on the research background of cascade utilization of repurposed batteries in power energy storage system, a series of analysis and research on RBES system based on CPET were carried out in this paper. e topological structure and control strategy of DC-AC stage and DC-DC stage in this system were analyzed.
CPET is composed of current fed dual active bridge modules and H-bridge module. AC side can be directly connected to the medium and high voltage power grid, and low-voltage side has an independent DC output end. e hierarchical control strategy of the RBES system based on CPET is composed of the energy layer, the power layer, and the SOC layer. e energy layer responds to the active and reactive scheduling instructions, the power layer controls the grid-connected current and tracks the grid voltage, and the SOC layer equates the charged state of repurposed batteries.

Data Availability
e data used to support the findings of this study were supplied by Jiaxun Teng under license and so cannot be made freely available. Requests for access to these data should be made to Jiaxun Teng (tengjiaxun@qq.com).

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this study.