摘要：相較其他國際工程項目，跨國油氣管道項目的政治風(fēng)險更加突出。以中緬天然氣管道項目為例，以復(fù)雜網(wǎng)絡(luò)理論和平均場理論為基礎(chǔ)，引入信息風(fēng)險并運用改進(jìn)SIS模型對中緬天然氣管道項目的政治風(fēng)險擴(kuò)散特性進(jìn)行分析，為其他跨國油氣管道項目政治風(fēng)險管理和風(fēng)險擴(kuò)散應(yīng)對提供有益參考。

關(guān)鍵詞：跨國油氣管道；政治風(fēng)險；復(fù)雜網(wǎng)絡(luò)；改進(jìn)SIS模型

0 引言

隨著我國“一帶一路”倡議的深化和對外開放合作步伐的加快，我國海外油氣項目呈現(xiàn)井噴式增長，涉及石油、天然氣、頁巖氣等多個產(chǎn)業(yè)。其中，跨國油氣管道項目跨越多個國家，其政治、經(jīng)濟(jì)、自然環(huán)境較為復(fù)雜，不可控性較多，如何做好油氣管道項目政治風(fēng)險管理是企業(yè)和項目團(tuán)隊急需解決的問題。油氣管道政治風(fēng)險還會進(jìn)一步擴(kuò)散，對其他行為主體產(chǎn)生風(fēng)險影響。本文以中緬天然氣管道項目為研究對象，在改進(jìn)SIS模型基礎(chǔ)上引入多層級網(wǎng)絡(luò)，運用ORIGIN軟件對項目政治風(fēng)險擴(kuò)散進(jìn)行定性分析，為下一步油氣管道項目政治風(fēng)險擴(kuò)散應(yīng)對提供思路和方法，為項目安全建設(shè)提供有益參考。

目前，SIS模型在研究傳染病、信息及風(fēng)險擴(kuò)散傳播機(jī)制方面已經(jīng)進(jìn)行了運用。例如，Zhang等^[1]運用SIS反應(yīng)擴(kuò)散對流模型的漸近行為研究飽和和自發(fā)感染機(jī)制。周琦萍等^[2]運用SIS模型對輿情預(yù)警機(jī)制進(jìn)行研究分析。宋之杰等^[3]運用優(yōu)化SIS模型對共享經(jīng)濟(jì)消費傳染的內(nèi)在機(jī)理、擴(kuò)散過程和變化規(guī)律進(jìn)行研究。目前，缺乏運用SIS模型及多層級網(wǎng)絡(luò)對油氣管道項目政治風(fēng)險擴(kuò)散進(jìn)行研究的文獻(xiàn)。

1 中緬天然氣管道項目概況與政治風(fēng)險識別

1.1 項目概況

中緬天然氣管道項目起點位于緬甸若開邦皎漂蘭里島，從云南瑞麗入境我國，延伸至寶山、大理、昆明歸入貴州，最終達(dá)到廣西貴港，全長2520km，緬甸境內(nèi)管道全長793km，輸氣能力為120億m³/a。該項目于2013年7月28日正式投產(chǎn)，由6方組成合資公司進(jìn)行建設(shè)和運營。

1.2 政治風(fēng)險識別

以中緬天然氣管道項目為例進(jìn)行政治風(fēng)險擴(kuò)散研究，研究周期為2013年7月—2022年1月。這一時期為中緬天然氣管道項目的運營期，項目的政治風(fēng)險比較復(fù)雜和突出。

為確保準(zhǔn)確識別中緬天然氣管道項目的政治風(fēng)險，采用三類組合方法進(jìn)行分析。第一類是中緬天然氣管道項目內(nèi)部資料，中國石油集團(tuán)東南亞管道公司作為中緬天然管道項目建設(shè)運營主體擁有第一手風(fēng)險管理資料，包括政治風(fēng)險分析和處置等風(fēng)險研究資料；第二類是通過中英文數(shù)據(jù)庫檢索得到的研究資料^[4-11]；第三類是通過對各參建方，包括設(shè)計單位、施工單位、監(jiān)理單位、政府部門的溝通調(diào)研獲得的資料。通過上述方法列出該時期中緬天然氣管道項目面臨的主要政治風(fēng)險。中緬天然氣管道項目主要政治風(fēng)險辨識圖如圖1所示。

2 基于改進(jìn)SIS模型的中緬天然氣管道項目政治風(fēng)險擴(kuò)散研究

2.1 政治風(fēng)險擴(kuò)散模型構(gòu)建

SIS模型是一種典型的流行病傳播模型，其中S表示易感者，I表示感染者。SIS模型可以描述如下：最初，所有的節(jié)點都處于易感染狀態(tài)。然后，部分節(jié)點接觸到風(fēng)險后，變?yōu)楦腥緺顟B(tài)。這些感染狀態(tài)的節(jié)點繼續(xù)感染其他易感染狀態(tài)的節(jié)點，或者進(jìn)入恢復(fù)狀態(tài)；恢復(fù)狀態(tài)后，仍然有再次被感染進(jìn)而繼續(xù)參與風(fēng)險的傳播的機(jī)會^[12]。

本文將SIS模型引入中緬天然氣管道項目政治風(fēng)險擴(kuò)散模型，研究范圍為多層級能源環(huán)境網(wǎng)絡(luò)。將能源環(huán)境網(wǎng)絡(luò)分為三個層級，依次為項目層、企業(yè)層和政府層，每個層級都有不同的主體，風(fēng)險擴(kuò)散除了在同層進(jìn)行，也會在不同層進(jìn)行。同時，政治風(fēng)險通過信息媒介傳播也是一個關(guān)鍵的風(fēng)險擴(kuò)散模式。而信息導(dǎo)致的風(fēng)險能量增加也是政治風(fēng)險擴(kuò)散的重要載體和手段。本文基于改進(jìn)SIS模型的中緬天然氣管道項目政治風(fēng)險擴(kuò)散研究，不僅考慮擴(kuò)散網(wǎng)絡(luò)的層次性，也把政治信息風(fēng)險納入其中^[13-14]。

中緬天然氣管道項目政治風(fēng)險擴(kuò)散示意圖如圖2所示。假設(shè)政治風(fēng)險在項目層L₁開始擴(kuò)散，即政治風(fēng)險從中緬天然氣管道項目開始，經(jīng)過t時刻后，通過自身政治風(fēng)險的耦合、放大等效應(yīng)開始擴(kuò)散，使項目層L₁的能源項目被激活和感染的數(shù)量越來越多，且積聚的風(fēng)險能量越來越大，不再只局限于在項目層L₁進(jìn)行風(fēng)險擴(kuò)散，而會通過層間的風(fēng)險關(guān)聯(lián)的連邊擴(kuò)散至其他層級的風(fēng)險節(jié)點，如擴(kuò)散到企業(yè)層L₂，即可視為中緬天然氣管道項目政治風(fēng)險擴(kuò)散到了企業(yè)節(jié)點。與此同時，企業(yè)層L₂的風(fēng)險節(jié)點被激活后，也會在本層傳播擴(kuò)散并進(jìn)入到其他層如政府層L₃，同時也可能回到項目層L₁。

本文用感染率、感染密度來表征中緬天然氣管道項目的政治風(fēng)險擴(kuò)散過程。風(fēng)險能量除通過關(guān)聯(lián)關(guān)系（連邊）進(jìn)行擴(kuò)散，還會通過信息流進(jìn)行擴(kuò)散。信息流主要是通過網(wǎng)絡(luò)等媒介，將風(fēng)險信息（如負(fù)面信息、謠言、模糊信息、失真信息等）通過相應(yīng)載體傳遞到其他節(jié)點，這也是政治風(fēng)險被動放大的主要原因。因中緬天然氣管道項目為初始風(fēng)險源，故假定風(fēng)險源在L_P層，政治風(fēng)險擴(kuò)散動力學(xué)模型公式如下

式中，L_n表示中緬天然氣管道項目政治風(fēng)險擴(kuò)散網(wǎng)絡(luò)中第n層網(wǎng)絡(luò)；N_n表示第L_n層的節(jié)點總數(shù)；ρ^L_n表示第L_n層中被感染節(jié)點的密度；N^r_n表示第L_n層中感染節(jié)點的數(shù)量；P_nv表示第L_n層節(jié)點與第L_v層節(jié)點的連接概率；α_n表示第L_n層穩(wěn)定節(jié)點被感染的可能性；γ_n表示第L_n層風(fēng)險節(jié)點信息流感染的可能性，主要指信息渠道和信息資源本身的負(fù)效應(yīng)；β_n表示第L_n層被感染節(jié)點重新回到穩(wěn)定狀態(tài)的可能性；

多層級風(fēng)險網(wǎng)絡(luò)綜合擴(kuò)散速率μ=（α+γ）/β=α/β+γ/β；接觸擴(kuò)散速率v=α/β；信息擴(kuò)散速率ψ=γ/β。因此，在第L_n層的風(fēng)險擴(kuò)散速率μ_n=（α_n+γ_n）/β_n。

2.2 政治風(fēng)險擴(kuò)散模型分析

本文選取的中緬天然氣管道項目政治風(fēng)險擴(kuò)散范圍限定為項目所處的能源環(huán)境體系，涉及能源勘察、開發(fā)、生產(chǎn)、加工、轉(zhuǎn)換、運輸、分配、儲備、使用等各環(huán)節(jié)，由與中緬天然氣管道項目相關(guān)的不同能源國、能源企業(yè)、能源項目各主體組成，是相互間政治、經(jīng)濟(jì)、社會等關(guān)系要素的宏觀體現(xiàn)。中緬天然氣管道項目所在能源環(huán)境體系也處于整個國民經(jīng)濟(jì)大環(huán)境中，會不斷與其他環(huán)境進(jìn)行物質(zhì)、信息、資源等交換、整合和傳遞。

結(jié)合中緬天然氣管道項目資料、檢索文獻(xiàn)和專家溝通，將政治風(fēng)險擴(kuò)散網(wǎng)絡(luò)分為三個層級，項目層選取22個節(jié)點、企業(yè)層選取18個節(jié)點、政府層選取8個節(jié)點?；趯＜覇柧矸椒?，對中緬天然氣管道項目風(fēng)險關(guān)聯(lián)關(guān)系進(jìn)行分析確認(rèn)。所選專家為跨國油氣管道項目的專業(yè)人員和管理人員。若超過50%的專家認(rèn)為中緬天然氣管道項目與某一主體節(jié)點存在直接風(fēng)險關(guān)聯(lián)，則認(rèn)為兩個主體存在連邊；若低于50%，則認(rèn)為不存在直接風(fēng)險關(guān)聯(lián)，但是會通過中間節(jié)點存在風(fēng)險關(guān)聯(lián)^[15]。按照上述方法，發(fā)放10份專家調(diào)查問卷，對中緬天然氣管道項目各個節(jié)點的風(fēng)險關(guān)聯(lián)關(guān)系進(jìn)行確認(rèn)。若存在直接風(fēng)險關(guān)聯(lián)，則選1；若不存在直接風(fēng)險關(guān)聯(lián)，則選0。若超過50%的專家認(rèn)為存在關(guān)聯(lián)節(jié)點，則認(rèn)為節(jié)點存在直接風(fēng)險關(guān)聯(lián)；若低于50%，則認(rèn)為不存在直接風(fēng)險關(guān)聯(lián)。最后得到各層級風(fēng)險關(guān)聯(lián)矩陣。

根據(jù)改進(jìn)的SIS模型，求得穩(wěn)態(tài)解，公式如下

2.2.1 項目層政治風(fēng)險擴(kuò)散分析

以中緬天然氣管道項目為政治風(fēng)險擴(kuò)散的源頭，即項目層為擴(kuò)散的起點。假設(shè)信息擴(kuò)散的層內(nèi)擴(kuò)散概率均為0.1，風(fēng)險擴(kuò)散速率范圍為0～5，通過選取不同的項目層風(fēng)險鄰居節(jié)點數(shù)，進(jìn)行ORIGIN軟件數(shù)據(jù)分析。項目層政治風(fēng)險擴(kuò)散速率與風(fēng)險節(jié)點密度關(guān)系圖如圖3所示。

當(dāng)中緬天然氣管道項目發(fā)生政治風(fēng)險擴(kuò)散時，風(fēng)險擴(kuò)散速率與信息擴(kuò)散速率分別用x和y軸表示，項目層的風(fēng)險節(jié)點密度用z軸表示。由圖3可知，當(dāng)初始階段項目層的風(fēng)險擴(kuò)散速率與信息風(fēng)險擴(kuò)散速率都相較較小時，項目層的風(fēng)險節(jié)點密度也較??；隨著風(fēng)險擴(kuò)散速率增加，項目層的風(fēng)險節(jié)點密度呈現(xiàn)快速增長勢頭。當(dāng)項目層的風(fēng)險節(jié)點密度達(dá)到0.8時，風(fēng)險擴(kuò)散速率達(dá)到約4。隨著風(fēng)險擴(kuò)散速率增加，項目層的風(fēng)險節(jié)點密度變化相對平緩。

2.2.2 企業(yè)層、政府層風(fēng)險擴(kuò)散分析

由于中緬天然氣管道項目為政治風(fēng)險擴(kuò)散的初始節(jié)點，首先在項目層進(jìn)行擴(kuò)散，導(dǎo)致項目層風(fēng)險節(jié)點數(shù)量增加。當(dāng)項目層分別與企業(yè)層、政府層的最大風(fēng)險鄰居節(jié)點數(shù)為1～5時，項目層中的風(fēng)險感染節(jié)點通過層間連接和信息關(guān)系繼續(xù)擴(kuò)散到企業(yè)層和政府層節(jié)點，并導(dǎo)致其他層節(jié)點的風(fēng)險感染。

運用ORIGIN軟件對企業(yè)層和政府層的風(fēng)險節(jié)點密度與該層的政治風(fēng)險擴(kuò)散速率的關(guān)系進(jìn)行分析。層內(nèi)的政治風(fēng)險擴(kuò)散速率范圍為0～5，項目層與兩個層級的風(fēng)險鄰居節(jié)點數(shù)量1～5，風(fēng)險節(jié)點密度范圍為0～1。

