地震资料叠前反演技术应用研究

地震资料叠前反演技术应用研究中油集团东方地球物理公司研究院2007年8月报告人：冉建斌中油集团东方地球物理公司研究院处理解释技术研讨会汇报内容前言地震油气检测技术发展地震叠前反演关键技术结论及建议中油集团东方地球物理公司研究院处理解释技术研讨会前言随着计算机运算能力不断提高、叠前地震资料处理技术的发展及弹性阻抗概念的提出（Ptrick，1999），推动了以叠后地震资料为基础的含油气储层的描述向以叠前叠前资料处理为基础的含油气储层描述的发展。叠后和叠前地震资料差异？中油集团东方地球物理公司研究院处理解释技术研讨会前言地震油气检测技术发展地震叠前反演关键技术结论及建议汇报内容中油集团东方地球物理公司研究院处理解释技术研讨会S1R1R2R3R4H实际观测S2多次覆盖S1S2S3R1R1R1地震资料油气检测技术的发展中油集团东方地球物理公司研究院处理解释技术研讨会),,,,,,(222111SpSpPPVVVVfR1212)(AIAIAIAIR垂直入射垂直入射及随偏移距变化反射系数与弹性参数的关系)sin(tansin)(:sinsin1)(:&2ln)(:22212)1(024cos212421)/ln()/ln(22122sincoscos2122222212121211122PPsspsPPPSsspppVVPPVVVVVVVVVVssVVVNIANIRCShueyRCRichardsAkiVVRCBortfeld近中远角度关系岩性关系流体项刚性项地震资料油气检测技术的发展中油集团东方地球物理公司研究院处理解释技术研讨会气水界面含气砂岩模型含气砂岩AVO正演Ostrander(1982)萨克拉曼多(美国加州首府)储层在6700ft,1,75s白垩系的深海扇沉积最大净产层厚度95ft。AVO效应被证实中油集团东方地球物理公司研究院处理解释技术研讨会零偏移距合成道零偏移距真实数据道偏移距振幅AVO被广泛应用中油集团东方地球物理公司研究院处理解释技术研讨会AVO截距AVO梯度2sinPPPRRGAVO被广泛应用中油集团东方地球物理公司研究院处理解释技术研讨会-0.1+0.1+0.2-0.210入射角反射系数20301类3类2类高阻抗油藏接近零阻抗油藏低阻抗油藏入射角亮点极性反转暗点实际CRP道集合成道集AVO类型图解油藏顶面油藏顶面油藏顶面油藏顶面油藏顶面小大亮点暗点小大小大AVO异常被用于检测油气的存在AVO被广泛应用中油集团东方地球物理公司研究院处理解释技术研讨会弹性阻抗反演是BP公司Connolly.Ptrick，1999年提出来的将波阻抗的概念推广到非零入射角，通过推导得到了非零入射角的反射系数与垂直入射角时反射系数公式具有相同的形式，此时对应波阻抗的物理量称为弹性阻抗，其表达形式如下：)()()()()()(12121212EIEIEIEIRVAIAIAIAIAIRP设K=Vs2∕Vp2≈0.21EI=AI=声阻抗＝纵波速度×密度弹性阻抗的提出(1999年)中油集团东方地球物理公司研究院处理解释技术研讨会弹性阻抗反演的改进中油集团东方地球物理公司研究院处理解释技术研讨会弹性阻抗反演及应用－岩性流体识别弹性阻抗对含油气层响应明显，而纵波阻抗没有反映中油集团东方地球物理公司研究院处理解释技术研讨会纵波阻抗对岩性和流体的区分弹性阻抗对岩性和流体的区分弹性阻抗反演岩性流体识别中油集团东方地球物理公司研究院处理解释技术研讨会低角度叠加声阻抗标定高角度叠加弹性阻抗标定地震的AVO效应能够被标定成为现实中油集团东方地球物理公司研究院处理解释技术研讨会波阻抗反演结果弹性阻抗反演结果弹性阻抗反演－岩性流体识别中油集团东方地球物理公司研究院处理解释技术研讨会合成数据由反射系数和子波褶积来模拟。AVA反射系数和弹性参数之间的关系可以由Knott-Zoeppritz方程或Aki-Richards近似式来求解。Knott-Zoeppritz方程：Aki-Richards近似式：叠前弹性参数反演中油集团东方地球物理公司研究院处理解释技术研讨会反演算法AVA约束的稀疏脉冲反演模块是基于这样的一种算法，对每一道来求解下列约束的优化问题：最小目标函数：F=Freflectivty＋Fcontrast＋Fseismic＋Ftrend＋Fspatial＋FGardner＋FMudrock所给的约束：Vpi（下限）＜Vpi＜Vpi（上限）Vsi（下限）＜Vsi＜Vsi（上限）ρi（下限）＜ρi＜ρi（上限）这里的Fx是各项误差的范数叠前弹性参数反演中油集团东方地球物理公司研究院处理解释技术研讨会0-10010-20020-300部分叠加角道集叠加中油集团东方地球物理公司研究院处理解释技术研讨会19-31°弹性阻抗与纵波阻抗曲线叠合图19-31°弹性阻抗与纵波阻抗反射系数叠合图用弹性阻抗和Aki计算反射系数交会图弹性阻抗及反射系数的求取中油集团东方地球物理公司研究院处理解释技术研讨会不同角道集AVO子波提取及层位标定老堡南1井全叠加资料子波及层位标定老堡南1井0-11°资料子波及层位标定老堡南1井9-21°资料子波及层位标定老堡南1井19-31°资料子波及层位标定中油集团东方地球物理公司研究院处理解释技术研讨会老堡南1井近、中、远子波叠合图研究区5口井近、中、远子波叠合图不同角道集AVA子波提取分析中油集团东方地球物理公司研究院处理解释技术研讨会信噪比、井相关和Lambda迭代次数等优化信噪比、井相关、标准化、井标准偏差系数脉冲数及综合误差与Lamda关系中油集团东方地球物理公司研究院处理解释技术研讨会泊松比纵横波速度比横波速度叠后波阻抗泊松比横波速度25°弹性阻抗纵横波速度比反演结果分析中油集团东方地球物理公司研究院处理解释技术研讨会1、储层的岩性和流体的识别储层表征（碎屑岩和碳酸岩）2、用4-D地震进行油藏压力和饱和度的估算，油藏监测3、应用叠前弹性阻抗反演和井的合成记录进行含气识别02040608010012014016018020019992001200320052007Joson:2000Paradigm2001Hampson-RussellEPTCGG...叠前反演应用及商业软件推出中油集团东方地球物理公司研究院处理解释技术研讨会要作好叠前地震资料岩性及流体预测，我们需要怎样工作？地震资料处理要求？如何一体化？需要发展哪些技术？...？国内－复杂含油气储层中油集团东方地球物理公司研究院处理解释技术研讨会汇报内容前言地震油气检测技术发展地震叠前反演关键技术结论及建议中油集团东方地球物理公司研究院处理解释技术研讨会南堡5号南堡2号南堡3号南堡4号据冀东油田公司海岸线（水深＜20m）125Km2南堡1号基于AVO分析三高处理监控研究区位置及地表条件工区概况及研究难点中油集团东方地球物理公司研究院处理解释技术研讨会Ed12715-2847m厚度：53.6m自喷99.61吨天Ng2116.2-2183.411m排液59.04吨天Ng和Ed1厚度20.4m自喷260.91吨天断裂复杂、火成岩发育和砂泥岩薄互层TRACE576SN火成岩中油集团东方地球物理公司研究院处理解释技术研讨会地震资料采集情况1中油集团东方地球物理公司研究院处理解释技术研讨会预处理几何扩散补偿道间能量均衡废坏道剔除野外静校正处理反Q滤波Q值试验带通滤波地表一致性反褶积第一次速度分析动校正确定切除参数地表一致性剩余静校正第二次速度分析地表一致性振幅校正动校正叠前去噪处理叠前剩余振幅补偿叠前偏移处理OBC处理地震资料处理及质量监控中油集团东方地球物理公司研究院处理解释技术研讨会陆检水检双检合并OBC地震资料处理及监控1一中油集团东方地球物理公司研究院处理解释技术研讨会陆检水检双检求和OBC地震资料处理及监控1中油集团东方地球物理公司研究院处理解释技术研讨会NP2-2LPN-1NP2-2LPN-1NP2-2LPN-1NP2-2LPN-12叠前配套处理技术应用及效果对比地表一致性反褶积地表一致性剩余静校正压制异常振幅、面波和多次波叠前三维随机噪音衰减中油集团东方地球物理公司研究院处理解释技术研讨会地表一致性反褶积后叠加叠前高信噪比处理及叠加效果对比Inline230地表一致性静校正后叠加叠前三维随机噪音衰减后叠加叠前压制多次波后叠加2中油集团东方地球物理公司研究院处理解释技术研讨会叠前去多次波后叠前时间偏移叠加三维去多次波＋随机噪音衰减后叠前时间偏移叠加3叠前高信噪比处理及偏移叠加效果对比中油集团东方地球物理公司研究院处理解释技术研讨会相干体技术用作地震资料品质监控手段1800ms4中油集团东方地球物理公司研究院处理解释技术研讨会三维随机噪音衰减有效带宽压制多次波后有效带宽CMP分选后的有效带宽地表一致性剩余剩余静校正有效带宽属性地震信噪比监控多道倾角滤波后叠前时间偏移三维随机噪音衰减后叠前时间偏移5中油集团东方地球物理公司研究院处理解释技术研讨会不同偏移距动校时差分析－速度精度及密度过NP203-NP2-3-NP2-2-NP2-1-LPN1井任意线9-21°－19-31°部分叠加叠合剖面6中油集团东方地球物理公司研究院处理解释技术研讨会中偏移距－远偏移距地震信息相关系数中偏移距－远偏移距地震信息可预测性中偏移距－远偏移距馆陶组底界时差分布中偏移距－远偏移距地震信息最小均方根方差不同偏移距地震资料的可预测性监控标准均方根残差相关值的平方7中油集团东方地球物理公司研究院处理解释技术研讨会地震资料的振幅相对保持剖面监控80-10°10-20°20-30°0-31°中油集团东方地球物理公司研究院处理解释技术研讨会CMP分选后最大波峰振幅地表一致性静校正后最大波峰振幅多道倾角滤波后最大波峰振幅随机噪音后最大波峰振幅叠前地震资料振幅相对保持平面、空间监控多道倾角滤波叠前时间偏移后最大波峰振幅9中油集团东方地球物理公司研究院处理解释技术研讨会精细处理后处理前处理前、后叠加地震资料的对比9中油集团东方地球物理公司研究院处理解释技术研讨会叠前地震资料提高分辨率处理监控10提高分辨率处理前及合成记录标定提高分辨率后处理及合成记录标定中油集团东方地球物理公司研究院处理解释技术研讨会全角道集扫描31°切除后的角道集叠前时间偏移CRP道集近、中、远角道集叠加角度范围扫描CRP-角度域弯曲射线模型同一中心角，角度重叠范围试验叠前部分角道集叠加处理0.000.200.400.600.801.001.200.0010.0020.0030.0040.0050.0060.0070.0080.0090.00入射角反射系数中油集团东方地球物理公司研究院处理解释技术研讨会在最大入射角内地震的部分角道集叠加0-10°10-20°20-30°中油集团东方地球物理公司研究院处理解释技术研讨会孔隙流体Kf基质Kma水孔隙流体Kf基质Kma油气岩石物理模拟技术2模拟含油气储层纵横波速度和密度中油集团东方地球物理公司研究院处理解释技术研讨会岩石物理模拟技术及研究流程泥岩体积岩石总孔隙度目的层的温度、压力、流体类型和密度含盐水储层纵横波模拟含油气储层纵横波模拟岩石孔隙的调整和岩性约束岩石含油气岩石参数的估算研究区目标层段的纵、横波速度、Vp/Vs输出岩石物理模型的选择：Xu&White(DEM)Self-consitentsDispersiveClayGrainSupportedMatrixSupportedCriticalPorosity)//1)(/1()/1(00200fdddddsKKKKKKKKKKKds/)/()/(3/4skmVskmVSKp纵横波吻合程度满意中油集团东方地球物理公司研究院处理解释技术研讨会岩石物理模拟技术中油集团东方地球物理公司研究院处理解释技术研讨会泥－砂SwФVpVsρSo流体置换模拟中油集团东方地球物理