A Finite-Volume Method with Hexahedral Multiblock Grids for Modeling Flow in Porous Media

This paper presents a finite-volume method for hexahedral multiblock grids to calculate multiphase flow in geologically complex reservoirs. Accommodating complex geologic and geometric features in a reservoir model (e.g., faults) entails non-orthogonal and/or unstructured grids in place of conventional (globally structured) Cartesian grids. To obtain flexibility in gridding as well as efficient flow computation, we use hexahedral multiblock grids. These grids are locally structured, but globally unstructured. One major advantage of these grids over fully unstructured tetrahedral grids is that most numerical methods developed for structured grids can be directly used for dealing with the local problems. We present several challenging examples, generated via a commercially available tool, that demonstrate the capabilities of hexahedral multiblock gridding. Grid quality is discussed in terms of uniformity and orthogonality. The presence of non-orthogonal grid and full permeability tensors requires the use of multi-point discretization methods. A flux-continuous finite-difference (FCFD) scheme, previously developed for stratigraphic hexahedral grid with full-tensor permeability, is employed for numerical flow computation. We extend the FCFD scheme to handle exceptional configurations (i.e. three- or five-cell connections as opposed to the regular four), which result from employing multiblock gridding of certain complex objects. In order to perform flow simulation efficiently, we employ a two-level preconditioner for solving the linear equations that results from the wide stencil of the FCFD scheme. The individual block, composed of cells that form a structured grid, serves as the local level; the higher level operates on the global block configuration (i.e. unstructured component). The implementation uses an efficient data structure where each block is wrapped with a layer of neighboring cells. We also examine splitting techniques [14] for the linear systems associated with the wide stencils of our FCFD operator. We present three numerical examples that demonstrate the method: (1) a pinchout, (2) a faulted reservoir model with internal surfaces and (3) a real reservoir model with multiple faults and internal surfaces.     

Control-Volume Discretization Method for Quadrilateral Grids with Faults and Local Refinements

This paper extends the multipoint flux-approximation (MPFA) control-volume method to quadrilateral grids for which the adjacent cells do not necessarily share corners. Examples are grids with faults and locally refined grids. This paper gives a derivation of the method for such grids. The difference between two-point flux-approximation (TPFA) results and MPFA results for faults and local grid refinements is demonstrated for synthetic problems. Further, the results are compared with results from uniform fine-grid simulations. The effect of repeated fault patterns as well as anisotropy is investigated. Large errors may be found for the TPFA method for flow through a series of faults in an anisotropic medium. Finally, a comparison is done for a reservoir field application.     

Modeling branched and intersecting faults in reservoir-geomechanics models with the extended finite element method

Faults are geological entities with thickness several orders of magnitude smaller than the grid blocks typically used to discretize geological formations. On using the extended finite element method (X-FEM), a structured mesh suffices and the faults can arbitrarily cut the elements in the mesh. Modeling branched and intersecting faults is a challenge, in particular when the faults work as internal fluid flow conduits that allow fluid flow in the faults as well as to enter/leave the faults. By appropriately selecting the enrichment function and the nodes to be enriched, we are able to capture the special characteristics of the solution in the vicinity of the fault. We compare different enrichment schemes for strong discontinuities and develop new continuous enrichment functions with discontinuous derivatives to model branched and intersecting weak discontinuities. Symmetric fluid flows within the regions embedded by branched, coplanar intersecting, and noncoplanar intersecting faults are considered to verify the effectiveness of the proposed enrichment scheme. For a complex fault consisting of branched and intersecting faults, the accuracy and efficiency of the X-FEM is compared with the FEM. We also demonstrate different slipping scenarios for branched and intersecting faults with the same friction coefficient. In addition, fault slipping triggered by an injection pressure and three-dimensional fluid flows are modeled to show the versatility of the proposed enrichment scheme for branched and intersecting weak discontinuities.     

G23FM: a tool for meshing complex geological media

We present G23FM, a mesh generation tool for discretizing two- and three-dimensional complex fractured geological media. G23FM includes different techniques to generate finite element grids that maintain the geometric integrity of input surfaces, and geologic data and produce optimal triangular/tetrahedral grids for flow and transport simulations. G23FM generates grid for two-dimensional cross-sections, represents faults and fractures, for three-dimensional fractured media, and has the capability of including finer grids. Different examples are presented to illustrate some of the main features of G23FM.     

Model reduction and discretization using hybrid finite volumes for flow in porous media containing faults

