用户名: 密码: 验证码:
裂缝性油藏离散裂缝网络模型数值模拟研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
基于双重连续介质理论的裂缝性油藏数值模拟方法,当裂缝系统发育较差时其连续性与实际相差较大,因而这类方法的适应性受到限制。实际上,裂缝系统是一个分隔基质岩块系统的网络系统,可视为离散的裂缝网络系统,更能准确描述裂缝性油藏的非均质性。论文以离散裂缝网络模型为基础,利用有限元方法,突破了只能求解不可压缩流体的限定,成功解决了二维离散裂缝网络模型的微可压缩油水两相流动的数值模拟问题。首先,在对裂缝发育特征分析的基础上研究了天然裂缝、人工裂缝及其开启、半开启、充填、闭合(潜在)裂缝的识别特征,提出了裂缝参数的表征方法,基于体积法和Parsons的平板流动理论形成了裂缝的孔隙度和渗透率的计算方法,为裂缝性油藏的数值模拟奠定了基础;第二,考虑地质上所描述的裂缝发育分布情况以及裂缝的形态、充填性、渗流特性等因素建立了离散裂缝网络物理模型和微可压缩油水两相流动的数学模型;第三,提出了离散裂缝网络数学模型的有限元求解方法,分别对基岩和裂缝系统采用三角形单元和线单元进行离散并推导出了其有限元插值函数,通过采用加权余量法得到了离散裂缝网络数值模拟的压力和饱和度方程的等效积分“弱”形式,形成了离散裂缝网络数值模拟的单元特性矩阵和列阵,建立了离散裂缝网络的有限元数值模型;第四,通过梯度法、共轭梯度法和预条件共轭梯度法的对比分析,依据收敛速度、需求存储空间等特点优选了预条件共轭梯度法,并形成了离散裂缝网络模型的有限元数值模型的预条件共轭梯度求解方法;最后,基于Windows XP操作系统和多物理场耦合分析软件(COMSOL Multiphysics)编制了裂缝性油藏离散裂缝网络数值模拟器,并对各种典型的裂缝性油藏概念模型进行了数值模拟研究,表明离散裂缝网络模型数值模拟方法能够准确描述裂缝属性对流动的影响,其油藏动态更加符合真实情况,具有较强的实用性和广阔的应用前景。
When fracture development is poor, its continuity has a very great difference compared with the reality. The numerical simulation method of fractured reservoir based on double porosity medium theory is limited when applying. In fact, fracture system is a network system separating matrix block and it can be regarded as discrete fracture network (DFN) system, which describe the heterogeneity of fractured reservoir accurately. On the basis of DFN model, the thesis overcomes the limit of only solving incompressible fluid and solves numerical simulation problems of slightly compressible oil-water two-phase flow with finite element method (FEM) successfully. Firstly, the fracture identification characteristics of natural fracture, artificial fracture, unlocked and half-unlocked fracture, fracture filling and fracture closure was studied based on the analysis of field data, such as FMI imaging logging etc. The characterization method of fracture parameters is proposed. And the calculation method of porosity and permeability is formed respectively using Volume Method and Parsons Plate Flow Theory, which provides basis for fractured reservoir numerical simulation. Secondly, according to fracture development and distribution in geological description mathematical model, fracture morphology, fracture filling, fracture seepage characteristics and other factors, the thesis has built DFN physical model and slightly compressible oil-water two-phase flow mathematical model. Thirdly, the solving method of mathematical model on DFN was proposed with FEM. Matrix is discretized using linear triangular elements; fractures are discretized by linear line elements. And finite element shape function is derived. The equivalent integral weak form of pressure and saturation equation of DFN numerical simulation is established using weighted residual method. Element characteristic matrix and array of DFN numerical simulation is formed, and finite element numerical model of DFN is built. Fourthly, according to convergence rate, demanding storage space and so on, preconditioned conjugate gradient method (PCG) is the optimal method after the comparison and analysis with gradient method and conjugate gradient method. And then PCG solving method of DFN finite element numerical model is formed. Finally, DFN numerical simulator of carbonate fractured reservoir was drawn up by COMSOL Multiphysics software at the circumstance of Windows XP. Then the numerical simulation on different typical conceptual models of fractured reservoir is studied. The results show that effects of fracture attribute for the fluid flow is accurately accounted by DFN numerical simulation method, and its reservoir performance is more close to realities. This method has strong practicability and wide application prospect.
引文
[1] M.Karimi-Fard, RERI, A.Firoozabadi. Numerical simulation of water injection in 2d fractured media using discrete-fracture model. the 2001 SPE Annual Technical Conference and Conference [C]. New Orleans, Louisiana: SPE 71615, 2001.
    [2] Arnaud Lange, Remy Basquet, Bernard Bourbiaux. Hydraulic characterization of faults and fractures using a dual medium discrete fracture network simulator. the 11th Abu Dhabi International Petroleum Exhibition and Conference[C]. Abu Dhabi, U.A.E: SPE 88675, 2004.
    [3] Bill Dershowitz, Paul Lapointe, Thorsten Eiben et al. Intergration of discrete feature network methods with conventional simulator approaches. the 1998 Annual Technical Conference and Exhibition[C]. New Orleans: SPE 62489,1998.
    [4] Jong Gyun Kim, Milind D. Deo. Comparison of the performance of a discrete fracture multiphase model with those using conventional methods. the 1999 Reservoir Simulation Symposium, Houston[C]. Texas: SPE 51928, 1999.
    [5] S. H. Lee, L. J. Durlofsky, M. F. Lough et al. Finite difference simulation of geologically complex reservoirs with tensor permeabilities. the 1997 SPE Reservoir Simulation Symposium[C]. Dallas: SPE 52637, 1997.
    [6] Ahamed Ouenes, Lee J.Hartley. Integrated fractured reservoir modeling using both discrete and continuum approaches. the 2000 SPE Annual Technical Conference and Exhibition[C]. Dallas, Texas: SPE 62939, 2000.
    [7] S. H. Lee, C. L. Jensen, M. F. Lough. An efficient finite difference model for flow in a reservoir with multiple length-scale fracures. the SPE Annual Technical Conference and Exhibition[C]. Houston, U.S.A: SPE 56752, 1999.
    [8] F. M. Verga, G. Giglio, Di Torino et al.: Calibration of fractured reservoirs with dynamic data. the SPE Reservoir Simulation Symposium[C]. Houston, Texas: SPE 66395, 2001.
    [9] S. Sarda, L. Jeannin, B. Bourbiaux. Hydraulic characterization of fractured reservoirs: simulation on discrete fracture models. the SPE Reservoirs Simulation Symposium[C]. Houston, Texas: SPE 66398, 2001.
    [10] M. Ziad Saghir, Mohamed A. Edali, M. R. Islam. Rigorous modeling of fractured in a porous medium. the SPE Latin American and Caribbean Petroleum Engineering Conference[C]. Buenos Aires, Argentina: SPE 69574, 2001.
    [11] Bernard Bourbiaux, Remy Basquet, Marie-Christine. An integrated workflow toaccount for multi-scale fractures in reservoir simulation models: implementation and benefits. the 10th Abu Dhabi International Exhibition and Conference[C]. Abu Dhabi: SPE 78489, 2002.
    [12] M. Karimi-Fard, L. J. Durlofsky, K. Aziz. An efficient discrete-fracture model applicable for general-purpose reservoir simulators. the 2003 SPE Reservoir Simulation Symposium[C]. Houston: SPE 88812, 2003.
    [13] Hugo Araujo, Pablo Lacentre, Tomas Zapata. Dynamic behavoir of discrete fracture network models. the 2004 SPE International Petroleum Conference[C]. Puebla, Mexico: SPE 91940, 2004.
    [14] J. E. P. Monteagudo, A. Firoozabadi. Numerical simulation of water injection in disconnected and connected fractured media using jacobian-free fully implicit control volume method. the 2004 SPE/DOE 14th Symposium on Improved Oil Recovery[C]. Tulsa, Oklahoma, U.S.A.: SPE 89449, 2004.
    [15]王洪涛,王恩志,金宜英.裂隙岩体渗流耦合数值分析方法及其工程应用[J].南京水利科学研究院水利水运科学研究,1998,(4):320-328.
    [16]王洪涛,王恩志.岩体主干裂隙系统三维非稳定渗流分析模型[J].水动力学研究与进展,1998,13(2):206-213.
    [17] J. E. Warren, P.J. Root. The behavior of naturally fractured reservoir. the Fall Meeting of the Society of Petroleum Engineers[C]. Los Angeles: SPE 426, 1962.
    [18] Roy H. Yamamoto, J. B. Padgett, Wayne T. Ford. Compositional reservoir simulator for fissured systems-the single-block model. SPE 44th Annual Fall Meeting[C]. Denver, Colo: SPE 2666, 1969.
    [19] H. Kazemi, L. S. Merrill. JR, K. L. Porterfield. Numerical simulation of water-oil flow in naturally fractured reservoirs. Simulation of Reservoir Performance[C]. Loa Angeies, 1976.
    [20] Rossen. Simulation of naturally fractured reservoirs with semi-implicit source terms. the SPE-AIME Fo~h SymPOsiVm on Numerical Simulation of Reservoir Performance[C]. Los Angeles, 1976.
    [21] Saidi. Simulations of naturally fractured reservoirs. the Reservoir Simulatlon Symposium[C]. San Francleco, CA, 1983.
    [22] Kent Thomas, Thomas N. Dixon, Ray G. Pierson. Fractured reservoir simulation. the 1980 SPE Annual Technical Reference and Exibition[C]. Dallas, 1980.
    [23] James R. Gilman, Hossein Kazemi. Improvements in simulation of naturallyfractured reservoirs. the 1982 SPE Reservoir Simulation Symposium[C]. New Orleans: SPE 10511, 1982.
    [24] Dimitrie Bossie-Codreanu, Paul R. Bia, Jean-Claude Sabathier. The“checker model”an improvement in modeling naturally fractured reservoirs with a tridimensional, triphasic, black-oil numerical model. the 1982 SPE Annual Technical Conference and Exhibition[C]. New Orleans: SPE 10977, 1982.
    [25] R. H. Dean, L. L. Lo. Simulation of naturally fractured reservoirs. the 1986 SPE Intl Meeting on Petroleum Engineering[C]. Beijing: SPE 14110, 1986.
    [26] Aanda M. Wijesinghe, Wilmer E. Culham. Single-well pressure testing solutions for naturally fractured reservoirs with arbitrary fracture connectivity. the 59th Annual Technical Conference and Exhibition[C]. Houston, Texas: SPE 13055, 1984.
    [27] A. C. Hill, G. W. Thomas. A new approach for simulating complex fractured reservoirs. the SPE 1985 Middle East Oil Technical Conference and Exhibition[C]. Bahrain: SPE 13537, 1985.
    [28] D. Waldren, A. F. Corrigan. An engineering and geological review of the problems encountered in simulating naturally fractured reservoirs. the SPE 1985 Middle East Oil Technical Conference and Exhibition[C]. Bahrain: SPE 13717, 1985.
    [29] F. Sonier, P. Souillard, F. T. Blaskovich. Numerical simulation of naturally fractured reservoirs. the 1986 SPE Annual Technical Conference and Exhibition[C]. New Orlean: SPE 15627, 1986.
    [30] Ishimoto. Improved matrix, fracture fluid transfer function in dual porosity models. the SPE International Meeting on Petroleum Engineering[C]. Tianjin, China: Paper SPE 17599, 1988.
    [31] T. A. Jones, S. O. Wooten, T. J. Kaluza. Single-phase flow through natural fractures. the 63rd Annual Technical Conference and Exhibition[C]. Houston, TX: SPE 18175, 1968.
    [32] Abbas Firoozabadi, L. Kent Thomas. Sixth spe comparative solution project: dual-porosity simulators. the 1989 SPE Symposium on Reservoir Simulation[C]. Houston: SPE 18741, 1989.
    [33] J. Kamath, S. H. Lee, C. L. Jensen et al. Modeling fluid flow in complex naturally fractured reservoirs. the SPE India Oil and Gas Conference and Exhibition[C]. New Delhi, India: SPE 39547, 1998.
    [34] Gurpinar, C. A. Kossack. Realistic numerical models for fractured reservoirs. the 2000 SPE International Petroleum Conference and Exhibition[C]. Villahermosa, Mexico: SPE 59041, 2000.
    [35] Lingli Wei, Jerry Hadwin, Keith Rawnsley ec al. Discriminating fracture patterns in fractured reservoirs by pressure transient tests. the 1998 SPE Annual Technical Conference and Exhibition[C]. New Orleans, Louisiana: SPE 49233, 1998.
    [36] Richard O. Baker, Frank Kuppe. Reservoir characterization for naturally fractured reservoirs. the 2000 SPE Annual Technical Conference and Exhibition[C]. Dallas, Texas: SPE 63286, 2000.
    [37] Javandel I, Witherspoon P A. Application of the finite element method to transient flow in porous media [J]. SPEJ, 1968.9:241-252 (SPE 2052, 1967.9).
    [38] McMichael C L, Thomas G W. Reservoir simulation by galerkin's method[J]. SPE 3558, 1973.
    [39] Langsrud Qyvind, Computas A S. Simulation of two-phase flow by finite element methods[J]. SPE 5725, 1976.
    [40] Spivak A, Price H 5, Settari A. Solution of the equations for multidimensional, two-phase, immiscible flow by variational methods [J]. SPEJ, 1977: 27-41 (SPE 5723, 1976).
    [41] Dalen Vilgeir. Simplified finite-element models for reservoir flow problems [J]. SPE 7196, 1979.
    [42] Young Larry C. A finite-element method for reservoir simulation [J]. SPE 7413, 1981.
    [43] Masato Tauchia, Kazuo Kohda, Yutaka Suzukawa. A well model for finite element reservoir simulators [J]. SPE 12936, 1984.
    [44] Bhatia K S, Advani S H, Lee J K. Finite element representation of two-phase fluid flow through a naturally fractured reservoir[J]. SPE 19069, 1989.
    [45] Nobuo Morita. Transient finite element code: a versatile tool for well performance analyses[J]. SPE 25878, 1993.
    [46] Murzxnko V. Finite element method for steady fluid flow in reservoir containing a fracture[J]. SPE 26306, 1994.
    [47] Jong Gyun Kim, Milind D Deo. Finite element, discrete-fracture model for multiphase flow in porous media[J]. AIChE Journal, 2000, 46(6):1120-1130.
    [48] Lewis R W, Verner E A, Zienkiewicz O C. A finite element approach to two-phase flow in porous media[M]. Finite Element in Fluids, John Wiley&Sons Ltd, London, 1975.
    [49] Kebaili A L, Thomas G W. A two-phases coning model using alternating direction galerkin procedure[J]. SPE 5724, 1976.
    [50]袁益让,王文洽.关于油水两项渗流平面弹性驱动问题的有限元方法[J].石油学报,1980,1(4):65-76.
    [51]袁益让.二相渗流数值模拟的一类有限元格式及其理论分析[J].科学通报,1984,4:193-197.
    [52]袁益让.油、水二相渗流驱动问题的变网格有限元方法及其理论分析[J].中国科学,1986,A(2):135-148.
    [53]梁栋.两项渗流问题特征混合元数值模拟方法[J].中国科学,1991,A(5):478-485.
    [54]程林松,郎兆新.水平井油—水两相渗流的有限元方法[J].水动力学研究与进展,1995,10(3):309-314.
    [55]周竹眉,郎兆新.水平井油藏数值模拟的有限元方法[J].水动力学研究与进展,1996,11(3):261-270.
    [56]程林松,李春兰,郎兆新.裂缝性底水油藏水平井三维油水两相有限元数值模拟方法[J].石油勘探与开发,1998,25(2):41-45.
    [57]程林松,李春兰,郎兆新.裂缝油藏三维油、水两相渗流的有限元方法[J].大庆石油地质与开发,1998,17(1):26-27, 33.
    [58]张子香,宋洪才,范江.两维两相不稳定渗流变分有限元数值模拟方法[J].石油学报,1992,13(3):79-86.
    [59]宋洪才,张继芬,王江.油藏数值模拟中的变分有限元方法[J].大庆石油学院学报,1996,20(3):97-101.
    [60]谢海兵,马远乐,桓冠仁等.非结构网格油藏数值模拟方法研究[J].石油学报,2001,22(1):63-66.
    [61]蒋廷学,郎兆新,单文文.低渗透油藏压裂井动态预测的有限元方法[J].石油学报,2002,23(5):53-58.
    [62]向晓丽,张玺.有限元技术在预测储层孔隙度、渗透率中的应用[J].油气采收率技术,1999,6(2):49-60.
    [63]袁奕群,袁庆峰.黑油模型在油田开发中的应用[M].北京:石油工业出版社,1995.
    [64]宋考平,宋洪才,吴文样.油藏数值模拟理论基础[M].北京:石油工业出版社,1996.
    [65]陈月明.油藏数值模拟基础[M].北京:石油工业出版社,1989.
    [66]韩大匡,陈钦雷,闰存章.油藏数值模拟基础[M].北京:石油工业出版社,1993.
    [67]王鸿勋,张士诚.水力压裂设计数值计算方法[M].北京:石油工业出版社,1998.
    [68]范家齐,樊俊才.实用有限元法[M].哈尔滨:黑龙江科学技术出版社,1993.
    [69] Cheung Y K, Yeo M F.实用有限元分析导论[M].北京:人民交通出版社,1982.
    [70]蒋孝煌.有限元法基础[M].北京:清华大学出版社,1992.
    [71]署恒木,仝兴华.工程有限元法[M].东营:石油大学出版社,2003.
    [72]张鸿庆,王鸣.有限元的数学理论[M].北京:科学出版社,1991.
    [73] Nobuo Morita. Transient finite element code: a versatile tool for well performance analyses[J].SPE 25878,1993.
    [74] Kuland B.R. Dean S.L., Ward dr B.J. Fracture core analysis: interpretation, logging and use of natural and induced fractures in core[J]. AAPG Methods in Exploration Series, Tulsa, Oklahoma, 1990, (8).
    [75]裘怿楠,陈子琪.中国油藏管理技术手册.油藏描述[M].北京:石油工业出版社,1996.
    [76] Ramsay J.G., Huber M.I. The techniques of modern structural geology, volume 2: folds and fractures[M]. London: Academic Press, 1987.
    [77] Zhao G., Johoson A.M.. Sequential and incremental formation of conjugate sets of faults[J]. Journal of Structural Geology, 1991, 13(8): 887~895.
    [78] Nelson R A. Geological analysis of naturally fractured reservoirs[M]. Houston: Gulf Publishing Company, 1985.
    [79] NARR W, LERCHE L. A method for estimating subsurface fracture density in core[J]. AAPG Bull, 1984, 68(5): 637-648.
    [80] Brown, A., Davies, M., Nicholson, H. et al. The machar field, unlocking the potential of a north sea chalk field. the Offshore Europe Conference[C]. Aberdeen, Scotland: SPE56974, 1999.
    [81] T. Tamagawa, T. Matsuura, T. Anraku, K. Tezuka, T. Namikawa. Construction of fracture network model using static and dynamic data. Society of Petroleum Engineers Annual Technical Conference and Exhibition[C]. Texas, USA: SPE77741, 2002.
    [82] Wu H. Pollard D D. Propagation of a set of opening-mode fractures in layered brittle materials under uniaxial strain cycling[J]. J. Geophs. Res., 1992, 97(B3): 3381~3396.
    [83] Wu H. Pollard D D. Effect of strain rate on a set of fractures[J]. J. Rock. Mech. Min. Sci.&Geomech, Abstr. 1993, 30(7): 869~872.
    [84] [美]Von Golf-Racht著.陈钟祥译.石油科学进展.裂缝油藏工程基础[M].北京:石油工业出版社,1986.
    [85] TWISS RJ, MOORES E M. Structural geology[M]. New York: Freeman & Co, 1992.
    [86] Nelson R A. Geologic analysis of naturally fractured reservoirs[M]. Houston: Gulf Publishing Company, 1985.
    [87] S.H.Lee, C.L.Jensen, M.F.Lough. An efficient finite difference model for flow in a reservoir with multiple length-scale fractures. the 1999 SPE Annual Technical Conference and Exhibition[C]. Houston, Texas, U.S.A.: SPE 56752, 2004.
    [88] Nanqun He, Seong H.Lee, Clair L.Jensen. Combination of analytical, numerical and geostatistical methods to model naturally fractured reservoirs. the SPE Western Regional Meeting[C]. Bakersfield, CA: SPE 68832, 2001.
    [89] Arnaud Lange, Ramey Basquet, Bernard Bourbiaux. Hydraulic characterization of faults and fractures using a dual medium discrete fracture network simulator. the 11th ADIPEC[C]. Abu Dhabi: SPE 88675, 2004.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700