石油学报 ›› 2022, Vol. 43 ›› Issue (10): 1474-1486.DOI: 10.7623/syxb202210010

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

基于扩展有限元的离散缝网渗流数值模拟方法

邓英豪, 夏阳, 金衍   

  1. 中国石油大学(北京)石油工程学院 油气资源与探测国家重点实验室 北京 102249
  • 收稿日期:2021-03-21 修回日期:2022-02-26 出版日期:2022-10-25 发布日期:2022-11-05
  • 通讯作者: 金衍,男,1972年8月生,2001年获石油大学(北京)博士学位,现为中国石油大学(北京)教授、博士生导师,主要从事岩石力学等石油工程方向的研究。Email:jiny@cup.edu.cn
  • 作者简介:邓英豪,男,1996年12月生,2020年获中国石油大学(北京)石油工程专业学士学位,现为中国石油大学(北京)石油工程学院石油与天然气工程专业博士研究生,主要从事岩石力学数值模拟等方向的研究。Email:dengyinghao2020@163.com
  • 基金资助:
    国家自然科学基金企业创新发展联合基金项目(No.U19B6003-05)和国家自然科学基金青年科学基金项目(No.51904318)资助。

Numerical simulation method of discrete fracture network flow based on the extended finite element method

Deng Yinghao, Xia Yang, Jin Yan   

  1. College of Petroleum Engineering, China University of Petroleum;State Key Laboratory of Petroleum Resources and Prospecting, Beijing 102249, China
  • Received:2021-03-21 Revised:2022-02-26 Online:2022-10-25 Published:2022-11-05

摘要: 裂缝性储层具有强非均质性和多尺度特征,传统方法在离散裂缝群数值模拟中对网格设置有较高的要求,通常难以平衡计算效率和精度,在高效且精细模拟多尺度无序交叉裂缝渗流方面的局限性日益凸显。基于离散裂缝模型显示表征每一条裂缝,针对不同尺度裂缝构造基于渐近解析解的加强形函数,建立了一种裂缝和网格解耦的多尺度扩展有限单元法用于模拟无序交叉裂缝群的渗流特征。该方法利用渐近解析解捕捉网格内裂缝局部渗流特征,避免使用非结构网格在裂缝壁面和尖端的局部加密,因此能有效降低计算量,提高计算精度。建立了单裂缝模型和大规模裂隙群模型用于模拟,新方法相较于传统有限单元法将裂缝壁面法向流速的计算精度提高80 %以上;在采用粗网格模拟大规模裂隙群时,在计算量(网格数和结点数)降低约90 %的情况下,单井产量误差保持在仅约5 %,验证了新方法的准确性和有效性。新方法为裂缝性油气藏流体流动模拟以及后续解决两相流、流固耦合、裂缝扩展等一系列复杂问题提供了一种新的思路和方法。

关键词: 多尺度, 扩展有限单元法, 离散裂缝, 加强形函数, 渗流

Abstract: Fractured reservoirs have strong inhomogeneity and multi-scale characteristics, and the traditional methods have high requirements on the grid setting during the numerical simulation of discrete fracture groups, which are usually difficult to balance the computational efficiency and accuracy; moreover, the limitations in efficient and elaborate simulation of multi-scale disordered cross-fracture seepage are becoming more and more prominent. Based on the discrete fracture model, each crack is explicitly described. Enrichment functions based on the asymptotic analytical solution are constructed for different scales of factures. And the multi-scale extended finite element method is constructed to simulate the flow characteristics of disordered intersected fracture groups with grid-crack decoupling. The method uses the asymptotic analytical solution to capture the local seepage characteristics of cracks in the grid, avoiding the local refinement at the crack walls and tips by use of unstructured grids, thus effectively reducing the computational amount and improving the computational accuracy. A single fracture model and a large-scale fracture group model are established for simulation. Using this method, the calculation accuracy of normal flux along the crack wall can be improved by more than 80 % when compared with the traditional finite element method; when using coarse grids to simulated a large-scale fracture cluster, the calculation complexity (number of grids and nodes) is reduced by about 90 %, and the error of single well production is maintained at only 5 %, which verifies the accuracy and effectiveness of this method. The method proposed in this paper provides a new idea and method for the simulation of fluid flow in fractured reservoirs and the subsequent solution of a series of complex problems such as two-phase flow, fluid-solid coupling, and fracture propagation.

Key words: multi-scale, extended finite element method, discrete fractures, enrichment functions, flow

中图分类号: