Approximate traveltime inversion in downhole microseismic monitoring

In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media.     

横观各向同性场地中埋置衬砌隧道对qP波的散射

采用间接边界元法(IBEM)结合"分区契合"技术,研究qP波入射下横观各向同性(TI)场地中衬砌隧道的动力响应问题。方法充分利用半空间和全空间动力格林函数在分别构造含孔半无限空间域和闭合域内散射波场方面的优势,将含有衬砌隧道的层状弹性半空间分解为含孔半无限空间域和一个环形衬砌闭合域来分别进行波场构造,有效地降低了求解时间和存储量。文中验证方法的正确性,并以均匀TI半空间和基岩上单一TI土层为例,计算分析弹性半空间场地中隧道衬砌内表面动应力放大问题。结果表明qP波入射下,TI介质与各向同性介质场地中埋置衬砌隧道的动力响应差异显著,TI介质参数的变化导致场地动力特性的改变,进而改变场地与衬砌隧道的动力相互作用机制,显著影响着衬砌内表面动应力的大小及其空间分布。     

Elastic solutions for a transversely isotropic half-space subjected to buried asymmetric-loads

Elastic closed-form solutions for the displacements and stresses in a transversely isotropic half-space subjected to various buried loading types are presented. The loading types include finite line loads and asymmetric loads (such as uniform and linearly varying rectangular loads, or trapezoidal loads). The planes of transverse isotropy are assumed to be parallel to its horizontal surface. These solutions are directly obtained from integrating the point load solutions in a transversely isotropic half-space, which were derived using the principle of superposition, Fourier and Hankel transformation techniques. The solutions for the displacements and stresses in transversely isotropic half-spaces subjected to linearly variable loads on a rectangular region are never mentioned in literature. These exact solutions indicate that the displacements and stresses are influenced by several factors, such as the buried depth, the loading types, and the degree and type of rock anisotropy. Two illustrative examples, a vertical uniform and a vertical linearly varying rectangular load acting on the surface of transversely isotropic rock masses, are presented to show the effect of various parameters on the vertical surface displacement and vertical stress. The results indicate that the displacement and stress distributions accounted for rock anisotropy are quite different for those calculated from isotropic solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用

由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性（VTI）介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。     

Inversion of multicomponent seismic time shifts for reservoir pressure and length: a feasibility study

Pressure drops associated with reservoir production generate excess stress and strain that cause travel‐time shifts of reflected waves. Here, we invert time shifts of P‐, S‐, and PS‐waves measured between baseline and monitor surveys for pressure reduction and reservoir length. The inversion results can be used to estimate compaction‐induced stress and strain changes around the reservoir. We implement a hybrid inversion algorithm that incorporates elements of gradient, global/genetic, and nearest neighbour methods and permits exploration of the parameter space while simultaneously following local misfit gradients. Our synthetic examples indicate that optimal estimates of reservoir pressure from P‐wave data can be obtained using the reflections from the reservoir top. For S‐waves, time shifts from the top of the reservoir can be accurately inverted for pressure if the noise level is low. However, if noise contamination is significant, it is preferable to use S‐wave data (or combined shifts of all three modes) from reflectors beneath the reservoir. Joint wave type inversions demonstrate improvements over any single pure mode. Reservoir length can be estimated using the time shifts of any mode from the reservoir top or deeper reflectors. We also evaluate the differences between the actual strain field and those corresponding to the best‐case inversion results obtained using P‐ and S‐wave data. Another series of tests addresses the inversion of the time shifts for the pressure drops in two‐compartment reservoirs, as well as for the associated strain field. Numerical testing shows that a potentially serious source of error in the inversion is a distortion in the strain‐sensitivity coefficients, which govern the magnitude of stiffness changes. This feasibility study suggests which wave types and reflector locations may provide the most accurate estimates of reservoir parameters from compaction‐induced time shifts.     

Dynamic analysis of slab track on multi-layered transversely isotropic saturated soils subjected to train loads

The dynamic responses of a slab track on transversely isotropic saturated soils subjected to moving train loads are investigated by a semi-analytical approach. The track model is described as an upper Euler beam to simulate the rails and a lower Euler beam to model the slab. Rail pads between the rails and slab are represented by a continuous layer of springs and dashpots. A series of point loads are formulated to describe the moving train loads. The governing equations of track-ground systems are solved using the double Fourier transform, and the dynamic responses in the time domain are obtained by the inverse Fourier transform. The results show that a train load with high velocity will generate a larger response in transversely isotropic saturated soil than the lower velocity load, and special attention should be paid on the pore pressure in the vicinity of the ground surface. The anisotropic parameters of a surface soil layer will have greater influence on the displacement and excess pore water pressure than those of the subsoil layer. The traditional design method taking ground soil as homogeneous isotropic soil is unsafe for the case of RE 1 and RG 1, so a transversely isotropic foundation model is of great significance to the design for high train velocities.     

考虑各向异性的层土-盾构隧道地震反应数值模拟

在层状各向异性土体-盾构隧道地震反应分析中,引入了横观各向同性弹塑性模型理论,建立了适合于横观各向同性介质的双渐近-多向透射边界条件。针对地铁区间盾构隧道抗震设计的特点,基于横观各向同性弹塑性模型,研制了考虑层状土体各向异性和施工开挖效应,适合于盾构隧道动力计算的各向异性弹塑性动力有限元程序。在程序中对于不同的材料采用了不同的本构关系和不同的单元形式,并采用了关联流动法则和多种屈服准则,可同时进行各向异性土体与地下结构的二维平面应力、平面应变和轴对称问题的静力和动力数值分析。最后利用所研制的程序进行了上海地铁二号线石门一路站附近区间隧道在不同超越概率地震动输入下的隧道反应计算。结果表明,在层状土体-地铁区间隧道的抗震分析中考虑土体各向异性的影响是必要的,所提出的计算模型是可行的。     

无限域横观各向同性岩体介质中折线裂纹边界元方法

岩体材料的各向异性导致其中的光滑裂纹在扩展后会变为折线裂纹。对于折线裂纹,其折线处裂纹面外法线方向不惟一,不能用连续单元离散。为此,在光滑裂纹问题方法的基础上,引入4种新的不连续单元来离散折线处裂纹面,建立了适用于折线裂纹问题的对偶边界元方法,该方法基于横观各向同性基本解。算例验证表明,该方法具有较好的精度。最后用该方法分析了横观各向同性岩体中的折线裂纹,得到了该类裂纹的应力强度因子。当裂纹面上作用法向均布力,横观各向同性岩体介质中的矩形光滑裂纹发生弯折时,折线两侧的裂纹面在张开时存在抑制效应,从而导致折线裂纹裂尖应力强度因子小于原光滑裂纹。同时还发现,随裂纹面的逐渐弯折,其裂纹面对各向同性面的倾角发生变化,因此,其裂尖断裂特性还受到岩体各向异性的影响。     

柱状节理坝基岩体三维各向异性数值分析

在分析柱状节理不连续性和各向异性的工程特性基础上,对某拱坝及坝基的柱状节理岩体进行三维数值模拟分析。考虑河谷初始地应力场和岩体各向异性力学性质,采用各向同性本构和各向异性弹塑性本构进行对比分析。计算表明,坝基柱状节理岩体各向异性特性对工程整体应力、变形的影响比较明显,采用各向异性弹塑性模型计算更能反映柱状节理的特殊性。其结果为该坝区岩体工程各向异性力学分析提供了示例。     

横观各向同性岩体中圆形巷道反分析的惟一性

惟一性研究是位移反分析的基础工作之一。推导了横观各向同性岩体中圆形巷道的位移解析解,利用参数可辨识条件对横观各向同性岩体中圆形巷道位移反分析的惟一性进行了探讨。结果表明,无论量测多少个点的位移也不能惟一地反演出所有6个参数;必须至少已知3个参数时,才有可能惟一地反演其他参数;2个地应力分量是否相等对反分析结果有明显的影响;2个地应力分量的可辨识性最好,各向同性面上的弹性模量和泊松比次之,垂直各向同性面方向的弹性模量和泊松比最差。