Pg回折波的上地壳三维层析成像

提出利用人工地震Ｐｇ回折波走时重建区域上地壳三维速度分布的方法。正问题利用经典最小二乘迭代求解垂向速度梯度模型的参数值及相应走时残差。反问题采用泛函空间的最小二乘准则来实现。对速度不需要预先进行参数化处理，可以计算反演区内任意点的速度值。     

利用三维弯曲界面的反演方法重建唐山震区莫霍界面三维构造形态

研究利用反射波的走时反演三维弯曲界面和介质层速度的计算方法．各层界面利用分段非完全三次多项式描述．正问题采用一种快速的三维射线追踪方法．反演过程采用变阻尼最小二乘法．数值模拟结果表明，解很快收敛到真模型．处理了通过唐山震区的实测资料，重建了该震区莫霍界面三维构造形态，并揭示了该区域构造与地震活动的关系．      

联合利用走时与波形反演技术研究地壳三维速度结构(Ⅱ)——应用

使用阻尼最小二乘法进行震源参数和地壳三维速度结构的走时联合反演．所用资料为Ｓ波和Ｐ波到时差,并用人工地震资料的二维解释结果作为三维速度模型的特定约束条件．为建立初始模型,又利用天然地震构成了准二维剖面．在走时反演基础上,利用遗传算法进行了几个地震事件的波形反演尝试,并对走时反演获得的地壳速度结构模型的局部进行了修正．以３４°～４２°Ｎ,９４°～１１２°Ｅ作为研究区域,在该区域中收集了１９８６年以来大量地震的Ｓ波和Ｐ波到时差资料,７条人工地震二维速度剖面资料和２个数字化地震台的几个地震的三分向记录资料．对这些资料进行了处理,最后得出了０～２５ｋｍ深度不同截面的速度分布,并对所得结果进行了分析．     

二维反射地震成象—界面与速度的联合反演

介绍了利用反射走时联合反演地震界面起伏形态和速度横向分布的方法。正问题采用了一种简单而有效的二维介质中两点射线追踪方法，反问题在广义最小二乘意义下求解，描述界面和速度的两类参数不必要先进行分离，数值模拟的结果表明方法是有效的。处理了通过唐山震区的各庄－丰南－丰宁深地震测深剖面，得到了意义的结果。     

有限差分二,三维速度层析成象

介绍了利用有限差分方法进行地震走时资料速度反的基本原理与方法,并对理论模型进行了正和反演。然后对唐山滦县的Ｐｇ波三维资料和菏泽－林县－－长治的二维Ｐｇ资料分别进行了三维和二维反演,反演结果与浅层地质资料较好吻合。由于反演是采用不用矩阵运算的反投影方法,可以实现稠密采样,这为我们提供了一种较好的得到高精度浅层速度的办法,为深部结果的反演提供了可靠的浅层速度资料。     

三维地质模型中的地震定位初步研究

快速精确地确定震源位置和发震时刻是地震速报及地震预警的前提和基础,可以为地震应急、震后救援及地震趋势预测等提供重要科学依据。传统的定位方法由于通常基于简单一维地球模型,对地下结构复杂区域地震事件定位误差较大。本研究采用块状结构建模方案来描述三维复杂地质体,发展了相应的逐段迭代射线追踪方法实现走时的快速正演计算;采用阻尼最小二乘方法,计算了走时对震源三维空间位置的偏导数,为实现三维地质模型中的地震定位研究奠定了基础。     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.     

界面和速度的两步反演-唐山震区三维细结构研究a

研究区域界面和速度的两步反演计算方法．三维界面采用分段非完全多项式描述,三维速度的重建采用泛函空间的最小二乘原理．计算采用两步进行:第1步反演三维界面形态;第2步将剩余走时残差按加权分配的方式进行三维速度反演．处理了唐山滦县震区地震测深资料,获得了唐山震区深部三维构造形态及唐山、滦县震区三维地壳速度分布．结果表明:唐山震区深部三维构造总体为北东走向,北东向的丰台——野鸡坨断裂与北东向的唐山断裂所夹的构造为莫霍界面隆起区. 该区中下地壳三维速度结构沿北东向的唐山断裂存在明显低速异常带,其位置与唐山地震活动带的位置相一致．沙河驿附近较大的低速异常块体对应于较密集的地震分布．该区域存在一条北西向的高速异常带,可能是一条隐伏深断裂．唐山7.8级地震震中区下地壳为北东向的低速异常带与北西向的高速异常带相交处．这两组构造对控制唐山7.8级地震的孕育和发生起了重要作用．      

