CN116362622B - Power grid power supply capacity assessment method and device under extreme meteorological conditions - Google Patents
Power grid power supply capacity assessment method and device under extreme meteorological conditions Download PDFInfo
- Publication number
- CN116362622B CN116362622B CN202310441512.0A CN202310441512A CN116362622B CN 116362622 B CN116362622 B CN 116362622B CN 202310441512 A CN202310441512 A CN 202310441512A CN 116362622 B CN116362622 B CN 116362622B
- Authority
- CN
- China
- Prior art keywords
- representing
- power grid
- fault rate
- load
- follows
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 14
- 238000012546 transfer Methods 0.000 claims abstract description 57
- 238000011156 evaluation Methods 0.000 claims abstract description 25
- 238000013461 design Methods 0.000 claims description 5
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 abstract 1
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Power Engineering (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention belongs to the field of safe operation of power grids, and particularly relates to a power supply capacity assessment method and device for a power grid under extreme meteorological conditions. The evaluation method comprises the following steps: analyzing the running state of the power grid under extreme meteorological conditions, and establishing a fault rate hybrid model; calculating a mixed fault rate, judging whether important loads need to be transferred or not, and determining transfer priority; and constructing a power grid evaluation function under extreme meteorological conditions, taking the minimum electric energy shortage of the important load as a target, establishing an objective function, and carrying out electric energy transfer on the important load according to the transfer priority. The invention not only can determine the mixed fault rate of the important load, but also can realize the electric energy supply of the important load, ensure the normal electricity utilization of the important load in the power grid, provide a novel guiding method for the safe operation of the power grid and improve the operation reliability of the power grid.
Description
Technical Field
The invention belongs to the field of safe operation of power grids, and particularly relates to a power supply capacity assessment method and device for a power grid under extreme meteorological conditions.
Background
The safe operation of the power grid can not only improve the operation efficiency of the power system and ensure the sharing of high-quality power users, but also promote the economic development and maximize the benefit. However, in recent years, various small-probability extreme weather seriously threatens the normal operation of the power grid, and causes great loss to national economy, so that the operation capability of the power grid under the extreme weather condition is improved, and the method has very important significance for the safe operation of a power system and the development of the national economy. The probability of the outage of the important load close to the extreme weather is far greater than that of the important load far from the extreme weather, and how to measure the outage probability of the important load under the extreme weather is an important premise for ensuring the power supply of the important load.
Disclosure of Invention
Aiming at the defects of the existing method, the invention provides a power supply capacity evaluation method and device for a power grid under extreme meteorological conditions.
In order to achieve the above purpose, the invention adopts the following technical scheme:
a power supply capacity evaluation method for a power grid under extreme meteorological conditions comprises the following steps:
analyzing the running state of the power grid under extreme meteorological conditions, and establishing a fault rate hybrid model;
calculating a mixed fault rate, judging whether important loads need to be transferred or not, and determining transfer priority;
and constructing a power grid evaluation function under extreme meteorological conditions, taking the minimum electric energy shortage of the important load as a target, establishing an objective function, and carrying out electric energy transfer on the important load according to the transfer priority.
Preferably, the fault rate hybrid model includes a wind speed impact power grid outage fault rate model, an overload fault rate model, a line fault rate model and other fault rate models, and the definitions are respectively as follows:
the wind speed influences the power grid to shut down the fault rate model:
wherein f w The wind speed influences the outage fault rate of the power grid, V represents the wind speed, and V represents the design wind speed;
the overload fault rate model is as follows:
f g =f g1 +f g2
wherein f g Indicating overload failure rate, f g1 Indicating the current overload rate of the branch, f g2 Representing a voltage threshold crossing rate;
the line fault rate model is as follows:
f x =f x1 +f x2
wherein f x Representing the failure rate of the line, f x1 Represents the broken line fault rate, f x2 Representing the reverse tower fault rate;
other failure rate models are:
wherein f q Representing other failure rate, x q Representing the variable, mu q Mean value, sigma q Representing standard deviation.
Preferably, the calculating the hybrid fault rate, determining whether the important load needs to be diverted, includes:
the mixed failure rate P f The calculation is performed according to the following formula:
P f =1-(1-f w )(1-f g )(1-f x )(1-f q )
determining a mixed failure rate P f Whether C is satisfied or not, C represents a critical value, and if so, important loads are determined to need transferring; if not, it indicates that the important load does not need to be transferred.
Preferably, the determining the priority of forwarding includes:
all disaster scenes are simulated in advance by using a Monte Carlo method, all possible transfer success paths are analyzed, and the running time of each possible transfer success path is recorded as (l, t) = { (l) i ,t i ) I=1, 2, …, n, where L represents the path set, T represents the run-time set, and transfer priority is determined from the run-time.
Preferably, the power grid evaluation function under the extreme meteorological conditions is:
F(λ,P,f)=λ w P w f w +λ g P g f g +λ x P x f x +λ q P q f q
wherein F (λ, P, F) represents an evaluation function, λ w Representing the coefficient of failure rate of wind speed influencing power grid outage, P w Indicating that important loads lack of electric energy value when wind speed influences power grid outage, f w The wind speed is expressed to influence the outage fault rate of the power grid; lambda (lambda) g Representing the overload failure rate coefficient, P g Indicating that the important load lacks an electric energy value when in overload fault, f g Indicating an overload failure rate; lambda (lambda) x Representing the line fault rate coefficient, P x Indicating the lack of electric energy value of important load during line fault, f x Representing line faultsA rate; lambda (lambda) q Representing other failure rate coefficients, P q Indicating that the important load lacks power values at other faults, and fq indicates other fault rates.
Preferably, the objective function is:
minF(λ,P,f)
s.t.λ i ∈[0,1]and (2) and
P i ∈[0,P i max ]
f i ∈[0,1]
wherein F (lambda, P, F) represents an evaluation function, lambda i represents a power grid outage fault rate coefficient caused by an ith scene, nF represents the total number of scenes, pi represents an important load lack electric energy value when the fault of the ith scene occurs, and P i max Indicating the lack of the maximum power value of the important load in the ith scene, f i Indicating the i-th scene failure rate.
An apparatus for evaluating power supply capacity of a power grid under extreme meteorological conditions, comprising:
the model construction module is used for analyzing the running state of the power grid under the extreme meteorological conditions and establishing a fault rate hybrid model;
the transfer supply determining module calculates the mixed fault rate, judges whether important loads need transfer supply or not, and determines transfer supply priority;
and the electric energy transfer module is used for constructing a power grid evaluation function under extreme weather conditions, taking the minimum electric energy shortage of the important load as a target, establishing an objective function, and carrying out electric energy transfer on the important load according to transfer priority.
The invention has the positive beneficial effects that:
1. although extreme weather is a small probability event, the probability of occurrence still exists, and once the occurrence happens, the safety operation of the power grid is seriously influenced, the normal development of national economy is influenced, and how to measure the probability of the shutdown of an important load under the extreme weather condition is an important premise for ensuring the power supply of the important load. The invention provides a power supply capacity assessment method of a power grid under extreme weather conditions, which comprises the steps of analyzing the running state of the power grid under the extreme weather conditions, establishing a fault rate mixing model, calculating each fault rate according to the fault rate mixing model, calculating to obtain a mixed fault rate, judging whether important loads in the power grid under the extreme weather conditions need to be supplied or not according to the magnitude of the mixed fault rate, if the mixed fault rate is larger than a critical value, indicating that the possibility of the important loads to be shut down is larger, and if the mixed fault rate is larger than the critical value, considering that the important loads need to be supplied, otherwise, indicating that the possibility of the important loads to be shut down is smaller, and then further determining the priority of the supply. According to the invention, through constructing an evaluation function, the outage probability of important loads in the power grid when extreme weather conditions occur is evaluated, an objective function which aims at minimizing the electric energy shortage of the important loads is established, and electric energy transfer is carried out on the important loads according to transfer priority.
The method can not only determine the mixed fault rate of the important load, but also realize the electric energy supply of the important load, ensure the normal operation of the important load in the power grid, provide a new method for evaluating the fault rate of the power grid under extreme meteorological conditions, and improve the operation reliability of the power grid.
Drawings
FIG. 1 is a flow chart of a method for evaluating power supply capacity of a power grid under extreme meteorological conditions.
Fig. 2 is a 14-node grid system diagram.
FIG. 3 is a graph of the mixed failure rate of each significant load node shutdown as typhoons pass.
FIG. 4 is a graph of power supplied from each load node.
Detailed Description
The invention will be further illustrated with reference to a few specific examples.
Example 1
The power supply capacity evaluation method for the power grid under the extreme meteorological conditions, referring to fig. 1, comprises the following steps:
step one: analyzing the running state of the power grid under extreme meteorological conditions, and establishing a fault rate hybrid model;
the fault rate hybrid model comprises a wind speed influence power grid outage fault rate model, an overload fault rate model, a line fault rate model and other fault rate models, and is defined as follows:
the wind speed influences the power grid to shut down the fault rate model:
wherein f w The wind speed affects the failure rate of the power grid, V represents the wind speed, and V represents the design wind speed.
The overload fault rate model is as follows:
f g =f g1 +f g2
wherein f g Indicating overload failure rate, f g1 Indicating the current overload rate of the branch, f g2 Represents the voltage out-of-limit rate, n I And n U Respectively representing the number of overload branches and the number of out-of-limit voltage buses, I i Representing the current in the branch i,represents the maximum value, U j The j-th bus voltage value is shown.
The line fault rate model is as follows:
f x =f x1 +f x2
wherein f x Representing the failure rate of the line, f x1 Represents the broken line fault rate, f x2 Indicating the reverse tower failure rate.
During meteorological disasters, the main reason for the disconnection of the line is that the acting force of disaster causing factors on the line exceeds the self-bearing limit. The wind load born by the wire in the horizontal direction corresponding to any wind speed v acting on the wire is as follows:
H s =0.61275ημ z μ sc β c dLbv 2 sin 2 α
wherein eta represents the wind pressure non-uniformity coefficient of the wire and mu z Represents the coefficient of variation, mu, of the wind pressure altitude sc Representing the body form factor, beta c The wind load adjusting coefficient of the wire is represented by d, the outer diameter of the wire is represented by L, the horizontal span of the tower is represented by b, the wind load increasing coefficient during icing is represented by v, the horizontal wind speed is represented by v, and the included angle between the wire and the wind direction is represented by alpha.
The dead weight of the wire is as follows:
H g =Ψm g Lg
wherein ψ represents the wire division number, m g The mass of the wire per unit length is shown, and g is the gravitational acceleration.
The total load of the lead is as follows:
the specific load of the lead is as follows:
the broken line fault rate is as follows:
wherein K is fx1 Representing the fault parameters of broken lines, F hs Represents the horizontal stress of the wire under wind load, l represents the line span,representing the overhead line elevation angle.
The unbalanced force born by the two sides of the pole tower is as follows:
wherein n is D Indicating the number of wires on the tower, F 1 And F 2 Respectively representing horizontal stress of wires at two ends of the tower, wherein theta represents an included angle of the wires at two ends of the tower, and for a right-angle pole theta=0.
The wind load born by the tower in the horizontal direction is as follows:
wherein A is gs Representing the body form factor of the tower and ρ gs Representing the wind vibration coefficient.
The failure rate of the inverted tower is as follows:
other failure rate models are:
wherein f q Representing other failure rate, x q Representing the variable, mu q Mean value, sigma q Representing standard deviation.
Step two: calculating a mixed fault rate, judging whether important loads need to be transferred or not, and determining transfer priority;
calculating a hybrid failure rate P f :
P f =1-(1-f w )(1-f g )(1-f x )(1-f q )
Determining a mixed failure rate P f Whether the value of C is more than C is met or not, wherein C represents a critical value, if the value is met, the important load is determined to need to be transferred according to experience, and if the value is not met, the important load is not required to be transferred;
simulating all disaster scenes in advance by using a Monte Carlo method, analyzing all possible transfer success paths and each possible transfer success pathRadial run time, noted as (l, t) = { (l) i ,t i )/l i E L, ti e T, i=1, 2, …, n, where L represents the path set and T represents the run-time set; the transfer priority is determined according to the running time, and the smaller t is, the higher the priority is.
Thirdly, constructing a power grid evaluation function under extreme meteorological conditions, taking the minimum of important load electric energy deficiency as a target, establishing an objective function, and carrying out electric energy transfer on important loads according to transfer priority;
the power grid evaluation function is:
F(λ,P,f)=λ w P w f w +λ g P g f g +λ x P x f x +λ q P q f q
wherein F (λ, P, F) represents the constructed evaluation function, λ w Representing the coefficient of failure rate of wind speed influencing power grid outage, P w Indicating that important loads lack of electric energy value when wind speed influences power grid outage, f w The wind speed is expressed to influence the outage fault rate of the power grid; lambda (lambda) g Representing the overload failure rate coefficient, P g The method comprises the steps that when an overload fault occurs, an important load lacks an electric energy value, and fg represents an overload fault rate; lambda (lambda) x Representing the line fault rate coefficient, P x The method comprises the steps that when a line fails, an important load lacks an electric energy value, and fx represents the line failure rate; lambda (lambda) q Representing other failure rate coefficients, P q Indicating the lack of power value of important load in other faults, f q Indicating other failure rates.
The objective function is:
minF(λ,P,f)
s.t.λ i ∈[0,1]and (2) and
P i ∈[0,P i max ]
f i ∈[0,1]
wherein F (λ, P, F) represents an evaluation function, λ i Representing the power failure rate coefficient of the power grid caused by the ith scene, n F Representing the total number of scenes, P i Indicating the lack of electric energy value, P, of important load when the ith scene fault occurs i max Indicating the lack of the maximum power value of the important load in the ith scene, f i Indicating the i-th scene failure rate.
An evaluation device for power supply capacity of a power grid under extreme meteorological conditions, which realizes the above evaluation method, comprises the following steps:
the model construction module is used for analyzing the running state of the power grid under the extreme meteorological conditions and establishing a fault rate hybrid model;
the transfer supply determining module calculates the mixed fault rate, judges whether important loads need transfer supply or not, and determines transfer supply priority;
and the electric energy transfer module is used for constructing a power grid evaluation function under extreme weather conditions, taking the minimum electric energy shortage of the important load as a target, establishing an objective function, and carrying out electric energy transfer on the important load according to the transfer priority.
The following description will take a 14-node power grid system diagram shown in fig. 2 as an example, where node 1, node 8 and node 10 are transformer nodes; node 12, node 6, node 4 are important load nodes. The typhoon is taken as an extreme weather to explain, when the typhoon moves according to a blue double solid line, the wind speed is set to be 25km/s when the typhoon arrives at a power grid, the wind speed is set to be 15km/s when the typhoon leaves the power grid, the design wind speed of the power grid is 14km/s, C takes 0.5, the total load value of a node 12 is 12MW, the total load of a node 6 is 10MW, the total load of a node 4 is 8MW, and the probability that the typhoon stops running through important load nodes 12, 6 and 4 of the power grid is obtained by constructing a fault rate hybrid model as shown in figure 3.
As can be seen from fig. 3, the outage probability of the important load node 12 is 0.7017, the outage probability of the important load node 6 is 0.5506, and the values of the two nodes are both greater than C, which indicates that the outage probability of the two nodes is relatively high; the outage probability of the important load node 4 is 0.3395, which is smaller than C, and the outage probability is smaller. Because the outage probability of the important load node 12 and the important load node 6 is high, two important loads need to be transferred; the transfer priority is then determined and the results are shown in table 1.
Table 1 priority of transfer for important load node 12 and node 6
Important load node | Priority of transfer 1 | Priority of transfer 2 | Priority of transfer 3 | Priority of transfer 4 |
6 | L3-L4 | L11-L4 | L8-L6-L5 | L1-L2-L7-L6-L5 |
12 | L10-L14 | L9-L14 | L3-L11-L12-L13 | L1-L2-L7-L8-L9-L14 |
As can be seen from table 1, four possible transfer paths are obtained for each of the critical load node 6 and the critical load node 12. For the important load node 6, the preferred transfer path is L3-L4, then L11-L4, again L8-L6-L5, and finally L1-L2-L7-L6-L5; for the critical load node 12, the preferred transfer paths are L10-L14, followed by L9-L14, followed by L3-L11-L12-L13, and finally followed by L1-L2-L7-L8-L9-L14.
By utilizing the evaluation function and the objective function constructed by the invention, the electric energy needed to be transferred by the important load node 12 and the node 6 is obtained as shown in the figure 4, and the node 12 can start to break down when typhoons pass through the power grid for 3.3 hours, and the electric energy minimum value is only 3.25MW; the node 6 fails when typhoons pass through the power grid for 5 hours, the minimum value of electric energy is 4.19MW, the important load nodes 12 are influenced by typhoons most, the electric energy needed to be transferred is the most, and the electric energy transfer is carried out on the important load nodes according to the transfer priority.
Finally, it is noted that the above-mentioned embodiments are merely for illustrating the technical solution of the present invention, and that other modifications and equivalents thereof by those skilled in the art should be included in the scope of the claims of the present invention without departing from the spirit and scope of the technical solution of the present invention.
Claims (4)
1. The power supply capacity evaluation method for the power grid under the extreme meteorological conditions is characterized by comprising the following steps of:
analyzing the running state of the power grid under extreme meteorological conditions, and establishing a fault rate hybrid model;
calculating a mixed fault rate, judging whether important loads need to be transferred or not, and determining transfer priority;
constructing a power grid evaluation function under extreme meteorological conditions, taking the minimum electric energy shortage of an important load as a target, establishing an objective function, and carrying out electric energy transfer on the important load according to transfer priority;
the fault rate hybrid model comprises a wind speed influence power grid outage fault rate model, an overload fault rate model, a line fault rate model and other fault rate models, and is defined as follows:
the wind speed influences the power grid to shut down the fault rate model:
wherein f w The wind speed influences the outage fault rate of the power grid, V represents the wind speed, and V represents the design wind speed;
the overload fault rate model is as follows:
f g =f g1 +f g2
wherein f g Indicating overload failure rate, f g1 Indicating the current overload rate of the branch, f g2 Represents the voltage out-of-limit rate, n I And n U Respectively representing the number of overload branches and the number of out-of-limit voltage buses, I i Representing the current in the branch i,represents the maximum value, U j Representing the voltage value of the j-th bus;
the line fault rate model is as follows:
f x =f x1 +f x2
wherein f x Representing the failure rate of the line, f x1 Represents the broken line fault rate, f x2 Representing the reverse tower fault rate;
the wind load born by the wire in the horizontal direction corresponding to any wind speed v acting on the wire is as follows:
H s =0.61275ημ z μ sc β c dLbv 2 sin 2 α
wherein eta represents the wind pressure non-uniformity coefficient of the wire and mu z Represents the coefficient of variation, mu, of the wind pressure altitude sc Representing the body form factor, beta c The wind load adjusting coefficient of the wire is represented by d, the outer diameter of the wire is represented by L, the horizontal span of the tower is represented by b, the wind load increasing coefficient during icing is represented by v, the horizontal wind speed is represented by v, and the included angle between the wire and the wind direction is represented by alpha;
the dead weight of the wire is as follows:
H g =Ψm g Lg
wherein ψ represents the wire division number, m g Representing the mass of a wire per unit length, g representing the gravitational acceleration;
the total load of the lead is as follows:
the specific load of the lead is as follows:
the broken line fault rate is as follows:
wherein K is fx1 Representing the fault parameters of broken lines, F hs Represents the horizontal stress of the wire under wind load, l represents the line span,representing the overhead line height difference angle;
the unbalanced force born by the two sides of the pole tower is as follows:
wherein n is D Indicating the number of wires on the tower, F 1 And F 2 Respectively representing horizontal stress of wires at two ends of the tower, wherein theta represents an included angle of the wires at two ends of the tower, and for a right-angle pole theta=0;
the wind load born by the tower in the horizontal direction is as follows:
wherein A is gs Representing the body form factor of the tower and ρ gs Representing wind vibration coefficients;
the failure rate of the inverted tower is as follows:
other failure rate models are:
wherein f q Representing other failure rate, x q Representing the variable, mu q Mean value, sigma q Representing standard deviation;
the calculating the mixed fault rate and judging whether the important load needs to be transferred or not comprises the following steps:
the mixed failure rate P f The calculation is performed according to the following formula:
P f =1-(1-f w )(1-f g )(1-f x )(1-f q )
determining a mixed failure rate P f Whether C is satisfied or not, C represents a critical value, and if so, important loads are determined to need transferring; if the load is not satisfied, indicating that the important load is not required to be transferred;
the power grid evaluation function under the extreme meteorological conditions is as follows:
F(λ,P,f)=λ w P w f w +λ g P g f g +λ x P x f x +λ q P q f q
wherein F (λ, P, F) represents an evaluation function, λ w Representing the coefficient of failure rate of wind speed influencing power grid outage, P w Indicating that important loads lack of electric energy value when wind speed influences power grid outage, f w The wind speed is expressed to influence the outage fault rate of the power grid; lambda (lambda) g Representing the overload failure rate coefficient, P g Representing the critical load in the event of overload failureLack of electric energy value, f g Indicating an overload failure rate; lambda (lambda) x Representing the line fault rate coefficient, P x Indicating the lack of electric energy value of important load during line fault, f x Representing a line failure rate; lambda (lambda) q Representing other failure rate coefficients, P q Indicating the lack of power value of important load in other faults, f q Representing other failure rates;
the objective function is:
minF(λ,P,f)
s.t.λ i ∈[0,1]and (2) and
P i ∈[0,P i max ]
f i ∈[0,1]
wherein F (λ, P, F) represents an evaluation function, λ i Representing the power failure rate coefficient of the power grid caused by the ith scene, n F Representing the total number of scenes, P i Indicating the lack of electric energy value, P, of important load when the ith scene fault occurs i max Indicating the lack of the maximum power value of the important load in the ith scene, f i Indicating the i-th scene failure rate.
2. The method for evaluating power supply capacity of a power grid under extreme weather conditions as set forth in claim 1, wherein said determining a transfer priority comprises:
all disaster scenes are simulated in advance by using a Monte Carlo method, all possible transfer success paths are analyzed, and the running time of each possible transfer success path is recorded as (l, t) = { (l) i ,t i )/l i ∈L,t i E T }, i=1, 2, …, n, where L represents the set of paths, T represents the set of runtimes, and the transfer priority is determined from the runtimes.
3. An apparatus for evaluating power supply capacity of a power grid under extreme meteorological conditions, comprising:
the model construction module is used for analyzing the running state of the power grid under the extreme meteorological conditions and establishing a fault rate hybrid model;
the transfer supply determining module calculates the mixed fault rate, judges whether important loads need transfer supply or not, and determines transfer supply priority;
the power transfer module is used for constructing a power grid evaluation function under extreme weather conditions, taking the minimum power shortage of an important load as a target, establishing an objective function, and carrying out power transfer on the important load according to transfer priority;
the fault rate hybrid model comprises a wind speed influence power grid outage fault rate model, an overload fault rate model, a line fault rate model and other fault rate models, and is defined as follows:
the wind speed influences the power grid to shut down the fault rate model:
wherein f w The wind speed influences the outage fault rate of the power grid, V represents the wind speed, and V represents the design wind speed;
the overload fault rate model is as follows:
f g =f g1 +f g2
wherein f g Indicating overload failure rate, f g1 Indicating the current overload rate of the branch, f g2 Represents the voltage out-of-limit rate, n I And n U Respectively representing the number of overload branches and the number of out-of-limit voltage buses, I i Representing the current in the branch i,represents the maximum value, U j Representing the voltage value of the j-th bus;
the line fault rate model is as follows:
f x =f x1 +f x2
wherein f x Representing the failure rate of the line, f x1 Represents the broken line fault rate, f x2 Representing the reverse tower fault rate;
the wind load born by the wire in the horizontal direction corresponding to any wind speed v acting on the wire is as follows:
H s =0.61275ημ z μ sc β c dLbv 2 sin 2 α
wherein eta represents the wind pressure non-uniformity coefficient of the wire and mu z Represents the coefficient of variation, mu, of the wind pressure altitude sc Representing the body form factor, beta c The wind load adjusting coefficient of the wire is represented by d, the outer diameter of the wire is represented by L, the horizontal span of the tower is represented by b, the wind load increasing coefficient during icing is represented by v, the horizontal wind speed is represented by v, and the included angle between the wire and the wind direction is represented by alpha;
the dead weight of the wire is as follows:
H g =Ψm g Lg
wherein ψ represents the wire division number, m g Representing the mass of a wire per unit length, g representing the gravitational acceleration;
the total load of the lead is as follows:
the specific load of the lead is as follows:
the broken line fault rate is as follows:
wherein K is fx1 Representing the fault parameters of broken lines, F hs Represents the horizontal stress of the wire under wind load, l represents the line span,representing the overhead line height difference angle;
the unbalanced force born by the two sides of the pole tower is as follows:
wherein n is D Indicating the number of wires on the tower, F 1 And F 2 Respectively representing horizontal stress of wires at two ends of the tower, wherein theta represents an included angle of the wires at two ends of the tower, and for a right-angle pole theta=0;
the wind load born by the tower in the horizontal direction is as follows:
wherein A is gs Representing the body form factor of the tower and ρ gs Representing wind vibration coefficients;
the failure rate of the inverted tower is as follows:
other failure rate models are:
wherein f q Representing other failure rate, x q Representing the variable, mu q Mean value, sigma q Representing standard deviation;
the calculating the mixed fault rate and judging whether the important load needs to be transferred or not comprises the following steps:
the saidHybrid failure rate P f The calculation is performed according to the following formula:
P f =1-(1-f w )(1-f g )(1-f x )(1-f q )
determining a mixed failure rate P f Whether C is satisfied or not, C represents a critical value, and if so, important loads are determined to need transferring; if not, it indicates that the important load does not need to be transferred.
4. A grid power supply capability assessment apparatus under extreme weather conditions according to claim 3, wherein said determining a transfer priority comprises:
all disaster scenes are simulated in advance by using a Monte Carlo method, all possible transfer success paths are analyzed, and the running time of each possible transfer success path is recorded as (l, t) = { (l) i ,t i )/l i ∈L,t i E T }, i=1, 2, …, n, where L represents the set of paths, T represents the set of runtimes, and the transfer priority is determined from the runtimes.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310441512.0A CN116362622B (en) | 2023-04-23 | 2023-04-23 | Power grid power supply capacity assessment method and device under extreme meteorological conditions |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310441512.0A CN116362622B (en) | 2023-04-23 | 2023-04-23 | Power grid power supply capacity assessment method and device under extreme meteorological conditions |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116362622A CN116362622A (en) | 2023-06-30 |
CN116362622B true CN116362622B (en) | 2023-11-21 |
Family
ID=86933523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310441512.0A Active CN116362622B (en) | 2023-04-23 | 2023-04-23 | Power grid power supply capacity assessment method and device under extreme meteorological conditions |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116362622B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107230015A (en) * | 2017-05-25 | 2017-10-03 | 天津大学 | A kind of power distribution network toughness appraisal procedure based on system information entropy |
CN108876194A (en) * | 2018-07-16 | 2018-11-23 | 国网福建省电力有限公司 | Power distribution network methods of risk assessment under a kind of typhoon disaster scene |
CN109617085A (en) * | 2019-01-28 | 2019-04-12 | 深圳供电局有限公司 | Rapid load shedding optimization control method in fault state |
CN111126841A (en) * | 2019-12-24 | 2020-05-08 | 东南大学 | Power distribution network load to energy conversion power evaluation method based on random power flow |
WO2022022101A1 (en) * | 2020-07-30 | 2022-02-03 | 国网甘肃省电力公司电力科学研究院 | Method for controlling coordinated peak regulation of emergency source network in sending-end grid fault state |
CN115587649A (en) * | 2022-09-16 | 2023-01-10 | 国网黑龙江省电力有限公司 | Power transmission line maintenance sequence optimization decision method under extreme ice and snow disasters |
-
2023
- 2023-04-23 CN CN202310441512.0A patent/CN116362622B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107230015A (en) * | 2017-05-25 | 2017-10-03 | 天津大学 | A kind of power distribution network toughness appraisal procedure based on system information entropy |
CN108876194A (en) * | 2018-07-16 | 2018-11-23 | 国网福建省电力有限公司 | Power distribution network methods of risk assessment under a kind of typhoon disaster scene |
CN109617085A (en) * | 2019-01-28 | 2019-04-12 | 深圳供电局有限公司 | Rapid load shedding optimization control method in fault state |
CN111126841A (en) * | 2019-12-24 | 2020-05-08 | 东南大学 | Power distribution network load to energy conversion power evaluation method based on random power flow |
WO2022022101A1 (en) * | 2020-07-30 | 2022-02-03 | 国网甘肃省电力公司电力科学研究院 | Method for controlling coordinated peak regulation of emergency source network in sending-end grid fault state |
CN115587649A (en) * | 2022-09-16 | 2023-01-10 | 国网黑龙江省电力有限公司 | Power transmission line maintenance sequence optimization decision method under extreme ice and snow disasters |
Non-Patent Citations (3)
Title |
---|
何嘉兴 等.极端灾害下考虑应急转供的配网停电过程仿真方法.智慧电力.2020,第48卷(第4期),第104-111页. * |
极端灾害下考虑应急转供的配网停电过程仿真方法;何嘉兴 等;智慧电力;第48卷(第4期);第104-111页 * |
计及转供电概率的配电网可靠性评估;李志铿 等;南方能源建设;第1卷(第01期);30-33 * |
Also Published As
Publication number | Publication date |
---|---|
CN116362622A (en) | 2023-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108470231B (en) | Power distribution network distributed energy storage site selection and volume fixing method considering energy storage system quantization characteristics | |
WO2023184700A1 (en) | Battery system charging and discharging control method based on dynamic reconfigurable battery network | |
CN110535144A (en) | The intelligent distribution network toughness quantitative analysis method of the load containing polymorphic type under dusty wind weather | |
CN113489003B (en) | Source network coordination planning method considering wind-light-water integrated complementary operation | |
CN115456304A (en) | Offshore wind farm reliability index calculation method and device considering typhoon influence | |
CN111552925A (en) | Method for constructing comprehensive evaluation index system of alternating current-direct current hybrid power distribution network | |
CN111313420B (en) | Power distribution network elastic lifting method and system based on multi-energy coordination in extreme weather | |
CN117578706A (en) | Multi-level load management system and method based on power grid | |
Zhang et al. | Boosting the power grid resilience under typhoon disasters by coordinated scheduling of wind energy and conventional generators | |
CN116362622B (en) | Power grid power supply capacity assessment method and device under extreme meteorological conditions | |
CN113541211A (en) | Method and system for determining toughness of alternating current-direct current power distribution network | |
CN115809836B (en) | Method for planning toughness of power distribution network by considering distributed energy storage emergency power supply capacity | |
CN117117918A (en) | Source network load storage optimization operation method and device | |
CN115765001A (en) | Power supply reliability calculation method for novel power system | |
CN116545009A (en) | Distributed wind-light resource admission capacity planning method and system | |
CN114268172B (en) | Multi-type energy storage operation control method | |
CN107886227B (en) | Method for evaluating disaster resistance improvement degree of power distribution network | |
CN114357671A (en) | Distribution network elastic lifting planning configuration method and system based on weak link identification | |
CN109149636A (en) | Wind power integration point transient stability based on Dynamic Security Region influences research method | |
CN115276051A (en) | Receiving-end urban power grid elasticity evaluation method considering new energy and energy storage response characteristics | |
CN107834542B (en) | Effectiveness analysis method for access of extra-high voltage power grid to receiving-end power grid | |
CN113315140A (en) | Offshore wind power reactive compensation coordination control method | |
CN113554272B (en) | Full life cycle supplement method for electrochemical energy storage power station | |
CN110445128B (en) | Sensitivity-based real-time evaluation method for new energy consumption capability of regional power grid | |
CN113900031B (en) | SOC safety verification method after energy storage system is accessed to AGC |
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 |