Automatic sampling of stream water during storm events in small remote catchments

In order to study the relationship between water composition and stream flow rate, it is desirable to sample at a frequency related to flow rate, especially during storm events. In a rural catchment of 18 ha near Oxford, the rate of rainfall was found to be linearly related to discharge on the rising limb of the stream hydrograph. A sampling system was therefore designed in which electrical pulses from a tipping-bucket raingauge were used to initiate and control the action of an automatic water sampler. A threshold rainfall intensity is set above which sampling commences. Sampling then continues at regular increments of rainfall until the intensity drops below the threshold, after which sampling occurs at regular intervals during the period that the stream flow reverts to normal. The CMOS electrical circuits which control the sampling also operate a cassette tape recorder which records the time of each tip of the raingauge and operation of the sampler. Since the sytem is designed to impose very little additional load on the battery which powers the water sampler, and can operate unattended for at least a fortnight, it is ideal for use in small, remote catchments. The system has been extended to include measurements of water temperature and could provide other measurements as well.     

Small channels: Their experimental formation,nature, and significance

The manner in which small channels are generated, from plane beds beneath sheet flows, has been experimentally elucidated. On plane, erodible, sand beds, the transition from thin, supercritical sheet flows to the channelled condition was studied over ranges of discharge, slope, and temperature. Secondary flow of the second kind, its action facilitated by steep vertical velocity gradients in the primary flows, caused sheet-flow instability. Along junctions between neighbouring secondary cells, both either raised or lowered elements of the primary flow. In the latter case, fast surface water was lowered to the bed, causing relatively intense, local, longitudinal scour. Dislodged grains were moved divergently to either side, leaving straight, central trenches. Development of positive feedback between cells and trenches led to rapid enlargement of the latter and concomitant growth of paired levees. The resulting structures, ‘protochannels’, were themselves ephemeral, developing two types of instability associated with secondary flow of the first kind. Firstly, small deviations from bilateral symmetry were enhanced, causing evolution into meandering channels. Secondly, headcutting led to multiple tributary development and, at resulting confluences, the action of strong pairs of secondary cells led to the development of braiding channels. Because they are shortlived, protochannels are but rarely seen in nature. Their seeding is markedly temperature-sensitive, reflecting their frictional origin. The erosive power of shallow overland flow largely depends on flow-energy concentration by secondary flow, firstly into channels, then within the channels themselves. Suppression of secondary flow, as by intense raindrop bombardment, can stabilize sheet flows. In deeper water, the effects of secondary flow appear relatively less dramatic. However, even if such motion is weak, bedload divergence from attachment lines can favour entrainment locally and thus affect bed geometry. Analogy between our results and river behaviour appears close and. on continental shelves where water must often flow as sheets, structures resembling giant protochannels evidently persist.     

地下水开采—地面沉降模型研究

本文建立三维流支－－一维非线性固结地面沉降模型,在混合井流,降雨滞后补给,初始水头形成,人为边界刻画,水流－固结耦合及软土层固结滞后于地下水开采层水头变化等方面有所改进,所建模型用于苏州市,模拟出若干地面沉降重要特征,地面沉降中心偏离地下水漏斗中心,地面沉降动态滞后于地下水的水头动态,软土层渗透系数随固结过程的变化及头在软土层中的传递特征。     

Hysteresis‐based flood‐wave analysis using the concept of strain

Hysteresis represents a loop in a rating curve and is a phenomenon which closely resembles that occurring in stress–strain curves used for studying the elastic properties of solid substances in engineering mechanics. Earlier hysteresis‐based studies used for defining floodwave propagation in open channels have qualitatively shown that hysteresis is an index of energy loss during floodwave propagation. Using the concept of elasticity, this paper introduces a new term called flow strain (defined as the ratio of change in discharge to the initial discharge) for investigating hysteresis. The usefulness of this new term is evaluated with use of four dam‐break studies. The study reveals that:

• 1 flow strain is a function of three wave speeds, Seddon speed, Lagrange speed, and elastic speed;
• 2 a single linear reservoir concept frequently used in flood routing is a specific variant of the Seddon speed formula;
• 3 the non‐linear storage–discharge relationship, widely used in overland flow modelling, is a variant of the kinematic wave representation;
• 4 the discharge ordinates on the recession part of a hydrograph follow a simple first‐order autoregressive form;
• 5 the hysteresis, phase difference and logarithmic decrement all define attenuation and are indices of energy loss during floodwave propagation.

Copyright © 2001 John Wiley & Sons, Ltd.     

Risk Management for Petroleum Reservoir Production: A Simulation-Based Study of Prediction

We consider numerical solutions of the Darcy and Buckley–Leverett equations for flow in porous media. These solutions depend on a realization of a random field that describes the reservoir permeability. The main content of this paper is to formulate and analyze a probability model for the numerical coarse grid solution error. We explore the extent to which the coarse grid oil production rate is sufficient to predict future oil production rates. We find that very early oil production data is sufficient to reduce the prediction error in oil production by about 30%, relative to the prior probability prediction.     

图们江流域中更新世火山泥石流分布及其特征

图们江流域中更新世火山泥石流仅分布在残留的Ⅲ阶地上，碎屑物的粒级范围很广，具正粒序层理，分选差，是在冰川作用或冰雪融化作用下快速搬运、沉积形成的，具有浊积岩的沉积模式。     

直线补给边界非稳定井流数值模拟研究

针对解析方法求解半无限含水层井流问题,对直线补给边界附近的井流采用数值方法进行模拟计算,并将两种方法进行比较,数值计算结果表明,存在补给边界时,抽水一定时间以后降深能达到稳定,但仅仅代表观测点处降深达到稳定,而整个渗流场并未达到稳定状态.     

SEAWAT 2000: modelling unstable flow and sensitivity to discretization levels and numerical schemes

A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60% of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.     

Post-drilling destabilization of temperature profile in borehole Yaxcopoil-1, Mexico

As part of the International Continental Scientific Drilling Program (ICDP), the 1.5-km-deep borehole Yaxcopoil-1, located in the Chixculub meteor impact structure in Mexico, has undergone further study after drilling operations ceased. Temperature logs were repeated ten times at intervals 0.3–0.8, 15, 24 and 34 months after borehole shut-in. The logs bear a distinct signature of transient heat transfer by groundwater flow manifested by a gradual distortion of the linear temperature profile when a cold wave of 0.8–1.6°C amplitude was detected propagating downward from 145 to 312 m at a rate of 4–6 m/month. To understand the nature of this moving anomaly, a 20-day monitoring of the cold wave was carried out at a depth of 307 m that showed further cooling of 0.6°C during the first 16 days of the passage followed by temperature stabilisation. As an explanation of this unusual phenomenon, a theory is proposed, whereby the drilling mud has accumulated within the overlying and cooler highly porous and permeable karstic rocks during the drilling and migrates downward. The observed migration rate suggests a permeability higher than 10^?11 m². This indicates a high vulnerability to contamination of the only freshwater aquifer in the Yucatan region.