CN113629766B - Automatic voltage control method and device for reducing voltage fluctuation of new energy collection area - Google Patents
Automatic voltage control method and device for reducing voltage fluctuation of new energy collection area Download PDFInfo
- Publication number
- CN113629766B CN113629766B CN202110925892.6A CN202110925892A CN113629766B CN 113629766 B CN113629766 B CN 113629766B CN 202110925892 A CN202110925892 A CN 202110925892A CN 113629766 B CN113629766 B CN 113629766B
- Authority
- CN
- China
- Prior art keywords
- voltage
- new energy
- bus
- period
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000010248 power generation Methods 0.000 claims abstract description 51
- 238000004364 calculation method Methods 0.000 claims abstract description 45
- 230000035945 sensitivity Effects 0.000 claims abstract description 29
- 238000002347 injection Methods 0.000 claims abstract description 6
- 239000007924 injection Substances 0.000 claims abstract description 6
- 230000008859 change Effects 0.000 claims description 20
- 230000007423 decrease Effects 0.000 claims description 12
- 230000005283 ground state Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 6
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims 1
- 230000001105 regulatory effect Effects 0.000 abstract description 3
- 238000011217 control strategy Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 230000002829 reductive effect Effects 0.000 description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 4
- 239000003990 capacitor Substances 0.000 description 4
- 229910052799 carbon Inorganic materials 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000033228 biological regulation Effects 0.000 description 3
- 230000003449 preventive effect Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/16—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The disclosure relates to an automatic voltage control method and device for reducing voltage fluctuation of a new energy collection area, and belongs to the technical field of automatic voltage control of power systems. The method comprises the following steps: temporarily acquiring a base state power flow calculation result of a new energy collection area in a power grid in each automatic voltage control period, and calculating voltage sensitivity of injection active power of a high-voltage side bus of each new energy power station in the new energy collection area to each bus; calculating the predicted value of each busbar voltage in the new energy collecting area at each moment in the predicted period by using the sensitivity; and calculating the upper limit of the voltage safety area and the lower limit of the voltage safety area of each bus in the prediction period according to the voltage prediction value, so as to be used for automatic voltage control of the prediction period. The method and the device are suitable for the concentrated grid-connected new energy power station collecting regional power grid, reactive power equipment of the power grid can be regulated prophylactically before the new energy power generation is changed rapidly, and the safety and stability of the voltage of the new energy collecting regional power grid are guaranteed.
Description
Technical Field
The disclosure belongs to the technical field of automatic voltage control of power systems, and particularly relates to an automatic voltage control method and device for reducing voltage fluctuation of a new energy collection area.
Background
An automatic voltage control (AVC, automatic Voltage Control) system is an important means for realizing safe (voltage stability margin improvement), economical (network loss reduction) and high-quality (voltage qualification rate improvement) operation of a power transmission network. The AVC system is constructed on a power grid Energy Management System (EMS), can utilize real-time operation data of a power grid, scientifically decides an optimal reactive voltage adjustment scheme from the perspective of global optimization of the power grid, and automatically distributes the optimal reactive voltage adjustment scheme to a power plant, a transformer substation and a lower power grid dispatching mechanism for execution. Sun Hong, zhang Baming, guo Qinglai describes the architecture of large grid automatic voltage control in soft partition based global voltage optimization control system design (power system automation, 2003, volume 27, 8 th, pages 16-20).
The master station part of the AVC system is realized based on software in a power system control center, and the voltage control strategy of the AVC system to the power transmission network mainly comprises reactive power control strategies to all generators of a power plant and reactive power equipment control strategies of a transformer substation in class 2. The reactive power control strategy for each generator of the power plant adopts the following main modes at present: after reactive power adjustment quantity of each unit of the power plant is obtained through reactive power optimization calculation by the AVC master station system of the dispatching center, the reactive power adjustment quantity is sent to an AVC substation system of the power plant through a data communication channel, and after the AVC substation of the power plant receives the reactive power adjustment quantity of the generator, the reactive power generated by the generator is adjusted in a stepping mode according to the running state of each generator in the current power plant until the reactive power reaches the adjustment quantity issued by the AVC master station. The reactive power equipment control strategy of the transformer substation is a switching instruction of reactive power compensation equipment, wherein the reactive power equipment mainly comprises a capacitor and a reactor, and when the capacitor is switched in or the reactor is cut off, the bus voltage is increased; when the capacitor is cut off or the reactor is put in, the bus voltage decreases. The AVC master station issues an instruction for switching in or switching off the reactive equipment, and an automatic monitoring system in the transformer substation finds a breaker connected with the reactive equipment according to the received instruction and switches on or off the breaker so as to complete switching in or switching off of the reactive equipment.
In order to achieve the strategic goals of national carbon peak in 2030 and carbon neutralization in 2060, the construction of a novel power system taking new energy as a main body is a necessary way for the power industry to promote self carbon emission reduction and support carbon emission reduction of the whole society. New energy has become the preferred scheme of the clean energy replacement of the regional propulsion of the lack of conditions, formed new energy and gathered the district, its grid-connected capacity increases fast, bring new challenges for the dispatch operation of electric wire netting. On the one hand, in some areas closer to the load center, such as offshore areas in coastal developed cities, and with wind all the year round, the active requirements of the electric load center are very suitable; on the other hand, a large-scale centralized development mode is adopted, reactive voltage support of conventional water and fire power plant units is lacking in a new energy centralized grid-connected area, the short-circuit capacity of the system is small, and the inherent intermittent change of active power generation of the new energy can cause larger voltage fluctuation, so that larger difficulty is brought to voltage regulation.
In a new energy power generation area which is concentrated and connected into a power grid in a large scale, the new energy power station cascading off-grid fault caused by overlarge voltage fluctuation easily occurs, and the safe and stable operation of the whole power grid is further affected. In order to solve the problem, on one hand, the reactive power regulation capability of the new energy station needs to be fully utilized to provide voltage support for new energy power generation; on the other hand, the method also reasonably regulates and controls various reactive resources of the new energy station and the sending-out channel to realize coordination control.
Automatic Voltage Control (AVC) has been widely used in various levels of grid dispatching centers. In recent years, voltage control of an access power grid is concentrated around large-scale new energy, and many research results are available. Guo Qinglai, wang Bin and Sun Hong propose an autonomous collaborative voltage control technology system for supporting centralized access of large-scale wind power in an autonomous collaborative voltage control technology (automatization of a power system, 2015, volume 39, 1 st phase, pages 88-93) for supporting centralized access of large-scale wind power, and realize autonomous control at a new energy power station level, and utilize a new energy power station substation to coordinately control regulating equipment with different time constants such as a static reactive compensator, a static reactive generator, a wind turbine generator, a capacitor reactor and the like so as to inhibit voltage fluctuation induced by intermittent wind power output; the cooperative control is realized at the system level, voltage fluctuation is reduced under normal conditions through agile secondary voltage control which can be self-adapted to wind power change, and when the off-grid risk is large, the normal and safe operation state of the collecting area is ensured by utilizing the preventive control based on the Safety Constraint Optimal Power Flow (SCOPF).
The automatic voltage control system (AVC) running in each stage of grid dispatching center is mostly controlled by periodic control (usually 5 minutes) and calculation for data of one data section of the grid. In the period of high active power of new energy, a large amount of active power change is likely to be generated in the interval time of two controls, so that the AVC system cannot respond timely, the bus voltage is caused to generate larger fluctuation, and even short-time out-of-limit conditions are caused.
Because the large-scale new energy collection area has stronger intermittence, the voltage fluctuation of the grid-connected area caused by active power change or fault disturbance is more severe, and the linkage off-grid of the new energy station in the area can be possibly caused. Therefore, the voltage safety operation domain of the new energy collection area needs to be determined through online voltage safety domain evaluation so as to meet the requirement that the voltage safety can be ensured after the new energy power generation trend changes in the future and the risk faults occur. Compared with the traditional AVC, the existing partial AVC is added with a new energy collection area voltage safety evaluation module, the module considers the power generation trend change of the new energy collection area in the future for a period of time, predicts the future voltage change trend according to the active power generation trend change, calculates the voltage static safety domain of the large-scale new energy collection area, and adjusts the bus voltage limit range according to the voltage change trend prediction. The adjusted limit value range is generally a subset of the planned limit value range which is originally set manually and is used as constraint conditions of three-level optimization and two-level voltage control of AVC, so that voltage safety prevention control considering future voltage change trend and risk condition of a large-scale new energy collection area is realized.
Therefore, in order to reduce the larger fluctuation of the bus voltage of the power grid caused by the inherent intermittent change of active power generation of the new energy unit, when automatic voltage control is performed, the change amplitude of the bus voltage of the near zone, which possibly occurs due to the change of active power, needs to be predicted according to the active power prediction of the new energy in the current new energy collecting zone, then the upper limit of the voltage safety domain and the lower limit of the voltage safety zone of each bus can be further calculated, and the voltage prevention control is performed on the power grid of the near zone according to the requirement, so that the fluctuation of the bus voltage entering the new energy collecting zone is reduced, and the condition that the power grid of the new energy collecting zone is out of limit in voltage is avoided.
Disclosure of Invention
The purpose of the present disclosure is to overcome the shortcomings of the prior art, and to provide an automatic voltage control method and device for reducing voltage fluctuation in a new energy collection area. The method and the device are suitable for the concentrated grid-connected new energy power station collecting regional power grid, reactive power equipment of the power grid can be regulated prophylactically before the new energy power generation is changed rapidly, and the safety and stability of the voltage of the new energy collecting regional power grid are guaranteed.
An embodiment of a first aspect of the present disclosure provides an automatic voltage control method for reducing voltage fluctuation in a new energy collection area, including:
Temporarily acquiring a base state power flow calculation result of a new energy collection area in the power grid in each automatic voltage control period;
calculating the voltage sensitivity of the injection active power of the high-voltage side bus of each new energy power station in the new energy collection area to each bus by using the ground state power flow calculation result;
calculating the predicted value of each busbar voltage in the new energy collecting area at each moment in the predicted period by using the sensitivity;
and calculating the upper limit of the voltage safety area and the lower limit of the voltage safety area of each bus in the prediction period according to the predicted value of each bus voltage, so as to be used for automatic voltage control of the prediction period.
In one embodiment of the present disclosure, the base state power flow calculation result of the new energy collection area in the power grid includes:
active power of each new energy power station and each bus voltage value in the new energy collection area;
wherein the active power of each new energy power station forms an active power set U of each new energy power station wstp ＝(P 1 ，P 2 ，…，P I Each busbar voltage value forms a busbar voltage set U wu ＝{U 1 ,U 2 ,…,U J -a }; i is the total number of new energy power stations contained in the new energy collecting area, and J is the total number of buses contained in the new energy collecting area.
In one embodiment of the present disclosure, calculating, using the sensitivity, each bus voltage predicted value in the new energy collection area at each time in a predicted period includes:
1) Let the prediction period be t 0 To t 0 +T 1 Wherein t is 0 Representing the starting time of the predicted period, T 1 Is the length of the prediction period;
2) Let the current prediction time be x, and the initial value of x be t 0 ；
3) For each new energy power station i in the new energy collection area, determining an active power generation predicted value P of the new energy power station i at the moment x i,x The method is characterized by comprising the following steps:
if x coincides with the corresponding time of any current day active power generation prediction initial value of the new energy power station i, taking the current day active power generation prediction initial value of the new energy power station i corresponding to the x time as the active power generation prediction value of the new energy power station i at the x time; otherwise, linear interpolation of the current day active power generation prediction initial values of two adjacent new energy power stations i before and after the x moment is adopted as the active power generation prediction value of the new energy power station i at the x moment, and the expression is as follows:
wherein, the current day active power generation prediction initial value D of the new energy power station i i,n Obtained from a grid energy management system, n=1..n, N is new energy electricityThe sequence number of the station active power generation prediction initial value, and v is used as a time interval between adjacent prediction initial values; x/v represents the integer result of dividing x by v, and x% v represents the remainder result of dividing x by v;
4) Utilizing the result of the step 3) to establish an active power generation predictive value vector F of all new energy power stations at the moment x x ＝[P 1,x P 2,x …P I,x ]；
According to set U wstp Form the corresponding vector f= [ P ] 1 P 2 …P I ]Calculating a change value vector delta F of active power of each new energy power station at the moment x x ＝F x -F＝[P 1,x -P 1 P 2,x -P 2 …P I,x -P I ]；
5) Calculating voltage change value vector delta U of each bus at x time x ＝ΔF x *S cv ＝[ΔU 1,x ΔU 2,x …ΔU J,x ]According to set U wu ＝{U 1 ,U 2 ,…,U J Voltage vector U formed 0 ＝[U 1 U 2 …U J ]The vector formed by the predicted values of the voltages of all the buses at the moment x is obtained to be U x ＝U 0 +ΔU x Wherein U is x ＝[U 1 +ΔU 1,x U 2 +ΔU 2,x …U J +ΔU J,x ]＝[U 1,x U 2,x …U J,x ]；
Wherein S is cv For the sensitivity matrix:
wherein S is ij The voltage sensitivity of the j-th bus is injected on the high-voltage side bus of the i-th new energy power station;
6) Let x=x+1, then return to step 3); until the current prediction time is t 0 +T 1 Obtaining the predicted values of the voltages of all buses in the new energy collecting area at all times of the predicted period, and constructing the predicted values at (t) 0 ,t 0 +T 1 ) New time periodEach busbar voltage predicted value in the energy collection area is assembled:
in one embodiment of the present disclosure, the calculating, according to the predicted value of each bus voltage, an upper voltage safety region limit and a lower voltage safety region limit of each bus in the predicted period for automatic voltage control of the predicted period includes:
1) Let t be the prediction period (t 0 ,t 0 +T 1 ) Any time in the process;
for each busbar j, an initial value DeltaU is set t.inc,j ＝0，ΔU t.dec,j =0; wherein DeltaU t.inc,j Represents the maximum amplitude of the continuous increase in voltage of bus j from time t over the predicted period, deltaU t.dec,j Representing the maximum amplitude of continuous decrease of the voltage of the busbar j in the prediction period from the moment t;
2) From the time x=t+1, Δu is calculated t.inc,j The specific steps are as follows:
2-1) from U j ＝{U j,x ,x＝t 0 ,...,t 0 +T 1 Respectively obtaining voltage predicted values U of buses j at x and x-1 moments j,x 、U j,x-1 Calculating the voltage variation delta U between the two predicted values j,x ＝U j,x -U j,x-1 The method comprises the steps of carrying out a first treatment on the surface of the Wherein U is j Representative bus j is in the prediction period (t 0 ,t 0 +T 1 ) A set of voltage predictors at each time;
2-2) vs. DeltaU j,x And (3) judging:
2-2-1) if DeltaU j,x > 0, let DeltaU t.inc,j ＝ΔU t.inc,j +ΔU j,x ；
And judging x: if x<t 0 +T 1 Let x=x+1, return to step 2-1); if x is greater than or equal to t 0 +T 1 Then the calculation is finished, the current delta U t.inc,j I.e. starting from time t and generating line jA final value of the maximum amplitude at which the voltage continuously increases during the prediction period;
2-2-2) if DeltaU j，x If the current delta U is less than or equal to 0, ending the calculation t.inc,j Namely, the final value of the maximum amplitude of the voltage of the bus j which continuously increases in the prediction period from the time t;
3) From the time x=t+1, Δu is calculated t.dec,j The specific steps are as follows:
3-1) from U j ＝{U j,x ,x＝t 0 ,...,t 0 +T 1 Respectively obtaining voltage predicted values U of buses j at x and x-1 moments j,x 、U j,x-1 Calculating the voltage variation delta U between the two predicted values j,x ＝U j,x -U j,x-1 The method comprises the steps of carrying out a first treatment on the surface of the Wherein U is j Representative bus j is in the prediction period (t 0 ,t 0 +T 1 ) A set of voltage predictors at each time;
3-2) vs. DeltaU j,x And (3) judging:
3-2-1) if DeltaU j,x < 0, let DeltaU t.dec,j ＝ΔU t.dec,j +ΔU j,x ；
And judging x: if x<t 0 +T 1 Let x=x+1, return to step 3-1); if x is greater than or equal to t 0 +T 1 Then the calculation is finished, the current delta U t.dec,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude value at which the predicted period of (a) decreases continuously;
3-2-2) if DeltaU j，x Not less than 0, the calculation is finished, and the current delta U is equal to or greater than 0 t.dec,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude value at which the predicted period of (a) decreases continuously;
4) Repeating the steps 2) -3) to obtain (t) 0 ,t 0 +T 1 ) Delta U corresponding to each time t in time period t.inc,j And DeltaU t.dec,j Each busbar j is calculated to be (t 0 ,t 0 +T 1 ) The upper voltage safety domain limit and the lower voltage safety domain limit in the period are as follows:
in the method, in the process of the invention,the upper and lower voltage plan limits for busbar j are respectively provided.
An embodiment of a second aspect of the present disclosure provides an automatic voltage control apparatus for reducing voltage fluctuation in a new energy collection area, including:
the base state power flow calculation module is used for temporarily obtaining a base state power flow calculation result of a new energy collection area in the power grid in each automatic voltage control period;
The sensitivity calculation module is used for calculating the voltage sensitivity of the injection active power of the high-voltage side bus of each new energy power station in the new energy collection area to each bus by using the ground state power flow calculation result;
the voltage prediction module is used for calculating the predicted value of each bus voltage in the new energy collection area at each moment in the prediction period by using the sensitivity;
and the voltage safety area calculation module is used for calculating the upper limit and the lower limit of the voltage safety area of each bus in the prediction period according to the voltage prediction value of each bus so as to be used for automatic voltage control of the prediction period.
An embodiment of a third aspect of the present disclosure proposes an electronic device, including:
at least one processor; and a memory communicatively coupled to the at least one processor;
wherein the memory stores instructions executable by the at least one processor, the instructions configured to perform an automatic voltage control method for reducing new energy collection zone voltage fluctuations as described above.
An embodiment of a fourth aspect of the present disclosure proposes a computer-readable storage medium storing computer instructions for causing the computer to execute the above-described automatic voltage control method of reducing voltage fluctuation in a new energy collection region.
The characteristics and beneficial effects of this disclosure lie in:
the method and the system are suitable for collecting regional power grids of the new energy power station in a centralized grid connection mode, the power grid voltage in the regions can generate larger fluctuation along with the intermittent characteristic of new energy power generation, in order to solve the problem, active power generation power prediction data of the new energy power station are obtained when each automatic voltage control period arrives, further voltage prediction values of all buses in the new energy collection region in a future period are calculated, then the upper limit of the voltage safety domain and the lower limit of the voltage safety domain of all buses can be further calculated, and the upper limit and the lower limit of the voltage safety domain are input into an automatic voltage control system to realize preventive control, so that reactive equipment of the power grid can be subjected to preventive regulation before the new energy power generation is rapidly changed, the fluctuation range of the voltage is reduced, the risk of power grid voltage out-of-limit is reduced, and the safety stability of the power grid voltage in the new energy collection region is guaranteed.
Drawings
Fig. 1 is an overall flowchart of an automatic voltage control method for reducing voltage fluctuation in a new energy collection area in an embodiment of the disclosure.
Detailed Description
An embodiment of a first aspect of the present disclosure provides an automatic voltage control method for reducing voltage fluctuation in a new energy collection area, where the overall flow is shown in fig. 1, and the method includes the following steps:
(1) Setting an automatic voltage control period T (T is usually 5 minutes) of the power grid;
(2) And acquiring a current ground state power flow calculation result of the power grid from an Energy Management System (EMS) system of the power grid temporarily in each control period T.
Recording that the new energy collection area contains I new energy power stations, and collecting the new energy power stations U wst ＝{W 1 ,W 2 ,…,W I W, where W i Representing an ith new energy power station, i=1, 2, …, I; note that the area contains J busbars, then busbar set U wbs ＝{B 1 ,B 2 ,…,B J A booster station bus, a switching station bus and a collecting station bus of each new energy power station, wherein B j Represents the J-th busbar, j=1, 2, …, J. Active power of each new energy power station in the ground state tide is obtained to form an active power set U of each new energy power station wstp ＝[P 1 ，P 2 ，…，P I }, wherein P i Represents the ith new energy power station W i Is a power source of the power source. Obtaining each bus voltage value in the ground state tide to form a bus voltage set U wu ＝{U 1 ,U 2 ,…,U J U, where j Representative bus B j Is a voltage value of (a).
(3) Calculating and obtaining voltage sensitivity S of high-voltage side buses of each new energy power station in the new energy collecting area by injecting active power into each bus in the collecting area based on the ground state power flow calculation result obtained in the step (2) ij ，S ij And the voltage variation of the j-th bus in the corresponding new energy collecting area is shown by injecting unit active power on the high-voltage side bus of the i-th new energy power station. Based on all new energy power stations and collecting substations in the new energy collecting area, all S is obtained ij And form an I x J order sensitivity matrix as follows:
wherein is S cv Sensitivity matrix, S ij And the voltage sensitivity of the j-th bus is actively injected on the high-voltage side bus of the i-th new energy power station, wherein the unit is kV/WM.
(4) Obtaining the current day active power generation prediction initial value D of any new energy power station i from EMS system i,n N=1, where N, N is a sequence number of an active power generation prediction initial value of the new energy power station, where v is a time interval between adjacent prediction initial values, generally v=15 minutes, and n=96; i represents the I new energy power station, i=1, 2, …, I. The prediction initial value can be from short-term new energy power generation prediction before the day, or from ultra-short-term new energy power generation prediction in the day.
(5) Calculating the predicted value of each busbar voltage in the new energy collecting area at each moment in the predicted period;
let the prediction period be t 0 To t 0 +T 1 Wherein t is 0 Representing the starting time of the predicted period, T 1 Is the length of the prediction period; in this embodiment, t is set 0 T is set to indicate the number of minutes from the current day 0 to the current time 1 The number of minutes for the rolling predicted future period (T1 typically takes 30 minutes); calculating the current time t 0 To t 0 +T 1 Predicted value U of each busbar voltage in new energy collection area in time interval per minute i,t The method comprises the following specific steps:
(5-1) setting the current prediction time to x, x being at (t) 0 ,t 0 +T 1 ) Within the range, let the initial value of x be t 0 ；
(5-2) determining an active power generation predicted value P of the new energy power station i at the time x for each new energy power station i in the new energy collection region i,x The method comprises the following steps:
if x coincides with the corresponding time of any current day active power generation prediction initial value of the new energy power station i, taking the current day active power generation prediction initial value of the new energy power station i corresponding to the x time as the active power generation prediction value of the new energy power station i at the x time; otherwise, linear interpolation of the current day active power generation prediction initial values of two adjacent new energy power stations i before and after the x moment is adopted as the active power generation prediction value of the new energy power station i at the x moment, and the expression is as follows:
where x/v denotes an integer result obtained by dividing x by v, and x% v denotes a remainder result obtained by dividing x by v.
(5-3) obtaining the active power generation predicted values of all the new energy power stations at the moment x by using the result of the step (5-2) to form an active power generation predicted value vector F of all the new energy power stations at the moment x x ＝[P 1,x P 2,x …P I,x ]；
According to set U wstp Form the corresponding vector f= [ P ] 1 P 2 …P I ]Calculating a change value vector delta F of active power of each new energy power station at the moment x x ＝F x -F＝[P 1,x -P 1 P 2,x -P 2 …P I,x -P I ]。
(5-4) according to the vector ΔF x And sensitivity matrix S cv Calculating the voltage change value vector DeltaU of each bus at the time x x ＝ΔF x *S cv ＝[ΔU 1,x ΔU 2,x …ΔU J,x ]According to set U wu ＝{U 1 ,U 2 ,…,U J Voltage vector U formed 0 ＝[U 1 U 2 …U J ]The vector formed by the predicted values of the voltages of all the buses at the moment x is obtained to be U x ＝U 0 +ΔU x Wherein U is x ＝[U 1 +ΔU 1,x U 2 +ΔU 2,x …U J +ΔU J,x ]＝[U 1,x U 2,x …U J,x ]
(5-5) making x=x+1, then returning to the step (5-2) again, and calculating the predicted value of each busbar voltage in the new energy collecting area at the next predicted time until the current predicted time is t 0 +T 1 Obtaining the predicted values of the voltages of all buses in the new energy collecting area at all times of the predicted period, and constructing the predicted values at (t) 0 ,t 0 +T 1 ) The set of each bus voltage predicted value in the new energy collection area of the time period:
(6) Calculating the upper limit and the lower limit of a voltage safety area of each bus in the new energy collection area in a prediction period;
in the embodiment, step (5) acquires the voltage predicted value U of each busbar j in the new energy collection region at each moment in the predicted period j ＝{U j,x ,x＝t 0 ,...,t 0 +T 1 The voltage of each busbar j can be further calculatedUpper limit of safety areaAnd a lower limit of the voltage safety region->The method comprises the following specific steps:
(6-1) setting t to (t) 0 ,t 0 +T 1 ) At any time within the range, an initial value DeltaU is set t.inc,j ＝0，ΔU t.dec,j =0; wherein DeltaU t.inc,j Represents the maximum amplitude of the continuous increase in voltage of bus j from time t over the predicted period, deltaU t.dec,j Representing the maximum amplitude at which the voltage of bus j decreases continuously during the predicted period from time t.
(6-2) calculating Δu from the time x=t+1 t.inc,j The specific steps are as follows:
(6-2-1) from U j Respectively obtaining voltage predicted values U of buses j at x and x-1 moments j,x 、U j,x-1 Calculating the voltage variation delta U between the two predicted values j,x ＝U j,x -U j,x-1 ；
(6-2-2) vs. DeltaU j,x And (3) judging:
(6-2-2-1) if DeltaU j,x > 0, let DeltaU t.inc,j ＝ΔU t.inc,j +ΔU j,x ；
Then, judging x: if x<t 0 +T 1 Let x=x+1, return to step (6-2-1); if x is greater than or equal to t 0 +T 1 Then the calculation is finished, the current delta U t.inc,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude that occurs continuously increasing for the predicted period of time;
(6-2-2-2) if DeltaU j，x If the current delta U is less than or equal to 0, ending the calculation t.inc,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude that occurs continuously increasing for the predicted period of time; (6-3) calculating Δu from the time x=t+1 t.dec,j Final of (2)The specific steps are as follows:
(6-3-1) from U j ＝{U j,x ,x＝t 0 ,...,t 0 +T 1 Respectively obtaining voltage predicted values U of buses j at x and x-1 moments j,x 、U j,x-1 Calculating the voltage variation delta U between the two predicted values j,x ＝U j,x -U j,x-1 The method comprises the steps of carrying out a first treatment on the surface of the Wherein U is j Representative bus j is in the prediction period (t 0 ,t 0 +T 1 ) A set of voltage predictors at each time;
(6-3-2) pair DeltaU j,x And (3) judging:
(6-3-2-1) if DeltaU j,x < 0, let DeltaU t.dec,j ＝ΔU t.dec,j +ΔU j,x ；
Then, judging x: if x<t 0 +T 1 Let x=x+1, repeat step (6-3-1) to update Δu j,x The method comprises the steps of carrying out a first treatment on the surface of the If x is greater than or equal to t 0 +T 1 Then the calculation is finished, the current delta U t.dec,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude value at which the predicted period of (a) decreases continuously;
(6-3-2-2) if DeltaU j，x Not less than 0, the calculation is finished, and the current delta U is equal to or greater than 0 t.dec,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude value at which the predicted period of (a) decreases continuously;
(6-4) repeating the steps (6-2) to (6-3) to obtain (t) 0 ,t 0 +T 1 ) Delta U corresponding to each time t in time period t.inc,j And DeltaU t.dec,j Each bus j considering the voltage trend is calculated at (t 0 ,t 0 +T 1 ) The upper voltage safety domain limit and the lower voltage safety domain limit in the period are as follows:
wherein t=t 0 ,…,t 0 +T 1 For each predicted minute time, the time of day,the upper limit value and the lower limit value of the voltage plan of the busbar j which are manually compiled in advance are respectively.
The upper limit and the lower limit of the voltage safety domain of each bus obtained by the calculation are only applicable to the prediction period (t 0 ,t 0 +T 1 ) When the next prediction period comes, the upper voltage safety domain limit and the lower voltage safety domain limit of each bus bar need to be recalculated.
(7) And (3) inputting the upper limit and the lower limit of the voltage safety domain of each bus in the new energy collection area in the prediction period obtained by the calculation in the step (6) into an automatic voltage control module of a dispatching monitoring system of a power grid dispatching center, and issuing corresponding control instructions to a near-area power plant and a transformer substation by the automatic voltage control module to execute, so that the voltage fluctuation of the new energy collection area in the prediction period is reduced, and the safe and stable operation of the voltage is ensured.
To achieve the above embodiments, an embodiment of a second aspect of the present disclosure provides an automatic voltage control device for reducing voltage fluctuation in a new energy collection area, including:
the base state power flow calculation module is used for temporarily obtaining a base state power flow calculation result of a new energy collection area in the power grid in each automatic voltage control period;
the sensitivity calculation module is used for calculating the voltage sensitivity of the injection active power of the high-voltage side bus of each new energy power station in the new energy collection area to each bus by using the ground state power flow calculation result;
the voltage prediction module is used for calculating the predicted value of each bus voltage in the new energy collection area at each moment in the prediction period by using the sensitivity;
and the voltage safety area calculation module is used for calculating the upper limit and the lower limit of the voltage safety area of each bus in the prediction period according to the voltage prediction value of each bus so as to be used for automatic voltage control of the prediction period.
To achieve the above embodiments, an embodiment of a third aspect of the present disclosure proposes an electronic device, including:
at least one processor; and a memory communicatively coupled to the at least one processor;
Wherein the memory stores instructions executable by the at least one processor, the instructions configured to perform an automatic voltage control method for reducing new energy collection zone voltage fluctuations as described above.
To achieve the above embodiments, an embodiment of a fourth aspect of the present disclosure proposes a computer-readable storage medium storing computer instructions for causing the computer to execute the above-described automatic voltage control method for reducing voltage fluctuation in a new energy collection area.
It should be noted that the computer readable medium described in the present disclosure may be a computer readable signal medium or a computer readable storage medium, or any combination of the two. The computer readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any of the foregoing. More specific examples of the computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this disclosure, a computer-readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. In the present disclosure, however, the computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, with the computer-readable program code embodied therein. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: electrical wires, fiber optic cables, RF (radio frequency), and the like, or any suitable combination of the foregoing.
The computer readable medium may be contained in the electronic device; or may exist alone without being incorporated into the electronic device. The computer-readable medium carries one or more programs which, when executed by the electronic device, cause the electronic device to perform an automatic voltage control method of reducing voltage fluctuations in a new energy collection region of the above-described embodiment.
Computer program code for carrying out operations of the present disclosure may be written in one or more programming languages, including an object oriented programming language such as Java, smalltalk, C ++ and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computer (for example, through the Internet using an Internet service provider).
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present application. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, the different embodiments or examples described in this specification and the features of the different embodiments or examples may be combined and combined by those skilled in the art without contradiction.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present application, the meaning of "plurality" is at least two, such as two, three, etc., unless explicitly defined otherwise.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and further implementations are included within the scope of the preferred embodiment of the present application in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present application.
Logic and/or steps represented in the flowcharts or otherwise described herein, e.g., a ordered listing of executable instructions for implementing logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium may even be paper or other suitable medium upon which the program is printed, as the program may be electronically captured, via, for instance, optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner, if necessary, and then stored in a computer memory.
It is to be understood that portions of the present application may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
Those of ordinary skill in the art will appreciate that all or part of the steps carried out in the method of the above-described embodiments may be implemented by a program to instruct related hardware, and the program may be stored in a computer readable storage medium, where the program when executed includes one or a combination of the steps of the method embodiments.
In addition, each functional unit in each embodiment of the present application may be integrated in one processing module, or each unit may exist alone physically, or two or more units may be integrated in one module. The integrated modules may be implemented in hardware or in software functional modules. The integrated modules may also be stored in a computer readable storage medium if implemented as software functional modules and sold or used as a stand-alone product.
The above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, or the like. Although embodiments of the present application have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the application, and that variations, modifications, alternatives, and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the application.
Claims (4)
1. An automatic voltage control method for reducing voltage fluctuation of a new energy collection area is characterized by comprising the following steps:
temporarily acquiring a base state power flow calculation result of a new energy collection area in the power grid in each automatic voltage control period;
calculating the voltage sensitivity of the injection active power of the high-voltage side bus of each new energy power station in the new energy collection area to each bus by using the ground state power flow calculation result;
calculating the predicted value of each busbar voltage in the new energy collecting area at each moment in the predicted period by using the sensitivity;
calculating the upper limit and the lower limit of a voltage safety area of each bus in the prediction period according to the predicted value of each bus voltage, so as to be used for automatic voltage control of the prediction period;
The base state power flow calculation result of the new energy collection area in the power grid comprises the following steps:
active power of each new energy power station and each bus voltage value in the new energy collection area;
wherein the active power of each new energy power station forms an active power set U of each new energy power station wstp ＝{P 1 ,P 2 ,…,P I Each busbar voltage value forms a busbar voltage set U wu ＝{U 1 ,U 2 ,…,U J -a }; i is new energy contained in the new energy collection areaThe total number of the power stations, J is the total number of the buses contained in the new energy collecting area;
calculating the predicted value of each bus voltage in the new energy collection area at each moment in the predicted period by using the sensitivity, wherein the method comprises the following steps:
1) Let the prediction period be t 0 To t 0 +T 1 Wherein t is 0 Representing the starting time of the predicted period, T 1 Is the length of the prediction period;
2) Let the current prediction time be x, and the initial value of x be t 0 ；
3) For each new energy power station i in the new energy collection area, determining an active power generation predicted value P of the new energy power station i at the moment x i,x The method is characterized by comprising the following steps:
if x coincides with the corresponding time of any current day active power generation prediction initial value of the new energy power station i, taking the current day active power generation prediction initial value of the new energy power station i corresponding to the x time as the active power generation prediction value of the new energy power station i at the x time; otherwise, linear interpolation of current day active power generation prediction initial values of two adjacent new energy power stations before and after the x moment is adopted as the active power generation prediction value of the new energy power station i at the x moment, and the expression is as follows:
Wherein, the current day active power generation prediction initial value D of the new energy power station i i,n Acquiring n=1 from a power grid energy management system, wherein N, N is a serial number of an active power generation prediction initial value of a new energy power station, and v is a time interval between adjacent prediction initial values; x/v represents the integer result of dividing x by v, and x% v represents the remainder result of dividing x by v;
4) Utilizing the result of the step 3) to establish an active power generation predictive value vector F of all new energy power stations at the moment x x ＝[P 1,x ，P 2,x ，…，P I,x ]；
According to set U wstp Form the corresponding vector f= [ P ] 1 ，P 2 ，…，P I ]Calculating a change value vector delta F of active power of each new energy power station at the moment x x ＝F x -F＝[P 1,x -P 1 ，P 2,x -P 2 ，…，P I,x -P I ]；
5) Calculating voltage change value vector delta U of each bus at x time x ＝ΔF x *S cv ＝[ΔU 1,x ，ΔU 2,x ，…，ΔU J,x ]According to set U wu ＝{U 1 ,U 2 ,…,U J Voltage vector U formed 0 ＝[U 1 ，U 2 ，…，U J ]The vector formed by the predicted values of the voltages of all the buses at the moment x is obtained to be U x ＝U 0 +ΔU x Wherein U is x ＝[U 1 +ΔU 1,x ，U 2 +ΔU 2,x ，…，U J +ΔU J,x ]＝[U 1,x ，U 2,x ，…，U J,x ]；
Wherein S is cv For the sensitivity matrix:
wherein S is ij The voltage sensitivity of the j-th bus is injected on the high-voltage side bus of the i-th new energy power station;
6) Let x=x+1, then return to step 3); until the current prediction time is t 0 +T 1 Obtaining the predicted values of the voltages of all buses in the new energy collecting area at all times of the predicted period, and constructing the predicted values at (t) 0 ,t 0 +T 1 ) The set of each bus voltage predicted value in the new energy collection area of the time period:
According to the predicted value of each bus voltage, calculating the upper limit of the voltage safety area and the lower limit of the voltage safety area of each bus in the predicted period, for automatic voltage control in the predicted period, including:
1) Let t be the prediction period (t 0 ,t 0 +T 1 ) Any time in the process;
for each busbar j, an initial value DeltaU is set t.inc,j ＝0，ΔU t.dec,j =0; wherein DeltaU t.inc,j Represents the maximum amplitude of the continuous increase in voltage of bus j from time t over the predicted period, deltaU t.dec,j Representing the maximum amplitude of continuous decrease of the voltage of the busbar j in the prediction period from the moment t;
2) From the time x=t+1, Δu is calculated t.inc,j The specific steps are as follows:
2-1) from U j ＝{U j,x ,x＝t 0 ,...,t 0 +T 1 Respectively obtaining voltage predicted values U of buses j at x and x-1 moments j,x 、U j,x-1 Calculating the voltage variation delta U between two predicted values j,x ＝U j,x -U j,x-1 The method comprises the steps of carrying out a first treatment on the surface of the Wherein U is j Representative bus j is in the prediction period (t 0 ,t 0 +T 1 ) A set of voltage predictors at each time;
2-2) vs. DeltaU j,x And (3) judging:
2-2-1) if DeltaU j,x > 0, let DeltaU t.inc,j ＝ΔU t.inc,j +ΔU j,x ；
And judging x: if x<t 0 +T 1 Let x=x+1, return to step 2-1); if x is greater than or equal to t 0 +T 1 Then the calculation is finished, the current delta U t.inc,j Namely, the final value of the maximum amplitude of the voltage of the bus j which continuously increases in the prediction period from the time t;
2-2-2) if DeltaU j,x If the current delta U is less than or equal to 0, ending the calculation t.inc,j Namely, the final value of the maximum amplitude of the voltage of the bus j which continuously increases in the prediction period from the time t;
3) From the time x=t+1, Δu is calculated t.dec,j The specific steps are as follows:
3-1) from U j ＝{U j,x ,x＝t 0 ,...,t 0 +T 1 Respectively obtaining voltage predicted values U of buses j at x and x-1 moments j,x 、U j,x-1 Calculating the voltage variation delta U between the two predicted values j,x ＝U j,x -U j,x-1 The method comprises the steps of carrying out a first treatment on the surface of the Wherein U is j Representative bus j is in the prediction period (t 0 ,t 0 +T 1 ) A set of voltage predictors at each time;
3-2) vs. DeltaU j,x And (3) judging:
3-2-1) if DeltaU j,x < 0, let DeltaU t.dec,j ＝ΔU t.dec,j +ΔU j,x ；
And judging x: if x<t 0 +T 1 Let x=x+1, return to step 3-1); if x is greater than or equal to t 0 +T 1 Then the calculation is finished, the current delta U t.dec,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude value at which the predicted period of (a) decreases continuously;
3-2-2) if DeltaU j,x Not less than 0, the calculation is finished, and the current delta U is equal to or greater than 0 t.dec,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude value at which the predicted period of (a) decreases continuously;
4) Repeating the steps 2) -3) to obtain (t) 0 ,t 0 +T 1 ) Delta U corresponding to each time t in time period t.inc,j And DeltaU t.dec,j Each busbar j is calculated to be (t 0 ,t 0 +T 1 ) The upper voltage safety domain limit and the lower voltage safety domain limit in the period are as follows:
In the method, in the process of the invention,the upper and lower voltage plan limits for busbar j are respectively provided.
2. An automatic voltage control device for reducing voltage fluctuation in a new energy collection area, comprising:
the base state power flow calculation module is used for temporarily obtaining a base state power flow calculation result of a new energy collection area in the power grid in each automatic voltage control period;
the sensitivity calculation module is used for calculating the voltage sensitivity of the injection active power of the high-voltage side bus of each new energy power station in the new energy collection area to each bus by using the ground state power flow calculation result;
the voltage prediction module is used for calculating the predicted value of each bus voltage in the new energy collection area at each moment in the prediction period by using the sensitivity;
the voltage safety area calculation module is used for calculating the upper limit and the lower limit of the voltage safety area of each bus in the prediction period according to the voltage prediction value of each bus so as to be used for automatic voltage control of the prediction period;
the base state power flow calculation result of the new energy collection area in the power grid comprises the following steps:
active power of each new energy power station and each bus voltage value in the new energy collection area;
Wherein the active power of each new energy power station forms an active power set U of each new energy power station wstp ＝{P 1 ,P 2 ,…,P I Each busbar voltage value forms a busbar voltage set U wu ＝{U 1 ,U 2 ,…,U J -a }; i is the total number of new energy power stations contained in the new energy collecting area, and J is the total number of buses contained in the new energy collecting area;
calculating the predicted value of each bus voltage in the new energy collection area at each moment in the predicted period by using the sensitivity, wherein the method comprises the following steps:
1) Let the prediction period be t 0 To t 0 +T 1 Wherein t is 0 Representing the starting time of the predicted period, T 1 Is the length of the prediction period;
2) Let the current prediction time be x, the initial of xWith value t 0 ；
3) For each new energy power station i in the new energy collection area, determining an active power generation predicted value P of the new energy power station i at the moment x i,x The method is characterized by comprising the following steps:
if x coincides with the corresponding time of any current day active power generation prediction initial value of the new energy power station i, taking the current day active power generation prediction initial value of the new energy power station i corresponding to the x time as the active power generation prediction value of the new energy power station i at the x time; otherwise, linear interpolation of current day active power generation prediction initial values of two adjacent new energy power stations before and after the x moment is adopted as the active power generation prediction value of the new energy power station i at the x moment, and the expression is as follows:
Wherein, the current day active power generation prediction initial value D of the new energy power station i i,n Acquiring n=1 from a power grid energy management system, wherein N, N is a serial number of an active power generation prediction initial value of a new energy power station, and v is a time interval between adjacent prediction initial values; x/v represents the integer result of dividing x by v, and x% v represents the remainder result of dividing x by v;
4) Utilizing the result of the step 3) to establish an active power generation predictive value vector F of all new energy power stations at the moment x x ＝[P 1,x ，P 2,x ，…，P I,x ]；
According to set U wstp Form the corresponding vector f= [ P ] 1 ，P 2 ，…，P I ]Calculating a change value vector delta F of active power of each new energy power station at the moment x x ＝F x -F＝[P 1,x -P 1 ，P 2,x -P 2 ，…，P I,x -P I ]；
5) Calculating voltage change value vector delta U of each bus at x time x ＝ΔF x *S cv ＝[ΔU 1,x ，ΔU 2,x ，…，ΔU J,x ]According to set U wu ＝{U 1 ,U 2 ,…,U J Formed by }Voltage vector U 0 ＝[U 1 ，U 2 ，…，U J ]The vector formed by the predicted values of the voltages of all the buses at the moment x is obtained to be U x ＝U 0 +ΔU x Wherein U is x ＝[U 1 +ΔU 1,x ，U 2 +ΔU 2,x ，…，U J +ΔU J,x ]＝[U 1,x ，U 2,x ，…，U J,x ]；
Wherein S is cv For the sensitivity matrix:
wherein S is ij The voltage sensitivity of the j-th bus is injected on the high-voltage side bus of the i-th new energy power station;
6) Let x=x+1, then return to step 3); until the current prediction time is t 0 +T 1 Obtaining the predicted values of the voltages of all buses in the new energy collecting area at all times of the predicted period, and constructing the predicted values at (t) 0 ,t 0 +T 1 ) The set of each bus voltage predicted value in the new energy collection area of the time period:
According to the predicted value of each bus voltage, calculating the upper limit of the voltage safety area and the lower limit of the voltage safety area of each bus in the predicted period, for automatic voltage control in the predicted period, including:
1) Let t be the prediction period (t 0 ,t 0 +T 1 ) Any time in the process;
for each busbar j, an initial value DeltaU is set t.inc,j ＝0，ΔU t.dec,j =0; wherein DeltaU t.inc,j Represents the maximum amplitude of the continuous increase in voltage of bus j from time t over the predicted period, deltaU t.dec,j The voltage representing the bus j starting from time t occurs in the predicted periodA continuously decreasing maximum amplitude;
2) From the time x=t+1, Δu is calculated t.inc,j The specific steps are as follows:
2-1) from U j ＝{U j,x ,x＝t 0 ,...,t 0 +T 1 Respectively obtaining voltage predicted values U of buses j at x and x-1 moments j,x 、U j,x-1 Calculating the voltage variation delta U between two predicted values j,x ＝U j,x -U j,x-1 The method comprises the steps of carrying out a first treatment on the surface of the Wherein U is j Representative bus j is in the prediction period (t 0 ,t 0 +T 1 ) A set of voltage predictors at each time;
2-2) vs. DeltaU j,x And (3) judging:
2-2-1) if DeltaU j,x > 0, let DeltaU t.inc,j ＝ΔU t.inc,j +ΔU j,x ；
And judging x: if x<t 0 +T 1 Let x=x+1, return to step 2-1); if x is greater than or equal to t 0 +T 1 Then the calculation is finished, the current delta U t.inc,j Namely, the final value of the maximum amplitude of the voltage of the bus j which continuously increases in the prediction period from the time t;
2-2-2) if DeltaU j,x If the current delta U is less than or equal to 0, ending the calculation t.inc,j Namely, the final value of the maximum amplitude of the voltage of the bus j which continuously increases in the prediction period from the time t;
3) From the time x=t+1, Δu is calculated t.dec,j The specific steps are as follows:
3-1) from U j ＝{U j,x ,x＝t 0 ,...,t 0 +T 1 Respectively obtaining voltage predicted values U of buses j at x and x-1 moments j,x 、U j,x-1 Calculating the voltage variation delta U between the two predicted values j,x ＝U j,x -U j,x-1 The method comprises the steps of carrying out a first treatment on the surface of the Wherein U is j Representative bus j is in the prediction period (t 0 ,t 0 +T 1 ) A set of voltage predictors at each time;
3-2) vs. DeltaU j,x And (3) judging:
3-2-1) if DeltaU j,x < 0, let DeltaU t.dec,j ＝ΔU t.dec,j +ΔU j,x ；
And judging x: if x<t 0 +T 1 Let x=x+1, return to step 3-1); if x is greater than or equal to t 0 +T 1 Then the calculation is finished, the current delta U t.dec,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude value at which the predicted period of (a) decreases continuously;
3-2-2) if DeltaU j,x Not less than 0, the calculation is finished, and the current delta U is equal to or greater than 0 t.dec,j Namely, the voltage of the bus j from the time t is within (t 0 ,t 0 +T 1 ) A final value of the maximum amplitude value at which the predicted period of (a) decreases continuously;
4) Repeating the steps 2) -3) to obtain (t) 0 ,t 0 +T 1 ) Delta U corresponding to each time t in time period t.inc,j And DeltaU t.dec,j Each busbar j is calculated to be (t 0 ,t 0 +T 1 ) The upper voltage safety domain limit and the lower voltage safety domain limit in the period are as follows:
In the method, in the process of the invention,the upper and lower voltage plan limits for busbar j are respectively provided.
3. An electronic device, comprising:
at least one processor; and a memory communicatively coupled to the at least one processor;
wherein the memory stores instructions executable by the at least one processor, the instructions configured to perform the method of claim 1.
4. A computer readable storage medium storing computer instructions for causing the computer to perform the method of claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110925892.6A CN113629766B (en) | 2021-08-12 | 2021-08-12 | Automatic voltage control method and device for reducing voltage fluctuation of new energy collection area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110925892.6A CN113629766B (en) | 2021-08-12 | 2021-08-12 | Automatic voltage control method and device for reducing voltage fluctuation of new energy collection area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113629766A CN113629766A (en) | 2021-11-09 |
| CN113629766B true CN113629766B (en) | 2024-01-26 |
Family
ID=78385036
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110925892.6A Active CN113629766B (en) | 2021-08-12 | 2021-08-12 | Automatic voltage control method and device for reducing voltage fluctuation of new energy collection area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113629766B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114374211B (en) * | 2021-12-10 | 2024-01-30 | 国网重庆市电力公司市区供电分公司 | Automatic voltage control method and device considering energy storage active plan |
| CN116799878B (en) * | 2023-08-18 | 2023-11-07 | 国网山西省电力公司运城供电公司 | Dynamic reactive power coordination control method for new energy collection area based on power generation trend prediction |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103218531A (en) * | 2013-04-15 | 2013-07-24 | 河海大学 | Remote voltage control mode accounted node type extended continuous load flow algorithm |
| CN103545814A (en) * | 2013-09-24 | 2014-01-29 | 国家电网公司 | Coordinated voltage control method for alternating-current and direct-current hybrid power system |
| CN105977988A (en) * | 2016-06-08 | 2016-09-28 | 清华大学 | Active and reactive power coordination control method for improving new energy influx area voltage security |
| CN112865108A (en) * | 2021-01-11 | 2021-05-28 | 国网山西省电力公司忻州供电公司 | Power grid automatic voltage control simulation method based on continuous power flow simulation |
-
2021
- 2021-08-12 CN CN202110925892.6A patent/CN113629766B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103218531A (en) * | 2013-04-15 | 2013-07-24 | 河海大学 | Remote voltage control mode accounted node type extended continuous load flow algorithm |
| CN103545814A (en) * | 2013-09-24 | 2014-01-29 | 国家电网公司 | Coordinated voltage control method for alternating-current and direct-current hybrid power system |
| CN105977988A (en) * | 2016-06-08 | 2016-09-28 | 清华大学 | Active and reactive power coordination control method for improving new energy influx area voltage security |
| CN112865108A (en) * | 2021-01-11 | 2021-05-28 | 国网山西省电力公司忻州供电公司 | Power grid automatic voltage control simulation method based on continuous power flow simulation |
Non-Patent Citations (4)
| Title |
|---|
| A novel method of power flow analysis with UPFC considering limit violations of variable;Qiang Tong .et al;2011 International Conference on Electric Information and Control Engineering;全文 * |
| Automatic voltage control of wind farms based on model predictive control theory;Xu Fengda.et al;Automation of Electric Power Systems;第39卷(第7期);全文 * |
| 基于模型预测控制理论的风电场自动电压控制;徐峰达 等;电力系统自动化;第39卷(第7期);全文 * |
| 肖浩 等.分布式电源对配电网电压的影响分析及其 优化控制策略.2016,第31卷(第S1期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113629766A (en) | 2021-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Weckx et al. | Optimal local reactive power control by PV inverters | |
| Bidgoli et al. | Combined local and centralized voltage control in active distribution networks | |
| Zeraati et al. | Distributed control of battery energy storage systems for voltage regulation in distribution networks with high PV penetration | |
| CN113629766B (en) | Automatic voltage control method and device for reducing voltage fluctuation of new energy collection area | |
| Glinkowski et al. | Advances in wind energy technologies in the context of smart grid | |
| Sakamuri et al. | Coordinated voltage control in offshore HVDC connected cluster of wind power plants | |
| US20170317498A1 (en) | Resiliency Controller for Voltage Regulation in Microgrids | |
| Xie et al. | Robust model predictive control based voltage regulation method for a distribution system with renewable energy sources and energy storage systems | |
| Rashid et al. | Combining Volt/Var & Volt/Watt modes to increase PV hosting capacity in LV distribution networks | |
| Parthasarathy et al. | Integration and control of lithium-ion BESSs for active network management in smart grids: Sundom smart grid backup feeding case | |
| Radosavljević | Voltage regulation in LV distribution networks with PV generation and battery storage | |
| Heidari Yazdi et al. | Over-voltage regulation of distribution networks by coordinated operation of PV inverters and demand side management program | |
| CN117374974A (en) | Distribution network scheduling method, system, medium and equipment | |
| Andreev et al. | Novel method for setting up the relay protection of power systems containing renewable energy sources and hydrogen energy storage systems | |
| CN114374211B (en) | Automatic voltage control method and device considering energy storage active plan | |
| Coppo et al. | Sliding time windows assessment of storage systems capability for providing ancillary services to transmission and distribution grids | |
| Abaci et al. | Voltage stability improvement with coordinated ULTC–STATCOM controller and VSC-HVDC in high wind penetration cases | |
| CN113872238B (en) | Automatic voltage control method and device for power system, electronic equipment and storage medium | |
| Bosisio et al. | A procedure for the optimal management of medium-voltage AC networks with distributed generation and storage devices | |
| Bagheri | Energy storage system coordinated with phase-shifting transformer and dynamic rating equipment for optimal operation of wind-rich smart power networks | |
| Kim et al. | Economic analysis on multi-terminal VSC HVDC systems with wind farms based on hierarchical optimal power flow with stability constraint | |
| Zhang et al. | Optimal Planning of Hybrid AC/DC Low-Voltage Distribution Networks Considering DC Conversion of Three-Phase Four-Wire Low-Voltage AC Systems | |
| Cheah | Offshore wind integration through high voltage direct current systems | |
| Yang et al. | An adaptive control for supporting village power grid integrating residential PV power generation | |
| Marten et al. | Kriegers Flak Combined Grid Solution-novel double use of offshore equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |