面波法勘探在工程勘察中的应用Word格式.doc
- 文档编号:8636606
- 上传时间:2023-05-12
- 格式:DOC
- 页数:18
- 大小:712KB
面波法勘探在工程勘察中的应用Word格式.doc
《面波法勘探在工程勘察中的应用Word格式.doc》由会员分享,可在线阅读,更多相关《面波法勘探在工程勘察中的应用Word格式.doc(18页珍藏版)》请在冰点文库上搜索。
第一章地震面波简介
地震波是地震震源在地球介质中产生的扰动。
在有介质分界面存在时,地震波除了像反射波和折射波那样在整个介质体内传播的体波外,还存在一类沿介质自由界面传播的面波,当它沿着自由表面传播时,其能量主要集中在自由表面附近,并随着深度的增加能量迅速衰减。
面波按其类型主要有瑞雷面波和勒夫面波两大类。
瑞雷面波是1887年由英国学者瑞雷首先在理论上确定的,这种面波分布在自由表面上,或者表面为疏松的覆盖层内。
当介质为均匀各向同性介质时,瑞雷面波的相速度和群速度将一致,否则瑞雷波的相速度将不一致,出现频散现象,当介质具有水平层状性质时,瑞雷面波的频散规律与介质的分层结构紧密相关。
瑞雷面波既有P波成分也有S波成分,而无S波成分。
瑞雷面波在天然地震中常常可以观测到,它对建筑物的破坏性极大。
在地震勘探中,瑞雷波已由过去的干扰波变成了可以利用的信号。
勒夫面波产生于介质表面的低速覆盖层以及该层与下面介质的分界面上。
勒夫面波面波是一种S型波,具有频散现象。
假定存在一均匀完全弹性的半无限空间,不均匀平面纵波与不均匀平面横波沿自由表面传播时相互叠加就产生了瑞利面波。
在各向均匀半无限空间弹性介质表面上,当一个圆形基础上下运动时,由它产生的弹性波入射能量的分配率已由Miller(1955年)计算出来,即P波占7%、S波占26%、R波占67%,亦就是说,R波的能量占全部激振能量的2/3,因此利用面波作为勘探方法,其信噪比会大大提高。
第二章瑞利波勘察原理及现场工作方法
2.1瑞利波勘察原理
瑞利波沿地表面传播,其穿透能力仅有一个波长,也就是说,可以达到距表层一个波长λ的深度范围。
如果能在水平方向的测线上记录同一波长不同点的V值,就可以反应地质界面在水平方向的变化特征。
若记录不同λ的V值,也就可以反映出不同深度的地层分布和特征。
瑞利波和反、折射波一样都是沿测线方向传播的。
在测线上以一定道间距ΔX设置N+1个检波器,就可以观测到瑞利波在NΔX长度范围内传播的过程。
设瑞利波的频率为fi,相邻两各检波器的瑞利波的到时差为Δt或相位差为Δφ,则相邻两道ΔX长度的范围内,瑞利波的传播速度,可以记为:
(2-1)式中(2-2)
测量范围NΔX内地层
的平均速度为
(2-3)
在同一地段测量出一定频率的值就可以得到一条(V-f)曲线,即所谓的频散特性曲线或把(V-f)曲线转化为(V-λ)曲线,可用下式表示:
(2-4)
由于(V-f)与(V-λ)曲线的变化规律与地层地质条件存在着内在的联系,因此通过对频散曲线的反演解释,可以得到地下某一定深度范围内的地质构造,也可以得到不同深度地层的V值。
2.2多道瞬态面波数据采集方法
2.2.1仪器及配件
一套完整的多道瞬态面波采集系统至少应该配备以下仪器和配件:
地震仪:
用于处理和存贮地震波信号。
一般使用SWS多通道工程地震仪,也可以采用其它通用多通道数字地震仪,数据通道不应低于六道。
数据传输线:
用于检波器和地震仪间的数据传输,其长度不应小于最大测线长度。
检波器:
用于接收地震波信号,面波采集时应采用低频检波器。
触发开关:
触发开关通过导线连接震源和地震仪,以保证震源的激发的同时地震开始记录数据,使所采集的地震信号具有时间特性。
震源:
一般的浅层面波勘探常采用锤击震源,也可采用落重或炸药震源。
电源:
根据不同的仪器要求,配备相适应电源。
2.2.2数据采集
1、侧线布置
野外数据采集时使用低频面波检波器在震源纵向方向等间距排列,如图2-1所示,排列长度应大于预期探测深度,排列线附近地面尽量避免有沟、坎、墙等能产生反射或散射的障碍物。
图2-1
2、参数设置
仪器开启进入面波采集系统后,需要对以下各种参数进行设置:
存盘路径:
用来指定数据的存放地点,方便以后调用。
文件名:
由于野外数据采集时往往数据量大,文件多,所以必须正确设置文件名,否则极易搞混淆。
文件名应该包括代号和代码两部分,代号部分一般用工程名称的拼音简写,代码部分用来表示数据采集的先后顺序,可以由仪器自动生成。
采样间隔:
常采用0.20或0.25ms
每道采样数:
常取1024、2048、4096等。
道数:
根据实际情况设置,常采用12道或24道。
道间距:
道间距由测线长度和道数控制,设测线长度为L,道数为n,则道间距为L/(n-1)。
道间距的设置还应该考虑分辨率要求,不能大于欲探测的最薄地层的厚度。
偏移距:
视具体情况而定,取值范围一般为2~10米。
3、震源激发
多道瞬态面波震源激发位置必须位于检波器排列的纵向方向,可置于前端也可置于后端,最小偏移距不宜低于2m。
震源能量视预期勘探深度而定,当预期勘探深度小于30m时,可用人工锤击震源,预期勘探深度在30~80m可采用落重
震源,预期勘探深度大于80m时一般应使用炸药震源。
瞬态面波测深要求采用的震源在时间上是单个脉冲的冲击。
在锤击或落重操作中往往会产生连击,甚至在爆炸时,由于围岩的影响,也能出现反冲。
如果两个脉冲的时间间隔小于期望获得的面波最长周期,就不可能用时间-空间窗口加以清除,而会在频率波数谱上出现周期性的能量强弱起伏,严重时甚至会导致相位的周期性扭曲。
4、数据检查与保存
当震源激发后,地震仪会将所接受到的地震波形记录显示在屏幕上,可以通过增益控制键调整波形幅度,观察有无缺道以及干扰的大小等,然后确定记录信号的质量是否合乎要求。
有的工区,由于客观条件的限制而无法避开干扰,这时就需要采用多次叠加技术来压制干扰,叠加的次数视具体情况而定。
当确定接收的信号合乎要求后,即可存盘。
第三章瑞利波资料整理与解释
3.1面波频散曲线的深度解释
要利用面波频散曲线进行地层划分,首先要确定面波波长与深度的转换系数β,以便将面波f-V曲线转换为H-V曲线。
瑞雷波的能量随深度按指数规律衰减,通常定义当振幅比µ
/µ
衰减到1/e时的深度为穿透深度,其中µ
为横向振动的振幅,µ
为纵向振动的振幅。
针对不同的岩土介质,我们可以计算出穿透深度与振幅能量之间的关系,如表3-1所示,从而确定出比较合理的深度H与波长λ的转换系数值β值。
从表3-1中可以看出,对于所有的介质,瑞雷波的穿透深度为0.55λ~0.875λ。
对于土体而言,泊松比σ=0.4-0.45,则穿透深度H≈(0.79-0.84)λ。
对于淤泥质软塑土层,穿透深度可取0.85λ。
对于一般土层穿透深度可采用:
(3-1)
实际应用中,由于各测区地层条件一般不会相同,所以应该根据现场对比试验来确定合适的深度的转换系数β。
一般来说,以上述β值绘出V—βλ曲线中的传播速度能够代表βλ深度以上的平均速度,其变化规律与V—λ曲线一致。
表3-1不同介质中瑞利波的穿透深度
3.2层厚度的计算方法
在实际勘察工作中,以V为横坐标,以H=βλ为纵坐标,绘制V-H曲线(如图3-2),曲线的纵坐标就可近似代表勘探深度。
分析V-H曲线的形态和变化规律,可以初步确定地层界面深度以及各层速度的大概范围。
精确确定地层段划分主要有以下两种方法:
图3-2
1、一阶导数极值点法
根据曲线的极值点对应的分层位置,求出其波长,并根据确定分层深度。
2、拐点法
根据V-βλ曲线上的拐点的位置,计算出拐点处所对应的波长,同样根据确定出分层深度。
需要注意的是瑞利波速度代表着深度以上介质的平均速度。
对于多层介质,深度计算公式需要作适当的修正,即:
,以消除层间的影响。
3.3层速度的计算方法
1、瑞利波速度
层速度计算一般采用近似计算方法,即近似的认为瑞雷波传播速度代表某一深度内各层波速的加权平均值。
分层速度由下式计算:
速度随H增大而增大时:
(3-2)
速度随深度减小时:
(3-3)
式中:
为第n点深度;
为第n-1点深度;
为n点深度以上平均波速;
为第n-1点深度以上平均波速;
为—深度间隔的波速。
2、横波速度
由瑞雷波的基本原理可知,在均匀各项同性半无限弹性介质中,瑞雷波速
度和横波速度的近似关系为。
对于不同的岩土介质,经理论计算,可得到不同泊松比σ与的关系见表3-3。
从表中可以看出,随着泊松比的增大,相对急剧增大,而与的值则趋于一致。
一般岩石泊松比在0.25左右,第四系地层泊松比为0.4~0.49,可以认为对土体而言,与基本相等。
其误差只有5%左右,即。
由此,在计算岩土的力学参数时,可利用代替进行近似计算。
表3-3瑞利波速度与横波速度的关系表
第四章工程实例
4.1工程概述
2007年春季,陕西延长县某河岸黄土陡边坡出现滑动迹象,严重威胁穿行而过的输气管道的安全。
滑坡区位于陕北黄土高原中部地区,属于典型的继承型和继承、侵蚀混合型的黄土沟壑地貌。
黄土广覆于下伏起伏的基岩古地形之上,在新构造运动的影响下,由于河流的长期侵蚀切割,形成沟谷深切的地形。
该处滑坡受河流的侵蚀切割作用,坡体前部形成陡坡临空面,为滑坡提供了滑动的空间条件。
该滑坡危及输气管线的安全,所以必须尽快勘察和治理。
由于时间紧迫,按常规的勘察方法很难满足工期要求,急需一种快速、高效、准确的勘探方法。
经分析和研究,最后确定了以多道瞬态面波勘探为主,配以少量钻孔加以验证的勘察方案。
面波勘察的任务是要划分地层并确定滑动面的位置。
4.2数据采集和处理
1、野外数据采集
仪器采用北京水电物探研究所SWS工程地震仪,检波器使用4赫兹面波检波器,震源采用18磅铁锤锤击合金垫板。
采集时参数设置如下:
0.25ms
2048
12道
2m
各点取4m和8m各采集一组数据
观测系统
测试中沿河堤走向布置测线,等间距设置6个测点(CD1~CD6),总覆盖长度约140m。
2、数据分析和处理
对所采集的数据按本文所述方法在时间-空间域和频率-波数域进行初步处理和分析,发现测区面波数据呈现两类不同的特征:
一类是基阶面波型号很强,高阶信号较弱,,如图4-1所示,类似Ⅰ类(波速由表层向底层逐层增高的横向均匀介质)地层面波频散特性。
二类是基阶面波和高阶面波都很强,但是能够区分,如图4-2所示,类似Ⅱ类(底层波速最高,中部含低速层的横向均匀介质)地层面波频散特性。
通过进一步分析发现,测区面波频散曲线形状非常类似,都具有两处明显拐折,较深处拐折出现速度倒转现象,说明该深度处具有低速夹层。
而测区部分面波信号(CD3-CD6)在时间-空间域和频率-波数域出现类似Ⅰ类地层的面波频散特征,则是由于该处低速层深度较大,减弱了高速波导效应。
事实上,我们从图上可以看出,高阶面波还是具有一定能量,要比正常Ⅰ类地层中高阶面波信号稍强。
由此可以得出,测区内地层都可以归为Ⅱ类地层。
将验证孔的地质和面波资料进行对比分析,经多次计算和比较,确定β取0.8最为合适。
按此换算系数得出的面波H-VR曲线上两处拐折很好的对应了两处地层分界面,如图4-3所示。
由此也确定了层段划分标准。
图4-1测点4实测面波S-T域波形记录及(右)F-K域图谱
图4-2测点1实测面波S-T域波形记录及(右)F-K域图谱
对所有资料进行时间空间域和频率-波数域处理后得到各测点基阶面波频散曲线,然后按上面得出的波长-深度转换系数将面波频散曲线转换为H-VR曲线。
各曲线按其测点相对位置显示在图4-4中。
4.3地层划分及滑动面确定
利用所得VR-H曲线反演横波速度结构,反演方法采用人机联合反演。
首先根据频散曲线的形状特点确定地层界面位置,然后再反演其横波速度结构。
反演结果见表4-1。
表4-1测点1~6横波速度反演结果表
CD1
CD2
CD3
层号
深度(s)
(m/s)
1
1.66
97
1.76
130
1.63
120
2
6.96
166
6.72
173
7.38
165
3
10.03
256
17.92
227
17.14
4
12.17
152
21.12
150
20.56
138
5
13.84
272
24.82
269
26.78
324
6
18.49
409
29.81
427
41.65
779
CD4
CD5
CD6
3.8
149
1.79
136
3.20
148
8.8
188
11.63
174
8.74
177
13.7
207
23.84
263
11.92
159
20.3
252
26.86
169
17.10
249
23.6
199
31.12
403
21.15
198
28.8
290
40.69
762
31.02
7
38.8
456
40.29
542
可以看出,地层横波速度总体随深度加深呈增大趋势,深度在20m左右出现低速夹层,横波速度值不超过200m/s。
归纳各点横波速度结构并结合验证孔地层结构情况,确定测线下方从上向下依次分布四层地层:
1、地表至8米深度左右,横波速度97~207(m/s),为黄土:
浅黄色,结构疏松,孔隙发育,天然状态下垂直节理十分发育,遇水显湿陷性,强度很低。
该层在陡坡处出现相对滑动。
2、黄土层下方,波速范围227~263(m/s),为黄土质砂粘土:
黄褐色,结构致密,垂直节理发育,无层理,无湿陷性,抗剪强度较高。
3、黄土质砂粘土下方,横波速度降低到138~198(m/s),为泥岩:
灰白色,强风化,多为土状,位于地下水位以下,遇水软化,抗剪强度极低。
该层即为滑动带。
4、最底层的横波速度明显增大,超过400(m/s),为砂岩:
灰色-灰白色,坚硬,抗剪强度高,顶部具有3-5米的风化带,具有垂直节理,但其空隙和节理中往往充填了粘土物质而形成相对的隔水层,易于形成滑坡面。
图4-5和图4-6为根据以上反演和解释所得的地层剖面图。
图4-5为该滑坡沿河堤走向的地层解释剖面,图4-6为该滑坡沿滑坡轴方向的地层解释剖面。
瞬态面波方法在该滑坡勘察中的应用取得了良好的效果:
为滑坡治理方案的设计提供了地层分层参数;
节约了勘察成本,缩短了勘察周期,为滑坡的治理赢得了宝贵的时间;
采用面波勘探方法,避免了大量钻探工作对土体的扰动而加速滑坡体的滑动。
图4-3实测面波频散曲线与钻孔柱资料对比图
图4-4测点1~6的实测面波频散曲线(β=0.8)
图4-5沿河堤走向的地层解释剖面
图4-6沿滑坡轴方向的地层解释剖面
第五章结论
本文结合工程实例,对瑞利波勘察的方法可以得出以下结论和认识:
1、瞬态瑞雷波法用于浅层地质勘探和工程检测具有高效,环保,经济,抗干扰强等优点,在中间夹软弱低速层的复杂地层条件下也能使用,能够弥补传勘探方法的不足。
2、多道瞬态面波数据处理过程中,采用频率-波数谱分析法求取频散曲线效率高,结果很稳定,但其只能反映地层在横向变化的平均效应。
3、不同结构类型地层的面波频散特征不同,相应的面波采集处理方法也应不一样。
4、在实际面波勘探工作中,应该根据测区面波资料与地质资料对比,确定合适的瑞雷波长与深度转换系数。
对于一般的黄土覆盖层,波长深度转换系数取0.8比较合适。
参考文献
[1]朱介寿等.地震学中的计算方法[M].北京:
地震出版社,1988.
[2]李录明,李正文.地震勘探原理、方法及解释[M].成都:
成都理工大学内部教材,2005.
[3]杨文采.地球物理反演的理论与方法[M].北京:
地质出版社,1996.
[4]袁学诚.中国地球物理图集[M].北京:
[5]左兆荣,房明山.震相与地球内部结构[M].北京:
地质出版社,1986.
[6]傅淑芳,刘宝诚.地震学教程[M].北京:
地震出版社,1991.
[7]曾融生.固体地球物理学导论[M].北京 科学出版社,1984.
[8]雷宛,肖宏跃,邓一谦.工程与环境物探教程[M].北京:
地质出版社,2006.
[9]杨成林等,瑞雷波勘探[M].1993年,地质出版社。
[10]吴如山,安艺敬一.主编.地震波的散射与衰减[M],李裕澈,卢寿德等译,吴建春,刘洪校,北京:
地震出版社,1993.
[11]安艺敬一,P.G.理查兹.定量地震学[M],李钦祖,邹其嘉等译,傅承义校,北京:
[12]G.Nolet,edit,SeismicTomography[M].Reidel,Dordrecht,1987.
[13]朱介寿等,用瑞雷面波研究东亚及西太平洋地壳上地幔速维结构[J].物探化探计算技术2005,27(3):
185-192.
[14]朱介寿,曹家敏,蔡学林等.东亚及西太平洋边缘海高分辨率面波层析成像[J].地球物理学报,2002,45(5):
646-664.
[15]黄绪德.由面波求表层横波速度[J].勘探地球物理进展[J].2004,27
(1),9-15.
[16]徐华全,顾勤平,刘远志.多道瞬态面波法在覆盖层调查中的应用[J].内蒙古石油化工.2008,17:
26-27.
[17]董海文.瑞雷面波法在公路质量检测中的理论与应用研究[D].长沙:
湖南大学土木工程学院,2005.
[18]李晶.面波在地震波场中的特性研究及应用[D].成都:
成都理工大学信息工程学院,2006.
[19]贾辉.瞬态瑞雷波法在隐伏断层探测中的应用研究[D].北京:
中国地震局地球物理研究所,2007.
[20]肖伯勋.高模式瑞雷面波及其正反演研究[D].长沙:
中南大学,2000.
[21]娄晓龙.瑞雷面波法在路基压实度检测中的试验研究[J].公路交通科技.2006,5:
66-68.
-16-
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 面波法 勘探 工程 勘察 中的 应用
![提示](https://static.bingdoc.com/images/bang_tan.gif)