王飛，阿斯?jié)h，王衛(wèi)強(qiáng)，柴多

基于LS-SVM的天然氣水合物生成條件預(yù)測(cè)模型建立

王飛，阿斯?jié)h，王衛(wèi)強(qiáng)，柴多

（遼寧石油化工大學(xué), 遼寧 撫順 113001）

在天然氣管線內(nèi)生成的水合物會(huì)嚴(yán)重影響天然氣的開(kāi)采、運(yùn)輸，因而天然氣水合物的預(yù)測(cè)方法和防治措施備受重視。針對(duì)天然氣水合物生成條件，考慮天然氣組分對(duì)水合物生成的影響，為簡(jiǎn)化計(jì)算、提高預(yù)測(cè)精度，引入一種能夠很好解決復(fù)雜物理問(wèn)題的最小二乘支持向量機(jī)（LS-SVM），并且通過(guò)Matlab語(yǔ)言編程，建立了一種包含CH4濃度、CO2濃度、H2S濃度以及水合物生成溫度為輸入，水合物生成壓強(qiáng)為輸出的天然氣水合物生成條件預(yù)測(cè)模型，同時(shí)將實(shí)驗(yàn)數(shù)據(jù)作為最小二乘支持向量機(jī)訓(xùn)練數(shù)據(jù)并進(jìn)行預(yù)測(cè)分析。結(jié)果表明，該預(yù)測(cè)模型不僅擁有較高的預(yù)測(cè)精度，而且方法簡(jiǎn)單、可行，為天然氣水合物生成條件預(yù)測(cè)提供了一種新的解決方法。

天然氣；水合物；生成條件；最小二乘支持向量機(jī)

天然氣在運(yùn)輸?shù)倪^(guò)程中當(dāng)溫度和壓力達(dá)到一定值的條件下水和天然氣組分中的小分子氣體會(huì)形成一種冰雪狀復(fù)合物也就是天然氣水合物，而且水合物的聚合累積會(huì)導(dǎo)致管道堵塞等不良后果，甚至事故的發(fā)生[1]。提高天然氣水合物的預(yù)測(cè)精度，可以為防治水合物提供一定理論支持，從而可以更好地處理由水合物造成的不良后果。

傳統(tǒng)的天然氣水合物生成條件預(yù)測(cè)模型都是以VDW—P模型為基礎(chǔ)，通過(guò)熱力學(xué)和動(dòng)力學(xué)規(guī)律建立起來(lái)的，雖然理論基礎(chǔ)強(qiáng)，但卻存在有計(jì)算復(fù)雜、推廣性不好的缺點(diǎn)[2]。而目前所采用智能算法，例如基于 BP網(wǎng)絡(luò)的預(yù)測(cè)模型[3,4]，雖然有解決了復(fù)雜性并且具有較好的推廣性，但是BP網(wǎng)絡(luò)模型必須是在有大量訓(xùn)練樣本的情況下才有效，在實(shí)際現(xiàn)場(chǎng)數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)都為小樣本的前提下會(huì)對(duì)預(yù)測(cè)精度產(chǎn)生一定的負(fù)面影響。

由Vapnik等人提出了一種在解決小樣本、非線性問(wèn)題方面具有優(yōu)勢(shì)的[5]，并且數(shù)學(xué)理論嚴(yán)密的機(jī)器學(xué)習(xí)算法支持向量機(jī)(SVM)[6,7]。近幾年，SVM憑借著其特有的優(yōu)勢(shì)和極強(qiáng)的泛化能力，已經(jīng)成為了一種新的建模熱點(diǎn)[8]，而且在解決實(shí)際問(wèn)題中得到了成功應(yīng)用[9,10]。在此基礎(chǔ)上由 Sukens等人改進(jìn)并提出了一種新的支持向量機(jī)方法最小二乘支持向量機(jī)（LS-SVM）[11]，其優(yōu)勢(shì)在于，將約束條件更改為等式約束，并且將經(jīng)驗(yàn)風(fēng)險(xiǎn)變?yōu)槎畏?，使用求解線性方程組的辦法得出最優(yōu)化結(jié)果，這樣做不僅降低了復(fù)雜度，而且還加快了收斂精度和求解速度，特別是應(yīng)用在非線性預(yù)測(cè)方面具有很大優(yōu)勢(shì)。

