Investigating the effects of intersection flow localization in equivalent-continuum-based upscaling of flow in discrete fracture networks

Author:

Kottwitz Maximilian O.ORCID,Popov Anton A.,Abe Steffen,Kaus Boris J. P.ORCID

Abstract

Abstract. Predicting effective permeabilities of fractured rock masses is a crucial component of reservoir modeling. Its often realized with the discrete fracture network (DFN) method, whereby single-phase incompressible fluid flow is modeled in discrete representations of individual fractures in a network. Depending on the overall number of fractures, this can result in high computational costs. Equivalent continuum models (ECMs) provide an alternative approach by subdividing the fracture network into a grid of continuous medium cells, over which hydraulic properties are averaged for fluid flow simulations. While continuum methods have the advantage of lower computational costs and the possibility of including matrix properties, choosing the right cell size to discretize the fracture network into an ECM is crucial to provide accurate flow results and conserve anisotropic flow properties. Whereas several techniques exist to map a fracture network onto a grid of continuum cells, the complexity related to flow in fracture intersections is often ignored. Here, numerical simulations of Stokes flow in simple fracture intersections are utilized to analyze their effect on permeability. It is demonstrated that intersection lineaments oriented parallel to the principal direction of flow increase permeability in a process we term intersection flow localization (IFL). We propose a new method to generate ECMs that includes this effect with a directional pipe flow parameterization: the fracture-and-pipe model. Our approach is compared against an ECM method that does not take IFL into account by performing ECM-based upscaling with a massively parallelized Darcy flow solver capable of representing permeability anisotropy for individual grid cells. While IFL results in an increase in permeability at the local scale of the ECM cell (fracture scale), its effects on network-scale flow are minor. We investigated the effects of IFL for test cases with orthogonal fracture formations for various scales, fracture lengths, hydraulic apertures, and fracture densities. Only for global fracture porosities above 30 % does IFL start to increase the systems permeability. For lower fracture densities, the effects of IFL are smeared out in the upscaling process. However, we noticed a strong dependency of ECM-based upscaling on its grid resolution. Resolution tests suggests that, as long as the cell size is smaller than the minimal fracture length and larger than the maximal hydraulic aperture of the considered fracture network, the resulting effective permeabilities and anisotropies are resolution-independent. Within that range, ECMs are applicable to upscale flow in fracture networks.

Funder

Bundesministerium für Bildung und Forschung

Publisher

Copernicus GmbH

Subject

Paleontology,Stratigraphy,Earth-Surface Processes,Geochemistry and Petrology,Geology,Geophysics,Soil Science

Reference82 articles.

1. Alghalandis, Y. F.: ADFNE: Open source software for discrete fracture network engineering, two and three dimensional applications, Comput. Geosci., 102, 1–11, https://doi.org/10.1016/j.cageo.2017.02.002, 2017. a, b, c, d

2. Andrä, H., Combaret, N., Dvorkin, J., Glatt, E., Han, J., Kabel, M., Keehm, Y., Krzikalla, F., Lee, M., Madonna, C., Marsh, M., Mukerji, T., Saenger, E. H., Sain, R., Saxena, N., Ricker, S., Wiegmann, A., and Zhan, X.: Digital rock physics benchmarks-Part I: Imaging and segmentation, Comput. Geosci., 50, 25–32, https://doi.org/10.1016/j.cageo.2012.09.005, 2013a. a

3. Andrä, H., Combaret, N., Dvorkin, J., Glatt, E., Han, J., Kabel, M., Keehm, Y., Krzikalla, F., Lee, M., Madonna, C., Marsh, M., Mukerji, T., Saenger, E. H., Sain, R., Saxena, N., Ricker, S., Wiegmann, A., and Zhan, X.: Digital rock physics benchmarks-part II: Computing effective properties, Comput. Geosci., 50, 33–43, https://doi.org/10.1016/j.cageo.2012.09.008, 2013b. a

4. Bai, T., Maerten, L., Gross, M. R., and Aydin, A.: Orthogonal cross joints: Do they imply a regional stress rotation?, J. Struct. Geol., 24, 77–88, https://doi.org/10.1016/S0191-8141(01)00050-5, 2002. a

5. Balay, S., Abhyankar, S., Adams, M., Brown, J., Brune, P., Buschelman, K., Dalcin, L., Dener, A., Eijkhout, V., Gropp, W., Karpeyev, D., Kaushik, D., Knepley, M., May, D., Curfman McInnes, L., Mills, R., Munson, T., Rupp, K., Sanan, P., Smith, B., Zampini, S., Zhang, H., and Zhang, H.: PETSc Users Manual: Revision 3.10, Tech. rep., Argonne National Lab.(ANL), Argonne, IL (US), 2018. a, b

Cited by 10 articles. 订阅此论文施引文献 订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3