CN111827966A - 一种多井声波测井曲线一致性处理方法、装置及存储介质 - Google Patents

一种多井声波测井曲线一致性处理方法、装置及存储介质 Download PDF

Info

Publication number
CN111827966A
CN111827966A CN202010216043.9A CN202010216043A CN111827966A CN 111827966 A CN111827966 A CN 111827966A CN 202010216043 A CN202010216043 A CN 202010216043A CN 111827966 A CN111827966 A CN 111827966A
Authority
CN
China
Prior art keywords
well
wells
acoustic logging
inversion
frequency
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.)
Granted
Application number
CN202010216043.9A
Other languages
English (en)
Other versions
CN111827966B (zh
Inventor
于承业
韦学锐
陈鸿安
于浩忙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Daqing Oilfield Co Ltd
Original Assignee
Daqing Oilfield Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Daqing Oilfield Co Ltd filed Critical Daqing Oilfield Co Ltd
Priority to CN202010216043.9A priority Critical patent/CN111827966B/zh
Publication of CN111827966A publication Critical patent/CN111827966A/zh
Application granted granted Critical
Publication of CN111827966B publication Critical patent/CN111827966B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Mining & Mineral Resources (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Remote Sensing (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种多井声波测井曲线一致性处理方法、装置及存储介质，涉及油田领域，其中，一种多井声波测井曲线一致性处理方法，包括：获取所有井的声波测井曲线；根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线；根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；所述每口井反演结果均达到第一设定反演结果符合率，则根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果；利用验证井检验多井反演结果，所述验证井的反演结果达到第二设定反演结果符合率，则完成多井声波测井曲线一致性处理。以解决多井测井曲线的一致性差，多井约束反演及多井建模结果失真的问题。

Description

一种多井声波测井曲线一致性处理方法、装置及存储介质
技术领域
本发明涉及油田领域，具体说是一种多井声波测井曲线一致性处理方法、装置及存储介质。
背景技术
多井反演及多井建模的关键之一是多井声波曲线一致性处理方法，在利用实际资料来做地震波阻抗反演时，波阻抗结果是根据声波和密度测井资料来求取的。因此，测井数据的相对保真及多井一致性是十分重要的，如果测井数据不准确，不仅会导致的到得合成地震记录不精确，与井旁道的匹配程度较差，影响地震、地质层位标定，地层对比等，当然导致经约束反演结果失真；另一方面，即使每口井测井曲线相对保真，做多井约束反演时，由于测量时间跨度大、所选仪器性能与型号的不同，使用泥浆不一致、井径不一致等也会造成多井测井曲线的一致性差，不同井约束的反演剖面差异较大，多井约束反演及多井建模结果失真。
所以，多井约束反演及建模，首先要对测井资料对其进行预必要的预处理，来消除非地质因素对测井曲线数据的影响，使得测井资料和地震资料相关性更高。
发明内容
有鉴于此，本发明提供一种多井声波测井曲线一致性处理方法、装置及存储介质，以解决多井测井曲线的一致性差，不同井约束的反演剖面差异较大，多井约束反演及多井建模结果失真的问题。
第一方面，本发明提供一种多井声波测井曲线一致性处理方法，包括：
获取所有井的声波测井曲线；
根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线；
根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
所述每口井反演结果均达到第一设定反演结果符合率，则根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果；
利用验证井检验所述所有井反演结果，所述验证井的反演结果达到第二设定反演结果符合率，则完成多井声波测井曲线一致性处理。
优选地，所述根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线的方法，包括：
分别根据所述所有井的声波测井曲线得到每口井的低频分量以及高频分量；
根据所述所有井的低频分量确定所述所有井的平均低频分量；
分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度；
分别根据所述所有井的声波测井曲线、所述高频幅度、所述高频分量平均幅度以及所述平均低频分量得到所述所有井的一致性处理后的声波测井曲线。
优选地，在所述根据所述所有井的低频分量确定所述所有井的平均低频分量，以及分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度之前，还包括：
确定所述所有井的低频分量取值范围；
若所述所有井中的某一口井的低频分量偏离其他井的低频分量取值范围的偏离量达到设定偏离量，则剔除此井的声波测井曲线。
优选地，在所述根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果之前，还包括：
获取所有井的实际测量地震道；
利用所述一致性处理后的声波测井曲线得到测井曲线合成的地震道；
根据所述测井曲线合成的地震道以及所述实际测量地震道的相关系数判断所述一致性处理后的声波测井曲线是否合格；
若合格，则根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
若不合格，则进行每口井的环境校正，得到校正后的声波测井曲线，根据所述校正后的声波测井曲线得到一致性处理后的声波测井曲线。
优选地，所述进行每口井的环境校正，得到校正后的声波测井曲线的方法，包括：
分别获取每口井的多个声波测井曲线；
分别将所述多个声波测井曲线拟合成一条所述声波测井曲线，得到所述所有井的声波测井曲线。
第二方面，本发明提供一种多井声波测井曲线一致性处理装置，包括：
获取单元，用于获取所有井的声波测井曲线；
生成单元，用于根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线；
单井反演单元，用于根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
多井反演单元，判断所述每口井反演结果均是否达到第一设定反演结果符合率，若达到所述第一设定反演结果符合率，则根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果；
检验单元，用于利用验证井检验多井反演结果，若所述验证井的反演结果达到第二设定反演结果符合率，则完成多井声波测井曲线一致性处理。
优选地，所述生成单元，包括：
分频单元，用于分别根据所述所有井的声波测井曲线得到每口井的低频分量以及高频分量；
平均低频分量确定单元，用于根据所述所有井的低频分量确定所述所有井的平均低频分量；
高频分量平均幅度确定单元，用于分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度；
声波测井曲线生产单元，用于分别根据所述所有井的声波测井曲线、所述高频幅度、所述高频分量平均幅度以及所述平均低频分量得到所述所有井的一致性处理后的声波测井曲线。
优选地，所述的装置，还包括：剔除单元；
所述剔除单元，用于确定所述所有井的低频分量取值范围，若所述所有井中的某一口井的低频分量偏离其他井的低频分量取值范围的偏离量达到设定偏离量，则剔除此井的声波测井曲线；
所述获取单元，获取经过所述剔除单元得到的所有井的声波测井曲线。
优选地，所述生成单元，还包括：判断单元；
所述判断单元，分别与所述生成单元以及所述单井反演单元连接，包括：
实际测量地震道获取单元，用于获取所有井的实际测量地震道；
测井曲线合成的地震道生成单元，用于利用所述一致性处理后的声波测井曲线得到测井曲线合成的地震道；
声波测井曲线判断单元，用于根据所述测井曲线合成的地震道以及所述实际测量地震道的相关系数判断所述一致性处理后的声波测井曲线是否合格；
若合格，则根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
若不合格，则进行每口井的环境校正，得到校正后的声波测井曲线，根据所述校正后的声波测井曲线得到一致性处理后的声波测井曲线。
第三方面，本发明提供一种存储介质，包括：处理器；
用于存储处理器可执行指令的存储器；
其中，所述处理器被配置为调用所述存储器存储的指令，以执行上述的方法。
本发明至少具有如下有益效果：
本发明提供一种多井声波测井曲线一致性处理方法、装置及存储介质，以解决多井测井曲线的一致性差，不同井约束的反演剖面差异较大，多井约束反演及多井建模结果失真的问题，消除非地质因素对测井曲线数据的影响，使得测井资料和地震资料相关性更高。
附图说明
通过以下参考附图对本发明实施例的描述，本发明的上述以及其它目的、特征和优点更为清楚，在附图中：
图1是本发明实施例一种多井声波测井曲线一致性处理方法的流程示意图；
图2是本发明实施例声波测井曲线高、低频分解过程示意图；
图3是本发明实施例多井声波测井曲线高低频分解结果图；
图4是本发明实施例多井声波测井曲线一致性处理后结果图；
图5是本发明实施例的合成地震记录对比示意图；其中，图5(a)为应用本发明前的合成地震记录示意图，图5(b)为应用本发明后的合成地震记录示意图；
图6是本发明实施例本发明实施例的地震反演对比示意图；其中，图6(a)为应用本发明前的地震反演示意图，图6(b)为应用本发明后的地震反演示意图。
具体实施方式
以下基于实施例对本发明进行描述，但是值得说明的是，本发明并不限于这些实施例。在下文对本发明的细节描述中，详尽描述了一些特定的细节部分。然而，对于没有详尽描述的部分，本领域技术人员也可以完全理解本发明。
此外，本领域普通技术人员应当理解，所提供的附图只是为了说明本发明的目的、特征和优点，附图并不是实际按照比例绘制的。
同时，除非上下文明确要求，否则整个说明书和权利要求书中的“包括”、“包含”等类似词语应当解释为包含的含义而不是排他或穷举的含义；也就是说，是“包含但不限于”的含义。
本公开实施例提供的一种多井声波测井曲线一致性处理方法的执行主体可以为任意的图像处理装置，例如多井声波测井曲线一致性处理方法可以由终端设备或服务器执行，其中，终端设备可以为用户设备(User Equipment，UE)、移动设备、用户终端、终端、蜂窝电话、无绳电话、个人数字处理(Personal Digital Assistant，PDA)、手持设备、计算设备、车载设备、可穿戴设备等。服务器可以为本地服务器或者云端服务器。在一些可能的实现方式中，该多井声波测井曲线一致性处理方法可以通过处理器调用存储器中存储的计算机可读指令的方式来实现。
现有商业软件多井声波曲线一致性处理方法要求的标准井和标准层，每一口约束井都具有不可替代的作用，对于测井曲线质量的判别是以合成地震记录最佳及反演效果最佳为判断准则。
本方法中，多井声波测井曲线一致性处理基本原理是：测井曲线的低频部分代表了相对大的沉积旋回，而测井曲线的高频部分代表了相对小的沉积旋回和岩性相对变化，对于坳陷沉积盆地，尤其是河流相沉积砂泥岩沉积环境条件下，在同一地震工区内，大的沉积旋回一般相同，也就是长周期的沉降环境相同，因此工区内所有井的低频分量应相近或相同，高频分量的幅度代表了不同岩性的速度差异，而不同井同一层位的声波测井曲线，尽管岩性不同，但泥岩的平均速度和砂岩的平均的速度及平均幅度差基本一致，所以我们通过多井声波测井曲线的高低频分解技术，将声波测井曲线分解，对高低频分量分别处理后再将其合并，获得多井一致性处理后的声波测井曲线。
图1是本发明实施例一种多井声波测井曲线一致性处理方法的流程示意图。如图1所示，一种多井声波测井曲线一致性处理方法，包括：步骤S101：获取所有井的声波测井曲线；步骤S102：根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线；步骤S103：根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；步骤S104：所述每口井反演结果均达到第一设定反演结果符合率，则根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果；步骤S105：利用验证井检验所述所有井反演结果，所述验证井的反演结果达到第二设定反演结果符合率，则完成多井声波测井曲线一致性处理。
若所述验证井的反演结果不能达到第二设定反演结果符合率，则从新选取所有井的声波测井曲线，继续S102-步骤S105的步骤，直到所述验证井的反演结果达到第二设定反演结果符合率为止。
本发明提供一种多井声波测井曲线一致性处理方法，以解决多井测井曲线的一致性差，不同井约束的反演剖面差异较大，多井约束反演及多井建模结果失真的问题，消除非地质因素对测井曲线数据的影响，使得测井资料和地震资料相关性更高。
步骤S101：获取所有井的声波测井曲线。
在本发明的实施例中，声波测井曲线为通过测井仪器在井下测量得到，可测井地面系统上得到需要的测井数据，可以用来计算储层的各种参数。
在步骤S101：获取所有井的声波测井曲线之前，可以对所有井的声波测井曲线进行预处理，完成单井测井曲线的环境影响校正。
在本发明的实施例中，测井公司或者现有商业测井解释软件，都可以完成单井测井曲线的环境影响校正，主要完成井径、泥浆浸入带、泥浆密度与矿化度、地层温度与压力、围岩以及仪器外径、间隙等非地质因素造成的误差。
例如，由泥岩段井壁垮塌造成的声波时差值增大的异常、在测井的顶底界面处产生的测井曲线的畸变、两次测井曲线进行拼接的时候出现的异常以及由于膏盐层段的存在，使得声波曲线缺失或出现声波异常值等，经过这些校正以后，以提高合成记录制作的精度为判别准则，或以同一口井，多种测井曲线一致解释效果最好，单井测井曲线环境校正完成后，才开始做多井声波曲线一致性处理。
步骤S102：根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线。
在本发明中，所述根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线的方法，包括：分别根据所述所有井的声波测井曲线得到每口井的低频分量以及高频分量；根据所述所有井的低频分量确定所述所有井的平均低频分量；分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度；分别根据所述所有井的声波测井曲线、所述高频幅度、所述高频分量平均幅度以及所述平均低频分量得到所述所有井的一致性处理后的声波测井曲线。
其中，每口井的低频分量以及高频分量为时间域的分量，即时间域的低频分量以及高频分量。
在本发明的实施例中，所述分别根据所述所有井的声波测井曲线得到每口井的低频分量以及高频分量的方法，包括：将所述所有井的声波测井曲线进行傅里叶变换，得到所述所有井的声波测井曲线的频域数据；将所述频域数据分解成低频频域分量以及高频频域分量；分别对低频频域分量以及高频频域分量做傅里叶逆变换，得到时间域的低频分量以及高频分量。其中，低频频域的范围可为3Hz-5Hz，高频频域的范围可为5Hz-400Hz。
在本发明的实施例中，对工区内所有井的声波测井曲线xij(深度域函数)做傅里叶变换，变换到频率域，如采用如下公式：
Figure BDA0002424454640000061
公式中，xij是声波测井曲线(深度域函数)，i代表井序号，j对应测井曲线深度样点，Xi(ω)是第i井的傅里叶变换结果(声波测井曲线的频域数据)，ω对应频率采样点，FFT[]为傅里叶变换。
在本发明的实施例中，将第i井的傅里叶变换结果(声波测井曲线的频域数据)Xi(ω)分解成低频频域分量以及高频频域分量，如：
Xi(ω)＝Xi低(ω)+Xi高(ω)。
公式中，低频频域分量Xi低(ω)，高频频域分量Xi高(ω)，低频的范围可为3Hz-5Hz，高频的范围可为5Hz-400Hz。
在本发明的实施例中，对低频频域分量和高频频域分量别做反傅里叶变换，变换回深度域(时间域)，完成了声波测井曲线的高低频分解，如：
Figure BDA0002424454640000062
Figure BDA0002424454640000063
公式中，时间域的低频分量xij低，高频分量xij高，FFT-1[]为傅里叶变换。
在所述根据所述所有井的低频分量确定所述所有井的平均低频分量，以及分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度之前，还包括：确定所述所有井的低频分量取值范围；若所述所有井中的某一口井的低频分量偏离其他井的低频分量取值范围的偏离量达到设定偏离量，则剔除此井的声波测井曲线。
在本发明的实施例中，若所述所有井中的某一口井的低频分量偏离其他井的低频分量取值范围的偏离量达到设定偏离量，则剔除此井的声波测井曲线的具体方法，包括：确定低频分量取值范围的最大值和最小值，若所述某一口井的低频分量偏离其他井的低频分量取值范围的最大值或最小值的设定偏离量，则剔除此井的声波测井曲线。
在本发明的实施例中，分析各井低频分量(时间域)，如果某井的低频分量偏离其他井的低频分量取值范围的20％(设定偏离量)以上，剔除该井，然后确定剔除该井后的所有井的平均低频分量以及根据剔除该井后所述每口井的高频分量确定每口井的针对目标层段的高频幅度。也就是说，低频分量的取值范围，包括低频分量的最大值和低频分量的最小值；例如，低频分量的最大值可为30Hz，低频分量的最小值可为3Hz-5Hz。如果某井的低频分量偏离其他井的低频分量取值范围的20％(设定偏离量)以上或者小于低频分量的最小值的20％，则剔除此井的声波测井曲线。
其中，根据所述所有井的低频分量确定所述所有井的平均低频分量，采用下公式计算。
Figure BDA0002424454640000071
公式中，平均低频分量为
Figure BDA0002424454640000072
i代表井序号，j对应测井曲线深度样点，n为工区内参与声波测井曲线一致性处理的井的个数，即所有井的数量。但是，如果剔除此井的声波测井曲线或剔除井的个数为1，则n为所有井的数量-1。
其中，分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度。
计算每口井的针对目标层段的高频幅度，采用下公式计算。
Figure BDA0002424454640000073
公式中，每口井的针对目标层段的高频幅度为σi高，j代表对应测井曲线深度样点，m为针对目标层段声波测井曲线样点个数。
其中，根据所述所有井的高频幅度确定所有井的高频分量平均幅度。
计算所有井的高频分量平均幅度，采用下公式计算。
Figure BDA0002424454640000081
公式中，n为工区内参与声波测井曲线一致性处理的井的个数，即所有井的数量。但是，如果剔除此井的声波测井曲线或剔除井的个数为1，则n为所有井的数量-1。
其中，分别根据所述所有井的声波测井曲线、所述高频幅度、所述高频分量平均幅度以及所述平均低频分量得到所述所有井的一致性处理后的声波测井曲线，采用下公式计算。用经过平均幅度调整后的高频和统一的低频分量相加求得一致性处理后的声波测井曲线。
Figure BDA0002424454640000082
公式中，x′ij为经过多井一致性处理后的声波测井曲线，即所有井的一致性处理后的声波测井曲线。
步骤S103：根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果。
在本发明的实施例中，反演是指波阻抗反演，反演的具体过程为：通过Z反演、Jason等具有地震反演功能的市场化软件。Z反演、Jason等具有地震反演功能的市场化软件可以从每口井的地震记录中估算地层反射系数，地层反射系数结合所述一致性处理后的声波测井曲线进行运算得到波阻抗信息(即，每口井的反演结果)，使其能够直观地与实际测井曲线计算出的波阻抗值进行对比。
步骤S104：所述每口井反演结果均达到第一设定反演结果符合率，则根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果。
在本发明的实施例中，第一设定反演结果符合率可为80％。所述每口井反演结果与该井实际测量值进行对比值一致的层段占该井总层段的80％(第一设定反演结果符合率)算达到标准。
在本发明的实施例中，根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果。即，Z反演、Jason等具有地震反演功能的市场化软件可以从所有井的地震记录中估算地层反射系数，地层反射系数结合所述一致性处理后的声波测井曲线进行运算得到波阻抗信息(即，所有井的反演结果)。
步骤S105：利用验证井检验所述所有井反演结果，所述验证井的反演结果达到第二设定反演结果符合率，则完成多井声波测井曲线一致性处理。
在本发明的实施例中，验证井不在前面参与运算的井，第二设定反演结果符合率可为80。
在本发明的实施例中，Z反演、Jason等具有地震反演功能的市场化软件可以从验证井的地震记录中估算地层反射系数，地层反射系数结合所述一致性处理后的声波测井曲线进行运算得到波阻抗信息(即，验证井的反演结果)。
在本发明的实施例中，验证井对反演结果验证：用验证井所在目的层的反演结果与该井实际测量值进行对比值一致的层段占该井总层段的80％(第二设定反演结果符合率)算达到标准。
在本发明中，在所述根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果之前，还包括：获取所有井的实际测量地震道；利用所述一致性处理后的声波测井曲线得到测井曲线合成的地震道；
根据所述测井曲线合成的地震道以及所述实际测量地震道的相关系数判断所述一致性处理后的声波测井曲线是否合格；
若合格，则根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
若不合格，则进行每口井的环境校正，得到校正后的声波测井曲线，根据所述校正后的声波测井曲线得到一致性处理后的声波测井曲线。
在本发明的实施例中，所述测井曲线合成的地震道，即用所述一致性处理后的声波测井曲线做合成地震记录检验。合成地震记录：由所述一致性处理后的声波测井曲线和测井密度曲线求出波阻抗曲线，再由波阻抗曲线转换为反射系数曲线，与地震子波褶积运算得到合成地震道，该方法称为合成地震记录。
地震记录的检验：由所述一致性处理后的声波测井曲线合成的地震道与实际测量地震道做相关性对比，检验两者的相似程度，得到相关系数，然后根据所述测井曲线合成的地震道以及所述实际测量地震道的相关系数判断所述一致性处理后的声波测井曲线是否合格。
在本发明的实施例中，所述据所述测井曲线合成的地震道以及所述实际测量地震道的相关系数判断所述一致性处理后的声波测井曲线是否合格的方法，包括：确定设定相关系数；所述测井曲线合成的地震道以及所述实际测量地震道的相关系数大于所述设定相关系数，则所述一致性处理后的声波测井曲线合格，然后根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果。其中，相关系数可为75％。的实施例中，用一致性处理后的声波曲线做合成地震记录检验，相关系数达到75％，处理后的声波测井曲线视为合格；相关系数未达到75％，一致性处理后的声波测井曲线视为不合格，则进行每口井的环境校正，得到校正后的声波测井曲线，根据所述校正后的声波测井曲线得到一致性处理后的声波测井曲线，即：分别根据所述校正后所有井的声波测井曲线得到每口井的低频分量以及高频分量；根据所述所有井的低频分量确定所述所有井的平均低频分量；分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度；分别根据所述所有井的声波测井曲线、所述高频幅度、所述高频分量平均幅度以及所述平均低频分量得到所述所有井的一致性处理后的声波测井曲线。具体可见上述的详细描述。
在本发明中，所述进行每口井的环境校正，得到校正后的声波测井曲线的方法，包括：分别获取每口井的多个声波测井曲线；分别将所述多个声波测井曲线拟合成一条所述声波测井曲线，得到所述所有井的声波测井曲线。
例如，每口井的单井环境校正的方法为一口井的多个声波测井曲线拟合成成一条声波测井曲线。N个井，即得到N条声波测井曲线。步骤S105中，若所述验证井的反演结果不能达到第二设定反演结果符合率，则可进行每口井的环境校正，得到校正后的声波测井曲线，继续S102-步骤S105的步骤，直到所述验证井的反演结果达到第二设定反演结果符合率为止。
本发明还提出了一种多井声波测井曲线一致性处理装置，包括：获取单元，用于获取所有井的声波测井曲线；生成单元，用于根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线；单井反演单元，用于根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；多井反演单元，判断所述每口井反演结果均是否达到第一设定反演结果符合率，若达到所述第一设定反演结果符合率，则根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果；检验单元，用于利用验证井检验多井反演结果，若所述验证井的反演结果达到第二设定反演结果符合率，则完成多井声波测井曲线一致性处理。具体可详见一种多井声波测井曲线一致性处理方法中的详细说明。
在本发明中，所述生成单元，包括：分频单元，用于分别根据所述所有井的声波测井曲线得到每口井的低频分量以及高频分量；平均低频分量确定单元，用于根据所述所有井的低频分量确定所述所有井的平均低频分量；高频分量平均幅度确定单元，用于分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度；声波测井曲线生产单元，用于分别根据所述所有井的声波测井曲线、所述高频幅度、所述高频分量平均幅度以及所述平均低频分量得到所述所有井的一致性处理后的声波测井曲线。
在本发明中，所述的装置，还包括：剔除单元；所述剔除单元，用于确定所述所有井的低频分量取值范围，若所述所有井中的某一口井的低频分量偏离其他井的低频分量取值范围的偏离量达到设定偏离量，则剔除此井的声波测井曲线；所述获取单元，获取经过所述剔除单元得到的所有井的声波测井曲线。
在本发明中，所述生成单元，还包括：判断单元；所述判断单元，分别与所述生成单元以及所述单井反演单元连接，包括：实际测量地震道获取单元，用于获取所有井的实际测量地震道；测井曲线合成的地震道生成单元，用于利用所述一致性处理后的声波测井曲线得到测井曲线合成的地震道；声波测井曲线判断单元，用于根据所述测井曲线合成的地震道以及所述实际测量地震道的相关系数判断所述一致性处理后的声波测井曲线是否合格；若合格，则根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；若不合格，则进行每口井的环境校正，得到校正后的声波测井曲线，根据所述校正后的声波测井曲线得到一致性处理后的声波测井曲线。
发明还提出了一种存储介质，包括：处理器；用于存储处理器可执行指令的存储器；其中，所述处理器被配置为调用所述存储器存储的指令以执行上述的方法。具体可详见一种多井声波测井曲线一致性处理方法中的详细说明。
图2是本发明实施例声波测井曲线高、低频分解过程示意图。如图2所示，1(实线)为原始声波测井曲线，即所有井的声波测井曲线。2(右图左侧)为声波测井曲线的高频部分，3(右图右侧)为声波测井曲线的低频部分，4(左图虚线)为声波测井曲线的高频部分和低频部分相加结果，表明高、低频分解后再合并后完全与原声波曲线一致。(注：为了把原始声波测井曲线与合成的声波测井曲线绘在一起，又能看出是两条曲线，绘图时人为将合成声波做了一个微小移位，否则两条完全覆盖，无法区分)。
图3是本发明实施例多井声波测井曲线高低频分解结果图。如图3所示，可见各井的低频分量不一致(如图右上方所示)，高频幅度分布范围比较大(如图右下方所示)图4是本发明实施例多井声波测井曲线一致性处理后结果图。如图4所示，可见高频幅度分布范围比较集中。
图5是本发明实施例的合成地震记录对比示意图；其中，图5(a)为应用本发明前的合成地震记录示意图，图5(b)为应用本发明后的合成地震记录示意图。如图5所示，图5(a)为本方法处理前合成地震记录，与井旁道相似系数0.795，图5(b)为本方法处理后合成地震记录，与井旁道相似系数0.859，与本方法处理前合成地震记录与井旁道相似系数有所提高。
图6是本发明实施例本发明实施例的地震反演对比示意图；其中，图6(a)为应用本发明前的地震反演示意图，图6(b)为应用本发明后的地震反演示意图。如图6所示，图6(a)为本方法处理前地震反演，图6(b)为本方法处理后地震反演结果，经本方法处理后，反演结果主要目标层(含油)更加清楚(见椭圆标注区域)
以上所述实施例仅为表达本发明的实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形、同等替换、改进等，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种多井声波测井曲线一致性处理方法，其特征在于，包括：
获取所有井的声波测井曲线；
根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线；
根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
所述每口井反演结果均达到第一设定反演结果符合率，则根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果；
利用验证井检验所述所有井反演结果，所述验证井的反演结果达到第二设定反演结果符合率，则完成多井声波测井曲线一致性处理。
2.根据权利要求1所述的方法，其特征在于，所述根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线的方法，包括：
分别根据所述所有井的声波测井曲线得到每口井的低频分量以及高频分量；
根据所述所有井的低频分量确定所述所有井的平均低频分量 ；
分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度；
分别根据所述所有井的声波测井曲线、所述高频幅度、所述高频分量平均幅度以及所述平均低频分量得到所述所有井的一致性处理后的声波测井曲线。
3.根据权利要求2所述的方法，其特征在于，在所述根据所述所有井的低频分量确定所述所有井的平均低频分量，以及分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度之前，还包括：
确定所述所有井的低频分量取值范围；
若所述所有井中的某一口井的低频分量偏离其他井的低频分量取值范围的偏离量达到设定偏离量，则剔除此井的声波测井曲线。
4.根据权利要求1-3任一项所述的方法，其特征在于，在所述根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果之前，还包括：
获取所有井的实际测量地震道；
利用所述一致性处理后的声波测井曲线得到测井曲线合成的地震道 ；
根据所述测井曲线合成的地震道以及所述实际测量地震道的相关系数判断所述一致性处理后的声波测井曲线是否合格；
若合格，则根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
若不合格，则进行每口井的环境校正，得到校正后的声波测井曲线，根据所述校正后的声波测井曲线得到一致性处理后的声波测井曲线。
5.根据权利要求4所述的方法，其特征在于，所述进行每口井的环境校正，得到校正后的声波测井曲线的方法，包括：
分别获取每口井的多个声波测井曲线；
分别将所述多个声波测井曲线拟合成一条所述声波测井曲线，得到所述所有井的声波测井曲线。
6.一种多井声波测井曲线一致性处理装置，其特征在于，包括：
获取单元，用于获取所有井的声波测井曲线；
生成单元，用于根据所述所有井的声波测井曲线得到一致性处理后的声波测井曲线；
单井反演单元，用于根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
多井反演单元，判断所述每口井反演结果均是否达到第一设定反演结果符合率，若达到所述第一设定反演结果符合率，则根据所述一致性处理后的声波测井曲线进行所有井的反演，得到所述所有井反演结果；
检验单元，用于利用验证井检验所述所有井反演结果，若所述验证井的反演结果达到第二设定反演结果符合率，则完成多井声波测井曲线一致性处理。
7.根据权利要求6所述的装置，其特征在于，所述生成单元，包括：
分频单元，用于分别根据所述所有井的声波测井曲线得到每口井的低频分量以及高频分量；
平均低频分量确定单元，用于根据所述所有井的低频分量确定所述所有井的平均低频分量 ；
高频分量平均幅度确定单元，用于分别根据所述每口井的高频分量确定每口井的针对目标层段的高频幅度，根据所述所有井的高频幅度确定所有井的高频分量平均幅度；
声波测井曲线生产单元，用于分别根据所述所有井的声波测井曲线、所述高频幅度、所述高频分量平均幅度以及所述平均低频分量得到所述所有井的一致性处理后的声波测井曲线。
8.根据权利要求7所述的装置，其特征在于，还包括：剔除单元；
所述剔除单元，用于确定所述所有井的低频分量取值范围，若所述所有井中的某一口井的低频分量偏离其他井的低频分量取值范围的偏离量达到设定偏离量，则剔除此井的声波测井曲线；
所述获取单元，获取经过所述剔除单元得到的所有井的声波测井曲线。
9.根据权利要求6-8任一项所述的装置，其特征在于，所述生成单元，还包括：判断单元；
所述判断单元，分别与所述生成单元以及所述单井反演单元连接，包括：
实际测量地震道获取单元，用于获取所有井的实际测量地震道；
测井曲线合成的地震道生成单元，用于利用所述一致性处理后的声波测井曲线得到测井曲线合成的地震道 ；
声波测井曲线判断单元，用于根据所述测井曲线合成的地震道以及所述实际测量地震道的相关系数判断所述一致性处理后的声波测井曲线是否合格；
若合格，则根据所述一致性处理后的声波测井曲线进行每口井的反演，得到所述每口井反演结果；
若不合格，则进行每口井的环境校正，得到校正后的声波测井曲线，根据所述校正后的声波测井曲线得到一致性处理后的声波测井曲线。
10.一种存储介质，其特征在于，包括：处理器；
用于存储处理器可执行指令的存储器；
其中，所述处理器被配置为调用所述存储器存储的指令，以执行权利要求1至5中任意一项所述的方法。
CN202010216043.9A 2020-03-25 2020-03-25 一种多井声波测井曲线一致性处理方法、装置及存储介质 Active CN111827966B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010216043.9A CN111827966B (zh) 2020-03-25 2020-03-25 一种多井声波测井曲线一致性处理方法、装置及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010216043.9A CN111827966B (zh) 2020-03-25 2020-03-25 一种多井声波测井曲线一致性处理方法、装置及存储介质

Publications (2)

Publication Number Publication Date
CN111827966A true CN111827966A (zh) 2020-10-27
CN111827966B CN111827966B (zh) 2022-04-15

Family

ID=72913532

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010216043.9A Active CN111827966B (zh) 2020-03-25 2020-03-25 一种多井声波测井曲线一致性处理方法、装置及存储介质

Country Status (1)

Country Link
CN (1) CN111827966B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4310887A (en) * 1972-08-28 1982-01-12 Schlumberger Technology Corporation Verification and calibration of well logs and reconstruction of logs
US5587966A (en) * 1994-10-13 1996-12-24 Schlumberger Technology Corporation Sonic well logging methods and apparatus for processing flexural wave in optimal frequency band
US20050261835A1 (en) * 2004-05-18 2005-11-24 Schlumberger Technology Corporation Sonic well logging methods and apparatus utilizing parametric inversion dispersive wave processing
US20070280050A1 (en) * 2006-06-05 2007-12-06 Stefano Volterrani Low frequency model estimation
CN102353985A (zh) * 2011-06-01 2012-02-15 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于非下采样Contourlet变换的拟声波曲线构建方法
CN103018778A (zh) * 2012-12-12 2013-04-03 中国石油天然气股份有限公司 校正声波测井曲线进行储层预测的方法及设备
CN103850679A (zh) * 2014-04-01 2014-06-11 北京师范大学 一种利用多种测井曲线对声波时差曲线进行重构的方法
CN104047597A (zh) * 2014-01-26 2014-09-17 中国石油大学(华东) 含油气泥页岩地层测井曲线标准化方法
CN105093280A (zh) * 2014-05-15 2015-11-25 中国石油化工股份有限公司 表层模型对地震数据影响的低频与高频成分的分解方法
CN107065011A (zh) * 2017-06-22 2017-08-18 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法
CN109343120A (zh) * 2018-10-17 2019-02-15 吉林大学 融入约束稀疏脉冲反演低频补偿的声波曲线重构方法
US20190169986A1 (en) * 2014-06-06 2019-06-06 Quantico Energy Solutions, Llc Real-time synthetic logging for optimization of drilling, steering, and stimulation

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4310887A (en) * 1972-08-28 1982-01-12 Schlumberger Technology Corporation Verification and calibration of well logs and reconstruction of logs
US5587966A (en) * 1994-10-13 1996-12-24 Schlumberger Technology Corporation Sonic well logging methods and apparatus for processing flexural wave in optimal frequency band
US20050261835A1 (en) * 2004-05-18 2005-11-24 Schlumberger Technology Corporation Sonic well logging methods and apparatus utilizing parametric inversion dispersive wave processing
US20070280050A1 (en) * 2006-06-05 2007-12-06 Stefano Volterrani Low frequency model estimation
CN102353985A (zh) * 2011-06-01 2012-02-15 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于非下采样Contourlet变换的拟声波曲线构建方法
CN103018778A (zh) * 2012-12-12 2013-04-03 中国石油天然气股份有限公司 校正声波测井曲线进行储层预测的方法及设备
CN104047597A (zh) * 2014-01-26 2014-09-17 中国石油大学(华东) 含油气泥页岩地层测井曲线标准化方法
CN103850679A (zh) * 2014-04-01 2014-06-11 北京师范大学 一种利用多种测井曲线对声波时差曲线进行重构的方法
CN105093280A (zh) * 2014-05-15 2015-11-25 中国石油化工股份有限公司 表层模型对地震数据影响的低频与高频成分的分解方法
US20190169986A1 (en) * 2014-06-06 2019-06-06 Quantico Energy Solutions, Llc Real-time synthetic logging for optimization of drilling, steering, and stimulation
CN107065011A (zh) * 2017-06-22 2017-08-18 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法
CN109343120A (zh) * 2018-10-17 2019-02-15 吉林大学 融入约束稀疏脉冲反演低频补偿的声波曲线重构方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
MOHAMAD等: "Seismic velocity deviation log: An effective method for evaluating spatial distribution of reservoir pore types", 《JOURNAL OF APPLIED GEOPHYSICS》 *
姜传金等: "拟声波曲线构建的意义及应用", 《大庆石油地质与开发》 *
宋维琪等: "测井声波时差反演重构技术研究及应用", 《地震地质》 *
李博南: "基于岩石物理分析的碳酸盐岩储层声波速度测井校正方法", 《科学技术与工程》 *
梁宇: "密井网区井震匹配性改进措施及应用", 《北京2018国际地球物理会议暨展览电子论文集》 *
殷积峰等: "合成声波测井及在储层研究中的应用", 《石油学报》 *
辛连发: "合成声波测井剖面在阿尔善油田储层横向预测中的应用", 《石油地球物理勘探》 *
陈志德等: "长垣油田密度曲线的重构与应用", 《大庆石油地质与开发》 *

Also Published As

Publication number Publication date
CN111827966B (zh) 2022-04-15

Similar Documents

Publication Publication Date Title
CN109425896B (zh) 白云岩油气储层分布预测方法及装置
CN106405651B (zh) 一种基于测井匹配的全波形反演初始速度模型构建方法
US11808143B2 (en) Methods and apparatus to measure formation features
CN111980676B (zh) 一种阵列声波测井评价固井质量的方法及处理装置
US11230922B2 (en) Fracture interpretation with resistivity and sonic logs in biaxial anisotropic formations
CN105588883B (zh) 三维岩石力学参数获取方法和系统
US7447114B2 (en) Non-linear seismic trace matching to well logs
WO2020215170A1 (zh) 基于地震岩石物理实验分析的测井与地震速度匹配方法
CN109143345B (zh) 基于模拟退火的品质因子q非线性反演方法及系统
CN111827966B (zh) 一种多井声波测井曲线一致性处理方法、装置及存储介质
CN111273349B (zh) 一种用于海底浅部沉积层的横波速度提取方法及处理终端
CN110007349B (zh) 一种弹性参数反演方法
CN111435174B (zh) 强反射地区地震资料振幅补偿方法及装置
CN115453629A (zh) 一种基于时频域自适应滤波的波形校正方法及装置
CN112180464B (zh) 一种储层物性的识别方法
CN113050191B (zh) 基于双参数的页岩油toc预测方法及装置
CN111458751B (zh) 一种综合时深标定方法及装置
CN111352157B (zh) 一种横波静校正方法及系统
CN109061737B (zh) 一种基于合成地震记录的储层预测方法及装置
CN109884701B (zh) 地质体散射角度导向深度成像方法
CN105178950B (zh) 一种确定细粒沉积岩中总有机碳的方法
CN113009579B (zh) 地震数据反演方法及装置
CN113219530B (zh) 非稳态盲反褶积方法及装置
CN114721044B (zh) 一种多频率接收函数和振幅比联合反演地壳结构的方法及系统
CN114114406B (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
TA01 Transfer of patent application right

Effective date of registration: 20220107

Address after: 163453 Heilongjiang Province, Daqing City Ranghulu District No. 233 South Central Avenue

Applicant after: Daqing Oilfield Co.,Ltd.

Applicant after: PetroChina Company Limited

Address before: 163453 Heilongjiang Province, Daqing City Ranghulu District No. 233 South Central Avenue

Applicant before: Daqing Oilfield Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant