CN115640995A - 一种隧洞超前地质预报风险评价方法及装置 - Google Patents
一种隧洞超前地质预报风险评价方法及装置 Download PDFInfo
- Publication number
- CN115640995A CN115640995A CN202211212264.4A CN202211212264A CN115640995A CN 115640995 A CN115640995 A CN 115640995A CN 202211212264 A CN202211212264 A CN 202211212264A CN 115640995 A CN115640995 A CN 115640995A
- Authority
- CN
- China
- Prior art keywords
- risk
- geological
- physical parameter
- tunnel
- distribution signal
- 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
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种隧洞超前地质预报风险评价方法及装置,属于隧洞超前地质预报技术领域,所述方法包括:获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。本发明提供的隧洞超前地质预报风险评价方法,操作简单并且能够实现对隧洞地质风险的快速定量评价。
Description
技术领域
本发明涉及隧洞超前地质预报技术领域,尤其涉及一种隧洞超前地质预报风险评价方法及装置。
背景技术
如今,随着隧洞工程的大量建设,地下工程面临的地质条件越来越复杂,因此,隧洞超前地质预报技术被广泛地应用于隧洞开挖过程中,作为前期地质勘察的补充,用于提前探明掌子面前方一定范围内的不良地质体,为隧洞施工过程中提前发现地质风险从而针对性进行风险处置提供技术依据,能够有效指导隧洞施工,保障隧洞施工安全。
目前国内外常用的隧洞超前地质预报技术方法种类较多,根据技术原理,超前地质预报技术可以主要分为地震波法、电磁法、电法等等,每种技术方法的特点与适用条件都不尽相同,由于方法技术的局限性、地质条件的复杂性、以及物探反演的多解性,利用单一的预报方法难以实现对地质风险的精确预报与识别,更加难以实现对地质风险进行分级评价。
虽然综合超前地质预报手段越来越被广泛的应用于隧洞工程中,但是目前行业内仅限于利用不同的物探预报手段相互验证、相互补充去发现地质异常,对于发现的地质异常到底能造成什么不良后果,造成的不良后果严重程度能够到什么级别,并没有建立一套切实有效的隧洞超前地质预报风险评价方法。
近年来,国内外有些学者也提出了一些风险评价的方法,但是大多数方法仅能实现地质风险的定性评价,并未真正通过算法或定量参数实现地质风险的定量风险评价;有些学者提出的评价方法需要采用的参数过多,包括岩石力学参数、物理参数、设计参数、支护参数、监测参数等等,虽然参数多有利于风险评价更全面,但是其中很多参数不是仅仅通过常规超前地质预报手段能够获取的,还需要进行室内试验以及多方配合,需要大量的参数搜集、试验与分析时间,与超前地质预报过程强调的时效性相违背,实际操作过程中不切实际,操作难度极大。
发明内容
本发明提供一种隧洞超前地质预报风险评价方法、装置、电子设备及存储介质,用以解决现有技术中的缺陷,实现了高效的对隧洞地质风险的定量风险评价。
第一方面,本发明提供一种隧洞超前地质预报风险评价方法,包括:获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
根据本发明提供的一种隧洞超前地质预报风险评价方法,根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号,包括:在所述任一目标物理参数与围岩质量呈正相关的情况下,根据所述分布信号与正相关参考平均值的比值,确定所述归一化分布信号;在所述任一目标物理参数与围岩质量呈负相关的情况下,对所述分布信号进行倒数处理,获取中间分布信号;根据所述中间分布信号与正相关参考平均值的比值,确定所述归一化分布信号。
根据本发明提供的一种隧洞超前地质预报风险评价方法,在根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号之前,还包括:在所述任一目标物理参数与围岩质量呈正相关的情况下,将所述任一目标物理参数的参考平均值作为所述正相关参考平均值;在所述任一目标物理参数与围岩质量呈负相关的情况下,将所述任一目标物理参数的参考平均值的倒数作为所述正相关参考平均值。
根据本发明提供的一种隧洞超前地质预报风险评价方法,在目标物理参数为弹性物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:利用地震波法,获取在掌子面前方的预设范围内围岩的弹性物理参数的分布信号;所述弹性物理参数至少包括以下物理参数中的一种:纵波波速、横波波速、泊松比、以及密度。
根据本发明提供的一种隧洞超前地质预报风险评价方法,在目标物理参数为电性物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:在所述电性物理参数为电阻率的情况下,利用瞬变电磁法获取所述电阻率的分布信号;在所述电性物理参数为电阻率和/或极化率的情况下,利用激发极化法获取所述电性物理参数的分布信号。
根据本发明提供的一种隧洞超前地质预报风险评价方法,在目标物理参数为电磁响应物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:利用地质雷达法,获取在掌子面前方的预设范围内围岩的电磁响应物理参数的分布信号;所述电磁响应物理参数至少包括以下物理参数中的一种:磁波振幅、电磁波频率以及电磁波同相轴形态。
根据本发明提供的一种隧洞超前地质预报风险评价方法,根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况,包括:根据所述综合权重评分值的分布状况和预设评分阈值,确定所述预设范围内的地质风险潜在区域;根据所述地质风险潜在区域内的综合权重评分值的分布状况以及预先设置的定量化风险评价与风险等级划分标准,确定所述地质风险潜在区域内的地质风险区域,以及所述地质风险区域的定量风险评价与风险等级;所述地质风险区域的综合权重评分平均值小于预设评分阈值,所述地质风险区域的覆盖范围大于预设范围阈值;所述综合权重评分平均值为所述地质风险区域内的不同位置区域的综合权重评分值的平均值。
第二方面,本发明还提供一种隧洞超前地质预报风险评价装置,包括:第一模块,第二模块,第三模块,第四模块。
第一模块,用于获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;
第二模块,用于根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;
第三模块,用于根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;
第四模块,用于根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
第三方面,本发明提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述任一种所述隧洞超前地质预报风险评价方法的步骤。
第四方面,本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如上述任一种所述隧洞超前地质预报风险评价方法的步骤。
本发明提供的隧洞超前地质预报风险评价方法及装置,通过获取隧洞的围岩的对应的目标物理参数的分布信号,并且基于加权分析计算预设范围内的不同位置区域的综合权重评分值,根据综合权重评分值的分布情况,实现了对隧洞地质风险状况的评价,该方法操作简单并且能够实现对隧洞地质风险的快速定量评价。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的隧洞超前地质预报风险评价方法的流程示意图之一;
图2是本发明提供的隧洞超前地质预报风险评价方法的流程示意图之二;
图3是本发明提供的综合权重评分值直方图;
图4是本发明提供的隧洞超前地质预报风险评价装置的结构示意图;
图5是本发明提供的电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,在本发明实施例的描述中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
本申请中的术语“第一”、“第二”等是用于区别类似的对象,而不用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施,且“第一”、“第二”等所区分的对象通常为一类,并不限定对象的个数,例如第一对象可以是一个,也可以是多个。此外,“和/或”表示所连接对象的至少其中之一,字符“/”,一般表示前后关联对象是一种“或”的关系。
下面结合图1-图5描述本发明实施例所提供的隧洞超前地质预报风险评价方法和装置。
图1是本发明提供的隧洞超前地质预报风险评价方法的流程示意图之一,如图1所示,包括但不限于以下步骤:
步骤101:获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号。
其中,所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成。预设范围内可以为掌子面的前方N米内。
可选地,根据目标地质灾害类型及其对应不良地质体的主要物性特征,确定用于目标地质灾害风险评价的主要的目标物理参数。
其中,所述不良地质体,主要包括断层、破碎带、软弱夹层、岩溶、地下水等,每种不良地质体的物性特征与结构形态都不尽相同,从而表现在各种不良地质体在某几种物理参数上表现最为突出。
目标物理参数的种类包括但不限于:纵波波速(Vp)、横波波速(Vs)、泊松比(V)、密度(ρ)、电阻率(ρs)、极化率(η)、电磁波振幅(A)、电磁波频率(f)、电磁波同相轴形态(S)等。
步骤102:根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号。
在所述任一目标物理参数与围岩质量呈正相关的情况下,根据所述分布信号与正相关参考平均值的比值,确定所述归一化分布信号;
在所述任一目标物理参数与围岩质量呈负相关的情况下,对所述分布信号进行倒数处理,获取中间分布信号;根据所述中间分布信号与正相关参考平均值的比值,确定所述归一化分布信号。
相应的归一化分布信号的表达式如下:
步骤103:根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值。
可选地,针对不同的地质灾害类型,采用加权分析法,确定各项目标物理参数对应的加权系数WP。
所述加权系数WP,指目标物理参数P对应的加权系数,代表该目标物理参数在可能导致对应地质灾害中的重要程度,赋值范围为0到1。
式中,k为参与风险评价的目标物理参数的数量,Pj为第j个目标物理参数P。
相应地,根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值的具体公式为:
Qi∈(0,1]
Qi代表掌子面前方第i米范围内(即不同位置区域)围岩的综合权重评分值,为常数,无单位,其数值越小代表围岩存在地质风险的可能性越大,反之则越小。
步骤104:根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
可选地,根据在预设范围内的不同区域的综合权重评分值Qi的分布情况,对预设范围内的隧洞地质风险状况进行评价。可选地,综合权重评分值越低的区域代表存在地质风险的可能性越大,综合权重评分值越高的区域代表存在地质风险的可能性越小。从而,本发明可以根据不同区域的综合权重评分值,对预设范围内的地质风险进行评价。
需要说明的是,为了可以更加方便的执行本发明的方法,在实施上述步骤之前,还可以预先建立本发明中的地质风险评价体系。具体方式如下。
根据前期勘察和施工地质资料,初步判断隧洞在施工过程中可能面临的主要地质灾害类型,针对目标地质灾害,建立对应的风险评价体系,并确定风险评价与风险等级划分标准。
所述主要地质灾害类型,主要包括:涌水突泥,破碎岩体塌方、软岩大变形、岩溶空洞、岩溶涌水等等。
表1是风险评价与风险等级划分标准表,下面结合表1对风险评价与风险等级划分标准进行说明。
对于风险等级,可根据工程需求分为N个级别,本发明以分为5个风险等级为例进行说明,风险等级分为A、B、C、D、E,其中A代表风险等级最高,E代表风险等级最低。
可选地,本发明从风险程度和风险规模两个方面进行风险等级的划分。
其中,所述风险程度代表不良地质体可能引发的地质灾害程度,与不良地质体的目标物理参数数值相关,目标物理参数对应的围岩条件越差,其对应的风险程度越高。
所述风险规模代表不良地质体在隧洞轴线范围内的尺寸,在相同的风险程度情况下,风险规模越大,其对应的风险等级越高。
所述风险规模临界值划分不同风险等级规模的阈值。
表1风险评价与风险等级划分标准表
本发明提供的隧洞超前地质预报风险评价方法,通过获取隧洞的围岩的对应的目标物理参数的分布信号,并且基于加权分析计算预设范围内的不同位置区域的综合权重评分值,根据综合权重评分值的分布情况,实现了对隧洞地质风险状况的评价,该方法操作简单并且能够实现对隧洞地质风险的快速定量评价。
基于上述实施例的内容,作为一种可选的实施例,本发明提供的隧洞超前地质预报风险评价方法,针对不同岩性的围岩,确定各项目标物理参数P在该岩性围岩质量较好情况下的参考平均值并对所有用于风险评价的目标物理参数进行正相关处理,形成正相关参考平均值使目标物理参数的数值与围岩质量RQ的好坏成正相关关系,并对隧洞模型一定范围内的围岩进行正相关参考平均值赋值。
具体地,正相关处理是指在所述任一目标物理参数与围岩质量呈正相关的情况下,将所述任一目标物理参数的参考平均值作为所述正相关参考平均值;在所述任一目标物理参数与围岩质量呈负相关的情况下,将所述任一目标物理参数的参考平均值的倒数作为所述正相关参考平均值。
可选地,当该物理参数数值大于正相关参考平均值时,可以认为围岩在该物理参数方面为较好的围岩,且数值越大,代表围岩质量在该方面越好;反之,当该物理参数数值小于正相关参考平均值时,可以认为围岩在该物理参数方面为较差的围岩,且数值越小,代表围岩质量在该方面越差。具体公式为:
基于上述实施例的内容,作为一种可选的实施例,本发明提供一种隧洞超前地质预报风险评价方法,在目标物理参数为弹性物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:利用地震波法,获取在掌子面前方的预设范围内围岩的弹性物理参数的分布信号;所述弹性物理参数至少包括以下物理参数中的一种:纵波波速、横波波速、泊松比、以及密度。
可选地,采用地震波法开展超前地质预报工作,获取掌子面前方N米范围内围岩的纵波波速(Vp)、横波波速(Vs)、泊松比(V)、以及密度(ρ)等物理参数的分布信号。
所述地震波法,包括目前国内外主流的地震波法预报系统,所述地震波法,包括目前国内外主流的地震波法预报系统,例如瑞士研发生产的TSP系统、俄罗斯研发生产的TGS系统、美国研发生产的TRT系统、中国研发生产的TGP系统等。
Vpi代表掌子面前方第i米范围内围岩的纵波波速,单位为m/s,其数值越大代表围岩在该目标物理参数方面越好,反之则越差。
{Vpi}={Vp1,Vp2,……,VpN},其中i=1,2,……,N。
Vsi代表掌子面前方第i米范围内围岩的横波波速,单位为m/s,其数值越大代表围岩在该目标物理参数方面越好,反之则越差。
{Vsi}={Vs1,Vs2,……,VsN},其中i=1,2,……,N。
Vi代表掌子面前方第i米范围内围岩的泊松比,为常数,无单位,其数值越大代表围岩在该目标物理参数方面越差,反之则越好。
{Vi}={V1,V2,……,VN},其中i=1,2,……,N。
ρi代表掌子面前方第i米范围内围岩的密度,单位为g/cm3,其数值越大代表围岩在该目标物理参数方面越好,反之则越差。
{ρi}={ρ1,ρ2,……,ρN},其中i=1,2,……,N。
基于上述实施例的内容,作为一种可选的实施例,本发明提供一种隧洞超前地质预报风险评价方法,在目标物理参数为电性物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:在所述电性物理参数为电阻率的情况下,利用瞬变电磁法获取所述电阻率的分布信号;在所述电性物理参数为电阻率和/或极化率的情况下,利用激发极化法获取所述电性物理参数的分布信号。
ηi代表掌子面前方第i米范围内围岩的极化率,单位为%,其数值越大代表围岩含水率越高,反之则越低。
{ηi}={η1,η2,……,ηN},其中i=1,2,……,N。
基于上述实施例的内容,作为一种可选的实施例,本发明提供一种隧洞超前地质预报风险评价方法,在目标物理参数为电磁响应物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:利用地质雷达法,获取在掌子面前方的预设范围内围岩的电磁响应物理参数的分布信号;所述电磁响应物理参数至少包括以下物理参数中的一种:磁波振幅、电磁波频率以及电磁波同相轴形态。
采用地质雷达法开展超前地质预报工作,获取掌子面前方N米范围内围岩的电磁波振幅(A)、电磁波频率(f)、电磁波同相轴形态(S)等物理参数的分布信号。
Ai代表掌子面前方第i米范围内围岩的电磁波振幅,单位为V,其数值越大代表围岩完整性越差,反之则越好。
{Ai}={A1,A2,……,AN},其中i=1,2,……,N。
fi代表掌子面前方第i米范围内围岩的电磁波频率,单位为MHz,其数值越大代表围岩质量越好,反之则越差。
{fi}={f1,f2,……,fN},其中i=1,2,……,N。
Si代表掌子面前方第i米范围内围岩的电磁波同相轴形态,例如水平状、倾斜状、抛物线状、交叉状、错乱状等。当围岩质量较好、岩质较均匀时,一般呈“水平状”,其它形态代表着岩体构造或者不良地质体的形态分布特征。
基于上述实施例的内容,作为一种可选的实施例,本发明提供一种隧洞超前地质预报风险评价方法,在所述任一目标物理参数与围岩质量呈正相关的情况下,根据所述分布信号与正相关参考平均值的比值,确定所述归一化分布信号;在所述任一目标物理参数与围岩质量呈负相关的情况下,对所述分布信号进行倒数处理,获取中间分布信号;根据所述中间分布信号与正相关参考平均值的比值,确定所述归一化分布信号。
当MP∝RQ时,
进一步的,根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值的具体公式为:
Qi∈(0,1]
Qi代表掌子面前方第i米范围内(即不同位置区域)围岩的综合权重评分值,为常数,无单位,其数值越小代表围岩存在地质风险的可能性越大,反之则越小。
基于上述实施例的内容,作为一种可选的实施例,本发明提供一种隧洞超前地质预报风险评价方法,根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况,包括:根据所述综合权重评分值的分布状况和预设评分阈值,确定所述预设范围内的地质风险潜在区域;根据所述地质风险潜在区域内的综合权重评分值的分布状况以及预先设置的定量化风险评价与风险等级划分标准,确定所述地质风险潜在区域内的地质风险区域,以及所述地质风险区域的定量风险评价与风险等级;所述地质风险区域的综合权重评分平均值小于预设评分阈值,所述地质风险区域的覆盖范围大于预设范围阈值;所述综合权重评分平均值为所述地质风险区域内的不同位置区域的综合权重评分值的平均值。
可选地,绘制超前预报风险评价范围内的综合权重评分值Qi直方图,与基于正相关物理平均值建立的经验模型对应的综合权重评分值直方图进行比对分析,提取综合权重评分值偏低的区域,即综合权重评分值Qi<1的区域,于是得到地质风险潜在区域。其中,预设评分阈值为1。
假设第m米至m+n米的综合权重评分值Qi均小于1,形成地质风险潜在区域对应的数列{LQl}。
{LQl}={Qm,Qm+1,……,Qm+n},其中l=1,2,……,n+1。
进一步地,制定定量化风险评价与风险等级划分标准,如表2所示出定量化风险评价与风险等级划分标准。
表2定量化风险评价与风险等级划分标准
表2中,LQX代表综合权重评分平均值对应X风险等级的临界值;FCX代表地质风险潜在区域对应X风险等级的临界值。X可以取A、B、C、D、E。
参照预先建立的风险评价与风险等级划分标准,引入综合权重评分平均值概念,对风险程度进行定量划分;并在地质风险潜在区域中筛选能够满足更高风险等级的地质风险区域,地质风险区域既能满足综合权重评分平均值的更高要求,也能满足地质风险区域规模的定量要求。
进一步地,本发明可以参照已建立的定量化风险评价与风险等级划分标准,综合考虑“综合权重评分平均值”和“地质风险区域”两个重要参数,确定隧洞超前地质预报的风险区域及其对应的风险等级。
为了更好的对本发明提供的隧洞超前地质预报风险评价方法进行说明,下面以一个完整的实施例对本发明的技术方案进行说明。
某省某引调水工程,在某隧洞掌子面开展超前地质预报工作时,掌子面岩性为灰岩,以弱风化为主,岩质坚硬,岩体较完整到完整性差,岩体自稳能力总体较好,地下水状态以渗水为主,综合判定围岩类别为Ⅲ类,此时掌子面揭露的围岩质量总体较好,从掌子面揭露的围岩情况看,目前存在地质风险的可能性较小。
但是,由于隧洞沿线地质条件极为复杂,特别表现在岩溶涌水问题,不能保证开挖过程中一直不存在风险,本次依然开展了综合物探超前地质预报工作,对前方30m范围内的围岩地质情况进行了地质风险评价与风险等级划分。
图2是本发明提供的隧洞超前地质预报风险评价方法的流程示意图之二,如图2所示,步骤如下:
步骤201:建立地质风险评价体系。
由于该隧洞工程沿线穿越强岩溶区域,且地下水发育,且施工过程中多次揭露溶洞、溶隙等岩溶结构,分析判断该隧洞在施工过程中面临的主要地质灾害类型为岩溶涌水突泥问题,针对该地质问题,建立涌水突泥风险评价体系,并确定风险等级的划分标准,如表3所示。
步骤202:确定用于风险评价的目标物理参数。
由于本次风险评价主要针对岩溶涌水风险,对应的不良地质体主要包括岩溶和地下水,根据岩溶和地下水的物性特征,确定本次参与风险评价的主要目标物理参数包括:横波波速(Vs)、泊松比(V)、电阻率(ρs)、极化率(η)、电磁波频率(f)。
表3风险等级划分标准表
步骤203:建立正相关物理参数经验模型。
首先,根据本隧洞的岩性灰岩,确定参与本次风险分析的物理参数横波波速(Vs)、泊松比(V)、电阻率(ρs)、极化率(η)、电磁波频率(f)在围岩质量较好情况下的参考平均值表4示出了目标物理参数平均值统计表。
表4目标物理参数参考平均值统计表
其中,横波波速(Vs)、电阻率(ρs)、电磁波频率(f)已与围岩质量RQ的好坏成正相关关系,因此:
其中,泊松比(V)、极化率(η)与围岩质量RQ的好坏成负相关关系,需要进行正相关处理,处理后:
正相关参考平均值的计算结果,可以参见表5。
表5目标物理参数的正相关参考平均值统计表
步骤204:采用地震波法预报获取相关弹性物理参数。
首先采用目前国内外主流的TSP地震波法开展超前地质预报工作,获取掌子面前方100m范围内围岩的纵波波速(Vp)、横波波速(Vs)、泊松比(V)、以及密度(ρ)等物理参数的分布情况。
由于本次风险评价主要针对隧洞前方30m范围内的围岩,并根据需要参与风险评价的目标物理参数,提取前30m范围内横波波速(Vs)、泊松比(V)的数值,形成对应的数列{Vsi}和{Vi},其中i=1,2,……,30。
步骤205:采用瞬变电磁法或激发极化法预报获取相关电性物理参数。
步骤206:采用地质雷达法预报获取相关电磁响应物理参数。
采用地质雷达法开展超前地质预报工作,获取掌子面前方30m范围内围岩的电磁波振幅(A)、电磁波频率(f)、电磁波同相轴形态(S)等物理参数的分布情况,并根据需要参与风险评价的物理参数,形成电磁波频率(f)物理参数对应的数列{fi},其中i=1,2,……,30。
步骤207:将实测值与参考平均值进行归一化比值计算。
其中,横波波速(Vs)、电阻率(ρs)、电磁波频率(f)已与围岩质量RQ的好坏成正相关关系,满足以下算法:
统计计算可得,表6所示的横波波速归一化比值计算表;表7所示的电阻率归一化比值计算表;表8所示的电磁波频率归一化比值计算表。
其中,泊松比(V)、极化率(η)与围岩质量RQ的好坏成负相关关系,满足以下算法:
统计计算可得,表9所示的泊松比归一化比值计算表,表10所示的极化率归一化比值计算表。
表6横波波速归一化比值计算表
第i米 | 实测值(m/s) | 正相关参考平均值(m/s) | 归一化比值 |
1 | 2750 | 2400 | 1.00 |
2 | 2720 | 2400 | 1.00 |
3 | 2700 | 2400 | 1.00 |
4 | 2740 | 2400 | 1.00 |
5 | 2690 | 2400 | 1.00 |
6 | 2700 | 2400 | 1.00 |
7 | 2630 | 2400 | 1.00 |
8 | 2510 | 2400 | 1.00 |
9 | 2440 | 2400 | 1.00 |
10 | 2420 | 2400 | 1.00 |
11 | 2220 | 2400 | 0.93 |
12 | 2110 | 2400 | 0.88 |
13 | 2100 | 2400 | 0.88 |
14 | 2130 | 2400 | 0.89 |
15 | 2050 | 2400 | 0.85 |
16 | 2000 | 2400 | 0.83 |
17 | 1950 | 2400 | 0.81 |
18 | 1970 | 2400 | 0.82 |
19 | 2040 | 2400 | 0.85 |
20 | 2090 | 2400 | 0.87 |
21 | 2140 | 2400 | 0.89 |
22 | 2200 | 2400 | 0.92 |
23 | 2400 | 2400 | 1.00 |
24 | 2470 | 2400 | 1.00 |
25 | 2690 | 2400 | 1.00 |
26 | 2640 | 2400 | 1.00 |
27 | 2750 | 2400 | 1.00 |
28 | 2800 | 2400 | 1.00 |
29 | 2780 | 2400 | 1.00 |
30 | 2770 | 2400 | 1.00 |
表7电阻率归一化比值计算表
第i米 | 实测值(Ω·m) | 正相关参考平均值(Ω·m) | 归一化比值 |
1 | 700 | 500 | 1.00 |
2 | 720 | 500 | 1.00 |
3 | 700 | 500 | 1.00 |
4 | 650 | 500 | 1.00 |
5 | 620 | 500 | 1.00 |
6 | 510 | 500 | 1.00 |
7 | 480 | 500 | 0.96 |
8 | 460 | 500 | 0.92 |
9 | 440 | 500 | 0.88 |
10 | 440 | 500 | 0.88 |
11 | 310 | 500 | 0.62 |
12 | 180 | 500 | 0.36 |
13 | 190 | 500 | 0.38 |
14 | 150 | 500 | 0.30 |
15 | 130 | 500 | 0.26 |
16 | 120 | 500 | 0.24 |
17 | 100 | 500 | 0.20 |
18 | 90 | 500 | 0.18 |
19 | 100 | 500 | 0.20 |
20 | 140 | 500 | 0.28 |
21 | 190 | 500 | 0.38 |
22 | 260 | 500 | 0.52 |
23 | 390 | 500 | 0.78 |
24 | 460 | 500 | 0.92 |
25 | 550 | 500 | 1.00 |
26 | 530 | 500 | 1.00 |
27 | 640 | 500 | 1.00 |
28 | 690 | 500 | 1.00 |
29 | 710 | 500 | 1.00 |
30 | 740 | 500 | 1.00 |
表8电磁波频率归一化比值计算表
第i米 | 实测值(MHz) | 正相关参考平均值(MHz) | 归一化比值 |
1 | 95 | 80 | 1.00 |
2 | 95 | 80 | 1.00 |
3 | 93 | 80 | 1.00 |
4 | 94 | 80 | 1.00 |
5 | 90 | 80 | 1.00 |
6 | 91 | 80 | 1.00 |
7 | 88 | 80 | 1.00 |
8 | 85 | 80 | 1.00 |
9 | 75 | 80 | 0.94 |
10 | 69 | 80 | 0.86 |
11 | 46 | 80 | 0.58 |
12 | 44 | 80 | 0.55 |
13 | 43 | 80 | 0.54 |
14 | 39 | 80 | 0.49 |
15 | 34 | 80 | 0.43 |
16 | 32 | 80 | 0.40 |
17 | 30 | 80 | 0.38 |
18 | 31 | 80 | 0.39 |
19 | 39 | 80 | 0.49 |
20 | 44 | 80 | 0.55 |
21 | 48 | 80 | 0.60 |
22 | 62 | 80 | 0.78 |
23 | 70 | 80 | 0.88 |
24 | 79 | 80 | 0.99 |
25 | 80 | 80 | 1.00 |
26 | 80 | 80 | 1.00 |
27 | 84 | 80 | 1.00 |
28 | 89 | 80 | 1.00 |
29 | 94 | 80 | 1.00 |
30 | 97 | 80 | 1.00 |
表9泊松比归一化比值计算表
第i米 | 实测值 | 正相关参考平均值 | 归一化比值 |
1 | 0.26 | 3.33 | 1.00 |
2 | 0.26 | 3.33 | 1.00 |
3 | 0.26 | 3.33 | 1.00 |
4 | 0.26 | 3.33 | 1.00 |
5 | 0.27 | 3.33 | 1.00 |
6 | 0.26 | 3.33 | 1.00 |
7 | 0.28 | 3.33 | 1.00 |
8 | 0.28 | 3.33 | 1.00 |
9 | 0.29 | 3.33 | 1.00 |
10 | 0.3 | 3.33 | 1.00 |
11 | 0.34 | 3.33 | 0.88 |
12 | 0.36 | 3.33 | 0.83 |
13 | 0.36 | 3.33 | 0.83 |
14 | 0.37 | 3.33 | 0.81 |
15 | 0.37 | 3.33 | 0.81 |
16 | 0.38 | 3.33 | 0.79 |
17 | 0.39 | 3.33 | 0.77 |
18 | 0.39 | 3.33 | 0.77 |
19 | 0.37 | 3.33 | 0.81 |
20 | 0.36 | 3.33 | 0.83 |
21 | 0.35 | 3.33 | 0.86 |
22 | 0.34 | 3.33 | 0.88 |
23 | 0.29 | 3.33 | 1.00 |
24 | 0.28 | 3.33 | 1.00 |
25 | 0.26 | 3.33 | 1.00 |
26 | 0.26 | 3.33 | 1.00 |
27 | 0.25 | 3.33 | 1.00 |
28 | 0.24 | 3.33 | 1.00 |
29 | 0.25 | 3.33 | 1.00 |
30 | 0.25 | 3.33 | 1.00 |
表10极化率归一化比值计算表
第i米 | 实测值 | 正相关参考平均值 | 归一化比值 |
1 | 1.50% | 50 | 1.00 |
2 | 1.30% | 50 | 1.00 |
3 | 1.30% | 50 | 1.00 |
4 | 1.40% | 50 | 1.00 |
5 | 1.50% | 50 | 1.00 |
6 | 1.60% | 50 | 1.00 |
7 | 1.70% | 50 | 1.00 |
8 | 1.80% | 50 | 1.00 |
9 | 2.20% | 50 | 0.91 |
10 | 2.60% | 50 | 0.77 |
11 | 4.50% | 50 | 0.44 |
12 | 6.20% | 50 | 0.32 |
13 | 6.70% | 50 | 0.30 |
14 | 7.80% | 50 | 0.26 |
15 | 8.80% | 50 | 0.23 |
16 | 9.20% | 50 | 0.22 |
17 | 10.70% | 50 | 0.19 |
18 | 10.10% | 50 | 0.20 |
19 | 8.60% | 50 | 0.23 |
20 | 8.30% | 50 | 0.24 |
21 | 7.40% | 50 | 0.27 |
22 | 3.80% | 50 | 0.53 |
23 | 2.50% | 50 | 0.80 |
24 | 2.10% | 50 | 0.95 |
25 | 2.10% | 50 | 0.95 |
26 | 1.70% | 50 | 1.00 |
27 | 1.60% | 50 | 1.00 |
28 | 1.30% | 50 | 1.00 |
29 | 1.40% | 50 | 1.00 |
30 | 1.40% | 50 | 1.00 |
步骤208:确定参与风险评价的各项目标物理参数的加权系数。针对岩溶涌水突泥问题,根据大量工程实践经验,采用加权分析法,确定各项目标物理参数对应的加权系数WP,本次风险评价中采用的加权系数如下:
WVs=0.15
WV=0.20
Wρs=0.30
Wη=0.20
Wf=0.15
步骤209:采用加权分析法计算综合权重评分值。
将加权系数分别乘以对应各项物理参数的归一化比值,计算获得综合权重评分值Qi。
统计计算可得,表11示出的综合权重评分值计算表。
步骤210:提取地质风险潜在区域。
绘制超前预报风险评价范围内的综合权重评分值Qi直方图。图3是本发明提供的综合权重评分值直方图,如图3所示,横坐标为距离单位为米,纵坐标为综合权重评分值。
与正相关物理参数经验模型对应的综合权重评分值直方图进行比对分析,提取综合权重评分值偏低的区域,即综合权重评分值Qi<1的区域,于是得到地质风险潜在区域,即第7m至第25m。
地质风险潜在区域数列:{LQl}={Q7,Q8,……,Q25}。
步骤211:定量化风险评价与风险等级划分标准。
参照步骤201建立的风险评价与风险等级划分标准,引入综合权重评分平均值概念,对风险程度进行定量划分;并在地质风险潜在区域中筛选能够满足更高风险等级的地质风险区域,地质风险区域既能满足的更高要求,也能满足地质风险区域规模的定量要求。请参见表12示出的定量化风险评价与风险等级划分标准表。
表11综合权重评分值计算表
步骤212:确定地质风险区域及其风险等级。
参照已建立的定量化风险评价与风险等级划分标准,综合考虑“综合权重评分平均值”和“地质风险区域”两个重要参数,从而得到了本次风险评价的结论:地质风险区域为第12m至第21m,风险类型为岩溶涌水,其风险等级为A级。
表12定量化风险评价与风险等级划分标准表
综上所述,本发明具有以下有益效果:
本发明提出了一种基于多元物理参数经验模型的隧洞超前地质预报风险评价方法,区别于仅能定性推断隧洞地质风险的现状,本发明通过对参与风险分析的物理参数首先进行正相关和归一化预处理,进一步进行加权分析计算和风险区域提取等步骤,创新实现了隧洞地质风险的定量评价与等级划分。
本发明提出了一种基于多元物理参数经验模型的隧洞超前地质预报风险评价方法,区别于现有需要采用较多类型、较多数量参数的评价方法,本发明仅需要使用常规隧洞超前地质预报技术获取的物理参数,且省去了繁琐的室内试验与相关资料搜集过程,明显提高了工作效率。
本发明提出了一种基于多元物理参数经验模型的隧洞超前地质预报风险评价方法,区别于仅能对某单一类型地质风险进行评价的方法,本发明针对各类不良地质体的物性特征与结构形态,针对性选择风险评价物理参数,基本能够适用于隧洞施工过程中面临的各类地质灾害类型,适用面更广。
图4是本发明提供的隧洞超前地质预报风险评价装置的结构示意图,如图4所示,所述装置包括:第一模块401,第二模块402,第三模块403,第四模块404。
第一模块401,用于获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;
第二模块402,用于根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;
第三模块403,用于根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;
第四模块404,用于根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
本发明提供的隧洞超前地质预报风险评价装置,通过获取隧洞的围岩的对应的目标物理参数的分布信号,并且基于加权分析计算预设范围内的不同位置区域的综合权重评分值,根据综合权重评分值的分布情况,实现了对隧洞地质风险状况的评价,该方法操作简单并且能够实现对隧洞地质风险的快速定量评价。
需要说明的是,本发明实施例提供的隧洞超前地质预报风险评价装置,在具体运行时,可以执行上述任一实施例所述的隧洞超前地质预报风险评价方法,对此本实施例不作赘述。
图5是本发明提供的电子设备的结构示意图,如图5所示,该电子设备可以包括:处理器(processor)510、通信接口(Communications Interface)520、存储器(memory)530和通信总线540,其中,处理器510,通信接口520,存储器530通过通信总线540完成相互间的通信。处理器510可以调用存储器530中的逻辑指令,以执行隧洞超前地质预报风险评价方法,该方法包括:获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
此外,上述的存储器530中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面,本发明还提供一种计算机程序产品,所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机执行时,计算机能够执行上述各实施例所提供的隧洞超前地质预报风险评价方法,该方法包括:获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
又一方面,本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各实施例提供的隧洞超前地质预报风险评价方法,该方法包括:获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种隧洞超前地质预报风险评价方法,其特征在于,包括:
获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;
根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;
根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;
根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
2.根据权利要求1所述的隧洞超前地质预报风险评价方法,其特征在于,根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号,包括:
在所述任一目标物理参数与围岩质量呈正相关的情况下,根据所述分布信号与正相关参考平均值的比值,确定所述归一化分布信号;
在所述任一目标物理参数与围岩质量呈负相关的情况下,对所述分布信号进行倒数处理,获取中间分布信号;根据所述中间分布信号与正相关参考平均值的比值,确定所述归一化分布信号。
3.根据权利要求2所述的隧洞超前地质预报风险评价方法,其特征在于,在根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号之前,还包括:
在所述任一目标物理参数与围岩质量呈正相关的情况下,将所述任一目标物理参数的参考平均值作为所述正相关参考平均值;
在所述任一目标物理参数与围岩质量呈负相关的情况下,将所述任一目标物理参数的参考平均值的倒数作为所述正相关参考平均值。
4.根据权利要求1所述的隧洞超前地质预报风险评价方法,其特征在于,在目标物理参数为弹性物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:
利用地震波法,获取在掌子面前方的预设范围内围岩的弹性物理参数的分布信号;
所述弹性物理参数至少包括以下物理参数中的一种:纵波波速、横波波速、泊松比、以及密度。
5.根据权利要求1所述的隧洞超前地质预报风险评价方法,其特征在于,在目标物理参数为电性物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:
在所述电性物理参数为电阻率的情况下,利用瞬变电磁法获取所述电阻率的分布信号;
在所述电性物理参数为电阻率和/或极化率的情况下,利用激发极化法获取所述电性物理参数的分布信号。
6.根据权利要求1所述的隧洞超前地质预报风险评价方法,其特征在于,在目标物理参数为电磁响应物理参数的情况下,所述获取在掌子面前方的预设范围内围岩的任一目标物理参数的分布信号,包括:
利用地质雷达法,获取在掌子面前方的预设范围内围岩的电磁响应物理参数的分布信号;
所述电磁响应物理参数至少包括以下物理参数中的一种:磁波振幅、电磁波频率以及电磁波同相轴形态。
7.根据权利要求1所述的隧洞超前地质预报风险评价方法,其特征在于,根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况,包括:
根据所述综合权重评分值的分布状况和预设评分阈值,确定所述预设范围内的地质风险潜在区域;
根据所述地质风险潜在区域内的综合权重评分值的分布状况以及预先设置的定量化风险评价与风险等级划分标准,确定所述地质风险潜在区域内的地质风险区域,以及所述地质风险区域的定量风险评价与风险等级;
所述地质风险区域的综合权重评分平均值小于预设评分阈值,所述地质风险区域的覆盖范围大于预设范围阈值;所述综合权重评分平均值为所述地质风险区域内的不同位置区域的综合权重评分值的平均值。
8.一种隧洞超前地质预报风险评价装置,其特征在于,包括:
第一模块,用于获取在掌子面前方预设范围内的围岩的任一目标物理参数的分布信号;所述分布信号由位于不同位置区域的围岩的任一目标物理参数的实测值组成;
第二模块,用于根据所述任一目标物理参数与围岩质量的相关关系,对所述分布信号进行归一化处理,获取归一化分布信号;
第三模块,用于根据每个目标物理参数的归一化分布信号与预设加权系数,确定在所述围岩的在不同位置区域的综合权重评分值;
第四模块,用于根据所述综合权重评分值的分布状况,评价所述预设范围内的隧洞地质风险状况。
9.一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至7任一项所述隧洞超前地质预报风险评价方法的步骤。
10.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述隧洞超前地质预报风险评价方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211212264.4A CN115640995A (zh) | 2022-09-30 | 2022-09-30 | 一种隧洞超前地质预报风险评价方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211212264.4A CN115640995A (zh) | 2022-09-30 | 2022-09-30 | 一种隧洞超前地质预报风险评价方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115640995A true CN115640995A (zh) | 2023-01-24 |
Family
ID=84941666
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211212264.4A Pending CN115640995A (zh) | 2022-09-30 | 2022-09-30 | 一种隧洞超前地质预报风险评价方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115640995A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117272213A (zh) * | 2023-11-21 | 2023-12-22 | 中南大学 | 地下污染物的地物化综合参数扫面方法、装置、设备及介质 |
-
2022
- 2022-09-30 CN CN202211212264.4A patent/CN115640995A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117272213A (zh) * | 2023-11-21 | 2023-12-22 | 中南大学 | 地下污染物的地物化综合参数扫面方法、装置、设备及介质 |
CN117272213B (zh) * | 2023-11-21 | 2024-02-02 | 中南大学 | 地下污染物的地物化综合参数扫面方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2016141598A1 (zh) | 共姿态道集地表一致性振幅补偿方法及装置 | |
CN111126471A (zh) | 微地震事件检测方法及系统 | |
CN106066493A (zh) | 贝叶斯岩相判别方法及装置 | |
CN102254095A (zh) | 基于多维分形克里格方法的成矿异常提取方法 | |
US10761229B2 (en) | Microseismic sensitivity analysis and scenario modelling | |
CN113222347B (zh) | 露天爆破安全风险灰色系统评估方法 | |
CN115061219B (zh) | 基于石油天然气探测的裂缝型储层预测识别方法及系统 | |
CN115640995A (zh) | 一种隧洞超前地质预报风险评价方法及装置 | |
CN113534240A (zh) | 微地震事件检测与定位方法及系统 | |
CN105740204A (zh) | 不规则地形下低频段大地电导率快速反演方法 | |
CN107507168B (zh) | 一种用于岩体节理裂隙模型检验的迹线图相似度判别方法 | |
CN110988985B (zh) | 基于波形特征的地震信号检测方法 | |
CN107632010A (zh) | 一种结合激光诱导击穿光谱对钢铁样品的定量方法 | |
Wang et al. | Fragmentation calculation method for blast muck piles in open-pit copper mines based on three-dimensional laser point cloud data | |
Liang et al. | Hybrid support vector machine optimization model for inversion of tunnel transient electromagnetic method | |
CN104155703A (zh) | 评价三维观测系统的方法和设备 | |
CN102709207A (zh) | 质量评估设备、方法和使计算机执行质量评估方法的程序 | |
CN114236624B (zh) | 基于电磁法估算压裂改造空间体积的方法和系统 | |
CN115220098A (zh) | 碳酸盐岩断控缝洞体自动识别方法及装置 | |
CN112395667A (zh) | 一种盾构隧道土体的随机场参数生成方法和装置 | |
CN108875163B (zh) | 一种评价三维裂缝网络连通性的方法和系统 | |
CN112800664A (zh) | 一种基于探地雷达A-scan数据估算树根直径的方法 | |
CN105528657A (zh) | 基于北斗和向量机的建筑物震害预测方法 | |
CN115079248B (zh) | 基于裂纹密度及声波阈值反演爆破损伤范围的方法及设备 | |
Liu et al. | Analyzing spatial variability of geologic profiles for four sites in Hong Kong |
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 |