饱和岩石的各向异性及非线性黏弹性响应

在MTS伺服压机上对不同饱和状态的砂岩、大理岩标本进行了垂直层理和平行层理两个方向的正弦波加载试验,研究了饱和岩石的各向异性及非线性黏弹性行为. 在最大载荷低于岩石的屈服点时,获得的衰减、杨氏模量、泊松比、波速都表现出显著的各向异性,同时它们还具有较明显的应变振幅效应和频率效应. 随着应变振幅的增大,衰减近似呈线性增长,杨氏模量和泊松比近似呈线性下降. 在0.005-4Hz频段内除衰减不明显依赖频率以外,杨氏模量、泊松比、波速对频率的依赖性都较强; 当频段扩大到15Hz时,这4个参数都表现出较强的频率效应. 饱和岩石的衰减、杨氏模量、泊松比、波速均随饱和液体的黏滞系数增大而增大.     

Frequency dependence of elastic wave speeds at high temperature: a direct experimental demonstration

At sufficiently high temperatures and/or long periods, the elastic behaviour of crystalline material gives way progressively to viscoelastic behaviour associated with the stress-induced migration of crystal defects. This transition is marked by the onset of appreciable strain energy dissipation accompanied by frequency dependence (dispersion) of the shear modulus and elastic wave speeds. Ultrasonic interferometry and torsional forced-oscillation techniques can be used to probe the low-amplitude stress-strain behaviour of fine-grained polycrystalline material in two very different frequency ranges, respectively 10-100 MHz and 1 mHz-1 Hz. Here we demonstrate and apply these two complementary methods in a study of the high-temperature mechanical behaviour of a fine-grained synthetic olivine polycrystal. At the high frequencies of ultrasonic interferometry, the shear wave speed varies linearly with temperature between room temperature and the highest experimentally accessible temperature (1300 °C) in close accord with expectations based on similarly high-frequency studies of the elastic behaviour of single-crystal olivine. However, at teleseismic frequencies (<1 Hz) and temperatures >900 °C, the shear wave speed becomes much more strongly temperature-sensitive reflecting markedly viscoelastic behaviour. Newly emerging laboratory-derived constraints on this viscoelastic enhancement of the temperature sensitivity of seismic wave speeds and its grain-size dependence will provide a more robust interpretation of seismological models for the variation of wave speeds and attenuation within the Earth's interior.     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。     

不同检测方法对混凝土波速的影响研究

为建立超声波法和低应变反射波法检测所得混凝土基桩波速的相互转换关系,本文给出了标准线性固体模型下粘弹性介质的三维和一维纵波速度之比及其实用简化公式,并且通过实验得出,混凝土构件的超声波速与1.054倍低应变动测波速之比处于1.03~1.12,具体数值会随着混凝土龄期、弹性模量和强度等级的增加而有所减小。     

基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明:本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础.     

常Q滞弹性介质地震波动数值模拟——时域本构优化逼近

在地震波动模拟中计入常Q滞弹性阻尼,可有效降低模拟波形的误差.就时域有限差分和有限元模拟而言,常基于广义标准线性体建立阻尼介质的时域本构逼近.广义标准线性体由若干标准线性体并联得到,增加标准线性体个数能有效提高模拟精度,但计算量及计算存储将成倍增长.目前尚未有普适的标准线性体个数优化取值方案.本文基于广义标准线性体参数的非线性优化拟合方法,详细分析了时域本构逼近误差的影响因素,清楚揭示了逼近误差仅取决于频带宽度,与频带上下限取值无关这一特性,阐明了构建具有普适性标准线性体个数优化取值方案的可行性.论证了波形模拟精度主要取决于波传播距离与模拟波长的比值以及标准线性体的个数取值.综合考虑上述两个控制因素,结合在波动正反演问题中广为采纳的波形时频误差衡量准则,对不同Q值介质给出了标准线性体个数优化取值表.进一步,本文提出采用不同个数标准线性体以近似不同Q值的阻尼介质时域本构,解决了以往波动数值模拟中统一采用相同个数标准线性体而导致的计算量及计算存储浪费或模拟精度低下等问题,并基于数值实验验证了这一方法的精度.本文工作对推进滞弹性介质波动数值模拟及其在全波形反演问题中的应用具有理论价值和实践意义.     

黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

有能量吸收的一维波动方程的阻抗反演

本文讨论了一维粘弹性波阻抗反演的迭代法。根据一维波动方程和ＫｅｌｖｉｎＶｏｉｇｔ粘弹模型导出了一个适于数值计算的第二类Ｖｏｌｔｅｒｒａ型积分方程，用迭代法求解。通过对已有方法的改进和完善，数值计算表明了该反演方法具有一定的精度、抗干扰能力和分辨率以及应用于实际水平层状地层反射地震记录反演阻抗的可能性．     

青藏高原更新世晚期冰川的均衡效应

本文基于Ｋｕｈｌｅ关于青藏更新世晚期存在统一大冰盖的假说，构制冰负荷史模型，利用Ｐｅｌｔｉｅｒ（１９７４，１９７６）粘弹地球表面负荷形变理论，计算出冰盖负荷引起的地表垂直位移和重力变化，并将其作为该区冰种均衡效应的最大估计，数值计算结果表明，冰盖负荷的均衡作用对高原隆升有一定的影响。     

Comparison between thermal–viscous coupling and frictional sliding

We investigated the similarity between thermal–viscous coupling (TVC) and frictional sliding, proposed by Kameyama and Kaneda [Pure Appl. Geophys. 159 (2002) 2011]. We consider a one-dimensional layer composed of viscous material, which is sandwiched and sheared by two thick elastic layers. The rate of viscous deformation depends on the temperature T_c in the viscous layer as well as shear stress τ. The temperature T_c changes owing to heating by viscous dissipation and conductive cooling. We carried out velocity-stepping tests for the steady-state deformation both numerically and analytically, and compared the temporal evolution of small perturbations with that of the spring-block model with rate- and state-dependent friction (RSF). We found that, as is the case of frictional slip stability, the manner of temporal evolution is classified into four regimes depending on whether it is stable or not and whether it is monotonous or oscillatory with time. By further interpreting TVC in terms of general RSF theory by Ruina [J. Geophys. Res. 88 (1983) 10359], we obtained the relations between the parameters appearing in the phenomenological RSF law and the nondimensional parameters which characterize the nature of TVC. A further improvement of this approach might be important for estimating the actual values of frictional constitutive parameters at the deeper portion of seismogenic faults of interplate or inland earthquakes where a ductile deformation is expected to be significant.