楊 威，薛 釗，蔡文慧，王 婷，田東海，梁昌晶

1.華北石油管理局有限公司蘇里格勘探開發(fā)分公司，內(nèi)蒙古鄂爾多斯150626

2.中國(guó)石油華北油田公司友信勘探開發(fā)服務(wù)有限公司，河北任丘062552

3.中國(guó)石油華北油田公司第二采油廠，河北霸州065700

4.中國(guó)石油華北油田公司二連分公司，內(nèi)蒙古錫林浩特026000

目前，水合物生成條件的預(yù)測(cè)方法主要有經(jīng)驗(yàn)圖解法、關(guān)聯(lián)公式法和熱力學(xué)模型法[1-3]，其中經(jīng)驗(yàn)圖解法和關(guān)聯(lián)公式法在油氣田現(xiàn)場(chǎng)應(yīng)用較多，但預(yù)測(cè)精度及效果不好，熱力學(xué)模型法是基于室內(nèi)試驗(yàn)的結(jié)果，通過(guò)引入狀態(tài)方程以提高預(yù)測(cè)結(jié)果精度，但涉及的參數(shù)較多，專業(yè)性較強(qiáng)，不能很好地推廣使用。

隨著機(jī)器學(xué)習(xí)技術(shù)的不斷發(fā)展，越來(lái)越多的學(xué)者采用人工智能算法對(duì)水合物的生成條件進(jìn)行預(yù)測(cè)[4]。卞小強(qiáng)等[5]通過(guò)引入CO2和H2S 的貢獻(xiàn)因子，采用SVM 對(duì)含酸性氣體的水合物生成條件進(jìn)行了預(yù)測(cè)，其平均相對(duì)偏差為5.7%；Mesbah 等[6]為了解決SVM 算法產(chǎn)生局部最優(yōu)問(wèn)題，引入了最小二乘向量機(jī)，可有效預(yù)測(cè)高酸天然氣水合物的生成；唐永紅等[7]將小波分析和神經(jīng)網(wǎng)絡(luò)相結(jié)合，通過(guò)求解權(quán)值向量、平滑因子、伸縮因子等變量，對(duì)I 型水合物的生成溫度進(jìn)行了預(yù)測(cè)；馬貴陽(yáng)等[8]將遺傳算法與SVM 相結(jié)合，通過(guò)數(shù)值計(jì)算擬合了相平衡曲線，其精度與熱力學(xué)模型的精度接近。綜上所述，以上研究均基于大量的試驗(yàn)數(shù)據(jù)，本身容易陷入局部最優(yōu)，出現(xiàn)超參數(shù)選取困難的現(xiàn)象，且未考慮抑制劑、水中鹽類的影響，對(duì)水合物結(jié)構(gòu)進(jìn)行智能算法分類也未見(jiàn)報(bào)道。在此，采用具有分類和回歸功能的隨機(jī)森林算法[9-10]，以不同的氣體組分為輸入變量，對(duì)水合物結(jié)構(gòu)進(jìn)行分類；以氣體組分、壓力、抑制劑、鹽類等為輸入變量，對(duì)水合物的生成溫度進(jìn)行預(yù)測(cè)，以期為水合物的相關(guān)研究提供理論依據(jù)和實(shí)際參考。

1 水合物相平衡影響因素

影響水合物相平衡的因素可分為氣體組分和外部因素兩方面，外部因素包括抑制劑、鹽類、溫度、壓力等[11]。

（1） 氣體組分。首先，水合物結(jié)構(gòu)的形成主要與氣體組分相關(guān)，CH4、C2H6等小分子氣體和非烴類氣體可形成I 型水合物；C3H8、i-C4H10和非烴類氣體可形成II 型水合物；H 型水合物的形成需要2，2-二甲基丁烷、甲基環(huán)己烷、甲基環(huán)戊烷等大分子配合CH4、N2等小分子[12]，在現(xiàn)場(chǎng)未檢測(cè)到有關(guān)組分，在此不予討論。其次，氣體組分的差異化可改變相平衡曲線，在純CH4中分別加入C2H6和C3H8可以促進(jìn)水合物的生成，且C3H8與十六面體大空腔的結(jié)合作用更強(qiáng)；加入H2S 和CO2也會(huì)使相平衡曲線右移，由于H2S 的溶解度更高，對(duì)水合物的促進(jìn)作用更強(qiáng)。

（2） 水合物抑制劑。常用的水合物抑制劑有甲醇、乙醇、乙二醇、二甘醇、三甘醇[13]，不同抑制劑的抑制效果不同。隨著抑制劑摩爾分?jǐn)?shù)的增加，甲醇的抑制效果基本不變，乙醇和二甘醇的抑制效果變小，乙二醇的抑制效果變大，三甘醇的抑制效果先變小后變大，主要與氣- 水兩相體系中非電荷基團(tuán)的數(shù)量和性能有關(guān)。

（3） 水中鹽類。對(duì)于氣田采出水，水中鹽類的成分和含量對(duì)水合物相平衡影響較大[14]。由于水中溶解的離子會(huì)形成強(qiáng)弱不同的電場(chǎng)，而水合物形成籠型結(jié)構(gòu)需要額外的能量來(lái)破壞這種電場(chǎng)，根據(jù)分子間氫鍵的破壞程度，Cl-對(duì)水合物的抑制作用比S要強(qiáng)[15]。

（4） 溫度、壓力。高壓、低溫環(huán)境下容易生成水合物，這是水合物生成的最重要熱力學(xué)因素[16]。低壓下相平衡曲線較敏感，高壓下敏感程度降低；低溫下相平衡曲線不敏感，高溫下敏感程度上升很快。

2 數(shù)據(jù)來(lái)源及研究方法

2.1 數(shù)據(jù)來(lái)源

分類問(wèn)題方面，以不同的氣體組分為輸入變量，其中C1+C2的摩爾分?jǐn)?shù)為第一變量，C3+C4+N2的摩爾分?jǐn)?shù)為第二變量，C5+摩爾分?jǐn)?shù)為第三變量，CO2摩爾分?jǐn)?shù)為第四變量，H2S 摩爾分?jǐn)?shù)為第五變量，以水合物的結(jié)構(gòu)類型為輸出值，對(duì)結(jié)構(gòu)類型進(jìn)行數(shù)字化標(biāo)簽轉(zhuǎn)化，I 型水合物定義為1，II 型水合物定義為2?；貧w問(wèn)題方面，以氣體組分、壓力、抑制劑含量、鹽類含量等為輸入變量，其中氣體組分變量輸入與回歸模型一致，抑制劑中以甲醇和乙二醇應(yīng)用最為廣泛，且一般不同時(shí)使用，以甲醇+乙二醇摩爾分?jǐn)?shù)為第六變量，對(duì)水合物形成影響最大的孔隙水中主要含有NaCl 和MgCl2，以NaCl+MgCl2的摩爾分?jǐn)?shù)為第七變量，以井口壓力為第八變量，以Du-Guo 公式計(jì)算的水合物形成溫度為輸出變量。

采用華北油田2018—2019 年氣井的水合物數(shù)據(jù)，包括氣質(zhì)組分、水質(zhì)組分、抑制劑注入量、壓力、溫度等，這些數(shù)據(jù)充分考慮了不同影響因素對(duì)水合物生成的貢獻(xiàn)程度，按照現(xiàn)場(chǎng)工況，CH4摩爾分?jǐn)?shù)為80%～96%，C2H6摩爾分?jǐn)?shù)為0%～5%，C3H8摩爾分?jǐn)?shù)為0%～5%，H2S 摩爾分?jǐn)?shù)為0%～5%，CO2摩爾分?jǐn)?shù)為0%～10%，N2摩爾分?jǐn)?shù)為0%～6%，抑制劑在水相中的摩爾分?jǐn)?shù)為0% ～30% ，NaCl+MgCl2在水相中的摩爾分?jǐn)?shù)為3%～8%，井口經(jīng)節(jié)流后壓力不超過(guò)5 MPa。取100 組數(shù)據(jù)，按照4∶1 的比例，取其中80 組作為訓(xùn)練集用于建立模型，取其中20 組作為測(cè)試集用于評(píng)價(jià)模型，部分?jǐn)?shù)據(jù)見(jiàn)表1。

表1 水合物生成條件部分?jǐn)?shù)據(jù)

2.2 隨機(jī)森林算法

隨機(jī)森林算法RF（Random Forest） 由Leo Breiman 在2001 年提出，根據(jù)bootstrap 重抽樣方法，從原始數(shù)據(jù)樣本中隨機(jī)抽取M 個(gè)觀測(cè)值，同時(shí)再隨機(jī)抽取k 個(gè)自變量作為分類樹的節(jié)點(diǎn)，產(chǎn)生成百上千個(gè)分類樹，是決策樹的隨機(jī)集成[17]。對(duì)于分類問(wèn)題，可根據(jù)不同分類樹的投票結(jié)果，按照特征判定標(biāo)準(zhǔn)決定分類結(jié)果，模型為RFC（Random Forest Classification）；對(duì)于回歸問(wèn)題，可根據(jù)預(yù)測(cè)均值，按照待回歸屬性決定回歸結(jié)果，模型為RFR（Random Forest Regression）。隨機(jī)森林算法中最為重要的兩個(gè)超參數(shù)為樹的數(shù)量（Nt） 和候選特征子集的數(shù)量（Mt）[18]，同時(shí)隨機(jī)森林算法本身對(duì)節(jié)點(diǎn)最小樣本數(shù)（Nodesize） 不敏感，在此不予考慮。采用運(yùn)算速度最快的網(wǎng)格搜索法GS（Grid Search），通過(guò)對(duì)不同固定范圍內(nèi)的超參數(shù)進(jìn)行尋優(yōu)，得到最優(yōu)解。

2.3 評(píng)價(jià)指標(biāo)

為了防止模型出現(xiàn)高方差或高偏差問(wèn)題，需要對(duì)模型的適應(yīng)性進(jìn)行評(píng)估，在此采用K 折交叉驗(yàn)證的方式，K 取5。采用二維混淆矩陣對(duì)分類模型的預(yù)測(cè)結(jié)果進(jìn)行評(píng)價(jià)，評(píng)價(jià)指標(biāo)為召回率Recall、精確率Precision 以及召回率和精確率的調(diào)和均值Fβ，為了加大Recall 的相對(duì)重要程度，β 值取2，公式如下：

采用相對(duì)誤差百分比、平均絕對(duì)百分比誤差MAPE、均方根誤差RMSE 對(duì)回歸模型的預(yù)測(cè)結(jié)果進(jìn)行評(píng)價(jià)，公式如下：

式中：n 為測(cè)試集樣本的個(gè)數(shù)，n=16；y 為實(shí)際值；y'為預(yù)測(cè)值，即經(jīng)隨機(jī)森林算法計(jì)算的預(yù)測(cè)值。

3 結(jié)果與討論

3.1 水合物結(jié)構(gòu)類型預(yù)測(cè)

（1） 超參數(shù)Mt預(yù)選。采用GS 對(duì)超參數(shù)Mtry進(jìn)行預(yù)選，采用殘差均方值（mean of squared residuals） 和擬合優(yōu)度（var explained） 作為判定依據(jù)，見(jiàn)表2。當(dāng)Mt=3 時(shí)，殘差均方值最小，擬合優(yōu)度最大，模型在計(jì)算精度、運(yùn)算時(shí)間、擬合程度上最優(yōu)。

表2 不同Mt 的殘差均方值和擬合優(yōu)度

當(dāng)Mt=3 時(shí)，對(duì)不同決策樹的數(shù)量Nt的殘差均方值變化進(jìn)行驗(yàn)算，見(jiàn)圖1。當(dāng)Nt＞100 時(shí)，模型誤差較小，但Nt在100～200 之間時(shí)，仍有小幅波動(dòng)。Nt數(shù)量過(guò)小，訓(xùn)練不穩(wěn)定，容易過(guò)擬合，過(guò)大則影響運(yùn)算速度，綜合各方因素，Nt取300。

圖1 Nt 數(shù)量與殘差均方值的關(guān)系

（2） 模型的適應(yīng)性對(duì)比。為了進(jìn)一步評(píng)價(jià)模型的適應(yīng)性，對(duì)比不同的超參數(shù)尋優(yōu)方法和分類模型，選取經(jīng)網(wǎng)格搜索優(yōu)化后的支持向量機(jī)分類模型（GS-SVC）、經(jīng)遺傳算法優(yōu)化后的支持向量機(jī)分類模型（GA-SVC）、神經(jīng)網(wǎng)絡(luò)分類模型（BP） 與經(jīng)網(wǎng)格搜索優(yōu)化后的隨機(jī)森林分類模型（GS-RFC）進(jìn)行對(duì)比，見(jiàn)表3。

表3 不同超參數(shù)尋優(yōu)方法和分類模型的交叉驗(yàn)證結(jié)果

GA-SVC 和BP 算法的訓(xùn)練集和測(cè)試集的準(zhǔn)確率均較低，訓(xùn)練集的準(zhǔn)確率不超過(guò)88%，測(cè)試集的準(zhǔn)確率不超過(guò)83%，這是由于SVC 求解的是二次規(guī)劃問(wèn)題，計(jì)算量較大，雖然經(jīng)過(guò)GA 算法進(jìn)行了超參數(shù)尋優(yōu)，但適應(yīng)性較差；BP 算法由于不同層數(shù)神經(jīng)元的連接，容易造成參數(shù)膨脹，出現(xiàn)過(guò)擬合，且網(wǎng)絡(luò)深度太淺，數(shù)據(jù)區(qū)分度不高，因此適應(yīng)性最差。GS-SVC 算法的訓(xùn)練集準(zhǔn)確率為92.8%，測(cè)試集準(zhǔn)確率為87.6%，較GA-SVC 有所提高，證明GS 比GA 的尋優(yōu)方式更好，模型適應(yīng)性有所提高。GS-RFC 算法訓(xùn)練集和測(cè)試集的準(zhǔn)確率均超過(guò)了92%，對(duì)于不同水合物結(jié)構(gòu)類型，雖然不同氣體組分形成的水合物類型具有明顯的指向性，但具體輕微的含量差別對(duì)于水合物結(jié)構(gòu)的內(nèi)在影響不易區(qū)分，不同組分的區(qū)分度不大，隨機(jī)森林算法可以處理離散型或連續(xù)性數(shù)據(jù)，也可以處理異常數(shù)據(jù)的缺失，因此模型的適應(yīng)型最佳。

