A distribution network dispatch optimization method based on aggregatable resources

Abstract

The present invention proposes a distribution network dispatch optimization method with aggregatable resources as the main body, which belongs to the field of smart grid technology and comprises the following steps: S1: constructing an aggregated power quantification model of aggregatable resources based on the operating characteristics of the aggregatable resources, dynamically quantifying the response power of the aggregatable resources, and realizing the coordinated response among the aggregatable resources; S2: constructing a distribution network Distflow branch flow model by considering three aspects: node power balance, node voltage correlation and line power; S3: considering the network loss of the optimized system operation, the operating cost of each unit in the optimized system operation, the maximum absorption of local distributed renewable energy and the reduction of abandoned power, and establishing an optimization objective function; S4: converging constraint conditions, solving the model constructed in S1-S3, and after calculating the results, sending them to the resource main body under the distribution network for execution. The method can realize multi-objective optimization of low network loss, low operating cost and maximum absorption of renewable energy in the local power grid.

Description

Technical Field

[0001] The present invention belongs to the technical field of smart grids, and in particular relates to a distribution network dispatching optimization method based on aggregatable resources. Background Art

[0002] Under the background of energy structure transformation, renewable energy represented by wind power and photovoltaic power is connected to modern power system in large quantities. Unlike traditional power generation units, renewable energy is constrained by natural factors, and its output is intermittent, volatile and uncertain. At the same time, with the increase of renewable energy penetration, the proportion of synchronous generators in operation decreases. The traditional mode of relying solely on the spare capacity of the power generation side for regulation is difficult to cope with the development needs of the future power generation system with a high proportion of renewable energy.

[0003] In addition, it is worth noting that with the large-scale construction and application of distributed power generation, energy storage, electric vehicle charging facilities and microgrids, source load nodes have appeared in the distribution network. Moreover, with the development of advanced measurement systems and Internet of Things technology, various loads can also adjust their own electricity demand according to incentives or electricity prices, especially electric vehicles, heat pumps, electric water heaters, etc. The aggregable resources on the demand side have become an effective dispatching and control object of the distribution network. Therefore, under the premise of comprehensively considering the renewable energy on the power generation side and the aggregable resources on the load side, a new distribution network operation and dispatching method is proposed to achieve multi-objective optimization of low network loss, low operating cost and maximum absorption of renewable energy in the local power grid, which has important research significance. Summary of the invention

[0004] In order to solve the above problems existing in the prior art, a distribution network scheduling optimization method based on aggregatable resources is provided.

[0005] The technical solution adopted by the present invention to solve its technical problem is:

[0006] The present invention proposes a distribution network dispatch optimization method based on aggregatable resources, which is characterized by

[0007] In, comprising the following steps:

[0008] S1: Based on the operating characteristics of aggregatable resources, an aggregation power quantification model of aggregatable resources is constructed to dynamically quantify the response power of aggregatable resources and realize coordinated response among aggregatable resources;

[0009] S2: Considering the node power balance, node voltage correlation and line power, the Distflow branch power flow model of the distribution network is constructed;

[0010] S3: Consider optimizing the network loss of the system operation, optimizing the operating cost of each unit in the system operation, maximizing the absorption of local distributed renewable energy and reducing power abandonment, and establishing the optimization objective function;

[0011] S4: Aggregate constraints, solve the model constructed in S1-S3, and after calculating the results, send them to the resource entities under the distribution network for execution.

[0012] Furthermore, in S1, the aggregated resources include heat pumps, electric water heaters, electric vehicles,

[0013] Transfer loads, translatable loads, and interruptible loads.

[0014] Furthermore, in S1, the modeling process of the heat pump model includes the following steps:

[0015] Considering the cooling service demand, the maximum operating power during the evaluation period corresponds to the power consumed by the current temperature state to reach the lower limit of the controlled temperature, as follows:

[0016] ;

[0017] In the formula, Indicates that the lower temperature limit was reached during the evaluation period Required cooling power, represents the equivalent heat capacity of the building walls and room air, Indicates the lower limit of the comfortable temperature allowed in the room. Indicates the room temperature. represents the time step, represents the radiant heat power of the room, Indicates the thermal power of the heat source in the room, Indicates the outdoor temperature, Indicates the equivalent thermal resistance of the building wall and air;

[0018] Considering the limitation of heat pump operating power, the maximum heat power allowed in the heat pump hot zone is calculated as follows:

[0019] (2);

[0020] In the formula, It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the upper limit of the cooling power of the heat pump equipment. Indicates that the lower temperature limit was reached during the evaluation period Required cooling power;

[0021] The minimum operating power of the heat pump during the evaluation period corresponds to the power consumed by the current temperature state to reach the controlled temperature upper limit, and the formula is as follows:

[0022] ;

[0023] In the formula, Indicates that the upper temperature limit was reached during the evaluation period Required cooling power, represents the equivalent heat capacity of the building walls and room air, Indicates the upper limit of the comfortable temperature allowed in the room. Indicates the room temperature. represents the time step, represents the radiant heat power of the room, Indicates the thermal power of the heat source in the room, Indicates the outdoor temperature, Indicates the equivalent thermal resistance of the building wall and air;

[0024] Considering that the operating power of the heat pump equipment cannot be negative, the minimum thermal power allowed in the heat pump thermal zone is calculated as follows:

[0025] ;

[0026] In the formula, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. Indicates that the upper temperature limit was reached during the evaluation period Required cooling power;

[0027] Calculate the aggregate flexible power range of the fixed-frequency heat pump cluster during the evaluation period using the following formula:

[0028] ;

[0029] In the formula, Indicates the lower limit of the aggregate power of the fixed-frequency heat pump cluster, Indicates the aggregate power upper limit of the fixed-frequency heat pump cluster. represents the scale of the fixed-frequency heat pump cluster, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the operating efficiency coefficient of the fixed-frequency heat pump in the room;

[0030] Calculate the aggregate flexible power range of the variable frequency heat pump cluster during the evaluation period using the following formula:

[0031] ;

[0032] In the formula, Indicates the lower limit of the aggregate power of the variable frequency heat pump cluster; Indicates the aggregate power upper limit of the variable frequency heat pump cluster; represents the size of the variable frequency heat pump cluster, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the operating efficiency coefficient of the variable frequency heat pump in the room.

[0033] Furthermore, in S1, the modeling process of the electric water heater model includes the following steps:

[0034] During the evaluation period, at the current hot water temperature, the hot water temperature cannot exceed the upper limit of the water temperature under the maximum operating power of the electric water heater. Taking into account the power limit of the equipment, the maximum thermal power allowed by the electric water heater is calculated. The formula is as follows:

[0035] ;

[0036] In the formula, Indicates that the water temperature of the electric water heater reaches The average operating power required, Indicates the upper limit of water temperature. Indicates the water temperature of the electric water heater. Indicates the heat capacity of the electric water heater. represents the time step, represents the thermal disturbance power caused by hot water use, Indicates the indoor temperature. represents the thermal resistance of the water heater, Indicates the operating efficiency of the electric water heater. It represents the maximum operating power allowed for electric water heater unit i during the evaluation period. Indicates the upper limit of the operating power of the electric water heater;

[0037] During the evaluation period, at the current hot water temperature, the minimum operating power of the electric water heater must meet the user's hot water demand at the set time, while taking into account the equipment operating power limit, and calculate the minimum thermal power allowed by the electric water heater. The formula is as follows:

[0038] ;

[0039] (13);

[0040] In the formula, Indicates that the water temperature of the electric water heater reaches the set time In the corresponding period The average operating power within represents the time step, Indicates that the water temperature of the electric water heater reaches the set time In the corresponding period The average operating power within Indicates the desired water temperature. Indicates the water temperature of the electric water heater. Indicates the heat capacity of the electric water heater. Indicates the set time for the water heater to work. represents the thermal disturbance power caused by hot water use, Indicates the indoor temperature. represents the thermal resistance of the water heater, Indicates the operating efficiency of the electric water heater. Indicates the upper limit of the operating power of the electric water heater. represents the minimum operating power of electric water heater unit i during the evaluation period;

[0041] The aggregate flexible power range of the electric water heater cluster during the evaluation period is calculated as follows:

[0042] (14);

[0043] (15);

[0044] In the formula, Indicates the lower limit of the aggregate power of the electric water heater cluster. Indicates the aggregate power upper limit of the electric water heater cluster. represents the size of the electric water heater cluster, represents the minimum operating power of electric water heater unit i during the evaluation period, It represents the maximum allowed operating power of electric water heater unit i during the evaluation period.

[0045] Furthermore, in S1, the modeling process of the electric vehicle model includes the following steps:

[0046] During the evaluation period, under the current SOC of the electric vehicle, the maximum charging power of the electric vehicle meets the SOC constraint. Considering the power limit of the charging equipment, the maximum charging power of the electric vehicle is calculated. The formula is as follows:

[0047] ;

[0048] In the formula, Indicates that the SOC of the electric vehicle reaches The average charging power required, represents the upper limit of charging power of electric vehicle cell i, It represents the maximum charging power allowed for electric vehicle single i during the evaluation period, Indicates the maximum value of the battery state of charge. Indicates the battery charge state of an electric vehicle. Indicates the battery capacity. represents the charging efficiency of electric vehicles, represents the time step;

[0049] During the evaluation period, under the current SOC of the electric vehicle, the maximum discharge power of the electric vehicle during the evaluation period should meet the SOC requirement of the user at the time of departure. The maximum discharge power allowed by the electric vehicle is calculated as follows:

[0050] ;

[0051] (19);

[0052] In the formula, Indicates that the SOC of the electric vehicle can meet the requirements in the subsequent period The average discharge power allowed during the evaluation period is: represents the upper limit of charging power of electric vehicle cell i, represents the charging efficiency of electric vehicles, Indicates the time when the electric vehicle leaves. represents the time step, represents the expected state of charge of the electric vehicle when it drives away, Indicates the battery charge state of an electric vehicle. Indicates the battery capacity. represents the discharge efficiency of electric vehicles, It represents the maximum discharge power allowed by electric vehicle cell i during the evaluation period, It represents the upper limit of discharge power of electric vehicle cell i;

[0053] The aggregated flexible power range of the electric vehicle cluster during the evaluation period is calculated as follows:

[0054] (20);

[0055] (twenty one);

[0056] In the formula, represents the lower limit of the aggregated power of the electric vehicle cluster, represents the aggregate power cap of the electric vehicle cluster, represents the size of the electric vehicle cluster, It represents the maximum discharge power allowed by electric vehicle cell i during the evaluation period, The maximum charging power allowed for electric vehicle cell i during the assessment period.

[0057] Furthermore, in S1, the modeling process of the transferable load includes the following steps:

[0058] Taking into account the flexibility of system operation regulation, the transferable load can flexibly deploy the load power within the allowed time period, and the operating power regulation equation formula is as follows:

[0059] (twenty two);

[0060] ;

[0061] In the formula, is the load after transfer, is the load before transfer, For transferable time periods, is a binary variable, if , ,otherwise, , is a set of scheduling time periods, and These are the lower and upper limits of the adjustable range of TFL.

[0062] Furthermore, in S1, the modeling process of the translatable load includes the following steps:

[0063] Taking into account industrial production processes, shiftable loads provide operational flexibility by shifting the load curve.

[0064] The operating power regulation equation is as follows:

[0065] Calculate the load after translation using the following formula:

[0066] (twenty four);

[0067] In the formula, represents the load after translation, represents the decision vector for translation, represents the translation scheme matrix;

[0068] Decision vector for translation The calculation formula is as follows:

[0069] (25);

[0070] (26);

[0071] In the formula, is a binary variable, Indicates that the corresponding translation scheme is selected, otherwise it is not selected;

[0072] Translation Matrix The calculation formula is as follows:

[0073] (27);

[0074] In the formula, The number of translation options for the translatable load, Indicates different scheduling time periods.

[0075] Furthermore, in S1, the modeling process of the interruptible load includes the following steps:

[0076] By giving economic compensation, the interruptible load can be interrupted or reduced according to the user's wishes.

[0077] The operating power regulation equation is as follows:

[0078] ;

[0079] In the formula, represents the load before the interruption, represents the load after the interruption, and are the lower and upper limits of the range in which interrupts are allowed.

[0080] Furthermore, in S2, the modeling process of the distribution network Distflow branch flow model includes:

[0081] Next steps:

[0082] Establish a radial distribution network model;

[0083] Considering node power balance, node voltage correlation and line power, the distribution network Distflow is constructed

[0084] Branch power flow model, where the node power balance equation is as follows:

[0085] ;

[0086] In the formula, the subscript They all represent line nodes. Representation Node The injected active power, Representation Node The injected reactive power, Representation Node , The active power of the line between Representation Node , The active power of the line between Representation Node , The reactive power of the line between Representation Node , The reactive power of the line between Representation Node , The square of the line current amplitude between Representation Node , The line resistance between Representation Node , The line reactance between

[0087] Calculate the node voltage correlation, the formula is as follows:

[0088] ;

[0089] In the formula, Representation Node The square of the voltage amplitude, Representation Node The square of the voltage amplitude, Represents the line resistance, Representation Node , The active power of the line between Representation Node , The line reactance between Representation Node , The reactive power of the line between Representation Node The square of the line current amplitude between them;

[0090] Calculate the line power using the following formula:

[0091] (32);

[0092] In the formula, Representation Node , The active power of the line between Representation Node , The reactive power of the line between Representation Node The square of the voltage amplitude, Representation Node , The square of the line current amplitude between them.

[0093] Furthermore, in S3, the calculation process of optimizing the objective function includes the following steps:

[0094] The network loss cost of the distribution network is calculated using the following formula:

[0095] ;

[0096] In the formula, represents the network loss cost target item, represents the network loss cost coefficient of the distribution network, represents the calculation time period, represents the set of distribution network lines, Representation Node , The square of the line current amplitude between Representation Node , The line resistance between Indicates the optimization scheduling time interval;

[0097] The operating cost of the distribution network is calculated as follows:

[0098] (34);

[0099] In the formula, represents the operating cost target item, represents the tie line power cost coefficient, is the tie line power, represents the static VAR generator (SVG) cost factor, Indicates SVG power, represents the cost coefficient of controllable distributed generator (DG), represents the DG power, represents the energy storage (ES) cost coefficient, Indicates ES charging power, represents the ES discharge power, represents the load response cost coefficient, Indicates load response regulation power, Indicates the optimization scheduling time interval, Indicates the calculation time period;

[0100] Load response power regulation The calculation formula is as follows:

[0101] (35);

[0102] In the formula, represents the aggregated constant-frequency heat pump (HPss) operating power, Indicates the load before the constant frequency heat pump (HPss) load is controlled. represents the aggregated variable speed heat pump (HPiv) operating power, represents the load before the aggregated variable speed heat pump (HPiv) type load is controlled, represents the aggregated electric water heater (EWH) operating power, represents the load before the aggregated electric water heater (EWH) type load is controlled, represents the aggregated electric vehicle (EV) operating power, represents the load before the electric vehicle (EV) load is controlled, represents the aggregated transferable load (TFL) operating power, represents the load before the transferable load (TFL) type load is controlled, represents the aggregate Translatable Load (TLL) operating power, represents the load before the Translatable Load (TLL) type load is controlled, represents the aggregated interruptible load (IL) operating power, It indicates the load before the interruptible load (IL) type load is controlled;

