Correlation and covariance of runoff

The application of objective methods for interpolation of stochastic fields is based on the assumption of homogeneity with respect to the correlation function, i.e. only the relative distance between two points is of importance. This is not the case for runoff data which is demonstrated in this paper. Taking into consideration the structure of the river network and the related drainage basin supporting areas theoretical expressions are derived for the correlation function for flow along a river from its outlet and upstream. The results are exact for a rectangular drainage basin. For more complex basin geometry a grid approximation is suggested. The found relations are demonstrated on a real world example with a good agreement between the theoretically calculated correlation functions and empirical data.     

A generalization of the timoshenko beam model for coupled vibration analysis of thin-walled beams

A generalization of the Timoshenko beam model is presented which accounts for the influence of the shear strains, due to non-uniform bending and torsion, on the flexural-torsional vibrations of thin-walled cores with open or closed cross-section. The axial displacement field incorporates the torsion secondary warping as well as the warping terms depending on the shear resultants. It is shown that exact solutions for the interior domain problem can be obtained under proper load conditions. Moreover, a discrete model for the free-vibration analysis is derived by adopting a linear interpolation of the unknown functions and a reduced integration in order to avoid locking phenomena. Various applications are developed, including the case of the coupled vibrations of a shear-core-steel-frame building.     

The influence of the degree of saturation on dynamic response of a cylindrical lined cavity in a nearly saturated medium

In most previous studies on the dynamic response of a long cylindrical cavity subjected to internal transient dynamic loads, the porous medium was usually assumed to be completely saturated by ground water. In practice, however, the full saturation condition does not always exist. In this paper the surrounding soil and the lining of the cavity are respectively treated as a nearly saturated porous medium and an elastic material, and the governing equations for the dynamic problem are derived. A set of exact solutions are obtained in the Laplace transform domain for three types of transient loads, i.e. suddenly applied constant load, gradually applied step load and triangular pulse load. By utilizing a reliable numerical method of inverse Laplace transforms, the time-domain solutions are then presented. The influence of the degree of saturation of the surrounding soil on the dynamic response of the lined cavity is examined for numerical examples.     

成层粘弹性土中桩土耦合纵向振动时域响应研究

从三维轴对称角度出发,采用粘性阻尼粘弹性连续土介质模型,考虑桩土相互作用效应,对成层土中桩土纵向耦合振动时的桩顶时域响应进行了解析研究。求解时,首先建立定解问题,然后利用拉氏变换先对底部土层进行求解得到其振动位移形式解,然后利用桩土接触界面连续条件来考虑桩土耦合作用,分析底层土中桩段的动力反应,然后利用桩段阻抗函数的传递性,进行逐层递推求解,最终得到桩顶时域和频域响应的半解析解。通过参数影响分析和与工程实测曲线的对比,讨论分析了成层土中桩土耦合振动的响应特性,验证了本文解。基于本文研究可为桩基抗震、防震设计、桩基动力检测提供新的理论支持。     

Following the curve? Reviewing the physical basis of the SCS curve number method for estimating storm runoff

Less attention has been paid to runoff generation from semi-arid than from humid-temperate catchments. The SCS curve number approach is simple to apply and widely used, but lacks physical underpinning. Here output from a runoff generation models is compared with data from field measurements, making use of 11 years data from rainfall and runoff events at the Sierra de Enguera Soil Erosion Experimental Station in Eastern Spain. Runoff from natural rainfall events was monitored for 10 years on bare plots of 1–16 m length. The largest storm event was of 142 mm, generating runoff of up to 115 mm on the smallest plots. The model presented simulates overland storm flow on a sloping rough and unvegetated surface, representing an area of 320 × 320 m. Green-Ampt infiltration constants are randomly assigned to each cell in a 128 × 128 grid, and rectangular storms applied at a range of total amounts and intensities to simulate runoff at each transect across the area. A simple algebraic expression is developed to estimate total runoff and storage in terms of storm size and duration, and plot length, with parameters that reflect infiltration behaviour, and this expression is compared with the SCS curve number approach. For the very largest storms, both expressions converge asymptotically towards 100% runoff, but the revised expression greatly improves estimates of runoff from smaller events. Output of these simulations is compared with measured storm runoff data on bare runoff plots at the Sierra de Enguera experimental Station in SE Spain and gives further support to the proposed expression for storm runoff.     

Absorbing boundary conditions for dynamic analysis of fluid-saturated porous media

