CN104832161B - 基于双尺度相关对比的自动深度校正方法及系统 - Google Patents
基于双尺度相关对比的自动深度校正方法及系统 Download PDFInfo
- Publication number
- CN104832161B CN104832161B CN201510188038.0A CN201510188038A CN104832161B CN 104832161 B CN104832161 B CN 104832161B CN 201510188038 A CN201510188038 A CN 201510188038A CN 104832161 B CN104832161 B CN 104832161B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- curve
- logging
- mover
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/04—Measuring depth or liquid level
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Measurement Of Radiation (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
一种基于双尺度相关对比的自动深度校正方法，包括：S1、加载电成像测井资料；并对电成像测井资料进行预处理；S2、按照常规测井资料中的常规测井自然伽马曲线的采样间隔，对成像测井自然伽马曲线进行欠采样插值；S3、在第一尺度窗口下对成像测井自然伽马的校正曲线与常规测井曲线通过第一次对比进行粗校正；S4、按照成像测井自然伽马曲线的采样间隔，对常规测井自然伽马曲线重采样；S5、在第二尺度窗口中对成像测井自然伽马曲线与常规测井自然伽马曲线在第一尺度窗口下得到粗校正的深度范围内通过第二次对比进行精细校正；S6、在得到双尺度下成像测井自然伽马的校正曲线后，其他成像测井曲线按照成像测井自然伽马的校正曲线规律进行校正。
Description
技术领域
本发明涉及测井技术的数据处理领域，特别涉及一种基于双尺度相关对比的自动深度校正方法及系统。
背景技术
地球物理测井通过多次组合测量来获取井剖面的各种地球物理参数，并以此来研究井剖面的地质问题和工程问题。然而，对于不同次组合测井测得的测井曲线，由于种种原因会有深度错动。所以，在对测井资料进行数字处理和综合解释时需要以某一次组合测井的深度为准，把其他各次的测井曲线与之深度对齐(即进行相对深度校正)。
测井曲线的相对深度校正可以用手工、人机交互和计算机自动处理来完成。其中，用计算机对测井曲线进行自动深度校正，是测井解释中需要解决的一个重要课题。
目前，自动深度校正方法主要采用以测井曲线间存在的相似性为依据进行的。众所周知，各类测井(如孔隙度测井、泥质指示测井、放射性测井、电阻率测井系列)曲线间通常具有一定的相似性，这才使利用某些数学方法由计算机实现测井曲线的自动校深成为可能。
实际应用中经常需要结合常规测井资料与电成像测井资料进行综合分析，而常规测井资料与电成像测井资料往往存在较大的深度偏差。相关函数法、图形识别法和频率分析法是目前较为成熟的几种常规测井资料与电成像测井资料之间自动深度校正方法，但它们在复杂条件下，尤其常规测井与电成像测井资料之间深度差较大时，存在校正精度不理想的问题。
发明内容
为了解决现有技术存在的在复杂条件下，尤其常规测井与电成像测井资料之间深度差较大时，校正精度低的缺点，提出一种基于双尺度相关对比的自动深度校正方法及系统，能够实现在常规测井资料与电成像测井资料之间精确校正。
一种基于双尺度相关对比的自动深度校正方法，其包括如下步骤：
S1、加载电成像测井资料；并对电成像测井资料进行预处理，得到电成像测井静态与动态图像；
S2、按照常规测井资料中的常规测井自然伽马曲线的采样间隔，对电成像测井静态与动态图像中的成像测井自然伽马曲线进行欠采样插值，得到第一尺度窗口下成像测井自然伽马曲线
S3、在第一尺度窗口下对成像测井自然伽马的校正曲线与常规测井曲线通过第一次对比进行粗校正，确定成像测井自然伽马曲线在第一尺度窗口下粗校正的深度范围；
S4、按照成像测井自然伽马曲线的采样间隔，对常规测井自然伽马曲线重采样，得到第二尺度窗口下常规自然伽马曲线第一尺度窗口大于第二尺度窗口；
S5、在第二尺度窗口中对成像测井自然伽马曲线与常规测井自然伽马曲线在第一尺度窗口下得到粗校正的深度范围内通过第二次对比进行精细校正，得到双尺度下成像测井自然伽马的校正曲线
S6、在得到双尺度下成像测井自然伽马的校正曲线后，其他成像测井曲线按照成像测井自然伽马的校正曲线规律进行校正，得到深度校正绘图。
一种基于双尺度相关对比的自动深度校正系统，其包括如下模块：
预处理模块，用于加载电成像测井资料；并对电成像测井资料进行预处理，得到电成像测井静态与动态图像；
第一采样模块，用于按照常规测井资料中的常规测井自然伽马曲线的采样间隔，对电成像测井静态与动态图像中的成像测井自然伽马曲线进行欠采样插值，得到第一尺度窗口下成像测井自然伽马曲线
粗校正模块，用于在第一尺度窗口下对成像测井自然伽马的校正曲线与常规测井曲线通过第一次对比进行粗校正，确定成像测井自然伽马曲线在第一尺度窗口下粗校正的深度范围；
第二采样模块，用于按照成像测井自然伽马曲线的采样间隔，对常规测井自然伽马曲线重采样，得到第二尺度窗口下常规自然伽马曲线第一尺度窗口大于第二尺度窗口；
精细校正模块，用于在第二尺度窗口中对成像测井自然伽马曲线与常规测井自然伽马曲线在第一尺度窗口下得到粗校正的深度范围内通过第二次对比进行精细校正，得到双尺度下成像测井自然伽马的校正曲线
其他成像测井曲线深度校正模块，用于在得到双尺度下成像测井自然伽马的校正曲线后，其他成像测井曲线按照成像测井自然伽马的校正曲线规律进行校正，得到深度校正绘图。
本发明提供的基于双尺度相关对比的自动深度校正方法及系统通过先在第一尺度窗口中确定了校正的大概深度范围，然后在第二尺度窗口中在粗校正深度范围内进行精细校正，第一尺度窗口大于第二尺度窗口，解决了复杂条件下深度准确校正问题，具有精度高，可靠性强的优点。
附图说明
图1是本发明实施的基于双尺度相关对比的自动深度校正方法流程图；
图2为基于双尺度相关对比的自动深度校正方法的示意图。
图3是本发明实施的基于双尺度相关对比的自动深度校正系统的结构框图。
具体实施方式
本发明实施例不限于应用的平台，比如可以通过Visual C++6.0开发平台实现以下实施例的内容。如图1所示，一种基于双尺度相关对比的自动深度校正方法，其包括如下步骤：
S1、加载电成像测井资料；并对电成像测井资料进行预处理，得到电成像测井静态与动态图像。
S2、按照常规测井资料中的常规测井自然伽马曲线的采样间隔，对电成像测井静态与动态图像中的成像测井自然伽马曲线进行欠采样插值，得到第一尺度窗口下成像测井自然伽马曲线
S3、在第一尺度窗口下对成像测井自然伽马的校正曲线与常规测井曲线通过第一次对比进行粗校正，确定成像测井自然伽马曲线在第一尺度窗口下粗校正的深度范围。如图2所示，图2中左边为第一尺度窗口下采样进行粗校正的示意图。
S4、按照成像测井自然伽马曲线的采样间隔，对常规测井自然伽马曲线重采样，得到第二尺度窗口下常规自然伽马曲线第一尺度窗口大于第二尺度窗口。
S5、在第二尺度窗口中对成像测井自然伽马曲线与常规测井自然伽马曲线在第一尺度窗口下得到粗校正的深度范围内通过第二次对比进行精细校正，得到双尺度下成像测井自然伽马的校正曲线如图2所示，图2中右边为第二尺度窗口下采样进行粗校正的示意图。
S6、在得到双尺度下成像测井自然伽马的校正曲线后，其他成像测井曲线按照成像测井自然伽马的校正曲线规律进行校正，得到深度校正绘图。
可选地，所述步骤S1中的预处理包括极板对齐处理以及均衡化处理。
可选地，所述步骤S3中第一次对比具体如下：
曲线与曲线间的相关函数为：
其中，是曲线上的第i个采样点的测井值，是曲线上的第i个采样点的测井值，t为相对于位移所对应的采样点数目，n为窗长内的采样点数，k为1/2探索长度内的采样点数,为曲线上进行相关对比曲线段上各采样点测井值的平均值，为曲线上进行对比井段内各采样点测井值的平均值。
可选地，所述步骤S5中第二次对比具体如下：
曲线与曲线间的相关函数为：
其中，是曲线上的第i个采样点的测井值，是曲线上的第i个采样点的测井值，t为相对于位移所对应的采样点数目，n为窗长内的采样点数，k为1/2探索长度内的采样点数，为曲线上进行相关对比曲线段上各采样点测井值的平均值,为曲线上进行对比井段内各采样点测井值的平均值。
如图3所示，一种基于双尺度相关对比的自动深度校正系统，其包括如下模块：
预处理模块10，用于加载电成像测井资料；并对电成像测井资料进行预处理，得到电成像测井静态与动态图像。
第一采样模块20，用于按照常规测井资料中的常规测井自然伽马曲线的采样间隔，对电成像测井静态与动态图像中的成像测井自然伽马曲线进行欠采样插值，得到第一尺度窗口下成像测井自然伽马曲线
粗校正模块30，用于在第一尺度窗口下对成像测井自然伽马的校正曲线与常规测井曲线通过第一次对比进行粗校正，确定成像测井自然伽马曲线在第一尺度窗口下粗校正的深度范围。
第二采样模块40，用于按照成像测井自然伽马曲线的采样间隔，对常规测井自然伽马曲线重采样，得到第二尺度窗口下常规自然伽马曲线第一尺度窗口大于第二尺度窗口。
精细校正模块50，用于在第二尺度窗口中对成像测井自然伽马曲线与常规测井自然伽马曲线在第一尺度窗口下得到粗校正的深度范围内通过第二次对比进行精细校正，得到双尺度下成像测井自然伽马的校正曲线
其他成像测井曲线深度校正模块60，用于在得到双尺度下成像测井自然伽马的校正曲线后，其他成像测井曲线按照成像测井自然伽马的校正曲线规律进行校正，得到深度校正绘图。
可选地，所述预处理模块10中的预处理包括极板对齐处理以及均衡化处理。
可选地，所述粗校正模块30中第一次对比具体如下：
曲线与曲线间的相关函数为：
其中，是曲线上的第i个采样点的测井值，是曲线上的第i个采样点的测井值，t为相对于位移所对应的采样点数目，n为窗长内的采样点数，k为1/2探索长度内的采样点数,为曲线上进行相关对比曲线段上各采样点测井值的平均值，为曲线上进行对比井段内各采样点测井值的平均值。
可选地，所述精细校正模块50中第二次对比具体如下：
曲线与曲线间的相关函数为：
其中，是曲线上的第i个采样点的测井值，是曲线上的第i个采样点的测井值，t为相对于位移所对应的采样点数目，n为窗长内的采样点数，k为1/2探索长度内的采样点数，为曲线上进行相关对比曲线段上各采样点测井值的平均值,为曲线上进行对比井段内各采样点测井值的平均值。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机储存器、内存、只读存储器、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其他形式的存储介质中。
可以理解的是，对于本领域的普通技术人员来说，可以根据本发明的技术构思做出其它各种相应的改变与变形，而所有这些改变与变形都应属于本发明权利要求的保护范围。
Claims (8)
1.一种基于双尺度相关对比的自动深度校正方法，其特征在于，其包括如下步骤：
S1、加载电成像测井资料；并对电成像测井资料进行预处理，得到电成像测井静态与动态图像；
S2、按照常规测井资料中的常规测井自然伽马曲线的采样间隔，对电成像测井静态与动态图像中的成像测井自然伽马曲线进行欠采样插值，得到第一尺度窗口下成像测井自然伽马曲线
S3、在第一尺度窗口下对成像测井自然伽马的校正曲线与常规测井曲线通过第一次对比进行粗校正，确定成像测井自然伽马曲线在第一尺度窗口下粗校正的深度范围；
S4、按照成像测井自然伽马曲线的采样间隔，对常规测井自然伽马曲线重采样，得到第二尺度窗口下常规自然伽马曲线第一尺度窗口大于第二尺度窗口；
S5、在第二尺度窗口中对成像测井自然伽马曲线与常规测井自然伽马曲线在第一尺度窗口下得到粗校正的深度范围内通过第二次对比进行精细校正，得到双尺度下成像测井自然伽马的校正曲线
S6、在得到双尺度下成像测井自然伽马的校正曲线后，其他成像测井曲线按照成像测井自然伽马的校正曲线规律进行校正，得到深度校正绘图。
2.根据权利要求1所述的基于双尺度相关对比的自动深度校正方法，其特征在于，所述步骤S1中的预处理包括极板对齐处理以及均衡化处理。
3.根据权利要求1所述的基于双尺度相关对比的自动深度校正方法，其特征在于，所述步骤S3中第一次对比具体如下：
曲线与曲线间的相关函数为：
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>i</mi>
<mo>+</mo>
<mi>t</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>i</mi>
<mo>+</mo>
<mi>t</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
</mrow>
其中，是曲线上的第i个采样点的测井值，是曲线上的第i个采样点的测井值，t为相对于位移所对应的采样点数目，n为窗长内的采样点数，k为1/2探索长度内的采样点数,为曲线上进行相关对比曲线段上各采样点测井值的平均值，为曲线上进行对比井段内各采样点测井值的平均值。
4.根据权利要求3所述的基于双尺度相关对比的自动深度校正方法，其特征在于，所述步骤S5中第二次对比具体如下：
曲线与曲线间的相关函数为：
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>i</mi>
<mo>+</mo>
<mi>t</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>i</mi>
<mo>+</mo>
<mi>t</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
</mrow>
其中，是曲线上的第i个采样点的测井值，是曲线上的第i个采样点的测井值，t为相对于位移所对应的采样点数目，n为窗长内的采样点数，k为1/2探索长度内的采样点数，为曲线上进行相关对比曲线段上各采样点测井值的平均值,为曲线上进行对比井段内各采样点测井值的平均值。
5.一种基于双尺度相关对比的自动深度校正系统，其特征在于，其包括如下模块：
预处理模块，用于加载电成像测井资料；并对电成像测井资料进行预处理，得到电成像测井静态与动态图像；
第一采样模块，用于按照常规测井资料中的常规测井自然伽马曲线的采样间隔，对电成像测井静态与动态图像中的成像测井自然伽马曲线进行欠采样插值，得到第一尺度窗口下成像测井自然伽马曲线
粗校正模块，用于在第一尺度窗口下对成像测井自然伽马的校正曲线与常规测井曲线通过第一次对比进行粗校正，确定成像测井自然伽马曲线在第一尺度窗口下粗校正的深度范围；
第二采样模块，用于按照成像测井自然伽马曲线的采样间隔，对常规测井自然伽马曲线重采样，得到第二尺度窗口下常规自然伽马曲线第一尺度窗口大于第二尺度窗口；
精细校正模块，用于在第二尺度窗口中对成像测井自然伽马曲线与常规测井自然伽马曲线在第一尺度窗口下得到粗校正的深度范围内通过第二次对比进行精细校正，得到双尺度下成像测井自然伽马的校正曲线
其他成像测井曲线深度校正模块，用于在得到双尺度下成像测井自然伽马的校正曲线后，其他成像测井曲线按照成像测井自然伽马的校正曲线规律进行校正，得到深度校正绘图。
6.根据权利要求5所述的基于双尺度相关对比的自动深度校正系统，其特征在于，所述预处理模块中的预处理包括极板对齐处理以及均衡化处理。
7.根据权利要求6所述的基于双尺度相关对比的自动深度校正系统，其特征在于，所述粗校正模块中第一次对比具体如下：
曲线与曲线间的相关函数为：
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>i</mi>
<mo>+</mo>
<mi>t</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>i</mi>
<mo>+</mo>
<mi>t</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
</mrow>
其中，是曲线上的第i个采样点的测井值，是曲线上的第i个采样点的测井值，t为相对于位移所对应的采样点数目，n为窗长内的采样点数，k为1/2探索长度内的采样点数,为曲线上进行相关对比曲线段上各采样点测井值的平均值，为曲线上进行对比井段内各采样点测井值的平均值。
8.根据权利要求7所述的基于双尺度相关对比的自动深度校正系统，其特征在于，所述精细校正模块中第二次对比具体如下：
曲线与曲线间的相关函数为：
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>i</mi>
<mo>+</mo>
<mi>t</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mi>i</mi>
</mrow>
<mn>1</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mi>i</mi>
<mo>+</mo>
<mi>t</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>-</mo>
<mover>
<mrow>
<msubsup>
<mi>GR</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
</mrow>
其中，是曲线上的第i个采样点的测井值，是曲线上的第i个采样点的测井值，t为相对于位移所对应的采样点数目，n为窗长内的采样点数，k为1/2探索长度内的采样点数，为曲线上进行相关对比曲线段上各采样点测井值的平均值,为曲线上进行对比井段内各采样点测井值的平均值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510188038.0A CN104832161B (zh) | 2015-04-20 | 2015-04-20 | 基于双尺度相关对比的自动深度校正方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510188038.0A CN104832161B (zh) | 2015-04-20 | 2015-04-20 | 基于双尺度相关对比的自动深度校正方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104832161A CN104832161A (zh) | 2015-08-12 |
| CN104832161B true CN104832161B (zh) | 2017-12-29 |
Family
ID=53810304
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510188038.0A Active CN104832161B (zh) | 2015-04-20 | 2015-04-20 | 基于双尺度相关对比的自动深度校正方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104832161B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105484732B (zh) * | 2015-12-28 | 2018-10-12 | 上海神开石油设备有限公司 | 用于水平井钻井地质导向施工过程井深的处理方法 |
| CN108571316B (zh) * | 2017-03-14 | 2023-01-13 | 中石化石油工程技术服务有限公司 | 一种无缆测井深度校正方法及装置 |
| CN108828687B (zh) * | 2018-08-09 | 2020-07-24 | 中国海洋石油集团有限公司 | 一种基于电成像面孔率的渗透率计算方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2591758B1 (fr) * | 1985-12-16 | 1988-05-20 | Total Petroles | Procede de correction de la profondeur de diagraphies |
| US6618675B2 (en) * | 2001-02-27 | 2003-09-09 | Halliburton Energy Services, Inc. | Speed correction using cable tension |
| CN102003172B (zh) * | 2010-10-11 | 2013-04-03 | 青岛杰瑞自动化有限公司 | 测井系统中的测井深度校正计算方法 |
| CN202370504U (zh) * | 2011-12-23 | 2012-08-08 | 东北石油大学 | 一种射孔深度自动校正装置 |
| CN104121009B (zh) * | 2013-04-26 | 2017-02-01 | 中国石油天然气股份有限公司 | 一种温度、压力曲线的生成方法及系统 |
| CN104156624B (zh) * | 2014-08-27 | 2017-08-01 | 中国电子科技集团公司第二十二研究所 | 测井设备的时深数据的处理方法和装置 |
-
2015
- 2015-04-20 CN CN201510188038.0A patent/CN104832161B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104832161A (zh) | 2015-08-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sokos et al. | Evaluating centroid‐moment‐tensor uncertainty in the new version of ISOLA software | |
| JP6610678B2 (ja) | ピーク検出方法及びデータ処理装置 | |
| EP1918877A2 (en) | Correlation peak finding method for image correlation displacement sensing | |
| CN105116448B (zh) | 一种转换波方位各向异性校正方法及装置 | |
| CN104832161B (zh) | 基于双尺度相关对比的自动深度校正方法及系统 | |
| US20230401691A1 (en) | Image defect detection method, electronic device and readable storage medium | |
| CN113536574B (zh) | 一种地质导向建模方法、系统及计算机可读存储介质 | |
| Xu et al. | Systemic comparison of seismometer horizontal orientations based on teleseismic earthquakes and ambient‐noise data | |
| CN107075567B (zh) | 用于荧光数据校正的方法 | |
| CN112859192B (zh) | 岩心饱和度预测模型构建方法与岩心饱和度预测方法 | |
| CN106199693A (zh) | 地震数据速度谱自动拾取方法和装置 | |
| CN111596379A (zh) | 一种地震观测系统异常识别方法和装置 | |
| CN104278990B (zh) | 测井数据质量恢复方法及装置 | |
| CN110907984A (zh) | 一种基于自回归移动平均模型的震前红外长波辐射异常信息探测方法 | |
| CN111898314B (zh) | 湖泊水体参数检验方法、装置、电子设备及存储介质 | |
| de Souza et al. | Using parity decomposition for interpreting magnetic anomalies from dikes having arbitrary dip angles, induced and remanent magnetization | |
| CN103364023A (zh) | 一种基于总体最小二乘的仪表图像中心点提取方法 | |
| CN108376123A (zh) | 一种卫星重力观测数据异常值提取方法、装置及电子设备 | |
| Lyubushin et al. | Prognostic Properties of Instantaneous Amplitudes Maxima of Earth Surface Tremor | |
| CN115541021A (zh) | 拉曼光谱特征峰的定位方法、电子设备及存储介质 | |
| CN116359925B (zh) | 基于水深匹配的浅剖位置坐标校正方法 | |
| CN113296157B (zh) | 一种利用广义高斯分布进行储层的预测方法及装置 | |
| Zheleznyak et al. | Marine gravity measurements without referencing to the coastal base stations | |
| CN111177886B (zh) | 基于物探数据分析的海上测距规划及土厚预测的方法 | |
| Al Ibrahim et al. | Automated Verification of Sidewall Core Recovery Depth using Borehole Image Logs |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |