水平井多裂缝同步扩展数值模拟

doi:10.7623/syxb201512012

石油学报 ›› 2015, Vol. 36 ›› Issue (12): 1571-1579.DOI: 10.7623/syxb201512012

水平井多裂缝同步扩展数值模拟

曾青冬, 姚军

中国石油大学石油工程学院山东青岛 266580

收稿日期:2015-05-23 修回日期:2015-09-26 出版日期:2015-12-25 发布日期:2016-01-06
通讯作者: 姚军,男,1964年3月生,1984年获华东石油学院学士学位,2000年获石油大学(华东)博士学位,现为中国石油大学(华东)石油工程学院教授、博士生导师,主要从事油气田开发工程的教学和科研工作。Email:yaojunhdpu@126.com
作者简介:曾青冬,男,1987年11月生,2010年获中国石油大学(华东)学士学位,现为中国石油大学(华东)石油工程学院博士研究生,主要从事非常规油气藏压裂改造、油藏数值模拟方面的研究。Email:upc.zengqd@163.com
基金资助:
国家自然科学基金项目(No.51234007)、国家自然科学基金重大项目(No.51490654)、国家自然科学基金青年科学基金项目(No.51404291)和中央高校基本科研业务费专项资金项目(13CX02052A)资助。

Numerical simulation of multiple fractures simultaneous propagation in horizontal wells

Zeng Qingdong, Yao Jun

School of Petroleum Engineering, China University of Petroleum, Shandong Qingdao 266580, China

Received:2015-05-23 Revised:2015-09-26 Online:2015-12-25 Published:2016-01-06

摘要/Abstract

摘要：

为探究页岩气藏水力压裂复杂裂缝网络的形成机理,开展了水平井多裂缝同步扩展的数值模拟研究。考虑井筒和裂缝内流体流动,耦合岩体变形,建立了多裂缝同步扩展数学模型,利用裂缝尖端渐近解隐式求解裂缝扩展步长,并采用扩展有限元方法对模型进行了求解。通过与经典模型计算结果比较,验证了数学模型与算法的正确性,并将数值算法运用于2个算例。多裂缝同步扩展时,两侧裂缝的长度和宽度比中间裂缝的大。裂缝间距越小时,两侧裂缝偏离初始最大主地应力方向角度越大,中间裂缝宽度越小;地应力差越小时,两侧裂缝偏离初始最大主地应力方向角度越大。研究结果表明,应力阴影作用对多裂缝同步扩展轨迹和几何形态具有重要作用,且裂缝间距越小,地应力差越小,应力阴影作用越强,裂缝形态越偏离单裂缝扩展形态。

关键词: 分段压裂, 同步扩展, 扩展有限元, 尖端渐进解, 水平集方法

Abstract:

To investigate the formation mechanism of complex fracture network during hydraulic fracturing in shale gas reservoirs, numerical simulation was carried out on multiple fractures simultaneous propagation in horizontal wells. When considering fluid flow in wellbore and fractures, the mathematical model of multiple fractures simultaneous propagation was established based on coupled rock deformation. The fracture-tip asymptotic solution and implicit solution were obtained for calculating the step length of fracture propagation, and the extended finite element method was adopted for model solution. The correctness of mathematical model and algorithm were validated through a comparison with the calculation results of classic models, and numerical algorithm was applied in calculation cases. During multiple fractures simultaneous propagation, the length and width of the fractures on both sides were larger than those of central fractures. The smaller the fracture spacing was, the larger the angle of the fractures on both sides deviated from initial maximum principle stress would be, and the smaller the width of central fracture would be. The greater the geostress difference was, the larger the angle of the fractures on both sides deviated from initial maximum principle stress would be. Research results show that stress shadowing has significant influences on the trace and geometry of multiple fractures simultaneous propagation. The smaller the fracture spacing is, the smaller the geostress difference will be, the stronger the stress shadowing effect will be, and the more deviated from single-fracture propagation the fracture morphology will be.

Key words: staged fracturing, simultaneous propagation, extended finite element method, tip asymptote solutions, level set method

中图分类号:

TE357.1

曾青冬, 姚军. 水平井多裂缝同步扩展数值模拟[J]. 石油学报, 2015, 36(12): 1571-1579.

Zeng Qingdong, Yao Jun. Numerical simulation of multiple fractures simultaneous propagation in horizontal wells[J]. Acta Petrolei Sinica, 2015, 36(12): 1571-1579.

参考文献

[1] Fisher M K,Heinze J R,Harris C D,et al.Optimizing horizontal completions in the Barnett shale with microseismic fracture mapping[J].Journal of Petroleum Technology,2005,57(3):41-43.
[2] 姚军,孙海,黄朝琴,等.页岩气藏开发中的关键力学问题[J].中国科学:物理学力学天文学,2013,43(12):1527-1547.
Yao Jun,Sun Hai,Huang Zhaoqin,et al.Key mechanical problems in the development of shale gas reservoirs[J].Scientia Sinica Physica,Mechanica & Astronomica, 2013,43(12):1527-1547.
[3] 郭建春,苟波,任山,等.川西页岩-砂岩交互水平井压裂参数优化设计[J].石油学报,2014,35(3):511-518.
Guo Jianchun,Gou Bo,Ren Shan,et al.Fracturing parameters optimization design for horizontal well with shale and sandstone interbeds in western Sichuan[J].Acta Petrolei Sinica,2014,35(3):511-518.
[4] Geertsma J,Klerk F D.A rapid method of predicting width and extent of hydraulically induced fractures[J].Journal of Petroleum Technology,1969,21(12):1571-1581.
[5] Perkins T K,Kern L R.Width of hydraulic fractures[J].Journal of Petroleum Technology,1961,13(9):937-949.
[6] Nordgren R P.Propagation of a vertical hydraulic fracture[J].Society of Petroleum Engineers Journal,1972,12(4):306-314.
[7] 时贤,程远方,蒋恕,等.页岩储层裂缝网络延伸模型及其应用[J].石油学报,2014,35(6):1130-1137.
Shi Xian,Cheng Yuanfang,Jiang Shu,et al.Simulation of complex fracture network propagation and its application for shale gas reservoir[J].Acta Petrolei Sinica,2014,35(6):1130-1137.
[8] Cheng Y.Mechanical interaction of multiple fractures-exploring impacts on selection of the spacing/number of perforation clusters on horizontal shale gas wells[J].Society of Petroleum Engineers Journal,2012,17(4):992-1001.
[9] Geilikman M,Xu G S,Wong S W.Interaction of multiple hydraulic fractures in horizontal wells[R].SPE 163982,2013.
[10] Wu K,Olson J E.Simultaneous multifracture treatments:fully coupled fluid flow and fracture mechanics for horizontal wells[J].SPE Journal,2015,20(2):337-346.
[11] 曾青冬,姚军.基于扩展有限元的页岩水力压裂数值模拟[J].应用数学和力学,2014,35(11):1239-1248.
Zeng Qingdong,Yao Jun.Numerical simulation of shale hydraulic fracturing based on the extended finite element method[J].Applied Mathematics and Mechanics,2014,35(11):1239-1248.
[12] Gordeliy E,Peirce A.Coupling schemes for modeling hydraulic fracture propagation using the XFEM[J].Computer Methods in Applied Mechanics and Engineering,2013,253:305-322.
[13] Taleghani A D.Analysis of hydraulic fracture propagation in fractured reservoirs:an improved model for the interaction between induced and natural fractures[D] .Austin:University of Texas at Austin,2009.
[14] Ren Q W,Dong Y W,Yu T T.Numerical modeling of concrete hydraulic fracturing with extended finite element method[J].Science in China Series E:Technological Sciences,2009,52(3):559-565.
[15] Lecampion B.An extended finite element method for hydraulic fracture problems[J].Communications in Numerical Methods in Engineering,2009,25(2):121-133.
[16] Gordeliy E,Peirce A.Implicit level set schemes for modeling hydraulic fractures using the XFEM[J].Computer Methods in Applied Mechanics and Engineering,2013,266:125-143.
[17] Linkov A M.Universal asymptotic umbrella for hydraulic fracture modeling[EB/OL].[2015-01-10].http://de.arxiv.org/ftp/arxiv/papers/1404/1404.4165.pdf.
[18] Zilian A,Fries T P.A localized mixed-hybrid method for imposing interfacial constraints in the extended finite element method (XFEM)[J].International Journal for Numerical Methods in Engineering,2009,79(6):733-752.
[19] Peirce A.A Hermite cubic collocation scheme for plane strain hydraulic fractures[J].Computer Methods in Applied Mechanics and Engineering,2010,199(29/32):1949-1962.
[20] Mousavi S E,Grinspun E,Sukumar N.Harmonic enrichment functions:a unified treatment of multiple,intersecting and branched cracks in the extended finite element method[J].International Journal for Numerical Methods in Engineering,2011,85(10):1306-1322.
[21] Spence D A,Sharp P.Self-similar solutions for elastohydrodynamic cavity flow[J].Proceedings of the Royal Society of London.Series A,Mathematical and Physical Sciences,1985,400(1819):289-313.

水平井多裂缝同步扩展数值模拟

Numerical simulation of multiple fractures simultaneous propagation in horizontal wells

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 4

编辑推荐

Metrics

本文评价

[1]	王志刚. 涪陵大型海相页岩气田成藏条件及高效勘探开发关键技术[J]. 石油学报, 2019, 40(3): 370-382.
[2]	杨正明, 刘学伟, 张仲宏, 周拓, 赵淑霞. 致密油藏分段压裂水平井注二氧化碳吞吐物理模拟[J]. 石油学报, 2015, 36(6): 724-729.
[3]	杨正明, 张仲宏, 刘学伟, 田文博, 徐庆岩. 低渗/致密油藏分段压裂水平井渗流特征的物理模拟及数值模拟[J]. 石油学报, 2014, 35(1): 85-92.
[4]	曾凡辉, 郭建春, 刘恒, 尹建. 致密砂岩气藏水平井分段压裂优化设计与应用[J]. 石油学报, 2013, 34(5): 959-968.