中国青藏高原深部地球物理探测与地球动力学研究(1958-2004)

经过半个多世纪的努力,中国学者在青藏高原进行的地壳上地幔地球物理探测工作取得了巨大进展. 累计完成长约45000km的深部地球物理探测工作,取得了许多科学数据,为探讨高原地壳上地幔结构、隆升机制和动力学研究奠定了基础. 为比较全面反映中国学者多年来的工作成果,作者广泛收集资料,总结了中国学者在青藏高原地壳上地幔地球物理探测工作程度,并按照方法分类绘制了系列工作程度图. 文中分别对地壳结构、上地幔的横向不均匀性、岩石圈的电性结构、青藏高原隆升机制、青藏高原地球动力学模型等几个方面所取得的主要成果做了概略的评述. 已有资料表明：青藏高原的莫霍面埋深变化较大,且在几条重要缝合带莫霍面两侧都有断错;根据目前的探测结果,高原在20±5km埋深范围内普遍存在壳内低速高导层,速度一般为56～58km/s,电阻率约为1～10Ω·m,厚度一般为5～10km,但横向分布不连续. 低速层与高导层的深度、厚度在趋势上一致,但不十分吻合.天然地震的研究结果表明,组成高原各个地块内部的地震各向异性方向大致相同,各地块的分界处各向异性方向往往有明显的变化;虽然对高原隆升机制还存在不同的看法,但至少认为高原是多期隆升、多种机制共同作用的结果这一点已达成共识. 综合已有的地球物理调查成果,结合地质地球化学资料建立的高原地球动力学模型,形象地表达出青藏高原岩石圈的双向挤压变形模式. 这些工作为研究青藏高原隆升和变形机制提供有价值的信息.     

ELECTRICAL STRUCTURE OF UPPER CRUST IN THE SOURCE REGION OF JINGGU YUNNAN MS6.6 EARTHQUAKE AND THE SEISMOGENIC ENVIRONMENT

The October 7, 2014 M_S6.6 earthquake in southwest of Jinggu in the southwestern Yunnan Province occurred as the result of shallow strike-slip faulting within the crust of the Eurasia plate in the broad plate boundary region between the India and Eurasia plates. The strike of fault plane is 140°, and the aftershock distribution shows that the rupture plane is also NNW-trending. Tectonics of the region are controlled by the convergence of the India plate with Eurasia, which has driven the uplift of the Himalayas to the west of this earthquake, and has caused the formation of numerous intraplate continental transform structures in the surrounding region. The pattern of elastic-wave radiation from the earthquake is consistent with the shock occurring either as the result of right-lateral faulting on a northwest-trending fault or as the result of left-lateral faulting on a northeast trending fault. Faults of both types have been mapped in southwestern Yunnan, and it is unclear at this time which type of fault hosted this event. Magnetotelluric survey line is across Jinggu earthquake zone. The advanced data processing and analysis technology of MT is employed and the quantitative data from field surveys are analyzed to acquire the reliable electrical model. The MT data are inverted using nonlinear conjugate gradient (NLCG) inversion algorithm. At last, the interpretation of the electrical model is performed considering the geology and the other geophysical data. Based on the final inversion model of the target profile, it is found that:(1) Electrical structure of the source region can be divided into four layers:The surface is relatively low resistivity layer(0~5km), consisting mainly of Mesozoic and Cenozoic Basin sedimentary rocks, the value of resistivity is 100Ω·m; The high resistivity layer(5~10km) in upper crust mainly consists of Proterozoic metamorphic rocks, with resistivity higher than 1 000Ω·m; there are the upper crust high-conductivity layer(15~25km) and crust-mantle transition zone(blow 25km); (2) The focal depth of the Jinggu earthquake is about 10km, which locates in the interface between high resistivity layer and high-conductivity layer; (3) Most of the focal depths of the aftershocks are in the range of 5km and 10km, and the two depths(5km & 10km) are corresponding to the resistivity gradient belt.     

滇西三江构造带电性结构特征——以福贡—巧家剖面为例

为查明滇西三江构造带及邻区复杂的构造特征,并揭示该区深部电性结构,沿福贡—巧家布设了一条长约410km的大地电磁剖面.共观测到61个物理点,其中宽频大地电磁测点41个,长周期大地电磁测点20个.通过对采集到的数据进行一系列的处理、反演,得到了沿剖面的壳幔电性结构模型.并结合研究区内区域地质资料及其他地球物理资料,对剖面所经过的各个主要地质构造单元及主要断裂带进行了综合解释.电性结构模型揭示沿剖面地壳电性层次复杂,深部电性结构由西往东呈分块展布,横向变化大,壳内广泛发育低阻异常.在中甸构造带(香格里拉地块)和盐源—永胜构造带深部壳幔存在大规模低阻异常,这可能与地下局部熔融体和地热流有关;康滇构造带壳幔存在大规模高阻异常,表明地壳中曾经有地幔物质侵入;在大凉山构造带地下10～50km深处存在一呈横向"半月形"展布的低阻体,电阻率值不满10Ωm,结合地质资料与前人的研究成果,推测该低阻体成因应与青藏高原东南缘"地壳管道流"有一定关联.     

Electrical Struture of the Crust in Manas Earthquake Area,Xinjiang,China

Magnetotelluric sounding data obtalned recently in Manas earthquake area were processed. Inthe result, curves of apparent resistivity, impedance Phase, skewness and optimum rotationangle versus period and the real magnetic induction vectors were obtained. Then the data ofall measuring points were interpreted by 2D automatic inversion. The result indicates thatalong the sounding profile the shallow crust can be divided into 5 segments and the deep crustcan be divided into 3 segments, with faults or deep-seated fault zones as the contactboundaries between them. The sedimentary cover along the profile extents down to depthabout 12 km in maximum and a low-resistivity body exists in the crust in southern section ofthe profile. The interpretation results are well consistent with geological and othergeophysical data. The Manas M7. 7 earthquake occurred near a contact zone where theelectrical structure of the crust sharply changes.     

利用接收函数和大地电磁数据联合反演南迦巴瓦构造结中部地区壳幔结构

利用2010年布设在西藏南迦巴瓦构造结的郎嘎、崩嘎、直白和拉格四个宽频地震台所观测到的近5个月的地震记录,采用时间域迭代反褶积技术处理得到接收函数,通过筛选多条相近震中距和反方位角的高质量接收函数求取其叠加平均.对大地电磁数据做Rhoplus分析处理得到视电阻率和相位曲线.利用单台接收函数和相同位置的大地电磁视电阻率和相位联合反演地下一维壳幔结构.联合反演采用遗传算法,并通过权衡图分析大地电磁和地震数据的兼容性.理论值和实测值的对比显示两种数据能同时得到较好拟合.联合反演结果表明:(1)中上地壳为9 km至14 km厚的高阻高速层覆盖于低阻低速层之上的结构,中地壳低阻低速层可能与深部流体和局部熔融共同作用有关.(2)下地壳存在最厚达20 km的高导的壳幔过渡层,波速在4 km/s左右;上地幔约130 km至150 km以下存在软流圈.(3)上地壳的高阻高速层解释为多雄拉组混合岩化角闪岩相变质岩,而直白台所显示的低阻低速层与高压麻粒岩的少量部分熔融有关,可能源于壳幔过渡带镁铁质岩石的相变或更深处幔源岩浆底侵作用的产物.     

Magnetotelluric Studies of Active Continent–Continent Collisions

Continent–continent collisions are an important tectonic process and have played a fundamental role in the evolution of the modern continents. A combination of geological and geophysical data has provided new constraints on the structure and temporal evolution of these orogens. Magnetotelluric (MT) studies have been an important part of these studies since they can constrain the fluid content and thermal structure which are key parameters for defining the rheology of the crust and upper mantle. MT studies of the Himalaya have defined the geometry of active faults associated with continued plateau growth. Orogen scale MT studies have shown that both the India–Asia collision (Tibetan Plateau and Himalaya) and the Arabia–Eurasia collision (Eastern Anatolia) have developed a low resistivity mid-crustal layer with upper surface at 10–20 km that is likely due to a combination of partial melt and associated aqueous fluids. The properties of this layer are consistent with a strength contrast that permits crustal flow over geological timescales. The upper mantle from the Moho to at least 100 km beneath both Northern Tibet and the Anatolian Plateau is characterized by low resistivity values (10–30 Ωm) indicating the presence of shallow asthenosphere. Future integrated seismic and MT studies of collision zones are needed fully to explore the 3D structures associated with deformation and further constrain geodynamic models.     

长白山天池火山区及邻近地区壳幔结构探测研究

对长白-敦化深地震测深剖面资料利用二维射线追踪程序包进行走时拟合及地震图计算,得到了长白山天池火山区及邻近地区地壳上地幔速度结构和深部构造. 结果表明,以C2界面为标志,研究区地壳可分为上部地壳和下部地壳. 上部地壳厚1-23km,P波速度为6.00-6.25km/s;下部地壳厚12-17km,它是由一个较均匀的速度层和一个厚6-km的壳幔过渡层构成. 地壳厚度由敦化一带31-33km向东南逐渐增厚,至天池火山区最深达3km. 在天池火山区地壳存在低速体,其速度较周围介质低约为0.15km/s. 利用地震剖面探测、地震CT和大地电磁测深等结果显示,在天池火山区地壳内存在低速、低密度及低阻异常体,该异常体可能表明壳内岩浆囊的存在.     

江苏及邻区地震重新定位和构造特征分析

本文将双差定位和遗传算法定位两种方法相结合,对江苏及邻区1980～2005年的地震进行重新定位研究.重新定位后大大改善了原地震定位的精度,地震在空间上更加集中分布在某些区域,地震震源深度分布结果更加合理.江苏及邻区地震主要发生在上地壳和中地壳;研究区不同构造单元的震源深度统计分析表明,下扬子断块江苏段、冀鲁断块和豫皖断块的地震震源深度特点相似,在10～11 km、15～17 km深处存在二个明显的地震优势分布,推测分别在上地壳底面和中地壳;在25 km深处也存在一小的地震优势分布面,但地震频次较低.大别山地区地震的震源深度与下扬子断块、冀鲁断块和豫皖断块内的地震震源深度存在明显的差异,主要差异为大别山地区10 km以上的浅源地震十分发育,在6～7 km的深处有一地震优势分布,该深度附近地震波速度可能较高,而在10 km以下差异不大.     

Deep massive sulphide exploration using 2D and 3D geoelectrical and induced polarization data in Skellefte mining district,northern Sweden

Geoelectrical and induced polarization data from measurements along three profiles and from one 3D survey are acquired and processed in the central Skellefte District, northern Sweden. The data were collected during two field campaigns in 2009 and 2010 in order to delineate the structures related to volcanogenic massive sulphide deposits and to model lithological contacts down to a maximum depth of 1.5 km. The 2009 data were inverted previously, and their joint interpretation with potential field data indicated several anomalous zones. The 2010 data not only provide additional information from greater depths compared with the 2009 data but also cover a larger surface area. Several high‐chargeability low‐resistivity zones, interpreted as possible massive sulphide mineralization and associated hydrothermal alteration, are revealed. The 3D survey data provide a detailed high‐resolution image of the top ～450 m of the upper crust around the Maurliden East, North, and Central deposits. Several anomalies are interpreted as new potential prospects in the Maurliden area, which are mainly concentrated in the central conductive zone. In addition, the contact relationship between the major geological units, e.g., the contact between the Skellefte Group and the Jörn Intrusive Complex, is better understood with the help of 2010 deep‐resistivity/chargeability data. The bottommost part of the Vargfors basin is imaged using the 2010 geoelectrical and induced polarization data down to ～1‐km depth.     

天山中段的深浅构造特征

天山起源于古生代的陆-陆碰撞作用,又经历了中新生代的典型陆内造山过程,其深浅构造结构和活动性一直是众多学者关注的热点。文中通过多种地球物理探测和综合地质构造分析, 以地学断面形式对其深浅构造进行填图,揭示了天山中段复杂的深浅构造特色。结果表明:沿古生代的陆-陆碰撞缝合带两侧分别呈现出主要构造地质单元由老到新的对称性,并伴有相应深部结构的复杂性,反映了碰撞过程及后期的构造演化特点;天山中部的上地幔顶部存在厚近 10km、宽近 200km、几乎涵盖整个天山的低速高导层,可能是中新生代以来天山的陆内再造山作用引起的壳幔拆沉作用形成的残留下地壳     

我国西北部分地区地壳-上地幔电性横向变化特征与地震活动的关系

由50余个测点的大地电磁测深资料,讨论了该区的地壳-上地幔电性横向变化特征。按照上地幔第一低阻层顶面埋深,将测区划分为四类地区:浅埋深(55-90公里)、中浅埋深(90-110公里)、中深埋深(110-160公里)和深埋深(160-250公里)。讨论了本区六次大震例的深部电性背景。指出了上地幔顶部的梯度带地区、地壳内电性横向变化剧烈的地区和地壳内存在局部特殊增厚的低阻层地区将可以作为潜在震源区的深部电性判据之一     

唐山地震区地壳结构和构造:深地震反射剖面结果

1976年7月28日,在唐山地区发生了7.8级大地震.为了研究该区的地壳结构和断裂的深浅构造关系,2009年,我们在唐山市南部的丰南地区,跨唐山断裂带完成了1条道间距40m、炮间距200m、50次覆盖的深地震反射探测剖面.结果表明:研究区的地壳厚度为32 ～ 34km,莫霍面自东向西逐渐加深,在丰南县和宣庄镇之间,中-...     

云南遮放-宾川和孟连-马龙宽角地震剖面的层析成像研究

依据有限差分反演和射线反演的方法,利用走时、振幅比和重力布格异常数据对云南地区遮放—宾川和孟连—马龙宽角地震剖面的地壳上地幔结构进行层析成像研究. 遮放—宾川测线地壳厚度为35～46km,孟连—马龙测线地壳厚度33～44.5km. 局部位置上地幔顶部P波速度值偏低,速度值变化范围大,反映出云南地区是典型的构造活动区的特点. 地壳上地幔结构存在一定程度的深浅构造一致性,意味着浅部物质活动存在深部背景. 遮放—宾川剖面速度结构显示地表怒江断裂东侧存在贯穿地壳的大规模低速异常,可能与深部物质上涌活动有关. 遮放—宾川线和孟连—马龙线速度结构显示作为一级构造单元分界线的红河断裂是超壳断裂,怒江断裂深及上地幔. 而昌宁—双江断裂显示为低角度铲式断层,意味着该断裂切割并不深. 云南地区强震的发生往往与延伸到上地幔的深大断裂有关,且一些浅源地震经常位于中上地壳深大断裂与其他断裂的交汇部位、高速块体与低速块体接触带上速度等值线弯曲的位置. 推测这样的位置有利于能量和区域应力的积累.     

用转换函数方法研究腾冲—临沧地区地壳结构

根据流动数字地震台网提供的三分量地震波形记录资料,应用转换函数及快速模拟退火算法对腾冲－临沧地区30个地震台站下的地壳横波速度结构进行了反演.反演结果说明,研究区壳幔边界清晰、莫霍界面附近速度跳跃明显,由此得出该区地壳厚度在40 km左右、并具有从南向北增厚趋势.一个普遍的现象是,在腾冲-宝山地块下地壳存在明显的低速带,低速带的厚度在10～20 km间.研究结果进一步表明各台站下方上地幔速度结构复杂.这些结果为探讨青藏高原东南缘下地壳的侧向黏性流动、碰撞板块边界处壳幔物质交换等均提供了重要的地球物理证据,为探讨印-藏汇聚过程中青藏高原东构造结岩石圈变形、高原隆升及其深部动力学有一定理论意义.     

北京南部地壳精细结构深地震反射探测研究

