波动方程模拟PML吸收影响因素分析

地震数值模拟中,完全匹配层(PML)边界能有效地吸收衰减地震波,得到无边界反射干扰的波场快照和地震记录。在前人研究的基础上,进一步分析了分裂完全匹配层(SPML)的衰减机制,通过波动方程模拟,分析了震源主频、空间网格间距、介质速度等参数对SPML边界吸收衰减特征的影响。得到了不同条件下,PML边界对地震波的吸收效果。通过分析对比,得出了震源主频对PML吸收效果无直接影响,空间网格间距与PML吸收效果成反比,高速层PML吸收效果缓慢等结论。并通过Marmousi2模型测算,对结论进行了验证。最后给出了复杂模型、不同尺度数值模拟中SPML参数的选取方法,为地震数值模拟中SPML边界参数的定量选取奠定了基础。     

共炮检距道集波动方程保幅叠前深度偏移方法

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

多尺度快速匹配追踪多域联合地震反演方法

多尺度快速匹配追踪多域联合地震反演是一种通过地震数据多尺度分解的迭代反演方法.与此同时,在快速匹配追踪算法中引入低频模型约束,有效提高了收敛精度,使反演结果具有丰富的高低频信息.首先通过对大尺度地震资料进行反演得到低频背景.在此基础上,采用中尺度与小尺度地震数据进行逐级迭代用以获得高频数据,因而有效缓解了常规反演方法对于初始模型精度的依赖.最后利用理论模型与实际地震数据进行测试,通过与常规时间域反演方法的反演结果进行对比可以看出,本文方法在地层连续变化处依然可以对变化地层进行精确刻画,且在纵向分辨率提升的同时保持了较好的横向连续性.     

TTI介质刚度矩阵的分解

刚度矩阵是表征传播介质特性的重要参数.TTI介质的刚度矩阵由VTI/HTI介质的刚度矩阵转换而得,因而TTI介质与VTI/HTI介质之间存在着某种联系.本文将TTI介质的刚度矩阵表示成VTI/HTI介质的刚度矩阵与一个‘扰动'矩阵之和的形式,从另一个角度揭示了TTI介质与VTI/HTI介质的关系.本文在TTI介质刚度矩阵分解的基础上,将TTI介质的波动方程分解成VTI/HTI介质波动方程和‘扰动介质'波动方程,然后对VTI/HTI介质波动方程和‘扰动'介质波动方程分别进行了模拟.将两个模拟结果进行合成并与TTI介质的波动方程模拟结果进行了对比,验证了方法的准确性和有效性.理论分析与数值模拟结果表明,TTI介质可认为是VTI/HTI介质的一种延伸.     

天然气水合物地层岩石物理模型构建

天然气水合物地震反射特征的识别和物性信息的定量提取存在多解性与不确定性,岩石物理模型是将水合物地层地球物理观测值地震信息定量转化为物性信息的桥梁,而传统的岩石物理模型少有对微观孔隙结构和孔隙充填水合物剪切性质的描述.文章针对孔隙充填和颗粒支撑两种微观分布模式的水合物地层,重点考虑水合物地层的微观孔隙结构以及水合物的剪切性质,基于等效介质理论定量描述地层矿物组分特征和孔隙连通性及形状,利用斑块饱和理论和广义Gassmann理论定量描述孔隙充填水合物的剪切性质,在此基础上构建了两种模式水合物地层的岩石物理模型,揭示了水合物地层宏观弹性性质与微观物理性质之间的定量关系.数值研究发现,纵横波速度比随孔隙度和水合物饱和度增大而减小;颗粒支撑模式地层的纵横波速度对水合物含量较为敏感,且孔隙越狭长,敏感性越高;孔隙充填模式地层的纵横波速度比在高水合物饱和度时对水合物组分剪切性质的敏感性更高.实验数据和神狐海域的实际资料应用结果表明,岩石物理模型可有效地计算水合物地层宏观弹性性质与微观物性特征之间的定量关系,提供常规测井中缺乏的横波速度信息,确定对水合物含量指示能力较强的弹性参数,可以根据实际数据求取水合物饱和度、孔隙纵横比等物性参数,为天然气水合物地层的定量解释和资源评价提供理论依据和数据支持.     

基于逆算子估计的AVO反演方法研究

传统反演算法以优化算法为主,而基于逆算子估计的AVO反演算法则利用了直接求逆的思路.算法的关键在于寻找存在逆函数的子域,进而可以在子域内直接求逆,这种解决反问题的思路不同于一般的优化类算法所采用的直接搜索解的方式,具有更高的效率.AVO反演利用了振幅随着偏移距的变化特征,反演的精度受到地震资料质量的影响,通过加入L1范数约束以及合理的初始模型有助于提高反演的稳定性以及准确度.模型测算和实际应用表明,基于逆算子估计的AVO反演方法具有较高的精确程度和可靠性.     

基于FFT-MA谱模拟的快速随机反演方法研究

虽然基于地质统计学的随机反演方法能够有效融合测井资料中的高频信息,但计算效率低,占用内存大,限制了它在实际资料中的应用.本文在保留传统随机反演方法优点的基础上,创造性地引入傅里叶滑动平均(Fast Fourier Transform-Moving Average,FFT-MA)谱模拟进行频率域的地质统计模拟,并利用逐步变形算法(Gradual Deformation Method,GDM)确保模拟结果与实际地震数据的匹配,构建了基于FFT-MA谱模拟的新的快速随机反演方法.与常规随机反演相比,新方法不仅分辨率高,而且能够使反演解得到快速收敛,有效提高计算效率,减少内存占用.模型试算获得了与理论模型吻合度较好的高分辨率反演结果.实际资料分析也表明新方法所得到的高分辨率反演结果能够对薄互储层进行良好的展示,为薄储层的识别提供高效可靠的技术支持.     

一种基于全局优化的交错网格有限差分法

在地震波场数值模拟中, 交错网格有限差分技术得到了广泛的应用, 但是在弹性模量变化较大时, 通常会因插值而导致模拟误差增大. 旋转交错网格可以很好地克服这个缺点, 因而适合于各向异性介质正演模拟. 但是对于同样大小的网格单元, 旋转交错网格需要的步长比常规交错网格要大, 这会使梯度和散度算子的误差增大因而更易产生空间数值频散. 针对这些问题, 本文提出了旋转交错网格与紧致有限差分相结合的方法, 并基于模拟退火算法进行全局优化, 压制数值频散, 拓宽波数范围. 数值模拟结果表明, 此方法可以有效地压制数值频散, 且具有较高的模拟精度.      

基于方位各向异性弹性阻抗的裂缝岩石物理参数反演方法研究

裂缝储层岩石物理参数的准确获得对地下裂缝预测具有重要意义,而叠前地震反演是获得裂缝岩石物理参数的有效手段.本文从裂缝岩石物理等效模型的构建出发,从测井数据上估测了裂缝岩石物理参数,通过推导含裂缝岩石物理参数的方位各向异性弹性阻抗公式,探讨了基于方位各向异性弹性阻抗的裂缝岩石物理参数地震反演方法.实际工区地震数据应用表明,基于方位各向异性弹性阻抗的裂缝岩石物理参数反演方法合理、可靠,可以降低裂缝岩石物理参数估测的不确定性,为地下裂缝预测提供有力的依据.