CN107154627B - A kind of small power station's group of planes Dynamic Equivalence suitable for bulk power grid analysis - Google Patents
A kind of small power station's group of planes Dynamic Equivalence suitable for bulk power grid analysis Download PDFInfo
- Publication number
- CN107154627B CN107154627B CN201710446143.9A CN201710446143A CN107154627B CN 107154627 B CN107154627 B CN 107154627B CN 201710446143 A CN201710446143 A CN 201710446143A CN 107154627 B CN107154627 B CN 107154627B
- Authority
- CN
- China
- Prior art keywords
- parameter
- group
- power station
- small power
- planes
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 14
- 238000010206 sensitivity analysis Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 23
- 230000035945 sensitivity Effects 0.000 claims description 19
- 238000004088 simulation Methods 0.000 claims description 9
- 239000002131 composite material Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 230000004069 differentiation Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 3
- 230000004888 barrier function Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 230000001360 synchronised effect Effects 0.000 description 7
- 230000003068 static effect Effects 0.000 description 5
- 230000001052 transient effect Effects 0.000 description 5
- 230000005611 electricity Effects 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000009434 installation Methods 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
- 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]
-
- H02J3/382—
-
- 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
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- 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
-
- 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/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Power Engineering (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of small power station's group of planes Dynamic Equivalences suitable for bulk power grid analysis, the small power station's unit strict equifrom that ensure that in choropleth i.e. small power station's group of planes is differentiated by same tone, key parameter and non-key parameter are distinguished using sensitivity analysis, and make full use of the system features for including in more PMU multiple faults information, dynamic equivalent model is improved to the interpretability of true small power station's group of planes, it treats identified parameters and carries out otherness identification, effectively eliminate dynamic equivalent model there are problems that solving more, to effectively improve the accuracy and robust ability of small power station's group of planes dynamic equivalent model.
Description
Technical field
The invention belongs to electric system simulation modeling technique fields, more specifically, are related to a kind of suitable for bulk power grid
Small power station's group of planes Dynamic Equivalence of analysis.
Background technique
In recent years, renewable including hydroelectric resources with the continuous aggravation of global energy crisis and environmental pollution
The energy has obtained quick development.Cut-off 2015, hydropower installed capacity accounts for 16.6%, Zhan Ke of world's total electricity installed capacity
The 70% of renewable sources of energy installed capacity, it is contemplated that future 25 years will be increased with annual 3.1% ratio.
Due to being influenced by geographical environment, in selected areas of China, there are a large amount of small power station's units, (single station capacity is less than
The phenomenon grid-connected with net side (35kV) 25MW) is concentrated in a certain region.A large amount of small power station's group of planes are flanked into power grid, no in distribution
Power distribution network can only be impacted, while the stability of bulk power grid can also be impacted.During bulk power grid modeling analysis,
Due to the features such as small power station's group of planes is large number of, and structure is complicated, and parameter is opaque, it is difficult to be built in detail in simulation model
Mould.Meanwhile bulk power grid analysis is carried out using detailed model, cumulative errors can be not only introduced, but also easily cause " dimension calamity ".
Dynamic equivalent (PMU based Dynamic Equivalent, PDE) based on PMU metrical information can be very good
Solve this problem.In the case where guaranteeing choropleth dynamic characteristic unanimous circumstances, PDE carries out choropleth using reduced-order model
Simplify modeling, the cumulative errors in detailed modeling process are effectively eliminated while improving modeling efficiency, avoid " dimension calamity ".So
And find in practical applications, although the small power station's group of planes Equivalent Model obtained based on traditional PDE method can be very good to reappear
It is faulty, but be difficult to the unknown failure of Accurate Prediction, that is, there is a problem of that Equivalent Model robust ability is not strong.By the people having the same aspiration and interest
Property it is theoretical it is found that only when the completely same timing of unit in choropleth dynamic equivalent operation could be carried out to it；Simultaneously as
The strong nonlinearity feature of electric system, PDE of the tradition based on single failure information are difficult to reflect the dynamic characteristic of system completely；Most
Afterwards, in traditional PDE method, since parameter to be identified is excessive, Equivalent Model is easy to appear the problem of solving more.These are all caused
The not strong major reason of small power station's group of planes Equivalent Model robust ability.Therefore, continue to increase renewable energy access electricity in various countries
Under the overall background of net ratio, establishes and be suitable for bulk power grid analysis, and the small power station with high accuracy and stronger robust ability
Group of planes dynamic equivalent model has important practical significance.
Summary of the invention
It is an object of the invention to overcome the deficiencies in the prior art, propose a kind of small power station's machine suitable for bulk power grid analysis
Group's Dynamic Equivalence, to effectively improve the accuracy and robust ability of small power station's group of planes dynamic equivalent model.
For achieving the above object, the present invention is suitable for small power station's group of planes Dynamic Equivalence of bulk power grid analysis,
It is characterized in that, comprising the following steps:
(1), using time-domain-simulation method of discrimination, same tone differentiation is carried out to needing to carry out equivalent small power station's group of planes, such as
The not exclusively same timing of small power station's unit in a fruit small power station group of planes, then according to same tone criterion to the small power station in small power station's group of planes
Unit carries out a point group, until the unit in each point of group meets same tone requirement；
(2), to the small power station's group of planes for meeting same tone requirement, small power station's group of planes dynamic equivalent model is established；
(3), sensitivity analysis is carried out to dynamic equivalent model parameter using relevant path sensitivity, selected to small power station's machine
For the parameter to be identified that group motion step response is affected as key parameter, remaining parameter to be identified is non-key parameter；
(4), it is based on more PMU multiple faults information, applied statistics theory and comentropy method, extraction can more reflect small power station
More PMU multiple faults features of group of planes dynamic characteristic；
(5), otherness identification is carried out to key parameter and non-key parameter based on more PMU multiple faults features, establishes dynamic
Equivalent Model.
The object of the present invention is achieved like this.
The present invention is suitable for small power station's group of planes Dynamic Equivalence of bulk power grid analysis, ensure that by same tone differentiation
The small power station's unit strict equifrom being worth in region, that is, small power station's group of planes distinguishes key parameter and non-key ginseng using sensitivity analysis
Number, and the system features for including in more PMU multiple faults information are made full use of, dynamic equivalent model is improved to true small power station's group of planes
Interpretability, treat identified parameters and carry out otherness identification, effectively eliminating dynamic equivalent model has more solutions, thus
Effectively improve the accuracy and robust ability of small power station's group of planes dynamic equivalent model.
Detailed description of the invention
Fig. 1 is the flow chart suitable for small power station's group of planes Dynamic Equivalence of bulk power grid analysis；
Fig. 2 is the true small power station's group of planes wiring diagram in Sichuan Province；
Fig. 3 is the power-angle curve in the case of a certain typical fault of small power station's group of planes shown in Fig. 2；
Fig. 4 is the active power curves at equivalent bus；
Fig. 5 is the active power curves at external system bus 20.
Specific embodiment
A specific embodiment of the invention is described with reference to the accompanying drawing, preferably so as to those skilled in the art
Understand the present invention.Requiring particular attention is that in the following description, when known function and the detailed description of design perhaps
When can desalinate main contents of the invention, these descriptions will be ignored herein.
Fig. 1 is the flow chart suitable for small power station's group of planes Dynamic Equivalence of bulk power grid analysis.
In the present embodiment, as shown in Figure 1, the present invention is suitable for small power station's group of planes Dynamic Equivalence of bulk power grid analysis
The following steps are included:
Step S1: same tone differentiates
Using time-domain-simulation method of discrimination, same tone differentiation is carried out to small power station's group of planes, if small in small power station's group of planes
The not exclusively same timing of Hydropower Unit then carries out a point group to small power station's unit in small power station's group of planes according to same tone criterion, until
Unit in each point of group meets same tone requirement.
In the present embodiment, described to use time-domain-simulation method of discrimination, it is carried out to needing to carry out equivalent small power station's group of planes
Same tone differentiates are as follows:
1.1), under a certain fault condition, record is differentiated each small power station's unit power-angle curve；
If 1.2), after the failure occurred in [0, τ] time range, any Liang Tai small power station unit at any time opposite
Generator rotor angle difference is not more than given value ε, then differentiates small power station's group of planes people having the same aspiration and interest, it may be assumed that
Wherein, Δ σmn(t) poor in the generator rotor angle of t moment for small power station's unit m and n, usual ε is 5~10 ° desirable, and τ desirable 1~
3s；
1.3) aforesaid operations, that is, step 1.1), 1.2), is repeated in different typical faults, if in all typical faults
In the case of, small power station's group of planes is all satisfied same tone requirement, then all motor people having the same aspiration and interest in small power station's group of planes, only when small power station's group of planes
Meet same tone criterion, dynamic equivalent operation could be carried out to it.
Step S2: small power station's group of planes dynamic equivalent model is established
To the small power station's group of planes for meeting same tone requirement, small power station's group of planes dynamic equivalent model is established.In the present embodiment,
The small power station's group of planes for meeting same tone is reduced to one 5 rank synchronous generators and a static state ZIP load parallel model, in which: 5 ranks
Synchronous generator mathematic(al) representation is as follows:
D axis electric parameter:
Q axis electric parameter:
Equation of rotor motion:
Wherein, Xd,XqFor d, q axis synchronous reactance；Xd' it is the transient state reactance of d axis；Xd”,Xq" it is d, q axis subtranient reactance；
Td0' it is d axis transient state open circuit time constant；Td0”,Tq0" it is time transient state open circuit time constant；Ed',Eq' it is d, q axis transient potential；
Ed”,Eq" it is d, q axis time transient potential；EfFor excitation potential；ud,uqFor d, q axis stator voltage；id,iqFor d, q axis stator electricity
Stream；
Wherein, static ZIP load model mathematic(al) representation is as follows:
Wherein, P, Q are respectively the active and reactive power of load consumption when busbar voltage is U；U0For voltage rating；P0,
Q0Active and reactive power consumed by load when for voltage rating；PZ,PI,PSAnd QZ,QI,QSThe respectively active and reactive function of load
The proportionality coefficient of constant-impedance, constant current, invariable power component in rate, and meet PZ+PI+PS=QZ+QI+QS=1.
Step S3: key parameter and non-key parameter are selected
Sensitivity analysis is carried out to dynamic equivalent model parameter using relevant path sensitivity, is selected dynamic to small power station's group of planes
For the parameter to be identified that step response is affected as key parameter, remaining parameter to be identified is non-key parameter.In this implementation
In example, specifically include:
3.1), defining the bus that small power station's group of planes dynamic equivalent model is connect with bulk power grid is system equivalent bus；
3.2), the active and reactive power at selecting system equivalence bus seeks parameter to be identified as observation track variable
θhActive and reactive relevant path sensitivity TSPh、TSQh:
Wherein,For active power track yPIn k-th of sampled point for parameter θ to be identifiedhPartial derivative, yP0
For active power track yPSteady-state value,For reactive power track yQIn k-th of sampled point for parameter θ to be identifiedh
Partial derivative, yQ0For reactive power track yQSteady-state value, θh0For parameter θ to be identifiedhStandard value, K is sampling number, and h is
The serial number of parameter to be identified, h=1,2 ..., H, H are the quantity of parameter to be identified；
By parameter θ to be identifiedhActive and reactive relevant path sensitivity TSPh、TSQhIt is added, obtains parameter θ to be identifiedh
Composite trochoid sensitivity；Key parameter selection threshold value is set, if composite trochoid sensitivity is greater than given threshold wait distinguish
Know parameter and be set to key parameter, remaining parameter is set to non-key parameter.
Step S4: more PMU multiple faults features are extracted
Based on more PMU multiple faults information, applied statistics theory and comentropy method, extraction can more reflect small power station's group of planes
More PMU multiple faults features of dynamic characteristic.In the present embodiment, comprising the following steps:
4.1), more PMU multiple faults information are chosen
Selected more PMU multiple faults information refer to as the error assessment in dynamic equivalent identification of Model Parameters objective function
It marks (Error Evaluation Indicator, EEI), should be able to more comprehensively reflect the dynamic characteristic of bulk power grid.In this reality
It applies in example, the important bus active and reactive power of system under the conditions of selection typical fault including equivalent bus is as small water
The EEI of motor group dynamic equivalent；
4.2), applied statistics theory and comentropy method extract more PMU multiple faults features
More PMU multiple faults feature extractions are weight of the selected EEI in Equivalent Model parameter identification objective function
Determination process:
4.2.1), according to small power station's group of planes gross rated capacity and the gross rated capacity of load, it is to be identified to provide Equivalent Model
The theoretical value of parameter；
In certain confidence interval (the 200% of the 10% of theoretical value to theoretical value) range, it is to be identified that multiple groups are randomly generated
Parameter combination, different faults condition reappeared in Equivalent Model PMU record historical failure event, accordingly, can obtain it is different to
Identified parameters combination under the conditions of different faults with the error matrix E of PMU measured data are as follows:
Wherein, i is failure number, and i=1,2 ..., I, I is number of faults, and j is under same fault condition by PMU record
EEI number, j=1,2 ..., J, J are EEI number, and l is the parameter group # to be identified generated at random, L parameter group number to be identified,For in l sample, error (Equivalent Model simulation value and the PMU measured data root mean square of j-th of EEI under i fault condition
Error),For j-th of EEI in first of sample under i fault condition equivalent model emulation value and PMU measured data in t
The difference of (t sampled point) is carved, T is sampling number；
In error matrix E, each column represent an EEI, and every a line represents used in primary parameter identification process
Whole EEI；
4.2.2), each EEI is subjected to negative sense criterion:
Wherein, l=1 ..., L；
4.2.3), the specific gravity of i-th of failure, j-th of EEI in sample l is calculated:
4.2.4), i-th of failure, the comentropy of j-th EEI are calculated:
Wherein, g=1/ln (L) is a constant related with sample number L, to make comentropy Entm∈[0,1]；
4.2.5), the redundancy of comentropy is calculated:
dij=1-Entij (12)
4.2.6), i-th of failure, the weight of j-th EEI are calculated:
Accordingly, it is based on Principle of Statistics and information entropy theory, obtains each EEI in Equivalent Model parameter identification objective function
Weighted value, by eijMultiplied by respective weights omegaijAnd be added, it can be obtained the more PMU that can more reflect small power station's group of planes dynamic characteristic
Multiple faults feature C:
Step S5: otherness identification is carried out to key parameter and non-key parameter based on more PMU multiple faults features, is established dynamic
State Equivalent Model.In the present embodiment, it specifically includes:
5.1), preliminary identification
Successively all parameters to be identified are recognized using more PMU information under the conditions of single failure, at this point, parameter is distinguished
The objective function of knowledge are as follows:
Wherein, MSEMPSF_iFor based on the comprehensive mean square error of EEI under the conditions of more PMU single faults；
5.2), emphasis recognizes
According to preliminary identification result, non-key parameter is assigned to average identifier when preliminary identification, while using more PMU
Multiple faults information carries out emphasis identification to key parameter；At this point, the objective function of parameter identification are as follows:
Wherein, MSEMPMFFor based on the comprehensive mean square error of EEI under the conditions of more PMU multiple faults.
In addition, dividing group to be reduced to one each small power station if small power station's group of planes is divided into n groups according to same tone criterion
Group dynamic equivalent model；
Successively group is divided to carry out dynamic equivalent operation each small power station according to step S2 to step S5；When to 1st point
When group carries out dynamic equivalent, remaining divides group using detailed model, when carrying out dynamic equivalent to i-th (2≤i≤n-1) a point of group,
Divide group to use Equivalent Model before i-th point of group, divides group using detailed model after i-th point of group, when to n-th point of group
When carrying out dynamic equivalent, remaining divides group using Equivalent Model.
Compared with prior art, the present invention having the advantage that
1, according to people having the same aspiration and interest Theory of Equivalence, only when the same timing of generating set, dynamic equivalent operation could be carried out to it.This hair
It is bright that same tone verifying is carried out to small power station's unit in choropleth before carrying out dynamic equivalent to small power station's group of planes, guarantee every
Small power station's group of planes strict equifrom of inferior Value Operations is eliminated not influence of the people having the same aspiration and interest unit to Equivalent Model dynamic characteristic, is effectively mentioned
High dynamic equivalence precision.
2, parameter to be identified is divided by key parameter and non-key parameter using relevant path sensitivity, and difference is carried out to it
Opposite sex identification.Effectively reduce the dimension of Equivalent Model solution space, eliminate parameter identification solve more problem to equivalent model accuracy and
The influence of robust ability.
3, the more PMU multiple faults data that can reflect system dynamic characteristic comprehensively are introduced during parameter identification, effectively
Equivalent Model is improved to the interpretability of real system, further increases the accuracy and robust ability of Equivalent Model.
Example
Fig. 2 show the true small power station's group of planes wiring diagram in Sichuan Province.9 Hydropower Unit G are shared in small power station's group of planes1~
G9, it is directly accessed 35kV distribution network.Total installation of generating capacity 75MVA in choropleth, total load 14MW+36MVAr, inner wire
Lu overall length 185km, small power station's group of planes with certain 220kV bus at access external bulk power grid.
1, using time-domain-simulation method of discrimination, equivalent, shown in Fig. 2 small power station's group of planes is carried out to needs and carries out same tone
Differentiate.
Under the conditions of all typical faults of external system, small power station's group of planes is all satisfied same tone requirement.A certain typical fault
In the case of, small power station's group of planes power-angle curve is as shown in Figure 3, wherein G9Generator rotor angle is with reference to 0 °.Generator rotor angle difference is up to 8.2 degree, meets
It is required that.
2, to the small power station's group of planes for meeting same tone requirement, small power station's group of planes dynamic equivalent model is established.
Small power station's group of planes is reduced to one 5 rank synchronous generators and a static state ZIP load parallel model.In the present embodiment
In, synchronous generator parameter to be identified are as follows: Xd, Xd', Xd", Xq, Xq", Td0', Td0", Tq0", H, D.ZIP load model is to be identified
Parameter is P0, Q0, PZ, PI, PS, QZ, QI, QS。
3, sensitivity analysis is carried out to dynamic equivalent model parameter using relevant path sensitivity, selected to small power station's group of planes
For the parameter to be identified that dynamic characteristic is affected as key parameter, remaining parameter to be identified is non-key parameter.
In this example, the active and reactive power at selecting system equivalence bus is sought as observation track variable wait distinguish
Know parameter θhActive and reactive relevant path sensitivity TSPh、TSQh, by parameter θ to be identifiedhActive and reactive relevant path spirit
Sensitivity TSPh、TSQhIt is added, obtains parameter θ to be identifiedhComposite trochoid sensitivity.
In this example, the relevant path sensitivity of parameter to be identified is as shown in table 1:
Table 1
To guarantee representativeness of the key parameter in parameter to be identified, 5 rank synchronous generators are respectively set and static state ZIP is negative
Lotus model key parameter selects threshold value.The 5 rank synchronous generator key parameter composite trochoid thresholds of sensitivity are set as 0.04, static state
The ZIP load model key parameter composite trochoid threshold of sensitivity is set as 0.4.Therefore, key parameter in parameter to be identified are as follows:
Xd, Xd', Xd", Xq, Xq", H, P0, Q0, QZ, QI, QS.Non-key parameter are as follows: Td0', Td0", Tq0", D, PZ, PI, PS。
4, it is based on more PMU multiple faults information, applied statistics theory and comentropy method, extraction can more reflect system dynamic
More PMU multiple faults features of characteristic:
4.1, as shown in Fig. 2, the active and reactive power at the equivalent bus of selection, bus 5 and bus 17 is as Equivalent Model
EEI in parameter identification objective function.In addition, the typical fault of three PMU record is also selected for EEI weight calculation, point
It is not failure 1: 1 singlephase earth fault of bus.Failure 2: 11 double earthfault of bus.Failure 3: 23 three-phase ground of bus event
Barrier.
4.2, according to the gross rated capacity of small power station's unit gross rated capacity and load in choropleth, Equivalent Model is provided
The theoretical value of parameter to be identified.In certain confidence interval (the 200% of the 10% of theoretical value to theoretical value) range, it is randomly generated
Different faults condition is arranged in multiple groups parameter combination to be identified, and the historical failure event of PMU record is reappeared in Equivalent Model.According to
This, can obtain different parameter combinations to be identified under the conditions of different faults with the error matrix E of PMU measured data.
Each EEI weight is calculated under the conditions of failure 1, failure 2 and failure 3 respectively, the results are shown in Table 2:
| EEIP1 | EEIP2 | EEIP3 | EEIQ1 | EEIQ2 | EEIQ3 | |
| Failure 1 | 0.1912 | 0.0855 | 0.0856 | 0.3712 | 0.0834 | 0.1831 |
| Failure 2 | 0.1777 | 0.0703 | 0.0704 | 0.3860 | 0.0896 | 0.2060 |
| Failure 3 | 0.1902 | 0.0861 | 0.0862 | 0.3717 | 0.0835 | 0.1822 |
Table 2
Wherein, EEIPiAnd EEIQiEEI active power, reactive power weight at respectively selected bus i.
Each EEI weight is calculated under three fault conditions, the results are shown in Table 3:
| EEIP11 | EEIP21 | EEIP31 | EEIQ11 | EEIQ21 | EEIQ31 |
| 0.0812 | 0.0367 | 0.0368 | 0.1580 | 0.0356 | 0.0777 |
| EEIP12 | EEIP22 | EEIP32 | EEIQ12 | EEIQ22 | EEIQ32 |
| 0.0237 | 0.0122 | 0.0122 | 0.0573 | 0.0140 | 0.0310 |
| EEIP12 | EEIP22 | EEIP32 | EEIQ12 | EEIQ22 | EEIQ32 |
| 0.0806 | 0.0366 | 0.0366 | 0.1572 | 0.0354 | 0.0772 |
Table 3
Wherein, EEIPijAnd EEIQijActive power, reactive power of the EEI under the conditions of failure j at respectively selected bus i
Weight.
5, otherness identification is carried out to key parameter and non-key parameter based on more PMU multiple faults features, establishes equivalent essence
Degree is high and has the dynamic equivalent model of stronger robust ability.
5.1, preliminary identification.
Respectively using the active and reactive function of equivalent bus, bus 5 and bus 17 under the conditions of failure 1, failure 2 and failure 3
Rate recognizes all parameters to be identified.Under the conditions of each single fault, EEI weight setting such as table 2 in parameter identification objective function
It is shown.
5.2, emphasis recognizes.According to preliminary identification result, non-key parameter is assigned to average identifier when preliminary identification,
Emphasis identification is carried out to key parameter using three fault messages simultaneously.Under the conditions of multiple faults, EEI is weighed in parameter identification objective function
It resets and sets as shown in table 3.
To illustrate accuracy and robust ability using present invention gained Equivalent Model, it is respectively adopted equivalent obtained by the present invention
(under the conditions of certain single failure (failure 1), dynamic equivalent objective function is active at equivalent bus for model and traditional equivalence method
The sum of idle root-mean-square error) gained Equivalent Model for playback system break down 4 (10 double earthfault of bus) when
Dynamic characteristic.Active power curves at equivalent bus are as shown in figure 4, the active power curves at external system bus 20 are such as schemed
Shown in 5.By Fig. 4 and Fig. 5 it is found that after using the mentioned method of the present invention, small power station's group of planes dynamic equivalent model accuracy and robust
Ability significantly improves.
Although the illustrative specific embodiment of the present invention is described above, in order to the technology of the art
Personnel understand the present invention, it should be apparent that the present invention is not limited to the range of specific embodiment, to the common skill of the art
For art personnel, if various change the attached claims limit and determine the spirit and scope of the present invention in, these
Variation is it will be apparent that all utilize the innovation and creation of present inventive concept in the column of protection.
Claims (3)
1. a kind of small power station's group of planes Dynamic Equivalence suitable for bulk power grid analysis, which comprises the following steps:
(1), using time-domain-simulation method of discrimination, same tone differentiation is carried out to needing to carry out equivalent small power station's group of planes, if small
The not exclusively same timing of small power station's unit in a water power group of planes, then according to same tone criterion to small power station's unit in small power station's group of planes
A point group is carried out, until the unit in each point of group meets same tone requirement；
(2), to the small power station's group of planes for meeting same tone requirement, small power station's group of planes dynamic equivalent model is established；
(3), sensitivity analysis is carried out to dynamic equivalent model parameter using relevant path sensitivity, selected dynamic to small power station's group of planes
The influential parameter to be identified of step response is non-key parameter as key parameter, remaining parameter to be identified:
3.1), defining the bus that small power station's group of planes dynamic equivalent model is connect with bulk power grid is system equivalent bus；
3.2), the active and reactive power at selecting system equivalence bus seeks parameter θ to be identified as observation track variableh's
Active and reactive relevant path sensitivity TSPh、TSQh:
Wherein,For active power track yPIn k-th of sampled point for parameter θ to be identifiedhPartial derivative, yP0To have
Function power track yPSteady-state value,For reactive power track yQIn k-th of sampled point for parameter θ to be identifiedhIt is inclined
Derivative, yQ0For reactive power track yQSteady-state value, θh0For parameter θ to be identifiedhStandard value, K is sampling number, and h is wait distinguish
Know the serial number of parameter, h=1,2 ..., H, H is the quantity of parameter to be identified；
By parameter θ to be identifiedhActive and reactive relevant path sensitivity TSPh、TSQhIt is added, obtains parameter θ to be identifiedhSynthesis
Trace sensitivity；Key parameter selection threshold value is set, if composite trochoid sensitivity is greater than the parameter to be identified of given threshold
It is set to key parameter, remaining parameter is set to non-key parameter；
(4), it is based on more PMU multiple faults information, applied statistics theory and comentropy method, extraction can more reflect small power station's group of planes
More PMU multiple faults features of dynamic characteristic:
4.1), more PMU multiple faults information are chosen
The important bus active and reactive power of system under the conditions of selection typical fault including equivalent bus is as small power station's machine
Error assessment index, that is, EEI of group's dynamic equivalent；
4.2), applied statistics theory and comentropy method extract more PMU multiple faults features
More PMU multiple faults feature extractions are that weight of the selected EEI in Equivalent Model parameter identification objective function determines
Process:
4.2.1), according to small power station's group of planes gross rated capacity and the gross rated capacity of load, Equivalent Model parameter to be identified is provided
Theoretical value；
In 10% 200% range to theoretical value of theoretical value, multiple groups parameter combination to be identified is randomly generated, in different faults
Condition reappears the historical failure event of PMU record in Equivalent Model, accordingly, can obtain different parameter combinations to be identified in different events
Under the conditions of barrier with the error matrix E of PMU measured data are as follows:
Wherein, i is failure number, and i=1,2 ..., I, I is number of faults, and j is to be compiled under same fault condition by the EEI of PMU record
Number, j=1,2 ..., J, J are EEI number, and l is the parameter group # to be identified generated at random, L parameter group number to be identified,For
In l sample, (Equivalent Model simulation value and PMU measured data root mean square miss the error of j-th of EEI under i fault condition
Difference),For j-th of EEI in first of sample under i fault condition equivalent model emulation value and PMU measured data in t moment
The difference of (t sampled point), T are sampling number；
In error matrix E, each column represent an EEI, and every a line represents used complete in primary parameter identification process
Portion EEI；
4.2.2), each EEI is subjected to negative sense criterion:
Wherein, l=1 ..., L；
4.2.3), the specific gravity of i-th of failure, j-th of EEI in sample l is calculated:
4.2.4), i-th of failure, the comentropy of j-th EEI are calculated:
Wherein, g=1/ln (L) is a constant related with sample number L, to make comentropy Entm∈[0,1]；
4.2.5), the redundancy of comentropy is calculated:
dij=1-Entij
4.2.6), i-th of failure, the weight of j-th EEI are calculated:
By eijMultiplied by respective weights omegaijAnd be added, it can be obtained the more PMU multiple faults features that can more reflect bulk power grid dynamic characteristic
C:
(5), otherness identification is carried out to key parameter and non-key parameter based on more PMU multiple faults features, establishes dynamic equivalent
Model.
2. small power station's group of planes Dynamic Equivalence according to claim 1, which is characterized in that in step (1), the use
Time-domain-simulation method of discrimination is judged to the small power station's group of planes progress same tone for needing to carry out equivalent:
1.1), under a certain fault condition, record is differentiated each small power station's unit power-angle curve；
If 1.2), after the failure occurred in [0, τ] time range, the opposite generator rotor angle of any Liang Tai small power station unit at any time
Difference is not more than given value ε, then differentiates small power station's group of planes people having the same aspiration and interest, it may be assumed that
Wherein, Δ σmn(t) poor in the generator rotor angle of t moment for small power station's unit m and n, ε takes range are as follows: 5≤ε≤10 °, τ take range
Are as follows: 1≤τ≤3s.
3. small power station's group of planes Dynamic Equivalence according to claim 1, which is characterized in that described to be based in step (5)
More PMU multiple faults features carry out otherness identification to key parameter and non-key parameter, establish dynamic equivalent model are as follows:
5.1), preliminary identification
Successively all parameters to be identified are recognized using more PMU information under the conditions of single failure, at this point, parameter identification
Objective function are as follows:
Wherein, MSEMPSF_iFor based on the comprehensive mean square error of EEI under the conditions of more PMU single faults；
5.2), emphasis recognizes
According to preliminary identification result, non-key parameter is assigned to average identifier when preliminary identification, while using more PMU mostly event
Hinder information and emphasis identification is carried out to key parameter；At this point, the objective function of parameter identification are as follows:
Wherein, MSEMPMFFor based on the comprehensive mean square error of EEI under the conditions of more PMU multiple faults.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710446143.9A CN107154627B (en) | 2017-06-14 | 2017-06-14 | A kind of small power station's group of planes Dynamic Equivalence suitable for bulk power grid analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710446143.9A CN107154627B (en) | 2017-06-14 | 2017-06-14 | A kind of small power station's group of planes Dynamic Equivalence suitable for bulk power grid analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107154627A CN107154627A (en) | 2017-09-12 |
| CN107154627B true CN107154627B (en) | 2019-09-24 |
Family
ID=59795260
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710446143.9A Active CN107154627B (en) | 2017-06-14 | 2017-06-14 | A kind of small power station's group of planes Dynamic Equivalence suitable for bulk power grid analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107154627B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107658865B (en) * | 2017-09-21 | 2019-08-20 | 电子科技大学 | A kind of multimachine PSS parameter coordination optimizing method under multi-operating condition |
| CN107832264A (en) * | 2017-11-20 | 2018-03-23 | 宁德师范学院 | The overall evaluation system and evaluation method of a kind of design of electrical motor |
| CN117559518B (en) * | 2023-10-16 | 2024-06-21 | 国家电网有限公司华东分部 | Equivalent model construction method and device for wind power generation system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104090985A (en) * | 2014-07-25 | 2014-10-08 | 武汉大学 | Active disconnection optimum fracture surface searching method based on electrical distance |
| CN106250640A (en) * | 2016-08-04 | 2016-12-21 | 山东大学 | A kind of layering Dynamic Equivalence being applicable to area power grid |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10097240B2 (en) * | 2013-02-19 | 2018-10-09 | Astrolink International, Llc | System and method for inferring schematic and topological properties of an electrical distribution grid |
-
2017
- 2017-06-14 CN CN201710446143.9A patent/CN107154627B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104090985A (en) * | 2014-07-25 | 2014-10-08 | 武汉大学 | Active disconnection optimum fracture surface searching method based on electrical distance |
| CN106250640A (en) * | 2016-08-04 | 2016-12-21 | 山东大学 | A kind of layering Dynamic Equivalence being applicable to area power grid |
Non-Patent Citations (3)
| Title |
|---|
| Comparison of Ensemble Decision Tree Methods for On-line Identification of Power System Dynamic Signature Considering Availability of PMU Measurements;Tingyan Guo，P. Papadopoulos, P. Mohammed, and J. V. Milanović;《2015 IEEE Eindhoven PowerTech》;20150903;1-6 * |
| 基于系统同调性的P M U 最优布点;许剑冰,薛禹胜;《电力系统自动化》;20041010;第28卷(第19期);22-26 * |
| 小水电机群对特高压电网功率振荡特性的影响研究;李阳海，刘 巨，文劲宇;《高电压技术》;20150331;第41卷(第3期);762-769 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107154627A (en) | 2017-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102377180B (en) | Power system load modeling method based on electric energy quality monitoring system | |
| CN107154627B (en) | A kind of small power station's group of planes Dynamic Equivalence suitable for bulk power grid analysis | |
| CN102539822B (en) | Short-term wind speed forecasting method | |
| Wu et al. | Review on transient stability prediction methods based on real time wide-area phasor measurements | |
| CN103279639A (en) | Receiving-end network voltage stabilization overall process situation assessment and prevention and control method based on responses | |
| CN111679158A (en) | Power distribution network fault identification method based on synchronous measurement data similarity | |
| CN108667069A (en) | A kind of short-term wind power forecast method returned based on Partial Least Squares | |
| CN110348114B (en) | Non-precise fault identification method for power grid completeness state information reconstruction | |
| CN107145707A (en) | It is a kind of to count and photovoltaic is exerted oneself the power distribution network transformer planing method of uncertain and overall life cycle cost | |
| CN103995948B (en) | A kind of voltage of oscillation center Forecasting Methodology based on multinomial model | |
| CN104362622A (en) | Load model of regional power grid with small hydroelectric generating set and modeling method | |
| CN108074198B (en) | Power system PMU bad data identification method | |
| CN111262238B (en) | Machine learning-based method for predicting short-circuit current of power distribution network containing IIDG | |
| CN108074048A (en) | It is included in the wind-electricity integration power system security methods of risk assessment of wind speed correlation properties | |
| CN104502800A (en) | Extraction method of fault factor characteristics of electrical power system | |
| He et al. | Power system frequency situation prediction method based on transfer learning | |
| CN102522759A (en) | Method for distinguishing transient state instability of electric power system in real time based on voltage track after disturbance | |
| CN105701265A (en) | Double-fed wind generator modeling method and apparatus | |
| Li et al. | Augmented state estimation method for fault location based on on-line parameter identification of PMU measurement data | |
| CN113093052B (en) | Method for judging turn number of stator winding fault of turn-to-turn short circuit motor | |
| CN105896525B (en) | Power grid large disturbances recognition methods based on WAMS | |
| CN106410836A (en) | Evaluation method of transient disturbance interactive influence states of MMC (multilevel modular converter)-HVDC (high voltage direct current transmission) | |
| Zhao et al. | A novel real-time transient stability prediction method based on post-disturbance voltage trajectories | |
| Li et al. | Research concept of pattern features of generalized load curves in future power systems | |
| CN109978715A (en) | User side distributed generation resource Data Reduction method and device |
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 |