A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

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.     

地下洞库勘察中的深斜孔钻探和压水试验技术

青岛市黄岛区拟建设海底地下油气储藏库,前期地质勘察工作要求施工倾角45°、深220m的取心钻孔4个,并在孔内分段进行压水试验.本文介绍了如何利用XY-4型金刚石小口径钻机施工基岩深斜孔及压水试验技术.     

变形分析的Bayes估计与检验

由于观测条件的变化，观测值方差将产生随机波动。本文根据Ｂａｙｅｓ估计理论，推导了线性模型待估参数的分布函数及观测值方差估计公式，并将其用于位移检验。计算结果表明，在同样的置信度下，本方法更符合实际。文中还对参数估值的统计性质进行了讨论。     

塔中4号油藏渗透率分布的非均质性及渗流方向研究

研究油藏内部油水渗流方向时,主要考虑渗透率影响,引入两个水平方向渗透率比、水平与垂向渗透率比参数来分析。前者反映样品所在储层平面上两个水平方向渗流能力的差异,后者反映样品所在储层水平与垂向渗流能力的差异,并经小柱样测试和全岩心测试两种方法确定两个参数。通过对塔中4号油藏CⅠ油组9口井2306个样品两种渗透率比值分布直方图分析,研究了CⅠ油组渗透率的非均质性及其渗流方向:CI油组储层在水平方向的渗流能力有很大差异,具有明显的方向性。储层垂向渗流能力明显小于水平方向的渗流能力。     

Post-solidification cooling and the age of Io's lava flows

The modeling of thermal emission from active lava flows must account for the cooling of the lava after solidification. Models of lava cooling applied to data collected by the Galileo spacecraft have, until now, not taken this into consideration. This is a flaw as lava flows on Io are thought to be relatively thin with a range in thickness from ∼1 to 13 m. Once a flow is completely solidified (a rapid process on a geological time scale), the surface cools faster than the surface of a partially molten flow. Cooling via the base of the lava flow is also important and accelerates the solidification of the flow compared to the rate for the ‘semi-infinite’ approximation (which is only valid for very deep lava bodies). We introduce a new model which incorporates the solidification and basal cooling features. This model gives a superior reproduction of the cooling of the 1997 Pillan lava flows on Io observed by the Galileo spacecraft. We also use the new model to determine what observations are necessary to constrain lava emplacement style at Loki Patera. Flows exhibit different cooling profiles from that expected from a lava lake. We model cooling with a finite-element code and make quantitative predictions for the behavior of lava flows and other lava bodies that can be tested against observations both on Io and Earth. For example, a 10-m-thick ultramafic flow, like those emplaced at Pillan Patera in 1997, solidifies in ∼450 days (at which point the surface temperature has cooled to ∼280 K) and takes another 390 days to cool to 249 K. Observations over a sufficient period of time reveal divergent cooling trends for different lava bodies [examples: lava flows and lava lakes have different cooling trends after the flow has solidified (flows cool faster)]. Thin flows solidify and cool faster than flows of greater thickness. The model can therefore be used as a diagnostic tool for constraining possible emplacement mechanisms and compositions of bodies of lava in remote-sensing data.     

Predicting ground-motion variations at the Turkey-flat test site,California

Results of a single group participating in an international experiment are analyzed. The experiment served to verify computational predictions of the ground-motion variations due to near-surface geological effects at a site established for that purpose by the California Department of Conservation. Based on an acceleration record at a rock location, and geotechnical model of medium, records at the other locations of a nearby sedimentary deposit were predicted. A 2-D finite-difference sensitivity analysis suggested that the lateral wave-propagation effects are negligibly small, and locally 1-D computations are sufficient for the present site. Those computations are compared with observations not available to the authors during the blind prediction. Peak accelerations, peak velocities and RMS accelerations were predicted with errors less than 159%, 114% and 62%, respectively. Maxima of the response spectra were fitted within a factor of 2. The predicted and observed Husid's plots (i.e., the normalized cumulative plots of the acceleration squared) have the correlation coefficients 0.98. The detected misfits do not show any simple relation to the instrument location, component, frequency, or time.     

Integrated approach for identification of potential groundwater zones in Seethanagaram Mandal of Vizianagaram District,Andhra Pradesh,India

Identifying a good site for groundwater exploration in hard rock terrain is a challenging task. In hard rocks, groundwater occurs in secondary porosity developed due to weathering, fracturing, faulting, etc., which is highly variable within short distance and contributing to near-surface inhomogeneity. In such situations topographic, hydrogeological and geomorphological features provide useful clues for the selection of suitable sites. Initially, based on satellite imagery, topographical, geomorphological and hydrogeological features, an area of about 149 km² was demarcated as a promising zone for groundwater exploration in the hard rock tract of Seethanagaram Mandal, Vizianagaram District, Andhra Pradesh, India. A total of 50 Vertical Electrical Soundings (VES) were carried out using Wenner electrode configuration. An interactive interpretation of the VES data sharpened the information inferred from geomorphological and hydrogeological reconnaissance. Ten sites were recommended for drilling. Drilling with Down-The-Hole Hammer (DTH) was carried out at the recommended sites down to 50 to 70 m depths. The interpreted VES results matched well with the drilled bore well lithologs. The yields of bore wells vary from 900 to 9000 liters per hour (lph).     

山东省济宁强磁异常区深部铁矿初步验证及其意义

山东省济宁磁异常是一个重、磁同源体,面积大于100 km2,磁异常峰值为3800nT。钻探验证在孔深1041.57～1796.54m位置发现铁矿体,矿体总厚度74.04～220m,磁性铁平均品位15.89～25.19%。矿石类型有条带状方解磁铁石英岩和条带状磁铁石英大理岩,矿石的主要组成矿物为石英、方解石、磁铁矿、磁赤铁矿、菱铁矿。矿体产于济宁岩群浅变质岩系中,矿床特征与条带状铁建造（BIF）铁矿或鞍山式铁矿有明显区别,铁矿成因类型属与千枚岩、变质中酸性火山岩、大理岩有关的沉积变质型铁矿床。该区铁矿资源潜力巨大。     

Applications of particle-tracking techniques to bank infiltration: a case study from El Paso,Texas, USA

This paper presents results of a small scale study that utilized particle-tracking techniques to evaluate transport of river water through an alluvial aquifer in a bank infiltration testing site in El Paso, Texas, USA. The particle-tracking survey was used to better define filtration parameters. Several simulations were generated to allow visualization of the effects of well placement and pumping rate on flow paths, travel time, the size of the pumping influence zone, and proportion of river-derived water and groundwater mixing in the pumping well. Simulations indicate that migration of river water into the aquifer is generally slow. Most water does not arrive at the well by the end of an 18-day pumping period at 0.54 m³/min pumping rate for a well located 18 m from the river. Forty-four percent of the water pumped from the well was river water. The models provided important information needed to design appropriate sampling schedules for bank filtration practices and ensured meeting adequate soil-retention times. The pumping rate has more effect on river water travel time than the location of the pumping well from the river. The examples presented in this paper indicate that operating the pumping well at a doubled distance from the river increased the time required for the water to travel to the well, but did not greatly change the capture zone.