季佳圓，趙伶玲，李偲宇

（東南大學(xué)能源熱轉(zhuǎn)換及其過程測控教育部重點(diǎn)實驗室，能源與環(huán)境學(xué)院，江蘇 南京 210096）

溫度及壓強(qiáng)對CO2-鹽水系統(tǒng)界面張力的影響

季佳圓，趙伶玲，李偲宇

（東南大學(xué)能源熱轉(zhuǎn)換及其過程測控教育部重點(diǎn)實驗室，能源與環(huán)境學(xué)院，江蘇 南京 210096）

在超臨界態(tài)CO2封存于深部鹽水層過程中，溫度、壓強(qiáng)等控制條件是影響封存效率和封存量的重要因素。應(yīng)用分子動力學(xué)模擬的方法對343～373 K和6～35 MPa范圍內(nèi)的CO2-NaCl鹽水系統(tǒng)進(jìn)行了界面張力（IFT）及界面特性的研究，分析了IFT隨溫度及壓強(qiáng)的變化關(guān)系，并觀測到了壓力平衡點(diǎn)pplateau；從分子尺度（物質(zhì)密度、界面過余量、界面水合物密度）分析了IFT隨壓強(qiáng)、溫度的變化，以及pplateau產(chǎn)生的原因。結(jié)果表明，pplateau前壓強(qiáng)升高或溫度降低將導(dǎo)致CO2密度升高，IFT下降，而pplateau后IFT趨于穩(wěn)定且受溫度影響較?。籆O2的界面過余量及界面處水合物數(shù)量隨壓強(qiáng)及溫度變化，與IFT的變化相反；高壓下界面水合物密度的飽和現(xiàn)象可能是pplateau產(chǎn)生的重要原因。

CO2-鹽水系統(tǒng)；界面張力；分子模擬；溫度；壓強(qiáng)；水合物

引 言

深部鹽水層CO2地質(zhì)封存所需的注射能耗及最大地質(zhì)埋存深度與 CO2-鹽水之間的界面張力（interfacial tension，IFT）直接相關(guān)[1-4]，并受溫度和壓強(qiáng)的制約。開展溫度、壓強(qiáng)對CO2-鹽水間界面張力的影響研究，不僅可以分析IFT隨溫度、壓強(qiáng)等控制參數(shù)的變化規(guī)律，還能闡述溫度、壓強(qiáng)對IFT產(chǎn)生影響的內(nèi)在機(jī)理進(jìn)而對指導(dǎo)不同環(huán)境條件（溫度、壓強(qiáng)）下的CO2地質(zhì)封存設(shè)計，提高注射安全性及存儲容量具有重大意義。

目前實驗[1,5-8]已測定相關(guān)儲層條件下 CO2-水和CO2-鹽水系統(tǒng)的IFT值，并觀測到IFT在定溫條件下會隨著壓強(qiáng)升高而降低，并在壓力平衡點(diǎn)pplateau之后趨于穩(wěn)定值。實驗還發(fā)現(xiàn)pplateau的大小與鹽的種類及鹽度無明顯聯(lián)系，僅隨溫度升高而上升。Chalbaud等[6]將pplateau的存在歸因于CO2溶解度的影響，但尚未展開深入分析。

分子動力學(xué)模擬（molecular dynamics simulation，MD模擬）可以研究多相界面系統(tǒng)的微觀特性，目前該方法已成功模擬了 CO2-水及 CO2-鹽水系統(tǒng)，可獲得與實驗一致的IFT值[9-13]，并能觀測界面的微觀現(xiàn)象，是一種有效的研究手段。

本文應(yīng)用 MD模擬方法，對 343～373 K和6～35 MPa范圍內(nèi)的CO2-NaCl系統(tǒng)進(jìn)行計算，分析了體相及界面各物質(zhì)性質(zhì)隨環(huán)境條件的變化規(guī)律，包括CO2的密度、CO2的界面過余量、界面處CO2水合物數(shù)量等，探討了IFT對溫度及壓強(qiáng)依賴關(guān)系的物理機(jī)理，尤其對pplateau現(xiàn)象的產(chǎn)生原因進(jìn)行了分析，可為IFT的控制和預(yù)測提供理論依據(jù)。

1 研究對象及方法

1.1 對象及模型

本文根據(jù)Chalbaud等[6]對CO2-NaCl系統(tǒng)大范圍溫度及壓強(qiáng)下的IFT實驗研究結(jié)果，選擇了溫度及壓強(qiáng)范圍為 343～373 K和 6～35 MPa的CO2-NaCl系統(tǒng)為研究對象。具體工況參數(shù)列于表1。

表1 CO2-NaCl系統(tǒng)的溫度和壓強(qiáng)條件Table 1 Temperature and pressure conditions of CO2-NaCl system

在計算過程中，綜合考慮了系統(tǒng)內(nèi)分子間非鍵結(jié)作用力（范德華力、庫倫靜電力）及分子內(nèi)鍵結(jié)作用力（鍵拉伸和鍵彎曲）。分別采用Lennard-Jones勢能函數(shù)[14]模擬范德華力，庫侖定律模擬庫侖靜電力，具體分子間勢能函數(shù)如式（1）所示

其中，rij為原子i與j之間的距離；εij為勢能阱的深度，εij(εiiεjj)1/2；σij為兩體互相作用的勢能為零時的距離，σij(σiiσjj)1/2；ε0為真空介電常數(shù)；qi及 qj為原子 i與j所帶電荷量。其中采用PME技術(shù)[15]模擬分子間長程庫侖作用力，范德華作用截距設(shè)定為0.9 nm[16-18]。此外，本文采用諧波勢能函數(shù)[5]模擬鍵拉伸和鍵角彎曲等分子內(nèi)鍵結(jié)作用力。

本文計算中，水分子選擇柔性F3C模型[19]，CO2選擇柔性EPM2模型[20]，鹽離子采用Chandrasekhar等[21]開發(fā)的模型。應(yīng)用MD軟件Gromacs4.5[22]并采用周期性邊界條件進(jìn)行計算，所建立的橫截面4 nm×4 nm的計算域示于圖1。計算域中間區(qū)域為鹽水，包括4323個水分子、147個Na+和147個Cl-，對應(yīng)鹽度為1.89 mol·L-1，兩側(cè)分別為 732個CO2分子。利用 Berendsen方法[23-24]來實現(xiàn)溫度和壓強(qiáng)的設(shè)定。由于系統(tǒng)在NPzT系綜下20 ns達(dá)平衡態(tài)，故模擬時間運(yùn)行30 ns，選取最后5 ns為有效數(shù)據(jù)進(jìn)行分析。

圖1 CO2-NaCl系統(tǒng)平衡狀態(tài)Fig.1 CO2-NaCl solution system in equilibrium state

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

計算系統(tǒng)中包括兩個界面，根據(jù)界面張力γ定義[25]，系統(tǒng)中γ可由式（2）計算

式中，pxx, pyy, pzz分別為沿x，y，z方向壓強(qiáng)張量對z向的角分量。

吉布斯分界面（Gibbs dividing surface，GDS）為以體相為參照，垂直于界面方向上CO2相多余水分子與溶液相缺乏水分子相等的位置，其厚度采用0.1～0.9倍水密度之間的距離表示[26]。

界面過余量[27]表征物質(zhì)在界面和體相中的量差異，本文中CO2-NaCl系統(tǒng)界面處i分量的界面總過余量可用式（3）表示

2 結(jié)果與討論

2.1 IFT值的變化

本文分別模擬了343 K和373 K的各壓強(qiáng)下鹽濃度為1.89 mol·L-1的CO2-NaCl系統(tǒng)，IFT的計算結(jié)果示于圖2。由圖2可以看出，模擬所得的 IFT值與Chalbaud等[6]的實驗數(shù)據(jù)吻合較好。在溫度為343 K時[圖2(a)]，IFT值隨著壓強(qiáng)的不斷升高，其減幅逐漸下降，最終在壓力平衡點(diǎn)pplateau=15 MPa之后達(dá)到穩(wěn)定值34 mN·m-1。373 K時[圖2(b)]的IFT值變化與343 K時極為相似，IFT值同樣隨著壓強(qiáng)的升高而下降，且下降速度不斷趨于平緩，由于高壓情況下實驗設(shè)備等因素受限制，實驗方法尚未得到25 MPa以上的數(shù)據(jù)點(diǎn)，故在實驗壓強(qiáng)范圍內(nèi)并未觀測到壓力平衡點(diǎn)。而本文應(yīng)用分子模擬的方法[28]，成功模擬了25～35 MPa下的CO2-NaCl系統(tǒng)，觀測到壓強(qiáng)在pplateau=25 MPa之后達(dá)到穩(wěn)定值33 mN·m-1。

圖2 343 K和373 K下CO2-NaCl系統(tǒng)的IFT值Fig.2 IFT values of CO2-NaCl system at 343 K and 373 K

對比圖2(a)與圖 2(b)還可以發(fā)現(xiàn)，當(dāng)鹽濃度一定（1.89 mol·L-1）時，pplateau之前的 CO2-NaCl系統(tǒng)任意壓強(qiáng)點(diǎn)的IFT值在溫度為343 K時均小于373 K；同時隨著壓強(qiáng)的不斷升高，兩者差值不斷縮小；而至壓力平衡點(diǎn)之后兩溫度下的 IFT穩(wěn)定值相接近，故溫度對高壓系統(tǒng)IFT穩(wěn)定值的影響不大。此外，由于343 K下壓力平衡點(diǎn)為15 MPa，373 K壓力平衡點(diǎn)為25 MPa，故可推測壓力平衡點(diǎn)會隨溫度升高而增大。

2.2 密度分布圖

圖3 343 K和373 K下界面CO2、H2O、鹽離子密度Fig.3 Density profiles of CO2, water and salts at 343 K and 373 K

為從分子角度更加詳細(xì)地描述界面性質(zhì)的變化，本文選取并計算了界面處 2 nm厚度區(qū)域的CO2、H2O、Na+及Cl-的密度，結(jié)果示于圖3。由圖3可以觀測到：當(dāng)從 CO2相向溶液相過渡時，CO2密度逐漸下降，H2O密度則逐漸上升；而鹽離子密度則從界面中心向溶液相，由0逐漸達(dá)到某一穩(wěn)定值。由圖3可以看出，溫度一定時，CO2在CO2相的密度隨著壓強(qiáng)的升高逐漸升高。由[圖 3(a)]可以看出，溫度為343 K時，當(dāng)壓強(qiáng)由6.5 MPa升高至22 MPa，CO2相的CO2密度從580 kg·m-3升高至800 kg·m-3左右；而溫度為 373 K 時[圖 3(b)]，當(dāng)壓強(qiáng)由11 MPa升高至17 MPa，CO2相的CO2密度從約420 kg·m-3升高至 600 kg·m-3左右。此外，對比圖3(a)和(b)還可知，當(dāng)壓強(qiáng)恒定在11 MPa時，溫度由343 K升高到373 K，CO2相的CO2密度則由700 kg·m-3降至 420 kg·m-3左右。

綜上所述，壓強(qiáng)的升高和溫度的降低均對界面處CO2相的CO2密度產(chǎn)生較大的影響，而CO2密度的變化與IFT值的變化聯(lián)系密切。具體而言，壓強(qiáng)升高或溫度降低時，CO2分子排列緊湊，界面分子受力的不均勻性有所好轉(zhuǎn)，增大了其與水分子之間的吸引力。高密度下CO2分子更易進(jìn)入水相，從而破壞了部分水分子間的氫鍵作用，界面上水分子受到內(nèi)部分子的引力有所降低。以上兩方面因素導(dǎo)致了界面水分子所受內(nèi)部吸引力的減小，進(jìn)而使得IFT值下降。

2.3 CO2界面過余量

由上述有關(guān)界面處密度的討論可知，溫度及壓強(qiáng)對CO2相的CO2密度有較大的影響，為了更加全面地分析界面處CO2分子對界面張力的影響，本文進(jìn)一步分析了 343～373 K和 6～35 MPa的CO2-NaCl系統(tǒng)的 CO2界面過余量，計算結(jié)果示于圖4。由圖4可知，343 K和373 K時，CO2界面過余量均隨著壓強(qiáng)的升高而逐漸升高，直至到達(dá)穩(wěn)定值。具體來說，在343 K時，CO2界面過余量隨著壓強(qiáng)的升高不斷升高，從壓強(qiáng) 6.5 MPa 的 1.24 μmol·m-2，直至壓力平衡點(diǎn) 15 MPa之后達(dá)到并穩(wěn)定在 1.50 μmol·m-2左右；373 K 時 CO2界面過余量變化規(guī)律與343 K時基本一致，界面過余量從壓強(qiáng)8 MPa 時的 0.93 μmol·m-2開始逐漸升高，直至壓力平衡點(diǎn) 25 MPa 后穩(wěn)定在約 1.57 μmol·m-2。

圖4 343 K和373 K下的界面張力和CO2界面過余量Fig.4 IFT values and surface excess of CO2at 343 K and 373 K

由圖4中CO2界面過余量及IFT和溫度、壓強(qiáng)的變化關(guān)系可以發(fā)現(xiàn)，在 pplateau之前，CO2界面過余量隨著壓強(qiáng)的升高或溫度的降低而升高，而pplateau后，CO2界面過余量趨于穩(wěn)定且受溫度影響較小。CO2界面過余量受溫度、壓強(qiáng)的影響與IFT受溫度、壓強(qiáng)的影響結(jié)果恰好相反。這表明IFT受溫度、壓強(qiáng)變化的影響可理解為：壓力平衡點(diǎn)之前的壓強(qiáng)的升高或溫度的降低均使得CO2分子在界面處的累積量增加，界面處CO2數(shù)量的增加可能是導(dǎo)致IFT降低的原因之一。

2.4 CO2水合物密度

由于CO2界面過余量隨溫度、壓強(qiáng)變化的關(guān)系與IFT值變化規(guī)律相反，界面處CO2分子對IFT有著較大影響?？紤]到界面附近CO2分子與水分子作用距離較近，極易形成CO2水合物[29-30]，故本文進(jìn)一步研究了343～373 K和6～35 MPa的CO2-NaCl系統(tǒng)界面處的水合物數(shù)量及密度。其中部分溫度、壓強(qiáng)條件下的界面處水合物分析結(jié)果示于圖5。

圖5 343 K和373 K下界面處水合物數(shù)量Fig.5 Hydrates quantities at 343 K and 373 K

由圖5可知，當(dāng)溫度恒定時，界面處水合物數(shù)量隨著壓強(qiáng)的升高而增加，并逐漸趨于穩(wěn)定值。當(dāng)溫度為343 K時，水合物數(shù)量由90不斷升高并穩(wěn)定至118（6.5～24 MPa）；而373 K時，水合物數(shù)量則由74升高并穩(wěn)定至130（11～35 MPa）。

為排除模型尺寸對水合物數(shù)量的影響，本文將各溫度、壓強(qiáng)條件下的水合物數(shù)量除以截面積，計算出水合物面密度，其值示于圖6。

由圖6可知，水合物密度隨壓強(qiáng)升高而逐漸降低，下降速率逐漸趨于平緩，最終在IFT的壓力平衡點(diǎn)pplateau處達(dá)到飽和，密度值恒定。水合物密度與IFT值隨溫度、壓強(qiáng)變化呈現(xiàn)相反趨勢。因此可推測，溫度及壓強(qiáng)直接影響了界面處水合物密度，而界面處水合物密度的變化可能是導(dǎo)致界面張力變化的重要因素。

圖6 343 K和373 K下的界面張力和界面處水合物密度Fig.6 IFI values and hydrates densities at 343 K and 373 K

3 結(jié) 論

本文應(yīng)用MD模擬的方法進(jìn)行了溫度和壓強(qiáng)在343～373 K 和 6～35 MPa范圍內(nèi)鹽濃度為 1.89 mol·L-1的CO2-NaCl鹽水系統(tǒng)的界面特性研究，探討了IFT隨溫度、壓強(qiáng)變化的微觀機(jī)理及壓力平衡點(diǎn)pplateau的存在原因，得出以下結(jié)論。

（1）CO2密度隨著壓強(qiáng)增加或溫度降低而升高，其密度升高將增大CO2分子與界面水分子的引力，降低界面水分子所受到的水相內(nèi)部的引力，最終導(dǎo)致壓強(qiáng)增加或溫度降低時IFT下降的趨勢。

（2）恒溫時 CO2界面過余量隨壓強(qiáng)升高而增加，并且在pplateau之后增加變緩。該值隨溫度、壓強(qiáng)變化的趨勢與IFT的變化趨勢相反，CO2界面過余量對于IFT存在顯著的負(fù)作用。

（3）界面處 CO2水合物的密度隨溫度、壓強(qiáng)變化的趨勢與IFT的變化趨勢相反。溫度恒定時，界面處水合物數(shù)量隨著壓強(qiáng)的升高而升高，升高速率逐漸降低，隨后趨于定值。界面處水合物數(shù)量在高壓下的飽和現(xiàn)象可能是導(dǎo)致壓力平衡點(diǎn)pplateau產(chǎn)生的主要因素。

符 號 說 明

A ——界面截面積，m2

p ——壓強(qiáng)，MPa

pplateau——壓力平衡點(diǎn)，MPa

pxx——x方向壓強(qiáng)張量對z向的角分量，Pa

pyy——y方向壓強(qiáng)張量對z向的角分量，Pa

pzz——z方向壓強(qiáng)張量對z向的角分量，Pa

qi——i離子的電荷數(shù)，e

qj——j離子的電荷數(shù)，e

rij——原子i和j之間的距離，nm

Si——i物質(zhì)的界面過余量，μmol·m-2

T——溫度，K

γ——界面張力，mN·m-1

εij——蘭納-瓊斯勢能阱深度，kJ·mol-1

ε0——真空介電常數(shù)

σij——蘭納-瓊斯勢能為零時距離，nm

[1] CHIQUET P, DARIDON J L, BROSETA D, et al. CO2/Water interfacial tensions under pressure and temperature conditions of CO2geological storage [J]. Energy Convers. Manage., 2007, 48 (3)：736-744.

[2] JR F M O. Onshore geologic storage of CO2[J]. Science, 2009, 325(5948)： 1656-1658.

[3] LI Z, DONG M, LI S, et al. CO2sequestration in depleted oil and gas reservoirs—caprock characterization and storage capacity [J]. Energy Convers. Manage., 2006, 47 (11/12)： 1372-1382.

[4] JUNG J W, WAN J. Supercritical CO2and ionic strength effects on wettability of silica surfaces： equilibrium contact angle measurements[J]. Energy Fuels, 2012, 26 (9)： 6053-6059.

[5] AGGELOPOULOS C A, ROBIN M, VIZIKA O. Interfacial tension between CO2and brine (NaCl+CaCl2) at elevated pressures and temperatures： the additive effect of different salts [J]. Adv. Water Resour., 2011, 34 (4)： 505-511.

[6] CHALBAUD C, ROBIN M, LOMBARD J M, et al. Interfacial tension measurements and wettability evaluation for geological CO2storage [J]. Adv. Water Resour., 2009, 32 (1)： 98-109.

[7] GEORGIADIS A, MAITLAND G, TRUSLER J P M, et al. Interfacial tension measurements of the (H2O+CO2) system at elevated pressures and temperatures [J]. J. Chem. Eng. Data, 2010, 55 (10)： 4168-4175.

[8] HEBACH A, OBERHOF A, DAHMEN N, et al. Interfacial tension at elevated pressures-measurements and correlations in the water+carbon dioxide system [J]. J. Chem. Eng. Data, 2002, 47 (6)： 1540-1546.

[9] ZHAO L L, LIN S C, MENDENHALL J D, et al. Molecular dynamics investigation of the various atomic force contributions to the interfacial tension at the supercritical CO2-water interface [J]. J.Phys. Chem. B, 2011, 115 (19)：6076-6087.

[10] ZHAO L L, JI J Y, TAO L, et al. Ionic effects on supercritical CO2-brine interfacial tensions： molecular dynamics simulations and a universal correlation with ionic strength, temperature, and pressure [J].Langmuir, 2016, 32 (36)： 9188-9196.

[11] DA-ROCHA S R P, JOHNSTON K P, WESTACOTT R E, et al.Molecular structure of the water-supercritical CO2interface [J]. J.Phys. Chem. B, 2001, 105 (48)： 12092-12104.

[12] MüLLER E A, MEJíA A. Resolving discrepancies in the measurements of the interfacial tension for the CO2+H2O mixture by computer simulation [J]. J. Phys. Chem. Lett., 2014, 5 (7)： 1267-1271.

[13] NIELSEN L C, BOURG I C, SPOSITO G. Predicting CO2-waterinterfacial tension under pressure and temperature conditions of geologic CO2storage [J]. Geochim. Cosmic. Acta, 2012, 81： 28-38.

[14] 范康年. 物理化學(xué) [M]. 2版. 北京：高等教育社, 2005： 222.FAN K N. Physical Chemistry [M]. 2nd ed. Beijing： Higher Education Press, 2005： 222.

[15] AGGELOPOULOS C A, ROBIN M, PERFETTI E, et al. CO2/CaCl2solution interfacial tensions under CO2geological storage conditions：influence of cation valence on interfacial tension [J]. Adv. Water Resour., 2010, 33 (6)： 691-697.

[16] KLAUDA J B, WU X, PASTOR R W, et al. Long-range Lennard-Jones and electrostatic interactions in interfaces： application of the isotropic periodic sum method [J]. J. Phys. Chem. B, 2007, 111(17)： 4393-4400.

[17] DE LARA L S, MICHELON M F, MIRANDA C R. Molecular dynamics studies of fluid/oil interfaces for improved oil recovery processes [J]. J. Phys. Chem. B, 2012, 116 (50)： 14667-14676.

