A controllable power resource incremental model establishment method

Abstract

This invention provides a method for establishing an incremental model of controllable power resources, belonging to the field of power regulation; it solves the problem of low efficiency in power resource regulation; specifically as follows: calculating the expected power consumption of each municipal district; constructing the 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 the expected power consumption of the power consumption area and constructing the admittance matrix; obtaining the actual power consumption of the power consumption area, and correcting the admittance matrix according to the expected power consumption, adjusting the supply voltage, phase, and power from each municipal district to the power consumption area; this invention, by acquiring, analyzing, and mathematically modeling relevant electrical data of the target area, simulates the power and voltage changes from transmission to consumption in each power consumption area of ​​the target area, which is helpful for the regulation of power resources in the target area.

Description

Technical Field

[0001] This invention relates to a method for establishing an incremental model of controllable power resources, and pertains to the field of power regulation. Background Technology

[0002] Existing methods for establishing incremental power resource models have the following shortcomings:

[0003] Model building is challenging: Existing methods for building incremental power resource models mostly involve the deep integration of multiple disciplines such as power electronics, communication 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 optimizing the admittance matrix depends on the analysis of complex power grid topologies. This interdisciplinary nature results in extremely high technical implementation thresholds, requiring the R&D team to have diverse professional backgrounds, making model building complex and data processing difficult.

[0004] Mismatch between static topology and dynamic power grid: Most existing methods for establishing incremental power resource models do not integrate dynamic topology identification modules, which leads to the inability to update the admittance matrix in real time during power grid reconfiguration (such as line faults or load switching); when a transmission line is disconnected due to a fault, but the model still assumes its existence, it may lead to power misallocation.

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

[0006] In view of the shortcomings of existing technologies, the purpose of this invention is to provide a method for establishing an incremental model of controllable power resources, which aims to solve the problem of low efficiency in power resource control.

[0007] To achieve the above objectives, the present invention provides a method for establishing an incremental model of controllable power resources, comprising:

[0008] Obtain the number of municipal districts in the target area, obtain the historical power consumption of each municipal district, perform time-series analysis on the historical power consumption, and calculate the expected power consumption of each municipal district.

[0009] Obtain the number of energy storage devices in the target area and construct the discharge characteristic function for each energy storage device; obtain the transmission power of the power supply end in the target area 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 municipality; if sufficient, adjust the discharge time of each energy storage device according to the discharge characteristic function; if insufficient, adjust the transmission power.

[0010] A time-series analysis of the historical power consumption of each municipal district to the power application area is performed to obtain the expected power consumption of each municipal district to the power application area, and an admittance matrix is ​​constructed. The actual power consumption of each municipal district to the power application area is obtained in real time and compared with the expected power consumption. Based on 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 power application area are adjusted.

[0011] Furthermore, the specific steps for constructing the discharge characteristic function of each energy storage device are as follows:

[0012] Count the number of energy storage devices ac, and construct the discharge characteristic function Qe of the first energy storage device. （1） ;

[0013] The average discharge power pd of the first energy storage device was obtained from day 1 to day mo of the past month. （1） ~pd （mo） Discharge time td （1） , td （2） ~td （mo） ; calories qd （1） ~qd （mo） Construct matrices X and Y;

[0014] Let the function Qe （1） The coefficient is β （0） ~β （5） Construct matrix Z;

[0015] Calculate β （0） ~β （5） Value: ;

[0016] Let the first discharge power be pe and the discharge time be te, and define the function Qe. （1） The mathematical expression:

[0017] ;

[0018] δ represents the correction value;

[0019] pd （1） ~pd （mo） As PE, td （1） ~td （mo）As te, qd （1） ~qd （mo） Substitute into function Qe （1） Calculate δ from the values;

[0020] Construct the discharge characteristic functions of the 2nd to acth energy storage devices to obtain the function Qe. （2） ~function Qe （ac） .

[0021] Furthermore, the specific steps for adjusting the discharge time and discharge duration of each energy storage device based on the discharge characteristic function are as follows:

[0022] Obtain and calculate the sum of expected electricity consumption (ael) for all municipal districts;

[0023] Obtain and calculate the sum of the remaining electricity of all energy storage devices;

[0024] Obtain the power transmission power Pu from the power supply end;

[0025] judge Is it valid?

[0026] If true, then no action is taken;

[0027] If not, then compare and Adjust the size of the energy storage device to control its discharge power and discharge time, or the power transmission power and reactive power at the power supply end.

[0028] like ≥ Then calculate the energy replenishment amount ree of the energy storage device; based on ree and the discharge characteristic function of each energy storage device, adjust the discharge power and discharge time of the energy storage device;

[0029] Calculate the discharge amount ec of the first to the acth energy storage devices. （1） ~ec （ac） ;

[0030] Let the discharge characteristic function of the i-th energy storage device be 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:

[0031] Let Qe （i） The mathematical expression is:

[0032] ;

[0033] β （0）（i） ~β （5） （i） Representation function Qe （i） The coefficient, δ （i） Representation function Qe （i） The correction value; the range of i is: 1 to ac;

[0034] Introducing the Lagrange multiplier λ, we construct equation L:

[0035] ;

[0036] Calculate equation L for pe （i） and te （i） The partial derivatives are then set to zero to eliminate λ, resulting in formula A-1:

[0037] ;

[0038] Simplifying formula A-1, we get formula A-2:

[0039] ;

[0040] Will Substituting into formula A-2, we get formula A-3:

[0041] ;

[0042] PE （i） As p, te （i） As t, will ec （i） As q, As a, As b, As c, As d, calculate pe （i） With te （i） The analytical expression.

[0043] Further, calculate pe （i） With te （i） The specific steps for obtaining the analytical expression are as follows:

[0044] Rearranging formula A-3, we get formula B-1:

[0045] ;

[0046] Introduce an auxiliary parameter o, let This leads to formula B-2:

[0047] ;

[0048] Where Aa, Bb, and Cc represent new parameters of p:

[0049] ;

[0050] By introducing the parameter k a second time, we can rewrite formula B-2 and obtain formula B-3:

[0051] ;

[0052] Define formula B-4:

[0053] ;

[0054] Calculate the smallest real root km of formula B-4;

[0055] Based on km, we obtain formula B-5:

[0056] ;

[0057] Based on formula B-5, the solution for o is obtained:

[0058] ;

[0059] Based on the solution for o, we obtain pe. （i） The analytical expression:

[0060] ;

[0061] According to PE （i） The analytical expression is obtained from te. （i） Analytical expression:

[0062] ;

[0063] Functions based on PE （i） and te （i） The analytical expression is used to calculate the power and discharge time of the first to the acth energy storage devices.

[0064] Furthermore, the step of adjusting the discharge time and discharge duration of each energy storage device according to the discharge characteristic function also includes:

[0065] like < Then, obtain the reactive power Qu at the power supply end of the target area, and adjust the transmission power and reactive power at the power supply end;

[0066] Calculate the additional power ΔPu at the power supply end:

[0067] ;

[0068] Calculate the power factor σ of the target region, σ = Pu / Qu;

[0069] Assuming σ remains constant, the new reactive power Qu is the result of an increase of ΔPu in transmission power: ;

[0070] Assuming Qu remains constant, the new power factor σu is the result of an increase in transmission power of ΔPu: ;

[0071] Calculate the reactive power variation coefficient bQ: ;

[0072] The power factor variation coefficient bσ: ;

[0073] Compare the magnitudes of bQ and bσ, and adjust the reactive power and voltage phase angle at the power supply end.

[0074] If bQ ≥ bσ, then the voltage phase angle at the power supply terminal will be increased by (1 + bσ) times;

[0075] If bQ < bσ, then the reactive power at the power supply end will be adjusted to Quu.

[0076] Furthermore, the specific steps for constructing the admittance matrix are as follows:

[0077] Count the number of electricity-consuming areas (u_a) within the first municipal district, obtain the topological relationship between the first municipal district and the electricity-consuming areas, and construct the initial admittance matrix of the first municipal district; obtain the expected power consumption rate and reactive power of the i-th to u_th electricity-consuming areas within the first municipal district, and obtain the complex power (se) of the first electricity-consuming area within the first municipal district. （1） ~se （ual） ;

[0078] Obtain the power supply voltage ug of the first to the ualth power consumption area. （1） ~ug （mn） ;

[0079] Calculate the equivalent admittance Y of the first power consumption area. （1,1） ~Y （1,ual） ;

[0080] Determine whether the first power consumption area is directly connected to the second power consumption area;

[0081] If connected, then the equivalent admittance Y of the first power consumption area with respect to the second power consumption area... （1,2） : ;

[0082] If they are not connected, then the equivalent admittance Y of the first power consumption area with respect to the second power consumption area... （1,2） =0;

[0083] Calculate the equivalent admittance Y of the first power consumption zone with respect to the third to the 1st power consumption zones. （1,3） ~Y（1,ual） ;

[0084] Calculate Y （1,3） ~Y （1,ual） The sum of these values, as the self-inductance Y of the first power consumption area. （1,1） .

[0085] Furthermore, the steps for constructing the admittance matrix also include:

[0086] Let the complex power of the m-th power consumption area be se. （m） The power supply voltage is ug. （m） The complex power of the nth power consumption area is se. （n） The power supply voltage is ug. （n） ;

[0087] Let Y be the equivalent admittance of the m-th power consumption area with respect to the n-th power consumption area. （m,n） Define the calculation formula C-1-1:

[0088] ;

[0089] Let Y be the equivalent admittance of the m-th power consumption area with respect to the l-th power consumption area. （m,l） Let the self-admittance of the m-th power consumption area be Y. （m,m） Define the calculation formula C-1-2:

[0090] ;

[0091] Calculate the equivalent admittance and self-admittance Y of the second to the 1ath power consumption areas. （2,l） ~Y （ual,ual） ;

[0092] Construct the initial admittance matrix YY （1） :

[0093] ;

[0094] Extract Y （1,l） ~Y （ual,ual） The real part is used to obtain the equivalent conductance G from the first to the 1st power consumption area. （1,l） ~G （ual,ual） Extract the imaginary part to obtain the equivalent susceptance B of the first to the 1st power consumption areas. （1,l） ~B （ual,ual） ;

[0095] Construct the initial admittance matrices for the 2nd to mnth municipal districts; based on matrix YY （1） The calculation equations for the initial power flow calculation of the first municipal district are constructed, and the initial phase of the power supply voltage of each power consumption area in the first municipal district is calculated.

[0096] Furthermore, the specific steps for calculating the initial phase of the supply voltage are as follows:

[0097] Let the power supply voltage of the w-th power consumption area be Uz. （w） The voltage phase is θ （w） The expected power consumption is Pz （w） The 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;

[0098] The power supply voltage for the vth power consumption area is Ug （v） The voltage phase is θ （v） ;

[0099] The equivalent admittance of the w-th power consumption area with respect to the v-th 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） ;

[0100] Let the conjugate of the power supply voltage of the w-th power consumption area be (Uz) （w） ) * Define formula C-2-1: ;

[0101] The conjugate of the power supply voltage for the v-th power consumption area is (Uz) （v） ) * , ;

[0102] Define formula C-2-2:

[0103] ;

[0104] Expanding formula C-2-2, we get formula C-2-3:

[0105] ;

[0106] Where, θ （w，v） This represents the phase difference between the power supply voltage of the w-th power consumption area and the v-th power consumption area;

[0107] Construct matrix I （1） Sum matrix U （1） ; Define formula C-3:

[0108] ;

[0109] Let the initial phase of the power supply voltage for the first to the 1st power consumption area be θz. （1） ~θz（mn） ;

[0110] Rearranging formula C-3, we get formula D. （1） Formula D （ual） :

[0111] Formula D （1） for:

[0112] ;

[0113] Similarly, formula D （ual） for:

[0114] ;

[0115] Calculate θz using Newton-Raphson's algorithm. （1） , θz （2） ~θz （mn） The value of is used to obtain the initial phase of the power supply voltage for the first to the ualth power consumption areas;

[0116] Calculate the initial phase of the application power zone for the 2nd to mnth municipal districts; obtain the actual power consumption of all power zones in the 1st district, and adjust the power supply of the 1st to ualth power zones.

[0117] Furthermore, the specific steps for adjusting the power supply of the first to the 11th power consumption areas are as follows:

[0118] Obtain the expected electricity consumption rate po for the first to the 1st electricity consumption area. （1） ~po （mn） Actual electricity consumption rate (pa) （1） ~pa （mn） ;

[0119] Determine whether the expected electricity consumption rate is equal to the actual electricity consumption rate, and adjust the power supply voltage and phase of each municipal district to the corresponding electricity consumption area based on the relationship between the actual power consumption and the expected power consumption.

[0120] In the first to the 11th power consumption area, if only one power consumption area has an expected power consumption rate that is not equal to the actual power consumption rate, then the power consumption area with the expected power consumption rate that is not equal to the actual power consumption rate is regarded as an abnormal area.

[0121] Obtain the expected power consumption rate pox, actual power consumption rate pax, actual power supply voltage ux, and actual power supply voltage phase θx for the abnormal area;

[0122] Obtain the expected supply voltage uy in the abnormal region and the initial supply voltage phase θy;

[0123] Determine whether ux and uy, as well as θx and θy, are equal, and adjust the power supply voltage and phase in abnormal areas or other areas accordingly.

[0124] If ux and uy are equal, it indicates that the reactive power in the abnormal region has changed;

[0125] Compare the values ​​of POX and PAX, and adjust the power supply voltage for other areas;

[0126] If pox > pax, then reduce the power supply voltage in other areas. times;

[0127] If pox < pax, then increase the power supply voltage in other areas. times;

[0128] If θx and θy are equal, it indicates that the active power in the abnormal region has changed;

[0129] Compare the values ​​of POX and PAX, and adjust the phase of the power supply voltage in other areas;

[0130] If pox > pax, then increase the phase of the power supply voltage in other zones. times;

[0131] If pox < pax, then reduce the phase of the power supply voltage in other zones. times;

[0132] If ux and uy, as well as θx and θy, are all unequal, then the outlier region is used as the balancing node to reconstruct matrix YY. （1） Adjust the actual power consumption and actual reactive power of all power consumption areas.

[0133] Furthermore, the specific steps for adjusting the actual power consumption and actual reactive power of all power consumption areas are as follows:

[0134] Let the abnormal zone be numbered r; recalculate the equivalent admittance rY of the abnormal zone with respect to the 1st to the 1st power consumption zone. （r，1） ~rY （r，ual） Self-guided nanorY （r，r） ;

[0135] Calculate the change in self-admittance ΔY in the abnormal region. （r，r） The change in equivalent admittance ΔY （r，1） ;

[0136] Calculate the changes in equivalent admittance and self-admittance for the first to the 11th power plants to obtain ΔY. （1，1） ~ΔY （ual，ual） , thus obtaining matrix YY （2） ;

[0137] Calculate the voltage change Δu at the slack node. Phase change Δθ ;

[0138] Calculate the change in actual power supply voltage Δu for the first to the 1st power station. （1） ~Δu （ual） The actual change in phase of the supply voltage Δθ （1） ~Δθ （ual） ;

[0139] Let ΔP be the change in actual power consumption of the first to the 1st power station. （1） ~ΔP （ual） The actual change in reactive power is ΔQ （1） ~ΔQ （ual） ;

[0140] Construct matrix Zz （1） : ; matrix Zz （2） : ;

[0141] Define formula C-4-2: ;

[0142] Based on formula C-4-2, the Newton-Raphson algorithm is reused to calculate ΔP. （1） ~ΔP （ual） and ΔQ （1） ~ΔQ （ual） The value;

[0143] The actual power consumption of the first to the 11th power stations will be increased by ΔP sequentially. （1） ~ΔP （ual） The actual reactive power increases by ΔQ sequentially. （1） ~ΔQ （ual） ;

[0144] If, in the first to the 11th power consumption zones, the expected power consumption rate of multiple power consumption zones is not equal to the actual power consumption rate, then the steps of adjusting the actual power consumption and actual reactive power of all power consumption zones are repeated.

[0145] Compared with the prior art, the beneficial effects of the present invention are:

[0146] Multi-dimensional data fusion: This invention incorporates heterogeneous data such as historical power consumption in the city's jurisdiction, discharge characteristics of energy storage devices, and topological relationships of power consumption areas into a unified framework, breaking through the limitations of traditional power regulation that relies on a single dimension; by comparing the expected power consumption with the 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; and realize bidirectional regulation between the transmission end and the power consumption end.

[0147] The hierarchical control architecture adopts a three-level control structure of "city-municipal district-electricity consumption area", which takes into account both global optimization and local fine control. At the city level, the energy storage device and the power supply end are coordinated. At the municipal district level, the demand is predicted by the ARIMA model. At the electricity consumption area level, the power supply parameters are adjusted by the admittance matrix, forming a hierarchical decision-making system. At the same time, this invention realizes quantitative control of the discharge process by constructing the discharge power-time-heat generation relationship (Qe(i) function) of the energy storage device. The discharge power (pe(i)) and time (te(i)) are optimized by using the Lagrange multiplier method, so as to minimize heat loss while meeting the discharge quantity (ec(i)) requirement.

[0148] Flexibility and robustness of power consumption-side regulation: By monitoring the deviation between the actual power consumption and the expected value in a power consumption area in real time, this invention can construct and adjust the admittance matrix parameters to optimize the power supply voltage and phase. When a power shortage occurs in a power consumption area, cross-regional power support can be achieved by reconstructing the matrix balance node. By adopting the "balance node" strategy, abnormal power consumption areas are isolated and reconstructed 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. Attached Figure Description

[0149] Other features, objects, and advantages of the invention will become more apparent from the following detailed description of non-limiting embodiments with reference to the accompanying drawings:

[0150] Figure 1 This is a schematic diagram of the method of the present invention;

[0151] Figure 2 This is a schematic diagram of the power consumption area topology of the present invention;

[0152] Figure 3 This is a schematic diagram of the power consumption area topology of the present invention. Detailed Implementation

[0153] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0154] Please see Figure 1 A method for establishing an incremental model of controllable power resources includes:

[0155] Step S1: Obtain the number of municipal districts in the target area, obtain the historical power consumption of each municipal district, perform time-series analysis on the historical power consumption, and calculate the expected power consumption of each municipal district;

[0156] The historical electricity consumption data for each municipal district was obtained by collecting hourly data for the past three days. Time series analysis was performed on the historical electricity consumption data to calculate the difference order, autoregressive parameters, and moving average parameters for each municipal district. Based on the difference order, autoregressive parameters, and moving average parameters for each municipal district, the autocorrelation coefficient and moving average coefficient for each municipal district were calculated. Finally, the ARIMA model was used to predict the total electricity consumption for each municipal district in the next day, based on the autocorrelation coefficient and moving average coefficient, as the expected electricity consumption for each municipal district.

[0157] It should be noted that, in this invention, "target area" refers to the city-level region where incremental modeling and regulation of power resources are carried out using this invention (a method for establishing an incremental model of controllable power resources).

[0158] Step S2: (Power Transmission Control) Obtain the number of energy storage devices in the target area and construct the discharge characteristic function for each energy storage device; obtain the power transmission capacity at the power supply end of the target area and the remaining power of each energy storage device, and determine whether the power transmission capacity is sufficient based on the expected power consumption of each municipality; if sufficient, adjust the discharge time of each energy storage device according to the discharge characteristic function; if insufficient, adjust the power transmission capacity.

[0159] It should be noted that in this invention, "energy storage device" refers to an electrical energy storage device; "power supply terminal" refers to a power plant or national grid that supplies power to the target area.

[0160] The specific steps of step S2 are as follows:

[0161] Count the number of energy storage devices (ac) within the target area, construct the discharge characteristic functions of the 1st to acth energy storage devices, and obtain the function Qe. （1） ~function Qe （ac） ;

[0162] Construct the discharge characteristic function Qe of the first energy storage device （1） ;

[0163] The average discharge power pd of the first energy storage device was obtained on the first, second, and mo days of the past month. （1） ,pd （2） ~pd （mo） ;

[0164] Discharge time td （1） , td （2） ~td （mo） ;

[0165] The heat qd generated by the first energy storage device due to discharge. （1） qd （2） ~qd （mo）Where mo represents the number of days in the past month;