为了研究北京平原区的地壳结构特征、断裂的空间展布、断裂活动性以及深浅构造关系,在北京平原区的南部完成了1条长90 km的深地震反射剖面.探测结果表明,该区地壳以双程反射时间(TWT)6～7 s的强反射带Tc为界分为上地壳和下地壳,上地壳厚约18～19 km,下地壳厚约16～17 km,Moho界面深度约为34～35 km.该区结晶基底起伏变化较大,上、下地壳分界面和Moho界面都是一个具有一定厚度的过渡带.上地壳反射层位丰富,断裂构造发育,构造形态清晰.在夏垫断裂西北,剖面揭示了4～5组能量较强的反射震相,表现为典型的隆起区特征;在夏垫断裂东南,上部为一套向东南倾伏的密集强反射层,下部为一套形态各异、结构复杂的强反射层,这些反射具有典型的沉积盆地特征,盆地最深处约为11 km.剖面揭示的地壳深断裂倾角较陡,向上切穿了上、下地壳分界面,延伸到上地壳沉积盆地的底部,向下切穿了壳幔过渡带,与上部断裂和沉积盆地构成了独特的组合关系.     

内蒙古锡林浩特-东乌旗剖面壳幔电性结构研究

为研究二连-东乌旗贺根山一带成矿构造环境,提供矿产资源勘查、预测、评价的地质背景依据,跨贺根山和锡林浩特板块缝合带一线布设了26个超宽频带长周期大地电磁测深点,点距3～6 km,剖面长度100 km,在对获取的资料采用Robust变换、互参考处理的基础上,定性分析了视电阻率和相位曲线、二维偏离度、电性主轴,并采用二维共...     

青藏高原东北缘祁连山造山带至阿拉善地块壳幔电性结构研究

在青藏高原东北缘祁连山造山带至阿拉善地块之间完成了一条372km的大地电磁剖面,通过二维反演计算,获得了沿剖面180km深的壳幔电性结构模型,结合研究区地质和地球物理资料开展综合分析,研究结果表明:(1)剖面自南向北所经过的祁连山造山带、走廊过渡带和阿拉善地块对应3种壳幔电性结构模型:东祁连壳幔高-低-高阻似层状电性结构、河西走廊壳幔低阻带状电性结构和阿拉善南缘壳幔高-低-高阻层状电性结构.(2)剖面所经过的主要断裂带在电性结构上表现为低阻异常带或电性梯度带,并且止于中上地壳或消失于下地壳低阻层中.除这些分布于中上地壳的断裂系统以外,在下地壳至上地幔顶部还存在两条切割莫霍面的壳幔韧性剪切带:西华山北缘壳幔韧性剪切带和阿拉善南缘壳幔韧性剪切带.其中,西华山北缘壳幔韧性剪切带可能是1920年海原8.6级地震发生的深部背景之一;而阿拉善南缘壳幔剪切带可能是卫宁北山燕山晚期和喜山期幔源岩浆上升到地壳浅部或喷出到地表的通道,为在该区域寻找晚中生代至新生代含矿隐伏岩体提供了深部电性结构依据.(3)由若干形状不规则、彼此不相连的"碎块状"极高阻块体组成的中上地壳与"似层状"的中下地壳低阻层共同构成的地壳电性结构,是引起青藏高原东北缘强烈破坏性地震最佳的地壳电性结构组合之一.印度板块向欧亚板块俯冲碰撞楔入引起青藏高原块体向北东方向运移与阿拉善地块向南的俯冲碰撞楔入,是青藏高原东北缘强震活动带产生的动力学背景.     

汶川MS8.0地震：地壳上地幔S波速度结构的初步研究

2008年5月12日我国四川省汶川地区发生了震惊世界的M_S8.0地震.历史上,同类地震在大陆内部极为罕见.该地震深部构造背景的研究对理解其成因极为重要.本文利用中国地震局地质研究所地震动力学国家重点实验室在川西地区布设的大规模密集流动宽频带地震台阵记录的远震P波波形数据和接收函数非线性反演方法,得到了沿北纬31°线的19个台站下方120 km深度范围内的S波速度结构及台站下方地壳的平均泊松比.该观测剖面穿越了主震区,总长度约为420 km. 我们的结果揭示了川滇地块、松潘-甘孜地块和四川盆地三个不同地块构造差异.上述三个地块的地壳结构特征可以概括为：（1）四川盆地前陆壳幔界面向西侧倾斜并有较为明显的横向变形,地壳厚度存在46～52 km的横向变化,中下地壳S波速度存在横向变化,地壳平均泊松比值较高（0.28～0.31）,但在龙门山断裂带附近,显示了坚硬地壳的特征,地壳平均泊松比仅为0.2;（2）松潘-甘孜地块地壳厚度由西侧靠近鲜水河断裂的60 km,向东减薄为52 km,在14～50 km深度范围内存在S波速度2.75～3.15 km/s的楔状低速区,其厚度由西侧的～30 km向东逐渐减薄为～15 km,相应区域的地壳平均泊松比高达0.29～0.31; （3）鲜水河断裂西侧,川滇地块地壳结构相对简单,地壳厚度为58 km,并在26 km深度存在约10 km厚度的高速层,地壳内平均泊松比约为0.25;（4）汶川大震区在12～23 km深度上具有近乎4.0 km/s的S波高速结构,而其下方的地壳为低速结构,地壳平均泊松比0.31～0.32,汶川大震的余震序列主要分布在高速介质区域内. 本文的结果表明松潘-甘孜地块的地壳相对软弱;而且并不存在四川盆地向西侧的俯冲.我们认为在青藏高原东向挤压的长期作用下,四川盆地强硬地壳的阻挡作用可导致松潘-甘孜地块内部蓄积很大的应变能量以及上、下地壳在壳内低速层顶部边界的解耦,在龙门山断裂带附近形成上地壳的铲形逆冲推覆.汶川大地震及其邻近区域所具有的坚硬上地壳和四川盆地的阻挡作用为低应变率下的高强度应力积累创造了必要条件,而松潘-甘孜地块长期变形积累的高应变能构成了孕育汶川大地震的动力来源.     

3D upper crustal density structure of the Deccan Syneclise,Central India

A constrained 3D density model of the upper crust along a part of the Deccan Syneclise is carried out based on the complete Bouguer anomaly data. Spectral analysis of the complete Bouguer gravity anomaly map of the study region suggests two major sources: short wavelength anomalies (<100 km) caused primarily due to the density inhomogeneities at shallow crustal level and long wavelength anomalies (>100 km) produced due to the sources deeper than the upper crust. A residual map of the short wavelength anomalies is prepared from the complete Bouguer anomaly using Butterworth high‐pass filter (100 km cut‐off wavelength). Utilizing the constraints from deep resistivity sounding, magnetotellurics and deep seismic sounding studies, 2.5D density models have been generated along 39 profiles of this region. The mismatch between the calculated response of the a priori 2.5D model with the residual (short wavelength) gravity anomalies is minimized by introducing high‐density intrusive bodies (≥2.81 g/cm³) in the basement. With these 2.5D density models, the initial geometry of our 3D density model, which includes alluvium, Deccan trap, Mesozoic sediment and high‐density intrusive bodies in the basement up to a depth of 7 km (upper crust), is generated. In the final 3D model, Deccan trap extends from 200 m to nearly 1700 m below the 90–150 m thick Quaternary sediment. Further down, the sub‐trappean Mesozoic sediment is present at a depth range of 600–3000 m followed by the basement. The derived 3D density model also indicates six intrusive bodies of density 2.83 g/cm³ in the basement at an average depth of about 4–7 km that best fits the residual gravity anomaly of the study area.     

内蒙、甘、宁某些地区地球上部电导率分布的横向变化特征

这篇文章给出了内蒙、甘肃、宁夏某些地区的大地电磁测深结果。这些结果显示出,在地震带内地壳和上地幔电导率分布的横向变化是很明显的,而在比较稳定的鄂尔多斯地台内电导卑分布的横向变化较小,地壳和上地幔的电导率结构比较简单。基于这些结果发现,上部地壳(10—15公里)低阻层在大震震中区较为发育,这可能与大震有关。我国西北地区已有的大地电磁测深结果显示出,莫霍面並不是一个清楚的电性分界面,而某些地区的地壳中部低阻层似乎与康拉界面相应,这就使得后者成为一个明显的电性分界面。