CN106452238B - A kind of Siemens PSS3B model parameter setting methods - Google Patents
A kind of Siemens PSS3B model parameter setting methods Download PDFInfo
- Publication number
- CN106452238B CN106452238B CN201610884501.XA CN201610884501A CN106452238B CN 106452238 B CN106452238 B CN 106452238B CN 201610884501 A CN201610884501 A CN 201610884501A CN 106452238 B CN106452238 B CN 106452238B
- Authority
- CN
- China
- Prior art keywords
- phase
- pss3b
- siemens
- pss
- frequency
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000005540 biological transmission Effects 0.000 claims abstract description 22
- 230000005284 excitation Effects 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000010363 phase shift Effects 0.000 claims description 9
- 230000000903 blocking effect Effects 0.000 claims description 5
- 230000003321 amplification Effects 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 238000013016 damping Methods 0.000 abstract description 7
- 230000010355 oscillation Effects 0.000 abstract description 7
- 239000003381 stabilizer Substances 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P9/00—Arrangements for controlling electric generators for the purpose of obtaining a desired output
- H02P9/14—Arrangements for controlling electric generators for the purpose of obtaining a desired output by variation of field
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of Siemens PSS3B model parameter setting methods, include the following steps：1, using Mason's formula, transmission function expression formula is obtained according to Siemens's PSS3B model framework charts；2, Siemens's PSS3B transmission functions are converted to lead-lag link series connection transmission function form；3, the uncompensated phase-frequency characteristic φ of actual measurement generator excited systemΔEq’, the uncompensated phase-frequency characteristic of excitation system being applied in engineering practice is obtained after carrying out " artificial smooth " to the phase significantly change near intrinsic point；4, it is based on Non-linear least-square curve fitting φPSS(Δ Pe) phase frequency function.5, processing fitting result obtains PSS3B setting parameters.The PSS3B parameters adjusted using this method can effectively be played power system stabilizer, PSS and improve power grid dynamic stability, play the role of power oscillation damping.
Description
Technical field
The invention belongs to net source coordination technical field, it is related to a kind of using the method for undetermined coefficients and nonlinear least square method
Siemens's PSS3B parameter tuning methods.
Background technology
Due to the delayed phase characteristic of field regulator, excitation system and generator field winding, field regulator generates
Phase lags behind generator rotor angle and the negative damping torque with rotation speed change antiphase, so as to cause low-frequency oscillation.Power system stability
Device (PSS) is a kind of additional excitation control technology studied for power oscillation damping.Power system stabilizer, PSS introducing is led over
The additional signal of axle speed overcomes the negative damping torque generated in former excitation voltage regulation device, is synthetically produced positive damping torque, solution
Certainly low-frequency oscillation problem is one of the important measures for improving power system dynamic stability.
Siemens's PSS3B models use unconventional form, not common lead-lag link cascade, Duo Gecan
Numerical digit causes PSS3B forms complicated in multiple concatenated feedback elements.Have to PSS phase compensations in power industry standard quantitative
It is required that, it is therefore desirable to Siemens's PSS3B model conversations are subjected to phase compensation for common lead-lag link.
When configuring PSS using phase compensation principle, need to determine the phase of PSS according to the uncompensated phase-frequency characteristic of excitation system
Frequency characteristic, and then adjust PSS parameter.If PSS phase-frequency characteristics areExcitation system is uncompensated, and phase-frequency characteristic isThen
It is theoretically rightRequirement be：
1) using Δ ω as the PSS of input signal, have
2) using-Δ Pe as the PSS of input signal, have
When being adjusted to PSS3B progress scene, the uncompensated phase-frequency characteristic of generator excited system is surveyedCalculate PSS
Phase-frequency characteristic isIt is mended using the method design PSS phases of fitting PSS phase frequency curves
Repay parameter.
Invention content
The technical problem to be solved by the present invention is to, in view of the shortcomings of the prior art, provide it is a kind of using the method for undetermined coefficients and
Siemens's PSS3B parameter tuning methods of nonlinear least square method, to meet the needs of field test.
To achieve the above object, the method applied in the present invention is：
A kind of Siemens PSS3B model parameter setting methods, include the following steps：
Step 1 obtains transmission function expression formula using Mason's formula according to Siemens's PSS3B model framework charts；
Step 2 converts the transmission function in step 1 to lead-lag link series connection transmission function form；
Step 3 surveys the uncompensated phase-frequency characteristic of generator excited system, and to the phase significantly change near intrinsic point
Carry out it is artificial smooth after obtain the uncompensated phase-frequency characteristic of excitation system being applied in engineering practice, then obtain PSS phases accordingly
Frequency characteristic；
Step 4, based on the phase frequency function for expressing PSS phase-frequency characteristics in Non-linear least-square curve fitting step 3, then
Series connection lead-lag phase compensation link, finally obtains the phase frequency function after compensation；
Step 5 obtains PSS3B according to the Relationship of Coefficients of fitting result in step 4 and conversion process in step 2 and adjusts
Parameter.
The method, in the step one, transmission function expression formula is：
Wherein,
s：Transmission function independent variable；
K：Gain coefficient；
TD:Blocking link time constant；
K0~K4:Phase shift link proportionality coefficient；
T1~T4:Phase shift link time constant.
The method, in the step two, the function expression after conversion is：
K in conversion0=1 and there are Relationship of Coefficients：
Wherein T1'~T4':The phase shift link time constant of transmission function after conversion；
A1~A4:Intermediate link coefficient；
B1~B4:Intermediate link coefficient.
The method, in the step three, generator excited system is uncompensated, and phase-frequency characteristic isPSS phases
Frequency characteristicFor：
The method, in the step four, the phase frequency function expression after compensation is：
F (f)=90 °-arctan (2 π fTD)+2[arctan(2πfT1')-arctan(2πfT2')+arctan(2π
fT3')-arctan(2πfT4')]。
The method, in the step five, PSS3B setting parameters are TD, T1、T2、T3、T4、K1、K2、K3、K4、K0,
Its large deviations amplification coefficient K is determined after being tested by critical gain.
The technical effects of the invention are that, it is proposed that a method of it being directed to Siemens's PSS3B model parameters and adjusts, energy
The phase compensation for providing 0.1~2Hz frequency ranges effectively improves the ability for inhibiting unit low-frequency oscillation.
Description of the drawings
Fig. 1 is implementing procedure figure of the present invention；
Fig. 2 is Siemens's PSS3B illustratons of model；
Fig. 3, which is PSS, should compensate phase and its matched curve figure；
Fig. 4 is that whether there is or not PSS microvariations active power response curves；
Fig. 5 is artificial smooth schematic diagram.
Specific implementation mode
As shown in Figure 1,1) utilize Mason's formula, transmission function expression formula is obtained according to Siemens's PSS3B model framework charts.
s：Transmission function independent variable；
K：Gain coefficient；
TD:Blocking link time constant；
K0~K4:Phase shift link proportionality coefficient；
T1~T4:Phase shift link time constant；
2) Siemens's PSS3B transmission functions are converted to lead-lag link series connection transmission function form.
The present embodiment is subtracted by the form of artificially specified lead-lag series connection link in the case where compensation effect is constant
The number for having lacked time constant undetermined, simplifies transfer process.
3) the uncompensated phase-frequency characteristic of generator excited system is surveyedTo the phase significantly change near intrinsic point into
The uncompensated phase-frequency characteristic of excitation system being applied in engineering practice is obtained after row " artificial smooth "；Calculate PSS phase-frequency characteristicsWherein artificial refers to smoothly that the SPA sudden phase anomalies point in phase-frequency characteristic is not examined
Consider, it is carried out thinking to be fitted according to trend before, for details, reference can be made to Fig. 5, the phase of the neighbouring significantly changes of 1.2Hz in figure,
The curve of oscillation point both sides is joined together using thickened portion curve is manually smoothly obtained.In addition, although transmission function is s
Function, but need to replace s to be substituted into jw when considering the transmission function of phasor relation, w=2 π f herein, wherein frequency f's
Ranging from 0.1~2.0Hz.
4) it is based on Non-linear least-square curve fittingTransformed series connection lead-lag phase compensation ring
The phase frequency function of section is
F (f)=90 °-arctan (2 π fTD)+2[arctan(2πfT1')-arctan(2πfT2')+arctan(2π
fT3')-arctan(2πfT4')] utilize the lsqcurvefit function pairs in matlabCarry out non-linear minimum
Square law is fitted.It is selected there are many ways to being wherein fitted, for the considerations of simplifying, the present embodiment directly uses
Least square function inside matlab is fitted.
Blocking link phase is contained in the present embodiment fitting function, that is, considers the time constant of blocking link.
5) it is PSS3B model parameters by lead-lag link Parameter Switch.Transformational relation is
T1'~T4':The phase shift link time constant of transmission function after conversion；
A1~A4:Intermediate link coefficient；
B1~B4:Intermediate link coefficient；
The present embodiment, which exempts from solution nonlinear equation, can carry out Parameter Switch.
Deviation amplification coefficient K is determined after being tested by critical gain.
Calculated examples：
The theoretical value for the uncompensated phase-frequency characteristic of excitation system that one one-machine infinite-bus power system of structure calculates its generator,
Uncompensated phase-frequency characteristic curve after " artificial smooth " within the scope of 0.1~2Hz of frequency take a little as shown in table 1,
1 uncompensated phase-frequency characteristic of table
| f/Hz | Phase (°) | f/Hz | Phase (°) |
| 0.1 | -17 | 1.1 | -77 |
| 0.2 | -32 | 1.2 | -79 |
| 0.3 | -43 | 1.3 | -80 |
| 0.4 | -52 | 1.4 | -80 |
| 0.5 | -58 | 1.5 | -80 |
| 0.6 | -64 | 1.6 | -80 |
| 0.7 | -68 | 1.7 | -80 |
| 0.8 | -71 | 1.8 | -80 |
| 0.9 | -73 | 1.9 | -80 |
| 1 | -75 | 2 | -80 |
According to formulaIt is as shown in table 2 that calculating PSS should compensate phase
2 PSS of table should compensate phase
| f/Hz | Phase (°) | f/Hz | Phase (°) |
| 0.1 | -73 | 1.1 | -13 |
| 0.2 | -58 | 1.2 | -11 |
| 0.3 | -47 | 1.3 | -10 |
| 0.4 | -38 | 1.4 | -10 |
| 0.5 | -32 | 1.5 | -10 |
| 0.6 | -26 | 1.6 | -10 |
| 0.7 | -22 | 1.7 | -10 |
| 0.8 | -19 | 1.8 | -10 |
| 0.9 | -17 | 1.9 | -10 |
| 1 | -15 | 2 | -10 |
Based on Non-linear least-square curve fittingTransformed series connection lead-lag phase compensation link
Phase frequency function be
F (f)=90 °-arctan (2 π fTD)+2[arctan(2πfT1')-arctan(2πfT2')+arctan(2π
fT3')-arctan(2πfT4')] utilize the lsqcurvefit function pairs in matlabCarry out non-linear minimum
Square law is fitted.Matched curve is as shown in Figure 3.
Fitting result is handled, it is as shown in table 3 to obtain PSS3B model parameters.Deviation amplification coefficient K passes through critical increasing
1.4 are finally taken after beneficial test simulation.
3 setting parameter of table
| Parameter | Numerical value | Parameter | Numerical value |
| TD | 5 | T4 | 10.835 |
| T1’ | 0.8437 | K0 | 1 |
| T2’ | 5.2324 | K1 | 1.1692 |
| T3’ | 0.1868 | K2 | 5.7092 |
| T4’ | 0.1851 | K3 | 7.1871 |
| T1 | 0.0894 | K4 | 37.782 |
| T2 | 0.3342 | K | 1.4 |
| T3 | 2.8991 |
It according to data setting PSS3B models in table 3 and carries out whether there is or not the microvariations of PSS emulation, microvariations form is excitation
4% step of reference voltage, simulated effect are as shown in Figure 4.Table 4 is that the small interference frequency domain of the corresponding one-machine infinite-bus power systems of Fig. 4 is special
Value indicative result of calculation.The result shows that positive damping increases after input PSS, play the role of inhibition to low-frequency oscillation.
Whether there is or not PSS unit active power analysis results for table 4
| State | Damping ratio | Frequency | Real part | Imaginary part |
| Without PSS steps | 0.071 | 1.3261 | -0.5933 | 8.3318 |
| There are PSS steps | 0.252 | 1.924 | -3.148 | 12.088 |
Claims (4)
1. a kind of Siemens PSS3B model parameter setting methods, which is characterized in that include the following steps：
Step 1 obtains transmission function expression formula using Mason's formula according to Siemens's PSS3B model framework charts；
Step 2 converts the transmission function in step 1 to lead-lag link series connection transmission function form；
Step 3 is surveyed the uncompensated phase-frequency characteristic of generator excited system, and is carried out to the phase significantly change near intrinsic point
The uncompensated phase-frequency characteristic of excitation system being applied in engineering practice is obtained after artificial smooth, and it is special then to obtain PSS phase frequencies accordingly
Property；
Then step 4 is connected based on the phase frequency function for expressing PSS phase-frequency characteristics in Non-linear least-square curve fitting step 3
Lead-lag phase compensation link finally obtains the phase frequency function after compensation；
Step 5 obtains PSS3B setting parameters according to the Relationship of Coefficients of conversion process in fitting result in step 4 and step 2；
In the step one, transmission function expression formula is：
Wherein,
s：Transmission function independent variable；
K：Gain coefficient；
TD:Blocking link time constant；
K0~K4:Phase shift link proportionality coefficient；
T1~T4:Phase shift link time constant；
In the step two, the function expression after conversion is：
K in conversion0=1 and there are Relationship of Coefficients：
Wherein T1'~T4':The phase shift link time constant of transmission function after conversion；
A1~A4:Intermediate link coefficient；
B1~B4:Intermediate link coefficient.
2. according to the method described in claim 1, it is characterized in that, in the step three, generator excited system is uncompensated
Phase-frequency characteristic isPSS phase-frequency characteristicsFor：
3. according to the method described in claim 1, it is characterized in that, in the step four, the phase frequency function representation after compensation
Formula is：
F (f)=90 °-arctan (2 π fTD)+2[arctan(2πfT1')-arctan(2πfT2')+arctan(2πfT3')-
arctan(2πfT4')]。
4. according to the method described in claim 1, it is characterized in that, in the step five, PSS3B setting parameters are TD, T1、
T2、T3、T4、K1、K2、K3、K4、K0, determination after large deviations amplification coefficient K is tested by critical gain.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610884501.XA CN106452238B (en) | 2016-10-10 | 2016-10-10 | A kind of Siemens PSS3B model parameter setting methods |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610884501.XA CN106452238B (en) | 2016-10-10 | 2016-10-10 | A kind of Siemens PSS3B model parameter setting methods |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106452238A CN106452238A (en) | 2017-02-22 |
| CN106452238B true CN106452238B (en) | 2018-08-10 |
Family
ID=58172339
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610884501.XA Active CN106452238B (en) | 2016-10-10 | 2016-10-10 | A kind of Siemens PSS3B model parameter setting methods |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106452238B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110610024B (en) * | 2019-08-09 | 2023-04-07 | 国网河北省电力有限公司电力科学研究院 | PSS parameter damping effect optimization method and device |
| CN111585276B (en) * | 2020-05-18 | 2021-08-17 | 国网重庆市电力公司电力科学研究院 | PSS parameter online setting method and device and readable storage medium |
| CN111682532A (en) * | 2020-06-12 | 2020-09-18 | 国网重庆市电力公司电力科学研究院 | Excitation system uncompensated phase-frequency characteristic online modeling method and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101447670A (en) * | 2008-05-27 | 2009-06-03 | 中国电力科学研究院 | Setting algorithm of parameters of power system stabilizer |
| CN101447679A (en) * | 2008-09-17 | 2009-06-03 | 中国电力科学研究院 | Method for implementing parallel power system stabilizer |
| CN102801175A (en) * | 2012-08-09 | 2012-11-28 | 国电南瑞科技股份有限公司 | PSS (Power System Stabilizer) phase compensation link time constant calculation method |
| CN103532130A (en) * | 2013-09-06 | 2014-01-22 | 国家电网公司 | III-type PSS (Power System Stabilizer) parameter setting method |
| WO2014154305A1 (en) * | 2013-03-28 | 2014-10-02 | Siemens Aktiengesellschaft | Method and apparatus for a load step robust unit commitment of power generation units in a power supply system |
-
2016
- 2016-10-10 CN CN201610884501.XA patent/CN106452238B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101447670A (en) * | 2008-05-27 | 2009-06-03 | 中国电力科学研究院 | Setting algorithm of parameters of power system stabilizer |
| CN101447679A (en) * | 2008-09-17 | 2009-06-03 | 中国电力科学研究院 | Method for implementing parallel power system stabilizer |
| CN102801175A (en) * | 2012-08-09 | 2012-11-28 | 国电南瑞科技股份有限公司 | PSS (Power System Stabilizer) phase compensation link time constant calculation method |
| WO2014154305A1 (en) * | 2013-03-28 | 2014-10-02 | Siemens Aktiengesellschaft | Method and apparatus for a load step robust unit commitment of power generation units in a power supply system |
| CN103532130A (en) * | 2013-09-06 | 2014-01-22 | 国家电网公司 | III-type PSS (Power System Stabilizer) parameter setting method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106452238A (en) | 2017-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fan | Modeling type-4 wind in weak grids | |
| CN107679769B (en) | Method for establishing frequency response model of power system containing wind power and method for calculating frequency characteristic index | |
| CN106452238B (en) | A kind of Siemens PSS3B model parameter setting methods | |
| CN106849135B (en) | The power/current quality Synergistic method of microgrid inverter and active filter | |
| CN106532685A (en) | Generalized impedance criterion calculation method for stability analysis of grid-connected inverter and application | |
| CN105470979B (en) | A kind of method and system inhibiting double-fed fan motor unit sub-synchronous oscillation | |
| CN108880300B (en) | Double-fed fan rectifier impedance calculation method based on double closed-loop control | |
| CN106786674B (en) | Double-fed blower-compensated transmission system subsynchronous resonance suppressing method and device | |
| CN106786673B (en) | The suppressing method and device of double-fed blower compensated transmission system subsynchronous resonance | |
| CN106712055A (en) | Power system stabilizer (PSS) configuration method in coordination with low-excitation limiting function | |
| CN115313524A (en) | Photovoltaic power generation grid-connected control method and system based on network-building type converter | |
| CN109217371A (en) | Consider voltage converter grid-connected system method for analyzing stability, apparatus and system that phaselocked loop influences | |
| Zhou et al. | Mitigation of subsynchronous oscillation in a VSC-HVDC connected offshore wind farm integrated to grid | |
| CN102522751A (en) | Control method for actively suppressing power grid harmonics | |
| CN107579529A (en) | A kind of subsynchronous suppressing method of synchronous machine based on the optimization of grid-connection converter phaselocked loop | |
| CN110970925A (en) | Double-fed fan based damping and modeling method for improving system through fast active power modulation | |
| CN112202186A (en) | Machine network coordination feedback control method for restraining subsynchronous oscillation of direct-drive fan | |
| CN104201953B (en) | Steam turbine output mechanical power control method in a kind of generating set | |
| CN109062274A (en) | A kind of magnetic bearing oscillatory torque suppressing method based on complex variable finite-dimensional repetitive control | |
| CN106532785A (en) | Load frequency control method of considering novel cloud model | |
| CN106817054A (en) | A kind of PMSG control methods of the mechanical elastic energy storage based on parameter identification | |
| Özdemir et al. | Stability analysis of the automatic voltage regulation system with PI controller | |
| CN106849130A (en) | A kind of setting method of parameters of power system stabilizer | |
| Chatterjee et al. | On the equality of modal damping power and the average rate of transient energy dissipation in a multimachine power system | |
| CN116073437A (en) | Three-phase grid-connected inverter droop control strategy based on Lyapunov equation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |