CN107134811A - Network load regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution - Google Patents
Network load regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution Download PDFInfo
- Publication number
- CN107134811A CN107134811A CN201710540775.1A CN201710540775A CN107134811A CN 107134811 A CN107134811 A CN 107134811A CN 201710540775 A CN201710540775 A CN 201710540775A CN 107134811 A CN107134811 A CN 107134811A
- Authority
- CN
- China
- Prior art keywords
- frequency
- load
- spare capacity
- value
- annual
- 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 80
- 230000033228 biological regulation Effects 0.000 title claims abstract description 54
- 230000004044 response Effects 0.000 claims description 16
- 230000001186 cumulative effect Effects 0.000 claims description 15
- 230000005619 thermoelectricity Effects 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000011068 loading method Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000000205 computational method Methods 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 239000000470 constituent Substances 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 241001269238 Data Species 0.000 claims description 2
- 238000009434 installation Methods 0.000 claims description 2
- 238000009825 accumulation Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 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
- 238000004458 analytical method Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method 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
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- 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
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Power Engineering (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
What regulation spare capacity of the present invention based on system was tackled is the power deviation that load fluctuation is caused, and proposes the regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution, assesses the reserve capacity for load variation in power size of balanced load fluctuation.The present invention comprises the following steps：The active power deviation that unit output and load fluctuation are caused causes system frequency excursion rated value, the structure regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution；Load fluctuation generation system frequency offseting value, according to power system resultant frequency characteristics, tries to achieve system active power deviation.
Description
Technical field
The present invention relates to reserve capacity of power grid field, more particularly to a kind of bulk power grid based on frequency shift (FS) probability distribution
Load Regulation spare capacity appraisal procedure.
Background technology
With the traditional calculated load based on generation load percentage of the improvement of sustainable development and power technology of power network
The need for spare capacity method does not adapt to Power System Analysis, while power system reform deepens continuously and extensive new
The access of the energy increases the uncertain fixed of Real-time Load, selection of the load fluctuation caused to power system load spare capacity
The influence of generation is increasing.
In actual motion, power system load is fluctuated real-time, and the change that generating set is exerted oneself needs
Certain time, generating set can not balanced load in real time, this allows for occurring active power between load and unit output
Deviation, causes system frequency excursion rated value 50Hz.Therefore, system frequency excursion reflects system active power deviation, and has
Work(power deviation reflects the size of load fluctuation again.The present invention proposes a kind of new Load Regulation based on system frequency excursion
Spare capacity appraisal procedure, has reference value to improving Power System Planning, operation standard.
The content of the invention
What the regulation spare capacity of system was tackled is the power deviation that load fluctuation is caused, from this angle, this hair
Bright regulation spare capacity appraisal procedure of the proposition based on frequency shift (FS) probability distribution, assesses the standby appearance of load of balanced load fluctuation
Measure size.
To achieve the above object, the present invention is adopted the following technical scheme that：A kind of big electricity based on frequency shift (FS) probability distribution
Net Load Regulation spare capacity appraisal procedure, comprises the following steps：
Step S1：The active power deviation of unit output and load causes system frequency excursion rated value, inclined based on frequency
Move the structure regulation spare capacity appraisal procedure of probability distribution；
Step S2：Load fluctuation generation system frequency offseting value Δ f, can be in the hope of according to power system resultant frequency characteristics
Obtaining system active power deviation is：
Δ P=β × Δ f (1), as adjust spare capacity,
Wherein：β is power system frequency response characteristic coefficient, and unit is MW/0.1Hz.
Further, in the step S1, power system is once adjusted, causes system frequency excursion, specific system change
Process is as follows：
When a certain moment A system loadings increase suddenly, cause generator output to be less than load, negative active power occur inclined
Poor Δ P, at this moment frequency decline rapidly, power system, which is once adjusted, to play a role rapidly, and X has arrived maximum at the time of after certain time
Frequency offset, a certain moment B after an event occurs, once adjustment terminates, and then proceeds by the secondary of frequency from moment B
Adjustment, AGC control unit increases are exerted oneself, to Time To Event after at the time of C, system frequency starts slowly to 50Hz to recover.
Further, the system frequency excursion value Δ f algorithms in the step S2 are as follows：
Power system is per second after being filtered in real time to system frequency interference signal by SCADA system in actual motion
Clock is once sampled, to meet nyquist sampling theorem, as sample frequency fs.maxMore than or equal to highest frequency f in signalmax
2 times when (fs.max≥2fmax), a frequency values are obtained using a stepped-frequency signal per minute, remembered using a large amount of frequency shift (FS)s
Data are recorded, can be in the hope of active power deviation delta P probability distribution situation can be obtained.
Further, frequency response characteristic factor beta is time-varying and nonlinear, and the β algorithms in the step S2 are as follows：
Method 1, FR coefficient is taken as the certain percentage of annual peak load：
In existing achievement in research, it is maximum that FR characteristic coefficient is all taken as the year by most of power system
The certain percentage of load, i.e. peak load percentage, the FR coefficient of most of control zones are estimated negative in annual highest
Between the 1%~1.5% of lotus；
Method 2, FR coefficient is taken as the certain percentage of actual load：
Consider the time-varying characteristics of FR factor beta, certain percentage is multiplied by using actual load and obtained；
Method 3, three-stage FR Y-factor method Y：
Different primary frequency modulation dead bands are generally set in view of thermoelectricity, Hydropower Unit in power system, with reference to pertinent literature,
Beta coefficient can be configured for three sections according to fired power generating unit, Hydropower Unit Regulation dead-band as boundary point, thermoelectricity, Hydropower Unit
The dead band of primary frequency modulation is respectively ± fFire, ± fWater, three-stage FR factor beta establishing method is as follows：
1) when frequency departure is less than ± fFireWhen, according to system loading and the frequency character of load coefficient during disturbance, calculate with
The related natural frequency coefficient of load, obtains first paragraph β value：
Wherein：βL*The active power change caused by frequency change 1%, according to system load level and constituent etc. really
It is fixed, typically take 1%~3%；PLe、feRespectively rated load and system nominal frequency；
2) frequency departure is between ± fFireWith ± fWaterBetween when, according to thermoelectricity and total installed capacity of hydropower ratio, approximately obtain second segment β
Value：
β2≈KGngh+KL≈(β-KL)ngh+KL (3)
Wherein：KG、KLThe respectively frequency characteristic coefficient of generator and load；nghTotal installed capacity is accounted for for thermoelectricity installed capacity to hold
The ratio of amount；
3) frequency departure exceedes ± fWaterWhen, " mains frequency is responded and natural frequency computational methods " recommended using NERC,
Calculate the 3rd section of β value：
Wherein：ΔPA、ΔfAPower offset value and exemplary frequency deviation values respectively before system disturbance；ΔPB、ΔfBRespectively
After system disturbance speed regulator follow closely action make frequency tend towards stability and AGC not yet act before power offset value and frequency departure
Value；
Further, the computational methods that spare capacity is adjusted in the step S2 are as follows：
Method 1：Based on the regulation spare capacity acquiring method under annual peak load percentage β value,
Certain percentage is multiplied by with annual peak load and obtains FR coefficient β value, according to annual frequency per minute
Rate offset data, active power deviation delta P probability distribution is obtained using formula Δ P=β × Δ f, adjusted standby by R=- Δs P
Capacity probability distribution, so as to try to achieve the regulation spare capacity under certain cumulative probability；
Method 2：Based on the regulation spare capacity acquiring method under annual actual load percentage β value,
Certain percentage is multiplied by with annual actual load and obtains FR coefficient β value, according to annual frequency per minute
Rate offset data, active power deviation delta P probability distribution is obtained using formula Δ P=β × Δ f, so as to try to achieve general in certain accumulation
Regulation spare capacity under rate；
Method 3：Spare capacity acquiring method is adjusted based on the three-stage under annual peak load percentage β value,
Certain percentage is multiplied by with annual peak load and obtains every band frequency characteristic response factor beta1、β2、β3, according to difference
Unit primary frequency modulation dead band is boundary, and using different FR coefficients, active power deviation delta is calculated using formula (5)
P probability distribution；
Method 4：Spare capacity acquiring method is adjusted based on the three-stage under annual actual load percentage β value,
Certain percentage is multiplied by with annual actual load and obtains every band frequency characteristic response factor beta1*、β2*、β3*, according to not
It is boundary with unit primary frequency modulation dead band, using different FR coefficients, active power deviation is calculated using formula (5)
Δ P probability distribution.
Method 5：Three-stage regulation spare capacity acquiring method is mixed,
Certain percentage is actually multiplied by with whole year and obtains first paragraph FR factor beta1*, with annual peak load
It is multiplied by certain percentage and obtains second and third band frequency characteristic response factor beta2、β3, it is boundary according to different unit primary frequency modulation dead bands
Limit, using different FR coefficients, active power deviation delta P probability distribution is calculated using formula (7).
Five kinds of methods of the above obtain active power deviation probability distribution, with reference to the result of calculation under different condition, probability theory
In 3 σ principles, so as to try to achieve the regulation spare capacity under certain cumulative probability.
By using the regulation spare capacity appraisal procedure proposed by the invention based on frequency shift (FS), can effectively it instruct
Power System Planning, runing adjustment spare capacity are chosen.The present invention is to adapting to new electric reshaping gesture and new energy development, being to meeting
System needs and power network future development has universality and practicality.
Brief description of the drawings
Fig. 1 is the flow chart of the present invention.
Fig. 2 is frequency fluctuation typical curve after generating set failure.
Fig. 3 is bulk power grid frequency shift (FS) probability distribution graph.
Fig. 4 is bulk power grid frequency shift (FS) cumulative probability curve.
Embodiment
Accompanying drawing being given for example only property explanation, it is impossible to be interpreted as the limitation to this patent；It is attached in order to more preferably illustrate the present embodiment
Scheme some parts to have omission, zoom in or out, do not represent the size of actual product；To those skilled in the art,
Some known features and its explanation may be omitted and will be understood by accompanying drawing.Being given for example only property of position relationship described in accompanying drawing
Explanation, it is impossible to be interpreted as the limitation to this patent.
As shown in figure 1, a kind of standby regulation spare capacity assessment side of bulk power grid load based on frequency shift (FS) probability distribution
Method, comprises the following steps：
Step S1：Under actual electric network running situation, unit output can not possibly be all the time all equal with load, the two it
Between always there is active power deviation, so as to cause system frequency excursion rated value, and build based on this inclined based on frequency
Move the regulation spare capacity appraisal procedure of probability distribution；
Step S2：Load fluctuation causes system frequency excursion value to be Δ f, according to power system resultant frequency characteristics, can be with
Trying to achieve system active power deviation is：
Δ P=β × Δ f (1), as adjust spare capacity,
Wherein：β is power system frequency response characteristic coefficient, and unit is MW/0.1Hz；
Wherein, step S1 power systems, which are once adjusted, causes system frequency excursion, and specific system change process is as follows：
As shown in Fig. 2 Fig. 2 is power system frequency fluctuation typical curve after certain generating set failure, it is contemplated that load
The reduction exerted oneself in generating set of increase Approximate Equivalent, system frequency fluctuation situation after load increase suddenly can be described.When certain
One moment (A points in figure) system loading increases suddenly, causes generator output to be less than load, negative active power deviation delta occurs
P, at this moment frequency decline rapidly, power system, which is once adjusted, to play a role rapidly.(X points in figure) peak frequency of arriving is inclined after 10 seconds
Shifting amount, after 30s occurs for event (B points in figure), once adjustment terminates.The secondary adjustment of frequency is then proceeded by from B points,
AGC control unit increases are exerted oneself, and are occurred to event after 60s (C points in figure), system frequency starts slowly to 50Hz to recover；
Wherein, for system frequency excursion value Δ f acquiring methods, system frequency excursion value-based algorithm in step s 2 is such as
Under：
Power system is per second after being filtered in real time to system frequency interference signal by SCADA system in actual motion
Clock is once sampled, to meet nyquist sampling theorem (as sample frequency fs.maxMore than highest frequency f in signalmax2
Times when (fs.max≥2fmax), the data signal after sampling just intactly remains the information in primary signal.), this research is adopted
A frequency values are obtained with a stepped-frequency signal per minute, can be in the hope of can obtain using a large amount of frequency shift (FS) record datas
To active power deviation delta P probability distribution situation；
Wherein, for power system frequency response characteristic coefficient acquiring method, time-varying and non-during frequency response characteristic factor beta
Linear, the β algorithms in step S2 are as follows：
Method 1, FR coefficient is taken as the certain percentage of annual peak load：
In existing achievement in research, it is maximum that FR characteristic coefficient is all taken as the year by most of power system
The certain percentage of load, i.e. peak load % (MW/0.1Hz).This method is typically all (as sent out by typical frequencies surge events
Group of motors catastrophic failure, large-disturbance test), FR factor beta is asked for according to frequency variation curve.Most of control zones
FR coefficient between 1%~1.5% (MW/0.1Hz) of annual highest forecast demand；
Method 2, FR coefficient is taken as the certain percentage of actual load：
Consider the time-varying characteristics of FR factor beta, certain percentage is multiplied by using actual load and obtained；
Method 3, three-stage FR Y-factor method Y
Different primary frequency modulation dead bands are generally set in view of thermoelectricity, Hydropower Unit in power system, with reference to pertinent literature,
Beta coefficient can be configured for three sections according to fired power generating unit, Hydropower Unit Regulation dead-band as boundary point.Thermoelectricity, Hydropower Unit
The dead band of primary frequency modulation is respectively ± fFire, ± fWater, three-stage FR factor beta establishing method is as follows：
1) when frequency departure is less than ± fFireWhen, according to system loading and the frequency character of load coefficient during disturbance, calculate with
The related natural frequency coefficient of load, obtains first paragraph β value：
Wherein：βL*The active power change caused by frequency change 1%, according to system load level and constituent etc. really
It is fixed, typically take 1%~3%；PLe、feRespectively rated load and system nominal frequency；
2) frequency departure is between ± fFireWith ± fWaterBetween when, according to thermoelectricity and total installed capacity of hydropower ratio, approximately obtain second segment β
Value：
β2≈KGngh+KL≈(β-KL)ngh+KL (3)
Wherein：KG、KLThe respectively frequency characteristic coefficient of generator and load；nghTotal installed capacity is accounted for for thermoelectricity installed capacity
The ratio of (thermoelectricity and water power) capacity.
3) frequency departure exceedes ± fWaterWhen, " mains frequency is responded and natural frequency computational methods " recommended using NERC,
Calculate the 3rd section of β value：
Wherein：ΔPA、ΔfAPower offset value and exemplary frequency deviation values respectively before system disturbance；ΔPB、ΔfBRespectively
After system disturbance speed regulator follow closely action make frequency tend towards stability and AGC not yet act before power offset value and frequency departure
Value；
5. the power system frequency response described in regulation spare capacity according to claim 1 and claim 4 is special
Property coefficient acquiring method, it is characterised in that：It is characterized in that：The computational methods that spare capacity is adjusted in the step S2 are as follows：
Method 1：Based on the regulation spare capacity acquiring method under annual peak load percentage β value
Certain percentage is multiplied by with annual peak load and obtains FR coefficient β value, according to annual frequency per minute
Rate offset data, active power deviation delta P probability distribution is obtained using formula Δ P=β × Δ f, adjusted standby by R=- Δs P
Capacity probability distribution, so as to try to achieve the regulation spare capacity under certain cumulative probability；
Method 2：Based on the regulation spare capacity acquiring method under annual actual load percentage β value
Certain percentage is multiplied by with annual actual load and obtains FR coefficient β value, according to annual frequency per minute
Rate offset data, active power deviation delta P probability distribution is obtained using formula Δ P=β × Δ f, so as to try to achieve general in certain accumulation
Regulation spare capacity under rate；
Method 3：Spare capacity acquiring method is adjusted based on the three-stage under annual peak load percentage β value
Certain percentage is multiplied by with annual peak load and obtains every band frequency characteristic response factor beta1、β2、β3.According to difference
Unit primary frequency modulation dead band is boundary, and using different FR coefficients, active power deviation delta is calculated using formula (5)
P probability distribution；
Method 4：Spare capacity acquiring method is adjusted based on the three-stage under annual actual load percentage β value
Certain percentage is multiplied by with annual actual load and obtains every band frequency characteristic response factor beta1*、β2*、β3*.According to not
It is boundary with unit primary frequency modulation dead band, using different FR coefficients, active power deviation is calculated using formula (5)
Δ P probability distribution.
Method 5：Mix three-stage regulation spare capacity acquiring method
Certain percentage is multiplied by with annual actual load and obtains first paragraph FR factor beta1*, with annual maximum
Load is multiplied by certain percentage and obtains second and third band frequency characteristic response factor beta2、β3.According to different unit primary frequency modulation dead bands
For boundary, using different FR coefficients, active power deviation delta P probability distribution is calculated using formula (7).
Five kinds of methods of the above obtain active power deviation probability distribution, with reference to the result of calculation under different condition, probability theory
In 3 σ principles, so as to try to achieve the regulation spare capacity under certain cumulative probability.
By taking certain bulk power grid as an example, the regulation spare capacity appraisal procedure based on frequency shift (FS) comprises the following steps：
(1) frequency shift (FS) cumulative probability distribution situation.
According to bulk power grid frequency log per minute, certain time bulk power grid frequency shift (FS) probability distribution is as shown in Figure 3.
Bulk power grid frequency accumulation probability distribution graph is obtained by frequency shift (FS) probability distribution as shown in Figure 4：When cumulative probability is
When 0.90, as frequency offseting value Δ f=0.034Hz；When cumulative probability is 0.95, frequency offseting value Δ f=0.037Hz；When
When cumulative frequency is 0.99, frequency offseting value Δ f=0.045Hz；When cumulative frequency is 0.9999, frequency offseting value Δ f=
0.067Hz。
(2) parameter beta is asked for.
1) a typical frequencies surge events, the tripping operation of certain DC bipolar, removal of load 1800MW occur for power network known to；A, B point
At intervals of 18 seconds, Δ P=-1800MW, Δ f=0.12Hz, the whole network β=1500MW/0.1Hz；Tripping operation moment the whole network load be
72791MW, can calculate β value and account for load percentage for 1500/72791 × 100%=2% (MW/0.1Hz)；
2) its time variation is considered, actual electric network FR factor beta is the 20% of load per minute, therefore by every point
Clock load data is multiplied by 20%；
3) three sections are configured
β1：Under accident conditions, the percentage that Static Load frequency characteristic coefficient accounts for load is 3.29%, then has：β1≈
0.329% × Load (MW/0.1Hz)；
β2：Thermoelectricity accounts for the 67% of total installation of generating capacity, i.e. ngh=67%, try to achieve：β21.4456% × Load of ≈ (MW/
0.1Hz)；
β3：β32% × Load of ≈ (MW/0.1Hz)；
2) power network peak load Load=119445MW, therefore, β1≈393MW/0.1Hz、β2≈1728MW/0.1Hz、β3
≈2389MW/0.1Hz。
(3) network load regulation spare capacity is assessed.
1) partial electric grid peak load is 80050MW, and it is 67.02% to account for bulk power grid peak load percentage, from system scale and frequency
Adjust from the point of view of responsibility, partial electric grid regulation spare capacity value should be the 67.02% of bulk power grid regulation spare capacity；
2) the regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution is passed through：
Method 1：Asked for based on the regulation spare capacity under annual peak load percentage β value
β value is：Network system is adjusted under β=119445 × 2%MW/0.1Hz=2389MW/0.1Hz, certain cumulative probability
Spare capacity value is as shown in table 1：
Method 2：Asked for based on the regulation spare capacity under actual load percentage β value per minute
β value is：The network system regulation spare capacity value such as institute of table 2 under the 20% of load per minute, certain cumulative probability
Show：
Method 3-5：Regulation spare capacity under three-stage β value is asked for, as shown in table 3：
Comprehensive five kinds of method result of calculations, the regulation spare capacity value of reply load fluctuation part, result of calculation collects
As shown in table 4：
(4) Load Regulation spare capacity recommends value
With reference to 3 σ principles in the result of calculation under different condition, probability theory, local power system loading undulating value can take
The 1.2% of peak load is most adjusted wholly to be annual.
Obviously, the above embodiment of the present invention is only intended to clearly illustrate example of the present invention, and is not pair
The restriction of embodiments of the present invention.For those of ordinary skill in the field, may be used also on the basis of the above description
To make other changes in different forms.There is no necessity and possibility to exhaust all the enbodiments.It is all this
Any modifications, equivalent substitutions and improvements made within the spirit and principle of invention etc., should be included in the claims in the present invention
Protection domain within.
Claims (5)
1. a kind of bulk power grid Load Regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution, it is characterised in that bag
Include following steps：
Step S1：The active power deviation of unit output and load causes system frequency excursion rated value, general based on frequency shift (FS)
The structure regulation spare capacity appraisal procedure of rate distribution；
Step S2：Load fluctuation generation system frequency offseting value Δ f, can be in the hope of being according to power system resultant frequency characteristics
System active power deviation be：
Δ P=β × Δ f (1), as adjust spare capacity,
Wherein：β is power system frequency response characteristic coefficient, and unit is MW/0.1Hz.
2. the regulation spare capacity appraisal procedure according to claim 1 based on frequency shift (FS), it is characterised in that：The step
In rapid S1, power system is once adjusted, causes system frequency excursion, specific system change process is as follows：
When a certain moment A system loadings increase suddenly, cause generator output to be less than load, negative active power deviation delta occur
P, at this moment frequency decline rapidly, power system, which is once adjusted, to play a role rapidly, and X has arrived peak frequency at the time of after certain time
Offset, a certain moment B after an event occurs, once adjustment terminates, and the secondary tune of frequency is then proceeded by from moment B
Whole, AGC control unit increases are exerted oneself, to Time To Event after at the time of C, system frequency starts slowly to 50Hz to recover.
3. the regulation spare capacity appraisal procedure according to claim 1 based on frequency shift (FS), it is characterised in that：The step
System frequency excursion value Δ f algorithms in rapid S2 are as follows：
Power system is in actual motion, and after being filtered in real time to system frequency interference signal by SCADA system, each second enters
Row is once sampled, to meet nyquist sampling theorem, as sample frequency fs.maxMore than or equal to highest frequency f in signalmax2
Times when, a frequency values are obtained using a stepped-frequency signal per minute, can be in the hope of using a large amount of frequency shift (FS) record datas
Active power deviation delta P probability distribution situation can be obtained.
4. the regulation spare capacity appraisal procedure according to claim 3 based on frequency shift (FS), it is characterised in that：Frequency is rung
It is time-varying and nonlinear to answer characteristic coefficient β, and the β algorithms in the step S2 are as follows：
Method 1, FR coefficient is taken as the certain percentage of annual peak load：
In existing achievement in research, FR characteristic coefficient is all taken as the annual peak load by most of power system
Certain percentage, i.e. peak load percentage, the FR coefficient of most of control zones is in annual highest forecast demand
Between 1%~1.5%；
Method 2, FR coefficient is taken as the certain percentage of actual load：
Consider the time-varying characteristics of FR factor beta, certain percentage is multiplied by using actual load and obtained；
Method 3, three-stage FR Y-factor method Y：
Different primary frequency modulation dead bands are generally set in view of thermoelectricity, Hydropower Unit in power system, can be with reference to pertinent literature
Beta coefficient is configured for three sections according to fired power generating unit, Hydropower Unit Regulation dead-band as boundary point, thermoelectricity, hydroelectric units primary
The dead band of frequency modulation is respectively ± fFire, ± fWater, three-stage FR factor beta establishing method is as follows：
1) when frequency departure is less than ± fFireWhen, according to the system loading and frequency character of load coefficient during disturbance, calculate and load
Related natural frequency coefficient, obtains first paragraph β value：
Wherein：βL*Active power change, is determined according to system load level and constituent etc. caused by frequency change 1%,
Typically take 1%~3%；PLe、feRespectively rated load and system nominal frequency；
2) frequency departure is between ± fFireWith ± fWaterBetween when, according to thermoelectricity and total installed capacity of hydropower ratio, approximately obtain second segment β value：
β2≈KGngh+KL≈(β-KL)ngh+KL (3)
Wherein：KG、KLThe respectively frequency characteristic coefficient of generator and load；nghTotal installation of generating capacity is accounted for for thermoelectricity installed capacity
Ratio；
3) frequency departure exceedes ± fWaterWhen, " mains frequency is responded and natural frequency computational methods " recommended using NERC is calculated
3rd section of β value：
Wherein：ΔPA、ΔfAPower offset value and exemplary frequency deviation values respectively before system disturbance；ΔPB、ΔfBRespectively system
After disturbance speed regulator follow closely action make frequency tend towards stability and AGC not yet act before power offset value and exemplary frequency deviation values.
5. the regulation spare capacity appraisal procedure according to claim 4 based on frequency shift (FS), it is characterised in that the step
The computational methods of regulation spare capacity are as follows in rapid S2：
Method 1：Based on the regulation spare capacity acquiring method under annual peak load percentage β value,
Certain percentage is multiplied by with annual peak load and obtains FR coefficient β value, it is inclined according to annual frequency per minute
Data are moved, active power deviation delta P probability distribution are obtained using formula Δ P=β × Δ f, by the adjusted spare capacities of R=- Δs P
Probability distribution, so as to try to achieve the regulation spare capacity under certain cumulative probability；
Method 2：Based on the regulation spare capacity acquiring method under annual actual load percentage β value,
Certain percentage is multiplied by with annual actual load and obtains FR coefficient β value, it is inclined according to annual frequency per minute
Data are moved, active power deviation delta P probability distribution are obtained using formula Δ P=β × Δ f, so as to try to achieve under certain cumulative probability
Regulation spare capacity；
Method 3：Spare capacity acquiring method is adjusted based on the three-stage under annual peak load percentage β value,
Certain percentage is multiplied by with annual peak load and obtains every band frequency characteristic response factor beta1、β2、β3, according to different units
Primary frequency modulation dead band is boundary, using different FR coefficients, calculates active power deviation delta P using formula (5) general
Rate is distributed；
Method 4：Spare capacity acquiring method is adjusted based on the three-stage under annual actual load percentage β value,
Certain percentage is multiplied by with annual actual load and obtains every band frequency characteristic response factor beta1*、β2*、β3*, according to different machines
Group primary frequency modulation dead band is boundary, and using different FR coefficients, active power deviation delta P is calculated using formula (5)
Probability distribution.
Method 5：Three-stage regulation spare capacity acquiring method is mixed,
Certain percentage is actually multiplied by with whole year and obtains first paragraph FR factor beta1*, one is multiplied by with annual peak load
Determine percentage and obtain second and third band frequency characteristic response factor beta2、β3, it is boundary according to different unit primary frequency modulation dead bands, uses
Different FR coefficients, active power deviation delta P probability distribution is calculated using formula (7).
Five kinds of methods of the above obtain active power deviation probability distribution, with reference to 3 σ in the result of calculation under different condition, probability theory
Principle, so as to try to achieve the regulation spare capacity under certain cumulative probability.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710540775.1A CN107134811A (en) | 2017-07-05 | 2017-07-05 | Network load regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710540775.1A CN107134811A (en) | 2017-07-05 | 2017-07-05 | Network load regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107134811A true CN107134811A (en) | 2017-09-05 |
Family
ID=59736108
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710540775.1A Pending CN107134811A (en) | 2017-07-05 | 2017-07-05 | Network load regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107134811A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110556844A (en) * | 2019-09-29 | 2019-12-10 | 云南电网有限责任公司电力科学研究院 | automatic frequency modulation method for asynchronous transmitting-end power grid |
| CN112308303A (en) * | 2020-10-22 | 2021-02-02 | 新奥数能科技有限公司 | High fault-tolerant energy supply group load scheduling method and device based on deviation distribution and terminal equipment |
| CN112564134A (en) * | 2020-12-22 | 2021-03-26 | 南方电网科学研究院有限责任公司 | Method, device, equipment and medium for configuring primary frequency modulation reserve capacity of power grid |
| CN112769147A (en) * | 2020-12-30 | 2021-05-07 | 国家电网有限公司 | Frequency stability early warning method based on power grid perception data |
| CN113890066A (en) * | 2021-11-18 | 2022-01-04 | 广东电网有限责任公司 | Frequency modulation method and device of multi-direct-current feed-in system based on energy storage system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102510060A (en) * | 2011-11-29 | 2012-06-20 | 武汉大学 | Computing method of frequency characteristic coefficient of electric system |
| CN104300535A (en) * | 2014-10-17 | 2015-01-21 | 华中电网有限公司 | In-network quick standby coordination and optimization scheduling method used after high-power vacancy accidents |
| CN104881570A (en) * | 2015-05-06 | 2015-09-02 | 东南大学 | Engineering calculating method used for simulating power system secondary frequency modulation dynamic process |
| CN104915723A (en) * | 2015-03-24 | 2015-09-16 | 国电南瑞科技股份有限公司 | Interruptible load participation rotation standby coordination optimization method considering frequency recovery effect |
| CN105098979A (en) * | 2014-12-19 | 2015-11-25 | 国网山东泗水县供电公司 | Automatic electric power scheduling system and method |
-
2017
- 2017-07-05 CN CN201710540775.1A patent/CN107134811A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102510060A (en) * | 2011-11-29 | 2012-06-20 | 武汉大学 | Computing method of frequency characteristic coefficient of electric system |
| CN104300535A (en) * | 2014-10-17 | 2015-01-21 | 华中电网有限公司 | In-network quick standby coordination and optimization scheduling method used after high-power vacancy accidents |
| CN105098979A (en) * | 2014-12-19 | 2015-11-25 | 国网山东泗水县供电公司 | Automatic electric power scheduling system and method |
| CN104915723A (en) * | 2015-03-24 | 2015-09-16 | 国电南瑞科技股份有限公司 | Interruptible load participation rotation standby coordination optimization method considering frequency recovery effect |
| CN104881570A (en) * | 2015-05-06 | 2015-09-02 | 东南大学 | Engineering calculating method used for simulating power system secondary frequency modulation dynamic process |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110556844A (en) * | 2019-09-29 | 2019-12-10 | 云南电网有限责任公司电力科学研究院 | automatic frequency modulation method for asynchronous transmitting-end power grid |
| CN112308303A (en) * | 2020-10-22 | 2021-02-02 | 新奥数能科技有限公司 | High fault-tolerant energy supply group load scheduling method and device based on deviation distribution and terminal equipment |
| CN112308303B (en) * | 2020-10-22 | 2024-03-08 | 新奥数能科技有限公司 | High fault tolerance energy supply group load scheduling method and device based on deviation distribution and terminal equipment |
| CN112564134A (en) * | 2020-12-22 | 2021-03-26 | 南方电网科学研究院有限责任公司 | Method, device, equipment and medium for configuring primary frequency modulation reserve capacity of power grid |
| CN112564134B (en) * | 2020-12-22 | 2023-02-24 | 南方电网科学研究院有限责任公司 | Method, device, equipment and medium for configuring primary frequency modulation reserve capacity of power grid |
| CN112769147A (en) * | 2020-12-30 | 2021-05-07 | 国家电网有限公司 | Frequency stability early warning method based on power grid perception data |
| CN112769147B (en) * | 2020-12-30 | 2022-04-15 | 国家电网有限公司 | Frequency stability early warning method based on power grid perception data |
| CN113890066A (en) * | 2021-11-18 | 2022-01-04 | 广东电网有限责任公司 | Frequency modulation method and device of multi-direct-current feed-in system based on energy storage system |
| CN113890066B (en) * | 2021-11-18 | 2023-07-28 | 广东电网有限责任公司 | Frequency modulation method and device for multi-DC feed-in system based on energy storage system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107134811A (en) | Network load regulation spare capacity appraisal procedure based on frequency shift (FS) probability distribution | |
| Seneviratne et al. | Frequency response due to a large generator loss with the increasing penetration of wind/PV generation–A literature review | |
| Bian et al. | Demand side contributions for system inertia in the GB power system | |
| Concordia et al. | Load shedding on an isolated system | |
| Ashton et al. | Application of phasor measurement units to estimate power system inertial frequency response | |
| CN104505855B (en) | A kind of wind farm grid-connected characteristic online rating system | |
| CN104242273A (en) | Self-adaptation positive sequence current quick breaking protection method for power grid containing inverter-based DG (distributed generation) | |
| Tan et al. | Multi-area automatic generation control scheme considering frequency quality in Southwest China grid: Challenges and solutions | |
| CN104793030A (en) | Supervision method of distributed photovoltaic power generation electric larceny | |
| CN104638638B (en) | A kind of safety on line stable tendency analysis method for bulk power grid | |
| US10544778B2 (en) | Method of operating a DFIG wind turbine under SSR | |
| Zhang et al. | Security assessment method for inertia and frequency stability of high proportional renewable energy system | |
| Shun et al. | Dynamic optimization of adaptive under-frequency load shedding based on WAMS | |
| CN110854852B (en) | Configuration method of key parameters of AGC main station in high-proportion hydropower area | |
| CN102684210A (en) | Method for calculating steady state power fluctuation amplitude of AC tie line of interconnected systems in two regions | |
| Adeen et al. | Statistical correlation between wind penetration and grid frequency variations in the Irish network | |
| Kim et al. | Conditional value-at-risk-based method for evaluating the economic risk of superconducting fault current limiter installation | |
| Zhang et al. | The emergency control strategies of short-run isolated island wind farm | |
| Yang et al. | Centralized adaptive under frequency load shedding schemes for smart grid using synchronous phase measurement unit | |
| Gbadamosi et al. | Evaluation of operational efficiency of Shiroro hydro-electric plant in Nigeria | |
| Dyśko | Risk assessment analysis to find optimum ROCOF protection settings | |
| Karki et al. | Reliability evaluation of a wind power delivery system using an approximate wind model | |
| CN109274112B (en) | Energy storage power station scheduling method considering regional stability | |
| Zhao et al. | Real-time Inertia Estimation Using Probing Method and PMU Measurement | |
| Das et al. | Data-Driven Inertia Estimation of Power Systems using Autoencoder-based Unsupervised Feature Learning |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170905 |