CN107147112B - A kind of power system modeling method and system based on weighting polymerization - Google Patents
A kind of power system modeling method and system based on weighting polymerization Download PDFInfo
- Publication number
- CN107147112B CN107147112B CN201710373744.1A CN201710373744A CN107147112B CN 107147112 B CN107147112 B CN 107147112B CN 201710373744 A CN201710373744 A CN 201710373744A CN 107147112 B CN107147112 B CN 107147112B
- Authority
- CN
- China
- Prior art keywords
- model
- discrete
- coefficient
- power grid
- subsystem
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000006116 polymerization reaction Methods 0.000 title claims abstract description 20
- 238000005094 computer simulation Methods 0.000 title claims abstract description 12
- 238000004458 analytical method Methods 0.000 claims abstract description 20
- 238000005070 sampling Methods 0.000 claims description 30
- 230000005611 electricity Effects 0.000 claims description 9
- 230000005284 excitation Effects 0.000 claims description 8
- 230000001052 transient effect Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000011426 transformation method Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract description 4
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 102100029390 CMRF35-like molecule 1 Human genes 0.000 description 1
- 102100029380 CMRF35-like molecule 2 Human genes 0.000 description 1
- 101000990055 Homo sapiens CMRF35-like molecule 1 Proteins 0.000 description 1
- 101000870135 Mus musculus Cytohesin-1 Proteins 0.000 description 1
- 101000855847 Mus musculus Cytohesin-3 Proteins 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 108010036356 cytohesin-2 Proteins 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000013316 zoning Methods 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- 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]
-
- 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)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Power Engineering (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The present invention discloses a kind of power system modeling method and system based on weighting polymerization, and step is to obtain each subsystem model under regional power grid same level to be modeled first；Then it is derived by corresponding discrete model, analyze to obtain the relationship of each model coefficient and each component parameters in regional power grid in discrete model again, and then according to the physical relation feature between model coefficient each in discrete model, judge whether each subsystem discrete model is effective；It is then based on each effective discrete model of subsystem under same level, the unified discrete model under same level is established using weighting polymerization；Ultimate analysis obtains the physical relation feature between each coefficient of unified discrete model of S4 foundation, and judge whether unified discrete model is correct, it is effective to unify discrete model for established regional power grid if correctly, can be used for subsequent model splicing etc. and applies.Present invention can apply to different compositions, the regional power grid of differing complexity, to facilitate the foundation of simulation model and the implementation of model splicing.
Description
Technical field
The present invention relates to power system modeling technical field, especially a kind of power system modeling side based on weighting polymerization
Method and system.
Background technique
China carries out the power grid control management mode of layering and zoning, and the Energy Management System of control centres at different levels is generally only built
The detailed model of power grid in this control centre's compass of competency is found, and adjacent area power grid generallys use Equivalent Model, it is complete to improve
The calculating speed of power grid.
The equivalence of regional power grid can be simple duty value or equivalent generator, also may include interconnection, etc.
The equivalent network being worth including plant stand, equivalent generator, duty value.In currently used equivalence method, Dynamic Equivalence with
The actual physics problem studied after equivalence to system is closely coupled, and practical Dynamic Equivalence can be divided into three classes: the people having the same aspiration and interest
The generator that same frequency and similar angle are waved is turned to one that is, according to the transient stability analysis of electric system by succusion
Group；Modal method based on characteristic value reduces the order of dynamic model using the number of compression mould, so that reduction is to calculating
Demand；Evaluation method method, the main online survey that equivalence is exported by estimating using interconnection information are distinguished.
For the regional power grid of different compositions, by Dynamic Equivalence, equivalent regional power grid model structure obtained
It may be different.In model splicing, if unified mould all can be used for different compositions, the regional power grid of differing complexity
Type structure describes, and is beneficial to the foundation of simulation model and the implementation of model splicing.
Explanation of nouns
CLM (Classic Load Models, classical load model), CLM model structure is refering to what is shown in Fig. 1, it is electricity
Master pattern in Force system comprehensive stability calculation procedure.Do not consider directly in use, the equivalent impedance of power distribution network is practical,
But be added on the equivalent stator impedance of motor, so being a kind of model structure for considering power distribution network indirectly.
Summary of the invention
The technical problem to be solved in the present invention are as follows: for different compositions, the regional power grid of differing complexity, saved with boundary
The voltage and current of point is state variable, describes regional power grid using unified discrete model, with facilitate simulation model foundation and
The implementation of model splicing.
A kind of the technical solution used in the present invention are as follows: power system modeling method based on weighting polymerization, comprising:
S1 obtains each subsystem model under regional power grid same level to be modeled；
S2 is derived by the model system of corresponding discrete model and discrete load model based on the model that S1 is obtained
Number；
S3, analysis obtain the relationship of each model coefficient and each component parameters in regional power grid in discrete model, and then basis
Physical relation feature in discrete model between each model coefficient judges whether each subsystem discrete model is effective；
S4 is established under same level based on each effective discrete model of subsystem under same level using weighting polymerization
Unified discrete model containing multiple subsystems；
S5, analysis obtain the physical relation feature between each coefficient of unified discrete model of S4 foundation, unified discrete with judgement
Whether model is correct, and it is effective to unify discrete model for established regional power grid if correct.
Load in power grid according to the present invention mainly includes induction-motor load and static load.Thus the present invention first
Based on classical load model, its discrete model is derived, obtains analytic parameter (the i.e. mould of discrete model for synchronizing classical load model
Type coefficient)；On the basis of classical load discrete model, the inner link of parameter of analytic model and each element of electric system, with
And the physical features between each model coefficient；Based on the discrete model of classical load, established under same level using weighting polymerization
Unified discrete model；On the basis of unified discrete model, the physical features between unified each coefficient of discrete model are analyzed；It will
Model coefficient meets the unified discrete model of relationship characteristic requirement as regional power grid valid model, is used for actual model splicing
Deng in application.
Further, each subsystem model is defined for classical load model, in classical load modelFor transient potential；
For voltage；For electric current；ω is rotor velocity；T′d0Time constant is wound for rotor；RsFor stator resistance；RrFor rotor electricity
Resistance；X is steady-state reactance；X ' is transient state reactance；XsFor stator reactance；XrFor rotor reactance；XmFor excitation reactance；
Then classical load model indicates are as follows:
In formula:X=Xs+Xm；X '=XsXm/(Xs+Xm)；Td0′
=(Xr+Xm)/Rr；
Based on formula (1), it is derived by using voltage, electric current as the load model of input/output relation are as follows:
Wherein: Ar=-(1+B Δ X)/T 'd0；Br=ω -1-G Δ X/T 'd0；Cr=G/T 'd0-B(ω-1)；Aj=-ω+1+
GΔX/T′d0；Bj=-(1+B Δ X)/T 'd0；Cj=[[B/T 'd0+G(ω-1)]]；Δ X=X-X '；G=RS/(RS 2+X′2)；B=
X′/(RS 2+X′2)；
I1rFor CLM model current real part；I1jFor CLM model current imaginary part；
The formula (2) of incremental form is subjected to Laplace transformation, obtains the load model under frequency domain are as follows:
In formula (3): S is laplace operator；
According to formula (3), the transmission function of the real part and imaginary part of electric current relative to port voltage is obtained:
Using bilinear transformation method, formula (4) and (5) are transformed to difference equation model:
ΔI1r(k+2)=θ11ΔI1r(k+1)+θ12ΔI1r(k)+θ13ΔU(k+2)+θ14ΔU(k+1)+θ15ΔU(k)
(6)
ΔI1j(k+2)=θ16ΔI1j(k+1)+θ17ΔI1j(k)+θ18ΔU(k+2)+θ19ΔU(k+1)+θ110ΔU(k)
(7)
Formula (6) and formula (7) are the discrete model of classical load model；In formula, k refer to discrete time (k=1,2,3,
4……)；The real part coefficient of discrete model are as follows:
The imaginary part coefficient of discrete model are as follows:
Wherein h is sampling step length, and z is the transform factor.
From the above equation, we can see that the real part coefficient (θ of the discrete model of CLM model1, θ12, θ13, θ14, θ15) or imaginary part coefficient
(θ16, θ17, θ18, θ19, θ110) and rotor velocity ω, rotor winding time constant T 'd0, stator resistance Rs, rotor resistance Rr, fixed
Sub- reactance Xs, rotor reactance Xr, excitation reactance XmAnd sampling step length h is related, illustrates that discrete model and each element of system exist
Connection, has certain mechanism meaning.
The port voltage for selecting each CLM model boundary node is input quantity, and the port current of model boundary node is output
Amount；In the corresponding discrete model of each classics load model, if when sampling step length h level off to 0 when, the coefficient of discrete model output item
The sum of level off to 1, the sum of coefficient of input item levels off to 0, then discrete model is effective.
I.e. when sampling step length h very little to level off to 0 when, there are following relationships for each coefficient:
As can be seen from the above equation when sampling step length h is less than or equal to 0.01s, the sum of coefficient of CLM model output item is close
Approximately equal to 1；The sum of the coefficient of CLM mode input item is approximately equal to 0.This conclusion can be the parameter identification and model of model
Whether effectively judgement provides reference frame.
The regional power grid containing classical load model is derived below unifies discrete model.
S4 comprising steps of
S41 seeks the incremental form expression formula of each subsystem classics load model median generatrix electric current:
DefinitionFor bus current, IirFor CLM model current real part, IijFor CLM model current imaginary part, n is same level
The quantity of lower subsystem；According to bus current formula:
Obtain incremental form are as follows:
S42, it is assumed that for regional power grid, in steady state equilibrium point, the electric current real and imaginary parts of classical load model are deposited
In following weight relationship:
In formula: KrFor the weight factor of subsystem electric current real part each under same level；KjFor subsystem electricity each under same level
Flow the weight factor of imaginary part；
S43, then convolution (6) and formula (7), the unification of formula (12) and formula (13) and formula (14) attainable region domain power grid from
Dissipate model are as follows:
ΔIr(k+2)=θ1ΔIr(k+1)+θ2ΔIr(k)+θ3ΔU(k+2)+θ4ΔU(k+1)+θ5ΔU(k) (8)
ΔIj(k+2)=θ6ΔIj(k+1)+θ7ΔIj(k)+θ8ΔU(k+2)+θ9ΔU(k+1)+θ10ΔU(k) (9)
In formula, θ1~θ5It respectively represents regional power grid and unifies discrete model real part coefficient, θ6~θ10Respectively represent regional power grid
Unified discrete model imaginary part coefficient；
As it can be seen that the model coefficient and rotor velocity ω of unified discrete model, rotor wind time constant T 'd0, stator electricity
Hinder Rs, rotor resistance Rr, stator reactance Xs, rotor reactance Xr, excitation reactance Xm, sampling step length and weight factor are related.
The port voltage of selection region power grid boundary node is input quantity, and the port current of regional power grid boundary node is defeated
Output；Then when sampling step length level off to 0 when, the sum of the coefficient of unified discrete model output item is approximately equal to 1, the coefficient of input item
The sum of be approximately equal to 0.
That is convolution (6)-(7) and formula (13)-(14) and formula (15)-(16), when sampling step length h is less than or equal to
When 0.01s, under certain condition, being derived from each coefficient, there are following relationships:
θ1+θ2≈1 (10)
θ3+θ4+θ5≈0 (11)
θ6+θ7≈1 (12)
θ8+θ9+θ10≈0 (13)
Comprehensive classics load model compositing area power grid unifies the equation of discrete model, it will thus be seen that 1) regional power grid mould
The model equation structure of type and classical load model is unified；2) model coefficient is all related with sampling step length, therefore sampling step length is
Key factor in modeling process；3) when sampling step length h very little, relationship that model coefficient still has are as follows: the sum of output item
It is approximately equal to 1, the sum of input item is approximately equal to 0.
Invention additionally discloses a kind of power system modeling systems based on weighting polymerization comprising:
Subsystem model obtains module, obtains each subsystem model under regional power grid same level to be modeled；
Subsystem discrete model obtains module, is based on each subsystem model, is derived by corresponding discrete model, Yi Jili
Dissipate the model coefficient of model；
Subsystem discrete model efficiency analysis module, analysis obtain in discrete model in each model coefficient and regional power grid
The relationship of each component parameters, and then according to the physical relation feature between model coefficient each in discrete model, judge each subsystem
Whether discrete model is effective；
Regional power grid unifies discrete model and establishes module, based on each effective discrete model of subsystem under same level, adopts
The unified discrete model under same level is established with weighting polymerization；
Unified discrete model efficiency analysis module, analysis obtain the physics between established unified each coefficient of discrete model
Relationship characteristic, to judge whether unified discrete model is correct, and it is effective to unify discrete model for established regional power grid if correct.
Beneficial effect
Compared with prior art, the present invention has following progress:
Firstly, being related using technical solution of the present invention model built with system component parameter, anticipate with certain mechanism
Justice, while picking out the parameter come and can be very good to verify curve itself, have certain adaptability；
Secondly, the physical features between discrete model coefficient can provide reference frame for the parameter identification result of system model；
By being described using unified model structure in power system modeling, can facilitate the foundation of simulation model with
And the implementation of model splicing.
Detailed description of the invention
Fig. 1 is classical load model figure；
Fig. 2 is the structural block diagram containing multiple classical load models under same level；
Fig. 3 is inventive algorithm schematic illustration；
Fig. 4 is CEPRI-36 node example of the present invention；
Fig. 5 is 3% curve graph of Voltage Drop；
Fig. 6 is 30% curve graph of Voltage Drop；
Fig. 7 is 3% electric current real part curve matching figure of Voltage Drop；
Fig. 8 is 3% electric current imaginary part curve matching figure of Voltage Drop；
Fig. 9 is 30% electric current real part curve matching figure of Voltage Drop；
Figure 10 is 30% electric current imaginary part curve matching figure of Voltage Drop.
Specific embodiment
It is further described below in conjunction with the drawings and specific embodiments.
Power system modeling method based on weighting polymerization of the invention, comprising:
S1 obtains each subsystem model under regional power grid same level to be modeled；
S2 obtains model based on S1, is derived by the model coefficient of corresponding discrete model and discrete model；
S3, analysis obtain the relationship of each model coefficient and each component parameters in regional power grid in discrete model, and then basis
Physical relation feature in discrete model between each model coefficient judges whether each subsystem discrete model is correct；
S4 is established under same level based on each correct discrete model of subsystem under same level using weighting polymerization
Unified discrete model；
S5, analysis obtain the physical relation feature between each coefficient of unified discrete model of S4 foundation, unified discrete with judgement
Whether model is correct, and it is effective to unify discrete model for established regional power grid if correct.
Load in power grid according to the present invention mainly includes induction-motor load and static load.Thus the present invention first
Based on classical load model, its discrete model is derived, obtains the analytic parameter for synchronizing classical load model；It is discrete in classical load
On the basis of model, the physics between parameter of analytic model and the inner link and each model coefficient of each element of electric system is special
Sign；Based on the discrete model of classical load, the unified discrete model under same level is established using weighting polymerization；It is unified from
On the basis of dissipating model, the physical features between unified each coefficient of discrete model are analyzed；Model coefficient is met relationship characteristic to want
The unified discrete model asked is as regional power grid valid model, in the application such as actual model splicing.
Each subsystem model is defined for classical load model, in classical load modelFor transient potential；For voltage；
For electric current；ω is rotor velocity；T′d0Time constant is wound for rotor；RsFor stator resistance；RrFor rotor resistance；X is stable state
Reactance；X ' is transient state reactance；XsFor stator reactance；XrFor rotor reactance；XmFor excitation reactance；
Then the classical load model of each subsystem indicates are as follows:
In formula:X=Xs+Xm；X '=XsXm/(Xs+Xm)；Td0'=
(Xr+Xm)/Rr；
Based on formula (1), it is derived by using voltage, electric current as the load model of input/output relation are as follows:
Wherein: Ar=-(1+B Δ X)/T 'd0；Br=ω -1-G Δ X/T 'd0；Cr=G/T 'd0-B(ω-1)；Aj=-ω+1+
GΔX/T′d0；Bj=-(1+B Δ X)/T 'd0；Cj=[B/T 'd0+G(ω-1)]；Δ X=X-X '；G=RS/(RS 2+X′2)；B=
X′/(RS 2+X′2)；
I1rFor CLM model current real part；I1jFor CLM model current imaginary part；
The formula (2) of incremental form is subjected to Laplace transformation, obtains the load model under frequency domain are as follows:
In formula (3): S is laplace operator；
According to formula (3), the transmission function of the real part and imaginary part of electric current relative to port voltage is obtained:
Using bilinear transformation method, formula (4) and (5) are transformed to difference equation model:
ΔI1r(k+2)=θ11ΔI1r(k+1)+θ12ΔI1r(k)+θ13ΔU(k+2)+θ14ΔU(k+1)+θ15ΔU(k)
(6)
ΔI1j(k+2)=θ16ΔI1j(k+1)+θ17ΔI1j(k)+θ18ΔU(k+2)+θ19ΔU(k+1)+θ110ΔU(k)
(7)
Formula (6) and formula (7) are the discrete model of classical load model；In formula, k refer to discrete time (k=1,2,3,
4……)；The real part coefficient of discrete model are as follows:
The imaginary part coefficient of discrete model are as follows:
Wherein h is sampling step length, and z is the transform factor.
From the above equation, we can see that the real part coefficient (θ of the discrete model of CLM model1, θ12, θ13, θ14, θ15) or imaginary part coefficient
(θ16, θ17, θ18, θ19, θ110) and rotor velocity ω, rotor winding time constant T 'd0, stator resistance Rs, rotor resistance Rr, fixed
Sub- reactance Xs, rotor reactance Xr, excitation reactance XmAnd sampling step length h is related, illustrates that discrete model and each element of system exist
Connection, has certain mechanism meaning.
The port voltage for selecting each CLM model boundary node is input quantity, and the port current of model boundary node is output
Amount；In the corresponding discrete model of each classics load model, if when sampling step length h level off to 0 when, the coefficient of discrete model output item
The sum of level off to 1, the sum of coefficient of input item levels off to 0, then discrete model is effective.
I.e. when sampling step length h very little to level off to 0 when, there are following relationships for each coefficient:
As can be seen from the above equation when sampling step length h is less than or equal to 0.01s, the sum of coefficient of CLM model output item is close
Approximately equal to 1；The sum of the coefficient of CLM mode input item is approximately equal to 0.This conclusion can be the parameter identification and model of model
Whether effectively judgement provides reference frame.
The regional power grid containing classical load model is derived below unifies discrete model.
S4 comprising steps of
S41 seeks the incremental form expression formula of each subsystem classics load model median generatrix electric current:
DefinitionFor bus current, IirFor CLM model current real part, IijFor CLM model current imaginary part, n is same level
The quantity of lower subsystem；According to bus current formula:
Obtain incremental form are as follows:
S42, it is assumed that for regional power grid, in steady state equilibrium point, the electric current real and imaginary parts of classical load model are deposited
In following weight relationship:
In formula: KrFor the weight factor of subsystem electric current real part each under same level；KjFor subsystem electricity each under same level
Flow the weight factor of imaginary part；
S43, then convolution (6) and formula (7), the unification of formula (12) and formula (13) and formula (14) attainable region domain power grid from
Dissipate model are as follows:
ΔIr(k+2)=θ1ΔIr(k+1)+θ2ΔIr(k)+θ3ΔU(k+2)+θ4ΔU(k+1)+θ5ΔU(k) (21)
ΔIj(k+2)=θ6ΔIj(k+1)+θ7ΔIj(k)+θ8ΔU(k+2)+θ9ΔU(k+1)+θ10ΔU(k) (22)
In formula, θ1~θ5It respectively represents regional power grid and unifies discrete model real part coefficient, θ6~θ10Respectively represent regional power grid
Unified discrete model imaginary part coefficient；
As it can be seen that the model coefficient and rotor velocity ω of unified discrete model, rotor wind time constant T 'd0, stator electricity
Hinder Rs, rotor resistance Rr, stator reactance Xs, rotor reactance Xr, excitation reactance Xm, sampling step length and weight factor are related.
The port voltage of selection region power grid boundary node is input quantity, and the port current of regional power grid boundary node is defeated
Output；Then when sampling step length level off to 0 when, the sum of the coefficient of unified discrete model output item is approximately equal to 1, the coefficient of input item
The sum of be approximately equal to 0.
That is convolution (6)-(7) and formula (13)-(14) and formula (15)-(16), when sampling step length h is less than or equal to
When 0.01s, under certain condition, being derived from each coefficient, there are following relationships:
θ1+θ2≈1 (23)
θ3+θ4+θ5≈0 (24)
θ6+θ7≈1 (25)
θ8+θ9+θ10≈0 (26)
Comprehensive classics load model compositing area power grid unifies the equation of discrete model, it will thus be seen that 1) regional power grid mould
The model equation structure of type and classical load model is unified；2) model coefficient is all related with sampling step length, therefore sampling step length is
Key factor in modeling process；3) when sampling step length h very little, relationship that model coefficient still has are as follows: the sum of output item
It is approximately equal to 1, the sum of input item is approximately equal to 0.
Embodiment
The present embodiment verifies the validity of the method for the present invention by CEPRI-36 node example system.
It chooses and is implemented for containing only the regional power grid of classical load model below.Fig. 4 is CEPRI-36 node example,
Assuming that accessing classical load model in regional power grid at BUS50, classical load model is made of induction-motor load and static load.
CLM1 Selection Model parameter are as follows:
Rs=0, Rr=0.02, Xs=0.18, Xr=0.12, Xm=3.5, h=0.01.
CLM2 Selection Model parameter are as follows:
Rs=0, Rr=0.02, Xs=0.181, Xr=0.12, Xm=3.5, h=0.01.
CLM3 Selection Model parameter are as follows:
Rs=0, Rr=0.02, Xs=0.295, Xr=0.12, Xm=3.5, h=0.01.
Required data are modeled in order to obtain, apply disturbance outside regional power grid: at 1 second, BUS16-BUS29 route
Single phase grounding fault is arranged in the side BUS16, cuts off after 0.2s.Setting sampling step length is 0.01s, and soft using PSASP
Input voltage, active power and the reactive power of part emulation record CLM model.
Using method of the invention, follow the steps below:
S1 obtains three subsystems model under same level；
S2 is derived by the discrete model of each subsystem；
S3 analyzes the relationship in discrete model between model coefficient and model parameter, that is, electric network element parameter, and then obtains
Physical relation feature between each discrete model model coefficient, to judge whether each discrete model is effective；
The port voltage for selecting each CLM model boundary node is input quantity, and the port current of model boundary node is output
Amount；In the corresponding discrete model of each classics load model, if the model coefficient of discrete model meets physical relation feature: working as sampling
Step-length h level off to 0 when, the sum of the coefficient of discrete model output item levels off to 1, and the sum of coefficient of input item levels off to 0, then discrete
Model is effective.
S4 is based on effective three subsystems discrete model, using weighting polymerization to the region electricity under same level
Net carries out the modeling of unified discrete model；
S5 analyzes the physical relation feature of model coefficient in unified discrete model, the end of selection region power grid boundary node
Mouth voltage is input quantity, and the port current of regional power grid boundary node is output quantity；Then model coefficient meets physical relation feature:
When sampling step length level off to 0 when, the sum of the coefficient of unified discrete model output item is approximately equal to 1, and the sum of coefficient of input item is close
Approximately equal to 0, then it is effective to unify discrete model.
The following table 1 is under different Voltage Drops, between the parameter and discrete model coefficient of regional power grid discrete model
Physical features.
1 regional power grid discrete model results of model parameter identification IrIj-U of table
Physical features between identification of Model Parameters and discrete model coefficient all demonstrate the validity of this patent.
The curve matching of Fig. 5 to Figure 10 also demonstrates the practicability and validity of the method for the present invention.
The present invention is based on the power system modeling systems for weighting polymerization to include:
Subsystem model obtains module, obtains each subsystem model under regional power grid same level to be modeled；
Subsystem discrete model obtains module, is based on each subsystem model, is derived by corresponding discrete model, Yi Jili
Dissipate the model coefficient of model；
Subsystem discrete model efficiency analysis module, analysis obtain in discrete model in each model coefficient and regional power grid
The relationship of each component parameters, and then according to the physical relation feature between model coefficient each in discrete model, judge each subsystem
Whether discrete model is effective；
Regional power grid unifies discrete model and establishes module, based on each effective discrete model of subsystem under same level, adopts
The unified discrete model under same level is established with weighting polymerization；
Unified discrete model efficiency analysis module, analysis obtain the physics between established unified each coefficient of discrete model
Relationship characteristic, to judge whether unified discrete model is correct, and it is effective to unify discrete model for established regional power grid if correct.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Claims (6)
1. a kind of power system modeling method based on weighting polymerization, characterized in that include:
S1 obtains each subsystem model under regional power grid same level to be modeled；
S2 is derived by the model coefficient of corresponding discrete model and discrete load model based on the model that S1 is obtained；
S3, analysis obtain the relationship of each model coefficient and each component parameters in regional power grid in discrete model, and then according to discrete
Physical relation feature in model between each model coefficient judges whether each subsystem discrete model is effective: preference pattern boundary
The port voltage of node is input quantity, and the port current of model boundary node is output quantity；It is each classics load model it is corresponding from
Dissipate in model, if when sampling step length h level off to 0 when, the sum of the coefficient of discrete model output item levels off to 1, the coefficient of input item
The sum of level off to 0, then discrete model is effective；
S4 is established under same level using weighting polymerization and is contained based on each effective discrete model of subsystem under same level
The unified discrete model of multiple subsystems；
S5, analysis obtain the physical relation feature between each coefficient of unified discrete model of S4 foundation, to judge unified discrete model
Whether correct, it is effective unify discrete model for established regional power grid if correct.
2. according to the method described in claim 1, it is characterized in that, based on S1 obtain model, be derived by corresponding discrete model
Step includes:
Each subsystem model is defined for classical load model, in classical load modelFor transient potential；For voltage；For electricity
Stream；ω is rotor velocity；T′d0Time constant is wound for rotor；RsFor stator resistance；RrFor rotor resistance；X is stable state electricity
It is anti-；X ' is transient state reactance；XsFor stator reactance；XrFor rotor reactance；XmFor excitation reactance；
Then classical load model indicates are as follows:
In formula:X=Xs+Xm；X '=XsXm/(Xs+Xm)；Td0'=(Xr+
Xm)/Rr；
Based on formula (1), it is derived by using voltage, electric current as the load model of input/output relation are as follows:
Wherein: Ar=-(1+B Δ X)/T 'd0；Br=ω -1-G Δ X/T 'd0；Cr=G/T 'd0-B(ω-1)；Aj=-ω+1+G Δ X/
T′d0；Bj=-(1+B Δ X)/T 'd0；Cj=[B/T 'd0+G(ω-1)]；Δ X=X-X '；G=RS/(RS 2+X′2)；B=X '/(RS 2+
X′2)；
I1rFor CLM model current real part；I1jFor CLM model current imaginary part；
The formula (2) of incremental form is subjected to Laplace transformation, obtains the load model under frequency domain are as follows:
In formula (3): S is laplace operator；
According to formula (3), the transmission function of the real part and imaginary part of electric current relative to port voltage is obtained:
Using bilinear transformation method, formula (4) and (5) are transformed to difference equation model:
ΔI1r(k+2)=θ11ΔI1r(k+1)+θ12ΔI1r(k)+θ13ΔU(k+2)+θ14ΔU(k+1)+θ15ΔU(k) (6)
ΔI1j(k+2)=θ16ΔI1j(k+1)+θ17ΔI1j(k)+θ18ΔU(k+2)+θ19ΔU(k+1)+θ110ΔU(k) (7)
Formula (6) and formula (7) are the discrete model of classical load model；In formula, k refers to discrete time (k=1,2,3,4 ...)；
The real part coefficient of discrete model are as follows:
The imaginary part coefficient of discrete model are as follows:
Wherein h is sampling step length, and z is the transform factor.
3. according to the method described in claim 2, it is characterized in that, S4 comprising steps of
S41 seeks the incremental form expression formula of each subsystem classics load model median generatrix electric current:
DefinitionFor bus current, IirFor CLM model current real part, IijFor CLM model current imaginary part, n is subsystem under same level
The quantity of system；According to bus current formula:
Obtain incremental form are as follows:
S42, it is assumed that for regional power grid, in steady state equilibrium point, the electric current real and imaginary parts of classical load model exist such as
Lower weight relationship:
In formula: KrFor the weight factor of subsystem electric current real part each under same level；KjIt is empty for subsystem electric current each under same level
The weight factor in portion；
S43, then convolution (6) and formula (7), the unified walk-off-mode of formula (12) and formula (13) and formula (14) attainable region domain power grid
Type are as follows:
ΔIr(k+2)=θ1ΔIr(k+1)+θ2ΔIr(k)+θ3ΔU(k+2)+θ4ΔU(k+1)+θ5ΔU(k) (15)
ΔIj(k+2)=θ6ΔIj(k+1)+θ7ΔIj(k)+θ8ΔU(k+2)+θ9ΔU(k+1)+θ10ΔU(k) (16)
In formula, θ1~θ5It respectively represents regional power grid and unifies discrete model real part coefficient, θ6~θ10It is unified to respectively represent regional power grid
Discrete model imaginary part coefficient；
4. method according to claim 1 or 3, characterized in that the model coefficient and rotor velocity of unified discrete model
ω, rotor wind time constant T 'd0, stator resistance Rs, rotor resistance Rr, stator reactance Xs, rotor reactance Xr, excitation reactance Xm,
Sampling step length and weight factor are related.
5. according to the method described in claim 4, it is characterized in that, the port voltage of selection region power grid boundary node is input
Amount, the port current of regional power grid boundary node are output quantity；Then when sampling step length level off to 0 when, unified discrete model output
The sum of the coefficient of item is approximately equal to 1, and the sum of coefficient of input item is approximately equal to 0.
6. a kind of power system modeling system based on weighting polymerization, characterized in that include:
Subsystem model obtains module, obtains each subsystem model under regional power grid same level to be modeled；
Subsystem discrete model obtains module, is based on each subsystem model, is derived by corresponding discrete model and walk-off-mode
The model coefficient of type；
Subsystem discrete model efficiency analysis module, analysis obtain each model coefficient and each member in regional power grid in discrete model
The relationship of part parameter, and then according to the physical relation feature between model coefficient each in discrete model, judge that each subsystem is discrete
Whether model is effective: the port voltage of preference pattern boundary node is input quantity, and the port current of model boundary node is output
Amount；In the corresponding discrete model of each classics load model, if when sampling step length h level off to 0 when, the coefficient of discrete model output item
The sum of level off to 1, the sum of coefficient of input item levels off to 0, then discrete model is effective；
Regional power grid unifies discrete model and establishes module, based on each effective discrete model of subsystem under same level, using adding
Power polymerization establishes the unified discrete model under same level；
Unified discrete model efficiency analysis module, analysis obtain the physical relation between established unified each coefficient of discrete model
Feature, to judge whether unified discrete model is correct, and it is effective to unify discrete model for established regional power grid if correct.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710373744.1A CN107147112B (en) | 2017-05-24 | 2017-05-24 | A kind of power system modeling method and system based on weighting polymerization |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710373744.1A CN107147112B (en) | 2017-05-24 | 2017-05-24 | A kind of power system modeling method and system based on weighting polymerization |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107147112A CN107147112A (en) | 2017-09-08 |
| CN107147112B true CN107147112B (en) | 2019-09-24 |
Family
ID=59779782
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710373744.1A Expired - Fee Related CN107147112B (en) | 2017-05-24 | 2017-05-24 | A kind of power system modeling method and system based on weighting polymerization |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107147112B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113949059A (en) * | 2021-09-28 | 2022-01-18 | 昆明理工大学 | Generalized discrete equivalent model building method of comprehensive load |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010246289A (en) * | 2009-04-07 | 2010-10-28 | Hiroshima Univ | Critical failure removing period calculation method, and program |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100444494C (en) * | 2006-09-29 | 2008-12-17 | 天津大学 | Coherency generating sets of electric power system and polymeric method of asynchronous motor group parameter |
| CN101882167B (en) * | 2009-11-25 | 2013-02-06 | 中国电力科学研究院 | Wind power station equivalent modeling method of large-scale wind power concentration access power grid |
| CN101789598B (en) * | 2010-03-05 | 2012-05-30 | 湖北省电力试验研究院 | Power system load modelling method |
| CN102760191B (en) * | 2012-06-26 | 2014-09-03 | 西安交通大学 | Double-fed machine set wind power station equivalent modeling system and method based on rotation speed grouping |
| US9563722B2 (en) * | 2012-11-13 | 2017-02-07 | Gridquant, Inc. | Sigma algebraic approximants as a diagnostic tool in power networks |
| CN105356454B (en) * | 2015-11-05 | 2017-11-21 | 上海电力学院 | A kind of power system load modeling method based on typical load group |
-
2017
- 2017-05-24 CN CN201710373744.1A patent/CN107147112B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010246289A (en) * | 2009-04-07 | 2010-10-28 | Hiroshima Univ | Critical failure removing period calculation method, and program |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107147112A (en) | 2017-09-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gurrala et al. | Parareal in time for fast power system dynamic simulations | |
| Osipov et al. | Adaptive nonlinear model reduction for fast power system simulation | |
| Milano et al. | A systematic method to model power systems as stochastic differential algebraic equations | |
| Shu et al. | A novel interfacing technique for distributed hybrid simulations combining EMT and transient stability models | |
| CN112636341B (en) | Power system inertia spatial distribution estimation method and device based on multiple innovation identification | |
| Allen et al. | Interactive object-oriented simulation of interconnected power systems using SIMULINK | |
| Jalili-Marandi et al. | Real-time simulation of grid-connected wind farms using physical aggregation | |
| CN103810646B (en) | Improved projection integral algorithm based active power distribution system dynamic simulation method | |
| CN110135031A (en) | Electric power system transient stability calculation method based on half implicit runge kutta method | |
| Tomim et al. | Parallel transient stability simulation based on multi-area Thévenin equivalents | |
| CN109145362A (en) | A kind of power network modeling method and system | |
| US11569682B2 (en) | System and method for a fast power network simulator | |
| Jardim et al. | A methodology for computing robust dynamic equivalents of large power systems | |
| CN112380670A (en) | Virtual rotor-based modeling method and system for sectional power supply linear induction motor | |
| Abusalah et al. | CPU based parallel computation of electromagnetic transients for large power grids | |
| Irwin et al. | Parallel processing and hybrid simulation for HVDC/VSC PSCAD studies | |
| CN107147112B (en) | A kind of power system modeling method and system based on weighting polymerization | |
| Plumier | Co-simulation of electromagnetic transients and phasor models of electric power systems | |
| Kyesswa et al. | A Matlab-based dynamic simulation module for power system transients analysis in the eASiMOV framework | |
| Ourari et al. | Generating units aggregation for dynamic equivalent of large power systems | |
| Almeida et al. | A software tool for the determination of dynamic equivalents of power systems | |
| Salam et al. | Transient Stability Analysis of a Three-machine Nine Bus Power System Network. | |
| Montier et al. | Proper generalized decomposition applied on a rotating electrical machine | |
| US20140088894A1 (en) | Method and device for extracting skeleton topology structure of electric power grid | |
| Teja et al. | Power network layout generation using force directed graph technique |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190924 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |