黃 濤，程林松，曹仁義，王 鵬，賈志豪

(中國石油大學(xué)(北京) 石油工程學(xué)院，北京 102249)

引 言

頁巖油藏礦物組成復(fù)雜，常見的無機(jī)礦物成分有石英、長石、碳酸鹽礦物和黏土礦物等[1-3]。不同礦物的親油性不同，會對頁巖油的賦存和運移產(chǎn)生不同程度的影響，因此大量學(xué)者研究了不同礦物對頁巖油賦存運移的影響。寧方興等[4]運用掃描電鏡等手段分析了濟(jì)陽坳陷主力巖相石英、方解石和黏土礦物等礦物中頁巖油的賦存狀態(tài)；王民等[5]通過有機(jī)地球化學(xué)等方法研究了不同礦物組成下吸附油和游離油的含量；CUI等[6]基于高斯混合理論建立了頁巖油賦存狀態(tài)的模型，研究了有機(jī)質(zhì)和無機(jī)質(zhì)礦物對頁巖油吸附狀態(tài)的影響；ZHANG等[7]通過實驗研究了溫度和壓力對不同礦物表面潤濕性的影響。上述學(xué)者多基于實驗和理論推導(dǎo)礦物組成對頁巖油賦存狀態(tài)的影響，缺少了對微觀現(xiàn)象的描述和機(jī)理解釋。

分子動力學(xué)模擬可以從納米級別上模擬油氣的賦存及運移，研究頁巖油與儲層的微觀作用機(jī)制，因此也有學(xué)者運用此方法研究頁巖油與不同礦物表面的相互作用機(jī)理，ZHAO等[8-9]運用分子動力學(xué)的方法研究了不同礦化度下方解石不同表面對油分子運動的影響；WANG等[10]利用分子動力學(xué)模擬研究了石英納米孔中正辛烷的靜態(tài)性質(zhì)和壓力驅(qū)動下的流動行為；吳春正等[11]運用分子動力學(xué)模擬方法研究了混合烴在二氧化硅和石墨狹縫中的密度分布特征及溫度對有機(jī)質(zhì)狹縫內(nèi)烴分子吸附的影響；LUAN等[12]運用分子動力學(xué)模擬研究了石英礦物中二氧化碳驅(qū)油的機(jī)理；TIAN等[13]利用分子動力學(xué)模擬研究了多組分油混合物在高嶺石孔隙中的吸附特性。但這些學(xué)者多研究一到兩種礦物中原油的賦存運移行為，不能完整地體現(xiàn)頁巖儲層無機(jī)礦物中的頁巖油的賦存運移規(guī)律。

本文首先通過構(gòu)建石英、鈉長石、方解石和高嶺石4種礦物模型來分別表征頁巖儲層中的陸源碎屑巖、碳酸鹽巖及黏土礦物，然后通過構(gòu)建甲苯和己烷兩種組分代表頁巖油，最后運用分子動力學(xué)模擬了頁巖油在不同礦物表面的吸附擴(kuò)散，從微觀上描述了不同礦物對頁巖油賦存運移的規(guī)律并解釋了機(jī)理，對頁巖油藏的勘探開發(fā)有重要的意義。

1 模型的構(gòu)建

1.1 巖石和流體模型的構(gòu)建

根據(jù)文獻(xiàn)中的數(shù)據(jù)，高嶺石和鈉長石屬于三斜晶系，高嶺石單位晶胞的空間群為P1[14]，鈉長石單位晶胞的空間群為C1[15]；方解石和石英屬于三方晶系[16]，方解石單位晶胞的空間群為R3C，石英單位晶胞的空間群為P3221。礦物模型的具體晶胞參數(shù)見表1。根據(jù)表1中的數(shù)據(jù)，運用Material Studio 8.0軟件建立晶胞模型。4種礦物的晶胞如圖1所示。

表1 礦物模型的晶胞參數(shù)

圖1 4種礦物晶胞圖

礦物存在不同的晶面，不同晶面的穩(wěn)定性不同，并對孔隙內(nèi)流體有著不同的影響。迄今為止對方解石的{104}晶面[17 -18]和其他3種礦物的{001}晶面[19]開展的研究較多，所以本文研究基于方解石的{104}面和其他3種礦物的{001}晶面進(jìn)行。

頁巖油中的飽和烴和芳香烴的占比較高，一般大于80 %，因此建立以甲苯為代表的芳香烴模型和以己烷為代表的飽和烴模型，并且由于甲苯和己烷結(jié)構(gòu)的對稱性不同，也可以代表頁巖油中的非極性物質(zhì)與極性物質(zhì)[20]。甲苯和己烷的模型如圖2所示。

圖2 流體模型

1.2 計算模型的構(gòu)建

通過構(gòu)建平板狀的計算模型來模擬石油組分在礦物表面的潤濕，計算模型構(gòu)建如下：(1)首先為了消除模型邊界效應(yīng)，采用周期性邊界條件；(2)其次在空間群為P1下建立4a×2b×1c的礦物超晶胞模型；(3)最后將30個甲苯或者己烷分子組成的液滴模型放在不同礦物的超晶胞模型上，構(gòu)成計算模型。計算模型如圖3所示。

圖3 計算模型

2 模擬過程

本文中礦物模型和流體模型的力場均采用COMPASS力場，這是一個將從頭算法和實驗參數(shù)相結(jié)合的力場，能夠準(zhǔn)確地預(yù)測無機(jī)物和有機(jī)物的相關(guān)性質(zhì)[21]。其總勢能的函數(shù)公式為：

(1)

式中：前5項為成鍵形式的能量，分別表示鍵伸縮、鍵角變化、二面角變化、交叉項和離平面的彎曲勢變化引起的能量改變的加和，b，θ，φ，χ分別為鍵長、鍵角、二面角和離平面振動的角度；最后2項表示非鍵相互作用，分別為靜電相互作用和范德華力相互作用，其中采用L-J9-6勢能模型計算范德華力和靜電力的相互作用[22]，其計算方法為：

(2)

(3)

計算模型采用Smart算法計算15 000步找到最優(yōu)結(jié)構(gòu)。靜電作用勢能計算選擇Ewald求和法，分子之間的范德華力計算采用Atom based算法，其截斷半徑為1.55 nm。

在動力學(xué)計算中，采用NVT系綜，控溫方法采用Nose-Hoover擴(kuò)展系統(tǒng)法，總模擬時長500 ps，步長為1 fs。其他參數(shù)與前文一致。

3 結(jié)果與討論

利用軟件模擬分析在溫度為298 K條件下甲苯和己烷在4種礦物表面運移擴(kuò)散的規(guī)律。

3.1 潤濕接觸角

結(jié)構(gòu)優(yōu)化后的甲苯和己烷在4種礦物表面的潤濕角如圖4和5所示，從圖中可以發(fā)現(xiàn)，甲苯在除高嶺石之外的3種礦物表面也會向球形轉(zhuǎn)變，甲苯在高嶺石表面會平鋪開來。通過擬合甲苯在4種礦物表面的接觸角可得石英的接觸角約為92.7°，方解石的接觸角約為80.7°，鈉長石的接觸角約為68.5°。己烷在4種礦物的表面都會聚集起來向球形轉(zhuǎn)變。通過擬合己烷在4種礦物表面的接觸角可得石英的接觸角約為89.7°，方解石的接觸角約為78.5°，鈉長石的接觸角約為63.5°，高嶺石的接觸角約為61.3°。接觸角越小，礦物對該液體的潤濕性越強，親油性越好，因此4種礦物的親油性排序由高到低為：高嶺石>鈉長石>方解石>石英。這與張永超[23]和吳春正[24]的結(jié)果基本類似，證明了模型的正確性。

圖4 甲苯在4種礦物表面的潤濕接觸角

圖5 己烷在4種礦物表面的潤濕接觸角

3.2 相對密度分布