[0103] The formula for calculating the RE curtailment cost of the distribution network is as follows:

[0104] (36);

[0105] In the formula, represents the target item of power curtailment cost, represents the calculation time period, represents the cost coefficient of RE power abandonment in the distribution network, represents the abandoned power of photovoltaic (PV), represents the abandoned power of the wind turbine (WT), Indicates the optimization scheduling time interval;

[0106] Calculate the optimal scheduling target, the formula is as follows:

[0107] (37);

[0108] In the formula, represents the optimization scheduling target item, represents the network loss cost target item, represents the operating cost target item, Represents the curtailment cost objective item.

[0109] Furthermore, in S4, the constraints consist of the following parts: load response constraints, tie line power constraints, static VAR generator constraints, distributed power supply regulation constraints, energy storage response constraints and distribution network system operation constraints, a total of six parts. Among them, the load response constraints include fixed frequency heat pump constraints, variable frequency heat pump constraints, electric water heater constraints and electric vehicle constraints.

[0110] Furthermore, the load response constraint formula is as follows:

[0111] Fixed frequency heat pump: (38);

[0112] Inverter heat pump:

[0113] (39);

[0114] Electric water heater:

[0115] (40);

[0116] Electric Vehicles:

[0117] (41);

[0118] In the formula, and They represent the lower and upper limits of the aggregated power of the fixed-frequency heat pump cluster, represents the aggregate power of the fixed-frequency heat pump cluster, and They represent the lower and upper limits of the aggregated power of the variable frequency heat pump cluster, represents the aggregate power of the variable frequency heat pump cluster, and They represent the lower and upper limits of the aggregated power of the electric water heater cluster, represents the aggregate power of the electric water heater cluster, and They represent the lower and upper limits of the aggregated power of the electric vehicle cluster, represents the aggregate power of the electric vehicle cluster.

[0119] Furthermore, the tie line power constraint formula is as follows:

[0120] (42);

[0121] (43);

[0122] In the formula, and They represent the lower and upper active limits of the interconnection line exchange power, represents the active power of the tie line exchange power, and They represent the reactive lower bound and reactive upper bound of the tie line exchange power, Reactive power representing the exchange power of the tie line.

[0123] Furthermore, the static VAR generator (SVG) constraint formula is as follows:

[0124] (44);

[0125] In the formula, and Indicates the lower and upper limits of the reactive power generated by the static VAR generator. Indicates the reactive power of the static VAR generator.

[0126] Furthermore, the distributed generation (DG) regulation constraint formula is as follows:

[0127] (45);

[0128] (46);

[0129] ;

[0130] In the formula, represents the reactive power of DG, and They represent the reactive power lower bound and reactive power upper bound of DG respectively, Indicates the active power of DG during this period, Indicates the active power of the next counting period, and Indicates the lower and upper limits of DG output. represents the power factor of DG, and They represent the lower and upper limits of DG's climbing ability respectively.

[0131] Furthermore, the energy storage response (ES) constraint formula is as follows:

[0132] ;

[0133] (49);

[0134] (50);

[0135] (51);

[0136] In the formula, It represents the energy of ES in this period, Indicates the energy of the next counting period of ES, Indicates ES charging power, Indicates ES charging efficiency, represents the ES discharge power, represents the ES discharge efficiency, Indicates the optimization scheduling time interval, and denote the lower and upper energy limits of ES, respectively. Indicates the upper limit of ES charging power. Indicates the upper limit of ES discharge power, Indicates the charging or discharging state of energy storage. It is a binary variable, 1 indicates charging and 0 indicates discharging.

[0137] Furthermore, the distribution network system constraint comprises the following steps:

[0138] The Distflow branch flow model is converted into a feasible solution mode by performing second-order cone relaxation. The formula is as follows:

[0139] (52);

[0140] (53);

[0141] (54);

[0142] In the formula, Representation Node , The active power of the line between Representation Node , The reactive power of the line between Indicates line , The square of the current between nodes, Representation Node The square of the voltage amplitude;

[0143] Considering the node injection power coupling constraint of load response, the formula is as follows:

[0144] (55);

[0145] (56);

[0146] In the formula, Representation Node The injected active power, Representation Node The active power of the interconnection line exchange power, Representation Node The active power of distributed generation, Representation Node The active power of photovoltaic Representation Node The active power of the fan, represents the ES discharge power, Indicates ES charging power, Representation Node The operating power of the fixed frequency heat pump (HPss) aggregated at Representation Node The aggregated variable speed heat pump (HPiv) operating power at Representation Node The operating power of the electric water heater (EWH) aggregated at Representation Node Aggregate electric vehicle (EV) operating power at Representation Node Aggregate transferable load (TFL) operating power at Representation Node Aggregate Translatable Load (TLL) operating power at Representation Node The aggregated interruptible load (IL) operating power at represents the active power of the node load, Representation Node The injected reactive power, Representation Node The reactive power of the interconnection line exchange power, Representation Node The reactive power of distributed generation, Representation Node SVG power, Represents the reactive power of the node load;

[0147] Considering the node voltage and line current constraints to ensure safe operation of the system, the formula is as follows:

[0148] (57);

[0149] (58);

[0150] In the formula, , They represent the square of the lower limit of the node voltage constraint and the square of the upper limit of the node voltage constraint, respectively. represents the square of the voltage at node j, , They represent the square of the lower limit of the line current constraint and the square of the upper limit of the line current constraint, respectively. express , The square of the line current between nodes. Compared with the prior art, the present invention has the following advantages:

[0151] The present application discloses a method for optimizing the dispatching of a distribution network with aggregatable resources as the main body. The aggregatable resources specifically include fixed-frequency heat pumps, variable-frequency heat pumps, electric water heaters, electric vehicles, transferable loads, shiftable loads and interruptible loads. The present application conducts targeted quantitative modeling of the aggregatable resources and simultaneously establishes a Distflow branch flow model of the distribution network. On the premise of comprehensively considering the renewable energy on the power generation side and the aggregatable resources on the load side, an optimization objective function is established, and multi-dimensional constraints are combined to achieve multi-objective optimization of low network loss, low operating cost and maximum absorption of renewable energy in the local power grid.

[0152] The present application discloses a Distflow branch power flow model of a distribution network, which is a multi-dimensional power flow model for radial distribution networks. It performs power flow modeling from three aspects: node power balance, node voltage correlation, and line power calculation, so as to maximize the stability and reliability of the distribution network during the dispatching and operation process.

[0153] The present application discloses an optimization objective function, which comprehensively considers the power generation side, the power grid side and the power consumption side, wherein the power generation side considers photovoltaic and wind power resources, the power grid side considers the network loss of the overall operation of the distribution network, and the power consumption side considers the operating cost of each power user. In this way, the operation and scheduling of the distribution network can be optimized in three aspects: optimizing the network loss of the system operation, the operating cost of each unit in the system operation, and maximizing the absorption of local distributed renewable energy and reducing power abandonment, and the optimization is more comprehensive.

[0154] This application discloses multi-dimensional constraints, including load response constraints, interconnection line power constraints, static VAR generator constraints, distributed power supply regulation constraints, energy storage response constraints and distribution network system operation constraints, covering a wide range. In addition, a second-order cone relaxation algorithm is used for the distribution network system operation constraints to improve the calculation accuracy and response speed of the distribution network operation constraints. BRIEF DESCRIPTION OF THE DRAWINGS

