油储地球物理会议-三维逆时深度偏移成像及3D VSP应用(NXPowerLite).ppt

三维逆时深度偏移成像及在3D VSP 数据上的应用油储地球物理研讨会暨祝贺李幼铭老师从事地球物理事业50周年 王华忠 同济大学海洋与地球科学学院 2010年6月21号 内容 诚挚感谢李幼铭老师 给李幼铭老师的简单工作结果汇报 逆时深度偏移方法 逆时深度偏移在三维地面地震数据上的应用 逆时深度偏移在三维VSP数据上的应用 结论与讨论 祝贺李幼铭老师从事地球物理事业五十周年。诚挚感谢李幼铭老师对年轻地球物理工作者的培养。油储地球物理自然科学基金重大项目十几年间培养了目前勘探地球物理研究领域的中坚力量。李幼铭老师功不可没！ 给李幼铭老师的简单工作结果汇报 千年之交，油储地球物理重大项目正处在“十五”攻关时期。 李老师提出：中科院地球物理所团队重点做积分法偏移成像；同济大学团队重点做波动方程叠前深度成像。 给李幼铭老师的简单工作结果汇报 十年左右的时间，同济大学团队在叠前成像方面主要做的工作体现在同济大学和胜利油田物探研究院共同开发的“神通STseis”叠前地震处理成像系统和同济大学与南京石油物探研究院共同开发的“iCluster”叠前地震成像系统中。 给李幼铭老师的简单工作结果汇报 主要工作包括： 三维叠后时间与深度偏移 Kirchhoff积分法 波动方程 时间空间差分 频率空间差分 Kirchhoff三维叠前时间与深度偏移 水平地表 非水平地表 各向异性介质 给李幼铭老师的简单工作结果汇报 主要工作包括： 波动方程叠前时间与深度偏移---单向波方程 （粘）声波方程叠前时间偏移成像---有限差分方法 （粘）声波方程单炮道集叠前深度偏移成像 各向异性介质拟单波单炮道集叠前深度偏移成像 双平方根算子CMP道集、炮集共方位角数据及全三维深度成像 最小二乘深度偏移成像 三维逆时深度偏移成像 数值试验 体偏移结果与Omega结果 体偏移输出成像道集 体偏移与Omega效率对比 体偏移技术扩展---射线束偏移 单个射线束偏移： 体偏移技术扩展---射线束偏移 数值试验 数值试验 体偏移技术扩展---弯曲射线偏移 数值试验 数值试验 偏移距平面波时间偏移 三维偏移距平面波分解 三维偏移距平面波偏移（声波） 三维偏移距平面波偏移（粘声） Feng, B., Wang, H., Liu, S., 2009, 3D offset plane-wave finite-difference pre-stack time migration: 79th Annual International Meeting, SEG, Expanded Abstracts, 2934-2938. 偏移距平面波分解 二维合成数据实验 子波频谱对比 计算效率对比 数值结果 数值结果 数值结果 吸收衰减补偿全三维叠前深度偏移 再次感谢李幼铭老师对年轻地球物理研究人员的培养。李老师在十多年的时间里为年轻地球物理研究人员搭建了一个安心发展的平台。 三维逆时深度---内容 概述 波场外推算子---三维差分格式 成像条件与低频噪音 计算策略、边界条件与存储 三维VSP逆时深度偏移 结论与讨论 逆时偏移成像---概述 逆时偏移成像---概述 它包括三部分关键步骤： 1）沿时间正向延拓震源波场； 2）沿时间反方向延拓接收点波场； 3）应用适当成像条件产生偏移成像结果。 逆时偏移成像---概述 逆时偏移成像---概述 目前存在的主要问题： 对于复杂的构造和散射体成像缺乏实用的成像条件。 常规的单向波成像条件：下行波的到达时等于上行波的出发时。Claerbout（1971）提出。 很容易实现，但存在很强的低频噪音。 双向波成像条件：反射界面处入射波场到达时等于反射波场出发时。修改后的激励时间成像条件 。 不容易实现，但可以很好地压制低频噪音。 逆时偏移成像---概述 逆时偏移成像---概述 内容 概述 波场外推算子---三维差分格式 成像条件与低频噪音 计算策略、边界条件与存储 三维VSP逆时深度偏移 结论与讨论 高阶差分波场外推算子 逆时偏移的主要计算问题就是双向波方程的差分格式和边界条件。 三维声波方程： 高阶差分波场外推算子 高阶差分格式 炮点波场外推 正向时间外推 高阶差分波场外推算子 波场外推算子主要关注精度和效率的平衡： 高精度要求网格密度大；网格密度大计算效率低。 显式差分格式不是无条件稳定的。稳定的格式对网格密度有要求。 目前，精度更高的是紧致差分格式；粘滞差分格式等。 内容 概述 波场外推算子---三维差分格式 成像条件与低频噪音 计算策略、边界条件与存储 三维VSP逆时深度偏移 结论与讨论 成像条件与低频噪音 波场分解成像条件 Liu, et al（2007） 坡印廷矢量成像条件 Yoon, et al. (2006)