合作客戶/
拜耳公司 |
同濟(jì)大學(xué) |
聯(lián)合大學(xué) |
美國保潔 |
美國強生 |
瑞士羅氏 |
相關(guān)新聞Info
-
> 熱力學(xué)模型計算MgO-B2O3-SiO2-CaOAl2O3富硼渣表面張力（一）
> 人工沖洗升級為超聲波清洗，可改善新能源電池沖壓配件的表面張力
> 聚氧乙烯鏈長度調(diào)控非離子Gemini表面活性劑的表面張力、接觸角（四）
> 納米生物質(zhì)體系性能評價及驅(qū)油特性實驗研究
> 表面張力的測定方法有哪些？
> 兩種烷基咪唑亞磷酸酯離子液體熱穩(wěn)定性、表面張力測定（二）
> 生物表面活性劑產(chǎn)生菌菌體密度、細(xì)胞疏水性與發(fā)酵液pH及表面張力的關(guān)系（二）
> 基于界面張力和表面張力測試評估商用UV油墨對不同承印紙張的表面浸潤性差異（二）
> 堿、表面活性劑濃度變化對三元復(fù)合體系界面張力的影響
> 液滴撞擊固體表面過程中黏度與表面張力系數(shù)變化（一）
推薦新聞Info
-
> FYXF-3煤粉懸浮劑潤濕吸附性能、?傷害性能及在煤層氣壓裂改造現(xiàn)場的實施方案（三）
> FYXF-3煤粉懸浮劑潤濕吸附性能、?傷害性能及在煤層氣壓裂改造現(xiàn)場的實施方案（二）
> FYXF-3煤粉懸浮劑潤濕吸附性能、?傷害性能及在煤層氣壓裂改造現(xiàn)場的實施方案（一）
> 基于水煤漿流變性和動態(tài)表面張力觀察水煤漿的微觀破裂特性（三）
> 基于水煤漿流變性和動態(tài)表面張力觀察水煤漿的微觀破裂特性（二）
> 基于水煤漿流變性和動態(tài)表面張力觀察水煤漿的微觀破裂特性（一）
> 免罩光水性素色面漆配方、制備方法及步驟
> 固體、鹽溶液表面張力測量及與其在潔凈硅橡膠表面接觸角的關(guān)系研究（三）
> 固體、鹽溶液表面張力測量及與其在潔凈硅橡膠表面接觸角的關(guān)系研究（二）
> 固體、鹽溶液表面張力測量及與其在潔凈硅橡膠表面接觸角的關(guān)系研究（一）
溫度及壓強對CO2-NaCl鹽水系統(tǒng)界面張力的影響（一）
來源：化工學(xué)報 瀏覽 253 次 發(fā)布時間:2025-05-13
在超臨界態(tài)CO2封存于深部鹽水層過程中，溫度、壓強等控制條件是影響封存效率和封存量的重要因素。應(yīng)用分子動力學(xué)模擬的方法對343～373 K和6～35 MPa范圍內(nèi)的CO2-NaCl鹽水系統(tǒng)進(jìn)行了界面張力（IFT）及界面特性的研究，分析了IFT隨溫度及壓強的變化關(guān)系，并觀測到了壓力平衡點pplateau；從分子尺度（物質(zhì)密度、界面過余量、界面水合物密度）分析了IFT隨壓強、溫度的變化，以及pplateau產(chǎn)生的原因。結(jié)果表明，pplateau前壓強升高或溫度降低將導(dǎo)致CO2密度升高，IFT下降，而pplateau后IFT趨于穩(wěn)定且受溫度影響較?。籆O2的界面過余量及界面處水合物數(shù)量隨壓強及溫度變化，與IFT的變化相反；高壓下界面水合物密度的飽和現(xiàn)象可能是pplateau產(chǎn)生的重要原因。
引言
深部鹽水層CO2地質(zhì)封存所需的注射能耗及最大地質(zhì)埋存深度與CO2-鹽水之間的界面張力（interfacial tension，IFT）直接相關(guān)，并受溫度和壓強的制約。開展溫度、壓強對CO2-鹽水間界面張力的影響研究，不僅可以分析IFT隨溫度、壓強等控制參數(shù)的變化規(guī)律，還能闡述溫度、壓強對IFT產(chǎn)生影響的內(nèi)在機理進(jìn)而對指導(dǎo)不同環(huán)境條件（溫度、壓強）下的CO2地質(zhì)封存設(shè)計，提高注射安全性及存儲容量具有重大意義。
目前實驗已測定相關(guān)儲層條件下CO2-水和CO2-鹽水系統(tǒng)的IFT值，并觀測到IFT在定溫條件下會隨著壓強升高而降低，并在壓力平衡點pplateau之后趨于穩(wěn)定值。實驗還發(fā)現(xiàn)pplateau的大小與鹽的種類及鹽度無明顯聯(lián)系，僅隨溫度升高而上升。Chalbaud等將pplateau的存在歸因于CO2溶解度的影響，但尚未展開深入分析。
分子動力學(xué)模擬（molecular dynamics simulation，MD模擬）可以研究多相界面系統(tǒng)的微觀特性，目前該方法已成功模擬了CO2-水及CO2-鹽水系統(tǒng)，可獲得與實驗一致的IFT值，并能觀測界面的微觀現(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對溫度及壓強依賴關(guān)系的物理機理，尤其對pplateau現(xiàn)象的產(chǎn)生原因進(jìn)行了分析，可為IFT的控制和預(yù)測提供理論依據(jù)。
1研究對象及方法
1.1對象及模型
本文根據(jù)Chalbaud等對CO2-NaCl系統(tǒng)大范圍溫度及壓強下的IFT實驗研究結(jié)果，選擇了溫度及壓強范圍為343～373 K和6～35 MPa的CO2-NaCl系統(tǒng)為研究對象。具體工況參數(shù)列于表1。
表1 CO2-NaCl系統(tǒng)的溫度和壓強條件
在計算過程中，綜合考慮了系統(tǒng)內(nèi)分子間非鍵結(jié)作用力（范德華力、庫倫靜電力）及分子內(nèi)鍵結(jié)作用力（鍵拉伸和鍵彎曲）。分別采用Lennard-Jones勢能函數(shù)模擬范德華力，庫侖定律模擬庫侖靜電力，具體分子間勢能函數(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ù)模擬分子間長程庫侖作用力，范德華作用截距設(shè)定為0.9 nm。此外，本文采用諧波勢能函數(shù)模擬鍵拉伸和鍵角彎曲等分子內(nèi)鍵結(jié)作用力。
圖1 CO2-NaCl系統(tǒng)平衡狀態(tài)
本文計算中，水分子選擇柔性F3C模型，CO2選擇柔性EPM2模型，鹽離子采用Chandrasekhar等開發(fā)的模型。應(yīng)用MD軟件Gromacs4.5并采用周期性邊界條件進(jìn)行計算，所建立的橫截面4 nm×4 nm的計算域示于圖1。計算域中間區(qū)域為鹽水，包括4323個水分子、147個Na+和147個Cl-，對應(yīng)鹽度為1.89 mol·L-1，兩側(cè)分別為732個CO2分子。利用Berendsen方法來實現(xiàn)溫度和壓強的設(shè)定。由于系統(tǒng)在NPzT系綜下20 ns達(dá)平衡態(tài)，故模擬時間運行30 ns，選取最后5 ns為有效數(shù)據(jù)進(jìn)行分析。