Distribution control cooperation method and system for multiple transformer areas of power distribution network

Abstract

The invention provides a multi-zone-area distribution control cooperation method and system for a power distribution network, and the method comprises the steps: constructing a power flow equation, and calculating a first fitness value of a distribution point in the power distribution network; calculating a second fitness value, a stationary scale and an intra-region matching scale between one sub-distribution point and other sub-distribution points, and calculating a comprehensive fitness value; fitness difference values are obtained through calculation, the fitness difference values between the other sub-distribution points and the initial sub-distribution points are ranked from large to small, and the fitness difference values are combined from large to small to obtain a to-be-combined partition; sequentially absorbing the remaining sub-distribution points to the to-be-combined partition until a preset condition is met; based on the to-be-combined partitions, calculating carbon assets and equipment investment prices, and optimizing the to-be-combined partitions to obtain new partitions; calculating a supply index of each time period of the power distribution network based on the new partition so as to calculate a sharing index; and obtaining official electricity purchasing and selling cost and power exchange cost, and obtaining an optimal planning scheme. According to the invention, the partitioning is more reasonable, and the subsequent operation cost is lower.

Description

Technical Field

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

[0002] Due to the development of new energy vehicles and electrical appliances, a large number of access nodes will be newly connected in a large-scale distribution network, which will in turn lead to the imbalance of the distribution system. In order to reduce the impact of this imbalance on the safety and lifespan of the distribution system, it is necessary to perform distributed combined control on nodes with high electrical coupling degree and strong correlation.

[0003] In the prior art, the research on the correlation between sub-distribution points in a distribution network generally predicts future operating parameters based on historical data, generally only using the influence degree between sub-distribution points in the distribution network as a reference standard. The sample data is less and the prediction algorithm is not stable enough. Also, the influence brought by constraint conditions to nodes is not considered, and the shared coordinability of the distribution network after zoning is not fully utilized to optimize the implementation plan. Summary of the Invention

[0004] Aiming at the deficiencies of the prior art, the purpose of the present invention is to provide a method for coordinated distribution control of multiple substations in a distribution network, aiming to solve the technical problems mentioned in the background art.

[0005] To achieve the above purpose, the present invention is implemented through the following technical solutions: A method for coordinated distribution control of multiple substations in a distribution network includes the following steps: Construct a power flow equation of the distribution network, and calculate the first fitness value of each sub-distribution point in the distribution network based on the power flow equation; Traverse and calculate the second fitness value, steady scale, and intra-region matching scale between one sub-distribution point and the remaining sub-distribution points, and calculate the comprehensive fitness value by synthesizing the second fitness value, steady scale, and intra-region matching scale; According to the comprehensive fitness value and the first fitness value, obtain the fitness difference, 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 the initial sub-distribution point to form a to-be-merged partition; Take the to-be-merged partition as a temporary sub-distribution point, and sequentially absorb the remaining sub-distribution points into the to-be-merged partition until a preset condition is met; Calculate the to-be-merged partition results of the remaining sub-distribution points in sequence; Based on the divided to-be-merged partitions, calculate the carbon assets and equipment investment prices, and optimize the to-be-merged partitions with the highest carbon asset acquisition rate as the goal to obtain new partitions; Based on the divided new partitions, obtain the electricity demand and calculate the supply indicators for each time period in the distribution network to calculate the sharing indicators; Obtain the official electricity purchase and sale costs and power exchange costs, and take the lowest operating cost as the goal to obtain the optimal planning scheme.

[0006] 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: Based on the power flow parameters of the distribution network, calculate the voltage magnitude change when all the sub-distribution points are connected to active and reactive devices respectively and ; ; Among them, represents the change in the voltage phase angle in the sub-distribution point, represents the change in the voltage magnitude in the sub-distribution point, represents the change in the active power in the sub-distribution point, represents the change in the reactive power in the sub-distribution point, represents the change in the voltage phase angle when the sub-distribution point is connected to an active device, represents the change in the voltage phase angle when the sub-distribution point is connected to a reactive device; According to the voltage magnitude changes and when the devices in the sub-distribution point are active and reactive, calculate the distance weight of this sub-distribution point; ; ; Among them, represents the influence index of the voltage magnitude of the sub-distribution point i when the other sub-distribution point j connected to the sub-distribution point i accesses reactive power, represents the influence index of the voltage magnitude of the other sub-distribution point j connected to the sub-distribution point i when the sub-distribution point i accesses reactive power, represents the influence index of the voltage magnitude of the sub-distribution point i when the other sub-distribution point j connected to the sub-distribution point i accesses active power, represents the influence index of the voltage magnitude of the other sub-distribution point j connected to the sub-distribution point i when the sub-distribution point i accesses active power, represents the reactive distance weight between two sub-distribution points, represents the active distance weight between two sub-distribution points; Set the reactive power influence factor ω and the active power influence factor (1 - ω), and combine the active power distance weight and the reactive power distance weight to obtain the overall weight; ; Normalize the overall weight to obtain the edge weight of the sub - distribution point; Calculate the first fitness value of the sub - distribution point through the edge weight of the sub - distribution point; ; ; Among them, and respectively represent the total weights of sub - distribution point i and sub - distribution point j, m V represents the total edge weight of the sub - distribution point. If the two sub - distribution points are connected, then ρ(i, j)=1, otherwise ρ(i, j)=0, represents the edge weight.