In this paper, we study two different model reduction strategies for solving problems involving single phase flow in a porous medium containing faults or fractures whose location and properties are known. These faults are represented as interfaces of dimension N ? 1 immersed in an N dimensional domain. Both approaches can handle various configurations of position and permeability of the faults, and one can handle different fracture permeabilities on the two inner sides of the fracture. For the numerical discretization, we use the hybrid finite volume scheme as it is known to be well suited to simulating subsurface flow. Some results, which may be of use in the implementation of the proposed methods in industrial codes, are demonstrated.     

Development of a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method

In the numerical modeling of fluid flow in heterogeneous geological media, large material contrasts associated with complexly intersected material interfaces are challenging, not only related to mesh discretization but also for the accurate realization of the corresponding boundary constraints. To address these challenges, we developed a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method (NMM) and the Lagrange multiplier method (LMM) for modeling boundary constraints. The advantages of NMM include meshing efficiency with fixed mathematical grids (covers), the convenience of increasing the approximation precision, and the high integration precision provided by simplex integration. In this discontinuous approach, the elements intersected by material interfaces are divided into different elements and linked together using the LMM. We derive and compare different forms of LMMs and arrive at a new LMM that is efficient in terms of not requiring additional Lagrange multiplier topology, yet stringently derived by physical principles, and accurate in numerical performance. To demonstrate the accuracy and efficiency of the NMM with the developed LMM for boundary constraints, we simulate a number of verification and demonstration examples, involving a Dirichlet boundary condition and dense and intersected material interfaces. Last, we applied the developed model for modeling fluid flow in heterogeneous media with several material zones containing a fault and an opening. We show that the developed discontinuous approach is very suitable for modeling fluid flow in strongly heterogeneous media with good accuracy for large material contrasts, complex Dirichlet boundary conditions, or complexly intersected material interfaces. Copyright © 2015 John Wiley & Sons, Ltd.     

基于地质剖面构建三维地质模型的方法研究

地质剖面是三维地质建模的重要数据源,运用地质剖面构建三维地质模型的方法应用较为广泛。在模型构建之前,统一确定模数据的坐标系和比例足,建立原始资料数据库。构建三维地质模型的关键是不同地质界面。本文详细介绍了模型边界面、断层面、地层界面、岩体界面等4种主要地质界面的构建流程与方法,尤其对褶皱构造、地层界面的断层效应、复杂岩体界面等的构建进行了重点阐述按模型边界面(模型的底界面和四周边界面)、DEM面、断层面、其他地质界面的顺序依次构建地质界面,构建断层面和其他地质界面时严格按先新后老的顺序。运用已构建好的地质界面按先新后老的顺序逐个、依次建立单个地质体,再将所有地质体的面模型组合成整个模型的面模型。通过对面模型进行网格(实体)填充和对网格赋予相应的属性值,最终构建三维地质模型。     

Open-source MATLAB implementation of consistent discretisations on complex grids

Accurate geological modelling of features such as faults, fractures or erosion requires grids that are flexible with respect to geometry. Such grids generally contain polyhedral cells and complex grid-cell connectivities. The grid representation for polyhedral grids in turn affects the efficient implementation of numerical methods for subsurface flow simulations. It is well known that conventional two-point flux-approximation methods are only consistent for K-orthogonal grids and will, therefore, not converge in the general case. In recent years, there has been significant research into consistent and convergent methods, including mixed, multipoint and mimetic discretisation methods. Likewise, the so-called multiscale methods based upon hierarchically coarsened grids have received a lot of attention. The paper does not propose novel mathematical methods but instead presents an open-source Matlab? toolkit that can be used as an efficient test platform for (new) discretisation and solution methods in reservoir simulation. The aim of the toolkit is to support reproducible research and simplify the development, verification and validation and testing and comparison of new discretisation and solution methods on general unstructured grids, including in particular corner point and 2.5D PEBI grids. The toolkit consists of a set of data structures and routines for creating, manipulating and visualising petrophysical data, fluid models and (unstructured) grids, including support for industry standard input formats, as well as routines for computing single and multiphase (incompressible) flow. We review key features of the toolkit and discuss a generic mimetic formulation that includes many known discretisation methods, including both the standard two-point method as well as consistent and convergent multipoint and mimetic methods. Apart from the core routines and data structures, the toolkit contains add-on modules that implement more advanced solvers and functionality. Herein, we show examples of multiscale methods and adjoint methods for use in optimisation of rates and placement of wells.     

