首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 171 毫秒
1.
将局部余弦基小波束波场分解、传播与观测系统沉降法叠前深度偏移相结合,推导了源-检波器观测系统沉降法传播算子.本算法中,先对频率域的共点源和共点检波器道集做局部余弦小波束分解,然后分别沿共小波束源和共小波束检波器在深度方向延拓得到下一层波场.每个深度的波场,都等效于把源和检波器放在该层后所能接收到的地震记录,每点的像值由炮点和检波点重合时的零时刻波场值给出.通过二维SEG/EAGE盐丘模型和Marmousi模型的偏移成像结果验证该方法理论推导的正确性.另外,结果显示该方法继承了小波束域波场延拓在速度扰动较大情况下波传播及成像精度高的优点.  相似文献   

2.
Dreamlet偏移成像目的是探索一类能够对地震波场和单程波传播算子同时分解和压缩的理论和方法,也即实现在压缩域的传播与成像、地震数据在时间和空间的非平稳性质,决定了要实现地震数据的有效稀疏表示,分解方法必须在时间和空间上同时具有局域化性质.Dreamlet由时间和空间局部分解原子的张量积构成,可以看作一种脉冲-小波束形式的波场分解原子.时空局域化的dreamlet单程波传播算子在对波场沿深度方向延拓时,地震数据在时间轴上总是向同一方向流动.随着深度的增加,部分用于成像浅层结构的数据归位至其空间位置后被dreamlet算子丢弃,波场的有效记录时间变短,每一步用于波场延拓的计算量也相应下降.为了充分发挥这一优势,本文介绍dreamlet传播算子与观测系统沉降法偏移相结合的理论与方法.观测系统沉降法偏移每一步都将记录到的所有数据向下延拓,沉降后的波场等效于把源和检波器都放置在目标深度所能接收的反射数据.Dreamlet观测系统沉降过程只保留用于成像观测系统下部地质结构的有效数据,自动丢弃已经用于成像观测系统上部而对下部成像没有贡献的信号.本文通过二维SEG/EAGE叠后和Marmousi叠前数据算例展示了dreamlet传播算子应用于观测系统沉降法偏移的这一特点.数值算例结果显示,在不影响成像质量的前提下,该偏移方法能够有效减少传播数据量,为发展一种快速高效的偏移方法提供了新的思路.  相似文献   

3.
角度域弹性波Kirchhoff叠前深度偏移速度分析方法   总被引:3,自引:3,他引:0       下载免费PDF全文
杜启振  李芳    秦童  毕丽飞 《地球物理学报》2011,54(5):1327-1339
为提高地震成像结果的准确性并真实反映实际地震波场在介质中的传播特性,应该充分利用多分量地震数据的矢量特征进行弹性波成像,其中,最为棘手的问题是纵横波偏移速度场的确定,为此,本文提出了直接利用多分量地震数据进行弹性波角度域偏移速度分析的方法.基于空移成像条件的弹性波Kirchhoff偏移方程提取了弹性波局部偏移距域共成像...  相似文献   

4.
基于Gabor框架合成小波束叠前深度偏移   总被引:1,自引:1,他引:0  
利用小波束域波场分解和传播在空间及方向上的双重局域性,通过Fourier传播算子进行深度外推波场,完成小波束的传播,最终应用成像条件求取成像值.合成小波束和小束源记录,进行部分小束源叠前偏移成像的数值试验,获得的成像结果表明,小波束叠前深度成像是切实有效的.  相似文献   

5.
散射系数矩阵的计算及其应用   总被引:6,自引:2,他引:6       下载免费PDF全文
利用Gabor Daubechies(G? D)小波束域波场分解和传播在空间和方向上的双重局域性,提出了基于G? D小波束域叠前深度偏移进行角度域成像和计算局部散射系数矩阵的方法. 以简单分层模型为例,对不同探测系统的局部散射系数矩阵分布特征进行分析. 分析结果表明,在一定的探测系统几何布局下,由本文方法得到的局部散射系数矩阵能够较真实地反映局部结构的散射(或反射)特性. 通过局部散射系数矩阵进一步外推具有不同速度反差的水平界面随角度变化的反射系数,并估计界面的空间位置和倾角等说明局部散射系数矩阵的潜在应用.  相似文献   

6.
利用GaborDaubechies(G?D)小波束域波场分解和传播在空间和方向上的双重局域性,提出了基于G?D小波束域叠前深度偏移进行角度域成像和计算局部散射系数矩阵的方法. 以简单分层模型为例,对不同探测系统的局部散射系数矩阵分布特征进行分析. 分析结果表明,在一定的探测系统几何布局下,由本文方法得到的局部散射系数矩阵能够较真实地反映局部结构的散射(或反射)特性. 通过局部散射系数矩阵进一步外推具有不同速度反差的水平界面随角度变化的反射系数,并估计界面的空间位置和倾角等说明局部散射系数矩阵的潜在应用.  相似文献   

7.
方位角度域共成像点道集能够客观反映地下介质的速度、各向异性参数异常以及振幅随角度变化(AVA)和裂缝信息。传统Kirchhoff PSTM通常输出偏移距域共成像点道集,对于速度分析、各向异性分析、AVA分析、裂缝识别等均存在诸多不便。本文提出了基于走时梯度的Kirchhoff叠前时间偏移全方位角度集输出方法并提出工业上切实可行的实现方案。通过走时场梯度计算波场传播方向矢量,形成能够反映观测系统参数和波场传播情况的全方位角度域共成像点道集。为了在大规模地震数据Kirchhoff积分叠前时间偏移中输出全方位角度道集,本文给出基于输入道方式的偏移实现方法,采用逐条inline线进行线偏移成像,从而大大降低了全方位角度道集输出对计算机内存的压力,显著提高了Kirchhoff积分时间偏移输出全方位角度道集的可行性。三维盐丘模型测试和海上某区块三维实际资料试验证明了本文方法的正确性。   相似文献   

8.
起伏的地表条件限制了采集孔径范围并造成深层地震照明不足,为改善该类地区的成像质量,本文提出了一种起伏地表条件下的照明补偿方法.首先,基于小波束波场延拓算子和逐步累加的外推方法在波场延拓过程中解决起伏地表面的影响,并引入空间滤波函数压制虚拟层内的偏移噪音;其次,利用局部指数标架对上、下行波场分解,得到局部角度域成像和照明补偿因子.再次,利用计算出的成像值和照明补偿因子,在局部倾角域完成照明补偿.SEG起伏地表模型测试证明了本方法的有效性,深层构造照明度明显加强,不同角度成像振幅更加均衡,该技术为提高起伏地表地区的成像品质提供了新的手段.  相似文献   

9.
窄方位地震数据双平方根方程偏移方法探讨   总被引:13,自引:2,他引:11       下载免费PDF全文
简单回顾双平方根(DSR)方程偏移一般理论之后,介绍了适应偏移速度分析的角度域成像快速算法.针对地震勘探中窄方位三维地震数据成像问题,探讨了限定输入数据体或限定波场传播方向的DSR方程偏移方法.数值试验表明,共方位角偏移等DSR方程叠前深度偏移方法在盐体或礁体、古潜山等强横向变速地区具有很好的应用前景.  相似文献   

10.
时-空局域化地震波传播方法:Dreamlet叠前深度偏移   总被引:1,自引:1,他引:0       下载免费PDF全文
提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(dreamlet算子)沿深度延拓.本文利用多维局部余弦变换实现dreamlet算法,分解后的波场系数和传播算子不仅有很好的稀疏性,且均为实数,也即波的传播和成像过程完全在实数域实现.文中推导了局部余弦基dreamlet波场分解和传播算子理论公式并将其应用于叠前深度偏移.在dreamlet相空间波的传播过程为稀疏矩阵相乘,而且延拓后的地表数据波场的有效时间长度随深度的增加不断减小,从而可以减少需要传播的波场系数.二维SEG/EAGE盐丘和SIGSBEE模型算例验证了理论推导的正确性,成像结果显示该方法在横向速度变化剧烈情况下有很好的精度.  相似文献   

11.
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.  相似文献   

12.
叠前地震数据的平面波深度偏移法   总被引:7,自引:3,他引:7       下载免费PDF全文
提出了一套基于平面波分解的波动方程叠前地震数据深度偏移方法. 通过对共炮点道集和共偏移距道集地震数据的平面波分解,分别得到适用于单平方根波场外推方程和双平方根波场外推方程的共ps(炮点坐标平面波参数)平面波道集和共ph(偏移距坐标平面波参数)平面波道集. 在对共炮点道集和共偏移距道集地震数据的平面波分解时,不需要进行通常意义下的τ p变换计算. 通过对共ps平面波道集和共ph平面波道集的偏移效果对比,我们认为在速度弱横向变化介质中,两种平面波道集偏移方法的效果相当,但对于速度强横向变化介质,共ps平面波道集偏移方法的效果要优于共ph平面波道集偏移方法. 在计算效率方面,共ps平面波道集偏移方法与共ph平面波道集偏移方法基本相同.  相似文献   

13.
Reverse‐time migration has become an industry standard for imaging in complex geological areas. We present an approach for increasing its imaging resolution by employing time‐shift gathers. The method consists of two steps: (i) migrating seismic data with the extended imaging condition to get time‐shift gathers and (ii) accumulating the information from time‐shift gathers after they are transformed to zero‐lag time‐shift by a post‐stack depth migration on a finer grid. The final image is generated on a grid, which is denser than that of the original image, thus improving the resolution of the migrated images. Our method is based on the observation that non‐zero‐lag time‐shift images recorded on the regular computing grid contain the information of zero‐lag time‐shift image on a denser grid, and such information can be continued to zero‐lag time‐shift and refocused at the correct locations on the denser grid. The extra computational cost of the proposed method amounts to the computational cost of zero‐offset migration and is almost negligible compared with the cost of pre‐stack shot‐record reverse‐time migration. Numerical tests on synthetic models demonstrate that the method can effectively improve reverse‐time migration resolution. It can also be regarded as an approach to improve the efficiency of reverse‐time migration by performing wavefield extrapolation on a coarse grid and by generating the final image on the desired fine grid.  相似文献   

14.
双平方根方程三维叠前深度偏移   总被引:10,自引:6,他引:10       下载免费PDF全文
从双平方根(DSR)形式的波动方程出发,基于沉降观测概念和地震波扰动理论,介绍了深度域的DSR全偏移算子及共成像道集的生成方法. 根据三维地震数据的方位角特征,通过对全偏移算子的稳相近似,依次导出了适应于零方位角道集、Cross line共偏移距道集以及共偏移距矢量道集的偏移算子. 理论分析与合成数据的数值试验表明,DSR全偏移算子、共方位角偏移算子对介质速度变化的适应性很强,而其余两种偏移算子仅适用于缓变速情况.  相似文献   

15.
三维VSP数据高效偏移成像的超道集方法   总被引:1,自引:1,他引:0       下载免费PDF全文
当前的三维VSP地震数据偏移成像实现都是在共炮点道集或共检波点道集中逐个道集循环进行的,计算效率相对较低.根据三维VSP观测系统中炮点和检波点布置的特殊性和地震波场满足线性叠加的特性,本文提出了一种三维VSP数据的高效偏移成像方法,即首先通过对三维VSP共接收点道集进行地震数据的广义合成得到一种超道集,然后在共接收点道集的波场深度外推过程中逐步应用多震源波场对超道集进行偏移成像,即利用一次波场深度外推循环完成对所有共检波点道集数据的偏移成像.通过三维VSP模型数据与实际地震数据的偏移成像试验验证了这种高效的超道集偏移成像方法可取得与常规共检波点道集相当的偏移成像效果,还具有极高的计算效率,其计算量与单个共检波点道集的偏移成像计算量相当.  相似文献   

16.
井间地震物理模拟研究   总被引:5,自引:3,他引:2       下载免费PDF全文
井间地震资料具有极高的分辨率,但是,其波场十分复杂,特别是有效一次反射波场能量较弱,且往往被管波等强相干干扰淹没,波场的识别与分离难度较大.地震物理模拟技术是认识复杂地震波场的有效手段.我们通过单一地质体模型、复杂地质体模型和真实井间模拟模型等多种逼真地质模型的物理模拟,揭示了井间地震观测下直达波、透射波、反射波、折射波,多次波和导波等各类波型的特征.物理模拟试验表明,在炮检域平面上分析直达波的信噪比与能量分布,可以使处理人员对旅行时的误差有更直观的认识.所以识别初至与拾取旅行时间时,不仅要在共震源点道集(CSG)上分析拾取,而且还要在共炮检距道集(COG)上分析拾取,这样就提高了直达波识别的可靠性与拾取时间的精度.  相似文献   

17.
Extracting true amplitude versus angle common image gathers is one of the key objectives in seismic processing and imaging. This is achievable to different degrees using different migration techniques (e.g., Kirchhoff, wavefield extrapolation, and reverse time migration techniques) and is a common tool in exploration, but the costs can vary depending on the selected migration algorithm and the desired accuracy. Here, we investigate the possibility of combining the local‐shift imaging condition, specifically the time‐shift extended imaging condition, for angle gathers with a Kirchhoff migration. The aims are not to replace the more accurate full‐wavefield migration but to offer a cheaper alternative where ray‐based methods are applicable and to use Kirchhoff time‐lag common image gathers to help bridge the gap between the traditional offset common image gathers and reverse time migration angle gathers; finally, given the higher level of summation inside the extended imaging migration, we wish to understand the impact on the amplitude versus angle response. The implementation of the time‐shift imaging condition along with the computational cost is discussed, and results of four different datasets are presented. The four example datasets, two synthetic, one land acquisition, and a marine dataset, have been migrated using a Kirchhoff offset method, a Kirchhoff time‐shift method, and, for comparison, a reverse time migration algorithm. The results show that the time‐shift imaging condition at zero time lag is equivalent to the full offset stack as expected. The output gathers are cleaner and more consistent in the time‐lag‐derived angle gathers, but the conversion from time lag to angle can be considered a post‐processing step. The main difference arises in the amplitude versus offset/angle distribution where the responses are different and dramatically so for the land data. The results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers. The same disadvantages with a ray‐based approach will apply using the extended imaging condition relative to a wave equation angle gather solution. Nevertheless, using this approach allows one to explore the relationship between the velocity model and focusing of the reflected energy, to use the Radon transformation to remove noise and multiples, and to generate consistent products from a ray‐based migration and a full‐wave equation migration, which can then be interchanged depending on the process under study.  相似文献   

18.
冯波  王华忠  冯伟 《地球物理学报》2019,62(4):1471-1479
地震波的运动学信息(走时、斜率等)通常用于宏观速度建模.针对走时反演方法,一个基本问题是走时拾取或反射时差的估计.对于成像域反演方法,可以通过成像道集的剩余深度差近似计算反射波时差.在数据域中,反射地震观测数据是有限频带信号,如果不能准确地确定子波的起跳时间,难以精确地确定反射波的到达时间.另一方面,如果缺乏关于模型的先验信息,则很难精确测量自地下同一个反射界面的观测数据同相轴和模拟数据同相轴之间的时差.针对走时定义及时差测量问题,首先从叠前地震数据的稀疏表达出发,利用特征波场分解方法,提取反射子波并估计局部平面波的入射和出射射线参数.进一步,为了实现自动和稳定的走时拾取,用震相的包络极值对应的时间定义反射波的到达时,实现了立体数据中间的自动生成.理论上讲,利用包络极值定义的走时大于真实的反射波走时,除非观测信号具有无限带宽(即delta脉冲).然而,走时反演的目的是估计中-大尺度的背景速度结构,因此走时误差导致的速度误差仍然在可以接受的误差范围内.利用局部化传播算子及特征波聚焦成像条件将特征波数据直接投影到地下虚拟反射点,提出了一种新的反射时差估计方法.既避免了周期跳跃现象以及串层等可能性,又消除了振幅因素对时差测量的影响.最后,在上述工作基础之上,提出了一种基于特征波场分解的新型全自动反射走时反演方法(CWRTI).通过对泛函梯度的线性化近似,并用全变差正则化方法提取梯度的低波数部分,实现了背景速度迭代反演.在理论上,无需长偏移距观测数据或低频信息、对初始模型依赖性低且计算效率高,可以为后续的全波形反演提供可靠的初始速度模型.理论和实际资料的测试结果证明了本文方法的有效性.  相似文献   

19.
角度域共成像点道集是衔接叠前地震数据与储层特征的重要桥梁,对地震偏移成像与储层描述具有重要意义.与克希霍夫偏移和单向波动方程偏移相比,逆时偏移是复杂地区最精确的成像方法.高效稳健地生成逆时偏移角度道集目前仍然是一个挑战.本文主要讨论如何采用光学流方法高效、高质量地提取角度道集.在逆时偏移波场外推过程中,光学流方法可以估计波场传播方向.其中Lucas-Kanade(LK)和Horn-Schunck(HS)方法是光学流方法中两种典型的方法.LK光学流方法是一种局部方法,该方法依赖于局部点的梯度值,但是容易出现奇异现象,HS光学流方法属于全局方法,波场方向估计依赖于整个波场,易受噪声影响,对异常值比较敏感,导致整体波场方向计算精度不高.本文提出采用局部和整体结合(Combining Local and Global,CLG)的光学流方法估计波场传播方向.该方法可以有效地提高波场方向的精度,并且简单高效,便于并行处理.对比HS光学流方法,CLG光学流方法几乎不增加额外的计算量.另外,为了弱化光学流方法无法处理波前重叠问题,本文利用解析波场和方向滤波对波场进行方向分解,仅需波场的空间傅里叶变换即可实现任意波场方向分解,将分解后的波场分别估计波场反向,提取成像结果.进一步地,在估计反射张角和方位角时,本文提出有效的归一化方法和改进的最小二乘除法,提高角度估计的精度和稳定性.最后,理论和实际资料例证了本文提出方法的有效性.  相似文献   

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

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