CN113376695B - 一种适用于煤层底板复杂陷落柱的全波形反演方法 - Google Patents
一种适用于煤层底板复杂陷落柱的全波形反演方法 Download PDFInfo
- Publication number
- CN113376695B CN113376695B CN202110650586.6A CN202110650586A CN113376695B CN 113376695 B CN113376695 B CN 113376695B CN 202110650586 A CN202110650586 A CN 202110650586A CN 113376695 B CN113376695 B CN 113376695B
- Authority
- CN
- China
- Prior art keywords
- field model
- velocity field
- coal seam
- velocity
- 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.)
- Active
Links
- 239000003245 coal Substances 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000002939 conjugate gradient method Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000005065 mining Methods 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000003384 imaging method Methods 0.000 description 10
- 238000001514 detection method Methods 0.000 description 8
- 238000003325 tomography Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000012795 verification 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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- 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/20—Arrangements of receiving elements, e.g. geophone pattern
-
- 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/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1429—Subsurface, e.g. in borehole or below weathering layer or mud line
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种适用于煤层底板复杂陷落柱的全波形反演方法,先获取震源的观测数据,然后经过计算获得震源的正传波场和伴随源的反传波场,然后确定速度场模型的梯度表达式,接着根据正传波场和反传波场及梯度计算表达式,计算得到速度场模型参数的梯度,有利于获得准确的速度场模型参数的更新方向,最终通过对速度场模型参数进行迭代更新,直至数据残差达到设定阈值或迭代更新次数达到预设值,停止迭代更新确定速度场模型参数,此时采用该参数的速度场模型进行全波形反演。因此本发明能对煤层底板下方条件复杂的陷落柱进行精确定位及形态探测,从而为后续开采提供数据支撑。
Description
技术领域
本发明涉及一种适用于煤层底板复杂陷落柱的全波形反演方法,属于地震勘探技术领域。
背景技术
陷落柱是煤系下伏灰岩中溶洞顶部覆岩发生塌陷所形成的陷落体,是我国华北型煤田中一种典型的隐伏垂向构造。因具有随机性及不确定性,其诱发的突水隐蔽性强、速度快、水量大、破坏性强,对煤矿安全生产的危害极大。地震三维勘探在绿色煤炭资源丰富的中西部地区受地表地形以及浅部上组煤的开采和多层采空区及地表沉陷的影响,陷落柱的探测精度受到了很大限制,矿井地震勘探具有距离目标体近、地震波能量和高频成分衰减少、地震波传播路径干扰少、分辨率高等特点。因此,开展井下针对陷落柱探测具有重要意义。
目前工作面回采过程中针对煤层陷落柱的探测以槽波最为典型。当陷落柱发育至煤层时,会破坏煤层的连续性,利用槽波的地震波属性进行探测,典型的探测方式是利用工作面其中一条巷道布置检波器,另外一条巷道进行震源激发进行CT探测。自Evison在新西兰首次于煤层中激发并且记录到了槽波,随后众多学者开展了大量研究,目前,槽波勘探已成为煤层工作面回采面内构造探测的常规工作。
在槽波发育的地区,槽波勘探作为主要的煤层勘探方法,但是,若煤层相对较薄或者煤层槽波不发育时,槽波则不能作为探测手段,此时发现沿着顶底板传播的折射波具有陷落柱探测可能性。因此,当单一类型的地震波成像精度受限,开展全波场信息的综合使用可以提升成像分辨率。为了获得煤田高精度的速度模型,全波形反演方法在煤田中的应用已经有了一些初探。有学者将全波形反演方法应用到煤田多参数反演中,并且全波形反演结果与测井曲线的吻合度比较高。还有一些学者在煤矿井下针对煤层内地质构造高精度成像问题,发展了槽波全波形反演方法,相比于常规的槽波成像方法,利用全波形反演可以获得更高精度的反演结果。但是目前研究仅限于合成数据的参数反演,同时现场陷落柱条件复杂,目前的全波形反演方法无法对条件复杂多个陷落柱的精确反演,因此如何能对条件复杂的陷落柱进行精确的全波形反演是亟需解决的问题。
发明内容
针对上述现有技术存在的问题,本发明提供一种适用于煤层底板复杂陷落柱的全波形反演方法,能对煤层底板下方条件复杂的陷落柱进行精确定位及形态探测,从而为后续开采提供数据支撑。
为了实现上述目的,本发明采用的技术方案是:一种适用于煤层底板复杂陷落柱的全波形反演方法,该方法的具体步骤为:
A、在煤层工作面一侧的巷道壁相互等间距布设多个震源;在煤层工作面另一侧的巷道壁相互等间距布设多个检波器,多个检波器均通过连接主线与采集基站及采集主机连接,使检波器、采集基站及采集主机形成地震观测系统;激发多个震源,通过地震观测系统接收观测数据并获得各个地震子波;
B、先构建基于卷积型的目标函数,具体定义如下:
其中合成数据中的正传波场通过弹性波动方程合成,具体公式如下:
L[m(x)]u(x,t;m)=s(xs,t) (2)
其中,xs是震源坐标,L[]是弹性波正演模拟算子,u(x,t;m)是正传波场,s(xs,t)为震源;
然后将步骤A获得地震子波及震源坐标代入公式(2)计算得出震源的正传波场;
C、利用拉格朗日伴随方法获得模型参数的梯度,弹性波动方程的伴随方程为:
根据公式(4)计算得出基于卷积型目标函数对应的伴随源,然后将求取的伴随源代入公式(3)中进行数据模拟,最终获得伴随源的反转波场;
D、先确定速度场模型参数m=(vp,vs)的梯度表达式如下:
将步骤B计算出的正传波场和步骤C计算出的反传波场代入到公式(5)中,最终计算得到速度场模型参数的梯度;
E、采用共轭梯度方法求取更新方向,根据得到的梯度对速度场模型参数进行迭代更新,直至数据残差达到设定阈值或迭代更新次数达到预设值,停止迭代更新确定速度场模型参数,此时采用该参数的速度场模型进行全波形反演。
进一步,所述步骤E的具体过程为:
采用共轭梯度方法求取更新方向Δmk+1,
其中下标k表示迭代次数,sk是用于计算共轭方向的一个中间变量,上标T表示矩阵转置,更新后的速度场模型为:
其中,tk表示选取的更新步长,tkΔmk+1即为此次迭代过程中速度场模型参数的更新量;
通过对速度场模型参数进行迭代更新,直至数据残差达到设定阈值或迭代更新次数达到预设值,停止迭代更新确定速度场模型参数,此时采用该参数的速度场模型进行全波形反演。
进一步,所述步骤A中在煤层工作面的切眼处等间距设有一排检波器,使其与煤层工作面另一侧的巷道壁内设置的检波器形成L型观测系统。在切眼处增设检波器能提高观测系统采集地震数据的精确性,从而使后续反演的精度更高。
与现有技术相比,本发明先获取震源的观测数据,然后经过计算获得震源的正传波场和伴随源的反传波场,然后确定速度场模型的梯度表达式,接着根据正传波场和反传波场及梯度计算表达式,计算得到速度场模型参数的梯度,有利于获得准确的速度场模型参数的更新方向,最终通过对速度场模型参数进行迭代更新,直至数据残差达到设定阈值或迭代更新次数达到预设值,停止迭代更新确定速度场模型参数,此时采用该参数的速度场模型进行全波形反演。因此本发明能对煤层底板下方条件复杂的陷落柱进行精确定位及形态探测,从而为后续开采提供数据支撑。
附图说明
图1是某煤矿工作面已知的陷落柱位置及形态示意图;
图2是采用本发明实施例对图1煤矿工作面进行纵波反演成像图;
图3是采用本发明实施例对图1煤矿工作面进行横波反演成像图;
图4是采用现有的纵波层析反演对图1煤矿工作面进行反演成像图;
图5是本发明实施例中L型地震观测系统的位置示意图;
其中,星号表示震源;三角形表示检波器。
具体实施方式
下面将结合本文实施例中的附图,对本文实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本文一部分实施例,而不是全部的实施例。基于本文中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本文保护的范围。
本发明实施例的反演方法,具体步骤为:
A、在煤层工作面一侧的巷道壁相互等间距布设多个震源;在煤层工作面另一侧的巷道壁相互等间距5m布设多个检波器,在煤层工作面的切眼处等间距5m设有一排检波器,各个检波器均通过连接主线与采集基站及采集主机连接,如图5所示,使检波器、采集基站及采集主机形成L型地震观测系统;激发多个震源,通过地震观测系统接收观测数据并获得各个地震子波;
B、先构建基于卷积型的目标函数,具体定义如下:
其中合成数据中的正传波场通过弹性波动方程合成,具体公式如下:
L[m(x)]u(x,t;m)=s(xs,t) (2)
其中,xs是震源坐标,L[]是弹性波正演模拟算子,u(x,t;m)是正传波场,s(xs,t)为震源;
然后将步骤A获得地震子波及震源坐标代入公式(2)计算得出震源的正传波场;
C、利用拉格朗日伴随方法获得模型参数的梯度,弹性波动方程的伴随方程为:
根据公式(4)计算得出基于卷积型目标函数对应的伴随源,然后将求取的伴随源代入公式(3)中进行数据模拟,最终获得伴随源的反转波场;
D、先确定速度场模型参数m=(vp,vs)的梯度表达式如下:
将步骤B计算出的正传波场和步骤C计算出的反传波场代入到公式(5)中,最终计算得到速度场模型参数的梯度;
E、采用共轭梯度方法求取更新方向Δmk+1,
其中下标k表示迭代次数,sk是用于计算共轭方向的一个中间变量,上标T表示矩阵转置,更新后的速度场模型为:
其中,tk表示选取的更新步长,tkΔmk+1即为此次迭代过程中速度场模型参数的更新量;
通过对速度场模型参数进行迭代更新,直至数据残差达到设定阈值或迭代更新次数达到预设值,停止迭代更新确定速度场模型参数,此时采用该参数的速度场模型进行全波形反演。
效果验证:
在某煤矿其工作面走向长度为500m,切眼长200m,如图1所示,通过现有的多种方法先确定其工作面内存在五个大小形态不一的陷落柱,陷落柱的位置及速度场模型参数如表1所示,五个陷落柱的尺寸差异较大,最小的陷落柱长轴约15m,短轴约10m,最大的陷落柱长轴约56m,短轴约26.6m。此外,陷落柱呈现不同的产状,长轴与水平方向的夹角在-37.5°到37°之间变化。由于位于同一工作面的陷落柱波速特征差异不大,确定五个陷落柱的纵波速度为2600m/s,横波速度为1501m/s,底板围岩的纵波速度为3200m/s,横波速度为1847m/s;
表1
然后在相同的观测数据和观测系统条件下,分别采用本发明实施例和现有的纵波层析反演方法对该工作面进行反演成像,其中本发明实施例获得的反演成像如图2和3所示,现有的纵波层析反演成像如图4所示;
将两者成像图分别与图1已知的现有位置及形态示意图进行比较可知,现有的层析成像结果分辨率明显低,并且无法准确地指示出尺寸陷落柱的形态以及物性参数,同时,对小尺度的陷落柱,不能有效地指示。而采用本发明实施例的方法,可以得到高分辨率的成像结果,并且显示的陷落柱位置及形态与实际情况吻合度较高,因此本发明具有较高的精度。
Claims (3)
1.一种适用于煤层底板复杂陷落柱的全波形反演方法,其特征在于,该方法的具体步骤为:
A、在煤层工作面一侧的巷道壁相互等间距布设多个震源;在煤层工作面另一侧的巷道壁相互等间距布设多个检波器,多个检波器均通过连接主线与采集基站及采集主机连接,使检波器、采集基站及采集主机形成地震观测系统;激发多个震源,通过地震观测系统接收观测数据并获得各个地震子波;
B、先构建基于卷积型的目标函数,具体定义如下:
其中合成数据中的正传波场通过弹性波动方程合成,具体公式如下:
L[m(x)]u(x,t;m)=s(xs,t) (2)
其中,xs是震源坐标,L[]是弹性波正演模拟算子,u(x,t;m)是正传波场,s(xs,t)为震源;
然后将步骤A中获得的地震子波及震源坐标代入公式(2),计算得出震源的正传波场;
C、利用拉格朗日伴随方法获得模型参数的梯度,弹性波动方程的伴随方程为:
根据公式(4)计算得出基于卷积型的目标函数对应的伴随源,然后将求取的伴随源代入公式(3)中进行数据模拟,最终获得伴随源的反传波场;
D、先确定速度场模型参数m=(vp,vs)的梯度表达式如下:
将步骤B计算出的正传波场和步骤C计算出的反传波场代入到公式(5)中,最终计算得到速度场模型参数的梯度;
E、采用共轭梯度方法求取更新方向,根据得到的梯度对速度场模型参数进行迭代更新,直至数据残差达到设定阈值或迭代更新次数达到预设值,停止迭代更新确定速度场模型参数,此时采用该速度场模型参数的速度场模型进行全波形反演。
3.根据权利要求1所述的一种适用于煤层底板复杂陷落柱的全波形反演方法,其特征在于,所述步骤A中在煤层工作面的切眼处等间距设有一排检波器,使其与煤层工作面另一侧的巷道壁内设置的检波器形成L型观测系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110650586.6A CN113376695B (zh) | 2021-06-11 | 2021-06-11 | 一种适用于煤层底板复杂陷落柱的全波形反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110650586.6A CN113376695B (zh) | 2021-06-11 | 2021-06-11 | 一种适用于煤层底板复杂陷落柱的全波形反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113376695A CN113376695A (zh) | 2021-09-10 |
CN113376695B true CN113376695B (zh) | 2022-07-05 |
Family
ID=77573787
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110650586.6A Active CN113376695B (zh) | 2021-06-11 | 2021-06-11 | 一种适用于煤层底板复杂陷落柱的全波形反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113376695B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114660659B (zh) * | 2022-03-30 | 2023-03-31 | 中国矿业大学 | 一种纵横波解耦的高精度微震定位方法 |
CN116626751B (zh) * | 2023-05-26 | 2024-03-19 | 北京中矿大地地球探测工程技术有限公司 | 基于多目标函数的黏弹性参数同步反演方法、装置和设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105445798A (zh) * | 2014-08-21 | 2016-03-30 | 中国石油化工股份有限公司 | 一种基于梯度处理的全波形反演方法和系统 |
CN106526674A (zh) * | 2016-11-14 | 2017-03-22 | 中国石油化工股份有限公司 | 一种三维全波形反演能量加权梯度预处理方法 |
CN107505654A (zh) * | 2017-06-23 | 2017-12-22 | 中国海洋大学 | 基于地震记录积分的全波形反演方法 |
CN108845351A (zh) * | 2018-06-26 | 2018-11-20 | 中国石油大学(华东) | 一种vsp地震资料转换波全波形反演方法 |
CN111766628A (zh) * | 2020-07-29 | 2020-10-13 | 浪潮云信息技术股份公司 | 一种预条件的时间域弹性介质多参数全波形反演方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9910174B2 (en) * | 2014-07-25 | 2018-03-06 | Seoul National University R&Db Foundation | Seismic imaging apparatus and method for performing iterative application of direct waveform inversion |
EP3168653B1 (en) * | 2015-11-05 | 2019-07-17 | CGG Services SAS | Device and method for full waveform inversion |
-
2021
- 2021-06-11 CN CN202110650586.6A patent/CN113376695B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105445798A (zh) * | 2014-08-21 | 2016-03-30 | 中国石油化工股份有限公司 | 一种基于梯度处理的全波形反演方法和系统 |
CN106526674A (zh) * | 2016-11-14 | 2017-03-22 | 中国石油化工股份有限公司 | 一种三维全波形反演能量加权梯度预处理方法 |
CN107505654A (zh) * | 2017-06-23 | 2017-12-22 | 中国海洋大学 | 基于地震记录积分的全波形反演方法 |
CN108845351A (zh) * | 2018-06-26 | 2018-11-20 | 中国石油大学(华东) | 一种vsp地震资料转换波全波形反演方法 |
CN111766628A (zh) * | 2020-07-29 | 2020-10-13 | 浪潮云信息技术股份公司 | 一种预条件的时间域弹性介质多参数全波形反演方法 |
Non-Patent Citations (2)
Title |
---|
Data-driven low-frequency signal recovery using deep learning predictions in full-waveform inversion;Jinwei Fang et al.;《GEOPHYSICS》;20200831;第85卷(第6期);第1-40页 * |
地震波形分类技术在地质异常体解释中的应用;程彦 等;《煤田地质与勘探》;20201231;第48卷(第6期);第87-92、102页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113376695A (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101930083B (zh) | 井间电磁波层析成像多孔对联合反演方法 | |
CN108680964A (zh) | 一种基于结构约束的归一化重磁电震联合反演方法 | |
CN103809216B (zh) | 一种电阻率数据与地震数据联合速度建场方法 | |
CN113376695B (zh) | 一种适用于煤层底板复杂陷落柱的全波形反演方法 | |
CN107132571A (zh) | 一种用于隧道超前地质预报的多源地震干涉法 | |
CN112485823B (zh) | 高效综合超前地质预报方法 | |
CN105549087A (zh) | 一种煤矿井下槽波地震勘探的走时和振幅联合反演方法 | |
CN106610503A (zh) | 一种煤矿采掘过程中的全方位槽波地震探测方法 | |
US20230384470A1 (en) | Method for three-dimensional velocity geological modeling with structures and velocities randomly arranged | |
CN109765606A (zh) | 一种基于反射槽波的回采工作面隐蔽断层性质的探测方法 | |
CN107490808A (zh) | 一种高可靠性地震勘探观测系统的建立方法 | |
CN109459787B (zh) | 基于地震槽波全波形反演的煤矿井下构造成像方法及系统 | |
CN111812707A (zh) | 一种基于微动和面波勘探的地铁沿线场地内孤石探测方法 | |
CN102877828A (zh) | 一种三维多井联合井地ct成像方法 | |
Hanson et al. | Advanced techniques in site characterization and mining hazard detection for the underground coal industry | |
CN107942372A (zh) | 井间及井地联合地震波ct成像方法及装置 | |
CN117192615A (zh) | 基于透射地震尾波的采煤工作面内隐伏地质构造探测方法 | |
CN107817524A (zh) | 三维地震层析成像的方法和装置 | |
CN111413747A (zh) | 低山地貌区隐伏塌陷的预测方法 | |
CN106772637A (zh) | 一种基于槽波超前探测tvsp扫描方法 | |
CN109991660A (zh) | 基于槽波地震和克里金插值的煤层厚度预测方法及装置 | |
Jiskani et al. | Fault orientation modeling of Sonda-Jherruck coalfield, Pakistan | |
CN110687591B (zh) | 基于先验数据的波形匹配确定煤层及围岩物性参数的方法 | |
Ren et al. | Pre-stack elastic reverse time migration in tunnels based on cylindrical coordinates | |
Wang et al. | A polarization migration velocity model building method for geological prediction ahead of the tunnel face |
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 |