劉順喜,樊坤雨,金 毅,于振鋒,董佳斌,王長征
(1.河南理工大學 資源環(huán)境學院,河南 焦作 454003;2.煤與煤層氣共采國家重點實驗室,山西 晉城 048012;3.煤炭安全生產與清潔高效利用省部共建協(xié)同創(chuàng)新中心,河南 焦作 454003)
當前我國煤層氣勘探開發(fā)正逐漸由淺層向深層發(fā)展,僅沁水盆地和鄂爾多斯盆地埋深1 000~3 000 m 的煤層氣資源量就達到33 萬億m3,勘探開發(fā)前景十分可觀[1-3]。深部煤層氣資源由于埋深大、蓋層厚,保存條件較好,其含氣量往往較高。但是埋深較高的同時也存在地應力高、滲透率低的不利因素,而滲透率是決定煤層氣開采成敗的關鍵參數之一,其高低直接影響煤層氣井產能和開發(fā)效益[4-6]。因此,探討滲透率的關鍵影響因素成為深部煤層氣資源勘探開發(fā)的熱點問題。
眾多學者從不同角度開展了滲透率的影響因素研究,認為應力是影響煤層滲透性的關鍵因素[7-10]。排水降壓初期,儲層壓力逐漸降低(導致有效應力增加),裂縫孔隙閉合,煤儲層滲透率隨著有效應力的增加而降低。M.A.Biot[11]1941 年首次提出滲透率隨有效應力呈指數下降。隨后,Zeng Kaihua 等[12]討論和分析了煤在三軸壓縮條件下的滲透率變化,認為裂隙的演化是影響煤滲透率變化的重要因素;陳世達[13]、黃強[14]等探究有效應力對中高階煤儲層滲透率的控制作用及其應力敏感性;馬如英等[15]研究了低階煤儲層孔隙率和滲透率特征及儲層應力敏感性對煤層氣開采的影響。程鳴[16]、孟雅[17]、孟召平[18]等開展煤系“三氣”覆壓孔滲對比實驗研究,分析不同煤系的應力敏感性特征,得出煤層應力敏感性隨不同有效應力階段而發(fā)生動態(tài)變化。Li Jianhua[19]、王攀[20]等基于煤體的應力應變關系及煤體裂縫與基質的相互作用,構建了考慮不同應力條件的各向異性滲透率模型。同時,也有部分學者開展了有效應力與不同方向煤巖滲透率的關系研究[21-24],結果表明,割理中部分充填了方解石和黏土礦物,水平層理面方向的滲透率比垂直層理面方向的滲透率大約高2.5 倍,但上述主要通過物理模擬實驗測試溫度、圍壓、飽和度及有效應力等對滲透率的總體影響趨勢,并未詳細闡述各向異性隨有效應力增大應力敏感性的變化和相應階段孔裂隙變化特征,且研究對象主要為中、淺部煤層。
因此,筆者選取深部煤層氣勘探開發(fā)的熱點區(qū)域-沁水盆地東北部橫嶺區(qū)塊,以15 號煤層為研究對象,通過覆壓孔滲試驗,對目標煤層在有效應力增加過程中,應力敏感性變化趨勢開展階段劃分,并分析平行層理面和垂直層理面方向上滲透率以及應力敏感性的各向異性特征,以期揭示深部煤儲層應力敏感性變化規(guī)律,為研究區(qū)煤層氣高效勘探開發(fā)提供理論支撐。
沁水盆地為中生代以來形成的構造盆地,先后歷經海西、印支、燕山和喜馬拉雅4 期構造運動疊加影響,盆地主體為一軸向呈NNE 展布的大型寬緩復向斜。盆地構造應力場以東西向擠壓應力為主,構造變形強度由盆緣向內逐漸減弱,邊緣多發(fā)育斷裂構造。橫嶺區(qū)塊位于沁水盆地東北部,沁水復向斜東翼沾尚-武鄉(xiāng)-陽城褶皺帶內,區(qū)內以正斷層為主,傾角大于70°,褶皺兩翼寬緩,傾角小于6°,構造線方向受區(qū)域構造應力控制,呈NNE 向(圖1)。樣品采自橫嶺區(qū)塊3 口煤層氣井的15 號煤。15 號煤形成于晚石炭-早二疊世太原組潮坪環(huán)境中,為區(qū)內穩(wěn)定可采煤層,埋深大于1 200 m,煤厚2.02~9.90 m,平均5.20 m,含1~3 層夾矸,煤階為貧煤,宏觀煤巖類型為半亮煤。
圖1 沁水盆地橫嶺區(qū)塊構造及采樣位置Fig.1 General tectonic characteristics of the study region and sampling location
為研究有效應力對滲透率各向異性的影響,分別對3 塊煤樣沿平行層理和垂直層理方向鉆取直徑為25 mm 的柱狀樣品(圖2)。平行層理面樣品編號為1-H、2-H 和3-H,垂直層理面樣品編號為1-V、2-V 和3-V,為保證樣品原始孔裂隙結構,未對柱狀樣進行特殊的烘干處理(表1)。
表1 煤樣基本特征Table 1 Basic properties of coal samples
圖2 煤樣制備Fig.2 Schematic diagram of coal sample preparation
實驗采用覆壓孔隙率-滲透率測試儀AP-608,非穩(wěn)態(tài)脈沖衰減法,執(zhí)行中國石油天然氣行業(yè)標準SY/T 6385-2016《覆壓下巖石孔隙度和滲透率測定方法》,以He 為測試介質。實驗過程中保持進口壓力不變,利用平流泵逐步增加樣品的圍壓值,依據研究區(qū)煤層埋深情況共設計7 個壓力點(3.5~15.0 MPa),壓力點平衡時間控制在30 min 以上,分別測定各壓力點的滲透率。
滲透率是評價煤層氣可采性的關鍵參數之一,反映煤儲層的滲流能力,是非常規(guī)天然氣資源開發(fā)的關鍵影響因素。對實驗數據進行回歸分析,獲得煤樣的滲透率隨有效應力的增大呈冪指數函數規(guī)律衰減,滿足下式:
式中:ke為 已知應力條件下的滲透率,10-3μm2;k0為初始有效應力為0 時的滲透率,10-3μm2;p為有效應力,MPa;c為衰減系數,MPa-1。
初始滲透率k0是煤樣脫離原位儲層,不受任何外部應力和不考慮基質吸附變形時的滲透率,無法通過實驗測得,可以通過式(1)進行回歸分析擬合得出。6 個樣品的初始滲透率如圖3 所示(相關系數R2>0.960 4)。由圖中可知,平行層理面初始滲透率大小為1-H>3-H>2-H,分別為2.778 5×10-3、0.214 4×10-3和0.111 8×10-3μm2;垂直層理面初始滲透率3-V>2-V>1-V,分別為1.332 3×10-3、0.085 2×10-3和0.052 3×10-3μm2。因此依據現有數據分析,研究區(qū)1 號井位附近平行層理面滲透率具有優(yōu)勢,而3 號井位附近垂直層理面滲透率具有優(yōu)勢。對比分析同一樣品平行層理面和垂直層理面滲透率發(fā)現,除3 號樣品外,其余2 個樣品平行層理面滲透率均大于垂直層理面滲透率。表明沉積物在沉積過程中形成的層理面是滲透率大小的優(yōu)勢方向。
圖3 滲透率與有效應力之間的關系Fig.3 Relationship between permeability and effective stress
進一步分析平行層理面和垂直層理面滲透率隨有效應力變化規(guī)律,并進行線性擬合,二者呈正相關性(R2>0.974 6,圖4)。從圖中可以看出,在不同方向上滲透率隨有效應力的增加,呈一定比例系數增大。樣品1-樣品3 平行層理面與垂直層理面滲透率比例系數分別為0.01,0.77 和3.90。從比例系數來看,由于煤樣個體存在差異性導致變化系數差別較大,并無特定比例關系,但比例系數提供了預測不同方向上滲透率的方法,對煤層氣勘探開發(fā)具有較強的現實指導意義。
圖4 平行層理面與垂直層理面滲透率變化關系Fig.4 Permeability variation relationship between parallel bedding plane and vertical bedding plane
煤層氣開采過程中隨著孔隙流體的排出,有效應力不斷增加,使得儲層巖體發(fā)生變形,產生應力敏感現象,導致滲透率降低。本文選取儲層孔裂隙壓縮系數、滲透率損害率和滲透率曲率3 個參數開展深部煤儲層應力敏感性階段劃分和各向異性特征研究。
2.3.1 孔裂隙壓縮系數
煤儲層孔裂隙隨有效應力的增加逐漸被壓縮,導致滲透率急劇下降,其變化范圍可達2~3 個數量級??琢严秹嚎s系數是定量表征滲透率隨有效應力變化的程度,是最直觀反映儲層應力敏感性參數。C.R.Mckee等[25](1988)給出了滲透率、有效應力與孔裂隙壓縮系數的關系式,即:
式中:a為孔裂隙壓縮系數,MPa-1;p0為初始有效應力,MPa。
由式(2)可推導出孔裂隙壓縮系數:
由式(3)可知,a是參數p-p0和 -[ln(ke/k0)]/3的線性函數系數,以二者為橫坐標和縱坐標進行線性擬合,直線的斜率即為孔裂隙壓縮系數(圖5)。
孔裂隙壓縮系數的差異導致應力敏感性具有各向異性。由圖5 可知,垂直層理面孔裂隙壓縮系數為0.073 1~0.100 9 MPa-1,平均0.082 6 MPa-1;平行層理面孔裂隙壓縮系數為0.068 4~0.074 9 MPa-1,平均0.070 9 MPa-1。無論是單個樣品結果還是平均值,在相同壓差的條件下,垂直層理面壓縮系數均大于平行層理面,表明儲層在垂直層理面方向上具有更高的應力敏感性。研究區(qū)15 號煤層兩個方向上的應力敏感性比例平均值為1.17∶1。
圖5 不同方向煤樣孔裂隙壓縮系數計算曲線Fig.5 Equating curves of the pore-fracture compression coefficient of coal samples in different directions
進一步分析孔裂隙壓縮系數與有效應力差值之間關系,發(fā)現壓縮系數以壓差5 MPa 為界限分為2段(圖6),前段呈正相關性,后段則表現出顯著的負相關性。從圖中可知,孔裂隙壓縮系數為一動態(tài)參數,并非定值,在有效應力增加初期,由于煤中宏觀裂隙具有大量可壓縮的空間,其壓縮系數逐漸增大;之后,隨著壓差不斷增加,儲層轉入顯微孔裂隙壓縮階段,可壓縮空間不斷降低,壓縮難度越來越大,因此壓縮系數呈現降低趨勢。
圖6 孔裂隙壓縮系數與壓差變化趨勢Fig.6 Change trend of pore-fracture compression coefficient and pressure difference
2.3.2 滲透率損害率
滲透率損害率是反映儲層應力敏感性的又一重要參數,可表征滲透率在有效應力增加情況下樣品受損傷的程度,用百分數表示。該值越大,表明儲層滲透率受有效應力損害程度越嚴重,對應的應力敏感性越高,其表達式為:
式中:Is為滲透率損害率,%;I1為第1 個應力點下的滲透率,10-3μm2;Ii為第i個應力點下的滲透率,10-3μm2。
由圖7 可知,滲透率損害率與有效應力之間呈指數關系(R2>0.985 5),有3 個變化區(qū)間:當有效應力達到5 MPa 時,滲透率損害率為68.38%~90.06%,平均76.13%,滲透率損害率近似呈線性關系快速增加;當有效應力達到12 MPa 時,滲透率損害率為91.81%~98.87%,平均94.99%,樣品的滲透率損害率增加程度放緩,與上一階段平均值差為19%,整體程度低于第一階段;當壓力達到15 MPa 時,滲透率損害率達到最大值,介于95.01%~99.43%,平均96.91%,滲透率降低程度有限,與上一階段平均值差僅為2%。
圖7 滲透率損害率與有效應力關系Fig.7 Relationship between permeability stress damage rates and effective stress
2.3.3 滲透率曲率
在數學上通常用曲率來描述物體的彎曲程度,曲率越大,表明物體在該點處彎曲變形程度越高。滲透率曲率可以更有效地表征滲透率隨有效應力增大而降低的速率,即滲透率曲率越大,滲透率的應力敏感性就越強;反之越弱[22]。因此,本次通過引入滲透率曲率來對比分析滲透率受應力的影響程度,進一步揭示滲透率對應力的敏感程度。
式中:kc為滲透率曲率,% ;k′、k′′分別為滲透率對有效應力的一階、二階導數。
滲透率曲率為:
式中:c為應力敏感系數,MPa-1。
由圖8 可知,6 個樣品滲透率曲率隨有效應力的增加呈負指數規(guī)律降低(R2>0.999 9),當有效應力小于5 MPa 時,滲透率曲率隨有效應力的增加快速下降,表明該階段儲層中孔裂隙的壓縮速率下降極快,滲透率快速降低;當有效應力為5~12 MPa 時,隨有效應力的增大,滲透率曲率下降趨勢逐漸變緩,表明孔裂隙的可壓縮空間已經降低,滲透率下降程度變緩;當有效應力大于12 MPa 時,滲透率變化程度微弱。
圖8 滲透率曲率與有效應力關系Fig.8 Relationship between permeability curvature and effective stress
綜上所述,同時結合滲透率損害率和滲透率曲率的變化程度,可將研究區(qū)樣品相應劃分為宏觀裂隙壓縮、顯微裂隙壓縮和孔隙壓縮3 個階段(圖7)。
平行層理面和垂直層理面上,其滲透率曲率具有顯著差異(表2):其中,平行層理面樣品其滲透率曲率大小為3-H>2-H>1-H,而垂直層理面其滲透率曲率大小為1-V>3-V>2-V,由此可知,3 號樣品在平行層理方向上、1 號樣品在垂直層理面方向上具有較高的應力敏感性。
表2 儲層應力敏感性評價參數Table 2 Stress sensitivity evaluation parameters of coal reservoirs
進一步對每個樣品在兩個方向上的滲透率曲率開展分析發(fā)現,研究區(qū)垂直層理面滲透率曲率均大于平行層理面方向,表明儲層的滲透率在垂向上受有效應力影響更大,具有更強的應力敏感性。
煤層氣排采經歷解吸-擴散-滲流3 個階段,包括排水、降壓和產氣過程。假定儲層水平等厚,煤層氣以達西平面徑向流的方式流向井底,考慮應力敏感性但不考慮表皮系數的影響,標準條件下氣井的產能公式[26]:
式中:Qsc為標準狀態(tài)下的產氣量,m3/d;h為煤層的有效厚度,m;μ為氣體黏度,MPa·s;Z為壓縮因子;T為煤儲層溫度,K;pwf為井底流壓,MPa;pe為原始儲層壓力,MPa;Le為泄壓半徑,m;Lw為井筒直徑,m;其中,標準狀態(tài)下的壓力、溫度和偏差系數依次?。?.101 328 MPa、293 K 和1。
當完全不考慮煤儲層應力敏感性的影響時,則氣體的滲流是滲透率為k0條件下的達西滲流[26],其氣井產氣量qsc為:
為分析煤儲層應力敏感性對本研究區(qū)產氣井的影響,對采樣點3 口井在不同生產壓差、分別考慮和不考慮應力敏感性條件下,通過式(7)和式(8)進行產能理論對比分析,煤層氣井相關參數見表3。由圖9 可知,煤儲層不考慮應力敏感性的氣井產量要高于考慮應力敏感性的氣井產量,且隨著生產壓差不斷增大,受到的有效應力逐漸增加,煤儲層滲透率指數降低,影響其氣井產能。應力敏感性對氣井產能有顯著影響,煤層氣排采過程中,隨著生產壓力的增加,氣井的產量增加幅度變小,并逐漸趨向穩(wěn)定,說明單純依靠增大生產壓差并不能獲得最大產量[25]。
表3 煤層氣井參數Table 3 Parameters of CBM wells
圖9 煤層氣井產能隨生產壓差的變化Fig.9 Changes of CBM well productivity with production pressure difference
為了反映應力敏感性對煤層氣產量的影響程度,可定義產量降低幅度為δ:
取井位1、2、3 的6 個樣品裂隙壓縮系數作為應力敏感性系數。為便于對比分析,統(tǒng)一生產壓差為5 MPa,表征應力敏感性對氣井產能的影響,由式(9)定義生產壓差為5 MPa 的產量降低幅度為δ5。
計算結果如圖10 所示,在同一生產壓差下可以看出,氣井的產量降低幅度δ5隨應力敏感性系數a的增大整體呈增高趨勢。究其原因,排水降壓階段隨著孔隙流體不斷采出,煤儲層的有效應力逐漸增大,造成孔裂隙開度降低甚至發(fā)生閉合,而排采過程中的壓力波動會造成孔裂隙的永久縮小,導致其滲透率急劇降低及滲透率損害率逐漸增大,進而影響煤層氣產量。
圖10 產能降低幅度與應力敏感性系數的關系(壓力差5 MPa)Fig.10 Relationship between stress sensitivity coefficient and productivity reduction rate (p-p0=5 MPa)
有效應力增加的實質是由于排水降壓導致儲層壓力降低引起,從而使?jié)B透率發(fā)生動態(tài)變化。因此,滲透率與有效應力之間關系反映了煤層氣井排采過程中儲層滲透率與生產壓差之間的關系[27]。生產壓差是儲層壓力與井底流壓之間的差值,因而煤層氣井排采需要制定合理的動液面高度(生產流量)來改變生產壓差。同時既要保證較高產量,又要避免生產壓差過大導致煤儲層產生明顯的應力敏感性,減少其對煤儲層傷害程度,確保氣井的高產和穩(wěn)產。因此,通過對研究區(qū)儲層應力敏感性的階段劃分,可以把握關鍵壓力點,更有效控制煤層氣的產出過程。
結合研究區(qū)3 口煤層氣井的應力敏感性特征,給出以下建議:
(1) 當有效應力小于5 MPa 時,滲透率隨有效應力的增加迅速下降,應力敏感性最強,橫嶺區(qū)塊平均儲層壓力梯度為0.7 MPa/hm,因此,生產中動液面下降區(qū)間為0~500 m,此階段應采用小排量來控制生產流量。
(2) 當有效應力為5~12 MPa 時,生產中對應的動液面下降區(qū)值大于500 m,滲透率隨有效應力的增加降低速度變緩,應力敏感性減弱,應采取中等排量方案來控制流量。
(3) 當有效應力大于12 MPa 時,滲透率下降速度極其緩慢,應力敏感性最弱,可將有效應力12 MPa 作為研究區(qū)滲透率動態(tài)變化的理論臨界值,大于此臨界值時采用大排量來控制生產流量。
a.深部煤儲層滲透率隨有效應力的增大呈指數函數規(guī)律降低,平行層理面的滲透率與垂直層理面的呈正相關性。對比分析同一樣品2 個方向上的滲透率,除3 號樣品外,其余2 個樣品平行層理面滲透率均大于垂直層理面滲透率,表明沉積物在沉積過程中形成的層理面是滲透率大小的優(yōu)勢方向。
b.儲層應力敏感性各向異性特征研究表明:在相同壓差的條件下,垂直層理面方向上孔裂隙壓縮系數要高于平行層理面方向;滲透率損害率隨有效應力增加呈指數函數規(guī)律增加,相反,滲透率曲率相應減少;同時垂直層理面方向上的滲透率損害率和滲透率曲率均大于平行層理面方向,表明研究區(qū)垂直層理面方向上應力敏感性更高。
c.基于應力敏感性參數研究,通過煤層氣氣井產能模型,分析得出,不考慮應力敏感性條件下氣井產量偏高,結合應力敏感性對煤層氣產量影響程度,可將研究區(qū)15 號煤儲層應力敏感性變化劃分為3 個階段,相應地在煤層氣排采過程中可采取小-中-大的排量方案來控制生產流量。
d.沁水盆地橫嶺區(qū)塊目前為煤層氣勘探開發(fā)的新區(qū)塊,3 號樣品所在位置儲層滲透率要高于另外兩個,因此,下一步進行深部煤層氣勘探開發(fā)過程中,可以考慮將新井位部署在滲透率相對較高的3 號井區(qū)附近。同時,由于垂直層理面方向應力敏感性高于平行層理面,且平行層理面為滲透率的優(yōu)勢方向,因此,在研究區(qū)煤儲層壓裂改造過程中,應注重孔裂隙在平行層理面方向上的擴展,從而激發(fā)其主導作用。