1、東北石油大學(xué)本科生畢業(yè)設(shè)計(jì)（論文）文獻(xiàn)翻譯PAGE PAGE 16模型預(yù)測(cè)油田水中溶解的碳酸鈣含量：壓力(yl)和溫度的影響 XXX 譯摘要(zhiyo)：油田(yutin)中水垢沉積會(huì)對(duì)儲(chǔ)層造成傷害、堵塞地層孔道、表面以及注入設(shè)備。碳酸鈣是水中最常見(jiàn)的結(jié)垢化合物之一，儲(chǔ)層產(chǎn)生的鹽水會(huì)使壓力和溫度降低，儲(chǔ)層壓力降低會(huì)使CaCO3的溶解度降低，進(jìn)而提高體系中碳酸鈣的飽和速率，而溫度下降會(huì)產(chǎn)生相反的結(jié)果。因此溫度和壓力一起作用的結(jié)果可能增加或減小CaCO3溶解度，用體系溫度的變化來(lái)指定其壓力的變化。因此，在石油生產(chǎn)系統(tǒng)中精確的預(yù)測(cè)方法的應(yīng)用備受關(guān)注。目前的研究重點(diǎn)是運(yùn)用基于最小二乘支持向量機(jī)(LS

2、SVM)預(yù)測(cè)模型來(lái)估計(jì)油田水中溶解碳酸鈣濃度的大小。用超優(yōu)化參數(shù)（r和C2）的遺傳算法（GA）嵌入到LSSVM模型，這種方法可簡(jiǎn)單準(zhǔn)確的預(yù)測(cè)油田鹵水中溶解碳酸鈣濃度的最小量。1.引言隨著油田鹵水壓力和溫度變化，氣體可能會(huì)從儲(chǔ)層到地表的運(yùn)動(dòng)，導(dǎo)致某些固體沉淀。為了保持注水井壓力平衡并將油運(yùn)移到生產(chǎn)井，有時(shí)需要將鹵水注入到儲(chǔ)層中，因此，過(guò)量的鹽垢可以沉積在儲(chǔ)層或井眼內(nèi)。對(duì)于大部分油田結(jié)垢多會(huì)發(fā)生在此過(guò)程中。碳酸鈣沉積通常是一個(gè)自發(fā)的過(guò)程，沉積形成的主要原因是二氧化碳從水相逸出，導(dǎo)致油氣層的壓力下降，該過(guò)程會(huì)除去了水中的碳酸，直到方解石溶解完全。在恒定二氧化碳分壓下，方解石的溶解性隨溫度的降低而降低

3、1-4。根據(jù)公式（1），碳酸鈣沉積垢來(lái)自碳酸鈣沉淀：Ca2+ + CO32- CaCO3 下面的公式為碳酸的電離式57：CO2 + H2O H2CO3 H2CO3 H+ + HCO3-HCO3- H+ + CO32-若要形成(xngchng)碳酸氫根離子和氫離子，碳酸要電離，因?yàn)樘妓岬牡谝浑婋x常數(shù)遠(yuǎn)大于它的第二(d r)電離常數(shù)，從碳酸第一電離離子化的氫離子與水中自由的碳酸根離子結(jié)合。此外，碳酸鈣沉淀的方程式可以(ky)說(shuō)明810：Ca(HCO3)2 CaCO3+ CO2+ H2O碳酸鈣的溶解度很大程度上取決于二氧化碳在水中的含量（即二氧化碳?xì)怏w逸出時(shí)所需最小的分壓）1012 。二氧化碳?xì)怏w注

4、入到儲(chǔ)層鹵水可產(chǎn)生酸性溶液，能溶解地層中的方解石，壓力的下降導(dǎo)致二氧化碳?xì)怏w從酸性溶液中逸出，并且在近井地帶的地層空隙和射孔眼處沉積碳酸鹽垢。而且，在井筒附近產(chǎn)生的垢會(huì)進(jìn)一步導(dǎo)致氣井產(chǎn)能的下降和更多結(jié)垢的產(chǎn)生13,14 。Bahadori借助范德蒙德行列式的相關(guān)性來(lái)研究估算油田鹵水溶解的碳酸鈣濃度15。本論文中涉及到的最小二乘支持向量機(jī)（LSSVM），是作為一種傳統(tǒng)簡(jiǎn)化的支持向量機(jī)（SVM）來(lái)預(yù)測(cè)油田鹵水中溶解碳酸鈣濃度量的大小。遺傳算法（GA）偶合來(lái)優(yōu)化主題，用于確定超參數(shù)。據(jù)我們所了解的，目前還沒(méi)有文獻(xiàn)記錄采用支持向量方法進(jìn)行天然氣產(chǎn)能建模。2.原理理論2.1 最小二乘支持向量機(jī)（LSSV

5、M）支持向量機(jī)（SVM）最初形成的是一種智能化的系統(tǒng)方法16。支持向量機(jī)已被廣泛應(yīng)用于分類(lèi)，回歸分析和模式識(shí)別，是一種非常有效的而方法17。支持向量機(jī)的原理是對(duì)于線(xiàn)性不可分的情況，通過(guò)使用非線(xiàn)性 HYPERLINK /view/21249.htm t _blank 映射算法將低維輸入空間線(xiàn)性不可分的樣本轉(zhuǎn)化為 HYPERLINK /view/2019275.htm t _blank 高維特征空間使其線(xiàn)性可分，從而使得高維特征空間采用線(xiàn)性算法對(duì)樣本的非線(xiàn)性特征進(jìn)行線(xiàn)性分析成為可能，并且基于結(jié)構(gòu)風(fēng)險(xiǎn)最小化理論之上在特征空間中建構(gòu)最優(yōu)分割超平面18。這種新方法基于結(jié)垢風(fēng)險(xiǎn)最小化理論（SRM）和統(tǒng)計(jì)學(xué)

6、理論（SLT）概念19。支持向量機(jī)通過(guò)求解二次規(guī)劃（QP）問(wèn)題來(lái)發(fā)現(xiàn)目標(biāo)函數(shù)的全局最小值，而不是獲得局部最優(yōu)解。與其他回歸技術(shù)不同，二次規(guī)劃問(wèn)題是一個(gè)凸函數(shù)20。然而，這種做法可能會(huì)耗時(shí)并且困難，因?yàn)樗枰鞣N非線(xiàn)性公式（QP），Suyken和Vandewalle 21 建議最小二乘支持向量機(jī)（LSSVM）作為支持向量機(jī)（SVM）的一個(gè)備用方法。SVM也具有LSSVM類(lèi)似的優(yōu)點(diǎn)。此外，需要一個(gè)二次規(guī)劃（QP）可代替多個(gè)非線(xiàn)性公式（線(xiàn)性規(guī)劃），這使得計(jì)算問(wèn)題簡(jiǎn)單且容易處理。LSSVM公式中的非線(xiàn)性函數(shù)可表示如下，給出訓(xùn)練集Xk，Yk，K=1，2，N，其中XkRn，是第i個(gè)輸入到輸入集合空間，Yk

7、R輸出值為特定輸入變量的對(duì)應(yīng)值（即Xk），N是指訓(xùn)練樣本的數(shù)目，通過(guò)一個(gè)非線(xiàn)性 HYPERLINK /view/21249.htm t _blank 映射，把 HYPERLINK /view/562162.htm t _blank 樣本空間 HYPERLINK /view/21249.htm t _blank 映射到一個(gè) HYPERLINK /view/2019275.htm t _blank 高維乃至無(wú)窮維的特征空間中，構(gòu)造的回歸模型如下：其中(qzhng)w表示權(quán)重向量，而b是一個(gè)偏置(pin zh)術(shù)語(yǔ)，上表“n”是數(shù)據(jù)(shj)區(qū)域的范圍，“nh”表示未知區(qū)域（高維和無(wú)窮維）的范圍20

