收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

华北地区瑞利面波噪声层析成像研究

房立华  
【摘要】: 华北具有我国大陆最古老的克拉通构造,自中生代以来,华北克拉通活化,发生大规模的伸展和裂解,岩石圈减薄,并成为中新生代构造强烈活动的地区。华北地区不但是我国的政治、经济和文化中心,也是地震灾害多发区,历史上曾发生多次破坏性地震。近代发生的唐山、滦县、宁河强震造成了巨大的经济损失和人员伤亡,因此该区一直是我国的地震重点监视防御区。区内断层分布纵横交错,地质构造非常复杂。有鉴于此,许多学者利用多种地球物理方法对华北地区的地壳上地幔速度结构进行了研究,取得了许多重要成果,使人们对华北地区的地质构造背景,地壳上地幔深部构造环境,地震活动性及破坏性地震的成因和机理有了一定程度的了解和认识。 尽管前人就华北地区的地质构造和地壳上地幔速度结构进行了很多研究,但由于地震定位精度、台站分布的不均匀性、台站密度和数据质量等多方面因素的限制,以前的研究也存在一些不足。如:P波速度结构研究较多,而S波速度结构的研究相对较少;早期面波层析成像研究多采用1°×1°以上的网格划分,只能揭示较大尺度的构造框架和主要构造单元的基本性质,分辨率较低;射线平均路径长,短周期面波信号不足,对地壳浅部速度结构的约束较差。 为研究华北克拉通的破坏和岩石圈减薄,深入探讨华北地区强震的孕震背景,获取高分辨率的地壳上地幔三维速度结构,为该区域中长期破坏性地震预测和华北地区的防震减灾提供科学依据,中国地震局地球物理研究所于2006年11月份在华北地区开展了大规模流动地震台阵观测,在该区域布设了250套流动地震观测设备。作为华北地区宽频带流动地震台阵综合研究的一部分,本文对噪声层析成像的理论和方法进行了研究,并将该方法成功应用于华北地区。利用布设在华北地区的190个宽频带、10个甚宽带台的2007年1月份至2008年2月份垂直分量的连续波形数据,采用互相关方法提取出了瑞利面波的格林函数,利用多重滤波和相位匹配滤波方法提取了基阶瑞利面波的群速度频散曲线,得到5630条高质量的频散曲线,继而对研究区进行0.25°×0.25°网格剖分求得了每个网格的群速度纯路径频散,分辨率比以往的研究大有提高。利用遗传算法反演得到了每个网格下方0-50km深度范围的横波速度分布,构建了华北地区的三维S波速度结构。迄今为止,这是中国大陆分辨率最高的群速度层析成像结果。 主要研究内容及结果概述如下: 一、对从背景噪声提取格林函数的理论,频散曲线的测量和反演,以及面波层析成像的理论和方法进行了研究。开发了噪声层析成像的数据处理程序和群速度频散曲线的自动测量程序。测量由地震背景噪声提取出的瑞利面波的频散曲线时,使用多重滤波方法大多能获得较为满意的结果,对于信噪比较差的格林函数,可配合使用相位匹配滤波方法。 二、对提取出的格林函数分析后发现,从地震背景噪声中提取的频散曲线的最大周期与台阵孔径有关,台阵孔径越大,提取的频散曲线频带越宽;通过引入最大可信周期可使测得的频散数据更可靠;发现格林函数的信噪比与观测时间的平方根成正比,使用对称分量可以提高信噪比;噪声源的不均匀分布,导致格林函数的不对称性;使用一年以上的观测数据,可以降低噪声源分布不均匀性对格林函数的影响,使提取出的格林函数更加对称,信噪比更高。 三、对噪声特征进行了研究,发现不同周期段的噪声特征是不一样的。4-10s周期段,互相关函数在零时刻附近出现一组振幅强、震相连贯的波组。由于短周期面波强烈的衰减和吸收效应,提取出的这一周期段的瑞利面波格林函数的信噪比较低。在10-20s周期段,噪声源存在季节性变化,噪声强度与全球海浪高度较一致。在20-50s周期段,瑞利面波格林函数基本对称且不随季节性变化,表明这一周期段的噪声源分布较均匀。在4-20s周期段,尽管互相关函数的正负分支在振幅上存在很明显的不对称性,但是到时却基本一致,说明噪声源分布对互相关函数振幅的影响较大,对到时的影响较小。 四、将研究区域划分为0.25°×0.25°的网格,用面波层析成像方法,反演得到了4s-40s周期的群速度分布图像。层析成像结果较好地揭示了地壳内部、尤其是浅部地壳的横向速度变化。短周期的群速度分布同地表地质结构、地形密切相关;华北盆地及山间的沉积盆地,如延庆-怀来盆地,大同盆地,显示低速异常,而基岩广泛出露的太行山和燕山隆起区,呈现高速异常;我们计算了9s-12s的群速度分布图的速度梯度场,将速度梯度划分为三个区间:<0.05km/s,0.05-0.15km/s,>0.15km/s,并同时给出了公元前780年以来发生的M≥6.0级的强震的震中分布,发现约2/3的强震(M≥6.0级)都发生在中等速度梯度带(0.05-0.15km/s)附近。 五、使用光滑约束的遗传算法反演了432个节点的瑞利面波群速度纯路径频散曲线,得到了华北地区0-50km深度的三维S波速度结构。反演结果较好地揭示了华北地区地壳内部的S波速度结构、结晶基底埋深和莫霍面埋深的横向变化。反演结果表明太行山、燕山隆起区结晶基底埋深一般约<2km,延庆-怀来盆地约5km,大同盆地约3km。华北盆地内部结晶基底厚度一般>6km,且不同构造单元存在差异,沧县隆起约8km,冀中坳陷和黄骅坳陷约11km。研究区域的地壳厚度横向变化较大,基本形态是西部厚、东部薄,由东南向西北逐渐加深,莫霍面的区域构造方向呈北东或北北东向展布。地壳厚度从东部渤海湾附近的28km逐渐向西北方向变厚,进入山西高原、张家口地区达到44km,其横向变化达16km。在太行山山前断裂附近,莫霍面厚度由东部的36km左右,至西部陡变为40km左右。燕山隆起带地壳厚度由渤海边的32km,向西逐渐加深至约44km。 六、S波速度横向分布图揭示:在0-12km深度范围,太行山山前断裂以东表现为低速异常,而以西表现为高速异常,13km以下,东西两部分的高低速速度异常开始变得模糊,仅在部分区域较明显。据此推测,太行山山前断裂可能仅在部分区域延伸至莫霍面。在0-8km深处,燕山隆起带为大范围的高速异常;在燕山隆起带南侧、与张渤地震带相邻的地方,高低速异常交替分布,可能是受北北东-北东向断裂的切割所致;沿张渤地震带附近,在10km深处有一条较为连续的北西向低速异常,这一低速异常带和燕山高速区的南缘勾画出了张渤断裂带的位置和北部边界。北京-天津-唐山区域在中下地壳(15km-32km)存在显著的低速异常,可能是由地幔热物质侵入造成的。 七、对唐山、滦县、宁河强震区的S波速度结构分析发现,唐山地震和滦县地震发生在高速体内部,且震源下方存在低速体,宁河地震发生在S波低速体与高速体的过渡带。震源区上地幔顶部的S波速度较低,莫霍面在宁河西南和滦县存在局部隆起。根据本文反演得到的高分辨率的S波速度结构,我们认为这三个强震的成因与上地幔的垂直形变、壳幔物质的交换和地壳中的流体有关。幔源岩浆沿壳幔边界的断错侵入,导致上地幔顶部的速度较低。同时岩浆的侵入使下地壳增热,粘度下降,某些矿物发生脱水作用,脱出的水上移并存储于中地壳内。这些流体长时间的存在影响了断层带的结构和组成,进而改变了断层带的应力状态,造成了孕震区的弱化,从而触发地震。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 王在民,李莉,朱鹏飞,张仲祜,郭宏宪,蒋立;用平均能量法分频识别和压制面波[J];新疆石油地质;2005年06期
2 陈爱萍,邹文;基于S-变换的面波压制技术[J];世界地质;2005年01期
3 毛海波;王在民;朱鹏飞;姚华;;平均能量-迭代法分频压制面波研究[J];天然气工业;2007年S1期
4 李杨;董守华;杜庆顺;;基于广义S变换的面波压制技术研究[J];能源技术与管理;2009年03期
5 鲁来玉;丁志峰;何正勤;;浅层有限频率面波成像中的3D灵敏度核分析[J];地球物理学报;2011年01期
6 申振华;于光明;杨振;颜中辉;王赟;;面波压制方法与道间距关系的数值分析[J];煤田地质与勘探;2012年01期
7 王振东;;微动应用技术讲座 第二讲:由微动数据提取面波的方法[J];国外地质勘探技术;1990年05期
8 陈定国,李衍达;一种新的消除面波的方法[J];石油物探;1992年04期
9 罗国安,杜世通;小波变换及信号重建在压制面波中的应用[J];石油地球物理勘探;1996年03期
10 蒋连斌,陈道宏,高振山,廉国芬,刘学伟,何樵登;一种去面波的新方法[J];石油地球物理勘探;1997年S2期
11 辛小春;面波新技术在工程地质勘察场地评价中的应用[J];北京水利;2003年01期
12 蔡宗文,庄进耀,张清香,郭铁栓;瑞雷面波在填海造地强夯检测工程中的应用[J];福建地震;2003年01期
13 蔡宗文,庄进耀,张清香,郭铁栓;瑞雷面波在填海造地强夯检测工程中的应用[J];华南地震;2003年04期
14 Laurent Duval;唐祥功;;丘陵地区面波衰减的两种革新方法[J];油气地球物理;2006年02期
15 梁志强;邓道静;;基于多阶模式的面波技术应用探讨[J];勘探地球物理进展;2008年04期
16 陈伟;徐义贤;张煜;;瑞雷面波在弹性楔状体中的传播特性研究[J];工程地球物理学报;2009年02期
17 陈文超;高静怀;包乾宗;;基于连续小波变换的自适应面波压制方法[J];地球物理学报;2009年11期
18 马见青;李庆春;;S变换和TT变换联合压制地震面波[J];物探化探计算技术;2011年03期
19 段俊;张白林;潘树林;;广义S变换去除面波技术及应用[J];物探化探计算技术;2011年03期
20 石颖;裴江云;林春华;何恕;赵旭;刘宇;;十字交叉排列面波压制方法及应用[J];东北石油大学学报;2012年05期
中国重要会议论文全文数据库 前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年
中国博士学位论文全文数据库 前6条
1 张智;川滇地区地壳及上地幔结构面波层析成像[D];吉林大学;2007年
2 李亚峻;地震勘探子波估计、面波消减新技术研究[D];吉林大学;2007年
3 张煜;部分饱和多相介质中Rayleigh面波传播特性研究[D];中国地质大学;2012年
4 房立华;华北地区瑞利面波噪声层析成像研究[D];中国地震局地球物理研究所;2009年
5 王未来;中国西南地区的接收函数和面波层析成像研究[D];中国地震局地球物理研究所;2014年
6 苏伟;青藏高原及邻区壳幔速度结构及面波方位各向异性[D];中国地震局地球物理研究所;2012年
中国硕士学位论文全文数据库 前10条
1 耿亮;基于分数阶傅里叶变换的地震面波压制方法研究[D];吉林大学;2010年
2 桂丽华;基追踪在面波分离与压制中的应用[D];西南交通大学;2010年
3 刘魁;面波有限频率层析成像理论及应用[D];成都理工大学;2007年
4 李媛媛;小波变换去除面波的方法研究[D];长安大学;2004年
5 梁丹丹;基于分频处理的叠前地震资料中面波压制[D];西南交通大学;2011年
6 张少梅;地震资料面波正演模拟及去噪技术研究[D];中国石油大学;2011年
7 李俊邦;基于多子波地震道分解技术消除面波[D];中国石油大学;2011年
8 蒋科植;利用面波方法研究欧亚大陆及周边海域岩石圈结构[D];成都理工大学;2012年
9 徐宏斌;欧亚大陆及西太平洋边缘海瑞利面波层析成像[D];成都理工大学;2006年
10 刘朝;基于Shearlet变换的面波压制[D];哈尔滨工业大学;2013年
中国重要报纸全文数据库 前2条
1 本刊记者 陈欣 谷玉强;10″,可贵的攻关[N];秦皇岛日报;2010年
2 ;我国首创电梯地震监测系统[N];政府采购信息报;2010年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978