CN111980685B - 测井曲线处理方法及装置 - Google Patents
测井曲线处理方法及装置 Download PDFInfo
- Publication number
- CN111980685B CN111980685B CN201910427496.3A CN201910427496A CN111980685B CN 111980685 B CN111980685 B CN 111980685B CN 201910427496 A CN201910427496 A CN 201910427496A CN 111980685 B CN111980685 B CN 111980685B
- Authority
- CN
- China
- Prior art keywords
- reflecting
- logging curve
- log
- logging
- lithology
- 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
- 238000003672 processing method Methods 0.000 title abstract description 7
- 230000008859 change Effects 0.000 claims abstract description 85
- 230000008021 deposition Effects 0.000 claims abstract description 46
- 238000013517 stratification Methods 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 17
- 239000013049 sediment Substances 0.000 claims abstract description 15
- 238000004590 computer program Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 11
- 230000004927 fusion Effects 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 11
- 238000004519 manufacturing process Methods 0.000 abstract description 11
- 239000011435 rock Substances 0.000 abstract description 9
- 230000015572 biosynthetic process Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 230000000630 rising effect Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000032798 delamination Effects 0.000 description 3
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035939 shock 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种测井曲线处理方法及装置，其中方法包括：获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。本发明将反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，融合后的结果能够全面完整地体现地层变化的轮廓与细节特征，在利用融合后的结果与岩心进行比对时，有利于建立更加可靠的匹配关系，进而有效提高地质分层结果的准确性，为油田开发生产奠定基础。
Description
技术领域
本发明涉及石油测井工程技术领域，尤其涉及测井曲线处理方法及装置。
背景技术
随着油田勘探开发进入中后期阶段，研究目标从简单构造油藏向复杂断块、复杂地层和复杂岩性方向转变，而这些油藏受到多次构造活动、纵向岩性混杂、横向地层尖灭等影响，基础地质分层工作难度较大。现有技术中通常直接利用测井曲线与岩心进行比对，确定测井曲线幅度、形态等地层特征与岩性、沉积旋回的匹配关系，从而利用匹配关系对地质进行分层。但是对于一些地质条件复杂的油田区块，测井曲线的地层特征不明显，获得的地质分层结果准确性较低，不利于油田开发生产。
发明内容
本发明实施例提供一种测井曲线处理方法，用以提高地质分层结果的准确性，为油田开发生产奠定基础，该方法包括：
获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；
对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。
本发明实施例提供一种测井曲线处理装置，用以提高地质分层结果的准确性，为油田开发生产奠定基础，该装置包括：
获取模块，获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；
融合模块，对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如下方法：
获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；
对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行如下方法的计算机程序：
获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；
对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。
相对于现有技术中直接利用测井曲线与岩心进行比对，确定测井曲线幅度、形态等地层特征与岩性、沉积旋回的匹配关系，利用匹配关系对地质进行分层的方案而言，本发明实施例通过获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线，对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。本发明实施例将反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，融合后的结果能够全面完整地体现地层变化的轮廓与细节特征，在利用融合后的结果与岩心进行比对时，有利于建立更加可靠的匹配关系，进而有效提高地质分层结果的准确性，为油田开发生产奠定基础。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中测井曲线处理方法示意图；
图2为本发明实施例中曲线形态定义示意图；
图3为本发明实施例中特征曲线与分层曲线结果示意图；
图4为本发明实施例中测井曲线处理装置结构图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
发明人发现，测井曲线可以分为反映整体沉积背景的测井曲线和反映岩性变化的测井曲线，反映整体沉积背景的测井曲线体现地层变化的轮廓特征，反映岩性变化的测井曲线体现地层变化的细节特征，轮廓特征和细节特征则共同构成了地层发生变化时的特征。在现有技术中，通常是将单一测井曲线与岩心进行比对，确定测井曲线幅度、形态等地层特征与岩性、沉积旋回的匹配关系，从而利用匹配关系对地质进行分层。而单一测井曲线无法全面完整地体现地层变化的轮廓与细节特征，因此在与岩心进行比对时无法建立可靠的匹配关系，进而获得的地质分层结果准确性较低，不利于油田开发生产。但是，如果将各测井曲线与岩心进行逐一比对，工作量繁重且工作效率低。
为了提高地质分层结果的准确性，为油田开发生产奠定基础，本发明实施例提供一种测井曲线处理方法，如图1所示，该方法可以包括：
步骤101、获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；
步骤102、对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。
由图1所示可以得知，本发明实施例通过获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线，对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。本发明实施例将反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，融合后的结果能够全面完整地体现地层变化的轮廓与细节特征，在利用融合后的结果与岩心进行比对时，有利于建立更加可靠的匹配关系，进而有效提高地质分层结果的准确性，为油田开发生产奠定基础。
具体实施时，获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线。发明人发现，测井曲线可以分为反映整体沉积背景的测井曲线和反映岩性变化的测井曲线，反映整体沉积背景的测井曲线体现地层变化的轮廓特征，反映岩性变化的测井曲线体现地层变化的细节特征，轮廓特征和细节特征则共同构成了地层发生变化时的特征。在现有技术中，通常是将单一测井曲线与岩心进行比对，确定测井曲线幅度、形态等地层特征与岩性、沉积旋回的匹配关系，从而利用匹配关系对地质进行分层。而单一测井曲线无法全面完整地体现地层变化的轮廓与细节特征，因此在与岩心进行比对时无法建立可靠的匹配关系，进而获得的地质分层结果准确性较低，不利于油田开发生产。因此，本发明实施例获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线，避免遗漏地层变化的轮廓或细节特征。
实施例中，反映整体沉积背景的测井曲线包括：时差测井曲线，密度测井曲线和自然电位测井曲线其中之一或任意组合。
实施例中，反映岩性变化的测井曲线包括：深侧向测井曲线，浅侧向测井曲线和自然伽马测井曲线其中之一或任意组合。
实施例中，获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线之后，分别从反映整体沉积背景的测井曲线和反映岩性变化的测井曲线中选取至少一个测井曲线，然后在对所述测井曲线进行融合时，对选取的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合。发明人发现，在反映整体沉积背景的测井曲线和反映岩性变化的测井曲线中，存在一些对地质分层没有贡献作用的测井曲线，也即测井曲线在两不同地层的值相同，如果将这些测井曲线一并进行融合，不仅不能提高地质分层的准确率，而且在融合结果中增加噪声，影响分层结果的准确性。因此，本发明实施例从反映整体沉积背景的测井曲线和反映岩性变化的测井曲线中进行选取，对选取的测井曲线进行融合，保证进行融合的测井曲线对地质分层有贡献作用，提高分层准确率。
实施例中，按如下方式选取至少一个测井曲线：首先，获取已知区块的测井曲线，然后根据岩心相资料，确定已知区块的地质分层结果，根据所述地质分层结果，确定所述已知区块的测井曲线在地质分层界线处的变化趋势，最后根据所述变化趋势，从所述已知区块的测井曲线中选取至少一个测井曲线作为待测区块的测井曲线。其中，岩心相资料包括钻井取心资料、古生物分析化验资料、同位素检验资料、井震对比结果。变化趋势包括：幅度、斜率、曲率和形态的变化趋势。在确定已知区块的地质分层结果之后，可以找到地质分层界线，从而对应找到已知区块的测井曲线的地质分层界线处。通过分析已知区块的测井曲线在地质分层界线处的幅度、斜率、曲率和形态的变化趋势，从所述已知区块的测井曲线中选取至少一个测井曲线作为待测区块的测井曲线，例如，可以选取已知区块的测井曲线在地质分层界线处的幅度、斜率、曲率或形态具有差异的测井曲线作为待测区块的测井曲线。
实施例中，获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线之后，按如下方式对测井曲线进行处理：对测井曲线进行去噪，对去噪后的测井曲线进行归一化处理，在对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合时，对归一化处理后的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合。
实施例中，可以利用Savitzky-Golay(S-G)滤波算法按如下公式对曲线进行滤波，消除测井曲线的噪声：
其中，Y为去噪前的测井曲线，为第j个去噪后的测井曲线，Ci为S-G多项式拟合系数，表示第i个测井曲线的权重，N为低通滤波器的窗口长度，其长度为2m+1。
实施例中，在对测井曲线进行去噪之后，对于不同测井数据量纲不同的情况，按如下公式对去噪后的测井曲线进行归一化处理：
其中，Y'为归一化处理后的结果，Y*为归一化处理之前的测井曲线，为归一化处理之前测井曲线的最小值，/>为归一化处理之前测井曲线的最大值。
具体实施时，对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。
实施例中，可以按如下方式对反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合：首先对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行频域正变换，然后对反映整体沉积背景的测井曲线对应的频域信号进行低通滤波，对反映岩性变化的测井曲线对应的频域信号进行高通滤波，最后将低通滤波后的结果与高通滤波滤波后的结果相加，对相加后的结果进行频域反变换，获得融合后的测井曲线。发明人发现，反映整体沉积背景的测井曲线体现地层变化的轮廓特征，有用信号集中在低频，噪声主要集中在高频，而反映岩性变化的测井曲线体现地层变化的细节特征，有用信号集中在高频，噪声主要集中在低频。因此本发明实施例对反映整体沉积背景的测井曲线对应的频域信号进行低通滤波，对反映岩性变化的测井曲线对应的频域信号进行高通滤波，可以有效去除两类测井曲线中的噪声，从而有利于提高分层准确率。
实施例中，在进行融合之后，对融合后的曲线进行符号化，根据曲线相邻值的幅度差赋予符号值。首先定义曲线值幅度差异阈值为ε，若设定区间内曲线值幅度差Δx＞ε，则这一区间的曲线定义为上升；若设定区间内曲线值幅度差Δx＜-ε，则这一区间的曲线定义为下降；若|Δx|≤ε，则这一区间的曲线定义为平稳。将相邻区间的相同符号进行合并，将合并后的相邻符号按照曲线形态进行定义，如图2所示，上升用①表示，下降用②表示，平稳用③表示，曲线形态分为6类特征，其中：将“上升→下降”(“①→②”)定义为峰，如图2中的A所示；将“下降→上升”(“②→①”)定义为谷，如图2中的B所示；将“上升→平稳”(“①→③”)定义为上升-平直，如图2中的C所示；将“平稳→下降”(“③→②”)定义为平直-下降，如图2中的D所示；将“平稳→上升”(“③→①”)定义为平直-上升，如图2中的E所示；将“下降→平稳”(“②→③”)定义为下降-平直，如图2中的F所示，从而获得特征曲线。基于获得的特征，获取用于地质分层的分层曲线。通过符号曲线对应深度值，读取符号对应特征曲线上的数值。在相邻符号对应的深度区间内，按如下公式以6类符号对应深度区间在特征曲线上响应数据点的平均值充填：
其中，为填充的平均值，X(a,b)为深度区间内的特征曲线值之和，S(a,b)为深度区间内的采样点数。结果如图3所示，其中图3中的A为特征曲线结果，图3中的B为分层曲线结果。
基于同一发明构思，本发明实施例还提供了一种测井曲线处理装置，如下面的实施例所述。由于这些解决问题的原理与测井曲线处理方法相似，因此装置的实施可以参见方法的实施，重复之处不再赘述。
图4为本发明实施例中测井曲线处理装置的结构图，如图4所示，该装置包括：
获取模块401，获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；
融合模块402，对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。
一个实施例中，获取模块401进一步用于：
获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线之后，分别从反映整体沉积背景的测井曲线和反映岩性变化的测井曲线中选取至少一个测井曲线；
所述融合模块进一步用于：对选取的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合。
一个实施例中，获取模块401进一步用于，按如下方式选取至少一个测井曲线：
获取已知区块的测井曲线；
根据岩心相资料，确定已知区块的地质分层结果，根据所述地质分层结果，确定所述已知区块的测井曲线在地质分层界线处的变化趋势；
根据所述变化趋势，从所述已知区块的测井曲线中选取至少一个测井曲线作为待测区块的测井曲线。
一个实施例中，融合模块402进一步用于，按如下方式对反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合：
对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行频域正变换；
对反映整体沉积背景的测井曲线对应的频域信号进行低通滤波，对反映岩性变化的测井曲线对应的频域信号进行高通滤波；
将低通滤波后的结果与高通滤波滤波后的结果相加，对相加后的结果进行频域反变换，获得融合后的测井曲线。
综上所述，本发明实施例通过获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线，对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层。本发明实施例将反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，融合后的结果能够全面完整地体现地层变化的轮廓与细节特征，在利用融合后的结果与岩心进行比对时，有利于建立更加可靠的匹配关系，进而有效提高地质分层结果的准确性，为油田开发生产奠定基础。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (12)
1.一种测井曲线处理方法，其特征在于，包括：
获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；
对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层；
按如下方式对反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合：
对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行频域正变换；
对反映整体沉积背景的测井曲线对应的频域信号进行低通滤波，对反映岩性变化的测井曲线对应的频域信号进行高通滤波；
将低通滤波后的结果与高通滤波滤波后的结果相加，对相加后的结果进行频域反变换，获得融合后的测井曲线。
2.如权利要求1所述的方法，其特征在于，所述反映整体沉积背景的测井曲线包括：时差测井曲线，密度测井曲线和自然电位测井曲线其中之一或任意组合。
3.如权利要求1所述的方法，其特征在于，所述反映岩性变化的测井曲线包括：深侧向测井曲线，浅侧向测井曲线和自然伽马测井曲线其中之一或任意组合。
4.如权利要求1所述的方法，其特征在于，获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线之后，分别从反映整体沉积背景的测井曲线和反映岩性变化的测井曲线中选取至少一个测井曲线；
对所述测井曲线进行融合包括：对选取的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合。
5.如权利要求4所述的方法，其特征在于，按如下方式选取至少一个测井曲线：
获取已知区块的测井曲线；
根据岩心相资料，确定已知区块的地质分层结果，根据所述地质分层结果，确定所述已知区块的测井曲线在地质分层界线处的变化趋势；
根据所述变化趋势，从所述已知区块的测井曲线中选取至少一个测井曲线作为待测区块的测井曲线。
6.如权利要求5所述的方法，其特征在于，所述变化趋势包括：幅度，斜率，曲率和形态的变化趋势。
7.如权利要求1所述的方法，其特征在于，获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线之后，按如下方式对测井曲线进行处理：对测井曲线进行去噪，对去噪后的测井曲线进行归一化处理；
对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合包括：对归一化处理后的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合。
8.一种测井曲线处理装置，其特征在于，包括：
获取模块，获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线；
融合模块，对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合，根据融合后的结果对所述待测区块进行地质分层；
所述融合模块进一步用于，按如下方式对反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合：
对所述反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行频域正变换；
对反映整体沉积背景的测井曲线对应的频域信号进行低通滤波，对反映岩性变化的测井曲线对应的频域信号进行高通滤波；
将低通滤波后的结果与高通滤波滤波后的结果相加，对相加后的结果进行频域反变换，获得融合后的测井曲线。
9.如权利要求8所述的装置，其特征在于，所述获取模块进一步用于：
获取待测区块的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线之后，分别从反映整体沉积背景的测井曲线和反映岩性变化的测井曲线中选取至少一个测井曲线；
所述融合模块进一步用于：对选取的反映整体沉积背景的测井曲线和反映岩性变化的测井曲线进行融合。
10.如权利要求9所述的装置，其特征在于，所述获取模块进一步用于，按如下方式选取至少一个测井曲线：
获取已知区块的测井曲线；
根据岩心相资料，确定已知区块的地质分层结果，根据所述地质分层结果，确定所述已知区块的测井曲线在地质分层界线处的变化趋势；
根据所述变化趋势，从所述已知区块的测井曲线中选取至少一个测井曲线作为待测区块的测井曲线。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至7任一所述方法。
12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至7任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910427496.3A CN111980685B (zh) | 2019-05-22 | 2019-05-22 | 测井曲线处理方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910427496.3A CN111980685B (zh) | 2019-05-22 | 2019-05-22 | 测井曲线处理方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111980685A CN111980685A (zh) | 2020-11-24 |
| CN111980685B true CN111980685B (zh) | 2023-08-22 |
Family
ID=73436319
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910427496.3A Active CN111980685B (zh) | 2019-05-22 | 2019-05-22 | 测井曲线处理方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111980685B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102353985A (zh) * | 2011-06-01 | 2012-02-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于非下采样Contourlet变换的拟声波曲线构建方法 |
| CN104295293A (zh) * | 2014-10-23 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种获取测井密度曲线的方法 |
| CN104462792A (zh) * | 2014-11-24 | 2015-03-25 | 中国海洋石油总公司 | 一种测井数据岩性层数值归约方法 |
| CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105426612A (zh) * | 2015-11-18 | 2016-03-23 | 中国石油天然气股份有限公司 | 一种地层组分最优化确定方法及装置 |
| CN109113729A (zh) * | 2018-06-20 | 2019-01-01 | 中国石油天然气集团有限公司 | 基于测井曲线的岩性识别方法及装置 |
-
2019
- 2019-05-22 CN CN201910427496.3A patent/CN111980685B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102353985A (zh) * | 2011-06-01 | 2012-02-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于非下采样Contourlet变换的拟声波曲线构建方法 |
| CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN104295293A (zh) * | 2014-10-23 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种获取测井密度曲线的方法 |
| CN104462792A (zh) * | 2014-11-24 | 2015-03-25 | 中国海洋石油总公司 | 一种测井数据岩性层数值归约方法 |
| CN105426612A (zh) * | 2015-11-18 | 2016-03-23 | 中国石油天然气股份有限公司 | 一种地层组分最优化确定方法及装置 |
| CN109113729A (zh) * | 2018-06-20 | 2019-01-01 | 中国石油天然气集团有限公司 | 基于测井曲线的岩性识别方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111980685A (zh) | 2020-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106951660B (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
| CN109709607B (zh) | 一种预测薄层砂岩储层厚度方法及装置 | |
| CN107748399B (zh) | 利用重力界面反演识别山前带深部构造层方法 | |
| CN111856566B (zh) | 湖相滩坝砂体中薄储层预测方法及装置 | |
| CN111624651B (zh) | 基于古地貌约束的储层预测方法及装置 | |
| EP2616978A1 (en) | System and method for sweet zone identification in shale gas reservoirs | |
| CN108629459B (zh) | 储层含烃孔隙的检测方法及装置 | |
| CN108873065B (zh) | 砂岩优质储层预测方法及装置 | |
| CN105240006A (zh) | 一种适用于火山岩储层的油水层识别方法 | |
| CN111983671B (zh) | 基于微古地貌恢复的浅水湖盆储层预测方法和装置 | |
| KR101292659B1 (ko) | 물리 검층 데이터의 해석 방법 | |
| EP2926171B1 (en) | System and method for seismic velocity anomaly analysis | |
| NO20121473A1 (no) | System for modellering av geologiske strukturer | |
| CN103969685A (zh) | 一种薄互层地震信号的处理方法 | |
| CN110632654B (zh) | 断块圈闭含油边界确定方法及装置 | |
| CN111980685B (zh) | 测井曲线处理方法及装置 | |
| CN112379462A (zh) | 电磁地震数据联合处理方法及装置 | |
| CN111239809B (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
| CN112485841A (zh) | 深层地层岩性识别方法及装置 | |
| CN112649887A (zh) | 基于钻井资料定量划分层序地层的方法及装置 | |
| Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis | |
| CN115113268A (zh) | 一种基于米氏旋回的三级层序划分方法 | |
| CN117270032A (zh) | 基于三维地震的砂岩型铀矿铀特征参数反演方法及装置 | |
| CN105467447B (zh) | 相控趋势能量匹配的地震保幅评价方法 | |
| CN110954947B (zh) | 时域线性低频融合方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |