利用背景噪声进行地震成像的新方法

利用背景地震噪声提取格林函数，即通过一对台站记录的地震噪声进行互相关计算，来近似台站间的格林函数，并进一步通过地震成像获取地下结构的认识，成为最近跨学科研究的热点问题之一.本文首先叙述了其发展背景及过程，然后分别从四个角度，简要介绍了目前对其物理原理的解释.如果该方法能得到进一步的发展，将不再依赖于天然地震以及人工地震，而仅仅利用地震台站记录的地震噪声便可获取高分辨率的地下成像.该方法将会大大促进地震学的发展.     

一种新的地震危险性表达方法研究

为了寻找一种精度较高、超越概率范围较广并且便于应用的方式来表达地震危险性,本文回顾了当前常用的几种地震危险性表达方法,提出基于一个新函数来拟合地震危险性曲线的“特征系数法”,并使用《中国地震动参数区划图》（GB 18306—2015）的基础数据对该函数的拟合效果进行了验证。结果表明,新函数与地震危险性曲线拟合良好,与极值函数相比有明显的提升,能够充分地表达一个场点的地震危险性。另外,本文结果还显示该函数中表征曲线形状的参数k （文中称为特征系数）与场点面临的地震环境有关,k值较低的场点危险性贡献基本来自近场,而k值较高的场点中远距离的贡献是不能忽视的。      

On the 650-km seismic discontinuity

The pressure-temperature conditions and the variations of both density and bulk sound velocity in the vicinity of the 650-km discontinuity have been compared with those calculated for the phase transitions in both the olivine and the pyroxene-garnet components of the mantle material. These studies suggest that the mantle below about 650 km is composed primarily of perovskite phase, as distinct from the olivine-rich upper mantle. Thus, the “650-km” discontinuity is not likely to be associated with any of the equilibrium phase boundaries observed in olivine, pyroxene, and garnet, and is proposed instead to be a chemical change. It is suggested that the following factors may be responsible for chemical separation: the pyroxene-garnet component transforms to much denser phases possessing the ilmenite and perovskite structures before the breakdown of the spinel phase into a mixture of perovskite plus rocksalt phases. The perovskite phase is also much denser than the rocksalt phase and the two phases may not form a gravitationally stable mixture. Thus, the denser phases may tend to sink to or stay at the deep part of the mantle, causing chemical separation. Possible separation processes are discussed and the supporting observations are presented.     

Station corrections for the Kaapvaal seismic network: Statistical properties and relation to lithospheric structure

Station corrections for body wave travel times are required to compensate for lateral variations in the crust and uppermost mantle in the analysis of seismic travel times that are used to determine deep Earth structure by various methods, including tomography. Station corrections to be applied to P wave arrival times from teleseismic earthquakes recorded by the Kaapvaal seismic network were estimated by five different methods: (1) averaging, (2) computing the median, and (3) weighted averaging of residuals; (4) least-squares regression, and (5) weighted least-squares regression. The corrections display variations that are related to the tectonic features of southern Africa inferred from surface geology, clearly delineating the southern and central areas of both the Kaapvaal and Zimbabwe cratons as regions of early arrivals, and the area around the Bushveld complex by later arrivals. Use of a simple ray method for generating synthetic station corrections suggests that lateral variations in the top 230 km of the Earth can explain the observed pattern of variations in station corrections. A satisfactory way of compensating for the biasing effects of outliers in the individual estimates of station corrections is through adaptation of a method originally developed by Jeffreys, which involves ascribing weights to the observations that reduce the standard deviation on a single estimate of a station correction from 0.123 to 0.096 s. Methods (2), (3) and (5) avoid serious bias by outliers, although methods (3) and (5) are preferred, because they also provide information on the causes of outliers. The presence of some outliers cannot be explained by errors in the measurement process, but must be caused by timing errors at the stations during recording, and/or errors introduced during the process of constructing the archived data files from the field data.     

Implication of seismic noise for determining the structure of the upper Earth Rock Mass

An assertion that the cross-correlation function of seismic noise, considered as a result of the superposition of the surface waves, excited by the sources, randomly distributed over the Earth’s surface, determines the Green function of the surface wave is verified by numerical modeling. The maximum wave periods, for which this assertion is correct, are estimated and the errors in determination of the phase and group velocity of the surface waves are evaluated. The procedure for the determination of the correlation function and estimation from it of the group velocity are tested thoroughly based on the example of the pair of the BJT and TLY stations in Asia. This procedure is used for obtaining the group and phase velocities of the Rayleigh waves on the traces between the OBN-ARU and PUL-ARU stations. The velocity sections of the transverse waves are built based on the dispersion curves of the phase and group velocities of averages along these traces. The region of the lowered velocity in the upper mantle at depths of 150–300 km is revealed on both traces. From the analysis of correlation functions, which are subjected to narrow-band filtering, it is shown that the frequency composition of noise varies from the East and from the West from the profiles between the stations: in the East (Siberia) the noise has an appreciably lower-frequency than in the West (Western Europe).     

A new method of migration imaging for MT data

Introduction The migrationof electromagnetic field is like that of elastic wave in seismic prospecting, butthere are some differences between them (Lee et al, 1989). The elastic wave field satisfies thegeneral wave equation, while the electromagnetic field in our study satisfies the diffusion equation.Both the electromagnetic field and the elastic wave field have a similar mechanism. The imagingfunction named map can characterize the interface of different mediums (MA, 1989), that is …     

A lithospheric seismic profile along the axis of the Alps, 1975

Summary From 8 to 20 September 1975 refraction seismic measurements were carried out in close European cooperation on a long range profile along the strike of the Alps between France and Hungary. The execution and first results of the Alpine Longitudinal Profile 1975 are presented in this paper, which is the first of a series. 20 shots from 9 different shotpoints were recorded by 193 mobile stations along a main line of a length of 850 km as well as on a number of fans and additional shorter profiles. The recordings were subsequently digitized and a number of computer generated record sections are presented to illustrate the quality of the data. First results are given in the form of a simple crustal cross section along the main profile and of two velocity depth functions, which indicate a substantial difference in type between the westernmost part and the eastern part of the profile.G. Angenheister, J. Ansorge, K. Aric, D. Bamford, R. Cassinis, H. Gebrande, I. Guerra, R. Gutdeutsch, W. Kaminski, R. King, H. Miller, C. Morelli, St. Mueller, R. Nicolich, G. perrier, K. Posgay, C. Prodehl, S. Scarascia, E. Schmedes, P. Steinhauser, F. Thouvenot.     

一种新的针对反射波的逆时偏移真振幅成像条件

真振幅成像是一种代表性的定量估计模型参数扰动高波数部分的地震波成像方法.经典的真振幅成像方法在高频近似和理想照明假设条件下求取显式对角Hessian逆矩阵作为偏移振幅加权算子,用以校正波传播过程中的几何扩散效应,得到模型参数扰动的带限估计.真振幅保真成像方法在利用逆时偏移（RTM）框架实现时会产生低波数噪声,影响对高波数参数估计的精度.本文给出了一种新的基于RTM框架的真振幅保真成像条件,该成像条件针对反射波数据,在高频近似下散射模式对应正问题及Bayes反问题框架下导出.与传统基于高频渐进反演的波动方程成像方法类似,利用本文提出RTM成像条件能够保证计算结果与高频近似下反演结果的一致性.同时,利用本文提出RTM真振幅成像条件能够在成像过程中自动保真的消除传统真振幅RTM算法中存在低波数噪声,模型数值实验结果验证了本文方法的正确性和有效性.     

一种新的地震子波提取与层速度反演方法

粒子群优化算法是近十年发展起来的一种基于群智能的非线性全局最优化新方法.本文详细介绍了粒子群优化算法的基本原理,并将其应用到子波提取与层速度反演中.通过模拟数值算例,从不同角度研究了粒子群优化算法的可行性及其效率问题.试算结果表明,粒子群优化算法在不同分辨率、不同信噪比、不同相位子波合成的地震记录反演中效果明显.     

IMS体系抗震性能研究的Push-over方法

应用三维有限元实体单元计算模型,通过Push-over方法对云南傣族新民居IMS体系结构进行了地震破坏历程分析,并与模型振动台试验结果进行了比较。三维实体单元Push-over方法结果能较好的反映结构破坏表象,对结构抗震性能研究具有很好的辅助作用。     

一种提高地震勘探解释精度、分辨率的新方法-DDM方法

DDM(Depth—Density Method)方法是将地震时间剖面直接反演成深度、密度剖面的一种地震资料处理方法,在几千米的勘探深度范围内,该方法可以使地震勘探的精度由相对精度用绝对精度来表示;使分辨率由十几米提高到2～3米;从深度、密度剖面中,地质人员可以解释出地层的岩性,时代、砂体的厚度、沉积旋回、沉积粒序,几米、十几米的断层和褶皱,本文主要简介该方法的工作原理、地质解释的原则和在生产中的应用实例。     

The pyroxene-garnet transformation and its implication for the 200-km seismic discontinuity

The P-T conditions and the density (and velocity) change associated with the reaction from an alumina-rich pyroxene to an alumina-poor pyroxene plus garnet, or the lower-pressure boundary for the two-phase region of the pyroxene-garnet transition, is thoroughly examined. The 200-km seismic discontinuity has not been satisfactorily explained previously, and it is proposed here that it may be related to the above phase transformation.     

地幔岩石圈热结构差异与中国大陆岩石圈均衡分析

观测表明,大陆地区存在很多传统均衡模型无法解释的现象,其根本原因在于传统均衡理论中没有考虑地幔岩石圈部分由于热结构差异导致密度差异的影响.本文基于岩石圈尺度的质量平衡模型研究了中国大陆20个构造单元地壳及地幔岩石圈对地形海拔的贡献,以及各块体的均衡状态.计算结果表明,在一些地区,如塔里木盆地、北山和柴达木盆地,尽管岩石圈均衡模型和Airy模型得到了一致的海拔值,但岩石圈均衡模型更能体现均衡过程的物理本质;除青藏高原造山带外的多数块体,岩石圈均衡模型的计算结果更接近观测海拔和地表垂直运动状态;总体上,考虑地幔岩石圈热结构影响后,中国大陆各地区的均衡结果普遍优于传统的均衡模型.通过对均衡状态分析,我们得到以下主要结论：(1)构造稳定地区均衡程度较高;(2) 青藏高原及周边造山带现今地壳运动主要为区域构造过程及深部动力学过程所控制,均衡调整过程不是主要控制因素;(3) 现今地壳垂直运动比较明显的块体处于均衡调整阶段,地表垂直运动的大小反映了该区所受的均衡力作用的程度;(4) 构造稳定地区基于岩石层均衡的理论计算海拔与观测海拔之差值和现今地壳垂直运动速率有较好的相关性,据此我们可以通过均衡分析研究构造块体的运动趋势和动力学性质;(5) 地幔热结构对现今地形、海拔及地壳垂直运动有显著影响,在处理均衡问题时,地幔岩石圈热结构是我们必须考虑的重要因素.     

New imaging method for seismic reflection wave and its theoretical basis

Some new imaging formulas for seismic reflection wave and their theoretical basis are given. Phenomena of wave propagation should be characterized by instantaneous spectrum and expressed by complex function of three variables (time, space and frequency) in mathematics. Various physical parameters of medium are also complex functions of two variables (space and frequency). The relationship between reflection coefficient of medium and spectrum of reflected wave is given. Multi-reflection and filter of formations are considered in inversion formulas. Problems in classical convolution model and wave equation are illustrated. All these inversion formulas can be used to image underground medium by wavelet transform and method of “3-basic colors”. Different colors mean different media.     

A discontinuous Galerkin method for seismic soil–structure interaction analysis in the time domain

A technique for modeling transient wave propagation in unbounded media is extended and applied to seismic soil–structure interaction analysis in the time domain. The technique, based on the discontinuous Galerkin method, requires lower computational cost and less storage than the boundary element method, and the time‐stepping scheme resulting from Newmark's method in conjunction with the technique is unconditionally stable, allowing for efficient and robust time‐domain computations. To extend the technique to cases characterized by seismic excitation, the free‐field motion is used to compute effective forces, which are introduced on the boundary of the computational domain containing the structure and the soil in the vicinity of the structure. A numerical example on a dam–foundation system subjected to seismic excitation demonstrates the performance of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

A magnetic method for determining the geometry of hydraulic fractures

Summary We propose a method that may be used to determine the spatial orientation of the fracture plane developed during hydraulic fracture. In the method, magnetic particles are injected into the crack with the fracturing fluid so as to generate a sheet of magnetized material. Since the magnetization of a body with extreme dimension ratios, such as a crack, exceeds that of an equidimensional body and since this magnetization is sensitive both to orientation and geometry, this could be used to obtain information about the crack. By measuring the vertical and horizontal components of the magnetic field and field gradients at the earth's surface surrounding the injection well with superconducting magnetometers having 10^–4 gamma sensitivity and also by measuring field direction within the well itself, it should be possible to calculate the orientation and perhaps infer the approximate geometry of the fracture surface. Experiments on electric field potential operated in conjunction with this experiment could further constrain estimates of shape and orientation.     

A replacement for the srss method in seismic analysis

It is well-known that the application of the Square-Root-of-Sum-of-Squares (SRSS) method in seismic analysis for combining modal maxima can cause significant errors. Nevertheless, this method continues to be used by the profession for significant buildings. The purpose of this note is to present an improved technique to be used in place of the SRSS method in seismic analysis. A Complete Quadratic Combination (CQC) method is proposed which reduces errors in modal combination in all examples studied. The CQC method degenerates into the SRSS method for systems with well-spaced natural frequencies. Since the CQC method only involves a small increase in numerical effort, it is recommended that the new approach be used as a replacement for the SRSS method in all response spectrum calculations.     

A new simple method to substantially increase the seismic stability of reinforced soil structures

A preloading and prestressing (PLPS) method has been proposed to substantially decrease the transient and residual vertical compression of geosynthetic-reinforced soil (GRS) structures subjected to long-term traffic load. It is shown that by using a newly developed device (called the ratchet system) in addition to the PLPS procedure, the seismic stability of PLPS GRS structures becomes very high. The ratchet system can not only maintain high prestress when the backfill tends to contract but also prevent the expansion of the backfill, both effectively restraining the shear and bending deformation of the structure subjected to seismic load.     

基于误差准则和循环迭代的时移地震匹配滤波方法

分析了用于时移地震互均衡处理中的传统匹配滤波算法适用范围,结果表明该方法在滤波算子的反算子为最小相位时效果最好.对该方法难以解决的混合相位问题,推导出通用公式,并提出基于最小平方误差准则和循环迭代的求解方法.理论和实际数据都验证了该方法较已有方法的优越性.     

地震干涉技术被动源地震成像

被动源地震成像是基于地下反射波响应和透射波响应之间的关系,通过在地表接收由地下非人工震源发出的透射波响应,利用互相关运算来合成反射波响应的方法.合成的反射波响应,又称虚炮集记录等价于地表地震剖面记录.本文研究了由地下随机分布的被动源发出的透射波记录来合成虚炮集记录的基本原理,推导了反射波记录和透射波记录的关系公式;并通过数值模拟方法,分析了被动源成像的可行性.结果表明:利用被动源地震信息进行成像是可行的.