Multiscale mixed/mimetic methods on corner-point grids

Multiscale simulation is a promising approach to facilitate direct simulation of large and complex grid models for highly heterogeneous petroleum reservoirs. Unlike traditional simulation, approaches based on upscaling/downscaling, multiscale methods seek to solve the full flow problem by incorporating subscale heterogeneities into local discrete approximation spaces. We consider a multiscale formulation based on a hierarchical grid approach, where basis functions with subgrid resolution are computed numerically to correctly and accurately account for subscale variations from an underlying (fine-scale) geomodel when solving the global flow equations on a coarse grid. By using multiscale basis functions to discretise the global flow equations on a (moderately sized) coarse grid, one can retain the efficiency of an upscaling method and, at the same time, produce detailed and conservative velocity fields on the underlying fine grid. For pressure equations, the multiscale mixed finite-element method (MsMFEM) has been shown to be a particularly versatile approach. In this paper, we extend the method to corner-point grids, which is the industry standard for modelling complex reservoir geology. To implement MsMFEM, one needs a discretisation method for solving local flow problems on the underlying fine grids. In principle, any stable and conservative method can be used. Here, we use a mimetic discretisation, which is a generalisation of mixed finite elements that gives a discrete inner product, allows for polyhedral elements, and can (easily) be extended to curved grid faces. The coarse grid can, in principle, be any partition of the subgrid, where each coarse block is a connected collection of subgrid cells. However, we argue that, when generating coarse grids, one should follow certain simple guidelines to achieve improved accuracy. We discuss partitioning in both index space and physical space and suggest simple processing techniques. The versatility and accuracy of the new multiscale mixed methodology is demonstrated on two corner-point models: a small Y-shaped sector model and a complex model of a layered sedimentary bed. A variety of coarse grids, both violating and obeying the above mentioned guidelines, are employed. The MsMFEM solutions are compared with a reference solution obtained by direct simulation on the subgrid.     

小尺度断层定量预测及对注水开发的影响——以渤中34油田为例

大尺度断层往往控制了沉积盆地的形成和油气成藏,而小尺度断层则影响着注水开发效果和剩余油分布.大尺度断层可以通过二维或三维地震资料识别,而小尺度断层的识别则特别困难.本文提出了一种基于断层分形生长模式和三维地质力学模拟相结合来确定小尺度断层的数量、发育位置和方位的方法,并根据油田开发动态资料来确定小尺度断层对注水开发和剩余油分布的影响.将地震上识别的大尺度断层引入到三维数值力学模型中,模拟大尺度断层形成时期断裂带附近的应力扰动作用,然后结合破裂准则来建立最易发生破裂的方位和最大库伦剪切应力网格,以这两套网格和断层尺度的幂律分布确定的小尺度断层数量为约束条件来确定随机模型,对小尺度断层的密度、产状和发育位置进行定量预测.研究表明:利用分形理论和三维地质力学模拟可以对大尺度断层伴生小尺度断层进行有效预测;小尺度断层对注水开发效果和剩余油分布的影响取决于小尺度断层的规模(断距)以及小尺度断层方位和注采方向的关系.     

Mortar Upscaling for Multiphase Flow in Porous Media

In mortar space upscaling methods, a reservoir is decomposed into a series of subdomains (blocks) in which independently constructed numerical grids and possibly different physical models and discretization techniques can be employed in each block. Physically meaningful matching conditions are imposed on block interfaces in a numerically stable and accurate way using mortar finite element spaces. Coarse mortar grids and fine subdomain grids provide two-scale approximations. In the resulting effective solution flow is computed in subdomains on the fine scale while fluxes are matched on the coarse scale. In addition the flexibility to vary adaptively the number of interface degrees of freedom leads to more accurate multiscale approximations. This methodology has been implemented in the Center for Subsurface Modeling's multiphysics multiblock simulator IPARS (Integrated Parallel Accurate reservoir Simulator). Computational experiments demonstrate that this approach is scalable in parallel and it can be applied to non-matching grids across the interface, multinumerics and multiphysics models, and mortar adaptivity. Moreover unlike most upscaling approaches the underlying systems can be treated fully implicitly.     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Enhanced Velocity Mixed Finite Element Methods for Flow in Multiblock Domains

The paper presents a new approach to discretizing flow in porous media via mixed finite element methods on non-matching multiblock grids. The velocity space along the interfaces is enhanced to give flux-continuous approximation. No additional matching conditions need to be imposed. The computational complexity of the resulting algebraic problem is comparable to the one for the single-block case. A priori error estimates for the pressure and the velocity and numerical experiments confirming the theory are presented.     

A phase-field method for the direct simulation of two-phase flows in pore-scale media using a non-equilibrium wetting boundary condition

Advances in pore-scale imaging (e.g., μ-CT scanning), increasing availability of computational resources, and recent developments in numerical algorithms have started rendering direct pore-scale numerical simulations of multi-phase flow on pore structures feasible. Quasi-static methods, where the viscous and the capillary limit are iterated sequentially, fall short in rigorously capturing crucial flow phenomena at the pore scale. Direct simulation techniques are needed that account for the full coupling between capillary and viscous flow phenomena. Consequently, there is a strong demand for robust and effective numerical methods that can deliver high-accuracy, high-resolution solutions of pore-scale flow in a computationally efficient manner. Direct simulations of pore-scale flow on imaged volumes can yield important insights about physical phenomena taking place during multi-phase, multi-component displacements. Such simulations can be utilized for optimizing various enhanced oil recovery (EOR) schemes and permit the computation of effective properties for Darcy-scale multi-phase flows.We implement a phase-field model for the direct pore-scale simulation of incompressible flow of two immiscible fluids. The model naturally lends itself to the transport of fluids with large density and viscosity ratios. In the phase-field approach, the fluid-phase interfaces are expressed in terms of thin transition regions, the so-called diffuse interfaces, for increased computational efficiency. The conservation law of mass for binary mixtures leads to the advective Cahn–Hilliard equation and the condition that the velocity field is divergence free. Momentum balance, on the other hand, leads to the Navier–Stokes equations for Newtonian fluids modified for two-phase flow and coupled to the advective Cahn–Hilliard equation. Unlike the volume of fluid (VoF) and level-set methods, which rely on regularization techniques to describe the phase interfaces, the phase-field method facilitates a thermodynamic treatment of the phase interfaces, rendering it more physically consistent for the direct simulations of two-phase pore-scale flow. A novel geometric wetting (wall) boundary condition is implemented as part of the phase-field method for the simulation of two-fluid flows with moving contact lines. The geometric boundary condition accurately replicates the prescribed equilibrium contact angle and is extended to account for dynamic (non-equilibrium) effects. The coupled advective Cahn–Hilliard and modified Navier–Stokes (phase-field) system is solved by using a robust and accurate semi-implicit finite volume method. An extension of the momentum balance equations is also implemented for Herschel–Bulkley (non-Newtonian) fluids. Non-equilibrium-induced two-phase flow problems and dynamic two-phase flows in simple two-dimensional (2-D) and three-dimensional (3-D) geometries are investigated to validate the model and its numerical implementation. Quantitative comparisons are made for cases with analytical solutions. Two-phase flow in an idealized 2-D pore-scale conduit is simulated to demonstrate the viability of the proposed direct numerical simulation approach.     

Modeling fractures as interfaces: a model for Forchheimer fractures

In this paper we are concerned with modeling single phase flow in a medium with known fractures. In particular we are interested in the case in which the flow rate in the fractures is large enough to make it appropriate to use Forchheimer’s law for modeling the flow in the fractures even though the flow in the surrounding domain is such that Darcy’s law is adequate. We describe a model in which the fractures are treated as interfaces. We also consider the case of intersecting fractures and the case of nonconforming meshes.     

基于水平集的复杂三维地层模型建模及在地下工程中的应用研究

由于地质体存在断层、尖灭、出露等复杂地质现象,在三维地层建模时,为了表达这些现象,无论是面模型还是体模型或混合模型,都存在空间分割或曲面求交线的问题。由于地质体拓扑关系的复杂性、数据误差以及计算机精度问题,使得这些模型在实际建模过程中常常失效。运用水平集理论可以有效解决这一问题,水平集用隐函数表达曲面（或超曲面）,可以实现复杂地质体的表达及并、交、差等拓扑运算。在三维地层建模中,插值生成各种地质界面后,用水平集表达这些地质界面,利用水平理论完成各种复杂的拓扑操作,建立以水平集表达的三维地层模型。在此基础上,插入水平集表达的各种工程活动界面,利用水平集理论进行拓扑操作,可构建各种工程活动后的地层模型。再利用Marching cube算法抽取各种地质界面或工程活动界面,构建可用于实时可视化或用于工程评估（如有限元计算）等的NMTINF-BR地层模型或工程活动后的NMTINF-BR地层模型。     

Mortar coupling and upscaling of pore-scale models

Pore-scale models are becoming increasingly useful as predictive tools for modeling flow and transport in porous media. These models can accurately represent the 3D pore-structure of real media. Currently first-principles modeling methods are being employed for obtaining qualitative and quantitative behavior. Generally, artificial, simple boundary conditions are imposed on a model that is used as a stand-alone tool for extracting macroscopic parameters. However, realistic boundary conditions, reflecting flow and transport in surrounding media, may be necessary for behavior that occurs over larger length scales or including pore-scale models in a multiscale setting. Here, pore-scale network models are coupled to adjacent media (additional pore-scale or continuum-scale models) using mortars. Mortars are 2D finite-element spaces employed to couple independent subdomains by enforcing continuity of pressure and flux at shared boundary interfaces. While mortars have been used in the past to couple subdomains of different models, physics, and meshes, they are extended here for the first time to pore-scale models. The approach is demonstrated by modeling single-phase flow in coupled pore-scale models, but the methodology can be utilized to model dynamic processes and perform multiscale modeling in 3D continuum simulators for flow and transport.     

Control‐volume mixed finite element methods

A key ingredient in simulation of flow in porous media is accurate determination of the velocities that drive the flow. Large‐scale irregularities of the geology (faults, fractures, and layers) suggest the use of irregular grids in simulation. This paper presents a control‐volume mixed finite element method that provides a simple, systematic, easily implemented procedure for obtaining accurate velocity approximations on irregular (i.e., distorted logically rectangular) block‐centered quadrilateral grids. The control‐volume formulation of Darcy’s law can be viewed as a discretization into element‐sized “tanks” with imposed pressures at the ends, giving a local discrete Darcy law analogous to the block‐by‐block conservation in the usual mixed discretization of the mass‐conservation equation. Numerical results in two dimensions show second‐order convergence in the velocity, even with discontinuous anisotropic permeability on an irregular grid. The method extends readily to three dimensions.     

Upscaled modeling of well singularity for simulating flow in heterogeneous formations

Subsurface flows are affected by geological variability over a range of length scales. The modeling of well singularity in heterogeneous formations is important for simulating flow in aquifers and petroleum reservoirs. In this paper, two approaches in calculating the upscaled well index to capture the effects of fine scale heterogeneity in near-well regions are presented and applied. We first develop a flow-based near-well upscaling procedure for geometrically flexible grids. This approach entails solving local well-driven flows and requires the treatment of geometric effects due to the nonalignment between fine and coarse scale grids. An approximate coarse scale well model based on a well singularity analysis is also proposed. This model, referred to as near-well arithmetic averaging, uses only the fine scale permeabilities at well locations to compute the coarse scale well index; it does not require solving any flow problems. These two methods are systematically tested on three-dimensional models with a variety of permeability distributions. It is shown that both approaches provide considerable improvement over a simple (arithmetic) averaging approach to compute the coarse scale well index. The flow-based approach shows close agreement to the fine scale reference model, and the near-well arithmetic averaging also offers accuracy for an appropriate range of parameters. The interaction between global flow and near-well upscaling is also investigated through the use of global fine scale solutions in near-well scale-up calculations.     

青藏公路及铁路沿线的活动构造与其次生灾害

基于对青藏公路及铁路沿线区域约63400km2范围内1∶10万地质环境遥感的调查研究结果,采用以ETM为主、重点地区以IKONOS为辅的遥感数据、遥感解译与地面验证相结合的技术方法在研究区解译出数百条断层。根据这些断层的延伸方向和分布位置可归纳为20条断裂带(RF),它们与以往以地面调查为主的断裂带(GF)基本吻合。研究区的断裂带基本上均为活动断裂带,且都存在发震断层,将已发生地震的震级和频率较高(Ms≥5至>8)的断裂带定为主要地震构造带。青藏线共有5个主要地震构造带,其中东昆仑地震构造带及当雄—羊八井地震构造带是研究区内可能对铁路工程建设及公路、铁路运行安全带来严重影响的最重要的地震构造带。滑坡、泥石流是青藏线活动构造的主要次生灾害。该区的活动断裂对滑坡分布有一定的影响,但并没有控制作用;地震构造不但控制了泥石流的分布,且地震活动对泥石流活动的触发作用也非常明显。