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.     

A traveltime reciprocity discrepancy in the Podvin & Lecomte time3d finite difference algorithm

We have found that the extensively used finite difference scheme time3d produces time fields which are 'asymmetric' in the sense that traveltimes computed to the right of the source are faster than traveltimes computed to the left. All finite difference schemes will, as they are approximations to the wave equation, to some extent fail to obey reciprocity perfectly. We show, however, that the errors in time3d may be significant—and unnecessarily large. An asymmetry in the point source initialization has been identified, and after correction time3d produces time fields with an improved reciprocity.     

利用三维地震在戈壁沙漠区解决复杂地质构造问题

在分析戈壁沙漠地区地球物理特征的基础上，针对性的提出了解决煤田地质构造问题的方法，包括三维地震野外施工方法的选取、施工因素的选择、资料处理手段；并以获得的实际地质成果证实三维地震在解决复杂地质构造问题的可行性。     

Seismic Hazards Assessment for Radioactive Waste Disposal Sites in Regions of Low Seismic Activity

A comparative tectonic quiescence and lack of earthquakes make the stable centres of continents attractive for siting long-term radioactive waste storage facilities. The low rates of deformation in such regions, however, make it difficult to characterize their long-term seismotectonic behaviour, leading to uncertain estimates for the very low probability hazard estimates required by society. In an attempt to overcome the deficiency of both contemporary seismicity and paleoseismic data in central Canada, we have used earthquake histories from regions with similar seismotectonic characteristics from around the world. Substituting space for time, we estimate a long-term rate per 10⁶ km² of 0.004 magnitude ≥6 earthquakes per annum, of which 33–100% might rupture to the surface.     

陆相油藏中砂泥岩互层时移地震遮掩效应及其对时移地震资料解释的影响

从差异波形特征的角度出发,分析时移地震差异波形特征与砂泥岩互层油藏中油层速度变化的关系。通过建立由5层地层组成的楔形地质模型,运用褶积公式及时移地震求差公式,得到泥岩速度为2 400 m/s和砂岩速度为2 700m/s条件下,互层中砂岩速度变化±150 m/s的时移地震响应。研究结果表明:砂泥岩互层的差异波形是一个多峰复合波,泥岩夹层的厚度是决定差异波形特征的一个主要因素,并且存在一个临界值,即四分之一主波长;如果泥岩夹层的厚度小于这一临界值,就会发生遮掩效应,即顶部油层的时移地震响应可以掩盖中、下部油层的时移地震响应。     

地震信号的相位恢复及误差分析

常见的几种地震子波有:零相位子波、最小相位子波、最大相位子波和混合相位子波。以不同类型的信号为激发子波,对不同地层模型合成地震道。仅用原始合成地震道记录的振幅谱来恢复原始记录,讨论不同方法下,由于激发子波的不同,对比分析原始记录与恢复记录的差异,及其恢复效果。     

深层地震勘探激发井深探讨

深层地震勘探野外资料采集的关键是在保证原始资料一定信噪比的基础上,同时兼顾分辨率,激发井深是决定资料品质的重要因素之一。通过对虚反射滤波作用分析,提出了充分利用虚反射能量设计最佳激发井深,增强下传地震波能量的采集方法。理论分析和试验资料表明,在高速层下4～7m的最佳激发岩性深度是深层勘探的理想井深。     

基于复小波变换的SAR图像斑点滤波方法

斑点噪声的存在使得SAR图像的应用受到较大的限制.文中首先介绍了SAR图像斑点噪声模型以及几种常用的SAR斑点噪声抑制方法;然后对小波变换斑点噪声滤波方法进行分析,提出了基于二元树复小波变换的SAR图像滤波方法.实验比较表明,二元树复小波变换抑制SAR斑点噪声效果明显,较小波变换更具有优越性.     

地震波混合阶褶积算法模拟

基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。