8、。當(dāng)應(yīng)用LSSVM時(shí)會(huì)產(chǎn)生一個(gè)新的優(yōu)化問(wèn)題，該方法涉及的優(yōu)化問(wèn)題如下：根據(jù)下列等式進(jìn)行約束：其中r表示正則化參數(shù)的平衡訓(xùn)練誤差和模型闡述19，ek是回歸誤差。拉格朗日乘子理論構(gòu)造如下，以找到優(yōu)化問(wèn)題的約束條件：其中，ak拉格朗日乘子或支集值，為了獲得上述方程的解，對(duì)上述方程中的變量w，b，ek，ak微分：在除去變量w和e后，獲得Karush-Kuhn-Trucker系統(tǒng)如下：其中y=y1，yN T，1N=1，1 T，a=a1，aN T，I是一個(gè)單位矩陣，是內(nèi)核函數(shù)，必須滿(mǎn)足Mercerd的條件內(nèi)核函數(shù)三種典型的選擇方法如下所示：（1）線(xiàn)性 HYPERLINK /view/8075712.htm

9、 t _blank 核函數(shù)K(x,y)=xy；（2）多項(xiàng)式 HYPERLINK /view/8075712.htm t _blank 核函數(shù)K(x,y)=(xy)+1d；（3） HYPERLINK /view/1683917.htm t _blank 徑向(jn xin)基函數(shù)K(x,y)=exp(-|x-y|2/d2）由此LSSVM模型函數(shù)估計(jì)(gj)如下：其中(qzhng)（b，a）是式（14）線(xiàn)性系統(tǒng)的解。文獻(xiàn)中，一些綜合性描述支持向量機(jī)內(nèi)容是可用的16,21,23，并且對(duì)最小二乘支持向量機(jī)理論的解釋很詳細(xì)21,24，此外，劉等人2527，提供了LSSVM理論價(jià)值的解釋?zhuān)约癝VM和LS

10、SVM在有關(guān)情況下的區(qū)別，所以，感興趣的讀者可以參考上述文獻(xiàn)以了解更多的細(xì)節(jié)。2.2 遺傳算法遺傳算法（Genetic Algorithm）是模擬 HYPERLINK /view/6739.htm t _blank 達(dá)爾文 HYPERLINK /view/183146.htm t _blank 生物進(jìn)化論的自然選擇和 HYPERLINK /subview/8907/11410208.htm t _blank 遺傳學(xué)機(jī)理的生物進(jìn)化過(guò)程的計(jì)算 HYPERLINK /view/96500.htm t _blank 模型，是一種通過(guò)模擬自然進(jìn)化過(guò)程搜索 HYPERLINK /view/1009692.

11、htm t _blank 最優(yōu)解的方法。遺傳算法是從代表問(wèn)題可能潛在的解集的一個(gè) HYPERLINK /view/43530.htm t _blank 種群（population）開(kāi)始的，而一個(gè)種群則由經(jīng)過(guò) HYPERLINK /view/8563.htm t _blank 基因（gene）編碼的一定數(shù)目的個(gè)體(individual)組成。每個(gè)個(gè)體實(shí)際上是 HYPERLINK /view/6633.htm t _blank 染色體(chromosome)帶有特征的實(shí)體。染色體作為 HYPERLINK /view/232169.htm t _blank 遺傳物質(zhì)的主要載體，即多個(gè)基因的 HYPE

12、RLINK /view/15216.htm t _blank 集合，其內(nèi)部表現(xiàn)（即 HYPERLINK /view/142644.htm t _blank 基因型）是某種基因組合，它決定了個(gè)體的形狀的外部表現(xiàn)。因此，在一開(kāi)始需要實(shí)現(xiàn)從 HYPERLINK /view/281787.htm t _blank 表現(xiàn)型到基因型的 HYPERLINK /view/21249.htm t _blank 映射即 HYPERLINK /view/237708.htm t _blank 編碼工作。由于仿照基因編碼的工作很復(fù)雜，我們往往進(jìn)行簡(jiǎn)化，如 HYPERLINK /view/18536.htm t _bl

13、ank 二進(jìn)制編碼，初代種群產(chǎn)生之后，按照適者生存和優(yōu)勝劣汰的原理，逐代（generation）演化產(chǎn)生出越來(lái)越好的近似解，在每一代，根據(jù)問(wèn)題域中個(gè)體的 HYPERLINK /view/2167684.htm t _blank 適應(yīng)度（fitness）大小選擇（selection）個(gè)體，并借助于自然遺傳學(xué)的遺傳 HYPERLINK /view/53313.htm t _blank 算子（genetic operators）進(jìn)行組合交叉（crossover）和變異（mutation），產(chǎn)生出代表新的解集的種群。這個(gè)過(guò)程將導(dǎo)致種群像自然進(jìn)化一樣的后生代種群比前代更加適應(yīng)于環(huán)境，末代種群中的最優(yōu)個(gè)體

14、經(jīng)過(guò) HYPERLINK /view/657399.htm t _blank 解碼（decoding），可以作為問(wèn)題近似最優(yōu)解。1. 初始化：此階段由r和C2超參數(shù)來(lái)創(chuàng)造染色體，C2表示正規(guī)化因素，r代表RBF核函數(shù)的寬度，染色體A的定義為A=a1,a2，其中a1和a2與r和C2相對(duì)應(yīng)的，此外，選定種群數(shù)量大小為20的，最大的代分配給100，初始化代分配給0。2. 評(píng)價(jià)適應(yīng)度: 此階段，對(duì)每一個(gè)解(染色體)指定一個(gè)適應(yīng)度的值，根據(jù) HYPERLINK /view/1099373.htm t _blank 問(wèn)題求解的實(shí)際接近程度來(lái)指定(以便逼近求解問(wèn)題的答案)。測(cè)試數(shù)據(jù)的平均相對(duì)偏差為適應(yīng)度函數(shù)

15、。其中(qzhng)Caiactual表示(biosh)實(shí)際碳酸鈣濃度的大小，Caipredicted表示(biosh)用LSSVM模型預(yù)測(cè)碳酸鈣濃度的級(jí)數(shù)。3. 選擇：選擇的目的是把優(yōu)化的個(gè)體(或解)直接遺傳到下一代或通過(guò)配對(duì)交叉產(chǎn)生新的個(gè)體再遺傳到下一代。輪盤(pán)賭選擇法 （roulette wheel selection）是最簡(jiǎn)單也是最常用的選擇方法。圖1 遺傳算法的超參數(shù)選擇流程圖4. 交叉：所謂交叉是指把兩個(gè)父代個(gè)體的部分結(jié)構(gòu)加以替換重組而生成新個(gè)體的操作。生成每一對(duì)新的染色體的可能性為0.8，且執(zhí)行單點(diǎn)交叉，新染色體結(jié)合為新的個(gè)體后代。5. 變異：變異操作遵循交叉操作，并確認(rèn)是否有染色

16、體變異來(lái)進(jìn)行進(jìn)一步創(chuàng)造，在目前的進(jìn)展中，均勻變異概率為0.05，此進(jìn)化方法設(shè)計(jì)并執(zhí)行。6. 下一代：一代一代演化下去，直到達(dá)到期望的解為止。7. 終止條件：當(dāng)最優(yōu)個(gè)體的 HYPERLINK /view/2167684.htm t _blank 適應(yīng)度達(dá)到給定的 HYPERLINK /view/409216.htm t _blank 閾值，或者最優(yōu)個(gè)體的適應(yīng)度和群體適應(yīng)度不再上升時(shí)，或者迭代次數(shù)達(dá)到預(yù)設(shè)的代數(shù)時(shí)，算法終止，預(yù)設(shè)代數(shù)為100。3.方法論GA-LSSVM是由實(shí)驗(yàn)室1.8免費(fèi)工具箱和遺傳算法工具箱MATLABR2009a對(duì)LSSVM變量進(jìn)行優(yōu)化。在識(shí)別和收集數(shù)據(jù)之后，再定義LSSVM模

17、型(mxng)的輸入和輸出參數(shù)，在文獻(xiàn)數(shù)據(jù)庫(kù)的基礎(chǔ)上，油田鹵水中溶解碳酸鈣濃度量作為L(zhǎng)SSVM預(yù)測(cè)模型的輸出參數(shù)，該模型有兩個(gè)輸入?yún)?shù)：壓力和溫度。識(shí)別收集的實(shí)際數(shù)據(jù)庫(kù)分為兩個(gè)子集。首先在預(yù)測(cè)培訓(xùn)的過(guò)程中，利用整個(gè)數(shù)據(jù)的80%（63個(gè)數(shù)據(jù)點(diǎn)），50個(gè)等于數(shù)據(jù)線(xiàn)，剩余的13個(gè)組成樣本設(shè)置，來(lái)檢驗(yàn)上述建議(jiny)方法的評(píng)估能力和完整性。與其他內(nèi)核函數(shù)比較，RBF核的參數(shù)少并且有最佳的整體性能，因此RBF核函數(shù)在文獻(xiàn)中被視為一種合適(hsh)有效的內(nèi)核函數(shù)30,31。在利用LSSVM模型和RBF核函數(shù)開(kāi)發(fā)模型，在上述方程中超參數(shù)r2和C2規(guī)范很重要，因?yàn)檫@兩個(gè)參數(shù)在獲得高估計(jì)完整性和概括性的LS

18、SVM模型中有很重要的作用，其中r2為正則化因子，C為樣本方差20。根據(jù)Ahmadi等，其中包括Levenberg-Marquardt（LM）和Simplex SlmelatedAnnealingAlgvrithm(M-SIMPSON)不倡導(dǎo)不是基于逐代優(yōu)化技術(shù)的方法，因?yàn)橹С窒蛄繖C(jī)無(wú)法處理高非線(xiàn)性問(wèn)題?；谥鸫鷥?yōu)化技術(shù)即遺傳算法，已經(jīng)用于對(duì)超參數(shù)（r2和C2）的評(píng)估，AARD的測(cè)試數(shù)據(jù)為遺傳算法的適應(yīng)度，重復(fù)優(yōu)化過(guò)程才肯能確定，相應(yīng)的ARRD變量如圖2所示.在圖2中，36代后ARRD穩(wěn)定，并達(dá)到最佳態(tài)，通過(guò)這種方式，從油田鹵水溶解碳酸鈣濃度量獲得r2和C2級(jí)數(shù)值，分別為20.96137079

19、4和1072577.9326。圖2-ARRD變異隨代數(shù)的變化4.結(jié)果(ji gu)和討論在本研究中，確定(qudng)參數(shù)r2和均方誤差（MSE）評(píng)價(jià)預(yù)測(cè)(yc)模型有效性和準(zhǔn)確性的指標(biāo)，計(jì)算MSE和r2的表達(dá)式如下：其中，N表示數(shù)據(jù)點(diǎn)的數(shù)量，yiactual表示第i個(gè)數(shù)值，yipredicted表示第i個(gè)輸出模型值，yiactual平均表示實(shí)驗(yàn)觀(guān)察數(shù)據(jù)的平均值。表1中有MSE和r2值，用于預(yù)測(cè)訓(xùn)練，測(cè)試和數(shù)據(jù)分析。最后，如果得到r2和MSE值分別接近于1和0，GA-LSSVM預(yù)測(cè)估計(jì)結(jié)果令人滿(mǎn)意。表1-用于預(yù)測(cè)油田鹵水中溶解碳酸鈣的GA-LSSVM方法性能和優(yōu)化超參數(shù)模型中得到的結(jié)果，油田鹵

20、水中溶解碳酸鈣濃度量大小隨溫度的實(shí)際變化，如圖3所示，此外，油田鹵水中溶解碳酸鈣濃度量大小隨壓力的實(shí)際變化，如圖4所示，此外，油田鹵水中溶解碳酸鈣濃度量大小隨壓力和溫度變化的三維圖見(jiàn)附錄A，GA-LSSVM模型預(yù)測(cè)和實(shí)際油田鹵水中溶解碳酸鈣濃度量大小的比較，如圖5所示，此圖所表示的GA-LSSVM模型預(yù)測(cè)和實(shí)際油田鹵水中溶解碳酸鈣濃度量大小值相互覆蓋的現(xiàn)象。圖3-油田鹵水中溶解碳酸鈣濃度量大小(dxio)隨溫度的實(shí)際變化1,15.圖4-油田鹵水中溶解碳酸鈣濃度(nngd)量大小隨壓力的實(shí)際變化1,15圖5- GA-LSSVM模型預(yù)測(cè)和實(shí)際油田鹵水中溶解(rngji)碳酸鈣濃度量大小的比較油田鹵

