CN103984988A - Method for correcting super-short-term prediction of photovoltaic power of ARMA module in real time through light metering network - Google Patents

Abstract

The invention discloses a method for correcting super-short-term prediction of the photovoltaic power of an ARMA module in real time through a light metering network. The method includes the steps of inputting data to obtain auto-regression sliding average model parameters, inputting light resource monitoring system data and running monitoring system data, correcting the startup capacity in real time according to the running monitoring system data, setting up an auto-regression sliding average model so as to obtain a photovoltaic power super-short-term prediction result, and introducing real-time light metering station data so as to correct the photovoltaic power super-short-term prediction result in real time. By introducing the real-time light metering station data to correct the photovoltaic power super-short-term prediction result in real time, the defect that in the existing ARMA technology, photovoltaic generation power super-short-term prediction accuracy is low is overcome, and high-precision photovoltaic power generation super-short-term prediction is achieved.

Description

Photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology

Technical field

The present invention relates to photovoltaic power electric powder prediction in generation of electricity by new energy process, relate to particularly a kind of arma modeling photovoltaic power ultra-short term Forecasting Methodology of photometry network real time correction.

Background technology

The large-scale new forms of energy base majority that China's photovoltaic generation produces after entering the large-scale development stage is positioned at " three northern areas of China " (northwest, northeast, North China); large-scale new forms of energy base is generally away from load center, and its electric power need to be transported to load center and dissolve through long-distance, high voltage.Due to intermittence, randomness and the undulatory property of wind, light resources, cause wind-powered electricity generation, the photovoltaic generation in extensive new forms of energy base to be exerted oneself fluctuation in a big way can occur thereupon, further cause the fluctuation of power transmission network charge power, bring series of problems to safe operation of electric network.

By in April, 2014, photovoltaic generation installed capacity has reached 4,350,000 kilowatts, accounts for 13% of Gansu Power Grid total installation of generating capacity, and simultaneously Gansu becomes China's photovoltaic generation largest province of installing.At present, Gansu Power Grid wind-powered electricity generation, photovoltaic generation installation exceed 1/3 of Gansu Power Grid total installation of generating capacity.Along with improving constantly of new-energy grid-connected scale, photovoltaic generation uncertainty and uncontrollability are brought problems to the safety and stability economical operation of electrical network.Accurately estimating available generating light resources is the basis to large-scale photovoltaic generating Optimized Operation.Photovoltaic generation power in photovoltaic generation process is predicted, be can be that generation of electricity by new energy Real-Time Scheduling, generation of electricity by new energy are planned, generation of electricity by new energy monthly plan, generation of electricity by new energy capability evaluation and abandon optical quantum and estimate to provide key message a few days ago.

ARMA (autoregressive moving-average model) is widely used in the prediction of photovoltaic generation power ultra-short term as a kind of machine learning method of maturation.Arma modeling is made up of autoregressive model (AR) and moving average model (MA), adopts to historical power is carried out autoregression computing and white noise sequence carried out to the photovoltaic generation that running mean comes in predict future 0-4 hour and exert oneself.ARMA method has many good qualities, therefore be widely used in the prediction of photovoltaic generation power ultra-short term, but the shortcoming of ARMA maximum is exactly the hysteresis quality of its prediction---, when photovoltaic generation is exerted oneself while changing, the pace of change of the result of ARMA prediction is generally slower than actual photovoltaic generated output pace of change.Therefore, have a strong impact on the precision of prediction of ARMA.

Summary of the invention

The object of the invention is to, for the problems referred to above, propose a kind of photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology, to realize the advantage of high precision photovoltaic generation power ultra-short term prediction.

For achieving the above object, the technical solution used in the present invention is:

A kind of photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology, comprises that input data obtain autoregressive moving-average model parameter;

Input light resources monitoring system data and operation monitoring system data, and according to operational monitoring real-time correction start capacity;

Thereby setting up autoregressive moving-average model obtains photovoltaic power ultra-short term and predicts the outcome;

Introducing light simultaneous measurement station data predicts the outcome and carries out real time correction photovoltaic power ultra-short term.

According to a preferred embodiment of the invention, described input data obtain autoregressive moving-average model parameter and comprise, input model training basic data;

Model is determined rank;

Adopt square method of estimation to estimate ARMA (p, the q) model parameter of determining rank.

According to a preferred embodiment of the invention, described input model training basic data, input data comprise, historical radiation data and historical power data.

According to a preferred embodiment of the invention, described model is determined rank and is specially:

Adopt residual error variogram method to carry out model and determine rank, be specially and establish x _tfor the item that needs are estimated, x _t-1, x _t-2..., x _t-nfor known historical power sequence, for ARMA (p, q) model, the value of Model Parameter p and q determined rank and determines by model;

The models fitting original series increasing progressively gradually with serial exponent number all calculates residual sum of squares (RSS) at every turn then draw exponent number and figure, when exponent number is during by little increase, can significantly decline, reach after true exponent number value can tend towards stability gradually, even on the contrary increase,

The observed value item number of actual observed value number actual use while referring to model of fit, for the sequence with N observed value, matching AR (p) model, the actual observed value using mostly is N-p most, model parameter number refers to the actual number of parameters comprising in set up model, and for the model that contains average, model parameter number is that model order adds 1, for the sequence of N observed reading, the residual error estimator of arma modeling is:

Wherein, the sum of squares function that Q is error of fitting, and θ _j(1≤j≤q) is model coefficient, and N is observation sequence length, it is the constant term in model parameter.

According to a preferred embodiment of the invention, described employing square method of estimation estimates that to determining ARMA (p, the q) model parameter on rank concrete steps are:

Historical photovoltaic plant power data is utilized to data sequence x ₁, x ₂..., x _trepresent, its sample autocovariance is defined as

${\widehat{\mathrm{\γ}}}_k=\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k},$

Wherein, k=0,1,2 ..., n-1, x _tand x _t-kbe data sequence x ₁, x ₂..., x _tin numerical value;

? ${\widehat{\mathrm{\γ}}}_0=\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2$

Historical power data sample autocorrelation function is:

${\widehat{\mathrm{\ρ}}}_k=\frac{{\widehat{\mathrm{\γ}}}_k}{{\widehat{\mathrm{\γ}}}_0}=\frac{\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2}=\frac{\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2},$

Wherein, k=0,1,2 ..., n-1;

The square of AR part is estimated as,

Order

Covariance function is

With estimation replace γ _k,

Can obtain parameter

To MA (q) model coefficient θ ₁, θ ₂..., θ _qemploying square estimates at

${\mathrm{\γ}}_0\left(y_t\right)=(1+{\mathrm{\θ}}_1^2+{\mathrm{\θ}}_2^2+...+{\mathrm{\θ}}_q^2){\mathrm{\σ}}_a^2$ Until

${\mathrm{\γ}}_k\left(y_t\right)=(-{\mathrm{\θ}}_k+{{\mathrm{\θ}}_1\mathrm{\θ}}_{k+1}+...+{\mathrm{\θ}}_{q-k}{\mathrm{\θ}}_q){\mathrm{\σ}}_a^2$

Wherein k=1,2 ..., m,

Above m+1 equation nonlinear equation, adopts process of iteration to solve and obtains autoregressive moving-average model parameter.

According to a preferred embodiment of the invention,

Described light resources monitoring system data comprise the light simultaneous measurement data that monitor at the photometry station relevant to photovoltaic plant to be predicted and the photovoltaic plant average radiation of numerical weather forecast data prediction, described operation monitoring system data are photovoltaic components in photovoltaic plant Real-Time Monitoring information to be predicted, comprise that photovoltaic DC-to-AC converter stops machine status information in real time.

According to a preferred embodiment of the invention, also comprise,

To predict the outcome and export in database, and show by chart and curve the contrast that predicts the outcome and show prediction and measured result.

According to a preferred embodiment of the invention, described autoregressive moving-average model is:

Wherein, and θ _j(1≤j≤q) is coefficient, α _tit is white noise sequence.

According to a preferred embodiment of the invention, described introducing light simultaneous measurement station data predict the outcome and carry out real time correction and be specially photovoltaic power ultra-short term:

If t ₁in the moment, it is I that the photovoltaic plant average irradiance obtaining is monitored at photometry station ₁, the photovoltaic plant average irradiance of data of weather forecast prediction is J ₁, the actual of photovoltaic plant exerted oneself as p ₁; Next time point t ₂in the moment, the photovoltaic plant average irradiance of data of weather forecast prediction is J ₂, the actual irradiance I of photovoltaic plant ₂for,

