Multi-Time Scale Economic Optimization Dispatch of the Park Integrated Energy System

The park integrated energy system (PIES) plays an important role in realizing sustainable energy development and carbon neutral. Furthermore, its optimization dispatch can improve the energy utilization efficiency and reduce energy systems operation cost. However, the randomness and volatility of renewable energy and the instability of load all bring challenges to its optimal operation. An optimal dispatch framework of PIES is proposed, which constructs the operation models under three different time scales, including day-ahead, intra-day and real-time. Demand response is also divided into three levels considering its response characteristics and cost composition under different time scales. The example analysis shows that the multi-time scale optimization dispatch model can not only meet the supply and demand balance of PIES, diminish the fluctuation of renewable energy and flatten load curves, but also reduce the operation cost and improve the reliability of energy systems.


INTRODUCTION
The park integrated energy system (PIES) is one of the main forms of distributed renewable energy utilization, which plays an important role in promoting the absorption of renewable energy and establishing a new power system under lower-carbon (Zahedi, 1994;Tang et al., 2018). Due to the fluctuation of renewable energy and the increase of load prediction error caused by the increase of forecast time scale, it is difficult to satisfy the balance of supply and demand through energy systems dispatch under a single time scale. Thus, the multi-time scale optimization dispatch of PIES is proposed. The dispatch plan is updated and adjusted according to the coordination of energy equipment and the use of more accurate forecast information of renewable energy and load, so as to ensure the economic operation, reduce the volatility of renewable energy and load, and improve the feasibility of the dispatch plan (Fan et al., 2021).
Many progresses have been made in the research on the coordinated planning of multiple types energy equipment in PIES and the improvement of overall operation economy. Literature (Zhang et al., 2019;Hui and Wei, 2020) constructed an integrated electricity-gas-thermal system planning model considering demand response, energy storage and energy conversion. Literature Jiang and Ai (2019) put forward an integrated energy system optimization operation model considering the new energy incentive mechanism to clarify the influence of different incentive mechanisms on the economic benefits. According to the quasi-dynamic characteristics of the thermal system, Literature Cen et al. (2020) established a building virtual energy storage model and an integrated energy utilization optimization model under the TOU price. Literature Lin and Fang (2019) established a regional integrated energy management system model based on multi-agent method to optimize operator profit and user cost. Literature Zeng et al. (2017) introduced a control optimization strategy considering economic performance and response performance to improve the multi-energy optimization of integrated energy system. Literature Zhang and Wei (2020) proposed a two-layer capacity allocation optimization method for IES based on the combined gas cooling, heating and power supply scheme.
However, the above literatures mainly studied the principle of economic optimal operation and the framework of system model. The integrated demand response is not considered or just regarded as a constraint condition of the model and the operation method and response strategy of the integrated demand response are lack of in-depth exploration. In Literature (Pu et al., 2020), an optimal operation method of power system considering power demand response was proposed to transform the energy supply mode of port integrated energy system. Literature Lv et al. (2021) established an optimized operation model of IES and analyzed the flexibility of demand response under different coupling degrees. Literature Zhao et al. (2019a) introduced an electric-gas-thermal system and proposed a rolling operation optimization strategy considering the multienergy demand response and wind penetration. Literature Qi et al. (2020) proposed a two-layer collaborative optimization method of IES considering DR and ES. In Literature , a multi-energy flow calculation model and a dayahead optimal scheduling model for the power-gas-heat-cooling integrated energy system were established considering the integrated demand-side response, wind power output prediction and energy storage. Literature Wang C. et al. (2020) established a demand response model based on compensation mechanism to analyze the scheduling elasticity of different load types. Literature Wang J. et al. (2020) proposed a two-stage optimal scheduling model for integrated energy system from the perspective of user demand-side load control.
The time scale characteristics of each energy flow in the comprehensive energy are different. In terms of operation dispatch of integrated demand response, it is necessary to discuss the response speed of different energy flows to price change and external demand, analyze the influence of start-stop time and climbing speed of equipment on the complementarity of energy sources, study the time requirements of different energy flows on the balance between supply and demand (Duan et al., 2017). The response strategy and cost of load under different time scales also need to be studied. In Literature (Yang et al., 2020), a day-ahead and intra-day optimal scheduling model for integrated energy systems considering demand side response was established to suppress the intra-day fluctuations of renewable energy and loads. Literature Ouyang et al. (2019) constructed a multi-time scale optimization model with three different time response characteristics of electricity, heat and cold to optimize the economic operation of energy systems. Literature Xue et al. (2020) presented a real-time energy management method for integrated heating and power system based on approximate dynamic programming algorithm. Although the above studies introduced integrated demand response to improve the flexibility of load dispatch and increase the economic benefits of energy systems, they did not consider the uncertainty of load response at a small time scale and the impact on the operation dispatch of energy systems, so they could not meet the real-time adjustment requirements of energy systems under a small time scale.
To sum up, this paper shows the comparisons of related works in Table 1, and summarizes the advantages and disadvantages of several significant papers in Table 2. Hence this paper proposes an economic optimization dispatch method for PIES based on multi-time scale demand response. Aiming at minimizing the operation cost of PIES, this paper constructs an optimization dispatch model under three-time scales, including day-ahead, intra-day and real-time. Based on the above model, this paper proposes a dispatch strategy and uses MATLAB to solve the mixed integer linear programming problem. Finally, the example analysis verifies the effectiveness and feasibility of the proposed method.
The key contributions of this paper are listed as follows.
1) Constructing a three-layer dispatch framework of PIES considering demand response, including the day-ahead, intra-day and real-time optimal dispatch. 2) Proposing the multi-time scale economic optimization strategy of PIES based on the multi-energy coordination and the effect of the fluctuation of renewable energy and load, including giving the objective functions, constraints and corresponding solution strategies in consideration of demand response. 3) Conducting case studies under a PIES scenario to verify the effectiveness of the proposed model and strategy, and analyzing the dispatch scheme and economic benefits of PIES under different time scales.
The organization of this study is as follows. PIES Dispatch Model introduces the structure and dispatch framework of PIES, and presents the dispatch model. Economic Optimal Dispatch of PIES Considering Demand Response proposes the dispatch strategies of PIES in consideration of demand response and solution scheme. Case Study describes the feasibility of PIES multi-time scale optimal dispatch based on the analysis of the calculation results. Conclusion summarizes the contributions and extensions of this paper.

Structure Diagram of PIES
According to the geographical factors and the characteristics of energy generation/transmission/distribution/use, the integrated energy system can be divided into trans-regional level, regional level, and user level, while PIES is the most typical representative of the regional integrated energy system (RIES). PIES takes active distribution network, hybrid energy storage and energy conversion technologies as the core, including intelligent distribution system, heating/cooling system, natural gas system and other multi-energy coupling networks . Compared with the traditional integrated energy system, PIES has a shorter energy transmission distance, more accurate and fast energy system operation, and a higher degree of coupling between energy systems.
The structure of PIES is shown in Figure 1. Energy systems are coupled by four energy flows: electricity, gas, heating and cooling, and them contain a variety of energy production, conversion and storage equipment to meet different loads. However, only CCHP is selected as the energy coupling equipment in the system structure, which consists of gas turbine, waste heat boiler, absorption refrigerator, and electricity chiller. On one hand, for general energy systems, CCHP can not only undertake the full function of energy conversion , but also is the most effective unit to improve the economic benefits (Ming et al., 2015;Zeng et al., 2017). On the other hand, this measure  - can simplify the calculation, which can avoid discussing the calculation of the start-stop state of other energy conversion equipment and the impact of the start-stop cost on the economy.
Since there are many studies on the steady state mathematical model of CCHP, this paper will not elaborate on it (Cheng et al., 2016;Zhao et al., 2018).

Dispatch Framework of PIES
In the optimization operation of PIES, day-ahead dispatch needs to formulate the operation plan of energy systems 24 h in advance, while intra-day and real-time dispatch should be scheduled according to energy production and coupling equipment output at various periods. Firstly, a day-ahead optimization dispatch can realize the optimal configuration under a long time scale and improve the economy of system operation. However, considering that the forecasting errors caused by renewable energy and load gradually decrease with the approaching of running time point (Drauz et al., 2018), dayahead dispatch cannot meet the actual operation requirements. Meanwhile, intra-day dispatch and real-time dispatch can make use of the forecasting information of load for dispatch under a short time scale, but they cannot consider the global coordination and collocation in the operation plan of energy systems. On the other hand, integrated demand response faces the following problems in the multi-time scale optimization dispatch (Duan et al., 2017). For one thing, from the perspective of time scales, it is necessary to adapt to the prediction error caused by different time dimensions, the adjustment of resource scope and the difference of cost. For another, it is essential to adapt to the difference of demand response and loads with different characteristics when they participate in the power grid dispatching operation.
Based on the above problems, this paper proposes a combinatorial optimization dispatch framework for PIES, which consists of day-ahead optimal dispatch model, intra-day optimal dispatch model and real-time dispatch model, as shown in Figure 2. First, the day-ahead optimization dispatch plan is formulated every 24 h and the time scale is set as 1 h. The interactive power between energy systems and power grid, the output of dispatchable energy equipment and flexible load response can be determined according to the short-term forecast output of day-ahead load and the predicted power of photovoltaic and wind turbines (Zhao et al., 2019b). Second, on the basis of day-ahead optimal dispatch, the intra-day optimization dispatch plan is made every 1 h to adjust the output of multiple equipment in the future 4 h and the time scale is set as 15 min. Third, real-time dispatch takes into account the intraday optimization dispatch and makes a plan every 15 min. It determines the next time scheduling plan according to ultrashort-term forecast load output.
There are three kinds of demand response in the optimal scheduling model of PIES, which are all incentive-oriented demand responses (IDR). According to the time of advance notification to users, the IDR is divided into three categories: 1 day in advance (A-type demand response), 4 hours in advance (B-type demand response), and 15 minutes in advance (C-type demand response).

Dispatch Model of PIES
This paper takes 24 h forecast power of load and renewable energy generation as the input of model, while dispatchable output of the combination of energy systems and equipment as the output of model, so as to construct an optimization dispatch model of PIES, which is expressed as follows: It is assumed that the dispatch plan of energy systems is divided into n time periods. In terms of input, P load is a column vector of length n that indicates forecasting load and its element P load i indicates the load forecasting power in the i-th period. P WT and P PV are column vectors indicating prediction power of wind turbine and photoelectric and their length are both n. Their elements P WT i and P PV i indicate the predicted power output of wind turbine and photovoltaic in the i-th period.
In terms of output, energy purchase of energy systems mainly includes electricity and gas, so P buy is a matrix of n × 2, whose element P buy ij indicates the purchased power of the j-th type energy in the i-th period. The output of energy storage system P SOC includes three parts: electricity, heat and cool. Therefore, P SOC is a matrix of n × 3, whose element P SOC ij indicates the storage and discharge power of the j-th type energy storage in the i-th period. Here positive values indicate charging and negative values indicate discharging. P CCHP is an n × 3 matrix indicating the power generation, heating and refrigeration of CCHP equipment (combined cooling heating and power) and its element P CCHP ij indicates the output power of the j-th type energy generation in the i-th period. The length of load response vector P DR is n and its element P DR i indicates the dispatchable load of energy systems in the i-th period under the corresponding time scale.
According to Eq. 1, the dispatch model of PIES can realize optimal dispatch through reasonable arrangement of energy purchase, energy storage equipment, CCHP equipment, and dispatch strategies of flexible loads under the condition that the energy supply demand is different between the predicted output of energy systems load and that of the renewable energy.

ECONOMIC OPTIMAL DISPATCH OF PIES CONSIDERING DEMAND RESPONSE
The Objective Function of Economic Optimal Dispatch

The Objective Function of Day-Ahead Dispatch
In the day-ahead optimization dispatch model, this paper takes the cost of energy storage and CCHP equipment operation, energy purchase and load response as decision variables to solve the minimum operating cost of energy systems in the future 24 h considering the short-term forecast load: where F 1 indicates the expected total cost of the day-ahead system; F s,x t , F CCHP t , F buy t and F sub,A t respectively indicate the operation and maintenance cost of energy storage equipment, the operation and maintenance cost of CCHP equipment, energy purchase cost and A-type demand response compensation in the t-th period; F grid t and F gas t indicate electricity purchase cost and gas purchase cost. c x,in and c x,out indicate charge and discharge maintenance price coefficients of energy storage equipment; c CCHP indicates the price coefficient of CCHP equipment operation and maintenance; c grid t and c gas indicate the electricity price purchased from power grid in the t-th period and unit natural gas price; c sub,A indicates the subsidy standard of A-type demand response. P x t indicates the charge and discharge power of the energy storage equipment in the t-th period, in which charging is positive and discharging is negative; P grid t and G gas t indicate the electricity purchase power and the gas purchase capacity in the t-th period; P CCHP t indicates the electricity power of CCHP equipment; P DR,A t indicates the regulation power of A-type demand response. η sub indicates the response coefficient of demand response. U x,in t and U x,out t are 0-1 variables, which indicate the charging and discharging states of the energy storage equipment. x indicates energy type and x is a member of the set Z, which contains three elements: E, H and C, respectively indicating electricity, heat and cool energy.

The Objective Function of Intra-day Dispatch
If intra-day dispatch plan is re-formulated in each period considering the day-ahead dispatch, there may be a large difference between them (Zhao et al., 2021), which obviously does not satisfy the actual operation requirements. Therefore, intra-day optimization dispatch adjusts the output of each equipment every 1 h on the basis of day-ahead optimization dispatch, so as to minimize the total adjustment cost of energy systems in the next 4 h: where F 2 indicates the total adjustment cost of the intra-day system; ΔF s,x 1,2 , ΔF CCHP 1,2 , ΔF buy 1,2 indicate the adjustment costs of energy storage equipment, CCHP equipment and energy purchasing plan of PIES between day-ahead and intra-day schedules in the t-th period; F sub,B t indicates B-type demand response compensation; c sub,B indicates the subsidy standard of B-type demand response; P DR,B t indicates the regulation power of B-type demand response; subscripts 1 and 2 indicate day-ahead dispatch and intra-day dispatch.

The Objective Function of Real-Time Dispatch
Real-time dispatch adjusts the output of each equipment every 15 min at the next moment considering the ultra-short-term predicted output of load on the basis of intra-day optimization dispatch, so to minimize the adjustment cost in a single period of PIES as: where F 3 indicates the total adjustment cost of the real-time system; ΔF s,x 2,3 , ΔF CCHP 2,3 and ΔF buy 2,3 indicate the adjustment costs of energy storage equipment, CCHP equipment and energy purchasing plan of PIES between intra-day and real-time schedules in the t-th period; F sub,C t indicates C-type demand response compensation; subscripts 2 and 3 indicate intra-day dispatch and real-time dispatch. where P PV t and P WT t indicate the power generation of photovoltaic and wind turbine in the t-th period; P ES t indicates the charge and discharge power of the storage device, in which charging is positive and discharging is negative; P DR,A t is positive when the load is transferred in and negative when the load is transferred out.

System Thermal Equilibrium Constraint
where H CCHP t indicates the heating power of CCHP equipment; P HS t indicates the charge and discharge power of heat storage equipment, in which charging is positive and discharging is negative; H LD t indicates the thermal load of energy systems in the t-th period; η heat indicates the ratio of heat energy utilization in heating grid.

System Cool Equilibrium Constraint
where C CCHP t indicates the refrigeration power of CCHP equipment; P CS t indicates the charge and discharge power of cool storage equipment, in which charging is positive and discharging is negative; C LD t indicates the cool load of energy systems in the t-th period; η cool indicates the ratio of cool energy utilization in cooling grid.

Energy Storage Equipment Constraint
where S x t indicates the capacity of various energy storage devices in the t-th period; S x, min and S x, max indicate the upper and lower limits of the capacity of energy storage equipment; η s,x indicates the self-consumption rate of energy storage equipment; S x 0 and S x T indicate the energy storage capacity at the beginning and the end of dispatch.

CCHP Operating Constraint
where P CCHP, min and P CCHP, max indicate the minimum and maximum generating power of CCHP equipment; H CCHP, min and H CCHP, max indicate the minimum and maximum heating power of CCHP equipment; C CCHP, min and C CCHP, max indicate the minimum and maximum refrigeration power of CCHP equipment.  (26) where μ indicates the constraint constant; subscripts 1, 2 and 3 indicate the day-ahead dispatch, intra-day dispatch and real-time dispatch. The deviation constraint is only used for intra-day and real-time dispatch.

Multi-Time Scale Economic Optimization Dispatching Strategy for PIES
In order to solve the optimal dispatch scheme of PIES, this paper schedules the cost coefficient and flexible load based on demand response under different time scales and the establishment of objective functions and constraint conditions in day-ahead, intraday and real-time dispatch. The specific solution process is shown in Figure 3.
In the above solving process, the multi-time scale economic optimal dispatch strategy of PIES constitutes the mixed integer linear programming (MILP) problem. Based on MATLAB software, this paper calls the CPLEX solver in the YALMIP installation package to solve the above MILP problem.

CASE STUDY Setting
In order to verify the effectiveness of the multi-time scale optimization dispatch strategy of PIES, two typical cases are designed to simulate the multi-constraint dynamic linear programming problem including CCHP equipment and energy storage equipment. This paper selects a PIES for analysis, and the main basic parameters are shown in Table 3. The unit price of natural gas c gas is 2.4 yuan/ m 3 , the low calorific value of natural gas is 36 MJ/ m 3 and the time-of-use electricity price is as follows: 0.821 yuan/ kWh in peak period, 0.485 yuan/ kWh in normal period and 0.368 yuan/ kWh in valley period; the operating range of the energy storage equipment is [5%,95%] and the initial capacity is 50%; the proportion of thermal energy utilization η heat is 0.9 and the proportion of cool energy utilization η cool is 0.9; the maximum capacity that can be invoked for A-type, B-type and C-type demand response shall not exceed 6% of the maximum of electric load. The response coefficient and subsidy standard of demand response are shown in Table 4.

Typical Day in Summer
This paper analyzes the daily output of each equipment on a typical day in summer. The power, heating and cooling load, and renewable energy curves of a typical day in summer are shown in Figure 4.

Day-Ahead Optimization Operation Result and Analysis
According to the above example parameters, this paper runs the day-ahead optimal dispatch model of PIES and obtains the Frontiers in Energy Research | www.frontiersin.org September 2021 | Volume 9 | Article 743619 6 FIGURE 3 | The optimal dispatch solution process of PIES. day-ahead optimal dispatch scheme on a typical day in summer, as shown in Figure 5.
It can be seen from Figure 5 that the day-ahead optimized operation plan can meet the balance between supply and demand of electricity, heat and cool load by reasonably arranging the output of each equipment.
When the electricity price is at the valley from 23:00 to 07:00 the next day, the electrical load is low and the thermal load and cool load are mainly powered by CCHP equipment. The electricity generated is enough to satisfy the demand of electrical load, so the electricity is less purchased and the storaging energy of electric storage equipment increases steadily.
Then the electrical load of energy systems reaches the first peak between 7:00 and 11:00. PIES increases the purchased electric power, adjusts the storage equipment to discharge and improves the storaging energy of cool storage equipment.
From 11:00 to 15:00, the first peak of electrical load ends and power generation of photovoltaic and wind turbine gradually rises, so the power supply of energy systems is relatively abundant. System reduces power purchase, while the energy storage of the electric storage equipment increases. Because of the continuous increase of cool load, CCHP equipment can only increase the output to maintain the cool balance, causing the cool storage energy further increases.
During 15:00 to 23:00, PIES ushers in the second peak of electrical load and increases the purchased electric power and adjusts the storage equipment for discharging. At about 19:00, the storage equipment consumes up the stored electric energy, so energy systems increase the purchased electric power again. As the cool load of energy systems also appears the second peak and the output of CCHP equipment changes little, energy systems arrange the cooling storage equipment to discharge energy to satisfy the cooling demand.

Typical Day in Winter
This paper analyzes the daily output of each equipment on a typical day in winter. The power, heating and cooling load, and renewable energy curves of a typical day in winter are shown in Figure 6.

Day-Ahead Optimization Operation Result and Analysis
According to the above example parameters, this paper runs the day-ahead optimal dispatch model of PIES and obtains the dayahead optimal dispatch scheme on a typical day in winter, as shown in Figure 7.
When the electricity price is at the valley from 22:00 to 07: 00 the next day, the power load is low while the thermal load is keeping at a high level. In order to meet the heating demand, the electricity generated by CCHP exceeds the power load, so that energy systems sell electricity, and electric energy storage increases. On the contrary, thermal energy storage is decreasing.  When it comes to the first peak from 7:00 to 11:00, as the heating load decreases, the output of CCHP also decreases. The insufficient power generation is supplemented by electric energy storage and purchased power.
From 11:00 to 15:00, the first peak of electrical load ends. During this period, the thermal load continues to decrease and reaches the lowest value at 13:00, while the output of CCHP remains basically unchanged in order to meet the power supply demand. As a result, thermal and electric energy storage increases.
When it comes to the second peak from 15:00 to 22:00, PIES increases the purchased electric power and adjusts the storage equipment for discharging. CCHP reaches maximum capacity at 21:00 and lasts for 2 hours. At the same time, the thermal load also reaches the peak, and the thermal energy storage has to operate at maximum discharging power.
In particular, cooling energy storage remains roughly the same regardless of how CCHP output changes. It is possible that energy systems release excess cooling energy in order to reduce the maintenance cost of the cooling storage equipment, when the output of CCHP is sufficient to satisfy the cooling load demand.

