Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data     

Zhiwei Xu  James Irving  Kyle Lindsay  John Bradford  Peimin Zhu  Klaus Holliger 《Geophysical Prospecting》2020,68(2):678-689

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.  相似文献   

与高岭石固相平衡的酸性天然地表水中铝形态分布的计算机模拟     

刘锋  毕树平  汤伟 《湖泊科学》2000,12(4):311-320

基于热力学平衡计算,用计算机模拟研究了与高岭石固相平衡的酸性天然地表水中铝的形态分布,讨论了影响铝莆态分布的各种因素,同三水铝石固相平衡相比,与高岭石固相平衡的酸性天然地表水中,铝的形态分布与溶解硅酸的浓度密度相关,聚合态铝浓度较低,总溶解铝浓度明湿减小约２个数量级,模型被用于实际水样的分析,所得结果与实验值符合较好,最后用Ｍｏｎｔｅ Ｃａｒｌｏ模拟考察了由于平衡常数的选择而带来的不确定性。  相似文献   

与高岭石固相平衡的酸性天然地表水中铝的缓冲作用的计算机模拟     

刘锋  毕树平 《湖泊科学》2000,12(3):233-239

基于质量守恒以及电中性条件,采用化学平衡计算法,计算模拟研究了高岭石固相平衡的酸性天然地表水中铝的缓冲作用对ｐＨ值,酸中和容量ＡＮＣ和缓冲强度β的影响,讨论了影响酸性地表水ｐＨ估算的各种因素,如溶解硅酸的浓度,溶解总无机碳的浓度,溶解总氟和总有机碳的浓度及温度等。模型被用于实际水样的分析,所得结果与实验值符合较好,用ＭｏｎｔｅＣａｒｌｏ模拟考察了由于平衡常数的选择以及水化学参数测量误差而带来的不确  相似文献   

基于Monte Carlo模拟的三维剖面地质界线不确定性分析     

侯卫生  杨翘楚  杨亮  崔婵婕 《吉林大学学报(地球科学版)》2017,47(3)

三维剖面地质界线是构建三维地质结构模型的重要基础数据,其不确定性会影响三维模型的几何形态和属性分布。以单一分布为假设前提的统计学不确定性分析方法掩盖了其他概率分布特征对模型的影响。突破单一误差分布条件的假设前提,本文使用Monte Carlo方法模拟了不同概率分布情况下地质剖面数据中地质界线的抽样采集,以及地质界线空间分布的不确定性;依托地质界线空间位置与地质属性的耦合关系,提出了用地质属性概率分布实现地质界线空间不确定性的定量可视化,并结合实际地质剖面探讨了多种概率分布条件下地质界线的空间不确定性。实例研究表明,基于Monte Carlo模拟的不确定性分析方法可以突破单一误差分布假设条件,结合地质属性概率可充分揭示出建模数据的内在不确定性与模型外在要素形态之间的耦合关系。  相似文献   

裂隙岩体三维柔度张量的数值试验方法     

向前  何吉  陈胜宏  吴锦亮 《岩土力学》2015,36(12):3625-3632

在裂隙岩体二维柔度张量数值试验的基础上,建立了裂隙岩体三维柔度张量及其表征单元体积（REV）尺度研究的简化数值试验方法。首先根据岩体裂隙的统计参数及分布规律,运用蒙特卡罗方法在研究域内获得与实际岩体裂隙同分布的三维随机裂隙网络,提取位于不同方位的岩体试件,运用二维柔度张量的数值试验方法求出各个平面方向上的二维柔度张量,然后根据二维与三维柔度张量的拓扑关系,用数学方法求解裂隙岩体的三维柔度张量。对于含3组正交全贯通裂隙的岩体,通过数值试验获得的柔度张量与理论解相比,其主对角线上各参数的误差在5%以内,表明该数值试验方法较为可靠。最后将此方法运用到小湾水电站工程中,确定左岸坝区裂隙岩体的应力REV为11 m×11 m×11 m,并获得该区域裂隙岩体的三维柔度张量。  相似文献   

路堤边坡稳定可靠度计算中的模型不确定性分析     

吴兴正  蒋良潍  罗强  孔德惠  张良 《岩土力学》2015,36(Z2):665-672

基于均质路堤边坡Monte Carlo法的稳定可靠度计算,分析了临界滑面搜索策略和稳定分析方法两类模型不确定性对边坡可靠度的影响特性,讨论了边坡失效概率随土工参数变异性的变化规律。研究表明,选用不同的临界滑面搜索策略所得可靠度结果差异不大,参数滑面法（overall slope）的失效概率略大于均值滑面法（global minimum）,但差别对边坡稳定性分析没有实质性影响;土性参数变异水平是影响边坡可靠度的最重要因素,边坡在相同设计参数安全系数下的可靠度指标随参数变异性增大而急剧降低;不同稳定性分析方法对应的安全系数概率密度函数曲线形态基本一致,但失效概率差异明显,因此目标可靠度指标取值应与稳定性分析方法相适应。提出的考虑土工参数变异水平的安全系数取值修正原则,对改进确定性设计的边坡稳定分析技术有积极意义。  相似文献   

Stochastic dynamic stiffness of a surface footing for offshore wind turbines: Implementing a subset simulation method to estimate rare events     

《Soil Dynamics and Earthquake Engineering》2014

The purpose of this study, which concerns the stochastic dynamic stiffness of foundations for large offshore wind turbines, is to quantify uncertainties related to the first natural frequency of a turbine supported by a surface footing and to estimate the low event probabilities. Herein, a simple model of a wind turbine structure with equivalent coupled springs at the base is calibrated with the mean soil property values. A semianalytical solution, based on the Green׳s function for a layered half-space is utilized for estimation of foundation responses. Soil elastic modulus and layer depth are considered as random variables with lognormal distributions. The uncertainties are quantified, and the estimation of rare events of the first natural frequency is discussed through an advanced reliability approach based on subset simulation. This analysis represents a first step in the estimation of the safety with respect to the failure of a turbine in the fatigue limit state.  相似文献   

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

《Advances in water resources》2014

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.  相似文献   

塔克拉玛干沙漠周边地区极端弱降水的时空变化特征   总被引：3，自引：1，他引：2  

王新萍  杨青 《中国沙漠》2014,34(5):1376-1385

选用塔克拉玛干沙漠周边40个气象站1961-2009年日降水资料和4个极端弱降水指标,分析该地区极端弱降水的时空变化特征。采用M-K法和F检验对各站点降水指标的变化趋势及变化率进行检验和计算,并利用Monte Carlo模拟进行区域显著性检验。由 Copula函数得到两降水指标的联合分布,计算两降水指标的联合重现期。结果表明：（1）年最长连续无降水日数（CDD）多为80～100 d,呈显著减少趋势;出现频率最高的每年日降水量小于降水日序列25%分位数的日数（D₂₅）为0～10 d,呈显著增加趋势;每年日降水量小于降水日序列25%分位数的总降水量（P₂₅）的值集中在0～1.5 mm,2～3 mm的 P₂₅从2000年才开始出现;出现频率最高的每年日降水量小于降水日序列25%分位数的日平均降水量（I₂₅）为0.1～0.3 mm,I₂₅超过0.4 mm的情况极少出现。（2）CDD与D₂₅和P₂₅各自五年一遇值的空间分布相反。除CDD和D₂₅均大于各自五年一遇值的联合重现期较长外,其余各类型联合重现期较短。沙漠周边地区发生不同类型极端弱降水事件的概率不同。  相似文献   

Robustness of the holistic seismic risk evaluation in urban centers using the USRi   总被引：2，自引：2，他引：0  

Mabel-Cristina Marulanda  Omar D. Cardona  Alex H. Barbat 《Natural Hazards》2009,49(3):501-516

The Urban Seismic Risk index (USRi) published in a previous article (Carreño et al., Nat Hazards 40:137–172, 2007) is a composite indicator that measures risk from an integrated perspective and guides decision-making for identifying the main interdisciplinary factors of vulnerability to be reduced or intervened. The first step of the method is the evaluation of the potential physical damage (hard approach) as a result of the convolution of the seismic hazard with the physical vulnerability of buildings and infrastructure. Subsequently, a set of social context conditions that aggravate the physical effects is also considered (soft approach). According to this procedure, the physical risk index is evaluated for each unit of analysis from existing loss scenarios, whereas the total risk index is obtained by multiplying the former index by an impact factor using an aggravating coefficient, based on variables associated with the socio-economic conditions of each unit of analysis. The USRi has been developed using the underlying holistic and multi-hazard approach of the Urban Risk Index framework proposed for the evaluation of disaster risk in different megacities worldwide. This article presents the sensitivity analysis of the index to different parameters such as input data, weights and transformation functions used for the scaling or normalization of variables. This analysis has been performed using the Monte Carlo simulation to validate the robustness of this composite indicator, understanding as robustness how the cities maintain the ranking as well as predefined risk level ranges, when compared with the deterministic results of risk. Results are shown for different cities of the world.  相似文献