# Optimal Control Strategy of Platform Load Oriented to Network and Load Cooperation

Authors: Xuesong Shao; Gaoying Cui; Xiao Chen; Xinrong Ji; Yongxian Yi
Journal: Complexity
 (2021)

---

## Abstract

In recent years, with the continuous growth of China’s power peak load and the rapid development of renewable energy, a large number of renewable energy sources are connected to the power grid, increasing the uncertainty of power grid operation and posing new major challenges to the power system regulation capacity. Flexible load has the characteristics of wide distribution, fast response, and high economy, which is an important control resource for the future power system. Based on the flexible load of commercial buildings and residential users, this paper studies the resource characteristics and response characteristics, clarifies the resource characteristics and demand response characteristic indexes of commercial and residential users, and establishes the response characteristics model of commercial buildings and residential users. Considering the influence of weather, holidays, incentive mechanism, and other factors on the response of flexible load, the quantitative analysis method of flexible load resource regulation potential for regional power grid dispatching was studied, and the feasibility of flexible load resources directly participating in the load control system was analyzed. Based on the uncertainty and mathematical characterization method of the active response of flexible loads, the optimal combination control strategy of demand response resources was proposed to eliminate the problems of heavy load and overload of regional power grid equipment by using the active response ability of flexible loads. Finally, the IEEE 14-node system is selected for simulation verification, which provides a theoretical basis for alleviating the power grid operation pressure in the peak load period of the power grid in the urban core area, improving the safety and economic operation level of regional power grid dispatching and the utilization rate of power grid equipment assets.


---

## Body

## 1. Introduction

With the rapid development of China’s economy and the increasing load, the contradiction between power supply and demand is further deepened. Flexible load has the advantages of fast response and high economy and reflects different users’ willingness to use electricity. Flexible load resources are important interactive resources in the framework of smart grid. Although a single load is not controllable at any time, a large number of small- and medium-sized loads will show a certain overall controllability when gathered together.In order to make the flexible loads of commercial and residential users better participate in the system scheduling through the demand response project, it is necessary to master the response characteristics when they participate in the demand response project. Domestic and foreign scholars have done a lot of research on this. The literature [1] studies the load regulation technology of a large number of large apartment air conditioning load clusters, studies the regulation characteristics of the air conditioning load cluster, and proposes a regulation method of state cohort based on this characteristic. The literature [2] proposes a new temperature setting adjustment method based on the load regulation characteristics of air conditioning regulated by the state queue method proposed by Lee et al., in order to improve problems such as state absence and load fluctuation in regulation. The literature [3] establishes the first-order model of a single air conditioner and establishes the power aggregation model of air conditioning load by using the Latin hypervertical method. After the establishment of the aggregation model, the controllable capacity of air conditioning load is evaluated, and the expression between the controllable capacity and the control temperature and outdoor temperature was obtained. The literature [4] compares three different ways to contract design in which demand response participates in power system regulation for incentive payments. The literature [5] mainly analyzes the relationship between the installation capacity of residents’ air conditioning load and residents’ characteristic variables. In the literature [6], the status quo of energy consumption of public buildings at home and abroad and the development process of load calculation methods are systematically reviewed and analyzed, the methods of investigation and field measurement of energy consumption of public buildings and the possible problems are discussed, and some suggestions are given. The literature [7] divides the air conditioning load into three industries, which are commercial users, administrative and public institution users, and residential users, and obtains the industry load through the ratio relationship between the industry electric quantity and the number of hours of the maximum load utilization. The literature [8] proposes a demand response scheme using a hopping scheme for consumption scheduling of appliances which can protect the privacy of electricity consumers. The literature [9] proposes a novel demand-side response model, which can be integrated into system scheduling model straightway and enables the optimal scheduling of energy/power consumption of a heterogeneous population of TCLs. The literature [10] analyzes and studies the load control technology of central air conditioning based on demand response. The literature [11] analyzes the response characteristic that the diversity of the load group was destroyed after the air conditioning load group participated in the regulation, proposes the control strategy to maintain the diversity of the air conditioning load group, analyzes the phenomenon of “secondary peak load” after the air conditioning load group was controlled, and puts forward the time-divided load recovery strategy. The literature [12] studies the response characteristics of refrigerated and nonrefrigerated central air conditioning systems used in public buildings. The literature [13] proposes a simple dynamic model on the basis of the analysis, estimation, and control of a variable-output compressor air conditioner (VOCAC).The literature [14] shows that with the development of smart power grid, the interaction ability of “source-network-load” is improved, and the participation of flexible loads in power grid scheduling makes the management of line overload more flexible and economical. Scholars at home and abroad have studied how to mobilize flexible loads in courts to eliminate heavy load on lines and formed a heavy-load elimination mechanism based on electricity price and incentive demand responses. The literature [15] sorts out and analyzes the lessons of the restructuring of the American energy industry based on the natural gas and electric power accidents in California. In the literature [16], a model predictive control (MPC) of distributed air conditioning loads is proposed in order to balance fluctuations in solar power. The literature [17–19] considers the impact of load demand elasticity on power grid power flow from the perspective of users and aims at maximizing users’ benefits by establishing reasonable electricity price to guide load transfer from peak period to trough period to eliminate line overload. The literature [20–23] invents interruptible loads to alleviate the problem of overloaded lines. In order to minimize the interruption load, interruption compensation price, and power purchase cost of users, the power flow of heavy-duty lines can be reduced by making reasonable load reduction plan. The literature [24] firstly analyzes different types of power loads, builds user load characteristic quantity that reflects the power consumption characteristics of loads, and then obtains user demand response load curve after the implementation of time-of use electricity price. Then, according to the principle of economics, the demand response behavior of users under the environment of TOU electricity price is analyzed, and the corresponding load reduction and transfer model is constructed. The literature [25] encourages users to respond by sharing the benefits brought by reducing grid congestion with users. However, these results have only been verified in a single-time cross-section and have yet to be verified in a continuous-time cross-section. At the same time, simple electricity price or incentive mechanism is difficult to fully mobilize the active response potential of flexible loads, while the coordinated application of electricity price and incentive mechanism can increase the diversity of choices on the user side, which is conducive to the maximum utilization of the demand response potential of flexible loads and effectively reduce the power flow of heavy-load line.The agglomeration effect caused by a large number of flexible load access will affect the existing dispatching operation mode of network in different degrees. At present, there are some research studies on flexible load participating in dispatching operation control of power grid at home and abroad which mainly focusing on the active response characteristics of flexible loads, load characteristics of electric vehicles, demand response scheduling strategies participating in wind power consumption and the influence on rotating reserve capacity, etc. However, there is little research on the model and algorithm of regional power grid optimization reconfiguration with flexible load. Therefore, it is urgent to carry out research on the key technologies of flexible load participating grid dispatching operation to support the integration of demand response dispatching and current dispatching methods.Based on the load response characteristics of commercial buildings and residential users, this paper analyzes the potential of flexible load regulation. On this basis, the flexible load optimal combination control model based on the maximum comprehensive benefit is established, and the medium- and long-term control, day-ahead control, and real-time control methods of flexible load are studied. And the feasibility of flexible load resources directly participating in the load control system was analyzed. Based on the uncertainty and mathematical characterization method of the active response of flexible loads, the optimal combination control strategy of demand response resources was proposed to eliminate the problems of heavy load and overload of regional power grid equipment by using the active response ability of flexible loads. In this paper, the complementary potential of the hybrid energy system and the response priority mechanism of the integrated energy system equipment are analyzed from the aspects of technical and economic factors. The design model of the integrated energy coordinated control is optimized. Finally, the IEEE 14-node system is selected for simulation verification. It provides a theoretical basis for improving the dynamic balance ability of the power grid, alleviating the operating pressure of the power grid and improving the economic operating level of the power grid.


## 2. An Analysis of Load Response Characteristics of Commercial Buildings and Residential Users Oriented to Courts

### 2.1. An Analysis of Load Response Characteristics of Commercial Building

In the central air conditioning system, the operating power of the refrigerator takes the largest proportion of the operating power of the central air conditioning system. Therefore, only the refrigerators in the central air conditioning system are regulated and mathematically modeled.Objective function:(1)maxPpot=maxPair⋅num_pott=Te_start,…,Te_end.In the above formula,Te_start and Te_end represent the start and end moments of emergency demand response events, respectively. During the duration of emergency demand response events, the load reduction of central air conditioning load group may be different at each moment, and the minimum load reduction over this period is defined as the regulation potential Ppot. Pair is the operating power of a single central air conditioning refrigerator; num_pot is a decision variable to be optimized, and its significance is the number of central air conditioning refrigerators actually reduced corresponding to the regulation potential. Therefore, the significance of the objective function is to maximize the regulation potential of the central air conditioning load group during the emergency demand response event.The constraints are as follows:(1)
Load reduction constraints during emergency demand response events(2)Pbaset−∑k=1nPair⋅Sk,t≥Ppott=Te_start,…,Te_end.In equation (2), Pbaset represents the original overall operating power of the central air conditioning group at time t, that is, the load baseline of the central air conditioning load group; ∑k=1nPair⋅Sk,t represents the overall actual operating power after the central air conditioning load group is controlled at time t, and n represents the number of controlled central air conditioning in the central air conditioning load group; Sk,t represents the operating state of the refrigerator of the k-th central air conditioning at time t, Sk,t = 1 means that the refrigerator of the k-th central air conditioning is in operation at time t, and Sk,t = 0 means that the refrigerator of the k-th central air conditioning is in shutdown state at time t.(2)
Room temperature time-varying equationAccording to the energy conservation principle, during the refrigeration period of the central air conditioning system, the indoor energy conservation equation can be obtained:(3)ρa⋅Vk⋅ca⋅dTinont=Qt1⋅dt+Qer⋅dt+Qnew⋅dt−Qt2−Qzon⋅dt,whereρa is air density; Vk is indoor volume; ca is the air’s specific heat capacity; Tinont is the indoor temperature within the period of t in the refrigeration period; and ρa⋅Vk⋅ca⋅dTinontρ represents the heat obtained in the room during the dt period. By solving the first-order differential equation, the time-varying equation of indoor temperature in the refrigeration period can be obtained, shown as follows:(4)Tinont=α3−Qpα2−α3−Qpα2−Tinon0⋅e−α2/α1⋅t.When the central air conditioning system is in shutdown period, its indoor energy conservation equation is(5)ρa⋅Vk⋅ca⋅dTinofft=Qt1⋅dt+Qer⋅dt+Qnew⋅dt−Qt2−Qzoff⋅dt.By solving the above first-order differential equation, the indoor temperature change equation during the shutdown period can be obtained, as follows:(6)Tinofft=α3α2+α4a⋅α1−α2⋅e−a⋅t−α4a⋅α1−α2+α3α2−Tinoff0⋅e−α2/α1⋅t.Specific expressions of some variables in the above equation are as follows:(7)α1=ρa⋅Vk⋅ca+ks⋅Swallin,α2=1.8⋅ktop⋅Stop+kwall⋅Swall,α3=α2⋅Tout+Qer,α4=a⋅mz⋅cw⋅Tw−out−Tw−in,whereVk is the indoor volume; ca is the specific heat capacity of air; ks is the heat storage coefficient of the inner wall; Swallin is the area of the inner wall; ktop and kwall are the coefficients of heat conduction of the roof and wall, respectively; Stop and Swall are the area of the roof and wall; mz is the quality of frozen water; cw is the specific heat capacity of frozen water; Tw−in and Tw−out are, respectively, the inlet and outlet water temperature of frozen water; and Qer is the total cooling load of indoor equipment, lighting, and personnel.(3)
Upper and lower limits on room temperature changes(8)Tink,t=Tinonk,t⋅Sk,t+Tinoffk,t⋅1−Sk,t,Tmint≤Tink,t≤Tmaxt,t=Tp_cool,…,Te_start,…,Te_end,whereTmaxt and Tmint represent the upper and lower limits of the temperature range at each time after the load group of central air conditioning is regulated by load aggregators. Tp_cool represents the start time of the advance notice period. Therefore, if the load aggregators use precooling control strategy, the advance notice period is t=Tp_cool,…,Te_start−1; otherwise, no advance notice is required, and at this time, Tp_cool=Te_start.(4)
Minimum operating time constraint of refrigerator(9)∑t=Tp_coolTp_cool+Gk−11−Sk,t=0,∑t=mm+UTk−1Sk,t≥UTkSk,n−Sk,n−1∀m=Tp_cool+Gk,…,Te_end−UTk+1,∑t=nTe_endSk,t−Sk,t−Sk,t−1≥0∀n=Te_end−UTk+2,…,Te_end,Gk=minTe_end−Tp_cool+1,UTk−Uk0Sk,0.In equation (9), UTk represents the minimum operating time of the k-th central air conditioning refrigerating machine; Uk0 represents the initial operating time of the k-th central air conditioning refrigerating machine; Sk,0 represents the initial operating state of the k-th central air conditioning refrigerating machine; and Gk refers to the required operating time at least of the refrigerating machine in order to keep continuous with the operating state before control during the initial period after the start of the control. The third expression in equation (9) is to ensure that the operation state of the central air conditioning unit remains continuous before the end of the controlled state.(5)
Minimum downtime constraint of the refrigerator(10)∑t=Tp_coolTp_cool+Lk−1Sk,t=0,∑t=mm+DTk−11−Sk,t≥DTkSk,n−1−Sk,n∀m=Tp_cool+Lk,…,Te_end−DTk+1,∑t=nTe_end1−Sk,t−Sk,t−1−Sk,t≥0∀n=Te_end−DTk+2,…,Te_end,Lk=minTe_end−Tp_cool+1,DTk−Qk01−Sk,0.In the above formula,DTk represents the minimum downtime of the k-th central air conditioning refrigerator; Qk0 represents the initial shutdown time of the k-th central air conditioning refrigerating machine; Lk refers to the required shutdown time at least of the refrigerating machine in order to keep continuous with the operating state before control during the initial period after the start of the control. The third expression in equation (10) is to ensure that the operation state of the central air conditioning units remains continuous before the end of the controlled state.(6)
Constraints on the value range of decision variables(11)Sk,t∈0,1t=Tpre−cooling,…,Te_start,…,Te_end,0≤num_pot≤minPbasetPairt=Te_start,…,Te_end.The goal of the above central air conditioning load group regulation potential optimization model is to find the maximum numbernum_pot of central air conditioning refrigerators to be cut corresponding to the regulation potential in the objective function (1) and the operating state Sk,t of each central air conditioning refrigerator at every moment controlled by the load aggregator under the condition that the constraints are satisfied.The refrigerator’s start-stop state matrixSk,t contains the advance notice period of all central air conditioning refrigerators when the precooling regulation strategy is adopted and the state information of the emergency demand response event duration period. According to the state matrix obtained by optimization, the load aggregator can realize the regulation potential by adjusting the operating state of each central air conditioning refrigerator in these two periods. In addition, the room temperature will drop to the lower limit of room temperature at the end of the advance notice period to realize prerefrigeration; during the duration of emergency demand response events, reduce the air conditioning load according to the regulation potential obtained by optimization.After the emergency demand response event, the central air conditioning refrigerators will be artificially and orderly controlled to gradually enter the operation state, so that the room temperature will return to the original room temperature range. When the room temperature enters the original room temperature range, the artificial control of the running state of the refrigerator will be cancelled.In the model established above, only the relevant constraints of the refrigerator are considered; that is, the refrigerator is mainly used to achieve load reduction, which is a reasonable control method based on the consideration of communication and control complexity. If greater reduction is required, the air supply unit can be considered for regulation. When optimizing the air supply unit, the model is as follows.The objective function:(12)maxPpot=maxminPbaset−Pair⋅∑k=1nSairk,t+Pfan⋅∑k=1nSfank,tt=Te_start,…,Te_end.where Pbaset is the baseline load taking into account the power of the forced draught blower; Pair and Pfan are, respectively, the power of refrigerator and air supply system; Sairk,t and Sfank,t are, respectively, the operating states of the refrigerating machine and air supply system in the k-th building at time t. Therefore, the operating power of commercial buildings’ air conditioning load group in time period t is jointly determined by the operating states of refrigerators and air supply systems in each building, as shown in the following formula:(13)Ptotal=Pair⋅∑k=1nSairk,t+Pfan⋅∑k=1nSfank,t.Taking into account the regulation potentialPpot of air conditioning load group of commercial buildings controlled by the air supply system, the expression is as follows:(14)Ppot=minPbaset−Pair⋅∑k=1nSairk,t+Pfan⋅∑k=1nSfank,t.It is the minimum value of load reduction of commercial building air conditioning load group during emergency demand response event.The constraints are as follows:(1)
Room temperature time-varying equationIndoor energy conservation equation:(15)ρa⋅Vk⋅ca⋅dTinofft=Qt1⋅dt+Qer⋅dt+Qnew⋅dt−Qt2.By solving the first-order differential equation, the room temperature time-varying equation under this condition can be obtained:(16)Tinoff2t=α3α2−α3α2−Tinoff0⋅e−α2/α1⋅t.(2)
Coupling constraints of the operating state of refrigerating unit and air supply system(17)Sairk,t−Sfank,t≤0.It can be seen from the above formula that when the refrigerator is in operation stateSairk,t=1, the air supply system must also be in operation state Sfank,t=1; when the refrigerator is off Sairk,t=0, the air supply system can be in either running state Sfank,t=1 or off state Sfank,t=0.(3)
Constraint of running time of the air supply system(18)∑t=Tpre−coolTendSfank,t−Sairk,t≤Tfan_on⋅Timeson−off.The operating time of the air supply system is related to the change of the residual cooling capacity of chilled water. For the convenience of modeling, the running time of the air supply system is simplified to a fixed timeTfan_on; that is, the running time of the air supply system shall not exceed Tfan_on during the shutdown period of the central air conditioning refrigerator.


### 2.2. Analysis of Load Response Characteristics of Residential Users

The objective function is consistent with that of air conditioning load in the commercial buildings. By controlling the air conditioning load of residential users, the maximum load reduction goal is achieved under the condition of satisfying the corresponding constraints, so that the potential of the air conditioning cluster can be fully developed.In addition to one of the above objective functions, due to the variety and quantity of residents’ air conditioning loads, in order to reduce the impact on users, the corresponding number of air conditioning units should be reduced as little as possible under the constraints of the maximum load reduction. The specific objective function is shown as follows:(19)minxrun=min∑k=1nSk,tt=Te_statrt,…,Te_end.The meanings of the above parameter are consistent with those in the air conditioning load model of commercial buildings wherexrun is the newly set variable representing the number of running machines at each time, and the objective function is added to make it affect the air conditioning load as little as possible under the condition of the maximum load reduction.For constraint, it also includes the above constraints, namely, load reduction constraints during emergency demand response events, room temperature time-varying equations, constraints on upper and lower limits of room temperature variation, constraints on minimum operating time of air conditioning, constraints on minimum downtime of air conditioning, constraints on the value range of decision variables, etc. Most of the constraints contain the same contents as the air conditioning cluster load in commercial buildings. The main difference is that the start-stop state matrixSk,t of the air conditioning load only contains the state information of all residents’ air conditioners during the duration of emergency demand response events but does not include the state in the advance notice period when the prerefrigeration control strategy is used. This is mainly because prerefrigeration measures for residential air conditioning load do not have its practicality. Similarly, according to the state matrix obtained by optimization, the load aggregator regulates the operating state of each resident’s air conditioner during this period, so as to realize the regulation potential and reduce the air conditioning load according to the regulation potential obtained by optimization during the duration of emergency demand response events.


## 2.1. An Analysis of Load Response Characteristics of Commercial Building

In the central air conditioning system, the operating power of the refrigerator takes the largest proportion of the operating power of the central air conditioning system. Therefore, only the refrigerators in the central air conditioning system are regulated and mathematically modeled.Objective function:(1)maxPpot=maxPair⋅num_pott=Te_start,…,Te_end.In the above formula,Te_start and Te_end represent the start and end moments of emergency demand response events, respectively. During the duration of emergency demand response events, the load reduction of central air conditioning load group may be different at each moment, and the minimum load reduction over this period is defined as the regulation potential Ppot. Pair is the operating power of a single central air conditioning refrigerator; num_pot is a decision variable to be optimized, and its significance is the number of central air conditioning refrigerators actually reduced corresponding to the regulation potential. Therefore, the significance of the objective function is to maximize the regulation potential of the central air conditioning load group during the emergency demand response event.The constraints are as follows:(1)
Load reduction constraints during emergency demand response events(2)Pbaset−∑k=1nPair⋅Sk,t≥Ppott=Te_start,…,Te_end.In equation (2), Pbaset represents the original overall operating power of the central air conditioning group at time t, that is, the load baseline of the central air conditioning load group; ∑k=1nPair⋅Sk,t represents the overall actual operating power after the central air conditioning load group is controlled at time t, and n represents the number of controlled central air conditioning in the central air conditioning load group; Sk,t represents the operating state of the refrigerator of the k-th central air conditioning at time t, Sk,t = 1 means that the refrigerator of the k-th central air conditioning is in operation at time t, and Sk,t = 0 means that the refrigerator of the k-th central air conditioning is in shutdown state at time t.(2)
Room temperature time-varying equationAccording to the energy conservation principle, during the refrigeration period of the central air conditioning system, the indoor energy conservation equation can be obtained:(3)ρa⋅Vk⋅ca⋅dTinont=Qt1⋅dt+Qer⋅dt+Qnew⋅dt−Qt2−Qzon⋅dt,whereρa is air density; Vk is indoor volume; ca is the air’s specific heat capacity; Tinont is the indoor temperature within the period of t in the refrigeration period; and ρa⋅Vk⋅ca⋅dTinontρ represents the heat obtained in the room during the dt period. By solving the first-order differential equation, the time-varying equation of indoor temperature in the refrigeration period can be obtained, shown as follows:(4)Tinont=α3−Qpα2−α3−Qpα2−Tinon0⋅e−α2/α1⋅t.When the central air conditioning system is in shutdown period, its indoor energy conservation equation is(5)ρa⋅Vk⋅ca⋅dTinofft=Qt1⋅dt+Qer⋅dt+Qnew⋅dt−Qt2−Qzoff⋅dt.By solving the above first-order differential equation, the indoor temperature change equation during the shutdown period can be obtained, as follows:(6)Tinofft=α3α2+α4a⋅α1−α2⋅e−a⋅t−α4a⋅α1−α2+α3α2−Tinoff0⋅e−α2/α1⋅t.Specific expressions of some variables in the above equation are as follows:(7)α1=ρa⋅Vk⋅ca+ks⋅Swallin,α2=1.8⋅ktop⋅Stop+kwall⋅Swall,α3=α2⋅Tout+Qer,α4=a⋅mz⋅cw⋅Tw−out−Tw−in,whereVk is the indoor volume; ca is the specific heat capacity of air; ks is the heat storage coefficient of the inner wall; Swallin is the area of the inner wall; ktop and kwall are the coefficients of heat conduction of the roof and wall, respectively; Stop and Swall are the area of the roof and wall; mz is the quality of frozen water; cw is the specific heat capacity of frozen water; Tw−in and Tw−out are, respectively, the inlet and outlet water temperature of frozen water; and Qer is the total cooling load of indoor equipment, lighting, and personnel.(3)
Upper and lower limits on room temperature changes(8)Tink,t=Tinonk,t⋅Sk,t+Tinoffk,t⋅1−Sk,t,Tmint≤Tink,t≤Tmaxt,t=Tp_cool,…,Te_start,…,Te_end,whereTmaxt and Tmint represent the upper and lower limits of the temperature range at each time after the load group of central air conditioning is regulated by load aggregators. Tp_cool represents the start time of the advance notice period. Therefore, if the load aggregators use precooling control strategy, the advance notice period is t=Tp_cool,…,Te_start−1; otherwise, no advance notice is required, and at this time, Tp_cool=Te_start.(4)
Minimum operating time constraint of refrigerator(9)∑t=Tp_coolTp_cool+Gk−11−Sk,t=0,∑t=mm+UTk−1Sk,t≥UTkSk,n−Sk,n−1∀m=Tp_cool+Gk,…,Te_end−UTk+1,∑t=nTe_endSk,t−Sk,t−Sk,t−1≥0∀n=Te_end−UTk+2,…,Te_end,Gk=minTe_end−Tp_cool+1,UTk−Uk0Sk,0.In equation (9), UTk represents the minimum operating time of the k-th central air conditioning refrigerating machine; Uk0 represents the initial operating time of the k-th central air conditioning refrigerating machine; Sk,0 represents the initial operating state of the k-th central air conditioning refrigerating machine; and Gk refers to the required operating time at least of the refrigerating machine in order to keep continuous with the operating state before control during the initial period after the start of the control. The third expression in equation (9) is to ensure that the operation state of the central air conditioning unit remains continuous before the end of the controlled state.(5)
Minimum downtime constraint of the refrigerator(10)∑t=Tp_coolTp_cool+Lk−1Sk,t=0,∑t=mm+DTk−11−Sk,t≥DTkSk,n−1−Sk,n∀m=Tp_cool+Lk,…,Te_end−DTk+1,∑t=nTe_end1−Sk,t−Sk,t−1−Sk,t≥0∀n=Te_end−DTk+2,…,Te_end,Lk=minTe_end−Tp_cool+1,DTk−Qk01−Sk,0.In the above formula,DTk represents the minimum downtime of the k-th central air conditioning refrigerator; Qk0 represents the initial shutdown time of the k-th central air conditioning refrigerating machine; Lk refers to the required shutdown time at least of the refrigerating machine in order to keep continuous with the operating state before control during the initial period after the start of the control. The third expression in equation (10) is to ensure that the operation state of the central air conditioning units remains continuous before the end of the controlled state.(6)
Constraints on the value range of decision variables(11)Sk,t∈0,1t=Tpre−cooling,…,Te_start,…,Te_end,0≤num_pot≤minPbasetPairt=Te_start,…,Te_end.The goal of the above central air conditioning load group regulation potential optimization model is to find the maximum numbernum_pot of central air conditioning refrigerators to be cut corresponding to the regulation potential in the objective function (1) and the operating state Sk,t of each central air conditioning refrigerator at every moment controlled by the load aggregator under the condition that the constraints are satisfied.The refrigerator’s start-stop state matrixSk,t contains the advance notice period of all central air conditioning refrigerators when the precooling regulation strategy is adopted and the state information of the emergency demand response event duration period. According to the state matrix obtained by optimization, the load aggregator can realize the regulation potential by adjusting the operating state of each central air conditioning refrigerator in these two periods. In addition, the room temperature will drop to the lower limit of room temperature at the end of the advance notice period to realize prerefrigeration; during the duration of emergency demand response events, reduce the air conditioning load according to the regulation potential obtained by optimization.After the emergency demand response event, the central air conditioning refrigerators will be artificially and orderly controlled to gradually enter the operation state, so that the room temperature will return to the original room temperature range. When the room temperature enters the original room temperature range, the artificial control of the running state of the refrigerator will be cancelled.In the model established above, only the relevant constraints of the refrigerator are considered; that is, the refrigerator is mainly used to achieve load reduction, which is a reasonable control method based on the consideration of communication and control complexity. If greater reduction is required, the air supply unit can be considered for regulation. When optimizing the air supply unit, the model is as follows.The objective function:(12)maxPpot=maxminPbaset−Pair⋅∑k=1nSairk,t+Pfan⋅∑k=1nSfank,tt=Te_start,…,Te_end.where Pbaset is the baseline load taking into account the power of the forced draught blower; Pair and Pfan are, respectively, the power of refrigerator and air supply system; Sairk,t and Sfank,t are, respectively, the operating states of the refrigerating machine and air supply system in the k-th building at time t. Therefore, the operating power of commercial buildings’ air conditioning load group in time period t is jointly determined by the operating states of refrigerators and air supply systems in each building, as shown in the following formula:(13)Ptotal=Pair⋅∑k=1nSairk,t+Pfan⋅∑k=1nSfank,t.Taking into account the regulation potentialPpot of air conditioning load group of commercial buildings controlled by the air supply system, the expression is as follows:(14)Ppot=minPbaset−Pair⋅∑k=1nSairk,t+Pfan⋅∑k=1nSfank,t.It is the minimum value of load reduction of commercial building air conditioning load group during emergency demand response event.The constraints are as follows:(1)
Room temperature time-varying equationIndoor energy conservation equation:(15)ρa⋅Vk⋅ca⋅dTinofft=Qt1⋅dt+Qer⋅dt+Qnew⋅dt−Qt2.By solving the first-order differential equation, the room temperature time-varying equation under this condition can be obtained:(16)Tinoff2t=α3α2−α3α2−Tinoff0⋅e−α2/α1⋅t.(2)
Coupling constraints of the operating state of refrigerating unit and air supply system(17)Sairk,t−Sfank,t≤0.It can be seen from the above formula that when the refrigerator is in operation stateSairk,t=1, the air supply system must also be in operation state Sfank,t=1; when the refrigerator is off Sairk,t=0, the air supply system can be in either running state Sfank,t=1 or off state Sfank,t=0.(3)
Constraint of running time of the air supply system(18)∑t=Tpre−coolTendSfank,t−Sairk,t≤Tfan_on⋅Timeson−off.The operating time of the air supply system is related to the change of the residual cooling capacity of chilled water. For the convenience of modeling, the running time of the air supply system is simplified to a fixed timeTfan_on; that is, the running time of the air supply system shall not exceed Tfan_on during the shutdown period of the central air conditioning refrigerator.


## 2.2. Analysis of Load Response Characteristics of Residential Users

The objective function is consistent with that of air conditioning load in the commercial buildings. By controlling the air conditioning load of residential users, the maximum load reduction goal is achieved under the condition of satisfying the corresponding constraints, so that the potential of the air conditioning cluster can be fully developed.In addition to one of the above objective functions, due to the variety and quantity of residents’ air conditioning loads, in order to reduce the impact on users, the corresponding number of air conditioning units should be reduced as little as possible under the constraints of the maximum load reduction. The specific objective function is shown as follows:(19)minxrun=min∑k=1nSk,tt=Te_statrt,…,Te_end.The meanings of the above parameter are consistent with those in the air conditioning load model of commercial buildings wherexrun is the newly set variable representing the number of running machines at each time, and the objective function is added to make it affect the air conditioning load as little as possible under the condition of the maximum load reduction.For constraint, it also includes the above constraints, namely, load reduction constraints during emergency demand response events, room temperature time-varying equations, constraints on upper and lower limits of room temperature variation, constraints on minimum operating time of air conditioning, constraints on minimum downtime of air conditioning, constraints on the value range of decision variables, etc. Most of the constraints contain the same contents as the air conditioning cluster load in commercial buildings. The main difference is that the start-stop state matrixSk,t of the air conditioning load only contains the state information of all residents’ air conditioners during the duration of emergency demand response events but does not include the state in the advance notice period when the prerefrigeration control strategy is used. This is mainly because prerefrigeration measures for residential air conditioning load do not have its practicality. Similarly, according to the state matrix obtained by optimization, the load aggregator regulates the operating state of each resident’s air conditioner during this period, so as to realize the regulation potential and reduce the air conditioning load according to the regulation potential obtained by optimization during the duration of emergency demand response events.


## 3. Establish an Optimal Control Model for the Flexible Load Participating in the Power Grid Dispatching

### 3.1. Analysis Method of Flexible Load Regulation Potential

In this paper, Xinjiekou Zhuanjinying Line# 1(hereinafter referred to as Jinying # 1) is selected as a typical representative of the power characteristics of the courts to analyze the load characteristics.Table1 shows the linear regression model summary of the load under the meteorological factor Jinying # 1 feeder in summer (June, July, and August), focusing on the value of the determination coefficient R squared. This coefficient indicates the percentage of the variance explained by the independent variable in the total variance. The greater the value, the greater the role of this factor. In the following table, the determination coefficient R reaches 71.7%, indicating that the regression relationship can explain the 71.7% variation of the dependent variable and has a good fitting degree.Table 1
Summary table of the linear regression model of Jinying# 1 line load in summer.
ModelRAdjustedRError in standard estimates10.7170.70318.423An analysis of variance was conducted on the model. The result of the analysis of variance showed that theF value was 54.974, and the P value was less than 0.05. Therefore, the model was statistically significant or valuable for predicting variables.The Table2 gives the estimated value and test results of the constant term and regression coefficient in the regression equation. It can be seen that the constant term = −31.001, and the four regression coefficients are 139.900, 30.094, 12.826, and 7.430, respectively. The regression equation can be written as follows:(20)load value=−31.001+139.900×temperature factor+30.094×humidity factor+12.826×precipitation factor+7.430×wind speed factor.Table 2
Table of coefficient of the linear regression model of Jinying# 1 line load in summer.
ModelNonstandard coefficientNormalized coefficientTSigBStandard errorBeta1(Constant)−31.00112.083—−2.5660.0122Temperature factor139.9009.8091.39714.2620.0003Humidity factor30.0943.7190.5648.0920.0004Precipitation factor12.8261.5000.7058.5500.0005Wind speed factor7.4302.2950.2073.2380.002Table3 shows a summary of the linear regression model of winter meteorological factor Jinying # 1 load. The R value in the model abstract is 0.859, indicating that the regression relationship can explain 85.9% variation of the dependent variable. In Table 4, Sig value is 0, less than 0.05, indicating that most independent variables in the model are correlated with dependent variables, indicating that the model has statistical significance.Table 3
Summary of the linear regression model of Jinying# 1 line load in winter.
ModelRAdjustedRError in standard estimates10.8590.81711.033Table 4
Variance diagram of the line load linear regression model of Jinying# 1 in winter.
ModelQuadratic sumDfMean squareFSig (b)Regression6484.93541621.23413.3190.000Residual5112.45642121.725Total11597.39146Table5 gives the estimated value and test results of the constant term and regression coefficient in the regression equation. The regression equation can be written as follows:(21)load value=103.295−32.514×temperature factor−2.414×humidity factor−5.539×precipitation factor−0.884×wind speed factor.Table 5
Table of coefficient of the linear regression model of Jinying# 1 line load in winter.
ModelNonstandard coefficientNormalized coefficientTSigBStandard errorBeta1(Constant)103.2952.247—45.9780.0002Temperature factor−32.5145.042−0.698−6.4480.0003Humidity factor−2.4141.520−0.170−1.5880.1204Precipitation factor−5.5393.026−0.191−1.8300.0045Wind speed factor−0.8842.205−0.042−0.4010.691According to the above construction method of the multifactor load characteristic regression model, each feeder can be independently analyzed, and the corresponding regression model of each feeder can be obtained, so that the feeder load value can be accurately predicted according to meteorological factors. According to different models in spring and summer, the load baseline mean values of feeders in spring and summer can be obtained. According to the difference between the two, the theoretical maximum load reduction potential can be obtained. According to this value, different load reduction potential targets can be formulated. After a large number of simulation calculations, the load reduction potential of each feeder under different prerefrigeration duration, different control duration, and different room temperature can be obtained. Finally, the flexible load regulation potential table can be obtained, as shown in Table6.Table 6
Flexible load regulation potential table.
Outdoor temperature 31–33°CPrecooling duration (min)Duration of event (min)Upper limit of room temperature (°C)Load reduction potential (MW)Compensation fees (yuan)03026552574027995405428115.57207229115.5810814026442956827776036828115.510348829115.5116424103026115.559459.427115.569369.328115.579279.229115.589189.1402660.544721.6279985377.628115.5113836.829115.5128066.4


### 3.2. Establish an Optimal Combination Control Model of Flexible Load Based on the Maximum Comprehensive Benefit

