CN105281360B - A kind of distributed photovoltaic automatic power generation control method based on sensitivity - Google Patents

A kind of distributed photovoltaic automatic power generation control method based on sensitivity Download PDF

Info

Publication number
CN105281360B
CN105281360B CN201510583592.9A CN201510583592A CN105281360B CN 105281360 B CN105281360 B CN 105281360B CN 201510583592 A CN201510583592 A CN 201510583592A CN 105281360 B CN105281360 B CN 105281360B
Authority
CN
China
Prior art keywords
mrow
mtd
node
mtr
power
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510583592.9A
Other languages
Chinese (zh)
Other versions
CN105281360A (en
Inventor
韩志军
王坚敏
严耀良
朱伟
陈国恩
王跃强
张磊
任志翔
罗里志
殷华
施晓光
张仲孝
王波
王波一
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Jiaxing Power Supply Co of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Jiaxing Power Supply Co of State Grid Zhejiang Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Jiaxing Power Supply Co of State Grid Zhejiang Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510583592.9A priority Critical patent/CN105281360B/en
Publication of CN105281360A publication Critical patent/CN105281360A/en
Application granted granted Critical
Publication of CN105281360B publication Critical patent/CN105281360B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Abstract

A kind of distributed photovoltaic automatic power generation control method based on sensitivity, the control method include the following steps:a)By to the local linearization near distribution steady-state operation operating point, the mathematical relationship that power distribution network node voltage disturbance quantity injects active reactive disturbance amount with node is solved, i.e., variation relation is injected by Sensitivity Analysis Method analysis node voltage magnitude and node active reactive;b)The sensitivity relation injected according to node voltage and node, choose respectively distribution interior joint it is active idle injection is out-of-limit to node voltage or the node of voltage sensibility maximum of the violent node of fluctuation be it is active Reactive-power control node, the regulation and control of You Gong idle work optimizations are participated in, the photovoltaic generating system for being not involved in adjusting always is operate on maximal power tracing state;c)With in distribution distributed photovoltaic contribute maximize and node voltage offrating be minimised as target, using it is active Reactive-power control node it is active idle power generating value for solution variable, establish distributed photovoltaic Optimized model.

Description

A kind of distributed photovoltaic automatic power generation control method based on sensitivity
Technical field
The present invention relates to a kind of distributed photovoltaic automatic power generation control method based on sensitivity, belong to photovoltaic generation Technical field.
Background technology
Increasingly serious with energy crisis, regenerative resource is greatly developed and becomes trend, and photovoltaic generation is it In a kind of important form.With the development of distributed photovoltaic power generation technology, energy policy and the further of electricity market open, The operation of power distribution network needs to meet the high degree of compatibility to distribution type renewable energy power generation.The distribution of access power distribution network at present Photovoltaic generating system is usually to run on maximal power tracing (MPPT) pattern as a uncontrollable power supply, but works as power distribution network When middle distributed photovoltaic permeability is continuously improved, due to the intrinsic randomness of photovoltaic generation and intermittence so that distribution power flow With supply voltage frequent fluctuation, power quality deteriorates, and power supply reliability reduces, and abandons optical phenomenon getting worse, causes serious money Source wastes, and lot of challenges is brought to power distribution network safety, economy, stable operation.Photovoltaic DC-to-AC converter has considerable active nothing in itself Work(regulating power, vehicle air-conditioning is carried out by being exported to photovoltaic DC-to-AC converter active reactive, and grid-connected band is eliminated from source The uncertainty come, can improve the access level of distributed photovoltaic in distribution, improve resource utilization, realize regenerative resource Close friend's access, has broad prospects and huge economic implications.
The voltage fluctuation problem brought currently for distributed photovoltaic access power distribution network, mainly using adjusting substation Load tap changer, switching shunt capacitor and the installation Continuous Var Compensation device (SVC) in power distribution network, are adjusted.But It is load tap changer, and shunt capacitor is frequently adjusted can substantially reduce service life of equipment, and be difficult to power distribution network Voltage realizes smooth adjustment, easily causes to impact to power distribution network, and the reactive-load compensation equipment installation cost such as SVC is higher, is not suitable for big rule Mould application.And when traditional Reactive-power control means in distribution cannot effectively suppress voltage fluctuation, and it be easy to cause extensive Optical phenomenon is abandoned, easily causes the waste of resource.For such a situation, research it is a kind of based on photovoltaic DC-to-AC converter itself it is active-nothing The Optimal Control Strategy of work(regulating power has huge economic value and broad prospect of application.
The content of the invention
It is an object of the invention to overcome the above-mentioned problems of the prior art, there is provided a kind of distribution based on sensitivity Photovoltaic automatic power generation control method;It derives node voltage amplitude disturbance and active-idle injection of node by sensitivity analysis The relation of disturbance, then choose influences maximum node active injection to the voltage magnitude of voltage out-of-limit node is injected to idle Regulated variable, photovoltaic active power output is maximum in a distributed manner is distributed with the minimum target of node voltage amplitude offrating, foundation Formula photovoltaic coordinates control Optimized model, then by the PSO Algorithm model, tries to achieve the active adjusting for participating in adjusting respectively The active output of node and the idle output of Reactive-power control node, ensure to maintain while distributed photovoltaic power generation is fully dissolved Power distribution network node voltage level is within normal operation range.
The technical solution adopted by the present invention is:A kind of distributed photovoltaic automatic power generation control method based on sensitivity, its It is characterized in that the control method includes the following steps:
A) by the local linearization near distribution steady-state operation operating point, solve power distribution network node voltage disturbance quantity with Node injects the mathematical relationship of active-reactive disturbance amount, i.e., has by Sensitivity Analysis Method analysis node voltage magnitude with node Work(is idle injection variation relation;
B) sensitivity relation injected according to node voltage and node, choose respectively distribution interior joint it is active idle injection The node of out-of-limit to node voltage or the violent node of fluctuation voltage sensibility maximum be it is active Reactive-power control node, participated in Gong idle work optimizations regulate and control, and the photovoltaic generating system for being not involved in adjusting always is operate on maximal power tracing (MPPT) state;
C) with distributed photovoltaic output maximizes in distribution and node voltage offrating is minimised as target, to have Work(Reactive-power control node it is active idle power generating value to solve variable, establish distributed photovoltaic Optimized model, and using a kind of Immune bifurcation PSO Algorithm is under conditions of meeting with network operation constraint constraint, active-idle output of adjustment node Optimal value.
Preferably, in step a) of the present invention, it is described by the office near distribution steady-state operation operating point Portion linearizes, and solves the mathematical relationship that power distribution network node voltage disturbance quantity injects active-reactive disturbance amount with node, specifically such as Under:
For the electric power networks with N number of node, n=N-1 is made, takes balance nodes to increase PV node as node is referred to Extensively into flow equation, the Newton Power Flow update equation formula of 2n polar form can be obtained:
Wherein:For Jacobian matrix, Δ θ and Δ V are respectively node voltage phase angle and amplitude disturbance vector.Then:
Wherein:V ties up node voltage amplitude diagonal matrix, i.e. V=diag (v for n1,v2,...,vn).B and G is respectively that node is led Receive the real and imaginary parts of battle array, Bcos θ are a kind of simplified literary style of matrix, it and B have identical structure, its each several part element is Corresponding element B in BijWith cos θijProduct, other are similar, in addition, above-mentioned P and Q are n rank diagonal matrixs, its diagonal element Respectively Pi/Vi 2And Qi/Vi 2.Under normal circumstances, θijIt is very small, therefore cos θ can be madeij=1, sin θij=0, above formula can be done into one Step abbreviation obtains:
In view of normal operation lower node voltage magnitude perunit value near 1.0pu, above formula Blocked portion is carried out Gaussian elimination, sensitivity of the voltage to node injecting power, which is calculated, is:
Δ V=((B+Q) (G-P)-1(B-Q)+(G+P))-1ΔP-((G-P)(B+Q)-1(G+P)+(B-Q))-1ΔQ
As preferred:In step c) of the present invention, contributed and maximized with distributed photovoltaic in distribution, distribution stable state fortune Row active power loss minimize and node voltage offrating be minimised as target, with it is active Reactive-power control node it is active without Work(power generating value is solution variable, establishes distributed photovoltaic Optimized model, specific as follows:
(1) object function:
In formula:DGiFor photovoltaic access node in distribution,By the photovoltaic active power output predicted value mounted in node i, Pi,DGFor distributed photovoltaic active power output control targe value in node i;vjFor node j voltage magnitudes,For node j voltage magnitudes Rated value (being perunit value), α are weight coefficient and 0 < α < 1;
(2) distribution power-balance constraint:
In formula:PiWith QiThe respectively active reactive injection of node i,For node i voltage phasor, YijFor system admittance square Battle array corresponding element;
(3) node voltage constrains:
In formula:viFor node i voltage magnitude,WithThe maximal and minmal value that respectively node voltage amplitude allows.
(4) critical point exchanges power constraint:
To suppress influence of the photovoltaic generation power fluctuation to higher level's power grid, it need to consider that power distribution network root node exchanges power limit System, i.e.,:
In formula:P1And Q1The active reactive power of power distribution network is respectively flowed into from root node;WithRespectively dispatch Centrally disposed critical point is active to exchange power bound,WithThe critical point reactive power exchange work(that respectively control centre is set Rate bound;
(5) distributed photovoltaic module operation constraint:
PQ models are used during photovoltaic module steady-state operation, by grid-connected inverters, grid-connected power can realize that active reactive is only Vertical control;In formula, Pi,DGFor node i photovoltaic module active power output value,For node i photovoltaic active power output predicted value, Qi,DGFor Node i photovoltaic module is idle power generating value, Si,DGFor node i photovoltaic module capacity.
The present invention has the following advantages compared with prior art:
The fact that do not decouple is injected with idle in view of power distribution network interior joint is active, comprehensively utilizes distributed photovoltaic power generation system Unite active-Reactive-power control ability, contribute than single adjusting photovoltaic or the idle delivery efficiency higher of reactive-load compensation equipment;Utilize The characteristic that photovoltaic DC-to-AC converter obtains in itself suppresses the fluctuation of power distribution network node voltage, it is not necessary to extra investment, is adapted to big Sizable application;Method based on sensitivity analysis chooses active-Reactive-power control node to key node voltage influence maximum, often Secondary adjusting process avoids the adjusting to all photovoltaic modules, while the solution variable in an optimization cycle will not be with distribution The increase of middle photovoltaic access point improves, and is not in that model is difficult to when distributed photovoltaic permeability gradually steps up in power distribution network Solve problems, are particularly suitable for distributed photovoltaic and access situation on a large scale.
Brief description of the drawings
Fig. 1 is the method for the present invention flow chart.
Fig. 2 is utilizes the flow chart of PSO Algorithm institute's established model of the present invention.
Embodiment
The embodiment of the present invention is further described below in conjunction with the accompanying drawings.A kind of point based on sensitivity Cloth photovoltaic automatic power generation control method, the control method include the following steps:
A) by the local linearization near distribution steady-state operation operating point, solve power distribution network node voltage disturbance quantity with Node injects the mathematical relationship of active-reactive disturbance amount, i.e., has by Sensitivity Analysis Method analysis node voltage magnitude with node Work(is idle injection variation relation;
B) sensitivity relation injected according to node voltage and node, choose respectively distribution interior joint it is active idle injection The node of out-of-limit to node voltage or the violent node of fluctuation voltage sensibility maximum be it is active Reactive-power control node, participated in Gong idle work optimizations regulate and control, and the photovoltaic generating system for being not involved in adjusting always is operate on maximal power tracing (MPPT) state;
C) with distributed photovoltaic output maximizes in distribution and node voltage offrating is minimised as target, to have Work(Reactive-power control node it is active idle power generating value to solve variable, establish distributed photovoltaic Optimized model, and using a kind of Immune bifurcation PSO Algorithm is under conditions of meeting with network operation constraint constraint, active-idle output of adjustment node Optimal value.
It is described by matching somebody with somebody to the local linearization near distribution steady-state operation operating point, solution in the step a) Grid nodes voltage disturbance amount injects the mathematical relationship of active-reactive disturbance amount with node, specific as follows:
For the electric power networks with N number of node, n=N-1 is made, takes balance nodes to increase PV node as node is referred to Extensively into flow equation, the Newton Power Flow update equation formula of 2n polar form can be obtained:
Wherein:For Jacobian matrix, Δ θ and Δ V are respectively node voltage phase angle and amplitude disturbance vector.Then:
Wherein:V ties up node voltage amplitude diagonal matrix, i.e. V=diag (v for n1,v2,...,vn).B and G is respectively that node is led Receive the real and imaginary parts of battle array, Bcos θ are a kind of simplified literary style of matrix, it and B have identical structure, its each several part element is Corresponding element B in BijWith cos θijProduct, other are similar, in addition, above-mentioned P and Q are n rank diagonal matrixs, its diagonal element Respectively Pi/Vi 2And Qi/Vi 2.Under normal circumstances, θijIt is very small, therefore cos θ can be madeij=1, sin θij=0, above formula can be done into one Step abbreviation obtains:
In view of normal operation lower node voltage magnitude perunit value near 1.0pu, above formula Blocked portion is carried out Gaussian elimination, sensitivity of the voltage to node injecting power, which is calculated, is:
Δ V=((B+Q) (G-P)-1(B-Q)+(G+P))-1ΔP-((G-P)(B+Q)-1(G+P)+(B-Q))-1ΔQ
In the step c), contributed and maximized with distributed photovoltaic in distribution, distribution steady-state operation active power loss is minimum Change and node voltage offrating be minimised as target, using it is active Reactive-power control node it is active idle power generating value for solution Variable, establishes distributed photovoltaic Optimized model, specific as follows:
(1) object function:
In formula:DGiFor photovoltaic access node in distribution,By the photovoltaic active power output predicted value mounted in node i, Pi,DGFor distributed photovoltaic active power output control targe value in node i;vjFor node j voltage magnitudes,For node j voltage magnitudes Rated value (being perunit value), α are weight coefficient and 0 < α < 1;
(3) distribution power-balance constraint:
In formula:PiWith QiThe respectively active reactive injection of node i,For node i voltage phasor, YijFor system admittance square Battle array corresponding element;
(3) node voltage constrains:
In formula:viFor node i voltage magnitude,WithThe maximal and minmal value that respectively node voltage amplitude allows.
(4) critical point exchanges power constraint:
To suppress influence of the photovoltaic generation power fluctuation to higher level's power grid, it need to consider that power distribution network root node exchanges power limit System, i.e.,:
In formula:P1And Q1The active reactive power of power distribution network is respectively flowed into from root node;P1 minWith P1 maxRespectively dispatch Centrally disposed critical point is active to exchange power bound,WithThe critical point reactive power exchange work(that respectively control centre is set Rate bound;
(6) distributed photovoltaic module operation constraint:
PQ models are used during photovoltaic module steady-state operation, by grid-connected inverters, grid-connected power can realize that active reactive is only Vertical control;In formula, Pi,DGFor node i photovoltaic module active power output value,For node i photovoltaic active power output predicted value, Qi,DGFor Node i photovoltaic module is idle power generating value, Si,DGFor node i photovoltaic module capacity.
Embodiment:
1) contributed according to the photovoltaic of subsequent time and predicted load calculates next moment distribution trend and is distributed;The present invention adopts Trend is calculated with Newton-Laphson method, obtains the node voltage situation at distribution each node next moment;
2) required according to distribution network voltage amplitude, the node of next moment voltage out-of-limit is searched, if node number is k;
3) sensitivity of calculate node voltage injection active to node-idle, calculation formula are as follows:
Δ V=((B+Q) (G-P)-1(B-Q)+(G+P))-1ΔP-((G-P)(B+Q)-1(G+P)+(B-Q))-1ΔQ
B and G is respectively the real and imaginary parts of node admittance battle array in above formula, and P and Q are n rank diagonal matrixs, its diagonal element divides Wei not PiAnd Qi, for active-idle injection value of node i;Then:
Psense=((B+Q) (G-P)-1(B-Q)+(G+P))-1
Qsense=((G-P) (B+Q)-1(G+P)+(B-Q))-1
Psense、QsenseIt is n rank square formations, wherein i rows j column elements represent, and the active and reactive injection of j-th of node is to i-th The sensitivity of a node voltage amplitude, and if node does not access distributed photovoltaic, the active reactive for representing j nodes can not Adjust, then j column elements are assigned a value of zero;
4) this is searched to participate in adjusting active node and participate in adjusting idle node;
Result more than searches the node for participating in optimization, respectively in Psense、QsenseSearched in matrix in row k absolutely To the element that value is maximum, the corresponding node number difference i, j of maximum absolute value element.Then represent the active output of i-node, j nodes Influence of the idle output to node k voltages it is the most notable, then choose the photovoltaic active power output P of i-nodei,DG, the photovoltaic of j nodes Idle output Qi,DGDistributed photovoltaic Coordination and Optimization Model is established to solve variable;
5) active-idle Coordination and Optimization Model of distributed photovoltaic is as follows:
(1) object function:
In formula:DGiFor photovoltaic access node in distribution,By the photovoltaic active power output predicted value mounted in node i, Pi,DGFor distributed photovoltaic active power output control targe value in node i.vjFor node j voltage magnitudes,For node j voltage amplitudes Value rated value (being perunit value), α are weight coefficient and 0 < α < 1;
(2) distribution power-balance constraint
In formula:PiWith QiThe respectively active reactive injection of node i,For node i voltage phasor, YijFor system admittance square Battle array corresponding element;
(3) node voltage constrains
In formula:Vi is node i voltage magnitude,WithThe maximal and minmal value that respectively node voltage amplitude allows;
(4) critical point exchanges power constraint
To suppress influence of the photovoltaic generation power fluctuation to higher level's power grid, it need to consider that power distribution network root node exchanges power limit System, i.e.,:
In formula:P1And Q1The active reactive power of power distribution network is respectively flowed into from root node.P1 minWith P1 maxRespectively dispatch Centrally disposed critical point is active to exchange power bound,WithThe critical point reactive power exchange work(that respectively control centre is set Rate bound;
(5) distributed photovoltaic module operation constraint
PQ models are used during photovoltaic module steady-state operation, by grid-connected inverters, grid-connected power can realize that active reactive is only Vertical control;In formula, Pi,DGFor node i photovoltaic module active power output value,For node i photovoltaic active power output predicted value, Qi,DGFor Node i photovoltaic module is idle power generating value, Si,DGFor node i photovoltaic module capacity.
It is that subsequent time distributed photovoltaic is active and idle export is joined by the above-mentioned Optimized model of PSO Algorithm Examine value.
According to the magnitude of voltage of measurement system Real-time Feedback whether in defined power supply voltage range, if condition is met etc. Treat next period optimization, if still having voltage out-of-limit node, reduce object function weight coefficient α in Optimized model =0.5 α, again solving-optimizing model obtain distributed photovoltaic active reactive output reference value, until meet voltage magnitude will Ask.

Claims (2)

1. a kind of distributed photovoltaic automatic power generation control method based on sensitivity, it is characterised in that the control method includes as follows Step:
A) by the local linearization near distribution steady-state operation operating point, solving power distribution network node voltage disturbance quantity and node The mathematical relationship of active-reactive disturbance amount is injected, that is, passes through Sensitivity Analysis Method analysis node voltage magnitude and the active nothing of node Work(injects variation relation;
B) sensitivity relation injected according to node voltage and node, choose respectively distribution interior joint it is active idle injection to section The node of voltage sensibility maximum of point voltage out-of-limit or the violent node of fluctuation be it is active Reactive-power control node, participate in it is active without Work(Optimum Regulation, the photovoltaic generating system for being not involved in adjusting always are operate on maximal power tracing state;
C) with distribution distributed photovoltaic contribute maximize and node voltage offrating be minimised as target, with it is active without Work(adjustment node it is active idle power generating value to solve variable, establish distributed photovoltaic Optimized model, and using a kind of immune double State PSO Algorithm is under conditions of meeting with network operation constraint constraint, active-idle output optimal value of adjustment node;
It is described by the local linearization near distribution steady-state operation operating point, solving power distribution network section in the step a) Point voltage disturbance amount injects the mathematical relationship of active-reactive disturbance amount with node, specific as follows:
For the electric power networks with N number of node, n=N-1 is made, takes balance nodes to arrive PV node augmentation as node is referred to In flow equation, the Newton Power Flow update equation formula of 2n polar form can be obtained:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>H</mi> </mtd> <mtd> <mi>N</mi> </mtd> </mtr> <mtr> <mtd> <mi>M</mi> </mtd> <mtd> <mi>L</mi> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>&amp;Delta;</mi> <mi>&amp;theta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&amp;Delta;</mi> <mi>V</mi> <mo>/</mo> <mi>V</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>&amp;Delta;</mi> <mi>P</mi> </mtd> </mtr> <mtr> <mtd> <mi>&amp;Delta;</mi> <mi>Q</mi> </mtd> </mtr> </mtable> </mfenced> </mrow>
Wherein:For Jacobian matrix, Δ θ and Δ V are respectively node voltage phase angle and amplitude disturbance vector;Then:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>H</mi> </mtd> <mtd> <mi>N</mi> </mtd> </mtr> <mtr> <mtd> <mi>M</mi> </mtd> <mtd> <mi>L</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>V</mi> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mi>V</mi> </mtd> </mtr> </mtable> </mfenced> <mo>{</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>B</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>G</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>G</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> </mrow> </mtd> <mtd> <mrow> <mi>B</mi> <mi>cos</mi> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>G</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> </mtd> <mtd> <mrow> <mi>B</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>B</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> </mtd> <mtd> <mrow> <mi>G</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mi>Q</mi> </mrow> </mtd> <mtd> <mi>P</mi> </mtd> </mtr> <mtr> <mtd> <mi>P</mi> </mtd> <mtd> <mi>Q</mi> </mtd> </mtr> </mtable> </mfenced> <mo>}</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>V</mi> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mi>V</mi> </mtd> </mtr> </mtable> </mfenced> </mrow>
Wherein:V ties up node voltage amplitude diagonal matrix, i.e. V=diag (V for n1, V2..., Vn);B and G is respectively node admittance battle array Real and imaginary parts, Bcos θ are a kind of simplified literary style of matrix, it and B have identical structure, its each several part element is in B Corresponding element BijWith cos θijProduct, in addition, above-mentioned P and Q are n rank diagonal matrixs, its diagonal element is respectivelyWithUnder normal circumstances, θijIt is very small, therefore cos θ can be madeij=1, sin θij=0, above formula can do further abbreviation and obtain:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>B</mi> <mo>+</mo> <mi>Q</mi> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>G</mi> <mo>-</mo> <mi>P</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>G</mi> <mo>-</mo> <mi>P</mi> </mrow> </mtd> <mtd> <mrow> <mi>B</mi> <mo>-</mo> <mi>Q</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>V</mi> <mi>&amp;Delta;</mi> <mi>&amp;theta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&amp;Delta;</mi> <mi>V</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>&amp;Delta;</mi> <mi>P</mi> <mo>/</mo> <mi>V</mi> </mtd> </mtr> <mtr> <mtd> <mi>&amp;Delta;</mi> <mi>Q</mi> <mo>/</mo> <mi>V</mi> </mtd> </mtr> </mtable> </mfenced> </mrow>
In view of normal operation lower node voltage magnitude perunit value near 1.0pu, gaussian elimination is carried out to above formula, is calculated Obtaining sensitivity of the voltage to node injecting power is:
Δ V=((B+Q) (G-P)-1(B-Q)+(G+P))-1ΔP
-((G-P)(B+Q)-1(G+P)+(B-Q))-1ΔQ
2. the distributed photovoltaic automatic power generation control method according to claim 1 based on sensitivity, it is characterised in that institute In the step c) stated, contributed and maximized with distributed photovoltaic in distribution, distribution steady-state operation active power loss minimizes and node electricity Pressure offrating be minimised as target, using it is active Reactive-power control node it is active idle power generating value for solve variable, foundation Distributed photovoltaic Optimized model, it is specific as follows:
(1) object function:
<mrow> <mi>min</mi> <mi> </mi> <mi>f</mi> <mo>=</mo> <mi>&amp;alpha;</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <msub> <mi>DG</mi> <mn>1</mn> </msub> </mrow> <mrow> <msub> <mi>DG</mi> <mi>N</mi> </msub> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>D</mi> <mi>G</mi> </mrow> <mrow> <mi>Pr</mi> <mi>e</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>D</mi> <mi>G</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mi>j</mi> </msub> <mo>-</mo> <msubsup> <mi>v</mi> <mi>j</mi> <mrow> <mi>r</mi> <mi>a</mi> <mi>t</mi> <mi>e</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>
In formula:DGiFor photovoltaic access node in distribution,By the photovoltaic active power output predicted value mounted in node i, Pi,DG For distributed photovoltaic active power output control targe value in node i;VjFor node j voltage magnitudes, Vj rateFor node j voltage magnitude volumes Definite value, α are weight coefficient and 0 < α < 1;
(2) distribution power-balance constraint:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>Re</mi> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>Y</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>V</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>Im</mi> <mrow> <mo>(</mo> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>Y</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>V</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
In formula:PiWith QiThe respectively active reactive injection of node i,For node i voltage phasor, YijFor system admittance matrix pair Answer element;
(3) node voltage constrains:
<mrow> <msubsup> <mi>v</mi> <mi>i</mi> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>&amp;le;</mo> <msubsup> <mi>v</mi> <mi>i</mi> <mi>max</mi> </msubsup> </mrow>
In formula:ViFor node i voltage magnitude,WithThe maximal and minmal value that respectively node voltage amplitude allows;
(4) critical point exchanges power constraint:
To suppress influence of the photovoltaic generation power fluctuation to higher level's power grid, it need to consider that power distribution network root node exchanges power limit, i.e.,:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msubsup> <mi>P</mi> <mn>1</mn> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mn>1</mn> <mi>max</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>Q</mi> <mn>1</mn> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>Q</mi> <mn>1</mn> </msub> <mo>&amp;le;</mo> <msubsup> <mi>Q</mi> <mn>1</mn> <mi>max</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced>
In formula:P1And Q1The active reactive power of power distribution network is respectively flowed into from root node;P1 minWith P1 maxRespectively control centre The critical point of setting is active to exchange power bound, Q1 minWith Q1 maxOn the critical point reactive power exchange power that respectively control centre is set Lower bound;
(5) distributed photovoltaic module operation constraint:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>0</mn> <mo>&amp;le;</mo> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>D</mi> <mi>G</mi> </mrow> </msub> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>D</mi> <mi>G</mi> </mrow> <mrow> <mi>Pr</mi> <mi>e</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>P</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>D</mi> <mi>G</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>Q</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>D</mi> <mi>G</mi> </mrow> <mn>2</mn> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>S</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>D</mi> <mi>G</mi> </mrow> <mn>2</mn> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced>
PQ models are used during photovoltaic module steady-state operation, by grid-connected inverters, grid-connected power can realize that active reactive is independently controlled System;In formula, Pi,DGFor node i photovoltaic module active power output value,For node i photovoltaic active power output predicted value, Qi,DGFor node I photovoltaic modules are idle power generating value, Si,DGFor node i photovoltaic module capacity.
CN201510583592.9A 2015-09-14 2015-09-14 A kind of distributed photovoltaic automatic power generation control method based on sensitivity Active CN105281360B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510583592.9A CN105281360B (en) 2015-09-14 2015-09-14 A kind of distributed photovoltaic automatic power generation control method based on sensitivity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510583592.9A CN105281360B (en) 2015-09-14 2015-09-14 A kind of distributed photovoltaic automatic power generation control method based on sensitivity

Publications (2)

Publication Number Publication Date
CN105281360A CN105281360A (en) 2016-01-27
CN105281360B true CN105281360B (en) 2018-05-11

Family

ID=55149906

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510583592.9A Active CN105281360B (en) 2015-09-14 2015-09-14 A kind of distributed photovoltaic automatic power generation control method based on sensitivity

Country Status (1)

Country Link
CN (1) CN105281360B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102014113262B4 (en) * 2014-09-15 2016-09-15 Sma Solar Technology Ag Method and device for operating a power plant with fluctuating capacity connected to a network generator and at least one load to a limited AC network
CN106058887B (en) * 2016-07-08 2018-10-16 燕山大学 Improve the idle work optimization method of the qualification rate of distribution network voltage containing distributed photovoltaic power
CN106340897B (en) * 2016-09-20 2018-11-06 山东省建设项目环境评审服务中心 A kind of determination method of power distribution network distributed generation resource allowed capacity
CN107317353B (en) * 2017-06-26 2020-09-25 中国电力科学研究院 Voltage control method and system for distribution network containing distributed photovoltaic power generation
CN109217284A (en) * 2017-07-05 2019-01-15 南京理工大学 A kind of reconstruction method of power distribution network based on immune binary particle swarm algorithm
CN107994587A (en) * 2017-12-06 2018-05-04 国网江苏省电力有限公司淮安供电分公司 A kind of high permeability photovoltaic power distribution network voltage control method on the spot
CN108494023A (en) * 2018-04-06 2018-09-04 冯辉 A kind of active output control system of photovoltaic plant and method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101373905A (en) * 2008-05-16 2009-02-25 清华大学 Visualization method for voltage stability domain of electric power system
CN102611118A (en) * 2012-03-14 2012-07-25 清华大学 Method for comprehensively controlling reactive voltage of wind farm with imported prediction method
CN104319775A (en) * 2014-09-27 2015-01-28 国家电网公司 Reactive power voltage control method applicable to large-scale wind power integration
CN104348188A (en) * 2014-11-21 2015-02-11 四川慧盈科技有限责任公司 Distributed generation running and monitoring method
CN104361200A (en) * 2014-10-10 2015-02-18 国家电网公司 Sensitivity-based reactive compensation and location optimization method
CN104753061A (en) * 2015-03-05 2015-07-01 中国农业大学 Distributed type power supply accessed into power distribution network and microgrid group zone control method and microgrid group zone control system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101373905A (en) * 2008-05-16 2009-02-25 清华大学 Visualization method for voltage stability domain of electric power system
CN102611118A (en) * 2012-03-14 2012-07-25 清华大学 Method for comprehensively controlling reactive voltage of wind farm with imported prediction method
CN104319775A (en) * 2014-09-27 2015-01-28 国家电网公司 Reactive power voltage control method applicable to large-scale wind power integration
CN104361200A (en) * 2014-10-10 2015-02-18 国家电网公司 Sensitivity-based reactive compensation and location optimization method
CN104348188A (en) * 2014-11-21 2015-02-11 四川慧盈科技有限责任公司 Distributed generation running and monitoring method
CN104753061A (en) * 2015-03-05 2015-07-01 中国农业大学 Distributed type power supply accessed into power distribution network and microgrid group zone control method and microgrid group zone control system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
兆瓦级光伏电站并网技术的研究;王佳;《CNKI优秀硕士论文全文数据库》;20101231;第13至53页 *

Also Published As

Publication number Publication date
CN105281360A (en) 2016-01-27

Similar Documents

Publication Publication Date Title
Wong et al. Review on the optimal placement, sizing and control of an energy storage system in the distribution network
Shivashankar et al. Mitigating methods of power fluctuation of photovoltaic (PV) sources–A review
Badwawi et al. A review of hybrid solar PV and wind energy system
Sahu et al. A hybrid DE–PS algorithm for load frequency control under deregulated power system with UPFC and RFB
US9436200B2 (en) Use of distributed generator (DG) inverters as statcoms for decreasing line losses
CN103150606B (en) A kind of distributed power source optimal load flow optimization method
CN104242337B (en) The real time coordination control method of photovoltaic microgrid system
Iqbal et al. Optimal placement of DG and DSTATCOM for loss reduction and voltage profile improvement
Zhao et al. Cooperative optimal control of battery energy storage system under wind uncertainties in a microgrid
CN103151797B (en) Multi-objective dispatching model-based microgrid energy control method under grid-connected operation mode
CN100578911C (en) Electric voltage idle-work fast control method of speed-changing constant frequency wind generator group wind power station
Enslin Network impacts of high penetration of photovoltaic solar power systems
CN102522781B (en) Method for participating area control error (ACE) control by uniformly modeling wind power plant and thermal power plant
CN103580030B (en) Parallel networking type photovoltaic power station reactive voltage control method and system
Guan et al. Distribution system feeder reconfiguration considering different model of DG sources
Sheng et al. Research and practice on typical modes and optimal allocation method for PV-Wind-ES in Microgrid
Han et al. Hierarchical energy management for PV/hydrogen/battery island DC microgrid
García-Triviño et al. Long-term optimization based on PSO of a grid-connected renewable energy/battery/hydrogen hybrid system
Satapathy et al. Stability improvement of PV-BESS diesel generator-based microgrid with a new modified harmony search-based hybrid firefly algorithm
CN103441506B (en) Method for multi-target coordination reactive power optimization control of distributed wind farm in different time scales
CN102611118A (en) Method for comprehensively controlling reactive voltage of wind farm with imported prediction method
CN106058887B (en) Improve the idle work optimization method of the qualification rate of distribution network voltage containing distributed photovoltaic power
CN105186556B (en) Based on the large-sized photovoltaic power station idle work optimization method for improving immunity particle cluster algorithm
CN103280821A (en) Multi-period dynamic reactive power optimization method of intelligent power distribution system
CN104362648B (en) A kind of idle phase modulation method of photovoltaic plant

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
GR01 Patent grant