特种油气藏 ›› 2026, Vol. 33 ›› Issue (4): 122-131.DOI: 10.3969/j.issn.1006-6535.2026.04.016

• 钻采工程 • 上一篇    下一篇

基于精细地质模型的水平井压裂区采气过程渗流-应力耦合有限元模拟

刘厚彬1, 杜爽1, 张立双2, 王帅1, 崔帅3   

  1. 1.西南石油大学,四川 成都 610500;
    2.中国石油吉林油田分公司地球物理勘探研究院,吉林 松原 138000;
    3.中国石油西南油气田分公司页岩气研究院,四川 成都 610051
  • 收稿日期:2025-12-30 修回日期:2026-05-28 出版日期:2026-08-25 发布日期:2026-10-15
  • 作者简介:刘厚彬(1980—),男,教授,2003年毕业于西南石油学院勘查技术与工程专业,2009年毕业于西南石油大学油气井工程专业,获博士学位,现主要从事岩石力学与井壁稳定等方面的教学与研究工作。
  • 基金资助:
    国家自然科学基金“深层页岩表面水化诱发细观结构损伤与宏观力学性能劣化机制研究”(52474008);四川省自然科学基金“泥岩水化学微观作用机理与物理化学耦合井壁稳定研究”(24NSFSC2397)

Flow in porous media-stress coupled finite-element simulation of gas production in the fractured zone of horizontal wells based on a refined geological model

Liu Houbin1, Du Shuang1, Zhang Lishuang2, Wang Shuai1, Cui Shuai3   

  1. 1. Southwest Petroleum University,Chengdu,Sichuan 610500,China;
    2. Research Institute of Geophysical Exploration,PetroChina Jilin Oilfield Company,Songyuan,Jilin 138000,China;
    3. Shale Gas Research Institute,PetroChina Southwest Oil & Gasfield Company,Chengdu,Sichuan 610051,China
  • Received:2025-12-30 Revised:2026-05-28 Online:2026-08-25 Published:2026-10-15

摘要: 为研究页岩气水平井压裂区采气过程的渗流-应力耦合机制与孔隙压力、地应力时空演化规律,建立了一种基于精细地质模型的有限元耦合模拟方法。该方法通过裂缝面嵌入地质网格及三维属性插值函数,将多尺度裂缝精确映射至有限元网格。在此基础上,构建了考虑非达西渗流及孔隙度、渗透率动态演化的流-固耦合模型,根据模型建立了生产井模拟压力与实际压力误差梯度下降的生产历史拟合方法。以四川泸州区块页岩气X井平台为例,在初始孔隙压力和渗透率动态标定及生产历史拟合的基础上,分析了采气驱动及裂缝系统影响下孔隙压力与应力的耦合响应。结果表明:历史拟合平均相对误差为4.6%,模型可准确反映实际生产动态;在82 d的短期采气过程中,井周孔隙压力下降约1.5 MPa,有效应力相应增加约0.9 MPa;井周压降漏斗与应力集中区同步向外扩散,影响范围主要限制在井周50 m以内;孔隙压力下降及岩石体积应变导致孔隙度、渗透率缓慢降低,降幅小于0.5%,表明短期开采对储层物性的影响有限;裂缝系统作为优势渗流通道加速压力传播,同时在力学响应中形成应力集中,使流-固耦合呈现显著的各向异性和非均质性特征。研究成果可为页岩气开发的产能评价、井间干扰分析、套管完整性评估等工程问题提供定量研究手段。

关键词: 页岩气, 水平井, 压裂区, 渗流-应力耦合, 裂缝, 孔隙压力, 有限元方法

Abstract: To investigate the flow in porous media-stress coupling mechanism during gas production from the fractured zone of shale-gas horizontal wells and the spatiotemporal evolution of pore pressure and in-situ stress,a finite-element coupled simulation method based on a refined geological model was established.By embedding fracture surfaces into geological grids and using a three-dimensional property interpolation function,this method accurately maps multiscale fractures onto finite-element grids.On this basis,a fluid-solid coupling model considering non-Darcy flow and the dynamic evolution of porosity and permeability was constructed.According to the model,a production-history matching method was established using gradient descent of the error between simulated and actual pressures of production wells.With the Well X platform in the Luzhou Block,Sichuan,as an example,the coupled responses of pore pressure and stress under gas-production drive and the influence of the fracture system were analyzed on the basis of dynamic calibration of initial pore pressure and permeability and production-history matching.The results show that the average relative error of history matching was 4.6%,indicating that the model accurately reflects actual production performance.During short-term gas production over 82 d,the near-wellbore pore pressure decreased by approximately 1.5 MPa,and the effective stress increased correspondingly by approximately 0.9 MPa.The near-wellbore pressure-depletion funnel and stress-concentration zone expanded outward synchronously,with the affected range mainly confined within 50 m of the well.The decrease in pore pressure and the volumetric strain of the rock led to slow reductions in porosity and permeability,with decreases of less than 0.5%,indicating that short-term production had limited effects on reservoir physical properties.As preferential seepage pathways,the fracture system accelerated pressure propagation and simultaneously formed stress concentration in the mechanical response,causing fluid-solid coupling to exhibit pronounced anisotropy and heterogeneity.The results provide a quantitative research tool for engineering issues such as productivity evaluation,interwell interference analysis,and casing-integrity assessment in shale-gas development.

Key words: shale gas, horizontal well, fractured zone, flow in porous media-stress coupling, fracture, pore pressure, finite-element method

中图分类号: