基于叠前拟横波反射率属性的混合阻抗砂体识别方法技术

本发明专利技术公开了一种基于叠前拟横波反射率属性的混合阻抗砂体识别方法，属于石油勘探领域。所述混合阻抗砂体识别方法包括：确定砂体的阻抗类型；基于砂体的阻抗类型，通过井震标定确定砂体的地震响应特征；基于砂体的地震响应特征，应用叠前地震角道集数据进行AVO属性反演求取P

【技术实现步骤摘要】
基于叠前拟横波反射率属性的混合阻抗砂体识别方法


[0001]本专利技术属于石油勘探领域，涉及石油勘探中地震资料解释领域的砂体识别刻画技术，特别涉及涉及一种基于叠前拟横波反射率属性的混合阻抗砂体识别方法。

技术介绍

[0002]三维地震识别砂体展布特征以砂体与围岩的阻抗差异为基础，砂体较围岩阻抗值低时，砂体表现为“顶谷底峰亮点反射”，砂体较围岩阻抗值高时，砂体为“顶峰底谷亮点反射”，基于砂体此种地震响应特征通过均方根振幅地层切片提取完成砂体平面刻画。
[0003]现有的使用叠后地震进行属性分析砂体刻画技术存在以下的问题和缺点：
[0004]1.常规叠后地震数据仅能刻画与围岩阻抗差异较大的高阻抗及低阻抗砂体，对于与围岩阻抗差异较小，地震反射亮点不突出的砂体识别能力较弱。
[0005]2.受埋藏深度、物源方向、成岩演化阶段等因素影响同一期次砂体可能表现为：高纵波阻抗、中纵波阻抗、低纵波阻抗特征，常规地震属性不能有效刻画砂体展布。
[0006]3.常规地层切片为等比例切片，相对等时，当沉积厚度发生变化时等比例切片等时性较差，砂体刻画准确率降低。

技术实现思路

[0007]本专利技术的目的在于克服现有技术的不足，提供一种基于叠前拟横波反射率属性的混合阻抗砂体识别方法。
[0008]本专利技术的目的是通过以下技术方案来实现的：基于叠前拟横波反射率属性的混合阻抗砂体识别方法，包括：
[0009]确定砂体的阻抗类型；
[0010]基于砂体的阻抗类型，通过井震标定确定砂体的地震响应特征；
[0011]基于砂体的地震响应特征，应用叠前地震角道集数据进行AVO属性反演，求取P
‑
G伪横波反射率属性；
[0012]分别对P
‑
G伪横波反射率属性和原始地震数据求取均方根振幅，并对P
‑
G伪横波反射率属性的均方根振幅和原始地震数据的均方根振幅进行振幅归一化；
[0013]对振归一化后的P
‑
G伪横波反射率属性的均方根振幅和原始地震数据的均方根振幅进行振幅相加计算，得到融合属性体；
[0014]采用基于全局自动体追踪的高密度等时切片提取融合属性体的振幅值，并根据融合属性体的振幅值刻画砂体的平面展布。
[0015]进一步地，所述P
‑
G伪横波反射率属性的求取方法为：
[0016]基于叠前地震角道集数，应用Aki
‑
Richards方程，得到AVO属性；
[0017]根据AVO属性得到截距和梯度，并将截距和梯度做插值得到P
‑
G伪横波反射率属性。
[0018]进一步地，所述原始地震数据的均方根振幅的求取方法为：
[0019]在一个给定的时窗内，计算振原始地震数据的幅平方和的平均值；
[0020]将原始地震数据的幅平方和的平均值开平方得到原始地震数据的均方根振幅。
[0021]进一步地，所述P
‑
G伪横波反射率属性的均方根振幅的求取方法为：
[0022]在一个给定的时窗内，计算P
‑
G伪横波反射率属性的振幅平方和的平均值；
[0023]将P
‑
G伪横波反射率属性的幅平方和的平均值开平方得到P
‑
G伪横波反射率属性的均方根振幅。
[0024]进一步地，对P
‑
G伪横波反射率属性的均方根振幅和原始地震数据的均方根振幅进行振幅归一化，包括：
[0025]分别将P
‑
G伪横波反射率属性的均方根振幅和原始地震数据的均方根振幅的值域范围通过数学算法规整为同一值域空间[1
‑
10000]。
[0026]进一步地，所述高密度等时切片的提取方法为：
[0027]基于全局自动体追踪建立相对等时框架模型；
[0028]基于相对等时框架模型，按照所需的目的层顶底窗口提取高密度的等时解释层位切片，得到高密度等时切片。
[0029]本专利技术的有益效果是：
[0030](1)本专利技术应用P
‑
G伪横波反射率属性弥补了使用常规地震属性进行中高阻抗砂体刻画精度较低的问题，将中高阻抗砂体刻画方法从叠后引入叠前；
[0031](2)本专利技术将伪横波反射率属性及常规地震属性进行归一化融合，可实现在单一属性上同时进行高、中、低纵波阻抗砂体识别；
[0032](3)本专利技术应用全局优化自动追踪结果建立等时地层切片较等比例切片提高了地层切片的等时性，进而提升了砂体识别的精度。
附图说明
[0033]图1为本专利技术中混合阻抗砂体识别方法的一种实施例的流程图；
[0034]图2为一个实施例中的合成地震记录剖面；
[0035]图3为一个实施例中的原始地震剖面；
[0036]图4为一个实施例中的伪横波反射率剖面；
[0037]图5为一个实施例中的融合后的均方根属性剖面；
[0038]图6为一个实施例中的原始地震均方根振幅等时地层切片；
[0039]图7为一个实施例中的伪横波反射率均方根振幅等时地层切片。
具体实施方式
[0040]下面将结合实施例，对本专利技术的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0041]参阅图1至图7，本专利技术提供一种基于叠前拟横波反射率属性的混合阻抗砂体识别方法：
[0042]如图1所示，基于叠前拟横波反射率属性的混合阻抗砂体识别方法，包括：
[0043]步骤S100.确定砂体的阻抗类型。
[0044]例如，由于钻井后砂体都有测井曲线，因此可以基于预设的分类标准，根据砂体的阻抗特征分类，得到高、中、低阻抗分类。
[0045]步骤S200.基于砂体的阻抗类型，通过井震标定确定砂体的地震响应特征。
[0046]所述井震标定为通过子波和地震反射系数褶积求取模拟的地震道过程。震褶积模型数学形式为：
[0047]X(t)＝W(t)*R(t)+N(t)
[0048]式中，X(t)表示地震记录；W(t)为地震子波；R(t)为地层反射系数；N(t)为随机环境噪音。
[0049]具体的，通过井震标定得到时深表，完成深度域井与时间域地震匹配，其中经标定后河道砂体所对应地震反射的振幅、相位和频率特征即为砂体的地震响应特征。
[0050]步骤S300.基于砂体的地震响应特征，应用叠前地震角道集数据进行AVO属性反演，求取P
‑
G伪横波反射率属性。
[0051]在一个实施例中，所述P
‑
G伪横波反射率属性的求取方法为：基于叠前地震角道集数，应用Aki
‑
Richards本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于叠前拟横波反射率属性的混合阻抗砂体识别方法，其特征在于，包括：确定砂体的阻抗类型；基于砂体的阻抗类型，通过井震标定确定砂体的地震响应特征；基于砂体的地震响应特征，应用叠前地震角道集数据进行AVO属性反演，求取P
‑
G伪横波反射率属性；分别对P
‑
G伪横波反射率属性和原始地震数据求取均方根振幅，并对P
‑
G伪横波反射率属性的均方根振幅和原始地震数据的均方根振幅进行振幅归一化；对振归一化后的P
‑
G伪横波反射率属性的均方根振幅和原始地震数据的均方根振幅进行振幅相加计算，得到融合属性体；采用基于全局自动体追踪的高密度等时切片提取融合属性体的振幅值，并根据融合属性体的振幅值刻画砂体的平面展布。2.根据权利要求1所述的基于叠前拟横波反射率属性的混合阻抗砂体识别方法，其特征在于，所述P
‑
G伪横波反射率属性的求取方法为：基于叠前地震角道集数，应用Aki
‑
Richards方程，得到AVO属性；根据AVO属性得到截距和梯度，并将截距和梯度做插值得到P
‑
G伪横波反射率属性。3.根据权利要求1所述的基于叠前拟横波反射率属性的混合阻抗砂体识别方法，其特征在于，所述原始地震数据的均方根振幅的求取...

【专利技术属性】
技术研发人员：李建海，邱文霜，
申请(专利权)人：福瑞升成都科技有限公司，
类型：发明
国别省市：