CN110489801A - Utilize the generator shafting multiple mass parameter predigesting method of particle swarm optimization algorithm - Google Patents
Utilize the generator shafting multiple mass parameter predigesting method of particle swarm optimization algorithm Download PDFInfo
- Publication number
- CN110489801A CN110489801A CN201910653770.9A CN201910653770A CN110489801A CN 110489801 A CN110489801 A CN 110489801A CN 201910653770 A CN201910653770 A CN 201910653770A CN 110489801 A CN110489801 A CN 110489801A
- Authority
- CN
- China
- Prior art keywords
- mass
- frequency
- shafting
- formula
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000002245 particle Substances 0.000 title claims abstract description 34
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 27
- 238000005457 optimization Methods 0.000 title claims abstract description 23
- 230000010355 oscillation Effects 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 230000000717 retained effect Effects 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000001351 cycling effect Effects 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 10
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000005452 bending Methods 0.000 description 5
- 239000012141 concentrate Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Biophysics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Primary Health Care (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a kind of generator shafting multiple mass parameter predigesting methods using particle swarm optimization algorithm for belonging to power system modeling and analysis field.Specifically relating to the analysis of fired power generating unit sub-synchronous oscillation carries out parameter predigesting with shafting multiple mass model, according to gradually generator mass is merged, mass section where judging subsynchronous frequency mode, and remaining mass is merged, and then obtain the simplified model without high frequency mode component, and model parameter is corrected based on particle swarm optimization algorithm, the shafting multiple mass model and simulation calculation being simplified, the foundation of multiple mass model in power system simulation software is simplified, and is able to reflect the dynamic characteristic of electric system fired power generating unit sub-synchronous oscillation.Make the high fdrequency component for not including shafting torsional oscillation in simplified shafting multiple mass model.The confidence level of Efficiency and result of study will be greatly improved.It can satisfy the requirement of electric system fired power generating unit sub-synchronous oscillation specificity analysis.
Description
Technical field
The invention belongs to power system modeling and analysis field, in particular to a kind of power generation using particle swarm optimization algorithm
Arbor system multiple mass parameter predigesting method, it is specifically simple with particle swarm optimization algorithm adjustment shafting multiple mass based on gradually merging
Change method.
Background technique
Generator shafting multiple mass model determines the natural torsion frequency of shafting, therefore establishes reasonable shaft system of unit mould
Type is to study the basis of sub-synchronous oscillation problem.And the shafting length of large sized unit is very long in practice, and is continuous unlimited
First shafting model, is generally not easy to analyze.It mainly includes simply concentrating that shaft system of unit torsional oscillation inherent characteristic, which calculates the model used,
Model, multistage lumped-mass model and continuum mass model.Continuum mass model computational solution precision is high, although can be accurately anti-
Higher-order torsional vibration characteristic and mechanical characteristics are reflected, but calculation amount is larger, the main influence for calculating shafting partial structurtes, for electromechanical coupling
The dynamic analysis of collaboration system, this method are simultaneously not suitable for.Multistage concentrates mass model calculation amount smaller compared to the former, calculates essence
Degree it is also ensured that, but shafting is divided into tens sections to several hundred sections by multistage lumped-mass model in practice, it is contemplated that the tool of shafting
Body details, explication de texte is crucial with non-key section, and calculation amount is still very huge.It is simple to concentrate mass quality model simplification side
Method has good effect when carrying out the qualitative analysis of shaft system of unit torsional oscillation inherent characteristic, although simplified simple concentration matter
Block models can bring some errors, but as long as the subsynchronous oscillation frequency for guaranteeing that the two calculates is consistent, and the simplified vibration shape
Overall trend it is consistent, that ensures that simplified reasonability.Some documents, which are computed, to be pointed out, certain adjacent rigid wheel discs close
The torsional vibration characteristic of model is not influenced substantially after and, therefore shaft system of unit modelling is letter by cylinder position by many researchers
It is single to concentrate mass model.In the research for the subsynchronous problem of fired power generating unit, the multimass block shafting torsional oscillation mould of producer's offer
It may include high fdrequency component in type, however only focus on the torsional mode of subsynchronous frequency component in sub-synchronous oscillation research, therefore,
It is calculated to simplify Analysis of Torsional Vibration, then its further progress is effectively merged, obtained low without multiple mass model high frequency mode
The simple mass model of rank.In the analysis for carrying out sub-synchronous oscillation problem, shafting design frequency and practical frequency that power plant provides
Generally there is small gap, along with the model is simplified Equivalent Model, calculated nature torsion frequency may be with
There are deviations for field measurement frequency, it is therefore desirable to simplified shafting model parameter is modified, from the point of view of in practice, and inertia
Time constant depends mainly on the quality of shafting mass, is easy to accurately obtain its value, and precision is higher, therefore corrects each mass block
Between stiffness coefficient more rationally, after corrected parameter also need carry out the vibration shape correction.
Summary of the invention
The purpose of the present invention is to propose to a kind of generator shafting multiple mass parameter predigesting sides using particle swarm optimization algorithm
Method；It is characterized in that, the generator shafting multiple mass parameter predigesting method is based on gradually Unite principle, by Steam Turbine axis
It is that multiple mass shafting model is further simplified as method containing only subsynchronous oscillation frequency component shafting equivalent model, and utilizes grain
Swarm optimization carries out parameter correction, and detailed process includes:
(1) the stiffness coefficient K between the inertia time constant M and each mass block of generator multiple mass is inputted；
(2) gradually generator multiple mass is merged, between the inertia time constant M and each mass block after merging
Stiffness coefficient K calculates new torsion frequency, and is compared with the subsynchronous oscillation frequency component of reservation；
(3) it obtains containing only subsynchronous oscillation frequency component shafting equivalent model, and calculates the inertia time constant M after merging
And the stiffness coefficient K between each mass block；
(4) stiffness coefficient between each mass block is corrected using particle swarm algorithm；
(5) meet precision:
If meeting required precision, then next step (6) are carried out；If it is not, then return step (4)；
(6) parameter after output optimization.
Stiffness coefficient K between the inertia time constant M and each mass block of step (1) the input generator multiple mass,
The torsional vibration frequency of Steam Turbine refers to that unit after being disturbed, produces between each mass block of rotor of steam turbo generator shafting
Raw natural frequency, they are only related with the inherent characteristic of mechanical system；If if therefore two neighboring mass will be merged,
It can then reduce by a natural frequency, section where subsynchronous frequency component is obtained by the value that judgement reduces natural frequency
Face；Wherein, the inertia time constant of new mass is to merge the sum of the former two, and the stiffness coefficient between each mass block after merging is pressed
Pro rate, specific formula for calculation are as follows:
Mj=Mi+Mi+1 (1-1)
In above formula, M indicates mass inertia time constant, and K indicates the stiffness coefficient parameter between each mass block, and i indicates matter
Mass positional number before merged block, j indicate the mass positional number after merging, i.e. M1For the inertia time constant of first mass, K1
For the stiffness coefficient between mass 1 and mass 2.
The subsynchronous frequency component that reservation can be selected in the torsion frequency of the offer provided by power plant, is denoted as:
Rest_f=[rest_f1, rest_f2..., rest_fk] (1-4)
In formula, rest_f is the subsynchronous frequency retained, rest_fkFor k-th of frequency of reservation.
The step (2), (3) are specially to merge two adjacent mass, based on formula (1-1), (1-2), (1-3)
Calculate the inertia time constant after merging and the stiffness coefficient between each mass block.When two neighboring mass mass merges,
Natural frequency after then merging, note:
F=[f1,f2,…,fh] (1-5)
In above formula, fhFor h-th of frequency of reservation.
The frequency component that each of formula (1-4) retains is made the difference with h frequency in (1-5) respectively, and takes it absolutely
To the minimum value of value, it is denoted as T:
T=[T1,T2,…,Tk] (1-6)
Wherein TkAre as follows:
Tk=min | rest_fk-f1|,|rest_fk-f2|,…,|rest_fk-fh|} (1-7)
Rest_f in above formulakK-th of the frequency component as retained, fhSame formula (1-5).
At this point, if TkValue be then to represent the frequency component close to zero and be not merged, conversely, then representing the frequency point
Amount is present between the adjacent mass, and then determines that the position needs to be retained；And according to the maximum value in T whether close to zero
Come judge merge after mass between whether comprising think reservation low frequency component, if the maximum value of T is not wrapped close to zero
Contain；
Here note mass quantity is n, then it is n-1 times total to merge number, then can accordingly obtain n-1 maxT, be denoted as S:
S=[S1,S2,…,Sn-1] (1-8)
S in above formulan-1Represent the maximum value corresponded in formula (1-6) after merging each time.
It can determine retained subsynchronous frequency component institute according to the position of the maximum value of k reserve frequency number preceding in S
Position, and remaining mass is merged, while calculating corresponding inertia time constant and each matter by formula 1-1~1-3
Stiffness coefficient between gauge block.
The step (4) is for it is excellent to be now based on population in view of the parameter after tentatively merging can generate certain error
Change algorithm to be corrected parameter, it is contemplated that the value of inertia time constant depends mainly on the quality of shafting mass, is easy to calibrated
Its value really is obtained, therefore the stiffness coefficient corrected between each mass block is more reasonable, therefore retains the value of inertia time constant, correction
Stiffness coefficient between each mass block accordingly checks equivalent shafting model, to guarantee the consistency of its vibration shape trend；
In iterative process each time, particle updates speed and the position of itself by individual extreme value and group's extreme value, more
New formula is as follows:
In formula, meaning representated by each coefficient is as follows:
--- the d of the t times iteration particle e flying speed ties up component；
--- the suitable d in the t times position iteration particle e ties up component；
c1,c2--- acceleration constant adjusts study maximum step-length；
r1,r2--- two random functions, value range [0,1], to increase search randomness；
ω --- inertia weight, nonnegative number adjust the search range to solution space；
For the rapidity for guaranteeing search, formula 1-11 is set by inertia weight ω, guarantees that inertia weight successively decreases, can be protected in this way
Card reduces search range when being optimal near solution, (as shown in Figure 3).
In formula, ωstart=0.9, ωend=0.4, maxgen are population scale size.
The concrete operations of the step (5) are with practical frequency and to calculate the sum of the absolute value of frequency difference as fitness letter
Numerical value, i.e.,
Fitness=| f1-real_f1|+|f2-real_f2|+…+|fr-real_fs| (1-12)
In formula, real_fsFor practical frequency value, frTo calculate frequency.
Cycling condition is jumped out in setting:
Fitness=ε < 10-4 (1-13)
It is to follow the process of current optimal solution that particle swarm algorithm, which entirely searches for renewal process, therefore, under normal circumstances, the calculation
The convergence rate of method quickly, is reaching maximum number of iterations, maxgen=50 is arranged；Or meets to jump out after formula (1-13) and follow
Ring obtains current optimal solution, i.e., the stiffness coefficient between each mass block.
The beneficial effects of the invention are as follows the present invention according to analysis shafting parameter in electric system fired power generating unit sub-synchronous oscillation
The requirement of stability analysis proposes and is further simplified as the mass model containing high fdrequency component to contain only subsynchronous oscillation frequency
The shafting model of rate component, and parameter correction is carried out to simplified model based on particle swarm optimization algorithm, simplify power train
The shafting model foundation of system simulation software PSCAD has obtained the simplification turbine shafting equivalent model without high fdrequency component and has imitated
It is true to calculate, it can satisfy the requirement of electric system fired power generating unit sub-synchronous oscillation specificity analysis.
Detailed description of the invention
Fig. 1 is to simplify mass flow chart.
Fig. 2 is that mass merges schematic diagram.
Fig. 3 is the relation curve of inertia weight and the number of iterations.
Fig. 4 is typical large turbine-generator set shafting structure schematic diagram.
Fig. 5 is that mass tentatively merges flow chart.
Fig. 6 is that particle swarm optimization algorithm corrects the stiffness coefficient flow chart between each mass block.
Fig. 7 is northern two power plant's axis multiple mass tying compositions.
Fig. 8-1 provides 38.30Hz mode bending vibation mode picture for northern two power plant.
Fig. 8-2 is 38.30Hz mode bending vibation mode picture after simplifying.
Fig. 9-1 provides 29.91Hz mode bending vibation mode picture for northern two power plant.
Fig. 9-2 is 29.91Hz mode bending vibation mode picture after simplifying.
Figure 10-1 provides the 17.49Hz mode vibration shape for northern two power plant.
Figure 10-2 is 17.49Hz mode bending vibation mode picture after simplifying.
Specific embodiment
The invention proposes a kind of generator shafting multiple mass parameter predigesting methods.The present invention is given with reference to the accompanying drawing
It further illustrates.
It is as shown in Figure 1 simplified mass flow chart.The generator shafting multiple mass parameter predigesting method is based on gradually
Balancing of Steam Turbine Shaft multiple mass shafting model is further simplified as containing only subsynchronous oscillation frequency component shafting etc. by Unite principle
The method for imitating model, and parameter correction is carried out using particle swarm algorithm, detailed process includes
(1) the stiffness coefficient K between the inertia time constant M and each mass block of generator multiple mass is inputted；
(2) gradually generator multiple mass is merged, between the inertia time constant M and each mass block after merging
Stiffness coefficient K calculates new torsion frequency, and is compared with the subsynchronous oscillation frequency component of reservation；
(3) it obtains containing only subsynchronous oscillation frequency component shafting equivalent model, and calculates the inertia time constant M after merging
And the stiffness coefficient K between each mass block；
(4) the stiffness coefficient K between each mass block is corrected using particle swarm algorithm；
(5) meet precision:
If meeting required precision, then next step (6) are carried out；If it is not, then return step (4)；
(7) parameter after output optimization.
Stiffness coefficient K between the inertia time constant M and each mass block of step (1) the input generator multiple mass,
Typical large turbine-generator set shafting structure schematic diagram (as shown in Figure 4), the torsional vibration frequency of Steam Turbine refer to that unit exists
After being disturbed, the natural frequency generated between each mass block of rotor of steam turbo generator shafting, they only and mechanical system
Inherent characteristic it is related；If two neighboring mass is therefore merged (as shown in Figure 2), it can reduce by a natural frequency,
The section where subsynchronous frequency component is obtained by the value that judgement reduces natural frequency；Wherein, the inertial time of new mass
Between constant be to merge the sum of the former two, the stiffness coefficient proportional assignment between each mass block after merging, specific formula for calculation
Are as follows:
Mj=Mi+Mi+1 (1-1)
In above formula, M indicates mass inertia time constant, and K indicates the stiffness coefficient parameter between each mass block, and i indicates matter
Mass positional number before merged block, j indicate the mass positional number after merging, (i.e. M1For the inertia time constant of first mass,
K1For the stiffness coefficient between mass 1 and mass 2).
The subsynchronous frequency component that reservation can be selected in the torsion frequency of the offer provided by power plant, is denoted as:
Rest_f=[rest_f1, rest_f2..., rest_fk] (1-4)
In formula, rest_f is the subsynchronous frequency retained, rest_fkFor k-th of frequency of reservation.
The step (2), (3) are specially to merge two adjacent mass, based on formula (1-1), (1-2), (1-3)
Calculate the shafting shafting inertia time constant after merging and the stiffness coefficient between each mass block.When two neighboring mass merges
When, then the natural frequency after merging, note:
F=[f1,f2,…,fh] (1-5)
In above formula, fhFor h-th of frequency of reservation.
The frequency component that each of formula (1-4) retains is made the difference with h frequency in (1-5) respectively, and takes it absolutely
To the minimum value of value, it is denoted as T:
T=[T1,T2,…,Tk] (1-6)
T in above formulakAre as follows:
Tk=min | rest_fk-f1|,|rest_fk-f2|,…,|rest_fk-fh|} (1-7)
Rest_f in above formulakK-th of the frequency component as retained, fhSame formula (1-5).
At this point, if TkValue be then to represent the frequency component close to zero and be not merged, conversely, then representing the frequency point
Amount is present between the adjacent mass, and then determines that the position needs to be retained；And according to the maximum value in T whether close to zero
Come judge merge after mass between whether comprising think reservation low frequency component, if the maximum value of T is not wrapped close to zero
Contain；
Here note mass quantity is n, then it is n-1 times total to merge number, then can accordingly obtain n-1 maxT, be denoted as S:
S=[S1,S2,…,Sn-1] (1-8)
S in above formulan-1Represent the maximum value corresponded in formula (1-6) after merging each time.
It can determine retained subsynchronous frequency component according to the position of k preceding in S (number of reserve frequency) maximum values
The position at place, and remaining mass is merged, while by the corresponding inertia time constant of formula 1-1~1-3 calculating and respectively
Stiffness coefficient between mass block.
The step (4) is for it is excellent to be now based on population in view of the parameter after tentatively merging can generate certain error
Change algorithm to be corrected parameter, it is contemplated that the value of inertia time constant depends mainly on the quality of shafting mass, is easy to calibrated
Its value really is obtained, therefore the stiffness coefficient corrected between each mass block is more reasonable, therefore retains the value of inertia time constant, correction
Stiffness coefficient between each mass block accordingly checks equivalent shafting model, to guarantee the consistency of its vibration shape trend；
In iterative process each time, particle updates speed and the position of itself by individual extreme value and group's extreme value, more
New formula is as follows:
In formula, meaning representated by each coefficient is as follows:
--- the d of the t times iteration particle e flying speed ties up component；
--- the suitable d in the t times position iteration particle e ties up component；
c1,c2--- acceleration constant adjusts study maximum step-length；
r1,r2--- two random functions, value range [0,1], to increase search randomness；
ω --- inertia weight, nonnegative number adjust the search range to solution space；
For the rapidity for guaranteeing search, formula 1-11 is set by inertia weight ω, guarantees that inertia weight successively decreases, can be protected in this way
Card reduces search range (as shown in Figure 3) when being optimal near solution.
In formula, ωstart=0.9, ωend=0.4, maxgen are population scale size.
The concrete operations of the step (5) are with practical frequency and to calculate the sum of the absolute value of frequency difference as fitness letter
Numerical value, i.e.,
Fitness=| f1-real_f1|+|f2-real_f2|+…+|fr-real_fs| (1-12)
In formula, real_fsFor practical frequency value, frTo calculate frequency.
Cycling condition is jumped out in setting:
Fitness=ε < 10-4 (1-13)
It is to follow the process of current optimal solution that particle swarm algorithm, which entirely searches for renewal process, therefore, under normal circumstances, the calculation
The convergence rate of method quickly, is reaching maximum number of iterations, maxgen=50 is arranged；Or circulation is jumped out after meeting formula 1-13,
Obtain current optimal solution, i.e., the stiffness coefficient between each mass block.
The simplification turbine shafting equivalent model without high fdrequency component has been obtained according to above method, can be used for simplifying
Generator shafting multiple mass model foundation and simulation calculation in PSCAD.
Typical large turbine-generator set shafting structure schematic diagram according to Fig.4, column are write shaft system of unit segmentation and are concentrated
Mass-spring model is mathematically represented as
In formula, meaning representated by each coefficient is as follows.(point above variable represents the differential equation)
δ is electrical torsional angle vector:
δ=[δ1 δ2 … δN]Τ
Wherein, δNElectrical torsional angle for n-th mass relative to synchronous rotary reference axis, unit are radian (rad).
ω is electrical angular speed: ω=[ω1 ω2 … ωN]Τ, wherein ωNIt is single for the electrical angular speed of n-th mass
Position is radian per second (rad/s).
Δ T is the difference of machine torque and electromagnetic torque.
M is inertia time constant matrix: M=diag [M1 M2 … MN]；Wherein, MNFor the inertia time of n-th mass
Constant, unit are kilogram-meter2(kg·m2)。
Stiffness coefficient of the K between each mass block, i.e. elastic coefficient matrix:
Wherein, Kp,p+1Stiffness coefficient between mass, n are mass number, and unit is ox rice/radian (Nm/rad).
D is damped coefficient matrix:
Wherein, DqqWith Dq-1,qSelf-damping and mutual damping coefficient between steam turbine mass, n are mass number, and unit is ox
Rice second/radian (Nms/rad).
By calculating M-1The available mode frequency of the characteristic value of K.
It is illustrated in figure 5 mass and tentatively merges flow chart.Include:
1) selection needs the subsynchronous frequency component retained, and is calculated as follows:
The parameter of 1.1 input generator multiple mass: the stiffness coefficient between inertia time constant and each mass block；
1.2 merge two adjacent mass, the inertia time after merging by the method in formula 1-1~1-3
Stiffness coefficient parameter between constant and each mass block calculates；
Parameter after merging is carried out mode frequency calculating by 1.3；
1.4 will need the frequency in the subsynchronous frequency component and 1.3 retained to make the difference respectively, and minimize, i.e.,
Tk=min | rest_fk-f1|,|rest_fk-f2|,…,|rest_fk-fh|} (1-15)
1.5 seek T=[T1,T2,…,Tk] maximum value, be denoted as S；
1.6 remember that mass quantity is that n obtains S=[S after n-1 times merges here1S2…Sn-1], by judging in S preceding k
The position of (number of reserve frequency) maximum value judges the section where subsynchronous frequency, and by the method in formula 1-1~1-3
Position without subsynchronous component is merged, the stiffness coefficient between preliminary inertia time constant and each mass block is calculated；
2) parameter correction (as shown in Figure 6) is carried out based on particle swarm optimization algorithm:
2.1, according to the stiffness coefficient in step (3) between obtained inertia time constant and each mass block, initialize excellent
Change the value (± 25% variation) of the stiffness coefficient between each mass block, setting Population Size is 100, and the number of iterations is 50 times, and
Initialize iteration speed weights omega；
For fitness function value, i.e., 2.2 with practical frequency and calculate the sum of absolute value of frequency difference
Fitness=| f1-real_f1|+|f2-real_f2|+…+|fr-real_fs| (1-16)
In formula, real_fsFor practical frequency value, frTo calculate frequency.
2.3 calculate individual extreme value and group's extreme value optimal solution at present, i.e., the stiffness coefficient between each mass block；
2.4 are iterated, and jump out circulation after meeting fitness function value condition or after reaching maximum number of iterations；
2.5 export optimal solutions, i.e., the stiffness coefficient between each mass block.
Embodiment
In order to verify the correctness of the above shafting equivalent model, it is directed at two power plant's shafting model of north of eastern direct current sending end power grid
Simplified, carries out vibration shape verifying after correcting the stiffness coefficient between each mass block according to practical frequency, guarantee simplified model
Reasonability.
The shafting structure figure that two power plant of north as shown in Figure 7 provides, inertia time constant that host manufacturer provides and
Stiffness coefficient between each mass block, as shown in table 1.
The northern two METHOD FOR TURBOGENERATOR SET shafting parameter models of table 1
Northern two power plant concentration can be calculated according to above-mentioned offer data, and mass --- the natural mode shape of spring model is shown in
Shown in table 2.
Northern two power plant of table 2 concentrate mass --- the natural mode shape of spring model
1) preliminary simplified model is obtained
From the point of view of intrinsic frequency, the subsynchronous frequency mode component that should retain is 17.49Hz, 29.91Hz and
38.30Hz because the multiple mass that the power plant provides shares 9, therefore gradually carries out 8 merging, time that will be retained respectively respectively
Synchronizing frequency component makes the difference with 7 frequencies that 8 times merge, and the S in step (2) can be obtained, and the result is shown in shown in table 3, is easy
The position for obtaining first three maximum value is 3,4,5, therefore mass 1,2,3 is merged into a mass, and mass 4,5 is respectively as one
A independent mass, 6,7,8 merge into the 4th mass, and merging method is detailed in summary of the invention, inertia time constant after merging, each
Stiffness coefficient and natural torsion frequency between mass block are shown in Table 4.
Table 3 judges that preliminary mass merges foundation according to the size of S value
Inertia time constant, stiffness coefficient and natural torsion frequency after the tentatively merging of table 4
2) stiffness coefficient between each mass block is corrected
The corresponding natural torsion frequency of mass model and realistic model after preliminary merging has certain deviation, and gained is each
Stiffness coefficient between mass block is optimized with particle swarm optimization algorithm, and the parameter and torsion frequency after optimization are shown in Table 5 institutes
Show.
The stiffness coefficient and torsion frequency between each mass block after the optimization of table 5
3) model corrects
To prove simplified model validity, because the vibration shape of simple mass model is by the line segment structure of a few point connection
At, and the vibration shape of multiple mass is the curve of relative smooth, so the former can regard the simplification curve of the latter as, the entirety of the two becomes
Gesture should be consistent, and the vibration shape of the three subsynchronous frequency components just retained separately below provides the vibration shape with power plant and is compared
(as shown in Fig. 8-1,8-2,9-1,9-2,10-1,10-2), it can be clearly seen that the vibration shape trend for simplifying front and back is consistent, to demonstrate,prove
The validity of simplified model is illustrated.
Claims (5)
1. a kind of generator shafting multiple mass parameter predigesting method using particle swarm optimization algorithm is it is characterized in that, the power generation
Arbor system multiple mass parameter predigesting method is to be based on gradually Unite principle, and Balancing of Steam Turbine Shaft multiple mass shafting model is further
It is reduced to the method containing only subsynchronous oscillation frequency component shafting equivalent model, and carries out parameter correction using particle swarm algorithm,
Detailed process includes:
(1) parameter of generator multiple mass: the stiffness coefficient K between inertia time constant M and each mass block is inputted；
(2) gradually generator multiple mass is merged, with the rigidity between the inertia time constant M and each mass block after merging
COEFFICIENT K calculates new torsion frequency, and is compared with the subsynchronous oscillation frequency component of reservation；
(3) it obtains containing only subsynchronous oscillation frequency component shafting equivalent model, and calculates the inertia time constant M after merging and each
Stiffness coefficient K between mass block；
(4) the stiffness coefficient K between mass block is corrected using particle swarm algorithm；
(5) meet precision:
If meeting required precision, then next step (6) are carried out；If it is not, then return step (4)；
(6) parameter after output optimization.
2. utilizing its feature of the generator shafting multiple mass parameter predigesting method of particle swarm optimization algorithm according to claim 1
It is, the stiffness coefficient K between the inertia time constant M and each mass block of step (1) the input generator multiple mass, steamer
The torsional vibration frequency of unit refers to unit after being disturbed, generated between each mass of rotor of steam turbo generator shafting from
Right frequency, they are only related with the inherent characteristic of mechanical system；If therefore merging two neighboring mass, one can be reduced
A natural frequency obtains the section where subsynchronous frequency component by the value that judgement reduces natural frequency；Wherein, new matter
The inertia time constant of block is to merge the sum of the former two, the stiffness coefficient proportional assignment between each mass block after merging, tool
Body calculation formula are as follows:
Mj=Mi+Mi+1 (1-1)
In above formula, M indicates the inertia time constant of mass, and K indicates that the stiffness coefficient between each mass block, i indicate that mass merges
Preceding mass positional number, j indicate the mass positional number after merging, M1For the inertia time constant of first mass, K1For mass 1
Stiffness coefficient between mass 2；
The subsynchronous frequency component that reservation can be selected in the torsion frequency of the offer provided by power plant, is denoted as:
Rest_f=[rest_f1, rest_f2..., rest_fk] (1-4)
In formula, rest_f is the subsynchronous frequency retained, rest_fkFor k-th of frequency of reservation.
3. utilizing its feature of the generator shafting multiple mass parameter predigesting method of particle swarm optimization algorithm according to claim 1
It is, the step (2), (3) are specially to merge two adjacent mass, are calculated by formula (1-1), (1-2), (1-3)
The stiffness coefficient between shafting inertia time constant and each mass block after merging；When two neighboring mass merges, then
Natural frequency after merging, note:
F=[f1,f2,…,fh] (1-5)
In above formula, fhFor h-th of frequency of reservation；
The frequency component that each of formula (1-4) retains is made the difference with h frequency in (1-5) respectively, and takes its absolute value
Minimum value, be denoted as T:
T=[T1,T2,…,Tk] (1-6)
T in above formulakAre as follows:
Tk=min | rest_fk-f1|,|rest_fk-f2|,…,|rest_fk-fh| rest_f in (1-7) above formulakI.e.
For k-th of frequency component of reservation, fhSame formula (1-5)；
At this point, if TkValue be close to zero, then represent the frequency component and be not merged, conversely, then represent the frequency component presence
Between the adjacent mass, and then determine that the position needs to be retained；And whether judged close to zero according to the maximum value in T
It whether include the low frequency component for wanting to retain between mass after merging, if the maximum value of T does not include close to zero and wants to protect
The low frequency component stayed；Here note mass quantity is n, then it is n-1 times total to merge number, then can accordingly obtain n-1 maxT, be denoted as
S:
S=[S1,S2,…,Sn-1] (1-8)
S in above formulan-1Represent the maximum value corresponded in formula (1-6) after merging each time；
It is where can determine retained subsynchronous frequency component according to the position of the maximum value of k reserve frequency number preceding in S
Position, and remaining mass is merged, while calculating corresponding inertia time constant and each mass block by formula 1-1~1-3
Between stiffness coefficient.
4. utilizing its feature of the generator shafting multiple mass parameter predigesting method of particle swarm optimization algorithm according to claim 1
It is, the step (4) is to be now based on Particle Swarm Optimization in view of the parameter after tentatively merging can generate certain error
Method is corrected parameter, it is contemplated that the value of inertia time constant depends mainly on the quality of shafting mass, is easy to accurately
Its value is obtained, therefore the stiffness coefficient corrected between each mass block is more reasonable, therefore retains the value of inertia time constant, corrects each matter
Stiffness coefficient between gauge block accordingly checks equivalent shafting model, to guarantee the consistency of its vibration shape trend；
In iterative process each time, particle updates speed and the position of itself by individual extreme value and group's extreme value, updates public
Formula is as follows:
In formula, meaning representated by each coefficient is as follows:
--- the d of the t times iteration particle e flying speed ties up component；
--- the suitable d in the t times position iteration particle e ties up component；
c1,c2--- acceleration constant adjusts study maximum step-length；
r1,r2--- two random functions, value range [0,1], to increase search randomness；
ω --- inertia weight, nonnegative number adjust the search range to solution space；
For the rapidity for guaranteeing search, formula 1-11 is set by inertia weight ω, guarantees that inertia weight successively decreases, can guarantee in this way
Search range is reduced when being optimal near solution,
In formula, ωstart=0.9, ωend=0.4, maxgen are population scale size.
5. the generator shafting multiple mass parameter predigesting method of particle swarm optimization algorithm is utilized according to claim 1, it is special
Sign is that the concrete operations of the step (5) are with practical frequency and to calculate the sum of the absolute value of frequency difference as fitness letter
Numerical value, i.e.,
Fitness=| f1-real_f1|+|f2-real_f2|+…+|fr-real_fs| (1-12)
In formula, real_fsFor practical frequency value, frTo calculate frequency；
Cycling condition is jumped out in setting:
ε < 10 fitness=-4 (1-13)
It is to follow the process of current optimal solution that particle swarm algorithm, which entirely searches for renewal process, therefore, under normal circumstances, the algorithm
Convergence rate quickly, is reaching maximum number of iterations, maxgen=50 is arranged；Or circulation is jumped out after meeting formula 1-13, it obtains
Current optimal solution, i.e., the stiffness coefficient between each mass block.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910653770.9A CN110489801B (en) | 2019-07-19 | 2019-07-19 | Generator shaft system multi-mass-block parameter simplification method by utilizing particle swarm optimization algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910653770.9A CN110489801B (en) | 2019-07-19 | 2019-07-19 | Generator shaft system multi-mass-block parameter simplification method by utilizing particle swarm optimization algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110489801A true CN110489801A (en) | 2019-11-22 |
| CN110489801B CN110489801B (en) | 2022-10-21 |
Family
ID=68547470
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910653770.9A Active CN110489801B (en) | 2019-07-19 | 2019-07-19 | Generator shaft system multi-mass-block parameter simplification method by utilizing particle swarm optimization algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110489801B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111597680A (en) * | 2020-04-09 | 2020-08-28 | 华北电力大学 | Subsynchronous resonance risk quantitative evaluation method based on electrical damping at modal frequency |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101515719A (en) * | 2009-04-02 | 2009-08-26 | 华北电力大学 | Method for establishing generator shafting multimass block damping variable model and device thereof |
| CN101807223A (en) * | 2010-03-16 | 2010-08-18 | 浙江大学 | Method for simulating mechanical damping of steam turbine shafting of power system |
| CN103336874A (en) * | 2013-07-15 | 2013-10-02 | 国家电网公司 | Electrical power system subsynchronous oscillation on-line analyzing and early-warning method based on time-domain simulation |
| CN103605819A (en) * | 2013-09-05 | 2014-02-26 | 昆明理工大学 | Simulation method for water turbine generator set shafting vibration transient state |
| CN104598989A (en) * | 2014-12-26 | 2015-05-06 | 神华国华（北京）电力研究院有限公司 | Method for optimizing parameters of elastic multi-mass-block module of shaft system of steam turbine generator unit |
| US20160357893A1 (en) * | 2016-08-15 | 2016-12-08 | Xianwu Ling | Contact stiffness estimation based on structural frequency responses |
| CN106599526A (en) * | 2017-02-17 | 2017-04-26 | 华北电力大学（保定） | Impedance analyzing method for sub-synchronous oscillation of thermal power unit |
| CN108733917A (en) * | 2018-05-17 | 2018-11-02 | 桂林电子科技大学 | A kind of cargo vehicle dynamic modeling method based on transfer matrix |
| CN109145464A (en) * | 2018-08-28 | 2019-01-04 | 暨南大学 | Merge the Structural Damage Identification of multiple target ant lion optimization and the sparse regularization of mark |
-
2019
- 2019-07-19 CN CN201910653770.9A patent/CN110489801B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101515719A (en) * | 2009-04-02 | 2009-08-26 | 华北电力大学 | Method for establishing generator shafting multimass block damping variable model and device thereof |
| CN101807223A (en) * | 2010-03-16 | 2010-08-18 | 浙江大学 | Method for simulating mechanical damping of steam turbine shafting of power system |
| CN103336874A (en) * | 2013-07-15 | 2013-10-02 | 国家电网公司 | Electrical power system subsynchronous oscillation on-line analyzing and early-warning method based on time-domain simulation |
| CN103605819A (en) * | 2013-09-05 | 2014-02-26 | 昆明理工大学 | Simulation method for water turbine generator set shafting vibration transient state |
| CN104598989A (en) * | 2014-12-26 | 2015-05-06 | 神华国华（北京）电力研究院有限公司 | Method for optimizing parameters of elastic multi-mass-block module of shaft system of steam turbine generator unit |
| US20160357893A1 (en) * | 2016-08-15 | 2016-12-08 | Xianwu Ling | Contact stiffness estimation based on structural frequency responses |
| CN106599526A (en) * | 2017-02-17 | 2017-04-26 | 华北电力大学（保定） | Impedance analyzing method for sub-synchronous oscillation of thermal power unit |
| CN108733917A (en) * | 2018-05-17 | 2018-11-02 | 桂林电子科技大学 | A kind of cargo vehicle dynamic modeling method based on transfer matrix |
| CN109145464A (en) * | 2018-08-28 | 2019-01-04 | 暨南大学 | Merge the Structural Damage Identification of multiple target ant lion optimization and the sparse regularization of mark |
Non-Patent Citations (2)
| Title |
|---|
| 和鹏 等: ""电力系统次同步振荡轴系模型研究"", 《电力系统保护与控制》 * |
| 贾锋 等: ""电网故障下风电机组轴系扭振抑制方法"", 《电力自动化设备》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111597680A (en) * | 2020-04-09 | 2020-08-28 | 华北电力大学 | Subsynchronous resonance risk quantitative evaluation method based on electrical damping at modal frequency |
| CN111597680B (en) * | 2020-04-09 | 2024-03-22 | 华北电力大学 | Subsynchronous resonance risk quantitative evaluation method based on electric damping at modal frequency |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110489801B (en) | 2022-10-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lee et al. | Multi‐flexible‐body dynamic analysis of horizontal axis wind turbines | |
| CN101908088B (en) | Time domain bidirectional iteration-based turbine vane flutter stress forecasting method | |
| CN103592081B (en) | A kind of rotor being applicable to distributed amount of unbalance is without test mass dynamic balance method | |
| CN104410107B (en) | A kind of passive integral sliding mode control method of double-fed wind power system | |
| Dandeno et al. | Practical application of eigenvalue techniques in the analysis of power system dynamic stability problems | |
| Voormeeren | Dynamic substructuring methodologies for integrated dydnamic analysis of wind turbines | |
| CN108804853B (en) | Elastic support lower torsional shoulder blade dynamic modeling method based on variable cross-section beam | |
| CN110492479B (en) | Method for identifying rotational inertia and damping of distributed grid-connected equipment | |
| Skjoldan | Aeroelastic modal dynamics of wind turbines including anisotropic effects | |
| CN105117539B (en) | Pneumatic equipment bladess modal frequency and its bimodal Gauss Distribution Fitting method | |
| CN106877363A (en) | A kind of SSSC suppression system sub-synchronous oscillation method and device | |
| CN108923464A (en) | The building of doubly fed induction generator three-phase imbalance steady-state model, emulation mode and device | |
| Ke et al. | Wind field simulation and wind‐induced responses of large wind turbine tower‐blade coupled structure | |
| CN103838141B (en) | Control-orientated large antenna modeling method | |
| CN110489801A (en) | Utilize the generator shafting multiple mass parameter predigesting method of particle swarm optimization algorithm | |
| Bottasso et al. | Wind tunnel testing of wind turbines and farms | |
| CN104408223B (en) | A kind of tracking optimal control method of Wind turbines | |
| Larsen et al. | Generator dynamics in aeroelastic analysis and simulations | |
| CN104795812B (en) | Slow-coherency partitioning method considering stator damping effect and adopting variable substitution | |
| CN103199524B (en) | A kind of tidal current computing method adapting to the access of multiple distributed power source | |
| Larwood | Dynamic analysis tool development for advanced geometry wind turbine blades | |
| CN106026182A (en) | Doubly-fed wind turbine generator set simulation system considering parameter randomness and construction method thereof | |
| Thomas et al. | Event analysis and real‐time validation of doubly fed induction generator‐based wind energy system with grid reactive power exchange under sub‐synchronous and super‐synchronous modes | |
| Rodriguez | Stability and dynamic properties of tip vortices shed from flexible rotors of floating offshore wind turbines | |
| CN104598989A (en) | Method for optimizing parameters of elastic multi-mass-block module of shaft system of steam turbine generator unit |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |