張路鋒，周福建，張士誠，汪杰，王晉

（中國石油大學(xué)（北京）非常規(guī)天然氣研究院·油氣資源與探測國家重點實驗室，北京102200）

0 引言

全球分布的致密氣可采儲量達114萬億立方米，國內(nèi)致密砂巖氣藏資源也極為豐富，可采儲量可達9～13萬億立方米，約占全國天然氣資源可采資源量的22%。致密砂巖氣藏是未來重要的油氣勘探區(qū)域，具有良好的發(fā)展前景。由于致密砂巖氣藏具有低孔、低滲、高毛管力以及低含水飽和度的儲層特性，鉆井過程中會引起鉆井液侵入到儲層，可能引發(fā)儲層中有效孔隙的堵塞以及各種敏感性傷害，從而導(dǎo)致儲層滲透率降低，進而降低儲層采收率，造成油氣藏開發(fā)效率低下[1-3]。塔里木克深區(qū)塊現(xiàn)場使用的鉆井液分為油基和水基兩類。與水基鉆井液相比，油基鉆井液具有抗高溫穩(wěn)定性強、抑制性強、潤滑性好、防塌效果好、抗污染能力強等優(yōu)點。但是，油基鉆井液也存在以下不足：油相進入地層使得有效流動面積減小，表面活性劑使油層發(fā)生潤濕反轉(zhuǎn)，降低油相滲透率、井眼凈化能力較差、配制成本較高等[4-6]。

目前，針對致密砂巖氣藏的鉆井液傷害評價研究還不成熟，國內(nèi)還是沿用常規(guī)儲層的評價方法，通過室內(nèi)測定直徑為2.54 cm、長度6 cm左右的小巖心滲透率在污染前后的變化來判斷儲層的傷害程度[7-8]。由于致密砂巖滲透率較低，常規(guī)巖心流動實驗驅(qū)替壓差大、穩(wěn)定時間長，并且實驗過程易受到隨時間引起的溫度變化以及測量系統(tǒng)自身流量不精確而造成較大誤差[9]。從而，許多學(xué)者使用劈裂的巖心模擬鉆井液對儲層裂縫的傷害，不考慮鉆井液對基質(zhì)的傷害[10-11]。但是Cipolla等[12]研究表明，當滲透率小于0.0001 mD，裂縫網(wǎng)絡(luò)對產(chǎn)能的極限貢獻在80%左右。也就是說，對于致密砂巖鉆井液對基質(zhì)的傷害不能忽略。因此，有必要開展鉆井液對基質(zhì)傷害的研究。

目前應(yīng)用最廣泛的超低滲巖心滲透率測量方法是非穩(wěn)態(tài)方法，包括壓力衰減法、壓力傳導(dǎo)法和周期震蕩法。壓力衰減法是由游利軍和康毅力等提出，并使用此方法評價了大牛地氣田盒3段儲層敏感性傷害[9]。周期振蕩法最早由Kranz等[13]借鑒熱擴散系數(shù)測量而提出，被運用到測量低滲巖石的滲透率上。Fischer等[14]進一步詳細闡述了周期振蕩法的理論背景，試驗設(shè)計以及數(shù)據(jù)處理。壓力傳導(dǎo)法是由 W.F .Brace等[15]在 1968 年首次提出，并用于測量Westerly花崗巖滲透系數(shù)，此外還給出了該方法的半解析解，但是其假設(shè)巖石孔隙度為零與實際巖心差別較大。隨后經(jīng)Hsieh和Dicker等人[16-18]修正改進，使得此方法在石油行業(yè)得以廣泛應(yīng)用。

考慮到研究區(qū)塊儲層滲透率較低，基于學(xué)者們對壓力傳導(dǎo)方法研究的基礎(chǔ)上，利用自主設(shè)計研發(fā)的壓力傳導(dǎo)儀，定量分析了塔里木克深區(qū)塊現(xiàn)場使用的油基、水基鉆井液濾液對基質(zhì)的傷害規(guī)律，建立了低滲儲層基質(zhì)傷害評價的新方法。

1 數(shù)學(xué)模型

1.1 模型建立

圖1是瞬態(tài)壓力傳導(dǎo)方法測量巖心滲透率的物理模型示意圖，初始時刻巖樣下方容器中的壓力為Po，實驗時巖樣上方始終以恒壓Pm驅(qū)替流體流動。顯而易見，該物理模型是沿垂直方向的一維滲流問題[19-23]，據(jù)此可以建立一個用于求取致密低滲巖樣滲透率的數(shù)學(xué)模型[24-26]。

圖1 壓力傳導(dǎo)滲透率測定示意圖

根據(jù)滲流力學(xué)建立模型的方法與步驟，可知，巖心內(nèi)一維擴散方程：

由于本實驗系統(tǒng)并不需要測量巖樣x=L處的流速（出口流速）。因此，通過壓縮系數(shù)對公式（5）進行修正：

1.2 模型求解

利用已有的邊界條件和初始條件，基于拉普拉斯變換，求解上述擴散方程。

首先，將擴散方程（1）變換到拉氏空間：

易知方程（7）為二階非齊次方程，因此，其通解為：

由于方程（8）中的系數(shù)C1和C2的值取決于邊界條件，為了求解系數(shù)C1和C2，必須將邊界條件變換成拉普拉斯形式。

邊界條件：

將邊界條件（9）帶入到方程（8），聯(lián)立求解，可以得到系數(shù)C1和C2：

因此，擴散方程（7）的通解可表示成式（11）：

針對方程（11），使用Carslaw等[24]給出的解，如下所示：

由于參數(shù)φn極大地取決于巖樣的孔隙體積與下游容器體積之比，當上述比值較小，方程在x=L時可簡化成如下方程：

式中，P（L，t）是下游壓力隨時間變化，其值等于巖樣在x=L處的壓力，即巖樣下端面的壓力，MPa與時間曲線的斜率；L為巖樣的高度，cm。

2 目標儲層

目標儲層位于中國西北部塔里木盆地，埋藏深度超過6000 m，是典型的超高溫、 高壓儲層。目標儲層厚度約為300 m，為分析其孔、 滲特征，鉆取273根巖心，進行室內(nèi)實驗（見表1）。巖心孔隙度分布在1.5%～7.6%，平均孔隙度為4.1%；巖心滲透率分布為小于 0.000 1 mD～2.905 mD，平均滲透率為 0.055 mD。

表1 目標儲層巖心孔隙度、滲透率分布頻率

3 室內(nèi)實驗

3.1 壓力傳導(dǎo)滲透率儀

壓力傳導(dǎo)滲透率儀是由中國石油大學(xué)（北京）儲層改造中心自主研發(fā)，由Core-Lab公司制作，其主要由 ISCO 泵（排量范圍為 0.000 01～60 mL/min）、儲液罐、三孔巖心釜（巖心夾持器）、恒溫箱以及數(shù)據(jù)采集傳輸系統(tǒng)等組成，見圖1。

圖2 壓力傳導(dǎo)滲透率儀（左）和三孔巖心釜（右）

3.2 實驗用巖樣

實驗中所使用的巖樣為塔里木克深儲層致密砂巖巖心，實驗巖樣信息如表2所示。

表2 實驗巖樣

實驗巖樣制備步驟如下：①使用線切割將從現(xiàn)場取回來的全直徑巖心，加工成3.6 cm×3.6 cm×15 cm的長方體巖心柱；②將加工好的巖心柱放入聚碳酸脂塑料管正中間，使巖心柱與塑料管的幾何中心一致，倒入環(huán)氧樹脂膠，靜止24 h；③待環(huán)氧樹脂膠干后，將整個體系放入110 ℃的恒溫箱，加熱1 h，確保巖心柱與環(huán)氧樹脂膠膠結(jié)為一體；④最后，使用線切割將其切為一個個0.65 cm厚的薄片。巖樣實物圖見圖3。

圖3 巖樣實物圖

3.3 實驗用鉆井液

塔里木克深致密砂巖儲層現(xiàn)場所用的鉆井液配方如下。

油基鉆井液 2%主乳化劑+3%輔乳化劑+2%有機土+2%氧化鈣+3%瀝青類降濾失劑+240 mL 0#柴油+ 60 mL氯化鈣溶液（20%）+加重材料

水基鉆井液 4%膨潤土漿+9%SMP-2+3%SPNH+1%PAC+0.5%AP-220+4%FT-1+3%碳酸鈣+2% 氫氧化鈉+加重材料

3.4 實驗步驟

①將煤油、鉆井液、標準鹽水裝入儲液罐中，巖樣放入三孔巖心釜中，連接好儀器，檢漏；②先上下游同時抽真空1 h，之后以一恒定壓力向下游注入煤油，壓力穩(wěn)定后記為下游初始壓力Po；③下游壓力穩(wěn)定后，以大于下游初始壓力的恒定壓力Pm向上游注入煤油，并以Pm的流動壓力持續(xù)流動；④監(jiān)測下游壓力隨時間的變化，即P（L，t）；⑤繪與t的關(guān)系曲線，并求取曲線的斜率ζ，利用公式（13）計算初始滲透率k1；⑥使用氮氣瓶恒壓驅(qū)替工作液，使工作液在巖樣上游持續(xù)流動3 h；⑦工作液傷害結(jié)束后，重復(fù)步驟②～⑤，可計算出傷害后的滲透率k2；⑧通過公式λd=（k1-k2）/k1×100%計算鉆井液對巖樣的傷害程度。

4 結(jié)果與討論

4.1 油基鉆井液傷害

裂縫性致密砂巖氣藏具有低孔低滲、裂縫發(fā)育、局部超低含水飽和度、高毛管壓力、地層壓力異常、高傷害潛力等特征，塔里木油田克深區(qū)塊具有以上所有特征，是典型的裂縫性致密砂巖氣藏。油基鉆井液在該區(qū)塊廣泛應(yīng)用，使用現(xiàn)場的油基鉆井液濾液開展了基質(zhì)傷害實驗。由于實驗巖心致密且黏土含量較高，容易受到水鎖、水敏等傷害。為了避免滲透率測量過程中，測試流體造成實驗巖心滲透率降低，上游流體選擇煤油，為了消除砂巖的滲析作用，下游流體也選擇煤油。圖4、圖5是實驗巖樣油基鉆井液傷害的壓力曲線，由于上下游存在滲透壓差以及水力作用，下游壓力會隨著時間逐漸上升，最終達到穩(wěn)定值，同時由于砂巖的半透膜效應(yīng)不明顯，下游壓力最終達到上游壓力水平。從圖4、圖5中可以看出，砂巖巖樣傷害前下游壓力上升速度較快，較短的時間壓力達到穩(wěn)定，傷害后下游壓力上升慢，下游壓力達到穩(wěn)定需要的時間較傷害前長，說明鉆井液對砂巖基質(zhì)存在一定的傷害。

圖6、圖7是砂巖實驗巖樣鉆井液傷害前后滲透率計算曲線。

圖4 油基鉆井液對KS-2#巖心傷害壓力傳導(dǎo)曲線

圖5 油基鉆井液對KS-2#巖心傷害壓力傳導(dǎo)曲線

圖6 油基鉆井液對KS-2#巖樣傷害滲透率計算曲線

圖7 油基鉆井液對KS-3#巖樣傷害滲透率計算曲線

通過擬合上下游實時壓差與初始壓差之比的對數(shù)值與時間的關(guān)系，可以得到斜率，利用公式（13）可以求出滲透率。易知KS-2#巖樣傷害前后滲透率分別為 0.65×10-3mD、0.45×10-3mD，傷害率30.77%；KS-3#巖樣傷害前后滲透率分別為1.025×10-3mD、0.7×10-3mD，傷害率為 31.71%。

4.2 水基鉆井液傷害

圖8、 圖9是克深砂巖巖樣受水基鉆井液傷害的壓力曲線，由于上下游存在滲透壓差以及水力作用，下游壓力會隨著時間逐漸上升，最終達到穩(wěn)定值。從圖8、圖9可以看出，與油基鉆井液傷害類似，克深致密砂巖傷害前下游壓力上升速度較快，較短的時間壓力達到穩(wěn)定，傷害后下游壓力上升慢，下游壓力達到穩(wěn)定需要的時間較傷害前長，說明水基鉆井液對砂巖基質(zhì)也存在一定的傷害。

圖8 水基鉆井液對KS-4#巖樣壓力傳導(dǎo)曲線

圖9 水基鉆井液對KS-5#巖樣壓力傳導(dǎo)曲線

圖10、圖11是克深砂巖實驗巖樣水基鉆井液傷害前后滲透率計算曲線。

圖10 水基鉆井液對KS-4#巖樣傷害滲透率計算曲線

圖11 水基鉆井液對KS-5#巖樣傷害滲透率計算曲線

由圖10、圖11可知，KS-4#巖樣傷害前后滲透率分別為 1.28×10-3、0.98×10-3mD，傷害率23.53%；KS-5#巖樣傷害前后滲透率分別為0.63×10-3mD、0.48×10-3mD，傷害率 23.81%。

5 結(jié)論

1.克深致密砂巖油基鉆井液平均傷害率為31.24%，水基鉆井液平均傷害率為23.67%。水基鉆井液對基質(zhì)傷害程度略低于油基鉆井液，但是現(xiàn)場使用什么類型的鉆井液需要綜合考慮。

2.實驗中上下游壓力最終趨于一致，說明砂巖膜效應(yīng)不明顯。

3.僅僅評價了油基/水基鉆井液對儲層基質(zhì)的傷害，但瞬態(tài)壓力傳導(dǎo)方法并不局限與此。該方法還可以用來評價致密低滲儲層敏感性傷害。