【摘要】：面波,通常指瑞雷波(Rayleigh wave)和勒夫波(Love wave),具有低速度、低频率、高振幅以及频散(相速度随频率变化)的特点。瑞雷波由P波和SV波相互干涉形成,它同时具有纵向(振动方向平行于波的传播方向)及横向粒子运动(振动方向垂直于波的传播方向)；勒夫波由不同形式的SH波干涉而形成,它仅包含横向粒子运动。在浅地表地震炮集记录中,高频面波信号占据了浅地表地震波场能量的主要部分,并包含有丰富的地下剪切波(S波)速度信息。基于对高频面波信号的分析,高频面波(包括瑞雷波和勒夫波)方法已广泛应用于求取浅地表S波速度结构。通过对多道地震信号的处理与分析,面波(包括瑞雷波和勒夫波)多道分析方法(Multichannel Analysis of Surface Waves, MASW,这里的Surface Waves代指瑞雷波;Multichannel Analysis of Love Waves, MALW)是目前最为广泛应用的高频面波方法。由于具有无损、高精度及高效的特点,面波多道分析方法越来越受到浅地表地球物理及地质工程学界的重视。瑞雷波和勒夫波多道分析方法具有相似的数据处理及反演步骤,具体包括：多道瑞雷波和/或勒夫波野外数据采集；野外数据向频率速度(f-v)域的变换；根据连续的频散能量峰值拾取面波相速度(又称为频散曲线)；反演拾取的相速度获得地下一维S波速度结构；重复前面的工作获得拟二维S波速度剖面图。大量的实例证明了面波多道分析方法的有效性。然而,在遇到一些较为复杂的地球模型时(例如含有低速夹层,断层,空洞等),现行的面波多道分析方法会遇到问题而无法给出准确地下速度结构。本博士学位论文的主要目的是提高高频面波方法对任意非均匀二维模型的成像能力,具体包括改进现行面波多道分析方法中频散曲线的正演算法,使其适用于含低速层的地球模型；改进勒夫波多道分析方法的野外采集系统,尝试新的勒夫波激发方式,并将勒夫波多道分析方法推广至三维数据采集系统；提出反演面波波形获取地下S波速度结构。计算地球模型对应的理论频散曲线是面波多道分析方法(MASW)的一个重要基础。然而,当地下出现低速层时,用传统方法计算瑞雷波频散曲线会出现问题,该问题具体可分为以下两类：当半空间为低速层时,由于它的频散方程中出现了复数矩阵,使其无法在实数域中获得方程的根,从而难以在某些频率范围内计算获得面波理论相速度值；当地层中含有低速夹层时且该层的S波速度低于表层S波速度时,计算得到的频散曲线在其高频部分会趋近该低速夹层的S波速度而非某个与表层S波速度相关的值。然而,基于波动方程的数值模拟结果显示其频散能量的高频部分趋近于表层S波速度的大约91%,即现行的频散曲线计算方法对于该类含有低速夹层的模型无法在高频范围获得正确的相速度。对于第一类问题,本文提出利用频散方程的实部求解瑞雷波相速度；对于第二类问题,本文基于瑞雷波的穿透深度对现行算法进行改进：我们建立一个仅含低速夹层及其上方地层的替代模型,在计算频散曲线过程中,当瑞雷波波长无法穿透该低速夹层时,我们采用替代模型计算其频散曲线。本文通过基于波动方程的数值模拟方法在频率-速度域中验证了这两种解决方案的正确性,通过数值模拟反演算例证明了这两类改进的频散计算方案在频散曲线反演过程中的稳定性。勒夫波仅由SH波形成,而与P波无关,这使得勒夫波频散曲线比瑞雷波更加简单,从而使勒夫波多道分析方法成为一种很有吸引力的获取浅地表S波速度的方法。在勒夫波多道分析方法中通常采用横向震源激发勒夫波信号。然而,横向震源相比垂向震源而言激发难度更大,且不适用于硬地表等实际情况,在一定程度上限制了勒夫波多道分析方法的适用范围。本文提出在勒夫波多道分析方法中采用垂向震源激发勒夫波。本文对垂向震源与横向震源激发勒夫波的效果进行了四组对比试验。通过对比勒夫波波形,野外数据振幅及其频散能量图,本文发现垂向震源可以激发出与横向震源能量相当、质量相似的勒夫波信号。这为将来在采用勒夫波多道分析方法时的震源方式提供了新的选择,拓宽了勒夫波多道分析方法的适用范围。现行的勒夫波多道分析方法基于二维地震采集系统(检波器与震源沿同一测线排列),无法适用于三维地震采集系统。由于勒夫波的粒子振动方向平行于地表且与传播方向垂直,这使得我们很难在三维地震采集系统中直接观测和采集到勒夫波信号。本文提出在地震数据采集时,在每个检波点同时采集两个相互垂直的水平分量(沿线方向和垂线方向水平分量)。通过这两个分量的旋转与合成,我们分别恢复每个采集点的勒夫波分量。在将二维空间域的勒夫波信号按照真实炮检距转换到一维空间域后,我们可以对重获的勒夫波信号进行频散分析及频散曲线反演,从而获得地下S波速度结构。本文采用数值模拟算例和一组实测数据证明了该方法的正确性。瑞雷波和勒夫波多道分析方法(MASW和MALW)均是基于频散曲线的反演获得地下S波速度结构。由于频散曲线正演基于水平层状模型假设,瑞雷波和勒夫波多道分析方法均仅适用于水平层状模型及横向S波速度缓变的地球模型。这被视为目前高频面波方法发展的一个重要局限。另一个目前高频面波方法遇到的问题是频散曲线及其模式的正确识别,特别是对含有软夹层、起伏地表等较为复杂的地球模型。波形反演直接拟合实测地震波波形,不需要拾取及反演频散曲线,适用于大多数地球模型。波形反演不受层状模型假设约束,同时也不依赖于频散曲线的拾取,而从突破了现行高频方法的局限及困难。本文提出在时间域通过反演面波(包括瑞雷波和勒夫波)波形获取浅地表S波速度结构。本文通过检测板模型测试了面波对地下介质的分辨率,并根据该分辨率将反演网格由浅至深依次增大,提出变反演网格。本文采用有限差分作为面波(包括瑞雷波和勒夫波)波形正演计算方法。通过共轭梯度算法更新变反演网格中每个网格节点的S波(或和P波)速度使目标函数收敛,从而获得地下S波速度结构。在地震波波形反演中地震子波效应的处理十分重要。本文提出两种方式来处理这个问题。第一种方式是利用虚实源对比法(Virtual Real Source)从两炮(或多炮)具有相同检波器排列的炮集记录中提取地震子波。第二种方式是对实测数据进行反卷积处理消除数据中的地震子波效应,从而实现不依赖地震子波的波形反演。本文利用一些数值模拟算例和野外实测数据证明了这两种方法的有效性。
[Abstract]:......
【学位授予单位】：中国地质大学
【学位级别】：博士
【学位授予年份】：2016
【分类号】：P631.4
【相似文献】
相关期刊论文 前10条
1 蔚贵宏;瞬态面波在西藏水库勘察中的应用[J];工程勘察;2000年02期
2 王在民,李莉,朱鹏飞,张仲祜,郭宏宪,蒋立;用平均能量法分频识别和压制面波[J];新疆石油地质;2005年06期
3 陈爱萍,邹文;基于S-变换的面波压制技术[J];世界地质;2005年01期
4 张华;潘冬明;刘松;;小波包变换在面波分离中的应用[J];物探与化探;2007年02期
5 刘魁;程先琼;朱介寿;孔春玉;;面波有限频率法3-D灵敏度核函数的计算[J];成都理工大学学报(自然科学版);2007年03期
6 ;面波、反射处理软件综合技术服务[J];水利技术监督;2007年03期
7 毛海波;王在民;朱鹏飞;姚华;;平均能量-迭代法分频压制面波研究[J];天然气工业;2007年S1期
8 孙学文;陈文超;高静怀;包乾宗;;小波域强面波衰减方法研究[J];石油地球物理勘探;2008年01期
9 徐鑫;张学强;徐涛;张晓敏;;小波变换在压制面波中的应用[J];工程地球物理学报;2008年02期
10 李文杰;魏修成;宁俊瑞;张建伟;;基于频率衰减特性的面波压制方法[J];石油物探;2008年03期
相关会议论文 前10条
1 谢金娥;肖鹏飞;;数学形态学在压制面波中的应用[A];中国地球物理·2009[C];2009年
2 董烈乾;李振春;;第二代曲波变换压制面波方法研究[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
3 王伟涛;王宝善;;噪声互相关函数中面波前驱信号的来源分析[A];中国地球物理学会第二十七届年会论文集[C];2011年
4 靳洪晓;赵永贵;王超凡;;单道瞬态面波的硬件及软件实现[A];中国地球物理第二十一届年会论文集[C];2005年
5 刘魁;朱介寿;;面波有限频率法3-D灵敏度核函数的计算[A];中国地球物理学会第22届年会论文集[C];2006年
6 王德营;李振春;;二维S变换域自适应面波压制[A];中国地球物理2013——第十八专题论文集[C];2013年
7 张在金;徐辉;于海铖;徐苗苗;张军华;;基于图像特征转换的面波去噪方法研究[A];2014年中国地球科学联合学术年会——专题14：地下介质结构及其变化的地震面波、背景噪声及尾波研究论文集[C];2014年
8 高磊;陈文超;高静怀;;基于盲源分离的面波压制方法[A];中国地球物理·2009[C];2009年
9 鲁来玉;曾融生;丁志峰;;三维各向同性结构体对面波的多次散射[A];中国地球物理学会第22届年会论文集[C];2006年
10 熊章强;张大洲;宁刚;刘彦华;;瑞雷面波有限差分法数值模拟[A];中国地球物理学会第22届年会论文集[C];2006年
相关重要报纸文章 前2条
1 本刊记者 陈欣 谷玉强;10″，可贵的攻关[N];秦皇岛日报;2010年
2 ;我国首创电梯地震监测系统[N];政府采购信息报;2010年
相关博士学位论文 前10条
1 吴勃;面波和地球自由振荡的本征模理论[D];中国科学技术大学;2016年
2 郑晨;利用接收函数与面波频散研究青藏高原东南缘壳幔结构[D];中国地震局地球物理研究所;2016年
3 高玲利;高频面波中不同模态的特性及应用[D];中国地质大学;2016年
4 潘雨迪;二维弹性非均匀介质高频面波反演[D];中国地质大学;2016年
5 于彦彦;三维沉积盆地地震效应研究[D];中国地震局工程力学研究所;2016年
6 路拓;矿井巷道面波—转换横波超前探测技术研究[D];中国矿业大学;2016年
7 张智;川滇地区地壳及上地幔结构面波层析成像[D];吉林大学;2007年
8 李亚峻;地震勘探子波估计、面波消减新技术研究[D];吉林大学;2007年
9 张煜;部分饱和多相介质中Rayleigh面波传播特性研究[D];中国地质大学;2012年
10 房立华;华北地区瑞利面波噪声层析成像研究[D];中国地震局地球物理研究所;2009年
相关硕士学位论文 前10条
1 吴春梅;噪声源瑞利面波波场的物理模拟研究[D];西南交通大学;2015年
2 耿亮;基于分数阶傅里叶变换的地震面波压制方法研究[D];吉林大学;2010年
3 桂丽华;基追踪在面波分离与压制中的应用[D];西南交通大学;2010年
4 刘魁;面波有限频率层析成像理论及应用[D];成都理工大学;2007年
5 李媛媛;小波变换去除面波的方法研究[D];长安大学;2004年
6 梁丹丹;基于分频处理的叠前地震资料中面波压制[D];西南交通大学;2011年
7 张少梅;地震资料面波正演模拟及去噪技术研究[D];中国石油大学;2011年
8 李俊邦;基于多子波地震道分解技术消除面波[D];中国石油大学;2011年
9 蒋科植;利用面波方法研究欧亚大陆及周边海域岩石圈结构[D];成都理工大学;2012年
10 徐宏斌;欧亚大陆及西太平洋边缘海瑞利面波层析成像[D];成都理工大学;2006年
，
本文编号：
2456345
