CN112379439A - 地震数据中纵波与横波匹配的方法及装置 - Google Patents
地震数据中纵波与横波匹配的方法及装置 Download PDFInfo
- Publication number
- CN112379439A CN112379439A CN202011202647.4A CN202011202647A CN112379439A CN 112379439 A CN112379439 A CN 112379439A CN 202011202647 A CN202011202647 A CN 202011202647A CN 112379439 A CN112379439 A CN 112379439A
- Authority
- CN
- China
- Prior art keywords
- waves
- longitudinal
- transverse
- wave
- pure
- 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 49
- 238000011160 research Methods 0.000 claims abstract description 107
- 238000012545 processing Methods 0.000 claims abstract description 35
- 238000005553 drilling Methods 0.000 claims abstract description 20
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 8
- 238000001914 filtration Methods 0.000 claims description 26
- 238000012795 verification Methods 0.000 claims description 24
- 238000004458 analytical method Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 16
- 230000015572 biosynthetic process Effects 0.000 claims description 12
- 238000010183 spectrum analysis Methods 0.000 claims description 11
- 238000001228 spectrum Methods 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 9
- 239000011435 rock Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 16
- 238000010586 diagram Methods 0.000 description 18
- 230000008859 change Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004519 manufacturing process 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
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种地震数据中纵波与横波匹配的方法及装置，其中该方法包括：根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。本发明可实现将纵波匹配到纯横波的时间域，可避免现有技术中气云区位置会产生同相轴下拉的畸变以及丢失原始数据信息的问题，改善了地震数据中纵波与横波匹配效果。
Description
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及地震数据中纵波与横波匹配的方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
在多波资料应用过程中，纵横波数据的匹配是纵横波联合解释、联合属性分析、叠后与叠前纵横波联合反演及岩性识别、储层预测和含气性识别的重要基础，是体现多波多分量地震勘探技术实际勘探开发应用价值的关键。
但纵观国内外纵波横波匹配技术的发展，几乎所有的方法均基于纵波和转换横波同一地层反射同相轴相似性最大的原则求取纵横波的速度比，进而实现纵波和转换横波在时间域的匹配，将转换横波匹配到纵波时间域。
对于地震数据中纵波与纯横波的匹配，由于纵波数据存在气云区构造下拉的特征，而纯横波无下拉特征，现阶段在实现中纵波与纯横波的匹配时，一般是将纯横波配到纵波的时间域。
但在将纯横波配到纵波的时间域后，会造成气云区位置产生同相轴下拉的畸变的问题，这不符合应用纯横波恢复气云区构造的初衷，而且匹配后的横波会丢失原始数据的部分信息，影响了地震数据中纵波与横波的匹配效果。
发明内容
本发明实施例提供一种地震数据中纵波与横波匹配的方法，用以改善地震数据中纵波与横波匹配效果，该方法包括：
根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；
对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；
根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；
根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。
本发明实施例还提供一种地震数据中纵波与横波匹配的装置，用以改善地震数据中纵波与横波匹配效果，该装置包括：
合成记录标定模块，用于根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；
地震层位解释处理模块，用于对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；
空间纵横波速度比预测模块，用于根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；
匹配模块，用于根据研究区的空间纵横波速度比，根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地震数据中纵波与横波匹配的方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述地震数据中纵波与横波匹配的方法的计算机程序。
本发明实施例中，根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，与现有技术中将纯横波匹配到纵波的时间域的技术方案对比，可实现将纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，可避免现有技术中气云区位置会产生同相轴下拉的畸变以及丢失原始数据信息的问题，改善了地震数据中纵波与横波匹配效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中一种地震数据中纵波与横波匹配的方法的流程示意图；
图2为本发明实施例中一种地震数据中纵波与横波匹配的方法的具体示例图；
图3为本发明实施例中原始纵波、横波地震数据剖面的示意图；
图4a为本发明实施例中横波数据匹配到纵波时间域的成果剖面的示意图；
图4b为本发明实施例中纵波数据匹配到横波时间域的成果剖面的示意图；
图5为本发明实施例中纵波数据与匹配后纵波数据(SS域)、横波数据与匹配后横波数据(PP域)瞬时振幅变化的对比图；
图6为本发明实施例中纵波数据、匹配后纵波数据(SS域)、横波数据的频谱分析示意图；
图7为本发明实施例中原始横波数据、匹配后纵波数据(SS域)与分频段带通滤波结果的对比图；
图8a，图8b，图8c，图8d为本发明实施例中应用地震数据中纵波与横波匹配的方法进行层位调整前后的匹配效果的对比图；
图9a为本发明实施例中将横波数据匹配到纵波时间域后联合反演结果的剖面图；
图9b为本发明实施例中将纵波数据匹配到横波时间域后联合反演结果的剖面图；
图10为本发明实施例中一种地震数据中纵波与横波匹配的装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1为本发明实施例中一种地震数据中纵波与横波匹配的方法的流程示意图,如图1所示，本发明实施例提供的地震数据中纵波与横波匹配的方法，可以包括：
步骤101：根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；
步骤102：对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；
步骤103：根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；
步骤104：根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。
本发明实施例中，根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，与现有技术中将纯横波匹配到纵波的时间域的技术方案对比，可实现将纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，可避免现有技术中气云区位置会产生同相轴下拉的畸变以及丢失原始数据信息的问题，改善了地震数据中纵波与横波匹配效果。
具体实施时，首先根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定。
实施例中，根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定，可以包括：对研究区测井资料和钻井资料的纵波数据和横波数据进行分析，确定地层纵波和横波的速度特征；根据地层纵波和横波的速度特征，对研究区地震数据中纵波和纯横波进行合成记录标定。
在上述实施例中，首先可以开展对研究区内有关纵波横波数据的测井资料的分析，明确地层纵波横波的速度特征，完成对纵波和纯横波地震数据的合成记录标定。
具体实施时，根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定，可以包括：在研究区的测井资料中不含有横波的测井数据时，以预设置的弹性参数，通过岩石物理建模，进行横波曲线估算；根据估算的横波曲线，对研究区地震数据中纯横波进行合成记录标定。
在上述实施例中，在研究区的测井资料中不含有横波的测井数据时(即无横波测井情况)，基于岩石物理建模的横波数据标定技术，通过预设置合理的弹性参数以及预选择的横波预测方法，开展岩石物理建模进行横波曲线的估算；应用估算的横波曲线进行纯横波地震数据的合成记录标定。
具体实施时，在根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定后，对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据。
实施例中，对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据，可以包括：根据进行合成记录标定后的纵波和纯横波，确定目的层的反射特征；根据目的层的反射特征，对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据。
在上述实施例中，通过确定目的层的反射特征，确定地震资料解释方案(即选择地震层位解释处理的具体处理方式)，分别完成纵波数据、纯横波数据的初步层位解释，得到纵波和纯横波的解释层位数据。
具体实施时，本发明实施例提供的地震数据中纵波与横波匹配的方法，还可以包括：根据研究区地震数据中纵波和纯横波，预计算井点位置目的层纵横波速度比。
实施例中，可按如下公式预计算井点位置目的层纵横波速度比：
其中，γ为井点位置目的层纵横波速度比；TSS为研究区地震数据中横波数据反射时长；TPP为研究区地震数据中纵波数据反射时长；h为研究区的反射地层深度；Vpp为研究区地震数据中纵波速度；VSS为研究区地震数据中横波速度。
在上述实施例中，通过预计算井点位置目的层纵横波速度比，有助于在后续步骤中预测研究区的空间纵横波速度比。
具体实施时，在对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据后，根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比。
在上述实施例中，可以叠后地震数据中纯横波数据为基础，首先确定井点位置目的层的纵横波速度比，通过纵横波解释层位数据的匹配，预测空间上的纵横波速度比。
具体实施时，根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。
实施例中，匹配过程是将纵波数据匹配到纯横波数据时间域，匹配方法是以叠后地震数据中纯横波数据为基础，将空间纵横波速度比的数据体，应用于纵波数据，实现将纵波地震数据匹配到横波数据时间域的目的，由此得到与横波匹配的横波时间域纵波。
在上述实施例中，可将纵波数据匹配到纯横波的时间域，从而基本解决了现有技术下将横波数据匹配到纵波时间域存在的问题，在此基础上实现了纵波横波联合解释和联合反演，进而可实现根据纵波横波进行岩性识别、储层预测和含气性识别的目的。
现有技术中下，纯横波数据与纵波数据在同一地层波组特征具有明显的差异，对于大套地层，波组强弱对应关系较好，但是纯横波数据的同相轴数据多于纵波数据，因而同一套地层在纵波、纯横波的剖面上并不具备明显相似的波形、振幅、频率等特征，难以进行对比分析，难以应用同一地层反射同相轴相似性最大的原则进行评价。因此纯横波数据与纵波数据的匹配，目前没有适用的匹配评价方法。本发明实施例，与现有技术中将纯横波匹配到纵波的时间域的技术方案对比，可实现将纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，可实现对横波纵波的匹配效果的验证与评价。
具体实施时，本发明实施例提供的地震数据中纵波与横波匹配的方法，还可以包括：根据地震数据中纯横波和预设置的滤波参数，对与横波匹配的横波时间域纵波进行验证；若验证通过，则输出与横波匹配的横波时间域纵波。
实施例中，根据地震数据中纯横波和预设置的滤波参数，对与横波匹配的横波时间域纵波进行验证，可以包括：分别对地震数据中纯横波和与横波匹配的横波时间域纵波进行频谱分析；分别对经频谱分析后的地震数据中纯横波和与横波匹配的横波时间域纵波，以预设置的滤波参数进行分频段带通滤波分析，得到纯横波和纵波的分频段分析结果；所述纯横波和纵波的分频段分析结果含有纯横波和纵波对应不同频段的频谱特征，波组强弱特征和同相轴；对纯横波和纵波的分频段分析结果中不同频段进行验证，在频段满足在该频段下纯横波和纵波的频谱特征的相似程度大于预设值、波组强弱特征的相似程度大于预设值且同相轴一一对应时，验证通过。
实施例中，还可以包括：若验证不通过，则重新对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位更新数据；根据纵波和纯横波的解释层位更新数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比的更新数据；根据研究区的空间纵横波速度比的更新数据，将地震数据中纵波与纯横波进行匹配，得到与横波匹配的横波时间域纵波的更新数据；重复执行上述步骤，直到在根据地震数据中纯横波，对与横波匹配的横波时间域纵波进行验证时验证通过。
在上述实施例中，可根据地震数据中纯横波和预设置的滤波参数，对与横波匹配的横波时间域纵波进行验证，实现对横波纵波的匹配效果的评价。将原始横波数据、匹配后纵波数据分别进行频谱分析，然后进行分频段带通滤波，选择两者频谱特征最接近的频段，要求在两者在该频段下经带通滤波后的两组数据波组的强弱特征相似，且同相轴一一对应；如果找不到符合上述条件的滤波频段，则需要重新调整解释层位，重新对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，直到根据地震数据中的纯横波，对与横波匹配的横波时间域纵波进行验证时验证通过。在验证通过时，可确定实现地震数据中纯横波，和与横波匹配的横波时间域纵波，二者的优势频段内同相轴一一对应。上述步骤中，重新进行地震反射层位解释并进行纵横波的匹配是一个迭代的过程。本发明实施例最终可实现基于优势频带同相轴相似性，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。
下面给出一个具体实施例，来说明本发明的方法的具体应用。在本实施例中，将与纵波匹配的纵波时间域横波用纯横波数据(PP域)表示；将与横波匹配的横波时间域纵波用纵波数据(SS域)表示。如图2所示，该具体实施例可以包括如下步骤：
1、开展对研究区内纵波横波测井资料和钻井资料的分析，明确地层纵波横波的速度特征；当研究区没有可以应用的横波测井数据，通过对试油和解释成果数据的分析，应用岩石物理建模技术进行横波速度估算，通过正演模拟分析波组特征，应用估算的横波速度结合波组特征进行横波数据的合成记录标定；结合区域地质资料对纵波数据进行合成记录标定，综合纵、横波数据标定进行连井地层对比，确定纵波数据和横波数据的解释层位数据，建立井-震地层统一格架并完成初步纵横波数据的联合解释；
2、根据纵横波资料解释的层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比，据此进行纵波数据和横波数据的匹配。图3展示了原始纵波和横波地震数据；纵波数据存在气云区构造下拉的特征，而纯横波数据无下拉，恢复了构造形态。借鉴转换横波与纵波的匹配方法，将纯横波数据匹配到纵波的时间域的结果，将与纵波匹配的纵波时间域横波用纯横波数据(PP域)表示，如图4a所示，可以看出气云区的位置同样会产生同相轴下拉的畸变，这不符合应用纯横波数据恢复气云区构造的初衷；而将纵波数据匹配到纯横波时间域，将与横波匹配的横波时间域纵波用纵波数据(SS域)表示，如图4b所示，可看出与横波匹配的横波时间域纵波数据的构造形态与纯横波数据一致；
3、对比两种匹配方法，选择不受气云区影响的斜坡位置进行对比，图5展示了纵波数据与纵波数据(SS域)、纯横波数据与纯横波数据(PP域)的瞬时振幅变化对比，可以看出纵波数据(SS域)的振幅与原始纵波数据相比变化很小，基本可以忽略；而纯横波数据(PP域)的振幅与原始纯横波数据相比，数值范围及分布形态均变化较大，也就是说纯横波数据(PP域)丢失了原始纯横波数据的部分信息；因此可确定本发明实施例提供的方法对纵波与横波匹配效果更好；
4、匹配结果是否合理的判断方法，由于纵波与转换横波的匹配是基于纵波和转换横波同一地层反射同相轴相似性最大的原则，这一方法不适用于纯横波数据，因此本发明实施例创新了基于优势频带同相轴相似的方法。通过将原始纯横波数据、与横波匹配的横波时间域纵波分别进行频谱分析，然后进行带通滤波分析，其中，选取的带通滤波参数要一致。带通滤波频带范围可分别预设置为3-13Hz、12-22Hz、22-32Hz三个频段，如图6所示；
图7为原始地震数据、3-13Hz、12-22Hz、22-32Hz三个频段的滤波结果的示意图，图7中每一部分的图示中黑色竖线左侧均为横波数据、右侧均为纵波数据(SS域)。从图7可以看出，与原始数据相比，滤波后三个频段的结果横波数据与纵波数据(SS域)其振幅能量强弱变化规律更为接近，尤其是12-22Hz这一频段两者的波组的对应关系基本一致，且同相轴一一对应，因此选择12-22Hz滤波参数进行匹配精度分析；
5、在匹配的过程中有两个关键点，一是滤波参数的选择。需要注意，同一测线不同位置频谱特征不同，同一滤波参数不能保证这一剖面目的层内所有的同相轴都达到一一对应的关系。二是纵波与横波解释层位一一对应关系是正确的，如图8a和图8b，由于横波数据K5和K9层位整体追踪靠上，使得滤波后同相轴不能一一对应；图8c和图8d为调整层位后滤波效果，对应关系较好。因此地震反射层位解释与纵横波匹配是迭代的过程，直至滤波后对应关系良好，就可以应用匹配的纵横波数据联合进行储层预测、含气检测和构造成图；
6、应用效果：对比纵波横波联合反演结果，如图9a所示，图9a展示了横波数据匹配到纵波域后的纵横波的联合反演结果，将纵波数据匹配到纯横波的时间域后，联合反演结果的分辨率与纯横波数据分辨能力一致；如图9b所示，图9b展示了纵波数据匹配到横波域后的纵横波的联合反演结果，将纯横波数据匹配到纵波时间域后，联合反演结果的分辨率低于横波数据的分辨率，高于原始纵波分辨率。
本发明实施例中，根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，与现有技术中将纯横波匹配到纵波的时间域的技术方案对比，可实现将纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，可避免现有技术中气云区位置会产生同相轴下拉的畸变以及丢失原始数据信息的问题，改善了地震数据中纵波与横波匹配效果。
本发明实施例中还提供了一种地震数据中纵波与横波匹配的装置，如下面的实施例所述。由于该装置解决问题的原理与地震数据中纵波与横波匹配的方法相似，因此该装置的实施可以参见地震数据中纵波与横波匹配的方法的实施，重复之处不再赘述。
图10为本发明实施例中一种地震数据中纵波与横波匹配的装置的结构示意图，如图10所示，本发明实施例提供的地震数据中纵波与横波匹配的装置，可以包括：
合成记录标定模块01，用于根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；
地震层位解释处理模块02，用于对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；
空间纵横波速度比预测模块03，用于根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；
匹配模块04，用于根据研究区的空间纵横波速度比，根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。
在一个实施例中，合成记录标定模块，具体用于：对研究区测井资料和钻井资料的纵波数据和横波数据进行分析，确定地层纵波和横波的速度特征；
根据地层纵波和横波的速度特征，对研究区地震数据中纵波和纯横波进行合成记录标定。
在一个实施例中，合成记录标定模块，具体用于：在研究区的测井资料中不含有横波的测井数据时，以预设置的弹性参数，通过岩石物理建模，进行横波曲线估算；
根据估算的横波曲线，对研究区地震数据中纯横波进行合成记录标定。
在一个实施例中，地震层位解释处理模块，具体用于：
根据进行合成记录标定后的纵波和纯横波，确定目的层的反射特征；
根据目的层的反射特征，对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据。
在一个实施例中，还包括：井点位置目的层纵横波速度比计算模块，用于：根据研究区地震数据中纵波和纯横波，预计算井点位置目的层纵横波速度比。
在一个实施例中，井点位置目的层纵横波速度比计算模块，具体用于：按如下公式预计算井点位置目的层纵横波速度比：
其中，γ为井点位置目的层纵横波速度比；TSS为研究区地震数据中横波数据反射时长；TPP为研究区地震数据中纵波数据反射时长；h为研究区的反射地层深度；Vpp为研究区地震数据中纵波速度；VSS为研究区地震数据中横波速度。
在一个实施例中，还包括：验证模块，用于：根据地震数据中纯横波和预设置的滤波参数，对与横波匹配的横波时间域纵波进行验证；若验证通过，则输出与横波匹配的横波时间域纵波。
在一个实施例中，验证模块，具体用于：
分别对地震数据中纯横波和与横波匹配的横波时间域纵波进行频谱分析；
分别对经频谱分析后的地震数据中纯横波和与横波匹配的横波时间域纵波，以预设置的滤波参数进行分频段带通滤波分析，得到纯横波和纵波的分频段分析结果；所述纯横波和纵波的分频段分析结果含有纯横波和纵波对应不同频段的频谱特征，波组强弱特征和同相轴；
对纯横波和纵波的分频段分析结果中不同频段进行验证，在频段满足在该频段下纯横波和纵波的频谱特征的相似程度大于预设值、波组强弱特征的相似程度大于预设值且同相轴一一对应时，验证通过。
在一个实施例中，验证模块，具体用于：
若验证不通过，则重新对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位更新数据；根据纵波和纯横波的解释层位更新数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比的更新数据；根据研究区的空间纵横波速度比的更新数据，将地震数据中纵波与纯横波进行匹配，得到与横波匹配的横波时间域纵波的的更新数据；重复执行上述步骤，直到在根据地震数据中纯横波，对与横波匹配的横波时间域纵波进行验证时验证通过。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地震数据中纵波与横波匹配的方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述地震数据中纵波与横波匹配的方法的计算机程序。
本发明实施例中，根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，与现有技术中将纯横波匹配到纵波的时间域的技术方案对比，可实现将纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波，可避免现有技术中气云区位置会产生同相轴下拉的畸变以及丢失原始数据信息的问题，改善了地震数据中纵波与横波匹配效果。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (20)
1.一种地震数据中纵波与横波匹配的方法，其特征在于，包括：
根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；
对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；
根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；
根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。
2.如权利要求1所述的方法，其特征在于，根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定，包括：
对研究区测井资料和钻井资料的纵波数据和横波数据进行分析，确定地层纵波和横波的速度特征；
根据地层纵波和横波的速度特征，对研究区地震数据中纵波和纯横波进行合成记录标定。
3.如权利要求1所述的方法，其特征在于，根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定，包括：
在研究区的测井资料中不含有横波的测井数据时，以预设置的弹性参数，通过岩石物理建模，进行横波曲线估算；
根据估算的横波曲线，对研究区地震数据中纯横波进行合成记录标定。
4.如权利要求1所述的方法，其特征在于，对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据，包括：
根据进行合成记录标定后的纵波和纯横波，确定目的层的反射特征；
根据目的层的反射特征，对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据。
5.如权利要求1所述的方法，其特征在于，还包括：根据研究区地震数据中纵波和纯横波，预计算井点位置目的层纵横波速度比。
6.如权利要求5所述的方法，其特征在于，按如下公式预计算井点位置目的层纵横波速度比：
其中，γ为井点位置目的层纵横波速度比；TSS为研究区地震数据中横波数据反射时长；TPP为研究区地震数据中纵波数据反射时长；h为研究区的反射地层深度；Vpp为研究区地震数据中纵波速度；VSS为研究区地震数据中横波速度。
7.如权利要求1所述的方法，其特征在于，还包括：根据地震数据中纯横波和预设置的滤波参数，对与横波匹配的横波时间域纵波进行验证；若验证通过，则输出与横波匹配的横波时间域纵波。
8.如权利要求7所述的方法，其特征在于，根据地震数据中纯横波和预设置的滤波参数，对与横波匹配的横波时间域纵波进行验证，包括：
分别对地震数据中纯横波和与横波匹配的横波时间域纵波进行频谱分析；
分别对经频谱分析后的地震数据中纯横波和与横波匹配的横波时间域纵波，以预设置的滤波参数进行分频段带通滤波分析，得到纯横波和纵波的分频段分析结果；所述纯横波和纵波的分频段分析结果含有纯横波和纵波对应不同频段的频谱特征，波组强弱特征和同相轴；
对纯横波和纵波的分频段分析结果中不同频段进行验证，在频段满足在该频段下纯横波和纵波的频谱特征的相似程度大于预设值、波组强弱特征的相似程度大于预设值且同相轴一一对应时，验证通过。
9.如权利要求8所述的方法，其特征在于，若验证不通过，则重新对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位更新数据；根据纵波和纯横波的解释层位更新数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比的更新数据；根据研究区的空间纵横波速度比的更新数据，将地震数据中纵波与纯横波进行匹配，得到与横波匹配的横波时间域纵波的更新数据；重复执行上述步骤，直到在根据地震数据中纯横波，对与横波匹配的横波时间域纵波进行验证时验证通过。
10.一种地震数据中纵波与横波匹配的装置，其特征在于，包括：
合成记录标定模块，用于根据研究区的测井资料和钻井资料，对研究区地震数据中纵波和纯横波进行合成记录标定；
地震层位解释处理模块，用于对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据；
空间纵横波速度比预测模块，用于根据纵波和纯横波的解释层位数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比；
匹配模块，用于根据研究区的空间纵横波速度比，根据研究区的空间纵横波速度比，将地震数据中纵波匹配到纯横波的时间域，得到与横波匹配的横波时间域纵波。
11.如权利要求10所述的装置，其特征在于，合成记录标定模块，具体用于：对研究区测井资料和钻井资料的纵波数据和横波数据进行分析，确定地层纵波和横波的速度特征；
根据地层纵波和横波的速度特征，对研究区地震数据中纵波和纯横波进行合成记录标定。
12.如权利要求10所述的装置，其特征在于，合成记录标定模块，具体用于：在研究区的测井资料中不含有横波的测井数据时，以预设置的弹性参数，通过岩石物理建模，进行横波曲线估算；
根据估算的横波曲线，对研究区地震数据中纯横波进行合成记录标定。
13.如权利要求10所述的装置，其特征在于，地震层位解释处理模块，具体用于：
根据进行合成记录标定后的纵波和纯横波，确定目的层的反射特征；
根据目的层的反射特征，对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位数据。
14.如权利要求10所述的装置，其特征在于，还包括：井点位置目的层纵横波速度比计算模块，用于：根据研究区地震数据中纵波和纯横波，预计算井点位置目的层纵横波速度比。
15.如权利要求14所述的装置，其特征在于，井点位置目的层纵横波速度比计算模块，具体用于：按如下公式预计算井点位置目的层纵横波速度比：
其中，γ为井点位置目的层纵横波速度比；TSS为研究区地震数据中横波数据反射时长；TPP为研究区地震数据中纵波数据反射时长；h为研究区的反射地层深度；Vpp为研究区地震数据中纵波速度；VSS为研究区地震数据中横波速度。
16.如权利要求10所述的装置，其特征在于，还包括：验证模块，用于：根据地震数据中纯横波和预设置的滤波参数，对与横波匹配的横波时间域纵波进行验证；若验证通过，则输出与横波匹配的横波时间域纵波。
17.如权利要求16所述的装置，其特征在于，验证模块，具体用于：
分别对地震数据中纯横波和与横波匹配的横波时间域纵波进行频谱分析；
分别对经频谱分析后的地震数据中纯横波和与横波匹配的横波时间域纵波，以预设置的滤波参数进行分频段带通滤波分析，得到纯横波和纵波的分频段分析结果；所述纯横波和纵波的分频段分析结果含有纯横波和纵波对应不同频段的频谱特征，波组强弱特征和同相轴；
对纯横波和纵波的分频段分析结果中不同频段进行验证，在频段满足在该频段下纯横波和纵波的频谱特征的相似程度大于预设值、波组强弱特征的相似程度大于预设值且同相轴一一对应时，验证通过。
18.如权利要求17所述的装置，其特征在于，验证模块，具体用于：
若验证不通过，则重新对进行合成记录标定后的纵波和纯横波进行地震层位解释处理，得到纵波和纯横波的解释层位更新数据；根据纵波和纯横波的解释层位更新数据，和预计算的井点位置目的层纵横波速度比，预测研究区的空间纵横波速度比的更新数据；根据研究区的空间纵横波速度比的更新数据，将地震数据中纵波与纯横波进行匹配，得到与横波匹配的横波时间域纵波的的更新数据；重复执行上述步骤，直到在根据地震数据中纯横波，对与横波匹配的横波时间域纵波进行验证时验证通过。
19.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至9任一所述方法。
20.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至9任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011202647.4A CN112379439A (zh) | 2020-11-02 | 2020-11-02 | 地震数据中纵波与横波匹配的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011202647.4A CN112379439A (zh) | 2020-11-02 | 2020-11-02 | 地震数据中纵波与横波匹配的方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112379439A true CN112379439A (zh) | 2021-02-19 |
Family
ID=74576509
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011202647.4A Pending CN112379439A (zh) | 2020-11-02 | 2020-11-02 | 地震数据中纵波与横波匹配的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112379439A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112859172A (zh) * | 2021-03-05 | 2021-05-28 | 中国石油天然气集团有限公司 | 纵横波微测井资料处理方法及装置 |
| CN113868979A (zh) * | 2021-09-29 | 2021-12-31 | 中国石油化工股份有限公司 | 浅层气钻井风险的确定方法、装置和服务器 |
| CN113917534A (zh) * | 2021-09-01 | 2022-01-11 | 展少辉 | 一种地球物理勘探用纵横波匹配方法和装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4686657A (en) * | 1984-09-24 | 1987-08-11 | Mobil Oil Corporation | Matching P wave and shear wave sections in seismic exploration |
| CN104570079A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种纵波、转换横波地震资料的时间匹配方法 |
| US20150234061A1 (en) * | 2006-04-17 | 2015-08-20 | Soundblast Technologies, Llc | System and Method for Harnessing Pressure Produced by a Detonation |
| CN106646611A (zh) * | 2016-12-08 | 2017-05-10 | 中国石油天然气集团公司 | 一种纵波和转换波地震数据时间匹配的方法及装置 |
| CN107589452A (zh) * | 2017-07-21 | 2018-01-16 | 中国石油天然气集团公司 | 纵波和转换波的数据匹配方法和装置 |
| CN110780351A (zh) * | 2018-07-31 | 2020-02-11 | 中国石油化工股份有限公司 | 纵波和转换波叠前联合反演方法及系统 |
-
2020
- 2020-11-02 CN CN202011202647.4A patent/CN112379439A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4686657A (en) * | 1984-09-24 | 1987-08-11 | Mobil Oil Corporation | Matching P wave and shear wave sections in seismic exploration |
| US20150234061A1 (en) * | 2006-04-17 | 2015-08-20 | Soundblast Technologies, Llc | System and Method for Harnessing Pressure Produced by a Detonation |
| CN104570079A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种纵波、转换横波地震资料的时间匹配方法 |
| CN106646611A (zh) * | 2016-12-08 | 2017-05-10 | 中国石油天然气集团公司 | 一种纵波和转换波地震数据时间匹配的方法及装置 |
| CN107589452A (zh) * | 2017-07-21 | 2018-01-16 | 中国石油天然气集团公司 | 纵波和转换波的数据匹配方法和装置 |
| CN110780351A (zh) * | 2018-07-31 | 2020-02-11 | 中国石油化工股份有限公司 | 纵波和转换波叠前联合反演方法及系统 |
Non-Patent Citations (2)
| Title |
|---|
| 王九拴 等: "多波地震资料在三湖地区油气预测中的应用", 石油地球物理勘探, vol. 47, no. 4, 31 August 2012 (2012-08-31), pages 606 * |
| 陈国文等: "纵横波联合解释技术在气云区的应用", 岩性油气藏, vol. 31, no. 6, 31 December 2019 (2019-12-31), pages 81 - 82 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112859172A (zh) * | 2021-03-05 | 2021-05-28 | 中国石油天然气集团有限公司 | 纵横波微测井资料处理方法及装置 |
| CN113917534A (zh) * | 2021-09-01 | 2022-01-11 | 展少辉 | 一种地球物理勘探用纵横波匹配方法和装置 |
| CN113868979A (zh) * | 2021-09-29 | 2021-12-31 | 中国石油化工股份有限公司 | 浅层气钻井风险的确定方法、装置和服务器 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112379439A (zh) | 地震数据中纵波与横波匹配的方法及装置 | |
| US8923093B2 (en) | Determining the quality of a seismic inversion | |
| CN102109612B (zh) | 一种地震波吸收衰减补偿方法 | |
| CN110703354B (zh) | 一种基于井震结合的岩性边界识别方法 | |
| CN107678064B (zh) | 一种声波时差实时提取方法 | |
| CN105425289A (zh) | 确定低频波阻抗的方法和装置 | |
| US11237288B2 (en) | Verifying measurements of elastic anisotropy parameters in an anisotropic wellbore environment | |
| CN111856566A (zh) | 湖相滩坝砂体中薄储层预测方法及装置 | |
| CN112946743B (zh) | 区分储层类型的方法 | |
| CN110118994B (zh) | 一种基于地震反演和机器学习的陆相烃源岩定量预测方法 | |
| CN111239809B (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
| CN112649887A (zh) | 基于钻井资料定量划分层序地层的方法及装置 | |
| CN112327357A (zh) | 基于三维地震资料的烃源岩成熟度预测方法及装置 | |
| CN107679614B (zh) | 一种基于粒子群优化的声波时差实时提取方法 | |
| CN112394392B (zh) | 对烃源岩分布情况进行评价的方法和装置 | |
| Torii et al. | Shear slowness evaluation by parametric inversion of dipole waveforms | |
| CN114114421A (zh) | 基于深度学习的导向自学习地震数据去噪方法及装置 | |
| CN112147687B (zh) | 一种储层含气性预测方法及预测系统 | |
| CN112130209A (zh) | 岩溶储层预测方法及装置 | |
| CN113589366B (zh) | 基于全波形反演的宽频融合建模方法 | |
| Wang et al. | Accurately estimating shear slowness using data-driven quadrupole sonic logging-while-drilling data processing | |
| CN112147700A (zh) | 速度异常区的低频模型构建方法及系统 | |
| CN113009579B (zh) | 地震数据反演方法及装置 | |
| CN113514904B (zh) | 地层参数模型建立方法及装置 | |
| CN112363244B (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 |