原油快速評價技術(shù)的應(yīng)用研究

李敬巖, 褚小立, 田松柏

(中國石化 石油化工科學(xué)研究院, 北京 100083)

原油快速評價技術(shù)的應(yīng)用研究

李敬巖, 褚小立, 田松柏

(中國石化 石油化工科學(xué)研究院, 北京 100083)

為了快速準(zhǔn)確測量原油的密度、酸值、硫含量等重要性質(zhì),采用紅外光譜技術(shù)結(jié)合化學(xué)計量學(xué)定量校正算法建立校正模型。結(jié)果表明,采用偏最小二乘算法(PLS)建模預(yù)測原油密度、酸值和硫含量的標(biāo)準(zhǔn)偏差分別為0.0085 g/cm3、0.31 mgKOH/g和0.35%,預(yù)測的重復(fù)性與再現(xiàn)性等同或優(yōu)于標(biāo)準(zhǔn)方法,具有測量快速、操作簡單、無需預(yù)處理、重復(fù)性好等優(yōu)點。

原油; 偏最小二乘(PLS); 紅外光譜(MIR); 快速評價

原油評價在原油開采、原油貿(mào)易和原油加工等方面發(fā)揮著十分重要的作用。盡管現(xiàn)已建立了一套較為完整的原油評價方法,但這些方法分析時間長、工作量大、成本高,遠(yuǎn)不能滿足實際應(yīng)用的需要。目前,也存在較為完善的原油評價數(shù)據(jù)庫,并不斷得到更新和擴(kuò)充,但原油的性質(zhì)會隨時間發(fā)生較大變化,且在儲運(yùn)過程中經(jīng)常發(fā)生不同種類原油混合的情況。此外,機(jī)會原油交易逐漸頻繁,進(jìn)口原油普遍含硫量較高,高硫原油在管輸及裝置進(jìn)料時[1-2]會產(chǎn)生明顯腐蝕。因此,快速檢測原油性質(zhì),對于優(yōu)化進(jìn)料、調(diào)和、提高企業(yè)效率具有重要意義。針對以上情況,國內(nèi)外大型石化企業(yè)都正在基于多種現(xiàn)代儀器分析手段研發(fā)原油快評技術(shù),包括色譜、質(zhì)譜、核磁共振和紅外/近紅外光譜等[3-4]。

光譜分析技術(shù)是一種無損的快速檢測技術(shù)。由于其分析過程中無需預(yù)處理,大大降低了分析成本,一次光譜數(shù)據(jù)的采集可以同時分析多個指標(biāo)。近些年來,隨著計算機(jī)技術(shù)與化學(xué)計量學(xué)的發(fā)展,紅外、近紅外光譜技術(shù)與多元分析方法的結(jié)合已經(jīng)被廣泛用于各個領(lǐng)域[5-6]。紅外、近紅外光譜技術(shù)具有分析快、精密度高、操作簡單等優(yōu)點,非常適合原油及油品的定量和定性分析[7],并已經(jīng)越來越多地用于測量石油產(chǎn)品的物性參數(shù),如汽油辛烷值[8]、烴族組成,預(yù)測生物柴油主要成分[9]等。建立穩(wěn)健的定量校正模型是紅外、近紅外光譜分析的核心之一。常用的多元校正方法有多元線性回歸[10](MLR) 、偏最小二乘(PLS)[11-13]、人工神經(jīng)網(wǎng)絡(luò)(ANN)[14]等。Aske等[15]曾用紅外和近紅外光譜結(jié)合多元校正方法分析了原油的四組分組成。而非線性校正算法有最小二乘支持向量機(jī)算法[16](LSSVM)、核偏最小二乘[17](KPLS)等。

筆者采用透射中紅外光譜結(jié)合PLS算法,建立原油中密度、酸值、硫含量等8種原油評價數(shù)據(jù)的預(yù)測模型,以期實現(xiàn)原油的快速評價,為確定原油加工方案和優(yōu)化生產(chǎn)決策及時獲得評價數(shù)據(jù)提供一種簡捷的方法。

1 實驗部分

1.1 樣品

