专利名称：Vti介质准p波方程各向异性逆时偏移方法
技术领域：
本发明属于勘探地球物理学中地震资料处理的技术领域，是一种基于各向异性理论的叠前VTI介质准P波方程各向异性逆时偏移方法。
背景技术：
近年来，随着地震勘探精度和要求的不断提高，对复杂地质构造的精确成像受到业界越来越多的关注，复杂地质构造成像的问题主要表现在地层倾角大，横向速度变化剧烈，断层发育、埋藏较深等；地震偏移技术正是为了满足这一需求而不断发展和完善起来的。现有技术中，叠前深度偏移方法主要有Kirchhoff积分法、基于单程波方程的有限差分偏移方法、频率波数域偏移方法及基于双程波方程的逆时偏移方法；其中，逆时偏移方法基于精确的双程波动方程，且允许波向各个方向传播，具有无倾角限制，能适应强横向速度变化及保幅性好等优点；因此，该项技术受到了业界的广泛重视。目前，基于各向同性理论的偏移方法已经比较成熟，在生产实践中也得到了广泛应用，产生了巨大效益。然而，各向异性在地下介质中是广泛存在的，越来越多的研究表明，基于各向同性理论的常规地震偏移方法不能有效解决各向异性介质中地震波的偏移归位问题，成像精度难以保证，因此有必要发展基于各向异性理论的偏移方法。当前，研究最为广泛的各向异性介质是具有垂直对称轴的横向各向同性(Transverse Isotropy with aVertical Axis of Symmetry，简称VTI)介质，一般由周期性薄互层组成。
发明内容
本发明的目的在于克服现有技术中基于各向同性理论常规偏移方法精度低、有倾角限制和不能适应强横向速度变化等缺点，研发了一种基于各向异性理论的叠前逆时深度偏移方法，以解决各向异性介质中地震波的偏移归位问题，获取地下真实构造形态，并为后续处理如AVO分析、属性分析等技术提供振幅和相位信息。实现本发明目的的解决技术方案是一种VTI介质准P波方程各向异性逆时偏移方法，该方法包括如下步骤(1)、采用交错网格对二维VTI介质中的一阶准P波方程及其PML吸收层边界方程进行差分离散，得到两个方程正向延拓和逆时延拓的高阶差分格式；(2)、通过正演模拟计算模型内的炮点正向延拓波�。�(3)、将检波点接收到的共炮点数据中的直达波切除掉，再将每炮数据输入到程序中，通过逆时延拓实现模型内地震波场的逆时传播，得到检波点逆时延拓波�。籊)、将炮点正向延拓波场与检波点逆时延拓波场进行归一化互相关运算，得到模型内各成像点的偏移成像结果；(5)、通过步骤(1)、O)，经计算得到每一炮数据对应的偏移结果，从偏移结果中抽取共成像点道集，进行偏移画弧切除和带通滤波后直接叠加，最后得到偏移剖面。
本发明相对于现有技术，其显著效果是(1)、本发明能够克服基于各向同性理论的常规偏移方法精度低、有倾角限制和不能适应强横向速度变化等缺点，使地震波准确地偏移归位，为解释人员提供更真实的地下构造。O)、本发明能处理横向速度变化强烈和陡倾角地层的偏移成像问题，对复杂模型具有良好的成像能力。(3)、本发明考虑了介质各向异性的影响，在各向异性地区采集的纵波数据应用本发明提供的方法能得到更好的成像效果。
表1为本发明P波波场分量和弹性参数的空间位置表图1是本发明交错网格示意2是本发明完全匹配层吸收边界示意3a是本发明二维VTI Marmousi纵波速度模型3b是本发明二维VTI Marmousi η模型4是本发明震源位于(3000m，4m)处时的合成记录5是本发明二维VTI介质一阶准P波方程逆时偏移剖面6是本发明Kirchhoff叠前时间偏移剖面7是本发明显式有限差分法叠前深度偏移剖面8是本发明叠前逆时深度偏移剖面9是本发明各向同性介质声波方程逆时偏移剖面图
具体实施例方式下面结合附图及实施例对本发明作详细描述一种VTI介质准P波方程各向异性逆时偏移方法，该方法实现过程概括有如下几个步骤1)、采用交错网格对二维VTI介质中的一阶准P波方程及其PML(Perfect MatchedLayer)吸收层边界方程进行差分离散，得到两个方程正向延拓和逆时延拓的高阶差分格式。2)、通过正演模拟计算模型内的炮点正向延拓波场。3)、首先将检波点接收到的共炮点数据中的直达波切除掉，然后将每炮数据输入到程序中，通过逆时延拓实现模型内地震波场的逆时传播，得到检波点逆时延拓波场。4)、将炮点正向延拓波场与检波点逆时延拓波场进行归一化互相关运算，得到模型内各成像点的偏移成像结果。5)、通过所述的步骤1) ,2)计算得到每一炮数据对应的偏移结果，从偏移结果中抽取共成像点道集，进行偏移画弧切除和带通滤波后直接叠加，最后得到偏移剖面。模型算例本发明利用二维各向异性Marmousi模型来检验VTI介质中一阶准P波方程叠前逆时深度偏移算法对复杂模型的成像效果；二维各向异性Marmousi模型如图3所示。
逆时偏移所用的地震记录按以下观测系统通过正演获得，纵波源激发，震源为主频50Hz的Ricker子波，炮点埋深細，800道接收，道间距細，炮间距40m，从模型左端开始放炮，当炮点位于模型左端1600m范围内时，模型左端800道接收，当炮点位于中间6000m范围内时，炮点两侧各400道接收，当炮点位于模型右端1600m范围内时，模型右端800道接收，记录长度3. 15s，共得231炮合成记录，如图4所示为震源位于(3000m，4m)处时的合成记录。应用本发明方法对去除直达波之后的所有炮集进行偏移成像得到的偏移剖面如图5所示。如图6所示是采用Kirchhoff叠前时间偏移得到的偏移结果；如图7所示是采用显式有限差分法叠前深度偏移得到的偏移结果；如图8所示是截取图5右侧3. 0-9. Okm得到的结果；对比图6、图7、图8可看到，Kirchhoff叠前时间偏移对模型浅部地层的成像效果较好，但对模型深部的背斜构造不能正确成像(可参见图6中白色椭圆内区域)，位于背斜顶部的油气储层无法得到正确显示；显式有限差分法叠前深度偏移对深部地层的成像效果有所改善，但对模型的成像精度不高，部分层位�：磺�(如图7中箭头所指处)；采用叠前逆时深度偏移得到的偏移结果中模型的层位成像准确，三条断层、断块、深部高速体、背斜构造、油气储层(如图8中箭头所指处)以及模型的细节特征在偏移剖面上都得到了较好的反映，说明本发明能处理横向速度变化强烈和陡倾角地层的偏移成像问题，对复杂模型具有良好的成像能力。为了分析各向异性对偏移效果的影响，采用各向同性介质中声波方程逆时偏移算法对上述地震记录进行偏移成像，得到的偏移剖面如图9所示。对比图5和图9可看到，采用各向同性偏移算法得到的偏移结果，其效果明显不如采用各向异性偏移算法得到的偏移结果，主要表现在(1)、前者存在较严重的干扰现象(如图9中白色箭头所指处)，且部分层位�：磺澹员洳�(如图9中黑色箭头所指处)，成像精度明显低于后者，在模型左上部各向异性较强的区域这种现象尤为突出(参见图9中白色椭圆内区域)J2)、对模型中部复杂构造的成像效果变差(参见图9中红色椭圆内区域)，陡倾角地层和断层面的连续性变差，精度降低，甚至出现失真现象(如图9中红色箭头所指处)。由上述分析可知，对地震资料的偏移处理不能忽略介质各向异性的影响，在各向异性地区采集的纵波数据用基于各向异性理论的偏移方法能得到更好的成像效果。本发明的方法所依据的原理是
⑴、.二维VTI介质中的一阶准P波方程
St1 dP ρ Sx(la)
dt1 dP ρ dz(lb)
k = -dx(Ic)
δψ dtQk S(Id)
权利要求
1. 一种VTI介质准P波方程各向异性逆时偏移方法，其特征在于该方法包括如下步骤(1)、采用交错网格对二维VTI介质中的一阶准P波方程及其PML吸收层边界方程进行差分离散，得到两个方程正向延拓和逆时延拓的高阶差分格式；O)、通过正演模拟计算模型内的炮点正向延拓波�。�(3)、将检波点接收到的共炮点数据中的直达波切除掉，再将每炮数据输入到程序中，通过逆时延拓实现模型内地震波场的逆时传播，得到检波点逆时延拓波�。籊)、将炮点正向延拓波场与检波点逆时延拓波场进行归一化互相关运算，得到模型内各成像点的偏移成像结果；(5)、通过所述的步骤(1)、(2)，经计算得到每一炮数据对应的偏移结果，从偏移结果中抽取共成像点道集，进行偏移画弧切除和带通滤波后直接叠加，最后得到偏移剖面。
全文摘要
本发明涉及一种VTI介质准P波方程各向异性逆时偏移方法，其步骤包括①采用交错网格对二维VTI介质中的一阶准P波方程及其PML吸收层边界方程进行差分离散，得到两个方程正向延拓和逆时延拓的高阶差分格式；②通过数值计算得到炮点的正向延拓波场和检波点逆时延拓波场，将两者进行归一化互相关运算，得到模型类各成像点的偏移成像结果；③从偏移结果中抽取共成像点道集得到最终的偏移剖面；本发明方法能处理横向速度变化强烈和陡倾角地层的偏移成像问题，也考虑了介质各向异性的影响，在各向异性地区采集的纵波数据用基于各向异性理论的偏移方法能得到更好的成像效果。
文档编号G01V1/28GK102590859SQ20111045943
公开日2012年7月18日 申请日期2011年12月31日 优先权日2011年12月31日
发明者吕彬, 李斐, 王宇超, 胡自多, 韩令贺 申请人:中国石油集团西北地质研究所
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
