地下水模型的Neumann展开Monte-Carlo随机有限元法

讨论了求解地下水模型(水流模型和水质模型)的Neumann展开Monte-Carlo随机有限元法。从基本的随机变量入手,避免了过程中随机变量的增多问题,给出结点水头(浓度)的均值、方差和水头在某区间的概率计算方法;改进了矩阵求逆的效率,对输入随机变量较多、随机变量变异较大的非稳定地下水问题特别有效。同时选取二维承压地下水水流问题(有解析解)作为例子,进行了随机数值模拟实验。     

A review of methods used to determine the fractal dimension of linear features

An in-depth review of the more commonly applied methods used in the determination of the fractal dimension of one-dimensional curves is presented. Many often conflicting opinions about the different methods have been collected and are contrasted with each other. In addition, several little known but potentially useful techniques are also reviewed. General recommendations which should be considered whenever applying any method are made.     

利用爆破事件分析不同速度模型下常规定位方法定位结果

地震定位对速度模型的依赖性很强。四川地区地形复杂,常规工作中可选取多种速度模型进行定位。川西龙门山断裂带为东南部四川盆地和西北部青藏高原东部山区的明显分界线,近年在此断裂带上发生多次较大地震。对发生在该断裂带附近的6个爆破事件和15个天然地震重新定位,并对比结果。研究表明,相同台站包围情况下,川滇3D速度模型稳定性最好,但对浅表爆破不太准确。相比HypoSat(一维速度模型)组合,台站分布对Hypo2000(一维速度模型)和Hypo2000(赵珠速度模型)组合的定位结果影响较大。     

Discussion of ‘Modelling viscous damping in nonlinear response history analysis of buildings for earthquake excitation’ by Anil K. Chopra and Frank McKenna

This discussion deals with recommendations in the paper on appropriate damping formulations for use in nonlinear response history analysis of buildings. Concern over potentially excessive damping forces and moments should extend beyond the damping moments produced by the stiffness proportional part of Rayleigh damping that corresponds to rotational springs used to explicitly model plastic hinges. The key to an appropriate damping formulation for nonlinear analysis is a realistic mechanism that allows all damping forces and moments to be meaningfully assessed. Then features can be added to keep these forces and moments within reasonable bounds. Copyright © 2016 John Wiley & Sons, Ltd.     

Prediction of river water temperature: a comparison between a new family of hybrid models and statistical approaches

River water temperature is a key physical variable controlling several chemical, biological and ecological processes. Its reliable prediction is a main issue in many environmental applications, which however is hampered by data scarcity, when using data‐demanding deterministic models, and modelling limitations, when using simpler statistical models. In this work we test a suite of models belonging to air2stream family, which are characterized by a hybrid formulation that combines a physical derivation of the key equation with a stochastic calibration of parameters. The air2stream models rely solely on air temperature and streamflow, and are of similar complexity as standard statistical models. The performances of the different versions of air2stream in predicting river water temperature are compared with those of the most common statistical models typically used in the literature. To this aim, a dataset of 38 Swiss rivers is used, which includes rivers classified into four different categories according to their hydrological characteristics: low‐land natural rivers, lake outlets, snow‐fed rivers and regulated rivers. The results of the analysis provide practical indications regarding the type of model that is most suitable to simulate river water temperature across different time scales (from daily to seasonal) and for different hydrological regimes. A model intercomparison exercise suggests that the family of air2stream hybrid models generally outperforms statistical models, while cross‐validation conducted over a 30‐year period indicates that they can be suitably adopted for long‐term analyses. Copyright © 2016 John Wiley & Sons, Ltd.     

Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.     

Evaluation of remotely‐sensed DEMs and modification based on plausibility rules and initial sediment budgets of an artificially‐created catchment

