2. 中国科学院地球科学研究院, 北京 100029;
3. 中国科学院大学, 北京 100049;
4. 同济大学土木工程防灾国家重点实验室, 上海 200092;
5. 成都理工大学地质灾害防治与地质环境保护国家重点实验室, 成都 610059
2. Institutes of Earth Science, Chinese Academy of Sciences, Beijing 100029, China;
3. University of Chinese Academy of Sciences, Beijing 100049, China;
4. State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China;
5. State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China
面波是一种沿地表传播的地震弹性波,其中的Rayleigh波成分是面波勘探的主要有效信号.Rayleigh(1885)在自由表面约束下解波动方程时,从理论上证明了Rayleigh波的存在:波速比体波小,随距离衰减比体波慢,但随深度呈指数迅速衰减,其波前是一个圆柱面,质点振动为逆时针旋转的椭圆.Haskell(1953)提出了Rayleigh波在多层介质中的频散曲线的计算方法,成为利用Rayleigh波获取地下介质岩土物性的理论基础.
按Rayleigh波信号来源的不同,Rayleigh波勘探方法分为主动源法和被动源法(Dai et al., 2015; 刘庆华等,2015).主动源方法是利用人工可控震源在地表激发,产生稳态或瞬态的地震波场,从地震记录中选取有效的Rayleigh波信号,进而提取频散曲线.受限于主动源波场的频带宽度,该方法主要应用于浅层面波勘探.较早被提出的方法是基于双道检波器的稳态Rayleigh波互相关法(Jones, 1958)和瞬态面波谱分析法(Stokoe and Nazarian, 1983),后被发展成多道观测(multichannel analysis of surface waves method,MASW)(刘云祯和王振东,1996;Park et al., 1999a).MASW法因其对体波干扰不敏感、线性排列和识别低速层等优势,被更普遍地应用(Miller et al., 1999; Yilmaz and Eser, 2002; Tian et al., 2003; Lin et al., 2004).被动源Rayleigh波信号频带相对较宽,可用于深部构造分析.因其非侵入式和低成本的优势,也常被用于城市浅层地质勘探和地质灾害调查.被动源方法利用非特定震源产生的地震波和大地噪声来提取Rayleigh波频散曲线.例如Aki(1957)提出的空间自相关法(SPAC)和Capon (1970)提出的频率-波数法(F-K);后来,Nakamura (1989)提出了水平-垂向分量谱比法(HVSM),Okada和Suto(2003)提出了扩展空间自相关法(ESPAC),提高了SPAC法的实用性(徐佩芬等,2009;张宝龙等,2016).Louie(2001)结合MASW法和谱分析法提出了一种折射微动法(ReMi),与此类似的,Park和Miller(2008)也在MASW法的基础上提出了路边被动MASW法(roadside passive MASW method),可用于滑坡、公路、机场等特殊环境中的面波勘探(Coccia et al., 2010).由此可见,随着Rayleigh波勘探技术的发展,特别是MASW法,处理主动源和被动源数据的技术隔阂越来越小,适用范围越来越广(夏江海等,2015).
上述方法利用大地噪声或表面震源激发的面波来获取介质信息,而对于浅层地下震源激发的面波及其近场效应的探讨鲜有报道.随着城市地下轨道交通的快速发展,地铁振动对周围建筑和居民带来结构震损、噪音干扰等隐患的同时(罗桂纯等,2011;杨立志等,2011),也为基于埋入源的MASW方法提供了潜在应用场景.地铁轮轨之间的动力学作用激发的振动,以地震波的形式传到地表,携带了地下介质的本征信息(剪切波速、基频、振幅和介质泊松比等),这些地震波具有被用于地质勘探的潜力.不过到目前为止,很少有学者对埋入源(如地铁)引发的面波进行地质、地球物理研究.就浅层勘探而言,考虑到成本问题,一般不会舍易就难去人工设置埋入源,因此如果能够利用既有的地下震源,例如普遍存在的地铁,会是一个新的勘探思路.特别是考虑到这种方法的非侵入式、低成本的优势,在城市环境中的应用意义更加值得关注.
根据前人对地表震源的MASW方法的研究,发现虽然MASW方法观测系统为线性排列,施工方便,数据采集快速高效,其中一个必须考虑的问题就是如何设置最小偏移距.由于面波在震中一定距离外才会产生,应当远离震源一段距离后开始布设检波器.但是偏移距的选取存在很强的经验性.因为MASW方法所采用的线性排列的参数(如最小偏移距、道间距和排列长度等)直接决定了数据质量,进而影响了频散图像的分辨率、探测深度和反演结果的可信度.具体而言,最小偏移距与近场效应直接相关,后者会造成低频段波速难以估计.道间距和排列长度的大小则分别与频散图的高频段和低频段分辨率有关.Stokoe等(1994)建议最优的偏移距应为1.5倍的Rayleigh波长;Park等(2001)根据合成MASW地震记录的频散分析,主张使用尽可能长的排列.Zhang等(2004)通过推导双层介质中Rayleigh波的频谱方程,基于介质的横波速度范围和界面深度的先验信息,提出了最优偏移距的理论公式.Xu等(2006)也基于双层介质提出了最小偏移距的理论公式.Dikmen等(2010)计算了不同偏移距和排列参数的频散图,通过对比频散图的分辨率,给出了以3~4倍道间距或1/3排列长度作为最小偏移距的结论.
以上对MASW方法的排列参数所做的研究为成功实施基于表面震源的MASW野外勘探提供了许多重要的指导资料.但对于埋入源的MASW方法,尚无可参考的公开信息.博采前人研究,为推动被动源的多道面波勘探,本文发展了一种最小偏移距的经验估计方法,并以城市环境中的地铁振动为实例,探讨该方法的有效性.
1 埋入源激发Rayleigh波的传播机制 1.1 Rayleigh波形成的波动地震学依据
Rayleigh波本质上并不是新的波,而是P波和SV波的在自由表面相互耦合的叠加波,波速小于体波(Rayleigh,1885),如图 1a所示.SV波以大于临界角(θc)入射到自由表面时,会发生全反射,如图 1b所示(Ewing et al., 1957;Lay and Wallace, 1995).由于横波速度(vs)比纵波速度(vp)小,设SV波入射角为α,P波反射角为β,根据斯奈尔定律(Snell′s Law),显然有α<β,且sinβ=(vp/vs)sinα>(vp/vs)sinθc=1.由波动地震学理论可知,此时会产生一个沿自由表面传播的、波速
根据几何地震学中的射线理论,本文分析了埋入源激发的Rayleigh波在层状介质中的传播机制,以点震源为例,如图 2所示,介质参数满足vsi>vsi-1,vpi>vpi-1,第n层为均匀无限半空间.
由于P波和SV波在透过界面时都可以产生SV波,因此入射波可能是它们两者之一.根据斯奈尔定律,只有以震源SV波入射和各层SV波透射的路径才可以在自由表面产生不均匀P波,进而产生Rayleigh波,且第i层的入射角θi满足
(1) |
记x0为Rayleigh波产生时的震中距,d表示界面深度.在MASW法中,检波器布设在离震源不足x0的距离以内是无法接收到Rayleigh波信号的.因此,x0可被视作“最小偏移距”.在n层层状介质中,不难发现x0和每层介质的厚度以及震源深度都是正切关系.x0可由下式表达:
(2) |
式中di是第i层介质的厚度,dn是震源到最后一层界面的距离,vsi是第i层介质的剪切波速,要比vp1小.
根据弹性波中vs、vp和泊松比σ的关系,
(3) |
式(2)可以进一步表示成
(4) |
式中σ1是第一层介质的泊松比,κn是剪切波速比vsn/vs1,d为震源深度.
依据实际情况,埋入源一般在25 m以浅.对于浅表地层而言,泊松比σ1的值一般在0.35到0.45之间.通过对第四纪浅表地层的剪切波速统计,分别以5 m和25 m地层的剪切波速作为vs1和vsn,剪切速度比值κn的统计结果(见表 1)大多在1.22到1.72之间.因此式(4)中的最小偏移距x0可以根据埋入源的深度d进行经验估计:
(5) |
基于σ1和κn的平均值,即σ1=0.4,κn=1.47,等式(5)可以简化成
(6) |
为了抑制近场效应,第一个检波器到震中的距离应该超过式(6)所计算出的经验估计值.
式(6)给出了在多层介质中计算最小偏移距的经验公式.这里需要说明的是,由于射线理论是对波场的高阶近似(Sheriff and Geldart, 1995),如果介质中存在断层、不整合面、尖灭点或采空区,那么Rayleigh波在复杂的近地表介质中的传播无法用封闭式的函数表达(解析解).因此利用式(6)的方法计算的结果会不准确.但在实际的MASW方法应用中,上述公式可用作最小偏移距的近似估计.
2 埋入源激发Rayleigh波的有限元数值分析为了验证本文提出的最小偏移距公式,进而发展基于埋入源的MASW方法,需要对Rayleigh波在实际的复杂介质中的传播情况进行探讨.如前文所述,Rayleigh波在复杂的浅层介质中传播的解析解并不存在.因此,可以利用数值模拟法对其进行研究(Gucunski and Woods, 1992),本文采用有限单元法(finite element method,FEM)(王妙月和郭亚曦,1987),选取了上海某地铁隧道及邻区作为研究区,从数值分析的角度探讨基于埋入源的Rayleigh波传播机制,验证本文提出的最小偏移距公式,以期为MASW法的实际应用提供指导.
2.1 最小偏移距的经验估计本文利用ABAQUS(v6.13)有限元模拟软件,建立了二维的地铁隧道及邻区的地质模型.由于地铁隧道呈左右环状对称,因此在进行数值模拟时只对其中的右半边进行建模,如图 3所示.通过调研,获得了研究区的钻孔资料——土层介质的物性和几何参数(见表 2).地铁的铁轨道床位于地下16.15 m处,与左边界距离2.3 m.
模型的几何尺寸设为100 m长,50 m深.为了更合理地模拟研究区的波场信息(Gardien and Stuit, 2003),模型右边界和下边界设为无限边界(Herrmann, 2013),上边界为自由边界,左边界则采用对称边界(Yang et al., 2003; Kouroussis et al., 2011)以满足原地质模型的对称条件.
对于本文所用的复杂多层介质,由于隧道的存在,式(2)和式(6)给出的最小偏移距的计算结果并不严格精确.假如不考虑隧道对土体的影响,即假设介质是层状连续的,震源位于道床上,则式(2)仍可用作最小偏移距的理论估算.考虑到道床与模型左边界的距离为2.3 m,理论的最小偏移距如式(7)所示:
(7) |
根据式(6),由于道床位于地下16.15 m,可以得到最小偏移距的经验值:
(8) |
由于隧道空洞会造成波场的衍射,因此最小偏移距应大于式(7)和(8)的计算结果.保守估计,为了能有效地接收到埋入源产生的Rayleigh波的信号,最小偏移距应大于15 m.
2.2 有限元分析的模型参数根据Zerwer等(2002)的Rayleigh波有限元建模方法,单元尺寸(Le)由最小波速vmin和最高频率fmax决定.由于Rayleigh波是地震波中波速最小的,因此最小波速可取vR,而最高频率fmax可由物理模型测定或参考其他先验信息.
(9) |
地铁引起的地表振动频率通常在0~90 Hz(Lai et al., 2005;楼梦麟等,2009),因此单元尺寸应不大于0.25×96.2/90=0.267 m.本文采用0.25 m×0.25 m的平面应变有限单元(CPE4R),无限边界设为平面应变无限单元(CINPE4).
本文基于模态分析的结果,确定了模型的Rayleigh阻尼系数,其中前8个振型计算结果如表 3所示.Rayleigh阻尼C=αM +βK (M是质量矩阵,K是刚度矩阵,α和β是阻尼系数).若阻尼比ξ在频率ωi~ωj处保持不变,根据振型的正定性,阻尼参数可以由下式确定:
(10) |
式中,ωi和ωj是对振动模态响应贡献最高(最大的有效质量)的第i和j振型的圆频率.
由表 3可知,垂向分量的1、2振型和水平分量的3、6振型的有效质量最大.将ω1和ω6代入式(10),并令ξ = 5%(Hung et al., 2004; He and Cui, 2015),可得阻尼系数α=0.1538,β= 0.0060.
2.3 最小偏移距的验证本文选用雷克子波作为震源加载到道床上.震源函数为
(11) |
其中A表示位移振幅,t是时间,持续30 ms,t0=15 ms,主频f0=80 Hz.
ABAQUS自适应的稳定时间步长约为2.5×10-2 ms.本文提取了地表上每1 m、共101个节点的位移垂向分量.从地震记录(图 4a)中可以明显的区分出P波、S波和Rayleigh波,结合频谱特征(图 4b)可以直观地发现:近场多高频成分,体波在近场振幅更大,能量更强,且衰减很快;而远场以低频成分为主,Rayleigh波在近场并不发育,随着偏移距增大,Rayleigh波的能量相对于体波逐渐增强.
根据经验公式(6)估计的最小偏移距为15 m.为了验证Rayleigh波最小偏移距计算公式,同时进一步研究最小偏移距对MASW方法数据质量的影响,本文对101道MASW数据进行了抽道选排,计算不同偏移距条件下的频散图,通过对比频散图的分辨率来分析最小偏移距的效应(图 5).利用F-K法(Capon, 1970)可以绘制出频率-相速度(f-v)图谱.当使用0~8 m范围的数据时,图上很显然是没有频散特征的(图 5a);当使用0~12 m的数据时,图上的频散分辨非常差(图 5b).而当偏移距超过15 m后,例如图 5c—5e所示,可以观察到较为明显的频散特征.这与本文提出的最小偏移距的计算公式所得的结果相吻合.另外,考虑到40 m之后的Rayleigh波较为发育,以40~70 m的数据所计算的频散图作为参照(图 5h),通过对比图 5e—5g,可以发现,排列长度相同时(均为25 m),0~25 m、15~40 m和25~50 m的数据所计算的频散图分辨率逐渐增强,越来越接近图 5h所示的稳定频散特征;图 5e中(0~25 m)虽有可被明显识别的频散能量团,但形态与稳定频散特征相差较大,而图 5d(15~40 m)的频散能量团已经与稳定频散特征趋于一致.
综上所述,基于有限元方法获得的数值结果,以15~25 m作为最小偏移距是比较合理的.这能够较好的验证本文所提出的最小偏移距估计方法.同时也说明在MASW法的实际应用中,最小偏移距是一个至关重要的参数,本文提出的最小偏移距的计算公式是比较可靠的.
3 实测分析在有限元分析的基础上,本文在该地铁线旁边进行了基于埋入源的MASW野外实测,按照本文提出的最小偏移距确定原则布设观测系统,对采集到的地震记录进行截取和频散分析,提取频散曲线,反演研究区的一维剪切波速度结构.通过与钻孔资料的对比,验证该方法的可行性.
3.1 野外观测系统与数据预处理本次试验共使用了10个检波器,由西向东呈线性排列,测线垂直于地铁行进方向,如图 6所示.遵循式(6)的理论推导和经验估计,本次试验的最小偏移距应为15 m.考虑到试验现场的空间限制(避开行车道和绿化区)、检波器的数量以及预期的探测深度等主观客观因素,最终决定在距离地铁线路垂直距离20 m处开始放置检波器,即最小偏移距为20 m,道间距为5 m,排列长度为45 m.每道布设3个单轴低频速度传感器(传感器相关指标详见表 4),分别记录x、y和z三个相互垂直的速度分量(x向朝北,y向朝东,z向垂直于地面向下).由于路面交通会对实验数据造成干扰,故选在车流量较小的清晨开展实验.观测小组在最近的地铁站厅中详细记录了地铁的进出站时刻,在测线西侧的人行道上记录了路面车辆的行驶情况,以方便后期能从连续的地震数据中截选出有效的地铁振动信号.
检波器的采样率设为500 Hz,每15 min记录一组数据,实验共记录了3组数据.图 7a展示了其中一组数据的归一化地震波形.由于地铁班次是等时间间隔的,每组数据均包含了6次地铁振动信号(图 7a中Ⅰ—Ⅵ),其他较强振幅的记录是路面交通引起的振动,剩余时段的数据可视为背景噪音.将图 7a的第Ⅰ个地铁振动信号放大后(图 7b)可以看出,地铁引起的振动信号每次大约持续15 s.
由于地铁与测线的相对位置是由远到近再到远的,所以在这15 s的地震记录中,前段和末段数据存在两个主要的问题,可能会导致后续的频散分析失效:
(1) 远距离信号Rayleigh波能量弱,易产生远场效应;
(2) 宽方位角接收的Rayleigh波视速度大于真速度(赵东,2010),低频段的频散曲线与实际不符.
因此,本文设置了一个3 s的时窗(图 7b中虚线框),截取15 s地铁振动信号的中段,切除一部分较远距离和较宽方位角的地震信号,以提取相对稳定可靠的数据进行频散分析.最终本文共提取了18组时长3 s的有效MASW数据集,均为10道.
3.2 频散分析与联合反演本文综合采用F-K法和ESPAC法进行频散分析.逐一对这18组数据进行频散分析后,将得到的频散谱图进行叠加,最终可得两种频散谱图,其中图 8a是F-K的结果,图 8b是ESPAC的结果.由于Rayleigh波信号不是沿测线接收,而是存在一定方位角(虽然做了截取,但无法完全消除),因此基于F-K法得到的频散图,应当从能量团的下沿提取频散曲线;而ESPAC法由于对方位角不敏感,仍可从能量团峰值处提取频散曲线.从两幅频散图上可以看到明显的基阶面波和高阶面波,因此在拾取频散曲线时应注意区分.图 8c和图 8d分别是从F-K频散图和ESPAC频散图中提取的频散曲线,二者区别很小.
本文利用遗传算法,对提取的两组频散曲线进行基阶面波与高阶面波的联合反演(罗银河等,2008).遗传算法反演的参数设置为:迭代次数60次,交叉概率0.90,变异概率为0.02.本文将钻孔资料给出的速度结构作为初始模型(图 8e,蓝线),控制地层深度,通过比较反演结果与钻孔资料的剪切波速,来验证基于埋入源的MASW方法的可行性.通常,MASW所记录的有效波长与排列长度相当(Park et al., 1999b; Martin, 2011).因此速度反演时没有使用波长大于45 m的频散数据.
对F-K法提取的频散曲线进行反演,结果如图 8e中红色折线所示,反演的均方根偏差(RMS Misfit)为2.02%;基于ESPAC法提取的频散曲线的反演结果在图 8e中用绿色折线表示,均方根偏差为1.87%.两种结果的误差都很小.
通常认为,初始值选择对于反演结果的收敛性至关重要(Xia et al., 1999).为了分析反演模型对初始剪切波速的敏感性,本文将初始模型(M)的速度减小10%(模型A)和增大10%(模型B),以F-K法提取的频散曲线为例,分别对模型A和B进行反演(图 9).从图 9可以看到,初始剪切波速增大和减小后,反演结果在视觉上的整体精度是相同的,模型A和B的均方根偏差分别是1.80%和1.69%,能够从两个方向收敛,特别是在中浅层,0~35 m范围内,收敛性很明显,结果较为可信.而在40 m以深的地层收敛性变差,这是低频段频散曲线反演的客观现象,但M、A、B模型各自反演的结果,仍然能够识别出40 m下方的低速变化.
表 5列出了钻孔资料和两种反演结果.从图 8e和表 5可以看出:
(1) 在0~19 m的浅层介质中,两种反演结果都比钻孔资料的速度大,这是由于高频Rayleigh波信号对方位角更加敏感;
(2) 在19 m以深的介质中,反演结果与钻孔资料的速度接近,ESPAC的结果变动幅度相对更大一些;
(3) 反演的层速度相对变化趋势与地层实际速度的变化趋势是一致的:例如从反演结果中可以识别7~15 m和40~45 m的低速层;而在2 m、15 m、19 m、25 m和35 m的正常层界面上,反演的层速度为逐渐增大,与实际速度吻合.
综上所述,本文采用MASW方法反演得到的剪切波速度可以稳定收敛,不会陷入局部极值,虽然反演所得的地铁隧道及邻区的剪切波速在各层介质中存在一定偏差,但反演的均方根偏差较小,且反演的地层波速随深度的变化规律(如低速层、正常层)与钻孔资料是一致的.这说明基于埋入源的MASW方法是可行的,也验证了本文提出的最小偏移距的经验公式.
4 结论多道面波分析(MASW)法以Rayleigh波的频散特征为基础,通过对频散曲线的反演得到地层的剪切波速度结构,被广泛应用于地球物理、土木工程、无损测试等多个领域中.根据前人和本文的研究,认为最小偏移距是关乎MASW方法应用效果好坏的一个重要参数.
本文引入了埋入源,根据波动地震学和几何地震学的基本原理,研究了埋入源激发的Rayleigh波在介质中的传播规律,基于射线理论提出了计算最小偏移距的经验公式x0≈0.75 d,进而为基于埋入源的MASW法应用提供指导.
以上海城市环境下的某地铁振动为研究实例,利用有限元方法进行数值分析,建立了地铁隧道及邻区的地质地球物理模型.通过对不同偏移距的道集的频散分析,验证了最小偏移距公式(x0≈0.75 d)的有效性.
按照最小偏移距的估算方法,在上海市某地铁附近设计了一次现场试验.利用F-K法和ESPAC法进行频散分析,并提取到基阶面波和高阶面波的频散曲线,并利用遗传算法反演得到了研究区的剪切波速结构.结合钻孔资料,发现反演结果能较好地与之吻合.
致谢同济大学唐和生副教授和李晓峰硕士为本文的野外试验提供了帮助,中国科学院上海高等研究院宋学行博士参与了有限元数值建模工作,加拿大骄佳(Geogiga)技术公司为本文的反演提供了软件支持,在此一并对上述专家学者表示感谢.
Aki K. 1957. Space and time spectra of stationary stochastic waves, with special reference to micro tremors. Bull. Earthq. Res. Inst. Tokyo, 35: 415-456. |
Capon J. 1970. Applications of detection and estimation theory to large array seismology. Proc. IEEE, 58(5): 760-770. DOI:10.1109/PROC.1970.7730 |
Chen G X, Xu J L, Yuan C Q. 1998. Relation between depth and shear wave velosity of soil and bedrock in Nanjing City. J. Nanjing Archit. Civ. Eng. Inst., (2): 32-37. |
Coccia S, Del Gaudio V, Venisti N, et al. 2010. Application of Refraction Microtremor (ReMi) technique for determination of 1-D shear wave velocity in a landslide area. J. Appl. Geophys., 71(2-3): 71-89. DOI:10.1016/j.jappgeo.2010.05.001 |
Dai K S, Song X H, Li X F, et al. 2015. Application of the surface wave survey method on multi-scale engineering problems:laboratory and field testing case studies. J. Test. Eval., 43(2): 1-9. DOI:10.1520/JTE20140083 |
Dikmen V, Arısoy M Ö, Akkaya İ. 2010. Offset and linear spread geometry in the MASW method. J. Geophys. Eng., 7(2): 211-222. DOI:10.1088/1742-2132/7/2/S07 |
Ewing W M, Jardetzky W S, Press F. 1957. Elastic Waves in Layered Media. New York: McGraw-Hill, doi: 10.1063/1.3060203.
|
Gao F, Sun X G. 2005. Characteristic analysis of Shear wave velocity of foundation ground in Shanghai region. Shanghai Geology, 26(2): 27-29. DOI:10.3969/j.issn.2095-1329.2005.02.008 |
Gardien W, Stuit H G. 2003. Modelling of soil vibrations from railway tunnels. J Sound Vib., 267(3): 605-619. DOI:10.1016/S0022-460X(03)00727-2 |
Gucunski N, Woods R D. 1992. Numerical simulation of the SASW test. Soil Dyn. Earthq. Eng., 11(4): 213-227. DOI:10.1016/0267-7261(92)90036-D |
Haskell N A. 1953. The dispersion of surface waves on multilayered media. Bull. Seismol. Soc. Am., 43(1): 17-34. |
He P P, Cui Z D. 2015. Dynamic response of a thawing soil around the tunnel under the vibration load of subway. Environ. Earth Sci., 73(5): 2473-2482. DOI:10.1007/s12665-014-3596-8 |
Herrmann R B. 2013. Computer programs in seismology:An evolving tool for instruction and research. Seismol. Res. Lett., 84(6): 1081-1088. DOI:10.1785/0220110096 |
Hung H H, Yang Y B, Chang D W. 2004. Wave barriers for reduction of train-induced vibrations in soils. J Geotech Geoenviron Eng., 130(12): 1283-1291. DOI:10.1061/(ASCE)1090-0241(2004)130:12(1283) |
Jones R. 1958. In-situ measurement of the dynamic properties of soil by vibration methods. Géotechnique, 8(1): 1-21. DOI:10.1680/geot.1958.8.1.1 |
Kouroussis G, Verlinden O, Conti C. 2011. Finite-dynamic model for infinite media:corrected solution of viscous boundary efficiency. J. Eng. Mech., 137(7): 509-511. DOI:10.1061/(ASCE)EM.1943-7889.0000250 |
Lai C G, Callerio A, Faccioli E, et al. 2005. Prediction of railway-induced ground vibrations in tunnels. J. Vib. Acoust., 127(5): 503-514. DOI:10.1115/1.201330 |
Lay T, Wallace T C. 1995. Modern Global Seismology. California: Academic Press.
|
Lin C P, Chang C C, Chang T S. 2004. The use of MASW method in the assessment of soil liquefaction potential. Soil Dyn. Earthq. Eng., 24(9-10): 689-698. DOI:10.1016/j.soildyn.2004.06.012 |
Liu Q H, Lu L Y, Wang K M. 2015. Review on the active and passive surface wave exploration method for the near-surface structure. Progress in Geophys., 30(6): 2906-2922. DOI:10.6038/pg20150660 |
Liu Y Z, Wang Z D. 1996. Data collection and processing system of transient surface wave method and examples of its application. Geophys. Geochem. Explor., 20(1): 28-34. |
Lou M L, Jia X P, Yu J Q. 2009. Field measurement and analysis of ground vibration induced by subway trains. J. Disaster Prev. Mitig. Eng., 29(3): 282-288. |
Louie J N. 2001. Faster, better:shear-wave velocity to 100 meters depth from refraction microtremor arrays. Bull. Seismol. Soc. Am., 91(2): 347-364. DOI:10.1785/0120000098 |
Luo G C, Liu L B, Qi C, et al. 2011. Structural response analysis of a reinforced concrete building based on excitation of microtremors and passing subway trains. Chinese J. Geophys., 54(10): 2708-2715. DOI:10.3969/j.issn.0001-5733.2011.10.028 |
Luo Y H, Xia J H, Liu J P, et al. 2008. Joint inversion of fundamental and higher mode Rayleigh waves. Chinese J. Geophys., 51(1): 242-249. |
Martin A J. 2011. Analysis of SASW, MASW, and passive surface wave data collected at the national geotechnical experimentation site at Texas A & M University. //Proceedings of the Geo-Risk 2011: Risk Assessment and Management. Atlanta, GA, 886-893, doi: 10.1061/41183(418)94
|
Miller R D, Xia J H, Park C B, et al. 1999. Multichannel analysis of surface waves to map bedrock. The Leading Edge, 18(12): 1392-1396. DOI:10.1190/1.1438226 |
Nakamura Y. 1989. A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Quart. Rep. Railw. Tech. Res. Inst., 30(1): 25-33. |
Nastev M, Parent M, Ross M, et al. 2016. Geospatial modelling of shear-wave velocity and fundamental site period of Quaternary marine and glacial sediments in the Ottawa and St. Lawrence Valleys, Canada. Soil Dyn. Earthq. Eng., 85: 103-116. DOI:10.1016/j.soildyn.2016.03.006 |
Okada H, Suto K. 2003. The Microtremor Survey Method. Geophysical Monograph Series No. 12. Tulsa, Oklahoma: Society of Exploration Geophysicists.
|
Park C B, Miller R D, Xia J H. 1999a. Multichannel analysis of surface waves. Geophysics, 64(3): 800-808. DOI:10.1190/1.1444590 |
Park C B, Miller R D, Xia J H, et al. 1999b. Higher mode observation by the MASW method. //Proceedings of the SEG Technical Program Expanded Abstracts. Society of Exploration Geophysicists, 524-527, doi: 10.1190/1.1821070.
|
Park C B, Miller R D, Xia J H. 2001. Offset and resolution of dispersion curve in multichannel analysis of surface waves (MASW). //Proceedings of the Symposium on the Application of Geophysics to Engineering and Environmental Problems (SAGEEP). Environment and Engineering Geophysical Society, SSM4, doi: 10.4133/1.2922953.
|
Park C B, Miller R D. 2008. Roadside passive multichannel analysis of surface waves (MASW). J. Environ. Eng. Geophys., 13(1): 1-11. DOI:10.2113/JEEG13.1.1 |
Peng Y J, Lv Y J, Shi C H, et al. 2011. Site classification based on Quaternary Geology and shear wave velocity in Beijing plain region. //Proceedings of the 1st National Conference on Disaster Prevention and Mitigation Engineering(in Chinese). Nanjing: China Association for Disaster Prevention, 108-114.
|
Qi W H, Liu D D, Lan J Y, et al. 2008. Statistical analysis of soil layers' shear wave in Yanliang, Xi'an. J. Inst. Disaster Prev. Sci. Technol., 10(4): 10-12. DOI:10.3969/j.issn.1673-8047.2008.04.003 |
Rayleigh L. 1885. On waves propagated along the plane surface of an elastic solid. Proc. Lond. Math. Soc., 17(1): 4-11. |
Sheriff R E, Geldart L P. 1995. Exploration seismology. 2nd ed. New York: Cambridge University Press.
|
Stokoe K H, Nazarian S. 1983. Effectiveness of ground improvement from spectral analysis of surface waves. //Proceeding of the 8th European Conference on Soil Mechanics and Foundation Engineering. Helsinki, Finland, 91-94.
|
Stokoe K H, Wright S G, Bay J A, et al. 1994. Characterization of geotechnical sites by SASW method. //Geophysical Characterization of Sites. New Delhi: Oxford Publishers, 15-25.
|
Tian G, Steeples D W, Xia J H, et al. 2003. Multichannel analysis of surface wave method with the autojuggie. Soil Dyn. Earthqu. Eng., 23(3): 61-65. DOI:10.1016/S0267-7261(02)00214-2 |
Wang M Y, Guo Y X. 1987. Finite element modelling of elastic wave propagating in 2-D media and its preliminary applications. Acta Geophysica Sinica, 30(3): 292-306. |
Wang S Y, Wang H Y. 2016. Site-dependent shear-wave velocity equations versus depth in California and Japan. Soil Dyn. Earthq. Eng., 88: 8-14. DOI:10.1016/j.soildyn.2016.05.001 |
Xia J H, Miller R D, Park C B. 1999. Estimation of near-surface shear-wave velocity by inversion of Rayleigh waves. Geophysics, 64(3): 691-700. DOI:10.1190/1.1444578 |
Xia J H, Gao L L, Pan Y D, et al. 2015. New findings in high-frequency surface wave method. Chinese J. Geophys., 58(8): 2591-2605. DOI:10.6038/cjg20150801 |
Xu P F, Li C J, Ling S Q, et al. 2009. Mapping collapsed columns in coal mines utilizing Microtremor Survey Methods. Chinese J. Geophys., 52(7): 1923-1930. DOI:10.3969/j.issn.0001-5733.2009.07.028 |
Xu Y X, Xia J H, Miller R D. 2006. Quantitative estimation of minimum offset for multichannel surface-wave survey with actively exciting source. J. Appl. Geophys., 59(2): 117-125. DOI:10.1016/j.jappgeo.2005.08.002 |
Yang L Z, He X L, Chi Q L. 2011. Vibration measurements of city environment. Progress in Geophys., 26(1): 380-388. DOI:10.3969/j.issn.1004-2903.2011.01.047 |
Yang Y B, Hung H H, Chang D W. 2003. Train-induced wave propagation in layered soils using finite/infinite element simulation. Soil Dyn. Earthq. Eng., 23(4): 263-278. DOI:10.1016/S0267-7261(03)00003-4 |
Yilmaz O, Eser M. 2002. A unified workflow for engineering seismology. //SEG Technical Program Expanded Abstracts 2002. Salt Lake City, UT: SEG, 1496-1499, doi: 10.1190/1.1816948.
|
Zerwer A, Cascante G, Hutchinson J. 2002. Parameter estimation in finite element simulations of Rayleigh waves. J. Geotech. Geoenviron. Eng., 128(3): 250-261. DOI:10.1061/(ASCE)1090-0241(2002)128:3(250) |
Zhan J Y, Chen G X, Liu J D. 2009. Empirical relationship between shear wave velocity and soil depth on deep soft sites in urban area of Suzhou city. World Earthq. Eng., 25(2): 11-17. |
Zhang B L, Li Z W, Bao F, et al. 2016. Shallow shear-wave velocity structures under the Weishan volcanic cone in Wudalianchi volcano field by microtremor survey. Chinese J. Geophys., 59(10): 3662-3673. DOI:10.6038/cjg20161013 |
Zhang S X, Chan L S, Xia J H. 2004. The selection of field acquisition parameters for dispersion images from multichannel surface wave data. Pure Appl. Geophys., 161: 185-201. DOI:10.1007/s00024-003-2428-7 |
Zhang Z L. 2007. Statistical analysis of shear wave velocity in construction sites [Master′s thesis]. Beijing: Institute of Engineering Mechanics, China Earthquake Administration.
|
Zhao D. 2010. Passive surface waves:methods and applications. Geophys. Geochem. Explor., 34(6): 759-764. |
陈国兴, 徐建龙, 袁灿勤. 1998. 南京城区岩土体剪切波速与土层深度的关系. 南京建筑工程学院学报, (2): 32-37. |
高飞, 孙小刚. 2005. 上海地区场地剪切波速的特征分析. 上海国土资源, 26(2): 27-29. DOI:10.3969/j.issn.2095-1329.2005.02.008 |
刘庆华, 鲁来玉, 王凯明. 2015. 主动源和被动源面波浅勘方法综述. 地球物理学进展, 30(6): 2906-2922. DOI:10.6038/pg20150660 |
刘云祯, 王振东. 1996. 瞬态面波法的数据采集处理系统及其应用实例. 物探与化探, 20(1): 28-34. |
楼梦麟, 贾旭鹏, 俞洁勤. 2009. 地铁运行引起的地面振动实测及传播规律分析. 防灾减灾工程学报, 29(3): 282-288. |
罗桂纯, 刘澜波, 齐诚, 等. 2011. 基于地脉动和地铁振动的钢筋混凝土建筑结构响应分析. 地球物理学报, 54(10): 2708-2715. DOI:10.3969/j.issn.0001-5733.2011.10.028 |
罗银河, 夏江海, 刘江平, 等. 2008. 基阶与高阶瑞利波联合反演研究. 地球物理学报, 51(1): 242-249. |
彭艳菊, 吕悅军, 施春花等. (2011). 北京平原区基于第四纪地质与剪切波速的场地类别划分. //全国第一届防灾减灾工程学术研讨会论文集. 南京: 中国灾害防御协会, 108-114.
|
齐文浩, 刘德东, 兰景岩, 等. 2008. 西安阎良区土层剪切波速统计分析. 防灾科技学院学报, 10(4): 10-12. DOI:10.3969/j.issn.1673-8047.2008.04.003 |
王妙月, 郭亚曦. 1987. 二维弹性波的有限元模拟及其初步实践. 地球物理学报, 30(3): 292-306. |
夏江海, 高玲利, 潘雨迪, 等. 2015. 高频面波方法的若干新进展. 地球物理学报, 58(8): 2591-2605. DOI:10.6038/cjg20150801 |
徐佩芬, 李传金, 凌甦群, 等. 2009. 利用微动勘察方法探测煤矿陷落柱. 地球物理学报, 52(7): 1923-1930. DOI:10.3969/j.issn.0001-5733.2009.07.028 |
杨立志, 何先龙, 匙庆磊. 2011. 城市环境振动的测量. 地球物理学进展, 26(1): 380-388. DOI:10.3969/j.issn.1004-2903.2011.01.047 |
战吉艳, 陈国兴, 刘建达. 2009. 苏州城区深软场地土剪切波速与土层深度的经验关系. 世界地震工程, 25(2): 11-17. |
张宝龙, 李志伟, 包丰, 等. 2016. 基于微动方法研究五大连池火山区尾山火山锥浅层剪切波速度结构. 地球物理学报, 59(10): 3662-3673. DOI:10.6038/cjg20161013 |
张忠利. 2007. 建设场地剪切波速的统计分析[硕士论文]. 北京: 中国地震局工程力学研究所.
|
赵东. 2010. 被动源面波勘探方法与应用. 物探与化探, 34(6): 759-764. |