When the grid dispatching needs to regulate an interactive resource to complete a certain interactive task, the grid needs to pay the corresponding cost to the user, that is, the interactive cost of flexible load dispatching on grid side, no matter whether the price type load or the incentive load responds. At this time, with the change of the relationship between electricity price and power, the income of power grid will also change.When the interactive resources on the flexible load side respond according to the price signal or incentive measures sent by the power grid, the electricity expense of the flexible load side will be different while the relationship between the electricity price and power changes. The variation of electricity expense is the interactive benefit of the flexible load side load response. Accordingly, the inverse value of the interaction benefit of load response represents the interactive cost of load response on the flexible load side. At the same time, the mode of electricity consumption will also change accordingly, which will affect the satisfaction of the mode of electricity consumption.Cinter, the comprehensive interaction cost of flexible load demand response, can be expressed as(22)Cinter=∑t∈nT∑i∈npCi,tp+∑j∈nICj,tI+∑m∈nGCm,tg,Cm,tg=amPm,tg2+bmPm,tg+λm,tcm,where Cinter is the comprehensive interaction cost of flexible load demand response; Ci,tp is the scheduling cost of power-price type flexible load i in time period t; Cj,tI is the dispatching cost of incentive flexible load j in time period t; np is the quantity of flexible load of electricity price type; nI is the quantity of excitation type flexible load; nG is the amount of generators in the system; Cm,tg is the generating cost of the generator m in the system during the period of time t; Pm,tg is the active power output of the generator m during the period of time t; am,bm,andcm are the power generation cost model coefficient of generator m; and λm,t is a variable of 0-1, which is 1 when the generator is put into operation in time period t, and 0 otherwise.In this paper, the optimal dispatching is carried out aiming at maximizing the comprehensive benefits of “source-network-load.” Since the numerical value of the comprehensive satisfaction degree of electricity consumptionR and the comprehensive interaction cost Cinter do not belong to the same order of magnitude, it is necessary to normalize the Cinter in order to ensure the optimization effect.The objective function is as follows:(23)MaxF=Max1np∑i∈npRi−Cinter∑t∈nT∑k∈nfPk,t0×p0,where Ri is the comprehensive electricity consumption satisfaction of power-price type load i; nf is the total number of flexible loads; p0 is the original unified power price of the system; and Pk,t0 is the original load of the flexible load k during the period of time t. There are mainly three constraints, which are unit constraint, system constraint, and load constraint, and the details are as follows:(1)
Unit constraintsThe upper and lower limits of unit output:(24)Pg,mmin≤Pm,tg≤Pg,mmax,m∈nG,wherePg,mmin and Pg,mmax indicate the lower and upper limits of the output of generator mseparately.Unit climbing constraint:(25)λm,tPm,tg−λm,t−1Pm,t−1g≤Ug,m,λm,t−1Pm,t−1g−λm,tPm,tg≤Dg,m,whereUg,m and Dg,m, respectively, represent the upward and downward climbing constraint values of generator m.Minimum start-stop time constraint of unit:(26)λm,t−λm,t−1⋅Tt−1−Ton_m≥0,λm,t−λm,t−1⋅Tt−1−Toff_m≤0,whereTt−1 is the duration of period t−1 and Ton_m and Toff_m are the minimum start-stop time of generator m, respectively.(2)
System constraintsActive power balance constraint:(27)∑k∈nfPk,t=∑m∈nGPm,tg,wherePk,t is the load of the flexible load k during the period of time t after participating in the demand response.Network constraints:(28)−bxy×θx,t−θy,t≤Lxy,∀t∈nT,∀x,y∈1,2,…,N,∐x≠y,wherebxy is the electrical susceptance of line x−y; θx,t and θy,t separately represent the phase angles of nodes x and y; Lxy is the maximum active power transmission capacity for line x−y; and N is the amount of nodes in the system.(3)
Load constraintThere are mainly two load constraints which are load maximum adjustment constraints and excitation response invocation time constraints, and the details are as follows:(29)Pkmin≤Pk,t≤Pkmax,k∈nf,0≤numj≤nI_jmax,j∈nI,where Pkmin and Pkmax are the minimum and maximum of the incentive flexible load j, respectively; numj is the number of times that the incentive flexible load j is actually called; and nI_jmax is the maximum number of times that the incentive flexible load j is called.The dispatcher can get the optimal nodal power price, load reduction period and load reduction amount of each incentive flexible load, and the output arrangement of each power generating unit in each period through optimization, so as to issue scheduling instructions.


## 3.1. Analysis Method of Flexible Load Regulation Potential

In this paper, Xinjiekou Zhuanjinying Line# 1(hereinafter referred to as Jinying # 1) is selected as a typical representative of the power characteristics of the courts to analyze the load characteristics.Table1 shows the linear regression model summary of the load under the meteorological factor Jinying # 1 feeder in summer (June, July, and August), focusing on the value of the determination coefficient R squared. This coefficient indicates the percentage of the variance explained by the independent variable in the total variance. The greater the value, the greater the role of this factor. In the following table, the determination coefficient R reaches 71.7%, indicating that the regression relationship can explain the 71.7% variation of the dependent variable and has a good fitting degree.Table 1
Summary table of the linear regression model of Jinying# 1 line load in summer.
ModelRAdjustedRError in standard estimates10.7170.70318.423An analysis of variance was conducted on the model. The result of the analysis of variance showed that theF value was 54.974, and the P value was less than 0.05. Therefore, the model was statistically significant or valuable for predicting variables.The Table2 gives the estimated value and test results of the constant term and regression coefficient in the regression equation. It can be seen that the constant term = −31.001, and the four regression coefficients are 139.900, 30.094, 12.826, and 7.430, respectively. The regression equation can be written as follows:(20)load value=−31.001+139.900×temperature factor+30.094×humidity factor+12.826×precipitation factor+7.430×wind speed factor.Table 2
Table of coefficient of the linear regression model of Jinying# 1 line load in summer.
ModelNonstandard coefficientNormalized coefficientTSigBStandard errorBeta1(Constant)−31.00112.083—−2.5660.0122Temperature factor139.9009.8091.39714.2620.0003Humidity factor30.0943.7190.5648.0920.0004Precipitation factor12.8261.5000.7058.5500.0005Wind speed factor7.4302.2950.2073.2380.002Table3 shows a summary of the linear regression model of winter meteorological factor Jinying # 1 load. The R value in the model abstract is 0.859, indicating that the regression relationship can explain 85.9% variation of the dependent variable. In Table 4, Sig value is 0, less than 0.05, indicating that most independent variables in the model are correlated with dependent variables, indicating that the model has statistical significance.Table 3
Summary of the linear regression model of Jinying# 1 line load in winter.
ModelRAdjustedRError in standard estimates10.8590.81711.033Table 4
Variance diagram of the line load linear regression model of Jinying# 1 in winter.
ModelQuadratic sumDfMean squareFSig (b)Regression6484.93541621.23413.3190.000Residual5112.45642121.725Total11597.39146Table5 gives the estimated value and test results of the constant term and regression coefficient in the regression equation. The regression equation can be written as follows:(21)load value=103.295−32.514×temperature factor−2.414×humidity factor−5.539×precipitation factor−0.884×wind speed factor.Table 5
Table of coefficient of the linear regression model of Jinying# 1 line load in winter.
ModelNonstandard coefficientNormalized coefficientTSigBStandard errorBeta1(Constant)103.2952.247—45.9780.0002Temperature factor−32.5145.042−0.698−6.4480.0003Humidity factor−2.4141.520−0.170−1.5880.1204Precipitation factor−5.5393.026−0.191−1.8300.0045Wind speed factor−0.8842.205−0.042−0.4010.691According to the above construction method of the multifactor load characteristic regression model, each feeder can be independently analyzed, and the corresponding regression model of each feeder can be obtained, so that the feeder load value can be accurately predicted according to meteorological factors. According to different models in spring and summer, the load baseline mean values of feeders in spring and summer can be obtained. According to the difference between the two, the theoretical maximum load reduction potential can be obtained. According to this value, different load reduction potential targets can be formulated. After a large number of simulation calculations, the load reduction potential of each feeder under different prerefrigeration duration, different control duration, and different room temperature can be obtained. Finally, the flexible load regulation potential table can be obtained, as shown in Table6.Table 6
Flexible load regulation potential table.
Outdoor temperature 31–33°CPrecooling duration (min)Duration of event (min)Upper limit of room temperature (°C)Load reduction potential (MW)Compensation fees (yuan)03026552574027995405428115.57207229115.5810814026442956827776036828115.510348829115.5116424103026115.559459.427115.569369.328115.579279.229115.589189.1402660.544721.6279985377.628115.5113836.829115.5128066.4


## 3.2. Establish an Optimal Combination Control Model of Flexible Load Based on the Maximum Comprehensive Benefit

When the grid dispatching needs to regulate an interactive resource to complete a certain interactive task, the grid needs to pay the corresponding cost to the user, that is, the interactive cost of flexible load dispatching on grid side, no matter whether the price type load or the incentive load responds. At this time, with the change of the relationship between electricity price and power, the income of power grid will also change.When the interactive resources on the flexible load side respond according to the price signal or incentive measures sent by the power grid, the electricity expense of the flexible load side will be different while the relationship between the electricity price and power changes. The variation of electricity expense is the interactive benefit of the flexible load side load response. Accordingly, the inverse value of the interaction benefit of load response represents the interactive cost of load response on the flexible load side. At the same time, the mode of electricity consumption will also change accordingly, which will affect the satisfaction of the mode of electricity consumption.Cinter, the comprehensive interaction cost of flexible load demand response, can be expressed as(22)Cinter=∑t∈nT∑i∈npCi,tp+∑j∈nICj,tI+∑m∈nGCm,tg,Cm,tg=amPm,tg2+bmPm,tg+λm,tcm,where Cinter is the comprehensive interaction cost of flexible load demand response; Ci,tp is the scheduling cost of power-price type flexible load i in time period t; Cj,tI is the dispatching cost of incentive flexible load j in time period t; np is the quantity of flexible load of electricity price type; nI is the quantity of excitation type flexible load; nG is the amount of generators in the system; Cm,tg is the generating cost of the generator m in the system during the period of time t; Pm,tg is the active power output of the generator m during the period of time t; am,bm,andcm are the power generation cost model coefficient of generator m; and λm,t is a variable of 0-1, which is 1 when the generator is put into operation in time period t, and 0 otherwise.In this paper, the optimal dispatching is carried out aiming at maximizing the comprehensive benefits of “source-network-load.” Since the numerical value of the comprehensive satisfaction degree of electricity consumptionR and the comprehensive interaction cost Cinter do not belong to the same order of magnitude, it is necessary to normalize the Cinter in order to ensure the optimization effect.The objective function is as follows:(23)MaxF=Max1np∑i∈npRi−Cinter∑t∈nT∑k∈nfPk,t0×p0,where Ri is the comprehensive electricity consumption satisfaction of power-price type load i; nf is the total number of flexible loads; p0 is the original unified power price of the system; and Pk,t0 is the original load of the flexible load k during the period of time t. There are mainly three constraints, which are unit constraint, system constraint, and load constraint, and the details are as follows:(1)
Unit constraintsThe upper and lower limits of unit output:(24)Pg,mmin≤Pm,tg≤Pg,mmax,m∈nG,wherePg,mmin and Pg,mmax indicate the lower and upper limits of the output of generator mseparately.Unit climbing constraint:(25)λm,tPm,tg−λm,t−1Pm,t−1g≤Ug,m,λm,t−1Pm,t−1g−λm,tPm,tg≤Dg,m,whereUg,m and Dg,m, respectively, represent the upward and downward climbing constraint values of generator m.Minimum start-stop time constraint of unit:(26)λm,t−λm,t−1⋅Tt−1−Ton_m≥0,λm,t−λm,t−1⋅Tt−1−Toff_m≤0,whereTt−1 is the duration of period t−1 and Ton_m and Toff_m are the minimum start-stop time of generator m, respectively.(2)
System constraintsActive power balance constraint:(27)∑k∈nfPk,t=∑m∈nGPm,tg,wherePk,t is the load of the flexible load k during the period of time t after participating in the demand response.Network constraints:(28)−bxy×θx,t−θy,t≤Lxy,∀t∈nT,∀x,y∈1,2,…,N,∐x≠y,wherebxy is the electrical susceptance of line x−y; θx,t and θy,t separately represent the phase angles of nodes x and y; Lxy is the maximum active power transmission capacity for line x−y; and N is the amount of nodes in the system.(3)
Load constraintThere are mainly two load constraints which are load maximum adjustment constraints and excitation response invocation time constraints, and the details are as follows:(29)Pkmin≤Pk,t≤Pkmax,k∈nf,0≤numj≤nI_jmax,j∈nI,where Pkmin and Pkmax are the minimum and maximum of the incentive flexible load j, respectively; numj is the number of times that the incentive flexible load j is actually called; and nI_jmax is the maximum number of times that the incentive flexible load j is called.The dispatcher can get the optimal nodal power price, load reduction period and load reduction amount of each incentive flexible load, and the output arrangement of each power generating unit in each period through optimization, so as to issue scheduling instructions.


## 4. A Comprehensive Optimization Control Strategy of Flexible Load for Courts

Through the above analysis, the comprehensive optimization control strategy of flexible load under different time scales can be obtained to eliminate the line overload problem.


### 4.1. Medium- and Long-Term Regulation of Flexible Loads

For medium- and long-term load dispatching planning, the main consideration is the demand response model of load under the electricity price incentive of time-of-use electricity price and interruptible load. Its model is mainly manifested as response modeling of cumulative electric quantity at various periods of peak and valley level:(30)qAp=∑i=1nqAip,qAm=∑i=1nqAim,qAv=∑i=1nqAiv,qAip=qAip0+εppiqAip0pTOUip−pTOUip0pTOUip0+εpviqAip0pTOUiv−pTOUiv0pTOUiv0−θiΔLILit2i−t1i,qAim=qAim0,qAiv=qAiv0+εvviqAiv0pTOUiv−pTOUiv0pTOUiv0+εvpiqAiv0pTOUip−pTOUip0pTOUip0,where n represents the number of loads under a certain node; θi is the probability of the i-th load signing an interruptible contract; ΔLILi is the amount of interruptible load of load i; and t1i∼t2i represents the start-stop time of the interruptible load for load i.


### 4.2. Day-Ahead Regulation of Flexible Load

For day-ahead load dispatching planning, the main consideration is the demand response model under day-ahead real-time power price and interruptible load:(31)PCt=∑i=1nPCit,PCit=PCitpSPit−ΔPILi.


### 4.3. Real-Time Regulation of Flexible Load

For real-time load dispatching planning, the demand response model of load under the incentive of real-time power price is considered. Its main features are shown as follows:(1)
Price elasticity of demand will be different while power price changes.(2)
The load at the next moment is the result of superposition of load response components based on the load at the previous moment.(32)PDt+1=∑i=1nPDit+1,PDit+1=PDit+εipSPit+1PDitpSPit+1−pSPitpSPit,εipSPit+1=∂PipSPit+1∂pSPit+1.According to the above analysis, the expectation and variance of load under different time scale responses can be obtained, as shown in Table7.Table 7
Load expectation and variance at different time scales.
Time scaleElectricity price incentivesExpectationVarianceMedium and long termTime-of-use electricity price + interruptible loadqApqAmqAvσAp2=σAMC2σAm2=0.264lnΔqAm+1σAv2=0.264lnΔqAv+1Day-aheadDay-ahead real-time electricity price + interruptible loadPCtσC2=σCMCt2t1≤t≤t20.264lnΔLCt+1elseReal-timeReal-time electricityPDtσD2=0.264lnΔLDt+1Note:σAMCt2, σBMCt2, and σCMCt2 represent the load variance values obtained by Monte Carlo simulation.


## 4.1. Medium- and Long-Term Regulation of Flexible Loads

For medium- and long-term load dispatching planning, the main consideration is the demand response model of load under the electricity price incentive of time-of-use electricity price and interruptible load. Its model is mainly manifested as response modeling of cumulative electric quantity at various periods of peak and valley level:(30)qAp=∑i=1nqAip,qAm=∑i=1nqAim,qAv=∑i=1nqAiv,qAip=qAip0+εppiqAip0pTOUip−pTOUip0pTOUip0+εpviqAip0pTOUiv−pTOUiv0pTOUiv0−θiΔLILit2i−t1i,qAim=qAim0,qAiv=qAiv0+εvviqAiv0pTOUiv−pTOUiv0pTOUiv0+εvpiqAiv0pTOUip−pTOUip0pTOUip0,where n represents the number of loads under a certain node; θi is the probability of the i-th load signing an interruptible contract; ΔLILi is the amount of interruptible load of load i; and t1i∼t2i represents the start-stop time of the interruptible load for load i.


## 4.2. Day-Ahead Regulation of Flexible Load

For day-ahead load dispatching planning, the main consideration is the demand response model under day-ahead real-time power price and interruptible load:(31)PCt=∑i=1nPCit,PCit=PCitpSPit−ΔPILi.


## 4.3. Real-Time Regulation of Flexible Load

For real-time load dispatching planning, the demand response model of load under the incentive of real-time power price is considered. Its main features are shown as follows:(1)
Price elasticity of demand will be different while power price changes.(2)
The load at the next moment is the result of superposition of load response components based on the load at the previous moment.(32)PDt+1=∑i=1nPDit+1,PDit+1=PDit+εipSPit+1PDitpSPit+1−pSPitpSPit,εipSPit+1=∂PipSPit+1∂pSPit+1.According to the above analysis, the expectation and variance of load under different time scale responses can be obtained, as shown in Table7.Table 7
Load expectation and variance at different time scales.
Time scaleElectricity price incentivesExpectationVarianceMedium and long termTime-of-use electricity price + interruptible loadqApqAmqAvσAp2=σAMC2σAm2=0.264lnΔqAm+1σAv2=0.264lnΔqAv+1Day-aheadDay-ahead real-time electricity price + interruptible loadPCtσC2=σCMCt2t1≤t≤t20.264lnΔLCt+1elseReal-timeReal-time electricityPDtσD2=0.264lnΔLDt+1Note:σAMCt2, σBMCt2, and σCMCt2 represent the load variance values obtained by Monte Carlo simulation.


## 5. Example Analysis

