井间层析成像混合最短射线追踪旅行时反演

传统的曲射线追踪方法存在一定的缺陷,弯曲法得到的旅行时可能收敛到局部最小而不是全局最小。为了使其得到较好的射线路径,本文介绍一种方法对弯曲法进行修改。其起始猜测路径采用由SPR得到的曲折射线,这样既可以获得全局最小走时,又避免了SPR方法的射线路径不够光滑,相当于对SPR所得射线路径做平滑处理。射线追踪路径实例和井间地震层析成像模型反演表明,这种混合方法效果很好。     

逆时偏移技术发展现状与趋势

逆时偏移比单程波波动方程偏移和克希霍夫积分偏移更适于复杂构造成像。基于双程波波动方程进行波场逆时外推,并应用成像条件提取成像值,可以得到逆时偏移数据体。实例分析表明,逆时偏移对盐体侧翼及盐下构造的成像效果要优于单程波波动方程偏移。然而,由于采用有限差分算法以及需要保存大量的震源波场数据,因而逆时偏移的计算成本很高,使其无法适应大数据量的地震数据处理;而由互相关成像条件引入的低频噪声也是困扰逆时偏移实用化的一个问题。在总结和分析逆时偏移技术的进展以及国外地球物理公司关于逆时偏移技术发展策略的基础上,结合研究现状,阐明了逆时偏移未来的发展趋势和研究重点。分析认为,进一步提高逆时偏移的计算效率仍然是未来几年逆时偏移技术发展的重要方向,弹性波逆时偏移和各向异性逆时偏移将会越来越受到重视。     

起伏地表逆时偏移在复杂山前带地震成像中的应用

针对山前带"复杂地表、复杂构造"双复杂地震地质条件造成的资料信噪比低、静校正问题严重、地震波场复杂等一系列地球物理难题,在"真地表"地震成像面的确定及高频静校正的基础上,以基于起伏地表的深度域速度分析与建模为重点,以起伏地表逆时叠前深度偏移为核心建立了一套高精度地震成像处理流程,将长波长静校正问题隐含在偏移成像过程中,直接从起伏地表进行波场延拓和偏移成像,以便更好地应对复杂山前带的地震成像问题。针对性模型试算和实际资料处理表明,该技术在应对双复杂地震资料偏移成像时具有更高的精度,是复杂山前带资料高精度地震成像更理想的技术手段。     

应用反射理论的最小二乘逆时偏移成像

最小二乘偏移是在线性反演理论下求解模型空间的精确解，相对于常规偏移具有更高成像分辨率、保幅性，且可减少成像假象。经典最小二乘偏移以散射波线性化表达为基础，成像实质是估计介质的散射强度。由于地下实际地层以层状介质为主，反射成像是实际应用中最小二乘偏移技术的目标。文中重点讨论散射理论与反射理论两种线性化表达方法的不同，根据实际应用中偏移成像的需求，选择忽略角度信息的反射理论线性化表达反偏移方法，建立了一套估计平均反射系数的最小二乘逆时偏移流程。Sigsbee2a模型验证了最小二乘逆时偏移方法的成像效果；实际三维探区应用结果表明，与常规逆时偏移方法相比，该方法“串珠”成像收敛效果更好，对深层大断裂、层间小断层的刻画更清晰。     

基于APSO-MCMC的叠前三参数同步随机反演方法研究

针对确定性叠前三参数反演存在的精度低、稳定性差、抗噪性差以及过度依赖初始模型的问题,虽然利用随机反演方法实现叠前三参数同步反演有助于解决这些问题,但此类方法反演效率较低。为此,提出了利用自适应粒子群优化-马尔科夫链蒙特卡洛(APSO-MCMC)算法以提高叠前三参数反演的计算效率。结合纵波速度、横波速度以及密度的变化统计扰动关系,建立模型扰动的高斯分布,根据三参数扰动先验分布抽样获取初始粒子群。在每一步迭代过程中,利用跃迁矩阵改变粒子的演化,根据粒子之间的距离计算演化因子,判定粒子群的收敛状态。此外,利用精英策略,帮助粒子实现跳出局部极小值的困境。由于传统叠前近似公式在宽角度区域与非近似公式误差较大,故利用Zoeppritz公式来保证叠前三参数反演在各个角度的精度。NS地区二维海洋实际数据应用表明,该方法估算纵波速度、横波速度和密度有效,稳定收敛,计算效率高于传统的随机方法,具有较好的精度和抗噪性。     

应用GPU的傅里叶有限差分逆时偏移

逆时偏移是针对复杂构造的主流叠前深度偏移成像技术,常规逆时偏移技术大多采用高阶有限差分和成像后滤波算法。由于实际偏移的计算网格较大造成空间频散现象和成像后的去低频噪声滤波会损害有效低频信息,因此成像频宽受到一定影响。为此,提出基于GPU的傅里叶有限差分逆时偏移方法。即采用傅里叶有限差分算法进行波场传播计算,波场外推中直接在时间—波数域做显式波场分解和互相关成像;同时采用基于GPU的算法,显著提高了傅里叶有限差分逆时偏移的计算速度。实际资料应用结果表明,该技术能实现对高、低频有效信息的充分保护和利用,从而提高了地震成像分辨率。