Wavepath traveltime tomography

The elastic-wave equation is used to construct sensitivity kernels relating perturbations in elastic parameters to traveltime deviations. Computation of the functions requires a correlation of the forward-propagating seismic wavefield with a backward propagation of the residual wavefield. The computation of the wavefields is accomplished using a finite difference algorithm and is efficiently executed on a CM-2 parallel processor. The source and receiver locations have maximum sensitivity to velocity structure. The sensitivity kernels or wavepaths are well suited for transmission traveltime inversion such as cross-borehole tomography and vertical seismic profiling. Conventional ray tomography and wavepath tomography are applied to a set of P -wave arrival times, from a cross-borehole experiment at Kesterson, California. Because the wavepaths have increased sensitivity near the source and receiver there are differences in resolution of the velocity structure. Both techniques recover the same relative variations in velocity where the coverage is adequate. The wavepath solution is more laterally continuous and the dominant variation is vertical, as is expected for the layered sediments in this region.     

Inference of core-mantle boundary topography from ISC PcP and PKP traveltimes

We invert ISC PcP and PKP absolute and differential traveltimes in an attempt to infer the long-wavelength topography of the core-mantle boundary (CMB). The data selection and processing methods are described and evaluated. These travel-time data are very noisy and the geographic distribution of the data is highly non-uniform, inhibiting reliable inference of CMB topography. Spatial averaging enhances the coherent component of the residual variance (related to heterogeneity), however, the random component of the variance is much larger than the coherent component. We show that for PcP data the coherent signal due to mantle heterogeneity overshadows that arising from the CMB, and that the effects of mantle heterogeneity are mapped into our inferred CMB solutions. The PcP data are not correlated across the spatial averaging bins and seem to have a strong bias due to small-scale structure and/or noise. The non-uniform geographic sampling of the data plays a role in the mapping of mantle heterogeneity onto the CMB. Spatial patterns of CMB models inferred from different phases do not agree. Amplitudes of seismically inferred CMB undulations vary greatly. The sensitivity of inferred CMB models to the processing, spatial averaging procedure, and inversion techniques are investigated. Topographic amplitudes increase strongly with increasing input residual variance. The power spectrum of inferred topography indicates that there are unmodelled heterogeneities that must be described with spherical harmonics of degree 6 and higher. Based on this work, we conclude that reliable inference of long-wavelength CMB topography is not likely with the current ISC data set or with a spherical harmonic expansion truncated to degree and order 6.     

Seismic recognition of local inhomogeneities

A simple method to contour local inhomogeneities using seismic data is proposed. It formalizes an approximate inversion method which is based on the interpretation of local inhomogeneities as making the differences between an actual seismic data set and a previous reference model. It uses the optimal statistical criteria of parameter estimation and recognition and the ray representation of the waves spreading. Any combination of direct, reflected and/or other types of waves may be used as the database. Inhomogeneities, having a size two times above the wavelength of the seismic waves, can be resolved. Laboratory experiments, using ultrasonic waves and analysis of data from field experiments, confirmed the theoretical results. The method can be used to search for ore bodies, kimberlite cubes, oiltraps, etc.     

阿德朗达克-西魁北克地震带地壳三维速度结构成像

根据在阿德朗达克-西魁北克地震带及其邻近地区的59台短周期数字地震仪记录的直达P波走时观测值,用地震层析成像方法反演了该地区地壳三维速度结构.结果表明,地壳上层（0-5km）在阿德朗达克穹隆山的中南部地区出现正速度异常;第二层（5-10km）及第三层（10-15km）的速度横向变化较小,介质相对比较均匀;第四层（15-25km）出现显著的速度异常,一个是位于地震带中部的正速度异常（+4％）,它与布格重力正异常一致;另一个是位于阿德朗达克穹隆山下的负速度异常（-4％）.结合已有的地质及地球物理资料进行对比,可以认为阿德朗达克是一个正在发展的穹隆上升山,但其热源前锋尚未到达地面的大陆热点.     

唐山滦县震区的三维地震透射研究──中、上地壳速度层析成像

利用布设在唐山滦县震区３０ｋｍ×４０ｋｍ范围内的由８８台数字和模拟地震仪组成的临时台阵，接收来自不同方向６个炮点激发产生的莫霍界面反射波走时，重建了台阵下１-９ｋｍ深度的Ｐ，Ｓ波速度和Ｖ_ｐ／Ｖ_ｓ图像．结果表明。中、上地壳结构存在明显的横向不均匀性．北北东走向的滦县-卢龙断裂自地表向下至少延伸到８ｋｍ深度处，并向北西方向倾斜．中、上地壳存在着近东西向的低速异常条带．这两组构造控制着该区的地震活动．     

岩样中包体构造的CT实验探测

CT技术图象重建方法可分为变换法和代数法两大类。本文讨论了几种代数重建法的特点,开展了地球物理CT技术在岩石实验中应用的探讨性工作。研究表明,实验数据的测量误差直接影响着反演结果,当测量误差比较大时,即使数据量很大,也难以提高反演结果的分辨。岩样的吸收、换能器与样品的耦合、岩样均匀性等因素对实验数据的测量精度都有着重要的影响。文中还针对实验条件进行了数值模拟,以检验本方案的可行性,并对当前流行的BPT、ART和SIRT算法及其组合算法的成象能力结合本实验的条件进行了性能的比较。结果表明:BPT方法计算简单,但所给出的图象偏粗糙,特别是在井间观测时水平方向的分辨率较差;ART方法收敛速度快,图象的反差大,反演结果的高频成分比较丰富,但边界异常大;SIRT方法收敛速度慢,图象的低频比较好,结构平滑,边界异常小,反演受测量误差的影响比较小。综合了ART和SIRT两种反演方法特点的ART-SIRT联合反演方法比较理想,发挥出了各自的长处。     

Tomographic determination of 3-D crustal structure — Joint inversion of explosion and earthquake data

In this paper, the theory and method, obtaining the tomographic determination of three-dimensional velocity structure of the crust by use of the joint inversion of explosion and earthquake data, are given. The velocity distribution of the crust is regarded as a continuous function of the spatial coordinates without parametrization of the velocity model ahead, so that the inversion solution would not be influenced by different parametrization procedures. The expressions of integration kernels, which relates the two kinds of data sets, are also given. The authors have processed the observed data in Tangshan earthquake region by the method proposed in this paper, and obtained the tomographic results of the middle and upper crust structures in this region. The comparison of these results with the result obtained only by the explosion data, has also been made.     

渭河断陷及邻近地区莫霍界面速度图象

根据陕西省及邻省地震台网的地震记录，利用地震及工矿爆破的首波走时，采用地震层析技术，重建了渭河断陷及邻近地区的莫霍界面速度图象。并对速度横向不均匀性及其与地质构造、地震活动的关系进行了探讨