Specific embodiment
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
In the normal operation work of electric system, there are factors that will will affect the variation of electric load, generally comprise
Controllable factor and uncontrollable factor, policy, economic season etc., which all can generate no small influence to power grid, to be influenced, and seems different
There are many connections between factor again.
Wherein major influence factors include:
(1) economic development and the influence of local policy
During national economy rapid growth, people's daily life and all trades and professions electricity consumption demand are also constantly
Soaring, corresponding change can also occur for load.The economy for determining that an area is even national of one economic policy of country
Developing direction can also generate highly important influence to electric system.
(2) influence of power grid electricity price.
Power grid electricity price also results in the variation of load, and power grid can realize the adjusting of load by adjusting electricity price, for example divide
When electricity price setting, tou power price can by reduce peak of power consumption when load, guidance user time-sharing section using electricity, play
The effect of peak load shifting.
(3) influence of meteorological condition and weather.
There are many loads more sensitive for weather and temperature change, such as refrigerator, furnace, sky for electric power demand side
Tune machine etc..These loads are the variations according to weather and change, for example, when time in summer weather is warmmer, especially in China south
A just band, the use time interval of air-conditioning greatly increases in one day, and at this moment power load aggravates;When weather is colder, heating is set
Standby investment increases, and also results in power load exacerbation.
(4) influence of time factor.
Time is also one of factor for influencing load variations, for load prediction, main time factor
Can be divided into following three kinds: two-day weekend with it is workaday intersect, replacement, the appearance of great festivals or holidays of seasonal climate transformation.For
For short term, the variation of load weekly causes short term to change, and working day is different from the load of two-day weekend
, on Monday belong to regular working day to Friday, commercial power biggish for power demand, workload demand amount is larger, in week
Six arrive Zhou Tian, and factory, which has a holiday, causes load to reduce.Great festivals or holidays are often due to the form of expression of business model increases power supply
Perhaps, demand, such as the Spring Festival, International Labour Day, the Dragon Boat Festival, National Day etc. can make electricity change.For term Load,
The replacement in season throughout the year leads to the variation of term Load, and the temperature in summer is higher, and the temperature in winter is lower, and weather is sensitive
Load using more, is become more for the demand of electricity.
(5) enchancement factor etc.
Enchancement factor, such as natural calamity, affect electric system, also result in load and generate unexpected variation.
Perhaps, some natural calamities also will affect the acquisition of load data, and load data is caused abnormal point occur, this number to load prediction
It is had an impact according to acquisition, so as to reduce the accuracy of prediction result.
Load Prediction In Power Systems are combined with society's electricity consumption situation, economy, weather, politics and other factors, to load
Supply and demand development trend carries out reasonable prediction.According to the correlation of the development trend of analysis of history load data influence factor associated therewith
Property, the prediction of science is made to the following electricity consumption.It, can be by load prediction point according to the difference of prediction purpose and predicted time range
For three classes:
Short-term load forecasting (STLF), usually 1 hour to 1 week, short-term load forecasting for generating control in real time, safety
Analysis and power exchange plan etc. operate extremely important.Compared to Mid-long term load forecasting, carrying out short-term load forecasting is in order to short
When control scheduling various regions power supply volume, on the one hand can help to formulate the fuel of generator in a short time for power supply side and supply
Plan, improves economy and reasonability;On the other hand facilitate power grid and formulate interruption maintenance plan.It can be also used for pacifying in real time
Complete analysis, Real-time Economic Dispatch etc..The characteristics of short-term load forecasting is to predict day and a few days ago with the rule of the transient change of period
Restrain it is more similar, different maturity periods it is interior have similar rule.The influence factor of short-term load forecasting mainly has week type, electricity
Valence, meteorologic factor etc..
Medium term load forecasting (MTLF), a usually Zhou Zhiyi, medium term load forecasting with from a few days to a few weeks or several
The period of the moon is related, and for meeting the load requirement in summer or winter peak period.The purpose of medium term load forecasting be for
The interruption maintenance plan of route and substation is formulated, control centre adjusts the method for operation etc. of electric system.Term Load is pre-
The load variations rule that is mainly characterized by surveyed has periodically, for example, the load variations trend in observation several years, can obtain and exist every year
The load trend in identical month has similitude.
It is more than 1 year that long term load forecasting (LTLF), which is often referred to the time,.Long term load forecasting is mainly used for Electric Power Network Planning, is
Power supply expands the macro-level policy-makings such as capacity, network structure upgrading, system scale expansion and provides strong foundation.Long term load forecasting
Feature is that data are substantially to be monotonically changed, and generally shows incremental, aperiodicity.The influence factor of long-term phase load prediction
Mainly there are national economic development situation, the variation of population, the adjustment of the industrial structure, variation of electrovalence policy etc..
Above-mentioned short-term, medium and long term load prediction, each type have important role, influence the reliable of electric system
Property and economy.In smart grid, electricity market is opened to electricity provider and consumer, for instance it can be possible that a company can
Using the company as normal product processing electric power, and weigh and consider the balance of demand response.Accurate medium term load forecasting and
Short-term load forecasting, which formulates aid decision making person in electricity market, suitably to be planned and to obtain maximum return most important.
It should be noted that medium term load forecasting is concentrated mainly on 1-5 years, the prediction generally as unit of year is main
Act on mid-term balance of electric power and ener, power supply and substation's constant volume addressing, Net Frame of Electric Network planning etc..Long term load forecasting refers generally to
Even longer-term, main function is to plan that electric system distant view develops within 10-20 years.The master of Mid-long term load forecasting
Feature is wanted to have:
(1) forecast sample is few: the time major part of Electric Power Network Planning load prediction was usually no more than 15 years in 10-15 years,
Shorter is possibly less than 10 years, this leads to have limitation in selection main reasons is that complete load sample negligible amounts.Its
His reason mainly includes that historical statistics means are insufficient, statistical error problem is larger and economic society transition problem, available
Sample size is also limited.Based on the few feature of load prediction sample, should be avoided as far as possible when carrying out the selection of load forecasting model
To the demanding model of sample size, and select efficiently excavate the model for containing effective information in small sample.
(2) predicted time span is larger: in Mid-long term load forecasting, the range of time span at least more than 1 year, this
Challenge is proposed to prediction precision.Impact factor is broadly divided into short-term factor and long term factor, and short-term factor is main
Comprising temperature change, rainfall situation, humidity condition etc., long term factor has financial crisis, industry restructuring, regional industrial to turn
Shifting, high energy-consuming industry development, income level of resident, electrovalence policy, subsidy policy stimulation etc., short-term factor to Mid-long Term Load not
Excessive influence can be generated, and long term factor is affected to Mid-long Term Load variation.It follows that economic development is in middle length
Phase load forecast has apparent influence, therefore should focus on collecting the related data information of economic development when being predicted,
The prediction model for establishing multivariable improves the precision of prediction of load prediction.
(3) demand of the prediction to section is more urgent: Mid-long term load forecasting mainly macro-plan and in long-term electric power
Obvious effect is planned, policymaker is desirable with effective information, obtains effectively as a result, providing the upper limit of predicted load under
Limit provides reasonable suggest to planning.
(4) Electric Power Network Planning, traffic control experience accurately play an important role to long term load forecasting.Mid-long term load forecasting has
The features such as time span is big, sample size is few, influence factor is more, therefore work as and possess management and running work expertise pair abundant
Its precision predicted has a preferable supplement, therefore more focuses on qualitative subjective analysis in system modelling and have with objective quantitative
Machine combines, and enhances the accuracy, reasonability and stability of forecast result of model.
Currently, experts and scholars both domestic and external propose the method for many load predictions, most of Predicting Techniques are negative
It carries prediction aspect to attempt, achieves different degrees of success.Since the time of load prediction divides the difference of classification, may lead
Cause identical method that may generate different results to different types of load prediction.These especially popular technologies in recent years
It include: to have econometrics model and statistics for long term load forecasting.There is neural network for short-term load forecasting,
Fuzzy logic, expert system, homing method, time series, similar day method and support vector machines.Due to influencing medium term load forecasting
Factor between short-term load forecasting and long term load forecasting, therefore be wherein suitable for it is long-term and short-term load forecasting perhaps
Multi-method can be applied to medium term load forecasting.
Load is divided into several types according to the length of time, for Short Term load Forecasting Technique, short term
The grasped data information of prediction is compared with horn of plenty, and the variation of load in a short time will not be particularly big, and influence factor does not have
Mid-long term load forecasting is so much, therefore the research of current short-term load forecasting method is more more mature, and artificial intelligence
With the technological innovation in terms of information processing, short-term load forecasting method has been got back further supplement and improvement.Compared to short
For phase load prediction, the span time of Mid-long term load forecasting is big, and influence factor is more and uncontrollable, and data information is collected relatively tired
It is more difficult will to carry out in this respect accurate quantitative study for difficulty.
In view of the above technical problems, referring to Fig.1, a kind of prediction technique of electric system Mid-long Term Load of the present invention is shown
Step flow chart, which comprises
Step S101: the Mid-long Term Load data in electric system are collected, the Mid-long Term Load data are located in advance
Reason, obtains load data sample;
Characteristics of Electric Load index is the quantity performance of part throttle characteristics, that is, the characteristic value of load variations.In last generation
" index explanation of power industry production statistics " file that Ministry of Energy, discipline the eighties China issues is proposed about part throttle characteristics
14 indexs, Guo Wang company has carried out addition supplement to it after 21st century, wherein most importantly increasing one
New index --- peak-valley ratio.The calculating of Load characteristics index and analysis are non-for the load variations characteristic for describing electric system
It is often necessary.By to part throttle characteristics carry out analysis obtain peak load using the development trend of hour be part throttle characteristics prediction base
This content.Load curve characteristics index is to carry out load Analysis, grasp rule and characteristic that load changes over time, is primarily determined
The method of operation in all kinds of power plant power stations and the parameter of installed capacity.
The index of Load characteristics index has very much, both at home and abroad also not no index system of unified specification, considers practicability,
Main loads characteristic index is divided into a year Load characteristics index, moon Load characteristics index and daily load characteristic index according to the time period.
The selection of Load characteristics index can largely influence the quality of load Analysis, and the method for load prediction also can be directly to pre-
Result is surveyed to have an impact.Load characteristics index can be divided into description class substantially according to content, compare class and class of a curve, referring to figure
2, show the classification situation of Load characteristics index in the prior art.As can be seen from Figure 2, description class mainly includes maximum load, most
Underload, average load, peak-valley difference.Comparing class mainly includes daily load rate, average daily load rate, ratio of minimum load to maximum load, peak-valley difference
Rate, maximum load utilize hourage etc..Class of a curve includes daily load curve, monthly load curve, yearly load curve, year continuous loading
Curve.
Mid-long Term Load data collected by the present invention refer to above-mentioned Load characteristics index.As the basis of load prediction,
Historical data will directly influence the reliability of prediction result, and load prediction is that basis carries out analysis expansion to historical data.?
It is general using SCADA system acquisition initial data (the Mid-long Term Load data i.e. in the present invention) in electric system.
In actual motion, the links of data collection system may be because of some external factors such as: communication meets with
To interference, losing occur in data in data transmission procedure, may result in the missing and exception of data acquisition.In addition, having
When some special event cutting loads have a power failure, line maintenance has a power failure, and large user, major issue impact etc. can cause load data to be dashed forward
Become.As lacked, being mutated these abnormal data, it is pre- that these exceptional values are directly applied into load without being pocessed if directly
Among survey, it will so that the accounting of load prediction technology and these abnormal datas, then, when being fitted emulation, perhaps
Normal match value can be deviateed because of these exceptional values, obtain less appropriate predicted value, cause forecasting accuracy lower.Cause
This, the present invention pre-processes initial data, to guarantee the foundation of correctly predicted model and the raising of precision of prediction.
In a preferred embodiment, the Mid-long Term Load data collected in electric system are proposed, in described
The step of long-term load data are pre-processed, and load data sample is obtained include:
Sub-step 1: the Mid-long Term Load data in electric system are collected;
Sub-step 2: the first pretreatment is carried out to the Mid-long Term Load data, obtains the first Mid-long Term Load data;
Sub-step 3: data smoothing processing is carried out to the first Mid-long Term Load data, obtains load data sample.
The embodiment of the present invention carries out the first pretreatment to collected Mid-long Term Load data first, then again to pretreatment
Data afterwards carry out data smoothing processing and remove some improper datas or special data by handling twice, reduce data
Fluctuation, to improve precision of prediction to a certain extent.
Data across comparison method, data longitudinal direction method of comparison, curve displacement method, probability statistics can be used in above-mentioned first pretreatment
One of method and the supplement of missing data or any several method, in which:
(1) data across comparison method
The across comparison method of data is using the load value of adjacent moment as reference, by the load value and ginseng of specified time point
Examine value to be compared, if the difference of the two has been more than a certain threshold value, then it is assumed that the load value at the moment be exceptional value, need to its into
Row repairing treatment.
(2) data longitudinal direction method of comparison
Electric load has certain periodicity, therefore the data variation of front and back mutually in the same time is little, they have certain
Similitude, if the difference of a certain time data and front and back load value in the same time has been more than some threshold value, then it is assumed that the data are
Bad data, using the load average value of front and back in the same time as the correction value at the moment.
(3) curve displacement method
Obvious abnormal and major break down day the load of load curve is rejected or replaced.
(4) probabilistic method
The primary election and amendment for rule of thumb carrying out Power system load data first, then further according to needing to be arranged appropriate confidence
Section reaches the identification and amendment to bad data in data acquired by this two step.
(5) supplement of missing data
Containing scarce in the data that the problems such as having a power failure as caused by information acquisition system reliability or special circumstances occurs
The case where losing data, the processing to this kind of data is that several days load datas in front and back are taken mean value to fill up.
Current data smoothing processing method includes " the simple rolling average of 2n+1 point " smothing filtering method, " weighting is mobile flat
" exponential smoothing, seasonal index smoothness etc..
In practical operation power, the production activity of people, the continuous variation of life style and meteorological condition, electric power
Supply can change with these variation, and this variation has diversity and complexity.Monthly load is come
It says, changing rule can be unfolded in terms of annual developmental sequence and monthly developmental sequence two.On the one hand, by the sequence of year development
On from the point of view of, with the development of social economy, the living standard of the people significantly improved, and continuous increasing occurs in the monthly load of electric system
Long situation;On the other hand, in each moon load in 1 year, with certain fluctuation tendency for being showed of variation in season.
Attribute annual period of load mainly has the following aspects:
The growth property of trend, two continuous time loads have same distribution character;
The fluctuation in season, general area one, February load level it is relatively low under, but to seven Augusts of summer
Its load level of part, which just has, to be obviously improved;
There is low ebb in the attribute of festivals or holidays, the load when general trend of usual regional load is grand festivals or holidays,
His period is in higher horizontal fluctuating change up and down.
The moon cyclic attributes of load are mainly:
A month the inside, load was substantially being changed using week as the period.It is workaday in a week
Load level is higher than the load level of two-day weekend.This is mainly due to the compositions of load and people's rule of life habit to be determined
's.The main constituents for the daily load that works are industrial load, but the electricity consumption proportion in two-day weekend resident living can be bright
Aobvious to rise, because of the specific gravity that network load shared by industrial load is very big, two-day weekend load should be able to be significantly lower than working day
Load.In addition to this, because industrial load on weekdays during generally lie in stable operating, workaday load
Variation has similitude.It is so analyzed from attribute annual period of load, similar load distribution performance can be to monthly load
Prediction provides accurate power load distributing trend.It is analyzed by the moon cyclic attributes of load, the prediction for the daily load that works will conduct
Influence the mostly important reason of the accuracy of the prediction modeling of monthly load.
Patentee considers the fluctuation problem of monthly load, in the present invention preferably using seasonal index smoothness to institute
It states the first Mid-long Term Load data and carries out data smoothing processing, to obtain accurate load data sample.
In a preferred embodiment, it shows according to seasonal index smoothness, to first Mid-long Term Load
Data carry out data smoothing processing, and the sub-step for obtaining load data sample further comprises:
Choose the first Mid-long Term Load data in multiple years;
The first Mid-long Term Load data are calculated in the average value of the same target time period in the multiple year, Yi Ji
The overall average in all months in the multiple year;
The average value of the same target time period is divided by with the overall average, the season for obtaining the target time period refers to
Number;
The seasonal index number of the average value and the target time period of the same target time period is divided by, when obtaining the target
The load data sample of phase.
After obtaining the load data sample, need to predict the load data sample.Current prediction technique
It is broadly divided into traditional prediction method and modern prediction technique.
Traditional prediction method is to build high-precision data model by collecting research data, and basic theories is derived from
The mathematical theories such as Probability Theory and Math Statistics.Such as: time series method, regression analysis and gray forecast approach prediction technique etc..
Modern prediction technique includes intelligent algorithm, fuzzy prediction method, wavelet analysis method and combinatorial forecast etc..Its
In, intelligent algorithm includes expert system, artificial neural network method etc..The method of load prediction is more, but each side
Method has the advantage and disadvantage and the scope of application of oneself, and inventor compares several common prediction techniques, and comparison result is such as
Shown in table 1.
Table 1
As known from Table 1, gray forecast approach prediction technique has that algorithm is simple, method is mature, calculation amount is small, speed is fast, sample
The advantages that requirement is small, and characteristic quantity need not be considered in modeling, the index prediction of nonlinear change is applied also for, and it is medium-term and long-term
The Property comparison of load data is close, so the present invention uses in gray forecast approach prediction technique as prediction side of the invention
Method, the prediction model with gray model GM (1,1) as load data sample.
Step S102: the load data sample is input in gray model GM (1,1), obtains the first predicted value;
In a preferred embodiment, it shows and the load data sample is input to gray model GM (1,1)
In, the step of obtaining the first predicted value includes:
Sub-step 4: the load data sample is input in GM (1,1) model;
Sub-step 5: in the GM (1,1) model, single order Accumulating generation is carried out to the load data sample, obtains one
Rank load sequence;
Sub-step 6: being fitted the changing rule of the single order load sequence using first-order equation, obtains the GM (1,1) mould
The time response function model of type;
Sub-step 7: the time response function model is subjected to the reduction of single order regressive, obtains the load data sample
Gray prediction result;
Sub-step 8: the gray prediction result is subjected to average generation processing, obtains the first predicted value.
It in embodiments of the present invention, can be by separation and war, irregular by carrying out single order Accumulating generation to load data sample
Over-and-over addition is generated new cumulative data sequence by initial data, after accumulated generation, so that it becomes having exponential increase rule
The ascending series of rule.
The specific implementation step of sub-step 6 includes:
The changing rule of the single order load sequence is fitted using first-order equation;
Use the undetermined parameter of single order load sequence described in least square solution;
The undetermined parameter is substituted into the first-order equation, obtains the time response function model of the GM (1,1) model.
Above-mentioned first-order equation, least square method and undetermined parameter belong to the common algorithms and data of Grey Model method, herein
Seldom repeat.
Next, by time response function model progress single order regressive reduction, at the result after Accumulating generation prediction
Reason, so that ascending series are reduced to normal ordered series of numbers.
Sub-step 8 is using average generation processing method to carry out school to the gray prediction result by the reduction of single order regressive
Just, the first predicted value of final gray model is obtained.
Above-mentioned average generation includes adjacent average generation, non-neighboring average generation.When adjacent average generation has equal ordered series of numbers
Away from the average value of adjacent two data, can subtract by the data handled in this way under using even time interval when generating new data
It is irregular under weak anomaly and emergency case, while can effective polishing defect value;With adjacent average generation on the contrary, non-neighboring average generation
Ordered series of numbers when away from be it is unequal, generally use and delete the exceptional value that occurs in ordered series of numbers, contain to be formed in entire ordered series of numbers
Vacancy so that ordered series of numbers when away from equal, after abnormal data is deleted, deleted data are often asked using the data on both sides
Average value replaces, thus the position that fills up the vacancy.
Step S103: compared with first predicted value is made difference with the Mid-long Term Load data, residual sequence is obtained;Its
In, the residual sequence takes absolute value;
Since load prediction work is the judgement to load future development rule, to predict following certain periods
Load value and load variations rule.Some uncontrollable factors are often accompanied by the prediction of load, these factors may cause to pre-
The reliability and accuracy for surveying result have an impact.If there is no selection suitable appropriate prediction model when being predicted, it is this
The accident changed at random all can make predicted value that can have a certain difference with actual value.Therefore, in step s 103, of the invention
Compared with first predicted value is made difference with the Mid-long Term Load data, school is carried out to the prediction result of above-mentioned prediction technique
It tests, to evaluate the precision of above-mentioned grey forecasting model.Since residual sequence has positive value also to have negative value, and gray model requires data
Sequence symbol is consistent, cannot directly use Grey Model, so making it meet ash herein by the absolute value of residual sequence
The requirement of color model.
Specific steps are as follows:
It is assumed that the sequence of history Mid-long Term Load data and the first predicted value is respectively X(0)(k)、Wherein, k=
At the time of 1,2 ..., n indicate record.ε (k) is represented by the residual error at k moment, i.e. actual value X(0)(k) and predicted valueIt
Between difference.The average value of residual epsilon (k)It may be expressed as:
The wherein average value of history Mid-long Term Load dataAre as follows:
History Mid-long Term Load data varianceExpression formula are as follows:
Residual varianceExpression formula are as follows:
Posteriority difference ratio C are as follows:
Small error possibility P are as follows:
Equation left side two indexes C, P value range in formula (5) and (6) are as follows: C > 0,0≤P≤1.C index value is smaller, indicates
The precision of predicted value is higher.When the residual error fitting degree that history Mid-long Term Load data are more discrete, and predicted value and actual value generate
Higher, C is smaller.When P is bigger, the difference of residual error and residual error average value is less than 0.6745S1Point it is more, thus, it can be known that with index C
On the contrary, the bigger precision of prediction of P is higher.According to above-mentioned two index parameter, establishes such as table 2 and evaluate obtained by grey forecasting model
Predicted value precision.
Table 2
Precision of prediction grade |
P |
C |
Precision of prediction grade |
P |
C |
Good (level-one) |
>0.95 |
<0.35 |
(three-level) reluctantly |
>0.7 |
<0.45 |
Qualified (second level) |
>0.8 |
<0.5 |
Unqualified (level Four) |
≥0.7 |
≥0.65 |
Based on precision evaluation standard described in table 2, in a preferred embodiment, showing evaluation method includes:
According to the residual sequence, judge whether first predicted value meets preset requirement;
If meeting, using first predicted value as the prediction result of the Mid-long Term Load data.
Above-mentioned preset requirement refers to using after residual sequence progress P, C metrics evaluation, whether judges the first predicted value
It is qualified.
Grey GM (1,1) model can simulate the basic trend of load sample data, and obtain predicted value, but such as
Fruit only uses grey forecasting model only come if being modeled, the prediction result that cannot be more satisfied with, trace it to its cause for,
Gray model cannot simulate the random information occurred in load sequence well.So the first predicted value be likely to occur it is qualified with not
Qualified two kinds of situations.
When the first predicted value qualification, directly using first predicted value as the prediction knot of the Mid-long Term Load data
Fruit.
Step S104: when first predicted value does not meet preset requirement, the residual error is corrected using Markov Chain
Sequence obtains the second predicted value of the residual sequence.
But when the first predicted value is unqualified, illustrate that there are certain precision problem, inventors by this gray model GM (1,1)
The residual sequence is modified using Markov Chain, the prediction that gray model GM (1,1) is obtained can be further decreased
The fluctuation of value.
Markov Chain mentioned by the present invention is the method for describing Stochastic Dynamic Process, it refers to system in each time
State in which be it is random, and unrelated with state before, that is, have the characteristics that markov property.Markov chain model is logical
It crosses and determines the current status of variable to predict the development tendency in following a period of time.If the transfer of system generating state is
Refer to the case where system variable is transferred to the value range of another state from the value range of a state.
Before formally being corrected using the Markov Chain amendment residual sequence, need gray model and Markov Chain
Be combined together, the method specifically includes the following steps:
The residual sequence is input in the GM (1,1) model, the gray prediction sequence of the residual sequence is obtained;
The residual sequence and the gray prediction sequence are divided by, the grey degree of fitting index of the residual sequence is obtained;
According to preset state classification condition, the grey degree of fitting index is divided into several state intervals;
According to several described state intervals, the Markov state transition probability square of the gray prediction sequence is constructed
Battle array.
In a preferred embodiment, before correcting the residual sequence using Markov Chain, also to the horse
Er Kefu state transition probability matrix carries out geneva inspection.
The implementation method of step S104 includes:
According to the Markov state transition probability matrix, obtains the gray prediction sequence and correspond to the state interval
Interval range;
Preset value is taken in the interval range, and the preset value is multiplied with the residual sequence, is obtained described residual
Second predicted value of difference sequence.
Above-mentioned steps are illustrated using formula below.
By obtaining a residual sequence after GM (1,1) model prediction, q can be denoted as(0)(k), it may be assumed that
Wherein x(0)(k) true Mid-long Term Load data are represented,Represent the predicted value of GM (1,1) model.
After the absolute value of residual sequence, then residual sequence ε(0)It can indicate are as follows:
ε(0)=(ε(0)(1),ε(0)(2),ε(0)(3),…,ε(0)(n)) (8)
By residual sequence ε(0)Be input in grey GM (1,1) model and predicted, establish about residual sequence GM (1,
1) model obtains the gray prediction sequence of residual sequence:
Parameter thereinWithIt can be obtained by least square method.
Calculate the grey degree of fitting index of residual sequence:
The gray model matched curve of residual sequence is rendered as the characteristic of index, and grey degree of fitting index Y (k) is mainly used
Reflect that original data surround the degree of fluctuation of matched curve, that is, reflects the degree of fluctuation of residual sequence.Due to residual
The variation of difference sequence shows a kind of the characteristics of non-stationary randomness, it is possible to be corrected by markovian application residual
The result of the prediction of GM (1,1) model of difference sequence.The basic model of Markov Model about Forecasting method are as follows:
X (k+1)=x (k) × P (11)
Wherein, x (k) indicates that quantity of state of the prediction object at the t=k moment, x (k+1) indicate prediction object in t=k+1
The quantity of state at quarter, P indicate the step transition probability matrix in Markov Chain.
1, the grey degree of fitting index of residual sequence is divided
Since grey fitting accuracy index is the non-stationary process of a random fluctuation, thus need to consider a state classification
Method, so that Y (k) is divided into multiple states.According to the size of Y (k), state interval can be divided:
Ei=[Yi1,Yi2], i=1,2 ..., k (12)
Wherein, Yi1That indicate is the lower limit of Y (k), Yi2What is indicated is the upper limit of Y (k).Under normal circumstances, the shape of division
State is more, and the precision of model is higher.
2, structural regime transition probability matrix
The prediction result of residual sequence is divided into several states:
Pij=nij(k)/ni (13)
Above formula indicates system by state EiBy k step to state EjThe number of appearance is nij(k), EiThe number of appearance is ni。
Then state transition probability matrix are as follows:
3, geneva is examined
It is (n-1) that statistic, which obeys freedom degree,2Distribution.Selected confidence level is α, can be by tabling look-up to obtain χa 2((n-1)2),
It enables:
IfThen it is with geneva.
4, the amendment of predicted value
It, can be with the corresponding section of grey degree of fitting index of prediction residual sequence by Markov state transition probability matrix
Range [Yi1,Yi2], preset value is then taken in the interval range, and the preset value is multiplied with the residual sequence, obtain
Obtain the second predicted value of the residual sequence.It is preferred that take state interval intermediate value and raw residual sequence product as residual error sequence
The correction value of column predicted value, i.e. the second predicted value.
Step S105: judge whether second predicted value meets preset requirement;
Step S106: if meeting, using second predicted value as the prediction result of the Mid-long Term Load data.
The verification method of the first predicted value of above-mentioned preset requirement reference, not described here any more.It, will when meeting preset requirement
Prediction result of second predicted value as the Mid-long Term Load data.
In addition, when second predicted value does not meet preset requirement, which comprises
When judging that second predicted value does not meet preset requirement, second predicted value is added to described medium-term and long-term
In load data, and second predicted value and the Mid-long Term Load data are pre-processed.
After second predicted value is added in the Mid-long Term Load data, subsequent processing steps are referring to centering of the present invention
The processing step of long-term load data, not described here any more.
To sum up, grey forecasting model has the characteristics that can handle " small sample ", " poor information ", and the present invention is based on its characteristics
Applied among long-medium term power load forecasting, the forecast sample to reduce Mid-long term load forecasting is few, predicted time across
Spend biggish problem.But for the load with fluctuation, the prediction of gray model it is ineffective, trace it to its cause for
Since gray model is the characteristics of input data is made accumulation process, makes it have exponential increase, once the sample data of input
When fluctuating larger, the precision of prediction will be reduced.To solve the problems, such as data fluctuations caused by Accumulating generation, centering of the present invention is long
Phase load data is pre-processed, and can be reduced to the fluctuation of Mid-long Term Load data, wherein especially utilizes seasonal index number
Exponential smoothing reduces fluctuation caused by Accumulating generation, and then by treated, sequence is predicted again, is finally gone back prediction result
Original, the prediction result which is replaced.In addition, occurring since grey forecasting model cannot simulate well in load sequence
Random information, and demand of the Mid-long term load forecasting to section is more urgent at present, the present invention using Markov Chain come
State interval is divided, the fluctuation tendency of load is established by probability transfer matrix, by grey forecasting model and Markov Chain phase
In conjunction with, the fluctuation of Mid-long Term Load data can be reduced again, finally realize a kind of prediction effect it is accurate, it is reasonable and stable in
Long term load forecasting method, to overcome existing technical problem.
Next, technical solution of the present invention is described in detail using a specific example.
Zhangye Prefecture's in January, 2011 in December, 2016 is chosen, 72 totally months Analyzing Total Electricity Consumption data are as original
Data are emulated.Specific initial data is shown in Table shown in 3.1, and the load of 2011 in monthly, 2016 electricity datas is at any time
Distribution map is as shown in Fig. 4 .1.
Table 3.1
The monthly electricity data of analysis 2011 to 2016 it can be concluded that, the song of the data of Zhangye Prefecture and Non-smooth surface
Line has trend growth property and cyclic swing.For monthly electricity, with the replacement in season, people's lives habit
Change, can all possess identical trend every year, for water power, since seasonal law has the area of dry season and wet season
Not;For people's household electricity, hot and winter the cold in summer leads to the use of the equipment of the adjusting temperature such as air-conditioning
Increase, these phenomenons are but also the characteristic of the load of every month is different.It can use gray model and solve the problems, such as trend growth,
The fluctuation pattern of load is solved using Markov Chain.
It is modeled using the monthly electricity data in Zhangye Prefecture in January, 2011 in December, 2016 as initial data,
Using the monthly electricity data in January, 2017 to December as test data, modeled using MATLAB R2016a platform.
It, can be by grey with two kinds of trend of fluctuation and trend growth property from Fig. 4 .1 it is found that for monthly
Model predicts its trend growth property, if but directly original loads data be input among gray model predict, ash
Color model is, then vulnerable to the influence of fluctuation, to cause prediction result undesirable with exponential growth.So carrying out grey
It before the prediction of model, needs to carry out initial data the smoothing processing of data, the fluctuation of initial data is reduced, then will place
Data after reason, which are input in gray model, to be predicted.
(1) 2011 to 2016 corresponding average value of each month is calculated first;
aiIndicate i-th month corresponding average value, bijRepresent the electricity data of the i month in jth year.
(2) the average electricity consumption in all months, i.e. overall average electricity consumption are calculated;
(3) electricity consumption that is averaged of each phase divided by overall average electricity consumption is obtained into the seasonal index number of each phase;
(4) it is modified monthly to be obtained into seasonal index number divided by the seasonal index number in corresponding month for monthly electricity consumption data sequence
Electricity consumption sequence.
2011 to 2016 seasonal index numbers calculated according to formula (1.1)~(1.3) are as shown in table 3.2.Table 3.2
For 2011 in monthly, the 2016 electricity data seasonal index numbers in Zhangye Prefecture.
Table 3.2
Same month sequence of average as shown in table 3.2, in data smoothing treatment process are as follows:
A=[16.458 14.225 14.667 15.378 15.339 16.328 ... 16.815 18.253
17.169 15.973 16.774 16.634]
Overall average electricity consumption is 16.168 hundred million KWh;Seasonal index number sequence are as follows:
ρ=[1.018 0.880 0.907 0.951 0.949 1.010 ... 1.040 1.129 1.062 0.988
1.038 1.029]
Through the value of revised 2011 in monthly, 2016 electricity of seasonal index number as shown in table 3.3, seasonal index number
Smoothed out load sequence is as shown in Fig. 4 .2.Table 3.3 is revised 2011 in monthly, 2016 electricity of seasonal index number
Value.
Table 3.3
From Fig. 4 .2 it is found that the smoothed out load sequence of seasonal index number compared with the load sequence before amendment for, fluctuation does not have
So obvious, only some tiny fluctuations are whole to have the trend increased upwards, can be by the smoothed out load sequence of seasonal index number
Column, which are input in gray model, to be predicted.
Establish the gray model of seasonal index number Orders Corrected:
It can be obtained by MATLAB platform emulation, parameter a and parameter u are as follows: a=-0.0027, u=14.671
A and u can be substituting in the time response function model of the GM (1,1) model, seasonal index number Orders Corrected
Gray model are as follows:
Using the gray model of foundation, in January, 2011 in December, 2016 is simulated, the value of fitting is 3.4 institute of table
Show, table 3.4 is 2011 to 2016 fitting results.The match value and the smoothed out load sequence pair of seasonal index number of gray model
Than scheming as shown in Fig. 4 .3.
Table 3.4
Utilize the monthly electricity in January, 2017 to December obtained in the smoothed out sequence inputting to gray model of seasonal index number
Predicted value is measured as shown in table 3.5.Table 3.5 is the monthly power quantity predicting value in January, 2017 to December.
Table 3.5
Analyze fitting result, table 3.5 it can be concluded that, the average relative error of gray model match value is 3.05%, effect
Preferably, but highest relative error has reached 9.20%, needs to carry out some amendments to it.
Next, calculating the error of gray model, residual sequence is constituted:
In formula:For the predicted value of gray model,For actual value.
Residual sequence figure is as shown in Fig. 4 .4 constructed by the present invention.
Can be seen that residual sequence from Fig. 4 .4 has positive value also to have negative value, and gray model requires data sequence symbol one
It causes, Grey Model cannot be directly used, so making it meet wanting for gray model the absolute value of residual sequence herein
It asks.
Will | ξ(0)(k) | then sequence uses GM (1,1) model, residual absolute value is calculated as original data sequence
Forecasting sequence
It can be obtained by MATLAB platform emulation, parameterAnd parameterAre as follows:
The match value that residual sequence can be calculated by above formula, as shown in table 3.4.Residual sequence match value calculates
Later, ash discharge precision index is calculated, several dynamical states are divided according to grey precision index, according to the point for falling into each state interval,
Markov Transition Probabilities matrix is calculated, the correction value of residual error is obtained.
Analyzing result to can be seen that the relative error after residual GM by table 3.4 is 1.72%.It is repaired compared to residual error
2.81% before just, precision of prediction has obtained better promotion.Gray model match value and seasonal index number are flat after residual GM
Match value and actual comparison figure such as Fig. 4 .6 and Fig. 4 .7 institute of the load sequence after cunning as shown in Fig. 4 .5, before and after residual GM
Show.
It is predicted using data of the gray model after residual GM in January, 2017 to December, obtained predicted value
Final predicted value should can be just reduced to multiplied by seasonal index number, final prediction result is as shown in table 3.6.Table 3.6 is
The monthly power quantity predicting value in January, 2017 to December.
Table 3.6
It can analyze and obtain from the predicted value in table 3.6, minimal error reaches 0.06% accuracy, average relative error
It is 1.86%.The result illustrates the gray model based on Markov Chain residual GM in the moon in January, 2017 to December
Degree power quantity predicting achieves preferable effect.
In order to prove the validity of this method, the present invention is by the predicted value of this method and wavelet neural network, ARIMA model
It is compared, as shown in Fig. 4 .8.
Can be seen that the method for the present invention predicted value from Fig. 4 .8 to be close with actual curve, from residual plot it can also be seen that
Compared with wavelet neural network and ARIMA model, the prediction effect of the method for the present invention is preferable.
It should be noted that for simple description, therefore, it is stated as a series of action groups for embodiment of the method
It closes, but those skilled in the art should understand that, embodiment of that present invention are not limited by the describe sequence of actions, because according to
According to the embodiment of the present invention, some steps may be performed in other sequences or simultaneously.Secondly, those skilled in the art also should
Know, the embodiments described in the specification are all preferred embodiments, and the related movement not necessarily present invention is implemented
Necessary to example.
In view of the above technical problems, referring to Fig. 3, a kind of prediction meanss of electric system Mid-long Term Load of the present invention are shown
Structural schematic diagram, described device comprises the following modules:
Data preprocessing module 301, for collecting the Mid-long Term Load data in electric system, to the Mid-long Term Load
Data are pre-processed, and load data sample is obtained;
Gray model processing module 302 is obtained for the load data sample to be input in gray model GM (1,1)
First predicted value out;
Residual sequence obtains module 303, compared with first predicted value is made difference with the Mid-long Term Load data,
Obtain residual sequence;Wherein, the residual sequence takes absolute value;
Markov Chain correction module 304, for utilizing Ma Erke when first predicted value does not meet preset requirement
Husband's chain corrects the residual sequence, obtains the second predicted value of the residual sequence;
Preset requirement judgment module 305, for judging whether second predicted value meets preset requirement;
Prediction result correction verification module 306, for when judging that second predicted value meets preset requirement, by described second
Prediction result of the predicted value as the Mid-long Term Load data.
For system embodiments, since it is basically similar to the method embodiment, related so being described relatively simple
Place illustrates referring to the part of embodiment of the method.
All the embodiments in this specification are described in a progressive manner, the highlights of each of the examples are with
The difference of other embodiments, the same or similar parts between the embodiments can be referred to each other.
Above to a kind of prediction technique and device of electric system Mid-long Term Load provided by the present invention, carry out in detail
It introduces, specific examples are applied in the present invention, and principle and implementation of the present invention are described, and above embodiments are said
It is bright to be merely used to help understand method and its core concept of the invention;At the same time, for those skilled in the art, foundation
Thought of the invention, there will be changes in the specific implementation manner and application range, in conclusion the content of the present specification is not
It is interpreted as limitation of the present invention.