模拟退火方法在三维速度模型地震波走时反演中的应用

采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。     

Travel time computations using a compact eikonal equation for vertical transverse isotropic media

Eikonal solvers often have stability problems if the velocity model is mildly heterogeneous. We derive a stable and compact form of the eikonal equation for P‐wave propagation in vertical transverse isotropic media. The obtained formulation is more compact than other formulations and therefore computationally attractive. We implemented ray shooting for this new equation through a Hamiltonian formalism. Ray tracing based on this new equation is tested on both simple as well as more realistic mildly heterogeneous velocity models. We show through examples that the new equation gives travel times that coincide with the travel time picks from wave equation modelling for anisotropic wave propagation.     

Near-bankfull floods in an Alpine stream: Effects on the sediment mobility and bedload magnitude

In a mountain environment, the transport of coarse material is a key factor for many fields such as geomorphology, ecology, hazard assessment, and reservoir management. Despite this, there have been only a few field investigations of bedload, in particular using multiple monitoring methods. In this sense, attention has frequently focused on the effects of “high magnitude/low frequency floods” rather than on “ordinary events”. This study aims to analyze the sediment dynamics triggered by three high-frequency floods (recurrence interval “RI” between 1.1 and 1.7 yr) that occurred in the Rio Cordon basin during 2014. The flood events were investigated in terms of both sediment mobility and bedload magnitude. The Rio Cordon is an Alpine basin located in northeastern Italy. The catchment has a surface area of 5 km², ranging between 1763 and 2763 m above sea level. The Rio Cordon flows on an armored streambed layer, with a stable step-pool configuration and large boulders. Since 1986, the basin has been equipped with a permanent station to continuously monitor water discharge and sediment flux. To investigate sediment mobility, 250 PIT-tags were installed in the streambed in 2012. The 2014 floods showed a clear difference in terms of tracer displacement. The near-bankfull events showed equal mobility conditions, with mean travel distance one order of magnitude higher than the below-bankfull event. Furthermore, only the near-bankfull events transported coarse material to the monitoring station. Both events had a peak discharge up to 2.06 m³ s^-1, but the bedload transport rates differed by more than one order of magnitude, proving that under the current supply-limited condition, the bedload appears more related to the sediment supply than to the magnitude of the hydrological features. In this sense, the results demonstrated that near-bankfull events can mobilize large amounts of material for long distances, and that floods of apparently similar magnitude may lead to different sediment dynamics, depending on the type and amount of sediment supply.     

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

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

蒸发波导对GNSS海面反射信号有效散射区域的影响

蒸发波导在海洋低空的发生概率高达90%,对舰船雷达、通信等电磁系统具有重要影响.为了分析利用GNSS卫星海面反射信号的时延-相关功率波形反演蒸发波导的可行性,本文提出了GNSS卫星海面反射信号的有效散射区域的概念,并将有效散射区域内的GNSS反射信号区分为GNSS标准反射信号和GNSS波导反射信号;然后,利用射线追踪方法,仿真分析了GNSS卫星海面反射信号的有效散射区域大小对蒸发波导的关键参数——波导高度的敏感性,并分析了在时延-相关功率波形上利用反射信号的传播时延将二者分离开的可行性.结果表明,GNSS卫星海面反射信号的有效散射区域对蒸发波导高度非常敏感,对于2～25m高的GNSS反射信号接收天线,当蒸发波导的高度由0m增加至20m时,GNSS反射信号有效散射区域半径的均值可由约14km迅速扩展至约160km;采用高度角足够大的GNSS卫星可以将有效散射区内的GNSS标准反射信号与GNSS波导反射信号在时延-相关功率波形上分离开.     

地层格架正则化约束下的二维立体层析反演

通过把地层格架信息作用于立体层析Fréchet导数矩阵,使得更新后的速度模型呈现出符合地质规律的块状特征.地层格架信息基于立体层析反演中得到的反射点位置进行非规则B样条插值拟合得到,因此在反演中它将会随着反射点位置的更新自然得到更新.与前人提出的保边缘层析算法或多层立体层析算法相比,本文提出的地层格架正则化无需引入混合正则化项或定义某种复杂的混合速度格式,更为直接也更容易实现.理论和实际数据算例证实了该正则化技巧的稳健性和可靠性,能够得到与实际地质构造特征更为一致的地质一致性反演结果.     

西准噶尔红山地区晚古生代赞岐岩锆石U Pb年代学、地球化学特征及其地质意义

