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

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

Info

Publication number
CN112099090B
CN112099090B CN201910520238.XA CN201910520238A CN112099090B CN 112099090 B CN112099090 B CN 112099090B CN 201910520238 A CN201910520238 A CN 201910520238A CN 112099090 B CN112099090 B CN 112099090B
Authority
CN
China
Prior art keywords
domain
static correction
seismic
apparent velocity
shot
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.)
Expired - Fee Related
Application number
CN201910520238.XA
Other languages
English (en)
Other versions
CN112099090A (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
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/362—Effecting static or dynamic corrections; Stacking
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/53—Statics 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)，其中x为水平方向网格点坐标，z为垂直方向网格点坐标；
S4：将共炮点地震记录变换到时间-视速度域，包括如下小步：
S41：设定变换参数vmin,vmax,nv,dv，其中，vmin为变换至视速度域后的最小视速度，vmax为变换至视速度域后的最大视速度，nv为变换至视速度域后视速度域中的道数，dv为变换至视速度域时两道之间的视速度间隔；
S42：按照设定好的参数对地震记录进行变换，得到视速度域数据g(va,τ)，变换公式如下：
Figure GDA0003032909680000021
其中，g(va,τ)为视速度域中数据；τ为视速度域中时间；va为视速度，范围为(vmin,vmax)，视速度间隔与变换后视速度个数分别为dv和nv；s(x,t)为原始时间域地震记录，h为偏移距，t为时间； R{·}表示从时间域到视速度域的映射关系；(h0,t0)为映射原点，取地震记录中最接近炮点的地震道对应的偏移距为h0,该道对应的地震波初至时间为t0；
S5：计算炮点端的静校正量Δt，并在视速度域中进行炮点端非一致性静校正，包括如下小步：
S51：对视速度域中每一视速度va，根据地表速度模型，自炮点出发，沿当前视速度对应角度方向计算对应的时差Δtd：
Figure GDA0003032909680000022
其中，zn为低速层底界面对应的网格点深度；zi为当前计算的网格点深度；dz为垂直方向网格点间隔；xs为该炮对应的炮点水平方向坐标；θ为地震波出射角度；v(xs+zi*tanθ,zi)为沿当前射线方向的第i个网格点的速度；v(xs,0)为炮点处地表速度；
S52：对视速度域中每一视速度对应的角度θ，计算自低速层底界面延拓到基准面的时差Δ tu，并与Δtd相加得到炮点端非一致性静校正量Δt：
Figure GDA0003032909680000023
Δt＝Δtd+Δtu (4)
其中，zh为基准面高程，vp为替换速度；
S53：用计算得到的炮点端非一致性静校正量Δt对相应视速度域的地震记录进行校正，然后计算下一炮对应的炮点端非一致性静校正量，重复步骤S51-S53，直至所有炮点均计算完毕；
S6：将校正后的视速度域地震记录反变换至时间-空间域，得到炮点端非一致性静校正后的地震记录s'(h,t)，变换公式如下：
Figure GDA0003032909680000031
其中R'{·}为反变换映射关系，由于时间-空间域中偏移距信息h是由观测系统决定的，因此需对反变换后得到的地震记录进行双线性插值处理得到观测系统与原始地震记录相同的记录；
S7：对进行过炮点端非一致性静校正后的地震记录s'(h,t)抽取共检波点道集，得到检波点域地震记录r(h,t)；
S8：根据互换原理，检波点域地震记录可等效视为自检波点位置激发，炮点位置接收的炮集，对检波点域记录重复步骤3到5，得到进行检波点端非一致性静校正后的地震记录r'(h,t)；
S9:将检波点域地震记录变换到炮点域，得到非一致性长波长静校正后的地震记录。
优选地，步骤二中使用菲涅尔体层析方法建立离散近地表速度模型。
优选地，步骤五中使用双线性插值方法获取观测系统与原始地震记录相同的记录。
本发明的有益效果是：本发明所述的视速度域非一致性长波长静校正方法，在视速度域中分别进行炮点与检波点端静校正，考虑到了地震波在近地表低速层中传播路径非垂直时的影响，解决了基于地表一致性假设的静校正方法无法用于复杂地区地震资料静校正的问题；本发明易于计算机自动化实现，相比现有波动方程基准面校正等其他非一致性静校正方法，具有计算速度快的优点，为近地表复杂时的地震资料处理提供了一种有效的非一致性静校正方法。
附图说明
图1为本发明的流程图；
图2为理论模型示意图；
图3为用理论模型合成的第一炮地震记录；
图4为对第一炮地震记录进行一致性静校正结果与理论结果对比；
图5为对第一炮地震记录进行非一致性静校正结果与理论结果对比；
具体实施方式
为了使本发明目的、技术方案更加清楚明白，下面结合实施例，对本发明作进一步详细说明。
实施例1：
S1：对原始地震资料的预处理，主要包括对原始数据的噪声压制、振幅补偿、废炮废道的剔除；
S2:从处理后的地震记录中拾取初至波走时信息；
S3：根据拾取的地震初至走时数据，使用菲涅尔体层析反演方法建立离散近地表速度模型 v(x,z)，其中x为水平方向网格点坐标，z为垂直方向网格点坐标；
S4：将共炮点地震记录变换到时间-视速度域，包括如下小步：
S41：设定变换参数vmin,vmax,nv,dv，其中，vmin为变换至视速度域后的最小视速度，vmax为变换至视速度域后的最大视速度，nv为变换至视速度域后视速度域中的道数，dv为变换至视速度域时两道之间的视速度间隔；
S42：按照设定好的参数对地震记录进行变换，得到视速度域数据g(va,τ)，变换公式如下：
Figure GDA0003032909680000041
其中，g(va,τ)为视速度域中数据；τ为视速度域中时间；va为视速度，范围为(vmin,vmax)，视速度间隔与变换后视速度个数分别为dv和nv；s(x,t)为原始时间域地震记录，h为偏移距，t为时间； R{·}表示从时间域到视速度域的映射关系；(h0,t0)为映射原点，一般设置h0＝0,t0＝0；
S5：计算炮点端的静校正量Δt，并在视速度域中进行炮点端非一致性静校正，包括如下小步：
S51：对视速度域中每一视速度va，根据地表速度模型，自炮点出发，沿当前视速度对应角度方向计算对应的时差Δtd：
Figure GDA0003032909680000042
其中，zn为低速层底界面对应的网格点深度；zi为当前计算的网格点深度；dz为垂直方向网格点间隔；xs为该炮对应的炮点水平方向坐标；θ为地震波出射角度；v(xs+zi*tanθ,zi)为沿当前射线方向的第i个网格点的速度；v(xs,0)为炮点处地表速度；
S52：对视速度域中每一视速度对应的角度θ，计算自低速层底界面延拓到基准面的时差Δ tu，并与Δtd相加得到炮点端非一致性静校正量Δt：
Figure GDA0003032909680000043
Δt＝Δtd+Δtu (4)
其中，zh为基准面高程，vp为替换速度；
S53：用计算得到的炮点端非一致性静校正量Δt对相应视速度域的地震记录进行校正，然后计算下一炮对应的炮点端非一致性静校正量，重复步骤S51-S53，直至所有炮点均计算完毕；
S6：将校正后的视速度域地震记录反变换至时间-空间域，得到炮点端非一致性静校正后的地震记录s'(h,t)，变换公式如下：
Figure GDA0003032909680000051
其中R'{·}为反变换映射关系，由于时间-空间域中偏移距信息h是由观测系统决定的，因此需对反变换后得到的地震记录进行双线性插值处理得到观测系统与原始地震记录相同的记录；
S7：对进行过炮点端非一致性静校正后的地震记录s'(h,t)抽取共检波点道集，得到检波点域地震记录r(h,t)；
S8：根据互换原理，检波点域地震记录可等效视为自检波点位置激发，炮点位置接收的炮集，对检波点域记录重复步骤3到5，得到进行检波点端非一致性静校正后的地震记录r'(h,t)；
S9:将检波点域地震记录变换到炮点域，得到非一致性长波长静校正后的地震记录。
实施例2：
下面结合具体实施方式，对本发明进行解释和说明。
为了进一步说明该方法的实现思路、实现过程及有效性，用理论模型进行了测试，并与地表一致性长波长静校正的结果进行比较。
S1：由于使用理论模型进行测试，因此不需要进行资料预处理与速度建模，直接使用理论模型进行数值计算，得到该模型的地震记录；
S2：将地震记录变换至视速度域；
S3：计算炮点端非一致性静校正量，并进行炮点端静校正；
S4：将校正后的视速度域地震记录反变换到时间域；
S5：对变换到时间域的地震记录，抽取共检波点道集；
S6：将共检波点道集地震记录再变换至视速度域；
S7：计算检波点端非一致性静校正量，并进行检波点端静校正；
S8：将视速度域地震记录变换到时间域，再变换回共炮点道集，即可获得非一致性静校正后的地震记录。
图2是本实施例所用速度模型，图3是未做静校正前第一炮地震记录，图4和图5分别是对第一炮地震记录进行一致性静校正与非一致性静校正后的效果图，图中黑线为无地表低速层时地震记录同相轴形态。从图可以看出，进行一致性静校正后，地震记录同相轴与理论同相轴形态有较大差别，而进行非一致性静校正后，地震记录同相轴更接近于理论形态，静校正效果得到了明显改善。这表明本发明对地表低速层不满足地表一致性假设的地区的适应性和良好的应用效果。
本发明可广泛运用于地震资料处理的静校正环节。
以上所述仅为本发明的较佳实施例而己，并不以本发明为限制，凡在本发明的精神和原则之内所作的均等修改、等同替换和改进等，均应包含在本发明的专利涵盖范围内。

Claims (4)

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

Family

ID=73748292

Family Applications (1)

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

Country Status (1)

Country Link
CN (1) CN112099090B (zh)

Families Citing this family (1)

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

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
共接收面元三维折射旅行时反演方法;金昌昆等;《石油地球物理勘探》;20160815(第04期);全文 *
地震静校正的相对时延法;杨华臣等;《中国海洋大学学报(自然科学版)》;20170315(第03期);全文 *

Also Published As

Publication number Publication date
CN112099090A (zh) 2020-12-18

Similar Documents

Publication Publication Date Title
CN113625337B (zh) 一种极浅水高精度地震资料快速成像方法
CN108196305B (zh) 一种山地静校正方法
CN103293554B (zh) 一种基于地震数据的静校正方法
CN109765616B (zh) 一种保幅波场延拓校正方法及系统
CN112327362B (zh) 速度域的海底多次波预测与追踪衰减方法
CN112099090B (zh) 地震资料视速度域非一致性长波长静校正方法
CN109490961B (zh) 起伏地表无射线追踪回折波层析成像方法
CN113075732B (zh) 一种消除高分辨率小多道地震地层异常起伏的方法
CN113406700B (zh) 一种地表主动源反射波干涉成像方法
CN111650645B (zh) 一种变偏移距vsp弯线校正处理方法和装置
CN112213784B (zh) 复杂地表地震数据一次处理快速静校正方法
CN111538088B (zh) 一种海上斜缆波场校正方法
CN109283582B (zh) 一种静校正的补偿方法
CN111983685B (zh) τ-p域地表非一致性静校正方法
CN111751874B (zh) 一种变偏移距vsp叠后变覆盖次数校正方法和装置
CN111929731B (zh) 地表一致性和非一致性联合静校正方法
CN114371507A (zh) 一种基于一致性谱相关成像的震源定位方法及系统
CN112630838B (zh) 一种时间域浮动面转换深度域计算方法及计算系统
CN113917539B (zh) 火山岩覆盖区地震数据叠前道集处理方法、系统和装置
CN112526607B (zh) 一种基于最小二乘反演的高频静校正量提取方法
CN113075734B (zh) 一种基于信噪比约束的剩余曲率谱计算方法及装置
CN116840888A (zh) 一种采集资料插值的理论观测系统生成方法
CN112269211A (zh) τ-p域地表非一致性长波长静校正方法
CN114397699A (zh) 一种起伏地表q偏移地震成像方法、装置及终端设备
CN116106970A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210730