CN109873447B - Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network - Google Patents

Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network Download PDF

Info

Publication number
CN109873447B
CN109873447B CN201910124910.3A CN201910124910A CN109873447B CN 109873447 B CN109873447 B CN 109873447B CN 201910124910 A CN201910124910 A CN 201910124910A CN 109873447 B CN109873447 B CN 109873447B
Authority
CN
China
Prior art keywords
power
active
time
scheduling
ess
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
Application number
CN201910124910.3A
Other languages
Chinese (zh)
Other versions
CN109873447A (en
Inventor
胡伟
邓振立
许晓慧
汪春
马洲俊
王勇
梁硕
夏俊荣
吴奕
荆江平
杨梓俊
许洪华
陈曦
陈逸如
李小荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Electric Power Research Institute Co Ltd CEPRI
State Grid Jiangsu Electric Power Co Ltd
Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
State Grid Jiangsu Electric Power Co Ltd
Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Electric Power Research Institute Co Ltd CEPRI, State Grid Jiangsu Electric Power Co Ltd, Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201910124910.3A priority Critical patent/CN109873447B/en
Publication of CN109873447A publication Critical patent/CN109873447A/en
Application granted granted Critical
Publication of CN109873447B publication Critical patent/CN109873447B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a multi-source cooperative active-reactive power regulation and control method for multiple time levels of an active power distribution network, which comprises the following steps: (1) Predicting wind power, photovoltaic and load power sequences; (2) Carrying out uncertain sampling of source-load prediction errors to obtain a scene set covering the whole sampling space; (3) Scene reduction is carried out by utilizing a synchronous back-substitution reduction technology to obtain a typical scene set; (4) In a long-time scale scheduling stage, obtaining an adjusting scheme of adjusting equipment with adjusting frequency limitation and an energy storage charging and discharging state; (5) In a short time scale scheduling stage, a scheduling instruction of the rapid continuous adjustable equipment is obtained through solving; (6) And in the real-time control stage, each adjustable device is sequentially adjusted based on the sensitivity analysis of the node voltage, so that the safe and stable operation of the power grid is guaranteed. The method of the invention can realize the coordinated optimization operation of each adjustable device in the active power distribution network, and the full consumption and the efficient utilization of renewable energy sources.

Description

Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network
Technical Field
The invention relates to a multi-source cooperative active power distribution network multi-time level active-reactive power regulation and control method, and belongs to the field of active power distribution network optimized operation.
Background
At present, in the face of the double pressure of energy exhaustion and environmental deterioration, distributed Generation (DG) having renewable green clean characteristics is rapidly moving to the stage of the world electric power industry. Meanwhile, the energy storage technology changes the problem of simultaneity of transmission and distribution of the traditional power system due to the time sequence energy regulation function of the energy storage technology, and has wide development prospect in the power system. The permeability of DG, ESS, SVC equipment in the distribution network improves day by day, and traditional distribution network is gradually evolving into the initiative distribution network that has numerous adjustable controllable resources, its core characterized in that is active and initiative: the active power distribution network contains distributed energy such as a fan (WT), a Photovoltaic (PV) and an ESS; the initiative refers to that the active power distribution network has active control capability and can perform coordination control and active management on controllable resources in the network.
The uncertainty of wind power and photovoltaic power brings great challenges to the safe operation and optimal scheduling of the power system. The DG grid connection with high permeability can cause voltage fluctuation or overvoltage to cause grid disconnection, the capacity of an active power distribution network for absorbing renewable energy sources for power generation is seriously restricted, and power grid resources and renewable energy sources are wasted. Core technologies capable of better solving the problems of overvoltage, wind abandonment, light abandonment, power blockage, line loss and the like of the active power distribution network are urgently needed. A multi-source cooperative active-reactive power regulation and control method for multiple time levels of an active power distribution network is an effective tool for solving the problems.
Disclosure of Invention
The invention provides a multi-time-level active-reactive power regulation and control method for a multi-source collaborative active power distribution network. The method adopts the principles of time sequence progressive and progressive refinement, takes the operation characteristics of adjustable resources into consideration, arranges the scheduling strategy of each device from a prospective view, finely manages and controls the problems occurring in real-time operation, can solve the key problems of overvoltage, wind abandonment, light abandonment, power blockage, line loss and the like of the active power distribution network, and ensures the safe, stable, efficient and green operation of the power grid.
The invention provides a multi-time-level active-reactive power regulation and control method for a multi-source collaborative active power distribution network, which comprises the following steps:
s1: predicting wind power, photovoltaic and load power based on a time sequence algorithm of a wavelet-BP neural network, decomposing a wind power output sequence, a photovoltaic output sequence and a load power sequence on different scales by utilizing wavelet transformation, predicting components of different frequencies by using a plurality of BP neural networks, and reconstructing each prediction result to obtain a complete prediction result; the method comprises the following specific steps:
s1.1: 3-scale wavelet decomposition is carried out on the original power sequence x to obtain a low-frequency part l 3 And a high-frequency part h 1 、h 2 And h 3
S1.2: using BP neural network to measure low frequency part l 3 And a high-frequency part h 1 、h 2 、h 3 Performing prediction to obtain a low-frequency part L of the prediction component 3 And a high frequency part H 1 、H 2 And H 3
S1.3: reconstructing the low-frequency prediction result and the high-frequency prediction result to obtain a complete prediction power sequence x * =L 3 +H 1 +H 2 +H 3
S2: based on a wind power, photovoltaic and load prediction error probability model, adopting Latin hypercube sampling to sample source-load uncertainty, and obtaining a scene set capable of covering all sampling spaces;
in consideration of the multi-dimensional randomness of wind power, photovoltaic and load which need to be described simultaneously in the generation of the source-load probability scene, the method utilizes a Latin hypercube sampling method to carry out multi-dimensional layered sampling by means of a multi-dimensional sampling theory and based on a wind power, photovoltaic and load prediction error probability model, and effectively avoids the phenomenon that the number of scenes is still too much after reduction caused by independent sampling of the wind power, the photovoltaic and the load.
The probability error model of the wind-solar load is as follows:
s2.1: the prediction error probability of the wind power output power meets Beta distribution, and the probability density function is as follows:
Figure BDA0001973206170000021
in the formula, alpha and beta are respectively a shape parameter and a scale parameter of a wind power output power prediction error probability density function; Δ P is the power prediction error.
S2.2: the load and photovoltaic contribution prediction errors are generally considered to satisfy a normal distribution, i.e., Δ P-N (μ, σ) 2 ) The probability density function is:
Figure BDA0001973206170000022
where μ and σ are the expected and standard deviation of the prediction error, respectively.
S3: scene reduction is carried out by utilizing a synchronous back-substitution reduction technology, a typical scene which can better reflect an original scene set is obtained, and corresponding probability of the typical scene is determined;
s4: in a long-time scale scheduling stage, a long-time scale active-reactive scheduling model of the active power distribution network facing the full consumption of renewable energy sources is built; obtaining a regulation scheme of regulation equipment with regulation frequency limit, such as OLTC, CB and the like, and the charge-discharge state of the ESS;
in multi-time scale scheduling, a long-time scale optimization scheduling model facing renewable energy sources to be fully consumed for the active power distribution network is built in a long-time scale scheduling stage, the action times of discrete reactive power adjusting equipment are strictly limited in a scheduling period, and in the long-time scale scheduling stage, adjusting schemes of adjusting equipment with adjusting time limit, such as OLTC (on-line communication technology), CB (circuit board) and the like are determined; considering the influence of frequent charge-discharge conversion on the ESS, determining the charge-discharge state of the ESS at the same time, and keeping the charge-discharge state unchanged in short-time scale scheduling, wherein an objective function of an active power distribution network long-time scale optimization scheduling model for full consumption of renewable energy sources is as follows:
Figure BDA0001973206170000031
in the formula (13), C T And C CB The calculation method of (c) is shown in equation (14).
Figure BDA0001973206170000032
In formula (13), P g,t ,P ess,t ,B ess,t ,Q g,t ,Q SVC,t ,Q CB,t ,Q dis,t ,k T,t Actual active power of renewable energy, charging and discharging power of ESS, charging and discharging states of ESS, reactive output power of renewable energy, reactive compensation power of SVC, discrete compensation power of CB, reactive output power of ESS and switching gear of OLTC at the moment t respectively; t is a regulation and control period; delta T is a long time scale regulation time interval; e is a system branch set; i all right angle ij,t The square of the current amplitude flowing through the branch ij at the moment t; r is a radical of hydrogen ij Resistance for branch ij; c loss 、C T And C CB The unit cost of the system operation grid loss, the gear adjusting cost in the OLTC dispatching cycle and the gear adjusting cost in the dispatching cycle of the capacitor group CB are respectively. The last term of the formula (13) is a power reduction penalty term of the renewable energy source, and lambda is an active power reduction penalty coefficient of the renewable energy source; phi is a re A node set accessed by renewable energy sources in the power distribution network;
Figure BDA0001973206170000033
and predicting the active power output of the renewable energy source at the moment t.
In the formula (14), p T And ρ C Unit adjustment costs of OLTC and CB, respectively; Δ K T And Δ K CB The adjustment times of the scheduling periods of the OLTC and the CB are respectively;
s5: on the basis of S4, in a short-time scale scheduling stage, a scheduling scheme of OLTC and CB obtained in a long-time scale scheduling stage is fixed, meanwhile, the influence of frequent charge-discharge conversion on the service life of the ESS is considered, in the short-time scale scheduling, the charge-discharge state of the ESS is consistent with the long-time scale, and on the basis of the prediction information of short-time scale wind, light and load, the minimum correction of a control variable is taken as an optimization target; solving to obtain a scheduling instruction of adjustable equipment such as wind power, photovoltaic, ESS and SVC;
the objective function of the short-time scale scheduling model for the efficient utilization of renewable energy sources is as follows:
Figure BDA0001973206170000041
wherein the content of the first and second substances,
S(k+i|k)=[P g,t (k+i|k),P ess,t (k+i|k),Q g,t (k+i|k),Q SVC,t (k+i|k),Q dis,t (k+i|k)] T (16)
Figure BDA0001973206170000042
in the formula (15), S (k + i | k) represents a scheduling plan vector of the short-time-scale-adjustable device for predicting the future k + i time at the time k; p g,t (k+i|k),P ess,t (k+i|k),Q g,t (k+i|k),Q SVC,t (k+i|k),Q dis,t (k + i | k) predicting the active output of the renewable energy source, the charging and discharging power of the ESS, the reactive output of the renewable energy source and the reactive output power scheduling sequence of the ESS at the future time k + i for the time k respectively;
Figure BDA0001973206170000043
scheduling plan vectors of all adjustable devices obtained by long-time scale scheduling are used as reference values of short-time scale scheduling;
Figure BDA0001973206170000044
Figure BDA0001973206170000045
the active output power of the renewable energy source, the charging and discharging power of the ESS, the reactive output power of the renewable energy source and the reactive output power of the ESS at the k + i moment of the long-time scale scheduling stage are scheduled in a plan sequence; and N is a short time scale scheduling step length.
S6: and in the real-time control stage, monitoring the voltage of each node, and when node voltage is out of limit, sequentially adjusting each adjustable device based on the sensitivity of the node where each adjustable device is located to the voltage out-of-limit node, so as to ensure the safe and stable operation of the power grid. The method for calculating the sensitivity of each adjustable device to the voltage out-of-limit node comprises the following steps:
s6.1: based on the alternating current power flow equation, linearizing the nonlinear power flow equation at the steady state solution, and obtaining the following matrix expression:
Figure BDA0001973206170000046
in the formula, Δ P and Δ Q are respectively an active and reactive variable matrix of each node injection system; delta theta and delta V are respectively variable matrixes of voltage phase angle and amplitude of each node; j is Jacobian matrix.
S6.2: inverting equation (18) in S6.1 yields:
Figure BDA0001973206170000051
s6.3: the relation between the voltage amplitude of each node in the power distribution network and the active power-reactive power can be expressed as follows:
ΔV=S 21 ΔP+S 22 ΔQ (20)
wherein, in the formula (20), S 21 、S 22 The voltage-active sensitivity factor and the voltage-reactive sensitivity factor are respectively, and the sensitivity of the node where the adjustable equipment is located to the voltage out-of-limit node can be measured by S 21 And S 22 And (4) calculating.
The invention has the beneficial effects that:
(1) Considering that the line resistance and the reactance value in the power distribution network are close to each other, so that the active-reactive power coupling is stronger, the method calculates the operating characteristics of each adjustable device from the perspective of active-reactive coordination optimization regulation and control, and performs comprehensive coordination optimization regulation and control on the adjustable devices in the active power distribution network;
(2) Compared with the prior art, the multi-time-level active-reactive power regulation and control method for the multi-source collaborative active power distribution network has the advantages that scheduling problems and control problems in optimized operation of the power distribution network are considered cooperatively, so that full consumption and efficient operation of renewable energy sources in the power distribution network are facilitated, and the safe, stable and efficient operation of the system is guaranteed;
(3) Compared with the prior art, the invention fully utilizes the PCS to provide reactive support for the active Power distribution network, improves the utilization rate of the ESS and exerts the value of the ESS to a greater extent;
(4) According to the multi-time-level active-reactive power regulation and control method, in a long-time scale scheduling stage, the regulation scheme of regulation equipment with regulation frequency limit, such as OLTC and CB, and the charging and discharging state of the ESS are determined, and the influence of frequent regulation of the OLTC and CB and frequent charging and discharging conversion of the ESS on the service life of the ESS is solved; meanwhile, on the basis of the prediction information of wind, light and load on the short time scale, the scheduling instruction of equipment which can be rapidly and frequently adjusted such as wind power, photovoltaic, ESS and SVC is corrected, the influence of inaccurate prediction of wind, light and load on the charging and discharging state scheduling instruction of OLTC, CB and ESS on the long time scale is made up, and the reliability of system operation is enhanced;
(5) In the real-time control stage, when the node voltage is out of limit, each adjustable device is adjusted in sequence based on the sensitivity of the node where each adjustable device is located to the voltage out-of-limit node, so that the impact of simultaneous action of multiple devices on the safety and stability of the system is reduced to a greater extent, and the safe and stable operation of a power grid can be effectively ensured.
Drawings
Fig. 1 is an overall flow chart of a multi-time-level active-reactive power regulation and control method of a multi-source collaborative active power distribution network of the invention;
FIG. 2 is a flow chart of prediction of wind power, photovoltaic and load power based on a wavelet-BP neural network time series algorithm.
Fig. 3 is a schematic diagram of a 3-scale wavelet decomposition.
Detailed Description
The present invention will be described in detail with reference to the accompanying drawings.
Fig. 1 is a general method flowchart of the present invention, and the present invention provides a multi-time-level active-reactive power regulation method for a multi-source cooperative active power distribution network, which includes the following steps:
s1: predicting wind power, photovoltaic and load power based on a time series algorithm of a wavelet-BP neural network, decomposing a wind power output sequence, a photovoltaic output sequence and a load power sequence on different scales by utilizing wavelet transformation, predicting components of different frequencies by using a plurality of BP neural networks, and reconstructing each prediction result to obtain a complete prediction result; as shown in the attached figure 2, the method comprises the following specific steps:
s1.1: 3-scale wavelet decomposition is carried out on the original power sequence x to obtain a low-frequency part l 3 And a high-frequency part h 1 、h 2 And h 3 As shown in fig. 3;
s1.2: using BP neural network to measure low frequency part l 3 And a high-frequency part h 1 、h 2 、h 3 Performing prediction to obtain a low-frequency part L of the prediction component 3 And a high frequency part H 1 、H 2 And H 3
S1.3: reconstructing the low-frequency prediction result and the high-frequency prediction result to obtain a complete prediction power sequence x * =L 3 +H 1 +H 2 +H 3
S2: based on a wind power, photovoltaic and load prediction error probability model, adopting Latin hypercube sampling to sample source-load uncertainty, and obtaining a scene set capable of covering all sampling spaces;
in consideration of the multi-dimensional randomness of wind power, photovoltaic and load which need to be described simultaneously in the generation of the source-load probability scene, the method utilizes a Latin hypercube sampling method to carry out multi-dimensional layered sampling by means of a multi-dimensional sampling theory and based on a wind power, photovoltaic and load prediction error probability model, and effectively avoids the phenomenon that the number of scenes is still too much after reduction caused by independent sampling of the wind power, the photovoltaic and the load.
The probability error model of the wind-solar load is as follows:
s2.1: the prediction error probability of the wind power output power meets Beta distribution, and the probability density function is as follows:
Figure BDA0001973206170000071
in the formula, alpha and beta are respectively a shape parameter and a scale parameter of a wind power output power prediction error probability density function; Δ P is the power prediction error.
S2.2: the load and photovoltaic contribution prediction errors are generally considered to satisfy a normal distribution, i.e., Δ P-N (μ, σ) 2 ) The probability density function is:
Figure BDA0001973206170000072
where μ and σ are the expected and standard deviation of the prediction error, respectively.
S3: scene reduction is carried out by utilizing a synchronous back-substitution reduction technology, a typical scene which can better reflect an original scene set is obtained, and corresponding probability of the typical scene is determined;
s4: in a long-time scale scheduling stage, a long-time scale active-reactive scheduling model of the active power distribution network facing the full consumption of renewable energy sources is built; obtaining an adjustment scheme of adjustment equipment with adjustment frequency limit, such as an on-load tap changing transformer, a capacitor bank and the like, and a charge-discharge state of an energy storage system;
in multi-time scale scheduling, a long-time scale optimization scheduling model facing renewable energy sources to be fully consumed for the active power distribution network is built in a long-time scale scheduling stage, the action times of discrete reactive power adjusting equipment are strictly limited in a scheduling period, and in the long-time scale scheduling stage, adjusting schemes of adjusting equipment with adjusting time limit, such as OLTC (on-line communication technology), CB (circuit board) and the like are determined; considering the influence of frequent charge-discharge conversion on the ESS, determining the charge-discharge state of the ESS, and keeping the charge-discharge state unchanged in short-time scale scheduling, wherein an objective function of a long-time scale optimization scheduling model of the active power distribution network for fully consuming renewable energy sources is as follows:
Figure BDA0001973206170000073
in the formula (23), C T And C CB The calculation method of (2) is shown in formula (24).
Figure BDA0001973206170000074
In formula (23), P g,t ,P ess,t ,B ess,t ,Q g,t ,Q SVC,t ,Q CB,t ,Q dis,t ,k T,t Actual active power of renewable energy sources, charging and discharging power of an ESS, charging and discharging states of the ESS, reactive output power of the renewable energy sources, reactive compensation power of an SVC, discrete compensation power of a CB, reactive output power of the ESS and switching gears of an OLTC at the moment t respectively; t is a regulation and control period; delta T is a long time scale regulation time interval; e is a system branch set; i all right angle ij , t The square of the current amplitude flowing through the branch ij at the moment t; r is ij Resistance for branch ij; c loss 、C T And C CB The unit cost of the system operation network loss, the gear adjusting cost in the OLTC dispatching cycle and the gear adjusting cost in the dispatching cycle of the capacitor group CB are respectively. The last item of the formula (23) is a power reduction penalty item of the renewable energy source, and lambda is an active power reduction penalty coefficient of the renewable energy source; phi is a unit of re A node set accessed by renewable energy sources in the power distribution network;
Figure BDA0001973206170000084
the predicted active power output of the renewable energy source at time t.
In the formula (24), p T And ρ C Unit adjustment costs of OLTC and CB, respectively; Δ K T And Δ K CB The adjustment times of the scheduling periods of the OLTC and the CB are respectively;
s5: on the basis of S4, in short time scale scheduling, on the basis of more accurate short-term prediction information of wind power, photovoltaic and load obtained from a time sequence algorithm of a wavelet-BP neural network, the switching capacity Q of the capacitor group CB obtained in a long time scale is kept CB,t Switching gear k of OLTC T,t And charging and discharging state B of ESS ess,t The minimum correction quantity of the other controlled variables is optimized without changingThe method comprises the following steps of aiming at building a short time scale scheduling model for efficient utilization of renewable energy, and solving to obtain scheduling instructions of adjustable equipment such as wind power, photovoltaic, ESS and SVC; the objective function of the short-time scale scheduling model for efficient utilization of renewable energy sources is as follows:
Figure BDA0001973206170000081
wherein the content of the first and second substances,
S(k+i|k)=[P g,t (k+i|k),P ess,t (k+i|k),Q g,t (k+i|k),Q SVC,t (k+i|k),Q dis,t (k+i|k)] T (26)
Figure BDA0001973206170000082
in the formula (25), S (k + i | k) represents a scheduling plan vector of the short-time-scale-adjustable device for predicting the k + i time in the future at the time k; p g,t (k+i|k),P ess,t (k+i|k),Q g,t (k+i|k),Q SVC,t (k+i|k),Q dis,t (k + i | k) predicting the active output of the renewable energy source, the charging and discharging power of the ESS, the reactive output of the renewable energy source and the reactive output power scheduling sequence of the ESS at the future time k + i for the time k respectively;
Figure BDA0001973206170000083
scheduling plan vectors of all adjustable devices obtained by long-time scale scheduling are used as reference values of short-time scale scheduling;
Figure BDA0001973206170000091
Figure BDA0001973206170000092
the active output power of the renewable energy source, the charging and discharging power of the ESS, the reactive output power of the renewable energy source and the reactive output power of the ESS at the k + i moment of the long-time scale scheduling stage are scheduled in a plan sequence; and N is a short time scale scheduling step length.
S6: and in the real-time control stage, monitoring the voltage of each node, and when node voltage is out of limit, sequentially adjusting each adjustable device based on the sensitivity of the node where each adjustable device is located to the voltage out-of-limit node, so as to ensure the safe and stable operation of the power grid. The method for calculating the sensitivity of each adjustable device to the voltage out-of-limit node comprises the following steps:
s6.1: based on the alternating current power flow equation, the nonlinear power flow equation is linearized at the steady state solution, and the following matrix expression can be obtained:
Figure BDA0001973206170000093
in the formula, Δ P and Δ Q are respectively an active and reactive variable matrix of each node injection system; delta theta and delta V are respectively variable matrixes of voltage phase angle and amplitude of each node; j is the Jacobian matrix.
S6.2: inversion of equation (28) in S6.1 yields:
Figure BDA0001973206170000094
s6.3: the relation between the voltage amplitude of each node in the power distribution network and the active power-reactive power can be expressed as follows:
ΔV=S 21 ΔP+S 22 ΔQ (30)
wherein, in the formula (30), S 21 、S 22 The voltage-active sensitivity factor and the voltage-reactive sensitivity factor are respectively, and the sensitivity of the node where the adjustable equipment is located to the voltage out-of-limit node can be measured by S 21 And S 22 And (4) calculating.
The foregoing describes specific embodiments of the present invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are intended as illustrations of the strategy of the invention, but that various changes and modifications may be made without departing from the spirit and scope of the invention, and the invention is to be construed as falling within the scope of the invention as claimed. The scope of the invention is defined by the appended claims.

Claims (5)

1. A multi-time-level active-reactive power regulation and control method for a multi-source collaborative active power distribution network comprises the following steps:
step S1: predicting wind power, photovoltaic and load power sequences based on a time sequence algorithm of a wavelet-BP neural network;
step S2: based on a wind power, photovoltaic and load prediction error probability model, performing source-load uncertainty sampling by adopting Latin hypercube sampling to obtain a scene set capable of covering all sampling spaces;
and step S3: scene reduction is carried out by utilizing a synchronous back-substitution reduction technology, a typical scene capable of reflecting an original scene set is obtained, and the corresponding probability of the typical scene is determined;
and step S4: in a long-time scale scheduling stage, building a day-ahead active-reactive power scheduling model of the active power distribution network facing the full consumption of renewable energy; obtaining a regulation scheme of discrete regulation equipment with regulation frequency limit of an on-load tap changer (OLTC) and a Capacitor Bank (CB) and a charging and discharging state of an Energy Storage System (ESS);
step S5: in a short-time scale scheduling stage, on the basis of the predicted power of short-time scale wind power, photovoltaic and load, the minimum correction quantity of a control variable is taken as an optimization target; solving to obtain a scheduling instruction of the wind power, photovoltaic and energy storage system and SVC adjustable equipment;
step S6: in the real-time control stage, monitoring the voltage of each node, and when node voltage is out of limit, sequentially adjusting each adjustable device based on the sensitivity of the node where each adjustable device is located to the voltage out-of-limit node to ensure the safe and stable operation of the power grid;
in the step S4, in the long time scale scheduling stage, the adjusting scheme of the adjusting equipment with the adjusting frequency limit and the charging and discharging state of the stored energy are determined and are kept unchanged in the short time scale scheduling; the objective function of the active power distribution network long-time scale optimization scheduling model for fully absorbing renewable energy sources is as follows:
Figure FDA0003726181130000011
in the formula (3), C T And C CB The calculation method of (2) is shown in formula (4):
Figure FDA0003726181130000012
in the formula (3), P g,t ,P ess,t ,B ess,t ,Q g,t ,Q SVC,t ,Q CB,t ,Q dis,t ,k T,t Actual active power of renewable energy sources, charging and discharging power of the ESS, charging and discharging states of the ESS, reactive output power of the renewable energy sources, reactive compensation power of a Static Var Compensator (SVC), discrete compensation power of a CB, reactive output power of the ESS and switching gears of an OLTC at the moment t respectively; t is a regulation and control period; delta T is a long time scale regulation time interval; e is a system branch set; i.e. i ij,t The square of the current amplitude flowing through the branch ij at the moment t; r is ij Resistance for branch ij; c loss 、C T And C CB Unit cost of system operation network loss, gear adjusting cost in an OLTC dispatching cycle and gear adjusting cost in a dispatching cycle of a capacitor group CB are respectively set; the last item of the formula (3) is a power reduction penalty item of the renewable energy source, and lambda is an active power reduction penalty coefficient of the renewable energy source; phi is a unit of re A node set accessed by renewable energy sources in the power distribution network;
Figure FDA0003726181130000021
predicting the active power output of the renewable energy source at the time t;
in the formula (4) (. Rho) T And ρ C Unit adjustment costs of OLTC and CB, respectively; Δ K T And Δ K CB The adjustment times of the scheduling periods of the OLTC and CB, respectively.
2. The multi-source cooperative active power distribution network multi-time level active-reactive power regulation and control method of claim 1, wherein: in step S1, the specific steps are as follows:
s1.1: 3-scale wavelet decomposition is carried out on the original power sequence x to obtain a low-frequency part l 3 And a high-frequency part h 1 、h 2 And h 3
S1.2: using BP neural network to measure low frequency part l 3 And a high-frequency part h 1 、h 2 、h 3 Performing prediction to obtain a low-frequency part L of the prediction component 3 And a high frequency part H 1 、H 2 And H 3
S1.3: reconstructing the low-frequency prediction result and the high-frequency prediction result to obtain a complete prediction power sequence x * =L 3 +H 1 +H 2 +H 3
3. The multi-source cooperative active power distribution network multi-time-level active-reactive power regulation and control method of claim 1, wherein: in the step S2, the process is carried out,
the probability error model of the wind-solar load is as follows:
s2.1: the prediction error probability of the wind power meets Beta distribution, and the probability density function is as follows:
Figure FDA0003726181130000022
in the formula, alpha and beta are respectively a shape parameter and a scale parameter of a wind power output power prediction error probability density function; Δ P is the power prediction error;
s2.2: the photovoltaic and load prediction errors satisfy a normal distribution, and the DeltaP follows the normal distribution, i.e., deltaP-N (mu, sigma) 2 ) The probability density function is:
Figure FDA0003726181130000031
where μ and σ are the expected and standard deviation of the prediction error, respectively.
4. The multi-source cooperative active power distribution network multi-time level active-reactive power regulation and control method of claim 1, wherein: in the step S5, on the basis of the predicted power of the short-time-scale wind power, the photovoltaic and the load, the minimum correction quantity of the control variable is taken as an optimization target; solving to obtain a scheduling instruction of the wind power, photovoltaic, ESS and SVC adjustable equipment; the objective function of the short-time scale scheduling model for the efficient utilization of renewable energy sources is as follows:
Figure FDA0003726181130000032
wherein the content of the first and second substances,
S(k+i|k)=[P g,t (k+i|k),P ess,t (k+i|k),Q g,t (k+i|k),Q SVC,t (k+i|k),Q dis,t (k+i|k)] T (6)
Figure FDA0003726181130000033
in the formula (5), S (k + i | k) represents a scheduling plan vector of the short-time-scale adjustable device for predicting the future k + i time at the time k; p g,t (k+i|k),P ess,t (k+i|k),Q g,t (k+i|k),Q SVC,t (k+i|k),Q dis,t (k + i | k) predicting a renewable energy active output, ESS charge and discharge power, renewable energy reactive output and ESS reactive output power scheduling sequence of the ESS at the future time k respectively;
Figure FDA0003726181130000034
scheduling plan vectors of all adjustable devices obtained by long-time scale scheduling are used as reference values of short-time scale scheduling;
Figure FDA0003726181130000035
Figure FDA0003726181130000036
are respectively longA renewable energy active power output, ESS charge and discharge power, renewable energy reactive output power and ESS reactive output power scheduling plan sequence at the time of k + i in the time scale scheduling stage; and N is a short time scale scheduling step size.
5. The multi-source cooperative active power distribution network multi-time level active-reactive power regulation and control method of claim 1, wherein:
in step S6, the method for calculating the sensitivity of the node where each adjustable device is located to the voltage out-of-limit node is as follows:
s6.1: based on the alternating current power flow equation, linearizing the nonlinear power flow equation at the steady state solution to obtain the following matrix expression:
Figure FDA0003726181130000041
in the formula, Δ P and Δ Q are respectively an active and reactive variable matrix of each node injection system; delta theta and delta V are respectively variable matrixes of voltage phase angle and amplitude of each node; j is Jacobian matrix;
s6.2: inversion of equation (8) in S6.1 yields:
Figure FDA0003726181130000042
s6.3: the relation between the voltage amplitude of each node in the power distribution network and the active power-reactive power can be expressed as follows:
ΔV=S 21 ΔP+S 22 ΔQ (10)
wherein, in formula (10), S 21 、S 22 The voltage-active and voltage-reactive sensitivity factors are respectively, and the sensitivity of the node where each adjustable device is located to the voltage out-of-limit node can be measured by S 21 And S 22 And (4) calculating.
CN201910124910.3A 2019-02-19 2019-02-19 Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network Active CN109873447B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910124910.3A CN109873447B (en) 2019-02-19 2019-02-19 Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910124910.3A CN109873447B (en) 2019-02-19 2019-02-19 Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network

Publications (2)

Publication Number Publication Date
CN109873447A CN109873447A (en) 2019-06-11
CN109873447B true CN109873447B (en) 2022-10-04

Family

ID=66918949

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910124910.3A Active CN109873447B (en) 2019-02-19 2019-02-19 Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network

Country Status (1)

Country Link
CN (1) CN109873447B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110601211B (en) * 2019-09-23 2022-12-02 国网辽宁省电力有限公司鞍山供电公司 Method for adjusting distribution network voltage and improving distributed power consumption based on SVC
CN111030125B (en) * 2019-12-23 2022-04-12 东南大学 Active power distribution network optimization regulation and control and distributed scheduling system
CN111756075B (en) * 2020-06-29 2021-12-28 国网经济技术研究院有限公司 Method for designing and testing power distribution system examples containing distributed power supply
CN111799847B (en) * 2020-07-16 2022-11-22 国网北京市电力公司 Predictive control method of risk-considering two-stage random model of active power distribution network
CN112260327B (en) * 2020-10-22 2022-09-30 国网河南省电力公司电力科学研究院 New energy reactive power support analysis method based on extra-high voltage alternating current-direct current power grid
CN113516278B (en) * 2021-04-26 2023-08-22 山东大学 Active power distribution network multi-time scale active and reactive power coordination optimization scheduling method and system
CN114118532A (en) * 2021-11-06 2022-03-01 深圳供电局有限公司 Scheduling method and device for island microgrid, computer equipment and storage medium
CN114421474B (en) * 2022-03-30 2022-06-07 南京易司拓电力科技股份有限公司 Power-voltage sensitivity estimation method between distribution network nodes

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106786806A (en) * 2016-12-15 2017-05-31 国网江苏省电力公司南京供电公司 A kind of power distribution network active reactive based on Model Predictive Control coordinates regulation and control method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106786806A (en) * 2016-12-15 2017-05-31 国网江苏省电力公司南京供电公司 A kind of power distribution network active reactive based on Model Predictive Control coordinates regulation and control method

Also Published As

Publication number Publication date
CN109873447A (en) 2019-06-11

Similar Documents

Publication Publication Date Title
CN109873447B (en) Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network
Ma et al. Optimal allocation of hybrid energy storage systems for smoothing photovoltaic power fluctuations considering the active power curtailment of photovoltaic
CN109687510B (en) Uncertainty-considered power distribution network multi-time scale optimization operation method
CN105846461B (en) Control method and system for large-scale energy storage power station self-adaptive dynamic planning
CN110690732B (en) Photovoltaic reactive power partition pricing power distribution network reactive power optimization method
CN110581571A (en) dynamic optimization scheduling method for active power distribution network
CN109378861B (en) Robust optimization scheduling method for active power distribution network considering time-space correlation
CN114597969B (en) Power distribution network double-layer optimization method considering intelligent soft switch and virtual power plant technology
CN103490426A (en) Comprehensive control method for voltage stabilization of large-scale wind electricity switching-in weak sending end electrical power system
CN109378864B (en) Source-network-load coordination optimization control method based on new energy consumption
CN107332252B (en) Comprehensive low-voltage treatment method for power distribution network considering generalized reactive power source
CN108933448B (en) Coordination control method and system for medium and low voltage distribution network containing photovoltaic power supply
CN108599259A (en) A kind of micro-capacitance sensor active operational decisions method based on sensitivity analysis
Guo et al. Two‐stage optimal MPC for hybrid energy storage operation to enable smooth wind power integration
CN116365506A (en) Energy-saving and loss-reducing optimization method and system for active power distribution network containing distributed photovoltaic
CN116316834A (en) Robust optimization method, system, equipment and storage medium for photovoltaic output spatial correlation power distribution network
CN111614110A (en) Receiving-end power grid energy storage optimization configuration method based on improved multi-target particle swarm optimization
Zhang et al. Coordinated voltage control of active distribution networks with photovoltaic and power to hydrogen
Xu et al. A hybrid energy storage strategy based on multivariable fuzzy coordinated control of photovoltaic grid‐connected power fluctuations
Pippi et al. Assessing the provision of ancillary services considering BES capacity degradation
CN114498740B (en) Energy storage configuration optimization method for normal and fault voltage fluctuation suppression of photovoltaic access power grid
CN115663791A (en) Intelligent power distribution network multi-target adaptive scheduling method based on time-varying property of operating environment
CN104836242A (en) Reactive power optimization method for regional power grid
CN115036931A (en) Active power grid reactive voltage affine adjustable robust optimization method and device
Eissa et al. A novel approach for optimum allocation of Flexible AC Transmission Systems using Harmony Search technique

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