CN114336598B - Regional power grid automatic voltage control method based on parallel simulation - Google Patents
Regional power grid automatic voltage control method based on parallel simulation Download PDFInfo
- Publication number
- CN114336598B CN114336598B CN202111517126.2A CN202111517126A CN114336598B CN 114336598 B CN114336598 B CN 114336598B CN 202111517126 A CN202111517126 A CN 202111517126A CN 114336598 B CN114336598 B CN 114336598B
- Authority
- CN
- China
- Prior art keywords
- bus
- power grid
- simulation
- voltage
- monitored
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 196
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000004364 calculation method Methods 0.000 claims abstract description 59
- 238000011217 control strategy Methods 0.000 claims description 26
- 239000003990 capacitor Substances 0.000 claims description 22
- 101150104728 GPR88 gene Proteins 0.000 claims description 18
- 238000005457 optimization Methods 0.000 claims description 10
- 238000012544 monitoring process Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000012797 qualification Methods 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 241000219061 Rheum Species 0.000 description 1
- 235000009411 Rheum rhabarbarum Nutrition 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- 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
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a regional power grid automatic voltage control method based on parallel simulation, which is technically characterized by comprising the following steps of: reading an initial section from a historical section automatically controlled by regional power grid voltage; setting a simulation task set corresponding to the regional power grid simulation case, and performing simulation calculation and index statistics on the simulation task set; outputting reactive voltage control indexes of the power grid; analyzing and deciding the reactive voltage control index data of the power grid to obtain upper and lower limit value parameters of voltage operation of the bus needing to be monitored; and updating the original states of the upper limit value and the lower limit value of the voltage operation of the bus to be monitored in the automatic voltage control system for real-time operation of the regional power grid by utilizing the upper limit value and the lower limit value parameters of the voltage operation of the bus to be monitored. The invention has reasonable design, reduces the active loss of the power grid while ensuring the qualification of the 10kV bus voltage of the regional power grid, thereby improving the automatic voltage control level of the regional power grid and the stability of the power system.
Description
Technical Field
The invention belongs to the technical field of power grid dispatching, and particularly relates to an automatic voltage control method for a regional power grid based on parallel simulation.
Background
The automatic voltage control (Automatic Voltage Control, AVC for short) system is an important means for realizing safe, economical and high-quality operation of regional power grids. The AVC control parameters are set to have important effects in improving the voltage qualification level of the regional power grid and reducing the active loss of the regional power grid. In practical application, with the promotion of regional power grid construction and the rapid change of electricity consumption, the power grid has the characteristics of nonlinear load, strong uncertainty, strong coupling and the like, and the power grid dispatching can not timely adjust AVC control parameters according to regional power grid load change, so that the problem of unreasonable voltage reactive power control exists.
Because 220kV and 110kV buses of a transformer substation in the regional power grid model are used as high-voltage sides to acquire electric energy from an external power grid, and the electric energy is provided for regional load users through low-voltage side 10kV buses after voltage transformation of a transformer, the voltage qualification level of the 10kV buses of a user side of the transformer substation is more concerned in the operation of the regional power grid.
In conventional regional power grid AVC control, the limit value parameter range of a 10kV bus at a user end of a transformer substation is relatively fixed, the limit value parameter ranges of 220kV and 110kV buses at a power receiving end of the transformer substation are generally set in a certain range according to the experience of scheduling, and the support of relevant data is lacking in the aspects of considering regional power grid safety, economy and high-quality operation, so that the obtained AVC control parameters often cannot meet the operation requirement of power grid voltage stability.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides an automatic voltage control method for a regional power grid based on parallel simulation, which can reduce the active loss of the power grid and improve the automatic voltage control level of the regional power grid.
The invention solves the technical problems in the prior art by adopting the following technical scheme:
an automatic voltage control method for a regional power grid based on parallel simulation comprises the following steps:
step 1, setting an automatic voltage control period T of a regional power grid c ；
Step 2, setting the date T of the current clock of the regional power grid m Reading the date T from a historical section of regional power grid voltage automatic control m Corresponding power grid model M and initial section F at time 0:00 m,0 ；
Step 3, setting the regional power grid clock date as T m Load prediction information Ld for +1 day fre ；
Step 4, setting a power grid model in the simulation calculation at T m Controllable switch state Cb of +1 day m,i ：
Step 5, setting scene information F corresponding to regional power grid simulation cases sch ：
Step 6, setting a simulation task set T corresponding to the regional power grid simulation case sch According to the scene information F of step 5 sch Bus limit value parameter variable setting information to be monitored to generate a simulation task set T sch ；
Step 7, step-by-stepSimulation task set T of step 6 sch Performing simulation calculation and index statistics on one task of the system;
step 8, repeating the step 7 by adopting a task parallel method to finish the simulation task set T of the step 6 sch The simulation calculation of all simulation tasks and the statistics of reactive power voltage control indexes of the power grid output the reactive power voltage control index K10 of the power grid p,i ；
Step 9, controlling the reactive voltage control index K10 of the power grid in step 8 p,i The data are analyzed and decided to obtain a group of voltage operation upper and lower limit value parameters which can be used for the regional power grid real-time operation automatic voltage control system and need to monitor buses;
step 10, the clock date of the regional power grid is T m At the time of 0:00 of +1 day, the original state Bs of the upper limit value and the lower limit value of the voltage operation of the bus to be monitored in the automatic voltage control system for real-time operation of the regional power grid is obtained by utilizing the upper limit value and the lower limit value parameters of the voltage operation of the bus to be monitored ora The array is updated.
The invention has the advantages and positive effects that:
1. according to the invention, a model and data of real-time operation of the regional power grid are utilized, a reactive voltage control effect index data set of the power grid is rapidly formed by using a parallel simulation calculation method from the perspective of the whole system, and an AVC control parameter set value of an automatic voltage control system for real-time operation of the regional power grid is obtained through decision analysis of the index data set, so that the active loss of the power grid is reduced while the voltage of a 10kV bus of the regional power grid is qualified, and the automatic voltage control level of the regional power grid is improved.
2. The method is used for forming the transformer substation model in the automatic voltage control, automatically generates the related transformer substation control equipment model and specific parameters and plan curves through the power grid model at the beginning of the control according to the defined model generation rule, and can reduce the manual maintenance workload and maintenance errors of the transformer substation model in the automatic voltage control process by applying the control model in the control so as to improve the automatic voltage control level and the stability of the power system.
Drawings
Fig. 1 is a schematic diagram of an automatic voltage control method for a regional power grid according to the present invention.
Detailed Description
Embodiments of the present invention are described in further detail below with reference to the accompanying drawings.
An automatic voltage control method for a regional power grid based on parallel simulation, as shown in fig. 1, comprises the following steps:
step 1, setting an automatic voltage control period of a regional power grid as T c ；
Step 2, setting the current clock date of the regional power grid as T m Reading the date T from a historical section of regional power grid voltage automatic control m Corresponding power grid model M and initial section F at time 0:00 m,0 ，F m,0 The method comprises the steps of carrying out information such as load active power, load reactive power, switch states, bus voltage limit states needing to be monitored and the like in a regional power grid:
F m,0 ＝{Ld 0,p ,Ld 0,q ,Cb 0 ,Bs 0 ,Loss 0 }；
Wherein the subscript m corresponds to the date T m ，Loss 0 Is the active loss of the initial section.
The specific implementation method of the steps is as follows:
step 2.1, setting the serial number of the adjustable load in the regional power grid as L, wherein l= … … L, L represents the total number of the adjustable loads in the power grid model M, and Ld 0,p And Ld 0,q For an initial active and initial reactive of the adjustable load at time 0:00:
Ld 0,p ＝[l,l＝1,..L]{P 0,l }；Ld 0,q ＝[l,l＝1,..L]{Q 0,l }；
wherein P is 0,l And Q 0,l The initial active power and the initial reactive power of the first adjustable load at the moment 0:00 respectively;
step 2.2, setting the serial number of the controllable switches in the regional power grid as C, c= … … C, wherein C represents the total number of the controllable switches in the power grid model M, cb 0 An initial state array for the controllable switch:
Cb 0 ＝[c,c＝1,..C]{Open 0,c }；
wherein Open is 0,c The switching state of the c-th controllable switch at the moment 0:00 is as follows:
step 2.3, setting the serial number of the bus needing to be monitored in the regional power grid as S, wherein s= … … S, and S represents the total number of the bus needing to be monitored in the power grid model M, and Bs 0 Initial state array for bus voltage limit value to be monitored:
Bs 0 ＝[s,s＝1,..S]{K s,i ,A s,i ,U s,i,val ,U s,i,min ,U s,i,max }(0≤i≤288)；
wherein K is s,i The voltage grade information corresponding to the s-th adjustable mother is obtained; a is that s,i The information of the region corresponding to the bus needing to be monitored is the s-th bus; u (U) s,i,val The voltage value information of the bus needing to be monitored at the ith moment is the s th bus; u (U) s,i,min And U s,i,max The operation upper limit value and the operation lower limit value of the bus needing to be monitored at the ith moment are respectively the s th bus.
Step 3, setting the regional power grid clock date as T m Load prediction information of +1 day is Ld fre ：
Wherein the subscript i is the load prediction time,respectively, the active and reactive information of the adjustable load corresponding to the ith prediction moment in the load prediction,/for the load prediction>And->The active power and the reactive power of the first load device at the ith prediction moment are respectively shown.
Step 4, setting a power grid model in the simulation calculation at T m The controllable switch state of +1 day is Cb m,i ：
Cb m,i ＝[i,i＝0,…288]{Cb m,i,name ,Cb m,i,val }；
Wherein Cb is m,i,name 、Cb m,i,val The switch name and the switch state at the ith simulation time are respectively.
The specific implementation method of the steps is as follows:
step 4.1, at T m Cb at time 0:00 of +1 day simulation calculation m,0 State and step 2 controllable switch initial state array Cb 0 The consistency is that:
Cb m,0 ＝Cb 0 ；
step 4.2, at T m At the subsequent time of +1 day simulation calculation, the switch state of the automatic voltage control capacitor reactor in the power grid model is determined according to a simulation control strategy given by AVC in the simulation process, and the state sources of other switches (the switches of the non-automatic voltage control capacitor reactor) are as follows:
step 4.2.1, set at T m The time of day start simulation calculation is T s Corresponding time is I s 。
Step 4.2.2, date T is set m From the i=1 th moment to the i=i s The state of a controllable switch (the switch of the non-automatic voltage control capacitor reactor) of the power grid model at moment is Cb m,val '：
Cb m,val ′＝[i＝1,..I s ]{Cb m,i,val }；
Step 4.2.3, date T is set m From i=i in-1 s The state of the power grid model controllable switch (the switch of the non-automatic voltage control capacitor reactor) from the moment to the i=288 moment is Cb m,val ″：
Cb m,val ″＝[i＝I s ,..288]{Cb m,i,val }；
Step 4.2.4, T is obtained according to steps 4.2.2 and 4.2.3 m +1 day switch (switch of non-automatic voltage controlled capacitive reactor) state information Cb m,i,val The method comprises the following steps:
Cb m,i,val ＝Cb m,val ′(1≤i≤I s )；
Cb m,i,val ＝Cb m,val ″(I s <i≤288)；
step 5, setting the scene information corresponding to the regional power grid simulation case as F sch ：
Let the serial number of the voltage class and the limit value parameter corresponding to the bus bar to be monitored in the power grid model be K, k= … … K, K represents the total number of the voltage class and the limit value parameter corresponding to the bus bar to be monitored in the power grid model.
F sch ＝[k,k＝1,..K]{bus k ,Step k ,Area k ,Type k ,U k,nom ,U k,min ,U k,max }；
Wherein, bus k Step for voltage class and limit value parameter information corresponding to the k-th bus needing to be monitored k Step value, area, of bus limit parameter variable to be monitored for k-th class k Bus limit value parameter effective area to be monitored for k-th class, type k Bus limit value parameter variable type, U, to be monitored for the kth class k,nom Default parameters for the voltage limit parameters of the bus to be monitored for class k, U k,min Voltage lower limit value of bus limit value parameter variable to be monitored for k-th class, U k,max The upper voltage limit value of the bus limit value parameter variable to be monitored for the k-th class;
step 6, setting a simulation task set corresponding to the regional power grid simulation case as T sch According to the scene information F of step 5 sch Bus limit value parameter variable setting information to be monitored to generate a simulation task set T sch . The specific implementation method of the step comprises the following steps:
and 6, step 6.1. Setting task serial numbers P, p= … … P in regional power grid simulation task sets, wherein P is the total number of tasks of the regional power grid simulation task sets, and bus limit value parameter information to be monitored of the kth class corresponding to the kth task in the power grid simulation task sets is Bs k,p ；
Step 6.2, type of bus limit parameter variable to be monitored if the kth class k Is constant, then Bs k,p Medium parameter value Val k,p The default value of the bus limit value parameter variable corresponding to the voltage limit value to be monitored for the k-th class is that:
Val k,p ＝U k,nom ；
step 6.3, type of bus limit parameter variable to be monitored if the kth class k Is a variable, then Bs k,p Medium parameter value Val k,p The combination of the corresponding voltage limits for the bus limit parameter variables that need to be monitored for class k, namely:
Val k,p ＝{U k,min +Step k }(U k,min ≤Val k,p ≤U k,max )；
step 6.4, generating a simulation task set T corresponding to the simulation case according to the steps 6.2 and 6.3 sch :
Setting a sample number p of a voltage limit value parameter combination of a k-th type bus needing to be monitored k ,p k ＝1…P k ，P k The total number of samples of the voltage limit parameter combination of the bus needs to be monitored for class k.
Wherein P is k The total number of voltage limit samples of the bus needs to be monitored for the k-th class, and P is the total number of tasks of the simulation task set.
T sch ＝[p,p＝1,..P]{Bs k,p ,A k,p }；
Wherein Bs is k,p The limit value parameter information of the bus needs to be monitored for the kth class corresponding to the kth task, A k,p The k-th class corresponding to the p-th task requires monitoring of the bus limit parametersEffective area.
Step 7, using the simulation task set T of step 6 sch The first task in (task number p=1) is illustrative of the simulation calculation and index statistics process. The specific implementation method of the step comprises the following steps:
step 7.1, at the 0:00 th moment of the start of the simulation task (task number p=1) calculation, according to the T of step 3 m Regional power grid load prediction information Ld corresponding to +1 day fre From Ld fre Extract the i=0 prediction timeAnd->Data, for the initial array Ld of the power grid load in the step 2.1 0,p And Ld 0,q Updating:
obtaining updated power grid load initial array Ld 0,p 1 And Ld 0,q 1 ：
Ld 0,p 1 ＝[l,l＝1,..L]{P 0,l 1 }；Ld 0,q 1 ＝[l,l＝1,..L]{Q 0,l 1 }；
Step 7.2, date T according to step 4 m Switch state information Cb corresponding to +1 m,i From Cb m,i Extracting the switch state Cb at the i=0 moment m,0 For the initial array Cb of the power grid controllable switch in the step 2.2 0 Updating:
using the switch name Cb m,0,name Retrieving initial array Cb of controllable switch of power grid 0 For matched Cb 0 Switch record uses Cb m,0,val Updating its switch state Open 0,c ：
Open 0,c 1 ＝Cb m,0,val ；
Obtain a betterNew grid switch state array Cb 0 1 ：
Cb 0 1 ＝[c,c＝1,..C]{Open 0,c 1 }
Step 7.3, obtaining an updated initial power grid model F through the steps 7.1 to 7.2 m,0 1 And starting one-time power flow calculation on the basis to obtain a power flow calculation result.
Step 7.4, bus limit parameter information Bs to be monitored according to step 6 k,1 For the initial state array Bs of the voltage limit value of the bus to be monitored in step 2.3 0 Updating:
step 7.4.1 according to Bs k,1 Corresponding scene information bus of step 5 k And Area k Retrieving the bus initial state array Bs to be monitored 0 K in (B) s,i ,A s,i Recording;
step 7.4.2 for Bs 0 Is matched with the bus to be monitored, if the corresponding voltage class and limit value parameter information bus k Belonging to the upper voltage limit, bs are used k,1 Parameter value Val of (2) k,1 Updating array Bs 0 U of corresponding ith point moment in time s,i,max Recording U s,i,min And U s,i,val Remain unchanged, namely:
U s,i,max 1 ＝Val k,1 ；U s,i,min 1 ＝U s,i,min ；U s,i,val 1 ＝U s,i,val ；
step 7.4.3 for Bs 0 Is matched with the bus to be monitored, if the corresponding voltage class and limit value parameter information bus k Belonging to the lower voltage limit, bs are used k,1 Parameter value Val of (2) k,1 Updating array Bs 0 U of corresponding ith point moment in time s,i,min Recording U s,i,max And U s,i,val Remain unchanged, namely:
U s,i,max 1 ＝U s,i,max ；U s,i,min 1 ＝Val k,1 ；U s,i,val 1 ＝U s,i,val ；
step 7.4.4, obtaining an updated voltage limit value array Bs of the bus needing to be monitored through step 7.4.2 and step 7.4.3 0 1 ：
Bs 0 1 ＝[s,s＝1,..S]{K s,i ,A s,i ,U s,i,val 1 ,U s,i,min 1 ,U s,i,max 1 }(0≤i≤288)；
Step 7.5, performing an automatic voltage control simulation calculation according to the tide calculation result in step 7.3 and the voltage limit value information of the bus needing to be monitored updated in step 7.4, and outputting a simulation control strategy A 0,strg ：
A 0,strg ＝[x,x＝0,..X]{Cb 0,x,strg }；
Wherein X is the serial number of the capacitive reactance analog control strategy, x= … … X, X represents the total number of the capacitive reactance analog control strategies, cb 0,x,strg And simulating switching information corresponding to the control strategy for the xth capacitive reactance.
Step 7.6, automatic voltage simulation control strategy A according to step 7.5 0,strg For the initial power grid model F of step 2 m,0 1 Update and utilize updated grid model F m,0 1 ' perform load flow calculation, comprising the steps of:
step 7.6.1, automatic voltage simulation control strategy A from step 7.5 0,strg X < th > of medium reading 0 Simulation control strategy for strip capacitive reactance deviceInitial grid model F at step 2 m,0 1 Middle search switch->In the initial state of the controllable switch, the array Cb 0 1 And the position is marked as c 0 For the initial power grid model F m,0 1 Initial state Cb of medium controllable switch 0 1 Updating: />
Step 7.6.2, traversing automatic voltage simulation control strategy A 0,strg Repeating step 7.6.1 to obtain an updated grid model F m,0 1 ' carrying out one-time power flow calculation to obtain a power flow calculation result F m,0 1″ Recording the load flow calculation result F m,0 1″ Capacitive reactance belonging to analog automatic voltage control model and corresponding switch state of Cb strg,cp An array.
F m,0 1″ ＝{Ld 0,p 1″ ,Ld 0,q 1″ ,Cb 0 1″ ,Bs 0 1″ ,Loss 0 1″ }；
Step 7.7, calculating the result F according to the tide of step 7.6 m,0 1″ The statistics of the regional power grid reactive voltage control index corresponding to the task p=1 in the simulation task set in the step 6 comprises the following steps:
step 7.7.1, setting the serial number of the bus which needs to be monitored for 10kV in the regional power grid model as Y, wherein y= … … Y, Y represents the total number of the bus which needs to be monitored for 10kV in the power grid model, and K10 p,i The voltage out-of-limit statistics of the bus needs to be monitored for 10 kV.
K10 p,i ＝[p＝1,…P]{Up i ,Dn i ,PLoss i }(0≤i≤288)；
Wherein, the subscript p is the task sequence number of the simulation task set in the step 6, the subscript i is the simulation time, up i ,Dn i The upper limit rate and the lower limit rate of the voltage of the bus to be monitored are respectively 10kV corresponding to the ith simulation time, and Loss is calculated i And (5) corresponding to the active loss of the power grid for the ith simulation moment.
Here, the task in the simulation task set p=1 corresponds to the K10 at the i=0 simulation time p,i The array is K10 1,0 。
Analyzing the load flow calculation result F of the step 7.6 m,0 1″ Data, bs according to step 7.4 0 1 Information, obtain Bs 0 1″ Array:
K s,i 1″ ＝K s,i 1 '；A s,i 1″ ＝A s,i 1 '；
U s,i,val 1″ ＝V s,i,val 1″ ；U s,i,min 1″ ＝U s,i,min 1 '；U s,i,max 1″ ＝U s,i,max 1 '；
Bs 0 1″ ＝[s,s＝1,..S]{K s,i 1″ ,A s,i 1″ ,U s,i,val 1″ ,U s,i,min 1″ ,U s,i,max 1″ }(0≤i≤288)；
wherein V is s,i,val 1″ Is F m,0 1″ The voltage value of the bus to be monitored in the s-th bus at the i=0 simulation time.
Step 7.7.2, traversing Bs of step 7.7.1 0 1″ Recording the 10kV bus to be monitored in the array, and counting the upper limit rate information of the voltage of the 10kV bus to be monitored:
wherein,,is Bs 0 1″ The upper limit of the voltage in the bus which needs to be monitored for 10kV in the array is added up.
Step 7.7.3, traversing Bs of step 7.7.1 0 1″ Recording the 10kV bus to be monitored in the array, and counting lower limit rate information of the voltage of the 10kV bus to be monitored:
wherein,,is Bs 0 1″ The lower limit of the voltage in the bus which needs to be monitored by 10kV in the array is added up.
Step 7.7.4 from the initial grid model F m,0 1″ Acquiring power grid active loss information at the i=0 simulation moment:
PLoss 0 ＝Loss 0 1″ ；
step 7.8, simulation task (task number p=1) i=1 simulation time, and load flow calculation result F at step 7.6.2 m,0 1″ Then, model update and simulation calculation at the i=1 time are performed with reference to steps 7.1 to 7.7. Wherein, the power flow calculation result F is calculated at the i=1 simulation moment m,0 1″ When updating, the controllable switch state Cb 0 1″ The switches connected to the capacitive reactors in the array are partially switched to maintain Cb of step 7.6.2 strg,cp And (3) in the running state, monitoring the voltage of the bus and keeping the voltage value obtained by the completion of the tide calculation in the step 7.6.
Step 7.9, traversing all simulation moments in the simulation task (task serial number p=1), repeating the steps 7.1-7.8 to complete simulation calculation and data statistics of the simulation task, and obtaining a power grid reactive voltage control index K10 corresponding to the simulation task 1,i The following are provided:
K10 1,i ＝{Up i ,Dn i ,PLoss i }(0≤i≤288)；
step 8, repeating the step 7 to complete the simulation task set T of the step 6 by adopting a task parallel method through a parallel simulation computing platform sch The simulation calculation of all simulation tasks and the reactive voltage control index statistics of the power grid output the following data:
K10 p,i ＝[p＝1,…P]{Up i ,Dn i ,PLoss i }(0≤i≤288)；
step 9, controlling the reactive voltage control index K10 of the power grid in step 8 p,i The data is analyzed and decided to obtain a group of voltage operation upper and lower limit value parameters of the bus needing to be monitored, which can be used for the regional power grid real-time operation automatic voltage control system, and the method comprises the following steps:
step 9.1, setting a simulation task set T sch The average value of the upper limit rate of the voltage of the 10kV bus corresponding to a single task is Up p,avg ：
Step 9.2, setting a simulation task set T sch The average value of the lower limit rate of the voltage of the 10kV bus corresponding to a single task is Dn p,avg ：
Step 9.3, setting a simulation task set T sch The average value of the active power loss of the power grid corresponding to a single task in the network is PLoss p,avg ：
Step 9.4, obtaining a simulation task set T according to the steps 9.1, 9.2 and 9.3 sch Grid reactive voltage control index average value K10 corresponding to all tasks in the system p,avg Array:
K10 p,avg ＝[p＝1,…P]{Up p,avg ,Dn p,avg ,PLoss p,avg }；
step 9.5, setting a simulation task set T sch The average value of the upper limit rate of the 10kV bus voltage corresponding to all tasks is Up sch,avg ：
Step 9.6, setting a simulation task set T sch The average value of the lower limit rate of the 10kV bus voltage corresponding to all tasks is Dn sch,avg ：
Step 9.7, setting a simulation task set T sch The corresponding optimization task sequence number of the 10kV bus voltage is U, u= … … U, and U represents a simulation task set T sch The total number of optimization tasks of the upper limit rate of the corresponding 10kV bus voltage is recorded as K10 u,opf Array:
K10 u,opf ＝[u＝0,…U]{Up u,opf ,Dn u,opf ,PLoss u,opf }；
retrieving the average value K10 of the reactive power voltage control indexes of the power grid in the step 9.4 p,avg Up in array p,avg Record if Up p,avg Up less than step 9.5 sch,avg Let in K10 u,opf Array:
Up u,opf ＝Up p,avg ；Dn u,opf ＝Dn p,avg ；PLoss u,opf ＝PLoss p,avg ；
step 9.8, setting a simulation task set T sch The corresponding optimization task sequence number of the lower limit rate of the 10kV bus voltage is D, d= … … D, and D represents a simulation task set T sch The total number of optimization tasks with lower limit rate of the corresponding 10kV bus voltage is recorded as K10 d,opf Array:
K10 d,opf ＝[d＝0,…D]{Up d,opf ,Dn d,opf ,PLoss d,opf }；
retrieving the average value K10 of the reactive power voltage control indexes of the power grid in the step 9.4 p,avg Dn in array p,avg Recording, if Dn p,avg Dn less than step 9.5 sch,avg Let in K10 d,opf Array:
Up d,opf ＝Up p,avg ；Dn d,opf ＝Dn p,avg ；PLoss d,opf ＝PLoss p,avg ；
step 9.9, K10 according to step 9.7 u,opf Array and K10 of step 9.8 d,opf An array, obtaining a simulation task set T after taking the intersection sch Corresponding 10kV bus voltage out-of-limit rate optimization task K10 h,opf Array:
K10 h,opf ＝[h＝0,…H]∩{K10 u,opf ,K10 d,opf }；
K10 h,opf ＝[h＝0,…H]{Up h,opf ,Dn h,opf ,PLoss h,opf }；
wherein the subscript h is K10 h,opf Task number H is K10 h,opf Including the total number of tasks.
Step 9.10, K10 from step 9.9 h,opf Retrieving PLoss in array h,opf The optimized task sequence number of the out-of-limit rate of the voltage of the 10kV bus corresponding to the minimum value is h new According to step 6, the power grid simulation task set T sch Information matchingThe simulation task serial number corresponding to the task is p o ,k o From T sch The analysis sequence number in the information is p o ,k o Obtaining an upper limit value array Bs and a lower limit value array Bs of voltage operation of a bus to be monitored opf ：
Step 10, the clock date of the regional power grid is T m At time 0:00 of +1 day, referring to step 7.4, the upper and lower limit value arrays Bs are run using the bus voltage to be monitored of step 9.10 opf Bus voltage operation upper and lower limit value original state Bs to be monitored in automatic voltage control system for real-time operation of regional power grid ora Updating the array to obtain Bs ora new ：
Bs ora new ＝[s＝1,..S]{U i,min new ,U i,max new }(0≤i≤288)；
Step 10.1, measuring a voltage value U of an ith bus needing to be monitored in an ith moment in a real-time operation system of a local power grid i,val Greater than U i,max new When the automatic voltage control system outputs the input of the transformer substation where the s-th bus needing to be monitored isA reactor or capacitor-stripping control strategy to eliminate the upper limit problem of the bus voltage to be monitored in the s th strip;
step 10.2, measuring a voltage value U of an ith bus needing to be monitored in an ith moment in a real-time operation system of a local power grid i,val Less than U i,min new When the voltage of the bus to be monitored is lower than the lower limit, the automatic voltage control system outputs a reactor-removing or capacitor-throwing control strategy of a transformer substation where the bus to be monitored is positioned, so that the problem that the voltage of the bus to be monitored is lower than the lower limit is solved;
step 10.3, strategy of outputting by the automatic voltage control system of step 10.1 and step 10.2, so that the upper and lower limit value Bs of the voltage operation of the bus needs to be monitored opf The array plays a role in automatic voltage control of regional power grid real-time operation. The voltage of the bus needs to be monitored to control the regional power grid to run in a reasonable interval, and meanwhile, the requirements of ensuring the voltage qualification level of the 10kV bus of the regional power grid and reducing the active loss of the power grid are also met.
The following description is made in connection with a specific example in the above-described manner:
(1) Acquiring an initial section of the power grid voltage of a certain region, wherein the initial section is automatically controlled at 2021-07-13 days: in this embodiment, the regional power grid is provided with the following files at time 0:00 of the 2021-07-13 day voltage automatic control history section: avcmodel20210713000000.Dgz.
(2) Load prediction information of a regional power grid in 2021-07-14 days is obtained: in this embodiment, the load forecast file for the regional power grid 2021-07-14 days is: dq_load_20210714.Dat, including information of type, active, reactive etc. corresponding to each load device in the power grid at the predicted time, is exemplified by partial load predicted data at time 0:00 on days 2021-07-14, as follows:
table 2-1 regional grid load prediction information
(3) Acquiring controllable switch state information of a regional power grid on days 2021-07-14 of simulation calculation: in this embodiment, the switching state file of the regional power grid on days 2021-07-14 is cb_stat_20210714.Dat, including information such as types and switching states corresponding to switching devices in the power grid at simulation time, and the partial switching states at the times 0:00 on days 2021-07-14 are illustrated as the following table:
TABLE 3-1 regional grid switch status information
(4) Setting simulation scene information of a regional power grid: in the embodiment, the simulation scene selects 4 limit parameter variables of the upper limit value of the 220kV bus voltage, the lower limit value of the 220kV bus voltage, the upper limit value of the 110kV bus voltage and the lower limit value of the 110kV bus voltage of the regional power grid, and the following table is adopted:
TABLE 4-1 simulation scenario information
(5) Generating a simulation task set corresponding to a power grid simulation scene in a certain region: according to the limit parameter variable information contained in the simulation scene definition, a simulation task set of the regional power grid is generated through permutation and combination of limit parameter variable value samples, and the simulation task set generated in the embodiment contains 169 simulation tasks.
220kV bus voltage upper limit parameter value: ranges [228-234], steps 0.5, totaling 13.
The lower limit parameter of the 220kV bus voltage is valued: ranges [224-224], constant, 1 total.
The upper limit parameter of the voltage of the 110kV bus is as follows: ranges [118-124], step sizes 0.5, totaling 13.
The lower limit parameter of the voltage of the 110kV bus is valued: the range [110-110], constant, total 1.
In summary, the total number of the valued combinations of 4 parameters in the simulation scene definition is as follows: 13 x 1 x 13 x 1 = 169, the resulting simulation task set is as follows:
table 5-1 simulation task set information (Unit kV)
/>
/>
Starting a parallel simulation platform to start simulation calculation and index statistics: after the simulation task set is generated, the parallel simulation platform is utilized to start simulation calculation and index statistics, and in the embodiment, the calculation resource node of the parallel simulation platform server is 25, so that the simulation calculation of 25 tasks can be simultaneously executed.
(6) Taking a 110kV east gateway station corresponding to a first task in a simulation task set as an example, a data updating flow in simulation calculation is described. At the 0:00 moment of simulation calculation, the 110kV east-note station in the regional power grid model needs to update the following data:
(6-1) updating the load data of the 110kV east-note station at the time of 0:00 according to the load prediction data of the step 2.
Table 6-1 eastern gateway station 0:00 time load prediction information
| Time | Load name | Type(s) | Active power | Reactive power |
| 2021/7/140:00 | 110kV Dongguan station Dongaluminum wire-394 | 4 | 0 | 0 |
| 2021/7/140:00 | 110kV east-closing station east combustion line-397 | 4 | 0 | 0 |
| 2021/7/140:00 | 110kV east-note line-398 | 4 | 0 | 0 |
| 2021/7/140:00 | 110kV Dongguan station Dong Sun Xian-195 | 4 | 0 | 0 |
| 2021/7/140:00 | Town IV return-592 at 110kV east Guangdong station | 4 | 2.049 | 0.4685 |
| 2021/7/140:00 | Town III return-596 | 4 | 0 | 0 |
| 2021/7/140:00 | Town I return-597 | 4 | 0.7116 | 0.187 |
| 2021/7/140:00 | 110kV east Guangdong station town II return-891 | 4 | 0 | 0 |
| 2021/7/140:00 | 110kV Dongguan station Yangjia Bay-892 | 4 | 0.6738 | 0.2219 |
| 2021/7/140:00 | 110kV east Guangdong station.10 kV town V return-893 | 4 | 0 | 0 |
| 2021/7/140:00 | 110kV east-Guangdong station, hometown line-895 | 4 | 3.05 | 1.0497 |
| 2021/7/140:00 | 110kV Dongguan station rhubarb slope line-896 | 4 | 0 | 0 |
| 2021/7/140:00 | Xinbao high-speed line-591 | 4 | 0 | 0 |
| 2021/7/140:00 | 110kV east Guangdong station, new City I loop-593 | 4 | 0.2598 | 0.0942 |
| 2021/7/140:00 | 110kV east Guangdong station New City II loop-894 | 4 | 0.6523 | 0.0291 |
(6-2) updating the switching state of the 110kV east-note station at the moment of 0:00 according to the switching state data of the step 3.
TABLE 6-2 east-note station 0:00 moment switch status information
/>
/>
(6-3) according to the simulation task set data of the step 5, the middle bus voltage limit value corresponding to the task one is as follows:
Table 6-3 simulation task A corresponding bus voltage limit information (Unit kV)
| 220kV bus voltage upper limit | 220kV bus voltage lower limit | 110kV bus voltage upper limit | 110kV bus voltage lower limit |
| 234 | 224 | 124 | 110 |
According to bus voltage limit information of the table (6-3), the lower voltage limit and the upper voltage limit of the 110kV bus of the 110kV east-west station need to be updated are respectively as follows: 110kV and 124kV.
And (6-4) updating data corresponding to simulation time by referring to a method of a 110kV east-close station by other stations in the regional power grid, starting one-time power flow calculation and automatic voltage control simulation calculation, adjusting the corresponding capacitive reactor switch state in the power grid model according to the output simulation control strategy, and counting reactive voltage control indexes corresponding to the simulation time.
And (6-5) repeating the processes from the step 6-1) to the step 6-4) at the subsequent simulation moment of the first task to update the power grid data, calculate the tide, simulate and calculate the automatic voltage control and calculate the reactive voltage control index until the simulation calculation at all moments is completed.
(6-6) referring to the task I, completing simulation calculation of all tasks in a simulation task set by using a parallel simulation platform, wherein the simulation calculation is as follows:
TABLE 6-4 parallel simulation calculation task State
/>
/>
/>
/>
Obtaining a power grid reactive voltage statistical index corresponding to the simulation task set:
In this embodiment, the reactive voltage statistics index of the power grid is: the upper limit rate of the 10kV bus voltage, the lower limit rate of the 10kV bus voltage and the active loss of the system are calculated, and after each simulation task in the simulation task set is completed, the corresponding statistical indexes are output as follows:
table 7-1 simulation task set grid reactive voltage statistical index
/>
/>
Analysis decision for reactive voltage statistical index of power grid
In this embodiment, according to the data analysis simulation task set in step 7, the average value of the overall 10kV busbar voltage over the upper limit rate is: 0.040; the average value of the upper limit rate of the 10kV bus voltage is as follows: 0.082; the average value of the active loss of the system is as follows: 8.045MW.
Obtaining voltage limit parameter information for a regional power grid 2021-07-14 days, where bus monitoring is required
In this embodiment, on the basis of the analysis result in step 8, the task number in step 7, which satisfies two conditions that the voltage threshold ratio of the 10kV bus is lower than the overall average value and the active loss of the system is at least, is 25, and the simulation task set list in step 5 is searched to obtain the voltage operation limit value parameters of the bus to be monitored, which are:
220kV bus voltage upper operation limit: 233.5kV; the range is as follows: a full network;
220kV bus voltage operation lower limit: 224.0kV; the range is as follows: a full network;
Upper limit of 110kV bus voltage operation: 118.5kV; the range is as follows: a full network;
110kV bus voltage operation lower limit: 110.0kV; the range is as follows: a full network;
AVC control strategy after 2021-07-14 days of updating bus voltage limit value of power grid in certain region
In this embodiment, at time 0:00 of regional power grid 2021-07-14, the voltage operation limit parameter obtained in step 9) is updated to an automatic voltage control system operating in real time, and the generation process and execution effect of AVC control strategy are illustrated by taking the power grid operation state at time 12:32:24 of 2021-07-14 as an example.
At this time, the lower limit (lower limit value [224.0 ]) of the 220kV south parent voltage measurement value of 220kV in the county-level station is 223.14kV, the AVC needs to be subjected to reactor withdrawal or capacitor throwing to [ increase reactive power ], and after the running states of the reactor and the capacitor equipment in the transformer substation are searched, the 220kV county-level station [ 774- - -capacitor ] is selected to be thrown to solve the problem of the 220kV south parent voltage out-of-limit of the 220kV county-level station, and the specific strategy is as follows:
control unit state: bus [220kV code county station.220 kV south bus ]: [223.06/223.20/223.14] _ upper limit [233.5] _ lower limit [224.0]; bus [220kV code county station.110 kV south bus ]: [114.27/114.27/114.26] _ upper limit [118.5] _ lower limit [110.0]; bus [220kV code county station.35 kV II section bus ]: [36.49/36.50/36.49] _ upper limit [38.5] _ lower limit [35.0].
The [1] st scanning, the reactor runs [00] and the capacitor can be alternatively switched [04 ].
[220kV Dynasty station.774- - -capacitor ] out of limit [ general 0] [ serious 0] after input, voltage estimation is that [115.04] in high [226.37] is low [36.99], and can be input.
[220kV Dynasty station.775 capacitor ] is out of limit [ general 0] [ serious 0] after being put into operation, and the voltage is estimated to be [115.04] low [36.99] in high [226.37], so that the capacitor can be put into operation.
[220kV Dynasty station.771- - -capacitor ] out of limit [ general 0] [ serious 0] after input, voltage estimation is that [115.04] in high [226.37] is low [36.99], and can be input.
[220kV Dynasty station.773- - -capacitor ] out of limit [ general 0] [ serious 0] after input, voltage estimation is that [115.04] in high [226.37] is low [36.99], and can be input.
[220kV county station.774- - -capacitor ] operating factor= [12199999], [220kV county station.775 capacitor ] operating factor= [12199999], [220kV county station.771- - -capacitor ] operating factor= [12199999], [220kV county station.773- - -capacitor ] operating factor= [12199999].
And selecting and putting in [220kV county station.774- - -capacitor ], and generating a control strategy [1] of the station [220kV county station ].
After the AVC strategy of the 220kV code county station is successfully executed, the voltage measurement of the high-middle-low three-side bus of the 220kV code county station is as follows: 220kV for the measurement of the south parent voltage of 220kV for the county station: 226.19kV;220kV code county station.110 kV south bus voltage measurement: 114.46kV;220kV code county station.35 kV II section busbar voltage measurement: 36.57kV.
It should be emphasized that the examples described herein are illustrative rather than limiting, and therefore the invention includes, but is not limited to, the examples described in the detailed description, as other embodiments derived from the technical solutions of the invention by a person skilled in the art are equally within the scope of the invention.
Claims (9)
1. An automatic voltage control method for regional power grids based on parallel simulation is characterized by comprising the following steps of: the method comprises the following steps:
step 1, setting an automatic voltage control period T of a regional power grid c ；
Step 2, setting the date T of the current clock of the regional power grid m Reading the date T from a historical section of regional power grid voltage automatic control m Corresponding power grid model M and initial section F at time 0:00 m,0 ；
Step 3, setting the regional power grid clock date as T m Load prediction information Ld for +1 day fre ；
Step 4, setting a power grid model in the simulation calculation at T m Controllable switch state Cb of +1 day m,i ：
Step 5, setting scene information F corresponding to regional power grid simulation cases sch ：
Step 6, setting a simulation task set T corresponding to the regional power grid simulation case sch According to the scene information F of step 5 sch Bus limit value parameter variable setting information to be monitored to generate a simulation task set T sch ；
Step 7, for the simulation task set T of step 6 sch Performing simulation calculation and index statistics on one task of the system;
step 8, repeating the step 7 by adopting a task parallel method to finish the step6 simulation task set T sch The simulation calculation of all simulation tasks and the statistics of reactive power voltage control indexes of the power grid output the reactive power voltage control index K10 of the power grid p,i ；
Step 9, controlling the reactive voltage control index K10 of the power grid in step 8 p,i The data are analyzed and decided to obtain a group of voltage operation upper and lower limit value parameters which can be used for the regional power grid real-time operation automatic voltage control system and need to monitor buses;
step 10, the clock date of the regional power grid is T m At the time of 0:00 of +1 day, the original state Bs of the upper limit value and the lower limit value of the voltage operation of the bus to be monitored in the automatic voltage control system for real-time operation of the regional power grid is obtained by utilizing the upper limit value and the lower limit value parameters of the voltage operation of the bus to be monitored ora Updating the array;
the specific implementation method of the step 7 is as follows:
step 7.1, calculating and starting the simulation task with the task serial number p=1 at the time of 0:00, and according to T in the step 3 m Regional power grid load prediction information Ld corresponding to +1 day fre From Ld fre Extract the i=0 prediction timeAnd->Data, for the initial array Ld of the power grid load in the step 2.1 0,p And Ld 0,q Updating:
Obtaining updated power grid load initial array Ld 0,p 1 And Ld 0,q 1 ：
Ld 0,p 1 ＝l,l＝1,..LP 0,l 1 ；Ld 0,q 1 ＝l,l＝1,..LQ 0,l 1 ；
Step 7.2 according to the steps of4 date T m Switch state information Cb corresponding to +1 m,i From Cb m,i Extracting the switch state Cb at the i=0 moment m,0 For the initial array Cb of the power grid controllable switch in the step 2.2 0 Updating:
using the switch name Cb m,0,name Retrieving initial array Cb of controllable switch of power grid 0 For matched Cb 0 Switch record uses Cb m,0,val Updating its switch state Open 0,c ：
Open 0,c 1 ＝Cb m,0,val ；
Obtaining updated power grid switch state array Cb 0 1 ：
Cb 0 1 ＝[c,c＝1,..C]{Open 0,c 1 }
Step 7.3, obtaining an updated initial power grid model F through the steps 7.1 to 7.2 m,0 1 Starting one-time power flow calculation on the basis to obtain a power flow calculation result;
step 7.4, bus limit parameter information Bs to be monitored according to step 6 k,1 For the initial state array Bs of the voltage limit value of the bus to be monitored in step 2.3 0 Updating;
step 7.4.1 according to Bs k,1 Corresponding scene information bus of step 5 k And Area k Retrieving the bus initial state array Bs to be monitored 0 K in (B) s,i ,A s,i Recording;
step 7.4.2 for Bs 0 Is matched with the bus to be monitored, if the corresponding voltage class and limit value parameter information bus k Belonging to the upper voltage limit, bs are used k,1 Parameter value Val of (2) k,1 Updating array Bs 0 U of corresponding ith point moment in time s,i,max Recording U s,i,min And U s,i,val Remain unchanged, namely:
U s,i,max 1 ＝Val k,1 ；U s,i,min 1 ＝U s,i,min ；U s,i,val 1 ＝U s,i,val ；
step 7.4.3 for Bs 0 Is matched with the bus to be monitored, if the corresponding voltage class and limit value parameter information bus k Belonging to the lower voltage limit, bs are used k,1 Parameter value Val of (2) k,1 Updating array Bs 0 U of corresponding ith point moment in time s,i,min Recording U s,i,max And U s,i,val Remain unchanged, namely:
U s,i,max 1 ＝U s,i,max ；U s,i,min 1 ＝Val k,1 ；U s,i,val 1 ＝U s,i,val ；
step 7.4.4, obtaining an updated voltage limit value array Bs of the bus needing to be monitored through step 7.4.2 and step 7.4.3 0 1 ：
Bs 0 1 ＝[s,s＝1,..S]{K s,i ,A s,i ,U s,i,val 1 ,U s,i,min 1 ,U s,i,max 1 }(0≤i≤288)；
Step 7.5, performing an automatic voltage control simulation calculation according to the tide calculation result in step 7.3 and the voltage limit value information of the bus needing to be monitored updated in step 7.4, and outputting a simulation control strategy A 0,strg ：
A 0,strg ＝[x,x＝0,..X]{Cb 0,x,strg }；
Wherein X is the serial number of the capacitive reactance analog control strategy, x= … … X, X represents the total number of the capacitive reactance analog control strategies, cb 0,x,strg Switching information corresponding to the simulation control strategy of the capacitor x;
step 7.6, automatic voltage simulation control strategy A according to step 7.5 0,strg For the initial power grid model F of step 2 m,0 1 Update and utilize updated grid model F m,0 1 ' perform load flow calculation, comprising the steps of:
step 7.6.1, automatic voltage simulation control strategy A from step 7.5 0,strg X < th > of medium reading 0 Simulation control strategy for strip capacitive reactance deviceInitial grid model F at step 2 m,0 1 Middle search switch->In the initial state of the controllable switch, the array Cb 0 1 And the position is marked as c 0 For the initial power grid model F m,0 1 Initial state Cb of medium controllable switch 0 1 Updating:
step 7.6.2, traversing automatic voltage simulation control strategy A 0,strg Repeating step 7.6.1 to obtain an updated grid model F m,0 1 ' carrying out one-time power flow calculation to obtain a power flow calculation result F m,0 1 ", record the load flow calculation result F m,0 1 Capacitive reactance belonging to analog automatic voltage control model and corresponding switch state of capacitive reactance being Cb strg,cp An array;
F m,0 1 ”＝{Ld 0,p 1” ,Ld 0,q 1” ,Cb 0 1” ,Bs 0 1” ,Loss 0 1” }；
step 7.7, calculating the result F according to the tide of step 7.6 m,0 1 The regional power grid reactive voltage control index corresponding to the task of p=1 in the simulation task set of the step 6 is counted, and the method comprises the following steps:
step 7.7.1, setting the serial number of the bus which needs to be monitored for 10kV in the regional power grid model as Y, wherein y= … … Y, Y represents the total number of the bus which needs to be monitored for 10kV in the power grid model, and K10 p,i Voltage out-of-limit statistics of the bus needs to be monitored for 10 kV:
K10 p,i ＝[p＝1,…P]{Up i ,Dn i ,PLoss i }(0≤i≤288)；
wherein, the subscript p is the task sequence number of the simulation task set in the step 6, the subscript i is the simulation time, up i ,Dn i The upper limit rate and the lower limit rate of the voltage of the bus to be monitored are respectively 10kV corresponding to the ith simulation time, and Loss is calculated i The i simulation time corresponds to the active loss of the power grid;
corresponding to the task in the simulation task set p=1, at the i=0 simulation time, K10 p,i The array is K10 1,0 ；
Analyzing the load flow calculation result F of the step 7.6 m,0 1 "data, bs according to step 7.4 0 1 Information, obtain Bs 0 1” Array:
K s,i 1” ＝K s,i 1' ；
A s,i 1” ＝A s,i 1' ；
U s,i,val 1” ＝V s,i,val 1” ；
U s,i,min 1” ＝U s,i,min 1' ；
U s,i,max 1” ＝U s,i,max 1' ；
Bs 0 1” ＝[s,s＝1,..S]{K s,i 1” ,A s,i 1” ,U s,i,val 1” ,U s,i,min 1” ,U s,i,max 1” }(0≤i≤288)；
wherein V is s,i,val 1 "is F m,0 1 "the voltage value of the bus to be monitored in the s-th strip at the i=0 simulation moment;
step 7.7.2, traversing Bs of step 7.7.1 0 1 Recording of 10kV to monitor buses in an array, and counting upper limit rate information of the voltage of the 10kV to monitor buses:
wherein,,is Bs 0 1 The total sum of the upper limit voltage in the buses which are 10kV in the array and need to be monitored is calculated;
step 7.7.3, traversing Bs of step 7.7.1 0 1 Recording of 10kV required monitoring buses in an array, and counting lower limit rate information of the voltage of the 10kV required monitoring buses:
wherein,,is Bs 0 1 The total sum of the lower limit voltage in the buses which are 10kV in the array and need to be monitored is calculated;
step 7.7.4 from the initial grid model F m,0 1 "acquiring power grid active loss information at the i=0 simulation moment:
PLoss 0 ＝Loss 0 1 ”；
step 7.8, simulation task with task number p=1, i=1 simulation time, and load flow calculation result F in step 7.6.2 m,0 1 On the basis, model updating and simulation calculation at the i=1 moment are carried out according to the steps 7.1 to 7.7; wherein, the power flow calculation result F is calculated at the i=1 simulation moment m,0 1 "when updating, controllable switch state Cb 0 1 "the switches connected to the capacitive reactors in the array part need to maintain Cb of step 7.6.2 strg,cp The operation state, the voltage of the monitoring bus keeps the voltage value of the completion of the tide calculation in the previous step 7.6;
step 7.9, traversing all simulation moments in the simulation task with the task sequence number p=1, and repeating the steps 7.1-7.8 to complete simulation calculation and data statistics of the simulation task to obtain a power grid reactive voltage control index K10 corresponding to the simulation task 1,i The following are provided:
K10 1,i ＝Up i ,Dn i ,PLoss i 0≤i≤288。
2. the regional power grid automatic voltage control method based on parallel simulation according to claim 1, wherein the method comprises the following steps of: the initial section F read in the step 2 m,0 The method comprises the following steps:
F m,0 ＝Ld 0,p ,Ld 0,q ,Cb 0 ,Bs 0 ,Loss 0 ；
wherein the subscript m corresponds to the date T m ，Loss 0 Active loss of the initial section;
initial section F m,0 The setting method is as follows:
step 2.1, setting the serial number of the adjustable load in the regional power grid as L, wherein l= … … L, L represents the total number of the adjustable loads in the power grid model M, and Ld 0,p And Ld 0,q For an initial active and initial reactive of the adjustable load at time 0:00:
Ld 0,p ＝l,l＝1,..LP 0,l ；
Ld 0,q ＝l,l＝1,..LQ 0,l ；
wherein P is 0,l And Q 0,l The initial active power and the initial reactive power of the first adjustable load at the moment 0:00 respectively;
step 2.2, setting the serial number of the controllable switches in the regional power grid as C, c= … … C, wherein C represents the total number of the controllable switches in the power grid model M, cb 0 An initial state array for the controllable switch:
Cb 0 ＝c,c＝1,..COpen 0,c ；
wherein Open is 0,c The switching state of the c-th controllable switch at the moment 0:00 is as follows:
step 2.3, setting the serial number of the bus needing to be monitored in the regional power grid as S, wherein s= … … S, and S represents the total number of the bus needing to be monitored in the power grid model M, and Bs 0 Initial state array for bus voltage limit value to be monitored:
Bs 0 ＝s,s＝1,..SK s,i ,A s,i ,U s,i,val ,U s,i,min ,U s,i,max 0≤i≤288；
wherein K is s,i The voltage grade information corresponding to the s-th adjustable mother is obtained; a is that s,i The information of the region corresponding to the bus needing to be monitored is the s-th bus; u (U) s,i,val The voltage value information of the bus needing to be monitored at the ith moment is the s th bus; u (U) s,i,min And U s,i,max The operation upper limit value and the operation lower limit value of the bus needing to be monitored at the ith moment are respectively the s th bus.
3. The regional power grid automatic voltage control method based on parallel simulation according to claim 1, wherein the method comprises the following steps of: the load prediction information set in the step 3 is Ld fre The method comprises the following steps:
wherein the subscript i is the load prediction time,respectively, the active and reactive information of the adjustable load corresponding to the ith prediction moment in the load prediction,/for the load prediction>And->The active power and the reactive power of the first load device at the ith prediction moment are respectively shown.
4. According to claimThe regional power grid automatic voltage control method based on parallel simulation described in claim 1 is characterized in that: the controllable switch state Cb set in the step 4 m,i The method comprises the following steps:
Cb m,i ＝i,i＝0,…288Cb m,i,name ,Cb m,i,val ；
wherein Cb is m,i,name 、Cb m,i,val The switch name and the switch state at the ith simulation moment are respectively;
the specific setting method comprises the following steps:
step 4.1, at T m Cb at time 0:00 of +1 day simulation calculation m,0 State and step 2 controllable switch initial state array Cb 0 And (3) coincidence: cb (Cb) m,0 ＝Cb 0 ；
Step 4.2, at T m At the subsequent time of +1 day simulation calculation, the switching state of the automatic voltage control capacitor reactor in the power grid model is determined according to a simulation control strategy given by AVC in the simulation process, and the switching state sources of the non-automatic voltage control capacitor reactor are as follows:
step 4.2.1, set at T m The time of day start simulation calculation is T s Corresponding time is I s ；
Step 4.2.2, date T is set m From the i=1 th moment to the i=i s The controllable switch state of the power grid model at moment is Cb m,val ′：
Cb m,val ′＝[i＝1,..I s }{Cb m,i,val }；
Step 4.2.3, date T is set m From i=i in-1 s The controllable switch state of the power grid model from moment to i=288 moment is Cb m,val ”：
Cb m,val ″＝[i＝I s ,..288]{Cb m,i,val }；
Step 4.2.4, T is obtained according to steps 4.2.2 and 4.2.3 m +1 day switch state information Cb m,i,val The method comprises the following steps:
Cb m,i,val ＝Cb m,val ′(1≤i≤I s )；
Cb m,i,val ＝Cb m,val ″(I s <i≤288)。
5. the regional power grid automatic voltage control method based on parallel simulation according to claim 1, wherein the method comprises the following steps of: the specific implementation method of the step 5 is as follows: let the serial number of the voltage class and the limit value parameter corresponding to the bus to be monitored in the power grid model be K, k= … … K, where K represents the total number of the voltage class and the limit value parameter corresponding to the bus to be monitored in the power grid model, and then the scene information F corresponding to the regional power grid simulation case sch The method comprises the following steps:
F sch ＝[k,k＝1,..K]{bus k ,Step k ,Area k ,Type k ,U k,nom ,U k,min ,U k,max }；
wherein, bus k Step for voltage class and limit value parameter information corresponding to the k-th bus needing to be monitored k Step value, area, of bus limit parameter variable to be monitored for k-th class k Bus limit value parameter effective area to be monitored for k-th class, type k Bus limit value parameter variable type, U, to be monitored for the kth class k,nom Default parameters for the voltage limit parameters of the bus to be monitored for class k, U k,min Voltage lower limit value of bus limit value parameter variable to be monitored for k-th class, U k,max The upper voltage limit of the bus limit parameter variable that needs to be monitored for class k.
6. The regional power grid automatic voltage control method based on parallel simulation according to claim 1, wherein the method comprises the following steps of: the specific implementation method of the step 6 is as follows:
step 6.1, setting task serial numbers P in regional power grid simulation task sets, wherein p= … … P, P is the total number of tasks in the regional power grid simulation task sets, and the k-th bus limit value parameter information to be monitored corresponding to the P-th task in the power grid simulation task sets is Bs k,p ；
Step 6.2, type of bus limit parameter variable to be monitored if the kth class k Is constant, then Bs k,p Medium parameter value Val k,p The default value of the bus limit value parameter variable corresponding to the voltage limit value to be monitored for the k-th class is that:
Val k,p ＝U k,nom ；
Step 6.3, type of bus limit parameter variable to be monitored if the kth class k Is a variable, then Bs k,p Medium parameter value Val k,p The combination of the corresponding voltage limits for the bus limit parameter variables that need to be monitored for class k, namely:
Val k,p ＝U k,min +Step k U k,min ≤Val k,p ≤U k,max ；
step 6.4, generating a simulation task set T corresponding to the simulation case according to the steps 6.2 and 6.3 sch ：
Setting a sample number p of a voltage limit value parameter combination of a k-th type bus needing to be monitored k ,p k ＝1…P k ，P k The total number of samples of the voltage limit parameter combination of the bus bar that need to be monitored for the k-th class;
wherein P is k The total number of voltage limit value samples of the bus to be monitored is k, and P is the total number of tasks of the simulation task set;
T sch ＝p,p＝1,..PBs k,p ,A k,p ；
wherein Bs is k,p The limit value parameter information of the bus needs to be monitored for the kth class corresponding to the kth task, A k,p The k-th class corresponding to the p-th task needs to monitor the effective area of the limit value parameters of the bus.
7. The regional power grid automatic voltage control method based on parallel simulation according to claim 1, wherein the method comprises the following steps of: the reactive voltage control index K10 of the power grid output in the step 8 p,i The method comprises the following steps:
K10 p,i ＝p＝1,…PUp i ,Dn i ,PLoss i 0≤i≤288。
8. the regional power grid automatic voltage control method based on parallel simulation according to claim 1, wherein the method comprises the following steps of: the specific implementation method of the step 9 comprises the following steps:
Step 9.1, setting a simulation task set T sch The average value of the upper limit rate of the voltage of the 10kV bus corresponding to a single task is Up p,avg ：
Step 9.2, setting a simulation task set T sch The average value of the lower limit rate of the voltage of the 10kV bus corresponding to a single task is Dn p,avg ：
Step 9.3, setting a simulation task set T sch The average value of the active power loss of the power grid corresponding to a single task in the network is PLoss p,avg ：
Step 9.4, obtaining a simulation task set T according to the steps 9.1, 9.2 and 9.3 sch Grid reactive voltage control index average value K10 corresponding to all tasks in the system p,avg Array:
K10 p,avg ＝p＝1,…PUp p,avg ,Dn p,avg ,PLoss p,avg ；
step 9.5, setting a simulation task set T sch The average value of the upper limit rate of the 10kV bus voltage corresponding to all tasks is Up sch,avg ：
Step 9.6, setting a simulation task set T sch The average value of the lower limit rate of the 10kV bus voltage corresponding to all tasks is Dn sch,avg ：
Step 9.7, setting a simulation task set T sch The corresponding optimization task sequence number of the 10kV bus voltage is U, u= … … U, and U represents a simulation task set T sch The total number of optimization tasks of the upper limit rate of the corresponding 10kV bus voltage is recorded as K10 u,opf Array:
K10 u,opf ＝u＝0,…UUp u,opf ,Dn u,opf ,PLoss u,opf ；
retrieving the average value K10 of the reactive power voltage control indexes of the power grid in the step 9.4 p,avg Up in array p,avg Record if Up p,avg Up less than step 9.5 sch,avg Let in K10 u,opf Array:
Up u,opf ＝Up p,avg ；
Dn u,opf ＝Dn p,avg ；
PLoss u,opf ＝PLoss p,avg ；
Step 9.8, setting a simulation task set T sch The corresponding optimization task sequence number of the lower limit rate of the 10kV bus voltage is D, d= … … D, and D represents a simulation task set T sch The total number of optimization tasks with lower limit rate of the corresponding 10kV bus voltage is recorded as K10 d,opf Array:
K10 d,opf ＝d＝0,…DUp d,opf ,Dn d,opf ,PLoss d,opf ；
retrieving the average value K10 of the reactive power voltage control indexes of the power grid in the step 9.4 p,avg Dn in array p,avg Recording, if Dn p,avg Dn less than step 9.5 sch,avg Let in K10 d,opf Array:
Up d,opf ＝Up p,avg ；
Dn d,opf ＝Dn p,avg ；
PLoss d,opf ＝PLoss p,avg ；
step 9.9, K10 according to step 9.7 u,opf Array and K10 of step 9.8 d,opf An array, obtaining a simulation task set T after taking the intersection sch Corresponding 10kV bus voltage out-of-limit rate optimization task K10 h,opf Array:
K10 h,opf ＝h＝0,…H∩K10 u,opf ,K10 d,opf ；
K10 h,opf ＝h＝0,…HUp h,opf ,Dn h,opf ,PLoss h,opf ；
wherein the subscript h is K10 h,opf Task number H is K10 h,opf The total number of tasks included;
step 9.10, K10 from step 9.9 h,opf Retrieving PLoss in array h,opf The optimized task sequence number of the out-of-limit rate of the voltage of the 10kV bus corresponding to the minimum value is h new According to step 6, the power grid simulation task set T sch Information matchingThe simulation task serial number corresponding to the task is p o ,k o From T sch The analysis sequence number in the information is p o ,k o Obtaining an upper limit value array Bs and a lower limit value array Bs of voltage operation of a bus to be monitored opf ：
9. The regional power grid automatic voltage control method based on parallel simulation according to claim 1, wherein the method comprises the following steps of: the step 10 is updated to obtain new Bs ora new The array is as follows:
Bs ora new ＝s＝1,..SU i,min new ,U i,max new 0≤i≤288；
the updating method comprises the following steps:
step 10.1, measuring a voltage value U of an ith bus needing to be monitored in an ith moment in a real-time operation system of a local power grid i,val Greater than U i,max new When the voltage of the bus to be monitored exceeds the upper limit, the automatic voltage control system outputs a reactor throwing or capacitor withdrawing control strategy of a transformer substation where the bus to be monitored is positioned so as to eliminate the problem that the voltage of the bus to be monitored exceeds the upper limit;
step 10.2, measuring a voltage value U of an ith bus needing to be monitored in an ith moment in a real-time operation system of a local power grid i,val Less than U i,min new When the voltage of the bus to be monitored is lower than the lower limit, the automatic voltage control system outputs a reactor-removing or capacitor-throwing control strategy of a transformer substation where the bus to be monitored is positioned, so that the problem that the voltage of the bus to be monitored is lower than the lower limit is solved;
step 10.3, strategy of outputting by the automatic voltage control system of step 10.1 and step 10.2, so that the upper and lower limit value Bs of the voltage operation of the bus needs to be monitored opf The array plays a role in automatic voltage control of regional power grid real-time operation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111517126.2A CN114336598B (en) | 2021-12-13 | 2021-12-13 | Regional power grid automatic voltage control method based on parallel simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111517126.2A CN114336598B (en) | 2021-12-13 | 2021-12-13 | Regional power grid automatic voltage control method based on parallel simulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114336598A CN114336598A (en) | 2022-04-12 |
| CN114336598B true CN114336598B (en) | 2023-07-28 |
Family
ID=81051563
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111517126.2A Active CN114336598B (en) | 2021-12-13 | 2021-12-13 | Regional power grid automatic voltage control method based on parallel simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114336598B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1610541A1 (en) * | 1988-08-09 | 1990-11-30 | Омский Институт Инженеров Железнодорожного Транспорта | Method of controlling capacitive compensation power |
| CN111555293A (en) * | 2020-04-24 | 2020-08-18 | 国网河北省电力有限公司电力科学研究院 | Reactive power equipment type selection system and method based on data model analysis |
| CN112865108A (en) * | 2021-01-11 | 2021-05-28 | 国网山西省电力公司忻州供电公司 | Power grid automatic voltage control simulation method based on continuous power flow simulation |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111564849B (en) * | 2020-05-15 | 2021-11-02 | 清华大学 | Two-stage deep reinforcement learning-based power grid reactive voltage control method |
-
2021
- 2021-12-13 CN CN202111517126.2A patent/CN114336598B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1610541A1 (en) * | 1988-08-09 | 1990-11-30 | Омский Институт Инженеров Железнодорожного Транспорта | Method of controlling capacitive compensation power |
| CN111555293A (en) * | 2020-04-24 | 2020-08-18 | 国网河北省电力有限公司电力科学研究院 | Reactive power equipment type selection system and method based on data model analysis |
| CN112865108A (en) * | 2021-01-11 | 2021-05-28 | 国网山西省电力公司忻州供电公司 | Power grid automatic voltage control simulation method based on continuous power flow simulation |
Non-Patent Citations (3)
| Title |
|---|
| 主动负荷参与配电网分布式智能电压控制;徐弢;李天楚;郭凌旭;林军;张世宇;王笑雪;;电力系统及其自动化学报(第01期);全文 * |
| 基于实时仿真系统的地区电网AVC系统闭环检测方法;唐永红;徐琳;范宏;丁会凯;;电网技术(第09期);全文 * |
| 大规模城市电网基于统一无功优化的自动电压控制;鄂志君;张宝栋;刘伟;汤磊;郭章喜;;电力系统及其自动化学报(第S1期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114336598A (en) | 2022-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108876163B (en) | Transient state power angle stability rapid evaluation method integrating causal analysis and machine learning | |
| CN109409770A (en) | A kind of rural power grids level of intelligence evaluation method neural network based | |
| CN114142530A (en) | Risk scheduling method considering N-1 security constraint based on near-end policy optimization algorithm | |
| CN104008279B (en) | Method for solving power network static security domain | |
| CN108270216B (en) | Multi-target-considered complex power distribution network fault recovery system and method | |
| CN115688993A (en) | Short-term power load prediction method suitable for power distribution station area | |
| CN114336690B (en) | Intelligent terminal-based three-phase imbalance management method, system and storage medium for transformer area | |
| CN106296002A (en) | A kind of power distribution network energy-saving benefit evaluating method based on life cycle management model | |
| CN111193295A (en) | Distribution network flexibility improvement robust optimization scheduling method considering dynamic reconfiguration | |
| Yin et al. | Sequential reconfiguration of unbalanced distribution network with soft open points based on deep reinforcement learning | |
| CN114336598B (en) | Regional power grid automatic voltage control method based on parallel simulation | |
| CN108694475B (en) | Short-time-scale photovoltaic cell power generation capacity prediction method based on hybrid model | |
| CN103606931A (en) | AVC (Automatic Voltage Control) system constant value regulation method based on probability statistics characteristic | |
| CN112052985B (en) | Middle-short-term low-voltage prediction method based on lightgbm | |
| CN117172097A (en) | Power distribution network dispatching operation method based on cloud edge cooperation and multi-agent deep learning | |
| CN105207255B (en) | A kind of power system peak regulation computational methods suitable for wind power output | |
| CN116029491B (en) | Power dispatching management system and control method thereof | |
| CN114552653B (en) | Grid-connected inverter system control parameter multi-objective optimization method based on intelligent algorithm | |
| CN115425650A (en) | Power supply station microgrid configuration method, device, equipment and medium | |
| CN112701700B (en) | Multi-objective optimization-based three-phase imbalance treatment method and system for transformer area | |
| CN114595543A (en) | Power distribution network topology identification method based on voltage correlation coefficient | |
| Zhang et al. | Optimal scheduling strategy for active distribution network considering demand side response | |
| CN108736483B (en) | Method and system for determining power grid operation mode | |
| CN111191948A (en) | Power grid operation state evaluation method and system for multi-provincial union intra-day connecting line | |
| Li et al. | K-Means Clustering Algorithm and LSTM based Short-term Load Forecasting for Distribution Transformer |
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 |