CN103913776B - 盾构施工中球状孤石的探测方法 - Google Patents
盾构施工中球状孤石的探测方法 Download PDFInfo
- Publication number
- CN103913776B CN103913776B CN201410164382.1A CN201410164382A CN103913776B CN 103913776 B CN103913776 B CN 103913776B CN 201410164382 A CN201410164382 A CN 201410164382A CN 103913776 B CN103913776 B CN 103913776B
- Authority
- CN
- China
- Prior art keywords
- boulder
- signal
- spherical
- shield
- spherical boulder
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种盾构施工中球状孤石的探测方法,要解决的技术问题是提高盾构区间球状孤石的探测能力。本发明的方法包括以下步骤:用井中地质雷达采集到沿钻孔深度方向的地质雷达回波信号,初步圈定球状孤石可能存在的空间方位,在两钻孔之间采集电位差和电流,得到两孔之间地层的真电阻率,将真电阻率绘制成等值线图,圈定得出两孔之间地层中球状孤石的大小和分布情况。本发明与现有技术相比,利用井中地质雷达和跨孔超高密度电法,联合探测盾构隧洞范围内球状孤石,能最大限度地探查隧洞范围内球状孤石的分布,同时还能减少钻孔数量,探测成本较小,同时提高了复杂地层盾构区间发现球状孤石的几率,保证了盾构机的顺利推进。
Description
技术领域
本发明涉及一种隧道施工的地质勘察方法,特别是一种探测盾构区间球状孤石分布的方法。
背景技术
中、微风化花岗岩球状孤石探测是地铁盾构隧道施工建设亟待解决的关键问题。现有技术用于花岗岩球状孤石探测的地质调查方法主要有钻芯取样和工程物探两种。其中,钻芯取样方法最为直接,往往作为地铁设计阶段最主要的地质勘查手段。但钻芯取样只是“一孔之见”,仅能反映钻孔及其周围有限范围内的球状孤石分布,通过钻芯取样揭露球状孤石的信息十分有限。
工程物探利用球状孤石与围岩介质之间的物性差异,如电阻率、磁化率、介电常数、弹性波速度等,来选择合适的物探方法对球状孤石进行探测。用地表观测系统的二维电阻率层析成像和地震波反射能探测到盾构区间10~30米深度内的球状孤石目标,但得到的数据质量受采集环境影响较大,从电阻率剖面与地震映像或叠加剖面上很难圈定球状孤石的大小。地质雷达虽探测精度高,但其穿透能力有限。在富水地层中,100MHz的天线仅能穿透5~10米,远未达到盾构区间10~30米的要求。地微动方法直接利用天然源来探测地层结构的S波速度,该方法能直观显示围岩的纵、横向变化,但其要求天然源为随机源,在交通繁忙地段其适用能力受限,且其探测细粒径球状孤石的能力还需要进一步研究。
现有技术用于山岭隧道地质超前探测的TSP、HSP和TRT地震散射成像方法,震源由炸药或破岩振动产生,地震波接收需要在隧洞壁上钻炸药孔和检波器埋设孔,对于盾构隧道管片衬砌一般不可取。在盾构机刀盘上安装地质雷达天线,同济大学做过类似实验且取得了一些进展,但其真正实用还尚待时日,而且盾构机刀盘结构改造和地质雷达天线保护的结构依旧具有挑战。在盾构机上安装声波软土探测系统或聚焦电流激化极化装置BEAM是国外盾构机常用的做法,这两类方法能实时了解盾构前方地层情况,包括土或岩的强度、均质性以及其它物理特性,特别适合地下工程施工,但由于受国外技术出口限制、采取只租不卖的形式或价格特别昂贵的原因,因此在国内盾构施工现场使用并不常见。
城市的工程物探因受地表条件、金属管线和居民/工业用电等的干扰,勘探精度大大降低。井中物探不但可以避开这些干扰因素,且人工或自然地球物理场更接近勘探目标体,测量数据能有效反映异常体信息,提高了原始数据的信噪比,反演更为准确,勘探精度高。广州地铁建设中针对球状孤石探测问题,选用过多达10余种工程物探方法,先后在3号线和6号线二期工程多次开展球状孤石工程物探方法试验和专题研究,结果表明:瞬变电磁、地面地质雷达、地震映像等地面物探方法均达不到理想效果,但钻孔物探优势明显,不同方法的效果也不尽相同。跨孔超高密度电法和跨孔地震CT虽然能取得较好效果,但其探测范围只是二维剖面,三维探测成本高且技术还不成熟。
发明内容
本发明的目的是提供一种盾构施工中球状孤石的探测方法,要解决的技术问题是提高盾构区间球状孤石的探测能力。
本发明采用以下技术方案:一种盾构施工中球状孤石的探测方法,包括以下步骤:
一、在确定探测盾构隧道区域的上方,沿盾构隧道轴线,在既有地质补勘钻孔里,用井中地质雷达对孔壁周围地质情况进行探测,采集到沿钻孔深度方向的地质雷达回波信号,对地质雷达回波信号做数据处理:零时调整、漂移处理与水平同相轴消除、时域滤波和增益处理,得到整个钻孔深度、沿孔壁周向距离的以回波信号双程走时为横坐标,以距离地表深度为纵坐标的井中地质雷达反射剖面图;
二、在地质雷达反射剖面图中找出代表球状孤石的弧状反射信号线,拾取弧状反射信号线中的弧状的顶点的双程走时t,根据地质补勘中测量的土壤孔隙率和含水饱和度,通过复折射系数(CRIM)公式:
非衰减或弱衰减介质中速度v与相对介电常数ε的关系:
计算雷达电磁波在该地层中的传播速度v;
式(6)中,φ和Sw分别为土壤孔隙率和含水饱和度,εs、εw和εa分别为土颗粒、水和空气的相对介电常数,式(7)中,c为电磁波在真空中传播速度;
通过r=vt/2计算球状孤石的几何中心与孔壁之间的距离r,初步圈定球状孤石可能存在的空间方位r,z,其中z表示球状孤石的深度;
三、在井中地质雷达反射剖面上存在弧状反射信号的两钻孔之间,用分别设置在两孔中的供电电极和测量电极,采集供电电极和测量电极之间电位差ΔU和电流I,再通过最小二乘反演最优化迭代算法,得到两孔之间地层的真电阻率;用于最小二乘反演最优化迭代算法最小二乘反演的目标函数S(m)为:
式(8)中,为二范数意义下的加权残差,m为模型参数,模型参数每一个元素mj和对应电阻率ρj的关系为mj=logρj,j=1,2,...,M,其中M为模型m中元素的个数;m0为先验模型,根据测井或室内试验得电阻率ρ0后再根据对数关系m0=logρ0求出;d为实测视电阻率数据对数矢量,d=log(KΔU/I);K为与电极排列有关的装置系数,满足1/K=I/4π(1/rs-1/rs′),其中rs和rs′分别为供电电极和测量电极所在的位置与测量点所在的空间位置的方向矢量;g(m)为正演模拟视电阻率对数矢量;Wd为数据权重矩阵,Wd=diag(1/εi),其中εi为数据误差;Wm为模型权重矩阵;λ为正则化因子;
采用高斯牛顿迭代法求解上述最小二乘问题:
式(9)中,Δmi为第i次迭代模型的修正量,mi为更新i次后的电阻率对数模型;
以垂直于地表的钻孔深度方向为纵坐标,以两孔连线距离为横坐标,将真电阻率绘制成等值线图,再根据高阻区域内球状孤石电阻率阈值,从等值线图上圈定得出两孔之间地层中球状孤石的大小和分布情况,球状孤石电阻率阈值范围为600~100000Ω·m。
本发明的步骤一零时调整为将每道地质雷达回波信号具有最大振幅波形的起跳时间作为零时刻点;
所述漂移处理与水平同相轴消除为去除地质雷达回波信号中的直流漂移量和水平干扰:
式(1)中,Xs′(i,t)为去除回波信号中直流漂移量后第i道地质雷达的回波信号在时刻t时的数值,X(i,t)为去除回波信号中直流漂移量前第i道地质雷达的回波信号在时刻t时的数值,i为当前地质雷达剖面道号,X(i,k)为去除回波信号中直流漂移量前第i道地质雷达的回波信号在t=kΔt时刻的数值,Δt为时间采样间隔,k为时间采样点,N为总的时间采样点数;式(2)中,X′(i,t)表示去除直流漂移和水平同相轴干扰后的波形数据,Xs′(n,t)为第n道地质雷达回波信号在时刻t时的数值,M为当前回波信号总的信号道数;
所述时域滤波为将X′(i,t)与带通滤波器h(t)进行卷积计算:
X′′(i,t)=X′(i,t)*h(t) (3)
式(3)中,X′′(i,t)为时域滤波后的波形数据,h(t)为时域带通滤波器;
时域理想带通滤波器h(t)的频率特征响应为:
(4)
式(4)中,FT为傅立叶变换算符,f1和f2分别为带通滤波器的低通频率和高通频率;
所述增益处理将每道波形数据X′′(i,t)信号分别与增益函数g(t)相乘:
X′′′(i,t)=X′′(i,t)g(t) (5)
式(5)中,g(t)=(1+at)ebt,a为:a=a′/Lp,b为:b=b′v/8.69,Lp为地质雷达回波信号脉冲宽度,v为雷达波在介质中的传播速度0.1m/nsec,a′=1,b′=5,t为时间。
本发明的步骤一井中地质雷达沿孔壁的探测半径不小于3m,能分辨的球状孤石粒径不小于0.5m,井中地质雷达的竿状天线的发射频率不高于300MHz,时间采样率不低于512,采样窗口长不少于200ns。
本发明的步骤一用于井中地质雷达探测钻孔的孔径不小于90mm,孔底位于隧洞底板以下3~5m。
本发明的步骤一在富水砂层或泥浆难以护住孔壁的松软地层环境中,将内径不小于60mm聚氯乙烯(PVC)管插入孔中构成护壁。
本发明的步骤二在地质雷达反射剖面图中无弧状反射信号,进行下一个地质钻孔探测,直至所有补勘钻孔探测完成。
本发明的步骤三在一个孔地质雷达反射剖面上发现弧状反射信号后,再根据该孔左右两边的孔的地质雷达反射剖面来确定在哪两个孔中之间进行跨孔超高密度电法探测,或直接在与该孔相邻的两孔中进行跨孔超高密度电法探测。
本发明的步骤三对于地下水位较深,孔中最深和最浅两处有效电极水位之间的距离D不满足D≥2L时,L为两孔之间的孔间距,需往钻孔中注水至满足该条件的水位。
本发明的步骤三在富水砂层或泥浆难以护壁的松软地层环境中,采用透水无纺布包裹的开有3~8mm的小孔、孔间距不超过20cm的聚氯乙烯(PVC)管构成护壁。
本发明的步骤三不能圈定球状孤石位置的情况,在两孔连线的中间位置的两侧增加地质钻孔,重复步骤一、步骤二、步骤三,若球状孤石位于原有钻孔与新增钻孔的剖面上,探测完成;
若球状孤石不位于原有钻孔与新增钻孔的剖面上,根据雷达反射三点定位方法确定球状孤石在钻孔构成的三角区域内的空间位置,再用地质钻孔验证。
本发明与现有技术相比,利用井中地质雷达和跨孔超高密度电法,联合探测盾构隧洞范围内球状孤石,能最大限度地探查隧洞范围内球状孤石的分布,同时还能减少钻孔数量,探测成本较小,同时提高了复杂地层盾构区间发现球状孤石的几率,保证了盾构机的顺利推进。
附图说明
图1为本发明的井中地质雷达和跨孔超高密度电法联测示意图。
图2为本发明的探测方法流程图。
图3为本发明实施例的井中地质雷达反射剖面图。
图4为本发明实施例用跨孔超高密度电法得到的真电阻率等值线图。
图5为本发明实施例增加钻孔与雷达反射三点定位方法图。
图6-1为本发明实施例的地质雷达回波单道信号图。
图6-2为本发明实施例零时调整后的地质雷达回波单道信号图。
图6-3为本发明实施例去除直流偏移量和水平干扰后的地质雷达回波单道信号图。
图6-4为本发明实施例时域滤波后的地质雷达回波单道信号图。
图6-5为本发明实施例的增益处理后的地质雷达回波单道信号图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。如图2所示,本发明的盾构施工中球状孤石的探测方法,采用井中地质雷达和跨孔超高密度电法共同测量,包括以下步骤:
一、在分析沿盾构隧道的地质钻孔(孔)资料球状孤石赋存地质条件的基础上,根据勘察设计阶段地质钻孔资料得到的球状孤石区段作为探测的区域,确定需要探测的盾构隧道区域。球状孤石通常赋存在自稳能力不好的残积土中,为花岗岩成岩后期构造作用下风化不均或石英富集部分不容易风化所致。
在确定探测盾构隧道区域的上方,沿盾构隧道轴线,在既有地质补勘钻孔里,用井中地质雷达对孔壁周围地质情况进行探测,采集到沿钻孔深度方向的地质雷达回波信号(回波信号)。对地质雷达回波信号做数据处理,包括零时调整、漂移处理与水平同相轴消除、时域滤波和增益处理,得到整个钻孔深度、沿孔壁周向距离的以回波信号双程走时为横坐标,以距离地表深度为纵坐标的井中地质雷达反射剖面图(反射剖面图,反射剖面)。对地质雷达回波信号做数据处理的具体流程为:
1、对采集到的地质雷达回波信号做零时调整。针对回波信号,将每道回波信号具有最大振幅波形(直达波)的起跳时间t0作为零时刻点。
2、以零时调整后的地质雷达回波信号作为输入信号,进行漂移处理和水平同相轴消除,去除回波信号中的直流漂移量和水平干扰。直流漂移量为回波信号上出现全正或全负、或正负半周不对称的情况所对应的数值。水平干扰为影响特征信号识别的水平同相轴,如直达波或地表反射波。
漂移处理和水平同相轴消除的表达式分别为(《探地雷达原理与应用》第157页第2段、8.1-1和8.1-2式,曾昭发,刘四新,冯晅等编著,电子工业出版社,2010):
式(1)中,Xs′(i,t)为去除回波信号中直流漂移量后第i道地质雷达的回波信号在时刻t时的数值,X(i,t)为去除回波信号中直流漂移量前第i道地质雷达的回波信号在时刻t时的数值,即步骤1中的输出信号,i为当前地质雷达剖面道号,X(i,k)为去除回波信号中直流漂移量前第i道地质雷达的回波信号在t=kΔt时刻的数值,Δt为时间采样间隔,k为时间采样点,N为总的时间采样点数。式(2)中,X′(i,t)表示去除直流漂移和水平同相轴干扰后的波形数据,Xs′(n,t)为第n道地质雷达回波信号在时刻t时的数值,M为当前回波信号总的信号道数。
3、以去除直流漂移和水平同相轴干扰后的波形数据X′(i,t)作为输入进行时域滤波,时域滤波为将X′(i,t)与带通滤波器h(t)进行卷积计算,其表达式为:
X′′(i,t)=X′(i,t)*h(t) (3)
式(3)中,X′′(i,t)为时域滤波后的波形数据,h(t)为时域带通滤波器,“*”为卷积算符。
时域理想带通滤波器h(t)的频率特征响应为:
(4)
式(4)中,FT为傅立叶变换算符,f1和f2分别为带通滤波器的低通频率和高通频率,,高、低通频率定义了信号的通带范围为只允许大于低通频率f1和小于高通频率f2的信号通过。所有理想滤波器的频率响应函数在通带频率处都是间断的,需要有效的镶边函数来避免信号由理想滤波器滤波后产生的吉布斯效应。时域带通滤波器及镶边函数按现有技术方法设计(《信号数字处理的数学原理》第二版第237-239页,程乾生著,石油工业出版社,1993)。
4、将时域滤波后的波形数据X′′(i,t)作为输入进行增益处理,即将每道波形数据X′′(i,t)信号分别与增益函数g(t)相乘,其表达式为:
X′′′(i,t)=X′′(i,t)g(t) (5)
式(5)中,增益函数g(t)为:g(t)=(1+at)ebt,包含线性部分和指数部分,其中,参数a为:a=a′/Lp,参数b为:b=b′v/8.69,Lp为地质雷达回波信号脉冲宽度,可直接通过雷达发射信号或直达波波形宽带直接量得,v为雷达波在介质中的传播速度,通常令v=0.1m/nsec,a′和b′为线性增益和指数增益的控制参数,可根据雷达回波信号的频带范围进行调节,一般令a′=1和b′=5(Karl Josef Sandmeier,Reflex2D-Quick,Karlsruhe,2012.卡尔.约瑟夫.沙德麦,Reflex2D快速参考手册,卡布斯鲁厄,2012,第13页,http://www.sandmeier-geo.de/Download/reflex2dquick.pdf),t为时间。
经增益处理后,便可得到在整个钻孔深度,沿时间或经过r=vt/2时间距离转换后沿孔壁周向距离的井中地质雷达反射剖面图,以回波信号双程走时为横坐标,以距离地表深度为纵坐标。
在确定探测的盾构隧道区域范围内对球状孤石进行探测,对井中地质雷达的要求为:沿孔壁的探测半径不小于3m,能分辨的球状孤石粒径不小于0.5m,井中地质雷达的竿状天线的发射频率不高于300MHz,时间采样率不低于512,采样窗口长不少于200ns。竿状天线设置在孔内,以天线接收、发射天线连线的中点作为测量深度的零点,从地表到孔底的深度测量由地质雷达测距轮完成。
用于井中地质雷达探测钻孔的终孔孔径不小于90mm,孔底位于隧洞底板以下3~5m。隧洞底板指圆形或椭圆形隧道的弧形底部砌体的上表面。在富水砂层或泥浆难以护住孔壁的松软地层环境中,将聚氯乙烯PVC管插入孔中构成护壁,以减少泥浆进入钻孔内,防止孔壁坍塌导致测量时竿状天线或连接竿状天线的电缆被卡。PVC管的内径根据所采用的井中地质雷达天线外径来选择,一般不小于60mm。
二、识别井中地质雷达反射剖面上的弧状反射信号,可以结合地质钻孔资料来判断球状孤石的存在。对反射剖面进行分析,如图3所示,在图中找出代表球状孤石的弧状反射信号线。在非均匀介质中弧状反射信号线为非线段,类似抛物线形状,在复杂介质中弧状反射信号线为无规则曲线。拾取弧状反射信号线中的弧状的顶点(顶点类似抛物线或曲线的顶点)的双程走时t(横坐标)。根据地质补勘中测量的土壤孔隙率和含水饱和度,通过复折射系数CRIM公式(《探地雷达原理与应用》第50页3.1-13式,曾昭发,刘四新,冯晅等编著,电子工业出版社,2010):
非衰减或弱衰减介质中速度v与相对介电常数ε的关系:
计算雷达电磁波在该地层中的传播速度v。
式(6)中,φ和Sw分别为土壤孔隙率和含水饱和度,由地质补勘测量得到,εs、εw和εa分别为土颗粒、水和空气的相对介电常数,可参考《探地雷达原理与应用》第54页表3.5选取经验值(《探地雷达原理与应用》,曾昭发,刘四新,冯晅等编著,电子工业出版社,2010)。式(7)中,c为电磁波在真空中传播速度。
根据雷达电磁波在该地层中的传播速度v,就可通过r=vt/2来计算球状孤石的几何中心与孔壁之间的距离r,初步圈定球状孤石可能存在的空间方位(r,z),其中z表示球状孤石的深度,为相对于地表钻孔口的距离,可以在雷达反射剖面坐标轴上直接读取。
若该孔周围无球状孤石响应(无弧状反射信号),则进行下一个地质钻孔探测,直至所有补勘钻孔探测完成。
由于球状孤石与周围地层存在电性差异,在地质雷达反射剖面上会造成同相轴的不连续,球状孤石多为不规则形状,在地质雷达反射剖面上表现为类似抛物线的弧状反射信号。
三、在井中地质雷达反射剖面上存在弧状反射信号的两钻孔之间,如图1所示的A、B两孔,用跨孔超高密度电法对两孔之间的地层电阻率情况进行探测。用分别设置在A、B两孔中的供电电极和测量电极,采集供电电极和测量电极之间电位差ΔU和电流I,再通过最小二乘反演最优化迭代算法,得到两孔之间地层的真电阻率。也可以在A孔反射剖面上发现弧状反射信号后,再根据A孔左右两边的孔的反射剖面来确定在哪两个孔中之间进行跨孔超高密度电法探测,或直接在与A孔相邻的两孔中进行跨孔超高密度电法探测。
根据文献1、Gunther Thomas,Inversion Methods and Resolution Analysisfor the2D/3D Reconstruction of Resistivity Structures from DCMeasurements.Freiberg University of Mining and Technology,2004,汤玛斯.甘特,利用直流电阻率测量方法重建二维和三维电阻率结构的反演方法及分辨率分析,弗赖堡矿业技术大学,2004第17页;文献2、Rücker Carsten,Advanced Electrical ResistivityModeling and Inversion Using Unstructured Discretization,University ofLeipzig,Dissertation,2011,卡斯顿.拉克,无结构离散的先进电阻率模拟和反演方法研究,博士论文,莱比锡大学,2011第63页;文献3、Kalscheuer Thomas,Improvement andAssessment of Two-Dimensional Resistivity Models Derived from Radiomagnetotelluric and Direct-Current Resistivity Data,Uppsala Uiversity,2008,汤玛斯.卡尔朔伊尔,大地电磁和直流电阻率数据反演二维电阻率模型的改进和评估方法,乌普萨拉大学,2008,第21-22页,得到用于最小二乘反演最优化迭代算法的目标函数S(m)为:
式(8)中,为二范数意义下的加权残差,m为模型参数,模型参数每一个元素mj和对应电阻率ρj的关系为mj=logρj,j=1,2,...,M,其中M为模型m中元素的个数(文献2,第62页);m0为先验模型,根据测井或室内试验得电阻率ρ0后再根据对数关系m0=logρ0求出;d为实测视电阻率数据对数矢量,d=log(KΔU/I);K为与电极排列有关的装置系数,满足1/K=I/4π(1/rs-1/rs′),其中rs和rs′分别为供电电极和测量电极所在的位置与测量点所在的空间位置的方向矢量;g(m)为正演模拟视电阻率对数矢量,可通过求解文献2中2.1式位场方程求得(文献2,第16页),求解过程见文献2第16-56页;Wd为数据权重矩阵,存在关系Wd=diag(1/εi),其中εi为数据误差,可通过仪器测得或者通过估算求得,估算方法见文献2第63页;Wm为模型权重矩阵,其设置方法见文献2第63-64页;λ为正则化因子,通过L曲线法求得(文献2,第66-67页)。
采用高斯牛顿迭代法求解上述最小二乘问题,其相应的迭代格式为:
式(9)中,Ji为雅可比矩阵,为供电电极和测量电极之间的视电阻率对数矢量g(mi)关于电阻率对数模型mi的偏导数Δmi为第i次迭代模型的修正量,mi为更新i次后的电阻率对数模型。
(1)根据文献1第17页公式2.13-2.17;文献2第64页公式4.5;文献3第22页公式4.1-4.6,式(9)的推导如下:将目标函数S(mi+1)在第i次迭代mi处进行泰勒展开,舍去高阶项得:
(1)为了使S(mi+1)逼近二次型(9)最小,那么存在:
(2)求目标函数S(mi)关于模型参数mi的梯度和二次导数分别为:
(3)舍去(12)式中二次导数2的高阶项就可以从式(10)推导出式(9):
将(13)式中的换成雅可比矩阵Ji,便可推导出(9)式。
以垂直于地表的钻孔深度方向为纵坐标,以两孔连线距离为横坐标,将真电阻率绘制成等值线图,如图4所示,得到两孔之间地层的高阻区域,高阻区域为较大等值线值围成的封闭区域。再根据高阻区域内球状孤石电阻率阈值,从等值线图上圈定得出两孔之间地层中球状孤石的大小和分布情况。球状孤石电阻率阈值为花岗岩球状孤石的电阻率,参考范围为600~100000Ω·m。
通过对已知花岗岩球状孤石进行跨孔超高密度电阻率法测量后知,深圳9号线某工区内花岗岩球状孤石阈值为2000Ω·m。
对于地下水位较深(地下水位较低),孔中最深和最浅两处有效电极水位之间的距离D不满足D≥2L时,其中L为两孔之间的孔间距,需往钻孔中注水至满足该条件的水位,在保证测量有效的情况下,通常直接注水到孔口即地表,以保证电极与其周围的介质形成导电通路,且保证孔中水位在测量时不降至从地表往下数的第一个有效电极以下。有效电极为电极与地层接触良好能形成导电通路的电极。有效电极的深度为从地表往下第一个有效电极与最后一个电极或多芯电缆末端电极之间的深度差。本实施例中,在一条电缆上等距离布设32个电极,最后一个电极或末端电极就是多芯电缆末端电极。
用于跨孔超高密度电法的孔的规格与井中地质雷达探测的孔相同,在盾构隧道上方,沿盾构隧道轴线布设的两孔轴线平行,孔间距不大于孔深的1/2。在富水砂层或泥浆难以护壁的松软地层环境中,采用透水无纺布包裹的开有3~8mm的小孔、孔间距不超过20cm的PVC管构成护壁,以减少泥浆进入孔内和防止孔壁坍塌,导致的测量时电缆被卡。用于保护井中地质雷达天线和跨孔超高密度电法电缆的护壁可为同一PVC管,只是用于跨孔超高密度电法的PVC管需要在管壁上钻孔,且用透水无纺布包裹PVC管。
四、参考地质钻孔资料的球状孤石的赋存地质条件,根据井中地质雷达的反射剖面和跨孔超高密度电法得到的等值线图,结合球状孤石的赋存地质条件,根据高阻区域内球状孤石电阻率阈值,则可准确圈定球状孤石在两孔之间的空间方位和大小。球状孤石通常赋存在自稳能力不好的残积土中,为花岗岩成岩后期构造作用下风化不均或石英富集部分不容易风化所致。
对步骤一中探测到球状孤石响应,步骤二找出反射剖面中代表球状孤石的弧状反射信号,而步骤三用跨孔超高密度电法不能圈定球状孤石位置的情况,如图5所示,在盾构隧道上方地面沿盾构隧道轴线,在两孔A、B连线的中间位置的两侧增加地质钻孔C、D,重复上述步骤一、步骤二、步骤三,若球状孤石位于原有钻孔A与新增钻孔C或D,或原有钻孔B与新增钻孔C或D的连线的剖面上,则探测完成。
若球状孤石不位于原有钻孔A与新增钻孔C或D,或原有钻孔B与新增钻孔C或D的连线的剖面上,根据雷达反射三点定位方法确定球状孤石在钻孔A、C、D或B、C、D构成的三角区域内的空间位置,再用地质钻孔验证。地质钻孔验证的具体方法为:直接在定位点上方地表进行地质钻探,根据取芯记录球状孤石上、下底界面距离地表的深度。
如图5所示,雷达反射三点定位方法确定球状孤石空间位置的方法具体为:选取钻孔位置B为坐标原点,分别拾取与其相邻的三钻孔A、C、D的井中地质雷达反射剖面图上针对同一疑似球状孤石的弧状反射信号线中弧状的顶点的电磁波发射和回波的双程走时ti和深度zi,i=A,C,D,计算孔壁与孤石之间的距离ri,通过A、C和D三钻孔空间坐标xi和yi,连同孤石弧状反射深度zi,建立三个三元二次方程:
式(14)中,x,y,z为球状孤石坐标,用现有技术非线性最小二乘求解,便可确定球状孤石存在的空间方位。非线性最小二乘求解为:先将三个三元二次方程写成目标函数形式:
然后求目标函数E关于空间位置p=(x,y,z)的一阶导数和二阶导数,分别得到梯度G和Hessian矩阵H,构成线性方程组:
HΔp=-G (16)
给定一初始位置p0=(x0,y0,z0),由牛顿法迭代求解上式得:
pi+1=pi+αΔpi,i=1,2,... (17)
则可得到球状孤石空间位置坐标(x,y,z),其中α为步长因子,可由线搜索求得。
实施例1,为探测盾构隧道区间范围内球状孤石赋存情况,在深圳地铁9号线某地质补勘区域内用本发明的方法,采用井中地质雷达和跨孔超高密度电法共同测量,对可能存在球状孤石的区域进行勘察。具体步骤为:
一、钻孔孔位布设如图1所示,用井中地质雷达对孔壁周围地质情况进行探测,将采集到沿钻孔深度方向的地质雷达回波信号做数据处理,选取第i=300道数据为例,如图6-1所示。
1、零时调整,将回波信号上该道信号具有最大振幅波形的起跳时间t0=4ns作为零时刻点,将零时刻点以上的数据进行切除,切除处理后的结果如图6-2所示;
2、以零时调整后的回波信号X(i,t)作为输入信号,按照式(1-1)和(1-2)去除回波信号中的直流漂移量和水平干扰:
式(1-1)中,N=512为时间采样点数;式(1-2)中,M=1028为同一测线总的信号道数。去除直流偏移量和水平干扰后的结果如图6-3所示。
3、进行时域滤波,时域滤波为将X′(i,t)与带通滤波器h(t)按照式(1-3)进行卷积计算,对于150MHz井中天线,带通滤波器的高、低通频率分别为f1=75MHz和f2=300MHz。
X′′(i,t)=X′(i,t)*h(t) (1-3)
其中带通滤波器的镶边函数采用余弦函数,镶边宽度为△f=10MHz。选择的h(t)其频率响应满足:
(1-4)
对H(f)做反傅立叶变换就可得到h(t)。滤波后的结果如图6-4所示。
4、将波形数据X′′(i,t)作为输入按照式(1-5)进行增益处理,g(t)=(1+at)ebt,其中,a=0.04ns-1,b=0.0575m/ns。
X′′′(i,t)=X′′(i,t)g(t) (1-5)
增益处理后的结果如图6-5所示。
将地质雷达反射剖面每一道数据按照上述步骤1-4处理便可得到在整个钻孔深度,沿时间的井中地质雷达反射剖面,如图3所示,以回波双程走时为横坐标,以距离地表深度为纵坐标。
实施例1对井中地质雷达的要求为:沿孔壁的探测半径不小于3m,能分辨的球状孤石粒径不小于0.5m,井中地质雷达的竿状天线的发射频率150MHz,采样时间0.25ns,采样窗口长200ns。用于井中地质雷达探测钻孔的终孔孔径90mm,孔底位于隧洞底板以下5m。
二、对反射剖面进行分析,如图3所示,图中箭头所示处代表球状孤石的弧状反射信号线,拾取弧状反射信号线中的弧状的顶点的双程走时t=97ns。根据地质补勘中测量的土壤孔隙率φ=45.9%和含水饱和度Sw=86.2%,以及在150MHz天线频率下实验测量的土颗粒相对介电常数参考值εs≈3、水的相对介电常数εw≈81和空气的相对介电常数εa=1,通过复折射系数CRIM公式:
和非衰减或弱衰减介质中速度v和相对介电常数ε的关系(1-7):
计算出雷达电磁波在该地层中的传播速度v为0.066m/ns。
根据该地层中的传播速度v,通过r=vt/2来计算球状孤石的几何中心与孔壁之间的距离r=3.2m,初步圈定球状孤石可能存在的空间方位为距离井壁3.2m且深度为18m。
三、在井中地质雷达反射剖面上存在弧状反射信号的两钻孔之间,如图1所示的A、B两孔,用跨孔超高密度电法对两孔之间的地层情况进行探测。用分别设置在A、B两孔的供电电极和测量电极,供电电极和测量电极末端深度为24米,采集两电极之间电位差ΔU和电流I,再通过关系d=log(KΔU/I)和目标函数(1-8):
采用高斯牛顿迭代法求解上述最小二乘问题可得第i次迭代模型的修正量Δmi,更新后的模型参数为:
mi+1=mi+Δmi,i=1,2,... (1-9)
重复上面的求解过程直到满足条件|S(mi+1)-S(mi)|/S(mi)<10-6,迭代终止,做指数变换后便可以得到两孔之间地层的真电阻率。
如图4所示,以垂直于地表的钻孔深度方向为纵坐标,以两孔连线距离为横坐标,将真电阻率绘制成等值线图,得到两孔之间地层的高阻区域,高阻区域为较大等值线值围成的封闭区域,如图4中的星号处。再根据高阻区域内球状孤石电阻率阈值,从等值线图上分析得出两孔之间地层中球状孤石的大小和分布情况。深圳9号线某工区内花岗岩球状孤石阈值为2000Ω·m。
根据高阻区域内球状孤石电阻率阈值,准确圈定球状孤石在两孔之间的空间方位为距离地表18m深度、中心距离右井壁3.2m,其长度在深度方向2m,在水平方向粒径为0.5m。
实施例1探测球状孤石效果明显,后经钻孔验证探测准确。
本发明利用既有补勘地质钻孔,不受地表场地限制,探测电极或天线接近不良地质目标体,受环境干扰小,能保证探测准确性。借助单孔中地质雷达无极化和方向的特点,探测孔壁周围360度范围内球状孤石的分布情况,并能根据地质雷达反射剖面图判断出球状孤石离孔壁距离,能准确界定球状孤石是否在隧洞区间内,具有较高的分辨率。采用跨孔超高密度电法进行跨孔探测,能弥补单孔中地质雷达探测无法定位方向的缺点,并能准确成像两孔之间真电阻率变化,根据真电阻率变化能准确判断球状孤石的位置和大小。本发明采用井中地质雷达和跨孔超高密度电法相结合,结合了两种方法的优势,能最大限度地探查隧道范围内球状孤石的分布,同时还能减少钻孔数量,探测成本较小,具有较好的经济效益。
Claims (10)
1.一种盾构施工中球状孤石的探测方法,包括以下步骤:
一、在确定探测盾构隧道区域的上方,沿盾构隧道轴线,在既有地质补勘钻孔里,用井中地质雷达对孔壁周围地质情况进行探测,采集到沿钻孔深度方向的地质雷达回波信号,对地质雷达回波信号做数据处理:零时调整、漂移处理与水平同相轴消除、时域滤波和增益处理,得到整个钻孔深度、沿孔壁周向距离的以回波信号双程走时为横坐标,以距离地表深度为纵坐标的井中地质雷达反射剖面图;
二、在地质雷达反射剖面图中找出代表球状孤石的弧状反射信号线,拾取弧状反射信号线中的弧状的顶点的双程走时t,根据地质补勘中测量的土壤孔隙率和含水饱和度,通过复折射系数(CRIM)公式:
非衰减或弱衰减介质中速度v与相对介电常数ε的关系:
计算雷达电磁波在地层中的传播速度v;
式(1)中,φ和Sw分别为土壤孔隙率和含水饱和度,εs、εw和εa分别为土颗粒、水和空气的相对介电常数,式(2)中,c为电磁波在真空中传播速度;
通过r=vt/2计算球状孤石的几何中心与孔壁之间的距离r,初步圈定球状孤石可能存在的空间方位r,z,其中z表示球状孤石的深度;
三、在井中地质雷达反射剖面上存在弧状反射信号的两钻孔之间,用分别设置在两孔中的供电电极和测量电极,采集供电电极和测量电极之间电位差ΔU和电流I,再通过最小二乘反演最优化迭代算法,得到两孔之间地层的真电阻率;用于最小二乘反演最优化迭代算法最小二乘反演的目标函数S(m)为:
式(3)中,为二范数意义下的加权残差,m为模型参数,模型参数每一个元素mj和对应电阻率ρj的关系为mj=logρj,j=1,2,…,M,其中M为模型m中元素的个数;m0为先验模型,根据测井或室内试验得电阻率ρ0后再根据对数关系m0=logρ0求出;d为实测视电阻率数据对数矢量,d=log(KΔU/I);K为与电极排列有关的装置系数,满足1/K=I/4π(1/rs-1/r′s),其中rs和r′s分别为供电电极和测量电极所在的位置与测量点所在的空间位置的方向矢量;g(m)为正演模拟视电阻率对数矢量;Wd为数据权重矩阵,Wd=diag(1/εi),其中εi为数据误差;Wm为模型权重矩阵;λ为正则化因子;
采用高斯牛顿迭代法求解上述最小二乘问题:
mi+1=mi+Δmi,i=1,2,... (4)
式(4)中,Δmi为第i次迭代模型的修正量,mi为更新i次后的电阻率对数模型;
以垂直于地表的钻孔深度方向为纵坐标,以两孔连线距离为横坐标,将真电阻率绘制成等值线图,再根据高阻区域内球状孤石电阻率阈值,从等值线图上圈定得出两孔之间地层中球状孤石的大小和分布情况,球状孤石电阻率阈值范围为600~100000Ω·m。
2.根据权利要求1所述的盾构施工中球状孤石的探测方法,其特征在于:步骤一中的零时调整为将每道地质雷达回波信号具有最大振幅波形的起跳时间作为零时刻点;
所述漂移处理与水平同相轴消除为去除地质雷达回波信号中的直流漂移量和水平干扰:
式(5)中,X′s(i,t)为去除回波信号中直流漂移量后第i道地质雷达的回波信号在时刻t时的数值,X(i,t)为去除回波信号中直流漂移量前第i道地质雷达的回波信号在时刻t时的数值,i为当前地质雷达剖面道号,X(i,k)为去除回波信号中直流漂移量前第i道地质雷达的回波信号在t=kΔt时刻的数值,Δt为时间采样间隔,k为时间采样点,N为总的时间采样点数;式(6)中,X′(i,t)表示去除直流漂移和水平同相轴干扰后的波形数据,X′s(n,t)为第n道地质雷达回波信号在时刻t时的数值,M为当前回波信号总的信号道数;
所述时域滤波为将X′(i,t)与时域带通滤波器h(t)进行卷积计算:
X″(i,t)=X′(i,t)*h(t) (7)
式(7)中,X″(i,t)为时域滤波后的波形数据,h(t)为时域带通滤波器;
时域带通滤波器h(t)的频率特征响应为:
式(8)中,FT为傅立叶变换算符,f1和f2分别为带通滤波器的低通频率和高通频率;
所述增益处理将每道波形数据X″(i,t)信号分别与增益函数g(t)相乘:
X″′(i,t)=X″(i,t)g(t) (9)
式(9)中,g(t)=(1+at)ebt,a为:a=a′/Lp,b为:b=b′v/8.69,Lp为地质雷达回波信号脉冲宽度,v为雷达波在介质中的传播速度0.1m/nsec,a′=1,b′=5,t为时间。
3.根据权利要求2所述的盾构施工中球状孤石的探测方法,其特征在于:井中地质雷达沿孔壁的探测半径不小于3m,能分辨的球状孤石粒径不小于0.5m,井中地质雷达的竿状天线的发射频率不高于300MHz,时间采样率不低于512,采样窗口长不少于200ns。
4.根据权利要求3所述的盾构施工中球状孤石的探测方法,其特征在于:用于井中地质雷达探测钻孔的孔径不小于90mm,孔底位于隧洞底板以下3~5m。
5.根据权利要求4所述的盾构施工中球状孤石的探测方法,其特征在于:所述步骤一在富水砂层或泥浆难以护住孔壁的松软地层环境中,将内径不小于60mm聚氯乙烯(PVC)管插入孔中构成护壁。
6.根据权利要求1所述的盾构施工中球状孤石的探测方法,其特征在于:在地质雷达反射剖面图中无弧状反射信号时,进行下一个地质钻孔探测,直至所有补勘钻孔探测完成。
7.根据权利要求1所述的盾构施工中球状孤石的探测方法,其特征在于:当在一个孔的地质雷达反射剖面上发现弧状反射信号后,再根据该孔左右两边的孔的地质雷达反射剖面来确定在哪两个孔中之间进行跨孔超高密度电法探测,或直接在与该孔相邻的两孔中进行跨孔超高密度电法探测。
8.根据权利要求1所述的盾构施工中球状孤石的探测方法,其特征在于:当地下水位较深、且孔中最深和最浅两处有效电极水位之间的距离D不满足D≥2L时,需往钻孔中注水至满足该条件的水,其中L为两孔之间的孔间距。
9.根据权利要求1所述的盾构施工中球状孤石的探测方法,其特征在于:在富水砂层或泥浆难以护壁的松软地层环境中,采用透水无纺布包裹的开有3~8mm的小孔、孔间距不超过20cm的聚氯乙烯(PVC)管构成护壁。
10.根据权利要求1所述的盾构施工中球状孤石的探测方法,其特征在于:当出现不能圈定球状孤石位置的情况时,在两孔连线中间位置的两侧增加地质钻孔,重复步骤一、步骤二、步骤三;若球状孤石位于原有钻孔与新增钻孔的剖面上,则探测完成;
若球状孤石不位于原有钻孔与新增钻孔的剖面上,根据雷达反射三点定位方法确定球状孤石在钻孔构成的三角区域内的空间位置,再用地质钻孔验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410164382.1A CN103913776B (zh) | 2014-04-23 | 2014-04-23 | 盾构施工中球状孤石的探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410164382.1A CN103913776B (zh) | 2014-04-23 | 2014-04-23 | 盾构施工中球状孤石的探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103913776A CN103913776A (zh) | 2014-07-09 |
CN103913776B true CN103913776B (zh) | 2016-08-10 |
Family
ID=51039600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410164382.1A Active CN103913776B (zh) | 2014-04-23 | 2014-04-23 | 盾构施工中球状孤石的探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103913776B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104407394B (zh) * | 2014-11-26 | 2015-07-22 | 山东大学 | 基于电阻率跨孔ct的地铁盾构区间孤石精细化探测方法 |
CN104391334B (zh) * | 2014-11-26 | 2017-01-04 | 山东大学 | 用于地下水运移过程监测的电阻率时间推移反演成像方法 |
CN104678443A (zh) * | 2015-01-22 | 2015-06-03 | 深圳市市政设计研究院有限公司 | 一种孤石探测仪及其控制方法 |
US10330618B2 (en) * | 2015-04-30 | 2019-06-25 | Schlumberger Technology Corporation | Method to estimate water saturation in electromagnetic measurements |
CN105604066B (zh) * | 2015-12-25 | 2017-11-03 | 天津市勘察院 | 电阻率剖面法在建筑基坑围护结构渗漏水检测中的应用 |
CN105804763B (zh) * | 2016-03-11 | 2017-03-29 | 山东大学 | 盾构掘进机搭载的超前三维电阻率跨孔ct孤石探测系统及探测方法 |
CN105807336A (zh) * | 2016-03-11 | 2016-07-27 | 山东大学 | 适用于多地球物理探测方法的综合孤石探测模型试验装置及试验方法 |
CN106324683B (zh) * | 2016-09-08 | 2019-07-12 | 山东大学 | 用于地铁盾构隧道前方孤石探测的声波装置及方法 |
CN106932772A (zh) * | 2017-03-15 | 2017-07-07 | 华北计算技术研究所(中国电子科技集团公司第十五研究所) | 一种面向数字地球的受地形影响的雷达探测范围显示方法 |
CN107436435A (zh) * | 2017-08-07 | 2017-12-05 | 广州地铁设计研究院有限公司 | 一种单孔探测孤石的探地雷达装置及其探测方法 |
CN108415100A (zh) * | 2018-03-01 | 2018-08-17 | 广州轨道交通建设监理有限公司 | 一种用于花岗岩球状风化体的高精度水平综合勘探方法 |
CN109188383B (zh) * | 2018-08-16 | 2023-03-14 | 广西有色勘察设计研究院 | 富水圆砾层盾构管片壁后注浆地质雷达探测方案 |
CN109252868A (zh) * | 2018-10-23 | 2019-01-22 | 中铁十二局集团有限公司 | 软硬不均地层盾构施工孤石及基岩凸起处理方法 |
CN110161576A (zh) * | 2019-06-13 | 2019-08-23 | 中国地质大学(北京) | 一种矿区复垦地表层土壤非均质程度的确定方法及系统 |
CN110375913B (zh) * | 2019-07-23 | 2021-01-12 | 中铁第四勘察设计院集团有限公司 | 一种盾构法隧道非预埋式结构健康监测方法 |
CN110469334B (zh) * | 2019-09-06 | 2020-10-02 | 广西大学 | 一种用于地铁盾构孤石和突出岩石爆破的爆破系统 |
CN111929347A (zh) * | 2020-07-31 | 2020-11-13 | 中铁九局集团第四工程有限公司 | 一种盾构机施工孤石探测装置及探测方法 |
CN115880482B (zh) * | 2023-02-17 | 2023-05-09 | 中海油田服务股份有限公司 | 一种测井图像遇卡识别与倾角校正方法、装置及计算设备 |
CN116360001A (zh) * | 2023-04-04 | 2023-06-30 | 中国地质科学院地球物理地球化学勘查研究所 | 消除冻土区低频非屏蔽探地雷达多次波信号干扰的方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5052166B2 (ja) * | 2007-03-06 | 2012-10-17 | 株式会社小松製作所 | 埋設物検出装置、管理装置、および建設機械監視システムの埋設物探索方法 |
CN102508308B (zh) * | 2011-11-03 | 2013-06-05 | 山东大学 | 一种隧道瞬变电磁多点阵列式探测方法 |
CN102518442B (zh) * | 2011-12-14 | 2014-04-23 | 中铁十二局集团第二工程有限公司 | 一种盾构掘进溶沟地段孤石群的处理方法 |
-
2014
- 2014-04-23 CN CN201410164382.1A patent/CN103913776B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103913776A (zh) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103913776B (zh) | 盾构施工中球状孤石的探测方法 | |
Martínez-Moreno et al. | Integrated geophysical methods for studying the karst system of Gruta de las Maravillas (Aracena, Southwest Spain) | |
Gómez-Ortiz et al. | Assessing the risk of subsidence of a sinkhole collapse using ground penetrating radar and electrical resistivity tomography | |
Martínez-Pagán et al. | The electrical resistivity tomography method in the detection of shallow mining cavities. A case study on the Victoria Cave, Cartagena (SE Spain) | |
dos Reis Jr et al. | Characterization of collapsed paleocave systems using GPR attributes | |
CN105116453B (zh) | 一种冻土带天然气水合物的瞬变电磁勘探方法及装置 | |
Cueto et al. | Karst‐induced sinkhole detection using an integrated geophysical survey: a case study along the Riyadh Metro Line 3 (Saudi Arabia) | |
CN107462924B (zh) | 一种不依赖于测井资料的绝对波阻抗反演方法 | |
CN105204078B (zh) | 一种计算三维大地电磁测深主阻抗的方法 | |
Verhegge et al. | Exploring integrated geophysics and geotechnics as a paleolandscape reconstruction tool: archaeological prospection of (prehistoric) sites buried deeply below the Scheldt Polders (NW Belgium) | |
CN107065014A (zh) | 一种基于微动探测结合地质钻探的地铁区间孤石群探测方法 | |
Hazreek et al. | Soil identification using field electrical resistivity method | |
Ramazi et al. | Contribution of geophysical inversion theory and geostatistical simulation to determine geoelectrical anomalies | |
Tao et al. | Evaluating the joint use of GPR and ERT on mapping shallow subsurface features of karst critical zone in southwest China | |
Zhao et al. | Application of 3D GPR attribute technology in archaeological investigations | |
Abdul-Nafiu et al. | Effects of electrode spacing and inversion techniques on the efficacy of 2D resistivity imaging to delineate subsurface features | |
Ganiyu et al. | Combined electrical resistivity imaging and ground penetrating radar study for detection of buried utilities in Federal University of Agriculture, Abeokuta, Nigeria | |
Durdağ et al. | Combined application of electrical resistivity tomography and multi-channel analysis of surface waves methods in the tunnel detection: A case study from Kocaeli University Campus Site, Turkey | |
CN116084927A (zh) | 地铁线路采空区的三维微动探测方法、装置及介质 | |
Danilov et al. | Study of deep structure of the kimberlite pipe named after M. Lomonosov of the Arkhangelsk diamondiferous province obtained by joint using of passive seismic and radiometric methods | |
Al-Zubedi et al. | Use of 2D azimuthal resistivity imaging in delineation of the fracture characteristics in Dammam aquifer within and out of Abu-Jir fault zone, central Iraq | |
Gobashy et al. | High resolution ground magnetic survey (HRGM) for determining the optimum location of subsurface dam in Wadi Nu’man, Makkah Al Mukarammah, KSA | |
CN104360390A (zh) | 基于cmpcc二维面波的吸收散射综合分析法 | |
Monahan | Investigating fault structure using electrical resistivity tomography | |
Abd El-Aal et al. | Geophysical and borehole investigations of subsurface hazardous karsts in Ar-Riyadh city, Kingdom of Saudi Arabia |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |