共查询到20条相似文献，搜索用时 31 毫秒
1.
The authors present a method for estimation of interval velocities using the downward continuation of the wavefield to perform
layer-stripping migration velocity analysis. The generalized, phase-shift migration MG(F-K) in wavenumber-frequency domain
was used for fulltime downward extrapolation of the wavefield. Such downward depth extrapolation accounts for strong changes
of velocity in lateral and vertical directions and helps in correct positioning of the wavefield image in complex structures.
Determination of velocity is the recursive process which means that the wavefield on depth level z
n−1 (n = 0, 1, ...) is an input data-set for determination of velocity on level z
n
. The velocity ν [x, z
n
− z
n−1] can be thus treated as interval velocity in Δz
n
= z
n
− z
n−1 step. This method was tested on synthetic Marmousi data-set and showed satisfactory results for complex, inhomogeneous media. 相似文献
2.
垂直地震剖面(Vertical Seismic Profiling,VSP)资料处理中波场分离是关键问题之一.随着属性提取技术的发展,新的属性参数(例如Q值)提取技术对波场分离的保真性要求越来越高.本文改进了传统奇异值分解(Singular Value Decomposition,SVD)法,给出了一种对波场的动力学特征具有更好的保真性,可以作为Q值提取的预处理步骤的零偏VSP资料上下行波场分离方法.该方法通过两步奇异值分解变换实现:第一步,排齐下行波同相轴,利用SVD变换压制部分下行波能量;第二步,在剩余波场中排齐上行波同相轴,使用SVD变换提取上行波场.在该方法的实现过程中,压制部分下行波能量后的剩余波场中仍然存在较强的下行波干扰,使得上行波同相轴的排齐比较困难.本文给出了一种通过极大化多道数据线性相关程度(Maximize Coherence,MC)排齐同相轴的算法,在一定程度上解决了低信噪比下排齐同相轴的问题.将本文提出的方法用于合成数据和实际资料的处理,并与传统SVD法的处理结果进行对比,结果表明本文提出的波场分离方法具有良好的保真性,得到波场的质量明显优于传统SVD法.通过对本文方法和传统SVD法处理合成数据得到的下行波场提取Q值,然后进行对比可知,本文方法可以有效提高所提取Q值的准确性,适合作为Q值提取的预处理步骤. 相似文献
3.
In the case of 3D multilayered structures the 2D interval velocity analysis may be inaccurate. This fact is illustrated by synthetic examples. The method proposed solves the 3D inverse problem within the scope of the ray approach. The solution, i.e. the interval velocities and the reflection interface position, is obtained using data from conventional 2D line profiles arbitrarily located and from normal incidence time maps. Although the input information is essentially limited, the method presented reveals only minor biased velocity estimates. In order to implement the proposed 3D inversion method, we developed a processing procedure. The procedure performs the evaluation of reflection time and ray parameters along line profiles, 3D interval velocity estimation, and time-to-depth map migration. Tools to stabilize the 3D inversion are investigated. The application of the 3D inversion technique to synthetic and real data is compared with results of the 2D inversion. 相似文献
4.
Multicomponent seismic data are acquired by orthogonal geophones that record a vectorial wavefield. Since the single component recordings are not independent, the processing should be performed jointly for all the components. A way to achieve this goal is to exploit quaternions, hyper‐complex numbers that due to their very nature are apt to represent multidimensional data. In fact, quaternion algebra allows us to extend coherence functionals used for scalar observations to multicomponent data. Therefore by means of quaternions we implement semblance and other methods based on matched filtering and on the data covariance properties. As an application we show the results from a quaternion velocity analysis carried out combining information from the geophones and from the hydrophones of an ocean bottom cable (OBC) survey, and thus recognizing the true vectorial nature of the incoming wavefield. This also allows one to relax, at least partially, vector fidelity constraints. We demonstrate that quaternion velocity analysis yields an improved resolution with respect to the single component velocity analysis for any coherence functional chosen and that it simultaneously evidences velocity trends pertaining to different wave modes. This facilitates the interpreter in the estimation of interval Vp/Vs by means of event correlation, and in making use of a priori information from VSP and well logs. It also speeds up the velocity picking that can be performed in a single pass on a multicomponent velocity panel, rather than once for each single component velocity panel. 相似文献
5.
E. LOINGER 《Geophysical Prospecting》1983,31(1):98-118
Variations of seismic interval velocities within the cable length cause anomalies in the stacking velocity analyses. Utilizing the approximation of rectilinear ray propagation, i.e. supposing that the velocity changes cause time delays only, it is shown that the stacking velocity anomalies are linearly related to the interval velocity variations. In particular, the stacking velocity anomaly is calculated when the interval velocity of an intermediate layer undergoes a stepwise variation. The amplitude of the anomaly increases with the ratio between horizon depth and cable length. From the forward model, a program for the inversion is derived in order to identify lateral changes of interval velocities from unsmoothed stacking velocity analyses. Some examples of the application of this technique to synthetic and real data are presented. 相似文献
6.
针对在矿山微地震监测过程中纵横波及其传播速度难以确定的难点,在球坐标系中以波动方程为基础,推导出了矿山微地震全波自动定位的层析成像的计算公式,并给出了参数选择准则,藉此可以获得微地震成像能量最大值及其对应参数。通过判断能量最大值在风险区内的位置,可以将风险区内外的微地震事件区分出来。对于能量最大值位于风险区域内的微地震事件,其震源位置和震动时间以及传播速度与能量最大值对应的坐标参数相同,能量最大值对应于震源函数相对强度;当能量最大值位于风险区域边界上时,微地震事件不在风险区域范围之内。理论模型试算和矿山实际微地震资料的应用表明:层析成像方法无需纵横波识别和传播速度假定,只要给定事故风险区域范围及其传播速度范围,便可完成微地震事件定位,是一种全波全自动定位方法,具有进一步研究价值。 相似文献
7.
Elastic full waveform inversion of seismic reflection data represents a data‐driven form of analysis leading to quantification of sub‐surface parameters in depth. In previous studies attention has been given to P‐wave data recorded in the marine environment, using either acoustic or elastic inversion schemes. In this paper we exploit both P‐waves and mode‐converted S‐waves in the marine environment in the inversion for both P‐ and S‐wave velocities by using wide‐angle, multi‐component, ocean‐bottom cable seismic data. An elastic waveform inversion scheme operating in the time domain was used, allowing accurate modelling of the full wavefield, including the elastic amplitude variation with offset response of reflected arrivals and mode‐converted events. A series of one‐ and two‐dimensional synthetic examples are presented, demonstrating the ability to invert for and thereby to quantify both P‐ and S‐wave velocities for different velocity models. In particular, for more realistic low velocity models, including a typically soft seabed, an effective strategy for inversion is proposed to exploit both P‐ and mode‐converted PS‐waves. Whilst P‐wave events are exploited for inversion for P‐wave velocity, examples show the contribution of both P‐ and PS‐waves to the successful recovery of S‐wave velocity. 相似文献
8.
The estimation of velocity and depth is an important stage in seismic data processing and interpretation. We present a method for velocity-depth model estimation from unstacked data. This method is formulated as an iterative algorithm producing a model which maximizes some measure of coherency computed along traveltimes generated by tracing rays through the model. In the model the interfaces are represented as cubic splines and it is assumed that the velocity in each layer is constant. The inversion includes the determination of the velocities in all the layers and the location of the spline knots. The process input consists of unstacked seismic data and an initial velocity-depth model. This model is often based on nearby well information and an interpretation of the stacked section. Inversion is performed iteratively layer after layer; during each iteration synthetic travel-time curves are calculated for the interface under consideration. A functional characterizing the main correlation properties of the wavefield is then formed along the synthetic arrival times. It is assumed that the functional reaches a maximum value when the synthetic arrival time curves match the arrival times of the events on the field gathers. The maximum value of the functional is obtained by an effective algorithm of non-linear programming. The present inversion algorithm has the advantages that event picking on the unstacked data is not required and is not based on curve fitting of hyperbolic approximations of the arrival times. The method has been successfully applied to both synthetic and field data. 相似文献
9.
10.
In this paper, regionalP-wave upper mantle structure is investigated using slant-stack velocity analysis of short-period earthquake data recorded at station MAJO (Matsushiro, Japan). Shallow earthquakes from 1980–1986 within 35° of MAJO are used to construct a common receiver gather. Processing of the wavefield data includes focal depth and static time corrections, as well as deterministic deconvolution, in order to equalize pulse shapes and align wavelets on the first arrivals. The processed wavefield data are slant stacked and interatively downward continued to obtain a regional upper mantle velocity model. The model includes a low velocity zone between 107 and 220 km. Beneath the LVZ, the velocity increases smoothly down to the discontinuity at 401 km. In the transition zone, the velocity model again increases linearly, although there is some suggestion of further complexity in the downward continued wavefield data. At the base of the transition zone, a second velocity discontinuity occurs at 660 km, with a linear velocity gradient below. In addition to slant-stack analysis, travel times and synthetic seismograms are computed and compared with the processed and unprocessed wavefield data. 相似文献
11.
Subsalt imaging is strongly dependent on the quality of the velocity model. However, rugose salt bodies complicate wavefield propagation and lead to subsalt multipathing, illumination gaps and shadow zones, which cannot be handled correctly by conventional traveltime‐based migration velocity analysis (MVA). We overcome these limitations by the wave‐equation MVA technique, introduced in a companion paper, and demonstrate the methodology on a realistic synthetic data set simulating a salt‐dome environment and a Gulf of Mexico data set. We model subsalt propagation using wave paths created by one‐way wavefield extrapolation. Those wave paths are much more accurate and robust than broadband rays, since they inherit the frequency dependence and multipathing of the underlying wavefield. We formulate an objective function for optimization in the image space by relating an image perturbation to a perturbation of the velocity model. The image perturbations are defined using linearized prestack residual migration, thus ensuring stability, relative to the first‐order Born approximation assumptions. Synthetic and real data examples demonstrate that wave‐equation MVA is an effective tool for subsalt velocity analysis, even when shadows and illumination gaps are present. 相似文献
12.
13.
Up–down wavefield decomposition is effectuated by a scaled addition or subtraction of the pressure and vertical particle velocity, generally on horizontal or vertical surfaces, and works well for data given on such surfaces. The method, however, is not applicable to decomposing a wavefield when it is given at one instance in time, i.e. on snapshots. Such situations occur when a wavefield is modelled with methods like finite-difference techniques, for the purpose of, for example, reverse time migration, where the entire wavefield is determined per time instance. We present an alternative decomposition method that is exact when working on snapshots of an acoustic wavefield in a homogeneous medium, but can easily be approximated to heterogeneous media, and allows the wavefield to be decomposed in arbitrary directions. Such a directional snapshot wavefield decomposition is achieved by recasting the acoustic system in terms of the time derivative of the pressure and the vertical particle velocity, as opposed to the vertical derivative in up–down decomposition for data given on a horizontal surface. As in up–down decomposition of data given at a horizontal surface, the system can be eigenvalue decomposed and the inverse of the eigenvector matrix decomposes the wavefield snapshot into fields of opposite directions, including up–down decomposition. As the vertical particle velocity can be rotated at will, this allows for decomposition of the wavefield into any spatial direction; even spatially varying directions are possible. We show the power and effectiveness of the method by synthetic examples and models of increasing complexity. 相似文献
14.
Xu Zhaotao Wang Huazhong 《应用地球物理》2006,3(4):226-233
Based on perturbation theory, the wave equation extrapolation operator with mixed domains has the ability to deal with lateral velocity variations. It is the image method that has undergone much research in seismology. All extrapolation operators face the problem of choosing the reference velocity due to continuation in depth. The wavefield extrapolation operator with a single reference velocity is suitable for media with weak lateral variation. The multi-reference velocity extrapolation operator can cope with severe lateral velocity variations and improve image accuracy. However, the calculation cost is large. We present a self-adaptive approach to automatically determine the number of selected reference velocities according to the complexity of structure and the given velocity threshold value. The approach can be used to construct the SSF, FFD, WXFD, and GSP multi-reference velocity wavefield extrapolation image algorithms. The result of a salt-dome model data test demonstrates that the self-adoptive multi-reference wavefield extrapolation algorithm has the ability to deal with severe lateral velocity variations and can also be used for structure edge detection. The method is flexible and computationally cost-effective. 相似文献
15.
基于扰动理论的混合域波动方程波场外推算子,具有一定介质横向速度变化适应能力,是反射地震学中研究较多的成像方法。此类波场外推算子沿深度层进行波场外推,都存在参考速度选择问题。单参考速度波场外推算子,适应地下介质弱横向变速,而多参考速度波场外推算子可以提高横向变速的适应能力和成像精度,但要以大量计算为代价。本文提出的自适应多参考速度选择策略,根据外推层地质构造的复杂度和给定的速度门槛值自动选择参考速度个数,利用该策略构造混合域SSF、FFD、WXFD和GSP等多参考速度波场外推成像算法。盐丘模型理论数据测试结果表明,自适应多参考速度波场外推算法具有强横向变速适应能力和较高成像精度。 相似文献
16.
S. BOISSE 《Geophysical Prospecting》1978,26(1):163-174
It is not possible to determine accurate geological velocities from seismic velocity analysis for thin layers or complex structural features, especially under an unconformity. Instead, we can approach the problem of interval velocity with seismic amplitudes analysis and compute the reflection coefficient along the unconformity surface. An error estimation has been made on a model to test the possibility of such a method and to choose the best parameters to be used. The method has been applied on an actual case: the computed interval velocities show good correlation with the values obtained by a sonic log. 相似文献
17.
Tariq Alkhalifah 《Geophysical Prospecting》2011,59(4):627-634
The wavefield dependence on a virtual shift in the source location can provide information helpful in velocity estimation and interpolation. However, the second‐order partial differential equation (PDE) that relates changes in the wavefield form (or shape) to lateral perturbations in the source location depends explicitly on lateral derivatives of the velocity field. For velocity models that include lateral velocity discontinuities this is problematic as such derivatives in their classical definition do not exist. As a result, I derive perturbation partial differential wave equations that are independent of direct velocity derivatives and thus, provide possibilities for wavefield shape extrapolation in complex media. These PDEs have the same structure as the wave equation with a source function that depends on the background (original source) wavefield. The solutions of the perturbation equations provide the coefficients of a Taylor's series type expansion for the wavefield. The new formulas introduce changes to the background wavefield only in the presence of lateral velocity variation or in general terms velocity variations in the perturbation direction. The accuracy of the representation, as demonstrated on the Marmousi model, is generally good. 相似文献
18.
Jürgen Mann Rainer Jger Thilo Müller German Hcht Peter Hubral 《Journal of Applied Geophysics》1999,42(3-4)
The simulation of a zero-offset (ZO) stack section from multi-coverage reflection data is a standard imaging method in seismic processing. It significantly reduces the amount of data and increases the signal-to-noise ratio due to constructive interference of correlated events. Conventional imaging methods, e.g., normal moveout (NMO)/dip moveout (DMO)/stack or pre-stack migration, require a sufficiently accurate macro-velocity model to yield appropriate results, whereas the recently introduced common-reflection-surface stack does not depend on a macro-velocity model. For two-dimensional seismic acquisition, its stacking operator depends on three wavefield attributes and approximates the kinematic multi-coverage reflection response of curved interfaces in laterally inhomogeneous media. The common-reflection-surface stack moveout formula defines a stacking surface for each particular sample in the ZO section to be simulated. The stacking surfaces that fit best to actual events in the multi-coverage data set are determined by means of coherency analysis. In this way, we obtain a coherency section and a section of each of the three wavefield attributes defining the stacking operator. These wavefield attributes characterize the curved interfaces and, thus, can be used for a subsequent inversion. In this paper, we focus on an application to a real land data set acquired over a salt dome. We propose three separate one-parametric search and coherency analyses to determine initial common-reflection-surface stack parameters. Optionally, a subsequent optimization algorithm can be performed to refine these initial parameters. The simulated ZO section obtained by the common-reflection-surface stack is compared to the result of a conventional NMO/DMO/stack processing sequence. We observe an increased signal-to-noise ratio and an improved continuity along the events for our proposed method — without loss of lateral resolution. 相似文献
19.
The seismic K-Horizon is the key to gaining understanding on the deep supercritical geothermal rocks in Southern Tuscany. The K-Horizon is hosted in metamorphic rocks, which cause strong seismic wavefield scattering resulting in a poor signal-to-noise ratio. Our study aims to reveal high-resolution seismic images of the K-Horizon below a geothermal field in Southern Tuscany, using an advanced three-dimensional seismic depth imaging approach. The key seismic pre-processing steps in the time domain include muting a large amount of persistent noise based on the statistical analysis of the seismic amplitudes, and tomostatics technique to correct for static effects. We carried out seismic depth imaging using Kirchhoff Pre-Stack Depth Migration and Fresnel Volume Migration techniques. Each migration technique was tested with constant and heterogeneous three-dimensional velocity models. Due to the difficulties in determining emergent angles for this low signal-to-noise ratio data set, the migration results with the heterogeneous three-dimensional velocity model show less coherent reflections compared to the migration results using the constant velocity model. Both velocity models however lead to relatively the same structure and depth of the K-Horizon, indicating the similarity of the average velocities along the wave propagation paths in both velocity models. With both velocity models Fresnel Volume Migration yields the K-Horizon with better reflection coherency and higher signal-to-noise ratio than standard Kirchhoff Pre-Stack Depth Migration. Nevertheless, both migration techniques have been able to reveal the K-Horizon with relatively high resolution and provide a reliable basis for geothermal rock characterization as well as steering of the first geothermal well penetrating the K-Horizon. 相似文献
20.
Ali A. Nowroozi 《Pure and Applied Geophysics》1990,133(1):103-115
In reflection surveys and velocity analysis, calculations of interval velocities and layer-thicknesses of a multilayered horizontal structure are often based on Dix's equation which requires the travel times at zero offsets and a prior estimate of the root mean squared velocities.In this paper a method is presented which requires only the reflection travel-time data. A set of equations are derived which relate the interval velocity and thickness of a layer to the reflection travel time from the top and the bottom of that layer, the offset distances and the ray parameter. It is shown that the difference of the offset distances and the difference of the picked travel times of any reflected rays with the same value of ray parameter from the top and the bottom of a horizontal layer can be used to calculate the interval velocity and thickness of that layer. 相似文献