Kirchhoff叠前时间偏移角度道集

第５３卷第５期　２０１０年５月　地球物理　学　报　Ｖｏｌ＿５３，Ｎｏ．５　Ｍａｙ，２０１０　ＣＨＩＮＥＳＥ　ＪＯＵＲＮＡＬ　ＯＦ　ＧＥＯＰＨＹＳＩＣＳ　邹振，刘洪，刘红伟．Ｋｉｒｃｈｈｏｆｆ叠前时间偏移角度道集．地球物理学报，２０１０，５３（５）：１２０７～１２１４，ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．　０００１—５７３３．２Ｏ１０．０５．０２３　Ｚｏｕ　Ｚ，Ｌｉｕ　Ｈ，Ｌｉｕ　Ｈ　Ｗ．Ｃｏｍｍｏｎ—ａｎｇｌｅ　ｇａｔｈｅｒｓ　ｂａｓｅｄ　ｏｎ　Ｋｉｒｃｈｈｏｆｆ　ｐｒｅ－ｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ．Ｃｈｉｎｅｓｅ　Ｊ．Ｇｅｏｐｈｙｓ．（ｉｎ　Ｃｈｉｎｅｓｅ），２０１０，５３（５）：１２０７￣１２１４，ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．０００１—５７３３．２０１０．０５．０２３　Ｋｉｒｃｈｈｏｆｆ叠前时间偏移角度道集　邹　振　，刘　洪　，刘红伟　１中国科学院地质与地球物理研究所中国科学院油气资源研究重点实验室，北京１０００２９　２中国科学院研究生院，北京１０００３９　摘要提出三维Ｋｉｒｃｈｈｏｆｆ叠前时间偏移角度域共像点道集的改进算法，克服传统角度求取算法局限，可计算相　对倾斜地层法线入射角；与Ｋｉｒｃｈｈｏｆｆ直射线叠前时间偏移求角度算法相比，本文方法考虑射线弯曲效应，包含层　速度，角度范围加大，更接近真实入射角；计算走时采取弯曲射线或者适应线性横向变速介质的非对称走时等算　法，角度道集在大角度处得到拉平；采用相对保幅的权因子以及覆盖次数校正技术，有利于叠前ＡＶＡ反演．模型测　试结果表明：叠前时间偏移角度道集，相对ＣＭＰ、ＣＲＰ所转化角度道集，更准确反应ＡＶＡ效应；实际三维数据测试　表明本文方法可以提供品质优良的角度道集，适用于ＡＶＡ分析、反演，提高叠前反演分辨率．　关键词Ｋｉｒｅｈｈｏｆｆ叠前时间偏移，角度道集，角度域共像点道集　１３０Ｉ：１０．３９６９／ｊ．ｉｓｓｒ￣０００１—５７３３．２０１０．０５．０２３　中图分类号Ｐ６３１　收稿日期２００９—０８—３１，２０１０－０２—０３收修定稿　Ｃｏｍｍｏｎ—ａｎｇｌｅ　ｇａｔｈｅｒｓ　ｂａｓｅｄ　ｏｎ　Ｋｉｒｃｈｈｏｆｆ　ｐｒｅ－ｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ　ＺＯＵ　Ｚｈｅｎ　，ＬＩＵ　Ｈｏｎｇ　，ＬＩＵ　Ｈｏｎｇ—Ｗｅｉ　’　１　Ｋｅｙ　Ｌａｂｏｒａｔｏｒｙ　ｏｆ　Ｐｅｔｒｏｌｅｕｍ　Ｒｅｓｏｕｒｃｅｓ　Ｒｅｓｅａｒｃｈ，Ｉｎｓｔｉｔｕｔｅ　ｏｆ　Ｇｅｏｌｏｇｙ　ａｎｄ　Ｇｅｏｐｈｙｓｉｃｓ，Ｃｈｉｎｅｓｅ　Ａｃａｄｅｍｙ　ｏｆ　Ｓｃｉｅｎｃｅｓ　Ｂｅｉｊｉｎｇ　１０００２９，Ｃｈｉｎａ　２　Ｇｒａｄｕａｔｅ　Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　Ｃｈｉｎｅｓｅ　Ａｃａｄｅｍｙ　ｏｆ　Ｓｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ　１０００３９，Ｃｈｉｎａ　Ａｂｓｔｒａｃｔ　Ｉｎ　ｔｈｉｓ　ｐａｐｅｒ　ｗｅ　ｐｒｅｓｅｎｔ　ｔｈｅ　ｉｍｐｒｏｖｅｄ　ａｌｇｏｒｉｔｈｍ　ｆｏｒ　ｃｏｍｍｏｎ—ａｎｇｌｅ　ｇａｔｈｅｒｓ　ｂａｓｅｄ　ｏｎ　Ｋｉｒｃｈｈｏｆｆ　ｐｒｅ—ｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ．Ｔｈｅ　ａｌｇｏｒｉｔｈｍ　ｃａｎ　ｃａｌｃｕｌａｔｅ　ｒｅｌａｔｉｖｅｌｙ　ｔｒｕｅ　ｉｎｃｉｄｅｎｔ　ａｎｇｌｅ　ｕｎｄｅｒ　ｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａ１　ｃｉｒｃｕｍｓｔａｎｃｅ。ｗｈｉｌｅ　ｔｈｅ　ｔｒａｄｉｔｉｏｎａｌ　ｒａｙ—ｐａｒａｍｅｔｅｒ　ｍｅｔｈｏｄ　ｃｏｕｌｄ　ｎｏｔ　ｈａｎｄｌｅ　ｔｈｅ　ｉｎｃｉｄｅｎｔ　ａｎｇｌｅ　ｆｏｒ　ｔｈｅ　ｄｉｐｐｉｎｇ　ｅｖｅｎｔｓ．Ｃｏｍｐａｒｅｄ　ｗｉｔｈ　ｔｈｅ　ｓｔｒａｉｇｈｔ—ｒａｙ　Ｋｉｒｃｈｈｏｆｆ　ａｎｇｌｅ　ｄｏｍａｉｎ　ｇａｔｈｅｒｓ　ｍｅｔｈｏｄ　ｐｒｏｐｏｓｅｄ　ｉｎ　ｒｅｃｅｎｔ　ｙｅａｒｓ，ｏｕｒ　ｍｅｔｈｏｄ　ｉｓ　ａｂｌｅ　ｔｏ　ｅｘｔｒａｃｔ　ｍｕｃｈ　ｗｉｄｅｒ　ａｎｇｌｅ　ｆｏｒ　ｔａｋｉｎｇ　ｉｎｔｅｒｖａｌ　ｖｅｌｏｃｉｔｙ　ａｎｄ　ｒａｙ　ｂｅｎｄｉｎｇ　ｅｆｆｅｃｔ　ｉｎｔｏ　ａｃｃｏｕｎｔ．Ｍｏｒｅｏｖｅｒ，ｗｅ　ａｐｐｌｙ　ｎｏｎ—　ｈｙｐｅｒｂｏｌｉｃ　ｍｏｖｅｏｕｔ　ｅｑｕａｔｉｏｎ　ｔＯ　ｃａｌｃｕｌａｔｅ　ｔｈｅ　ｔｒａｖｅｌ　ｔｉｍｅ，ｔｈｕｓ　ｔｈｅ　ｃｏｍｍｏｎ—ａｎｇｌｅ　ｇａｔｈｅｒｓ　ａｒｅ　ｆｌａｔｔｅｒ　ａｔ　ｔｈｅ　ｌａｒｇｅ　ａｎｇｌｅ　ｔｈａｎ　ｓｔｒａｉｇｈｔ　ｒａｙ　ｍｅｔｈｏｄ．Ｅｆｆｉｃｉｅｎｔ　ｔｒｕｅ—ａｍｐｌｉｔｕｄｅ　ｗｅｉｇｈｔｉｎｇ　ｆｕｎｃｔｉｏｎ　ａｓ　ｗｅｌｌ　ａｓ　ｈｉｔ—ｃｏｕｎｔ　ｃｏｒｒｅｃｔｉｏｎ　ｉｓ　ｂｅｎｅｆｉｃｉａｌ　ｆｏｒ　ＡＶＡ　ｉｎｖｅｒｓｉｏｎ．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｏｆ　ｓｙｎｔｈｅｔｉｃ　ｍｏｄｅｌ　ｄａｔａ　ｈａｖｅ　ｖｅｒｉｆｉｅｄ　ｔｈａｔ　ｃｏｍｐａｒｅｄ　ｗｉｔｈ　ａｎｇｌｅ　ｇａｔｈｅｒｓ　ｔｒａｎｓｆｏｒｍｅｄ　ｆｒｏｍ　ｏｆｆｓｅｔ　ｔｏ　ａｎｇｌｅ　ｄｏｍａｉｎ　ｗｉｔｈ　ＣＭＰ　ｇａｔｈｅｒｓ　ａｎｄ　ＣＲＰ　ｇａｔｈｅｒｓ　ａｆｔｅｒ　ｍｉｇｒａｔｉｏｎ，ｔｈｅ　ｃｏｍｍｏｎ—ａｎｇｌｅ　ｇａｔｈｅｒｓ　ｂａｓｅｄ　ｏｎ　Ｋｉｒｃｈｈｏｆｆ　ｐｒｅ—ｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ　ｈａｖｅ　ｔｈｅ　ｆｏｌｌｏｗｉｎｇ　ａｄｖａｎｔａｇｅｓ：ｍｕｃｈ　ｍｏｒｅ　ａｃｃｕｒａｔｅ　ｉｎｃｉｄｅｎｔ　ａｎｇｌｅ，ｔｈｅ　ｅｆｆｅｃｔ　ｏｆ　ａｍｐｌｉｔｕｄｅ　ｖｅｒｓｕｓ　ａｎｇｌｅ（ＡＶＡ）ｉＳ　ｒｏｂｕｓｔ　ａｎｄ　ｃｌｅａｒ．Ｔｈｅ　ｒｅａｌ　３－Ｄ　ｄａｔａ　ｔｅｓｔ　ｓｈｏｗｅｄ　ｔｈａｔ　基金项目　国家“８６３”高技术研究发展计划项目（２００６ＡＡ０９Ａ１０２—０８），国家“９７３”重点基础研究发展计划项目（２００７ＣＢ２０９６０３）联合资助　作者简介邹振，男，１９８２年生，中国科学院地质与地球物理研究所在读博士生，主要从事地震成像及反演方法的研究．　Ｅ—ｍａｉｌ：ｚｏｕｚｈｅｎ＠ｍａｉｌ．ｉｇｇｃａｓ．ａｃ．ｃｎ　地球物理学报（Ｃｈｉｎｅｓｅ　Ｊ．Ｇｅｏｐｈｙｓ．）　５３卷　ｔｈｅ　ａｎｇｌｅ　ｇａｔｈｅｒｓ　ｐｒｏｖｉｄｅｄ　ｗｉｔｈ　ｏｕｒ　ｍｅｔｈｏｄ　ｗｅｒｅ　ｂｅｎｅｆｉｃｉａｌ　ｆｏｒ　ＡＶＡ　ａｎａｌｙｓｉｓ　ａｎｄ　ＡＶＡ　ｉｎｖｅｒｓｉｏｎ，ａｎｄ　ｃｏｕｌｄ　ｉｍｐｒｏｖｅ　ｐｒｅ—ｓｔａｃｋ　ｉｎｖｅｒｓｉｏｎ　ｒｅｓｏｌｕｔｉｏｎ．　Ｋｅｙｗｏｒｄｓ　Ｋｉｒｃｈｈｏｆｆ　ＰＳＴＭ，Ｃｏｍｍｏｎ—ａｎｇｌｅ　ｇａｔｈｅｒｓ（ＣＡＧ），　Ａｎｇｌｅ—ｄｏｍａｉｎ　ｃｏｍｍｏｎ　ｉｍａｇｅ　ｇａｔｈｅｒｓ（ＡＤＣＩＧ）　基于平面波分解叠前时间偏移所产生的射线参数Ｐ　１　引　言　振幅随偏移距变化（ＡＶＯ）技术广泛应用于储　道集，并用于ＡＶＯ分析，但推广到三维尚存在难度．　本文首先分析传统算法偏移距域转换角度域所　存在问题；提出新的三维Ｋｉｒｃｈｈｏｆｆ叠前时间偏移　层反演、预测岩石物性，提高了储层预测准确性．　ＡＶＯ技术理论基础是振幅随入射角变化（Ａ、，Ａ），二　者仅在地下反射层为水平界面情况下等价¨１］．偏移　距向入射角／反射角转换是ＡＶＯ关键步骤．国内工　业界俗称的角度道集（ＣＡＧ）是将ＮＭＯ校正后共　中心点道集（ＣＭＰ）或叠前偏移后偏移距域共像点　道集（ＣＩＧ）／共反射点道集（ＣＲＰ）映射到角度域．偏　移距域到角度域映射较少采用对层速度要求较高的　射线追踪方法，通常利用Ｗａｌｄｅｎ近似　］，该公式基　于水平层状介质假设，无法得到相对于倾斜地层法　线的入射角．Ｒｅｓｎｉｃｋ［３］等分析了ＮＭＯ后、ＤＭＯ后　以及叠前时间偏移后ＣＲＰ道集ＡＶＯ效应，建议　Ｋｉｒｃｈｈｏｆｆ偏移同ＡＶＯ分析一起使用；Ｍｏｓｈｅｒ　等［４］导出平面波叠前时间偏移射线参数Ｐ道集，并　与ＣＭＰ道集相比较，展示了叠前时间偏移对于改　善地层分辨率、ＡＶＯ目标局部精度等方面的价值．　近年来，用于速度以及ＡＶＡ分析的角度域共　像点道集（ＡＤＣＩＧ）引起学者广泛关注　，ＡＤＣＩＧ　相对于偏移距共像点道集（ＯＤＣＩＧ）可以很大程度　上减少假象．这些研究多基于叠前深度偏移（ＰＳＤＭ），　与叠前时间偏移（ＰＳＴＭ）地表偏移距近似得到的入　射角相比，前者在速度精确情况下求得入射角为地　下局部入射角；但由于ＰＳＤＭ速度建模困难、计算　成本昂贵等原因，基于此的ＡＤＣＩＧ尚未见到应用　于三维ＡＶＡ反演实例．Ｋｉｒｃｈｈｏｆｆ积分法叠前时间　偏移对观测系统适应性强、计算成本低等特点使得　该方法一直在实际应用中占主导地位，以此输出的　角度道集适用于ＡＶＡ分析．ＺｈｅｎｇＥｎ３、Ｐｅｒｅｚ＿】幻提　出一种基于直射线近似Ｋｉｒｃｈｈｏｆｆ叠前时间偏移角　度域成像方法，并在高分辨成像和裂缝分析（方位　ＡＶＡｚ）得到应用；程玖兵提出走时计算基于射线追　踪的表驱Ｋｉｒｃｈｈｏｆｆ叠前时间偏移角度域算法ｌ】３＿以　及非双曲时差方位／角度域叠前时间偏移口　，走时　计算考虑射线弯曲效应以及等效横向各向同性，可　以获得相对准确的入射角、走时．Ｗａｎｇ等口　分析了　角度道集提取算法，计算走时采用适应长偏移距（大　入射角）优化六阶ＮＭＯ方程［】　，或者对线性横向　变速介质有更好聚焦性的非对称走时算法＿】　，采　用相对保真的加权因子以及覆盖次数校正技术；理　论模型和实际数据表明Ｋｉｒｃｈｈｏｆｆ叠前时间偏移角　度道集，可更准确获得偏移距与入射角映射关系，其　振幅相对保真，适于ＡＶＡ反演．　２　叠前时间偏移角度道集与ＣＭＰ、　ＣＲＰ道集对比　实际生产中，ＣＭＰ、ＣＲＰ／ＣＩＧ向角度道集转化　较多采用公式（１）Ｅ２３：　ａ＝＝＝ｓｉｎ－　（　）．　㈩　如图１所示，ａ为入射射线与时间轴之间夹角，　。　为层速度，ＶＲｍ　为均方根速度，ｔ为双程走时，ｚ为炮　检距，该公式仅在二维水平层状介质下成立，地层倾　斜时转化的角度不准确．图１中，地下界面水平时，　共深度点Ｍ即为反射点，但当地层倾斜时，反射点　在Ｎ处，分选ＣＭＰ道集将来自Ｎ处反射信息归为　Ｍ处ＣＭＰ道集，Ｎ点入射角　与水平层入射角ａ　不相同，随着地层倾角７值增大，二者偏离增大，Ｍ　点振幅值也将受到影响．Ｍｕｅｒｄｔｅｒ等　¨　数值分析　了不同地层倾角下偏移距与反射角变化，表明　ＡＶＡ分析必须考虑地层倾斜．除了角度因素外，在　图１　ＣＭＰ、ＣＲＰ、ＣＡＧ道集示意图　Ｆｉｇ．１　Ｃｏｎｔｒａｓｔ　ｂｅｔｗｅｅｎ　ＣＭＰ，ＣＲＰ　ａｎｄ　ＣＡＧ　ｇａｔｈｅｒｓ　５期　邹　振等：Ｋｉｒｃｈｈｏｆｆ叠前时间偏移角度道集　ＣＭＰ道集上，断层等横向不连续所产生断面波、绕　射波与目的层反射波互相干涉，以此基础上的　ＡＶＯ／ＡＶＡ分析不可信．因此，基于叠前偏移的　ＡＶＡ分析变得十分必要．叠前偏移使倾斜反射归　位到真正地下界面的位置，并使绕射波收敛，改善了　地震资料的横向分辨率和信噪比，消除界面曲率、倾　角对地震波相位与振幅的影响，为属性处理和反演　提供振幅相对保真的基础资料＿２　．　ＣＲＰ道集转化角度道集存在如下问题：地下构　造较复杂时会出现射线阴影区，阴影区中ＣＲＰ道集　存在假象　；ＣＲＰ道集不记录炮检点坐标，分选叠　加不考虑方位因素，当地层倾角较大时，根据　ｗａｌｄｅｎ公式转换入射角与真实值偏离较大；ＣＲＰ　道集按不同偏移距进行分选，近、中、远偏移距参与　叠加次数、偏移距分选问隔处理不当，都将造成振幅　异常假象，使ＡＶＡ走向误区．角度域叠前偏移算法　在偏移同时记录入射角／反射角信息，既保留ＣＲＰ　道集优势，又避开ＣＲＰ道集转化角度道集对水平层　状假设的依赖，合适的振幅加权因子、角度域覆盖次　数纠正技术使得偏移后的角度道集可应用于ＡＶＡ　反演；另外，角度道集上拉伸因子不变，可以较好地　消除影响成像分辨率、影响ＡＶＡ分析效果的远道　拉伸效应＿２　．　３叠前时间偏移角度道集提取算法　叠前时间偏移提取角度道集就是将偏移距改为　入射角／Ｋ射角作为分选单位，同一偏移距地震数据　被映射到不同角度中．图２ａ为直射线ＰＳＴＭ射线　路径图，计算走时采取ＤＳＲ双曲走时：　…　＝＝＝√（号）　＋熹＋√（号）。＋　，　（２）　ｒｓ－一￣／（　－一ｚｓ）　＋（　－一　ｓ）　（３）　ｒｉｇ一√（ｘｉ—ｚｇ）　＋（　。一　ｇ）　，　其中，ｔ。为零炮检距双程走时，　分别为炮点、　检波点到成像点水平距离，炮点ｓ坐标为（　，　。，　ｏ），检波点Ｒ坐标为（　，　，ｏ），成像点Ｉ坐标为　（　”　ｉ，ｔ　）．　入射角　为人射射线Ｓ　与反射射线ＪＲ夹角一　半，二者夹角可由余弦定理求出口引，该算法克服公　式（１）对水平层状介质假设依赖，可求得入射射线与　倾斜地层法线之间夹角．直射线法ＰＳＴＭ有如下不　足：采用双曲走时公式（２）仅保留偏移距ｒ的２次　项，随着ｒ增加，所计算走时与真实值偏差增大，使　得大偏移距／大入射角无法拉平；由于未考虑射线弯　曲效应，该算法所计算入射角范围对目的层照明不　足，严重影响大倾角成像以及ＡＶＡ反演精度．工业　界多采用大偏移距采集，陆上最大偏移距一般可达　到５０００￣６０００　Ｉｎ，甚至更高，直射线法在走时以及　入射角计算精度上都不能达到要求．Ｒｏｓｓｌ２引、Ｂａｌｅ　等＿２　提出使用非双曲走时进行ＮＭＯ校正与ＡＶＯ　相结合，可以提高ＡＶＯ反演的精度．我们将以上思　路结合，拓展到Ｋｉｒｃｈｈｏｆｆ叠前时间偏移中，计算走　时采用非双曲走时方程，如优化６阶ＮＭＯ方程，适　应线性横向变速的非对称走时算法等，根据射线参　数法并考虑偏移距高阶项计算入射角，因此大角度　处理更精确．　我们以优化６阶ＮＭｏ方程为例口　，说明非双　曲走时（弯曲射线）角度道集提取策略，公式（４）为双　平方根（ＤＳＲ）优化６阶ＮＭＯ走时方程：　一　＋￡　≈　３　４－∞　＋　３ｇ＋ｃｃ　Ｃ＿４Ｙｉｇ，（４）　３ｓ　３ｇ　其中，　。　一￣／ｃ　＋ｃ　，．　＋　。ｒ　，ｔ　一￣／ｃ　＋ｃ　ｒ　＋４ｆ。ｒ　，　ｃｌ，Ｃ２，Ｃ３，Ｃ４与Ｔａｎｅｒ等＿２　给出的前四项系数相　同，ＣＣ为常数．根据Ｓｎｅｌｌ定理，即水平射线参数　不变原理，对炮点到像点走时ｔ　以及检波点到像点　走时ｔ　作ｒ方向求导得：　ｊ‘　’　　Ｔｓ　＋８ｃ３ｒ。３　４－６ｃｃ＊ｆ４　５）　Ｕ，ｓｉ　一＿１（３　ｃｃ　ｃ４ｒ６　（ｃ２ｒｓｉ＋８ｃ３ｒｓ３。，）　；。　’　一去　＋８ｃｓｒ　＋６　ｒ　）　一　！！　！！　堡±　！！　，　（５）　ｊ　一ｓ　ｔ　），　㈤　ａｚ＝ｓｉ一（　ｔ　）．　实际计算Ｏｔ　、ａ　时，公式（５）保留到ｒ的５次项即可，　可节约计算量，而该方法对以度为计量单位的入射　角精度影响甚微，本文模型以及实际数据都采用这　种截断，角度之间采用交替叠加以提高信噪比．如图　２ｂ所示，由公式（６）分别求得射线ＭＩ与时间轴ｔ。　夹角ａ　以及射线Ｎ　与时间轴夹角ａ。，入射角　可　由空间向量法则求得　：　１２１４　地球物理学报（Ｃｈｉｎｅｓｅ　Ｊ．Ｇｅｏｐｈｙｓ．）　５３卷　ｍｉｇｒａｔｉｏｎ　ｏｎ　ＡＶｏ．Ｇｅｏｐｈｙｓｉｃｓ，１９９６，６１（６）：１６０３～１６１５　［５］　Ｐｒｕｃｈａ　Ｍ　Ｌ，Ｂｉｏｎｄｉ　Ｂ　Ｌ，Ｓｙｍｅｓ　Ｗ　Ｗ．Ａｎｇｌｅ—ｄｏｍａｉｎ　ｉｍａｇｅ　ｇａｔｈｅｒｓ　ｂｙ　ｗａｖｅ—ｅｑｕａｔｉｏｎ　ｍｉｇｒａｔｉｏｎ：Ａｎｇｌｅ　ｄｏｍａｉｎ　ｃｏｍｍｏｎ　ｉｍａｇｅ　ｇａｔｈｅｒｓ　ｂｙ　ｗａｖｅ　ｅｑｕａｔｉｏｎ．６９　Ａｎｎｕａｌ　Ｍｅｅｔｉｎｇ，ＳＥＧ，　Ｅｘｐａｎｄｅｄ　Ａｂｓｔｒａｃｔｓ．１９９９．８２４～８２７　［６］　Ｍｏｓｈｅｒ　Ｃ　Ｃ，Ｊｉｎ　Ｓ，Ｆｏｓｔｅ　Ｄ　Ｊ．Ｍｉｇｒａｔｉｏｎ　ｖｅｌｏｃｉｔｙ　ａｎａｌｙｓｉｓ　ｕｓｉｎｇ　ｃｏｍｍｏｎ　ａｎｇｌｅ　ｉｍａｇｅ　ｇａｔｈｅｒｓ．　７１ｓｔ　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｍｅｅｔｉｎｇ，ＳＥＧ，Ｅｘｐａｎｄｅｄ　Ａｂｓｔｒａｃｔｓ，２００１．　８８９～８９２　［７］　Ｘｕ　Ｓ，Ｃｈａｕｒｉｓ　Ｈ，Ｌａｍｂａｒｅ　Ｇ，ｅｔ　ａ１．Ｃｏｍｍｏｎ－ａｎｇｌｅ　ｍｉｇｒａｔｉｏｎ：　Ａ　ｓｔｒａｔｅｇｙ　ｆｏｒ　ｉｍａｇｉｎｇ　ｃｏｍｐｌｅｘ　ｍｅｄｉａ．　Ｇｅｏｐｈｙｓｉｃｓ，２００１，６６（６）：１８７７——１８９４　［８］　Ｓａｖａ　Ｐ　Ｃ，Ｆｏｍｅｌ　Ｓ．Ａｎｇｌｅ—ｄｏｍａｉｎ　ｃｏｍｍｏｎ—ｉｍａｇｅ　ｇａｔｈｅｒｓ　ｂｙ　ｗａｖｅｆｉｅｌｄ　ｃｏｎｔｉｎｕａｔｉｏｎ　ｍｅｔｈｏｄｓ．Ｇｅｏｐｈｙｓｉｃｓ，２００３，６８（３）：　１０６５～１０７４　［９］　Ｓｔｏｌｋ　Ｃ　Ｃ，ｄｅ　Ｈｏｏｐ　Ｍ　Ｖ，Ｓｙｍｅｓ　Ｗ　Ｗ．Ｋｉｎｅｍａｔｉｃｓ　ｏｆ　ｐｒｅｓｔａｃｋ　ｓｈｏｔ—ｇｅｏｐｈｏｎｅ　ｍｉｇｒａｔｉｏｎ．７５ｔｈ　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　ｍｅｅｔｉｎｇ，ＳＥＧ，Ｅｘｐａｎｄｅｄ　ａｂｓｔｒａｃｔｓ，２００５．１８６６－－１８６９　［１０］　陈凌，吴如山，王伟君．基于Ｇａｂｏｒ　Ｄａｕｂｅｃｈｉｅｓ小波束叠前　深度偏移的角度域共成像道集．地球物理学报，２００４，４７（５）：　８７６～８８５　Ｃｈｅｎ　Ｌ，Ｗｕ　Ｒ　Ｓ，Ｗａｎｇ　Ｗ　Ｊ．Ｃｏｍｍｏｎ　ａｎｇｌｅ　ｉｍａｇｅ　ｇａｔｈｅｒｓ　ｏｂｔａｉｎｅｄ　ｆｒｏｍ　Ｇａｂｏｒ　ｄａｕｂｅｃｈｉｅｓ　ｂｅａｍｌｅａｔ　ｐｒｅｓｔａｃｋ　ｄｅｐｔｈ　ｍｉｇｒａｔｉｏｎ．Ｃｈｉｎｅｓｅ　Ｊ．Ｇｅｏｐｈｙｓ．（ｉｎ　Ｃｈｉｎｅｓｅ），２００４，４７（５）：　８７６～８８５　［１１］　Ｚｈｅｎｇ　Ｙ．Ｓｅｉｓｍｉｃ　ａｚｉｍｕｔｈａｌ　ａｎｉｓｏｔｒｏｐｙ　ａｎｄ　ｆｒａｃｔｕｒｅ　ａｎａｌｙｓｉｓ　ｆｒｏｍ　ＰＰ　ｒｅｆｌｅｃｔｉｏｎ　ｄａｔａ［Ｐｈ．Ｄ．ｔｈｅｓｉｓ］．Ｃａｌｇａｒｙ：Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　Ｃａｌｇａｒｙ，２００６　［１２］　Ｐｅｒｅｚ　Ｇ，Ｍａｒｆｕｒｔ　Ｋ　Ｊ．Ｉｍｐｒｏｖｉｎｇ　ｌａｔｅｒａｌ　ａｎｄ　ｖｅｒｔｉｃａｌ　ｒｅｓｏｌｕｔｉｏｎ　ｏｆ　ｓｅｉｓｍｉｃ　ｉｍａｇｅｓ　ｂｙ　ｃｏｒｒｅｃｔｉｎｇ　ｆｏｒ　ｗａｖｅｌｅｔ　ｓｔｒｅｔｃｈ　ｉｎ　ｃｏｍｍｏｎ－ａｎｇｌｅ　ｍｉｇｒａｔｉｏｎ．Ｇｅｏｐｈｙｓｉｃｓ，２００７，７２（６）：９４～１０４　［１３］　程玖兵，王楠，马在田．表驱三维角度域Ｋｉｒｃｈｈｏｆｆ叠前时　间偏移成像方法．地球物理学报，２００９，５２（３）：７９２～８００　Ｃｈｅｎｇ　Ｊ　Ｂ，Ｗａｎｇ　Ｎ，Ｍａ　Ｚ　Ｔ．ｒｒａｂ１ｅ　ｄｒｉｖｅｎ　３　Ｄ　ａｎｇｌｅ－ｄｏｍａｉｎ　ｉｍａｇｉｎｇ　ａｐｐｒｏａｃｈ　ｆｏｒ　Ｋｉｒｃｈｈｏｆｆ　ｐｒｅｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ．Ｃｈｉｎｅｓｅ　Ｊ．Ｇｅｏｐｈｙｓ．（ｉｎ　Ｃｈｉｎｅｓｅ），２００９，５２（３）：７９２￣８００　［１４］　程玖兵，王　楠．非双曲时差方位／角度域叠前时间偏移．　ＣＰＳ／ＳＥＧ　Ｂｅｉｊｉｎｇ　２００９　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｇｅｏｐｈｙｓｉｃａｌ　Ｃｏｎｆｅｒｅｎｃｅ　＆Ｅｘｐｏｓｉｔｉｏｎ．２００９．ＩＤ１　１　７９　Ｃｈｅｎｇ　Ｊ　Ｂ，Ｗａｎｇ　Ｎ．Ａｚｉｍｕｔｈ／ａｎｇｌｅ　ｄｏｍａｉｎ　ｉｍａｇｉｎｇ　ｉｎ　ｐｒｅｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ　ｕｓｉｎｇ　ｎｏｎ—ｈｙｐｅｒｂｏｌｉｃ　ｍｏｖｅｏｕｔ．ＣＰＳ／　ＳＥＧ　Ｂｅｉｊｉｎｇ　２００９　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｇｅｏｐｈｙｓｉｃａｌ　Ｃｏｎｆｅｒｅｎｃｅ　８Ｌ　Ｅｘｐｏｓｉｔｉｏｎ．２００９．ＩＤ１１７９　［１５］　Ｗａｎｇ　Ｄ，Ｚｈｅｎｇ　Ｘ　Ｄ，Ｃｈｅｎ　Ｊ　Ｂ．Ａｍｐｌｉｔｕｄｅ—ｐｒｅｓｅｒｖｉｎｇ　ｐｌａｎｅ—　ｗａｖｅ　ｐｒｅｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ　ｆｏｒ　ＡＶＯ　ａｎａｌｙｓｉｓ．Ａｐｐｌｉｅｄ　Ｇｅｏｐｈｙｓｉｃｓ，２００８，５（３）：２１２～２１８　［Ｉ６］　Ｓｕｎ　Ｃ　Ｗ．Ｏｐｔｉｍｉｚｅｄ　６ｔｈ　ＮＭＯ　ｃｏｒｒｅｃｔｉｏｎ　ｆｏｒ　ｌｏｎｇ—ｏｆｆｓｅｔ　ｓｅｉｓｍｉｃ　ｄａｔａ．７２ｔｈ　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｍｅｅｔｉｎｇ，ＳＥＧ，　Ｅｘｐａｎｄｅｄ　Ａｂｓｔｒａｃｔｓ，２００２　［Ｉ７］　刘洪，王秀闽，曾　锐等．单程波算子积分解的象征表示．　地球物理学进展，２００７，２２（２）：４６３～４７１　Ｌｉｕ　Ｈ，Ｗａｎｇ　Ｘ　Ｍ，Ｚｅｎｇ　Ｒ，ｅｔ　ａ１．Ｓｙｍｂｏｌ　ｄｅｓｃｒｉｐｔｉｏｎ　ｔｏ　ｉｎｔｅｇｒａｌ　ｓｏｌｕｔｉｏｎ　ｏｆ　ｏｎｅ—ｗａｙ　ｗａｖｅ　ｏｐｅｒａｔｏｒ．Ｐｒｏｇｒｅｓｓ　ｉｎ　Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ　Ｃｈｉｎｅｓｅ），２００７，２２（２）：４６３～４７１　［１８］李博，刘国峰，刘洪．地震叠前时间偏移的一种图形提速　实现方法．地球物理学报，２００９，５２（１）：１～８　Ｌｉ　Ｂ，Ｌｉｕ　Ｇ　Ｆ，Ｌｉｕ　Ｈ．Ａ　ｍｅｔｈｏｄ　ｏｆ　ｕｓｉｎｇ　ＧＰＵ　ｔｏ　ａｃｃｅｌｅｒａｔｅ　ｓｅｉｓｍｉｃ　ｐｒｅ—ｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ．Ｃｈｉｎｅｓｅ　Ｊ．Ｇｅｏｐｈｙｓ．（ｉｎ　Ｃｈｉｎｅｓｅ），２００９，５２（１）：１～８　［１９］Ｍｕｅｒｄｔｅｒ　Ｄ，ＬＬＣ　Ｌ，Ｋｅｌｌｙ　Ｍ，ｅｔ　ａ１．Ｔｈｅ　ｅｆｆｅｃｔ　ｏｆ　ｒｅｆｌｅｃｔｏｒ　ｄｉｐ　ｏｎ　ＡＶ０　ａｎａｌｙｓｉｓ．７５　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｍｅｅｔｉｎｇ，　ＳＥＧ，Ｅｘｐａｎｄｅｄ　Ａｂｓｔｒａｃｔｓ，２００５　［２Ｏ］刘洪林，朱秋影．基于叠前深度偏移的ＡＶＯ反演及解释．地　球物理学进展，２００７，２２（３）：９０５～９１３　Ｌｉｕ　Ｈ　Ｌ，Ｚｈｕ　Ｑ　Ｙ．ＡＶＯ　ｉｎｖｅｒｓｉｏｎ　ａｎｄ　ｉｎｔｅｒｐｒｅｔａｔｉｏｎ　ｂａｓｅｄ　ｏｎ　ｐｒｅｓｔａｅｋ　ｄｅｐｔｈ　ｍｉｇｒａｔｉｏｎ．Ｐｒｏｇｒｅｓｓ　ｉｎ　Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ　Ｃｈｉｎｅｓｅ），２００７，２２（３）：９０５～９１３　［２１］胡贤根，尚新民，石林光等．基于叠前时间偏移资料的ＡＶＯ　处理技术．石油物探，２００２，４１（３）：３４３～３４６　Ｈｕ　Ｙ　Ｇ，Ｓｈａｎｇ　Ｘ　Ｍ，Ｓｈｉ　Ｌ　Ｇ，ｅｔ　ａ１．ＡＶＯ　ｐｒｏｃｅｓｓ　ｔｅｃｈｎｉｑｕｅ　ｂａｓｅｄ　ｏｎ　ｐｒｅｓｔａｃｋ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ　ｄａｔａ．Ｇｅｏｐｈｙｓｉｃａｌ　Ｐｒｏｓｐｅｃｔｉｎｇ　＿厂０ｒ　Ｐｅｔｒｏｌｅｕｍ（ｉｎ　Ｃｈｉｎｅｓｅ），２００２，４１（３）：３４３￣３４６　［２２］ｒｒｙｇｅｌ　Ｍ，Ｓａｎｔｏｎｓ　Ｌ　Ｔ，Ｓｃｈｌｅｉｃｈｅｒ　Ｊ．Ｋｉｒｃｈｈｏｆｆ　ｉｍａｇｉｎｇ　ａｓ　ａ　ｔｏｏｌ　ｆｏｒ　ＡＶＯ／ＡＶＡ．Ｔｈｅ　Ｌｅａｄｉｎｇ　Ｅｄｇｅ，１９９９，１８（８）：９４０～９４５　［２３］Ａｖｓｅｔｈ　Ｐ，Ｍｕｋｅｒｊｉ　Ｔ，Ｍａｙｋｏ　Ｇ，Ｑｕａｎｔｉｔａｔｉｖｅ　Ｓｅｉｓｍｉｃ　Ｉｎｔｅｒｐｒｅｔａｔｉｏｎ．Ｃａｍｂｒｉｄｇｅ　Ｕｎｉｖｅｒｓｉｔｙ　Ｐｒｅｓｓ，２００５．１９Ｏ～２００　［２４］Ｒｏｙ　Ｂ，ＡＮｍｏ　Ｐ　Ｄ，Ｂａｕｍｅｌ　Ｒ，ｅｔ　ａ１．Ａｎａｌｙｔｉｃ　ｃｏｒｒｅｃｔｉｏｎ　ｆｏｒ　ｗａｖｅｌｅｔ　ｓｔｒｅｔｃｈ　ｄｕｅ　ｔｏ　ｉｍａｇｉｎｇ．７５　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｍｅｅｔｉｎｇ，ＳＥＧ，Ｅｘｐａｎｄｅｄ　Ａｂｓｔｒａｃｔｓ，２００５．２３４～２３７　［２５］Ｒｏｓｓ　Ｃ．ＡＶＯ　ａｎｄ　ｎｏｎｈｙｐｅｒｂｏｌｉｃ　ｍｏｖｅｏｕｔ：Ａ　ｐｒａｃｔｉｃａｌ　ｅｘａｍｐｌｅ．Ｆｉｒｓｔ　Ｂｒｅａｋ，１９９７，１５：４３～４８　［２６］Ｂａｌｅ　Ｒ，Ｌｅａｎｅｙ　Ｓ，Ｄｕｍｉｔｒｕ　Ｇ．Ｏｆｆｓｅｔ—ｔｏ－ａｎｇｌｅ　ｔｒａｎｓｆｏｒｍａｔｉｏｎｓ　ｆｏｒ　ＰＰ　ａｎｄ　ＰＳ　ＡＶＯ　ａｎａｌｙｓｉｓ．７１　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｍｅｅｔｉｎｇ，　ＳＥＧ，Ｅｘｐａｎｄｅｄ　Ａｂｓｔａｃｔｓ，２００１．２３５～２３８　［２７］Ｔａｎｅｒ　Ｍ　Ｔ，Ｋｏｅｈｌｅｒ　Ｆ．Ｖｅｌｏｃｉｔｙ　ｓｐｅｃｔｒａ　ｄｉｇｉｔａｌ　ｃｏｍｐｕｔｅｒ　ｄｅｒｉｖａｔｉｏｎ　ａｎｄ　ａｐｐｌｉｃａｔｉｏｎｓ　ｏｆ　ｖｅｌｏｃｉｔｙ　ｆｕｎｃｔｉｏｎｓ．Ｇｅｏｐｈｙｓｉｃｓ，　１９６９，３４（６）：８５９～８８１　［２８］Ｃａｓｔａｌ　Ｒ　Ｊ．Ａ　ｔｈｅｏｒｙ　ｏｆ　ｎｏｒｍａｌ　ｍｏｖｅｏｕｔ．Ｇｅｏｐｈｙｓｉｃｓ，１９９４，　５９（６）：９８３～９９９　［２９］Ｚｏｕ　Ｚ，Ｌｉｕ　Ｈ，Ｌｉｕ　Ｗ　Ｈ．３　Ｄ　ａｎｇｌｅ—ｄｏｍａｉｎ　ｃｏｍｍｏｎ　ｉｍａｇｉｎｇ　ｇａｔｈｅｒｓ　ｄｕｒｉｎｇ　Ｋｉｒｃｈｈｏｆｆ　ＰＳＴＭ．７９　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｍｅｅｔ　ｎｇ，ＳＥＧ，Ｅｘｐａｎｄｅｄ　Ａｂｓｔｒａｃｔｓ，２００９．３２６～３３０　［３Ｏ］Ｚｈａｎｇ　Ｙ，Ｇｒａｙ　Ｓ，Ｙｏｎｇ　Ｊ．Ｅｘａｃｔ　ａｎｄ　ａｐｐｒｏｘｉｍａｔｅ　ｗｅｉｇｈｔｓ　ｆｏｒ　Ｋｉｒｃｈｈｏｆｆ　ｍｉｇｒａｔｉｏｎ．７０　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｍｅｅｔｉｎｇ，　ＳＥＧ，Ｅｘｐａｎｄｅｄ　Ａｂｓｔａｃｔｓ，２０００．１０３９～１０３９　［３１］Ｌｅｅ　Ｓ，Ｋｉｎｇ　Ｄ，Ｌｉｎ　Ｓ．Ｅｆｆｉｃｉｅｎｔ　ｔｒｕｅ—ａｍｐｌｉｔｕｄｅ　ｗｅｉｇｈｔｓ　ｉｎ　Ｋｉｒｃｈｈｏｆｆ　ｔｉｍｅ　ｍｉｇｒａｔｉｏｎ．７４　Ａｎｎｕａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｍｅｅｔｉｎｇ，　ＳＥＧ。Ｅｘｐａｎｄｅｄ　Ａｂｓｔｒａｃｔｓ，２００４．１０８９～１０９２　［３２］Ｄｅｌｌｉｎｇｅｒ　Ｊ　Ａ，Ｇｒａｙ　Ｓ　Ｈ，Ｍｕｒｐｈｙ　Ｇ　Ｅ，ｅｔ　ａ１．Ｅｆｆｉｃｉｅｎｔ　２．５－Ｄ　ｔｒｕｅ－ａｍｐｌｉｔｕｄｅ　ｍｉｇｒａｔｉｏｎ．Ｇｅｏｐｈｙｓｉｃｓ，２０００，６５（３）：９４３～９５０　［３３］Ａｕｄｅｂｅｒｔ　Ｆ，Ｎｉｃｏｌｅｔｉｓ　Ｌ，Ｆｒｏｉｄｅｖａｕｘ　Ｐ．Ｒｅｇｕｌａｒｉｚａｔｉｏｎ　ｏｆ　ｉｌｌｕｍｉｎａｔｉｏｎ　ｉｎ　ａｎｇｌｅ　ｄｏｍａｉｎｓ－－Ａ　ｋｅｙ　ｔｏ　ｔｒｕｅ　ａｍｐｌｉｔｕｄｅ　ｍｉｇｒａｔｉｏｎ．ＴｈｅＬｅａｄｉｎｇ　Ｅｄｇｅ，２００５，２４（６）：６４３￣６５４　（本文编辑胡素芳）