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

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

河北省南部及邻区地壳P波速度结构三维反演

1980—2012年河北省及邻区测震台网地震记录,使用了河北省南部及邻区（34.0°—38.0°N,112.0°—118.0°E）63个固定地震台站和4 540个地震事件,得到27 709条P波到时数据,采用速度结构与地震位置联合反演的方法,获得研究区内地壳P波三维速度结构,重新确定中小地震震源位置。速度结构揭示：研究区域内地壳的P波速度结构存在明显的横向不均匀性,在10—25 km深度上横向不均匀性更加显著;大地震基本发生在速度异常体或高低速交界区域。地震重新定位结果显示：地震P波走时均方根残差（RMS）从1.68 s降到0.82 s;地震呈明显条带状分布,震源深度与地质构造年代具有一定负相关性。     

滇西北地区主要断裂带及邻区三维P波速度结构的双差地震层析成像

利用区域固定台站和中国地震科学探测台阵记录的7349个近震事件的60471个Pg波绝对走时和196465个相对走时数据,采用双差地震层析成像获得滇西北地区(25°～28.2°N,99.5°～101.5°E)横向分辨率为0.2°的中、上地壳的三维P波速度模型,重点分析了区域内各主要断裂带及其邻区的速度结构特征.结果表明:金沙江—红河断裂带北段15 km以上的P波速度较低,重定位后中、小地震的震中主要位于低速异常的内部,且震源深度在断裂带两侧相似;推测金沙江—红河断裂带作为川滇菱形块体西南边界的剪切控制作用已弱化,分界能力局部减弱,并且断裂带下方主体的低速异常可能为跨越断裂的动力传递提供条件.丽江—小金河断裂带(西南段)两侧存在大范围的P波低速异常,推测此低速体可能是青藏高原东南向挤出的物质,而程海断裂带以东从近地表至25 km深处明显的P波高速异常体则可能会阻挡高原物质东南向的逃逸.     

青藏高原东北缘地壳及上地幔顶部速度结构研究

本文利用青藏高原东北缘71个固定台站与418个流动台站记录到的天然地震事件资料,采用双差层析成像方法对近震走时数据进行反演,获得了研究区高分辨率的三维P、S波速度结构和地震重定位结果.研究结果表明,本文给出的P、S波速度模型较已有的全球模型能更好的解释体波走时与面波相速度观测资料.松潘—甘孜和祁连构造带下方20～40 km深度范围表现为显著的P、S波低速异常,其中松潘—甘孜地块的壳内低速层可能与地壳部分熔融有关,而祁连构造带的壳内低速层则可能与地壳增厚有关.精定位后的岷漳6.7级地震和九寨沟7.0级地震震源深度都位于脆性的上地壳.两个地震的震源区地处不同块体的边界,均处在高、低速过渡带.震源区的壳内低速层可能处于部分熔融或易于蠕变的状态,脆性上地壳更容易积累应变能,从而导致地震的发生.     

利用拟牛顿法和信赖域法联合反演震中分布与一维速度结构

本文将拟牛顿法与信赖域法相结合的反演方法用于地震震中和一维速度结构的联合反演. 通过数值试验比较发现：用拟牛顿法构造Hessian矩阵,虽然仅利用目标函数的一阶导数信息, 但由于保证了正定性和拟牛顿条件, 比解析求得Hessian矩阵更加高效稳定, 通过与信赖域方法的结合,使得该联合反演方法与经典的联合反演方法相比不仅具有较高的计算效率和较高的反演精度,同时也具有更强的总体收敛性.利用本文方法对台湾集集地震重新定位结果（纬度2386°N, 经度12080°E, 深度10km, 发震时刻1999年9月21日17:47:160）与前人主流结果一致.此外,该方法对不同初始速度模型以及不同的地震到时数据集均可给出较一致的结果,进一步表明该方法的优越性.本文还利用1998～2001年间华北遥测地震台网40个台站记录到的64次地震的971条P直达波和首波到时资料,采用本文提出的方法吩北京及其邻区的一维地壳P波速度结构进行联合反演, 结果表明该地区浅层地壳速度较低,中部地壳速度较高, 下部地壳速度逐渐递增.     

环渤海地区的地震层析成像与地壳上地幔结构

利用环渤海地区的天然地震P波到时资料,采用纬度和经度方向分别为05°×06°的网格划分,反演了该地区地壳上地幔的三维P波速度结构.初步结果表明,环渤海地区地壳上地幔的速度结构具有明显的横向不均匀性：京津唐地区地壳中上部的速度异常反映了浅表层的地质构造特征,造山带和隆起区对应于高速异常,坳陷区和沉积盆地对应于低速异常;地壳下部出现大规模的低速异常与华北地区广泛存在的高导层相对应,估计与壳内的滑脱层和局部熔融、岩浆活动有关;莫霍面附近的速度异常反映了地壳厚度的变化及壳幔边界附近热状态的差异;上地幔顶部大范围的低速异常可能是上地幔软流层热物质大规模上涌所致.     

澳大利亚大陆上地幔顶部的速度结构:Pn走时层析成像

澳大利亚大陆的地震速度结构多以面波层析成像、地震折射、接收函数和噪声层析成像为主。澳大利亚大陆的三维速度结构已被总结成澳大利亚地震参考模型(AuSREM)。Pn波走时层析成像可提供整个该大陆更强约束的上地幔结构信息。本文调查了自1993年以来发生在澳大利亚大陆内部大的地震事件的波形和地震震相报告,这些地震事件被固定台站和流动台站所记录。为克服澳大利亚东西部大的横向速度对比和保证震源参数的一致性,所有的地震事件在三维澳大利亚参考速度模型中进行了地震重定位。地震重定位后发现,地震震相的走时残差特征具有良好的一致性。定义Pn波为震中距1.8~15°范围内的初至波。利用FMTOMO方法反演走时差以获得上地幔顶部的P波速度结构,其分辨率为3°×3°。反演结果表明,整个澳大利亚大陆上地幔顶部P波速度呈强烈的横向不均匀性。P波的最大速度为8.36km/s,位于澳大利亚西部和中部的前寒武纪克拉通下方,其比全球ak135模型的P波速度快3.99%;最慢的P波速度为7.66km/s,其位于澳大利亚大陆东缘显生宙造山带之下,较全球ak135模型慢4.74%,这可能与晚第三纪火山活动有关。     

沧州地区地壳上地幔的远震P波三维速度结构成像

利用架设在沧州地区的20台流动地震台站于2006年12月至2010年7月间记录的271个远震事件,读取的2308个P波到时数据,采用地震层析成像方法反演得到沧州及其邻区(38.0°N～39.0°N,116.5°E～117.5°E)的地壳上地幔P波三维速度结构。层析成像结果表明,沧东断裂两侧的地壳介质的速度分布表现出明显的横向差异,浅层速度分布同地表地质结构分布相一致。沧东断裂西北侧沧县隆起的地壳速度较高,表明其基底抬升;断裂带东南边的黄骅拗陷速度较低,说明基底埋藏较深。本文的远震层析成像研究结果和前人使用重力、电磁和人工地震的探测结果都表明,沧东断裂带两侧的地质构造和地球物理性质有明显的变化,这种构造差异在整个地壳中都有体现。     

中国边缘海及邻区地壳和上地幔速度结构特征

本文利用中国地震台网和世界地震台网1 222个地震台记录到的27 909个研究区内地震的479 707条P波到时数据和4 995个研究区外远震事件的90 530条P波到时数据,采用多波多震相地震层析成像法重建了中国边缘海及其邻近地区（10°~45°N,105°~145°E,深至1 061 km）的三维速度图像,揭示了中国边缘海及其邻近地区地壳和上地幔三维速度结构特征.     

京津唐地区中上地壳三维P波速度结构

根据华北遥测台网和延庆-怀来数字地震台网记录的京津唐地区(113.5～119.5E，38.0～41.5N)1992～1999年共16 048条P波首波到时资料，利用近震层析成像方法反演了该地区中、上地壳三维P波速度结构，并用所得到的地震波三维P波速度模型同时得到了震级ＭＬ1.7～6.2的1 132次地震的精确震源参数．反演结果揭示了京津唐地区地震P波速度结构明显的横向不均匀性．      

An approach of refraction seismology for processing and modeling of local earthquake seismogram sections of virtual sources at multiple depths in seismogenic regions—application to Koyna-Warna region,India, for upper crustal P and S velocity structure