The IEEE 33-node system in the literature [26] was selected for example analysis. The IEEE 33-node system is shown in Figure 1. There are 5 generators in the system, which are located on nodes 8, 12, 17, 19, and 22. Tables 8–10 show the system’s nodes, active power output, circuit, and load parameters. The unit of cost coefficient a,b,c are $/MW2, $/MW, and $, respectively, and the units of impedance and ground capacitance are Ω.Figure 1
IEEE 33-node system.Table 8
Active power output data of the IEEE 33-node system.
Node numberUpper limit of active power output/(MW)Lower limit of active power output/(MW)Cost coefficientabc830000130012300002300173500033501935000635022500008500Table 9
Line parameters of the IEEE 33-node system.
Branch circuitBranch impedanceGround capacitance0-10.0922 +j0.0470.04121-20.01938 +j0.059170.05282-30.04699 +j0.197970.04383-40.05811 +j0.176320.03744-50.05403 +j0.223040.04925-60.05695 +j0.173880.03406-70.06701 +j0.171030.03467-80.00000 +j0.176150.00007-90.00000 +j0.110010.00008-91.0440 +j0.74000.00009-100.1966 +j006500.000010-110.3744 +j0.12380.000011-121.4680 +j1.5500.000012-130.5416 +j0.71290.000013-140.5910 +j0.56200.000014-150.7463 +j0.54500.000015-161.2890 +j1.72100.00001-180.09498 +j0.19890.000018-190.12291 +j0.255810.000019-200.06615 +j0.130270.000020-210.12711 +j0.270380.00002-220.08205 +j0.192070.000022-230.22092 +j0.199880.000023-240.17093 +j0.348020.00005-250.2030 +j0.10340.000025-260.2842 +j0.14470.000026-271.0590 +j0.93370.000027-280.8042 +j0.70060.000028-290.5075 +j0.25850.000029-300.9744 + 0.96300.000030-310.3105 +j0.36190.000031-320.3410 +j0.53620.0000Table 10
Load demand response attributes of each user.
Load numberMaximum /average loadLoad start timeLoad end timeMinimum running timeMinimum /average loadMinimum downturn timeIncentive scheduling cost (first stage)Incentive scheduling cost (second stage)Incentive scheduling cost (third stage)Stand-by costTotal limit calls number of stimulated DRMaximum reserve capacity of incentive DRPDR/IDRLoad type11.1512440.7545101538200PDR&IDRTranslational21.1512430.753510203830PDR&IDRTransferable31.1512420.752510253830PDR&IDRReducible41.1512420.752515303840PDR&IDRReducible51.1512410.75050505020180PDR&IDRReducible61.1512420.752491438200PDR&IDRTranslational71.1512430.75349192830PDR&IDRTransferable81.1512420.75249242830PDR&IDRReducible91.1512420.752414292840PDR&IDRReducible101.1512410.75050505020180PDR&IDRReducible111.1512400.7501000100010001000160PDR121.1512400.750100010001000100080PDR131.1512400.750100010001000100080PDR141.1512400.750100010001000100080PDR151.1512400.750100010001000100010PDRIn Table10, all users can participate in price-based demand response. Users 1 to 10 participate in incentive demand response, which is divided into translational, reducible, and transferable classes. Take μ1=μ2=0.5, take K=1,α=0.15, and the original price p0=42.7 in the electricity price response model. When the breaker breaking capacity constraint is not taken into account, the optimal power flow is adopted to calculate and obtain the active power output data of each generator in 24 periods, as shown in Table 11. The active power flow of each branch is shown in Table 12.Table 11
Active power output of each generator without considering line overload constraint.
Generator123456789101112131415161718192021222324825.00.00.00.00.00.00.00.075.0150.0225.0289.4232.9157.982.97.90.075.0150.0225.0225.0150.075.00.0120.00.00.00.00.00.00.024.591.5166.5241.5221.3296.3225.0150.075.00.074.7149.7100.2166.291.216.20.01762.50.087.50.02.890.3177.8265.3350.0340.8350.0350.0350.0350.0262.5175.0247.1320.1350.0350.0350.0324.2236.7149.219225.1255.4276.0188.5149.8175.0262.5350.0322.7328.1350.0350.0350.0350.0332.8297.0350.0262.5350.0350.0350.0262.5175.087.522125.0250.0375.0476.8500.0407.2500.0500.0500.0500.0500.0500.0500.0500.0500.0500.0500.0500.0500.0500.0500.0500.0382.4396.7Table 12
Active power flow in 24 periods of each line without considering line capacity constraint.
Circuit1234567891011121314151617181920212223240-19.218.623.18.92.3−13.6−4.57.4−24.0−35.4−77.3−89.631.238.549.457.371.792.4115.2131.691.766.362.945.21-2−44.410.037.930.626.128.645.255.543.611.6−21.9−43.1−35.81.330.456.657.654.825.99.743.473.439.626.72-3−33.3−7.3−0.35.89.68.32.72.6−19.2−56.3−94.7−114.8−105.5−67.9−28.03.24.3−1.7−40.8−64.5−49.4−10.85.68.03-437.870.895.866.041.154.191.0118.4117.488.962.648.656.285.5105.3117.6118.0124.3102.881.795.6114.889.154.54-5−34.0−10.3−3.713.526.820.910.36.3−16.9−57.8−99.6−120.5−109.6−69.4−26.33.95.41.7−39.8−68.4−56.3−15.113.319.35-6−27.5−12.8−13.3−3.03.50.9−11.9−15.5−39.5−76.1−114.0−132.6−123.1−87.6−46.2−15.2−14.0−20.8−61.2−86.1−78.1−38.8−6.31.16-7−78.8−92.6−114.6−62.5−20.0−41.1−94.8−130.6−153.2−162.4−175.1−180.5−178.0−170.5−149.0−132.0−131.1−142.2−160.0−165.0−168.6−148.5−89.5−43.37-8−30.986.4198.7323.6449.9430.5426.5417.5414.3407.6386.7379.9386.7400.7401.6404.7405.3416.5405.7396.7406.4413.3421.0425.87-9−21.9−75.2−146.3−195.2−252.8−250.7−272.0−283.0−300.7−324.1−349.3−357.6−354.1−340.7−298.0−275.1−275.5−286.2−317.2−331.2−326.6−301.1−262.0−247.98-928.8−9.5−39.6−68.7−97.7−83.8−92.2−90.5−92.5−72.3−53.4−43.1−49.3−71.0−80.7−79.0−80.4−92.7−86.3−69.2−87.2−96.8−81.8−75.99-10−75.2−28.55.580.0149.3129.8100.676.450.414.3−24.0−46.3−31.82.641.075.973.972.135.05.622.155.6101.9129.710-11−33.0−26.2−27.8−9.55.1−1.9−17.4−28.8−42.3−60.7−80.4−90.6−84.3−68.6−45.3−27.3−28.2−31.1−51.5−65.5−58.3−40.4−14.8−1.411-1210.821.922.12.0−18.1−6.1−5.54.1−25.0−64.6−104.9−129.2−114.8−77.2−29.715.113.31.1−46.9−76.0−64.7−21.27.72.512-13−65.2−63.8−74.1−37.4−7.8−24.7−53.7−77.0−86.5−96.2−108.2−110.6−110.9−103.2−88.2−81.4−80.1−82.0−93.7−97.9−96.7−85.4−55.6−28.613-14−30.2−34.5−45.5−38.5−33.5−37.7−51.9−65.3−71.6−79.0−85.7−91.5−88.4−81.5−71.4−61.4−63.7−66.3−75.1−81.5−78.9−71.0−49.7−36.114-15−71.5−78.1−101.3−78.8−62.5−75.5−108.1−139.2−153.7−170.3−188.8−196.4−194.3−179.9−153.7−132.4−137.0−142.6−163.1−174.7−170.1−152.2−104.5−73.115-160.0−13.5−38.7−64.3−91.9−83.2−90.0−84.2−102.9−113.9−127.5−131.3−133.6−120.7−99.6−79.1−76.1−92.3−118.4−117.8−128.3−102.4−79.7−75.516-17−30.4−43.9−76.2−88.6−106.0−105.2−123.9−139.6−153.5−170.0−189.3−196.0−195.5−180.3−150.6−124.7−130.5−139.7−162.2−175.2−170.8−152.4−113.3−96.51-1844.638.036.53.5−25.4−9.55.823.223.426.129.730.029.428.625.628.428.423.923.126.021.722.810.5−6.118-19−7.9−6.8−7.8−2.12.2−0.9−4.0−7.2−8.5−9.8−13.4−10.9−13.5−12.9−8.8−7.7−7.8−8.2−10.1−9.6−9.9−8.4−4.5−1.419-20−58.8−59.1−74.4−47.0−27.0−42.1−67.8−92.7−99.1−110.1−123.7−126.6−126.3−118.2−99.8−90.3−93.1−92.7−102.5−112.4−105.1−98.0−67.2−42.220-21−32.6−57.2−63.1−71.0−93.4−112.2−134.5−142.8−167.3−172.1−181.6−194.5−52.1−66.3−68.0−76.2−80.2−86.3−91.2−99.4−83.8−78.4−63.2−55.62-2222.411.27.62.1−19.6−11.611.227.432.138.944.751.3−130.2−113.5−101.1−92.6−77.3−62.1−57.9−48.5−39.2−27.9−25.3−22.122-23−52.914.026.937.255.468.248.132.422.39.2−17.3−35.2−178.2−152.4−139.5−122.9−107.1−97.4−81.8−63.5−71.6−77.4−55.2−38.523-2444.351.265.782.996.4116.3137.2132.6117.478.372.162.9−31.5−29.3−27.8−22.9−19.6−17.4−15.3−14.0−11.9−8.2−6.9−3.05-25−32.1−28.9−22.4−13.6−9.4−3.85.610.5−17.2−44.3−59.3−88.9−152.6−134.7−98.0−92.3−81.0−114.9−108.6−97.2−84.2−62.9−55.8−42.625-26−99.7−103.2−119.3−88.4−54.7−45.9−66.9−78.2−92.4−127.3−130.6−164.2−128.9−110.2−77.3−48.6−18.311.827.1−36.0−44.7−57.0−23.6−16.726-274.222.839.699.7118.4132.6100.783.479.354.221.913.449.255.663.959.178.294.371.399.882.538.443.237.127-28−5.6−23.4−11.97.925.739.610.63.2−17.7−48.1−66.2−98.9−112.8−97.6−88.3−47.0−13.923.713.1−9.239.564.762.644.128-29−46.5−22.1−3.45.914.6−13.2−19.4−31.1−39.3−58.3−72.0−81.922.141.055.271.989.2118.446.927.410.954.984.689.729-3023.213.97.8−3.4−16.7−20.6−41.2−65.9−72.8−99.3−121.7−128.4−33.6−27.4−19.6−9.48.724.933.542.946.437.136.826.530-31−44.2−34.1−22.8−13.2−3.28.3−4.9−27.9−34.1−48.3−63.9−87.9−182.0−178.6−150.4−144.9−127.3−92.9−75.3−62.7−77.9−65.7−56.2−42.631-3232.830.524.916.78.739.648.961.388.692.5111.9126.531.728.426.122.919.316.115.712.58.87.45.0−5.9The maximum breaking capacity of circuit breakers on each line is 300 MW. It can be seen from Table12 that the power flow of lines 7-8 and 7–9 exceeds the limit in some time periods. To eliminate the network congestion of lines 7-8 and 7–9, four different scenarios are constructed in the literature [26]. Scenario settings are shown in Table 13.Table 13
Four different scenarios for eliminating network congestion.
Scene sequenceScene settingsScenario 1Eliminating network congestion only by generation side dispatchingScenario 2Adjust generation output and implement electricity price dispatch on load side to eliminate network congestionScenario 3Adjust generation output and implement incentive dispatch on load side to eliminate network congestionScenario 4Adjust generation output and implement electricity price and incentive dispatch on load side to eliminate network congestion (integrated dispatch of generation and consumption)In the four scenarios, MATLAB was used for simulation calculation. The computer CPU was Intel(R) Core(TM) I5-4590, the main frequency was 3.3 GHz, and the memory was 8 GB. The optimized simulation time in each scenario is 10.57 s, 24.07 s, 28.1 s, and 57.03 s, respectively.The optimized active power flow curves of lines 7-8 and 7–9 are shown in Figures2 and 3.Figure 2
Active power flow curve of line 7-8 after optimized scheduling in different scenarios.Figure 3
Active power flow curve of line 7–9 after optimized scheduling in different scenarios.In Scenario 1, the output of the generator is adjusted without considering the demand response. The optimized output of each generator node is shown in Table14.Table 14
The variation of active power output of each generator when only the output regulation of the generator is considered.
Generator12345678910111213141516171819202122232480.00.00.00.00.00.05.527.227.227.227.210.663.567.167.167.10.00.00.00.00.00.00.00.0120.00.00.00.00.00.00.0−24.5−27.2−27.2−27.2−10.6−63.5−67.1−67.1−67.10.00.00.00.00.00.00.00.017112.6230.4199.4199.4187.284.784.784.7−7.60.00.00.00.00.00.06.321.70.00.00.00.0−61.7−61.7−61.719−112.6−230.4−199.4−188.5−149.8−87.5−87.5−87.57.60.00.00.00.00.00.0−6.3−21.70.00.00.00.061.761.761.7220.00.00.0−10.9−37.42.8−2.70.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.0As can be seen from Figures2 and 3, in Scenario 1, the power flow of the blocked line can only be reduced to 340 MW by adjusting the output distribution of the power generation side to the minimum, and the problem of power flow beyond the limit cannot be eliminated. At the same time, it can be seen from Table 14 that when the overload of line active power flow is eliminated, the variation of active power output of generator nodes 17 and 19 is too large, which is difficult to be realized in the actual power grid.However, in Scenarios 2, 3, and 4, due to the implementation of integrated dispatching for power generation and consumption side at the same time, the off-limit power flow is under ideal control, and under this dispatching mode, the blocked line still has room for power flow to fall.


### 5.1. Interactive Cost Analysis of Flexible Loads

The interaction costs of flexible loads in each scenario are shown in Table15.Table 15
Comparison of load interaction costs under different scenarios.
ScenarioTotal cost of interactionTotal cost of generationElectricity price response costIncentive response costScenario 173343073343000Scenario 2749350733359159910Scenario 375966075153908121Scenario 474319073192994621799The comparison of unit power generation costs in each scenario is shown in Table16:Table 16
Comparison of unit power generation cost after response scheduling in different scenarios.
ScenarioTotal cost of generationGross generationUnit generation costScenario 173343027440.4326.73Scenario 273335927419.2126.75Scenario 375153927506.6427.32Scenario 473192927593.0826.53Since the network congestion cannot be eliminated only through generation side dispatching in Scenario 1, the interaction cost in Scenario 1 is not analyzed.As can be seen from Table15, in Scenario 4, when electricity price and incentive comprehensive response scheduling are adopted, the total interactive cost is the minimum. This is mainly because the collaborative scheduling on the generating and consuming side can effectively reduce the output of high-cost units on the basis of ensuring the power consumption and increase the output of the low-cost units, so as to effectively reduce the total cost of power generation and effectively reduces the cost of electricity price response and incentive response of load under the comprehensive action of electricity price and incentive at the same time.And in Scenario 3, although the interaction response cost on the load side (electricity price response cost + incentive response cost) is the lowest, it calls more high-cost units (the total output of units 8 and 12 is 4756 MW in Scenario 3, only 3961 MW in Scenario 4), making its unit power generation cost (from Table16) higher, pushing up the total cost of the interaction. In addition, it can be found from the data that, under the same dispatching adjustment effect, the total interactive cost of electricity price response is obviously smaller than that of incentive response.


### 5.2. Comparison of Comprehensive Satisfaction with Electricity Consumption

Under different scenarios, the comprehensive satisfaction of electricity consumption is shown in Table17:Table 17
Comparison of comprehensive satisfaction with electricity consumption after response scheduling in different scenarios.
ScenarioSatisfaction of power consumption modeInteractive benefit satisfactionComprehensive satisfactionScenario 1111Scenario 20.8151.0430.929Scenario 3111Scenario 40.9870.9960.992As can be seen from the chart, since there is no electricity response in Scenarios 1 and 3, electricity comprehensive satisfaction is 1, and the comprehensive satisfaction with electricity consumption in Scenario 4 is greater than that in Scenario 2. This is mainly due to the addition of incentive response mechanism in Scenario 4, and the incentive scheduling plan is formulated according to the load interruptible plan independently reported by users. In this case, the autonomous power consumption willingness of flexible load can be guaranteed to the maximum extent. Leading to a high degree of satisfaction with its power consumption mode. In Scenario 2, although the user obtains a certain income from electricity fee through the response of electricity price, the mode of electricity consumption changes greatly, which leads to a low degree of comprehensive satisfaction with electricity consumption.It can be found from the comparison of comprehensive economic benefits and satisfaction that when network congestion occurs, in Scenario 4, the optimal scheduling mode of adjusting power generation output and simultaneously implementing electricity price and incentive scheduling on the load side to eliminate network congestion (integrated comprehensive scheduling of power generation and electricity consumption) is relatively the most appropriate.


## 5.1. Interactive Cost Analysis of Flexible Loads

The interaction costs of flexible loads in each scenario are shown in Table15.Table 15
Comparison of load interaction costs under different scenarios.
ScenarioTotal cost of interactionTotal cost of generationElectricity price response costIncentive response costScenario 173343073343000Scenario 2749350733359159910Scenario 375966075153908121Scenario 474319073192994621799The comparison of unit power generation costs in each scenario is shown in Table16:Table 16
Comparison of unit power generation cost after response scheduling in different scenarios.
ScenarioTotal cost of generationGross generationUnit generation costScenario 173343027440.4326.73Scenario 273335927419.2126.75Scenario 375153927506.6427.32Scenario 473192927593.0826.53Since the network congestion cannot be eliminated only through generation side dispatching in Scenario 1, the interaction cost in Scenario 1 is not analyzed.As can be seen from Table15, in Scenario 4, when electricity price and incentive comprehensive response scheduling are adopted, the total interactive cost is the minimum. This is mainly because the collaborative scheduling on the generating and consuming side can effectively reduce the output of high-cost units on the basis of ensuring the power consumption and increase the output of the low-cost units, so as to effectively reduce the total cost of power generation and effectively reduces the cost of electricity price response and incentive response of load under the comprehensive action of electricity price and incentive at the same time.And in Scenario 3, although the interaction response cost on the load side (electricity price response cost + incentive response cost) is the lowest, it calls more high-cost units (the total output of units 8 and 12 is 4756 MW in Scenario 3, only 3961 MW in Scenario 4), making its unit power generation cost (from Table16) higher, pushing up the total cost of the interaction. In addition, it can be found from the data that, under the same dispatching adjustment effect, the total interactive cost of electricity price response is obviously smaller than that of incentive response.


## 5.2. Comparison of Comprehensive Satisfaction with Electricity Consumption

Under different scenarios, the comprehensive satisfaction of electricity consumption is shown in Table17:Table 17
Comparison of comprehensive satisfaction with electricity consumption after response scheduling in different scenarios.
ScenarioSatisfaction of power consumption modeInteractive benefit satisfactionComprehensive satisfactionScenario 1111Scenario 20.8151.0430.929Scenario 3111Scenario 40.9870.9960.992As can be seen from the chart, since there is no electricity response in Scenarios 1 and 3, electricity comprehensive satisfaction is 1, and the comprehensive satisfaction with electricity consumption in Scenario 4 is greater than that in Scenario 2. This is mainly due to the addition of incentive response mechanism in Scenario 4, and the incentive scheduling plan is formulated according to the load interruptible plan independently reported by users. In this case, the autonomous power consumption willingness of flexible load can be guaranteed to the maximum extent. Leading to a high degree of satisfaction with its power consumption mode. In Scenario 2, although the user obtains a certain income from electricity fee through the response of electricity price, the mode of electricity consumption changes greatly, which leads to a low degree of comprehensive satisfaction with electricity consumption.It can be found from the comparison of comprehensive economic benefits and satisfaction that when network congestion occurs, in Scenario 4, the optimal scheduling mode of adjusting power generation output and simultaneously implementing electricity price and incentive scheduling on the load side to eliminate network congestion (integrated comprehensive scheduling of power generation and electricity consumption) is relatively the most appropriate.


## 6. Conclusions

In this paper, the mathematical model of air conditioning load in commercial buildings and residential users is established, and the response characteristics of commercial buildings and residential users are analyzed. On this basis, Xinjiekou Zhuanjinying Line# is selected as the typical representative of the power consumption characteristics of the courts to analyze the load characteristics, and the flexible load regulation potential table is obtained. A flexible load optimal combination control model based on the maximum comprehensive benefit was established, and through this model, the dispatcher can obtain the optimal node electricity price, load reduction period and load reduction amount of each incentive flexible load, and the output arrangement of each generator set in each period, so as to issue scheduling instructions. Finally, the comprehensive optimization control strategy of the flexible load oriented to the courts is proposed, and the medium- and long term-regulation, day-ahead regulation, and real-time regulation methods of the flexible load are studied, and the IEEE 14-node system is selected for example analysis.In this paper, the complementary potential of the hybrid energy system and the response priority mechanism of the integrated energy system equipment are analyzed from the aspects of technical and economic factors. The design model of the integrated energy coordinated control is optimized, and the solution analysis is carried out by using the multitime scale cooperative control strategy. Through the analysis of the example, the optimal scheduling results of multitime scale of source storage load are obtained. This study is of great significance for the construction of the hybrid energy system, the promotion of new energy consumption, the increase of energy storage application in power grid, and the realization of more sustainable energy management.



---

*Source: 9976473-2021-06-08.xml*