中亚造山带西部西准噶尔地区红山花岗岩体内部发育多期似岩墙状安山质暗色条带,LA-ICP-MS锆石U-Pb同位素定年和岩石化学与同位素分析表明,它们具有不同的形成年龄、相同的源区和相似的形成过程。其中,具有不规则状或环状形态的暗色条带,其锆石U-Pb年龄为319.1±2.9 Ma和313.3±2.4 Ma,远大于红山岩体花岗岩锆石结晶年龄(305～301Ma),可能是红山岩体侵位过程中所捕掳的围岩;具有线性展布特征的安山质暗色条带,其锆石U-Pb年龄为295±2Ma,形成于红山岩体侵位之后,构成伸展岩墙群。红山岩体中的安山质暗色条带和线状岩墙群具有相似的岩石化学组成,富SiO_2(56.48%～63.09%)、MgO(3.56%～6.31%),具有高的Mg#值(51.74～62.40)及Na_2O/K_2O值(1.34～3.43);球粒陨石标准化稀土元素(REE)配分模式呈明显的右倾型,富集轻稀土元素(LREE)和大离子亲石元素(LILE)Rb、K、Ba、U、Sr,亏损高场强元素(HFSE)Th、Nb、Ce、P和重稀土元素,具较弱的负铕异常;其同位素组成特征为(~(87 )Sr/~(86 )Sr)i=0.703295～0.703620,(~(143) Nd/~(144) Nd)i=0.512612～0.512618,εNd(t)=6.91～7.62,(~(206) Pb/~(204) Pb)t为17.6883～17.9876,(~(207) Pb/~(204) Pb)t为15.5313～15.5686,(~(208)Pb/~(204)Pb)t为37.4460～38.0581。它们具有与赞岐岩类似的地球化学特征,总体表现出与弧岩浆作用相关的地球化学特征,可能具有共同的物质来源,为准噶尔洋板片俯冲消减后同一地幔源区在不同阶段的产物。其中,形成于弧岩浆作用时期的赞岐岩(319～313 Ma),构成与岛弧花岗岩类伴生的环状似岩墙状安山质暗色条带;形成于后造山伸展岩浆作用晚期的赞岐岩(～295Ma),构成与达拉布特左行走滑作用相关的陆内伸展岩墙群。安山质暗色条带(岛弧火山作用)、红山岩体(后造山伸展)和线状岩墙群(陆内伸展)记录了西准噶尔红山地区洋陆转换的全过程。     

Bedrock erosion due to hoeing as tillage technique in a hilly agricultural landscape,southwest China

Tillage on hillslopes may not only induce severe soil erosion, but may also cause bedrock erosion under certain conditions. Yet, little is known about bedrock erosion by tillage in a hilly agricultural landscape, southwest China. The aim of this study is to quantify the translocation of rock fragments derived from bedrock fragmentation by hoeing under different conditions, including slope gradient, hoeing depth and soil-covered thickness using a gravel tracing method. The reliability of the gravel tracing method was confirmed by the bedrock dyeing tracing method. Hoeing depth is a significant factor affecting the translocation rate of rock fragments (Q_r ). Meanwhile, under the condition of overlying soil layers (0.06−0.10 m thick), the values of Q_r were significantly smaller with a reduction of 20.7−25.6%, compared with rock fragmentation by hoeing for bare bedrock. However, slope gradient was found to have insignificant effects on Q_r . Fractured bedrock moved as individual small fragments, which was mainly controlled by the hitting force of the hoe, while soil moved in the shape of lumps, which was dominated by both drag force of the hoe and gravity. This study suggests that hoeing into soil-covered bedrock can diminish bedrock erosion while providing soil matrix for shallow soil layers. Our work presents a quantitative assessment of bedrock erosion by hoeing and an underlying insight into characteristics of bedrock erosion by tillage operations in hilly agricultural regions with mudstone and shale, southwest China. © 2020 John Wiley & Sons, Ltd.     

地球自由振荡的保结构模拟

大地震等诸多激励均能激发全球自由振荡现象,通常表现为驻波形式的全球整体振荡.现有的地震波数值模拟方法多为非保结构方法,无法压制长时程计算中的积累误差.本文采用优化的三阶辛格式谱元法,对地球自由振荡及全球尺度的地震波传播进行了长时程模拟.通过与传统的基于Newmark算法的谱元方法结果对比分析,明确验证了本文所得优化的三阶辛格式谱元法在模拟地球自由振荡等大规模长时程问题上的优越性和准确性.上述进展在方法论层面为今后探测、刻画全球尺度地球非均匀结构的驻波数值方法奠定了部分基础,并为相关研究领域提供了新的选择.     

复杂山地随机介质GMM-ULTI法射线追踪

对复杂山地介质的非均质性以及介质中地震波运动学特征进行深入研究,对于提高复杂山地区域地震勘探的效果有着重要的理论意义和实际价值.为了研究复杂山地非均质性和该介质中地震波的一些运动特性,提出了一种复杂山地随机介质的建模方法和一种新的射线追踪算法.与常规算法相比,复杂山地随机介质的生成方法采用更贴近实际介质特点的梯度介质作为背景介质,并在模型生成过程中加入地形修正步骤;新提出的GMM-ULTI射线追踪算法,充分融合群推进法、迎风思想、走时插值法的优势,采用先计算走时后追踪射线路径的两步策略完成射线追踪.算法分析与计算实例表明:复杂山地随机介质的生成方法能灵活、精细且更贴近实际地刻画复杂山地介质的非均质特点;新射线追踪算法兼顾精度和效率、能无条件稳定且灵活地适应复杂山地随机介质的特点;同时基于对几个模型试算结果的分析也得出了复杂山地随机介质中的地震波的一些传播规律.