I ₂＝I ₁+(J ₂-J ₁)

The parameter correction of predicting power of photovoltaic plant is

$k=(\frac{I_2-I_1}{I_1}\×100\%)\×p_1$ (I ₁≠ 0 o'clock).

According to a preferred embodiment of the invention, output finally predict the outcome into:

Wherein, X _tthe prediction of exerting oneself of t moment photovoltaic plant, and θ _j(1≤j≤q) is coefficient, α _tbe white noise sequence, λ is weighting coefficient, I _tit is the average irradiance of t moment photovoltaic plant.

Technical scheme of the present invention has following beneficial effect:

Technical scheme of the present invention is predicted by the photovoltaic generation power in photovoltaic generation process, and by introducing light simultaneous measurement station data, photovoltaic generation power ultra-short term is predicted the outcome and carries out real time correction, overcome the low defect of photovoltaic generation power ultra-short term precision of prediction in existing ARMA technology, reach the object of high-precision photovoltaic generation power ultra-short term prediction.

Below by drawings and Examples, technical scheme of the present invention is described in further detail.

Brief description of the drawings

Fig. 1 is the photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology theory diagram described in the embodiment of the present invention.

Embodiment

Below in conjunction with accompanying drawing, the preferred embodiments of the present invention are described, should be appreciated that preferred embodiment described herein, only for description and interpretation the present invention, is not intended to limit the present invention.

A kind of photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology, comprises that input data obtain autoregressive moving-average model parameter;

Input light resources monitoring system data and operation monitoring system data, and according to operational monitoring real-time correction start capacity;

Thereby setting up autoregressive moving-average model obtains photovoltaic power ultra-short term and predicts the outcome;

Introducing light simultaneous measurement station data predicts the outcome and carries out real time correction photovoltaic power ultra-short term.

As shown in Figure 1, the photovoltaic generation power ultra-short term prediction that technical solution of the present invention proposes can be divided into two stages: model training stage and power prediction stage.

Stage 1: model training

Step 1.1: model training basic data input

Photovoltaic generation power prediction system model training required input data comprise historical radiation data, historical power data etc.Basic data is input to and in forecast model, carries out model training.

Step 1.2: model is determined rank

Need to set up estimation function with the item of how many known time sequences owing to cannot determining in advance, so need to carry out determining rank judgement to model.

If x _tfor the item that needs are estimated, x _t-1, x _t-2..., x _t-nfor known historical power sequence, for ARMA (p, q) model, it is exactly the value of determining Model Parameter p and q that model is determined rank.

Adopt residual error variogram method to carry out model and determine rank.Hypothetical model is limited rank autoregressive models, if the exponent number arranging is less than true exponent number, be a kind of not enough matching, thereby matching residual sum of squares (RSS) must be bigger than normal, now can significantly reduce residual sum of squares (RSS) by improving exponent number.Otherwise, if exponent number has reached actual value, increase so again exponent number, be exactly overfitting, now increase exponent number and can not make residual sum of squares (RSS) significantly reduce, even can slightly increase.

Carry out matching original series with the model that a series of exponent numbers increase progressively gradually like this, all calculate residual sum of squares (RSS) at every turn then draw exponent number and figure.When exponent number is during by little increase, can significantly decline, reach after true exponent number value can tend towards stability gradually, sometimes even on the contrary increase.The estimator of residual error variance is:

The observed value item number of actual use when " actual observed value number " refers to model of fit, for the sequence with N observed value, matching AR (p) model, the actual observed value using mostly is N-p most.

" model parameter number " refers to the actual number of parameters comprising in set up model, and for the model that contains average, model parameter number is that model order adds 1.For the sequence of N observed reading, the residual error estimator of corresponding arma modeling is:

(formula 1)

Step 1.3: model parameter estimation

Adopt square method of estimation to estimate the model parameter of ARMA (p, q).First, historical photovoltaic plant power data is utilized to data sequence x ₁, x ₂..., x _trepresent, its sample autocovariance is defined as

${\widehat{\mathrm{\γ}}}_k=\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}$ (formula 2)

Wherein, k=0,1,2 ..., n-1, x _tand x _t-kbe data sequence x ₁, x ₂..., x _tin numerical value.

Especially,

${\widehat{\mathrm{\γ}}}_0=\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2$ (formula 3)

