CN101943599A - 一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法 - Google Patents

一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法 Download PDF

Info

Publication number
CN101943599A
CN101943599A CN 201010290417 CN201010290417A CN101943599A CN 101943599 A CN101943599 A CN 101943599A CN 201010290417 CN201010290417 CN 201010290417 CN 201010290417 A CN201010290417 A CN 201010290417A CN 101943599 A CN101943599 A CN 101943599A
Authority
CN
China
Prior art keywords
wave
measuring
velocity
tunnel
face
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
CN 201010290417
Other languages
English (en)
Other versions
CN101943599B (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN 201010290417 priority Critical patent/CN101943599B/zh
Publication of CN101943599A publication Critical patent/CN101943599A/zh
Application granted granted Critical
Publication of CN101943599B publication Critical patent/CN101943599B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法,是在掌子面上设一条水平和一条铅垂向测线,测线上每隔20~40cm设一测点,在掌子面后方边墙设铅垂向短测线,测线上设4~8个测点,采用5~40cm震-检距测出各测点的时间曲线,然后通过归一化将一条测线各测点时间曲线汇成时间剖面,在时间剖面上得到测点到反射面反射时间。边墙测量是在激震点钻孔、孔中插入短钢筋,锤击短钢筋激震,分别测取纵波或横波。设定多个虚拟波速,计算边墙测点在虚拟波速情况下对虚拟反射面的反射时间与实测反射时间作对比,求出各测点的波速差值,绘出差值-波速曲线,曲线上差值为零的点对应的波速即为实际波速。

Description

一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法
技术领域
一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法,本发明涉及在隧道和其它地下工程掌子面探查前方地质时测取前方岩体波速的方法,属工程地质勘查技术领域。
背景技术
目前投入使用的探测隧道掌子面前方地质的方法中,属于反射类的方法有地震反射负视速度法(隧道VSP)、瑞士的TSP、日本的HSP、陆地声纳等弹性波反射类方法,以及探地雷达等电磁类方法。它们均是通过测量地质体的反射面的反射时间t和岩体波速V求取反射体空间位置:
距离L=Vt/2
其中L-距离(米),
V-岩体波速(米/秒),
t-反射时间(秒)。
但是,岩体波速测定,目前仍存在难题。以弹性波反射类方法为例,都是通过在掌子面后面边墙或隧道底部等其它部位设测线,用一点激震、多点接收或一点接收、多点激震的方式测得直达波或岩体内部松弛带与压密区等界面的折射波,测得第一层岩体波速。当反射面与掌子面平行时,可在隧道侧墙或隧道底部等设测线,在距掌子面远端设激震点,近端设若干检波点,用一点激震,多点接收方式,可以测出第二层,第三层等岩体的波速。但当反射面不平行掌子面(反射面不垂直测线)时,因为反射面的产状不知道,而是要根据波速及其它参数求取的,故尚不能实测岩体波速。专利申请2010102666019给出了一种知道波速后根据掌子面上设的一条水平测线、一条测铅垂向测线,根据测线上测点测得的反射波反射时间求取地质体反射面空间位置的方法,但该专利申请未涉及如何测取波速的技术。
发明内容
本发明的目的是为克服上述缺点而提供一种测取掌子面前方岩体弹性波速的方法。
本发明的技术方案是:在掌子面上设定一条水平和一条铅垂且相互交叉的测线,在测线上每隔20~40cm间距设一测点,采用5~40cm震-检距用锤击激震测出各测点时间曲线,然后将一条测线各测点的时间曲线汇成时间剖面,在时间剖面上得到测点对反射面的反射时间。在掌子面后方边墙或隧道底或拱部等部位设一测线,在测线上设4~8个测点,点间距为20~40cm,震-检距为5~40cm,测得时间曲线并汇成时间剖面,在时间剖面上确定与掌子面得到的时间剖面上相同的反射面,并得到测点对反射面的反射时间。
在测量纵波波速时,掌子面上的测点用纵波检波器测量沿隧道纵向运动质点的振动,在测量横波波速时用横波检波器测量垂直隧道纵向运动的质点的振动。在边墙或隧道底部、拱部等部位的测点,则要向掌子面方向锤击激震,并用定向的检波器接收来自掌子面方向的反射波,测得各测点对反射面的实测反射时间,得到其时间曲线。在边墙或隧道底部、拱部施测的方法是:在激震点处钻深10~20cm的小孔,孔径略大于所用钢筋的的直径,孔中插入钢筋,向掌子面方向锤击钢筋激震。
用假设岩体波速的方法,根据掌子面上的测点测得的各测点对反射面的反射时间,按专利申请2010102666019提供的方法,求出虚拟的反射面的空间位置,然后计算边墙测点在假设的波速情况下对虚拟的反射面的反射时间,与实测反射时间比较,得到实测反射时间与虚拟反射时间的差值。假设多个波速,得到多个实测反射时间与虚拟反射时间的差值,绘出差值-波速曲线,则差值为零的点即为实际波速。或将反射时间的差值取绝对值后绘|差值|-波速曲线,曲线极小值点对应的波速即为实际波速。
附图说明:
图1为在隧道掌子面上布置一条水平测线和一条铅垂向的测线示意图。
图中:M-隧道掌子面,A-水平测线,B-铅垂向的测线,
1、2…14和I、II…VII-测点号。
图2为在隧道边墙布置的测点及其向掌子面前方反射面的入射、反射射线示意图。
图中:M-隧道掌子面,b-隧道边墙,(1)、(2)、(3)、(4)-测点号,
S-入射及反射射线,S′-入射及反射射线,
1-反射面1,面2-反射面2。
图3为一张实测的时间剖面图。
图中:①、②、③、④、⑤-分别为1~5反射面反射波同相轴。
图4为用以确定V0的|Δt|-V曲线示意图。
图中:V-所假设的波速值,|Δt|-反射时间差值的绝对值,
V0-所求出的实际波速值
min|Δt|-|Δt|的极小点,它所对应的V值,即为V0
图5为用以确定V0的Δt-V曲线示意图。
图中:V-所假设的波速值,Δt-反射时间差值,
V0-所求出的实际波速值
具体实施方式:
(1)在隧道掌子面上设一条水平测线A、一条铅垂向的测线B,在测线上每隔20~40cm设一测点,在每一测点按5~40cm震-检距用锤击锤和激震杆激震,用时在接收点接收,然后通过计算机作归一化处理将各测点的时间曲线汇成时间剖面,根据同相轴确定各反射面;
(2)在隧道边墙或隧底、拱部等部位的4~8个测点上,用5~40cm震-检距的方式,激震得到某反射面P到这些测点(1)、(2)、(3)…(g)的反射时间tP1、tP2、tP3…(P=1、2、3、4…,下同);
(3)在边墙或隧底或拱部的测量方法:
在已设每一个锤击点处打一深10-20cm的小孔,其直径略大于所用钢筋的直径,孔中插入短钢筋,钢筋插入后外露15cm左右,锤击钢筋背向掌子面的侧面,采用2分量或3分量检波器,若测纵波速度,则用朝向隧道纵向的检波器,测取平行隧道纵向运动的质点振动的纵波;若测量横波波速,则采用2分量或3分量检波器测取掌子面前方反射的横波。
(4)在A、B测线的时间剖面上,根据同相轴确定各反射面,并选定3个或3个以上不在同一测线上的测点1、2、3、…n(n≥3)和I、II…m(m≥3);
假设第P层岩体均方波速VσP设,根据VσP设,计算出第P个反射面面设P的空间位置;
(5)根据面设P的位置,计算点(1)、(2)、(3)…的反射时间tP1设、tP2设、…并计算出ΔtP1、ΔtP2、…ΔtPg
ΔtPg=tPg-tPg设
若ΔtP1、ΔtP2…ΔtPg=0,
则所选定的波速Vσ设P等于实际波速VσP
若ΔtP1、ΔtP2…ΔtPg>0,
则所选定的波速VσP设小于实际波速VσP
若ΔtP1、ΔtP2…<0,
则所选定的波速VσP设大于实际波速VσP
假设多个波速,得到多个实测反射时间与虚拟反射时间的差值,绘出|差值|-波速曲线(|差值|即差值的绝对值),则|差值|为零的点或曲线极小值点对应的波速即为实际波速。或绘出差值-波速曲线,差值为零的点对应的波速即为实际波速。
作为一个实例,对本发明进一步说明如下:
(1)在隧道掌子面上设一条水平测线A,一条铅垂向测线B,每隔30cm设一测点,即1、2、3、4、5…(A线),I、III、…(B线),采用震-检为30cm,测出各测点时间曲线,通过计算机归一化计算汇成时间剖面(图3),从时间剖面上追踪出第一个反射面面1,第二个反射面面2,第三个反射面3…的同相轴;
(2)在掌子面后方边墙设一铅垂方向测线,以30cm点距设4个测点,以20-40cm的震-检距测量,并由此4个测点组成的测线的时间剖面分辨出面1,面2,面3…的同相轴,读出面1的反射时间t11、t12、t13、t14,面2的反射时间t21、t22、t23、t24,和面3、面4的反射时间。
(3)计算各层岩体的均方根波速VσP,以第2层均方根波速Vσ2为例:
设定一组V2,V21设=2000m/s、V22设=2050m/s、V23设=2100m/s…V2N设=V22N-1设+50m/s,
计算出4个边墙测点数组虚拟的反射时间:
t121设、t122设、t123设
t221设、t222设、t223设
t321设、t322设、t323设
t421设、t422设、t423设
再计算出各点计算的反射时间与实测反射时间之差
Δt121、Δt122、Δt123
Δt221、Δt222、Δt223
Δt321、Δt322、Δt323
Δt421、Δt422、Δt423
其中Δt121=t21-t121设、Δt122=t21-t122设、Δt123=t21-t123设、…Δt12N=t21-t12N设
Δt221=t22-t221设、Δt222=t22-t222设、Δt223=t22-t223设、…Δt22N=t22-t22N设
Δt321=t23-t321设、Δt322=t23-t322设、Δt323=t23-t323设、…Δt32N=t23-t32N设
Δt421=t24-t421设、Δt422=t24-t422设、Δt423=t24-t423设、…Δt42N=t24-t42N设
然后用以上4组数据绘出4条|Δt|-V曲线(图4),其纵坐标为|Δt|,横坐标为Vσ;从曲线上找到极小点,即min|Δt|的点。min|Δt|的点对应的|Δt|应等于0或接近于0,它所对应的V值,即为Vσ2,取上述4个测点所得的Vσ2的平均值,即得第2层岩体均方根波速。
照此,可计算出各层岩体的均方根波速。
均方根波速与各层岩体的波速符合下式
V 0 = ( Σ i = 1 n t i V i 2 / Σ i = 1 n t i ) 1 / 2
其中ti为相应反射面的反射时间,
Vi为相应层的岩体波速,
可据此式计算各层的岩体波速。

Claims (2)

1.一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法,其特征在于:在掌子面上设定一条水平、一条铅垂且相互交叉的测线,在测线上每隔20~40cm设一测点,采用5~40cm震-检距测出各测点时间曲线,然后将一条测线各点时间曲线汇成时间剖面,在时间剖面上得到测点对反射面的反射时间。在掌子面后方边墙设铅垂向测线,上设4-8个测点,点间距为用20~40cm,震-检距为5~40cm,测得时间曲线并汇成时间剖面,在时间剖面上确定于掌子面得到的时间剖面上相同的反射面,并得到测点对反射面的反射时间;用锤击锤及激震杆在掌子面激震,测量纵波速度时,采用测量质点沿隧道纵向运动的纵波检波器,测量横波波速时,采用测量质点垂直隧道纵向运动的横波检波器;在边墙激震点处打一深10~20cm的小孔,孔中插入钢筋,钢筋插入后外露15cm左右,锤击钢筋背向掌子面的侧面激震,分别测量质点沿隧道纵向运动的纵波和垂直隧道纵向运动的横波;采用先假设岩体波速的方法,根据掌子面上的测点测得的对反射面的反射时间求出虚拟的反射面的空间位置,然后计算边墙测点在假设的波速情况下对反射面的反射时间,与实测反射时间比较,假设多个波速值,分别求出与实测值差值,最后选定真波速值。
2.根据权利要求1所述的在隧道中测取掌子面前方岩体波速的方法,其特征在于:先假设某反射面P到掌子面岩体的均方根波速VσP设,由掌子面上的测线上不少于3个不在同一测线上的相应测点的反射时,计算出一个虚拟的反射面设P,然后计算出边墙各测点对面设P的相应反射时tP1设、tP2设、tP3设…,与它们各自的实测反射时相比较,得到其差值Δ,设定多个波速,得到多个差值,然后绘出|差值|-波速曲线,则其极小值即为实测波速值;或绘出差值-波速曲线,其差值为零的点对应的波速即为实测波速值;所得到的是层面到掌子面的岩体的均方根速度。
CN 201010290417 2010-09-25 2010-09-25 一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法 Expired - Fee Related CN101943599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010290417 CN101943599B (zh) 2010-09-25 2010-09-25 一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010290417 CN101943599B (zh) 2010-09-25 2010-09-25 一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法

Publications (2)

Publication Number Publication Date
CN101943599A true CN101943599A (zh) 2011-01-12
CN101943599B CN101943599B (zh) 2013-08-21

Family

ID=43435713

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010290417 Expired - Fee Related CN101943599B (zh) 2010-09-25 2010-09-25 一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法

Country Status (1)

Country Link
CN (1) CN101943599B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103697999A (zh) * 2013-12-30 2014-04-02 中国科学院武汉岩土力学研究所 一种高应力硬岩tbm施工隧道微震波速实时获取方法
CN110174694A (zh) * 2019-02-26 2019-08-27 成都畅达通检测技术股份有限公司 一种超前地质预报数据采集及分析方法
CN112130204A (zh) * 2020-09-02 2020-12-25 武汉大学 一种适用于深埋隧洞围岩波速的现场测试方法
CN113267806A (zh) * 2021-05-28 2021-08-17 长江水利委员会长江科学院 Tbm刀盘破岩噪声源的多波采集系统及超前探测方法
CN113311480A (zh) * 2021-05-28 2021-08-27 山东大学 一种基于高频弹性波的隧道不良地质超前探测系统及方法
CN114674418A (zh) * 2022-02-18 2022-06-28 深圳大学 隧道竖井爆破施工围岩内部质点振动的监测方法及装置
CN117233840A (zh) * 2023-11-15 2023-12-15 中国电建集团西北勘测设计研究院有限公司 一种基于地震波的地下洞室全空间超前地质预报方法
CN117233850A (zh) * 2023-10-19 2023-12-15 中国地质调查局成都地质调查中心(西南地质科技创新中心) 一种大地电磁信号处理方法与系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1967177A (zh) * 2005-11-16 2007-05-23 钟世航 一种在隧道中测取掌子面前方岩体波速的方法
JP2008267843A (ja) * 2007-04-17 2008-11-06 Tobishima Corp トンネル切羽面の測量システム
CN101672929A (zh) * 2009-09-23 2010-03-17 山东大学 在隧道中探查掌子面前方含水层并估算其出水量的方法
CN101839994A (zh) * 2010-05-07 2010-09-22 钟世航 锤击激震杆

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1967177A (zh) * 2005-11-16 2007-05-23 钟世航 一种在隧道中测取掌子面前方岩体波速的方法
JP2008267843A (ja) * 2007-04-17 2008-11-06 Tobishima Corp トンネル切羽面の測量システム
CN101672929A (zh) * 2009-09-23 2010-03-17 山东大学 在隧道中探查掌子面前方含水层并估算其出水量的方法
CN101839994A (zh) * 2010-05-07 2010-09-22 钟世航 锤击激震杆

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《物探与化探》 19970630 钟世航 陆地声纳法及其应用效果 172-179 1-2 第21卷, 第3期 2 *
《隧道建设》 20070831 钟世航等 隧道掌子面前方地质预报的进展 7-11 1-2 第27卷, 第4期 2 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103697999A (zh) * 2013-12-30 2014-04-02 中国科学院武汉岩土力学研究所 一种高应力硬岩tbm施工隧道微震波速实时获取方法
CN110174694A (zh) * 2019-02-26 2019-08-27 成都畅达通检测技术股份有限公司 一种超前地质预报数据采集及分析方法
CN112130204A (zh) * 2020-09-02 2020-12-25 武汉大学 一种适用于深埋隧洞围岩波速的现场测试方法
CN113267806A (zh) * 2021-05-28 2021-08-17 长江水利委员会长江科学院 Tbm刀盘破岩噪声源的多波采集系统及超前探测方法
CN113311480A (zh) * 2021-05-28 2021-08-27 山东大学 一种基于高频弹性波的隧道不良地质超前探测系统及方法
CN114674418A (zh) * 2022-02-18 2022-06-28 深圳大学 隧道竖井爆破施工围岩内部质点振动的监测方法及装置
CN114674418B (zh) * 2022-02-18 2022-11-22 深圳大学 隧道竖井爆破施工围岩内部质点振动的监测方法及装置
CN117233850A (zh) * 2023-10-19 2023-12-15 中国地质调查局成都地质调查中心(西南地质科技创新中心) 一种大地电磁信号处理方法与系统
CN117233850B (zh) * 2023-10-19 2024-01-26 中国地质调查局成都地质调查中心(西南地质科技创新中心) 一种大地电磁信号处理方法与系统
CN117233840A (zh) * 2023-11-15 2023-12-15 中国电建集团西北勘测设计研究院有限公司 一种基于地震波的地下洞室全空间超前地质预报方法
CN117233840B (zh) * 2023-11-15 2024-02-20 中国电建集团西北勘测设计研究院有限公司 一种基于地震波的地下洞室全空间超前地质预报方法

Also Published As

Publication number Publication date
CN101943599B (zh) 2013-08-21

Similar Documents

Publication Publication Date Title
CN101943599B (zh) 一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法
Bakulin et al. Smart DAS upholes for simultaneous land near-surface characterization and subsurface imaging
DONG et al. Three-dimensional analytical solution of acoustic emission or microseismic source location under cube monitoring network
CN101975820B (zh) 海底沉积声学参数原位测量装置
CN103534615A (zh) 使用分布式声传感器的地下监控
JP2007231729A (ja) トンネル建設の際に事前探査するための方法および装置
WO2012044481A2 (en) Method for offset timing of simultaneous seismic source firing
CN111856551A (zh) 浅层横向高分辨率瑞雷波勘探方法及系统
CN108802193A (zh) 一种巷道围岩松动圈的探测设备及探测方法
CN106443776A (zh) 一种基于时间切片法的海底地震仪重定位方法
KR20120076952A (ko) 굴절법 해양탄성파 탐사를 위한 obc 타입 스트리머 장치
CN104122582A (zh) 利用叠加速度求取高精度地震波速度的方法
CN112130204B (zh) 一种适用于深埋隧洞围岩波速的现场测试方法
Živor et al. Measurement of P-and S-wave velocities in a rock massif and its use in estimation elastic moduli
US9423521B2 (en) Method of determining the position of a detector disposed at the bottom of the sea
CN102087368B (zh) 一种近于零震-检距的反射法测取波速的方法
US20120287757A1 (en) Forward Looking Seismics From Drill-Bit
JP5186538B2 (ja) 地山探査方法
GB2442244A (en) Determining the position and orientation of electromagnetic receivers
CN104502913B (zh) 揭煤巷道与煤层距离的测量方法及装置
RU2046376C1 (ru) Способ оценки напряженного состояния горных пород
CN110780345A (zh) 一种隧道超前地震探测地震数据的三维速度分析方法
Du et al. Coal mine geological hazardous body detection using surface ground penetrating radar velocity tomography
JP2004138447A (ja) 岩盤の物性評価方法
CN104406681B (zh) 一种实时确定微震波速的测试方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB02 Change of applicant information

Address after: 100081 Beijing city Haidian District Xizhimen Daliushu Road No. 17 Fuhai Center Building 5, room 1603

Applicant after: Zhong Shihang

Applicant after: Wang Rong

Address before: 100081 Beijing city Haidian District Daliushu Road No. 2 Building 1 North Railway Science Research Institute 5121

Applicant before: Zhong Shihang

Applicant before: Wang Rong

C14 Grant of patent or utility model
GR01 Patent grant
C56 Change in the name or address of the patentee
CP02 Change in the address of a patent holder

Address after: 100044 Beijing city Haidian District West Park Village Hospital No. 3 of Changyuan No. 6 Building 1408 room

Patentee after: Zhong Shihang

Patentee after: Wang Rong

Address before: No. 5 building 100081 Beijing city Haidian District Xizhimen Daliushu Road No. 17 Fuhai center room 1603

Patentee before: Zhong Shihang

Patentee before: Wang Rong

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130821

CF01 Termination of patent right due to non-payment of annual fee