共查询到20条相似文献，搜索用时 562 毫秒
1.
2.
深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构. 相似文献
3.
井间地震资料具有极高的分辨率,但是,其波场十分复杂,特别是有效一次反射波场能量较弱,且往往被管波等强相干干扰淹没,波场的识别与分离难度较大.地震物理模拟技术是认识复杂地震波场的有效手段.我们通过单一地质体模型、复杂地质体模型和真实井间模拟模型等多种逼真地质模型的物理模拟,揭示了井间地震观测下直达波、透射波、反射波、折射波,多次波和导波等各类波型的特征.物理模拟试验表明,在炮检域平面上分析直达波的信噪比与能量分布,可以使处理人员对旅行时的误差有更直观的认识.所以识别初至与拾取旅行时间时,不仅要在共震源点道集(CSG)上分析拾取,而且还要在共炮检距道集(COG)上分析拾取,这样就提高了直达波识别的可靠性与拾取时间的精度. 相似文献
4.
The interpretation of deep seismic sounding (DSS) data has been made on the basis of a two-dimensional inhomogeneous model. The refracted first arrivals as well as reflected and diffracted waves on the seismic records have been utilized. The seismic section was modeled in the iso-veolcity lines v(x, y) = const, taking into account the zones of diffraction associated with deep faults. Gravity observations have been used to construct a block model of the Earth's crust with vertical boundaries. It is suggested to define the base of the crust as the zone with velocities between 7.8 and 8.2 km/s. The reflecting boundaries of different length occurring in this zone can be conformal or unconformal with the iso-velocity lines near the base of the crust. As an example of our approach to the interpretation of DSS data the folded-blocky structure of the crust with horizontal inhomogeneities of velocity and density is shown in the Kzyl-Orda-Dzheskazgan profile in Central Kazakhstan. 相似文献
5.
根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震（OBS）共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。 相似文献
6.
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. 相似文献
7.
2009年,中国地质科学院地质研究所与美国俄克拉荷马大学合作实施了一条长453 km的深地震反射、宽角反射与折射、三分量反射地震联合探测剖面. 剖面南起怀来盆地,向北依次穿过燕山造山带西缘、内蒙地轴、白乃庙弧带、温都尔庙杂岩带,到达索伦缝合带. 其中,宽角反射与折射剖面采用8个0.5~1.5 t炸药震源激发,使用300套Texan单分量数字检波器接收,获得了高质量的地震资料. 通过资料分析和处理,识别出沉积层及结晶基底的折射波（Pg）、来自上地壳底界面的反射波（Pcp）,中地壳底界面的反射波（Plp）,莫霍界面的反射波（Pmp）及上地幔顶部的折射波（Pn）等5个震相. 分别采用Hole有限差分层析成像和Rayinvr算法对华北克拉通北缘及中亚造山带南部进行了上地壳P波速度结构成像和全地壳二维射线追踪反演成像. 结果显示：（1）中亚造山带地壳厚度~40 km,变化平缓,低于全球平均造山带地壳平均厚度,可能为造山后区域伸展的结果. 阴山-燕山带附近莫霍明显加深,推测其为燕山期造山过程形成的山根,但该山根很可能在后期被改造. （2）测线中部地壳上部速度较高,对应地表大面积花岗岩出露,而下地壳速度较低,速度梯度低,呈通道状,推测其可能曾为古亚洲洋向南俯冲消亡的主动陆缘,并在碰撞后演变为伸展环境下岩浆侵入的通道. （3）华北克拉通北缘与中亚造山带显示出不同速度变化特征,前者变化相对缓而后者则变化剧烈,二者的分界出现在赤峰-白云鄂博断裂附近. 相似文献
8.
The construction of S-wave velocity models of marine sediments down to hundreds of meters below the seafloor is important
in a number of disciplines. One of the most significant trends in marine geophysics is to use interface waves to estimate
shallow shear velocities which play an important role in determining the shallow crustal structure. In marine settings, the
waves trapped near the fluid–solid interface are called Scholte waves, and this is the subject of the study. In 1998, there
were experiments on the Ninetyeast Ridge (Central Indian Ocean) to study the shallow seismic structure at the drilled site.
The data were acquired by both ocean bottom seismometer and ocean bottom hydrophone. A new type of seafloor implosion sources
has been used in this experiment, which successfully excited fast and high frequency (>500 Hz) body waves and slow, intermediate
frequency (<20 Hz) Scholte waves. The fundamental and first higher mode Scholte waves have both been excited by the implosion
source. Here, the Scholte waves are investigated with a full waveform modeling and a group velocity inversion approach. Shear
wave velocities for the uppermost layers of the region are inferred and results from the different methods are compared. We
find that the full waveform modeling is important to understand the intrinsic attenuation of the Scholte waves between 1 and
20 Hz. The modeling shows that the S-wave velocity varies from 195 to 350 m/s in the first 16 m of the uppermost layer. Depths
levels of high S-wave impedance contrasts compare well to the layer depth derived from a P-wave analysis as well as from drilling
data. As expected, the P- to S-wave velocity ratio is very high in the uppermost 16 m of the seafloor and the Poisson ratio
is nearly 0.5. Depth levels of high S-wave impedance contrasts are comparable to the layer depth derived from drilling data. 相似文献
9.
The interpretation of seismic refraction and wide angle reflection data usually involves the creation of a velocity model based on an inverse or forward modelling of the travel times of crustal and mantle phases using the ray theory approach. The modelling codes differ in terms of model parameterization, data used for modelling, regularization of the result, etc. It is helpful to know the capabilities, advantages and limitations of the code used compared to others.This work compares some popular 2D seismic modelling codes using the dataset collected along the seismic wide-angle profile DOBRE-4, where quite peculiar/uncommon reflected phases were observed in the wavefield.The ~505 km long profile was realized in southern Ukraine in 2009, using 13 shot points and 230 recording stations. Double PMP phases with a different reduced time (7.5–11 s) and a different apparent velocity, intersecting each other, are observed in the seismic wavefield. This is the most striking feature of the data. They are interpreted as reflections from strongly dipping Moho segments with an opposite dip. Two steps were used for the modelling. In the previous work by Starostenko et al. (2013), the trial-and-error forward model based on refracted and reflected phases (SEIS83 code) was published. The interesting feature is the high-amplitude (8–17 km) variability of the Moho depth in the form of downward and upward bends. This model is compared with results from other seismic inversion methods: the first arrivals tomography package FAST based on first arrivals; the JIVE3D code, which can also use later refracted arrivals and reflections; and the forward and inversion code RAYINVR using both refracted and reflected phases. Modelling with all the codes tested showed substantial variability of the Moho depth along the DOBRE-4 profile. However, SEIS83 and RAYINVR packages seem to give the most coincident results. 相似文献
10.
11.
A new algorithm for tomographic inversion of traveltimes of reflected and refracted seismic waves is developed. The inversion gives interface configurations and velocity distributions in layers. The important features of the algorithm are: (a) the inclusion of shot time delays in the list of unknown parameters; (b) the regularization is applied in such a way that the most probable model is characterized by the similarity of neighbouring interfaces. As the problem under consideration is non-linear, several iterations are necessary in order to obtain the final model. In the case of a very inexact initial model, a 'layer-by-layer' inversion strategy is recommended as a first inversion step. The inversion program is supplied with a user interface, thanks to which one can: (a) pick interactively and identify seismic traveltimes; (b) build and edit depth/velocity models; and (c) display calculated traveltime curves and compare them with picked traveltimes as well as with the original seismic sections. The efficiency of the inversion software developed is illustrated by a numerical example and a field example in which shallow seismic data are considered. Application to wide-aperture reflection/refraction profiling (WARRP) data is also possible. 相似文献
12.
13.
A least squares estimation procedure is used to estimate the pulse-shape, amplitude function, and arrival time of multiple reflected signals. The estimates of the multiple reflections are subtracted from the data which are subsequently processed by standard methods. The estimation algorithm is applied continuously along the seismic line for each shot point or common datum point. In some cases it is advantageous to apply a pulse-shaping filter prior to using the estimation algorithm. The effectiveness of the technique is demonstrated by studying common shot point gathers, velocity analyses, and stacked sections derived from field data. 相似文献
14.
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. 相似文献
15.
The main problem in seismic prospecting is to infer from the observed reflection response the distribution of density and seismic velocity with depth. This process is generally called the inversion of the reflection data. For plane waves propagating through plane parallel stratification, it can be shown that at any depth the ratio between the amplitude of the transmitted and reflected wave satisfies the Riccati equation. Based on this equation we have formulated an iterative inversion method, which is found to be suitable for numerical computations. We have applied this method on synthetic reflection data, and found that it provides a very fast and accurate inversion. 相似文献
16.
G. NEUMANN 《Geophysical Prospecting》1981,29(2):161-177
Lateral inhomogeneities generate fluctuations in the traveltime of seismic waves. By evaluation of these traveltime fluctuations from different source and receiver positions, lateral inhomogeneities can be located using a pseudo inverse matrix method (Aki, Christoffersson and Husebye 1977). The formulation of the problem is possible for transmitted waves as well as for reflected and refracted waves. In reflection seismics this method is of importance, if no reflections from the inhomogeneities themselves, but only reflections from lower boundaries can be observed. The basic assumptions for the mathematical formulation are (1) the average velocities and depths of the reflecting horizons are known already from standard processing methods, and (2) the traveltime residuals are due to lateral velocity changes between different reflectors or between reflectors and the surface. The area of the earth to be considered is divided into layers and the layers into rectangular blocks. The parallel displacement of a ray after passing a disturbed block is neglected, only the traveltime residual is taken into account. In this paper the method and its application to data obtained with two-dimensional models are described. 相似文献
17.
B. J. RAFISON 《Geophysical Prospecting》1988,36(8):847-856
Determination of petrography and pore fluid content is an ultimate goal of an integrated seismic-petrophysical study. For lack of a general inversion technique, forward modelling is useful in studying the relations between lithology, stratigraphy, pore fluid content and the seismic response. This report describes a study of two clastic sequences in Utah, from which 32 rock samples were analysed. A detailed petrographic study was done. Laboratory measurements were made of ultrasonic compressional- and shear-wave velocity as a function of pressure. We computed the velocities at seismic frequencies for the samples when dry, over-pressured, brine saturated, and oil saturated. The velocities were sensitive to the porosity, carbonate cementation and the depositional facies. We generated velocity profiles for hypothetical reservoirs for a range of saturation states. The velocity profiles were used to generate synthetic seismic shot gathers to study the seismic response of these clastic reservoirs. The fluid-saturation strongly affects the seismic respone, as does the presence of a coal seam. An amplitude change with offset is often observed. However, stratigraphy appears to have a stronger effect on the seismic response. 相似文献
18.
多震相走时和波形三维地震成象方法及在祁连山地壳结构研究中的应用 总被引：1，自引：0，他引：1
本文的理论方法是以几何射线理论为基础发展起来的、天然地震走时反演技术及天然地震层析成像技术。它存在着震源函数与介质参数的解耦问题。本研究采取了五种方法来改善反演结果。包括,利用,Pg,Pn等震相增大约束条件;用已有精度较高的人工地震测深结果作速度约柬：用波形反演来修改模型,把诸多物理量开发出来互为约束,以修改后的模型再作反演,使解的稳定性大大提高：采用最优化过程,选择遗传算法。可以进行震源定位,走时反演,波形反演：得到任意深度的速度分布及从地表到Moho面的速度剖面。用于在祁连山地区的结果表明,这些层析剖面对认识大地构造、重大深部事件动力学是很有益的。 相似文献
19.
以三维高分辨地震与海底高频地震仪（OBS）联合勘探数据为基础,获得海底之下沉积层的地震反射成像剖面及多波信息,并以此确定研究区含天然气水合物沉积层的纵、横波速度的变化特征.根据走时反演获得的横波速度与纵波速度对比分析发现,研究区海底之下500 m深度范围内的某些沉积层具有较高的纵横波速度,这一纵波速度升高区域与水合物稳定带对应,而纵波速度下降并且横波速度变化较小的区域,可能与游离气的存在相关.游离气的可能存在与基于这一区域2007年钻探测井结果的普遍认识不完全相符. 相似文献
20.
Three-component seismic and geoelectrical in-mine surveys were carried out in Lyukobanya colliery near Miskolc, Hungary to determine the in situ petrophysical parameter distributions and to detect inhomogeneities in the coal seam. The seismic measurements comprise an underground vertical seismic profile, using body waves, and an in-seam seismic amplitude-depth distribution and transmission survey, using channel waves. The geoelectrical measurements are based on the drift- and seam-sounding method. Interval traveltime-, amplitude-, multiple-filter- and polarization analysis methods are applied to the seismic data. They lead to a five-layer model for the strata including the coal seam. The coal seam and two underlying beds act as a seismic waveguide. The layer sequence supports the propagation of both normal and leaky mode channel waves of the Love- and Rayleigh type. A calculation of the total reflected energy for each interface using Knott's energy coefficients shows that the velocity ranges of high reflection energy and of normal and leaky mode wavegroups coincide. The excitation of wavegroups strongly depends on the seismic source. A simultaneous inversion of a geoelectrical drift- and seam-sounding survey prevents misinterpretations of the seismic data by clearly identifying the low-velocity coal seam as a high-resistivity bed. Calculations of dispersion and sounding curves improve the resolution of the slowness and resistivity in each layer. Both diminished amplitudes and distortions in the polarization of transmission seismo-grams and decreasing resistivities in a geoelectrical pseudosection of the coal seam are related to an inhomogeneity. A calculation of synthetic seismograms for Love and Rayleigh channel waves with the finite-difference and the Alekseev-Mikhailenko method agrees well with the field data for the main features, i.e., particular arrivals in the wave train, wavegroups, velocities and symmetries or asymmetries. This in-mine experiment demonstrates that the simultaneous acquisition, processing and interpretation of seismic and geoelectrical data improve the lithological interpretation of petrophysical parameter distributions. Coal seam inhomogeneities can also be detected more reliably by the two independent surveys than by one alone. 相似文献