陳圓圓 于炳松

(中國地質(zhì)大學地球科學與資源學院 北京 100083)

碳酸鹽巖溶解-沉淀熱力學模型及其在塔北地區(qū)的應用①

陳圓圓 于炳松

(中國地質(zhì)大學地球科學與資源學院 北京 100083)

碳酸鹽巖作為一種化學巖類其物性受成巖作用的影響十分明顯,溶解作用是否發(fā)育直接關系到碳酸鹽巖儲層質(zhì)量的好壞。碳酸鹽礦物在埋藏環(huán)境中的溶解度與地層水的化學成分、環(huán)境的溫度、pH等有著密切的關系。利用地層水化學資料,根據(jù)化學熱力學的理論,通過方解石及白云石與地層水作用的平衡反應,計算吉布斯自由能變化ΔG,以此來判斷方解石及白云石的溶解趨勢,并建立了封閉系統(tǒng)中這兩種礦物與水作用的熱力學模型。根據(jù)這個理論模型,對塔北地區(qū)159口井奧陶系及6口井的寒武系地層水進行了ΔG的計算,結果顯示,塔北地區(qū)碳酸鹽巖與地層水反應的ΔG基本為負值,表明該地區(qū)總體處于碳酸鹽巖溶解的環(huán)境,寒武系白云巖反應的自由能小于奧陶系灰?guī)r,即在這種深埋條件下白云巖比灰?guī)r更易溶解,這一點在實際巖芯資料上及前人的實驗模擬中均有顯示。塔北于奇地區(qū)自由能的計算結果也為負值,說明該地區(qū)強烈的充填作用并非發(fā)生在現(xiàn)今條件下,于奇地區(qū)現(xiàn)今的埋藏環(huán)境是有利于碳酸鹽礦物溶解的,這一結論也在該地區(qū)奧陶系及寒武系巖芯上得到了驗證。

碳酸鹽礦物 化學熱力學 溶解作用 吉布斯自由能

0 引言

碳酸鹽巖在地殼中分布極廣,不僅是有用的礦產(chǎn),被廣泛用于建筑、裝飾等工業(yè),更是重要的油氣儲集體。長期以來,國內(nèi)外對于碳酸鹽巖儲層的大量研究成果已證實了巖溶作用對于儲層形成的重要性[1,2]。古巖溶作用形成的各種溶蝕孔洞在后期埋藏中的保存機理以及現(xiàn)代埋藏條件下巖溶發(fā)育趨勢的研究,仍是學術界的一個難題。而這方面的研究,將是決定溶蝕形成的各種孔洞能否在埋藏條件下成為有效儲集空間的關鍵。碳酸鹽巖的溶解與沉淀的實質(zhì)是水-巖-氣之間的相互物理化學作用,其作用的結果主要表現(xiàn)在:碳酸鹽巖礦物的溶解形成孔洞,而沉淀則會造成孔洞的充填。以化學熱力學為基礎的相平衡研究目前是巖石學發(fā)展的前沿,在推動巖石學由定性到定量的研究中起著重要作用。碳酸鹽礦物在一定物理化學條件下是溶解還是沉淀可以用化學熱力學參數(shù)定量表示。熱力學是研究能量相互轉換過程中應遵循的規(guī)律的科學,主要研究一定條件下某種過程能否自發(fā)進行以及進行程度的問題。熱力學的一切結論基本是建立在熱力學三大定律的基礎上,這三個熱力學定律雖不能用理論方法加以證明,但其正確性已由無數(shù)次的實驗事實所證實。應用熱力學基本原理解決在化學過程及與化學有關的物理過程中的能量效應、反應可行性及被研究物質(zhì)的穩(wěn)定性等一系列問題上具有重要意義。在熱力學上,化學反應的吉布斯自由能增量(ΔG)可作為熱力學過程的方向和限度的判據(jù),還可作為過程不可逆性大小的量度。當ΔG>0,過程不可能自動發(fā)生;ΔG=0,過程平衡;ΔG<0,過程自動發(fā)生(不可逆),ΔG值越低,說明自動過程越易發(fā)生,而且反應越快[3].elgeson等[4]1978年對造巖礦物的熱力學性質(zhì)進行了總結,給出各種礦物及離子的熱力學數(shù)據(jù),成為相平衡研究的重要數(shù)據(jù)庫.ribble等[5]1995年研究了白云石等礦物的晶體化學及熱力學性質(zhì).ierre Genthon等[6]1997年研究了碳酸鹽礦物在碎屑巖儲層中對次生孔隙發(fā)育的貢獻.ohn W Morse[7]等2002年探討了幾種碳酸鹽礦物的溶解動力學。國內(nèi),林傳仙等1985年建立了礦物及有關化合物的熱力學數(shù)據(jù)庫[8],王高尚等1991年將Helgeson等人提出的估算高溫高壓條件下水溶液中電解質(zhì)和離子的標準偏摩爾性質(zhì)的HKF模型做了詳細介紹[9]。前人的眾多成果使得各種化學反應的熱力學計算有了數(shù)據(jù)支持。近年,于炳松、劉建清等研究了地層水條件下碳酸鹽礦物熱力學平衡條件并將其應用于油氣田的儲層分析,取得了良好的效果[10～14]。

本文建立了碳酸鹽巖體系中水-巖作用的化學熱力學模型,通過計算現(xiàn)代埋藏條件下方解石和白云石化學反應的吉布斯自由能ΔG值,提供一種判別現(xiàn)代埋藏條件下碳酸鹽巖巖溶發(fā)育與保存趨勢的方法,并通過在塔北地區(qū)的實際應用,對該方法進行了初步的應用效果驗證,旨在探索碳酸鹽巖儲層研究的新思路。

1 化學熱力學模型

熱力學在解決問題時所用的方法是嚴格的數(shù)理邏輯推理方法,具有以下特點:首先,熱力學的研究對象是具有足夠大量質(zhì)點的系統(tǒng),熱力學只研究物質(zhì)的宏觀性質(zhì),對于物質(zhì)的微觀性質(zhì)即個別或少數(shù)分子、原子的行為無法解釋;其次,熱力學只需知道系統(tǒng)的起始狀態(tài)和最終狀態(tài)以及過程進行的外界條件,就可以進行相應的計算,它不依賴于物質(zhì)結構的知識,也無需了解過程進行的機理,由于這個原因,對過程自發(fā)性的判斷,只能是表面了解而不知其內(nèi)在原因;在熱力學研究的變量中,沒有時間的概念,所以它不涉及過程進行的時間問題,所以只能說明過程能不能自發(fā)進行以及進行到什么程度,至于過程發(fā)生的時間及具體的速率,熱力學無法預測。這些既是熱力學方法的優(yōu)點也是它的局限性[3]。

熱力學模型的建立通常涉及到系統(tǒng)、環(huán)境、狀態(tài)、相、熱力學平衡幾個概念。研究的對象物質(zhì)統(tǒng)稱為系統(tǒng),將系統(tǒng)之外的但和系統(tǒng)之間關系密切的其余物質(zhì)和空間稱為環(huán)境。系統(tǒng)的物理性質(zhì)和化學性質(zhì)的綜合表現(xiàn)稱為狀態(tài),系統(tǒng)狀態(tài)的性質(zhì)稱為狀態(tài)性質(zhì)或狀態(tài)函數(shù),如溫度、壓力、焓、熵等,應注意系統(tǒng)的熱力學狀態(tài)性質(zhì)只說明系統(tǒng)當時所處的狀態(tài),而不能說明系統(tǒng)以前的狀態(tài)[15]。反應處于熱力學平衡態(tài)時,系統(tǒng)與環(huán)境之間沒有任何物質(zhì)和能量交換,系統(tǒng)中各個狀態(tài)性質(zhì)均不隨時間變化。

考慮到熱力學研究的特點,本文選擇已知溫度的封閉系統(tǒng)作為平衡體系,不考慮系統(tǒng)的物質(zhì)遷移轉化、熱量的進出,對碳酸鹽巖在現(xiàn)今埋藏環(huán)境中溶解反應的發(fā)生方向及反應程度進行了預測。

1.1 方解石-水相互作用熱力學模型

方解石在地層水中主要存在如下電離方程式:

任意化學反應的吉布斯自由能ΔG可表示為式中ΔrG0為反應的標準吉布斯自由能變化,Qa為活度商,表示為Qa=∏B(aB)υB,B表示反應式中的任一物質(zhì),aB為溶液中離子的活度,υB為該物質(zhì)的化學計量數(shù),Keq為反應的平衡常數(shù)。因此對于方解石溶解沉淀的反應有如下方程式:

其中aCa2+為現(xiàn)今實際溶液中Ca2活度積,計算公式如下:

其中[Ca2+]為現(xiàn)今溶液中Ca2+濃度,γCa2+為Ca2+活度系數(shù)。

Debye和Huckel于1923年引入了離子云的概念,用數(shù)學方法導出了離子溶液中活度系數(shù)的計算公式:

式中γi為第i種組分的活度系數(shù);A、B為主要取決于水的溫度T(℃)的常數(shù),其函數(shù)關系可近似表示為A=0.0009T+0.4869,B=0.0002T+0.3241;ai為與離子水化學半徑有關的常數(shù);I為溶液的離子強度,其定義式為I=mizi,M為水溶液中離子組分的個數(shù),mi為第i種離子組分的濃度(mol/L),zi為第i種離子所帶的電荷數(shù)[16]。大量實驗表明 Debye-Huckel公式僅適用于離子強度小于0.1的溶液.resdell和Jones(1974)修改了Debye-Huckel公式使計算結果在離子強度不大于6的范圍內(nèi)均是可信的[16],公式如下:

表1 碳酸鹽巖地層水中常出現(xiàn)離子的活度系數(shù)參數(shù)Table1 Activity coefficients parameters of common ions in carbonate strata water

(2)式中aCa2+為反應達到平衡時Ca2+活度積,T為水樣所在深度的溫度。稀溶液中物質(zhì)的活度近似于濃度(mol.L-1),即ai=[i]。水的活度可近似為1。即反應的進行程度主要取決于現(xiàn)今地下水中的Ca2+活度、反應平衡態(tài)時Ca2+的理論活度以及實際地下水的溫度。

反應平衡態(tài)時Ca2+活度可以根據(jù)現(xiàn)今地層水溶液中離子濃度及一定溫度下反應平衡常數(shù)求得。反應前后整個固液體系遵從質(zhì)量守恒原則,溶液中Ca2+變化的數(shù)量應與碳酸總量的變化相等,即

由上述各反應平衡常數(shù)KCc、K1、K2與反應式中各離子活度的函數(shù)關系可以推導出

將式(7)代入式(4)中可得:

根據(jù)Van't Hoff公式可計算化學反應平衡常數(shù)隨溫度的變化,在恒壓條件下有

表2 碳酸鹽溶解反應中常見組分的標準焓變Δr H 0mTable2 The standard enthalpy changeΔr H0m of common components in carbonate dissolve reaction

將(3)式及(8)式代入(2)式中即可求得ΔG。

1.2 白云石-水相互作用熱力學模型

埋藏成巖作用中白云石(CaMg(CO3)2)的溶解和沉淀取決于埋藏環(huán)境中流體與礦物之間發(fā)生的各種物理化學反應。對于白云巖與地層水的封閉體系,存在如下電離方程式:

與求解方解石溶解沉淀反應的自由能原理相同,白云石與水作用過程的自由能變化可由下面的公式表示:

其中aCa2+aMg2+為實際溶液中Ca2+和Mg2+活度積,可用Debye-Huckle模型計算得出。

γCa2+和γMg2+為活度系數(shù),可由Debye-Huckel公式求出。

aCa2+eqaMg2+eq為反應達到平衡時 Ca2+和 Mg2+活度積,T為水樣所在深度的溫度。稀溶液中物質(zhì)的活度近似于濃度(mol/L),即ai=[i]。水的活度可近似為1。

由白云石在水中反應的平衡常數(shù)推導可得:

將(13)(14)(15)代入(10)式可得:

反應達到平衡時,反應混合物中各組分的濃度保持不變,反應式中生成物的濃度之比等于化學計量數(shù)之比,因此有

將(16)(17)式代入(9)式可得:

最后將(12)(18)代入(11)式即可求出ΔG。

1.3 模型適用條件

本文中方解石及白云石溶解沉淀的熱力學模型均建立在封閉系統(tǒng)內(nèi),本質(zhì)為計算恒溫定壓條件下方解石、白云石與含CO2的地層水發(fā)生化學反應過程中做非體積功的最大限度,即吉布斯自由能的變化。本文模型是在已知離子濃度及反應平衡常數(shù)K的前提下進行計算的,K主要隨溫度的變化而變化,大量實驗證明壓力對其也有一定影響,因此系統(tǒng)必須是恒溫恒壓的環(huán)境,且溫度、壓力可求。

計算離子活度時應用的 Debye-Huckel模型,是建立在非締合電解質(zhì)溶液中僅存在庫倫力、水為連續(xù)介質(zhì)等假設的基礎上,認為離子服從Boltzman分布,并采用Possion方程處理離子氛電位,從而推導出離子的平均活度系數(shù)方程,該理論是第一個嚴格導出的有解析式的電解質(zhì)溶液理論,適用于離子強度I< 0.1 mol/L的極稀電解質(zhì)溶液[18],Tresdell和Jones在此基礎上添加新參數(shù)而引出的修正模型,使之在I< 6 mol/L的溶液中均適用[16]。

2 熱力學模型在塔北地區(qū)的應用

采用化學熱力學方法計算吉布斯自由能ΔG正負值可以判別方解石及白云石溶解與沉淀趨勢,理論的正確性毋庸置疑,但這種方法在巖溶儲層預測中是否具有可操作性或者說能否切合實際,并最終取得良好的效果,需要用實際的資料與計算結果進行吻合驗證。塔北地區(qū)鉆探程度較高,地層水資料豐富,碳酸鹽巖地層在該區(qū)埋深較大。而在深埋條件下,流體運移緩慢,可作為一個封閉的系統(tǒng),在一定時間段內(nèi)特定的深度區(qū)域之間其溫度及壓力基本保持不變,地層溫度及壓力均可根據(jù)鉆井資料求出,因此塔北地區(qū)碳酸鹽巖儲層具備本文熱力學模型的應用條件。

2.1 塔北地區(qū)地質(zhì)概況

塔河油田位于塔里木盆地塔北隆起區(qū)阿克庫勒凸起中段的西南斜坡,西鄰哈拉哈塘凹陷,東靠草湖凹陷,南接順托果勒低凸起和滿加爾凹陷,北部為雅克拉斷凸[19]。塔里木盆地北部存在加里東中晚期-海西早期、海西晚期-印支期、燕山期-喜馬拉雅早期及喜馬拉雅晚期四大水文地質(zhì)旋回。塔河油區(qū)處于阿克庫勒凸起南斜坡,具有長期相對穩(wěn)定的有利區(qū)域水動力背景。加里東-海西早期水文地質(zhì)環(huán)境造就了奧陶系主要是下奧陶統(tǒng)規(guī)模巨大的碳酸鹽巖巖溶-縫洞型儲集體;海西晚期以來,滿加爾、草湖烴源巖區(qū)的沉積壓實水為油氣運移提供了動力與載體;喜馬拉雅晚期由于塔北上古生界-中新生界地層反轉,造成部分中生界成藏油氣的調(diào)整,但未改變下古生界宏觀南傾的格局;晚期油氣依然沿南傾的T、T等不整合面由南部及東部烴源區(qū)向塔河運移聚集[20]。地層水資料統(tǒng)計結果表明,塔河油田奧陶系地層水以CaCl2水型為主,偶見Na2SO4水型,礦化度在150～ 300 g/L之間。

塔河油田奧陶系的易溶巖層主要位于中-下奧陶統(tǒng),為一套較為純凈的碳酸鹽巖;而上奧陶統(tǒng)則混雜了較多的泥質(zhì),不利于巖溶作用的進行[25](圖1)。研究區(qū)大部分鉆井都揭示了一間房組及鷹山組,從巖性上看,主體都為顆?；?guī)r及泥晶灰?guī)r,其中一間房組泥晶灰?guī)r-生物礁灰?guī)r-顆粒灰?guī)r-藻粘結灰?guī)r的旋回式發(fā)育特征較為明顯,而鷹山組向下白云石化呈增強趨勢。塔河油田目前鉆遇寒武系的井數(shù)量較少,且基本停留在上寒武丘里塔格組,發(fā)育一套微晶-細晶白云巖。

2.2 化學熱力學分析

2.2.1 方解石溶解-沉淀化學熱力學分析

為了解塔河油田奧陶系深埋藏碳酸鹽巖儲集物性與方解石溶解-沉淀、孔隙水化學之間的成因聯(lián)系,本文選取了塔河油田159口井,利用埋藏條件下地下水-方解石相互作用的熱力學原理,對這些井深5 000～5 500 m的奧陶系碳酸鹽巖儲集層在現(xiàn)今埋藏地下水作用下的巖溶趨勢進行研究。研究區(qū)奧陶系儲集巖主要為灰?guī)r和白云質(zhì)灰?guī)r,孔隙度及滲透率均較高,地下水成分為塔里木油田分公司實驗中心用分光光度計測得。

圖1 塔河地區(qū)寒武-奧陶系地層柱狀圖(據(jù)嚴威等,2011)[25]Fig.1 Synthetic columnar of the Cambrian-Ordovician formations in Tahe area

根據(jù)鉆井測溫資料,塔河地區(qū)平均地溫梯度約2.2℃/100m,塔河油田奧陶系地層溫度在120～ 130℃之間。地層壓力采用以下公式計算[12]:

塔河油田奧陶系泥漿當量密度大約為1.08～ 1.12 g/cm3,前人通過大量的數(shù)據(jù)點計算認為塔河油田奧陶系儲層的地層壓力約在50～100 MPa之間。根據(jù)塔河油田的溫壓條件,選擇在100℃及100 MPa條件下的各反應平衡常數(shù)(表3),以此來計算反應平衡狀態(tài)下的Ca2+eq活度。經(jīng)過初步計算,地層水中離子強度在2～5之間,因此現(xiàn)今地層水中的Ca2+活度選用Tresdell和Jones修改后的Debye-Huckel的活度模型來計算。根據(jù)塔河油田基本覆蓋全區(qū)的159口井的地層水各項值計算得到了現(xiàn)代條件下方解石在這些井中的巖溶趨勢(△G),其中鷹山組96口井,一間房組63口井,并應用Surfer軟件繪制了相應的溶解趨勢等值線圖(圖2、3)。

表3 碳酸鹽巖反應的平衡常數(shù)Table3 The equilibrium constants of the carbonate reactions

2.2.2 白云石溶解-沉淀化學熱力學分析

為了了解白云石溶解-沉淀與碳酸鹽巖儲集物性的關系,本文選擇了塔北地區(qū)寒武系白云巖作為研究對象。塔北地區(qū)寒武系埋藏深度大,多在6 000 m以下,鉆揭寒武系的探井也很少。據(jù)現(xiàn)有資料,選取了6口揭示寒武系丘里塔格組的鉆井進行分析。丘里塔格組以大套結晶白云巖及殘余顆粒白云巖為主,是寒武系中最好的碳酸鹽儲層[21]。該組泥漿當量密度大約為1.02～1.16 g/cm3,與奧陶系大致相同,本文選擇在100℃及100 MPa條件下的各反應平衡常數(shù)來進行計算。經(jīng)過初步計算,寒武系地層水中離子強度在1～3之間,因此現(xiàn)今地層水中的 Ca2+及Mg2+活度也選用Tresdell和Jones修改后的 Debye-Huckel的活度模型來計算。

圖2 塔河地區(qū)鷹山組埋藏巖溶趨勢圖Fig.2 Burial karst trend contourmap of Yingshan group in Tahe area

白云石是一種較為特殊的碳酸鹽礦物,其活度的計算要考慮晶體離子占位的問題。理想白云石晶體的化學式為CaMg(CO3)2,其中,Ca和Mg離子分別主要分布在M2和M1位上。基于Bragg-Williams理論模型,有s=2XM2-1=2XM1-1。其中,XM2、XM1分別指M 2位和M 1位上Ca和Mg的摩爾分數(shù).=0時為完全無序白云石,s=1時為有序白云石,s為白云石的有序度。白云石晶體中有序白云石組分(aoD)和無序白云石組分(adD)的活度分別為[12]

圖3 塔河地區(qū)一間房組埋藏巖溶趨勢圖Fig.3 Burial karst trend contourmap of Yijianfang group in Tahe area

對于接近理想成分的白云石晶體來說,其s主要決定于溫度.oldsmith和Heard(1961)作出了0～ 1200℃之間白云石有序度與溫度的關系圖(圖4)[4]。研究區(qū)寒武系地層溫度在400K左右,有序度s≈1,因此白云石活度aDol=1。將其代入白云石-水相互作用的熱力學模型,計算結果如表4所示。

2.3 討論

2.3.1 奧陶系灰?guī)r埋藏巖溶趨勢分析

根據(jù)塔河地區(qū)奧陶系鷹山組及一間房組現(xiàn)今埋藏條件下方解石溶解-沉淀反應的ΔG平面分布情況(圖2、3)可知,塔河奧陶系方解石溶解反應的ΔG基本都小于0,表明該地區(qū)現(xiàn)今的埋藏條件整體是有利于方解石溶解的環(huán)境。如圖5所示,T704井、T705井縫合線溶蝕及膠結物溶解的現(xiàn)象,說明埋藏溶蝕作用在該區(qū)較為發(fā)育。局部ΔG接近0或大于0的區(qū)域位于塔河油田北部、于奇地區(qū)的南部,呈北東向展布。

圖4 白云石有序度與溫度的函數(shù)關系[4]Fig.4 Function relationship between dolomite order degree and temperature

表4 寒武系地層水離子強度、離子活度與吉布斯自由能Table4 The ionic intensity,ionic activity and Gibbs free energy of Cam brian oilfield water

對于奇地區(qū)鉆井資料進行分析可知,該地區(qū)大量的縫洞均被充填,但是在巖芯觀察中,可以發(fā)現(xiàn)原來被充填的縫或洞又發(fā)生了溶蝕現(xiàn)象,如圖6所示,說明盡管地質(zhì)歷史時期于奇地區(qū)大量的方解石發(fā)生沉淀,但在現(xiàn)今埋藏條件下,方解石與地層水的反應向發(fā)生溶解的方向進行。在鷹山組埋藏巖溶趨勢圖中,于奇地區(qū)自由能在-3 kJ/mol左右,這也表明計算結果與實際資料吻合性較好。

圖5 塔河地區(qū)鉆井埋藏巖溶現(xiàn)象a.亮晶砂屑灰?guī)r,膠結物發(fā)生溶蝕.T705井,5661.55m,X60;b.沿縫合線溶蝕.T704井,5740.31m,X60Fig.5 Buried karst phenomenon in Tahe area

圖6 于奇地區(qū)巖芯縫洞溶蝕特征a.于奇1井O1-2 y淺黃灰色泥晶灰?guī)r,充填溶蝕縫內(nèi)發(fā)育溶孔;b.于奇6井O1-2 y黃灰色泥晶砂屑灰?guī)r,被方解石全充填的溶洞發(fā)育溶孔,未充填.Fig.6 Dissolution characteristics on cores in Yuqi area

巖溶作用及其孔洞縫網(wǎng)絡體系的發(fā)育與巖溶古地貌關系密切。巖溶高地是巖溶大氣水的區(qū)域補給區(qū),地下水動力以垂向滲入滲流帶為主,水平潛流帶欠發(fā)育,決定了儲滲空間以溶蝕裂縫、孤立孔洞發(fā)育為特征,潛水面洞穴層和洞穴層儲層欠發(fā)育;而巖溶谷地的地下水水動力分帶不明顯,地下水的動力和化學溶蝕能力均已減弱;巖溶斜坡相帶是連接巖溶高地和巖溶谷地的過渡帶,地下滲流和水平潛流均發(fā)育,決定了該帶是巖溶儲層發(fā)育的最好區(qū)帶。在塔河地區(qū),海西早期阿克庫勒凸起沖斷褶皺強烈,下奧陶統(tǒng)頂部及上覆地層剝蝕夷平,塔河地區(qū)巖溶斜坡發(fā)育。在奧陶紀末至早石炭世,阿克庫勒凸起以北部高地、東西兩側的谷地以及中部的丘叢-洼地相間的古地貌形態(tài)出現(xiàn),使得塔河地區(qū)奧陶系碳酸鹽巖發(fā)育大規(guī)模、多旋回巖溶帶。由圖7可知,塔河地區(qū)奧陶系巖溶高地呈北東向延伸,這與圖2、圖3中的ΔG相對高值區(qū)分布情況有著良好的對應關系。

鷹山組和一間房組在塔河地區(qū)大部分鉆井中均有揭示,而于奇地區(qū)則缺失一間房組。鷹山組屬于典型的碳酸鹽臺地相沉積,巖性以泥晶灰?guī)r為主,夾雜部分顆?；?guī)r。一間房組主要為臺地淺灘-臺內(nèi)礁相沉積,巖性表現(xiàn)為泥晶灰?guī)r、生物丘灰?guī)r、顆?；?guī)r均有發(fā)育,其中砂屑灰?guī)r、生屑灰?guī)r等顆?；?guī)r占主要優(yōu)勢。受不同地區(qū)沉積微相、沉積后成巖作用及構造運動的影響,奧陶系中下統(tǒng)灰?guī)r段非均質(zhì)性較為顯著[21]。在圖2、圖3上表現(xiàn)為,塔河油田南部地區(qū)鷹山組的ΔG明顯小于一間房組,而部分地區(qū)(T727-T705井區(qū))鷹山組的ΔG則大于一間房組。但整體上顯示出鷹山組巖溶趨勢略強于一間房組,這在鷹山組及一間房組巖芯溶蝕孔洞發(fā)育程度上也有所體現(xiàn)。

圖7 塔河地區(qū)奧陶系巖溶地貌Fig.7 Ordovician karst landform in Tahe area

2.3.2 寒武系白云巖埋藏巖溶趨勢分析

根據(jù)表4可知,寒武系丘里塔格組吉布斯自由能在-6～-8 kJ/mol之間,屬于典型的白云石溶蝕環(huán)境。分析表中6口井資料也可以發(fā)現(xiàn),這些井極為發(fā)育溶蝕孔洞,且基本未充填,如圖7所示。

以YQ6井為例,如前文所述,該地區(qū)大量的縫洞被充填,但是這種充填作用并不是發(fā)生在現(xiàn)今條件下,根據(jù)YQ6井寒武系的巖芯觀察,在全充填的縫洞上發(fā)育了許多未充填或半充填的溶孔(圖9a),此外如圖8c,掃描電鏡下,白云石礦物發(fā)生溶蝕,這表明現(xiàn)今條件該地區(qū)地層水是有利于白云石溶解的。塔深1井(TS1)是目前鉆探深度最大的井,計算得出該井上寒武統(tǒng)吉布斯自由能為-8.14 kJ/mol,根據(jù)熱力學原理,吉布斯自由能越小反應越容易進行,而該井在埋深8 000多米的上寒武統(tǒng)白云巖中揭示了蜂窩狀溶蝕孔洞的存在(圖9b),而且在溶洞中還發(fā)現(xiàn)了黃褐色的液態(tài)烴,恰好驗證了這一原理。

與奧陶系灰?guī)r溶解的吉布斯自由能相比,寒武系白云巖自由能值明顯較低,這就表示在壓力100 MPa、溫度100℃條件下白云巖比灰?guī)r更容易發(fā)生溶解。從圖6b及圖9a上即可看出,同一口井的寒武系白云巖次生孔隙明顯多于奧陶系灰?guī)r。楊俊杰、黃思靜等1995年至2003年通過大量實驗模擬了表生及埋藏條件下碳酸鹽巖溶蝕過程,實驗證明表生與相對淺埋藏的溫壓條件下,石灰?guī)r的巖溶作用較白云巖發(fā)育;但在深埋藏階段,由溶解作用造成的白云巖次生孔隙應比方解石更為發(fā)育[22～24]。本文理論計算判別結果與實際資料溶蝕孔洞發(fā)育情況及實驗模擬結果都具有較好的一致性。由此可知,應用化學熱力學方法判別巖溶的發(fā)育與保存是可靠而有效的。

CaMg(CO3)2在現(xiàn)代海水中呈過飽和狀態(tài),但在常溫常壓下無法從海水中直接沉淀出白云石,所以白云石的溶解沉淀不僅僅是熱力學問題,更是動力學問題。根據(jù)化學動力學理論,在任何反應中,并不是所有的分子都能參加反應,而是具有一定能量水平的分子(離子)才能參加反應,這些分子稱為活化分子,活化分子的平均能量與所有分子的平均能量的差叫做活化能。反應的發(fā)生首先需要克服活化能,因此即使根據(jù)熱力學原理能夠自發(fā)進行的反應也未必能夠發(fā)生。克服活化能的方法包括升高溫度及提高反應物的濃度,從而提高分子有效碰撞次數(shù)[26]。塔北地區(qū)寒武系埋深基本都大于5 500 m,由地溫梯度換算其地層溫度都在120℃以上,地層水為高礦化度鹵水。這種深埋藏的物理化學環(huán)境有利于白云石溶解沉淀反應克服活化能,因此在符合該反應發(fā)生的動力學條件下,應用熱力學手段判斷反應進行的方向也是可行的。

圖8 白云石溶蝕現(xiàn)象a.XH2井鑄體薄片5 516.27 m白云巖,X40;b.DG1井鑄體薄片6 258.81 m白云巖,X40; c.YQ6井7 118.58 m,白云石溶蝕,X500;d.DG1井6 256.31 m,晶間溶孔,X692Fig.8 Dissolution phenomenon of dolomite

圖9 寒武系巖芯照片a.YQ6井第4筒上寒武統(tǒng)7 115.32 m淺灰色中細晶白云巖;b.塔深1井第4筒上寒武統(tǒng)7 873.98 m白云巖溶蝕孔洞Fig.9 Cambrian core photo pictures

3 結論

(1)本文熱力學模型建立在溫度壓力可知的封閉系統(tǒng)內(nèi),不考慮物質(zhì)轉化遷移及熱量傳遞,是一個較為簡單的地球化學模型。推導過程主要應用了質(zhì)量守恒定律及Debye-Huckel模型,質(zhì)量守恒定律是自然界普遍存在的基本定律之一,而Debye-Huckel模型則不斷的被修正使其適用范圍更廣,因此了解模型的適用條件是十分必要的。

(2)根據(jù)熱力學模型的各種參數(shù)可以看出,埋藏條件下碳酸鹽巖溶解反應的吉布斯自由能主要受溫度、壓力、Ca2+和Mg2+離子濃度及孔隙流體pH的控制。由計算結果可知,塔北地區(qū)埋藏環(huán)境總體上有利于碳酸鹽巖溶解的,寒武系白云巖溶蝕反應比奧陶系灰?guī)r更易發(fā)生。

(3)于奇地區(qū)在地質(zhì)歷史時期發(fā)生過嚴重的充填作用,但現(xiàn)今埋藏條件對方解石及白云石以溶蝕作用為主,理論上該區(qū)吉布斯自由能表現(xiàn)為負值,在巖心及掃描電鏡分析中可見曾經(jīng)充填的孔、縫內(nèi)發(fā)生進一步溶蝕的現(xiàn)象,說明了計算結果的正確性。

(4)熱力學模型計算法在判別塔北地區(qū)碳酸鹽巖儲層次生孔隙發(fā)育中能很好的切合實際,說明化學熱力學可以作為一種新的輔助手段應用于油氣儲層預測中。

References)

1 梅博文譯。儲層地球化學譯文集[M]。西安:西北大學出版社, 1991:1-241[Translated by Mei Bowen.Reservoir Geochemistry (Translated Papers)[M].i'an:Northwestern University Press, 1991:1-241]

2 蔡春芳。沉積盆地流體-巖石相互作用研究的現(xiàn)狀[J]。地球科學進展,1996,11(6):575-578[CaiChunfang.Advancements of fluidrock interactions in sedimentary basin[J].dvance in Earth Science, 1996,11(6):575-578]

3 印永嘉,奚正楷,李大珍。物理化學簡明教程[M]。北京:高等教育出版社,1992:1-580[Yin Yongjia,Xi Zhengkai,Li Dazhen.The Concise Course of Physicochemistry[M].eijing:Higher Education Press,1992:1-580]

4 Helgeson H C,Delany JM,NesbittHW,etal.Summary and critique of the thermodynamic properties of rock-forming minerals[J].merican Journal of Science,1978,278-A:1-229

5 Tribble JS,Arvidson R S,Lane M III,et al.Crystal chemistry,and thermodynamic and kinetic properties of calcite,dolomite,apatite and biogenic silica:Applications to petrologic problems[J].hemical Geology,1995,95:11-37

6 Pierre Genthon,Jacques Schott,Jean-Louis Dandurand.Carbonate diagenesis during thermo-convection:Application to secondary porosity generation in elastic reservoirs[J].hemical Geology,1997,142:41-61

7 John W.Morse,Rolf S.Arvidson.The dissolution kinetics of major sedimentary carbonate minerals[J].arth-Science Reviews,2002, 58:51-84

8 王高尚,李鵬九。熱液平衡計算數(shù)據(jù)手冊[M]。北京:地質(zhì)出版社,1991:1-100[Wang Gaoshang,Li Pengjiu.The Data Handbook of Hydrothermal Equilibrium Calculation[M].eijing:Geological Publishing House,1991:1-100]

9 林傳仙,白正華,張哲儒。礦物及有關化合物熱力學數(shù)據(jù)手冊[M]。北京:科技出版社,1985:1-354[Lin Chuanxian,Bai Zhenghua,Zhang Zheru.The Thermodynamic Data Handbook ofMineral and Related Compounds[M].eijing:Science and Technology Press, 1985:1-354]

10 賴興運,于炳松,陳軍元,等。碎屑巖骨架顆粒溶解的熱力學條件及其在克拉2氣田的應用[J]。中國科學:D輯,2004,42:45-53[Lai Xingyun,Yu Bingsong,Chen Junyuan,etal.The thermodynamic condition of rock fragmentdissolution and its application in Kela 2 gas field[J].cience in China:Series D,2004,42:45-53]

11 于炳松,賴興運?？死?氣田儲集巖中方解石膠結物的溶解及其對次生孔隙的貢獻[J]。礦物巖石,2006,26(2):74-79[Yu Bingsong,Lai Xingyun.Dissolution of calcite cementand its contribution to the secondary pores of reservoir in the Kela 2 gas field in the Tarim Basin[J].ournal ofMineralogy and Petrology,2006,26(2): 74-79]

12 劉建清,賴興運,于炳松。地層水條件下碳酸鹽礦物熱力學平衡條件及其在克拉2氣田的應用[J]。沉積學報,2006,24(5): 636-640[Liu Jianqing,Lai Xingyun,Yu Bingsong.The thermodynamics equilibrium condition of carbonate minerals in strata environment and its application in Kela 2 gas field[J].cta Sedimentologica Sinica,2006,24(5):636-640]

13 于炳松,賴興運。成巖作用中的地下水碳酸體系與方解石溶解度[J]。沉積學報,2006,24(5):627-635[Yu Bingsong,Lai Xingyun.Carbonic acid system of groundwater and the solubility of calcite during diagenesis[J].cta Sedimentologica Sinica,2006,24 (5):627-635]

14 于炳松,賴興運,高志前。克拉2氣田砂巖儲層中成巖方解石-白云石的平衡及其對儲層質(zhì)量的影響[J]。自然科學進展,2007, 3:339-345[Yu Bingsong,Lai Xingyun,Gao Zhiqian.The equilibrium between calcite and dolomite during diagenesis in sandstone reservoir in Kela 2 gas field and the influence to quality of reservoir[J].rogress in Natural Science,2007,3:339-345]

15 陳家瑋,鮑征宇。地球化學反應平衡模型的建模方法[J]。地質(zhì)科技情報,2001,20(3):41-46[Chen Jiawei,Bao Zhengyu.Approaches to geochemical reaction equilibrium modeling[J].eological Science and Technology Information,2001,20(3):41-46]

16 錢會,馬致遠。水文地球化學[M]。北京:地質(zhì)出版社,2005:1-235[Qian Hui,Ma Zhiyuan.Hydrogeochemistry[M].eijing:Geological Publishing House,2005:1-235]

17 于炳松,林暢松。油氣儲層埋藏成巖過程中的地球化學熱力學[J]。沉積學報,2009,27(5):896-903[Yu Bingsong,Lin Changsong.Geochemical thermodynamics of diagenesis in reservoirs for oil and gas[J].cta Sedimentologica Sinica,2009,27(5):896-903]

18 覃東棉,劉桂華,李小斌,等.ebye-Hückel理論的研究進展[J]。材料導報,2010,24(13):107-112[Qin Dongmian,Liu Guihua,Li Xiaobin,etal.Research progress of Debye-Hückel theory [J].aterials Review,2010,24(13):107-112]

19 朱蓉,樓章華,金愛民,等。塔河油田S48縫洞單元流體分布及開發(fā)動態(tài)響應[J]。浙江大學學報:工學版,2009,43(7):1344-1348[Zhu Rong,Lou Zhanghua,Jin Aimin,etal.Fluid distribution and dynamic responses to exploitation in fracture-cave unit S48 in Tahe Oilfield[J].ournal of Zhejiang University:Engineering Science,2009,43(7):1344-1348]

20 顧憶,黃繼文,馬紅強。塔河油區(qū)油氣分布特點及其控制因素[J]。中國西部油氣地質(zhì),2006,2(1):19-25[Gu Yi,Huang Jiwen,Ma Hongqiang.Petroleum distribution and controlling factors in the Tahe Oilfield[J].est China Petroleum Geosciences,2006,2 (1):19-25]

21 葉德勝,王根長,林忠民,等。塔里木盆地北部寒武-奧陶系碳酸鹽巖儲層特征及油氣前景[M]。成都:四川大學出版社,2000: 1-151[Ye Desheng,Wang Gengchang,Lin Zhongmin,et al.Characteristics of Cambrian-Ordovician Carbonate Reservoir and Its Gas and Oil Foreground in the Northern Tarim Basin[M].hengdu:Sichuan University Press,2000:1-151]

22 楊俊杰,黃思靜,張文正,等。表生和埋藏成巖作用的溫壓條件下不同組成碳酸鹽巖溶蝕成巖過程的實驗模擬[J]。沉積學報, 1995,13(4):49-54[Yang Junjie,Huang Sijing,ZhangWenzheng, et al.Experimental simulation of dissolution for carbonate with different composition under the conditions from epigenesis to burial diagenesis environment[J].cta Sedimentologica Sinica,1995,13(4):49-54]

23 楊俊杰,張文正,黃思靜,等。埋藏成巖作用的溫壓條件下白云巖溶解過程的實驗模擬研究[J]。沉積學報,1995,13(3):83-87 [Yang Junjie,Zhang Wenzheng,Huang Sijing,et al.Experimental simulation for dolomite dissolution under the conditions of burial temperature and pressure[J].cta Sedimentologica Sinica,1995,13 (3):83-87]

24 肖林萍,黃思靜。方解石和白云石溶蝕實驗熱力學模型及地質(zhì)意義[J]。礦物巖石,2003,23(1):113-116[Xiao Linping,Huang Sijing.Model of thermodynamics for dissolution of carbonate and its geological significances[J].ournal of Mineralogy and Petrology, 2003,23(1):113-116]

25 嚴威,王興志,張廷山,等。塔河油田加里東中期第Ⅲ幕巖溶作用[J]。石油學報,2011,32(3):411-416[Yan Wei,Wang Xingzhi,Zhang Tingshan,etal.Investigations of the karst during the episodeⅢof the mid-Caledonian in the Tahe Oilfield[J].cta Petrolei Sinica,2011,32(3):411-416]

26 陳永權,周新源,楊文靜。白云石形成過程中的熱力學與動力學基礎及白云巖形成環(huán)境劃分[J]。海相油氣地質(zhì),2009,14(1): 21-25[Chen Yongquan,Zhou Xinyuan,Yang Wenjing.Thermodynamic and kinetic principle in dolomitizition process and the division of dolostone diagenetic environment[J].arine Origin Petroleum Geology,2009,14(1):21-25]

Dissolution-Precipitation Thermodynam ic M odels of Carbonate Rock and the Application in the Northern Part of Tarim Basin

CHEN Yuan-yuan YU Bing-song
(School of Earth Sciences and Resources,China University of Geosciences,Beijing 100083)

As a kind of chemical rock,the porosity and permeability of carbonate rock is influenced obviously by diagenesis.Dissolution is directly related to the quality of carbonate reservoir.The solubility of carbonateminerals in the buried environments is closely controlled by the chemistry of groundwater,the temperature,pH value and so on.According to the chemical thermodynamics theory,Gibbs free-energy changeΔG of the equilibrium reaction of calcite and dolomite can be calculated to judge the dissolution trend of calcite and dolomite.Thermodynamic models of the reactions between theminerals and groundwater were established in a closed system.This theoreticalmodel is used in the northern part of Tarim Basin to evaluate the dissolution trend of the Cambrian-Ordovician under the buried condition.The calculated results ofΔG from the reactions between carbonate rock and formation water are almost negative in the study area,which indicates that the buried environment is favorite for the dissolution of carbonate rock in the Cambrian-Ordovician.Gibbs free-energy of Cambrian dolomite is less than Ordovician limestone,and thatmeans dolomite is easier to dissolve than limestone in this deeply buried condition.This conclusion is also supported by cores and previous experimental simulation.Calculation results in Yuqiarea are also negative,which implies the strong precipitation of calcite in the caves in this region does not happen in present condition.This conclusion also can be proved from the cores of the Cambrian-Ordovician.

carbonateminerals;chemical thermodynamics;dissolution;Gibbs free energy

陳圓圓 女 1988年出生 碩士 儲層地質(zhì)學 E-mail:517371457@qq.com

P588.24+5

A

1000-0550(2012)02-0219-12

①國家973項目(編號:2011CB201100-03,2012CB214802)和國家油氣專項(編號:2011ZX05005-004-HZ06,2011ZX05009-002)聯(lián)合資助。

2011-01-20;收修改稿日期:2011-04-21