微动探测方法在城市地铁盾构施工“孤石”探测中的应用——以福州地铁1号线为例

刘宏岳， 黄佳坤， 孙智勇， 宗全兵

(1. 福建省建筑设计研究院， 福建 福州 350001； 2. 福州市城市地铁有限责任公司， 福建 福州 350001)

微动探测方法在城市地铁盾构施工“孤石”探测中的应用
——以福州地铁1号线为例

刘宏岳1， 黄佳坤1， 孙智勇2， 宗全兵2

(1. 福建省建筑设计研究院， 福建 福州 350001； 2. 福州市城市地铁有限责任公司， 福建 福州 350001)

孤石在我国南方花岗岩地区普遍存在，其分布无明显规律，形状各异，强度可达到100 MPa以上。福州地铁1号线在施工过程中遇到了孤石问题，盾构区间范围内的孤石如果没有提前查明并处理，会给盾构施工造成重大不利影响和安全隐患。采用微动探测方法有效解决了福州地铁1号线的孤石问题，运用微动探测成果结合少量钻孔验证进行盾构施工的安全性分析，为盾构施工提供指导。简要阐述了微动探测方法的基本原理，介绍其在福州地铁1号线多个盾构区间孤石探测的成功案例，说明微动探测方法适用于城市复杂的环境条件，在盾构施工孤石探测方面具有良好的应用效果。

微动探测； 福州地铁； 盾构施工； 孤石

0 引言

花岗岩风化土中的球状风化核，俗称“孤石”，主要由不易风化的石英矿物形成石英角砾质残留核(球状风化核)组成，在我国南方燕山期花岗岩发育地区普遍存在，其形状各异，大小从几十厘米到几米不等，强度可以达到100 MPa以上[1-2]。孤石会给地铁施工带来极大风险，盾构掘进中若遇到未探明的孤石，经常会严重损坏盾构，甚至会造成喷涌、塌方等意外情况；被动地处理孤石可能会造成环境破坏、工期延误、路面交通堵塞和投资控制的不确定性增大等问题，造成的损失往往是巨大的。因此，在盾构施工开始之前探明孤石非常有必要。

孤石探测方法一般分为钻探和地球物理2大类。由于孤石的发育无规律性，埋藏及分布较为随机，很难通过地质钻探查明其分布情况。虽然通过加密钻孔可以提高发现孤石的概率，降低工程风险，但受成本和场地条件等限制，钻探常常难于实施；同时，在城市道路中往往存在密集分布的电力、电信、雨水、污水、燃气和路灯的地下复杂管道等障碍物，钻探本身就具有很高风险[3]。因此，寻求一种适应城市复杂环境的地球物理方法来探测孤石势在必行。

城市地铁盾构施工孤石探测是公认的技术难题，如广州地铁建设中针对孤石的专题研究，采用十余种物探方法，由于地铁建设所处的城市综合环境复杂，电磁波类方法受城市复杂电磁背景干扰严重，常规地震方法及电法类受城市狭小地面空间和复杂地表情况影响较大，各类跨孔CT方法虽有一定效果，但也受到孔间距等影响，诸多物探方法均达不到预期效果[4-10]。

福州地铁1号线的盾构施工同样遇到了孤石问题[11]，对此，采用微动方法解决了孤石探测这一技术难题。这里介绍福州地铁1号线几个盾构区间孤石探测的成功案例，希望可对类似工程提供参考。

在世界上不同的地区和国家，微动有不同的名称。在美洲，称为被动源面波(Passive Surface Wave)；在日本，称为微动(Microtremor)；而在欧洲则是环境随机振动(Ambient Vibration)；在我国，有时又叫作天然源。尽管名称各异，但实质相同。如今，在美洲、欧洲和日本等地区和国家，微动探测方法在无损检测和场地地震评价中的应用日益广泛，但运用于浅地表工程领域的孤石探测基本没有。国内徐佩芬等[12-14]在这方面开展了比较多的研究工作，成功将微动方法应用于探测地热断层、煤炭陷落柱等中深部勘探以及岩土工程尺度的城市地铁勘察；在深圳地铁7号线采用微动方法得到2个剖面共19个勘探点的速度，并对可能存在的花岗岩孤石进行划分，但未见钻探验证资料。总体而言，微动探测方法在岩土工程方面的应用较少，在国内尚处于起步阶段。

1 微动与微动勘探

地球表面存在一种微弱波动，它源于自然界和人类的各种活动。自然界中的风、潮汐、气压变化和火山活动等都会产生振动，而人类活动产生的振动包括车辆移动、工厂机械运行，甚至人的行走等。所有这些振动的能量都以波的形式向远处传播，其中含有各种体波，但能量传播的主要形式是面波，因此微动方法也可以称为被动源面波法。

微动震源通常分为2大类——自然现象和人类活动，以对应震源的不同频段。低频段信号反映深部地层信息，而岩土工程尺度的探测主要使用高频段的信号，如孤石探测要求的深度较浅，其震源主要为测点附近数百米内的交通和机器振动。当震源距离台阵较近时，微动波场包含体波和面波；当震源距离台阵较远时，面波则为主要传播形式。

微动是一种由体波(P波和S波)和面波(Rayleigh波和Love波)组成的复杂振动，并且面波的能量占信号总能量的70%以上，大部分的能量以基阶模式传播。垂直分量检波器仅能探测到Rayleigh面波，水平分量的检波器能同时探测Rayleigh面波和Love面波。尽管微动信号的振幅和形态随时空变化而发生变化，但在一定时空范围内具有统计稳定性，可用时间和空间上的平稳随机过程描述。微动探测就是以平稳随机过程理论为依据，从微动信号中提取面波(Rayleigh波)的频散曲线，通过对频散曲线的反演，获得地下介质的横波速度结构[15]。

2 微动探测工作原理简介

微动探测是一种基于微动台阵探测的地球物理探测方法。微动探测工作原理如图1所示，采用台阵(圆形台阵或者内嵌三角形台阵等)布置方式进行信号采集，运用空间自相关(SPAC)算法[16]或频率波数域(F-K)算法[17]从各测点微动台阵记录中提取Rayleigh波频散曲线，然后直接绘制相速度等值线图，或者计算视S波速度vx，再经插值光滑计算获得二维视S波速度剖面。相速度等值线图或者视S波速度剖面均能客观且直观地反映地层岩性变化，是地质解释的基本依据。本文中案例采用高分辨率频率波数域(HRFK)算法。

图1 速度剖面获取流程图

在实际计算中，将台阵微动信号划分为多个时间窗，用HRFK法计算各个时间窗和不同频率对应的慢度(速度的倒数)，统计所有时间窗的慢度概率分布，舍弃误差大且明显偏离频散曲线形态的点，求出各个频率的慢度平均值，即频散曲线。

福州地铁1号线隧洞直径为6.2m，一般采用半径为2.5～3m的正五边形阵列[18]，如图1观测台阵，每个圆形阵列由放置于正五边形顶点和中心点的6个检波器和数据采集系统组成，按5m测点间距逐点进行，以形成二维剖面观测，观测阵列可基本覆盖整个隧洞范围。现场测试仪器可采用地震仪与2Hz三分量检波器通过电缆连接方式同步采集振动信号，也可以采用一体化数字地震仪采集，后者采用无线连接方式，通过GPS授时功能实现各台地震仪的同步信号采集。

微动方法探测“孤石”以孤石与周围风化土体的速度差异为地球物理前提，在测试成果形成的速度剖面上圈定高速异常区域，并通过钻孔验证直接揭露孤石或者起到排除孤石的作用。

3 福州地铁1号线微动探测案例

福州地铁1号线采用微动方法探测盾构施工不良地质体，共完成探测长度超过8km，1 600多个台阵测点，根据微动探测成果共提出验证钻孔96个，实际完成验证钻孔78个，其中揭露孤石20个，揭露基岩凸起16个，发现岩脉6处，验证钻孔揭露的最小孤石40cm。根据钻孔的验证结果，进行盾构掘进施工的安全性分析，划分出盾构施工的危险区、警示区和安全区，对盾构施工起到指导作用。从探测成果来看，微动探测取得了理想的效果。

在验证的78个钻孔中，41个孔发现不良地质体(包括孤石、基岩凸起、岩脉和注浆混凝土块)，准确率53%，若扣除茶瓶站—达道站区间18个孔，其他区间的不良地质验证准确率高达68%。茶瓶站—达道站区间异常实际验证的结果主要为红黏土地层，该地层速度高，出现地层速度倒转。

福州地铁1号线目前微动探测区域已全部完成洞通，福州地铁公司反馈微动划分的安全区100%准确，微动成果未漏过直径大于1m的孤石，危险区域经施工单位的提前处理和施工预案，未出现盾构卡住的现象。

以下介绍福州地铁1号线新店出入段、屏山站—东街口站区间和达道站—上藤站区间3个盾构区间的不良地质体微动探测案例。

3.1 案例1： 新店出入段

图2为新店出入段入段线RDK1+700～+500视S波速度剖面图，其中平行双实线为隧道洞身范围，结合详勘资料，隧道洞身上方速度明显偏高的土体推断为碎卵石层，厚度为5m到7m多不等。新店出入段为福州地铁1号线最早采用微动技术探测的一个试验段，盾构在前期施工中发现一些异常，在盾构前方200m区域进行探测。洞身范围内RC13—RC21测点表现为较大范围土体速度明显偏高，推断为碎卵石层深入到隧道洞身范围。由于施工工期紧，且场地条件不允许，未进行钻孔验证和处理，盾构掘进过程中在RC13—RC20遇到孤石(大粒径碎卵石)被卡住。图3为该段盾构掘进中所取孤石岩芯样本(照片由施工方提供)，孤石已被盾构搅碎，能够顺利取出，取出的孤石最大直径为53cm。

图2 新店入段线RDK1+700～+500视S波速度剖面图

(a)(b)

图3 孤石岩芯样本

Fig. 3 Photos of boulder core samples

3.2 案例2： 屏山站—东街口站区间

该区域地质资料显示场地基岩岩性主要为花岗岩，局部地段辉绿岩呈脉岩状产出。图4中屏山站—东街口站区间下行线里程XK7+840～XK8+610完成探测长度770 m，台阵测点共155个，图4中红色点线为隧道的洞身范围(根据各频点的面波半波长穿透深度推算)。根据地质资料和验证钻孔揭露岩性，推断图中速度较大的4处地方应为后期侵入的岩脉，但风化程度有所不同。

图4 屏山站—东街口站XK7+840～XK8+610面波相速度等值线图

根据探测成果剖面可知，该剖面共布置6个异常验证钻孔，实际完成5个，验证结果见表1。

表1 钻孔验证结果

其中,X4揭露厚度达4 m的孤石;X7揭露基岩侵入隧道洞身；X34揭露碎块状；X13验证孔由于交通原因，钻孔位置偏离较远。值得注意的是测点X100验证孔洞身范围内未揭露异常，但在同一个微动台阵范围内，靠近X99测点补充验证孔洞身范围揭露到孤石，孤石厚6.5 m。所布钻孔均位于4处岩脉中，其余区段无明显速度异常，直接列为盾构掘进的安全区，该段盾构隧道现已贯通，隧道其余位置未出现妨碍盾构掘进的孤石，安全区的准确度达100%。

3.3 案例3： 达道站—上藤站区间

图5为达道站—上藤站区间下行线XK12+415.81～+609.81面波相速度等值线图，图5中红色点为根据各频点的面波半波长穿透深度推算的隧道顶底深度范围，根据该剖面成果，共布置5处验证钻孔。

其中X2和X8验证钻孔揭露花岗岩孤石，厚度分别为0.4 m和1.9 m，见图6，其余钻孔未见异常。

图5 达道站—上藤站XK12+415～+710面波相速度等值线图

(a)X2(b)X8

图6 X2和X8验证孔岩芯照片

Fig. 6 Testifying borehole core samples of X2 and X8

综合探测成果剖面和钻孔验证结果，对该段进行盾构掘进的安全性分析，针对揭露孤石的X2和X8位置附近应进一步查明孤石的大小和分布位置，XK12+415.81～+455列为盾构掘进的危险区，剩余区段XK12+455～+609.81，其中X1、X14和X22测点位置未揭露孤石，但能起到排除孤石可能性的作用，钻孔验证到的2处孤石位于强风化花岗岩地层中，X2孔孤石埋深17.0～17.8 m，在洞顶附近;X8孔孤石埋深20.6～22.5 m，在洞身范围。孤石的尺寸相对埋深比例较小，前者为1∶21，后者1∶11，突破了常规物探方法径深比1∶10的限制，分辨率也远超常规地震勘探四分之一波长的限制。达道站—上藤站区间在下行线探测到的异常验证到的最小孤石为40 cm。

4 结论与讨论

微动探测有效利用环境噪声，并从中获取面波的频散特性以推断地下速度结构，这使得微动探测技术非常适用于城市的复杂环境，市区繁忙的交通不仅不影响观测，还为浅层微动勘探提供了丰富的高频信号源。微动探测不受场地条件限制，可在交通繁忙、建筑物密集、钻探难于实施的闹市区或地质信息盲区进行有效探测，是一种无损、经济、高效的地球物理探测手段。微动探测与少量钻孔结合，可以得到较精确的地下构造二维剖面，面波相速度等值线图和视S波速度剖面能直观显示岩土层的纵、横向变化，横向速度对异常情况有明显反映，可推断划分出盾构掘进位置正常地层与可疑孤石地层的界限，经少量钻孔校正与验证，可辨识出盾构施工掘进的安全区、警示区和危险区位置，以便提早处理孤石或准备施工方案。

微动探测对孤石异常有明显的反映，孤石埋深估算较准确，一般说来孤石越大，速度异常也越明显，但对孤石的大小、三维形态无法辨识。对于探测验证到的孤石，目前施工单位常规做法仍采用钻探以密间距确定边界范围，并结合孤石爆破一起处理，也有采用钻孔孔间CT穿透确定异常孤石边界[19]。由于孤石的尺寸相对台阵面积一般较小，下一步的研究应考虑微动台阵采用更多的检波器确定孤石异常在台阵中的位置。

微动探测成果划分盾构施工的危险区、警示区和安全区综合利用了相速度和H/V曲线2个参数，H/V曲线可以反映测点位置从浅部到深部的波阻抗变化规律，该部分的内容本文暂未涉及，以后将另文阐述。

微动探测对一定埋深的小尺度孤石异常有较好的反映，分辨率相当高，与常规地震勘探相比提高了一个数量级，笔者认为与微动探测利用众多的震动源有关系，国内目前尚未见到这方面的文章，期待更多的同行进行这方面的研究。

[1] 鲍晓东.深圳地区花岗岩残积土工程特性的研究[J].铁道勘察，2004 (2)： 72-74.(BAO Xiaodong. Engineering geological properties of the remainder soil from weathering granite in Shenzhen[J]. Railway Investigation and Surveying, 2004(2): 72-74.(in Chinese))

[2] 刘宏岳，梁奎生，段建庄.地震反射波CDP叠加技术在台山核电海域花岗岩孤石探测中的应用[J].隧道建设，2011，31(6)： 657-661.(LIU Hongyue, LIANG Kuisheng, DUAN Jianzhuang. Application of marine seismic reflection CDP stacking technology in detection of weathered granite boulder: Case study of water intake tunnel of Taishan Nuclear Power Station[J]. Tunnel Construction, 2011, 31(6): 657-661. (in Chinese))

[3] 古力.盾构机破碎孤石条件及预处理方法[J].隧道建设，2006，26(增刊2): 12-13,22.(GU Li. Qualifications for boulder broke by shield machines and pre-treatment of boulders [J]. Tunnel Construction, 2006, 26(S2): 12-13,22. (in Chinese))

[4] 王典.地铁工程球状风化专题勘察方案研究[J].广州建筑，2011，39(2): 42-46.(WANG Dian. Discussion about special investigation programmes of spheroidal weathering in subway project[J]. Guangzhou Architecture, 2011, 39(2): 42-46. (in Chinese))

[5] 王俊超，师学明，万方方，等.探测孤石高阻体的跨孔电阻率CT水槽物理模拟实验研究[J]. CT理论与应用研究，2012，21(4): 647-657．(WANG Junchao, SHI Xueming, WAN Fangfang, et al. The physical experiment research on laboratory of cross-hole electric resistivity tomography for detecting boulders[J]. CT Theory and Applications, 2012, 21(4): 647-657. (in Chinese))

[6] 党如姣.孤石的物探探测方法[J].隧道建设，2012，32(增刊2): 56-60.(DANG Rujiao. Geophysical prosfecting method for boulder[J]. Tunnel Construction,2012, 32 (S2): 56-60. (in Chinese))

[7] 李洋，党如姣，朱培民，等.花岗岩风化残留体地震法探测应用研究[J].隧道建设，2015，35(5)： 435-438.(LI Yang, DANG Rujiao, ZHU Peimin, et al. Application of seismic method in detection of weathered-granite boulders [J]. Tunnel Construction, 2015, 35 (5): 435-438. (in Chinese))

[8] 刘成禹，林毅鹏，林超群，等.球状孤石在探地雷达探测成果中的表现特征[J].物探与化探，2015，39(4): 860-866.(LIU Chengyu, LIN Yipeng, LIN Chaoqun, et al. The performance characteristic of spherical boulder in the georadar detection[J]. Geophysical and Geochemical Exploration, 2015, 39(4): 860-866. (in Chinese))

[9] 姚金.高密度电法在某地铁工程孤石勘察中的应用[J].勘察科学技术，2015(4): 57-59,64.(YAO Jin. Application of high density resistivity method in boulder detecting in a Metro engineering[J]. Site Investigation Science and Technology, 2015(4): 57-59,64. (in Chinese))

[10] 朱亚军，龙昌东，杨道煌，等.井间电磁波CT技术在地下孤石探测中的应用[J].工程地球物理学报，2015，12(6): 721-725.(ZHU Yajun, LONG Changdong, YANG Daohuang, et al. The application of cross borehole electromagnetic wave CT technology to subsurface isolated stone detection[J]. Chinese Journal of Engineering Geophysics, 2015, 12 (6): 721-725. (in Chinese))

[11] 张庆林，詹晨曦，陈晔翔.福州地铁1号线工程特性及地质风险研究[J].福建地质，2014，33(2): 142-148.(ZHANG Qinglin, ZHAN Chenxi, CHEN Yexiang. Study of engineering properties and geological risk of Fuzhou Metro Line 1 [J]. Geology of Fujian, 2014, 33(2): 142-148. (in Chinese))

[12] 徐佩芬，李传金，凌甦群，等.利用微动勘察方法探测煤矿陷落柱[J].地球物理学报，2009，52(7): 1923-1930.(XU Peifen, LI Chuanjin, LING Suqun, et al. Mapping collapsed columns in coal mines utilizing microtremor survey methods [J]. Chinese Journal of Geophysics, 2009, 52(7): 1923-1930. (in Chinese))

[13]XU Peifen, LING Suqun, LI Chuanjin, et al. Mapping deeply-buried geothermal faults using microtremor array analysis[J]. Geophysical Journal International, 2012, 188(1): 115-122.

[14] 徐佩芬，侍文，凌甦群，等.二维微动剖面探测“孤石”： 以深圳地铁7号线为例[J].地球物理学报，2012，55(6): 2120-2128.(XU Peifeng, SHI Wen, LING Suqun, et al. Mapping spherically weathered “Boulders” using 2D microtremor profiling method: A case study along subway Line 7 in Shenzhen[J]. Chinese Journal of Geophysics, 2012, 55(6): 2120-2128. (in Chinese))

[15] 赵东.被动源面波勘探方法与应用[J].物探与化探，2010，34(6): 759-764.(ZHAO Dong. Passive surface waves: Methods and applications[J]. Geophysical & Geochemical Exploration, 2010, 34(6): 759-764. (in Chinese))

[16] Aki Keiiti. Space and time spectra of stationary stochastic waves, with special reference to microtremors[J]. Bull. Earthq. Res. Inst, 1957, 35: 415-457.

[17]Capon J. High-resolution frequency-wavenumber spectrum analysis[J]. Proceedings of the IEEE, 1969, 57(8): 1408-1418.

[18]Hiroshi Okada. Theory of efficient array observations of microtremors with special reference to the SPAC method [J]. Exploration Geophysics, 2006, 37(1): 73-85.

[19] 郭云峰，胡晓娟，游敬密.地震波跨孔CT探测孤石中直线与曲线追踪法的应用[J].工程地球物理学报，2015，12(3): 361-366.(GUO Yunfeng, HU Xiaojuan, YOU Jingmi. The exploitation of straight and curve tracing method in boulders prospected by seismic cross-hole tomography[J]. Chinese Journal of Engineering Geophysics, 2015, 12 (3): 361-366. (in Chinese))

Application of Microtremor Method to Boulders Detection in Urban Metro Shield Construction: Case Study of Fuzhou Metro Line No.1

LIU Hongyue1, HUANG Jiakun1, SUN Zhiyong2, ZONG Quanbing2

(1. Fujian Provincial Institute of Architectural Design and Research, Fuzhou 350001, Fujian, China;2.FuzhouUrbanMetroCo.,Ltd.,Fuzhou350001,Fujian,China)

The boulders (spherical weathered granite) exist commonly in south China; and the compression strength of boulder can reach 100 MPa. The boulders occurred during shield tunneling of Fuzhou Metro Line No. 1, which would lead to potential safety risk. The microtremor method is adopted to detect the boulders of Fuzhou Metro Line No. 1. The application results are testified by borehole drilling. The working principle of microtremor detection method is introduced. The application results of the microtremor detection method to Fuzhou Metro Line No. 1 show that the above-mentioned method is reliable and effective in complex urban area.

microtremor detection; Fuzhou Metro; shield tunneling; boulder

2016-05-03;

2016-07-18

刘宏岳(1967—)，男，福建永泰人，1990年毕业于同济大学，勘查地球物理专业，本科，教授级高级工程师，主要从事工程物探、地震数据处理研究及微动技术研究工作。E-mail: 475068460@qq.com。

10.3973/j.issn.1672-741X.2016.12.014

U 45

A

1672-741X(2016)12-1500-07