[0155] The above and / or additional aspects and advantages of the present invention will become apparent and easily understood from the description of the embodiments in conjunction with the following drawings, in which:

[0156] Figure 1 It is the overall flow chart of the present invention;

[0157] Figure 2It is a schematic diagram of power flow of the power distribution network when the present invention is implemented;

[0158] Figure 3 It is a parameter setting diagram when the present invention is implemented;

[0159] Figure 4 is the original load diagram of the transferable load, the translatable load and the interruptible load when the present invention is implemented;

[0160] Figure 5 It is a data diagram of outside temperature, renewable energy output and time-of-use electricity price when the present invention is implemented;

[0161] Figure 6 It is an operational effect diagram after the implementation of the present invention. DETAILED DESCRIPTION

[0162] Embodiments of the present invention are described in detail below, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain the present invention, and cannot be understood as limiting the present invention.

[0163] like Figure 1 As shown, this embodiment proposes a distribution network dispatch optimization method based on aggregatable resources, wherein:

[0164] Aggregable resources refer to a variety of flexible and controllable resources on the load side of the smart grid, including fixed-frequency heat pumps, variable-frequency heat pumps, electric water heaters, electric vehicles, transferable loads, shiftable loads, and interruptible loads.

[0165] The following steps are involved:

[0166] S1: Based on the operating characteristics of aggregatable resources, an aggregation power quantification model of aggregatable resources is constructed to dynamically quantify the response power of aggregatable resources and realize coordinated response among aggregatable resources;

[0167] S2: Considering the node power balance, node voltage correlation and line power, the Distflow branch power flow model of the distribution network is constructed;

[0168] S3: Consider optimizing the network loss of the system operation, optimizing the operating cost of each unit in the system operation, maximizing the absorption of local distributed renewable energy and reducing power abandonment, and establishing the optimization objective function;

[0169] S4: Aggregate constraints, solve the model constructed in S1-S3, and after calculating the results, send them to the resource entities under the distribution network for execution.

[0170] This application conducts targeted quantitative modeling of aggregatable resources and establishes a Distflow branch flow model for the distribution network. On the premise of comprehensively considering renewable energy on the power generation side and aggregatable resources on the load side, an optimization objective function is established, combined with multi-dimensional constraints, to achieve multi-objective optimization of low network losses, low operating costs and maximum absorption of renewable energy in the local power grid.

[0171] The heat pump (HP) model includes a fixed-frequency heat pump and a variable-frequency heat pump. A heat pump is a device that can obtain low-level heat energy from the air, water or soil in nature, and provide high-level heat energy that can be used by people through electrical energy. In S1, the modeling process of the heat pump model includes the following steps:

[0172] Considering the cooling service demand, the maximum operating power during the evaluation period corresponds to the power consumed by the current temperature state to reach the lower limit of the controlled temperature, as follows:

[0173] ;

[0174] In the formula, Indicates that the lower temperature limit was reached during the evaluation period Required cooling power, represents the equivalent heat capacity of the building walls and room air, Indicates the lower limit of the comfortable temperature allowed in the room. Indicates the room temperature. represents the time step, represents the radiant heat power of the room, It represents the thermal power of the heat source in the room, that is, the rate at which heat is generated in the room. Indicates the outdoor temperature, Indicates the equivalent thermal resistance of the building wall and air;

[0175] Considering the limitation of heat pump operating power, the maximum thermal power allowed in the heat pump thermal zone is calculated as follows:

[0176] (2);

[0177] In the formula, It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the upper limit of the cooling power of the heat pump equipment. Indicates that the lower temperature limit was reached during the evaluation period Required cooling power;

[0178] The minimum operating power of the heat pump during the evaluation period corresponds to the power consumed by the current temperature state to reach the controlled temperature upper limit, and the formula is as follows:

[0179] ;

[0180] In the formula, Indicates that the upper temperature limit was reached during the evaluation period Required cooling power, represents the equivalent heat capacity of the building walls and room air, Indicates the upper temperature limit, Indicates the room temperature. represents the time step, represents the radiant heat power of the room, It represents the thermal power of the heat source in the room, that is, the rate at which heat is generated in the room. Indicates the outdoor temperature, Indicates the equivalent thermal resistance of the building wall and air;

[0181] Considering that the operating power of the heat pump equipment cannot be negative, the minimum thermal power allowed in the heat pump thermal zone is calculated as follows:

[0182] ;

[0183] In the formula, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. Indicates that the upper temperature limit was reached during the evaluation period Required cooling power;

[0184] Calculate the aggregate flexible power range of the fixed-frequency heat pump cluster during the evaluation period using the following formula:

[0185] ;

[0186] In the formula, Indicates the lower limit of the aggregate power of the fixed-frequency heat pump cluster. Indicates the aggregate power upper limit of the fixed-frequency heat pump cluster. represents the scale of the fixed-frequency heat pump cluster, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the operating efficiency coefficient of the fixed-frequency heat pump in the room;

[0187] Calculate the aggregate flexible power range of the variable frequency heat pump cluster during the evaluation period using the following formula:

[0188] ;

[0189] In the formula, Indicates the lower limit of the aggregate power of the variable frequency heat pump cluster; Indicates the aggregate power upper limit of the variable frequency heat pump cluster; represents the size of the variable frequency heat pump cluster, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the operating efficiency coefficient of the variable frequency heat pump in the room.

[0190] In this model, the calculation is carried out based on the cooling demand, and the state limit value under the heating demand is the opposite. When calculating the aggregate flexibility of cluster-scale HP, the adjustment power of HPss is treated as continuous here, and its discrete variables can be continuous when analyzing the response capacity of large-scale HPss.

[0191] An electric water heater (EWH) is a water heater that uses electricity as energy for heating. In S1, the modeling process of the electric water heater model includes the following steps:

[0192] During the evaluation period, at the current hot water temperature, the hot water temperature cannot exceed the upper limit of the water temperature under the maximum operating power of the electric water heater. Taking into account the power limit of the equipment, the maximum thermal power allowed by the electric water heater is calculated. The formula is as follows:

[0193] ;

[0194] In the formula, Indicates that the water temperature of the electric water heater reaches The average operating power required, Indicates the upper limit of water temperature. Indicates the water temperature of the electric water heater. Indicates the heat capacity of the electric water heater. represents the time step, represents the thermal disturbance power caused by hot water use, Indicates the indoor temperature. represents the thermal resistance of the water heater, Indicates the operating efficiency of the electric water heater. It represents the maximum operating power allowed for electric water heater unit i during the evaluation period. Indicates the upper limit of the operating power of the electric water heater;

[0195] During the evaluation period, at the current hot water temperature, the minimum operating power of the electric water heater must meet the user's hot water demand at the set time, while taking into account the equipment operating power limit, and calculate the minimum thermal power allowed by the electric water heater. The formula is as follows:

[0196] ;

[0197] (13);

[0198] In the formula, Indicates that the water temperature of the electric water heater reaches the set time In the corresponding period The average operating power within represents the time step, Indicates that the water temperature of the electric water heater reaches the set time In the corresponding period The average operating power within Indicates the desired water temperature. Indicates the water temperature of the electric water heater. Indicates the heat capacity of the electric water heater. Indicates the set time for the water heater to work. represents the thermal disturbance power caused by hot water use, Indicates the indoor temperature. represents the thermal resistance of the water heater, Indicates the operating efficiency of the electric water heater. Indicates the upper limit of the operating power of the electric water heater. represents the minimum operating power of electric water heater unit i during the evaluation period;

[0199] The aggregate flexible power range of the electric water heater cluster during the evaluation period is calculated as follows:

[0200] (14);

[0201] (15);

[0202] In the formula, Indicates the lower limit of the aggregate power of the electric water heater cluster. Indicates the aggregate power upper limit of the electric water heater cluster. represents the size of the electric water heater cluster, represents the minimum operating power of electric water heater unit i during the evaluation period, It represents the maximum allowed operating power of electric water heater unit i during the evaluation period.

[0203] Electric vehicles (EVs) are vehicles that use an onboard power source as a power source, use an electric motor to drive the wheels, and comply with various requirements of road traffic and safety regulations. In S1, the modeling process of the electric vehicle model includes the following steps:

[0204] During the evaluation period, under the current SOC of the electric vehicle, the maximum charging power of the electric vehicle meets the SOC constraint. Considering the power limit of the charging equipment, the maximum charging power of the electric vehicle is calculated. The formula is as follows:

[0205] ;

[0206] In the formula, Indicates that the SOC of the electric vehicle reaches The average charging power required, represents the upper limit of charging power of electric vehicle cell i, It represents the maximum charging power allowed for electric vehicle single i during the evaluation period, Indicates the maximum value of the battery state of charge. Indicates the battery charge state of an electric vehicle. Indicates the battery capacity. represents the charging efficiency of electric vehicles, represents the time step;

[0207] During the evaluation period, under the current SOC of the electric vehicle, the maximum discharge power of the electric vehicle during the evaluation period should meet the SOC requirement of the user at the time of departure. The maximum discharge power allowed by the electric vehicle is calculated as follows:

[0208] ;

[0209] (19);

[0210] In the formula, Indicates that the SOC of electric vehicles can meet the requirements in the subsequent period The average discharge power allowed during the evaluation period is: represents the upper limit of charging power of electric vehicle cell i, represents the charging efficiency of electric vehicles, Indicates the time when the electric vehicle leaves. represents the time step, represents the expected state of charge of the electric vehicle when it drives away, Indicates the battery charge state of an electric vehicle. Indicates the battery capacity. represents the discharge efficiency of electric vehicles, It represents the maximum discharge power allowed by electric vehicle cell i during the evaluation period, It represents the upper limit of discharge power of electric vehicle cell i;

[0211] The aggregated flexible power range of the electric vehicle cluster during the evaluation period is calculated as follows:

[0212] (20);

[0213] (twenty one);

[0214] In the formula, represents the lower limit of the aggregated power of the electric vehicle cluster, represents the aggregate power cap of the electric vehicle cluster, represents the size of the electric vehicle cluster, It represents the maximum discharge power allowed by electric vehicle cell i during the evaluation period, The maximum charging power allowed for electric vehicle cell i during the assessment period.

[0215] Shiftable load (TFL) refers to the load that users can change during their electricity consumption period according to electricity prices or incentives.

[0216] This type of load is usually sensitive to electricity prices or incentives, and can achieve load transfer during this period. Transferable load is a significant means of load management, which is not only conducive to the economic operation of microgrids, but also can regulate the elasticity of load demand and effectively achieve the purpose of peak shaving and valley filling.

[0217] The transferable load can be reasonably modified and adjusted according to the actual situation. At the same time, the load compensation costs required for transferring loads of different time periods and power sizes must be considered.

[0218] In S1, the modeling process of transferable loads includes the following steps:

[0219] Taking into account the flexibility of system operation regulation, the transferable load can flexibly deploy the load power within the allowed time period, and the operating power regulation equation formula is as follows:

[0220] (twenty two);

[0221] ;

[0222] In the formula, is the load after transfer, is the load before transfer, For transferable time periods, is a binary variable, if , ,otherwise, , is a set of scheduling time periods, and is the adjustable range of TFL.

[0223] The controllability of the shiftable load (TLL) is very strong, and the load power supply time can be changed according to the plan (the total power supply is unchanged, the power supply period is variable, the compensation is very small, or even no compensation is required from the power supply department). It is different from the interruptible load (using the interruptible electricity price management, directly cutting off the load and providing certain compensation). This type of load mainly includes washing machines, disinfection cabinets, etc.

[0224] In S1, the modeling process of the translatable load includes the following steps:

[0225] Considering the industrial production process, the shiftable load provides operation flexibility by shifting the load curve. The operating power adjustment equation is as follows:

[0226] Calculate the load after translation using the following formula:

[0227] (twenty four);

[0228] In the formula, represents the load after translation, represents the decision vector for translation, represents the translation scheme matrix;

[0229] Decision vector for translation The calculation formula is as follows:

[0230] (25);

[0231] (26);

[0232] In the formula, is a binary variable, Indicates that the corresponding translation scheme is selected, otherwise it is not selected;

[0233] Translation Matrix The calculation formula is as follows:

[0234] (27);

[0235] In the formula, The number of translation options for the translatable load, Indicates different scheduling time periods.

[0236] Interruptible load (IL) refers to the part of the user load center that can be interrupted during the peak period of the power grid or in an emergency. In S1, the modeling process of interruptible load includes the following steps:

[0237] By giving economic compensation, IL can be interrupted or reduced according to the user's wishes, and its operating power regulation equation is as follows:

[0238] ;

[0239] In the formula, represents the load before the interruption, represents the load after the interruption, and This is the range in which interrupts are allowed.

[0240] In S2, Figure 2As shown in the figure, the modeling process of the distribution network Distflow branch flow model includes the following steps:

[0241] Establish a radial distribution network model;

[0242] Considering the node power balance, node voltage correlation and line power, the Distflow branch power flow model of the distribution network is constructed. The node power balance equation is as follows:

[0243] ;

[0244] In the formula, the subscript They all represent line nodes. Representation Node The injected active power, Representation Node The injected reactive power, Representation Node , The active power of the line between Representation Node , The active power of the line between Representation Node , The reactive power of the line between Representation Node , The reactive power of the line between Representation Node , The square of the line current amplitude between Representation Node , The line resistance between Representation Node , The line reactance between

[0245] The node voltage association equation formula is as follows:

[0246] ;

[0247] In the formula, Representation Node The square of the voltage amplitude, Representation Node The square of the voltage amplitude, Represents the line resistance, Representation Node , The active power of the line between Representation Node , The line reactance between Representation Node , The reactive power of the line between Representation Node , The square of the line current amplitude between them;

[0248] The line power equation formula is as follows:

[0249] (32);

[0250] In the formula, Representation Node , The active power of the line between Representation Node , The reactive power of the line between Representation Node The square of the voltage amplitude, Representation Node , The square of the line current amplitude between them.

[0251] The Distflow branch power flow model is a multi-dimensional power flow model for radial distribution networks. It models power flow from three aspects: node power balance, node voltage correlation and line power calculation, so as to maximize the stability and reliability of the distribution network during the dispatching and operation process.

[0252] In S3, the calculation process of optimizing the objective function includes the following steps:

[0253] The network loss cost of the distribution network is calculated using the following formula:

[0254] ;

[0255] In the formula, represents the network loss cost target item, represents the network loss cost coefficient of the distribution network, represents the calculation time period, represents the set of distribution network lines, Representation Node , The square of the line current amplitude between Representation Node , The line resistance between Indicates the optimization scheduling time interval;

[0256] The operating cost of the distribution network is calculated as follows:

[0257] (34);

[0258] In the formula, represents the operating cost target item, represents the tie line power cost coefficient, is the tie line power, represents the static VAR generator (SVG) cost factor, Indicates SVG power, represents the cost coefficient of controllable distributed generator (DG), represents the DG power, represents the energy storage (ES) cost coefficient, Indicates ES charging power, represents the ES discharge power, represents the load response cost coefficient, Indicates load response regulation power, Indicates the optimization scheduling time interval, Indicates the calculation time period;

[0259] The calculation formula is as follows:

[0260] (35);

[0261] In the formula, Indicates the aggregated HPss operating power, Indicates the load before the HPss class load is controlled. represents the aggregated HPiv operating power, represents the load before the aggregated HPiv class load is controlled, represents the aggregated EWH operating power, represents the load before the aggregated EWH load is controlled, represents the aggregated EV operating power, Indicates the load before the EV load is controlled. represents the aggregated TFL operating power, represents the load before the TFL type load is controlled, represents the aggregated TLL operating power, Indicates the load before the TLL load is controlled. represents the aggregated IL operating power, It represents the load before the IL type load is controlled;

[0262] The formula for calculating the RE curtailment cost of the distribution network is as follows:

[0263] (36);

[0264] In the formula, represents the target item of power curtailment cost, represents the calculation time period, represents the cost coefficient of RE power abandonment in the distribution network, represents the abandoned power of photovoltaic (PV), represents the abandoned power of the wind turbine (WT), Indicates the optimization scheduling interval.

[0265] Calculate the optimal scheduling target, the formula is as follows

[0266] (37);

[0267] In the formula, represents the optimization scheduling target item, represents the network loss cost target item, represents the operating cost target item, Represents the curtailment cost objective item.

[0268] It should be specially noted that , , as well as The data is calculated periodically and updated in real time. , , as well as The data will be saved in sequence to form a set of optimization scheduling targets to be selected. When the optimization scheduling target needs to be determined, the smallest one is selected from the currently saved optimization scheduling target set. The value is used as the final optimization dispatch target and sent to the resource entities under the distribution network for execution.

[0269] The optimization objective function comprehensively considers the power generation side, the grid side and the power consumption side. The power generation side considers photovoltaic and wind power resources, the grid side considers the network loss of the overall operation of the distribution network, and the power consumption side considers the operating cost of each power user. It can optimize the distribution network operation and scheduling in three aspects: optimizing the network loss of system operation, the operating cost of each unit in the system operation, and maximizing the absorption of local distributed renewable energy and reducing power abandonment. The optimization is more comprehensive.

[0270] In S4, the constraints consist of the following parts: load response constraints, tie line power constraints, static VAR generator constraints, distributed power generation regulation constraints, energy storage response constraints and distribution network system operation constraints, a total of six parts. Among them, the load response constraints include fixed-frequency heat pump constraints, variable-frequency heat pump constraints, electric water heater constraints and electric vehicle constraints.

[0271] Specifically, the load response constraint formula is as follows:

[0272] Fixed frequency heat pump:

[0273] (38);

[0274] Inverter heat pump:

[0275] (39);

[0276] Electric water heater:

[0277] (40);

[0278] Electric Vehicles:

[0279] (41);

[0280] In the formula, and They represent the lower and upper limits of the aggregated power of the fixed-frequency heat pump cluster, represents the aggregate power of the fixed-frequency heat pump cluster, and They represent the lower and upper limits of the aggregated power of the variable frequency heat pump cluster, represents the aggregate power of the variable frequency heat pump cluster, and They represent the lower and upper limits of the aggregated power of the electric water heater cluster, represents the aggregate power of the electric water heater cluster, and They represent the lower and upper limits of the aggregated power of the electric vehicle cluster, represents the aggregate power of the electric vehicle cluster.

[0281] The tie line power constraint formula is as follows:

[0282] (42);

[0283] (43);

[0284] In the formula, and They represent the lower and upper active limits of the interconnection line exchange power, represents the active power of the tie line exchange power, and They represent the reactive lower bound and reactive upper bound of the tie line exchange power, Reactive power representing the exchange power of the tie line.

[0285] The static VAR generator (SVG) constraint formula is as follows:

[0286] (44);

[0287] In the formula, and Indicates the lower and upper limits of the reactive power generated by the static VAR generator. Indicates the reactive power of the static VAR generator.

[0288] The distributed generation (DG) regulation constraint formula is as follows:

[0289] (45);

[0290] (46);

[0291] ;

[0292] In the formula, represents the reactive power of DG, and They represent the reactive power lower bound and reactive power upper bound of DG respectively, Indicates the active power of DG during this period, Indicates the active power of the next counting period, and Indicates the lower and upper limits of DG output. represents the power factor of DG, and They represent the lower and upper limits of DG's climbing ability respectively.

[0293] The energy storage response (ES) constraint formula is as follows:

[0294] ;

[0295] (49);

[0296] (50);

[0297] (51);

[0298] In the formula, It represents the energy of ES in this period, Indicates the energy of the next counting period of ES, Indicates ES charging power, Indicates ES charging efficiency, represents the ES discharge power, represents the ES discharge efficiency, Indicates the optimization scheduling time interval, and denote the lower and upper energy limits of ES, respectively. Indicates the upper limit of ES charging power. Indicates the upper limit of ES discharge power, Indicates the charging or discharging state of energy storage. It is a binary variable, 1 indicates charging and 0 indicates discharging.

[0299] The distribution network system constraints include the following steps:

[0300] The Distflow branch flow model is converted into a feasible solution mode by performing second-order cone relaxation. The formula is as follows:

[0301] (52);

[0302] (53);

[0303] (54);

[0304] Considering the node injection power coupling constraint of load response, the formula is as follows:

[0305] (55);

[0306] (56);

[0307] In the formula, Representation Node The injected active power, Representation Node The active power of the interconnection line exchange power, Representation Node The active power of distributed generation, Representation Node The active power of photovoltaic Representation Node The active power of the fan, represents the ES discharge power, Indicates ES charging power, Representation Node The operating power of the fixed frequency heat pump (HPss) aggregated at Representation Node The aggregated variable speed heat pump (HPiv) operating power at Representation Node The operating power of the electric water heater (EWH) aggregated at Representation Node Aggregated electric vehicle (EV) operating power, Representation Node Aggregate transferable load (TFL) operating power at Representation Node Aggregate Translatable Load (TLL) operating power at Representation Node The aggregated interruptible load (IL) operating power at represents the active power of the node load, Representation Node The injected reactive power, Representation Node The reactive power of the interconnection line exchange power, Representation Node The reactive power of distributed generation, Indicates the SVG power at the node, Indicates the reactive power of the node load.

[0308] Considering the node voltage and line current constraints to ensure safe operation of the system, the formula is as follows:

[0309] (57);

[0310] (58);

[0311] In the formula, , They represent the square of the lower limit of the node voltage constraint and the square of the upper limit of the node voltage constraint, respectively. represents the square of the voltage at node j, , They represent the square of the lower limit of the line current constraint and the square of the upper limit of the line current constraint, respectively. express , The square of the line current between nodes.

[0312] In this application, multi-dimensional constraints are set, including load response constraints, interconnection line power constraints, static VAR generator constraints, distributed power supply regulation constraints, energy storage response constraints and distribution network system operation constraints, covering a wide range. In addition, a second-order cone relaxation algorithm is used for the distribution network system operation constraints to improve the calculation accuracy and response speed of the distribution network operation constraints.

[0313] In S4, the model constructed in steps S1-S3 is solved, and after the results are calculated, they are sent to the resource entities under the distribution network for execution, thereby achieving multi-objective optimization of low network losses, low operating costs and maximum absorption of renewable energy in the local power grid.

[0314] The present invention also provides a system for optimizing the dispatching of distribution networks taking into account the response of aggregatable resources. The system implements a method for optimizing the dispatching of distribution networks taking into account the response of aggregatable resources. Specifically, the system includes a modeling module, a building module and an aggregation module. The modeling module is used to complete the modeling of the distribution network, the building module is used to complete the construction of the objective function, and the aggregation module is used to complete the aggregation of constraint conditions. The modeling module, the building module and the aggregation module are integrated into a local optimization software platform, which can interact with the superior dispatching organization for data and respond to instructions.

[0315] When the distribution network has a dispatching demand, the dispatching instruction is first sent to the local optimization software platform. The software platform calls three functional modules to perform calculations and solutions. After giving the calculation results, it is sent to the resource entity under the jurisdiction of the distribution network for execution.

[0316] It should be noted that in the process of calculation and solution, within the same calculation cycle, the modeling module is started first and then the aggregation module is started for calculation. There is no requirement for the order of the construction module and other modules, and they can be started before or after, earlier or later.

[0317] In addition, the above calculation process is a dynamic optimization process, that is, a specific duration can be set as a time period, and the parameters in each time period change dynamically, but the solution model and algorithm are the same; and the n+1th calculation period will use the data generated by the nth calculation period to update the parameters, that is, there is a correlation between the previous and next calculation periods, rather than being independent of each other, in order to achieve the optimal effect.

[0318] Please refer to the attached Figure 3 To Attachment Figure 5 In the specific implementation, an example of IEEE 33-node power distribution system was constructed for analysis. Various types of resource elements were designed in the system environment, including residential HPss, residential EWH, commercial HPiv, EV charged at commercial locations (EVc), and EV charged at residential locations (EVr). EVc and EVr have the same parameters, but their driving-in and driving-out rate characteristics are different. The charging power of electric vehicles in commercial areas is designed to be 30kW, corresponding to fast charging. The residential area is designed to be 10kW, corresponding to slow charging. The thermal parameters of HPss, HPiv and EWH all obey normal distribution. The temperature dead zone of HP is set to , the allowable water temperature of EWH is , the original load curve of the large-capacity load elastic resource is as follows Figure 4 The ambient temperature, the total predicted output of photovoltaic (PV) and wind turbine (WT), and the time-of-use electricity price of the tie line power are shown in Figure 5 As shown in the figure, the resource response cost on the load side is set to 40% of the grid electricity price. TLL, IL and TFL are distributed at nodes 6, 7 and 29, EVr and EVc are distributed at nodes 23 and 24, EWH is distributed at node 7, and HPss and HPiv are distributed at nodes 30 and 31.

[0319] According to the attached Figure 6 It can be seen that after the optimization of the distribution network dispatching optimization method based on aggregatable resources of the present invention, the operation cost, network loss, RE absorption rate and RE power reduction of the distribution network have been well improved, and the optimization operation level of the distribution network has been greatly improved. Compared with before optimization, the system operation cost has been reduced by 7.0%, the network loss of the distribution network has been reduced by 9.3%, RE has achieved a high proportion of absorption, and the absorption rate has been increased by 8.7%.

[0320] Although the embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the claims and their equivalents.

Claims

1. A distribution network dispatch optimization method based on aggregatable resources, characterized in that: The following steps are involved: S1: Based on the operating characteristics of aggregatable resources, an aggregation power quantification model of aggregatable resources is constructed to dynamically quantify the response power of aggregatable resources and realize coordinated response among aggregatable resources; S2: Considering the node power balance, node voltage correlation and line power, the Distflow branch power flow model of the distribution network is constructed; S3: Consider optimizing the network loss of the system operation, optimizing the operating cost of each unit in the system operation, maximizing the absorption of local distributed renewable energy and reducing power abandonment, and establishing the optimization objective function; S4: Aggregate constraints, solve the model constructed in S1-S3, and after calculating the results, send them to the resource entities under the distribution network for execution; Aggregable resources include heat pumps, electric water heaters, electric vehicles, transferable loads, shiftable loads, and interruptible loads; The modeling process of the heat pump model includes the following steps: Considering the cooling service demand, the maximum operating power during the evaluation period corresponds to the power consumed by the current temperature state to reach the lower limit of the controlled temperature, as follows: ； In the formula, Indicates that the lower temperature limit was reached during the evaluation period Required cooling power, represents the equivalent heat capacity of the building walls and room air, Indicates the lower limit of the comfortable temperature allowed in the room. Indicates the room temperature. represents the time step, represents the radiant heat power of the room, Indicates the thermal power of the heat source in the room, Indicates the outdoor temperature, Indicates the equivalent thermal resistance of the building wall and air; Considering the limitation of heat pump operating power, the maximum thermal power allowed in the heat pump thermal zone is calculated as follows: （2）； In the formula, It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the upper limit of the cooling power of the heat pump equipment. Indicates that the lower temperature limit was reached during the evaluation period Required cooling power; The minimum operating power of the heat pump during the evaluation period corresponds to the power consumed by the current temperature state to reach the controlled temperature upper limit, and the formula is as follows: ； In the formula, Indicates that the upper temperature limit was reached during the evaluation period Required cooling power, represents the equivalent heat capacity of the building walls and room air, Indicates the upper limit of the comfortable temperature allowed in the room. Indicates the room temperature. represents the time step, represents the radiant heat power of the room, Indicates the thermal power of the heat source in the room, Indicates the outdoor temperature, Indicates the equivalent thermal resistance of the building wall and air; Considering that the operating power of the heat pump equipment cannot be negative, the minimum thermal power allowed in the heat pump thermal zone is calculated as follows: ； In the formula, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. Indicates that the upper temperature limit was reached during the evaluation period Required cooling power; Calculate the aggregate flexible power range of the fixed-frequency heat pump cluster during the evaluation period using the following formula: ； In the formula, Indicates the lower limit of the aggregate power of the fixed-frequency heat pump cluster. Indicates the aggregate power upper limit of the fixed-frequency heat pump cluster. represents the scale of the fixed-frequency heat pump cluster, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the operating efficiency coefficient of the fixed-frequency heat pump in the room; Calculate the aggregate flexible power range of the variable frequency heat pump cluster during the evaluation period using the following formula: ； In the formula, Indicates the lower limit of the aggregate power of the variable frequency heat pump cluster; Indicates the aggregate power upper limit of the variable frequency heat pump cluster; represents the size of the variable frequency heat pump cluster, It represents the minimum cooling power of the house when the upper limit of the allowable comfort temperature is met. It indicates the maximum cooling power of the house when the lower limit of the allowable comfort temperature is met. Indicates the operating efficiency coefficient of the variable frequency heat pump in the room; The modeling process of the electric water heater model includes the following steps: During the evaluation period, at the current hot water temperature, the hot water temperature cannot exceed the upper limit of the water temperature under the maximum operating power of the electric water heater. Taking into account the power limit of the equipment, the maximum thermal power allowed by the electric water heater is calculated. The formula is as follows: ； In the formula, Indicates that the water temperature of the electric water heater reaches The average operating power required, Indicates the upper limit of water temperature. Indicates the water temperature of the electric water heater. Indicates the heat capacity of the electric water heater. represents the time step, represents the thermal disturbance power caused by hot water use, Indicates the indoor temperature. represents the thermal resistance of the water heater, Indicates the operating efficiency of the electric water heater. It represents the maximum operating power allowed for electric water heater unit i during the evaluation period. Indicates the upper limit of the operating power of the electric water heater; During the evaluation period, at the current hot water temperature, the minimum operating power of the electric water heater must meet the user's hot water demand at the set time, while taking into account the equipment operating power limit, and calculate the minimum thermal power allowed by the electric water heater. The formula is as follows: ； （13）； In the formula, Indicates that the water temperature of the electric water heater reaches the set time In the corresponding period The average operating power within represents the time step, Indicates that the water temperature of the electric water heater reaches the set time In the corresponding period The average operating power within Indicates the desired water temperature. Indicates the water temperature of the electric water heater. Indicates the heat capacity of the electric water heater. Indicates the set time for the water heater to work. represents the thermal disturbance power caused by hot water use, Indicates the indoor temperature. represents the thermal resistance of the water heater, Indicates the operating efficiency of the electric water heater. Indicates the upper limit of the operating power of the electric water heater. represents the minimum operating power of electric water heater unit i during the evaluation period; The aggregate flexible power range of the electric water heater cluster during the evaluation period is calculated as follows: （14）； （15）； In the formula, Indicates the lower limit of the aggregate power of the electric water heater cluster. Indicates the aggregate power upper limit of the electric water heater cluster. represents the size of the electric water heater cluster, represents the minimum operating power of electric water heater unit i during the evaluation period, It represents the maximum operating power allowed for electric water heater unit i during the evaluation period; The modeling process of the electric vehicle model includes the following steps: During the evaluation period, under the current SOC of the electric vehicle, the maximum charging power of the electric vehicle meets the SOC constraint. Considering the power limit of the charging equipment, the maximum charging power of the electric vehicle is calculated. The formula is as follows: ； In the formula, Indicates that the SOC of the electric vehicle reaches The average charging power required, represents the upper limit of charging power of electric vehicle cell i, It represents the maximum charging power allowed for electric vehicle single i during the evaluation period, Indicates the maximum value of the battery state of charge. Indicates the battery charge state of an electric vehicle. Indicates the battery capacity. represents the charging efficiency of electric vehicles, represents the time step; During the evaluation period, under the current SOC of the electric vehicle, the maximum discharge power of the electric vehicle during the evaluation period should meet the SOC requirement of the user at the time of departure. The maximum discharge power allowed by the electric vehicle is calculated as follows: ； （19）； In the formula, Indicates that the SOC of the electric vehicle can meet the requirements in the subsequent period The average discharge power allowed during the evaluation period is: represents the upper limit of charging power of electric vehicle cell i, represents the charging efficiency of electric vehicles, Indicates the time when the electric vehicle leaves. represents the time step, represents the expected state of charge of the electric vehicle when it drives away, Indicates the battery charge state of an electric vehicle. Indicates the battery capacity. represents the discharge efficiency of electric vehicles, It represents the maximum discharge power allowed by electric vehicle cell i during the evaluation period, It represents the upper limit of discharge power of electric vehicle cell i; The aggregated flexible power range of the electric vehicle cluster during the evaluation period is calculated as follows: （20）； （21）； In the formula, represents the lower limit of the aggregated power of the electric vehicle cluster, represents the aggregate power cap of the electric vehicle cluster, represents the size of the electric vehicle cluster, It represents the maximum discharge power allowed by electric vehicle cell i during the evaluation period, The maximum charging power allowed for electric vehicle unit i during the assessment period; The process of modeling transferable loads includes the following steps: Taking into account the flexibility of system operation regulation, the transferable load can flexibly deploy the load power within the allowed time period, and the operating power regulation equation formula is as follows: （22）； ； In the formula, is the load after transfer, is the load before transfer, For transferable time periods, is a binary variable, if , ,otherwise, , is a set of scheduling time periods, and are the lower and upper limits of the adjustable range of TFL; The process of modeling a translatable load consists of the following steps: Considering the industrial production process, the shiftable load provides operation flexibility by shifting the load curve. The operating power adjustment equation is as follows: Calculate the load after translation using the following formula: （24）； In the formula, represents the load after translation, represents the decision vector for translation, represents the translation scheme matrix; Decision vector for translation The calculation formula is as follows: （25）； （26）； In the formula, is a binary variable, Indicates that the corresponding translation scheme is selected, otherwise it is not selected; Translation Matrix The calculation formula is as follows: （27）； In the formula, The number of translation options for the translatable load, Indicates different scheduling time periods; The process of modeling interruptible loads includes the following steps: By giving economic compensation, the interruptible load can be interrupted or reduced according to the user's wishes, and its operating power adjustment equation is as follows: ； In the formula, represents the load before the interruption, represents the load after the interruption, and are the lower and upper limits of the range where interruption is allowed.

