CN111342501B - Reactive power control method for microgrid with distributed power supply - Google Patents

Reactive power control method for microgrid with distributed power supply Download PDF

Info

Publication number
CN111342501B
CN111342501B CN202010112982.9A CN202010112982A CN111342501B CN 111342501 B CN111342501 B CN 111342501B CN 202010112982 A CN202010112982 A CN 202010112982A CN 111342501 B CN111342501 B CN 111342501B
Authority
CN
China
Prior art keywords
wind
power
value
output
station
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
CN202010112982.9A
Other languages
Chinese (zh)
Other versions
CN111342501A (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.)
Jiashan County Power Supply Co Of State Grid Zhejiang Electric Power Co ltd
Jiaxing Power Supply Co of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Jiashan County Power Supply Co Of State Grid Zhejiang Electric Power Co ltd
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 Jiashan County Power Supply Co Of State Grid Zhejiang Electric Power Co ltd, Jiaxing Power Supply Co of State Grid Zhejiang Electric Power Co Ltd filed Critical Jiashan County Power Supply Co Of State Grid Zhejiang Electric Power Co ltd
Priority to CN202010112982.9A priority Critical patent/CN111342501B/en
Publication of CN111342501A publication Critical patent/CN111342501A/en
Application granted granted Critical
Publication of CN111342501B publication Critical patent/CN111342501B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/50Controlling the sharing of the out-of-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/18Arrangements for adjusting, eliminating or compensating reactive power in networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • 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
    • 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
    • Y02E70/00Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30Systems combining energy storage with energy generation of non-fossil origin

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to the technical field of power distribution networks, in particular to a reactive power control method of a microgrid with distributed power sources, which comprises the following steps: A) reading an active power predicted value and a reactive power predicted value, and an output predicted value of the wind and light power station; B) calculating a characteristic value, and if the characteristic value is smaller than a set threshold value, outputting active power only by the wind-solar power station, otherwise, outputting the active power and reactive power; C) calculating the deviation probability of the wind and light power station; D) and monitoring the output of the wind and light power station in real time, and scheduling the energy storage station and the wind and light power station according to the deviation probability. The substantial effects of the invention are as follows: by reasonably scheduling reactive power output, the stability of the microgrid can be improved, the cost of newly adding traditional energy sources for adapting to the access of renewable energy sources is obviously reduced, the access of the renewable energy sources is favorably expanded, the microgrid with high permeability and renewable energy sources is suitable, and the running stability of the power distribution network is improved.

Description

