A computational method for approximating a Darcy–Stokes system governing a vuggy porous medium

We develop and analyze a mixed finite element method for the solution of an elliptic system modeling a porous medium with large cavities, called vugs. It consists of a second-order elliptic (i.e., Darcy) equation on part of the domain coupled to a Stokes equation on the rest of the domain, and a slip boundary condition (due to Beavers–Joseph–Saffman) on the interface between them. The tangential velocity is not continuous on the interface. We consider a 2-D vuggy porous medium with many small cavities throughout its extent, so the interface is not isolated. We use a certain conforming Stokes element on rectangles, slightly modified near the interface to account for the tangential discontinuity. This gives a mixed finite element method for the entire Darcy–Stokes system with a regular sparsity pattern that is easy to implement, independent of the vug geometry, as long as it aligns with the grid. We prove optimal global first-order L ² convergence of the velocity and pressure, as well as the velocity gradient in the Stokes domain. Numerical results verify these rates of convergence and even suggest somewhat better convergence in certain situations. Finally, we present a lower dimensional space that uses Raviart–Thomas elements in the Darcy domain and uses our new modified elements near the interface in transition to the Stokes elements.     

Comparison of Galerkin-type discretization techniques for two-phase flow in heterogeneous porous media

A comparison of Standard Galerkin, Petrov-Galerkin, and Fully-Upwind Galerkin methods for the simulation of two-phase flow in heterogeneous porous media is presented. On the basis of the coupled pressure-saturation equations, a generalized formulation for all three finite element methods is derived and analysed. For flow in homogeneous media, the Petrov-Galerkin method gives excellent results. But this method fails miserably for problems with heterogeneous media. This is because it is not able to capture correctly processes that take place at interfaces when, for instance, the capillary pressure-saturation relationship after Brooks and Corey is assumed. The Fully-Upwind Galerkin method is superior to the Petrov-Galerkin approach because it is able to give correct results for flow in homogeneous and heterogeneous media for the two models of van Genuchten and Brooks-Corey. The widely used formulation which is correct for the homogeneous case cannot be used for heterogeneous media. Instead the straightforward approach of gradp_c in combination with a chord-slope technique must be utilized.     

General absorbing boundary conditions for dynamic analysis of fluid-saturated porous media

General absorbing boundary conditions based on Biot's two-phase mixture theory and paraxial approximation is presented for the dynamic analysis of fluid-saturated porous media with isotropic, transverse isotropic, and anisotropic properties. For the last two cases, the equivalent Lame's constants, under conditions of uniqueness, are introduced to facilitate the analytical solutions. The numerical results show that the proposed absorbing boundary can greatly suppress spuriously reflected waves and efficiently model the far field of the system with sufficient accuracy.     

Adaptive hierarchical upscaling of flow in heterogeneous reservoirs based on an a posteriori error estimate

This paper treats the upscaling of the absolute permeability in a heterogeneous reservoir. By replacing the fine scale permeability tensor with an upscaled, or effective permeability tensor, a modelling error is introduced. An a posteriori error estimate on this modelling error is formulated and tested. An implementation of the theory, based on domain decomposition coupled with a hierarchical representation of the absolute permeability field, is given. As hierarchical basis functions we have chosen the Haar system, which leads to a wavelet representation of the permeability. The wavelet representation offers a natural upscaling technique which resembles the highcut filters commonly used in signal analysis. This procedure represents an adaptive upscaling method. The numerical results show that this method conserves both the dissipation and the mean velocity in the problem fairly well. The a posteriori error estimate on the modelling error coupled with domain decomposition methods constitutes a powerful modelling tool.     

Velocities of elastic waves andv P/v S ratios in dry and water-saturated rock samples at high pressure

Elastic wave velocities and the ratiov _P/v _S were studied for dry and initially saturated samples of carbonate and crystalline rocks at pressures to 2 kbar. In initially saturated samplesv _P increases in crystalline rock, whereas in sedimentary rock it may either increase or decrease with increasing pressure. Under the same conditionsv _S remains approximately constant in crystalline rocks and decreases in sedimentary samples. The ratiov _P/v _S as a function of pressure increases in dry rocks and decreases in saturated ones. Saturated samples always have higherv _P/v _S ratios than dry ones. It appears that the geometry of the pore space influences the acoustic properties of saturated crystalline rocks.     

横向非均匀介质中弹性地震波传播的数值模拟

应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法，将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质，给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度。     

Sediment transport patterns at Trafalgar offshore windfarm

A simple analytical model of wave propagation has been developed in order to study the potential sediment transport patterns due to the action of currents and waves in the neighborhood of cylindrical structures as well inside a group of these structures. The attention is focused on the study of Trafalgar offshore windfarm, a case in which it has been necessary to analyze the flow trough porous structures in order to model fish growing-cages planned to be installed at each aerogenerator structure. The results are obtained by averaging over one period of wave in order to evaluate the net potential sediment transport. The analysis of the results reveals how the processes of wave diffraction and reflection give rise to periodic patterns of sediment transport around and between the structures.     

雪宝顶流域钙华景观生态风险评价的研究探讨

景观生态风险评价能够有效地对区域景观格局进行优化并为风险防范提供理论基础。文献计量学分析表明,国内岩溶景观生态风险评价研究主要从“研究对象”、“评价模型”、“评价方法”和“评价指标”4个方面分析;而钙华作为特殊的岩溶体,目前对它的生态风险评价几乎未开展。文章从岩溶景观生态风险评价入手,基于模糊数学综合评判法和层次分析法,提出以地表钙华体与地表水循环作为景观主体,并对主体因素逐步细化,形成多层次综合评价模型,展开对雪宝顶流域钙华景观的生态风险评价。     

VTI介质中的弹性阻抗与参数提取

弹性阻抗是波阻抗的扩展,可以更好的反映储层和流体性质.常规弹性阻抗方程忽略了各向异性的影响,不能正确地描述各向异性介质中弹性阻抗的性质.本文基于Ruger给出的VTI介质反射系数近似公式,结合Connolly推导弹性阻抗的思想,推导出了一种VTI介质中的弹性阻抗方程.通过引入归一化常数,对方程进行了标准化,消除了方程量纲随角度变化而变化的不足.结合常规弹性阻抗反演的流程,得出了VTI介质中弹性阻抗反演的流程.最后,开展了从VTI介质弹性阻抗数据体中提取岩石物性参数的方法研究.数值测试表明,所讨论的参数提取方法具有一定的理论和应用价值.