Large to intermediate-scale aquifer heterogeneity in fine-grain dominated alluvial fans (Cenozoic As Pontes Basin,northwestern Spain): insight based on three-dimensional geostatistical reconstruction     

O. Falivene  L. Cabrera  A. Sáez 《Hydrogeology Journal》2007,15(5):861-876

Facies reconstructions are used in hydrogeology to improve the interpretation of aquifer permeability distribution. In the absence of sufficient data to define the heterogeneity due to geological processes, uncertainties in the distribution of aquifer hydrofacies and characteristics may appear. Geometric and geostatistical methods are used to understand and model aquifer hydrofacies distribution, providing models to improve comprehension and development of aquifers. However, these models require some input statistical parameters that can be difficult to infer from the study site. A three-dimensional reconstruction of a kilometer scale fine-grain dominated Cenozoic alluvial fan derived from more than 200 continuously cored, closely spaced, and regularly distributed wells is presented. The facies distributions were reconstructed using a genetic stratigraphic subdivision and a deterministic geostatistical algorithm. The reconstruction is only slightly affected by variations in the geostatistical input parameters because of the high-density data set. Analysis of the reconstruction allowed identification in the proximal to medial alluvial fan zones of several laterally extensive sand bodies with relatively higher permeability; these sand bodies were quantified in terms of volume, mean thickness, maximum area, and maximum equivalent diameter. These quantifications provide trends and geological scenarios for input statistical parameters to model aquifer systems in similar alluvial fan depositional settings.  相似文献   

Predicting DNAPL entrapment and recovery: the influence of hydraulic property correlation     

L. D.?LemkeEmail author  L. M.?Abriola 《Stochastic Environmental Research and Risk Assessment (SERRA)》2003,17(6):408-418

The influence of aquifer property correlation on multiphase fluid migration, entrapment and recovery was explored by incorporating correlated and uncorrelated porosity, permeability, and capillary pressure-saturation (P_c-Sat) parameter fields in a cross-sectional numerical multiphase flow model. Comparison of two-dimensional entrapped organic saturation distributions for a simulated tetrachloroethylene (PCE) spill in ensembles of aquifer realizations suggests that the degree of spatial correlation in P_c-Sat parameters exerts a controlling influence on dense nonaqueous phase liquid (DNAPL) spreading and redistribution in saturated aquifers. The predicted evolution of DNAPL source zones and resultant remediation efficiency under surfactant enhanced aquifer remediation (SEAR) also appear to be strongly influenced by the spatial correlation of aquifer parameters and multiphase flow constitutive relationships. Results for a limited number of realizations selected from each ensemble showed that removal of 60% to 99% of entrapped PCE could reduce dissolved contaminant concentration and mass flux by approximately two orders of magnitude under natural gradient conditions. Aqueous phase contaminant mass flux did not vary uniformly as a function of % DNAPL removed, however, and notable differences in behavior were observed for models incorporating correlated versus uncorrelated P_c-Sat and permeability fields. Although these results must be confirmed through analysis of additional realizations, it is likely that similar or larger differences between correlated and uncorrelated system behavior will be observed in aquifers with greater spatially variability than that of the nonuniform, homogeneous sand aquifer studied here. Funding for this research was provided by the United States Environmental Protection Agency, Great Lakes and Mid-Atlantic Center for Hazardous Substance Research under Grant No. R-825540, the Michigan Department of Environmental Quality under Contract No. Y80011, and the Strategic Environmental Research and Development Program under Project No. CU-1293. The content of this publication does not necessarily represent the views of these agencies and has not been subject to agency review.  相似文献   

Structured flow-based gridding and upscaling for modeling subsurface flow   总被引：2，自引：0，他引：2  

C. He  L.J. Durlofsky 《Advances in water resources》2006,29(12):1876-1892

  相似文献   

The role of surface heterogeneity in modelling the stable boundary layer     

A. McCabe  A. R. Brown 《Boundary-Layer Meteorology》2007,122(3):517-534

High-resolution numerical simulations are performed for three nights for each of two areas of the United Kingdom. Area-averaging techniques are then used to calculate effective stability functions for turbulence parametrizations in models using typical mesoscale and global spatial resolutions. Comparisons are made with parametrizations commonly used in numerical weather prediction models. The present results do not suggest that significant enhancement of the stability functions above 50 m is justified. Closer to the surface, more significant enhancement is observed in some regions. It is shown that the amount of enhancement is related to the variability of the orography.  相似文献   

A dense array experiment for the observation of waveform perturbations   总被引：1，自引：0，他引：1  

Stphane Gaffet  Christophe Larroque  Anne Deschamps  Florence Tressols 《Soil Dynamics and Earthquake Engineering》1998,17(7-8):475-484

We analyze observed three-component seismograms recorded by a dense array of 21 short-period receivers (f>2 Hz) installed on a homogeneous Quaternary plain with 1.5×3.5 km² horizontal dimensions. The experiment was conducted in the district of Caille, in SE France for 3 months. The observed waveform data are related to diffraction and refraction of the direct wavefield by different geological structures and topographic morphologies surrounding the array. The origin of the coda energy following the energetic Pg or Sg arrivals is either influenced by heterogeneities close to and inside of the valley, for instance thrust of Jurassic limestones over Cretaceous marls in the northern and southern limits of the valley, or by heterogeneities located outside of the valley such as the interface unconformity between crystalline basement and sedimentary cover in the northern border of the Maures–Tanneron ranges. The main energetic arrivals of the seismograms arrive at the array with a definite range of azimuths that is not correlated with the event back-azimuth. Moreover, these arrivals do not appear with azimuths related to the characteristic NS and EW valley orientations, and thus may be related to heterogeneities located outside of the valley.  相似文献   

岩体多尺度异质性及其力学行为     

兰恒星  包含  孙巍锋  刘世杰 《工程地质学报》2022,30(1):37-52

异质性是岩体的重要属性之一,从岩体的物质组成到结构特征,再到微宏观力学行为表现,异质性现象普遍存在。本文在梳理已有研究成果的基础上,从岩体的细微观异质性、结构面异质性和宏观异质性3个不同的尺度,对异质性的产生原因、表现特征、力学效应和时空演化规律等进行阐述,并介绍了相关研究方法以及异质性在不同尺度上的差异表现。岩体异质性的研究,是岩体工程地质力学理论体系的延伸,而多尺度的研究方法也成为全面揭示岩体变形破坏机制的一种有效手段。  相似文献   

Boundary-Layer Adjustment Over Small-Scale Changes of Surface Heat Flux   总被引：1，自引：0，他引：1  

L.?MahrtEmail author  Dean?Vickers 《Boundary-Layer Meteorology》2005,116(2):313-330

Four months of eddy correlation data collected over a grass field and a nearby sage brush community are analyzed to examine the adjustment of the boundary-layer structure as it flows from the heated brush to the snow-covered grass. The grass site includes a 34-m tower with seven levels of eddy correlation data. The midday heat flux over the snow-covered grass and bare ground surfaces is often downward particularly with melting conditions, while the corresponding heat flux over the brush is almost always upward. For most of these cases, a stable internal boundary layer over the snow is well defined in terms of vertical profiles of the buoyancy flux over the snow-covered grass. The stable internal boundary layer is generally embedded within a deeper layer of flux divergence corresponding to increasing upward heat flux with height above the internal boundary layer. With thin snow cover, the surface heat flux over the grass is weak upward due to heating of grass protruding above the snow so that the flow adjusts to a decrease of the upward surface heat flux in the downwind direction. This common case of an adjusting boundary layer contrasts with the formation of an internal boundary layer due to a change of sign of the surface heat in flux the downwind direction. The adjustment of the boundary layer to the decrease of the surface heat flux leads to vertical divergence of the upward heat flux in contrast to the usual heated boundary layer over homogeneous surfaces. The consequences of the cooling due to the vertical divergence of the heat flux are discussed in terms of the heat budget of the adjusting and internal boundary layers.  相似文献   

Eulerian spatial moments for solute transport in three-dimensional heterogeneous,dual-permeability media     

J.?Xu  B.X.?HuEmail author 《Stochastic Environmental Research and Risk Assessment (SERRA)》2004,18(2):47-60