We present an approach based on controlled source seismology (CSS) methods, especially developed for processing and modeling of the local earthquake seismograms. Record sections of the local earthquake seismograms generated for multiple source depths illuminate the upper crustal velocity structure in the region. Extensive travel times and synthetic seismograms modeling of the observed record sections reveal the P and S velocity structure in the region. The strength of this approach essentially lies with the possibility of validating the upper crustal velocity models inferred in various subregions of the seismogenic region. A redundant and significantly large number of virtual source local earthquake seismogram sections, gathered for multiple source depths and varying source mechanisms in each of the subregions, validate the same set of P and S velocity models in that region. Further, those models are found to generate the synthetic seismograms consistent with the observed sections. The proposed approach effectively utilizes a reliable dataset from a great volume of well-located local earthquake recordings of a state-of-the-art digital seismograph network. Such a dataset of local earthquake seismograms in the Koyna-Warna active earthquake zone is used here to demonstrate this approach and obtained subregion-specific models of upper crustal P and S velocity structure in the epicentral region. The results indicate that the technique presented here is efficient for processing and modeling the local earthquake seismograms and deriving upper crustal velocity models in the seismogenic regions.     

Crustal and upper mantle structure of the Kozani-Grevena area obtained by non-linear inversion of P and S travel times

The present study focuses on the P and S crustal and uppermost mantle velocity structure in the broader Kozani-Grevena area. The velocity structure is derived from the inversion of travel times of local events. The main data source is the travel times from the aftershock sequence of the large event of 13 of May 1995 (M_w = 6.6) which occurred in the study area. An appropriate preconditioning of the final linearized system is used to reduce ray density effects on the results. An attempt is made to interpret the features and details of the crustal structure in terms of the geotectonic setting of the area. The observed features of the deeper crustal and uppermost mantle structure are in very good agreement with previous results. Specifically, a crustal thickening is observed along a ENE-WSW direction, perpendicular to the well-known Dinaric trend (NNW-SSE) of the geological formations of the area, in accordance with the theoretical expectation of a thicker crust under the accretion prism which starts at the SW edge of the study area.     

上海及其邻近地区三维地壳结构层析成像

利用上海及其邻近地区地震台网的地震记录及人工爆破资料,采用地震层析成像方法反演研究区地壳三维P波、S波速度结构模型。成像结果表明:不同深度的P波、S波速度扰动呈北西向展布特征;地震大多数分布在低速区的块体内,说明地震不仅与断裂活动有关,有可能还与物质结构和性质有关;上海及其邻近地区上地壳速度结构,断裂活动和物质性质三者之间存在密切关系。     

Path-specific Velocity Structure of Western China from Surface-wave Dispersion

We develop one-dimensional (1-D) path-specific velocity models in western China using new Rayleigh and Love wave group and phase velocity dispersion measurements for 20 events in the region. The earthquakes were grouped into three geographic clusters from which we compute the average phase and group velocity dispersion. We invert the average dispersion curves simultaneously for 1-D shear-velocity models appropriate for the three central Asian paths, using three previous shear-velocity models as initial models. The models are validated by forward modeling waveforms of recent events. The crustal thickness beneath western China in the vicinity of the Lop Nor test site is 50–60 km and our velocity models are consistent with major geologic features (e.g., basins and mountain ranges) and previous structural models for this region.     

Finite-frequency tomography in a crustal environment: Application to the western part of the Gulf of Corinth

In this paper we investigate finite-frequency effects in crustal tomography. We developed an inversion procedure based on an exact numerical computation of the sensitivity kernels. In this approach we compute the 3D travel-time sensitivity kernels by using (1) graph theory and an additional bending to estimate accurately both rays and travel-times between source/receiver and diffraction points and (2) paraxial ray theory to estimate the amplitude along theses rays. We invert both the velocity and the hypocentre parameters, using these so-called banana-doughnut kernels and the LSQR iterative solver. We compare the ray-theoretical and the finite-frequency tomography to image the intermediate structures beneath the Gulf of Corinth (Greece), which has long been recognized as the most active continental rifting zone in the Mediterranean region. Our dataset consists of 451 local events with 9233 P- first-arrival times recorded in the western part of the Gulf (Aigion area) in the framework of the 3F-Corinth European project. Previous tomographic images showed a complex velocity crustal model and a low-dip surface that may accommodate the deformation. Accurate velocity models will help to better constrain the rifting process, which is still a subject of debate. The main results of this study show that finite-frequency tomography improves crustal tomographic images by providing better resolved images of the 3D complicated velocity structure. Because the kernels spread the information over a volume, finite-frequency tomography results in a sharpening of layer boundaries as we observed for the shallower part of the crust (down to 5 km depth) beneath the Gulf of Corinth.     

利用双差层析成像方法反演青藏高原东南缘地壳速度结构

本文利用云南及周边区域地震台网2010—2016年记录到的近震资料,采用双差层析成像方法进行地震重定位并获得了青藏高原东南缘的三维地壳速度结构。结果显示:重定位后的震源位置精度得到明显提高,震源主要分布于20 km深度以上的中上地壳;地震分布与速度结构存在一定的相关性,大多数地震发生在中上地壳的低速异常区内以及高、低速异常区域之间;研究区上地壳速度结构存在明显的横向不均匀性,其速度异常与地表地形及地质特征密切相关;中下地壳分布着两条主要的低速带,一条沿着安宁河断裂、小江断裂分布在川滇菱形地块的东侧;另一条主要分布在川西北次级地块内,并穿过丽江断裂向南延伸,推测这两条低速带可能是青藏高原中下地壳物质向南逃逸的两条通道。      

UK 1-D regional velocity models by analysis of variance of P-wave travel times from local earthquakes

A method is presented for deriving 1-D velocity depth models from earthquake bulletin data. The models can be used as initial models for more advanced modelling techniques such as tomographic inversion. The method is useful when there is little or no refraction and long-range reflection survey data. The bulletin travel times are subjected to an analysis of variance, where they are separated into source, distance, and receiving station terms. The distance terms describe the variation of travel time with distance, and the associated trend lines allow 1-D velocity models for the crustal layers to be determined. The velocity models provide an average crustal model for the region derived from local data. This does not include superficial layers which are necessarily poorly determined. Earthquake bulletin P-wave data from propagation paths across three different regions of the UK are employed to illustrate the use of the technique.     

曲线坐标系下的块体运动与应变模型研究

本文建立了顾及地球扁率和局部切标架随点变化特性的椭球坐标系下的刚体运动模型和块体运动与应变模型,以及球坐标系下顾及局部切标架随点变化特性的严密的块体运动与应变模型,分析了球坐标系下块体运动与应变模型及椭球坐标系下的块体运动与应变模型间的差异;通过计算具体讨论了地球扁率和曲线坐标系的局部切标架随点变化特性对欧拉矢量与应变张量的影响.结果表明：地球扁率对刚体欧拉矢量和应变参数的影响甚小,具体计算时可以不予考虑,但曲线坐标系的局部切标架随点变化特性对两者的影响较大,在建模过程中需要顾及,常用的Savage模型需要修正.     

Tomographic imaging of P wave velocity structure beneath the region around Beijing

The three-dimension crustal and upper mantle structures in the region around Beijing were studied by seismic tomography. We used the P wave arrival times from local and teleseismic events. These events were recorded by 250 stations of the North China Seismic Array and 108 stations of the Beijing Telemetry Seismic Network. 118 869 P wave arrivals from 10 285 local events and 12 189 P wave arrivals from 107 teleseismic events were used in the inversion. We obtained the 3-D P wave velocity structure of the cru...     

利用双差成像方法反演青藏高原东北缘及其邻区地壳速度结构

利用区域尺度双差层析成像方法,使用2009年1月至2017年2月的近震资料,对青藏高原东北缘及其邻区内记录到的地震事件进行震源位置和三维速度结构的联合反演.重新定位后震源空间位置得到明显改善,浅层的地震波速与地形和沉积层厚度对应较好,研究区地震主要发生在河西走廊过渡带的低泊松比区域.本文将研究区分为五个区域并分别对其层析成像结果进行了讨论,结果显示研究区不同地块之间地壳结构变化明显,地壳物质整体呈酸性,青藏高原东北缘地壳增厚可能主要发生在中下地壳.