三维地质模型中地震波共轭梯度非线性走时反演

地震体波走时层析成像是探测地球内部速度结构的重要方法之一。基于三维块状建模以及三角形拼接的界面描述方式,结合快速高效的逐段迭代射线追踪方法,获得三维复杂地质模型中的地震射线路径与走时信息,采用共轭梯度非线性反演算法,进行地震波走时反演。实验结果表明共轭梯度反演算法在三维层状模型中具有较高的有效性。     

三维复杂层状介质中多震相走时联合反演成像

采用新近提出的改进型不规则最短路径多次波射线追踪正演技术,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,讨论了三维复杂层状模型中利用多震相走时资料进行联合反演成像的技术方法.考虑到不同震相种类走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数;另外,考虑到同时反演速度模型和反射界面起伏中不同参数变化对...     

二维复杂层状介质中地震多波走时联合反演成像

采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法.     

三角网格剖分下速度与反射界面的同时反演

为了解决复杂速度模型中的走时正、反演问题,例如:含不规则起伏地表、不规则地下波阻抗界面、以及不规则速度异常体的复杂地学模型,本文采用三角网格单元模型参数化下的分区多步改进型最短路径算法,实现了多震相地震射线的追踪计算,结合共轭梯度法求解带约束的阻尼最小二乘反演问题,实现了多震相走时联合同时反演成像的方法技术.当界面起伏较大时出现散射,从而造成散射点所在区域射线密度过密,导致该区域内速度和界面的过度更新.为了克服上述问题,我们在同时反演中引入了射线密度的概念,从而有效地解决了上述过度更新问题.数值模拟实验表明:采用三角网格单元进行模型参数化,可保证在复杂模型中的正演计算具有较高的计算精度;同时反演中可以准确地刻画不规则异常体和不规则反射界面.因此,本文提出的走时成像方法技术具有较广的实用价值.     

基于多震相联合反演江苏地区一维速度模型研究

利用地震走时数据,采用联合反演方法获取了江苏地区的一维P波速度模型。与仅采用初至波走时的传统天然地震走时获取方法相比,该方法充分利用了大量存在的续至波参与反演,能有效改进中下地壳的反演能力。针对地震震相目录中常存在震相标识错误的问题,采用的自动判别筛选震相方法能最大限度提高数据走时的精度,可以对不同震相进行有效区分。与其他常用一维速度模型相比,本文反演的模型对Pg、Pn震相走时拟合效果最佳,残差最小。当所用走时数据拥有较高定位精度时,该反演方法能为研究区三维速度结构成像和地震定位提供较可靠的一维速度模型。     

三维介质中速度结构和界面的联合成像

根据波逆行原理推导了三维介质中地震波射线走时对界面偏导数的完整基本关系式,并对基本关系式进行简化,得到其在二维和一维介质中的关系式. 给出了任意多个复杂界面情况下,反演时所需的走时对界面偏导数系数矩阵. 为了检验三维介质中速度和界面联合成像理论的有效性,进行了数值模拟计算,很好地对三维速度结构和界面进行了重建.     

近震资料确定三维地震速度结构的一种方法

我们研究了由近震得到的地震走时资料估计地震波速度三维结构的反演方法。这种方法由同时确定速度界面的二维深度分布、三维速度分布、台站校正和震源来表征。速度界面的深度分布和慢度扰动的分布以纬度、经度和深度的幂级数来拟合。这种三维地震结构的表示不仅有助于减少未知参数的个数,也有助于解析地求取界面深度、速度和走时关于未知参数偏导数的值。连接震源和台站之间的射线轨迹,是在具有确切界面深度分布的简化的速度分布条件下由称为试射法的射线追踪方法确定的。数值试验证明,我们的反演方法可成功地估计三维速度结构。     

起伏地形下的高精度反射波走时层析成像方法

全球造山带及中国大陆中西部普遍具有强烈起伏的地形条件.复杂地形条件下的地壳结构成像问题像一面旗帜引领了当前矿产资源勘探和地球动力学研究的一个重要方向.深地震测深记录中反射波的有效探测深度可达全地壳乃至上地幔顶部,而初至波通常仅能探测上地壳浅部.为克服和弥补初至波探测深度的不足,本文基于前人对复杂地形条件下初至波成像的已有研究成果,采用数学变换手段将笛卡尔坐标系的不规则模型映射到曲线坐标系的规则模型,并将快速扫描方法与分区多步技术相结合,发展了反射波走时计算和射线追踪的方法.进而利用反射波走时反演,实现起伏地形下高精度的速度结构成像,从而为起伏地形下利用反射波数据高精度重建全地壳速度结构提供了一种全新方案.数值算例从正演计算精度、反演中初始模型依赖性、反演精度、纵横向分辨率以及抗噪性等方面验证了算法的正确性和可靠性.     

各向异性TI介质qP反射波走时层析成像

地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力.     

地壳和上地幔三维速度与界面结构层析成像——理论部分

提出了一种三维模型参数化方法,使得正问题和反问题可以用解析式表达．地壳和上地幔速度与界面结构分别由一组非常光滑的函数描述．在正问题和反问题中,使用了初值射线追踪方法,走时对模型网格节点参数的偏微商在射线追踪的同时进行了解析计算．由于速度和界面函数具有二阶连续偏微商,由射线路径的发散和聚焦在地表引起的影区可以避免．射线追踪之后,可以获得由矩阵和矢量组成的反问题．解决反问题使用了具有阻尼因子的奇异值分解法．一组由几条纵剖面组成的合成数据被用来进行数值模拟计算,结果表明方法是健全的．与二维方法比较,在地壳模型非常复杂的情况下,三维反演方法可以给出准确的速度结构和界面位置．      

中国西部地区地壳结构特征与强震活动相关性研究

中国西部地区是地震活动十分强烈的地区,天山、阿尔泰、帕米尔和西昆仑都是著名的地震构造带,在这些地震构造带和周边地区发生了多次震级大于5级的强震.本文通过分析西部地区的重力场特征,根据重力数据结合地震剖面、应用Parker-Oldenburg方法反演得到了研究区莫霍面深度,通过对比地震层析成像的反演结果,分析了研究区的地壳结构特征.计算结果表明,研究区地壳结构不均匀特征明显,在造山带地区一般是莫霍面坳陷区,盆地则是莫霍面隆起区,主要造山带地壳速度结构表现为高速区,盆地和主要凹陷区为低速区.根据计算结果和以往强震震中位置分析了地壳构造与强震活动的相关性,西部地区的地震活动与地壳结构的横向不均匀密切相关,强震主要发生在地壳速度变化带附近和地壳速度结构差异较大的地区,在构造应力作用下,这些地壳介质非均匀地区易发生强震,这是中国西部造山带和盆-山边界附近频发强震的构造原因之一.     

JOINT INVERSION OF SURFACE WAVE DISPERSION AND RECEIVER FUNCTIONS FOR CRUSTAL AND UPPERMOST MANTLE STRUCTURE BENEATH CHINESE TIENSHAN AND ITS ADJACENT AREAS

