Incremental model establishment method capable of regulating and controlling electric power resources

Abstract

The invention provides an incremental model establishment method capable of regulating and controlling electric power resources, and belongs to the field of electric power regulation and control. The problem of low power resource regulation and control efficiency is solved; the method specifically comprises the following steps: calculating expected power consumption of each municipal district; constructing a discharge characteristic function of the energy storage device, and adjusting the discharge time of each energy storage device according to the discharge characteristic function; calculating expected electricity utilization power of the electricity utilization area, and constructing an admittance matrix; obtaining the actual power utilization power of the power utilization area, correcting the admittance matrix according to the expected power utilization power, and adjusting the power supply voltage, phase and power of each municipal district to the corresponding power utilization area; according to the invention, by obtaining, analyzing and mathematically modeling the related electrical data of the target area, the power and voltage changes of each power utilization area of the target area from power transmission to power utilization are simulated, and the power resource regulation and control of the target area are facilitated.

Description

Technical Field

[0001] The present invention discloses a method for establishing an incremental model of controllable electric power resources, and relates to the field of electric power control. Background Art

[0002] The existing methods for establishing power resource increment models have the following deficiencies: Difficulty in model establishment: Existing methods for establishing incremental power resource models mostly involve the deep integration of multiple disciplines such as power electronics, communications technology, computer science, and economics. For example, constructing the discharge characteristic function of an energy storage device requires combining electrochemical models with thermodynamic analysis, while the optimization of the admittance matrix relies on the analysis of complex grid topologies. This interdisciplinary nature leads to an extremely high threshold for technical implementation, requiring the R&D team to have a diverse professional background, complex model establishment, and difficult data processing.

[0003] Mismatch between static topology and dynamic grids: Most existing methods for establishing incremental power resource models do not integrate dynamic topology identification modules, resulting in the inability to update the admittance matrix in real time during grid reconstruction (such as line failures and load switching). When a transmission line is disconnected due to a fault, but the model still assumes its existence, power may be misallocated.

[0004] Ignoring the impact of distributed energy: Most existing methods for establishing incremental power resource models only consider the transmission relationship from the national grid to the municipal grid, without considering the actual power resource allocation issues of the municipal grid; at the same time, existing methods have not considered the impact of distributed power resources on the overall power supply relationship in the target area; in the scenario of high penetration of new energy, traditional control strategies may not be able to cope with sudden changes in output, resulting in frequency or voltage exceeding the limit. Summary of the Invention

[0005] In view of the shortcomings of the existing technology, the purpose of the present invention is to provide a method for establishing an incremental model of controllable power resources, aiming to solve the problem of low efficiency in power resource control.

[0006] In order to achieve the above-mentioned object, the present invention is implemented through the following technical solutions: A method for establishing an incremental model of controllable power resources includes: Obtain the number of districts in the target area, obtain the historical power consumption of each district, perform time series analysis on the historical power consumption, and calculate the expected power consumption of each district; Obtain the number of energy storage devices in the target area and construct a discharge characteristic function for each energy storage device; obtain the transmission power of the target area's power supply terminal and the remaining power of each energy storage device, and determine whether the transmission power is sufficient based on the expected power consumption of each municipal district; if sufficient, adjust the discharge time and discharge duration of each energy storage device according to the discharge characteristic function; if not sufficient, adjust the transmission power; A time series analysis is performed on the historical power consumption of each municipal district to the application power area to obtain the expected power consumption of each municipal district to the application power area, and an admittance matrix is ​​constructed; the actual power consumption of each municipal district to the application power area is obtained in real time and compared with the expected power consumption. According to the relationship between the actual power consumption and the expected power consumption, the admittance matrix is ​​corrected, and the power supply voltage, phase and power of each municipal district to the application power area are adjusted.

[0007] Furthermore, the specific steps for constructing the discharge characteristic function of each energy storage device are as follows: Count the number of energy storage devices ac and construct the discharge characteristic function Qe of the first energy storage device （1） ; Get the average discharge power pd of the first energy storage device from the first to the moth day in the past month （1） ~pd （mo） ;Discharge time td （1） 、td （2） ~td （mo） ; Calories qd （1） ~qd （mo） ;Construct matrix X and matrix Y; Let function Qe （1） The coefficient is β （0） ~β （5） , construct matrix Z; Calculating β （0） ~β （5） Value: ; Assume the first discharge power is pe, discharge time is te, define the function Qe （1） The mathematical expression is: ; δ represents the correction value; PD （1） ~pd （mo） As pe,td （1） ~td （mo） As te,qd （1） ~qd （mo） Substitute into the function Qe （1） Calculate δ from the value of Construct the discharge characteristic function of the 2nd to acth energy storage devices and obtain the function Qe （2） ~Function Qe （ac） .

[0008] Furthermore, the specific steps of adjusting the discharge time of each energy storage device according to the discharge characteristic function and the discharge time are as follows: Obtain and calculate the expected electricity consumption and ael of all city districts; Obtain and calculate the sum of the remaining power of all energy storage devices; Obtain the transmission power Pu of the power supply end; judge whether it is established; If established, no processing will be done; If not, compare and The size of the energy storage device is used to adjust the discharge power and discharge time or the transmission power and reactive power at the power supply end; like ≥ , then calculate the electric energy replenishment amount ree of the energy storage device; according to ree and the discharge characteristic function of each energy storage device, adjust the discharge power and discharge time of the energy storage device; Calculate the discharge capacity ec of the 1st to acth energy storage devices （1） ~ec （ac） ; Assume that the discharge characteristic function of the i-th energy storage device is function Qe （i） , the discharge power of the i-th energy storage device is pe （i） , the discharge time of the i-th energy storage device is te （i） , define pe （i） and te （i） The analytical expression of : Let Qe （i） The mathematical expression is: ; β （0） （i） ~β （5） （i） Represents the function Qe （i） The coefficient of δ （i） Represents the function Qe （i） The correction value of i is in the range of 1 to ac. Introducing the Lagrange multiplier λ, the equation L is constructed: ; Calculation equation L versus pe （i） and te （i） The partial derivative of , and set the partial derivative to zero to eliminate λ, we get formula A-1: ; Simplifying formula A-1, we get formula A-2: ; Will Substituting into formula A-2, we get formula A-3: ; will pe （i）As p,te （i） As t, ec （i） As q, As a, As b, As c, As d, calculate pe （i） With te （i） The analytical expression of .

[0009] Furthermore, calculate pe （i） With te （i） The specific steps of the analytical expression are as follows: Arranging formula A-3, we get formula B-1: ; Introduce auxiliary parameter o, let , we get formula B-2: ; Among them, Aa, Bb and Cc represent the new parameters of p: ; Introducing parameter k a second time, rewriting formula B-2, we get formula B-3: ; Define formula B-4: ; Calculate the smallest real root km of Formula B-4; Based on km, we get formula B-5: ; Based on formula B-5, we get the solution for o: ; According to the solution of o, we get pe （i） The analytical expression of : ; According to pe （i） The analytical expression of te （i） Analytical formula: ; Function according to pe （i） and te （i） , calculate the power and discharge time of the 1st to acth energy storage devices.

