CN110336299A - A kind of Distribution system method considering integrated energy system small interference stability - Google Patents
A kind of Distribution system method considering integrated energy system small interference stability Download PDFInfo
- Publication number
- CN110336299A CN110336299A CN201910514943.9A CN201910514943A CN110336299A CN 110336299 A CN110336299 A CN 110336299A CN 201910514943 A CN201910514943 A CN 201910514943A CN 110336299 A CN110336299 A CN 110336299A
- Authority
- CN
- China
- Prior art keywords
- small interference
- power supply
- state
- integrated energy
- interference stability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000005611 electricity Effects 0.000 claims abstract description 24
- 230000001360 synchronised effect Effects 0.000 claims description 32
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 30
- 230000008569 process Effects 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 16
- 239000003345 natural gas Substances 0.000 claims description 15
- 238000005516 engineering process Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 10
- 230000003534 oscillatory effect Effects 0.000 claims description 10
- 239000007789 gas Substances 0.000 claims description 6
- 230000006698 induction Effects 0.000 claims description 6
- 238000010438 heat treatment Methods 0.000 claims description 5
- 239000003990 capacitor Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 230000001052 transient effect Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000001816 cooling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000004804 winding 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
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- 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)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of Distribution system methods for considering integrated energy system small interference stability, mainly in region class integrated energy system, influence of the Distribution system to the small interference stability of power supplys all kinds of in system and whole system, by building respective small interference stability model to power supplys various types of in integrated energy system, and final integrated energy system model is formed by coordinate transform, and small interference stability is added in the constraint of distribution reconstruct, simultaneously by wind-powered electricity generation, the power output of photovoltaic is arranged in system and runs under most unstable scene, ensure integrated energy system on the basis of meeting small interference stability, more economical operation.
Description
Technical field:
The present invention relates to a kind of power transmission and transformation fields, consider matching for integrated energy system small interference stability more particularly to a kind of
Net reconstructing method.
Background technique:
With the development and the improvement of people's living standards of society, demand of the people to the energy is gradually increased.Continue out
Hair utilizes limited fossil energy, not only results in environmental pollution, also results in lack of energy, problem of energy crisis occurs.It opens
Hair is become using clean renewable energy, the utilization rate for improving the energy to solve the effective of environmental pollution and problem of energy crisis
Means.In recent years, researcher is to the research and application of the renewable energy power generations technology such as wind, light so that the above problem is had
The alleviation of effect.In addition to this, the proposition of integrated energy system is to the permeability and efficiency of energy utilization for improving renewable energy
Important in inhibiting, it is intended to break the barrier between each energy supply system, promotes the coupling of multifunctional system and coordinates excellent
Change operation.
Integrated energy system can be divided into transregional grade, region class and user class according to different utilization scenes.For by
The region class integrated energy system that distribution system, natural gas system and heat supply/cooling system are constituted.Wherein, distribution system contains
Renewable energy based on double-fed wind generating technology and photovoltaic DC-to-AC converter generation technology, natural gas system and heating system difference
It is coupled by gas turbine and cogeneration units with distribution system.Therefore, synchronization is contained in entire integrated energy system
Generator, asynchronous generator and inverter power generation apparatus, the dynamic process time scale span in system are big.Its is stable
Property can be divided into small interference stability and Large disturbance stability according to the size of disturbance, wherein the network of the characteristics of all kinds of power supplys, power distribution network
Structure is its major influence factors.
Under conditions of high proportion renewable energy is grid-connected, the small interference stability problem values of integrated energy system must cause to close
Note.For different types of power supply (such as photovoltaic, wind-powered electricity generation, gas turbine), due to its operation mechanism and the difference of control mode,
Its small signal stability is also not quite similar.It in addition to this, generally will not be as synchronous in power transmission network since these power supply capacities are smaller
Machine equally configures power system stabilizer, PSS, along with their power output has very strong uncertainty, so as to the small interference of system
Stable influence becomes increasingly complex.
Distribution system can reduce the loss of system, but the operation side of system is also changed while changing structure
Formula affects the small signal stability of system.When carrying out Distribution system, focus is placed on economy by most researchers
On, the dynamic process of all kinds of power supplys and system is had ignored, that is, has ignored the stability of system.
For Distribution system generally to reduce network loss, raising network voltage level etc. as objective function, constraint condition includes trend
Constraint, the radial constraint of network, Branch Power Flow constraint and active and reactive units limits.After network structure changes, respectively
The small signal stability of class power supply operation is necessarily influenced to different extents.If Distribution system only considers economy, then very
Be likely to cause certain power supplys can not stable operation the case where, lead to off-grid.In addition to this, wind-powered electricity generation, photovoltaic power output have
The uncertainty of randomness, power output inevitably results in the variation of trend and node voltage in system, to influence unit
Electromagnetic power so that unit is constantly subjected to the test of small interference stability in distribution, due to all kinds of power supplys small signal stability not
Together, the network structure after Distribution system can cause positive or negative effect to the small signal stability of system.The prior art
In be present to ensure that the Distribution system method of distributed generation resource small interference stability, but all regard all power supplys as synchronous machine, do not have
Have and various types of power supplys are distinguished.
Summary of the invention:
The technical problems to be solved by the present invention are: overcome the deficiencies in the prior art, to various types of in integrated energy system
The power supply of type establishes dynamic model, and small interference stability is added in the constraint of distribution reconstruct, it is ensured that integrated energy system is full
On the basis of sufficient small interference stability, the Distribution system method of more economical the considerations of running integrated energy system small interference stability.
The technical scheme is that a kind of Distribution system method for considering integrated energy system small interference stability, to comprehensive
It closes various types of power supplys in energy resource system and establishes dynamic model, and small interference stability is added in the constraint of distribution reconstruct, really
Integrated energy system is protected on the basis of meeting small interference stability, it is more economical to run, it the steps include: Step 1: to based on double
Present wind-powered electricity generation power supply, the photo-voltaic power supply based on inverter generation technology, the natural pneumoelectric based on gas turbine of wind generating technology
Respective small interference stability model is established as object in source respectively；
Step 2: the dq coordinate system of each model is transformed into common coordinate system DQ coordinate system by coordinate transform, carry out
Model integration；
Step 3: calculate initial network structure under each power supply operating status, and calculate under this condition the characteristic value of system,
Small interference stability is constrained and is linearized by feature vector；
Step 4: using loss minimization as objective function, and corresponding constraint condition is introduced, establish the reconstruct for considering small interference
Model；
Step 5: the small interference stability constraint after linearisation is added in reconstruction model, Distribution system scheme is solved；
Step 6: calculating the operating status of each power supply under the reconfiguration scheme, and the dominant characteristics value of computing system again；
Step 7: judging whether the dominant characteristics value of system meets small interference stability constraint；Meet, then exports reconstruct side
Case is unsatisfactory for, then using the reconfiguration scheme as primary condition, then return step three.
Further, in the step 1, the small interference of wind-powered electricity generation power supply is established with legacy system, induction machine and filter
Stability analysis model, and linearization process is carried out, the state-space expression of wind-powered electricity generation power supply can be obtained:
In formula, xDFIGFor the quantity of state of double-fed blower, uDFIGFor the control amount of double-fed blower, yDFIGFor the defeated of double-fed blower
Output, ADFIGFor the state matrix of double-fed blower, element therein by blower under the parameter and the operating status of blower state
Amount is constituted.
Further, it in the step 1, is established with voltage controller, current controller, filter, capacitor, route
Photo-voltaic power supply Small signal stability analysis model, and linearization process is carried out, the state-space expression of photo-voltaic power supply can be obtained:
In formula, xpvFor the quantity of state of photovoltaic, upvFor the control amount of photovoltaic, ypvFor the output quantity of photovoltaic, ApvFor the shape of photovoltaic
State matrix, element therein are made of the quantity of state of photovoltaic under the parameter of photovoltaic and the operating status.
Further, in the step 1, synchronous generator connects electrical system with natural gas system as heating system
Tie, using traditional third-order model as natural gas power supply Small signal stability analysis model, and carry out linearization process, can obtain
The state-space expression of photovoltaic plant:
In formula, xsynFor the quantity of state of synchronous machine, usynFor the control amount of synchronous machine, AsynFor the state matrix of synchronous machine,
In element be made of the quantity of state of synchronous machine under the parameter of synchronous machine and the operating status.
Further, in the step 2, the electric current that fan power is injected into system is stator current and net side inversion
Device exports the sum of electric current；The electric current that photo-voltaic power supply is injected to system is that filter exports electric current；Natural gas power supply is infused into system
The electric current entered is stator current；
Respective dq coordinate system is transformed into common coordinate system DQ coordinate system according to coordinate transform formula, coordinate transform is such as
Under:
Wherein, DQ coordinate system represents common coordinate system, diqiCoordinate system represents the coordinate system of each power supply, δiIt is two
Angle between a coordinate system；
Increased according to the voltage and current of photovoltaic in double-fed blower in synchronous machine in natural gas power supply, wind-powered electricity generation power supply, photo-voltaic power supply
Amount, and the state of the state-space expression of the state-space expression of convolution wind-powered electricity generation power supply, photo-voltaic power supply, photovoltaic plant is empty
Between expression formula can obtain the state-space expression of the integrated energy system:
In formula, the state matrix A of system is made of network, double-fed blower, photovoltaic, synchronous machine and their operating status.
Further, in the step 3, the shape of the system can be obtained by the state-space expression of integrated energy system
State space matrix A, the characteristic value of calculating matrix A can obtain the system by the Oscillatory mode shape after small interference；
It is leading Oscillatory mode shape with Oscillatory mode shape corresponding to the characteristic value near the imaginary axis, characteristic value is leading characteristic value
λ1；Dominant characteristics value is constrained to λref, then small interference stability constrains are as follows:
λ1≤λref (42)
According to the calculation formula of characteristic value:
| A- λ E |=0 (43)
Small interference stability is constrained and is linearized.
Further, in the step 4, with the objective function of loss minimization are as follows:
The constraint condition includes the radial constraint of network, power-balance constraint, voltage constrains, Branch Power Flow constrains and small
Interference stability constraint.
The beneficial effects of the present invention are:
1, the present invention is mainly in region class integrated energy system, and Distribution system is to power supplys all kinds of in system and entire system
The influence of the small interference stability of system, by establishing dynamic model to power supplys various types of in integrated energy system, and will be small dry
It disturbs to stablize and be added in the constraint of distribution reconstruct, it is ensured that integrated energy system is on the basis of meeting small interference stability, more economically
Operation.
2, small interference stability is constrained and is added in the constraint of Distribution system by the present invention, while the power output of wind-powered electricity generation, photovoltaic being set
It sets and is run under most unstable scene in system, so that the result of Distribution system can be in the base for guaranteeing system small interference stability
It is more economical to run on plinth.
Detailed description of the invention:
Fig. 1 is double-fed blower fan structure schematic diagram.
Fig. 2 is photovoltaic plant structural schematic diagram.
Fig. 3 is voltage, current controller schematic diagram.
Fig. 4 is the transformation of coordinate system.
Fig. 5 is generating set power characteristic.
Specific embodiment:
Embodiment: referring to Fig. 1, Fig. 2, Fig. 3, Fig. 4 and Fig. 5；In Fig. 1,1- blade, 2- gear-box, 3- transmission system, 4- sense
Answer motor, 5- current transformer, 6- filter inductance, 7- transformer and route, 8- bulk power grid.
The Distribution system method for considering integrated energy system small interference stability, to electricity various types of in integrated energy system
Dynamic model is established in source, and small interference stability is added in the constraint of distribution reconstruct, it is ensured that integrated energy system is meeting small do
It disturbs on the basis of stablizing, it is more economical to run.
Embodiment: according to different utilization scenes, integrated energy system is divided into transregional grade, region class and user class.
The embodiment of the present application uses the region class integrated energy system being made of distribution system, natural gas system and heat supply/cooling system.Its
In, distribution system contains the renewable energy based on double-fed wind generating technology and photovoltaic DC-to-AC converter generation technology, natural gas
System and heating system pass through gas turbine respectively and cogeneration units are coupled with distribution system.It the steps include: Step 1: right
Wind-powered electricity generation power supply based on double-fed wind generating technology, the photo-voltaic power supply based on inverter generation technology, the day based on gas turbine
Respective small interference stability model is established as object in right pneumoelectric source respectively；
Step 2: the dq coordinate system of each model is transformed into common coordinate system DQ coordinate system by coordinate transform, carry out
Model integration；
Step 3: calculate initial network structure under each power supply operating status, and calculate under this condition the characteristic value of system,
Small interference stability is constrained and is linearized by feature vector；
Step 4: using loss minimization as objective function, and corresponding constraint condition is introduced, establish the reconstruct for considering small interference
Model；
Step 5: the small interference stability constraint after linearisation is added in reconstruction model, Distribution system scheme is solved；
Step 6: calculating the operating status of each power supply under the reconfiguration scheme, and the dominant characteristics value of computing system again；
Step 7: judging whether the dominant characteristics value of system meets small interference stability constraint；Meet, then exports reconstruct side
Case is unsatisfactory for, then using the reconfiguration scheme as primary condition, then return step three.
Wind-powered electricity generation power supply Small signal stability analysis model foundation process is as follows:
Double-fed blower relies on its flexible control mode and lower operating cost, is widely used in each wind field,
Structure is as shown in Figure 1.Double-fed blower is made of transmission system, induction machine, current transformer, filter inductance, passes through transformer and line
Road is accessed in bulk power grid.The dynamic process of double-fed blower is complicated, and time scale difference is obvious, the time scale of current transformer is smaller,
The time scale of transmission system is larger, since double-fed blower unstable state appears in transmission system and induction machine part, because
This, the application has ignored the dynamic process of current transformer and its controller in double-fed blower, establishes comprising legacy system, induced electricity
The Small signal stability analysis model of machine and filter.
Transmission system uses simple substance block models, and dynamical equation is as follows:
In formula, ωrFor rotor speed, Tm,TeRespectively driving torque, electromagnetic torque, HgFor rotary inertia.
Induction machine uses influence generator model, and introduces transient potential, and dynamical equation is as follows:
In formula, ids,iqsFor stator current, idr,iqrFor rotor current, v'd,v'qFor induction machine transient voltage, vdr,vqr
For rotor voltage, Lss,LrrRespectively stator and rotor reactance, ωs,ωBThe respectively per unit value and a reference value of synchronous speed, Xs,Xs' point
Not Wei stator reactance and transient state reactance, RsFor stator resistance, Lm,Lrr,LssRespectively magnetizing inductance, inductor rotor and stator electricity
Sense.
Inverter outlet filter model:
In formula, idg,iqgFor filter current, Rg,LgRespectively filter resistance and inductance, vds,vqsFor stator voltage.
Stator current and inverter current are finally injected into power grid by double-fed blower.Formula (1)-(7) are attached in equalization point
It is close to carry out linearization process, state-space expression of the double-fed blower under a certain operating status can be obtained:
In formula, xDFIGFor the quantity of state of double-fed blower, uDFIGFor the control amount of double-fed blower, yDFIGFor the defeated of double-fed blower
Output, ADFIGFor the state matrix of double-fed blower, element therein by blower under the parameter and the operating status of blower state
Amount is constituted.
The establishment process of photo-voltaic power supply Small signal stability analysis model is as follows:
Compared with double-fed blower, the photovoltaic plant based on inverter have lesser time constant, dynamic process variation compared with
Fastly, therefore the dynamic process of its controller should be also taken into account.The application regards photovoltaic plant as constant voltage source, passes through line
Road is connected on bulk power grid, and structural schematic diagram is as shown in Fig. 2, controller schematic diagram is as shown in Figure 3.In figure, Lf,CfRespectively
Filter inductance, capacitor, LcFor line inductance, vbFor bulk power grid.Voltage controller: its input value is filter exit potential vod,
voq, electric current iod,ioq, output valveAs the target value of current controller, dynamical equation is as follows:
In formula, xvd,xvqFor voltage controller state variable.
Current controller: the target value and inverter outlet electric current i that voltage controller is given are receivedld,ilq, export inversion
The control voltage value of deviceDynamical equation is as follows:
Filter, capacitor, circuit model are as follows:
In formula, ω is mains frequency, vid,viqFor inverter outlet voltage.
In the current implantation systems that photovoltaic plant finally exports filter, therefore its output quantity is filter outlet electricity
Stream.Formula (9)-(20) are subjected to linearization process, state-space expression of the photovoltaic plant under a certain operating status can be obtained:
In formula, xpvFor the quantity of state of photovoltaic, upvFor the control amount of photovoltaic, ypvFor the output quantity of photovoltaic, ApvFor the shape of photovoltaic
State matrix, element therein are made of the quantity of state of photovoltaic under the parameter of photovoltaic and the operating status.
The establishment process of natural gas power supply Small signal stability analysis model is as follows:
In integrated energy system, synchronous generator connects the knob of electrical system as heating system with natural gas system
Band.Synchronous machine uses traditional third-order model, as follows:
Ud=XqIq (22)
Uq=E'q-X'dId (23)
P δ=ω -1 (24)
In formula, δ is phase angle difference, and ω is rotor speed, TJFor rotary inertia, Tm,TeRespectively machine torque, electromagnetic torque,
Ef,E'qRespectively excitation voltage, motor q axis transient potential, Xd,X'dThe respectively reactance of d axis and transient state reactance, idFor d axis stator
Electric current, T 'd0For d axis open circuit time constant.
Formula (22)-(24) are subjected to linearization process, state space table of the photovoltaic plant under a certain operating status can be obtained
Up to formula:
In formula, xsynFor the quantity of state of synchronous machine, usynFor the control amount of synchronous machine, AsynFor the state matrix of synchronous machine,
In element be made of the quantity of state of synchronous machine under the parameter of synchronous machine and the operating status.
By coordinate transform, the dq coordinate system of each model is transformed into common coordinate system DQ coordinate system, it is whole to carry out model
It closes, detailed process is as follows:
The electric current that double-fed blower is injected into system is the sum of stator current and net side inverter outlet electric current, i.e. (ids+
idg)+j(iqs+iqg)；The electric current that photovoltaic plant is injected to system is that filter exports electric current, i.e. iod+jioq；Synchronous machine is to system
The electric current of middle injection is stator current, i.e. id+jiq.It is worth noting that, above-mentioned model is all built upon one machine infinity bus system
Under, these Injection Currents are built upon entirely under respective dq coordinate system.If it is intended to the dynamic characteristic to whole system is ground
Study carefully, need to be coordinately transformed, respective dq coordinate system is transformed into common coordinate system DQ coordinate system, coordinate transform is as follows:
Wherein, DQ coordinate system represents common coordinate system, diqiCoordinate system represents the coordinate system of each power supply, δiIt is two
Angle between a coordinate system.
For photovoltaic plant, there is no the rotors of rotation, therefore the not equilibrium equation of torque；For synchronous machine and double-fed
Blower can obtain the electromagnetic power increment of unit if ignoring the internal resistance of stator winding are as follows:
ΔPe=Ux0ΔIx+Uy0ΔIy+Ix0ΔUx+Iy0ΔUy (29)
In formula, Δ Ix, Δ IyIt is unit to network Injection Current increment, Δ Ux,ΔUyFor grid entry point voltage increment.In the following,
Calculate separately the voltage and current increment of synchronous machine, double-fed blower, photovoltaic.Synchronous machine can be obtained by formula (22), (23), (28)
Double-fed blower can be obtained by formula (8), (28)
Photovoltaic can be obtained by formula (21), (28)
Meanwhile eliminating network equationIn contact node, and substitute into formula (30)-(32) and can obtain,
The increment expression formula of electromagnetic torque, and convolution (8), (21), (27) can be obtained by bringing formula (30)-(33) into formula (29)
The state-space expression of the integrated energy system can be obtained:
In formula, the state matrix A of system is made of network, double-fed blower, photovoltaic, synchronous machine and their operating status,
Model is complex.
Using loss minimization as objective function, and corresponding constraint condition is introduced, establishes the reconstruction model for considering small interference,
Establishment process is as follows:
For Distribution system generally to reduce network loss, raising network voltage level etc. as objective function, constraint condition includes trend
Constraint, the radial constraint of network, Branch Power Flow constraint and active and reactive units limits.After network structure changes, respectively
The small signal stability of class power supply operation is necessarily influenced to different extents.If Distribution system only considers economy, then very
Be likely to cause certain power supplys can not stable operation the case where, lead to off-grid.In addition to this, wind-powered electricity generation, photovoltaic power output have
The uncertainty of randomness, power output inevitably results in the variation of trend and node voltage in system, to influence unit
Electromagnetic power, so that unit is constantly subjected to the test of small interference stability in distribution.
Power distribution network reconfiguration should ensure that the small signal stability of the system operation before reconstruct next time.Go out in view of wind-powered electricity generation, photovoltaic
The uncertainty of power needs to guarantee in the worst cases to be still small interference stability.It, can be with simply for blower and synchronous machine
The power characteristic of electric system analyzes its analysis for running most unstable scene, as shown in Figure 5.
Wherein, PmFor mechanical output, PeFor electromagnetic power, there are two intersection points a, b for two curves, when unit operates in point a
When, by can be settled out after small interference；Conversely, when unit operates in point b, by being unstable after small interference.Cause
This, when unit mechanical output improves, and electromagnetic power reduces, the operating status of unit is very easy to be transitioned into point b from point a, i.e.,
It is unstable to be changed into small interference from small interference stability.Therefore, it is larger to appear in unit output to scene most unstable in system, but
In the lesser situation of network electromagnetic torque.
The application is proposed to constrain small interference stability and is added in the constraint of Distribution system, while by the power output of wind-powered electricity generation, photovoltaic
System is arranged in run under most unstable scene, so that the result of Distribution system can guarantee system small interference stability
On the basis of, it is more economical to run.It models as follows:
1) using loss minimization as objective function:
2) constraint condition are as follows:
The radial constraint of network
Power-balance constraint
Voltage constraint
Vkmin≤Vk≤Vkmax (39)
Branch Power Flow constraint:
Pijmin≤Pij×sij≤Pijmax (40)
Qijmin≤Qij×sij≤Qijmax (41)
Small interference stability constraint: the state that can obtain the system by the state-space expression (34) of integrated energy system is empty
Between matrix A, the characteristic value of calculating matrix A can obtain the system by the Oscillatory mode shape after small interference, in all characteristic values, most
It is minimum close to the corresponding Oscillatory mode shape stability margin of the characteristic value of the imaginary axis, and it is most slow to decay, therefore the Oscillatory mode shape is referred to as master
Oscillatory mode shape is led, characteristic value is also referred to as dominant characteristics value, and the application is denoted as λ1.For the small signal stability for guaranteeing system, need
Will the dominant characteristics value to system use restraint, be set as λref, then small interference stability is constrained to
λ1≤λref (42)
In above-mentioned model, how small interference stability is constrained into linearisation and expressed, become the work of a great challenge.
According to the calculation formula of characteristic value
| A- λ E |=0 (43)
It is found that characteristic value is obtained by solving Equation of Higher Degree with One Unknown, it is difficult directly to acquire the analytical expression of characteristic value,
It can only seek its numerical solution.In addition to this, which belongs to nonlinearity equation, and includes the operation of system
State, non trivial solution can with the variation of system running state rather than linearly change.
For these reasons, present applicant proposes the method for solving of small interference stability verification, this method is intended to provide
The fast solution method of one consideration small interference stability Distribution system, it cannot be guaranteed that method eventually converges to globally optimal solution.
(1) operating status of each power supply under initial network structure is calculated, and calculates the characteristic value of system, feature under this condition
Vector, and small interference stability is constrained linearize accordingly；
(2) the small interference stability constraint after linearisation is added in model, solves Distribution system scheme；
(3) operating status of each power supply under the reconfiguration scheme, and the dominant characteristics value of computing system again are calculated；
(4) judge whether the dominant characteristics value of system meets small interference stability constraint.If it is satisfied, then output reconstruct side
Case, if conditions are not met, then using the reconfiguration scheme as primary condition, return step (1) is recalculated.
Although the program is calculated using original state always in an iterative process, again to small dry after reconstruct
It disturbs stability to be verified, can guarantee that result calculated meets small interference stability constraint.
The above described is only a preferred embodiment of the present invention, be not intended to limit the present invention in any form, it is all
It is any simple modification, equivalent change and modification to the above embodiments according to the technical essence of the invention, still falls within
In the range of technical solution of the present invention.
Claims (7)
1. a kind of Distribution system method for considering integrated energy system small interference stability, to various types of in integrated energy system
Power supply establishes dynamic model, and small interference stability is added in the constraint of distribution reconstruct, it is ensured that integrated energy system meet it is small
It is more economical to run on the basis of interference stability, it the steps include: Step 1: with the wind-powered electricity generation electricity based on double-fed wind generating technology
Source, the photo-voltaic power supply based on inverter generation technology, the natural gas power supply based on gas turbine are established respective respectively as object
Small interference stability model；
Step 2: the dq coordinate system of each model is transformed into common coordinate system DQ coordinate system by coordinate transform, model is carried out
Integration；
Step 3: calculating the operating status of each power supply under initial network structure, and calculate the characteristic value of system, feature under this condition
Small interference stability is constrained and is linearized by vector；
Step 4: using loss minimization as objective function, and corresponding constraint condition is introduced, establish the reconstruct mould for considering small interference
Type；
Step 5: the small interference stability constraint after linearisation is added in reconstruction model, Distribution system scheme is solved；
Step 6: calculating the operating status of each power supply under the reconfiguration scheme, and the dominant characteristics value of computing system again；
Step 7: judging whether the dominant characteristics value of system meets small interference stability constraint；Meet, then exports reconfiguration scheme, no
Meet, then using the reconfiguration scheme as primary condition, then return step three.
2. the Distribution system method according to claim 1 for considering integrated energy system small interference stability, it is characterized in that: institute
It states in step 1, wind-powered electricity generation power supply Small signal stability analysis model is established with legacy system, induction machine and filter, and carry out
Linearization process can obtain the state-space expression of wind-powered electricity generation power supply:
In formula, xDFIGFor the quantity of state of double-fed blower, uDFIGFor the control amount of double-fed blower, yDFIGFor the output quantity of double-fed blower,
ADFIGFor the state matrix of double-fed blower, element therein by blower under the parameter and the operating status of blower quantity of state institute structure
At.
3. the Distribution system method according to claim 1 for considering integrated energy system small interference stability, it is characterized in that: institute
It states in step 1, photo-voltaic power supply small interference stability point is established with voltage controller, current controller, filter, capacitor, route
Model is analysed, and carries out linearization process, the state-space expression of photo-voltaic power supply can be obtained:
In formula, xpvFor the quantity of state of photovoltaic, upvFor the control amount of photovoltaic, ypvFor the output quantity of photovoltaic, ApvFor the state square of photovoltaic
Battle array, element therein are made of the quantity of state of photovoltaic under the parameter of photovoltaic and the operating status.
4. the Distribution system method according to claim 1 for considering integrated energy system small interference stability, it is characterized in that: institute
It states in step 1, synchronous generator connects the tie of electrical system as heating system with natural gas system, with three traditional ranks
Model carries out linearization process as natural gas power supply Small signal stability analysis model, can obtain the state space of photovoltaic plant
Expression formula:
In formula, xsynFor the quantity of state of synchronous machine, usynFor the control amount of synchronous machine, AsynIt is therein for the state matrix of synchronous machine
Element is made of the quantity of state of synchronous machine under the parameter of synchronous machine and the operating status.
5. the Distribution system method according to claim 1 for considering integrated energy system small interference stability, it is characterized in that: institute
It states in step 2, the electric current that fan power is injected into system is the sum of stator current and net side inverter outlet electric current；Photovoltaic
The electric current that power supply is injected to system is that filter exports electric current；The electric current that natural gas power supply is injected into system is stator current；
Respective dq coordinate system is transformed into common coordinate system DQ coordinate system according to coordinate transform formula, coordinate transform is as follows:
Wherein, DQ coordinate system represents common coordinate system, diqiCoordinate system represents the coordinate system of each power supply, δiIt is sat for two
Angle between mark system；
According to the voltage and current increment of photovoltaic in double-fed blower in synchronous machine in natural gas power supply, wind-powered electricity generation power supply, photo-voltaic power supply, and
The state-space expression of convolution wind-powered electricity generation power supply, the state-space expression of photo-voltaic power supply, photovoltaic plant state space table
The state-space expression of the integrated energy system can be obtained up to formula:
In formula, the state matrix A of system is made of network, double-fed blower, photovoltaic, synchronous machine and their operating status.
6. the Distribution system method according to claim 1 for considering integrated energy system small interference stability, it is characterized in that: institute
It states in step 3, the state space matrices A of the system, calculating matrix can be obtained by the state-space expression of integrated energy system
The characteristic value of A can obtain the system by the Oscillatory mode shape after small interference；
It is leading Oscillatory mode shape with Oscillatory mode shape corresponding to the characteristic value near the imaginary axis, characteristic value is leading eigenvalue λ1；It is main
That leads characteristic value is constrained to λref, then small interference stability constrains are as follows:
λ1≤λref (42)
According to the calculation formula of characteristic value:
| A- λ E |=0 (43)
Small interference stability is constrained and is linearized.
7. the Distribution system method according to claim 1 for considering integrated energy system small interference stability, it is characterized in that: institute
It states in step 4, with the objective function of loss minimization are as follows:
The constraint condition includes the radial constraint of network, power-balance constraint, voltage constraint, Branch Power Flow constraint and small interference
Scleronomic constraint.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910514943.9A CN110336299B (en) | 2019-06-14 | 2019-06-14 | Distribution network reconstruction method considering small interference stability of comprehensive energy system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910514943.9A CN110336299B (en) | 2019-06-14 | 2019-06-14 | Distribution network reconstruction method considering small interference stability of comprehensive energy system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110336299A true CN110336299A (en) | 2019-10-15 |
| CN110336299B CN110336299B (en) | 2023-12-15 |
Family
ID=68142045
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910514943.9A Active CN110336299B (en) | 2019-06-14 | 2019-06-14 | Distribution network reconstruction method considering small interference stability of comprehensive energy system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110336299B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113300704A (en) * | 2021-06-01 | 2021-08-24 | 上海交通大学 | Parameter design method and system for phase-locked loop of complex filter |
| CN114285037A (en) * | 2022-03-07 | 2022-04-05 | 天津大学 | Method for determining control parameter stability region of regional electricity-gas integrated energy system |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103050992A (en) * | 2012-12-10 | 2013-04-17 | 杜文娟 | Analyzing method for damping torque having small interfering stable influence on power system caused by wind power integration |
| CN103904649A (en) * | 2014-03-28 | 2014-07-02 | 南方电网科学研究院有限责任公司 | Aggregation method suitable for generators different in order model |
| CN104933639A (en) * | 2015-07-02 | 2015-09-23 | 浙江大学 | A small-interference stability rapid analysis method targeted at a large scale electric power system |
| US20150318697A1 (en) * | 2013-03-29 | 2015-11-05 | Gansu Electric Power Corporation Wind Power Technology Center | A method for improving small disturbance stability after double-fed unit gets access to the system |
| CN105914783A (en) * | 2016-05-10 | 2016-08-31 | 西安交通大学 | Island type microgrid small interference stability analysis method |
| CN106203883A (en) * | 2016-07-21 | 2016-12-07 | 中国电力科学研究院 | A kind of based on interval power distribution network small-signal stability computational methods |
| CN107834564A (en) * | 2017-11-02 | 2018-03-23 | 上海电力学院 | A kind of control method of the Small disturbance voltage stability of micro-grid system |
| CN108063461A (en) * | 2018-01-04 | 2018-05-22 | 华北电力大学（保定） | A kind of active power dispatch method of electric system containing wind-powered electricity generation for considering small interference stability risk |
| CN108565863A (en) * | 2018-04-13 | 2018-09-21 | 国网浙江省电力有限公司电力科学研究院 | A kind of regional complex energy resource system multiple target tide optimization method considering randomness |
-
2019
- 2019-06-14 CN CN201910514943.9A patent/CN110336299B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103050992A (en) * | 2012-12-10 | 2013-04-17 | 杜文娟 | Analyzing method for damping torque having small interfering stable influence on power system caused by wind power integration |
| US20150318697A1 (en) * | 2013-03-29 | 2015-11-05 | Gansu Electric Power Corporation Wind Power Technology Center | A method for improving small disturbance stability after double-fed unit gets access to the system |
| CN103904649A (en) * | 2014-03-28 | 2014-07-02 | 南方电网科学研究院有限责任公司 | Aggregation method suitable for generators different in order model |
| CN104933639A (en) * | 2015-07-02 | 2015-09-23 | 浙江大学 | A small-interference stability rapid analysis method targeted at a large scale electric power system |
| CN105914783A (en) * | 2016-05-10 | 2016-08-31 | 西安交通大学 | Island type microgrid small interference stability analysis method |
| CN106203883A (en) * | 2016-07-21 | 2016-12-07 | 中国电力科学研究院 | A kind of based on interval power distribution network small-signal stability computational methods |
| CN107834564A (en) * | 2017-11-02 | 2018-03-23 | 上海电力学院 | A kind of control method of the Small disturbance voltage stability of micro-grid system |
| CN108063461A (en) * | 2018-01-04 | 2018-05-22 | 华北电力大学（保定） | A kind of active power dispatch method of electric system containing wind-powered electricity generation for considering small interference stability risk |
| CN108565863A (en) * | 2018-04-13 | 2018-09-21 | 国网浙江省电力有限公司电力科学研究院 | A kind of regional complex energy resource system multiple target tide optimization method considering randomness |
Non-Patent Citations (1)
| Title |
|---|
| 邓年俊: "光伏并网发电系统小信号建模与分析", 《电气开关》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113300704A (en) * | 2021-06-01 | 2021-08-24 | 上海交通大学 | Parameter design method and system for phase-locked loop of complex filter |
| CN114285037A (en) * | 2022-03-07 | 2022-04-05 | 天津大学 | Method for determining control parameter stability region of regional electricity-gas integrated energy system |
| CN114285037B (en) * | 2022-03-07 | 2022-05-27 | 天津大学 | Method for determining control parameter stability region of regional electricity-gas integrated energy system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110336299B (en) | 2023-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xiong et al. | Backstepping-based DPC strategy of a wind turbine-driven DFIG under normal and harmonic grid voltage | |
| CN102856917B (en) | Reactive power optimization method of power distribution network | |
| Fathabadi | Novel high-efficient unified maximum power point tracking controller for hybrid fuel cell/wind systems | |
| CN105429170B (en) | A kind of microgrid inverter droop control method based on adjustable virtual impedance | |
| CN106374528A (en) | Microgrid dispersing type control policy model based on improved droop control | |
| CN103414212B (en) | Consider the analogy method of the distribution network systems of double-fed wind generator machine access | |
| Bezza et al. | Sensorless MPPT fuzzy controller for DFIG wind turbine | |
| CN105680479A (en) | Virtual synchronous generator control method and system by considering dynamic characteristics of photovoltaic power supply | |
| Bouzelata et al. | Exploration of optimal design and performance of a hybrid wind-solar energy system | |
| CN105608244B (en) | Distributed photovoltaic power generation system based on PSCAD/EMTDC | |
| CN108400593A (en) | Active distribution network electrical model method for building up based on layering multi-agent technology | |
| Shao et al. | Adequacy of the single-generator equivalent model for stability analysis in wind farms with VSC-HVDC systems | |
| CN110206686A (en) | A kind of adaptive maximum power tracking and controlling method for wind power generating set | |
| CN106451539B (en) | It is a kind of meter and permanent magnet direct-drive wind turbine group dynamic characteristic wind farm grid-connected Method of Stability Analysis | |
| CN110336299A (en) | A kind of Distribution system method considering integrated energy system small interference stability | |
| CN106294959B (en) | The modeling and simulating method that model reference self-adapting control is combined with double-fed blower reduced-order model | |
| Zhang et al. | Enhanced wind turbine maximum wind-energy capture based on the inverse-system method | |
| CN106712032A (en) | Optimal power flow model construction method considering active power voltage regulation capacity of wind turbine generator set | |
| CN104795842B (en) | Symmetrical fault ride-through control method for hybrid wind farm group containing double-fed wind farm and permanent magnetic direct drive wind farm | |
| CN106786766B (en) | A method of the raising wind-powered electricity generation maximum grid connection capacity based on P2G technology | |
| Coath et al. | Effect of steady-state wind turbine generator models on power flow convergence and voltage stability limit | |
| Eminoglu et al. | Incorporation of a new wind turbine generating system model into distribution systems load flow analysis | |
| Li et al. | Structure preserving aggregation method for doubly-fed induction generators in wind power conversion | |
| Xuefeng et al. | A brief analysis of distributed generation connected to distribution network | |
| CN109546687A (en) | Using the new-energy grid-connected method for improving droop control technology based on virtual power |
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 |