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.     

Semi-Analytic Modelling of Subsidence

This paper presents a forward model for subsidence prediction caused by extraction of hydrocarbons. The model uses combinations of analytic solutions to the visco-elastic equations, which approximate the boundary conditions. There are only a few unknown parameters to be estimated, and, consequently, calculations are very fast. The semi-analytic model is applicable to a uniform and layer-cake stratigraphy, with visco-elastic parameters changing per layer, and an arbitrary depletion pattern. By its capabilities to handle a multi-layered visco-elastic subsurface, the semi-analytic model fills the gap between the analytic single-layered elastic models available to date and the more elaborate numerical, e.g. finite element, models.     

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

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

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

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

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

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