龐 銘，陳華興，唐洪明，盧 浩，王 昭

(1.中海石油(中國)有限公司天津分公司，天津 300450；2.西南石油大學(xué)，四川 成都 610500)

0 引 言

中國碳酸鹽巖油氣藏油氣資源極其豐富，占全國油氣資源總量的40%以上[1]。應(yīng)力敏感性損害是碳酸鹽巖儲層的主要的儲層損害方式之一[2-4]。前人針對有效應(yīng)力、儲層類型、裂縫充填程度與充填類型對應(yīng)力敏感性損害的影響已進行大量的巖心實驗、有限元數(shù)值模擬及離散元數(shù)值模擬。研究普遍認為：溶蝕孔隙、溶蝕孔洞的發(fā)育弱化了裂縫型碳酸鹽巖儲層的應(yīng)力敏感性損害[5-6]；裂縫充填程度越大，應(yīng)力敏感性損害越弱[7]；縫面接觸關(guān)系是影響應(yīng)力敏感性損害的關(guān)鍵因素[8-11]。碳酸鹽巖儲層通常需要進行酸化解除固相堵塞以提高油氣生產(chǎn)效率[12]，甚至需要進行酸化壓裂產(chǎn)生人工裂縫以溝通儲集體[13-14]。油氣開采過程中，不同有效應(yīng)力下的酸蝕裂縫導(dǎo)流能力(長期導(dǎo)流能力)是油藏開采效率的重要指標，酸蝕前后裂縫型碳酸鹽巖儲層的應(yīng)力敏感性損害存在差異，而目前針對酸蝕對裂縫型碳酸鹽巖應(yīng)力敏感性損害的影響研究基本處于空白。由于應(yīng)力敏感性實驗流程(SY/T 5338—2010)對巖心具有不可逆破壞作用，目前尚未建立一套針對酸蝕前后的裂縫型碳酸鹽巖儲層應(yīng)力敏感性損害系統(tǒng)評價方法。因此，選用LH油田奧陶系裂縫型碳酸鹽巖儲層巖心，利用裂縫三維特征掃描獲取縫面微凸體及縫寬的分布、分形特征，在裂縫重構(gòu)的基礎(chǔ)上，利用有限元方法模擬裂縫閉合過程，對不同閉合階段的裂縫流動空間特征和滲透率變化進行量化分析，明確酸蝕前后裂縫的應(yīng)力敏感性變化，對指導(dǎo)裂縫型碳酸鹽巖儲層高效開發(fā)具有現(xiàn)實意義。

1 實驗方法

1.1 實驗材料

實驗樣品取自LH油田奧陶系，屬于裂縫-孔隙(孔洞)型儲層，巖心孔隙度為0.11%～8.97%，平均為1.38%；滲透率為0.002～320.254 mD，平均為1.390 mD。裂縫受構(gòu)造過程控制并經(jīng)過成巖作用改造。該區(qū)裂縫以構(gòu)造縫為主，裂縫多呈高角度縫特征，主要為微縫(縫寬小于0.10 mm)和小縫(縫寬為0.10～1.00 mm)；局部地區(qū)發(fā)育大縫、中縫，占總裂縫的17.2%，是重要的油氣滲流通道。全充填縫約占總裂縫的30.0%，半充填縫約占總裂縫的10.0%，非充填縫占總裂縫的60.0%。將巖心樣品制成長度為50 mm，直徑為100 mm的圓柱狀，裂縫內(nèi)未充填。

實驗所使用的光柵掃描系統(tǒng)為江蘇先臨公司生產(chǎn)的shinning 3D系列，掃描精度為0.01 mm，掃描時CCD相機分辨率可達1 280×1 024像素，10 s內(nèi)的一次測量可以得到數(shù)量為130×104的點云集合，能夠達到表征微凸體特征的精度。

1.2 實驗過程