Reactive power control method for microgrid with distributed power supply
Technical Field
The invention relates to the technical field of power distribution networks, in particular to a reactive power control method for a microgrid with a distributed power supply.
Background
The distributed power supply refers to a power supply with a voltage grade of 35kV or below, which is not directly connected with a centralized power transmission system, and mainly comprises power generation equipment and an energy storage device. The distributed power generation can widely utilize clean renewable energy sources, and reduce the consumption of fossil energy and the emission of harmful gases. And the distributed power generation is positioned at the user side and close to the load center, so that the construction cost and the loss of the transmission and distribution network are greatly reduced. The novel energy production system has the advantages of low pollution emission, flexibility, convenience, high reliability and high efficiency. Distributed Generation (DG) devices can be classified into cogeneration, internal combustion engine set Generation, gas turbine Generation, small hydroelectric power Generation, wind power Generation, solar photovoltaic power Generation, fuel cells, and the like, depending on the technology used. At present, with the rapid development of wind power and photovoltaic power generation, the power distribution network in partial areas has the trend of gradually increasing the density of wind and photovoltaic station type distributed power supplies, and the proportion of the distributed power supplies in the power supply of the power distribution network is large. Wind and photoelectric stations are clean energy sources, and cannot damage the environment during power supply, so that the wind and photoelectric stations are mainly used for outputting active power. However, when the permeability is high, the power factor of the conventional energy equipment such as the transformer works in an unreasonable range, and the working efficiency and the working stability of the equipment such as the transformer are reduced. At present, the stable operation of the power grid is a main factor for limiting the further expansion of the scale of the wind and photovoltaic station. Therefore, a control method suitable for a high-penetration distribution network needs to be researched.
For example, chinese patent CN105470991B, published 2018, 8, 28, an output current-limiting control method for an inverter-type distributed power supply, which can ensure the maximization of the output active power of the inverter-type distributed power supply both under symmetric faults and asymmetric faults, and at the same time, improve the positive-sequence voltage amplitude of the output of the inverter-type distributed power supply, reduce the negative-sequence voltage amplitude, and reduce the asymmetry of the output voltage of the inverter-type distributed power supply, under different current-limiting methods: when the asymmetric fault occurs, the high voltage output level of the inverter type distributed power supply is ensured, so that the low-voltage protection action is not easy to trigger. Meanwhile, for the negative sequence overcurrent protection, the current limiting method provided by the method limits the output of the negative sequence current, so that the negative sequence overcurrent protection is not easy to trigger. But it can only solve the problem of grid stability degradation due to voltage and cannot solve the grid operation problem due to power factor.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the technical problem that the permeability of renewable energy sources in the existing power distribution network is not high is solved. A reactive power control method of a microgrid with distributed power supplies is provided, and equipment in the microgrid can work at a more reasonable power factor.
To solve the problems mentioned aboveThe technical problem is that the technical scheme adopted by the invention is as follows: a reactive power control method of a microgrid with a distributed power supply, wherein the distributed power supply comprises an energy storage station and a wind-solar power station, and the method comprises the following steps: A) dividing a day into n periods t j,j∈[1,n] Reading the next time interval t j+1 Active power prediction value of load
Figure GDA0003614820290000021
And reactive power prediction
Figure GDA0003614820290000022
Obtaining wind-solar power station G i,i∈[1,l] In the next time period t j+1 Predicted value of output
Figure GDA0003614820290000023
l is the number of wind and light power stations; B) calculating the next time period t j+1 Characteristic value
Figure GDA0003614820290000024
Figure GDA0003614820290000025
If the characteristic value
Figure GDA0003614820290000026
Less than a first set threshold lambda thr Then t is j+1 Wind-solar power station G in time interval i,i∈[1,l] Only active power is output, otherwise, calculation is carried out
Figure GDA0003614820290000027
Wind-light power station G i,i∈[1,l] Output active power
Figure GDA0003614820290000028
And reactive power
Figure GDA0003614820290000029
And satisfy
Figure GDA00036148202900000210
Figure GDA00036148202900000211
t j+1 In time period, wind-light power station G i Reactive power of
Figure GDA00036148202900000212
Keeping the same; C) wind and light power station G i,i∈[1,l] Real-time force of
Figure GDA00036148202900000213
Deviation from predicted value
Figure GDA00036148202900000214
Probability of deviation of
Figure GDA00036148202900000215
Tau is the rate of deviation and is the rate of deviation,
Figure GDA00036148202900000216
σ is the deviation rate τ | t j+1 The probability of occurrence; D) time period t j+1 Internal and real-time monitoring wind-light power station G i,i∈[1,l] Real-time force of
Figure GDA00036148202900000217
Deviation rate of (τ | t) j+1 If deviation probability
Figure GDA00036148202900000218
Then the energy storage station E is added i,i∈[1,m] Charging power or reducing energy storage station E i,i∈[1,m] M is the number of energy storage stations, otherwise, the wind-light power station G is increased according to a set range i,i∈[1,l] Output reactive power
Figure GDA00036148202900000219
Until real-time output
Figure GDA00036148202900000220
Deviation rate τ | t of (1) j+1 Corresponding deviation probability
Figure GDA00036148202900000221
Fall back to sigma thr The following. At present, a wind and light power station only outputs active power, so that the operation of a distribution network is easy to be unstable, the upper limit of the output of the wind and light power station is often limited, wind and light abandon is caused, and the access scale of the wind and light power station is also limited. By reasonably scheduling reactive power output, the method can improve the stability of the microgrid, obviously reduce the cost of newly adding traditional energy sources for adapting to the access of renewable energy sources, and is suitable for the microgrid with renewable energy sources with high permeability.
Preferably, in step B), the wind-solar power station G is predicted i The next time period t i+1 Output prediction
Figure GDA00036148202900000222
The method comprises the following steps: B11) enumerating influences on wind-solar power station G i,i∈[1,l] The factors of the output, constitute factor vectors; B12) selecting a wind-solar power station G i The output interval is set
Figure GDA00036148202900000223
Is equally divided into N small intervals, and the cell size is equal to the cell size,
Figure GDA00036148202900000224
for wind-solar power station G i Maximum output of (c); B13) wind-light power station G i Operating for a period of time to acquire the values of the factors in the factor vector and the wind-solar power station G at a set frequency i Output, collection of factor vectors and wind-light power station G i The output of the system is associated with the cells to form sample data; B14) establishing a neural network model, wherein the input of the neural network model is a factor vector, the output is the probability of N small intervals, and the neural network model is trained and verified by using the sample data obtained in the step B13); B15) in the next time period t i+1 Setting time before starting, measuring wind-light power station G i The value of the factor in the factor vector is input into the neural network model, and the mean value among the cells output by the neural network model is used as the output prediction
Figure GDA0003614820290000031
The value of (c). The ultra-short term output prediction of the wind and light power station has various published prediction models, the ultra-short term output prediction is obtained by the models according to the commonality of the wind and light power stations, the evaluation index can reflect the deviation between the model prediction result and the actual result under the actual condition, the difference between the wind and light power stations and the difference between the same wind and light power station in different time periods is reflected, the error of the model prediction is corrected, and the accuracy of the output prediction is improved.
Preferably, step B13) further comprises factor screening, comprising: B131) associating the factor vector with the corresponding wind-solar power station G i Form a correlation vector V g,g∈[1,M] =<a 1 ,a 2 ,…,a k ,S>M is the number of the collected factor vector values, k is the number of the factors, and the difference value of each factor and the output force is calculated among the associated vectors to form a difference value vector V g-k =<Δa 1 ,Δa 2 ,…,Δa k ,ΔS>Wherein g is [1, M ]],k∈[1,M]G ≠ k, difference vector V g-k Each element in the vector is an association vector V g And V k Difference values of corresponding elements to obtain at least M difference value vectors V g-k (ii) a B132) The difference vector V g-k Dividing the value of each element by the value corresponding to the output force, converting the output force into a unit 1, and obtaining a corrected difference vector
Figure GDA0003614820290000032
B133) Calculating a corrected difference vector V' g-k Standard deviation of each element in
Figure GDA0003614820290000033
If it is
Figure GDA0003614820290000034
If the value is larger than the first set threshold value, the corresponding factor a is used i And (4) screening. When the output of the wind and light power station is not changed greatly, irrelevant factors can reflect the characteristic of random fluctuation, and screening can be performed through standard deviation.
Preferably, in step D), the energy storage station E is controlled i,i∈[1,m] The method for increasing the charging power comprises the following steps: D11) computing energy storage station E i,i∈[1,m] Total increased charging power of
Figure GDA0003614820290000035
k is a set margin coefficient, k>1; D12) will be of period t j+1 Dividing the evaluation function into N small periods
Figure GDA0003614820290000036
Where z represents the number of small cycles,
Figure GDA0003614820290000037
for feeder i during time period t j+1 The average load of the z-th small period of time,
Figure GDA0003614820290000038
h is the number of feeders, i is the upper limit of the load of the feeder i; D13) at a set time before the beginning of the z-th small period, starting from the (z-1) th small period to the wind-light power station G within the current time i Real-time force of
Figure GDA0003614820290000039
Average value calculation of
Figure GDA00036148202900000310
A value of (d); D14) using an optimization algorithm, an evaluation function is obtained
Figure GDA00036148202900000311
Energy storage station E with the smallest value i,i∈[1,m] The charging power is increased in real time, and the process returns to step D13) at a predetermined timing before the start of the next small cycle. The energy storage equipment is scheduled and coordinated on site in a small period, and the energy utilization efficiency is improved.
Preferably, in step B), the wind-solar power station G i,i∈[1,l] Active power of
Figure GDA0003614820290000041
And reactive power
Figure GDA0003614820290000042
Distribution method packageComprises the following steps: B11) establishing an evaluation function
Figure GDA0003614820290000043
Figure GDA0003614820290000044
For feeder i during time period t j+1 The average active power transferred in-between,
Figure GDA0003614820290000045
for feeder i during time period t j+1 The average reactive power transferred in-between,
Figure GDA0003614820290000046
h is the number of feeders, i is the upper limit of the load of the feeder i; B12) using an optimization algorithm, an evaluation function is obtained
Figure GDA0003614820290000047
Minimum value wind-solar power station G i,i∈[1,l] Reactive power of
Figure GDA0003614820290000048
Value, active power
Figure GDA0003614820290000049
Preferably, in step B), the wind-solar power station is calculated
Figure GDA00036148202900000410
Operating power factor of
Figure GDA00036148202900000411
Figure GDA00036148202900000412
Time period t j+1 Internal active power
Figure GDA00036148202900000413
And reactive power
Figure GDA00036148202900000414
Is distributed to
Figure GDA00036148202900000415
Figure GDA00036148202900000416
Always true, λ' thr Is a second set threshold value, λ' thrthr . The working state of the wind and light power station is ensured to be at a better level, and excessive renewable energy waste is avoided.
Preferably, in step C), the deviation probability is calculated
Figure GDA00036148202900000417
The method comprises the following steps: C11) will be of period t j Dividing into M small periods, and counting wind-light power station G in each small period i,i∈[1,l] Mean value of real-time forces of
Figure GDA00036148202900000418
r∈[1,M](ii) a C12) Will be of period t j+1 The former w time periods execute the step C11), and the wind and light power station G in each small period in each time period is obtained i,i∈[1,l] Mean value of real-time forces of
Figure GDA00036148202900000419
u∈[j-w,j],r∈[1,M](ii) a C13) Statistics of
Figure GDA00036148202900000420
Maximum and minimum values of
Figure GDA00036148202900000421
Section of will
Figure GDA00036148202900000422
Equally dividing the data into a plurality of value intervals, and respectively counting the data falling into each interval
Figure GDA00036148202900000423
Of each value interval
Figure GDA00036148202900000424
The ratio of the number of the wind-solar power station G to the number w.M is used as the wind-solar power station G i,i∈[1,l] Real-time force of
Figure GDA00036148202900000425
Deviation probability corresponding to falling value interval
Figure GDA00036148202900000426
Through deviation probability, the possible duration of the output fluctuation of the current wind and light power station is obtained, the deviation probability corresponding to the output level is large, the duration of the output level is longer, therefore, a scheduling mode more suitable for a long time needs to be adopted for the fluctuation output level, if the possibility that the duration of the output level is short is large, the influence of the output fluctuation can be eliminated through short-term compensation of an energy storage device, and flexible scheduling is achieved.
The substantial effects of the invention are as follows: by reasonably scheduling reactive power output, the stability of the microgrid can be improved, the newly increased cost of the traditional energy sources for adapting to the access of renewable energy sources is obviously reduced, the access of the renewable energy sources is favorably expanded, the microgrid with high permeability and renewable energy sources is suitable, and the running stability of the power distribution network is improved.
Drawings
Fig. 1 is a flowchart illustrating a method for controlling a reactive power according to an embodiment.
FIG. 2 is a flow chart of a wind-solar power plant output prediction method according to an embodiment.
Fig. 3 is a flowchart illustrating a method for increasing charging power of a power storage station according to an embodiment of the invention.
Detailed Description
The following provides a more detailed description of the present invention, with reference to the accompanying drawings.
The first embodiment is as follows:
a reactive power control method for a microgrid with distributed power sources, wherein the distributed power sources include an energy storage station and a wind-light power station, as shown in fig. 1, the embodiment includes the following steps:
A) dividing a day into n periods t j,j∈[1,n] Reading the next time interval t j+1 Active power prediction value of load
Figure GDA0003614820290000051
And reactive power prediction
Figure GDA0003614820290000052
Obtaining wind-solar power station G i,i∈[1,l] In the next time period t j+1 Predicted value of output
Figure GDA0003614820290000053
And l is the number of the wind-solar power stations.
B) Calculating the next time period t j+1 Characteristic value
Figure GDA0003614820290000054
If the characteristic value
Figure GDA0003614820290000055
Less than a first set threshold lambda thr Then t is j+1 Wind-solar power station G in time interval i,i∈[1,l] Only active power is output, otherwise, calculation is carried out
Figure GDA0003614820290000056
Wind-light power station G i,i∈[1,l] Output active power
Figure GDA0003614820290000057
And reactive power
Figure GDA0003614820290000058
And satisfy
Figure GDA0003614820290000059
t j+1 In time period, wind-light power station G i Of reactive power
Figure GDA00036148202900000510
Remain unchanged.
Wind-light power station G i,i∈[1,l] Active power of
Figure GDA00036148202900000511
And reactive power
Figure GDA00036148202900000512
The distribution method comprises the following steps: B11) establishing an evaluation function
Figure GDA00036148202900000513
Figure GDA00036148202900000514
For feeder i at time period t j+1 The average active power transferred in-between,
Figure GDA00036148202900000515
for feeder i at time period t j+1 The average reactive power transferred in-between,
Figure GDA00036148202900000516
h is the number of feeders, i is the upper limit of the load of the feeder i; B12) using an optimization algorithm, an evaluation function is obtained
Figure GDA00036148202900000517
Minimum value wind-solar power station G i,i∈[1,l] Of reactive power
Figure GDA00036148202900000518
Value, active power
Figure GDA00036148202900000519
Wind-solar power station G i,i∈[1,l] Operating power factor of
Figure GDA00036148202900000520
Time period t j+1 Internal active power
Figure GDA00036148202900000521
And reactive power
Figure GDA00036148202900000522
Is distributed to
Figure GDA00036148202900000523
Always true, λ' thr Is a second set threshold value, λ' thrthr . The working state of the wind and light power station is ensured to be at a better level, and excessive renewable energy waste is avoided. As shown in fig. 2, predictive wind-solar plant G i The next time period t i+1 Output prediction
Figure GDA00036148202900000524
The method comprises the following steps: B11) enumerating influences on wind-solar power station G i,i∈[1,l] The factors of the output, constitute factor vectors; B12) selecting a wind-solar power station G i The output interval is set
Figure GDA0003614820290000061
Is equally divided into N small intervals,
Figure GDA0003614820290000062
for wind-solar power station G i The maximum output of (c); B13) wind-light power station G i Operating for a period of time to acquire values of factors in the factor vector and the wind and light power station G at a set frequency i The output of (1), collecting the factor vector and the wind-light power station G i The output of (2) is associated with the cells to form sample data; B14) establishing a neural network model, wherein the input of the neural network model is a factor vector, the output is the probability of N small intervals, and the neural network model is trained and verified by using the sample data obtained in the step B13); B15) in the next time period t i+1 Setting time before starting, measuring wind-light power station G i The value of the factor in the factor vector is input into the neural network model, and the mean value among the cells output by the neural network model is used as the output prediction
Figure GDA0003614820290000063
The value of (c). The ultra-short term output prediction of the wind and light power station has various published prediction models, but the ultra-short term output prediction is obtained by the models according to the commonality of the wind and light power station, and the evaluation index can reflect the actual outputAnd the deviation of the model prediction result and the actual result under the condition reflects the difference between the wind and light power stations and the difference of the same wind and light power station in different time periods, so that the error of the model prediction is corrected, and the accuracy of the output prediction is improved.
Step B13) further comprises factor screening comprising: B131) associating the factor vector with the corresponding wind-solar power station G i Form a correlation vector V g,g∈[1,M] =<a 1 ,a 2 ,…,a k ,S>M is the number of the collected factor vector values, k is the number of the factors, and the difference value of each factor and the output force is calculated among the associated vectors to form a difference value vector V g-k =<Δa 1 ,Δa 2 ,…,Δa k ,ΔS>Wherein g is [1, M ]],k∈[1,M]G ≠ k, difference vector V g-k Each element in the vector is an association vector V g And V k Difference values of corresponding elements to obtain at least M difference value vectors V g-k (ii) a B132) The difference vector V g-k Dividing the value of each element by the value corresponding to the output force, converting the output force into a unit 1, and obtaining a corrected difference vector
Figure GDA0003614820290000064
B133) Calculating a corrected difference vector V g-k Standard deviation of each element in
Figure GDA0003614820290000065
If it is
Figure GDA0003614820290000066
If the value is larger than the first set threshold value, the corresponding factor a is used i And (6) screening. When the output of the wind and light power station is not changed greatly, irrelevant factors can reflect the characteristic of random fluctuation, and screening can be performed through standard deviation.
C) Wind and light power station G i,i∈[1,l] Real time force output
Figure GDA0003614820290000067
Deviation from predicted value
Figure GDA0003614820290000068
Deviation probability of (2)
Figure GDA0003614820290000069
Tau is the rate of deviation and is the rate of deviation,
Figure GDA00036148202900000610
σ is the deviation rate τ | t j+1 The probability of occurrence. Calculating the probability of departure
Figure GDA00036148202900000611
The method comprises the following steps: C11) will be of period t j Dividing into M small periods, and counting wind-light power station G in each small period i,i∈[1,l] Mean value of real-time forces
Figure GDA00036148202900000612
r∈[1,M](ii) a C12) Will be of period t j+1 The former w time periods execute the step C11), and the wind-solar power station G in each small period in each time period is obtained i,i∈[1,l] Mean value of real-time forces of
Figure GDA00036148202900000613
u∈[j-w,j],r∈[1,M](ii) a C13) Statistics of
Figure GDA00036148202900000614
Maximum and minimum values of
Figure GDA00036148202900000615
Section of will
Figure GDA00036148202900000616
Equally dividing the data into a plurality of value intervals, and respectively counting the data falling into each interval
Figure GDA0003614820290000071
Of each value interval
Figure GDA0003614820290000072
The ratio of the number of the wind-solar power station G to the number w.M is used as the wind-solar power station G i,i∈[1,l] Real-time force of
Figure GDA0003614820290000073
Deviation probability corresponding to falling value interval
Figure GDA0003614820290000074
Through deviation probability, the possible duration of the output fluctuation of the current wind and light power station is obtained, the deviation probability corresponding to the output level is large, the duration of the output level is longer, therefore, a scheduling mode more suitable for a long time needs to be adopted for the fluctuation output level, if the possibility that the duration of the output level is short is large, the influence of the output fluctuation can be eliminated through short-term compensation of an energy storage device, and flexible scheduling is achieved.
D) Time period t j+1 Internal and real-time monitoring wind and light power station G i,i∈[1,l] Real time force output
Figure GDA0003614820290000075
Deviation rate τ | t of (1) j+1 If the probability of deviation is high
Figure GDA0003614820290000076
Then add the energy storage station E i,i∈[1,m] Charging power or reducing energy storage station E i,i∈[1,m] M is the number of energy storage stations, otherwise, the wind-light power station G is increased according to a set range i,i∈[1,l] Reactive power of output
Figure GDA0003614820290000077
Until real-time output
Figure GDA0003614820290000078
Deviation rate of (τ | t) j+1 Corresponding deviation probability
Figure GDA0003614820290000079
Fall back to sigma thr The following. As shown in fig. 3, the energy storage station E is controlled i,i∈[1,m] The method for increasing the charging power comprises the following steps: D11) computing energy storage station E i,i∈[1,m] Total increased charging power of
Figure GDA00036148202900000710
Figure GDA00036148202900000711
k is a set margin coefficient, k>1; D12) will be of period t j+1 Dividing the evaluation function into N small periods
Figure GDA00036148202900000712
Where z represents a small cycle number and,
Figure GDA00036148202900000713
for feeder i during time period t j+1 The average load of the z-th small period of (c),
Figure GDA00036148202900000714
h is the number of feeders, i is the upper limit of the load of the feeder i; D13) at a set time before the beginning of the z-th small period, the wind-light power station G within the (z-1) th small period to the current time i Real-time force of
Figure GDA00036148202900000715
Mean value calculation of
Figure GDA00036148202900000716
Figure GDA00036148202900000717
A value of (d); D14) using an optimization algorithm, an evaluation function is obtained
Figure GDA00036148202900000718
Energy storage station E with the smallest value i,i∈[1,m] The charging power is increased in real time, and the process returns to step D13) at a predetermined timing before the start of the next small cycle. The energy storage equipment is scheduled and coordinated on site in a small period, and the energy utilization efficiency is improved.
At present, the wind and light power station only outputs active power, so that the operation of a distribution network is easy to be unstable, the upper limit of the output of the wind and light power station is often limited, wind and light abandon is caused, and the access scale of the wind and light power station is also limited. By reasonably scheduling reactive power output, the method can improve the stability of the microgrid, obviously reduce the cost of newly adding traditional energy sources for adapting to the access of renewable energy sources, and is suitable for the microgrid with renewable energy sources with high permeability.
The above-described embodiments are only preferred embodiments of the present invention, and are not intended to limit the present invention in any way, and other variations and modifications may be made without departing from the spirit of the invention as set forth in the claims.

Claims (7)

1. A reactive power control method of a microgrid with distributed power supplies, wherein the distributed power supplies comprise an energy storage station and a wind-solar power station,
the method comprises the following steps:
A) dividing a day into n periods t j,j∈[1,n] Reading the next time interval t j+1 Active power prediction value of load
Figure FDA0003614820280000011
And reactive power prediction
Figure FDA0003614820280000012
Obtaining a wind-solar power station G i,i∈[1,l] In the next time period t j+1 Predicted value of output
Figure FDA0003614820280000013
l is the number of wind and light power stations;
B) calculating the next time period t j+1 Characteristic value
Figure FDA0003614820280000014
If the characteristic value
Figure FDA0003614820280000015
Less than a first set threshold lambda thr Then t is j+1 Wind-solar power station G in time interval i,i∈[1,l] Only active power is output, otherwise, calculation is carried out
Figure FDA0003614820280000016
Wind-light power station G i,i∈[1,l] Output active power
Figure FDA0003614820280000017
And reactive power
Figure FDA0003614820280000018
And satisfy
Figure FDA0003614820280000019
t j+1 In time period, wind-light power station G i Reactive power of
Figure FDA00036148202800000110
Keeping the same;
C) wind and light power station G i,i∈[1,l] Real time force output
Figure FDA00036148202800000111
Deviation from predicted value
Figure FDA00036148202800000112
Probability of deviation of
Figure FDA00036148202800000113
Tau is the rate of deviation and is the rate of deviation,
Figure FDA00036148202800000114
σ is the deviation rate τ | t j+1 The probability of occurrence;
D) time period t j+1 Internal and real-time monitoring wind and light power station G i,i∈[1,l] Real time force output
Figure FDA00036148202800000115
Deviation rate of (τ | t) j+1 If deviation probability
Figure FDA00036148202800000116
Then the energy storage station E is added i,i∈[1,m] Charging power or reducing energy storage station E i,i∈[1,m] M is the number of energy storage stations, otherwise, the wind-light power station G is increased according to a set range i,i∈[1,l] Reactive power of output
Figure FDA00036148202800000117
Until real-time output
Figure FDA00036148202800000118
Deviation rate of (τ | t) j+1 Corresponding deviation probability
Figure FDA00036148202800000119
Fall back to sigma thr The following.
2. The reactive power control method of claim 1, wherein the reactive power control method comprises the steps of,
in step B), the wind and light power station G is predicted i The next time period t i+1 Output prediction
Figure FDA00036148202800000120
The method comprises the following steps:
B11) enumerating influences on wind-solar power station G i,i∈[1,l] The factors of the output, constitute factor vectors;
B12) selecting a wind-solar power station G i Applying a force to the interval
Figure FDA00036148202800000121
Is equally divided into N small intervals,
Figure FDA00036148202800000122
for wind-solar power station G i The maximum output of (c);
B13) wind-light power station G i Operating for a period of time to acquire the values of the factors in the factor vector and the wind-solar power station G at a set frequency i Output, collection of factor vectors and wind and lightPower station G i The output of (2) is associated with the cells to form sample data;
B14) establishing a neural network model, wherein the input of the neural network model is a factor vector, the output is the probability of N small intervals, and the neural network model is trained and verified by using the sample data obtained in the step B13);
B15) in the next time period t i+1 Setting the time before starting, measuring the wind-light power station G i The value of the factor in the factor vector is input into the neural network model, and the mean value among the cells output by the neural network model is used as the output prediction
Figure FDA0003614820280000021
The value of (c).
3. The method of claim 2, wherein the reactive power control of the microgrid comprising a distributed power source is carried out,
step B13) further comprises factor screening comprising:
B131) associating the factor vector with the corresponding wind-solar power station G i Form a correlation vector V g,g∈[1,M] =<a 1 ,a 2 ,…,a k ,S>M is the number of the collected factor vector values, k is the number of the factors, and the difference value of each factor and the output force is calculated among the associated vectors to form a difference value vector V g-k =<Δa 1 ,Δa 2 ,…,Δa k ,ΔS>Wherein g is [1, M ]],k∈[1,M]G ≠ k, difference vector V g-k Each element in the vector is an association vector V g And V k Difference values of corresponding elements to obtain at least M difference value vectors V g-k
B132) The difference vector V g-k Dividing the value of each element by the value corresponding to the output force, converting the output force into a unit 1, and obtaining a corrected difference vector
Figure FDA0003614820280000022
B133) Calculating a corrected difference vector V' g-k Standard deviation of each element in
Figure FDA0003614820280000023
If it is
Figure FDA0003614820280000024
If the value is larger than the first set threshold value, the corresponding factor a is used i And (6) screening.
4. The reactive power control method comprising a distributed power microgrid of claim 1, 2 or 3,
in step D), the energy storage station E is controlled i,i∈[1,m] The method for increasing the charging power comprises the following steps:
D11) computing energy storage station E i,i∈[1,m] Total increased charging power of
Figure FDA0003614820280000025
k is a set margin coefficient, k>1;
D12) Will be of period t j+1 Dividing into N small periods, and establishing evaluation function
Figure FDA0003614820280000026
Where z represents a small cycle number and,
Figure FDA0003614820280000027
for feeder i during time period t j+1 The average load of the z-th small period of (c),
Figure FDA0003614820280000028
h is the number of feeders, i is the upper limit of the load of the feeder i;
D13) at a set time before the beginning of the z-th small period, the wind-light power station G within the (z-1) th small period to the current time i Real time force output
Figure FDA0003614820280000029
Mean value calculation of
Figure FDA00036148202800000210
A value of (d);
D14) using an optimization algorithm, an evaluation function is obtained
Figure FDA0003614820280000031
Energy storage station E with the smallest value i,i∈[1,m] The charging power is increased in real time, and the process returns to step D13) at a predetermined timing before the start of the next small cycle.
5. The reactive power control method comprising a distributed power microgrid of claim 1, 2 or 3,
in step B), wind-solar power station G i,i∈[1,l] Active power of
Figure FDA0003614820280000032
And reactive power
Figure FDA0003614820280000033
The distribution method comprises the following steps:
B11) establishing an evaluation function
Figure FDA0003614820280000034
Figure FDA0003614820280000035
For feeder i at time period t j+1 The average active power transferred in-between,
Figure FDA0003614820280000036
for feeder i at time period t j+1 The average reactive power transferred in-between,
Figure FDA0003614820280000037
h is the number of feeders, i is the upper limit of the load of the feeder i;
B12) using an optimization algorithm, an evaluation function is obtained
Figure FDA0003614820280000038
Minimum value wind-solar power station G i,i∈[1,l] Reactive power of
Figure FDA0003614820280000039
Value, active power
Figure FDA00036148202800000310
6. The reactive power control method comprising a distributed power microgrid of claim 1, 2 or 3,
in step B), calculating the wind-light power station G i,i∈[1,l] Operating power factor of
Figure FDA00036148202800000311
Time period t j+1 Internal active power
Figure FDA00036148202800000312
And reactive power
Figure FDA00036148202800000313
Is distributed to
Figure FDA00036148202800000314
Always true, λ' thr Is a second set threshold value, λ' thrthr
7. The reactive power control method comprising a distributed power microgrid of claim 1, 2 or 3,
in step C), the deviation probability is calculated
Figure FDA00036148202800000315
The method comprises the following steps:
C11) will be of period t j Dividing into M small periods, and counting wind-light power station G in each small period i,i∈[1,l] Mean value of real-time forces
Figure FDA00036148202800000316
C12) Will be of period t j+1 The former w time periods execute the step C11), and the wind-solar power station G in each small period in each time period is obtained i,i∈[1,l] Mean value of real-time forces of
Figure FDA00036148202800000317
C13) Statistics of
Figure FDA00036148202800000318
Maximum and minimum values of
Figure FDA00036148202800000319
Section of will
Figure FDA00036148202800000320
Equally dividing the data into a plurality of value intervals, and respectively counting the data falling into each interval
Figure FDA00036148202800000321
Of each value interval
Figure FDA00036148202800000322
The ratio of the number of the wind-solar power station G to the number w.M is used as the wind-solar power station G i,i∈[1,l] Real-time force of
Figure FDA00036148202800000323
Deviation probability corresponding to falling value interval
Figure FDA00036148202800000324
CN202010112982.9A 2020-02-24 2020-02-24 Reactive power control method for microgrid with distributed power supply Active CN111342501B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010112982.9A CN111342501B (en) 2020-02-24 2020-02-24 Reactive power control method for microgrid with distributed power supply

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010112982.9A CN111342501B (en) 2020-02-24 2020-02-24 Reactive power control method for microgrid with distributed power supply

Publications (2)

Publication Number Publication Date
CN111342501A CN111342501A (en) 2020-06-26
CN111342501B true CN111342501B (en) 2022-09-27

Family

ID=71185558

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010112982.9A Active CN111342501B (en) 2020-02-24 2020-02-24 Reactive power control method for microgrid with distributed power supply

Country Status (1)

Country Link
CN (1) CN111342501B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111404195B (en) * 2020-02-24 2021-08-27 国网浙江嘉善县供电有限公司 Intelligent gateway-based scheduling method for microgrid with distributed power supply
CN111969432A (en) * 2020-07-10 2020-11-20 国网浙江省电力有限公司湖州供电公司 Wind-solar complementary power supply device for ring main unit dehumidification system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104638772A (en) * 2013-11-14 2015-05-20 国家电网公司 Battery energy storage power station energy management method based on wind power prediction
CN107069807A (en) * 2017-03-22 2017-08-18 广西大学 Containing uncertain budget adjust without balance nodes microgrid Robust Scheduling method
KR101786904B1 (en) * 2016-08-17 2017-11-15 인천대학교 산학협력단 Model Predictive Control Apparatus and Method to BESS for Microgrid Control
CN109066744A (en) * 2018-08-08 2018-12-21 中国电力科学研究院有限公司 One kind coordinated scheduling method of power distribution network containing energy storage and system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104638772A (en) * 2013-11-14 2015-05-20 国家电网公司 Battery energy storage power station energy management method based on wind power prediction
KR101786904B1 (en) * 2016-08-17 2017-11-15 인천대학교 산학협력단 Model Predictive Control Apparatus and Method to BESS for Microgrid Control
CN107069807A (en) * 2017-03-22 2017-08-18 广西大学 Containing uncertain budget adjust without balance nodes microgrid Robust Scheduling method
CN109066744A (en) * 2018-08-08 2018-12-21 中国电力科学研究院有限公司 One kind coordinated scheduling method of power distribution network containing energy storage and system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Impacts_of_Stochastic_Wind_Power_and_Storage_Participation_on_Economic_Dispatch_in_Distribution_Systems;Abebe W. Bizuayehu;《IEEE TRANSACTIONS ON SUSTAINABLE ENERGY》;20160731;第1336-1345页 *
考虑电压稳定约束的无功补偿优化配置;程振龙;《电力系统及其自动化学报》;20150430;第86-91页 *
考虑风电消纳的配电网多目标无功协同优化;金尚婷;《广东电力》;20190131;第52-59页 *

Also Published As

Publication number Publication date
CN111342501A (en) 2020-06-26

Similar Documents

Publication Publication Date Title
CN109325621B (en) Park energy internet two-stage optimal scheduling control method
CN110245794B (en) Flexibility-considered double-layer optimization method for central fire storage capacity in multi-energy convergence
CN111342501B (en) Reactive power control method for microgrid with distributed power supply
CN110137952A (en) It is a kind of to realize source-lotus-storage coordinated operation Optimized model and method
CN112671035A (en) Virtual power plant energy storage capacity configuration method based on wind power prediction
CN107359611B (en) Power distribution network equivalence method considering various random factors
CN107392791B (en) Distributed photovoltaic and gas-electricity hybrid capacity planning method and system for multi-energy complementary system
CN114050609B (en) Adaptive robust day-ahead optimization scheduling method for high-proportion new energy power system
CN114301081B (en) Micro-grid optimization method considering storage battery energy storage life loss and demand response
CN111130145A (en) Wind and light installed capacity optimization planning method based on reasonable wind and light abandonment
Han et al. Analysis of economic operation model for virtual power plants considering the uncertainties of renewable energy power generation
CN115719967A (en) Active power distribution network energy storage device optimal configuration method for improving power supply reliability
CN111985844A (en) Day-ahead economic dispatching method for wind power and light energy comprehensive energy system
Zhang et al. Study on the microgrid pluripotent complementary and comprehensive economic optimization planning based on virtual energy theory
Zhang et al. Analysis of peak regulation strategy with considering renewable energy injection and power from outside
CN115528752B (en) Micro-grid group control method and device
Xie et al. Multi-objective stochastic optimization planning method based on microgrid
CN115021333A (en) New energy consumption capability assessment method based on time sequence state
Lu et al. Optimal modeling of integrated energy system demand response operation considering wind power absorption
Hu et al. Multi-Objective Configuration Method of Energy Storage Capacity of New Energy Microgrid Based on Multi-Functional Application
Liu et al. Research on optimal scheduling of wind-solar-storage system based on improved multi-objective particle swarm optimization algorithm
Yu et al. Optimal Allocation of Electrochemical Energy Storage of Source-Grid-Load Sides in Power System Considering Social Benefit
Zhang et al. Multi-objective optimal scheduling of wind-photoelectric-storage virtual power plants considering demand response
CN115882452A (en) New energy consumption capability analysis and evaluation method considering source load uncertainty
Zhang et al. Research on multi-time scale optimal scheduling of integrated energy system based on digital twinning

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
TA01 Transfer of patent application right

Effective date of registration: 20220908

Address after: 314100 No. 858, Chezhan South Road, Luoxing street, Jiashan County, Jiaxing City, Zhejiang Province

Applicant after: Jiashan County power supply company of State Grid Zhejiang Electric Power Co.,Ltd.

Applicant after: JIAXING POWER SUPPLY COMPANY OF STATE GRID ZHEJIANG ELECTRIC POWER Co.,Ltd.

Address before: 314100 No. 858, Chezhan South Road, Luoxing street, Jiashan County, Jiaxing City, Zhejiang Province

Applicant before: STATE GRID ZHEJIANG JIASHAN COUNTY POWER SUPPLY Co.,Ltd.

Applicant before: JIAXING POWER SUPPLY COMPANY OF STATE GRID ZHEJIANG ELECTRIC POWER Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant