Summary of the invention
The purpose of this invention is to provide a kind of method that adopts sparse coefficient autoregressive model prediction gas, concrete is exactly to carry out the gas density early warning prediction that suddenlys change unusually according to the field monitoring data characteristic of gas density in the gas tunnel, overcomes the influence of the unsafe factor that ex post forecasting brings.
The method of the employing sparse coefficient autoregressive model prediction gas among the present invention, be to lay computer in ground central station, data communication interface and substation, face in the tunnel, lining cutting, add the broadband and monitoring point, return air inlet place is laid the face sensor respectively, the lining cutting sensor, add broadband sensor and return air inlet sensor, described four monitoring points are laid left arch springing sensor respectively, right arch springing sensor and dome sensors, all the sensors is connected with substation by communication cable, substation is connected with data communication interface by the transmission data wire, data communication interface is connected with computer, it is characterized in that, comprise the steps:
(1) in described computer, sets up the sparse coefficient autoregression model;
(2) described various sensor is transformed into the signal of telecommunication to the gas density that is monitored and is transferred to described substation;
(3) described substation carries out analyzing and processing respectively to the gas density data that difference detects the position, simultaneously gas density is passed to described computer;
(4) computer is drawn the time graph figure according to gas density, determines 3 important abnormity point positions in the curvilinear figure: unusual starting point, slowly rise to the slow rising of fast rise/be climbed to/slowly drop to fast rise/drop to fast turning point, the peak point of slow decline;
(5) computer is determined the parameter of sparse coefficient autoregression technology model according to described time graph figure and abnormity point position, calls sparse coefficient autoregression technology model trend prediction is carried out in the unusual sudden change of gas density in future.
At gas Forecasting Methodology sparse coefficient autoregression model of the present invention:
In the formula
Subscript be 1,2,3,4 ... the subclass of P}, P makes BIC criterion function BIC
P=log (σ
ε(P))
2+ [Plog (N-n)
2The value of]/(N-n) (n 〉=P 〉=0) extreme value minimum, this moment, model was best forecasting model.
The present invention is by carrying out the gas prediction to the gas monitor data, by sparse coefficient autoregression forecast model the unusual early warning of greatly may carrying out constantly of the outstanding preceding appearance of gas is predicted, provide reference frame for formulating control gas technical measures, simultaneously, guarantee constructor's safety, therefore, the present invention is significant to the gas tunnel construction safety.
The specific embodiment
Disclosed all features in this manual, or the step in disclosed all methods or the process except mutually exclusive feature and/or step, all can make up by any way.
Disclosed arbitrary feature in this manual (comprising any accessory claim, summary and accompanying drawing) is unless special narration all can be replaced by other equivalences or the alternative features with similar purpose.That is, unless special narration, each feature is an example in a series of equivalences or the similar characteristics.
Preferred embodiment: sparse coefficient autoregression forecast model
The principle of this model is:
If temporal model is
Y
t=α
0+α
1Y
t-1+α
2Y
t-2+…+α
pY
t-p+ε
t(1)
In the formula: p---model order; ε
t---the model residual error.
In order to try to achieve the α in (1) formula
0, α
1, α
2..., α
p, must at first set up matrix and find the solution.Concrete grammar is as follows:
At first make p=n (n is the maximum amount of lagging, i.e. the model top step number).Can get following equation group according to the principle of least square:
In the formula:
Wherein:
After solving Al from (2) formula, related parameter arranged in must (1) formula:
α
0=Y-(α
1Y
1+α
2Y
2+…+α
nY
n)
Selecting variable in order to the equation (1) that last method is obtained without BIC criterion, also is not the time series forecasting equation of an optimum.For trying to achieve optimal models, must set up bordered matrix S.Make S
Ij=B
Ij, X
i=B
I0(1≤i, j≤n)
X herein
0=1.Following formula is carried out stepwise regression analysis, and find out with BIC criterion and to make the minimum model of criterion function value as best forecasting model.
So-called BIC criterion is exactly a value of calculating following formula, makes it the standard of accepting or rejecting as model.
BIC
P=log(σ
ε(P))
2+[Plog(N-n)
2]/(N-n)(n≥P≥0)
(δ in the formula
ε(P))
2Match residual error variance for P rank model.Concrete treatment step to (3) formula is as follows:
At first (3) formula is carried out continuous n cancellation conversion, can get preceding n element X of n+1 row after the conversion
1, X
2... X
n, these are exactly the parameter of n rank standardized model.Final reduction formula is:
(σ
ε(n))
2=A·Y
0/(N-n)
BIC
n=log(σ
ε(n))
2+[n·log(N-n)
2]/(N-n)
A
JjThe i row of expression s-matrix to the A of strain capacity j in the A matrix
JjValue, Y
0It is the lower right corner element of the s-matrix after the conversion.
In order to delete the parameter item less to regression effect contribution, existing with following formula as the standard of whether deleting:
Get
Following formula has illustrated that the regression effect of in variable is the poorest, should delete it.The method of deletion is that the in the matrix capable (row) is transferred to outermost layer, drag then with (n=1, n=1) cancellation conversion, result of variations with regard to cancellation the variable of desire deletion, and preceding n-1 element X being listed as of matrix n
1, X
2..., X
N-1Be exactly the parameter of that model of the residual sum of squares (RSS) minimum of n-1 rank standardized model, (n, n) element is exactly the residual sum of squares (RSS) of this model to S.Calculate the BIC value of n-1 rank model, and calculate P
N-2 (i)(1≤i≤n-1), the decision deletion is what variable, and then carries out the adjustment of matrix ranks, carries out (n, n) cancellation conversion.And the like, until the zeroth order model.Relatively the BIC value of each rank model is minimum as the BIC value of P rank model, just determine that so the P rank are best model, relevant parameters is
And model is
In the formula
Subscript normally 1,2,3,4 ... the subclass of P}.
Here the monitored data with certain tunnel tunnel face sensor is an example, utilizes the dynamic data processing method of sparse coefficient autoregression forecast that its " the unusual very big possibility that occurs constantly " carried out analyzing and processing.This face sensor passes to substation to the gas density data that monitored, in substation, gas density is carried out analyzing and processing with analysis software, simultaneously analysis processing result is formed gas density time graph figure in computer, computer gathers the initial data in the curvilinear figure as follows automatically:
Computer gathers the initial data in the curvilinear figure as follows automatically:
The sequence data of after original data processing, setting up:
i |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
Y
i |
6 |
12 |
2 |
7 |
9 |
12 |
9 |
8 |
4 |
5 |
6 |
15 |
i |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
21 |
22 |
|
|
Y
i |
11 |
6 |
2 |
6 |
15 |
12 |
11 |
2 |
3 |
4 |
|
|
In order to make time series data meet the characteristics of sparse coefficient model, time interval Y
iThe unit of sequence is the sky.
Get n=9, by computer program to Y
iTime series is carried out analyzing and processing, model order P and corresponding BIC value result:
P |
BIC |
P |
BIC |
9 |
0.3584E+01 |
4 |
0.2466E+01 |
8 |
0.3207E+01 |
3 |
0.2657E+01 |
7 |
0.2831E+01 |
2 |
0.2754E+01 |
6 |
0.2459E+01 |
1 |
0.2494E+01 |
5 |
0.2178E+01 |
0 |
0.2942E+01 |
What the BIC value was minimum is the autoregression model of P=5.Determine model parameter
Be respectively: 41.702 ,-0.953 ,-1.167 ,-1.659 ,-0.792,0.334; Variable of selecting and corresponding order: 3,9,6,8,1.So forecasting model is:
Final prognostic equation is:
In order to check the effect of prognostic equation, below done concrete prediction:
According to predicted value, the 23rd time gas density takes place is unusually after 7.121 days, and (be Y this 25 days October in 2007 with reality
23=6) the gas density maximum value is 1.34% (Y
24=15) compare, precision of forecasting model is reliable.Therefore be feasible unusually with sparse coefficient autoregressive model prediction gas.
The present invention is not limited to the aforesaid specific embodiment.The present invention expands to any new feature or any new combination that discloses in this manual, and the arbitrary new method that discloses or step or any new combination of process.