多震相菲涅耳体射线走时同时反演成像

高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线附近一定范围(称其为第一菲涅耳体)内的速度异常分布有关.鉴于此,本文提出了计算多震相地震波菲涅耳体有限频射线的方法,并定义了走时敏感核函数,同时给出了利用多震相菲涅耳体有限频射线进行速度模型和反射界面同时反演成像的公式.利用多震相走时资料,使用传统射线层析成像方法与有限频射线层析成像方法进行了速度和界面的同时反演成像.结果表明,当射线密度较小时,无论是对速度模型的重建还是对反射界面几何形状的更新,有限频射线层析成像方法均优于传统射线层析成像方法, 而变频有限频射线层析成像则是实际地震层析成像的首选反演算法.      

Simultaneous Inversion for Velocity and Reflector Geometry Using Multi-phase Fresnel Volume Rays

Traditional ray tomography methods based on the high frequency assumption are sometimes unable to obtain a high resolution tomographic picture due to a deficient coverage of ray paths in real applications, especially for low velocity anomalous regions. In contrast, finite-frequency ray theory is more suitable for handling real seismic propagation problems because the travel time depends not only on the velocity distribution along a central ray (or traditional geometric ray), but also on the velocity values within a region (referred to as the first Fresnel Volume) which incorporates the central ray. In this study, we develop an algorithm to calculate multi-phase Fresnel Volume finite-frequency rays, and then present an inversion method to simultaneous invert for both velocity and reflector geometry by using these multi-phase Fresnel Volume finite-frequency rays. Using synthetic data examples, we compare the reconstructions of the velocity field and the reflector orientation using the Fresnel Volume ray tomographic methods and the traditional ray tomography approach. Results show that the former is advantageous over the latter, especially when the ray density is relatively low. An additional benefit of the Fresnel Volume finite-frequency ray tomographic method is that it can start with a low frequency to capture the coarse velocity structure, thereby mitigating the local minimum trapping problem, and then be tuned to a high frequency for delineating the fine velocity structure.     

基于网格走时计算的弯曲界面下菲涅尔体研究

为了研究弯曲界面曲率变化对分辨率的影响,首先推导了垂直入射下来自弯曲界面的反射波和透射波界面菲涅尔带近似的解析公式,证明了公式中曲率为零恰好对应已经被推导的平界面菲涅尔带的解析公式,然后给出了利用网格走时计算方法计算弯曲界面下反射波和透射波菲涅尔体的数值实现策略,这一实现策略同时保证了网格走时计算的精度和菲涅尔体计算的精度,最后对比了不同弯曲界面（不同曲率）下的菲涅尔体相对于平界面（曲率为零）下菲涅尔体的变化.研究结果表明,界面下高速时,向斜弯曲造成菲涅尔体在界面附近变宽,使得分辨率降低;背斜弯曲造成菲涅尔体在界面附近变窄,使得分辨率提高.并且向斜弯曲对分辨率的影响程度要明显大于背斜弯曲.而界面下低速时,结论正好相反.     

Discussion of the Features of Synthetic Amplitude-Distance Curves of Teleseismic P Waves for Global Earth Models

The reason why the synthetic amplitude distance curves of P-waves for models IASP91 and PREM, observed for shallow sources up to epicentral distance of 28°, oscillate is illustrated by means of synthetic seismograms. Furthermore the position of the beginning of the diffraction of the P-wave at the CMB, depending on the prevailing signal period, is discussed in connection with the extension of Fresnel volumes.     

Analytical expression for fresnel volumes and interface fresnel zones of seismic body waves. Part 1: Direct and unconverted reflected waves

Summary Fresnel volumes, plane-sectional Fresnel zones and interface Fresnel zones of direct and unconverted reflected waves are studied. Exact analytical expressions for various parameters of Fresnel volumes and Fresnel zones are derived and discussed. Among others, these expressions are related to semi-axes of Fresnel zones, to overshooting and penetration distances, and to the off-ray shifts of Fresnel zones.     

初至波相位走时层析

近地表速度结构通常是利用射线走时层析或菲涅尔体走时层析等反演方法得到的,但它们的目标函数仍利用射线走时残差构建,导致反演精度不高.为此,本文提出了基于散射积分算法的初至波相位走时层析成像方法.该方法的核心是:(1)提出了依赖于频率的相位走时概念;(2)利用依赖于频率的相位走时信息,而非单一的无限频率射线走时;(3)发展了一种改进的相位展开方法,即通过监测相位不连续性和2π周期判定来消除相位折叠现象;(4)考虑了地震波传播的有限频特征,即基于波动理论而非传统的射线路径或有限空间的菲涅尔体构建核函数.通过利用Overthrust模型的数值实验及与传统射线走时层析和菲涅尔体走时层析的对比表明:本文提出的方法是一种有效的初至波走时反演方法.同时,基于Overthrust模型的数值试验还证明了下列结论,即通过挖掘更多的走时信息的确可以获得更高的反演精度和分辨率.     

Analytical expressions for Fresnel volumes and interface fresnel zones of seismic body waves. Part 2: Transmitted and converted waves. Head waves

Summary Fresnel volumes and interface Fresnel zones of transmitted and head waves are studied. The relation derived for transmitted waves may also be used for converted reflected waves. Considerable attention is devoted to the penetration of Fresnel volumes across structural interfaces, particularly for head waves.     

Wave equation tomographic velocity inversion method based on the Born/Rytov approximation

This paper discusses Born/Rytov approximation tomographic velocity inversion methods constrained by the Fresnel zone. Calculations of the sensitivity kernel function and traveltime residuals are critical in tomographic velocity inversion. Based on the Born/Rytov approximation of the frequency-domain wave equation, we derive the traveltime sensitivity kernels of the wave equation on the band-limited wave field and simultaneously obtain the traveltime residuals based on the Rytov approximation. In contrast to single-ray tomography, the modified velocity inversion method improves the inversion stability. Tests of the near-surface velocity model and field data prove that the proposed method has higher accuracy and Computational efficiency than ray theory tomography and full waveform inversion methods.     

Lithospheric structure of Southeast Anatolia from joint inversion of local and teleseismic data

A joint tomographic inversion of local and teleseismic arrival times recorded at 41 seismic stations in southeast Anatolia is conducted to study the 3-D lithospheric velocity structure and its relation to the prevailing tectonic processes. A total of 21300 arrivals from local and teleseismic events are used in the final inversion. The tomographic model reveals prominent lower crustal/uppermost mantle low-velocity anomalies. High-velocity zones are imaged in the western part of the study area. The background seismic activity occurs mainly at the low-velocity areas and to a lesser extent in some high-velocity zones. Large crustal earthquakes occur in average velocity zones, but not in high-velocity areas that can resist stress. Results of the checkerboard resolution test indicate the reliability of the obtained images; while the large hit counts at most depth slices denote reasonable ray-path coverage for most parts of the study area. The obtained velocity anomalies are generally consistent with many previous geophysical measurements and give much deeper understanding of the current seismotectonic processes occurring in the region.     

Potential and limits of double-difference tomographic methods

Double-difference tomographic methods use directly accurate time delays computed between similar signals. Such methods are designed to image very heterogeneous media, such as volcanoes or fault zones. In seismological applications, similar signals are recorded at a given station from earthquakes sharing similar and close-by sources. In seismic exploration experiments, similar signals are often recorded at neighbouring receivers. After a brief presentation of the tomographic algorithm used, a seismological application is summarized. The potential and limits of double-difference tomographic methods are explored using various numerical experiments. They show that two effects are competing in double-difference tomography: (i) the degradation of the stability of the inversion due to the geometrical proximity of the rays used in the differentiation and (ii) the decrease in modelling error, which allows improving the stability of the inversion and using smaller quantities of a priori information when data are sufficiently accurate. The best resolution is obtained for an optimal value of the inter-source or inter-receiver distance. For optimal values of these distances and a priori information, tomography using traveltime differences provides significantly better resolved results than using traveltimes.     

二维菲涅耳带共反射面元叠加方法研究

大量研究证明CRS叠加能提高地震勘探的信噪比,但是能否提高分辨率的关键在于如何确定CRS叠加孔径.本文详细探讨了地震波反射过程中菲涅耳带的特征,认为起伏地形下菲涅耳带可以采用椭圆予以近似,在此基础上提出了一种通过菲涅耳带来确定CRS叠加孔径的方法,并应用于泌阳凹陷陡坡带的地震剖面.结果表明,由于菲涅耳带确定的叠加范围使地震信号的能量达到最佳,相对于CMP叠加,菲涅耳带CRS叠加同时提高了地震资料的信噪比和分辨率,特别显示了中深部较弱的地震信息,而常规的CRS叠加则只在于提高资料的信噪比和改善浅部的地震信息.     

复杂介质小波多尺度井间地震层析成像方法研究

复杂介质井间地震层析成像是一个很复杂的非线性反 演问题，常规的线性化反演方法无法得到好的解. 采用基于图形的弯曲射线追踪方法， 并将小波多尺度思想引入到井间层析成像，提出了小波多尺度井间地震层析成像方法，很好 地解决了非线性成像的难题，提高了图像的质量和分辨率. 物理模型实验结果表明， 该方法适合于复杂介质成像，并具有良好的实用性和效果.     

菲涅耳体地震层析成像分辨率研究

层析成像分辨率的研究,不仅可以帮助分析层析方法的反演能力,评价层析反演的效果,还可以帮助指导层析参数设置,优化观测系统设计等.本文对比研究了前人提出的两种菲涅耳体层析成像分辨率的计算方法,并针对其存在的问题进行了优化.文中通过对二维理论模型的定量计算,总结了菲涅耳体地震层析成像分辨率的一些规律,并将其与射线层析的分辨率...     

THE RESOLUTION POSSIBLE IN IMAGING WITH DIFFRACTED SEISMIC WAVES1

The determination of the vertical and lateral extent of discontinuities is an important aspect of interpreting seismic reflection data. The Common Fault Point (CFP) stacking method appears to be promising in imaging discontinuities in acoustic impedance by making use of diffracted energy from a spatial array of receivers. The problems of vertical and lateral resolution in the method are most important when carrying out an interpretation. Source signature, subsurface velocities and the depth of the discontinuity are the most important parameters affecting the resolution. We use, for a perfectly coherent source, the first derivative of the Gaussian function which is an antisymmetric band-limited wavelet. Rayleigh's, Ricker's and Widess' criteria are also applicable to this wavelet. The limits of vertical and lateral resolution are illustrated by using a step fault and a dike model respectively. The vertical resolution of the CFP method is found to be of the order of λ/16 which is half the theoretically predicted value for a single receiver. The lateral resolution is still limited by the size of the Fresnel zone which depends upon the velocity, two-way time and the dominant frequency of the wavelet. The resolution limits of the CFP method are compared with that of the CDP method, prestack migration and post-stack migration. Obtaining high resolution with real data is limited by the extent to which it is possible to generate a coherent source or to simulate one during computer processing with before stack seismic data. The CFP method is an artificial intelligence approach to imaging diffracting points as it localizes parts of the structure that scatter acoustic waves.     

The European Upper Mantle as Seen by Surface Waves

We derive a global, three-dimensional tomographic model of horizontally and vertically polarized shear velocities in the upper mantle. The model is based on a recently updated global database of Love- and Rayleigh-wave fundamental-mode phase-anomaly observations, with a good global coverage and a particularly dense coverage over Europe and the Mediterranean basin (broadband stations from the Swiss and German seismic networks). The model parameterization is accordingly finer within this region than over the rest of the globe. The large-scale, global structure of our model is very well correlated with that of earlier shear-velocity tomography models, based both on body- and surface-wave observations. At the regional scale, within the region of interest, correlation is complicated by the different resolution limits associated to different databases (surface waves, compressional waves, shear waves), and, accordingly, to different models; while a certain agreement appears to exist for what concerns the grand tectonic features in the area, heterogeneities of smaller scale are less robustly determined. Our new model is only one step towards the identification of a consensus model of European/Mediterranean upper-mantle structure: on the basis of the findings discussed here, we expect that important improvements will soon result from the combination, in new tomographic inversions, of fundamental-mode phase-anomaly data like ours with observations of surface-wave overtones, of body-wave travel times, of ambient “noise”, and by accounting for an a-priori model of crustal structure more highly resolved than the one employed here.     

Tomographic Imaging of the Seismic Structure Beneath the East Anatolian Plateau, Eastern Turkey

The high level of seismic activity in eastern Turkey is thought to be mainly associated with the continuing collision of the Arabian and Eurasian tectonic plates. The determination of a detailed three-dimensional (3D) structure is crucial for a better understanding of this on-going collision or subduction process; therefore, a body wave tomographic inversion technique was performed on the region. The tomographic inversion used high quality arrival times from earthquakes occurring in the region from 1999 to 2001 recorded by a temporary 29 station broadband IRIS-PASSCAL array operated by research groups from the Universities of Bo?azi?i (Turkey) and Cornell (USA). The data was inverted and consisted of 3,114 P- and 2,298 S-wave arrival times from 252 local events with magnitudes (M _D) ranging from 2.5 to 4.8. The stability and resolution of the results were qualitatively assessed by two synthetic tests: a spike test and checkerboard resolution test and it was found that the models were well resolved for most parts of the imaged domain. The tomographic inversion results reveal significant lateral heterogeneities in the study area to a depth of ~20?km. The P- and S-wave velocity models are consistent with each other and provide evidence for marked heterogeneities in the upper crustal structure beneath eastern Turkey. One of the most important features in the acquired tomographic images is the high velocity anomalies, which are generally parallel to the main tectonic units in the region, existing at shallow depths. This may relate to the existence of ophiolitic units at shallow depths. The other feature is that low velocities are widely dispersed through the 3D structure beneath the region at deeper crustal depths. This feature can be an indicator of the mantle upwelling or support the hypothesis that the Anatolian Plateau is underlain by a partially molten uppermost mantle.     

A numerical comparison of time and frequency‐domain marine electromagnetic methods for hydrocarbon exploration in shallow water

In shallow water the frequency domain controlled source electromagnetic method is subject to airwave saturation that strongly limits the sensitivity to resistive hydrocarbon targets at depth. It has been suggested that time‐domain CSEM may offer an improved sensitivity and resolution of these deep targets in the presence of the airwave. In order to examine and test these claims, this work presents a side‐by‐side investigation of both methods with a main focus on practical considerations, and how these effect the resolution of a hydrocarbon reservoir. Synthetic noisy data for both time‐domain and frequency domain methods are simulated using a realistic frequency dependent noise model and frequency dependent scaling for representative source waveforms. The synthetic data studied here include the frequency domain response from a compact broadband waveform, the time‐domain step‐response from a low‐frequency square wave and the time‐domain impulse response obtained from pseudo‐random binary sequences. These data are used in a systematic resolution study of each method as a function of water‐depth, relative noise and stacking length. The results indicate that the broadband frequency domain data have the best resolution for a given stacking time, whereas the time‐domain data require prohibitively longer stacking times to achieve similar resolution.     

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.     

Optimization of Cell Parameterizations for Tomographic Inverse Problems

—?We develop algorithms for the construction of irregular cell (block) models for parameterization of tomographic inverse problems. The forward problem is defined on a regular basic grid of non-overlapping cells. The basic cells are used as building blocks for construction of non-overlapping irregular cells. The construction algorithms are not computationally intensive and not particularly complex, and, in general, allow for grid optimization where cell size is determined from scalar functions, e.g., measures of model sampling or a priori estimates of model resolution. The link between a particular cell j in the regular basic grid and its host cell k in the irregular grid is provided by a pointer array which implicitly defines the irregular cell model. The complex geometrical aspects of irregular cell models are not needed in the forward or in the inverse problem. The matrix system of tomographic equations is computed once on the regular basic cell model. After grid construction, the basic matrix equation is mapped using the pointer array on a new matrix equation in which the model vector relates directly to cells in the irregular model. Next, the mapped system can be solved on the irregular grid. This approach avoids forward computation on the complex geometry of irregular grids. Generally, grid optimization can aim at reducing the number of model parameters in volumes poorly sampled by the data while elsewhere retaining the power to resolve the smallest scales warranted by the data. Unnecessary overparameterization of the model space can be avoided and grid construction can aim at improving the conditioning of the inverse problem. We present simple theory and optimization algorithms in the context of seismic tomography and apply the methods to Rayleigh-wave group velocity inversion and global travel-time tomography.     

Estimating Crustal Heterogeneity from Double-difference Tomography

Seismic velocity parameters in limited, but heterogeneous volumes can be inferred using a double-difference tomographic algorithm, but to obtain meaningful results accuracy must be maintained at every step of the computation. MONTEILLER et al. (2005) have devised a double-difference tomographic algorithm that takes full advantage of the accuracy of cross-spectral time-delays of large correlated event sets. This algorithm performs an accurate computation of theoretical travel-time delays in heterogeneous media and applies a suitable inversion scheme based on optimization theory. When applied to Kilauea Volcano, in Hawaii, the double-difference tomography approach shows significant and coherent changes to the velocity model in the well-resolved volumes beneath the Kilauea caldera and the upper east rift. In this paper, we first compare the results obtained using MONTEILLER et al.'s algorithm with those obtained using the classic travel-time tomographic approach. Then, we evaluated the effect of using data series of different accuracies, such as handpicked arrival-time differences (``picking differences'), on the results produced by double-difference tomographic algorithms. We show that picking differences have a non-Gaussian probability density function (pdf). Using a hyperbolic secant pdf instead of a Gaussian pdf allows improvement of the double-difference tomographic result when using picking difference data. We completed our study by investigating the use of spatially discontinuous time-delay data.