[0010] Furthermore, the step of adjusting the discharge time of each energy storage device and the discharge time according to the discharge characteristic function further includes: like < , then obtain the reactive power Qu of the power supply end in the target area and adjust the transmission power and reactive power of the power supply end; Calculate the additional power ΔPu of the power transmission at the power supply end: ; Calculate the power factor σ of the target area, σ=Pu / Qu; Assuming σ remains unchanged, the new reactive power when the transmission power increases by ΔPu is Quu: ; Assuming that Qu remains unchanged, the new power factor σu of the transmission power increased by ΔPu is: ; Calculate the coefficient of variation of reactive power bQ: ; The coefficient of variation of power factor bσ: ; Compare the size of bQ and bσ, and adjust the reactive power and voltage phase angle of the power supply end; If bQ ≥ bσ, then increase the voltage phase angle at the power supply end by (1 + bσ) times; If bQ<bσ, the reactive power at the power supply end is adjusted to Quu.

[0011] Furthermore, the specific steps of constructing the admittance matrix are as follows: Count the number of power consumption areas ual in the first municipal district, obtain the topological relationship between the first municipal district and the power consumption areas, and construct the initial admittance matrix of the first municipal district; obtain the expected power consumption rate and reactive power of the ualth to ualth power consumption areas in the first municipal district, and obtain the complex power se of the first power consumption area in the first municipal district （1） ~se （ual） ; Get the power supply voltage ug from the 1st to the ualth power consumption area （1） ~ug （mn） ; Calculate the equivalent admittance Y of the first power consumption area （1,1） ~Y （1,ual） ; Determine whether the first power consumption area is directly connected to the second power consumption area; If they are connected, the equivalent admittance Y of the first power consumption area with respect to the second power consumption area is （1,2） : ; If they are not connected, the equivalent admittance Y of the first power consumption area with respect to the second power consumption area is （1,2） =0; Calculate the equivalent admittance Y of the first power consumption area with respect to the third to the ual power consumption areas （1,3） ~Y （1,ual） ; Calculate Y （1,3）~Y （1,ual） The sum of the self-admittance Y of the first power consumption area （1,1） .

[0012] Furthermore, the step of constructing the admittance matrix also includes: Assume that the complex power of the mth electricity consumption area is se （m） , the supply voltage is ug （m） ; The complex power of the nth electricity consumption area is se （n） , the supply voltage is ug （n） ; Assume that the equivalent admittance of the mth power consumption area with respect to the nth power consumption area is Y （m,n） , define calculation formula C-1-1: ; Assume that the equivalent admittance of the mth power consumption area with respect to the lth power consumption area is Y （m,l） ; Let the self-admittance of the mth power consumption area be Y （m,m） , define calculation formula C-1-2: ; Calculate the equivalent admittance and self-admittance Y from the second power consumption area to the ual power consumption area （2,l） ~Y （ual,ual） ; Construct the initial admittance matrix YY （1） : ; Extract Y （1,l） ~Y （ual,ual） The real part of the first power consumption area to the ual power consumption area is obtained. （1,l） ~G （ual,ual） ; Extract the imaginary part and get the equivalent susceptance B of the first to the ual power consumption area （1,l） ~B （ual,ual） ; Construct the initial admittance matrix of the 2nd to mnth districts; according to the matrix YY （1） , construct the calculation equation for the initial power flow calculation of the first urban district, and calculate the initial phase of the power supply voltage of each power consumption area in the first urban district.

[0013] Furthermore, the specific steps for calculating the initial phase of the supply voltage are as follows: Assume that the supply voltage of the wth power consumption area is Uz （w） , the voltage phase is θ （w） , the expected power consumption is Pz （w） , reactive power is Qz （w） ; The complex power of the w-th power consumption area is Sz （w） , Sz （w） =Pz （w） +Qz（w） ×j; The supply voltage of the vth power consumption area is Ug （v） , the voltage phase is θ （v） ; The equivalent admittance of the wth power consumption area with respect to the vth power consumption area is Y （w，v） ; Y （w，v） The corresponding equivalent conductance is G （w,v） , the equivalent reactance is B （w,v） ; Assume that the conjugate of the supply voltage of the wth power consumption area is (Uz （w） ) * , define formula C-2-1: ; The conjugate of the supply voltage of the vth power consumption area is (Uz （v） ) * , ; Define formula C-2-2: ; Expanding formula C-2-2, we get formula C-2-3: ; Among them, θ （w，v） Indicates the phase difference between the supply voltage of the wth power consumption area and the vth power consumption area; Construct Matrix I （1） and matrix U （1） ;Define formula C-3: ; Assume that the initial phase of the supply voltage from the 1st to the ualth power consumption area is θz （1） ~θz （mn） ; Arrange formula C-3 to get formula D （1） ~Formula D （ual） : Formula D （1） for: ; Similarly, Formula D （ual） for: ; Using the Newton-Raphson algorithm, calculate θz （1） ,θz （2） ~θz （mn） The value of is used to obtain the initial phase of the supply voltage of the first to the ual power consumption area; Calculate the initial phase of the corresponding power consumption areas of the 2nd to mnth municipal districts; obtain the actual power consumption of all power consumption areas in the first district, and adjust the power supply power of the 1st to ualth power consumption areas.

[0014] Furthermore, the specific steps for adjusting the power supply of the first to the ualth power consumption areas are as follows: Get the expected power consumption rate po of the first to the ual power consumption area （1） ~po （mn） , actual power consumption rate pa （1） ~pa （mn） ; Determine whether the expected power consumption rate is equal to the actual power consumption rate, and adjust the power supply voltage and phase from each municipal district to the corresponding power consumption district based on the relationship between the actual power consumption and the expected power consumption; In the first to the ualth power consumption zones, if there is only one power consumption zone whose expected power consumption rate is not equal to the actual power consumption rate, the power consumption zone with the different expected power consumption rate and the actual power consumption rate is regarded as an abnormal zone; Obtain the expected power consumption rate pox, actual power consumption rate pax, actual supply voltage ux, and actual supply voltage phase θx of the abnormal area; Obtain the expected power supply voltage uy and the initial power supply voltage phase θy in the abnormal area; Determine whether ux and uy, as well as θx and θy, are equal, and adjust the power supply voltage and power supply voltage phase in the abnormal area or other areas; If ux and uy are equal, it means that the reactive power in the abnormal area has changed; Compare the sizes of pox and pax, and adjust the power supply voltage of other areas; If pox>pax, then reduce the power supply voltage of other areas times; If pox < pax, then increase the supply voltage of other areas times; If θx and θy are equal, it means that the active power in the abnormal area has changed; Compare the size of pox and pax, and adjust the phase of the power supply voltage in other areas; If pox>pax, then increase the phase of the power supply voltage in other areas times; If pox<pax, then reduce the phase of the power supply voltage in other areas times; If ux and uy as well as θx and θy are not equal, the abnormal area is used as the balance node and the matrix YY is reconstructed. （1） The actual power consumption and actual reactive power of all power consumption areas shall be adjusted.

[0015] Furthermore, the specific steps for adjusting the actual power consumption and actual reactive power of all power consumption areas are as follows: Assume that the abnormal area is numbered as r; recalculate the equivalent admittance rY of the abnormal area with respect to the 1st to ualth power consumption areas （r，1） ~rY （r，ual） , self-admittance rY （r，r） ; Calculate the change in the self-admittance of the abnormal area ΔY （r，r） ; Change in equivalent admittance ΔY （r，1） ; Calculate the changes in the equivalent admittance and self-admittance of the first to the ual power station, and obtain ΔY （1，1） ~ΔY （ual，ual） , get the matrix YY （2） ; Calculate the voltage change Δu at the balance node, ; Phase change Δθ, ; Calculate the change in actual power supply voltage from the 1st to the ualth power station Δu （1） ~Δu （ual） , the actual supply voltage phase change Δθ （1） ~Δθ （ual） ; Assume that the change in actual power consumption from the 1st to the ualth power station is ΔP （1） ~ΔP （ual） , the actual change in reactive power is ΔQ （1） ~ΔQ （ual） ; Construct matrix Zz （1） : Matrix Zz （2） : ; Define formula C-4-2: ; Based on Formula C-4-2, use the Newton-Raphson algorithm to calculate ΔP （1） ~ΔP （ual） and ΔQ （1） ~ΔQ （ual） The value of Increase the actual power consumption of the first to the ual power stations by ΔP （1） ~ΔP （ual） , the actual reactive power increases in turn by ΔQ （1） ~ΔQ （ual） ; In the first to the ualth power consumption areas, if the expected power consumption rates of multiple power consumption areas are not equal to the actual power consumption rates, the step of adjusting the actual power consumption and the actual reactive power of all power consumption areas is repeated.

[0016] Compared with the prior art, the present invention has the following beneficial effects: Multi-dimensional data fusion: This invention integrates heterogeneous data such as historical power consumption in the city, discharge characteristics of energy storage devices, and topological relationships of power consumption areas into a unified framework, breaking through the limitation of traditional power regulation relying on a single dimension. By comparing expected power consumption with actual power supply capacity in real time, when the power supply is insufficient, the system can automatically adjust the discharge strategy of the energy storage device or the parameters of the power supply end to ensure the stability of the power grid, thus realizing two-way regulation between the transmission end and the power consumption end.

[0017] Hierarchical control architecture: A three-level control structure of "city-municipal district-power consumption zone" is adopted, taking into account both global optimization and local refined control. The city level coordinates energy storage devices and power supply terminals, the municipal district level uses the ARIMA model to predict demand, and the power supply parameters are adjusted using the admittance matrix at the power consumption zone level, forming a hierarchical decision-making system. At the same time, the present invention realizes quantitative control of the discharge process by constructing the discharge power-time-heating relationship (Qe(i) function) of the energy storage device. The Lagrange multiplier method is used to optimize the discharge power (pe(i)) and time (te(i)), minimizing heat loss while meeting the discharge capacity (ec(i)) requirements.

[0018] Flexibility and robustness of power consumption-side regulation: By real-time monitoring of the deviation between the actual power consumption of the power consumption area and the expected value, the present invention can construct and adjust the admittance matrix parameters to optimize the power supply voltage and phase. When a power gap occurs in a power consumption area, cross-regional power support is achieved by reconstructing the matrix balancing node. The "balancing node" strategy is adopted to isolate and reconstruct abnormal power consumption areas to prevent the spread of local faults. Optimizing the power supply voltage and phase can reduce line losses and improve the utilization rate of power resources. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] Other features, objects and advantages of the present invention will become more apparent upon reading the detailed description of non-limiting embodiments with reference to the following drawings: Figure 1 Schematic diagram of the method of the present invention; Figure 2 This is a schematic diagram of the topology of the power consumption area of ​​the present invention; Figure 3 This is a schematic diagram of the power consumption area topology of the present invention. DETAILED DESCRIPTION