21、水中溶解碳酸鈣濃度量大小的實(shí)際值與GA-LSSVM預(yù)測(cè)估計(jì)值關(guān)系，如圖6所示，油田鹵水中溶解碳酸鈣濃度量大小預(yù)測(cè)值分散在1:1線(xiàn)出，表明GA-optimized模型可以很好的預(yù)測(cè)油田鹵水中溶解碳酸鈣濃度量，描述比較LSSVM模型預(yù)測(cè)值和油田鹵水中溶解碳酸鈣濃度量實(shí)際值與溫度關(guān)系圖，如圖7所示，從向量模型所得輸出值和實(shí)際數(shù)據(jù)值有很高的相似度，為了證明，油田鹵水中溶解碳酸鈣濃度量預(yù)測(cè)值和實(shí)際值的相對(duì)偏差，如圖8所示，相對(duì)偏差在4.264%到4.576%，最低相對(duì)偏差為0%，平均相對(duì)偏差為1.7121%，這也證明了油田鹵水中溶解碳酸鈣濃度量預(yù)測(cè)值和實(shí)際值很接近。圖6-油田鹵水中溶解(rngji)碳酸

22、鈣濃度量大小的實(shí)際值與GA-LSSVM預(yù)測(cè)估計(jì)值關(guān)系圖7-LSSVM模型預(yù)測(cè)值和油田鹵水中溶解碳酸鈣濃度量實(shí)際(shj)值與溫度關(guān)系圖1,15圖8-油田(yutin)鹵水中溶解碳酸鈣濃度量預(yù)測(cè)值和實(shí)際值的相對(duì)偏差確定異常值的數(shù)學(xué)方法在解決(jiju)實(shí)用性問(wèn)題問(wèn)題中起著至關(guān)重要的作用，離群值表明單個(gè)的數(shù)據(jù)有別于數(shù)據(jù)庫(kù)中部分主要數(shù)據(jù)36,37。結(jié)果表明，需要一個(gè)可行的實(shí)際數(shù)據(jù)用于生產(chǎn)分析(fnx)，因?yàn)閮?yōu)化方法估計(jì)能力具有不確定性。為了得到這個(gè)結(jié)果，可應(yīng)用杠桿統(tǒng)計(jì)的方法37,38.異常或可疑數(shù)據(jù)圖形通過(guò)Williams plot 獲得，決定H值，在之前的文獻(xiàn)中可以找到詳細(xì)的數(shù)學(xué)背景和計(jì)算過(guò)程36

23、38. 如圖9所示的為Williams plot 與向量機(jī)方法結(jié)果的比對(duì)，在所有的數(shù)據(jù)樣本中，-3R3和0H0.18，對(duì)進(jìn)化方法認(rèn)證確認(rèn)，此外，說(shuō)明了整個(gè)數(shù)據(jù)集的數(shù)據(jù)分布在合理范圍內(nèi)。圖9-選擇(xunz)出在測(cè)量油田(yutin)鹵水中溶解碳酸鈣濃度量并且應(yīng)用主要的建議方法分析生產(chǎn)油井的水，pH值為0.185.結(jié)論(jiln)根據(jù)研究，可得到以下結(jié)論：1. 據(jù)已經(jīng)獲得的實(shí)際數(shù)據(jù)庫(kù)，使用LSSVM技術(shù)和RBF核函數(shù)來(lái)預(yù)測(cè)油田鹵水中溶解碳酸鈣濃度。2. 適合的超參數(shù)r2和C2對(duì)LSSVM預(yù)測(cè)和泛化能力很重要，使用遺傳算法（GA），從油田鹵水溶解碳酸鈣濃度獲得r2和C2級(jí)數(shù)值，分別為20.9613

