中俄輸油管道泄露過(guò)程中原油于土壤中的遷移研究

王學(xué)力，王榕，吳官生，高鵬，王博，王東坡

(1.天津大學(xué)材料科學(xué)與工程學(xué)院，天津 300072；2.中國(guó)石油管道公司大慶(加格達(dá)奇)輸油氣分公司，黑龍江 加格達(dá)奇 165000；3.四川農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院，四川 成都 611130)

?

中俄輸油管道泄露過(guò)程中原油于土壤中的遷移研究

王學(xué)力1,2，王榕3，吳官生2，高鵬2，王博2，王東坡1

(1.天津大學(xué)材料科學(xué)與工程學(xué)院，天津 300072；2.中國(guó)石油管道公司大慶(加格達(dá)奇)輸油氣分公司，黑龍江 加格達(dá)奇 165000；3.四川農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院，四川 成都 611130)

摘要：為了探尋中俄輸油管道泄露原油過(guò)程中原油在土壤中的遷移規(guī)律，采用箱體動(dòng)態(tài)模擬和球形積分模型對(duì)此進(jìn)行研究，實(shí)驗(yàn)土壤為沿線的粉質(zhì)粘土和粘土混碎石.結(jié)果表明，原油在不同土體各方向遷移速率的大小皆為：垂直向下方向(Z負(fù)軸)>管線方向(Y正軸/Y負(fù)軸)>水平垂直管線方向(X正軸/X負(fù)軸)>垂直向上方向(Z正軸).分析土壤含水率對(duì)原油遷移距離的影響，結(jié)果表明，含水率在6.09%15.02%范圍內(nèi)時(shí)，原油的遷移距離隨著土壤含水率的增大先增大后減小，含水率為10.33%時(shí)遷移距離最大.箱體動(dòng)態(tài)模擬實(shí)驗(yàn)結(jié)果和球形積分模型的預(yù)測(cè)結(jié)果一致，原油的遷移速率大小為：粘土混碎石>粉質(zhì)粘土，說(shuō)明了球形積分模型的合理性.

關(guān)鍵詞：管道泄漏；遷移；原油；箱體動(dòng)態(tài)模擬

0引言

近年來(lái)，輸油管道原油泄漏事件的頻繁發(fā)生，不僅造成了嚴(yán)重的資源浪費(fèi)和經(jīng)濟(jì)損失，并且對(duì)生態(tài)環(huán)境帶來(lái)了巨大的威脅.中俄輸油管道是一條埋藏式的年輕管道，且該管道穿越的地區(qū)地形復(fù)雜、生態(tài)環(huán)境脆弱，因此研究中俄輸油管道原油在土壤中的遷移規(guī)律至關(guān)重要.目前，國(guó)內(nèi)外研究原油及其污染物在土壤中遷移規(guī)律的方法主要有：現(xiàn)場(chǎng)取樣分析法[1-4]、土壤過(guò)濾法[5-7]、淋濾法[8-11]、離心模型法[12-13]等，這些方法主要針對(duì)已經(jīng)脫離管道存在于地表上或土壤中的原油遷移，而對(duì)于在管道泄漏過(guò)程中原油在土壤中遷移規(guī)律的研究方法較少.為了在發(fā)現(xiàn)原油泄漏時(shí)能及時(shí)正確地采取處理措施減少損失和危害，本文中采用箱體動(dòng)態(tài)模擬實(shí)驗(yàn)研究中俄輸油管道泄漏原油過(guò)程中原油在自然土體中的遷移規(guī)律.

1材料方法

1.1實(shí)驗(yàn)材料1)油品.實(shí)驗(yàn)用油采用中俄輸油管道林源進(jìn)站俄原油.經(jīng)檢測(cè)分析，原油的密度為0.843 3 g/cm3，凝點(diǎn)小于-25 ℃，水含量為0.064%，蠟含量為4.35%，膠質(zhì)和瀝青含量為6.67%.2)土樣.實(shí)驗(yàn)土樣粉質(zhì)粘土和粘土混碎石采集于漠大線加格達(dá)奇地段，其中粘土混碎石由花崗巖與砂巖發(fā)育而來(lái).土樣的基本理化參數(shù)見(jiàn)表1，表中的容重及含水率為在實(shí)驗(yàn)箱體中采樣并測(cè)定的結(jié)果.3)熒光燈.采用LGY JL-101型手電筒式熒光燈，燈光顏色為藍(lán)色；經(jīng)檢測(cè)，原油在藍(lán)色燈光的照射下顯現(xiàn)出鮮明的暗黃色，原油落在土壤中后可以明顯地區(qū)分土壤與油，并且可以清晰分辨出原油的遷移路徑.

表1　實(shí)驗(yàn)土樣的基本理化性質(zhì)

1.2實(shí)驗(yàn)方法

1.2.1箱體動(dòng)態(tài)模擬實(shí)驗(yàn)采用邊長(zhǎng)為20 cm的立方體木箱進(jìn)行箱體動(dòng)態(tài)模擬實(shí)驗(yàn).將土樣每2 cm一層壓實(shí)，把土壤一層一層裝入箱體，同時(shí)將輸油管鋪設(shè)在箱體中心固定，裝好的土壤應(yīng)達(dá)到均勻密實(shí)、接近自然土體.箱體裝置見(jiàn)圖1.然后將原油從豎直輸油管連續(xù)地勻速加入，在實(shí)驗(yàn)時(shí)間結(jié)束時(shí)停止加入，立馬打開(kāi)排油閥門(mén)，待輸油管內(nèi)殘留的油全部流出后，馬上拔出輸油管打開(kāi)箱體；用削土刀由進(jìn)油口方向往排油口方向根據(jù)一定的距離迅速切下土層，用熒光手電筒觀察切下的每個(gè)剖面上原油滲透的痕跡，用刻度尺量取每個(gè)剖面上原油斑邊際距離剖面中心的距離并記錄.經(jīng)過(guò)預(yù)實(shí)驗(yàn)，設(shè)置粉質(zhì)粘土的實(shí)驗(yàn)時(shí)間梯度為：1、2、3、4、5 min；粘土混碎石的實(shí)驗(yàn)時(shí)間梯度為：15、30、60、120 s；每一時(shí)間梯度進(jìn)行3次重復(fù)試驗(yàn)并取其平均值.

圖1　箱體模擬實(shí)驗(yàn)裝置

1.2.2原油遷移距離預(yù)測(cè)根據(jù)在土體剖面XOZ(X、Y、Z軸方向見(jiàn)圖1)上采集的實(shí)驗(yàn)數(shù)據(jù)，利用Origin 8.5軟件，采用高斯方程(Gauss)進(jìn)行數(shù)據(jù)擬合，得到每個(gè)XOZ平面的橢圓面積曲線方程，然后對(duì)該曲線沿輸油管埋藏方向(Y軸)進(jìn)行積分，得到原油遷移污染土體的體積；接著對(duì)體積—時(shí)間進(jìn)行異速增長(zhǎng)曲線(allometric)擬合，得出原油在土體中的遷移體積與時(shí)間關(guān)系，即V-t關(guān)系.根據(jù)兩種土壤的V-t擬合公式與球體體積變形公式求出以漏油點(diǎn)為中心的遷移半徑r，遷移半徑r即為類似遷移距離L.該預(yù)測(cè)模型稱為球體積分模型.該模型的建立基于以下幾個(gè)假設(shè)：1) 原油以勻速連續(xù)泄漏；2)XOZ平面原油滲透斑跡為橢圓；3) 污染土體為球體.

1.2.3含水率對(duì)原油遷移的影響將采集的新鮮粉質(zhì)粘土經(jīng)過(guò)風(fēng)干或加水?dāng)嚢璧姆椒ㄖ苽湟唤M含水率不同的土樣，該組土樣含水率為：6.09%、8.55%、8.86%、10.33%、11.56%、12.3%、15.02%；將這一系列含水率的土壤分別進(jìn)行3 min的箱體動(dòng)態(tài)模擬實(shí)驗(yàn).

2結(jié)果與討論

2.1原油在土體各方向的遷移原油在土體中X軸的正負(fù)方向的遷移距離很接近，Y軸也如此，因此僅需X正軸、Y正軸、Z正軸和Z負(fù)軸方向的數(shù)據(jù)，分析原油在粉質(zhì)粘土和粘土混碎石土體中各方向遷移的時(shí)間距離關(guān)系.

由圖2、圖3可知，在粉質(zhì)粘土和粘土混碎石土體中，原油在各方向遷移距離的大小皆為：Z負(fù)軸>Y正軸=Y負(fù)軸>X正軸=X負(fù)軸>Z正軸.原油在土壤中遷移的作用力主要為重力和毛細(xì)管力[14]，在Z負(fù)軸，重力作為驅(qū)動(dòng)力促進(jìn)原油的移動(dòng)；在X軸、Y軸、Z正軸，重力阻礙原油的移動(dòng)，在Z正軸的阻礙作用遠(yuǎn)遠(yuǎn)大于X、Y軸.而原油在Y軸遷移的距離大于X軸，可能是因?yàn)?，Y軸方向?yàn)檩斢凸苈癫氐姆较颍寥琅c輸油管不能完全融合而留下一條縫隙，這條縫隙使原油更易在Y軸方向運(yùn)移.

圖2　原油在粉質(zhì)粘土中個(gè)各方位的遷移距

圖3　原油在粘土混碎石中各方位的遷移距

土壤粒徑越小，比表面積越大，吸附截留原油的性能越強(qiáng)[15]，原油在其中的遷移速率受到牽制；并且，粘土混碎石中含有粒徑大于10 mm的碎石顆粒，碎石顆粒幾乎沒(méi)有吸附性，能與周圍土壤顆粒間產(chǎn)生較大的縫隙，減小了對(duì)原油遷移的阻礙.有機(jī)質(zhì)雖然對(duì)土壤吸附原油有較大影響[16]，但兩種土壤的有機(jī)質(zhì)含量相差太小，可忽略其影響.因此，在該實(shí)驗(yàn)中，原油的遷移速率主要受土壤顆粒的影響，且原油在粉質(zhì)粘土中的遷移速率小于在粘土混碎石中的遷移速率.

2.2原油遷移距離預(yù)測(cè)結(jié)果根據(jù)1.2中的球體積分模型，將計(jì)算得到的土體被污染體積與泄漏時(shí)間(t)進(jìn)行異速增長(zhǎng)曲線(allometric)進(jìn)行擬合分析，得到被污染土體體積(V)與泄漏時(shí)間(t)之間的關(guān)系方程式.如下：

粉質(zhì)粘土：

V=17.93t2.06,R2=0.991 5

(1)

粘土混碎石：

V=0.16t2.02,R2=0.969 9

(2)

根據(jù)兩種土壤的V-t擬合公式與球體體積變形公式：

r=(3V/4π)1/3

(3)得到類似遷移距離L(即r)，進(jìn)一步圖像擬合，得到時(shí)間-遷移距離預(yù)測(cè)圖(見(jiàn)圖4).由圖4可明顯看出原油在粉質(zhì)粘土、粘土混碎石中遷移的距離隨時(shí)間的增加而增大，粘土混碎石的遷移速率大于粉質(zhì)粘土，且其遷移速率為粉質(zhì)粘土的兩倍多，這與2.1中的實(shí)驗(yàn)結(jié)果相同，說(shuō)明了球體積分模型的合理性.

圖4　原油在土壤中遷移距離預(yù)

圖5　含水率對(duì)遷移距離的影

3結(jié)論

1) 原油在粉質(zhì)粘土和粘土混碎石土體中各方向遷移距離的大小皆為：Z負(fù)軸>Y正軸=Y負(fù)軸>X正軸=X負(fù)軸>Z正軸.管道泄漏時(shí)，可以加強(qiáng)對(duì)Z負(fù)軸和Y軸方向的原油截留，防止原油遷移帶來(lái)更嚴(yán)重的污染.

2) 原油在粉質(zhì)粘土中的遷移距離隨著土壤含水率的增大先增大后減小，在含水率為10.33%時(shí)遷移距離達(dá)到最大.管道泄漏時(shí)，可以提高土壤含水率的方式減小原油的遷移速率，為清理回收原油爭(zhēng)取時(shí)間.

3) 箱體動(dòng)態(tài)模擬實(shí)驗(yàn)結(jié)果和球形積分模型預(yù)測(cè)結(jié)果一致，原油在粘土混碎石中的遷移速率大于在粉質(zhì)粘土中的遷移速率.

4參考文獻(xiàn)

[1] Akbari A,Ardestani M,Shayegan J.Distribution and mobility of petroleum hydrocarbons in soil:case study of the south pars gas complex,southern Iran[J].Iranian Journal of Science,2012,36(C2):265-275.

[2] Lu Zhe,Zeng Fangang,Xue Nandong,et al.Occurrence and distribution of polycyclic aromatic hydrocarbons in organo-mineral particles of alluvial sandy soil profiles at a petroleum-contaminated site[J].Science of the Total Environment,2012,433:50-57.

[3] Pan Feng,Ma Jinzhu,Wang Yunquan,et al.Simulation of the migration and transformation of petroleum pollutants in the soils of the Loess Plateau:a case study in the Maling oil field of northwestern China[J].Environmental Monitoring and Assessment,2013,185(10):8023-8034.

[4]李龍媛,羅澤嬌,彭輝.石油類污染物在黏性土壤中的垂直分布規(guī)律[J].安全與環(huán)境工程,2014,21(02):57-62.

