祁连山-河西走廊地壳速度结构及速度与电性的联合解释

联合利用天然地震和人工地震测深资料对祁连山-河西走廊地区进行地壳三维速度结构反演,获得了该区地壳三维速度结构图像．给出了按5km间隔的中、上地壳速度切面和地壳速度结构与莫霍界面等深线图．综合大地热流实测结果和深部岩石物性等,建立了下地壳部分熔融的韧性模型．给出了计算地震波速度和电阻率的方法与公式,进行了部分熔融状态下深部电性与地震波速度的联合数值模拟和联合解释,其结果与实测的基本一致．分析了各构造单元之间下地壳电阻率、速度与热流的关系,对下地壳部分熔融程度进行了估计．     

南北地震带北段速度结构与电阻率特性

建立了南北地震带北段下地壳部分熔融几何形状和分布的模型。熔融的几何分布和连通性由以下3个因素决定：①熔融的位置;②熔融所占总体积的比例;③固－固和固－液界面的最小自由能。熔融形状有薄版形、管道形和椭球包体形。用此模型计算出含有部分溶融岩石的体积、剪切弹性模量和电导率的一般表示式。选择参数并计算出南北地震带北段下地壳的Vp、Vs和电阻率ρs值,即在鄂尔多斯地台下地壳的速度为Vp ≈6.71km/s,Vs≈3.96km/s,电阻率为上百至几百欧姆米。而在祁连、秦岭褶皱系下地壳速率Vp约为6．4～6.6km/s,Vs约为3．72～3．89km/s,电阻率为几至几十欧姆米。数值模拟的结果与地球物理反演结果基本相符,并进行了定性的构造解释。用一个模型对下地壳中的速度特征与电性特性进行了联合数值模拟,说明了此模型在解释下地壳中的Vp、Vs和ρs的有效性。     

朝鲜半岛南部三维地壳速度结构成像

对多震相走时层析成像方法进行改进和优化, 使用朝鲜半岛陆地和近海的地震记录, 对朝鲜半岛南部(345°N～39°N，125°E～130°E)进行地震波走时反演， 获得了本研究区的三维速度结构成像结果.分析了沉积层、基底、上地壳、康腊面、下地壳和莫霍面的特征.研究了不同地质区的深部结构的主要差别、主要断裂特征、 海陆交界的地壳特征等.可以发现京畿山地与沃川褶皱带、太白山褶皱带、临津江褶皱带深部结构的差异尤为明显.     

三峡水库坝址及邻区中上地壳P波速度结构

利用三峡地震台网的到时资料,采用地震层析成像技术重建了三峡水库坝址及邻区中上地壳P波三维速度图像,深部构造研究结果是：（1）黄陵背斜岩性均匀完整、固结程度高,高速异常到20km深度依然可见,形状随深度变化．（2）秭归盆地沉积厚度不均匀,最厚处约6km．（3）结晶基底顶面在0～7km、底面在14～16km．（4）黄陵背斜与秭归盆地之间及附近区域和培石附近区域深部存在上地幔部分熔融物质,现今构造发展与演化与上地幔熔融物质冷却固化过程密切相关．（5）九湾溪断裂东西两侧地区地壳深部构造的速度结构不同,是三峡地区重力梯度带成因需要考虑的因素．（6）库水渗透作用对地壳浅层速度结构产生影响．     

青藏高原东北缘岩石圈三维密度结构

综合重力观测资料和地震波走时资料反演了青藏高原东北缘岩石圈三维密度结构,并对该区岩石圈结构及动力学特征进行了讨论.首先利用收集到的P波近震和远震走时数据进行地震层析成像,得到研究区岩石圈三维P波速度结构.然后利用速度-密度经验关系式,将速度扰动转化为密度扰动建立研究区三维初始密度模型.最后利用分离的布格重力异常反演得到了岩石圈三维密度结构.反演结果表明:青藏高原东北缘地壳密度结构特征有利于地震孕育发生和地壳物质侧向流动;地壳内,密度异常等值线走向与地表断裂走向基本一致,进入地幔后,密度异常等值线走向发生了顺时针旋转,这表明青藏高原东北缘地壳和地幔具有不同的构造运动模式,暗示该区可能发生了壳幔解耦;80~100 km深度上,P波速度异常较密度异常明显偏低,推测该区可能发生了部分熔融或者岩石含水量的增加;印度板块俯冲和周围坚硬块体阻挡联合作用,使得青藏高原东北缘形成了强大的区域构造应力场,并导致深部软流圈热物质上涌,为该区壳幔解耦、部分熔融和P波速度降低创造了条件.     

青藏高原地震波三维速度结构的研究

根据青藏高原及其邻区的模拟地震台站和宽频带数字地震记录资料,采用区域体波层析成你和瑞利面波层析成像,反演得到了青藏高原地区的三维地震波速度结构。两种层析成像方法得到的地壳上地幔P波和S波速度结构的结果非常相似,它们显示,青藏高原南部的拉萨块体的上地壳存在明显的低速区,青藏高原北部的羌塘地区的下地壳和上地幔顶部相对速度较低。这些结果与青藏高原南、北部处于碰撞过程中的不同阶段有关。     

藏南岩石圈导电性结构与流变性——超宽频带大地电磁测深研究结果

为了比较全面、客观地认识青藏高原岩石圈深部结构,探讨高原岩石圈形变特征、应变状态、热结构、板块（或地体）运动和壳、幔物质流变等重要科学问题,必须研究高原内部各个地质构造区域壳、幔电性结构沿东西方向的变化;INDEPTH—MT在喜马拉雅．西藏南部地区完成了6条超宽频带大地电磁深探测剖面研究．通过这6条剖面的电性结构成像,讨论了研究区地壳和上地幔导电性三维结构特点;发现西藏南部沿东西方向超出1000km范围,较普遍存在中、下地壳高导层,这高导层并不完全是连续的,向雅鲁藏布江大拐弯处高导层变薄、变浅、电阻率升高．讨论了藏南岩石圈的流变性问题,认为藏南中、下地壳具有良导电性,可以证明西藏巨厚的地壳中确实存在部分“熔融体”和“热流体”,藏南巨厚的中、下地壳的物质状态是热的、软弱的、塑性的,甚至可能是“流变”的;结合岩石物理实验结果的讨论认为,与藏南大地电磁资料相适应的地壳部分熔融百分比应能达到5％～14％;对于地壳中的细晶岩来说,在这个熔融百分比下引起黏度的降低量有可能达到引起地壳“流变性”的要求：但对于花岗岩来说,也许不足以引起地壳产生“流变”．     

利用福建爆破实验记录检验华南速度模型

基于2010~2012年"福建及台湾海峡地壳深部结构陆海联合探测"项目中18次爆破实验的准确发震时间、发震位置所确定的多震相走时数据,进行爆破震相特征分析,采用华南速度模型对单纯型、Hyposat、Hypo2000及Locsat等4种定位方法进行比较,得出最优定位方法。利用14炮爆破记录反演得到的速度模型与华南速度模型作比较,结果表明,反演模型的上地壳厚度小于华南模型的,下地壳厚度大于华南模型的,上地壳速度与华南模型的相当,下地壳速度大于华南模型的。利用本研究的反演结果对其他4次爆破进行重新定位,定位精度明显提高。     

云南腾冲火山区上部地壳三维地震速度层析成像

根据1999年腾冲火山区人工地震探测资料，用层析成像方法反演了腾冲火山区上地壳三维P波速度结构．层析成像结果显示出，在腾冲火山区下方的上地壳内存在明显的P波低速区，该低速区在地下7～8km深度，可能与上地壳内的岩浆囊或部分熔融体有关．层析结果还显示出低速区主体在热海地区的东北侧．地表附近的低速区可能是构造破裂发育地区形成了有利于深部热流体上升的通道，因此，在热海形成最活跃的地热显示．      

台湾地区的壳幔P波速度和Vp／Vs波速比结构研究

利用中国台湾地震台网1992-2004年间记录到的P波和S波到时数据,使用地震层析成像方法联合反演了台湾及邻近地区地壳上地幔的P波速度和Vp/Vs波速比结构。与已有结果相比,采用的球坐标系有限差分算法适合于地壳波速结构强烈不均匀的情况,能够提高走时计算和三维射线追踪的准确程度,而联合反演得到的P波速度、Vp/Vs波速比及其重新定位后的地震分布为揭示壳幔结构的横向变化和深部特征提供了更多约束。结果表明：在地壳浅表层,沉积盆地和造山带具有明显不同的波速特征,沉积盆地P波速度偏低、Vp/Vs波速比偏高,造山带P波速度偏高、Vp/Vs波速比偏低;台东纵谷作为欧亚大陆和菲律宾海板块的碰撞缝合带,Vp/Vs波速比明显偏高,并延伸至中、下地壳和上地幔,反映了岩石破裂、部分熔融和流体活动对壳幔深部物质的影响;台湾东北部琉球岛弧附近的P波速度、Vp/Vs波速比和重新定位后的地震分布勾画出了菲律宾海板块俯冲至欧亚大陆下方的形态,高速的海洋俯冲板片及板片上方P波低速和高Vp/Vs波速比的特征一直延伸至80km;沿着东西向的纵向波速剖面,地壳厚度在中央山脉东侧达到60km左右,并呈现向东倾的趋势,反映出欧亚大陆板块的向东俯冲和挤压使造山带中下地壳物质变形增厚。     

青藏高原东缘中下地壳流与地壳变形

地壳缩短导致青藏隆升造山是普遍的认识.然而,在青藏东部,越来越多的观测数据和研究支持了中下地壳流与隆升造山的关系.目前,地壳缩短造山机制和中下地壳流造山机制仍然处于争论之中.本文建立了二维黏弹塑性有限元模型,模拟了龙门山断层带的多个地震循环的应变与变形,探讨了无与有中下地壳流情况下,地壳地表的位移、速度与变形的分布和演化;以及有中下地壳流情况下,不同流动范围、速度与黏度对模型结果的影响;并结合地形变观测数据的约束,推测了青藏东缘中下地壳流的流动状态.模拟结果显示,通过对比有和无中下地壳流的模拟结果,发现青藏东部震间的地表垂向速度在变形样式及数值上存在较大差异,即存在地壳流的地表垂向抬升速率显著大于无地壳流;震间在龙门山断层西侧附近产生的垂向凸状隆起随中下地壳流的速度、黏度及通道长度的变化而变化.此外,本文研究结果对青藏其他地区可能存在的地壳流的研究也具有一定的参考意义.     

青藏高原东部和周边地区地壳速度结构的背景噪声层析成像

利用连续地震背景噪声记录和互相关技术获得瑞利面波格林函数,进而反演获得了青藏高原东部和周边地区的地壳三维速度结构.地震数据源于北京大学宽频带流动观测地震台阵,国家数字测震台网数据备分中心提供的部分固定台站的连续记录及INDEPTH IV宽频带流动观测地震台阵.首先对观测数据进行处理和分析取得所有可能台站对的面波经验格林函数和瑞利波相速度频散曲线,反演得到了观测台阵下方周期从6~60s的瑞利波相速度异常分布图像.并且进一步反演获得研究区域三维剪切波速度结构和莫霍面深度分布.短周期(6~14s)相速度异常分布与地表地质构造特征吻合较好,在青藏高原和四川盆地之间存在一个明显的南北向转换带.而本文最重要的结果是周期大于25s的相速度异常分布图像显示,以昆仑断裂带为界,柴达木盆地和祁连山脉地区呈现与青藏高原截然不同的中地壳速度结构,反而与青藏高原东缘地区和川滇菱形块体速度结构相似.反演获得的剪切波速度在27.5~45km深度的切片也明显地揭示:青藏高原的松潘—甘孜地块和羌塘地块呈现均一的低速层;然而,柴达木盆地和祁连山脉地区则呈现较强的横向不均匀性,尤其是柴达木盆地的高速异常和四川盆地的高速异常相对应.这些结果为前人提出的青藏高原东北向台阶式增长模式提供了重要的地震学观测证据.与全球一维平均速度模型(AK135)相比较发现,本文测量和反演获得的研究区域内平均相速度和剪切波速度都比AK135模型慢很多,尤其是青藏高原的中地壳(25~40km)剪切波速度显著低于全球平均速度模型.进一步的层析成像反演证实松潘—甘孜和羌塘地块中地壳(27.5~45km)呈现大范围均一的低速层,为青藏高原可能存在大规模中下地壳"层流"提供地震学观测证据.在祁连山脉的27.5~45km深度观测到的明显低速异常体可能对应于该造山带下地幔岩浆活动导致的底侵作用,表明引起该地区地壳增厚的主要机制可能是来自地幔岩浆的底侵作用.     

Azimuthal anisotropy of Rayleigh waves beneath the Tibetan Plateau and adjacent areas

The crustal and upper mantle azimuthal anisotropy of the Tibetan Plateau and adjacent areas was studied by Rayleigh wave tomography. We collected sufficient broadband digital seismograms trav-ersing the Tibetan Plateau and adjacent areas from available stations, including especially some data from the temporary stations newly deployed in Yunnan, eastern Tibet, and western Sichuan. They made an adequate path coverage in most regions to achieve a reasonable resolution for the inversion. The model resolution tests show that the anisotropic features of scope greater than 400 km and strength greater than 2% are reliable. The azimuthal anisotropy pattern inside the Tibetan Plateau was similar to the characteristic of tectonic partition. The crustal anisotropy strength is greater than 2% in most re-gions of East Tibet, and the anisotropy shows clockwise rotation surrounding the eastern Himalayan syntaxis. Vertically, the anisotropy direction indicates a coherent pattern within the upper crust, lower crust, and lithosphere mantle of the Tibetan Plateau, which also is consistent with GPS velocity field and SKS fast polarization directions. The result supports that the crust-mantle deformation beneath the Tibetan Plateau is vertically coherent. The anisotropy strength of crust and lithospheric upper mantle in Yunnan outside the Tibetan Plateau is lower than 2%, so SKS splitting from core-mantle boundary to station should largely be attributed to the anisotropy of asthenosphere.     

Analysis of the rising and sinking movement of the crust in Qinghai－Xizang Plateau，China

ＡｎａｌｙｓｉｓｏｆｔｈｅｒｉｓｉｎｇａｎｄｓｉｎｋｉｎｇｍｏｖｅｍｅｎｔｏｆｔｈｅｃｒｕｓｔｉｎＱｉｎｇｈａｉ－ＸｉｚａｎｇＰｌａｔｅａｕ，ＣｈｉｎａＺｉ－ＬｉｎｇＬＵ（吕梓龄），ＺｈｅｎｇＺＨＵＡＮＧ（庄真），Ｚｈｕ－ＷｕＦＵ（傅竹武），Ｊｉａ－...     

南北地震带岩石圈S波速度结构面波层析成像

本文利用天然地震面波记录和层析成像方法,研究了南北地震带及邻近区域的岩石圈S波速度结构和各向异性特征.结果表明南北地震带的东边界不但是地壳厚度剧变带,也是地壳速度的显著分界.其西侧中下地壳的S波速度显著低于东侧,强震大多发生在低速区内部和边界.青藏高原东缘中下地壳速度显著低于正常大陆地壳,在松潘甘孜地块和川滇地块西部大约25～45 km深度存在壳内低速层;这些低速特征与高原主体的低速区相连,有利于下地壳物质的侧向流动.地壳的各向异性图像与下地壳流动模式相符,即下地壳物质绕喜马拉雅东构造结运动,东向的运动遇到扬子坚硬地壳阻挡而变为向南和向北东的运动.面波层析成像结果支持青藏高原地壳运动的下地壳流动模型.南北地震带的岩石圈厚度与其东侧的扬子和鄂尔多斯地块相似但速度较低.川滇西部地块上地幔顶部(莫霍面至88 km左右)异常低速;松潘甘孜地块上地幔盖层中有低速夹层(约90～130 km深度).岩石圈上地幔的速度分布图像与地壳显著不同,在高原主体与川滇之间存在北北东向高速带,可能会阻挡地幔物质的东向运动.上地幔各向异性较弱且与地壳的分布图像显然不同.因此青藏高原岩石圈地幔的构造运动具有与地壳不同的模式,软弱的下地壳提供了壳幔运动解耦的条件.     

Inversion ofQ value structure beneath the Tibet Plateau

A total of 11 earthquakes with 15 Rayleigh wave paths, recorded at 11 broadband digital PASSCAL seismometers installed in the Tibet Plateau by the Sino-U.S. joint research group, were used to determine the phase velocity and attenuation coefficient of surface waves in periods of 10–130 s. The average shear wave velocity and quality factor {ie271-1} structures in the crust and upper mantle were obtained in this region. The result shows the average {ie271-2} is low and there exists a high attenuation ({ie271-3}=93–141) layer in the crust. The depth range of the low {ie271-4} value layer (16–42 km) is consistent with the range of low velocity layer (21–51 km) in the crust. Below 63 km in the lower crust, {ie271-5} decreases with depth from 114 to 34 at depth of 180 km. The low shear wave velocity and low value of {ie271-6} at the same depth range in the crust indicate that the rocks in the range is probably melted or partially melted. According to the shear wave velocity structure, the average thickness of the crust is about 71 km and a clear velocity discontiniuty appears at the depth of 51 km. The low-velocity zone (4. 26 km/s) at depth of 96–180 km may be corresponding to the asthenosphere.     

Research on the 3—D Seismic Structures in Qinghai—Xizang Plateau

Based on the recording data from the analogue and broadband digital seismic stations in and around Qinghai-Xizang (Tibet)Platean,the three dimensiomal 3-D) seismic velocity stroctures in Qinghai-Xizang Plateau were obtained by using the regional body wave tomography and surface wave tomography.The results from these two tomography methods have similar characteristics for P-and S-wave velocity structures in crust and upper mantle.They show that there are remarkahle low velocity zones in the upper crust of L hasa block in the southern Qinghai-Xizang Plateau and the lower crust and upper mantle of Qiangtang block in the northern Qinghai-Xizang Plateau.These phenomena may be related to the different steps of collision process in southern and northern Qinghai-Xizang Plateau.     

Inversion of Q value structure beneath the Tibet Plateau

A total of 11 earthquakes with 15 Rayleigh wave paths, recorded at 11 broadband digital PASSCAL seismometers installed in the Tibet Plateau by the Sino-U.S. joint research group, were used to determine the phase velocity and attenuation coefficient of surface waves in periods of 10–130 s. The average shear wave velocity and quality factor {ie271-1} structures in the crust and upper mantle were obtained in this region. The result shows the average {ie271-2} is low and there exists a high attenuation ({ie271-3}=93–141) layer in the crust. The depth range of the low {ie271-4} value layer (16–42 km) is consistent with the range of low velocity layer (21–51 km) in the crust. Below 63 km in the lower crust, {ie271-5} decreases with depth from 114 to 34 at depth of 180 km. The low shear wave velocity and low value of {ie271-6} at the same depth range in the crust indicate that the rocks in the range is probably melted or partially melted. According to the shear wave velocity structure, the average thickness of the crust is about 71 km and a clear velocity discontiniuty appears at the depth of 51 km. The low-velocity zone (4. 26 km/s) at depth of 96–180 km may be corresponding to the asthenosphere. Contribution No. 96A0047, Institute of Geophysics, SSB, China. This study was supported by the National Natural Science Foundation of China.     

Preliminary study of crust-upper mantle structure of the Tibetan Plateau by using broadband teleseismic body waveforms

As part of a joint Sino-U.S. research project to study the deep structure of the Tibetan Plateau, 11 broadband digital seismic recorders were deployed on the Plateau for one year of passive seismic recording. In this report we use teleseimic P waveforms to study the seismic velocity structure of crust and upper mantle under three stations by receiver function inversion. The receiver function is obtained by first rotating two horizontal components of seismic records into radial and tangential components and then deconvolving the vertical component from them. The receiver function depends only on the structure near the station because the source and path effects have been removed by the deconvolution. To suppress noise, receiver functions calculated from events clustered in a small range of back-azimuths and epicentral distances are stacked. Using a matrix formalism describing the propagation of elastic waves in laterally homogeneous stratified medium, a synthetic receiver function and differential receiver functions for the parameters in each layer can be calculated to establish a linearized inversion for one-dimensional velocity structure. Preliminary results of three stations, Wen-quan, Golmud and Xigatze (Coded as WNDO, TUNL and XIGA), located in central, northern and southern Plateau are given in this paper. The receiver functions of all three stations show clear P-S converted phases. The time delays of these converted phases relative to direct P arrivals are: WNDO 7.9s (for NE direction) and 8.3s (for SE direction), TUNL 8.2s, XIGA 9.0s. Such long time delays indicate the great thickness of crust under the Plateau. The differences between receiver function of these three station shows the tectonic difference between southern and north-central Plateau. The waveforms of the receiver functions for WNDO and TUNL are very simple, while the receiver function of XIGA has an additional midcrustal converted phase. The S wave velocity structures at these three stations are estimated from inversions of the receiver function. The crustal shear wave velocities at WNDO and TUNL are vertically homogeneous, with value between 3.5–3.6 km/s down to Moho. This value in the lower crust is lower than the normal value for the lower crust of continents, which is consistent with the observed strong Sn attenuation in this region. The velocity structure at XIGA shows a velocity discontinuity at depth of 20 km and high velocity value of 4.0 km/s in the midcrust between 20–30 km depth. Similar results are obtained from a DSS profile in southern Tibet. The velocity under XIGA decreases below a depth of 30 km, reaching the lowest value of 3.2 km/s between 50–55 km. depth. This may imply that the Indian crust underthrusts the low part of Tibetan crust in the southern Plateau, forming a “double crust”. The crustal thickness at each of these sites is: WNDO, 68 km; TUNL, 70 km; XI-GA, 80 km. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, Supp., 581–592, 1992.