CN109471169A - 一种基于地震折射波法的岩体完整性评价方法 - Google Patents

一种基于地震折射波法的岩体完整性评价方法 Download PDF

Info

Publication number
CN109471169A
CN109471169A CN201811320037.7A CN201811320037A CN109471169A CN 109471169 A CN109471169 A CN 109471169A CN 201811320037 A CN201811320037 A CN 201811320037A CN 109471169 A CN109471169 A CN 109471169A
Authority
CN
China
Prior art keywords
velocity
rock mass
longitudinal wave
rock
refraction
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
Application number
CN201811320037.7A
Other languages
English (en)
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.)
China Railway Design Corp
Original Assignee
China Railway Design Corp
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 China Railway Design Corp filed Critical China Railway Design Corp
Priority to CN201811320037.7A priority Critical patent/CN109471169A/zh
Publication of CN109471169A publication Critical patent/CN109471169A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface

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

本发明涉及一种基于地震折射波法的岩体完整性评价方法,包括以下步骤:多炮点的精细化数据采集地震折射波,利用折射层析成像技术获取洞身纵波速度值;新鲜岩体的纵波速度根据规范进行取值;提取洞身纵波速度,结合新鲜岩体纵波速度计算岩体完整性系数Kv;根据隧道里程、岩性、新鲜岩体纵波速度、岩体完整性系数Kv的数据给出岩体完整性评价表;绘制纵波速度断面图,按照洞身地层岩性及速度异常进行分段解释说明纵波速度分布情况、纵波速度范围、地层、岩体完整性、裂隙发育情况,实现岩体完整性评价。本发明能够解决以往地层纵波速度无法连续获取,岩体完整性系数Kv无法连续定量计算等问题,有效提高岩体完整性评价的精度及效率,降低施工风险。

Description

