深度域储层地震反演方法及装置制造方法及图纸

本发明专利技术涉及一种深度域储层地震反演方法及装置，包括：根据时间域叠后地震数据，获得深度域叠后地震数据；将测井曲线和深度域叠后地震数据进行标定，获得深度域测井曲线模型；对深度域叠后地震数据进行地层层位解释，并将解释结果数据网格化，获得深度域地质框架模型；分别计算单井的深度域测井曲线模型和深度域叠后地震数据的分形参数，将每一个地震道的分形参数转换为测井曲线的分形参数体；根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果。本发明专利技术基于分形算法，在避免了深度域子波的求取的同时有效地保证了地震反演结果的高分辨率。

【技术实现步骤摘要】
深度域储层地震反演方法及装置
本专利技术涉及油气勘探
，尤其涉及深度域储层地震反演方法及装置。
技术介绍
在石油天然气工业中，地震反演作为储层预测的核心技术，能够将测井纵向上的高分辨率优势和地震横向上的密集采样优势紧密结合起来，估算出地层横向上的岩性、物性等特征信息，目前已经被广泛应用于石油勘探开发的各个阶段。但是常规的时间域地震反演方法对于薄储层的识别能力较差，其测井曲线从深度域换算到时间域的过程中丢失了很多宝贵的高层信息，存在很大的局限性。近年来，对于地震反演储层预测的最为常见的做法是基于褶积模型的深度域地震资料反演。该方法的首要问题是子波的获取，即通过时深转换将时间域子波转换为深度域子波，然后根据地震道褶积和反子波，消除子波的褶积效应，从而获得一个近似的地层反射系数系列，以便用于反演高分辨率的岩层型剖面。专利技术人在实现本专利技术的过程中，发现上述现有技术存在如下不足：该方法虽然可以保留测井资料的高频信息，但是在操作过程中难以建立进行时深转换所需要的精确三维速度场和提取深度域地震子波，且在地震子波确定时，相位谱的确定容易出现误差，不适宜大规模应用。
技术实现思路
本专利技术实施例提供一种深度域储层地震反演方法，用以克服测井曲线在时深转换过程中高频信息易丢失的缺点，且避免了在反演过程中求取深度域子波，包括：根据时间域叠后地震数据，获得深度域叠后地震数据；将测井曲线和深度域叠后地震数据进行标定，获得深度域测井曲线模型；对深度域叠后地震数据进行地层层位解释，并将解释结果数据网格化，获得深度域地质框架模型；分别计算单井的深度域测井曲线模型和深度域叠后地震数据的分形参数，将每一个地震道的分形参数转换为测井曲线的分形参数体；根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果。一个实施例中，所述根据时间域叠后地震数据，获得深度域叠后地震数据，包括：建立空变速度场，根据空变速度场将时间域叠后地震数据转换为深度域叠后地震数据。一个实施例中，所述深度域测井曲线模型包括：波阻抗测井曲线模型、自然伽马测井曲线模型和孔隙度测井曲线模型其中之一或任意组合。一个实施例中，所述根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果，根据如下公式进行计算：其中，S(a)是插值点a的深度域储层地震反演结果，fline是分形随机插值的线性分量，funine是分形随机插值的随机分量，K是标定系数，H是反应深度域叠后地震数据细节和粗糙度的分维参数，σ是深度域叠后地震数据正态分布的方差，G是服从标准正态分布的Gauss随机变量，ΔX是已知点和插值点之间的距离。一个实施例中，所述分形随机插值的线性分量fline采用方位距离加权法实现，其中，在研究区域内，网格节点z(x,y)的估计值为：其中，D是以网格节点z(x,y)为原点构成的四个象限内总的等分角数，l、i、j、k是已知点的序号，n是已知点的总数，λi是在D范围内第i个已知点的权系数，zi是在D范围内第i个已知点的参数值，rj、rl分别是第j个及第l个等分角内离待估计点最近已知点的距离。一个实施例中，深度域储层地震反演方法还包括：在获得深度域储层地震反演结果后，对深度域储层地震反演结果进行精度检验。一个实施例中，所述在获得深度域储层反演结果后，对深度域储层地震反演结果进行精度检验，包括：分别抽取全部反演井的每一口井作为验证井，对未被抽取的反演井进行深度域储层地震反演，将未被抽取的反演井的反演结果与验证井的测井曲线进行对比，计算出未被抽取的反演井的反演结果和验证井的测井曲线之间的反演误差。本专利技术实施例中还提供了一种深度域储层地震反演装置，包括：时深转换模块，用于根据时间域叠后地震数据，获得深度域叠后地震数据；井震联合标定模块，用于将测井曲线和深度域叠后地震数据进行标定，获得深度域测井曲线模型；地质框架建立模块，用于对深度域叠后地震数据进行地层层位解释，并将解释结果数据网格化，获得深度域地质框架模型；分形参数计算模块，用于分别计算单井深度域测井曲线模型和深度域叠后地震数据的分形参数，将每一个地震道地震分形参数转换为测井曲线的分形参数体；分形随机插值模块，用于根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果。一个实施例中，所述时深转换模块具体用于：建立空变速度场，根据空变速度场将时间域叠后地震数据转换为深度域叠后地震数据。一个实施例中，所述深度域测井曲线模型具体包括：波阻抗测井曲线模型、自然伽马测井曲线模型和孔隙度测井曲线模型其中之一或任意组合。一个实施例中，所述分形随机插值模块，根据如下公式进行计算：其中，S(a)是插值点a的深度域储层地震反演结果，fline是分形随机插值的线性分量，funine是分形随机插值的随机分量，K是标定系数，H是反应深度域叠后地震数据细节和粗糙度的分维参数，σ是深度域叠后地震数据正态分布的方差，G是服从标准正态分布的Gauss随机变量，ΔX是已知点和插值点之间的距离。一个实施例中，所述分形随机插值的线性分量fline采用方位距离加权法实现，其中，在研究区域内，网格节点z(x,y)的估计值为：其中，D是以网格节点z(x,y)为原点构成的四个象限内总的等分角数，l、i、j、k是已知点的序号，n是已知点的总数，λi是在D范围内第i个已知点的权系数，zi是在D范围内第i个已知点的参数值，rj、rl分别是第j个及第l个等分角内离待估计点最近已知点的距离。一个实施例中，深度域储层地震反演装置还包括：结果检验模块，用于在获得深度域储地震层反演结果后，对深度域储层地震反演结果进行精度检验。一个实施例中，所述结果检验模块，具体用于：分别抽取全部反演井的每一口井作为验证井，对未被抽取的反演井进行深度域储层地震反演，将未被抽取的反演井的反演结果与验证井的测井曲线进行对比，计算出未被抽取的反演井的反演结果和验证井的测井曲线之间的反演误差。本专利技术实施例的深度域储层地震反演方法及装置，将常规时间域叠后地震数据进行时深转换得到深度域叠后地震数据；然后将测井曲线和地震数据在深度域进行标定，在地质框架约束下，采用分形随机插值算法重构深度域每一个地震道的储层参数，完成高分辨率地震反演。由于本专利技术实施例不是基于褶积模型，因此避免了深度域子波的求取，且本专利技术实施例同样克服了测井曲线时深转换过程中由于重采样而易造成高频信息丢失的缺点，有效地保证了地震反演结果的高分辨率特征。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：图1为本专利技术实施例中深度域储层地震反演方法的实施流程示意图。图2为本专利技术实施例中时间域地震剖面图。图3为本专利技术实施例中深度域地震剖面图。图4为本专利技术实施例中深度域储层地震反演装置的结构示意图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本专利技术实施例做进一步详细说本文档来自技高网...

【技术保护点】
一种深度域储层地震反演方法，其特征在在于，包括：根据时间域叠后地震数据，获得深度域叠后地震数据；将测井曲线和深度域叠后地震数据进行标定，获得深度域测井曲线模型；对深度域叠后地震数据进行地层层位解释，并将解释结果数据网格化，获得深度域地质框架模型；分别计算单井的深度域测井曲线模型和深度域叠后地震数据的分形参数，将每一个地震道的分形参数转换为测井曲线的分形参数体；根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果。

【技术特征摘要】
1.一种深度域储层地震反演方法，其特征在在于，包括：根据时间域叠后地震数据，获得深度域叠后地震数据；将测井曲线和深度域叠后地震数据进行标定，获得深度域测井曲线模型；对深度域叠后地震数据进行地层层位解释，并将解释结果数据网格化，获得深度域地质框架模型；分别计算单井的深度域测井曲线模型和深度域叠后地震数据的分形参数，将每一个地震道的分形参数转换为测井曲线的分形参数体；根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果。2.如权利要求1所述的深度域储层地震反演方法，其特征在于，所述根据时间域叠后地震数据，获得深度域叠后地震数据，包括：建立空变速度场，根据空变速度场将时间域叠后地震数据转换为深度域叠后地震数据。3.如权利要求1所述的深度域储层地震反演方法，其特征在于，所述深度域测井曲线模型包括：波阻抗测井曲线模型、自然伽马测井曲线模型和孔隙度测井曲线模型其中之一或任意组合。4.如权利要求1所述的深度域储层地震反演方法，其特征在于，所述根据深度域测井曲线模型和深度域地质框架模型，对测井曲线的分形参数体进行分形随机插值，获得深度域储层地震反演结果，根据如下公式进行计算：其中，S(a)是插值点a的深度域储层地震反演结果，fline是分形随机插值的线性分量，funine是分形随机插值的随机分量，K是标定系数，H是反应深度域叠后地震数据细节和粗糙度的分维参数，σ是深度域叠后地震数据正态分布的方差，G是服从标准正态分布的Gauss随机变量，ΔX是已知点和插值点之间的距离。5.如权利要求4所述的深度域储层地震反演方法，其特征在于，所述分形随机插值的线性分量fline采用方位距离加权法实现，其中，在研究区域内，网格节点z(x,y)的估计值为：其中，D是以网格节点z(x,y)为原点构成的四个象限内总的等分角数，l、i、j、k是已知点的序号，n是已知点的总数，λi是在D范围内第i个已知点的权系数，zi是在D范围内第i个已知点的参数值，rj、rl分别是第j个及第l个等分角内离待估计点最近已知点的距离。6.如权利要求1所述的深度域储层地震反演方法，其特征在于，还包括：在获得深度域储层地震反演结果后，对深度域储层地震反演结果进行精度检验。7.如权利要求6所述的深度域储层地震反演方法，其特征在于，所述在获得深度域储层地震反演结果后，对深度域储层地震反演结果进行精度检验，包括：分别抽取全部反演井的每一口井作为验证井，对未被抽取的反演井进行深度域储层地震反演，将未被抽取的反演井的反演结果与验证井的测井曲线进行对比，计算出未被抽取的反演井的反演结果和验证井的测井曲线之间...

【专利技术属性】
技术研发人员：张静，杨丽，韩乾凤，郝涛，沙雪梅，郝晋进，张亚军，牛慧赟，雷明，郑长龙，何巍巍，邓国鑫，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11