VTI介质qP波方程高精度有限差分算子   总被引：10，自引：4，他引：6       下载免费PDF全文

吴国忱  梁锴 《地球物理学进展》2007,22(3):896-904

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题,会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散,本文针对VTI介质地震波数值模拟问题,构造了频率-空间域qP波波动方程高精度有限差分优化算子,根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析,并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟,通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度,有效压制了传统有限差分算子数值模拟中的数值频散现象,提高了有限差分算子精度,为VTI介质频率-空间域qP波正演模拟奠定了基础.  相似文献   

One‐dimensional dynamics of saturated incompressible porous media: analytical solutions and influence of inertia terms          下载免费PDF全文

Yunpeng Zhang  Dorival M. Pedroso  Wolfgang Ehlers 《国际地质力学数值与分析法杂志》2016,40(18):2489-2513

The complexity of formulations for the hydromechanical coupled mechanics of porous media is typically minimised by simplifying assumptions such as neglecting the effect of inertia terms. For example, three formulations commonly employed to model practical problems are classified as fully dynamic, simplified dynamic and quasi‐static. Thus, depending on the porous media conditions, each formulation will have advantages and limitations. This paper presents a comprehensive analysis of these limitations when solving one‐dimensional fully saturated porous media problems in addition to a new solution that considers a more general loading situation. A phase diagram is developed to assist on the selection of which formulation is more appropriate and convenient regarding particular cases of porosity and hydraulic conductivity values. Non‐dimensional formulations are proposed to achieve this goal. Results using the analytical solutions are compared against numerical values obtained with the finite element method, and the effect of porosity is investigated. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Numerical analysis of frost effects in porous media. Benefits and limits of the finite element poroelasticity formulation     

Stéphane Multon  Alain Sellier  Bernard Perrin 《国际地质力学数值与分析法杂志》2012,36(4):438-458

The aim of this paper is to analyse the performance of a finite element formulation usable for predicting the mechanical consequences of frost effects on porous media. It considers the characteristics of porous media and how the frost action can be assessed. The problem is then separated into two parts: thermal and poromechanical calculations. The constitutive equations developed in the framework of poromechanics are presented and the implementation in a usual finite element poroelasticity formulation based on Zuber's method is adopted. An analysis of the time‐step influence on the convergence rate is given and leads us to propose a simple method in order to obtain objectivity of the finite element response and avoid over‐long calculations. Frost effect simulations are carried out on real porous media (two fired clays) as a case study. Although the experimental behaviour of the porous media subjected to frost action is in accordance with some observations, the calculated strains appear to be overestimated compared with measurements. The problem could be largely attributable to the difficulty of assessing permeability evolution during frost development. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Internal erosion of chemically reinforced granular materials: a mathematical modeling approach     

Tony Khalil  Nadia Saiyouri  Bogdan Muresan  Pierre‐Yves Hicher 《国际地质力学数值与分析法杂志》2013,37(5):491-502

Internal erosion (IE) affects the stability of natural and reinforced materials by causing instability within their granular structure. The dislodgement and transport of eroded particles affect both the particulate concentration of eroding fluid and the pore network of eroded material. In this study, we examined these modifications using a transport model with a finite element code. First, IE tests on chemically reinforced sand columns were performed to obtain information about eroded material loss of mass, particulate concentration of effluent, porosity and permeability modifications, and existing IE stages. Second, based on experimental results, a mathematical one‐dimensional model has been formulated to monitor the evolution and spatial distribution of erodible solids, fluidized particles, porosity, permeability, and seepage stresses. The model consists of a set of coupled nonlinear differential equations solved in sequence. It provides valuable information about the extent and the dynamics of structural changes, which can be used to estimate an IE time for the hydraulic work to reach failure. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Dynamics of unsaturated poroelastic solids at finite strain     

Ryosuke Uzuoka  Ronaldo I. Borja 《国际地质力学数值与分析法杂志》2012,36(13):1535-1573

We derive the governing equations for the dynamic response of unsaturated poroelastic solids at finite strain. We obtain simplified governing equations from the complete coupled formulation by neglecting the material time derivative of the relative velocities and the advection terms of the pore fluids relative to the solid skeleton, leading to a so‐called u^s ? p^w ? p^a formulation. We impose the weak forms of the momentum and mass balance equations at the current configuration and implement the framework numerically using a mixed finite element formulation. We verify the proposed method through comparison with analytical solutions and experiments of quasi‐static processes. We use a neo‐Hookean hyperelastic constitutive model for the solid matrix and demonstrate, through numerical examples, the impact of large deformation on the dynamic response of unsaturated poroelastic solids under a variety of loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Model for degradation‐induced settlements as part of a coupled landfill model          下载免费PDF全文

S. Bente  V. Krase  U. Kowalsky  D. Dinkler 《国际地质力学数值与分析法杂志》2017,41(12):1390-1410

  相似文献   

二维不均匀介质中点源P-SV波响应的有限差分近似算法   总被引：1，自引：0，他引：1       下载免费PDF全文

谢小碧  姚振兴 《地球物理学报》1988,31(5):540-555

本文从一阶方程组形式的波动方程出发,发展了一种计算二维不均匀介质中点源P-SV波响应的近似方法。该方法通过引入线分布的应力作为震源,利用二维有限差分方法计算出线源响应,然后再经过波形校正和几何扩散校正得出相应的近似点源响应。通过把波形和振幅与精确解比较表明,该方法具有较好的精度。由于有限差分方法对于介质中速度和密度的分布没有特殊要求,另一方面,本文所给出的震源可以适用于位错点源、爆炸源或集中力源,因此上述方法十分适合于研究横向不均匀介质中的近场强地运动、爆炸振动或地震勘探等问题。  相似文献   

Seismic wave attenuation in fluid-saturated porous media     

James G. Berryman 《Pure and Applied Geophysics》1988,128(1-2):423-432

Contrary to the traditional view, seismic attenuation in Biot's theory of fluid-saturated porous media is due to viscous damping of local (not global) pore-fluid motion. Since substantial inhomogeneities in fluid permeability of porous geological materials are to be expected, the regions of highest local permeability contribute most to the wave energy dissipation while those of lowest permeability dominate the fluid flow rate if they are uniformly distributed. This dichotomy can explain some of the observed discrepancies between computed and measured attenuation of compressional and shear waves in porous earth. One unfortunate consequence of this result is the fact that measured seismic wave attenuation in fluid-filled geological materials cannot be used directly as a diagnostic of the global fluid-flow permeability.  相似文献   

Overall Properties of A Material With Inclusions Or Cavities   总被引：1，自引：0，他引：1  

J. A. Hudson 《Geophysical Journal International》1994,117(2):555-561

  相似文献   

Acoustic reflection calculation using transfer matrices     

G. W. Farnell 《Geophysical Journal International》1994,118(2):358-368

  相似文献