[5] Ali Nuhad M,Saad Sarra A M,Elamin A.Mobility of total petroleum hydrocarbons in Shambat soil in Sudan[J].World Journal of Science,Technology and Sustainable Development,2014,11(2):134-143.

[7] Amro Mohammed Musa,Benzagouta Mohammed Said,Karnanda Wimpy.Investigation on crude oil penetration depth into soils[J].Arabian Journal of Geosciences,2013,6(3):873-880.

[8] Sun Xiaonan,Liu Anping,Chen Xiurong,et al.Research on the migration of petroleum hydrocarbon contamination in the soil in different leaching amount[J].Advanced Materials Research,2012,414:121-125.

[9] 潘峰,陳麗華,付素靜,等.石油類污染物在隴東黃土塬區(qū)土壤中遷移的模擬試驗(yàn)研究[J].環(huán)境科學(xué)學(xué)報(bào),2012,32(2):410-418.

[10] Yang M,Yang YS,Du X,et al.Fate and transport of petroleum hydrocarbons in Vadose Zone:compound-specific natural attenuation[J].Water Air Soil Pollut,2013,224(3):1-14.

[11] Sang Yuquan,Zheng Jingtang,Ma Xiaolei.Study on migration model of petroleum pollutant in soil[J].Advances in Information Sciences & Service Sciences,2013,5(2):661-669.

[12] 胡黎明,邢巍巍,周小文.非飽和土中多相流動(dòng)的試驗(yàn)研究和數(shù)值模擬[J].工程力學(xué)，2008,25(11)：162-166.

[13] Kererat Chusak,Sasanakul Inthuorn,Soralump Suttisak.Centrifuge modeling of LNAPL migration with a soil-cement barrier in case of flow and no-flow conditions[J].Bridges,2014,10:3513-3522.

[14 ] 李興柏.多年凍土區(qū)石油污染物遷移實(shí)驗(yàn)研究[M].甘肅:蘭州理工大學(xué),2012:38-46.

[15] 劉曉燕,李英麗,朱謙雅,等.石油類污染物在土壤中的吸附/解吸機(jī)理研究及展望[J].礦物巖石地球化學(xué)通報(bào),2007,26(1):82-87.

[16] 李永霞,鄭西來(lái),馬艷飛.石油污染物在土壤中的環(huán)境行為研究進(jìn)展[J].安全與環(huán)境工程.2011,18(4):43-47.

(責(zé)任編輯游俊)

The migration of crude oil in soil during the leakage of theChina-Russia oil pipeline

WANG Xueli1,2,WANG Rong3,WU Guansheng2,GAO Peng2,WANG Bo2,WANG Dongpo1

(1.School of Materials Science and Engineering,Tianjin University,Tianjin 300072,China;2.Daqing Oil and Gas Pipelines Branch Company of PetroChina Pipeline Company,Jiagedaqi 165000,China;3.College of Resources & Environment,Sichuan Agricultural University,Chengdu,611130,China)

Abstract:In order to explore the migration patterns of crude oil in soil during China-Russia oil pipeline leak,the method of box body dynamic simulation and integral sphere model had been applied with silty clay and clay mixed macadam as the experimental soil.The results showed that the migration distance of crude oil in each direction of soil decreased in the following order:vertically downward direction (negative Z axis) > direction of the pipeline (positive Y axis and negative Y axis) > horizontal vertical pipeline direction (positive X axis and negative X axis) > vertical upward direction (positive Z axis).The soil moisture content and oil pressure effect on crude oil migration distance were analyzed.As results showed,crude oil migration distance along with the increase of soil moisture content increased first,then decreased,in the case of soil moisture content within the scope of 6.09%15.02%.In addition,the moisture content of 10.33% migration distance was the largest.Meanwhile,the result predicted by spherical integral model that the migration rate of crude oil in soil decreased in the order:clay mixed macadam > silty clay.It was consistent with the result of box body dynamic simulation experimental,which proved the rationality of the spherical integral model.

Key words:pipeline leak; migration; crude oil; box body dynamic simulation

中圖分類號(hào):X53

文獻(xiàn)標(biāo)志碼:A

DOI:10.3969/j.issn.1000-2375.2016.03.013

文章編號(hào):1000-2375(2016)03-0242-04

作者簡(jiǎn)介：王學(xué)力(1974-)，男，高級(jí)工程師，E-mail：kjwxl@petrochina.com.cn

基金項(xiàng)目：校企合作項(xiàng)目(漠大線油品遷移規(guī)律及控制技術(shù)研究)資助

收稿日期:2015-08-18