基于二维三次卷积插值算法的辛几何射线追踪

射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法（SAM-Symplectic Algorithm Method）及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.     

Multiridge Euler deconvolution

Potential field interpretation can be carried out using multiscale methods. This class of methods analyses a multiscale data set, which is built by upward continuation of the original data to a number of altitudes conveniently chosen. Euler deconvolution can be cast into this multiscale environment by analysing data along ridges of potential fields, e.g., at those points along lines across scales where the field or its horizontal or vertical derivative respectively is zero. Previous work has shown that Euler equations are notably simplified along any of these ridges. Since a given anomaly may generate one or more ridges we describe in this paper how Euler deconvolution may be used to jointly invert data along all of them, so performing a multiridge Euler deconvolution. The method enjoys the stable and high‐resolution properties of multiscale methods, due to the composite upward continuation/vertical differentiation filter used. Such a physically‐based field transformation can have a positive effect on reducing both high‐wavenumber noise and interference or regional field effects. Multiridge Euler deconvolution can also be applied to the modulus of an analytic signal, gravity/magnetic gradient tensor components or Hilbert transform components. The advantages of using multiridge Euler deconvolution compared to single ridge Euler deconvolution include improved solution clustering, increased number of solutions, improvement of accuracy of the results obtainable from some types of ridges and greater ease in the selection of ridges to invert. The multiscale approach is particularly well suited to deal with non‐ideal sources. In these cases, our strategy is to find the optimal combination of upward continuation altitude range and data differentiation order, such that the field could be sensed as approximately homogeneous and then characterized by a structural index close to an integer value. This allows us to estimate depths related to the top or the centre of the structure.     

Statistical models for shear strength of RC beam‐column joints using machine‐learning techniques

This paper proposes a new set of probabilistic joint shear strength models using the conventional multiple linear regression method, and advanced machine‐learning methods of multivariate adaptive regression splines (MARS) and symbolic regression (SR). In order to achieve high‐fidelity regression models with reduced model errors and bias, this study constructs extensive experimental databases for reinforced and unreinforced concrete joints by collecting existing beam‐column joint subassemblage tests from multiple sources. Various influential parameters that affect joint shear strength such as material properties, design parameters, and joint configuration are investigated through tests of statistical significance. After performing a set of regression analyses, the comparison of simulation results indicates that MARS approach is the best estimation method. Moreover, the accuracy of analytical predictions of the derived MARS model is compared with that of existing joint shear strength relationships. The comparison results show that the proposed model is more accurate compared to existing relationships. This joint shear strength prediction model can be readily implemented into joint response models for evaluation of earthquake performance and inelastic responses of building frames. Copyright © 2014 John Wiley & Sons, Ltd.     

中国西部地区利用烈度数据估计地震动参数的方法

通过最近收集到的中国西部地区的6次地震的地震动参数和烈度数据,建立了该地区适用于贝叶斯方法的各烈度档P(I|GM)的经验分布。利用2011年发生的四川炉霍地震和新疆伽师地震资料验证了依据烈度数据用贝叶斯方法估算峰值加速度的可行性。研究结果还表明,用贝叶斯方法估计的参数精度与选择的先验概率(衰减关系)显著相关。针对这2次地震,用各烈度档地震动参数的均值法、地震动衰减关系法及贝叶斯法3种方法估算的峰值加速度值与峰值加速度观测值之间的均方根比较表明,用贝叶斯方法估算的参数精度优于用另2种方法所获结果。     

Sampling efficiency in Monte Carlo based uncertainty propagation strategies: Application in seawater intrusion simulations

The implementation of Monte Carlo simulations (MCSs) for the propagation of uncertainty in real-world seawater intrusion (SWI) numerical models often becomes computationally prohibitive due to the large number of deterministic solves needed to achieve an acceptable level of accuracy. Previous studies have mostly relied on parallelization and grid computing to decrease the computational time of MCSs. However, another approach which has received less attention in the literature is to decrease the number of deterministic simulations by using more efficient sampling strategies. Sampling efficiency is a measure of the optimality of a sampling strategy. A more efficient sampling strategy requires fewer simulations and less computational time to reach a certain level of accuracy. The efficiency of a sampling strategy is highly related to its space-filling characteristics.This paper illustrates that the use of optimized Latin hypercube sampling (OLHS) strategies instead of the widely employed simple random sampling (SRS) and Latin hypercube sampling (LHS) strategies, can significantly improve sampling efficiency and hence decrease the simulation time of MCSs. Nine OLHS strategies are evaluated including: improved Latin hypercube sampling (IHS); optimum Latin hypercube (OLH) sampling; genetic optimum Latin hypercube (GOLH) sampling; three sampling strategies based on the enhanced stochastic evolutionary (ESE) algorithm namely φ_p-ESE which employs the φ_p space-filling criterion, CLD-ESE which utilizes the centered L₂-discrepancy (CLD) space-filling criterion, and SLD-ESE which uses the star L₂-discrepancy (SLD) space-filling criterion; and three sampling strategies based on the simulated annealing (SA) algorithm namely φ_p-SA which employs the φ_p criterion, CLD-SA which uses the CLD criterion, and SLD-SA which utilizes the SLD criterion. The study applies SRS, LHS and the nine OLHS strategies to MCSs of two synthetic test cases of SWI. The two test cases are the Henry problem and a two-dimensional radial representation of SWI in a circular island. The comparison demonstrates that the CLD-ESE strategy is the most efficient among the evaluated strategies. This paper also demonstrates how the space-filling characteristics of different OLHS designs change with variations in the input arguments of their optimization algorithms.     

云南地区近期地壳活动特征

基于云南省的GNSS连续观测资料,通过最小二乘、克里金插值等方法进行速度场、应变参数的求解,获得了自2010年起至今云南地区的地壳活动特征。结果表明,云南地区最大剪应变变化速率存在西强东弱,北强南弱的特征,目前滇西北地区最大剪应变增强最快,滇南至滇西南等区域挤压应变积累最明显,且2011年缅甸7.2级地震的发生未改变这一地区的应变积累趋势。     

Sediment settling in the Latian Dam in Iran

Among the difficulties that influence future dam operations,reservoir sedimentation is the most problematic for engineers.This study predicted the amount and pattern of sedimentation for use in estimation of the useful lifespan of reservoirs and identification of optimal locations for outlets and intakes at the initial stages of dam design.Hydrographic surveys of different dams can provide better insight into this phenomenon.Latian Dam in Iran has conducted hydrographic surveys during 7 time periods.The amount and process of sedimentation in this reservoir were determined,and predictions of distribution of sediments were validated by well-known,common methods.The formation of a delta in the reservoir was investigated for different time periods after operation.Future problems due to the impacts of sedimentation on dam operation and the useful lifespan of the reservoir were predicted.In addition,the study results may be used for developing empirical methods to predict sedimentation patterns in other reservoirs.     

泉州盆地地震效应特征的一维等效线性和二维非线性分析的比较

采用等效线性动粘弹性模型描述土的动力非线性特性,基于一维等效线性波传法,对泉州盆地地震效应进行了分析;同时,采用修正Martin-Seed-Davidenkov动粘弹塑性模型描述土的动力非线性特性,对泉州盆地非线性地震效应进行了大尺度二维精细化有限元分析,研究了地形地貌和土层横向不均匀性对地震效应的影响。将两种分析结果进行对比,结果表明：①随着基岩输入地震动强度增大,地表峰值加速度PGA放大效应总体呈现减小趋势,中震与小震、大震与小震的地表PGA放大系数之比依次为0.83~0.99、0.72~0.97;②该盆地Ⅲ类场地处,基岩、地表起伏不大,且土层横向分布较均匀,两种方法计算得到的地震效应特征类似;基岩或地表起伏剧烈、土层横向分布明显不均匀的Ⅱ类场地上,二维非线性分析给出的地表PGA放大系数明显大于一维等效线性结果,两种方法得到的地表加速度反应谱及PGA随土层深度的变化特征存在显著差异,二维非线性分析给出的地表加速度反应谱大多呈现双峰甚至多峰现象,且PGA在土层特定深度处存在聚集效应,使PGA随土层深度的变化呈现非单调性。