Historical power data sample autocorrelation function is:

${\widehat{\mathrm{\ρ}}}_k=\frac{{\widehat{\mathrm{\γ}}}_k}{{\widehat{\mathrm{\γ}}}_0}=\frac{\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2}=\frac{\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2}$ (formula 4)

Wherein, k=0,1,2 ..., n-1.

The square of AR part is estimated as

(formula 5)

Order

(formula 6)

Covariance function is

(formula 7)

With estimation replace γ _k, have

(formula 8)

Can obtain parameter

To MA (q) model coefficient θ ₁, θ ₂..., θ _qemploying square estimates at

${\mathrm{\γ}}_0\left(y_t\right)=(1+{\mathrm{\θ}}_1^2+{\mathrm{\θ}}_2^2+...+{\mathrm{\θ}}_q^2){\mathrm{\σ}}_a^2$ (formula 9)

……

……

……

${\mathrm{\γ}}_k\left(y_t\right)=(-{\mathrm{\θ}}_k+{{\mathrm{\θ}}_1\mathrm{\θ}}_{k+1}+...+{\mathrm{\θ}}_{q-k}{\mathrm{\θ}}_q){\mathrm{\σ}}_a^2$ (formula 10)

Wherein k=1,2 ..., m.

Below comprise altogether m+1 equation, for its parameter, equation is non-linear, adopts process of iteration to solve.

Concrete steps are as follows, and equation is deformed into:

${\mathrm{\σ}}_a^2{=\mathrm{\γ}}_0/(1+{\mathrm{\θ}}_1^2+{\mathrm{\θ}}_2^2+...+{\mathrm{\θ}}_q^2)$ (formula 11)

${\mathrm{\θ}}_k=-\frac{{\mathrm{\γ}}_k}{{\mathrm{\σ}}_a^2}+{\mathrm{\θ}}_1{\mathrm{\θ}}_{k+1}+...+{\mathrm{\θ}}_{q-k}{\mathrm{\θ}}_q,$ K=1,2 ..., m (formula 12)

Given θ ₁, θ ₂..., θ _qwith one group of initial value, as

${\mathrm{\θ}}_1={\mathrm{\θ}}_2=...={\mathrm{\θ}}_q=0,{\mathrm{\σ}}_a^2={\mathrm{\γ}}_0$ (formula 13)

The above two formula the right of substitution, the value that the left side obtains is first step iterative value, is designated as by the right side of two formulas in this value successively substitution, just obtain second step iterative value again, the like, until adjacent twice iteration result while being less than given threshold value, got the result of gained as the approximate solution of parameter.

Find by above-mentioned solution procedure, solve the exponent number of time series models, will obtain seasonal effect in time series predicted value; Obtain seasonal effect in time series predicted value, must first set up concrete anticipation function; Set up concrete anticipation function, must know the exponent number of model.

According to experiment, time series models exponent number is generally no more than 5 rank.So in the time of this algorithm specific implementation, first hypothesized model is 1 rank, utilize method for parameter estimation in step 1.3 to obtain the parameter of first order modeling, and then set up estimation function and just can estimate to obtain each predicted value in the hope of first order modeling time series models, thereby try to achieve the residual error variance of first order modeling; Afterwards, hypothesized model is second order, tries to achieve the residual error of second-order model with said method; By that analogy, can obtain the residual error of 1 to 5 rank model, select the exponent number of model of residual error minimum as the exponent number of final mask.Determine after model order, just can calculate parameter θ ₁, θ ₂..., θ _qvalue.

Stage 2: power prediction

Step 2.1: light resources monitoring system data input

Light resources monitoring system data mainly comprise the photovoltaic plant average radiation of the light simultaneous measurement data that monitor at the photometry station relevant to photovoltaic plant to be predicted and NWP (numerical weather forecast data) prediction.

Step 2.2: operation monitoring system data input

Operation monitoring system data refer to photovoltaic components in photovoltaic plant Real-Time Monitoring information to be predicted, mainly comprise that photovoltaic DC-to-AC converter stops the status informations such as machine situation in real time.

Step 2.3: operational monitoring real-time correction start capacity

In photovoltaic plant operational process; because a variety of causes can cause shutdown situation; such as the photoelectricity station of typical 50,000 kilowatts of installations; 50,000 kilowatts of the average common less thaies of start capacity; therefore can know the start capacity of photovoltaic plant reality by real-time photovoltaic operational monitoring data; predict but not carry out photovoltaic generation power ultra-short term with the installed capacity of photovoltaic plant, thereby obtain higher precision of prediction.