[0166] Construct a matrix X of size (mo×6):

[0167] ;

[0168] Construct a (mo×1) matrix Y: ;

[0169] Let the function Qe （1） The coefficient is β （0） β （1） ~β （5） ;

[0170] Construct a (6×1) matrix Z: ;

[0171] Calculate β （0） ~β （5） Value: Where T represents the transpose of the matrix, and -1 represents the inverse of the matrix;

[0172] Let the first discharge power be pe and the discharge time be te, and define the function Qe. （1） The mathematical expression:

[0173] ;

[0174] Where δ represents the function Qe （1） The correction value (the initial value of δ is 0);

[0175] pd （1） ~pd （mo） As PE, td （1） ~td （mo） As te, qd （1） ~qd （mo） As a function Qe （1） Substitute the value into the function Qe （1） In the mathematical expression, the value of δ is calculated to obtain the function Qe. （1） ;

[0176] Repeatedly construct function Qe （1） Following the same steps, construct the discharge characteristic functions of the 2nd to acth energy storage devices to obtain the function Qe. （2） ~function Qe （ac） ;

[0177] Obtain the expected electricity consumption of the first to the mnth municipal districts, and get el （1） ~el （mn） Where mn represents the number of municipal districts; calculate el （1） ~el（mn） and ael;

[0178] Obtain the remaining power of the 1st, 2nd, and up to the acth energy storage devices, and get re. （1） re （2） ~re （ac） ; calculate re （1） ~re （ac） and are;

[0179] Obtain the transmission power Pu at the power supply end of the (target area);

[0180] judge Is it valid?

[0181] If true, then no action is taken;

[0182] If not, then compare and Adjust the size of the energy storage device to control its discharge power and discharge time, or the power transmission power and reactive power at the power supply end.

[0183] like ≥ Then calculate the electrical energy replenishment amount ree of the energy storage device: ;

[0184] Based on ree and the discharge characteristic function of each energy storage device, adjust the discharge power and discharge time of the energy storage device;

[0185] Calculate the discharge quantity ec of the 1st, 2nd, up to the acth energy storage device. （1） ec （2） ~ec （ac） ;

[0186] Among them, ec （1） The formula for calculation is: ;

[0187] ec （2） The formula for calculation is: ;

[0188] And so on, ec （ac） The formula for calculation is: ;

[0189] Let the discharge characteristic function of the i-th energy storage device be 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:

[0190] Let Qe（i） The mathematical expression is:

[0191] ;

[0192] Where, β （0） （i） ~β （5） （i） Representation function Qe （i） The coefficient, δ （i） Representation function Qe （i） The correction value; the range of i is: 1 to ac;

[0193] Define formula A-1: ;ec （i） This represents the discharge amount of the i-th energy storage device;

[0194] Introducing the Lagrange multiplier λ, we construct equation L:

[0195] ;

[0196] Calculate equation L for pe （i） and te （i） The partial derivatives of and set them to zero, yield formula A-2:

[0197] ;

[0198] By combining equations A-2 and eliminating λ, we obtain equation A-3:

[0199] ;

[0200] Simplifying formula A-3, we get formula A-4:

[0201] ;

[0202] Based on formula A-1, it can be seen that if and only if At that time, the i-th energy storage device has the best discharge effect (i.e., the lowest discharge heat generation).

[0203] Will Substituting into formula A-4, we get formula A-5:

[0204] ;

[0205] PE （i） As p, te （i） As t, will ec （i） As q, As a, As b, As c, As d;

[0206] PE is calculated based on formula A-5. （i） With te （i） The analytical expression:

[0207] Rearranging formula A-5, we get formula B-1:

[0208] ;

[0209] Introduce an auxiliary parameter o, let Eliminating the cubic term of p in formula B-1 yields formula B-2:

[0210] ;

[0211] Where Aa, Bb, and Cc represent new parameters of p:

[0212] ;

[0213] By introducing the parameter k a second time, formula B-2 is rewritten as a perfect square, resulting in formula B-3:

[0214] ;

[0215] For the right-hand side of formula B-3 to also satisfy a perfect square, the discriminant of the right-hand side of formula B-3 must be zero, resulting in formula B-4-1:

[0216] ;

[0217] Expanding formula B-4-1, we obtain a cubic equation in k, which serves as formula B-4-2:

[0218] ;

[0219] Let kz be the smallest real root of formula B-4-2, and define formula S to test the value of kz; the expression of formula S is:

[0220] ;

[0221] Using Cardano's formula (i.e., the formula for finding the roots of a cubic equation), calculate formula B-4-2 and extract the smallest real root km that satisfies the S test of formula;

[0222] Based on km, we can rearrange formula B-3 to obtain formula B-5:

[0223] ;

[0224] Based on formula B-5, the solution for o is obtained:

[0225] ;

[0226] Based on the solution for o, we obtain pe. （i） The analytical expression:

[0227] ;

[0228] According to PE （i） The analytical expression is obtained from te. （i） Analytical expression:

[0229] ;

[0230] Functions based on PE （i） and te （i） Using the analytical expression, calculate the power and discharge time of the first to the acth energy storage devices;

[0231] like < Then, obtain the reactive power Qu at the power supply end of the target area, and adjust the transmission power and reactive power at the power supply end;

[0232] Calculate the additional power ΔPu at the power supply end:

[0233] ;

[0234] Calculate the power factor σ of the target region, σ = Pu / Qu;

[0235] Assuming σ remains constant, the new reactive power Qu is the result of an increase of ΔPu in transmission power: ;

[0236] Assuming Qu remains constant, the new power factor σu is the result of an increase in transmission power of ΔPu: ;

[0237] Calculate the reactive power variation coefficient bQ: ;

[0238] The power factor variation coefficient bσ: ;

[0239] Compare the magnitudes of bQ and bσ, and adjust the reactive power and voltage phase angle at the power supply end.

[0240] If bQ ≥ bσ, then the voltage phase angle at the power supply terminal will be increased by (1 + bσ) times;

[0241] If bQ < bσ, then the reactive power at the power supply end will be adjusted to Quu.

[0242] Step S3: (Power Consumption Control) Perform time-series analysis on 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 the 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, correct the admittance matrix and adjust the power supply voltage, phase and power of each municipal district to the power application area;

[0243] It should be noted that, in this invention, "electricity consumption area" refers to the area within each municipal district that consumes electricity, such as residential areas, commercial areas, industrial areas, public facilities areas, etc.

[0244] The specific steps of step S3 are as follows:

[0245] Step S31: Count the number of electricity-consuming areas in the 1st, 2nd, up to the mnth municipal district, and denote it as ua. （1） ua （2） ~ua （mn） ;

[0246] Obtain the hourly power consumption of each municipal district in the power application area for the past three days, perform time-series analysis on the historical power consumption of each municipal district in the power application area, and obtain the expected power consumption of each municipal district in the power application area.