收集了國內(nèi)外具有代表性的原油樣本200種,基本覆蓋了世界主要原油產(chǎn)區(qū)的原油品種,其中石蠟基原油50個、中間基原油82個、環(huán)烷基原油32個、環(huán)烷-中間基原油15個、石蠟-中間基原油9個、中間-石蠟基原油14個,密度分布范圍0.7687～1.009 g/cm3。分別采用標(biāo)準(zhǔn)方法測定所收集原油的密度、酸值以及硫質(zhì)量分?jǐn)?shù)。

1.2 儀器與光譜采集

采用某型手持式紅外光譜儀,附件為100 μm ZnSe透射液體池。光譜采集范圍4750～950 cm-1,陣列檢測器,累積掃描次數(shù)20 次。所有樣本的紅外光譜均在常溫下采集,并實時扣除H2O和CO2的干擾。測量前,原油樣品需要鋪滿整個樣品池表面,對于黏度較大、不易鋪展的原油可以用棉簽將其涂抹在晶體表面。測量結(jié)束后,先用脫脂棉將樣品池中的原油大部分抹去,然后再用石油醚清洗數(shù)次,直至通過背景測試。

1.3 數(shù)據(jù)處理

采用“RIPP化學(xué)計量學(xué)光譜分析軟件3.0”對以標(biāo)準(zhǔn)方法測得的原油樣品的紅外光譜及其密度、酸值和硫質(zhì)量分?jǐn)?shù)進(jìn)行編輯,生成標(biāo)準(zhǔn)矩陣式數(shù)據(jù)庫。用Kennard-Stone(K-S)方法將原油樣品分為校正集和驗證集,其中校正集樣本150個,驗證集樣本50個。所用的化學(xué)計量學(xué)方法均采用MATLAB 6.5編寫,程序運(yùn)行平臺為ThinkPad T440p,i5(2.50 GHz),4 GB RAM。

2 結(jié)果與討論

2.1 紅外光譜預(yù)處理及參數(shù)選擇

圖1為RIPP原油庫中106#原油的原始紅外光譜。圖1中,1800～2400 cm-1范圍內(nèi)因為衍射效果比較嚴(yán)重,無法利用而予以剔除,波數(shù)4000 cm-1以后超出中紅外傳統(tǒng)的區(qū)間范圍也一并剔除,余下的作為最終建模選擇的光譜范圍,包含93個光譜變量。

圖1 106#原油的紅外光譜

在建立校正模型前需要對光譜進(jìn)行預(yù)處理,常用的算法有微分、歸一化等。一般使用微分處理以消除樣品顏色、溫度、基線漂移等因素的影響,但由于本研究得到的光譜變量數(shù)目較少,且缺乏連續(xù)性,微分不能直接用于離散波長光譜消除基線漂移。為此,首先采用3次樣條插值法對離散波長光譜進(jìn)行擬合,然后對擬合光譜進(jìn)行一階微分處理,再從微分光譜中取出對應(yīng)于原離散波長光譜數(shù)值的數(shù)值,構(gòu)成原光譜的微分光譜,從而實現(xiàn)離散波長光譜的基線漂移校正。

模型建立前,選取原油樣本分別構(gòu)成校正集和驗證集。所有樣本均在室溫下采集光譜,通過校正標(biāo)準(zhǔn)偏差(RMSECV)和預(yù)測標(biāo)準(zhǔn)偏差(RMSEP)來評價模型的準(zhǔn)確性。RMSECV和 RMSEP分別由式(1)和式(2)計算。

(1)

(2)

式(1)、(2)中,n為校正集的樣本總數(shù);m為預(yù)測集的樣本數(shù);yi,actual、yj,actual分別為通過標(biāo)準(zhǔn)方法得到校正集樣本和預(yù)測集樣本的實測值;yi,predicted為交互驗證預(yù)測值;yj,predicted為紅外光譜預(yù)測值。

2.2 校正集樣本的訓(xùn)練

原油成分極其復(fù)雜,其在中紅外光譜范圍內(nèi)(4000～400 cm-1),尤其是指紋區(qū)內(nèi)信息非常豐富,但譜帶重疊嚴(yán)重。筆者通過優(yōu)選后的區(qū)域參與建立校正模型。校正集樣本的選擇與分布對模型影響較大,校正集樣本需覆蓋預(yù)測集樣本的濃度范圍,同時分布盡量均勻。將經(jīng)過優(yōu)化后的光譜區(qū)間進(jìn)行一階微分處理后的吸光度值形成吸光度矩陣(X),相應(yīng)樣品用標(biāo)準(zhǔn)方法測定密度、酸值和硫含量等分別組成濃度矩陣(Y),然后用PLS法分別建立原油性質(zhì)的校正模型。

使用留一法對校正集樣本進(jìn)行交互驗證,PLS通過交互驗證的預(yù)測殘差平方和與主因子數(shù)的關(guān)系,確定最佳主因子并建立校正模型。圖2為酸值模型交互驗證的預(yù)測殘差平方和(PRESS)曲線。由圖2可見,當(dāng)主因子數(shù)為8時,PLS交互驗證標(biāo)準(zhǔn)偏差(RMSECV)達(dá)到最小并趨于穩(wěn)定,因此確定光譜最佳主因子數(shù)為8。其他如殘?zhí)恐怠⒌?、蠟含量、膠質(zhì)和瀝青質(zhì)含量的模型參數(shù)的最佳主因子數(shù)和相關(guān)系數(shù)列于表1。

圖2 原油酸值模型交互驗證的PRESS曲線

2.3 原油性質(zhì)的PLS預(yù)測分析

為了保證紅外光譜方法的準(zhǔn)確性,用建立的校正模型對驗證集的50個原油樣本進(jìn)行預(yù)測。圖3為原油性質(zhì)的PLS預(yù)測值和實際值的相關(guān)性,其中紅色為校正集樣本,黑色為預(yù)測集樣本。從圖3(a)可以看出,對模型外的樣本依然有很好的預(yù)測能力。從

表1 原油性質(zhì)的PLS模型參數(shù)的最佳主因子數(shù)和相關(guān)系數(shù)

圖3 中紅外光譜方法與標(biāo)準(zhǔn)方法測定200種原油樣品性質(zhì)的相關(guān)性比較

原油密度的交互驗證模型結(jié)果可以得出,相關(guān)系數(shù)R2為0.9541,交互驗證得到的RMSECV為0.0109 g/cm3,預(yù)測標(biāo)準(zhǔn)偏差(RMSEP)為0.0090 g/cm3,達(dá)到了標(biāo)準(zhǔn)方法的誤差要求;從原油酸值的交互驗證模型結(jié)果可以得出,相關(guān)系數(shù)R2為0.9562,交互驗證得到的RMSECV為0.39 mgKOH/g,RMSEP為0.31 mgKOH/g;從原油硫含量的交互驗證結(jié)果可以得出,相關(guān)系數(shù)R2為0.9721,交互驗證得到的RMSECV為0.31%,RMSEP為0.35%。因此,采用PLS算法結(jié)合紅外光譜預(yù)測原油密度、酸值和硫含量等主要性質(zhì)的校正模型完全可行,預(yù)測結(jié)果基本滿足了快速分析和過程分析的要求。其他如殘?zhí)恐怠⒌?、蠟含量、膠質(zhì)和瀝青質(zhì)含量的模型訓(xùn)練與預(yù)測結(jié)果列于表2。

2.5 紅外光譜分析方法測定原油性質(zhì)的重復(fù)性

使用透射光譜測量方式具有較好的穩(wěn)定性,為了檢驗紅外光譜分析方法的重復(fù)性,隨機(jī)選取3個原油樣本重復(fù)測定8次紅外光譜,并調(diào)用建立的原油密度、酸值和硫含量PLS預(yù)測模型計算,結(jié)果列于表3。由表3可知,紅外光譜分析方法測量原油密度、酸值和硫含量的標(biāo)準(zhǔn)偏差均小于相對應(yīng)的標(biāo)準(zhǔn)方法的,說明采用透射紅外光譜方法具有較高的精密度,可以用于原油的快速評價。

