CN111337982A - 绕射波叠前炮集的分离方法、装置和电子设备 - Google Patents
绕射波叠前炮集的分离方法、装置和电子设备 Download PDFInfo
- Publication number
- CN111337982A CN111337982A CN202010309275.9A CN202010309275A CN111337982A CN 111337982 A CN111337982 A CN 111337982A CN 202010309275 A CN202010309275 A CN 202010309275A CN 111337982 A CN111337982 A CN 111337982A
- Authority
- CN
- China
- Prior art keywords
- data
- gather
- shot
- filtering
- seismic
- 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 67
- 238000001914 filtration Methods 0.000 claims abstract description 92
- 238000001514 detection method Methods 0.000 claims abstract description 27
- 238000000926 separation method Methods 0.000 claims abstract description 24
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 6
- 239000000523 sample Substances 0.000 claims description 6
- 239000007983 Tris buffer Substances 0.000 claims description 3
- 230000002146 bilateral effect Effects 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 8
- 230000008569 process Effects 0.000 description 16
- 230000002159 abnormal effect Effects 0.000 description 8
- 238000004891 communication Methods 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- HWJPHQNEWARZLH-UHFFFAOYSA-N 1,1,2,2,3,3,4,4,5,5-decafluoro-6,6-bis(trifluoromethyl)cyclohexane Chemical compound FC(F)(F)C1(C(F)(F)F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C1(F)F HWJPHQNEWARZLH-UHFFFAOYSA-N 0.000 description 3
- 239000003245 coal Substances 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004432 carbon atom Chemical group C* 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 238000013507 mapping 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
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
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
技术领域
本发明涉及地震勘探技术领域,尤其是涉及一种绕射波叠前炮集的分离方法、装置和电子设备。
背景技术
在地质勘探中,对于地质不连续性的准确预测对于探测精度以及生产安全具有重要意义。在进行人工震源操作过程中,形成的地震波在地面传播时,当遇到地质异常区域会形成绕射波。在绕射波中包含该地质异常体的有效数据,可通过分析得到该异常体的详细属性,对于开采相关行业具有指导意义。
绕射波具有揭示小于地震波波长的小尺度物体的特性,然而由于绕射相对于镜面反射较弱,存在于绕射中的细节信息通常被镜面反射所掩盖。此外,在常规处理过程中,特别是在常规动校正和叠加过程中,绕射会遭到破坏。因此,现有技术中对于绕射波的分离方式,还有进一步提升空间。
发明内容
有鉴于此,本发明的目的在于提供一种绕射波叠前炮集的分离方法、装置和电子设备,利用反射波和绕射波的运动学差异构建共虚震源道集,通过滤波的方式去除共虚震源道集中的反射波,实现了对绕射波叠前炮集的有效分离,提升了叠前炮集中绕射波的分离效果。
第一方面,本发明实施例提供了一种绕射波叠前炮集的分离方法,该方法包括:
获取叠前炮集地震数据;叠前炮集地震数据中包含反射波数据以及绕射波数据;
将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集;
对共虚震源道集中的数据沿检波点方向进行滤波操作,滤除共虚震源道集中包含的反射波数据后得到滤波结果;
将滤波结果转换至炮集中,得到已完成分离的绕射波叠前炮集。
在一些实施方式中,上述将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集的步骤,包括:
遍历叠前炮集地震数据,获得具有同一同相轴的叠前炮集地震数据;
根据转换算式将叠前炮集地震数据转化为共虚震源道集;转换算式为:
V(s,sv,r)=∫U(s,r,t)δ(t-tr(s,sv,r))dt,
式中,s为炮点信息;sv为虚震源位置信息;t为双程旅行时间;r为检波点;tr为由虚震源sv到检波点r的旅行时间;U(s,r,t)为叠前炮集地震数据;V(s,sv,r)为共虚震源道集。
在一些实施方式中,上述虚震源sv到检波点r的旅行时间tr由以下算式计算得到:
式中,x为炮点的水平位置;xs为虚震源的水平位置;xr为检波点的水平位置;z为虚震源深度;v(o)为虚震源与炮点中点o处的均方根速度;v0为地表速度。
在一些实施方式中,上述对共虚震源道集中的数据沿检波点方向进行滤波操作,滤除共虚震源道集中包含的反射波数据后得到滤波结果的步骤,包括:
遍历共虚震源道集中的数据,获取共虚震源道集中的炮点数据以及虚震源点数据;
对共虚震源道集依次沿着每个炮点以及虚震源点中的检波点方向,进行滤波操作,得到滤波计算结果;
将共虚震源道集中的数据与对应的滤波计算结果进行差值运算,得到滤波结果。
在一些实施方式中,上述滤波操作为中值滤波、双边滤波或高斯滤波中的任意一种。
在一些实施方式中,上述对共虚震源道集依次沿着每个炮点以及虚震源点中的检波点方向,进行滤波操作,得到滤波计算结果的步骤,通过以下算式得以实现:
在一些实施方式中,上述将滤波结果转换至炮集中,得到已完成分离的绕射波叠前炮集的步骤,通过以下算式得以实现:
Ud(s,r,t)=∫Vd(s,sv,r)δ(t-tr(s,sv,r))d sv
式中,Ud(s,r,t)为已完成分离的绕射波叠前炮集;Vd(s,sv,r)为仅含绕射波的共虚震源道集。
第二方面,本发明实施例提供了一种绕射波叠前炮集的分离装置,该装置包括:
数据获取模块,用于获取叠前炮集地震数据;叠前炮集地震数据中包含反射波数据以及绕射波数据;
共虚震源道集转换模块,将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集;
滤波操作模块,用于对共虚震源道集中的数据沿检波点方向进行滤波操作,滤除共虚震源道集中包含的反射波数据后得到滤波结果;
获取模块,用于将滤波结果转换至炮集中,得到已完成分离的绕射波叠前炮集。
第三方面,本发明实施例还提供一种电子设备,包括存储器、处理器,存储器中存储有可在处理器上运行的计算机程序,其中,处理器执行计算机程序时实现上述第一方面所述的方法的步骤。
第四方面,本发明实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质,其中,程序代码使处理器执行上述第一方面所述方法。
本发明实施例带来了以下有益效果:
本发明提供了一种绕射波叠前炮集的分离方法、装置和电子设备,该方法首先获取叠前炮集地震数据,在叠前炮集地震数据中包含反射波数据以及绕射波数据;然后将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集;再对共虚震源道集中的数据沿检波点方向进行滤波操作,滤除共虚震源道集中包含的反射波数据后得到滤波结果;最后将滤波结果转换至炮集中,得到已完成分离的绕射波叠前炮集。该方法利用反射波和绕射波的运动学差异构建共虚震源道集,通过滤波的方式去除共虚震源道集中的反射波,实现了对绕射波叠前炮集的有效分离,提升了叠前炮集中绕射波的分离效果。
本发明的其他特征和优点将在随后的说明书中阐述,或者,部分特征和优点可以从说明书推知或毫无疑义地确定,或者通过实施本发明的上述技术即可得知。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施方式,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的绕射波叠前炮集的分离方法的流程图;
图2为本发明实施例提供的绕射波叠前炮集的分离方法中步骤S102的流程图;
图3为本发明实施例提供的绕射波叠前炮集的分离方法中步骤S103的流程图;
图4为本发明实施例提供的使用绕射波叠前炮集的分离方法所得到的炮集图像;
图5为本发明实施例提供的分离前的原始炮集图;
图6为本发明实施例提供的使用绕射波叠前炮集的分离方法得到的断层成像图;
图7为本发明实施例提供的常规绕射波成像方法所得到的断层成像图;
图8为本发明实施例提供的一种绕射波叠前炮集的分离装置的结构示意图;
图9为本发明实施例提供的一种电子设备的结构示意图。
图标:
810-数据获取模块;820-共虚震源道集转换模块;830-滤波操作模块;840-获取模块;101-处理器;102-存储器;103-总线;104-通信接口。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在地质勘探中,对于地质不连续性的准确预测对于探测精度以及生产安全具有重要意义,例如在散射体、小尺度断层等地下地质不连续性的探测过程中,碳酸盐岩储层中含有大量的裂缝、孔洞等细小的散射体,会有较大的安全生产隐患;煤炭勘探中,断层和陷落柱经常诱发水或气通道。实际生产过程中需要对这些地质不连续性进行准确预测,最大程度的避免开采过程中潜在的突水和瓦斯突出等事故。在油气勘探和生产中具有重要作用。
在进行人工震源操作过程中,震源也称为炮点。施工过程中,采用单炮激发,多检波器接收的方式。形成的地震波在地面传播时,当遇到地质异常区域会形成绕射波。在绕射波中包含该地质异常体的有效数据,可通过分析得到该异常体的详细属性,对于开采相关行业具有指导意义。
绕射波具有揭示小于地震波波长的小尺度物体的特性,然而由于绕射相对于镜面反射较弱,存在于绕射中的细节信息通常被镜面反射所掩盖。此外,在常规处理过程中,特别是在常规动校正和叠加过程中,绕射会遭到破坏。因此,现有技术中对于绕射波的分离方式,还有进一步提升空间。基于此,本发明实施例提供的一种绕射波叠前炮集的分离方法、装置和电子设备,可以利用反射波和绕射波的运动学差异构建共虚震源道集,通过滤波的方式去除共虚震源道集中的反射波,实现了对绕射波叠前炮集的有效分离,提升了叠前炮集中绕射波的分离效果。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种绕射波叠前炮集的分离方法进行详细介绍。
参见图1所示的一种绕射波叠前炮集的分离方法的流程图,其中,该方法具体步骤包括:
步骤S101,获取叠前炮集地震数据;叠前炮集地震数据中包含反射波数据以及绕射波数据。
地震数据的采集,采用人工震源的方式在野外进行激发进而收集地震数据。上述过程中的震源也称为炮点,采用单炮激发、多检波器接收的方式获取数据,这些数据叫做炮集地震数据。叠前炮集地震数据是指地震波还未叠加时的初始数据,包括反射波数据以及绕射波数据,同时还包含了炮点信息以及检波点信息。
步骤S102,将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集。
获取共虚震源道集的目的,是为了分离弱绕射和镜面反射,由于绕射相对于镜面反射较弱,存在于绕射中的细节信息通常被镜面反射所掩盖,因此为获得高分辨率的信息,有必要单独分离或成像绕射波。
对于同一反射器所产生的反射可认为是来自同一虚震源,根据这一特点,反射可以转换到相应的虚震源位置,建立共虚震源道集。与反射不同的是,绕射不能被认为来自虚震源。因此,绕射具有弯曲的形状,而反射在共虚震源道集中表现出平坦特征。转换后的反射波具有同一深度点(即虚震源深度)。该转换过程称为正变换,即将炮集转换至共虚震源道集。
共虚震源道集和炮集有很好的映射,可以互相转换,转换过程几乎没有能量损失,并且在虚震源道集中去除水平反射波比在炮集中去除双曲反射波更容易,因此将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集。
步骤S103,对共虚震源道集中的数据沿检波点方向进行滤波操作,滤除共虚震源道集中包含的反射波数据后得到滤波结果。
将反射波转换到共虚震源位置生成共虚震源道集后,该反射会表现为水平特征;因此,可沿着水平方向对共虚震源道集进行滤波操作,将沿着水平方向振幅连续性较好的反射滤除,将绕射波保留。
在具体实现过程中,可对共虚震源道集的数据依次沿着每个炮点和虚震源点中的检波点方向进行滤波操作,滤波结束后得到的结果与共虚震源道集原始数据进行差值运算,得到仅含绕射波的共虚震源道集。
步骤S104,将滤波结果转换至炮集中,得到已完成分离的绕射波叠前炮集。
对应于上述将共虚震源道集转换至炮集转换过程的正变换,该步骤将炮集转换至共虚震源道集,可称为反变换。由于共虚震源道集和炮集可以互相转换,转换过程几乎没有能量损失,因此将滤波结果反转换至炮集中,即可得到已完成分离的绕射波叠前炮集。
通过上述实施例可知,本实施例中的绕射波叠前炮集的分离方法,首先获取叠前炮集地震数据,在叠前炮集地震数据中包含反射波数据以及绕射波数据;然后将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集;再对共虚震源道集中的数据沿检波点方向进行滤波操作,滤除共虚震源道集中包含的反射波数据后得到滤波结果;最后将滤波结果转换至炮集中,得到已完成分离的绕射波叠前炮集。该方法利用反射波和绕射波的运动学差异构建共虚震源道集,通过滤波的方式去除共虚震源道集中的反射波,实现了对绕射波叠前炮集的有效分离,提升了叠前炮集中绕射波的分离效果,分离后的绕射波波场可用于地下小尺度构造检测,如断层断点,煤层尖灭点,陷落柱等,用于实现这些地下地质异常体的高分辨率成像。
在一些实施方式中,上述将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集的步骤S102,如图2所示,包括:
步骤S201,遍历叠前炮集地震数据,获得具有同一同相轴的叠前炮集地震数据。
由于共虚震源道集的建立是建立在同一反射器所产生反射的基础上,转换后的反射波具有同一深度点(即虚震源深度),因此在具体实施过程中,需要获得具有同一同相轴的叠前炮集地震数据,然后再执行后续的转换过程。
步骤S202,根据转换算式将叠前炮集地震数据转化为共虚震源道集。
转换算式为:
V(s,sv,r)=∫U(s,r,t)δ(t-tr(s,sv,r))dt
上式中,s为炮点信息;sv为虚震源位置信息;t为双程旅行时间;r为检波点;tr为由虚震源sv到检波点r的旅行时间;U(s,r,t)为叠前炮集地震数据;V(s,sv,r)为共虚震源道集。
上述虚震源sv到检波点r的旅行时间tr由以下算式计算得到:
上式中,x为炮点的水平位置;xs为虚震源的水平位置;xr为检波点的水平位置;z为虚震源深度;v(o)为虚震源与炮点中点o处的均方根速度;v0为地表速度。
对共虚震源道集中的数据沿检波点方向进行滤波操作,滤除共虚震源道集中包含的反射波数据后得到滤波结果的步骤S103,如图3所示,包括:
步骤S301,遍历共虚震源道集中的数据,获取共虚震源道集中的炮点数据以及虚震源点数据。
该步骤认为是初始化步骤,通过对共虚震源道集中的数据进行遍历,获取共虚震源道集中的炮点数据以及虚震源点数据,具体实现过程中可一定范围的数据集进行遍历。
步骤S302,对共虚震源道集依次沿着每个炮点以及虚震源点中的检波点方向,进行滤波操作,得到滤波计算结果。
滤波的过程是沿着检波点的方向得以实现,通常是沿着炮点以及虚震源点的方向进行滤波操作,在虚震源道集中去除水平反射波比在炮集中去除双曲反射波更容易,因此该步骤可滤除共虚震源道集中的水平反射波,保留具有曲线特征的绕射波信号。
具体的,上述步骤通过以下算式得以实现:
在一些实施方式中,该步骤中的滤波操作为中值滤波、双边滤波或高斯滤波中的任意一种。
步骤S303,将共虚震源道集中的数据与对应的滤波计算结果进行差值运算,得到滤波结果。
具体实现过程可将共虚震源道集中的数据与对应的滤波计算结果相减,最终得到滤波结果中已将反射波进行去除,只含绕射波的共虚震源道集。
在一些实施方式中,上述将滤波结果转换至炮集中,得到已完成分离的绕射波叠前炮集的步骤S104,通过以下算式得以实现:
Ud(s,r,t)=∫Vd(s,sv,r)δ(t-tr(s,sv,r))d sv
上式中,Ud(s,r,t)为已完成分离的绕射波叠前炮集;Vd(s,sv,r)为仅含绕射波的共虚震源道集。
对应于上述将共虚震源道集转换至炮集转换过程的正变换,上述算式为可称为反变换,可将炮集转换至共虚震源道集。由于共虚震源道集和炮集可以互相转换,转换过程几乎没有能量损失,因此将滤波结果反转换至炮集中,即可得到已完成分离的绕射波叠前炮集。
上述实施例中的所有时间单位为秒,速度单位为米/秒,距离或位置单位为米,均为国际单位制。
通过上述实施例提到的绕射波叠前炮集的分离方法,得到绕射波分离后的炮集图像,如图4所示;与分离前的原始炮集图像图5进行对比,从箭头指引的反射波成像可以发现,通过使用本实施例中的绕射波叠前炮集的分离方法后,可对箭头所指的反射波进行了有效抑制。
通过上述实施例提到的绕射波叠前炮集的分离方法最终得到的断层成像结果如图6所示,与常规的断层成像图7进行对比,可发现,图6中箭头所指的区域就是断层边界区域,而图7中相同位置的成像效果远不如图6中箭头所指区域的成像效果。
上述实施例可知,通过利用反射波和绕射波的运动学差异构建共虚震源道集,通过滤波的方式去除共虚震源道集中的反射波,实现了对绕射波叠前炮集的有效分离,提升了叠前炮集中绕射波的分离效果,分离后的绕射波波场可用于地下小尺度构造检测,如断层断点,煤层尖灭点,陷落柱等,用于实现这些地下地质异常体的高分辨率成像。
对应于上述方法实施例,本发明实施例还提供了一种绕射波叠前炮集的分离装置,其结构示意图如图8所示,其中,该装置包括:
数据获取模块810,用于获取叠前炮集地震数据;叠前炮集地震数据中包含反射波数据以及绕射波数据;
共虚震源道集转换模块820,将同相轴相同的叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集;
滤波操作模块830,用于对共虚震源道集中的数据沿检波点方向进行滤波操作,滤除共虚震源道集中包含的反射波数据后得到滤波结果;
获取模块840,用于将滤波结果转换至炮集中,得到已完成分离的绕射波叠前炮集。
本发明实施例提供的绕射波叠前炮集的分离装置,与上述实施例提供的绕射波叠前炮集的分离方法具有相同的技术特征,所以也能解决相同的技术问题,达到相同的技术效果。为简要描述,实施例部分未提及之处,可参考前述方法实施例中相应内容。
本实施例还提供一种电子设备,为该电子设备的结构示意图如图9所示,该设备包括处理器101和存储器102;其中,存储器102用于存储一条或多条计算机指令,一条或多条计算机指令被处理器执行,以实现上述绕射波叠前炮集的分离方法。
图9所示的电子设备还包括总线103和通信接口104,处理器101、通信接口104和存储器102通过总线103连接。
其中,存储器102可能包含高速随机存取存储器(RAM,Random Access Memory),也可能还包括非不稳定的存储器(non-volatile memory),例如至少一个磁盘存储器。总线103可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示,图9中仅用一个双向箭头表示,但并不表示仅有一根总线或一种类型的总线。
通信接口104用于通过网络接口与至少一个用户终端及其它网络单元连接,将封装好的IPv4报文或IPv4报文通过网络接口发送至用户终端。
处理器101可能是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器101中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器101可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(DigitalSignal Processor,简称DSP)、专用集成电路(Application Specific IntegratedCircuit,简称ASIC)、现场可编程门阵列(Field-Programmable Gate Array,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本公开实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本公开实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器102,处理器101读取存储器102中的信息,结合其硬件完成前述实施例的方法的步骤。
本发明实施例还提供了一种计算机可读存储介质,该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器运行时执行前述实施例的方法的步骤。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统、设备和方法,可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,又例如,多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口,设备或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以用软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种绕射波叠前炮集的分离方法,其特征在于,所述方法包括:
获取叠前炮集地震数据;所述叠前炮集地震数据中包含反射波数据以及绕射波数据;
将同相轴相同的所述叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集;
对所述共虚震源道集中的数据沿检波点方向进行滤波操作,滤除所述共虚震源道集中包含的反射波数据后得到滤波结果;
将所述滤波结果转换至炮集中,得到已完成分离的所述绕射波叠前炮集。
2.根据权利要求1所述的方法,其特征在于,将同相轴相同的所述叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集的步骤,包括:
遍历所述叠前炮集地震数据,获得具有同一同相轴的所述叠前炮集地震数据;
根据转换算式将所述叠前炮集地震数据转化为所述共虚震源道集;所述转换算式为:
V(s,sv,r)=∫U(s,r,t)δ(t-tr(s,sv,r))dt
式中,s为炮点信息;sv为虚震源位置信息;t为双程旅行时间;r为检波点;tr为由虚震源sv到检波点r的旅行时间;U(s,r,t)为叠前炮集地震数据;V(s,sv,r)为共虚震源道集。
4.根据权利要求2所述的方法,其特征在于,对所述共虚震源道集中的数据沿检波点方向进行滤波操作,滤除所述共虚震源道集中包含的反射波数据后得到滤波结果的步骤,包括:
遍历所述共虚震源道集中的数据,获取所述共虚震源道集中的炮点数据以及虚震源点数据;
对所述共虚震源道集依次沿着每个炮点以及虚震源点中的检波点方向,进行滤波操作,得到滤波计算结果;
将所述共虚震源道集中的数据与对应的所述滤波计算结果进行差值运算,得到所述滤波结果。
5.根据权利要求4所述的方法,其特征在于,所述滤波操作为中值滤波、双边滤波或高斯滤波中的任意一种。
7.根据权利要求2所述的方法,其特征在于,将所述滤波结果转换至炮集中,得到已完成分离的所述绕射波叠前炮集的步骤,通过以下算式得以实现:
Ud(s,r,t)=∫Vd(s,sv,r)δ(t-tr(s,sv,r))dsv
式中,Ud(s,r,t)为已完成分离的所述绕射波叠前炮集;Vd(s,sv,r)为仅含绕射波的共虚震源道集。
8.一种绕射波叠前炮集的分离装置,其特征在于,所述装置包括:
数据获取模块,用于获取叠前炮集地震数据;所述叠前炮集地震数据中包含反射波数据以及绕射波数据;
共虚震源道集转换模块,将同相轴相同的所述叠前炮集地震数据的反射波数据转换至相应的虚震源处,得到共虚震源道集;
滤波操作模块,用于对所述共虚震源道集中的数据沿检波点方向进行滤波操作,滤除所述共虚震源道集中包含的反射波数据后得到滤波结果;
获取模块,用于将所述滤波结果转换至炮集中,得到已完成分离的所述绕射波叠前炮集。
9.一种电子设备,包括存储器、处理器,所述存储器中存储有可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述权利要求1至7任一项所述的方法的步骤。
10.一种具有处理器可执行的非易失的程序代码的计算机可读介质,其特征在于,所述程序代码使所述处理器执行所述权利要求1至7任一所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010309275.9A CN111337982A (zh) | 2020-04-17 | 2020-04-17 | 绕射波叠前炮集的分离方法、装置和电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010309275.9A CN111337982A (zh) | 2020-04-17 | 2020-04-17 | 绕射波叠前炮集的分离方法、装置和电子设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111337982A true CN111337982A (zh) | 2020-06-26 |
Family
ID=71182841
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010309275.9A Pending CN111337982A (zh) | 2020-04-17 | 2020-04-17 | 绕射波叠前炮集的分离方法、装置和电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111337982A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113640872A (zh) * | 2021-08-12 | 2021-11-12 | 中国矿业大学(北京) | 绕射波分离方法、装置和电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102455439A (zh) * | 2010-11-02 | 2012-05-16 | 中国石油大学(北京) | 基于克希霍夫积分法的绕射波场分离方法 |
CN108845358A (zh) * | 2018-09-26 | 2018-11-20 | 中国矿业大学(北京) | 断层及构造异常体识别方法及装置 |
CN108919350A (zh) * | 2018-09-26 | 2018-11-30 | 中国矿业大学(北京) | 绕射波分离方法及装置 |
US20190187313A1 (en) * | 2017-12-15 | 2019-06-20 | Pgs Geophysical As | Processes and systems to enhance illumination and resolution of seismic images using multiple reflected wavefields |
-
2020
- 2020-04-17 CN CN202010309275.9A patent/CN111337982A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102455439A (zh) * | 2010-11-02 | 2012-05-16 | 中国石油大学(北京) | 基于克希霍夫积分法的绕射波场分离方法 |
US20190187313A1 (en) * | 2017-12-15 | 2019-06-20 | Pgs Geophysical As | Processes and systems to enhance illumination and resolution of seismic images using multiple reflected wavefields |
CN108845358A (zh) * | 2018-09-26 | 2018-11-20 | 中国矿业大学(北京) | 断层及构造异常体识别方法及装置 |
CN108919350A (zh) * | 2018-09-26 | 2018-11-30 | 中国矿业大学(北京) | 绕射波分离方法及装置 |
Non-Patent Citations (8)
Title |
---|
CHUANGJIAN LI ET AL.: "Diffraction imaging using an adaptive phase filter", 《GEOPHYSICAL PROSPECTING》 * |
CHUANGJIAN LI ET AL.: "Enhancing Subsurface Diffractions Using Demigration Method", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
JINGTAO ZHAO ET AL.: "Separating prestack diffractions with SVMF in the flattened shot domain", 《JOURNAL OF GEOPHYSICS AND ENGINEERING》 * |
T.J. MOSER ET AL.: "Diffraction imaging in depth", 《GEOPHYSICAL PROSPECTING》 * |
V. KHAIDUKOV ET AL.: "Diffraction imaging by focusing-defocusing:An outlook on seismic superresolution", 《GEOPHYSICS》 * |
李清娟等: "《人工智能与产业变革》", 31 January 2020, 上海:上海财经大学出版社 * |
王志辉等: "地震绕射波分离方法研究进展", 《地球物理学进展》 * |
黄建平等: "《绕射波分离与成像方法》", 31 July 2015, 东营:中国石油大学出版社 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113640872A (zh) * | 2021-08-12 | 2021-11-12 | 中国矿业大学(北京) | 绕射波分离方法、装置和电子设备 |
CN113640872B (zh) * | 2021-08-12 | 2022-03-08 | 中国矿业大学(北京) | 绕射波分离方法、装置和电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111337981B (zh) | 绕射波成像方法、装置和电子设备 | |
CN104536044B (zh) | 一种地震数据的插值去噪方法以及系统 | |
CN106855640B (zh) | 一种海底电缆地震数据初至计算方法及装置 | |
JP5874039B2 (ja) | ステレオ画像処理装置及びステレオ画像処理方法 | |
CN111337982A (zh) | 绕射波叠前炮集的分离方法、装置和电子设备 | |
WO2005017564A1 (en) | Dip value in seismic images | |
CN110118993B (zh) | 绕射波成像方法及装置 | |
CN103412325B (zh) | 一种地震纵、横波波场分离与去噪的方法 | |
CN105301648A (zh) | 一种获取共反射面元叠加参数的方法 | |
CN109143343A (zh) | 绕射波成像方法、装置和电子设备 | |
CN109116423B (zh) | 一种绕射多次波压制方法及装置 | |
CN103941289A (zh) | 一种绕射点成像的方法和设备 | |
CN104765069A (zh) | 一种压制同步激发采集邻炮干扰的方法 | |
CN111399038B (zh) | 斜率参数提取方法、装置及计算机可读存储介质 | |
Wang et al. | Ground roll wave suppression based on wavelet frequency division and radial trace transform | |
CN109212590B (zh) | 一种保幅的逆时偏移低频噪音压制方法及系统 | |
CN112379420B (zh) | 高精度弯曲测线叠前时间域成像方法及装置 | |
CN104954704A (zh) | 一种用于拉曼光谱仪ccd信号降噪方法 | |
CN115184987B (zh) | 地层倾角信息的提取方法、装置、服务器及存储介质 | |
CN116184413B (zh) | 用于多波束测深系统的底检测方法及装置 | |
CN115903043A (zh) | 一种绕射波分离方法和装置 | |
CN116027412A (zh) | 一种基于包络误差函数的反射波波形反演方法及电子设备 | |
CN111142155B (zh) | 基于矢量面元的复数域弱信号恢复方法及电子设备 | |
CN116381786B (zh) | 绕射波成像方法和装置 | |
CN113391348B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200626 |