[0247] Obtain the topological relationship between the first municipal district and the first, second, up to the ualth power consumption area, and construct the initial admittance matrix of the first municipal district;

[0248] ua （1） As ual, obtain the expected electricity consumption rate po of the 1st, 2nd, up to ualth electricity consumption areas within the 1st municipal district. （1） ,po （2） ~po （ual） ;

[0249] Obtain the (average) reactive power qo of the first, second, up to the ualth electricity consumption area in the first municipal district yesterday. （1） ,po （2） ~po （ual） ;

[0250] Based on the expected electricity consumption rate and reactive power, combine the complex power se of the first electricity consumption area within the first municipal district. （1） ,se （1） =po （1） +qo （1） ×j;

[0251] The combined power of the second power consumption area is se （2） ,se （2） =po （2） +qo（2） ×j;

[0252] And so on, the complex power of the ualth power consumption area se （ual） ,se （ual） =po （ual） +qo （ual） ×j; where j represents the imaginary unit;

[0253] Obtain the power supply voltage ug from the first municipal district (yesterday) to the first, second, and so on, up to the ualth power consumption area. （1） ug （2） ~ug （mn） ;

[0254] Calculate the equivalent admittance Y of the first power consumption area. （1,1） ~Y （1,ual） ;

[0255] Determine whether the first power consumption area is directly connected to the second power consumption area;

[0256] Please see 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） : ;

[0257] Please see 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;

[0258] (Repeat Y) （1,2） (Calculation process) Calculate the equivalent admittance Y of the first power consumption area with respect to the third to 1ath power consumption areas. （1,3） ~Y （1,ual） ;

[0259] Calculate Y （1,3） ~Y （1,ual） The sum of these values, as the self-inductance Y of the first power consumption area. （1,1） ;

[0260] Let the complex power of the m-th power consumption area be se. （m） The supply voltage (for the first municipal district to supply power to the m-th power consumption area) is ug. （m） The complex power of the nth power consumption area is se. （n） The supply voltage (for the first municipal district to supply power to the nth power-consuming area) is ug. （n） ;

[0261] Let Y be the equivalent admittance of the m-th power consumption area with respect to the n-th power consumption area. （m,n） Define the calculation formula C-1-1:

[0262] ;

[0263] The range of values ​​for m is 2 to (ual-1).

[0264] The range of n is (m+1) to ual;

[0265] Let Y be the equivalent admittance of the m-th power consumption area with respect to the l-th power consumption area. （m,l） The value range of l is 1 to (m-1); let the self-admittance of the m-th power consumption area be Y. （m,m） Define the calculation formula C-1-2:

[0266] ;

[0267] Based on formulas C-1-1 and C-1-2, calculate the equivalent admittance Y of the second power consumption zone with respect to the first to the 1st power consumption zone. （2,l） ~Y （2,ual） and self-guided Y （2,2） ;

[0268] The equivalent admittance Y of the third power consumption zone with respect to the first to the last ual （3,l） ~Y （3,ual） and self-guided Y （3,3） ;

[0269] Similarly, the equivalent admittance Y of the ual-th power consumption area with respect to the 1st to (ual-1)th power consumption area is... （ual,l） ~Y （ual,（ual-1）） and self-guided Y （ual,ual） ;

[0270] Construct a (ual×ual) zero matrix and fill it into Y. （1,l） ~Y （ual,ual） The zero matrix yields the initial admittance matrix YY of the first municipal district. （1） :

[0271] ;

[0272] Among them, matrix YY （1） The diagonal elements are the self-admittances of the 1st to the 1st power consumption region, and the off-diagonal elements are equal in value to the "elements symmetric about the diagonal" (similar to (-Y)). （1,2） ) = (-Y （2,1） (-Y) （ual,1） ) = (-Y （1,ual） ));

[0273] Extract Y （1,l） Y （1,2） ~Y （1,ual） The real part is used to obtain the equivalent conductance G of the first power consumption area. （1,l） G （1,2） ~G（1,ual） Extract Y （1,l） Y （1,2） ~Y （2,ual） The imaginary part is used to obtain the equivalent susceptance B of the first power consumption area. （1,l） B （1,2） ~B （1,ual） ;

[0274] Extract Y （1,l） Y （1,2） ~Y （1,ual） The real part is used to obtain the equivalent conductance G of the first power consumption area. （1,l） G （1,2） ~G （1,ual） Extract Y （1,l） Y （1,2） ~Y （2,ual） The imaginary part is used to obtain the equivalent susceptance B of the first power consumption area. （1,l） B （1,2） ~B （1,ual） ;

[0275] Extract Y （2,l） Y （2,2） ~Y （1,ual） The real part is used to obtain the equivalent conductance G of the second power consumption area. （2,l） G （2,2） ~G （2,ual） Extract Y （2,l） Y （2,2） ~Y （2,ual） The imaginary part is used to obtain the equivalent susceptance B of the first power consumption area. （2,l） B （2,2） ~B （2,ual） ;

[0276] And so on, extract Y （ual,l） Y （ual,2） ~Y （ual,ual） The real part is used to obtain the equivalent conductance G of the first power consumption area. （ual,l） G （ual,2） ~G （ual,ual） Extract Y （ual,l） Y （ual,2） ~Y （ual,ual） The imaginary part is used to obtain the equivalent susceptance B of the first power consumption area. （ual,l） B （ual,2） ~B （ual,ual） ;

[0277] Repeating matrix YY （1） The construction steps are as follows: construct the admittance matrix of the 2nd to mnth municipal districts;

[0278] According to matrix YY （1）Construct the calculation equations for the initial power flow calculation of the first municipal district, and calculate the initial phase of the power supply voltage for each power consumption area within the first municipal district;

[0279] Let the power supply voltage of the w-th power consumption area (within the first municipal district) be Uz. （w） The voltage phase is θ （w） The expected power consumption is Pz （w） The 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;

[0280] The power supply voltage for the vth power consumption area is Ug （v） The voltage phase is θ （v） ;

[0281] The equivalent admittance of the w-th power consumption area with respect to the v-th 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） ; where the values ​​of w and v are both in the range of 1 to ual;

[0282] Let the conjugate of the power supply voltage of the w-th power consumption area be (Uz) （w） ) * Define formula C-2-1: ;

[0283] The conjugate of the power supply voltage for the v-th power consumption area is (Uz) （v） ) * , ;

[0284] Define formula C-2-2:

[0285] ;

[0286] Expanding formula C-2-2, we get formula C-2-3:

[0287] ;

[0288] Where, θ （w，v） This represents the phase difference between the power supply voltage of the w-th power consumption area and the v-th power consumption area;

[0289] Construct a (ual×1) matrix I （1） : ;

[0290] Construct a (ual×1) matrix U （1） : ;

[0291] Define formula C-3-1: ;

[0292] Expanding and rearranging formula C-3-1, we get formula C-3-2:

[0293] ;

[0294] Let the initial phase of the power supply voltage for the 1st, 2nd, and up to the ualth power consumption area be θz. （1） , θz （2） ~θz （mn） ;

[0295] po （1） ~po （mn） and qo （1） ~po （mn） Substituting into formula C-3-2 and rearranging, we obtain formula D. （1） Formula D （2） Formula D （ual） :

[0296] Formula D （1） for:

[0297] ;

[0298] Formula D （2） for:

[0299] ;

[0300] And so on, formula D （ual） for:

[0301] ;

[0302] Calculate θz using Newton-Raphson's algorithm. （1） , θz （2） ~θz （mn） The value of is used to obtain the initial phase of the power supply voltage for the 1st, 2nd, and up to the ualth power consumption area (within the 1st municipal district);

[0303] Repeat θz （1） ~θz （mn） The calculation process calculates the initial phase of the application power zone for the 2nd to mnth municipal districts;

[0304] Obtain the actual power consumption of all power consumption areas within the first jurisdiction, and adjust the power supply of the first to the 1st power consumption areas according to the expected power consumption of each power consumption area;

[0305] Obtain the expected electricity consumption rate po for the 1st, 2nd, up to the ualth electricity consumption area within the 1st municipal district. （1） ,po （2） ~po （mn） ;

[0306] Obtain the actual electricity consumption rate pa of the first, second, up to the ualth electricity consumption area within the first municipal district. （1） pa （2） ~pa （mn） ;

[0307] Determine whether the expected power consumption rate and the actual power consumption rate of the first to the 11th power consumption areas are equal, and adjust the power supply voltage and phase of each municipal district to the corresponding power consumption area based on the relationship between the actual power consumption and the expected power consumption.

[0308] In the first to the 11th power consumption area, if only one power consumption area has an expected power consumption rate that is not equal to the actual power consumption rate, then the power consumption area with the expected power consumption rate that is not equal to the actual power consumption rate is regarded as an abnormal area.

[0309] Obtain the expected power consumption rate pox in the abnormal region, where pox∈{po （1） ,po （2） ~po （ual）}; Actual electricity consumption rate pax, pax∈{pa （1） pa （2） ~pa （mn）};

[0310] Obtain the actual power supply voltage ux in the abnormal region, ux∈{ug （1） ug （2） ~ug （ual）}; Actual supply voltage phase θx, θx∈{θz} （1） , θz （2） ~θz （ual）};

[0311] Obtain the expected supply voltage uy in the abnormal region and the initial supply voltage phase θy;

[0312] Determine whether ux and uy and θ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 power consumption areas other than the abnormal area in the first municipal district from the first to the ualth power consumption area).

[0313] If ux and uy are equal, it indicates that the reactive power in the abnormal region has changed;

[0314] Compare the values ​​of POX and PAX, and adjust the power supply voltage for other areas;

[0315] If pox > pax, then reduce the power supply voltage in other areas. times;

[0316] If pox < pax, then increase the power supply voltage in other areas. times;

[0317] If θx and θy are equal, it indicates that the active power in the abnormal region has changed;

[0318] Compare the values ​​of POX and PAX, and adjust the phase of the power supply voltage in other areas;

[0319] If pox > pax, then increase the phase of the power supply voltage in other zones. times;

[0320] If pox < pax, then reduce the phase of the power supply voltage in other zones. times;

[0321] If ux and uy, as well as θx and θy, are all unequal, then the outlier region is used as the balancing node to reconstruct matrix YY. （1） Adjust the actual power consumption and actual reactive power of all power consumption areas;

[0322] Let the number of the abnormal region be r, and the value range of r be 1 to u1;

[0323] Based on 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 1st, 2nd, and up to the ualth power consumption area. （r，1） ,rY （r，2） ~rY （r，ual） Self-guided nanorY （r，r） ;

[0324] Get matrix YY （1） The equivalent admittance tY of the first, second, and up to the ualth power consumption areas in the abnormal region. （r，1） ,tY （r，2） ~tY （r，ual） Self-guided nanorY （r，r） ;

[0325] Calculate the change in self-admittance ΔY in the abnormal region. （r，r） , ;

[0326] Calculate the change in equivalent admittance in the abnormal region to obtain ΔY. （r，1） ΔY （r，2） ~ΔY （r，ual） ;

[0327] in, , And so on, ;

[0328] The equivalent admittance of the w-th power consumption area with respect to the abnormal area is rY （w，r） The self-admittance of the wth power consumption area is rY. （w，w） ;

[0329] Let ΔY be the change in self-admittance of the w-th power consumption area. （w，w） Equivalent admittance Y （w，v） The change is ΔY （w，v） Define the calculation formula C-4-1:

[0330] ;

[0331] According to formula C-4-1, calculate the changes in equivalent admittance and self-admittance of the 1st, 2nd, and up to the ualth power station, and obtain ΔY. （1，1） ΔY （1，2） ~ΔY （ual，ual） ;

[0332] matrix YY （1） Replace the middle element with ΔY （1，1） ΔY （1，2） ~ΔY （ual，ual） , thus obtaining matrix YY （2） ;

[0333] Treating the abnormal region as a balancing node, calculate the voltage change Δu at the balancing node. Phase change Δθ ;

[0334] Obtain the actual power supply voltage uv of the 1st, 2nd, up to the ualth power station (excluding abnormal areas). （1） uv （2） ~uv （ual） Phase θv of the actual supply voltage （1） θv （2） ~θv （ual） ;

[0335] Calculate the change in actual power supply voltage Δu for the 1st, 2nd, and up to the ualth power station (excluding abnormal areas). （1） , Δu （2） ~Δu （ual） The actual change in phase of the supply voltage Δθ （1） , Δθ （2） ~Δθ （ual） ;

[0336] in, , And so on. ;

[0337] Let ΔP be the change in actual power consumption of the 1st, 2nd, and up to the ualth power station (excluding abnormal areas). （1） ΔP （2） ~ΔP （ual） The actual change in reactive power is ΔQ （1） ΔQ （2） ~ΔQ （ual） ;

[0338] Construct two (ual×2) matrices to obtain matrix Zz. （1） And matrix Zz （2） ;

[0339] Matrix Zz （1） : ; matrix Zz （2） : ;

[0340] Define formula C-4-2: ;

[0341] Based on formula C-4-2, repeat formula D. （1） Formula D （ual） The expansion method is described, and the Newton-Raphson algorithm is used to calculate ΔP. （1） ΔP （2） ~ΔP （ual） and ΔQ （1） ΔQ （2） ~ΔQ （ual） The value;

[0342] The actual power consumption of the 1st, 2nd, and so on up to the ualth power station (excluding abnormal areas) will be increased by ΔP sequentially. （1） ΔP （2） ~ΔP （ual） The actual reactive power increases by ΔQ sequentially. （1） ΔQ （2） ~ΔQ （ual） ;

