共查询到19条相似文献，搜索用时 127 毫秒
1.
相比地面地震勘探,垂直地震剖面(VSP)资料信噪比高,从VSP偏移结果中能够准确地识别井旁小的地质构造.在VSP偏移方法中,VSP逆时偏移方法基于双程波动方程,具有成像精确且无倾角限制、适应任意复杂速度模型等诸多优点,受到人们广泛关注.传统的VSP逆时偏移方法会在成像剖面上产生低频噪声,同时井中检波点分布位置附近会产生井轨迹噪声,严重影响了成像剖面的质量.在本文中,我们利用优化时空域频散关系求解差分系数,提高数值模拟精度;采用混合吸收边界条件压制边界反射,减少边界反射对成像的影响;通过Hilbert变换来构建复数波场,利用波场分解成像条件进行VSP逆时偏移成像,来压制低频噪声和井轨迹噪声.数值模型和实际资料测试表明,基于复数波场分解成像方法能够实现高精度的VSP逆时偏移,并且能够有效地压制VSP逆时偏移中的低频成像噪声和井轨迹噪声;相比地面地震逆时偏移,本文方法能够对井旁构造和盐丘下构造取得更清晰的成像结果. 相似文献
2.
相对保幅的角度域VSP逆时偏移(英文) 总被引：3，自引：1，他引：2
本文介绍了一种改进的角度域VSP逆时偏移方法。对VSP逆时偏移中的逆推公式进行了改进,为方便数值计算出相对保幅的角度域共成像点道集(ADCIGs)。此外VSP记录到的波场信息丰富,包括上行波场、下行波场和直达波场等,本文分析了这些波场的响应特征,发现直达波和下行波在角度域共成像点道集(ADCIGs)上都产生了成像噪音,直达波产生的噪音尤为严重。把该方法用于我国西部地区实际观测的VSP资料,不仅获得相对保幅角度域共成像点道集(ADCIGs),而且压制了成像噪音。通过数值模型试算,实际资料的应用验证了该方法的实用性与有效性,从而为VSP偏移速度分析、VSP AVA/AVO分析和反演等提供可靠的基础资料。 相似文献
3.
与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性. 相似文献
4.
静校正和近地表速度建模是地震资料处理的重要组成部分.在近地表复杂区,静校正问题和近地表速度建模问题表现尤为突出.本文利用地震干涉测量在波场重构方面的优势,给出了一种同时避免静校正和近地表速度建模过程的方法,直接由重构的SWP波场来成像高陡构造.文中首先阐述了地震干涉测量的原理,给出了VSP波场向SWP波场重构的数学表达式;接着,采用几种不同的近地表模型,对比VSP波场向SWP波场重构的效果;最后,建立局部速度模型,直接利用重构的SWP波场进行偏移处理,完成高陡构造成像.模型试验结果表明该方法的可行性,有利于实现复杂地表区深层高陡构造的成像. 相似文献
5.
当前的三维VSP地震数据偏移成像实现都是在共炮点道集或共检波点道集中逐个道集循环进行的,计算效率相对较低.根据三维VSP观测系统中炮点和检波点布置的特殊性和地震波场满足线性叠加的特性,本文提出了一种三维VSP数据的高效偏移成像方法,即首先通过对三维VSP共接收点道集进行地震数据的广义合成得到一种超道集,然后在共接收点道集的波场深度外推过程中逐步应用多震源波场对超道集进行偏移成像,即利用一次波场深度外推循环完成对所有共检波点道集数据的偏移成像.通过三维VSP模型数据与实际地震数据的偏移成像试验验证了这种高效的超道集偏移成像方法可取得与常规共检波点道集相当的偏移成像效果,还具有极高的计算效率,其计算量与单个共检波点道集的偏移成像计算量相当. 相似文献
6.
《应用地球物理》2019,(4)
本文从反射声波测井逆时偏移成像方法出发,考虑逆时偏移本身存在海量波场数据存储问题,探讨了波场存储优化策略,提出并实现了激发振幅成像条件逆时偏移和随机边界法逆时偏移两种避免波场数据存储的优化方法。首先,激发振幅成像条件方法利用计算激发时间作为成像时间判断,只需存储小部分波场,如几十个时间点的波场数据,通过实例证明甚至仅需两个时间点波场即可完成成像。传统逆时偏移通常需存储数千时间点的波场数据,与之相比数据存储和读写量可缩小几十倍甚至数千倍;其次,随机边界方法利用波场在随机介质中发生散射而非反射进行震源声场反传重构,从而避免正传波场数据的存储。数值计算表明两种方法均可有效减少反射声波逆时偏移的波场数据存储,与其它偏移算法和常规逆时偏移相比,在保证成像精度的同时提高了数据处理效率,为井旁缝洞体高效、高精度成像提供了重要的技术手段。 相似文献
7.
逆时偏移方法利用双程波算子模拟波场的正向和反向传播,通常采用互相关成像条件获得偏移剖面,是一种高精度的成像方法.但是传统的互相关成像条件会在偏移结果中产生低频噪声;此外,如果偏移速度中存在剧烈速度变化还可能进一步产生偏移假象.为了提高逆时偏移的成像质量,可在成像过程中先对震源波场和检波点波场分别进行波场分离,然后选择合适的波场成分进行互相关成像.本文基于Hilbert变换,推导了可在偏移过程中进行上下行和左右行波场分离的高效波场分离公式以及相应的成像条件,结合Sigsbee 2B合成数据,给出了不同波场成分的互相关成像结果.数值算例结果表明,采用本文提出的高效波场分离算法以及合理的波场成分互相关成像条件可以获得高信噪比的成像结果. 相似文献
8.
逆时偏移成像是基于双程波波动方程,通过对波场正向和反向延拓来实现偏移的高精度成像方法.逆时偏移能有效利用全波场信息,且相对于单程波偏移方法而言,无传播方向限制,不需要对波场进行上、下行波的分离;将波场看作矢量波的叠加,避免了倾角限制,并能较好地实现回转波、多次波成像.在给定速度场足够准确的情况下,逆时偏移理论上能对任意复杂介质精准成像,是现在最为精确的偏移方法之一.但是,低频噪声是制约其发展的因素之一.本文在分析各种成像条件优缺点的基础上,选用归一化互相关成像条件,对多个典型模型进行逆时偏移成像测试,以验证程序的正确性;选取波场分离、Poynting矢量、拉普拉斯算子滤波等三种去噪方法,分别对组合模型和Marmousi2模型进行测试,对比三种去噪方法在叠前逆时偏移成像中的去噪效果,为实际地震资料处理提供方法借鉴. 相似文献
9.
逆时偏移作为一种先进的地震偏移成像方法,其成像结果的好坏取决于很多因素,其中成像奈件是关键的一个因素.成像条件一般采用互相关成像条件,它是将源波场沿时间进行正推,接收波场沿时间进行逆推,然后将源波场和接收波场进行互相关,从而得到偏移成像结果.源波场正推和接收波场逆推都是采用同一个波动方程,这里采用波动方程一阶应力-速度形式,数值计算方法采用交错网格有限差分方法.采用振幅补偿拉普拉斯滤波方法压制逆时偏移成像中的低频噪声.采用三次样条插值方法解决逆时偏移成像中的波形不光滑问题.在接收波场逆推过程中,是将地震记录作为边界条件.而在源波场正推过程中,需要给定震源子波作为初值条件.理论上来讲,震源子波为脉冲时是最佳的,但在实际当中难以实现.在源波场正推时,一般采用子波函数.子波函数形式非常多,最常用的是雷克子波.不同主频震源子波产生的源波场频率也不同,从而导致逆时偏移互相关成像效果也不同.那么震源子波主频对逆时偏移成像影响到底有多大呢?为了简便,本文基于不同主频雷克子波,初步分析了相应的逆时偏移成像效果,同时对速度模型存在不同误差时的情形也进行了试验分析.结果 表明:在不同速度误差情况下,震源子波主频越低,成像同相轴连续性越好,散射点聚焦能量越强;震源子波主频越高,虽然分辨率提高,但成像同相轴连续性变差,散射点聚焦能量变弱,同时背景噪声增强.本文研究结果对逆时偏移成像技术在实践中的应用具有一定的参考价值. 相似文献
10.
探讨了井巷隧道工程反射波法超前地质预报成像问题.从一阶速度－应力弹性波方程出发,推导了二维各向同性介质情况下弹性波逆时传播的高阶差分格式,实现了弹性波在数值空间中的逆时延拓.构建反射波法隧道超前探测中的断层、软弱夹层等介质模型,以反射波探测的正演记录作为初始条件,并从程函方程出发,采用逆时差分格式求取介质模型网格空间中各点的直达波旅行时作为弹性波逆时偏移的成像条件,实现多波多分量资料的逆时偏移.偏移结果表明,逆时偏移能够使隧道壁接收到的波场准确归位,提高隧道反射波超前探测的资料处理的精度. 相似文献
11.
偏移成像是VSP数据处理中的一个重要环节,常规的VSP成像方法通常利用VSP-CDP转换或Kirchhoff偏移,均存在保幅性差及成像精度低等问题,而波动方程叠前深度偏移被认为是对地下复杂构造进行成像的精确偏移方法.任意广角波动方程作为一种高精度的空间域单程波波动方程,同时由于只含有二阶偏导数项,易于数值实现,与其他单程波波动方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文将AWWE推广应用到VSP数据成像中,实现了VSP时空域高角度单程波方程偏移.首先从三维标量任意广角波动方程出发,推导了完全匹配层吸收边界条件,在基本不增加计算量的前提下有效地压制了边界反射成像噪音,同时利用非线性反演算法优选参考速度来提高平方根算子的近似程度,从而提高高角度地层的成像精度.模型数值模拟实验验证了该方法的有效性,同时表明该方法在陡倾角构造情况下能取得很好的成像效果.最后对某地区实际观测的VSP资料进行了偏移成像,并与地面地震偏移结果进行了对比,显示出VSP波动方程偏移在成像分辨率上的优势. 相似文献
12.
研究不同偏移方法的成像机理是实现复杂构造高精度成像的前提,研究了两类共炮点域的相关型叠前深度偏移成像方法：基于单程波动方程的叠前深度成像方法和基于双程波动方程的叠前深度成像方法,同时对比了它们在计算效率、数据存储量、成像精度、成像机理、速度敏感性等方面的差异及其共性。以复杂构造模型为例,采用了傅里叶有限差分法（FFD）和逆时成像法（RTM）,这两种方法实现了121个共炮点道集的叠前深度偏移成像处理。计算结果表明,当速度准确时,两种深度域波动方程成像方法均可以恢复出各个地质反射界面,其中逆时偏移对陡倾角成像效果显著,当速度存在百分比误差和随机扰动情况时,逆时成像结果要差于单程波方法,因此,逆时偏移方法对速度的敏感性较大,且低频噪声较为严重。 相似文献
13.
VSP资料上下行波场发育丰富.本文在分析VSP直达波、上行反射波、下行反射波传播路径及其照明范围的基础上,指出了常规VSP波动方程偏移方法缺陷,进而通过修改波场延拓方式,提出了上下行反射波联合成像方法,并在高频近似下分析了该方法的成像原理.该方法不需要进行VSP上下行反射波场分离,能够同时对VSP资料中的一次反射波、自由表面多次波、层间多次波进行成像,比常规成像剖面具有更宽的成像范围和更好的成像效果.该方法能够对下行一次反射波进行成像,从而可以实现常规偏移方法难以处理的高陡倾角构造成像.模拟资料和实际资料处理证明了本文方法的正确性. 相似文献
14.
逆时偏移成像(RTM)常用来处理复杂速度模型,包括陡倾角及横向速度变化剧烈的模型.与常规偏移成像方法(如Kirchhoff偏移)相比,逆时偏移成像能提供更好的偏移成像结果,近些年逆时偏移成像越来越广泛地应用到勘探地震中,它逐渐成为石油地震勘探中的一种行业标准.电磁波和弹性波在动力学和运动学上存在相似性,故本文开发了基于麦克斯韦方程组的电磁波逆时偏移成像算法,并将其应用到探地雷达数据处理中.时间域有限差分(FDTD)用于模拟电磁波正向和逆向传播过程,互相关成像条件用于获得最终偏移结果.逆时偏移成像算法中,偏移成像结果受初始模型影响较大,而其中决定电磁波传播速度的介电常数的影响尤为重要.本文基于时间域全波形反演(FWI)算法反演获得了更为精确的地下介电常数模型,并将其反演结果作为逆时偏移成像的初始介电常数模型.为了验证此算法的有效性,首先构建了一个复杂地质结构模型,合成了共偏移距及共炮点探地雷达数据,分别应用常规Kirchhoff偏移算法及逆时偏移成像算法进行偏移处理,成像结果显示由逆时偏移成像算法得到的偏移结果与实际模型具有较高的一致性;此外本文在室内沙槽中进行了相关的物理模拟实验,采集了共偏移距及共炮点探地雷达数据,分别应用Kirchhoff和叠前逆时偏移成像算法进行处理,结果表明叠前逆时偏移成像在实际应用中能获得更好的成像效果. 相似文献
15.
基于双程波动方程的逆时偏移被认为是目前最好的偏移成像技术,更适合于复杂构造成像.然而,大计算量和大存储量使得逆时偏移的计算成本很高而无法用于大数据量的地震成像.本文分析了目前常用的存储策略,并分别在空间和时间上对存储策略进行了研究:空间上,根据有限差分格式,在边界存储策略的基础上通过修改波场逆向传播的边界条件,提出了有效边界存储策略.该策略可在不增加任何计算量的情况下大幅降低逆时偏移对存储量的需求;在时间上,使用checkpointing技术对有效边界存储策略进行了改进,使叠前逆时偏移在增加少量计算量的情况下进一步降低存储量需求.Marmousi模型测试结果表明了有效边界存储策略的有效性和优越性. 相似文献
16.
Quantitative imaging of complex structures from dense wide-aperture seismic data by multiscale traveltime and waveform inversions: a case study 总被引：3，自引：0，他引：3
S. Operto C. Ravaut L. Improta J. Virieux A. Herrero P. Dell'Aversana 《Geophysical Prospecting》2004,52(6):625-651
An integrated multiscale seismic imaging flow is applied to dense onshore wide‐aperture seismic data recorded in a complex geological setting (thrust belt). An initial P‐wave velocity macromodel is first developed by first‐arrival traveltime tomography. This model is used as an initial guess for subsequent full‐waveform tomography, which leads to greatly improved spatial resolution of the P‐wave velocity model. However, the application of full‐waveform tomography to the high‐frequency part of the source bandwidth is difficult, due to the non‐linearity of this kind of method. Moreover, it is computationally expensive at high frequencies since a finite‐difference method is used to model the wave propagation. Hence, full‐waveform tomography was complemented by asymptotic prestack depth migration to process the full‐source bandwidth and develop a sharp image of the short wavelengths. The final traveltime tomography model and two smoothed versions of the final full‐waveform tomography model were used as a macromodel for the prestack depth migration. In this study, wide‐aperture multifold seismic data are used. After specific preprocessing of the data, 16 frequency components ranging from 5.4 Hz to 20 Hz were inverted in cascade by the full‐waveform tomography algorithm. The full‐waveform tomography successfully imaged SW‐dipping structures previously identified as high‐resistivity bodies. The relevance of the full‐waveform tomography models is demonstrated locally by comparison with a coincident vertical seismic profiling (VSP) log available on the profile. The prestack depth‐migrated images, inferred from the traveltime, and the smoothed full‐waveform tomography macromodels are shown to be, on the whole, consistent with the final full‐waveform tomography model. A more detailed analysis, based on common‐image gather computations, and local comparison with the VSP log revealed that the most accurate migrated sections are those obtained from the full‐waveform tomography macromodels. A resolution analysis suggests that the asymptotic prestack depth migration successfully migrated the wide‐aperture components of the data, allowing medium wavelengths in addition to the short wavelengths of the structure to be imaged. The processing flow that we applied to dense wide‐aperture seismic data is shown to provide a promising approach, complementary to more classical seismic reflection data processing, to quantitative imaging of complex geological structures. 相似文献
17.
本文深入分析了多次波逆时偏移的成像原理和串扰假象产生机制,并提出了多次波分阶逆时偏移策略,即首先应用自由界面多次波衰减方法对原始炮集记录中的多次波进行剔除得到一次波记录,然后应用一次波记录预测获得各阶多次波记录,最后将各阶多次波记录分别进行逆时偏移成像.模型实验结果表明,多次波分阶逆时偏移其各阶多次波成像剖面均可对真实界面正确成像且其能够有效压制串扰假象,其中一阶多次波逆时偏移剖面的成像精度最高,其深层成像质量明显优于常规多次波逆时偏移. 相似文献
18.
传统有限差分方法在处理起伏地表时存在一些困难,而坐标变换法可将起伏地表映射为水平地表以克服此缺点.但同时,地下构造被变换得更加复杂,导致了波传播和成像的不准确.本文提出了一种分层的坐标变换方法,并应用到了弹性波逆时偏移中,此方法既可以克服起伏地表的影响,又可以不破坏地下构造.波场正向延拓、逆时延拓和分离是在辅助坐标系下完成的,而成像是在笛卡尔坐标系下完成的.通过对简单起伏模型和中原起伏模型的试算证明了本文提出方法的准确性.同时,对两种极端起伏地层高程不准确的情况进行测试可以看出:分层坐标变换逆时偏移方法的成像效果远好于传统坐标变换方法. 相似文献
19.
弹性波矢量波场逆时偏移可以综合利用纵横波场信息，对地下空间进行清晰成像，且对于成像介质没有角度限制，可以对复杂构造进行更清晰的成像.而弹性波逆时偏移中最重要的就是求解波动方程的算法，其直接影响成像的精度以及效率.本文引进电力系统谐波分析中常用的余弦组合窗函数，并通过一种新的优化算法得到了改进的余弦组合窗函数从而得到优化后的有限差分算子.并将此算子应用于解耦方程的矢量波场分离算法从而提高了成像精度.数值测试表明基于新算法的逆时偏移的成像精度和清晰度得到了明显的提高. 相似文献