CN112230285B - 地震数据激发延迟时间的校正方法及校正装置 - Google Patents
地震数据激发延迟时间的校正方法及校正装置 Download PDFInfo
- Publication number
- CN112230285B CN112230285B CN201910634891.9A CN201910634891A CN112230285B CN 112230285 B CN112230285 B CN 112230285B CN 201910634891 A CN201910634891 A CN 201910634891A CN 112230285 B CN112230285 B CN 112230285B
- Authority
- CN
- China
- Prior art keywords
- data
- seismic data
- post
- processing
- delay time
- 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
- 230000005284 excitation Effects 0.000 title claims abstract description 89
- 238000012937 correction Methods 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000012545 processing Methods 0.000 claims abstract description 101
- 238000004590 computer program Methods 0.000 claims description 10
- 239000013535 sea water Substances 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 abstract description 8
- 238000004891 communication Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 17
- 230000009471 action Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000010304 firing Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000013480 data collection Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000006467 substitution reaction 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/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/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/51—Migration
- G01V2210/514—Post-stack
-
- 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
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
技术领域
本发明涉及地震数据处理技术领域,具体涉及一种地震数据激发延迟时间的校正方法及校正装置。
背景技术
海上石油勘探利用拖缆进行地震数据采集时,为了施工的便利,在检波器开始记录数据后,作为震源的气枪一般都会延迟一段时间再进行激发。这就会造成检波器的原始记录时间与地震波实际旅行时间之间存在时间差值,该时间差值即为激发延迟时间。
在拖缆地震数据处理过程中,为了得到地震波实际旅行时间,在对地震数据进行处理时都会把激发延迟时间从原始记录时间中减去。采集拖缆地震数据时,检波器开始记录数据后间隔一激发延迟时间,气枪再进行激发。记录下该激发延迟时间并提供给后续的拖缆地震数据处理。
但是,由于气枪和检波器的时间同步误差、检波器陈放深度误差等会引起实际激发延迟时间与施工设计的激发延迟时间不同,存在实际的激发延迟时间与设计的激发延迟时间之间误差大、拖缆地震数据精度低的问题。
发明内容
针对现有技术中的问题,本发明提供一种地震数据激发延迟时间的校正方法及校正装置,能够有效提高激发延迟时间的精确性,进而对拖缆地震数据进行校正,提高拖缆地震数据精度以及提高基于拖缆地震数据的成像质量。
为解决上述技术问题,本发明提供以下技术方案:
第一方面,本发明提供一种地震数据激发延迟时间的校正方法,包括:
对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;
根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;
将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;
对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;
基于所述道数据中最大值所对应的时间对地震数据的激发延迟时间进行校正。
进一步的,所述对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据,包括:
抽取所述地震数据中的共中心点道集,采用第一速度对该共中心点道集进行动校正;
对所述共中心点道集进行叠加处理得到叠后地震数据。
进一步的,所述第一速度为地震波在海水中的传播速度。
进一步的,所述基于所述道数据中最大值所对应的时间对地震数据的激发延迟时间进行校正,包括:
地震数据的记录时间中减去所述道数据中最大值所对应的时间,完成对地震数据激发延迟时间的校正。
第二方面,本发明提供一种地震数据激发延迟时间的校正装置,包括:
动校正单元,用于对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;
褶积单元,用于根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;
处理单元,用于将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;
叠加单元,用于对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;
校正单元,用于基于所述道数据中最大值所对应的时间对地震数据的激发延迟时间进行校正。
进一步的,所述动校正单元包括:
动校正子单元,用于抽取所述地震数据中的共中心点道集,采用第一速度对该共中心点道集进行动校正;
叠加处理子单元,用于对所述共中心点道集进行叠加处理得到叠后地震数据。
进一步的,所述第一速度为地震波在海水中的传播速度。
进一步的,所述校正单元包括:
校正子单元,用于地震数据的记录时间中减去所述道数据中最大值所对应的时间,完成对地震数据激发延迟时间的校正。
第三方面,本发明提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现所述的地震数据激发延迟时间的校正方法的步骤。
第四方面,本发明提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现所述的地震数据激发延迟时间的校正方法的步骤。
由上述技术方案可知,本发明提供一种地震数据激发延迟时间的校正方法及校正装置,通过对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;基于所述道数据中最大值所对应的时间对地震数据激发延迟时间进行校正,能够有效提高激发延迟时间的精确性,为后续地震数据处理中提供更加精确的地震波实际旅行时间,进而对拖缆地震数据进行校正,提高拖缆地震数据精度以及提高基于拖缆地震数据的成像质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的地震数据激发延迟时间的校正装置的一种通信结构示意图。
图2为本发明的地震数据激发延迟时间的校正装置的另一种通信结构示意图。
图3为本发明实施例中的地震数据激发延迟时间的校正方法的流程示意图。
图4为本发明实施例中的地震数据激发延迟时间的校正装置的结构示意图。
图5为本发明实施例中的电子设备的结构示意图。
图6为本发明具体应用实例中的激发延迟时间校正前的炮记录图。
图7为本发明具体应用实例中的动校正和叠加后的地震数据图。
图8为本发明具体应用实例中地震数据自褶积的叠后多次波数据图。
图9为本发明具体应用实例中叠后地震数据和多次波数据进行互相关的互相关数据图。
图10为本发明具体应用实例中的互相关数据进行叠加的数据图。
图11为本发明具体应用实例中的激发延迟时间校正后的炮记录图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
考虑到现有的对地震波实际旅行时间进行校正的方法中存在的校正准确性低而无法保证地震波实际旅行时间的精度的问题。本发明提供一种拖缆地震数据激发延迟时间的校正方法、地震数据激发延迟时间的校正装置、电子设备及计算机可读存储介质,通过对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;基于所述道数据中最大值所对应的时间对地震数据激发延迟时间进行校正,能够有效提高激发延迟时间的精确性,为后续地震数据处理中提供更加精确的地震波实际旅行时间,进而对拖缆地震数据进行校正,提高拖缆地震数据精度以及提高基于拖缆地震数据的成像质量。
基于上述内容,本发明还提供一种地震数据激发延迟时间的校正装置,该装置可以为一种服务器A1,参见图1,该服务器A1可以与客户端设备B1通信连接,用户可以将地震资料及其他相关数据输入所述客户端设备B1,所述客户端设备B1可以在线将地震资料及其他相关数据发送至服务器A1,所述服务器A1可以在线接收所述客户端设备B1发送的地震资料及其他相关数据,而后离线或在线根据地震资料获取对应的地震数据,对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;基于所述道数据中最大值所对应的时间对地震数据激发延迟时间进行校正。而后,所述服务器A1在线将目标地震数据发送至所述客户端设备B1,使得用户经由所述客户端设备B1获知最终校正得到的目标地震数据。
进一步来说,上述服务器A1还可以与地震资料采集设备C1通信连接,参见图2,该地震资料采集设备C1可以直接自目标区域获取地震资料及其他相关数据,也可以与数据库D1通信连接,自该数据库D1中获取对应的地震资料及其他相关数据。而后,所述地震资料采集设备C1将地震资料及其他相关数据发送至所述服务器A1。
可以理解的是,所述客户端设备B1可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中,所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。
在实际应用中,进行地震数据激发延迟时间的校正的部分可以在如上述内容所述的服务器A1侧执行,即,如图1所示的架构,也可以所有的操作都在所述客户端设备B1中完成。具体可以根据所述客户端设备B1的处理能力,以及用户使用场景的限制等进行选择。本发明对此不作限定。若所有的操作都在所述客户端设备B1中完成,所述客户端设备B1还可以包括处理器,用于进行地震数据激发延迟时间的校正的具体处理。
上述的客户端设备可以具有通信模块(即通信单元),可以与远程的服务器进行通信连接,实现与所述服务器的数据传输。例如,通信单元可以将上述用户输入的地震资料及其他相关数据发送至服务器,以便服务器根据这些地震资料及其他相关数据进行地震数据激发延迟时间的校正。通信单元还可以接收服务器返回的校正结果。所述服务器可以包括任务调度中心一侧的服务器,其他的实施场景中也可以包括中间平台的服务器,例如与任务调度中心服务器有通信链接的第三方服务器平台的服务器。所述的服务器可以包括单台计算机设备,也可以包括多个服务器组成的服务器集群,或者分布式装置的服务器结构。
所述服务器与所述客户端设备之间可以使用任何合适的网络协议进行通信,包括在本发明提交日尚未开发出的网络协议。所述网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然,所述网络协议例如还可以包括在上述协议之上使用的RPC协议(Remote Procedure Call Protocol,远程过程调用协议)、REST协议(Representational State Transfer,表述性状态转移协议)等。
为了能够有效提高激发延迟时间的精确性,为后续地震数据处理中提供更加精确的地震波实际旅行时间并对拖缆地震数据进行校正。本发明提供一种地震数据激发延迟时间的校正方法的实施例,参见图3,该地震数据激发延迟时间的校正方法具体包含有如下内容:
S101:对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;
可以理解的是,该地震数据是通过接收的方式或者采集的方式得到的海上拖缆地震勘探数据,该地震勘探数据为叠前地震数据;
该叠前地震数据表示为d0(x0,ti),其中,x0是叠前数据对应的位置,ti是数据对应的采样时间,i是地震数据采样的序列,满足i=1,2,…,n。n为单道数据总的采样数目。
在本步骤中,抽取地震数据中的共中心点道集,采用第一速度对该共中心点道集进行动校正处理;其中,该第一速度为地震波在海水中的传播速度,在本实施例中,地震波在海水中的传播速度为1500m/s。
对该动校正处理之后的共中心点道集进行叠加处理得到叠后地震数据,该叠后地震数据表示为d(x,ti),其中,x是叠后数据对应的位置,ti是数据对应的采样时间,i是地震数据采样的序列,满足i=1,2,…,n。
需要说明的是,叠加处理得到的对应的叠后地震数据即为地震数据进行动校正处理和叠加处理后得到的叠后地震数据。
S102:根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;
在本步骤中,根据步骤S101确定的叠后地震数据d(x,ti),并对该叠后地震数据d(x,ti)进行自褶积处理得到叠后多次波数据f(x,tj)。
其中,按照下式对该叠后地震数据d(x,ti)进行自褶积处理得到叠后多次波数据f(x,tj):
其中,j=1,2,…,n。
需要说明的是,自褶积处理得到的对应的叠后多次波数据即为对叠后地震数据进行自褶积处理得到的叠后多次波数据。
S103:将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;
在本步骤中,叠后地震数据d(x,ti)和叠后多次波数据f(x,tj)进行互相关处理得到互相关数据r(x,tk);可以理解的是。
按照下式对叠后地震数据d(x,ti)和叠后多次波数据f(x,tj)进行互相关,得到互相关数据r(x,tk):
其中,k=1,2,…,n。
需要说明的是,互相关处理得到的对应的互相关数据即为叠后地震数据d(x,ti)和叠后多次波数据f(x,tj)之间的互相关数据。
S104:对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;
在本步骤中,该预设范围选取为海底水平或者起伏平缓的一段,海底水平或者起伏平缓的一段,就是在叠后数据上选择海底水平或者起伏平缓的部分。通过该段内的有关的数据确定的激发延迟时间,能够提高激发延迟时间的精准度。
因此确定海底水平或者起伏平缓的一段后并确定该段的起始位置xmin和终止位置xmax,按照下式将选取的预设范围内的互相关数据r(x,tk)叠加在一起,得到道数据D(ti):
其中,k=1,2,…,n。
通过该道数据D(ti)来对地震数据激发延迟时间进行校正,能够提高激发延迟时间的精准度。
需要说明的是,叠加处理得到的对应的道数据即为对确定的互相关数据进行叠加处理得到的道数据。
S105:基于所述道数据中最大值所对应的时间对地震数据激发延迟时间进行校正。
在本步骤中,根据步骤S104中确定的道数据D(ti),确定该道数据D(ti)中的最大值以及该最大值所对应的时间t0,该最大值所对应的时间t0即为检波器开始记录时间与震源的激发时间之间的差值。通过该最大值所对应的时间t0对地震数据激发延迟时间进行校正,能够有效提高激发延迟时间的精确性,进而对拖缆地震数据进行校正,提高拖缆地震数据精度以及提高基于拖缆地震数据的成像质量。
进一步的,地震数据的记录时间中减去道数据中最大值所对应的时间,完成对地震数据激发延迟时间的校正。
从上述描述可知,本发明提供一种地震数据激发延迟时间的校正方法,通过对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;基于所述道数据中最大值所对应的时间对地震数据激发延迟时间进行校正,能够有效提高激发延迟时间的精确性,为后续地震数据处理中提供更加精确的地震波实际旅行时间,进而对拖缆地震数据进行校正,提高拖缆地震数据精度以及提高基于拖缆地震数据的成像质量。
从软件层面来看,为了能够有效提高激发延迟时间的精确性,为后续地震数据处理中提供更加精确的地震波实际旅行时间并对拖缆地震数据进行校正。本发明提供一种能够实现所述地震数据激发延迟时间的校正方法中全部内容的地震数据激发延迟时间的校正装置的具体实施方式,参见图4,所述地震数据激发延迟时间的校正装置具体包括如下内容:
动校正单元10,用于对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;
褶积单元20,用于根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;
处理单元30,用于将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;
叠加单元40,用于对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;
校正单元50,用于基于所述道数据中最大值所对应的时间对地震数据的激发延迟时间进行校正。
进一步的,所述动校正单元10包括:
动校正子单元,用于抽取所述地震数据中的共中心点道集,采用第一速度对该共中心点道集进行动校正;
叠加处理子单元,用于对所述共中心点道集进行叠加处理得到叠后地震数据。
进一步的,所述第一速度为地震波在海水中的传播速度。
进一步的,所述校正单元50包括:
校正子单元,用于地震数据的记录时间中减去所述道数据中最大值所对应的时间,完成对地震数据激发延迟时间的校正。
本发明提供的地震数据激发延迟时间的校正装置的实施例具体可以用于执行上述实施例中的地震数据激发延迟时间的校正方法的实施例的处理流程,其功能在此不再赘述,可以参照上述方法实施例的详细描述。
从上述描述可知,本发明提供一种地震数据激发延迟时间的校正装置,通过对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;基于所述道数据中最大值所对应的时间对地震数据激发延迟时间进行校正,能够有效提高激发延迟时间的精确性,为后续地震数据处理中提供更加精确的地震波实际旅行时间,进而对拖缆地震数据进行校正,提高拖缆地震数据精度以及提高基于拖缆地震数据的成像质量。
本发明的实施例还提供能够实现上述实施例中的地震数据激发延迟时间的校正方法中全部步骤的一种电子设备的具体实施方式,参见图5,所述电子设备具体包括如下内容:
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604;
其中,所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信;所述处理器601用于调用所述存储器602中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的地震数据激发延迟时间的校正方法中的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:通过对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;基于所述道数据中最大值所对应的时间对地震数据激发延迟时间进行校正。
本发明的实施例还提供能够实现上述实施例中的地震数据激发延迟时间的校正方法中全部步骤的一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的地震数据激发延迟时间的校正方法的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:通过对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;基于所述道数据中最大值所对应的时间对地震数据激发延迟时间进行校正。
基于上述内容,本发明还给出一种应用所述地震数据中强反射信号分离装置、电子设备或存储介质实现的地震数据中强反射信号分离方法的具体应用实例,具体包含有如下内容:
(1)采集海上地震勘探数据,取得叠前地震数据d0(x0,ti);
如图6所示,是没有做激发延迟时间校正的炮记录。该数据野外施工设计的激发延迟时间为25ms。
(2)抽取共中心点道集,做动校正后进行叠加,得到叠后地震数据d(x,ti);
如图7所示,是基于地震波在海水中的传播速度做动校正和叠加后的地震数据。
(3)对叠后地震数据d(x,ti)进行自褶积运算,得到预测的叠后多次波数据f(x,tj):
如图8所示,是由叠加后的地震数据自褶积预测的多次波。
(4)对叠后地震数据d(x,ti)和多次波数据f(x,tj)进行互相关,得到互相关数据r(x,tk)。
如图9所示,是多次波数据与叠后地震数据相关后的数据。
(5)从叠后地震数据d(x,ti)中选取海底水平或者起伏平缓的一段,把选取的海底水平或者起伏平缓的一段所对应的互相关数据r(x,tk)叠加在一起,得到对应的道数据,然后该道数据的最大值,并确定最大值对应的时间。
如图10所示,是互相关数据的叠加结果,从中可以找到其最大值对应的时间为-25.8ms,则激发延迟时间为25.8ms,这与野外施工设计的激发延迟时间25ms接近,但更精确。
(6)按照步骤(5)确定的时间(25.8ms),将叠前地震数据d0(x0,ti)的记录时间中减去激发延迟时间(25.8ms),完成激发延迟时间校正。
如图11所示,是激发延迟时间校正后的炮记录。
从图6和图11的对比中可以看出,经过校正后,地震数据的旅行时减少,激发延迟引起的旅行时误差被弥补,激发延迟效应得到了有效消除。
从上述描述可知,根据地震叠加数据的自褶积得到叠后预测多次波,将叠后多次波与地震叠加数据的互相关,确定互相关的最大值对应的时间,互相关最大值对应的时间就是激发延迟时间,然后利用求得的激发延迟时间校正地震数据,来消除激发延迟引起的地震波旅行时误差,得到符合满足地震波传播规律的地震数据,为后续地震处理提供旅行时更加精确的数据,从而提高成像质量。
虽然本发明提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
本领域技术人员应明白,本说明书的实施例可提供为方法、系统或计算机程序产品。因此,本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
本发明的说明书中,说明了大量具体细节。然而能够理解的是,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。类似地,应当理解,为了精简本发明公开并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释呈反映如下意图:即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。本发明并不局限于任何单一的方面,也不局限于任何单一的实施例,也不局限于这些方面和/或实施例的任意组合和/或置换。而且,可以单独使用本发明的每个方面和/或实施例或者与一个或更多其他方面和/或其实施例结合使用。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围,其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (10)
1.一种地震数据激发延迟时间的校正方法,其特征在于,包括:
对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;
根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;
将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;
对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;
基于所述道数据中最大值所对应的时间对地震数据的激发延迟时间进行校正。
2.根据权利要求1所述的地震数据激发延迟时间的校正方法,其特征在于,所述对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据,包括:
抽取所述地震数据中的共中心点道集,采用第一速度对该共中心点道集进行动校正;
对该动校正处理之后的共中心点道集进行叠加处理得到叠后地震数据。
3.根据权利要求2所述的地震数据激发延迟时间的校正方法,其特征在于,所述第一速度为地震波在海水中的传播速度。
4.根据权利要求1所述的地震数据激发延迟时间的校正方法,其特征在于,所述基于所述道数据中最大值所对应的时间对地震数据的激发延迟时间进行校正,包括:
地震数据的记录时间中减去所述道数据中最大值所对应的时间,完成对地震数据激发延迟时间的校正。
5.一种地震数据激发延迟时间的校正装置,其特征在于,包括:
动校正单元,用于对地震数据进行动校正处理和叠加处理得到对应的叠后地震数据;
褶积单元,用于根据所述叠后地震数据进行自褶积处理得到对应的叠后多次波数据;
处理单元,用于将所述叠后地震数据和所述叠后多次波数据进行互相关处理得到对应的互相关数据;
叠加单元,用于对预设范围内的所述互相关数据进行叠加处理,得到对应的道数据;
校正单元,用于基于所述道数据中最大值所对应的时间对地震数据的激发延迟时间进行校正。
6.根据权利要求5所述的地震数据激发延迟时间的校正装置,其特征在于,所述动校正单元包括:
动校正子单元,用于抽取所述地震数据中的共中心点道集,采用第一速度对该共中心点道集进行动校正;
叠加处理子单元,用于对该动校正处理之后的共中心点道集进行叠加处理得到叠后地震数据。
7.根据权利要求6所述的地震数据激发延迟时间的校正装置,其特征在于,所述第一速度为地震波在海水中的传播速度。
8.根据权利要求5所述的地震数据激发延迟时间的校正装置,其特征在于,所述校正单元包括:
校正子单元,用于地震数据的记录时间中减去所述道数据中最大值所对应的时间,完成对地震数据激发延迟时间的校正。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至4任一项所述的地震数据激发延迟时间的校正方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至4任一项所述的地震数据激发延迟时间的校正方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910634891.9A CN112230285B (zh) | 2019-07-15 | 2019-07-15 | 地震数据激发延迟时间的校正方法及校正装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910634891.9A CN112230285B (zh) | 2019-07-15 | 2019-07-15 | 地震数据激发延迟时间的校正方法及校正装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112230285A CN112230285A (zh) | 2021-01-15 |
CN112230285B true CN112230285B (zh) | 2023-09-26 |
Family
ID=74111123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910634891.9A Active CN112230285B (zh) | 2019-07-15 | 2019-07-15 | 地震数据激发延迟时间的校正方法及校正装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112230285B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201409A (zh) * | 2006-12-14 | 2008-06-18 | 中国石油天然气集团公司 | 一种地震数据变相位校正方法 |
CN101813786A (zh) * | 2010-04-02 | 2010-08-25 | 中国石油集团西北地质研究所 | 子波处理二步法反褶积方法 |
CN102636811A (zh) * | 2012-04-10 | 2012-08-15 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种海上二维地震资料中多次波的消除方法 |
GB201416562D0 (en) * | 2013-09-19 | 2014-11-05 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
CN104459794A (zh) * | 2014-12-01 | 2015-03-25 | 郑鸿明 | 共反射点道集时变时间差值的校正方法及装置 |
US9075162B2 (en) * | 2011-11-10 | 2015-07-07 | Pgs Geophysical As | Method and system for separating seismic sources in marine simultaneous shooting acquisition |
CN105954798A (zh) * | 2016-04-26 | 2016-09-21 | 中国石油天然气集团公司 | 一种确定起伏地表的叠前时间偏移速度的方法和装置 |
US9671511B2 (en) * | 2012-08-31 | 2017-06-06 | Cgg Services Sas | Horizontal streamer broadband marine seismic acquisition configuration and processing |
WO2017106127A1 (en) * | 2015-12-18 | 2017-06-22 | Saudi Arabian Oil Company | Automated near surface analysis by surface-consistent refraction methods |
CN106932824A (zh) * | 2017-03-24 | 2017-07-07 | 北京大学 | 陆地地震勘探资料的降维自适应层间多次波压制方法 |
CN107193045A (zh) * | 2017-07-07 | 2017-09-22 | 中国石油天然气集团公司 | 一种地震数据处理方法及装置 |
WO2017167191A1 (zh) * | 2016-04-01 | 2017-10-05 | 中国石油天然气股份有限公司 | 地震数据处理方法和装置 |
CN107942386A (zh) * | 2017-11-15 | 2018-04-20 | 河南理工大学 | 一种计算槽波地震数据激发延迟时间的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3260887A1 (en) * | 2016-06-23 | 2017-12-27 | CGG Services SAS | Methods and data processing apparatus for seismic signal separation |
-
2019
- 2019-07-15 CN CN201910634891.9A patent/CN112230285B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201409A (zh) * | 2006-12-14 | 2008-06-18 | 中国石油天然气集团公司 | 一种地震数据变相位校正方法 |
CN101813786A (zh) * | 2010-04-02 | 2010-08-25 | 中国石油集团西北地质研究所 | 子波处理二步法反褶积方法 |
US9075162B2 (en) * | 2011-11-10 | 2015-07-07 | Pgs Geophysical As | Method and system for separating seismic sources in marine simultaneous shooting acquisition |
CN102636811A (zh) * | 2012-04-10 | 2012-08-15 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种海上二维地震资料中多次波的消除方法 |
US9671511B2 (en) * | 2012-08-31 | 2017-06-06 | Cgg Services Sas | Horizontal streamer broadband marine seismic acquisition configuration and processing |
GB201416562D0 (en) * | 2013-09-19 | 2014-11-05 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
CN104459794A (zh) * | 2014-12-01 | 2015-03-25 | 郑鸿明 | 共反射点道集时变时间差值的校正方法及装置 |
WO2017106127A1 (en) * | 2015-12-18 | 2017-06-22 | Saudi Arabian Oil Company | Automated near surface analysis by surface-consistent refraction methods |
WO2017167191A1 (zh) * | 2016-04-01 | 2017-10-05 | 中国石油天然气股份有限公司 | 地震数据处理方法和装置 |
CN105954798A (zh) * | 2016-04-26 | 2016-09-21 | 中国石油天然气集团公司 | 一种确定起伏地表的叠前时间偏移速度的方法和装置 |
CN106932824A (zh) * | 2017-03-24 | 2017-07-07 | 北京大学 | 陆地地震勘探资料的降维自适应层间多次波压制方法 |
CN107193045A (zh) * | 2017-07-07 | 2017-09-22 | 中国石油天然气集团公司 | 一种地震数据处理方法及装置 |
CN107942386A (zh) * | 2017-11-15 | 2018-04-20 | 河南理工大学 | 一种计算槽波地震数据激发延迟时间的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112230285A (zh) | 2021-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2974636B1 (fr) | Dispositif et procede d'extrapolation angulaire adaptative pour des rassemblements angulaires tridimensionnels de migration inverse dans le temps | |
CN101630019B (zh) | 一种校正预测的海底多次波大时差的方法 | |
FR2787201A1 (fr) | Methode et dispositif d'acquisition synchronisee de signaux sismiques | |
WO2016156530A1 (en) | Method and system of processing seismic data by providing surface aperture common image gathers | |
WO2013093063A4 (en) | Device and method for deghosting variable depth streamer data | |
CN111221037A (zh) | 解耦弹性逆时偏移成像方法和装置 | |
Barak et al. | Six-component seismic land data acquired with geophones and rotation sensors: Wave-mode selectivity by application of multicomponent polarization filtering | |
MX2013005640A (es) | Metodo y aparato para la eliminacion de fantasma de pre-pila de datos sismicos. | |
US20160054465A1 (en) | Seismic data spectrum restoring and broadening | |
CN111781635B (zh) | 海底四分量弹性波高斯束深度偏移方法和装置 | |
FR2974187B1 (fr) | Dispositif et procede pour calculer la migration reverse time dans des milieux orthorhombiques inclines en 3d | |
CN112230285B (zh) | 地震数据激发延迟时间的校正方法及校正装置 | |
CN110907989A (zh) | 重建拟地面地震反射波成像方法及系统 | |
CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
CN113721296A (zh) | 远震数据处理方法及装置 | |
CN111337972B (zh) | 用于地震数据重建的规则化域选取方法、系统 | |
CN111781644B (zh) | 浅中地层地震数据的线性干扰衰减方法及装置 | |
AU2011248910B2 (en) | System and method for compensating time and offset varying near-surface effects in seismic data background | |
Liu et al. | Multiple attenuation using multichannel predictive deconvolution in radial domain | |
CN112180448A (zh) | 剩余静校正方法、装置和电子设备 | |
CN111551993B (zh) | 压制鸣震的方法及装置 | |
FR2966253A1 (fr) | Proceder et dispositif pour retirer des fantomes de donnees sismiques | |
CN111610562B (zh) | 地震数据中强反射信号分离方法及装置 | |
CN113031071B (zh) | 地震波长波长静校正方法及装置 | |
US20150219778A1 (en) | Method for Obtaining an Adapted Sequence of Shot Predictions From a Raw Sequence of Shot Predictions |
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 |