一种基于地震折射波法的岩体完整性评价方法
技术领域
本发明属于隧道岩体完整程度测量的技术领域,涉及一种基于地震折射波法的岩体完整性评价方法。
背景技术
岩体完整性评价是隧道围岩分级的重要指标,隧道围岩分级是隧道设计最重要的依据。不同围岩等级的隧道,工程造价差异极大,围岩等级不准确可能导致项目资金预算缺口使得项目无法顺利进行;不同的围岩等级需要采取不同的施工工艺,围岩等级判断不准确将带来极大的施工风险。
目前岩体完整性评价的测试方法现状为:
1.地震折射波法测试
地震折射波法是获取纵波速度一种较为常用的方法,可以测试土层和岩石的波速。采用地面激发以产生弹性波,地面由检波器接收弹性波。地震折射波法的解释采用常规的表层剥去法及T0法,只能对土石分界进行划分,无法提取洞身高程的纵波速度。
2.岩芯速度测试
需要取芯测试,完整岩石速度准确,破损的岩石速度不准确,不能测土体速度。只能进行单孔测试,无法进行剖面测试。该方法利用钻孔岩芯,成本高、周期长,其评价的主体在钻孔周边围岩,不能连续、完整的进行隧道围岩评价。
3.岩体完整性
现阶段岩体完整性的计算只能通过室内测试岩芯的纵波速度,并结合新鲜岩体的经验值获取,获取值单一,评价简单。
因此,现有技术中的岩体完整性评价方法存在着不同的缺陷,急需一种简便快捷的实用方法来解决这种难题。
发明内容
本发明的目的是提供一种基于地震折射波法的岩体完整性评价方法。本发明利用多炮点的地震波精细化数据采集技术获得丰富的地层信息,采用折射层析成像技术获取地层岩体速度,进行精细成像和网格剖分,提取洞身纵波速度值,根据规范和地区经验给定新鲜岩体纵波速度,计算岩体完整性系数Kv,给出隧道洞身完整、连续的岩体完整性评价,可提高生产效率,确保项目的预算与实施,并极大的减少工程施工风险。
为解决上述问题,本发明采用如下技术方案:
一种基于地震折射波法的岩体完整性评价方法,包括以下步骤:
(1)岩体完整性评价表:
①多炮点的精细化数据采集地震折射波,利用折射层析成像技术获取洞身纵波速度值;
②新鲜岩体的纵波速度根据规范和地区经验进行取值;
③提取洞身纵波速度,结合新鲜岩体纵波速度计算岩体完整性系数Kv;
④根据隧道里程、岩性、新鲜岩体纵波速度、岩体完整性系数Kv的数据给出岩体完整性评价表;
(2)绘制纵波速度断面图,按照洞身地层岩性及速度异常进行分段解释说明纵波速度分布情况、纵波速度范围、地层、岩体完整性、裂隙发育情况,实现岩体的完整性评价。
优选的,多炮点的数目为13~15炮,每20米一炮,具体炮数依据现场采集情况确定,炮点设置在地面上即可。
优选的,步骤①中,对地震折射波法采用RefractionWav程序进行表层剥去法/T0法解释,即采用RefractionWav程序对地震折射资料进行处理,对土石分界进行划分;同时要对时距曲线进行细致划分。
优选的,步骤①中,对每一炮的初至采用Plotrefa程序进行层析成像处理;即采用Plotrefa程序对地震折射资料进行层析成像反演,绘制洞身纵波速度断面图。
优选的,步骤①中,采用SUFER软件提取隧道洞身纵波速度。
优选的,步骤②中,新鲜岩体的纵波速度具体的取值规则为根据《工程岩土学》中公开的规范和地区经验取新鲜岩体对应的范围的最大值,具体见(见表2)
优选的,步骤②中,岩体完整性系数Kv计算方法为Kv=(Vpm/Vpr)*(Vpm/Vpr);其中,Vpm为提取的洞身纵波速度,Vpr为新鲜岩体纵波速度。
优选的,隧道里程实测得到,岩性根据钻孔资料获得
有益效果
本发明能够解决以往地层纵波速度无法连续获取,岩体完整性系数Kv无法连续定量计算等问题,将有效提高岩体完整性评价的精度及效率,降低施工风险,有良好的应用前景和经济效益。
附图说明
图1为某隧道物探解释成果图。
具体实施方式
下面结合具体实施方式,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
实施例1
(1)岩体完整性评价表:
①多炮点的精细化数据采集地震折射波,利用折射层析成像技术获取洞身纵波速度值;常规地震折射波法外业采集需5炮可以达到土石分界的目的,现采用加炮采集(13~15炮),达到足够的射线密度,对地震折射波法采用RefractionWav程序进行表层剥去法/T0法解释;
其中,表层剥去法/T0法:采用RefractionWav程序对地震折射资料进行处理,对土石分界进行划分;同时要对时距曲线进行细致划分;
层析成像:对每一炮的初至采用Plotrefa程序进行层析成像处理,即采用Plotrefa程序对地震折射资料进行层析成像反演,绘制纵波速度断面图;层析成像反演结果网格化后提取洞身顶部速度,采用SUFER软件提取隧道洞身纵波速度,sufer克里格网格化(横纵)10*1;
②新鲜岩体的纵波速度根据规范和地区经验进行取平均速度;
③提取洞身纵波速度,结合新鲜岩体纵波速度计算岩体完整性系数Kv;岩体完整性系数Kv计算方法为Kv=(Vpm/Vpr)*(Vpm/Vpr);其中,Vpm为提取的洞身纵波速度,Vpr为新鲜岩体纵波速度。
④根据隧道里程、岩性、新鲜岩体纵波速度、岩体完整性系数Kv的数据给出岩体完整性评价表;通过岩体完整性系数Kv确定岩体的完整程度,具体的对比关系见表1;通过表2确定新鲜岩体纵波速度;隧道里程实测得到,岩性根据钻孔资料获得。最后根据隧道里程、岩性、新鲜岩体纵波速度、岩体完整性系数Kv的数据给出岩体完整性评价表,具体见表3。
绘制纵波速度断面图,按照洞身地层岩性及速度异常进行分段解释说明纵波速度分布情况、纵波速度范围、地层、岩体完整性、裂隙发育情况,实现岩体的完整性评价,具体见图1,其中,图中数字及等值线为纵波速度,代表地层及岩体的纵波速度,纵波速度分布情况、纵波速度范围从图中可以看出范围区间,或从反演数据列表中提取、地层根据地质专业提供钻孔资料或现场调查得出、岩体完整性为计算的Kv与定性划分的岩体完整程度的对应关系表中对比得出、裂隙发育情况裂隙发育情况根据岩体完整性及图中速度、等值线变化推断解释;下部推断解释为常规物探推断解释:根据图中等值线变化、速度变化区间、计算的岩体完整性系数、洞身埋深等做推断解释。从图1可以看出该段洞身纵波速度呈低俗分布,纵波速度在1.0km/s~1.5km/s范围内,推断洞身在闪长岩(强风化)地层中,穿过,岩体完整性为极破碎,局部存在裂隙,该段洞身埋深较浅,设计施工时应注意塌方、涌水等地质灾害。
表1 Kv与定性划分的岩体完整程度的对应关系(GB50021-2001岩土工程勘察规范)
Kv >0.75 0.75~0.55 0.55~0.35 0.35~0.15 <0.15
完整程度 完整 较完整 较破碎 破碎 极破碎
表2常见岩体的纵波速度(km/s)(据唐大雄等,1987)
表3岩体完整性评价表

Claims (9)

1.一种基于地震折射波法的岩体完整性评价方法,其特征是,包括以下步骤:
(1)岩体完整性评价表:
①多炮点的精细化数据采集地震折射波,利用折射层析成像技术获取洞身纵波速度值;
②新鲜岩体的纵波速度根据规范和地区经验进行取值;
③提取洞身纵波速度,结合新鲜岩体纵波速度计算岩体完整性系数Kv;
④根据隧道里程、岩性、新鲜岩体纵波速度、岩体完整性系数Kv的数据给出岩体完整性评价表;
(2)绘制纵波速度断面图,按照洞身地层岩性及速度异常进行分段解释说明纵波速度分布情况、纵波速度范围、地层、岩体完整性、裂隙发育情况,实现岩体的完整性评价。
2.根据权利要求1所述的基于地震折射波法的岩体完整性评价方法,其特征在于,多炮点设置时的炮点间的距离为20米。
3.根据权利要求2所述的基于地震折射波法的岩体完整性评价方法,其特征在于,多炮点的数目为13~15炮。
4.根据权利要求1所述的基于地震折射波法的岩体完整性评价方法,其特征在于,步骤①中,对地震折射波法采用RefractionWav程序进行表层剥去法/T0法解释,即采用RefractionWav程序对地震折射资料进行处理,对土石分界进行划分;同时要对时距曲线进行细致划分。
5.根据权利要求1所述的基于地震折射波法的岩体完整性评价方法,其特征在于,步骤①中,对每一炮的初至采用Plotrefa程序进行层析成像处理;即采用Plotrefa程序对地震折射资料进行层析成像反演,绘制洞身纵波速度断面图。
6.根据权利要求1所述的基于地震折射波法的岩体完整性评价方法,其特征在于,步骤①中,采用SUFER软件提取隧道洞身纵波速度。
7.根据权利要求1所述的基于地震折射波法的岩体完整性评价方法,其特征在于,步骤②中,新鲜岩体的纵波速度具体的取值规则为根据《工程岩土学》中公开的规范和地区经验取新鲜岩体对应的范围的最大值。
8.根据权利要求1所述的基于地震折射波法的岩体完整性评价方法,其特征在于,步骤②中,岩体完整性系数Kv计算方法为Kv=(Vpm/Vpr)*(Vpm/Vpr);其中,Vpm为提取的洞身纵波速度,Vpr为新鲜岩体纵波速度。
9.根据权利要求1所述的基于地震折射波法的岩体完整性评价方法,其特征在于,步骤②中,隧道里程实测得到,岩性根据钻孔资料获得。
CN201811320037.7A 2018-11-07 2018-11-07 一种基于地震折射波法的岩体完整性评价方法 Pending CN109471169A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811320037.7A CN109471169A (zh) 2018-11-07 2018-11-07 一种基于地震折射波法的岩体完整性评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811320037.7A CN109471169A (zh) 2018-11-07 2018-11-07 一种基于地震折射波法的岩体完整性评价方法

Publications (1)

Publication Number Publication Date
CN109471169A true CN109471169A (zh) 2019-03-15

Family

ID=65671907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811320037.7A Pending CN109471169A (zh) 2018-11-07 2018-11-07 一种基于地震折射波法的岩体完整性评价方法

Country Status (1)

Country Link
CN (1) CN109471169A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110794039A (zh) * 2019-11-06 2020-02-14 长江勘测规划设计研究有限责任公司 利用岩体波速计算帷幕灌浆岩体裂隙充填率的方法
CN112923879A (zh) * 2021-02-04 2021-06-08 中水北方勘测设计研究有限责任公司 一种利用声波速度评价洞室围岩松弛厚度及等级的方法
CN113010942A (zh) * 2021-02-25 2021-06-22 中国铁路设计集团有限公司 隧道开挖风险预警及围岩分级评价方法
CN114355447A (zh) * 2022-01-13 2022-04-15 中国安能集团第三工程局有限公司 一种隧道工程围岩等级快速划分方法及装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102087367A (zh) * 2010-11-19 2011-06-08 中冶集团武汉勘察研究院有限公司 纵波-面波速度比值法划分基岩风化带的方法
CN102435295A (zh) * 2011-09-07 2012-05-02 中铁第四勘察设计院集团有限公司 隧道围岩岩体弹性波速度的探测方法
CN102680577A (zh) * 2010-12-23 2012-09-19 中冶集团武汉勘察研究院有限公司 工程围岩完整性系数Kv的地应力修正取值方法
CN106290576A (zh) * 2016-08-12 2017-01-04 中国电建集团成都勘测设计研究院有限公司 一种复合型蚀变岩/带的工程地质分类方法
CN107091623A (zh) * 2017-03-20 2017-08-25 西安工业大学 隧道围岩松动圈厚度计算方法
CN107506556A (zh) * 2017-08-31 2017-12-22 中水北方勘测设计研究有限责任公司 一种确定新鲜完整岩块声波纵波速度值的简便方法
CN107741488A (zh) * 2017-09-20 2018-02-27 西安工业大学 隧道施工期岩体完整性系数Kv的计算方法
CN108267778A (zh) * 2018-01-04 2018-07-10 中国铁路设计集团有限公司 一种地层速度测试系统及测试方法
CN108562648A (zh) * 2018-05-24 2018-09-21 山东大学 用于判断破碎岩体完整度的装置及方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102087367A (zh) * 2010-11-19 2011-06-08 中冶集团武汉勘察研究院有限公司 纵波-面波速度比值法划分基岩风化带的方法
CN102680577A (zh) * 2010-12-23 2012-09-19 中冶集团武汉勘察研究院有限公司 工程围岩完整性系数Kv的地应力修正取值方法
CN102435295A (zh) * 2011-09-07 2012-05-02 中铁第四勘察设计院集团有限公司 隧道围岩岩体弹性波速度的探测方法
CN106290576A (zh) * 2016-08-12 2017-01-04 中国电建集团成都勘测设计研究院有限公司 一种复合型蚀变岩/带的工程地质分类方法
CN107091623A (zh) * 2017-03-20 2017-08-25 西安工业大学 隧道围岩松动圈厚度计算方法
CN107506556A (zh) * 2017-08-31 2017-12-22 中水北方勘测设计研究有限责任公司 一种确定新鲜完整岩块声波纵波速度值的简便方法
CN107741488A (zh) * 2017-09-20 2018-02-27 西安工业大学 隧道施工期岩体完整性系数Kv的计算方法
CN108267778A (zh) * 2018-01-04 2018-07-10 中国铁路设计集团有限公司 一种地层速度测试系统及测试方法
CN108562648A (zh) * 2018-05-24 2018-09-21 山东大学 用于判断破碎岩体完整度的装置及方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
交通部人事劳动司 编: "《交通勘测 (四)》", 31 December 1994, 人民交通出版社 *
吕宝辉 等: ""综合物探技术在某大型隧道岩土工程勘察及地址灾害勘察中的应用"", 《工程勘察》 *
孙毅 等: ""测定围岩松动圈的折射波法探讨"", 《铁道建筑技术》 *
王奇 等: ""浅层地震折射波在某隧道勘察中的应用"", 《内江科技》 *
许广春 等: ""中浅埋隧道综合物探勘察模式研究"", 《铁道工程学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110794039A (zh) * 2019-11-06 2020-02-14 长江勘测规划设计研究有限责任公司 利用岩体波速计算帷幕灌浆岩体裂隙充填率的方法
CN110794039B (zh) * 2019-11-06 2021-12-21 长江勘测规划设计研究有限责任公司 利用岩体波速计算帷幕灌浆岩体裂隙充填率的方法
CN112923879A (zh) * 2021-02-04 2021-06-08 中水北方勘测设计研究有限责任公司 一种利用声波速度评价洞室围岩松弛厚度及等级的方法
CN112923879B (zh) * 2021-02-04 2022-06-21 中水北方勘测设计研究有限责任公司 一种利用声波速度评价洞室围岩松弛厚度及等级的方法
CN113010942A (zh) * 2021-02-25 2021-06-22 中国铁路设计集团有限公司 隧道开挖风险预警及围岩分级评价方法
CN114355447A (zh) * 2022-01-13 2022-04-15 中国安能集团第三工程局有限公司 一种隧道工程围岩等级快速划分方法及装置

Similar Documents

Publication Publication Date Title
CN109471169A (zh) 一种基于地震折射波法的岩体完整性评价方法
CN104133245B (zh) 一种地震资料的静校正方法及系统
CN105277982B (zh) 一种泥页岩总有机碳含量地震预测方法
CN105607146B (zh) 一种曲流河砂体规模的定量表征方法
Seers et al. Comparison of digital outcrop and conventional data collection approaches for the characterization of naturally fractured reservoir analogues
CN103088803B (zh) 鉴别岩溶塌陷致塌因素的方法
CN103728659A (zh) 一种提高地下岩溶探测精度的方法
CN109884710B (zh) 针对激发井深设计的微测井层析成像方法
Putz-Perrier et al. Distribution of faults and extensional strain in fractured carbonates of the North Malta Graben
Choi et al. Variation of rock quality designation (RQD) with scanline orientation and length: a case study in Korea
Johnson et al. Statistical comparison of methods for estimating sediment thickness from horizontal-to-vertical spectral ratio (HVSR) seismic methods: An example from Tylerville, Connecticut, USA
CN110886604A (zh) 一种基于计算机模拟技术的高效地热资源勘察方法
KR20180055518A (ko) 3차원 잠재광산 부존지 예측 시스템 및 방법
CN103678778A (zh) 一种放射性物化探信息综合方法
CN106650192A (zh) 一种火山岩型铀矿床磁性界面反演方法
CN109799540A (zh) 基于地质信息约束下的火山岩型铀矿床磁化率反演方法
CN102928871B (zh) 一种基于旋转菱形的属性提取与断层描述方法
CN105487106B (zh) 一种基于高斯射线束目的层能量照明的补炮方法
Berrino et al. Repeated absolute gravity measurements on a dense network at Campi Flegrei–a reliable tool for volcano monitoring
CN104280774A (zh) 一种单频地震散射噪声的定量分析方法
CN105259577B (zh) 一种确定地层界面的角度信息的方法及装置
CN104076391A (zh) 基于tti介质四阶旅行时方程的局部角度域各向异性偏移方法
Yu et al. Research on site classification method based on BP neural network
Smith et al. Near-surface seismic imaging of tunnels using 3D elastic full-waveform inversion
CN107765299B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190315

RJ01 Rejection of invention patent application after publication