转换波波动方程叠前深度偏移

克希霍夫叠前时间偏移被广泛运用于转换波地震成像领域,该方法具有计算效率高、速度模型易于获取的优点。但是,纵横波速度差使纵波和转换波的成像在时间域不匹配,且克希霍夫偏移也无法较好地处理多次到达和保幅等问题。这里提出一种基于波动方程的转换波叠前深度偏移方法,以改善克希霍夫时间偏移的不足。     

共偏移距可分表示法叠前深度偏移

基于双平方根方程的共偏移距可分表示法叠前深度偏移用于复杂介质成像, 该方法在中点-偏移距坐标中同时向下延拓炮点和检波点波场, 实现方式采用正反傅立叶变换.构造的双平方根方程波场延拓算子能够使波数域变量与空间(速度) 域变量分离, 波数域内进行相移计算, 在空间域对因介质横向变速引起的时移作修正.地震数据偏移不需要逐炮计算, 具有较高的效率.对Marmousi模型数据的偏移成像结果显示, 该方法较好地成像强横向变速介质中的复杂构造.      

偏移距域/角度域共成像点道集与偏移速度的关系

为了得到波动方程偏移距域共成像点道集(ODCIGs) 和角度域共成像点道集(ADCIGs) 与偏移速度的关系, 在匀速单层水平反射情况下, 对ODCIGs和ADCIGs随偏移速度的变化进行了定量的推导, 得到偏移速度偏小时ODCIGs与半偏移距成双曲关系, 反之成椭圆关系; 结合深度聚焦分析, 得到偏移速度偏大时, ADCIGs中视入射角大于真入射角, 反之, 小于真入射角; 在速度偏小时, ADCIGs上剩余时差(RMO) 与视入射角的正切值成椭圆关系, 反之成双曲关系.在地震有效入射角度范围内, ADCIGs对于偏大的偏移速度更加敏感.在偏移速度分析(MVA) 中, 从稍大的速度开始分析会对速度分析更加有利.      

两种不同求和路径的叠前时间偏移方法

介绍了等效偏移距偏移和非成像炮集偏移的基本原理,二者都是用一个新的偏移距变量代替常规的偏移距,然后对以新的偏移距变量定义的旅行时双曲线进行动校正和叠加,从而实现叠前时间偏移。叠前时间偏移是对双平方根方程定义的旅行时曲面的绕射求和,而这两种方法的区别就在于其在该曲面上定义的求和路径不同。在金属矿勘探中,金属矿脉由于其成因受构造剧烈运动等因素的影响,往往倾角陡,规模小,大多呈非水平层分布,而这两种方法可以为金属矿成像提供较好的成像方法。     

叠前时间偏移在华北地区的应用

从应用的角度较详细地论述了三维叠前克希霍夫时间偏移的过程。在给定合理的偏移孔径下,三维叠前克希霍夫时间偏移利用均方根速度对叠前地震数据进行逐道处理,从而得到准确的归位数据。叠前时间偏移相对叠后时间偏移而言,前者能够进一步提高地下信息的成像质量及归位准确性,信噪比和分辨率也得到了改善,能够满足精细的地震资料的解释。因此,三维叠前克希霍夫时间偏移对于速度变化平缓、地下构造较复杂的地区具有广泛的应用价值。     

叠前时间偏移成像技术及其应用

零炮检距剖面是地震反射成像过程中的重要成果。通过弯曲反射界面模型,分析了共中心点(CMP)道集与共反射点(CRP)道集的差别,认为当反射界面倾斜时,常规处理中的CMP叠加不能得到准确的零偏移距剖面,而沿CRP道集叠加的叠前时间偏移能得到准确的零偏移距剖面,另外,还介绍了叠前时间偏移的实现方法和步骤,包括预处理、偏移速度场的建立和偏移后CRP道集的处理等,并结合实际资料进行了效果分析。     

压缩三维走时表提高克希霍夫偏移计算效率

地震波走时是克希霍夫叠前深度偏移的重要参数。三维克希霍夫偏移成像由于需频繁读取大数据量的走时表文件,所以计算效率不高。这里提出一种通过压缩三维射线追踪走时表,来提高克希霍夫偏移计算效率的方法：首先规划一个具有规则网格控制点并覆盖所有走时表点集的最小长方体区域,然后以三维三次B样条函数为插值基函数进行最小二乘法曲面体拟合,求出规则网格控制点的数值并以数组形式存储入内存,采用稀疏化存储进一步节省了内存空间。在偏移成像时,再由这些规则网格控制点的数值,使用线性插值公式解编出走时表。实际资料算例验证了该走时表压缩方法不仅近似精度高,计算稳定度高,计算效率高,而且由于省去了频繁进行大数据量走时表文件的读写操作,所以克希霍夫偏移的计算效率提高了二倍以上。     

波动方程角度域共成像道集

针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础.     

等效偏移距在垂直地震剖面中的研究及应用

基于散射理论的等效偏移距方法在地震资料数据处理中具有独特的优势,尤其在速度分析和偏移成像方面。笔者将该理论应用到垂直地震剖面(VSP)数据处理中,研究了基于该理论的VSP数据共散射点道集的形成方法,并试验了VSP共散射点道集在速度分析和成像方面的应用效果。理论模型和实际资料测试结果表明,该方法能方便有效地为VSP提供不受初始速度和地层倾角影响的有效偏移速度,除了能在常规成像方法的成像区域内准确成像外,还能在常规成像方法的盲区外有效成像,在一定程度上扩大了VSP方法的成像范围,具有有较好的研究意义和应用前景。     

真振幅共炮点偏移归位

在克希霍夫偏移中,为了计算用真振幅恢复震源脉冲所必需的复变的加权函数,需要进行两步动学射线追踪计算。一是从激发点成像点,另一是从接点到成像点。     

地震波场叠前零偏移距化算法研究

在能源地震勘探中,为了更好地实现叠前倾斜地层的时差校正,提高速度分析和叠加的精度,提出了一种快速实现叠前部分偏移的算法,即零偏移距化（MZO）算法。该算法应用波场外推的技术及波动方程偏移的有限差分算法,可在共偏移道集上实现,对于3D资料须应用同方位角共偏移距道集。经实际资料测试,其计算速度比同类的DMO方法快一倍,且处理效果相当。      

共散射点道集与角道集串级优化叠前偏移速度分析

推导了基于角度域共成像点道集的叠前深度偏移层析速度分析公式,提出一种共散射点（CSP）道集与角道集串级优化叠前偏移的速度分析方法。该方法通过基于CSP道集的叠前时间偏移速度分析获取初始速度,利用基于角度域共成像点道集（ADCIGs）的叠前深度偏移速度分析进行速度更新。实现步骤概括为：将叠前地震数据映射为CSP道集,利用CSP道集叠加速度谱拾取能量团获取均方根（RMS）速度场;通过Dix公式将RMS速度转换为层速度作为层析的初始输入速度,基于ADCIGs实现叠前深度偏移层析速度反演,最终得到高精度的叠前偏移速度场。断层模型和实际资料试算结果验证了该方法的正确性和有效性。     

基于F-K域和Curvelet-中值滤波联合去噪的混采数据分离方法

由于炮分离过程会不可避免地损失部分有效信息,为解决这一问题,笔者建立模型试验,将混采数据转至频率域进行F-K滤波预处理,在保护地震同相轴的同时压制大部分混合噪声;并变换回时间域共偏移距道集进行小时窗中值滤波,同时设计迭代公式在共检波点道集进行Curvelet阈值迭代去噪;通过不断优化得到最后分离结果。模拟数据和实际数据处理结果证明:经本文方法分离后的单震源记录清晰,具有较好的实际应用价值。     

共聚焦点道集偏移速度建模

偏移速度建模可在多种道集实现.基于等时原理和差异时移分析,提出了共聚焦点道集偏移速度建模方法.偏移速度的更新是利用约束参数迭代反演实现,模型的参数化主要依据实际情况而定.为适应横向变速,选用了垂直速度梯度参数.模型试算表明:偏移速度相对误差在0.25%范围内,反射层深度误差小于10m.     

共反射面叠加在实际地震资料处理中的应用

共反射面叠加方法运用3个地震波场参数来描述地下反射面元时距关系,不对地下反射界面形状做任何假设,并且是一种不依赖于宏观速度模型的反射波成像方法。共反射面叠加利用了邻近多个共反射点道集的相似性,提高了资料的覆盖次数,从而压制了随机噪声,增强反射波同相轴的连续性,大幅提高了地震资料的信噪比。通过在实际地震资料处理中的应用,证明了该方法的这些特性。     

The Offset-Domain Prestack Depth Migration with Optimal Separable Approximation

The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.