CN103544362A - Harmonic medium and long term prediction method based on two-dimensional curve prediction - Google Patents
Harmonic medium and long term prediction method based on two-dimensional curve prediction Download PDFInfo
- Publication number
- CN103544362A CN103544362A CN201310539260.1A CN201310539260A CN103544362A CN 103544362 A CN103544362 A CN 103544362A CN 201310539260 A CN201310539260 A CN 201310539260A CN 103544362 A CN103544362 A CN 103544362A
- Authority
- CN
- China
- Prior art keywords
- sequence
- year
- prediction
- predicted
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a harmonic medium and long term prediction method based on two-dimensional curve prediction. The method includes the following steps: 1) acquiring historical electric energy quality index data of at least one monitoring point in a set time period in a prediction interval and conducting data preprocessing to obtain a historical data series; 2) judging whether a year to be predicted has an innovation value, executing the step 3) on yes judgment, adjusting the historical data series on no judgment and executing the step 3) with the adjusted historical data series; 3) conducting time series curve prediction on a monthly development series of monthly quantity according to the historical data series to obtain a transverse prediction series of the year to be predicted; 4) conducting time series curve prediction on a yearly development series of the monthly quantity according to the historical data series to obtain a vertical prediction series of the year to be predicted; 5) conducting weighted average calculation to obtain the two-dimensional prediction series of the year to be predicted. Compared with the prior art, the method has the advantages of being high in prediction accuracy, anti-jamming, simple in principle and the like.
Description
Technical field
The present invention relates to a kind of electrical network quality of power supply Forecasting Methodology, especially relate to a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction.
Background technology
The quality of power supply is being related to safety and stability, the economical operation of electrical network, along with the transformation of modern power network to high reliability and quality supply, improve the electrical network quality of power supply become guarantee power system safety and stability operation in the urgent need to.
The electrical network quality of power supply on the one hand because of non-linear, large capacity and the strong load equipment of impact constantly access aggravate, on the other hand because its control measures and the continuous extension to increase capacity of electrical network improve.Also there are some uncertain factors, as the automatic protection of spontaneous phenomenon and power equipment and device and the change of normal operating mode etc. have all caused serious interference to the trend of the quality of power supply simultaneously.
Power transmission network, as the chief component of electric system, is being born jumbo electric power transfer task, and its power quality problem influence area is large, therefore most important to the safety of whole electrical network.If can the quality of power supply be predicted, analyze its development trend, find that early the quality of power supply is maybe by the problem worsening, take corresponding measure to be improved and administer, thereby the loss causing thus is even avoided in reduction, for network optimization provides strong foundation and decision support, there is important theory value and realistic meaning.Therefore the quality of power supply is predicted to be very necessary, also more and more higher to the requirement of quality of power supply prediction, medium-and long-term forecasting is also had to urgent demand.Medium-term forecast refers in 1~2 year, month or the prediction in season, long-term forecasting refers to moon, the season of 1~10 year, the prediction in year.
Summary of the invention
Object of the present invention is exactly to provide in order to overcome the defect of above-mentioned prior art existence a kind of harmonic wave medium-and long-term forecasting method that precision of prediction is high, anti-interference, principle is simply predicted based on two-dimensional curve.
Object of the present invention can be achieved through the following technical solutions:
A harmonic wave medium-and long-term forecasting method for two-dimensional curve prediction, the method comprises the following steps:
1) obtain the historical power quality index data of at least one monitoring point in setting-up time section in forecast interval, the line number of going forward side by side Data preprocess, obtains historical data sequence;
2) whether have new breath value, whether have known electric energy quality index data, if so, perform step 3 if judging to be predicted year if judging to be predicted year), if not, historical data sequence is adjusted, and with the historical data sequence execution step 3 after adjusting);
3) according to historical data sequence, to the moon, the monthly developmental sequence of tolerance carries out timing curve prediction, obtains the lateral prediction sequence of to be predicted year
4) according to historical data sequence, to the moon, the annual developmental sequence of tolerance carries out timing curve prediction, obtains longitudinal forecasting sequence of to be predicted year
5) lateral prediction sequence and longitudinal forecasting sequence are weighted on average, obtain the two-dimensional prediction sequence of to be predicted year
and add in historical data sequence:
In formula, w
1, w
2the weight of horizontal, the longitudinal predicted value of difference.
Described setting-up time section is greater than 1 year.
Described data pre-service is specially:
The power quality index data that each monitoring point is obtained are weighted the comprehensive historical data sequence that on average obtains forecast interval, and computing formula is as follows:
Y
region=wy
point
w
i=p
i/p
sum,i=1,2,...,N
Wherein, w
ifor i the element of weight vectors w, represent the weight of i monitoring point, p
ibe the power of i monitoring point, p
sumfor the general power of each monitoring point in forecast interval, Y
tegionfor the comprehensive historical data sequence of forecast interval, y
pointpower quality index matrix for each monitoring point in forecast interval.
Described step 3) be specially:
301) characteristic parameter of the horizontal curve in year to be predicted is predicted, described characteristic parameter comprises average perunit value ρ and minimum perunit value β;
302) take on the time close to time of to be predicted year be standard year, according to monthly perunit value sequence and the step 301 of standard year) characteristic parameter recording in advance obtains the monthly perunit value sequence of to be predicted year;
303) according to the monthly perunit value sequence of to be predicted year and the new breath value of to be predicted year, carry out newly ceasing famousization processing, obtain lateral prediction value.
Described step 301), in, the characteristic parameter Forecasting Methodology of employing comprises the moving method of average, regression analysis or exponential smoothing.
Described step 302) comprising:
Will
after descending sequence, become sequence
if the monthly perunit value sequence d of to be predicted year
tsequence postscript is
before and after sequence, under corresponding original of sequence, be labeled as h
j, per days, number scale was T, had following relation:
X
t, iwith the relational expression of characteristic parameter ρ, β be:
B) set up mathematical model:
Wherein,
C) to step b) in mathematical model carry out iterative, iterative process is specially:
C1) introduce Lagrange multiplier w
t=[w
1, w
2..., w
t-1] and v
t=[v
1, v
2], and remember W
0=diag{w
i, W
0e=w, e
t=[1,1 ..., 1]; Put initial value W
0=0, iterations q=1, given condition of convergence ε (ε > 0);
C2) calculate v:v=(AA
t)
-1[b-A (x
(0)+ W
0e)];
C3) calculate x
(*)=x
(0)+ W
0e+A
tv, judgement x
(*)in each component
if
put w
i=0; Otherwise, order
obtain thus x
(*), W
0;
C4) judgement || Ax
(*)||
2/ || b||
2whether < ε sets up, and if so, stops iteration, obtains optimum solution x
(*); If not, put q=q+1, return to step c2);
D) to step c) the optimum solution x that obtains
(*)carrying out contrary generation processes:
Described step 303), in, newly cease famousization processing and be specially:
If the power quality index data of the front m of to be predicted year month are known, the new breath value sequence of to be predicted year is { y
t, 1, y
t, 2, y
t, m, the lateral prediction value in to be predicted year residue month is
Wherein, d
t, j,d
t, kbe respectively the perunit value of the to be predicted year j month, the k month, y
t, kmeet { y
t, k| min (v
(k)), v
(k)computing formula as follows
Described step 4), in, when known annual developmental sequence is greater than three values, longitudinally predicted value is identical with the computing method of lateral prediction value.
Described step 4) in, when known annual developmental sequence only has two values, adopt growth ratio method to calculate longitudinal predicted value, its formula is as follows:
Wherein, k=1,2 ..., m, j=m+1, m+2 ..., T.
Described step 5) in, the weight w of horizontal, longitudinal predicted value
1, w
2meet
Compared with prior art, the present invention has the following advantages:
1. the inventive method has adopted the large value of moon tolerance 95% probability of Detecting Power Harmonics point to predict, has weakened thus harmonic wave random factor to the harmful effect predicting the outcome.
2. the inventive method adopts the interval integrated value of harmonic wave to predict, has weakened the impact that each monitoring point random variation produces.
3. the inventive method has considered harmonic wave moon tolerance the same year (monthly development trend is also horizontal trend) and trend that the same period (annual development trend, i.e. longitudinal development trend), two dimension developed over the years month by month, has used for reference the method for load prediction.Because measuring two-way trend the moon, formed the netted development relation in its space, within each month, tolerance is on the point of crossing of the netted comprehensive development trend in space, therefore take into account the two while predicting, takes full advantage of its natural law.Month tolerance annual developmental sequence point between be spaced apart 1 year, embodied the development and change rule under its overall background improving constantly at social development levels; And interval between its monthly developmental sequence point is 1 month, has embodied its rule with seasonal variations.The present invention adopts the method for curve prediction to utilize respectively these two kinds of rules to predict it, according to the weight of asking for, Two-way measured value is weighted and on average obtains two-dimensional prediction result.
4. the inventive method can accurately provide the development trend of the quality of power supply, and principle is simple.
Accompanying drawing explanation
Fig. 1 is the principle schematic of correction sequence of the present invention;
Fig. 2 is schematic flow sheet of the present invention;
Fig. 3 is per-unit curve model solution process schematic diagram of the present invention;
Fig. 4 is the inventive method and error comparison diagram horizontal, longitudinally independent Forecasting Methodology.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.The present embodiment be take technical solution of the present invention and is implemented as prerequisite, provided detailed embodiment and concrete operating process, but protection scope of the present invention is not limited to following embodiment.
A harmonic wave medium-and long-term forecasting method for two-dimensional curve prediction, for harmonic wave control, the deterioration of the prevention quality of power supply provide strong foundation.The method take the Monitoring Data of harmonic wave as basis, take into account its development trend same period over the years (annual development trend) and the same year month by month development trend (monthly development trend) carry out two-dimensional prediction.First for harmonic wave year developmental sequence and monthly developmental sequence predict respectively, then according to Least Square Theory, ask for weight, finally respectively annual trend and monthly trend prediction value are weighted and on average obtain two-dimensional prediction result.
The method is theoretical based on normal distribution model.Least square method is to determine the common method of forecast model function expression unknown parameter.Its model is:
aw
T=1 (2)
1≥w≥0 (3)
Wherein, z, a, Y is known quantity,
w is unknown quantity.This model is constrained linear least-squares problem.In the situation that data volume is less, can adopt process of iteration.
While predicting, if there is no known power quality index data in year to be predicted, original series does not meet predicted condition, needs to adjust.For famousization of the new breath of the inventive method prediction step, proposed the concept of correction sequence, adjusted original series and make it meet the realization that Forecasting Methodology requires to facilitate this step.
Known array is y
j=[y
i, 1y
i, 2y
i, T], i=1,2 ..., t-1 (t>=2).Forecasting sequence is: y
t=[y
t, 1y
t, 2 ... y
t, T], y wherein
t, k(k=1,2 ..., m) known.Initial month of original series is respectively: January and Dec, the i.e. natural division in year.When known 2 years and above data (being t >=2) and while within to be predicted year, there is no new breath value (being m=0), need first adjust sequence, redefine start-stop month, make m ≠ 0, thereby obtain a new sequence, be called correction sequence, its schematic diagram is as shown in Figure 1.
The key step of the inventive method is the foundation of transverse and longitudinal two-dimensional curve forecast model and solves and the asking for of weight.Curve prediction model mainly adopts the foundation of existing Day Load Curve Forecasting model and method for solving to obtain, then carry out famousization as the case may be.Because existing, the index amount of the quality of power supply take the feature that year as the cycle changes, therefore required minimum data should be the power quality index value of a year.The following forecasting process of first introducing known 1 year above historical data.
As shown in Figure 2, the harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction of the present invention comprises the following steps:
1) obtain the historical power quality index data of at least one monitoring point in setting-up time section (being greater than a year) in forecast interval, the line number of going forward side by side Data preprocess, obtains historical data sequence;
2) whether have new breath value, if so, perform step 3 if judging to be predicted year), if not, historical data sequence is adjusted, and with the historical data sequence execution step 3 after adjusting);
3) according to historical data sequence, to the moon, the monthly developmental sequence of tolerance carries out timing curve prediction, obtains the lateral prediction sequence of to be predicted year
4) according to historical data sequence, to the moon, the annual developmental sequence of tolerance carries out timing curve prediction, obtains longitudinal forecasting sequence of to be predicted year
5) lateral prediction sequence and longitudinal forecasting sequence are weighted on average, obtain the two-dimensional prediction sequence of to be predicted year
and add in historical data sequence:
In formula, w
1, w
2the weight of horizontal, the longitudinal predicted value of difference.
Below this method is specifically introduced.
1 transverse and longitudinal curve prediction
To the moon, the monthly developmental sequence of tolerance carries out timing curve prediction, the overall trend of sequence is predicted, is called lateral prediction.To the moon, the annual developmental sequence of tolerance carries out timing curve prediction, is called longitudinal prediction.Transverse and longitudinal prediction all adopts the method for curve prediction.Mainly the process with horizontal curve prediction is introduced.Because existing, the index amount of harmonic wave take the feature that year as the cycle changes, therefore can think that its per-unit curve each year is basically identical.
1.1 data pre-service
If carry out interval comprehensive electric energy quality prediction, need to carry out to raw data the pre-service of following steps.Single monitoring point does not need to carry out this step while predicting.
Because each electric energy quality monitoring point is proportionate on the ratio of interval impact and its power and general power.Therefore account for the ratio of general power is weighted each data of monitoring point of forecast interval on average to obtain interval comprehensive electric energy quality curve according to the power of monitoring point.The power quality index data that each monitoring point is obtained are weighted the comprehensive historical data sequence that on average obtains forecast interval, and computing formula is as follows:
Y
region=wy
point (4)
w
i=p
i/p
sum,i=1,2,...,N (5)
Wherein, w
ifor i the element of weight vectors w, represent the weight of i monitoring point, p
ibe the power of i monitoring point, p
sumfor the general power of each monitoring point in forecast interval, Y
regionfor the comprehensive historical data sequence of forecast interval, y
pointpower quality index matrix for each monitoring point in forecast interval.Interval censored data is compared with one point data, and enchancement factor is weakened, regular enhancing.With interval comprehensive electric energy quality data, predict the accuracy that can improve prediction.
1.2 curve prediction
Curve prediction is divided into two steps: curvilinear characteristic parameter prediction and per-unit curve prediction.Its lateral prediction result is a month horizontal per-unit curve for tolerance.
(1) prediction of curvilinear characteristic parameter
The prediction of curvilinear characteristic parameter can adopt the moving method of average, regretional analysis, exponential smoothing etc.Here adopt the moving method of average to obtain the year average perunit value ρ of characteristic parameter to be predicted and minimum perunit value β (0 < β < ρ < 1).
Make T=12 (being annual moon number), first with the annual laterally maximal value y of sequence
0, ito y
i(i=1,2 ..., t-1) carry out standardization, obtain corresponding year per-unit curve d
i(i=1,2 ..., t-1), there is following relation:
d
t,j=y
i,j/y
0,i,j=1,2,…,T (8)
Feature and shape that the average perunit value ρ of main characteristic parameters of monthly developmental sequence per-unit curve and minimum perunit value β can reflect curve, it changes the variation that has substantially reflected monthly developmental sequence curve.The curve of known t-1, its average perunit value and minimum perunit value are respectively ρ
iand β
i(i=1,2 ..., t-1).
The prediction of curvilinear characteristic parameter can adopt the moving method of average, regretional analysis, exponential smoothing etc.Here adopt the moving method of average.The characteristic parameter of predicting the monthly developmental sequence curve of to be predicted year is respectively:
(2) per-unit curve prediction
First per-unit curve prediction will determine datum curve.Can select historical each year curve to make generalized analysis, as weighted comprehensive (near big and far smaller principle), determine and represent curve.Also can select certain year actual curve with typicalness as datum curve.On the employing time of the present invention close to certain year interval comprehensive harmonic curve of to be predicted year as datum curve.Foundation and method for solving according to existing Day Load Curve Forecasting model obtain the quality of power supply per-unit curve of to be predicted year.
Suppose that monthly developmental sequence curve post youngest of known reference year value sequence is
with to be predicted year characteristic parameter ρ, in the situation of β (0 < β < ρ < 1) (being obtained by upper step prediction), carry out the prediction of this year curve.Suppose to be predicted year curve post the one value sequence d
t, j with there is similar shape.
1. Raw Data Generation is processed
In order to weaken the randomness of raw data, and provide intermediate information for setting up mathematical model, introduce the thought of gray system, to raw data
generate processing:
A. sequence is processed
Will
after descending sequence, become sequence
if the monthly perunit value sequence d of to be predicted year
isequence postscript is
before and after sequence, under corresponding original of sequence, be labeled as h
j, per days, number scale was T, T=12, has following relation:
B. difference is processed
X
t, iwith the relational expression of characteristic parameter ρ, β be:
2. mathematical model
By generation, process, problem is converted into and makes sequence x
twith
difference as far as possible little, mathematical model is:
x
t,i≥0,i=1,2,…,T-1 (26)
Order
The matrix form of problem is:
s.t.Ax=b (28)
x≥0 (29)
3. model solution
This model is a typical quadratic programming problem, can adopt the method for solving of quadratic programming.Given this problem has following feature: the gloomy matrix in sea of objective function is unit matrix, and equality constraint is linear restriction.Succinct method for solving is as follows.
Introduce Lagrange multiplier w
t=[w
1, w
2..., w
t-1] and v
t=[v
1, v
2], and remember W
0=diag{w
i, make e
t=[1,1 ..., 1],
W
0e=w (30)
Set up following Lagrangian function:
Quadratic programming is as the special case of convex programming, and K-T condition, as sufficient and necessary condition, can be expressed as, at optimum point x
(*)place:
x
(*)-x
(0)-W
0e-A
Tv=0 (32)
Ax
(*)-b=0 (33)
W
0x
(*)=0 (34)
x
(*)≥0,W
0≥0 (35)
By above-mentioned, obtained:
v=(AA
T)
-1□[b-A(x
(0)+W
0e)] (36)
(AA wherein
t)
-1for constant matrices.
As shown in Figure 3, the iterative process of model is:
A. put initial value W
0=0, iterations q=1, given condition of convergence ε (ε > 0);
B. according to formula (36), calculate v;
C. calculate x
(*)=x
(0)+ W
0e+A
tv, judgement x
(*)in each component
if
put w
i=0; Otherwise, order
obtain thus x
(*)w
0;
D. whether judgement (33) formula is set up, and is converted into the judgement condition of convergence || Ax
(*)||
2/ || b||
2whether < ε sets up, and if so, stops iteration, obtains optimum solution x
(*); If not, put q=q+1, return to step a.
4. contrary generation of result processed
A. unfavourable balance number is processed
B. contrary sequence is processed
1.3 new famousization of breath
According to the value in the known month (1~m month) of to be predicted year, new breath value is carried out the famous value that minimum variance estimate obtains prediction curve, is called famousization of information.Sequence { the y that the known month in time to be predicted is worth
t, 1, y
t, 2..., y
t, m,
Wherein, d
t,
j, d
t, kbe respectively the t j month, the perunit value of the k month, y
t, kmeet { y
t, k| min (v
(k)), v
(k)for take k month numerical value m before famousization of benchmark curve mean square deviation of individual month, by formula (42), calculated.
The lateral prediction value of t m+1~T month is
this method takes full advantage of new breath value.
The data that obtain by virtual year method by the value prediction of whole year calendar year not out, if desired can according to the estimation of prediction in virtual year newly breath value re-start prediction and famousization of calendar year per-unit curve, draw annual value.
In next the method that can not adopt of situation that only known a year and a day, data were predicted the data of Second Year to famousization of prediction curve.Need to predict and make famousization of prediction curve the maximal value in prediction year.Can rule of thumb determine maximal value, prediction accuracy in the case reduces.
While only knowing the data of a year, longitudinally prediction cannot realize.When if known longitudinal sequence only has two values, because the longitudinal iterative process of model solution does not restrain, should not adopt curve prediction model.Now adopt growth ratio method, its formula is as follows:
Wherein, k=1,2 ..., m, j=m+1, m+2 ..., T.
2 two-dimensional predictions
Given data is carried out two-dimensional prediction by weighted mean during more than 1 year, makes full use of its month to measure horizontal trend and longitudinal trend correlativity separately.
The criterion that weight is chosen is to make predicated error meet least-squares estimation.Make two-dimensional estimation value be
, weight is ww=(w
1, w
2), have:
Wherein, i=1,2 ..., m, the actual value of i month and the error of its estimated value are:
Criterion of least squares is exactly to wish that required weight can make the quadratic sum of evaluated error reach minimum, even performance index:
Reach minimum weight, wherein, z
ibe the actual value of i month,
for sequence
in i element,
for sequence
in i element.According to the process of iteration of Constrained least square method, solve.Given initial value: w (0)=[0 1], select suitable iteration step length to carry out iterative.
Therefore the estimated value for the t m+1 month is:
When only knowing the data of a year, longitudinally prediction cannot be carried out, therefore directly using lateral prediction result as two-dimensional prediction result.
Below with instantiation explanation the inventive method.
The interval integrated voltage total percent harmonic distortion (VTHD) of four monitoring point 2009-2011s and the data of total harmonic current of table 1 for calculating according to formula (4), its weight is as shown in table 2.What the inventive method adopted is the large value of monthly data 95% probability, to weaken the impact of quality of power supply undulatory property.Predicting the outcome of correspondingly obtaining is also 95% probable value.
Integrated data between table 1 monitoring occupied area
Each monitoring point weight of table 2
The known No. 2 monitoring points value of 2009 and the total percent harmonic distortion data of voltage of 4 months before 2010, predict the value in residue month in 2010.The result of prediction is as shown in table 3, and error as shown in Figure 4.
Table 3 predicts the outcome
In table 4, list absolute average and the variance of each method predicated error, evaluated respectively its prediction accuracy and stability.From predicting the outcome and Error Graph and table 4, the result of two-dimensional prediction is with respect to laterally or longitudinally predicting the outcome separately, and its error is less than both, and the poor minimum of mean square of error illustrates that two-dimensional prediction accuracy and stability are better than single directional prediction.
Table 4 comparison that predicts the outcome
Claims (10)
1. the harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction, is characterized in that, the method comprises the following steps:
1) obtain the historical power quality index data of at least one monitoring point in setting-up time section in forecast interval, the line number of going forward side by side Data preprocess, obtains historical data sequence;
2) whether have new breath value, whether have known electric energy quality index data, if so, perform step 3 if judging to be predicted year if judging to be predicted year), if not, historical data sequence is adjusted, and with the historical data sequence execution step 3 after adjusting);
3) according to historical data sequence, to the moon, the monthly developmental sequence of tolerance carries out timing curve prediction, obtains the lateral prediction sequence of to be predicted year
4) according to historical data sequence, to the moon, the annual developmental sequence of tolerance carries out timing curve prediction, obtains longitudinal forecasting sequence of to be predicted year
5) lateral prediction sequence and longitudinal forecasting sequence are weighted on average, obtain the two-dimensional prediction sequence of to be predicted year
and add in historical data sequence:
In formula, w
1, w
2the weight of horizontal, the longitudinal predicted value of difference.
2. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 1, is characterized in that, described setting-up time section is greater than 1 year.
3. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 1, is characterized in that, described data pre-service is specially:
The power quality index data that each monitoring point is obtained are weighted the comprehensive historical data sequence that on average obtains forecast interval, and computing formula is as follows:
Y
region=wy
point
w
i=p
i/p
sum,i=1,2,...,N
Wherein, w
ifor i the element of weight vectors w, represent the weight of i monitoring point, p
ibe the power of i monitoring point, p
sumfor the general power of each monitoring point in forecast interval, Y
tegionfor the comprehensive historical data sequence of forecast interval, y
pointpower quality index matrix for each monitoring point in forecast interval.
4. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 1, is characterized in that described step 3) be specially:
301) characteristic parameter of the horizontal curve in year to be predicted is predicted, described characteristic parameter comprises average perunit value ρ and minimum perunit value β;
302) take on the time close to time of to be predicted year be standard year, according to monthly perunit value sequence and the step 301 of standard year) characteristic parameter recording in advance obtains the monthly perunit value sequence of to be predicted year;
303) according to the monthly perunit value sequence of to be predicted year and the new breath value of to be predicted year, carry out newly ceasing famousization processing, obtain lateral prediction value.
5. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 4, is characterized in that described step 301) in, the characteristic parameter Forecasting Methodology of employing comprises the moving method of average, regression analysis or exponential smoothing.
6. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 4, is characterized in that described step 302) comprising:
Will
after descending sequence, become sequence
if the monthly perunit value sequence d of to be predicted year
tsequence postscript is
before and after sequence, under corresponding original of sequence, be labeled as h
j, per days, number scale was T, had following relation:
X
i, jwith the relational expression of characteristic parameter ρ, β be:
B) set up mathematical model:
Wherein,
C) to step b) in mathematical model carry out iterative, iterative process is specially:
C1) introduce Lagrange multiplier w
t=[w
1, w
2..., w
t-1] and v
t=[v
1, v
2], and remember W
0=diag{w
j, W
0e=w, e
t=[1,1 ..., 1]; Put initial value W
0=0, iterations q=1, given condition of convergence ε (ε > 0);
C2) calculate v:v=(AA
t)
-1[b-A (x
(0)+ W
0e)];
C3) calculate x
(*)=x
(0)+ W
0e+A
tv, judgement x
(*)in each component
if
put w
i=0; Otherwise, order
obtain thus x
(*), W
0;
C4) judgement || Ax
(*)||
2/ || b||
2whether < ε sets up, and if so, stops iteration, obtains optimum solution x
(*); If not, put q=q+1, return to step c2);
D) to step c) the optimum solution x that obtains
(*)carrying out contrary generation processes:
7. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 6, is characterized in that described step 303) in, newly cease famousization processing and be specially:
If the power quality index data of the front m of to be predicted year month are known, the new breath value sequence of to be predicted year is { y
t, 1, y
t, 2..., y
t, m, the lateral prediction value in to be predicted year residue month is
Wherein, d
t, j, d
t, kbe respectively the perunit value of the to be predicted year j month, the k month, y
t, kmeet { y
t, k| min (v
(k)), v
(k)computing formula as follows
Obtain thus the lateral prediction sequence of to be predicted year
8. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 7, is characterized in that described step 4) in, when known annual developmental sequence is greater than three values, longitudinally predicted value is identical with the computing method of lateral prediction value.
9. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 1, it is characterized in that described step 4) in, when known annual developmental sequence only has two values, adopt growth ratio method to calculate longitudinal predicted value, its formula is as follows:
Wherein, k=1,2 ..., m, j=m+1, m+2 ..., T.
10. a kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction according to claim 1, is characterized in that described step 5) in, the weight w of horizontal, longitudinal predicted value
1, w
2meet
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310539260.1A CN103544362B (en) | 2013-11-04 | 2013-11-04 | A kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310539260.1A CN103544362B (en) | 2013-11-04 | 2013-11-04 | A kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103544362A true CN103544362A (en) | 2014-01-29 |
CN103544362B CN103544362B (en) | 2016-08-17 |
Family
ID=49967810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310539260.1A Active CN103544362B (en) | 2013-11-04 | 2013-11-04 | A kind of harmonic wave medium-and long-term forecasting method based on two-dimensional curve prediction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103544362B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104010029A (en) * | 2014-05-12 | 2014-08-27 | 上海交通大学 | Distributed computing environment performance predicting method based on transverse and longitudinal information integration |
CN105184425A (en) * | 2015-10-21 | 2015-12-23 | 国网河南省电力公司电力科学研究院 | Reference sample optimizing method for feeder short term load forecast |
CN107977986A (en) * | 2016-10-21 | 2018-05-01 | 北京君正集成电路股份有限公司 | The Forecasting Methodology and device of a kind of movement locus |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622657A (en) * | 2012-04-26 | 2012-08-01 | 上海市电力公司 | Annual maximum power load prediction method |
WO2012114974A1 (en) * | 2011-02-23 | 2012-08-30 | 株式会社日立製作所 | Apparatus for forecasting electric power load of industrial park, system and method |
CN102867221A (en) * | 2012-07-17 | 2013-01-09 | 中国电力科学研究院 | Method for dynamically and transversely correcting medium and long-term electric quantity forecasting |
CN103036231A (en) * | 2012-12-07 | 2013-04-10 | 温州电力局 | Forecasting method, device, and upper computer of power load |
-
2013
- 2013-11-04 CN CN201310539260.1A patent/CN103544362B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012114974A1 (en) * | 2011-02-23 | 2012-08-30 | 株式会社日立製作所 | Apparatus for forecasting electric power load of industrial park, system and method |
CN102622657A (en) * | 2012-04-26 | 2012-08-01 | 上海市电力公司 | Annual maximum power load prediction method |
CN102867221A (en) * | 2012-07-17 | 2013-01-09 | 中国电力科学研究院 | Method for dynamically and transversely correcting medium and long-term electric quantity forecasting |
CN103036231A (en) * | 2012-12-07 | 2013-04-10 | 温州电力局 | Forecasting method, device, and upper computer of power load |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104010029A (en) * | 2014-05-12 | 2014-08-27 | 上海交通大学 | Distributed computing environment performance predicting method based on transverse and longitudinal information integration |
CN104010029B (en) * | 2014-05-12 | 2017-09-08 | 上海交通大学 | DCE performance prediction method based on laterally longitudinal information integration |
CN105184425A (en) * | 2015-10-21 | 2015-12-23 | 国网河南省电力公司电力科学研究院 | Reference sample optimizing method for feeder short term load forecast |
CN107977986A (en) * | 2016-10-21 | 2018-05-01 | 北京君正集成电路股份有限公司 | The Forecasting Methodology and device of a kind of movement locus |
CN107977986B (en) * | 2016-10-21 | 2021-11-23 | 北京君正集成电路股份有限公司 | Method and device for predicting motion trail |
Also Published As
Publication number | Publication date |
---|---|
CN103544362B (en) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Han et al. | Mid-to-long term wind and photovoltaic power generation prediction based on copula function and long short term memory network | |
Costantini et al. | A hierarchical procedure for the combination of forecasts | |
CN104636823B (en) | A kind of wind power forecasting method | |
CN105678403A (en) | Region saturation load prediction method based on model family decomposition and integration technology | |
CN102682218A (en) | Method for evaluating electricity energy efficiency of industrial user | |
CN106682763B (en) | Power load optimization prediction method for large amount of sample data | |
CN104616079A (en) | Temperature change based power grid daily electricity consumption prediction method | |
CN104933629A (en) | Power user equipment evaluation method based on interval level analysis and interval entropy combination | |
CN106446571A (en) | Error evaluation method of wind power forecasting | |
CN104573857A (en) | Power grid load rate prediction method based on intelligent algorithm optimization and combination | |
CN106779277A (en) | The classification appraisal procedure and device of a kind of distribution network loss | |
Sevlian et al. | Value of aggregation in smart grids | |
CN105069536A (en) | Electricity demand predication method based on temperature and economic growth | |
CN102509173A (en) | Markov chain based method for accurately forecasting power system loads | |
CN108876114A (en) | Consider the power grid efficiency benefit appraisal procedure of new energy access | |
CN107293118A (en) | A kind of traffic speed motion interval Forecasting Approach for Short-term | |
CN106934500A (en) | A kind of method being predicted to regional saturation electricity based on nonparametric model | |
CN103530700B (en) | Urban distribution network saturation loading Comprehensive Prediction Method | |
CN101899820A (en) | Method for determining amount of available surface water of river basin facing to protection of river ecosystem | |
CN103544362A (en) | Harmonic medium and long term prediction method based on two-dimensional curve prediction | |
CN106096779A (en) | A kind of power transmission and distribution cost analysis, risk evaluation model and risk structure optimization method | |
CN104134103A (en) | Method for predicting energy consumption of hot oil pipeline by using modified BP neural network model | |
CN107679666A (en) | A kind of power consumption prediction method based on Shapley values and economic development | |
CN104112204A (en) | Evaluation method for efficient operation of power supply quality | |
CN105373865A (en) | Industrial structure based electricity consumption demand prediction method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |