共查询到20条相似文献，搜索用时 46 毫秒
1.
The improvement in accuracy and efficiency of wave-equation migration techniques is an ongoing topic of research. The main problem is the correct imaging of steeply dipping reflectors in media with strong lateral velocity variations. We propose an improved migration method which is based on cascading phase-shift and finite-difference operators for downward continuation. Due to these cascaded operators we call this method‘Fourier finite-difference migration’(FFD migration). In our approach we try to generalize and improve the split-step Fourier migration method for strong lateral velocity variations using an additional finite-difference correction term. Like most of the current migration methods in use today, our method is based on the one-way wave equation. It is solved by first applying the square-root operator but using a constant velocity at each depth step which has to be the minimum velocity. In a second step, the approximate difference between the correct square-root operator and this constant-velocity squareroot operator (the error made in the first step) is implemented as an implicit FD migration scheme, part of which is the split-step Fourier correction term. Some practical aspects of the new FFD method are discussed. Its performance is compared with that of split-step and standard FD migration schemes. First applications to synthetic and real data sets are presented. They show that the superiority of FFD migration becomes evident by migrating steeply dipping reflectors with complex overburden having strong lateral velocity variations. If velocity is laterally constant, FFD migration has the accuracy of the phase-shift method. The maximum migration angle is velocity adaptive, in contrast to conventional FD migration schemes. It varies laterally depending on the local level of velocity variation. FFD migration is more efficient than higher-order implicit FD schemes. These schemes use two cascaded downward-continuation steps in order to attain comparable migration performance. 相似文献
2.
Conventional two‐way splitting Fourier finite‐difference migration for 3D complex media yields azimuthal anisotropy where an additional phase correction is needed with much increase of computational cost. We incorporate the alternating‐direction‐implicit plus interpolation scheme into the conventional Fourier finite‐difference method to reduce azimuthal anisotropy. This scheme retains the high‐order remnants ignored by the two‐way splitting in the form of a wavefield interpolation in the wavenumber domain. The wavefield interpolation for each step of downward extrapolation is implemented between the wavefields before and after the conventional Fourier finite‐difference extrapolation. As the Fourier finite‐difference migration is implemented in the space and wavenumber dual space, the Fourier transforms between space and wavenumber domain that were needed for the alternating‐direction‐implicit plus interpolation in frequency domain (FD) migration are saved in Fourier finite‐difference migration. Since the azimuth anisotropy in Fourier finite‐difference is much less than that in FD, the application of the alternating‐direction‐implicit plus interpolation scheme in Fourier finite‐difference migration is superior to that in FD migration in handling complex media with large velocity contrasts and steep dips. Impulse responses show that the presented method reduces the azimuthal anisotropy at almost no extra cost. 相似文献
3.
Andrey V. Masjukov 《Geophysical Prospecting》2003,51(4):333-345
A 3D F–K dip-moveout (DMO) is developed, which is applicable to data acquired in an elementary single-fold cross-spread. The key idea is that a 3D log-stretch transform and the inherent regularity of the cross-spread geometry make it possible to transform 3D Fourier DMO. The derived theory generalizes the 2D Fourier shot-gather DMO in the log-stretch domain; 2D turns out to be a special case. Similarly to 2D, the cross-spread DMO becomes convolutional after multidimensional logarithmic stretch. The proposed method works for orthogonal and slanted acquisition geometries; the cross-spread DMO relationships are found to be independent of the intersection angle of the shot and receiver lines. In contrast to integral (Kirchhoff-style) methods, the cross-spread F–K DMO does not degrade from the inevitable irregularity in 3D sampling of offsets in a CMP gather. The newly derived F–K DMO operator can be approximated by finite-difference (FD) schemes; the low-order FD cross-spread DMO equation is shown to be the 3D extension of the Bolondi and Rocca offset continuation. It is shown that F–K and low-order FD operators are effective in a synthetic case. 相似文献
4.
本文提出了复杂构造地区的目标导向观测系统的设计方法.使用波动方程正演模拟来指导并在二维声波方程的一阶速度-应力方程中应用交错网格有限差分法实现.使用了四阶精度的差分算子和完全匹配层吸收边界条件.通过分析理论模型的模拟结果,展示了如何将地面地震响应与地下目标构造匹配.通过分析桥口地区实际地质模型的模拟结果,指出波动方程正演模拟在小断块、小背斜增生的复杂地区中相对于传统方法更精确,图像更清晰,更利于分析和指导观测系统设计.跟踪目标区的反射同相轴并得到其到达地面的接收范围来约束炮检距的范围,通过比较不同道距的模拟记录并综合考虑成本和任务目标来选取最佳道距.最终获得了实际效果达到设计要求的观测系统参数. 相似文献
5.
Reverse-time migration (RTM) is based on seismic numerical modeling algorithms, and the accuracy and efficiency of RTM strongly depend on the algorithm used for numerical solution of wave equations. Finite-difference (FD) methods have been widely used to solve the wave equation in seismic numerical modeling and RTM. In this paper, we derive a series of time–space domain staggered-grid FD coefficients for acoustic vertical transversely isotropic (VTI) equations, and adopt these difference coefficients to solve the equations, then analyze the numerical dispersion and stability, and compare the time–space domain staggered-grid FD method with the conventional method. The numerical analysis results demonstrate that the time–space domain staggered-grid FD method has greater accuracy and better stability than the conventional method under the same discretizations. Moreover, we implement the pre-stack acoustic VTI RTM by the conventional and time–space domain high-order staggered-grid FD methods, respectively. The migration results reveal that the time–space domain staggered-grid FD method can provide clearer and more accurate image with little influence on computational efficiency, and the new FD method can adopt a larger time step to reduce the computation time and preserve the imaging accuracy as well in RTM. Meanwhile, when considering the anisotropy in RTM for the VTI model, the imaging quality of the acoustic VTI RTM is better than that of the acoustic isotropic RTM. 相似文献
6.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献
7.
常规伪谱方法二阶时间差分格式时间精度较低,且对于大步长时间采样间隔,常规伪谱方法不稳定.拟解析方法对于速度变化剧烈的模型,在时间和空间上均有较大误差.本文提出了一种基于解耦的二阶位移弹性波方程波场模拟及矢量波场分解的优化拟解析方法,将归一化的拟拉普拉斯算子分别应用于P波和S波波场延拓,延拓矢量波场的同时,可分解并延拓纯纵波和纯横波波场.利用弹性波优化拟微分算子表示拟拉普拉斯算子,该拟微分算子不仅包括原始微分算子的谱估计,而且还包含一个时间补偿项,其可在波数-空间域精确地补偿波动方程在时间方向上采用二阶有限差分引起的误差.利用低秩分解近似求解弹性波优化拟微分算子,可有效提高计算效率.2D均匀模型、层状模型以及部分盐丘模型数值正演模拟结果表明:相比较于常规的伪谱法和拟解析法,本文方法在时间和空间上均有很高的精度,并且稳定性条件比较宽松. 相似文献
8.
基于共聚焦点技术的叠前AVP（振幅随射线参数变化）分析与常规叠后反演方法相比优势明显,但传统通过褶积和互相关运算来实现的方法依赖于聚焦算子,而在复杂构造区走时计算困难且子波难以精确提取,从而导致了聚焦算子不准确,而且褶积和互相关运算会影响信噪比和分辨率,基于此,本文提出了基于保真振幅单程波延拓算法的叠前AVP成像方法.该方法利用保真振幅傅里叶有限差分延拓算法实现两步聚焦,分别生成共聚焦点道集和网格点道集,既充分利用了保真振幅延拓算法在振幅保持方面的优势,也可以发挥傅里叶有限差分方法对复杂构造区横向变速适应性强的优势,而且两步聚焦过程都不需要聚焦算子,从而解决了传统方法中走时计算和子波提取的问题.模型试算结果表明了方法的正确性和可行性,而针对实际地震资料的试处理结果与传统方法相比具有更高的信噪比和分辨率,表明了方法的有效性.该方法为复杂构造区油气检测提供了一种新的地球物理依据. 相似文献
9.
与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性. 相似文献
10.
Ｔｈｅｃｏｎｖｏｌｕｔｉｏｎａｌｄｉｆｅｒｅｎｔｉａｔｏｒｍｅｔｈｏｄｆｏｒｎｕ┐ｍｅｒｉｃａｌｍｏｄｅｌｉｎｇｏｆａｃｏｕｓｔｉｃａｎｄｅｌａｓｔｉｃｗａｖｅ┐ｆｉｅｌｄｓＺＨＯＮＧ－ＪＩＥＺＨＡＮＧ（张中杰），ＪＩ－ＷＥＮＴＥＮＧ（滕吉文）ａｎ... 相似文献
11.
传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性. 相似文献
12.
Acoustic finite-difference modeling beyond conventi-onal Courant-Friedrichs-Lewy stability limit: App-roach based on variable-length temporal and spatial operators 
下载免费PDF全文
下载免费PDF全文 Conventional finite-difference (FD) methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy (CFL) numbers 0.707 and 0.577 for two-dimensional (2D) and three-dimensional (3D) equal spacing cases, respectively, thereby limiting time step selection. Based on the definition of temporal and spatial FD operators, we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy. First, to simulate wave propagation beyond the conventional CFL stability limit, the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones. Second, to preserve the modeling accuracy, the velocity-dependent lengths of the temporal and spatial operators are adaptively varied. The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples, respec-tively. We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media. 相似文献
13.
O. HOLBERG 《Geophysical Prospecting》1988,36(2):99-114
Numerical wavefield extrapolation represents the backbone of any algorithm for depth migration pre- or post-stack. For such depth imaging techniques to yield reliable and interpretable results, the underlying wavefield extrapolation algorithm must propagate the waves through inhomogeneous media with a minimum of numerically induced distortion, over a range of frequencies and angles of propagation. A review of finite-difference (FD) approximations to the acoustic one-way wave equation in the space-frequency domain is presented. A straightforward generalization of the conventional FD formulation leads to an algorithm where the wavefield is continued downwards with space-variant symmetric convolutional operators. The operators can be precomputed and made accessible in tables such that the ratio between the temporal frequency and the local velocity is used to determine the correct operator at each grid point during the downward continuation. Convolutional operators are designed to fit the desired dispersion relation over a range of frequencies and angles of propagation such that the resulting numerical distortion is minimized. The optimization is constrained to ensure that evanescent energy and waves propagating at angles higher than the maximum design angle are attenuated in each extrapolation step. The resulting operators may be viewed as optimally truncated and bandlimited spatial versions of the familiar phase shift operator. They are unconditionally stable and can be applied explicitly. This results in a simple wave propagation algorithm, eminently suited for implementation on pipelined computers and on large parallel computing systems. 相似文献
14.
15.
傅里叶有限差分法三维波动方程正演模拟 总被引：4，自引：6，他引：4
傅里叶有限差分(FFD)法兼有相位屏法和隐式有限差分法二者的优势，能够处理复杂地质构造中的波传播问题，但在三维情形下，算子的双向分裂会引起明显的方位各向异性误差.本文用Fourier变换计算双向分裂过程中的高阶交叉项，消除了方位各向异性误差.该方法充分利用了FFD法在双域实现的算法结构，明显减少了由于引入误差校正所带来的计算量.将该方法应用于修改后的三维French模型的地震正演问题，并将得到的叠后记录、单炮记录同全波有限差分法的模拟结果进行对比，结果证实了该方法对一次反射波具有较高的模拟精度，在内存需求和计算效率方面则具有更大的优势. 相似文献
16.
叠前逆时偏移在理论上是现行偏移方法中最为精确的一种成像方法,其实现过程中的核心步骤之一是波动方程的波场延拓,而波场延拓的本质是求解波动方程,所以精确、快速地求解波动方程对逆时偏移至关重要.本文采用一种基于时空域频散关系的有限差分方法来求解声波方程,分析其频散和稳定性,实现波场数值模拟,并将分析和模拟结果与传统有限差分法进行对比.分析结果和模型数值模拟结果都表明时空域有限差分法模拟精度更高、稳定性更好.将时空域高阶有限差分法应用到叠前逆时偏移波场延拓的方程求解中,然后再利用归一化互相关成像条件成像,理论模型数据偏移处理获得了精度更高的成像.同时,在逆时偏移波场延拓的实现中,采用自适应变长度的空间差分算子求解空间导数的有限差分策略,在不影响数值模拟和成像精度的前提下,有效地提高了计算效率. 相似文献
17.
Various numerical methods have been used in the literature to simulate single and multiphase flow in fractured media. A promising approach is the use of the discrete-fracture model where the fracture entities in the permeable media are described explicitly in the computational grid. In this work, we present a critical review of the main conventional methods for multiphase flow in fractured media including the finite difference (FD), finite volume (FV), and finite element (FE) methods, that are coupled with the discrete-fracture model. All the conventional methods have inherent limitations in accuracy and applications. The FD method, for example, is restricted to horizontal and vertical fractures. The accuracy of the vertex-centered FV method depends on the size of the matrix gridcells next to the fractures; for an acceptable accuracy the matrix gridcells next to the fractures should be small. The FE method cannot describe properly the saturation discontinuity at the matrix–fracture interface. In this work, we introduce a new approach that is free from the limitations of the conventional methods. Our proposed approach is applicable in 2D and 3D unstructured griddings with low mesh orientation effect; it captures the saturation discontinuity from the contrast in capillary pressure between the rock matrix and fractures. The matrix–fracture and fracture–fracture fluxes are calculated based on powerful features of the mixed finite element (MFE) method which provides, in addition to the gridcell pressures, the pressures at the gridcell interfaces and can readily model the pressure discontinuities at impermeable faults in a simple way. To reduce the numerical dispersion, we use the discontinuous Galerkin (DG) method to approximate the saturation equation. We take advantage of a hybrid time scheme to alleviate the restrictions on the size of the time step in the fracture network. Several numerical examples in 2D and 3D demonstrate the robustness of the proposed model. Results show the significance of capillary pressure and orders of magnitude increase in computational speed compared to previous works. 相似文献
18.
Broadband constant-coefficient propagators 总被引：4，自引：1，他引：4
Li-Yun Fu 《Geophysical Prospecting》2005,53(3):299-310
The phase error between the real phase shift and the Gazdag background phase shift, due to lateral velocity variations about a reference velocity, can be decomposed into axial and paraxial phase errors. The axial phase error depends only on velocity perturbations and hence can be completely removed by the split‐step Fourier method. The paraxial phase error is a cross function of velocity perturbations and propagation angles. The cross function can be approximated with various differential operators by allowing the coefficients to vary with velocity perturbations and propagation angles. These variable‐coefficient operators require finite‐difference numerical implementation. Broadband constant‐coefficient operators may provide an efficient alternative that approximates the cross function within the split‐step framework and allows implementation using Fourier transforms alone. The resulting migration accuracy depends on the localization of the constant‐coefficient operators. A simple broadband constant‐coefficient operator has been designed and is tested with the SEG/EAEG salt model. Compared with the split‐step Fourier method that applies to either weak‐contrast media or at small propagation angles, this operator improves wavefield extrapolation for large to strong lateral heterogeneities, except within the weak‐contrast region. Incorporating the split‐step Fourier operator into a hybrid implementation can eliminate the poor performance of the broadband constant‐coefficient operator in the weak‐contrast region. This study may indicate a direction of improving the split‐step Fourier method, with little loss of efficiency, while allowing it to remain faster than more precise methods such as the Fourier finite‐difference method. 相似文献
19.
本文在前人工作的基础上,建立了一种基于Shannon奇异核的交错网格褶积微分算子方法.文中不仅详细讨论了影响算子精度的各种因素,同时也着重分析了其在弹性波模拟中的频散关系和稳定性条件.通过和交错网格有限差分算子比较,发现该算子即使在高波数域也具有较高的精度.均匀介质中的数值试验也表明,该方法9点格式就基本上达到了解析解精度.而分层均匀介质和复杂介质中的地震波数值模拟也同时证实了该方法精度高,稳定性好,是一种研究复杂介质中地震波传播的有效数值方法. 相似文献
20.
We present a two‐dimensional (2D) gradient operator that produces more accurate results than known traditional operators such as Ando, Sobel and the so‐called Isotropic operator. We further extend the derivation to three‐dimensional (3D), a powerful feature missing in all conventional operators. We start by constructing a parameterized formula that generically represents all 2D numerical gradient operators. We then solve for the required parameter by equating this numerical gradient with that obtained analytically from a single Fourier harmonic (or, equivalently here, a stationary plane wave). As this parameter is frequency‐ and direction‐dependent (by virtue of the underlying Fourier harmonic), we construe a pragmatic version of it that is independent of these two variables yet capable of significantly reducing the error associated with traditional operators. Extension to 3D is achieved similarly; it requires dealing with two parameters as opposed to only one in the 2D case. Synthetic and real‐data results confirm higher accuracy from this operator than from traditional ones. 相似文献