CN117556969A - Flexible power distribution network distributed reactive power optimization method based on probability scene driving - Google Patents
Flexible power distribution network distributed reactive power optimization method based on probability scene driving Download PDFInfo
- Publication number
- CN117556969A CN117556969A CN202410042537.8A CN202410042537A CN117556969A CN 117556969 A CN117556969 A CN 117556969A CN 202410042537 A CN202410042537 A CN 202410042537A CN 117556969 A CN117556969 A CN 117556969A
- Authority
- CN
- China
- Prior art keywords
- reactive power
- distribution network
- node
- scene
- 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.)
- Granted
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 92
- 238000000034 method Methods 0.000 title claims abstract description 85
- 230000001133 acceleration Effects 0.000 claims abstract description 26
- 238000010168 coupling process Methods 0.000 claims description 19
- 238000005859 coupling reaction Methods 0.000 claims description 19
- 230000008878 coupling Effects 0.000 claims description 16
- 230000009977 dual effect Effects 0.000 claims description 15
- 239000003990 capacitor Substances 0.000 claims description 11
- 238000004146 energy storage Methods 0.000 claims description 11
- 230000001105 regulatory effect Effects 0.000 claims description 10
- 230000003068 static effect Effects 0.000 claims description 9
- 238000011478 gradient descent method Methods 0.000 claims description 6
- 238000003064 k means clustering Methods 0.000 claims description 6
- 238000007599 discharging Methods 0.000 claims description 4
- 230000003190 augmentative effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 230000033228 biological regulation Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 230000007423 decrease Effects 0.000 claims description 2
- 230000006872 improvement Effects 0.000 claims description 2
- 238000012804 iterative process Methods 0.000 claims description 2
- 230000010355 oscillation Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
- Operations Research (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a distributed reactive power optimization method of a flexible power distribution network based on probabilistic scene driving, which comprises the steps of firstly, establishing a reactive power optimization model of the flexible power distribution network with the aim of minimizing system loss; secondly, comprehensively considering confidence constraints of 1-norm and ++norm, and constructing a flexible power distribution network distribution robust reactive power optimization model based on a probability scene fuzzy set; on the basis, a distributed optimization model is taken as an external framework, a consistent acceleration gradient Alternating Direction Multiplier Method (ADMM) is adopted to carry out overall coordination and updating iterative solution, and each sub-region distribution robust optimization model is taken as an internal framework, and a Column and Constraint Generation (CCG) algorithm is adopted to carry out solution; and finally, through system example simulation, the distributed reactive power optimization method for the flexible power distribution network has better convergence, and balance between economy and robustness is considered.
Description
Technical Field
The invention belongs to the technical field of reactive power optimization of a flexible power distribution network, and particularly relates to a distributed reactive power optimization method of the flexible power distribution network based on probabilistic scene driving.
Background
The intelligent soft switch is used as a novel power electronic device to replace a traditional interconnection switch to be widely interconnected in a traditional power distribution network, so that the intelligent soft switch gradually evolves into a flexible power distribution network. Compared with the traditional power distribution network, the flexible power distribution network can flexibly adjust the power flow distribution among the feed lines, improves the power supply reliability and reduces the operation loss. Meanwhile, as the distributed power supply (DG) is largely integrated into the flexible power distribution network, uncertainty of DG output and load demand makes reactive power optimization of the flexible power distribution network face a great challenge, so that how to coordinate the reactive power compensation device and SOP to improve reactive power distribution and power supply quality has very important significance.
According to the different uncertainty modeling forms of DG output and load, the uncertainty reactive power optimization method of the flexible power distribution network is mainly divided into a random optimization method and a robust optimization method. Random optimization models uncertainty of DG output force and load by adopting a discrete scene; and the robust optimization models the uncertainty of the DG output and the load by adopting a set, and the probability distribution of the DG output and the load does not need to be acquired, so that the method is easier to realize. The probability scene driving-based distributed robust optimization method combines a random optimization method and robust optimization, and searches for a solution under the most-averagely probability distribution under the fuzzy set of the probability scene, so that the accuracy of the solution can be ensured, and the robustness is good.
With the continuous development of the scale and the structure of the flexible power distribution network, the adoption of centralized optimization can lead to the problems of rapid expansion of the model solving scale, slow solving speed and the like. Thus, more and more students began to combine the alternating direction multiplier method and the DRO method to achieve distributed optimization.
Disclosure of Invention
The invention aims to solve the technical problems that: the distributed reactive power optimization method for the flexible power distribution network based on probability scene driving is used for reducing uncertainty introduced in reactive power optimization by accessing massive distributed resources into the flexible power distribution network in a layered distribution mode.
The technical scheme adopted by the invention for solving the technical problems is as follows: a distributed reactive power optimization method of a flexible power distribution network based on probability scene driving comprises the following steps:
s1: initializing, namely inputting active power distribution network line parameters, load prediction reference values, distributed photovoltaic parameters, energy storage system parameters, grouping switching capacitor parameters, static reactive power compensator parameters, on-load voltage regulating transformer parameters, intelligent soft switching parameters, ADMM algorithm parameters and CCG algorithm parameters;
s2: establishing a reactive power optimization model of the flexible power distribution network, wherein the reactive power optimization model comprises an objective function and operation constraint;
s3: clustering historical data of distributed photovoltaic output and load to obtain a typical scene set, comprehensively considering 1-norm and ++norm confidence constraint of probability distribution of each scene, and constructing a flexible distribution network distribution robust reactive power optimization model based on probability scene driving;
s4: introducing a gradient descent algorithm and a consistency theory improvement ADMM algorithm, and providing a distributed solving method of a flexible power distribution network distributed robust reactive power optimization model based on a consistency acceleration gradient ADMM algorithm;
s5: taking a flexible power distribution network distribution robust reactive power optimization model as an external frame, and adopting a consistency accelerating gradient ADMM algorithm to carry out global coordination and updating iterative solution; and taking each sub-region distribution robust model as an internal framework, solving by adopting a CCG algorithm, and outputting a result.
According to the above scheme, in the step S2, the specific steps are as follows:
s21: is provided withIs thattTime branch->The magnitude of the current flowing upwards, +.>For branch->Resistance of->Is thattTime SOP at nodeiThe active power loss at the point(s),Tfor optimizing the total number of time periods of the scheduled operation +.>And->Respectively a load node set and an SOP node set; establishing a reactive power optimization objective function of the flexible power distribution network by taking minimum network loss and SOP loss as targetsfThe method comprises the following steps:
;
s22: reactive power optimization operation constraint of the flexible power distribution network comprises distributed photovoltaic output constraint, system power flow constraint, energy storage system constraint, on-load voltage regulating transformer constraint, grouping switching capacitor constraint, static reactive compensator constraint and intelligent soft switch constraint.
Further, in the step S22, it is set thatAnd->Respectively istTime nodejActive power and reactive power emitted by the PV, < >>Is thattTime nodejPredicted power of PV, +.>Is a nodejThe capacity of the PV; the distributed photovoltaic output constraint is:
,
;
is provided withAnd->Respectively istTime inflow nodejActive power and reactive power, +.>And->Respectively istTime slave nodeiInflow branch->Active power and reactive power, +.>For branch->Reactance of->And->Respectively istTime nodejActive power and reactive power required for the load, +.>Is thattTime nodeiSquare of voltage magnitude +.>Andrespectively->Upper and lower limits of>Is thattTime branch->Square of the magnitude of the current flowing upwards, +.>;/>Is->Upper limit of->And->Respectively istTime nodejActive power when charging ESS and ESS to nodejActive power during discharge, +.>Is thattTime CB to nodejReactive power of compensation->Is thattTime SVC nodejThe reactive power of the compensation is calculated,and->Respectively istTime-of-day SOP inflow nodeiActive power and reactive power of (a); the system power flow constraint is:
,
,
,
;
set 0-1 variableAnd->Respectively istTime nodejIn the charge-discharge state of the accessed ESS, +.>And->Respectively nodesjUpper and lower limits of ESS charging power at access, < ->And->Respectively nodesjUpper and lower limits of ESS discharge power at access, < ->Is thattTime nodejESS state of charge of access>And->Respectively nodesjCharging and discharging efficiency of ESS accessed at +.>For scheduling time intervals, +.>Is a nodejESS capacity size of access>And->Respectively nodesjThe upper and lower limits of the charge state of the ESS are accessed; the energy storage system constraints are:
,
,
,
;
by introducing auxiliary nodesoDividing on-load voltage regulating transformer branch into branchesAnd branch->Two parts are provided with,/>Is thattTime auxiliary nodeoSquare of the voltage amplitude>For branch->The number of gear positions of the upper transformer,for branch->Increment of each gear of upper transformer, +.>For branch->Upper transformer transformation ratio lower limit, +.>Is thattTime branch->The gear position of the upper transformer is changed,Mconstant introduced for Big-M method, < ->、/>And->Intermediate auxiliary variables introduced for linearization; the on-load tap changer constraint is:
,
,
;
is provided withIs thattTime nodejNumber of compensation groups at CB>Is a nodejSingle set of compensation powers at CB, +.>And->Respectively istTime nodejThe number of compensation groups where CB increases and decreases and are all non-negative, and +.>For a scheduled time periodTUpper limit of the number of actions of the inner CB>The maximum compensation group number of CB; the group switching capacitor constraint is:
,
;
is provided withAnd->Respectively nodesjAt the upper and lower limits of the SVC compensation power, the static var compensator is constrained as follows:
;
is provided withIs a nodeiLoss coefficient of SOP at->And->Is a nodeiAt the upper and lower limit of reactive power for SOP transmission, < ->Is a nodeiThe capacity of the SOP; the intelligent soft switch operation constraint is:
,
,
。
according to the above scheme, in the step S3, the specific steps are as follows:
s31: through K-means clustering pairsKScene clustering is carried out on the solar distributed photovoltaic output and load historical data to obtainA limited discrete typical scene and its initial scene probability distribution +.>The method comprises the steps of carrying out a first treatment on the surface of the Is provided with->For the actual typical scene probability distribution, +.>And->Maximum deviation value of probability of typical scene under 1-norm and + -norm limits, respectively, +.>Is a typical scenesIs a function of the probability of (1),for scenes obtained by K-means clusteringsIs>And->Confidence levels of 1-norm and + -norm respectively, 0-1 auxiliary variable->And->Respectively indicate->For->A positive offset state and a negative offset state of (a); building a distributed model based on 1-norms and ≡norms probability scene ambiguity set of photovoltaic and load +.>The method comprises the following steps:
,
;
s32: constructing a flexible power distribution network distribution robust reactive power optimization model based on probability scene driving;
according to the time scale of equipment regulation, setting the gear of the on-load voltage regulating transformer, the switching group number of the grouping switching capacitors and the charge and discharge state related discrete variables of the energy storage system as first-stage decision variablesOptimizing first, setting the rest continuous variable as the decision variable of the second stage +.>The following is shown:
,
;
is provided with、/>、/>、/>、/>、/>、/>、/>、/>、/>、/>、/>、/>And->Respectively corresponding coefficient matrixes in constraint conditions, < ->Constraint set for decision variables of second stage, +.>Constraint set for decision variables of second stage, +.>Is a scenesThe magnitude of the lower distributed photovoltaic output and load; converting the deterministic reactive power optimization model of the flexible power distribution network established in the step S1 into a flexible power distribution network two-stage distribution robust reactive power optimization model based on probabilistic scene driving:
,
,
。
further, in the step S32, the two-stage robust problem is decomposed into a main problem and a sub-problem based on the CCG algorithm, and the distributed robust reactive power optimization model is solved through iteration.
According to the above scheme, in the step S4, the specific steps are as follows:
s41: is provided withIs a regionaOptimized variable of->Coupling branch for region a->Related variable of->Is the area ofaAdjacent areasbCoupling branch->Related variable of->Penalty parameter for ADMM algorithm, +.>And->Respectively the firstkRegion at time of iterationaSum regionbMiddle coupling branch->Vector composed of global variables corresponding to related variables, +.>And->Respectively the firstkRegion at time of iterationaSum regionbMiddle coupling branch->A vector of dual variables; the traditional ADMM algorithm is improved by introducing the consistency theory, and the method is obtained:
;
s42: is provided withFor the acceleration step of the gradient descent method, +.>;/>And->The dual variables and the global variables after acceleration are respectively; nesterov-based gradient descent method for ∈K>And->And (5) performing acceleration update:
。
further, in the step S4, it is set thatFor the convergence accuracy based on the consistency acceleration gradient ADMM algorithm, the original residual error of the ADMM algorithm is used +.>And dual residual->The convergence criterion for the reference construction algorithm is as follows:
;
when (when)When in use, let->The method comprises the steps of carrying out a first treatment on the surface of the Acceleration step size by changing value of convergence accuracy in iterative process>Restarting to prevent the cause->Excessive large results in excessive oscillations of the objective function and convergence residual during algorithm convergence.
Further, in the step S5, the specific steps are as follows:
s51: initializing; in the subareaaIn which the iteration number is set,/>=0,/>,/>The method comprises the steps of carrying out a first treatment on the surface of the Given->Value of (2) and convergence accuracy->;
S52: solving the sub-problems of each region; solving the distribution robustness through each subarea to obtain the value of the related variable of each area coupling branch, providing an exchange variable for the next iteration of the external consistency acceleration gradient ADMM algorithm frame, and outputting the optimal solution of the objective function of each area;
constructing an augmented Lagrangian function of the objective function of each region and a sub-regionaThe objective function of (2) is:
;
s53: exchanging each area information and updating; each region receives the exchange variable of the adjacent region and updates the dual variable of each regionGlobal variable->And acceleration step +.>;
S54: judging whether convergence exists or not; computing the sum of the original residuals of the subareas in each sceneSum of sub-region dual residuals in each scene +.>Judging the maximum value of residual error in each scene +.>Whether or not to converge to +.>Outputting an optimal solution if convergence is achieved; otherwise, executing step S55;
s55: judgingWhether or not is greater than->If yes, restarting, and (E)>The method comprises the steps of carrying out a first treatment on the surface of the Whether or notThen it is not restarted; order thek=k+1, step S52 is performed.
A computer storage medium having stored therein a computer program executable by a computer processor for performing a probabilistic scene driven distributed reactive power optimization method for a flexible distribution network.
The beneficial effects of the invention are as follows:
1. according to the flexible power distribution network distributed reactive power optimization method based on probability scene driving, firstly, a flexible power distribution network reactive power optimization model is built with the aim of minimizing system loss, meanwhile, distributed photovoltaic output and load historical data are clustered to obtain a typical scene set, confidence constraints of 1-norm and ++norm are comprehensively considered, and a flexible power distribution network distribution robust DRO reactive power optimization model based on probability scene fuzzy set is built; then a distributed optimization framework is built according to the flexible power distribution network structure, a gradient descent algorithm and a consistency theory are introduced to improve the ADMM, and a distributed solving method of a flexible power distribution network DRO reactive power optimization model based on the consistency acceleration gradient ADMM is provided; finally, carrying out global coordination and updating iterative solution by taking a distributed optimization model as an external framework and adopting a consistent acceleration gradient Alternating Direction Multiplier Method (ADMM), solving each sub-region distribution robust DRO optimization model as an internal framework by adopting a Column and Constraint Generation (CCG) algorithm, and verifying the effectiveness of the provided flexible power distribution network distributed reactive power optimization method based on probability scene driving through system example simulation; the function of reducing uncertainty introduced to reactive power optimization by accessing the flexible power distribution network through layered distribution of mass distributed resources is realized, and the method has good convergence and balances economy and robustness.
2. According to the invention, information interaction between the subareas is carried out through the coupling branch and the SOP, so that the distributed optimization of the flexible power distribution network is realized, the data interaction burden is reduced, and the solving efficiency is accelerated.
3. The ADMM algorithm based on the consistency acceleration gradient can effectively solve the distributed robust reactive power optimization model, and has good convergence performance and stability.
4. The distributed photovoltaic and load uncertainty is considered, the flexible power distribution network distributed robust reactive power optimization model is established under the comprehensive constraint of the 1-norm and the infinity-norm, the problem that the robust optimization is too conservative is solved, and the balance of economy and robustness is realized to a certain extent.
Drawings
FIG. 1 is a schematic diagram of an OLTC branch model in accordance with an embodiment of the present invention.
Fig. 2 is a distributed reactive power optimization flow chart of the flexible power distribution network based on probabilistic scene driving in an embodiment of the invention.
Fig. 3 is a diagram of an improved 33 node test system in accordance with an embodiment of the present invention.
Fig. 4 is a graph of distributed photovoltaic output and load in a typical scenario of an embodiment of the present invention.
Fig. 5 is a flow chart of an embodiment of the present invention.
Detailed Description
The invention will be described in further detail with reference to the drawings and the detailed description.
Referring to fig. 3, taking an improved 33-node system as an example, the embodiment of the invention adopts a flexible power distribution network distributed reactive power optimization method based on probability scene driving, which comprises the following steps of:
the first step: initializing, namely inputting active power distribution network line parameters, load prediction reference values, distributed Photovoltaic (PV) parameters, energy Storage System (ESS) parameters, grouping switching Capacitor (CB) parameters, static Var Compensator (SVC) parameters, on-load voltage regulating transformer (OLTC) parameters, intelligent soft Switching (SOP) parameters, ADMM algorithm parameters and CCG algorithm parameters.
In the embodiment, 33 node power distribution network line parameters are input; the access locations and parameter settings of the devices in the system are shown in table 1. The reference voltage of the flexible power distribution network is 12.66kV, and the reference capacity is 10MVA. The amplitude range of the node voltage is set to be 0.95 pu-1.05 pu, and the upper limit of the amplitude of the branch current is set to be 0.2pu; scheduling periodsT=24; the convergence accuracy of the ADMM algorithm was set to 10 -4 The convergence accuracy of the CCG algorithm was set to 10 -6 ;;/>. And clustering the distributed photovoltaic output and load historical data by adopting a K-means clustering algorithm to obtain a distributed photovoltaic output and load curve in a typical scene shown in figure 4.
Table 1 system parameter settings
And a second step of: and establishing a reactive power optimization model of the flexible power distribution network, wherein the reactive power optimization model comprises a reactive power optimization objective function of the flexible power distribution network and reactive power optimization operation constraint of the flexible power distribution network.
(1) And (3) establishing a reactive power optimization objective function of the flexible power distribution network by taking minimum network loss and SOP loss as targets:
in the method, in the process of the invention,is thattTime branch->The magnitude of the current flowing upwards; />For branch->Resistance of (2); />Is thattTime SOP at nodeiActive loss at the site;Tthe total time period number of the operation is optimized; />And->Respectively a load node set and an SOP node set.
(2) Reactive power optimization operation constraints of the flexible power distribution network include distributed Photovoltaic (PV) output constraints, system power flow constraints, energy Storage System (ESS) constraints, on-load voltage regulating transformer (OLTC) constraints, group switched Capacitor (CB) constraints, static Var Compensator (SVC) constraints, and intelligent soft Switching (SOP) constraints.
1) Distributed Photovoltaic (PV) output constraint
In the method, in the process of the invention,and->Respectively istTime nodejActive power and reactive power emitted by the PV; />Is thattTime nodejPredicted power at PV; />Is a nodejAt the capacity of the PV.
2) System power flow constraint
In the method, in the process of the invention,and->Respectively istTime inflow nodejActive power and reactive power of (a); />And->Respectively istTime slave nodeiInflow branch->Active power and reactive power of (a); />For branch->Is a reactance of (2); />And->Respectively istTime nodejActive power and reactive power required by the load; />Is thattTime nodeiSquare of the magnitude of the voltage;and->Respectively->Upper and lower limits of (2); />Is thattTime branch->Square of the magnitude of the current flowing upwards, +.>;/>Is->Upper limit of (2); />And->Respectively istTime nodejActive power when charging ESS and ESS to nodejActive power during discharge; />Is thattTime CB to nodejCompensating reactive power; />Is thattTime SVC nodejCompensating reactive power; />And->Respectively istTime-of-day SOP inflow nodeiActive power and reactive power of (a).
3) Energy Storage System (ESS) constraints
Wherein the variables 0-1And->Respectively istTime nodejThe ESS is in a charging and discharging state; />And->Respectively nodesjUpper and lower limits of the charging power of the ESS accessed; />And->Respectively nodesjUpper and lower limits of the discharge power of the ESS accessed. />Is thattTime nodejThe state of charge of the ESS accessed; />And->Respectively nodesjCharging and discharging efficiency of the ESS accessed at the position; />For a scheduling time interval; />Is a nodejThe capacity of the ESS accessed at the location; />And->Respectively nodesjAnd the upper and lower limits of the charge state of the ESS are accessed.
4) On-load step-down transformer (OLTC) constraints
The schematic diagram of the branch model containing the OLTC is shown in figure 1, and the auxiliary node is introducedoDividing an OLTC branch into branchesAnd branch->Two parts, therefore OLTC operating constraints are:
in the method, in the process of the invention,;/>is thattTime auxiliary nodeoSquare of the magnitude of the voltage; />For branch->The number of gear positions of the upper transformer; />For branch->Incremental of each gear of the upper transformer; />For branch->The upper transformer ratio is lower. />Is thattTime branch->The gear position of the upper transformer is changed,Ma larger number; />、/>And->Intermediate auxiliary variables introduced for linearization.
5) Group switched Capacitor (CB) constraints
In the method, in the process of the invention,is thattTime nodejThe number of compensation groups of the CB; />Is a nodejSingle set of compensation powers at CB;and->Respectively istTime nodejThe number of compensation groups increased and decreased by CB is non-negative; />For a scheduled time periodTAn upper limit of the number of times of the inner CB actions; />The maximum compensation group number for CB.
6) Static Var Compensator (SVC) constraint
In the method, in the process of the invention,and->Respectively nodesjAt the upper and lower limits of the SVC compensation power.
7) Intelligent soft Switch (SOP) operational constraints
/>
In the method, in the process of the invention,is a nodeiA loss coefficient of the SOP; />And->Is a nodeiThe upper and lower limits of reactive power of SOP transmission are set; />Is a nodeiThe size of the SOP.
And a third step of: clustering the distributed photovoltaic output and load historical data to obtain a typical scene set, comprehensively considering the 1-norm and the infinity-norm confidence constraint of probability distribution of each scene, and constructing a flexible distribution network Distribution Robust (DRO) reactive power optimization model based on probability scene driving.
1) Probability scene ambiguity set for distributed photovoltaic output and load
Because the randomness of the distributed photovoltaic output and the load historical data is larger, it is difficult to obtain accurate scene probability distribution through a large amount of historical data, and the obtained scene probability distribution has a certain limitation; thus through K-means cluster pairsKScene clustering is carried out on the sky-distributed photovoltaic output and load history data curve to obtainA limited discrete typical scene and its initial scene probability distribution +.>. However, the actual scene probability distribution has an error with the initial scene probability obtained by the scene clustering. In contrast, the invention constructs the probability scene fuzzy set of distributed photovoltaic and load based on 1-norm and ++norm to restrict the error range of each typical scene probability to obtainProbability scene ambiguity set +.>:
Wherein, through K-means cluster pairsKScene clustering is carried out on the solar distributed photovoltaic output and load historical data to obtainA limited discrete typical scene and its initial scene probability distribution +.>;/>For the actual typical scene probability distribution, +.>And->Maximum deviation value of probability of typical scene under 1-norm and + -norm limits, respectively, +.>Is a typical scenesProbability of->For scenes obtained by K-means clusteringsIs>And->Confidence of 1-norm and + -norm respectivelyDegree, 0-1 auxiliary variable +.>And->Respectively indicate->For->A positive offset state and a negative offset state of (a); meanwhile, the fluctuation range of the probability of the typical scene is limited by adopting the 1-norm and the +_norm, so that the conservation of the probability distribution of the worst scene can be reduced by the distribution robustness searching.
2) Distributed robust reactive power optimization model solving method based on CCG algorithm
Setting the gear of the OLTC, the switching group number of the CB and the charge and discharge state related discrete variables of the ESS as first-stage decision variables according to the time scale of equipment adjustmentOptimizing first, setting the rest continuous variable as the decision variable of the second stage +.>The following is shown:
converting the established deterministic reactive power optimization model of the flexible power distribution network into a probability scene-driven flexible power distribution network two-stage distribution robust reactive power optimization model:
in the method, in the process of the invention,,/>,/>,/>,/>,/>,/>,/>,/>,/>,/>,/>,/>,/>the coefficient matrix is the corresponding coefficient matrix in the constraint condition; />Constraint sets for second-stage decision variables; />Constraint sets for second-stage decision variables; />Is a scenesLower distributed photovoltaic output and load size.
The flexible power distribution network distribution robust reactive power optimization model based on probability scene driving is a two-stage robust optimization problem, and iterative solution can be carried out through a CCG algorithm. The CCG algorithm decomposes the two-stage robust problem into a Main Problem (MP) and a sub-problem (SP) for iterative solution, when the optimal solution difference of the MP and the SP is smaller than a given convergence accuracyAnd then, the algorithm converges and an optimal solution is output.
Fourth step: and constructing a distributed optimization framework according to the flexible power distribution network structure, introducing a gradient descent algorithm and a consistency theory to improve the ADMM, and providing a DRO reactive power optimization model distributed solving method of the flexible power distribution network based on the consistency acceleration gradient ADMM.
Is provided withIs a regionaOptimized variable of->Coupling branch for region a->Related variable of->Is the area ofaAdjacent areasbCoupling branch->Related variable of->Penalty parameter for ADMM algorithm, +.>And->Respectively the firstkRegion at time of iterationaSum regionbMiddle coupling branch->A vector composed of global variables corresponding to the related variables; />And->Respectively the firstkRegion at time of iterationaSum regionbMiddle coupling branch->A vector of dual variables; the traditional ADMM algorithm is improved by introducing the consistency theory, and the method is obtained:
;
the ADMM algorithm based on the consistency theory realizes the synchronous ADMM algorithm by solving the problems of all the subareas in parallel and then carrying out coordinated updating through coupling branches among the areas to finally obtain the global solution of the original problem. In order to further improve the convergence of the ADMM algorithm, the gradient descent method pair based on NesterovAnd->And (5) performing acceleration update:
in the method, in the process of the invention,for the acceleration step of the gradient descent method, +.>;/>And->The dual variables and the global variables after acceleration are respectively.
The invention uses the original residual error of ADMM algorithmAnd dual residual->For reference, an algorithm convergence criterion is constructed as follows:
in the method, in the process of the invention,is the convergence accuracy of the gradient ADMM algorithm based on consistency acceleration.
In order to improve the stability of ADMM algorithm in the convergence process, the invention accelerates the step length by the change value of convergence precision in the iteration processRestarting to prevent the cause->Excessive size results in excessive target function and convergence residual during algorithm convergenceOscillation, i.e. when->When in use, let->。
Fifth step: and carrying out global coordination and updating iterative solution by taking the distributed optimization model as an external framework and adopting a consistency accelerating gradient ADMM, and carrying out CCG algorithm solution by taking each subregion DRO model as an internal framework to output a result.
The distributed reactive power optimization of the flexible power distribution network based on probability scene driving is established by taking the distributed optimization of a consistency accelerating gradient ADMM algorithm as an external framework, coordinating all subareas to update global optimal solutions, taking the distributed robust optimization of all areas as an internal framework, and providing the value of a coupling branch related variable for the external distributed optimization. Fig. 2 is a distributed reactive power optimization flow chart of a flexible power distribution network based on probabilistic scene driving, which comprises the following specific steps:
(1) initializing. Setting the iteration times;/>=0,/>;/>The method comprises the steps of carrying out a first treatment on the surface of the Given->Value of (2) and convergence accuracy。
(2) Solving the sub-problems of each region. Constructing an augmented Lagrangian function of the objective function of each region to obtain the regionaFor example, the sub-region objective function is:
,
in the method, in the process of the invention,iandjrepresenting a regionaAnd coupling leg nodes of adjacent regions. And solving the distribution robustness through each subarea to obtain the value of the related variable of each area coupling branch, providing the exchange variable for the next iteration of the external consistency accelerating gradient ADMM algorithm frame, and outputting the optimal solution of the objective function of each area.
(3) Exchanging the area information and updates. Each region receives the exchange variable of the adjacent region and updates the dual variable of each regionGlobal variable->And acceleration step +.>。
(4) And judging whether convergence is achieved. Computing the sum of the original residuals of the subareas in each sceneSum of sub-region dual residuals in each scene +.>Judging the maximum value of residual error in each scene +.>Whether or not to converge to +.>If the convergence is carried out, outputting an optimal solution; otherwise, continuing to step (5).
(5) JudgingWhether or not is greater than->If yes, restarting, and (E)>The method comprises the steps of carrying out a first treatment on the surface of the Otherwise, not restarting. Order thek=k+1 returns to step (2).
Table 2 gives the optimization results for the 3 ADMM algorithms at different initial penalty parameters. It can be seen from table 2 that the algorithm 3 proposed by the present invention converges most rapidly, least time-consuming and better under different initial penalty parameters.
Table 2 results data of each algorithm
The above embodiments are merely for illustrating the design concept and features of the present invention, and are intended to enable those skilled in the art to understand the content of the present invention and implement the same, the scope of the present invention is not limited to the above embodiments. Therefore, all equivalent changes or modifications according to the principles and design ideas of the present invention are within the scope of the present invention.
Claims (9)
1. A distributed reactive power optimization method of a flexible power distribution network based on probability scene driving is characterized by comprising the following steps of: the method comprises the following steps:
s1: initializing, namely inputting active power distribution network line parameters, load prediction reference values, distributed photovoltaic parameters, energy storage system parameters, grouping switching capacitor parameters, static reactive power compensator parameters, on-load voltage regulating transformer parameters, intelligent soft switching parameters, ADMM algorithm parameters and CCG algorithm parameters;
s2: establishing a reactive power optimization model of the flexible power distribution network, wherein the reactive power optimization model comprises an objective function and operation constraint;
s3: clustering historical data of distributed photovoltaic output and load to obtain a typical scene set, comprehensively considering 1-norm and ++norm confidence constraint of probability distribution of each scene, and constructing a flexible distribution network distribution robust reactive power optimization model based on probability scene driving;
s4: introducing a gradient descent algorithm and a consistency theory improvement ADMM algorithm, and providing a distributed solving method of a flexible power distribution network distributed robust reactive power optimization model based on a consistency acceleration gradient ADMM algorithm;
s5: taking a flexible power distribution network distribution robust reactive power optimization model as an external frame, and adopting a consistency accelerating gradient ADMM algorithm to carry out global coordination and updating iterative solution; and taking each sub-region distribution robust model as an internal framework, solving by adopting a CCG algorithm, and outputting a result.
2. The flexible power distribution network distributed reactive power optimization method based on probabilistic scene driving as claimed in claim 1, wherein the method comprises the following steps: in the step S2, the specific steps are as follows:
s21: is provided withIs thattTime branch->The magnitude of the current flowing upwards, +.>For branch->Resistance of->Is thattTime SOP at nodeiThe active power loss at the point(s),Tfor optimizing the total number of time periods of the scheduled operation +.>And->Respectively a load node set and an SOP node set; establishing a reactive power optimization objective function of the flexible power distribution network by taking minimum network loss and SOP loss as targetsfThe method comprises the following steps:
;
s22: reactive power optimization operation constraint of the flexible power distribution network comprises distributed photovoltaic output constraint, system power flow constraint, energy storage system constraint, on-load voltage regulating transformer constraint, grouping switching capacitor constraint, static reactive compensator constraint and intelligent soft switch constraint.
3. The flexible power distribution network distributed reactive power optimization method based on probabilistic scene driving as claimed in claim 2, wherein the method comprises the following steps: in the step S22, set upAnd->Respectively istTime nodejActive power and reactive power emitted by the PV, < >>Is thattTime nodejPredicted power of PV, +.>Is a nodejThe capacity of the PV; the distributed photovoltaic output constraint is:
,
;
is provided withAnd->Respectively istTime inflow nodejIs provided withActive and reactive power, +.>And->Respectively istTime slave nodeiInflow branch->Active power and reactive power, +.>For branch->Reactance of->And->Respectively istTime nodejActive power and reactive power required for the load, +.>Is thattTime nodeiSquare of voltage magnitude +.>And->Respectively->Upper and lower limits of>Is thattTime branch->Square of the magnitude of the current flowing upwards, +.>;/>Is->Is set at the upper limit of (c),and->Respectively istTime nodejActive power when charging ESS and ESS to nodejThe active power at the time of the discharge,is thattTime CB to nodejReactive power of compensation->Is thattTime SVC nodejReactive power of compensation->Andrespectively istTime-of-day SOP inflow nodeiActive power and reactive power of (a); the system power flow constraint is:
,
,
,
;
set 0-1 variableAnd->Respectively istTime nodejIn the charge-discharge state of the accessed ESS, +.>And->Respectively nodesjUpper and lower limits of ESS charging power at access, < ->And->Respectively nodesjUpper and lower limits of ESS discharge power at access, < ->Is thattTime nodejESS state of charge of access>And->Respectively nodesjCharging and discharging efficiency of ESS accessed at +.>For scheduling time intervals, +.>Is a nodejESS capacity size of access>And->Respectively nodesjThe upper and lower limits of the charge state of the ESS are accessed; the energy storage system constraints are:
,
,
,
;
by introducing auxiliary nodesoDividing on-load voltage regulating transformer branch into branchesAnd branch->Two parts are provided with,/>Is thattTime auxiliary nodeoSquare of the voltage amplitude>For branch->The number of gear positions of the upper transformer,for branch->Increment of each gear of upper transformer, +.>For branch->Upper transformer transformation ratio lower limit, +.>Is thattTime branch->The gear position of the upper transformer is changed,Mconstant introduced for Big-M method, < ->、/>And->Intermediate auxiliary variables introduced for linearization; the on-load tap changer constraint is:
,
,
;
is provided withIs thattTime nodejNumber of compensation groups at CB>Is a nodejSingle set of compensation powers at CB, +.>And->Respectively istTime nodejThe number of compensation groups where CB increases and decreases and are all non-negative, and +.>For a scheduled time periodTUpper limit of the number of actions of the inner CB>The maximum compensation group number of CB; the group switching capacitor constraint is:
,
;
is provided withAnd->Respectively nodesjAt the upper and lower limits of the SVC compensation power, the static var compensator is constrained as follows:
;
is provided withIs a nodeiLoss coefficient of SOP at->And->Is a nodeiAt the upper and lower limit of reactive power for SOP transmission, < ->Is a nodeiThe capacity of the SOP; the intelligent soft switch operation constraint is:
,
,
。
4. the flexible power distribution network distributed reactive power optimization method based on probabilistic scene driving as claimed in claim 1, wherein the method comprises the following steps: in the step S3, the specific steps are as follows:
s31: through K-means clustering pairsKScene clustering is carried out on the solar distributed photovoltaic output and load historical data to obtainA limited discrete typical scene and its initial scene probability distribution +.>The method comprises the steps of carrying out a first treatment on the surface of the Is provided with->For the actual typical scene probability distribution, +.>And->Maximum deviation value of probability of typical scene under 1-norm and + -norm limits, respectively, +.>Is a typical scenesProbability of->For scenes obtained by K-means clusteringsIs>And->Confidence levels of 1-norm and + -norm respectively, 0-1 auxiliary variable->And->Respectively indicate->For->Is of positive offset form (1)A state and a negative offset state; building a distributed model based on 1-norms and ≡norms probability scene ambiguity set of photovoltaic and load +.>The method comprises the following steps:
,
;
s32: constructing a flexible power distribution network distribution robust reactive power optimization model based on probability scene driving;
according to the time scale of equipment regulation, setting the gear of the on-load voltage regulating transformer, the switching group number of the grouping switching capacitors and the charge and discharge state related discrete variables of the energy storage system as first-stage decision variablesOptimizing first, setting the rest continuous variable as the decision variable of the second stage +.>The following is shown:
,
;
is provided with、/>、/>、/>、/>、/>、/>、/>、/>、/>、/>、/>、/>And->Respectively corresponding coefficient matrixes in the constraint conditions,constraint set for decision variables of second stage, +.>Constraint set for decision variables of second stage, +.>Is a scenesThe magnitude of the lower distributed photovoltaic output and load; converting the deterministic reactive power optimization model of the flexible power distribution network established in the step S1 into a flexible power distribution network two-stage distribution robust reactive power optimization model based on probabilistic scene driving:
,
,
。
5. the flexible power distribution network distributed reactive power optimization method based on probabilistic scene driving as claimed in claim 4, wherein the method comprises the following steps: in the step S32, the two-stage robust problem is decomposed into a main problem and a sub-problem based on the CCG algorithm, and the distributed robust reactive power optimization model is solved through iteration.
6. The flexible power distribution network distributed reactive power optimization method based on probabilistic scene driving as claimed in claim 1, wherein the method comprises the following steps: in the step S4, the specific steps are as follows:
s41: is provided withIs a regionaOptimized variable of->Coupling branch for region a->Related variable of->Is the area ofaAdjacent areasbCoupling branch->Related variable of->Penalty parameter for ADMM algorithm, +.>And->Respectively the firstkRegion at time of iterationaSum regionbMiddle coupling branch->Vector composed of global variables corresponding to related variables, +.>And->Respectively the firstkRegion at time of iterationaSum regionbMiddle coupling branch->A vector of dual variables; the traditional ADMM algorithm is improved by introducing the consistency theory, and the method is obtained:
;
s42: is provided withFor the acceleration step of the gradient descent method, +.>;/>And->The dual variables and the global variables after acceleration are respectively; nesterov-based gradient descent method for ∈K>And->And (5) performing acceleration update:
。
7. the probability scene driving-based distributed reactive power optimization method for the flexible power distribution network, as claimed in claim 6, is characterized by comprising the following steps: in the step S4, set upFor the convergence accuracy based on the consistency acceleration gradient ADMM algorithm, the original residual error of the ADMM algorithm is used +.>And dual residual->The convergence criterion for the reference construction algorithm is as follows:
;
when (when)When in use, let->The method comprises the steps of carrying out a first treatment on the surface of the Acceleration step size by changing value of convergence accuracy in iterative process>Restarting to prevent the cause->Excessive large results in excessive oscillations of the objective function and convergence residual during algorithm convergence.
8. The probability scene driving-based distributed reactive power optimization method for the flexible power distribution network, as claimed in claim 7, is characterized by comprising the following steps: in the step S5, the specific steps are as follows:
s51: initializing; in the subareaaIn which the iteration number is set,/>=0,/>,/>The method comprises the steps of carrying out a first treatment on the surface of the Given->Value of (2) and convergence accuracy->;
S52: solving the sub-problems of each region; solving the distribution robustness through each subarea to obtain the value of the related variable of each area coupling branch, providing an exchange variable for the next iteration of the external consistency acceleration gradient ADMM algorithm frame, and outputting the optimal solution of the objective function of each area;
constructing an augmented Lagrangian function of the objective function of each region and a sub-regionaThe objective function of (2) is:
;
s53: exchanging each area information and updating; each region receives the exchange variable of the adjacent region and updates the dual variable of each regionGlobal variable->And acceleration step +.>;
S54: judging whether convergence exists or not; computing the sum of the original residuals of the subareas in each sceneSum of sub-region dual residuals in each scene +.>Judging the maximum value of residual error in each scene +.>Whether or not to converge to +.>Outputting an optimal solution if convergence is achieved; otherwise, executing step S55;
s55: judgingWhether or not is greater than->If yes, restarting, and (E)>The method comprises the steps of carrying out a first treatment on the surface of the Otherwise, not restarting; order thek=k+1, step S52 is performed.
9. A computer storage medium, characterized by: a computer program executable by a computer processor is stored therein, the computer program executing a probability scene driven based distributed reactive power optimization method of a flexible power distribution network as claimed in any one of claims 1 to 8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410042537.8A CN117556969B (en) | 2024-01-11 | 2024-01-11 | Flexible power distribution network distributed reactive power optimization method based on probability scene driving |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410042537.8A CN117556969B (en) | 2024-01-11 | 2024-01-11 | Flexible power distribution network distributed reactive power optimization method based on probability scene driving |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117556969A true CN117556969A (en) | 2024-02-13 |
CN117556969B CN117556969B (en) | 2024-03-26 |
Family
ID=89820844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410042537.8A Active CN117556969B (en) | 2024-01-11 | 2024-01-11 | Flexible power distribution network distributed reactive power optimization method based on probability scene driving |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117556969B (en) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3261209A1 (en) * | 2016-06-24 | 2017-12-27 | Siemens Aktiengesellschaft | Method, control apparatus, system and computer program product for reactive power and voltage control in distribution grids |
US20180375336A1 (en) * | 2015-12-22 | 2018-12-27 | Sinewatts, Inc. | Systems and methods for advanced grid integration of distributed generators and energy resources |
CN109274134A (en) * | 2018-11-08 | 2019-01-25 | 东南大学 | A kind of active distribution network robust active reactive coordination optimizing method based on time series scene analysis |
CN110032828A (en) * | 2019-05-10 | 2019-07-19 | 四川大学 | It is a kind of meter and demand response soft readjustment power distribution network two stages distribution robust D G distribute linear method rationally |
CN110460036A (en) * | 2019-05-10 | 2019-11-15 | 四川大学 | A kind of probabilistic alternating current-direct current power distribution network distributed optimization method of consideration wind-powered electricity generation |
CN113708421A (en) * | 2021-08-23 | 2021-11-26 | 国网吉林省电力有限公司长春供电公司 | Improved two-stage robust operation optimization method and system for flexible power distribution network |
CN113937825A (en) * | 2020-06-29 | 2022-01-14 | 华北电力大学(保定) | DG double-layer optimization configuration method based on E-C-Kmeans clustering and SOP optimization |
CN115733166A (en) * | 2022-11-18 | 2023-03-03 | 国网天津市电力公司经济技术研究院 | Distributed robust operation method and system for flexible power distribution network |
CN115841187A (en) * | 2022-12-28 | 2023-03-24 | 国网河北省电力有限公司经济技术研究院 | Method, device, equipment and storage medium for optimizing operation strategy of flexible power distribution network |
CN116565830A (en) * | 2023-02-21 | 2023-08-08 | 国网浙江省电力有限公司丽水供电公司 | Soft switch distribution robust optimization method based on power scene probability driving |
CN116581752A (en) * | 2023-05-24 | 2023-08-11 | 国网湖北省电力有限公司电力科学研究院 | Flexible power distribution network two-stage power supply recovery robust optimization method considering DG different operation modes |
CN116599031A (en) * | 2023-04-06 | 2023-08-15 | 国网河北省电力有限公司衡水供电分公司 | Flexible power distribution network robust optimal scheduling method based on convex hull uncertain set |
CN116722536A (en) * | 2023-06-05 | 2023-09-08 | 上海交通大学 | Flexible power distribution network distribution robust extension planning method and system based on data driving |
CN117200244A (en) * | 2023-09-11 | 2023-12-08 | 南京理工大学 | Distributed robust reactive voltage optimization method considering spatial correlation of new energy |
-
2024
- 2024-01-11 CN CN202410042537.8A patent/CN117556969B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180375336A1 (en) * | 2015-12-22 | 2018-12-27 | Sinewatts, Inc. | Systems and methods for advanced grid integration of distributed generators and energy resources |
EP3261209A1 (en) * | 2016-06-24 | 2017-12-27 | Siemens Aktiengesellschaft | Method, control apparatus, system and computer program product for reactive power and voltage control in distribution grids |
CN109274134A (en) * | 2018-11-08 | 2019-01-25 | 东南大学 | A kind of active distribution network robust active reactive coordination optimizing method based on time series scene analysis |
CN110032828A (en) * | 2019-05-10 | 2019-07-19 | 四川大学 | It is a kind of meter and demand response soft readjustment power distribution network two stages distribution robust D G distribute linear method rationally |
CN110460036A (en) * | 2019-05-10 | 2019-11-15 | 四川大学 | A kind of probabilistic alternating current-direct current power distribution network distributed optimization method of consideration wind-powered electricity generation |
CN113937825A (en) * | 2020-06-29 | 2022-01-14 | 华北电力大学(保定) | DG double-layer optimization configuration method based on E-C-Kmeans clustering and SOP optimization |
CN113708421A (en) * | 2021-08-23 | 2021-11-26 | 国网吉林省电力有限公司长春供电公司 | Improved two-stage robust operation optimization method and system for flexible power distribution network |
CN115733166A (en) * | 2022-11-18 | 2023-03-03 | 国网天津市电力公司经济技术研究院 | Distributed robust operation method and system for flexible power distribution network |
CN115841187A (en) * | 2022-12-28 | 2023-03-24 | 国网河北省电力有限公司经济技术研究院 | Method, device, equipment and storage medium for optimizing operation strategy of flexible power distribution network |
CN116565830A (en) * | 2023-02-21 | 2023-08-08 | 国网浙江省电力有限公司丽水供电公司 | Soft switch distribution robust optimization method based on power scene probability driving |
CN116599031A (en) * | 2023-04-06 | 2023-08-15 | 国网河北省电力有限公司衡水供电分公司 | Flexible power distribution network robust optimal scheduling method based on convex hull uncertain set |
CN116581752A (en) * | 2023-05-24 | 2023-08-11 | 国网湖北省电力有限公司电力科学研究院 | Flexible power distribution network two-stage power supply recovery robust optimization method considering DG different operation modes |
CN116722536A (en) * | 2023-06-05 | 2023-09-08 | 上海交通大学 | Flexible power distribution network distribution robust extension planning method and system based on data driving |
CN117200244A (en) * | 2023-09-11 | 2023-12-08 | 南京理工大学 | Distributed robust reactive voltage optimization method considering spatial correlation of new energy |
Non-Patent Citations (4)
Title |
---|
QIULING YANG 等: "Deep Policy Gradient for Reactive Power Control in Distribution Systems", 《2020 IEEE INTERNATIONAL CONFERENCE ON COMMUNICATIONS, CONTROL, AND COMPUTING TECHNOLOGIES FOR SMART GRIDS (SMARTGRIDCOMM)>, 30 December 2020 (2020-12-30), pages 1 - 6 * |
张志荣 等: "基于ADMM 的交直流混合配电网分布式电压优化模型", 《电网技术》, vol. 45, no. 11, 30 November 2021 (2021-11-30), pages 4551 - 4559 * |
杨行: "主动配电网中分布式电源与智能软开关优化配置", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》, no. 2021, 15 February 2021 (2021-02-15), pages 042 - 1253 * |
王家怡 等: "考虑风电不确定性的交直流混合配电网分布式优化运行", 《中国电机工程学报》, vol. 40, no. 2, 20 January 2020 (2020-01-20), pages 550 - 562 * |
Also Published As
Publication number | Publication date |
---|---|
CN117556969B (en) | 2024-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105846461B (en) | Control method and system for large-scale energy storage power station self-adaptive dynamic planning | |
CN113285485B (en) | Power distribution network source network charge storage multi-terminal cooperative voltage regulation method under long, short and multi-time scales | |
CN114362196A (en) | Multi-time-scale active power distribution network voltage control method | |
CN111293682B (en) | Multi-microgrid energy management method based on cooperative model predictive control | |
CN113516278B (en) | Active power distribution network multi-time scale active and reactive power coordination optimization scheduling method and system | |
CN110768262A (en) | Active power distribution network reactive power supply configuration method based on node clustering partition | |
CN106712035A (en) | Economic dispatching method for power system | |
CN114662751B (en) | Garden multifunctional short-term load forecasting and optimizing method based on LSTM | |
CN113224769A (en) | Multi-time scale power distribution network voltage optimization method considering photovoltaic multi-state adjustment | |
CN116388302A (en) | Active-reactive power combined optimization method for power distribution network for coordinating network side resources | |
CN117039882A (en) | Resource aggregation regulation and control method and system based on binary consistency algorithm | |
CN115345380A (en) | New energy consumption electric power scheduling method based on artificial intelligence | |
CN115481856A (en) | Comprehensive energy system multi-scale scheduling method and system considering comprehensive demand response | |
CN114784831A (en) | Active power distribution network multi-objective reactive power optimization method based on mobile energy storage | |
Linlin et al. | Research on multi-objective reactive power optimization of power grid with high proportion of new energy | |
CN117556969B (en) | Flexible power distribution network distributed reactive power optimization method based on probability scene driving | |
Hao et al. | Solving economic load dispatch problem of power system based on differential evolution algorithm with different mutation strategies | |
CN113178868A (en) | Voltage optimization method for flexible interconnection power distribution system | |
CN111160808A (en) | Distributed event-triggered power system economic dispatching method with uncertain parameters | |
CN110556878A (en) | Distributed voltage control optimization method and system applied to wind power plant | |
Radziukynas et al. | Optimization methods application to optimal power flow in electric power systems | |
CN104408531B (en) | A kind of uniform dynamic programming method of multidimensional multistage complicated decision-making problems | |
CN115719113A (en) | Intelligent power grid economic dispatching distributed accelerated optimization method based on directed imbalance topology | |
CN116050576A (en) | Flexible resource coordination optimization method and system for active power distribution network | |
Nasir et al. | VSI improvement using SVC with aid of a modified lightning search algorithm |
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 |