油气藏评价与开发 ›› 2025, Vol. 15 ›› Issue (6): 1025-1033.doi: 10.13809/j.cnki.cn32-1825/te.2025.06.008

• 油气开发 • 上一篇    下一篇

碎裂煤中煤粉运移对孔渗影响的数值模拟研究

石慧1(), 谢天成1, 刘子亮1, 蒋志坤1, 魏迎春1,2()   

  1. 1.中国矿业大学(北京)地球科学与测绘工程学院,北京 100083
    2.中国矿业大学(北京)煤炭精细勘探与智能开发全国重点实验室,北京 100083
  • 收稿日期:2024-09-19 发布日期:2025-10-24 出版日期:2025-12-26
  • 通讯作者: 魏迎春(1977—),女,博士,教授,主要从事煤与煤系矿产资源地质勘探与开发方面的教学和科研工作。地址:北京市海淀区学院路丁11号,邮政编码:100083。E-mail: wyc@cumtb.edu.cn
  • 作者简介:石慧(2001—),女,在读硕士研究生,从事煤层气开发地质研究。地址:北京市海淀区学院路丁11号,邮政编码:100083。E-mail: shih_forward@163.com
  • 基金资助:
    国家自然科学基金项目“不同类型构造煤储层中煤粉运移规律及其储层意义”(41972174);国家自然科学基金项目“惰质组与镜质组大分子结构演化差异性及其动力学机制”(42042197)

Numerical simulation study on influence of coal fines migration on porosity and permeability in cataclastic coal

SHI Hui1(), XIE Tiancheng1, LIU Ziliang1, JIANG Zhikun1, WEI Yingchun1,2()   

  1. 1. College of Geoscience and Surveying Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China
    2. State Key Laboratory for Fine Exploration and Intelligent Development of Coal Resources, China University of Mining and Technology (Beijing), Beijing 100083, China
  • Received:2024-09-19 Online:2025-10-24 Published:2025-12-26

摘要:

在煤层气排采过程中,储层中煤粉会发生运移,可能造成孔隙和喉道堵塞,导致渗透率下降,从而影响煤层气的产量。为了探究煤粉运移对碎裂煤储层孔渗的影响,建立了煤粉在孔喉通道内的运移—沉降数值模型,基于python(高级程序设计语言)编写了煤粉运移的数值模拟程序,模拟了煤粉在碎裂煤基质孔隙网络中的运移,探讨了煤粉运移过程中的碎裂煤储层孔渗变化规律以及煤粉运移对孔渗的影响。数值模拟结果表明:压差和煤粉粒度是影响储层孔渗的2个关键因素。模拟初始阶段,煤粉运移引起储层渗透率迅速降低,不同粒径煤粉颗粒运移、沉降和排出进一步受到压差的影响,煤粉粒度越大,在低压差、低流速条件下越难以启动运移,而在压差较大流速较高的情况下,则会发生运移。对比小粒径煤粉,大粒径煤粉颗粒更加容易堵塞有效孔隙,导致渗透率迅速下降。压差增加会使煤粉沉积位置向出口端转移,导致煤粉沉降范围变大。在煤粉粒径小于喉道半径的条件下,驱替压力存在一个阈值,当小于这个阈值时,渗透率的下降速度随着驱替压力的增大而增大;当大于这个阈值时,渗透率的下降速度会随着驱替压力的增大而降低。将数值模拟实验数据与关于煤粉运移的物理模拟实验结果相结合表明:模型的数值模拟结果与物理模拟实验中的煤粉运移过程中碎裂煤的渗透率变化和煤粉产出沉降情况一致。

关键词: 煤层气, 碎裂煤, 煤粉运移, 孔渗, 煤粉粒径, 数值模拟

Abstract:

During coalbed methane (CBM) production, coal fines within the reservoir can migrate, potentially blocking pore throats and resulting in a significant reduction in reservoir permeability. This process adversely affects the final CBM yield. To investigate the influence of coal fines migration on the porosity and permeability of cataclastic coal reservoirs, this study focuses on the processes of fines initiation, migration, and deposition within reservoir channels. The pore size distribution characteristics of cataclastic coal were analyzed using low-field nuclear magnetic resonance and low-temperature liquid nitrogen adsorption experiments. Subsequently, a three-dimensional pore network model was constructed, and a numerical model for coal fines migration and deposition in pore throat channels was developed. By integrating with an existing mechanical model for coal fines initiation and a probabilistic model for particle deposition and throat blockage, the Monte Carlo method was used to simulate the migration and blockage of coal fines within reservoir pores. A numerical simulation program written in Python was developed to simulate coal fines migration within the pore network of the cataclastic coal matrix. The variations in porosity and permeability of the reservoir during coal fines migration, as well as the influence of this migration, were discussed. The analysis revealed the internal mechanisms by which pressure difference and coal fines particle size influenced coal fines output and model permeability. Both factors were significant, and their interactions were complicated. Specifically, the particle size of coal fines directly affected their migration, deposition, and output characteristics under different hydrodynamic conditions. Under low pressure difference and low flow velocity, large coal fines particles were difficult to mobilize and initiate migration. In contrast, under high pressure difference and high flow velocity, these particles became mobile but were more likely to block effective pores, resulting in a sharp decline in permeability. In addition, an increase in pressure difference had dual effects. It promoted coal fines output but also accelerated permeability decline rate. As the displacement pressure difference increased, the deposition location of coal fines shifted toward the outlet end, accompanied by a higher proportion of small throats. When the coal fines particle size was constant and smaller than the throat radius, a displacement pressure threshold was observed. On either side of this threshold, the relationship between the permeability decline rate and the displacement pressure showed distinct trends. During the drainage and depressurization stages of actual CBM production, the output characteristics and particle size distribution of coal fines served as important indicators for evaluating production efficiency and reservoir permeability changes. As drainage intensity gradually increased, the output intensity of coal fines experienced an initial slow growth followed by a rapid decline. Simultaneously, the particle size distribution of the produced coal fines reached its widest range, encompassing sizes from small to large, particularly during the initial drainage stage. When the drainage intensity was low, only small coal fines particles were mobilized and produced by fluid flow. To further investigate this, numerical simulations were conducted to replicate low drainage intensity conditions by setting a low pressure difference. The simulation results indicated that under low flow rates, small coal fines could indeed migrate and be produced, accompanied by a relatively gentle decline in reservoir permeability. These results were consistent with field observations under low drainage intensity, confirming the accuracy and reliability of the numerical simulation in predicting coal fines migration behavior. Furthermore, the numerical simulation results were compared with those from physical simulation experiments on coal fines migration. The model simulation results showed that as drainage continued, large coal fines particles settled preferentially within the pores. The deposition of these large particles formed channel barriers, blocking pore throats and significantly reducing permeability. Simultaneously, the deposition probability of small coal fines also increased rapidly. This made subsequent migration and production of coal fines increasingly difficult, resulting in a rapid decrease in coal fines production over time. In addition, physical simulation experiments of coal fines migration in cataclastic coal reservoirs provided valuable reference data. The experimental results showed that the permeability decrease of coal samples primarily occurred during the early stage of water flooding, and a higher fracture density corresponded to a higher average permeability. Pores and fissures with diameters greater than 1 000 nm served as the main channels for coal fines migration and blockage. The migrating coal fines further reduced permeability by plugging connected pores larger than 10 000 nm. These findings were consistent with the numerical simulation results, further confirming the alignment between permeability evolution and coal fines production during migration in cataclastic coal reservoirs. In simulations, high drainage intensity was simulated by setting a high pressure difference. When large coal fines were introduced, their deposition locations shifted compared with those under low pressure differences. Also, fines output decreased, and deposition became concentrated near the outlet end. When small coal fines were introduced, fines output increased significantly, and the permeability reduction was less severe than that caused by large coal fines. Overall, the model simulation results were consistent with coal fines output observed in actual production under increased drainage intensity. The numerical simulation results indicated that larger coal fines particles led to more concentrated deposition and blockage near the inlet end, with a smaller overall deposition range. Under low pressure difference, coal fines deposition was mainly concentrated near the inlet end. As the pressure difference increased, the deposition locations of coal fines shifted closer to the outlet end, and the deposition range expanded. In visual physical simulations of coal fines migration within fractures, the deposition area gradually decreased from the inlet toward both ends along and perpendicular to the main migration direction. In summary, the numerical simulation results were consistent with the permeability changes and coal fines deposition patterns observed during coal fines migration in physical simulation experiments.

Key words: coalbed methane, cataclastic coal, coal fines migration, porosity and permeability, coal fines particle size, numerical simulation

中图分类号: 

  • TE357