Missing precipitation data estimation using optimal proximity metric-based imputation,nearest-neighbour classification and cluster-based interpolation methods

Abstract

New optimal proximity-based imputation, K-nearest neighbour (K-NN) classification and K-means clustering methods are proposed and developed for estimation of missing daily precipitation records. Mathematical programming formulations are developed to optimize the weighting, classification and clustering schemes used in these methods. Ten different binary and real-valued distance metrics are used as proximity measures. Two climatic regions, Kentucky and Florida, (temperate and tropical) in the USA, with different gauge density and network structure, are used as case studies to evaluate the new methods. A comprehensive exercise is undertaken to compare the performances of the new methods with those of several deterministic and stochastic spatial interpolation methods. The results from these comparisons indicate that the proposed methods performed better than existing methods. Use of optimal proximity metrics as weights, spatial clustering of observation sites and classification of precipitation data resulted in improvement of missing data estimates.

Editor D. Koutsoyiannis; Associate editor C. Onof     

Beerkan multi-runs for characterizing water infiltration and spatial variability of soil hydraulic properties across scales

A method is presented for characterizing the spatial variability of water infiltration and soil hydraulic properties at the transect and field scales. The method involves monitoring a set of 10 Beerkan runs distributed over a 1-m length of soil, and running BEST (Beerkan estimation of soil transfer parameters) methods to derive hydraulic parameters. The Beerkan multi-runs (BMR) method provides a significant amount of data at the transect scale, allowing the determination of correlations between water infiltration variables and hydraulic parameters, and the detection of specific runs affected by preferential flow or water repellence. The realization of several BMRs at several transects on the same site allows comparison of the variation between locations (spatial variability at the field scale) and at the transect scale (spatial variability at the metre scale), using analysis of variance. From the results, we determined the spatial variability of water infiltration and hydraulic parameters as well as its characteristic scale (transect versus field).     

基于居民地不同烈度下的地震灾害评估系统设计与实现

现有地震灾害评估系统多以设定地震的方式进行震害评估,无法满足地震灾害风险评估中按照不同烈度对整个行政区域进行地震灾害评估的业务需求,本文以现有计算模型为框架、以ArcEngine 为GIS开发平台、采用C/S结构进行基于居民地的不同烈度下地震灾害评估系统设计。系统以居民地为评估单元,根据选择的评估区域,按照Ⅵ度、Ⅶ度、Ⅷ度、Ⅸ度、大于等于Ⅹ度5个级别分别进行地震灾害建筑物损坏情况和人员伤亡评估,结果可为各级政府及相关部门进行区域地震风险评估、地震灾害预评估提供技术支持。     

火山岩风化壳地层型油气藏评价预测方法研究——以新疆北部石炭系为例

针对火山岩风化壳地层型油气藏强非均质性,评价预测难度大,勘探成功率低等难题,本文采用重磁电剥层处理、沿层延拓信号增强反演方法有效顶测区域火山岩分布;建立风化淋滤剥蚀后不同岩石组合的不完整火山机构和形态识别模式,利用相干体和振幅分析等方法有效识别火山岩目标;在建立单次火山喷发岩石序列及储层分布模式基础上,开发了基于次生溶...     

Statistical facies classification from multiple seismic attributes: comparison between Bayesian classification and expectation–maximization method and application in petrophysical inversion

We present here a comparison between two statistical methods for facies classifications: Bayesian classification and expectation–maximization method. The classification can be performed using multiple seismic attributes and can be extended from well logs to three‐dimensional volumes. In this work, we propose, for both methods, a sensitivity study to investigate the impact of the choice of seismic attributes used to condition the classification. In the second part, we integrate the facies classification in a Bayesian inversion setting for the estimation of continuous rock properties, such as porosity and lithological fractions, from the same set of seismic attributes. The advantage of the expectation–maximization method is that this algorithm does not require a training dataset, which is instead required in a traditional Bayesian classifier and still provides similar results. We show the application, comparison, and analysis of these methods in a real case study in the North Sea, where eight sedimentological facies have been defined. The facies classification is computed at the well location and compared with the sedimentological profile and then extended to the 3D reservoir model using up to 14 seismic attributes.     

设计反应谱标定方法的研究现状与讨论

设计反应谱及其标定方法的研究是工程抗震研究领域的基本问题之一。在现有研究成果的基础上,归纳和总结了设计反应谱及其标定方法的研究历史和现状;评述了研究进程中的若干节点问题;介绍了设计反应谱的标定原理和当前几种有代表性的设计反应谱的标定方法;分析了设计反应谱标定参数的影响因素。在此基础上,对设计反应谱及其标定方法研究中的强震资料积累问题、谱形状问题、标定参数的确定问题以及标定方法等问题进行了讨论并提出了改进和进一步研究的建议。     

Permanent earthquake‐induced actions in buried pipelines: Numerical modeling and experimental verification

Buried pipelines are often constructed in seismic and other geohazard areas, where severe ground deformations may induce severe strains in the pipeline. Calculation of those strains is essential for assessing pipeline integrity, and therefore, the development of efficient models accounting for soil‐pipe interaction is required. The present paper is aiming at developing efficient tools for calculating ground‐induced deformation on buried pipelines, often triggered by earthquake action, in the form of fault rupture, liquefaction‐induced lateral spreading, soil subsidence, or landslide. Soil‐pipe interaction is investigated by using advanced numerical tools, which employ solid elements for the soil, shell elements for the pipe, and account for soil‐pipe interaction, supported by large‐scale experiments. Soil‐pipe interaction in axial and transverse directions is evaluated first, using results from special‐purpose experiments and finite element simulations. The comparison between experimental and numerical results offers valuable information on key material parameters, necessary for accurate simulation of soil‐pipe interaction. Furthermore, reference is made to relevant provisions of design recommendations. Using the finite element models, calibrated from these experiments, pipeline performance at seismic‐fault crossings is analyzed, emphasizing on soil‐pipe interaction effects in the axial direction. The second part refers to full‐scale experiments, performed on a unique testing device. These experiments are modeled with the finite element tools to verify their efficiency in simulating soil‐pipe response under landslide or strike‐slip fault movement. The large‐scale experimental results compare very well with the numerical predictions, verifying the capability of the finite element models for accurate prediction of pipeline response under permanent earthquake‐induced ground deformations.     

Inverse reconstructions of ecosystem flows in investigating regime shifts: impact of the choice of objective function

Inverse analysis is increasingly used in ecosystem modelling to objectively reconstruct a large number of unknown flows or interactions from a small number of observations. This type of analysis may be useful in relating observed regime shifts in ecosystem structure to underlying processes. Inversions of ecosystem flow networks currently use a constrained least-squares solution which at the same time minimizes the squared norm (the sum of squares) of the reconstructed flows. This minimum norm (MN) inversion is thought to be a parsimonious solution to the ecosystem flow inverse problem, but it may well not reflect how ecosystems are organised. It has been proposed instead that ecosystems evolve to maximize energy/mass flows or that they maximize the information content of the network weighted by ecosystem flows (ascendancy). We used simulated inverse experiments, where inverse analyses are applied to simulations of flow networks, to explore objective functions different than the MN generally used. We could not compute inverse solutions that maximize ascendancy because the objective function is unbounded. We could calculate inversions that maximize flows; however, these generally overestimated the simulated flows, even though the simulations were designed to maximize flows. It appears that the ecosystem flow inverse problem is too under-determined (too few data relative to the number of unknowns) to allow the use of these maximizing goal functions. We introduce a new minimization that simultaneously minimizes the squared flows and the squared differences between flows. This smoothing minimization makes the inverse flows as even as possible and it helps with some technical issues with MN inversions. The simulated inverse experiments indicated that this smoothed norm (SM) is the most robust in comparative analyses of contrasting ecosystem states, such as those that can be associated with regime shifts. Like the MN inversion, the SM inversion has no ecological basis. However, it is a conservative norm that is less likely to produce false differences between the dynamics of regimes.