基于回歸分析的北京市天然氣供應(yīng)規(guī)模預(yù)測(cè)

摘 要

摘要：通過(guò)對(duì)北京市1997年以來(lái)天然氣年供應(yīng)量關(guān)于GDP、年末房屋建筑面積的回歸分析，得出線性回歸方程，預(yù)測(cè)天然氣年供應(yīng)量。由建立的管道總長(zhǎng)度關(guān)于天然氣年供應(yīng)量的一元線性回
摘要：通過(guò)對(duì)北京市1997年以來(lái)天然氣年供應(yīng)量關(guān)于GDP、年末房屋建筑面積的回歸分析，得出線性回歸方程，預(yù)測(cè)天然氣年供應(yīng)量。由建立的管道總長(zhǎng)度關(guān)于天然氣年供應(yīng)量的一元線性回歸方程，預(yù)測(cè)管網(wǎng)規(guī)模。
關(guān)鍵詞：GDP；年末房屋建筑面積；回歸分析；負(fù)荷預(yù)測(cè)；年供氣量；管網(wǎng)規(guī)模
Forecasting of Natural Gas Supply Scale in Beijing City Based on Regression Analysis
LIU Yan，ZHANG Ying-hui，DU Xue-ping
Abstract：Based on the regression analysis between annual natural gas supply volume and GDP/year-end housing construction area since 1 997 in Beijing City，the linear regression equation is obtained to predict the annual natural gas supply volume. The scale of pipeline network is predicted by the simple linear regression equation between total pipeline length and annual natural gas supply volume.
Key words：gross domestic product(GDP)；year-end housing construction area； regression analysis；load forecasting；annual gas supply volume；scale of pipeline network
1 概述
    北京市天然氣市場(chǎng)仍將在未來(lái)一段時(shí)期內(nèi)保持快速增長(zhǎng)。2008年北京市燃?xì)饧瘓F(tuán)天然氣供應(yīng)量約為49×108m3，高峰流量為196×104m3/h，2020年有望突破200×108m3的供應(yīng)規(guī)模。陜京三線及唐山LNG正在規(guī)劃建設(shè)中，計(jì)劃給北京供氣。處于天然氣產(chǎn)業(yè)鏈中上游的長(zhǎng)輸管道和LNG氣源地建設(shè)，可能會(huì)與下游持續(xù)增長(zhǎng)的負(fù)荷需求不相適應(yīng)[1、2]，中上游配套設(shè)施規(guī)劃建設(shè)相對(duì)滯后，會(huì)導(dǎo)致城市天然氣氣源保障不足，類似于2004年冬季供氣緊張的局面可能再度出現(xiàn)。所以，采用多種方法準(zhǔn)確預(yù)測(cè)未來(lái)負(fù)荷需求對(duì)可靠供氣有重要意義，本文采用一元線性回歸方法預(yù)測(cè)天然氣的供應(yīng)規(guī)模。
2 線性回歸方法預(yù)測(cè)天然氣年供應(yīng)量
    城市天然氣輸配管網(wǎng)的年輸送量與當(dāng)年的GDP、年末房屋建筑面積等因素存在一定關(guān)系，可以對(duì)這種關(guān)系進(jìn)行回歸分析。回歸分析指從某一組數(shù)據(jù)出發(fā)，判斷變量之間是否存在相關(guān)關(guān)系，確定它們之間的數(shù)學(xué)表達(dá)式，此數(shù)學(xué)表達(dá)式即為回歸方程，并對(duì)其可信程度作統(tǒng)計(jì)檢驗(yàn)，最終利用確定的回歸方程進(jìn)行預(yù)測(cè)和控制[3、4]。將GDP、年末房屋建筑面積分別作為自變量，天然氣年供應(yīng)量作為因變量，通過(guò)回歸分析預(yù)測(cè)某一階段的供氣量及管道總長(zhǎng)度。
   ① 天然氣年供應(yīng)量關(guān)于GDP的回歸分析
   1997年—2008年北京市的GDP與天然氣年供應(yīng)量見(jiàn)表1，對(duì)應(yīng)關(guān)系見(jiàn)圖1，從圖1可知天然氣年供應(yīng)量與GDP存在線性關(guān)系。
    線性回歸理論一直被廣泛應(yīng)用，在許多實(shí)際問(wèn)題中，當(dāng)自變量局限于一定范圍時(shí)，可以取線性模型作為真實(shí)模型的近似，其誤差可忽略[3]。
表1 1997年—2008年北京市GDP與天然氣年供應(yīng)量
年份
1997
1998
1999
2000
2001
2002
GDP/元
2076×108
2376×108
2678×108
3161×108
3711×108
4330×108
天然氣年供應(yīng)量/m3
1.8×108
3.8×108
7.6×108
11.0×108
16.7×108
20.5×108
年份
2003
2004
2005
2006
2007
2008
GDP/元
5024×108
6060×108
6886×1088
7861×108
9353×108
10488×108
天然氣年供應(yīng)量/m3
23.8×108
27.0×108
31.7×108
38.4×108
43.4×108
60.1×108
 

   北京市天然氣年供應(yīng)量與GDP基本呈線性關(guān)系，可用式(1)表示。
    V=aX+b    (1)
式中X——年GDP，元
    V——年天然氣供應(yīng)量，m3
    a、b——回歸系數(shù)
   根據(jù)表1數(shù)據(jù)，得出回歸方程：
   V=0.0061422X-8.9436×108    (2)
   建立經(jīng)驗(yàn)回歸方程的目的在于揭示天然氣年供應(yīng)量與GDP之間的內(nèi)在規(guī)律，所建立的回歸方程是否有意義需要經(jīng)過(guò)回歸顯著性檢驗(yàn)。用F檢驗(yàn)對(duì)回歸方程作顯著性檢驗(yàn)的方法稱為方差分析，檢驗(yàn)線性回歸方程顯著性的方差分析[4]見(jiàn)表2。
表2 線性回歸方程顯著性方差分析
方差來(lái)源
偏差平方和
自由度
方差
F
Fα
顯著性
回歸
SR
1
VR=SR/1
F=VR/Ve
F0.05(1,n-2)
 
剩余
Se
n-2
Ve=Se(n-2)
 
F0.01(1,n-2)
 
總和
ST
n-1
 
   SR：回歸平方和，反映了X變化引起V波動(dòng)的大小，是通過(guò)X對(duì)V的相關(guān)性引起的。
   Se：剩余平方和，反映了除X對(duì)V的線性影響之外的剩余因素引起V波動(dòng)的大小。
    ST：總的偏差平方和，它反映V1，V2，…，Vn的離散程度，ST=SR+5。。
    若F＞F0.01(1,n-2)，則顯著性為高度顯著，用**表示；若F0.05(1,n-2)≤F≤F0.01(1,n-2)，則顯著性為顯著，用術(shù)表示；若F＜F0.05(1,n-2)，則顯著性為不顯著，不作標(biāo)記。
    若回歸方程有意義，即V波動(dòng)主要是由X變化引起的，其他因素是次要的，即要求SR盡可能大，而SR盡可能小。回歸直線對(duì)樣本數(shù)據(jù)點(diǎn)的擬合程度，通常用可決系數(shù)R2來(lái)表示，R2=SR/ST，R2越接近于1，表明回歸直線對(duì)樣本數(shù)據(jù)點(diǎn)的擬合程度越高[3、4]，這種檢驗(yàn)方法稱為R檢驗(yàn)。1997年—2008年北京市天然氣年供應(yīng)量與GDP方差分析見(jiàn)表3。
表3 1997年—2008年北京市天然氣年供應(yīng)量與GDP方差分析
方差來(lái)源
偏差平方和
自由度
方差
F
Fα
顯著性
回歸
3279.0
1
3279.0
408.8
5.0
**
剩余
80.2
10
8.0
 
10.0
 
總和
3359.2
11
 
    可決系數(shù)R2=0.97612，標(biāo)準(zhǔn)誤差Sy=2.8322。F＞F0.01(1,10)，故顯著性為高度顯著。
    天然氣年供應(yīng)量關(guān)于GDP的回歸方程通過(guò)了F檢驗(yàn)與R檢驗(yàn)，即可用來(lái)預(yù)測(cè)。預(yù)測(cè)是在確定自變量為某一個(gè)值時(shí)求因變量的估計(jì)值，分為點(diǎn)預(yù)測(cè)和區(qū)間預(yù)測(cè)。
    假設(shè)在2008年后12年內(nèi)，GDP保持每年10.9%的增幅，預(yù)測(cè)得到2009年—2020年天然氣的供應(yīng)量，見(jiàn)表4。經(jīng)驗(yàn)證，2009年的實(shí)際天然氣供應(yīng)量在預(yù)測(cè)天然氣供應(yīng)量的置信度為95%的預(yù)測(cè)區(qū)間內(nèi)。
    回歸方程表明GDP每增加1000×108元，天然氣年供應(yīng)量增加約6.14×108m3。隨著能源結(jié)構(gòu)的調(diào)整，今后燃?xì)獍l(fā)電的比例將有很大提高，雖然以GDP為自變量的回歸分析能反映未來(lái)的發(fā)展趨勢(shì)，但是對(duì)1997年—2008年基礎(chǔ)數(shù)據(jù)的分析并不能體現(xiàn)能源政策調(diào)整力度加大所帶來(lái)的結(jié)果。以GDP為自變量的分析表明：按1997年—2008年的發(fā)展模式，到2020年天然氣供應(yīng)量將達(dá)到約210×108m3。新清潔能源政策與既有發(fā)展模式的共同作用，將導(dǎo)致2020年天然氣供應(yīng)量為210×108m3與新增電廠用氣量之和。
表4 GDP年增幅為10.9%時(shí)的天然氣供應(yīng)量預(yù)測(cè)值
年份
預(yù)測(cè)天然氣供應(yīng)量/m3
預(yù)測(cè)天然氣供應(yīng)量的置信度為95%的預(yù)測(cè)區(qū)間/m3
2009
62.497×108
56.33×108～68.67×108
2010
70.284×108
64.11×108～76.46×108
2011
78.920×108
72.75×108～85.09×108
2012
88.497×108
82.33×108～94.67×108
2013
99.118×108
92.95×108～105.29×108
2014
110.897×108
104.73×108～117.07×108
2015
123.960×108
117.79×108～130.13×108
2020
213.999×108
207.83×108～220.17×108
    ② 天然氣年供應(yīng)量關(guān)于建筑面積的回歸分析
    天然氣年供應(yīng)量與GDP的相關(guān)性分析表明一元線性回歸分析可以用于預(yù)測(cè)天然氣供應(yīng)量。根據(jù)北京市近幾年的用氣結(jié)構(gòu)特點(diǎn)推斷北京市天然氣供應(yīng)量與年末房屋建筑面積也存在相關(guān)關(guān)系。1997年—2007年北京市天然氣年供應(yīng)量與年末房屋建筑面積見(jiàn)表5，并對(duì)這些數(shù)據(jù)進(jìn)行一元線性回歸分析。
表5 1997年—2007年北京市天然氣年供應(yīng)量與年末房屋建筑面積
年份
年末房屋建筑面積/m2
天然氣年供應(yīng)量/m3
1997
29 294×104
1.8×108
1998
31019×104
3.8×108
1999
33273×104
7.6×108
2000
35531×104
11.0×108
2001
37717×104
16.7×108
2002
40483×104
20.5×108
2003
43122×104
23.8×108
2004
46523×104
27.0×108
2005
50707×104
31.7×108
2006
54282×104
38.4×108
2007
57361×104
43.4×108
注：數(shù)據(jù)來(lái)自1998年—2008年北京統(tǒng)計(jì)年鑒。
回歸分析得到回歸方程：
    V=14.525A-40.1318×108    (3)
式中A——年末房屋建筑面積，m2
    可決系數(shù)R2=0.99214，標(biāo)準(zhǔn)誤差，Sy=1.2956，顯著性為高度顯著。
    年末房屋建筑面積為自變量的回歸分析結(jié)果表明天然氣年供應(yīng)量與年末房屋建筑面積存在線性關(guān)系，所得回歸方程通過(guò)了F檢驗(yàn)和R檢驗(yàn)，可用于預(yù)測(cè)。年末房屋建筑面積每增加1×108m2，天然氣年供應(yīng)量需增加約14.53×108m3，可決系數(shù)最接近于1，所以建筑面積應(yīng)作為天然氣需求預(yù)測(cè)的重點(diǎn)關(guān)注對(duì)象。若北京市年末房屋建筑面積達(dá)到15×108m3，則天然氣供應(yīng)量將達(dá)到178×108m3。與利用GDP的分析類似，以年末房屋建筑面積為自變量的回歸分析，部分反映了未來(lái)的發(fā)展趨勢(shì)，清潔能源政策與既有發(fā)展模式的共同作用將導(dǎo)致未來(lái)天然氣供應(yīng)量為178×108m3與新增電廠用氣量之和。以GDP和年末房屋建筑面積為自變量的回歸分析表明北京市2020年天然氣供應(yīng)量有望達(dá)到200×108m3。
    ③ 藍(lán)天數(shù)關(guān)于天然氣年供應(yīng)量的回歸分析
    我國(guó)發(fā)展城市燃?xì)獾闹饕康氖歉纳拼髿猸h(huán)境[5]，表6所示為2000年—2008年北京市天然氣年供應(yīng)量與藍(lán)天數(shù)，回歸分析表明天然氣年供應(yīng)量與藍(lán)天數(shù)具有相關(guān)關(guān)系。天然氣在北京的利用確實(shí)達(dá)到了改善大氣環(huán)境的目的。
    分析得到回歸方程：
    S=1.9528×10-8V+164.6477    (4)
式中S——藍(lán)天數(shù)，即空氣質(zhì)量二級(jí)和好于二級(jí)的天數(shù)，d
   可決系數(shù)R2=0.91285，標(biāo)準(zhǔn)誤差Sy=9.7066，顯著性為高度顯著。
表6 2000年—2008年北京市藍(lán)天數(shù)與天然氣年供應(yīng)量
年份
2000
2001
2002
2003
2004
天然氣年供應(yīng)量/m3
11×108
17×108
20×108
24×108
27×108
藍(lán)天數(shù)/d
177
185
203
224
229
年份
2005
2006
2007
2008
 
天然氣年供應(yīng)量/m3
32×108
38 ×108
43×108
60×108
 
藍(lán)天數(shù)/d
234
241
246
274
 
3 天然氣管網(wǎng)規(guī)模預(yù)測(cè)
    ① 管道總長(zhǎng)度的回歸分析
    研究表明，2030年之前我國(guó)天然氣產(chǎn)業(yè)將處于成長(zhǎng)期[6]，仍然以城市供氣為主要方向。管網(wǎng)建設(shè)與負(fù)荷的增加存在相關(guān)關(guān)系[7]，1997年—2007年北京市管道總長(zhǎng)度與天然氣年供應(yīng)量見(jiàn)表7，分析得到如下的回歸方程：
    L=149.8919×10-8V+2027.1271 (5)
式中L——天然氣管道總長(zhǎng)度，km
   可決系數(shù)R2=0.91092，標(biāo)準(zhǔn)誤差Sy=685.178，顯著性為高度顯著。
    回歸系數(shù)149.8919的意義是天然氣年供應(yīng)量每增加1×108m3，需要增建天然氣管道約149.9km。
    ② 中壓管道總長(zhǎng)度的回歸分析
    截至2008年，北京天然氣管網(wǎng)次高壓A以上管道長(zhǎng)度僅占管道總長(zhǎng)度的8.6%，且大多處于城市外圍，合理增建高壓管道具備條件。中壓管道長(zhǎng)度已占管道總長(zhǎng)度的40%，中壓管道幾乎遍布市區(qū)所有主干道路，形成了密集的環(huán)狀管網(wǎng)。中壓管道建設(shè)是連接用戶和高壓轉(zhuǎn)輸平臺(tái)的關(guān)鍵，除電廠以外的其他用戶幾乎都要經(jīng)過(guò)中壓管道的轉(zhuǎn)輸，中壓管道的完善是用戶發(fā)展的關(guān)鍵。1997年—2007年北京市中壓管道總長(zhǎng)度與天然氣年供應(yīng)量(不包括電廠用氣量)見(jiàn)表8。
表7 1997年—2007年北京市管道總長(zhǎng)度與天然氣年供應(yīng)量
年份
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
天然氣年供應(yīng)量/m3
1.8×108
3.8×108
7.6×108
11.0×108
16.7×108
20.5×108
23.8×108
27.0×108
31.7×108
38.4×108
43.4×108
管道總長(zhǎng)度/km
1222
2182
3277
3924
5019
5435
6464
6756
7004
7295
7551
表8 1997年—2007年北京市中壓管道總長(zhǎng)度與天然氣年供應(yīng)量(不包括電廠用氣量)
年份
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
不包括電廠用氣量的天然氣年供應(yīng)量/m3
1.8×108
3.8×108
7.6×108
11.0×108
16.7×108
20.5×108
23.8×108
27.0×108
31.7×108
35.7×108
40.7×108
中壓管道總長(zhǎng)度/km
609
698
839
1011
1212
1696
2045
2394
2639
2782
2912
分析得到回歸方程：
    L′=66.3613×10-8V+383.419    (6)
式中L′——天然氣中壓管道長(zhǎng)度，km
   可決系數(shù)R2=0.97209，標(biāo)準(zhǔn)誤差Sy=155.2057，顯著性為高度顯著。
    回歸系數(shù)66.3613的意義是通過(guò)中壓管道的天然氣年供應(yīng)量每增加1×108m3，需要增建中壓管道約66.4km。北京市平均每建設(shè)1km中壓管道的費(fèi)用約為150×104元，則每增加1×108m3的天然氣供應(yīng)量需要投資約1×108元。
4 結(jié)論
    隨著北京市城市建設(shè)的發(fā)展，郊區(qū)縣作為城市多中心進(jìn)行建設(shè)開(kāi)發(fā)，燃?xì)庑枨蠛腿細(xì)庠O(shè)施的建設(shè)必將達(dá)到新的水平。但是，各用氣點(diǎn)較為分散，管道建設(shè)量較大。由市中心區(qū)歷史數(shù)據(jù)得到的回歸方程可作為郊區(qū)縣管網(wǎng)規(guī)劃的參考。天然氣年供應(yīng)量與GDP、年末房屋建筑面積線性相關(guān)，可以得出一般性結(jié)論，用于預(yù)測(cè)臨近年份基本因素的變動(dòng)對(duì)總供氣量的作用程度。目前樣本數(shù)據(jù)有限，更精確的預(yù)測(cè)模型有待于樣本數(shù)據(jù)的增加完善。天然氣供應(yīng)規(guī)模受能源政策調(diào)整的影響，線性模型可以作為真實(shí)模型的近似，但要考慮能源政策的影響，對(duì)預(yù)測(cè)結(jié)果進(jìn)行修正。
參考文獻(xiàn)：
[1] 劉毅軍.天然氣產(chǎn)業(yè)鏈下游市場(chǎng)風(fēng)險(xiǎn)研究[M].北京：石油工業(yè)出版社，2007.
[2] 肖久明.城市燃?xì)庳?fù)荷的特點(diǎn)與預(yù)測(cè)模型的特征[J].煤氣與熱力，2005，25(1)：67-70.
[3] 何燦芝，羅漢.應(yīng)用數(shù)理統(tǒng)計(jì)[M].長(zhǎng)沙：湖南大學(xué)出版社，2004.
[4] 王巖，隋思連，王愛(ài)青.數(shù)理統(tǒng)計(jì)與MATLAB工程數(shù)據(jù)分析[M].北京：清華大學(xué)出版社，2006.
[5] 田貫三，付林.“西氣東輸”中天然氣合理應(yīng)用方式研究[M].北京：中國(guó)建筑工業(yè)出版社，2009.
[6] 章兆淇.國(guó)外天然氣發(fā)展的基本特點(diǎn)和啟示[J].國(guó)際石油經(jīng)濟(jì)，2001，(1)：20-25.
[7] 李猷嘉.燃?xì)廨斉湎到y(tǒng)的設(shè)計(jì)與實(shí)踐[M].北京：中國(guó)建筑工業(yè)出版社，2007.
 
(本文作者：劉燕1 張應(yīng)輝1 杜學(xué)平2 1.北京市燃?xì)饧瘓F(tuán)有限責(zé)任公司 北京 100035；2.北京建筑工程學(xué)院 北京 100044)