CN113376696A - 基于多子波的井震标定方法及装置 - Google Patents
基于多子波的井震标定方法及装置 Download PDFInfo
- Publication number
- CN113376696A CN113376696A CN202010115378.1A CN202010115378A CN113376696A CN 113376696 A CN113376696 A CN 113376696A CN 202010115378 A CN202010115378 A CN 202010115378A CN 113376696 A CN113376696 A CN 113376696A
- Authority
- CN
- China
- Prior art keywords
- seismic
- well
- wavelets
- data
- different
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000012937 correction Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 230000002159 abnormal effect Effects 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 21
- 230000008569 process Effects 0.000 description 10
- 238000010183 spectrum analysis Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect 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
- 230000000704 physical effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
- G01V1/366—Seismic filtering by correlation of seismic signals
-
- 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
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于多子波的井震标定方法及装置，该方法包括：获取目标井位置处采集的测井数据和井旁地震道数据；根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗；在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录；对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。本发明采用不同频率不同相位的地震子波进行井震标定，能够使得生成的合成地震记录与井旁地震道数据具有良好的匹配，从而实现更准确的层位标定。
Description
技术领域
本发明涉及石油勘探领域，尤其涉及一种基于多子波的井震标定方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
地震资料解释是地球物理勘探的一项重要工作，它利用地震勘探获得的二维或者三维地震资料进行分析，研究和确定地下地质特征，将地震资料中包含的地震信息转化为地质成果。具体地，根据地震资料中包含的有效波的运动学特征(旅行时、速度等)和动力学特征(振幅、频率、相位、波形等)，结合地质、钻井、测井等资料研究地下地质构造形态，推测地层岩性、厚度及地层接触关系，进行储层、油藏描述及综合解释，绘制相关成果图件确定各种类型圈闭并评价其含油气性，指出有利含油气区带，提供钻探井位和勘探部署意见。
地震-地质层位标定是地震资料解释过程中一项非常重要的研究内容，在收集所有测井数据(主要是声波时差曲线和密度曲线)后，首先对测井曲线进行预处理(包括环境校正及质量检查)，再对测井数据进行基准面校正(即从井位地面到地震基准面的高程校正)。通常，由声波测井资料和密度测井资料计算出地层界面的反射系数序列，将反射系数与提取的地震子波进行褶积得到初始合成地震记录。在初次完成合成地震记录后，将录井和测井资料解释的地质分层数据标定到合成地震记录上，再采用手动或自动方式将合成地震记录与井旁地震道进行波组和波形对比，找到相似系数最大时两者对应的同相轴，根据地质层位相当的地震层位标定在井旁地震道上，再根据地质层位深度和地震层位时间计算出平均速度。
需要注意的是，当人工合成地震记录与井旁地震道相关性不好的情况下，需要修改速度或者调整子波参数，通常在制作合成地震记录的过程中使用常相位单子波，子波可以是理论子波(如雷克子波、双余弦子波等)或者从井旁地震道数据中提取的实际地震子波，再调整速度从而得到和井旁地震道相关性较好的人工合成地震记录。
然而，对于一些地震资料频率、相位差别大的地区，使用上述方法很难达到人工合成地震记录和井旁道地震道匹配良好，这会导致地质层位标定到时间域地震剖面上产生一定的误差，从而对后期的地震层位追踪解释、成图等造成影响。
针对上述问题，目前尚未提出有效的解决方案。
发明内容
本发明实施例提供一种基于多子波的井震标定方法，用以解决现有技术中，由于实测地震道数据纵向分辨率不同，使得采用单一频率的地震子波对测井数据进行井震标定得到的合成地震记录，与井旁地震道数据匹配结果较差，影响到层位标定准确性的技术问题，该方法包括：获取目标井位置处采集的测井数据和井旁地震道数据；根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗；在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录；对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
本发明实施例还提供一种基于多子波的井震标定装置，用以解决现有技术中，由于实测地震道数据纵向分辨率不同，使得采用单一频率的地震子波对测井数据进行井震标定得到的合成地震记录，与井旁地震道数据匹配结果较差，影响到层位标定准确性的技术问题，该装置包括：数据集采集模块，用于获取目标井位置处采集的测井数据和井旁地震道数据；地震记录时窗划分模块，用于根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗；合成地震记录生成模块，用于在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录；合成地震记录修正模块，用于对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
本发明实施例还提供一种计算机设备，用以解决现有技术中，由于实测地震道数据纵向分辨率不同，使得采用单一频率的地震子波对测井数据进行井震标定得到的合成地震记录，与井旁地震道数据匹配结果较差，影响到层位标定准确性的技术问题，该计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述基于多子波的井震标定方法。
本发明实施例还提供一种计算机可读存储介质，用以解决现有技术中，由于实测地震道数据纵向分辨率不同，使得采用单一频率的地震子波对测井数据进行井震标定得到的合成地震记录，与井旁地震道数据匹配结果较差，影响到层位标定准确性的技术问题，该计算机可读存储介质存储有执行上述基于多子波的井震标定方法的计算机程序。
本发明实施例中，在获取到目标井位置处采集的测井数据和井旁地震道数据后，根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗，并在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录，进而对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
通过本发明实施例，采用不同频率不同相位的地震子波进行井震标定，能够使得生成的合成地震记录与井旁地震道数据具有良好的匹配，从而实现更准确的层位标定。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中提供的一种基于多子波的井震标定方法流程图；
图2为本发明实施例中提供的一种采用不同频率不同相位多子波对测井数据进行井震标定的方法流程图；
图3为本发明实施例中提供的对某一目标井的井旁地震道数据进行频谱分析的结果示意图；
图4为本发明实施例中提供的一种采用雷克子波进行井震标定的结果示意图；
图5为本发明实施例中提供的又一种采用雷克子波进行井震标定的结果示意图；
图6为本发明实施例中提供的在不同时窗内采用不同频率的子波对测井数据进行井震标定的结果示意图；
图7a为本发明实施例中提供的时窗Δt1对应的子波示意图；
图7b为本发明实施例中提供的时窗Δt2对应的子波示意图；
图7c为本发明实施例中提供的时窗Δt3对应的子波示意图；
图7d为本发明实施例中提供的时窗Δt1、Δt2和Δt3对应的子波比对结果示意图；
图8为本发明实施例中提供的在不同时窗内采用可变相位的子波对测井数据进行井震标定的结果；
图9为本发明实施例中提供的一种基于多子波的井震标定装置示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施，其中的步骤顺序不作限定，可根据需要作适当调整。
本发明实施例中提供了一种基于多子波的井震标定方法，图1为本发明实施例中提供的一种基于多子波的井震标定方法流程图，如图1所示，该方法包括如下步骤：
S101，获取目标井位置处采集的测井数据和井旁地震道数据。
需要说明的是，地震资料解释是根据地震剖面上的反射同相轴对地质界面的空间分布或地层构造进行描述，因而，准确地进行层位标定以建立地震反射同相轴与地质层位的对应关系对于地震资料解释十分重要。通过测井工作(即沿钻井轨迹，自浅至深对井壁附近地层的岩性和物性参数进行测量和记录)得到的测井资料中包含井点位置的地层特征和岩性特征，通过建立测井信息与地震信息之间的对应关系，即可将测井得到的地层信息对应到地震剖面上，从而实现对地震剖面上反射信号的层位标定。对于某一目标井，在井点位置，既有实际地震记录，也有实际测井数据，对测井数据进行井震标定，得到合成地震记录，将其与实际地震记录进行比对，即可建立地质层位与地震反射之间的对应关系。
本发明实施例上述S101中获取的测井数据为某一目标井的实际测井数据，井旁地震道数据可以是该目标井附近预设道数的地震道检测到的实际地震数据。
可选地，本发明实施中涉及的地震数据可以是二维地震数据，也可以是三维地震数据。
S102，根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗。
需要说明的是，无论是二维地震数据，还是三维地震数据，从地表到目的层位再到更深的地质层位，由于各种采集因素或地层本身特性的影响不同，使得不同深度对应的实际地震数据的分辨率和频率是不同的，实际上，地震波在传播的过程中，高频成分被逐渐吸收，使得层位越深，地震记录的主频会降低。如果按照现有技术采用单一固定频率(即相同频率)的地震子波进行井震标定，则生成的合成地震记录往往很难与井旁地震道数据具有较好的相关性，从而使得基于二者比对结果实现的层位标定结果不准确。
由此，本发明实施例通过上述S102根据地震数据纵向分辨率的不同，将井旁地震道数据(即实际地震数据)划分为多个时窗，以便在不同的时窗内采用与井旁地震道主频一致的子波进行井震标定，使得生成的合成地震记录与井旁地震道数据具有良好的相关性，从而实现准确的层位标定。
S103，在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录。
需要说明的是，本发明实施例中采用不同频率子波对测井数据进行井震标定的方法与现有技术中采用单频率子波对测井数据进行井震标定的过程类似，不同点在于，在不同的时窗内采用的子波的频率不同，即由声波和密度测井曲线计算得到反射系数序列，进而将反射系数序列与各个时窗内不同频率的地震子波进行褶积得到合成地震记录。因而，一种可选的实施方式中，上述S103具体可以通过如下步骤来实现：根据测井数据中包含的声波曲线和密度曲线，计算反射系数序列；计算反射系数序列与不同时窗内不同频率的子波的褶积，生成测井数据对应的合成地震记录。
可选地，为了确保用于井震标定段的制作的人工合成记录不受测井曲线数据的影响，在根据测井数据中包含的声波曲线和密度曲线，计算反射系数序列之前，本发明实施例提供的基于多子波的井震标定方法还可以包括如下步骤：对测井数据中包含的声波曲线和密度曲线进行环境校正和异常值处理。
需要注意的是，本发明实施例中采用的不同频率的子波可以是不同频率的典型子波(例如，雷克子波)，也可以是从地震数据中提取出来的不同频率的子波。
一种实施例中，在计算反射系数序列与不同时窗内不同频率的子波的褶积，生成测井数据对应的合成地震记录的时候，可以提取不同时窗对应的不同主频的雷克子波，进而计算反射系数序列与提取的不同主频的雷克子波的褶积，生成测井数据对应的合成地震记录。
另一种实施例中，在计算反射系数序列与不同时窗内不同频率的子波的褶积，生成测井数据对应的合成地震记录的时候，可以从井旁地震道数据中，提取不同时窗对应的不同频率的地震子波，进而计算反射系数序列与提取的不同频率的地震子波的褶积，生成测井数据对应的合成地震记录。
S104，对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
需要说明的是，由于在相同频带下零相位信号比最小相位信号、最大相位信号与混合相位信号对地层有更高的分辨率，因而，现有地震资料解释中，会将子波信号进行零相位处理，即在不同的时窗内，均采用零相位的子波。然而，实际地震资料通常并非是零相位的，所以从地震数据中提取的子波也是有零相位、最小相位、最大相位和混合相位等多种类型，因而，本发明实施例通过上述S103在不同时窗采用不同频率的子波对测井数据进行井震标定生成合成地震记录后，通过上述S104对不同时窗内的子波进行相位调整，来使得生成的合成地震记录与井旁地震道数据满足预设相关性条件，能够使得生成的合成地震记录与井旁地震道数据能达到最大的相关性。由于本发明实施例在对测井数据进行井震标定生成合成地震记录的过程中，在不同地震纵向分辨率的时间段(即不同的时窗)内使用不同频率不同相位的子波(雷克子波或者由地震资料提取的非常数相位地震子波)，因而，通过本发明实施例生成的合成地震记录与井旁地震道数据能达到最大的相关性。
可选地，本发明实施例中预设相关性条件可以是合成地震记录与井旁地震道数据之间的相关系数高于预设阈值。
由上可知，本发明实施例提供的基于多子波的井震标定方法，在获取到目标井位置处采集的测井数据和井旁地震道数据后，根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗，并在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录，进而对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
通过本发明实施例提供的基于多子波的井震标定方法，采用不同频率不同相位的地震子波进行井震标定，能够使得生成的合成地震记录与井旁地震道数据具有良好的匹配，从而实现更准确的层位标定。
图2为本发明实施例中提供的一种采用不同频率不同相位多子波对测井数据进行井震标定的方法流程图，如图2所示，包括如下步骤：
①分析井旁地震道附近地震资料的信噪比和分辨率，根据井旁地震道数据的分辨率，将地震数据从浅至深划分为多个连续的时窗，例如，n个连续的时窗。n的取值取决于地震数据的频率和地质层序界面标准层的数量，根据实际资料情况n≥2，需要注意的是，时窗不宜过多，但时窗最低应大于子波长度的3倍。
②对测井曲线进行环境校正和异常值处理，便于求取准确的反射系数。反射系数反映了上下地层之间的阻抗差，其表达式为：
R＝(P2-P1)/(P2+P1)＝(ρ2v2-ρ1v1)/(ρ2v2+ρ1v1) (1)
其中，R表示反射系数；P表示纵波阻抗；ρ表示密度；v表示速度。
由公式(1)可知，要得到精确的反射系数，必须要有准确的速度和密度，速度可通过声波曲线求取倒数获得。在测井过程中井壁垮塌和井壁噪声是影响声波曲线质量的重要因素，部分井声波曲线存在周波跳跃假象，使测井资料失去了真实性和可靠性。通过数理统计的方法分析测井曲线去除环境影响以及曲线中的异常值，有效提高计算反射系数的准确性。
③分时窗提取不同主频的雷克子波，采用不同主频的雷克子波对测井数据进行井震标定，生成对应的合成地震记录。由于地震剖面从浅层到深层的频率变化较大，浅层频率较高，深层频率较低，如果只用一个固定频率的子波进行标定效果不理想。在生成初始合成地震记录的时候，可以先提取与地震资料主频相当的雷克子波，即对各个不同的时窗，提取不同主频的雷克子波(相邻时窗子波频率不同)。由于雷克子波只有一个正峰，两侧旁瓣延续时间短，能够很快收敛。
④如果采用不同主频的雷克子波进行井震标定生成的合成地震记录与井旁地震道数据不符合预设相似条件，则可使用井旁地震道数据和测井数据提取不同频率的子波来对测井数据进行井震标定以生成合成地震记录。
需要说明的是，合成地震记录是地层反射系数序列和子波的褶积再加上噪音，噪音通常认为是随机的，与地震道数据互不相关，噪音在采集过程中对地震数据影响较小，可以忽略。则合成地震记录表达式为：
F(t)＝R(t)×S(t)+N(t) (2)
其中，F(t)表示合成地震记录；R(t)表示反射系数序列；S(t)表示子波；N(t)表示噪音。
需要注意的是，本发明实施例提供的基于多子波的井震标定方法，在具体实施过程中，可以先采用不同主频的雷克子波进行井震标定，再采用从井旁地震道数据中提取的不同频率的地震子波进行井震标定。
例如，对于一个四层序界面的地质模型，可以先使用不同主频的雷克子波与反射系数序列的褶积得到合成记录标定区域内一级、二级层序界面后，再使用从地震资料提取的不同频率的地震子波对三级、四级层序界面进行准确标定。使用井旁地震道数据中提取的实际子波再次制作合成记录进行地震-地质层位的标定，能够在雷克子波制作合成记录标定的基础上，将三级、四级层序界面标定的更加准确。
⑤针对研究的目的层段，在标志层对准的前提下对合成记录进行局部微调，主要从波形特征、波阻特征等方面对比。对单个时窗内的子波进行相位的小幅度调整，增加或者减小子波的相位，其目的是改变了局部的速度，以改变合成地震记录(也称人工合成地震记录)，这样便能将人工合成地震记录与井旁地震道合成地震记录的相关性进一步提高。
以某一研究地区的A1井为例，在获取到A1井位置处的三维地震资料、测井资料和地质分层数据后，对A1井附近地震道的地震数据进行频谱分析。图3为本发明实施例中提供的对A1井的井旁地震道数据进行频谱分析的结果示意图，如图3所示，在A1井位置的左侧和右侧各取20道地震道数据，目的层以1920-2650m左右的地层为主要目的层，从上到下分别选取Δt1(1920ms～2100ms)、Δt2(2100ms～2400ms)和Δt3(2400ms～2650ms)三个时窗进行地震资料的频谱分析，从图3的频谱分析结果可以看到，时窗Δt1内地震资料的主频为40hz，振幅能量最弱，最大振幅为100000；时窗Δt2内地震资料的主频明显降低，约为30hz，振幅能量增强，最大振幅为240000；时窗Δt3内地震资料的主频增加到33hz，振幅能量位于中间，最大振幅为130000。
分析可知，不同深度段的地震数据会根据岩性、温度、压力、噪音等多种因素造成在采集后的成果地震数据在不同地方的能量、频率、相位等的不同。在本发明实施例图3所示的频谱分析结果，时窗Δt1内地震数据的主频最高，振幅能量最弱；时窗Δt2内地震数据的主频最低，振幅能量最高，在地震剖面上显示出强振幅。
对采集到的测井数据中声波曲线和密度曲线进行环境校正和异常值处理，以确保用于井震标定段的制作的人工合成记录不受测井曲线数据的影响，通过使用相应的测井解释软件对曲线进行预处理。A1井附近的地震道的1920ms～2650ms的时窗内提取主频为33hz的雷克子波进行初次标定，结合全区的区域地质资料，在测井曲线上确定2925米为一个标志层，同时也是岩性分界面，在地震数据上表现为一套强反射波峰的顶界，经过分析认为它是一个三级层序界面(1938ms)标志层，将这样的大界面作为井震标定的一级标志层。
图4所示为采用雷克子波进行井震标定一种结果，分析可知，如果在1920ms～2100ms时窗范围内的一组强反射波峰和人工合成记录标定准确后，则在2400ms～2650ms左右的另外一套强反射同相轴(标志层)未能准确进行地震地质分层的标定，此时在目的层段井震标定相关系数只有0.65左右。图5所示为采用雷克子波进行井震标定又一种结果，分析可知，如果将2450ms-2650ms左右的地震标志层标定准确后，则在1920ms-2100ms时窗范围内标定结果也会出现偏差，相关系数也只有0.55。
通过对地震资料分析可知，由于不同深度的地震资料的主频、分辨率不同，而在实际标定过程中采用单一主频、单一相位的子波会造成图4或图5的结果。
为解决上述问题，本发明实施例在三个时窗分别先提取主频为40hz、30hz、33hz的理论子波(雷克子波)，对时窗Δt2强反射的顶底进行标定，生成合成地震记录的原理和单子波方法类似，不同的是在不同的深度采用了不同频率的子波，相应地生成的波阻抗也不一致，不同时间段用不同子波计算的人工地震道也会不同，标定结果如图6所示，此时井震标定的相关系数达到0.7，然后从井旁地震道数据提取三个时窗对应的可变相位的子波进行精细标定，图7a所示为时窗Δt1对应的地震子波、图7b所示为时窗Δt2对应的地震子波、图7c所示为时窗Δt3对应的地震子波，图7d为三个时窗Δt1、Δt2和Δt3对应三种地震子波的比对结果，采用三种提取地震子波进行井震标定的标定结果如图8所示，此时井震标定的相关系数达到了0.85，比零相位的雷克子波的标定结果更好，说明实际地震资料本身是有一定的相位变化，本次使用的地震资料在1920ms～2650ms时窗内有-50～+50的相位变化。可见相对于单一频率子波的标定，本发明实施例应用分频可变相位多子波分段标定方法显著提高了井震标定的精度，将其运用在各种地震资料解释中，能够有效提高地质层位地震标定的效果。
基于同一发明构思，本发明实施例中还提供了一种基于多子波的井震标定装置，如下面的实施例所述。由于该装置实施例解决问题的原理与基于多子波的井震标定方法相似，因此该装置实施例的实施可以参见方法的实施，重复之处不再赘述。
图9为本发明实施例中提供的一种基于多子波的井震标定装置示意图，如图9所示，该装置可以包括：数据集采集模块91、地震记录时窗划分模块92、合成地震记录生成模块93和合成地震记录修正模块94。
其中，数据集采集模块91，用于获取目标井位置处采集的测井数据和井旁地震道数据；地震记录时窗划分模块92，用于根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗；合成地震记录生成模块93，用于在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录；合成地震记录修正模块94，用于对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
由上可知，本发明实施例提供的基于多子波的井震标定装置，通过数据集采集模块91获取目标井位置处采集的测井数据和井旁地震道数据；通过地震记录时窗划分模块92根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗；通过合成地震记录生成模块93在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录；通过合成地震记录修正模块94对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
通过本发明实施例提供的基于多子波的井震标定装置，采用不同频率不同相位的地震子波进行井震标定，能够使得生成的合成地震记录与井旁地震道数据更好地匹配，从而实现更准确的层位标定。
在一种可选的实施例中，本发明实施例提供的基于多子波的井震标定装置中，合成地震记录生成模块93可以包括：第一计算子模块931，用于根据测井数据中包含的声波曲线和密度曲线，计算反射系数序列；第二计算子模块932，用于计算反射系数序列与不同时窗内不同频率的子波的褶积，生成测井数据对应的合成地震记录。
可选地，合成地震记录生成模块93还可以包括：雷克子波提取子模块933，用于提取不同时窗对应的不同主频的雷克子波；其中，第二计算子模块932还用于计算反射系数序列与提取的不同主频的雷克子波的褶积，生成测井数据对应的合成地震记录。
进一步地，合成地震记录生成模块93还可以包括：地震子波提取子模块934，用于从井旁地震道数据中，提取不同时窗对应的不同频率的地震子波；其中，第二计算子模块932还用于计算反射系数序列与提取的不同频率的地震子波的褶积，生成测井数据对应的合成地震记录。
在一种可选的实施例中，本发明实施例提供的基于多子波的井震标定装置还可以包括：测井数据预处理模块95，用于对测井数据中包含的声波曲线和密度曲线进行环境校正和异常值处理。
本发明实施例中还提供了一种计算机设备，用以解决现有技术中，由于实测地震道数据纵向分辨率不同，使得采用单一频率的地震子波对测井数据进行井震标定得到的合成地震记录，与井旁地震道数据匹配结果较差，影响到层位标定准确性的技术问题，该计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述任一项的基于多子波的井震标定方法。
本发明实施例中还提供了一种计算机可读存储介质，用以解决现有技术中，由于实测地震道数据纵向分辨率不同，使得采用单一频率的地震子波对测井数据进行井震标定得到的合成地震记录，与井旁地震道数据匹配结果较差，影响到层位标定准确性的技术问题，该计算机可读存储介质存储有执行上述任一项的基于多子波的井震标定方法的计算机程序。
综上所述，本发明实施例中提供了一种基于多子波的井震标定方法、装置、计算机设备及计算机可读存储介质，利用二维或三维地震数据以及常规测井数据(例如声波、密度测井数据)，根据声波曲线和密度曲线计算得到反射系数序列，同时分时窗提取多个雷克子波和井旁地震道子波，采用反射系数序列和多个子波褶积得到合成地震记录，再通过微调不同时窗内子波的相位，使得生成的合成地震记录与井旁地震道数据达到最大相关性，提高了地质层位和时间域地震反射界面的匹配准确度。
实验证明，仅采用单个固定频率的子波进行合成地震记录标定，会产生上部合成记录波形匹配好而中下部波形匹配不好，或者上下部合成记录波形对应良好而中部对应不好的现象，这都是因为地震数据在整个井段的信噪比和分辨率不同造成的。可见，现有技术采用单一固定频率的子波对测井数据进行井震标定得到的合成地震记录，不能反应实际的地震资料情况；而通过本发明提供的基于不同频率不同相位的多子波进行井震标定的方法，能够使得生成的合成地震记录与井旁地震道数据更好地匹配，从而实现更准确的层位标定，为下一步地质构造解释与成图与勘探目标评价奠定坚实的基础数据。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (12)
1.一种基于多子波的井震标定方法，其特征在于，包括：
获取目标井位置处采集的测井数据和井旁地震道数据；
根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗；
在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录；
对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
2.如权利要求1所述的方法，其特征在于，在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录，包括：
根据测井数据中包含的声波曲线和密度曲线，计算反射系数序列；
计算反射系数序列与不同时窗内不同频率的子波的褶积，生成测井数据对应的合成地震记录。
3.如权利要求2所述的方法，其特征在于，计算反射系数序列与不同时窗内不同频率的子波的褶积，生成测井数据对应的合成地震记录，包括：
提取不同时窗对应的不同主频的雷克子波；
计算反射系数序列与提取的不同主频的雷克子波的褶积，生成测井数据对应的合成地震记录。
4.如权利要求2所述的方法，其特征在于，计算反射系数序列与不同时窗内不同频率的子波的褶积，生成测井数据对应的合成地震记录，包括：
从井旁地震道数据中，提取不同时窗对应的不同频率的地震子波；
计算反射系数序列与提取的不同频率的地震子波的褶积，生成测井数据对应的合成地震记录。
5.如权利要求2所述的方法，其特征在于，在根据测井数据中包含的声波曲线和密度曲线，计算反射系数序列之前，所述方法还包括：
对测井数据中包含的声波曲线和密度曲线进行环境校正和异常值处理。
6.一种基于多子波的井震标定装置，其特征在于，包括：
数据集采集模块，用于获取目标井位置处采集的测井数据和井旁地震道数据；
地震记录时窗划分模块，用于根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗；
合成地震记录生成模块，用于在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录；
合成地震记录修正模块，用于对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。
7.如权利要求6所述的装置，其特征在于，所述合成地震记录生成模块包括：
第一计算子模块，用于根据测井数据中包含的声波曲线和密度曲线，计算反射系数序列；
第二计算子模块，用于计算反射系数序列与不同时窗内不同频率的子波的褶积，生成测井数据对应的合成地震记录。
8.如权利要求7所述的装置，其特征在于，所述合成地震记录生成模块还包括：
雷克子波提取子模块，用于提取不同时窗对应的不同主频的雷克子波；
所述第二计算子模块还用于计算反射系数序列与提取的不同主频的雷克子波的褶积，生成测井数据对应的合成地震记录。
9.如权利要求7所述的装置，其特征在于，所述合成地震记录生成模块包括：
地震子波提取子模块，用于从井旁地震道数据中，提取不同时窗对应的不同频率的地震子波；
所述第二计算子模块还用于计算反射系数序列与提取的不同频率的地震子波的褶积，生成测井数据对应的合成地震记录。
10.如权利要求7所述的装置，其特征在于，所述装置还包括：
测井数据预处理模块，用于对测井数据中包含的声波曲线和密度曲线进行环境校正和异常值处理。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一项所述基于多子波的井震标定方法。
12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一项所述基于多子波的井震标定方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010115378.1A CN113376696A (zh) | 2020-02-25 | 2020-02-25 | 基于多子波的井震标定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010115378.1A CN113376696A (zh) | 2020-02-25 | 2020-02-25 | 基于多子波的井震标定方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113376696A true CN113376696A (zh) | 2021-09-10 |
Family
ID=77568461
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010115378.1A Pending CN113376696A (zh) | 2020-02-25 | 2020-02-25 | 基于多子波的井震标定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113376696A (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101013161A (zh) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
| CN105182408A (zh) * | 2015-08-28 | 2015-12-23 | 中国石油天然气集团公司 | 一种合成地震记录的制作方法和装置 |
| CN105759312A (zh) * | 2016-03-01 | 2016-07-13 | 中国海洋石油总公司 | 一种储层特征法井震标定方法 |
-
2020
- 2020-02-25 CN CN202010115378.1A patent/CN113376696A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101013161A (zh) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
| CN105182408A (zh) * | 2015-08-28 | 2015-12-23 | 中国石油天然气集团公司 | 一种合成地震记录的制作方法和装置 |
| CN105759312A (zh) * | 2016-03-01 | 2016-07-13 | 中国海洋石油总公司 | 一种储层特征法井震标定方法 |
Non-Patent Citations (2)
| Title |
|---|
| 冯晅 等: "分时窗提取地震子波及在合成地震记录中的应用", 《地球物理学进展》, vol. 17, no. 1, pages 71 - 75 * |
| 李金磊 等: "分时窗多子波方法在制作高精度合成记录中地应用", 《物探化探计算技术》, vol. 33, no. 1, pages 40 - 44 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104656142B (zh) | 一种利用垂直地震剖面与测井联合的地震层位标定方法 | |
| CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
| CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
| CN108802812A (zh) | 一种井震融合的地层岩性反演方法 | |
| CN105044777B (zh) | 基于经验模态分解检测地震标志层强反射振幅消除的方法 | |
| CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
| CN104237945B (zh) | 一种地震资料自适应高分辨处理方法 | |
| CN109425900B (zh) | 一种地震储层预测方法 | |
| CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
| CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
| CN104536042A (zh) | 一种二维叠后地震资料振幅补偿方法及装置 | |
| CN104459794A (zh) | 共反射点道集时变时间差值的校正方法及装置 | |
| CN104570078A (zh) | 一种基于频率域倾角的相似性横向变化率的溶洞检测方法 | |
| CN111208564A (zh) | 一种深度域层位标定方法及装置 | |
| CN110244383B (zh) | 基于近地表数据的地质岩性综合模型创建方法 | |
| Chen et al. | Improving the Precision of Surface Seismic Data Processing by Walkaway VSP | |
| CN113376696A (zh) | 基于多子波的井震标定方法及装置 | |
| CN110703329A (zh) | 基于弱振幅地震反射形成机制的岩性油藏边界确定方法 | |
| CN109283577B (zh) | 一种地震层位标定方法 | |
| CN112147717B (zh) | 一种叠前叠后井震联合时深标定的方法 | |
| CN116500679A (zh) | 一种井震标定方法、装置和相关设备 | |
| CN113589365B (zh) | 基于时频域信息的储层尖灭线描述方法 | |
| CN111060968B (zh) | 一种用于球面扩散补偿精确井控速度校正方法 | |
| CN110568491B (zh) | 一种品质因子q的估算方法 |
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 |