共查询到20条相似文献，搜索用时 156 毫秒
1.
高斯束偏移是一种准确、灵活、高效的深度域成像方法,其不但具有接近于波动方程偏移的成像精度,还保留了Kirchhoff偏移灵活、高效的特点以及对复杂地表条件良好的适应性.本文提出了一种适用于复杂地表条件的且具有相对振幅保持特点的高斯束偏移方法.通过考虑地表高程、倾角以及实际的道间距等信息,推导了基于高斯束表示的波场反向延拓公式,并结合反褶积成像条件,得到了复杂地表条件下的共炮域保幅高斯束偏移公式.同原有方法相比,本文方法不但可以直接在起伏的地表面进行局部平面波的分解,具有更高的成像精度,而且可以得到反映地下随角度变化反射系数的成像结果.数值模型的试算验证了上述结论. 相似文献
2.
试图运用波场动力学途径,旨在增强观测系统分析评估,本文提出基于单程波方程的角度域照明分析方法.本文方法不仅可用于观测系统设计,对叠前成像资料的考量、叠前振幅补偿和AVA（Amplitude Versus Angle,振幅随角度变化）反演也必具重要的意义.文中通过目的层邻域单个绕射点的正、反传波场分析,给出地下目标的水平和倾斜界面随角度变化的成像照明强度.单程波方法可用于模拟地震波在复杂速度构造下的传播,并得到较准确的幅值和多次到达的波场,用以对绕射点正、反传分析时可在照明分析中简明地利用幅值信息和多次到达的波场.本文方法可显现复杂盖层下目标的照明情况.本文中也对照明能量随传播距离的几何扩散进行补偿,突显了角度域的照明均匀度和范围,从而使得所分析的结果可直接应用于AVA的分析及补偿.针对照明分析的特点,文中建议可用一种快速的单程波波场延拓策略——频率速度相关变步长波场深度延拓.数值实验结果表明了本文方法是可行且有效的. 相似文献
3.
油气勘探的重点正转向复杂地表条件和复杂地质条件的区域.双复杂条件下的叠前深度偏移是解决复杂地表条件和复杂地质构造成像的有效手段.基于“逐步累加”的“直接下延”法是解决复杂地表成像的有效手段,能够较好地消除地形起伏的影响.波动方程频率空间域有限差分(xwfd)叠前深度偏移对介质速度横向变化有较强的适应性,适宜于复杂构造的偏移成像,同其他常规波动方程深度偏移一样,常规的xwfd偏移方法,主要也是针对相位进行波场延拓,没有对振幅做任何处理.我们基于保幅单程波方程,推导出了基于xwfd的保幅波场延拓算子,针对xwfd求解时引入误差的影响,我们在xwfd保幅波场延拓过程中加入了误差补偿,实现了带误差补偿的xwfd保幅偏移.基于带误差补偿的xwfd保幅算子,应用适合起伏地表的直接下延法,对双复杂介质模型和实际资料进行了试算,改善了双复杂介质的成像效果.其中,误差补偿可以在若干个外推步长上进行,所以相对于保幅傅里叶有限差分(ffd)法偏移来说,该方法在改善成像质量的同时,也具有较高的运算效率. 相似文献
4.
《应用地球物理》2017,(4)
研究了基于Poynting矢量的角度域逆时偏移成像及成像幅度的校正方法。根据Poynting矢量进行波场角度分解,由此构建局部成像矩阵及局部照明矩阵。在局部成像矩阵中建立的角度域成像条件,有效地消除了低波数干扰,同时可在局部成像矩阵中进行角度域共成像点道集抽取、倾角估计等运算。利用局部照明矩阵进行了基于全波波动方程的时域照明分析,在局部照明矩阵中计算倾角域幅度校正因子。根据逆时偏移的像计算共倾角像,利用校正因子对各角度的像进行校正,进而实现对成像结果的校正,从而实现了一种高效的倾角域幅度校正方法。最后通过SEG/EAGE模型进行数值计算验证了文中所述的计算方法。 相似文献
5.
常规的转换波静校正的基本思想都是从地震波的运动学特征出发,基于地表一致性假设.在地表条件复杂和地表高程相差较大的地区,它不仅无法解决严重的静校正问题,反而会带来新的畸变.本文基于频率波数域波动方程偏移原理,采用波场延拓方法实现转换波静校正,其关键点在于时间空间域和频率波数域的相对应.文中通过坐标变换将起伏地表转化为新坐标系下的水平地表,把炮点和检波点映射到同一水平面上,然后在新坐标系下推导频率域波动方程延拓公式,接着对下行波P和上行转换波SV分别利用近地表速度向上延拓到基准面,恢复起伏地表到基准面之间的真实波场,最后转换到原始坐标系取出基准面数据完成转换波静校正.通过对模拟和实际数据处理,证明该方法是正确和有效的. 相似文献
6.
提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(dreamlet算子)沿深度延拓.本文利用多维局部余弦变换实现dreamlet算法,分解后的波场系数和传播算子不仅有很好的稀疏性,且均为实数,也即波的传播和成像过程完全在实数域实现.文中推导了局部余弦基dreamlet波场分解和传播算子理论公式并将其应用于叠前深度偏移.在dreamlet相空间波的传播过程为稀疏矩阵相乘,而且延拓后的地表数据波场的有效时间长度随深度的增加不断减小,从而可以减少需要传播的波场系数.二维SEG/EAGE盐丘和SIGSBEE模型算例验证了理论推导的正确性,成像结果显示该方法在横向速度变化剧烈情况下有很好的精度. 相似文献
7.
角度域成像道集是叠前深度偏移的重要输出结果,它是偏移速度分析、各向异性分析和AVA分析等研究工作的基础.目前存在的角度域成像道集的生成方法受计算效率或角度分辨率的影响,仍然满足不了实际生产的要求.角度域成像道集的生成方法可以大致分为直接法和间接法两大类.在直接方法中,波矢量方向计算和局部平面波分解是两个最重要的内容,它们共同决定角度域成像条件的实现效率和角度域成像道集的质量.为了完善现有的角度域成像道集生成方法,本文提出一种新的波矢量计算方法和局部平面波分解方法.本文先用波动方程任意时刻的柯西条件构造一个只含原波场负频率成分的柯西波场,然后根据柯西波场在时间波数域的振幅谱计算波场的波矢量方向.该方法仅在需要计算波矢量方向的时刻合成柯西波场,不需要增加额外的数据读写操作,是一种高效的波矢量计算方法.本文还以柯西波场为基础提出一种高效的局部平面波分解方法,保证角度域成像条件的实现效率.结合柯西波场和局部平面波分解方法,本文最后给出一种新的角度域成像方法.文中最后的数值实验证明该方法得到的角度域成像道集具有理想的角度分辨率,可以反映地下构造的角度照明情况. 相似文献
8.
弹性波逆时偏移不受倾角和偏移孔径的限制,能够实现任意复杂构造的高精度多波成像,是目前最精确的多分量资料偏移成像方法之一.逆时偏移算法的核心是波场延拓,传统波场延拓以水平基准面为边界条件,基于固定采样步长进行规则网格剖分,采用阶梯近似法处理起伏地表和复杂构造界面时会产生台阶散射,严重影响起伏地表复杂构造的成像精度.基于无网格节点模型,定量分析了弹性波模拟中径向基函数有限差分法的频散关系和稳定性条件.基于此,提出一种基于QR径向基函数的高精度有限差分方法,并提出一种优化的起伏地表自适应节点剖分方法,推导了精确的无网格自由边界条件和弹性波无网格混合吸收边界条件,形成了新的基于无网格的起伏地表弹性波数值模拟方法.此外,本文将此无网格径向基函数有限差分方法应用于精确的纵横波场矢量分解公式,实现了起伏地表弹性波逆时偏移成像.通过对高斯山丘模型,起伏凹陷模型和起伏地表Marmousi-2模型进行数值试算,验证了本文方法的有效性和可行性. 相似文献
9.
传统炮检距域共像集（CIG）在复杂介质中因波传播的多路径而存在反射体位置不确定的问题. 角度域CIG由于克服了这一缺陷而逐步成为速度分析、AVA以及振幅保真偏移成像等研究的主要手段. 以波动理论为基础的地震偏移成像方法的发展为获得高质量的角度域CIG提供了可靠的实现途径. 其中,基于波场局域化分解和传播的小波束域波场延拓和偏移成像方法,因其波场分解基本函数和传播算子在空间和方向上的双重局域特性,而成为角度相关分析研究的有效工具. 本文在采用Gabor Daubechies框架分解的小波束叠前角度域偏移成像基础上,利用不同的叠加方法由局部角度域像矩阵得到了反射角域CIG（CRAIG）和倾角域CIG（CDAIG）. 以SEG EAGE二维盐体模型为例,通过对CRAIG和CDAIG的对比,探讨了这两种角度域CIG的特点及其在地震偏移成像中的潜在应用. 相似文献
10.
随着能源和资源勘查开采工作的深入，地形强烈起伏的盆山耦合地区的地震资料处理解释技术正日益成为山地地震勘探面临的重要挑战.逆时偏移方法作为精确的地震偏移成像方法之一，能对地下结构进行高精度成像.逆时偏移的核心是地震波场延拓，由于传统的地震波场延拓技术往往基于水平地表条件，相应的方法在直接处理强地形起伏条件下的地震资料时往往存在一定的精度损失.本文引入一种精度无损的处理起伏边界的模型参数化方法：基于贴体网格的地形"平化"策略发展了与地形有关的地震波波动方程数值模拟方法，采用零延迟归一化互相关成像条件实现了起伏地表条件下的弹性波场逆时偏移成像.对工业界的标准Marmousi模型和盐丘模型进行改造，获得了相应起伏地形条件下的复杂几何模型，开展了起伏地表下的地震偏移成像数值试验.结果表明基于贴体网格"平化"策略的逆时偏移成像方法具有较高的灵活性，可适应不同类型起伏地表采集的地震资料，显示出该方法在地震勘探领域的良好应用前景. 相似文献
11.
12.
基于炮集数据的常规叠前深度偏移波场外推的计算量非常大,控制照明叠前深度偏移是基于快速的组合炮技术,通过在目的层上控制震源波场,高效实现对地质目标的高质量成像. 组合炮所需的合成算子,由在目的层定义的震源波场反向外推到地表建立. 本文提出一种更快速的地表旋转控制照明方法,直接在地表对震源波场进行面向目标的照明控制,实现目标区域高质量成像. 该方法只需对某一核合成算子进行简单的几何旋转建立合成算子函数,避免了复杂的波场外推. 这种新的控制照明方法与快速的一般组合炮偏移效率相同,但对复杂目标区域的成像能力显著增强. 应用本文方法对Marmousi模型进行计算,得到深层复杂构造非常理想的成像结果. 相似文献
13.
Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration 
下载免费PDF全文
下载免费PDF全文 The key objective of an imaging algorithm is to produce accurate and high‐resolution images of the subsurface geology. However, significant wavefield distortions occur due to wave propagation through complex structures and irregular acquisition geometries causing uneven wavefield illumination at the target. Therefore, conventional imaging conditions are unable to correctly compensate for variable illumination effects. We propose a generalised wave‐based imaging condition, which incorporates a weighting function based on energy illumination at each subsurface reflection and azimuth angles. Our proposed imaging kernel, named as the directional‐oriented wavefield imaging, compensates for illumination effects produced by possible surface obstructions during acquisition, sparse geometries employed in the field, and complex velocity models. An integral part of the directional‐oriented wavefield imaging condition is a methodology for applying down‐going/up‐going wavefield decomposition to both source and receiver extrapolated wavefields. This type of wavefield decomposition eliminates low‐frequency artefacts and scattering noise caused by the two‐way wave equation and can facilitate the robust estimation for energy fluxes of wavefields required for the seismic illumination analysis. Then, based on the estimation of the respective wavefield propagation vectors and associated directions, we evaluate the illumination energy for each subsurface location as a function of image depth point and subsurface azimuth and reflection angles. Thus, the final directional‐oriented wavefield imaging kernel is a cross‐correlation of the decomposed source and receiver wavefields weighted by the illuminated energy estimated at each depth location. The application of the directional‐oriented wavefield imaging condition can be employed during the generation of both depth‐stacked images and azimuth–reflection angle‐domain common image gathers. Numerical examples using synthetic and real data demonstrate that the new imaging condition can properly image complex wave paths and produce high‐fidelity depth sections. 相似文献
14.
提出了用局部指数标架小波束进行角度域分解的方法,解决了局部余弦基和局部正弦基缺乏单一方向性的问题.局部指数标架由局部余弦基和局部正弦基线性组合而成,是冗余度为2的紧标架.利用局部余弦变换和局部正弦变换的快速算法,能使基于局部指数标架进行方向照明分析的计算效率较常用的局部倾斜叠加和Gabor-Daubechies标架等方法具有更为明显的优势.通过计算二维SEG/EAGE模型和SIGSBEE模型的方向照明图以及采集系统倾角响应图证实了本文方法的有效性.该方法的高效性使其在三维模型的方向照明分析和大规模的工业应用中具有广阔前景. 相似文献
15.
采用弹性波全波形反演方法精确重建深部金属矿多参数模型,建模过程采用基于地震照明的反演策略.首先给出基于照明理论的观测系统可视性定义,利用可视性分析构建新的目标函数,对反演目标可视性较高的炮检对接收到的地震记录在波场匹配时占有更高的权重,确保了参与反演计算中的地震数据的有效性;其次将给定观测系统对地下介质的弹性波场照明强度作为优化因子,根据地震波在波阻抗界面处的能量分配特点,自适应补偿波场能量分布和优化速度梯度,以提高弹性波全波形反演过程的稳定性和反演结果的精度.理论模型和金属矿模型反演试验结果表明,基于可视性分析和能量补偿的反演策略可以使弹性波全波形反演更快地收敛到目标函数的全局极小值,获得适用于金属矿高分辨率地震偏移成像的多参数模型. 相似文献
16.
Unequal illumination of the subsurface highly impacts the quality of seismic imaging. Different image points receive different folds of reflection‐angle illumination, which can be caused by irregular acquisition or by wave propagation in complex media. Illumination problems can deteriorate amplitudes in migrated images. To address this problem, we present a method of stacking angle‐domain common‐image gathers, in which we use local similarity with soft thresholding to determine the folds of local illumination. Normalization by local similarity regularizes local illumination of reflection angles for each image point of the subsurface model. This approach compensates for irregular illumination by selective stacking in the image space, regardless of the cause of acquisition or propagation irregularities. Additional migration is not required because the methodology is implemented in the reflection angle domain after migration. We use two synthetic examples to demonstrate that our method can normalize migration amplitudes and effectively suppress migration artefacts. 相似文献
17.
We propose a method based on the Poynting vector that combines angle-domain imaging and image amplitude correction to overcome the shortcomings of reverse-time migration that cannot handle different angles during wave propagation. First, the local image matrix (LIM) and local illumination matrix are constructed, and the wavefield propagation directions are decomposed. The angle-domain imaging conditions are established in the local imaging matrix to remove low-wavenumber artifacts. Next, the angle-domain common image gathers are extracted and the dip angle is calculated, and the amplitude-corrected factors in the dip angle domain are calculated. The partial images are corrected by factors corresponding to the different angles and then are superimposed to perform the amplitude correction of the final image. Angle-domain imaging based on the Poynting vector improves the computation efficiency compared with local plane-wave decomposition. Finally, numerical simulations based on the SEG/EAGE velocity model are used to validate the proposed method. 相似文献
18.
受地下复杂构造和地震数据采集系统的影响，使地震波对地下目标的照明出现不均匀性，地震采集系统难以有效地获取地下某些目标的反射信息，进而使数据偏移成像在这些目标体上出现成像阴影. 根据波场和Green函数的窗口Fourier框架展开，利用角度域波动方程偏移成像和波动方程照明分析，并结合波动方程反演理论，提出一种角度域波动方程偏移成像阴影照明补偿方法. 这种补偿方法能同时考虑地震数据采集系统和波场传播路径对偏移成像的影响，消除复杂构造区的偏移成像阴影，改进波动方程叠前深度偏移成像在复杂构造区的成像效果. 相似文献
19.
Attenuation in seismic wave propagation is a common cause for poor illumination of subsurface structures. Attempts to compensate for amplitude loss in seismic images by amplifying the wavefield may boost high‐frequency components, such as noise, and create undesirable imaging artefacts. In this paper, rather than amplifying the wavefield directly, we develop a stable compensation operator using stable division. The operator relies on a constant‐Q wave equation with decoupled fractional Laplacians and compensates for the full attenuation phenomena by performing wave extrapolation twice. This leads to two new imaging conditions to compensate for attenuation in reverse‐time migration. A time‐dependent imaging condition is derived by applying Q‐compensation in the frequency domain, whereas a time‐independent imaging condition is formed in the image space by calculating image normalisation weights. We demonstrate the feasibility and robustness of the proposed methods using three synthetic examples. We found that the proposed methods are capable of properly compensating for attenuation without amplifying high‐frequency noise in the data. 相似文献
20.
现有的微震记录直接成像方法是将微震记录既当作入射记录，也当作散射记录，从而实现偏移成像.但此方法并不能突出透射波所携带的来自震源下方的深层散射波信息.本文在假设已知微震位置与子波的前提下，提出了对微震下方构造进行逆时偏移的成像方法.该方法类似于常规的逆时偏移，只是震源位置在地下.这使得在成像时，地下更深部的入射波场相比震源在地表时会更为精确，因此能够获得更加准确的成像结果.该方法会给成像结果带来一种尾波高频干扰：地下的震源发出的上行波与上方介质作用后，所产生的多级散射波会干扰反传波场.对此，在成像过程中，对入射场和散射场都进行左右行波分离，以压制该噪声.而在子波信息未知，无法重构入射场时，使用了激发时间成像条件，也能够实现同等效果的偏移成像，且不会出现尾波高频干扰.利用数值实验验证了本文方法的有效性. 相似文献