The Tienshan orogenic belt is one of the most active intracontinental orogenic belts in the world. Studying the deep crust-mantle structure in this area is of great significance for understanding the deep dynamics of the Tienshan orogen. The distribution of fixed seismic stations in the Tianshan orogenic belt is sparse. The low resolution of the existing tomographic results in the Tienshan orogenic belt has affected the in-depth understanding of the deep dynamics of the Tienshan orogenic belt. In this paper, the observation data of 52 mobile seismic stations in the Xinjiang Seismic Network and the 11 new seismic stations in the Tienshan area for one-year observations are used. The seismic ambient noise tomography method is used to obtain the Rayleigh surface wave velocity distribution image in the range of 10~50s beneath the Chinese Tienshan and its adjacent areas (41°~48° N, 79°~91° E). The joint inversion of surface wave and receiver function reveals the S-wave velocity structure of the crust and uppermost mantle and the crustal thickness below the station beneath the Chinese Tienshan area(41°~46° N, 79°~91° E). The use of observation data from mobile stations and new fixed seismic stations has improved the resolution of surface wave phase velocity imaging and S-wave velocity structure models in the study area.
The results show that there are many obvious low-velocity layers in the crust near the basin-bearing zone in the northern Tienshan Mountains and the southern Tienshan Mountains. There are significant differences in the structural characteristics and distribution range of the low-velocity zone in the northern margin and the southern margin. Combining previous research results on artificial seismic profiles, receiver function profiles, teleseismic tomography, and continental subduction simulation experiments, it is speculated that the subduction of the Tarim Basin and the Junggar Basin to the Tienshan orogenic belt mainly occurs in the middle of the Chinese Tienshan orogenic belt, and the subduction of the southern margin of the Tienshan Mountains is larger than that of the northern margin, and the subduction of the eastern crust is not obvious or in the early subduction stage. There are many low-velocity layers in the inner crust of the Tienshan orogenic belt, and most of them correspond to the strong uplifting areas that are currently occurring. The thickness of the crust below the Tienshan orogenic belt is between 55km and 63km. The thickness of the crust(about 63km)is the largest near the BLT seismic station in the Bazhou region of Xinjiang. The average crustal thickness of the Tarim Basin is about 45km, and that of the Junggar Basin is 47km. The S-wave velocity structure obtained in this study can provide a new deep basis for the study of the segmentation of the Tienshan orogenic belt and the difference of the basin-mountain coupling type.     

Constraining the anisotropy structure of the crust by joint inversion of seismic reflection travel times and wave polarizations

In the context of wide-angle seismic profiling, the determination of the physical properties of the Earth crust, such as the elastic layer depth and seismic velocity, is often performed by inversion of P- and/or S-phases propagation data supplying the geometry of the medium (reflector depths) or any other structural parameter (P- or S-wave velocity, density...). Moreover, the inversion for velocity structure and interfaces is commonly performed using only seismic reflection travel times and/or crustal phase amplitudes in isotropic media. But it is very important to utilize more available information to constrain the non-uniqueness of the solution. In this paper, we present a simultaneous inversion method of seismic reflection travel times and polarizations data of transient elastic waves in stratified media to reconstruct not only layer depth and vertical P-wave velocity but also the anisotropy feature of the crust based on the estimation of the Thomsen’s parameters. We carry out a checking with synthetic data, comparing the inversion results obtained by anisotropic travel-time inversion to the results derived by joint inversion of seismic reflection travel times and polarizations data. The comparison proves that the first procedure leads to biased anisotropic models, while the second one fits nearly the real model. This makes the joint inversion method feasible. Finally, we investigate the geometry, P-wave velocity structure and anisotropy of the crust beneath Southeastern China by applying the proposed inversion method to previously acquired wide-angle seismic data. In this case, the anisotropy signature provides clear evidence that the Jiangshan-Shaoxing fault is the natural boundary between the Yangtze and Cathaysia blocks.     

巴颜喀拉块体地壳结构多样性探测

青藏高原内部地壳岩性的改造、岩性随深度变化及形变构造是探索研究地壳增厚、物质运动问题的关键.巴颜喀拉块体位于青藏高原中北部,地域广袤,通过对块体内中、东部不同区域的深地震广角反射/折射震相的综合分析,利用反射率理论地震图方法对不同性质震相走时及振幅特征进行细致的模拟计算,进一步研究巴颜喀拉块体内部不同区域地壳精细结构.结果显示:巴颜喀拉块体地壳厚度50~60 km、整体向西逐渐增厚,结晶地壳平均速度6.07~6.18 km·s~(-1)、岩层速度大幅降低,壳内多强反射界面结构、但不同区域差异明显;东部若尔盖盆地地壳介质速度整体低速、壳幔边界较为清晰;中部玉树一玛多段下部地壳发现约6.8 km·s~(-1)的"高速度"介质结构,壳幔边界不清、被改造为2~4 km厚的高速度梯度层,显示了巴颜喀拉块体内部地壳增厚、介质岩性结构被改造的差异性.地壳内部多组强反射、低视速度走时震相揭示了介质岩性的低速破碎、弱化蠕变以及可能的壳内解耦构造.局部地区下地壳的高视速度震相特征显示了青藏高原地壳改造增厚大背景下可能存在稳定的"原始地壳"结构残留或是与上地幔物质的浸入交流.巴颜喀拉块体内不同区域地壳增厚、岩性结构、结晶基底及壳幔边界性质被改造的多样性为深入认识青藏高原地壳形变及动力学过程带来新的启迪.     

