CN109978717A - Distributed photovoltaic maximum permeability determines method and device - Google Patents
Distributed photovoltaic maximum permeability determines method and device Download PDFInfo
- Publication number
- CN109978717A CN109978717A CN201910198323.9A CN201910198323A CN109978717A CN 109978717 A CN109978717 A CN 109978717A CN 201910198323 A CN201910198323 A CN 201910198323A CN 109978717 A CN109978717 A CN 109978717A
- Authority
- CN
- China
- Prior art keywords
- distributed photovoltaic
- permeability
- actual load
- setting
- load data
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Abstract
The present invention provides a kind of distributed photovoltaic maximum permeabilities to determine method and device, this method comprises: circulation executes following steps, until target area load prediction stability and default load prediction stability difference within a preset range, determine newest distributed photovoltaic permeability be distributed photovoltaic maximum permeability: obtain setting distributed photovoltaic permeability under actual load data；Obtain the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability；Obtain the load prediction stability of target area；If the difference of the load prediction stability of target area and default load prediction stability within a preset range, does not adjust setting distributed photovoltaic permeability, by the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting.Present invention may determine that distributed photovoltaic maximum permeability, and load prediction stability is considered, it can effectively support load prediction and power balance arrangement.
Description
Technical field
The present invention relates to distributed photovoltaic power generation fields more particularly to a kind of distributed photovoltaic maximum permeability to determine method
And device.
Background technique
The distributed energy of " generating power for their own use " or " remaining electricity online " formula as multiple spot access power distribution network, distributed photovoltaic go out
For power by multifactor impacts such as weather conditions, power output has many characteristics, such as that randomness is big, fluctuation is strong, seasonal periodicity variation, it is difficult to
Etc. being all conventional generator group.Under permeability of the distributed photovoltaic in regional distribution network increasingly increased background, scale
Distributed photovoltaic cluster, or the small-scale terrestrial photovoltaic plant of similar access 35kV voltage class therewith adjusts area
Each time dimensions such as medium-term and long-term distribution planning, a few days ago power balance arrangement and the real time execution control of degree mechanism routine produce
It significantly affects.Although present various circles of society have carried out more research to distributed photovoltaic access distribution influence, study mostly all
It is to be influenced from distribution dynamic stability angle analysis distributed photovoltaic bring, and distributed photovoltaic maximum in area is assessed with this
Permeability, but influence of the distributed photovoltaic access to the load prediction a few days ago of regional scheduling institution, power balance arrangement is ignored, it is difficult
To support the power grid operation after large-scale distributed photovoltaic access.
Summary of the invention
The present invention proposes that a kind of distributed photovoltaic maximum permeability determines method, to determine the infiltration of distributed photovoltaic maximum
Rate, and load prediction stability is considered, it can effectively support load prediction and power balance arrangement, the maximum permeability is based on
Load prediction stability determines, comprising:
According to the distributed photovoltaic active power data and actual load data of target area, the original of target area is obtained
Load data；
Circulation executes following steps, until the difference of the load prediction stability of target area and default load prediction stability
Value within a preset range, determines that newest distributed photovoltaic permeability is distributed photovoltaic maximum permeability:
According to the distributed photovoltaic active power data and original loads data of target area, setting distributed photovoltaic is obtained
Actual load data under permeability；
Actual load data under setting distributed photovoltaic permeability are decomposed, setting distributed photovoltaic infiltration is obtained
The low frequency component and high fdrequency component of the sequence of actual load data under rate；
According to the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, obtain
Obtain the load prediction stability of target area；
If the difference of the load prediction stability of target area and default load prediction stability within a preset range, is not adjusted
Section setting distributed photovoltaic permeability, by the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting.
The present invention proposes a kind of distributed photovoltaic maximum permeability determining device, to determine the infiltration of distributed photovoltaic maximum
Rate, and load prediction stability is considered, it can effectively support load prediction and power balance arrangement, the maximum permeability is based on
Load prediction stability determines, comprising:
Original loads data obtaining module, for the distributed photovoltaic active power data and actual load according to somewhere
Data obtain the original loads data of target area；
Distributed photovoltaic maximum permeability determining module executes following steps for recycling, until the load of target area
The difference of prediction stability and default load prediction stability within a preset range, determines that newest distributed photovoltaic permeability is
Distributed photovoltaic maximum permeability:
According to the distributed photovoltaic active power data and original loads data of target area, setting distributed photovoltaic is obtained
Actual load data under permeability；
Actual load data under setting distributed photovoltaic permeability are decomposed, setting distributed photovoltaic infiltration is obtained
The low frequency component and high fdrequency component of the sequence of actual load data under rate；
According to the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, obtain
Obtain the load prediction stability of target area；
If the difference of the load prediction stability of target area and default load prediction stability within a preset range, is not adjusted
Section setting distributed photovoltaic permeability, by the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting.
The embodiment of the present invention also proposed a kind of computer equipment, including memory, processor and storage are on a memory
And the computer program that can be run on a processor, the processor are realized above-mentioned based on load when executing the computer program
The distributed photovoltaic maximum permeability of prediction stability determines method.
The embodiment of the present invention also proposed a kind of computer readable storage medium, the computer-readable recording medium storage
Have and executes the above-mentioned computer program for determining method based on the distributed photovoltaic maximum permeability of load prediction stability.
In embodiments of the present invention, according to the distributed photovoltaic active power data and actual load data of target area,
Obtain the original loads data of target area；Circulation executes following steps, until the load prediction stability of target area and pre-
If the difference of load prediction stability is within a preset range, determine that newest distributed photovoltaic permeability is that distributed photovoltaic is maximum
Permeability: according to the distributed photovoltaic active power data and original loads data of target area, setting distributed photovoltaic is obtained
Actual load data under permeability；Actual load data under setting distributed photovoltaic permeability are decomposed, are set
Determine the low frequency component and high fdrequency component of the sequence of the actual load data under distributed photovoltaic permeability；According to the distributed light of setting
The low frequency component and high fdrequency component for lying prostrate the sequence of the actual load data under permeability, the load prediction for obtaining target area are stablized
Degree；If not within a preset range, adjusting is set the difference of the load prediction stability of target area and default load prediction stability
Distributed photovoltaic permeability is determined, by the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting.In this hair
In bright embodiment, by setting the load prediction stability of distributed photovoltaic computing permeability target area, if target area
The difference of load prediction stability and default load prediction stability within a preset range, does not adjust setting distributed photovoltaic infiltration
Rate executes the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting repeatedly, until target area
Within a preset range, therefore, above procedure considers target to the difference of load prediction stability and default load prediction stability
The load prediction stability in region, and finally determining distributed photovoltaic maximum permeability have passed through successive ignition, precision is high, can
Effectively support power balance arrangement is executed repeatedly to the load prediction stability of target area and default load prediction stability
Difference within a preset range, therefore can effectively support load prediction.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.In the accompanying drawings:
Fig. 1 is the flow chart that distributed photovoltaic maximum permeability determines method in the embodiment of the present invention；
Fig. 2 is the schematic diagram data of Hebei province Langfang Prefecture in the embodiment of the present invention；
Fig. 3 is the schematic diagram of each component of the sequence of actual load data when setting distributed photovoltaic permeability as 0%；
Fig. 4 is showing for each component of the sequence of actual load data when setting distributed photovoltaic permeability as 12.77%
It is intended to；
Fig. 5 is the signal of each component of the sequence of actual load data when setting distributed photovoltaic permeability as 20%
Figure；
Fig. 6 is the signal of each component of the sequence of actual load data when setting distributed photovoltaic permeability as 40%
Figure；
Fig. 7 is the signal of each component of the sequence of actual load data when setting distributed photovoltaic permeability as 60%
Figure；
Fig. 8 is the signal of the low frequency component of the sequence of the actual load data under different set distributed photovoltaic permeability
Figure；
Fig. 9 is the signal of the high fdrequency component of the sequence of the actual load data under different set distributed photovoltaic permeability
Figure；
Figure 10 is that the actual load of different distributions formula photovoltaic permeability predicts the schematic diagram of stability；
Figure 11 is the structural schematic diagram of distributed photovoltaic maximum permeability determining device in the embodiment of the present invention.
Specific embodiment
Understand in order to make the object, technical scheme and advantages of the embodiment of the invention clearer, with reference to the accompanying drawing to this hair
Bright embodiment is described in further details.Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
Fig. 1 is the flow chart that distributed photovoltaic maximum permeability determines method in the embodiment of the present invention, as shown in Figure 1, institute
Maximum permeability is stated to determine based on load prediction stability, this method comprises:
Step 101, according to the distributed photovoltaic active power data and actual load data of target area, target area is obtained
The original loads data in domain；
Step 102, circulation executes following steps, until the load prediction stability of target area and default load prediction are steady
Surely the difference spent within a preset range, determines that newest distributed photovoltaic permeability is distributed photovoltaic maximum permeability:
Step 1021, it according to the distributed photovoltaic active power data and original loads data of target area, is set
Actual load data under distributed photovoltaic permeability；
Step 1022, the actual load data under setting distributed photovoltaic permeability are decomposed, obtains setting distribution
The low frequency component and high fdrequency component of the sequence of actual load data under formula photovoltaic permeability；
Step 1023, according to the low frequency component and height of the sequence of the actual load data under setting distributed photovoltaic permeability
Frequency component obtains the load prediction stability of target area；
Step 1024, if the difference of the load prediction stability of target area and default load prediction stability is not default
In range, setting distributed photovoltaic permeability is adjusted, by the distributed photovoltaic permeability replacement setting distributed photovoltaic after adjusting
Permeability.
In embodiments of the present invention, in embodiments of the present invention, by setting distributed photovoltaic computing permeability target area
The load prediction stability in domain, if the difference of the load prediction stability of target area and default load prediction stability is not pre-
If in range, adjusting setting distributed photovoltaic permeability, the distributed photovoltaic permeability replacement after adjusting is set into distributed light
Volt permeability, executes repeatedly, until the load prediction stability of target area and the difference of default load prediction stability are pre-
If in range, therefore, above procedure considers the load prediction stability of target area, and finally determining distributed photovoltaic is most
Big permeability have passed through successive ignition, and precision is high, can effectively support power balance arrangement, be executed repeatedly to the load of target area
The difference of prediction stability and default load prediction stability within a preset range, therefore can effectively support load prediction.
Firstly, it is necessary to obtain the distributed photovoltaic active power data and actual load data of target area, target area
Distributed photovoltaic active power data can indicate are as follows:
Ppv(t)={ Ppvt| t=1,2,3 ... } (1)
Wherein, PpvtFor the distributed photovoltaic active power data of target area；
Ppv(t) for target area distributed photovoltaic active power data sequence.
The actual load data of target area can indicate are as follows:
Pl(t)={ Plt| t=1,2,3 ... } (2)
Wherein, PltFor the actual load data of target area；
Pl(t) for target area actual load data sequence.
It is obtained using following formula according to the distributed photovoltaic active power data and actual load data of target area
The original loads data of target area:
Pel(t)={ Pelt| t=1,2,3 ... }=Pl(t)+Ppv(t) (3)
Wherein: PeltFor the original loads data of target area, Pel(t) for target area original loads data sequence.
In one embodiment, the actual load data under setting distributed photovoltaic permeability are decomposed, is set
The low frequency component and high fdrequency component of the sequence of actual load data under distributed photovoltaic permeability may include:
Fourier decomposition is carried out to the sequence of the actual load data under setting distributed photovoltaic permeability, obtains setting point
The actual load data of frequency domain under cloth photovoltaic permeability；
From the actual load data of the frequency domain under setting distributed photovoltaic permeability, extracts setting distributed photovoltaic and seep
The low frequency component and high fdrequency component of the actual load data of frequency domain under saturating rate；
The low frequency component and high fdrequency component of the actual load data of frequency domain under setting distributed photovoltaic permeability are carried out
Inverse transformation obtains the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability.
In one embodiment, according to the low frequency component of the sequence of the actual load data under setting distributed photovoltaic permeability
And high fdrequency component, the load prediction stability of target area is obtained, may include:
According to the high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, target area is obtained
Load prediction stability the upper limit；
According to the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, obtain
Obtain the lower limit of the load prediction stability of target area；
According to the upper and lower bound of the load prediction stability of target area, the load prediction for calculating target area is stablized
Degree.
In one embodiment, if the difference of the load prediction stability of target area and default load prediction stability does not exist
In preset range, setting distributed photovoltaic permeability is adjusted, may include:
If the load prediction stability of target area is greater than default load prediction stability, increases setting distributed photovoltaic and seep
Saturating rate；
If the load prediction stability of target area is less than default load prediction stability, reduces setting distributed photovoltaic and seep
Saturating rate.
In one embodiment, following formula can be used, according to the distributed photovoltaic active power data of target area and
Original loads data obtain the actual load data under setting distributed photovoltaic permeability:
P'l(t)={ P'lt| t=0,1,2,3 ... }=Pel(t)-P'pv(t) (6)
Wherein, MpvmTo set distributed photovoltaic permeability；
P'pvIt (t) is the sequence of the distributed photovoltaic active power data under setting distributed photovoltaic permeability；
P'lIt (t) is the sequence of the actual load data under setting distributed photovoltaic permeability；
MpvFor the distributed photovoltaic permeability of target area；
Ppv(t) for target area distributed photovoltaic active power data sequence；
Pel(t) for target area original loads data sequence.
In one embodiment, the actual load data under setting distributed photovoltaic permeability are decomposed, is set
The low frequency component and high fdrequency component of the sequence of actual load data under distributed photovoltaic permeability, comprising:
Using following formula, Fourier point is carried out to the sequence of the actual load data under setting distributed photovoltaic permeability
Solution obtains the actual load data of the frequency domain under setting distributed photovoltaic permeability:
X(ωi)=N (ai-jbi) (8)
Wherein, P'lIt (t) is the sequence of the actual load data under setting distributed photovoltaic permeability；
X(ωi) it is the actual load data for setting the frequency domain under distributed photovoltaic permeability；
N is the length for setting the sequence of the actual load data under distributed photovoltaic permeability；
ωiFor angular frequency；
a0+ D (t) is component diurnal periodicity for setting the actual load data under distributed photovoltaic permeability；
W (t) is the cycle component for setting the actual load data under distributed photovoltaic permeability；
L (t) is the low frequency component for setting the actual load data under distributed photovoltaic permeability；
H (t) is the high fdrequency component for setting the actual load data under distributed photovoltaic permeability；
The period of low frequency component L (t), the precision of prediction of low frequency component was higher, and load prediction stability is got over generally higher than for 24 hours
Good, high fdrequency component H (t) reflects the randomness fluctuation of actual load data, for short-term load forecasting, can not establish effective
Model prediction, therefore high fdrequency component belongs to unpredictable component.
It extracts and sets from the actual load data of the frequency domain under setting distributed photovoltaic permeability using following formula
Determine the low frequency component and high fdrequency component of the actual load data of the frequency domain under distributed photovoltaic permeability:
Wherein: component a diurnal periodicity of the sequence of actual load data0The actual load data of the corresponding frequency domain of+D (t)
Diurnal periodicity, the angular frequency of component wasFor example, by taking 14 light, 96 point sampling as an example,
The set of all possible value compositions of angular frequency subscript i are as follows: and i | mod (i, 14)=0, i=0,1 ..., N-1 } ∪ { 0 }；
The cycle of the actual load data of the corresponding frequency domain of the cycle component W (t) of the sequence of actual load data point
The angular frequency of amount isFor example, with 14 light, 96 point sampling
For, the collection of all possible values of angular frequency subscript i be combined into i | mod (i, 2)=0 and | mod (i, 14) ≠ 0, i=0,1 ...,
N-1}；
The low frequency component of the actual load data of the corresponding frequency domain of the low frequency component L (t) of the sequence of actual load data
Angular frequency isThe collection of all possible values of angular frequency subscript i be combined into i |
Mod (i, 2) ≠ 0, i=0,1 ..., 14 }；
The high fdrequency component of the actual load data of the corresponding frequency domain of the high fdrequency component H (t) of the sequence of actual load data
Angular frequency isThe collection of all possible values of angular frequency subscript i be combined into i |
Mod (i, 2) ≠ 0, i=0,1 ..., 14 }.
When it is implemented, the low frequency component and height of the actual load data to the frequency domain under setting distributed photovoltaic permeability
Frequency component carries out inverse transformation, obtains the low frequency component and height of the sequence of the actual load data under setting distributed photovoltaic permeability
Frequency component, i.e., to X'(ωi) inverse transformation is carried out, the actual load data under setting distributed photovoltaic permeability can be respectively obtained
Sequence Dx (t), Wx (t), Lx (t), Hx (t), take the real part of Lx (t), Hx (t), low frequency component L as calculated respectively
(t) and high fdrequency component H (t).
In one embodiment, following formula can be used, according to the actual load number under setting distributed photovoltaic permeability
According to sequence high fdrequency component, obtain the upper limit of the load prediction stability of target area:
Wherein, ξupperFor the upper limit of the load prediction stability of target area；
H (t) is the high fdrequency component for setting the sequence of the actual load data under distributed photovoltaic permeability；
P'lIt (t) is the sequence of the actual load data under setting distributed photovoltaic permeability；
N is the length for setting the sequence of the actual load data under distributed photovoltaic permeability；
The low frequency component of actual load data, is mainly influenced by meteorologic factor, models in bus short-term load forecasting
In can improve prediction stability in view of low frequency component, but either the photovoltaic generating system that electric load still accesses is gone out
Power, the power by meteorologic factor affecting laws is not consistent, and the precision of weather forecast (numerical weather forecast) is not also total
Energy Accurate Prediction, therefore, low frequency component, which belongs to part, can be predicted component, can use following formula, according to setting distribution light
The low frequency component and high fdrequency component for lying prostrate the sequence of the actual load data under permeability, the load prediction for obtaining target area are stablized
The lower limit of degree:
Wherein, ξlowerFor the lower limit of the load prediction stability of target area；
L (t) is the low frequency component for setting the sequence of the actual load data under distributed photovoltaic permeability；
Following formula can be used, according to the upper and lower bound of the load prediction stability of target area, calculates target area
The load prediction stability in domain:
Wherein, ξ0For the load prediction stability of target area.
A specific embodiment is given below, illustrates that distributed photovoltaic maximum permeability proposed by the present invention determines the tool of method
Body application.
The Hebei province Langfang Prefecture access 35kV power distribution network on May 1st, 2017 to May 31 is acquired, and in place platform area
The reality of the 40MW distributed photovoltaic active power data of consumption and the 220kV Long Hezhan bus of photovoltaic generating system access
The resolution ratio of border load data, distributed photovoltaic active power and actual load data is 5 minutes.
Firstly, according to the distributed photovoltaic active power data and actual load data of Hebei province Langfang Prefecture, using public affairs
Formula (3) obtains the original loads data of this area, and Fig. 2 is the schematic diagram data of Hebei province Langfang Prefecture in the embodiment of the present invention.
Then, circulation executes following steps, until the load prediction stability of target area and default load prediction are stablized
The difference of degree within a preset range, determines that newest distributed photovoltaic permeability is distributed photovoltaic maximum permeability:
Using formula (4), formula (5) and formula (6), according to the distributed photovoltaic active power data and original of target area
Beginning load data obtains the actual load data under setting distributed photovoltaic permeability.
Using formula (7), formula (8), formula (9) and formula (10), to the reality under setting distributed photovoltaic permeability
Load data is decomposed, and the low frequency component and height of the sequence of the actual load data under setting distributed photovoltaic permeability are obtained
Frequency component, Fig. 3-Fig. 7 are each of the sequence of the actual load data in the embodiment of the present invention under setting distributed photovoltaic permeability
The schematic diagram of component, wherein Fig. 3 is each point of the sequence of actual load data when setting distributed photovoltaic permeability as 0%
The schematic diagram of amount, Fig. 4 are each component of the sequence of actual load data when setting distributed photovoltaic permeability as 12.77%
Schematic diagram, Fig. 5 are the schematic diagram of each component of the sequence of actual load data when setting distributed photovoltaic permeability as 20%,
Fig. 6 is the schematic diagram of each component of the sequence of actual load data when setting distributed photovoltaic permeability as 40%, and Fig. 7 is to set
Determine the schematic diagram of each component of the sequence of actual load data when distributed photovoltaic permeability is 60%.
Fig. 8 is the signal of the low frequency component of the sequence of the actual load data under different set distributed photovoltaic permeability
Figure, as shown in figure 8, the low frequency component of the sequence of actual load data is presented with the raising of setting distributed photovoltaic permeability
The phenomenon that unordered variation, reason are relative to conventional electric power load, and photovoltaic power generation is with the variation of scale above for 24 hours, by continuous
The influence of the irradiation level condition of day.
Fig. 9 is the signal of the high fdrequency component of the sequence of the actual load data under different set distributed photovoltaic permeability
Figure, as shown in figure 9, with the raising of setting distributed photovoltaic permeability, the high fdrequency component of the sequence of actual load data is obvious
Increase, illustrates that, relative to conventional electric power load, photovoltaic power output is influenced to become apparent by meteorology.
By Fig. 8 and Fig. 9 it is found that after distributed photovoltaic permeability increase, low frequency component and height with stochastic volatility
Frequency component increases, so that the accuracy that conventional load prediction model carries out short-term forecast to bus load reduces.
Then, using formula (11), formula (12) and formula (13), according to the reality under setting distributed photovoltaic permeability
The low frequency component and high fdrequency component of the sequence of load data obtain the load prediction stability of target area,
Table 1 is that the actual load of different distributions formula photovoltaic permeability predicts the estimated value of stability, and Figure 10 is different distributions
The schematic diagram of the actual load prediction stability of formula photovoltaic permeability, can see in conjunction with table 1 and Figure 10, with dispersion access dragon
The photovoltaic permeability of river 220kV bus gradually increases, and the upper limit of load prediction stability, lower limit gradually decline, illustrate with point
After cloth photovoltaic is injected into power distribution network as a kind of negative sense load, due to photovoltaic power generation characteristic, to random wave caused by meteorology
Dynamic property and enchancement factor significantly much compared to conventional electric power load cause using conventional load prediction model to access photovoltaic
Power generation, especially under the application scenarios of high permeability photovoltaic power generation, precision of prediction will appear apparent reduction.
The estimated value of the actual load prediction stability of 1 different distributions formula photovoltaic permeability of table
The default load prediction stability of this area is 80%, obtains this area by iterative calculation and meets the default load
The distributed photovoltaic maximum permeability for predicting stability is 42%.
In embodiments of the present invention, according to the distributed photovoltaic active power data and actual load data of target area,
Obtain the original loads data of target area；Circulation executes following steps, until the load prediction stability of target area and pre-
If the difference of load prediction stability is within a preset range, determine that newest distributed photovoltaic permeability is that distributed photovoltaic is maximum
Permeability: according to the distributed photovoltaic active power data and original loads data of target area, setting distributed photovoltaic is obtained
Actual load data under permeability；Actual load data under setting distributed photovoltaic permeability are decomposed, are set
Determine the low frequency component and high fdrequency component of the sequence of the actual load data under distributed photovoltaic permeability；According to the distributed light of setting
The low frequency component and high fdrequency component for lying prostrate the sequence of the actual load data under permeability, the load prediction for obtaining target area are stablized
Degree；If not within a preset range, adjusting is set the difference of the load prediction stability of target area and default load prediction stability
Distributed photovoltaic permeability is determined, by the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting.In this hair
In bright embodiment, by setting the load prediction stability of distributed photovoltaic computing permeability target area, if target area
The difference of load prediction stability and default load prediction stability within a preset range, does not adjust setting distributed photovoltaic infiltration
Rate executes the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting repeatedly, until target area
Within a preset range, therefore, above procedure considers target to the difference of load prediction stability and default load prediction stability
The load prediction stability in region, and finally determining distributed photovoltaic maximum permeability have passed through successive ignition, precision is high, can
Effectively support power balance arrangement is executed repeatedly to the load prediction stability of target area and default load prediction stability
Difference within a preset range, therefore can effectively support load prediction.
Based on the same inventive concept, it determines and fills the embodiment of the invention also provides a kind of distributed photovoltaic maximum permeability
It sets, as described in following implementation.Since these principles solved the problems, such as determine that method is similar to distributed photovoltaic maximum permeability,
Therefore the implementation of device may refer to the implementation of method, repeats place and is not repeating.
Figure 11 is the structural schematic diagram of distributed photovoltaic maximum permeability determining device in the embodiment of the present invention, such as Figure 11 institute
Show, which includes:
Original loads data obtaining module 1101, for the distributed photovoltaic active power data and reality according to somewhere
Load data obtains the original loads data of target area；
Distributed photovoltaic maximum permeability determining module 1102 executes following steps for recycling, until target area
The difference of load prediction stability and default load prediction stability within a preset range, determines newest distributed photovoltaic infiltration
Rate is distributed photovoltaic maximum permeability:
According to the distributed photovoltaic active power data and original loads data of target area, setting distributed photovoltaic is obtained
Actual load data under permeability；
Actual load data under setting distributed photovoltaic permeability are decomposed, setting distributed photovoltaic infiltration is obtained
The low frequency component and high fdrequency component of the sequence of actual load data under rate；
According to the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, obtain
Obtain the load prediction stability of target area；
If the difference of the load prediction stability of target area and default load prediction stability within a preset range, is not adjusted
Section setting distributed photovoltaic permeability, by the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting.
In conclusion in embodiments of the present invention, according to the distributed photovoltaic active power data and reality of target area
Load data obtains the original loads data of target area；Circulation executes following steps, until the load prediction of target area is steady
The difference of fixed degree and default load prediction stability within a preset range, determines newest distributed photovoltaic permeability for distribution
Photovoltaic maximum permeability: according to the distributed photovoltaic active power data and original loads data of target area, setting point is obtained
Actual load data under cloth photovoltaic permeability；Actual load data under setting distributed photovoltaic permeability are divided
Solution obtains the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability；According to setting
The low frequency component and high fdrequency component for determining the sequence of the actual load data under distributed photovoltaic permeability obtain the negative of target area
Lotus predicts stability；If the difference of the load prediction stability of target area and default load prediction stability is not in preset range
It is interior, setting distributed photovoltaic permeability is adjusted, by the distributed photovoltaic permeability replacement setting distributed photovoltaic infiltration after adjusting
Rate.In embodiments of the present invention, by setting the load prediction stability of distributed photovoltaic computing permeability target area, if mesh
The difference of the load prediction stability and default load prediction stability of marking region within a preset range, it is distributed not adjust setting
Photovoltaic permeability executes the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting, until mesh repeatedly
The difference of the load prediction stability and default load prediction stability of marking region within a preset range, therefore, examine by above procedure
Consider the load prediction stability of target area, and finally determining distributed photovoltaic maximum permeability have passed through successive ignition,
Precision is high, can effectively support power balance arrangement, execute repeatedly pre- to the load prediction stability of target area and default load
The difference of survey stability within a preset range, therefore can effectively support load prediction.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the present invention, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection scope of invention.
Claims (10)
1. a kind of distributed photovoltaic maximum permeability determines method, which is characterized in that the maximum permeability is based on load prediction
Stability determines, comprising:
According to the distributed photovoltaic active power data and actual load data of target area, the original loads of target area are obtained
Data；
Circulation executes following steps, until the load prediction stability of target area and the difference of default load prediction stability exist
In preset range, determine that newest distributed photovoltaic permeability is distributed photovoltaic maximum permeability:
According to the distributed photovoltaic active power data and original loads data of target area, setting distributed photovoltaic infiltration is obtained
Actual load data under rate；
Actual load data under setting distributed photovoltaic permeability are decomposed, are obtained under setting distributed photovoltaic permeability
Actual load data sequence low frequency component and high fdrequency component；
According to the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, mesh is obtained
Mark the load prediction stability in region；
If not within a preset range, adjusting is set the difference of the load prediction stability of target area and default load prediction stability
Distributed photovoltaic permeability is determined, by the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting.
2. distributed photovoltaic maximum permeability as described in claim 1 determines method, which is characterized in that the distributed light of setting
Actual load data under volt permeability are decomposed, and the sequence of the actual load data under setting distributed photovoltaic permeability is obtained
The low frequency component and high fdrequency component of column, comprising:
Fourier decomposition is carried out to the sequence of the actual load data under setting distributed photovoltaic permeability, it is distributed to obtain setting
The actual load data of frequency domain under photovoltaic permeability；
From the actual load data of the frequency domain under setting distributed photovoltaic permeability, setting distributed photovoltaic permeability is extracted
Under frequency domain actual load data low frequency component and high fdrequency component；
Inversion is carried out to the low frequency component and high fdrequency component of the actual load data of the frequency domain under setting distributed photovoltaic permeability
It changes, obtains the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability.
3. distributed photovoltaic maximum permeability as described in claim 1 determines method, which is characterized in that distributed according to setting
The low frequency component and high fdrequency component of the sequence of actual load data under photovoltaic permeability, the load prediction for obtaining target area are steady
Fixed degree, comprising:
According to the high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, the negative of target area is obtained
The upper limit of lotus prediction stability；
According to the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, mesh is obtained
Mark the lower limit of the load prediction stability in region；
According to the upper and lower bound of the load prediction stability of target area, the load prediction stability of target area is calculated.
4. distributed photovoltaic maximum permeability as described in claim 1 determines method, which is characterized in that if target area is negative
The difference of lotus prediction stability and default load prediction stability within a preset range, does not adjust setting distributed photovoltaic infiltration
Rate, comprising:
If the load prediction stability of target area is greater than default load prediction stability, increase setting distributed photovoltaic infiltration
Rate；
If the load prediction stability of target area is less than default load prediction stability, reduce setting distributed photovoltaic infiltration
Rate.
5. distributed photovoltaic maximum permeability as described in claim 1 determines method, which is characterized in that following formula is used,
According to the distributed photovoltaic active power data and original loads data of target area, obtain under setting distributed photovoltaic permeability
Actual load data:
P′l(t)={ P 'lt| t=0,1,2,3 ... }=Pel(t)-P′pv(t)
Wherein, MpvmTo set distributed photovoltaic permeability；
P′pvIt (t) is the sequence of the distributed photovoltaic active power data under setting distributed photovoltaic permeability；
P′lIt (t) is the sequence of the actual load data under setting distributed photovoltaic permeability；
MpvFor the distributed photovoltaic permeability of target area；
Ppv(t) for target area distributed photovoltaic active power data sequence；
Pel(t) for target area original loads data sequence.
6. distributed photovoltaic maximum permeability as claimed in claim 2 determines method, which is characterized in that the distributed light of setting
Actual load data under volt permeability are decomposed, and the sequence of the actual load data under setting distributed photovoltaic permeability is obtained
The low frequency component and high fdrequency component of column, comprising:
Using following formula, Fourier decomposition is carried out to the sequence of the actual load data under setting distributed photovoltaic permeability,
Obtain the actual load data of the frequency domain under setting distributed photovoltaic permeability:
X(ωi)=N (ai-jbi)
Wherein, P 'lIt (t) is the sequence of the actual load data under setting distributed photovoltaic permeability；
X(ωi) it is the actual load data for setting the frequency domain under distributed photovoltaic permeability；
N is the length for setting the sequence of the actual load data under distributed photovoltaic permeability；
ωiFor angular frequency；
a0+ D (t) is component diurnal periodicity for setting the sequence of the actual load data under distributed photovoltaic permeability；
W (t) is the cycle component for setting the sequence of the actual load data under distributed photovoltaic permeability；
L (t) is the low frequency component for setting the sequence of the actual load data under distributed photovoltaic permeability；
H (t) is the high fdrequency component for setting the sequence of the actual load data under distributed photovoltaic permeability；
Setting point is extracted from the actual load data of the frequency domain under setting distributed photovoltaic permeability using following formula
The low frequency component and high fdrequency component of the actual load data of frequency domain under cloth photovoltaic permeability:
Wherein: component a diurnal periodicity of the sequence of actual load data0The diurnal periodicity of the actual load data of the corresponding frequency domain of+D (t)
The angular frequency of component is
The cycle component of the actual load data of the corresponding frequency domain of the cycle component W (t) of the sequence of actual load data
Angular frequency is
The angle of the low frequency component of the actual load data of the corresponding frequency domain of the low frequency component L (t) of the sequence of actual load data
Frequency is
The angle of the high fdrequency component of the actual load data of the corresponding frequency domain of the high fdrequency component H (t) of the sequence of actual load data
Frequency is
7. distributed photovoltaic maximum permeability as claimed in claim 3 determines method, which is characterized in that following formula is used,
According to the high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, the load for obtaining target area is pre-
Survey the upper limit of stability:
Wherein, ξupperFor the upper limit of the load prediction stability of target area；
H (t) is the high fdrequency component for setting the sequence of the actual load data under distributed photovoltaic permeability；
P′lIt (t) is the sequence of the actual load data under setting distributed photovoltaic permeability；
N is the length for setting the sequence of the actual load data under distributed photovoltaic permeability；
Using following formula, according to the low frequency component and height of the sequence of the actual load data under setting distributed photovoltaic permeability
Frequency component obtains the lower limit of the load prediction stability of target area:
Wherein, ξlowerFor the lower limit of the load prediction stability of target area；
L (t) is the low frequency component for setting the sequence of the actual load data under distributed photovoltaic permeability；
The load of target area is calculated according to the upper and lower bound of the load prediction stability of target area using following formula
Predict stability:
Wherein, ξ0For the load prediction stability of target area.
8. a kind of distributed photovoltaic maximum permeability determining device, which is characterized in that the maximum permeability is based on load prediction
Stability determines, comprising:
Original loads data obtaining module, for the distributed photovoltaic active power data and actual load number according to somewhere
According to obtaining the original loads data of target area；
Distributed photovoltaic maximum permeability determining module executes following steps for recycling, until the load prediction of target area
The difference of stability and default load prediction stability within a preset range, determines newest distributed photovoltaic permeability for distribution
Formula photovoltaic maximum permeability:
According to the distributed photovoltaic active power data and original loads data of target area, setting distributed photovoltaic infiltration is obtained
Actual load data under rate；
Actual load data under setting distributed photovoltaic permeability are decomposed, are obtained under setting distributed photovoltaic permeability
Actual load data sequence low frequency component and high fdrequency component；
According to the low frequency component and high fdrequency component of the sequence of the actual load data under setting distributed photovoltaic permeability, mesh is obtained
Mark the load prediction stability in region；
If not within a preset range, adjusting is set the difference of the load prediction stability of target area and default load prediction stability
Distributed photovoltaic permeability is determined, by the distributed photovoltaic permeability replacement setting distributed photovoltaic permeability after adjusting.
9. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor
Calculation machine program, which is characterized in that the processor realizes any side of claim 1 to 7 when executing the computer program
Method.
10. a kind of computer readable storage medium, which is characterized in that the computer-readable recording medium storage has perform claim
It is required that the computer program of 1 to 7 any the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910198323.9A CN109978717A (en) | 2019-03-15 | 2019-03-15 | Distributed photovoltaic maximum permeability determines method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910198323.9A CN109978717A (en) | 2019-03-15 | 2019-03-15 | Distributed photovoltaic maximum permeability determines method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109978717A true CN109978717A (en) | 2019-07-05 |
Family
ID=67079059
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910198323.9A Pending CN109978717A (en) | 2019-03-15 | 2019-03-15 | Distributed photovoltaic maximum permeability determines method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109978717A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103854067A (en) * | 2012-12-06 | 2014-06-11 | 国家电网公司 | Method and device for detecting power load forecasting precision |
| CN106600173A (en) * | 2017-01-25 | 2017-04-26 | 北京中电普华信息技术有限公司 | Distributed photovoltaic power load permeability estimating method and distributed photovoltaic power load permeability estimating device |
| WO2017155421A1 (en) * | 2016-03-07 | 2017-09-14 | Centro De Investigação Em Energia Ren - State Grid, S.A | Method and system for forecasting the power output of a group of photovoltaic power plants and managing the integration of said power output into a power grid |
| CN107482676A (en) * | 2017-08-11 | 2017-12-15 | 中国电力科学研究院 | Consider the determination method and apparatus of the distributed power source maximum permeability of fluctuation |
-
2019
- 2019-03-15 CN CN201910198323.9A patent/CN109978717A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103854067A (en) * | 2012-12-06 | 2014-06-11 | 国家电网公司 | Method and device for detecting power load forecasting precision |
| WO2017155421A1 (en) * | 2016-03-07 | 2017-09-14 | Centro De Investigação Em Energia Ren - State Grid, S.A | Method and system for forecasting the power output of a group of photovoltaic power plants and managing the integration of said power output into a power grid |
| CN106600173A (en) * | 2017-01-25 | 2017-04-26 | 北京中电普华信息技术有限公司 | Distributed photovoltaic power load permeability estimating method and distributed photovoltaic power load permeability estimating device |
| CN107482676A (en) * | 2017-08-11 | 2017-12-15 | 中国电力科学研究院 | Consider the determination method and apparatus of the distributed power source maximum permeability of fluctuation |
Non-Patent Citations (1)
| Title |
|---|
| 黄巍: "电能质量约束下主动配电网光伏最大渗透率研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105184391B (en) | Wind farm wind velocity and power forecasting method based on wavelet decomposition and support vector machines | |
| Al-Yahyai et al. | Nested ensemble NWP approach for wind energy assessment | |
| CN110659448B (en) | Non-orthogonal ensemble prediction initial value disturbance algorithm | |
| CN108985521A (en) | Output power predicting method, device, equipment and the medium of photovoltaic generating system | |
| Cai et al. | Bred vectors of the Zebiak–Cane model and their potential application to ENSO predictions | |
| Liu et al. | Application of particle swarm optimization and extreme learning machine forecasting models for regional groundwater depth using nonlinear prediction models as preprocessor | |
| Ye et al. | Combined approach for short-term wind power forecasting based on wave division and Seq2Seq model using deep learning | |
| Kaplan et al. | A novel method based on Weibull distribution for short-term wind speed prediction | |
| Zhou et al. | A hybrid grey prediction model for small oscillation sequence based on information decomposition | |
| CN103345585A (en) | Wind power prediction correction method and system based on support vector machine | |
| CN105846425A (en) | Economic dispatching method based on general wind power forecasting error model | |
| Cao et al. | Physics-informed graphical representation-enabled deep reinforcement learning for robust distribution system voltage control | |
| CN111832800A (en) | Photovoltaic power station power prediction method and device | |
| Ma et al. | Optimal economic operation of islanded microgrid by using a modified pso algorithm | |
| Sun et al. | Forecasting Carbon Dioxide Emissions in China Using Optimization Grey Model. | |
| CN111651896A (en) | Wind power plant flow field calculation method based on actual wind speed and actual complex terrain | |
| Kim et al. | A probabilistic approach to potential estimation of renewable energy resources based on augmented spatial interpolation | |
| CN109978277A (en) | Region online load prediction technique and device in photovoltaic power generation | |
| CN113191091A (en) | Wind speed prediction method, system and equipment based on hybrid deep learning mechanism | |
| CN117394310A (en) | Method and system for calculating area light Fu Dengxiao light resources affected by multiple weather types | |
| CN109978717A (en) | Distributed photovoltaic maximum permeability determines method and device | |
| Kuznetsova et al. | High resolution waves and weather forecasts using adapted WAVEWATCH III and WRF models | |
| Yang et al. | Improved nonlinear mapping network for wind power forecasting in renewable energy power system dispatch | |
| CN113779861B (en) | Photovoltaic Power Prediction Method and Terminal Equipment | |
| Dong et al. | Very short-term photovoltaic power forecasting using uncertain basis function |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190705 |
|
| RJ01 | Rejection of invention patent application after publication |