李洪兵 劉盈 唐紫涵 韓咪 羅雨 劉可

（1.四川師范大學(xué)工學(xué)院，四川 成都 610101；2.中國石油集團共享運營有限公司成都中心，四川 成都 610051；3.西南石油大學(xué)經(jīng)濟管理學(xué)院，四川 成都 610500)

0 引言

中國能源消耗量現(xiàn)居世界首位，碳排放量不斷增加［1］。2020 年，習(xí)近平總書記提出了“碳達(dá)峰、碳中和”（以下簡稱“雙碳”）目標(biāo)。隨著“雙碳”目標(biāo)的實施，天然氣已成為能源低碳轉(zhuǎn)型期最具消費潛力的重要化石能源。在實現(xiàn)碳達(dá)峰前，天然氣消費量持續(xù)上升［2-5］，未來天然氣消費快速反彈的現(xiàn)象將涌現(xiàn)。隨著“十四五”能源規(guī)劃的推出，加速了能源結(jié)構(gòu)優(yōu)化升級的進程，天然氣已成為中國能源低碳轉(zhuǎn)型期不可替代的重要能源［6］。天然氣相較于煤炭、石油，其碳排放量更低且比電能、煤炭、石油效率更高［7］。與此同時，在“一帶一路”和新型城鎮(zhèn)化的推動下，天然氣在能源結(jié)構(gòu)中的地位顯著提高，天然氣將成為未來能源發(fā)展的主流，天然氣行業(yè)具有巨大的發(fā)展空間［8］。天然氣在一次能源消費中占比持續(xù)提高，供不應(yīng)求現(xiàn)象明顯，致使中國天然氣供需矛盾加劇。全球天然氣供需市場格局正悄然發(fā)生變化，國際天然氣貿(mào)易的不確定性日漸凸顯。能源安全與社會穩(wěn)定和經(jīng)濟發(fā)展息息相關(guān)，天然氣安全又是能源安全不可或缺的部分，影響天然氣安全的核心是天然氣供需狀態(tài)的變化。因此為保障能源安全，維護社會穩(wěn)定和經(jīng)濟健康持續(xù)發(fā)展，開展天然氣供需預(yù)測預(yù)警研究迫在眉睫，基于灰色系統(tǒng)理論對中國天然氣供需進行預(yù)測預(yù)警研究十分必要。

1 理論方法與模型

1.1 GM（1，1）預(yù)測模型

GM（1，1）模型在灰色系統(tǒng)理論中不僅是最基本也是被運用最廣泛的一種預(yù)測模型。該預(yù)測模型主要應(yīng)用于解決存在指數(shù)變化規(guī)律的數(shù)據(jù)預(yù)測。優(yōu)點是不需要大量的數(shù)據(jù)就可以對中短期進行預(yù)測，可以解決因數(shù)據(jù)不全和可靠性較低的問題，且精確度高、易于檢驗、計算簡便。計算步驟如下：已知數(shù)據(jù)X(0)＝X(0)(1 )，X(0)( 2 )，…，X(0)(n)

（1）對已知數(shù)據(jù)進行級比檢驗

如果所有的級比λ(k)都落在內(nèi)，即符合模型灰色預(yù)測初等條件。反之則需要對序列x(0)做必要的變換處理，使其可以包含于θ內(nèi)。

（2）根據(jù)方程式x(0)(k)＋az(1)(k)＝b計算a和b的值

計算已知數(shù)據(jù)的一階方程

（4）得到非齊次函數(shù)

其中a為GM（1，1）模型的發(fā)展系數(shù)，b為灰作用量。

1.2 預(yù)警評價方法與模型

1.2.1 預(yù)警指標(biāo)權(quán)重的確定方法

為實現(xiàn)預(yù)警指標(biāo)權(quán)重值準(zhǔn)確性提高的目標(biāo)，分別運用熵值法和序關(guān)系分析法從客觀和主觀兩個角度來確定預(yù)警指標(biāo)權(quán)重。熵值法是根據(jù)歷年指標(biāo)值的改變趨勢來確定各項指標(biāo)的權(quán)重，用該方法可以避免由于人為因素而帶來的誤差。但由于熵值法會因指標(biāo)細(xì)微的變化導(dǎo)致指標(biāo)權(quán)重發(fā)生大的改變，因此將熵值法與序關(guān)系分析法相結(jié)合計算綜合權(quán)重可提高準(zhǔn)確性。序關(guān)系分析法對于熵值法而言，序關(guān)系分析法不會由于評價指標(biāo)的變化而影響指標(biāo)間相對重要性程度。相對于層次分析法優(yōu)勢在于不需要建構(gòu)判斷矩陣且無需一致性檢查，操作性強。

1）熵值法

①選取數(shù)據(jù)：共有m個樣本數(shù)量，n個因素指標(biāo)Yij表示第i個樣本中第j個指標(biāo)，其中i＝1，2，3，...m；j＝1，2，3，...，n。

②對數(shù)據(jù)標(biāo)準(zhǔn)化處理確保各指標(biāo)計量單位一致和方向統(tǒng)一。

正向指標(biāo)（越大越好的指標(biāo)）：

負(fù)向指標(biāo)（越小越好的指標(biāo)）

③計算樣本權(quán)重：

④計算第j項指標(biāo)的熵值：

⑤計算第j項指標(biāo)差異系數(shù)

⑥計算各項指標(biāo)

2）序關(guān)系分析法（G1法）

①根據(jù)熵值法計算的權(quán)重值，來確定各項指標(biāo)互相的序列關(guān)系：若指標(biāo)xi比xj重要則記作xi＞xj，按照這一原則，并結(jié)合行業(yè)專家提供的評價結(jié)果，確定評價指標(biāo)的重要性排序。

②確定相鄰指標(biāo)間的相對重要程度：設(shè)定專家對評價對象的評價指標(biāo)xk-1和xk的相對重要程度之比進行賦值：

其中n是天然氣安全的評價指標(biāo)總數(shù)，n＝9

③計算指標(biāo)權(quán)重：

最后通過G1法求得權(quán)重集合為wk＝(w1w2，…，wn)

3）綜合權(quán)重

綜合權(quán)重即是將熵值法所得權(quán)重與序關(guān)系分析法所得權(quán)重進行組合得出的綜合權(quán)重值。本研究將客觀權(quán)重與主觀權(quán)重視為同等重要，即各占50%，則通過加權(quán)計算出天然氣各預(yù)警指標(biāo)綜合權(quán)重值，計算公式為：wi＝0.5qj＋0.5wk。

1.2.2 灰色關(guān)聯(lián)綜合評價方法

灰色關(guān)聯(lián)分析法是通過子序列和母序列的幾何形狀類似程度來判斷其相互關(guān)系是否緊密，進而得到其安全度。灰色關(guān)聯(lián)分析法的優(yōu)點是對樣本或者樣本的數(shù)量是否具有規(guī)律性沒有特殊要求，且其量化結(jié)果和定性分析結(jié)果差異性較小，同時計算量小。計算步驟如下：

（1）確定子序列和母序列：有i個樣本數(shù)量和j個指標(biāo)， 每一個yi(j)(i＝1，2，3，…，m；j＝1，2，3，…，n)表示第i年第j個指標(biāo)值。y0(j)為每個指標(biāo)的最優(yōu)值。

子序列為：

母序列為：

（2）將各個指標(biāo)數(shù)據(jù)進行無量綱化處理：

無量綱化子序列：

無量綱化母序列為：

（3）計算比較序列中每個元素與參考序列相對應(yīng)元素的差的絕對值、最大值和最小值：

（4）確定各項指標(biāo)的權(quán)重值，已從客觀和主觀角度求得各項指標(biāo)的綜合權(quán)重。

（5）計算灰色關(guān)聯(lián)系數(shù)：

其中ξ為分辨系數(shù)，通常取值為0.5；δmin＝

（6）計算安全度：

利用灰色關(guān)聯(lián)系數(shù)與各預(yù)警指標(biāo)綜合權(quán)重的乘積表示測度，將所得預(yù)警測度值轉(zhuǎn)變到（0，10］范圍，計算公式和轉(zhuǎn)換公式分別如下。

2 中國天然氣供需預(yù)測預(yù)警分析

2.1 預(yù)警指標(biāo)體系建構(gòu)

在能源結(jié)構(gòu)轉(zhuǎn)型進程的大力推動下，中國天然氣能源供需便成為關(guān)注的焦點。為充分了解中國天然氣未來供需發(fā)展?fàn)顟B(tài)，通過對目前已有權(quán)威文獻(xiàn)的梳理分析［9-16］，并結(jié)合當(dāng)前中國國情，從供給、需求、市場、社會4個方面，選取儲采比、自給率、儲量指數(shù)、消費占比、消費強度、進口依存度、進口集中度、保障系數(shù)、碳強度及人均GDP 共10 個指標(biāo)，構(gòu)建中國天然氣供需預(yù)測預(yù)警指標(biāo)體系，如表1所示。

表1 中國天然氣供需預(yù)警指標(biāo)體系表

其中指標(biāo)包括儲量指數(shù)、保障系數(shù)的原始數(shù)據(jù)均來自《bp世界能源統(tǒng)計年鑒（2012-2021）》；消費占比以及人均GDP 均來源于國家統(tǒng)計局；國內(nèi)生產(chǎn)總值、天然氣產(chǎn)量、天然氣消費量及能源消費總量的數(shù)據(jù)均來自于國家統(tǒng)計局；碳年均排放總量、天然氣前3名國家進口量、天然氣總進口量的數(shù)據(jù)均來自《bp 世界能源統(tǒng)計年鑒（2022）》；中國年底可采儲量、天然氣采出量數(shù)據(jù)來源于《中國統(tǒng)計年鑒（2019-2021）》，自給率、消費強度、進口集中度、保障系數(shù)、碳強度由于未找到原始數(shù)據(jù)來源，將已知數(shù)據(jù)帶入表1中計算公式中計算得出數(shù)據(jù)。

2.2 預(yù)測分析

根據(jù)GM（1，1）模型的基本原理，運用matlab 軟件繪制出中國天然氣供需預(yù)測預(yù)警指標(biāo)的預(yù)測值，結(jié)果如圖1 所示。根據(jù)平均相對誤差絕對值（MAPE）公式：，可計算得出儲采比的MAPE 為6%≤10%、自給率的MAPE 為2%≤10%、儲量指數(shù)MAPE 為4%≤10%、消費占比MAPE 為3%≤10%、消費強度MAPE 為3%≤10%、進口依存度MAPE 為4%≤10%、進口集中度MAPE 為3%≤10%、保障系數(shù)的MAPE 為1%≤10%、碳強度MAPE為2%≤10%、人均GDP的MAPE為3%≤10%，表明GM（1，1）模型預(yù)測效果為優(yōu)［9］，預(yù)測結(jié)果可以用于預(yù)警綜合分析。

圖1 預(yù)警指標(biāo)趨勢變化圖

由圖1分析可知，中國天然氣儲采比和儲量指數(shù)2022—2030 年將呈指數(shù)增長；中國天然氣自給率在“十四五”和“十五五”期間均將表現(xiàn)為線性降低趨勢；中國天然氣消費強度及消費占比將有穩(wěn)定攀升的趨勢，在2030 年前中國天然氣消費量將持續(xù)增長。在新冠肺炎疫情席卷全球的背景下，2019 和2020 年中國天然氣進口依存度和進口集中度出現(xiàn)驟降點，造成全球貿(mào)易市場經(jīng)濟相對滯后從而影響天然氣市場安全。后新冠疫情時期，天然氣市場出現(xiàn)回暖現(xiàn)象。隨著國內(nèi)經(jīng)濟的穩(wěn)定發(fā)展、能源轉(zhuǎn)型的不斷推動、新型城鎮(zhèn)化的實行，中國人均GDP 和中國天然氣碳強度在“十五五”期間，前者表現(xiàn)為震蕩式增長后者為震蕩式下降，表明相關(guān)政策的實施將取得一定的成效。

2.3 預(yù)警分析

2.3.1 預(yù)警界限與等級劃分

根據(jù)熵值法對天然氣供需預(yù)測預(yù)警指標(biāo)進行客觀計算權(quán)重，采用序關(guān)系分析法求各指標(biāo)主觀權(quán)重，最后計算出綜合權(quán)重如表2所示：中國天然氣供需預(yù)測預(yù)警指標(biāo)權(quán)重較大的4 個因子依次為儲采比、自給率、儲量指數(shù)、進口集中度；意味著該四個因子對中國天然氣供需影響相較其余指標(biāo)影響更大。

表2 預(yù)警指標(biāo)權(quán)重表

在已有研究成果基礎(chǔ)上，結(jié)合本次研究的預(yù)警指標(biāo)權(quán)重，分析得出天然氣供需各預(yù)警指標(biāo)預(yù)警界線和預(yù)警等級如表3、表4所示。

表3 預(yù)警界線表

表4 預(yù)警等級標(biāo)準(zhǔn)

2.3.2 預(yù)警綜合評價分析

由圖2分析可知，中國天然氣供需安全度將呈緩慢上升態(tài)勢，將從臨界安全狀態(tài)變?yōu)榘踩珷顟B(tài)，預(yù)警等級將從Ⅲ級轉(zhuǎn)變?yōu)棰艏?。結(jié)果表明：在2022—2030 年，中國天然氣供給將從臨界安全狀態(tài)上升至安全狀態(tài)，隨著油氣增儲上產(chǎn)持續(xù)加強，中國天然氣供給預(yù)警等級將由Ⅲ級上升至Ⅳ級；隨著能源低碳轉(zhuǎn)型的持續(xù)推進，中國天然氣需求量將不斷增長，中國天然氣需求安全等級將維持在極度危險狀態(tài)，預(yù)警等級將處于Ⅱ級；中國天然氣市場安全度雖有所下降，但在未來仍處于臨界安全狀態(tài)，預(yù)警等級將達(dá)到Ⅲ級；中國天然氣社會將處于安全狀態(tài)，雖然安全度將呈現(xiàn)出平緩下降趨勢，但其預(yù)警等級仍保持在Ⅳ級。

圖2 天然氣供需安全度趨勢變化圖

3 結(jié)論

從天然氣供給、需求、市場、社會4個維度選擇10 個預(yù)警指標(biāo)構(gòu)建天然氣供需預(yù)測預(yù)警指標(biāo)體系，基于GM（1，1）模型對中國天然氣供需未來發(fā)展態(tài)勢進行預(yù)測；構(gòu)建灰色關(guān)聯(lián)預(yù)警模型對中國天然氣供需狀態(tài)進行預(yù)警分析，得出以下結(jié)論：

1）中國天然氣供需總體狀態(tài)發(fā)展水平較為平穩(wěn)，2012—2028 年中國天然氣供需狀況處于臨界安全狀態(tài)，預(yù)警等級處于Ⅲ級，2029—2030 年中國天然氣供需預(yù)警等級提升至Ⅳ級，表現(xiàn)為安全狀態(tài)。不斷促進清潔新能源和可再生能源的使用，將是提高中國天然氣供需預(yù)警等級的有效措施。

2）供給安全度穩(wěn)步上升，從臨界安全狀態(tài)發(fā)展到安全狀態(tài)。國家能源局提出油氣增儲上產(chǎn)七年行動計劃，中國天然氣探明儲量和產(chǎn)量將持續(xù)增長［17-20］，中國天然氣供應(yīng)能力將進一步提升。而天然氣消耗量的增長對天然氣供給能力會造成負(fù)面影響，天然氣消費力度持續(xù)增強促進了中國天然氣供需矛盾不斷加深。由此得出正確引導(dǎo)中國天然氣消費、增大天然氣開采和提高中國天然氣自給率是保障中國天然氣供給安全的重要措施。

3）需求安全度不斷下降，狀態(tài)屬性由臨界安全狀態(tài)轉(zhuǎn)變?yōu)椴话踩珷顟B(tài)。中國經(jīng)濟高質(zhì)量發(fā)展和“煤改氣”政策的實施，刺激了天然氣消費強度和消費占比持續(xù)升高。中國天然氣消費占比不斷提升，表明如今中國經(jīng)濟發(fā)展愈來愈依靠天然氣能源，也從側(cè)面反映中國天然氣的供應(yīng)壓力越來越大，致使2022—2030 年天然氣需求安全度也將不斷降低。因此提倡能源多元化、發(fā)掘其他新能源、降低天然氣消耗量減輕天然氣供應(yīng)壓力、提升能源供需比等將是提高中國天然氣需求預(yù)警等級和縮小天然氣供需缺口的有效方法。

4）受到國際地緣政治影響，國際局勢的變化擾亂了國際天然氣能源發(fā)展的穩(wěn)定狀態(tài)，致使市場安全形勢呈下降狀況，從安全狀態(tài)降至臨界安全狀態(tài)，預(yù)警等級由Ⅱ級轉(zhuǎn)變?yōu)棰蠹墶ｋm然中國天然氣進口集中度在降低，但其進口依存度仍在上升，造成中國天然氣市場總體預(yù)警等級處于較低狀態(tài)。進口依存度持續(xù)升高正面反映了中國天然氣進口依然受到其他國家的制約。所以，保障中國天然氣進口的多元化，選擇穩(wěn)定、安全的國家進口，降低中國天然氣進口集中度，改變天然氣進口鏈單一化現(xiàn)狀以及激勵性開采頁巖氣、煤層氣等資源是規(guī)避中國天然氣供需風(fēng)險的又一重要舉措。

（5）社會安全度有緩慢降低趨勢，維持在安全狀態(tài)，預(yù)警等級仍保持Ⅳ級。煤炭、石油等化石能源的大量使用，導(dǎo)致大氣污染現(xiàn)狀加重和碳排放量持續(xù)升高；人口老齡化問題不斷降低中國勞動力質(zhì)量進而影響國家經(jīng)濟發(fā)展，是導(dǎo)致社會安全度降低的重要因素。因此以減排為目的、改善人口質(zhì)量為宗旨、高質(zhì)量發(fā)展經(jīng)濟為核心，持續(xù)推進“氣改電”進程、不斷改善生育政策和高質(zhì)量穩(wěn)定發(fā)展中國經(jīng)濟以降低中國天然氣社會風(fēng)險；堅持“一帶一路”，不斷向“十四五”規(guī)劃和“雙碳”目標(biāo)奮進，將“藍(lán)天計劃”貫徹到底是保障天然氣社會安全的重要措施。