CN103323252B - Extraction method for power response characteristic parameter of steam turbine - Google Patents
Extraction method for power response characteristic parameter of steam turbine Download PDFInfo
- Publication number
- CN103323252B CN103323252B CN201310292709.9A CN201310292709A CN103323252B CN 103323252 B CN103323252 B CN 103323252B CN 201310292709 A CN201310292709 A CN 201310292709A CN 103323252 B CN103323252 B CN 103323252B
- Authority
- CN
- China
- Prior art keywords
- power
- time
- pressure cylinder
- high pressure
- steam turbine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000001105 regulatory effect Effects 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000009499 grossing Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Landscapes
- Control Of Turbines (AREA)
Abstract
The invention discloses an extraction method for a power response characteristic parameter of a steam turbine. The extraction method comprises the following steps of: obtaining a power change curve of a unit according to parameter actual measurement tests of the steam turbine and a regulating system of the steam turbine, and confirming a starting time T0 when the power of the unit of the steam turbine changes and an initial power P0 of the unit; taking data of the power change curve of the unit in a time period from the starting time T0 to T0+deltaT as effective data to extract a maximum power output power PHP of a high pressure cylinder of the steam turbine and a peak time THP of the high pressure cylinder; carrying out polynomial fit on the extracted effective data, so as to confirm the PHP and the THP. According to the method disclosed by the invention, the maximum power output power PHP of the high pressure cylinder of the steam turbine and the peak time THP of the high pressure cylinder can be accurately extracted, and the problem of relatively low accuracies of the obtained PHP and THP caused by strong manual examination subjectivity and noise interference is avoided.
Description
Technical field
The present invention relates to steam turbine and regulating system actual measurement with modeling field, particularly a kind of extracting method of steam turbine power response characteristic parameter.
Background technology
Steam turbine and regulating system actual measurement thereof are the important component parts that electric system " four large parameters are surveyed " works with modeling.For a long time, electric system situation about not being inconsistent in the characteristic that the steam turbine control system model carrying out using in system stability calculating exists model parameter and actual modeling object.The model of distortion will cause the deviation even mistake of Power System Stability Analysis result.The method establishment that steam turbine and regulating system actual measurement thereof and the object of modeling are by actual tests truly can reflect the realistic model of actual steam turbine and regulating system characteristic thereof.
" synchronous generator prime mover and regulating system parameter thereof are surveyed and modeling directive/guide " the i.e. Q/GDW748-2012 that industry is promulgated requires that the model set up should confirm its correctness through emulating the mode of checking.In emulation check process, the high pressure cylinder peak power of actual set power response curve is exerted oneself power P
hPwith high pressure cylinder time to peak T
hPthe calculating of checking index is used to as important Q factor.At present to P
hPand T
hPextraction mainly rely on testing crew to carry out judging and distinguish according to the variation tendency of actual set power response curve and obtain.This mode subjectivity is strong and noise in process of the test will add big error further, thus causes P
hPand T
hPprecision not high, the final confidence level reducing check result.
Summary of the invention
Technical matters to be solved by this invention is, not enough for prior art, provides a kind of extracting method of steam turbine power response characteristic parameter, to exert oneself power P to high pressure cylinder peak power
hPwith high pressure cylinder time to peak T
hPaccurately extract, avoid that artificial subjectivity of checking is strong, noise causes gained P
hPand T
hPthe lower problem of precision.
For solving the problems of the technologies described above, the technical solution adopted in the present invention is: a kind of extracting method of steam turbine power response characteristic parameter, and it is characterized in that, the method is:
1) data prediction: determine the initial time T that steam-turbine unit power changes
0with the initial power P of unit
0;
2) initial time T is got
0to T
0the data of the power of the assembling unit change curve in this time period of+△ T to be exerted oneself power P as extraction steam turbine high-pressure cylinder peak power
hPwith high pressure cylinder time to peak T
hPvalid data; Described power of the assembling unit change curve obtains according to steam turbine and regulating system parameter measured test thereof;
3) carry out fitting of a polynomial to the valid data of said extracted, obtaining fitting of a polynomial result is: P=f (t), and wherein t is time variable, and P represents the fitting power of change of corresponding time corresponding to t;
4) to P=f (t) differentiate, the derivative equation of fitting of a polynomial result is obtained
solve described derivative equation, obtain all solutions of described derivative equation;
5) high pressure cylinder peak power time T is determined
1: when described derivative equation a certain solution is greater than T
0+ △ T
1(△ T
1be a time regulatable interval, its value can adjust between 0.1s ~ 0.2s) minimum real solution time, then this root solution is high pressure cylinder peak power time T
1value;
6) high pressure cylinder peak power P is determined
1: according to T obtained above
1, according to time and power, the corresponding relation on power of the assembling unit change curve determines T
1time the power of the assembling unit value of inscribing, this value is high pressure cylinder peak power P
1numerical value;
7) calculate high pressure cylinder peak power to exert oneself power P
hP: according to described high pressure cylinder peak power P
1with initial power P
0calculate P
hP, P
hP=P
1-P
0;
8) high pressure cylinder time to peak T is calculated
hP: according to described high pressure cylinder peak power time T
1with initial time T
0calculate T
hP, T
hP=T
1-T
0.
Compared with prior art, the beneficial effect that the present invention has is: method of the present invention can be exerted oneself power P to high pressure cylinder peak power
hPwith high pressure cylinder time to peak T
hPaccurately extract, avoid that artificial subjectivity of checking is strong, noise causes gained P
hPand T
hPthe lower problem of precision; The computational accuracy of the inventive method is high, thus ensure that the accuracy and confidence of follow-up check.
Accompanying drawing explanation
Fig. 1 is one embodiment of the invention method flow diagram;
Fig. 2 is power of the assembling unit change curve;
The changing trend diagram of the frequency departure that Fig. 3 obtains according to steam turbine and regulating system parameter measured test thereof for one embodiment of the invention;
Fig. 4 is the power of the assembling unit change curve that one embodiment of the invention obtains according to steam turbine and regulating system parameter measured test thereof;
Fig. 5 is the matched curve figure that the present invention is based on embodiment polynomial fitting;
Fig. 6 is one embodiment of the invention time and the corresponding relation figure of power on power of the assembling unit change curve.
Embodiment
As shown in Figure 1, the extracting method of one embodiment of the invention steam turbine power response characteristic parameter comprises the steps:
(1) data prediction:
(1.1) the initial time T that the power of the assembling unit changes is determined
0: according to the variation tendency determination initial time T of the frequency departure that steam turbine and regulating system parameter measured test thereof obtain
0, the time corresponding to frequency departure generation Spline smoothing is initial time T
0.
(1.2) the initial power P of unit is determined
0: power of the assembling unit change curve (as Fig. 2) can be obtained according to steam turbine and regulating system parameter measured test thereof, by power of the assembling unit change curve determination initial power P
0.Get initial time T in (1.1)
0on power of the assembling unit change curve before, the mean value of all data is as initial power P
0.
(2) valid data are selected: get initial time T
0to T
0the data of the power of the assembling unit change curve in this time period of+△ T are as extraction P
hPand T
hPvalid data.△ T is an adjustable time interval, and when △ T is too little, then the valid data comprised in selected data very little thus cannot to P
hPand T
hPeffectively extract; When △ T is excessive, then the invalid data comprised in selected data too much thus increase and extract P
hPand T
hPcalculated amount and error.From actual effect, it is comparatively suitable that △ T was taken between 1 ~ 3 second.
(3) fitting of a polynomial and derivative equation rooting thereof:
(3.1) fitting of a polynomial: according to the valid data selected by (2) joint, fitting of a polynomial is carried out to it.Polynomial order is selected between 8 ~ 12, when order less or larger time, all may reduce the fitting precision of fitting of a polynomial.T
0to T
0the curve of the polynomial fitting in this time period of+△ T accurately can reflect the trend of same time period internal power response change curve, this curve is a level and smooth curve thus eliminates the impact of noise simultaneously, the fitting of a polynomial result finally obtained is following form: P=f (t), in formula, t is time variable, and P represents the fitting power corresponding to corresponding time variable t.
(3.2) derivative equation rooting: get fitting of a polynomial result P=f (t) that (3.1) obtain, first carries out differentiate to it, namely has
thus obtain the derivative equation of fitting of a polynomial result
derivative equation is solved to all solutions finally obtaining equation.
(4) high pressure cylinder peak power time T is determined
1with high pressure cylinder peak power P
1
(4.1) high pressure cylinder peak power time T is determined
1: the high pressure cylinder peak power time be first sharp wave on power of the assembling unit change curve in (1.2) crest corresponding to time.Mathematically understand, this point is the Local Extremum of power of the assembling unit change curve, therefore T
1one of them that in (3.2), all of derivative equation is separated.All of derivative equation solution in search (3.2), when the root of one of them root is greater than T
0+ △ T
1(△ T
1be a time regulatable interval, its value can adjust between 0.1s ~ 0.2s) real number and minimum time, then the root of this root is high pressure cylinder peak power time T
1numerical value.
(4.2) high pressure cylinder peak power P is determined
1: according to obtained T
1, according to time and power, the corresponding relation on power of the assembling unit change curve determines T
1time the power of the assembling unit value of inscribing, this value is high pressure cylinder peak power P
1numerical value.
(5) calculate high pressure cylinder peak power to exert oneself power P
hPwith high pressure cylinder time to peak T
hP:
(5.1) calculate high pressure cylinder peak power to exert oneself power P
hP: the high pressure cylinder peak power P determined according to (4.2)
1(1.2) determined initial power P
0calculate P
hP, P
hP=P
1-P
0.
(5.2) high pressure cylinder time to peak T is calculated
hP: according to (4.1) determined high pressure cylinder peak power time T
1(1.1) determined initial time T
0calculate T
hP, T
hP=T
1-T
0.
Method of the present invention is described in detail below for a certain steam-turbine unit:
(1) data prediction:
(1.1) the initial time T that the power of the assembling unit changes is determined
0: the variation tendency (as shown in Figure 3) of the frequency departure obtained according to steam turbine and regulating system parameter measured test thereof determines initial time T
0, the time corresponding to frequency departure generation Spline smoothing is initial time T
0, now T
0=21.45s.
(1.2) the initial power P of unit is determined
0: power of the assembling unit change curve (as shown in Figure 4) can be obtained according to steam turbine and regulating system parameter measured test thereof, by power of the assembling unit change curve determination initial power P
0.Get initial time T
0(T
0=21.45s) before power of the assembling unit change curve on the mean value of all data as initial power P
0, now P
0=496.74MW.
(2) valid data are selected: make △ T=1s, get initial time T
0to T
0the data of the power of the assembling unit change curve (as shown in Figure 2) in this time period of+△ T are as extraction P
hPand T
hPvalid data.
(3) fitting of a polynomial and derivative equation rooting thereof:
(3.1) fitting of a polynomial: according to the valid data (as shown in Figure 3) selected by (2) joint, fitting of a polynomial is carried out to it.Polynomial order is chosen as 8, and the fitting of a polynomial result finally obtained is following form:
P=f(t)=-8.8994t
8+849.6632t
7-1.0329e4t
6-1.9652e6t
5+1.2043e8t
4-3.2852e9t
3
4.8921e10t
2-3.8760e11t+1.2854e12, in formula, t is time variable, and P represents the fitting power corresponding to corresponding time variable t, based on this polynomial fitting matched curve as shown in Figure 5.
(3.2) derivative equation rooting: get fitting of a polynomial result P=f (t) that (3.1) obtain, first carries out differentiate to it, namely has
thus obtain the derivative equation of fitting of a polynomial result
derivative equation is solved to all solutions finally obtaining equation, its root solution set=-48.3464,22.767,22.4861,22.1719,21.4859 ± 0.1309i, 21.4894}.
(4) high pressure cylinder peak power time T is determined
1with high pressure cylinder peak power P
1
(4.1) high pressure cylinder peak power time T is determined
1: the high pressure cylinder peak power time be first sharp wave on power of the assembling unit change curve in (1.2) crest corresponding to time.Mathematically understand, this point is the Local Extremum of power of the assembling unit change curve, therefore T
1one of them that in (3.2), all of derivative equation is separated.All of derivative equation solution in search (3.2), when the root of one of them root is greater than T
0+ △ T
1(make △ T
1=0.2s, △ T
1be a time regulatable interval, its value can adjust between 0.1s ~ 0.2s) real number and minimum time, then the root of this root is high pressure cylinder peak power time T
1numerical value, now T
1=22.1719s.
(4.2) high pressure cylinder peak power P is determined
1: according to obtained T
1, according to time and power, the corresponding relation on power of the assembling unit change curve (as shown in Figure 6) determines T
1time the power of the assembling unit value of inscribing, this value is high pressure cylinder peak power P
1numerical value, now P
1=516.7MW.
(5) calculate high pressure cylinder peak power to exert oneself power P
hPwith high pressure cylinder time to peak T
hP:
(5.1) calculate high pressure cylinder peak power to exert oneself power P
hP: the high pressure cylinder peak power P determined according to (4.2)
1(1.2) determined initial power P
0calculate P
hP, P
hP=P
1-P
0=516.7MW-496.74MW=19.96MW.
(5.2) high pressure cylinder time to peak T is calculated
hP: according to (4.1) determined high pressure cylinder peak power time T
1(1.1) determined initial time T
0calculate T
hP, T
hP=T
1-T
0=22.1719s-21.45s=0.7219s.
Claims (3)
1. an extracting method for steam turbine power response characteristic parameter, is characterized in that, the method is:
1) data prediction: determine the initial time T that steam-turbine unit power changes
0with the initial power P of unit
0; Initial time T
0deterministic process for: according to the variation tendency determination initial time T of the frequency departure that steam turbine and regulating system parameter measured test thereof obtain
0, the time corresponding to frequency departure generation Spline smoothing is initial time T
0; Initial power P
0deterministic process be: obtain power of the assembling unit change curve according to steam turbine and regulating system parameter measured test thereof, get initial time T
0on power of the assembling unit change curve before, the mean value of all data is as initial power P
0;
2) initial time T is got
0to T
0the data of the power of the assembling unit change curve in this time period of+△ T to be exerted oneself power P as extraction steam turbine high-pressure cylinder peak power
hPwith high pressure cylinder time to peak T
hPvalid data; Described power of the assembling unit change curve obtains according to steam turbine operation data;
3) carry out fitting of a polynomial to the valid data of said extracted, obtaining fitting of a polynomial result is: P=f (t), and wherein t is time variable, and P represents the fitting power corresponding to corresponding time variable t;
4) to P=f (t) differentiate, the derivative equation of fitting of a polynomial result is obtained
solve described derivative equation, obtain all solutions of described derivative equation;
5) high pressure cylinder peak power time T is determined
1: when described derivative equation a certain solution is greater than T
0+ △ T
1minimum real solution time, then this root solution is high pressure cylinder peak power time T
1value, wherein △ T
1for 0.1s ~ 0.2s;
6) high pressure cylinder peak power P is determined
1: according to T obtained above
1, according to time and power, the corresponding relation on power of the assembling unit change curve determines T
1time the power of the assembling unit value of inscribing, this value is high pressure cylinder peak power P
1numerical value;
7) calculate high pressure cylinder peak power to exert oneself power P
hP: according to described high pressure cylinder peak power P
1with initial power P
0calculate P
hP, P
hP=P
1-P
0;
8) high pressure cylinder time to peak T is calculated
hP: according to described high pressure cylinder peak power time T
1with initial time T
0calculate T
hP, T
hP=T
1-T
0.
2. the extracting method of steam turbine power response characteristic parameter according to claim 1, is characterized in that, described step 2) in, △ T value is 1 ~ 3 second.
3. the extracting method of steam turbine power response characteristic parameter according to claim 1, is characterized in that, described step 3) in, when carrying out fitting of a polynomial, polynomial order is selected between 8 ~ 12.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310292709.9A CN103323252B (en) | 2013-07-12 | 2013-07-12 | Extraction method for power response characteristic parameter of steam turbine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310292709.9A CN103323252B (en) | 2013-07-12 | 2013-07-12 | Extraction method for power response characteristic parameter of steam turbine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103323252A CN103323252A (en) | 2013-09-25 |
| CN103323252B true CN103323252B (en) | 2015-06-10 |
Family
ID=49192141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310292709.9A Active CN103323252B (en) | 2013-07-12 | 2013-07-12 | Extraction method for power response characteristic parameter of steam turbine |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103323252B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103955134B (en) * | 2014-05-05 | 2016-06-22 | 国家电网公司 | Steam turbine model parameter identification method based on power response characteristic |
| CN108573113A (en) * | 2018-05-09 | 2018-09-25 | 南方电网科学研究院有限责任公司 | Method and system for automatically identifying quality parameters of turbine frequency step response |
| CN110595795B (en) * | 2019-09-05 | 2021-02-23 | 安徽江淮汽车集团股份有限公司 | Vehicle emission comparison test method, device, equipment and computer readable storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102882207A (en) * | 2012-09-18 | 2013-01-16 | 安徽省电力公司 | Method for controlling scheduling tidal current limit of power grid |
| CN102914633A (en) * | 2012-10-16 | 2013-02-06 | 重庆大学 | Method for determining optimum water content for soil compaction |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7461544B2 (en) * | 2006-02-24 | 2008-12-09 | General Electric Company | Methods for detecting water induction in steam turbines |
-
2013
- 2013-07-12 CN CN201310292709.9A patent/CN103323252B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102882207A (en) * | 2012-09-18 | 2013-01-16 | 安徽省电力公司 | Method for controlling scheduling tidal current limit of power grid |
| CN102914633A (en) * | 2012-10-16 | 2013-02-06 | 重庆大学 | Method for determining optimum water content for soil compaction |
Non-Patent Citations (4)
| Title |
|---|
| 《300MW汽轮机调速***模型与非线性分析》;王颖等;《汽轮机技术》;20070228;第49卷(第1期);第17-20页 * |
| 《基于响应特性的原动机及调速器建模与参数辨识》;武诚等;《电工技术学报》;20120630;第27卷(第6期);第226-231页 * |
| 晁晖等.《3 术语和定义》.《Q/GDW 748-2012 同步发电机原动机及其调节***参数实测与建模导则》.2012,第3页. * |
| 高葆新等.《2.2.2目标函数的极值 2.3.2多项式拟合》.《微波电路计算机辅助设计》.1988,第61页、第64-72页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103323252A (en) | 2013-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103955134B (en) | Steam turbine model parameter identification method based on power response characteristic | |
| CN109713685B (en) | Online positioning method suitable for VSC access-induced subsynchronous oscillation | |
| CN109698521B (en) | Low-penetration characteristic identification method of photovoltaic inverter based on measured data | |
| CN103323252B (en) | Extraction method for power response characteristic parameter of steam turbine | |
| CN110516275A (en) | Simulation parameters check method based on disturbance information | |
| CN111507637B (en) | Water turbine for stable calculation of electric power system and diversion system parameter modeling and testing method | |
| CN102052233A (en) | Water turbine regulating system module used for stability analysis of power system | |
| CN107989743A (en) | A kind of Kaplan turbine water-saving and synergistic adjusting method and regulating system | |
| CN101930494A (en) | Method for identifying aircraft model with undetermined order and parameters based on mode segmentation and genetic algorithm | |
| CN103807090A (en) | A impulse turbine governing system for electric power system stability analysis | |
| CN109800455B (en) | Transient reactive characteristic simulation method and device for doubly-fed wind turbine generator | |
| CN104236893A (en) | Performance parameter test system and performance parameter test method of hydraulic damper | |
| CN105305438A (en) | New energy power station model verification method based on variable impedance and controlled alternating-current voltage source | |
| CN104134013A (en) | Wind turbine blade modal analysis method | |
| CN108767905A (en) | Primary frequency modulation Commissioning Analysis method is carried out based on synchronized phase measurement device data | |
| CN105260548A (en) | Turbine model building method based on unit practical features | |
| CN105356498A (en) | Verification system and method for low-voltage ride-through simulation model of wind power unit | |
| CN110098610B (en) | Real-time identification method and system for oscillation leading mode of power system under fault disturbance | |
| CN105467348A (en) | Data acquisition system, virtual verification system and method for fiber current transformer | |
| CN110803297B (en) | Airworthiness test method for single-shot performance of CCAR29 helicopters | |
| CN204270138U (en) | Steam turbine governing system servo card controling parameters test macro | |
| CN105406788B (en) | A kind of generator dominant parameters discrimination method based on power grid dynamic characteristic | |
| CN110943485B (en) | Index evaluation method for simulation reliability of equivalent model of doubly-fed wind power plant | |
| CN106547985B (en) | Application of SignalCalc in excitation system static model identification | |
| CN109407543B (en) | Verification method and device for voltage response characteristics of electrical model of wind turbine generator |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180725 Address after: 410004 398 new Shao Dong Road, Tianxin District, Changsha, Hunan Co-patentee after: Power Science Research Institute of Hunan Electric Power Co., Ltd. Patentee after: State Grid Hunan Electric Power Co., Ltd. Co-patentee after: State Grid Corporation of China Address before: 100033 West Chang'an Avenue, Xicheng District, Xicheng District, Beijing Co-patentee before: Science Research Institute of Hunan Electric Power Co., Ltd. Patentee before: State Grid Corporation of China Co-patentee before: Hunan Xiangdian Electric Power Text & Research Co., Ltd. |
|
| TR01 | Transfer of patent right |