原位地应力约束下煤储层自封闭作用及其成藏效应

陈世达1,2,3,汤达祯1,2,3,陶 树1,2,3

(1.中国地质大学(北京) 能源学院,北京 100083; 2.煤层气开发利用国家工程研究中心 煤储层实验室,北京 100083; 3.非常规天然气地质评价与开发工程北京市重点实验室,北京 100083)

摘 要:关于煤层气叠置成藏效应的研究通常注重煤系地层层序地层格架的时空配置,对于原位地应力制约下储层的“自封闭效应”关注不足,应力场垂向转换诱导的煤储层渗透性的非单调性变化及其对储层压力、含气性等成藏特征参数的调控作用常被忽视。系统分析了黔西地区煤储层地应力场的垂向分布规律及其构造控制效应,揭示了渗透率随埋深的非指数变化规律及其在沁水、鄂东等含煤盆地的普适性,探讨了储层压力、压力系数垂向差异性分布及其与地应力-渗透率的匹配关系。黔西地区煤储层水平主应力(200~1 300 m)是埋深和构造综合作用的结果,含煤向斜轴部是水平主应力最为集中的区域。根据应力梯度垂向演变规律可以将其划分为应力挤压区(200~500 m,水平构造应力主导)、应力释放区(500~750 m,垂直应力主导)、应力过渡区(750~1 000 m,近向斜轴部)和构造集中区(>1 000 m,向斜轴部低点部位,高应力区)。埋深中段的应力释放区有利于相对高渗储层的形成(平均0.2×10-15 m2),在此深度区间上下(200~500 m,平均0.06×10-15 m2;>750 m,平均0.02×10-15 m2)渗透率普遍较低。渗透率随埋深的这种非单调性变化规律具有普遍性,黔西盘关—土城向斜、比德—三塘向斜埋深中段均存在渗透率相对高值区,沁南(650~800 m)、鄂东(800~950 m)、滇东(600~800 m)、准南(600~800 m)等含煤盆地内这一现象也并不鲜见。原位条件下,低渗储层(<0.1×10-15 m2)的“自封闭”作用使其可以不依赖于盖层等封堵条件,即可构成层间相对独立的流体单元。例如,黔西地区200~500 m,>750 m埋深段的低渗储层自封闭成藏作用显著,储层压力与埋深相关性较差,含气量与压力系数垂向呈“波动式”的无规律变化,常压、欠压、超压储层均有分布,流体压力系统叠置发育。500~750 m的相对高渗储层需外围形成致密的封堵盖层才能阻断层间的流体联系,相对于其他埋深区间更有利于形成统一的流体压力系统,储层压力随埋深增加或层位降低而单调递增,压力系数无明显波动,基本以常压储层为主,含气量随着埋深增大而增大。在同一煤层或气藏内部,储层相对致密部分也可能会封堵高孔渗部分的流体,导致流体压力系统横向上多段叠置。

关键词:煤层气;含气系统;地应力场;渗透率;自封闭作用

黔西—滇东煤层群发育区煤层气资源丰富,但单层资源量较低[1-2]。现行技术条件下,多煤层共采是实现区内煤层气资源高效开发的必要手段之一。然而,煤系地层内“隔水阻气关键层”的存在使得流体垂向联通性受限,煤层群甚至单一煤层的流体压力系统相互封闭和独立,储层压力、压力系数及含气量等储层物性参数呈非单调性变化,即“多层叠置独立含煤层气系统”普遍存在[3-4]。这种现象致使产层组合设计难度增加,排采过程中可能由于合采层位不兼容触发层间矛盾,最终导致合采效果差[5-7]

关于其成藏机制,目前更多的研究在于关注层序地层格架控制下的低渗透性岩层对含气单元间流体的封闭能力。杨兆彪等[8]、沈玉林等[9]先后针对黔西煤层群发育区,进一步阐明了层序地层格架的控藏效应,但关于隔水层的发育位置及其形成机理的结论差异较大。近年来,国内学者以地应力为切入点,发现地应力的垂向转化与含气系统、资源类型存在一定的关联性[10-12]。澳大利亚悉尼盆地也曾报道过水平主应力各向异性在含气性分带中的指示作用[13]。原位埋深条件下,煤系地层特别是煤储层中的割理-裂隙系统必然受到地应力的压缩而变形,从而影响其封隔性,进而使得含煤层气系统发生相应的调整。此前,徐宏杰等[14-15],XU等[16-17],JU等[18]先后对黔西煤储层地应力场的垂向展布及对渗透率的控制作用进行了研究。杨兆彪等[10]、郭晨和卢玲玲[11]进一步揭示了黔西地区原位应力场条件下煤层气的成藏特性和含气系统叠置特征,并认为临界转换深度界面以浅有利于形成统一的流体压力系统。然而,上述研究均未考虑地应力的非线性和渗透率的非指数性变化规律及其所承载的重要成藏(或煤层气富集)信息。

据此,笔者基于煤层气试井和原位地应力测量资料,系统分析了黔西地区含煤向斜构造内地应力随埋深的变化规律,揭示了应力场制约下煤储层渗透率、储层压力、压力系数及含气性的垂向差异性变化规律及其匹配关系,探讨了地层条件下煤储层的自封闭作用及其成藏效应,以期拓展煤层气富集机理和成藏理论认识。

1 地质背景与数据处理

黔西地区属于二叠世上扬子聚煤沉积盆地的一部分,含煤地层为上二叠统长兴组和龙潭组,受高频振荡海平面变化影响,泥炭沼泽稳定期短,造成煤层数量多但单层厚度薄的现状[8-9]。多期构造活动的改造使得原型含煤盆地肢解为众多北东向、北西向的构造盆地,煤系地层在这些规模大小不一的向斜内得以保存,外围则被剥蚀殆尽[19]

本次研究共搜集到黔西地区不同含煤构造32口井共计77层次的试井数据。根据地面垂直钻孔水力压裂测量地应力方法,闭合压力pc为最小水平主应力σh[20],即

σh=pc

(1)

最大水平主应力σH[20]

σH=3pc-pf-p0+T

(2)

式中,pf为破裂压力,MPa;p0为储层压力,MPa;T为煤或岩石的抗拉强度,可由破裂压力和重张压力的差值求取,MPa。

在任意深度,地层垂直主应力σv等于上覆岩石压力,可以通过上覆岩层的密度与埋深进行估算。本文中垂直主应力σv由经验公式[21]计算获得,即

σv≈0.027H

(3)

其中,H为垂直深度,m。水平主应力梯度(G,MPa/hm)可以根据计算得到的水平应力大小(σ)和对应埋深(H)求取,即

(4)

为保证数据的准确性,本次研究剔除了部分最大水平主应力<最小水平主应力的异常数据点(如普定县打磨冲煤矿503孔16号煤层、赫章县野马川507孔1号煤层、大方县对江南2602钻孔M18煤层、龙场煤矿检3孔18号煤层等);200 m以浅仅一口井2组数据,且位于破头山背斜,故也未作针对性分析。因此,共计保留了200~1 300 m埋深范围内的68组数据用以开展统计分析(表1)。

表1 注入/压降试井及原位地应力测量结果
Table 1 Results of injection/fall-off test and in-situ stress measurement

位置煤层深度H/mK/10-15 m2p0/MPaσh/MPaσH/MPaσv/MPa水城都格矿区3-2∗7号807.890.045 96.309.5611.5221.8113-2号869.480.043 46.9613.3319.0523.48普安县幸福井田J202∗18号341.950.000 73.938.2011.069.2319号366.940.000 24.1111.0118.299.91普安县地瓜坡J4105∗9号712.640.008 18.6612.8014.5019.2412号726.630.010 69.3714.2214.9319.6217-1号742.840.002 78.9718.1419.9720.0617-2号752.040.010 97.4013.1816.7620.3119号771.730.007 48.8612.8617.1020.84普定县打磨冲煤矿802孔∗6-1号269.990.060 72.544.254.807.296-2号280.870.051 12.594.255.297.5827号569.900.061 05.258.5210.8015.39赫章县野马川507孔∗4号795.700.031 95.4611.6315.1121.487号826.900.029 15.5413.7620.4522.33龙场煤矿检3孔∗20号215.250.013 32.297.199.815.81大方县对江南2602∗M29号321.960.008 84.569.1412.668.69M51号368.690.063 86.3112.5818.769.95M73号439.510.070 66.5111.4215.9111.87M78号464.270.039 15.3611.5615.8112.54威宁矿区疙瘩营井田J1302∗3号299.860.175 82.796.629.818.107号362.570.073 44.416.737.949.7910号383.680.092 84.687.178.8910.36

续 表

位置煤层深度H/mK/10-15 m2p0/MPaσh/MPaσH/MPaσv/MPa金沙县安洛煤矿7-3孔∗4号338.070.087 44.477.028.449.139号352.120.075 83.708.5612.689.5115号408.160.002 56.4511.1914.2211.02黔西向斜织1井9号216.500.008 71.184.664.495.85珠藏向斜织2井16号380.950.017 92.958.0210.3110.2923号431.380.000 23.0415.5926.2411.65珠藏向斜织4井27号510.000.074 45.5212.0111.6713.77比德向斜织6井32号899.000.000 15.2615.1523.1924.27比德向斜织7井6-1号1 027.000.003 59.1413.1215.9827.73化乐一矿1602孔2号464.040.500 22.978.0912.5312.535号502.260.322 84.418.7512.5913.566号523.350.300 04.689.3113.6514.13化乐一矿3603孔2号520.170.107 45.128.9012.1214.046号577.760.168 25.6911.7516.7615.60织金县中寨煤矿503孔6号220.770.010 01.045.048.175.9616号325.290.031 43.435.416.408.7827号399.940.085 72.246.638.6110.80盘关向斜黔红1井12号1 133.900.001 012.8927.2140.3930.6218号1 178.400.002 910.9627.3638.7931.8224号1 243.600.002 411.0128.3339.7333.58盘关向斜黔红2井12号1 061.900.001 08.6223.1933.2228.6715号1 080.000.009 612.2823.7629.8429.16盘关向斜黔红3井10号1 139.700.002 010.1323.9333.1730.77盘关向斜金佳矿区JV-1井3号558.050.324 05.6610.2314.1715.0710号611.100.259 06.0012.8018.7416.5013号641.580.015 76.1415.7924.2417.3222号706.930.313 07.0411.8516.1919.09土城向斜松河矿区GC-1井1+3号617.750.106 06.9410.6712.7016.689号659.620.093 47.0211.1113.2817.8116号737.100.210 010.1113.1413.5419.9027号917.790.059 912.3521.0125.4224.78盘关向斜月亮田矿区YV-1井6号674.450.192 06.2610.3713.8118.2112号722.090.573 06.3212.1716.3419.5018号774.390.049 26.9612.1815.4120.9124号832.890.057 810.8313.5716.6413.24盘关向斜贵煤1井7号560.000.416 06.6015.6822.5415.12盘关向斜贵煤2井3号359.090.017 33.9510.4015.849.709号408.530.004 45.2713.3317.8011.03青山向斜MY01井17号632.360.480 07.149.017.4117.18青山向斜MY03井17号292.220.007 83.546.288.197.8919号329.400.011 03.668.2910.688.8929号426.950.013 04.579.5613.9211.531206钻孔19号493.290.011 04.0911.5317.5013.3226号628.610.201 06.3711.5315.2216.97黔西向斜小林华矿区XC-1井4号890.410.026 88.0412.4112.6724.0415号964.820.027 68.1117.5223.0426.05

注:标注*的数据引自文献[10]。

2 原位地应力场垂向分布规律

2.1 含煤向斜构造水平主应力变化

为查明含煤向斜构造对水平主应力的影响,笔者对盘关向斜8井次24层次的试井数据进行了剖析(图1)。其中,黔红1,2,3井位于向斜轴部的亮山区块;JV-1,GC-1,YV-1和贵煤1,2井位于向斜翼部。实测结果表明,向斜两翼最大、最小水平主应力梯度随埋深增大逐渐降低,在埋深中段应力梯度大小达到最小值,应力场逐渐过渡为正断层应力场类型(σv>σH>σh)。位于向斜轴部的黔红1,2,3井的水平主应力大小及梯度较之埋深中段明显增高,属于高应力区(σh>18 MPa)的范畴,三向应力相对大小关系表现为σH>σv>σh,即走滑断层应力场类型。从统计结果来看,应力梯度随埋深变化整体上表现为“弧形”,浅部和靠近向斜轴部煤层构造应力占绝对优势,是水平主应力最为集中的区域,埋深中段构造应力集中现象相对不明显,以垂直主应力为主。

图1 盘关向斜最大、最小水平主应力及应力梯度垂向变化
Fig.1 Vertical variation of horizontal stresses and horizontal stresses gradient in the Panguan syncline

地应力是存在于地壳中的内应力,包括自重应力与水平应力。自重应力场比较简单,可用上覆岩层的容重与埋藏深度估算。水平应力场则相对较为复杂,一方面是在重力作用下煤岩层水平方向移动受限而产生;另一方面与构造运动及岩层地质构造有关,尤以水平方向构造应力影响最大[22]。因此,水平主应力大小是自重应力产生的泊松效应和构造运动产生的构造应力叠加的结果。在向斜轴部构造应力的挤压作用下,自重应力形成的水平应力基础上附加了较大构造应力,煤岩层内水平应变增大,导致水平应力迅速增加。HAN 等[23]也发现在距离平顶山向斜轴部的一定范围内,水平主应力的梯度将显著增高。图2给出了盘关向斜同一煤层(18号,24号)或邻近煤层(1+3号,3号,6号,7号,9号,10号,12号,13号,15号,16号)水平主应力梯度随埋深的变化趋势,其客观上反映了煤层水平主应力自浅部至深部、自向斜翼部至轴部的变化过程。可见,在靠近向斜轴部过程中煤层应力梯度先降后增,约750 m埋深处向斜轴部应力集中作用开始显现,水平主应力梯度随埋深增大而递增且逐渐大于垂直主应力梯度。

图2 相同或邻近煤层应力梯度随埋深变化
Fig.2 Changes in horizontal stress gradient of the same or the adjacent seams with increasing depth

2.2 原位地应力状态垂向转换规律

从实测结果来看,黔西地区平均最小水平主应力梯度(2.1 MPa/hm)较之我国其他含煤盆地略高,但由于不同含煤盆地应力测量数据在不同埋深区间出现的频率不同,应力大小或梯度的平均值难以反映煤岩所承受的现今地应力强度(表2)。从应力值大小及梯度与埋深的关系的角度,在给定埋深范围内,黔西地区最小水平主应力变化区间范围与我国其他含煤盆地基本相当,并非传统认识中典型的“高应力区”(图3)。依据地应力量级判定标准(σh>30 MPa为超高应力区;18~30 MPa为高应力区;10~18 MPa为中应力区;0~10 MPa为低应力区)[22],400 m以浅表现为为低应力区,400~1 000 m为中应力区,1 000~1 300 m才进入高应力区的范畴。从拟合趋势来看,1 400 m以深的煤层将承受超高地应力的影响。

表2 我国不同含煤盆地煤储层地应力实测结果对比
Table 2 Comparison of in-situ stress measurements results of coals in different coal-bearing basins

区域测点数量/个埋深/m最小主应力/MPa最小主应力梯度/(MPa·hm-1)最大主应力/MPa最大主应力梯度/(MPa·hm-1)黔西68215~1 2446044.25~28.3312.081.18~3.612.104.49~40.3916.141.17~6.082.80鄂东77488~1 7288365.8~28.415.201.00~3.261.883.94~33.3322.570.79~5.042.80沁南67353~1 2737156.24~29.0914.721.19~3.082.0410.48~44.5723.371.72~5.853.30滇东58130~1 1836033.16~19.0110.891.47~2.481.874.22~25.5514.801.77~3.622.55准南64278~1 5608284.67~29.1813.631.00~2.571.665.95~43.7718.860.63~4.122.32

注:数据格式

图3 不同含煤盆地煤储层最小水平主应力大小及梯度垂向变化
Fig.3 Vertical variation in minimum horizontal stress magnitude and gradient

在统计的68个测点中,逆断层应力场类型(σH>σh>σv)有9个,埋深介于215~560 m;正断层应力场类型(σv>σH>σh)有37个,埋深介于210~1 027 m;走滑断层应力场类型(σH>σv>σh)有22个,埋深介于221~1 244 m。水平主地应力梯度随埋深的变化规律与盘关向斜实测结果基本一致。500 m以浅和1 000 m以深共计进行了34层次的试井,逆断层和走滑断层应力场类型占76%,以水平主应力为主导,属于典型的构造应力场类型;500~1 000 m埋深范围内进行了34层次的试井,正断层应力场类型占85%,以垂直主应力为主导,属于大地静力场类型。

考虑应力梯度变化、应力的相对大小关系及向斜构造的影响,可以将煤储层地应力状态划分为4个深度区间(图4,表3):① 200~500 m埋深内煤储层呈应力挤压状态,表现为水平应力随埋深增大而增大,且增加幅度(或应力梯度)大于垂直应力,400 m附近正断层应力场类型消失,水平构造应力占据主导地位;② 500~750 m煤储层位于应力释放区,水平主应力梯度较之浅部有所收敛且明显降低,应力值变化趋于平缓,走滑、逆断层应力场类型基本不可见,应力集中现象不明显;③ 750~1 000 m地应力场仍以垂直主应力为主,但向斜轴部构造应力集中作用开始显现,水平主应力梯度变化趋势发生反转,最大水平主应力逐步接近垂直应力,属于应力过渡区;④ 1 000 m以深进入向斜轴部的构造应力集中区,水平主应力增加幅度远大于垂直主应力,最大水平主应力为最大主应力,最小水平主应力增高至所测量数据的最大值,属于高应力区的范畴。

表3 黔西地区煤储层地应力场垂向分带及其地质意义
Table 3 Vertical variation of in situ stress field in western Guizhou and its geology significance

埋深分段划分依据应力场类型地质意义应力挤压区或增高区(200~500 m)水平应力大小及梯度逐步增大(水平构造应力附加值增大)走滑和逆断层应力场类型渗透率<0.1×10-15 m2、煤层气自封闭成藏应力释放区或降低区(500~750 m)应力卸荷(或松弛)使水平主应力值明显降低(侧压系数<1)正断层应力场类型渗透率>0.1×10-15 m2、压力系统相对统一应力过渡区(750~1 000 m)轴部应力集中作用开始显现,水平应力梯度增大,由正断层应力场向走滑断层应力场类型过渡正断层应力场类型渗透率<0.1×10-15 m2、煤层气自封闭成藏应力集中区(1 000~1 300 m)向斜轴部低点部位,构造应力集中强度增高,形成高应力区走滑断层应力场类型渗透率<0.01×10-15 m2、煤层气自封闭成藏

图4 地应力状态划分
Fig.4 Division of in-situ stress regimes

3 原位地应力控渗、控藏效应

黔西地区煤储层渗透率在0.000 1×10-15~0.573×10-15 m2,平均仅0.09×10-15 m2。渗透率的垂向分布具有明显的阶段性,200~500,500~750,750~1 000 m以及1 000~1 300 m埋深范围分别对应低(平均0.05×10-15 m2)、高(平均0.2×10-15 m2)、低(平均0.03×10-15 m2)、极低(平均0.003×10-15 m2)渗透率(图5)。渗透率的这种垂向阶段分布特点并非个例,在盘关—土城向斜以及比德—三塘向斜均较为明显(图6)。在500~750 m埋深内,2者均存在一个渗透率相对高值区(> 0.1×10-15 m2),如盘关—土城向斜金佳1井(3号,10号,22号)、松参1井(1+3号,9号,16号)、月亮田1井(6号,12号)以及贵煤1井(7号),比德-三塘向斜的织4井(27号)、化乐一矿1602孔(2号,5号,6号)以及化乐一矿3603孔(2号,6号);该埋深区间上下煤层渗透率则普遍较低。

图5 黔西地区煤层试井渗透率、水平主应力差及侧压系数垂向变化
Fig.5 Vertical variation of permeability,horizontal stress difference and lateral pressure coefficient

图6 盘关—土城向斜、比德-三塘向斜渗透率垂向变化
Fig.6 Vertical variation of permeability in Panguan-Tucheng and Bide-Santang synclines

从应力场类型来看,正断层应力场类型相对有利于形成较高渗透的储层,其次为走滑应力场类型,逆断层应力场类型的煤层渗透率均较低。垂向上,储层裂隙的开合程度(即煤储层的渗流能力)不单纯与某一方向地应力绝对值相关,也与应力梯度变化趋势所反映出来的储层受力“状态”相关(图5):当应力梯度趋于降低时,应力状态相对拉张,渗透率则有所恢复;当应力梯度趋于增高时,则处于相对挤压的状态,渗透率普遍较低。水平主应力差也是控制裂缝开合程度和渗透率大小的一个关键参数,在较高的水平主应力差条件下,绝大多数天然裂缝通常处于闭合状态,导致其导流能力明显下降。就黔西地区而言,200~500 m埋深内应力梯度和水平主应力差明显增大,导致煤岩裂缝闭合,煤储层渗透率相应降低;500~750 m埋深区间内应力集中程度较低,水平主应力梯度和水平主应力差均有所下降,相对有利于煤岩裂缝保持张开的状态;750~1 000 m埋深段内,向斜轴部构造挤压作用开始显现,水平主应力大小、梯度及水平主应力差值再次升高,导致煤岩渗透率显著降低;1 000 m以深进入高地应力区的范畴,应力高度集中,渗透率极低(<0.01×10-15 m2)。原位地应力制约下煤储层渗透率的垂向差异性分带与储层压力、含气性之间存在较强的关联性。黔西地区煤层的储层压力为0.72~12.89 MPa,储层压力梯度介于0.28~1.70 MPa/hm。图7表明,储层压力与埋深的线性关系并不明显,特别是在200~500 m及>750 m埋深范围内存在极强的离散型,这种离散性反映在储层压力梯度上则更为显著。200~500 m和>750 m埋深范围内,储层压力梯度分布离散,常压、欠压以及超压储层均有分布且不存在明显的变化规律;500~750 m埋深范围内,储层压力与埋深相关性较好,压力系统统一程度高,平均储层压力梯度约为1 MPa/hm,与静水压力梯度大小相当。

图7 煤储层地应力梯度、渗透率、储层压力、压力系数及含气量垂向耦合关系
Fig.7 Vertical variation of stress gradient,permeability,pressure,pressure gradient and gas content

同一钻孔内,不同煤层的储层压力随埋深增加(或层位降低)的变化也基本符合上述规律(图8)。例如,盘关向斜金佳1井3号~22号煤层位于550~750 m埋深范围内,储层压力随埋深增加单调递增,压力梯度无明显波动,说明其属于统一的流体压力系统;松参1井1+3号~27号埋深在550~950 m,储层压力梯度在约750 m附近发生突变,由常压储层变为超压储层,类似的变化也发生在月亮田1井的6号~24号煤;黔红1井12号~24号煤层埋深大于1 100 m,渗透率极低,储层压力大小及梯度不增反降,说明层间压力系统叠置发育;比德向斜化乐1矿1602钻孔在460~540 m埋深段内,储层压力单调递增且至5号煤(约500 m)压力梯度开始趋于稳定;中寨煤矿503钻孔6号~27号煤的埋深介于200~450 m,随着埋深增大或层位的降低,储层压力大小及梯度呈波动变化,超压(16号煤)、欠压(6号,27号煤)储层均有发育地应力、渗透率、储层压力的垂向差异性变化与煤岩含气量分布存在较高的一致性。土城向斜松河矿区的松参1井进行了600~960 m埋深内共计30层次的现场含气量测试,其中部分煤层仅完成了解吸气量和损失气量的分析,未对残余气量进行进一步分析。尽管这部分样品的含气量较总含气量略低,但2者整体的变化趋势完全一致。由图9可知,600~750 m埋深内的煤层属于相对统一的含气系统(1+3号煤→12号煤),含气量随埋深的变化基本符合单调递增的一般性规律,含气量由5.93 m3/t增高至18.83 m3/t;当埋深>750 m,煤层含气量与埋深失去相关性,呈无明显规律的波动变化,说明该埋深段内含气系统叠置发育,相关成藏参数难以有效预测。含气量的这种垂向分带式变化规律在盘关向斜月亮田1井(YV-1井)及金佳1井(JV-1井)中也较为明显。

图8 同一井筒内不同煤层储层压力及压力梯度变化特征
Fig.8 Variation of reservoir pressure and pressure gradient of different coal seams in the same borehole

图9 GC-1井、YV-1井、JV-1井含气量(解吸气+残余气)垂向变化规律
Fig.9 Variation of coal seam gas content(desorbed gas remined gas) in Well GC-1,Well YV-1 and Well JV-1

4 地应力控藏机制探讨及研究展望

原位地应力场、储层渗透率、流体压力系统及含气性分带之间存在着必然的内在联系,其对于煤层气资源可采性的评价和甜点区段的优选具有重要的意义。黔西地区500~750 m埋深内应力集中现象不明显,垂直应力占主导地位,渗透率相对较高(>0.1×10-15 m2),储层压力与埋深存在良好的线性正相关性,压力梯度(或压力系数)波动幅度不大,含气量随埋深的增大而增高,流体压力系统统一程度较高。相反,在200~500 m和>750 m埋深内,水平主应力大小、梯度及水平主应力差均随埋深增加而增大,导致形成低渗-极低渗煤储层,储层压力、压力系数及含气量与埋深的相关性较差,欠压、常压、超压储层均有分布,流体压力系统普遍叠置发育。从这一点来看,500~750 m埋深区间多煤层合采的甜点层段,这与现阶段的高产煤层气井的产层埋深分布具有较高的一致性,如杨梅树向斜杨煤参1井的13-2号(657.07~659.00 m)、7号(635.87~637.7 m)、5-2号(596.35~599.01 m)[24],三塘向斜织3井14号(697.7~700.3 m)、16号(735.4~737.5 m)煤层,珠藏向斜小试验井组10口合采煤层气井的20号,23号,27号,30号煤层(埋深450~650 m)[25]

具体而言,原位地应力场控藏效应的实质是以控制储层渗透率为桥梁,进而决定储层自封闭能力的一种地质作用。自封闭(Self-sealing)最早由FACCA和TONANI[26-27]提出,常规油气地质领域多概指地层温压变化及地下流体共同作用下岩体发生重结晶和再胶结,最终封闭或封堵盖层之下储集岩体内流体的地质作用,它使地下流体与盖层之外失去联系后构成相对独立的流体单元或成藏单元。这类“自封闭”作用是目标岩体外围产生了相对致密部分,并对其高孔渗部分形成类似盖层的封堵能力。秦勇等[28]提出了高煤级煤层气成藏“弹性自封闭效应”的观点,认为这种效应是我国高煤级煤地区煤层气普遍富集的重要地质原因。贾承造等[27]将这一概念引入了非常规油气领域,并将其定义为非常规油气在沉积盆地内由于自身特殊的物理化学特性或在特殊储集层介质条件和特殊温压环境共同作用下,依赖油气自身内部或油气与储集层介质界面之间的分子间作用力,不依赖储集体之外的圈闭等上倾封堵条件,与外界隔离并独立成藏富集保存的地质作用。

不同于常规油气藏,非常规油气的自封闭是在一定的边界条件或地质门限条件下,发生在非常规油气藏内部而不是在其外部或边缘的一种成藏作用。对于煤层气而言,自封闭成藏作用需要致密储集条件。例如,黔西地区200~500 m埋深区间内煤层试井渗透率小于0.1×10-15 m2,750 m以深煤层渗透率<0.01×10-15 m2,这种低孔隙度、低渗透率条件使得煤层气藏不依赖于外部“封隔层”即可阻断储集层内流体的自由流出和外部流体的自由进入,形成独立流体压力系统垂向多层叠置的现象。这与杨兆彪等[29]的研究结论基本一致,即原位渗透率0.1×10-15 m2可以作为划分煤层是否具有统一含气系统的标志。对于500~750 m埋深范围内的相对高渗储层(>0.1×10-15 m2),需外围形成相对致密的封堵盖层(即考虑层序地层格架特点)阻断层间的流体联系才能构成独立的流体单元或成藏单元,因此相对于其他埋深区间更有利于形成统一的流体压力系统。

除黔西地区煤层群发育区外,地应力的非线性和渗透率的非指数型垂向变化在其他含煤盆地也并不鲜见,如滇东(600~800 m)、沁水盆地南部(650~800 m)、鄂尔多斯盆地东缘(800~950 m)及准噶尔盆地南缘(600~800 m)等地区埋深中段均存在水平主应力梯度低值区和渗透率相对高值区,且水平应力梯度随埋深的阶段式增减与渗透率变化均存在较好的对应关系(图10)。地应力在垂向上的这种复杂变化与盆地类型有关,我国含煤盆地多属于挤压型盆地,挤压构造力的传播伴随着压缩吸收能量和断层释放能量的过程,在区域应力背景下通常叠加了局部构造应力作用,因此地应力类型和构造特征具有极强的分带性[30-31]。从实测结果来看,我国几个主要含煤盆地应力梯度变化多表现为数字“3”形态:近地表靠近背斜轴部或近轴翼部等区域属于为拉张应力区,渗透率普遍较高;随着埋深增大,构造应力附加值逐步增大,过渡为走滑或逆断层应力场类型,煤体发生塑性变形,渗透率快速降低;随着埋深进一步增大,地应力状态发生转换,水平应力梯度和水平主应力差减小,渗透率恢复至相对较高的水平;局部(如向斜轴部)挤压应力的叠加会使深部煤层应力梯度将再次升高,同时伴随渗透率的骤降。

图10 滇东[32]、沁南、鄂东、准南地区应力梯度、渗透率随埋深变化
Fig.10 Vertical variation of permeability and horizontal stress gradient in the eastern Yunnan[32],the southern Qinshui basin,the eastern margin of Ordos basin,and the southern Junggar basin

地应力-渗透率的这种垂向非单调函数变化模式明显不同于埋深增大→地应力增高→渗透率负指数降低的传统认识,而不同埋深段储层的致密情况对束缚流体动力场、调控油气聚集成藏的作用不可忽视。对于我国西南地区的煤层群发育区,低渗储层的自封闭成藏效应势必会导致流体压力系统层间叠置发育的更加复杂化。同一储层在不同埋深或不同构造位置处所承受的现今地应力也显著不同,储集层内部的致密情况分布必然也有所差异,相对致密部分也可能会封堵高孔渗部分的流体,导致同一气藏内流体压力系统在横向上多段叠置,这对以单煤层单采为主的沁南、鄂东等含煤盆地内煤层气的富集机理、资源类型也具有重要意义,在煤层气勘探开发过程中应予以足够重视。

5 结 论

(1)黔西地区地应力量级与我国其他含煤沉积盆地基本相当,并非典型的“高应力区”,1 000 m以深才进入高应力区的范畴,向斜轴部是水平主应力最为集中的区域。

(2)根据应力梯度变化可将黔西地区地应力状态划分应力挤压区、应力释放区、应力过渡区和构造集中区。应力释放区有利于相对高渗储层的形成,在此深度区间上下渗透率普遍较低(<0.1×10-15 m2)。

(3)地应力控藏效应的实质是一定渗透率门限条件下的煤储层自封闭成藏作用。低渗储层(<0.1×10-15 m2)可阻断内部流体的流出和外部流体的进入,促使流体压力系统叠置发育。相对高渗储层(>0.1×10-15 m2)则需外围形成致密的封堵盖层才能构成独立的流体单元或成藏单元,相对有利于形成统一的流体压力系统。

(4)地应力的非线性和渗透率的非指数型垂向变化在沁水、鄂东、准南、滇东等含煤盆地也较为普遍,其所承载的重要成藏信息在煤层气勘探过程中应予以足够重视。

参考文献(References):

[1] 吴财芳,刘小磊,张莎莎.滇东黔西多煤层地区煤层气“层次递阶”地质选区指标体系构建[J].煤炭学报,2018,43(6):1647-1653.

WU Caifang,LIU Xiaolei,ZHANG Shasha.Construction of index system of “Hierarchical progressive” geological selection of coalbed methane in multiple seam area of eastern Yunnan and western Guizhou[J].Journal of the China Coal Society,2018,43(6):1647-1653.

[2] 桂宝林.煤层气勘探目标评价方法——以滇东黔西地区为例[J].天然气工业,2004,24(5):33-35.

GUI Baolin.Exploration target evaluation of coalbed methane in east Yunnan—west Guizhou region[J].Natural Gas Industry,2004,24(5):33-35.

[3] 秦勇,熊孟辉,易同生,等.论多层叠置独立含煤层气系统——以贵州织金—纳雍煤田水公河向斜为例[J].地质论评,2008,54(1):65-70.

QIN Yong,XIONG Menghui,YI Tongsheng,et al.On unattached multiple superposed coalbed-methane system:In a case of the Shuigonghe Syncline,Zhijin-Nayong Coalfield,Guizhou[J].Geological Review,2008,54(1):65-70.

[4] 郭晨,夏玉成,卢玲玲,等.黔西比德—三塘盆地多层叠置独立含煤层气系统发育规律与控制机理[J].天然气地球科学,2017,28(4):622-632.

GUO Chen,XIA Yucheng,LU Lingling,et al.Development features and mechanism of multi-layer superimposed independent CBM systems in Bide-Santang Basin,western Guizhou,South China[J].Natural Gas Geoscience,2017,28(4):622-632.

[5] 秦勇,吴建光,张争光,等.基于排采初期生产特征的煤层气合采地质条件分析[J].煤炭学报,2020,45(1):241-257.

QIN Yong,WU Jianguang,ZHANG Zhengguang,et al.Analysis of geological conditions for coalbed methane co-production based on production characteristics in early stage of drainage[J].Journal of the China Coal Society,2020,45(1):241-257.

[6] 许江,张超林,彭守建,等.多层叠置煤层气系统合采方式及其优化[J].煤炭学报,2018,43(6):1677-1686.

XU Jiang,ZHANG Chaolin,PENG Shoujian,et al.Multiple layers superposed CBM system commingled drainage schedule and its optimization[J].Journal of the China Coal Society,2018,43(6):1677-1686.

[7] 杨兆彪,李洋阳,秦勇,等.煤层气多层合采开发单元划分及有利区评价[J].石油勘探与开发,2019,46(3):559-568.

YANG Zhaobiao,LI Yangyang,QIN Yong,et al.Development unit division and favorable area evaluation for joint mining coalbed methane[J].Petroleum Exploration & Development,2019,46(3):559-568.

[8] 杨兆彪,秦勇,高弟,等.煤层群条件下的煤层气成藏特征[J].煤田地质与勘探,2011,39(5):22-26.

YANG Zhaobiao,QIN Yong,GAO Di,et al.Coalbed methane(CBM) reservoir-forming character under conditions of coal seam groups[J].Coal Geology & Exploration,2011,39(5):22-26.

[9] 沈玉林,秦勇,郭英海,等.黔西上二叠统含煤层气系统特征及其沉积控制[J].高校地质学报,2012,18(3):427-432.

SHEN Yulin,QIN Yong,GUO Yinghai,et al.The upper Permian coalbed methane bearing system and its sedimentary control in Western Guizhou,China[J].Geological Journal of China Universities,2012,18(3):427-432.

[10] 杨兆彪,秦勇.地应力条件下的多层叠置独立含气系统的调整研究[J].中国矿业大学学报,2015,44(1):70-75.

YANG Zhaobiao,QIN Yong.A study of the unattached multiple superposed coalbed-methane system under stress conditions[J].Journal of China University of Mining & Technology,2015,44(1):70-75.

[11] 郭晨,卢玲玲.黔西煤层气成藏特性空间分异及其对开发的启示[J].煤炭学报,2016,41(8):2006-2016.

GUO Chen,LU Lingling.Spatial distribution and variation of coalbed methane reservoir characteristics and its significance for CBM development in Western Guizhou[J].Journal of the China Coal Society,2016,41(8):2006-2016.

[12] 吴财芳,王肖,刘小磊,等.滇东老厂矿区多煤层条件下地应力特征及其影响研究[J].煤炭科学技术,2019,47(1):118-124.

WU Caifang,WANG Xiao,LIU Xiaolei,et al.Study on geostress features and influences under multi-seam condition in Laochang Mining Area of East Yunnan[J].Coal Science and Technology,2019,47(1):118-124.

[13] BURRA Agi,ESTERLE Joan,GOLDING Suzanne.Horizontal stress anisotropy and effective stress as regulator of coal seam gas zonation in the Sydney Basin,Australia[J].International Journal Coal Geology,2014,132:103-116.

[14] 徐宏杰,桑树勋,易同生,等.黔西地应力场特征及构造成因[J].中南大学学报:自然科学版,2014,45(6):1960-1966.

XU Hongjie,SANG Shuxun,YI Tongsheng,et al.Characteristics of in-situ stress field and its tectonic origin in Western Guizhou[J].Journal of Central South University(Science and Technology),2014,45(6):1960-1966.

[15] 徐宏杰,桑树勋,易同生,等.黔西地区煤层埋深与地应力对其渗透性控制机制[J].地球科学,2014,23(11):1507-1516.

XU Hongjie,SANG Shuxun,YI Tongsheng,et al.Control mechanism of buried depth and in-situ stress for coal reservoir permeability in western Guizhou[J].Earth Science,2014,23(11):1507-1516.

[16] XU Hongjie,SANG Shuxun,YANG Jingfen,et al.In-situ stress measurements by hydraulic fracturing and its implication on coalbed methane development in Western Guizhou,SW China[J].Journal of Unconventional Oil and Gas Resources,2016,15:1-10.

[17] XU Hongjie,SANG Shuxun,YANG Jingfen,et al.Selection of suitable engineering modes for CBM development in zones with multiple coalbeds:A case study in western Guizhou Province,Southwest China[J].Journal of Natural Gas Science and Engineering,2016,36:1264-1275.

[18] JU Wei,YANG Zhaobiao,QIN Yong,et al.Characteristics of in-situ stress state and prediction of the permeability in the Upper Permian coalbed methane reservoir,western Guizhou region,SW China[J].Journal of Petroleum Science and Engineering,2018,165:199-211.

[19] 窦新钊,姜波,秦勇,等.黔西地区构造演化及其对晚二叠世煤层的控制[J].煤炭科学技术,2012,40(3):109-114.

DOU Xinzhao,JIANG Bo,QIN Yong,et al.Structure evolution in west of Guizhou area and control to seam in late Permian[J].Coal Science and Technology,2012,40(3):109-114.

[20] HAIMSON B,FAIRHURST C.In-situ stress determination at great depth by means of hydraulic fracturing[A].The 11th US symposium on rock mechanics(USRMS)[C].American Rock Mechanics Association,1969.

[21] HOEK E,BROWN E T.Underground excavations in rock[D].London:The Institution of Mining Metallurgy,1980.

[22] KANG H,ZHANG X,SI L,et al.In-situ stress measurements and stress distribution characteristics in underground coal mines in China[J].Engineering Geology,2010,116(3-4):333-345.

[23] HAN Jun,ZHANG Hongwei,LIANG Bin,et al.Influence of large syncline on in situ stress field:A case study of the Kaiping coalfield,China[J].Rock Mechanics and Rock Engineering,2016,49(11):4423-4440.

[24] 单衍胜,毕彩芹,迟焕鹏,等.六盘水地区杨梅树向斜煤层气地质特征与有利开发层段优选[J].天然气地球科学,2018,29(1):122-129.

SHAN Yansheng,BI Caiqin,CHI Huanpeng,et al.Geological characteristics of coalbed methane and optimization for favorable productive intervals of Yangmeishu syncline in Liupanshui area[J].Natural Gas Geoscience,2018,29(1):122-129.

[25] 李伶俐.织金区块不同煤组煤层气产能评价及影响因素研究[J].煤炭技术,2020,39(7):59-61.

LI Lingli.Productivity evaluation and influencing factors of different coal groups in Zhijin block[J].Coal Technology,2020,39(7):59-61.

[26] FACCA G,TONANI F.The self-sealing geothermal field[J].Bulletin Volcanologique,1967,30(1):271-273.

[27] 贾承造,庞雄奇,宋岩.论非常规油气成藏机理:油气自封闭作用与分子间作用力[J].石油勘探与开发,2021,48(3):437-452.

JIA Chengzao,PANG Xiongqi,SONG Yan.The mechanism of unconventional hydrocarbon formation:Hydrocarbon self-containment and intermolecular forces[J].Petroleum Exploration and Development,2021,48(3):437-452.

[28] 秦勇,傅雪海,吴财芳,等.高煤级煤储层弹性自调节作用及其成藏效应[J].科学通报,2005,50(S1):82-86.

QIN Yong,FU Xuehai,WU Caifang,et al.Self-adjust elastic effect of high-rank coal reservoirs and its accumulation effect[J].Science Bulletin,2005,50(S1):82-86.

[29] 杨兆彪,秦勇,陈世悦,等.多煤层储层能量垂向分布特征及控制机理[J].地质学报,2013,87(1):139-144.

YANG Zhaobiao,QIN Yong,CHEN Shiyue,et al.Controlling mechanism and vertical distribution characteristics of reservoir energy of multi-coalbeds[J].Acta Geologica Sinica,2013,87(1):139-144.

[30] KANG Yongshang,JIANG Shanyu,DENG Ze,et al.Non-monotonic decrease of coal permeability with depth in highly compressional tectonic settings:A case study in the Bide-Santang Basin,SW China[J].Marine and Petroleum Geology,2020,116:104321.

[31] 王平.地质力学方法研究——不同构造力作用下地应力的类型和分布[J].石油学报,1992,13(1):1-12.

WANG Ping.A geomechanical technique—types and distribution of geostress under various tectonic forces[J].Acta Petrolei Sinica,1992,13(1):1-12.

[32] JU Wei,JIANG Bo,MIAO Qi,et al.Variation of in situ stress regime in coal reservoirs,eastern Yunnan region,South China:Implications for coalbed methane production[J].AAPG Bulletin,2018,102(11):2283-2303.

Self-sealing of coals and CBM accumulation effect constrained by in situ stress

CHEN Shida1,2,3,TANG Dazhen1,2,3,TAO Shu1,2,3

(1.School of Energy Resources,China University of Geosciences (Beijing),Beijing 100083,China; 2.Coal Reservoir Laboratory,National Engineering Research Center of CBM Development & Utilization,Beijing 100083,China; 3.Beijing Key Laboratory of Unconventional Natural Gas Geological Evaluation and Development Engineering,Beijing 100083,China)

Abstract:Research on the superimposed coalbed methane system usually pays attention to the sequence stratigraphic framework of the coal-measure stratigraphy.There is insufficient attention to the “elastic self-sealing effect” of the reservoir under the restriction of in-situ in situ stress.The non-monotonic changes in coal reservoir permeability and its control on reservoir pressure,gas content are often overlooked.Here,the vertical variation in in-situ stress and its structural control factors in western Guizhou were systematically analyzed,and the non-exponential variation in permeability with increasing depth in several coal-bearing basins were revealed.Finally,the matching relationship between reservoir pressure,pressure coefficient and permeability in the vertical were further discussed.In western Guizhou,the horizontal stress gradient changes with depth is non-uniform and is overprinted by the effect of syncline,and the stress gradient is remarkably high nearing the axis.These observations within 200-500,500-750,750-1 000 and >1 000 m depths are similar with the variations in horizontal stress difference and lateral pressure coefficient,corresponding to stress extrusion (horizontal stresses dominated),stress release (vertical stress dominated),stress transition (stress concentration in syncline axis begins to appear),and stress concentration zones (the low point of the syncline axis,high stress area),respectively.The stress release zone is favorable for a relatively high permeability reservoir (mean 0.2×10-15 m2).Below (750-1 000 m,mean 0.04×10-15 m2;>1 000 m,mean 0.003×10-15 m2) and above(200-500 m,mean 0.06×10-15 m2)this zone,the reservoirs have a very low permeability.The non-monotonic change law in permeability with depth is universal,the relatively high permeability zones at the middle depths can be found in the Panguan-Tucheng and Bide-Santang synclines,as well as in the southern Qinshui Basin (650-800 m) and the eastern margin of the Ordos Basin (800-950 m).Under in-situ conditions,the “self-sealing” effect of low-permeability reservoirs (<0.1×10-15 m2) makes it possible to form relatively independent fluid units or accumulation units without relying on caprocks and other sealing conditions.For seams at 200-750 m and >750 m in depth,penetration is restricted by reduced permeability,resulting in discontinuous gas-bearing systems with irregular gas content distributions and unpredictable reservoir pressure gradient.The higher permeability reservoirs at 500-750 m depths need tight sealing cap rock to block the fluid connection between layers,and thereby is favorable to form a unified gas-bearing system,where the gas content generally increase with depth to a“peak gas”horizon under hydrostatic pressure.The relatively tight part of the reservoir may also block the high porosity and permeability part,causing multiple section of fluid pressure system in the same gas reservoir.

Key words:CBM;gas-bearing system;insitu stress;permeability;self-sealing effect

中图分类号:P618.11

文献标志码:A

文章编号:0253-9993(2021)08-2466-13

移动阅读

收稿日期:2021-04-20

修回日期:2021-06-26

责任编辑:韩晋平

DOI:10.13225/j.cnki.jccs.CB21.0686

基金项目:国家科技重大专项资助项目(2016ZX05044-001);国家自然科学基金资助项目(41772132);中国博士后科学基金资助项目(2021M692998)

作者简介:陈世达(1991—),男,山东日照人,讲师。E-mail:Shida.Chen@cugb.edu.cn

通讯作者:汤达祯(1957—),男,江苏盐城人,教授,博士生导师。E-mail:tang@cugb.edu.cn

引用格式:陈世达,汤达祯,陶树.原位地应力约束下煤储层自封闭作用及其成藏效应[J].煤炭学报,2021,46(8):2466-2478.

CHEN Shida,TANG Dazhen,TAO Shu.Self-sealing of coals and CBM accumulation effect constrained by in situ stress[J].Journal of China Coal Society,2021,46(8):2466-2478.