首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 312 毫秒
1.
The neutron moisture probe is widely applicable to vadose zone monitoring problems which require measuring variable moisture contents. Neutron data are proportional to hydrogen density (modified by local chemistry) and sensitive to wetting fronts as well as changing volumes of hydrocarbon liquids. They cannot, however, be used to confirm contaminant chemistry, nor to detect steady-state flow. Neutron data are amenable to statistical analysis, providing a measure of the significance of data variations. Detection of incipient moisture changes at numerous monitoring locations is more practical using raw neutron data than data calibrated for moisture content because calibrations suffer from uncertainties associated with soil heterogeneities. When properly applied, the neutron probe is an effective monitoring tool as illustrated by three example applications described in this paper: (1) neutron moisture logs are used to detect subtle lithologic changes and identify monitoring horizons; (2) sequential neutron data are used to track induced saturation at a soil flushing pilot study; and (3) neutron logs from a horizontal access tube beneath a waste facility are used to pinpoint moisture anomalies.  相似文献   

2.
Preferred infiltration is mainly perceived as vertically down whereas subsurface storm flow is thought to occur parallel to slopes. The transition from vertical to lateral flow in a layered hillslope soil is the focus of the contribution. Transient flow is assumed to move as a wetting front. Three time‐domain reflectometry (TDR) wave‐guides, each 0·15 m long, were mounted in the shape of a truncated tetrahedron with its peak pointing down. Each wave‐guide focuses the front velocity along its axis. The three front‐velocity vectors are decomposed into their x, y and z components, which are then assembled to the resultant velocity vector. The volume density flux of preferred flow is the product of the front velocity and the mobile water content. The latter is the amplitude of transient soil moisture measured with each wave‐guide. The resultant vector of the volume flux density is computed similarly to the velocity vector. The experimental approach allows for the rapid assessment of transient flows without relying on the variation of water potentials. The experiments indicate that the directions of the resultant vectors of velocity and volume flux density can be estimated if the moisture variations of the three TDR wave‐guides are strongly correlated during the passing of the wetting front. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

3.
The saturated and unsaturated flow properties of a field soil under two tillage treatments were obtained with ponded rings and disc permeameters of dissimilar radii. No difference was observed between tillage treatments but the flow properties displayed a distinct macropore-matrix dichotomy, with K changing by an order of magnitude as ψ0 went from just - 30 mm to zero. Accurate prediction of time to incipient ponding was achieved using both numerical and analytical models calibrated with field hydraulic properties that were characteristic of the soil matrix. However, extension of the numerical model to the prediction of the wetting front development under non-ponded conditions was less accurate due to localized preferred wetting It is hypothesized that at this site, localized concentration of rainfall and hence, preferred wetting, May, occur by interception and stemflow associated with lines of standing stubble present in the original seeding slots.  相似文献   

4.
Interplay between capillary, gravity and viscous forces in unsaturated porous media gives rise to a range of complex flow phenomena affecting morphology, stability and dynamics of wetting and drainage fronts. Similar average phase contents may result in significantly different fluid distribution and patterns affecting macroscopic transport properties of the unsaturated medium. The formulation of general force balance within simplified pore spaces yields scaling relationships for motion of liquid elements in which gravitational force in excess of capillary pinning force scales linearly with viscous force. Displacement fluid front morphology is described using dimensionless force ratios expressed as Bond and Capillary numbers. The concise representations of a wide range of flow regimes with scaling relations, and predictive capabilities of front morphology based on dimensionless numbers lend support to certain generalizations. Considering available experimental data, we are able to define conditions for onset of unstable and intermittent flows leading to enhanced liquid and gas entrapment. These results provide a basis for delineation of a tentative value of Bo ∼ 0.05 as an upper limit of applicability of the Richards equation (at pore to sample scales) and related continuum-based flow models.  相似文献   

5.
砂岩核磁共振响应模拟及受限扩散   总被引:3,自引:1,他引:2       下载免费PDF全文
本文运用随机游走方法模拟了砂岩储层中流体的核磁共振(NMR)响应及其受限扩散现象.通过改变数字岩心的分辨率模拟生成不同孔隙尺寸的砂岩,研究了不同孔隙尺寸砂岩饱含水时流体扩散系数随扩散时间的变化关系,同时模拟了砂岩饱和单相流体和两相流体的NMR响应;研究了流体的受限扩散系数与横向弛豫时间T2的关系,分析了表面弛豫率和胶结指数对润湿相流体受限扩散系数线位置的影响,并将其用于解释砂岩储层的D-T2分布.结果表明:孔隙流体的扩散系数会随扩散时间的增加而逐渐减小并趋于定值.随着岩石孔隙尺寸的减小,受限扩散现象越明显,受限扩散对岩石NMR响应的影响也越大.润湿相流体受限扩散系数线的位置受岩石胶结指数和表面弛豫率的影响较大.由于润湿相流体扩散系数减小,导致D-T2分布中润湿相流体信号偏离其自由扩散系数线,需要利用流体的受限扩散系数线准确识别D-T2分布中的润湿相流体.  相似文献   

6.
The fate and transport of contaminants in the vicinity of septic fields remains poorly understood in many hydrogeomorphological environments. We report hydrometric data from an intensive hillslope‐scale experiment conducted between 29 August and 11 November 1998 at a residential leach field in New York State. The objective of our study was to characterize water flux within the vadose zone, understand the physical controls on the flux, and predict how this ultimately will affect subsurface water quality. Soil‐water flux was calculated using matric potential measurements from a network of 25 tensiometer nests, each nest consisting of three tensiometers installed to depths of 10, 50 and 130 cm. Unsaturated hydraulic conductivity curves were derived at each depth from field‐determined time‐domain reflectometry–tensiometry moisture‐release curves and borehole permeametry measurements. Flownets indicated that a strong upward flux of soil water occurred between rainstorms. Following the onset of (typically convective) rainfall, low near‐surface matric potentials were rapidly converted to near‐saturated and saturated conditions, promoting steep vertical gradients through the near‐surface horizons of the hillslope. Lateral hydraulic gradients were typically 10 times smaller than the vertical gradients. Resultant flow vectors showed that the flux was predominantly vertical through the vadose zone, and that the flux response to precipitation was short‐lived. The flux response was controlled primarily by the shape of the unsaturated hydraulic conductivity curves, which indicated a rapid loss of conductivity below saturation. Thus, soil water had a very high residence time in the vadose zone. The absence of rapid wetting at 130 cm and the delayed and small phreatic zone response to rainfall indicated that water movement through macropores did not occur on this hillslope. These results are consistent with a Cl tracing experiment, which demonstrated that the tracer was retained in the vadose zone for several months after injection to the system. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

7.
Flood irrigation is globally one of the most used irrigation methods. Typically, not all water that is applied during flood irrigation is consumed by plants or lost to evaporation. Return flow, the portion of applied water from flood irrigation that returns back to streams either via surface or subsurface flow, can constitute a large part of the water balance. Few studies have addressed the connection between vertical and lateral subsurface flows and its potential role in determining return flow pathways due to the difficulty in observing and quantifying these processes at plot or field scale. We employed a novel approach, combining induced polarization, time‐lapse electrical resistivity tomography, and time‐lapse borehole nuclear magnetic resonance, to identify flow paths and quantify changes in soil hydrological conditions under nonuniform application of flood irrigation water. We developed and tested a new method to track the wetting front in the subsurface using the full range of inverted resistivity values. Antecedent soil moisture conditions did not play an important role in preferential flow path activation. More importantly, boundaries between lithological zones in the soil profile were observed to control preferential flow pathways with subsurface run‐off occurring at these boundaries when saturation occurred. Using the new method to analyse time‐lapse resistivity measurements, we were able to track the wetting front and identify subsurface flow paths. Both uniform infiltration and preferential lateral flows were observed. Combining three geophysical methods, we documented the influence of lithology on subsurface flow processes. This study highlights the importance of characterizing the subsurface when the objective is to identify and quantify subsurface return flow pathways under flood irrigation.  相似文献   

8.
Modeling unsaturated flow in porous media requires constitutive relations that describe the soil water retention and soil hydraulic conductivity as a function of either potential or water content. Often, the hydraulic parameters that describe these relations are directly measured on small soil cores, and many cores are needed to upscale to the entire heterogeneous flow field. An alternative to the forward upscaling method using small samples are inverse upscaling methods that incorporate soft data from geophysical measurements observed directly on the larger flow field. In this paper, we demonstrate that the hydraulic parameters can be obtained from cross borehole ground penetrating radar by measuring the first arrival travel time of electromagnetic waves (represented by raypaths) from stationary antennae during a constant flux infiltration experiment. The formulation and coupling of the hydrological and geophysical models rely on a constant velocity wetting front that causes critical refraction at the edge of the front as it passes by the antennae. During this critical refraction period, the slope of the first arrival data can be used to calculate (1) the wetting velocity and (2) the hydraulic conductivity of the wet (or saturated) soil. If the soil is undersaturated during infiltration, then an estimate of the saturated water content is needed before calculating the saturated hydraulic conductivity. The hydraulic conductivity value is then used in a nonlinear global optimization scheme to estimate the remaining two parameters of a Broadbridge and White soil.  相似文献   

9.
This study demonstrates that positron emission tomography (PET) and neutron radiography (NR) techniques are complementary methods for determining the fluid pathway and porosity in crystalline rock. After preliminary injection of an organic solvent (e.g. isopropanol) front followed by the injection of the polymer solution (e.g. epoxy used for both techniques) and resin hardening, rock cutting may be performed. Flow pathway may be imaged by using a β+ emitter (e.g. 68Ga) in the resin. With a high-resolution PET camera, determination of the original water carrier features is possible in granodiorite pieces 20 cm in size and in simulated features with porosities of the order of 0.2. The use of a β+ tracer and the camera field, however, limit the lateral resolution of the technique (10 mm). Neutron radiography makes it possible to visualize the simulated porous phases by neutron transmission. The transmission process depends on the neutron scattering properties of the hydrogen-rich material (e.g. epoxy resin). Combination of 2D pictures may rebuild the 3D pattern. Lateral resolution may be in the range of 1 mm; however, the thickness of the rock sample must not exceed 10 cm. Complementarity of these techniques is discussed and they are compared with other methods used to determine porosity.  相似文献   

10.
Soils are often structured as fine-porous aggregates separated by large inter-aggregate pores. Under unsaturated conditions, water is mostly stored in the aggregates and water flow depends on the properties of the aggregates as well as on those of the contacts between aggregates. The goal of this study is to model and evaluate the hydraulic properties of the contacts. We used neutron radiography to monitor the infiltration of water through series of aggregates. The flow process was numerically simulated by considering the hydraulically conducting contact area between aggregates as a variable that depends on the capillary pressure. This contact area was evaluated by matching the observed and simulated water flow across aggregates. We determined the conductivity of the contacts assuming that it scales with the contact area. We also measured the equivalent conductivity of series of aggregates. We found that during drainage the hydraulically conducting contact area drastically decreases and the conductivity of the contacts becomes much smaller than that of the aggregates. We also found that the equivalent conductivity of the aggregate series decreases as the conductivity of the contacts. We concluded that the contacts control the flow: they are highly conductive when wet, but act as bottle-necks under drained conditions. The abrupt transition between these two limiting cases indicates that the contact region is more rapidly drained than the aggregate interior. Our findings might be extended to describe water transfer to an evaporating boundary, infiltration and storage of water in unsaturated aggregated soils.  相似文献   

11.
The loss of P in overland flow from most cultivated soils is controlled by erosion, and in‐turn soil moisture. We evaluated the effect of soil moisture on erosion and P transport in overland flow by applying rainfall (7 cm h?1) to packed soil boxes (1 m long and 0·15 m wide) and field plots (1 and 10 m long by 1 m wide) of silt loams in a central Pennsylvania (USA) catchment. Flow from packed soil boxes took longer to initiate as antecedent soil moisture decreased from field capacity (2 min) to air dried (8 to 9 min). Even in the more complex field plots (i.e. soil heterogeneity and topography), the wetter site (1 by 10 m plot; 70% field capacity) produced flow more quickly (3 min) and in greater volume (439 L) than the drier site (1 by 10 m plot; 40% field capacity, 15 min, and 214 L, respectively). However, less suspended sediment was transported from wetter soil boxes (1·6 to 2·5 g L?1) and field plots (0·9 g L?1) than drier boxes (2·9 to 4·2 g L?1) and plots (1·2 g L?1). Differences are attributed to their potential for soil aggregate breakdown, slaking and dispersion, which contribute to surface soil sealing and crusting, as dry soils are subject to rapid wetting (by rainfall). During flow, selective erosion and antecedent moisture conditions affected P transport. At field capacity, DRP and PP transport varied little during overland flow. Whereas P transport from previously dry soil decreased rapidly after the initiation of flow (6 to 1·5 mg TP L?1), owing to the greater slaking and dispersion of P‐rich particles into flow at the beginning than end of the flow event. These results indicate that soil moisture fluctuations greatly effect erosion and P transport potential and that management to decrease the potential for loss should consider practices such as conservation tillage and cover crops, particularly on areas where high soil P and erosion coincide. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

12.
Previous studies have shown that water retention curve (WRC) and the hydraulic conductivity vary because of changes of the void ratio or porosity of soil. However, limited documents pointed out the change of hydraulic properties of soil when compacted to different porosities while considering both of the drying and wetting processes of the WRC. This information is sometimes necessary for research like finger flow analysis or the occurrence of wetting and drying cycles as what would be seen in the field. Therefore, this study aims to examine the change of WRC characteristics with varied porosity considering both of the drying and wetting path in WRC by conducting a sand box experiment. Results show that the same type of sand compacted to various porosities have different hydraulic parameters. Hydraulic conductivities generally decrease with reduced porosities; shape parameter α of the van Genuchten equation (1980) linearly decreases with declining porosity and shape parameter n in a reversal manner for the sands of interest whether in the drying process or wetting process. The unsaturated properties of sand are further characterized by inspecting the variations of moisture content, matric suction and vertical displacement of soil body subject to periodic changes of the water level by another sand box experiment. The outcomes suggest that the saturated water content and residual water content are changing during the wetting–drying process, which can be an implication of the changed properties of WRC. The characteristics of volumetric deformation might be varied as well because of the observation of the dissimilar patterns of the changing vertical displacements among each wetting–drying process. Infiltration patterns of the sands also are identified through numerical modelling by introducing a constant infiltration flux from the surface followed by a no‐influx condition. Results indicate that less water accumulates in the sand near the surface for the sand compacted to higher porosity, but water can move deeper. Hydraulic conductivity is found as the prime factor dominating the evolvement of wetting fronts. However, shape parameters of water retention curves also affect the infiltration pattern to some extent. In addition, different sands with similar porosities can have quite different infiltrating characteristics. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

13.
The steady seepage rates at large times from flat-bottomed channels and channels of semi-circular cross-section were determined in laboratory tank experiments using various sands. Good agreement was found with the theoretical relationships which assume the flow to be confined to a saturated region bounded by a capillary-fringe surface and to be uniform and vertical at great depths. The steady large-time seepage rates were also obtained in laboratory sand-tank experiments for the three-dimensional cases of seepage from circular shallow ponds and hemispherical sources. These agreed with relationships obtained using an electrolytic tank analogue with approximate boundary conditions assumed for the flow region. A method of analysis of large-time seepage measurements from irrigation channels and infiltrometer rings is proposed, which yields the hydraulic conductivity and pressure head at the wetting front from experiments with different size channels or rings.  相似文献   

