CN111523229A - Subsynchronous oscillation joint simulation method based on multi-segment centralized quality model - Google Patents
Subsynchronous oscillation joint simulation method based on multi-segment centralized quality model Download PDFInfo
- Publication number
- CN111523229A CN111523229A CN202010323069.3A CN202010323069A CN111523229A CN 111523229 A CN111523229 A CN 111523229A CN 202010323069 A CN202010323069 A CN 202010323069A CN 111523229 A CN111523229 A CN 111523229A
- Authority
- CN
- China
- Prior art keywords
- model
- torsional
- mass
- section
- variation
- 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
Images
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
The invention discloses a subsynchronous oscillation joint simulation method based on a multi-section centralized quality model, and belongs to the technical field of subsynchronous oscillation of a power system. The method has the advantages that the simple shafting model is expanded into the multi-section concentrated mass model, the electric models of the generator and the power system and the data interfaces of the multi-section concentrated mass model and the electric models are established, electromechanical joint simulation under the same simulation platform is achieved, simulation efficiency is high, results are accurate, shafting local stress strain information can be flexibly obtained, the mechanism of the steam turbine generator unit for sub-synchronous oscillation is revealed, and research on the sub-synchronous oscillation mechanism is facilitated.
Description
Technical Field
The invention belongs to the technical field of subsynchronous oscillation of a power system, and particularly relates to a subsynchronous oscillation joint simulation method based on a multi-section centralized quality model.
Background
China's coal resources are mainly concentrated in western regions, and electric power consumption is mainly concentrated in the middle east, so that thermal power and electric energy need to be remotely transmitted to a load center, and in order to improve the electric power transmission capacity, a series compensation alternating current transmission mode or a high-voltage direct current transmission mode is often adopted. The series capacitor is additionally arranged in the alternating-current transmission system, so that partial inductive impedance of a transmission line can be compensated, the electrical distance between generator sets is shortened, the synchronous torque is increased, the voltage drop and the power angle difference caused by power transmission are reduced, the transmission capacity of a line is increased, and the risk of sub-synchronous resonance of a steam turbine generator unit is increased.
The time domain simulation method is a main method for analyzing the subsynchronous oscillation characteristics of the turbo generator set, a system differential equation is solved by adopting machine network coupling analysis software based on PSCAD/EMTDC at present, the system differential equation comprises an electrical system and a mechanical system, wherein a shafting model in the PSCAD/EMTDC adopts a simple quality model, when a centralized parameter model is adopted, if the inherent frequency error caused by aggregation is very small, the total node number and the required analysis order number are required to meet certain conditions, the subsynchronous frequency of the turbo generator set usually comprises the torsional inherent frequency of the first three orders, and if six sections or five sections are adopted, the calculation precision can obviously not be very high; in addition, the simple lumped parameter model can only obtain the torsional vibration responses of a plurality of specific sections, and if the local torque responses need to be analyzed in a detailed mode, simulation results need to be output to other analysis systems. The distribution parameter model is closer to the actual situation and can calculate the influence of the local structure of the shafting, but the calculation is difficult, the calculation amount is large, and the coupling oscillation characteristics between the computer nets are not generally used for analyzing.
Disclosure of Invention
In order to solve the above problems, the present invention aims to provide a sub-synchronous oscillation joint simulation method based on a multi-segment lumped mass model, which has high simulation efficiency and accurate result, can flexibly obtain shafting local stress strain information, reveals a mechanism of the sub-synchronous oscillation of a steam turbine generator unit, and is beneficial to developing a research on the sub-synchronous oscillation mechanism.
The invention is realized by the following technical scheme:
a subsynchronous oscillation joint simulation method based on a multi-segment centralized quality model comprises the following steps:
step 1: acquiring various parameters of a shaft system of a steam turbine generator unit, establishing a multi-section centralized quality model of the shaft system and forming a parameter initialization file of the shaft system model;
step 2: constructing a torsional stress response calculation module of the torsional vibration fault lower shaft system according to the multi-section concentrated mass model established in the step 1;
and step 3: acquiring various parameters of a steam turbine generator unit and an electric power system, and establishing electric models of a generator and the electric power system;
and 4, step 4: establishing a data exchange interface between the multi-section centralized quality model in the step 1 and the electric model in the step 3;
and 5: initializing an electromagnetic torque value in an electric model, and initializing torsion angles, angular velocities and angular acceleration values of each mass block in the multi-section concentrated mass model under the torsional vibration fault;
step 6: starting an electrical model, transmitting the electromagnetic torque value of the current step length obtained by calculation according to the initial value of the step 5 to a multi-section concentrated mass model, subtracting the electromagnetic torque value of the current step length from the electromagnetic torque value of the previous step length to obtain the electromagnetic torque variation and transmitting the electromagnetic torque variation to the torsional stress response calculation module of the step 2;
and 7: the torsional stress response calculation module calculates and obtains the variation of the torsional angle, the angular velocity and the angular acceleration of each mass unit in the multi-section concentrated mass model, and the variation is respectively added with the numerical value of the last step length to obtain the torsional angle, the angular velocity and the angular acceleration value of the current step length;
and 8: calculating the rotating speed of a mass unit where the generator rotor is located, transmitting the rotating speed to the electrical model, and starting the next long simulation; simultaneously outputting torsional stress and torsional strain results of each dangerous section;
and step 9: and (5) repeating the step 6 to the step 8 until the preset total simulation time is reached.
Preferably, the parameters of the steam turbine generator unit shafting comprise structural parameters and material parameters of the shafting.
Preferably, in step 2, the torsional stress response calculation module is constructed based on a transfer matrix method and in an incremental form of an equation, and the input of the torsional stress response calculation module is the electromagnetic torque variation and the output is the variation of the torsional angle, the angular velocity and the angular acceleration of each mass unit.
Preferably, in step 3, the electrical models of the generator and the power system are established by using PSCAD.
Preferably, the multi-segment lumped mass model and the electrical model are written in the same programming language.
Further preferably, the programming language is a Fortran language.
Compared with the prior art, the invention has the following beneficial technical effects:
the time domain simulation method is a main method for performing subsynchronous oscillation analysis on the steam turbine generator unit, currently, machine network coupling analysis software PSCAD/EMTDC is often adopted to carry out research, a shafting model for torsional oscillation analysis in PSCAD is a simple shafting model, and the simple shafting model cannot accurately simulate the coupling oscillation characteristics between machine networks and cannot obtain shafting local torsional oscillation response.
The invention discloses a subsynchronous oscillation joint simulation method based on a multi-section concentrated mass model, which expands a simple shafting model into the multi-section concentrated mass model, establishes an electric model of a generator and a power system and a data interface of the multi-section concentrated mass model and the electric model, realizes electromechanical joint simulation under the same simulation platform, greatly improves the calculation precision under torsional vibration response, can obtain torsional stress and torsional strain response of local details of the shafting, has higher simulation efficiency and more accurate result, has slower calculation speed and needs to occupy a large amount of calculation resources although the existing joint simulation model based on finite elements can also meet the precision requirement, and can greatly improve the calculation speed; and the shafting local stress strain information can be flexibly obtained, the mechanism of the sub-synchronous oscillation of the steam turbine generator unit is revealed, and the research on the sub-synchronous oscillation mechanism is favorably developed.
Furthermore, the multi-section centralized quality model and the electric model are written by adopting the same programming language, so that the data of the two models can be quickly communicated, and the calculation speed is increased.
Drawings
FIG. 1 is a flow chart of a sub-synchronous oscillation joint simulation method based on a multi-segment lumped mass model according to the invention;
FIG. 2 is a schematic structural diagram of a steam turbine generator unit shafting in the embodiment;
FIG. 3 is a multi-segment concentrated mass model diagram of the turbo generator set shafting established in the embodiment;
FIG. 4 is a graph of the change of the rotor speed of the generator with time obtained by the simulation method of the present invention and the conventional simulation method;
fig. 5 is a comparative enlarged view of fig. 4.
Detailed Description
The invention will now be described in further detail with reference to the drawings and specific examples, which are given by way of illustration and not by way of limitation.
Referring to fig. 1, a flow chart of the multistage centralized mass model-based subsynchronous oscillation joint simulation method of the present invention is shown, and a 600MW steam turbine in China is taken as an example for simulation, wherein the steam turbine is a subcritical, primary intermediate reheating, single-shaft, three-cylinder (high-middle pressing cylinder), four-steam-discharge and condensing steam turbine designed and produced by a harbin steam turbine plant, and has a unit model number of N600-16.7/538/538. The generator is a three-phase alternating-current hidden-pole type synchronous turbonator produced by Harbin motor factories, the cooling mode is water-hydrogen, namely the stator winding is cooled in water, and the rotor winding, the stator core and the end part structure are cooled by hydrogen, and the model is QFSN-600-2 YHG. The axis structure is shown in figure 2. The simulation method comprises the following specific steps:
step 1: acquiring structural parameters and material parameters of a shaft system of the steam turbine generator unit, establishing a multi-section centralized quality model of the shaft system as shown in figure 3, and being suitable for modeling the shaft system of a large steam turbine generator unit into dozens to hundreds of sections; forming a shafting model parameter initialization file;
step 2: constructing a torsional stress response calculation module of the torsional vibration fault lower shaft system according to the multi-section concentrated mass model established in the step 1, calculating the torsional stress response of the torsional vibration fault lower shaft system by using a transfer matrix method, wherein the algorithm adopts an increment form and is compiled by using Fortran language, the input of the torsional stress response calculation module is electromagnetic moment variation, and the output is the variation of the torsional angle, the angular velocity and the angular acceleration of each mass unit;
and step 3: acquiring various parameters of a steam turbine generator unit and an electric power system, and establishing an electric model of the generator and the electric power system by adopting PSCAD;
and 4, step 4: establishing a data exchange interface between the multi-section concentrated quality model in the step 1 and the electric model in the step 3, writing a calculation bottom program of the electric model by adopting a Fortran language, and writing a torsional stress response calculation module based on the Fortran language, so that the quick communication of the data of the two models can be realized;
and 5: initializing a data file Save _ te.temp for storing the electromagnetic torque of the previous step of simulation, and setting an initial value to be 0; initializing a file Save _ Res.temp for storing the torsional angle, the angular velocity and the angular acceleration under the torsional vibration action of each mass block in the previous step, and setting the initial value to be 0;
step 6: after receiving a project operation signal, starting an electrical model, obtaining a calculation result according to an initial value in the step 5, transmitting an electromagnetic torque value of the current step length to a multi-section concentrated mass model, reading an electromagnetic torque value of the previous step length stored in Save _ te.temp, subtracting the electromagnetic torque value of the previous step length stored in the Save _ te.temp to obtain an electromagnetic torque variation, transmitting the electromagnetic torque variation to a torsional stress response calculation module in the step 2, and storing the electromagnetic torque value of the step length in the Save _ te.temp to replace the original numerical value;
and 7: the torsional stress response calculation module calculates and obtains the torsional angle, the angular velocity and the variation of the angular acceleration of each mass unit in the multi-section concentrated mass model, adds the torsional angle, the angular velocity and the variation of the angular acceleration to the last step length response data stored in the Save _ res.temp to obtain the torsional angle, the angular velocity and the angular acceleration value of the current step length, and stores the result in the Save _ res.temp to replace the original data;
and 8: calculating the rotating speed of a mass unit where the generator rotor is located, transmitting the rotating speed to the electrical model, and starting the next long simulation; simultaneously outputting torsional stress and torsional strain results of each dangerous section;
and step 9: and (5) repeating the step 6 to the step 8 until the preset total simulation time is reached.
As shown in fig. 4 and 5, the torque value of 15s after a certain section of a shaft system is stabilized under the subsynchronous oscillation fault of a certain 600MW steam turbine generator unit is obtained through simulation, and compared with the conventional electrical simulation by using a simple mass block model, the torque value can be seen from the graph: the amplitude and the frequency of the result obtained by adopting the method are obviously different from those of the traditional simulation method, and the frequency of the result obtained by adopting the calculation method is closer to the second-order natural frequency of the shafting, so that the shafting oscillation characteristic can be accurately reflected.
When the PSCAD runs to a self-defined model, the model is compiled into a Fortran main program, and then the Fortran main program is called and executed through a Fortran compiler, so that the calculation of each step length of the PSCAD is serial and single-threaded, when the PSCAD runs to the self-defined model, the PSCAD enters the torsional vibration response calculation process of a shafting, the program is in a loop waiting state, the variables of other modules are kept unchanged, and after the self-defined model finishes the calculation and outputs related variables, the next long calculation is started, so that the joint simulation can be realized in an orderly and serial mode.
It should be noted that the above description is only a part of the embodiments of the present invention, and equivalent changes made to the system described in the present invention are included in the protection scope of the present invention. Persons skilled in the art to which this invention pertains may substitute similar alternatives for the specific embodiments described, all without departing from the scope of the invention as defined by the claims.
Claims (6)
1. A subsynchronous oscillation joint simulation method based on a multi-segment centralized quality model is characterized by comprising the following steps:
step 1: acquiring various parameters of a shaft system of a steam turbine generator unit, establishing a multi-section centralized quality model of the shaft system and forming a parameter initialization file of the shaft system model;
step 2: constructing a torsional stress response calculation module of the torsional vibration fault lower shaft system according to the multi-section concentrated mass model established in the step 1;
and step 3: acquiring various parameters of a steam turbine generator unit and an electric power system, and establishing electric models of a generator and the electric power system;
and 4, step 4: establishing a data exchange interface between the multi-section centralized quality model in the step 1 and the electric model in the step 3;
and 5: initializing an electromagnetic torque value in an electric model, and initializing torsion angles, angular velocities and angular acceleration values of each mass block in the multi-section concentrated mass model under the torsional vibration fault;
step 6: starting an electrical model, transmitting the electromagnetic torque value of the current step length obtained by calculation according to the initial value of the step 5 to a multi-section concentrated mass model, subtracting the electromagnetic torque value of the current step length from the electromagnetic torque value of the previous step length to obtain the electromagnetic torque variation and transmitting the electromagnetic torque variation to the torsional stress response calculation module of the step 2;
and 7: the torsional stress response calculation module calculates and obtains the variation of the torsional angle, the angular velocity and the angular acceleration of each mass unit in the multi-section concentrated mass model, and the variation is respectively added with the numerical value of the last step length to obtain the torsional angle, the angular velocity and the angular acceleration value of the current step length;
and 8: calculating the rotating speed of a mass unit where the generator rotor is located, transmitting the rotating speed to the electrical model, and starting the next long simulation; simultaneously outputting torsional stress and torsional strain results of each dangerous section;
and step 9: and (5) repeating the step 6 to the step 8 until the preset total simulation time is reached.
2. The multisegment lumped mass model-based subsynchronous oscillation joint simulation method of claim 1, wherein the parameters of the turbo generator set shafting comprise structural parameters and material parameters of the shafting.
3. The method for the subsynchronous oscillation joint simulation based on the multistage lumped mass model as recited in claim 1, wherein in the step 2, the torsional stress response calculation module is constructed by adopting an incremental form of an equation based on a transfer matrix method, and the input of the torsional stress response calculation module is the variation of the electromagnetic moment and the output is the variation of the torsional angle, the angular velocity and the angular acceleration of each mass unit.
4. The method for the subsynchronous oscillation joint simulation based on the multi-segment lumped mass model as recited in claim 1, wherein in the step 3, the electrical models of the generator and the power system are established by using PSCAD.
5. The method for subsynchronous oscillation joint simulation based on multi-segment lumped mass models as claimed in claim 1, wherein the multi-segment lumped mass models and the electrical models are written in the same programming language.
6. The method according to claim 5, wherein the programming language is Fortran's language.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010323069.3A CN111523229A (en) | 2020-04-22 | 2020-04-22 | Subsynchronous oscillation joint simulation method based on multi-segment centralized quality model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010323069.3A CN111523229A (en) | 2020-04-22 | 2020-04-22 | Subsynchronous oscillation joint simulation method based on multi-segment centralized quality model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111523229A true CN111523229A (en) | 2020-08-11 |
Family
ID=71910446
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010323069.3A Pending CN111523229A (en) | 2020-04-22 | 2020-04-22 | Subsynchronous oscillation joint simulation method based on multi-segment centralized quality model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111523229A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113779697A (en) * | 2021-09-13 | 2021-12-10 | 中国华能集团清洁能源技术研究院有限公司 | Turbogenerator unit torsional vibration joint simulation method and device and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104236705A (en) * | 2014-09-19 | 2014-12-24 | 华北电力大学 | Method and system for monitoring large disturbance transient impact torsional vibration of turbine-generator shaft system |
| CN104361172A (en) * | 2014-11-17 | 2015-02-18 | 清华大学 | Subsynchronous resonance electro-mechanical combination simulation method and system |
-
2020
- 2020-04-22 CN CN202010323069.3A patent/CN111523229A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104236705A (en) * | 2014-09-19 | 2014-12-24 | 华北电力大学 | Method and system for monitoring large disturbance transient impact torsional vibration of turbine-generator shaft system |
| CN104361172A (en) * | 2014-11-17 | 2015-02-18 | 清华大学 | Subsynchronous resonance electro-mechanical combination simulation method and system |
Non-Patent Citations (1)
| Title |
|---|
| 赵鹏程: "汽轮发电机组轴系扭振机理及安全性分析", 《中国优秀博硕士学位论文全文数据库(博士)》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113779697A (en) * | 2021-09-13 | 2021-12-10 | 中国华能集团清洁能源技术研究院有限公司 | Turbogenerator unit torsional vibration joint simulation method and device and storage medium |
| CN113779697B (en) * | 2021-09-13 | 2022-12-13 | 中国华能集团清洁能源技术研究院有限公司 | Turbogenerator unit torsional vibration joint simulation method and device and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107086606B (en) | Method for identifying load model parameters of comprehensive equivalent asynchronous motor of power distribution network | |
| CN106599526B (en) | Impedance analysis method for sub-synchronous oscillation of thermal power generating unit | |
| CN109766586A (en) | A kind of method and system automatically generating large scale electric network electromagnetic transient simulation model | |
| CN103678900B (en) | A kind of network decoupling calculation method for regional distribution network real-time simulation | |
| CN103810646B (en) | Improved projection integral algorithm based active power distribution system dynamic simulation method | |
| CN110854884B (en) | Wind power collection area subsynchronous oscillation risk online assessment and early warning method | |
| CN110705078A (en) | Simulation model development system and analysis method of nuclear power device | |
| CN108899895A (en) | A kind of electromechanics-electromagnetism the hybrid simulation method and system of the electric system of UPFC containing THE UPFC | |
| CN111597680A (en) | Subsynchronous resonance risk quantitative evaluation method based on electrical damping at modal frequency | |
| CN108418242B (en) | Doubly-fed wind turbine dynamic equivalence method based on similarity coherence | |
| CN113937764A (en) | Low-voltage distribution network high-frequency measurement data processing and topology identification method | |
| CN111523229A (en) | Subsynchronous oscillation joint simulation method based on multi-segment centralized quality model | |
| CN112307629A (en) | Comprehensive energy multi-energy flow unified real-time simulation system and method based on Thermolib software package | |
| CN108549758A (en) | A kind of modularization multi-level converter real-time emulation method suitable for FPGA | |
| CN115912330A (en) | Two-stage opportunity constraint optimization method and system of active power distribution network topology evolution model | |
| CN105529736B (en) | Crow bar puts into double-fed fan motor group of planes fault current Equivalent Model computational methods | |
| CN112072692A (en) | Impedance equivalence method and device for new energy power generation station | |
| CN111740449B (en) | Method for analyzing damping characteristics of wind turbine group to synchronous generator set | |
| CN106991229B (en) | Wind power plant equivalent modeling method for complex topology | |
| CN104361172A (en) | Subsynchronous resonance electro-mechanical combination simulation method and system | |
| CN116865318A (en) | Power transmission network and energy storage joint planning method and system based on two-stage random optimization | |
| CN114298478A (en) | Small disturbance stability identification method and system for wind power grid-connected system | |
| CN110489801B (en) | Generator shaft system multi-mass-block parameter simplification method by utilizing particle swarm optimization algorithm | |
| CN112688321B (en) | Large power grid subsynchronous oscillation path acquisition method and system | |
| CN101340121A (en) | Transient simulation simplified model for internal fault of generator, modeling method and application thereof |
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: 20200811 |