CN101950030B - 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法 - Google Patents

一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法 Download PDF

Info

Publication number
CN101950030B
CN101950030B CN2010102666019A CN201010266601A CN101950030B CN 101950030 B CN101950030 B CN 101950030B CN 2010102666019 A CN2010102666019 A CN 2010102666019A CN 201010266601 A CN201010266601 A CN 201010266601A CN 101950030 B CN101950030 B CN 101950030B
Authority
CN
China
Prior art keywords
line
point
reflecting surface
angle
survey
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
Application number
CN2010102666019A
Other languages
English (en)
Other versions
CN101950030A (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 CN2010102666019A priority Critical patent/CN101950030B/zh
Publication of CN101950030A publication Critical patent/CN101950030A/zh
Application granted granted Critical
Publication of CN101950030B publication Critical patent/CN101950030B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法。是在掌子面上设一条水平测线AB和一条与其交叉的铅垂测线CD,在测线上每隔20~40cm设一测点,采用20~40cm的震—检距激发弹性波并测出每个测点的时间曲线,然后将一条测线各点时间曲线汇成时间剖面,在时间剖面上得到测点到反射面的反射时间。在知道波速的情况下,计算出反射面的产状并确定他在在隧道掌子面前方的空间位置。计算辅助角α′=sin-1〔(BB′-AA′)/AB〕;α′比反射面的走向角α略大,反射面倾角γ在70°~90°和∠α在0°~30°的情况下,∠α≈∠α′。反射面倾角γ=90°-β,计算辅助角β=sin-1〔(DD′-CC′)/CD)〕。线段AAl≈AA′/cosβcosα′,线段BBl≈BB′/cosβcosα′。

Description