[0007] According to one aspect of the above - mentioned technical solution, traverse and calculate the second fitness value, the stability scale, and the in - area matching scale between one sub - distribution point and the other sub - distribution points, and calculate the comprehensive fitness value by combining the second fitness value, the stability scale, and the in - area matching scale. Specifically, it includes: Select a sub - distribution point, and traverse and calculate the second fitness value between this sub - distribution point i and the other sub - distribution points q ; ; Obtain all the devices among the sub - distribution point i and the other sub - distribution points q to obtain the reactive power supply quantity N Q ; Obtain the over - limit value of the voltage to obtain the minimum reactive power supply quantity N W ; Based on the reactive power supply quantity N Q , the minimum reactive power supply quantity N W , and the reactive power distance weight , calculate the stability scale between sub - distribution point i and sub - distribution point q; ; Based on the edge weight and the current partition situation, calculate the first in - area matching index and the second in - area matching index; ; ; Among them, C represents the serial number of the partition, G represents the total number of partitions, G C represents the C partition, N CLet \(C\) represent the total number of sub - distribution points in the \(C\) partition, and \(B\) represent the total number of sub - distribution points in this calculation. Let \(w_{ij}\) represent the reactive power distance weight between sub - distribution point \(i\) and sub - distribution point \(q\). Let \(e_{ij}\) represent the edge weight between sub - distribution point \(i\) and sub - distribution point \(q\). Calculate the in - area matching scale according to the matching index in the first area and the matching index in the second area. ; Calculate the comprehensive fitness value by synthesizing the second fitness value, the stability scale, and the in - area matching scale. .

[0008] According to one aspect of the above - mentioned technical solution, based on the comprehensive fitness value and the first fitness value, obtain the fitness difference, and arrange the fitness differences between each of 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 the initial sub - distribution point to form a to - be - merged partition. Specifically, it includes: Subtract the first fitness value of sub - distribution point \(i\) from the comprehensive fitness value between sub - distribution point \(i\) and the remaining sub - distribution points \(q\) to obtain the fitness difference. Arrange the fitness differences between each of the remaining sub - distribution points \(q\) and the initially selected sub - distribution point \(i\) from large to small. Select the sub - distribution point \(q\) corresponding to the largest fitness difference and the initial sub - distribution point \(i\) to form a to - be - merged partition.

[0009] According to one aspect of the above - mentioned technical solution, take the to - be - merged partition as a temporary sub - distribution point, and sequentially absorb the remaining sub - distribution points into the to - be - merged partition until a preset condition is met. Specifically, it includes: Take the to - be - merged partition as a temporary sub - distribution point \(iq\). Subtract the comprehensive fitness value of the temporary sub - distribution point \(iq\) from the comprehensive fitness value between the temporary sub - distribution point \(iq\) and the remaining un - merged sub - distribution points \(q\) to obtain the fitness difference. And sequentially absorb all the remaining un - merged sub - distribution points \(q\) into the to - be - merged partition until the fitness difference becomes negative.

[0010] According to one aspect of the above - mentioned technical solution, based on the divided to - be - merged partition, calculate the carbon assets and equipment investment prices, and optimize the to - be - merged partition with the goal of the highest carbon asset acquisition rate to obtain a new partition. Specifically, it includes: Calculate carbon assets based on market factors and planning schemes, and calculate the equipment investment price of the to-be-determined merger partition based on the access equipment situation of the sub-distribution points in the to-be-determined merger partition; ; ; Among them, C Y represents carbon assets, D represents the number of typical load calculation days in a year, W T represents the instant price of the carbon market, Tθ represents the number of time periods within the typical load calculation day of the year, M is the number of all sub-distribution points of the distribution network, W a,k,t represents the output power of the kth sub-distribution point in the t time period, γ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 trading cost, G l represents the annual network loss cost, G b represents the annual maintenance charge; Based on the carbon assets and the equipment investment price, with the goal of maximizing the carbon asset acquisition rate, construct the first objective function and set the first constraint conditions according to the requirements; ; Among them, P C represents the carbon asset acquisition rate; Solve the first objective function through a machine algorithm to obtain the first optimal solution and the first optimal planning scheme of the first objective function; Preset the first optimal solution threshold. If the first optimal solution is less than the first optimal solution threshold, remove the sub-distribution points in the reverse order of the construction of the to-be-determined merger partition, and calculate the first optimal solution based on the updated to-be-determined merger partition in a loop until the first optimal solution is greater than the first optimal solution threshold; When the first optimal solutions of the first objective function solved based on each to-be-determined merger partition converge, determine the final to-be-determined merger partition as the new partition.

[0011] According to one aspect of the above technical solution, based on the divided new partition, obtain the electricity demand and calculate the supply index of each time period in the distribution network to calculate the sharing index, specifically including: Based on the equipment access situation of each sub-distribution point in each new partition, obtain the electricity demand in each new partition; Based on the electricity demand in each new partition, calculate the supply index of each time period in the distribution network; ; Among them, represents the electricity sales volume in the Nth new partition during the t period, represents the electricity purchase volume in the Nth new partition during the t period, and tg(t) represents the supply index during the t period; Based on the current market carbon emission factor, calculate the shared median price; ; Among them, represents the purchase price of electricity from the superior power grid, represents the selling price of electricity to the superior power grid, G P represents the carbon emission index, G Pmax represents the maximum value of the carbon emission index, and ε is a constant, represents the shared median price; Optimize the shared median price, and calculate the shared electricity purchase price and the shared electricity selling price based on the optimized shared median price.

[0012] ; ; ; Among them, represents the shared correction amount, tc = 1 / tg, represents the optimized shared electricity purchase price, represents the optimized shared electricity selling price.

[0013] According to one aspect of the above technical solution, the obtaining of the official electricity purchase and selling costs and the electricity exchange cost, with the goal of the lowest operating cost, to obtain the optimal planning scheme, specifically includes: Obtain the official electricity purchase and selling costs and the electricity exchange cost; ; ; Among them, T Y represents the official electricity purchase and selling cost, represents the electricity purchase volume from the superior power grid, represents the electricity selling volume to the superior power grid, D c represents the moment in a day, T η represents the electricity exchange cost, η g represents the energy storage penalty coefficient, represents the energy storage charging power, represents the energy storage discharging 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, represents the second adjustable electric load power; With the goal of minimizing operating costs, a second objective function is constructed, and second constraint conditions are set according to requirements; ; Among them, T represents the operating cost; The second objective function after considering the shared power purchase price and shared power sales price is simulated and operated using a power system model to obtain the second optimal solution and the second best planning scheme of the second objective function.

[0014] The present invention also provides a distribution network multi-substation area distribution control cooperation system, including: The first calculation module: used to construct the power flow equation of the distribution network, and calculate the first fitness value of each sub-distribution point in the distribution network based on the power flow equation; The second calculation module: used to traverse and calculate the second fitness value, the smooth scale, and the in-area matching scale between one sub-distribution point and the remaining sub-distribution points, and calculate the comprehensive fitness value by synthesizing the second fitness value, the smooth scale, and the in-area matching scale; The first merging module: used to obtain the fitness difference according to the comprehensive fitness value and the first fitness value, 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 the initial sub-distribution point for combination to obtain a to-be-merged partition; The second merging module: used to use the to-be-merged partition as a temporary sub-distribution point, and sequentially absorb the remaining sub-distribution points into the to-be-merged partition until a preset condition is met; The third merging module: used to sequentially calculate the to-be-merged partition results of the remaining sub-distribution points; The optimization module: used to calculate the carbon assets and equipment investment prices based on the divided to-be-merged partition, and optimize the to-be-merged partition with the goal of maximizing the carbon asset acquisition rate to obtain a new partition; The third calculation module: used to obtain the electricity demand based on the divided new partition and calculate the supply index of each time period in the distribution network to calculate the sharing index; The fourth calculation module: used to obtain the official power purchase and sales costs and power exchange costs, and obtain the optimal planning scheme with the goal of minimizing the operating cost.

[0015] Compared with the prior art, the beneficial effects of the present invention are as follows: By calculating the first fitness value of each sub-distribution point, traversing and calculating the second fitness value, the steady scale, and the in-region matching scale between one sub-distribution point and the rest of the sub-distribution points, a comprehensive fitness value can be obtained. By calculating the comprehensive fitness value of two sub-distribution points and the first fitness value of the initial sub-distribution point itself, a fitness difference can be obtained. According to the fitness difference, the fitness relationship between the initial sub-distribution point and the rest of the sub-distribution points can be determined, and then the merging order can be determined. Each time the initial sub-distribution point merges with one of the remaining sub-distribution points, it is marked as a pending merging partition. The pending merging partition is used as a temporary sub-distribution point, and the remaining unmerged sub-distribution points are merged in turn until the preset conditions are met. Then, another initial sub-distribution point is selected and operated according to the above steps until all pending merging partitions are obtained. After this step, the initial merging of the sub-distribution points is completed; then, by calculating the carbon assets and equipment investment prices of the pending merging partitions, if they exceed the preset target, the sub-distribution points in the pending merging partitions are subtracted in turn, and the calculation is repeated until the requirements are met, and then new partitions are obtained; then, the sharing index is calculated according to the operation conditions of the sub-distribution points in the newly divided partitions. Finally, with the goal of the lowest operating cost, the optimal planning scheme between the sub-distribution points in the new partitions is calculated.

[0016] When partitioning in the present invention, the steady scale and the in-region matching scale between sub-distribution points are fully considered, and then the comprehensive fitness value is obtained, which can make the partitioning between sub-distribution points more reasonable, with a higher matching degree and stronger applicability. Secondly, by calculating the carbon assets and equipment investment prices of the pending merging partitions, a better balance can be achieved between the partitioning matching degree and the carbon emission cost, avoiding over-optimization in one aspect while ignoring the other aspect. The sharing index is optimally incorporated to calculate the optimal planning scheme between the sub-distribution points in the new partitions, realizing the coordinated control of the distribution network system. Brief Description of the Drawings

[0017] Figure 1 It is a flowchart of the coordinated method for multi-substation area distribution control in the first embodiment of the present invention; Figure 2 It is a structural block diagram of the coordinated system for multi-substation area distribution control in the second embodiment of the present invention; The following specific embodiments will further illustrate the present invention in combination with the above-mentioned drawings. Specific Embodiments

[0018] To facilitate the understanding of the present invention, the present invention will be described more comprehensively below with reference to the relevant drawings. Several embodiments of the present invention are given in the drawings. However, the present invention can be implemented in many different forms and is not limited to the embodiments described herein. On the contrary, the purpose of providing these embodiments is to make the disclosure of the present invention more thorough and comprehensive.

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

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

[0021] Please refer to Figure 1 , which shows the multi-substation area distribution control coordination method in the first embodiment of the present invention, including the following steps: S10. Construct a power flow equation of the distribution network, and calculate the first fitness value of each sub-distribution point in the distribution network based on the power flow equation; S20. Traverse and calculate the second fitness value, stability scale, and in-area matching scale between one sub-distribution point and the other sub-distribution points, and calculate the comprehensive fitness value by synthesizing the second fitness value, stability scale, and in-area matching scale; S30. Obtain a fitness difference according to the comprehensive fitness value and the first fitness value, arrange the fitness differences between the other 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 the initial sub-distribution point to form a to-be-merged partition; S40. Take the to-be-merged partition as a temporary sub-distribution point, and sequentially absorb the remaining sub-distribution points into the to-be-merged partition until a preset condition is met; S50. Sequentially calculate the to-be-merged partition results of the remaining sub-distribution points; S60. Based on the divided to-be-merged partition, calculate the carbon assets and equipment investment prices, and optimize the to-be-merged partition with the highest carbon asset acquisition rate as the goal to obtain a new partition; S70. Based on the divided new partition, obtain the power demand and calculate the supply index for each time period in the distribution network to calculate the sharing index; S80. Obtain the official power purchase and sale cost and power exchange cost, and obtain the optimal planning scheme with the lowest operating cost as the goal.

[0022] It can be understood that in the present invention, by calculating the first fitness value of each sub-distribution point and traversing to calculate the second fitness value, steady-state scale, and in-region matching scale between one sub-distribution point and the remaining sub-distribution points, a comprehensive fitness value can be obtained. By calculating the comprehensive fitness value between two sub-distribution points and the first fitness value of the initial sub-distribution point itself, a fitness difference can be obtained. Based on the fitness difference, the fitness relationship between the initial sub-distribution point and the remaining sub-distribution points can be determined, and then the merging order can be determined. Each time the initial sub-distribution point merges with one of the remaining sub-distribution points, it is marked as a pending merge partition. The pending merge partition is used as a temporary sub-distribution point, and the remaining unmerged sub-distribution points are merged in sequence until the preset conditions are met. Then, another initial sub-distribution point is selected and operated according to the above steps until all pending merge partitions are obtained. After this step, the initial merge of the sub-distribution points is completed; then, the carbon assets and equipment investment prices of the pending merge partition are calculated. If they exceed the preset target, the sub-distribution points in the pending merge partition are subtracted in sequence, and the calculation is repeated until the requirements are met, thereby obtaining a new partition; then, the sharing index is calculated based on the operation conditions of the sub-distribution points in the newly divided partition. Finally, with the goal of minimizing the operating cost, the optimal planning scheme between the sub-distribution points in the new partition is calculated.

[0023] When partitioning in the present invention, the steady-state scale and in-region matching scale between sub-distribution points are fully considered, and then the comprehensive fitness value is obtained, which can make the partitioning between sub-distribution points more reasonable, with a higher matching degree and stronger applicability. Secondly, by calculating the carbon assets and equipment investment prices of the pending merge partition, a better balance can be achieved between the partitioning matching degree and the carbon emission cost, avoiding one aspect being too excellent while ignoring the other aspect. Finally, the sharing index is incorporated to calculate the optimal planning scheme between the sub-distribution points in the new partition, realizing the coordinated control of the distribution network system.

[0024] Further, the step S10 specifically includes: Based on the power flow parameters of the distribution network, calculate the voltage magnitude change when all the sub-distribution points are connected to active devices and reactive devices and ; ; wherein, represents the change in voltage phase angle in the sub-distribution point, represents the change in voltage magnitude 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, represents the change in voltage phase angle when the sub-distribution point is connected to an active device, Indicates the change in voltage phase angle when the sub-distribution point accesses the reactive power device; According to the change in voltage amplitude value when there is active and reactive power in the devices at the sub-distribution point and , calculate the distance weight of this sub-distribution point; ; ; Among them, Indicates 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 access reactive power, Indicates the influence index of the voltage amplitude value of other sub-distribution points j connected to sub-distribution point i when sub-distribution point i accesses reactive power, Indicates the influence index of the voltage amplitude value of sub-distribution point i when other sub-distribution points j connected to the sub-distribution point i access active power, Indicates the influence index of the voltage amplitude value of other sub-distribution points j connected to sub-distribution point i when sub-distribution point i accesses active power, Indicates the reactive distance weight between two sub-distribution points, Indicates the active distance weight between two sub-distribution points; Set the reactive influence factor ω and the active influence factor (1 - ω), and combine the active distance weight and the reactive distance weight to obtain the overall weight; ; Normalize the overall weight to obtain the edge weight of the sub-distribution point; Calculate the first fitness value of the sub-distribution point through the edge weight of the sub-distribution point; ; ; Among them, and respectively represent the total weights of sub-distribution point i and sub-distribution point j, and m V represents the total marginal weight of the sub-distribution point. If the two sub-distribution points are connected, then ρ(i, j) = 1, otherwise ρ(i, j) = 0, represents the edge weight.

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

[0026] Further, the specific steps of step S20 include: 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 devices among the sub - distribution point i and the remaining sub - distribution points q to get the reactive power supply capacity N Q ; Obtain the over - limit value of the voltage to get the minimum reactive power supply quantity N W ; Based on the reactive power supply capacity N Q , the minimum reactive power supply quantity N W , and the reactive power distance weight , calculate the steady - state scale between the sub - distribution point i and the sub - distribution point q; ; Based on the edge weight and the current partition situation, calculate the matching index in the first area and the matching index in the second area; ; ; Among them, C represents the serial number of the partition, G represents the total number of partitions, G C represents the C - th partition, N C represents the total number of sub - distribution points in the C - th partition, B represents the total number of sub - distribution points in this calculation, represents the reactive power distance weight between the sub - distribution point i and the sub - distribution point q, represents the edge weight between the sub - distribution point i and the sub - distribution point q; According to the matching index in the first area and the matching index in the second area, calculate the intra - area matching scale; ; Comprehensively calculate the comprehensive fitness value based on the second fitness value, the steady - state scale, and the intra - area matching scale; .

[0027] It can be understood that the calculation method of the second fitness value is the same as that of the first fitness value. When calculating the first fitness value, it is necessary to calculate the distance weight between the initially selected sub-distribution points and the remaining adjacent sub-distribution points. When calculating the second fitness value, it is necessary to calculate the distance weight between the initially selected sub-distribution points and all the remaining sub-distribution points in the distribution network. Therefore, some steps are omitted in this step. For the omitted part, refer to the content of calculating the first fitness value. The second fitness value can reflect the matching degree between the initially selected sub-distribution points and the remaining sub-distribution points. The higher the second fitness value, the higher the matching degree. Then, according to the available reactive power supply N Q , the minimum reactive power supply N W Calculate the steady scale between two sub-distribution points. The steady scale reflects the reactive power balance ability between two sub-distribution points. Then, according to the edge weight and the current partition situation, calculate the in-zone matching scale. The in-zone matching scale reflects the coordination degree between each sub-distribution point in the partition under this partition situation, improving the accuracy of the partition. The above two indicators are two other important indicators for judging between two sub-distribution points, avoiding poor partition effects caused by only considering the voltage amplitude value influence index between two sub-distribution points.

[0028] Further, the step S30 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 value; Arrange the fitness difference values between the remaining sub-distribution points q and the initially selected sub-distribution point i from large to small; Select the sub-distribution point q corresponding to the largest fitness difference value and combine it with the initial sub-distribution point i to obtain a to-be-merged partition.

[0029] It can be understood that by taking the difference, the fitness difference value between the initially selected sub-distribution point and the remaining sub-distribution points can be calculated. The fitness difference value reflects "how much the second fitness value and the first fitness value will change if the initial sub-distribution point is combined with one of the remaining sub-distribution points". If the fitness difference value is positive and the value is larger, it means that the effect of this combination is very good. Therefore, the fitness difference value is used to determine the combination order between sub-distribution points.

[0030] Further, the step S40 specifically includes: Take the to-be-merged partition as a temporary sub-distribution point iq; Subtract the comprehensive fitness value between the temporary sub-distribution point iq and the remaining un-merged sub-distribution points q from the comprehensive fitness value of the temporary sub-distribution point iq to obtain a fitness difference value; And then absorb all the remaining unmerged sub - distribution points q into the to - be - determined merging partition in sequence until the fitness difference appears negative.

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

[0032] Further, step S60 specifically includes: Calculate carbon assets based on market factors and planning schemes, and calculate the equipment investment price of the to - be - determined merging partition based on the situation of the access equipment of the sub - distribution points in the to - be - determined merging partition; ; ; Among them, C Y represents carbon assets, D represents the number of days of the typical load calculation day in a year, W T represents the instant price of the carbon market, Tθ represents the number of time periods within the typical load calculation day in a year, M is the number of all sub - distribution points of the distribution network, W a,k,t represents the output power of the k - th sub - distribution point in the t - th time period, γ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 trading cost, G l represents the annual network loss cost, G b represents the annual maintenance charge; Specifically, ; ; ; ; Among them, q represents the discount rate, R a represents the planned years of the equipment in the sub - distribution point, e1 represents the investment cost of the equipment in the sub - distribution point, e oa represents the maintenance charge of the equipment in the sub - distribution point, e p represents the average purchase price of the transmission network, P s,k,t represents the output power of the thermal power grid at the i - th sub - distribution point in the t - th time period, e l represents the network loss cost coefficient, Lkc Indicates the current of line c of device k within the sub-distribution point, C kc Indicates the impedance of line c of device k within the sub-distribution point; Based on the carbon assets and the device investment price, with the goal of maximizing the carbon asset acquisition rate, construct a first objective function and set first constraint conditions according to requirements; ; Among them, P C Indicates the carbon asset acquisition rate; Solve the first objective function through a machine algorithm to obtain the first optimal solution and the first optimal planning scheme of the first objective function; Preset a first optimal solution threshold. If the first optimal solution is less than the first optimal solution threshold, remove the sub-distribution point in the reverse order of the construction of the to-be-determined merged partition, and loop to calculate the first optimal solution based on the updated to-be-determined merged partition until the first optimal solution is greater than the first optimal solution threshold; When the first optimal solutions of the first objective function solved based on each to-be-determined merged partition converge, determine the final to-be-determined merged partition as the new partition.

[0033] It is understandable that after all the to-be-merged partitions are determined, the carbon assets and equipment investment prices of the entire distribution network operation under this partition condition are calculated based on the operation conditions of the devices at the sub-distribution points in the to-be-merged partitions. Then, based on the carbon assets and equipment investment prices, with the goal of maximizing the carbon asset acquisition rate, a first objective function is constructed, and the first constraint conditions are set according to the requirements to solve the first objective function. The solution process can be implemented using machine algorithms, such as the SSA algorithm. It should be noted that the first constraint conditions are industry norms, such as distribution network power flow constraints, power balance constraints, and equipment capacities at sub-distribution points, which will not be elaborated here. After obtaining the first optimal solution and the first optimal planning scheme by solving the first objective function, the first optimal solution threshold can be determined according to the requirements. That is, if the obtained maximum carbon asset acquisition rate is less than the target, some nodes in the to-be-merged partitions need to be abandoned in the reverse order of merging, and the calculation is repeated until the first optimal solution meets the first optimal solution threshold. Then, the to-be-merged partitions are determined as new partitions, and the abandoned nodes need to be recalculated and incorporated into the appropriate to-be-merged partitions according to the steps of S30 and S40, and then the first optimal solutions of the remaining to-be-merged partitions after the update are solved to obtain all the new partitions. 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 so that all abandoned sub-distribution points can be incorporated into new partitions as much as possible in the remaining to-be-merged partitions. Machine algorithms can be used to achieve this (such as the particle swarm optimization algorithm).

[0034] Further, the step S70 specifically includes: Based on the device access conditions of each sub-distribution point in each of the new partitions, obtain the electricity demand in each of the new partitions; Based on the electricity demand in each of the new partitions, calculate the supply index for each time period in the distribution network; ; Wherein, represents the electricity sales volume in the Nth new partition at time t, represents the electricity purchase volume in the Nth new partition at time t, and tg(t) represents the supply index at time t; Based on the current market carbon emission factor, calculate the shared median price; ; Wherein, represents the purchase price of electricity from the superior power grid, represents the selling price of electricity to the superior power grid, G P represents the carbon emission index, G Pmax represents the maximum value of the carbon emission index, ε is a constant, Represents the shared median price; Optimize the shared median price, and calculate the shared purchase electricity price and the shared selling electricity price based on the optimized shared median price.

[0035] ; ; ; Among them, Represents the shared correction amount, tc = 1 / tg, Represents the optimized shared purchase electricity price, Represents the optimized shared selling electricity price.

[0036] It can be understood that after determining the new partition, determine the electricity demand according to the equipment access situation of the sub-distribution points in the new partition, and then calculate the supply index. The supply index reflects the internal supply and demand situation among the sub-distribution points in the new partition. Then, calculate the shared median price, the shared purchase electricity price, and the shared selling electricity price in sequence. The shared purchase electricity price and the shared selling electricity price can enable the sharing of electricity prices among the new partitions in the distribution network. Furthermore, in the actual application process, mobilize the energy interaction within the new partition, effectively carry out load reduction and load shifting, optimize the energy distribution in the distribution network, and play a guiding role in the coordinated control among the partitions in the distribution network.

[0037] Furthermore, the specific steps of step S80 include: Obtain the official purchase and sale electricity costs and the power exchange cost; ; ; Among them, T Y Represents the official purchase and sale electricity cost, Represents the purchased electricity quantity from the superior power grid, Represents the sold electricity quantity to the superior power grid, D c Represents the moment in a day, T η Represents the power exchange cost, η g Represents the energy storage penalty coefficient, Represents the energy storage charging power, Represents the energy storage discharging power, η k Represents the first adjustable electrical load penalty coefficient, Represents the first adjustable electrical load power, η l Represents the second adjustable electrical load penalty coefficient, Represents the second adjustable electrical load power; With the goal of minimizing the operating cost, construct the second objective function and set the second constraint conditions according to the requirements; ; Among them, T represents the operating cost; The second objective function after considering the shared power purchase price and the shared power selling price is simulated and operated by using a power system model to obtain the second optimal solution and the second best planning scheme of the second objective function.

[0038] It can be understood that the official power purchase and sale costs are calculated through the current market, and the power exchange cost is calculated according to the usage of the equipment in each new sub-region. Then, with the goal of minimizing the operating cost, a second objective function is constructed, and the second constraint condition is set according to the requirements. After considering the shared power purchase price and the shared power selling price calculated in step S70, the operation of the distribution network is simulated by using a power system model to solve the second objective function and obtain the second optimal solution and the second best planning scheme. Because the shared power purchase price and the shared power selling price are considered, the enthusiasm for energy interaction between the new sub-regions in the distribution network can be improved when calculating the power exchange cost, and the advantages of flexible regulation can be exerted, which will affect the power of the two indicators of the first adjustable power load and the second adjustable power load. Therefore, during the process of solving the second objective function, the operating cost is reduced, which provides a guiding role for the coordinated control among the sub-regions in the distribution network. It should be noted that the first adjustable power load is mainly the power load that can be adjusted according to time periods, and the second adjustable power load is mainly the power load that can be appropriately reduced for the stability of the distribution network. The second constraint condition is industry norms, such as the storage constraint of electric energy, the power balance constraint, etc., which will not be elaborated here. The power system model can adopt the commonly used IEEE33 in the industry, and after setting parameters, it can simulate the operation of the distribution network to solve the second objective function.

[0039] In summary, the distribution network multi-substation area distribution control coordination method in the above embodiments of the present invention makes the partitioning of each sub-distribution point more reasonable and the operating cost lower during subsequent use. Please refer to Figure 2 , which shows the distribution network multi-substation area distribution control coordination system in the second embodiment of the present invention, including: The first calculation module 11: constructing a 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; The second calculation module 12: traversing and calculating the second fitness value, the smooth scale, and the intra-region matching scale between one sub-distribution point and the rest of the sub-distribution points, and calculating the comprehensive fitness value by integrating the second fitness value, the smooth scale, and the intra-region matching scale; The first merging module 13: obtaining the fitness difference according to the comprehensive fitness value and the first fitness value, arranging the fitness differences between the rest of the sub-distribution points and the initially selected sub-distribution point from large to small, and selecting the sub-distribution point corresponding to the largest fitness difference to be combined with the initial sub-distribution point to obtain a to-be-merged partition; The second merging module 14: uses the to-be-merged partition as a temporary sub-distribution point, and sequentially absorbs the remaining sub-distribution points into the to-be-merged partition until a preset condition is met; The third merging module 15: sequentially calculates the to-be-merged partition results of the remaining sub-distribution points; The optimization module 16: based on the partitioned to-be-merged partition, calculates the carbon assets and equipment investment prices, and optimizes the to-be-merged partition with the goal of maximizing the carbon asset acquisition rate to obtain a new partition; The third calculation module 17: based on the partitioned new partition, obtains the electricity demand and calculates the supply indicators for each period in the distribution network to calculate the sharing indicator; The fourth calculation module 18: obtains the official power purchase and sale costs and power exchange costs, and obtains the optimal planning scheme with the goal of minimizing the operating cost.

[0040] In the description of this specification, the descriptions referring to terms such as "one embodiment", "some embodiments", "example", "specific example", or "some examples" mean that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0041] The above-described embodiments merely represent several implementation manners of the present invention. The descriptions are relatively specific and detailed, but should not be construed as limiting the scope of the patent of the present invention. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present invention, several modifications and improvements can still be made, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention shall be subject to the appended claims.

