Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.     

半无限空间中圆形孔洞周围SH波的散射

建立了求解在含有圆形孔洞的弹性半空间中SH波散射与圆形孔洞附近动应力集中问题的解析方法。利用SH波散射的对称性和多极坐标的方法，构造了一个可以预先满足半空间自由表面上应力自由边界条件的圆形孔洞对SH波散射的波函数。利用这一波函数，则可将该问题转化成对一个圆形孔洞散射的求解问题。该问题的解答最终又可归结为对一组无究代数方程组的求解问题，并可利用截断有限项的方法对其进行计算，最后给出了有关圆形也洞附近动应力集中问题的算例和数值结果，并讨论了波数与圆孔至自由边界距离变化对动应力集中的影响。     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

中国大陆及邻近海域的Rayleigh波群速度分布

利用Rayleigh波群速度资料反演得到中国大陆及其临近海域的(70°E-145°E,10°N-55°N)15-120s周期的群速度分布图像. 塔里木盆地在15s处清楚地显示为低速,在16-33s左右没有显示,但在36-5s显示为高速,说明塔里木盆地有较深的根. 青藏高原块体是44s至120s图像中最为突出的低速块体,南面与印度板块的分界线以及与北面的塔里木盆地、柴达木盆地的分界清晰,其块体中西部的速度低于东部. 泰国清迈附近存在一尺度为1000km左右的低速带,可能是青藏高原块体的物质向东南方向迁移造成上地幔物质上涌的结果. 南北地震带表现为强烈的速度梯度带,西面为低速,东面为高速. 中国南海的中央、日本海中央、菲律宾海表现为海洋性地壳. 菲律宾海的图像与地形及地震带完全吻合. 环绕菲律宾海及日本海存在400km左右宽的低速带,可能是岩浆活动带.     

克孜尔水库大坝的面波检测

1999年 3月 15日克孜尔地震之后对克孜尔水库坝体出现的裂缝进行了面波勘探。首先介绍了面波勘探法的基本原理 ,然后通过典型剖面检测与资料处理找出异常区、带并进行解释。试验表明面波勘探法适用于克孜尔水库大坝监测