2. A distribution network dispatch optimization method based on aggregatable resources according to claim 1, characterized in that: In S2, the modeling process of the distribution network Distflow branch power flow model includes the following steps: Establish a radial distribution network model; Considering the node power balance, node voltage correlation and line power, the Distflow branch power flow model of the distribution network is constructed. The node power balance equation is as follows: ； In the formula, the subscript They all represent line nodes. Representation Node The injected active power, Representation Node The injected reactive power, Representation Node , The active power of the line between Representation Node , The active power of the line between Representation Node , The reactive power of the line between Representation Node , The reactive power of the line between Representation Node , The square of the line current amplitude between Representation Node , The line resistance between Representation Node , The line reactance between Calculate the node voltage correlation, the formula is as follows: ； In the formula, Representation Node The square of the voltage amplitude, Representation Node The square of the voltage amplitude, Represents the line resistance, Representation Node , The active power of the line between Representation Node , The line reactance between Representation Node , The reactive power of the line between Representation Node , The square of the line current amplitude between them; Calculate the line power using the following formula: （32）； In the formula, Representation Node , The active power of the line between Representation Node , The reactive power of the line between Representation Node The square of the voltage amplitude, Representation Node , The square of the line current amplitude between them.

3. A distribution network dispatch optimization method based on aggregatable resources according to claim 1, characterized in that: In S3, the calculation process of optimizing the objective function includes the following steps: The network loss cost of the distribution network is calculated using the following formula: ； In the formula, represents the network loss cost target item, represents the network loss cost coefficient of the distribution network, represents the calculation time period, represents the set of distribution network lines, Representation Node , The square of the line current amplitude between Representation Node , The line resistance between Indicates the optimization scheduling time interval; The operating cost of the distribution network is calculated as follows: （34）； In the formula, represents the operating cost target item, represents the tie line power cost coefficient, is the tie line power, represents the static VAR generator (SVG) cost factor, Indicates SVG power, represents the cost coefficient of controllable distributed generator (DG), represents the DG power, represents the energy storage (ES) cost coefficient, Indicates ES charging power, represents the ES discharge power, represents the load response cost coefficient, Indicates load response regulation power, Indicates the optimization scheduling time interval, Indicates the calculation time period; Load response power regulation The calculation formula is as follows: （35）； In the formula, represents the aggregated constant-frequency heat pump (HPss) operating power, Indicates the load before the constant frequency heat pump (HPss) load is controlled. represents the aggregated variable speed heat pump (HPiv) operating power, represents the load before the aggregated variable speed heat pump (HPiv) type load is controlled, represents the aggregated electric water heater (EWH) operating power, represents the load before the aggregated electric water heater (EWH) type load is controlled, represents the aggregated electric vehicle (EV) operating power, represents the load before the electric vehicle (EV) load is controlled, represents the aggregated transferable load (TFL) operating power, represents the load before the transferable load (TFL) type load is controlled, represents the aggregate Translatable Load (TLL) operating power, represents the load before the Translatable Load (TLL) type load is controlled, represents the aggregated interruptible load (IL) operating power, It indicates the load before the interruptible load (IL) type load is controlled; The formula for calculating the RE curtailment cost of the distribution network is as follows: （36）； In the formula, represents the target item of power curtailment cost, represents the calculation time period, represents the cost coefficient of RE power abandonment in the distribution network, represents the abandoned power of photovoltaic (PV), Indicates the abandoned power of the wind turbine (WT), Indicates the optimization scheduling time interval; Calculate the optimal scheduling target, the formula is as follows: （37）； In the formula, represents the optimization scheduling target item, represents the network loss cost target item, represents the operating cost target item, Represents the curtailment cost objective item.

Images