CN112099090A - 地震资料视速度域非一致性长波长静校正方法 - Google Patents

地震资料视速度域非一致性长波长静校正方法 Download PDF

Info

Publication number
CN112099090A
CN112099090A CN201910520238.XA CN201910520238A CN112099090A CN 112099090 A CN112099090 A CN 112099090A CN 201910520238 A CN201910520238 A CN 201910520238A CN 112099090 A CN112099090 A CN 112099090A
Authority
CN
China
Prior art keywords
domain
static correction
seismic
apparent velocity
uniformity
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
Application number
CN201910520238.XA
Other languages
English (en)
Other versions
CN112099090B (zh
Inventor
张建中
任凯
杨华臣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ocean University of China
Original Assignee
Ocean University of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ocean University of China filed Critical Ocean University of China
Priority to CN201910520238.XA priority Critical patent/CN112099090B/zh
Publication of CN112099090A publication Critical patent/CN112099090A/zh
Application granted granted Critical
Publication of CN112099090B publication Critical patent/CN112099090B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/53Statics correction, e.g. weathering layer or transformation to a datum

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

地震资料视速度域非一致性长波长静校正方法
技术领域
本发明涉及一种视速度域非一致性长波长静校正方法,属于地震勘探资料处理的静校正技术领域。
背景技术
长波长静校正是地震资料处理过程中一种消除地表低速层起伏及速度变化对地震资料影响的技术。长波长静校正量主要来源于地表起伏及低速层横向速度变化,影响地震成像剖面上的构造形态。目前常用的长波长静校正技术多基于地表一致性假设,即假设地震波在低速层中的传播路径垂直于地表,这种假设在地表速度较低的地区是成立的。对于山地等复杂地区,其地表速度往往偏高,甚至有基岩出露现象,地震波在近地表层中的传播路径不符合地表一致性假设,使一致性长波长静校正不能取得令人满意的效果。
目前我国的油气勘探的重心正逐步向复杂地区转移,而静校正是获得高质量地震剖面的关键技术,地表一致性假设的误差为静校正工作带来了许多困难。波动方程基准面校正虽然可以解决一致性假设带来的问题,但是这种方法计算量大,难以应用到实际勘探过程中。目前,还缺乏易于应用的非一致性长波长静校正方法,严重制约了复杂地区地震勘探资料的处理效果。
发明内容
针对现有技术存在的上述缺陷,本发明提出了一种视速度域非一致性长波长静校正方法,其解决了地表一致性假设不适用于复杂地区的静校正问题。
本发明所述的地震资料视速度域非一致性长波长静校正方法,其特征在于,包括以下步骤:
S1:对原始地震资料的预处理,主要包括对原始数据的噪声压制、振幅补偿、废炮废道的剔除;
S2:从处理后的地震记录中拾取初至波走时信息;
S3:根据拾取的地震初至走时数据,使用菲涅尔体层析反演方法建立离散近地表速度模型v(x,z),其中 为水平方向网格点坐标, 为垂直方向网格点坐标;
S4:将共炮点地震记录变换到时间-视速度域,包括如下小步:
S41:设定变换参数vmin, vmax, nv, dv,其中,vmin为变换至视速度域后的最小视速度,vmax为变换至视速度域后的最大视速度,nv为变换至视速度域后视速度域中的道数,dv为变换至视速度域时两道之间的视速度间隔;
S42:按照设定好的参数对地震记录进行变换,得到视速度域数据g(va,τ),变换公式如下:
Figure DEST_PATH_IMAGE001
(1)
其中,g(va,τ)为视速度域中数据;τ为视速度域中时间;va为视速度,范围为(vmin,vmax),视速度间隔与变换后视速度个数分别为dv和nv;s(x,t)为原始时间域地震记录,x为偏移距,t为时间;R{•}表示从时间域到视速度域的映射关系;(h0,t0)为映射原点,一般设置h0= 0, t0 = 0;
S5:计算炮点端的静校正量Δt,并在视速度域中进行炮点端非一致性静校正,包括如下小步:
S51:对视速度域中每一视速度va,根据地表速度模型,自炮点出发,沿当前视速度对应角度方向计算对应的时差Δtd
Figure DEST_PATH_IMAGE002
(2)
其中,zn为低速层底界面对应的网格点深度;zi为当前计算的网格点深度;dz为垂直方向网格点间隔;xs为该炮对应的炮点水平方向坐标;θ为地震波出射角度;v(xs+zi*tanθ,zi)为沿当前射线方向的第i个网格点的速度;v(xs,0)为炮点处地表速度;
S52:对视速度域中每一视速度对应的角度θ,计算自低速层底界面延拓到基准面的时差 Δtu,并与Δtd相加得到炮点端非一致性静校正量Δt:
Figure DEST_PATH_IMAGE003
(3)
Figure DEST_PATH_IMAGE004
(4)
其中,zh为基准面高程,vp为替换速度;
S53:用计算得到的炮点端非一致性静校正量Δt对相应视速度域的地震记录进行校正,然后计算下一炮对应的炮点端非一致性静校正量,重复步骤S51-S53,直至所有炮点均计算完毕;
S6:将校正后的视速度域地震记录反变换至时间-空间域,得到炮点端非一致性静校正后的地震记录s'(h,t),变换公式如下:
Figure DEST_PATH_IMAGE005
(5)
其中,R'{•}为反变换映射关系,由于时间-空间域中偏移距信息h是由观测系统决定的,因此需对反变换后得到的地震记录进行双线性插值处理得到观测系统与原始地震记录相同的记录;
S7:对进行过炮点端非一致性静校正后的地震记录s'(h,t)抽取共检波点道集,得到检波点域地震记录r(h,t);
S8:根据互换原理,检波点域地震记录可等效视为自检波点位置激发,炮点位置接收的炮集,对检波点域记录重复步骤3到5,得到进行检波点端非一致性静校正后的地震记录 ;
S9:将检波点域地震记录变换到炮点域,得到非一致性长波长静校正后的地震记录。
优选地,步骤二中使用菲涅尔体层析方法建立离散近地表速度模型。
优选地,步骤五中使用双线性插值方法获取观测系统与原始地震记录相同的记录。
本发明的有益效果是:本发明所述的视速度域非一致性长波长静校正方法,在视速度域中分别进行炮点与检波点端静校正,考虑到了地震波在近地表低速层中传播路径非垂直时的影响,解决了基于地表一致性假设的静校正方法无法用于复杂地区地震资料静校正的问题;本发明易于计算机自动化实现,相比现有波动方程基准面校正等其他非一致性静校正方法,具有计算速度快的优点,为近地表复杂时的地震资料处理提供了一种有效的非一致性静校正方法。
附图说明
图1为本发明的流程图;
图2为理论模型示意图;
图3为用理论模型合成的第一炮地震记录;
图4为对第一炮地震记录进行一致性静校正结果与理论结果对比;
图5为对第一炮地震记录进行非一致性静校正结果与理论结果对比。
具体实施方式
为了使本发明目的、技术方案更加清楚明白,下面结合实施例,对本发明作进一步详细说明。
实施例1:
S1:对原始地震资料的预处理,主要包括对原始数据的噪声压制、振幅补偿、废炮废道的剔除;
S2:从处理后的地震记录中拾取初至波走时信息;
S3:根据拾取的地震初至走时数据,使用菲涅尔体层析反演方法建立离散近地表速度模型v(x,z),其中 为水平方向网格点坐标, 为垂直方向网格点坐标;
S4:将共炮点地震记录变换到时间-视速度域,包括如下小步:
S41:设定变换参数vmin, vmax, nv, dv,其中,vmin为变换至视速度域后的最小视速度,vmax为变换至视速度域后的最大视速度,nv为变换至视速度域后视速度域中的道数,dv为变换至视速度域时两道之间的视速度间隔;
S42:按照设定好的参数对地震记录进行变换,得到视速度域数据g(va,τ),变换公式如下:
Figure 163433DEST_PATH_IMAGE001
(1)
其中,g(va,τ)为视速度域中数据;τ为视速度域中时间;va为视速度,范围为(vmin,vmax),视速度间隔与变换后视速度个数分别为dv和nv;s(x,t)为原始时间域地震记录,x为偏移距,t为时间;R{•}表示从时间域到视速度域的映射关系;(h0,t0)为映射原点,一般设置h0= 0, t0 = 0;
S5:计算炮点端的静校正量Δt,并在视速度域中进行炮点端非一致性静校正,包括如下小步:
S51:对视速度域中每一视速度va,根据地表速度模型,自炮点出发,沿当前视速度对应角度方向计算对应的时差Δtd
Figure 533104DEST_PATH_IMAGE002
(2)
其中,zn为低速层底界面对应的网格点深度;zi为当前计算的网格点深度;dz为垂直方向网格点间隔;xs为该炮对应的炮点水平方向坐标;θ为地震波出射角度;v(xs+zi*tanθ,zi)为沿当前射线方向的第i个网格点的速度;v(xs,0)为炮点处地表速度;
S52:对视速度域中每一视速度对应的角度θ,计算自低速层底界面延拓到基准面的时差 Δtu,并与Δtd相加得到炮点端非一致性静校正量Δt:
Figure DEST_PATH_IMAGE006
(3)
Figure DEST_PATH_IMAGE007
(4)
其中,zh为基准面高程,vp为替换速度;
S53:用计算得到的炮点端非一致性静校正量Δt对相应视速度域的地震记录进行校正,然后计算下一炮对应的炮点端非一致性静校正量,重复步骤S51-S53,直至所有炮点均计算完毕;
S6:将校正后的视速度域地震记录反变换至时间-空间域,得到炮点端非一致性静校正后的地震记录s'(h,t),变换公式如下:
Figure 437475DEST_PATH_IMAGE005
(5)
其中,R'{•}为反变换映射关系,由于时间-空间域中偏移距信息h是由观测系统决定的,因此需对反变换后得到的地震记录进行双线性插值处理得到观测系统与原始地震记录相同的记录;
S7:对进行过炮点端非一致性静校正后的地震记录s'(h,t)抽取共检波点道集,得到检波点域地震记录r(h,t);
S8:根据互换原理,检波点域地震记录可等效视为自检波点位置激发,炮点位置接收的炮集,对检波点域记录重复步骤3到5,得到进行检波点端非一致性静校正后的地震记录 ;
S9:将检波点域地震记录变换到炮点域,得到非一致性长波长静校正后的地震记录。
实施例2:
下面结合具体实施方式,对本发明进行解释和说明。
为了进一步说明该方法的实现思路、实现过程及有效性,用理论模型进行了测试,并与地表一致性长波长静校正的结果进行比较。
S1:由于使用理论模型进行测试,因此不需要进行资料预处理与速度建模,直接使用理论模型进行数值计算,得到该模型的地震记录;
S2:将地震记录变换至视速度域;
S3:计算炮点端非一致性静校正量,并进行炮点端静校正;
S4:将校正后的视速度域地震记录反变换到时间域;
S5:对变换到时间域的地震记录,抽取共检波点道集;
S6:将共检波点道集地震记录再变换至视速度域;
S7:计算检波点端非一致性静校正量,并进行检波点端静校正;
S8:将视速度域地震记录变换到时间域,再变换回共炮点道集,即可获得非一致性静校正后的地震记录。
图2是本实施例所用速度模型,图3是未做静校正前第一炮地震记录,图4和图5分别是对第一炮地震记录进行一致性静校正与非一致性静校正后的效果图,图中黑色虚线为无地表低速层时地震记录同相轴形态,即静校正后理论同相轴。从图可以看出,进行一致性静校正后,地震记录同相轴与理论同相轴形态有较大差别,而进行非一致性静校正后,地震记录同相轴更接近于理论形态,静校正效果得到了明显改善。这表明本发明对地表低速层不满足地表一致性假设的地区的适应性和良好的应用效果。
本发明可广泛运用于地震资料处理的静校正环节。
以上所述仅为本发明的较佳实施例而己,并不以本发明为限制,凡在本发明的精神和原则之内所作的均等修改、等同替换和改进等,均应包含在本发明的专利涵盖范围内。

Claims (4)

1.一种地震资料视速度域非一致性长波长静校正方法,其特征在于,包括以下步骤:
S1:对原始地震资料的预处理,主要包括对原始数据的噪声压制、振幅补偿、废炮废道的剔除;
S2:从处理后的地震记录中拾取初至波走时信息;
S3:根据拾取的地震初至走时数据,使用菲涅尔体层析反演方法建立离散近地表速度模型v(x,z),其中 为水平方向网格点坐标, 为垂直方向网格点坐标;
S4:将共炮点地震记录变换到时间-视速度域,包括如下小步:
S41:设定变换参数vmin, vmax, nv, dv,其中,vmin为变换至视速度域后的最小视速度,vmax为变换至视速度域后的最大视速度,nv为变换至视速度域后视速度域中的道数,dv为变换至视速度域时两道之间的视速度间隔;
S42:按照设定好的参数对地震记录进行变换,得到视速度域数据g(va,τ),变换公式如下:
Figure 609810DEST_PATH_IMAGE001
(1)
其中,g(va,τ)为视速度域中数据;τ为视速度域中时间;va为视速度,范围为(vmin,vmax),视速度间隔与变换后视速度个数分别为dv和nv;s(x,t)为原始时间域地震记录,x为偏移距,t为时间;R{•}表示从时间域到视速度域的映射关系;(h0,t0)为映射原点,一般设置h0 = 0,t0 = 0;
S5:计算炮点端的静校正量Δt,并在视速度域中进行炮点端非一致性静校正,包括如下小步:
S51:对视速度域中每一视速度va,根据地表速度模型,自炮点出发,沿当前视速度对应角度方向计算对应的时差Δtd
Figure 663960DEST_PATH_IMAGE002
(2)
其中,zn为低速层底界面对应的网格点深度;zi为当前计算的网格点深度;dz为垂直方向网格点间隔;xs为该炮对应的炮点水平方向坐标;θ为地震波出射角度;v(xs+zi*tanθ,zi)为沿当前射线方向的第i个网格点的速度;v(xs,0)为炮点处地表速度;
S52:对视速度域中每一视速度对应的角度θ,计算自低速层底界面延拓到基准面的时差 Δtu,并与Δtd相加得到炮点端非一致性静校正量Δt:
Figure 230071DEST_PATH_IMAGE003
(3)
Figure 531739DEST_PATH_IMAGE004
(4)
其中,zh为基准面高程,vp为替换速度;
S53:用计算得到的炮点端非一致性静校正量Δt对相应视速度域的地震记录进行校正,然后计算下一炮对应的炮点端非一致性静校正量,重复步骤S51-S53,直至所有炮点均计算完毕;
S6:将校正后的视速度域地震记录反变换至时间-空间域,得到炮点端非一致性静校正后的地震记录s'(h,t),变换公式如下:
Figure 892313DEST_PATH_IMAGE005
(5)
其中,R'{•}为反变换映射关系,由于时间-空间域中偏移距信息h是由观测系统决定的,因此需对反变换后得到的地震记录进行双线性插值处理得到观测系统与原始地震记录相同的记录;
S7:对进行过炮点端非一致性静校正后的地震记录s'(h,t)抽取共检波点道集,得到检波点域地震记录r(h,t);
S8:根据互换原理,检波点域地震记录可等效视为自检波点位置激发,炮点位置接收的炮集,对检波点域记录重复步骤3到5,得到进行检波点端非一致性静校正后的地震记录;
S9:将检波点域地震记录变换到炮点域,得到非一致性长波长静校正后的地震记录。
2.根据权利要求1所述的视速度域非一致性长波长静校正方法,其特征在于,所述S4中,将原始地震记录变换至视速度域。
3.根据权利要求1所述的地震资料视速度域非一致性长波长静校正方法,其特征在于,所述S5中,需计算不同角度对应的静校正量。
4.根据权利要求1所述的地震资料视速度域非一致性长波长静校正方法,其特征在于,所述方法需在炮点端和检波点端分别计算非一致性静校正量并在视速度域进行静校正,再从视速度域反变换到时间域。
CN201910520238.XA 2019-06-17 2019-06-17 地震资料视速度域非一致性长波长静校正方法 Active CN112099090B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910520238.XA CN112099090B (zh) 2019-06-17 2019-06-17 地震资料视速度域非一致性长波长静校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910520238.XA CN112099090B (zh) 2019-06-17 2019-06-17 地震资料视速度域非一致性长波长静校正方法

Publications (2)

Publication Number Publication Date
CN112099090A true CN112099090A (zh) 2020-12-18
CN112099090B CN112099090B (zh) 2021-07-30

Family

ID=73748292

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910520238.XA Active CN112099090B (zh) 2019-06-17 2019-06-17 地震资料视速度域非一致性长波长静校正方法

Country Status (1)

Country Link
CN (1) CN112099090B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113640878A (zh) * 2021-08-12 2021-11-12 西南石油大学 利用虚拟震源扫描构建方位角-视速度雷达图的方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101630015A (zh) * 2008-07-16 2010-01-20 中国石油天然气集团公司 一种提高初至波拾取精度和效率的方法
CN102944896A (zh) * 2012-11-09 2013-02-27 中国石油集团川庆钻探工程有限公司地球物理勘探公司 表层调查数据的模型法静校正方法
CN104330826A (zh) * 2014-07-31 2015-02-04 中国石油化工股份有限公司 一种去除复杂地表条件下多种噪音的方法
KR20150123550A (ko) * 2014-04-25 2015-11-04 주식회사 코어밸런스 도시건축물에 적용되는 지진데이터 측정 및 감시 시스템 및 방법
CN105093320A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 针对高速结晶盐壳覆盖区层析静校正初至拾取方法
CN107561589A (zh) * 2017-10-25 2018-01-09 中国石油化工股份有限公司 一种近地表横波层速度模型建立方法
CN108279436A (zh) * 2017-12-20 2018-07-13 中国石油天然气股份有限公司 一种拾取地震波初至的方法及装置
CN108387936A (zh) * 2018-03-20 2018-08-10 中国石油集团东方地球物理勘探有限责任公司 一种可控震源资料初至高效拾取方法及层析静校正方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101630015A (zh) * 2008-07-16 2010-01-20 中国石油天然气集团公司 一种提高初至波拾取精度和效率的方法
CN102944896A (zh) * 2012-11-09 2013-02-27 中国石油集团川庆钻探工程有限公司地球物理勘探公司 表层调查数据的模型法静校正方法
KR20150123550A (ko) * 2014-04-25 2015-11-04 주식회사 코어밸런스 도시건축물에 적용되는 지진데이터 측정 및 감시 시스템 및 방법
CN105093320A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 针对高速结晶盐壳覆盖区层析静校正初至拾取方法
CN104330826A (zh) * 2014-07-31 2015-02-04 中国石油化工股份有限公司 一种去除复杂地表条件下多种噪音的方法
CN107561589A (zh) * 2017-10-25 2018-01-09 中国石油化工股份有限公司 一种近地表横波层速度模型建立方法
CN108279436A (zh) * 2017-12-20 2018-07-13 中国石油天然气股份有限公司 一种拾取地震波初至的方法及装置
CN108387936A (zh) * 2018-03-20 2018-08-10 中国石油集团东方地球物理勘探有限责任公司 一种可控震源资料初至高效拾取方法及层析静校正方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨华臣等: "地震静校正的相对时延法 ", 《中国海洋大学学报(自然科学版)》 *
金昌昆等: "共接收面元三维折射旅行时反演方法 ", 《石油地球物理勘探》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113640878A (zh) * 2021-08-12 2021-11-12 西南石油大学 利用虚拟震源扫描构建方位角-视速度雷达图的方法
CN113640878B (zh) * 2021-08-12 2024-03-29 西南石油大学 利用虚拟震源扫描构建方位角-视速度雷达图的方法

Also Published As

Publication number Publication date
CN112099090B (zh) 2021-07-30

Similar Documents

Publication Publication Date Title
CN113625337B (zh) 一种极浅水高精度地震资料快速成像方法
CN108196305B (zh) 一种山地静校正方法
CN109765616B (zh) 一种保幅波场延拓校正方法及系统
CN103293554B (zh) 一种基于地震数据的静校正方法
CN112327362B (zh) 速度域的海底多次波预测与追踪衰减方法
CN105204065A (zh) 一种拾取初至波的方法和装置
CN112099090B (zh) 地震资料视速度域非一致性长波长静校正方法
CN109490961B (zh) 起伏地表无射线追踪回折波层析成像方法
CN113075732B (zh) 一种消除高分辨率小多道地震地层异常起伏的方法
CN111650645B (zh) 一种变偏移距vsp弯线校正处理方法和装置
CN112213784B (zh) 复杂地表地震数据一次处理快速静校正方法
CN111983685B (zh) τ-p域地表非一致性静校正方法
CN112630838B (zh) 一种时间域浮动面转换深度域计算方法及计算系统
CN111751874B (zh) 一种变偏移距vsp叠后变覆盖次数校正方法和装置
CN114371507A (zh) 一种基于一致性谱相关成像的震源定位方法及系统
CN111929731B (zh) 地表一致性和非一致性联合静校正方法
CN116840888A (zh) 一种采集资料插值的理论观测系统生成方法
CN109283582B (zh) 一种静校正的补偿方法
CN112526607B (zh) 一种基于最小二乘反演的高频静校正量提取方法
CN113917539B (zh) 火山岩覆盖区地震数据叠前道集处理方法、系统和装置
CN112269211A (zh) τ-p域地表非一致性长波长静校正方法
CN112327361B (zh) 基于线性同相轴迭代追踪衰减的倾斜干扰剔除方法
CN113820748B (zh) 一种地震sv波发生全反射时的地面滑行波时程获取方法
CN112083489B (zh) 一种基于多信息约束的叠前深度偏移速度更新方法
CN114397699A (zh) 一种起伏地表q偏移地震成像方法、装置及终端设备

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