CN114493365A - Method for evaluating cascading failure vulnerability of power system including wind power plant - Google Patents
Method for evaluating cascading failure vulnerability of power system including wind power plant Download PDFInfo
- Publication number
- CN114493365A CN114493365A CN202210198717.6A CN202210198717A CN114493365A CN 114493365 A CN114493365 A CN 114493365A CN 202210198717 A CN202210198717 A CN 202210198717A CN 114493365 A CN114493365 A CN 114493365A
- Authority
- CN
- China
- Prior art keywords
- power
- line
- total
- cascading failure
- vulnerability
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000035515 penetration Effects 0.000 claims abstract description 29
- 238000004088 simulation Methods 0.000 claims abstract description 8
- 230000005540 biological transmission Effects 0.000 claims abstract description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000012905 input function Methods 0.000 claims description 3
- 241000764238 Isis Species 0.000 claims description 2
- 239000000654 additive Substances 0.000 claims description 2
- 230000000996 additive effect Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 abstract description 15
- 238000010248 power generation Methods 0.000 abstract description 6
- 230000035699 permeability Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 5
- 239000004020 conductor Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- 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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
-
- 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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/388—Islanding, i.e. disconnection of local power supply from the network
-
- 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]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention relates to an evaluation method for cascading failure vulnerability of a power system including a wind power plant, which comprises the following steps of: 1) constructing an optimal power flow model based on a line overload model; 2) under the power grid cascading failure, the line shutdown routing is carried out by combining a thermal stability model and a random model of the power transmission line, and the total trip occupation ratio of the line is obtained; 3) performing cascading failure simulation based on the mixed optimal power flow model and the random model, constructing a power grid vulnerability index considering the wind power uncertainty level and the penetration rate, and evaluating the influence of the power grid vulnerability through the total percentage of line tripping and the total load shedding percentage under the cascading failure. Compared with the prior art, the method effectively solves the problems of uncertainty and permeability of the power grid caused by popularization of renewable energy sources such as wind power generation and the like, quickly and accurately evaluates the vulnerability of the power grid to different uncertainty levels and penetration rates during cascading failures, and is helpful for technicians and decision-makers to more efficiently solve related practical problems.
Description
Technical Field
The invention relates to the technical field of cascading failure vulnerability assessment of a power system, in particular to a cascading failure vulnerability assessment method of a power system including a wind power plant.
Background
While promoting economic development, the global industrial revolution also uses a large amount of limited and environmentally unfriendly fossil energy, which leads to rapid consumption of energy and large carbon emission, and the threat is becoming serious, so that renewable energy sources such as wind energy and solar energy are receiving wide attention due to their environmental protection properties and contribution to saving fossil fuel resources. Meanwhile, the ratio of new energy power generation such as wind power generation in a power system has gradually increased, but the evaluation of cascading failure of the power system is also more challenging.
Large-scale blackouts caused by cascading failures, which originate from strong interdependencies within the grid, weaken the power system when they occur, and may lead to system instability and large-scale blackouts, represent a significant economic and social cost each year. The enormous economic and social impact of such events has prompted people to study the vulnerability of the grid to cascading failures and to find more efficient evaluation methods. In power systems, the penetration of new energy generation can bring the grid closer to its operating limits and introduce a large amount of uncertainty, which comes from the randomness of the new energy, thus changing the dynamic performance of the grid, which is first of all the increase of cascading failures related to wind farms. Therefore, it is becoming more and more important to evaluate the impact of different wind power uncertainty levels and penetration rates under cascading failures on grid vulnerability under cascading failures.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide an evaluation method for cascading failure vulnerability of a power system of a wind power plant.
The purpose of the invention can be realized by the following technical scheme:
a method for evaluating cascading failure vulnerability of a power system of a wind power plant comprises the following steps:
1) constructing an optimal power flow model based on a line overload model;
2) under the power grid cascading failure, a line shutdown path finding is carried out by combining a power transmission line thermal stability model and a random model, and the total trip occupation ratio of the line is obtained;
3) performing cascading failure simulation based on a mixed optimal power flow model and a random model, constructing a power grid vulnerability index considering wind power uncertainty level and penetration rate, identifying the most vulnerable line in the power grid under different wind power uncertainty levels and penetration rates, and evaluating the influence of the power grid vulnerability through the total line tripping percentage mu and the total load shedding percentage eta under cascading failure.
In the step 1), the total cost is taken as an optimization target, the initial scheduling output of the conventional generator is determined according to the load and the wind power prediction information, and an optimal power flow model is established through an online island detection algorithm and an automatic power balance algorithm, so that the method comprises the following steps:
wherein,as a polynomial cost function of the generator i,is the output power of generator i, ngIn order to be the total number of the generators,is the maximum value of the output power of the generator i, Pl jIs the demanded power of load j, nlFor the number of loads, G is the generator set, L is the load set, FijFor the power flow between the node where generator i is located and the node where load j is located,for the tidal current capacity between the node where generator i is located and the node where load j is located, presentation generation and consumptionBalancing the total power of (a).
In step 1), during a cascading failure, all islands generated due to line tripping are identified according to an online island detection algorithm, and an automatic power balance algorithm is adopted to reduce the total load to the minimum so as to limit the maximum capacity of a generator, maintain the power balance of each island, and change the electrical frequency of a system based on the power balance, then:
wherein, PGenFor total emitted power, PLoadFor total power consumption, H is total inertia, f is the electrical frequency of the system, and t is time.
In the step 2), the uncertainty of wind and load in the power flow process is simulated through the overload probability of the line, the line with the minimum overload distance is obtained, and the line is tripped in the case of cascading failure, so that the total tripping percentage of the line is obtained.
The uncertainty modeling specifically comprises the following steps:
21) adopting a sliding window method as an m multiplied by m matrix, and calculating a covariance matrix C of each time stepF(t, τ) there are:
22) The uncertainty of wind and load in the process of line overload probability simulation load flow is introduced into a stochastic model, and the load flow F of the first line is assumed that a power input function P (t) meets Gaussian distributionl(t) satisfies the Gaussian distribution, lineProbability of road overload ρl(t) calculated by the Q function, then there are:
wherein, alFor the standard overload distance of the l line, the expression of the Q function is For the tidal current capacity of the l-th line,respectively is the power flow mean value and the variance of the first line;
23) according to the standard overload distance and overload probability rho of each linelCalculating a load flow F at the time tlAverage overload time in (t)Then there are:
wherein, BWlThe equivalent bandwidth of the flow process of the first line is obtained;
24) when average overload timeGreater than the trip time ttrWhen it is, the step returns to step 21), when the average overload time isIs less than or equal to the tripping time ttrAnd then, tripping off the overload circuit after the tripping time, detecting a new island formed by the tripping circuit by adopting an island detection algorithm, balancing the generated energy and the load of each new island, and reducing the load loss to the maximum extent to obtain the total tripping percentage mu and the total load shedding percentage eta of the circuit.
In the step 24), the tripping time t of the overload line is obtained according to the line thermal stability modeltrThen, there are:
wherein F is the overload line current, FopFor the initial operating trend, FmaxAs a line flow threshold, TthAre thermal time constants related to the type of wire and environmental parameters, specifically wind speed and ambient temperature.
In the step 3), acquiring a vulnerability index psi of the wind power uncertainty level to the power grid under the cascading failure according to the total trip percentage mu and the total load shedding percentage eta of the line1Then, there are:
wherein, WuncIs a weight parameter of the wind power uncertainty level in the system under cascading failure.
In the step 3), the wind power uncertainty level represents the relative error of wind power prediction, and the uncertainty level of the wind power generator is dependent on factorsIs increased by the amount of the additive, wherein,as a new uncertainty index for wind power,For the initial uncertainty index of wind power, as gamma increases, the total trip proportion of a line during cascading failure increases, so that more islands and greater load shedding are formed, namely, as the uncertainty level of wind power increases, the vulnerability of the system is continuously increased.
In the step 3), acquiring a vulnerability index psi of wind power penetration rate to the power grid under cascading failure according to the total trip percentage mu and the total load shedding percentage eta of the line2Then, there are:
wherein, WpenThe weight parameter of the wind power penetration rate in the system under the cascading failure is shown.
In the step 3), the wind power penetration rate is defined as:
wherein,in order to be the total capacity of the wind turbine,for the total capacity of all generators, as alpha is increased, cascading failures are also upgraded, so that the total duty ratio of line tripping is increased, more islands are generated, and larger load shedding is caused, namely, when the wind power penetration rate is increased, the vulnerability of the system is increased.
Compared with the prior art, the invention has the following advantages:
the invention provides a more efficient evaluation method for the vulnerability evaluation of a power grid with cascading failures of a power system of a wind power plant, and the spreading time of the cascading failures is very short, and the uncertainty and permeability of the wind power also bring difficulty to the vulnerability evaluation, therefore, the invention firstly takes the total cost as an optimization target, combines a line overload model, an online island detection algorithm and an automatic power balance algorithm, establishes an optimal power flow model based on direct current power flow, greatly improves the power flow calculation time under the condition of ensuring the high fidelity of the system, then combines a thermal stability model of a power transmission line with a random model, simulates the uncertainty of the wind power and the load in the power flow process through the line overload probability, reduces the randomness in the actual process, provides reliable parameters for the calculation of evaluation indexes, and finally is based on an optimal power flow and random model mixing method, the vulnerability index under the cascading failure of the power grid is established, the wind power uncertainty level and the penetration rate under the cascading failure are calculated, the problem of evaluation of the vulnerability of the power grid under different wind power uncertainty levels and penetration rates is solved, and the evaluation efficiency is greatly improved.
Drawings
Fig. 1 is a flow chart of a power balancing algorithm for newly forming islands in a power grid.
Fig. 2 is a general flow chart of a cascading failure simulation program and vulnerability assessment of a power system based on an optimal power flow & stochastic model hybrid.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments.
Examples
The invention provides an evaluation method for cascading failure vulnerability of a power system including a wind power plant, which comprises the following three steps:
step 2, providing a line outage path-finding method combining a thermal stability model and a random model of the power transmission line under the power grid cascading failure, wherein the method provides a line overload probability to simulate the uncertainty of wind and load in a tidal current process, calculates a line with the minimum overload distance, trips the line in the cascading failure, and finally provides the total tripping percentage of the line for the step 3;
and 3, establishing an evaluation method for the wind power uncertainty level and penetration rate to the power grid vulnerability under the power grid cascading failure, establishing a power grid vulnerability index considering the wind power uncertainty level and the penetration rate based on the optimal power flow and random model mixing method in the step 1 and the step 2, identifying the most vulnerable line in the power grid under different wind power uncertainty levels and penetration rates, and evaluating the influence of the power grid vulnerability through the total line tripping percentage and the total load shedding percentage under the cascading failure.
In step 1, determining the steady-state operation condition of the power grid requires solving a complete power flow equation, so that a large number of solutions exist in the evolution process of cascading faults, repeated solving of the complete nonlinear power flow equation becomes very difficult to calculate, the response time of the cascading faults is short, and the evaluation benefit is also influenced by too long calculation time. Because the evaluation method only concerns the influence of the uncertainty and the penetration rate of the load and the wind power on the vulnerability of the power grid, a complete nonlinear network equation is not necessarily required, a line overload model can be established, direct current load flow calculation is carried out, and the total cost of the system is taken as an optimization target to establish an optimal load flow function:
wherein,as a polynomial cost function of the generator i,is the output power of generator i, ngIn order to be the total number of the generators,is the maximum value of the output power of the generator i, Pl jIs the demanded power of load j, nlFor the number of loads, G is the generator set, L is the load set, FijFor the power flow between the node where generator i is located and the node where load j is located,for the tidal current capacity between the node where generator i is located and the node where load j is located, representing the total power balance of production and consumption.
During a cascading failure, all islands generated due to line tripping are identified based on an online island detection algorithm. Reducing the total load to the minimum through an automatic power balance algorithm to limit the maximum capacity of the generator, maintaining the power balance of each island, and changing the electrical frequency of the system based on the power balance, wherein the specific formula is as follows:
in the formula, PGenFor total emitted power, PLoadFor total power consumption, H is total inertia, f is the electrical frequency of the system, and t is time.
Usually the frequency in the power system is considered as a global parameter and is not affected by a small part. However, when a part of the network becomes islanded, inertia and load balancing depend only on the generator and load within the islanded, where load shedding may be required. In the power balancing algorithm for any island, the total load and the total generation capacity are compared to each other and if the total demand exceeds the maximum available generation, some de-loading is required to maintain power balance. Likewise, if the total demand is less than the current power generation, one or more gensets should reduce their power generation. A flow chart of the automatic power balancing algorithm is shown in fig. 1. Let t be t in the grid0There are k independent islands if t ═ t0A line trip at + Δ t results in the formation of a new island, it is necessary to run an automatic power balancing algorithm on the newly formed island and its separate parent island. Thus, the automatic power balancing algorithm will balance the power generation and load of the two clusters for the next load flow calculation.
In step 2, for modeling of system uncertainty, a covariance matrix C at each time step is first calculated as an m × m matrix using a sliding window methodF(t,τ):
In the formula,is an uncertain item in the first line tide signal at the time t, and a matrix element
The line overload probability is introduced into a stochastic model to simulate the uncertainty of wind and load in the power flow process, and the power input function P (t) is assumed to meet the Gaussian distribution, so that the power flow F of the first linel(t) also satisfies the Gaussian distribution, line overload probabilityCan be calculated using the Q function as shown below:
in the formula,is the standard overload distance of the l-th line, and the Q function is For the tidal current capacity of the l-th line,respectively, the power flow mean value and the variance of the ith line.
Using the standard overload distance (a) of each linel) And overload probability (p)l) Calculating a load flow F at the time tlThe average overload time in (t) is shown as follows:
wherein, BWlIs the equivalent bandwidth of the flow process of the l-th line, which can be calculated using the Spectral Power Density (SPD) of the flow process.
In the line thermal stability model, the trip time of the thermal overload relay is determined according to the maximum allowable current flowing in the conductor, and thermal instability is not caused. Generally, overload protection for high voltage transmission lines has a time dependent trip characteristic, which is determined using a dynamic thermal balance between thermal gain and loss in the conductor. The trip time of the thermal relay is determined by the maximum point temperature and, in order to analyze the dc current flow more effectively, taking into account the initial operating current, the current is replaced by a current flow measured in units, assuming that the voltage in the entire network remains V1.0 p.u., so the trip time is given by:
wherein F is overload line current (p.u.), and FopFor the initial operating trend (p.u.), FmaxAs a line flow threshold, TthIs a thermal time constant associated with the type of conductor and environmental parameters such as wind speed and ambient temperature.
The optimal power flow & stochastic model hybrid method cascading failure simulation and the general flow chart for the vulnerability assessment of the power system are shown in fig. 2.
In step 3, establishing a grid vulnerability index considering wind power uncertainty level and penetration rate, identifying the most vulnerable line in the grid under different wind power uncertainty levels and penetration rates, and evaluating the influence of the grid vulnerability through the total line tripping percentage and the total load shedding percentage under cascading failure, wherein the specific evaluation method is as follows:
(1) wind power uncertainty level under cascading failures
The wind power uncertainty level represents the relative error of wind power prediction, and the uncertainty level of a wind power generator is dependent on factorsIs increased by an increase in whichIs a new uncertainty index of the wind power,is an initial uncertainty index of wind power and belongs towThe uncertainty index is determined by the randomness, the volatility and the intermittence of the output active power in the wind power cascading failure stage.
Based on optimal power flow&Performing cascading failure simulation on a system model established by a random model mixing method to obtain a total tripping proportion mu and a total load shedding percentage eta of a line, and regarding a vulnerability index psi of a wind power uncertainty level under cascading failure to a power grid1As follows:
in the formula, WuncIs a weight parameter of the wind power uncertainty level in the system under cascading failure.
As γ increases, the line trip total occupancy increases during cascading failures, creating more islands and greater load shedding, i.e., as the wind uncertainty level increases, the vulnerability of the system also increases.
(2) Penetration rate of wind power
The wind power penetration rate is defined asWhereinIn order to be the total capacity of the wind turbine,is the total capacity of all generators.
Based on optimal power flow&Performing cascading failure simulation on a system model established by a random model mixing method to obtain a total tripping proportion mu and a total load shedding percentage eta of a line, and regarding a vulnerability index psi of wind power penetration rate under cascading failure to a power grid2As follows:
in the formula, WpenThe weight parameter of the wind power penetration rate in the system under the cascading failure is shown.
As alpha increases, cascading failures will also escalate, causing the line trip total duty to climb, creating more islands and causing greater load shedding, i.e., as wind penetration increases, the vulnerability of the system will also increase.
Claims (10)
1. A method for evaluating cascading failure vulnerability of a power system of a wind power plant is characterized by comprising the following steps:
1) constructing an optimal power flow model based on a line overload model;
2) under the power grid cascading failure, a line shutdown path finding is carried out by combining a power transmission line thermal stability model and a random model, and the total trip occupation ratio of the line is obtained;
3) performing cascading failure simulation based on a mixed optimal power flow model and a random model, constructing a power grid vulnerability index considering wind power uncertainty level and penetration rate, identifying the most vulnerable line in the power grid under different wind power uncertainty levels and penetration rates, and evaluating the influence of the power grid vulnerability through the total line tripping percentage mu and the total load shedding percentage eta under cascading failure.
2. The method for evaluating the cascading failure vulnerability of the power system comprising the wind power plant according to claim 1, characterized in that in the step 1), the total cost is taken as an optimization target, the initial dispatching output of the conventional generator is determined according to the load and the wind power prediction information, and an optimal power flow model is established through an online island detection algorithm and an automatic power balance algorithm, and then:
wherein,as a polynomial cost function of the generator i,is the output power of generator i, ngIn order to be the total number of the generators,is the maximum value of the output power of the generator i,is the demanded power of load j, nlFor the number of loads, G is the generator set, L is the load set, FijFor the power flow between the node at which generator i is located and the node at which load j is located,for the power flow capacity between the node at which generator i is located and the node at which load j is located, representing the total power balance of production and consumption.
3. The method for evaluating cascading failure vulnerability of power system including wind farm according to claim 1, characterized in that in step 1), during cascading failure, all islands generated due to line tripping are identified according to online island detection algorithm, and total load is reduced to minimum by using automatic power balance algorithm to limit maximum capacity of generator, maintain power balance of each island, and change electrical frequency of system based on power balance, then:
wherein, PGenFor total emitted power, PLoadFor total power consumption, H is total inertia, f is the electrical frequency of the system, and t is time.
4. The method for evaluating the cascading failure vulnerability of the power system with the wind farm according to claim 1, characterized in that in the step 2), the line with the minimum overload distance is obtained by simulating the uncertainty of wind power and load in the power flow process through the line overload probability, and the line is tripped in the cascading failure, so that the total tripping proportion of the supply line is obtained.
5. The method for evaluating cascading failure vulnerability of wind farm-containing power system according to claim 4, wherein the uncertainty modeling specifically comprises the following steps:
21) adopting a sliding window method as an m multiplied by m matrix, and calculating a covariance matrix C of each time stepF(t, τ) there are:
22) The uncertainty of wind and load in the process of line overload probability simulation load flow is introduced into a stochastic model, and the load flow F of the first line is assumed that a power input function P (t) meets Gaussian distributionl(t) satisfying Gaussian distribution, line overload probability ρl(t) calculated by the Q function, then there are:
wherein, alFor the standard overload distance of the l line, the expression of the Q function is For the tidal current capacity of the l-th line,respectively is the power flow mean value and the variance of the first line;
23) according to the standard overload distance and overload probability rho of each linelCalculating a load flow F at the time tlAverage overload time in (t)Then there are:
wherein, BWlThe equivalent bandwidth of the flow process of the first line is obtained;
24) when average overload timeGreater than the trip time ttrWhen it is, the step returns to step 21), when the average overload time isIs less than or equal to the tripping time ttrAnd then, tripping off the overload circuit after the tripping time, detecting a new island formed by the tripping circuit by adopting an island detection algorithm, balancing the generated energy and the load of each new island, and reducing the load loss to the maximum extent to obtain the total tripping percentage mu and the total load shedding percentage eta of the circuit.
6. The method for evaluating cascading failure vulnerability of power system of wind power plant according to claim 5, wherein in step 24), the tripping time t of the overload line is obtained according to the line thermal stability modeltrThen, there are:
wherein F is the overload line current, FopFor the initial operating trend, FmaxAs a line flow threshold, TthAre thermal time constants related to the type of wire and environmental parameters, specifically wind speed and ambient temperature.
7. The method for evaluating cascading failure vulnerability of power system including wind farm according to claim 5, characterized in that in step 3), wind power uncertainty level under cascading failure is obtained according to total line trip ratio μ and total load shedding percentage η to obtain grid vulnerability index ψ1Then, there are:
wherein, WuncIs a weight parameter of the wind power uncertainty level in the system under cascading failure.
8. The method for evaluating cascading failure vulnerability of power system of wind power plant according to claim 7, wherein in step 3), wind power uncertainty level represents relative error of wind power prediction, and uncertainty level of wind power generator is dependent on factorsIs increased by the amount of the additive, wherein,is a new uncertainty index of the wind power,for the initial uncertainty index of wind power, as gamma increases, the total trip proportion of a line during cascading failure increases, so that more islands and greater load shedding are formed, namely, as the uncertainty level of wind power increases, the vulnerability of the system is continuously increased.
9. The method for evaluating cascading failure vulnerability of power system including wind farm according to claim 7, characterized in that in step 3), wind power penetration rate to grid vulnerability index psi under cascading failure is obtained according to total line trip ratio μ and total load shedding percentage η2Then, there are:
wherein, WpenThe weight parameter of the wind power penetration rate in the system under the cascading failure is shown.
10. The method for evaluating cascading failure vulnerability of wind power system according to claim 9, wherein in step 3), the wind power penetration rate is defined as:
wherein,in order to be the total capacity of the wind turbine,for the total capacity of all generators, as alpha is increased, cascading failures are also upgraded, so that the total duty ratio of line tripping is increased, more islands are generated, and larger load shedding is caused, namely, when the wind power penetration rate is increased, the vulnerability of the system is increased.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210198717.6A CN114493365A (en) | 2022-03-02 | 2022-03-02 | Method for evaluating cascading failure vulnerability of power system including wind power plant |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210198717.6A CN114493365A (en) | 2022-03-02 | 2022-03-02 | Method for evaluating cascading failure vulnerability of power system including wind power plant |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114493365A true CN114493365A (en) | 2022-05-13 |
Family
ID=81484158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210198717.6A Pending CN114493365A (en) | 2022-03-02 | 2022-03-02 | Method for evaluating cascading failure vulnerability of power system including wind power plant |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114493365A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114925993A (en) * | 2022-05-06 | 2022-08-19 | 国网上海市电力公司 | Cascading failure searching method and system for power system containing new energy |
CN115859630A (en) * | 2022-12-07 | 2023-03-28 | 南京师范大学 | Electric power traffic coupling network vulnerability assessment method based on probability map |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109118098A (en) * | 2018-08-21 | 2019-01-01 | 山东大学 | The cascading failure methods of risk assessment and system of high proportion wind-electricity integration |
US20210296897A1 (en) * | 2019-11-27 | 2021-09-23 | Robert F. Cruickshank, III | System method and apparatus for providing a load shape signal for power networks |
-
2022
- 2022-03-02 CN CN202210198717.6A patent/CN114493365A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109118098A (en) * | 2018-08-21 | 2019-01-01 | 山东大学 | The cascading failure methods of risk assessment and system of high proportion wind-electricity integration |
US20210296897A1 (en) * | 2019-11-27 | 2021-09-23 | Robert F. Cruickshank, III | System method and apparatus for providing a load shape signal for power networks |
Non-Patent Citations (1)
Title |
---|
MIR HADI ATHARI等: "Impacts of Wind Power Uncertainty on Grid Vulnerability to Cascading Overload Failures", 《IEEE TRANSACTIONS ON SUSTAINABLE ENERGYIND POWER UNCERTAINTY ON GRID VULNERABILITY TO CASCADING OVER》, vol. 9, no. 1, pages 128 - 137 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114925993A (en) * | 2022-05-06 | 2022-08-19 | 国网上海市电力公司 | Cascading failure searching method and system for power system containing new energy |
CN115859630A (en) * | 2022-12-07 | 2023-03-28 | 南京师范大学 | Electric power traffic coupling network vulnerability assessment method based on probability map |
CN115859630B (en) * | 2022-12-07 | 2023-06-16 | 南京师范大学 | Electric traffic coupling network vulnerability assessment method based on probability map |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Orfanos et al. | Transmission expansion planning of systems with increasing wind power integration | |
Aien et al. | Probabilistic optimal power flow in correlated hybrid wind–photovoltaic power systems | |
Ding et al. | Short-term and medium-term reliability evaluation for power systems with high penetration of wind power | |
Hozouri et al. | On the use of pumped storage for wind energy maximization in transmission-constrained power systems | |
Wen et al. | A review on reliability assessment for wind power | |
US8606416B2 (en) | Energy generating system and control thereof | |
CN107679658B (en) | Power transmission network planning method under high-proportion clean energy access | |
Mosadeghy et al. | A time-dependent approach to evaluate capacity value of wind and solar PV generation | |
CN114493365A (en) | Method for evaluating cascading failure vulnerability of power system including wind power plant | |
CN106354985B (en) | Power distribution system reliability assessment method considering distributed power supply | |
CN112288326B (en) | Fault scene set reduction method suitable for toughness evaluation of power transmission system | |
CN110146785B (en) | Method for identifying fragile line of power grid containing wind and solar power supply | |
CN106355308B (en) | A method of wind power integration system core equipment is recognized based on decision tree | |
EP2381094A1 (en) | Energy network and control thereof | |
Wang et al. | Integrated platform to design robust energy internet | |
CN109672215A (en) | Based on load can time shift characteristic distributed photovoltaic dissolve control method | |
McLaughlin et al. | Application of dynamic line rating to defer transmission network reinforcement due to wind generation | |
CN110661250A (en) | Reliability evaluation method and system for wind-solar energy storage and power generation power transmission system | |
Karki et al. | Reliability assessment of a wind integrated hydro-thermal power system | |
Chiodo et al. | Wind farm production estimation under multivariate wind speed distribution | |
McKenna et al. | Impact of wind curtailment and storage on the Irish power system 2020 renewable electricity targets: A free open-source electricity system balancing and market (ESBM) model | |
Li et al. | Composite Power System Reliability Evaluation Considering Space-time Characteristics of Wind Farm | |
Das et al. | A probabilistic approach to assess the adequacy of wind and solar energy | |
Zhang et al. | Adequacy evaluation of wind farm integration in power generation and transmission systems | |
Bagchi et al. | Composite system adequacy assessment incorporating virtual power plants |
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 |