CN113626922B - 基于地质雷达的隧道顶板富水区域结构涌水量预测方法 - Google Patents

基于地质雷达的隧道顶板富水区域结构涌水量预测方法 Download PDF

Info

Publication number
CN113626922B
CN113626922B CN202110935113.0A CN202110935113A CN113626922B CN 113626922 B CN113626922 B CN 113626922B CN 202110935113 A CN202110935113 A CN 202110935113A CN 113626922 B CN113626922 B CN 113626922B
Authority
CN
China
Prior art keywords
water
rich
rich area
point
water inflow
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
CN202110935113.0A
Other languages
English (en)
Other versions
CN113626922A (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.)
Yunnan Weishidun Technology Co ltd
Kunming University of Science and Technology
Original Assignee
Yunnan Weishidun Technology Co ltd
Kunming University of Science and Technology
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 Yunnan Weishidun Technology Co ltd, Kunming University of Science and Technology filed Critical Yunnan Weishidun Technology Co ltd
Priority to CN202110935113.0A priority Critical patent/CN113626922B/zh
Publication of CN113626922A publication Critical patent/CN113626922A/zh
Application granted granted Critical
Publication of CN113626922B publication Critical patent/CN113626922B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Business, Economics & Management (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Strategic Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Human Resources & Organizations (AREA)
  • Data Mining & Analysis (AREA)
  • Economics (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Evolutionary Computation (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Structural Engineering (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Civil Engineering (AREA)
  • Algebra (AREA)
  • Architecture (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及基于地质雷达的隧道顶板富水区域结构涌水量预测方法,属于顶板涌水量预测技术领域。本发明首先通过地质雷达超前探测收集富水区域结构资料,根据所获资料建立天然状态下的隧道顶板富水区域结构模型,根据富水区域结构参数,修正涌水量预测公式进行涌水量预测。本发明以地质雷达所获物探资料为基础,在此基础上划分富水区域建立天然状态下的隧道顶板富水区域结构模型,根据富水区域结构参数,修正涌水量预测公式,并建立达到富水区域的涌水量预测公式,其预测曲线更能反映实际工况,在短距离涌水量预测方面其预测精度更高。

Description

基于地质雷达的隧道顶板富水区域结构涌水量预测方法
技术领域
本发明涉及基于地质雷达的隧道顶板富水区域结构涌水量预测方法,属于顶板涌水量预测技术领域。
背景技术
随着我国经济的高速发展,人们对于铁路、公路交通的需求和要求不断提升,现在我国各领域运营、开挖隧道已接近5万千米,但在滇南地区,由于地质条件的复杂性,在隧道开挖过程中,会遇到溶洞、断层、破碎带、涌水等各种不良地质现象,面临的问题的挑战性也越来越高,其中涌水灾害,在国内外隧道特大施工中,其发生次数和死亡人数均在前列,影响隧道工期延误,危害施工人员生命安全,因此在施工前根据隧道勘测资料,对掌子面顶板前方富水不良地质体进行探测进行判断,根据其结构建立结构模型,修正涌水量预测公式,及时给出涌水量预测,进而对堵、排水及时给出方案,使设计人员根据所获信息及时更改设计方案,指导隧道施工,保障隧道施工安全。
发明内容
本发明提供了基于地质雷达的隧道顶板富水区域结构涌水量预测方法,以用于修正涌水量预测公式,及时给出涌水量预测,进而对堵、排水及时给出方案提供支撑。
本发明的技术方案是:基于地质雷达的隧道顶板富水区域结构涌水量预测方法,首先通过地质雷达超前探测收集富水区域结构资料,根据所获资料建立天然状态下的隧道顶板富水区域结构模型,根据富水区域结构参数,修正涌水量预测公式进行涌水量预测。
进一步地,所述方法的具体步骤如下:
Step1、通过在掌子面前布线进行雷达探测测取掌子面前0-30m富水区域雷达正演图,并针对其富水区域结构特征圈定富水区域规模;
Step2、根据上述所获地质雷达探测资料,划分富水区域建立天然状态下的隧道顶板富水区域结构模型;
Step3、根据建立的隧道顶板富水区域结构模型中的富水区域结构参数,修正涌水量预测公式,建立达到富水区域的涌水量预测公式,进行涌水量预测。
进一步地,所述Step2中,划分富水区域建立天然状态下的隧道顶板富水区域结构模型包括:
设在介电常数∈1的空间内存在一个介电常数为∈2的椭圆体,距离掌子面距离为h,富水区域长轴为a,富水区域短轴为b,电磁波发射点为P,O点为电磁波到达H点时PH发射点的中点,Q点为椭圆区域最远反射点,设PH=x1,PQ=x2,OP=y;
整理后得
式(2)为富水区域圆锥曲线模型,P为动点通过在Y轴上的移动使得到的a、b值更加的精确,式中x1、x2、y均能通过雷达探测可得。
进一步地,所述Step3中,修正后的涌水量预测公式为:
其中,q0、q1为隧道可能的最大初期涌水量(m3/d);K为岩体的渗透系数;H为含水层厚度(m);a为富水区域短轴;b为富水区域长轴;r为隧道半径,L为掌子面到富水区域的距离(m),N为富水区域厚度(m)。
本发明的有益效果是:
1、本发明以地质雷达所获物探资料为基础,在此基础上划分富水区域建立天然状态下的隧道顶板富水区域结构模型,根据富水区域结构参数,修正涌水量预测公式,并建立达到富水区域的涌水量预测公式,其预测曲线更能反映实际工况,在短距离涌水量预测方面其预测精度更高。
2、本富水区域结构涌水量预测方法为针对滇南地区特深长隧道顶板涌水频发现象所设立,为减少涌水事故所造成的工期延误,危害施工人员生命安全,从而为隧道掘进过程提供保障,根据涌水量预测值大小及时反馈业主。
附图说明
图1为本发明的结构流程;
图2为富水区域结构模型;
图3为测线布置图;
图4为测线1从右至左地质雷达成果图;
图5为测线2从左至右地质雷达成果图。
具体实施方式
实施例1:如图1-5所示,基于地质雷达的隧道顶板富水区域结构涌水量预测方法,首先通过地质雷达超前探测收集富水区域结构资料,根据所获资料建立天然状态下的隧道顶板富水区域结构模型,根据富水区域结构参数,修正涌水量预测公式进行涌水量预测。
进一步地,所述方法的具体步骤如下:
Step1、通过在掌子面前布线进行雷达探测测取掌子面前0-30m富水区域雷达正演图,并针对其富水区域结构特征圈定富水区域规模;布线方式如图3所示;得到的测线1、测线2两张地质雷达成果图,如图4、图5所示;通过对雷达图的分析可知,在6-10米区域内呈现双曲线波形,且该处的电磁信号频率较低,振幅较强,具有明显的双曲线状顶底界面反射,且多次反射,含水量较大,含水区域呈椭圆形分布;
Step2、根据上述所获地质雷达探测资料,划分富水区域建立天然状态下的隧道顶板富水区域结构模型;
进一步地,所述Step2中,划分富水区域建立天然状态下的隧道顶板富水区域结构模型包括:
设在介电常数∈1的空间内存在一个介电常数为∈2的椭圆体,距离掌子面距离为h,富水区域长轴为a,富水区域短轴为b,电磁波发射点为P,O点为电磁波到达H点时PH发射点的中点,Q点为椭圆区域最远反射点,设PH=x1,PQ=x2,OP=y;
整理后得
式(2)为富水区域圆锥曲线模型,P为动点通过在Y轴上的移动使得到的a、b值更加的精确,式中x1、x2、y均能通过雷达探测可得。划分富水区域建立天然状态下的隧道顶板富水区域结构模型,如图2所示;
具体的,此式中可以进行如下取值,x2=7、y=16.8,根据富水区域结构模型取到a=3.09,带入得b=20.08;
Step3、根据建立的隧道顶板富水区域结构模型中的富水区域结构参数,修正涌水量预测公式,建立达到富水区域的涌水量预测公式,进行涌水量预测。
进一步地,所述Step3中,修正后的涌水量预测公式为:
其中,q0、q1为隧道可能的最大初期涌水量(m3/d);K为岩体的渗透系数;H为含水层厚度(m);a为富水区域短轴;b为富水区域长轴;r为隧道半径,L为掌子面到富水区域的距离(m),N为富水区域厚度(m)。
具体的,K为岩体的渗透系数为5×10-5m·s-1
H为含水层厚度180m;a为富水区域短轴3.09m;b为富水区域长轴20.08m;r为隧道半径3m,L为掌子面到富水区域的距离3m,N为富水区域厚度5m;m为转换系数一般取0.86。
根据雷达图像所划富水区域厚度为5m当L为3m时带入得q0=907.606m3/d,当N为5m时q1=6433.437m3/d,q1为到达富水区域后隧道开挖后的可能的最大初期涌水量,q0值在未到达富水区域时其涌水量较少,q1值反映了在隧道挖掘过程中在接触富水区域后涌水量突增的情况,与实际开挖时遇到富水区域后涌水量突增的情况更加贴切,应将此值及时向业主反馈,做好排-堵的防护,若q1值过大应通知设计人员及时更改施工方案。
上面结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (1)

1.基于地质雷达的隧道顶板富水区域结构涌水量预测方法,其特征在于:
所述方法的具体步骤如下:
Step1、通过在掌子面前布线进行雷达探测测取掌子面前0-30m富水区域雷达正演图,并针对富水区域结构特征圈定富水区域规模;
Step2、通过地质雷达超前探测收集富水区域结构资料,根据上述所获地质雷达探测资料,划分富水区域建立天然状态下的隧道顶板富水区域结构模型;
Step3、根据建立的隧道顶板富水区域结构模型中的富水区域结构参数,修正涌水量预测公式,建立达到富水区域的涌水量预测公式,进行涌水量预测;
所述Step2中,划分富水区域建立天然状态下的隧道顶板富水区域结构模型包括:
设在介电常数∈1的空间内存在一个介电常数为∈2的椭圆体,距离掌子面距离为h,富水区域长轴为a,富水区域短轴为b,电磁波发射点为P,H点为椭圆上靠近P点曲率最大的点,O点为P点发射的电磁波到达H点时,发生反射,反射波落在Y轴的上半轴上的点与P点连线的中点,Q点为椭圆区域最远反射点,设PH=x1,PQ=x2,OP=y;
整理后得
式(2)为富水区域圆锥曲线模型,P为动点通过在Y轴上的移动使得到的a、b值更加的精确,式中x1、x2、y均能通过雷达探测可得。
CN202110935113.0A 2021-08-16 2021-08-16 基于地质雷达的隧道顶板富水区域结构涌水量预测方法 Active CN113626922B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110935113.0A CN113626922B (zh) 2021-08-16 2021-08-16 基于地质雷达的隧道顶板富水区域结构涌水量预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110935113.0A CN113626922B (zh) 2021-08-16 2021-08-16 基于地质雷达的隧道顶板富水区域结构涌水量预测方法

Publications (2)

Publication Number Publication Date
CN113626922A CN113626922A (zh) 2021-11-09
CN113626922B true CN113626922B (zh) 2023-07-28

Family

ID=78385491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110935113.0A Active CN113626922B (zh) 2021-08-16 2021-08-16 基于地质雷达的隧道顶板富水区域结构涌水量预测方法

Country Status (1)

Country Link
CN (1) CN113626922B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017181253A (ja) * 2016-03-30 2017-10-05 大成建設株式会社 地下水流の数値解析方法
CN107367772A (zh) * 2017-08-29 2017-11-21 西南石油大学 一种预测隧道前方不良地质的超前预报方法
CN110795878A (zh) * 2019-10-25 2020-02-14 中国铁路设计集团有限公司 一种隧道涌水量预测方法
CN111191849A (zh) * 2020-01-03 2020-05-22 中煤科工集团西安研究院有限公司 一种西部矿区深埋工作面涌水量预测方法
CN112147671A (zh) * 2020-09-23 2020-12-29 四川振通公路工程检测咨询有限公司 隧道含导水不良地质构造综合超前预报方法
CN112177617A (zh) * 2020-09-25 2021-01-05 中铁二十局集团有限公司 高压富水断层隧道施工超前地质预报预测方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102141552B1 (ko) * 2019-12-02 2020-08-06 대한민국 터널에서 사고 발생을 방지하기 위한 방법 및 장치

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017181253A (ja) * 2016-03-30 2017-10-05 大成建設株式会社 地下水流の数値解析方法
CN107367772A (zh) * 2017-08-29 2017-11-21 西南石油大学 一种预测隧道前方不良地质的超前预报方法
CN110795878A (zh) * 2019-10-25 2020-02-14 中国铁路设计集团有限公司 一种隧道涌水量预测方法
CN111191849A (zh) * 2020-01-03 2020-05-22 中煤科工集团西安研究院有限公司 一种西部矿区深埋工作面涌水量预测方法
CN112147671A (zh) * 2020-09-23 2020-12-29 四川振通公路工程检测咨询有限公司 隧道含导水不良地质构造综合超前预报方法
CN112177617A (zh) * 2020-09-25 2021-01-05 中铁二十局集团有限公司 高压富水断层隧道施工超前地质预报预测方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
地质雷达正演模拟及在 断层富水带超前地质预报的应用研究;蒋建国等;《铁道科学与工程学报》;第16卷(第11期);2801-2808 *
富水岩层隧道区域涌水量预测方法及工程应用;王健华;李术才;李利平;许振浩;石少帅;;人民长江(14);40-45 *

Also Published As

Publication number Publication date
CN113626922A (zh) 2021-11-09

Similar Documents

Publication Publication Date Title
CN109613616B (zh) 一种浅层岩溶隧道地质综合立体探测方法
CN101980053B (zh) 一种复杂礁滩储层预测方法
CN103033846B (zh) 地质相控制的地震反演系统和地震反演方法
CN111948645B (zh) 煤矿巷道和隧道随钻钻孔雷达超前探测装置与方法
CN104749617A (zh) 一种多尺度裂缝储层正演模型建立方法
CN113627080B (zh) 一种基于支持向量机的煤矿掘进工作面风量需求预测方法
CN103912270A (zh) 井间测井深度与速度同步控制定位系统及方法
US20230184983A1 (en) Vector-resistivity-based real-time advanced detection method for water-bearing hazard body
CN116030207A (zh) 一种岩溶地质公路隧道施工综合超前三维地质建模方法
CN103529486A (zh) 一种地球化学异常圈定方法
CN105005081A (zh) 煤机采动激励下综采面近场煤岩动态层析成像系统及方法
CN113626922B (zh) 基于地质雷达的隧道顶板富水区域结构涌水量预测方法
CN111830499A (zh) 基于电磁波技术的煤岩识别方法
CN112285802B (zh) 隧道地震和瞬变电磁联合探测方法
CN114839671B (zh) 地面和地巷瞬变电磁法精细识别煤系地层电性界面方法
CN110687591B (zh) 基于先验数据的波形匹配确定煤层及围岩物性参数的方法
CN105220694A (zh) 一种工程抛石量估算方法
CN115586577A (zh) 一种定源瞬变电磁非中心点观测数据全时转换方法
CN109683193B (zh) 一种双洞隧道施工前方不良地质体的主动探测方法
Zhao et al. Case study of ground penetrating radar for geological prediction in Qiyueshan Tunnel
CN111596372B (zh) 一种基于电磁波测量系统的含水孔隙度及矿化度反演方法
CN103576189B (zh) 一种随钻地震直达波组搜索方法
CN107894611B (zh) 一种水平井数据转直井数据的方法及用于地震反演的方法
Wang et al. Application of LTD-2100 ground penetrating radar to advanced forecast in diversion tunnel of jinping hydroelectric power station
CN116068627A (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