一种基于大地电磁测深电阻率判定隧围岩等级的方法技术

本发明专利技术公开了一种基于大地电磁测深电阻率判定隧道围岩等级的方法，包括以下步骤：（1）根据已知区域的隧道围岩等级与洞身电阻率数值样本，建立隧道围岩等级与电阻率区间的对应关系模型；（2）获取待分级区域的电阻率数值；（3）根据步骤2获得的待分级区域的电阻率以及步骤1建立的隧道围岩等级与电阻率区间的对应关系模型，得出待分级区域的隧道围岩等级。采用上述技术方案，具有较强的适用性与定量性，克服野外调查和人工判定的误差问题，克服钻孔波速确定的无法连续计算的缺陷，克服地震波速求取的困难性问题，可以根据大地测深电阻率判定围岩等级，为隧道设计提供围岩分级可靠的依据。

【技术实现步骤摘要】
一种基于大地电磁测深电阻率判定隧围岩等级的方法
本专利技术涉及一种基于大地电磁测深电阻率判定隧道围岩等级的方法，属于地质勘探施工
，尤其涉及铁路隧道围岩等级划分方法领域。
技术介绍
近年来随着高速铁路的发展，复杂地区的铁路建设越来越多，因此隧道的设计工作越多越来，隧道围岩分级是隧道设计最重要的依据。不同围岩等级的隧道，工程造价差异极大，围岩等级不准确可能导致项目资金预算缺口使得项目无法顺利进行。不同的围岩等级需要采取不同的施工工艺，围岩等级判断不准确将带来极大的施工风险。围岩分类是选择施工方法的依据、是进行科学管理及正确评价经济效益、确定结构上的荷载(松散荷载)、确定衬砌结构的类型及尺寸、制定劳动定额、材料消耗标准等的基础。隧道设计阶段，隧道围岩等级的确定依据根据野外调查、洞身处岩石的波速、物探异常、埋深等因素。野外调查是通过收集该地区的地质资料、现场踏勘、野外露头地质测绘，分析岩石的地质年代、地层岩性等确定岩体的隧道围岩等级。该方法仅仅根据露头岩石的特征分析其围岩等级，由于无法得到深部岩石的特性误差较大。确定洞身处岩石的波速主要有两种方法，一种是钻孔取芯测波速，该方法直观、准确，但仅仅代表该孔处的岩石波速，无法得到整个地区的岩石波速，且对环境有一定的破坏作用，一种是地震反射波波确定岩石波速，但该方法在复杂地区无法实施。物探异常、埋深等因素，一般来讲，隧道口都给5级，浅埋段给4级(埋深小于40～50m)，埋藏较深的一般给3级，很少给2级，如果有断裂等构造，物探异常、岩性接触等不利因素的段落，一般会提升一级。该方法人为因素较大，误差较大。因此，现有技术中的隧道围岩等级划分存在着不同的缺陷，急需一种简便快捷的实用方法来解决这种难题。
技术实现思路
因此，本专利技术的目的在于提供一种基于大地电磁测深电阻率判定隧道围岩分级的方法，根据大地电磁测深电阻率确定隧道围岩等级，为隧道设计提供可靠的依据。为了实现上述目的，本专利技术的一种基于大地电磁测深电阻率判定隧道围岩等级的方法，包括以下步骤：(1)根据已知区域的隧道围岩等级与洞身电阻率数值样本，建立隧道围岩等级与电阻率区间的对应关系模型；(2)获取待分级区域的电阻率数值；(3)根据步骤2获得的待分级区域的电阻率以及步骤1建立的隧道围岩等级与电阻率区间的对应关系模型，得出待分级区域的隧道围岩等级。所述电阻率数值通过以下步骤获取：获取标准成果为surfer格式的大地电磁数据；对GRID数据文件以及洞身坐标BLN数据文件进行提取；采用MapGen软件读入提取出的GRID文件和洞身坐标BLN数据文件，计算出洞身线与GRID数据交点处的横坐标、纵坐标和电阻率数值并输出。当通过MapGen软件提取隧道洞身电阻率时，如果隧道洞身线与计算的电阻率数值没有交点，采用电阻率数值加权平均得到。至少对电阻率数值进行连续两次正态分布概率密度函数计算，并将每次计算结果小于99.7％的数值去除。所述隧道围岩等级与电阻率区间的对应关系为不同岩性的隧道围岩等级与电阻率区间的对应关系。所述大地电磁数据的采集间距为5～10m。统计的电阻率数值个数大于20000个。根据大地电磁采集的影响因素，去除由于采集的原因导致的电阻率的异常值。绘制电阻率断面图，按照洞身地层岩性及电阻率进行分段解释说明电阻率分布情况、电阻率范围、隧道围岩等级，实现隧道的围岩等级评价。根据电阻率数值样本的来源地区确定基于电阻率区间判定隧道围岩等级的适用地区。采用上述技术方案，本专利技术的基于大地电磁测深电阻率判定隧道围岩等级的方法，与现有技术相比，具有较强的适用性与定量性，克服野外调查和人工判定的误差问题，克服钻孔波速确定的无法连续计算的缺陷，克服地震波速求取的困难性问题，可以根据大地测深电阻率判定围岩等级，为隧道设计提供围岩分级可靠的依据。附图说明图1为提取GRID数据文件的脚本代码截图。图2为对不同围岩等级中不同点的电阻率的原始数据示意图。图3为第一次正太分布后的数据示意图。图4为第二次正太分布后的数据示意图。图5为第三次正太分布后的数据示意图。图6为绘制电阻率剖面的示意图。具体实施方式以下通过附图和具体实施方式对本专利技术作进一步的详细说明。本专利技术提供一种基于大地电磁测深电阻率判定隧道围岩等级的方法，根据某地区已知的电阻率数值以及隧道围岩等级样本建立数据模型，通过获取该地区的待分级区域的大地电磁数据并提取电阻率数值，并通过已建立的数据模型得出隧道围岩等级。在本实施例中，以东北地区已经竣工的隧道为例，统计隧道围岩等级数值，收集该地区对应的大地电磁资料，利用MapGen软件提取洞身处的电阻率，此次统计共有8种岩性21055个样本点，统计过程中将明显的异常点去除。具体而言，通过大地电磁探测获取了统一处理平台下标准参数处理的成果后，为了和洞身围岩进行比对，需要提取洞身处电阻率的数值。获取的标准成果为surfer格式，为此，在surfer平台下进行洞身电阻率的提取，即在surfer平台下对成果图进行网格化，以洞身高程在网格中画线，洞身线和测点的交点定义为网格化的坐标，根据坐标找出电阻率属性，把电阻率数值读取出来，形成电阻率文件。Surfer图形中包含了用GRID生成的等值线，首先对GRID文件提取，但是Surfer操作接界面中却没有提供GRID数据文件的提取功能，通过图1所示的脚本代码，可以实现Surfer隐藏的GRID文件的提取功能，为数据的进一步处理和分析奠定了基础。其次，对洞身线bln文件提取。洞身线有的是用实际坐标通过导入bln白化边界文件实现的，有的是用多段线在Surfer画线来表示的洞身线。BLN白化文件一般是GRID数据文件的边界、特征线等，此处特指洞身线，为了进行后续的洞身线数据提取，需要把洞身线文件图提取出来。不同情况采用不同的处理方式，对于真实坐标的洞身线处理，可以直接使用Surfer平台进行提取，对于用多段线在Surfer画线来表示洞身线情况，可以通过鼠标进行点选提取。洞身数据的计算与获取采用专门的软件MapGen来完成。该软件读入提取出的GRID文件和洞身坐标bln文件，通过数学算法计算出洞身线与GRID数据交点处的横坐标、纵坐标和电阻率数值，并形成文件输出。如果隧道洞身线与计算的电阻率数值没有交点，采用电阻率数值加权平均得到。利用电阻率进行围岩等级划分，需要收集大量的电阻率和围岩等级数据。在收集的资料中有不同的单位和人实施了不同的方法，受现场人文噪声以及雷暴等二次场场源影响不同，不同地质人员认识的不同确定的围岩等级存在差异，不同岩性电阻率受岩性、水文条件、风化等影响的差异，这些差异对某些电阻率来说可能有比较大的影响，但放入大数据以后，可以把以上差异产生的误差认为是随机而离散的。在大数据的状态下，单一岩性电阻率和围岩分级应服从正态分布。正态分布也称“常态分布”，又名高斯分布，是一个在数学、物理及工程等领域都非常重本文档来自技高网...

【技术保护点】
1.一种基于大地电磁测深电阻率判定隧道围岩等级的方法，其特征在于，包括以下步骤：/n（1）根据已知区域的隧道围岩等级与洞身电阻率数值样本，建立隧道围岩等级与电阻率区间的对应关系模型；/n（2）获取待分级区域的电阻率数值；/n（3）根据步骤2获得的待分级区域的电阻率以及步骤1建立的隧道围岩等级与电阻率区间的对应关系模型，得出待分级区域的隧道围岩等级。/n

【技术特征摘要】
1.一种基于大地电磁测深电阻率判定隧道围岩等级的方法，其特征在于，包括以下步骤：
（1）根据已知区域的隧道围岩等级与洞身电阻率数值样本，建立隧道围岩等级与电阻率区间的对应关系模型；
（2）获取待分级区域的电阻率数值；
（3）根据步骤2获得的待分级区域的电阻率以及步骤1建立的隧道围岩等级与电阻率区间的对应关系模型，得出待分级区域的隧道围岩等级。


2.如权利要求1所述的基于大地电磁测深电阻率判定隧道围岩等级的方法，其特征在于，所述电阻率数值通过以下步骤获取：获取标准成果为surfer格式的大地电磁数据；对GRID数据文件以及洞身坐标BLN数据文件进行提取；采用MapGen软件读入提取出的GRID文件和洞身坐标BLN数据文件，计算出洞身线与GRID数据交点处的横坐标、纵坐标和电阻率数值并输出。


3.如权利要求2所述的基于大地电磁测深电阻率判定隧道围岩等级的方法，其特征在于：当通过MapGen软件提取隧道洞身电阻率时，如果隧道洞身线与计算的电阻率数值没有交点，采用电阻率数值加权平均得到。


4.如权利要求2所述的基于大地电磁测深电阻率判定隧道围岩等级的方法，其特征在于：至少对电阻率数值进行连续两次正态分布概率密度函数计算，并将每次计算结果小于99.7%...

【专利技术属性】
技术研发人员：许广春，秦海旭，祁晓雨，刘英，牛永效，
申请(专利权)人：中国铁路设计集团有限公司，
类型：发明
国别省市：天津;12