石油勘探与开发 >

2020 , Vol. 47 >Issue 1: 163 - 174

DOI: https://doi.org/10.11698/PED.2020.01.16

石油工程

水平井分段压裂平面三维多裂缝扩展模型求解算法

  • 陈铭 ,
  • 张士诚 ,
  • 胥云 ,
  • 马新仿 ,
  • 邹雨时
展开
  • 1. 中国石油大学（北京）,北京 102249;
    2. Texas A & M University,College Station 77840,USA;
    3. 中国石油勘探开发研究院,北京 100083;
    4. 中国石油油气藏改造重点实验室,河北廊坊 065007
陈铭（1990-）,男,山东泰安人,中国石油大学（北京）在读博士研究生,主要从事水力压裂理论与数值模拟研究工作。地址：北京市昌平区府学路18号,中国石油大学（北京）289信箱,邮政编码：102249。E-mail: xmcm0122@126.com

收稿日期: 2019-04-28

  网络出版日期: 2020-01-17

基金资助

国家科技重大专项“储层改造关键技术及装备”（2016ZX05023）; 国家重点基础研究发展计划（973）项目“陆相致密油高效开发基础研究”（2015CB250903）

收起

A numerical method for simulating planar 3D multi-fracture propagation in multi-stage fracturing of horizontal wells

  • CHEN Ming ,
  • ZHANG Shicheng ,
  • XU Yun ,
  • MA Xinfang ,
  • ZOU Yushi
Expand
  • 1. China University of Petroleum, Beijing 102249, China;
    2. Texas A&M University, College Station 77840, USA;
    3. PetroChina Research Institute of Petroleum Exploration & Development, Beijing 100083, China;
    4. The Key Laboratory of Research Stimulation, PetroChina, Langfang 065007, China

Received date: 2019-04-28

  Online published: 2020-01-17

Fold

摘要

针对多层油气藏水平井分段多簇压裂设计问题,提出平面三维多裂缝扩展模型高效解法。采用三维边界积分方程计算固体变形,考虑井筒-射孔-裂缝耦合流动及缝内流体滤失。利用显式积分算法求解流固耦合方程,根据尖端统一解析解和最短路径算法计算裂缝扩展边界。方法的准确性通过与径向裂缝解析解、隐式水平集算法和有机玻璃压裂实验对比得到全面验证。与隐式水平集算法相比,新算法计算速度大幅提高。以浙江油田页岩气水平井为例,重点分析了平面应力分布、射孔数分布等对多簇裂缝扩展与进液的影响。研究表明,减小单簇射孔数可以平衡簇间应力非均质分布的影响;调整各簇射孔数量可以实现均衡进液,各簇射孔数差别应控制在1~2孔;增加高应力簇的射孔数有利于均匀进液,但进液均匀并不等于裂缝形态一致,裂缝形态受应力干扰和层间应力剖面共同控制。图17参41

关键词： 水平井; 分段压裂; 多裂缝扩展; 三维边界元; 平面应力非均质; 射孔优化

本文引用格式

陈铭 , 张士诚 , 胥云 , 马新仿 , 邹雨时 . 水平井分段压裂平面三维多裂缝扩展模型求解算法[J]. 石油勘探与开发, 2020 , 47(1) : 163 -174 . DOI: 10.11698/PED.2020.01.16

Abstract

To resolve the issue of design for multi-stage and multi-cluster fracturing in multi-zone reservoirs, a new efficient algorithm for the planar 3D multi-fracture propagation model was proposed. The model considers fluid flow in the wellbore-perforation-fracture system and fluid leak-off into the rock matrix, and uses a 3D boundary integral equation to describe the solid deformation. The solid-fluid coupling equation is solved by an explicit integration algorithm, and the fracture front is determined by the uniform tip asymptotic solutions and shortest path algorithm. The accuracy of the algorithm is verified by the analytical solution of radial fracture, results of implicit level set algorithm and results of organic glass fracturing experiment. Compared with the implicit level set algorithm (ILSA), the new algorithm is much higher in computation speed. The numerical case study is conducted based on a horizontal well in shale gas formation of Zhejiang oilfield. The impact of stress heterogeneity among multiple clusters and perforation number distribution on multi-fracture growth and fluid distribution among multiple fractures are analyzed by numerical simulation. The results show that reducing perforation number in each cluster can counteract the effect of stress contrast among perforation clusters. Adjusting perforation number in each cluster can promote uniform flux among clusters, and the perforation number difference should better be 1-2 among clusters. Increasing perforation number in the cluster with high in situ stress is conducive to uniform fluid partitioning. However, uniform fluid partitioning is not equivalent to uniform fracture geometry. The fracture geometry is controlled by the stress interference and horizontal principal stress profile jointly.

Key words： horizontal well; multi-stage fracturing; multi-fracture growth; 3D boundary element; planar stress heterogeneity; perforation optimization

参考文献

[1] 吴奇, 胥云, 王腾飞, 等. 增产改造理念的重大变革: 体积改造技术概论[J]. 天然气工业, 2011, 31(4): 7-12.
WU Qi, XU Yun, WANG Tengfei, et al.The resolution of reservoir stimulation: An introduction of volume fracturing[J]. Natural Gas Industry, 2011, 31(4): 7-12.
[2] 吴奇, 胥云, 王晓泉, 等. 非常规油气藏体积改造技术: 内涵、优化设计与实现[J]. 石油勘探与开发, 2012, 39(3): 252-258.
WU Qi, XU Yun, WANG Xiaoquan, et al.Volume fracturing technology of unconventional reservoirs: Connotation, optimization design and implementation[J]. Petroleum Exploration and Development, 2012, 39(3): 252-258.
[3] 吴奇, 胥云, 张守良, 等. 非常规油气藏体积改造技术核心理论与优化设计关键[J]. 石油学报, 2014, 35(4): 706-714.
WU Qi, XU Yun, ZHANG Shouliang, et al.The core theories and key optimization designs of volume stimulation technology for unconventional reservoirs[J]. Acta Petrolei Sinica, 2014, 35(4): 706-714.
[4] 胥云, 雷群, 陈铭, 等. 体积改造技术理论研究进展与发展方向[J]. 石油勘探与开发, 2018, 45(5): 874-887.
XU Yun, LEI Qun, CHEN Ming, et al.Progress and development of volume stimulation techniques[J]. Petroleum Exploration and Development, 2018, 45(5): 874-887.
[5] SOMANCHI K, BREWER J, REYNOLDS A.Extreme limited entry design improves distribution efficiency in plug-n-perf completions: Insights from fiber-optic diagnostics[R]. SPE 184834, 2017.
[6] LECAMPION B, BUNGER A, ZHANG X.Numerical methods for hydraulic fracture propagation: A review of recent trends[J]. Journal of Natural Gas Science and Engineering, 2018, 49: 66-83.
[7] ADACHI J, SIEBRITS E, PEIRCE A, et al.Computer simulation of hydraulic fractures[J]. International Journal of Rock Mechanics and Mining Sciences, 2007, 44(5): 739-757.
[8] SETTARI A, CLEARY M P.Three-dimensional simulation of hydraulic fracturing[J]. Journal of Petroleum Technology, 1982, 36(7): 1177-1190.
[9] KRESSE O, COHEN C, WENG X W, et al.Numerical modeling of hydraulic fracturing in naturally fracture formations[R]. ARMA 11-363, 2011.
[10] 赵金洲, 陈曦宇, 李勇明, 等. 水平井分段多簇压裂模拟分析及射孔优化[J]. 石油勘探与开发, 2017, 44(1): 117-124.
ZHAO Jinzhou, CHEN Xiyu, LI Yongming, et al.Numerical simulation of multi-stage fracturing and optimization of perforation in a horizontal well[J]. Petroleum Exploration and Development, 2017, 44(1): 117-124.
[11] DONTSOV E V, PEIRCE A P.An enhanced pseudo-3D model for hydraulic fracturing accounting for viscous height growth, non-local elasticity, and lateral toughness[J]. Engineering Fracture Mechanics, 2015, 142: 116-139.
[12] ADVANI S H, LEE T S, LEE J K.Three-dimensional modeling of hydraulic fractures in layered media: part I: Finite element formulations[J]. Journal of Energy Resources Technology, 1990, 112(1): 1-9.
[13] BARREE R D.A practical numerical simulator for three-dimensional fracture propagation in heterogeneous media[R]. SPE 12273, 1983.
[14] PEIRCE A P, DETOURNAY E.An implicit level set method for modeling hydraulically driven fractures[J]. Computer Methods in Applied Mechanics and Engineering, 2008, 197(33): 2858-2885.
[15] DONTSOV E V, PEIRCE A P.A non-singular integral equation formulation to analyse multiscale behaviour in semi-infinite hydraulic fractures[J]. Journal of Fluid Mechanics, 2015, 781(3): 248-254.
[16] DONTSOV E V, PEIRCE A P.A multiscale implicit level set algorithm (ILSA) to model hydraulic fracture propagation incorporating combined viscous, toughness, and leak-off asymptotics[J]. Computational Methods in Applied Mechanic and Engineering, 2017, 313: 53-84.
[17] CHEN X Y, LI Y M, ZHAO J Z, et al.Numerical investigation for simultaneous growth of hydraulic fractures in multiple horizontal wells[J]. Journal of Natural Gas Science and Engineering, 2017, 51: 44-52.
[18] CARTER B J, DESROCHES J, INGRAFFEA A R, et al.Simulating fully 3D hydraulic fracturing[R]. Ithaca: Cornell University, 2000.
[19] XU G S, WONG S W.Interaction of multiple non-planar hydraulic fractures in horizontal wells[R]. IPTC 17043-MS, 2013.
[20] WANG J, REN L, XIE L Z, et al.Maximum mean principal stress criterion for three-dimensional brittle fracture[J]. International Journal of Solids and Structures, 2016, 102: 142-154.
[21] ZHANG X, JEFFREY R G, THIERCELIN M.Mechanics of fluid-driven fracture growth in naturally fractured reservoirs with simple network geometries[J]. Journal of Geophysical Research: Solid Earth, 2010, 114(B12): 1-16.
[22] WENG X, KRESSE O, COHEN C, et al.Modeling of hydraulic fracture network propagation in a naturally fractured formation[R]. SPE 140253, 2011.
[23] LI S B, ZHANG D X.A fully coupled model for hydraulic-fracture growth during multiwall-fracturing treatments: Enhancing fracture complexity[J]. SPE Journal, 2018, 33(2): 235-250.
[24] ZOU Y S, MA X F, ZHANG S C, et al.Numerical investigation into the influence of bedding plane on hydraulic fracture network propagation in shale formations[J]. Rock Mechanics & Rock Engineering, 2016, 49(9): 3597-3614.
[25] PEIRCE A P, BUNGER A P.Interference fracturing: Nonuniform distributions of perforation clusters that promote simultaneous growth of multiple hydraulic fractures[J]. SPE Journal, 2015, 20(2): 384-395.
[26] DONTSOV E V.An approximate solution for a penny-shaped hydraulic fracture that accounts for fracture toughness, fluid viscosity and leak-off[J]. Royal Society Open Science, 2016, 3: 1-18.
[27] WU R, BUNGER A P, JEFFREY R G, et al.A comparison of numerical and experiemental results of hydraulic fracture growth into a zone of lower confining stress[R]. ARMA 08-267, 2008.
[28] ZIA H, LECAMPION B.Explicit versus implicit front advancing schemes for the simulation of hydraulic fracture growth[J]. International Journal of Numerical and Analytical Methods in Geomechanics, 2019: 1-16.
[29] BUNGER A P, ZHANG X, ROBERT G J.Parameters affecting the interaction among closely spaced hydraulic fractures[J]. SPE Journal, 2012, 17(1): 292-306.
[30] TANG H Y, WINTERFELD P H, WU Y S, et al.Integrated simulation of multi-stage hydraulic fracturing in unconventional reservoirs[J]. Journal of Natural Gas Science and Engineering, 2016, 36: 875-892.
[31] 胥云, 陈铭, 吴奇, 等. 水平井体积改造应力干扰计算模型及其应用[J]. 石油勘探与开发, 2016, 43(5): 780-786.
XU Yun, CHEN Ming, WU Qi, et al.Stress interference calculation model and its application in volume stimulation of horizontal wells[J]. Petroleum Exploration and Development, 2016, 43(5): 780-786.
[32] WU K, OLSON J E.Mechanisms of simultaneous hydraulic-fracture propagation from multiple perforation clusters in horizontal wells[J]. SPE Journal, 2016, 21(3): 1000-1008.
[33] CRUMP J B, CONWAY M W.Effects of perforation-entry friction on bottomhole treating analysis[J]. Journal of Petroleum Technology, 1988, 40(8): 1041-1048.
[34] CROUCH S L, STARFIELD A M.Boundary element methods in solid mechanics: With applications in rock mechanics and geological engineering[M]. New South Wales: Allen & Unwin, 1982.
[35] PEIRCE A P.Localized Jacobian ILU preconditioners for hydraulic fractures[J]. International Journal for Numerical Methods in Engineering, 2006, 65: 1935-1946.
[36] MEYER C D, BALSARA D S, ASLAM T D.A stabilized Runger-Kutta-Legendre method for explicit super-time-stepping of parabolic and mixed equations[J]. Journal of Computational Physics, 2013, 257: 594-626.
[37] LECAMPION B, PEIRCE A, DETOURNAY E.The impact of the near-tip logic on the accuracy and convergence rate of hydraulic fracture simulations compared to reference solutions[R]. Brisbane, Australia: The International Conference for Effective and Sustainable Hydraulic Fracturing, 2013.
[38] WU K, OLSON J, BALHOFF M, et al.Numerical analysis for promoting uniform development of simultaneous multiple-fracture propagation in horizontal wells[J]. SPE Production & Operations, 2017, 32(1): 41-50.
[39] LECAMPION B, DESROCHES J, WENG X, et al.Can we engineer better multistage horizontal completions? Evidence of the importance of near-wellbore fracture geometry from theory, lab and field experiments[R]. SPE 173363, 2015.
[40] YANG Z Z, YI L P, LI X G, et al.Pseudo-three-dimensional numerical model and investigation of multi-cluster fracturing within a stage in a horizontal well[J]. Journal of Petroleum Science and Engineering, 2018, 162: 190-213.
[41] CRAMER D, FRIEHAUF K, ROBERTS G, et al.Integrating DAS, treatment pressure analysis and video-based perforation imaging to evaluate limited entry treatment effectiveness[R]. SPE 194334, 2019.
Options
文章导航

/

〈 〉

版权所有：《石油勘探与开发》编辑部；《石油勘探与开发（英文）》编辑部
地址：北京市海淀区学院路20号，中国石油勘探开发研究院数据中心605室　

京ICP备15044687号    技术支持：北京玛格泰克科技发展有限公司