一种花岗岩地层划分对比方法技术

本发明专利技术提供了一种花岗岩地层划分对比方法，其包括以下步骤：1)测井曲线数据准备；2)测井曲线最大熵谱分析；3)获取数据差异曲线；4)将数据差异曲线转换为频谱曲线；5)对重构的频谱曲线进行滤波处理；6)确定异常点；7)对异常点进行编号；8)进行地层划分与对比：按序连接骨干剖面上各井关键相似差异点，即完成花岗岩地层划分与对比。本发明专利技术是对各井中的异常点进行分析比较，与传统的相似地层对比方法相比，对花岗岩这种岩性差别小、垂向变化小、旋回性不突出的地层而言，具有更强的可操作性，较常规地层对比方法更具有实用性。本发明专利技术由于采用了最大熵谱分析方法，因此地层对比定量程度加强，可操作性更强，地层对比的效率大大加强。

【技术实现步骤摘要】

本专利技术涉及一种花岗岩地层划分对比方法，属于油田勘探领域。
技术介绍
地层划分与对比是油田勘探与开发工作中必不可少的一个环节，后续的地质建模、开发层系划分以及井网部署等都是在精确的地层划分对比的基础上进行的。目前基于常规碎屑岩和碳酸盐地层已经探索出了诸如地震地层学，层序地层学，年代地层学等多种地层划分对比方法，并形成了相应的理论体系及技术手段。但是花岗岩地层划分却鲜有涉及，其主要存在两方面原因，一方面是现有的理论体系由于是基于时间范畴的，强调的是等时理念，与沉积旋回概念密切相关，而花岗岩地层无沉积旋回概念、岩性大部分较为单一，地层内幕地震反射特征通常不清楚、单井垂向电性特征大多较均匀，因此不适用于岩浆侵入及后续构造、风化作用形成的储层；另一方面是花岗岩类型油田发现较少，因此相应的研究技术手段相对匮乏，传统等时的地层划分对比思路不适用于该类储层地层划分对比工作。
技术实现思路
针对上述问题，本专利技术的目的是提供一种基于异常点分析的花岗岩地层划分对比方法。为实现上述目的，本专利技术采取以下技术方案：一种花岗岩地层划分对比方法，其包括以下步骤：1)测井曲线数据准备：将待测油田各井的测井曲线的起止位置以及拼接位置的负异常值变为0；2)测井曲线最大熵谱分析：通过最大熵谱分析函数对各井的测井曲线数据进行处理，获取各井的测井频谱曲线；3)获取数据差异曲线：计算各井的测井频谱曲线的每一个深度点的最大熵值，并根据每一个深度点对应的测井曲线的实际数据，得到各井的数据差异曲线；4)将数据差异曲线转换为频谱曲线：将各井的数据差异曲线通过复化辛普森积分公式进行数学积分变换，重构出各井的频谱曲线；5)对重构的频谱曲线进行滤波处理：根据待测油田的每口井录井及井壁取心岩性统计结果获取风化地层厚度的上限、中值和下限，并根据有效波长与地层厚度的正相关关系获取风化地层厚度的上限、中值和下限的有效波长，依据有效波长的范围选取滤波参数，从而对重构的频谱曲线进行滤波处理，得到各井的上限、中值和下限的三条滤波曲线；6)确定异常点：根据滤波处理后得到的各井的三条滤波曲线，选取三条曲线中均存在较大变化幅度的点作为差异点；7)对异常点进行编号：对步骤6)中确定的各井的异常点从上到下进行垂向编号；8)进行地层划分与对比：按序连接骨干剖面上各井关键相似差异点，即完成花岗岩地层划分与对比。在所述步骤2)中，最大熵函数SX(ω)和相关函数rx(m)之间的关系由维纳-辛钦定理确定,即 S X ( ω ) = Σ - ∞ ∞ r x ( m ) e - j ω m ]]>其中，rx(m)是采样得到的离散时间序列，j是虚数单位，ω是功率谱信号密度，m是离散变量。在所述步骤3)中，数据差异曲线的计算公式为：数据差异＝实际数据-最大熵值。在所述步骤5)中，有效波长的计算公式为：W＝k×h其中，W为有效波长，k为相关系数，h为地层厚度，k＝0.8。在所述步骤5)中，保留频谱曲线中的有效波长，将无效的低波、中高波呈梯形状滤去，从而使曲线变得更加平滑；其中0.5W有效波长～1.5W有效波长为有效波长区间，0.25W有效波长～0.5W有效波长和1.5W有效波长～1.75W有效波长为无效的低波和中高波区间。在所述步骤6)中，以各井的上限滤波曲线来确定异常点发育的深度段，上限滤波曲线向左侧的最大振幅处对应的深度点的上下5米范围内作为异常点发育的深度段，再以异常点发育的深度段所对应的中值滤波曲线的振幅最大位置确定异常点。本专利技术由于采取以上技术方案，其具有以下优点：1、本专利技术是对各井中的异常点进行分析比较，与传统的相似地层对比方法相比，对花岗岩这种岩性差别小、垂向变化小、旋回性不突出的地层而言，具有更强的可操作性，较常规地层对比方法更具有实用性。2、本专利技术由于采用了最大熵谱分析方法，因此地层对比定量程度加强，可操作性更强，地层对比的效率大大加强。3、本专利技术由于采用对比方式将数据差异曲线转换为空间域的频谱曲线，因此更能适应花岗岩这一特殊储层的特点。附图说明图1为本专利技术划分对比地层的示意图。具体实施方式下面结合附图和实施例对本专利技术进行详细的描述。本专利技术提供的一种花岗岩地层划分对比方法，其包括以下步骤：1)测井曲线数据准备通常测井人员提供的曲线已经做过相应的预处理，但由于处理方法和习惯的不同在曲线的起止位置以及拼接位置会赋以固定的常值-999或-999.25。而这种负异常会严重影响本专利技术的后续处理步骤，因此需要先将各井的测井曲线的起止位置以及拼接位置的负异常值变为0；2)测井曲线最大熵谱分析通过最大熵谱分析函数对测井曲线数据进行处理，获取各井的测井频谱曲线。其中，最大熵函数SX(ω)和相关函数rx(m)之间的关系由维纳-辛钦定理确定,即 S X ( ω ) = Σ - ∞ ∞ r x ( m ) e - j ω m ]]>其中，rx(m)是采样得到的离散时间序列，j是虚数单位，ω是功率谱信号密度，m是离散变量；本实施例中，将步骤1)中预处理后的测井数据以las格式输入到商业软件cyclolog软件中，预先给定目的层段顶底深度和窗口宽度，然后调用最大熵谱分析函数对测井数据进行处理，即可得到各井的测井频谱曲线。3)获取数据差异曲线计算各井的测井频谱曲线每一个深度点的最大熵值，并根据每一个深度点对应的测井曲线的实际数据，得到各井的数据差异曲线；其中，数据差异＝实际数据-最大熵值。4)将数据差异曲线转换为频谱曲线将步骤3)中的各井的数据差异曲线通过复化辛普森积分公式进行数学积分变换，重构出一条新的各井的频谱曲线。其中，该过程可借用商业软件cyclolog软件的Inpefa模块来实现，具体步骤如下：在数据子窗口上点击鼠标右键，选择Processing→Spectral Attributes→INPEFA；在预测误差滤波窗口选择分析的目的层段(默认值是整口井)，并给出预测窗口的大小(默认值是10米)，即可完成数据差异曲线到频谱曲线的转换。5)对重构的频谱曲线进行滤波处理根据待测油田的每口井录井及井壁取心岩性统计结果获取风化地层厚度的上限、中值和下限，并根据有效波长与地层厚度的正相关关系获取风化地层厚度的上限、中值和下限的有效波长，即：W有效波长＝k相关系数×h地层本文档来自技高网...

【技术保护点】
一种花岗岩地层划分对比方法，其包括以下步骤：1)测井曲线数据准备：将待测油田各井的测井曲线的起止位置以及拼接位置的负异常值变为0；2)测井曲线最大熵谱分析：通过最大熵谱分析函数对各井的测井曲线数据进行处理，获取各井的测井频谱曲线；3)获取数据差异曲线：计算各井的测井频谱曲线的每一个深度点的最大熵值，并根据每一个深度点对应的测井曲线的实际数据，得到各井的数据差异曲线；4)将数据差异曲线转换为频谱曲线：将各井的数据差异曲线通过复化辛普森积分公式进行数学积分变换，重构出各井的频谱曲线；5)对重构的频谱曲线进行滤波处理：根据待测油田的每口井录井及井壁取心岩性统计结果获取风化地层厚度的上限、中值和下限，并根据有效波长与地层厚度的正相关关系获取风化地层厚度的上限、中值和下限的有效波长，依据有效波长的范围选取滤波参数，从而对重构的频谱曲线进行滤波处理，得到各井的上限、中值和下限的三条滤波曲线；6)确定异常点：根据滤波处理后得到的各井的三条滤波曲线，选取三条曲线中均存在较大变化幅度的点作为差异点；7)对异常点进行编号：对步骤6)中确定的各井的异常点从上到下进行垂向编号；8)进行地层划分与对比：按序连接骨干剖面上各井关键相似差异点，即完成花岗岩地层划分与对比。...

【技术特征摘要】
1.一种花岗岩地层划分对比方法，其包括以下步骤：1)测井曲线数据准备：将待测油田各井的测井曲线的起止位置以及拼接位置的负异常值变为0；2)测井曲线最大熵谱分析：通过最大熵谱分析函数对各井的测井曲线数据进行处理，获取各井的测井频谱曲线；3)获取数据差异曲线：计算各井的测井频谱曲线的每一个深度点的最大熵值，并根据每一个深度点对应的测井曲线的实际数据，得到各井的数据差异曲线；4)将数据差异曲线转换为频谱曲线：将各井的数据差异曲线通过复化辛普森积分公式进行数学积分变换，重构出各井的频谱曲线；5)对重构的频谱曲线进行滤波处理：根据待测油田的每口井录井及井壁取心岩性统计结果获取风化地层厚度的上限、中值和下限，并根据有效波长与地层厚度的正相关关系获取风化地层厚度的上限、中值和下限的有效波长，依据有效波长的范围选取滤波参数，从而对重构的频谱曲线进行滤波处理，得到各井的上限、中值和下限的三条滤波曲线；6)确定异常点：根据滤波处理后得到的各井的三条滤波曲线，选取三条曲线中均存在较大变化幅度的点作为差异点；7)对异常点进行编号：对步骤6)中确定的各井的异常点从上到下进行垂向编号；8)进行地层划分与对比：按序连接骨干剖面上各井关键相似差异点，即完成花岗岩地层划分与对比。2.如权利要求1所述的一种花岗岩地层划分对比方法，其特征在于：在所述步骤2)中，最大熵函数SX(ω)和相关函数rx(m)之间的关系由维纳-辛钦定理确定,即 S X ( ω ) = ...

【专利技术属性】
技术研发人员：范廷恩，梁旭，宋来明，胡光义，蔡文涛，王晖，周建楠，汤婧，张雨晴，潘晓庆，
申请(专利权)人：中国海洋石油总公司，中海油研究总院，
类型：发明
国别省市：北京;11