高级检索

改性褐煤对酸性矿山废水中Fe2+、Mn2+的吸附特性研究

狄军贞, 张思怡, 杨逾, 梁冰, 王显军, 孟凡康, 李增新, 郭建志, 李拓达, 张伟

狄军贞, 张思怡, 杨逾, 梁冰, 王显军, 孟凡康, 李增新, 郭建志, 李拓达, 张伟. 改性褐煤对酸性矿山废水中Fe2+、Mn2+的吸附特性研究[J]. 煤炭科学技术, 2022, 50(11): 243-250.
引用本文: 狄军贞, 张思怡, 杨逾, 梁冰, 王显军, 孟凡康, 李增新, 郭建志, 李拓达, 张伟. 改性褐煤对酸性矿山废水中Fe2+、Mn2+的吸附特性研究[J]. 煤炭科学技术, 2022, 50(11): 243-250.
DI Junzhen, ZHANG Siyi, YANG Yu, LIANG Bing, WANG Xianjun, MENG Fankang, LI Zengxin, GUO Jianzhi, LI Tuoda, ZHANG Wei. Study on adsorption characteristics of modified lignite for Fe2+, Mn2+ in acid mine drainage[J]. COAL SCIENCE AND TECHNOLOGY, 2022, 50(11): 243-250.
Citation: DI Junzhen, ZHANG Siyi, YANG Yu, LIANG Bing, WANG Xianjun, MENG Fankang, LI Zengxin, GUO Jianzhi, LI Tuoda, ZHANG Wei. Study on adsorption characteristics of modified lignite for Fe2+, Mn2+ in acid mine drainage[J]. COAL SCIENCE AND TECHNOLOGY, 2022, 50(11): 243-250.

改性褐煤对酸性矿山废水中Fe2+、Mn2+的吸附特性研究

基金项目: 

国家自然科学基金资助项目(41672247,41102157);辽宁省“兴辽英才”青年拔尖人才计划支持资助项目(XLYC1807159)

Study on adsorption characteristics of modified lignite for Fe2+, Mn2+ in acid mine drainage

Funds: 

National Natural Science Foundation of China (4167224741102157); Liaoning Province "Xingliao Talents" Youth Top notch Talents Program Support and Funding Project (XLYC1807159)

  • 摘要: 针对酸性矿山废水(Acid Mine Drainage, AMD)中Fe2+、Mn2+含量高,褐煤吸附饱和量有限等问题,以褐煤为试验材料,利用扫描电镜(SEM)、傅里叶红外光谱仪(FT-IR)、X射线衍射仪(XRD)等手段分析褐煤的吸附机理。同时,采用静态烧杯试验,并基于NaCl浸泡改性、超声波辅助NaCl改性褐煤的方法以及吸附等温线、吸附动力学原理,研究褐煤及改性褐煤吸附Fe2+、Mn2+特性和机理。结果表明:褐煤的表面孔隙结构发达,褐煤主要由高岭石、方解石、菱铁矿、石英等矿物组成,含有丰富的—OH、—C〖KG-*5〗〖FY=,1〗O、—CHO等活性基团。超声波辅助NaCl改性褐煤对Fe2+、Mn2+的吸附效果最佳,平衡时吸附饱和量分别为6.775 0 mg/g和8.092 5 mg/g,较未改性褐煤提高了28.44%、22.47%。褐煤及其改性褐煤对Fe2+、Mn2+的吸附过程均符合Langmuir等温模型。研究表明,褐煤及改性褐煤表面均匀并且为单分子层吸附,被吸附的Fe2+、Mn2+之间没有相互作用力。褐煤及改性褐煤的吸附能力大小为:超声波辅助NaCl改性褐煤>NaCl浸泡改性褐煤>褐煤。其中,吸附效果最好的超声波辅助NaCl〖JP2〗改性褐煤对Fe2+、Mn2+的Langmuir拟合方程分别为y=0.069 95x+0.139 24、y=0.070 73x+0.425 16。〖JP〗褐煤对AMD中Fe2+的吸附过程符合准二级动力学模型,说明褐煤对Fe2+的吸附机理以离子交换作用为主;NaCl浸泡改性褐煤、超声波辅助NaCl改性褐煤对AMD中Fe2+的吸附现象符合颗粒内扩散模型,表明该吸附过程是由膜扩散和颗粒内扩散共同控制。褐煤及改性褐煤对AMD中Mn2+的吸附均符合准二级动力学模型,说明褐煤及改性褐煤对AMD中Mn2+的吸附以化学吸附为主。其中吸附效果最好的超声波辅助NaCl改性褐煤对Fe2+的颗粒内扩散模型方程为y=0.480 24x+0.453 12;对Mn2+的准二级动力学拟合方程为y=0.103 44x+3.415 11。
    Abstract: In view of the high content of Fe2+ and Mn2+ in acid mine drainage (AMD) and the limited adsorption saturation of lignite, the adsorption mechanism of lignite was analyzed by scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FT-IR) and X-ray diffraction (XRD). At the same time, the static beaker experiment was used to study the characteristics and mechanism of adsorption of Fe2+ and Mn2+ by lignite and modified lignite based on the methods of NaCl soaking modification, ultrasonic assisted NaCl modification, adsorption isotherm and adsorption kinetics. The results show that the surface pore structure of lignite is developed. Lignite is mainly composed of kaolinite, calcite, siderite, quartz and other minerals, and rich in —OH, —C〖KG-*5〗〖FY=,1〗O, —CHO and other active groups. The results show that the adsorption of Fe2+ and Mn2+ by ultrasonic assisted NaCl modified lignite is the best, and the saturated adsorption capacity is 6.775 0 mg/g and 8.092 5 mg/g respectively, which is 28.44% and 22.47% higher than that of unmodified lignite. The adsorption process of Fe2+ and Mn2+ on lignite and its modified lignite conformed to Langmuir isotherm model. The results show that the surface of lignite and modified lignite is homogeneous and monolayer adsorption, and there is no interaction between Fe2+ and Mn2+ adsorbed. The adsorption capacity of lignite and modified lignite is as follows: ultrasonic assisted NaCl modified lignite>NaCl soaked modified lignite>lignite. 〖JP2〗Among them, the Langmuir fitting equations of Fe2+ and Mn2+ by ultrasonic assisted NaCl modified lignite were y=0.069 95x+ 0.139 24 〖JP〗and y=0.070 73x+0.425 16, respectively. The adsorption process of lignite for Fe2+ in AMD accords with the pseudo second order kinetic model, which indicates that the adsorption mechanism of lignite for Fe2+ is mainly ion exchange; the adsorption of Fe2+ in AMD by NaCl soaking modified lignite and ultrasonic assisted NaCl modified lignite conforms to the intra particle diffusion model, which indicates that the adsorption process is controlled by both membrane diffusion and intra particle diffusion. The adsorption of Mn2+ in AMD by lignite and modified lignite is consistent with the pseudo second order kinetic model, which indicates that the adsorption of Mn2+ in AMD by lignite and modified lignite is mainly chemical adsorption. Among them, the best adsorption effect of ultrasonic assisted NaCl modified lignite is y=0.480 24x+0.453 12 for Fe2+, and y=0.103 44x+3.415 11 for Mn2+.
  • 排水降压是中−浅部煤层气开发的主流方式。煤层气井排水−降压−产气的有效性,决定着气井是否能够获得高产[1-2],存在2个关键环节,即排水能否诱导降压与降压能否诱导解吸产气,前者涉及排采水源、煤储层及围岩含水性、地下水动力条件等问题[3-5],能够诱导储层降压的水源为弹性水,也称有效水,反之为水侵水,也称无效水[6-7];后者涉及资源量与资源丰度、含气饱和度、解吸效率、煤储层渗透性与压裂改造效果等问题[8-12]。针对低产低效井,如何将这2个环节有效剥离,是寻找产气受限关键因素,制定针对性提产改造措施的前提[13-15]。另一方面,现场丰富的排采数据隐含着排水、降压、产气等有效性信息,基于排采数据构建一种相对简易、操作性强,且具有较强实用性的煤层气井排采效果诊断与水源识别方法,并获取气藏相关属性参数,具有重要的实践指导意义。

    黔西织金区块煤层气勘探开发近年来取得重要进展,在笔者长期的跟踪研究中,基于产出水地球化学特征,提出表层水干扰深刻制约着该区煤层气高效开发,导致气井产水量大而产气量低[16-18]。原因在于,黔西地区石灰岩大面积出露,发育喀斯特地貌,且降雨量丰富,地表水与浅部地下水循环交替强烈,我国南方岩溶地区降水入渗系数可达0.8以上,导致上部煤层气系统排采容易受到浅层地下水影响,限制储层降压效果,从而降低开发效率[19-20]

    针对上述问题,笔者前期建立了基于生产特征曲线的外部水体侵入(水侵)判识方法,提出了煤系弹性产水能力极限值,并在黔西织金区块取得了应用成效[21]。经过实践检验,该方法可在3个方面进行优化:一是引入了生产压差这一在煤层气井排采控制重要的参数,在保持井口套压相对稳定的前提下,实际计算中若直接以井底流压为分析对象将更加简便;二是无法直接计算某些气藏参数,也未对产水量构成进行劈分;三是无法剥离排水降压与降压产气的有效性,忽略了产气量这一重要参数。优化的基本思路为:将生产压差替换为井底流压,并借鉴油井生产指示曲线分析方法[22],引入累计产气量,根据典型参数相关性曲线分析排水、降压、产气效果及其联动关系与约束环节。

    为此,笔者立足黔西织金区块煤层气井排采数据,以煤层气井排水—降压—产气有效性分析为目的,基于累计产水量与井底流压、累计产气量与井底流压、以及累计产水量与累计产气量3条关系曲线,对排水−降压有效性与降压−产气有效性进行剥离,寻找产能主要制约因素;提取曲线中隐含的气藏属性参数,包括启抽压力、弹性产水指数、井控动态储量等;探索弹性水(有效水)与水侵水(无效水)定量劈分方法,为促进煤层气高效开发提供技术支持。

    在常规排采曲线动态分析的基础上,提取井底流压、累计产水量、累计产气量3类数据,绘制两两关系曲线。累计产水量与井底流压关系曲线可以反映排水降压效率,累计产气量与井底流压关系曲线可以反映降压产气效率,累计产水量与累计产气量关系曲线可反映排水产气效率,是前两方面效果的综合体现。本文将上述3条曲线统称为排采有效性分析曲线,通过其耦合分析,为揭示影响产气效率的关键因素提供依据。

    对曲线中的近似线性段进行拟合是重要的分析手段,拟合直线的斜率、截距等信息可指示若干气藏属性。选段拟合的基本原则是线性、连续、稳定:拟合直线R2在0.6以上,能延伸一定长度(实质为较长的排采时间,至少60 d),无明显人为操作干扰,尤其不能有频繁的停井和修井作业,曲线较为平滑,代表稳定的排采阶段,保障压降规律分析效果。

    根据累计产水量与井底流压关系曲线,可获取:① 排采初期的拟合直线:其Y轴截距为初始井底流压,称为启抽压力,当压裂未沟通含水层,可视为原始储层压力,若沟通开放型含水层,会导致原始压力释放;斜率为单位产水量的降压幅度,代表排水降压效率;斜率的倒数为单位井底流压降幅的产水量,代表弹性产水指数,其值越高,储层含水性越强,水侵亦会导致该值升高;② 对于下凹型曲线:排采后期大量产气导致产水受限,使曲线呈下凹型,代表未发生水侵;③ 对于上凹型曲线(发生水侵):通过确定曲线拐点将曲线划分为2段,前面下降段代表弹性产水阶段,后面水平段代表水侵产水阶段。前期研究借鉴未饱和弹性驱动油藏生产指示曲线分析方法,以生产初期直线段斜率为弹性产液指数,代表单位压降的产液量,该值越高,水体弹性能越大。据此获得织金区块煤系极限弹性产水指数为200 m3/MPa,可区分含煤地层内源水与外源水2类水源。若下降段的弹性产水指数超过此极限值,说明发生水侵[19]。这种情况下,超过部分可视为水侵量,即200 m3/MPa这一极限值与实际弹性产水指数之比,为弹性水在弹性产水阶段总产水量中所占的比例,其余为水侵量比例;④ 水平段为一理想状态,实际情况多为近水平,若微弱上升,则代表完全水侵(产水的同时井底流压上升),若微弱下降则指示有一小部分弹性水产出,此时弹性水量与水侵量的劈分方法与下降段相同。

    基于累计产气量与井底流压关系曲线,拟合直线的斜率之倒数为单位井底流压降幅对应的产气量,X轴截距为井控动态储量。累计产水量与累计产气量关系曲线中,直线斜率能够反映排水产气的效率,即单位产水量对应的产气量。

    本文选择织金区块6口典型煤层气排采井开展分析,搜集产水、产气、流压、套压等生产数据,以及层位组合、埋深、跨度、厚度等产层信息(表1图1)。

    表  1  煤层气井数据信息
    Table  1.  CBM well data information
    井号 产层组合(煤层编号) 产层埋深/
    m
    最大跨度/
    m
    产层厚度/
    m
    平均日产水量 /
    (m3·d−1)
    平均日产气量 /
    (m3·d−1)
    峰值日产气量 /
    (m3·d−1)
    Z−1 9、10、炭质泥岩段、23、24 214~376 162 14.1 3.81 305 428
    Z−2−1 20、23 407~432 25 3.2 0.56 1650 2772
    Z−2−2 6−1、7、8、10、12、14、17 240~389 149 6.8 3.51 387 770
    Z−3 14、16 697~737 40 4.7 0.80 817 1137
    Z−4 6、7、17、20、23、27、30 283~531 248 11.7 1.98 798 2465
    Z−5 16、17、20、23、27 355~426 71 7.6 0.85 1433 2806
    Z−S1 8、9、10、14、15、16、17、19、
    20、21、22、23、26、27
    270~453 183 18.4 23.08 5 68
      注:Z−2−1代表Z−2井第1生产阶段;Z−2−2代表Z−2井第2生产阶段,下同。
    下载: 导出CSV 
    | 显示表格
    图  1  织金区块地质特征
    Figure  1.  Geological setting of the Zhijin Block

    黔西地区上二叠统海陆交互相含煤地层普遍发育多层叠置含煤层气系统,系指煤层(组)之间由于缺乏流体联系与能量交换,形成具有独立流体压力状态、水动力条件与含气性特征的含煤层气单元[23],不同含气单元垂向叠置,客观上造成了储层属性、流体能量与产气潜力的垂向波动,进而易诱发合采兼容性问题[24-25]。织金区块位于贵州省西部,构造上位于扬子板块黔中隆起,以岩脚向斜为主体构造单元,包括珠藏、阿弓、三塘、水公河、比德等含煤次向斜,具有煤层数量多、累计厚度大、含气量高、埋藏浅、煤体结构较完整、割理裂隙较发育、向斜控气等地面开发优势条件[26-27]。Z−2井、Z−4井、Z−5井与Z−S1井位于珠藏次向斜,分布较为集中,井间距500~600 m;Z−3井位于三塘次向斜,距离珠藏次向斜集中分布的煤层气井直线距离约12 km;Z−1井位于岩脚向斜以北的黔西向斜,距离珠藏次向斜煤层气井直线距离约40 km。6口井未受断层、陷落柱等影响,均位于NE向构造单元,煤体结构较为完整,地质条件相近,变量主要为排采层位及其埋深。构造煤主要发育于NW向构造,以比德次向斜为代表,导致煤层气开发效果不理想。珠藏次向斜位于织金区块东部,已被证实是织金区块乃至黔西地区最具煤层气开发潜力的构造单元,也是目前的主力建产区。

    根据层序地层格架及其控制下的流体赋存与能量分布特点,将该区上二叠统含煤地层划分为3套垂向叠置的含煤层气系统(简称含气系统),包括上部1~9号煤层含气系统,中部10~21号煤层含气系统,以及下部22~35号煤层含气系统[28]图1)。上部含气系统由于海侵增强发育K1~K4共4套石灰岩标志层,浅埋藏条件下岩溶发育,形成岩溶裂隙含水层,构成6~9号煤层开采的直接充水含水层[29]。抽水试验显示,当抽水层段包括灰岩层或以灰岩层为边界时,多显示较好的富水与导水条件[28]。例如,珠藏次向斜抽水试验钻孔6−4孔对长兴组顶部−6号煤底,6号煤底−16号煤底,16号煤底−龙潭组底部3个层段分别开展了抽水试验,获得的单位涌水量分别为0.01140.0038、0.0009 L/(s·m),渗透系数分别为0.05710.01020.0011 m/d,反映上部层段具有最强的水动力条件,源于K1、K2与K3灰岩含水层的贡献。根据水文柱状结构,上部含气系统含水层厚度约占地层总厚度的一半,强含水层主要为灰岩层,弱含水层主要为砂岩层,隔水层岩性主要为泥质岩类与煤层,三者互层产出,单层厚度多<10 m。另外,6号煤的顶板为K3灰岩,7号煤顶板为K4灰岩,导致煤层气开发过程中发生水侵的概率极高(图1)。

    开发实践表明,中部与下部含气系统具有相似的流体状态与产能特点,合采兼容性好[30]。因此,本文将传统含气系统概念向开发领域延伸,进一步提出产气系统概念,指具有相近流体动力条件、开发动态及生产规律的含气系统及其组合,鲜明特征是不同产气系统具有差异显著的气水产出特征与能量驱动机制,合采兼容性差。据此,将1~9号煤层定义为易受浅部含水层影响的上部产气系统,10~35号煤层为水文地质条件相对封闭的中下部产气系统。

    煤层含气量高但原位渗透率低,上部系统煤层渗透率主体位于0.01×10−3~0.1×10−3 μm2水平,含气量主体介于10~17 m3/t;中下部系统煤层渗透率主体介于0.0001×10−3~0.01×10−3 μm2,含气量主体介于14~22 m3/t[26-28,30]。排采总体遵循连续、缓慢、稳定原则,排水阶段液面降速<3 m/d,产气阶段液面降速<1 m/d,在不吐砂吐粉的情况下,见气前可适当加快排采速率[31-32]。采用活性水压裂液,配方为清水+0.2%助排剂+0.05%杀菌剂。施工加砂强度介于10~15 m3/m,加液强度>150 m3/m。破裂压力介于10~45 MPa,具有随埋深增加而增大的趋势,但储层改造效果与破裂压力呈反比[33]。研究区煤储层总体低渗,煤体结构接近,改造条件与改造工艺相似,本文相关分析均是以各井煤层物性参数(改造后)无明显差异为前提的。

    Z−2井分为2个生产阶段,第1阶段排采20号、23号煤,第2阶段在封堵下部煤层之后,排采上中部的6、7、8、10、12、14、17共7个煤层(图2)。

    图  2  煤层气井排采曲线
    Figure  2.  Production curves of CBM wells

    第1阶段(0~第263天):排采时间263 d,见气时间57 d,最高日产气量2772 m3,平均日产气量1650 m3,最高日产水量4.08 m3,平均日产水量0.56 m3,见气前平均日产水量1.86 m3,见气后平均日产水量0.20 m3。产气效果较好,见气后日产气量迅速提高,日产水量则迅速降低,产气量长期保持在1000 m3/d以上(图2a)。

    第2阶段(第266天~第814天):排采时间548 d,见气时间45 d,最高日产气量770 m3,平均日产气量387 m3,最高日产水量6.11 m3,平均日产水量3.51 m3。见气前平均日产水量1.12 m3,见气后平均日产水量3.72 m3,见气后的产水量高于见气前,属于异常状况。该阶段产能效果不佳,相比第1阶段,呈现产水量高而产气量低的特征(图2a)。产出水地球化学特征检测表明,该阶段产出水微量元素浓度与典型煤层水相比显著降低,表明排采过程发生了水侵,降低了产气效率[34]。第2阶段在排采煤层数量和产层总厚度增加的前提下,产气效果反而变差,源于产层组上返排采后遭遇浅部含水层干扰,产水量明显超过煤层实际供液能力[21],压降无法有效扩展。

    Z−1井合采9、10、炭质泥岩段、23号、24号煤层,排采时间543 d,最高日产气量428 m3,平均日产气量305 m3,最高日产水量5.88 m3,平均日产水量3.81 m3,产能效率低。排采至第59天见气,见气前平均日产水量2.83 m3,见气后平均日产水量3.93 m3,见气后产水量高于见气前。由于层位浅,跨度大,跨越2个产气系统,易受水侵影响(图2b)。

    Z−3井合采14号与16号煤层,排采时间272 d,最高日产气量1137 m3,平均日产气量817 m3,最高日产水量2.03 m3,平均日产水量0.80 m3。见气时间97 d,见气前平均日产水量1.7 m3,见气后产气量逐步上升,产水量逐步下降,平均仅0.29 m3/d。本井产气情况稳定,无明显水侵特征,产气量突破了1000 m3/d且持续超过3个月的储量报审标准,昭示了区内煤层气开发潜力(图2c)。

    Z−4井合采6、7、17、20、23、27、30号共7个煤层,排采时间351 d。最高日产气量2465 m3,平均日产气量798 m3,最高日产水量21.22 m3,出现于排采初始,平均日产水量1.98 m3。见气时间47 d,产气曲线呈现双峰型特征,第1个峰出现于210 d左右,产气量接近1000 m3/d,随后逐步降低;排采第306 d,产气量开始明显增加,第2个峰出现于315 d左右,产气量接近2500 m3/d。推测第2个峰源于下部产气系统的解吸,在排采300 d时生产出现波动,体现为井底流压、套压与产气量明显增高,昭示新的且强烈的解吸气源显现(图2 d)。

    具体而言,不同产气系统合采过程中由于生产接替,形成时间序列上具有明显差异的分段式气水产出特征,且接续过程伴随生产曲线的强烈波动,波动程度随产气系统流体能量与资源基础差异增强而增大。Z−4井合采2个产气系统,上部产气系统(6、7号煤层,埋深283~290 m)产水量高、产气量低;下部系统(17、20、23、27、30号煤层,埋深350~531 m)产水量低、产气量高。排采至300 d时,6、7号煤已处于暴露状态,此时产气量的剧增与压力的波动,极大概率源于下部系统的强烈解吸。

    井底流压在第300天开始出现显著的回升现象,代表下部系统开始产气,且直接进入敏感解吸阶段[8, 19],日产气量与套压急剧增高并达到峰值,产水量明显降低。本井出现产气系统接替现象,可归因于产层主体位于300 m以深且以下部产气系统为主,水文地质条件较为封闭,水侵能量有限。产出水地球化学分析也表明排采后期水源以煤层水为主[13]。若上部产气系统水侵能量强烈,可能全程屏蔽下部产气系统,使产能无法释放,比如Z−2井第2生产阶段与Z−1井。

    Z−5井合采16、17、20、23、27号煤层,排采时间314 d,最高日产气量2806 m3,平均日产气量1433 m3,最高日产水量7.01 m3,平均日产水量0.85 m3。排采87 d后见气,见气前日均产水量1.43 m3,见气后日均产水量0.63 m3。从产气曲线来看,存在3个产气峰,第1个峰出现于第160天,产气量为1675 m3/d,第2个峰出现于第185天,产气量2499 m3/d,第3个峰出现于第263天,产气量2806 m3/d。产气效果较好且多煤层接续产气特征明显(图2e)。相比Z−4井,本井产层组合未包含上部6号与7号煤层,表现出单一产气系统特征,流压总体平稳下降,产水产气情况与Z−4井下部产气系统相似。

    Z−S1井排采煤层数量多达14层,几乎涵盖了从上至下的所有可采煤层。最大跨度183 m,累计煤层厚度18.4 m,但产气效果极差。在整个259 d的排采周期内,最高日产气量仅68 m3,平均日产气量5 m3;产水量长期保持在高位水平,最高日产水量44.30 m3,平均日产水量23.08 m3,属于典型的高产水、低产气井。排采至第37天开始第1阶段产气,但产气量极低,持续47 d。排采至第189天开始第2阶段产气,平均日产气量31 m3,持续28 d。排采期间经历2次停抽,停抽后动液面快速恢复至初始水平(甚至高出),说明气井补给条件好,水侵明显,气体无法解吸产出(图2f)。该井油管直径73 mm,采用直径为56 mm的管式泵排采,冲程6 m,冲次根据产水量调整,主要介于1.5~3.0次/min,据此泵的理论排量介于31.90~63.81 m3/d,能够满足本井排水需要。产水量代表了地层的实际供液能力,未受到排采设备与管柱的限制。

    一个新勘探区块的煤层气井排水采气通用的控制原则是“连续、缓慢、长期、稳定”, 液降速率通常分阶段控制,在单相流即起始降液阶段降液速率可相对较快,在初期产气即气液两相流阶段需降低降液速率,在稳定产气阶段动液面基本保持稳定,主要靠控制井口套压控制产气速率。3个阶段的产液、产气特征各不相同。因此,结合研究区排采曲线形态,按照见气时间,产气高峰时间划分排采曲线,可基本对应上述3个阶段。原因在于2个方面:一是黔西织金区块见气后多直接进入敏感解吸阶段,产气高峰来临较早[8],以Z−2井第1阶段为代表,二是产气高峰来临后动液面基本保持稳定状态。通过对比,可获得以下信息:① Z−2−1井、Z−3井、Z−4井和Z−5井的3个阶段平均日产气量逐渐增加,平均日产水量逐渐降低,表现为正常生产规律,峰值日产气量均超过1000 m3/d;② Z−1井和Z−2−2井初期产气阶段日产水量高,而见气前日产水量低,不符合正常规律,峰值日产气量低于800 m3/d,平均日产气量低于500 m3/d;③ Z−1井和Z−2−2井均打开了上部产气系统,且产层组合顶界埋深浅于250 m,易受表层水干扰(表2)。Z−S1井由于强烈水侵,无明确生产规律,不再作阶段划分。

    表  2  煤层气井不同排采阶段产水产气特征
    Table  2.  Characteristics of water and gas production in different production stages of CBM wells
    井号 阶段划分 峰值日产气量/ (m3·d−1 平均日产气量/(m3·d−1 峰值日产水量/(m3·d−1 平均日产水量/(m3·d−1
    Z−1 见气前 0 0 5.60 2.83
    初期产气 427.51 305.88 5.88 4.62
    稳定产气 424.43 304.69 4.01 3.44
    Z−2−1 见气前 0 0 4.08 1.86
    初期产气 2771.84 1974.02 2.80 0.66
    稳定产气 2658.10 1548.40 0.07 0.06
    Z−2−2 见气前 0 0 3.34 1.12
    初期产气 769.92 320.54 6.11 4.18
    稳定产气 765.83 422.12 4.02 3.48
    Z−3 见气前 0 0 2.03 1.70
    初期产气 1137.26 412.79 1.56 0.68
    稳定产气 1076.36 978.71 0.32 0.14
    Z−4 见气前 0 0 21.22 5.69
    初期产气 2465.41 429.59 3.47 1.52
    稳定产气 2434.18 2001.89 1.42 0.76
    Z−5 见气前 0 0 7.01 1.43
    初期产气 2805.94 1356.00 1.90 0.78
    稳定产气 2798.43 1696.79 0.16 0.13
    下载: 导出CSV 
    | 显示表格

    图3为累计产水量与井底流压的关系曲线。

    图  3  累计产水量与井底流压的关系曲线
    Figure  3.  Relationship between cumulative water production and bottom-hole flow pressure

    Z−2井2个生产阶段的累计产水量与井底流压关系曲线呈现了不同特点,第1阶段呈下凹型(图3a),第2阶段呈上凹型(图3b)。前者反映排水能够持续诱导井壁储层压降,属于煤储层弹性水产出;后者反映排水后期不能有效诱导井壁储层压降(排采480 d之后),随累计产水量增加井底流压基本保持稳定,说明发生了水侵。具体而言,第1次降低冲次之后,补给量与产水量基本达到平衡,日产水量维持在3.5 m3/d左右,动液面稳定在350 m左右。因此,可根据累计产水量与井底流压关系曲线形态判断排水降压的有效性。

    对累计产水量与井底流压关系曲线进行分段线性拟合。由于井底流压随累计产水量增加总体呈降低趋势,在讨论拟合直线的斜率时忽略负号,针对斜率的绝对值进行分析。Z−2井第1阶段的第1直线段斜率为0.0192,第2直线段斜率为0.1018,代表每产100 m3水,井底流压分别降低1.92 MPa和10.18 MPa。见气后的排水降压效率显著高于见气前,曲线整体呈现斜率逐渐增大的下凹型,排水有效性逐步增强(图3a)。

    Z−2井第2阶段的3个直线段斜率依次为0.00430.00010.0002,代表每产100 m3水,井底流压分别降低0.43、0.01、0.02 MPa,曲线整体呈现斜率逐渐减小的上凹型,水侵导致排水降压效果逐步减弱(图3b)。

    Z−1井的关系曲线总体呈现上凹型,说明发生水侵。第1直线段每排100 m3水,井底流压下降0.42 MPa,在此区间见气前后降压效率没有发生明显变化。在降低冲次后,曲线发生转折,井底流压随累计产水量增加不再降低而趋于稳定,第2直线段呈水平趋势。经历停抽之后,快排使井底流压迅速下降,随后维持稳定。第2、3直线段排水未能形成有效压降,代表外源水补给量与产水量达到平衡,气井长期处于低产低效状态(图3c)。

    Z−3井的关系曲线具有下凹型特点,显示了较好的降压效果。具体而言,在见套前,累计产水量与井底流压线性关系显著,每产100 m3水,井底流压降低2.75 MPa。见气后,曲线发生弯曲下凹,第2直线段斜率增大,每产100 m3水井底流压降低6.45 MPa。见气后斜率的增大,源于气、水同产阶段气体产出的附加降压效应,以及气体产出后对产水的阻碍效应(图3 d)。

    Z−4井的累计产水量与井底流压关系曲线形态不够明朗。根据见气(套)点与产气高峰点,可识别3个典型的直线段进行拟合,拟合直线斜率逐渐增高,每产100 m3水井底流压分别降低0.43、0.92、1.83 MPa。产气高峰出现时导致的井底流压回升现象,与产气量大幅提高后的憋套压有关。在此阶段套压从0.6 MPa升至1.5 MPa,涨幅0.9 MPa,井底流压从3.0 MPa升至3.8 MPa,涨幅0.8 MPa,两者基本一致。前文述及,本井具有双产气系统特征,下部系统接替产气后,迎来产气高峰,并在井底流压下降阶段,直线斜率进一步增加。据此说明曲线总体呈现下凹型,代表较好的降压效果(图3e)。

    Z−5井的累计产水量与井底流压关系曲线具有下凹型特点。第1直线段位于见气之前,每产100 m3水井底流压降低1.12 MPa,第2拟合直线段发生于见气之后,产气高峰来临之前,每产100 m3水井底流压降低0.80 MPa;第3拟合直线段对应高产阶段,斜率显著增大,每产100 m3水井底流压降低4.07 MPa,该阶段平均日产气量1885 m3,峰值日产气量达2806 m3图3f)。

    Z−S1井的排采存在2次停抽、2次降冲次以及2次产气(极低效产气),导致累计产水量与井底流压关系曲线形态较为复杂。但在正常排采阶段,仍然存在3个较为明显的直线段。拟合结果显示,3个直线段每产100 m3水分别降压0.10、0.06、0.12 MPa,降压效率极低,且停抽后井底流压能够快速恢复,反映水侵强烈(图3 g)。

    Z−2井第1阶段拟合结果显示,斜率倒数指示井底流压每降低1 MPa的产气量为684932 m3图4a)。Z−2井第2阶段线性拟合结果为925926 m3图4b)。可见,单位井底流压降幅下,Z−2井第2阶段的产气量大于第1阶段。

    图  4  累计产气量与井底流压的关系曲线
    Figure  4.  Relationship between cumulative gas production and bottom-hole flow pressure

    Z−1井、Z−3井典型直线段拟合结果显示,单位压降产气量分别为1259450 m3649351 m3图4c图4 d)。Z−4井2个产气系统分别形成了2个典型直线段,单位井底流压降幅产气量分别为124378 m3246305 m3,后者大致为前者的2倍(图4e)。Z−5井存在2个典型直线段,指示单位井底流压降幅产气量依次为232558 m31257862 m3,呈现增高趋势,反映了较好的供气潜力(图4f)。Z−S1井产气量极低且产气天数少,不满足选段拟合标准,不再分析。

    综合而言,Z−1井与Z−5井的降压产气效率最高,但这2口井的实际产气效果与水侵程度差别巨大,说明降压诱导产气的效率不是关键制约因素。

    除产气初期的过渡阶段外,Z−2井累计产水量与累计产气量关系曲线总体呈现直线型特点。2个生产阶段差异显著,第1阶段每产1 m3水能产14207 m3煤层气,第2阶段仅能产122 m3煤层气,为前者的0.86%,反映了水侵导致的产气受限(图5a)。

    图  5  累计产水量与累计产气量的关系曲线
    Figure  5.  Relationship between cumulative water production and cumulative gas production

    Z−1井、Z−3井、Z−5井的累计产水量与累计产气量线性拟合显示,每产1 m3水产气量分别为91、735212043 m3,前者水侵干扰明显,后两者产气效率较高,无明显水侵特征(图5b图5c图5e)。

    Z−4井的累计产水量与累计产气量曲线存在两个明显的直线段,分别对应2个产气系统,验证了前文关于双产气系统的论述。每产1 m3水,第1直线段产气445 m3,第2直线段产气4160 m3,后者为前者的9.35倍,下部产气系统无明显水侵现象,产气效率高(图5 d)。

    Z−S1井的两个直线段显示产气效率极低,每产1 m3水分别产气0.12 m3和1.09 m3,气井变水井,水侵严重(图5f)。

    排水降压有效性:单位产水量的井底流压降幅以Z−2−1,Z−3,Z−4,Z−5井较高,且这些井不同阶段降压效率呈现逐渐增高的趋势(下凹型),与气体大量产出有关,可判识为非水侵井,高产潜力大。而Z−1井、Z−2−2井、Z−S1井的单位产水量井底流压降幅很低,均低于0.005 MPa/m3,且呈现逐渐降低的趋势(上凹型),为水侵井(图6a)。据此,可将排采初期(见气前)单位产水量井底流压降幅0.005 MPa/m3,视为水侵的临界值,低于此值,极大可能发生了水侵,若后续此值进一步降低,则可确定发生水侵。无水侵井的判识标志为排采初期单位产水量井底流压降幅大于0.005 MPa/m3,且后续有递增趋势。另外,产出水中的Cl、Na+离子浓度可作为分析滞留程度与排采水源的重要依据,3口水侵井的Cl+Na+离子质量浓度平均值介于265.66~483.86 mg/L,而非水侵井上述平均值介于906.94~1366.89 mg/L,研究区地表水样的Cl+Na+离子质量浓度介于6.95~12.07 mg/L,产出水地球化学封闭性的差异佐证了上述分类,也反映了水侵的客观存在[16]

    图  6  排水—降压—产气有效性分析
    Figure  6.  Effectiveness analysis of drainage, depressurization and gas production

    降压产气有效性:单位井底流压降幅的产气量总体处于较高水平,除Z−S1井外,其余井均高于10万m3/MPa,尤其Z−1井、Z−5井超过100万m3/MPa,Z−2−2井接近100万m3/MPa(图6b)。前文指出,Z−1井与Z−2−2井均为水侵井,产气效率低,但单位井底流压降幅的产气量却最高,说明限制产气的关键在于排水如何诱导降压。进一步而言,水侵导致的低效降压是制约织金地区煤层气高效开发的关键,一旦能够有效降压,实现高产的潜力就会很大,这与黔西地区煤层较高的含气量、较短的吸附时间密切相关[8,23]

    单位产水量的产气量井间差异与单位产水量的井底流压降幅井间差异高度一致(图6c),佐证了上述观点。如何限制水侵,将是该区低产井提产改造的基本出发点。织金区块煤层气勘探开发目标层位经历了从上部煤层群向下部煤层群的转换,水侵无疑是重要原因之一[3,21,35-37]

    基于排采有效性分析曲线提取气藏属性参数,结果显示(表3),3口水侵井均具有启抽压力低(<3 MPa),弹性产水指数高(>200 m3/MPa)的特点。考虑到水侵井均打开了上部产气系统,说明上部系统富水性强,且水动力系统开放,压裂沟通含水层导致储层卸压。非水侵井启抽压力高(>4 MPa),弹性产水指数低(<100 m3/MPa)。Z−4井的弹性产水指数为233 m3/MPa,但后续发生产气系统的接续,有效抑制了水侵。非水侵井产层主体隶属中下部产气系统,可见中下部系统富水性弱,水文地质开启程度低,开发过程不易发生水侵,能够保证降压效果。启抽压力实质反映排采初始静液面深度,水侵井初始静液面深度大,介于69~141 m,反映了漏失现象;非水侵井初始静液面深度小,介于0~19 m。

    表  3  气藏属性参数提取与产水量劈分结果
    Table  3.  Results of gas reservoir property parameters extraction and water production sources identification
    气井 气藏参数 产水量劈分
    启抽压力/MPa 弹性产水指数/(m3·MPa−1) 井控动态储量/106m3 弹性量/m3 水侵量/m3 总量/m3
    Z−1 2.65 238 0.19 500 1570 2070
    Z−2−1 4.47 52 0.84 147 0 147
    Z−2−2 2.87 233 0.70 709 1216 1925
    Z−3 7.22 36 0.44 203 0 203
    Z−4 5.73 233 1.06 659 38 697
    Z−5 4.13 89 0.57 265 0 265
    Z−S1 2.95 1000 5978
    下载: 导出CSV 
    | 显示表格

    井控动态储量而言,以Z−4井最高,说明产气系统接续后,其具有最高的开发潜力。水侵井的井控动态储量并未明显低于非水侵井(表3),说明资源条件不是关键制约因素,与黔西地区煤层气资源量大、资源丰度高的特点相符[32]

    产水量劈分结果显示,水侵井的水侵产水量大于弹性产水量,总产水量达到或超过2000 m3;非水侵井主要产出弹性水,总产水量低于700 m3表3)。两类井的产水能量驱动类型差异显著,水侵能量越高,降压产气效果越差。Z−S1井为强烈水侵井,生产波动极大,曲线极不稳定,不再进行产水量劈分。Z−S1井发生如此强烈水侵的原因需要依靠地震资料作进一步阐释,目前尚不能断定是否受到断层的影响,但该井排采煤层数量最多(14层),产层组顶部埋深270 m,平均埋深360 m,无疑增加了沟通含水层的机率与程度,构成高产水的潜在原因。

    1)构建了基于累计产水量与井底流压、累计产气量与井底流压、累计产水量与累计产气量关系曲线的煤层气井排采有效性分析方法,根据曲线形态与典型直线段拟合结果,能够识别水侵现象,获取气藏属性参数,劈分弹性产水量与水侵产水量,剥离排水降压有效性与降压产气有效性。该方法对揭示该区域煤层气井产能关键制约因素,促进开发方案优化调整具有指导意义。

    2)在黔西织金区块的应用表明,当累计产水量与井底流压的关系曲线呈现上凹型,指示发生水侵,对应低产井;下凹型则代表未发生水侵,气体大量产出导致曲线下凹,对应高产井。限制或避免水侵将是织金区块老井提产改造与新井开发设计的重要思路。

    3)水侵井具有启抽压力低,弹性产水指数高的特点,非水侵井与之相反。水侵井的水侵产水量大于弹性产水量,总产水量高;非水侵井主要产出弹性水,总产水量低。井控动态储量在两类井之间无规律性差异,资源条件不是主要制约因素。启抽压力低或排采初始静液面深度大、见气前单位产水量井底流压降幅低,可作为水侵井的识别标志。

    4)提出产气系统概念,不同产气系统合采兼容性差。黔西织金区块上部产气系统水文地质开启程度高,容易发生水侵,产气效率低,以Z−1、Z−2−2、Z−S1井为代表;中下部产气系统水文地质条件封闭,无明显水侵现象,产气效率高,以Z−2−1,Z−3,Z−5井为代表;Z−4井表现出产气系统接替特征,下部系统接替后具有高产潜力。研发针对表层水的降扰技术与不同产气系统的产能有效接替技术将是下一步攻关方向。

计量
  • 文章访问数:  56
  • HTML全文浏览量:  7
  • PDF下载量:  132
  • 被引次数: 0
出版历程
  • 网络出版日期:  2023-04-02
  • 发布日期:  2022-11-24

目录

/

返回文章
返回