基于深度降噪自编码神经网络的中国大陆地壳厚度反演

本文采用基于数据驱动的深度降噪自编码网络构建了瑞雷面波群速度、相速度频散特性与地壳厚度的正反演函数关系,并利用最新频散模型反演了中国大陆的地壳厚度.对于神经网络架构体系的评价,除了考虑传统意义上的测试误差、训练误差之外,本文还用已知物理原理的正演结果与网络预测结果进行比较;在设计网络构架时,同时考虑地球模型和面波频散的...     

2-D crustal structure of Changbaishan-Tianchi volcanic region determined by seismic traveltime inversion

IntroductionTrial-and-error forward modeling of wide-angle seismic reflection/refraction traveltimes for 2-D velocity structure is extremely time-consuming, even for experienced data interpreters. For wide-angle seismic reflection/refraction experiments that consist of numerous shots along a single line, it is quite difficult through repeated trial-and-error forward modeling to construct a 2-D model that fits the data within acceptable limits (Cerveny, et al, 1977; ZHANG, et al, 200 . In ad…     

远震反演中新的地壳异常体改正方法

近年来，远震走时层析成像方法有了长足的发展，在地下结构反演的研究中获得了众多成果。针对射线在台站下方覆盖率较差而导致远震反演方法对地壳速度约束不足的问题，本文提出了一种新的地壳异常体改正方法——对同一台站的相对走时残差进行求和平均去均值，消除了地壳中复杂的速度异常体对上地幔速度结构反演的干扰，反演结果具有更优的相对走时残差分布，同时反演得到的速度模型具有更小的数据方差。     

南海西沙地块岛屿地震观测和海陆联测初步结果

为了研究南海西沙地块下方的地壳结构,我们在岛屿区架设了流动地震台站进行天然地震观测和海底地震仪的人工地震探测.本文利用远震接收函数方法和射线追踪方法对琛航台的远震数据和海陆联测数据进行了初步处理和分析,建立了琛航岛下方的一维横波速度模型以及横穿琛航岛的二维地壳结构模型.琛航岛顶部存在2 km厚的新生代低速沉积层,下地壳...     

中国中南地区综合地质地球物理研究

利用重、磁和天然地震资料对中南地区进行了综合地质地球物理研究.根据研究区的重、磁异常的分布特征,可将其分为4个重力异常区、3个航磁异常区.根据重力资料反演计算得到的研究区的地壳厚度在295～41 km之间,总的趋势为西厚东薄,地壳厚度与地形起伏基本上呈镜像关系.根据磁力资料计算得到的研究区居里界面在12～40 km之间变化.地震层析成像结果表明研究区内的速度分布总体上体现了纵横交错的断块特征.在纵向尺度上,江汉－洞庭盆地以及周缘造山带的上地壳结构变化不大,中地壳和下地壳则普遍受到现今构造活动的改造,以致岩石的结构发生了一定的变化.它们主要表现为低速区域的扩大,尤其是在地壳下部尤为突出,这与断陷盆地的拉张以及造山带岩石层的底侵和拆沉作用密切相关.通过对研究区地球物理场的分析计算,在研究区共提取主要断裂带34条.根据岩石层板块大地构造理论,依据岩石层结构、地壳结构和结晶基底等深部结构的不同,将研究区中板内不同构造单元——块体作为一级构造单元,块体之间的深大断裂带作为块体的边界——块体结合带,据此原则在研究区中划分出两个一级构造单元,五个二级构造单元.     

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

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