浅部频率域电磁勘探方法综述

适用于近地表(2000m以内)勘探的频率域电磁法主要有音频大地电磁法(audio-frequency magnetotellurics,AMT),无线电大地电磁法(radio-magnetotellurics,RMT),可控源音频大地电磁法(controlled source audiofrequency magnetotellurics,CSAMT),广域电磁法(Wide Field Electromagnetic Method,WFEM).本文拟从最新的数据采集技术、数据处理技术、正反演算法、实例等四个方面,论述适用于浅部勘探的AMT,RMT,CSAMT和WFEM方法的国内外最新进展,总结目前AMT,RMT,CSAMT和WFEM方法遇到的困难,对潜在的发展方向提出建议.综述表明:(1)张量测量、多站阵列、多站叠加可提高AMT、RMT和CSAMT数据的质量.利用近区数据WFEM法可获得良好的效果.国产与国外仪器在质量方面的差距正在逐步缩小.(2)数学形态滤波技术、Hilbert-Huang变换等可有效分离出有用的数据,局部畸变仍然是亟待解决的难题,需要更为深入的研究.(3)矢量有限元与非结构网格的出现大幅度提高了有限元处理复杂电磁问题模拟的精度与应用范围,成为目前电磁正演的首选工具.完全非线性反演算法仍然局限于1D、2D问题,共轭梯度法和高斯牛顿算法等为解决3D问题的发展趋势.地质约束的引入和多数据联合反演可以减小反演的非唯一性.各向异性的反演为目前反演研究的热点之一.(4)野外数据解释的正确性严重依赖于对地下结构先期的维性判别,在2D特性不明显、3D特性明显时,需要采用3D进行反演解释.     

基于合理权重赋值方法选择的多因素综合评价模型——以JCR中70种人文地理期刊为例

针对地理学多因素综合评价中较少探讨指标权重赋值方法的选择,以及少有研究涉及人文地理期刊的多因素综合评价,本文设计了基于合理权重赋值方法选择的多因素综合评价模型,并将其应用于科学网数据库内的“期刊引证报告（JCR）”中分类为人文地理类的70种人文地理期刊的综合评价。该模型针对7种指标权重客观赋值方法,采用“异同比较”,分析不同方法在不同参数设置下和聚类与否情况下的权重值、评价值和评价序的变化特征与分布规律,并依据该变化特征与分布规律,实现了对各种方法的“优劣比较”和最优选择。应用该模型对JCR中70种人文地理期刊进行排序分等,其中排名第一的“全球环境变化和政策”等14种期刊为位列第一等分的高影响力期刊。本研究完善了多因素综合评价体系,也丰富了对人文地理期刊的认知。     

基于交通信息流的江苏省流空间网络结构研究

基于江苏省13个地级市之间的公路、普速列车、高铁（CRH）与信息网络数据,从流空间视角对江苏省空间结构特征及其发展趋势进行综合分析与判断。结果表明：① 公路和普速铁路网络所展现的区域空间结构较为无序且影响范围有限,高铁网络则在大幅度拓展城市有效联系范围的同时凸显了发达地区的重要地位,而信息网络所反映的空间结构更加完整、有序与协调,代表了未来江苏省城市联系网络的发展方向;② 在高铁与信息网络的冲击下,流空间视角下的江苏省省域空间结构研究显出了一些区别于传统认知的新特征与变化趋势,如跳跃性的与非相邻地区产生主要联系、“等级+网络”的变化趋势等;③ 基于交通信息流的江苏省流空间结构呈现出明显的“中心-外围”格局,区域联系线密度由“沪宁线”向“沿海线”逐级递减,空间结构指数的计算结果也基本支持这一论断;④ 以高铁与信息为代表的新要素流对江苏省省域空间格局网络化的效率与公平产生了积极与消极两种影响,前者体现在改善区际可达性促进区域发展,后者体现在核心城市的“虹吸效应”。     

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.     

Application of diffracted wave analysis to time‐lapse seismic data for CO2 leakage detection

Time‐lapse seismic analysis is utilized in CO₂ geosequestration to verify the CO₂ containment within a reservoir. A major risk associated with geosequestration is a possible leakage of CO₂ from the storage formation into overlaying formations. To mitigate this risk, the deployment of carbon capture and storage projects requires fast and reliable detection of relatively small volumes of CO₂ outside the storage formation. To do this, it is necessary to predict typical seepage scenarios and improve subsurface seepage detection methods. In this work we present a technique for CO₂ monitoring based on the detection of diffracted waves in time‐lapse seismic data. In the case of CO₂ seepage, the migrating plume might form small secondary accumulations that would produce diffracted, rather than reflected waves. From time‐lapse data analysis, we are able to separate the diffracted waves from the predominant reflections in order to image the small CO₂ plumes. To explore possibilities to detect relatively small amounts of CO₂, we performed synthetic time‐lapse seismic modelling based on the Cooperative Research Centre for Greenhouse Gas Technologies (CO2CRC) Otway project data. The detection method is based on defining the CO₂ location by measuring the coherency of the signal along diffraction offset‐traveltime curves. The technique is applied to a time‐lapse stacked section using a stacking velocity to construct offset‐traveltime curves. Given the amount of noise found in the surface seismic data, the predicted minimum detectable amount of CO₂ is 1000–2000 tonnes. This method was also applied to real data obtained from a time‐lapse seismic physical model. The use of diffractions rather than reflections for monitoring small amounts of CO₂ can enhance the capability of subsurface monitoring in CO₂ geosequestration projects.     

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.     

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.