A method and system for coordinated distributed control of multiple distribution areas in a distribution network

Abstract

The present invention provides a method and system for coordinated distributed control of multiple distribution zones in a distribution network. The method includes constructing a power flow equation to calculate a first fitness value of a distribution point in the distribution network; calculating a second fitness value, a stationary scale, and an intra-zone matching scale between a sub-distribution point and the remaining sub-distribution points, and calculating a comprehensive fitness value; calculating a fitness difference, and arranging the fitness differences between the remaining sub-distribution points and the initial sub-distribution point from large to small, and combining them from large to small to obtain a pending merged partition; sequentially absorbing the remaining sub-distribution points into the pending merged partition until a preset condition is met; calculating carbon assets and equipment investment prices based on the pending merged partition, optimizing the pending merged partition to obtain a new partition; calculating the supply index of the distribution network in each time period based on the new partition to calculate a sharing index; obtaining the official purchase and sale cost of electricity and the power exchange cost, and obtaining an optimal planning scheme. The present invention can make the partitioning more reasonable and reduce subsequent operating costs.

Description

Technical Field

[0001] The present invention relates to the technical field of distribution networks, and in particular to a method and system for coordinated distributed control of multiple distribution areas in a distribution network. Background Art

[0002] Due to the development of new energy vehicles and electrical appliances, a large number of new access nodes will be added to large-scale distribution networks, which will lead to imbalance in the distribution system. In order to reduce the impact of this imbalance on the safety and life of the distribution system, it is necessary to distribute and merge nodes with high electrical coupling and strong correlation.

[0003] In the existing technology, the research on the correlation between the sub-distribution points in the distribution network is generally based on historical data to predict future operating parameters. Generally, only the degree of influence between the sub-distribution points in the distribution network is used as a reference standard. The sample data is small and the prediction algorithm is not stable enough. The impact of the constraints on the nodes is not taken into account. At the same time, the shared coordination of the distribution network after partitioning is not fully utilized to optimize the implementation plan. Summary of the Invention

[0004] In view of the deficiencies in the prior art, the present invention aims to provide a method for coordinated distributed control of multiple substations in a distribution network, aiming to solve the technical problems mentioned in the background technology.

[0005] In order to achieve the above object, the present invention is implemented through the following technical solutions:

[0006] A method for coordinated distributed control of multiple distribution areas in a distribution network comprises the following steps:

[0007] Constructing a power flow equation of the distribution network, and calculating a first fitness value of each sub-distribution point in the distribution network based on the power flow equation;

[0008] traversing and calculating a second fitness value, a stable scale, and an intra-region matching scale between a sub-distribution point and the remaining sub-distribution points, and calculating a comprehensive fitness value by integrating the second fitness value, the stable scale, and the intra-region matching scale;

[0009] Obtaining a fitness difference according to the comprehensive fitness value and the first fitness value, and arranging the fitness differences between the remaining sub-distribution points and the initially selected sub-distribution point from largest to smallest, selecting the sub-distribution point corresponding to the largest fitness difference and combining it with the initial sub-distribution point to obtain a pending merged partition;

[0010] The pending merged partition is used as a temporary sub-distribution point, and the remaining sub-distribution points are absorbed into the pending merged partition in sequence until the preset conditions are met;

[0011] Calculate the pending merge partition results of the remaining sub-distribution points in sequence;

[0012] Based on the partitioned pending merged partitions, calculating carbon assets and equipment investment prices, optimizing the pending merged partitions with the goal of maximizing the carbon asset acquisition rate to obtain new partitions;

[0013] Based on the newly divided partitions, the power demand is obtained and the supply index of each time period in the distribution network is calculated to calculate the sharing index;

[0014] Obtain the official electricity purchase and sales costs and electricity exchange costs, and obtain the optimal planning solution with the goal of minimizing operating costs.

[0015] According to one aspect of the above technical solution, constructing the power flow equation of the distribution network and calculating the first fitness value of each sub-distribution point in the distribution network based on the power flow equation specifically includes:

[0016] Based on the power flow parameters of the distribution network, the voltage amplitude changes when all the sub-distribution points are connected to active equipment and reactive equipment are calculated respectively. and ;

[0017] ;

[0018] in, represents the change in the voltage phase angle in the sub-distribution point, represents the change in the voltage amplitude value in the sub-distribution point, represents the change in active power in the sub-distribution point, represents the change in reactive power in the sub-distribution point, Indicates the change in voltage phase angle when the sub-distribution point is connected to active equipment, Indicates the change in voltage phase angle when the sub-distribution point is connected to reactive equipment;

[0019] According to the change of the voltage amplitude value when the equipment in the sub-distribution point is active and reactive and , calculate the distance weight of the sub-distribution point;

[0020] ;

[0021] ;

[0022] in, It represents the influence index of the voltage amplitude value of sub-distribution point i when other sub-distribution points j connected to sub-distribution point i are connected to reactive power. It indicates the influence index of the voltage amplitude value of other connected sub-distribution points j when the sub-distribution point i is connected to reactive power. It represents the influence index of the voltage amplitude value of the sub-distribution point i when other sub-distribution points j connected to the sub-distribution point i are connected to the active power, It indicates the influence index of the voltage amplitude value of other connected sub-distribution points j when the sub-distribution point i is connected to the active power. represents the reactive distance weight between two sub-distribution points, Represents the active distance weight between two sub-distribution points;

[0023] The reactive power influence factor ω and the active power influence factor (1-ω) are set, and the active power distance weight and the reactive power distance weight are integrated to obtain the overall weight;

[0024] ;

[0025] Normalizing the overall weight to obtain the edge weight of the sub-distribution point;

[0026] Calculating a first fitness value of the sub-distribution point according to the edge weight of the sub-distribution point;

[0027] ;

[0028] ;

[0029] in, and Represents the total weight of sub-distribution point i and sub-distribution point j, m V Represents the total edge weight of the sub-distribution points. If the two sub-distribution points are connected, then ρ(i, j) = 1, otherwise ρ(i, j) = 0. represents the edge weight.

[0030] According to one aspect of the above technical solution, the traversal calculation of the second fitness value, the stationary scale, and the intra-region matching scale between the sub-distribution point and the remaining sub-distribution points, and the calculation of the comprehensive fitness value by integrating the second fitness value, the stationary scale, and the intra-region matching scale, specifically includes:

[0031] Select a sub-distribution point and traverse to calculate the second fitness value between the sub-distribution point i and the remaining sub-distribution points q ;

[0032] ;

[0033] Obtain all the equipment in the sub-distribution point i and the remaining sub-distribution points q to obtain the reactive power supply N Q ;

[0034] Get the voltage over-limit value and get the minimum reactive power supply N W ;

[0035] Based on the reactive power supply amount N Q , Minimum reactive power supply N W , and reactive distance weight , calculate the stationary scale between sub-distribution point i and sub-distribution point q;

[0036] ;

[0037] Calculating a first intra-region matching index and a second intra-region matching index based on the edge weights and the current partitioning situation;

[0038] ;

[0039] ;

[0040] Among them, C represents the serial number of the partition, G represents the total number of partitions, G C Indicates C partition, N C Indicates the total number of neutron distribution points in the C partition, B indicates the total number of neutron distribution points in this calculation, represents the reactive distance weight between sub-distribution point i and sub-distribution point q, represents the edge weight between sub-distribution point i and sub-distribution point q;

[0041] Calculating an intra-region matching scale based on the first intra-region matching index and the second intra-region matching index;

[0042] ;

[0043] Calculating a comprehensive fitness value by combining the second fitness value, the stability scale, and the intra-region matching scale;

[0044] .

[0045] According to one aspect of the above technical solution, the fitness difference is obtained based on the comprehensive fitness value and the first fitness value, and the fitness differences between the remaining sub-distribution points and the initially selected sub-distribution point are arranged from large to small, and the sub-distribution point corresponding to the largest fitness difference is selected and combined with the initial sub-distribution point to obtain a pending merged partition, which specifically includes:

[0046] Subtract the comprehensive fitness value between the sub-distribution point i and the remaining sub-distribution points q from the first fitness value of the sub-distribution point i to obtain a fitness difference;

[0047] Arrange the fitness differences between the remaining sub-distribution points q and the initially selected sub-distribution point i from large to small;

[0048] The sub-distribution point q corresponding to the largest fitness difference is selected and combined with the initial sub-distribution point i to obtain a pending merged partition.

[0049] According to one aspect of the above technical solution, the method of using the pending merged partition as a temporary sub-distribution point and sequentially absorbing the remaining sub-distribution points into the pending merged partition until a preset condition is met specifically includes:

[0050] The pending merged partition is used as a temporary sub-distribution point iq;

[0051] Subtract the comprehensive fitness value between the temporary sub-distribution point iq and the remaining sub-distribution points q that have not been merged from the comprehensive fitness value of the temporary sub-distribution point iq to obtain a fitness difference;

[0052] All other sub-distribution points q that have not been merged are absorbed into the pending merge partition in turn until the fitness difference becomes negative.

[0053] According to one aspect of the above technical solution, the calculation of carbon assets and equipment investment prices based on the partitioned pending merged partitions, and the optimization of the pending merged partitions with the goal of maximizing the carbon asset acquisition rate to obtain new partitions specifically includes:

[0054] Calculating carbon assets based on market factors and planning schemes, and calculating equipment investment prices for the to-be-determined merged subareas based on the status of access equipment at sub-distribution points in the to-be-determined merged subareas;

[0055] ;

[0056] ;

[0057] Among them, C Y represents carbon assets, D represents the number of typical load calculation days covered in the year, and W T represents the instantaneous price of carbon market, Tθ represents the number of time periods covering the typical annual load calculation day, M is the number of all sub-distribution points in the distribution network, and W a,k,t represents the output power of the kth sub-distribution point in time period t, γ1 represents the carbon dioxide emission index, C J represents the equipment investment price, G inv Represents the annual investment cost of the distribution network, G o represents the annual electricity transaction fee, G l Represents the annual network loss cost, G b Indicates annual maintenance fee;

[0058] Based on the carbon assets and the equipment investment price, with the goal of maximizing the carbon asset acquisition rate, a first objective function is constructed, and a first constraint condition is set according to demand;

[0059] ;

[0060] Among them, P C represents the carbon asset acquisition rate;

[0061] Solving the first objective function by a machine algorithm to obtain a first optimal solution and a first optimal planning scheme for the first objective function;

[0062] Preset a first optimal solution threshold. If the first optimal solution is smaller than the first optimal solution threshold, remove the sub-distribution points in the reverse order of the construction of the pending merged partitions, and cyclically calculate the first optimal solution based on the updated pending merged partitions until the first optimal solution is larger than the first optimal solution threshold.

[0063] When the first optimal solutions of the first objective function solved based on the partitions to be determined and merged converge, the final partition to be determined and merged is determined as a new partition.

[0064] According to one aspect of the above technical solution, based on the divided new partition, obtaining the power demand and calculating the supply index of each time period in the distribution network to calculate the sharing index specifically includes:

[0065] Based on the device access status of each sub-distribution point in each new partition, obtaining the power demand in each new partition;

[0066] Calculating the supply index of each time period in the distribution network based on the power demand in each of the new zones;

[0067] ;

[0068] in, represents the electricity sales in the Nth new partition during period t, represents the amount of electricity purchased in the Nth new partition during period t, and tg(t) represents the supply index during period t;

[0069] Calculate the shared median price based on the current market carbon emission factors;

[0070] ;

[0071] in, Indicates the electricity purchase price of the upper power grid, Indicates the electricity price of the upper power grid, G P Represents the carbon emission index, G Pmax represents the maximum value of carbon emission index, ε is a constant, represents the shared median price;

[0072] The shared median price is optimized, and the shared electricity purchase price and the shared electricity selling price are calculated based on the optimized shared median price.

[0073] ;

[0074] ;

[0075] ;

[0076] in, represents the shared correction, tc=1 / tg, Indicates the shared electricity purchase price after optimization, Indicates the shared electricity sales price after optimization.

[0077] According to one aspect of the above technical solution, obtaining the official electricity purchase and sales costs and the electricity exchange costs, with the goal of minimizing operating costs, to obtain the optimal planning solution specifically includes:

[0078] Obtain official electricity purchase and sales costs and electricity exchange costs;

[0079] ;

[0080] ;

[0081] Among them, T Y represents the official cost of electricity purchase and sale, Indicates the amount of electricity purchased by the upper power grid. Indicates the power sales of the upper power grid, D c Indicates the time of day, T η represents the electricity exchange cost, η g represents the energy storage penalty coefficient, Indicates the energy storage charging power, Represents the energy storage discharge power, η k represents the first adjustable electric load penalty coefficient, represents the first adjustable electric load power, η l represents the second adjustable electric load penalty coefficient, Indicates the second adjustable electric load power;

[0082] With the goal of minimizing operating costs, construct the second objective function and set the second constraint conditions according to requirements;

[0083] ;

[0084] Where T represents operating cost;

[0085] The second objective function after considering the shared electricity purchase price and the shared electricity sales price is simulated and run using a power system model to obtain a second optimal solution and a second optimal planning scheme for the second objective function.

[0086] The present invention also provides a distribution network multi-zone distributed control coordination system, comprising:

[0087] A first calculation module is used to construct a power flow equation of the distribution network, and calculate a first fitness value of each sub-distribution point in the distribution network based on the power flow equation;

[0088] A second calculation module is configured to traverse and calculate a second fitness value, a stationary scale, and an intra-region matching scale between a sub-distribution point and the remaining sub-distribution points, and calculate a comprehensive fitness value by integrating the second fitness value, the stationary scale, and the intra-region matching scale;

[0089] A first merging module is configured to obtain a fitness difference according to the comprehensive fitness value and the first fitness value, and arrange the fitness differences between the remaining sub-distribution points and the initially selected sub-distribution point from large to small, and select the sub-distribution point corresponding to the largest fitness difference and combine it with the initial sub-distribution point to obtain a pending merged partition;

[0090] The second merging module is configured to use the pending merged partition as a temporary sub-distribution point and sequentially absorb the remaining sub-distribution points into the pending merged partition until a preset condition is met;

[0091] A third merging module is used to sequentially calculate the pending merging partition results of the remaining sub-distribution points;

[0092] Optimization module: used for calculating carbon assets and equipment investment prices based on the partitioned pending merged partitions, optimizing the pending merged partitions with the goal of maximizing the carbon asset acquisition rate to obtain new partitions;

[0093] A third calculation module is configured to obtain the power demand based on the newly divided partitions and calculate the supply index of each time period in the distribution network to calculate the sharing index;

[0094] The fourth calculation module is used to obtain the official electricity purchase and sales costs and electricity exchange costs, with the goal of minimizing operating costs to obtain the optimal planning solution.

[0095] Compared with the prior art, the present invention has the following beneficial effects:

[0096] By calculating the first fitness value of each sub-distribution point, and traversing and calculating the second fitness value, stable scale, and intra-region matching scale between one of the sub-distribution points and the remaining sub-distribution points, the comprehensive fitness value can be obtained. By calculating the comprehensive fitness value of the two sub-distribution points and the first fitness value of the initial sub-distribution point itself, the fitness difference can be obtained. The fitness relationship between the initial sub-distribution point and the remaining sub-distribution points is determined according to the fitness difference, and then the merging order is determined. Each time the initial sub-distribution point merges with a remaining sub-distribution point, it is marked as a pending merge partition, and the pending merge partition is used as a temporary sub-distribution point, and the remaining The unmerged sub-distribution points are selected until the preset conditions are met, and then another initial sub-distribution point is selected and the above steps are followed until all pending merged partitions are obtained. After this step, the initial merger of the sub-distribution points is completed; then, by calculating the carbon assets and equipment investment prices of the pending merged partitions, if they exceed the preset targets, the sub-distribution points in the pending merged partitions are subtracted in turn, and the calculation is repeated until the requirements are met to obtain a new partition; then, the sharing indicators are calculated based on the operating conditions of the sub-distribution points in the divided new partitions, and finally, with the lowest operating cost as the goal, the optimal planning scheme between the sub-distribution points in the new partitions is calculated.

[0097] The present invention fully considers the stationary scale and intra-area matching scale between sub-distribution points when zoning, and then obtains a comprehensive fitness value, which can make the zoning between sub-distribution points more reasonable, with higher matching and stronger applicability. Secondly, by calculating the carbon assets and equipment investment prices of the pending merged partitions, a better balance can be achieved between partition matching and carbon emission costs, avoiding excessive optimization of one aspect while ignoring the other factor. The shared indicators are optimally incorporated to calculate the optimal planning scheme between each sub-distribution point in the new partition, thereby realizing coordinated control of the distribution network system. BRIEF DESCRIPTION OF THE DRAWINGS

[0098] Figure 1 Flowchart of the method for coordinated distributed control of multiple distribution areas in a distribution network according to the first embodiment of the present invention;

[0099] Figure 2 This is a structural block diagram of a multi-zone distributed control coordination system for a distribution network according to a second embodiment of the present invention;

[0100] The following specific embodiments will further illustrate the present invention in conjunction with the above-mentioned drawings. DETAILED DESCRIPTION

[0101] To facilitate understanding of the present invention, the present invention will be described more fully below with reference to the accompanying drawings. The drawings illustrate several embodiments of the present invention. However, the present invention may be implemented in many different forms and is not limited to the embodiments described herein. Rather, these embodiments are provided to provide a more thorough and comprehensive understanding of the present invention.

[0102] It should be noted that when an element is referred to as being "fixed to" another element, it may be directly on the other element or there may be an intermediate element. When an element is referred to as being "connected to" another element, it may be directly connected to the other element or there may be an intermediate element. The terms "vertical," "horizontal," "left," "right," and similar expressions used herein are for illustrative purposes only.

[0103] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one skilled in the art to which this invention pertains. The terms used in this specification of the present invention are for the purpose of describing specific embodiments only and are not intended to limit the present invention. The term "and / or" as used herein includes any and all combinations of one or more of the associated listed items.

[0104] See also Figure 1 , which shows a method for coordinated distributed control of multiple distribution areas in a distribution network according to a first embodiment of the present invention, comprising the following steps:

[0105] S10, constructing a power flow equation for the distribution network, and calculating a first fitness value of each sub-distribution point in the distribution network based on the power flow equation;

[0106] S20, traversing and calculating a second fitness value, a stable scale, and an intra-region matching scale between the sub-distribution point and the remaining sub-distribution points, and calculating a comprehensive fitness value by combining the second fitness value, the stable scale, and the intra-region matching scale;

[0107] S30, obtaining a fitness difference based on the comprehensive fitness value and the first fitness value, and arranging the fitness differences between the remaining sub-distribution points and the initially selected sub-distribution point from largest to smallest, selecting the sub-distribution point corresponding to the largest fitness difference and combining it with the initial sub-distribution point to obtain a pending merged partition;

[0108] S40, using the pending merged partition as a temporary sub-distribution point, and sequentially absorbing the remaining sub-distribution points into the pending merged partition until a preset condition is met;

[0109] S50, sequentially calculating pending merge partition results of the remaining sub-distribution points;

[0110] S60, based on the partitions to be merged, calculating carbon assets and equipment investment prices, optimizing the partitions to be merged with the goal of maximizing the carbon asset acquisition rate to obtain new partitions;

[0111] S70, based on the newly divided partition, obtaining power demand and calculating supply indicators for each time period in the distribution network to calculate a sharing indicator;

[0112] S80, obtains the official electricity purchase and sales costs and electricity exchange costs, and obtains the optimal planning scheme with the goal of minimizing operating costs.

[0113] It can be understood that the present invention calculates the first fitness value of each sub-distribution point and traverses and calculates the second fitness value, the stability scale, and the intra-region matching scale between one of the sub-distribution points and the remaining sub-distribution points to obtain a comprehensive fitness value. By calculating the comprehensive fitness value of the two sub-distribution points and the first fitness value of the initial sub-distribution point itself, a fitness difference can be obtained. The fitness relationship between the initial sub-distribution point and the remaining sub-distribution points is determined based on the fitness difference, and then the merging order is determined. Each time the initial sub-distribution point merges with a remaining sub-distribution point, it is marked as a pending merge partition, and the pending merge partition is used as a temporary sub-distribution point. Merge the remaining unmerged sub-distribution points in sequence until the preset conditions are met, and then select another initial sub-distribution point and follow the above steps until all the pending merged partitions are obtained. After this step, the initial merger of the sub-distribution points is completed; then calculate the carbon assets and equipment investment prices of the pending merged partitions. If they exceed the preset targets, subtract the sub-distribution points in the pending merged partitions in sequence, and repeat the calculation until the requirements are met to obtain a new partition; then calculate the sharing indicators based on the operating conditions of the sub-distribution points in the divided new partitions, and finally calculate the optimal planning scheme between the sub-distribution points in the new partition with the lowest operating cost as the goal.

[0114] The present invention fully considers the stationary scale and intra-area matching scale between sub-distribution points when zoning, and then obtains a comprehensive fitness value, which can make the zoning between sub-distribution points more reasonable, with higher matching and stronger applicability. Secondly, by calculating the carbon assets and equipment investment prices of the pending merged partitions, a better balance can be achieved between partition matching and carbon emission costs, avoiding excessive optimization of one aspect while ignoring the other aspect. Finally, shared indicators are incorporated to calculate the optimal planning scheme between each sub-distribution point in the new partition, thereby realizing coordinated control of the distribution network system.

[0115] Furthermore, the step S10 specifically includes:

[0116] Based on the power flow parameters of the distribution network, the voltage amplitude changes when all the sub-distribution points are connected to active equipment and reactive equipment are calculated respectively. and ;

[0117] ;

[0118] in, represents the change in the voltage phase angle in the sub-distribution point, represents the change in the voltage amplitude value in the sub-distribution point, represents the change in active power in the sub-distribution point, represents the change in reactive power in the sub-distribution point, Indicates the change in voltage phase angle when the sub-distribution point is connected to active equipment, Indicates the change in voltage phase angle when the sub-distribution point is connected to reactive equipment;

[0119] According to the change of the voltage amplitude value when the equipment in the sub-distribution point is active and reactive and , calculate the distance weight of the sub-distribution point;

[0120] ;

[0121] ;

[0122] in, It represents the influence index of the voltage amplitude value of sub-distribution point i when other sub-distribution points j connected to sub-distribution point i are connected to reactive power. It indicates the influence index of the voltage amplitude value of other connected sub-distribution points j when the sub-distribution point i is connected to reactive power. It represents the influence index of the voltage amplitude value of the sub-distribution point i when other sub-distribution points j connected to the sub-distribution point i are connected to the active power, It indicates the influence index of the voltage amplitude value of other connected sub-distribution points j when the sub-distribution point i is connected to the active power. represents the reactive distance weight between two sub-distribution points, Represents the active distance weight between two sub-distribution points;

[0123] The reactive power influence factor ω and the active power influence factor (1-ω) are set, and the active power distance weight and the reactive power distance weight are integrated to obtain the overall weight;

[0124] ;

[0125] Normalizing the overall weight to obtain the edge weight of the sub-distribution point;

[0126] Calculating a first fitness value of the sub-distribution point according to the edge weight of the sub-distribution point;

[0127] ;

[0128] ;

[0129] in, and Represents the total weight of sub-distribution point i and sub-distribution point j, m V Represents the total edge weight of the sub-distribution points. If the two sub-distribution points are connected, then ρ(i, j) = 1, otherwise ρ(i, j) = 0. represents the edge weight.

[0130] It can be understood that the first fitness value of the sub-distribution point itself can be calculated by the above method. The main idea is to select a sub-distribution point and calculate the voltage amplitude value influence index of the remaining adjacent sub-distribution points, and then obtain the distance weight and the first fitness value. The higher the first fitness value, the higher the degree of matching between the sub-distribution point and the remaining adjacent sub-distribution points.

[0131] Furthermore, the step S20 specifically includes:

[0132] Select a sub-distribution point and traverse to calculate the second fitness value between the sub-distribution point i and the remaining sub-distribution points q ;

[0133] ;

[0134] Obtain all the equipment in the sub-distribution point i and the remaining sub-distribution points q to obtain the reactive power supply N Q ;

[0135] Get the voltage over-limit value and get the minimum reactive power supply N W ;

[0136] Based on the reactive power supply amount N Q , Minimum reactive power supply N W , and reactive distance weight , calculate the stationary scale between sub-distribution point i and sub-distribution point q;

[0137] ;

[0138] Calculating a first intra-region matching index and a second intra-region matching index based on the edge weights and the current partitioning situation;

[0139] ;

[0140] ;

[0141] Among them, C represents the serial number of the partition, G represents the total number of partitions, G C Indicates C partition, N C Indicates the total number of neutron distribution points in the C partition, B indicates the total number of neutron distribution points in this calculation, represents the reactive distance weight between sub-distribution point i and sub-distribution point q, represents the edge weight between sub-distribution point i and sub-distribution point q;

[0142] Calculating an intra-region matching scale based on the first intra-region matching index and the second intra-region matching index;

[0143] ;

[0144] Calculating a comprehensive fitness value by combining the second fitness value, the stability scale, and the intra-region matching scale;

[0145] .

[0146] It can be understood that the calculation method of the second fitness value is consistent with the calculation method of the first fitness value, except that when calculating the first fitness value, it is necessary to calculate the distance weight between the initially selected sub-distribution point and the remaining adjacent sub-distribution points, while when calculating the second fitness value, it is necessary to calculate the distance weight between the initially selected sub-distribution point and all the remaining sub-distribution points in the distribution network. Therefore, this step omits some steps and omits some reference to the content of calculating the first fitness value. The second fitness value can reflect the matching degree between the initially selected sub-distribution point and the remaining sub-distribution points. The higher the second fitness value, the higher the matching degree. Then, according to the reactive power supply amount N Q , Minimum reactive power supply N W The stability scale between the two sub-distribution points is calculated. The stability scale reflects the reactive balancing ability between the two sub-distribution points. Then, based on the edge weight and the current partitioning situation, the intra-region matching scale is calculated. The intra-region matching scale reflects the degree of coordination between the sub-distribution points in this partitioning situation, thereby improving the accuracy of the partitioning. The above two indicators are two more important indicators for judging the relationship between the two sub-distribution points, avoiding the poor partitioning effect caused by only considering the voltage amplitude value between the two sub-distribution points.

[0147] Furthermore, the step S30 specifically includes:

[0148] Subtract the comprehensive fitness value between the sub-distribution point i and the remaining sub-distribution points q from the first fitness value of the sub-distribution point i to obtain a fitness difference;

[0149] Arrange the fitness differences between the remaining sub-distribution points q and the initially selected sub-distribution point i from large to small;

[0150] The sub-distribution point q corresponding to the largest fitness difference is selected and combined with the initial sub-distribution point i to obtain a pending merged partition.

[0151] It can be understood that the fitness difference between the initially selected sub-distribution point and the remaining sub-distribution points can be calculated by taking the difference. The fitness difference reflects "how much the second fitness value will change from the first fitness value if the initial sub-distribution point is combined with the remaining sub-distribution point". If the fitness difference is positive, and the larger the value is, the better the effect of the merger is. Therefore, the order of merging sub-distribution points is determined by the fitness difference.

[0152] Furthermore, the step S40 specifically includes:

[0153] The pending merged partition is used as a temporary sub-distribution point iq;

[0154] Subtract the comprehensive fitness value between the temporary sub-distribution point iq and the remaining sub-distribution points q that have not been merged from the comprehensive fitness value of the temporary sub-distribution point iq to obtain a fitness difference;

[0155] All other sub-distribution points q that have not been merged are absorbed into the pending merge partition in turn until the fitness difference becomes negative.

[0156] It can be understood that when the pending merge partition absorbs sub-distribution points in turn, the updated pending merge partition is regarded as a new sub-distribution point i each time, and is merged with the remaining unmerged sub-distribution points q according to steps S20 and S30. Therefore, when matching the scale in the calculation area in step S30, G and B will continue to change according to this calculation. When the merger ends depends on when the fitness difference is negative. If the fitness difference is negative, it means that the merger effect is very poor, or even less than the fitness value of the sub-distribution point itself.

[0157] Furthermore, the step S60 specifically includes:

[0158] Calculating carbon assets based on market factors and planning schemes, and calculating equipment investment prices for the to-be-determined merged subareas based on the status of access equipment at sub-distribution points in the to-be-determined merged subareas;

[0159] ;

[0160] ;

[0161] Among them, C Y represents carbon assets, D represents the number of typical load calculation days covered in the year, and W Trepresents the instantaneous price of carbon market, Tθ represents the number of time periods covering the typical annual load calculation day, M is the number of all sub-distribution points in the distribution network, and W a,k,t represents the output power of the kth sub-distribution point in time period t, γ1 represents the carbon dioxide emission index, C J represents the equipment investment price, G inv Represents the annual investment cost of the distribution network, G o represents the annual electricity transaction fee, G l Represents the annual network loss cost, G b Indicates annual maintenance fee;

[0162] Specifically,

[0163] ;

[0164] ;

[0165] ;

[0166] ;

[0167] Where q represents the discount rate, R a represents the planned number of years for the equipment in the sub-distribution point, e1 represents the investment cost of the equipment in the sub-distribution point, and e oa Indicates the maintenance fee for the equipment in the sub-distribution point, e p represents the average electricity purchase price of the transmission network, P s,k,t represents the output power of the thermal power grid at the ith sub-distribution point during period t, e l Indicates the network loss cost coefficient, L kc Indicates the current of the equipment k and c lines in the sub-distribution point, C kc Indicates the impedance of the k and c lines of the equipment within the sub-distribution point;

[0168] Based on the carbon assets and the equipment investment price, with the goal of maximizing the carbon asset acquisition rate, a first objective function is constructed, and a first constraint condition is set according to demand;

[0169] ;

[0170] Among them, P C represents the carbon asset acquisition rate;

[0171] Solving the first objective function by a machine algorithm to obtain a first optimal solution and a first optimal planning scheme for the first objective function;

[0172] Preset a first optimal solution threshold. If the first optimal solution is smaller than the first optimal solution threshold, remove the sub-distribution points in the reverse order of the construction of the pending merged partitions, and cyclically calculate the first optimal solution based on the updated pending merged partitions until the first optimal solution is larger than the first optimal solution threshold.

[0173] When the first optimal solutions of the first objective function solved based on the partitions to be determined and merged converge, the final partition to be determined and merged is determined as a new partition.

[0174] It is understandable that after all pending merged partitions are determined, the carbon assets and equipment investment prices of the entire distribution network operation under such partitioning conditions are calculated based on the operating conditions of the equipment at the sub-distribution points in the pending merged partitions. Then, based on the carbon assets and equipment investment prices, with the highest carbon asset acquisition rate as the goal, a first objective function is constructed, and the first constraint condition is set according to the requirements to solve the first objective function. The solution process can be implemented using a machine algorithm, such as the SSA algorithm. It should be noted that the first constraint condition is an industry standard. The industry standard is the distribution network flow constraint, power balance constraint, equipment capacity of the sub-distribution point, etc., which will not be described in detail here. After solving the first objective function to obtain the first optimal solution and the first optimal planning scheme, the first optimal solution threshold can be determined according to the requirements. That is, if the maximum carbon asset acquisition rate obtained is less than the target, it is necessary to abandon some nodes in the pending merged partitions, and abandon them in the reverse order of the merger, and the calculation is repeated until the first optimal solution meets the first optimal solution threshold. The pending merged partition is determined as a new partition, and the abandoned nodes need to be recalculated and assigned to the appropriate pending merged partition according to steps S30 and S40. Then, the first optimal solution of the remaining updated pending merged partitions is solved, and all new partitions are obtained in sequence. As an extension of this embodiment, during the calculation process, the first optimal solution threshold can be adjusted so that all sub-distribution points can reach an optimal solution between fitness and carbon asset acquisition rate, or all abandoned sub-distribution points can be divided into new partitions as much as possible in the remaining pending merged partitions. Machine algorithms (such as particle swarm optimization algorithm) can be used here to achieve this.

