Changes in hydrogeological properties of the River Choushui alluvial fan aquifer due to the 1999 Chi-Chi earthquake, Taiwan

Changes in hydrogeological properties of the River Choushui alluvial fan aquifer before and after the 1999 Chi-Chi earthquake, Taiwan, have been identified using pumping tests. Three wells, SH2, YL2 and SC2, located in a compressional zone with high coseismic groundwater levels, were tested. The threshold of the aquifer deformation with respect to transmissivity (T) is greater than that with respect to storage coefficient (S). Decreases in the post-earthquake S are approximately 60% at SH2 and SC2, indicating aquifer compression after the Chi-Chi earthquake. Changes in the post-earthquake T range from 61% increase to 0.8% decrease. Moreover, results from anisotropy analysis of T at SC2 further illustrate that normal stresses induced by the Chi-Chi earthquake have consolidated soil particles. Soil particles dilated laterally after the earthquake, resulting in an increase of the equivalent T. The changes in hydrogeological properties have a considerable influence on spatiotemporal fluid pressure and horizontal groundwater movement, resulting in different amounts of drawdown during post-earthquake pumping.     

DCT域遥感影像融合算法

提出一种基于离散余弦变换(Discrete Cosine Transform,DCT)与IHS变换的多光谱与全色遥感影像融合方法及其改进算法。本文方法依据DCT系数的特点在DCT域进行遥感影像融合,适合压缩格式的遥感影像快速融合。目视效果与客观评价表明,相比常用遥感影像融合方法,本文方法能在提高空间分辨率与保持光谱特性之间得到更好的折衷。     

The mass accommodation coefficient of ammonia on water

The mass accommodation coefficient of ammonia gas on water has been determined by measuring the absorption rate of 50–200 ppm NH₃ in one atm of air or helium into a liquid jet of 97 µm diameter as function of the exposed jet length, and comparing the results with numerical simulations which treat as the only free parameter. The model considers in detail transport of NH₃ by molecular diffusion, penetration of the gas/water interface, hydrolysis in the acidified water, and transport of the solutes from the surface into the jet. A correction is applied for the time evolution of the jet surface speed, using literature data on the fluid mechanics of liquid jets. The result of nine sets of independent measurements is

Trace gas exchange at the air/water interface: Measurements of mass accommodation coefficients

A liquid jet of 90 m diameter and variable length has been utilized to determine absorption rates and, hence, mass accommodation coefficients , of atmospheric trace gases. The compounds investigated are HCl (0.01), HNO₃ (0.01), N₂O₅ (0.005), peroxyacetyl nitrate (>0.001), and HONO (0.005). It is concluded that the absorption of these trace gases by liquid atmospheric water is not significantly retarded by interfacial mass transport. The strengths and limitations of the liquid jet technique for measuring mass accommodation coefficients are explored.     

Discrete numerical modelling of rockfill dams

The aim of this study is to obtain quantitative information on the behaviour of rockfill used in embankment dams, and particularly on the influence of block breakage on the displacement field, from a numerical analysis using the Distinct element method. A methodology is set up to define the resistance of the 2D particles so that the same probability of breaking blocks may be reproduced as in a 3D material. The model uses the discrete element code PFC^2D (Itasca Consulting Group Inc., PFC2D (Particle Flow Code in Two Dimensions), Version 3.0, 2002) and considers breakable clusters of 2D balls. The different parameters are determined from experimental data obtained from laboratory tests performed on rock blocks. The model is validated by comparing the results of the simulation of shearing tests with actual triaxial tests on rockfill material published in the literature. The numerical analysis of block crushing in an actual dam is proposed in the last part of this paper. Copyright © 2006 John Wiley & Sons, Ltd.     

A Connectivity Index for Discrete Fracture Networks

Connectivity is an important measure for assessing flow transport in rock, especially through fractures. In this paper, rock fracture systems are modelled by a discrete fracture model simulated by a marked point process. A connectivity index is then introduced to quantify the connectivity between any two points in space. Monte Carlo simulation is used to evaluate the connectivity index for stationary cases and relationships between the connectivity index and the parameters of the discrete fracture model are analysed. The average number of intersections per fracture, X_f, and the fracture intensity, P12 (P32), are calculated and the relationships between these parameters and the connectivity index are investigated, concluding that X_f is the more suitable parameter for the classification of rock mass flow properties. The relationships between the connectivity index and the percolation state of the fractured medium are also discussed. An edge correction is briefly discussed and a practical example is used to demonstrate the method of computing the connectivity index.     

Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling

Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.     

粗差发现和定位能力与相关系数的关系

根据粗差判断方程中的判断矩阵和两个统计检验量之间相关系数的函数式,论证了两种不同的研究方法所确定的粗差不能定位的数学模型实际上是相等的。通过算例,不仅说明两种研究方法对观测量不能定位粗差的判断是一致的,而且使用判断矩阵研究观测量的粗差发现和定位能力会更加方便简单。     

基于线性四叉树的全球离散格网索引

将四元三角剖分的三角形格网组织成四叉树剖分的菱形块格网,以菱形块单元来进行数据的组织索引,通过线性四叉树成熟的Morton编码作为关键字来标识查找菱形块,并在此基础上经过改进来索引三角形格网、进行邻域搜索。     

Laboratory Studies Of Wind Stress Over Surface Waves

Simultaneous laboratory observations of wind speed, wind stress, and surfacewind-wave spectra are made under a variety of wind forcing patterns using cleanwater as well as water containing an artificial surfactant. Under typical experimentalconditions, more than half of the total stress is supported by the wave-induced stressrather than by the surface viscous stress. When the surfactant reduces the shortwind-wave spectra, the wind stress also decreases by as much as 20–30% at agiven wind speed. When the wind forcing is modulated in time, the wind stresstends to be higher under decreasing wind than under increasing wind at a givenwind speed, mainly because the response of short wind-wave spectra to varyingwind forcing is delayed in time. These examples clearly demonstrate that therelationship between the wind speed and the wind stress can be significantlymodified if the surface wave field is not in equilibrium with the wind forcing.Next, we examine whether the wind stress is estimated accurately if the wave-inducedstress by all surface wave components is explicitly evaluated by linear superpositionand is added to the surface viscous stress. It is assumed that the surface viscous stressis uniquely related to the wind speed, and that the wind input rate is determined by thelocal, reduced turbulent stress rather than the total stress. Our wind stress estimatesincluding the wave contributions agree well with observed wind stress values, evenif the surface wave field is away from its equilibrium with the wind in the presenceof surface films and/or under time-transient wind forcing. These observations stronglysuggest that the wind stress is accurately evaluated as a sum of the wave-induced stressand the surface viscous stress. At very high winds, our stress estimates tend to be lowerthan the observations. We suspect that this is because of the enhancement of wind stressover very steep (or breaking) short wind-waves.