CN105977969B - Extensive stability of muilt-timelag electric power system method of discrimination based on SOD LMS - Google Patents
Extensive stability of muilt-timelag electric power system method of discrimination based on SOD LMS Download PDFInfo
- Publication number
- CN105977969B CN105977969B CN201610408812.9A CN201610408812A CN105977969B CN 105977969 B CN105977969 B CN 105977969B CN 201610408812 A CN201610408812 A CN 201610408812A CN 105977969 B CN105977969 B CN 105977969B
- Authority
- CN
- China
- Prior art keywords
- mrow
- power system
- time
- msub
- matrix
- 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 51
- 239000011159 matrix material Substances 0.000 claims abstract description 62
- 238000001228 spectrum Methods 0.000 claims abstract description 19
- 238000013507 mapping Methods 0.000 claims abstract description 16
- 238000012546 transfer Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000003416 augmentation Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000009897 systematic effect Effects 0.000 claims description 2
- 206010019133 Hangover Diseases 0.000 claims 1
- 238000004891 communication Methods 0.000 abstract description 7
- 230000010355 oscillation Effects 0.000 description 8
- 238000013016 damping Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 6
- 230000005611 electricity Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000001934 delay Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000003381 stabilizer Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000010921 in-depth analysis Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
-
- 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
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Power Engineering (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses the extensive stability of muilt-timelag electric power system method of discrimination based on SOD LMS, time-lag power system model is set up;According to the mapping relations between time-lag power system characteristic value and Solution operator characteristic value, the eigenvalue problem for solving time-lag power system is converted into the spectrum problem for solving Solution operator T (h);Discretization is carried out to Solution operator T (h) using linear multistep method, step-length h is shifted by reasonable selection, obtains accurately differentiating the discretization matrix T of time-lag power system stabilityN;Discretization matrix T is calculated using sequential method or subspace method (such as implicit restarted Arnoldi algorithm)NThe maximum setting number of modulus value characteristic value;According to spectrum mapping relations, by discretization matrix TNCharacteristic value be converted to time-lag power system be located at the complex plane rightmost side characteristic value, the small signal stability for differentiating extensive muilt-timelag electric power system.When SOD LMS methods proposed by the present invention are used for the stability for the critical eigenvalue and judgement system for calculating real system, taken into full account the scale of real system, and communication delay influence.
Description
Technical field
The present invention relates to stability of power system analysis field, and in particular to during a kind of extensive many based on SOD-LMS
Stagnant stability of power system method of discrimination, SOD-LMS is " Solution Operator Discretization-Linear
MultiStep " abbreviation, Chinese implication:Solution operator linear multi step discretization.
Background technology
With the continuous growth of the energy and electricity needs, set up since transregional, transnational large-scale interconnected power system is 21 century
The new stage of electric power development.But at the initial stage of interconnection, complicated electrical structure and weak transmission of electricity link cause interconnected network more
Low-frequency oscillation between easy region.For local oscillation, the coverage of inter-area oscillations is wide, participates in the machine of vibration
Electrical link between group is sufficiently complex, and the influence to the stability of extensive interconnected network is more prominent.Using local signal to be defeated
The power system stabilizer, PSS (Power System Stabilizer, PSS) entered can preferably calm local oscillation pattern, but
It is to be difficult to calm down inter-area oscillation mode, so as to be difficult to ensure that the stability of extensive interconnected network.
The appearance and development of WAMS (Wide-Area Measurement System, WAMS) so that use
The information realization electric power system stability control of global, distal end is possibly realized, to suppress the interregional low of extensive interconnected network
Frequency vibration, which is swung, brings opportunity.WAMS is mainly made up of three parts:Phasor measurement unit (Phasor Measurement Unit,
PMU), communication network, and monitoring system.PMU positioned at plant stand utilizes global positioning system (Global Positioning
System, GPS) high accurate clock signal, the quantity of state of each pivot point of synchro measure power network.Pass through the height in communication system
Target metric data is sent to monitoring system when fast information channel will be posted, and realizes the real-time monitoring, protection and control of power system
System.Metric data in the real-time data base handled by analysis, can not only be used for detecting the low-frequency oscillation letter in power system
Breath, moreover it is possible to which the wide area feedback signal of low frequency oscillation mode between effective reflecting regional is provided for damping controller.And then, optimal control
The performance of device, the ability of the gentle long-distance and large-capacity power transmission of damping water of lifting system.
However, wide area measurement data are in leading to for being made up of different communication medias (such as optical fiber, digital microwave, power line)
When transmitting and handle in letter system, there is the communication delay changed between tens to hundreds of milliseconds.Work as the complicated network structure, and transmission
When data volume is big, Wide-area Measurement Information is actual present in the analyzing and processing of PMU collection, the conveying of communication network and monitoring system
Time lag is often bigger than theoretical value.Time lag is to cause one kind weight of the failure of system control law, operation conditions deterioration and system unstability
Want inducement.When time lag constant is larger, relatively large deviation can occur for system features value, or even can change the small disturbance stability of system completely
Property.In summary, when modern extensive interconnected network carries out power system closed-loop control using wide area measurement information, it is necessary to count and
The influence of time lag.
Consider after time-delay, the state of power system is not only relevant with the state at current time, additionally depends on system mistake
The state gone.Therefore, the model of time-lag power system can be described with a Delay Differential-algebraic equation.With ordinary differential side
System described by journey is different, and the solution space of the power system described by differential equations with delay is infinite dimensional.In a frequency domain, when
The corresponding characteristic equation of stagnant power system, exists and surmounts item, there is infinite multiple solutions (characteristic value).Accordingly, it is considered to after time-delay,
The difficulty of stability of power system research is greatly increased.
At present, in the patent for being related to the differentiation of time-lag power system small signal stability, a kind of time lag electricity of Chinese invention patent
The stable method of discrimination .201010123345.8 [P] of Force system determines the time lag stable region of system using method for tracing is searched for.So
And, time lag characteristic value is solved in the above method and the amount of calculation of search procedure is larger.Chinese invention patent is based on Pad é approximate
Time-lag power system characteristic value is calculated and approached with Convenient stable criterion .201210271783.8 [P] by rational polynominal
Time Delay, and then the critical eigenvalue of the computing system rightmost side, and judge the time lag stability of system.However, the above method
Accuracy and result of calculation precision need through in-depth analysis with research.During extensive based on EIGD of Chinese invention patent
Stagnant power system feature value calculating method .201510055743.3 [P] proposes a kind of discrete based on explicit infinitesimal generator
Change the extensive time-lag power system of (Explicit Infinitesimal Generator Discretization, EIGD)
Feature value calculating method.The above method is needed by Multiple-Scan [0.1,2.5] Hz low-frequency oscillations frequency range, close to the imaginary axis
Critical eigenvalue, the time lag stability of system could be judged.Power System Delay of the Chinese invention patent based on SOD-PS is steady
Qualitative discrimination method .201510229738.X [P] proposes a kind of based on Solution operator puppet spectrum discretization (Solution
Operator Discretization-Pseudospectral, SOD-PS) Power System Delay Convenient stable criterion.Should
Method only needs to calculate a maximum characteristic value of Solution operator discretization matrix norm value or a pair of Con-eigenvalues, it is possible to sentence
Break the stability for system under fixed time lag.However, this method in calculating section critical eigenvalue, it is necessary to ask for discretization
Submatrix it is inverse, it is computationally intensive, calculate the time it is long.
The content of the invention
To solve the deficiency that prior art is present, the invention discloses the extensive multiple time delay power train based on SOD-LMS
System Convenient stable criterion, the small signal stability quick and precisely to differentiate extensive muilt-timelag electric power system.SOD-LMS side
Method is not related to matrix inversion operation during computing system Partial key characteristic value, and calculating speed is fast, computational efficiency is high.Together
When, this method need to only calculate the characteristic value of the maximum setting number of Solution operator discretization matrix norm value, it becomes possible to accurate to differentiate
The stability of extensive muilt-timelag electric power system.
To achieve the above object, concrete scheme of the invention is as follows:
Extensive stability of muilt-timelag electric power system method of discrimination based on SOD-LMS, comprises the following steps:
S1:Set up time-lag power system model;
S2:It is the ODE that Solution operator T (h) is represented by the time-lag power system model conversion obtained in step S1
Initial-value problem;
S3:Discretization is carried out to the Solution operator T (h) in step S2, the discretization matrix T of Solution operator is obtainedN;
S4:Discretization matrix T in calculation procedure S3NThe maximum setting number of modulus value characteristic value μ;
S5:According to spectrum mapping relations, by the discretization matrix T in step S4NCharacteristic value μ be converted to time-lag power system
Eigenvalue λ.
Further, in step one, the time-lag power system model of foundation:In the modeling of actual large-scale electrical power system
During introduce Time Delay, the system model suitable for time-lag power system analysis on Small Disturbance Stability is obtained, during with one group
Stagnant differential equation group is described.
Further, the characteristic equation of time-lag power system is obtained according to one group of time-delayed differential equations, and time lag is electric
The characteristic equation of Force system is expressed as augmentation form of equal value.
Further, in step 2, according to the mapping relations between time-lag power system characteristic value and Solution operator characteristic value,
The eigenvalue problem for solving time-lag power system is converted into the spectrum problem for solving Solution operator T (h).
Further, Solution operator T (h):X → X is defined as the original state at the θ moment in space XWhen being mapped to h+ θ
Etching system state ψ linear operator;Wherein, h is transfer step-length, 0≤h≤τmax;
Wherein:S is integration variable,WithThe state of respectively 0 and h+ θ moment time-lag power systems;0<τ1<τ2
<…<τi…<τmFor the time lag constant of Time Delay, maximum of which time lag is τm=τmax;For systematic observation matrix,
For a dense matrix,It is sparse matrix for system time lags state matrix;Δ x (s) is s moment system state variableses
Increment, Δ x (s- τi) it is s- τiThe increment of moment system state variables,For the increasing of s moment system state variables derivatives
Amount.
Further, in step 3, discretization is carried out to Solution operator T (h) using linear multistep method, obtains accurately sentencing
The discretization matrix T of other time-lag power system stabilityN。
Further, corresponding with Solution operator T (h), discretization matrix TNIt is expressed as follows:
TNLast block row Γ be polynomial eigenvalue problem coefficient matrix, be specifically represented by:
In formula:K is the step number of linear multistep method, αk, βkFor the coefficient of linear multistep method,For
Kronecker product operation,It is sparse matrix for system time lags state matrix,For middle and auxiliary vector.
Further, in the corresponding discretization matrix T of application Solution operatorNSolve the characteristic value of extensive time-lag power system
When, the characteristic value of the maximum setting number of its modulus value is calculated using sequential method or subspace method.
Further, in step 4, it is specially:If k-th Krylov vector representations are qK, then matrix-vector product qK+1
=TNqKFlow it is as follows:
Step (4-1):, will from the direction of rowBoil down to matrix
WhereinJ=1 ..., L+k;L is discrete points;
Step (4-2):qK+1(1:(L+k-1) n, 1)=qK((n+1):(L+k)n,1);
Step (4-3):Using the property of Kronecker product, calculate:
In formula:Vec () is by computing that matrix compression is column vector;
Step (4-4):Calculate qK+1((L+k-1)n:(L+k)n,1):
qK+1((L+k-1)n+1:(L+k) n, 1)=(R)-1z。
Further, discretization matrix TNCharacteristic value μ and time-lag power system eigenvalue λ mapping relations it is as follows:
In formula:H for transfer step-length, σ () represent spectrum, represent exclude.
Beneficial effects of the present invention:
Firstth, SOD-LMS proposed by the present invention is used to calculate the critical eigenvalue of real system and the stabilization of judgement system
During property, taken into full account the scale of real system, and communication delay influence.
Secondth, SOD-LMS proposed by the present invention is not related to Matrix Calculating during computing system Partial key characteristic value
Inverse operation, calculating speed is fast, computational efficiency is high.
3rd, the discretization matrix T that SOD-LMS proposed by the present invention is obtainedNStep-length h is shifted by reasonable selection, can
Ensure the accuracy of stability distinguishing result.
4th, SOD-LMS proposed by the present invention need to only calculate discretization matrix TNThe maximum setting number of modulus value spy
Value indicative, it is possible to differentiate the small signal stability of time-lag power system.
5th, SOD-LMS proposed by the present invention makes full use of the sparse characteristic of sytem matrix, and accurate calculate obtains Solution operator
The critical eigenvalue of modulus value maximum (successively decreasing), the stability of the extensive muilt-timelag electric power system of Quick in discretization matrix.
Brief description of the drawings
Fig. 1 is time-lag power system schematic diagram;
Fig. 2 (a) and Fig. 2 (b) is the graph-based of spectral mapping theorem;
Fig. 3 is the graph-based of transfer step-length h basis for selecting;
Fig. 4 is the flow chart of the extensive stability of muilt-timelag electric power system method of discrimination based on SOD-LMS.
Embodiment:
The present invention is described in detail below in conjunction with the accompanying drawings:
As shown in Figure 1:Time Delay is introduced in the modeling process of actual large-scale electrical power system.Time-lag power system bag
Containing without time-lag power system, wide area Feedback Delays, wide area damping control and wide area output time lag four part, between each several part
Annexation is as shown in the figure.In Fig. 1, yfFor the output without time-lag power system, ydfTo consider the wide area feedback after Feedback Delays
Signal is simultaneously used as the input of damping controller, ycFor the output of wide area damping control, ydcTo consider wide area delayed during output
The output of damping controller, while also serving as the control input of no time-lag power system.
Time lag system is located at the characteristic value of left half complex plane and is mapped in Fig. 2 (b) Solution operator and is located at unit in Fig. 2 (a)
Characteristic value within circle, and time lag system is mapped as Solution operator in Fig. 2 (b) positioned at the characteristic value of right half complex plane in Fig. 2 (a)
Modulus value is more than 1 characteristic value, and outside unit circle.Therefore, the characteristic value of Solution operator is utilized, it is possible to judge former time lag system
The stability of system.If Solution operator at least has the characteristic value that a modulus value is more than 1, it is unstable for may determine that former time lag system
Fixed, if the modulus value of all characteristic values of Solution operator is respectively less than 1, former time lag system is asymptotically stability.
Solid line represents the border that the absolute stability regions of 2 rank backward difference methods are formed after 1/h amplifies, border in Fig. 3
Interior is unstable region.Dotted line represents the characteristic value distributed areas of simple time-lag power system.According to the characteristic of linear multistep method,
Step-length h is shifted by reasonable selection so that the unstable region that solid line is surrounded, which is included, is located at complex plane in characteristic value distributed areas
The part of RHP.When the transfer step-length h of selection meets above-mentioned condition, the discretization matrix obtained by linear multistep method
It is capable of the small signal stability of accurate judgement system.
As shown in Figure 4:Extensive stability of muilt-timelag electric power system method of discrimination based on SOD-LMS, including following step
Suddenly:
S1:Set up time-lag power system model;
S2:The eigenvalue problem for solving time-lag power system is converted into the spectrum problem for solving Solution operator T (h);
S3:By linear multi step discretization, Solution operator T (h) finite dimension discretization matrix T is obtainedN;
S4:The Solution operator discretization matrix T obtained using implicit restarted Arnoldi algorithm come calculation procedure S3NMould
It is worth the characteristic value μ of maximum setting number;
S5:After calculating obtains μ, according to spectrum mapping relations, the eigenvalue λ of time-lag power system is obtained.
So far, the critical eigenvalue of the setting number of judgement system stability has been obtained.
In step S1, it is considered to after time-delay, power system can be described with one group of following time-delayed differential equations:
In formula:For the state variable vector of power system, n is system state variables sum.T is current time.0
<τ1<τ2<…<τi…<τmFor the time lag constant of Time Delay, maximum of which time lag is τm=τmax。For system shape
State matrix, is a dense matrix.It is sparse matrix for system time lags state matrix.Δ x (t) is
The increment of t system state variables, Δ x (t- τi) it is t- τiThe increment of moment system state variables,For t system
The increment of state variable derivative.Δ x (0) is the initial value (i.e. primary condition) of system state variables, and is abbreviated as
The characteristic equation for the time-lag power system that above formula is represented is:
In formula:λ is characterized value, and v is characterized the corresponding right characteristic vector of value.
The augmentation form of equal value of above formula is:
In formula:For middle and auxiliary vector.If InFor n rank unit matrixs, then A'(λ) and B'(λ) can specific table
It is shown as:
In formula:A0, B0, C0, D0, Ai, BiFor highly sparse linearisation coefficient matrix.
In step S2, Solution operator T (h):X → X is defined as the primary condition (state) at the θ moment in space XIt is mapped to
H+ θ moment system modes ψ linear operator;Wherein, h is transfer step-length, 0≤h≤τmax;
Wherein:S is integration variable,WithThe state of respectively 0 and h+ θ moment time-lag power systems.Time lag electricity
Relation between characteristic value and Solution operator T (h) spectrum of Force system model:
There is such as ShiShimonoseki between spectral mapping theorem, Solution operator T (h) characteristic value μ and the eigenvalue λ of time lag system
System:
In formula:H for transfer step-length, σ () represent spectrum, represent exclude.
T (h) is by the stable region of time-lag power system, i.e., left half s Planar Mappings are within the unit circle of z-plane.This causes
Several characteristic values that can be maximum only by calculating T (h) modulus value, you can judge the stability of time-lag power system.If | μ |>1,
Then corresponding λ real part is more than zero, i.e. Re (λ)>0, system is unstable.If | μ |<1, then Re (λ)<0, system asymptotically stability.If |
μ |=1, then Re (λ)=0, system neutrality.
Solution operator T (h) is the Infinite Dimension Linear operator for describing X → X mappings., can from the spectrum map feature of Solution operator
With the spectrum by calculating Solution operator, the Partial key characteristic value of time-lag power system is obtained.However, calculating T (h) characteristic value is
One Infinite-dimensional eigenvalue problem.Therefore, it is necessary to first to T (h) carry out discretization, then by calculate corresponding finite dimension from
The characteristic value of dispersion matrix, using the partial feature value as time lag system.
The step of step S3, is as follows:
, discretization matrix T corresponding with Solution operator T (h)NIt is expressed as follows:
TNLast block row Γ be polynomial eigenvalue problem coefficient matrix, be specifically represented by:
In formula:K is the step number of linear multistep method,For Kronecker product operation, especially, when
When system is integer only containing a time lag and N=τ/h, Γ can be explicitly expressed as follows:
Γ=[Γ0,0n×n(N-k-1),Γ1]
In formula:αj,βj(j=0 ..., it is k) coefficient of linear k footworks
In step S4, discretization matrix TNExponent number be (L+k) n.For large-scale electrical power system, matrix TNExponent number will
It is very huge.Therefore, in the corresponding discretization matrix T of application Solution operatorNWhen solving the time lag characteristic value of large-scale electrical power system,
The characteristic value of the maximum setting number of its modulus value must be calculated using sparse features value-based algorithm.
The step of step S4, is as follows:
If k-th Krylov vector representations are qK, then matrix-vector product qK+1=TNqKFlow it is as follows:
Step (4-1):, will from the direction of rowBoil down to matrix
WhereinJ=1 ..., L+k;L is discrete points.
Step (4-2):qK+1(1:(L+k-1) n, 1)=qK((n+1):(L+k)n,1);
Step (4-3):Using the property of Kronecker product, calculate:
In formula:Vec () is by computing that matrix compression is column vector.
Step (4-4):Calculate qK+1((L+k-1)n:(L+k)n,1):
qK+1((L+k-1)n+1:(L+k) n, 1)=(R)-1z
In step (4-3) it can be seen from above-mentioned flowIn step (44) (R)-1Z calculating
Amount is maximum, time-consuming most long.So as to deduce, the amount of calculation of whole flow process is approximately equal to L+k+1 general characteristics value point of progress
The amount of calculation of analysis.It is worth noting that, computationally intensive two that time-consuming can be with sparse reality in step (4-3) and step (4-4)
It is existing, so as to reduce amount of calculation, improve computational efficiency.
In step S5, after calculating obtains μ, time-lag power system is obtained positioned at the complex plane rightmost side by spectrum mapping
The eigenvalue λ of number is set, so as to judge the stability of system.Discretization matrix TNCharacteristic value μ and time-lag power system it is special
Value indicative λ mapping relations are as follows:
In formula:H for transfer step-length, σ () represent spectrum, represent exclude.
The present invention sets up the model of the electric power system model and wide area damping control that do not consider time-delay respectively, passes through
Time Delay is introduced, the mathematical modeling for the closed loop power system for considering time-delay is set up.By time-lag power system model at it
Linearized near steady-state operation point, obtain the system model available for time-lag power system analysis on Small Disturbance Stability.
By means of Solution operator T (h), by primary condition (state)It is mapped as Account Dept and decomposes (state) ψ, obtains sign and reflect
Penetrate the expression formula of relation.In order to obtain the spectrum of Solution operator, linear multi step is carried out to the Solution operator T (h) of the time-lag power system
Discretization, obtains the finite dimension Solution operator discretization matrix T with explicit expressionN。
Using implicit restarted Arnoldi algorithm, discretization matrix T is calculatedNThe characteristic value of the maximum setting number of modulus value
μ.After calculating obtains μ, the eigenvalue λ that time-lag power system is located at the complex plane rightmost side is obtained according to spectrum mapping relations, is used
In the stability of judgement system.
The SOD-LMS methods of the present invention need to only be calculated by a sub-eigenvalue, it is possible to quick to obtain being used for judgement system
The critical eigenvalue positioned at the complex plane rightmost side of stability.
Although above-mentioned the embodiment of the present invention is described with reference to accompanying drawing, not to present invention protection model
The limitation enclosed, one of ordinary skill in the art should be understood that on the basis of technical scheme those skilled in the art are not
Need to pay various modifications or deform still within protection scope of the present invention that creative work can make.
Claims (8)
1. the extensive stability of muilt-timelag electric power system method of discrimination based on SOD-LMS, it is characterized in that, comprise the following steps:
S1:Set up time-lag power system model;
S2:By the initial value that the time-lag power system model conversion obtained in step S1 is the ODE that Solution operator T (h) is represented
Problem;
S3:Discretization is carried out to the Solution operator T (h) in step S2, the discretization matrix T of Solution operator is obtainedN;
S4:Discretization matrix T in calculation procedure S3NThe maximum setting number of modulus value characteristic value μ;
S5:According to spectrum mapping relations, by the discretization matrix T in step S4NCharacteristic value μ be converted to the spy of time-lag power system
Value indicative λ;
In step 3, discretization is carried out to Solution operator T (h) using linear multistep method, obtains accurately differentiating time lag power train
The discretization matrix T for stability of unitingN;
, discretization matrix T corresponding with Solution operator T (h)NIt is expressed as follows:
TNLast block row Γ be polynomial eigenvalue problem coefficient matrix, be specifically represented by:
<mrow>
<mi>&Gamma;</mi>
<mo>=</mo>
<msup>
<mi>R</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
In formula:K is the step number of linear multistep method,For Kronecker product operation,For system
Hangover state matrix, is sparse matrix,For middle and auxiliary vector.
2. the extensive stability of muilt-timelag electric power system method of discrimination as claimed in claim 1 based on SOD-LMS, its feature
It is, in step one, the time-lag power system model of foundation:Time lag is introduced in the modeling process of actual large-scale electrical power system
Link, obtains the system model suitable for time-lag power system analysis on Small Disturbance Stability, is come with one group of time-delayed differential equations
Description.
3. the extensive stability of muilt-timelag electric power system method of discrimination as claimed in claim 2 based on SOD-LMS, its feature
It is that the characteristic equation of time-lag power system is obtained according to one group of time-delayed differential equations, and by the feature side of time-lag power system
Journey is expressed as augmentation form of equal value.
4. the extensive stability of muilt-timelag electric power system method of discrimination as claimed in claim 1 based on SOD-LMS, its feature
It is that in step 2, according to the mapping relations between time-lag power system characteristic value and Solution operator characteristic value, time lag electric power will be solved
The eigenvalue problem of system is converted into the spectrum problem for solving Solution operator T (h).
5. the extensive stability of muilt-timelag electric power system method of discrimination based on SOD-LMS as described in claim 1 or 4, its
It is characterized in, Solution operator T (h):X → X is defined as the original state at the θ moment in space XIt is mapped to h+ θ moment system modes ψ
Linear operator;Wherein, h is transfer step-length, 0≤h≤τmax;
Wherein:S is integration variable,WithThe state of respectively 0 and h+ θ moment time-lag power systems;0<τ1<τ2<…<
τi…<τmFor the time lag constant of Time Delay, maximum of which time lag is τm=τmax;It is one for systematic observation matrix
Individual dense matrix,It is sparse matrix for system time lags state matrix;Δ x (s) is the increasing of s moment system state variableses
Amount, Δ x (s- τi) it is s- τiThe increment of moment system state variables,For the increment of s moment system state variables derivatives.
6. the extensive stability of muilt-timelag electric power system method of discrimination as claimed in claim 1 based on SOD-LMS, its feature
It is, in the corresponding discretization matrix T of application Solution operatorNWhen solving the characteristic value of extensive time-lag power system, using sequential method
Or subspace method calculates the characteristic value of the maximum setting number of its modulus value.
7. the extensive stability of muilt-timelag electric power system method of discrimination based on SOD-LMS as described in claim 1 or 6, its
It is characterized in step 4, to be specially:If k-th Krylov vector representations are qK, then matrix-vector product qK+1=TNqKFlow
It is as follows:
Step (4-1):, will from the direction of rowBoil down to matrixWhereinJ=1 ..., L+k;L is discrete points;
Step (4-2):qK+1(1:(L+k-1) n, 1)=qK((n+1):(L+k)n,1);
Step (4-3):Using the property of Kronecker product, calculate
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>z</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>q</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>v</mi>
<mi>e</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>v</mi>
<mi>e</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<msubsup>
<mi>Pw</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>v</mi>
<mi>e</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>L</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<msubsup>
<mi>w</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
In formula:Vec () is by computing that matrix compression is column vector;
Step (4-4):Calculate qK+1((L+k-1)n:(L+k)n,1):
qK+1((L+k-1)n+1:(L+k) n, 1)=(R)-1z。
8. the extensive stability of muilt-timelag electric power system method of discrimination as claimed in claim 1 based on SOD-LMS, its feature
It is, discretization matrix TNCharacteristic value μ and time-lag power system eigenvalue λ mapping relations it is as follows:
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>h</mi>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mi>&mu;</mi>
<mo>,</mo>
<mi>&mu;</mi>
<mo>&Element;</mo>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>\</mo>
<mo>{</mo>
<mn>0</mn>
<mo>}</mo>
</mrow>
In formula:H for transfer step-length, σ () represent spectrum, represent exclude.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610408812.9A CN105977969B (en) | 2016-06-08 | 2016-06-08 | Extensive stability of muilt-timelag electric power system method of discrimination based on SOD LMS |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610408812.9A CN105977969B (en) | 2016-06-08 | 2016-06-08 | Extensive stability of muilt-timelag electric power system method of discrimination based on SOD LMS |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105977969A CN105977969A (en) | 2016-09-28 |
CN105977969B true CN105977969B (en) | 2017-10-03 |
Family
ID=57011087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610408812.9A Active CN105977969B (en) | 2016-06-08 | 2016-06-08 | Extensive stability of muilt-timelag electric power system method of discrimination based on SOD LMS |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105977969B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109685400B (en) * | 2018-02-24 | 2020-07-31 | 山东大学 | Time-lag power system stability discrimination method based on time integral IGD |
CN108923421B (en) * | 2018-07-13 | 2020-04-21 | 山东大学 | Time-lag power system key characteristic value calculation method based on low-order SOD-IRK algorithm |
CN108808705B (en) * | 2018-07-13 | 2020-05-22 | 山东大学 | Time-lag power system electromechanical oscillation mode calculation method based on low-order SOD-PS-II-R algorithm |
CN109033022B (en) * | 2018-07-13 | 2020-07-31 | 山东大学 | Time-lag power system characteristic value calculation method based on low-order SOD-L MS algorithm |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217896A (en) * | 2013-03-29 | 2013-07-24 | 国家电网公司 | Multi-FACTS (flexible AC transmission system) anti-delay coordination and control method based on free-weighting matrix method |
CN103227467A (en) * | 2013-04-19 | 2013-07-31 | 天津大学 | Lyapunov stability analysis method of time delay electric system |
CN103414245A (en) * | 2013-06-04 | 2013-11-27 | 浙江工业大学 | Quantization-based wide-area power system output feedback control method |
CN105468909A (en) * | 2015-11-24 | 2016-04-06 | 山东大学 | Time delay power system electromechanical oscillation mode computing method based on SOD-PS-R R algorithm |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9043040B2 (en) * | 2012-01-31 | 2015-05-26 | General Electric Company | Power system stabilization |
US9876352B2 (en) * | 2013-02-26 | 2018-01-23 | Washington State University | Voltage stability monitoring in power systems |
-
2016
- 2016-06-08 CN CN201610408812.9A patent/CN105977969B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217896A (en) * | 2013-03-29 | 2013-07-24 | 国家电网公司 | Multi-FACTS (flexible AC transmission system) anti-delay coordination and control method based on free-weighting matrix method |
CN103227467A (en) * | 2013-04-19 | 2013-07-31 | 天津大学 | Lyapunov stability analysis method of time delay electric system |
CN103414245A (en) * | 2013-06-04 | 2013-11-27 | 浙江工业大学 | Quantization-based wide-area power system output feedback control method |
CN105468909A (en) * | 2015-11-24 | 2016-04-06 | 山东大学 | Time delay power system electromechanical oscillation mode computing method based on SOD-PS-R R algorithm |
Non-Patent Citations (1)
Title |
---|
时变时滞电力***鲁棒稳定性的改进型判据;孙国强;《电力***自动化》;20150210;第39卷(第3期);第59-62页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105977969A (en) | 2016-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105977969B (en) | Extensive stability of muilt-timelag electric power system method of discrimination based on SOD LMS | |
Snyder et al. | Delayed-input wide-area stability control with synchronized phasor measurements and linear matrix inequalities | |
CN105449665B (en) | Time-lag power system stability method of discrimination based on SOD PS | |
CN105932666A (en) | Complex distribution network multi-time-scale digital-analog hybrid simulation system and simulation method therefor | |
CN101505061B (en) | Computation method capable of tracking Davinan equivalence parameter base on time domain simulation | |
CN105184027B (en) | A kind of power load modelling approach based on interacting multiple model algorithm | |
CN105574809B (en) | Electromagnetic transient simulation graphics processor parallel calculating method based on matrix exponetial | |
CN106250640A (en) | A kind of layering Dynamic Equivalence being applicable to area power grid | |
CN102801158A (en) | Method for calculating time-lag electric power system eigenvalue and discriminating stability based on Pade approximation | |
CN109274116A (en) | A kind of alternating current-direct current mixing micro-capacitance sensor continuous tide calculation method | |
CN105468909A (en) | Time delay power system electromechanical oscillation mode computing method based on SOD-PS-R R algorithm | |
CN108462192B (en) | Mode analysis method for broadband oscillation of power system | |
CN103258103A (en) | Thevenin equivalent parameter identification method based on partial least squares regression | |
CN104199302A (en) | Molding system and method of pump storage group speed regulating system | |
CN103984822A (en) | Three phase flow implementation method based on advanced digital power system simulator | |
CN106655190A (en) | Method for solving P-OPF (Probabilistic-Optimal Power Flow) of wind power stations | |
CN104917174A (en) | Static voltage stability judgment method for active power distribution network | |
CN109787250A (en) | A kind of low-frequency oscillation of electric power system pattern identification method based on polynary empirical mode decomposition | |
CN106300345A (en) | Based on the low-frequency oscillation parameter identification method improving Prony algorithm | |
CN103760400A (en) | Voltage stabilization on-line monitoring method based on wide area measuring system | |
CN107834593A (en) | A kind of droop control type isolated island micro-capacitance sensor static voltage stability probability evaluation method of failure | |
CN110460085A (en) | A method of consider wind-powered electricity generation and part throttle characteristics to effect on power system | |
CN104950261B (en) | The hardware-in-the-loop test method and system of battery | |
CN105244887A (en) | Closed-loop control method for transient power angle instability of power system | |
CN109149566A (en) | A kind of modeling method of the simulation model of the high-power minimum point prediction of missing lower frequency |
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 |