一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法
技术领域
一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法。本发明涉及在隧道或其他地下工程掌子面上用弹性波反射类方法探查前方岩体中断层等反射面的测量和计算方法,属隧道或地下工程施工地质预报技术领域。
背景技术
在隧道或其他地下工程中探查掌子面前方地质体时确定其产状和空间位置是一个重要内容,目前所采用的多种弹性波反射类方法、探地雷达等,都为全面而准确地确定这些地质体的空间位置而努力,突破的一种途径是力图用三维空间布置采集测点来解决。但非三维空间布置的采集测点的数据采集方法则遇到了困难。TSP、HSP等方法在隧道边墙布置排列的弹性波反射类方法仅能粗略地给出反射体位置,无法作出空间定位和给出断层等的产状,因而无法准确确定地质体反射面的空间位置。
发明内容:
本发明的目的是提供一种隧道施工地质预报时确定地质体位置的方法,即采用弹性波反射类方法测量、计算出掌子面前方地质体的反射面的产状,并确定其空间位置的方法。对溶洞等物体确定其距掌子面最近点的距离和位置,没有产状问题,较为简单。而对于断层等面状反射面,需要确定其走向、倾向、倾角的产状,并根据他们确定其空间位置,比较复杂。其做法是:
在掌子面上设定一条水平测线AB和一条与其相交叉的铅垂向测线CD,在测线上以20~40cm的距离设置测点。以20~40cm的震—检距激发弹性波并测出每个测点反射波的时间曲线,然后将一条测线上各点的时间曲线汇成时间剖面,在时间剖面上得到测点到反射面的反射时间。在知道波速的情况下,计算出地质体反射面的产状(走向、倾向、倾角)并确定它在隧道掌子面前方的空间位置。
结合附图来说明计算方法。
附图说明:
图1为在隧道掌子面上布置一条水平测线和一条铅垂向的测线示意图。
图中M—隧道掌子面;AB—水平测线;CD—铅垂向测线;
1、2…14—水平向测线AB上的测点号;
Ⅰ、Ⅱ…Ⅶ—铅垂向测线CD上的测点号;
图2为测线、射线及计算辅助线示意图。
图中M—隧道掌子面;AB—水平测线;CD—铅垂向测线;
N—反射面(断层面);
AA′⊥N、AA′为A点的入射波射线,A′为反射点;
BB′⊥N、BB′为B点的入射波射线,B′为反射点;
CC′⊥N、CC′为C点的入射波射线,C′为反射点;
DD′⊥N、DD′为D点的入射波射线,D′为反射点;
AAlBlB面为通过测线AB的水平面,垂直掌子面M,与反射面N交于AlBl线;
Al为通过测线AB的水平面上A点到反射面N并垂直掌子面的线AAl与N面的交点;
Bl为通过测线AB的水平面上B点到反射面N并垂直掌子面的线BBl与N面的交点;
CC′D′D面⊥M面并通过铅垂向测线CD,与面N交于线C′D′;
AAl⊥M面;
BBl⊥M面;
Ap点是通过A′点并与测线AB平行的线A′Ap与BB′线的交点;
Bp点是通过Al点并与测线AB平行的线AlBp与BBl线的交点;
D″是在DD′上、∠D′C′D″(∠β)的一边与线DD′的交点,C′D″⊥DD′;
面BB′B″⊥N面并与AlBl交于B″点;
BB″为面BB′B″和通过AB的水平面(AAlBlB面)与N面的交线;
A点的入射波射线AA′=V×TAA′/2,
B点的入射波射线BB′=V×TBB′/2,
C点的入射波射线CC′=V×TCC′/2,
D点的入射波射线DD′=V×TDD′/2;
式中:V—波速,
TAA′—测点A到反射面的反射时间,
TBB′—测点B到反射面的反射时间,
TCC′—测点C到反射面的反射时间,
TDD′—测点D到反射面的反射时间。
A l B p / / = AB / / = A ′ A p ;
∠B′A′Ap=α′,为平面ABB′A′上的一个角,为计算辅助角;
∠BlAlBp=α,为反射面N的走向角;
C′D″⊥DD′;
∠D′C′D″=β;β为C′D′与C′D″线的夹角;
γ=90°-β,为反射面N的倾角;
由于
Figure BSA00000248387600033
面,
sinβ=(DD′-CC′)/C′D″;
β=sin-1[(DD′-CC′)/C′D″];
由于
Figure BSA00000248387600041
以及A′A∥BB′,
所以sinα′=(BB′-AA′)/A′Ap=(BB′-AA′)/AB,
α′=sin-1[(BB′-AA′)/AB];
sinα=BlBp/AlBl
α=sin-1(BlBp/AlBl);
由于面A′B′Ap与面AAlBlB不平行(相互夹角β),
所以∠α′>∠α,
在工程实践中,当∠β在0°~20°和∠α在0°~30°的情况下,在计算反射面产状时可以接受∠α≈∠α′;
由于面BB′B″⊥面N∥面CC′D′D,
B′B″为面BB′B″和N面的交线,
所以B′B″=BB′/cosβ;
而面BB′B″与线BBl夹角为α,
即∠BlBB″=α,
又由于BBp⊥AlBp(∠α的一条边),
所以BB″也要垂直∠α的另一条边,
即BB″⊥AlBl
sinα′=(BB′-AA′)/AB,
α′=sin-1[(BB′-AA′)/AB],
所以BBl=BB″/cosα
       ≈BB″/cosα′;
结论:α′=sin-1[(BB′-AA′)/AB],
①反射面走向角∠α≈∠α′;
②β=sin-1[(DD′-CC′)/CD],
反射面倾角γ=90°-β;
③BB″=BB′/cosβ,
BBl≈BB′/cosβcosα′;
④AAl≈AA′/cosβcosα′;
⑤若DD′>CC′,则β>0,反射面倾向开挖方向,
若DD′<CC′,则β<0,反射面倾向反开挖方向。
计算完成。
这是通过在一定误差范围内的简化得出的。其条件是反射面倾角在70°~90°间,倾角小于70°则误差加大,此计算方法经过对超过100条断层的计算,在隧道开挖后实际丈量,误差在可容许范围内。
结合附图,本发明举例如下:
已知:波速V=3900m/s
A点到N面的反射时间TAA′=21.2ms;
B点到N面的反射时间TBB′=22.1ms;
C点到N面的反射时间TCC′=22.6ms;
D点到N面的反射时间TDD′=22.8ms;
AB=11.4m;
CD=2.1m;
AA′=V×TAA′=41.3m;
BB′=V×TBB′=43.1m;
CC′=V×TCC′=44.1m;
DD′=V×TDD′=44.5m;
得到反射面走向角α′=sin-1[(BB′-AA′)/AB]
                    =sin-1[(43.1m-41.3m)/11.4m]=8°;
β=sin-1[(DD′-CC′)/CD],
  =sin-1[(44.5m-44.1m)/2.1m]=11.2°;
反射面倾角γ=90°-β=78°;
AAl≈AA′/cosβcosα′
   =43.1m/cos11.2°cos8°=42.6m;
BBl≈BB′/cosβcosα′
   =41.3m/cos11.2°cos8°=44.5m;
在计算出地质体反射面的走向角和倾角、倾向就得到了地质体反射面的产状;计算出AAl和BBl后,即可得到地质体反射面与通过AB的水平面的交线的位置。从而确定反射面的空间位置。
若反射面不一定是平面,可计算AB测线上4~10个测点的AAl和BBl长度来确定地质体反射面与通过AB线的水平面的交线。

Claims (1)

1.一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法,是在掌子面上设定一条水平测线AB和一条与其交叉的铅垂向测线CD,在测线上每隔20~40cm设一测点,采用20~40cm的震-检距激震和接收弹性波测出每个测点的时间曲线,然后将一条测线各点时间曲线汇成时间剖面,在时间剖面上得到测点到反射面的反射时间T;在已知波速V的情况下,计算出反射面的走向、倾角和倾向,并确定它在隧道中的空间位置,其特征在于,计算掌子面前方地质体的产状按下列步骤进行:
第一步,在掌子面M上的水平测线AB和铅垂测线CD的基础上与反射面N之间建立计算模型图,标示出测线、射线、计算辅助线及反射面N的走向角α、倾角γ、辅助计算角α和β的相对位置,线AA′为A点的入射波射线,点A′为反射点;线BB′为B点的入射波射线,点B′为反射点;线CC′为C点的入射波射线,点C′为反射点;线DD′为D点的入射波射线,点D′为反射点;面
Figure FSB00000725025600011
为通过测线AB的水平面,垂直于掌子面M,与反射面N交于
Figure FSB00000725025600012
线;面CC′D′D为通过铅垂测线CD的垂直面,垂直于掌子面M,与反射面N交于C′D′线;点Ap是通过点A′并与测线AB平行的线A′Ap与线BB′的交点;点Bp是通过点
Figure FSB00000725025600013
并与测线AB平行的线
Figure FSB00000725025600014
与线BB′的交点;点D″是角D′C′D″的一边C′D″与线DD′的交点,C′D″垂直于DD″;点B″是面BB′B″垂直于反射面N的线BB″与线
Figure FSB00000725025600015
的交点;线BB″是面BB′B″和面
Figure FSB00000725025600016
在反射面N上的交线;
第二步,根据L=V·T/2的数学公式分别计算出各测点至反射面的距离,式中:L-两点间距,V-波速,t-测点到反射面的反射时间;
A点的入射波射线长度AA′=V·TAA/2,
B点的入射波射线长度BB′=V·TBB′/2,
C点的入射波射线长度CC′=V·TCC′/2,
D点的入射波射线长度DD′=V·TDD′/2,
第三步,反射面的走向角α、倾角γ的计算:
①反射面的走向角α的计算
∠B′A′Ap=α′是平面AA′B′B上的一个计算辅助角,
Figure FSB00000725025600017
是反射面N的走向角,
由于
Figure FSB00000725025600018
以及A′A//B′B,
所以sinα′=(BB′-AA′)/A′Ap=(BB′-AA′)/AB
α′=sin-1((BB′-AA′)/AB)
Figure FSB00000725025600021
Figure FSB00000725025600022
由于面A′B′Ap与面
Figure FSB00000725025600023
不平行,有一相互夹角β,
所以∠α′>∠α;
②反射面N的倾角γ的计算
γ=90°-β
β=∠D′D′D″,是C′C′线与C′D″线的夹角,
由于
Figure FSB00000725025600024
所以sinβ=((DD′-CC′)/C′D″;
又由于C′D″=CD,
所以β=sin-1[(DD′-CC′)/CD]
在工程实践中,当∠β在0°~20°和∠α在0°~30°的情况下,计算反射面产状时,可以接受∠α≈∠α′
第四步,反射面N倾向的判断
由于β=sin-1[(DD′-CC′)/CD]
所以若DD′>CC′,则β>0,反射面倾向开挖方向,
若DD′<CC′,则β<0,反射面倾向反开挖方向;
根据与反射面的走向角α、倾角γ和倾向即可确定反射面的空间位置;
若反射面不一定是平面,可计算测线AB上4~10个测点的
Figure FSB00000725025600026
Figure FSB00000725025600027
长度来确定地质体反射面与通过测线AB的水平面的交线,按照上述方法和步骤进行计算,判定地质体反射面的产状和空间位置。
CN2010102666019A 2010-08-30 2010-08-30 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法 Active CN101950030B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102666019A CN101950030B (zh) 2010-08-30 2010-08-30 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102666019A CN101950030B (zh) 2010-08-30 2010-08-30 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法

Publications (2)

Publication Number Publication Date
CN101950030A CN101950030A (zh) 2011-01-19
CN101950030B true CN101950030B (zh) 2012-05-23

Family

ID=43453593

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102666019A Active CN101950030B (zh) 2010-08-30 2010-08-30 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法

Country Status (1)

Country Link
CN (1) CN101950030B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105738471B (zh) * 2014-12-11 2018-06-19 安徽惠洲地质安全研究院股份有限公司 一种智能电子敲帮问顶便携式仪器
CN105891873B (zh) * 2016-06-17 2018-07-17 钟世航 一种采用弹性波反射法探查地下岩体中流体和气体的方法
CN106780741B (zh) * 2017-02-10 2020-05-12 中国电建集团成都勘测设计研究院有限公司 基于多点产状数据的地质曲面快速拟合插值方法
CN111337979B (zh) * 2020-03-19 2022-08-16 辽宁工程技术大学 一种用解析方法确定地质界面真倾角的方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101672929B (zh) * 2009-09-23 2011-04-20 山东大学 在隧道中探查掌子面前方含水层并估算其出水量的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10260167A (ja) * 1997-03-18 1998-09-29 Mitsubishi Heavy Ind Ltd 壁面検査装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101672929B (zh) * 2009-09-23 2011-04-20 山东大学 在隧道中探查掌子面前方含水层并估算其出水量的方法

Also Published As

Publication number Publication date
CN101950030A (zh) 2011-01-19

Similar Documents

Publication Publication Date Title
CN106088166B (zh) 一种管桩测斜装置及其测斜管居中定位组件
CN103148833B (zh) 岩层产状参数获取和计算方法
CN101950030B (zh) 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法
CN103697854A (zh) 一种非接触式结构面产状测量方法
WO2013067590A1 (en) Vibration analysis for blasting
CN102536196A (zh) 基于激光测距和加速度测量的井下姿态测量系统及方法
US9977124B2 (en) Method and system for characterization of subsurface cavities using joint inversion of gravity and ground penetrating radar data
CN103852059B (zh) 反铲挖掘机的铲斗定位装置及方法
CN102012509B (zh) 地质结构面产状的雷达探测方法
CN101943599B (zh) 一种在隧道中用弹性波反射法时测取掌子面前方岩体波速的方法
CN107165678A (zh) 一种深埋地下供水管线的爆破振动监测方法
CN104199109A (zh) 确定钻井目标层位的视倾角的方法和设备
CN103389523A (zh) 隧道地质的预报方法和系统
Puente et al. Monitoring of progressive damage in buildings using laser scan data
Jinyi et al. A health monitoring system for inverted arch of salt rock tunnel based on laser level deformation monitor and wFBG
CN103399356A (zh) 隧道地质的预报方法和系统
CN102749066A (zh) 一种电子激光靶及其测量装置
CN104535041B (zh) 观测站位置确定方法和高层建筑物倾斜变形测定方法
CN103837604B (zh) 一种岩石爆破损伤跨孔声波测试中跨距的修正方法
CN105783849B (zh) 一种岩层真倾角和真走向的测量系统及方法
CN108104822A (zh) 一种隧道穿越大型岩溶空洞施工地质灾害实时预测方法
Song et al. Determination of RVE size based on the 3D fracture persistence
CN107152964A (zh) 一种隧道台阶法施工中深埋地下供水管线的爆破振动监测方法
CN103389524B (zh) 隧道地质的预报方法和系统
JP7419119B2 (ja) 作業機械

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
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: No. 5 building 100081 Beijing city Haidian District Xizhimen Daliushu Road No. 17 Fuhai center room 1603

Patentee after: Zhong Shihang

Patentee after: Wang Rong

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

Patentee before: Zhong Shihang

Patentee before: Wang Rong

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 Road, No. 3 hospital of Jiaotong University Yuen Cheong Park Building 6, room 1408

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

DD01 Delivery of document by public notice

Addressee: Zhong Shihang

Document name: payment instructions

DD01 Delivery of document by public notice