共查询到20条相似文献，搜索用时 156 毫秒
1.
地震数据中发育的层间多次波是影响速度分析和偏移成像的精度和可靠性的关键.通常情况下，层间多次波的动校正量、叠加速度和频率与一次波并无明显差异，从而对识别、预测和压制多次波带来了极大挑战.传统虚同相轴方法基于物理图像和定性公式，其预测的层间多次波振幅和相位精度难以满足实际需求，造成了其对匹配算法的过度依赖.本文针对传统虚同相轴方法的理论缺陷和计算精度问题，通过理论推导得到了新的自适应虚同相轴方法.相比于传统方法，自适应虚同相轴方法能够显著提高压制多次波能力，同时减少对匹配算法的依赖.本文给出了自适应虚同相轴方法的推导过程，并运用一维和二维模型算例验证了方法相较于传统虚同相轴方法的多次波预测精度优势.通过在PLUTO模型和实际陆地地震数据上的应用实例，证明了本文新研究的自适应虚同相轴方法对去除层间多次波，恢复并突出目标储层同相轴，提高地震成像分辨率的显著作用. 相似文献
2.
3.
利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高. 相似文献
4.
《应用地球物理》2019,(4)
传统的Radon变换是利用一次波与多次波的速度差异将其聚焦在Radon域内的不同"点"或"直线"上。然而有限的"偏移孔径"、"离散采样"和"地震数据的AVO特性"会使Radon变换的聚焦变差,降低了该方法的分辨率,影响多次波压制精度。此外传统Radon变换没有考虑地震数据的AVO特性,使用L_0范数的近似形式L1范数来提高Radon域的聚焦特性,计算量极大。本文就上述问题,本文将SL_0范数和正交多项式同时引入λ-f域Radon变换,提出基于SL_0范数的高阶高分辨率λ-f域Radon变换。在考虑地震数据的AVO特性基础上,将正交多项式和Radon变换相结合。首先推导出时间域高阶Radon正反变换形式,随后使用傅里叶变换将其变换到频率域,再引入变量,将频率和曲率相结合,对其进行整体采样,具备了较快的计算效率,最后使用平滑高斯函数替代L_1范数来近似L_0范数(SL_0范数),通过最速下降和梯度投影原理求得Radon变换在SL_0范数下的稀疏解。理论模型和实际资料试算表明,同基于L_1范数的Radon变换比较,本文方法既在Radon域内具有更好的聚焦特性,同时在有效压制多次波(NMO)后存在剩余时差)的同时能更好地保存一次波的AVO特性。 相似文献
5.
在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义. 相似文献
6.
老爷庙地区位于渤海湾盆地南堡凹陷的西北部,是冀东油田的重要油气分布区,高质量的地震资料是该区油气勘探的重要依据.但是,由于受到中浅层的复杂断裂、中深层发育的三套火成岩和强横向非均质性的影响,采集的地震资料信噪比低,多次波发育.为了提高探区资料的成像品质,采用了针对性的处理策略,本文详细的分析了每种策略的基本原理并展示了应用效果,首先通过近地表Q补偿技术改善由于地表吸收造成的振幅、频率损失及相位畸变;其次,基于稳健地表一致性反褶积技术较好地解决了低信噪比资料的高分辨率处理问题,避免传统反褶积方法引起的噪声能量加强;再次,应用可适应非规则观测系统的广义多次波衰减技术(GSMP)压制了部分强火成岩引起的多次波.针对性处理对策克服了常规处理的不足,提高了该地区地震资料的成像品质,获得较高分辨率和信噪比的成果数据,为老爷庙地区油气勘探奠定了良好的基础. 相似文献
7.
8.
《地球物理学进展》2015,(1)
3D VSP地震数据具有高分辨率和高信噪比的特点,常规的VSP成像方法在大庆油田资料的处理中一直未见到理想的效果.依据大庆油田VSP资料的特点,设计了3D VSP资料的处理流程,主要包括三分量波场逐点定向方法和Q补偿的子波反褶积等预处理方法,所设计的流程具有高分辨率、高保真性的特点.大庆油田M1井的实际资料处理结果显示,所形成的方法可将VSP资料的频带展宽约20 Hz,同时对地面地震也进行了反褶积处理,反褶积后,地面地震的分辨率得到有效提高,波组特征更加清楚.实际数据计算表明,所形成的方法具有很好的实用性,可实现地下介质的精细刻画,可在具有同类地质条件的研究区推广应用. 相似文献
9.
多次波偏移中的假象主要来自于不同地震事件之间的互相关,由于这种互相关满足成像条件,很难直接在偏移过程中去除.但是对于准确的速度模型,真实的成像结果在角度域内应该是平直的.根据这个判断准则,可以在角度域内移除多次波偏移中的假象.本文以数据自相关偏移为例,提出了在单程波多次波偏移中移除假象的主要流程:首先在在单程波偏移过程中高效地提取角度域共成像点道集,然后对角度域共成像点道集应用高分辨率的抛物线型Radon变换,用合适的切除函数处理后,反变换回到角度域,最后叠加各个角度成分,得到偏移结果.Marmousi模型的合成数据测试表明,这种方法可以很好地压制多次波偏移过程中产生的假象,有效地提高成像结果的信噪比. 相似文献
10.
基于波动方程预测的表面多次波压制方法可处理复杂地下介质的地震资料,但计算成本较高.基于滤波的多次波压制方法计算效率较高,但其成功应用仅局限于一次波和多次波有明显时差差别的地震数据,对来自速度逆转等复杂介质数据则较难获得满意的压制效果.本文将波动方程预测的反馈迭代法和滤波法有效结合,采用GPU(图形处理器)和CPU协同并行加速计算粗略预测表面多次波,随后在双曲Radon域比较分析原始数据和预测的多次波,设计合理有效的Butterworth型自适应滤波器,滤出原始数据Radon域中的多次波能量,进行Radon反变换后,在时空域将多次波从原始数据中减去,得多次波压制结果.文中对理论模拟的单炮数据、复杂的SMAART模型以及实际地震数据进行了计算,结果表明,结合基于波动方程预测和双曲Radon变换的方法有效突破了两种方法各自的局限性,可高效高精度地压制复杂地下介质的表面多次波. 相似文献
11.
将多次波转换成反射波并按传统反射波偏移算法成像,是多次波成像的一种方法.聚焦变换能准确的将多次波转换为纵向分辨率更高的新波场记录,其中一阶多次波转换为反射波.本文对聚焦变换提出了两点改进:1)提出局部聚焦变换,以减小存储量和计算量,增强该方法对检波点随炮点移动的采集数据的适应性;2)引入加权矩阵,理论上证明原始记录的炮点比检波点稀疏时,共检波点道集域的局部聚焦变换可以将多次波准确转换成炮点与检波点有相同采样频率的新波场记录.本文在第一个数值实验中对比了对包含反射波与多次波的原始记录做局部聚焦变换和直接对预测的多次波做局部聚焦变换两种方案,验证了第二种方案转换得到的波场记录信噪比更高且避免了第一个方案中切聚焦点这项比较繁杂的工作.第二个数值实验表明:在炮点采样较为稀疏时,该方法能有效的将一阶多次波转换成反射波;转换的反射波能提供更丰富的波场信息,成像结果更均衡、在局部有更高的信噪比,以及较高的纵向分辨率. 相似文献
12.
提出一种基于Radon 变换和正交多项式变换的多方向正交多项式变换压制多次波方法.抛物Radon变换对不同曲率方向的同相轴叠加,根据速度差异区分一次波和多次波,但Radon反变换会损伤振幅特性,不利于AVO分析.多方向正交多项式变换在Radon变换(某一曲率方向的零阶特性)的基础上,利用正交多项式变换进一步分析同相轴的高阶多项式分布特性,用正交多项式谱表征同相轴AVO特性;根据一次波和多次波速度差异和同相轴能量分布特征实现多次波压制.该方法的优点是仅用一个曲率参数就可描述同相轴剩余时差参数,提高了一次波和多次波的剩余时差分辨率.实验结果表明,该方法可以有效压制多次波并保留一次波AVO特性. 相似文献
13.
一种改进的基于非高斯性最大化的预测反褶积算法 总被引：3，自引：1，他引：2
The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results. 相似文献
14.
The parabolic Radon transform has been widely used in multiple attenuation. To further improve the accuracy and efficiency of the Radon transform, we developed the 2- fdomain high-resolution Radon transform based on the fast and modified parabolic Radon transform presented by Abbad. The introduction of a new variable 2 makes the transform operator frequency-independent. Thus, we need to calculate the transform operator and its inverse operator only once, which greatly improves the computational efficiency. Besides, because the primaries and multiples are distributed on straight lines with different slopes in the 2-fdomain, we can easily choose the filtering operator to suppress the multiples. At the same time, the proposed method offers the advantage of high-resolution Radon transform, which can greatly improve the precision of attenuating the multiples. Numerical experiments suggest that the multiples are well suppressed and the amplitude versus offset characteristics of the primaries are well maintained. Real data processing results further verify the effectiveness and feasibility of the method. 相似文献
15.
本文深入分析了多次波逆时偏移的成像原理和串扰假象产生机制,并提出了多次波分阶逆时偏移策略,即首先应用自由界面多次波衰减方法对原始炮集记录中的多次波进行剔除得到一次波记录,然后应用一次波记录预测获得各阶多次波记录,最后将各阶多次波记录分别进行逆时偏移成像.模型实验结果表明,多次波分阶逆时偏移其各阶多次波成像剖面均可对真实界面正确成像且其能够有效压制串扰假象,其中一阶多次波逆时偏移剖面的成像精度最高,其深层成像质量明显优于常规多次波逆时偏移. 相似文献
16.
Multiples have longer propagation paths and smaller reflection angles than primaries for the same source–receiver combination, so they cover a larger illumination area. Therefore, multiples can be used to image shadow zones of primaries. Least-squares reverse-time migration of multiples can produce high-quality images with fewer artefacts, high resolution and balanced amplitudes. However, viscoelasticity exists widely in the earth, especially in the deep-sea environment, and the influence of Q attenuation on multiples is much more serious than primaries due to multiples have longer paths. To compensate for Q attenuation of multiples, Q-compensated least-squares reverse-time migration of different-order multiples is proposed by deriving viscoacoustic Born modelling operators, adjoint operators and demigration operators for different-order multiples. Based on inversion theory, this method compensates for Q attenuation along all the propagation paths of multiples. Examples of a simple four-layer model, a modified attenuating Sigsbee2B model and a field data set suggest that the proposed method can produce better imaging results than Q-compensated least-squares reverse-time migration of primaries and regular least-squares reverse-time migration of multiples. 相似文献
17.
Reverse time migration of multiples can be used to construct subsurface structures where primaries cannot illuminate well. However, the images generated using multiples suffer from severe artefacts due to the cross-talks created by interference among unrelated multiples. We developed a migration approach using water-bottom-related multiples to reduce these cross-talk artefacts. This approach first isolates primaries from the original data and predicts water-column primaries. The nth-order water-column multiples can be obtained by auto-convolving the water-column primaries n times, followed by convolving the nth-order water-column multiples with the primaries to extract the (n+1)th-order water-bottom-related multiples. The approach takes the nth-order water-column multiples as the secondary source and regards the (n+1)th-order water-bottom-related multiples as the receiver wavefield, followed by a cross-correlation imaging condition. Numerical examples from synthetic and field data sets reveal that our approach can provide images with substantially fewer cross-talk artefacts than conventional reverse time migration using multiples, as well as greatly improving shallow imaging compared with reverse time migration of primaries. 相似文献
18.
Short-period multiple reflections pose a particular problem in the North Sea where predictive deconvolution is often only partially successful. The targeted multiple attenuation (TMA) algorithm comprises computation of the covariance matrix of preflattened prestack or post-stack seismic data, the determination of the dominating eigenvectors of the covariance matrix, and subtraction of the related eigenimages followed by reverse flattening. The main assumption made is that the flattened multiple reflections may be represented by the first eigenimage(s) which implies that the spatial amplitude variations of primaries and associated multiples are similar. This assumption usually limits the method to short-period multiple reflections. TMA is applicable post-stack or prestack to common-offset gathers. It is computationally fast, robust towards random noise, irregular geometry and spatial aliasing, and it preserves the amplitudes of primaries provided they are not parallel to the targeted multiples. Application of TMA to 3D wavefields is preferable because this allows a better discrimination between primaries and multiples. Real data examples show that the danger of partially removing primary energy can be reduced by improving the raw multiple model that is based on eigenimages, for example by prediction filtering. 相似文献
19.
提出了以二进小波变换为基础的自适应Kalman滤波反褶积(AKFD)新方法，针对该方法的计算复杂程度，提出了一种快速实现方法.二进小波变换的AKFD抛弃了传统预测反褶积对信号平稳性的假设，克服了提高分辨率而信噪比明显降低的问题，具有很好的抗噪性能.在小波域进行的AKFD在压制假反射以及提高分辨率方面比时间域的AKFD好，克服了在时域内进行AKFD抬升低频成分的缺陷.利用二维地震数据的局部平稳性的假设提出了快速实现方法，通过分段求取自适应预测算子，分别于横向及纵向采用样条插值的方法进行插值，来减少求取自适应预测算子的计算量，达到快速实现的目的.经过大量实验表明计算速度提高数百倍，仍能保持原来的计算效果. 相似文献
20.
The key processes in marine seismic imaging include (i) removing from seismic data all seismic events (free-surface multiples and ghosts) which contain at least one reflection at the sea surface in their wave-propagation path, and leaving those with no reflection at the free surface (internal multiples and primaries), (ii) removing events with at least two reflections in the subsurface (internal multiples), and leaving events with only one reflection in the subsurface (primaries), and then (iii) locating the scattering points and reflectors inside the subsurface which are the sources of primaries and internal multiple events. All these processes are here explained, derived, and optimized via scattering diagrams (diagrammatica) in a way similar to the way the quantum field theory is often explained via Feynman diagrams. Our discussion of the removal of events with free-surface reflections from the data will be brief, as the diagrammatica of these events are now well understood.The main focus of this paper is the diagrammatica of internal multiples and primaries. Although these events do not contain any reflection at the sea surface, it is important to reconstruct them with scattering points near the sea surface, where seismic data are recorded. So our diagrammatica of primaries and internal multiples include events which are not directly recorded in seismic data but which can be constructed from seismic data. These events have allowed us to construct scattering diagrams of primaries and internal multiples with scattering points near the sea surface. Furthermore, these new diagrammatica of internal multiples and primaries can be used to remove internal multiples from the data. 相似文献