CN111965703A - 井中地震数据中单炮记录的下行波反褶积处理方法及系统 - Google Patents
井中地震数据中单炮记录的下行波反褶积处理方法及系统 Download PDFInfo
- Publication number
- CN111965703A CN111965703A CN201910418079.2A CN201910418079A CN111965703A CN 111965703 A CN111965703 A CN 111965703A CN 201910418079 A CN201910418079 A CN 201910418079A CN 111965703 A CN111965703 A CN 111965703A
- Authority
- CN
- China
- Prior art keywords
- downlink
- wavelet
- record
- single shot
- wave
- 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
Links
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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示出本申请实施例中一种井中地震数据中单炮记录的下行波反褶积处理方法的流程示意图。
图2a示出本申请实施例中井中地震数据单炮三分量记录。
图2b示出本申请实施例中井中地震数据P波初至拉平记录。
图2c示出本申请实施例中井中地震数据P波初至拉平记录频谱图。
图2d示出本申请实施例中反褶积前后的井中地震数据波场记录对比(左:反褶积前;中:下行子波长度为250毫秒的反褶积;右:下行子波长度为250毫秒的反褶积)。
图2e示出本申请实施例中反褶积后的井中地震数据频谱图。
图2f示出本申请实施例中反褶积后的井中地震数据走廊叠加嵌入地面地震后对比图(走廊叠加记录位置在横坐标道号为480处)。
图3示出本申请实施例中图1中步骤S200的具体流程示意图。
图4示出本申请实施例中图3的步骤S240的具体流程示意图。
图5示出本申请实施例中图3中步骤S300的具体流程示意图。
图6示出本申请实施例中一种井中地震数据中单炮记录的下行波反褶积处理系统结构示意图。
图7示出本申请实施例中图6中波场分离模块200的具体结构示意图。
图8示出本申请实施例中图7中下行子波获取单元240的具体结构示意图。
图9示出本申请实施例中图6中反褶积处理模块300的具体结构示意图。
图10示出适于用来实现本申请实施例的电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本申请主要基于井中地震数据下行子波随深度变化的假设,通过波场分离估算每个深度的下行子波,并设计合理的期望输出子波,达到最佳的反褶积效果,在实际生产应用中,对于压制井中地震数据的近地表多次波,提高数据的分辨率和一致性具有重要意义。
下面结合附图对本申请实施例进行详细说明。
本申请第一方面实施例提供一种井中地震数据中单炮记录的下行波反褶积处理方法,如图1所示,包括:
S100:获取井中地震波场数据中的任意一个单炮记录;所述单炮记录是通过设于井筒内的多个传感器采集震源单次激发的地震波得到,每个传感器处于不同深度;
S200:对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波;
S300:基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果。
本方面提供的井中地震数据中单炮记录的下行波反褶积处理方法,首先对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波,然后基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果,通过本方法的处理后井中地震数据近地表多次波可以得到很大程度的压制,资料分辨率明显提升,为后续的地震数据处理解释提供了保障。
下面具体对上述步骤进行详细说明。
在步骤S100中,首先利用三分量或单分量传感器下入井中接收,在地面或井中利用爆炸震源或人工可控震源激发,采集得到井中地震波场数据。
本申请中,所述井中地震波场数据,是指激发点在地表或井中布设、传感器在井中布设的观测方式下所采集得到的井中地震波场数据。
井中地震波场数据具体由多个单炮记录构成,单炮记录即为源单次激发的地震波被传感器记录的数据。
可以理解,本申请可以通过单分量传感器获得上述井中地震波场数据中的单炮记录,也可以通过三分量传感器获得,本申请对此不做限制,当本申请实施例中采用的是单分量传感器时,其方向固定,不会发生随机旋转,因此不需要确定P波的方向,其径向波场分量即为单分量传感器记录的分量,但对于三分量传感器,由于三分量传感器记录地震数据时会产生随机旋转,因此应对同一个深度的三分量按初至时窗振幅最大化且相位一致原则进行三分量旋转合成,得到每个传感器纵波初至能量最大化的径向分量。
一实施例中,第一下行子波和第二下行子波可以通过对径向波场分量进行归一化校正处理得到,具体的,如图3所示,步骤S200包括:
S210:从所述单炮记录中提取每个传感器记录的所有径向波场分量;
S220:分别对每个传感器对应的所有径向波场分量进行归一化校正处理,得到与每个传感器一一对应并且振幅统一的波形记录;
S230:去除每个所述波形记录中的非下行波,得到一一对应的下行波记录;
S240:基于每个下行波记录,获得对应每个传感器所处深度的第一下行子波和第二下行子波。
步骤S220中,若为三分量数据,按径向分量P初至时窗的均方根振幅对同一单炮记录的全部径向分量道进行归一化处理得到归一化后的记录W1i,并记下归一化校正因子Ai;若为单分量数据,则按单分量初至时窗的均方根振幅对同一单炮记录的全部道进行归一化处理得到归一化后的记录W1i,并记下归一化校正因子Ai。i表示第i个传感器的序号。
在优选实施例中,初至时窗长度的选择按初至波形的主周期长度确定,一般为1-3倍主周期长度。
此外,步骤S230可以通过中值滤波、压制等方式处理,以中值滤波作为示例,对归一化后的记录W1i沿水平方向进行中值滤波,得到拉平后的下行波记录W2i。
基于下行波记录得到两个下行子波,可以通过分级截取的方式,截取出时间长度不相同的两个下行子波,具体的,如图4所示,步骤S240包括:
S241:根据多次波的发育程度,对所述下行波记录执行截取操作,获得时长为初至波第一个周期时长的第一倍数的第一下行子波记录;
S242:对每个第一下行子波记录的最后若干样点进行镶边处理,使镶边处理的样点对应的样点值沿朝向每个第一下行子波记录尾部方向逐步减小至零,生成对应每个传感器的第一下行子波;
S243:对每个第一下行子波记录执行截取操作,获得时长为初至波第一个周期时长的第二倍数的第二下行子波记录;
S244:对每个第二下行子波记录的最后若干样点进行镶边处理,使镶边处理的样点对应的样点值沿朝向每个第二下行子波记录尾部方向逐渐减小至零,生成对应每个传感器的第二下行子波。
该实施例中,根据多次波的发育程度,若多次波数量较多(发育较好),对应地第一倍数即相应更大,一般地,第一倍数可以选取3-10,例如一具体实施例中,根据初至波形的主周期和近地表多次波程度对下行波截取Lw长度的数据记录并对最后10个样点值进行镶边处理得到下行随深度变化的第一下行子波Wi,即通过斜坡函数使其值逐步趋向于零。
由于执行截取操作后,得到的是不完整的下行子波,这时如果不进行镶边处理,无法输出正确的结果(输入需要是完整的子波),因此镶边处理即相当于使得每一个被截取的不完整波可以等效于完整的子波。
同样的理由可以适用于获得第二下行子波,一实施例中,第二倍数相较于第一倍数更小(第二下行子波通过截取第一下行子波记录得到),对W2i执行截取操作,得到的下行子波截取Ld长度的数据记录并对最后10个样点值进行镶边处理得到下行随深度变化的第二下行子波Wdi,第二下行子波Wdi又称为期望输出子波。一般Ld长度等于初至波形第一个主周期长度的1.2-1.5倍。即第二倍数为1.2-1.5。
进一步的,在一优选实施例中,上述方法还包括:
S001:对井中地震波场数据中的单炮记录进行预处理。
具体的,步骤S001包括如下至少之一:
S001a:去除所述单炮记录中的随机干扰;
S001b:设置所述单炮记录中的炮点接收点坐标。
预处理步骤的目的在于减少后续的频谱分析的处理量以及误差,经过炮点接收点坐标设置、去除随机干扰等处理,可以减少后续的数据处理量,减少处理负担。
预处理步骤还包括初至波拾取,即在对所述单炮记录进行波场分离之前,还包括:
S001c:获取所述单炮记录中对应每个传感器的初至波到达时间点;
S001d:根据初至波的到达每个传感器的时间点对所述单炮记录进行拉平处理,使所述单炮记录对应的下行波排列整齐。
本实施例中所述拉平处理,若初至时间值不能被采样率整除时,应通过数据内插确保时移准确性。
此外,如图5所示,步骤S300具体包括:
S301:根据每个第一下行子波和对应的第二下行子波,生成对应每个传感器所处深度的反褶积因子;
S302:根据每个所述反褶积因子对所述单炮记录的下行波进行反褶积处理,得到下行波随深度变化的反褶积结果。
本申请中的反褶积处理基于现有的反褶积技术,可以理解,本申请的发明构思在于对单炮记录进行波场分离,得到多个对应不同深度的下行子波,进而进行反褶积处理后能够得到考虑深度影响的处理结果,因此本申请对具体的反褶积处理方法不做限制。
对于每一个第一下行子波Wi和对应的第二下行子波Wdi,设计一对应的反褶积因子(又称反滤波因子)Fi。作为一种现有技术中的反褶积技术,反褶积因子可以通过最小平方反褶积技术、两步法预测误差反褶积等获得,本申请在此不做赘述。
下面通过一具体场景进行详细说明。
一具体场景中,首先利用三分量井中检波器下入井中接收,在地面或井中利用人工可控震源激发,采集得到井中地震波场数据。
然后对采集数据进行预处理,经过炮点接收点坐标设置、去除随机干扰、拾取纵波初至等三步处理,得到预处理后数据(如图2a)。
利用所得到的预处理后数据根据拾取的纵波初至时间进行拉平,即按初至时间值tp向上时移,使下行波排齐(如图2b)。从其频谱可以看出下行波频宽有明显的随深度变化且受多次波影响频谱有多个峰值(如图2c)。
对同一个深度的三分量按初至时窗振幅最大化且相位一致原则进行三分量旋转合成,得到每个传感器纵波初至能量最大化的径向分量P。
按径向分量P初至时窗的均方根振幅对同一单炮记录的全部径向分量道进行归一化处理得到归一化后的记录W1i,并记下归一化校正因子Ai。i表示第i个传感器的序号。
这里初至时窗长度的选择按初至波形的1个主周期长度确定。
对归一化后的记录W1i沿水平方向进行中值滤波,得到拉平后的下行波记录W2i,根据初至波形的主周期和近地表多次波程度对下行波截取Lw长度的数据记录并对最后10个样点值进行镶边处理得到下行随深度变化的子波Wi,即通过斜坡函数使其值逐步趋向于零。这里下行子波长度Lw等于初至波形第一个主周期长度的约8倍,即300毫秒。
对下行子波截取Ld长度的数据记录并对最后10个样点值进行镶边处理得到下行随深度变化的期望输出子波Wdi。这里期望输出子波长度Ld长度为50毫秒。
以下行子波Wi和期望输出子波Wdi设计反褶积因子Fi,利用常规反褶积方法实现对输入波场的随深度变化反褶积处理。
选择较浅层的单个期望输出子波Wd作为所有传感器记录反褶积的期望输出可以实现Q吸收的补偿,在压制近地表多次波的同时提高深层波场的分辨率(如图2d),频谱展平变宽(如图2e)。得到的走廊叠加记录与地面地震的镶嵌对比表明反褶积效果良好(如图2f)
可以知晓,本方面提供的井中地震数据中单炮记录的下行波反褶积处理方法,首先对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波,然后基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果,通过本方法的处理后井中地震数据近地表多次波可以得到很大程度的压制,资料分辨率明显提升,为后续的地震数据处理解释提供了保障。
基于相同的发明构思,本申请第二方面实施例提供一种井中地震数据中单炮记录的下行波反褶积处理系统,如图6所示,包括:
获取模块100,获取井中地震波场数据中的任意一个单炮记录;所述单炮记录是通过设于井筒内的多个传感器采集震源单次激发的地震波得到,每个传感器处于不同深度;
波场分离模块200,对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波;
反褶积处理模块300,基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果。
本方面提供的井中地震数据中单炮记录的下行波反褶积处理系统,首先对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波,然后基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果,通过本方法的处理后井中地震数据近地表多次波可以得到很大程度的压制,资料分辨率明显提升,为后续的地震数据处理解释提供了保障。
基于与上一方面相同的理由,井中地震波场数据具体由多个单炮记录构成,单炮记录即为源单次激发的地震波被传感器记录的数据。
基于与上一方面相同的理由,一实施例中,第一下行子波和第二下行子波可以通过对径向波场分量进行归一化校正处理得到,具体的,如图7所示,所述波场分离模块200,包括:
径向波场分离提取单元210,从所述单炮记录中提取每个传感器记录的所有径向波场分量;
归一化校正单元220,分别对每个传感器对应的所有径向波场分量进行归一化校正处理,得到与每个传感器一一对应并且振幅统一的波形记录;
非下行波去除单元230,去除每个所述波形记录中的非下行波,得到一一对应的下行波记录;
下行子波获取单元240,基于每个下行波记录,获得对应每个传感器所处深度的第一下行子波和第二下行子波。
基于与上一方面相同的理由,基于下行波记录得到两个下行子波,可以通过分级截取的方式,截取出时间长度不相同的两个下行子波,具体的,如图8所示,所述下行子波获取单元240,包括:
第一下行子波记录获取单元241,根据多次波的发育程度,对所述下行波记录执行截取操作,获得时长为初至波第一个周期时长的第一倍数的第一下行子波记录;
第一下行子波生成单元242,对每个第一下行子波记录的最后若干样点进行镶边处理,使镶边处理的样点对应的样点值沿朝向每个第一下行子波记录尾部方向逐步减小至零,生成对应每个传感器的第一下行子波;
第二下行子波记录获取单元243,对每个第一下行子波记录执行截取操作,获得时长为初至波第一个周期时长的第二倍数的第二下行子波记录;
第二下行子波生成单元244,对每个第二下行子波记录的最后若干样点进行镶边处理,使镶边处理的样点对应的样点值沿朝向每个第二下行子波记录尾部方向逐渐减小至零,生成对应每个传感器的第二下行子波。
在一优选实施例中,本系统还包括:
预处理模块,对井中地震波场数据中的单炮记录进行预处理。
具体的,所述预处理模块包括如下至少一个:
随机干扰去除单元,去除所述单炮记录中的随机干扰;
接收点坐标设置单元,设置所述单炮记录中的炮点接收点坐标。
预处理步骤的目的在于减少后续的频谱分析的处理量以及误差,经过炮点接收点坐标设置、去除随机干扰等处理,可以减少后续的数据处理量,减少处理负担。
一更优选实施例中,本系统还包括:
初至波拾取模块,获取所述单炮记录中对应每个传感器的初至波到达时间点;
单炮记录拉平模块,根据初至波的到达每个传感器的时间点对所述单炮记录进行拉平处理,使所述单炮记录对应的下行波排列整齐。
本实施例中所述拉平处理,若初至时间值不能被采样率整除时,应通过数据内插确保时移准确性。
进一步的,基于同样的理由,如图9所示,所述反褶积处理模块,包括:
反褶积因子生成单元301,根据每个第一下行子波和对应的第二下行子波,生成对应每个传感器所处深度的反褶积因子;
反褶积处理单元302,根据每个所述反褶积因子对所述单炮记录的下行波进行反褶积处理,得到下行波随深度变化的反褶积结果。
可以知晓,本方面提供的井中地震数据中单炮记录的下行波反褶积处理系统,首先对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波,然后基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果,通过本方法的处理后井中地震数据近地表多次波可以得到很大程度的压制,资料分辨率明显提升,为后续的地震数据处理解释提供了保障。
本申请的实施例还提供能够实现上述实施例中的方法中全部步骤的一种电子设备的具体实施方式,参见图10,所述电子设备具体包括如下内容:
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604;
其中,所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信;
所述处理器601用于调用所述存储器602中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的方法中的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
S100:获取井中地震波场数据中的任意一个单炮记录;所述单炮记录是通过设于井筒内的多个传感器采集震源单次激发的地震波得到,每个传感器处于不同深度;
S200:对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波;
S300:基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果。
从上述描述可知,本申请提供的电子设备处理后井中地震数据近地表多次波可以得到很大程度的压制,资料分辨率明显提升,为后续的地震数据处理解释提供了保障。
本申请的实施例还提供能够实现上述实施例中的方法中全部步骤的一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的方法的全部步骤,例如,所述处理器执行所述计算机程序时实现下述步骤:
S100:获取井中地震波场数据中的任意一个单炮记录;所述单炮记录是通过设于井筒内的多个传感器采集震源单次激发的地震波得到,每个传感器处于不同深度;
S200:对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波;
S300:基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果。
从上述描述可知,本申请提供的计算机可读存储介质处理后井中地震数据近地表多次波可以得到很大程度的压制,资料分辨率明显提升,为后续的地震数据处理解释提供了保障。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于硬件+程序类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或终端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。本领域技术人员应明白,本说明书的实施例可提供为方法、系统或计算机程序产品。因此,本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。以上所述仅为本说明书实施例的实施例而已,并不用于限制本说明书实施例。对于本领域技术人员来说,本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本说明书实施例的权利要求范围之内。
Claims (16)
1.一种井中地震数据中单炮记录的下行波反褶积处理方法,其特征在于,包括:
获取井中地震波场数据中的任意一个单炮记录;所述单炮记录是通过设于井筒内的多个传感器采集震源单次激发的地震波得到,每个传感器处于不同深度;
对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波;
基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果。
2.根据权利要求1所述的下行波反褶积处理方法,其特征在于,所述基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,包括:
根据每个第一下行子波和对应的第二下行子波,生成对应每个传感器所处深度的反褶积因子;
根据每个所述反褶积因子对所述单炮记录的下行波进行反褶积处理,得到下行波随深度变化的反褶积结果。
3.根据权利要求1所述的下行波反褶积处理方法,其特征在于,所述对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波,包括:
从所述单炮记录中提取每个传感器记录的所有径向波场分量;
分别对每个传感器对应的所有径向波场分量进行归一化校正处理,得到与每个传感器一一对应并且振幅统一的波形记录;
去除每个所述波形记录中的非下行波,得到一一对应的下行波记录;
基于每个下行波记录,获得对应每个传感器所处深度的第一下行子波和第二下行子波。
4.根据权利要求3所述的下行波反褶积处理方法,其特征在于,所述基于每个下行波记录,获得对应每个传感器所处深度的第一下行子波和第二下行子波,包括:
根据多次波的发育程度,对所述下行波记录执行截取操作,获得时长为初至波第一个周期时长的第一倍数的第一下行子波记录;
对每个第一下行子波记录的最后若干样点进行镶边处理,使镶边处理的样点对应的样点值沿朝向每个第一下行子波记录尾部方向逐步减小至零,生成对应每个传感器的第一下行子波;
对每个第一下行子波记录执行截取操作,获得时长为初至波第一个周期时长的第二倍数的第二下行子波记录;
对每个第二下行子波记录的最后若干样点进行镶边处理,使镶边处理的样点对应的样点值沿朝向每个第二下行子波记录尾部方向逐渐减小至零,生成对应每个传感器的第二下行子波。
5.根据权利要求1所述的下行波反褶积处理方法,其特征在于,在对所述单炮记录进行波场分离之前,还包括:
对井中地震波场数据中的单炮记录进行预处理。
6.根据权利要求5所述的下行波反褶积处理方法,其特征在于,所述井中地震波场数据中的单炮记录进行预处理,包括如下至少之一:
去除所述单炮记录中的随机干扰;
设置所述单炮记录中的炮点接收点坐标。
7.根据权利要求1所述的下行波反褶积处理方法,其特征在于,在对所述单炮记录进行波场分离之前,还包括:
获取所述单炮记录中对应每个传感器的初至波到达时间点;
根据初至波的到达每个传感器的时间点对所述单炮记录进行拉平处理,使所述单炮记录对应的下行波排列整齐。
8.一种井中地震数据中单炮记录的下行波反褶积处理系统,其特征在于,包括:
获取模块,获取井中地震波场数据中的任意一个单炮记录;所述单炮记录是通过设于井筒内的多个传感器采集震源单次激发的地震波得到,每个传感器处于不同深度;
波场分离模块,对所述单炮记录进行波场分离,生成对应每个传感器所处深度的第一下行子波和第二下行子波;
反褶积处理模块,基于每个第一下行子波和对应的第二下行子波,对所述单炮记录的下行波进行反褶积处理,生成下行波随深度变化的反褶积结果。
9.根据权利要求8所述的下行波反褶积处理系统,其特征在于,所述反褶积处理模块,包括:
反褶积因子生成单元,根据每个第一下行子波和对应的第二下行子波,生成对应每个传感器所处深度的反褶积因子;
反褶积处理单元,根据每个所述反褶积因子对所述单炮记录的下行波进行反褶积处理,得到下行波随深度变化的反褶积结果。
10.根据权利要求8所述的下行波反褶积处理系统,其特征在于,所述波场分离模块,包括:
径向波场分离提取单元,从所述单炮记录中提取每个传感器记录的所有径向波场分量;
归一化校正单元,分别对每个传感器对应的所有径向波场分量进行归一化校正处理,得到与每个传感器一一对应并且振幅统一的波形记录;
非下行波去除单元,去除每个所述波形记录中的非下行波,得到一一对应的下行波记录;
下行子波获取单元,基于每个下行波记录,获得对应每个传感器所处深度的第一下行子波和第二下行子波。
11.根据权利要求10所述的下行波反褶积处理系统,其特征在于,所述下行子波获取单元,包括:
第一下行子波记录获取单元,根据多次波的发育程度,对所述下行波记录执行截取操作,获得时长为初至波第一个周期时长的第一倍数的第一下行子波记录;
第一下行子波生成单元,对每个第一下行子波记录的最后若干样点进行镶边处理,使镶边处理的样点对应的样点值沿朝向每个第一下行子波记录尾部方向逐步减小至零,生成对应每个传感器的第一下行子波;
第二下行子波记录获取单元,对每个第一下行子波记录执行截取操作,获得时长为初至波第一个周期时长的第二倍数的第二下行子波记录;
第二下行子波生成单元,对每个第二下行子波记录的最后若干样点进行镶边处理,使镶边处理的样点对应的样点值沿朝向每个第二下行子波记录尾部方向逐渐减小至零,生成对应每个传感器的第二下行子波。
12.根据权利要求8所述的下行波反褶积处理系统,其特征在于,还包括:
预处理模块,对井中地震波场数据中的单炮记录进行预处理。
13.根据权利要求12所述的下行波反褶积处理系统,其特征在于,所述预处理模块包括如下至少一个:
随机干扰去除单元,去除所述单炮记录中的随机干扰;
接收点坐标设置单元,设置所述单炮记录中的炮点接收点坐标。
14.根据权利要求8所述的下行波反褶积处理系统,其特征在于,还包括:
初至波拾取模块,获取所述单炮记录中对应每个传感器的初至波到达时间点;
单炮记录拉平模块,根据初至波的到达每个传感器的时间点对所述单炮记录进行拉平处理,使所述单炮记录对应的下行波排列整齐。
15.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至7任一项所述的井中地震数据中单炮记录的下行波反褶积处理方法的步骤。
16.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至7任一项所述的井中地震数据中单炮记录的下行波反褶积处理方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418079.2A CN111965703B (zh) | 2019-05-20 | 2019-05-20 | 井中地震数据中单炮记录的下行波反褶积处理方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418079.2A CN111965703B (zh) | 2019-05-20 | 2019-05-20 | 井中地震数据中单炮记录的下行波反褶积处理方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111965703A true CN111965703A (zh) | 2020-11-20 |
CN111965703B CN111965703B (zh) | 2023-04-07 |
Family
ID=73358280
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910418079.2A Active CN111965703B (zh) | 2019-05-20 | 2019-05-20 | 井中地震数据中单炮记录的下行波反褶积处理方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111965703B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040122596A1 (en) * | 2002-12-19 | 2004-06-24 | Core Laboratories, Inc. | Method for high frequency restoration of seismic data |
CN103460074A (zh) * | 2011-03-31 | 2013-12-18 | 埃克森美孚上游研究公司 | 全波场反演中小波估计和多次波预测的方法 |
CN103913770A (zh) * | 2014-03-27 | 2014-07-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于vsp资料对地震数据进行处理的方法 |
-
2019
- 2019-05-20 CN CN201910418079.2A patent/CN111965703B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040122596A1 (en) * | 2002-12-19 | 2004-06-24 | Core Laboratories, Inc. | Method for high frequency restoration of seismic data |
CN103460074A (zh) * | 2011-03-31 | 2013-12-18 | 埃克森美孚上游研究公司 | 全波场反演中小波估计和多次波预测的方法 |
CN103913770A (zh) * | 2014-03-27 | 2014-07-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于vsp资料对地震数据进行处理的方法 |
Non-Patent Citations (2)
Title |
---|
李彦鹏 等: "VSP处理新技术及应用", 《石油地球物理勘探》 * |
李文杰 等: "利用VSP提高地震资料处理质量的新途径", 《新疆石油地质》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111965703B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9103935B2 (en) | Method of first arrival picking of seismic refraction data | |
CN1275051C (zh) | 用于对级联扫描数据的噪声去除方法 | |
WO2012158456A1 (en) | Frequency-varying filtering of simultaneous source seismic data | |
CN107966728A (zh) | 炸药震源多震源系统的地震采集方法及装置 | |
CN103645497B (zh) | 一种基于出射角的多分量波场分离方法 | |
CN110907995B (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
CN112394414B (zh) | 两步法地震绕射波场叠前分离的方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN1295672A (zh) | 地滚衰减方法 | |
CN104133244A (zh) | 一种震源横波勘探的近地表静校正处理方法 | |
CN111965703B (zh) | 井中地震数据中单炮记录的下行波反褶积处理方法及系统 | |
CN110907989A (zh) | 重建拟地面地震反射波成像方法及系统 | |
CN111123361B (zh) | 垂直地震剖面地震数据规则化重建方法及装置、存储介质 | |
CN105301656B (zh) | 一种压制虚反射信号的方法及装置 | |
CN100349005C (zh) | 地震数据多炮组合叠前深度偏移方法 | |
CN112198547A (zh) | 深层或超深层地震资料处理方法及装置 | |
CN105093275B (zh) | 基于速度模型的有效信号提取iss波场分离方法 | |
CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
CN113721296A (zh) | 远震数据处理方法及装置 | |
CN102798887A (zh) | 一种提高可控震源记录初至质量的滤波方法 | |
CN113945978A (zh) | 一种深度域地震数据优化方法及装置 | |
CN113138419B (zh) | 提取下行子波和衰减参数的方法、装置 | |
CN111965701B (zh) | 近地表结构反演方法及系统 | |
CN112799132B (zh) | 微局部线性噪声压制方法及装置 | |
CN113917539B (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 |