（3） 預(yù)測(cè)結(jié)果對(duì)比。不同分類模型的預(yù)測(cè)結(jié)果對(duì)比見(jiàn)表4。

表4 不同分類模型的預(yù)測(cè)結(jié)果

BP 算法的召回率、精確率和Fβ值最小，共出現(xiàn)了4 個(gè)分類錯(cuò)誤，BP 模型需要具有獨(dú)立的分類假設(shè)，因此分類效果不佳；GA-SVC 算法出現(xiàn)了3個(gè)分類錯(cuò)誤，數(shù)據(jù)集的數(shù)據(jù)具有離散非線性特點(diǎn)，對(duì)于除二元分類問(wèn)題外，表現(xiàn)效果不佳；GS-RFC算法的召回率、精確率和Fβ值分別為0.94、0.91和0.915，共出現(xiàn)了1 個(gè)分類錯(cuò)誤，在4 種算法模型中的準(zhǔn)確率最高，分類效果最好。

3.2 水合物生成溫度預(yù)測(cè)

水合物生成溫度預(yù)測(cè)模型與水合物結(jié)構(gòu)類型預(yù)測(cè)模型相比，輸入變量增加了壓力、抑制劑和鹽類等參數(shù)。超參數(shù)尋優(yōu)過(guò)程與分類模型相同，其中Mt=3，Nt=300。選取經(jīng)網(wǎng)格搜索優(yōu)化后的支持向量機(jī)回歸模型（GS-SVR） 和基于徑向基核函數(shù)的支持向量機(jī)模型（RBF-SVR） 與經(jīng)網(wǎng)格搜索優(yōu)化后的隨機(jī)森林回歸模型（GS-RFR） 進(jìn)行對(duì)比，預(yù)測(cè)結(jié)果見(jiàn)圖2、圖3、表5。其中GS-SVR 和RBF-SVR 的平均絕對(duì)誤差和均方根誤差較大，預(yù)測(cè)模型效果較差，而GS-RFR 預(yù)測(cè)模型的誤差最小，平均絕對(duì)誤差為1.91%，均方根誤差為0.35。主要是由于SVR 模型雖然是將非線性的數(shù)據(jù)映射到高維空間后再進(jìn)行分類處理，但仍然屬于淺層模型，對(duì)于區(qū)分度不大、維度較高的數(shù)據(jù)集預(yù)測(cè)效果不好，而隨機(jī)森林算法抗過(guò)擬合能力較強(qiáng)，對(duì)于預(yù)測(cè)不平衡的數(shù)據(jù)集誤差很小，屬于無(wú)偏估計(jì)數(shù)據(jù)驅(qū)動(dòng)的深層模型，更適合水合物生成條件的預(yù)測(cè)。從相關(guān)系數(shù)R 來(lái)看，GS-RFR預(yù)測(cè)模型的相關(guān)系數(shù)最大為0.987，證明預(yù)測(cè)結(jié)果與實(shí)際值更接近。

圖2 水合物生成溫度預(yù)測(cè)結(jié)果

圖3 水合物生成溫度預(yù)測(cè)結(jié)果相對(duì)誤差

表5 預(yù)測(cè)結(jié)果評(píng)價(jià)

4 結(jié)論

（1） 通過(guò)交叉驗(yàn)證，GS-RFC 算法訓(xùn)練集和測(cè)試集的準(zhǔn)確率均超過(guò)了92%，模型的適應(yīng)性最佳，召回率、精確率和Fβ值分別為0.94、0.91 和0.915，共出現(xiàn)了1 個(gè)分類錯(cuò)誤，分類效果最好。

（2） 對(duì)于水合物形成溫度的預(yù)測(cè)，GS-RFR預(yù)測(cè)模型的誤差最小，平均絕對(duì)誤差為1.91%，均方根誤差為0.35，相關(guān)系數(shù)為0.987，模型的準(zhǔn)確性和保守性最好，屬于無(wú)偏估計(jì)數(shù)據(jù)驅(qū)動(dòng)的深層模型，證明機(jī)器學(xué)習(xí)算法可為水合物的相關(guān)研究提供理論依據(jù)和實(shí)際參考。

（3） 研究表明，可通過(guò)預(yù)測(cè)結(jié)果判斷是否采取添加抑制劑或電加熱等措施。