24、70794和1072577.9326。3. LSSVM技術(shù)與遺傳算法（GA）的交叉在預(yù)測(cè)油田鹵水溶解碳酸鈣濃度取得了很好的結(jié)果，LSSVM模型預(yù)測(cè)值和油田鹵水中溶解碳酸鈣濃度實(shí)際值具有很好的一致性，r2和MSE值分別為0.944和0.0000208。4. LSSVM技術(shù)具有很多優(yōu)良的特性，比如效率，易于計(jì)算，這可以帶來(lái)一個(gè)很有吸引力的高度非線(xiàn)性模型，如油田海水淡化系統(tǒng)。附錄A本部分闡述油田鹵水中溶解碳酸鈣濃度量隨溫度和壓力的變化三維圖，如圖A。圖A-油田鹵水(lshu)中溶解碳酸鈣濃度量隨溫度和壓力的變化三維圖參考文獻(xiàn)1 Davies M, Scott PJB. Oilfield water

25、technology. Houston, TX: NACE International, 2006.2 Honsy R, Desouky SEM, Ramzi M, Abdel-Moghny T, El-Dars FMS, Farag AB.J Dispers Sci Technol, 2009, 30: 20412.3 Vralstad T, Fossen M, Sjolom J, Randhol P. J Dispers Sci Technol 2008, 29: 4406.4 Adenidi OD, Nwalor JU, Ako CT. J Dispers Sci Technol 200

26、8, 29: 3625.5 Bang J-H, Jang YN, Kim W, Song KS, Jeon CW, Chae SC, et al. Chem Eng J 2011, 174: 41320.6 Bahadori A. SPE Prod Facil Construct, 2011, 6: 15865.7 Zhang G, Taberner C, Cartwright L, Xu T. SPE J, 2011,16: 95967.8 Soong Y, Dilmore R, Hedges S, Griffith C. CO2 Capture=sequestration via brin

27、e and caustic material. In: Conference Proceedings, AIChE Annual Meeting, 2009.9 Berger C, Dandeu A, Carteret C, Humbert B, Muhr H, Plasari E, et al. Chem EngTrans, 2009, 17: 6816.10 Duguid A, Scherer GW. Int J Greenhouse Gas Control, 2010, 4: 54660.11 Ramstad K, Tydal T, Askvik KM, Fotland P. SPE J

28、, 2005, 10: 36373.12 Dyer SJ, Graham GM. J Petrol Sci Eng 2002, 35: 95107.13 Xiaoyan L, Jungang L, Qianya Z, Jinlai F, Yingli L, Jingxin S. J Petrol Sci Eng, 2009, 66: 1614.14 Al-Anezi K, Hilal N. Desalination 2007, 204: 385402.15 Bahadori A, Zendehboudi S. J Dispers Sci Technol 2013, 34(6): 7939.16

29、 Vapnik V. Statistical learning theory. New York. Wiley, 1998.17 Cortes C, Vapnik V. Machine Learn 1995, 20: 27397.18 Baylar A, Hanbay D, Batan M. Expert Syst Appl, 2009, 36: 836874.19 Mehdizadeh B, Movagharnejad K. Chem Eng Res Des. 2011, 89: 24207.20 Vong CM, Wong PK, Li YP. Eng Appl Artif Intell

30、2006, 19:27787.21 Suykens J, Vandewalle J. Neural Process Lett, 1999, 9: 293300.22 Li C, Zhu X, Cao G, Sui Sh, Hu M. J Power Sources. 2008, 175: 30316.23 Burges C. Data Mining Knowl Discuss. 1998, 2:12167.24 Suykens J, Gestel TV, Brabanter JD, Moor BD, Vandewalle J. Least squares support vector machines. Singapore: World Scientific, 2002.25 Liu H, Yao X, Zhang R, Liu M, Hu Z, Fan B. J Phys Chem B 2005, 109: 6571.26 Liu HX, Yao XJ, Zhang RS, Liu MC, Hu ZD, Fan BT. J Comput Aided Mol Des. 2005, 19: 499508.27 Li J, Liu H, Yao X, Liu M, Hu Z, Fan B. Chemomet Intell Lab Syst, 2007, 87: 13946.28 Ah