綜上，針對(duì)天然氣水合物這一多組分、多相的復(fù)雜系統(tǒng)，并且存在樣本數(shù)據(jù)較少的情況，選擇最小二乘支持向量機(jī)算法，來(lái)建立天然氣水合物的預(yù)測(cè)模型。

1 算法原理

最小二乘支持向量機(jī)原理如下：

式中：ω—特征空間權(quán)向量；

b∈R—偏置。

將這一優(yōu)化問(wèn)題轉(zhuǎn)化為其對(duì)偶空間中，采用拉格朗日法求解問(wèn)題，定義Lagrange函數(shù)為：

得到最小二乘支持向量機(jī)的估計(jì)函數(shù)為：

從式中能夠看出，最小二乘支持向量機(jī)只需要函數(shù)參數(shù)σ以及正則化參數(shù)c，不需要去選擇不敏感損失函數(shù)的ε值，從而簡(jiǎn)化了計(jì)算。

為了能夠提高LS-SVM算法的效率及效果，需要對(duì)訓(xùn)練數(shù)據(jù)和預(yù)測(cè)數(shù)據(jù)的輸入量和輸出量進(jìn)行歸一化處理：

式中：xmax—數(shù)據(jù)最大值；

xmin—數(shù)據(jù)最小值。

2 天然氣水合物預(yù)測(cè)模型建立

為了更好的驗(yàn)證上述算法的可靠性，所選仿真平臺(tái)為：2.5 GHz處理器，4 G內(nèi)存，Windows 7操作系統(tǒng)PC上的 Matlab2012a軟件平臺(tái)，并如圖1所示流程建立天然氣水合物生成條件預(yù)測(cè)模型。

圖1 天然氣水合物預(yù)測(cè)模型建模流程圖Fig.1 Natural gas hydrate forecast model

天然氣水合物是天然氣在管線內(nèi)運(yùn)輸時(shí)，當(dāng)溫度和壓強(qiáng)達(dá)到一定值的時(shí)候，天然氣組份內(nèi)的小分子氣體和水形成的，考慮到天然氣組分對(duì)水合物生成造成的影響，確定CH4濃度、CO2濃度、H2S濃度、水合物生成時(shí)溫度T作為L(zhǎng)S-SVM輸入變量，水合物生成時(shí)壓強(qiáng)P作為L(zhǎng)S-SVM輸出變量。

并且在世界范圍內(nèi)酸性天然氣普遍使用，所以引用相關(guān)文獻(xiàn)中所提供的59組（CH4+ CO2+ H2S）三元酸性天然氣在純水中的實(shí)驗(yàn)數(shù)據(jù)（xi,yi）（i=1，2，……59）作為樣本數(shù)據(jù)[12]，并將訓(xùn)練樣本數(shù)據(jù)帶入公式（6）進(jìn)行歸一化處理，將經(jīng)過(guò)處理之后的前50組實(shí)驗(yàn)數(shù)據(jù)作為模型的訓(xùn)練集，后9組數(shù)據(jù)作為模型的測(cè)試集，利用訓(xùn)練集建立天然氣水合物的預(yù)測(cè)模型，并調(diào)用Matlab中的LS-SVM工具箱，對(duì)數(shù)據(jù)進(jìn)交叉驗(yàn)證，從而確定核函數(shù)參數(shù)σ=20正則化參數(shù)c=12。

最后利用經(jīng)過(guò)歸一化處理的后9組實(shí)驗(yàn)數(shù)據(jù)進(jìn)行預(yù)測(cè)分析，從而驗(yàn)證出預(yù)測(cè)模型的預(yù)測(cè)精度以及推廣性質(zhì)。

3 預(yù)測(cè)結(jié)果分析

基于文獻(xiàn)數(shù)據(jù)[12]，采用Matlab語(yǔ)言編寫(xiě)最小二乘向量機(jī)程序，將經(jīng)過(guò)歸一化處理的后9組實(shí)驗(yàn)作為預(yù)測(cè)集，進(jìn)行預(yù)測(cè)分析，其中T為實(shí)驗(yàn)溫度，P為實(shí)驗(yàn)溫度下天然氣水合物生成時(shí)的壓強(qiáng)，Pc為預(yù)測(cè)模型預(yù)測(cè)出的天然氣水合物生成壓強(qiáng)，ΔP為預(yù)測(cè)值與實(shí)驗(yàn)值之間的誤差，結(jié)果見(jiàn)下圖2，詳細(xì)數(shù)據(jù)見(jiàn)下表1。

圖2 預(yù)測(cè)值與仿真實(shí)驗(yàn)對(duì)比圖Fig.2 The predicted values and the simulation contrast figure

圖2 中，星號(hào)代表在實(shí)驗(yàn)條件下天然氣水合物生成時(shí)的壓力值，圓圈代表預(yù)測(cè)的天然氣水合物生成壓力值，表1為預(yù)測(cè)數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)的詳細(xì)信息。分析可知，其平均誤差在 0.11%，可見(jiàn)預(yù)測(cè)模型的預(yù)測(cè)精度較好，可以應(yīng)用于天然氣水合物的生成預(yù)測(cè)，并且有著較好的推廣性。

表1 實(shí)驗(yàn)值與仿真預(yù)測(cè)實(shí)驗(yàn)誤差Table 1 The experimental data and simulation experimental error

4 結(jié) 論

（1）在分析現(xiàn)有的天然氣水合物生成條件預(yù)測(cè)模型的優(yōu)缺點(diǎn)基礎(chǔ)上，考慮天然氣組分對(duì)天然氣水合物生成的影響，引入了一種在解決小樣本數(shù)據(jù)上具有優(yōu)勢(shì)的最小二乘支持向量機(jī)（LS-SVM），并通過(guò)Matlab語(yǔ)言編程，建立了天然氣水合物生成條件預(yù)測(cè)模型。

（2）通過(guò)參閱一系列相關(guān)文獻(xiàn)，選擇了酸性天然氣在純水中的實(shí)驗(yàn)數(shù)據(jù)作為本模型的訓(xùn)練集和測(cè)試集，由于所選的輸入變量為實(shí)驗(yàn)數(shù)據(jù)，存在有一定的局限性，在實(shí)際應(yīng)用中會(huì)對(duì)誤差有所影響，以后應(yīng)該樣本數(shù)量更加豐富，并且優(yōu)化模型的數(shù)據(jù)處理，從而強(qiáng)本預(yù)測(cè)模型的推廣性。

（3）利用本模型對(duì)水合物生成條件進(jìn)行預(yù)測(cè)分析，從中可以看出實(shí)驗(yàn)數(shù)據(jù)與預(yù)結(jié)果具有很好的一致性，驗(yàn)證了本模型的有效、可靠。

（4）實(shí)際上影響天然氣水合物生成的因素還有很多，而且內(nèi)在關(guān)系復(fù)雜，運(yùn)用最小二乘支持向量機(jī)算法對(duì)天然氣水合物生成條件進(jìn)行預(yù)測(cè)，作為一種新的嘗試，為天然氣水合物預(yù)測(cè)提供了一種新的途徑，同時(shí)為后續(xù)研究奠定了基礎(chǔ)。

［1］王新，裘俊紅.籠型水合物研究和應(yīng)用[J].化工生產(chǎn)與技術(shù)，2004,11: 19-25.

［2］喻西崇，趙金州,等.天然氣水合物生成條件預(yù)測(cè)模型的比較[J].油氣儲(chǔ)運(yùn),2002,21(1):20-24.

［3］顧建新，韓小輝,等.用人工神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)天然氣水合物的生成[J].天然氣工業(yè),2003,23(5):102-104.

［4］彭遠(yuǎn)近,劉建儀,等.天然氣水合物生成條件的神經(jīng)網(wǎng)絡(luò)方法預(yù)測(cè)[J].天然氣工業(yè),2006,26(7):85-87.

［5］Cao L J, Tay F E H. Support vector machine with adaptive parameters in financial time forecasting[J]. IEEE Trans. on Nature Network, 2003, 14(6): 1506-1518.

［6］Vapnik V. The Nature of Statistical Learning Theory[M]. New York: Springer-Verlag, 1995.

［7］Vapnik V N. An Overview of statistical learning theory[J]. IEEE Trans. on Nature Network, 1999, 10(5): 988-999.

［8］Xiong Weili(熊偉麗), Yao Yue(姚樂(lè)), Xu Baoguo(徐保國(guó)). Chaos Least Squares Support Vector Machine and Its Application on Fermentation Process Modeling[J]. CIESC Journal(化工學(xué)報(bào)), 2013, 64(12): 4587-4590.

［9］Qian Xiaoshan(錢(qián)曉山), Yang Chunhua(陽(yáng)春華), Xu Lisha(徐麗莎). Soft sensor of sodium aluminate solution concentration based on improved differential evolution algorithm and LSSVM[J]. CIESC Journal(化工學(xué)報(bào)), 2013, 64(5): 1074-1079.

［10］Li Hao(李浩), Yang Min(楊敏). Soft sensing method based on principal curves and its application in distillation columns[J]. CIESC Journal(化工學(xué)報(bào)), 2012, 63(8): 2492-2499.

［11］Suykens J A K, Gestel T V, Brabanter J D. Least Squares Support Vector Machines[M].Singapore: World Scientific Publishers,2002.

［12］黃強(qiáng)，孫長(zhǎng)宇，陳光進(jìn),等.含( CH4+ CO2+ H2S) 酸性天然氣水合物形成條件實(shí)驗(yàn)與計(jì)算[J] . 化工學(xué)報(bào), 2005, 56 (7) : 1159- 1163.

Establishment of the Prediction Model of the Forming Conditions of Natural Gas Hydrate Based on Least Square Support Vector Machine

WANG Fei，A Si-han，WANG Wei-qiang，CHAI Duo
（Liaoning Shihua University, Liaoning Fushun 113001，China）

The hydrate generated in the natural gas pipeline will seriously affect the natural gas production, transportation, and the prediction methods and control measures of natural gas hydrate (NGH) have attracted great attention. For natural gas hydrate formation conditions, considering the influence of the components on the gas hydrate formation, in order to simplify the calculation and improve the prediction accuracy, through introducing least squares support vector machines (LS - SVM) and using Matlab language programming, a prediction model of natural gas hydrate formation conditions was established with the concentrations of CH4, H2S , CO2and hydrate formation temperature as the input, the pressure of hydrate formation as the output. At the same time, the experimental data were used as Least Squares Support Vector Machine training data to carry out the forecast analysis. The results show that the forecasting model has higher prediction accuracy, and the method is simple and feasible, can provide a new solution for prediction of natural gas hydrate formation conditions.

Natural gas; Hydrate; Formation conditions; Least squares support vector machine

TE 8

： A

： 1671-0460（2015）04-0789-03

2015-01-13

王飛（1988-），男，吉林長(zhǎng)春人，碩士研究生，2015年畢業(yè)于遼寧石油化工大學(xué)油氣儲(chǔ)運(yùn)工程專(zhuān)業(yè)，研究方向：管道輸運(yùn)技術(shù)研究。E-mail：tywanghk@126.com。