分子動力學(xué)平衡過后，甲苯和己烷在4種礦物表面的分布如圖6所示。從圖6中可知，甲苯和己烷會在礦物表面重新分布，但不同礦物表面的甲苯和己烷分子的密度不同，其中甲苯和己烷都會在高嶺石、鈉長石和方解石礦物表面形成一個明顯的聚集層，但甲苯則在遠(yuǎn)離石英礦物表面的真空層形成聚集，己烷會在石英礦物的表面和遠(yuǎn)離礦物表面的真空層中各形成一個吸附層。由前文的分析可知，這是由于不同礦物表面的親油性不同造成的，石英的親油性最差，因此礦物表面的油分子吸附得少，多余的油分子會運動到遠(yuǎn)離礦物的真空層中；高嶺石親油性最強，因此表面的油分子密度最大。

圖6 甲苯和己烷在4種礦物表面的相對密度分布

3.3 速度剖面

甲苯和己烷在4種礦物表面的速度剖面圖分別如圖7和圖8所示，其中速度分為沿X軸、Y軸和Z軸3個方向的分速度。從圖中可知，甲苯在鈉長石和方解石表面的速度峰值出現(xiàn)在遠(yuǎn)離礦物表面的真空層中，石英和高嶺石表面的甲苯運動速度較快，而在真空層中的甲苯運動的速度相對較小，尤其是石英真空層中甲苯運動速度幾乎為0；己烷在石英礦物的表面運動速度很快，在遠(yuǎn)離礦物表面的真空層中運動速度只有礦物表面速度的一半左右，而己烷在另外3種礦物表面的表現(xiàn)正好相反，其速度的峰值出現(xiàn)在真空層中。結(jié)合前文分析的相對密度分布可知，這是由于分子動力學(xué)平衡后甲苯和己烷在高嶺石、鈉長石和方解石表面形成聚集層，所以在真空層中的甲苯和己烷會加速運動到礦物表面，因此鈉長石和方解石真空層中的甲苯以及高嶺石、鈉長石和方解石真空層中的己烷的運動速度較大。甲苯在高嶺石礦物表面結(jié)構(gòu)優(yōu)化后就平鋪在礦物表面，并且表面的甲苯分子較多，所以會在高嶺石表面出現(xiàn)速度剖面的峰值；甲苯和己烷在石英礦物表面密度較小，真空層中密度較高，因此表面的甲苯和己烷加速運動到真空層中，所以石英礦物表面的甲苯和己烷速度較大。

圖7 甲苯在4種礦物表面的速度剖面圖

圖8 己烷在4種礦物表面的速度剖面圖

3.4 相互作用機(jī)制

模擬計算的4種礦物與甲苯和己烷之間的相互作用能如圖9所示，其中靜電能是指由于分子中的原子帶有部分電荷而引起靜電吸引或排斥作用，主要描述粒子之間的庫侖力。圖中能量正負(fù)值代表吸附和排斥作用，正值代表排斥，負(fù)值代表吸附。從圖中可知，甲苯和己烷與4種礦物之間的相互作用總能量均為負(fù)值，說明甲苯和己烷與4種礦物之間為吸附關(guān)系。甲苯和己烷與4種礦物的相互作用總能量按從大到小排序為高嶺石>鈉長石>方解石>石英，這與前文的4種礦物的親油性的關(guān)系一致，表明礦物表面與液體分子之間的作用力的大小決定了礦物的親油性，礦物表面與油分子的相互作用力越強，礦物越親油。4種礦物與甲苯之間的靜電能大于范德華能， 表明4種礦物與甲苯之間吸附時靜電力起主要作用。4種礦物吸附己烷時，范德華能的絕對值均大于靜電能的絕對值，說明4種礦物吸附己烷時范德華能起主要作用；除了高嶺石外其他3種礦物與己烷之間的靜電能均為正值，說明高嶺石與己烷相互作用時靜電能起吸附作用，其他3種礦物與己烷之間的靜電能起排斥作用。甲苯與礦物表面相互作用能均大于己烷與同種礦物表面的相互作用能，說明在相同條件下，礦物表面對甲苯的吸附作用更強。

圖9 甲苯和己烷與4種礦物的相互作用能

3.5 自擴(kuò)散系數(shù)

油分子在4種礦物表面自擴(kuò)散系數(shù)見表2。甲苯在4種礦物表面的自擴(kuò)散系數(shù)均小于己烷，結(jié)合前文的相互作用能分析認(rèn)為，甲苯作為一種極性分子與礦物表面的相互作用力大于己烷這種非極性分子，導(dǎo)致甲苯在礦物表面的吸附作用更強，更加不容易擴(kuò)散。甲苯和己烷在4種礦物表面的自擴(kuò)散系數(shù)的大小關(guān)系均為石英>方解石>鈉長石>高嶺石，與前文的相互作用能大小成反比，說明了油分子與礦物表面的相互作用力越強，油分子在礦物表面越不容易擴(kuò)散。

表2 甲苯和己烷在4種礦物中的自擴(kuò)散系數(shù)

3.6 徑向分布函數(shù)

徑向分布函數(shù)可以解釋為以空間中一個粒子為中心，去尋找周圍粒子的概率，即粒子在周期性邊界條件下的區(qū)域密度和全局密度的比值[25]。分子動力學(xué)計算徑向分布函數(shù)的方法為：

(4)

式中：g(r)為徑向分布函數(shù)；ρ為空間的密度；N為粒子的數(shù)目；T為計算的時間(步數(shù))；δr為設(shè)定的距離差；ΔN為介于r→r+δr間的粒子數(shù)目。

徑向分布函數(shù)的值為0時對應(yīng)的最大橫坐標(biāo)值為粒子之間的最小間距，徑向分布函數(shù)的峰值對應(yīng)的橫坐標(biāo)值為大部分粒子之間的間距，峰值越大，表示該粒子間距對應(yīng)的區(qū)域密度越大，粒子的聚集程度越高。

分子動力學(xué)平衡之后，4種礦物表面甲苯和己烷分子間的徑向分布函數(shù)如圖10所示。從圖中可知，當(dāng)分子之間的距離小于0.375 nm時，甲苯和己烷的徑向分布函數(shù)都為0，這表示油分子在4種礦物表面的最近的間距不會小于0.375 nm。甲苯在鈉長石和高嶺石表面的徑向分布函數(shù)都存在一個峰值，鈉長石峰值對應(yīng)的距離為0.725～0.825 nm，高嶺石峰值對應(yīng)的距離為0.675～0.775 nm，且鈉長石的峰值大于高嶺石的峰值，說明甲苯在鈉長石表面分布得更集中；甲苯在方解石和石英表面的徑向分布函數(shù)峰值相對不明顯，表示甲苯在這2種礦物表面的分布比較分散。己烷在4種礦物表面的徑向分布函數(shù)的曲線具有一致性，曲線比較平滑，沒有突出的峰值，說明己烷在4種礦物表面的排列比較分散。結(jié)合前文的分析可知，礦物的親油性越強，油分子的極性越高，油分子的排列越緊密。

圖10 甲苯和己烷在4種礦物表面的徑向分布函數(shù)

4 結(jié) 論

(1)通過分析甲苯和己烷在4種礦物的接觸潤濕角得到了4種礦物的親油性大小關(guān)系為高嶺石>鈉長石>方解石>石英。

(2)礦物的親油性越強，表面油分子的密度越大，排列越緊密，越不容易擴(kuò)散，并且真空層中的油分子會加速運動到表面。

(3)礦物表面與油分子的相互作用力越強，礦物越親油。4種礦物吸附甲苯時靜電力和范德華力都起吸附作用，但以靜電力為主。4種礦物吸附己烷時以范德華力為主，靜電力在鈉長石、方解石和石英中起排斥作用，只有在高嶺石中起吸附作用。

(4)以甲苯為代表的頁巖油極性物質(zhì)在同種礦物表面的潤濕接觸角更小，吸附密度更高，排列更緊密，相互作用力更強，吸附更穩(wěn)定，說明礦物對頁巖油中極性物質(zhì)的吸附作用強于非極性物質(zhì)。