分块三维速度模型生成及理论地震图的计算

本文提出了在计算机上实现分块三维地壳模型及利用加权最小二乘拟合生成平缓光滑的三维速度函数的方法,给出了适用于分块、块内速度连续变化的三维模型中Cauchy射线追踪的新算法,简介了基于上述方法反射线的基本理论所编制的合成三维理论地震图的程序包RSSGTD.给出的两个盆地状模型的算例表明,所使用的模型生成方法具有模拟复杂地壳结构的能力;与三维样条函数方法比较,最小二乘拟合方法能给出更加适合射线方法合成地震图计算的速度函数,并且内存小、计算速度快;所给出的Cauchy射线追踪算法能够适合块状模型中任何体波射线的追踪.     

黄土高原及邻区地壳P波速度结构

利用地震波走时联合反演算法(改进型最短路径算法)进行三维弯曲地震射线追踪正演,以及共轭梯度法求解带约束的阻尼最小二乘问题进行反演,同时更新速度模型和地震震中位置,结合地方震和区域地震走时资料得到了黄土高原(含汾渭断陷盆地)及邻区地壳三维P波速度结构.其横向变化结果表明,研究区地壳内的P波高速异常区与其内的地震活动构造带相一致,地震多发生在P波高速异常区的边缘或高、低速异常区的交汇处.秦岭山区和鄂尔多斯块体东南区为P波低速异常区.而垂向变化结果则表明研究区存在低速异常区.     

THE 3-D VELOCITY STRUCTURE OF CRUST AND UPPERMOST MANTLE AND ITS TECTONIC IMPLICATIONS IN FUJIAN PROVINCE

It is important to detect the fine velocity structures of the crust and uppermost mantle to understand the regional tectonic evolution, earthquake generation processes, and to conduct earthquake risk assessment. The inversion of uppermost mantle velocity and Moho depth are strongly influenced by crustal velocity heterogeneity. In this study, we collected first arrivals of Pg and Pn and secondary arrivals of Pg wave from the seismograms recorded at Fujian provincial seismic network stations. New 3-D P-wave velocities were inverted by multi-phase joint inversion method in Fujian Province. Our results show that the fault zones in Fujian Province have various velocity patterns. The shallow crust is characterized by high velocity that represents mountains, while the mid-lower crust shows low velocities. The anomalous velocities are correlated closely with tectonic faults in Fujian Province. Velocity anomalies mainly show NE-trending distribution, especially in the mid-lower crust and uppermost mantle, which is consistent with the NE-trending of the regional main fault zones. Meanwhile, a part of velocity patterns show NW trending, which is related to the secondary NW-oriented faults. Such velocity distribution also shows a geological structural pattern of "zoning in east-west direction and blocking in north-south direction" in Fujian area. In the crust, a low velocity zone is found along Zhenghe-Dapu fault zone as mentioned by previous study, however our result shows the low velocity exists at depth of 20~30km in mid-lower crust. Compared with previous study, this low velocity zone is larger and deeper both in range and depth. The crustal thickness of 28~35km from our joint inversion is similar to the results from the receiver functions of previous studies. The thinnest crust(28km)is observed at offshore in the north of Quanzhou; while the thickest crust(35km)is located west of Zhangzhou near the Zhenghe-Dapu fault zone. Generally, thinner crustal thickness is found in offshore of Fujian Province, and thicker crustal thickness is in the mainland. However, we also found that crustal thickness becomes thinner along the east side of Yongan-Jinjiang Fault. The values of Pn velocities in the region vary from 7.71 to 8.26km/s. The velocity distribution of the uppermost mantle presents a large inhomogeneity, which is correlated with the distribution of the fault zone. High Pn velocity anomalies are found mainly along the west side of the Zhenghe-Dapu fault zone(F₂), and the east side of the Shaowu-Heyuan fault zone(F₁), which is strip-shaped throughout the central part of Fujian. Low Pn velocity anomalies are observed along the coast and Taiwan Straits, including the Changle-Zhaoan fault zone, the coastal fault zone, and the Fuzhou Basin. We also found a low Pn velocity anomaly zone, which extends to the coast, in the Shaowu-Heyuan fault zone at the junction of the Fujian, Guangdong and Jiangxi Provinces. In the west of Taiwan Straits, both high and low Pn velocity anomalies are observed. Our results show that the historical strong earthquakes(larger than magnitude 6.0) are mainly distributed between positive and negative anomaly zones at different depth profiles of the crust, and similar anomalies distribution also exists at the uppermost mantle, suggesting that the occurrence of strong earthquakes in the region is not only related to the anomalous crustal velocity structure, but also affected by the velocity anomaly structure from the uppermost mantle.     

琼东北及雷州半岛深地震测深剖面Pg波资料处理结果——用地震波前成像和射线数分布分析法

本文用地震波前成像和射线数分布分析法对上世纪80年代在琼东北及雷州半岛地区完成的三条宽角反射/折射地震剖面的Pg波资料进行了处理,得到了三条测线下方的上地壳顶部P波速度结构及基底的空间展布形态特征.结果显示研究区域表现出较复杂的上部地壳结构特征,浅部P波速度结构的高、低速区域分别与地质构造上的一些凸起和凹陷相对应,用射线数分布分析法所显示的基底结构及形态特征与地震剖面穿过区域的地质构造有较好的对应关系,测线穿过的一些断裂下方P波速度结构表现为强烈的横向非均匀性.     

唐山震区上地壳平均波速比分析

收集河北地震台网记录的唐山震区2012-2016年ML≥1.0地震,选取震源深度5-10km、波形信噪比较高的波形数据,拾取直达P波、S波震相到时,利用和达法计算平均波速比,并分析研究区平均波速比水平分布特征与时变特征,结果表明:唐山震区上地壳平均波速比约为1.71,且水平分布呈区域化特征,反映了研究区复杂的地壳结构。研究区内3次ML≥4.0地震前,波速比曲线均呈较显著的下降-低值-回升形态,表明该区平均波速比变化可能具有一定地球物理异常特征。     

青藏高原中东部地壳和上地幔顶部P波层析成像

为获取青藏高原中东部地壳和上地幔顶部的精细结构,本文基于1万4 484条天然地震的P波（Pg和Pn）到时数据,对青藏高原中东部地壳和上地幔顶部进行P波三维速度结构层析成像,获取了该区域内地壳P波、上地幔顶部Pn波的速度结构和地壳厚度信息。层析成像结果显示,青藏高原中东部地壳P波速度范围为5.2—7.2 km/s,上地幔顶部Pn波速度范围为7.7—8.4 km/s,地壳厚度范围为48.0—68.6 km,地壳和上地幔顶部存在强烈的横向不均匀性,与地质块体分布有较好的对应关系。地壳P波速度结构显示,研究区中、下地壳分布有较大范围的低速区,上地壳与中下地壳P波分布存在明显的差异:羌塘地块和巴颜喀拉地块在上地壳主要表现为高速异常,随着深度增加逐渐表现为低速异常;而柴达木地块在上地壳主要表现为低速异常,下地壳则表现为高速异常;柴达木地块和拉萨地块在上地幔顶部表现为较高的Pn波速度,最高约为8.4 km/s,而在巴颜喀拉地块和羌塘地块东部,Pn波总体上表现为低速,最低约为7.7 km/s。研究区内地壳厚度的总体特征表现为南厚北薄,其中羌塘地块东部和拉萨地块的地壳较厚,而柴达木地块和巴颜喀拉地块东部的地壳相对较薄,羌塘地块西部存在局部的地壳变薄现象,反映了印度板块对欧亚板块北向俯冲作用下的岩石圈变形特征。      

Imaging 3-D crustal P-wave velocity structure of western Yunnan with bulletin data

Western Yunnan is a region with intensive tectonic activity and serious earthquake risk. It is of significant importance to study three dimensional crustal structure of this region to understand the tectonic setting and disaster mechanism. Densification and digitalization of seismic networks in this region provides an opportunity to study the velocity structure with bulletin data. In this study, we collect P-wave data of 10 403 regional earthquakes recorded by 79 seismic stations from January 2008 to December 2010. In addition to first arrivals data (Pg with epicentral distance less than 200 km and Pn), the Pg (or P) data with epicentral distance more than 200 km are also considered as later direct arrivals in the tomographic inversion. We also compare the quantity and the quality of the seismic data before 2010 and after 2010. The test results show that adding the follow-up Pg phase can effectively improve the inversion ability of crustal imaging, and quantity and the data quality are significantly improved since 2010. The tomographic results show that: (1) The Honghe fault zone, which is the major fault systems in this region, may cut through the entire crust, and the velocity contrasts between two sides at lower crust beneath the Honghe fault are estimated at higher than 10%, while the velocity difference below Nujiang fault zone extends only in the upper crust; (2) Most of the earthquakes in the region occurred at the interface of high-velocity media and low-velocity media, i.e., the areas with high velocity gradient, which has been validated in other areas.     

地壳三维构造反演和速度层析成像

本文给出利用空间深地震测深反射资料重建地壳三维构造和速度分布的方法，适用于任意分层并且每层可由若干断块组成的地壳模型。通过采用适当的反演方法获得地壳界面的三维构造，并在此结果基础之上，进一步将剩余走时残差归因于地壳速度的不均匀性，采用东分块的反演方法，重建地壳三维速度分布，数值模拟的结果表明算法是有效的。     

Tomographic determination of 3-D crustal structure — Joint inversion of explosion and earthquake data

In this paper, the theory and method, obtaining the tomographic determination of three-dimensional velocity structure of the crust by use of the joint inversion of explosion and earthquake data, are given. The velocity distribution of the crust is regarded as a continuous function of the spatial coordinates without parametrization of the velocity model ahead, so that the inversion solution would not be influenced by different parametrization procedures. The expressions of integration kernels, which relates the two kinds of data sets, are also given. The authors have processed the observed data in Tangshan earthquake region by the method proposed in this paper, and obtained the tomographic results of the middle and upper crust structures in this region. The comparison of these results with the result obtained only by the explosion data, has also been made.     

Tomographic determination of 3－D crustal structure──Joint inversion of explosion and earthquake data

Ｔｏｍｏｇｒａｐｈｉｃｄｅｔｅｒｍｉｎａｔｉｏｎｏｆ３－Ｄｃｒｕｓｔａｌｓｔｒｕｃｔｕｒｅ──Ｊｏｉｎｔｉｎｖｅｒｓｉｏｎｏｆｅｘｐｌｏｓｉｏｎａｎｄｅａｒｔｈｑｕａｋｅ　ｄａｔａＸｉａｎ－ＫａｎｇＺＨＡＮＧ（张先康）；Ｚｈｕｏ－ＸｉｎＹＡＮＧ（杨卓...     

应用层析成像技术研究华北地壳速度结构

利用华北地区１６９个地震台站记录的１９８４￣１９９５年的１６９６个地震的Ｐ波到时资料,重建了华北地区地壳三维速度图像,不同深度的速度图像表明：地壳中地震波速度的横向不均匀性显著,上地壳的速度分布与地有已知构造特征相联系;太行山山前断裂带是切割到地壳底部的深断裂带,２５ｋｍ深度普遍存在速度间断面;地壳底部起伏明显;强震孕育发生与其周围介质的不均匀性和土地幔的隆起或凹陷密切相关。     

河套地震带Pg波速度过渡区的地震空间相关长度幂律变化

基于Pg波速度反演和地震重新定位, 运用单键群算法对Pg波速度过渡带的地震空间相关长度进行了幂律拟合分析． Pg波速度反演结果表明, 其速度的横向变化表现出构造相依的特征, 速度高低与地壳厚度呈正相关, 并在包头—西山嘴凸起和岱海凹陷两个区域形成Pg波速度过渡带． 利用重新定位的地震数据计算了这两个速度过渡区的地震空间相关长度, 结果显示其幂律拟合曲线均呈一定的增长趋势, 表明2008年以来两个Pg波速度过渡区域的应力作用不断集聚、 增强, 断层有逐步进入协同化阶段的可能, 加之速度过渡带通常是地壳运动强烈区域, 未来将成为孕育中强地震的有利场所． 在有效控制定位误差的条件下, 重新定位可以明显减小地震空间相关长度的离散形态, 提高计算精度．