[18] PATEL S A, BROOKS C L. Revisiting the hexane-water interface via molecular dynamics simulations using nonadditive alkane-water potentials [J]. J. Chem. Phys., 2006, 124 (20)： 137-159.

[19] LEVITT M, HIRSHBERG M, SHARON R, et al. Calibration and testing of a water model for simulation of the molecular dynamics of proteins and nucleic acids in solution [J]. J. Phys. Chem. B, 1997, 101(25)： 5051-5061.

[20] NIETO-DRAGHI C, DE BRUIN T, P REZ-PELLITERO J, et al.Thermodynamic and transport properties of carbon dioxide from molecular simulation [J]. J. Chem. Phys., 2007, 126 (6)： 064509.

[21] CHANDRASEKHAR J, SPELLMEYER D C, JORGENSEN W L.Energy component analysis for dilute aqueous solutions of lithium(1+), sodium (1+), fluoride (1-), and chloride (1-) ions [J]. J. Am. Chem.Soc., 1984, 106 (4)： 903-910.

[22] PRONK S, PALL S, SCHULZ R, et al. GROMACS 4.5： a high-throughput and highly parallel open source molecular simulation toolkit [J]. Bioinformatics, 2013, 29 (7)： 845-854.

[23] BERENDSEN H J, POSTMA J P M, VAN GUNSTEREN W F, et al.Molecular dynamics with coupling to an external bath [J]. J. Chem.Phys., 1984, 81 (8)： 3684-3690.

[24] BUSSI G, DONADIO D, PARRINELLO M. Canonical sampling through velocity rescaling [J]. J. Chem. Phys., 2007, 126 (1)： 014101.

[25] ALEJANDRE J, TILDESLEY D J, CHAPELA G A. Molecular dynamics simulation of the orthobaric densities and surface tension of water [J]. J. Chem. Phys., 1995, 102 (11)： 4574-4583.

[26] YUET P K, BLANKSCHTEIN D. Molecular dynamics simulation study of water surfaces： comparison of flexible water models [J]. J.Phys. Chem. B, 2010, 114 (43)： 13786-13795.

[27] LI X, ROSS D A, TRUSLER J P, et al. Molecular dynamics simulations of CO2and brine interfacial tension at high temperatures and pressures [J]. J. Chem. Eng. Data, 2013, 117 (18)： 5647-5652.

[28] LEACH A R. Molecular Modelling： Principles and Applications [M].2nd ed. England, Harlow： Prentice Hall, 2001.

[29] SADLEJ J, MAKAREWICZ J, CHA?ASI?SKI G. Ab initio study of energy, structure and dynamics of the water-carbon dioxide complex[J]. J. Chem. Phys., 1998, 109 (10)： 3919.

[30] TEWES F, BOURY F. Thermodynamic and dynamic interfacial properties of binary carbon dioxide-water systems [J]. J. Phys. Chem.B, 2004, 108 (7)： 2405-2412.

Temperature and pressure effect on interfacial tensions of CO2-brine system

JI Jiayuan, ZHAO Lingling, LI Siyu
(Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy & Environment,Southeast University, Nanjing 210096, Jiangsu, China)

Temperature and pressure are important effect factors of the efficiency and quantity of CO2storage in the deep saline aquifers. Molecular dynamics (MD) simulation is applied to investigate the CO2-NaCl systems in 343—373 K and 6—35 MPa, the interfacial tensions (IFT) which are obtained from the simulations are consistent with experimental results and the pressure balance point pplateauare observed as well in this paper. Meanwhile, the interfacial tension variations with the temperature and pressure are analyzed and the reasons of pplateaufrom the molecular viewpoint are explained. The results show that the pressure rise and temperature decline will increase CO2density and decrease IFT before pplateau, but after pplateauthe IFT will be stable and less affected by temperature.In addition, the changes of CO2surface excess and the hydrate quantities with temperature and pressure showed the opposite trend compared with IFT variations, the saturation phenomena of hydrates at the interface under high pressure may be the fundamental reason of pplateau.

CO2-brine system; interfacial tension; molecular simulation; temperature; pressure; hydrates

date：2016-12-30.

Prof. ZHAO Lingling, zhao_lingling@seu.edu.cn

supported by the National Natural Science Foundation of China (51106027).

TQ 021.2

A

0438—1157（2017）07—2880—06

10.11949/j.issn.0438-1157.20161837

2016-12-30收到初稿，2017-03-14收到修改稿。

聯(lián)系人：趙伶玲。

季佳圓（1993—），女，碩士研究生。

國家自然科學(xué)基金項目（51106027）。