To quantify landscape change resulting from processes of erosion and deposition and to establish spatially distributed sediment budgets, ‘models of change’ can be established from a time series of digital elevation models (DEMs). However, resolution effects and measurement errors in DEMs may propagate to these models. This study aimed to evaluate and to modify remotely‐sensed DEMs for an improved quantification of initial sediment mass changes in an artificially‐created catchment. DEMs were constructed from photogrammetry‐based, airborne (ALS) and ground‐based laser scanning (TLS) data. Regions of differing morphological characteristics and vegetation cover were delineated. Three‐dimensional (3D) models of volume change were established and mass change was derived from these models. DEMs were modified region‐by‐region for rill, interrill and alluvial areas, based on logical and hydro‐geomorphological principles. Additional DEMs were constructed by combining multi‐source, modified data. Models were evaluated by comparison with d‐GPS reference data and by considering sediment budget plausibility. Comprehensive evaluation showed that DEM usability depends on a relation between the technique used to obtain elevation data, surface morphology and vegetation cover characteristics. Photogrammetry‐based DEMs were suited to quantification of change in interrill areas but strongly underestimated surface lowering in erosion rills. TLS DEMs were best suited to rill areas, while ALS DEMs performed best in vegetation‐covered alluvial areas. Agreement with reference data and budget plausibility were improved by modifications to photogrammetry‐ and TLS‐based DEMs. Results suggest that artefacts in DEMs can be reduced and hydro‐geomorphic surface structures can be better represented by applying region‐specific modifications. Photogrammetry‐based DEMs can be improved by combining higher and lower resolution data in defined structural units and applying modifications based on principles given by characteristic hydro‐geomorphic evolution. Results of the critical comparative evaluation of remotely‐sensed elevation data can help to better interpret DEM‐based quantifications of earth‐surface processes. Copyright © 2012 John Wiley & Sons, Ltd.     

Seismic demand models for probabilistic risk analysis of near fault vertical ground motion effects on ordinary highway bridges

The influence of vertical ground motions on the seismic response of highway bridges is not very well understood. Recent studies suggest that vertical ground motions can substantially increase force and moment demands on bridge columns and girders and cannot be overlooked in seismic design of bridge structures. For an evaluation of vertical ground motion effects on the response of single‐bent two‐span highway bridges, a systematic study combining the critical engineering demand parameters (EDPs) and ground motion intensity measures (IMs) is required. Results of a parametric study examining a range of highway bridge configurations subjected to selected sets of horizontal and vertical ground motions are used to determine the structural parameters that are significantly amplified by the vertical excitations. The amplification in these parameters is modeled using simple equations that are functions of horizontal and vertical spectral accelerations at the corresponding horizontal and vertical fundamental periods of the bridge. This paper describes the derivation of seismic demand models developed for typical highway overcrossings by incorporating critical EDPs and combined effects of horizontal and vertical ground motion IMs depending on the type of the parameter and the period of the structure. These models may be used individually as risk‐based design tools to determine the probability of exceeding the critical levels of EDP for pre‐determined levels of ground shaking or may be included explicitly in probabilistic seismic risk assessments. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

IM‐based and EDP‐based decision models for the verification of the seismic collapse safety of buildings

Decision models for the verification of seismic collapse safety of buildings are introduced. The derivations are based on the concept of the acceptable (target) annual probability of collapse, whereas the decision making involves comparisons between seismic demand and capacity, which is familiar to engineering practitioners. Seismic demand, which corresponds to the design seismic action associated with a selected return period, can be expressed either in terms of an intensity measure (IM) or an engineering demand parameter (EDP). Seismic capacity, on the other hand, is defined by dividing the near‐collapse limit‐state IM or EDP by an appropriate risk‐targeted safety factor (γ _im or γ _edp ), which is the only safety factor used in the proposed decision model. Consequently, the seismic performance assessment of a building should be based on the best possible estimate. For a case study, it is shown that if the target collapse risk is set to 10^?4 (0.5% over a period of 50 years), and if the seismic demand corresponds to a return period of 475 years (10% over a period of 50 years), then it can be demonstrated that γ _im is approximately equal to 2.5 for very stiff buildings, whereas for buildings with long periods the value of γ _im can increase up to a value of approximately 5. The model using γ _edp is equal to that using γ _im only if it can be assumed that displacements, with consideration of nonlinear behavior, are equal to displacements from linear elastic analysis.