CN108565854A - Analogy method and system for reappearing power system frequency characteristic - Google Patents

Analogy method and system for reappearing power system frequency characteristic Download PDF

Info

Publication number
CN108565854A
CN108565854A CN201810361848.5A CN201810361848A CN108565854A CN 108565854 A CN108565854 A CN 108565854A CN 201810361848 A CN201810361848 A CN 201810361848A CN 108565854 A CN108565854 A CN 108565854A
Authority
CN
China
Prior art keywords
power
system frequency
model
moment
value
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
Application number
CN201810361848.5A
Other languages
Chinese (zh)
Other versions
CN108565854B (en
Inventor
刘辉
孙大卫
田博
葛俊
宋鹏
江浩
程雪坤
刘海涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
North China Electric Power Research Institute Co Ltd
State Grid Jibei Electric Power Co Ltd
Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
North China Electric Power Research Institute Co Ltd
State Grid Jibei Electric Power Co Ltd
Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, North China Electric Power Research Institute Co Ltd, State Grid Jibei Electric Power Co Ltd, Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201810361848.5A priority Critical patent/CN108565854B/en
Publication of CN108565854A publication Critical patent/CN108565854A/en
Application granted granted Critical
Publication of CN108565854B publication Critical patent/CN108565854B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Turbines (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

The present invention provides a kind of analogy methods and system for reappearing power system frequency characteristic, are related to technical field of power systems.Method includes:According to the system frequency value of current time t, load active power is determined using pre-set load model;According to pre-set fuel reference signal and pitch aperture, air inlet magnitude is determined using pre-set main vapour pressure model;The mechanical output of steam turbine model is determined using pre-set steam turbine model;Electromagnetic power is generated using pre-set generating set model;According to electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element;According to load active power and fired power generating unit output power, the system frequency predicted value that equation determines the subsequent time t+1 of current time t is waved using synchronous machine, and using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment;Power system frequency characteristic is determined according to the system frequency predicted value at each moment.

Description

Analogy method and system for reappearing power system frequency characteristic
Technical field
The present invention relates to technical field of power systems more particularly to a kind of simulations for reappearing power system frequency characteristic Method and system.
Background technology
Currently, with the fast development of the new energy technologies such as wind-power electricity generation, photovoltaic generation, new energy technology generates electricity in power grid In accounting increasingly rise.As new energy technology power generation accounting in power grid is higher and higher, the frequency characteristic in electric system There is larger differences with the frequency characteristic of the electric system before based on fired power generating unit.
In order to study characteristic of the new energy technology power generation in power grid, need first to carry out existing power system frequency characteristic Reproduction research.However, existing power system digital simulation model is difficult to reflect the frequency characteristic of electric system of new generation, especially It is to be difficult to reappear frequency characteristic of the electric system after there is high-power vacancy failure.
Invention content
The embodiment of the present invention provides a kind of analogy method and system for reappearing power system frequency characteristic, to realize Accurate reproduction to current power system response has established base for characteristic of the research new energy technology power generation in power grid Plinth.
In order to achieve the above objectives, the present invention adopts the following technical scheme that:
A kind of analogy method for reappearing power system frequency characteristic, including:
According to the system frequency value of current time t, load active power is determined using pre-set load model;
According to pre-set fuel reference signal and pitch aperture, using pre-set main vapour pressure model determine into Tolerance value;
According to the air inlet magnitude, the mechanical output of steam turbine model is determined using pre-set steam turbine model;
According to the mechanical output, electromagnetic power is generated using pre-set generating set model;
According to the electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element;
According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and is determined currently The system frequency predicted value of the subsequent time t+1 of moment t, and using the system frequency predicted value of subsequent time t+1 as the t+1 moment System frequency value;
Power system frequency characteristic is determined according to the system frequency predicted value at each moment.
Specifically, the system frequency value according to current time t, determines that load has using pre-set load model Work(power, including:
System frequency value of the system frequency predicted value of t moment as t moment is obtained, or t moment is obtained by sampling System frequency value of the system frequency actual value as t moment;
Using pre-set load model:Pt=PN(1+(fN-fst)×LD) determine the load active-power P of t momentt
Wherein, PNFor load rated active power;fNFor system nominal frequency;fstFor the system frequency value of t moment;LDFor Load coefficient of frequency modulation.
Specifically, it is described according to pre-set fuel reference signal and pitch aperture, using pre-set main vapour pressure Power model determines air inlet magnitude, including:
According to pre-set fuel reference signal RrefWith a main vapour pressure initial value GZ0, obtained by pi controller To pi controller output signal RrefPI
According to pi controller output signal RrefPI, using formula:Determine the flow L of main vapour pressure model;Its In, TwFor the time constant of boiler oil release;S is the signal that frequency is indicated in Fourier transformation;L0For main vapour pressure model Flow initial value;TDFor drum accumulation of heat volume time constant;U is superheater and live steam piping discharge coefficient;
Pitch aperture is obtained from pre-set prime mover and investigation model, according to the flow L of the main vapour pressure model, Using formula:Determine air inlet magnitude GJ;Wherein, GJ0For the first of air inlet magnitude Value;TshFor superheater volume time constant;MDFor the pitch aperture.
Specifically, according to the air inlet magnitude, the machinery of steam turbine model is determined using pre-set steam turbine model Power, including:
According to the air inlet magnitude GJ, using pre-set steam turbine model:Determine steam turbine model Mechanical output Pm;Wherein, TchFor vapor volume time constant;FHPFor high pressure cylinder power proportions;FIPFor intermediate pressure cylinder power proportions; TrhFor reheater time constant;FLPFor low pressure (LP) cylinder power proportions;TcoCross pipe time constant.
Specifically, according to the load active power and the fired power generating unit output power, equation is waved using synchronous machine It determines the system frequency predicted value of the subsequent time t+1 of current time t, and the system frequency predicted value of subsequent time t+1 is made For the system frequency value at t+1 moment, including:
According to the load active-power P of t momenttWith the fired power generating unit output power, equation is waved using synchronous machine:The system for determining the subsequent time t+1 of current time t Frequency predication value fout,t+1;Wherein, JhFor fired power generating unit inertia coeffeicent;PhFor the fired power generating unit output power;KfFor thermal motor Group coefficient of frequency modulation;
Using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment.
A kind of simulation system for reappearing power system frequency characteristic, including:
Load model unit is determined for the system frequency value according to current time t using pre-set load model Load active power;
Main vapour pressure model unit, for according to pre-set fuel reference signal and pitch aperture, using setting in advance The main vapour pressure model set determines air inlet magnitude;
Steam turbine model unit, for according to the air inlet magnitude, steamer to be determined using pre-set steam turbine model The mechanical output of machine model;
Generating set model unit, for according to the mechanical output, being generated using pre-set generating set model Electromagnetic power;
Fired power generating unit output power generation unit is used for according to the electromagnetic power, by pre-set synchronous machine member Part generates fired power generating unit output power;
Power system frequency properties determining unit, for according to the load active power and the fired power generating unit output work Rate waves the system frequency predicted value that equation determines the subsequent time t+1 of current time t using synchronous machine, and by subsequent time t System frequency value of+1 system frequency predicted value as the t+1 moment;Electric power is determined according to the system frequency predicted value at each moment System response.
In addition, the load model unit, is specifically used for:
System frequency value of the system frequency predicted value of t moment as t moment is obtained, or t moment is obtained by sampling System frequency value of the system frequency actual value as t moment;
Using pre-set load model:Pt=PN(1+(fN-fst)×LD) determine the load active-power P of t momentt
Wherein, PNFor load rated active power;fNFor system nominal frequency;fstFor the system frequency value of t moment;LDFor Load coefficient of frequency modulation.
In addition, the main vapour pressure model unit, is specifically used for:
According to pre-set fuel reference signal RrefWith a main vapour pressure initial value GZ0, obtained by pi controller To pi controller output signal RrefPI
According to pi controller output signal RrefPI, using formula:Determine the flow L of main vapour pressure model;Its In, TwFor the time constant of boiler oil release;S is the signal that frequency is indicated in Fourier transformation;L0For main vapour pressure model Flow initial value;TDFor drum accumulation of heat volume time constant;U is superheater and live steam piping discharge coefficient;
Pitch aperture is obtained from pre-set prime mover and investigation model, according to the flow L of the main vapour pressure model, Using formula:Determine air inlet magnitude GJ;Wherein, GJ0For the first of air inlet magnitude Value;TshFor superheater volume time constant;MDFor the pitch aperture.
In addition, the steam turbine model unit, is specifically used for:
According to the air inlet magnitude GJ, using pre-set steam turbine model:Determine steam turbine model Mechanical output Pm;Wherein, TchFor vapor volume time constant;FHPFor high pressure cylinder power proportions;FIPFor intermediate pressure cylinder power proportions; TrhFor reheater time constant;FLPFor low pressure (LP) cylinder power proportions;TcoCross pipe time constant.
In addition, the power system frequency properties determining unit, is specifically used for;
According to the load active-power P of t momenttWith the fired power generating unit output power, equation is waved using synchronous machine:The system for determining the subsequent time t+1 of current time t Frequency predication value fout,t+1;Wherein, JhFor fired power generating unit inertia coeffeicent;PhFor the fired power generating unit output power;KfFor thermal motor Group coefficient of frequency modulation;
Using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment.
A kind of computer readable storage medium, is stored thereon with computer program, which realizes when being executed by processor Following steps:
According to the system frequency value of current time t, load active power is determined using pre-set load model;
According to pre-set fuel reference signal and pitch aperture, using pre-set main vapour pressure model determine into Tolerance value;
According to the air inlet magnitude, the mechanical output of steam turbine model is determined using pre-set steam turbine model;
According to the mechanical output, electromagnetic power is generated using pre-set generating set model;
According to the electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element;
According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and is determined currently The system frequency predicted value of the subsequent time t+1 of moment t, and using the system frequency predicted value of subsequent time t+1 as the t+1 moment System frequency value;
Power system frequency characteristic is determined according to the system frequency predicted value at each moment.
A kind of computer equipment, including memory, processor and the meter that is stored in storage above and can run on a processor Calculation machine program, the processor realize following steps when executing described program:
According to the system frequency value of current time t, load active power is determined using pre-set load model;
According to pre-set fuel reference signal and pitch aperture, using pre-set main vapour pressure model determine into Tolerance value;
According to the air inlet magnitude, the mechanical output of steam turbine model is determined using pre-set steam turbine model;
According to the mechanical output, electromagnetic power is generated using pre-set generating set model;
According to the electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element;
According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and is determined currently The system frequency predicted value of the subsequent time t+1 of moment t, and using the system frequency predicted value of subsequent time t+1 as the t+1 moment System frequency value;
Power system frequency characteristic is determined according to the system frequency predicted value at each moment.
Analogy method and system provided in an embodiment of the present invention for reappearing power system frequency characteristic, first basis are worked as The system frequency value of preceding moment t determines load active power using pre-set load model;Later according to pre-set Fuel reference signal and pitch aperture determine air inlet magnitude using pre-set main vapour pressure model;According to the air inflow Value, the mechanical output of steam turbine model is determined using pre-set steam turbine model;According to the mechanical output, using advance The generating set model of setting generates electromagnetic power;According to the electromagnetic power, according to the electromagnetic power, by pre-setting Synchronous machine element generate fired power generating unit output power;Generate fired power generating unit output power;According to the load active power and The fired power generating unit output power waves equation using synchronous machine and determines that the system frequency of the subsequent time t+1 of current time t is pre- Measured value, and using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment;According to the system at each moment Frequency predication value determines power system frequency characteristic.As it can be seen that the embodiment of the present invention may be implemented to current power system frequency spy Property accurate reproduction, lay a good foundation for characteristic of the research new energy technology power generation in power grid.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention without having to pay creative labor, may be used also for those of ordinary skill in the art With obtain other attached drawings according to these attached drawings.
Fig. 1 is a kind of flow chart for reappearing the analogy method of power system frequency characteristic provided in an embodiment of the present invention One;
Fig. 2 is the schematic diagram one of the simulation model system in the embodiment of the present invention;
Fig. 3 is a kind of flow chart for reappearing the analogy method of power system frequency characteristic provided in an embodiment of the present invention Two;
Fig. 4 is the logical schematic of the main vapour pressure model in the embodiment of the present invention;
Fig. 5 is the logical schematic of the steam turbine model in the embodiment of the present invention;
Fig. 6 is the schematic diagram two of the simulation model system in the embodiment of the present invention;
Fig. 7 is pair of the frequency characteristic of the power system frequency characteristic curve and real system in the embodiment of the present invention Compare schematic diagram;
Fig. 8 is the structural representation provided in an embodiment of the present invention for reappearing the simulation system of power system frequency characteristic Figure.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
As shown in Figure 1, the embodiment of the present invention provides a kind of analogy method for reappearing power system frequency characteristic, packet It includes:
Step 101, the system frequency value according to current time t, determine that load is active using pre-set load model Power.
Step 102, according to pre-set fuel reference signal and pitch aperture, using pre-set main vapour pressure mould Type determines air inlet magnitude.
Step 103, according to the air inlet magnitude, the machinery of steam turbine model is determined using pre-set steam turbine model Power.
Step 104, according to the mechanical output, electromagnetic power is generated using pre-set generating set model.
Step 105, according to the electromagnetic power, generate fired power generating unit output work by pre-set synchronous machine element Rate.
Step 106, according to the load active power and the fired power generating unit output power, equation is waved using synchronous machine It determines the system frequency predicted value of the subsequent time t+1 of current time t, and the system frequency predicted value of subsequent time t+1 is made For the system frequency value at t+1 moment.
101 can be returned to step after step 106, by the continuous iteration of step 101 to step 106, finally The system frequency predicted value at each moment can be formed, and then executes subsequent step 107.
Step 107 determines power system frequency characteristic according to the system frequency predicted value at each moment.
Analogy method provided in an embodiment of the present invention for reappearing power system frequency characteristic, first according to current time The system frequency value of t determines load active power using pre-set load model;Joined later according to pre-set fuel Signal and pitch aperture are examined, air inlet magnitude is determined using pre-set main vapour pressure model;According to the air inlet magnitude, use Pre-set steam turbine model determines the mechanical output of steam turbine model;According to the mechanical output, use is pre-set Generating set model generates electromagnetic power;According to the electromagnetic power, thermal motor is generated by pre-set synchronous machine element Group output power;According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and is determined The system frequency predicted value of the subsequent time t+1 of current time t, and using the system frequency predicted value of subsequent time t+1 as t+1 The system frequency value at moment;Power system frequency characteristic is determined according to the system frequency predicted value at each moment.As it can be seen that the present invention is real Example is applied by partly being modeled key equipment, component, the device etc. in electric system, it can be achieved that current power system frequency The accurate reproduction of rate characteristic, to study characteristic of the new energy technology power generation in power grid, such as frequency modulation measure effect, analysis system frequency Rate stability etc. is laid a good foundation.
For above-mentioned pre-set load model, main vapour pressure model, steam turbine model, generating set model and same Step machine element, as shown in Fig. 2, following simulation model system realization can be used in the embodiment of the present invention, including:Load model 21, master Steam pressure model 22, steam turbine model 23, generating set model 24 and synchronous machine element 25, in addition, the tune in above-mentioned steps 102 Door aperture can be given by a prime mover and investigation model 26, above-mentioned main vapour pressure model 22, steam turbine model 23, generator Group model 24, synchronous machine element 25 and prime mover and investigation model 26 constitute entire fired power generating unit model 20.In addition, in load The resistance and capacitor model 27 for transmission line simulation actual conditions are provided between model 21 and fired power generating unit model 20.It is right The correlation modeling method and model of prime mover and investigation model 26 and generating set model 24 in fired power generating unit model 20, with And synchronous machine element 25, correlative study is had existed in the prior art, will not repeated, and can specifically joined in embodiments of the present invention Examine following document:(military sincere, Xu Zheng, Ni Qiulong wait prime mover and governor modeling and parameter identification based on response characteristic [J] electrotechnics journals, 2012,27 (6):226-231.).For load model 21, main vapour pressure model 22, steam turbine mould The specific implementation of type 23, also see subsequent embodiment.
In order to make those skilled in the art be better understood by the present invention, a more detailed embodiment is set forth below, As shown in figure 3, the embodiment of the present invention provides a kind of analogy method for reappearing power system frequency characteristic, including:
Step 301 obtains system frequency value of the system frequency predicted value of t moment as t moment, or is obtained by sampling Take the system frequency actual value of t moment as the system frequency value of t moment.
Step 302, using pre-set load model:Pt=PN(1+(fN-fst)×LD) determine that the load of t moment has Work(power Pt
Wherein, PNFor load rated active power;fNFor system nominal frequency;fstFor the system frequency value of t moment;LDFor Load coefficient of frequency modulation.
Step 303, according to pre-set fuel reference signal RrefWith a main vapour pressure initial value GZ0, by proportional integration Controller obtains pi controller output signal RrefPI
Wherein, fuel reference signal RrefIt is generally set 1 (perunit value).
Step 304, according to pi controller output signal RrefPI, using formula:Determine the flow L of main vapour pressure model.
Wherein, TwFor the time constant of boiler oil release;S is the signal that frequency is indicated in Fourier transformation;L0For main vapour The flow initial value of pressure model;TDFor drum accumulation of heat volume time constant;U is superheater and live steam piping discharge coefficient.
Step 305 obtains pitch aperture from pre-set prime mover and investigation model, according to the main vapour pressure model Flow L, using formula:Determine air inlet magnitude GJ;Wherein, GJ0For air inflow The initial value of value;TshFor superheater volume time constant;MDFor the pitch aperture.
For above-mentioned steps 303 to step 305, which may be used realization method as shown in Figure 4, Wherein PI is pi controller.
Step 306, according to the air inlet magnitude GJ, using pre-set steam turbine model:Determine steam turbine model Mechanical output Pm
Wherein, TchFor vapor volume time constant;FHPFor high pressure cylinder power proportions;FIPFor intermediate pressure cylinder power proportions;TrhFor Reheater time constant;FLPFor low pressure (LP) cylinder power proportions;TcoCross pipe time constant.
For above-mentioned steps 306, which may be used realization method as shown in Figure 5.
Step 307, according to the mechanical output, electromagnetic power is generated using pre-set generating set model.
Step 308, according to the electromagnetic power, generate fired power generating unit output work by pre-set synchronous machine element Rate.
Step 309, the load active-power P according to t momenttWith the fired power generating unit output power, shaken using synchronous machine Put equation:Determine the subsequent time t+ of current time t 1 system frequency predicted value fout,t+1
Wherein, JhFor fired power generating unit inertia coeffeicent;PhFor the fired power generating unit output power;KfFor fired power generating unit frequency modulation system Number.
Step 310, using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment.
After step 310, it can return to and execute step 301, to which the iteration gradually Jing Guo step 301 to step 310 can With formed each moment (t, t+1, t+2 ..., t+n) system frequency predicted value.In this way when subsequently continuing to execute step 311, Have:
Step 311 determines power system frequency characteristic according to the system frequency predicted value at each moment.
Illustrate the above-mentioned steps 301 of the embodiment of the present invention to the effect of step 311 with reference to one group of detailed data:
Corresponding to above-mentioned simulation model system as shown in Figure 2, by combining 22 He of main vapour pressure model as shown in Figure 4 Steam turbine model 23 as shown in Figure 5, and combine already existing prime mover and investigation model, generating set mould in the prior art Type 24 and synchronous machine element 25, can form the specific example of simulation model system as shown in FIG. 6.
Each relevant parameter is as shown in table 1 in simulation model system wherein as shown in FIG. 6:
Table 1:
301 is as shown in Figure 7 to the finally formed power system frequency characteristic curve of step 311 through the above steps, it is seen that: The frequency characteristic of finally formed power system frequency indicatrix (simulation model curve) and real system of the embodiment of the present invention It is sufficiently close to, therefore the embodiment of the present invention can reappear power system frequency characteristic completely.
Analogy method provided in an embodiment of the present invention for reappearing power system frequency characteristic, by electric system The part such as key equipment, component, device modeled the accurate reproduction, it can be achieved that current power system response, be Characteristic of the new energy technology power generation in power grid is studied, such as frequency modulation measure effect, analysis of frequency stability have established base Plinth.
Corresponding to above-mentioned Fig. 1 and embodiment of the method shown in Fig. 3, it is used for as shown in figure 8, the embodiment of the present invention provides one kind Reappear the simulation system of power system frequency characteristic, including:
Load model unit 41, it is true using pre-set load model for the system frequency value according to current time t Constant load active power.
Main vapour pressure model unit 42 is used for according to pre-set fuel reference signal and pitch aperture, using advance The main vapour pressure model of setting determines air inlet magnitude.
Steam turbine model unit 43, for according to the air inlet magnitude, vapour to be determined using pre-set steam turbine model The mechanical output of turbine model.
Generating set model unit 44, for according to the mechanical output, being given birth to using pre-set generating set model At electromagnetic power.
Fired power generating unit output power generation unit 45 is used for according to the electromagnetic power, by pre-set synchronous machine Element generates fired power generating unit output power.
Power system frequency properties determining unit 46, for being exported according to the load active power and the fired power generating unit Power waves the system frequency predicted value that equation determines the subsequent time t+1 of current time t using synchronous machine, and by lower a period of time Carve system frequency value of the system frequency predicted value of t+1 as the t+1 moment;Electricity is determined according to the system frequency predicted value at each moment Force system frequency characteristic.
In addition, the load model unit 41, is specifically used for:
System frequency value of the system frequency predicted value of t moment as t moment is obtained, or t moment is obtained by sampling System frequency value of the system frequency actual value as t moment.
Using pre-set load model:Pt=PN(1+(fN-fst)×LD) determine the load active-power P of t momentt
Wherein, PNFor load rated active power;fNFor system nominal frequency;fstFor the system frequency value of t moment;LDFor Load coefficient of frequency modulation.
In addition, the main vapour pressure model unit 42, is specifically used for:
According to pre-set fuel reference signal RrefWith a main vapour pressure initial value GZ0, obtained by pi controller To pi controller output signal RrefPI
According to pi controller output signal RrefPI, using formula:Determine the flow L of main vapour pressure model;Its In, TwFor the time constant of boiler oil release;S is the signal that frequency is indicated in Fourier transformation;L0For main vapour pressure model Flow initial value;TDFor drum accumulation of heat volume time constant;U is superheater and live steam piping discharge coefficient.
Pitch aperture is obtained from pre-set prime mover and investigation model, according to the flow L of the main vapour pressure model, Using formula:Determine air inlet magnitude GJ;Wherein, GJ0For the first of air inlet magnitude Value;TshFor superheater volume time constant;MDFor the pitch aperture.
In addition, the steam turbine model unit 43, is specifically used for:
According to the air inlet magnitude GJ, using pre-set steam turbine model:Determine steam turbine model Mechanical output Pm;Wherein, TchFor vapor volume time constant;FHPFor high pressure cylinder power proportions;FIPFor intermediate pressure cylinder power proportions; TrhFor reheater time constant;FLPFor low pressure (LP) cylinder power proportions;TcoCross pipe time constant.
In addition, the power system frequency properties determining unit 46, is specifically used for;
According to the load active-power P of t momenttWith the fired power generating unit output power, equation is waved using synchronous machine:The system for determining the subsequent time t+1 of current time t Frequency predication value fout,t+1;Wherein, JhFor fired power generating unit inertia coeffeicent;PhFor the fired power generating unit output power;KfFor thermal motor Group coefficient of frequency modulation.
Using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment.
Simulation system provided in an embodiment of the present invention for reappearing power system frequency characteristic, by electric system The part such as key equipment, component, device modeled the accurate reproduction, it can be achieved that current power system response, be Characteristic of the new energy technology power generation in power grid is studied, such as frequency modulation measure effect, analysis of frequency stability have established base Plinth.
In addition, the embodiment of the present invention also provides a kind of computer readable storage medium, it is stored thereon with computer program, it should Following steps are realized when program is executed by processor:
According to the system frequency value of current time t, load active power is determined using pre-set load model.
According to pre-set fuel reference signal and pitch aperture, using pre-set main vapour pressure model determine into Tolerance value.
According to the air inlet magnitude, the mechanical output of steam turbine model is determined using pre-set steam turbine model.
According to the mechanical output, electromagnetic power is generated using pre-set generating set model.
According to the electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element.
According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and is determined currently The system frequency predicted value of the subsequent time t+1 of moment t, and using the system frequency predicted value of subsequent time t+1 as the t+1 moment System frequency value.
Power system frequency characteristic is determined according to the system frequency predicted value at each moment.
In addition, the embodiment of the present invention also provides a kind of computer equipment, including memory, processor and it is stored in storage And the computer program that can be run on a processor, the processor realize following steps when executing described program:
According to the system frequency value of current time t, load active power is determined using pre-set load model.
According to pre-set fuel reference signal and pitch aperture, using pre-set main vapour pressure model determine into Tolerance value.
According to the air inlet magnitude, the mechanical output of steam turbine model is determined using pre-set steam turbine model.
According to the mechanical output, electromagnetic power is generated using pre-set generating set model.
According to the electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element.
According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and is determined currently The system frequency predicted value of the subsequent time t+1 of moment t, and using the system frequency predicted value of subsequent time t+1 as the t+1 moment System frequency value.
Power system frequency characteristic is determined according to the system frequency predicted value at each moment.
It should be understood by those skilled in the art that, the embodiment of the present invention can be provided as method, system or computer program Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention Apply the form of example.Moreover, the present invention can be used in one or more wherein include computer usable program code computer The computer program production implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) The form of product.
The present invention be with reference to according to the method for the embodiment of the present invention, the flow of equipment (system) and computer program product Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions every first-class in flowchart and/or the block diagram The combination of flow and/or box in journey and/or box and flowchart and/or the block diagram.These computer programs can be provided Instruct the processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce A raw machine so that the instruction executed by computer or the processor of other programmable data processing devices is generated for real The device for the function of being specified in present one flow of flow chart or one box of multiple flows and/or block diagram or multiple boxes.
These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works so that instruction generation stored in the computer readable memory includes referring to Enable the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device so that count Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, in computer or The instruction executed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one The step of function of being specified in a box or multiple boxes.
Specific embodiment is applied in the present invention, and principle and implementation of the present invention are described, above example Explanation be merely used to help understand the present invention method and its core concept;Meanwhile for those of ordinary skill in the art, According to the thought of the present invention, there will be changes in the specific implementation manner and application range, in conclusion in this specification Appearance should not be construed as limiting the invention.

Claims (12)

1. a kind of analogy method for reappearing power system frequency characteristic, which is characterized in that including:
According to the system frequency value of current time t, load active power is determined using pre-set load model;
According to pre-set fuel reference signal and pitch aperture, air inflow is determined using pre-set main vapour pressure model Value;
According to the air inlet magnitude, the mechanical output of steam turbine model is determined using pre-set steam turbine model;
According to the mechanical output, electromagnetic power is generated using pre-set generating set model;
According to the electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element;
According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and determines current time The system frequency predicted value of the subsequent time t+1 of t, and be as the t+1 moment using the system frequency predicted value of subsequent time t+1 System frequency values;
Power system frequency characteristic is determined according to the system frequency predicted value at each moment.
2. the analogy method according to claim 1 for reappearing power system frequency characteristic, which is characterized in that described According to the system frequency value of current time t, load active power is determined using pre-set load model, including:
System frequency value of the system frequency predicted value of t moment as t moment is obtained, or is by what sampling obtained t moment System frequency value of the frequency actual value of uniting as t moment;
Using pre-set load model:Pt=PN(1+(fN-fst)×LD) determine the load active-power P of t momentt
Wherein, PNFor load rated active power;fNFor system nominal frequency;fstFor the system frequency value of t moment;LDFor load Coefficient of frequency modulation.
3. the analogy method according to claim 2 for reappearing power system frequency characteristic, which is characterized in that described According to pre-set fuel reference signal and pitch aperture, air inlet magnitude is determined using pre-set main vapour pressure model, is wrapped It includes:
According to pre-set fuel reference signal RrefWith a main vapour pressure initial value GZ0, compared by pi controller Example integral controller output signal RrefPI
According to pi controller output signal RrefPI, using formula:Determine the flow L of main vapour pressure model;Its In, TwFor the time constant of boiler oil release;S is the signal that frequency is indicated in Fourier transformation;L0For main vapour pressure model Flow initial value;TDFor drum accumulation of heat volume time constant;U is superheater and live steam piping discharge coefficient;
Pitch aperture is obtained from pre-set prime mover and investigation model, according to the flow L of the main vapour pressure model, is used Formula:Determine air inlet magnitude GJ;Wherein, GJ0For the initial value of air inlet magnitude;Tsh For superheater volume time constant;MDFor the pitch aperture.
4. the analogy method according to claim 3 for reappearing power system frequency characteristic, which is characterized in that according to institute Air inlet magnitude is stated, the mechanical output of steam turbine model is determined using pre-set steam turbine model, including:
According to the air inlet magnitude GJ, using pre-set steam turbine model:Determine steam turbine model Mechanical output Pm;Wherein, TchFor vapor volume time constant;FHPFor high pressure cylinder power proportions;FIPFor intermediate pressure cylinder power proportions; TrhFor reheater time constant;FLPFor low pressure (LP) cylinder power proportions;TcoCross pipe time constant.
5. the analogy method according to claim 4 for reappearing power system frequency characteristic, which is characterized in that according to institute Load active power and the fired power generating unit output power are stated, lower a period of time that equation determines current time t is waved using synchronous machine The system frequency predicted value of t+1 is carved, and using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment, Including:
According to the load active-power P of t momenttWith the fired power generating unit output power, equation is waved using synchronous machine:The system for determining the subsequent time t+1 of current time t Frequency predication value fout,t+1;Wherein, JhFor fired power generating unit inertia coeffeicent;PhFor the fired power generating unit output power;KfFor thermal motor Group coefficient of frequency modulation;
Using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment.
6. a kind of simulation system for reappearing power system frequency characteristic, which is characterized in that including:
Load model unit determines load for the system frequency value according to current time t using pre-set load model Active power;
Main vapour pressure model unit, for pre-set according to pre-set fuel reference signal and pitch aperture, use Main vapour pressure model determines air inlet magnitude;
Steam turbine model unit, for according to the air inlet magnitude, steam turbine mould to be determined using pre-set steam turbine model The mechanical output of type;
Generating set model unit, for according to the mechanical output, electromagnetism to be generated using pre-set generating set model Power;
Fired power generating unit output power generation unit, for according to the electromagnetic power, being given birth to by pre-set synchronous machine element At fired power generating unit output power;
Power system frequency properties determining unit is used for according to the load active power and the fired power generating unit output power, The system frequency predicted value that equation determines the subsequent time t+1 of current time t is waved using synchronous machine, and by subsequent time t+1 System frequency value of the system frequency predicted value as the t+1 moment;Power train is determined according to the system frequency predicted value at each moment System frequency characteristic.
7. the simulation system according to claim 6 for reappearing power system frequency characteristic, which is characterized in that described negative Lotus model unit, is specifically used for:
System frequency value of the system frequency predicted value of t moment as t moment is obtained, or is by what sampling obtained t moment System frequency value of the frequency actual value of uniting as t moment;
Using pre-set load model:Pt=PN(1+(fN-fst)×LD) determine the load active-power P of t momentt
Wherein, PNFor load rated active power;fNFor system nominal frequency;fstFor the system frequency value of t moment;LDFor load Coefficient of frequency modulation.
8. the simulation system according to claim 7 for reappearing power system frequency characteristic, which is characterized in that the master Steam pressure model unit, is specifically used for:
According to pre-set fuel reference signal RrefWith a main vapour pressure initial value GZ0, compared by pi controller Example integral controller output signal RrefPI
According to pi controller output signal RrefPI, using formula:Determine the flow L of main vapour pressure model;Its In, TwFor the time constant of boiler oil release;S is the signal that frequency is indicated in Fourier transformation;L0For main vapour pressure model Flow initial value;TDFor drum accumulation of heat volume time constant;U is superheater and live steam piping discharge coefficient;
Pitch aperture is obtained from pre-set prime mover and investigation model, according to the flow L of the main vapour pressure model, is used Formula:Determine air inlet magnitude GJ;Wherein, GJ0For the initial value of air inlet magnitude;Tsh For superheater volume time constant;MDFor the pitch aperture.
9. the simulation system according to claim 8 for reappearing power system frequency characteristic, which is characterized in that the vapour Turbine model unit, is specifically used for:
According to the air inlet magnitude GJ, using pre-set steam turbine model:Determine steam turbine model Mechanical output Pm;Wherein, TchFor vapor volume time constant;FHPFor high pressure cylinder power proportions;FIPFor intermediate pressure cylinder power proportions; TrhFor reheater time constant;FLPFor low pressure (LP) cylinder power proportions;TcoCross pipe time constant.
10. the simulation system according to claim 9 for reappearing power system frequency characteristic, which is characterized in that described Power system frequency properties determining unit, is specifically used for;
According to the load active-power P of t momenttWith the fired power generating unit output power, equation is waved using synchronous machine:The system for determining the subsequent time t+1 of current time t Frequency predication value fout,t+1;Wherein, JhFor fired power generating unit inertia coeffeicent;PhFor the fired power generating unit output power;KfFor thermal motor Group coefficient of frequency modulation;
Using the system frequency predicted value of subsequent time t+1 as the system frequency value at t+1 moment.
11. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is by processor Following steps are realized when execution:
According to the system frequency value of current time t, load active power is determined using pre-set load model;
According to pre-set fuel reference signal and pitch aperture, air inflow is determined using pre-set main vapour pressure model Value;
According to the air inlet magnitude, the mechanical output of steam turbine model is determined using pre-set steam turbine model;
According to the mechanical output, electromagnetic power is generated using pre-set generating set model;
According to the electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element;
According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and determines current time The system frequency predicted value of the subsequent time t+1 of t, and be as the t+1 moment using the system frequency predicted value of subsequent time t+1 System frequency values;
Power system frequency characteristic is determined according to the system frequency predicted value at each moment.
12. a kind of computer equipment, including memory, processor and the calculating that is stored in storage above and can run on a processor Machine program, which is characterized in that the processor realizes following steps when executing described program:
According to the system frequency value of current time t, load active power is determined using pre-set load model;
According to pre-set fuel reference signal and pitch aperture, air inflow is determined using pre-set main vapour pressure model Value;
According to the air inlet magnitude, the mechanical output of steam turbine model is determined using pre-set steam turbine model;
According to the mechanical output, electromagnetic power is generated using pre-set generating set model;
According to the electromagnetic power, fired power generating unit output power is generated by pre-set synchronous machine element;
According to the load active power and the fired power generating unit output power, equation is waved using synchronous machine and determines current time The system frequency predicted value of the subsequent time t+1 of t, and be as the t+1 moment using the system frequency predicted value of subsequent time t+1 System frequency values;
Power system frequency characteristic is determined according to the system frequency predicted value at each moment.
CN201810361848.5A 2018-04-20 2018-04-20 Simulation method and system for reproducing frequency characteristics of power system Active CN108565854B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810361848.5A CN108565854B (en) 2018-04-20 2018-04-20 Simulation method and system for reproducing frequency characteristics of power system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810361848.5A CN108565854B (en) 2018-04-20 2018-04-20 Simulation method and system for reproducing frequency characteristics of power system

Publications (2)

Publication Number Publication Date
CN108565854A true CN108565854A (en) 2018-09-21
CN108565854B CN108565854B (en) 2024-01-23

Family

ID=63536082

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810361848.5A Active CN108565854B (en) 2018-04-20 2018-04-20 Simulation method and system for reproducing frequency characteristics of power system

Country Status (1)

Country Link
CN (1) CN108565854B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020199644A1 (en) * 2019-04-04 2020-10-08 华北电力科学研究院有限责任公司 Prime mover simulation method and apparatus, and multi-machine synchronization system frequency characteristic simulation system

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU903574A1 (en) * 1980-05-22 1982-02-07 Белорусский Ордена Трудового Красного Знамени Политехнический Институт Power-generating unit load control system
WO2006097495A2 (en) * 2005-03-18 2006-09-21 Siemens Aktiengesellschaft Method and device for production of controlled power by means of a combined gas and steam turbine plant
CN101446807A (en) * 2008-08-15 2009-06-03 中国电力科学研究院 Realization method for heat-engine plant speed regulating system model in power system simulation
CN102946224A (en) * 2012-12-10 2013-02-27 贵州电网公司电网规划研究中心 Thermal power unit power prediction control system for ensuring safety and stability of sending-end electric grid
JP2013174223A (en) * 2012-02-27 2013-09-05 Toshiba Corp Speed governing controller for steam turbine, method for controlling the same and steam turbine
CN103336911A (en) * 2013-07-25 2013-10-02 哈尔滨工业大学 Mechanism-based mathematic model establishing method suitable for frequency modulation analysis of multi-machine power system
CN103670536A (en) * 2013-05-30 2014-03-26 甘肃大唐国际连城发电有限责任公司 Adjustment method for steam turbine control valve flows in thermal power plant
CN104794297A (en) * 2015-04-29 2015-07-22 国网上海市电力公司电力科学研究院电源技术中心 Method for eliminating distortion in steam turbine speed regulation system model
CN105275509A (en) * 2015-11-25 2016-01-27 广东电网有限责任公司电力科学研究院 Method and system for controlling opening degree of steam turbine governing valve of thermal power generating unit
CN105870943A (en) * 2016-05-18 2016-08-17 国网山东省电力公司电力科学研究院 Primary frequency modulation control system and method applied to thermal power generating unit and based on safe and stable operation
CN107341559A (en) * 2016-11-25 2017-11-10 国网山东省电力公司电力科学研究院 A kind of Primary frequency control ability Forecasting Methodology based on main vapour pressure tolerance
CN107482649A (en) * 2017-07-13 2017-12-15 上海电力学院 A kind of two domain interacted system LOAD FREQUENCY control methods based on frequency dividing control

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU903574A1 (en) * 1980-05-22 1982-02-07 Белорусский Ордена Трудового Красного Знамени Политехнический Институт Power-generating unit load control system
WO2006097495A2 (en) * 2005-03-18 2006-09-21 Siemens Aktiengesellschaft Method and device for production of controlled power by means of a combined gas and steam turbine plant
CN101446807A (en) * 2008-08-15 2009-06-03 中国电力科学研究院 Realization method for heat-engine plant speed regulating system model in power system simulation
JP2013174223A (en) * 2012-02-27 2013-09-05 Toshiba Corp Speed governing controller for steam turbine, method for controlling the same and steam turbine
CN102946224A (en) * 2012-12-10 2013-02-27 贵州电网公司电网规划研究中心 Thermal power unit power prediction control system for ensuring safety and stability of sending-end electric grid
CN103670536A (en) * 2013-05-30 2014-03-26 甘肃大唐国际连城发电有限责任公司 Adjustment method for steam turbine control valve flows in thermal power plant
CN103336911A (en) * 2013-07-25 2013-10-02 哈尔滨工业大学 Mechanism-based mathematic model establishing method suitable for frequency modulation analysis of multi-machine power system
CN104794297A (en) * 2015-04-29 2015-07-22 国网上海市电力公司电力科学研究院电源技术中心 Method for eliminating distortion in steam turbine speed regulation system model
CN105275509A (en) * 2015-11-25 2016-01-27 广东电网有限责任公司电力科学研究院 Method and system for controlling opening degree of steam turbine governing valve of thermal power generating unit
CN105870943A (en) * 2016-05-18 2016-08-17 国网山东省电力公司电力科学研究院 Primary frequency modulation control system and method applied to thermal power generating unit and based on safe and stable operation
CN107341559A (en) * 2016-11-25 2017-11-10 国网山东省电力公司电力科学研究院 A kind of Primary frequency control ability Forecasting Methodology based on main vapour pressure tolerance
CN107482649A (en) * 2017-07-13 2017-12-15 上海电力学院 A kind of two domain interacted system LOAD FREQUENCY control methods based on frequency dividing control

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
EN LU等: "Black-start strategy for power grids including fast cut thermal power units", 《2013 IEEE POWER & ENERGY SOCIETY GENERAL MEETING》 *
梅柏杉;张德兰;陈;: "火电机组一次调频的仿真研究", 电力科学与工程, no. 11 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020199644A1 (en) * 2019-04-04 2020-10-08 华北电力科学研究院有限责任公司 Prime mover simulation method and apparatus, and multi-machine synchronization system frequency characteristic simulation system

Also Published As

Publication number Publication date
CN108565854B (en) 2024-01-23

Similar Documents

Publication Publication Date Title
Liu et al. A method for sizing energy storage system to increase wind penetration as limited by grid frequency deviations
CN101969199B (en) Fault loss estimation method for risk assessment of transient power angle stability
CN106849088B (en) A kind of mains frequency characteristic computing method based on the control of wind-powered electricity generation pitch primary frequency modulation
CN102346442B (en) Method for carrying out parallel transient simulation computing on microgrid suitable for nodal analysis frameworks
Zhang et al. Synchrophasor measurement-based wind plant inertia estimation
CN110750882A (en) Wind power ratio limit value analytical calculation method considering frequency constraint
CN109449937B (en) Method for determining new energy bearing capacity of power grid through system frequency stability constraint
CN105303032A (en) Analysis method for objective factor affecting generator set energy efficiency
CN105162164A (en) Method of building low-order dynamic frequency response model with wind power integration system
CN109657309B (en) Method and device for simplifying and calculating long-process frequency response of power system
Suehiro et al. Decentralized control of smart grid by using overlapping information
CN107220414A (en) A kind of nuclear power generating sets cooling agent shielding main pump dynamic modeling and analysis method
Chen et al. Online emergency control to suppress frequency oscillations based on damping evaluation using dissipation energy flow
Dong et al. Spectrum-domain stability assessment and intrinsic oscillation for aggregated mobile energy storage in grid frequency regulation
Ju et al. Analytic assessment of the power system frequency security
Mondai et al. Analysis of limiting bounds for stalling of natural gas genset in the CERTS microgrid test bed
RU2621422C2 (en) System and method for steam turbine parameter testing
Mohamed et al. Mathematical modelling for coal fired supercritical power plants and model parameter identification using genetic algorithms
Vanfretti et al. A Modelica power system component library for model validation and parameter identification
CN108565854A (en) Analogy method and system for reappearing power system frequency characteristic
CN109002741B (en) Method and system for simulating transmission power of primary and secondary loop systems of pressurized water reactor nuclear motor unit
CN114254540A (en) Simulation method and simulation system of permanent magnet synchronous generator based on digital goldenrain tree generator
El-Guindy et al. Formal LPV control for transient stability of power systems
CA2829683C (en) Method for high fidelity modeling of an aircraft electrical power system
Wang et al. Development of wind-energy modeling technology and standards

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