Based on Biot's two-phase mixture theory and the paraxial approximation, the absorbing boundary conditions in the time domain for u-w, u-U and u-p formulations are presented for the dynamic analysis of fluid-saturated porous media. These absorbing boundaries are equivalent to the viscous boundaries in the fundamental mode. The expressions for the energy ratio and reflection coefficient between the reflected and incident waves along the absorbing boundary are given. The numerical results show that the proposed absorbing boundaries can greatly suppress the spuriously reflected wave and model the far field of the system. These results also dynamic analysis of infinite fluid-saturated porous media. for the transient dynamic analysis of infinite fluid-saturated porous media.     

Research Note: A workaround for the corner problem in numerically exact non-reflecting boundary conditions

Simulations of wave propagation in the Earth usually require truncation of a larger domain to the region of interest to keep computational cost acceptable. This introduces artificial boundaries that should not generate reflected waves. Most existing boundary conditions are not able to completely suppress all the reflected energy, but suffice in practice except when modelling subtle events such as interbed multiples. Exact boundary conditions promise better performance but are usually formulated in terms of the governing wave equation and, after discretization, still may produce unwanted artefacts. Numerically exact non-reflecting boundary conditions are instead formulated in terms of the discretized wave equation. They have the property that the numerical solution computed on a given domain is the same as one on a domain enlarged to the extent that waves reflected from the boundary do not have the time to reach the original truncated domain. With a second- or higher-order finite-difference scheme for the one-dimensional wave equation, these boundary conditions follow from a recurrence relation. In its generalization to two or three dimensions, a recurrence relation was only found for a single non-reflecting boundary on one side of the domain or two of them at opposing ends. The other boundaries should then be zero Dirichlet or Neumann. If two non-reflecting boundaries meet at a corner, translation invariance is lost and a simple recurrence relation could not be found. Here, a workaround is presented that restores translation invariance by imposing classic, approximately non-reflecting boundary conditions on the other sides and numerically exact ones on the two opposing sides that otherwise would create the strongest reflected waves with the classic condition. The exact ones can also be applied independently. As a proof of principle, the method is applied to the two-dimensional acoustic wave equation, discretized on a rectangular domain with a second-order finite-difference scheme and first-order Enquist–Majda boundary conditions as approximate ones. The method is computationally costly but has the advantage that it can be reused on a sequence of problems as long as the time step and the sound speed values next to the boundary are kept fixed.     

Stochastic analysis of velocity spatial variability in bounded rectangular heterogeneous aquifers

Transport of inert solutes in two-dimensional bounded heterogeneous porous media is investigated in a stochastic framework. After adopting a first-order approximation of the flow equations, analytical expressions are derived for the velocity covariances. Effects of the boundary conditions and aquifer size upon the statistical moments are analyzed. While the size of the domain is shown to have small influence on the covariances in most cases, the solutions are considerably modified by the boundaries. The results are compared with analytical solutions on infinite domains, and several discrepancies are demonstrated. For example, while the velocity variances on infinite domains are homogeneous, the present results are strongly non-stationary. Finally, the problem is solved numerically by the Monte Carlo simulation method. The results, including the behavior near the boundaries, are shown to be in close agreement with analytical solutions.     

单程波算子积分解的象征表示

单程波波场延拓算子在地震偏移成像中有重要应用.单程波波场延拓算子按其实现方式可分为Kirchhoff积分、空间隐式有限差分和Fourier变换方法,他们代表了算子的不同表示方法,当截断使用这些方法时会得到不同的精度.象征表示对这些方法的导出和精度分析有重要作用.算子作用于正弦波函数所得函数称为算子的象征.算子的象征是褶积算子Fourier变换的推广.Fourier变换方法则直接用象征函数的可分表示求出.空间隐式有限差分则可以用象征函数的Padè近似或部分分式导出.单程波算子在深度域的积分称为单程波算子积分解.本文推导了单程波算子积分解的象征表达式,给出了算子象征的代数运算的头几阶表达式,这些表达式还未在前人文献中发现.Kirchhoff积分所需格林函数可以通过象征函数和鞍点法导出.基于积分解的象征表达式给出了非对称走时公式,对改善Kirchhoff积分的聚焦性能有重要意义.     

Direct simulation of solute recycling in irrigated areas

Solute recycling from irrigation can be described as the process that occurs when the salt load that is extracted from irrigation wells and distributed on the fields is returned to the groundwater below irrigated surfaces by deep percolation. Unless the salt load leaves the system by means of drains or surface runoff, transfer to the groundwater will take place, sooner or later. This can lead to solute accumulation and thus to groundwater degradation, particularly in areas where extraction rates exceed infiltration rates (semi-arid and arid regions). Thus, considerable errors can occur in a predictive solute mass budget if the recycling process is not accounted for in the calculation. A method is proposed which allows direct simulation of solute recycling. The transient solute response at an extraction well is shown to be a superposition of solute mass flux contributions from n recycling cycles and is described as a function of the travel time distribution between a recycling point and a well. This leads to an expression for a transient ‘recycling source’ term in the advection–dispersion equation, which generates the effect of solute recycling. At long times, the ‘recycling source’ is a function of the local capture probability of the irrigation well and the solute mass flux captured by the well from the boundaries. The predicted concentration distribution at steady state reflects the maximum spatial concentration distribution in response to solute recycling and can thus be considered as the solute recycling potential or vulnerability of the entire domain for a given hydraulic setting and exploitation scheme. Simulation of the solute recycling potential is computationally undemanding and can therefore, for instance, be used for optimisation purposes. Also, the proposed method allows transient simulation of solute recycling with any standard flow and transport code.     

Generalized failure envelope for caisson foundations in cohesive soil: Static and dynamic loading

The response of massive caisson foundations to combined vertical (N), horizontal (Q) and moment (M) loading is investigated parametrically by a series of three-dimensional finite element analyses. The study considers foundations in cohesive soil, with due consideration to the caisson-soil contact interface conditions. The ultimate limit states are presented by failure envelopes in dimensionless and normalized forms and the effects of the embedment ratio, vertical load and interface friction on the bearing capacity are studied in detail. Particular emphasis is given on the physical and geometrical interpretation of the kinematic mechanisms that accompany failure, with respect to the loading ratio M/Q. Exploiting the numerical results, analytical expressions are derived for the capacities under pure horizontal, moment and vertical loading, for certain conditions. For the case of fully bonded interface conditions, comparison is given with upper bound limit equilibrium solutions based on Brinch Hansen theory for the ultimate lateral soil reaction. A generalized closed-form expression for the failure envelope in M–Q–N space is then proposed and validated for all cases examined. It is shown that the incremental displacement vector of the caisson at failure follows an associated flow rule, with respect to the envelope, irrespective of: (a) the caisson geometry, and (b) the interface conditions. A simplified geometrical explanation and physical interpretation of the associativity in M-Q load space is also provided. Finally, the derived failure envelope is validated against low (0.67 Hz) and high frequency (5 Hz) dynamic loading tests and the role of radiation damping on the response of the caisson at near failure conditions is unraveled.     

One-dimensional magnetotelluric inversion using an adaptation of Zohdy's resistivity method

An iterative refinement method for determining a layered resistivity model from a Schlumberger or Wenner sounding curve is adapted to determine a layered resistivity model by using apparent resistivity and phase derived from the magnetotelluric impedance. Magnetotelluric observations presented as a function of period are first converted to an approximate resistivity–depth profile using Schmucker's transformation and this is used to construct an initial guess (starting) model. A two-stage procedure is then invoked. Keeping resistivities constant, layer boundaries are first adjusted to give a minimum misfit between measured data and responses and this is followed by resistivity adjustments with fixed layer boundaries to reduce the misfit further. The method is illustrated by application to some synthetic data both exact and with added noise, to a real field data set and to some magnetotelluric profile data obtained in a survey over the Carnmenellis granites in south Cornwall. The method is validated by recovering conductivity models from the exact and noisy 1D synthetic data. For complicated three-dimensional data at a single site and along a profile of stations, the method is shown to produce acceptable solutions which may be used as starting models in further two- or three-dimensional studies.     

地面导线环的辐射电阻

地面导线环是变频测深发射天线的一种形式.Wait曾计算了离地面高度为h的小导线环的辐射电阻,但由于所用近似方法的限制,他得出的结果不能用在置于地面的导线环的情况（若令h=0,所给出的辐射电阻将等于∞）.在本文中,我们导出了具有有限半径的地面导线环的辐射电阻表达式（准确的以及近似的）,并给出了数值计算的结果.     

Accurate and efficient shot-gather dip moveout processing in the log-stretch domain1

An accurate analytical expression for shot-gather dip-moveout (DMO) in the timespace log-stretch domain has until now not been published. We present a simpler, alternative derivation of the exact DMO relationships of Black et al. which correctly take account of the repositioning of the midpoint. A new computationally efficient frequency-wavenumber (F-K) DMO operator for shot profiles is then derived, based on these DMO relationships in the time-space log-stretch domain. The newly derived DMO operator is, unlike most other log-stretch DMO operators) accurate for the full range of reflector dips. Along with other schemes which are performed in the log-stretch domain, it offers considerable time savings over conventional DMO processing. We have compared numerically the impulse response of the new operator with those of a number of other shot-gather DMO operators, and found it to be superior and well match to the theoretical elliptical DMO response.     

A simple model of hillslope response for overland flow generation

This paper deals with the derivation of the hydrological response of a hillslope on the assumption of quick runoff by surface runoff generation. By using the simple non‐linear storage based model, first proposed by Horton, an analytical solution of the overland flow equations over a plane hillslope was derived. This solution establishes a generalization for different flow regimes of Horton's original solution, which is valid for the transitional flow regime only. The solution proposed was compared successfully with that of Horton and, for the turbulent flow regime, to the one derived from kinematic wave theory. This solution can be applied easily to both stationary and non‐stationary rainfall excess events. An analytical solution for the instantaneous response function (IRF) was also derived. Finally, simple expressions to compute peak and time to peak of IRF are proposed. Copyright © 2001 John Wiley & Sons, Ltd.     

基于改进点估计法的结构整体概率抗震能力分析

确定能力中位值和能力离差值是结构整体概率抗震能力分析的两个关键问题,文中分析了现有方法存在的缺点。在Zhao-Ono点估计法的基础上,引入基于随机向量边缘概率分布信息的Nataf变换,提出了改进的点估计法。将改进点估计法与Pushover分析相结合,提出了评估结构整体概率抗震能力统计矩的随机Pushover分析方法。以某五层三跨钢筋混凝土框架结构为例,应用本方法,进行结构整体概率抗震能力分析,得到了结构整体抗震能力的易损性曲线。分析表明,所提方法是一种具有较高效率和较好精度的结构整体概率抗震能力的分析方法。     

Reliable yield of reservoirs and possible effects of climatic change

Abstract

Statistics of seasonal and daily rainfall in nine regions covering the United Kingdom are presented, together with the best currently available predictions of changes in rainfall and evaporation, based on various general circulation model (GCM) assessments of the climate 40 years ahead. On that basis, runoff sequences are generated for the region of northwest England with north Wales (NWE & NW, humid temperate) and for southeast England (SEE, temperate but less humid), both for the current climate and for the year 2030 estimate, the latter incorporating higher evaporation losses. The future annual runoff is reckoned as 8% less than the current average in the SEE region and 4% less than that in the NWE & NW region. Direct supply reservoir yield versus storage graphs are then derived in each case for a range of duration-based probability of emptiness, their yields (i.e. sustained draw-off rates) and storage capacities being expressed as proportions of current annual runoff. The year 2030 graphs fall systematically below those for the current climate, representing 8% to 15% loss of yield from existing storage in the SEE region and 4% to 25% loss of yield in the NWE & NW region.     

SWAT-based runoff and sediment simulation in a small watershed,the loessial hilly-gullied region of China:capabilities and challenges

Model calibration and validation are necessary before applying it for scenario assessment and watershed management.This study presented the methodology of evaluating Soil and Water Assessment Tool(SWAT) and tested the feasibility of SWAT on runoff and sediment load simulation in the Zhifanggou watershed located in hilly-gullied region of China.Daily runoff and sediment event data from 1998-2008 were used in this study;data from 1998-2003 were used for calibration and 2004-2008 for validation.The evaluation statistics for the daily runoff simulation showed that the model results were acceptable,but the model underestimated the runoff for high-flow events.For sediment load simulation,the SWAT performed well in capturing the trend of sediment load,while the model tended to underestimate sediment load during both the calibration and validation periods. The disparity between observed and simulated data most likely resulted from limitations of the existing SCS-CN and MUSLE methods in the model.This study indicated that the modification of SWAT components is needed to take rainfall intensity and its duration into account to enhance the model performance on peak flow and sediment load simulation during heavy rainfall season.     

3D viscous-spring artificial boundary in time domain

After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.     

Mapping average annual runoff: a hierarchical approach applying a stochastic interpolation scheme

Abstract

A novel approach for mapping river runoff is presented. It is based on a disaggregation of the mean annual streamflow measured at the outlet of a basin to estimate water depths on elements of an exact partition of this basin. The developed technique is based on geostatistical interpolation procedures to which a global constraint of water balance has been added. The methodology is illustrated by a case study from a tributary to the Rhône River, France. The results were compared to an established method-the nested approach, and a cross-validation was performed for each mapping technique. The disaggregation approach appears to give the most consistent results. Finally, two gridded maps were derived by applying the disaggregation twice to assess water depth on an increasingly finer grid mesh. The global constraint of water balance was applied to each element of the coarser mesh to give estimates for the finer one.