范 婕，許欣怡，周詩崠，周年勇

（常州大學 石油工程學院 江蘇省油氣儲運技術(shù)重點實驗室，江蘇 常州 213164）

天然氣水合物又稱籠型水合物，是一種類冰狀的結(jié)晶物質(zhì)，由天然氣（主要成分為甲烷）和一定量的水在低溫高壓的條件下形成[1]。作為一種清潔能源，天然氣水合物的應(yīng)用前景十分廣闊[2]。然而在天然氣集輸管道中，天然氣水合物的生成會影響輸送管線的安全運行。大量生成的天然氣水合物在流動過程中會集聚并堵塞管道，導致管內(nèi)壓差增大，嚴重時會損壞管道[3-4]。為確保天然氣在管道中的流動安全，提高天然氣水合物生成條件的預(yù)測精度是十分有必要的[5]。

在天然氣水合物形成條件的預(yù)測方面，國內(nèi)外學者取得了相對豐富的研究成果，形成了以vdW-P模型和Chen-Guo模型[6-7]為代表的熱力學計算模型。但在實際過程中，影響天然氣水合物生成的因素較為復雜，需要考慮的變量較多，計算復雜繁瑣，這限制了純理論模型的應(yīng)用。為進一步拓展水合物生成條件的預(yù)測方法，REBAI等[8]提出了一種預(yù)測混合組分體系水合物生成壓力的方法，使用人工神經(jīng)網(wǎng)絡(luò)（ANN），對基于vdW-P模型的熱力學模型近似地補充了壓力修正項，結(jié)果顯示，壓力的總體相對誤差從23.75%顯著降低到3.15%。SOROUSH等[9]建立了一種前饋式ANN用于預(yù)測糖類和酸性氣體體系的水合物形成，該模型與其他關(guān)聯(lián)式熱力學模型相比，尤其是對混酸氣體進行生成預(yù)測時，預(yù)測精度很高，總均方誤差為0.349。MEHRIZADEN[10]采用ANN和自適應(yīng)神經(jīng)模糊干擾系統(tǒng)（ANFIS）來估算不同天然氣系統(tǒng)的水合物生成壓力，將預(yù)測結(jié)果與實驗方程進行比較，結(jié)果表明，ANFIS模型在所有情況下都優(yōu)于ANN模型，且ANFIS模型比經(jīng)驗公式的預(yù)測精度更高。鄭秋海等[11]通過深度神經(jīng)網(wǎng)絡(luò)提取生產(chǎn)數(shù)據(jù)的網(wǎng)絡(luò)特征，將提取的特征融合生產(chǎn)數(shù)據(jù)來增強數(shù)據(jù)區(qū)分度，使用非線性支持向量機（SVM）對融合數(shù)據(jù)進行水合物生成預(yù)測，實驗表明深度神經(jīng)網(wǎng)絡(luò)結(jié)合支持向量機（FDNN-SVM）模型滿足海底管線生產(chǎn)運行需求。彭炎等[12]采用支持向量回歸方法進行天然氣水合物成藏預(yù)測研究，結(jié)果表明，支持向量回歸方法在凍土區(qū)天然氣水合物生成預(yù)測方面的應(yīng)用是有效的。徐小虎等[13]建立了基于粒子群算法（PSO）的最小二乘支持向量機（LSSVM）模型和誤差反向傳播（BP）神經(jīng)網(wǎng)絡(luò)，對含抑制劑體系下的天然氣水合物生成進行了預(yù)測，結(jié)果表明，前者更為精準。可見，支持向量機能夠結(jié)合不同的優(yōu)化算法，對水合物生成預(yù)測領(lǐng)域進行研究。

本文針對天然氣水合物樣本數(shù)少這一問題，將粒子群算法結(jié)合SVM，選用預(yù)測效果好的核函數(shù)構(gòu)建預(yù)測模型，預(yù)測純水體系中，純組分和多組分天然氣水合物的生成條件。PSO-SVM模型的輸入?yún)?shù)為天然氣組分及溫度，輸出參數(shù)為壓力。

1 模擬部分

1.1 PSO-SVM模型的算法原理

本文選用SVM和PSO算法建立天然氣水合物生成預(yù)測模型，并借助MATLAB語言編程實現(xiàn)。實驗中能夠得到的天然氣水合物生成數(shù)據(jù)并不多，支持小樣本數(shù)據(jù)下的訓練預(yù)測是SVM的特點之一，因此選用SVM構(gòu)建天然氣水合物生成預(yù)測模型。天然氣水合物生成數(shù)據(jù)所繪制的溫度壓力圖是近似平滑的曲線，SVM在曲線擬合方面也較為精準。PSO算法不需要額定的參數(shù)，靈活性強，與SVM相結(jié)合能解決SVM初始準確度不高和自身學習能力較弱的缺點。

1.1.1 支持向量機

SVM定義最優(yōu)線性超平面，并把尋找最優(yōu)線性超平面的算法歸納為求解一個凸規(guī)劃問題[14]?；贛ercer核展開定理，通過非線性映射Φ，把樣本空間映射到一個高維乃至無窮維的特征空間（Hilbert空間），使在特征空間中可以應(yīng)用線性學習機的方法，解決樣本空間中的高度非線性分類和預(yù)測等問題[15]。簡而言之，就是先將問題升維，然后線性化。SVM是用于分類和回歸的算法。天然氣水合物生成條件預(yù)測是一個低維、非線性問題，本文將對天然氣水合物生成條件構(gòu)建SVM回歸預(yù)測模型，使其變成能夠被解決的高維、線性回歸問題。

對于SVM的基本原理，給定訓練樣本D= {(xi,yi),i= 1, 2, …,N}，xi為自變量，yi為因變量，N為樣本量，回歸模型如下：

式中，f(x)為模型的輸出變量；ω為特征空間權(quán)向量；x為輸入變量；b為偏置向量。

SVM回歸中所使用的結(jié)構(gòu)風險函數(shù)表達式如下：

式中，C為懲罰因子；m為樣本量；Lε為ε不敏感損失函數(shù)；ε為損失函數(shù)的損失因子。

SVM回歸函數(shù)如下：

式中，αi為第i個樣本的Lagrange乘子，K(xi,x)為內(nèi)積，即核函數(shù)。

對于內(nèi)積問題，一般會選用性能較好、應(yīng)用較多的RBF核函數(shù)，其表達式如下：

式中，g為核函數(shù)的寬度因子。

1.1.2 粒子群算法

PSO是一種進化計算技術(shù)[16]，具有容易實現(xiàn)和不需要許多參數(shù)調(diào)節(jié)的特點。這使得PSO能夠被廣泛地應(yīng)用到函數(shù)優(yōu)化、神經(jīng)網(wǎng)絡(luò)及其參數(shù)優(yōu)化和模糊系統(tǒng)控制等應(yīng)用領(lǐng)域[17-18]。PSO通過設(shè)計一種無質(zhì)量的例子來模擬鳥群中的鳥，粒子僅具有速度和位置兩個屬性，用速度體現(xiàn)粒子移動的快慢，用位置體現(xiàn)粒子移動的方向。

粒子通過下式來更新自身的速度和位置：

式中，i= 1, 2, …,N，N為此群中粒子的總數(shù)；vi為粒子的速度；rand()為介于(0, 1)之間的隨機數(shù)；pbesti為每個粒子的個體最優(yōu)值；gbesti為全局粒子的最優(yōu)值；xi為粒子的當前位置；c1和c2為學習因子，通常取2；vi的最大值為Vmax（大于0），如果vi大于Vmax，則vi= Vmax。

1.1.3 PSO-SVM預(yù)測模型的實現(xiàn)方法

管線輸送天然氣的過程中，天然氣中各氣體組分含量（物質(zhì)的量分數(shù)）、溫度和壓力等都是天然氣水合物生成的影響因素。將天然氣的組分含量、溫度作為PSO-SVM模型的變量，將壓力作為因變量。為了減少迭代時間，保證算法的高效性，需要對SVM中的懲罰因子C和核函數(shù)參數(shù)g進行優(yōu)化。PSO中的個體最優(yōu)對應(yīng)SVM中的懲罰因子C，群體最優(yōu)對應(yīng)核函數(shù)參數(shù)g。