Claims

1. A cooperative method for multi-region distribution control of a distribution network, characterized in that, It includes the following steps: Construct the power flow equation of the distribution network, and calculate the first fitness value of each sub-distribution point in the distribution network based on the power flow equation; Traverse and calculate the second fitness value, steady scale, and in-region matching scale between one sub-distribution point and the rest of the sub-distribution points, and calculate the comprehensive fitness value by synthesizing the second fitness value, steady scale, and in-region matching scale; Obtain the fitness difference according to the comprehensive fitness value and the first fitness value, arrange the fitness differences between the rest of the 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 the initial sub-distribution point for combination to obtain the pending merger partition; Take the pending merger partition as a temporary sub-distribution point, and sequentially absorb the remaining sub-distribution points into the pending merger partition until the preset conditions are met; Calculate the pending merger partition results of the remaining sub-distribution points in sequence; Based on the divided pending merger partition, calculate the carbon assets and equipment investment prices, and optimize the pending merger partition with the highest carbon asset acquisition rate as the goal to obtain a new partition; Based on the divided new partition, obtain the power demand and calculate the supply index for each period in the distribution network to calculate the sharing index; Obtain the official power purchase and sale cost and power exchange cost, and obtain the optimal planning scheme with the lowest operating cost as the goal.

2. The multi-substation area distribution control collaboration method for a distribution network according to claim 1, wherein The construction of the power flow equation of the distribution network and the calculation of the first fitness value of each sub-distribution point in the distribution network based on the power flow equation specifically include: Based on the power flow parameters of the distribution network, calculate the voltage magnitude change when all the sub-distribution points are connected with active devices and reactive devices respectively and ; ； Among them, represents the change in the voltage phase angle among the sub-distribution points, represents the change in the voltage amplitude value among the sub-distribution points, represents the change in the active power among the sub-distribution points, represents the change in the reactive power among the sub-distribution points, represents the change in the voltage phase angle when an active device is connected to the sub-distribution point, represents the change in the voltage phase angle when a reactive device is connected to the sub-distribution point; According to the voltage amplitude value change amount at the time of active power and reactive power of the device in the sub-distribution point and , calculate the distance weight of the sub-distribution point; ； ； Among them, represents the influence index of the voltage amplitude value of the sub-distribution point i when the other sub-distribution point j connected to the sub-distribution point i accesses reactive power, represents the influence index of the voltage amplitude value of the other sub-distribution point j connected to the sub-distribution point i when the sub-distribution point i accesses reactive power, represents the influence index of the voltage amplitude value of the sub-distribution point i when the other sub-distribution point j connected to the sub-distribution point i accesses active power, represents the influence index of the voltage amplitude value of the other sub-distribution point j connected to the sub-distribution point i when the sub-distribution point i accesses active power, represents the reactive distance weight between two sub-distribution points, represents the active distance weight between two sub-distribution points; Set the reactive influence factor ω and the active influence factor (1 - ω), and synthesize the active distance weight and the reactive distance weight to obtain the overall weight; ； Normalize the overall weight to obtain the edge weight of the sub-distribution point; Calculate the first fitness value of the sub-distribution point through the edge weight of the sub-distribution point; ； ； Among them, and respectively represent the total weights of sub-distribution points i and sub-distribution point j, and m V represents the total marginal 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.

3. The distribution control coordination method for multiple substations in a distribution network according to claim 2, wherein The traversal calculation of the second fitness value, steady scale, and in-region matching scale between one sub-distribution point and the rest of the sub-distribution points, and the calculation of the comprehensive fitness value by synthesizing the second fitness value, steady scale, and in-region matching scale specifically include: 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 devices among the sub-distribution points i and the remaining sub-distribution points q to obtain the reactive power supply capacity N Q ; Obtain the overlimit value of the voltage to get the minimum reactive power supply quantity N W ; Based on the available reactive power supply amount N Q , the minimum reactive power supply amount N W , and the reactive power distance weight , calculate the stationary scale between sub-distribution points i and sub-distribution point q; ； Calculate the first in-region matching index σ1 and the second in-region matching index σ2 based on the edge weight and the current partition situation; ； ； Among them, C represents the serial number of the partition, G represents the total number of partitions, G C represents partition C, N C represents the total number of sub-distribution points in partition C, B represents the total number of sub-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; Calculate the in-region matching scale σ according to the first in-region matching index and the second in-region matching index; ； Calculate the comprehensive fitness value by synthesizing the second fitness value, steady scale, and in-region matching scale; 。 4. The multi-substation area distribution control and coordination method for a distribution network according to claim 3, wherein The obtaining of the fitness difference according to the comprehensive fitness value and the first fitness value, arranging the fitness differences between the rest of the sub-distribution points and the initially selected sub-distribution point from large to small, and selecting the sub-distribution point corresponding to the largest fitness difference and the initial sub-distribution point for combination to obtain the pending merger partition specifically includes: Subtract the first fitness value of the sub-distribution point i from the comprehensive fitness value between the sub-distribution point i and the rest of the sub-distribution points q to obtain the fitness difference; Arrange the fitness differences between the rest of the sub-distribution points q and the initially selected sub-distribution point i from large to small; Select the sub - distribution point q corresponding to the largest fitness difference and combine it with the initial sub - distribution point i to obtain a to - be - merged partition.

