基于入射角的avo近似公式进行属性提取的方法
【技术领域】
[0001] 本发明设及地震资料处理领域，是一种利用地震数据进行储层识别的方法。
【背景技术】
[0002] AV0技术是根据振幅随炮检距的变化规律所反映出的地下岩性及其孔隙流体的性 质来直接预测油气和估计地壳岩性参数的一项技术。其理论基础是描述平面波在水平分界 面上反射和折射的Zo巧pritz方程。完整的Zo巧pritz方程形式复杂，物理意义不明确， 无法直接应用于地震数据的处理和参数反演。尽管人们早就认识到了该一点，但在早期的 研究中，人们仅从数值计算的角度对Zoeppritz描述的反射系数进行了各种讨论和近似。 Muskat系统地计算了弹性介质中平面波的反射和透射系数；Koe化ed利用17组不同弹性参 数的模型，详细地研究了泊松比对两个各向同性介质之间的反射界面所产生的产生的反射 系数的影响，并第一次给出了将泊松比与反射系数直接联系起来的Zoeppritz近似方程。 Bcxrtfeld第一个给出了物理意义比较明确的反射系数近似计算公式，使得基于Zo巧pritz 方程的反射和透射问题真正应用到了勘探地球物理领域。在W上工作的基础上，不少人做 了大量的工作，分别从不同的方面对Zoeppritz方程进行了详细的讨论和简化，提出了不 同假设条件下的纵波反射振幅的近似表达式。但是所有的近似公式中角度都是应用的上层 入射角和下层透射角的平均，给进一步的参数估算带来一定的误差。
【发明内容】
[0003] 本发明的目的在于避免现有技术的不足而提供基于入射角的AV0近似公式进行 属性提取的方法。
[0004] 本发明的技术方案是：基于入射角的AV0近似公式进行属性提取的 方法，实现过程包括；1)基于入射角的AV0近似公式的提出；2)基于入射角 的AV0近似公式的理论推导；3)应用常规方法提取P、G，Rs的属性；4)将属性归 一化；5)由井数据或者研究区域的经验公式确定系数b值，然后应用公式（14)
对常规方法提取的属性进行修正；6)将 修正结果与井数据进行对比分析，根据分析结果对修正结果进行校正，得到最终的属性；其 中，P是由零炮检距构成的地震道，即纵波叠加道，代表对反射界面两侧的波阻抗变化的响 应；斜率G为梯度叠加道，代表对横波速度、纵波速度和体密度变化的响应，是由入射角近 似公式计算得到的，GO代表对横波速度、纵波速度和体密度变化的响应，由常规近似公式计 算的得到；Rs表示由入射角近似公式计算的横波阻抗反射系数，馬》为由常规近似公式计算 的横波阻抗反射系数； 具体的求取方法为： 1)由Aki和化U巧近似可得；
R( 0 )表示纵波反射系数，妍和巧分别为反射界面上、下地层介质的密度，和分 别是反射界面上下上、下地层的纵波速度，分别是上、下地层的横波速度，巧和巧分 别是地层的入射角和透射角 2)根据Snell定律有；
对公式（1)左边项进行泰勒展开，及
将公式（3)代入公式（2)可得：
令
4)Garner等人于1974提出的密度与纵波速度之间的经验公式的形式为[6];
其中，a和b为系数（一般取a=0. 31，b=0. 25) 由上式可得
即：
5) -般由化uey近似计算截距和梯度属性P和G(进一步可计算横波阻抗反射系数Rs 等属性）的公式为：
其中，巧是根据平均角计算出的梯度 6) 在实际应用公式（11)计算属性时，由于求取角度A不方便，常常用入射角S代替6进 行估算，即应用下式计算P和Rs:
比较式（11)和式（13)可得：
所W，计算横波阻抗反射系数Rs的公式为：
该里的P和Rs。由公式（13 )计算得到。
[0005] 与现有技术相比，本发明的有益效果是：计算方便、精确，与实际的参数误差小。
【附图说明】
[0006] 图 1 为Ostrander模型。
[0007] 图 2 为Goodway模型。
[0008] 图3为MarmousiII模型中的气层。
[0009] 图4为修改前后误差对比分析。
[0010] 图5为从mar2模型中抽取的一道数据。
[0011] 图6为MarmousiII模型的属性对比。
[0012] 图7MarmousiII模型横波合成地震记录的比较。其中（a)为垂直入射计算的 Rs与子波權积得到的合成地震记录。化）为常规方法计算的Rs与子波權积得到的合成地 震记录。（C)为修正之后的Rs与子波權积得到的合成地震记录。
[0013] 图8为实施例的叠前角度道集。其中（a)为某地区8°角度叠加道集。（b)为某 地区16°角度叠加道集。（C)为某地区24°角度叠加道集。
[0014] 图9为实施例基于修正公式计算的Rs属性。
【具体实施方式】
[0015] 基于入射角的AV0近似公式进行属性提取的方法，实现过程包括；1)基于入射角 的AV0近似公式的提出；2)基于入射角的AV0近似公式的理论推导；3)应用常规方法提取 P、G，Rs的属性；4)将属性归一化；5)由井数据或者研究区域的经验公式确定系数b值，然 后应用公式（14)
对常规方法提取的属性 进行修正；6)将修正结果与井数据进行对比分析，根据分析结果对修正结果进行校正，得 到最终的属性；其中，P是由零炮检距构成的地震道，即纵波叠加道，代表对反射界面两侧 的波阻抗变化的响应；斜率G为梯度叠加道，代表对横波速度、纵波速度和体密度变化的响 应，是由入射角近似公式计算得到的，GO代表对横波速度、纵波速度和体密度变化的响应， 由常规近似公式计算的得到；Rs表示由入射角近似公式计算的横波阻抗反射系数，馬C为由 常规近似公式计算的横波阻抗反射系数； 具体的求取方法为： 1) 由Aki和化U巧近似可得；
R( 0 )表示纵波反射系数，A和A分别为反射界面上、下地层介质的密度，?和分 别是反射界面上下上、下地层的纵波速度Ki和吃分别是上、下地层的横波速度，巧和巧分 别是地层的入射角和透射角 2) 根据Snell定律有；
对公式（1)左边项进行泰勒展开，及
3)该样可W得到基于入射角的AVO近似公式：
4) Ga