CN107800138A - A kind of deviation peak computational method and device based on mains frequency deviation variation rate - Google Patents
A kind of deviation peak computational method and device based on mains frequency deviation variation rate Download PDFInfo
- Publication number
- CN107800138A CN107800138A CN201711092192.3A CN201711092192A CN107800138A CN 107800138 A CN107800138 A CN 107800138A CN 201711092192 A CN201711092192 A CN 201711092192A CN 107800138 A CN107800138 A CN 107800138A
- Authority
- CN
- China
- Prior art keywords
- mains frequency
- frequency deviation
- variation rate
- deviation variation
- expression formula
- 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.)
- Granted
Links
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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/08—Synchronising of 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)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
A kind of deviation peak computational method based on mains frequency deviation variation rate provided in an embodiment of the present invention, including：Mains-power circuit is reduced to unit model system, reads corresponding circuit parameter, the expression formula of mains frequency deviation variation rate is calculated；Mains frequency deviation variation rate is set 0 and to calculate and obtaining proportionality coefficient α between mains frequency deviation peak value and grid disturbance amount；Time value in the expression formula of mains frequency deviation variation rate is set 0 and to calculate and obtaining proportionality coefficient β between the initial value of mains frequency deviation variation rate and grid disturbance amount；The ratio α/β obtained between mains frequency deviation peak value and the initial value of mains frequency deviation variation rate is calculated according to proportionality coefficient α and β；Detect the power network of initial time occurs for failure frequency departure rate of change and combination ratio α/β calculating acquisition mains frequency deviation peak value.The present invention can realize the quick calculating to mains frequency deviation peak value, and then provide support to improving mains frequency stability.
Description
Technical field
The present invention relates to electric power system stability control technical field, more particularly to one kind to be based on mains frequency deviation variation rate
Deviation peak computational method and device.
Background technology
China's power network Energy Base and load center are in contrary distribution, can be divided into sending end power network and receiving end power network.If
Alternating current interconnection between sending end power network and receiving end power network disconnects, and relies only on direct current and sends electric energy outside, is achieved that between regional power grid
Asynchronous interconnection.After asynchronous interconnection, the angle stability problem of sending end and receiving end is no longer closely related, even if sending direct current outside occurs event
Hinder, the change of the trend distribution between regional power grid will not cause power network generator rotor angle unstability of overall importance, for maintaining load center electricity
Net angle stability is significant.But after the asynchronous interconnection of power network, power shortage caused by DC Line Fault can not be by friendship arranged side by side
Stream interconnection is transported to receiving end power network, and the power shortage of receiving end power network can cause mains frequency to decline, be unfavorable for the frequency of power network
Stability.
South electric network is complexity most rambunctious in the world using remote, Large Copacity, ultrahigh voltage AC and DC mixing transmission of electricity
One of system.The Guangdong Power Grid inside south electric network realizes asynchronous interconnection, substantial amounts of power with Yunnan Power System by direct current at present
Guangdong Power Grid is transported to by Yunnan hydropower base by multiple-circuit line, Guangdong Power Grid is a typical receiving end power network.Due to straight
It is larger to flow the performance number of conveying, while without exchange interconnection, direct current hair between the Guangdong Power Grid of receiving end and the Yunnan Power System of sending end
Raw failure is stopped transport, and will cause the larger vacancy of Guangdong Power Grid internal power, causes Guangdong Power Grid frequency to decline.For Guangdong Power Grid
Frequency stabilization problem, the frequency departure peak value of power network is an important indicator after failure, influences the formulation of related steady prosecutor case.
Frequency stabilization influences the safe operation of power network, and Zhang Zhiqiang, Xu Youping, Yuan Rongxiang et al. are in entitled Large-Scale Interconnected area
Sending end mains frequency characteristic and high-frequency cutting machine scheme after the grid disconnection of domain (electric power network technique, 2015, volume 39,288-293)
The frequency departure analytic expression of one-of-a-kind system has been derived in document, have studied the high-frequency cutting machine scheme of Central China Power Grid；Zhao Qiang, Wang Li
It is quick, Liu Zhaoxu et al. entitled nationwide integrated power grid interacted system frequency characteristic and UFLS scheme (electric power network technique, 2009,
Volume 33,35-40) document in devise unit computation model, assay mains frequency transient stability characteristic and low frequency
The adaptability of off-load scheme.
But the studies above does not have the calculating for considering mains frequency deviation peak value, does not provide mains frequency deviation peak value
Quick calculation method.
The content of the invention
The embodiments of the invention provide a kind of deviation peak computational method and device based on mains frequency deviation variation rate,
Mains frequency deviation peak value can be quickly calculated, simplifies the calculating of mains frequency deviation peak value, solving prior art can not
The quick technical problem for calculating mains frequency deviation peak value.
A kind of deviation peak computational method based on mains frequency deviation variation rate provided in an embodiment of the present invention, including：
S1：Mains-power circuit is reduced to unit model system, corresponding circuit parameter is read, it is inclined that mains frequency is calculated
The expression formula of poor rate of change；
S2：Set mains frequency deviation variation rate and calculate acquisition as the expression formula of 0 and combination mains frequency deviation variation rate
Proportionality coefficient α between mains frequency deviation peak value and grid disturbance amount；
S3：The time value set in the expression formula of mains frequency deviation variation rate becomes as 0 and combination mains frequency deviation
The expression formula of rate calculates the proportionality coefficient β between the initial value for obtaining mains frequency deviation variation rate and grid disturbance amount；
S4：Calculated according to proportionality coefficient α and β and obtain the initial of mains frequency deviation peak value and mains frequency deviation variation rate
Ratio α/β between value；
S5：Detect the power network of initial time occurs for failure frequency departure rate of change and combination ratio α/β calculating acquisition electricity
Net frequency departure peak value.
Preferably, the step S1 is specifically included：
S101：Mains-power circuit is reduced to unit model system, the unit inertia constant in unit model system is read, is
Unite adjustment effect of load constant, generator work(frequency static characteristic coefficient and speed regulator generalized time constant this four circuit parameters；
S102：The grid disturbance of circuit parameter substitution unit model system is caused to the transmission function of mains frequency change,
Calculate the time solution for obtaining mains frequency deviation；
S103：Derived function is carried out to the time according to the time solution of mains frequency deviation, obtains mains frequency change of error
The expression formula of rate.
Preferably, the step S2 is specifically included：
S201：Mains frequency deviation variation rate is set as 0, it is inclined that the expression formula of mains frequency deviation variation rate is reduced to frequency
The expression formula of poor peak value；
S202：By the expression formula of frequency departure peak value divided by grid disturbance amount expression formula, mains frequency deviation peak value is obtained
The expression formula of proportionality coefficient α between grid disturbance amount；
S203：The expression formula for reading relevant parameter and passing ratio factor alpha calculates the concrete numerical value for obtaining proportionality coefficient α.
Preferably, the step S3 is specifically included：
S301：The time value in the expression formula of mains frequency deviation variation rate is set as 0, mains frequency deviation variation rate
Expression formula be reduced to mains frequency deviation variation rate initial value expression formula；
S302：By the expression formula of the initial value of mains frequency deviation variation rate divided by grid disturbance amount expression formula, electricity is obtained
The expression formula of proportionality coefficient β between the initial value and grid disturbance amount of net frequency departure rate of change；
S303：The expression formula for reading relevant parameter and passing ratio factor beta calculates the concrete numerical value for obtaining proportionality coefficient β.
Preferably, the expression formula of the mains frequency deviation variation rate is specially：
Wherein, M is unit inertia constant, and D is system loading mediating effect+6 constant, and K is generator work(frequency static characteristic
Coefficient, T are speed regulator generalized time constant.
Preferably, the grid disturbance cause mains frequency change transmission function be：
Wherein, M is unit inertia constant, and D is system loading mediating effect+6 constant, and K is generator work(frequency static characteristic coefficient,
T is speed regulator generalized time constant, Δ PeFor the active power variable quantity in system, that is, send the change of dc power, Δ P outsidemFor
Prime mover mechanical output variable quantity caused by active power change, Δ PDFor load power change caused by active power change
Amount.
Preferably, the step S102 is specifically included：
The change of initialization system active power turns to step disturbance, then Δ Pe(s)=Δ Pe/s；
The transmission function for causing mains frequency to change to grid disturbance solves to obtain the time solution of frequency departure：
Wherein,
Preferably, when setting mains frequency deviation variation rate as 0, Δ ω can be obtainedmax=α Δs Pe；
When setting the time value in the expression formula of mains frequency deviation variation rate as 0, it can obtain
Formula can be obtainedAnd frequency departure peak value is calculated according to this formula.
Preferably, when setting mains frequency deviation variation rate as 0, Δ ω can be obtainedmax=α Δs Pe；
When setting the time value in the expression formula of mains frequency deviation variation rate as 0, it can obtain
Formula can be obtainedAnd frequency departure peak value is calculated according to this formula.
A kind of deviation peak computational device based on mains frequency deviation variation rate provided by the invention, including a kind of calculating
Machine readable storage medium storing program for executing, is stored thereon with computer program, and side described above is realized when the computer program is executed by processor
The step of method.
As can be seen from the above technical solutions, the embodiment of the present invention has advantages below：
A kind of deviation peak computational method based on mains frequency deviation variation rate provided in an embodiment of the present invention, including：
Mains-power circuit is reduced to unit model system, reads corresponding circuit parameter, mains frequency deviation variation rate is calculated
Expression formula；Mains frequency deviation variation rate is set 0 and to calculate and obtaining between mains frequency deviation peak value and grid disturbance amount
Proportionality coefficient α；Time value in the expression formula of mains frequency deviation variation rate is set 0 and to calculate and obtaining mains frequency deviation
Proportionality coefficient β between the initial value and grid disturbance amount of rate of change；It is inclined that acquisition mains frequency is calculated according to proportionality coefficient α and β
Ratio α/β between the initial value of poor peak value and mains frequency deviation variation rate；Detect the power network of failure generation initial time
Frequency departure rate of change simultaneously combines ratio α/β calculating acquisition mains frequency deviation peak value.The present invention can be realized to mains frequency
The quick calculating of deviation peak value, and then provide support to improving mains frequency stability.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, without having to pay creative labor, may be used also
To obtain other accompanying drawings according to these accompanying drawings.
Fig. 1 is a kind of deviation peak computational method based on mains frequency deviation variation rate provided in an embodiment of the present invention
The schematic diagram of one embodiment；
Fig. 2 is a kind of deviation peak computational method based on mains frequency deviation variation rate provided in an embodiment of the present invention
The schematic diagram of another embodiment；
Fig. 3 is in a kind of deviation peak computational method based on mains frequency deviation variation rate provided in an embodiment of the present invention
The schematic diagram of unit model system；
Fig. 4 is the asynchronous interconnected network schematic diagram in twoth area of concrete application in application examples of the present invention；
Fig. 5 is mains frequency deviation schematic diagram in application examples of the present invention.
Embodiment
The embodiments of the invention provide a kind of deviation peak computational method and device based on mains frequency deviation variation rate,
Mains frequency deviation peak value can be quickly calculated, simplifies the calculating of mains frequency deviation peak value, solving prior art can not
The quick technical problem for calculating mains frequency deviation peak value.
To enable goal of the invention, feature, the advantage of the present invention more obvious and understandable, below in conjunction with the present invention
Accompanying drawing in embodiment, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that disclosed below
Embodiment be only part of the embodiment of the present invention, and not all embodiment.Based on the embodiment in the present invention, this area
All other embodiment that those of ordinary skill is obtained under the premise of creative work is not made, belongs to protection of the present invention
Scope.
A kind of referring to Fig. 1, deviation peak computational based on mains frequency deviation variation rate provided in an embodiment of the present invention
One embodiment of method, including：
101：Mains-power circuit is reduced to unit model system, corresponding circuit parameter is read, mains frequency is calculated
The expression formula of deviation variation rate；
102：Set mains frequency deviation variation rate and calculate acquisition as the expression formula of 0 and combination mains frequency deviation variation rate
Proportionality coefficient α between mains frequency deviation peak value and grid disturbance amount；
103：The time value set in the expression formula of mains frequency deviation variation rate becomes as 0 and combination mains frequency deviation
The expression formula of rate calculates the proportionality coefficient β between the initial value for obtaining mains frequency deviation variation rate and grid disturbance amount；
104：Calculated according to proportionality coefficient α and β and obtain mains frequency deviation peak value at the beginning of with mains frequency deviation variation rate
Ratio α/β between initial value；
105：Detect the power network of initial time occurs for failure frequency departure rate of change and combination ratio α/β calculating acquisition electricity
Net frequency departure peak value.
The present invention can realize the quick calculating to mains frequency deviation peak value, and then be carried to improving mains frequency stability
For supporting.
Above is to a kind of deviation peak computational side based on mains frequency deviation variation rate provided in an embodiment of the present invention
One embodiment of method is described in detail, below will be to provided in an embodiment of the present invention a kind of based on the change of mains frequency deviation
Another embodiment of the deviation peak computational method of rate is described in detail.
A kind of deviation peak computational method based on mains frequency deviation variation rate provided in an embodiment of the present invention it is another
Individual embodiment, including：
201：Mains-power circuit is reduced to unit model system, the unit inertia constant in unit model system is read, is
Unite adjustment effect of load constant, generator work(frequency static characteristic coefficient and speed regulator generalized time constant this four circuit parameters；
202：The grid disturbance of circuit parameter substitution unit model system is caused to the transmission function of mains frequency change, meter
Calculate the time solution for obtaining mains frequency deviation；
203：Derived function is carried out to the time according to the time solution of mains frequency deviation, obtains mains frequency deviation variation rate
Expression formula.
204：Mains frequency deviation variation rate is set as 0, it is inclined that the expression formula of mains frequency deviation variation rate is reduced to frequency
The expression formula of poor peak value；
Following formula (expression formula of mains frequency deviation variation rate) is reduced to multiplying for proportionality coefficient α and grid disturbance amount
Value, so as to draw mains frequency deviation variation rate be 0 when corresponding proportionality coefficient α；
205：By the expression formula of frequency departure peak value divided by grid disturbance amount expression formula, obtain mains frequency deviation peak value with
The expression formula of proportionality coefficient α between grid disturbance amount；
206：The expression formula for reading relevant parameter and passing ratio factor alpha calculates the concrete numerical value for obtaining proportionality coefficient α.
207：The time value in the expression formula of mains frequency deviation variation rate is set as 0, mains frequency deviation variation rate
Expression formula be reduced to mains frequency deviation variation rate initial value expression formula；
Multiply value by what the expression formula of mains frequency deviation variation rate was reduced to proportionality coefficient β and grid disturbance amountSo as to obtain the time value in the expression formula of mains frequency deviation variation rate for 0 when pair
The proportionality coefficient β answered.
208：By the expression formula of the initial value of mains frequency deviation variation rate divided by grid disturbance amount expression formula, power network is obtained
The expression formula of proportionality coefficient β between the initial value and grid disturbance amount of frequency departure rate of change；
209：The expression formula for reading relevant parameter and passing ratio factor beta calculates the concrete numerical value for obtaining proportionality coefficient β.
210：Calculated according to proportionality coefficient α and β and obtain mains frequency deviation peak value at the beginning of with mains frequency deviation variation rate
Ratio α/β between initial value；
211：Detect the power network of initial time occurs for failure frequency departure rate of change and combination ratio α/β calculating acquisition electricity
Net frequency departure peak value.
The expression formula of the mains frequency deviation variation rate is specially：
Wherein, M is unit inertia constant, and D is system loading mediating effect+6 constant, K
For generator work(frequency static characteristic coefficient, T is speed regulator generalized time constant.
The grid disturbance cause mains frequency change transmission function be：
Wherein, M is unit inertia constant, and D is system loading mediating effect+6 constant, and K is generator work(frequency static characteristic coefficient,
T is speed regulator generalized time constant, Δ PeFor the active power variable quantity in system, that is, send the change of dc power, Δ P outsidemFor
Prime mover mechanical output variable quantity caused by active power change, Δ PDFor load power change caused by active power change
Amount.
The change of initialization system active power turns to step disturbance, then Δ Pe(s)=Δ Pe/s；
The transmission function for causing mains frequency to change to grid disturbance solves to obtain the time solution of frequency departure：
Wherein,
Further, when setting mains frequency deviation variation rate as 0, Δ ωmax=α Δs Pe；
When setting the time value in the expression formula of mains frequency deviation variation rate as 0,
Formula can be obtained by merging two formula aboveIt is and inclined according to this formula calculating frequency
Poor peak value.
Further, when setting mains frequency deviation variation rate as 0, Δ ω can be obtainedmax=α Δs Pe；
When setting the time value in the expression formula of mains frequency deviation variation rate as 0, it can obtain
Formula can be obtained by merging two formula aboveAnd frequency departure is calculated according to this formula
Peak value.
The principle of the present invention will be explained below：
Present embodiment is illustrated using one-of-a-kind system as shown in Figure 3, can based on the one-of-a-kind system model shown in Fig. 3
Using obtain grid disturbance cause mains frequency change transmission function as：
In formula (1)：M is unit inertia constant, and D is system loading mediating effect+6 constant, and K is generator work(frequency static characteristic system
Number, T is speed regulator generalized time constant.ΔPeFor the active power variable quantity in system, that is, send the change of dc power outside.Δ
PmFor prime mover mechanical output variable quantity, Δ P caused by active power changeDBecome for load power caused by active power change
Change amount.
If the change of system active power turns to step disturbance, you can is expressed as Δ Pe(s)=Δ Pe/ s, to shown in formula (1)
Dynamic characteristic of power frequency transmission function is solved, and the time solution for obtaining frequency departure is：
In formula (2)：
Frequency departure in formula (2) can be obtained on the derivation of time t below：
It can be obtained by formula (3)：As t=0,
From formula (4), the rate of change that the system frequency deviation at moment occurs for disturbance is directly proportional to the size of disturbance quantity, with
System inertia constant is inversely proportional.
When system frequency deviation reaches peak value, the rate of change of system frequency deviation is zero in corresponding (3), can obtain peak value
Time meets：
N=0 in formula, 1,2
Therefore the time to peak of frequency departure is only relevant with power network inherent parameters, unrelated with disturbance quantity size.Substitution formula (2)
It is directly proportional to grid disturbance amount size that frequency departure peak value can be obtained, be represented by：
Δωmax=α Δs Pe (6)
When usual power network disturbs, disturbance quantity size delta PeIt not can determine that typically, so directly can not be determined by formula (6)
The frequency departure peak value of power network.But after disturbance occurs, the rate of change of mains frequency deviation but can be measured quickly, with reference to
Formula (4) can obtain with formula (6)：
Therefore from formula (7), the frequency departure peak value after disturbance occurs can be become by disturbing the frequency departure at moment
Rate is tried to achieve.In asynchronous interconnected network, DC Line Fault can be represented with step disturbance, so the low frequency peak value of receiving end power network
It can be tried to achieve by formula (7).
The application examples of the present invention will be described below：
The power system that present embodiment is studied is the asynchronous interconnected network in twoth area as shown in Figure 2, and research purpose is meter
Calculate the frequency departure peak value that section dc power changes receiving end power network in time domain 2.For the power system, using based on power network
The mains frequency deviation peak computational method flow of frequency departure rate of change is as follows：
(1) region 2 is equivalent with the one-of-a-kind system in Fig. 1, it is M=10, K=20, T=5, D to obtain one-of-a-kind system parameter
=2.；
(2) parameter tried to achieve in step (1) is updated to one-of-a-kind system transmission function, solution obtains mains frequency deviation
Time solution；
(3) time solution for the frequency departure tried to achieve in step (2) is subjected to derivation to the time, obtains the change of mains frequency deviation
The expression formula of rate；
(4) when mains frequency deviation variation rate is 0, mains frequency deviation reaches peak value, by the electricity obtained in step (3)
The expression formula of net frequency departure rate of change can be in the hope of the proportionality coefficient between mains frequency deviation peak value and grid disturbance amount
0.133；
(5) time in the expression formula of the mains frequency deviation variation rate obtained in step (3) is taken as 0, can be in the hope of
Proportionality coefficient between the initial value and grid disturbance amount of mains frequency deviation variation rate is 0.1；
(6), can be in the hope of mains frequency deviation peak value and mains frequency deviation variation rate according to step (4) and step (5)
Initial value between proportionality coefficient be 1.33.
(7) locking emulation, power network initial time frequency departure are carried out at the 1s moment to a direct current in power network shown in Fig. 2
Rate of change is -1.1Hz/s, and it is -1.463Hz that can try to achieve mains frequency deviation peak value with reference to the proportionality coefficient in step (6).It is actual
The mains frequency aberration curve of emulation is as shown in figure 3, simulation result is -1.464Hz.Proof is asked for based on frequency departure rate of change
The correctness of frequency departure peak method.
Mains frequency deviation computational methods of the present invention have advantages below：
(a) the inventive method has very strong theory analysis foundation, explicitly points out mains frequency deviation peak value and power network frequency
The inner link of rate deviation variation rate, the frequency departure rate of change in method occur for failure after mains frequency deviation variation rate
Initial value.
(b) present invention has practicality in actual electric network, because mains frequency deviation peak value and mains frequency deviation become
It is proportionate relationship between rate, can directly tries to achieve mains frequency deviation peak value and power network by l-G simulation test in actual electric network
Proportionality coefficient between frequency departure rate of change, quickly it can be calculated when an error occurs based on mains frequency deviation variation rate
To the frequency departure peak value of power network.
The aspect reason of summary two, it can be seen that：Compared with prior art, the present invention can realize inclined to mains frequency
The quick calculating of poor peak value, and then provide support to improving mains frequency stability.
In addition, a kind of deviation peak computational device based on mains frequency deviation variation rate provided by the invention, including one
Kind computer-readable recording medium, is stored thereon with computer program, is realized as above when the computer program is executed by processor
The step of method stated.
Described above, the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although with reference to before
Embodiment is stated the present invention is described in detail, it will be understood by those within the art that：It still can be to preceding
State the technical scheme described in each embodiment to modify, or equivalent substitution is carried out to which part technical characteristic；And these
Modification is replaced, and the essence of appropriate technical solution is departed from the spirit and scope of various embodiments of the present invention technical scheme.
Claims (10)
- A kind of 1. deviation peak computational method based on mains frequency deviation variation rate, it is characterised in that including：S1：Mains-power circuit is reduced to unit model system, reads corresponding circuit parameter, the change of mains frequency deviation is calculated The expression formula of rate；S2：Set mains frequency deviation variation rate and calculate acquisition power network as the expression formula of 0 and combination mains frequency deviation variation rate Proportionality coefficient α between frequency departure peak value and grid disturbance amount；S3：Time value in the expression formula of mains frequency deviation variation rate is set 0 and to combine mains frequency deviation variation rate Expression formula calculate obtain mains frequency deviation variation rate initial value and grid disturbance amount between proportionality coefficient β；S4：Calculated according to proportionality coefficient α and β obtain mains frequency deviation peak value and mains frequency deviation variation rate initial value it Between ratio α/β；S5：Detect the power network of initial time occurs for failure frequency departure rate of change and combination ratio α/β calculating acquisition power network frequency Rate deviation peak value.
- 2. a kind of deviation peak computational method based on mains frequency deviation variation rate according to claim 1, its feature It is, the step S1 is specifically included：S101：Mains-power circuit is reduced to unit model system, reads the unit inertia constant in unit model system, system is born This four circuit parameters of lotus mediating effect+6 constant, generator work(frequency static characteristic coefficient and speed regulator generalized time constant；S102：The grid disturbance of circuit parameter substitution unit model system is caused to the transmission function of mains frequency change, calculated Obtain the time solution of mains frequency deviation；S103：Derived function is carried out to the time according to the time solution of mains frequency deviation, obtains mains frequency deviation variation rate Expression formula.
- 3. a kind of deviation peak computational method based on mains frequency deviation variation rate according to claim 1, its feature It is, the step S2 is specifically included：S201：Mains frequency deviation variation rate is set as 0, the expression formula of mains frequency deviation variation rate is reduced to frequency departure peak The expression formula of value；S202：By the expression formula of frequency departure peak value divided by grid disturbance amount expression formula, mains frequency deviation peak value and electricity are obtained The expression formula of proportionality coefficient α between net disturbance quantity；S203：The expression formula for reading relevant parameter and passing ratio factor alpha calculates the concrete numerical value for obtaining proportionality coefficient α.
- 4. a kind of deviation peak computational method based on mains frequency deviation variation rate according to claim 1, its feature It is, the step S3 is specifically included：S301：The time value in the expression formula of mains frequency deviation variation rate is set as 0, the table of mains frequency deviation variation rate The expression formula of the initial value of mains frequency deviation variation rate is reduced to up to formula；S302：By the expression formula of the initial value of mains frequency deviation variation rate divided by grid disturbance amount expression formula, power network frequency is obtained The expression formula of proportionality coefficient β between the initial value and grid disturbance amount of rate deviation variation rate；S303：The expression formula for reading relevant parameter and passing ratio factor beta calculates the concrete numerical value for obtaining proportionality coefficient β.
- 5. a kind of deviation peak computational method based on mains frequency deviation variation rate according to claim 1, its feature It is, the expression formula of the mains frequency deviation variation rate is specially：Wherein, M is unit inertia constant, and D is system loading mediating effect+6 constant, and K is generator work(frequency static characteristic Coefficient, T are speed regulator generalized time constant.
- 6. a kind of deviation peak computational method based on mains frequency deviation variation rate according to claim 2, its feature It is, the transmission function that the grid disturbance causes mains frequency to change is：<mrow> <mfrac> <mrow> <mi>&Delta;</mi> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&Delta;P</mi> <mi>e</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mi>s</mi> <mi>T</mi> <mo>+</mo> <mn>1</mn> </mrow> <mrow> <msup> <mi>MTs</mi> <mn>2</mn> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mi>M</mi> <mo>+</mo> <mi>D</mi> <mi>T</mi> <mo>)</mo> </mrow> <mi>s</mi> <mo>+</mo> <mrow> <mo>(</mo> <mi>D</mi> <mo>+</mo> <mi>K</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>;</mo> </mrow>Wherein, M is unit inertia constant, and D is system loading mediating effect+6 constant, and K is generator work(frequency static characteristic coefficient, and T is Speed regulator generalized time constant, Δ PeFor the active power variable quantity in system, that is, send the change of dc power, Δ P outsidemTo have Prime mover mechanical output variable quantity caused by work(changed power, Δ PDFor load power variable quantity caused by active power change.
- 7. a kind of deviation peak computational method based on mains frequency deviation variation rate according to claim 6, its feature It is, the step S102 is specifically included：The change of initialization system active power turns to step disturbance, then Δ Pe(s)=Δ Pe/s；The transmission function for causing mains frequency to change to grid disturbance solves to obtain the time solution of frequency departure：Wherein,
- 8. a kind of deviation peak computational method based on mains frequency deviation variation rate according to claim 7, its feature It is, when setting mains frequency deviation variation rate as 0, Δ ωmax=α Δs Pe；When setting the time value in the expression formula of mains frequency deviation variation rate as 0,Formula can be obtainedAnd frequency departure peak value is calculated according to this formula.
- 9. a kind of deviation peak computational method based on mains frequency deviation variation rate according to claim 7, its feature It is, when setting mains frequency deviation variation rate as 0, Δ ωmax=α Δs Pe；When setting the time value in the expression formula of mains frequency deviation variation rate as 0,According to formulaCalculate frequency departure peak value.
- 10. a kind of deviation peak computational device based on mains frequency deviation variation rate, it is characterised in that including a kind of computer Readable storage medium storing program for executing, be stored thereon with computer program, when the computer program is executed by processor realize as claim 1 to In 9 the step of method described in any one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711092192.3A CN107800138B (en) | 2017-11-08 | 2017-11-08 | Deviation peak value calculation method and device based on power grid frequency deviation change rate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711092192.3A CN107800138B (en) | 2017-11-08 | 2017-11-08 | Deviation peak value calculation method and device based on power grid frequency deviation change rate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107800138A true CN107800138A (en) | 2018-03-13 |
| CN107800138B CN107800138B (en) | 2020-08-11 |
Family
ID=61549209
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711092192.3A Active CN107800138B (en) | 2017-11-08 | 2017-11-08 | Deviation peak value calculation method and device based on power grid frequency deviation change rate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107800138B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108664749A (en) * | 2018-05-18 | 2018-10-16 | 西安理工大学 | A kind of optimization method of dynamo governor frequency feedback coefficient |
| CN108988359A (en) * | 2018-08-01 | 2018-12-11 | 深圳供电局有限公司 | Non-linear automatic power generation control method, device, equipment and medium |
| 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 |
| CN115480122A (en) * | 2022-10-19 | 2022-12-16 | 云南电网有限责任公司电力科学研究院 | Method, device, equipment and medium for judging disturbance of power quality |
| CN117293860A (en) * | 2023-09-28 | 2023-12-26 | 暨南大学 | Frequency safety index evaluation method, system and storage medium of power system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103066620A (en) * | 2012-12-24 | 2013-04-24 | 中国电力科学研究院 | Design method of automatic generation control model under intermittent energy grid-connection |
| US8922043B1 (en) * | 2014-03-06 | 2014-12-30 | Industrial Cooperation Foundation Chonbuk National University | Time variant droop based inertial control method for wind generator |
| CN104993489A (en) * | 2015-07-09 | 2015-10-21 | 华北电力大学 | Method for determining optimal load shedding amount under circumstances of isolated grid emergency |
| CN105633988A (en) * | 2016-01-18 | 2016-06-01 | 华北电力科学研究院有限责任公司 | Method and device for energy storage system for participating in primary frequency control of power grid |
| CN105680483A (en) * | 2016-04-22 | 2016-06-15 | 广东电网有限责任公司电力科学研究院 | Virtual synchronous machine damping configuration method and device |
| CN106786714A (en) * | 2016-12-07 | 2017-05-31 | 中国南方电网有限责任公司 | A kind of asynchronous power network direct current locking dynamic frequency monitoring and Security Checking method |
-
2017
- 2017-11-08 CN CN201711092192.3A patent/CN107800138B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103066620A (en) * | 2012-12-24 | 2013-04-24 | 中国电力科学研究院 | Design method of automatic generation control model under intermittent energy grid-connection |
| US8922043B1 (en) * | 2014-03-06 | 2014-12-30 | Industrial Cooperation Foundation Chonbuk National University | Time variant droop based inertial control method for wind generator |
| CN104993489A (en) * | 2015-07-09 | 2015-10-21 | 华北电力大学 | Method for determining optimal load shedding amount under circumstances of isolated grid emergency |
| CN105633988A (en) * | 2016-01-18 | 2016-06-01 | 华北电力科学研究院有限责任公司 | Method and device for energy storage system for participating in primary frequency control of power grid |
| CN105680483A (en) * | 2016-04-22 | 2016-06-15 | 广东电网有限责任公司电力科学研究院 | Virtual synchronous machine damping configuration method and device |
| CN106786714A (en) * | 2016-12-07 | 2017-05-31 | 中国南方电网有限责任公司 | A kind of asynchronous power network direct current locking dynamic frequency monitoring and Security Checking method |
Non-Patent Citations (1)
| Title |
|---|
| 张建新等: "基于频率变化率的低频减载方案在南方电网应用的探讨", 《广东电力》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108664749A (en) * | 2018-05-18 | 2018-10-16 | 西安理工大学 | A kind of optimization method of dynamo governor frequency feedback coefficient |
| CN108664749B (en) * | 2018-05-18 | 2021-11-16 | 西安理工大学 | Optimization method for frequency feedback coefficient of generator speed regulator |
| CN108988359A (en) * | 2018-08-01 | 2018-12-11 | 深圳供电局有限公司 | Non-linear automatic power generation control method, device, equipment and medium |
| 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 |
| CN115480122A (en) * | 2022-10-19 | 2022-12-16 | 云南电网有限责任公司电力科学研究院 | Method, device, equipment and medium for judging disturbance of power quality |
| CN117293860A (en) * | 2023-09-28 | 2023-12-26 | 暨南大学 | Frequency safety index evaluation method, system and storage medium of power system |
| CN117293860B (en) * | 2023-09-28 | 2024-04-23 | 暨南大学 | Frequency safety index evaluation method, system and storage medium of power system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107800138B (en) | 2020-08-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107800138A (en) | A kind of deviation peak computational method and device based on mains frequency deviation variation rate | |
| Zheng et al. | Comprehensive control strategy of virtual synchronous generator under unbalanced voltage conditions | |
| ES2382786B1 (en) | METHOD AND VOLTAGE CONTROL SYSTEM OF AN ELECTRICAL GENERATION CENTER AND WIND FARM | |
| CN106712067B (en) | A kind of DC grid small signal stability determination method | |
| US20200401740A1 (en) | Aggregated model of large-scale wind farms for power system simulation software tools | |
| CN103760400B (en) | Voltage stabilization on-line monitoring method based on wide area measurement system | |
| CN103972899B (en) | A kind of STATCOM access point voltage compensating method | |
| CN109586337A (en) | VSC grid-connected system sub-synchronous oscillation methods of risk assessment based on Modeling In Frequency Domain | |
| CN102904266B (en) | Method for determining inactive compensation capacity network adaptability of wind power plant | |
| CN104979833B (en) | Consider the static electric voltage stability determination methods of area power grid reactive power support ability | |
| CN110137976A (en) | The stability judging method and device of voltage-source type converter station access multi-infeed systems | |
| CN109167380A (en) | A kind of judgment method for the multi-infeed systems stability accessing voltage-source type converter station | |
| CN103326354A (en) | Closed loop power distribution network tide estimation method based on load moment | |
| CN104485670B (en) | The control method of voltage sensitivity industrial load time-varying damping characteristic in island network | |
| CN109217302A (en) | The evaluation method of Traditional DC receiving end grid strength in a kind of mixing multi-infeed systems | |
| Bidram et al. | Capacitor design considering first swing stability of distributed generations | |
| CN110034567A (en) | Power distribution network short circuit type voltage dip characteristic quantity calculating method based on propagation characteristic | |
| Bhole et al. | Improvement of voltage stability in power system by using SVC and STATCOM | |
| CN103715704B (en) | A kind of micro-electrical network common bus Voltage unbalance inhibition method | |
| CN103529388B (en) | A kind of method and device measuring difference modulation rate of power generation set | |
| Zhang et al. | Research overview of sub-synchronous oscillation in DFIG-BASED wind farms connected to grid | |
| Wang et al. | A novel directional element for transmission line connecting inverter-interfaced renewable energy power plant | |
| Shao et al. | Medium-frequency and sub-synchronous oscillation analysis of direct-drive wind farms connected to the parallel-compensated AC grid | |
| CN102684210A (en) | Method for calculating steady state power fluctuation amplitude of AC tie line of interconnected systems in two regions | |
| CN102868177B (en) | Automatic identification and correction method for system impedance of wind power plant grid-connected point |
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 |