ITERATIVE TOMOGRAPHIC METHODS TO LOCATE SEISMIC LOW-VELOCITY ANOMALIES: A MODEL STUDY1

The possibilities for reconstructing seismic velocity distributions containing low-velocity anomalies by iterative tomographic methods are examined studying numerical and analogue 2D model data. The geometrical conditions of the model series were designed to generalize the geometrical characteristics of a typical cross-hole tomographic field case. Models with high (30%) and low (8%) velocity contrasts were realized. Traveltimes of 2D ultrasonic P-waves, determined for a dense net of raypaths across each model, form the analogue data set. The numerical data consists of traveltimes calculated along straight raypaths. Additionally, a set of curved-ray traveltimes was calculated for a smoothed version of the high-contrast model. The Simultaneous Iterative Reconstruction Technique (SIRT) was chosen from the various tomographic inversion methods. The abilities of this standard procedure are studied using the low-contrast model data. The investigations concentrate on the resolving power concerning geometry and velocity, and on the effects caused by erroneous data due to noise or a finite time precision. The grid spacing and the source and receiver patterns are modified. Smoothing and slowness constraints were tested. The inversion of high-contrast analogue model data shows that curved raypaths have to be considered. Hence, a ray-tracing algorithm using velocity gradients was developed, based on the grid structure of the tomographic inversion. This algorithm is included in the SIRT-process and the improvements concerning anomaly localization, resolution and velocity reconstruction are demonstrated. Since curved-ray tomography is time-consuming compared with straight-ray SIRT, it is necessary to consider the effects of grid spacing, ray density, slowness constraints and the     

Acoustic tomography in the atmospheric surface layer

Acoustic tomography is presented as a technique for remote monitoring of meteorological quantities. This method and a special algorithm of analysis can directly produce area-averaged values of meteorological parameters. As a result consistent data will be obtained for validation of numerical atmospheric micro-scale models. Such a measuring system can complement conventional point measurements over different surfaces. The procedure of acoustic tomography uses the horizontal propagation of sound waves in the atmospheric surface layer. Therefore, to provide a general overview of sound propagation under various atmospheric conditions a two-dimensional ray-tracing model according to a modified version of Snell’s law is used. The state of the crossed atmosphere can be estimated from measurements of acoustic travel time between sources and receivers at different points. Derivation of area-averaged values of the sound speed and furthermore of air temperature results from the inversion of travel time values for all acoustic paths. Thereby, the applied straight ray two-dimensional tomographic model using SIRT (simultaneous iterative reconstruction technique) is characterised as a method with small computational requirements, satisfactory convergence and stability properties as well as simple handling, especially, during online evaluation.     

跨孔雷达全波形反演成像方法的研究

跨孔雷达全波形反演是一种使用全波形信息反演两钻孔之间地下信息的层析成像技术.常规的层析成像反演大部分采用射线追踪方法,其中基于初至时的射线追踪方法可以反演出速度剖面（介电常数剖面）,基于最大振幅的层析成像可以反演出衰减剖面（电导率剖面）.常规射线追踪方法有许多不足,究其原因是该方法仅使用了小部分的信号信息.为了进一步提高成像分辨率,本文全面推导了全波形跨孔雷达层析成像反演方法,该方法利用雷达波全幅度相位信息能够反演出地下高分辨率的介电常数和电导率图像.本文通过基于局域网的分布式并行算法,有效地解决了巨量数据正演计算问题.文中首先建立了基于单轴各向异性介质完全匹配层的时间域有限差分二维正演算法,进而通过应用包括时间维度在内的全波场信息与残场逆向传播的全波场信息乘积来计算梯度方向,通过求取以步长为自变量的目标函数的极值确定步长公式,并提出以第一次介电常数反演作为同步反演的初始模型,能够有效提高收敛速度.本文对多组模型进行成像实验,取得了较好的反演效果.     

Geostatistical traveltime tomography in elliptically anisotropic media

In geological materials, anisotropy may arise due to different mechanisms and can be found at different scales. Neglecting anisotropy in traveltime tomographic reconstruction leads to artefacts that can obscure important subsurface features. In this paper, a geostatistical tomography algorithm to invert cross‐hole traveltime data in elliptically anisotropic media is presented. The advantages of geostatistical tomography are that the solution is regularized by the covariance of the model parameters, that known model parameters can be used as constraints and fitted exactly or within a prescribed variance and that stochastic simulations can be performed to appraise the variability of the solution space. The benefits of the algorithm to image anisotropic media are illustrated by two examples using synthetic georadar data and real seismic data. The first example confirms suspected electromagnetic anisotropy in the vadose zone caused by relatively rapid water content variations with respect to wavelength at georadar frequencies. The second presents how sonic log data can be used to constrain the inversion of cross‐well seismic data and how geostatistical simulations can be used to infer parameter uncertainty. Results of both examples show that considering anisotropy yields a better fit to the data at high ray angles and reduces reconstruction artefacts.     

A hybrid fast algorithm for first arrivals tomography

A hybrid algorithm, combining Monte-Carlo optimization with simultaneous iterative reconstructive technique (SIRT) tomography, is used to invert first arrival traveltimes from seismic data for building a velocity model. Stochastic algorithms may localize a point around the global minimum of the misfit function but are not suitable for identifying the precise solution. On the other hand, a tomographic model reconstruction, based on a local linearization, will only be successful if an initial model already close to the best solution is available. To overcome these problems, in the method proposed here, a first model obtained using a classical Monte Carlo-based optimization is used as a good initial guess for starting the local search with the SIRT tomographic reconstruction. In the forward problem, the first-break times are calculated by solving the eikonal equation through a velocity model with a fast finite-difference method instead of the traditional slow ray-tracing technique. In addition, for the SIRT tomography the seismic energy from sources to receivers is propagated by applying a fast Fresnel volume approach which when combined with turning rays can handle models with both positive and negative velocity gradients. The performance of this two-step optimization scheme has been tested on synthetic and field data for building a geologically plausible velocity model.This is an efficient and fast search mechanism, which permits insertion of geophysical, geological and geodynamic a priori constraints into the grid model and ray path is completed avoided. Extension of the technique to 3D data and also to the solution of 'static correction' problems is easily feasible.     

变阻尼约束层析成像及其在VSP资料中的应用(英文)

初至波走时层析成像已经取得了广泛的应用,然而,由于观测系统的限制,射线在模型中分布不均匀,导致层析结果的分辨能力不足。变阻尼约束方法应用不均匀的先验信息来匹配不均匀的数据分布,可以减小速度模型校正量与射线覆盖程度的相关性。本文将变阻尼约束方法应用于初至波旅行时层析成像中,并将平滑约束方法加入正则化方程组中来避免单独使用变阻尼约束带来的不稳定性,利用阿尔法滤波器对反演中间迭代结果进行平滑和去噪,采用LSQR算法求解线性方程组来提高收敛速度和压制误差传递。本文应用上述层析成像算法对VSP观测系统进行速度反演,分别应用于检测板速度模型数据和实际VSP资料速度反演中,结果表咀变阻尼约束层析成像可以改善射线不均匀覆盖带来的影响,从而提高速度反演结果的质量;VSP资料检波点附近的速度反演结果可靠性高。     

顾及电离层变化的层析反演新算法

区别于以往GPS电离层层析研究主要关注迭代模型的思路,本文从两方面入手提高GPS电离层层析迭代算法的反演精度:一方面,顾及传统电离层层析迭代模型仅与对电子密度误差起放大作用的GPS射线截距权重相关的不足,提出考虑层析像素格网中的电子密度对GPS TEC的贡献建立新的迭代模型,在不同电子密度像素格网内重新分配GPS TEC实测值与其反演值之间的差距;另一方面,顾及电离层层析迭代算法中松弛因子对反演结果的影响,提出考虑电子密度变化构造新的松弛因子,抑制传播噪声对电子密度反演精度的影响.实验结果显示,相对于传统代数重构算法(ART),新方法反演的电离层电子密度剖面更接近于电离层测高仪观测的电子密度剖面,提高了电子密度反演精度.     

井间电磁场的一维、二维联合反演方法

提出了针对轴对称二维井间模型的一维、二维联合反演方法.该方法基于纵向成层背景地层的Green函数,以一维反演结果作为联合反演的迭代初始值,同时反演一维背景地层和二维异常剖面的电导率.采用递推技术计算Green函数的系数,可以很方便地同时得到Green函数对空间坐标的积分和对背景层电导率的微分,大大加快了计算雅可比矩阵〖WTHX〗M〖WTBZ〗的速度并使计算方便、准确.反演实例说明了该方法的有效性.反演实例还显示,背景层电导率较二维成像剖面的电导率收敛速度要快,测量数据的精度对二维剖面电导率成像质量的影响大于对背景层电导率分布的影响.     

