A distributed power grid optimization scheduling method based on virtual power plant

Abstract

The present invention relates to the technical field of power grid dispatching, and specifically discloses a distributed power grid optimization dispatching method based on a virtual power plant, comprising the following steps: S1: using the distance of power transmission and the amount of power transmitted as independent variables, and the amount of power lost as a dependent variable, to obtain an amount of power loss function; obtaining a curve of power consumption over time in the power supply area; S2: obtaining the power generation power of the distributed power grid at the starting point of the dispatching cycle, and calculating the power difference; determining the positive power area and the negative power area, and determining the positive power value and the negative power value; S3: obtaining an amount of power loss set in the negative power area; selecting a lost amount of power from the amount of power loss set corresponding to the negative power area according to a preset rule, and calculating the total loss; dispatching power according to the minimum value of the total loss. The present invention can optimize power dispatching and balance power demand and power transmission losses.

Description

Technical Field

[0001] The present invention relates to the technical field of power grid dispatching, and in particular to a distributed power grid optimization dispatching method based on a virtual power plant. Background Art

[0002] A virtual power plant is a power coordination and management system that uses advanced information and communication technologies and software systems to aggregate and coordinate distributed energy resources such as distributed power sources, energy storage systems, controllable loads and electric vehicles, so that it can participate in the electricity market and power grid operation as a special power plant.

[0003] The core concept of distributed power grid is "decentralization". Traditional centralized power generation, such as coal-fired power plants, natural gas power plants, nuclear power plants, large hydropower plants and solar power plants, usually requires long-distance transmission of electricity. Distributed power grids are composed of multiple small power generation devices, which are close to the power consumption points, which can reduce the distance and loss of power transmission.

[0004] The purpose of designing a distributed power grid is to meet the power demand of a specific area or user. In actual situations, the power demand of different areas during the same period may show significant differences. For example, in residential areas and corporate areas, people mainly work in corporate areas during the day and use various electrical appliances, which leads to relatively high power demand in this area. In contrast, in residential areas, there are fewer people at home during the day, so the power demand is relatively low. In this case, the load of the entire distributed power grid can be balanced by transmitting power from areas with lower power demand to areas with higher power demand. However, there will be losses during power transmission, so how to balance the losses and power demand during power transmission has become an urgent problem to be solved. Summary of the invention

[0005] The purpose of the present invention is to provide a distributed power grid optimization scheduling method based on a virtual power plant to solve the above technical problems.

[0006] The purpose of the present invention can be achieved through the following technical solutions:

[0007] A distributed power grid optimization scheduling method based on a virtual power plant comprises the following steps:

[0008] S1: Take the power transmission distance D and the power K transmitted as independent variables, and the power loss Q as the dependent variable to obtain the power loss function f(D, K);

[0009] Set the collection area and obtain the curve of the power consumption of each power supply area in the collection area over time f '(t), the power supply area is the area where power is supplied through a distributed power grid;

[0010] S2: Set the scheduling period [T sta , T end ], and obtain the power generation power P of the distributed power grid at the starting point of the scheduling cycle, and calculate the power difference, the calculation formula is:

[0011] ;

[0012] According to the electric charge difference, the positive charge area and the negative charge area are determined, and the positive charge value K is determined. + and negative charge K - , the positive power area is the power supply area that needs to deliver electricity to the outside world during the current scheduling cycle, and the positive power value is the amount of electricity that the positive power area can deliver to the outside world. The negative power area is the power supply area that needs to deliver electricity to the outside world during the current scheduling cycle, and the negative power value is the amount of electricity that the negative power area needs to deliver to the outside world;

[0013] S3: Starting from the i-th negatively charged region, obtain the distance between the i-th negatively charged region and the positively charged region, and substitute the negative values ​​corresponding to the i-th negatively charged region into the power loss function to obtain the power loss set Qi of the i-th negatively charged region = (Q i1 , Q i2 , …, Q im ), where Q im It represents the power loss when the mth positively charged region transmits power to the ith negatively charged region, where m is the total number of positively charged regions;

[0014] Selecting a power loss from the power loss set corresponding to the negative power area according to a preset rule, and calculating the sum of the selected power losses as the total loss, wherein the preset rule is that the power loss corresponding to the same positive power area cannot be selected;

[0015] The minimum value of the total loss is obtained, and the element in each power loss set for which the minimum value is calculated is determined as a scheduling element, and within the current scheduling cycle, power is scheduled according to the negative power area and the positive power area corresponding to the scheduling element.

[0016] As a further solution of the present invention: in step S1, the process of obtaining the power loss function specifically includes:

[0017] Setting a distance range for power transmission as the distance range, and selecting a number of distance nodes within the distance range at preset distance intervals;

[0018] Setting a range of the amount of electricity to be transmitted as an amount of electricity range, and selecting a number of electricity nodes within the amount of electricity range at preset intervals of electricity;

[0019] While keeping the power transmission distance at the ath distance node, the power corresponding to each power node is transmitted respectively, and the power loss is recorded respectively to generate data points (D a , K b , Q ab ), where D a Indicates the distance corresponding to the a-th distance node, K b Indicates the power corresponding to the bth power node, Q ab Indicates the power loss corresponding to the a-th distance node and the b-th power node;

[0020] The data points are fitted through nonlinear regression analysis to obtain the power loss function.

[0021] As a further solution of the present invention: in the step S1, the process of obtaining the curve of the power consumption of each power supply area in the acquisition area over time specifically includes:

[0022] Set m collection cycles, where m is a preset value. Set several time nodes at preset time intervals within the collection cycle, obtain the average power consumption of the power supply area at the same time node, as the power consumption at the time node, and fit the average power consumption at all time nodes by the least squares method to obtain a curve of power consumption changing with time.

[0023] As a further solution of the present invention: in the step S3, the positive charge area and the negative charge area are determined according to the charge difference, and the positive charge value K is determined. + and negative charge K - The process specifically includes:

[0024] When the power difference is greater than or equal to 0, the corresponding power supply area is regarded as a positive power area, and the positive power value K is calculated. + =K'-K ret +∆K, where K' is the current amount of electricity stored in the energy storage device in the power supply area, K ret The preset guaranteed power is the power that should be retained without affecting the normal use in the power supply area;

[0025] When the power difference is less than 0, calculate the undetermined value K tbd =K'+∆K;

[0026] When the undetermined value is less than 0, the corresponding power supply area is regarded as a negative power area, and the undetermined value is regarded as a negative power value K. - ;

[0027] When the undetermined value is greater than or equal to 0, calculate the judgment value K jud =K'-K ret+∆K, and perform the following steps:

[0028] When the determination value is less than or equal to 0, the power supply area does not participate in subsequent steps in the current scheduling cycle;

[0029] When the judgment value is greater than 0, the power supply area is regarded as a positive power area, and the judgment value is regarded as a positive power value.

[0030] As a further solution of the present invention: in the step S3, the positive charge area and the negative charge area are determined according to the charge difference, and the positive charge value K is determined. + and negative charge K - The following steps are also performed during the process:

[0031] When the sum of the positive and negative electrical values ​​is less than 0, an early warning message is sent to inform the staff;

[0032] When the sum of the positive electrical value and the negative electrical value is greater than or equal to 0, the subsequent steps are performed normally.

[0033] As a further solution of the present invention: in the step S3, when the negative charge value corresponding to the i-th negative charge region is greater than the positive charge value of the j-th positive charge region, the integration step is performed:

[0034] Calculate the missing power K qs =K +j +K -i , where K +j represents the positive charge value of the jth positively charged region, K -i represents the negative charge value of the i-th negative charge region;

[0035] Obtain the distance between the positively charged area except the jth positively charged area and the ith negatively charged area, and substitute the distance and the absolute value of the missing power into the power loss function to obtain the power loss when the jth positively charged area and the remaining positively charged areas jointly supply power to the ith negatively charged area, and record it as Q ij , and update the positive value of the positive region c that supplies power to the i-th negative region together with the j-th positive region to K +c -K qs。

[0036] As a further solution of the present invention: when there is a positively charged region with a positive charge value less than the missing charge, the integration step is performed until the charge delivered by the positively charged region to the i-th negatively charged region reaches the negative charge value of the i-th negatively charged region.

[0037] As a further solution of the present invention: in the step S3, when the negative charge value corresponding to the i-th negative charge region is less than the positive charge value of the j-th positive charge region, the following steps are performed:

