Iteratively re‐weighted and refined least squares algorithm for robust inversion of geophysical data

A robust metric of data misfit such as the ?₁‐norm is required for geophysical parameter estimation when the data are contaminated by erratic noise. Recently, the iteratively re‐weighted and refined least‐squares algorithm was introduced for efficient solution of geophysical inverse problems in the presence of additive Gaussian noise in the data. We extend the algorithm in two practically important directions to make it applicable to data with non‐Gaussian noise and to make its regularisation parameter tuning more efficient and automatic. The regularisation parameter in iteratively reweighted and refined least‐squares algorithm varies with iteration, allowing the efficient solution of constrained problems. A technique is proposed based on the secant method for root finding to concentrate on finding a solution that satisfies the constraint, either fitting to a target misfit (if a bound on the noise is available) or having a target size (if a bound on the solution is available). This technique leads to an automatic update of the regularisation parameter at each and every iteration. We further propose a simple and efficient scheme that tunes the regularisation parameter without requiring target bounds. This is of great importance for the field data inversion where there is no information about the size of the noise and the solution. Numerical examples from non‐stationary seismic deconvolution and velocity‐stack inversion show that the proposed algorithm is efficient, stable, and robust and outperforms the conventional and state‐of‐the‐art methods.     

用接收函数方法研究汶川和芦山地震之间未破裂段的地壳结构

芦山与汶川地震之间存在约40 km的地震空区.震源区和地震空区的深部构造背景的研究对深入了解中强地震的深部孕育环境及地震空区的地震活动性具有重要科学意义.利用本小组布设的15个临时观测地震台以及21个芦山科考台站和21个四川省地震局固定台站记录的远震数据,用H-K叠加方法得到各个台站的地壳厚度和平均泊松比,并构建了接收函数共转换点(CCP)偏移叠加图像以及反演得到台站下方的S波速度模型.我们的结果揭示了震源区和地震空区地壳结构特征差异:(1)汶川震源区的地壳平均泊松比为~0.28;芦山震源区为~0.29;而地震空区处于泊松比变化剧烈的区域;(2)汶川地震与芦山地震的震源区以西下方的Moho面呈现深度上的突变(这与前人的研究成果基本一致),分别从~44 km突变到~59 km,~40 km突变到~50 km,而地震空区地壳平均厚度呈现渐变性变化;(3)地震空区Moho面下凹且具有低速的上地壳.综合一维S波速度结构和H-k以及CCP的初步结果,这可能显示汶川地震的发震断裂在深部方向上向西倾斜并形成切割整个地壳的大型断裂;芦山地震则可能是由于上、下地壳解耦引起的;而地震空区处于两种地震形成机制控制区域的过渡带中.     

中国东南部P波速度结构与构造分析

利用中国东南部地震台站的波形数据,通过远震P波层析成像方法反演了福建及台湾地区的上地幔P波速度结构,据此分析了华夏块体的构造属性并探讨了中国东南部上地幔的深部动力学机制.结果表明,研究区上地幔速度结构存在明显的横向非均匀性,它们与区域构造的深部动力学成因密切相关.宏观上,NE向展布的低速异常与东南沿海地区中新生代火成岩的条带状分布保持了较好的一致性,推测其上地幔可能存在幔源物质上涌;以政和—大埔断裂带为界,东、西华夏块体的速度结构存在差异,结合前人的研究结果分析,东、西华夏块体岩石圈减薄和中新生代岩浆活动的驱动机制有所不同,西华夏块体为“板块俯冲+岩浆底侵(岩石圈拆沉)”的动力学模式,东华夏块体是“玄武质岩浆迁移+岩浆底侵”的动力学模式.     

Establishing an effective offset-dependent velocity model by ray tracing and its application in Kirchhoff time migration

Conventional Kirchhoff prestack time migration based on the hyperbolic moveout can cause ambiguity in laterally inhomogeneous media, because the root mean square velocity corresponds to a one-dimensional model under the horizontal layer assumption; it does not include the lateral variations. The shot/receiver configuration with different offsets and azimuths should adopt different migration velocities as they contribute to a single image point. Therefore, we propose to use an offset-vector to describe the lateral variations through an offset-dependent velocity corresponding to the difference in offset from surface points to the image point. The offset-vector is decomposed into orthogonal directions along the in-line and cross-line directions so that the single velocity can be expressed as a series of actual velocities. We use a simple Snell's law-based ray tracing to calculate the travel time recorded at the image point and convert the travel time to an equivalent velocity corresponding to a pseudo-straight ray. The double-square-root equation using such an equivalent velocity in the offset-vector domain is non-hyperbolic and asymmetrical, which improves the accuracy of the migration. Numerical examples using the Marmousi model and a wide azimuth field data show that the proposed method can achieve reasonable accuracy and significantly enhances the imaging of complex structures.     

地震动残差分析盆地附加放大效应:日本关东盆地为例

我国地震动预测及地震危险性分析通常仅考虑局部场地浅层岩土层对地震动的放大效应,不能考虑较大范围的地质条件影响,如沉积盆地厚沉积层对地震动的附加放大效应通常被忽略,造成盆地内地震动及地震危险性预测结果普遍被低估。本文以地震动观测记录数据充足的日本关东盆地为例,采用地震动残差分析方法评估盆地附加放大效应,分析覆盖层厚度、盆地内空间位置、震级、震源距对地震动放大效应的影响,建立关东盆地附加放大效应经验评估模型。分析表明:关东盆地附加放大效应与反应谱周期相关,整体上从短周期的1.0逐渐增大至周期为5s时的1.5,附加放大效应与覆盖层厚度相关性较小,主要受盆地空间位置和震源距的影响;盆地北部边缘及西北部地区附加放大效应更强烈,盆地南部附加放大效应较小,这可能与盆地边缘效应密切相关。本文建立的关东盆地附加放大效应经验模型略高于BSSA14和ASK14模型的放大效应预测。相关研究结果可用于我国地震动预测、下一代地震动区划图修订等。     

基于三维模型中的射线追踪方法测定鲁甸M_s6.5地震震源深度

2014年8月3日云南鲁甸发生M_S6.5地震,造成重大人员伤亡和财产损失。此次地震发震构造复杂,引起了地震学界的广泛关注。本文基于三维建模方法,建立鲁甸地区三维层状非均匀速度模型,采用逐段迭代射线追踪方法对鲁甸M_S6.5地震进行三维射线追踪走时计算;利用震中距150km范围内的近震Pg波震相走时数据,通过射线追踪走时拟合获得鲁甸地震的震源深度约为12km,与前人研究成果基本一致,表明了本文采用的三维建模和射线追踪方法的有效性。     

Design response spectra for Aqaba City–Jordan

Strong motion records for earthquakes epicentred in the Gulf of Aqaba during swarms in August 1993 and November 1995 are studied. These recordings are made at two sites in the City of Aqaba, one near the seashore and the other inland. The former site has deep alluvial profile with low shear-wave velocity, while the latter's soil profile is constituted of shallow stiff alluvium overlaying rock. The magnitudes of these events range from 4·8 to 6·2 on the Richter scale. Ground motion on the seashore site is much greater than the inland site. Response spectra are derived for each of the available accelerograms. Thirteen records at the inland site are statistically treated to arrive at mean response spectrum using 5 per cent damping. The derived spectra are found to be sharp with small plateau, this is associated with small peak displacement. Empirical site-dependent response spectra are also derived based on recent research and on recommendations of the UBC'97. Comparison is made between the empirical and measured spectra. These were found to be quite close for the inland site but data are not sufficient for the seashore site to make a conclusive comparison. Copyright © 1999 John Wiley & Sons, Ltd.     

北极海冰减退引起的北极放大机理与全球气候效应

自20世纪70年代以来,全球气温持续增高,对北极产生了深刻的影响。21世纪以来,北极的气温变化是全球平均水平的2倍,被称为"北极放大"现象。北极海冰覆盖范围呈不断减小的趋势,2012年北极海冰已经不足原来的40%,如此大幅度的减退是过去1 450年以来独有的现象。科学家预测,不久的将来,将会出现夏季无冰的北冰洋。全球变暖背景下北极内部发生的正反馈过程是北极放大现象的关键,不仅使极区的气候发生显著变化,而且对全球气候产生非常显著的影响,导致很多极端天气气候现象的发生。北极科学的重要使命之一是揭示这些正反馈过程背后的机理。北极放大有关的重大科学问题主要与气—冰—海相互作用有关,海冰是北极放大中最活跃的因素,要明确海冰结构的变化,充分考虑融池、侧向融化、积雪和海冰漂移等因素,将海冰热力学特性的改变定量表达出来。海洋是北极变化获取能量的关键因素,是太阳能的转换器和储存器,要认识海洋热通量背后的能量分配问题,即能量储存与释放的联系机理,认识淡水和跃层结构变化对海气耦合的影响。全面认识北极气候系统的变化是研究北极放大的最终目的,要揭示气—冰—海相互作用过程、北极海洋与大气之间反馈的机理、北极变化过程中的气旋和阻塞过程、北极云雾对北极变化的影响。在对北极海冰、海洋和气候深入研究的基础上,重点研究极地涡旋罗斯贝波的核心作用,以及罗斯贝波变异的物理过程,深入研究北极变化对我国气候影响的主要渠道、关键过程和机理。     

OIQC技术在雷达反演VAD廓线资料退模糊中的应用研究

多普勒天气雷达在强风速天气条件下经常出现的速度模糊误差,严重限制了风场信息的广泛应用,如反演VAD廓线等。本文将最优插值客观分析质量控制技术(OIQC),应用于雷达反演VAD廓线资料中速度模糊的质量控制。对用于反演VAD风廓线的那些格点圈层上的资料,以简化的多变量统计插值技术(Cressman插值),为被检查的资料提供一个对比分析的参考值,从而确定折叠次数,修正模糊速度。根据多普勒雷达径向速度资料的特点,水平方向采用非各向同性的单变量背景场误差协方差函数,垂直方向协方差函数取为高斯型。修正模糊速度后运用VAD技术得到水平风的垂直廓线。对强龙卷天气过程的数值试验结果表明,OIQC技术应用于雷达资料质量控制,无需引入外界风场,有能力修正模糊速度。结合VAD方法,可以从存在严重折叠的雷达资料中反演出较为准确的风廓线。     

登陆台风维持和暴雨增幅实例的能量学分析

对登陆后迅速消亡的８１１６号台风和登陆后长期维持不消的８４０７号台风的动能平衡作对比分析，结果表明：对流层上部动能水平通量辐散和穿越等压线的运动使８１１６号台风迅速消亡，其中以正压过程为主；８４０７号台风则从环境大气输入动能和位能，这可能是台风暖性低压长久不消的重要原因。最后讨论了台风动能水平通量辐散与外围暴雨增幅的关系。