考虑土体固液二相性质的地下管线地震反应研究

本文将管线周围的土体视为固-液二相介质，采用考虑地震过程中土体内孔隙水压力的增长和消散的有效应力分析方法，并以非线性的土体本构模型模拟土体动力特性的改变以及在孔隙水压作用下的变形，以求接近土体的真实地震反应。研究表明随着土体软化程度的不同，受到土体约束的管线应力亦发生相应变化。文中详细研究了地震动峰值、管线直径、行波波速和入射角度等对管线应力的影响。     

Numerical simulation of seismic wave propagation in complex media by convolutional differentiator

We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.     

横观各向同性饱和两相介质动力反应计算分析

针对横观各向同性饱和两相介质的弹性波动方程组,应用基于显式逐步积分计算格式的时域显式有限元方法对其在输入地震波作用下的动力反应进行计算和分析,并将计算结果与完全各向同性饱和两相介质的计算结果进行对比研究。结果显示横观各向同性饱和两相介质与完全各向同性饱和两相介质的动力反应具有较为显著的差异。同时表明时域显式有限元方法是进行横观各向同性饱和两相介质动力反应计算分析的一种有效方法。     

Variational analysis of dynamics of fissured poroelastic rocks

Rocks can be modeled in a continuum framework as fissured, poroelastic materials, i.e., materials with two degrees of porosity, one due to the fissures and another one due to the pores. The governing equations of motion of fissured poroelastic rocks established by Beskos are rederived here by establishing a variational statement which also provides the boundary conditions of the problem. This is accomplished by considering strain, dissipation and kinetic energies as well as the work of external forces. The above statement is also derived here by employing the method of weighted residuals.     

Divergence of solutions to perturbation-based advection-dispersion moment equations

Moment equation methods are popular and powerful tools for modeling transport processes in randomly heterogeneous porous media, but the application of these methods to advection-dispersion equations often leads to erroneous oscillations. Perturbative methods, required to close systems of moment equations, become inaccurate for large perturbations; however, little quantitative theory exists for determining when this occurs for advection-dispersion equations. We consider three different methods (asymptotic approximation, Eulerian truncation, and iterative solution) for closing and solving advection-dispersion moment equations describing transport in stratified porous media with random permeability. We obtain approximate analytical expressions for time above which the asymptotic approximation to the mean diverges, in particular quantifying the impact that dispersion has on delaying—but not eliminating—divergence. We demonstrate that Eulerian truncation and iterative solution methods do not eliminate divergent behavior either. Our divergence criteria provide a priori estimates that signal a warning to the practitioner of stochastic advection-dispersion equations to carefully consider whether to apply perturbative approaches.     

DuMu: DUNE for multi-{phase, component, scale, physics, …} flow and transport in porous media

DuMu^x is a free and open-source simulator for flow and transport processes in porous media, based on the Distributed and Unified Numerics Environment (DUNE). Its main intention is to provide a sustainable and consistent framework for the implementation and application of model concepts, constitutive relations, discretizations, and solvers. The paper provides an overview of DuMu^x with the focus on software-related aspects. Selected examples highlight the multi-scale and the parallel capabilities.     

A semi-empirical approach to model pressure dependence of elastic moduli in granular media accounting for variations of coordination-number and Poisson-ratio

The effective medium theory based on the Hertz–Mindlin contact law is the most popular theory to relate dynamic elastic moduli (or elastic velocities) and confining pressure in dry granular media. However, many experimental results proved that the effective medium theory predicts pressure trends lower than experimental ones and over-predicts the shear modulus. To mitigate these mispredictions, several evolutions of the effective medium theory have been presented in the literature. Among these, the model named modified grain contact theory is an empirical approach in which three parametric curves are included in the effective medium theory model. Fitting the parameters of these curves permits to adjust the pressure trends of the Poisson ratio and the bulk modulus. In this paper, we present two variations of the modified grain contact theory model. First, we propose a minor modification in the fitting function for the porosity dependence of the calibration parameters that accounts for non-linearity in the vicinity of the critical porosity. Second, we propose a major modification that reduces the three-step modified grain contact theory model to a two-step model, by skipping the calibration parameter–porosity fit in the model and directly modelling the calibration parameter–pressure relation. In addition to an increased simplicity (the fitting parameters are reduced from 10 to 6), avoiding the porosity fit permits us to apply the model to laboratory data that are not provided with accurate porosity measurements. For this second model, we also estimate the uncertainty of the fitting parameters and the elastic velocities. We tested this model on dry core measurements from literature and we verified that it returns elastic velocity trends as good as the original modified grain contact theory model with a reduced number of fitting parameters. Possible developments of the new model to add predictive power are also discussed.     

地壳介质各向异性弹性本构关系讨论

论述了地壳介质中存在的各类向异性及其弹性本构关系,推导了应变能变换公式、ＥＤＡ薄层各向异性弹性常数及非对称面上ＥＤＡ构成的单斜对称弹性常数的计算公式;讨论了几种各向异性对称系统的本征方程并给出了弹性波的角散表达式,最后选取合理的参数对波的角散现象作了数值模拟。     

Seismic wave attenuation in fluid-saturated porous media

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.     

地震各向异性——多组裂隙对横波偏振的影响

通过对多分量地震资料的分析,我们发现随着频率的增加横波分裂时差减小.对于深部接收的VSP数据来说快横波的偏振方向保持不变,而对于浅层接收的VSP数据来说偏振方向却存在一个最大可以达到20°的旋转.尽管多尺度随机分布微裂隙岩石物理模型已经成功地模拟并解释了横波分裂时差随频率变化的现象,却不能解释与频率相关的横波分裂.据推测,如果微裂隙的排列方向和大裂隙的排列方向不同,利用低频信息获得的偏振方向将指示裂隙主方向,而利用高频信息获得的偏振方向则指示微裂隙方向.在背景多孔隙介质中存在多组裂隙的情况下,推导出垂直入射条件下横波偏振方向的解析式,给出了系统研究横波在介质中传播的方法.研究结果表明,横波偏振方向会随着频率的变化而变化,并且在入射方位、角度一定的条件下,是裂隙方位和密度的函数,这些认识可能有助于揭示观测到的、依赖频率变化的横波偏振现象.