龙门山及其邻区航磁、卫星磁测资料处理与解释.docVIP

龙门山及其邻区航磁、卫星磁测资料处理与解释 赵亚博 刘林静 王赛昕 任宇鹏 （中国地质大学 地球物理与空间信息学院 武汉 430074） 摘要：龙门山断裂带地处青藏高原东缘，大地构造位置特殊，断裂活动强烈.本文根据航磁和卫星磁测资料对龙门山及其邻近地区（102°E～107°E,29.8°N～33.8°N）的磁异常特征、断裂体系进行了分析并对居里等温面进行了反演与解释,揭示航磁、卫磁异常与断裂活动、深部物质流动的关系，为研究该区域复杂的地质构造背景，揭示研究区地震高发的深层原因提供帮助. 关键词：龙门山断裂;航磁;卫星磁测;小波多尺度分析;居里等温面 1 龙门山构造带的构造特征 图1 龙门山及其邻区断裂带分布及历史地震分布图（来源于国家地震科学数据共享中心） 龙门山构造带位于扬子地块的西缘,其东、西侧分别为松潘至甘孜褶皱带和四川盆地,北东向止于秦岭造山带,总体呈北东-南西向延伸,是中国最典型的逆冲推覆构造之一[1].大量地质学研究表明,龙门山断裂带主要由3条断裂组成,自东向西分别为龙门前山主边界(包括大川—双石断裂、灌县—二王庙—安县断裂和江油(或称马角坝)断裂)断裂;龙门山主中央(包括盐井—五龙断裂、北川—映秀断裂、茶坝—林庵寺断裂)断裂;后龙门山(汶川—茂汶断裂与青川—平武断裂)断裂.由于其所处的大地构造位置特殊,长期以来被地质学家所重视并在该区进行了大量的研究工作. 2 小波多尺度分解 小波变换是磁异常分解的有效工具,利用小波多尺度分析方法,可以将磁异常分解到不同的空间尺度上,尺度的大小决定了磁异常所反映的地质体规模和埋深的大小.根据小波多尺度分析方法的原理[2],地磁异常可以分解成不同空间尺度下逼近和细节异常,假设二维卫星磁测异常场为△T(x,y)= f(x,y),进行3阶小波分解后卫星磁测异常分解表达式可简记为： △T=A3G+D1G+D2G+D3G 式中,D1G～D3G为1～3阶小波变换异常细节,A3G为3阶小波变换异常逼近.通过对1～3阶小波变换细节进行径向对数功率谱分析,可以求取其对应的场源深度[3]. 采用以上方法首先对研究区航磁异常图进行小波多尺度分解，并对分解后的各阶磁异常细节进行了功率谱分析，最后对该地区的断裂及构造进行不同深度分析.同样可以利用小波变换来分析信号的奇异性,在信号突变的位置上,小波变换后的系数具有模量极大值,可用于断裂分析,来判断如断裂等构造形迹的走向及位置等. 3 研究区航磁异常特征分析 3.1航磁异常图的分区特征 随着航空磁测技术的发展，利用航空磁测对各种地质问题进行分析解释越来越受地球物理学家的重视，由于航空磁测接近磁异常的源,因而分辨率较高,对磁性体的几何形态可作出比较精确的计算和描画.本区航磁异常在磁场强度、分布范围、形态、梯度和走向等方面有较为显著的规律性,形态分布规律性较强,代表不同的地质构造单元,经分析可将研究区磁异常划分为以下三个区域（图2）： = 1 \* ROMAN I 川西坳陷高值异常区 图2 龙门山及其邻区磁异常图（单位：nT） 川西凹陷高值异常区为泸定—都江堰—绵阳—巴中一线，自西南-北东向背景磁异常升高，分布有两个近椭圆形的局部高值正异常,以南充为中心的异常最为明显，最大值490nT.推测该区升高的磁场背景是由元古代磁性基底引起[4]. = 2 \* ROMAN II 龙门山磁异常陡变带 龙门山磁异常陡变带在研究区呈北东—南西方向展布,斜穿整个研究区,该异常带是整个研究区内等值线较密集、磁场最为复杂、局部异常特别发育的地带.按地质构造划分,这一北东向异常带在龙门山冲断带的范畴之内,与本地区存在的三条主干断裂对应.航磁异常图上显示,在理县的西南侧显示一近圆形高值异常区,与地表出现的中生代花岗岩体群对应.在江油附近出现的北东向正磁异常可能与地表出现很厚的飞仙关组对应.在汶川到灌县有+150nT的局部磁异常，与地表的茂汶杂岩体对应[5].该带中部茂县—平武地区显示北东向负异常带,峰值为-105nT,推测这一异常反映的是志留系茂县群的一套变质程度不同的砂泥岩夹碳酸岩.该区东北部,在负异常背景上,出现了椭圆状的局部升高,长轴方向为东西方向,峰值-10nT,推测该异常与该处上三叠统须家河组碎屑砂岩岩屑中所含的大量花岗岩岩屑和火成岩岩屑相应. 研究区地震沿龙门山构造带南段密集成带分布,形成一条北东向的小地震密集活动条带,而龙门山构造带北段地震活动相对稀疏。 = 3 \* ROMAN III 红原—黑水低缓负异常区 红原—黑水低缓负异常区位于研究区西北部,为极平缓的低值负异常区,负磁场背景平稳缓变,负磁场峰值为-90nT,该区展布的负磁场背景是由台区上三叠统基底引起的[4]. 3.2小波分析结果分析 地名？其他北东向异常？ 图3 小波分析2阶细节图