A Eulerian analytical method is developed for nonreactive solute transport in heterogeneous, dual-permeability media where the hydraulic conductivities in fracture and matrix domains are both assumed to be stochastic processes. The analytical solution for the mean concentration is given explicitly in Fourier and Laplace transforms. Instead of using the fast fourier transform method to numerically invert the solution to real space (Hu et al., 2002), we apply the general relationship between spatial moments and concentration (Naff, 1990; Hu et al., 1997) to obtain the analytical solutions for the spatial moments up to the second for a pulse input of the solute. Owing to its accuracy and efficiency, the analytical method can be used to check the semi-analytical and Monte Carlo numerical methods before they are applied to more complicated studies. The analytical method can be also used during screening studies to identify the most significant transport parameters for further analysis. In this study, the analytical results have been compared with those obtained from the semi-analytical method (Hu et al., 2002) and the comparison shows that the semi-analytical method is robust. It is clearly shown from the analytical solution that the three factors, local dispersion, conductivity variation in each domain and velocity convection flow difference in the two domains, play different roles on the solute plume spreading in longitudinal and transverse directions. The calculation results also indicate that when the log-conductivity variance in matrix is 10 times less than its counterpart in fractures, it will hardly influence the solute transport, whether the conductivity field is matrix is treated as a homogeneous or random field.  相似文献   

降水气候强迫下非均匀地表区域平均径流的一种参数化方案     

刘晶淼  丁裕国  周秀骥  李云 《气象学报》2009,67(1):28-34

将任一中尺度区域的平均瞬间径流率考虑为区域平均降水量和地表土壤层水分渗透垦的余项.根据降水量在地理空间上分布的实测资料拟合其空间概率密度函数(PDF),并结合土壤入渗物理过程的数学描述及其经验公式,精确估计出地表土壤渗透率及其空间分布,由此建立区域地表径流率的统计-动力学估计方案.换言之,区域内地表产流率可视为区域平均降水量与区域平均的土壤下渗量之差值,而区域内土壤的平均下渗量又町分为非饱和区和饱和区两部分的下渗量来分别计算.就陆面水分循环的物理过程而言,地表入渗现象是在一定的下垫面特性基础上,由一定的水分供应源而形成的.根据大气降水向地表层输送水分的物理过程,在满足植被表层覆盖需水(截流水)和地表层土壤人渗水基础上,多余的降水量才会形成地表径流.凶此,推求地表产流率的主要关键在于地表土壤层需水量.为此奉文根据土壤水分通量方程推导出水分入渗公式.又从描述土壤水分和降水的空间PDF出发,推导出非均匀土壤含水量及降水气候强迫所形成的次网格尺度区域平均径流率计算公式.利用长江三角洲地区1996年降水量和土壤特性等实测资料建立区域平均地表径流率的估计公式,并对其影响凶素进行敏感性试验.结果表明,该方法与用Mosaic方法计算的区域径流率(或产流率)结果十分接近.由此可见,该文提出的降水气候强迫下非均匀地表区域平均径流的这种统计-动力参数化方案,具有相当的可靠性与可行性.  相似文献   

Analysis of small-scale heterogeneity in clastic rocks by using computerized X-ray tomography (CT)     

Janos Geiger  Zoltan Hunyadfalvi  Peter Bogner   《Engineering Geology》2009,103(3-4):112-118

The aim of this study is to recognize the basic clastic rock types on the basis of Hounsfield Units (HUs) by using statistical methods (hypothesis tests, distribution fitting, and confidence intervals). How does this recognition depend on depositional history of samples tested? Does the numerical pattern of HUs obtained from small-scale analyzing of some particular sedimentary structures coincide to the textural expectations of those sedimentary structures? Are the HU values provided by the CT measurements capable of evidencing micro-cycles belonging to a particular depositional history?For testing the numerical recognition, we analyze macroscopically homogeneous argillaceous marl, siltstone and fine sandstone samples with similar depositional background, age, and degree of diagenesis. The preliminary research shows that the Hounsfield Units can identify different sedimentary rocks considerably well. There are significant differences among the confidence intervals belonging to the different rock types. However, there also have been some slight overlaps among them.For analyzing the effects of the depositional processes in case of a particular rock type, samples from channel sandstone, distributary mouth-bar, massive (structureless) sandstone from a delta-fed turbidity fan, and channel sandstone with traction carpet origin from a delta-fed turbidity fan (all coming from Pannonian basin filling series) are used. The results prove again the significant differences in the terms of both confidence intervals and distribution-types; however, some overlapping also occurs.  相似文献