[0020] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0021] See also Figure 1 , a method for establishing an incremental model of controllable power resources includes: Step S1: Obtain the number of districts in the target area, obtain the historical power consumption of each district, perform time series analysis on the historical power consumption, and calculate the expected power consumption of each district; Obtain the hourly electricity consumption of each municipal district over the past three days as historical electricity consumption; perform time series analysis on the historical electricity consumption and calculate the differential order, autoregressive parameter, and moving average parameter corresponding to each municipal district's historical electricity consumption; calculate the autocorrelation coefficient and moving average coefficient corresponding to each municipal district based on the differential order, autoregressive parameter, and moving average parameter corresponding to each municipal district; then, use the ARIMA model, combined with the autocorrelation coefficient and moving average coefficient corresponding to each municipal district, to estimate the total electricity consumption of each municipal district in the next day as the expected electricity consumption of each municipal district; It should be noted that the "target area" in the present invention refers to: a city-level area where the present invention (a method for establishing an incremental model of controllable electric resources) is used to perform incremental modeling of electric resources and control of electric resources; Step S2: (Power Transmission Control) Obtain the number of energy storage devices in the target area and construct a discharge characteristic function for each energy storage device. Obtain the transmission power at the target area's power supply terminal and the remaining power of each energy storage device. Combined with the expected power consumption of each municipal district, determine whether the transmission power is sufficient. If sufficient, adjust the discharge time and discharge duration of each energy storage device based on the discharge characteristic function. If insufficient, adjust the transmission power. It should be noted that the "energy storage device" in the present invention refers to an electrical energy storage device; the "power supply end" refers to a power station or national power grid that supplies power to the target area; The specific steps of step S2 are as follows: Count the number of energy storage devices (in the target area) ac, construct the discharge characteristic function of the 1st to acth energy storage devices, and obtain the function Qe （1） ~Function Qe （ac） ; Construct the discharge characteristic function Qe of the first energy storage device （1） ; Get the average discharge power pd of the first energy storage device on the first, second, and mo days of the past month （1） 、pd （2） ~pd （mo） ; Discharge time td （1） 、td （2） ~td （mo） ; The heat generated by the first energy storage device due to discharge qd （1） 、qd （2） ~qd （mo） ; Among them, mo represents the number of days in the past month; Construct a (mo×6) matrix X: ; Construct a (mo×1) matrix Y: ; Let function Qe （1） The coefficient is β （0） , β （1） ~β （5） ; Construct a (6×1) matrix Z: ; Calculating β （0） ~β （5） Value: ; Where T represents the transpose of the matrix and -1 represents the inverse of the matrix; Assume the first discharge power is pe, discharge time is te, define the function Qe （1） The mathematical expression is: ; Where δ represents the function Qe （1） The correction value of (the initial value of δ is 0); PD （1） ~pd （mo） As pe,td （1） ~td （mo） As te,qd （1） ~qd （mo） As a function of Qe （1） Substitute the value of into the function Qe （1） In the mathematical expression of , calculate the value of δ and get the function Qe （1） ; Repeatedly construct function Qe （1） The same steps are used to construct the discharge characteristic functions of the 2nd to acth energy storage devices, and the function Qe is obtained. （2） ~Function Qe （ac） ; Get the expected electricity consumption of the 1st to mnth municipal districts and get el （1） ~el （mn） ; Where mn represents the number of urban districts; calculate el （1） ~el （mn） of and ael; Get the remaining power of the first, second, and ac energy storage devices, and get re （1） 、re （2） ~re （ac） ; Calculate re （1） ~re （ac） of and are; Obtain the transmission power Pu of the power supply end (in the target area); judge whether it is established; If established, no processing will be done; If not, compare and The size of the energy storage device is used to adjust the discharge power and discharge time or the transmission power and reactive power at the power supply end; like ≥ , then calculate the electric energy replenishment amount ree of the energy storage device: ; Adjust the discharge power and discharge time of the energy storage device according to ree and the discharge characteristic function of each energy storage device; Calculate the discharge capacity ec of the first, second, and ac energy storage devices （1） 、ec （2） ~ec （ac） ; Among them, ec （1） The calculation formula is: ; ec （2） The calculation formula is: ; By analogy, ec （ac） The calculation formula is: ; Assume that the discharge characteristic function of the i-th energy storage device is function Qe （i） , the discharge power of the i-th energy storage device is pe （i） , the discharge time of the i-th energy storage device is te （i） , define pe （i） and te （i） The analytical expression of : Let Qe （i） The mathematical expression is: ; Among them, β （0） （i） ~β （5） （i） Represents the function Qe （i） The coefficient of δ （i） Represents the function Qe （i） The correction value of i is in the range of 1 to ac. Define formula A-1: ;ec （i） represents the discharge capacity of the i-th energy storage device; Introducing the Lagrange multiplier λ, the equation L is constructed: ; Calculation equation L versus pe （i） and te （i）The partial derivative of , and set the partial derivative to zero, we get formula A-2: ; Eliminating λ from Formula A-2, we obtain Formula A-3: ; Simplifying Formula A-3, we get Formula A-4: ; Based on formula A-1, we know that if and only if When , the discharge effect of the i-th energy storage device is the best (i.e., the discharge heat generation is the lowest); Will Substituting into formula A-4, we get formula A-5: ; will pe （i） As p,te （i） As t, ec （i） As q, As a, As b, As c, as d; Calculate PE based on Formula A-5 （i） With te （i） The analytical expression of : Arranging formula A-5, we get formula B-1: ; Introduce auxiliary parameter o, let , eliminating the cubic term of p in formula B-1, we get formula B-2: ; Among them, Aa, Bb and Cc represent the new parameters of p: ; Introducing the parameter k twice, rewriting Formula B-2 into a perfect square form, we get Formula B-3: ; To make the right side of Formula B-3 also satisfy the perfect square form, the discriminant of the right side of Formula B-3 must be zero, and we get Formula B-4-1: ; Expanding formula B-4-1, we get the cubic equation for k, which is formula B-4-2: ; Assume that the smallest real root of formula B-4-2 is kz, and define formula S to test the value of kz; the expression of formula S is: ; Use Cardano's formula (the formula for finding roots of cubic equations) to calculate formula B-4-2 and extract the smallest real root km that satisfies the formula S test; Based on km, we can rearrange Formula B-3 to obtain Formula B-5: ; Based on formula B-5, we get the solution for o: ; According to the solution of o, we get pe （i） The analytical expression of : ; According to pe （i） The analytical expression of te （i） Analytical formula: ; Function according to pe （i） and te （i） Calculate the power and discharge time of the 1st to acth energy storage devices using the analytical formula; like < , then obtain the reactive power Qu of the power supply end in the target area and adjust the transmission power and reactive power of the power supply end; Calculate the additional power ΔPu of the power transmission at the power supply end: ; Calculate the power factor σ of the target area, σ=Pu / Qu; Assuming σ remains unchanged, the new reactive power when the transmission power increases by ΔPu is Quu: ; Assuming that Qu remains unchanged, the new power factor σu of the transmission power increased by ΔPu is: ; Calculate the coefficient of variation of reactive power bQ: ; The coefficient of variation of power factor bσ: ; Compare the size of bQ and bσ, and adjust the reactive power and voltage phase angle of the power supply end; If bQ ≥ bσ, then increase the voltage phase angle at the power supply end by (1 + bσ) times; If bQ<bσ, the reactive power at the power supply end is adjusted to Quu.

[0022] Step S3: (Power Consumption Control) Perform a time series analysis of the historical power consumption of each municipal district to the power application area to obtain the expected power consumption of each municipal district to the power application area and construct an admittance matrix. Obtain the actual power consumption of each municipal district to the power application area in real time and compare it with the expected power consumption. Based on the relationship between the actual power consumption and the expected power consumption, modify the admittance matrix and adjust the power supply voltage, phase, and power of each municipal district to the power application area. It should be noted that the “electricity consumption area” in the present invention means: the area that consumes electricity within each municipal area, such as residential areas, commercial areas, industrial areas, public facilities areas, etc.; The specific steps of step S3 are as follows: Step S31: Count the number of electricity consumption areas in the first, second, and mnth municipal districts, recorded as ua （1） 、ua （2） ~ua （mn） ; Obtain the hourly power consumption of each municipal district for the past three days, perform a time series analysis on the historical power consumption of each municipal district for the application area, and obtain the expected power consumption of each municipal district for the application area; Obtain the topological relationship between the first urban district and the first, second, and finally the ualth electricity consumption area, and construct the initial admittance matrix of the first urban district; ua （1） As ual, obtain the expected power consumption rate po of the first, second, and finally ual power consumption areas in the first municipal district. （1） 、po （2） ~po （ual） ; Get the (average) reactive power qo of the first, second, and third electricity consumption areas in the first municipal district yesterday （1） 、po （2） ~po （ual） ; According to the expected power consumption rate and reactive power, the complex power se of the first power consumption area in the first urban district is combined （1） ,se （1） =po （1） +qo （1） ×j; Complex power se of the second electricity consumption area （2） ,se （2） =po （2） +qo （2） ×j; By analogy, the complex power se of the ual power consumption area is （ual） ,se （ual） =po （ual） +qo （ual） ×j; where j represents an imaginary unit; Get the power supply voltage ug from the first municipal district (yesterday) to the first, second, and finally the ual power consumption area （1） ,ug （2） ~ug （mn） ; Calculate the equivalent admittance Y of the first power consumption area （1,1） ~Y （1,ual） ; Determine whether the first power consumption area is directly connected to the second power consumption area; See also Figure 2 , if connected, then the equivalent admittance Y of the first power consumption area with respect to the second power consumption area is （1,2） : ; See also Figure 3 , if they are not connected, then the equivalent admittance Y of the first power consumption area with respect to the second power consumption area is （1,2） =0; (Repeat Y （1,2） Calculate the equivalent admittance Y of the first power consumption area with respect to the third to the ual power consumption areas （1,3） ~Y （1,ual） ; Calculate Y （1,3） ~Y （1,ual） The sum of the self-admittance Y of the first power consumption area （1,1） ; Assume that the complex power of the mth electricity consumption area is se （m） , the power supply voltage (from the first municipal district to the mth power consumption area) is ug （m） ; The complex power of the nth electricity consumption area is se （n） , the power supply voltage (the first municipal district supplies power to the nth power consumption area) is ug （n） ; Assume that the equivalent admittance of the mth power consumption area with respect to the nth power consumption area is Y （m,n） , define calculation formula C-1-1: ; The value range of m is: 2~(ual-1); The value range of n is: (m+1)~ual; Assume that the equivalent admittance of the mth power consumption area with respect to the lth power consumption area is Y （m,l） , the value range of l is: 1~(m-1); let the self-admittance of the mth power consumption area be Y （m,m） , define calculation formula C-1-2: ; According to the calculation formula C-1-1 and C-1-2, calculate the equivalent admittance Y of the second power consumption area with respect to the first to the ual （2,l） ~Y（2,ual） and self-admittance Y （2,2） ; The equivalent admittance Y of the third power consumption area with respect to the first to the second （3,l） ~Y （3,ual） and self-admittance Y （3,3） ; Similarly, the equivalent admittance Y of the ualth power consumption area with respect to the 1st to (ual-1)th （ual,l） ~Y （ual,（ual-1）） and self-admittance Y （ual,ual） ; Construct a (ual×ual) zero matrix and fill it into Y （1,l） ~Y （ual,ual） Zero matrix, get the initial admittance matrix YY of the first district （1） : ; Among them, the matrix YY （1） The diagonal elements are the self-admittances of the first to the ual power consumption areas, and the non-diagonal elements are equal to the "elements symmetrical about the diagonal" (similar to (-Y （1,2） )=(-Y （2,1） ), (-Y （ual,1） )=(-Y （1,ual） )); Extract Y （1,l） 、Y （1,2） ~Y （1,ual） The real part of the first power consumption area is obtained by （1,l） , G （1,2） ~G （1,ual） ; Extract Y （1,l） 、Y （1,2） ~Y （2,ual） The imaginary part of the first power consumption area is obtained by （1,l） 、B （1,2） ~B （1,ual） ; Extract Y （1,l） 、Y （1,2） ~Y （1,ual） The real part of the first power consumption area is obtained by （1,l） , G （1,2） ~G （1,ual） ; Extract Y （1,l） 、Y （1,2） ~Y （2,ual） The imaginary part of the first power consumption area is obtained by （1,l） 、B （1,2） ~B （1,ual） ; Extract Y （2,l） 、Y （2,2） ~Y （1,ual）The real part of the second power consumption area is obtained by （2,l） , G （2,2） ~G （2,ual） ; Extract Y （2,l） 、Y （2,2） ~Y （2,ual） The imaginary part of the first power consumption area is obtained by （2,l） 、B （2,2） ~B （2,ual） ; Similarly, extract Y （ual,l） 、Y （ual,2） ~Y （ual,ual） The real part of the first power consumption area is obtained by （ual,l） , G （ual,2） ~G （ual,ual） ; Extract Y （ual,l） 、Y （ual,2） ~Y （ual,ual） The imaginary part of the first power consumption area is obtained by （ual,l） 、B （ual,2） ~B （ual,ual） ; Repeating matrix YY （1） The construction steps are to construct the admittance matrix of the 2nd to mnth urban districts; According to the matrix YY （1） , construct the calculation equation for the initial power flow calculation of the first urban district, and calculate the initial phase of the power supply voltage of each power consumption area in the first urban district; Assume that the supply voltage of the wth power consumption area (within the first city district) is Uz （w） , the voltage phase is θ （w） , the expected power consumption is Pz （w） , reactive power is Qz （w） ; The complex power of the w-th power consumption area is Sz （w） , Sz （w） =Pz （w） +Qz （w） ×j; The supply voltage of the vth power consumption area is Ug （v） , the voltage phase is θ （v） ; The equivalent admittance of the wth power consumption area with respect to the vth power consumption area is Y （w，v） ; Y （w，v） The corresponding equivalent conductance is G （w,v） , the equivalent reactance is B （w,v） ; Among them, the value range of w and v is: 1~ual; Assume that the conjugate of the supply voltage of the wth power consumption area is (Uz （w） ) * , define formula C-2-1: ; The conjugate of the supply voltage of the vth power consumption area is (Uz （v） ) * , ; Define formula C-2-2: ; Expanding formula C-2-2, we get formula C-2-3: ; Among them, θ （w，v） Indicates the phase difference between the supply voltage of the wth power consumption area and the vth power consumption area; Construct a (ual×1) matrix I （1） : ; Construct a (ual×1) matrix U （1） : ; Define formula C-3-1: ; Expand and organize formula C-3-1 to obtain formula C-3-2: ; Assume that the initial phase of the supply voltage of the first, second, and finally the ual power consumption area is θz （1） ,θz （2） ~θz （mn） ; will po （1） ~po （mn） and qo （1） ~po （mn） Substitute into formula C-3-2 and sort it out to get formula D （1） Formula D （2） ~Formula D （ual） : Formula D （1） for: ; Formula D （2） for: ; By analogy, formula D （ual） for: ; Using the Newton-Raphson algorithm, calculate θz （1） ,θz （2） ~θz （mn） The value of is used to obtain the initial phase of the power supply voltage of the first, second, and finally the ual power consumption area (within the first municipal district); Repeat θz （1） ~θz（mn） The calculation process is to calculate the initial phase of the power application area of ​​the 2nd to mnth municipal districts; Obtain the actual power consumption of all power consumption areas within the first jurisdiction, and adjust the power supply to the first to the ual power consumption areas based on the expected power consumption of each power consumption area; Get the expected power consumption rate po of the first, second, and third power consumption areas in the first city district （1） 、po （2） ~po （mn） ; Get the actual power consumption rate pa of the first, second, and third power consumption areas in the first city district （1） 、pa （2） ~pa （mn） ; Determine whether the expected power consumption rate and the actual power consumption rate of the first to the ual power consumption areas are equal, and adjust the power supply voltage and phase of each municipal district to the corresponding power consumption area according to the relationship between the actual power consumption and the expected power consumption; In the first to the ualth power consumption zones, if there is only one power consumption zone whose expected power consumption rate is not equal to the actual power consumption rate, the power consumption zone with the different expected power consumption rate and the actual power consumption rate is regarded as an abnormal zone; Get the expected power consumption rate pox in the abnormal area, pox∈{po （1） 、po （2） ~po （ual）}; actual power consumption rate pax, pax∈{pa （1） 、pa （2） ~pa （mn）}; Get the actual power supply voltage ux in the abnormal area, ux∈{ug （1） ,ug （2） ~ug （ual）}; Actual supply voltage phase θx, θx∈{θz （1） ,θz （2） ~θz （ual）}; Obtain the expected power supply voltage uy and the initial power supply voltage phase θy in the abnormal area; Determine whether ux and uy, as well as θx and θy, are equal, and adjust the power supply voltage and power supply voltage phase of the abnormal area or other areas (other areas refer to the other power consumption areas within the first municipal district, from the first to the ualth power consumption areas, excluding the abnormal area); If ux and uy are equal, it means that the reactive power in the abnormal area has changed; Compare the sizes of pox and pax, and adjust the power supply voltage of other areas; If pox>pax, then reduce the power supply voltage of other areas times; If pox < pax, then increase the supply voltage of other areas times; If θx and θy are equal, it means that the active power in the abnormal area has changed; Compare the size of pox and pax, and adjust the phase of the power supply voltage in other areas; If pox>pax, then increase the phase of the power supply voltage in other areas times; If pox<pax, then reduce the phase of the power supply voltage in other areas times; If ux and uy as well as θx and θy are not equal, the abnormal area is used as the balance node and the matrix YY is reconstructed. （1） If yes, adjust the actual power consumption and actual reactive power of all power consumption areas; Let the number of the abnormal area be r, and the value range of r is: 1~ual; According to the calculation formula C-1-1 and Y （1,1） ~Y （1,ual） The calculation formula is used to recalculate the equivalent admittance rY of the abnormal area with respect to the first, second, and finally the ual power consumption area. （r，1） 、rY （r，2） ~rY （r，ual） , self-admittance rY （r，r） ; Get matrix YY （1） The equivalent admittance tY of the abnormal area in the first, second, and third power consumption areas （r，1） ,tY （r，2） ~tY （r，ual） , self-admittance rY （r，r） ; Calculate the change in the self-admittance of the abnormal area ΔY （r，r） , ; Calculate the change in the equivalent admittance of the abnormal area and obtain ΔY （r，1） , ΔY （r，2） ~ΔY （r，ual） ; in, , And so on, ; The equivalent admittance of the wth power consumption area with respect to the abnormal area is rY （w，r） ; The self-admittance of the w-th power consumption area is rY （w，w） ; Assume that the change in self-admittance of the w-th power consumption area is ΔY （w，w） , equivalent admittance Y （w，v） The change is ΔY （w，v） ; Define calculation formula C-4-1: ; According to the calculation formula C-4-1, calculate the changes in the equivalent admittance and self-admittance of the first, second, and finally the ual power station, and obtain ΔY （1，1） , ΔY （1，2） ~ΔY （ual，ual） ; The matrix YY （1） The elements in the equation are replaced by ΔY （1，1） , ΔY （1，2） ~ΔY （ual，ual） , get the matrix YY （2） ; The abnormal area is regarded as the balance node, and the voltage change Δu of the balance node is calculated. ; Phase change Δθ, ; Get the actual power supply voltage uv of the first, second, and finally the ual power station (excluding abnormal areas) （1） 、uv （2） ~uv （ual） and the actual supply voltage phase θv （1） 、θv （2） ~θv （ual） ; Calculate the change in actual power supply voltage Δu of the first, second, and finally the ual power station (excluding abnormal areas) （1） , Δu （2） ~Δu （ual） , the actual supply voltage phase change Δθ （1） , Δθ （2） ~Δθ （ual） ; in, , , and so on, ; Assume that the change in actual power consumption of the first, second, and finally the ual power station (excluding abnormal areas) is ΔP （1） , ΔP （2） ~ΔP （ual） , the actual change in reactive power is ΔQ （1） , ΔQ （2） ~ΔQ （ual） ; Construct two (ual×2) matrices and get the matrix Zz （1） and matrix Zz （2） ; Matrix Zz （1） : Matrix Zz （2） : ; Define formula C-4-2: ; Repeat formula D based on formula C-4-2 （1） ~Formula D （ual） The expansion method and the Newton-Raphson algorithm are used to calculate ΔP （1） , ΔP （2） ~ΔP （ual） and ΔQ （1） , ΔQ （2） ~ΔQ （ual） The value of Increase the actual power consumption of the first, second, and third power stations (excluding abnormal areas) by ΔP （1） , ΔP （2） ~ΔP （ual） , the actual reactive power increases in turn by ΔQ （1） , ΔQ （2） ~ΔQ （ual） ; In the first to the ualth power consumption areas, if the expected power consumption rate and actual power consumption rate of multiple power consumption areas are not equal, then the power consumption area with the fewest directly connected power consumption areas and the smallest difference between the actual power consumption and the expected power consumption is selected as the balancing node, and the steps of adjusting the actual power consumption and actual reactive power of all power consumption areas are repeated to adjust the actual power consumption and actual reactive power of all power consumption areas; Repeat the steps of adjusting the power supply voltage, power supply voltage phase, actual power consumption and reactive power of all power consumption areas within the first municipal district, and adjust the power consumption areas corresponding to the second to mnth municipal districts.

[0023] The above formulas are all dimensionless and calculated by taking their numerical values. The formula is a formula for the most recent real situation obtained by collecting a large amount of data and performing software simulation. The preset parameters in the formula are set by technicians in this field according to actual conditions. For example, if there are weight coefficients and proportional coefficients, the size of the settings is to quantify each parameter to obtain a specific value, which is convenient for subsequent comparison. Regarding the size of the weight coefficient and the proportional coefficient, as long as it does not affect the proportional relationship between the parameter and the quantized value, it is fine.

[0024] Finally, it should be noted that the above-described embodiments are only specific implementations of the present invention, which are used to illustrate the technical solutions of the present invention, rather than to limit them. The scope of protection of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the above-described embodiments, those skilled in the art should understand that any person skilled in the art can modify or easily conceive of changes to the technical solutions described in the above-described embodiments within the technical scope disclosed by the present invention, or replace some of the technical features therein with equivalents. Such modifications, changes, or replacements do not deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be based on the scope of protection of the claims.

Claims

1. A method for establishing an incremental model of controllable power resources, characterized in that: The method comprises: Obtain the number of districts in the target area, obtain the historical power consumption of each district, perform time series analysis on the historical power consumption, and calculate the expected power consumption of each district; Obtain the number of energy storage devices in the target area and construct a discharge characteristic function for each energy storage device; obtain the transmission power of the target area's power supply terminal and the remaining power of each energy storage device, and determine whether the transmission power is sufficient based on the expected power consumption of each municipal district; if sufficient, adjust the discharge time and discharge duration of each energy storage device according to the discharge characteristic function; if not sufficient, adjust the transmission power; A time series analysis is performed on the historical power consumption of each municipal district to the application power area to obtain the expected power consumption of each municipal district to the application power area, and an admittance matrix is ​​constructed; the actual power consumption of each municipal district to the application power area is obtained in real time and compared with the expected power consumption. According to the relationship between the actual power consumption and the expected power consumption, the admittance matrix is ​​corrected, and the power supply voltage, phase and power of each municipal district to the application power area are adjusted.

2. The method for establishing an incremental model of controllable power resources according to claim 1, characterized in that: The specific steps for constructing the discharge characteristic function of each energy storage device are as follows: Count the number of energy storage devices ac and construct the discharge characteristic function Qe of the first energy storage device （1） ; Get the average discharge power pd of the first energy storage device from the first to the moth day in the past month （1） ~pd （mo） ;Discharge time td （1） ~td （mo） ; Calories qd （1） ~qd （mo） ;Construct matrix X and matrix Y; Let function Qe （1） The coefficient is β （0） ~β （5） , construct matrix Z; Calculating β （0） ~β （5） Value: ; Assume the first discharge power is pe, discharge time is te, define the function Qe （1） The mathematical expression is: ； δ represents the correction value; PD （1） ~pd （mo） As pe,td （1） ~td （mo） As te,qd （1） ~qd （mo） Substitute into the function Qe （1） Calculate δ from the value of Construct the discharge characteristic function of the 2nd to acth energy storage devices and obtain the function Qe （2） ~Function Qe （ac） .

3. The method for establishing an incremental model of controllable power resources according to claim 1, characterized in that: The specific steps for adjusting the discharge time and discharge time of each energy storage device according to the discharge characteristic function are as follows: Obtain and calculate the expected electricity consumption and ael of all city districts; Obtain and calculate the sum of the remaining power of all energy storage devices; Obtain the transmission power Pu of the power supply end; judge whether it is established; If established, no processing will be done; If not, compare and The size of the energy storage device is used to adjust the discharge power and discharge time or the transmission power and reactive power at the power supply end; like ≥ , then calculate the electric energy replenishment amount ree of the energy storage device; according to ree and the discharge characteristic function of each energy storage device, adjust the discharge power and discharge time of the energy storage device; Calculate the discharge capacity ec of the 1st to acth energy storage devices （1） ~ec （ac） ; Assume that the discharge characteristic function of the i-th energy storage device is function Qe （i） , the discharge power of the i-th energy storage device is pe （i） , the discharge time of the i-th energy storage device is te （i） , define pe （i） and te （i） The analytical expression of : Let Qe （i） The mathematical expression is: ； β （0） （i） ~β （5） （i） Represents the function Qe （i） The coefficient of δ （i） Represents the function Qe （i） The correction value of i is in the range of 1 to ac. Introducing the Lagrange multiplier λ, the equation L is constructed: ； Calculation equation L versus pe （i） and te （i） The partial derivative of , and set the partial derivative to zero to eliminate λ, we get formula A-1: ； Simplifying formula A-1, we get formula A-2: ； Will Substituting into formula A-2, we get formula A-3: ； will pe （i） As p,te （i） As t, ec （i） As q, As a, As b, As c, As d, calculate pe （i） With te （i） The analytical expression of .

4. The method for establishing an incremental model of controllable power resources according to claim 3, characterized in that: Calculate PE （i） With te （i） The specific steps of the analytical expression are as follows: Arranging formula A-3, we get formula B-1: ； Introduce auxiliary parameter o, let , we get formula B-2: ； Among them, Aa, Bb and Cc represent the new parameters of p: ； Introducing parameter k a second time, rewriting formula B-2, we get formula B-3: ； Define formula B-4: ； Calculate the smallest real root km of Formula B-4; Based on km, we get formula B-5: ； Based on formula B-5, we get the solution for o: ； According to the solution of o, we get pe （i） The analytical expression of : ； According to pe （i） The analytical expression of te （i） Analytical formula: ； Function according to pe （i） and te （i） , calculate the power and discharge time of the 1st to acth energy storage devices.

5. The method for establishing an incremental model of controllable power resources according to claim 3, characterized in that: The step of adjusting the discharge time and discharge time of each energy storage device according to the discharge characteristic function also includes: like < , then obtain the reactive power Qu of the power supply end in the target area and adjust the transmission power and reactive power of the power supply end; Calculate the additional power ΔPu of the power transmission at the power supply end: ； Calculate the power factor σ of the target area, σ=Pu / Qu; Assuming σ remains unchanged, the new reactive power when the transmission power increases by ΔPu is Quu: ; Assuming that Qu remains unchanged, the new power factor σu of the transmission power increased by ΔPu is: ; Calculate the coefficient of variation of reactive power bQ: ; The coefficient of variation of power factor bσ: ; Compare the size of bQ and bσ, and adjust the reactive power and voltage phase angle of the power supply end; If bQ ≥ bσ, then increase the voltage phase angle at the power supply end by (1 + bσ) times; If bQ<bσ, the reactive power at the power supply end is adjusted to Quu.

6. The method for establishing an incremental model of controllable power resources according to claim 1, characterized in that: The specific steps to construct the admittance matrix are as follows: Count the number of power consumption areas ual in the first municipal district, obtain the topological relationship between the first municipal district and the power consumption areas, and construct the initial admittance matrix of the first municipal district; obtain the expected power consumption rate and reactive power of the ualth to ualth power consumption areas in the first municipal district, and obtain the complex power se of the first power consumption area in the first municipal district （1） ~se （ual） ; Get the power supply voltage ug from the 1st to the ualth power consumption area （1） ~ug （mn） ; Calculate the equivalent admittance Y of the first power consumption area （1,1） ~Y （1,ual） ; Determine whether the first power consumption area is directly connected to the second power consumption area; If they are connected, the equivalent admittance Y of the first power consumption area with respect to the second power consumption area is （1,2） : ; If they are not connected, the equivalent admittance Y of the first power consumption area with respect to the second power consumption area is （1,2） =0; Calculate the equivalent admittance Y of the first power consumption area with respect to the third to the ual power consumption areas （1,3） ~Y （1,ual） ; Calculate Y （1,3） ~Y （1,ual） The sum of the self-admittance Y of the first power consumption area （1,1） .

7. The method for establishing an incremental model of controllable power resources according to claim 6, characterized in that: The steps of constructing the admittance matrix also include: Assume that the complex power of the mth electricity consumption area is se （m） , the supply voltage is ug （m） ; The complex power of the nth electricity consumption area is se （n） , the supply voltage is ug （n） ; Assume that the equivalent admittance of the mth power consumption area with respect to the nth power consumption area is Y （m,n） , define calculation formula C-1-1: ； Assume that the equivalent admittance of the mth power consumption area with respect to the lth power consumption area is Y （m,l） ; Let the self-admittance of the mth power consumption area be Y （m,m） , define calculation formula C-1-2: ； Calculate the equivalent admittance and self-admittance Y from the second power consumption area to the ual power consumption area （2,l） ~Y （ual,ual） ; Construct the initial admittance matrix YY （1） : ； Extract Y （1,l） ~Y （ual,ual） The real part of the first power consumption area to the ual power consumption area is obtained. （1,l） ~G （ual,ual） ; Extract the imaginary part and get the equivalent susceptance B of the first to the ual power consumption area （1,l） ~B （ual,ual） ; Construct the initial admittance matrix of the 2nd to mnth districts; according to the matrix YY （1） , construct the calculation equation for the initial power flow calculation of the first urban district, and calculate the initial phase of the power supply voltage of each power consumption area in the first urban district.

8. The method for establishing an incremental model of controllable power resources according to claim 7, characterized in that: The specific steps for calculating the initial phase of the supply voltage are as follows: Assume that the supply voltage of the wth power consumption area is Uz （w） , the voltage phase is θ （w） , the expected power consumption is Pz （w） , reactive power is Qz （w） ; The complex power of the w-th power consumption area is Sz （w） , Sz （w） =Pz （w） +Qz （w） ×j; The supply voltage of the vth power consumption area is Ug （v） , the voltage phase is θ （v） ; The equivalent admittance of the wth power consumption area with respect to the vth power consumption area is Y （w，v） ; Y （w，v） The corresponding equivalent conductance is G （w,v） , the equivalent reactance is B （w,v） ; Assume that the conjugate of the supply voltage of the wth power consumption area is (Uz （w） ) * , define formula C-2-1: ; The conjugate of the supply voltage of the vth power consumption area is (Uz （v） ) * , ; Define formula C-2-2: ； Expanding formula C-2-2, we get formula C-2-3: ； Among them, θ （w，v） Indicates the phase difference between the supply voltage of the wth power consumption area and the vth power consumption area; Construct Matrix I （1） and matrix U （1） ;Define formula C-3: ； Assume that the initial phase of the supply voltage from the 1st to the ualth power consumption area is θz （1） ~θz （mn） ; Arrange formula C-3 to get formula D （1） ~Formula D （ual） : Formula D （1） for: ； Similarly, Formula D （ual） for: ； Using the Newton-Raphson algorithm, calculate θz （1） ,θz （2） ~θz （mn） The value of is used to obtain the initial phase of the supply voltage of the first to the ual power consumption area; Calculate the initial phase of the corresponding power consumption areas of the 2nd to mnth municipal districts; obtain the actual power consumption of all power consumption areas in the first district, and adjust the power supply power of the 1st to ualth power consumption areas.

9. The method for establishing an incremental model of controllable power resources according to claim 8, characterized in that: The specific steps for adjusting the power supply to the first to the ualth power consumption area are as follows: Get the expected power consumption rate po of the first to the ual power consumption area （1） ~po （mn） , actual power consumption rate pa （1） ~pa （mn） ; Determine whether the expected power consumption rate is equal to the actual power consumption rate, and adjust the power supply voltage and phase from each municipal district to the corresponding power consumption district based on the relationship between the actual power consumption and the expected power consumption; In the first to the ualth power consumption zones, if there is only one power consumption zone whose expected power consumption rate is not equal to the actual power consumption rate, the power consumption zone with the different expected power consumption rate and the actual power consumption rate is regarded as an abnormal zone; Obtain the expected power consumption rate pox, actual power consumption rate pax, actual supply voltage ux, and actual supply voltage phase θx of the abnormal area; Obtain the expected power supply voltage uy and the initial power supply voltage phase θy in the abnormal area; Determine whether ux and uy, as well as θx and θy, are equal, and adjust the power supply voltage and power supply voltage phase in the abnormal area or other areas; If ux and uy are equal, it means that the reactive power in the abnormal area has changed; Compare the sizes of pox and pax, and adjust the power supply voltage of other areas; If pox>pax, then reduce the power supply voltage of other areas times; If pox < pax, then increase the supply voltage of other areas times; If θx and θy are equal, it means that the active power in the abnormal area has changed; Compare the size of pox and pax, and adjust the phase of the power supply voltage in other areas; If pox>pax, then increase the phase of the power supply voltage in other areas times; If pox<pax, then reduce the phase of the power supply voltage in other areas times; If ux and uy as well as θx and θy are not equal, the abnormal area is used as the balance node and the matrix YY is reconstructed. （1） The actual power consumption and actual reactive power of all power consumption areas shall be adjusted.

10. The method for establishing an incremental model of controllable power resources according to claim 9, characterized in that: The specific steps to adjust the actual power consumption and actual reactive power of all power consumption areas are as follows: Assume that the abnormal area is numbered as r; recalculate the equivalent admittance rY of the abnormal area with respect to the 1st to ualth power consumption areas （r，1） ~rY （r，ual） , self-admittance rY （r，r） ; Calculate the change in the self-admittance of the abnormal area ΔY （r，r） ; Change in equivalent admittance ΔY （r，1） ; Calculate the changes in the equivalent admittance and self-admittance of the first to the ual power station, and obtain ΔY （1，1） ~ΔY （ual，ual） , get the matrix YY （2） ; Calculate the voltage change Δu at the balance node, ; Phase change Δθ, ; Calculate the change in actual power supply voltage from the 1st to the ualth power station Δu （1） ~Δu （ual） , the actual supply voltage phase change Δθ （1） ~Δθ （ual） ; Assume that the change in actual power consumption from the 1st to the ualth power station is ΔP （1） ~ΔP （ual） , the actual change in reactive power is ΔQ （1） ~ΔQ （ual） ; Construct matrix Zz （1） Matrix Zz （2） ; Define formula C-4-2: ; Based on Formula C-4-2, use the Newton-Raphson algorithm to calculate ΔP （1） ~ΔP （ual） and ΔQ （1） ~ΔQ （ual） The value of Increase the actual power consumption of the first to the ual power stations by ΔP （1） ~ΔP （ual） , the actual reactive power increases in turn by ΔQ （1） ~ΔQ （ual） ; In the first to the ualth power consumption areas, if the expected power consumption rates of multiple power consumption areas are not equal to the actual power consumption rates, the step of adjusting the actual power consumption and the actual reactive power of all power consumption areas is repeated.

Images