CN114814931A - 一种观测***排列长度选择方法 - Google Patents
一种观测***排列长度选择方法 Download PDFInfo
- Publication number
- CN114814931A CN114814931A CN202210431768.9A CN202210431768A CN114814931A CN 114814931 A CN114814931 A CN 114814931A CN 202210431768 A CN202210431768 A CN 202210431768A CN 114814931 A CN114814931 A CN 114814931A
- Authority
- CN
- China
- Prior art keywords
- seismic
- point
- response
- formula
- reflection
- 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
- 238000010187 selection method Methods 0.000 title claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims description 18
- 238000012937 correction Methods 0.000 claims description 17
- 230000005284 excitation Effects 0.000 claims description 6
- 239000013535 sea water Substances 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 13
- 230000014509 gene expression Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 238000004445 quantitative analysis Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/168—Deployment of receiver elements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/133—Generating seismic energy using fluidic driving means, e.g. highly pressurised fluids; using implosion
- G01V1/137—Generating seismic energy using fluidic driving means, e.g. highly pressurised fluids; using implosion which fluid escapes from the generator in a pulsating manner, e.g. for generating bursts, airguns
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3817—Positioning of seismic devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3843—Deployment of seismic devices, e.g. of streamers
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
技术领域
本发明涉及野外地质调查和油气资源勘探技术领域,尤其涉及一种观测***排列长度选择方法。
背景技术
野外地质调查及油气资源勘探中观测***参数的选择是获取目标层位高质量成像的关键,特别是合理排列长度的选择对于潜山、凹陷、逆掩推覆等复杂构造区域的成像尤为重要。随着油气资源开发迈向中深层,地震地质条件的复杂程度提高,地震采集精度和分辨率要求随之提高,观测***排列长度的精准确定受到越来越多的关注。过小的排列长度使得有效信息丢失,过大的排列长度受动校拉伸的影响使得优势主频降低并且容易接收到更多的噪音,造成成像质量的降低。传统的观测***排列长度的选择是基于水平叠加理论,此方法的前提是假设地下为水平层状介质,但是该方法对于双复杂区勘探时地震波场会产生严重畸变。而面向地质目标的地震波正演模拟是解决该问题的有效手段,目前正演模拟技术主要基于两个方向,基于射线理论的正演模拟和基于波动理论的正演模拟,两种技术分别在计算效率和计算精度上各有优势,但是在计算精度和计算速度上存在问题。针对以上问题,发展起来了以高斯射线束方法为代表的地震波束理论。既有运动学特征,又有波动力学特征,可以在一定程度上克服射线的盲区效应,提高正演照明精度,同时具有高效性和灵活性,能够适应复杂的地质条件和地震采集***。但是该方法仍然存在问题,比如在参数优化和属性评价等方面不具普适性。
综上所述,以往观测***排列长度的选择方法和技术较复杂,在工程量和普适性方面存在问题。
发明内容
本发明的目的在于提供一种基于波场传播理论的观测***排列长度选择方法,该更简单、更具有普适性。
一种观测***排列长度选择方法,包括以下步骤:
步骤S101:基于波场传播理论获取地表任意一点的地震响应公式;
步骤S102:根据所述地震响应公式,从不同考虑因素确定最佳排列长度的选择准则;
所述不同考虑因素包括:目标层深度、速度分析精度、动校拉伸、反射波能量、AVO精度;
步骤S103:利用所述最佳排列长度的选择准则,分别得到所述不同考虑因素的排列长度;
进一步地,如上所述的观测***排列长度选择方法,所述步骤S101包括:
步骤1):考虑地震波在海水中***激发,经过地下某点反射/绕射后,传达到接收拖缆,该点的地震响应用地震绕射理论对该点进行定量分析,则:
水平方向某点的地震响应为:
式中,x,y,z为任意反射/绕射点的坐标,单位为m,t为地震波传播时间,单位为s,h为该点的垂直深度,单位为m,c为地震子波振幅,单位为m,单位为m/s,p为拉普拉斯变量,V为地震波速度,单位为m/s,r为震源激发点到反射/绕射点的距离,单位为m,S为反射点所在的反射界面;
步骤2)、根据反射界面位置,将该点响应分为反射波响应和绕射波响应:
式中,θ为反射点与地面的夹角,ξ为激发点到反射/绕射点距离的广义定义,单位为m;
步骤3)、根据所述反射波响应和绕射波响应,得到地表任意一点的地震响应公式为:
式中,f为地震子波主频,单位为Hz,j为虚数单位。
进一步地,如上所述的观测***排列长度选择方法,所述步骤S102包括:
(1)最大排列长度应该接近于目标层深度,需满足
|x-h|<ε (4)
式中,ε为任意无穷小的量。
(2)最大排列长度满足速度分析精度要求,即
(3)动校拉伸不大于12.5%,即排列长度小于地震波达到目标层的单程路径,即
式中,k为动校拉伸系数;
(4)考虑反射波能量和AVO精度。
进一步地,如上所述的观测***排列长度选择方法,所述步骤S103包括:
对所述公式(3)进行分析可知,当目标炮检距有响应时,该公式的最小值应大于等于0,此时应满足:
取所述公式(5)不等式中的等式部分带入公式(3),可得:
即振幅响应为与子波主频,目标层深度与速度有关的振幅为正的指数函数,恒大于0,满足速度分析精度要求;
将所述公式(6)带入目标方程(3)可得:
由于t0v<2h,得到F(f,x)>0,满足动校拉伸的要求;
反射界面入射角小于临界角时,反射能量稳定,同时考虑到保证AVO分析精度的要求,入射角为40°。
与现有技术相比,本发明可以获得以下技术效果:
本发明提供的基于波场传播理论的观测***排列长度选择方法,根据波场传播中的惠更斯原理,即任意一点的波场可以作为新的震源进行传播,可以定义一次、二次、多次震源,并对其进行定量分析,并根据互易定理,将目的层的地震响应与一次震源互换,形成理想观测***下目的层的有效响应,此时可以得到该观测***下的菲涅尔体,满足成像要求,通过以上分析,结合图1中所示的三角关系得到目标层位的最佳排列长度应为目的层深度的倍。
本发明提出的选择方法综合考虑了目标层深度,满足速度分析精度、满足动校拉伸要求、保证反射系数稳定的要求,本发明提出的排列长度的论证方法相对常规参数论证方法具有面向目标层位的更加准确、更具普适性的优点,对于海洋野外高分辨率、高信噪比三维地震勘探具有重要意义。
附图说明
图1为波场特征图;
图2为地震绕射理论对地震响应进行定量分析的原理图;
图3是本申请基于波场传播理论的观测***排列长度选择方法流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本申请提供一种基于波场传播理论的观测***目标层位最佳排列长度选择方法。所述确定观测***目标层位最佳排列长度的方法提供由地表任意点的震源激发的地震子波按照球面波场方式进行传播。
在本实施方式中,所述球面波场中任意一点按照惠更斯原理都可以作为新的震源进行传播,由此可以定义一次、二次、多次震源,波场特征。如附图1所示,震源处波场的传播和来自目的层的反射符合附图1中所示的实线和点画线关系,即一次震源激发的波场半径为一次震源到反射点的距离1,根据惠更斯原理,目的层位的反射点会形成二次震源,二次震源激发的波场半径为倍的一次震源到反射点的距离。根据互易定理,可以将目的层的地震响应与一次震源互换,形成理想观测***下目的层的有效地震响应,即附图1中所示的虚线。此时可以形成该观测***下的菲涅尔体,即附图1中所示的粗线段。根据地震绕射理论对图1中所示的理想观测***下目的层的有效地震响应进行定量分析,如附图2所示,得到地表任意位置的地震响应表达式。根据图1所示的波场特征,综合考虑目标层深度、满足速度分析精度要求(5%)、满足动校拉伸要求(12.5%)、考虑反射系数稳定等因素,通过对步骤三中地表任意位置的地震响应表达式的深入分析可以确定观测***最佳排列长度。
如图3所示,本发明实施例提供一种基于波场传播理论的观测***排列长度选择方法,包括以下步骤:
步骤S101:基于波场传播理论获取地表任意一点的地震响应公式;
具体地,该地震响应公式的获取包括以下步骤:
步骤1:考虑地震波在海水中***激发,经过地下某点反射/绕射后,传达到接收拖缆。该点的地震响应可以用地震绕射理论对该点进行定量分析。根据上述分析得到的波场特征,可以给出最大排列长度。
水平方向某点的地震响应为:
式中,x,y,z为任意反射/绕射点的坐标,单位为m,t为地震波传播时间,单位为s,h为该点的垂直深度,单位为m,c为地震子波振幅,单位为m,单位为m/s,p为拉普拉斯变量,V为地震波速度,单位为m/s,r为震源激发点到反射/绕射点的距离,单位为m,S为反射点所在的反射界面。
步骤2、根据反射界面位置,将该点响应分为反射波响应和绕射波响应:
式中,θ为反射点与地面的夹角,ξ为激发点到反射/绕射点距离的广义定义,单位为m。
步骤3、经过推导,可得地表任意位置的地震响应为:
式中,f为地震子波主频,单位为Hz,j为虚数单位。
即,基于惠更斯原理,对波场传播进行分析,可以得到地表任意一点的地震信号的特征,该特征与地下介质目标层深度,层速度和激发地震的子波主频,振幅有关。根据上述公式,我们可以综合各方面因素,更加深入、详细的给出最佳排列长度。
步骤S102:根据所述地震响应公式,从不同考虑因素确定最佳排列长度的选择准则;所述不同考虑因素包括:目标层深度、速度分析精度、动校拉伸、反射波能量、AVO精度。
具体地,从目标层深度、速度分析精度、动校拉伸、反射系数稳定等方面确定最佳排列长度的选择准则,分别为:
(1)最大排列长度应该接近于目标层深度,需满足
|x-h|<ε (4)
式中,ε为任意无穷小的量。
(2)最大排列长度满足速度分析精度要求,即
(3)动校拉伸不大于12.5%,即排列长度小于地震波达到目标层的单程路径,即
式中,k为动校拉伸系数。
(4)考虑反射波能量和AVO精度。
步骤S103:利用所述最佳排列长度的选择准则,分别得到所述不同考虑因素的排列长度;
利用所述的最佳排列长度的选择准则,对所述的地表任意位置的地震响应表达式(公式3)进行深入推导和分析,分别得到基于步骤S102所述四个选择准则的排列长度,分别为:
对所述公式(3)进行分析可知,当目标炮检距有响应时,该公式的最小值应大于等于0,此时应满足:
取所述公式(5)不等式中的等式部分带入公式(3),可得:
即振幅响应为与子波主频,目标层深度与速度有关的振幅为正的指数函数,恒大于0,满足速度分析精度要求。
将上述公式(6)带入目标方程(3)可得:
由于t0v>2h,得到F(f,x)>0,满足动校拉伸的要求。
反射界面入射角小于临界角时,反射能量稳定,同时考虑到保证AVO分析精度的要求,入射角最好为40°。
综合考虑以上所述各因素最终确定的观测***的最佳排列长度为目标层位深度的倍,该发明中确定的最佳排列长度由分析波场特征而来,基于的是惠更斯原理,即每一点都被当做新的震源进行传播,因此可以完全满足高分辨、高信噪比的地震采集***的要求。
具体地,首先通过惠更斯原理和震源互易定理得到震源的波场特征。如图1所示,根据波场波场传播理论以及三角函数关系,一次震源与目标层位的距离为1,二次震源形成的波场半径为倍的一次震源与目标层位的距离,根据震源互易定理,在一次震源处存在二次震源的互易震源,互易震源的半径为倍的一次震源与目标层位的距离,而该半径可以作为针对目标层位的最佳排列长度。然后根据地震绕射理论对波场特征进行定量分析得到地表任意位置地震响应的表达式。最后分别从速度分析精度、动校拉伸、AVO等方面对该地震响应的表达式进行具体分析,得到目标层位的排列长度均符合如图1所示的基于波场传播理论得到的目标层位的最佳排列长度。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (4)
2.根据权利要求1所述的观测***排列长度选择方法,其特征在于,所述步骤S101包括:
步骤1):考虑地震波在海水中***激发,经过地下某点反射/绕射后,传达到接收拖缆,该点的地震响应用地震绕射理论对该点进行定量分析,则:
水平方向某点的地震响应为:
式中,x,y,z为任意反射/绕射点的坐标,单位为m,t为地震波传播时间,单位为s,h为该点的垂直深度,单位为m,c为地震子波振幅,单位为m,单位为m/s,p为拉普拉斯变量,V为地震波速度,单位为m/s,r为震源激发点到反射/绕射点的距离,单位为m,S为反射点所在的反射界面;
步骤2)、根据反射界面位置,将该点响应分为反射波响应和绕射波响应:
式中,θ为反射点与地面的夹角,ξ为激发点到反射/绕射点距离的广义定义,单位为m;
步骤3)、根据所述反射波响应和绕射波响应,得到地表任意一点的地震响应公式为:
式中,f为地震子波主频,单位为Hz,j为虚数单位。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210431768.9A CN114814931A (zh) | 2022-04-23 | 2022-04-23 | 一种观测***排列长度选择方法 |
PCT/CN2022/099062 WO2023201865A1 (zh) | 2022-04-23 | 2022-06-16 | 一种观测***排列长度选择方法 |
US18/547,620 US20240094418A1 (en) | 2022-04-23 | 2022-06-16 | Selection method of array length of observation system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210431768.9A CN114814931A (zh) | 2022-04-23 | 2022-04-23 | 一种观测***排列长度选择方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114814931A true CN114814931A (zh) | 2022-07-29 |
Family
ID=82508120
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210431768.9A Pending CN114814931A (zh) | 2022-04-23 | 2022-04-23 | 一种观测***排列长度选择方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20240094418A1 (zh) |
CN (1) | CN114814931A (zh) |
WO (1) | WO2023201865A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118011471B (zh) * | 2024-01-30 | 2024-06-25 | 北京锐浪石油技术有限公司 | 一种面向多层系立体勘探的陆上三维地震采集观测***设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605158A (zh) * | 2013-10-30 | 2014-02-26 | 中国石油天然气集团公司 | 一种最大炮检距确定方法及装置 |
CN104536041A (zh) * | 2014-12-17 | 2015-04-22 | 中国石油天然气集团公司 | 一种地震观测***参数的优化方法 |
CN105223630A (zh) * | 2014-06-26 | 2016-01-06 | 中石化石油工程地球物理有限公司胜利分公司 | 基于地质模型的全方位观测***参数论证方法 |
CN107085236A (zh) * | 2017-03-23 | 2017-08-22 | 中国石油天然气集团公司 | 最大炮检距的确定方法和装置 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4611312A (en) * | 1983-02-09 | 1986-09-09 | Chevron Research Company | Method of seismic collection utilizing multicomponent receivers |
US5173880A (en) * | 1989-12-26 | 1992-12-22 | Exxon Production Research Company | Method of generating seismic wavelets using seismic range equation |
US8358562B2 (en) * | 2008-09-03 | 2013-01-22 | Christof Stork | Measuring and modifying directionality of seismic interferometry data |
CN104360388B (zh) * | 2014-10-31 | 2019-01-22 | 中国石油化工股份有限公司 | 一种三维地震观测***评价方法 |
CN107490808B (zh) * | 2017-06-14 | 2019-02-15 | 中国石油化工股份有限公司 | 一种高可靠性地震勘探观测***的建立方法 |
CN109212589A (zh) * | 2017-07-06 | 2019-01-15 | 中国石油化工股份有限公司 | 一种基于gpu/cpu协同并行的地震采集观测***设计方法 |
CN107526101B (zh) * | 2017-07-07 | 2018-12-21 | 中国地质调查局油气资源调查中心 | 一种获取地震反射波的采集和处理方法 |
US20190094397A1 (en) * | 2017-09-25 | 2019-03-28 | Real Time Geomechanics, Llc | Surface detection and location of microseismic events and earthquakes without the use of a velocity model |
CN107831533A (zh) * | 2017-11-17 | 2018-03-23 | 中国石油天然气集团公司 | 基于能量统计的采集方案检波器排列长度定量分析方法 |
-
2022
- 2022-04-23 CN CN202210431768.9A patent/CN114814931A/zh active Pending
- 2022-06-16 US US18/547,620 patent/US20240094418A1/en active Pending
- 2022-06-16 WO PCT/CN2022/099062 patent/WO2023201865A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605158A (zh) * | 2013-10-30 | 2014-02-26 | 中国石油天然气集团公司 | 一种最大炮检距确定方法及装置 |
CN105223630A (zh) * | 2014-06-26 | 2016-01-06 | 中石化石油工程地球物理有限公司胜利分公司 | 基于地质模型的全方位观测***参数论证方法 |
CN104536041A (zh) * | 2014-12-17 | 2015-04-22 | 中国石油天然气集团公司 | 一种地震观测***参数的优化方法 |
CN107085236A (zh) * | 2017-03-23 | 2017-08-22 | 中国石油天然气集团公司 | 最大炮检距的确定方法和装置 |
Non-Patent Citations (1)
Title |
---|
王超越: "煤田高密度三维地震勘探观测***面元属性评价与优选", 中国优秀硕士学位论文全文数据库 基础科学辑(月刊), 15 March 2022 (2022-03-15), pages 20 - 21 * |
Also Published As
Publication number | Publication date |
---|---|
WO2023201865A1 (zh) | 2023-10-26 |
US20240094418A1 (en) | 2024-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Duquet et al. | Kirchhoff modeling, inversion for reflectivity, and subsurface illumination | |
CN102298156B (zh) | 用于反虚反射地震数据的方法和装置 | |
Gazdag et al. | Migration of seismic data by phase shift plus interpolation | |
EP2673662B1 (en) | Method of analyzing seismic data | |
RU2570825C2 (ru) | Эффективное вычисление миграции по волновому уравнению угловых сейсмограмм | |
Farra et al. | Ray-theoretical modeling of secondary microseism P waves | |
Yokoi et al. | Consistency of the spatial autocorrelation methodwith seismic interferometry and its consequence | |
US6188964B1 (en) | Method for using global optimization to the estimation of surface-consistent residual statics | |
US8731838B2 (en) | Fresnel zone fat ray tomography | |
Wassermann et al. | Toward a single‐station approach for microzonation: Using vertical rotation rate to estimate Love‐wave dispersion curves and direction finding | |
CN106896409A (zh) | 一种基于波动方程边值反演的变深度缆鬼波压制方法 | |
CN114814931A (zh) | 一种观测***排列长度选择方法 | |
Colin et al. | Fine-scale velocity distribution revealed by datuming of very-high-resolution deep-towed seismic data: Example of a shallow-gas system from the western Black Sea | |
CN103675900A (zh) | 一种确定转换波叠前时间偏移最佳速度剖面的方法 | |
CN111103620A (zh) | 一种岩巷超前探测三维偏移成像方法 | |
Sun et al. | Reflection and diffraction separation in the dip-angle common-image gathers using convolutional neural network | |
CN105974463B (zh) | 基于源检双定向的勘探地震方法 | |
Yan et al. | Analysis of converted-wave extended images for migration velocity analysis | |
Liang et al. | Strong Ground Motion Recorded by High‐Rate GPS during the 2021 M s 6.4 Yangbi, China, Earthquake | |
US12000971B2 (en) | Method and system for seismic processing using virtual trace bins based on offset attributes and azimuthal attributes | |
Pei et al. | Data processing of the Kuiyang-ST2000 deep-towed high-resolution multichannel seismic system and application to South China Sea data | |
Wang et al. | Research on vertical cable seismic multiples imaging processing | |
Bauer et al. | Improving the resolution of wavefront tomography by utilizing diffractions | |
Duan et al. | 3D angle decomposition for elastic reverse time migration | |
Li et al. | Application of Snell's law in reflection raytracing using the multistage fast marching method |
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 |