薄储层的识别方法及装置制造方法及图纸

本发明专利技术公开了一种薄储层的识别方法及装置，属于油气物探工程领域。所述方法包括：获取待识别储层的原始地震道信号S；分别计算所述原始地震道信号S的第2阶至第n阶导数，得到n个导数：S

【技术实现步骤摘要】
薄储层的识别方法及装置
本专利技术涉及油气物探工程领域，特别涉及一种薄储层的识别方法及装置。
技术介绍
随着油气勘探程度的提高，油气勘探的目标更加隐蔽，油气储层的厚度更薄，特别是对于厚度小于10米的薄储层，例如砂泥岩薄储层，激发地震波后，从该薄储层反射的地震道信号较为微弱，根据该微弱的反射地震道信号识别薄储层的难度较高。相关技术中，一般采用谱白化算法、谱蓝化算法以及反Q滤波算法对待识别储层反射的地震道信号进行检测和处理，通过上述算法能够展宽信号的频谱，提高地震道信号的分辨率，进而能够更准确的识别出该待识别储层中的薄储层。但是，相关技术中的薄储层的识别算法计算时间较长，计算复杂度较高，对薄储层的识别效率低。
技术实现思路
为了解决相关技术中对薄储层的识别效率较低的问题，本专利技术实施例提供了一种薄储层的识别方法及装置。所述技术方案如下：一方面，提供了一种薄储层的识别方法，所述方法包括：获取待识别储层的原始地震道信号S；分别计算所述原始地震道信号S的第2阶至第n阶导数，得到n-1个导数：S2，...Si，...，Sn，其中，Si表示所述原始地震道信号S的第i阶导数，所述n为根据所述原始地震道信号S的导数确定的叠加阈值，所述i为大于等于2，且小于等于所述n的整数；根据所述原始地震道信号S和所述n-1个导数：S2，...Si，...，Sn，通过总体经验模态分解CEEMD算法确定目标地震道信号S’，所述CEEMD算法为：其中，Gau()表示对括号内的内容进行高斯变换，Nor()表示对括号内的内容进行归一化处理，M为预设的第一阈值，N为预设的第二阈值，所述第一阈值M为满足条件：4M-2≤n中的最大值，所述第二阈值其中表示向下取整；根据所述目标地震道信号S’，识别所述待识别储层中的薄储层。可选的，所述分别计算所述原始地震道信号S的第2阶至第n阶导数，得到n-1个导数：S2，...Si，...，Sn，包括：当j≥2时，判断待分析导数对应的主频Dj是否满足循环终止条件，所述待分析导数包括：所述原始地震道信号S的第j阶导数Sj，第j+1阶导数Sj+1和第j-1阶导数Sj-1，所述循环终止条件为：其中，ε为预设的终止参数，Dj，Dj-1和Dj+1分别为根据主频计算公式确定的所述原始地震道信号S的第j阶导数Sj的主频，所述原始地震道信号S的第j+1阶导数Sj+1的主频，以及所述原始地震道信号S的第j-1阶导数Sj-1的主频，所述主频计算公式为：其中，f为频率变量，Sj(f)为所述原始地震道信号S的第j阶导数Sj的频谱函数；当所述述待分析导数对应的主频Dj满足循环终止条件时，将j+2得到新的j，重复执行上述判断过程，直至所述待分析导数对应的主频Dj不满足所述循环终止条件，将不满足所述循环终止条件的待分析导数对应的j确定为所述叠加阈值。可选的，所述预设的终止参数ε的取值范围为0.2至0.3。可选的，所述Nor()表示对括号内的内容进行归一化处理后得到的数值在区域[-1，1]内。另一方面，提供了一种薄储层的识别装置，所述装置包括：获取模块，用于获取待识别储层的原始地震道信号S；计算模块，用于分别计算所述原始地震道信号S的第2阶至第n阶导数，得到n-1个导数：S2，...Si，...，Sn，其中，Si表示所述原始地震道信号S的第i阶导数，所述n为根据所述原始地震道信号S的导数确定的叠加阈值，所述i为大于等于2，且小于等于所述n的整数；确定模块，用于根据所述原始地震道信号S和所述n-1个导数：S2，...Si，...，Sn，通过CEEMD算法确定目标地震道信号S’，所述CEEMD算法为：其中，Gau()表示对括号内的内容进行高斯变换，Nor()表示对括号内的内容进行归一化处理，M为预设的第一阈值，N为预设的第二阈值，所述第一阈值M为满足条件：4M-2≤n中的最大值，所述第二阈值其中表示向下取整；识别模块，用于根据所述目标地震道信号S’，识别所述待识别储层中的薄储层。可选的，所述计算模块，还用于：当j≥2时，判断待分析导数对应的主频Dj是否满足循环终止条件，所述待分析导数包括：所述原始地震道信号S的第j阶导数Sj，第j+1阶导数Sj+1和第j-1阶导数Sj-1，所述循环终止条件为：其中，ε为预设的终止参数，Dj，Dj-1和Dj+1分别为根据主频计算公式确定的所述原始地震道信号S的第j阶导数Sj的主频，所述原始地震道信号S的第j+1阶导数Sj+1的主频，以及所述原始地震道信号S的第j-1阶导数Sj-1的主频，所述主频计算公式为：其中，f为频率变量，Sj(f)为所述原始地震道信号S的第j阶导数Sj的频谱函数；当所述述待分析导数对应的主频Dj满足循环终止条件时，将j+2得到新的j，重复执行上述判断过程，直至所述待分析导数对应的主频Dj不满足所述循环终止条件，将不满足所述循环终止条件的待分析导数对应的j确定为所述叠加阈值。可选的，所述预设的终止参数ε的取值范围为0.2至0.3。可选的，所述Nor()表示对括号内的内容进行归一化处理后得到的数值在区域[-1，1]内。本专利技术实施例提供的技术方案带来的有益效果是：本专利技术实施例提供了一种薄储层的识别方法及装置，该方法包括：获取待识别储层的原始地震道信号S；分别计算所述原始地震道信号S的第2阶至第n阶导数，得到n-1个导数：S2，...Si，...，Sn；根据所述原始地震道信号S和所述n-1个导数：S2，...Si，...，Sn，最后通过CEEMD算法确定目标地震道信号S’，并根据所述目标地震道信号S’，识别所述待识别储层中的薄储层。该薄储层的识别算法复杂度较低，计算用时较短，提高了薄储层的识别效率。附图说明为了更清楚地说明本专利技术实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是本专利技术实施例提供的一种薄储层的识别方法的流程图；图2是本专利技术实施例提供的另一种薄储层的识别方法的流程图；图3是本专利技术实施例提供的一种薄储层的识别装置的结构示意图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚，下面将结合附图对本专利技术实施方式作进一步地详细描述。本专利技术实施例提供了一种薄储层的识别方法，参见图1，该方法包括：步骤101、获取待识别储层的原始地震道信号S。步骤102、分别计算该原始地震道信号S的第2阶至第n阶导数，得到n-1个导数：S2，...Si，...，Sn，其中，Si表示该原始地震道信号S的第i阶导数，该n为根据该原始地震道信号S的导数确定的叠加阈值，该i为大于等于2，且小于等于该n的整数。步骤103、根据该原始地震道信号S和该n-1个导数：S2，...Si，...，Sn，通过总体经验模态分解(英文：ComplementaryEnsembleEmpiricalModeDecomposition；简称：CEEMD)算法确定目标地震道信号S’。该CEEMD算法为：其中，Gau()表示对括号内的内容进行高斯变换，Nor()表示对括号内的内容进行归一化处理，M为预设的第一阈值，N为预设的第二阈值，该第一阈值M为满足条件：4M-2≤n中的最大值，该第本文档来自技高网...

【技术保护点】
一种薄储层的识别方法，其特征在于，所述方法包括：获取待识别储层的原始地震道信号S；分别计算所述原始地震道信号S的第2阶至第n阶导数，得到n‑1个导数：S

【技术特征摘要】
1.一种薄储层的识别方法，其特征在于，所述方法包括：获取待识别储层的原始地震道信号S；分别计算所述原始地震道信号S的第2阶至第n阶导数，得到n-1个导数：S2，...Si，...，Sn，其中，Si表示所述原始地震道信号S的第i阶导数，所述n为根据所述原始地震道信号S的导数确定的叠加阈值，所述i为大于等于2，且小于等于所述n的整数；根据所述原始地震道信号S和所述n-1个导数：S2，...Si，...，Sn，通过总体经验模态分解CEEMD算法确定目标地震道信号S’，所述CEEMD算法为：其中，Gau()表示对括号内的内容进行高斯变换，Nor()表示对括号内的内容进行归一化处理，M为预设的第一阈值，N为预设的第二阈值，所述第一阈值M为满足条件：4M-2≤n中的最大值，所述第二阈值其中表示向下取整；根据所述目标地震道信号S’，识别所述待识别储层中的薄储层。2.根据权利要求1所述的方法，其特征在于，所述分别计算所述原始地震道信号S的第2阶至第n阶导数，得到n-1个导数：S2，...Si，...，Sn，包括：当j≥2时，判断待分析导数对应的主频Dj是否满足循环终止条件，所述待分析导数包括：所述原始地震道信号S的第j阶导数Sj，第j+1阶导数Sj+1和第j-1阶导数Sj-1，所述循环终止条件为：其中，ε为预设的终止参数，Dj，Dj-1和Dj+1分别为根据主频计算公式确定的所述原始地震道信号S的第j阶导数Sj的主频，所述原始地震道信号S的第j+1阶导数Sj+1的主频，以及所述原始地震道信号S的第j-1阶导数Sj-1的主频，所述主频计算公式为：其中，f为频率变量，Sj(f)为所述原始地震道信号S的第j阶导数Sj的频谱函数；当所述述待分析导数对应的主频Dj满足循环终止条件时，将j+2得到新的j，重复执行上述判断过程，直至所述待分析导数对应的主频Dj不满足所述循环终止条件，将不满足所述循环终止条件的待分析导数对应的j确定为所述叠加阈值。3.根据权利要求2所述的方法，其特征在于，所述预设的终止参数ε的取值范围为0.2至0.3。4.根据权利要求1至3任一所述的方法，其特征在于，所述Nor()表示对括号内的内容进行归一化处理后得到的数值在区域...

【专利技术属性】
技术研发人员：刘喜恒，易远元，秦凤启，陈伟，王亚，耿海燕，王旭峰，冯小英，叶秋焱，祖志勇，叶玮，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11