CN117826237A - 一种预测隧道内盾构机工作诱发地震波传播的方法 - Google Patents
一种预测隧道内盾构机工作诱发地震波传播的方法 Download PDFInfo
- Publication number
- CN117826237A CN117826237A CN202410010823.6A CN202410010823A CN117826237A CN 117826237 A CN117826237 A CN 117826237A CN 202410010823 A CN202410010823 A CN 202410010823A CN 117826237 A CN117826237 A CN 117826237A
- Authority
- CN
- China
- Prior art keywords
- tunnel
- shield machine
- model
- propagation
- vibration
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000011435 rock Substances 0.000 claims abstract description 25
- 230000005641 tunneling Effects 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 230000007246 mechanism Effects 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 208000005647 Mumps Diseases 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 3
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 230000005404 monopole Effects 0.000 claims description 3
- 208000010805 mumps infectious disease Diseases 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 7
- 238000009412 basement excavation Methods 0.000 abstract description 2
- 238000005293 physical law Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 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/02—Generating seismic energy
- G01V1/159—Generating seismic energy using piezoelectric or magnetostrictive driving means
-
- 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/143—Generating seismic energy using mechanical driving means, e.g. motor driven shaft
-
- 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
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
本发明属于隧道挖掘技术领域,具体涉及一种预测隧道内盾构机工作诱发地震波传播的方法,所述方法包括以下步骤:S1:构建盾构机工作模型,根据隧道盾构机工作机理,简化盾构机工作模型;S2:将盾构机模型网格化,在隧道横截面上构建四边形网格,以此为基础,沿着隧道方向进行网格扫描;S3:模拟震源震动,盾构机处于工作状态时,盾构机不断敲击岩壁,产生地震波;S4:当隧道充满空气时,利用时间域声波方程作为隧道的控制方程;S5:使用弹性波方程的时间显示形式描述物理规律,且岩壁中地震波的传播以弹性波传播为主。本发明能够更贴近实际工作状态,同时计算速度更快,计算效果更加良好,简化了模型,充分保留了其工作特征和预测精度。
Description
技术领域
本发明属于隧道挖掘技术领域,具体涉及一种预测隧道内盾构机工作诱发地震波传播的方法。
背景技术
盾构机是挖掘隧道的利器,但其在前进过程中无法探知前方可能存在的地质灾害。地震波传播距离远,对非均质构造有很好的勘探效果。为此,利用盾构机在工作时的震动作为震源,在岩壁上置入检波器,进而预测盾构机前方的非均质构造(如溶洞等)是研究者们长期以来致力的工作方向。
然而,由于环境的复杂性和地震信号的弱振幅特点,加之实际测试中存在非常多的噪音,根本无法利用测得的地震信号分析非均质构造,甚至无法分析其规律。为此,首先需构建模型,预测上述实际情况下地震波的特征。与此同时,通过模拟获得地震波的传播特征,进而指导研究人员在实际资料中寻找到有效信号。然而,目前几乎没有方法可以有效模拟盾构机工作时诱发的地震波的传播规律,因此,建立一个合适的有限元模型,并基于其预测隧道内盾构机工作状态是十分必要的。
发明内容
本发明的目的是提供一种预测隧道内盾构机工作诱发地震波传播的方法,能够更贴近实际工作状态,同时计算速度更快,计算效果更加良好,简化了模型,充分保留了其工作特征和预测精度。
本发明采取的技术方案具体如下:
一种预测隧道内盾构机工作诱发地震波传播的方法,所述方法包括以下步骤:
S1:构建盾构机工作模型,根据隧道盾构机工作机理,简化盾构机工作模型;
S2:将盾构机模型网格化,在隧道横截面上构建四边形网格,以此为基础,沿着隧道方向进行网格扫描;
S3:模拟震源震动,盾构机处于工作状态时,盾构机不断敲击岩壁,产生地震波;
S4:当隧道充满空气时,利用时间域声波方程作为隧道的控制方程;
S5:使用弹性波方程的时间显示形式描述物理规律,且岩壁中地震波的传播以弹性波传播为主;
S6:采用压电的震动方式模拟盾构机震动;
S7:隧道模型中岩壁为圆柱形,且使用完美匹配层边界条件,应用于岩壁外侧;
S8:使用MUMPS方法作为直接求解器。
所述S1中,模型构建直径为D1,长度为L1的隧道区域;内径为D1,外直径为D2,长度为L2的岩石区域;直径为D1,长度为L3的盾构机区域,并将模型沿中心方向切分。
所述S3中,构建震动函数V(t):
V(t)=V0*sin(2*pi*f0*t)*rect1(t)
其中:V0为震动幅度大小,f0为震动的主频,t为时间,rect1(t)为截断函数。
所述S4中,隧道控制方程具体为:
其中,ρ是密度,pt是整体声压,ut是声音速度场扰动,c是声波速度,是声压对时间的导数,/>是声音速度场的散度,Qn是单极域源,I是单位矩阵,/>是声压乘以单位矩阵的散度,qd是偶极域源,/>·是散度算子。
所述S5中,弹性波方程具体为:
S=C∶ε-E·e
其中,v是速度场,ρ是密度,S是应力张量,ε是应变张量,C是刚度张量,Fv是体力,是应变张量对时间的导数,/>是速度场的梯度,/>是速度场梯度的转置,E是电场,e是耦合矩阵。
所述S6中,盾构机推力产生的震动频率低,传播距离远,其控制方程具体为:
D=∈0∈rE
其中,·D为电位移的散度,D是电位移,∈0为空气介电常数,∈r为相对介电常数。
所述S7中,将对称边界条件应用到模型中心截面,用于增强计算效率。
所述S8中,使用牛顿迭代法获得最终解,即应力张量S,应变张量ε。
本发明取得的技术效果为:
本发明的一种预测隧道内盾构机工作诱发地震波传播的方法针对传统方法中相关策略的失误,通过到现场观察实际工作状态,提出流体泵的推拉作用是产生低频地震信号的主要原因,可以利用压电效应模拟其震动效果,从而避免考虑刀盘产生的高频影响。
本发明的一种预测隧道内盾构机工作诱发地震波传播的方法在模型控制方程和边界条件设置方面,使用多场耦合作为控制方程,充分模拟了实际情况;同时,使用完美匹配层模拟岩壁的无限空间;使用对称边界条件,提高计算效率。目前没有传统方法能与上述做法契合。
附图说明
图1是本发明实施例的流程框图;
图2是本发明实施例隧道模型的结构示意图;
其中,A、隧道;B、岩壁;C、刀盘;D、匹配边界区域;
图3是本发明实施例隧道模型的网格剖分图;
图4是本发明实施例矩形截断函数示意图;
图5是本发明实施例震源震动模拟示意图;
图6是本发明实施例求解后三个检波点位置的示意图。
具体实施方式
为了使本发明的目的及优点更加清楚明白,以下结合实施例对本发明进行具体说明。应当理解,以下文字仅仅用以描述本发明的一种或几种具体的实施方式,并不对本发明具体请求的保护范围进行严格限定。
实施例1:
如图1-6所示,一种预测隧道内盾构机工作诱发地震波传播的方法,所述方法包括以下步骤:
S1:构建盾构机工作模型,根据隧道盾构机工作机理,简化盾构机工作模型;模型构建直径为D1,长度为L1的隧道区域;内径为D1,外直径为D2,长度为L2的岩石区域;直径为D1,长度为L3的盾构机区域,并将模型沿中心方向切分。
S2:将盾构机模型网格化,在隧道横截面上构建四边形网格,以此为基础,沿着隧道方向进行网格扫描;
S3:模拟震源震动,盾构机处于工作状态时,为保证盾构机能够向前移动,需使用流体泵不断向前推,盾构机不断敲击岩壁,产生地震波;由于盾构机上的刀片同样会产生一些地震波,但是其频率过高,噪音非常大;所以,勘探以向前推进的震动为主,构建震动函数V(t):
V(t)=V0*sin(2*pi*f0*t)*rect1(t)
其中:V0为震动幅度大小,f0为震动的主频,t为时间,rect1(t)为截断函数。
S4:当隧道充满空气时,利用时间域声波方程作为隧道的控制方程;隧道控制方程具体为:
其中,ρ是密度,pt是整体声压,ut是声音速度场扰动,c是声波速度,是声压对时间的导数,/>是声音速度场的散度,Qm是单极域源,I是单位矩阵,/>是声压乘以单位矩阵的散度,qd是偶极域源,/>·是散度算子。
S5:使用弹性波方程的时间显示形式描述物理规律,且岩壁中地震波的传播以弹性波传播为主;弹性波方程具体为:
S=C∶ε-E·e
其中,v是速度场,ρ是密度,S是应力张量,ε是应变张量,C是刚度张量,Fv是体力,是应变张量对时间的导数,/>是速度场的梯度,/>是速度场梯度的转置,E是电场,e是耦合矩阵。
S6:采用压电的震动方式模拟盾构机震动;盾构机推力产生的震动频率低,传播距离远,其控制方程具体为:
D=∈0∈rE
其中,·D为电位移的散度,D是电位移,∈0为空气介电常数,∈r为相对介电常数。
S7:隧道模型中岩壁为圆柱形,且使用完美匹配层边界条件,应用于岩壁外侧;将对称边界条件应用到模型中心截面,用于增强计算效率。
S8:使用MUMPS方法作为直接求解器,使用牛顿迭代法获得最终解,即应力张量S,应变张量ε。
实施例2:
如图2-图3所示,基于实施例1的基础上,本实施例对具体实施方式做出清楚完整的说明:
依据S1:构建直径为5米,长度为15米的隧道;内径为5米,外径为10米,长度为7.5米的岩石区域;直径为5米,长度为2.5米的盾构机区域,其中隧道为A区,岩石区域为B区,盾构机区域为C区,D区为拓展区域;
依据S2:构建隧道模型的网格,隧道的横截面采用四边形网格,然后沿着隧道方向扫描;
依据S3:模拟盾构机的震动。其中V0取50福特,频率f0为20Hz,rect1(t)为矩形截断函数,效果如图4所示。将震动函数与矩形函数相乘获得的最终震源函数,如图5所示;
依据S4:将控制方程加载到隧道区域,即图2中所示A区,同时也包括拓展区域(图2中所示D区);
依据S5:将岩壁控制方程加载到岩壁区域,即图2所示B区;
依据S6:将盾构机控制方程加载到盾构机区域,即图2所示C区;
依据S7:将完美匹配层条件加载到岩壁的外圈,即图2所示B区的外层部分;
依据S8:求解,结果如图4所示,可以看到,除了为考虑到的噪音外,本发明的预测结果与真正测量结果基本一致。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本发明中未具体描述和解释说明的结构、装置以及操作方法,如无特别说明和限定,均按照本领域的常规手段进行实施。
Claims (8)
1.一种预测隧道内盾构机工作诱发地震波传播的方法,其特征在于:所述方法包括以下步骤:
S1:构建盾构机工作模型,根据隧道盾构机工作机理,简化盾构机工作模型;
S2:将盾构机模型网格化,在隧道横截面上构建四边形网格,以此为基础,沿着隧道方向进行网格扫描;
S3:模拟震源震动,盾构机处于工作状态时,盾构机不断敲击岩壁,产生地震波;
S4:当隧道充满空气时,利用时间域声波方程作为隧道的控制方程;
S5:使用弹性波方程的时间显示形式描述物理规律,且岩壁中地震波的传播以弹性波传播为主;
S6:采用压电的震动方式模拟盾构机震动;
S7:隧道模型中岩壁为圆柱形,且使用完美匹配层边界条件,应用于岩壁外侧;
S8:使用MUMPS方法作为直接求解器。
2.根据权利要求1所述的一种预测隧道内盾构机工作诱发地震波传播的方法,其特征在于:所述S1中,模型构建直径为D1,长度为L1的隧道区域;内径为D1,外直径为D2,长度为L2的岩石区域;直径为D1,长度为L3的盾构机区域,并将模型沿中心方向切分。
3.根据权利要求1所述的一种预测隧道内盾构机工作诱发地震波传播的方法,其特征在于:所述S3中,构建震动函数V(t):
V(t)=V0*sin(2*pi*f0*t)*rect1(t)
其中:V0为震动幅度大小,f0为震动的主频,t为时间,rect1(t)为截断函数。
4.根据权利要求1所述的一种预测隧道内盾构机工作诱发地震波传播的方法,其特征在于:所述S4中,隧道控制方程具体为:
其中,ρ是密度,pt是整体声压,ut是声音速度场扰动,c是声波速度,是声压对时间的导数,/>是声音速度场的散度,Qm是单极域源,i是单位矩阵,/>是声压乘以单位矩阵的散度,qd是偶极域源,/>是散度算子。
5.根据权利要求1所述的一种预测隧道内盾构机工作诱发地震波传播的方法,其特征在于:所述S5中,弹性波方程具体为:
S=C∶ε-E·e
其中,v是速度场,ρ是密度,S是应力张量,ε是应变张量,C是刚度张量,Fv是体力,是应变张量对时间的导数,/>是速度场的梯度,/>是速度场梯度的转置,E是电场,e是耦合矩阵。
6.根据权利要求1所述的一种预测隧道内盾构机工作诱发地震波传播的方法,其特征在于:所述S6中,盾构机推力产生的震动频率低,传播距离远,其控制方程具体为:
D=∈0∈rE
其中,为电位移的散度,D是电位移,∈0为空气介电常数,∈r为相对介电常数。
7.根据权利要求1所述的一种预测隧道内盾构机工作诱发地震波传播的方法,其特征在于:所述S7中,将对称边界条件应用到模型中心截面,用于增强计算效率。
8.根据权利要求1所述的一种预测隧道内盾构机工作诱发地震波传播的方法,其特征在于:所述S8中,使用牛顿迭代法获得最终解,即应力张量S,应变张量ε。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410010823.6A CN117826237A (zh) | 2024-01-04 | 2024-01-04 | 一种预测隧道内盾构机工作诱发地震波传播的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410010823.6A CN117826237A (zh) | 2024-01-04 | 2024-01-04 | 一种预测隧道内盾构机工作诱发地震波传播的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117826237A true CN117826237A (zh) | 2024-04-05 |
Family
ID=90518908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410010823.6A Pending CN117826237A (zh) | 2024-01-04 | 2024-01-04 | 一种预测隧道内盾构机工作诱发地震波传播的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117826237A (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1826588A2 (de) * | 2006-02-28 | 2007-08-29 | GeoForschungsZentrum Potsdam | Verfahren und Vorrichtung zur Vorauserkundung beim Tunnelbau |
CN102681004A (zh) * | 2012-05-14 | 2012-09-19 | 中国矿业大学(北京) | 以掘进机为震源的巷道随掘地震超前探测装置及方法 |
JP2014106075A (ja) * | 2012-11-27 | 2014-06-09 | Fujita Corp | トンネル掘削中の地質探査方法 |
KR101547508B1 (ko) * | 2014-09-04 | 2015-08-26 | (주)희송지오텍 | 전방 지질 예측 시스템 및 그 방법 |
US20170218757A1 (en) * | 2015-03-11 | 2017-08-03 | Shandong University | Rock breaking seismic source and active source three-dimensional seismic combined advanced detection system using tunnel boring machine |
CN109188512A (zh) * | 2018-09-17 | 2019-01-11 | 中国石油大学(华东) | 基于非规则扇形网格剖分的起伏隧道空间正演模拟***及方法 |
CN114969884A (zh) * | 2022-03-21 | 2022-08-30 | 武汉大学 | 一种盾构隧道开挖过程及地表变形三维有限差分数值模拟方法 |
US20220326403A1 (en) * | 2020-05-12 | 2022-10-13 | Shandong University | Multi-wavefield seismic detection method and system based on construction noise of shield machine |
CN115561810A (zh) * | 2022-09-29 | 2023-01-03 | 安徽理工大学 | 一种矿井随掘地震数据正演模拟方法 |
-
2024
- 2024-01-04 CN CN202410010823.6A patent/CN117826237A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1826588A2 (de) * | 2006-02-28 | 2007-08-29 | GeoForschungsZentrum Potsdam | Verfahren und Vorrichtung zur Vorauserkundung beim Tunnelbau |
CN102681004A (zh) * | 2012-05-14 | 2012-09-19 | 中国矿业大学(北京) | 以掘进机为震源的巷道随掘地震超前探测装置及方法 |
JP2014106075A (ja) * | 2012-11-27 | 2014-06-09 | Fujita Corp | トンネル掘削中の地質探査方法 |
KR101547508B1 (ko) * | 2014-09-04 | 2015-08-26 | (주)희송지오텍 | 전방 지질 예측 시스템 및 그 방법 |
US20170218757A1 (en) * | 2015-03-11 | 2017-08-03 | Shandong University | Rock breaking seismic source and active source three-dimensional seismic combined advanced detection system using tunnel boring machine |
CN109188512A (zh) * | 2018-09-17 | 2019-01-11 | 中国石油大学(华东) | 基于非规则扇形网格剖分的起伏隧道空间正演模拟***及方法 |
US20220326403A1 (en) * | 2020-05-12 | 2022-10-13 | Shandong University | Multi-wavefield seismic detection method and system based on construction noise of shield machine |
CN114969884A (zh) * | 2022-03-21 | 2022-08-30 | 武汉大学 | 一种盾构隧道开挖过程及地表变形三维有限差分数值模拟方法 |
CN115561810A (zh) * | 2022-09-29 | 2023-01-03 | 安徽理工大学 | 一种矿井随掘地震数据正演模拟方法 |
Non-Patent Citations (5)
Title |
---|
ANDRE LAMERT等: "Imaging disturbance zones ahead of a tunnel by elastic full-waveform inversion: Adjoint gradient based inversion vs. parameter space reduction using a level-set method", 《UNDERGROUND SPACE》, vol. 3, no. 1, 31 March 2018 (2018-03-31), pages 21 - 33 * |
刘绍宝等: "TBM施工地下水预报技术初探", 《工程地球物理学报》, vol. 6, no. 1, 31 December 2009 (2009-12-31), pages 117 - 120 * |
罗帅等: "基于柱坐标系的隧道空间溶洞探测全波场数值模拟", 《地球物理学进展》, vol. 35, no. 03, 30 June 2020 (2020-06-30), pages 1175 - 1182 * |
郭立全等: "巷道构造震波超前探测技术及其应用", 《煤炭科学技术》, vol. 36, no. 11, 30 November 2008 (2008-11-30), pages 99 - 101 * |
鲁光银等: "隧道反射波超前探测有限差分正演模拟与偏移处理", 《中南大学学报(自然科学版)》, vol. 42, no. 01, 31 January 2011 (2011-01-31), pages 136 - 141 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110109895B (zh) | 适用于tbm掘进隧道的围岩分级联合预测方法及应用 | |
Sharpe | The production of elastic waves by explosion pressures. I. Theory and empirical field observations | |
Cai et al. | FLAC/PFC coupled numerical simulation of AE in large-scale underground excavations | |
CN108957521B (zh) | 一种用于隧道长距离三维超前地质预报方法 | |
WO2017210231A1 (en) | Estimation of horizontal stresses and nonlinear constants in anisotropic formations such as interbedded carbonate layers in organic-shale reservoirs | |
CN100456035C (zh) | 基于***源的大坝无损检测*** | |
GB2472300A (en) | Using dipole compressional data obtained in a borehole to determined properties of a subterranean formation | |
Buffi et al. | Calibration of finite element models of concrete arch-gravity dams using dynamical measures: the case of Ridracoli | |
Gortsas et al. | Numerical modelling of micro-seismic and infrasound noise radiated by a wind turbine | |
Panji et al. | Seismic analysis of semi-sine shaped alluvial hills above subsurface circular cavity | |
Jiao et al. | Nonlinear seismic response of rock tunnels crossing inactive fault under obliquely incident seismic P waves | |
Gao et al. | Strength reduction model for jointed rock masses and peridynamics simulation of uniaxial compression testing | |
CN117826237A (zh) | 一种预测隧道内盾构机工作诱发地震波传播的方法 | |
CN117471532A (zh) | 一种基于***地震波速的深部岩体的地应力反演方法 | |
Zhuang et al. | The rotated Cartesian coordinate method to remove the axial singularity of cylindrical coordinates in finite‐difference schemes for elastic and viscoelastic waves | |
CN116227309B (zh) | 融合多源数据的盾构掘进数字孪生地层构建方法及*** | |
CHE et al. | Acoustic Field in Formation Generated by Linear Phased Array Transmitters in Fluid‐Filled Boreholes | |
Xu | Blasting vibration safety criterion of surrounding rock of a circular tunnel | |
Wu et al. | Research on dynamic response of shallow buried tunnel lining constructed by drilling and blasting method | |
Wang et al. | Analysis of blasting vibration effect of railway tunnel and determination of reasonable burial depth | |
Johari et al. | System reliability analysis of site slope using the conditional spatial variability of soil properties | |
JP3410965B2 (ja) | 地質探査方法、弾性波発生方法、及び弾性波発生装置 | |
Chen et al. | Advance grouting test using seismic prospecting method in TBM tunneling: A case study in the water supply project from Songhua River, Jilin, China | |
CN207144934U (zh) | 一种过钻头偶极子声波测井发射器及测井装置 | |
Yue et al. | Acoustic field simulations of logging while drilling by cylindrical finite difference with variable grids |
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 |