Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.     

Modelling of groundwater flow in heterogeneous porous media by finite element method

The HySuf‐FEM code (Hydrodynamic of Subsurface Flow by Finite Element Method) is proposed in this article in order to estimate the spatial variability of the transmissivity values of the Berrechid aquifer (Morocco). The calibration of the model is based on the hydraulic head, hydraulic conductivity and recharge. Three numerical tests are used to validate the model and verify its convergence. The first test case consists in using the steady analytical solution of the Poisson equation. In the second, the model has been compared with the hydrogeological system which is characterized by an unconfined monolayer (isotropic layer) and computed by using PMWIN‐MODFLOW software. The third test case is based on the comparison between the results of HySuf‐FEM and the multiple cell balance method in the aquifer system with natural boundaries case. Good agreement between the Hydrodynamic of Subsurface Flow, the numerical tests and the spatial distribution of the thickening of the hydrogeological system is deduced from the analysis and the interpretations of hydrogeological wells. Copyright © 2011 John Wiley & Sons, Ltd.     

A new finite element technique for the solution of two-phase flow through porous media

A new upstream weighting finite element technique is developed for improved solution of the two-phase immiscible flow equations. Unlike the upstream weighting technique used by previous investigators, the new technique does not employ finite difference concepts to achieve the required upstream weighting of relative permeabilities or mobilities. Instead, upstream weighting is achieved by (1) representing the relative permeabilities or mobilities as continuous functions expressed in terms of the shape functions and nodal values (2) using asymmetric weighting functions to weight the spatial terms in the flow equations. These weighting functions are constructed such that they are dependent on the flow direction along each side of an element.In conjunction with the proposed technique, two solution schemes for treating the resulting set of non-linear algebraic equations are presented. These are the fully-implicit chord slope incremental solution scheme and the Newton-Raphson solution scheme. Both schemes allow the use of large time steps without being unstable.The proposed numerical technique is applied to two problems (1) the one-dimensional Buckley-Leverett problem (2) the two-dimensional five-spot well flow problem. Results indicate that this technique is superior to not only earlier finite element schemes but also five-point upstream finite difference formulae.     

Temperature-driven fluid flow in porous media using a mixed finite element method and a finite volume method

We present a numerical scheme for the computation of conservative fluid velocity, pressure and temperature fields in a porous medium. For the velocity and pressure we use the primal–dual mixed finite element method of Trujillo and Thomas while for the temperature we use a cell-centered finite volume method. The motivation for this choice of discretization is to compute accurate conservative quantities. Since the variant of the mixed finite element method we use is not commonly used, the numerical schemes are presented in detail. We sketch the computational details and present numerical experiments that justify the accuracy predicted by the theory.     

Adaptive finite element simulation of Stokes flow in porous media

The Stokes problem describes flow of an incompressible constant-viscosity fluid when the Reynolds number is small so that inertial and transient-time effects are negligible. The numerical solution of the Stokes problem requires special care, since classical finite element discretization schemes, such as piecewise linear interpolation for both the velocity and the pressure, fail to perform. Even when an appropriate scheme is adopted, the grid must be selected so that the error is as small as possible. Much of the challenge in solving Stokes problems is how to account for complex geometry and to capture important features such as flow separation. This paper applies adaptive mesh techniques, using a posteriori error estimates, in the finite element solution of the Stokes equations that model flow at pore scales. Different selected numerical test cases associated with various porous geometrics are presented and discussed to demonstrate the accuracy and efficiency of our methodology.     

Analysis of Rayleigh waves in saturated porous elastic media by finite element method

A numerical procedure for the analysis of Rayleigh waves in saturated porous elastic media is proposed by use of the finite element method. The layer stiffness matrix, the layer mass matrix and the layer damping matrix in a layered system are presented for the discretized form of the solid-fluid equilibrium equation proposed by Biot. In order to consider the influence of the permeability coefficient on the behavior of Rayleigh waves, attention is focused on the following states: ‘drained’ state, ‘undrained’ state and the states between two extremes of ‘drained’ and ‘undrained’ states. It is found from computed results that the permeability coefficient exerts a significant effect on dispersion curves and displacement distributions of Rayleigh waves in saturated porous media.     

On the solution of the time-dependent convection-diffusion equation by the finite element method

A Petrov-Galerkin finite element method is presented for the time-dependent convection-diffusion equation. The scheme is based on bilinear time-space trial and quadratic in time-linear in space test functions, the latter being nonconforming. Second order in time and third order in space accuracy is obtained, and the schemes are free of numerical diffusion and disperson effects. Numerical results are presented which show excellent approximation properties.     

A mixed-dimensional finite volume method for two-phase flow in fractured porous media

We present a vertex-centered finite volume method for the fully coupled, fully implicit discretization of two-phase flow in fractured porous media. Fractures are discretely modeled as lower dimensional elements. The method works on unstructured, locally refined grids and on parallel computers with distributed memory. An implicit time discretization is employed and the nonlinear systems of equations are solved with a parallel Newton-multigrid method. Results from two-dimensional and three-dimensional simulations are presented.     

Lanczos method for the solution of groundwater flow in discretely fractured porous media

One of the more advanced approaches for simulating groundwater flow in fractured porous media is the discrete-fracture approach. This approach is limited by the large computational overheads associated with traditional modeling methods. In this work, we apply the Lanczos reduction method to the modeling of groundwater flow in fractured porous media using the discrete-fracture approach. The Lanczos reduction method reduces a finite element equation system to a much smaller tridiagonal system of first-order differential equations. The reduced system can be solved by a standard tridiagonal algorithm with little computational effort. Because solving the reduced system is more efficient compared to solving the original system, the simulation of groundwater flow in discretely fractured media using the reduction method is very efficient. The proposed method is especially suitable for the problem of large-scale and long-term simulation. In this paper, we develop an iterative version of Lanczos algorithm, in which the preconditioned conjugate gradient solver based on ORTHOMIN acceleration is employed within the Lanczos reduction process. Additional efficiency for the Lanczos method is achieved by applying an eigenvalue shift technique. The “shift” method can improve the Lanczos system convergence, by requiring fewer modes to achieve the same level of accuracy over the unshifted case. The developed model is verified by comparison with dual-porosity approach. The efficiency and accuracy of the method are demonstrated on a field-scale problem and compared to the performance of classic time marching method using an iterative solver on the original system. In spite of the advances, more theoretical work needs to be carried out to determine the optimal value of the shift before computations are actually carried out.     

Coupling multiscale finite element method for consolidation analysis of heterogeneous saturated porous media

The multiscale finite element method is developed for solving the coupling problems of consolidation of heterogeneous saturated porous media under external loading conditions. Two sets of multiscale base functions are constructed, respectively, for the pressure field of fluid flow and the displacement field of solid skeleton. The coupling problems are then solved with a multiscale numerical procedure in space and time domain. The heterogeneities induced by permeabilities and mechanical parameters of the saturated porous media are both taken into account. Numerical experiments are carried out for different cases in comparison with the standard finite element method. The numerical results show that the coupling multiscale finite element method can be successfully used for solving the complicated coupling problems. It reduces greatly the computing effort in both memory and time for transient problems.     

Coupling upscaling finite element method for consolidation analysis of heterogeneous saturated porous media

The coupling upscaling finite element method is developed for solving the coupling problems of deformation and consolidation of heterogeneous saturated porous media under external loading conditions. The method couples two kinds of fully developed methodologies together, i.e., the numerical techniques developed for calculating the apparent and effective physical properties of the heterogeneous media and the upscaling techniques developed for simulating the fluid flow and mass transport properties in heterogeneous porous media. Equivalent permeability tensors and equivalent elastic modulus tensors are calculated for every coarse grid block in the coarse-scale model of the heterogeneous saturated porous media. Moreover, an oversampling technique is introduced to improve the calculation accuracy of the equivalent elastic modulus tensors. A numerical integration process is performed over the fine mesh within every coarse grid element to capture the small scale information induced by non-uniform scalar field properties such as density, compressibility, etc. Numerical experiments are carried out to examine the accuracy of the developed method. It shows that the numerical results obtained by the coupling upscaling finite element method on the coarse-scale models fit fairly well with the reference solutions obtained by traditional finite element method on the fine-scale models. Moreover, this method gets more accurate coarse-scale results than the previously developed coupling multiscale finite element method for solving this kind of coupling problems though it cannot recover the fine-scale solutions. At the same time, the method developed reduces dramatically the computing effort in both CPU time and memory for solving the transient problems, and therefore more large and computational-demanding coupling problems can be solved by computers.     

Similarity solution of axisymmetric flow in porous media

Applications of the axisymmetric Boussinesq equation to groundwater hydrology and reservoir engineering have long been recognised. An archetypal example is invasion by drilling fluid into a permeable bed where there is initially no such fluid present, a circumstance of some importance in the oil industry. It is well known that the governing Boussinesq model can be reduced to a nonlinear ordinary differential equation using a similarity variable, a transformation that is valid for a certain time-dependent flux at the origin. Here, a new analytical approximation is obtained for this case. The new solution,, which has a simple form, is demonstrated to be highly accurate.     

流体饱和多孔隙介质二维弹性波方程正演模拟的小波有限元法

本文将小波有限元法引入到流体饱和多孔隙介质二维波动方程的正演模拟中,以二维Daubechies小波的尺度函数代替多项式函数作为插值函数,构造二维张量积小波单元.引入一类特征函数解决了Daubechies小波没有显式解析表达式所带来的基函数积分值计算问题,并推导出计算分数节点上Daubechies小波函数值的递推公式,从而构造出由小波系数空间到波场位移空间的快速小波变换.数值模拟结果表明该方法是有效的.     

Computer modeling of geoelectromagnetic fields in three-dimensional media by the finite element method

This work describes the basic approaches to the solution of three-dimensional (3D) problems of geoelectromagnetism with the use of the finite element method and the possibilities of the GeoEM program complex for its implementation. The methods of modeling geoelectromagnetic fields for the most widely known types of controllable sources and the mathematical apparatus for the solution of problems of magnetotelluric soundings are considered. Examples of the calculations of 3D fields for the horizontal and vertical electric lines in planning electric exploration works on the shelf are presented, as well as an example of the 3D interpretation of array transient electromagnetic sounding data during the search for deep-seated target objects in conditions of a mostly heterogeneous upper part of the section.     

On the applicability of the Biot solution to the problems of elastic propagation in saturated porous media

In this paper, the solution of the system of homogeneous Biot equations, which was derived by Biot for the displacement vectors of plane monochrome elastic waves propagating in a homogeneous infinite two-phase medium, is expanded to the case where the propagation area of the elastic waves is limited and the wavefront is a piecewise smooth curved surface. It is shown that the arbitrary system of homogeneous Biot equations for the displacement vectors of the solid and liquid phases can be reduced to three different equations pertaining to the class of Helmholtz equations. From this, irrespective of the geometry of the seismic wavefront and the boundaries of the studied two-phase medium, there is the following. (1) Each displacement vector (of the solid and liquid phase) splits into three independent vectors satisfying three different Helmholtz equations. Two of these vectors correspond to the two types of compressional waves, namely, fast waves (waves of the first kind) and slow waves (waves of the second kind). The third vector describes shear waves. (2) The similar (related to the same wave type) components of the displacement vector in the solid and liquid phases satisfy the same Helmholtz equation and are linked with each other through a corresponding scalar factor that is expressed in terms of the coefficients of the Biot equations. Taking into account the established properties of the displacement vectors in the solid and liquid phases seems to be helpful in the problems dealing with calculation of elastic fields of arbitrary sources in piecewise-homogeneous two-phase media.     

Directional interpolation infinite element for dynamic problems in saturated porous media

In this study, a directional interpolation infinite element suited to a saturated porous medium is presented to account for dynamic problems with semi-infi     

Free-surface flow in porous media and periodic solution of the shallow-flow approximation

The equations describing the flow of liquid in a porous medium with a free surface are expanded when the shallow-flow assumption holds. Second-order theory is used to describe the propagation of steady periodic motion in the medium, driven by the oscillating level of a reservoir in contact with it. A linearized solution of the second-order theory is compared with a numerical solution and is found adequate even when the amplitude of the motion is comparable to the mean depth of the liquid. The predictions of the analysis are found in good agreement with two laboratory experiments.     

An Hermitian finite element solution of the two-dimensional saturated-unsaturated flow equation

This paper describes a Galerkin-type finite element solution of the two-dimensional saturated-unsaturated flow equation. The numerical solution uses an incomplete (reduced) set of Hermitian cubic basis functions and is formulated in terms of normal and tangential coordinates. The formulation leads to continuous pressure gradients across interelement boundaries for a number of well-defined element configurations, such as for rectangular and circular elements. Other elements generally lead to discontinuous gradients; however, the gradients remain uniquely defined at the nodes. The method avoids calculation of second-order derivatives, yet retains many of the advantages associated with Hermitian elements. A nine-point Lobatto-type integration scheme is used to evaluate all local element integrals. This alternative scheme produces about the same accuracy as the usual 9- or 16-point Gaussian quadrature schemes, but is computationally more efficient.     

A new benchmark semi-analytical solution for density-driven flow in porous media

A new benchmark semi-analytical solution is proposed for the verification of density-driven flow codes. The problem deals with a synthetic square porous cavity subject to different salt concentrations at its vertical walls. A steady state semi-analytical solution is investigated using the Fourier–Galerkin method. Contrarily to the standard Henry problem, the cavity benchmark allows high truncation orders in the Fourier series and provides semi-analytical solutions for very small diffusion cases. The problem is also investigated numerically to validate the semi-analytical solution. The obtained results represent a set of new test case high quality data that can be effectively used for benchmarking density-driven flow codes.