14.
北京平原区西北部大地热流与深部地温分布特征   总被引:3,自引:0,他引:3       下载免费PDF全文
北京平原区蕴藏着丰富的中-低温水热型地热资源,其西北部分布着小汤山地热田和京西北地热田,两大地热田以南口—孙河断裂带为界.地热田及其外围地区基础的地热地质研究工作较少.为给地热学研究和地热资源精细勘探提供科学依据,本文基于前人23眼钻孔的温度测量数据以及近期完成的548件热导率和100件放射性生热率实测数据,研究了区域大地热流和0~4 km深部地温特征.结果表明:(1)研究区现今地温梯度为11.31~94.89℃·km-1,平均值为31.79℃·km-1;岩石热导率为0.895~5.111 W·(m·K)-1,放射性生热率为0.257~2.305 μW·m-3,大地热流为48.1~99.1 mW·m-2,平均值为68.3 mW·m-2,热流的分布受基底形态和断裂构造控制.研究区东部南口—孙河断裂带两侧小汤山和郑各庄地区为高热流异常区,中部马池口地区也存在局部高热流异常区.(2)在南口—孙河断裂带的不同位置,不同深度地层温度差异明显,体现出区域现今地温场不只受控于该活动断裂,更是多期次构造事件复合叠加的结果.(3)南口—孙河断裂带南侧存在两处有意义的较高地温异常区,分别为郑各庄异常区和马池口异常区,其中马池口异常区是未来地热开发利用有一定潜力的地区.  相似文献   

15.
Stormwater infiltration systems are a popular method for urban stormwater control. They are often designed using an assumption of one‐dimensional saturated outflow, although this is not very accurate for many typical designs where two‐dimensional (2D) flows into unsaturated soils occur. Available 2D variably saturated flow models are not commonly used for design because of their complexity and difficulties with the required boundary conditions. A purpose‐built stormwater infiltration system model was thus developed for the simulation of 2D flow from a porous storage. The model combines a soil moisture–based model for unsaturated soils with a ponded storage model and uses a wetting front‐tracking approach for saturated flows. The model represents the main physical processes while minimizing input data requirements. The model was calibrated and validated using data from laboratory 2D stormwater infiltration trench experiments. Calibrations were undertaken using five different combinations of calibration data to examine calibration data requirements. It was found that storage water levels could be satisfactorily predicted using parameters calibrated with either data from laboratory soils tests or observed water level data, whereas the prediction of soil moistures was improved through the addition of observed soil moisture data to the calibration data set. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

16.
There is a need to elucidate the impact of ethanol on the subsurface environment because of the application of ethanol as automotive fuel. This study quantifies the effects of changes in surface tension, viscosity, and density induced by ethanol on the transmission and retention of water in the vadose zone. The HYDRUS‐1D model was modified to simulate two different scenarios of flow in a sandy loam involving ponding (constant head) or spillage with a subsequent rainfall event (constant flux). Solutions containing different amounts of the highly miscible ethanol (10, 50, and 100% by weight) as well as pristine water were considered. During ponding, ethanol reduced the amount of fluid entering the soil and slowed down the advancement of the wetting front. Viscosity effects were predominant for this scenario, reducing the average depth of the infiltrating liquid up to 44%. The total amount of pure ethanol that entered the soil was 11.38 vs. 17.64 cm for pure water. For the spillage scenario, the results suggest that density has little impact on the liquid movement. Surface tension effects are predominant in the upper portion of the soil. The changes in hydraulic conductivity due to ethanol‐induced modifications of solution viscosity are responsible for the slower advancement of the moisture front. The 10% ethanol solution moved 43.1% faster than pure ethanol during the first 2 d because of viscosity and surface tension effects.  相似文献   

17.
本文利用藏北地区三口天然气水合物钻孔测温数据,在分析样品热导率测试结果基础上,计算了藏北地区的热流值.对于样品热导率值,首先根据样品孔隙度对实验室测试结果进行了饱水校正,计算热流时采用的是对应井段的岩石热导率饱水校正值的厚度加权平均值.地温梯度以三口钻孔48 h的测温数据为基础,回归三口井的地温梯度,计算时去除了浅部受地表温度和冻土带对温度影响的数值.A钻孔地温梯度分为200~438 m和438~882 m两段回归,分段热流的加权平均值作为钻孔热流值,计算结果为42.7 mW·m-2; B钻孔和C钻孔回归地温梯度时未分段,热流计算结果分别为58.3 mW·m-2、70 mW·m-2.综合分析认为,岩石圈断裂、地幔上涌、碰撞造山过程中的剪切生热等因素可能造成了班公湖—怒江缝合带以南热流值较高,而北部羌塘地块热流值相对较低.  相似文献   

18.
Rohit Salve 《水文研究》2011,25(18):2907-2915
As the scope of hydrologic investigations extend deeper into the subsurface profile, and increasingly include fractured rock, there is a growing need for techniques that can accurately monitor saturation changes at a high spatial and temporal resolution in this environment. We have developed a technique, the Electrical Resistance Sensor Array System (ERSAS), to track moisture dynamics in vadose zone regions that include both fractured rock and soil. The performance of ERSAS was compared with the time domain reflectometry (TDR) technique under controlled and field conditions. We found that ERSAS was effective in determining patterns of saturation changes along vertical soil/rock profiles. Because of the small size of individual sensors, it was able to resolve travel times associated with a wetting front and peak saturation better than TDR. In addition, ERSAS is significantly cheaper than the TDR system, and the sensor arrays are relatively easier to install in the subsurface profile. Published in 2011 by John Wiley & Sons, Ltd.  相似文献   

19.
两淮煤田大地热流分布及其构造控制   总被引:2,自引:1,他引:1       下载免费PDF全文
基于127块煤系地层岩石样品的热导率测定结果,并结合59个井田内可靠的系统测温数据,计算得出两淮煤田的大地热流值,并编制大地热流分布图,其结果表明:两淮煤田大地热流值变化范围为29.7~83.9mW·m-2,平均值为58.3mW·m-2,和其他沉积盆地存在一定的差异,且淮南煤田大地热流值(63.7mW·m-2)远大于淮北煤田(55.2mW·m-2).综合分析得出,两淮煤田大地热流与其他盆地的差异以及淮南煤田热流值高于淮北煤田的现象为构造演化和区域地质背景的控制结果;而研究区内热流的分布不均主要是由于受地质构造对地温场的影响所致,推覆构造上下盘现今热流值的差异尤为突出.  相似文献   

20.
In this paper, the feasibility of using magnetic resonance imaging (MRI) to study water infiltration into a heterogeneous soil is examined, together with its difficulties and limitations. MRI studies of ponded water infiltration into an undisturbed soil core show that the combination of one- and two-dimensional imaging techniques provides a visual and non-destructive means of monitoring the temporal changes of soil water content and the moisture profile, and the movement of the wetting front. Two-dimensional images show air entrapment in repetitive ponded infiltration experiments. During the early stages of infiltration, one-dimensional images of soil moisture profiles clearly indicate preferential flow phenomena. The observed advance of wetting fronts can be described by a linear relationship between the square root of infiltration time (√t) and the distance of the wetting front from the soil surface. Similarly, the cumulative infiltration is also directly proportional to √t. Furthermore, from the MRI infiltration moisture profiles, it is possible to estimate the parameters that feature in infiltration equations. © 1997 by John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号