基于有限体积法的裂缝性油藏两相流动模型

doi:10.7623/syxb201808008

石油学报 ›› 2018, Vol. 39 ›› Issue (8): 924-936.DOI: 10.7623/syxb201808008

基于有限体积法的裂缝性油藏两相流动模型

唐潮¹, 陈小凡¹, 杜志敏¹, 刘东晨¹, 汤连东¹, 魏嘉宝²

1. 西南石油大学油气藏地质及开发工程国家重点实验室四川成都 610500;
2. 中国石油化工股份有限公司胜利油田分公司河口采油厂山东东营 257015

收稿日期:2017-11-24 修回日期:2018-06-06 出版日期:2018-08-25 发布日期:2018-09-04
通讯作者: 陈小凡,男,1968年6月生,1988年获西南石油学院学士学位,1991年获西南石油学院硕士学位,现为西南石油大学石油与天然气工程学院副教授,主要从事油藏工程及数值模拟方面研究与教学工作。Email:382632351@qq.com
作者简介:唐潮,男,1988年4月生,2010年获西南石油大学学士学位,现为西南石油大学博士研究生,主要从事渗流力学及油气藏数值模拟研究。Email:304410863@qq.com
基金资助:
国家科技重大专项（2016ZX05010-002-002）、国家自然科学基金项目（No.51474179）资助。

Fractured reservoir two-phase flow model based on finite volume method

Tang Chao¹, Chen Xiaofan¹, Du Zhimin¹, Liu Dongchen¹, Tang Liandong¹, Wei Jiabao²

1. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Sichuan Chengdu 610500, China;
2. Hekou Oil Production Plant, Sinopec Shengli Oilfield Company, Shandong Dongying 257015, China

Received:2017-11-24 Revised:2018-06-06 Online:2018-08-25 Published:2018-09-04

摘要/Abstract

摘要：

针对裂缝性储层中裂缝分布的强非均质性和多尺度性，提出了一种多尺度离散裂缝两相流动模型的数值求解方法。该方法将n维模型中大尺度裂缝降维处理成n-1维实体，并基于格心形式的有限体积法对模型进行了离散，建立了数值求解格式，这些处理方法对降低求解时间及内存需求有显著影响。该方法具有严格的物理意义，保证了模型在局部范围内仍然遵循守恒定律。通过建立上游权值的有界高阶迎风格式，解决了一阶迎风格式所产生的假扩散现象以及大尺度裂缝与基质-微裂缝界面处饱和度突变引起的非物理振荡。使用IMPES顺序求解法对模型进行求解，通过结果对比证明了模型准确性。

关键词: 裂缝建模, 有限体积法, 有界迎风格式, 油水两相流, 多尺度裂缝

Abstract:

Aiming at the strong heterogeneity and multiscale distribution of fractures in fractured reservoirs, a numerical solution method is proposed for the two-phase flow model for multiscale discrete fractures. In this method, the large-scale fractures in n-dimensional model turn into n-1 dimensional entity after dimensionality reduction. Meanwhile, the cell-centered finite volume method is used to discrete model, and then the numerical solution format is established. This method has significant impacts on reducing solving time and memory requirements, characterized by strict physical significance, ensuring the model following the law of conservation on a regional scale. The bounded high-order upwind scheme of upstream weight is established to solve the pseudo-diffusion phenomenon caused by the first-order upwind scheme and the unphysical oscillation due to the sudden change of saturation at the interface between large-scale fractures and matrix-micro fractures. Moreover, the IMPES sequence solution method is adopted for model solving, and the model accuracy is proven by result contrast.

Key words: fracture modeling, finite volume method, bounded upwind scheme, oil-water two-phase flow, multi-scale facture

中图分类号:

TE312

唐潮, 陈小凡, 杜志敏, 刘东晨, 汤连东, 魏嘉宝. 基于有限体积法的裂缝性油藏两相流动模型[J]. 石油学报, 2018, 39(8): 924-936.

Tang Chao, Chen Xiaofan, Du Zhimin, Liu Dongchen, Tang Liandong, Wei Jiabao. Fractured reservoir two-phase flow model based on finite volume method[J]. Acta Petrolei Sinica, 2018, 39(8): 924-936.

参考文献

[1] GILMAN J R,KAZEMI H.Improvements in simulation of naturally fractured reservoirs[J].Society of Petroleum Engineers Journal,1983,23(4):695-707.
[2] DEAN R H,LO L L.Simulations of naturally fractured reservoirs[J].SPE Reservoir Engineering,1988,3(2):638-648.
[3] STALGOROVA E,BABADAGLI T.Field scale modeling of tracer injection in naturally fractured reservoirs using the random-walk particle-tracking simulation[J].SPE Journal,2012,17(2):580-592.
[4] 姚军,王子胜,张允,等.天然裂缝性油藏的离散裂缝网络数值模拟方法[J].石油学报,2010,31(2):284-288.
YAO Jun,WANG Zisheng,ZHANG Yun,et al.Numerical simulation method of discrete fracture network for naturally fractured reservoirs[J].Acta Petrolei Sinica,2010,31(2):284-288.
[5] YOUNG L C.A finite-element method for reservoir simulation[J].Society of Petroleum Engineers Journal,1981,21(1):115-128.
[6] KIM J G,DEO M D.Finite element,discrete-fracture model for multiphase flow in porous media[J].AIChE Journal, 2000,46(6):1120-1130.
[7] KARIMI-FARD M,FIROOZABADI A.Numerical simulation of water injection in fractured media using the discrete-fracture model and the Galerkin method[J].SPE Reservoir Evaluation & Engineering,2003,6(2):117-126.
[8] MONTEAGUDO J E P,FIROOZABADI A.Control-volume method for numerical simulation of two-phase immiscible flow in two-and three-dimensional discrete-fractured media[J].Water Resources Research,2004,40(7):W07405.
[9] BASTIAN P,CHEN Zhangxin,EWING R E,et al.Numerical simulation of multiphase flow in fractured porous media[M]//CHEN Zhangxin,EWING R E,SHI Zhongci.Numerical Treatment of Multiphase Flows in Porous Media.Berlin Heidelberg:Springer,2000:50-68.
[10] GEIGER S,ROBERTS S,MATTHÄI S K,et al.Combining finite volume and finite element methods to simulate fluid flow in geologic media[J].ANZIAM Journal,2003,44:C180-C201.
[11] CHEN Zhangxin,HUAN Guanren,MA Yuanle.Computational methods for multiphase flows in porous media[M].Philadelphia:Society for Industrial and Applied Mathematics,2006.
[12] JIANG Jiamin,YOUNIS R M.A multimechanistic multicontinuum model for simulating shale gas reservoir with complex fractured system[J].Fuel,2015,161:333-344.
[13] 张庆福,黄朝琴,姚军,等.基于多尺度混合有限元的离散裂缝两相渗流数值模拟[J].科学通报,2017,62(13):1392-1401.
ZHANG Qingfu,HUANG Zhaoqin,YAO Jun,et al.Two-phase numerical simulation of discrete fracture model based on multiscale mixed finite element method[J].Chinese Science Bulletin,2017,62(13):1392-1401.
[14] 任宗孝,吴晓东,屈展,等.倾斜裂缝水平井非稳态压力模型[J].石油学报,2017,38(9):1059-1065.
REN Zongxia,WU Xiaodong,QU Zhan,et al.Transient pressure model for horizontal wells with inclined fractures[J].Acta Petrolei Sinica,2017,38(9):1059-1065.
[15] 郭小哲,王晶,刘学锋.页岩气储层压裂水平井气-水两相渗流模型[J].石油学报,2016,37(9):1165-1170.
GUO Xiaozhe,WANG Jing,LIU Xuefeng.Gas-water two phase porous flow model of fractured horizontal well in shale gas reservoir[J].Acta Petrolei Sinica,2016,37(9):1165-1170.
[16] 刘志帆,刘志强,施安峰,等.水平井水力压裂影响参数的数值模拟[J].地球科学,2017,42(8):1394-1402.
LIU Zhifan,LIU Zhiqiang,SHI Anfeng,et al.Numerical simulation on influence parameters of horizontal well hydraulic fracturing[J].Earth Science,2017,42(8):1394-1402.
[17] 李道伦,杨景海,闫术,等.致密油大规模多段压裂水平试井解释及外区渗透率对试井曲线的影响[J].地球科学,2017,42(8):1324-1332.
LI Daolun,YANG Jinghai,YAN Shu,et al.Numerical well test interpretation of massive multistage fractured horizontal wells in tight oil reservoirs and effect of permeability of exterior region on well test curves[J].Earth Science,2017,42(8):1324-1332.
[18] 王濡岳,丁文龙,龚大建,等.渝东南-黔北地区下寒武统牛蹄塘组页岩裂缝发育特征与主控因素[J].石油学报,2016,37(7):832-845.
WANG Ruyue,DING Wenlong,GONG Dajian,et al.Development characteristics and major controlling factors of shale fractures in the Lower Cambrian Niutitang Formation,southeastern Chongqing-northern Guizhou area[J].Acta Petrolei Sinica,2016,37(7):832-845.
[19] 李颖川.采油工程[M].北京:石油工业出版社,2009.
LI Yingchuan.Oil recovery engineering[M].Beijing:Petroleum Industry Press,2009.
[20] 孔祥言.高等渗流力学[M].2版.北京:中国科学技术大学出版社,2010.
KONG Xiangyan.Advanced mechanics of flow in porous medias[M].2nd ed.Beijing:China University of Science and Technology Press,2010.
[21] WARREN J E,ROOT P J.The behavior of naturally fractured reservoirs[J].Society of Petroleum Engineers Journal,1963,3(3):245-255.
[22] MCCAIN W D.The properties of petroleum fluids[M].Tulsa,OK:Penn Well Books,1974.
[23] MONTEAGUDO J E P,FIROOZABADI A.Control-volume model for simulation of water injection in fractured media:incorporating matrix heterogeneity and reservoir wettability effects[J].SPE Journal,2007,12(3):355-366.
[24] KIM J G,DEO M D.Comparison of the performance of a discrete fracture multiphase model with those using conventional methods:Proceedings of the 199915th Symposium on Reservoir Simulation[C].Houston,TX,USA:SPE,1999.
[25] ZHANG Ruihan,ZHANG Liehui,LUO Jianxin,et al.Numerical simulation of water flooding in natural fractured reservoirs based on control volume finite element method[J].Journal of Petroleum Science and Engineering,2016,146:1211-1225.
[26] 李潮流,胡法龙,侯雨庭,等.基于有限元的致密砂岩储层电阻率特性模拟[J].石油学报,2016,37(6):787-795.
LI Chaoliu,HU Falong,HOU Yuting,et al.Simulation of electrical resistivity characteristics of tight sands reservoir by FEM[J].Acta Petrolei Sinica,2016,37(6):787-795.
[27] 李允.油藏模拟[M].东营:石油大学出版社,1998.
Li Yun.Reservoir simulation[M].Dongying:Petroleum University Press,1998.
[28] PATANKAR S V.Numerical heat transfer and fluid flow[M].Boca Raton:CRC Press,1980:125-126.
[29] DARWISH M,MOUKALLED F.A new approach for building bounded skew-upwind schemes[J].Computer Methods in Applied Mechanics and Engineering,1996,129(3):221-233.
[30] SPEKREIJSE S.Multigrid solution of monotone second-order discretizations of hyperbolic conservation laws[J].Mathematics of Computation, 1987,49(179):135-155.
[31] BARTH T,JESPERSEN D.The design and application of upwind schemes on unstructured meshes:27th Aerospace Sciences Meeting[C].Reno, NV:AIAA,1989.
[32] LEONARD B P.The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection[J].Computer Methods in Applied Mechanics and Engineering,1991,88(1):17-74.
[33] MOUKALLED F,MANGANI L,DARWISH M.The finite volume method in computational fluid dynamics:an advanced introduction with OpenFOAM○ R and matlab[M].Switzerland:Springer International Publishing,2016.

基于有限体积法的裂缝性油藏两相流动模型

Fractured reservoir two-phase flow model based on finite volume method

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 3

编辑推荐

Metrics

本文评价

[1]	苏皓, 雷征东, 李俊超, 龚斌, 赵文琪, 鞠斌山, 张荻萩, 秦鹤. 储集层多尺度裂缝高效数值模拟模型[J]. 石油学报, 2019, 40(5): 587-593,634.
[2]	熊兆洪　李振林　宫　敬　李　克. 埋地管道小泄漏模型及数值求解[J]. 石油学报, 2012, 33(3): 493-498.
[3]	姚军, 吴明录, 戴卫华, 朱玉国. 流线数值试井解释模型[J]. 石油学报, 2006, 27(3): 96-99.