酸蝕裂縫滲流能力受裂縫有效應(yīng)力、縫寬、滲透率等因素決定，縫寬的變化也會引起滲透率的變化，即縫寬對滲流能力的影響大于雙滲透率的影響。首先對比酸蝕前后裂縫縫寬的分布特征，然后施加有效應(yīng)力，得到應(yīng)力作用后的縫寬特征，最后評價不同有效應(yīng)力下的滲流特性。具體數(shù)值模擬流程：①使用3D光柵掃描系統(tǒng)測量酸蝕前裂縫的表面形貌，然后用10%的鹽酸對裂縫表面進行均勻刻蝕(酸化反應(yīng)時間2 h)，再次測量酸蝕后裂縫的表面形貌。②將酸蝕前后裂縫的點云數(shù)據(jù)進行逆向重構(gòu)并實體化，將實體化裂縫導(dǎo)入有限元軟件中并網(wǎng)格化處理，模擬裂縫閉合過程。③輸出每個閉合步長的裂縫流動空間，計算不同閉合程度下的裂縫滲透率。由于對每個裂縫面進行單獨掃描，為保證縫面耦合時的準確性，需要使縫面單獨掃描與樣品整體掃描時處于同一坐標系下。采用特征掃描技術(shù)可有效解決該問題，將樣品外部添加的標志點及表面已掃描點云數(shù)據(jù)作為特征點，先將樣品整體掃描，再單獨掃描2個縫面，掃描時進行特征識別，通過確定的特征點空間相對位置關(guān)系，重新建立與巖心樣品掃描時相同的局部坐標系[15]。

裂縫閉合有限元數(shù)值模擬采用ABAQUS6.4軟件進行，裂縫重構(gòu)模型網(wǎng)格劃分采用修正的二次四面體網(wǎng)格單元(C3D10M)。將裂縫的下表面巖塊設(shè)置為法向約束，上表面巖塊施加位移邊界條件，每個步長的位移量為0.05 mm，在閉合過程中縫面間施加硬接觸，進而分析裂縫閉合過程中裂縫面上微凸體的應(yīng)力應(yīng)變情況。在地層條件下，裂縫側(cè)面的變形特征與其周圍的巖石變形應(yīng)相同，因此，在側(cè)面施加對稱邊界條件，對稱邊界條件同時能夠約束裂縫巖塊受轉(zhuǎn)動產(chǎn)生的剛性位移，有利于計算收斂。根據(jù)巖心的實測巖石力學(xué)參數(shù)設(shè)置其他條件：彈性模量為30.6 GPa，泊松比為0.32，內(nèi)聚力為23.78 MPa，內(nèi)摩擦角為58.6 °。在每一步長后，提取裂縫上部巖塊頂表面上的平均應(yīng)力作為每個時間步的法向應(yīng)力，并輸出裂縫流動空間，分析裂縫滲透率變化。

將每個時間步的裂縫流動空間導(dǎo)入計算流體動力學(xué)軟件中，并在水平方向施加壓力梯度，模擬不同方向的滲透率變化。在滲流模擬過程中假設(shè)：①由于裂縫滲透率遠大于基質(zhì)滲透率，忽略基質(zhì)滲透率和基質(zhì)-裂縫的物質(zhì)交換；②流體為不可壓縮牛頓流體；③在雷諾數(shù)較小的情況下，忽略慣性效應(yīng)；④裂縫中水平方向上的速度梯度遠小于法向上的速度梯度。流體控制方程為[16-17]：

(1)

式中：p為裂縫中的流體壓力，MPa；w為局部裂縫寬度，m；x、y分別為裂縫縫面橫、縱坐標。

2 結(jié)果分析

2.1 酸蝕前后裂縫特征

研究裂縫時，通常將裂縫面假設(shè)為光滑的平面，裂縫寬度固定或漸變。然而在實際儲層中，裂縫依靠縫面上粗糙、不規(guī)則的微凸體自支撐從而保留一定的滲流空間，縫面微凸體之間的耦合規(guī)律復(fù)雜，極大地影響對裂縫閉合規(guī)律的研究。為定量研究儲層裂縫空間特性及其閉合規(guī)律，需要對縫面微凸體及其耦合特征進行測量及描述，而評價縫面微凸體分布的重要參數(shù)為裂縫粗糙度系數(shù)[18-24]。表1為裂縫縫面酸蝕前后的粗糙度。由表1可知，酸蝕后裂縫粗糙度明顯增加，這是由于礦物組分差異、儲集空間發(fā)育不均引起的非均勻酸蝕造成的。

表1 酸蝕前后裂縫縫面粗糙度對比

裂縫由一對粗糙裂縫縫面耦合而成，影響閉合規(guī)律的關(guān)鍵因素不是單一縫面微凸體分布特征，而是2個縫面合成之后的物理結(jié)構(gòu)。合成的裂縫表面反映了2個粗糙縫面之間的相關(guān)程度，體現(xiàn)了2個粗糙縫面微凸體之間的接觸特性，可以用來分析裂縫受力、微凸體接觸變形特征?？p面合成方法采用Brown & Scholz合成表面原理，將裂縫2個縫面的接觸過程簡化為合成表面與基準面之間的接觸關(guān)系[17]。裂縫縫寬可由下式計算：