為了便于比較項目層不同風(fēng)險鄰居節(jié)點對企業(yè)層和政府層的擴(kuò)散速率和層內(nèi)風(fēng)險節(jié)點密度的影響，假設(shè)企業(yè)層和政府層的信息擴(kuò)散速率與接觸擴(kuò)散速率相同，即""nbsp; ν=ψ""" ，將三維空間關(guān)系簡化為二維平面關(guān)系。企業(yè)層政治風(fēng)險擴(kuò)散速率與風(fēng)險節(jié)點密度關(guān)系圖如圖4所示，政府層政治風(fēng)險擴(kuò)散速率與風(fēng)險節(jié)點密度關(guān)系圖如圖5所示。

通過軟件計算發(fā)現(xiàn)，隨著企業(yè)層、政府層政治風(fēng)險擴(kuò)散速率增加，各層的風(fēng)險節(jié)點密度增加。當(dāng)層內(nèi)政治風(fēng)險擴(kuò)散速率接近4時，層內(nèi)風(fēng)險節(jié)點密度變化較小，趨于穩(wěn)定。當(dāng)企業(yè)層節(jié)點與項目層節(jié)點的最大風(fēng)險鄰居數(shù)量不斷增加時，各層的風(fēng)險節(jié)點密度不斷增加。項目層最大風(fēng)險節(jié)點越多，在相同的擴(kuò)散速率下，企業(yè)層和政府層的風(fēng)險節(jié)點密度也越大。隨著項目層最大風(fēng)險節(jié)點越大，企業(yè)層和政府層的風(fēng)險節(jié)點密度變化減小。

2.2.3 各層級政治風(fēng)險擴(kuò)散對比分析

結(jié)合ORIGIN軟件進(jìn)行數(shù)值模擬。當(dāng)項目層的最大風(fēng)險節(jié)點數(shù)為1、2、3時，三個層級政治風(fēng)險擴(kuò)散速率、風(fēng)險節(jié)點密度對比圖如圖6～圖8所示。

以中緬天然氣管道項目為風(fēng)險源，軟件計算結(jié)果顯示：當(dāng)各層風(fēng)險擴(kuò)散速率超過4時，三個層級的風(fēng)險節(jié)點密度變化趨于平穩(wěn)；當(dāng)項目層最大風(fēng)險節(jié)點增加時，三個層級在相同風(fēng)險擴(kuò)散速率下，層內(nèi)的風(fēng)險節(jié)點密度也相應(yīng)增加度。

在相同的項目層最大風(fēng)險節(jié)點和最大風(fēng)險擴(kuò)散速率下，三個層級中項目層的風(fēng)險節(jié)點密度最大，也即中緬天然氣管道項目率先對項目層節(jié)點進(jìn)行擴(kuò)散和感染，其次為政府層，風(fēng)險節(jié)點密度最小的為企業(yè)層。政府層如果與項目層的風(fēng)險連接概率越大，被政治風(fēng)險擴(kuò)散感染的概率也越高。而連接概率的高低與選取的節(jié)點數(shù)量和節(jié)點度有關(guān)。由本文選取的節(jié)點特性可知，項目層與企業(yè)層的連接概率高于項目層與政府層。

3 結(jié)語

基于改進(jìn)SIS模型對中緬天然氣管道項目風(fēng)險擴(kuò)散的研究，得出以下結(jié)論：

（1）風(fēng)險擴(kuò)散速率與風(fēng)險節(jié)點密度呈正相關(guān)關(guān)系。當(dāng)風(fēng)險擴(kuò)散速率和信息擴(kuò)散速率越大時，層內(nèi)的風(fēng)險節(jié)點密度隨之不斷增大。當(dāng)風(fēng)險擴(kuò)散速率和信息擴(kuò)散速率在4以上時，項目層、企業(yè)層、政府層的層內(nèi)風(fēng)險節(jié)點密度變化趨于平緩，即感染節(jié)點密度增長越來越慢。

（2）最大風(fēng)險鄰居節(jié)點與風(fēng)險節(jié)點密度呈正相關(guān)關(guān)系。當(dāng)項目層的風(fēng)險節(jié)點密度越大，即企業(yè)層、政府層的風(fēng)險鄰居節(jié)點越多時，企業(yè)層和政治層的風(fēng)險節(jié)點密度會不斷增長；而當(dāng)項目層的最大風(fēng)險鄰居節(jié)點很大時，企業(yè)層、政府層風(fēng)險節(jié)點密度趨于穩(wěn)定。即節(jié)點數(shù)越大，風(fēng)險節(jié)點密度的變化越小。

（3）層間連接度與風(fēng)險節(jié)點密度呈正相關(guān)關(guān)系。由于政府層與項目層連接概率比企業(yè)層與項目層更高，從平均場角度，政府層的風(fēng)險節(jié)點密度在相同條件下會大于企業(yè)層。

（4）信息擴(kuò)散風(fēng)險導(dǎo)致風(fēng)險節(jié)點密度增長更快。由于信息擴(kuò)散速率大于0，對于風(fēng)險節(jié)點密度公式而言是增函數(shù)，也就是在考慮了信息擴(kuò)散風(fēng)險的情況下，其風(fēng)險節(jié)點密度增長速率必然大于只考慮風(fēng)險接觸擴(kuò)散的情況。而信息擴(kuò)散速率越大，層內(nèi)風(fēng)險節(jié)點密度也會增長越快。

本文以多層級復(fù)雜網(wǎng)絡(luò)為中緬天然氣管道項目政治風(fēng)險擴(kuò)散的研究范圍，運用改進(jìn)SIS模型，引入信息擴(kuò)散風(fēng)險，以平均場理論為支撐，從宏觀視角對中緬天然氣管道項目的政治風(fēng)險擴(kuò)散特性進(jìn)行分析，厘清了政治風(fēng)險感染的相關(guān)影響因素，對天然氣管道項目政治風(fēng)險管理有一定的積極意義，也為開展類似項目政治風(fēng)險擴(kuò)散研究提供了有益參考。

參考文獻(xiàn)

[1]ZHANG J L，CUI R H.Asymptotic behavior of an SIS reaction-diffusion-advection model with saturation" and spontaneous infection mechanism[J].Zeitschrift Für Angewandte Mathematik Und Physik，2020，71（5）：1-21.

[2]周琦萍，楊芳.基于SIS模型的網(wǎng)絡(luò)輿情無監(jiān)督預(yù)警機(jī)制研究[J].情報科學(xué)，2019，37（8）：51-55.

[3]宋之杰，唐曉莉.基于SIS模型的共享經(jīng)濟(jì)消費傳染機(jī)制研究[J].數(shù)學(xué)的實踐與認(rèn)識，2020，50（1）：37-43.

[4]ZHAO H，KEEBENG O.The China-Myanmar energy pipelines：risks and benefits[J].Perspective，2013，5（30）：82-89.

[5]WU X，LIU X.The Myanmar-China oil and gas pipelines：is Myanmar dependable？[J].Geopolitics of Energy，2014，8（36）：2-5.

[6]TODD J.Energy and culture，the social construction of the Myanmar-China pipelines[D].Wellington：Victoria University of Wellington，2015.

[7]YANG Y，YIN Z Y，LEI S，et al.Impacts to Myanmar with China-Myanmar oil and gas pipeline project[C].Proceedings of the 2016 International Symposium on Business Cooperation and Development，2016.

[8]KULKARNI S S，NATHAN H S K.The elephant and the tiger：Energy security，geopolitics，and national strategy in China and India’s cross border gas pipelines[J].Energy Research amp; Social Science，2016（11）：183-194.

[9]ZHEN J，LU G S.Myanmar’s economic and political transition and chinese investment in Myanmar[J].Myanmar： Reintegrating into the International Community，2016，6（13）：245-265.

[10]DAWEI L，KENSUKE Y，HISASHI Y.Understanding the motivations behind the Myanmar-China energy pipeline：multiple streams and energy politics in China[J].Energy Policy，2017，1（107）：403-412.

[11]AUNG T S，LUAN S J，CONDON S.Evaluation of the environmental impact assessment （EIA） of Chinese EIA in Myanmar： Myitsone Dam，the lappadaung copper mine and the Sino-Myanmar oil and gas pipelines[J].Impact Assessment amp; Project Appraisal，2019，37（1）：71-85.

[12]劉偉，侯立娟，王亞舒.緬甸天然氣市場展望及中國企業(yè)投資合作建議[J].國際石油經(jīng)濟(jì)，2019，27（8）：90-98.

[13]陳湘球.海外企業(yè)在緬項目開發(fā)的創(chuàng)新思考[J].北京石油管理干部學(xué)院學(xué)報，2016（1）：29-32，37.

[14]趙剛，楊英寶，包旭.供應(yīng)鏈網(wǎng)絡(luò)風(fēng)險擴(kuò)散動力學(xué)模型及其應(yīng)用[J].系統(tǒng)工程理論與實踐，2015，35（8）：2014-2024.

[15]辜穗，余曉鐘.基于ISM模型的石油文化動態(tài)演化影響因素系統(tǒng)分析[J].北京石油管理干部學(xué)院學(xué)報，2016，23（3）：58-63.

收稿日期：2023-06-25

作者簡介：

邱湖淼（1986—），男，博士，高級工程師，PMP，研究方向：油氣項目風(fēng)險管理。

魏東（1981—），男，工程師，研究方向：能源管理。

宋振騏（1935—），男，中國科學(xué)院院士，研究方向：能源管理。