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
- mover
- imaging logging
- 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
- 230000000052 comparative effect Effects 0.000 title claims abstract description 29
- 238000003384 imaging method Methods 0.000 claims abstract description 100
- 230000005251 gamma ray Effects 0.000 claims abstract description 58
- 238000012937 correction Methods 0.000 claims abstract description 44
- 238000005070 sampling Methods 0.000 claims abstract description 44
- 238000011088 calibration curve Methods 0.000 claims abstract description 25
- 238000012952 Resampling Methods 0.000 claims abstract description 6
- 230000003068 static effect Effects 0.000 claims description 11
- 238000005314 correlation function Methods 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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 |
---|---|---|
Austermann et al. | Detection of a dynamic topography signal in last interglacial sea-level records | |
Sokos et al. | Evaluating centroid‐moment‐tensor uncertainty in the new version of ISOLA software | |
Yang et al. | Simultaneous earthquake detection on multiple stations via a convolutional neural network | |
CN104280777B (zh) | 一种压制陆上地震资料多次波干扰的方法 | |
CN105607122B (zh) | 一种基于全变分地震数据分解模型的地震纹理提取与增强方法 | |
CN104832161B (zh) | 基于双尺度相关对比的自动深度校正方法及系统 | |
MacLeod | Punctuated anagenesis and the importance of stratigraphy to paleobiology | |
CN113267822B (zh) | 基于地形约束因子权重优化提高海底地形反演精度的方法 | |
CN111381275A (zh) | 一种地震数据的初至拾取方法和装置 | |
Yu et al. | LPPN: A lightweight network for fast phase picking | |
Pilz et al. | Data-driven and machine learning identification of seismic reference stations in Europe | |
Spasojevic | Statistical analysis of ground‐based chorus observations during geomagnetic storms | |
CN108957554A (zh) | 一种地球物理勘探中地震反演方法 | |
Bhardwaj et al. | Multi-parameter algorithm for earthquake early warning | |
Wilson et al. | Detection and characterization of pulses in broadband seismometers | |
WO2023123971A1 (zh) | 基于vsp的深度域地震剖面层位标定方法及装置 | |
CN116482757A (zh) | 基于浅剖井震联合的地层结构识别分析方法 | |
US11927713B2 (en) | Method, device and medium for acquiring logging parameters | |
CN109857979A (zh) | 一种基于小波分析和bp神经网络的测井曲线拟合方法 | |
Castellaro et al. | A statistical low noise model of the earth | |
CN115113268A (zh) | 一种基于米氏旋回的三级层序划分方法 | |
CN116009080A (zh) | 一种地震波阻抗反演方法及系统、电子设备、存储介质 | |
Brandt | Performance of the South African National Seismograph Network from October 2012 to February 2017: spatially varying magnitude completeness | |
CN112666610B (zh) | 地层有机碳含量计算方法及系统 | |
Yamada | Estimation of fault rupture extent using near-source records for earthquake early warning |
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 |