首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 171 毫秒
1.
基于单程波方程的角度域保幅偏移(英文)   总被引:8,自引:6,他引:2  
传统叠前深度偏移只能够提供地下的构造信息,但工业界在需要构造信息的同时还要与地下界面反射系数成比例的振幅信息。最近几年,基于单程波方程的保幅叠前深度偏移算法有了一定的发展,但是,基于炮域、单程波的保幅型叠前深度偏移必须应用反褶积型的成像条件,这种成像条件在构造复杂、速度变化剧烈的地区会出现不稳定现象。基于角度域的保幅深度偏移克服了这一不稳定性缺点的同时,还域的保幅深度偏移,模型和实际资料的试算分析验证该思路方法的正确性和有效性。  相似文献   

2.
特征高斯波包叠前深度偏移方法   总被引:3,自引:3,他引:0       下载免费PDF全文
李辉  王华忠  冯波  胡英  张才 《地球物理学报》2014,57(7):2258-2268
高斯波包(Gaussian packet)传播算子可在局部时空域高效地计算局部波包的传播.高斯波包叠前深度偏移的基础是在Gabor变换域描述观测数据,再利用高斯波包传播算子计算炮点波场和检波点波场,两者相关即可得到偏移结果.利用炮道集的局部τ-p特征在Gabor变换域表达观测数据,可以仅关注部分高斯波包框架函数上的数据投影,这样既实现了波场的压缩存储,同时可利用高斯波包传播算子反传框架函数以实现整个炮道集的快速反传.这些综合了观测数据局部τ-p特征的高斯波包函数称为特征高斯波包(characteristic Gaussian packet,CGP),相应的波场反传称为特征高斯波包反传.理论及数值分析证明了上述特征高斯波包反传方法是有效且快速的.炮点正传波场也利用高斯波包传播算子模拟.利用互相关成像条件可实现特征高斯波包叠前深度偏移(characteristic Gaussian packet pre-stack depth migration,CGPM).由于高斯波包传播算子描述了局部方向及局部空间的波场,所以CGPM可以自然地提取角度域成像道集(ADCIG),并易于实现面向目标叠前深度偏移,从而作为偏移引擎为偏移速度分析(MVA)服务.数值实验证明了CGPM和面向目标CGPM的有效性和实用性.  相似文献   

3.
多分量地震资料叠前深度偏移技术可以对地下复杂地质构造进行更准确的成像,精确成像的前提是获取准确的纵横波偏移速度。本文采用高斯束偏移方法对多波地震数据进行偏移速度分析,首先分别给出纵波和转换波共偏移距域高斯束叠前深度偏移方法原理,在此基础上抽取纵波和转换波偏移距域共成像点道集;然后根据共成像点道集拉平准则,分别对纵波和横波速度进行更新;当两种波成像深度不一致时,对纵波和转换波成像剖面进行深度匹配,完成高精度的纵横波偏移速度分析。模型数据和实际资料试算表明,该方法是一种有效的多波偏移速度分析方法。  相似文献   

4.
Kirchhoff型真振幅反偏移的基本目的是通过等时线叠加得到偏移的输入场,是将深度偏移的结果再反映射到时间域去,而真振幅加权函数是Kirchhoff型真振幅反偏移的核心内容,它没有唯一的数学形式,目前发表在文献当中的加权函数公式是在物理直观下的最优表述.本文在前人研究的基础上推导了真振幅反偏移加权函数公式,对影响加权...  相似文献   

5.
τ-P域地震资料三维叠前深度偏移及其实用算法   总被引:3,自引:2,他引:3       下载免费PDF全文
本文提出了在τ-P域实现三维叠前深度偏移的计算方法,该方法把地震道集由时-空域变换到τ-P域,然后对共P数据体进行相移或相移加插值偏移,将偏移后的各共P数据体叠加得到叠前深度偏移结果.该算法能够实现陡倾地层的正确归位,降低偏移处理维数,减少计算机内存需求,易于并行处理,但是计算量和所花费的计算机时间仍然很大.为此我们又提出了针对目的层的分时(层)偏移方法,由于只考虑目的层,可节省大量的计算时间.文中以二维实际地震资料及三维数值模型资料的处理说明了上述方法的可行性与有效性.  相似文献   

6.
双复杂条件下的波动方程叠前深度偏移   总被引:1,自引:1,他引:0       下载免费PDF全文
复杂地表条件下的地震勘探越来越被人们所关注.双复杂条件下的叠前深度偏移方法是解决复杂地表条件和复杂地质构造成像的有效手段."波场上延"法能实现由非水平观测界面开始的偏移过程,解决复杂地表对地下构造成像的影响.复杂理论模型的试算以及实际资料处理表明,"波场上延"方法较好地克服了起伏地形对地下构造成像的影响,取得了令人满意的效果,实现了波动方程基准面校正和深度成像的有机结合.  相似文献   

7.
二维大地电磁资料的偏移   总被引:3,自引:2,他引:3       下载免费PDF全文
本文采用反射映像(即U/D成像)原理对大地电磁资料进行成像,将地面观测的波场分解成上行波和下行波,并分别向下延拓,利用上、下行波的时间一致性原理确定地下电性界面的位置.对TE和TM模式的MT响应分别处理可以得到两个深度剖面,两剖面的一致性便为实际反射界面的最佳估计.理论和实际资料的计算表明,MT偏移技术是比较有效的,它可获得地下界面直观的映像.和常规MT反演方法相比较,该方法具有算法简单等优点,并且还能获得表示地下地质构造真实映像的MT深度剖面.  相似文献   

8.
利用τ-p波场反演法确定东沙群岛的地壳速度模型   总被引:4,自引:0,他引:4       下载免费PDF全文
本文利用τ-p波场反演法直接对OBS原始数据进行转换,得到一维速度深度函数.该方法包含倾斜迭加和向下延拓两个线性转换:倾斜迭加将原始数据由时间-偏移距域(t-x)转换为时间截距-慢度域(τ-p);向下延拓将时间截距-慢度域(τ-p)波场转换为深度-慢度域(z-p)波场.然后根据五个OBS台站得到的一维速度深度函数,合成东沙群岛地壳的二维速度结构模型.为了检验该速度结构模型,我们用拾取走时和利用模型计算的理论走时进行拟合,得到了满意的结果,从而为下一步进行层析成像提供了可靠的初始速度结构模型.  相似文献   

9.
最小二乘傅立叶有限差分偏移   总被引:6,自引:1,他引:5       下载免费PDF全文
一般偏移算法是用反演算子通过解析方法求解.最小二乘偏移方法采用另一种思路,即采用数值方法,通过解一个线性离散反问题来索求解.这样我们试着寻找一个模型匹配地震数据并能表现出其某些特点来得到偏移图像.最小二乘法能减少偏移赝像,得到更精确的偏移效果.Kirchhoff算子在最小二乘偏移方法中应用较广,但需要较多的迭代次数,而且具有Kirchhoff偏移的缺点.本文把最小二乘方法运用到基于波长延拓的波动方程偏移方法中,为提高最小二乘偏移的效率,可采用效率较高的正传播算子和反传播算子.我们利用效率较高,能适应剧烈横向变速的傅立叶有限差分正传播和反传播算子来做叠后最小二乘偏移.数值实例表明,通过少数的共轭梯度法迭代,就能得到与真实模型差别很微小的偏移效果.对于傅式变换我们采用了数值软件FFTW,其变换速度比常规FFT算法一般要快六倍以上,进一步提高了效率.本文算法很容易在并行机上实现,这些特点在处理大型数据时大有裨益.  相似文献   

10.
多次波偏移中的假象主要来自于不同地震事件之间的互相关,由于这种互相关满足成像条件,很难直接在偏移过程中去除.但是对于准确的速度模型,真实的成像结果在角度域内应该是平直的.根据这个判断准则,可以在角度域内移除多次波偏移中的假象.本文以数据自相关偏移为例,提出了在单程波多次波偏移中移除假象的主要流程:首先在在单程波偏移过程中高效地提取角度域共成像点道集,然后对角度域共成像点道集应用高分辨率的抛物线型Radon变换,用合适的切除函数处理后,反变换回到角度域,最后叠加各个角度成分,得到偏移结果.Marmousi模型的合成数据测试表明,这种方法可以很好地压制多次波偏移过程中产生的假象,有效地提高成像结果的信噪比.  相似文献   

11.
线性同相轴波场分离的高分辨率τ-p变换法   总被引:7,自引:5,他引:7       下载免费PDF全文
基于最小二乘τ-p变换和τ-p域模型稀疏分布的假设,本文给出高分辨率τ-p变换的推导及其模型空间域的离散采样公式,同时给出了保振幅线性同相轴波场分离的算法流程.在求解本文给出的高分辨率τ-p正变换时,由于待求解的矩阵不具备最小平方法所具有的Toeplitz结构,故采用Cholesky分解法进行计算.本文模拟了井间地震和阵列声波测井中的Stoneley上下行波的分离算法过程,高分辨率正反τ-p变换且滤波所得结果显示本文算法误差小和保振幅的特点.对于在τ-p域距离很近或时间域同相轴近于水平的线性波场,高分辨率算法的聚焦作用使得所分离波场畸变小,体现本文算法精度高的优点.理论模型试算表明本文给出的高分辨率τ-p变换线性波场分离算法具有稳定性、精度高和保振幅的特点.  相似文献   

12.
多尺度τ-p谱及其应用   总被引:2,自引:1,他引:1       下载免费PDF全文
把小波变换与τ-p变换有机地结合起来,本文提出了多尺度τ-p谱的概念,给出了基于多尺度τ-p谱作滤波处理和波场分离的方法.与通常的τ-p变换相比,该方法将x-t域的地震记录变换到(τ,p,a)或(τ,p,f)空间(即多尺度τ-p谱).多尺度τ-p谱比通常τ-p谱增加了一维,因此用于波场分离或去噪,前者优于后者.若恰当地选择小波函数,该方法抗噪能力强,计算精度高.文中给出了模型及实际地震资料处理实例;证明了此方法的有效性.  相似文献   

13.
共炮检距道集波动方程保幅叠前深度偏移方法   总被引:9,自引:9,他引:0       下载免费PDF全文
本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括:(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率-波数域和频率-空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.  相似文献   

14.
Iterative migration velocity analysis is computationally expensive, where most of the computation time is used for generating prestack depth images. By using a reduced form of Kirchhoff migration, denoted as wave path migration, we can significantly speed up the depth imaging process and reduce the entire velocity analysis expense accordingly. Our results with 2D synthetic and field data show that wave path migration velocity analysis can efficiently improve the velocity model and the wave path migration velocity analysis updated velocity correlates well with that from the Kirchhoff migration velocity analysis. The central processing unit comparison shows that, for a 2D synthetic and field data set, wave path migration velocity analysis is six times faster than Kirchhoff migration velocity analysis. This efficiency should be even greater for 3D data.  相似文献   

15.
从波动方程出发,推导出平面波场传播基本公式,阐述了利用地表观测波场反演地球内部波速结构的理论关系.观测波场可通过-p 变换分解为地表平面波场,其最大振幅的轨迹能稳定地反映出地球内部波速随深度变化的趋势,这一特征可用来对反演解空间进行约束.波场延拓技术可以充分地利用观测波场中包含的丰富信息,在迭代反演中不仅能简单快速地得到反演解,而且所得解具有良好的稳定性,较少受主观因素影响,这是一种很好的反演技术.对波场延拓进行了理论分析与数值模拟,采用了同态反褶积等改善资料信噪比的措施,使解的分辨能力得到了提高.文中对南海北部一个声纳折射剖面进行了分析和计算.结果表明:该区1.4km 深处,存在一个从1.76km/s 到2.21km/s 的速度间断面.间断面上下两层的速度梯度分别为0.54kms-1/km,0.63kms-1/km.最后,从构造演化的角度对浅海构造特征进行了探讨.   相似文献   

16.
适于复杂介质的高精度波场延拓算子是叠前深度偏移研究的重要内容。本文采用最优可分表示方法,运用正反傅立叶变换构造了三维单程波场延拓算子,算子实现了波数域变量与空间(速度)域变量分离。波数域内进行相移计算,在空间域对因介质横向变速引起的时移作修正。脉冲响应显示在区域内各速度的脉冲计算值与理论值基本一致,说明最优可分表示法叠前深度偏移可适用于强变速条件下复杂介质的成像需求。SEG/EAGE模型和实测数据的成像结果验证了本文方法对复杂构造的成像能力。  相似文献   

17.
An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure.  相似文献   

18.
基于控制照明的合成震源记录交互剩余偏移速度分析   总被引:10,自引:0,他引:10       下载免费PDF全文
提出了一种新的偏移速度分析方法——基于控制照明的合成震源记录交互剩余偏移速度分析方法.与其他类似偏移速度分析方法的不同点在于:(1) 叠前深度偏移采用基于波动理论的快速合成震源记录算法;(2)偏移方法采用平面波震源,与速度分析方法一致;(3)应用控制照明技术,避免了因横向变速而导致的平面波震源波场在传播过程中的畸变,从而减小了速度分析的误差;(4)实用的速度谱设计,使交互偏移速度分析可行且易于操作.模型和新疆实际资料的试算表明该方法是一种有效和实用的偏移速度分析方法.  相似文献   

19.
Depth velocity model building remains a difficult step within the seismic depth imaging sequence. Stereotomography provides an efficient solution to this problem but was limited until now to a picking of seismic data in the prestack time un-migrated domain. We propose here a method for stereotomographic data picking in the depth migrated domain. Picking in the depth migrated domain exhibits the advantage of a better signal-to-noise ratio and of a more regular distribution of picked events in the model, leading to a better constrained tomographic inverse problem. Moreover, any improvement on the velocity model will improve the migrated results, again leading to improved picking. Our strategy for obtaining a stereotomographic dataset from a prestack depth migration is based on migration of attributes (and not on a kinematic demigration approach!). For any locally coherent event in the migrated image, migration of attributes allows one to compute ray parameter attributes corresponding to the specular reflection angle and dip. For application to stereotomography, the necessary attributes are the source/receiver locations, the traveltime and the data slopes. For the data slope, when the migration velocity model is erroneous, some additional corrections have to be applied to the result of migration of the attributes. Applying these corrections, our picking method is theoretically valid whatever the quality of the migration velocity model. We first present the theoretical aspects of the method and then validate it on 2D synthetic and real seismic reflection data sets.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号