CN116203633A - 基于天然源震动成像的隧道围岩完整性评价方法 - Google Patents
基于天然源震动成像的隧道围岩完整性评价方法 Download PDFInfo
- Publication number
- CN116203633A CN116203633A CN202310249797.8A CN202310249797A CN116203633A CN 116203633 A CN116203633 A CN 116203633A CN 202310249797 A CN202310249797 A CN 202310249797A CN 116203633 A CN116203633 A CN 116203633A
- Authority
- CN
- China
- Prior art keywords
- surrounding rock
- tunnel
- impedance value
- integrity
- natural source
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 101
- 238000003384 imaging method Methods 0.000 title claims abstract description 26
- 238000011156 evaluation Methods 0.000 title claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims description 10
- 238000000034 method Methods 0.000 abstract description 12
- 238000005553 drilling Methods 0.000 abstract description 5
- 239000011148 porous material Substances 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005070 sampling 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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
本发明公开了一种基于天然源震动成像的隧道围岩完整性方法,包括:采集天然源震动数据,提取特征频率,通过时深转换和成像数据处理得到隧道围岩阻抗值,确定隧道围岩的阻抗变化率、完整基岩基准阻抗值、破碎体基准阻抗值、归一化阻抗值、计算隧道围岩完整性系数,进而对隧道围岩完整性进行评价。该方法所用的勘探设备轻便,不受场地和人文干扰因素的影响,场地适应性强,无需成孔、计算方便快捷、勘察成本低,勘察效率高,在艰险复杂山区等调绘、钻探、波速实施困难的深埋隧道围岩评价方面优势明显,具有较大的使用推广价值。
Description
技术领域
本发明涉及工程地质勘察领域,尤其是涉及一种基于天然源震动成像的隧道围岩完整性评价方法。
背景技术
围岩的完整性评价是铁路、公路等隧道工程设计的关键环节,目前主要是通过调绘、钻探和波速测试等方法。调绘只能获取地层表面的地质信息,间接对隧道围岩完整性进行评价,准确度不高,同时难以对深部隧道完整性进行评价;钻探取样成本高、周期长,在地形条件复杂的山区实施困难,且是一孔之见,对横向变化大的隧道围岩误差较大。波速测试需要同时测试岩体和岩块的纵波速度,现场测试较为困难,尤其是深部隧道围岩的测试尤其困难,难以获得准确的岩体和岩块速度。
发明内容
针对现有技术存在的问题,本发明提供一种计算过程简单、快捷且计算精度较高的基于天然源震动成像的隧道围岩完整性评价方法。
为此,本发明采用以下技术方案:
一种基于天然源震动成像的隧道围岩完整性评价方法,包括以下步骤:
S1,场地天然源震动成像数据采集:将多个三分量节点式地震仪等间距线性排列布置,并同时采集一定时间内的天然源震动数据;
S2,天然源震动成像数据处理:分别对每个三分量节点式地震仪采集的天然源震动数据进行处理,提取特征频率,通过时深转换和成像,获得测试区域的天然源震动成像阻抗成果,并确定隧道围岩类型;
S3,计算步骤S2确定的隧道围岩类型中某一类隧道围岩的第i个计算点的阻抗变化率Si;
S4,确定该类隧道围岩的完整基岩基准阻抗值Iw:
其中:
Iij为该类围岩中阻抗变化率≥0.9的第j个阻抗值,单位为g/m2·s;
n1为该类围岩中阻抗变化率≥0.9的计算点数;
S5,确定该类隧道围岩的破碎体基准阻抗值Io:
其中,
Iik为该类围岩中阻抗变化率≤0.1的第k个阻抗值,单位为g/m2·s;
n2为该类围岩中阻抗变化率≤0.1的计算点数;
S6,利用步骤S3得到的隧道围岩阻抗变化率Si和步骤S5得到的基准阻抗值I0计算该类隧道围岩第i个计算点的归一化阻抗值INi;
S7,利用步骤S6得到的隧道围岩归一化阻抗值INi计算该类隧道围岩的完整性系数KIi:
KIi=INi/IW
其中:
INi的单位为g/m2·s;
IW的单位为g/m2·s;
S8,对该类隧道围岩的完整性进行评价;
S9,重复步骤S3-S8,完成对所有类型隧道围岩的完整性评价。
优选的是,步骤S1中所述的一定时间为2小时。
步骤S3中,所述阻抗变化率Si通过下式计算:
其中,Ii为该类围岩第i个计算点的阻抗值,单位为g/m2·s;
Imin为该类围岩的最小阻抗值,单位为g/m2·s;
Imax为该类围岩的最大阻抗值,单位为g/m2·s。
步骤S6中所述的归一化阻抗值INi的计算公式为:INi=Si·(IW-Io)+Io。
步骤S8的评价标准为:
当KIi≥0.9时,第i个计算点处围岩完整性为I级;
当0.65≤KIi<0.9时,第i个计算点处围岩完整性为II级;
当0.35≤KIi<0.65时,第i个计算点处围岩完整性为III级;
当0.1≤KIi<0.35时,第i个计算点处围岩完整性为IV级;
当KIi<0.1时,第i个计算点处围岩完整性为V级。
本发明的基于天然源震动成像的隧道围岩完整性评价方法,通过在地面设计合理的观测系统,采集天然源震动数据,提取特征频率,通过时深转换和成像数据处理得到隧道围岩阻抗值,确定隧道围岩阻抗变化率、完整基岩基准阻抗值、破碎体基准阻抗值、归一化阻抗值、计算隧道围岩完整性系数,进而对隧道围岩完整性进行评价。
与现有技术相比,本发明具有以下有益效果:
1.本发明直接应用天然源震动成像勘探结果进行隧道围岩完整性评价,计算所需参数来源于天然源震动成像勘探数据,无需调绘、钻探、波速测试等,计算过程简单、快捷,精度高,计算效率提高50%以上。
2.该方法所用勘探设备轻便,勘察成本低,勘探成本节约50-60%。
3.该方法测试过程中勘探设备放置在地面上即可,采集的是天然源震动信号,不受场地和人文干扰因素的影响,场地适应性强,在艰险复杂山区等调绘、钻探、波速实施困难的深埋隧道围岩完整性评价方面优势明显,具有较大的使用推广价值。
附图说明
图1为本发明的基于天然源震动成像的隧道围岩完整性评价方法的流程图。
具体实施方式
下面结合附图对本发明的评价方法进行详细说明。
如图1所示,本发明的基于天然源震动成像的隧道围岩完整性评价方法包括以下步骤:
S1,场地天然源震动成像数据采集:将多个三分量节点式地震仪等间距线性排列布置,并同时采集一定时间(例如2小时)内的天然源震动数据;
S2,天然源震动成像数据处理:分别对每个三分量节点式地震仪采集的天然源震动数据进行处理,提取特征频率,通过时深转换和成像,获得测试区域的天然源震动成像阻抗值成果,并确定隧道围岩类型;
S3,计算步骤S2确定的隧道围岩类型中某一类隧道围岩的第i个计算点的阻抗变化率Si,计算公式如下:
其中,Ii为该类围岩第i个计算点的阻抗值,单位为g/m2·s;
Imin为该类围岩的最小阻抗值,单位为g/m2·s;
Imax为该类围岩的最大阻抗值,单位为g/m2·s。
S4,确定该类隧道围岩的完整基岩基准阻抗值Iw:
其中:
Iij为该类围岩中阻抗变化率≥0.9的第j个阻抗值,单位为g/m2·s;
n1为该类围岩中阻抗变化率≥0.9的计算点数;
S5,确定该类隧道围岩的破碎体基准阻抗值Io:
其中,
Iik为该类围岩中阻抗变化率≤0.1的第k个阻抗值,单位为g/m2·s;
n2为该类围岩中阻抗变化率≤0.1的计算点数;
S6,利用步骤S3得到的隧道围岩阻抗变化率Si和步骤S5得到的基准阻抗值I0计算该类隧道围岩第i个计算点的归一化阻抗值INi:
INi=Si·(IW-Io)+Io。
S7,利用步骤S6得到的隧道围岩归一化阻抗值INi计算该类隧道围岩的完整性系数KIi:
KIi=INi/IW
其中:
INi的单位为g/m2·s;
IW的单位为g/m2·s;
S8,对该类隧道围岩的完整性进行评价:
当KIi≥0.9时,第i个计算点处围岩完整性为I级;
当0.65≤KIi<0.9时,第i个计算点处围岩完整性为II级;
当0.35≤KIi<0.65时,第i个计算点处围岩完整性为III级;
当0.1≤KIi<0.35时,第i个计算点处围岩完整性为IV级;
当KIi<0.1时,第i个计算点处围岩完整性为V级。
S9,重复步骤S3-S8,完成对所有类型隧道围岩的完整性评价。
Claims (5)
1.一种基于天然源震动成像的隧道围岩完整性评价方法,包括以下步骤:
S1,场地天然源震动成像数据采集:将多个三分量节点式地震仪等间距线性排列布置,并同时采集一定时间内的天然源震动数据;
S2,天然源震动成像数据处理:分别对每个三分量节点式地震仪采集的天然源震动数据进行处理,提取特征频率,通过时深转换和成像,获得测试区域的天然源震动成像阻抗成果,并确定隧道围岩类型;
S3,计算步骤S2确定的隧道围岩类型中某一类隧道围岩的第i个计算点的阻抗变化率Si;
S4,确定该类隧道围岩的完整基岩基准阻抗值Iw:
其中:
Iij为该类围岩中阻抗变化率≥0.9的第j个阻抗值,单位为g/m2·s;
n1为该类围岩中阻抗变化率≥0.9的计算点数;
S5,确定该类隧道围岩的破碎体基准阻抗值Io:
其中,
Iik为该类围岩中阻抗变化率≤0.1的第k个阻抗值,单位为g/m2·s;
n2为该类围岩中阻抗变化率≤0.1的计算点数;
S6,利用步骤S3得到的隧道围岩阻抗变化率Si和步骤S5得到的基准阻抗值I0计算该类隧道围岩第i个计算点的归一化阻抗值INi;
S7,利用步骤S6得到的隧道围岩归一化阻抗值INi计算该类隧道围岩的完整性系数KIi:
KIi=INi/IW,
其中:
INi的单位为g/m2·s;
IW的单位为g/m2·s;
S8,对该类隧道围岩的完整性进行评价;
S9,重复步骤S3-S8,完成对所有类型隧道围岩的完整性评价。
2.根据权利要求1所述的隧道围岩完整性评价方法,其特征在于:步骤S1中,所述一定时间为2小时。
4.根据权利要求1所述的隧道围岩完整性评价方法,其特征在于:步骤S6所述的归一化阻抗值INi的计算公式为:INi=Si·(IW-Io)+Io。
5.根据权利要求1所述的隧道围岩完整性评价方法,其特征在于:步骤S8的评价标准为:
当KIi≥0.9时,第i个计算点处围岩完整性为I级;
当0.65≤KIi<0.9时,第i个计算点处围岩完整性为II级;
当0.35≤KIi<0.65时,第i个计算点处围岩完整性为III级;
当0.1≤KIi<0.35时,第i个计算点处围岩完整性为IV级;
当KIi<0.1时,第i个计算点处围岩完整性为V级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310249797.8A CN116203633B (zh) | 2023-03-15 | 2023-03-15 | 基于天然源震动成像的隧道围岩完整性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310249797.8A CN116203633B (zh) | 2023-03-15 | 2023-03-15 | 基于天然源震动成像的隧道围岩完整性评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116203633A true CN116203633A (zh) | 2023-06-02 |
CN116203633B CN116203633B (zh) | 2024-03-22 |
Family
ID=86514596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310249797.8A Active CN116203633B (zh) | 2023-03-15 | 2023-03-15 | 基于天然源震动成像的隧道围岩完整性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116203633B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117075199A (zh) * | 2023-08-18 | 2023-11-17 | 中国铁路设计集团有限公司 | 基于天然源共振成像的既有铁路路基注浆效果评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110398784A (zh) * | 2019-08-07 | 2019-11-01 | 中铁四局集团第四工程有限公司 | 一种隧道穿越断层破碎带围岩稳定性评价方法 |
CN113010942A (zh) * | 2021-02-25 | 2021-06-22 | 中国铁路设计集团有限公司 | 隧道开挖风险预警及围岩分级评价方法 |
-
2023
- 2023-03-15 CN CN202310249797.8A patent/CN116203633B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110398784A (zh) * | 2019-08-07 | 2019-11-01 | 中铁四局集团第四工程有限公司 | 一种隧道穿越断层破碎带围岩稳定性评价方法 |
CN113010942A (zh) * | 2021-02-25 | 2021-06-22 | 中国铁路设计集团有限公司 | 隧道开挖风险预警及围岩分级评价方法 |
Non-Patent Citations (6)
Title |
---|
夏金儒;: "地震映像和瞬变电磁法在采区空洞勘查中的综合应用", 地质找矿论丛, no. 1, 30 October 2006 (2006-10-30) * |
彭功勋;赵旭光;张陆军;: "管波探测法在广州西北部复杂岩体桩基勘察中的应用", 广东土木与建筑, no. 11, 25 November 2018 (2018-11-25) * |
李光旭;顾汉明;詹凤林;刘振明;: "地质雷达结合地震映像法在铁路隧道底部隐伏岩溶探测中的应用", 工程地球物理学报, no. 06, 30 December 2009 (2009-12-30) * |
李志华;: "铁路煤矿采空区的综合物探研究", 物探与化探, no. 1, 15 October 2012 (2012-10-15) * |
李法滨;俞仁泉;何文勇;: "综合物探在复杂岩性地区断层勘察中的应用", 公路交通技术, no. 01, 25 February 2013 (2013-02-25) * |
郝勇;王建军;: "高密度电法及瞬变电磁法在工程地质勘察中的综合应用", 中国水运(下半月), no. 10, 15 October 2012 (2012-10-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117075199A (zh) * | 2023-08-18 | 2023-11-17 | 中国铁路设计集团有限公司 | 基于天然源共振成像的既有铁路路基注浆效果评价方法 |
CN117075199B (zh) * | 2023-08-18 | 2024-03-22 | 中国铁路设计集团有限公司 | 基于天然源共振成像的既有铁路路基注浆效果评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116203633B (zh) | 2024-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102373923B (zh) | 一种储层识别方法 | |
CN113295850A (zh) | 基于多源数据融合的隧道围岩定量快速分级方法与装置 | |
CN110823962B (zh) | 一种滑坡体的三维成像方法及系统 | |
CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
CN106157181B (zh) | 一种基于随钻参数对土质围岩快速实时分级的方法 | |
CN116203633B (zh) | 基于天然源震动成像的隧道围岩完整性评价方法 | |
CN112746837A (zh) | 基于分布式光纤传感的页岩油藏勘探数据采集系统及方法 | |
CN115291281B (zh) | 基于深度学习的实时微地震震级计算方法及装置 | |
CN112593922B (zh) | 一种阵列声波测井评价固井二界面胶结质量的方法及装置 | |
CN102073061B (zh) | 使用数字地听仪的地听信息高密度记录系统 | |
CN106908837A (zh) | 一种压裂裂缝形态和裂缝高度确定方法 | |
CN107506556B (zh) | 一种确定新鲜完整岩块声波纵波速度值的方法 | |
CN115327616A (zh) | 一种海量数据驱动的矿山微震震源自动定位方法 | |
CN107703538A (zh) | 地下不良地质勘测数据采集分析系统及方法 | |
CN111852467B (zh) | 一种砂岩铀矿矿体延伸范围的圈定方法及系统 | |
CN115373029B (zh) | 基于深度学习的实时微地震震源机制计算方法及系统 | |
CN106324661A (zh) | 一种利用天然源面波进行岩体风化分带的方法 | |
CN112379412B (zh) | 采集地震数据的质量监控方法及装置 | |
CN116413831A (zh) | 多井合成地震记录自动标定及构造解释和储层预测方法 | |
CN110579806B (zh) | 一种密井网条件下快速的井震标定方法 | |
CN109598049B (zh) | 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法 | |
CN107679614A (zh) | 一种基于粒子群优化的声波时差实时提取方法 | |
CN110780345A (zh) | 一种隧道超前地震探测地震数据的三维速度分析方法 | |
CN105929445B (zh) | 微破裂向量扫描方法 | |
CN112360447B (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 |