[0038] Calculate the excess power K dy =K +j +K -i , and split the jth positively charged region into two positively charged regions J1 and J2, the positive charge value of J1 is K -i The absolute value of J1 is the excess charge K. dy , and perform subsequent steps, and the two positively charged regions J1 and J2 are at the same position.

[0039] Beneficial effects of the present invention: In the present invention, firstly, the power loss function and the curve of power consumption changing with time are obtained, which is the basis for the subsequent balancing of power demand and loss in power transmission process; then the dispatching cycle is set, and the power difference is calculated according to the power generation power of the distributed power grid at the starting point of the dispatching cycle; it can be understood that the starting point of the dispatching cycle is the current time, and the reason for using the power generation power of the distributed power grid at the starting point of the dispatching cycle is that the power generation power at the current time reflects the current power production status, and provides initial conditions and reference benchmarks for power dispatch in the next period of time (i.e., the dispatching cycle); then, the positive power area and the negative power area, and determine the positive power value and the negative power value; it can be understood that this is In order to determine which power supply areas can supply power to the outside world and which power supply areas need external power supply to ensure normal use; afterwards, obtain the power loss set, and select a lost power according to the preset rules to calculate the total loss, and determine the scheduling element according to the minimum value of the total loss, that is, which negative power areas should receive power from which positive power areas, so as to perform effective power scheduling in the current scheduling cycle; it is worth noting that after the positive power area supplies power to a certain negative power area, the remaining power available for other negative power areas will be reduced accordingly. Therefore, the same positive power area cannot be selected to supply different negative power areas multiple times to ensure the accuracy of calculating the loss and avoid over-scheduling, which affects the normal use of the positive power area. The present invention can optimize power scheduling and balance the power demand and the loss in the power transmission process. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] The present invention will be further described below in conjunction with the accompanying drawings.

[0041] Figure 1 It is a flow chart of a distributed power grid optimization scheduling method based on a virtual power plant in the present invention. DETAILED DESCRIPTION

[0042] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0043] See also Figure 1 As shown, the present invention is a distributed power grid optimization scheduling method based on a virtual power plant, comprising the following steps:

[0044] S1: Take the power transmission distance D and the power K transmitted as independent variables, and the power loss Q as the dependent variable to obtain the power loss function f(D, K);

[0045] Set the collection area and obtain the curve of the power consumption of each power supply area in the collection area over time f '(t), the power supply area is the area where power is supplied through a distributed power grid;

[0046] S2: Set the scheduling period [T sta , T end ], and obtain the power generation power P of the distributed power grid at the starting point of the scheduling cycle, and calculate the power difference, the calculation formula is:

[0047] ;

[0048] According to the electric charge difference, the positive charge area and the negative charge area are determined, and the positive charge value K is determined. + and negative charge K - , the positive power area is the power supply area that needs to deliver electricity to the outside world during the current scheduling cycle, and the positive power value is the amount of electricity that the positive power area can deliver to the outside world. The negative power area is the power supply area that needs to deliver electricity to the outside world during the current scheduling cycle, and the negative power value is the amount of electricity that the negative power area needs to deliver to the outside world;

[0049] S3: Starting from the i-th negatively charged region, obtain the distance between the i-th negatively charged region and the positively charged region, and substitute the negative values ​​corresponding to the i-th negatively charged region into the power loss function to obtain the power loss set Qi of the i-th negatively charged region = (Q i1 , Q i2 , …, Q im ), where Q im It represents the power loss when the mth positively charged region transmits power to the ith negatively charged region, where m is the total number of positively charged regions;

[0050] Selecting a power loss from the power loss set corresponding to the negative power area according to a preset rule, and calculating the sum of the selected power losses as the total loss, wherein the preset rule is that the power loss corresponding to the same positive power area cannot be selected;

[0051] The minimum value of the total loss is obtained, and the element in each power loss set for which the minimum value is calculated is determined as a scheduling element, and within the current scheduling cycle, power is scheduled according to the negative power area and the positive power area corresponding to the scheduling element.

[0052] It should be noted that first, the power loss function and the curve of power consumption changing with time are obtained, which is the basis for the subsequent balance of power demand and loss in the power transmission process; then the dispatch cycle is set, and the power difference is calculated based on the power generation of the distributed power grid at the starting point of the dispatch cycle; it can be understood that the starting point of the dispatch cycle is the current time. The reason for using the power generation of the distributed power grid at the starting point of the dispatch cycle is that the power generation at the current time reflects the current power production status, which provides initial conditions and reference benchmarks for power dispatch in the next period of time (i.e., the dispatch cycle); then, the positive power area and the negative power area, and determine the positive power value and the negative power value; it can be understood that this is to determine Which power supply areas can supply power to the outside world, and which power supply areas need external power supply to ensure normal use; then, obtain the power loss set, and select a lost power according to the preset rules to calculate the total loss, and determine the scheduling element based on the minimum value of the total loss, that is, which negative power areas should receive power from which positive power areas, so as to carry out effective power scheduling within the current scheduling cycle; it is worth noting that after the positive power area supplies power to a negative power area, the remaining power available for other negative power areas will be reduced accordingly. Therefore, the same positive power area cannot be selected multiple times to supply different negative power areas to ensure the accuracy of calculating losses and avoid over-scheduling, which affects the normal use of the positive power areas.

[0053] In another preferred implementation of the present invention, in step S1, the process of obtaining the power loss function specifically includes:

[0054] Setting a distance range for power transmission as the distance range, and selecting a number of distance nodes within the distance range at preset distance intervals;

[0055] Setting a range of the amount of electricity to be transmitted as an amount of electricity range, and selecting a number of electricity nodes within the amount of electricity range at preset intervals of electricity;

[0056] While keeping the power transmission distance at the ath distance node, the power corresponding to each power node is transmitted respectively, and the power loss is recorded respectively to generate data points (D a , K b , Q ab ), where D a Indicates the distance corresponding to the a-th distance node, K b Indicates the power corresponding to the bth power node, Q abIndicates the power loss corresponding to the a-th distance node and the b-th power node;

[0057] The data points are fitted through nonlinear regression analysis to obtain the power loss function.

[0058] It is worth noting that the distance range and power range are set to reduce the amount of experimental data required and save computing resources.

[0059] In another preferred implementation of the present invention, in step S1, the process of obtaining the curve of power consumption of each power supply area in the acquisition area over time specifically includes:

[0060] Set m collection cycles, where m is a preset value. Set several time nodes at preset time intervals within the collection cycle, obtain the average power consumption of the power supply area at the same time node, as the power consumption at the time node, and fit the average power consumption at all time nodes by the least squares method to obtain a curve of power consumption changing with time.

[0061] In another preferred implementation of the present invention, in the step S3, the positive charge area and the negative charge area are determined according to the charge difference, and the positive charge value K is determined. + and negative charge K - The process specifically includes:

[0062] When the power difference is greater than or equal to 0, the corresponding power supply area is regarded as a positive power area, and the positive power value K is calculated. + =K'-K ret +∆K, where K' is the current amount of electricity stored in the energy storage device in the power supply area, K ret The preset guaranteed power is the power that should be retained without affecting the normal use in the power supply area;

[0063] When the power difference is less than 0, calculate the undetermined value K tbd =K'+∆K;

[0064] When the undetermined value is less than 0, the corresponding power supply area is regarded as a negative power area, and the undetermined value is regarded as a negative power value K. - ;

[0065] When the undetermined value is greater than or equal to 0, calculate the judgment value K jud =K'-K ret +∆K, and perform the following steps:

[0066] When the determination value is less than or equal to 0, the power supply area does not participate in subsequent steps in the current scheduling cycle;

[0067] When the judgment value is greater than 0, the power supply area is regarded as a positive power area, and the judgment value is regarded as a positive power value.

[0068] In another preferred implementation of the present invention, in the step S3, the positive charge area and the negative charge area are determined according to the charge difference, and the positive charge value K is determined. + and negative charge K - The following steps are also performed during the process:

[0069] When the sum of the positive and negative electrical values ​​is less than 0, an early warning message is sent to inform the staff;

[0070] When the sum of the positive electrical value and the negative electrical value is greater than or equal to 0, the subsequent steps are performed normally.

[0071] It should be noted that the preset guaranteed power is to ensure that each power supply area has sufficient power to guarantee basic operation and to avoid affecting normal power consumption in the power supply area due to excessive output.

[0072] In another preferred implementation of the present invention, in step S3, when the negative charge value corresponding to the i-th negative charge region is greater than the positive charge value of the j-th positive charge region, the integration step is performed:

[0073] Calculate the missing power K qs =K +j +K -i , where K +j represents the positive charge value of the jth positively charged region, K -i represents the negative charge value of the i-th negative charge region;

[0074] Obtain the distance between the positively charged area except the jth positively charged area and the ith negatively charged area, and substitute the distance and the absolute value of the missing power into the power loss function to obtain the power loss when the jth positively charged area and the remaining positively charged areas jointly supply power to the ith negatively charged area, and record it as Q ij , and update the positive value of the positive region c that supplies power to the i-th negative region together with the j-th positive region to K +c -K qs。

[0075] It is understandable that when the demand of a negatively charged area exceeds the amount of electricity that a single positively charged area can provide, power needs to be obtained from multiple positively charged areas. By calculating the missing electricity, it is possible to determine how much electricity is needed in addition to the jth positively charged area to meet the demand of the ith negatively charged area, thereby ensuring that the demand of the negatively charged areas is met while avoiding excessive output from the positively charged areas.

[0076] In another preferred implementation of the present invention, when there is a positively charged region whose positive charge value is less than the missing charge, the integration step is performed until the charge delivered by the positively charged region to the i-th negatively charged region reaches the negative charge value of the i-th negatively charged region.

[0077] In another preferred implementation of the present invention, in step S3, when the negative charge value corresponding to the i-th negative charge region is less than the positive charge value of the j-th positive charge region, the following steps are performed:

[0078] Calculate the excess power K dy =K +j +K -i , and split the jth positively charged region into two positively charged regions J1 and J2, the positive charge value of J1 is K -i The absolute value of J1 is the excess charge K. dy , and perform subsequent steps, and the two positively charged regions J1 and J2 are at the same position.

[0079] It is worth noting that in actual situations, the positive charge value of a positively charged area may simultaneously meet the negative charge values ​​of multiple negatively charged areas, and when the positively charged area supplies power to the negatively charged areas, the loss may also be relatively small. By splitting the positively charged areas, the excess power in the positively charged areas can be fully utilized.

[0080] The above is a detailed description of an embodiment of the present invention, but the content is only a preferred embodiment of the present invention and cannot be considered to limit the scope of implementation of the present invention. All equivalent changes and improvements made within the scope of the present invention should still fall within the scope of the patent coverage of the present invention.

Claims

1. A distributed power grid optimization scheduling method based on virtual power plant, characterized in that: The following steps are involved: S1: Take the power transmission distance D and the power K transmitted as independent variables, and the power loss Q as the dependent variable to obtain the power loss function f(D, K); Setting a collection area, and obtaining a curve f'(t) of the power consumption of each power supply area in the collection area over time, wherein the power supply area is an area powered by a distributed power grid; S2: Set the dispatching period [Tsta, Tend], obtain the power generation power P of the distributed power grid at the starting point of the dispatching period, and calculate the power difference, the calculation formula is: According to the electric charge difference, the positive charge area and the negative charge area are determined, and the positive charge value K is determined. + and negative power value K-, the positive power area is the power supply area that needs to deliver power to the outside world in the current scheduling cycle, the positive power value is the power that the positive power area can deliver to the outside world, the negative power area is the power supply area that needs to deliver power to the outside world in the current scheduling cycle, and the negative power value is the power that the negative power area needs to deliver to the outside world; S3: Starting from the ith negatively charged region, obtain the distance between the ith negatively charged region and the positively charged region, and substitute the negative power values ​​corresponding to the ith negatively charged region into the power loss function to obtain the power loss set Qi=(Qi1, Qi2, ..., Qim) of the ith negatively charged region, where Qim represents the power loss when the mth positively charged region transmits power to the ith negatively charged region, and m is the total number of positively charged regions; Selecting a power loss from the power loss set corresponding to the negative power area according to a preset rule, and calculating the sum of the selected power losses as the total loss, wherein the preset rule is that the power loss corresponding to the same positive power area cannot be selected; Obtaining the minimum value of the total loss, and determining the element in each power loss set for which the minimum value is calculated as the scheduling element, and scheduling power according to the negative power area and the positive power area corresponding to the scheduling element in the current scheduling period; According to the electric charge difference, the positive charge area and the negative charge area are determined, and the positive charge value K is determined. + The process of generating a negative charge value K- specifically includes: When the power difference is greater than or equal to 0, the corresponding power supply area is regarded as a positive power area, and the positive power value K is calculated. + =K'-Kret+ΔK, where K' is the current power of the energy storage device in the power supply area, and Kret is the preset guaranteed power, which is the power that should be retained without affecting the normal use in the power supply area; When the power difference is less than 0, the undetermined value Ktbd=K'+ΔK is calculated; When the undetermined value is less than 0, the corresponding power supply area is regarded as a negative power area, and the undetermined value is regarded as a negative power value K-; When the undetermined value is greater than or equal to 0, the judgment value Kjud=K'-Kret+ΔK is calculated, and the following steps are performed: When the determination value is less than or equal to 0, the power supply area does not participate in subsequent steps in the current scheduling cycle; When the judgment value is greater than 0, the power supply area is regarded as a positive power area, and the judgment value is regarded as a positive power value.