表2 原油性質(zhì)的PLS模型訓(xùn)練和預(yù)測效果

MIR—By MIR prediction

表3 紅外光譜分析方法測定原油性質(zhì)的重復(fù)性

SD—Standard deviation

3 結(jié) 論

(1)采用透射測量方式測定了200種原油的紅外光譜,并使用偏最小二乘方法建立了測定原油密度、酸值和硫含量等重要性質(zhì)的非線性定量校正模型。

(2) 對于原油特別是深色、黏稠的重質(zhì)原油,紅外光譜的測量存在進(jìn)樣和清洗的困難,而翻轉(zhuǎn)液體池附件可以有效解決,且光程長、重復(fù)性好。

(3)PLS方法可以準(zhǔn)確預(yù)測原油密度、酸值、硫含量、殘?zhí)恐?、氮含量、蠟含量、膠質(zhì)和瀝青質(zhì)含量,交互驗證標(biāo)準(zhǔn)偏差分別為0.0109 g/cm3、0.39 mgKOH/g、0.31/%、1.02%、0.05%、3.3%、2.5%和2.2%,預(yù)測標(biāo)準(zhǔn)差分別為0.0090 g/cm3、0.31 mgKOH/g、0.35%、0.88%、0.06%、2.9%、3.4%和1.8%,均接近于標(biāo)準(zhǔn)方法的誤差要求。

(4) 與標(biāo)準(zhǔn)方法相比,建立的方法具有測量快速、操作簡便等優(yōu)點,其結(jié)果滿足快速分析和過程控制分析的需要,在一定場合可以替代標(biāo)準(zhǔn)方法用于原油密度、酸值和元素硫含量的快速預(yù)測。

(5)小型紅外光譜儀雖然分辨率低、信息有限,但通過結(jié)合化學(xué)計量學(xué)算法也能滿足原油重要性質(zhì)快速評價的要求。

[1] 吳考民, 王可中, 石鑫, 等. 含硫原油對輸油管道的腐蝕性[J]. 油氣儲運(yùn), 2010, 29(8): 616-619. (WU Kaomin, WANG Kezhong, SHI Xin, et al.Corrosivity of sour crude on oil pipelines[J]. Oil & Gas Storage and Transportation, 2010, 29(8): 616-619.)

[2] 羅雄, 黃忠橋, 呂瑞典, 等. 加工高硫高酸原油引起的設(shè)備腐蝕問題探討[J]. 化學(xué)工程與裝備, 2010,(2): 59-62.(LUO Xiong, HUANG Zhongqiao, Lü Ruidian, et al. Research on corrosion in highly sulfur & acid crude oil processing[J]. Chemical Engineering & Equipment, 2010, (2): 59-62.)

[3] YUAN Hongfu, CHU Xiaoli, LI Haoran, et al. Determination of multi-properties of residual oils using mid-infrared attenuated total reflection spectroscopy[J]. Fuel, 2006, 85(12-13): 1720-1728.

[4] LI Jingyan, CHU Xiaoli, TIAN Songbai. Research on determination of total acid number of petroleum using mid-infrared attenuated total reflection spectroscopy[J]. Energy Fuels, 2012, 26(9):5633-5637.

[5] LIU F, ZHANG F, JIN Z L, et al. Determination of acetolactate synthase activity and protein content of oilseed rape(Brassica napus L.)leaves using visible/near infrared spectroscopy[J]. Analytica Chimica Acta,2008, 629(1-2): 56-65.

[6] KEMENY G J. Handbook of Near-Infrared Analysis[M]. New York: Marcel Dekker, 2001: 1-6.

[7] 褚小立, 許育鵬, 陸婉珍. 用于近紅外光譜分析的化學(xué)計量學(xué)方法研究與應(yīng)用進(jìn)展[J]. 分析化學(xué), 2008, 36(5): 702-709. (CHU Xiaoli, XU Yupeng, LU Wanzhen. Research and application progress of chemometrics methods in near infrared spectroscopic analysis[J]. Chinese Journal of Analytical Chemistry, 2008, 36(5): 702-709.)

[8] KELLYJ J. CALLIS J B. Nondestructive analytical procedure for simultaneous estimation of the major classes of hydrocarbon constituents finished gasolines[J]. Analytical Chemistry, 1990, 62(14): 1444-1451.

[9] PIMENTELA M F, RIBEIROB M G S, ROSENIRA S, et al. Determination of biodiesel content when blended with mineral diesel fuel using infrared spectroscopy and multivariate calibration[J]. Microchemical Journal, 2006, 82(2): 201-206.

[10] ANDREAS A K, NIKOS P. Autoregressive modeling of near-IR spectra and MLR to predict RON values of gasolines[J]. Fuel, 2010, 89(1): 158-161.

[11] DU Y P, LIANG Y Z, JIANG J H, et al. Spectral regions selection to improve prediction ability of PLS models by changeable size moving window partial least squares and searching combination moving window partial least squares[J]. Analytica Chimica Acta, 2004, 501(2): 183-191.

[12] 褚小立, 許育鵬, 陸婉珍. 偏最小二乘法方法在光譜定性分析中的應(yīng)用研究[J]. 現(xiàn)代儀器, 2007, (5): 13-15.(CHU Xiaoli, XU Yupeng, LU Wanzhen. The study of use of partial least squares in spectroscopy qualitative analysis[J]. Modern Instruments, 2007, (5): 13-15.)

[13] MARCELO M S, RONEI J P. N-way PLS applied to simultaneous spectrophotometric determination of acetylsalicylic acid, paracetamol and caffeine[J]. Journal of Pharmaceutical and Biomedical Analysis, 2004, 34(1): 27-34.

[14] 王艷斌, 劉偉, 袁洪福, 等. 人工神經(jīng)網(wǎng)絡(luò)在近紅外分析方法中的應(yīng)用及深色油品的分析[J]. 石油煉制與化工, 2002, 33(7): 62-66.(WANG Yanbin, LIU Wei, YUAN Hongfu, et al. Application of artificial neural network in near-infrared analysis and dark-colored oil analysis[J]. Petroleum Processing and Petrochemicals, 2002, 33(7): 62-66.)

[15] ASKE N, KALLEVIK H, SJOBLOM J. Determination of saturate, aromatic,resin, and asphaltenic (SARA) components in crude oils by means of infrared and near-infrared spectroscopy[J]. Energy Fuels, 2001, 15(5): 1304-1312.

[16] WU Di, HE Yong, FENG Shuijuan, et al. Study on infrared spectroscopy technique for fast measurement of protein content in milk powder based on LS-SVM[J]. Journal of Food Engineering, 2008, 84(1): 124-131.

[17] VAN DEN BROEK W, DERKS E, VAN DE VEN E W, et al. Plastic identification by remote sensing spectroscopic NIR imaging using kernel partial least squares (KPLS)[J]. Chemometrics and Intelligent Laboratory Systems, 1996, 35(2): 187-197.

Research on Crude Oil Fast Evaluation

LI Jingyan, CHU Xiaoli, TIAN Songbai

(ResearchInstituteofPetroleumProcessing,SINOPEC,Beijing100083,China)

For the fast and accurate measurement of density, TAN and sulfur mass fraction, the mid-infrared spectroscopy and chemometric quantitative calibration algorithm were used to establish the calibration models with the partial least square method (PLS). The prediction standard errors (SEP) of the density, acid value and sulfur mass fraction for crude oil were 0.0085 g/cm3,0.31 mgKOH/g and 0.35%, respectively, very close to those determined by standard methods. Compared with standards, this method was provided with the advantages of high-speed, simplicity, no-pretreatment and good-repeatability.

crude oil; partial least square(PLS); mid-infrared spectroscopy(MIR); fast-evaluation

2014-09-02

李敬巖，男，高級工程師，博士，從事紅外光譜分析與原油快速評價方面的研究 ；Tel：010-82368342；E-mail: lijy.ripp@sinopec.com

1001-8719(2015)06-1376-05

O657.33

A

10.3969/j.issn.1001-8719.2015.06.018