Scattering of plane SV waves by cylindrical canyons in saturated porous medium   总被引：4，自引：0，他引：4  

Li Weihua  Zhao Chenggang 《Soil Dynamics and Earthquake Engineering》2005,25(12):981-995

On the basis of Biot dynamic theory, an analytic solution of two-dimensional scattering and diffraction of plane SV waves by circular cylindrical canyons in a half space of saturated porous media is presented in this paper for the first time. The solution is obtained by employing the Fourier–Bessel series expansion technique. Parametric studies had been carried out, which includes: the angle of incidence, the frequency of the incident SV wave, the porosity of saturated porous medium and the stiffness and Poisson's ratio of the solid-skeleton. All the outcomes are useful for the seismic analysis of the surface topography conditions.  相似文献   

Penny-shaped fractures revisited     

V. Grechka 《Studia Geophysica et Geodaetica》2005,49(3):365-381

All methods of seismic characterization of fractured reservoirs are based on effective media theories that relate geometrical and material properties of fractures and surrounding rock to the effective stiffnesses. In exploration seismology, the first-order theory of Hudson is the most popular. It describes the effective model caused by the presence of a single set of thin, aligned vertical fractures in otherwise isotropic rock. This model is known to be transversely isotropic with a horizontal symmetry axis (HTI). Following the theory, one can invert the effective anisotropy for the crack density and type of fluid infill of fractures, the quantities of great importance for reservoir appraisal and management.Here I compute effective media numerically using the finite element method. I deliberately construct models that contain a single set of vertical, ellipsoidal, non-intersecting and non-interconnected fractures to check validity of the first-order Hudson’s theory and establish the limits of its applicability. Contrary to conventional wisdom that Hudson’s results are accurate up to crack density e ≈ 0.1, I show that they consistently overestimate the magnitudes of all effective anisotropic coefficients ε^(V), δ^(V), and γ^(V). Accuracy of theoretically derived anisotropy depends on the type of fluid infill and typically deteriorates as e grows. While the theory gives | ε^(V)|, |δ^(V)|, |γ^(V)| and close to the upper bound of the corresponding numerically obtained values for randomly distributed liquid-filled fractures, theoretical predictions of ε^(V), δ^(V) are not supported by numerical computations when the cracks are dry. This happens primarily because the first-order Hudson’s theory makes no attempt to account for fracture interaction which contributes to the final result much stronger for gas- than for liquid-filled cracks. I find that Mori-Tanaka’s theory is superior to Hudson’s for all examined crack densities and both types of fluid infill.The paper was presented at the 11th International Workshop on Seismic Anisotropy (11IWSA) held in St. John’s, Canada in 2004.  相似文献   

Modelling of multiple events using empirical Green''s functions: method, application to swarm earthquakes and implications for their rupture propagation     

Tomá&#; Fischer 《Geophysical Journal International》2005,163(3):991-1005

  相似文献   

General expressions for internal deformation fields due to a dislocation source in a multilayered elastic half-space     

Yukitoshi Fukahata  Mitsuhiro Matsu'ura 《Geophysical Journal International》2005,161(2):507-521

  相似文献   

Détermination expérimentale des paramètres équivalents d'un milieu poreux hétérogène en écoulements uniforme ou radial     

Charles Danquigny  Philippe Ackerer 《Comptes Rendus Geoscience》2005,337(6):563-570

Tracer tests are carried out in a heterogeneous porous medium that has a 3D correlated random distribution of the permeabilities. The fitting of numerical models provides the values of equivalent permeability and macrodispersivity characterizing a 2D homogeneous horizontal medium. Different flow configurations are studied: uniform, radial and pump and treat (doublet). The fitted parameter sets are independent of the flow type, except for the doublet. They are greater than the values predicted by stochastic theories, due to the small number of correlation lengths explored by the tracer and the limited extension of the experimental set-up. To cite this article: C. Danquigny, P. Ackerer, C. R. Geoscience 337 (2005).  相似文献   

Mise au point d'un milieu de culture pour l'étude de l'altération de silicates en présence de Pseudomonas aeruginosa     

Georges Aouad  Valérie Geoffroy  Jean-Marie Meyer  Jean-Louis Crovisier  Denis Damidot  Peter Stille 《Comptes Rendus Geoscience》2005,337(15):1340-1347

The influence of microorganisms on mineral alteration is not easy to determine in environmental conditions, because of the difficulty to raise for comparison purposes an identical but abiotic system. Another problem in this context is the choice of reliable tracers to evaluate the alteration rate of materials during in vitro experiments. To face such difficulties, we elaborated a defined medium allowing both the growth of Pseudomonas aeruginosa and a precise measurement of the elements solubilized from the minerals. Thanks to this medium, we were able to quantitatively determine the amounts of major elements solubilized from the materials in the presence of bacterial growth, compared to a sterile system. Moreover, the analysis by ICP-MS of trace elements was possible after a chromatographic treatment, which selectively eliminated 99% of the sodium content of the medium. To cite this article: G. Aouad et al., C. R. Geoscience 337 (2005).  相似文献   

Fracture mechanics analysis of multiple cracks in anisotropic media     

Chia‐Huei Tu  C.‐S. Chen  T.‐T. Yu 《国际地质力学数值与分析法杂志》2011,35(11):1226-1242

This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation          下载免费PDF全文

Chongbin Zhao  Thomas Poulet  Klaus Regenauer‐Lieb 《国际地质力学数值与分析法杂志》2015,39(6):571-593

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

A model for moisture and heat flow in fractured unsaturated porous media          下载免费PDF全文

A.M. Varnosfaderani  B. Gatmiri  E.M. Haghighi 《国际地质力学数值与分析法杂志》2017,41(6):828-858

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Hydrologic response of engineered media in living roofs and bioretention to large rainfalls: experiments and modeling          下载免费PDF全文

Ruifen Liu  Elizabeth Fassman‐Beck 《水文研究》2017,31(3):556-572

The hydrologic response of engineered media plays an important role in determining a stormwater control measure's ability to reduce runoff volume, flow rate, timing, and pollutant loads. Five engineered media, typical of living roof and bioretention stormwater control measures, were investigated in laboratory column experiments for their hydrologic responses to steady, large inflow rates. The inflow, medium water content response, and outflow were all measured. The water flow mechanism (uniform flow vs. preferential flow) was investigated by analyzing medium water content response in terms of timing, magnitude, and sequence with depth. Modeling the hydrologic process was conducted in the HYDRUS‐1D software, applying the Richards equation for uniform flow modeling, and a mobile–immobile model for preferential flow modeling. Uniform flow existed in most cases, including all initially dry living roof media with bimodal pore size distributions and one bioretention medium with unimodal pore size distribution. The Richards equation can predict the outflow hydrograph reasonably well for uniform flow conditions when medium hydraulic properties are adequately represented by appropriate functions. Preferential flow was found in two media with bimodal pore size distributions. The occurrence of preferential flow is more likely due to the interaction between the bimodal pore structure and the initial water content rather than the large inflow rate.  相似文献