[0343] In the first to the 11th power consumption area, if the expected power consumption rate of multiple power consumption areas is not equal to the actual power consumption rate, then select the power consumption area with fewer directly connected power consumption areas and a smaller difference between the actual power consumption and the expected power consumption as the balancing node, and repeat the steps of adjusting the actual power consumption and actual reactive power of all power consumption areas to adjust the actual power consumption and actual reactive power of all power consumption areas.

[0344] Repeat the adjustment steps for the power supply voltage, power supply voltage phase, actual power consumption and reactive power of all power consumption areas in the first municipal district, and then adjust the power consumption areas corresponding to the second to mn municipal districts.

[0345] The above formulas are all dimensionless calculations. The formulas are derived from software simulations using a large amount of collected data to obtain the most recent real-world results. The preset parameters in the formulas are set by those skilled in the art according to the actual situation. For example, there are weighting coefficients and proportional coefficients. The values ​​set are to quantify each parameter to obtain a specific value, which is convenient for subsequent comparison. The values ​​of the weighting coefficients and proportional coefficients are only required to not affect the proportional relationship between the parameters and the quantified values.

[0346] Finally, it should be noted that the above-described embodiments are merely specific implementations of the present invention, used to illustrate the technical solutions of the present invention, and not to limit it. 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 foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments within the technical scope disclosed in the present invention, or make equivalent substitutions for some of the technical features; and these modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be covered within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for establishing an incremental model of controllable power resources, characterized in that, The method includes: Obtain the number of municipal districts in the target area, obtain the historical electricity consumption of each municipal district, perform time-series analysis on the historical electricity consumption, and calculate the expected electricity consumption of each municipal district; Obtain the number of energy storage devices in the target area and construct the discharge characteristic function for each energy storage device; obtain the transmission power of the power supply end in the target area 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 power and discharge time of each energy storage device according to the discharge characteristic function; if insufficient, adjust the transmission power. Count the number of energy storage devices ac, construct the discharge characteristic functions 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） ; The average discharge power pd of the first energy storage device was obtained 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 qd generated by the first energy storage device due to discharge. （1） qd （2） ~qd （mo） Where mo represents the number of days in the past month; Construct a matrix X of size (mo×6): ； Construct a (mo×1) matrix Y: ; Let the function Qe （1） The coefficient is β （0） β （1） ~β （5） ; Construct a (6×1) matrix Z: ; Calculate β （0） ~β （5） Value: Where T represents the transpose of the matrix, and -1 represents the inverse of the matrix; Let the first discharge power be pe and the discharge time be te, and define the function Qe. （1） The mathematical expression: ； Where δ represents the function Qe （1） The correction value (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 Qe （1） Substitute the value into the function Qe （1） In the mathematical expression, the value of δ is calculated to obtain the function Qe. （1） ; Repeatedly construct function Qe （1） Following the same steps, construct the discharge characteristic functions of the 2nd to acth energy storage devices to obtain the function Qe. （2） ~function Qe （ac） ; A time-series analysis of the historical power consumption of each municipal district to the power application area is performed to obtain the expected power consumption of each municipal district to the power application area, and an admittance matrix is ​​constructed. The actual power consumption of each municipal district to the power application area is obtained in real time and compared with the expected power consumption. Based on 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 power application 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 adjusting the power and discharge time of each energy storage device based on the discharge characteristic function are as follows: Obtain and calculate the sum of expected electricity consumption (ael) for all municipal districts; Obtain and calculate the sum of the remaining electricity of all energy storage devices; Obtain the power transmission power Pu from the power supply end; judge Is it valid? If true, then no action is taken; If not, then compare and Adjust the size of the energy storage device to control its discharge power and discharge time, or the power transmission power and reactive power at the power supply end. like ≥ Then calculate the energy replenishment amount ree of the energy storage device; based on 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 amount ec of the first to the acth energy storage devices. （1） ~ec （ac） ; Let the discharge characteristic function of the i-th energy storage device be 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: Let Qe （i） The mathematical expression is: ； β （0） （i） ~β （5） （i） Representation function Qe （i） The coefficient, δ （i） Representation function Qe （i） The correction value; the range of i is: 1 to ac; Introducing the Lagrange multiplier λ, we construct equation L: ； Calculate equation L for pe （i） and te （i） The partial derivatives are then set to zero to eliminate λ, resulting in formula A-1: ； Simplifying formula A-1, we get formula A-2: ； Will Substituting into formula A-2, we get formula A-3: ； 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.

3. The method for establishing an incremental model of controllable power resources according to claim 2, characterized in that, Calculate PE （i） With te （i） The specific steps for obtaining the analytical expression are as follows: Rearranging formula A-3, we get formula B-1: ； Introduce an auxiliary parameter o, let This leads to formula B-2: ； Where Aa, Bb, and Cc represent new parameters of p: ； By introducing the parameter k a second time, we can rewrite formula B-2 and obtain formula B-3: ； Define formula B-4: ； Calculate the smallest real root km of formula B-4; Based on km, we obtain formula B-5: ； Based on formula B-5, the solution for o is obtained: ； Based on the solution for o, we obtain pe. （i） The analytical expression: ； According to PE （i） The analytical expression is obtained from te. （i） Analytical expression: ； Functions based on PE （i） and te （i） The analytical expression is used to calculate the power and discharge time of the first to the acth energy storage devices.

4. The method for establishing an incremental model of controllable power resources according to claim 2, characterized in that, The steps of adjusting the power and discharge time of each energy storage device based on the discharge characteristic function also include: like < Then, obtain the reactive power Qu at the power supply end of the target area, and adjust the transmission power and reactive power at the power supply end; Calculate the additional power ΔPu at the power supply end: ； Calculate the power factor σ of the target region, σ = Pu / Qu; Assuming σ remains constant, the new reactive power Qu is the result of an increase of ΔPu in transmission power: ; Assuming Qu remains constant, the new power factor σu is the result of an increase in transmission power of ΔPu: ; Calculate the reactive power variation coefficient bQ: ; The power factor variation coefficient bσ: ; Compare the magnitudes of bQ and bσ, and adjust the reactive power and voltage phase angle at the power supply end. If bQ ≥ bσ, then the voltage phase angle at the power supply terminal will be increased by (1 + bσ) times; If bQ < bσ, then the reactive power at the power supply end will be adjusted to Quu.

5. The method for establishing an incremental model of controllable power resources according to claim 1, characterized in that, The specific steps for constructing the admittance matrix are as follows: The number of electricity-consuming areas (u_a) within the first municipal district is counted, the topological relationship between the first municipal district and the electricity-consuming areas is obtained, and the initial admittance matrix of the first municipal district is constructed. The expected power consumption and reactive power of the first to u_a electricity-consuming areas within the first municipal district are obtained, and the complex power Se of the first to u_a electricity-consuming areas within the first municipal district is obtained. （1） ~Se （ual） ; Obtain the supply voltage Ug of the first to the ualth power consumption areas. （1） ~Ug （ual） ; 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 connected, then the equivalent admittance Y of the first power consumption area with respect to the second power consumption area... （1,2） : ; If they are not connected, then the equivalent admittance Y of the first power consumption area with respect to the second power consumption area... （1,2） =0; Calculate the equivalent admittance Y of the first power consumption zone with respect to the third to the 1st power consumption zones. （1,3） ~Y （1,ual） ; Calculate Y （1,3） ~Y （1,ual） The sum of these values, as the self-inductance Y of the first power consumption area. （1,1） .

6. The method for establishing an incremental model of controllable power resources according to claim 5, characterized in that, The steps for constructing the admittance matrix also include: Let the complex power of the m-th power consumption area be Se. （m） The power supply voltage is Ug （m） The complex power of the nth power consumption area is Se. （n） The power supply voltage is Ug （n） ; Let Y be the equivalent admittance of the m-th power consumption area with respect to the n-th power consumption area. （m,n） Define the calculation formula C-1-1: ； Let Y be the equivalent admittance of the m-th power consumption area with respect to the l-th power consumption area. （m,l） Let the self-admittance of the m-th power consumption area be Y. （m,m） Define the calculation formula C-1-2: ； Calculate the equivalent admittance and self-admittance Y of the second to the 1ath power consumption areas. （2,l） ~Y （ual,ual） ; Construct the initial admittance matrix YY （1） : ； Extract Y （1,l） ~Y （ual,ual） The real part is used to obtain the equivalent conductance G from the first to the 1st power consumption area. （1,l） ~G （ual,ual） Extract the imaginary part to obtain the equivalent susceptance B of the first to the 1st power consumption areas. （1,l） ~B （ual,ual） ; Construct the initial admittance matrices for the 2nd to mnth municipal districts; based on matrix YY （1） The calculation equations for the initial power flow calculation of the first municipal district are constructed, and the initial phase of the power supply voltage of each power consumption area in the first municipal district is calculated.

7. The method for establishing an incremental model of controllable power resources according to claim 6, characterized in that, The specific steps for calculating the initial phase of the supply voltage are as follows: Let the power supply voltage of the w-th power consumption area be Ug. （w） The voltage phase is θ （w） The expected power consumption is Pz （w） The reactive power is Qz （w） The complex power of the w-th power consumption area is Se. （w） Se （w） =Pz （w） +Qz （w） ×j; The power supply voltage for the vth power consumption area is Ug （v） The voltage phase is θ （v） ; The equivalent admittance of the w-th power consumption area with respect to the v-th 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） ; Let the conjugate of the power supply voltage of the w-th power consumption area be (Ug （w） ) Define formula C-2-1: ; The conjugate of the power supply voltage for the vth power consumption area is (Ug （v） ) , ; Define formula C-2-2: ； Expanding formula C-2-2, we get formula C-2-3: ； Where, θ （w，v） This represents the phase difference between the power supply voltage of the w-th power consumption area and the v-th power consumption area; Construct matrix I （1） Sum matrix U （1） ; Define formula C-3: ； Let the initial phase of the power supply voltage for the first to the 1st power consumption area be θz. （1） ~θz （ual） ; Rearranging formula C-3, we get formula D. （1） Formula D （ual） : Formula D （1） for: ； Similarly, formula D （ual） for: ； Calculate θz using Newton-Raphson's algorithm. （1） ~θz （ual） The value of is used to obtain the initial phase of the power supply voltage for the first to the ualth power consumption areas; Calculate the initial phase of the application power zone for the 2nd to mnth municipal districts; obtain the actual power consumption of all power zones in the 1st district, and adjust the power supply of the 1st to ualth power zones.

8. The method for establishing an incremental model of controllable power resources according to claim 7, characterized in that, The specific steps for adjusting the power supply of the first to the 11th power consumption areas are as follows: Obtain the expected power consumption po of the first to the 1st power consumption areas. （1） ~po （ual） Actual power consumption (Pa) （1） ~pa （ual） ; Determine whether the expected power consumption is equal to the actual power consumption, and adjust the power supply voltage and phase of each municipal district to the corresponding power consumption area based on the relationship between the actual power consumption and the expected power consumption. In the first to the 11th power consumption area, if the expected power consumption of only one power consumption area is not equal to the actual power consumption, then the power consumption area where the expected power consumption is not equal to the actual power consumption is designated as an abnormal area. Obtain the expected power consumption pox, actual power consumption pax, actual supply voltage ux, and actual supply voltage phase θx in the abnormal region; Obtain the expected supply voltage uy in the abnormal region and the initial supply voltage phase θy; Determine whether ux and uy, as well as θx and θy, are equal, and adjust the power supply voltage and phase in abnormal areas or other areas accordingly. If ux and uy are equal, it indicates that the reactive power in the abnormal region has changed; Compare the values ​​of POX and PAX, and adjust the power supply voltage for other areas; If pox > pax, then reduce the power supply voltage in other areas. times; If pox < pax, then increase the power supply voltage in other areas. times; If θx and θy are equal, it indicates that the active power in the abnormal region has changed; Compare the values ​​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 zones. times; If pox < pax, then reduce the phase of the power supply voltage in other zones. times; If ux and uy, as well as θx and θy, are all unequal, then the outlier region is used as the balancing node to reconstruct matrix YY. （1） Adjust the actual power consumption and actual reactive power of all 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 actual power consumption and actual reactive power of all power consumption areas are as follows: Let the abnormal zone be numbered r; recalculate the equivalent admittance rY of the abnormal zone with respect to the 1st to the 1st power consumption zone. （r，1） ~rY （r，ual） Self-guided nanorY （r，r） ; Calculate the change in self-admittance ΔY in the abnormal region. （r，r） The change in equivalent admittance ΔY （r，1） ; Calculate the changes in equivalent admittance and self-admittance for the first to the 11th power plants to obtain ΔY. （1，1） ~ΔY （ual，ual） Resulting in matrix YY （2） ; Calculate the voltage change Δu at the slack node. Phase change Δθ ; Calculate the change in actual power supply voltage Δu for the first to the 1st power station. （1） ~Δu （ual） The actual change in phase of the supply voltage Δθ （1） ~Δθ （ual） ; Let ΔP be the change in actual power consumption of the first to the 1st power station. （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, the Newton-Raphson algorithm is reused to calculate ΔP. （1） ~ΔP （ual） and ΔQ （1） ~ΔQ （ual） The value; The actual power consumption of the first to the 11th power stations will be increased by ΔP sequentially. （1） ~ΔP （ual） The actual reactive power increases by ΔQ sequentially. （1） ~ΔQ （ual） ; If, in the first to the 11th power consumption areas, the expected power consumption of multiple power consumption areas is not equal to the actual power consumption, then the steps of adjusting the actual power consumption and actual reactive power of all power consumption areas are repeated.

Images