CN109581477A - 预测地震反射界面的方法及*** - Google Patents
预测地震反射界面的方法及*** Download PDFInfo
- Publication number
- CN109581477A CN109581477A CN201710910942.7A CN201710910942A CN109581477A CN 109581477 A CN109581477 A CN 109581477A CN 201710910942 A CN201710910942 A CN 201710910942A CN 109581477 A CN109581477 A CN 109581477A
- Authority
- CN
- China
- Prior art keywords
- seismic
- function
- wavelet
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000006870 function Effects 0.000 claims abstract description 44
- 230000009131 signaling function Effects 0.000 claims abstract description 35
- 238000005070 sampling Methods 0.000 claims abstract description 24
- 238000001228 spectrum Methods 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 5
- 230000002123 temporal effect Effects 0.000 abstract description 6
- 238000012545 processing Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012360 testing method 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/01—Measuring or predicting earthquakes
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
本发明公开了一种预测地震反射界面的方法及***,包括:1)取地震数据的其中一道数据作为输入数据,建立地震信号函数;2)从步骤1)选取的道数据中提取子波,建立子波函数;3)基于地震信号函数和子波函数,获得地震信号的时频分析函数,在设定频段内检测每个时间样点处是否为地震反射点;4)取下一道地震数据进行检测,直到对所有地震道检测完成,获得地震反射界面。本发明通过信号时频分析理论,推导并预测出可能出现地震反射的位置,实现了从较低时间分辨率的地震剖面上获得反射界面。
Description
技术领域
本发明属于石油天然气地球领域,更具体地,涉及一种预测地震反射界面的方法及***。
背景技术
地震勘探的目的是通过地下地层的反射波信息查明地下地质构造进而为矿产普查、油气勘探服务,由于大地滤波的作用初始的地震脉冲信号会随着地层深度逐渐衰减,此种作用导致的后果之一就是使得较薄地层的界面在地震剖面上难以识别也即地震时间分辨率下降,要想从较低时间分辨率的地震剖面上得到地层反射界面不是容易的事情,需要对地震数据进行特殊的处理,这也是地球物理中的一个热门研究方向前人也做过很多的研究和探索,例如《提高分辨率技术在苏格里气田西区的应用》通过对实际资料进行反褶积和Q补偿方法试验,认为在苏里格气田西区应用地表一致性反褶积和Q补偿串联的方法能有效提高地震资料的分辨率,为后续薄砂层的储层预测提供了高品质的地震成果;同时通过参数优选总结出了苏里格气田西区品质因子Q的范围,该研究成果为同类地区开展高分辨率方法研究提供了借鉴。《井间地震提高分辨率技术》(孙渊,长安大学硕士论文)井间地震信号的主频通常是地面地震的数倍甚至更高,可以达到几百赫兹。通过反褶积处理、反Q滤波、反射系数反演等方法是提高分辨率的常用方法,通过这些特殊的地震资料处理以达到提高对薄层的识别能力。因此,我们提供一种方法,从相对低分辨率的地震数据中恢复被低分辨率子波掩盖的薄层。
发明内容
本发明通过信号时频分析理论,推导并预测出可能出现地震反射的位置,实现了从较低时间分辨率的地震剖面上获得反射界面。
根据本发明的一方面,提出了一种预测地震反射界面的方法,该方法可以包括:1)取地震数据的其中一道数据作为输入数据,建立地震信号函数;2)从步骤1)选取的道数据中提取子波,建立子波函数;3)基于所述地震信号函数和所述子波函数,获得地震信号的时频分析函数,在设定频段内检测每个时间样点处是否为地震反射点;4)取下一道地震数据进行检测,直到对所有地震道检测完成,获得地震反射界面。
优选地,所述地震信号函数s(t)表示为:
s(t)=r(y0,x0,t), (1)
其中,y0为线位置;x0为道位置;t为时间。
优选地,所述子波函数表示为:
其中,σ为能量衰减因子;t为时间。
优选地,所述地震信号的时频分析函数表示为:
s(t,f)=CWT(s(t),W(t)), (3)
其中,s(t)为地震信号函数;W(t)为子波函数;CWT为连续小波变换;f为频率。
优选地,所述时间样点t0满足以下公式时,将其判断为地震反射点,
其中,t0为子波延迟时间;s(t0,f)为t0时刻的地震信号频谱;f为频率。
优选地,如果判断所述时间样点t0处为地震反射点,则令r(y0,x0,t0)=1,否则r(y0,x0,t0)=0;对所有地震道检测完成后,获得的r(y,x,t)即为所述地震反射界面。
根据本发明的另一方面,提出了一种预测地震反射界面的***,所述***包括:存储器,存储有计算机可执行指令;处理器,运行所述存储器上的计算可执行指令时,所述处理器实现以下步骤:1)取地震数据的其中一道数据作为输入数据,建立地震信号函数;2)从步骤1)选取的道数据中提取子波,建立子波函数;3)基于所述地震信号函数和所述子波函数,获得地震信号的时频分析函数,在设定频段内检测每个时间样点处是否为地震反射点;4)取下一道地震数据进行检测,直到对所有地震道检测完成,获得地震反射界面。
优选地,所述地震信号函数s(t)表示为:
s(t)=r(y0,x0,t), (1)
其中,y0为线位置;x0为道位置;t为时间。
优选地,所述地震信号的时频分析函数表示为:
s(t,f)=CWT(s(t),W(t)), (3)
其中,s(t)为地震信号函数;W(t)为子波函数;CWT为连续小波变换;f为频率。
优选地,所述时间样点t0满足以下公式时,将其判断为地震反射点,
其中,t0为子波延迟时间;s(t0,f)为t0时刻的地震信号频谱;f为频率。
本发明的有益效果在于:通过建立选定地震道的地震信号函数和与其对应的子波函数,经过连续小波变换获得地震信号的时频分析函数,在设定频段内对每个时间样点进行检测是否为地震反射点并记录,然后取下一道地震数据进行检测,直到对所有地震道检测完成,从相对低分辨率地震数据中恢复被低分辨率子波掩盖的薄层。
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的预测地震反射界面的方法的流程图。
图2示出了根据本发明的一个实施例的原始地震数据。
图3示出了根据本发明的一个实施例的地层界面的解释结果。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
实施例1
在该实施例中,根据本发明的预测地震反射界面的方法可以包括:1)取地震数据的其中一道数据作为输入数据,建立地震信号函数;2)从步骤1)选取的道数据中提取子波,建立子波函数;3)基于地震信号函数和子波函数,获得地震信号的时频分析函数,在设定频段内检测每个时间样点处是否为地震反射点;4)取下一道地震数据进行检测,直到对所有地震道检测完成,获得地震反射界面。
该实施例通过信号时频分析理论,推导并预测出可能出现地震反射的位置,实现了从较低时间分辨率的地震剖面上获得反射界面。
图1示出了根据本发明的预测地震反射界面的方法的流程图。下面参考图1详细说明根据本发明的预测地震反射界面的方法的具体步骤。
步骤1,取地震数据的其中一道数据作为输入数据,建立地震信号函数。
在一个示例中,地震信号函数s(t)表示为:
s(t)=r(y0,x0,t), (1)
其中,y0为线位置;x0为道位置;t为时间。
步骤2,从步骤1)选取的道数据中提取子波,建立子波函数。
在一个示例中,子波函数表示为:
其中,σ为能量衰减因子;t为时间。
对子波做傅里叶变换,可得:
设子波时间延迟t0,则下面对时间延迟t0的子波做如下处理:
u(f)=FFT(W(t-t0)), (6)
在这里如果|σ|<<2πf,则
此时,
步骤3,基于地震信号函数和子波函数,获得地震信号的时频分析函数,在设定频段内检测每个时间样点处是否为地震反射点。
地震信号是典型的非平稳信号即信号的频率随着时间发生变化,常规的傅立叶变换难以满足地震信号的时频分析需要。
在一个示例中,地震信号的时频分析函数表示为:
s(t,f)=CWT(s(t),W(t)), (3)
其中,s(t)为地震信号函数;W(t)为子波函数;CWT为连续小波变换;f为频率。
设在t0时刻地震信号的频谱为:
u(f)=s(t0,f), (11)
在一个示例中,时间样点t0满足以下公式时,将其判断为地震反射点,
如果在t0时刻地震信号的频谱与分解的子波频谱高度一致,则频谱会满足式(10)的表达式即:
其中,t0为子波延迟时间;s(t0,f)为t0时刻的地震信号频谱;f为频率。
具体地,方法的使用前提条件是:2πf>>σ,在本实施例中我们取f的范围[60—120]单位是hz,2πf的范围为[376.8—753.6],单位是rad/s,σ取0.5。依次对每个地震信号的时间样点按照式检测,如果满足式(4)则认为此点是地震反射点,否则认为不是地震反射点。
步骤4,取下一道地震数据进行检测,直到对所有地震道检测完成,获得地震反射界面。
在一个示例中,如果判断时间样点t0处为地震反射点,则令r(y0,x0,t0)=1,否则r(y0,x0,t0)=0;对所有地震道检测完成后,获得的r(y,x,t)即为地震反射界面。
本实施例通过建立选定地震道的地震信号函数和与其对应的子波函数,经过连续小波变换获得地震信号的时频分析函数,在设定频段内对每个时间样点进行检测是否为地震反射点并记录,然后取下一道地震数据进行检测,直到对所有地震道检测完成,从相对低分辨率地震数据中恢复被低分辨率子波掩盖的薄层。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
首先,取三维地震数据,取地震的其中一道数据作为输入数据,建立如公式(1)所示的地震信号函数s(t)=r(y0,x0,t);
其次,从选取的道数据中提取子波,建立如公式(2)所示的子波函数W(t),在本应用实例中σ取0.5;
然后,对地震信号函数s(t)按照公式(3)做连续小波变换,在[60—120]频段内对每一个时间样点检测是否满足公式(3),如果满足,记三维数组r(y0,x0,t0)=1.0,否则,r(y0,x0,t0)=0.0,将此道地震数据的所有时间样点都检测完毕;
最后,取下一道地震数据进行检测,直到对所有地震道检测完成为止,最终得到的三维数组r(y,x,t)就是最后的结果。
图2示出了在上述应用示例中的原始地震数据,从图2可以看出,地震信号子波有明显的衰减,由于复合效应部分薄层难以识别;图3示出了在上述应用示例中的地层界面的解释结果,是应用本发明提出的方法处理后的结果,实现了从较低时间分辨率的地震剖面上获得反射界面。
本应用示例通过建立选定地震道的地震信号函数和与其对应的子波函数,经过连续小波变换获得地震信号的时频分析函数,在设定频段内对每个时间样点进行检测是否为地震反射点并记录,然后取下一道地震数据进行检测,直到对所有地震道检测完成,从相对低分辨率地震数据中恢复被低分辨率子波掩盖的薄层。
实施例2
根据本发明的实施例,提供了一种预测地震反射界面的***,***包括:存储器,存储有计算机可执行指令;处理器,运行存储器上的计算可执行指令时,处理器实现以下步骤:1)取地震数据的其中一道数据作为输入数据,建立地震信号函数;2)从步骤1)选取的道数据中提取子波,建立子波函数;3)基于地震信号函数和子波函数,获得地震信号的时频分析函数,在设定频段内检测每个时间样点处是否为地震反射点;4)取下一道地震数据进行检测,直到对所有地震道检测完成,获得地震反射界面。
该实施例通过信号时频分析理论,推导并预测出可能出现地震反射的位置,实现了从较低时间分辨率的地震剖面上获得反射界面。
在一个示例中,地震信号函数s(t)表示为:
s(t)=r(y0,x0,t), (1)
其中,y0为线位置;x0为道位置;t为时间。
在一个示例中,地震信号的时频分析函数表示为:
s(t,f)=CWT(s(t),W(t)), (3)
其中,s(t)为地震信号函数;W(t)为子波函数;CWT为连续小波变换;f为频率。
在一个示例中,时间样点t0满足以下公式时,将其判断为地震反射点,
其中,t0为子波延迟时间;s(t0,f)为t0时刻的地震信号频谱;f为频率。
本实施例通过建立选定地震道的地震信号函数和与其对应的子波函数,经过连续小波变换获得地震信号的时频分析函数,在设定频段内对每个时间样点进行检测是否为地震反射点并记录,然后取下一道地震数据进行检测,直到对所有地震道检测完成,从相对低分辨率地震数据中恢复被低分辨率子波掩盖的薄层。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种预测地震反射界面的方法,其特征在于,该方法包括:
1)取地震数据的其中一道数据作为输入数据,建立地震信号函数;
2)从步骤1)选取的道数据中提取子波,建立子波函数;
3)基于所述地震信号函数和所述子波函数,获得地震信号的时频分析函数,在设定频段内检测每个时间样点处是否为地震反射点;
4)取下一道地震数据进行检测,直到对所有地震道检测完成,获得地震反射界面。
2.根据权利要求1所述的预测地震反射界面的方法,其中,所述地震信号函数s(t)表示为:
s(t)=r(y0,x0,t), (1)
其中,y0为线位置;x0为道位置;t为时间。
3.根据权利要求1所述的预测地震反射界面的方法,其中,所述子波函数表示为:
其中,σ为能量衰减因子;t为时间。
4.根据权利要求1所述的预测地震反射界面的方法,其中,所述地震信号的时频分析函数表示为:
s(t,f)=CWT(s(t),W(t)), (3)
其中,s(t)为地震信号函数;W(t)为子波函数;CWT为连续小波变换;f为频率。
5.根据权利要求1所述的预测地震反射界面的方法,其中,所述时间样点t0满足以下公式时,将其判断为地震反射点,
其中,t0为子波延迟时间;s(t0,f)为t0时刻的地震信号频谱;f为频率。
6.根据权利要求5所述的预测地震反射界面的方法,其中,如果判断所述时间样点t0处为地震反射点,则令r(y0,x0,t0)=1,否则r(y0,x0,t0)=0;对所有地震道检测完成后,获得的r(y,x,t)即为所述地震反射界面。
7.一种预测地震反射界面的***,其特征在于,所述***包括:
存储器,存储有计算机可执行指令;
处理器,运行所述存储器上的计算可执行指令时,所述处理器实现以下步骤:
1)取地震数据的其中一道数据作为输入数据,建立地震信号函数;
2)从步骤1)选取的道数据中提取子波,建立子波函数;
3)基于所述地震信号函数和所述子波函数,获得地震信号的时频分析函数,在设定频段内检测每个时间样点处是否为地震反射点;
4)取下一道地震数据进行检测,直到对所有地震道检测完成,获得地震反射界面。
8.根据权利要求7所述的预测地震反射界面的***,其中,所述地震信号函数s(t)表示为:
s(t)=r(y0,x0,t), (1)
其中,y0为线位置;x0为道位置;t为时间。
9.根据权利要求7所述的预测地震反射界面的***,其中,所述地震信号的时频分析函数表示为:
s(t,f)=CWT(s(t),W(t)), (3)
其中,s(t)为地震信号函数;W(t)为子波函数;CWT为连续小波变换;f为频率。
10.根据权利要求7所述的预测地震反射界面的***,其中,所述时间样点t0满足以下公式时,将其判断为地震反射点,
其中,t0为子波延迟时间;s(t0,f)为t0时刻的地震信号频谱;f为频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710910942.7A CN109581477B (zh) | 2017-09-29 | 2017-09-29 | 预测地震反射界面的方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710910942.7A CN109581477B (zh) | 2017-09-29 | 2017-09-29 | 预测地震反射界面的方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109581477A true CN109581477A (zh) | 2019-04-05 |
CN109581477B CN109581477B (zh) | 2020-08-25 |
Family
ID=65919123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710910942.7A Active CN109581477B (zh) | 2017-09-29 | 2017-09-29 | 预测地震反射界面的方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109581477B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1086318A (zh) * | 1993-10-27 | 1994-05-04 | 地质矿产部西北石油地质局 | 地震油气检测方法 |
CN102305941A (zh) * | 2011-05-25 | 2012-01-04 | 东北石油大学 | 由叠前时间偏移直接扫描确定地层叠加品质因子方法 |
CN103235339A (zh) * | 2013-04-09 | 2013-08-07 | 中国石油大学(北京) | 一种时频分解地震流体识别方法 |
CN103364832A (zh) * | 2013-07-01 | 2013-10-23 | 西安交通大学 | 一种基于自适应最优核时频分布的地震衰减定性估计方法 |
CN103558635A (zh) * | 2013-10-30 | 2014-02-05 | 北京诺克斯达石油科技有限公司 | 基于偶函数地震响应以估算薄层厚度的方法及装置 |
CN104237941A (zh) * | 2013-06-18 | 2014-12-24 | 中国石油化工股份有限公司 | 一种基于频率衰减的煤层气预测方法 |
CN106291689A (zh) * | 2016-08-18 | 2017-01-04 | 中国石油天然气集团公司 | 一种提取地震数据频散属性的处理方法、装置及预测*** |
CN106610506A (zh) * | 2015-10-26 | 2017-05-03 | 中国石油天然气股份有限公司 | 地震勘探薄层识别方法 |
-
2017
- 2017-09-29 CN CN201710910942.7A patent/CN109581477B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1086318A (zh) * | 1993-10-27 | 1994-05-04 | 地质矿产部西北石油地质局 | 地震油气检测方法 |
CN102305941A (zh) * | 2011-05-25 | 2012-01-04 | 东北石油大学 | 由叠前时间偏移直接扫描确定地层叠加品质因子方法 |
CN103235339A (zh) * | 2013-04-09 | 2013-08-07 | 中国石油大学(北京) | 一种时频分解地震流体识别方法 |
CN104237941A (zh) * | 2013-06-18 | 2014-12-24 | 中国石油化工股份有限公司 | 一种基于频率衰减的煤层气预测方法 |
CN103364832A (zh) * | 2013-07-01 | 2013-10-23 | 西安交通大学 | 一种基于自适应最优核时频分布的地震衰减定性估计方法 |
CN103558635A (zh) * | 2013-10-30 | 2014-02-05 | 北京诺克斯达石油科技有限公司 | 基于偶函数地震响应以估算薄层厚度的方法及装置 |
CN106610506A (zh) * | 2015-10-26 | 2017-05-03 | 中国石油天然气股份有限公司 | 地震勘探薄层识别方法 |
CN106291689A (zh) * | 2016-08-18 | 2017-01-04 | 中国石油天然气集团公司 | 一种提取地震数据频散属性的处理方法、装置及预测*** |
Non-Patent Citations (3)
Title |
---|
万涛 等: "三参数小波及其在薄互层地震资料分析中的应用", 《煤田地质与勘探》 * |
高静怀 等: "三参数小波及其在地震资料分析中的应用", 《地球物理学报》 * |
高静怀 等: "广义S变换与薄互层地震响应分析", 《地球物理学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109581477B (zh) | 2020-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
de Ridder et al. | Ambient seismic noise eikonal tomography for near-surface imaging at Valhall | |
Wang | Multichannel matching pursuit for seismic trace decomposition | |
Wu et al. | Directional illumination analysis using beamlet decomposition and propagation | |
CN100383557C (zh) | 一种提高地震分辨率的方法 | |
US11880011B2 (en) | Surface wave prediction and removal from seismic data | |
CN107102356B (zh) | 基于ceemd的地震信号高分辨率处理方法 | |
US20050273266A1 (en) | Seismic event correlation and Vp-Vs estimation | |
Wang et al. | Reverse time migration of multiples: Eliminating migration artifacts in angle domain common image gathers | |
MX2012004856A (es) | Metodos y aparatos para procesar los datos de series de tiempo para propagar señales en formacion subterranea. | |
CN109738951A (zh) | 一种基于地震同相轴子波谱的时变反褶积方法 | |
CN101545981A (zh) | 可控震源地震数据零相位子波最小相位化方法 | |
CN105116449B (zh) | 一种弱反射储层的识别方法 | |
CA2774170A1 (en) | Seismic imaging systems and methods employing a 3d reverse time migration with tilted transverse isotropy | |
Wang et al. | Separation and imaging of seismic diffractions using a localized rank-reduction method with adaptively selected ranks | |
CN113805237B (zh) | 偏移陆地交叉排列地震的使用压缩感测模型的方法和*** | |
CN107390270B (zh) | 一种基于弹性波逆时偏移ADCIGs的AVA分析方法 | |
EP2537049A2 (en) | Estimating internal multiples in seismic data | |
Assimaki et al. | Attenuation and velocity structure for site response analyses via downhole seismogram inversion | |
Zhao et al. | Denoising of seismic data in desert environment based on a variational mode decomposition and a convolutional neural network | |
CN112255690A (zh) | 基于地震相位分解的自适应围岩强反射分离方法 | |
Xie et al. | Observation of higher‐mode Rayleigh waves from ambient noise in the Tarim Basin, China | |
CN109581477A (zh) | 预测地震反射界面的方法及*** | |
CN110673211B (zh) | 一种基于测井与地震数据的品质因子建模方法 | |
Spitzer et al. | Enhancing subbasalt reflections using parabolic τ-p transformation | |
Guo et al. | Denoising and wavefield separation method for DAS VSP via deep learning |
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 |