Three‐term inversion of prestack seismic data using a weighted l2, 1 mixed norm

We present a new inversion method to estimate, from prestack seismic data, blocky P‐ and S‐wave velocity and density images and the associated sparse reflectivity levels. The method uses the three‐term Aki and Richards approximation to linearise the seismic inversion problem. To this end, we adopt a weighted mixed l_{2, 1}‐norm that promotes structured forms of sparsity, thus leading to blocky solutions in time. In addition, our algorithm incorporates a covariance or scale matrix to simultaneously constrain P‐ and S‐wave velocities and density. This a priori information is obtained by nearby well‐log data. We also include a term containing a low‐frequency background model. The l_{2, 1} mixed norm leads to a convex objective function that can be minimised using proximal algorithms. In particular, we use the fast iterative shrinkage‐thresholding algorithm. A key advantage of this algorithm is that it only requires matrix–vector multiplications and no direct matrix inversion. The latter makes our algorithm numerically stable, easy to apply, and economical in terms of computational cost. Tests on synthetic and field data show that the proposed method, contrarily to conventional l₂‐ or l₁‐norm regularised solutions, is able to provide consistent blocky and/or sparse estimators of P‐ and S‐wave velocities and density from a noisy and limited number of observations.     

大地电磁三维正演并行算法研究

大地电磁法三维正演算法计算量大,用传统的串行程序计算相当耗时。而三维正演是逐个频率按顺序计算的,并行性好,适合并行运算。结合MPI自身的优越性,在深入分析大地电磁三维正演串行程序实现流程的基础上,确定了并行计算的思路,实现了三维正演的并行计算。通过三个理论模型对实现的三维正演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维正演并行程序运行结果正确,效率提高明显。此思路可为解决其它地球物理超大计算量问题所借鉴。     

A finite element model for seismic response analysis of deformable rocking frames

A new finite element model to analyze the seismic response of deformable rocking bodies and rocking structures is presented. The model comprises a set of beam elements to represent the rocking body and zero‐length fiber cross‐section elements at the ends of the rocking body to represent the rocking surfaces. The energy dissipation during rocking motion is modeled using a Hilber–Hughes–Taylor numerically dissipative time step integration scheme. The model is verified through correct prediction of the horizontal and vertical displacements of a rigid rocking block and validated against the analytical Housner model solution for the rocking response of rigid bodies subjected to ground motion excitation. The proposed model is augmented by a dissipative model of the ground under the rocking surface to facilitate modeling of the rocking response of deformable bodies and structures. The augmented model is used to compute the overturning and uplift rocking response spectra for a deformable rocking frame structure to symmetric and anti‐symmetric Ricker pulse ground motion excitation. It is found that the deformability of the columns of a rocking frame does not jeopardize its stability under Ricker pulse ground motion excitation. In fact, there are cases where a deformable rocking frame is more stable than its rigid counterpart. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluating a process‐based model for use in streambank stabilization: insights on the Bank Stability and Toe Erosion Model (BSTEM)

Streambank retreat is a complex cyclical process involving subaerial processes, fluvial erosion, seepage erosion, and geotechnical failures and is driven by several soil properties that themselves are temporally and spatially variable. Therefore, it can be extremely challenging to predict and model the erosion and consequent retreat of streambanks. However, modeling streambank retreat has many important applications, including the design and assessment of mitigation strategies for stream revitalization and stabilization. In order to highlight the current complexities of modeling streambank retreat and to suggest future research areas, this paper reviewed one of the most comprehensive streambank retreat models available, the Bank Stability and Toe Erosion Model (BSTEM), which has recently been integrated with several popular hydrodynamic and sediment transport models including the Hydrologic Engineering Center's River Analysis System (HEC‐RAS). The objectives of this paper were to: (i) comprehensively review studies that have utilized BSTEM and report their findings, (ii) address the limitations of the model so that it can be applied appropriately in its current form, and (iii) suggest directions of research that will help make the model a more useful tool in future applications. The paper includes an extensive overview of peer reviewed studies to guide future users of BSTEM. The review demonstrated that the model needs further testing and evaluation outside of the central United States. Also, further development is needed in terms of accounting for spatial and temporal variability in geotechnical and fluvial erodibility parameters, incorporating subaerial processes, and accounting for the influence of riparian vegetation on streambank pore‐water pressure dynamics, applied shear stress, and erodibility parameters. Copyright © 2016 John Wiley & Sons, Ltd.     

柴达木盆地东部地震地面运动放大效应

柴达木盆地是青藏高原东北部大型断陷山间盆地,该地区的流动观测记录了2008年11月10日发生于大柴旦附近的M_W6.3地震。和附近的基岩上的记录相比,盆地内部的记录显示出非常显著的地面运动放大效应,表现为峰值速度的增大、持续时间的延长,其呈现出长持续时间的后续震相。傅里叶频谱分析表明盆地内部显著的后续震相的频率和直达波相比较低,地面质点运动轨迹图显示后续震相为面波运动特征。为了解释地面运动的差异,构建二维模型,通过交错网格高阶有限差分方法计算了地震波在盆地内部的传播过程,结果显示盆地内部低速层的存在造成直达波的放大以及多次反射与转换,盆地边缘结构造成的波的相干叠加产生了强烈的次生面波,其低频、大振幅、长持续时间的特征是盆地内部地面运动放大的主要原因。     

土壤水环境中有机污染物生物处理过程模型研究

生物处理是一种经济有效处理土壤水环境中有机污染物的手段,本文在研究土地生物处理过程的基础上,建立了综合描述有机污染物在土壤－水－微生物系统中扩散、吸附／解吸、屏蔽和生物降解过程的数学模型。为确定模型中各参数在模型计算中的作用和相对重要性,进行了参数灵敏度分析,预计数学模型可以定量预测有机污染物进行土地生物处理所需的要时间和程度,为构建土地生物处理工程提供参考。     

Recent increase of river–floodplain suspended sediment exchange in a reach of the lower Amazon River

We analyzed variation of channel–floodplain suspended sediment exchange along a 140 km reach of the lower Amazon River for two decades (1995–2014). Daily sediment fluxes were determined by combining measured and estimated surface sediment concentrations with river–floodplain water exchanges computed with a two‐dimensional hydraulic model. The average annual inflow to the floodplain was 4088 ± 2017 Gg yr^?1 and the outflow was 2251 ± 471 Gg yr^?1, respectively. Prediction of average sediment accretion rate was twice the estimate from a previous study of this same reach and more than an order of magnitude lower than an estimate from an earlier regional scale study. The amount of water routed through the floodplain, which is sensitive to levee topography and increases exponentially with river discharge, was the main factor controlling the variation in total annual sediment inflow. Besides floodplain routing, the total annual sediment export depended on the increase in sediment concentration in lakes during floodplain drainage. The recent increasing amplitude of the Amazon River annual flood over two decades has caused a substantial shift in water and sediment river–floodplain exchanges. In the second decade (2005–2014), as the frequency of extreme floods increased, annual sediment inflow increased by 81% and net storage increased by 317% in relation to the previous decade (1995–2004). Copyright © 2017 John Wiley & Sons, Ltd.     

A one‐dimensional process‐based approach to study reservoir sediment dynamics during management operations

A comprehensive understanding of the dynamics of erosion and sedimentation in reservoirs under different management conditions is required to anticipate sedimentation issues and implement effective sediment management strategies. This paper describes a unique approach combining fluvial geomorphology tools and morphodynamic modeling for analyzing the sediment dynamics of an elongated hydropower reservoir subjected to management operations: the Génissiat Reservoir on the Rhône River. Functional sub‐reaches representative of the reservoir morphodynamics were delineated by adapting natural river segmentation methods to elongated reservoirs. The segmentation revealed the link between the spatial and temporal reservoir changes and the variability of longitudinal flow conditions during reservoir management operations. An innovative modeling strategy, incorporating the reservoir segmentation into two sediment transport codes, was implemented to simulate the dynamics of erosion and sedimentation at the reach scale during historic events. One code used a bedload approach, based on the Exner equation with a transport capacity formula, and the other used a suspended load approach based on the advection–dispersion equation. This strategy provided a fair quantification of the dynamics of erosion and sedimentation at the reach scale during different management operations. This study showed that the reservoir morphodynamics is controlled by bedload transport in upper reaches, graded suspended load transport of sand in middle reaches and suspended load transport of fine sediments in lower reaches. Eventually, it allowed a better understanding of the impact of dam management on sediment dynamics. Copyright © 2017 John Wiley & Sons, Ltd.     

应急演练中基于G-R关系的主余型地震余震震级模拟

在现有的日常地震演练过程中,与应急救援密切相关的地震余震信息产品较为缺乏,直接影响发震构造的判断以及影响场修正等关键环节。本文从计算机系统提供的均匀分布随机数出发,运用反函数法模拟生成余震序列,并进行系统检验,证实该方法产生的余震序列满足G-R频次关系。模拟生成的余震震级数据既有助于增强地震应急救援演练的现实性,也有助于丰富地震应急宣传产品,提升地震部门的履职能力。     

托素湖沉积物中重金属元素分析及地球化学特征研究

首次对托素湖沉积岩芯上部沉积物的重金属元素Cr、Cd、Cu、Zn和Pb的含量进行分析测试。精确测得沉积物中Cr,Cd,Cu,Zn和Pb的平均浓度分别为32.9 mg.kg-1,2.2 mg.kg-1,19.5 mg.kg-1,188.5 mg.kg-1和32.3 mg.kg-1。与区域背景值和标准值相比,托素湖沉积物样品中Cd和Zn都偏高。研究结果表明除了Cd污染程度比较严重,托素湖当前的Cr、Cu、Zn和Pb污染程度相对较轻。各元素经Ti校正后的分析结果显示,重金属元素含量在沉积剖面的10cm处出现波动性增大的明显转折,表明托素湖沉积物中Cd、Cr、Cu、Zn和Pb主要来源于流域侵蚀的造岩矿物元素,Cd污染可能源于人类活动所致。 关键词：托素湖;重金属元素;污染程度;地球化学特征