Comparative Analysis of Day-Ahead, Intra-day and Real-Time Optimal Dispatch
This paper substitutes the result of day-ahead optimization dispatch into the intra-day optimization dispatch model and the intra-day optimization dispatch plan can be obtained through calculation. Then this paper substitutes the intra-day dispatch result into the real-time dispatch model and also obtains the real-time dispatch result. Figure 8; Figure 9; Figure 10 show the curves of purchasing power from electric grid, gas purchasing capacity and demand response under multi-time scale.
As it can be seen from Figure 8, after the adjustment of intraday and real-time dispatch, the interactive power curve between energy systems and power grid fluctuates with the change of electricity price. The electricity purchase decreases during the peak period and the decrease are more obvious at the second peak period. It can be seen from Figure 9 that, in the peak period of electrical load, the gas purchase volume of energy systems under intra-day and real-time dispatching is improved compared with day-ahead dispatching. This is because electricity price is higher than the price of natural gas at the same energy consumption in the peak period, so adding power from CCHP equipment is more cost-effective than purchasing electricity.
As it can be seen from Figure 10, A-type demand response has a good function of peak clipping and valley filling, reducing the load in the peak period and increasing the load in the valley period. B-type demand response also has the function of peak clipping and valley filling, but the effect is weakened compared with A-type. C-type demand response has a small capacity and fluctuates, mainly because real-time dispatching needs to adjust the load response continuously in a short time to smooth down     the fluctuation of load and renewable energy output, so as to ensure the power balance of energy systems. However, it should be noted in particular that, since the output of CCHP reaches its maximum at about 21:00 on a typical day in winter, the demand response has no ability to adjust downward. Table 5 shows the Comparison of operating costs under different dispatch schemes. As it can be seen from Table 5, compared with day-ahead dispatch, the operation cost of multitime dispatch decreases by 2.23% in summer and 2.05% in winter. The above result verifies that the multi-time scale economic optimization dispatch of PIES considering the demand response can reduce the operation cost and improve economy efficiency.

CONCLUSION
In order to further improve the operation economy of PIES, this paper puts forward an economic optimization dispatch strategy under multiple time scales. Then energy systems optimize the dispatch plan by arranging the scheduling strategy of energy purchase, coupling energy equipment and flexible load reasonably. Finally, this paper analyzes the multi-time scale optimal dispatch results and economic benefits of energy systems. The specific contributions are as follows: 1) Adding the demand response conforming to different time scales into the optimization dispatch plan of PIES, which can realize the effect of peak clipping and valley filling, flatten the load curve, and reduce the operation cost of energy systems. 2) Scheduling electricity, heat and cool loads on different time scales by using the three-layer optimal dispatch model which is composed of day-ahead, intra-day and real-time, so that energy systems can balance the supply and demand relationship and ensure the stable operation.
Meanwhile, in view of the problems reflected in the research, this paper considers from the following aspects to further expand: 1) This paper only formulates the response strategy of electric load, without considering the demand response of cooling and heating load. 2) This article will consider coupling other types of energy conversion equipment, including gas Boiler, Power to Gas, and others. Moreover, when CCHP reaches the maximum rated output, the demand response loses the negative adjustment ability. It is considered to add equipment with repetitive functions in CCHP to make up for the deficiency of CCHP capacity. This paper will also consider using gas storage equipment in the future. 3) This paper does not consider the technical characteristics and cost composition of energy transmission and conversion in different PIES, such as industrial PIES, residential PIES and business PIES. It is necessary to verify whether the model is applicable to different types of energy systems.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
PL conceptualized the study, contributed to the study methodology, and wrote the original draft. FZ contributed to the study methodology, data curation and investigation. XM contributed to data analysis and investigation. SY contributed to software and formal analysis. ZLZ contributed to investigation and writing-original draft. PY contributed to supervision and writing-review & editing. ZLZ contributed to investigation and writing-review & editing. CSL contributed to data curation. LLL contributed to writing-review & editing. All authors have read and agreed to the published version of the manuscript.