本文中PSO-SVM模型的訓練過程如下：（1）導入原始數(shù)據(jù)，按比例隨機劃分訓練集和測試集，同時對數(shù)據(jù)進行歸一化處理；（2）對PSO進行初始化參數(shù)設(shè)置，將訓練集代入到PSO中，結(jié)合SVM迭代運算，更新權(quán)值，得到兩個最優(yōu)參數(shù)懲罰因子C和核函數(shù)參數(shù)g；（3）將經(jīng)過PSO計算后的最優(yōu)參數(shù)代入SVM模型中，用測試集進行運算，對得到的數(shù)據(jù)反歸一化；（4）計算實際值與預(yù)測值的誤差。

1.2 模型的準備工作

1.2.1 實驗數(shù)據(jù)準備及處理

通過調(diào)研現(xiàn)有文獻[19-22]數(shù)據(jù)，得到純水體系下水合物生成數(shù)據(jù)共132組。其中，純組分天然氣水合物生成數(shù)據(jù)69組，溫度范圍273.49～298.20 K，壓力范圍2.71～43.80 MPa；多組分天然氣水合物生成數(shù)據(jù) 63組，溫度范圍 273.59～298.00 K，壓力范圍0.82～24.85 MPa，氣體組分如表1。本研究以該數(shù)據(jù)，建立基于PSO-SVM模型的天然氣水合物生成預(yù)測模型，并以8:2的比例隨機劃分訓練集和測試集。

對數(shù)據(jù)進行訓練之前，需要先對數(shù)據(jù)進行歸一化處理，使得數(shù)據(jù)范圍穩(wěn)定，避免存在極端的最大最小值，同時避免數(shù)據(jù)中的異常值和極端值帶來的影響。歸一化的效果是將原數(shù)據(jù)規(guī)整到一個閉區(qū)間內(nèi)，通過Mapminmax函數(shù)來實現(xiàn)。該程序選用的區(qū)間是[1, 2]，即yi∈[0, 1](i= 1, 2, 3, …,N)，這種歸一化方式稱為[1, 2]區(qū)間歸一化。所用的歸一化映射如下式：

式中，x,y∈R；xmin= min(x)；xmax= max(x)。

1.2.2 SVM模型核函數(shù)選取

就SVM中無法處理的非線性可分問題而言，數(shù)據(jù)需要從低維向高維映射，核函數(shù)能夠?qū)⑻卣饔成涞礁呔S的空間，而計算過程在低維中運行，實際效果表現(xiàn)在高維上，從而解決了維度爆炸的問題。常用的核函數(shù)有線性核函數(shù)、多項式核函數(shù)、Sigmoid核函數(shù)和徑向基核函數(shù)[23]。本文在實際運算過程中，以PSO-SVM模型為基礎(chǔ)進行了試算，計算結(jié)果見2。由表2可知，徑向基核函數(shù)的結(jié)果最理想，訓練集平均平方誤差（MSE）和測試集MSE最小，因此選擇徑向基核函數(shù)。

表2 不同核函數(shù)下PSO-SVM模型的參數(shù)及誤差Table 2 Parameters and errors of PSO-SVM model with different kernel functions

1.2.3 運算參數(shù)

PSO參數(shù)設(shè)置對算法本身的運行有很大影響[24]。為進一步研究PSO-SVM模型在天然氣水合物生成預(yù)測領(lǐng)域中的精準度，對PSO初始化參數(shù)進行了優(yōu)化設(shè)計，見表3。PSO-SVM模型運算得到的最優(yōu)懲罰因子C為97.5331，最優(yōu)核函數(shù)參數(shù)g為0.6439。

表3 PSO初始參數(shù)Table 3 Initial parameters of PSO

2 模擬結(jié)果及分析

2.1 不適應(yīng)度