5. The multi-substation area distribution control and coordination method for a distribution network according to claim 4, characterized in that, Taking the to - be - merged partition as a temporary sub - distribution point, successively absorb the remaining sub - distribution points into the to - be - merged partition until a preset condition is met, which specifically includes: Take the to - be - merged partition as the temporary sub - distribution point iq; Calculate the difference between the comprehensive fitness value between the temporary sub - distribution point iq and the remaining unmerged sub - distribution points q and the comprehensive fitness value of the temporary sub - distribution point iq to obtain a fitness difference; And successively absorb all the remaining unmerged sub - distribution points q into the to - be - merged partition until the fitness difference becomes negative.

6. The multi-substation area distribution control and coordination method for a distribution network according to claim 1, wherein Based on the divided to - be - merged partition, calculate the carbon assets and equipment investment price, and optimize the to - be - merged partition with the highest carbon asset acquisition rate as the goal to obtain a new partition, which specifically includes: Calculate carbon assets based on market factors and planning schemes, and calculate the equipment investment price of the to - be - merged partition based on the equipment access situation of the sub - distribution points in the to - be - merged partition; ； ； Among them, C Y represents carbon assets, D represents the number of typical load calculation days in a year, W T represents the instantaneous price of the carbon market, Tθ represents the number of time periods within the typical load calculation days in a year, M is the number of all sub-distribution points of the distribution network, W a,k,t represents the output power of the k-th sub-distribution point in the t-th time period, γ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 trading cost, G l represents the annual network loss cost, G b represents the annual maintenance fee; Based on the carbon assets and the equipment investment price, construct a first objective function with the highest carbon asset acquisition rate as the goal, and set a first constraint condition according to requirements; ； Among them, P C represents the carbon asset acquisition rate; Solve the first objective function through a machine algorithm to obtain the first optimal solution and the first optimal planning scheme of the first objective function; Preset a first optimal solution threshold. If the first optimal solution is less than the first optimal solution threshold, remove the sub - distribution points in the reverse order of the construction of the to - be - merged partition, and circularly calculate the first optimal solution based on the updated to - be - merged partition until the first optimal solution is greater than the first optimal solution threshold; When the first optimal solutions of the first objective function solved based on each to - be - merged partition converge, determine the final to - be - merged partition as the new partition.

7. The multi-substation area distribution control and coordination method for a distribution network according to claim 1, wherein Based on the divided new partition, obtain the electricity demand and calculate the supply index for each period in the distribution network to calculate the sharing index, which specifically includes: Based on the equipment access situation of each sub - distribution point in each new partition, obtain the electricity demand in each new partition; Based on the electricity demand in each new partition, calculate the supply index for each period in the distribution network; ； Among them, represents the electricity sales volume in the Nth new sub-region at time t, represents the electricity purchase volume in the Nth new sub-region at time t, and tg(t) represents the supply index at time t; Based on the current market carbon emission factor, calculate the sharing median price; ； Among them, represents the purchase price of the upper-level power grid, represents the selling price of the upper-level power grid, G P represents the carbon emission index, G Pmax represents the maximum value of the carbon emission index, and ε is a constant, represents the shared median price; Optimize the sharing median price, and calculate the sharing purchase electricity price and sharing selling electricity price based on the optimized sharing median price; ； ； ； Among them, represents the shared correction amount, tc = 1 / tg, represents the optimized shared purchase electricity price, represents the optimized shared selling electricity price.

8. The distribution control collaboration method for multiple substations in a distribution network according to claim 7, wherein Obtain the official purchase and sale electricity costs and power exchange costs, and obtain the optimal planning scheme with the lowest operating cost as the goal, which specifically includes: Obtain the official purchase and sale electricity costs and power exchange costs; ； ； Among them, T Y represents the official electricity purchase and sale cost, represents the electricity purchase quantity from the superior power grid, represents the electricity sale quantity to the superior power grid, D c represents the moment in a day, T η represents the electricity exchange cost, η g represents the energy storage penalty coefficient, represents the energy storage charging power, represents the energy storage discharging power, η k represents the first adjustable electricity load penalty coefficient, represents the first adjustable electricity load power, η l represents the second adjustable electricity load penalty coefficient, represents the second adjustable electricity load power; Construct a second objective function with the lowest operating cost as the goal, and set a second constraint condition according to requirements; ； Where T represents the operating cost; Use the power system model to simulate the operation of the second objective function considering the sharing purchase electricity price and sharing selling electricity price to obtain the second optimal solution and the second best planning scheme of the second objective function.

9. A multi-region distribution control collaborative system for a distribution network, characterized in that Include: The first calculation module: used to construct the power flow equation of the distribution network and calculate the first fitness value of each sub-distribution point in the distribution network based on the power flow equation; The second calculation module: used to traverse and calculate the second fitness value, stationary scale, and in-region matching scale between one sub-distribution point and the rest of the sub-distribution points, and calculate the comprehensive fitness value by synthesizing the second fitness value, stationary scale, and in-region matching scale; The first merging module: used to obtain the fitness difference according to the comprehensive fitness value and the first fitness value, arrange the fitness differences between the rest of the 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 the initial sub-distribution point for combination to obtain a to-be-merged partition; The second merging module: used to take the to-be-merged partition as a temporary sub-distribution point and sequentially absorb the remaining sub-distribution points into the to-be-merged partition until a preset condition is met; The third merging module: used to sequentially calculate the to-be-merged partition results of the remaining sub-distribution points; The optimization module: used to calculate the carbon assets and equipment investment prices based on the divided to-be-merged partition, and optimize the to-be-merged partition with the highest carbon asset acquisition rate as the goal to obtain a new partition; The third calculation module: used to obtain the power demand based on the divided new partition and calculate the supply indicators for each period in the distribution network to calculate the sharing indicators; The fourth calculation module: used to obtain the official power purchase and sale costs and power exchange costs, and obtain the optimal planning scheme with the lowest operating cost as the goal.

Images