w(x，y)=z1(x，y)-z2(x，y)-min[z1(x，y)-z2(x，y)]

(2)

式中：z1和z2分別為裂縫上下表面相對于基準面的高度，mm。

裂縫合成后的縫寬分布特征如圖1所示，縫寬基本參數(shù)如表2所示，其中，分形維數(shù)采用三菱柱法進行計算。由表2可知：酸蝕后平均縫寬增加了近0.5 mm，最大縫寬增加了近2.0 mm，均方根值增加了0.2 mm，分形維數(shù)增加0.06，表明酸蝕后縫寬分布非均質(zhì)特征更加嚴重，裂縫面間耦合性變差。對比酸蝕前后縫寬(圖1)可知，酸蝕導(dǎo)致左上角產(chǎn)生了一條明顯的溝槽，右上角產(chǎn)生多個較大的酸蝕孔洞。通過對碳酸鹽巖裂縫酸蝕形貌特征的研究進行調(diào)研發(fā)現(xiàn)，由于難溶礦物的支撐、酸巖反應(yīng)不均一等因素，碳酸鹽巖裂縫酸蝕后普遍存在溝槽結(jié)構(gòu)，而有效應(yīng)力下裂縫的機械變形和水力性質(zhì)受縫面影響，其應(yīng)力敏感性損害存在較大差異。

圖1 酸蝕前后裂縫寬度分布特征

表2 酸蝕前后裂縫寬度參數(shù)

2.2 裂縫閉合過程中的變形機理

在地應(yīng)力作用下，微凸體之間的不完全耦合導(dǎo)致裂縫不能完全閉合，保持一定的滲流空間。在裂縫受到有效應(yīng)力而閉合的過程中，縫面上突起較高的微凸體(主要微凸體)首先支撐并抗拒閉合。隨有效應(yīng)力的繼續(xù)增大，主要微凸體發(fā)生壓縮形變，縫面間的距離變小，突起較低的微凸體(次級微凸體)開始起到支撐作用，進一步增加了裂縫面的抗閉合能力。

圖2為酸蝕前后的裂縫表面Mises應(yīng)力隨著有效應(yīng)力增加的變化，圖3展示了裂縫接觸面積和平均裂縫寬度隨有效應(yīng)力的變化。由圖2、3可知：隨有效應(yīng)力增加，縫面微凸體接觸面積明顯增加，裂縫逐漸閉合，裂縫縫寬減小。當(dāng)有效應(yīng)力小于10 MPa時，裂縫接觸面積較小，裂縫抗閉合能力較弱，裂縫快速閉合，裂縫縫寬快速減小。該階段酸蝕前裂縫的接觸面積略大，占比為2.82%，而酸蝕后裂縫的接觸面積較小，占比僅為1.03%，這是因為酸蝕導(dǎo)致裂縫耦合變差，進一步影響裂縫的抗閉合能力。當(dāng)有效應(yīng)力為5 MPa時，酸蝕前大部分裂縫均已閉合，而酸蝕后的裂縫仍然具有酸蝕形成的溝槽，保留了較好的滲流能力。當(dāng)有效應(yīng)力超過10 MPa時，縫面微凸體接觸面積明顯增加，此時次級微凸體逐漸開始對裂縫產(chǎn)生支撐作用，裂縫閉合速度減緩。由圖2c、d可知，該階段部分微凸體所承受的Mises應(yīng)力大于1 000 MPa，超過了一般巖石的彈性形變極限。這表明當(dāng)巖石裂縫發(fā)生閉合變形后，主要微凸體已經(jīng)發(fā)生不可恢復(fù)的損傷，即使儲層的壓力恢復(fù)，裂縫的寬度和流動能力也難以恢復(fù)。

圖2 酸蝕前后裂縫在閉合過程中的應(yīng)力分布對比

圖3 酸蝕前后裂縫的閉合規(guī)律

2.3 應(yīng)力敏感性損害程度評價

酸蝕導(dǎo)致的裂縫耦合度變化、有效應(yīng)力導(dǎo)致的裂縫閉合均會造成裂縫中滲流場發(fā)生改變。因此，對不同有效應(yīng)力下酸蝕前后裂縫的滲流場進行了計算(圖4)，進而利用數(shù)值模擬結(jié)果繪制酸蝕前后有效應(yīng)力與滲透率保留率的關(guān)系曲線(圖5)。由圖4可知：當(dāng)有效應(yīng)力為5 MPa時，酸蝕前后的裂縫滲流速度均有所減緩，但存在較大差異，酸蝕前的裂縫幾乎不可見優(yōu)勢滲流通道，而酸蝕后的裂縫依然具有優(yōu)勢滲流通道(即酸蝕形成的溝槽)，該通道可以有效減緩應(yīng)力敏感性損害。由圖5可知：當(dāng)有效應(yīng)力增至40 MPa時，酸蝕前裂縫滲透率損害率大于95%，而酸蝕后裂縫滲透率損害率小于75%，應(yīng)力敏感性損害減弱。這是因為裂縫經(jīng)過酸蝕后，縫面粗糙程度增加，縫面間耦合程度降低，形成了明顯的溝槽。隨著法向有效應(yīng)力的增加，酸蝕后的裂縫閉合量和縫面接觸率均較酸蝕前的裂縫小，造成滲透率保留率仍然相對較大，應(yīng)力敏感性損害程度相對較低。裂縫應(yīng)力敏感性損害受縫面接觸關(guān)系、巖石力學(xué)性質(zhì)和滲流能力影響，其中，縫面接觸關(guān)系是基礎(chǔ)，巖石力學(xué)性質(zhì)是直接誘因，滲流能力降低是最終表現(xiàn)。在碳酸鹽巖油氣藏開發(fā)過程中，對井筒進行酸洗、對儲層進行酸化壓裂，不僅有助于增加儲層人工裂縫滲透性，還能增加縫面間的非耦合性，降低生產(chǎn)過程中的應(yīng)力敏感性損害，防止裂縫閉合。

圖4 不同法向應(yīng)力下的酸蝕前后裂縫滲流差異

圖5 裂縫閉合過程中的滲流規(guī)律變化

3 結(jié)果驗證

為驗證數(shù)值模擬方法研究結(jié)果的準確性，采用石油天然行業(yè)標準SY/T 5358—2010中的內(nèi)壓變換方式，對研究區(qū)裂縫型碳酸鹽巖儲層相鄰的4塊巖心進行酸蝕應(yīng)力敏感性損害實驗。巖心基本參數(shù)如表3所示，實驗結(jié)果見圖6。

表3 巖心基本參數(shù)

圖6 應(yīng)力敏感性損害實驗中裂縫巖心滲透率變化

由圖6可知：隨著地層壓力的降低，滲透率逐漸減??；當(dāng)?shù)貙訅毫?8 MPa時，有效應(yīng)力增大，滲透率降低幅度逐漸變小，即滲透率隨著地層壓力的降低逐漸穩(wěn)定。當(dāng)?shù)貙訅毫抵?0 MPa時，酸蝕后的裂縫應(yīng)力敏感性損害為50%～70%，損害程度中等偏強；未酸蝕的裂縫滲透率損害率為76%～99%，損害程度強。與未酸蝕的裂縫相比，酸蝕后裂縫應(yīng)力敏感性損害變?nèi)?，這是因為酸蝕不均勻形成的蚓洞等優(yōu)勢滲流通道或非酸溶礦物(如半充填自生石英)對裂縫起到了一定的支撐作用，防止了裂縫完全閉合。實驗結(jié)果與數(shù)值模擬結(jié)果基本一致，說明該方法準確，結(jié)果可靠。

4 結(jié) 論

(1) 提出了一種基于裂縫重構(gòu)評價酸蝕對碳酸鹽巖應(yīng)力敏感損害影響的數(shù)值模擬方法，從縫面幾何形貌、縫面接觸關(guān)系、縫面接觸應(yīng)力、裂縫滲流特征多個指標對應(yīng)力敏感性損害進行評價。

(2) 當(dāng)有效應(yīng)力小于10 MPa時，裂縫寬度快速下降，滲透率急劇損失；當(dāng)有效應(yīng)力大于10 MPa后，裂縫寬度下降和滲透率損害程度逐漸平緩。當(dāng)有效應(yīng)力為40 MPa時，酸蝕前裂縫縫寬減小量大于55%，滲透率損害率大于90%；而酸蝕后縫寬減小量小于40%，滲透率損害率小于75%，應(yīng)力敏感性損害降低。

(3) 數(shù)值模擬及實驗結(jié)果表明，裂縫型儲層中縫面礦物組分存在差異，酸液對縫面產(chǎn)生非均勻酸蝕，酸蝕后縫面粗糙度增加，易產(chǎn)生溝槽，進而降低裂縫應(yīng)力敏感性損害程度?，F(xiàn)場生產(chǎn)實踐中，酸化壓裂、酸洗不僅有助于提高儲層滲透性，還能在一定程度上防止裂縫閉合。