CN102799765B - The real-time unfixed point interpolation computing method of switch in a kind of step-length - Google Patents
The real-time unfixed point interpolation computing method of switch in a kind of step-length Download PDFInfo
- Publication number
- CN102799765B CN102799765B CN201210210272.5A CN201210210272A CN102799765B CN 102799765 B CN102799765 B CN 102799765B CN 201210210272 A CN201210210272 A CN 201210210272A CN 102799765 B CN102799765 B CN 102799765B
- Authority
- CN
- China
- Prior art keywords
- point
- state value
- length
- calculating
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 36
- 230000033001 locomotion Effects 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 36
- 230000001052 transient effect Effects 0.000 claims abstract description 19
- 230000010354 integration Effects 0.000 claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 15
- 230000004048 modification Effects 0.000 claims abstract description 12
- 238000012986 modification Methods 0.000 claims abstract description 12
- 238000013016 damping Methods 0.000 claims abstract description 11
- 238000013213 extrapolation Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 9
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- BYACHAOCSIPLCM-UHFFFAOYSA-N 2-[2-[bis(2-hydroxyethyl)amino]ethyl-(2-hydroxyethyl)amino]ethanol Chemical group OCCN(CCO)CCN(CCO)CCO BYACHAOCSIPLCM-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Abstract
The present invention relates to switch real-time unfixed point interpolation computing method in a kind of step-length.The method comprises the steps：First step-length point-like state value of electro-magnetic transient is calculated, finds there is switch motion in step-length；Interpolated value calculating is carried out to switch motion point state value；Next beans-and bullets shooter long status value is calculated from switch motion point；From current step point extrapolated value to next calculating synchronizing point long.When the method carries out switching characteristic circuit counting, it is possible to increase the calculating speed of switch subnet, real-time is conducive to calculate.So as to provide time margin for the precision that switch is calculated, can be applied to real-time, super real-time, off-line calculation switching characteristic circuit electromagnetic transient simulation calculate in, wherein whole step size computation can be mutually compatible with various electro-magnetic transient switching algorithms and be not limited to such as trapezoidal integration, backward-Euler method, the trapezoidal integration with damping or its modification change combination, the calculating speed and precision of emulation can be effectively improved, optimization electro-magnetic transient calculates performance in real time.
Description
Technical field
The present invention relates to a kind of computational methods, and in particular to the real-time unfixed point interpolation computing method of switch in a kind of step-length.
Background technology
Time-domain-simulation has become the important tool of Power System Analysis, design and research, with Power Electronic Technique
Extensive use, flexible AC transmission, D.C. high voltage transmission and distributed power generation in interior power system links, especially
It is during a large amount of renewable energy power generation equipment are typically necessary and could access power network by electronic power convertor, and power electronics
This network topology structure being substantially continually changing with the time of switch proposes new to traditional Power System Time Domain Simulation Under
It is required that and challenge, especially electro-magnetic transient real-time simulation aspect, on the premise of how ensuring to calculate real-time, calculate frequent
The electronic power switch action put for a long time with non-synchronizing, is a thorny technological difficulties.
The content of the invention
In view of the shortcomings of the prior art, the present invention provides the real-time unfixed point interpolation computing method of switch in a kind of step-length, should
Method on the basis of existing electromagnetic transient in power system software etc., when carrying out switching characteristic circuit counting, it is possible to increase
The calculating speed of subnet is switched, is conducive to real-time to calculate.So as to provide time margin for the precision that switch is calculated, can answer
It is temporary for real-time, super real-time, off-line calculation the electromagnetism comprising switching characteristic circuits such as electronic power switch, high frequency electrical switches
In state simulation calculation, wherein whole step size computation can it is mutually compatible with various electro-magnetic transient switching algorithms and be not limited to as：Trapezoidal product
Point-score, backward-Euler method, the trapezoidal integration with damping or its modification change combination etc., can effectively improve the calculating speed of emulation
And precision, optimize electro-magnetic transient and calculate performance in real time.
The purpose of the present invention is realized using following technical proposals：
The real-time unfixed point interpolation computing method of switch in a kind of step-length, it is theed improvement is that, methods described includes following
Step：
A, first step-length point-like state value of electro-magnetic transient is calculated, find there is switch motion in step-length；
B, interpolated value calculating is carried out to switch motion point state value；
C, next beans-and bullets shooter long status value is calculated from the switch motion point；
D, from current step point extrapolated value to it is next calculating synchronizing point long.
Preferably, in the step B, the state value x based on previous step1With current step point-like state value x2, to the switch
The state value x at operating point k1 momentk1Carry out interpolated value calculating, the xk1Represented with following 1. formulas：
①；
In formula：xk1It is the state value at switch motion point k1 moment；tk1It is k1 time points at moment；t1It is state value x1Place
The whole material calculation moment；Δ t is a time span for the whole material calculation of standard.
Preferably, in the step C, from the switch motion point xk1Using electromagnetic switch algorithm to next beans-and bullets shooter shape long
State value xk1sCalculated.
More preferably, it is characterised in that the electromagnetic switch algorithm is using trapezoidal integration, backward-Euler method, band damping
Trapezoidal integration or its modification change combination.
Preferably, in the step D, from the state value x of step C current step-length pointk1sThe next calculating synchronizing of extrapolated value is long
Point-like state value x3, the x3Represented with following 2. formulas：
②；
In formula：x3It is next calculating synchronizing point-like state value long；t3It is state value x3The whole material calculation moment at place.
Preferably, methods described each step calculating after, either using backward-Euler method, hiding-trapezium integral method,
Trapezoidal integration with damping or its modification change combination, will carry out event searching and treatment.Event searching is to calculate one
After individual step-length it is necessary to detect whether have in this step-length calculate in do not considered switch motion, network structure change etc., if
Have it is necessary to be processed, so-called treatment is exactly to be calculated using method of the present invention interpolated value or extrapolated value.
The amount to be asked of interpolation algorithm is the state value at certain the calculating moment inside the current step-length cycle for calculating, and is adopted
Linear interpolation calculating is carried out to certain state value in step-length span with known state value；The amount to be asked of extrapolation value-based algorithm be
The state value at certain the calculating moment outside the current step-length cycle for calculating, using known state value to the shape outside step-length span
State value carries out linear interpolation calculating.
Compared with the prior art, the beneficial effect that reaches of the present invention is：
1st, the present invention proposes that the present invention provides the real-time unfixed point interpolation computing method of switch in a kind of step-length, to solve one
The power system circuit of generation multiple switch action can not meet calculating and want in electromagnetic transient simulation medium velocity and precision in step-length
Ask and propose the present invention, the present invention calculates method for designing using new switch, Multi- Switch occurs in one step-length of reply calculates
Situation, amount of calculation and the calculating time of the EM transient calculation containing switching characteristic circuit can be greatly reduced, can be applied in real time,
Super real-time, off-line calculation the electromagnetic transient simulation meter comprising switching characteristic circuits such as electronic power switch, high frequency electrical switches
In calculation, can be mutually compatible with various electro-magnetic transient switching algorithms, the calculating speed and precision of emulation are effectively improved, optimization electromagnetism is temporary
State calculates performance in real time.
2nd, the present invention on the basis of existing electromagnetic transient in power system software etc., carry out switching characteristic circuit counting when
Wait, it is possible to increase the calculating speed of switch subnet, be conducive to real-time to calculate.During so as to being provided for the precision that switch is calculated
Between nargin, can be applied to it is real-time, super it is real-time, off-line calculation comprising electronic power switch, high frequency electrical switch etc. switching characteristic
During the electromagnetic transient simulation of circuit is calculated, wherein whole step size computation can be mutually compatible with various electro-magnetic transient switching algorithms and limit
Yu Ru：Trapezoidal integration, backward-Euler method, the trapezoidal integration with damping or its modification change combination etc., can effectively improve imitative
Genuine calculating speed and precision, optimization electro-magnetic transient calculate performance in real time.
3rd, the present invention, can be quick using the algorithm of interpolated value in the case of the single or multiple switch motions of step-length
Calculate each switch quantity of state, due to interpolated value calculate speed quickly, it is possible to increase calculating speed, be conducive in real time calculating
Realization.
4th, the present invention reduces the number of times of the whole step size computation of multistep to the full extent, and no matter switch motion is more in a step-length
It is few secondary, on off state point interpolation calculation, method of the present invention design is for improving the speed of simulation calculation and imitating in real time
The true time reliability for calculating has very important meaning.
Brief description of the drawings
Fig. 1 is the calculation procedure figure of the interior real-time unfixed point interpolation computing method of switch of step-length that the present invention is provided, wherein, 1.
For fixed point fixed step size is calculated；2. for interpolated value is calculated；3. for unfixed point fixed step size is calculated；It is switch motion moment point；☆ is
Setpoint step size Connection Time point；6. for extrapolated value is calculated；
Fig. 2 is the workflow diagram of the interior real-time unfixed point interpolation computing method of switch of step-length that the present invention is provided；
Fig. 3 is that have the calculating in the case of more switch motions in a step-length of the specific embodiment that the present invention is provided
Block diagram, wherein：1. for fixed point fixed step size is calculated；2. and 4. for interpolated value is calculated；3. and 5. for unfixed point fixed step size is calculated；
It is switch motion moment point；☆ is setpoint step size Connection Time point；6. for extrapolated value is calculated.
Specific embodiment
Specific embodiment of the invention is described in further detail below in conjunction with the accompanying drawings.
The calculation procedure figure and workflow diagram of the real-time unfixed point interpolation computing method of switch in the step-length that the present invention is provided
Respectively as illustrated in fig. 1 and 2, the method comprises the steps：
A, first step-length point-like state value of electro-magnetic transient is calculated, find there is switch motion in step-length；
B, interpolated value calculating is carried out to switch motion point state value：State value x based on previous step1With current step point-like
State value x2, to the state value x at the switch motion point k1 momentk1Carry out interpolated value calculating, the xk1Represented with following 1. formulas：
①；
In formula：xk1It is the state value at switch motion point k1 moment；tk1It is k1 time points at moment；t1It is state value x1Place
The whole material calculation moment；Δ t is a time span for the whole material calculation of standard.
C, next beans-and bullets shooter long status value is calculated from the switch motion point：From switch motion point xk1Using electromagnetism
Switching algorithm is to next beans-and bullets shooter long status value xk1sCalculated.
Electromagnetic switch algorithm is using trapezoidal integration, backward-Euler method, the trapezoidal integration with damping or its modification change
Combination.
D, from current step point extrapolated value to it is next calculating synchronizing point long：From the state value x of step C current step-length pointk1s
Extrapolated value is next to calculate synchronizing point-like state value x long3, the x3Represented with following 2. formulas：
②；
In formula：x3It is next calculating synchronizing point-like state value long；t3It is state value x3The whole material calculation moment at place.
After the calculating of each step, either using backward-Euler method, hiding-trapezium integral method, the trapezoidal integration with damping
Method or its modification change combination, will carry out event searching and treatment.Event searching is it is necessary to detect after a step-length is calculated
Whether there is the switch motion not considered in calculating, network structure to change in this step-length etc., if it has, will be processed, institute
Meaning treatment is exactly to be calculated using method of the present invention interpolated value or extrapolated value.
With reference to specific embodiment, the present invention is described in further detail.
Embodiment
There is the calculation procedure in the case of more switch motions in one step-length of the specific embodiment that the present invention is provided
Figure is as shown in figure 3, comprise the following steps：
（1）First calculating of step-length is carried out, finds there is switch motion in step-length；
（2）Interpolated value calculating is carried out to first state value of switch motion point：State value x based on previous step-length point1
With the state value x of this step-length point2, to first state value x at switch motion point k1 momentk1Carry out interpolated value calculating, xk1With
Following 1. formulas are represented：
①；
In formula：xk1It is first state value at switch motion point k1 moment；tk1It is k1 time points at moment；t1It is state value x1
The whole material calculation moment at place；Δ t is a time span for the whole material calculation of standard.
（3）From first switch motion point xk1Carry out next step-length point xk1sCalculating；
（4）It was found that second switch motion point, based on first state value x of switch motion pointk1With step size computation point x2
State value, to second switch calculate point state value xk2Carry out interpolated value calculating, xk2Represented with following 3. formulas：
③；
In formula：xk2It is to second switch motion point k2 moment state value；tk2It is k2 time points at moment；t2It is state value x2
The whole material calculation moment at place.
（5）From second state value x of switch motion pointk2Using band damping hiding-trapezium integral method, backward-Euler method, band
The trapezoidal integration of damping or its modification change combination carry out next step-length point-like state value xk2sCalculating；
（6）From previous step current step point xk2sState value x of the extrapolated value to next calculating synchronizing point long3, x3With it is following 2.
Formula is represented：
②；
In formula：x3It is next calculating synchronizing point-like state value long；t3It is state value x3The whole material calculation moment at place.
After the calculating of each step, either using backward-Euler method, hiding-trapezium integral method, the trapezoidal integration with damping
Method or its modification change combination, will carry out event searching and treatment.Event searching is it is necessary to detect after a step-length is calculated
Whether there is the switch motion not considered in calculating, network structure to change in this step-length etc., if it has, will be processed, institute
Meaning treatment is exactly to be calculated using method of the present invention interpolated value or extrapolated value.
Finally it should be noted that：The above embodiments are merely illustrative of the technical solutions of the present invention rather than its limitations, to the greatest extent
Pipe has been described in detail with reference to above-described embodiment to the present invention, and those of ordinary skill in the art should be understood：Still
Specific embodiment of the invention can be modified or equivalent, and without departing from any of spirit and scope of the invention
Modification or equivalent, it all should cover in the middle of scope of the presently claimed invention.
Claims (1)
1. real-time unfixed point interpolation computing method is switched in a kind of step-length, it is characterised in that methods described comprises the steps：
A, first step-length point-like state value of electro-magnetic transient is calculated, find there is switch motion in step-length；
B, interpolated value calculating is carried out to switch motion point state value；
C, next step-length point-like state value is calculated from the switch motion point；
D, from current step point extrapolated value to it is next calculating synchronizing point long；
In the step B, the state value x based on previous step1With current step point-like state value x2, during to the switch motion point k1
The state value x at quarterk1Carry out interpolated value calculating, the xk1Represented with following 1. formulas：
In formula：xk1It is the state value at switch motion point k1 moment；tk1It is k1 time points at moment；t1It is state value x1The whole meter at place
Calculate the step-length moment；Δ t is a time span for the whole material calculation of standard；
In the step C, from the state value x at the switch motion point k1 momentk1Using electromagnetic switch algorithm to next step-length
Point-like state value xk1sCalculated；
The electromagnetic switch algorithm is using trapezoidal integration, backward-Euler method, the trapezoidal integration with damping or its modification change
Combination；
In the step D, the next step-length point-like state value x described in step Ck1sExtrapolated value is next to calculate synchronizing dotted state long
Value x3, the x3Represented with following 2. formulas：
In formula：x3It is next calculating synchronizing point-like state value long；t3It is state value x3The whole material calculation moment at place；
Methods described each step calculating after, either using backward-Euler method, hiding-trapezium integral method, band damp it is trapezoidal
Integration method or its modification change combination, will carry out event searching and treatment；Event searching be calculate a step-length after it is necessary to
Detect whether there is the switch motion not considered in calculating, network structure to change in this step-length, if it has, will be processed,
The treatment is exactly to be calculated using the interpolated value or extrapolated value of step A~D；
The amount to be asked of interpolation algorithm is the state value at certain the calculating moment inside the current step-length cycle for calculating, using
The state value known carries out linear interpolation calculating to certain state value in step-length span；The amount to be asked of extrapolation value-based algorithm is current
The state value at certain the calculating moment outside the step-length cycle of calculating, using known state value to the state value outside step-length span
Carry out linear interpolation calculating.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210210272.5A CN102799765B (en) | 2012-06-19 | 2012-06-19 | The real-time unfixed point interpolation computing method of switch in a kind of step-length |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210210272.5A CN102799765B (en) | 2012-06-19 | 2012-06-19 | The real-time unfixed point interpolation computing method of switch in a kind of step-length |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102799765A CN102799765A (en) | 2012-11-28 |
| CN102799765B true CN102799765B (en) | 2017-06-23 |
Family
ID=47198874
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210210272.5A Active CN102799765B (en) | 2012-06-19 | 2012-06-19 | The real-time unfixed point interpolation computing method of switch in a kind of step-length |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102799765B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104461672B (en) * | 2013-09-12 | 2018-03-13 | 南京南瑞继保电气有限公司 | A kind of method and device for extending power system digital simulation model emulation total duration |
| CN106611069B (en) * | 2015-10-26 | 2021-10-22 | 范圣韬 | Variable-step-size element adjoint model generation method for constant impedance |
| CN112711720A (en) * | 2020-12-31 | 2021-04-27 | 广州宸祺出行科技有限公司 | Boarding point recommendation method and system for reducing travel cost |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101382969A (en) * | 2008-10-31 | 2009-03-11 | 中国电力科学研究院 | Multi-step step length variable electromagnetical transient emulation method |
| CN101692608A (en) * | 2009-08-03 | 2010-04-07 | 和芯微电子(四川)有限公司 | High linearity variable bit width interpolator |
-
2012
- 2012-06-19 CN CN201210210272.5A patent/CN102799765B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101382969A (en) * | 2008-10-31 | 2009-03-11 | 中国电力科学研究院 | Multi-step step length variable electromagnetical transient emulation method |
| CN101692608A (en) * | 2009-08-03 | 2010-04-07 | 和芯微电子(四川)有限公司 | High linearity variable bit width interpolator |
Non-Patent Citations (1)
| Title |
|---|
| 一种计及多重开关的电力电子时域仿真插值算法;王成山等;《电工技术学报》;20100630;正文第1－2节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102799765A (en) | 2012-11-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104866665B (en) | Based on the equivalent hybrid simulation method containing power electronic equipment with interacting of interface | |
| CN102750416B (en) | Topological subnetting method of electromagnetic transient simulation containing switching characteristic circuit | |
| CN104462661B (en) | Electronic power switch interpolation real-time emulation method | |
| CN109271698A (en) | A kind of mode of resonance double active bridging parallel operation modelings, depression of order, design method, apparatus and system | |
| CN105610158B (en) | A kind of Distributed Power Flow controller and its control method | |
| CN104158405A (en) | Two-way DC/DC converter for micro-grid and magnetic bias digital suppression method | |
| CN105896587A (en) | Multi-port UPFC topology and appropriate configuration and control method thereof | |
| CN106096106B (en) | High-frequency high-voltage transformer for electrostatic dust collection optimum design method | |
| CN104809308B (en) | A kind of transverter switch function modeling method suitable for asymmetric operating state | |
| CN102799765B (en) | The real-time unfixed point interpolation computing method of switch in a kind of step-length | |
| CN103020385A (en) | RTDS-based 500-Kv major network modeling simulation systems for power network | |
| CN107612409A (en) | Simplify matrix converter model predictive control method with magnetic saturation control | |
| CN102916450A (en) | Hybrid real-time simulation method for three-phase asymmetrical alternating-current/direct-current power system | |
| CN106549377B (en) | A kind of handover method of alternating current-direct current combined hybrid system electromechanics-electromagnetic transient hybrid simulation | |
| CN102163845B (en) | Optimal configuration method of distributed generations (DG) based on power moment algorithm | |
| CN103729502B (en) | A kind of method improving electromagnetic transient in power system simulation velocity | |
| CN106229995B (en) | Based on the backup power source shunt reactor parameter optimization method under the Anti-Typhoon operational mode of wind power plant | |
| Yu et al. | Quadratic models of AC–DC power flow and optimal reactive power flow with HVDC and UPFC controls | |
| Miao et al. | Generalized steady-state model for energy router with applications in power flow calculation | |
| CN105119269A (en) | Random power flow calculation method taking regard of multiterminal unified power flow controller | |
| CN102570952B (en) | Demagnetizing control method for double-feeding type wind power generation system to LVRT (Low Voltage Ride Through) | |
| CN110445174A (en) | A kind of windfarm system transient stability appraisal procedure considering random air speed influence | |
| CN103378595B (en) | Consider the hybrid active filter parameter optimization configuration of resonance | |
| CN105305392A (en) | Symmetrical component method for short circuit calculation of voltage-controlled type IIDG included power distribution network | |
| CN105896558A (en) | VSC-based UPFC electromechanical transient modular modeling method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| ASS | Succession or assignment of patent right |
Owner name: STATE ELECTRIC NET CROP. Effective date: 20130424 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20130424 Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant after: China Electric Power Research Institute Applicant after: State Grid Corporation of China Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant before: China Electric Power Research Institute |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |