基于改进的L_1范数非因果滤波器多次波自适应减方法

多次波压制是地震勘探中的重要课题。基于匹配滤波器设计的多次波自适应减方法,对L1范数多次波自适应减方法进行改进,增加了一个门限值参数。同时针对常规自适应减方法中采用的因果滤波器只能沿时间轴正向进行多次波匹配的不足,采用了非因果滤波器进行多次波自适应减,增加了算子零时间的位置参数。通过引入新参数增强方法的适应性,改进方法提高了计算速度。通过对所设计的1D数据模型和合成地震记录的处理表明,在预测的多次波与原始记录中的多次波存在振幅和到达时间差别的情况下,非因果滤波器能给出较理想的结果;同时改进的自适应减方法可以减少计算时间。     

面向海上地震资料的τ-p域保幅处理方法

AVO分析和反演需要地震资料处理时最大限度地保持振幅相对关系,然而常规处理中球面扩散补偿、拉伸畸变和多次波均会影响振幅相对关系,海上地震资料尤为严重。提出一套面向海上地震资料的τ-p域保幅处理方法。该方法通过τ-p域无需球面扩散补偿且多次波周期性增强的特点,减少球面扩散补偿与多次波对振幅相对关系的破坏,通过τ-p域拉伸因子非时空变换的特点减少动校拉伸产生的频率和能量畸变。最后通过挪威海上实际资料的处理验证了方法的合理性与可行性。     

等时叠加波动方程叠前正演

非零炮检距地震正演模拟在复杂区叠前地震资料处理和解释中有着重要意义。通常使用射线追踪法和全程波动方程法来实现这一模拟。通过提出等时叠加原理和数学检波器表达式，形成了非零炮检距地震正演的新一类方法——单程波动方程法。该方法与射线追踪法相比能保持波的动力学特征，正演结果与实际地震记录相似性好。并且在计算过程中，人工干预的程度低，操作简便。与全程波动方程法相比，该方法计算效率高，且不产生直达波、多次波等干扰，信噪比高，与目的层相对应的一次波突出，波型识别难度不大。因此，单程波动方程法特别适合于复杂油气藏反射地震资料的正、反演研究。     

层状弹性介质中地震波运动学问题研究

完全弹性介质中，进行地震波射线追踪时需要确定射线的出射角；在非弹性介质中，除了确定此射线出射角之外，还需要确定另外一个角，即波前面与常振幅平面的夹角──衰减角。基于惠更斯原理，研究了确定射线出射角与衰减角的一种有效方法，进而对层状非弹性介质中地震波运动学问题进行了一定的理论研究与数值计算。     

层状弹性介质中地震波运动学问题研究

完全弹性介质中，进行地震波射线追踪时需要确定射线的出射角；在非弹性介质中，除了确定此射线出射角之外，还需要确定另外一个角，即波前面与常振幅平面的夹角－衰减角。基于惠更斯原理，研究了确定射线出射角与衰减角的一种有效方法，进而对层状非弹性介质中地震波运动学问题进行了一定的理论研究与数值计算。     

波场延拓法VSP正演模拟

波场延拓法模拟VSP,可以根据具体情况的需要同时对上、下行波场进行模拟,也可以只模拟上行波场或下行波场。对多次波能够有效地控制,通常在只关注一次波的情况下,可以不对多次波进行模拟。该方法相对其他波动方程正演模拟方法,有更高的计算效率。采用分步Fourier法延拓波场,既保证了延拓过程的稳定性,又不失较为复杂模型的模拟精度。模型试验表明,该方法能够有效地对VSP波场进行模拟,模拟的结果正确、可靠。     

鄂南地区地震储层预测存在问题与思考

通过分析鄂尔多斯盆地南部地区复杂黄土塬地表和延安组煤层对长8段储层描述的影响,认为受煤层多次波影响,地震振幅属性在反映该区储层发育情况时有很强的多解性,采用了基于地震波形特征的地震反射结构分析技术——透镜体描述河道砂体的方法。重点分析了延安组多套煤层形成的层间多次波对地震波形的改造问题,通过正演模拟技术分析了层间多次波的产生机理和发育特征,从资料解释方面提出了延安组煤层间多次波发育区的识别方法,完善了基于地震波形特征的砂体描述技术,提高了该区储层预测的精度。     

利用多次波和二维匹配滤波的近偏移距数据重构方法

在地震资料采集时,受地表条件、震源等因素的影响,几乎无法采集到近偏移距数据,不利于后续数据的处理。针对这一问题,研究了一种利用多次波和二维匹配滤波的近偏移距数据重构方法。首先给出了由表层多次波构建准一次波的互相关方法;然后在时间-空间滑动窗口内利用准一次波计算二维匹配滤波器,并用该滤波器重构近偏移距数据;最后使用均方根振幅校正方法校正重构数据的振幅。整个数据重构的过程在频率域进行,避免了基于维纳滤波方法求取滤波器的大量计算;同时使用快速傅里叶变换算法,提高了计算效率。与一维匹配滤波的近偏移距重构方法相比,采用二维匹配滤波方法的重构数据在空间上连续性更好。通过对Sigsbee 2B数据的处理及实际资料的测试,验证了该方法的正确性和有效性。     

Love型槽波透射记录的频域合成算法

利用射线追踪对网格化的二维模型进行走时计算,在频率域对Love型透射槽波记录进行了合成。在考虑了频散、吸收衰减等对波场振幅、相位影响的基础上,对含薄煤带模型的透射槽波理论记录进行了频域合成,并进行了波场特征分析。结果表明:频域合成的槽波透射记录在时间、频散特征、能量等属性上与理论分析一致,且与三维弹性波波动方程模拟的基阶Love型槽波波场特征相似,但频域合成算法的计算速度快1~2个数量级。     

反射率法波场正演模拟及其在惠州某区薄储层地震响应模拟中的应用

惠州某区薄储层发育,波场特征复杂,前人研究表明,常规正演模拟方法难以解决薄层模拟精度问题。鉴于此,采用了基于薄层的反射率法地震模拟,推导了相应的计算公式,并编制了计算程序,该方法模拟的地震记录中包含反射波、转换波、多次波信息,符合薄储层地震波传播规律。通过对理论模型正演模拟研究表明,地震波场中包含多次波、转换波信息,且得到了走时验证,证明该方法是可行的;通过应用该方法对惠州某区实际资料的正演模拟,揭示了薄砂体的展布及连通关系,得到了HZ-X井的钻井验证,取得了较好的应用效果。     

改进的地震模型初值射线追踪方法

初值射线追踪方法是现代射线追踪方法中的一个很重要的理论,它克服了两点法射线追踪方法耗时的计算效率问题。以程函方程为基础,对初值射线追踪方法进行改进,即利用平方慢度来替换模型中的速度参数,使得程函方程产生解析解,从而进一步导出当射线遇到界面时的反射和透射慢度向量的计算表达式,以及反射、透射系数的函数表达式。通过对简单的两层界面向斜模型及复杂的多层盐丘模型的射线追踪,表明该初值射线追踪方法的改进相比于以往的龙格库塔离散数值解法,不但使射线追踪效率得到了大幅度提高（10倍左右）,且也扩大了射线法使用范围。     

强反射能量屏蔽补偿在渤海C构造中的应用

在渤海海域西南部广泛分布着新生界的溢流相火成岩,高速层的能量屏蔽是造成火成岩发育区下覆地层地震波场信号弱的主要原因之一。传统的自动增益控制方法(AGC)在分时窗滑动补偿下覆地层被屏蔽的弱信号的同时,破坏了火成岩与下覆地层的能量关系,且时窗内火成岩层之间的弱信号也无法得到有效恢复。这里从反射系数与振幅补偿的关系出发,基于模型论证了利用射线追踪提取地震波传播路径上的反射系数来求取补偿因子,进而开展相对振幅补偿过程的有效性。在渤海油田多个构造的振幅补偿结果,较好地恢复了火成岩屏蔽地层的能量损失,且地层间的振幅关系也得到相对保持,取得了较好的实际应用效果。     

柴达木盆地自然伽玛曲线在古气候及沉积环境研究中的应用

应用自然伽玛测井（GR）曲线反演古气候和古环境的方法在国际上逐渐得到重视并取得了进展，柴达木盆地第四纪应用其反演2.85Ma B.P来的古气候变化也取得了成功。本文以柴达木盆地侏罗系冷科1井和雨33井自然伽玛曲线并结合孢粉等地质资料进行了古气候与沉积环境分析，证明GR曲线是有价值的古气候研究参数。GR值的大小及古气候的干湿变化和沉积水体的深浅变化均有内在联系，GR曲线振幅与频率的变化快慢反映出沉积环境动荡与稳定的相对特征，据此对两口井的沉积环境演变作了具体分析。结合古气候与沉积环境综合分析，GR曲线的异常波动还具有特殊的指相意义。GR值增大-古气候加湿-沉积水体加深的同步表现特征是一般性规律，构造-古地理背景的变动可改变这一对应关系，因而提供了研究古地理演变的信息，进一步丰富了GR曲线研究的地质涵义。说明在侏罗纪古气候和古环境研究上，GR曲线是一有价值的参考指标。     

多频去耦在双频激电数据去耦中的应用

在频率域激电法所用频段0.1~10Hz范围内,单纯由激电效应决定的视电阻率幅值的频率特性,可近似地表示为频率对数的线性函数;电磁耦合效应使视电阻率幅值发生的变化,可近似表示为频率的幂函数。因此,多频测量法被提出用于频率域激电法的电磁耦合校正。SQ-3型双频激电仪设有4个频组,频率从1/13~8Hz,正处于近似公式适用的范围。将前人多频去耦的思路应用于双频激电数据的处理中,根据双频激电的特点,提出双频仪的多频组电磁耦合校正,经野外实际应用,去耦效果良好。     

吸收介质中射线路径,走时与振幅计算——基于一阶扰动方法

实际地层由于波前扩散、地层吸收和散射等原因发生衰减,为得到更加详细的地下信息和分辨率更高的图像,基于波前构建法的射线追踪,实现了吸收介质中对射线路径,走时与振幅的计算。计算表明吸收介质中的波动属性是复值;当Q为5~20时,衰减是可见的,明显的;而当Q30时,衰减微弱,为弱衰减介质。对吸收介质中的射线追踪进行计算,可提高地震波的分辨率,得到更清晰的地震记录。     

粘弹介质中高斯束二维理论地震记录的合成

本文对无限空间二维变化粘弹介质中的高斯束及有关过界面的情况进行了讨论。在讨论中发现,利用高斯束法既可不用定点射线追踪,还可保证处处规则,较近年来发展较快的标准射线方法有很大的灵活性,适于处理各种复杂地质情况。并对理论地震记录的合成方法作了简要介绍。     

浅海海底表层速度建模技术及其应用研究

浅海海底起伏和速度变化对OBC资料成像质量产生较大影响。建立浅海海底表层速度模型不仅能够解决OBC资料的静校正问题,也可用于海底反射系数计算、双检资料合并、多次波压制等,但是,目前针对浅海海底表层速度建模的研究还不多。文中提出了针对浅海地区OBC资料的海底表层速度建模的三维初至走时反演技术,主要包括:(1)震源位置校正技术。根据地震波在海水中传播特征,把在海水中激发的震源位置校正至海底,使震源和接收点都位于海底,利于初至走时反演计算;(2)快速三维初至走时反演方法。利用回折波走时和射线方程,形成了高效率初至走时反演方法。将该技术应用于胜利油田莱州湾浅海区海底OBC资料的处理中,建立了三维海底表层速度模型,用此速度模型进行静校正,取得了良好的应用效果。     

利用线性旅行时插值射线追踪 计算近地表模型初至波走时

利用线性旅行时插值射线追踪对近地表模型进行正演计算,可以快捷、准确地获得初至波走时和射线路径。由于该算法计算的初至波不局限于折射波,因此很好地解决了浅层折射勘探中的低速"隐蔽层"问题;而且,由于该算法是基于网格划分和线性插值,因此它不仅可以追踪任意复杂介质的初至波,而且可以使得追踪的初至波射线路径逼于真实,避免了同类算法直接连接网络节点形成射线路径的缺陷(路径过于弯折,计算走时偏大)。将LTI算法同其他几种算法的追踪结果进行的对比和分析表明,LTI算法在计算初至波走时和射线路径方面较其他算法更为精准、稳定,是一种有效的射线追踪方法。     

Prediction of peak particle velocity using multi regression analysis: case studies

Ground vibrations produced from blasting operations cause structural vibrations, which may weaken structure if it occurs at the resonant frequency. Measurable parameters associated with ground vibrations are peak particle velocity (PPV), amplitude and dominant frequency (frequency of highest PPV amongst translational, vertical and horizontal vibrations). In this paper, an attempt is made to correlate measurable parameters associated with ground vibrations with scaled distance. Using the correlated data, it was found that a predictor equation can be determined for the amplitude and PPV, but not for dominant frequency as it is dynamic and depends upon infinitesimal changes that occur within a number of other parameters. Another analysis of the same is made using multiple linear regression analysis. This included predicting the PPV using scaled distance, maximum charge per delay, amplitude as predictors. A considerable improvement is seen in the prediction on adding the interaction of the predictors in multiple regressions. A comparison of different combination of predictors is made so as to assess the best combination giving the best R² value for the given mine. Frequency is also plotted using the aforementioned method. However, it was found that the dominant frequency cannot be predicted with high accuracy even with this method.