PSO-SVM模型模擬的適應(yīng)度較高，所以在程序編寫過程中，對適應(yīng)度的倒數(shù)進行曲線描述，如圖1所示。由圖1可知，隨著PSO迭代次數(shù)增加，適應(yīng)度倒數(shù)的值越小，該值也指均方差均值，誤差越小越好。說明PSO-SVM模型適用于天然氣水合物的生成預(yù)測。

圖1 適應(yīng)度曲線Fig.1 Fitness curve

2.2 純組分天然氣水合物生成預(yù)測

通過PSO-SVM模型，對收集到的數(shù)據(jù)組進行訓練預(yù)測。純組分天然氣水合物生成預(yù)測結(jié)果如圖2所示。圖2中方點代表天然氣水合物生成的真實壓力，圓點代表PSO-SVM模型計算得到的預(yù)測壓力。從圖2可知，兩條曲線的數(shù)據(jù)點很接近，但在高于295 K后，預(yù)測壓力誤差較大。說明在 273.49～295.00 K之間，PSO-SVM模型預(yù)測的可靠性很高，模型相對誤差的最大值為8.768%，最小值為0.015%。

圖2 純組分天然氣水合物生成預(yù)測Fig.2 Formation prediction of pure component natural gas hydrate

2.3 多組分天然氣水合物生成預(yù)測

4個不同體系多組分天然氣水合物生成預(yù)測結(jié)果如圖3所示。由圖3可知，體系1中，兩條曲線的數(shù)據(jù)點基本吻合，個別預(yù)測值誤差較大；體系2中，數(shù)據(jù)點在273.80～291.00 K之間的誤差較小，但溫度高于291 K后，數(shù)據(jù)的吻合度不高；體系3和體系4中，兩條曲線的數(shù)據(jù)點大體吻合?？偟目磥恚琍SO-SVM模型對多組分天然氣水合物生成預(yù)測的準確度也很高。

圖3 多組分天然氣水合物生成預(yù)測Fig.3 Formation prediction of multi-component gas hydrate

2.4 PSO-SVM模型誤差驗證

為進一步檢驗PSO-SVM模型預(yù)測的可行性和準確性，采用MSE、平方相關(guān)系數(shù)（R2）和平均絕對比例誤差（MAPE）等指標進行評價。評價指標定義如下：

式中，n為樣本數(shù)；Yi為實際值；為預(yù)測值。

MSE為所有樣本的樣本誤差絕對值的均值，該值越接近0，模型越準確；R2為因變量的方差能被自變量解釋的程度，該值越接近1，模型效果越好；MAPE為所有樣本誤差的絕對值占實際值的比值，該值越接近0，模型越準確。PSO-SVM模型對測試集預(yù)測的評價結(jié)果見表4。

表4 評價結(jié)果Table 4 Evaluation results

由表4可知，對于純組分和多組分天然氣水合物，模型的MSE和R2都很理想，并且MAPE也比較低。

3 結(jié)論

建立了預(yù)測純組分和多組分天然氣水合物生成條件的PSO-SVM模型，并對模型進行了分析與評價，得到如下結(jié)論。

（1）SVM的核函數(shù)是影響PSO-SVM模型預(yù)測結(jié)果的主要因素，其中徑向基核函數(shù)的預(yù)測效果最優(yōu)。預(yù)測天然氣水合物生成條件，模型最優(yōu)懲罰因子C為97.5331，最優(yōu)核函數(shù)參數(shù)g為0.6439，訓練集和測試集的MAPE分別為2.74%和2.84%。

（2）PSO-SVM模型中，純組分和多組分天然氣水合物的適用溫度分別為273.49～295.00 K和273.59～298.00 K。

（3）PSO-SVM模型對純水體系下，純組分和多組分天然氣水合物生成條件具有較好的預(yù)測精度，預(yù)測純組分天然氣水合物時，MSE、R2和MAPE分別為0.0003963、0.9996和2.84%；預(yù)測多組分天然氣水合物時，MSE、R2和MAPE分別為 0.0006870、0.9983和2.74%。