CN110275110A - A kind of synchronous generator q axis parameter testing and discrimination method - Google Patents
A kind of synchronous generator q axis parameter testing and discrimination method Download PDFInfo
- Publication number
- CN110275110A CN110275110A CN201910534761.8A CN201910534761A CN110275110A CN 110275110 A CN110275110 A CN 110275110A CN 201910534761 A CN201910534761 A CN 201910534761A CN 110275110 A CN110275110 A CN 110275110A
- Authority
- CN
- China
- Prior art keywords
- generator
- axis
- parameter
- data
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/34—Testing dynamo-electric machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
The present invention relates to Parameter Estimation of Synchronous Machines fields, more particularly to a kind of synchronous generator q axis parameter testing and discrimination method, step includes: adjustment generator operating condition, so that generator rotor angle δ is equal with power-factor angle ψ, load dump test is carried out under the specific operation, records stator voltage, excitation voltage, exciting current, active change curve, idle change curve at this time；It is established and the consistent simulation model of actual condition using BPA simulation software；Carry out cutting load using BPA simulation software to test, obtains the generator unit stator voltage data in simulation process；It adjusts separately q axis Steady-state Parameters, transient state parameter, secondary transient state parameter to be emulated, substep Curve fitting simulation curve and measured curve, finally recognizes generator q axis full set parameter.The present invention has fully considered the variation of the exciting current in load dump test dynamic process, and compared to traditional graphic-arts technique discrimination method, the parameter of identification is more accurate.
Description
Technical field
The present invention relates to Parameter Estimation of Synchronous Machines fields, and in particular to a kind of synchronous generator q axis parameter testing and
Discrimination method.
Background technique
Synchronous generator is most important one kind equipment in electric system.Whether reasonable, model parameter is generator model
It is no accurate, the credibility of Power System Stability Analysis result is not only directly affected, but also directly affect such as power system stability
Device (Power system stabilizer, PSS) etc. controls the control strategy design of equipment.
The parameter of synchronous generator can be determined by following three kinds of modes: calculate ginseng using the information in design of electrical motor stage
Number, choose handbook in typical data, under certain research technique, according to the motor response data System Discrimination of measurement
Method estimates parameter.Compared with first two determines the mode of parameter, the third mode has more realistic meaning.This is because: nothing
By being canonical parameter on design of electrical motor data or handbook, all there may be bigger difference with the parameter under practical operation situation,
So as to cause the result and actual conditions grave fault of simulation analysis；In service life, the parameter of generator is not to fix not
Become, due to factors such as physical material aging, overhauls, parameter will change, and with the means of identification, can relatively easily obtain
To the parameter value of continuous renewal；It, can with measurement data with the development of system identification theory, measuring technique and test method
Obtain accurate parameter value.
Industry and academia all take much count of the parameter identification work of synchronous generator.The North America electric reliability committee
(North American Electric Reliability Council, NERC) requires its industry member regularly to test
Generator parameter.International Electrotechnical Commission (International Electrotechnical Commission, IEC) and beauty
State's Institute of Electrical and Electronics Engineers (Institute of Electrical anq Electronics Engineers,
IEEE standard IEC 34-4 and the IEEE Stq 115 for generator parameter test) has been published respectively.The newest state that China promulgates
In family standard GB/T1029-2005, the method that alsies specify parameter testing.The end of the seventies in last century, U.S. professor qemello mentioned
Going out the method for carrying out generator parameter test using removal of load method, removal of load method clear physics conception, field test are easy, and
And consider saturation effect etc. and influence, a kind of method for practical of can yet be regarded as.Generator parameter identification based on removal of load method
Method it is most basic be graphic-arts technique, be common method in Practical Project.But for from shunt dynamo, although in removal of load
During test, field regulator is in permanent excitation con-trol mode, but exciting current will change during its test,
It is very difficult accurately to find out characteristic point from response curve using graphic-arts technique in this case.
Summary of the invention
To solve the above-mentioned problems, the present invention provides a kind of synchronous generator q axis parameter testing and discrimination methods, specifically
Technical solution is as follows:
A kind of synchronous generator q axis parameter testing and discrimination method, comprising the following steps:
S1: adjustment generator operating condition carries out removal of load so that generator rotor angle δ is equal with power-factor angle ψ under the specific operation
Test records stator voltage, excitation voltage, exciting current, active change curve, idle change curve at this time；
S2: it is established and the consistent simulation model of actual condition, including generator basic parameter, excitation system using BPA simulation software
System model, load model, cutting load model；
S3: carry out cutting load using BPA simulation software and test, obtain the generator unit stator voltage data in simulation process；
S4: adjusting separately q axis Steady-state Parameters, transient state parameter, secondary transient state parameter and emulated, substep Curve fitting simulation curve and actual measurement
Curve finally recognizes generator q axis full set parameter.
Preferably, the step S1 the following steps are included:
S11: generator is made to be in grid connection state；
S12: estimating operating condition of test, simplifies actual measurement power-factor angle and generator rotor angle；
S13: adjusting field regulator and be in MANUAL CONTROL mode, operating condition is adjusted by above-mentioned estimation result, so that generator
Excitation keeps power-factor angle and generator rotor angle equal；
S14: disconnecting generator outlet breaker, while recording stator voltage, excitation voltage, exciting current, active variation song
Line, idle change curve.
Preferably, shown in estimating that detailed process is as follows in the step S12:
The experimental condition of q axis parameter testing is generator rotor angle δ equal with power-factor angle ψ, according to the generator rotor angle equation of generator and is got rid of negative
Known to the basic principle of lotus method:
；(1)
；(2)
Work as U0It is when per unit value is 1, then available by two formulas above:
；(3)
The then relationship of reactive power Q and active-power P are as follows:
；(4)
In above formula,E q For electromotive force；U 0 For removal of load presteady state voltage；x q For q axis reactance；δFor generator rotor angle；U q0 After removal of load
Steady state voltage；ψFor power-factor angle.
Preferably, the step S2 the following steps are included:
S21: one machine infinity bus system trend is established using BPA Load Flow Program；
S22: and then simulation model identical with actual condition is established using BPA stability program；
With value nodes such as the locking adopted Infinite bus systems of BS in S23:BPA Load Flow Program, with the locking justice PQ node of B, with L card definition wires
Circuit-switched data fills in corresponding card content according to actual condition；
The locking adopted Infinite bus system of MC is used in S24:BPA stability program, with the locking adopted generator basic parameter of M, MF, excitation system
Model selection E* card is inputted live excitation voltage test data as a curve；
S25: defining load model using LB card, realizes that cutting load is operated with LS card, is then filled out according to live actual condition
Write corresponding card content.
Preferably, the step S4 the following steps are included:
S41: cutting load will be carried out using producer's given parameters and emulates resulting stator voltage data and the progress of field measurement data
Comparison fitting；
Then S42: q axis reactance xq in adjustment simulation model is emulated, until gained emulation data and field measurement data exist
The fitting of stable state section is consistent, then q axis reactance xq data adjusted are accurate parameters；
S43: similarly, adjusting q axis transient state reactance xq' in simulation model, then emulated, until gained emulation data and scene
Measured data is consistent in the fitting of transient state section, then q axis transient state reactance xq' data adjusted are accurate parameters；
S44: similarly, q axis subtranient reactance xq ", super transition open circuit time constant Tq0 " in simulation model is adjusted, is then imitated
Very, it is fitted unanimously up to gained emulates data and field measurement data in secondary transient state section, then q axis subtranient reactance adjusted
Xq ", super transition open circuit time constant Tq0 " data are accurate parameters；
S41: generator q axis full set parameter to sum up then can accurately be picked out.
The invention has the benefit that the present invention is by field test with emulating the method that combines to synchronous generator
Q axis parameter is recognized, and the variation of the exciting current in load dump test dynamic process has been fully considered, compared to traditional
The parameter of graphic-arts technique discrimination method, identification is more accurate, and traditional pattern recognition method can not accurately be found out on trial curve
Stable state, transient state, the characteristic point in secondary transient state each stage can not accurately calculate d axis relevant parameter, and use field test and emulation
The method combined can really reflect live actual operating mode.
Detailed description of the invention
Fig. 1 is flow diagram of the invention；
Fig. 2 is q axis load dump test scene stator voltage, excitation voltage waveform diagram；
Fig. 3 is q axis simulation recognition parameter, traditional graph method identified parameters and measured waveform comparison diagram in BPA.
Specific embodiment
In order to better understand the present invention, the present invention will be further explained below with reference to the attached drawings and specific examples:
As shown in Figure 1, a kind of synchronous generator q axis parameter testing and discrimination method, comprising the following steps:
S1: adjustment generator operating condition carries out removal of load so that generator rotor angle δ is equal with power-factor angle ψ under the specific operation
Test records stator voltage, excitation voltage, exciting current, active change curve, idle change curve at this time.Including with
Lower step:
S11: generator is made to be in grid connection state；
S12: estimating operating condition of test, simplifies actual measurement power-factor angle and generator rotor angle；Estimation detailed process is as follows institute
Show:
The experimental condition of q axis parameter testing is generator rotor angle δ equal with power-factor angle ψ, according to the generator rotor angle equation of generator and is got rid of negative
Known to the basic principle of lotus method:
；(1)
；(2)
Work as U0It is when per unit value is 1, then available by two formulas above:
；(3)
The then relationship of reactive power Q and active-power P are as follows:
；(4)
In above formula,E q For electromotive force；U 0 For removal of load presteady state voltage；x q For q axis reactance；δFor generator rotor angle；U q0 After removal of load
Steady state voltage；ψFor power-factor angle；
S13: adjusting field regulator and be in MANUAL CONTROL mode, operating condition is adjusted by above-mentioned estimation result, so that generator
Excitation keeps power-factor angle and generator rotor angle equal；
S14: disconnecting generator outlet breaker, while recording stator voltage, excitation voltage, exciting current, active variation song
Line, idle change curve.
S2: using BPA simulation software establish with the consistent simulation model of actual condition, including generator basic parameter, encourage
Magnetic system model, load model, cutting load model.The following steps are included:
S21: one machine infinity bus system trend is established using BPA Load Flow Program；
S22: and then simulation model identical with actual condition is established using BPA stability program；
With value nodes such as the locking adopted Infinite bus systems of BS in S23:BPA Load Flow Program, with the locking justice PQ node of B, with L card definition wires
Circuit-switched data fills in corresponding card content according to actual condition；
The locking adopted Infinite bus system of MC is used in S24:BPA stability program, with the locking adopted generator basic parameter of M, MF, excitation system
Model selection E* card is inputted live excitation voltage test data as a curve；Wherein, using the locking adopted generator of M, MF
Producer's given parameters are filled in when basic parameter first；
S25: defining load model using LB card, realizes that cutting load is operated with LS card, is then filled out according to live actual condition
Write corresponding card content.
S3: carry out cutting load using BPA simulation software and test, obtain the generator unit stator voltage data in simulation process.
S4: adjusting separately q axis Steady-state Parameters, transient state parameter, secondary transient state parameter and emulated, substep Curve fitting simulation curve with
Measured curve finally recognizes generator q axis full set parameter.The following steps are included:
S41: cutting load will be carried out using producer's given parameters and emulates resulting stator voltage data and the progress of field measurement data
Comparison fitting；
Then S42: q axis reactance xq in adjustment simulation model is emulated, until gained emulation data and field measurement data exist
The fitting of stable state section is consistent, then q axis reactance xq data adjusted are accurate parameters；
S43: similarly, adjusting q axis transient state reactance xq' in simulation model, then emulated, until gained emulation data and scene
Measured data is consistent in the fitting of transient state section, then q axis transient state reactance xq' data adjusted are accurate parameters；
S44: similarly, q axis subtranient reactance xq ", super transition open circuit time constant Tq0 " in simulation model is adjusted, is then imitated
Very, it is fitted unanimously up to gained emulates data and field measurement data in secondary transient state section, then q axis subtranient reactance adjusted
Xq ", super transition open circuit time constant Tq0 " data are accurate parameters；
S41: generator q axis full set parameter to sum up then can accurately be picked out.Emulate data and measured waveform fitting comparing result such as
Shown in Fig. 3.Q axis parameter identification is carried out with certain hydraulic turbine, the parameter picked out is as shown in table 1, can be seen that emulation is distinguished in conjunction with Fig. 3
Parameter that parameter is picked out compared to traditional graph method is known closer to field measurement waveform, illustrates that the parameter picked out is more smart
Really, more meet actual operating mode.
1 q axis identified parameters of table
| D/q axis parameter | The parameter that graphic-arts technique picks out | Simulation recognition parameter |
| x q | 0.731 | 0.76 |
| x q ' | 0.731 | 0.76 |
| x q " | 0.281 | 0.290 |
| T q0 ' | 0.07 | 0.09 |
| T q0 " | 0.132 | 0.143 |
The present invention is not limited to above-described specific embodiment, the foregoing is merely preferable case study on implementation of the invention
, it is not intended to limit the invention, any modifications, equivalent replacements, and improvements done within the spirit and principles of the present invention
Deng should all be included in the protection scope of the present invention.
Claims (5)
1. a kind of synchronous generator q axis parameter testing and discrimination method, it is characterised in that: the following steps are included:
S1: adjustment generator operating condition carries out removal of load so that generator rotor angle δ is equal with power-factor angle ψ under the specific operation
Test records stator voltage, excitation voltage, exciting current, active change curve, idle change curve at this time；
S2: it is established and the consistent simulation model of actual condition, including generator basic parameter, excitation system using BPA simulation software
System model, load model, cutting load model；
S3: carry out cutting load using BPA simulation software and test, obtain the generator unit stator voltage data in simulation process；
S4: adjusting separately q axis Steady-state Parameters, transient state parameter, secondary transient state parameter and emulated, substep Curve fitting simulation curve and actual measurement
Curve finally recognizes generator q axis full set parameter.
2. a kind of synchronous generator q axis parameter testing according to claim 1 and discrimination method, it is characterised in that: described
Step S1 the following steps are included:
S11: generator is made to be in grid connection state；
S12: estimating operating condition of test, simplifies actual measurement power-factor angle and generator rotor angle；
S13: adjusting field regulator and be in MANUAL CONTROL mode, operating condition is adjusted by above-mentioned estimation result, so that generator
Excitation keeps power-factor angle and generator rotor angle equal；
S14: disconnecting generator outlet breaker, while recording stator voltage, excitation voltage, exciting current, active variation song
Line, idle change curve.
3. a kind of synchronous generator q axis parameter testing according to claim 2 and discrimination method, it is characterised in that: described
Shown in estimating that detailed process is as follows in step S12:
The experimental condition of q axis parameter testing is generator rotor angle δ equal with power-factor angle ψ, according to the generator rotor angle equation of generator and is got rid of negative
Known to the basic principle of lotus method:
；(1)
；(2)
Work as U0It is when per unit value is 1, then available by two formulas above:
；(3)
The then relationship of reactive power Q and active-power P are as follows:
；(4)
In above formula,E q For electromotive force；U 0 For removal of load presteady state voltage；x q For q axis reactance；δFor generator rotor angle；U q0 After removal of load
Steady state voltage；ψFor power-factor angle.
4. a kind of synchronous generator q axis parameter testing according to claim 1 and discrimination method, it is characterised in that: described
Step S2 the following steps are included:
S21: one machine infinity bus system trend is established using BPA Load Flow Program；
S22: and then simulation model identical with actual condition is established using BPA stability program；
With value nodes such as the locking adopted Infinite bus systems of BS in S23:BPA Load Flow Program, with the locking justice PQ node of B, with L card definition wires
Circuit-switched data fills in corresponding card content according to actual condition；
The locking adopted Infinite bus system of MC is used in S24:BPA stability program, with the locking adopted generator basic parameter of M, MF, excitation system
Model selection E* card is inputted live excitation voltage test data as a curve；
S25: defining load model using LB card, realizes that cutting load is operated with LS card, is then filled out according to live actual condition
Write corresponding card content.
5. a kind of synchronous generator q axis parameter testing according to claim 1 and discrimination method, it is characterised in that: described
Step S4 the following steps are included:
S41: cutting load will be carried out using producer's given parameters and emulates resulting stator voltage data and the progress of field measurement data
Comparison fitting；
Then S42: q axis reactance xq in adjustment simulation model is emulated, until gained emulation data and field measurement data exist
The fitting of stable state section is consistent, then q axis reactance xq data adjusted are accurate parameters；
S43: similarly, adjusting q axis transient state reactance xq' in simulation model, then emulated, until gained emulation data and scene
Measured data is consistent in the fitting of transient state section, then q axis transient state reactance xq' data adjusted are accurate parameters；
S44: similarly, q axis subtranient reactance xq ", super transition open circuit time constant Tq0 " in simulation model is adjusted, is then imitated
Very, it is fitted unanimously up to gained emulates data and field measurement data in secondary transient state section, then q axis subtranient reactance adjusted
Xq ", super transition open circuit time constant Tq0 " data are accurate parameters；
S41: generator q axis full set parameter to sum up then can accurately be picked out.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910534761.8A CN110275110A (en) | 2019-06-20 | 2019-06-20 | A kind of synchronous generator q axis parameter testing and discrimination method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910534761.8A CN110275110A (en) | 2019-06-20 | 2019-06-20 | A kind of synchronous generator q axis parameter testing and discrimination method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110275110A true CN110275110A (en) | 2019-09-24 |
Family
ID=67961110
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910534761.8A Pending CN110275110A (en) | 2019-06-20 | 2019-06-20 | A kind of synchronous generator q axis parameter testing and discrimination method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110275110A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101430365A (en) * | 2008-12-12 | 2009-05-13 | 南京工程学院 | Identification system and method for actually measured electric parameter of synchronous generator |
| EP2270522A1 (en) * | 2009-07-03 | 2011-01-05 | ABB Oy | Induction motor parameter identification |
| CN102073012A (en) * | 2009-11-20 | 2011-05-25 | 华北电力科学研究院有限责任公司 | Method and system for obtaining parameters of synchronous generator |
| CN102073014A (en) * | 2011-01-07 | 2011-05-25 | 华北电网有限公司 | Method for identifying synchronous generator parameters |
| CN102510263A (en) * | 2011-10-25 | 2012-06-20 | 中国电力科学研究院 | Method for identifying practical parameters of synchronous generator on basis of load rejection test and numerical difference |
| CN102904518A (en) * | 2012-09-27 | 2013-01-30 | 安徽省电力公司 | Synchronous generator q shaft parameter on-line identification method |
-
2019
- 2019-06-20 CN CN201910534761.8A patent/CN110275110A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101430365A (en) * | 2008-12-12 | 2009-05-13 | 南京工程学院 | Identification system and method for actually measured electric parameter of synchronous generator |
| EP2270522A1 (en) * | 2009-07-03 | 2011-01-05 | ABB Oy | Induction motor parameter identification |
| CN102073012A (en) * | 2009-11-20 | 2011-05-25 | 华北电力科学研究院有限责任公司 | Method and system for obtaining parameters of synchronous generator |
| CN102073014A (en) * | 2011-01-07 | 2011-05-25 | 华北电网有限公司 | Method for identifying synchronous generator parameters |
| CN102510263A (en) * | 2011-10-25 | 2012-06-20 | 中国电力科学研究院 | Method for identifying practical parameters of synchronous generator on basis of load rejection test and numerical difference |
| CN102904518A (en) * | 2012-09-27 | 2013-01-30 | 安徽省电力公司 | Synchronous generator q shaft parameter on-line identification method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110346718A (en) | A kind of synchronous generator d axis parameter testing and discrimination method | |
| Amin et al. | A gray-box method for stability and controller parameter estimation in HVDC-connected wind farms based on nonparametric impedance | |
| CN106533290B (en) | A kind of the stator current method for limiting and device of automatic excitation adjustor of generator | |
| US11714934B2 (en) | Aggregated model of large-scale wind farms for power system simulation software tools | |
| CN104699996A (en) | Electric system asymmetric short circuit power frequency electric quantity calculation method based on double-fed induction wind generation set equivalence | |
| CN106503346B (en) | A kind of generating set electromagnetical transient emulation method of quick and smooth starting | |
| Jain et al. | Grid-supportive loads—a new approach to increasing renewable energy in power systems | |
| Alinejad et al. | PMU-based distribution network load modelling using Harmony Search Algorithm | |
| Tan et al. | Automatic generation and parameter verification of large-scale EMT simulation models based on TSP projects | |
| Ashourian et al. | Controlling and modeling power-electronic interface DERs in islanding mode operation micro grid | |
| CN110275110A (en) | A kind of synchronous generator q axis parameter testing and discrimination method | |
| CN105701265A (en) | Double-fed wind generator modeling method and apparatus | |
| Rengifo et al. | Full-load range in-situ efficiency estimation method for induction motors using only a direct start-up | |
| CN106897514A (en) | A kind of method for building up of the calculation of short-circuit current model of total power changing type new energy station | |
| CN115656658A (en) | Network-based converter inertia response testing device and testing and evaluating method thereof | |
| CN106055817B (en) | Double-fed asynchronous generator Double-working-condition model modelling approach for electromagnetic transient analysis | |
| Han et al. | Nonlinear transient mathematical model of large-capacity synchronous condenser based on time-varying reactance parameters | |
| Diao et al. | Simplified Modeling of an Inverter Air Conditioner based on Experimental Data | |
| Tao et al. | Transient stability analysis of AC/DC system considering electromagnetic transient model | |
| Wang et al. | Modeling of load frequency characteristics in the load model for power system digital simulation | |
| Wang et al. | Power emulator of variable speed drive with grid frequency support in multi-converter based power grid emulation system | |
| Subiabre et al. | Realistic loss modelling and minimisation in an air-cored permanent magnet generator for wind energy applications | |
| Gao et al. | Improved extended kalman filter based dynamic equivalent method of DFIG wind farm cluster | |
| Wamkeue et al. | Alternative approaches for linear analysis and prediction of a synchronous generator under partial-and full-load rejection tests | |
| Khan et al. | Performance comparison of single winding and double winding self-excited induction generators |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190924 |