共查询到18条相似文献，搜索用时 272 毫秒
1.
强横向变速和陡倾角介质的成像问题是地震偏移成像的难点.本文在频率-空间域利用有限差分高阶分裂法求解任意广角波动方程实现了地震偏移成像.试验表明,参考速度的选取是影响成像精度的关键因素.脉冲响应测试表明通过优化参数任意广角波动方程能够以较低阶的方程获得较大的偏移角度.Marmousi模型数据叠前深度偏移试验表明,该方法能... 相似文献
2.
为了处理横向强变速介质中的深度成像问题,本文提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,并在基于反射系数估算的成像条件下,可实现叠前深度偏移成像.该算子具有方程阶数低且能对陡倾角成像的特征,并采用有限差分法波场延拓,能适应速度场的任意变化.当在频率-空间域进行计算时,相对于纯粹的时间-空间域有限差分算法有计算效率高、成像方便的优点.脉冲响应测试和对Marmousi模型进行的叠前深度偏移结果表明,该偏移方法在强横向变速情况下具有非常好的成像效果. 相似文献
3.
《应用地球物理》2017,(1)
频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播,但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性,本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验,验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理,基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像,说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力,我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理,对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力,具有一定实际应用价值。 相似文献
4.
波动方程宽角抛物逼近得到的通常是非常系数的单程波传播算子，其系数是速度横向变化的函数，因此需要利用有限差分(FD)进行数值实施. 通过对Lippmann Schwinger单程波动积分方程的退化核逼近,本文研究了一类宽角退化算子的偏移成像. 这种退化偏移算子只用快速Fourier变换进行波场延拓，将常规的Fourier分裂步地震偏移方法(SSF)推广适应强速度横向变化介质和大角度传播波场. 退化的Fourier偏移算子通过在两个分裂步项之间作波数域线性插值来实现波场延拓，每延拓一层需要比常规的SSF地震偏移方法多一次快速Fourier变换(FFT). 通过SEG/EAGE盐丘模型和实际地震资料的应用表明，退化Fourier偏移算子能很好地对盐下的陡倾角断层和实际地震剖面上的复杂小断块和大断裂地质构造成像. 相似文献
5.
《地球物理学进展》2015,(3)
断层或盐体悬伸侧翼等陡倾角构造可以产生棱柱波,当地震波的速度随着深度线性增加时则会产生回折波,这两类地震特殊波包含了一次反射波所不能反映的地下构造的信息.逆时偏移是一种成像精度高和能够处理强横向速度变化的双程波动方程偏移方法,不对方程进行近似,能够对高陡构造进行准确成像,其应用越来越广泛.本文将逆时偏移和棱柱波、回折波相结合,充分利用棱柱波和回折波的信息,探索了逆时偏移成像特殊波的能力,分析了逆时偏移成像棱柱波时对速度场精度的要求,用棱柱波来描述盐丘侧翼的垂直边界.同时分析了不同成像条件下逆时偏移成像回折波的效果,并用回折波成像盐体内部和盐体遮挡下构造. 相似文献
6.
本文基于波场的局部平面波分解，得到空间和方向都具有局部性的局部平面波小波束，并由此对Fourier有限差分传播算子进行局部平面波分解.这种用小波束分解的单向传播算子，可被用来对采集系统进行定向照明分析和采集倾角响应估算，以分析采集系统布局和上覆结构对面向目标的叠前偏移成像质量的影响.进而，利用局部平面波成像矩阵及其成像相册，根据目标结构的主倾角方向，抽取并叠加部分局部平面波像，可以进行目标结构定向化特征成像.最后，对具有速度强横向变化上覆结构的陡倾角目标断层的复杂二维SEG-EAGE盐丘模型数据，进行采集系统目标照明分析和目标结构定向化特征成像. 相似文献
7.
地震波传播算子的计算效率和精度是制约三维叠前深度偏移的关键因素. 广义屏传播算子(GSP, Generalized Screen Propagator)是一种在双域中实现的广角单程波传播算子. 这一方法略去了在非均匀体之间发生的交混回响,但它可以正确处理包括聚焦、衍射、折射和干涉在内的各种多次前向散射现象. 通过背景速度下的相移和扰动速度下的陡倾角校正,广义屏算子能够适应地层速度的强烈横向变化. 这种算子可以直接应用于炮集叠前偏移,通过将广义屏算子作用于双平方根方程,还可以获得一种高效率、高精度的炮检距域叠前深度偏移方法,用于二维共炮检距道集和三维共方位角道集的深度域成像. 本文首先简述了炮检距域广义屏传播算子的理论,进而讨论了共照射角成像(CAI, Common Angle Imaging)条件,由此给出各个不同照射角(炮检距射线参数)下的成像结果,进而得到共照射角像集. 由于照射角和炮检距的对应关系,共照射角像集又为偏移速度分析和AVO(振幅随炮检距变化)分析等提供了有力工具. 相似文献
8.
傅里叶有限差分法三维波动方程正演模拟 总被引：4，自引：6，他引：4
傅里叶有限差分(FFD)法兼有相位屏法和隐式有限差分法二者的优势，能够处理复杂地质构造中的波传播问题，但在三维情形下，算子的双向分裂会引起明显的方位各向异性误差.本文用Fourier变换计算双向分裂过程中的高阶交叉项，消除了方位各向异性误差.该方法充分利用了FFD法在双域实现的算法结构，明显减少了由于引入误差校正所带来的计算量.将该方法应用于修改后的三维French模型的地震正演问题，并将得到的叠后记录、单炮记录同全波有限差分法的模拟结果进行对比，结果证实了该方法对一次反射波具有较高的模拟精度，在内存需求和计算效率方面则具有更大的优势. 相似文献
9.
基于任意广角波动方程的频率-空间域深度偏移方法研究 总被引：1，自引：1，他引：0
孙歧峰 《CT理论与应用研究》2010,19(2):17-24
传统的单程波动方程偏移算法对大倾角成像困难,本文基于时空域任意广角单程声波方程,通过对时间变量进行傅立叶变换得到其频率-空间域形式,利用有限差分进行离散化,设计并实现了频率-空间域有限差分叠后偏移成像算法。模型试算表明,该方法能够通过参数优化使得低阶数算子适应较大倾角地层的偏移成像。 相似文献
10.
同时震源（Simultaneous-source,SS）地震采集技术能有效地提高地震数据采集效率,但直接对SS混合数据偏移成像会在最后的成像剖面中引入很强的串扰噪声.将SS数据偏移成像看作一个反演问题,利用最小二乘（Least-squares,LS）求解是压制SS直接成像中串扰噪声的一种有效尝试.构造增强滤波（Structure-enhancing filter,SE）约束的最小二乘逆时偏移（LSRTM）方法可以有效地压制SS数据成像中的串扰噪声,但SE实质为低通滤波,会将成像中的陡倾角等细节信息平滑涂抹,降低成像分辨率.本文在利用SE对LSRTM约束的基础上,提出了基于加权构造增强约束的LSRTM方法（WSE-LSRTM）并应用于SS数据的反演成像中.该方法不仅能够有效地压制串扰噪声（cross-talk）、保留结构信息,而且可以保护成像中的陡倾角结构不被过度平滑而破坏.在对简单模型和复杂Marmousi模型的数值测试中,该方法都取得了良好的效果. 相似文献
11.
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. 相似文献
12.
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. 相似文献
13.
Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference (symplectic FFD) method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps. 相似文献
14.
研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎. 相似文献
15.
多次波偏移中的假象主要来自于不同地震事件之间的互相关,由于这种互相关满足成像条件,很难直接在偏移过程中去除.但是对于准确的速度模型,真实的成像结果在角度域内应该是平直的.根据这个判断准则,可以在角度域内移除多次波偏移中的假象.本文以数据自相关偏移为例,提出了在单程波多次波偏移中移除假象的主要流程:首先在在单程波偏移过程中高效地提取角度域共成像点道集,然后对角度域共成像点道集应用高分辨率的抛物线型Radon变换,用合适的切除函数处理后,反变换回到角度域,最后叠加各个角度成分,得到偏移结果.Marmousi模型的合成数据测试表明,这种方法可以很好地压制多次波偏移过程中产生的假象,有效地提高成像结果的信噪比. 相似文献
16.
单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合. 相似文献
17.
地震波在各向异性介质中以一个准P波(qP)和两个准S波(qS1和qS2)的形式传播.研究三种波的相速度、群速度以及偏振方向等传播性质能够为各向异性介质中的正反演问题提供有效支撑.具有比横向各向同性(TI)介质更一般对称性的正交各向异性介质通常需要9个独立参数对其进行描述,这使得对传播特征的计算更为复杂.当两个准S波速度相近时具有耦合性,从而令慢度的计算产生奇异性.因此,奇异点(慢度面的鞍点和交叉点)附近的反射与透射(R/T)系数的求解不稳定,会导致波场振幅不准确.本文首次通过结合耦合S波射线理论和基于迭代的各向异性相速度与偏振矢量的高阶近似解,得到了适用于正交各向异性介质以qP波入射所产生的二阶R/T系数的计算方法.与基于一阶近似的结果相比,基于二阶近似的方法提高了qP波R/T系数的精度,能得到一阶耦合近似无法表达的准确的qP-qS转换波的R/T系数解,且方法适用于较强的各向异性介质. 相似文献
18.
We propose an effective and reliable time-domain substructure technique which takes soil-structure interaction effects into account and uses the flexibility coefficients of unbounded soil obtained in the frequency domain. Compared with previous methods, the frequency points to calculate flexibility coefficients, and computational loads in the calculation of time-domain interaction forces, are reduced. In the formulation, we have assumed the flexibility coefficient to be a periodic function, obtained within the bandlimited frequency range, which also includes the predominant frequencies of the structure and incident wave. Then we simulate the periodic flexibility coefficients using discrete impulse responses in the time-domain analyses. However, the real and imaginary parts of the bandlimited flexibility coefficients do not form a Hilbert transform pair; the discrete impulse responses should be modified to be causal for the time-domain analyses. We present various discrete impulse responses which have been obtained from only the real part, only the imaginary part and from both the real and imaginary parts of the frequency-domain flexibility coefficients by FFT with causality conditions. Through a numerical example the relationship between the modified discrete impulse responses and the validity of the time-domain substructure method is presented. 相似文献