rner等人于1974提出的密度与纵波速度之间的经验公式的形式为[6];
其中，a和b为系数（一般取a=0. 31，b=0. 25) 由上式可得
5) -般由化uey近似计算截距和梯度属性P和G(进一步可计算横波阻抗反射系数Rs 等属性）的公式为：
其中，巧是根据平均角计算出的梯度 6)在实际应用公式（11)计算属性时，由于求取角度A不方便，常常用入射角巧代替0进 行估算，即应用下式计算P和Rs:
比较式（11)和式（13)可得：
所W，计算横波阻抗反射系数Rs的公式为：
该里的P和Rs。由公式（13 )计算得到。
[0016] 按照上面的方法原理，首先应用模型进行了分析，下面给出具体分析结果。
[0017] a、模型简介 本技术应用了立个模型，包括化trander所构造的含气砂岩模型w、G〇〇dway根据实测 资料给出的含气砂岩模型和含气砂岩的四种典型AV0模型W及MarmousiII中的气层 模型W。其参数分别见表1 - 3。
[0018] 表1 化trander(1984)的S层含油气砂岩与页岩模型
b、修巧刖后F-P汲化巧《甄团线卿此分机 图1-3中，实线为精确的zoeppritz方程计算结果，点线是将精确的shuey近似表示成 随入射角变换的函数得到的结果，"+ "线是应用修改之后的公式计算的结果。
[001引从图1-3中，可W看出当上、下储层差值比较大时，如果直接应用入射角代替平 均角存在比较大的误差。该样就为下一步的角道集提取，W及属性的计算带来很大的误 差。而修改之后的公式在角度比较小时（（伴<沒《3<r))，结果是非常准确的，与精确的 zoeppritz方程计算结果之间差别非常的小。修改前后，误差定量分析如图所示。
[0020] 在图4中，"1"代表Ostrander模型；"2"代表Goodway模型；"3"MarmousiII中 的气层模型。黑色为修正之前的相对误差绝对值，灰色为修正之后的相对误差绝对值。
[0021] 从图4中可W看出，经过修正之后误差大大降低。特别是气层引起上、下储层物性 差别比较大的时候。如图中的"1"和"3"所示。
[0022] C、MarmousiII模型分析 从MarmousiII模型中抽取一道进行方法验证，所取的道集中包含了气、油层W及水 层。
[0023] 图6中黑线是用阻抗差计算的结果，星号是利用常规的方法计算Rs属性，菱形是 用修正公式计算的结果。从图中可W看到在气层出现的地方，常规方法提取的Rs属性与用 阻抗差公式计算得到结果是完全相反的。而应用修正公式修正后巧日图中，菱形所示)，结果 能够比较好的与精确值吻合，特别是在气层出现的位置。图7是应用各种方法计算Rs然后 与40化雷克子波權积得到的合成地震记录。图7 (a)是由模型的横波速度和密度直接计算 Rs，然后再与40化雷克子波權积得到的合成地震记录。图7(b)是应用常规方法计算的Rs 与40化雷克子波權积得到的合成地震记录。图7(C)是应用修正公式计算的Rs与40化雷 克子波權积得到的合成地震记录。从图7的对比中可W看出，常规方法计算的Rs在气层位 置会产生很大的误差，出现了数值反转巧日图7 (b)所示)，经过修正公式修正之后气层处的 误差消失，并且数值与精确计算结果一致，如图7(C)所示。
【主权项】
1.基于入射角的AVO近似公式进行属性提取的方法，其特征在于：实现 过程包括：1)基于入射角的AVO近似公式的提出；2)基于入射角的AVO近 似公式的理论推导；3)应用常规方法提取P、G，Rs的属性；4)将属性归一 化；5)由井数据或者研宄区域的经验公式确定系数b值，然后应用公式（14)对常规方法提取的属性进行修正；6)将 修正结果与井数据进行对比分析，根据分析结果对修正结果进行校正，得到最终的属性；其 中，P是由零炮检距构成的地震道，即纵波叠加道，代表对反射界面两侧的波阻抗变化的响 应；斜率G为梯度叠加道，代表对横波速度、纵波速度和体密度变化的响应，是由入射角近 似公式计算得到的，GO代表对横波速度、纵波速度和体密度变化的响应，由常规近似公式计 算的得到；Rs表示由入射角近似公式计算的横波阻抗反射系数，4。为由常规近似公式计算 的横波阻抗反射系数； 具体的求取方法为： 1) 由Aki和Shuey近似可得：R ( Θ )表示纵波反射系数，A和A分别为反射界面上、下地层介质的密度，&和&3分 别是反射界面上下上、下地层的纵波速度，&和^分别是上、下地层的横波速度，幾和咚分 别是地层的入射角和透射角 2) 根据Snell定律有：对公式（1)左边项进行泰勒展开，及又因为：4. Gardner等人于1974提出的密度与纵波速度之间的经验公式的形式为[6]:其中，a和b为系数（一般取a=0. 31，b=0. 25) 由上式可得则5) -般由Shuey近似计算截距和梯度属性P和G (进一步可计算横波阻抗反射系数Rs 等属性）的公式为：其中，G是根据平均角计算出的梯度 6) 在实际应用公式（11)计算属性时，由于求取角度爲不方便，常常用入射角4代替进 行估算，即应用下式计算P和Rs:比较式（11)和式（13)可得：所以，计算横波阻抗反射系数Rs的公式为：这里的P和Rsci由公式（13 )计算得到。
【专利摘要】基于入射角的AVO近似公式进行属性提取的方法，其特征在于：实现过程包括：1）基于入射角的AVO近似公式的提出；2）基于入射角的AVO近似公式的理论推导；3）应用常规方法提取P、G，Rs的属性；4）将属性归一化；5）由井数据或者研究区域的经验公式确定系数b值，然后应用公式（14） 和公式（15）对常规方法提取的属性进行修正；6）将修正结果与井数据进行对比分析，根据分析结果对修正结果进行校正，得到最终的属性,其优点是：计算方便、精确，与实际的参数误差小。
【IPC分类】G01V1/28
【公开号】CN104977610
【申请号】CN201510348231
【发明人】张广智, 郑静静, 印兴耀
【申请人】中国石油大学(华东)
【公开日】2015年10月14日
【申请日】2015年6月23日