Step 2.4: the photovoltaic generation power ultra-short term prediction based on arma modeling

By model parameter estimation out after, in conjunction with the model order of having estimated, just can obtain the time series equation for photovoltaic generation power ultra-short term prediction.The p drawing according to above-mentioned steps and q value, and θ ₁, θ ₂..., θ _qvalue set up autoregressive moving-average model;

Autoregressive moving-average model is as follows:

(formula 14)

Wherein, and θ _j(1≤j≤q) is coefficient, α _tit is white noise sequence.

Step 2.5: monitoring resource real-time correction photovoltaic generation power ultra-short term predicts the outcome

Can find out by above-mentioned ARMA forecast model, for the real-time change of photovoltaic generation power, above-mentioned model always has hysteresis quality, and the present invention predicts the outcome and carries out real time correction photovoltaic generation power ultra-short term by introducing light simultaneous measurement station data.

If t ₁in the moment, it is I that the photovoltaic plant average irradiance obtaining is monitored at photometry station ₁, the photovoltaic plant average irradiance of NWP prediction is J ₁, the actual of photovoltaic plant exerted oneself as p ₁; Next time point t ₂in the moment, the photovoltaic plant average irradiance of NWP prediction is J ₂, the actual irradiance I of photovoltaic plant ₂for,

I ₂=I ₁+ (J ₂-J ₁) (formula 15)

The parameter correction of predicting power of photovoltaic plant is

$k=(\frac{I_2-I_1}{I_1}\×100\%)\×p_1$ (I ₁≠ 0 o'clock) (formula 16)

Step 2.6: output and displaying finally predict the outcome

The arma modeling photovoltaic generation power ultra-short term of photometry network real time correction predict the outcome into

(formula 17)

Wherein, X _tthe prediction of exerting oneself of t moment photovoltaic plant, and θ _j(1≤j≤q) is coefficient, α _tbe white noise sequence, λ is weighting coefficient, I _tit is the average irradiance of t moment photovoltaic plant.

By introducing the revised prediction irradiance of photometry station Real-time Monitoring Data, can make weighting adjustment to next step prediction of arma modeling, thereby solve the hysteresis quality problem of arma modeling prediction.

To predict the outcome and export in database, and show by chart and curve the contrast that predicts the outcome, shows prediction and measured result.

Finally it should be noted that: the foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, although the present invention is had been described in detail with reference to previous embodiment, for a person skilled in the art, its technical scheme that still can record aforementioned each embodiment is modified, or part technical characterictic is wherein equal to replacement.Within the spirit and principles in the present invention all, any amendment of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.

Claims (10)

1. a photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology, is characterized in that, comprises that input data obtain autoregressive moving-average model parameter;

Input light resources monitoring system data and operation monitoring system data, and according to operational monitoring real-time correction start capacity;

Thereby setting up autoregressive moving-average model obtains photovoltaic power ultra-short term and predicts the outcome;

Introducing light simultaneous measurement station data predicts the outcome and carries out real time correction photovoltaic power ultra-short term.

2. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 1, is characterized in that, described input data obtain autoregressive moving-average model parameter and comprise, input model training basic data;

Model is determined rank;

Adopt square method of estimation to estimate ARMA (p, the q) model parameter of determining rank.

3. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 2, is characterized in that, described input model training basic data, and input data comprise, historical radiation data and historical power data.

4. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 3, is characterized in that, described model is determined rank and is specially:

Adopt residual error variogram method to carry out model and determine rank, be specially and establish x _tfor the item that needs are estimated, x _t-1, x _t-2..., x _t-nfor known historical power sequence, for ARMA (p, q) model, the value of Model Parameter p and q determined rank and determines by model;

The models fitting original series increasing progressively gradually with serial exponent number all calculates residual sum of squares (RSS) at every turn then draw exponent number and figure, when exponent number is during by little increase, can significantly decline, reach after true exponent number value can tend towards stability gradually, even on the contrary increase,

The observed value item number of actual observed value number actual use while referring to model of fit, for the sequence with N observed value, matching AR (p) model, the actual observed value using mostly is N-p most, model parameter number refers to the actual number of parameters comprising in set up model, and for the model that contains average, model parameter number is that model order adds 1, for the sequence of N observed reading, the residual error estimator of arma modeling is:

Wherein, the sum of squares function that Q is error of fitting, and θ _j(1≤j≤q) is model coefficient, and N is observation sequence length, it is the constant term in model parameter.

5. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 4, is characterized in that, described employing square method of estimation estimates that to determining ARMA (p, the q) model parameter on rank concrete steps are:

Historical photovoltaic plant power data is utilized to data sequence x ₁, x ₂..., x _trepresent, its sample autocovariance is defined as

${\widehat{\mathrm{\γ}}}_k=\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k},$

Wherein, k=0,1,2 ..., n-1, x _tand x _t-kbe data sequence x ₁, x ₂..., x _tin numerical value;

? ${\widehat{\mathrm{\γ}}}_0=\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2$

Historical power data sample autocorrelation function is:

${\widehat{\mathrm{\ρ}}}_k=\frac{{\widehat{\mathrm{\γ}}}_k}{{\widehat{\mathrm{\γ}}}_0}=\frac{\frac{1}{n}\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\frac{1}{n}\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2}=\frac{\underset{t=k+1}{\overset{n}{\mathrm{\Σ}}}{x}_t{x}_{t-k}}{\underset{t=1}{\overset{n}{\mathrm{\Σ}}}{x}_t^2},$

Wherein, k=0,1,2 ..., n-1;

The square of AR part is estimated as,

Order

Covariance function is

With estimation replace γ _k,

Can obtain parameter

To MA (q) model coefficient θ ₁, θ ₂..., θ _qemploying square estimates at

${\mathrm{\γ}}_0\left(y_t\right)=(1+{\mathrm{\θ}}_1^2+{\mathrm{\θ}}_2^2+...+{\mathrm{\θ}}_q^2){\mathrm{\σ}}_a^2$ Until

${\mathrm{\γ}}_k\left(y_t\right)=(-{\mathrm{\θ}}_k+{{\mathrm{\θ}}_1\mathrm{\θ}}_{k+1}+...+{\mathrm{\θ}}_{q-k}{\mathrm{\θ}}_q){\mathrm{\σ}}_a^2$

Wherein k=1,2 ..., m,

Above m+1 equation nonlinear equation, adopts process of iteration to solve and obtains autoregressive moving-average model parameter.

6. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 5, is characterized in that,

Described light resources monitoring system data comprise the light simultaneous measurement data that monitor at the photometry station relevant to photovoltaic plant to be predicted and the photovoltaic plant average radiation of numerical weather forecast data prediction, described operation monitoring system data are photovoltaic components in photovoltaic plant Real-Time Monitoring information to be predicted, comprise that photovoltaic DC-to-AC converter stops machine status information in real time.

7. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 6, is characterized in that, also comprise,

To predict the outcome and export in database, and show by chart and curve the contrast that predicts the outcome and show prediction and measured result.

8. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 7, is characterized in that, described autoregressive moving-average model is:

Wherein, and θ _j(1≤j≤q) is coefficient, α _tit is white noise sequence.

9. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 8, is characterized in that, described introducing light simultaneous measurement station data predict the outcome and carry out real time correction and be specially photovoltaic power ultra-short term:

If the t1 moment, it is I1 that the photovoltaic plant average irradiance obtaining is monitored at photometry station, and the photovoltaic plant average irradiance of data of weather forecast prediction is J ₁, the actual of photovoltaic plant exerted oneself as p ₁; Next time point t ₂in the moment, the photovoltaic plant average irradiance of data of weather forecast prediction is J ₂, the actual irradiance I of photovoltaic plant ₂for,

I ₂＝I ₁+(J ₂-J ₁)

The parameter correction of predicting power of photovoltaic plant is

$k=(\frac{I_2-I_1}{I_1}\×100\%)\×p_1$ (I ₁≠ 0 o'clock).

10. photometry network real time correction arma modeling photovoltaic power ultra-short term Forecasting Methodology according to claim 9, is characterized in that, output finally predict the outcome into:

Wherein, X _tthe prediction of exerting oneself of t moment photovoltaic plant, and θ _j(1≤j≤q) is coefficient, α _tbe white noise sequence, λ is weighting coefficient, I _tit is the average irradiance of t moment photovoltaic plant.