共查询到20条相似文献，搜索用时 78 毫秒
1.
The applicability of seismic refraction profiling for the detection of velocity inversion, which is also known as a low-velocity layer (LVL), is investigated with the aid of synthetic seismogram computations for a range of models. Our computational models focus on the inherent ambiguities in the interpretation of first-arrival time delays or 'skips' in terms of LVL model parameters. The present modelling results reveal that neither the measure nor even the existence of a shadow zone and/or a time shift (skip) in first arrivals is necessarily indicative of an LVL. Besides attenuation effects, the cap-layer velocity gradient is a critical parameter, determining the termination point of the cap-layer diving wave and thus the time skip.
We suggest that shallow LVLs can be delineated more reliably by traveltime and amplitude modelling of coherent phases reflected from their top and bottom boundaries, often clearly observed in the pre- and near-critical ranges in seismogram sections of refraction profiling experiments with a close receiver spacing. We demonstrate the applicability of this approach for a field data set of a refraction profile in the West Bengal Basin, India. The inferred LVL corresponds to the Gondwana sediments underlying the higher-velocity layer of the Rajmahal Traps. This interpretation is consistent with the data from a nearby well in the region. 相似文献
We suggest that shallow LVLs can be delineated more reliably by traveltime and amplitude modelling of coherent phases reflected from their top and bottom boundaries, often clearly observed in the pre- and near-critical ranges in seismogram sections of refraction profiling experiments with a close receiver spacing. We demonstrate the applicability of this approach for a field data set of a refraction profile in the West Bengal Basin, India. The inferred LVL corresponds to the Gondwana sediments underlying the higher-velocity layer of the Rajmahal Traps. This interpretation is consistent with the data from a nearby well in the region. 相似文献
2.
Finite-Difference Modeling and Dispersion Analysis of High-Frequency Love Waves for Near-Surface Applications 总被引：3，自引：0，他引：3
Yinhe Luo Jianghai Xia Yixian Xu Chong Zeng Jiangping Liu 《Pure and Applied Geophysics》2010,167(12):1525-1536
Love-wave propagation has been a topic of interest to crustal, earthquake, and engineering seismologists for many years because it is independent of Poisson’s ratio and more sensitive to shear (S)-wave velocity changes and layer thickness changes than are Rayleigh waves. It is well known that Love-wave generation requires the existence of a low S-wave velocity layer in a multilayered earth model. In order to study numerically the propagation of Love waves in a layered earth model and dispersion characteristics for near-surface applications, we simulate high-frequency (>5 Hz) Love waves by the staggered-grid finite-difference (FD) method. The air–earth boundary (the shear stress above the free surface) is treated using the stress-imaging technique. We use a two-layer model to demonstrate the accuracy of the staggered-grid modeling scheme. We also simulate four-layer models including a low-velocity layer (LVL) or a high-velocity layer (HVL) to analyze dispersive energy characteristics for near-surface applications. Results demonstrate that: (1) the staggered-grid FD code and stress-imaging technique are suitable for treating the free-surface boundary conditions for Love-wave modeling, (2) Love-wave inversion should be treated with extra care when a LVL exists because of a lack of LVL information in dispersions aggravating uncertainties in the inversion procedure, and (3) energy of high modes in a low-frequency range is very weak, so that it is difficult to estimate the cutoff frequency accurately, and “mode-crossing” occurs between the second higher and third higher modes when a HVL exists. 相似文献
3.
本文介绍了根据反射地震数据进行波速成象的一种方法,其基础为多种反演技术的综合。由于要求的波速图象C（x,z）具有间断性,除利用走时数据T（x,t）外,在地层比较水平的情况下,还利用了均方根速度V（x,t）和统计子波W（t）的数据来成象。计算机层析成象过程分为三步:首先重做速度分析,取得与初次反射走时一致的均方根速度数据;然后用反射走时与均方根速度联合反演对应分析道的层速度和界面深度;最后由联合反演结果和反射面走时求波速图象函数的数字化版。文中还给出了波速成象方法在我国西北某沉积盆地上的应用及验证结果。 相似文献
4.
本文系统地求解了不连续走时曲线的反演问题.所得到的解答表明,除了在低速层(若存在)内,其余各处的速度分布以及低速层的厚度均可唯一地确定.本文指出Slichter(1932)关于低速层厚度上界的推导和Gerver-Markushevich(1966)的反演公式是不合理的. 相似文献
5.
在Ⅰ测线观测资料的基础上,本文结合测区特点给出了高程校正及沉积层校正方法,利用反射波与莫霍面折射波运动学特性进行一维、二维地壳结构的正、反演计算,获得了自甘肃夏河至灵台460公里范围内的分层地壳结构,认为西部地壳厚度约51公里,东部仅为44公里,已属鄂尔多斯地台区结构,与地槽区的分界线大致对应地表陇县固关镇。武山与张家川之间在15-20公里深处存在低速层。结合地震波的动力学特性,可判定本测线通过的超壳断层有鄂尔多斯西缘-六盘山断裂、西秦岭北缘断裂和武山南北向断裂。 相似文献
6.
初至波走时层析反演技术作为建立近地表速度模型的重要手段,是解决陆地资料复杂静校正问题的关键技术。而折射波广泛发育的海洋地震资料,对折射波信息的关注与运用并没有得到广泛的重视。本文首次将层析反演方法应用于海洋拖缆地震数据的近海底速度模型的建立。本文方法与陆地资料层析反演的主要区别在于:①在震源信号的最小相位化处理后进行初至时间的拾取,避免了混合相位子波初至拾取不准带来的误差;②以海水深度与海水速度作为反演约束条件,减小了迭代误差。实测二维资料的层析反演结果表明,本文方法可反演出较为精确的海洋地层速度结构。 相似文献
7.
Field static corrections in general need be applied to all onshore seismic reflection data to eliminate the disturbing effects a weathering layer or near-surface low velocity zone has on the continuity of deep seismic reflections. The traveltimes of waves refracted at the bottom of the low velocity zone (or intermediate refracting interfaces) can often be observed as first breaks on shot records and used to develop a laterally inhomogeneous velocity model for this layer, from which the field static corrections can then be obtained. A simple method is described for computing accurate field statics from first breaks. It is based on a linearization principal for traveltimes and leads to the algorithms that are widely and successfully applied within the framework of seismic tomography. We refine an initial model for the low velocity layer (estimated by a standard traveltime inversion technique) by minimizing the errors between the observed first arrivals on field records and those computed by ray theory through an initial model of the low velocity layer. Thus, one can include more lateral velocity variations within the low velocity layers, which are important to obtain good field static corrections. Traditional first break traveltime inversion methods cannot, in general, provide such refined velocity values. The technique is successfully applied to seismic data from the Amazon Basin. It is based on a simple model for the low velocity layer that consists of an undulating earth surface and one planar horizontal refractor overlain by a laterally changing velocity field. 相似文献
8.
地震波的运动学信息(走时、斜率等)通常用于宏观速度建模.针对走时反演方法,一个基本问题是走时拾取或反射时差的估计.对于成像域反演方法,可以通过成像道集的剩余深度差近似计算反射波时差.在数据域中,反射地震观测数据是有限频带信号,如果不能准确地确定子波的起跳时间,难以精确地确定反射波的到达时间.另一方面,如果缺乏关于模型的先验信息,则很难精确测量自地下同一个反射界面的观测数据同相轴和模拟数据同相轴之间的时差.针对走时定义及时差测量问题,首先从叠前地震数据的稀疏表达出发,利用特征波场分解方法,提取反射子波并估计局部平面波的入射和出射射线参数.进一步,为了实现自动和稳定的走时拾取,用震相的包络极值对应的时间定义反射波的到达时,实现了立体数据中间的自动生成.理论上讲,利用包络极值定义的走时大于真实的反射波走时,除非观测信号具有无限带宽(即delta脉冲).然而,走时反演的目的是估计中-大尺度的背景速度结构,因此走时误差导致的速度误差仍然在可以接受的误差范围内.利用局部化传播算子及特征波聚焦成像条件将特征波数据直接投影到地下虚拟反射点,提出了一种新的反射时差估计方法.既避免了周期跳跃现象以及串层等可能性,又消除了振幅因素对时差测量的影响.最后,在上述工作基础之上,提出了一种基于特征波场分解的新型全自动反射走时反演方法(CWRTI).通过对泛函梯度的线性化近似,并用全变差正则化方法提取梯度的低波数部分,实现了背景速度迭代反演.在理论上,无需长偏移距观测数据或低频信息、对初始模型依赖性低且计算效率高,可以为后续的全波形反演提供可靠的初始速度模型.理论和实际资料的测试结果证明了本文方法的有效性. 相似文献
9.
利用青海和甘肃地震台网2007—2009年记录的远震波形资料,提取多频段P波接收函数,反演得到了青藏高原东北缘及相邻地块下方0~100km深度的地壳和上地幔S波速度结构.结果表明:(1)青藏高原东北缘的上、下地壳之间普遍存在一个S波速度低速层,其深度由南端的约35km向北变浅约为20km,推测该低速层为一壳内滑脱层,表明东北缘地区的上地壳变形与下地壳解耦,从滑脱层的深度分布可以认为青藏高原东北缘的地壳缩短自南向北进行,现阶段以上地壳增厚为主;(2)昆仑—西秦岭造山带的下地壳厚度较北侧的祁连地块薄,一种推测是西秦岭造山带的下地壳抗变形能力更强,也可能这种差异在块体拼合前已经存在;(3)青藏高原东北缘及鄂尔多斯和阿拉善地块的下地壳S波速度随深度的增加而增加,这种正梯度增加的S波速度结构反映较高黏滞性的下地壳,推测青藏高原东北缘的地壳结构不利于下地壳流的发育. 相似文献
10.
Derecke Palmer 《Geophysical Prospecting》2010,58(2):203-218
Non‐uniqueness occurs with the 1D parametrization of refraction traveltime graphs in the vertical dimension and with the 2D lateral resolution of individual layers in the horizontal dimension. The most common source of non‐uniqueness is the inversion algorithm used to generate the starting model. This study applies 1D, 1.5D and 2D inversion algorithms to traveltime data for a syncline (2D) model, in order to generate starting models for wave path eikonal traveltime tomography. The 1D tau‐p algorithm produced a tomogram with an anticline rather than a syncline and an artefact with a high seismic velocity. The 2D generalized reciprocal method generated tomograms that accurately reproduced the syncline, together with narrow regions at the thalweg with seismic velocities that are less than and greater than the true seismic velocities as well as the true values. It is concluded that 2D inversion algorithms, which explicitly identify forward and reverse traveltime data, are required to generate useful starting models in the near‐surface where irregular refractors are common. The most likely tomogram can be selected as either the simplest model or with a priori information, such as head wave amplitudes. The determination of vertical velocity functions within individual layers is also subject to non‐uniqueness. Depths computed with vertical velocity gradients, which are the default with many tomography programs, are generally 50% greater than those computed with constant velocities for the same traveltime data. The average vertical velocity provides a more accurate measure of depth estimates, where it can be derived. Non‐uniqueness is a fundamental reality with the inversion of all near‐surface seismic refraction data. Unless specific measures are taken to explicitly address non‐uniqueness, then the production of a single refraction tomogram, which fits the traveltime data to sufficient accuracy, does not necessarily demonstrate that the result is either ‘correct’ or the most probable. 相似文献
11.
Cross double‐difference inversion for simultaneous velocity model update and microseismic event location 
下载免费PDF全文
下载免费PDF全文 Microseismic monitoring is an approach for mapping hydraulic fracturing. Detecting the accurate locations of microseismic events relies on an accurate velocity model. The one‐dimensional layered velocity model is generally obtained by model calibration from inverting perforation data. However, perforation shots may only illuminate the layers between the perforation shots and the recording receivers with limited raypath coverage in a downhole monitoring problem. Some of the microseismic events may occur outside of the depth range of these layers. To derive an accurate velocity model covering all of the microseismic events and locating events at the same time, we apply the cross double‐difference method for the simultaneous inversion of a velocity model and event locations using both perforation shots and microseismic data. The cross double‐difference method could provide accurate locations in both the relative and absolute sense, utilizing cross traveltime differences between P and S phases over different events. At the downhole monitoring scale, the number of cross traveltime differences is sufficiently large to constrain events locations and velocity model as well. In this study, we assume that the layer thickness is known, and velocities of P‐ and S‐wave are inverted. Different simultaneous inversion methods based on the Geiger's, double‐difference, and cross double‐difference algorithms have been compared with the same input data. Synthetic and field data experiments suggest that combining both perforation shots and microseismic data for the simultaneous cross double‐difference inversion of the velocity model and event locations is available for overcoming the trade‐offs in solutions and producing reliable results. 相似文献
12.
地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力. 相似文献
13.
14.
V.B. Piip 《Geophysical Prospecting》2001,49(4):461-482
A method using simple inversion of refraction traveltimes for the determination of 2D velocity and interface structure is presented. The method is applicable to data obtained from engineering seismics and from deep seismic investigations. The advantage of simple inversion, as opposed to ray‐tracing methods, is that it enables direct calculation of a 2D velocity distribution, including information about interfaces, thus eliminating the calculation of seismic rays at every step of the iteration process. The inversion method is based on a local approximation of the real velocity cross‐section by homogeneous functions of two coordinates. Homogeneous functions are very useful for the approximation of real geological media. Homogeneous velocity functions can include straight‐line seismic boundaries. The contour lines of homogeneous functions are arbitrary curves that are similar to one another. The traveltime curves recorded at the surface of media with homogeneous velocity functions are also similar to one another. This is true for both refraction and reflection traveltime curves. For two reverse traveltime curves, non‐linear transformations exist which continuously convert the direct traveltime curve to the reverse one and vice versa. This fact has enabled us to develop an automatic procedure for the identification of waves refracted at different seismic boundaries using reverse traveltime curves. Homogeneous functions of two coordinates can describe media where the velocity depends significantly on two coordinates. However, the rays and the traveltime fields corresponding to these velocity functions can be transformed to those for media where the velocity depends on one coordinate. The 2D inverse kinematic problem, i.e. the computation of an approximate homogeneous velocity function using the data from two reverse traveltime curves of the refracted first arrival, is thus resolved. Since the solution algorithm is stable, in the case of complex shooting geometry, the common‐velocity cross‐section can be constructed by applying a local approximation. This method enables the reconstruction of practically any arbitrary velocity function of two coordinates. The computer program, known as godograf , which is based on this theory, is a universal program for the interpretation of any system of refraction traveltime curves for any refraction method for both shallow and deep seismic studies of crust and mantle. Examples using synthetic data demonstrate the accuracy of the algorithm and its sensitivity to realistic noise levels. Inversions of the refraction traveltimes from the Salair ore deposit, the Moscow region and the Kamchatka volcano seismic profiles illustrate the methodology, practical considerations and capability of seismic imaging with the inversion method. 相似文献
15.
B. URSIN 《Geophysical Prospecting》1986,34(8):1213-1218
The one-dimensional seismic inverse problem consists of recovering the acoustic impedance (or reflectivity function) as a function of traveltime from the reflection response of a horizontally layered medium excited by a plane-wave impulsive source. Most seismic sources behave like point sources, and the data must be corrected for geometrical spreading before the inversion procedure is applied. This correction is usually not exact because the geometrical spreading is different for primary and multiple reflections. An improved algorithm is proposed which takes the geometrical spreading from a point source into account. The zero-offset reflection response from a stack of homogeneous layers of variable thickness is used to compute the thickness, velocity and density of each layer. This is possible because the geometrical spreading contains additional information about the velocities. 相似文献
16.
利用中国数字测震台网（CDSN）记录到的台湾地区两个地震事件6°~30°震中距范围的三重震相波形资料,基于观测与理论波形拟合法,获到华南地区上地幔P波和S波的最佳波形拟合速度模型及其VP/VS比值.与AK135模型相比,华南地区410 km深度上方存在明显低速层：S波低速区厚度约为70 km,速度降为2%~5%;而P波低速区厚度为70~230 km,速度降为5%~6%.另外,410 km间断面整体表现为一个梯度层,厚度约为10~40 km,VP跃增量为4.0%~5.4%,而VS跃增量为2.6%~11.7%.研究区内,低速层的VP和VS异常值大小和410 km间断面速度跃变量由北向南逐步减小.结合以往的接收函数和地震层析成像结果,华南地区410 km间断面上方的低速区可能与太平洋俯冲板块脱水有关. 相似文献
17.
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. 相似文献
18.
Alexey Stovas 《Geophysical Prospecting》2010,58(2):219-227
The factorized velocity model that incorporates both vertical heterogeneity and constant anisotropy is one of the complicated analytical models used in seismic data processing and interpretation. In this paper, I derive the analytic equations for offset, traveltime and relative geometrical spreading for the quasi‐compressional (qP‐) waves that can be used for modelling and inversion of the traveltime parameters. I show that the presence of anelliptic anisotropy usually dominates over the vertical heterogeneity with respect to the non‐hyperbolicity of the factorized velocity model. 相似文献
19.
3-D simultaneous inversion for velocity and reflector geometry using multiple classes of arrivals in a spherical coordinate frame 总被引：1，自引：0，他引：1
Guo-jiao Huang Chao-ying Bai Xing-wang Li Stewart Greenhalgh 《Journal of Seismology》2014,18(1):123-135
The two key requirements in conducting 3-D simultaneous traveltime tomography on real data at the regional and global scale with multiple classes of arrival time information are (1) it needs an efficient and accurate arrival tracking algorithm for multiply transmitted, reflected (or refracted) and converted waves in a 3-D variable velocity model with embedded velocity discontinuities (or subsurface interfaces), and (2) a subdimensional inversion solver is required which can easily search for different types of model parameters to balance the trade-off between the different types of model parameter updated in the simultaneous inversion process. For these purposes, we first extend a popular grid/cell-based wavefront expanding ray tracing algorithm (the multistage irregular shortest-path ray tracing method), which previously worked only in Cartesian coordinate at the local scale, to spherical coordinates appropriate to the regional or global scale. We then incorporated a fashionable inversion solver (the subspace method) to formulate a simultaneous inversion algorithm, in which the multiple classes of arrivals (including direct and reflected arrivals from different velocity discontinuities) can be used to simultaneously update both the velocity fields and the reflector geometries. Numerical tests indicate that the new inversion method is both applicable and flexible in terms of computational efficiency and solution accuracy, and is not sensitive to a modest level of noise in the traveltime data. It offers several potential benefits over existing schemes for real data seismic imaging. 相似文献
20.
本文对波动方程首先进行富里叶—贝赛尔积分变换,在波数k域内构成(z,t)的有限差分隐格式进行迭代,由此计算出纵向非均匀的层状模型的合成地震图。对含有低速层和高速薄层的几种模型做了对比计算,通过时间场与空间场的波动分析,揭示了几种主要震相的传播与形成过程。计算结果表明,无论高速层的厚薄如何。反射波始终很强烈。但初至首波在薄层构造中不清晰,一种属转换型的续至薄层首波震相值得注意;低速层的顶界面难以形成能量较强的上行波,因此在推断低速层埋深上存在不确定性。 相似文献