共查询到20条相似文献，搜索用时 31 毫秒
1.
Diganta Bhusan Das 《水文研究》2002,16(17):3393-3418
Hydrodynamic modelling for analysis of groundwater flow through permeable reactive barriers (PRBs) is addressed in this paper. Permeable reactive barriers constitute an emerging technology for in situ remediation of groundwater contamination and have many advantages over the traditional ex situ treatment methods. The transport domains during groundwater flow through PRBs often may involve free‐flow or non‐porous sections. To model the fluid mobility efficiently in such situations, the free and porous flow zones (PRBs) must be studied in conjunction with each other. The present paper is devoted to the analysis of groundwater flow through combined free flow domains and PRBs. The free‐flow regime is modelled using the Navier–Stokes equations whereas the permeable barriers are simulated by either the Darcy or the Brinkman equation. In order to couple the governing equations of motions, well‐posed mathematical formulations of matching boundary conditions are prescribed at the interface between the free‐groundwater‐flow zones and the permeable barriers. Combination of the Navier–Stokes equations with the Brinkman equation is more straightforward owing to their analogous forms. However, the Navier–Stokes and Darcy equations are incompatible mathematically and cannot be linked directly. The problem is resolved in this paper by invoking validated hydrodynamical expressions for describing the flow behaviour at the interfaces between free‐flow and porous zones. Three schemes for the analyses of fluid flow in combined domains are applied to the case of groundwater flow through permeable reactive barriers and different model results are compared. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
2.
3.
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. 相似文献
4.
Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with small-scale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference (FD) method of Biot’s poroelastic equations, with an arbitrary even-order (2L) accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer (CPML) absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultrasonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thicknesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Q sc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments. 相似文献
5.
An absorbing boundary for saturated porous media is developed that can be used for transient analyses in the time domain. The elastic constitutive equations for the saturated porous media follow Bowen's formulation. The method consists of applying viscous tractions along the artificial boundary. The absorbing boundary behaviour is assumed linear and isotropic. Hadamard's conditions provide the speeds of the dilatational and shear waves that propagate in saturated porous media. Since these expressions are frequency independent, the intensities of the viscous tractions are evaluated in the time domain, and the two dilatational waves are accounted for. The viscous tractions are defined from the drained characteristics, assuming an infinite permeability, at variance with the traditional ‘undrained’ method based on undrained characteristics and a null permeability. Solid media and materials with low permeability are also retrieved as subcases. The results show that, at no additional cost, this ‘drained’ method is more accurate for all permeabilities than the ‘undrained’ method, which disregards the existence of the second dilatational wave. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
6.
In subsurface porous media, the soil water retention curve (WRC) and unsaturated hydraulic conductivity curve (UHC) are two important soil hydraulic property curves. Spatial heterogeneity is ubiquitous in nature, which may significantly affect soil hydraulic property curves. The main theme of this paper is to investigate how spatial heterogeneities, including their arrangements and amounts in soil flumes, affect soil hydraulic property curves. This paper uses a two‐dimensional variably saturated flow and solute transport finite element model to simulate variations of pressure and moisture content in soil flumes under a constant head boundary condition. To investigate the behavior of soil hydraulic property curves owing to variations of heterogeneities and their arrangements as well, cases with different proportions of heterogeneities are carried out. A quantitative evaluation of parameter variations in the van Genuchten model (VG model) resulting from heterogeneity is presented. Results show that the soil hydraulic properties are strongly affected by variations of heterogeneities and their arrangements. If the pressure head remains at a specific value, the soil moisture increases when heterogeneities increase in the soil flumes. On the other hand, the unsaturated hydraulic conductivity decreases when heterogeneities increase in the soil flumes under a constant pressure head. Moreover, results reveal that parameters estimated from both WRC and UHC also are affected by shapes of heterogeneity; this indicates that the parameters obtained from the WRC are not suitable for predicting the UHC of different shapes in heterogeneous media. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
7.
《Advances in water resources》2004,27(6):657-667
In this article, the quadrupole method is implemented in order to simulate the effects of heterogeneities on one dimensional advective and diffusive transport of a passive solute in porous media. Theoretical studies of dispersion in heterogeneous stratified media can bring insight into transport artefacts linked to scale effects and apparent dispersion coefficients. The quadrupole method is an efficient method for the calculation of transient response of linear systems. It is based here on the Laplace transform technique. The analytical solutions that can be derived by this method assists understanding of upscaled parameters relevant to heterogeneous porous media.First, the method is developed for an infinite homogeneous porous medium. Then, it is adapted to a stratified medium where the fluid flow is perpendicular to the interfaces. The first heterogeneous medium studied is composed of two semi-infinite layers perpendicular to the flow direction each having different transport properties. The concentration response of the medium to a Dirac injection is evaluated. The case studied emphasises the importance in the choice of the boundary conditions.In the case of a periodic heterogeneous porous medium, the concentration response of the medium is evaluated for different numbers of unit-cells. When the number of unit cells is great enough, depending on the transport properties of each layer in the unit cell, an equivalent homogeneous behaviour is reached. An exact determination of the transport properties (equivalent dispersion coefficient) of the equivalent homogeneous porous medium is given. 相似文献
8.
Oversampling techniques are often used in porous media simulations to achieve high accuracy in multiscale simulations. These methods reduce the effect of artificial boundary conditions that are imposed in computing local quantities, such as upscaled permeabilities or basis functions. In the problems without scale separation and strong non-local effects, the oversampling region is taken to be the entire domain. The basis functions are computed using single-phase flow solutions which are further used in dynamic two-phase simulations. The standard oversampling approaches employ generic global boundary conditions which are not associated with actual flow boundary conditions. In this paper, we propose a flow based oversampling method where the actual two-phase flow boundary conditions are used in constructing oversampling auxiliary functions. Our numerical results show that the flow based oversampling approach is several times more accurate than the standard oversampling method. We provide partial theoretical explanation for these numerical observations. 相似文献
9.
Application of a perfectly matched layer in seismic wavefield simulation with an irregular free surface 总被引：2，自引：0，他引：2 
下载免费PDF全文
下载免费PDF全文 Haiqiang Lan Jingyi Chen Zhongjie Zhang Youshan Liu Jianguo Zhao Ruiqi Shi 《Geophysical Prospecting》2016,64(1):112-128
Recently, an effective and powerful approach for simulating seismic wave propagation in elastic media with an irregular free surface was proposed. However, in previous studies, researchers used the periodic condition and/or sponge boundary condition to attenuate artificial reflections at boundaries of a computational domain. As demonstrated in many literatures, either the periodic condition or sponge boundary condition is simple but much less effective than the well‐known perfectly matched layer boundary condition. In view of this, we intend to introduce a perfectly matched layer to simulate seismic wavefields in unbounded models with an irregular free surface. We first incorporate a perfectly matched layer into wave equations formulated in a frequency domain in Cartesian coordinates. We then transform them back into a time domain through inverse Fourier transformation. Afterwards, we use a boundary‐conforming grid and map a rectangular grid onto a curved one, which allows us to transform the equations and free surface boundary conditions from Cartesian coordinates to curvilinear coordinates. As numerical examples show, if free surface boundary conditions are imposed at the top border of a model, then it should also be incorporated into the perfectly matched layer imposed at the top‐left and top‐ right corners of a 2D model where the free surface boundary conditions and perfectly matched layer encounter; otherwise, reflections will occur at the intersections of the free surface and the perfectly matched layer, which is confirmed in this paper. So, by replacing normal second derivatives in wave equations in curvilinear coordinates with free surface boundary conditions, we successfully implement the free surface boundary conditions into the perfectly matched layer at the top‐left and top‐right corners of a 2D model at the surface. A number of numerical examples show that the perfectly matched layer constructed in this study is effective in simulating wave propagation in unbounded media and the algorithm for implementation of the perfectly matched layer and free surface boundary conditions is stable for long‐time wavefield simulation on models with an irregular free surface. 相似文献
10.
Different theoretical and laboratory studies on the propagation of elastic waves in layered hydrocarbon reservoir have shown characteristic velocity dispersion and attenuation of seismic waves. The wave‐induced fluid flow between mesoscopic‐scale heterogeneities (larger than the pore size but smaller than the predominant wavelengths) is the most important cause of attenuation for frequencies below 1 kHz. Most studies on mesoscopic wave‐induced fluid flow in the seismic frequency band are based on the representative elementary volume, which does not consider interaction of fluid flow due to the symmetrical structure of representative elementary volume. However, in strongly heterogeneous media with unsymmetrical structures, different courses of wave‐induced fluid flow may lead to the interaction of the fluid flux in the seismic band; this has not yet been explored. This paper analyses the interaction of different courses of wave‐induced fluid flow in layered porous media. We apply a one‐dimensional finite‐element numerical creep test based on Biot's theory of consolidation to obtain the fluid flux in the frequency domain. The characteristic frequency of the fluid flux and the strain rate tensor are introduced to characterise the interaction of different courses of fluid flux. We also compare the behaviours of characteristic frequencies and the strain rate tensor on two scales: the local scale and the global scale. It is shown that, at the local scale, the interaction between different courses of fluid flux is a dynamic process, and the weak fluid flux and corresponding characteristic frequencies contain detailed information about the interaction of the fluid flux. At the global scale, the averaged strain rate tensor can facilitate the identification of the interaction degree of the fluid flux for the porous medium with a random distribution of mesoscopic heterogeneities, and the characteristic frequency of the fluid flux is potentially related to that of the peak attenuation. The results are helpful for the prediction of the distribution of oil–gas patches based on the statistical properties of phase velocities and attenuation in layered porous media with random disorder. 相似文献
11.
Laboratory Study of Air Sparging: Air Flow Visualization 总被引：15，自引：0，他引：15
Wei Ji Amine Dahmani David P. Ahlfeld Jia Ding Lin Edward Hill III 《Ground Water Monitoring & Remediation》1993,13(4):115-126
Laboratory flow visualization experiments, using glass beads as the porous medium, were conducted to study air sparging, an innovative technology for subsurface contaminant remediation. The purpose of these experiments was to observe how air flows through saturated porous media and to obtain a basic understanding of air plume formation and medium heterogeneity effects. The experiments indicate that air flow occurring in discrete, stable channels is the most probable flow behavior in medium to fine grained water saturated porous media and that medium heterogeneity plays an important role in the development of air channels. Several simulated scales of heterogeneities, from pore to field, have been studied. The results suggest that air channel formation is sensitive to the various scales of heterogeneities. Site-specific hydrogeologic settings have to be carefully reviewed before air sparging is applied to remediate sites contaminated by volatile organic compounds. 相似文献
12.
13.
Saturation of porous rocks with a mixture of two fluids has a substantial effect on seismic‐wave propagation. In particular, partial saturation causes significant attenuation and dispersion of the propagating waves due to the mechanism of wave‐induced fluid‐flow. Such flow arises when a passing wave induces different fluid pressures in regions of rock saturated by different fluids. Most models of attenuation and dispersion due to mesoscopic heterogeneities imply that fluid heterogeneities are distributed in a regular way. However, recent experimental studies show that mesoscopic heterogeneities have less idealized distributions and that the distribution itself affects attenuation and dispersion. Based on an approximation for the coherent wavefield in random porous media, we develop a model which assumes a continuous distribution of fluid heterogeneities. As this continuous random media approach assumes that there will be a distribution of different patch sizes, it is expected to be better suited to modelling experimental data. We also show how to relate the random functions to experimentally measurable parameters. 相似文献
14.
In general, the accuracy of numerical simulations is determined by spatial and temporal discretization levels. In fractured porous media, the time step size is a key factor in controlling the solution accuracy for a given spatial discretization. If the time step size is restricted by the relatively rapid responses in the fracture domain to maintain an acceptable level of accuracy in the entire simulation domain, the matrix tends to be temporally over-discretized. Implicit sub-time stepping applies smaller sub-time steps only to the sub-domain where the accuracy requirements are less tolerant and is most suitable for problems where the response is high in only a small portion of the domain, such as within and near the fractures in fractured porous media. It is demonstrated with illustrative examples that implicit sub-time stepping can significantly improve the simulation efficiency with minimal loss in accuracy when simulating flow and transport in fractured porous media. The methodology is successfully applied to density-dependent flow and transport simulations in a Canadian Shield environment, where the flow and transport is dominated by discrete, highly conductive fracture zones. 相似文献
15.
T. AkiyoshiK. FuchidaH. L. Fang 《Soil Dynamics and Earthquake Engineering》1994,13(6):387-397
Based on Biot's two-phase mixture theory and the paraxial approximation, the absorbing boundary conditions in the time domain for u-w, u-U and u-p formulations are presented for the dynamic analysis of fluid-saturated porous media. These absorbing boundaries are equivalent to the viscous boundaries in the fundamental mode. The expressions for the energy ratio and reflection coefficient between the reflected and incident waves along the absorbing boundary are given. The numerical results show that the proposed absorbing boundaries can greatly suppress the spuriously reflected wave and model the far field of the system. These results also dynamic analysis of infinite fluid-saturated porous media. for the transient dynamic analysis of infinite fluid-saturated porous media. 相似文献
16.
Generation of three dimensional flow fields for statistically anisotropic heterogeneous porous media
M. A. Cushy A. Bellin Y. Rubin 《Stochastic Environmental Research and Risk Assessment (SERRA)》1995,9(1):89-104
A methodology for generating three dimensional (3D) flow fields for statistically anisotropic heterogeneous porous media is presented and demonstrated. The simulated flow fields are shown to exhibit the input spatial correlation structure and observe mass continuity. Sample flow fields are presented in the form of cross sectional slices of the 3D formation. These cross sections demonstrate visually the characteristics of subsurface flow. The method was found to be faster than traditional techniques in terms of its computational requirements. Given this method, it is possible to generate the large number of realizations of a velocity field necessary to compute high order statistics in transport problems. 相似文献
17.
地震波传播激发的不同尺度的流固相对运动(宏观、中观和微观)是许多沉积岩地层中地震波频散和衰减的主要原因,然而野外观测和试验测量都难以对非均匀多孔介质孔隙压力弛豫物理过程进行精细刻画.通过数字岩石物理技术,本文建立了三个典型的数字岩心分别用于表征孔隙结构、岩石骨架和斑状饱和流体引起的非均质性,利用动态应力应变模拟技术计算数字岩心的位移和孔隙流体增量图像.通过分析和比较三个数字岩心的位移和孔隙压力增量图像,细致刻画了发生于非均匀含流体多孔介质内的宏观、中观和微观尺度的流固相对运动:1)宏观尺度的波致孔隙流体流动导致波长尺度上数字岩心不同区域的孔隙压力和位移差异;2)中观尺度的流体流动发生在软层与硬层之间、气层与液层之间;3)微观尺度的流体流动发生在孔隙内部或相邻孔隙之间.数值模拟试验也证明基于数字岩心的动态应力应变模拟技术可以从微观尺度上更好的理解波致孔隙流体流动发生的物理机理,从而为建立岩石骨架、孔隙流体、孔隙结构非均质性和弹性波频散-衰减特征的映射关系奠定基础. 相似文献
18.
Vittorio Di Federico Marco Pinelli Rita Ugarelli 《Stochastic Environmental Research and Risk Assessment (SERRA)》2010,24(7):1067-1076
An understanding of the interplay between non-Newtonian effects in porous media flow and field-scale domain heterogeneity
is of great importance in several engineering and geological applications. Here we present a simplified approach to the derivation
of an effective permeability for flow of a purely viscous power–law fluid with flow behavior index n in a randomly heterogeneous porous domain subject to a uniform pressure gradient. A standard form of the flow law generalizing
the Darcy’s law to non-Newtonian fluids is adopted, with the permeability coefficient being the only source of randomness.
The natural logarithm of the permeability is considered a spatially homogeneous and correlated Gaussian random field. Under
the ergodic hypothesis, an effective permeability is first derived for two limit 1-D flow geometries: flow parallel to permeability
variation (serial-type layers), and flow transverse to permeability variation (parallel-type layers). The effective permeability
of a 2-D or 3-D isotropic domain is conjectured to be a power average of 1-D results, generalizing results valid for Newtonian
fluids under the validity of Darcy’s law; the conjecture is validated comparing our results with previous literature findings.
The conjecture is then extended, allowing the exponents of the power averaging to be functions of the flow behavior index.
For Newtonian flow, novel expressions for the effective permeability reduce to those derived in the past. The effective permeability
is shown to be a function of flow dimensionality, domain heterogeneity, and flow behavior index. The impact of heterogeneity
is significant, especially for shear-thinning fluids with a low flow behavior index, which tend to exhibit channeling behavior. 相似文献
19.
Clogging of saturated porous media by silt‐sized suspended solids under varying physical conditions during managed aquifer recharge 下载免费PDF全文
下载免费PDF全文 Managed aquifer recharge is an effective method for utilizing excess flood flows, but clogging of porous media is a limiting factor in the implementation of this water storage technique. In recent years, much research on the physical clogging of porous media during artificial recharge has been conducted. However, the understanding of clogging due to silt‐sized suspended solids (SS) is still inadequate, especially under varying physical conditions. Here, we subjected sand columns to controlled rates of flow and SS suspensions to investigate the influence of media size, SS size, SS concentration, and flow velocity on the clogging of porous media by silt‐sized SS. The results show that the diameter ratio of SS particles to sand grains is the dominant factor influencing the position of physical clogging. As pore velocity increased, the mobility of silt‐sized SS was enhanced and retention in the porous media decreased noticeably. The spatial retention profiles in the porous media were found to vary greatly at different flow velocities. The SS concentration of the infiltrating suspension also dramatically influenced the mobility and deposition of silt‐sized SS particles, such that high concentrations accelerated the clogging process. As the different physical factors changed, the breakthrough curves and retention profiles of silt‐sized SS particles changed obviously and the mechanisms of retention differed. On the whole, clogging position is mainly determined by particle size ratio, but clogging rate is dominated by a variety of factors including particle size ratio, SS concentration, and flow velocity. 相似文献
20.
The effects of ice cover on flow characteristics in meandering rivers are still not completely understood. Here, we quantify the effects of ice cover on flow velocity, the vertical and spatial flow distribution, and helical flow structure. Comparison with open‐channel low flow conditions is performed. An acoustic doppler current profiler (ADCP) is used to measure flow from up to three meander bends, depending on the year, in a small sandy meandering subarctic river (Pulmanki River) during two consecutive ice‐covered winters (2014 and 2015). Under ice, flow velocities and discharges were predominantly slower than during the preceding autumn open‐channel conditions. Velocity distribution was almost opposite to theoretical expectations. Under ice, velocities reduced when entering deeper water downstream of the apex in each meander bend. When entering the next bend, velocities increased again together with the shallower depths. The surface velocities were predominantly greater than bottom/riverbed velocities during open‐channel flow. The situation was the opposite in ice‐covered conditions, and the maximum velocities occurred in the middle layers of the water columns. High‐velocity core (HVC) locations varied under ice between consecutive cross‐sections. Whereas in ice‐free conditions the HVC was located next to the inner bank at the upstream cross‐sections, the HVC moved towards the outer bank around the apex and again followed the thalweg in the downstream cross‐sections. Two stacked counter‐rotating helical flow cells occurred under ice around the apex of symmetric and asymmetric bends: next to the outer bank, top‐ and bottom‐layer flows were towards the opposite direction to the middle layer flow. In the following winter, no clear counter‐rotating helical flow cells occurred due to the shallower depths and frictional disturbance by the ice cover. Most probably the flow depth was a limiting factor for the ice‐covered helical flow circulation, similarly, the shallow depths hinder secondary flow in open‐channel conditions. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献