CN112904422B - 一种拾取***震源地震数据初至时间起跳点的方法及*** - Google Patents
一种拾取***震源地震数据初至时间起跳点的方法及*** Download PDFInfo
- Publication number
- CN112904422B CN112904422B CN201911133018.8A CN201911133018A CN112904422B CN 112904422 B CN112904422 B CN 112904422B CN 201911133018 A CN201911133018 A CN 201911133018A CN 112904422 B CN112904422 B CN 112904422B
- Authority
- CN
- China
- Prior art keywords
- time
- arrival
- peak
- seismic data
- moment
- 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
- 239000002360 explosive Substances 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims abstract description 80
- 230000008569 process Effects 0.000 claims description 19
- 238000005070 sampling Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 18
- 238000009499 grossing Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 10
- 230000003247 decreasing effect Effects 0.000 claims description 4
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 claims description 2
- 238000012937 correction Methods 0.000 abstract description 25
- 230000003068 static effect Effects 0.000 abstract description 24
- 238000013508 migration Methods 0.000 abstract description 11
- 230000005012 migration Effects 0.000 abstract description 11
- 238000012545 processing Methods 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 230000005284 excitation Effects 0.000 description 5
- 238000004587 chromatography analysis Methods 0.000 description 4
- 238000003325 tomography Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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
技术领域
本发明涉及地震资料处理的初至拾取技术领域，尤指一种拾取***震源地震数据初至时间起跳点的方法及***。
背景技术
随着油田勘探开发的不断深入，大规模的构造型油藏已经被发现，为了发现更为隐蔽的油气藏，更加精确的地震成像成为目前地震勘探技术发展的主要方向。
在提高地震资料成像的处理技术中，静校正和叠前深度偏移速度建模技术是非常关键的两项技术。静校正量可以通过多种静校正方法获得：
如，高程静校正是最简单的计算静校正量的方法，但由于其未考虑表层低降速带厚度速度变化，精确程度较低；
野外模型静校正是利用近地表表层调查数据得到近地表低降速带模型，可以较好的恢复近地表速度模型，但由于表层调查点密度较为稀疏，其精度受到较大限制；
利用单炮初至反演近地表的速度、厚度模型、计算静校正量的方法，由于炮点密度远远大于近地表调查点密度，其建立的浅表层速度模型更加准确，利用该速度模型计算的静校正量精度更高，叠加成像效果更好。而在叠前深度偏移速度建模中，浅表层速度模型的精度至关重要，常规通过速度拾取得到浅层速度模型的方式不仅耗时耗力，且由于地震资料浅层覆盖次数低、信噪比低，难以求取相对准确的浅表层速度，而利用单炮初至反演得到的近地表的速度模型精度更高，在速度建模过程中通过嵌入浅表层速度模型来确立更为准确的浅表层速度场，提高了速度模型的精确程度，对叠前深度偏移成像精度的提高有着重要的作用。
地震资料单炮初至拾取在地震资料处理中是必不可少的环节，初至拾取的精度和效率决定了近地表速度模型和静校正的计算精度，也影响着地震资料的处理周期。目前常用的初至拾取算法有相关法、能量比法、分形分维法、人工神经网络法和基于图像处理的方法等。
地震资料单炮初至拾取是进行初至波层析反演的基础。对于常规的初至拾取，由于初至波的波峰振幅强度大，稳定性好，易于辨认，因此在初至拾取时往往是对初至的波峰位置进行拾取。对于可控震源的数据来说，由于其地震资料是由震源信号与接收信号相关获得，因此信号是零相位子波，其数据中的初至波波峰位置可以认为是初至时间；而对于***震源的数据，由于其震源特性，获得的地震子波为最小相位的，其波峰起跳位置与真实的初至时间有一定的差异，因此在初至时间拾取时，拾取波峰位置是不能准确代表初至达到时间。初至时间的准确拾取是保证精确反演近地表速度建模的必备条件，且进一步影响静校正量的计算精度以及后续深度偏移速度模型的精度。作为整个地震资料数据处理的起始过程，准确拾取初至时间的意义重大。
通过对实际地震记录中的地震子波进行研究分析，发现了***震源的地震子波起跳点位置与理论上的初至时间是吻合的，因此，可以考虑以地震子波起跳点作为***震源数据的初至时间。但是由于受环境噪音等各种噪音的影响，起跳点位置处的波形特征不易识别，直接拾取会大大降低拾取精度和工作效率，因此在拾取初至时大多拾取波峰位置。目前，还缺乏考虑***震源数据的起跳点拾取方法，这对于获得高精度的初至波层析反演数据是不利的，同时也进一步的影响了后续高精度地震剖面的获得。
因此，亟需一种可以准确拾取***震源地震数据初至时间起跳点的技术方案。
发明内容
为解决上述问题，本发明提出了一种拾取***震源地震数据初至时间起跳点的方法及***，结合实际情况中***震源激发产生最小相位子波波峰位置偏离初至时间的影响，利用起跳点位置代替波峰位置作为初至时间，可以准确拾取***震源地震数据初至时间起跳点，对后续的静校正处理以及叠前深度偏移的近地表速度模型建立提供更为有效的数据支持。
在本发明一实施例中，提出了一种拾取***震源地震数据初至时间起跳点的方法，该方法包括：
获取地震资料；
根据所述地震资料，拾取***震源地震数据初至波峰位置；
根据所述***震源地震数据初至波峰位置，截取一段地震记录，根据截取的所述地震记录重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal；
根据截取的所述地震记录，对相应的子波进行拟合，得到与所述子波近似的正弦波，并确定所述正弦波的波峰位置；
在拟合得到的所述正弦波中确定第一个零值的时刻，作为波形起跳点的位置；
根据所述波形起跳点的位置对应的时刻以及所述正弦波的波峰位置对应的时刻，得到波峰时刻与起跳时刻的时差Δt；
根据所述实际地震数据中的波峰初至时刻treal与时差Δt，得到初至起跳点时刻。
在本发明另一实施例中，还提出了一种拾取***震源地震数据初至时间起跳点的***，该***包括：
地震资料获取模块，用于获取地震资料；
初至波峰位置拾取模块，用于根据所述地震资料，拾取***震源地震数据初至波峰位置；
波峰初至时刻计算模块，用于根据所述***震源地震数据初至波峰位置，截取一段地震记录，根据截取的所述地震记录重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal；
正弦波拟合模块，用于根据截取的所述地震记录，对相应的子波进行拟合，得到与所述子波近似的正弦波，并确定所述正弦波的波峰位置；
起跳点位置获取模块，用于在拟合得到的所述正弦波中确定第一个零值的时刻，作为波形起跳点的位置；
时差计算模块，用于根据所述波形起跳点的位置对应的时刻以及所述正弦波的波峰位置对应的时刻，得到波峰时刻与起跳时刻的时差Δt；
初值起跳点时刻计算模块，用于根据所述实际地震数据中的波峰初至时刻treal与时差Δt，得到初至起跳点时刻。
在本发明另一实施例中，还提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现拾取***震源地震数据初至时间起跳点的方法。
在本发明另一实施例中，还提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现拾取***震源地震数据初至时间起跳点的方法。
本发明提出的拾取***震源地震数据初至时间起跳点的方法及***，可以结合实际情况中***震源激发产生最小相位子波波峰位置偏离初至时间的影响，利用起跳点位置代替波峰位置作为初至时间，准确拾取***震源地震数据初至时间起跳点，使初至时间更为准确，对后续的静校正处理以及叠前深度偏移的近地表速度模型建立提供更为有效的数据支持。
附图说明
图1是本发明一实施例的拾取***震源地震数据初至时间起跳点的方法流程示意图。
图2是本发明一实施例的获取实际地震数据中的波峰初至时刻的方法流程示意图。
图3是本发明一实施例的拾取***震源地震数据初至时间起跳点的***架构示意图。
图4是本发明一实施例的计算机设备结构示意图。
图5是本发明一具体实施例的波峰初至时间拾取示意图。
图6是本发明一具体实施例的单道中波峰初至时间附近确定的子波形态示意图。
图7是本发明一具体实施例的拟合得到与子波形态一致的正弦波形示意图。
图8是本发明一具体实施例的标记出起跳点位置的正弦波示意图。
图9是本发明一具体实施例的单道记录中的起跳点初至位置示意图。
图10A是本发明一具体实施例的波峰初至时间层析静校正的共炮检距剖面示意图。
图10B是本发明一具体实施例的起跳点初至时间层析静校正的共炮检距剖面示意图。
图11A是本发明一具体实施例的嵌入波峰初至层析近地表速度模型的偏移剖面示意图。
图11B是本发明一具体实施例的起跳点初至层析近地表模型的偏移剖面示意图。
具体实施方式
下面将参考若干示例性实施方式来描述本发明的原理和精神。应当理解，给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明，而并非以任何方式限制本发明的范围。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本领域技术人员知道，本发明的实施方式可以实现为一种***、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
首先，需要说明的术语有：
初至拾取：初至波为地震记录中最先到达接收点出的地震波信息，初至拾取是获得地震波初至的一种手段。对于***震源来说，其子波类型为最小相位子波，它的特性是能量集中在波形的前部，但是有一个从弱逐渐到强的过程。初至波决定着地震波在地下介质中的传播时间，进一步反应了地震波经过的介质的速度。对于最小相位子波，其能量最大的位置并不能代表地震波最先到达检波点的时间。拾取***震源的地震记录初至时间时，应该选择能够代表地震波最先到达检波点时间的起跳点位置作为初至点。
初至波层析：初至波层析是使用初至波作为输入的观测数据，通过模型正演，得到理论模型下的计算数据。通过修改理论模型逐步减小计算数据与观测数据的差。当计算数据与观测数据的差达到足够小的误差范围内时，则可认为理论模型与真实地下模型是一致的。
根据本发明的实施方式，提出了一种拾取***震源地震数据初至时间起跳点的方法及***，该方法及***在原有拾取波峰初至的基础上，将初至重新定义到更为合理的起跳点位置，进而提高初至时间的准确程度，为初至波层析提供更准确的初至时间数据，从而更准确的建立地表速度模型，进一步的提高静校正的质量以及后续偏移的浅表层速度模型质量，为地下精确成像提供有利的技术支撑。
下面参考本发明的若干代表性实施方式，详细阐释本发明的原理和精神。
图1是本发明一实施例的拾取***震源地震数据初至时间起跳点的方法流程示意图。如图1所示，该方法包括：
步骤S101，获取地震资料。
其中，地震资料可以是***震源原始单炮记录。
步骤S102，根据所述地震资料(即，***震源原始单炮记录)，拾取***震源地震数据初至波峰位置。在本步骤的基础上，继续进行后续的起跳点时刻计算过程。
步骤S103，根据所述***震源地震数据初至波峰位置，截取一段地震记录，该地震记录作为计算起跳点时刻的基础，根据截取的所述地震记录重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal。
本步骤的具体过程参考图2所示，包括：
步骤S1031，以所述***震源地震数据初至波峰位置的时刻为起始，向时间减少方向进行搜索并判断得到子波的负振幅位置，进一步得到向时间减少方向搜索过程中，距离所述***震源地震数据初至波峰位置对应的波峰时刻最近的振幅小于0的时刻t1；其中，
时刻t1的表达式为：A(t1)＜0；其中，A为振幅函数。
步骤S1032，以时刻t1对应的负振幅位置为起始，向时间减少方向进行搜索并判断得到子波的起始位置，进一步得到向时间减少方向搜索过程中，距离时刻t1最近的振幅大于0的时刻t2；其中，
时刻t2的表达式为：A(t2)＞0；其中，A为振幅函数。
步骤S1033，以所述***震源地震数据初至波峰位置的时刻为起始，向时间增加方向进行搜索并判断得到子波的终止位置，进一步得到向时间增加方向搜索过程中，距离所述***震源地震数据初至波峰位置对应的波峰时刻最近的振幅大于0的时刻t3；其中，
时刻t3的表达式为：A(t3)＜0；其中，A为振幅函数。
步骤S1034，根据时刻t2及时刻t3，分别向时间减少方向和时间增加方向延伸一定数量的采样点，截取当前地震道上的子波，并以此作为计算起跳点的基础。这里延伸的采样点数量可以由操作人员视地震资料品质而定。
步骤S1035，根据截取的所述当前地震道上的子波，重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal，计算公式如下：
treal＝t2+n×s；
其中，treal为实际地震数据中的波峰初至时刻；n为根据时刻t2向时间减少方向延伸的采样点数，n可以是负值；s为所述地震资料的采样率。
在完成步骤S103之后，进一步需要根据截取的地震记录，计算波形起跳点的位置对应的时刻t0与正弦波中的波峰对应的时刻t'real的时间差，具体过程如步骤S104至步骤S106所述。
步骤S104，根据截取的所述地震记录，对相应的子波进行拟合，得到与所述子波近似的正弦波，并确定所述正弦波的波峰位置。其中，
根据截取的所述地震记录，基于采样点振幅值对子波进行拟合，得到一个近似的正弦波：
其中，Areal为所述地震记录对应的采样点振幅值；以所述地震记录对应的子波的起始位置为0时刻，t为所在的时刻；为正弦波的相位；ω为正弦波的频率；A为拟合得到的正弦波振幅值。这里的未知量为拟合正弦波振幅值A、频率ω及相位/>已知数据为时刻t及采样点振幅值Areal。
步骤S105，在拟合得到的所述正弦波中确定第一个零值的时刻，作为波形起跳点的位置。
具体过程为：根据拟合得到的所述正弦波，选取第一个振幅值为零值的情况：
得到正弦波中第一个零值的时刻t0，作为波形起跳点的位置。
步骤S106，根据所述波形起跳点的位置对应的时刻t0以及所述正弦波的波峰位置对应的时刻t'real，得到波峰时刻与起跳时刻的时差Δt。
波峰时刻与起跳时刻的时差Δt是利用以下公式计算得到：
Δt＝t'real-t0；
其中，t'real是所述正弦波的波峰位置对应的时刻；t0为所述波形起跳点的位置对应的时刻。
步骤S107，根据所述实际地震数据中的波峰初至时刻treal与时差Δt，得到初至起跳点时刻。
在一具体实施例中，为准确获取初始起跳点时刻，在执行完步骤S106之后，可以对多个波峰时刻与起跳时刻的时差Δt进行平滑，得到平滑处理后的以此来保证得到稳定的初至时间；
平滑处理的公式如下：
其中，为平滑处理后的时差；Δti为第i个波峰时刻与起跳时刻的时差；trn为设定的时窗范围内参与计算的道数；mid表示将时窗内的时差按照从小到大进行排列，并选择位于/>向上取整的位置处的时差；可以根据设定的时窗范围选定式(1)或式(2)进行平滑处理。
具体的，时窗范围可以由操作人员给定，并选用合适的平滑公式进行平滑处理；根据所述地震资料中的初至时间，如果初至资料较好，例如，初至较为明显的资料，则可以选择局部选择中值滤波处理，即式(2)；如果初至资料较差，例如，资料中的初至时间跳动不稳定，则可以选择线性平滑处理，即式(1)。
结合步骤S107，可以利用以下公式计算得到初至起跳点时刻：
其中，tcal为初至起跳点时刻；是对多个波峰时刻与起跳时刻的时差Δt平滑处理后的时差。
利用上述方法拾取***震源地震数据初至时间起跳点时，考虑到了实际情况中***震源激发产生最小相位子波波峰位置偏离初至时间的影响，从而改用起跳点位置代替波峰位置作为初至时间，解决了原有的***震源初至时间拾取与真实初至时间有差异的问题，使初至时间更为准确；该方法易于计算机自动化实现，对于以***震源激发的地震资料有很好的效果，为后续的静校正以及进一步的叠前深度偏移等工作提供了良好的保证。
需要说明的是，尽管在上述实施例及附图中以特定顺序描述了本发明方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
在介绍了本发明示例性实施方式的方法之后，接下来，参考图3对本发明示例性实施方式的拾取***震源地震数据初至时间起跳点的***进行介绍。
拾取***震源地震数据初至时间起跳点的***的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的术语“模块”或者“单元”，可以是实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
基于同一发明构思，本发明还提出了一种拾取***震源地震数据初至时间起跳点的***，如图3所示，该***包括：
地震资料获取模块310，用于获取地震资料；
初至波峰位置拾取模块320，用于根据所述地震资料，拾取***震源地震数据初至波峰位置；
波峰初至时刻计算模块330，用于根据所述***震源地震数据初至波峰位置，截取一段地震记录，根据截取的所述地震记录重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal；
正弦波拟合模块340，用于根据截取的所述地震记录，对相应的子波进行拟合，得到与所述子波近似的正弦波，并确定所述正弦波的波峰位置；
起跳点位置获取模块350，用于在拟合得到的所述正弦波中确定第一个零值的时刻，作为波形起跳点的位置；
时差计算模块360，用于根据所述波形起跳点的位置对应的时刻以及所述正弦波的波峰位置对应的时刻，得到波峰时刻与起跳时刻的时差Δt；
初值起跳点时刻计算模块370，用于根据所述实际地震数据中的波峰初至时刻treal与时差Δt，得到初至起跳点时刻。
应当注意，尽管在上文详细描述中提及了拾取***震源地震数据初至时间起跳点的***的若干模块，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多模块的特征和功能可以在一个模块中具体化。反之，上文描述的一个模块的特征和功能可以进一步划分为由多个模块来具体化。
基于前述发明构思，如图4所示，本发明还提出了一种计算机设备400，包括存储器410、处理器420及存储在存储器410上并可在处理器420上运行的计算机程序430，所述处理器420执行所述计算机程序430时实现前述拾取***震源地震数据初至时间起跳点的方法。
基于前述发明构思，本发明还提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现拾取***震源地震数据初至时间起跳点的方法。
为了对上述拾取***震源地震数据初至时间起跳点的方法、***、计算机设备及计算机可读存储介质进行更为清楚的解释，下面结合一个具体的实施例来进行说明，然而值得注意的是该实施例仅是为了更好地说明本发明，并不构成对本发明不当的限定。
以塔里木盆地的某区域为例，采集该区域实际***震源激发的地震记录(***震源原始单炮记录)，运用本发明图1所示的方法对拾取的波峰初至进行调整，使初至时间变为起跳点初至。
具体流程如下：
步骤S1，首先，根据***震源原始单炮记录，如进行波峰初至时间的拾取，图5所示是波峰初至时间拾取示意图。
步骤S2，在单道记录中，根据波峰初至附近采样点的采样点值情况，截取一段包含能够代表该段处波形特征的记录，如图6所示，是本具体实施例的单道中波峰初至时间附近确定的子波形态示意图。
步骤S3，如图7所示，通过对子波701的采样点进行拟合，得到一个主频与地震子波相近的正弦波702，该正弦波形与子波形态一致。
步骤S4，基于图7拟合得到的正弦波形，得到第一个0点的位置即为起跳点的位置，该起跳点的位置如图8中的标记“×”所示。；
步骤S5，在单道中根据起跳点位置与波峰位置的时间差计算平移所需的时间，并根据多道的时间差对整体进行平滑处理，防止出现异常起跳点的出现。
步骤S6，根据计算得到的每一道时差，对单道中的波峰初至位置进行平移，得到单道记录中的起跳点初至位置，如图9所示，标记线901标示出了起跳点初至位置。
更进一步的，步骤S7，使用计算得到的起跳点初至进行初至波层析反演，得到静校正量，对地震数据进行静校正；如图10A所示，为波峰初至时间层析静校正的共炮检距剖面示意图；如图10B所示，为起跳点初至时间层析静校正的共炮检距剖面示意图。
步骤S8，将近地表速度模型嵌入整体的偏移速度模型并进行叠前深度偏移；如图11A所示，为嵌入波峰初至层析近地表速度模型的偏移剖面示意图；如图11B所示，为起跳点初至层析近地表模型的偏移剖面示意图。
通过上述步骤S7及步骤S8分别得到的共炮检距剖面及偏移剖面示意图来看，由于本发明改用起跳点位置代替波峰位置作为初至时间，解决了原有的***震源初至时间拾取与真实初至时间有差异的问题，使初至时间更准确，建立的地表速度模型也更准确，为后续的静校正以及进一步的叠前深度偏移等工作提供了良好的保证，提高了静校正的质量以及后续偏移的浅表层速度模型质量。
利用上述方法，在塔里木盆地探区进行推广应用，共完成总工区面积约5000平方公里的***震源起跳点初至时间拾取，为浅表层速度模型的建立提供了更为准确的初至信息，对后续的静校正处理以及叠前深度偏移的近地表速度模型建立提供了更为有效的数据，为地下精确成像提供有利的技术支撑。
本发明提出的拾取***震源地震数据初至时间起跳点的方法及***，可以结合实际情况中***震源激发产生最小相位子波波峰位置偏离初至时间的影响，利用起跳点位置代替波峰位置作为初至时间，准确拾取***震源地震数据初至时间起跳点，使初至时间更为准确，对后续的静校正处理以及叠前深度偏移的近地表速度模型建立提供更为有效的数据支持。
虽然已经参考若干具体实施方式描述了本发明的精神和原理，但是应该理解，本发明并不限于所公开的具体实施方式，对各方面的划分也不意味着这些方面中的特征不能组合以进行受益，这种划分仅是为了表述的方便。本发明旨在涵盖所附权利要求的精神和范围内所包括的各种修改和等同布置。
Claims (8)
1.一种拾取***震源地震数据初至时间起跳点的方法，其特征在于，该方法包括：
获取地震资料；
根据所述地震资料，拾取***震源地震数据初至波峰位置；
根据所述***震源地震数据初至波峰位置，截取一段地震记录，根据截取的所述地震记录重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal；
根据截取的所述地震记录，对相应的子波进行拟合，得到与所述子波近似的正弦波，并确定所述正弦波的波峰位置；
在拟合得到的所述正弦波中确定第一个零值的时刻，作为波形起跳点的位置；
根据所述波形起跳点的位置对应的时刻以及所述正弦波的波峰位置对应的时刻，得到波峰时刻与起跳时刻的时差Δt；
根据所述实际地震数据中的波峰初至时刻treal与时差Δt，得到初至起跳点时刻；
其中，根据所述***震源地震数据初至波峰位置，截取一段地震记录，根据截取的所述地震记录重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal，包括：
以所述***震源地震数据初至波峰位置的时刻为起始，向时间减少方向进行搜索并判断得到子波的负振幅位置，进一步得到向时间减少方向搜索过程中，距离所述***震源地震数据初至波峰位置对应的波峰时刻最近的振幅小于0的时刻t1；
以时刻t1对应的负振幅位置为起始，向时间减少方向进行搜索并判断得到子波的起始位置，进一步得到向时间减少方向搜索过程中，距离时刻t1最近的振幅大于0的时刻t2；
以所述***震源地震数据初至波峰位置的时刻为起始，向时间增加方向进行搜索并判断得到子波的终止位置，进一步得到向时间增加方向搜索过程中，距离所述***震源地震数据初至波峰位置对应的波峰时刻最近的振幅小于0的时刻t3；
根据时刻t2及时刻t3，分别向时间减少方向和时间增加方向延伸一定数量的采样点，截取当前地震道上的子波；
根据截取的所述当前地震道上的子波，重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal，计算公式如下：
treal＝t2+n×s；
其中，treal为实际地震数据中的波峰初至时刻；n为根据时刻t2向时间减少方向延伸的采样点数；s为所述地震资料的采样率；
其中，根据截取的所述地震记录，对相应的子波进行拟合，得到与所述子波近似的正弦波，并确定所述正弦波的波峰位置，包括：
根据采样点振幅值，对子波进行拟合，得到一个近似的正弦波：
其中，Areal为所述地震记录对应的采样点振幅值；以所述地震记录对应的子波的起始位置为0时刻，t为所在的时刻；为正弦波的相位；ω为正弦波的频率；A为拟合得到的正弦波振幅值；
其中，在拟合得到的所述正弦波中确定第一个零值的时刻，作为波形起跳点的位置，包括：
根据拟合得到的所述正弦波，选取第一个振幅值为零值的情况：
得到正弦波中第一个零值的时刻t0，作为波形起跳点的位置；
其中，根据所述波形起跳点的位置对应的时刻以及所述正弦波的波峰位置对应的时刻，得到波峰时刻与起跳时刻的时差Δt，包括：
利用以下公式计算得到波峰时刻与起跳时刻的时差Δt：
Δt＝t'real-t0；
其中，t'real是所述正弦波的波峰位置对应的时刻，即，波峰时刻；t0为所述波形起跳点的位置对应的时刻，即，起跳时刻。
2.根据权利要求1所述的拾取***震源地震数据初至时间起跳点的方法，其特征在于，获取地震资料，包括***震源原始单炮记录。
3.根据权利要求2所述的拾取***震源地震数据初至时间起跳点的方法，其特征在于，根据所述地震资料，拾取***震源地震数据初至波峰位置，包括：
在所述***震源原始单炮记录中，拾取***震源地震数据初至波峰位置。
4.根据权利要求1所述的拾取***震源地震数据初至时间起跳点的方法，其特征在于，该方法还包括：
对多个波峰时刻与起跳时刻的时差进行平滑，得到平滑处理后的Δt；平滑处理的公式如下：
其中，为平滑处理后的时差；Δti为第i个波峰时刻与起跳时刻的时差；trn为设定的时窗范围内参与计算的道数；mid表示将时窗内的时差按照从小到大进行排列，并选择位于向上取整的位置处的时差；
根据所述地震资料中的初至时间，如果初至时间跳动不稳定选择式(1)进行平滑处理，否则选择式(2)。
5.根据权利要求4所述的拾取***震源地震数据初至时间起跳点的方法，其特征在于，根据所述实际地震数据中的波峰初至时刻treal与时差Δt，得到初至起跳点时刻，还包括：
利用以下公式计算得到初至起跳点时刻：
其中，tcal为初至起跳点时刻；是对多个波峰时刻与起跳时刻的时差Δt平滑处理后的时差。
6.一种拾取***震源地震数据初至时间起跳点的***，其特征在于，该***包括：
地震资料获取模块，用于获取地震资料；
初至波峰位置拾取模块，用于根据所述地震资料，拾取***震源地震数据初至波峰位置；
波峰初至时刻计算模块，用于根据所述***震源地震数据初至波峰位置，截取一段地震记录，根据截取的所述地震记录重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal；
正弦波拟合模块，用于根据截取的所述地震记录，对相应的子波进行拟合，得到与所述子波近似的正弦波，并确定所述正弦波的波峰位置；
起跳点位置获取模块，用于在拟合得到的所述正弦波中确定第一个零值的时刻，作为波形起跳点的位置；
时差计算模块，用于根据所述波形起跳点的位置对应的时刻以及所述正弦波的波峰位置对应的时刻，得到波峰时刻与起跳时刻的时差Δt；
初值起跳点时刻计算模块，用于根据所述实际地震数据中的波峰初至时刻treal与时差Δt，得到初至起跳点时刻；
其中，根据所述***震源地震数据初至波峰位置，截取一段地震记录，根据截取的所述地震记录重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal，包括：
以所述***震源地震数据初至波峰位置的时刻为起始，向时间减少方向进行搜索并判断得到子波的负振幅位置，进一步得到向时间减少方向搜索过程中，距离所述***震源地震数据初至波峰位置对应的波峰时刻最近的振幅小于0的时刻t1；
以时刻t1对应的负振幅位置为起始，向时间减少方向进行搜索并判断得到子波的起始位置，进一步得到向时间减少方向搜索过程中，距离时刻t1最近的振幅大于0的时刻t2；
以所述***震源地震数据初至波峰位置的时刻为起始，向时间增加方向进行搜索并判断得到子波的终止位置，进一步得到向时间增加方向搜索过程中，距离所述***震源地震数据初至波峰位置对应的波峰时刻最近的振幅小于0的时刻t3；
根据时刻t2及时刻t3，分别向时间减少方向和时间增加方向延伸一定数量的采样点，截取当前地震道上的子波；
根据截取的所述当前地震道上的子波，重新定义波峰初至时刻，得到实际地震数据中的波峰初至时刻treal，计算公式如下：
treal＝t2+n×s；
其中，treal为实际地震数据中的波峰初至时刻；n为根据时刻t2向时间减少方向延伸的采样点数；s为所述地震资料的采样率；
其中，根据截取的所述地震记录，对相应的子波进行拟合，得到与所述子波近似的正弦波，并确定所述正弦波的波峰位置，包括：
根据采样点振幅值，对子波进行拟合，得到一个近似的正弦波：
其中，Areal为所述地震记录对应的采样点振幅值；以所述地震记录对应的子波的起始位置为0时刻，t为所在的时刻；为正弦波的相位；ω为正弦波的频率；A为拟合得到的正弦波振幅值；
其中，在拟合得到的所述正弦波中确定第一个零值的时刻，作为波形起跳点的位置，包括：
根据拟合得到的所述正弦波，选取第一个振幅值为零值的情况：
得到正弦波中第一个零值的时刻t0，作为波形起跳点的位置；
其中，根据所述波形起跳点的位置对应的时刻以及所述正弦波的波峰位置对应的时刻，得到波峰时刻与起跳时刻的时差Δt，包括：
利用以下公式计算得到波峰时刻与起跳时刻的时差Δt：
Δt＝t'real-t0；
其中，t'real是所述正弦波的波峰位置对应的时刻，即，波峰时刻；t0为所述波形起跳点的位置对应的时刻，即，起跳时刻。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1至5任一所述方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911133018.8A CN112904422B (zh) | 2019-11-19 | 2019-11-19 | 一种拾取***震源地震数据初至时间起跳点的方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911133018.8A CN112904422B (zh) | 2019-11-19 | 2019-11-19 | 一种拾取***震源地震数据初至时间起跳点的方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112904422A CN112904422A (zh) | 2021-06-04 |
| CN112904422B true CN112904422B (zh) | 2024-03-26 |
Family
ID=76103157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911133018.8A Active CN112904422B (zh) | 2019-11-19 | 2019-11-19 | 一种拾取***震源地震数据初至时间起跳点的方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112904422B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102243320A (zh) * | 2011-04-12 | 2011-11-16 | 中国石油天然气股份有限公司 | 一种地震波初至拾取方法 |
| CN103616722A (zh) * | 2013-11-28 | 2014-03-05 | 中国石油天然气股份有限公司 | 一种初至波旅行时拾取方法及装置 |
| CN108072896A (zh) * | 2016-11-18 | 2018-05-25 | 中国石油化工股份有限公司 | 一种全自动地震波初至拾取方法及*** |
| CN108279436A (zh) * | 2017-12-20 | 2018-07-13 | 中国石油天然气股份有限公司 | 一种拾取地震波初至的方法及装置 |
| CN109655885A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 一种自动初至拾取方法及*** |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1325030A (zh) * | 2000-05-18 | 2001-12-05 | 中国石油化工股份有限公司***分公司 | ***机精确井口时间获取方法与装置 |
| US9103935B2 (en) * | 2013-02-04 | 2015-08-11 | King Fahd University Of Petroleum And Minerals | Method of first arrival picking of seismic refraction data |
| US10345463B2 (en) * | 2014-02-21 | 2019-07-09 | Ion Geophysical Corporation | Methods and systems for using known source events in seismic data processing |
| US11372124B2 (en) * | 2018-04-18 | 2022-06-28 | GeoTomo, LLC | First-break picking of seismic data and generating a velocity model |
-
2019
- 2019-11-19 CN CN201911133018.8A patent/CN112904422B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102243320A (zh) * | 2011-04-12 | 2011-11-16 | 中国石油天然气股份有限公司 | 一种地震波初至拾取方法 |
| CN103616722A (zh) * | 2013-11-28 | 2014-03-05 | 中国石油天然气股份有限公司 | 一种初至波旅行时拾取方法及装置 |
| CN108072896A (zh) * | 2016-11-18 | 2018-05-25 | 中国石油化工股份有限公司 | 一种全自动地震波初至拾取方法及*** |
| CN109655885A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 一种自动初至拾取方法及*** |
| CN108279436A (zh) * | 2017-12-20 | 2018-07-13 | 中国石油天然气股份有限公司 | 一种拾取地震波初至的方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| Automatic first-arrival picking with Reinforcement Learning;Yue Ma等;《CPS/SEG Beijing 2018 International Geophysical Conference & Exposition Electronic papers》;第515-519页 * |
| 初至智能拾取技术;刘志成;《石油物探》;第46卷(第05期);第521-530页 * |
| 基于时频分析的初至拾取方法研究;岳龙;刘怀山;刘凯;徐秀刚;邢磊;《石油物探》;第54卷(第05期);第508-520页 * |
| 高精度初至自动拾取综合方法研究;徐钰;段卫星;徐维秀;曾维辉;徐雷良;《物探与化探》;第34卷(第05期);第595-599页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112904422A (zh) | 2021-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106353792B (zh) | 一种适用于水力压裂微震震源定位的方法 | |
| RU2694621C1 (ru) | Способ и устройство для обработки сейсмических данных | |
| CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
| CN104570076A (zh) | 一种基于二分法的地震波初至自动拾取方法 | |
| CN103616722A (zh) | 一种初至波旅行时拾取方法及装置 | |
| CN102937721A (zh) | 利用初至波走时的有限频层析成像方法 | |
| CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
| CN106772593B (zh) | 绕射波的成像方法及装置 | |
| CN105093318B (zh) | 一种自适应波动方程波场延拓静校正方法 | |
| CN108375789B (zh) | 联合采集地震数据的同步匹配方法 | |
| GB2520128A (en) | Method of Calibrating A geologic Model | |
| WO2021155754A1 (zh) | 光纤声波传感地震数据中井筒波干扰的去除方法及装置 | |
| CN112904422B (zh) | 一种拾取***震源地震数据初至时间起跳点的方法及*** | |
| CN105445793A (zh) | 一种确定坏道数据的方法及装置 | |
| CN108693560A (zh) | 一种基于互相关道的散射波成像方法及*** | |
| CN111538082B (zh) | 一种地震波时频域初至自动拾取方法 | |
| CN104278990A (zh) | 测井数据质量恢复方法及装置 | |
| CN110109184B (zh) | 一种基于多日变点的被动场源类三维电场勘探方法 | |
| CN109143364B (zh) | 初至时间拾取方法 | |
| CN113552623A (zh) | 曲流河沉积砂体边界的确定方法和确定装置 | |
| CN107356970B (zh) | 一种高精度的井震数据匹配方法 | |
| CN113433588B (zh) | 一种分偏移距扫描叠加的近地表速度分析方法 | |
| CN114859414B (zh) | 从地震数据中自动提取地层倾角信息的方法和装置 | |
| CN111123356A (zh) | 基于初至信息的异常道智能判识方法 | |
| CN111596357B (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 |