废弃石油井周围煤层中油层气扩散-渗流规律

索永录 刘颖凯 肖江

摘 要：为揭示煤矿废弃石油井周围煤层内油層气流动的压力与影响半径的变化规律，基于煤体孔隙裂隙二重介质假设，考虑孔隙率与渗透率动态变化，结合煤体变形特征，建立了油层气流动的扩散-渗流耦合模型，针对油层气组分的特殊性，运用COMSOL Multiphysics软件进行模拟计算。结果表明：流动过程中，煤层各处油层气压力的增加速度都先加快后减慢直至趋于稳定，距离废弃石油井越近，油层气压力短时间内变化越快，压力也越大;随着时间的推移，油层气影响半径逐渐增加，最后趋于稳定，井内压力恒定时，影响半径增加速度一直衰减，井内压力线性增加时，影响半径增加速度先增后减;第30年时，当井内压力线性增加至15 MPa时，影响半径为310 m.模拟结果与现场实测结果基本吻合，说明该耦合模型合理，该研究为受废弃石油井威胁的煤层的安全高效开采提供理论指导。

关键词：废弃石油井;油层气;扩散-渗流;数值模拟;影响半径

中图分类号：TD 712           文献标志码：A

文章编号：1672-9315（2019）05-0753-08

Abstract：In order to reveal the changing laws of the pressure and influence radius of the oilbed gas flowing in the coal seam around the abandoned oil wells in coal mines，based on the hypothesis of dualporous media，the diffusionseepage coupling model of oilbed gas flow was established with the dynamic changes of porosity and permeability as well as the deformation characteristics of coal in mind.The coupling model with the specificity of oilbed gas composition was simulated by COMSOL Multiphysics software.The results shows that：during the flow，the increase rate of the oilbed gas pressure in the coal seam accelerates first，then slows down until it stabilizes.The closer to the abandoned oil well，the faster the pressure changes in a short time and the greater the pressure is.As time went on，the influence radius gradually increases，and finally stabilizes.The increase speed of influence radius has attenuated when the pressure in the well is constant while it increases first and then decreases when the pressure in the well increases linearly.At 30 years，when the pressure in the well increases linearly to 15 MPa，the influence radius will be 310 m.The simulation results are roughly consistent with the reality.So this coupling model is proved to be right.This study provides a theoretical guidance for safe and efficient mining in the coal seam threatened by abandoned oil wells.

Key words：abandoned oil well;oilbed gas;diffusionseepage;numerical simulation;influence radius

0 引 言

我国有不少煤矿处于煤炭资源与石油资源重叠区[1-2]。宁东地区马家滩矿区的双马煤矿井田内有170口废弃石油井。目前井下位于I0104105综采工作面内的马探31号和马探30号石油井，在采用有效的地面和井下综合防治措施后，已分别于2017年11月22日和2018年1月24日安全通过。地面通井和井下揭露显示，马探31号无套管（可称裸眼井），马探30号有套管，但由于成井久远，套管破损严重。揭井时发现这两口井内都存在大量高压油层气，其对综采工作面安全生产构成巨大威胁[3]。

油层气[4]是指来自于含油地层的有害气体的总称，其主要成分为CH4，其余成分中H2S含量相对较高[5]。石油井将含油层和煤层导通后，井内的高压油层气会向临近煤层扩散渗流[6-7]，从而使井周围形成一定范围的油层气富聚区，严重制约综采工作面安全生产。关于煤层中气体扩散渗流规律，杨其銮等人依据扩散传质理论和试验研究，提出极限煤粒假说并建立了煤粒瓦斯扩散的微分方程[8-9];鲜学福等人用真实气体状态方程推导出了煤层瓦斯扩散渗流的微分方程[10-11];周世宁等人在煤层瓦斯流动符合线性达西定律等假设的基础上，建立了煤层瓦斯渗流的物理数学模型[12-13];赵阳升等人根据弹性变形和瓦斯渗流理论建立了煤体瓦斯耦合物理数学模型[14-15];吴世跃等人根据气固耦合理论建立了煤层气与煤层耦合运动的物理模型[16-18]。以上多为单一研究煤层瓦斯的纯扩散或纯渗流过程，对组分与瓦斯存在较大差异的油层气的扩散-渗流耦合运动的研究较少。

此文章在前人的基础上，综合分析煤层内游離、吸附态油层气运移机理，结合油层气流动过程中煤体渗透率的动态变化及煤体变形机理，建立了三维空间下的煤层内油层气流动模型，针对油层气特殊组分和宁夏双马煤矿煤层具体赋存条件，运用COMSOL Multiphysics软件对无套管或套管已破损的废弃石油井内的高压油层气向临近煤层扩散渗流过程进行模拟研究，掌握流动过程中油层气压力的变化规律，确定某一段时间后油层气对临近煤层的影响范围，得到井内不同压力值下具有时间效应的影响半径，为受废弃石油井威胁煤层的安全高效开采提供理论依据。

1 扩散-渗流耦合模型

1.1 模型假设

油层气的主要成分与煤层气相差不大，油层气在煤层中的扩散-渗流模型可借鉴煤层气的建立。参照文献[19]，可基于扩散渗流并存的孔隙裂隙二重介质假设建立煤层中油层气流动的扩散-渗流耦合模型。该假设认为：煤层是由含分子尺度孔隙的煤粒（基质）骨架及煤粒骨架之间间隙组成的孔隙裂隙二重介质，煤层内油层气有吸附态和游离态2种赋存状态，吸附态油层气主要吸附在基质的孔隙内，而游离态油层气则主要赋存在煤粒骨架之间的裂隙中，结构如图1所示。该模型还满足如下假设

1）煤层为均匀各项同性的线弹性体，并且其变形属于微小变形。

2）石油井内高压作用下，油层气在裂隙中渗流，裂隙内油层气浓度升高后，其在裂隙与孔隙浓度梯度下向基质孔隙吸附扩散，其对孔隙系统是流入，对裂隙系统为流出。

3）油层气为理想气体，其流动为等温流动，扩散过程遵循菲克定律，渗流过程遵循达西定律。

2.3 定解条件

初始条件：实测4-1煤层初始气体压力为0.21 MPa;通井时实测平均油层气压力为5 MPa.

边界条件：模型上部为压力边界，承受上覆岩层重力;下部为固定约束边界;四周为辊支撑约束边界;内部加载自重载荷;上下面和四周均为零通量不透气边界;石油井壁为裂隙流边界。

3 计算结果及分析

3.1 油层气压力变化规律

扩散渗流过程中煤层内油层气压力变化如图4所示，主要在3个时间段内表现出不同变化规律。第1个时间段：从t=0 d到t=30 d，压力变化不大，尤其是在t=0 d到t=15 d内，并且从图5可知，近处煤层压力不连续。第2个时间段：从t=30 d到t=60 d，压力增加的很多。第3个时间段：从t=60 d到t=120 d，压力继续增加，但增加量相对前一时段较少，并且在t=90 d到t=120 d内增加幅度也小于在t=60 d到t=90 d内增加幅度。

有效应力是煤体骨架之间传递的力，它与流体间传递的孔隙压力共同抵抗煤体受到的外力[22]。一般来说，煤体内孔隙压力越小有效应力就越大，骨架承受的力越大压缩就越严重，骨架之间的裂隙越小，从而煤体孔隙率、渗透率越小，因此孔隙压力和有效应力是影响煤体渗透率的重要因素[23-24]。在有效应力和吸附作用下，煤体会发生2种变形[16]：一是骨架变形（整体变形）;二是吸附变形。

煤体内较大的有效应力和较小的孔隙压力共同抵抗上覆岩层重力，因此煤层整体处于受压状态，煤体发生骨架压缩变形，煤体渗透率很小。从t=0 d到t=30 d，由石油井流向煤层的油层气会使石油井近处煤层的孔隙压力增大，有效应力减小，但由于处于扩散渗流初期，油层气压力较小，在较低压力水平下，基质内向吸附膨胀变形（也称内膨胀）占主导作用，油层气不断向基质孔隙内吸附，吸附膨胀会减小基质孔隙体积并挤压裂隙空间，煤体孔隙率减小，所以此时虽然有效应力减小，但煤体渗透率不但不增加，反而也进一步减小，从而导致油层气压力在该时间段内增加得缓慢，而同时在油层气压力不断缓慢增加的过程中，压缩的裂隙在裂隙内油层气压力作用下逐渐被导通，煤体内裂隙由不连续状态逐渐转变为连续状态，油层气压力相对之前增加变快，所以在t=0 d到t=15 d内增加的很少，而在t=15 d到t=30 d内增加的相对多一点。

从t=30 d到t=60 d，石油井近处煤层内裂隙基本被导通，随着井内油层气源源不断地向煤层内补充，煤层内孔隙压力越来越大，有效应力越来越小，在较高压力水平下，基质外向吸附膨胀变形（也称外膨胀）占主导作用，此时上覆岩层对内有较高气体压力的煤层的位移约束作用减小，吸附膨胀会在减小基质孔隙体积的同时又扩张裂隙空间，煤体内裂隙逐渐被打开，孔隙率、渗透率逐渐变大，油层气流动加快，压力增加的多。

从t=60 d到t=120 d，随着时间的推移，煤体外向吸附膨胀变形逐渐趋于稳定，同时石油井近处高压油层气在向远处流动时一直伴随着能量损耗，受到的阻力越来越大，所以煤体孔隙率、渗透率增加的越来越少，油层气压力增加量相对前一时段较小，并且在t=90 d到t=120 d内增加量同样也小于在t=60 d到t=90 d内增加量。

为了能更直观地看出流动过程中某一位置处油层气压力变化规律，在模型y=0 m，z=-2 m，x={10，20，30，40，50} m处各取一个观测点，在数值模拟过程中分别记录不同时间下各点处油层气压力值，如图6所示。

图6中各处煤层内油层气压力的增加速度都先加快后减慢，油层气压力在流动后期都缓慢增加直至趋于稳定;距离废弃石油井越近，油层气压力短时间内变化越快，压力也越大。

3.2 影响半径变化规律

由于废弃石油井内高压气体在煤层中的扩散渗流是一个较为漫长的过程，随着时间的推移，其影响范围会逐渐扩大，因此对原几何模型进行修改：将几何模型扩大至2 000 m×2 000 m×4 m.

瓦斯压力超过0.74 MPa的煤层为突出危险煤层，文章取煤层内油层气压力0.74 MPa为扩散渗流影响范围临界压力值，对不同油层气压力下、不同时间内影响半径的变化规律进行模拟研究。图7为模拟过程中影响半径变化图，图中井内油层气压力为10 MPa，当井内油层气向煤层流动至第10年时，影响半径为220 m.油层气分别在5，10，15和20 MPa井内恒定压力下流动过程中具有时间效应的影响半径，如图8所示。

图8中不同恒定压力下的扩散渗流影响半径随着时间的推移都在增大，并且在扩散渗流前期增大的速率大于后期增大的速率;井内气体压力越大，影响半径增大的速率越大，同一时刻影响半径值也越大;第30 a时，油层气在4种压力下的影响半径分别为188，361，465，537 m.

考虑到实际上许多废弃石油井内的高压力油层气都是随时间逐渐聚集起来的，油层气在这些变化压力下的流动与上述流动有较大差异，其影响半径也会有所不同，因此继续对变化压力下油层气流动进行模拟研究。假设井内压力从标准大气压随时间线性增加，30 a后压力分别达到5，10，15，20 MPa，即4种随时间变化的压力为

pi=kit+p0，i=1，2，3，4

式中 pi为第t年时井内气体压力，MPa;ki为线性系数，4种情况下取值依次为16，13，12，23.

4种变化压力下油层气流动过程中的扩散渗流影响半径，如图9所示。

图9中影响半径增加趋势大致与图8相似，不同的是每条变化曲线都呈现出不同程度的“S”型，并且井内气体压力变化越快，其对应曲线“S”程度越高。原因是恒定压力下，前期煤层内油层气压力很大，影响半径增加速度大，但随着油层气流动距离的加大，油层气流动受到的阻力不断增加，影响半径增加速度一直减小，所以曲线斜率一直在减小。变化压力下，井内油层气压力是随时间逐渐增加的，前期油层气压力很小，影响半径增加得慢，但随着时间的推移，井内油层气不断聚集，影响半径增加变快，曲线斜率增加，一段时间后，远距离流动的油层气受到的阻力不断增大，远处油层气对该线性增加下压力的增加已不再那么敏感，影响半径增加速度变慢，曲线斜率又开始减小，故整条曲线呈现“S”型。井内气体压力变化越快，前期油层气流动越快，影响半径增加越快，对应曲线的“S”型程度越高。

第30 a时，图9中4种变化压力下的影响半径依次为120，239，310，362 m，对比图8可以看出，尽管此时井内油层气压力相同，但两者影响半径之间存在的差异较大。由于许多废弃石油井内油层气压力不断增加，所以图9更能真实地反映影响半径的实际变化规律。

4 现场实测比较

在采掘工作面靠近废弃石油井过程中，双马煤矿进行了实测。已揭露的马探31号石油井已废弃30余年，地面通井时测得其井内气体最大压力13.992 MPa;在I0104105综采面巷道内靠近马探31号石油井布置了探孔，根据M2305，M2409和M3305等探孔实测，距离马探30号越近，煤层内气体压力越高;根据I0104105综采工作面实测，随着综采面靠近马探31号，煤层内气体压力逐渐升高，实测距离该井300 m的煤层气体压力达到0.74 MPa，现场据此划分出马探31号石油井明显影响范围为300 m，这与数值计算结果中压力变化规律及第30年井内油层气压力线性增加至15 MPa时的310 m油层气扩散渗流影响半径基本吻合，说明建立的模型合理，能反映出废弃石油井周围煤层内油层气的扩散渗流规律。

5 结 论

1）流动过程中，石油井附近煤层各处油层气压力的增加速度都先加快后变慢，油层气压力在流动后期都缓慢增加直至趋于稳定，距离井越近，油层气压力短时间内变化越快，压力也越大。

2）油层气影响半径随着时间逐渐增加，最后都趋于稳定，井内压力恒定时，影响半径的增加速度一直在衰减，井内压力线性增加时，影响半径的增加速度先增后减。

3）第30年时，当井内油层气压力线性增加至15 MPa时，油层气影响半径为310 m，数值模拟结果与现场实测结果基本吻合。

参考文献（References）：

[1] 张 健.子长矿区煤油地质特征及合理开发对策研究[D].西安：西安科技大学，2013.

ZHANG Jian.The research on the coal and petroleum geological features and the reasonable development countermeasures in Zichang mining area[D].Xi’an：Xi’an University of Science and Technology，2013.

[2] 張 荃.黄陵油型气涌出矿井CH4扩散规律数值模拟[D].西安：西安科技大学，2015.

ZHANG Quan.Numerical simulation of CH4 diffusion regularities in the coalmine of Huangling with oiltype gas emission[D].Xi’an：Xi’an University of Science and Technology，2015.

[3]孙四清.煤油气共存矿井掘进工作面底板油型气涌出机理探讨[J].矿业安全与环保，2017，44（4）：90-94.

SUN Siqing.Discussion on oiltype gas emission mechanism from heading face floor in coaloilgas coexisting mine[J].Mining Safety & Environmental Protection，2017，44（4）：90-94.

[4]张俭让，张 荃，董丁稳，等.掘进巷道油型气扩散规律数值模拟[J].西安科技大学学报，2015，35（3）：307-312.

ZHANG Jianrang，ZHANG Quan，DONG Dingwen，et al.Numerical simulation of oiltype gas diffusion regularities in roadway excavation[J].Journal of Xi’an University of Science and Technology，2015，35（3）：307-312.

[5]朱扬明，王积宝，郝 芳，等.川东宣汉地区天然气地球化学特征及成因[J].地质科学，2008（3）：518-532.

ZHU Yangming，WANG Jibao，HAO Fang，et al.Geochemical characteristics and origin of natural gases from Xuanhan area in Eastern Sichuan[J].Chinese Journal of Geology，2008（3）：518-532.

[6]金之钧，张发强.油气运移研究现状及主要进展[J].石油与天然气地质，2005（3）：263-270.

JIN Zhijun，ZHANG Faqiang.Current status and main progress of oil and gas migration[J].Oil & Gas Geology，2005（3）：263-270.

[7]李明诚.石油与天然气运移研究综述[J].石油勘探与开发，2000（4）：3-10，109-117.

LI Mingcheng.Summary of research on oil and gas migration[J].Petroleum Exploration and Development，2000（4）：3-10，109-117.

[8]杨其銮.煤屑瓦斯扩散理论及其应用[J].煤炭学报，1986，11（3）：67-70.

YANG Qiluan.Diffusion theory and application of methane in coal seam[J].Journal of China Coal Society，1986，11（3）：67-70.

[9]Wu D M，Wang H F，Ge C G，et al.Research on forced gas draining from coal seams by surface well drilling[J].Mining Science and Technology，2011，21（2）：229-232.

[10]鲜学福，张广祥，谭学术，等.煤层瓦斯运移数学模型[J].重庆大学学报，1994，17（4）：53-64.

XIAN Xuefu，ZHANG Guangxiang，TAN Xueshu，et al.Mathematical model of methane transport in coal seam[J].Journal of Chongqing University，1994，17（4）：53-64.

[11]李 波.受载含瓦斯煤渗流特性及其应用研究[D].北京：中国矿业大学（北京），2013.

LI Bo.Seepage characteristics of loaded coal containing gas and its application[D].Beijing：China University of Mining & Technology（Beijing），2013.

[12]周世宁.瓦斯在煤层中流动的机理[J].煤炭学报，1990，15（1）：61-74.

ZHOU Shining.Mechanism of methane flow in coal seam[J].Journal of China Coal Society，1990，15（1）：61-74.

[13]周世宁，孙辑正.煤层瓦斯流动理论及其应用[J].煤炭学报，1965，2（1）：24-36.

ZHOU Shining，SUN Jizheng.Theory and application of methane flow in coal seam[J].Journal of China Coal Society，1965，2（1）：24-36.

[14]趙阳升.煤体-瓦斯耦合数学模型及数值解法[J].岩石力学与工程学报，1994，13（3）：229-239.

ZHAO Yangsheng.Coupled mathematical model on coal massmethane and its numerical method[J].Chinese Journal of Rock Mechanics and Engineering，1994，13（3）：229-239.

[15]石军太，李相方，徐兵祥，等.煤层气解吸扩散渗流模型研究进展[J].中国科学：物理学天文学力学，2013，43（12）：1548-1557.

SHI Juntai，LI Xiangfang，XU Bingxiang，et al.Review on desorptiondiffusionflow model of coalbed methane[J].Sci SinPhys Mech Astron，2013，43（15）：1548-1557.

[16]吴世跃.煤层气与煤层耦合运动理论及其应用的研究[D].沈阳：东北大学，2006.

WU Shiyue.Research of methanecoalbed coupling movement theory and its application[D].Shenyang：Northeastern University，2006.

[17]郭 平，曹树刚，张遵国，等.煤体吸附膨胀变形模型理论研究[J].岩石力学，2014，35（12）：3467-3472.

GUO Ping，CAO Shugang，ZHANG Zunguo，et al.Theoretical study of deformation model of coal swell induced by gas adsorption[J].Rock and Soil Mechanics，2014，35（12）：3467-3472.

[18]Li B，Wei J P，Li P，Numerical simulation on gas drainage of boreholes in coal seam based on gassolid coupling model[J].Computer Modelling & New Technologies，2014，18（12A）：1-8.

[19]林柏泉，刘 厅，杨 威.基于动态扩散的煤层多场耦合模型建立及应用[J].中国矿业大学学报，2018，47（1）：32-39.

LIN Boquan，LIU Ting，YANG Wei.Solidgas coupling model for coalseams based on dynamic diffusion and its application[J].Journal of China University of Mining and Technology，2018，47（1）：32-39.

[20]張钧祥.基于扩散-渗流机理瓦斯抽采三维模拟研究[J].地下空间与工程学报，2018，2（1）：110-115.

ZHANG Junxiang.Research on 3D numerical simulation of gas extraction based on diffusionseepage mechanism[J].Chinese Journal of Underground Space and Engineering，2018，2（1）：110-115.

[21]魏凯丰，宋少英，张作群.天然气混合气体粘度和雷诺数计算研究[J].计量学报，2008（3）：248-250.

WEI Kaifeng，SONG Shaoying，ZHANG Zuoqun.Research on calculation of natural mix gas viscosity and reynolds number[J].Acta Metrologica Sinica，2008（3）：248-250.

[22]Jasinge D，Ranjith P G，Choi S K.Effects of effective stress changes on permeability of latrobe valley brown coal[J].Fuel，2011，90（3）：1292-1300.

[23]张景飞，郭 倩，朱同功，等.多场耦合下煤岩渗透率演化规律——以平煤十矿为例[J].西安科技大学学报，2018，38（5）：713-720.

ZHANG Jingfei，GUO Qian，ZHU Tonggong，et al.Evolution of coal rock permeability with multiphysics coupling：taking Pingmei No.10 mine as an example[J].Journal of Xi’an University of Science and Technology，2018，38（5）：713-720.

[24]孙光中，荆永滨.加卸载作用下构造煤峰值前渗透率变化特征试验研究[J].西安科技大学学报，2016，36（2）：193-199.

SUN Guangzhong，JING Yongbin.Experimental study on prepeak permeability of coal samples under axial stress loading and unloading[J].Journal of Xi’an University of Science and Technology，2016，36（2）：193-199.