[0175] Furthermore, the step S70 specifically includes:

[0176] Based on the device access status of each sub-distribution point in each new partition, obtaining the power demand in each new partition;

[0177] Calculating the supply index of each time period in the distribution network based on the power demand in each of the new zones;

[0178] ;

[0179] in, represents the electricity sales in the Nth new partition during period t, represents the amount of electricity purchased in the Nth new partition during period t, and tg(t) represents the supply index during period t;

[0180] Calculate the shared median price based on the current market carbon emission factors;

[0181] ;

[0182] in, Indicates the electricity purchase price of the upper power grid, Indicates the electricity price of the upper power grid, G P Represents the carbon emission index, G Pmax represents the maximum value of carbon emission index, ε is a constant, represents the shared median price;

[0183] The shared median price is optimized, and the shared electricity purchase price and the shared electricity selling price are calculated based on the optimized shared median price.

[0184] ;

[0185] ;

[0186] ;

[0187] in, represents the shared correction, tc=1 / tg, Indicates the shared electricity purchase price after optimization, Indicates the shared electricity sales price after optimization.

[0188] It can be understood that after the new partition is determined, the electricity demand is determined according to the equipment access status of the sub-distribution points in the new partition, and then the supply index is calculated. The supply index reflects the internal supply and demand situation between the sub-distribution points in the new partition. Then the shared median price, shared electricity purchase price and shared electricity selling price are calculated in turn. The shared electricity purchase price and shared electricity selling price can enable the new partitions in the distribution network to share electricity prices. In actual application, the energy interaction within the new partition is mobilized, load reduction and load shifting are effectively carried out, the energy distribution in the distribution network is optimized, and it plays a guiding role in the coordinated control between the partitions in the distribution network.

[0189] Furthermore, the step S80 specifically includes:

[0190] Obtain official electricity purchase and sales costs and electricity exchange costs;

[0191] ;

[0192] ;

[0193] Among them, T Y represents the official cost of electricity purchase and sale, Indicates the amount of electricity purchased by the upper power grid. Indicates the power sales of the upper power grid, D c Indicates the time of day, T η represents the electricity exchange cost, η g represents the energy storage penalty coefficient, Indicates the energy storage charging power, Represents the energy storage discharge power, η k represents the first adjustable electric load penalty coefficient, represents the first adjustable electric load power, η l represents the second adjustable electric load penalty coefficient, Indicates the second adjustable electric load power;

[0194] With the goal of minimizing operating costs, construct the second objective function and set the second constraint conditions according to requirements;

[0195] ;

[0196] Where T represents operating cost;

[0197] The second objective function after considering the shared electricity purchase price and the shared electricity sales price is simulated and run using a power system model to obtain a second optimal solution and a second optimal planning scheme for the second objective function.

[0198] It can be understood that the official electricity purchase and sales costs are calculated through the current market, and the electricity exchange costs are calculated based on the usage of equipment in each new partition. Then, with the lowest operating cost as the goal, the second objective function is constructed, and the second constraint is set according to the demand. After considering the shared electricity purchase price and shared electricity sales price calculated in step S70, the power system model is used to simulate the operation of the distribution network to solve the second objective function and obtain the second optimal solution and the second optimal planning scheme. Because the shared electricity purchase price and shared electricity sales price are taken into account, when calculating the electricity exchange cost, the enthusiasm for energy interaction between the new partitions in the distribution network can be improved, and the advantages of flexible regulation can be brought into play, which will affect the two indicator powers of the first adjustable electric load and the second adjustable electric load, thereby reducing the operating costs in the process of solving the second objective function, and providing guidance for the coordinated control between the partitions in the distribution network. It should be noted that the first adjustable electric load is mainly an electric load that can be adjusted according to the time period, the second adjustable electric load is mainly for the stability of the distribution network and can be appropriately reduced. The second constraint condition is the industry standard, such as the storage constraint of electric energy and the balance constraint of electric power, which will not be elaborated here. The power system model can adopt the IEEE33 commonly used in the industry. After setting the parameters, it can simulate the operation of the distribution network to solve the second objective function.

[0199] In summary, the distribution control coordination method for multiple distribution areas of a distribution network in the above-mentioned embodiment of the present invention makes the partitioning of each sub-distribution point more reasonable and reduces the operating cost during subsequent use.

[0200] Please refer to Figure 2 , which shows a distribution network multi-zone distributed control coordination system in a second embodiment of the present invention, including:

[0201] A first calculation module 11 is configured to construct a power flow equation for a distribution network, and calculate a first fitness value of each sub-distribution point in the distribution network based on the power flow equation;

[0202] The second calculation module 12 is configured to traverse and calculate the second fitness value, the stability scale, and the intra-region matching scale between the sub-distribution point and the remaining sub-distribution points, and calculate a comprehensive fitness value by integrating the second fitness value, the stability scale, and the intra-region matching scale;

[0203] A first merging module 13 is configured to obtain a fitness difference value based on the comprehensive fitness value and the first fitness value, and to arrange the fitness differences between the remaining sub-distribution points and the initially selected sub-distribution point from largest to smallest, and to select the sub-distribution point corresponding to the largest fitness difference value and combine it with the initial sub-distribution point to obtain a pending merged partition;

[0204] The second merging module 14 uses the pending merged partition as a temporary sub-distribution point and sequentially absorbs the remaining sub-distribution points into the pending merged partition until a preset condition is met;

[0205] The third merging module 15 is used to sequentially calculate the pending merging partition results of the remaining sub-distribution points;

[0206] Optimization module 16: Based on the partitioned pending merged partition, calculate the carbon asset and equipment investment price, optimize the pending merged partition with the goal of maximizing the carbon asset acquisition rate to obtain a new partition;

[0207] The third calculation module 17 is configured to obtain the power demand based on the newly divided partitions and calculate the supply index of each time period in the distribution network to calculate the sharing index;

[0208] The fourth calculation module 18: obtains the official electricity purchase and sales costs and the electricity exchange costs, and obtains the optimal planning scheme with the goal of minimizing operating costs.

[0209] Throughout this specification, reference to terms such as "one embodiment," "some embodiments," "examples," "specific examples," or "some examples" means that a specific feature, structure, material, or characteristic described in conjunction with that embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0210] The above-described embodiments merely illustrate several implementations of the present invention, and while their descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that a person skilled in the art would be able to make numerous variations and improvements without departing from the spirit of the present invention, all of which fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention shall be determined by the appended claims.

Claims

1. A method for coordinated distributed control of multiple distribution areas in a distribution network, characterized in that: The steps include: Constructing a power flow equation of the distribution network, and calculating a first fitness value of each sub-distribution point in the distribution network based on the power flow equation; traversing and calculating a second fitness value, a stable scale, and an intra-region matching scale between a sub-distribution point and the remaining sub-distribution points, and calculating a comprehensive fitness value by integrating the second fitness value, the stable scale, and the intra-region matching scale; Obtaining a fitness difference according to the comprehensive fitness value and the first fitness value, and arranging the fitness differences between the remaining sub-distribution points and the initially selected sub-distribution point from largest to smallest, selecting the sub-distribution point corresponding to the largest fitness difference and combining it with the initial sub-distribution point to obtain a pending merged partition; The pending merged partition is used as a temporary sub-distribution point, and the remaining sub-distribution points are absorbed into the pending merged partition in sequence until the preset conditions are met; Calculate the pending merge partition results of the remaining sub-distribution points in sequence; Based on the partitioned pending merged partitions, calculating carbon assets and equipment investment prices, optimizing the pending merged partitions with the goal of maximizing the carbon asset acquisition rate to obtain new partitions; Based on the newly divided partitions, the power demand is obtained and the supply index of each time period in the distribution network is calculated to calculate the sharing index; Obtain the official electricity purchase and sales costs and electricity exchange costs, and obtain the optimal planning solution with the goal of minimizing operating costs.

2. The method for coordinated distributed control of multiple distribution areas in a distribution network according to claim 1, characterized in that: The constructing of a power flow equation for the distribution network and calculating a first fitness value of each sub-distribution point in the distribution network based on the power flow equation specifically includes: Based on the power flow parameters of the distribution network, the voltage amplitude changes when all the sub-distribution points are connected to active equipment and reactive equipment are calculated respectively. and ; ； in, represents the change in the voltage phase angle in the sub-distribution point, represents the change in the voltage amplitude value in the sub-distribution point, represents the change in active power in the sub-distribution point, represents the change in reactive power in the sub-distribution point, Indicates the change in voltage phase angle when the sub-distribution point is connected to active equipment, Indicates the change in voltage phase angle when the sub-distribution point is connected to reactive equipment; According to the change of the voltage amplitude value when the equipment in the sub-distribution point is active and reactive and , calculate the distance weight of the sub-distribution point; ； ； in, It represents the influence index of the voltage amplitude value of sub-distribution point i when other sub-distribution points j connected to sub-distribution point i are connected to reactive power. It indicates the influence index of the voltage amplitude value of other connected sub-distribution points j when the sub-distribution point i is connected to reactive power. It represents the influence index of the voltage amplitude value of the sub-distribution point i when other sub-distribution points j connected to the sub-distribution point i are connected to the active power, It indicates the influence index of the voltage amplitude value of other connected sub-distribution points j when the sub-distribution point i is connected to the active power. represents the reactive distance weight between two sub-distribution points, Represents the active distance weight between two sub-distribution points; The reactive power influence factor ω and the active power influence factor (1-ω) are set, and the active power distance weight and the reactive power distance weight are integrated to obtain the overall weight; ； Normalizing the overall weight to obtain the edge weight of the sub-distribution point; Calculating a first fitness value of the sub-distribution point according to the edge weight of the sub-distribution point; ； ； in, and Represents the total weight of sub-distribution point i and sub-distribution point j, m V Represents the total edge weight of the sub-distribution points. If the two sub-distribution points are connected, then ρ(i, j) = 1, otherwise ρ(i, j) = 0. Represents edge weight.

3. The method for coordinated distributed control of multiple distribution areas in a distribution network according to claim 2, characterized in that: The traversal calculation of the second fitness value, the stable scale, and the intra-region matching scale between the sub-distribution point and the remaining sub-distribution points, and the calculation of the comprehensive fitness value by integrating the second fitness value, the stable scale, and the intra-region matching scale, specifically includes: Select a sub-distribution point and traverse to calculate the second fitness value between the sub-distribution point i and the remaining sub-distribution points q ; ； Obtain all the equipment in the sub-distribution point i and the remaining sub-distribution points q to obtain the reactive power supply N Q ; Get the voltage over-limit value and get the minimum reactive power supply N W ; Based on the reactive power supply amount N Q , Minimum reactive power supply N W , and reactive distance weight , calculate the stationary scale between sub-distribution point i and sub-distribution point q; ； Based on the edge weights and the current partitioning situation, calculate the first intra-region matching index σ1 and the second intra-region matching index σ2; ； ； Among them, C represents the serial number of the partition, G represents the total number of partitions, G C Indicates C partition, N C Indicates the total number of neutron distribution points in the C partition, B indicates the total number of neutron distribution points in this calculation, represents the reactive distance weight between sub-distribution point i and sub-distribution point q, represents the edge weight between sub-distribution point i and sub-distribution point q; Calculating an intra-region matching scale σ according to the first intra-region matching index and the second intra-region matching index; ； Calculating a comprehensive fitness value by combining the second fitness value, the stability scale, and the intra-region matching scale; 。 4. The method for coordinated distributed control of multiple distribution areas in a distribution network according to claim 3, characterized in that: The step of obtaining a fitness difference based on the comprehensive fitness value and the first fitness value, arranging the fitness differences between the remaining sub-distribution points and the initially selected sub-distribution point from large to small, selecting the sub-distribution point corresponding to the largest fitness difference and combining it with the initial sub-distribution point to obtain a pending merged partition, specifically includes: Subtract the comprehensive fitness value between the sub-distribution point i and the remaining sub-distribution points q from the first fitness value of the sub-distribution point i to obtain a fitness difference; Arrange the fitness differences between the remaining sub-distribution points q and the initially selected sub-distribution point i from large to small; The sub-distribution point q corresponding to the largest fitness difference is selected and combined with the initial sub-distribution point i to obtain a pending merged partition.

5. The method for coordinated distributed control of multiple distribution areas in a distribution network according to claim 4, characterized in that: The step of using the pending merged partition as a temporary sub-distribution point and sequentially absorbing the remaining sub-distribution points into the pending merged partition until a preset condition is met specifically includes: The pending merged partition is used as a temporary sub-distribution point iq; Subtract the comprehensive fitness value between the temporary sub-distribution point iq and the remaining sub-distribution points q that have not been merged from the comprehensive fitness value of the temporary sub-distribution point iq to obtain a fitness difference; All other sub-distribution points q that have not been merged are absorbed into the pending merge partition in turn until the fitness difference becomes negative.

6. The method for coordinated distributed control of multiple distribution areas in a distribution network according to claim 1, characterized in that: The method of calculating the carbon asset and equipment investment price based on the partitioned pending merged partition, optimizing the pending merged partition with the goal of maximizing the carbon asset acquisition rate to obtain a new partition, specifically includes: Calculating carbon assets based on market factors and planning schemes, and calculating equipment investment prices for the to-be-determined merged subareas based on the status of access equipment at sub-distribution points in the to-be-determined merged subareas; ； ； Among them, C Y represents carbon assets, D represents the number of typical load calculation days covered in the year, and W T represents the instantaneous price of carbon market, Tθ represents the number of time periods covering the typical annual load calculation day, M is the number of all sub-distribution points in the distribution network, and W a,k,t represents the output power of the kth sub-distribution point in time period t, γ1 represents the carbon dioxide emission index, C J represents the equipment investment price, G inv Represents the annual investment cost of the distribution network, G o represents the annual electricity transaction fee, G l Represents the annual network loss cost, G b Indicates annual maintenance fee; Based on the carbon assets and the equipment investment price, with the goal of maximizing the carbon asset acquisition rate, a first objective function is constructed, and a first constraint condition is set according to demand; ； Among them, P C represents the carbon asset acquisition rate; Solving the first objective function by a machine algorithm to obtain a first optimal solution and a first optimal planning scheme for the first objective function; Preset a first optimal solution threshold. If the first optimal solution is smaller than the first optimal solution threshold, remove the sub-distribution points in the reverse order of the construction of the pending merged partitions, and cyclically calculate the first optimal solution based on the updated pending merged partitions until the first optimal solution is larger than the first optimal solution threshold. When the first optimal solutions of the first objective function solved based on the partitions to be determined and merged converge, the final partition to be determined and merged is determined as a new partition.

7. The method for coordinated distributed control of multiple distribution areas in a distribution network according to claim 1, characterized in that: The method of obtaining the power demand based on the newly divided partition and calculating the supply index of each time period in the distribution network to calculate the sharing index specifically includes: Based on the device access status of each sub-distribution point in each new partition, obtaining the power demand in each new partition; Calculating the supply index of each time period in the distribution network based on the power demand in each of the new zones; ； in, represents the electricity sales in the Nth new partition during period t, represents the amount of electricity purchased in the Nth new partition during period t, and tg(t) represents the supply index during period t; Calculate the shared median price based on the current market carbon emission factors; ； in, Indicates the electricity purchase price of the upper power grid, Indicates the electricity price of the upper power grid, G P Represents the carbon emission index, G Pmax represents the maximum value of carbon emission index, ε is a constant, represents the shared median price; Optimizing the shared median price, and calculating a shared electricity purchase price and a shared electricity sales price based on the optimized shared median price; ； ； ； in, represents the shared correction, tc=1 / tg, Indicates the shared electricity purchase price after optimization, Indicates the shared electricity sales price after optimization.

8. The method for coordinated distributed control of multiple distribution areas in a distribution network according to claim 7, characterized in that: The above-mentioned acquisition of official electricity purchase and sales costs and electricity exchange costs, with the goal of minimizing operating costs, results in the optimal planning scheme, specifically including: Obtain official electricity purchase and sales costs and electricity exchange costs; ； ； Among them, T Y represents the official cost of electricity purchase and sale, Indicates the amount of electricity purchased by the upper power grid. Indicates the power sales of the upper power grid, D c Indicates the time of day, T η represents the electricity exchange cost, η g represents the energy storage penalty coefficient, Indicates the energy storage charging power, Represents the energy storage discharge power, η k represents the first adjustable electric load penalty coefficient, represents the first adjustable electric load power, η l represents the second adjustable electric load penalty coefficient, Indicates the second adjustable electric load power; With the goal of minimizing operating costs, construct the second objective function and set the second constraint conditions according to requirements; ； Where T represents operating cost; The second objective function after considering the shared electricity purchase price and the shared electricity sales price is simulated and run using a power system model to obtain a second optimal solution and a second optimal planning scheme for the second objective function.

9. A distribution network multi-zone distributed control coordination system, characterized in that: include: A first calculation module is used to construct a power flow equation of the distribution network, and calculate a first fitness value of each sub-distribution point in the distribution network based on the power flow equation; A second calculation module is configured to traverse and calculate a second fitness value, a stationary scale, and an intra-region matching scale between a sub-distribution point and the remaining sub-distribution points, and calculate a comprehensive fitness value by integrating the second fitness value, the stationary scale, and the intra-region matching scale; A first merging module is configured to obtain a fitness difference according to the comprehensive fitness value and the first fitness value, and arrange the fitness differences between the remaining sub-distribution points and the initially selected sub-distribution point from large to small, and select the sub-distribution point corresponding to the largest fitness difference and combine it with the initial sub-distribution point to obtain a pending merged partition; The second merging module is configured to use the pending merged partition as a temporary sub-distribution point and sequentially absorb the remaining sub-distribution points into the pending merged partition until a preset condition is met; A third merging module is used to sequentially calculate the pending merging partition results of the remaining sub-distribution points; Optimization module: used for calculating carbon assets and equipment investment prices based on the partitioned pending merged partitions, optimizing the pending merged partitions with the goal of maximizing the carbon asset acquisition rate to obtain new partitions; A third calculation module is configured to obtain the power demand based on the newly divided partitions and calculate the supply index of each time period in the distribution network to calculate the sharing index; The fourth calculation module is used to obtain the official electricity purchase and sales costs and electricity exchange costs, with the goal of minimizing operating costs to obtain the optimal planning solution.

Images