海洋拖缆地震资料初至波走时层析反演

初至波走时层析反演技术作为建立近地表速度模型的重要手段,是解决陆地资料复杂静校正问题的关键技术。而折射波广泛发育的海洋地震资料,对折射波信息的关注与运用并没有得到广泛的重视。本文首次将层析反演方法应用于海洋拖缆地震数据的近海底速度模型的建立。本文方法与陆地资料层析反演的主要区别在于:①在震源信号的最小相位化处理后进行初至时间的拾取,避免了混合相位子波初至拾取不准带来的误差;②以海水深度与海水速度作为反演约束条件,减小了迭代误差。实测二维资料的层析反演结果表明,本文方法可反演出较为精确的海洋地层速度结构。      

三维地震与地面微地震联合校正方法

由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景.     

Application of tomographic inversion in studying airglow in the mesopause region

It is pointed out that observations of periodic nightglow structures give excellent information on atmospheric gravity waves in the mesosphere and lower thermosphere. The periods, the horizontal wavelengths and the phase speeds of the waves can be determined from airglow images and, using several cameras, the approximate altitude of the luminous layer can also be determined by triangulation. In this paper the possibility of applying tomographic methods for reconstructing the airglow structures is investigated using numerical simulations. A ground-based chain of cameras is assumed, two-dimensional airglow models in the vertical plane above the chain are constructed, and simulated data are calculated by integrating the models along a great number of rays with different elevation angles for each camera. After addition of random noise, these data are then inverted to obtain reconstructions of the models. A tomographic analysis package originally designed for satellite radiotomography is used in the inversion. The package is based on a formulation of stochastic inversion which allows the input of a priori information to the solver in terms of regularization variances. The reconstruction is carried out in two stages. In the first inversion, constant regularization variances are used within a wide altitude range. The results are used in determining the approximate altitude range of the airglow structures. Then, in the second inversion, constant non-zero regularization variances are used inside this region and zero variances outside it. With this method reliable reconstructions of the models are obtained. The number of cameras as well as their separations are varied in order to find out the limitations of the method.     

跨孔地震层析成像的级联方法

对厚度小于1/4波长的超薄波速干扰体进行高分辨率成像,作者曾采用走时反演和波场反演相结合的方法,体现了由低分辨率向高分辨率逐步逼近的思想.级联算法是这种思想的进一步发展,我们将具有不同分辨率的算法串联起来,以达到高分辨率成像的目的.本文介绍一种三级串联的算法,并进行了算法分析,数值计算的例子说明这种级联算法分辨率高、稳定性好,只需要地震资料而不要求其它先验信息,因此能较好地满足实际应用的要求.     

3-D multi-parameter type traveltime tomography in a spherical coordinate frame: comparison of double and triple class simultaneous inversions

It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later(primary reflected) arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm(in a spherical coordinate system) with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data(noise free and added noise) indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.     

用于速度图象重建的层析成象法

本文提出了一种用于速度图象重建的层析成象法。特别注意到地震学上的ST与医学上的CT的不尽相同之处,并处理了由此引起的困难。与速度图象重建有关的正问题的分析表明,不同频率范围的资料其分辨能力是不同的,因此,区域地震与远震资料的解释应予注意。本文有关速度图象重建的讨论不仅涉及到与成象方法有关的反演理论和方法,还给出了进行数值计算的内存量和运算量的分析。关于重建图象的质量评价,还给出了反演解的可靠性分析,为分辨特性的描述提供了一种既简便又形象直观的方案。 本文的方法与现行主要方法的比较表明:1.本文的方法允许存在速度间断面,有利于揭示不同构造区域的地壳厚度的差异;2.在常用的方法中,作为正问题假定的常速度块同反演结果解释时非常速度块之间不协调。由于我们把给定网格点的速度值的插值函数作为速度的空间函数,这一困难已被解决;3.相对于Aki等的一般方法,本文的内存量节省约一个量级、运算量减少约一半;4.相对于ART类方法本身不能给出成象结果的可靠性分析,本文的方法则具有显著的优点。 作者在分析了天然地震资料的状况之后相信,对充分利用我国地震台网的现有观测资料,本文的方法是有效的。     

Two-dimensional joint inversion of radiomagnetotelluric and direct current resistivity data

An algorithm for the two-dimensional (2D) joint inversion of radiomagnetotelluric and direct current resistivity data was developed. This algorithm can be used for the 2D inversion of apparent resistivity data sets collected by multi-electrode direct current resistivity systems for various classical electrode arrays (Wenner, Schlumberger, dipole-diplole, pole-dipole) and radiomagnetotelluric measurements jointly. We use a finite difference technique to solve the Helmoltz and Poisson equations for radiomagnetotelluric and direct current resistivity methods respectively. A regularized inversion with a smoothness constrained stabilizer was employed to invert both data sets. The radiomagnetotelluric method is not particularly sensitive when attempting to resolve near-surface resistivity blocks because it uses a limited range of frequencies. On the other hand, the direct current resistivity method can resolve these near-surface blocks with relatively greater accuracy. Initially, individual and joint inversions of synthetic radiomagnetotelluric and direct current resistivity data were compared and we demonstrated that the joint inversion result based on this synthetic data simulates the real model more accurately than the inversion results of each individual method. The developed 2D joint inversion algorithm was also applied on a field data set observed across an active fault located close to the city of Kerpen in Germany. The location and depth of this fault were successfully determined by the 2D joint inversion of the radiomagnetotelluric and direct current resistivity data. This inversion result from the field data further validated the synthetic data inversion results.     

Determination of a shallow velocity–depth model from seismic refraction data by coherence inversion1

Seismic refractions have different applications in seismic prospecting. The travel- times of refracted waves can be observed as first breaks on shot records and used for field static calculation. A new method for constructing a near-surface model from refraction events is described. It does not require event picking on prestack records and is not based on any approximation of arrival times. It consists of the maximization of the semblance coherence measure computed using shot gathers in a time window along refraction traveltimes. Time curves are generated by ray tracing through the model. The initial model for the inversion was constructed by the intercept-time method. Apparent velocities and intercept times were taken from a refraction stacked section. Such a section can be obtained by appling linea moveout corrections to common-shot records. The technique is tested successfully on synthetic and real data. An important application of the proposed method for solving the statics problem is demonstrated.     

采用射线寻迹的井间层析成像

本文采用二维井间观测走时反演地震波速,用迭代射线寻迹走时成像方法进行层析成像.作者导出三次多项式插值基函数计算离散网格内的折射率分布,还给出了具有自己特点的二步扫描法,不但加快了正演计算时间,同时方便地建立走时残差与模型修正量的关系.考虑到实际应用时,震源与接收器常不位于同一平面,为此进行了校正,并给出了校正方法,提高了重建精度. 数值的与实际应用的结果表明,本文给出的层析成像算法是快速的、稳定的、适用的.     

采用射线寻迹的井间层析成像

本文采用二维井间观测走时反演地震波速,用迭代射线寻迹走时成像方法进行层析成像.作者导出三次多项式插值基函数计算离散网格内的折射率分布,还给出了具有自己特点的二步扫描法,不但加快了正演计算时间,同时方便地建立走时残差与模型修正量的关系.考虑到实际应用时,震源与接收器常不位于同一平面,为此进行了校正,并给出了校正方法,提高了重建精度. 数值的与实际应用的结果表明,本文给出的层析成像算法是快速的、稳定的、适用的.     

Iterative tomographic analysis based on automatic refined picking

The ever-growing size of data sets for active and passive seismic imaging makes the availability of automatic procedures for rapid analysis more and more valuable. Such procedures are especially important for time-critical applications like emergency decisions or re-orienting of ongoing seismic surveys. In this paper a new, iterative scheme for 3D traveltime tomography is presented. The technique, based on a tool originally developed for earthquake data, uses cross-correlation to examine waveform similarity and to adjust arrival times on seismic sections. A preliminary set of reference arrival times is first corrected by the cross-correlation lag and then used to build an initial 3D tomographic velocity model through a standard inversion code; traveltimes calculated from this model are then taken as new reference arrivals and the process of pick adjustment is repeated. The result is a tomographic image, upgraded and refined at each iteration of the procedure. The test performed on the waveform data set recorded during the 2001 SERAPIS active seismic survey in the gulfs of Naples and Pozzuoli (Southern Italy) shows that the 3D iterative tomography scheme produces a velocity image of the structure of the Campi Flegrei caldera which is consistent with the results from previous studies, employing just a fraction of the time needed by a human analyst to identify first breaks. We believe that this technique can be effectively employed for rapid analysis of large data-sets within time-critical or time-dependent tasks and for automatic 4D tomographic investigations.