基于聚类有效性函数的面状地理实体聚类

为解决聚类数未知条件下面状地理实体的聚类问题,文中提出了一种基于聚类有效性函数的聚类方法。给出了适合面状地理实体k-中心点聚类算法的聚类有效性函数;将该有效性函数改写为适应度函数,设计了基于遗传算法的面状地理实体聚类算法。该算法在计算聚类数的同时能得到划分聚类结果。实验结果从一定程度上反映了数据集的结构信息特征。     

Landform‐oriented flow‐routing algorithm for the dual‐structure loess terrain based on digital elevation models

The loess landform in the Loess Plateau of China is with typical dual structure, namely, the upper smooth positive terrain and the lower cliffy negative terrain (P–N terrain for short). Obvious differences in their morphological feature, geomorphological mechanism, and hydrological process could be found in the both areas. Based on the differences, a flow‐routing algorithm that separately addresses the dual‐structure terrain would be necessary to encompass this spatial variation in their hydrological behaviour. This paper proposes a mixed flow‐routing algorithm to address aforementioned problems. First, the loess landform surface is divided into P–N terrains based on digital elevation models. Then, specific catchment area is calculated with the new algorithm to simulate the water flows in both positive and negative terrain areas. The mixed algorithm consists of the multiple flow‐routing algorithm (multiple‐flow direction) for positive areas and the D8 algorithm for negative areas, respectively. The approach is validated in two typical geomorphologic areas with low hills and dense gullies in the northern Shaanxi Loess Plateau. Four indices are used to examine the results, which show that the new algorithm is more suitable for loess terrain in simulating the spatial distribution of water accumulation, as well as in modeling the flow characteristics of the true surface by considering the morphological structures of the terrain. Copyright © 2013 John Wiley & Sons, Ltd.     

Scalar advection schemes for ocean modelling on unstructured triangular grids

Several schemes for scalar advection on unstructured triangular grids are assessed for possible use in ocean modelling applications. Finite element, finite volume and finite volume–element approaches are evaluated. A series of tests, including a numerical order of convergence analysis, idealized rotating cone and cylinder experiments, and transport of a tracer through the Stommel Gyre representation of ocean basin-scale circulation, are carried out. Volume element Eulerian–Lagrangian and third-order Runge-Kutta discontinuous Galerkin schemes are recommended for use in tracer studies. Taylor–Galerkin and second-order Runge–Kutta discontinuous Galerkin are found to be robust and accurate second-order schemes. When positivity is required, a fluctuation redistribution scheme was found to be an easily implemented, accurate, and computationally efficient approach. Responsible editor: Phil Dyke     

A parallel evolutionary strategy based simulation–optimization approach for solving groundwater source identification problems

Groundwater characterization involves the resolution of unknown system characteristics from observation data, and is often classified as an inverse problem. Inverse problems are difficult to solve due to natural ill-posedness and computational intractability. Here we adopt the use of a simulation–optimization approach that couples a numerical pollutant-transport simulation model with evolutionary search algorithms for solution of the inverse problem. In this approach, the numerical transport model is solved iteratively during the evolutionary search. This process can be computationally intensive since several hundreds to thousands of forward model evaluations are typically required for solution. Given the potential computational intractability of such a simulation–optimization approach, parallel computation is employed to ease and enable the solution of such problems. In this paper, several variations of a groundwater source identification problem is examined in terms of solution quality and computational performance. The computational experiments were performed on the TeraGrid cluster available at the National Center for Supercomputing Applications. The results demonstrate the performance of the parallel simulation–optimization approach in terms of solution quality and computational performance.     

关于-遗传算法收敛性的注记

遗传算法是一种受到广泛注意的全局优化算法,已经在包括地震工程的很多领域中获得应用.本文将结合这一算法的实际操作步骤,对简单遗传算法的不收敛性和包含最优个体保护策略的遗传算法的收敛性给出一个简要的证明.     

3D Gravity Inversion by Growing Bodies and Shaping Layers at Mt. Vesuvius (Southern Italy)

To improve our knowledge of the structural pattern of Mt. Vesuvius and its magmatic system, which represents one of the three volcanoes located in the Neapolitan area (together with Campi Flegrei and Ischia; southern Italy), we analyze here the Bouguer gravity map that is already available through its interpretation by means of 2.5-dimensional modelling. We have carried out a three-dimensional interpretation using a new and original algorithm, known as ‘Layers’, that has been especially processed for this purpose. Layers works in an automatic and non-subjective way, and allows the definition of the structural settings in terms of several layers, each representing a specific geological formation. The same data are also interpreted in terms of isolated and shallow anomalous density bodies using a well tested algorithm known as ‘Growth’. We focus our inversions on the Mt. Vesuvius volcano, while globally analyzing the entire Neapolitan area, in order to investigate the deep structures, and in particular the deep extended ‘sill’ that has been revealed by seismic tomography. The final models generally confirm the global setting of the area as outlined by previous investigations, mainly for the shape and depth of the carbonate basement below Mt. Vesuvius. The presence of lateral density contrasts inside the volcano edifice is also shown, which was only hypothesized in the 2.5-dimensional inversion. Moreover, the models allow us to note a high density body that rises from the top of the carbonate basement and further elongates above sea level. This probably represents an uprising of the same basement, which is just below the volcano and which coincides with the V_P and V_P/V_S anomalies detected under the crater. The three-dimensional results also reveal that the two inversion methods provide very similar models, where the high density isolated body in the Growth model can be associated with the rising high density anomaly in the Layers model. Taking into account the density of these modelled bodies, we would also suggest that they represent solidified magma bodies, as suggested by other studies. Finally, we did not clearly detect any deep anomalous body that can be associated with the sill that was suggested by seismic tomography.     

Are modern metaheuristics successful in calibrating simple conceptual rainfall–runoff models?

In recent years sampling approaches have been used more widely than optimization algorithms to find parameters of conceptual rainfall–runoff models, but the difficulty of calibration of such models remains in dispute. The problem of finding a set of optimal parameters for conceptual rainfall–runoff models is interpreted differently in various studies, ranging from simple to relatively complex and difficult. In many papers, it is claimed that novel calibration approaches, so-called metaheuristics, outperform the older ones when applied to this task, but contradictory opinions are also plentiful. The present study aims at calibration of two simple lumped conceptual hydrological models, HBV and GR4J, by means of a large number of metaheuristic algorithms. The tests are performed on four catchments located in regions with relatively similar climatic conditions, but on different continents. The comparison shows that, although parameters found may somehow differ, the performance criteria achieved with simple lumped models calibrated by various metaheuristics are very similar and differences are insignificant from the hydrological point of view. However, occasionally some algorithms find slightly better solutions than those found by the vast majority of methods. This means that the problem of calibration of simple lumped HBV or GR4J models may be deceptive from the optimization perspective, as the vast majority of algorithms that follow a common evolutionary principle of survival of the fittest lead to sub-optimal solutions.     

Evolutionary aseismic design and retrofit of structures with passive energy dissipation

A new computational framework is developed for the design and retrofit of building structures by considering aseismic design as a complex adaptive process. For the initial phase of the development within this framework, genetic algorithms are employed for the discrete optimization of passively damped structural systems. The passive elements may include metallic plate dampers, viscous fluid dampers and viscoelastic solid dampers. The primary objective is to determine robust designs, including both the non‐linearity of the structural system and the uncertainty of the seismic environment. Within the present paper, this computational design approach is applied to a series of model problems, involving sizing and placement of passive dampers for energy dissipation. In order to facilitate our investigations and provide a baseline for further study, we introduce several simplifications for these initial examples. In particular, we employ deterministic lumped parameter structural models, memoryless fitness function definitions and hypothetical seismic environments. Despite these restrictions, some interesting results are obtained from the simulations and we are able to gain an understanding of the potential for the proposed evolutionary aseismic design methodology. Copyright © 2005 John Wiley & Sons, Ltd.     

Accelerating the convergence of coupled geomechanical‐reservoir simulations

The pressure variations during the production of petroleum reservoir induce stress changes in and around the reservoir. Such changes of the stress state can induce marked deformation of geological structures for stress sensitive reservoirs as chalk or unconsolidated sand reservoirs. The compaction of those reservoirs during depletion affects the pressure field and so the reservoir productivity. Therefore, the evaluation of the geomechanical effects requires to solve in a coupling way the geomechanical problem and the reservoir multiphase fluid flow problem. In this paper, we formulate the coupled geomechanical‐reservoir problem as a non‐linear fixed point problem and improve the resolution of the coupling problem by comparing in terms of robustness and convergence different algorithms. We study two accelerated algorithms which are much more robust and faster than the conventional staggered algorithm and we conclude that they should be used for the iterative resolution of coupled reservoir‐geomechanical problem. Copyright © 2006 John Wiley & Sons, Ltd.     

Quantifying flows along hydrological pathways by applying a new filtering algorithm in conjunction with master recession curve analysis

Quantifying the proportion of the river hydrograph derived from the different hydrological pathways is essential for understanding the behaviour of a catchment. This paper describes a new approach using the output from master recession curve analysis to inform a new algorithm based on the Lyne and Hollick ‘one‐parameter’ signal analysis filtering algorithm. This approach was applied to six catchments (including two subcatchments of these) in Ireland. The conceptual model for each catchment consists of four main flow pathways: overland flow, interflow, shallow groundwater and deep groundwater. The results were compared with those of the master recession curve analysis, a recharge coefficient approach developed in Ireland and the semi‐distributed, lumped and deterministic hydrological model Nedbør‐Afstrømings‐Model. The new algorithm removes the ‘free variable’ aspect that is typically associated with filtering algorithms and provides a means of estimating the contribution of each pathway that is consistent with the results of hydrograph separation in catchments that are dominated by quick response pathways. These types of catchments are underlain by poorly productive aquifers that are not capable of providing large baseflows in the river. Such aquifers underlie over 73% of Ireland, ensuring that this new algorithm is applicable in the majority of catchments in Ireland and potentially in those catchments internationally that are strongly influenced by the quick‐responding hydrological pathways. Copyright © 2013 John Wiley & Sons, Ltd.