2. A distributed power grid optimization scheduling method based on a virtual power plant according to claim 1, characterized in that: In step S1, the process of obtaining the power loss function specifically includes: Setting a distance range for power transmission as the distance range, and selecting a number of distance nodes within the distance range at preset distance intervals; Setting a range of the amount of electricity to be transmitted as an amount of electricity range, and selecting a number of electricity nodes within the amount of electricity range at preset intervals of electricity; While keeping the power transmission distance at the ath distance node, transmit the power corresponding to each power node respectively, and record the power loss respectively, and generate data points (Da, Kb, Qab), where Da represents the distance corresponding to the ath distance node, Kb represents the power corresponding to the bth power node, and Qab represents the power loss corresponding to the ath distance node and the bth power node; The data points are fitted through nonlinear regression analysis to obtain the power loss function.

3. The distributed power grid optimization scheduling method based on virtual power plant according to claim 1 is characterized in that: In the step S1, the process of obtaining the curve of the power consumption of each power supply area in the acquisition area over time specifically includes: Set m collection cycles, where m is a preset value. Set several time nodes at preset time intervals within the collection cycle, obtain the average power consumption of the power supply area at the same time node, as the power consumption at the time node, and fit the average power consumption at all time nodes by the least squares method to obtain a curve of power consumption changing with time.

4. The distributed power grid optimization scheduling method based on virtual power plant according to claim 1 is characterized in that: In the step S3, the positive charge area and the negative charge area are determined according to the charge difference, and the positive charge value K is determined. + And negative value K- process also performs the following steps: When the sum of the positive and negative electrical values ​​is less than 0, an early warning message is sent to inform the staff; When the sum of the positive electrical value and the negative electrical value is greater than or equal to 0, the subsequent steps are performed normally.

5. The distributed power grid optimization scheduling method based on virtual power plant according to claim 1 is characterized in that: In the step S3, when the negative charge value corresponding to the i-th negative charge region is greater than the positive charge value of the j-th positive charge region, the integration step is performed: Calculate the missing power Kqs = K + j+Ki, where K + j represents the positive charge value of the jth positively charged region, and Ki represents the negative charge value of the ith negatively charged region; Obtain the distance between the positive areas except the jth positive area and the ith negative area, and substitute the distance and the absolute value of the missing electricity into the electricity loss function to obtain the electricity loss when the jth positive area and the remaining positive areas jointly supply power to the ith negative area, and record it as Qij, and update the positive value of the positive area c that supplies power to the ith negative area together with the jth positive area to K + c-Kqs.

6. A distributed power grid optimization scheduling method based on virtual power plant according to claim 5, characterized in that: When there is a positively charged region whose positive charge value is less than the missing charge, the integration step is performed until the charge delivered by the positively charged region to the i-th negatively charged region reaches the negative charge value of the i-th negatively charged region.

7. A distributed power grid optimization scheduling method based on virtual power plant according to claim 5, characterized in that: In the step S3, when the negative charge value corresponding to the i-th negative charge region is smaller than the positive charge value of the j-th positive charge region, the following steps are performed: Calculate the excess power Kdy = K + j+Ki, and split the j-th positively charged area into two positively charged areas J1 and J2, the positive charge value of J1 is the absolute value of Ki, and power is supplied to the i-th negatively charged area through J1, the positive charge value of J2 is the excess power Kdy, and the two positively charged areas J1 and J2 are in the same position, and the subsequent steps are performed.

Images