首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 859 毫秒
1.
In porous aquifers, groundwater flow and solute transport strongly depend on the sedimentary facies distribution at fine scale, which determines the heterogeneity of the conductivity field; in particular, connected permeable sediments could form preferential flow paths. Therefore, properly defined statistics, e.g. total and intrinsic facies connectivity, should be correlated with transport features. In order to improve the assessment of the relevance of this relationship, some tests are conducted on two ensembles of equiprobable realizations, obtained with two different geostatistical simulation methods—sequential indicator simulation and multiple point simulation (MPS)—from the same dataset, which refers to an aquifer analogue of sediments deposited in a fluvial point-bar/channel association. The ensembles show different features; simulations with MPS are more structured and characterised by preferential flow paths. This is confirmed by the analysis of transport connectivities and by the interpretation of data from numerical experiments of conservative solute transport with single and dual domain models. The use of two ensembles permits (1) previous results obtained for single realizations to be consolidated on a more firm statistical basis and (2) the application of principal component analysis to assess which quantities are statistically the most relevant for the relationship between connectivity indicators and flow and transport properties.  相似文献   

2.
The fine-scale heterogeneity of porous media affects the large-scale transport of solutes and contaminants in groundwater and it can be reproduced by means of several geostatistical simulation tools. However, including the available geological information in these tools is often cumbersome. A hierarchical simulation procedure based on a binary tree is proposed and tested on two real-world blocks of alluvial sediments, of a few cubic meters volume, that represent small-scale aquifer analogs. The procedure is implemented using the sequential indicator simulation, but it is so general that it can be adapted to various geostatistical simulation tools, improving their capability to incorporate geological information, i.e., the sedimentological and architectural characterization of heterogeneity. When compared with a standard sequential indicator approach on bi-dimensional simulations, in terms of proportions and connectivity indicators, the proposed procedure yields reliable results, closer to the reference observations. Different ensembles of three-dimensional simulations based on different hierarchical sequences are used to perform numerical experiments of conservative solute transport and to obtain ensembles of equivalent pore velocity and dispersion coefficient at the scale length of the blocks (meter). Their statistics are used to estimate the impact of the variability of the transport properties of the simulated blocks on contaminant transport modeled on bigger domains (hectometer). This is investigated with a one-dimensional transport modeling based on the Kolmogorov-Dmitriev theory of branching stochastic processes. Applying the proposed approach with diverse binary trees and different simulation settings provides a great flexibility, which is revealed by the differences in the breakthrough curves.  相似文献   

3.
张嘉  王明玉 《地学前缘》2010,17(6):152-158
在地下水污染模拟预报中,弥散参数是很难确定的一个模型参数。因实验室小尺度弥散规律一般不能用于大尺度弥散过程,而野外示踪试验却耗资大、周期长,限制了其实用性。文中利用随机数值模拟手段、基于随机理论的蒙特卡罗方法及序贯高斯模拟技术来生成渗透系数随机场,并研究渗透系数对数场的方差、相关长度以及变异函数类型在不同尺度上对纵向弥散度的影响,进而建立纵向弥散度与随机分布渗透系数场的方差和相关长度的统计定量关系,并与Gelhar理论计算结果进行比较。数值模拟结果表明,经过一定迁移距离后纵向弥散度与随机分布渗透系数对数场的方差和相关长度具有良好的线性统计关系,与Gelhar理论公式表达的关系类型类似。但对于较大的方差,纵向弥散度模拟结果明显大于Gelhar理论计算值,而对于较大相关长度在迁移距离不很大时,纵向弥散度模拟结果明显小于Gelhar理论计算值。本研究可为野外大尺度地下水污染预报模型中水动力弥散参数的确定提供方法借鉴。  相似文献   

4.
非均质含水层中渗流与溶质运移研究进展   总被引:1,自引:0,他引:1  
Natural aquifer heterogeneity controls the groundwater flow and solute transport, and how to accurately quantify the flow and solute transport in heterogeneous aquifers has received wide attention by many scholars, and has become a hot research topic in earth science. Theoretically, a systematic review is given by the following aspect: flow and solute transport model, moment analysis, multi scale analysis. The resolved and remained issues for scale conversion in hydrogeology research are pointed out. Secondly, recent advances of heterogeneous field test, uncertainty and velocity connectivity are analyzed. Finally, the geophysical inversion of aquifer heterogeneity, stochastic theory and development of stochastic simulation software, scale conversion and uncertainty of velocity connectivity, and the relationship between heterogeneity and hydrogeological condition on the major four aspects of the future research direction is pointed out.  相似文献   

5.
6.
Saturated flow takes place in geological formations of spatially variable permeability which is regarded as a stationary random space function of given statistical moments. The flow is assumed to be uniform in the mean and the Eulerian velocity field has stationary fluctuations. Water carries solutes that react according to the nonlinear equilibrium Freundlich isotherm. Neglecting pore scale dispersion (high Peclet number), we study the behavior of an initially finite pulse injection of constant concentration. Mean flux-averaged concentration is derived in a simple manner by using the previously determined solution of transport in a homogeneous one-dimensional medium and the Lagrangian methodology developed by Cvetkovic and Dagan [5] to model reactive transport in a three-dimensional flow field. The mean breakthrough curves are computed and the combined effect of reactive parameters and heterogeneity upon reduction of the concentration peak is investigated. Furthermore, with the aid of temporal moments, we determine equivalent reaction and macrodispersion coefficients pertinent to a homogeneous medium. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

7.
根据渗透系数在多相岩层中的协方差函数,文中首先推导了用以描述溶质在三维模型中运移特征的宏观弥散系数方程,然后导出了纵向弥散系数用来表达该参数的尺度效应。利用一个例子,讨论了宏观弥散随时间的变化趋势,以及同一岩相内转移项和越相转移项对弥散的相对贡献。灵敏度分析指出纵向弥散系数是与岩相平均长度正相关,也与不同岩相的渗透系数的差值正相关,其取值随渗透系数的总方差的增加而增大。  相似文献   

8.
Ma  Lei  Liao  Huasheng  Qian  Jiazhong  Zhao  Weidong  Li  Shuguang 《Hydrogeology Journal》2023,31(4):985-1004

The heterogeneity of hydrofacies is represented as spatial variability on different scales, and it has a significant impact on the behavior of groundwater flow and pollutant transport. However, effectively characterizing hydrofacies heterogeneity on different scales remains one of the most challenging problems in hydrogeology. In this study, an upscaling hydrofacies simulation (UHS) framework is proposed by integrating the upscaling borehole generalization (UBG) approach and transition probability geostatistics (TPG). A new UBG approach for generating virtual boreholes with equivalent hydrofacies information based on relatively high-density borehole lithological data is proposed, and the TPG is used to delineate the multiscale facies distribution. The results show that the UBG approach can significantly reduce borehole data volume while retaining the key equivalent hydrofacies information on a coarser scale. The UHS method can well characterize the overall distribution of equivalent hydrofacies on coarser scales, with the minor-component hydrofacies underestimated and the major-component hydrofacies overestimated to a lesser extent, and more equivalent facies appearing in strong heterogeneous areas. These results demonstrate that the UHS method can provide valuable capacity insights and advantages in characterizing hydrofacies heterogeneity on different scales using such high-density borehole lithological data.

  相似文献   

9.
Advances in photogrammetry have eased the acquisition of high-resolution digital information from outcrops, enabling faster, non-destructive data capturing and improved reservoir modeling. Geocellular models for flow dynamics with in the virtual outcrop in siliciclastic deposits at different sets of sandstone facies architecture remain, however, a challenge. Digital maps of bedding, lithological contrast, spatial-temporal variations of bedding and permeability characteristics make it more easy to understand flow tortuosity in a particular architecture. An ability to precisely model these properties can improve reservoir characterization and flow modeling at different scales. Here we demonstrate the construction of realistic 2 D sandstone facies based models for a pragmatic simulation of flow dynamics using a combination of digital point clouds dataset acquired from LiDAR and field investigation of the Sandakan Formation, Sabah, Borneo.Additionally, we present methods for enhancing the accuracy of outcrop digital datasets for producing high resolution flow simulation. A well-exposed outcrop from the Sandakan Formation, Sabah, northwest Borneo having a lateral extent of 750 m was chosen in order to implement our research approach. Sandstone facies and its connectivity are well constrained by outcrop observations, data from air-permeability measurements, bilinear interpolation of permeability, grid construction and water vector analysis for flow dynamics.These proportions were then enumerated in terms of static digital outcrop model(DOM) and facies model based on sandstone facies bedding characteristics. Flow simulation of water vector analysis through each of the four sandstone facies types show persistent spatial correlation of permeability that align with either cross-bedded orientation or straight with more dispersion high quality sandstone(porosity 21.25%-41.2%and permeability 1265.20-5986.25 mD) and moderate quality sandstone(porosity 10.44%-28.75% and permeability 21.44-1023.33 mD). Whereas, in more heterolithic sandstone(wavy-to flaser-bedded and bioturbated sandstone), lateral variations in permeability show spatially non-correlated patterns over centimeters to tens of meters with mostly of low quality sandstone(porosity 3.4%-12.31% and permeability < 1 mD to 3.21 mD). These variations reflect the lateral juxtaposition in flow dynamics. It has also been resulted that the vertical connectivity and heterogeneities in terms of flow are mostly pragmatic due to the interconnected sandstone rather than the quality of sandstone.  相似文献   

10.
11.
The approach of dynamic stratigraphy aims to understand genetic processes that form stratigraphic units in a hierarchy of spatial and temporal scales. This approach was used to investigate Quaternary gravel deposits in terms of their sedimentology and in order to characterize the various sedimentary units in terms of their hydrogeological properties. Facies analysis within 62 gravel pits, laboratory permeability measurements of field samples and geophysical surveys (3-D georadar, 2-D seismic reflection) led to the detection and classification of sedimentary heterogeneity according to the following six scales whereby each scale can be translated into defined hydrostratigraphic units. (1) Particles and pores (micro scale) that reflect depositional and diagenetic fluid dynamics as well as source material behaviour (e.g. grain-size, roundness, lithological composition). This was found to be important for the hydrogeochemistry of groundwater in gravel aquifers (e.g. higher sorption capacity of carbon-rich limestone particles for organic pollutants). (2) Strata (meso scale) contain the recognition of sorting, fabric, texture and stratinomic features, which can give an indication of transport and depositional dynamics. Five major lithofacies groups, for example, were distinguished within fluvial gravel-bed deposits. Their variable hydraulic properties led to their subdivision into 12 hydrofacies types. They form the smallest mappable hydrostratigraphic units, which may result in either preferred pathways for fluid flow or flow barriers. (3) Depositional elements (macro scale) enable reconstruction of sedimentary/geomorphic elements and their dynamics within a depositional system (e.g. gravel-bed braided river systems are dominated by gravel sheet, gravel dunes and scour pool depositional elements). Hydrostratigraphically, the architecture of depositional elements influences the hydraulic connectivity and local permeability structure/distribution within an aquifer body. Five types of depositional elements in fluvial gravel-bed deposits were distinguished and their geometries/dimensions quantified. (4) Facies bodies (mega scale) composed of a stack of depositional elements and strata recording distinct environmental systems and their dynamics (e.g. a coarse-grained prograding delta system). Hydrostratigraphically, facies bodies represent major compartments of an aquifer. Six major types of meltwater-controlled facies bodies were identified in the study area. (5) Genetic sequences (mega scale) reflect the shifts of depositional environments caused by allocyclic changes (e.g. glacial advance recorded by a coarsening upward sequence) or autocyclic changes of landscape shaping events. These sequences may form separate hydrostratigraphic units or aquifer storeys. (6) Basin fill (giga scale) comprising the lateral and vertical stacking of facies bodies and genetic sequences controlled by either long-term glacier dynamics or short term flood events. The regional distribution of permeable gravel units and, for example, less permeable diamicts builds the larger scale hydrostratigraphy.  相似文献   

12.
13.
Analysis of groundwater chemistry can yield important insights about subsurface conditions, and provide an alternative and complementary method for characterizing basin hydrogeology, especially in areas where hydraulic data are limited. More specifically, hydrochemical facies have been used for decades to help understand basin flow and transport, and a set of facies were developed for the Roswell Artesian Basin (RAB) in a semi-arid part of New Mexico, USA. The RAB is an important agricultural water source, and is an excellent example of a rechargeable artesian system. However, substantial uncertainties about the RAB hydrogeology and groundwater chemistry exist. The RAB was a great opportunity to explore hydrochemcial facies definition. A set of facies, derived from fingerprint diagrams (graphical approach), existed as a basis for testing and for comparison to principal components, factor analysis, and cluster analyses (statistical approaches). Geochemical data from over 300 RAB wells in the central basin were examined. The statistical testing of fingerprint-diagram-based facies was useful in terms of quantitatively evaluating differences between facies, and for understanding potential controls on basin groundwater chemistry. This study suggests the presence of three hydrochemical facies in the shallower part of the RAB (mostly unconfined conditions) and three in the deeper artesian system of the RAB. These facies reflect significant spatial differences in chemistry in the basin that are associated with specific stratigraphic intervals as well as structural features. Substantial chemical variability across faults and within fault blocks was also observed.  相似文献   

14.
油气储层相控随机建模技术的约束方法   总被引:31,自引:4,他引:27  
提出根据相控随机建模技术描述油气储层特征,并用已知地质数据和沉积微相的研究成果,如平面展布方向、宽厚比等来指导相控建模,文章提出“多层分级控制,同级套合管理”的相控随机建模策略,并将相序、概率、定量知识库或变差函数三个方面相结合来表征油气储层的非均质特性。具体则是从沉积形成与演化的成因角度指导沉积储层随机建模的过程,应用多参数协同、分层次约束的方法,运用沉积相带的平面展布和垂向演化来控制建模的结果,并用宝力格油田的实际地质和钻井数据验证优选得到的多个模型实现。结果表明:优选模型真实地反映了地下油气储层的非均质特性和连通性展布特征,并经新钻井数据验证确实有效。  相似文献   

15.
Combining groundwater flow models with solute transport models represents a common challenge in groundwater resources assessments and contaminant transport modeling. Groundwater flow models are usually constructed at somewhat larger scales (involving a coarser discretization) to include natural boundary conditions. They are commonly calibrated using observed groundwater levels and flows (if available). The groundwater solute transport models may be constructed at a smaller scale with finer discretization than the flow models in order to accurately delineate the solute source and the modeled target, to capture any heterogeneity that may affect contaminant migration, and to minimize numerical dispersion while still maintaining a reasonable computing time. The solution that is explored here is based on defining a finer grid subdomain within a larger coarser domain. The local-grid refinement (LGR) implemented in the Modular 3D finite-difference ground-water flow model (MODFLOW) code has such a provision to simulate groundwater flow in two nested grids: a higher-resolution sub-grid within a coarse grid. Under the premise that the interface between both models was well defined, a comprehensive sensitivity and uncertainty analysis was performed whereby the effect of a parameter perturbation in a coarser-grid model on transport predictions using a higher-resolution grid was quantified. This approach was tested for a groundwater flow and solute transport analysis in support of a safety evaluation of the future Belgian near-surface radioactive waste disposal facility. Our reference coarse-grid groundwater flow model was coupled with a smaller fine sub-grid model in two different ways. While the reference flow model was calibrated using observed groundwater levels at a scale commensurate with that of the coarse-grid model, the fine sub-grid model was used to run a solute transport simulation quantifying concentrations in a hypothetical well nearby the disposal facility. When LGR coupling was compared to a one-way coupling, LGR was found to provide a smoother flow solution resulting in a more CPU-efficient transport solution. Parameter sensitivities performed with the groundwater flow model resulted in sensitivities at the head observation locations. These sensitivities identified the recharge as the most sensitive parameter, with the hydraulic conductivity of the upper aquifer as the second most sensitive parameter in regard to calculated groundwater heads. Based on one-percent sensitivity maps, the spatial distribution of the observations with the highest sensitivities is slightly different for the upper aquifer hydraulic conductivity than for recharge. Sensitivity analyses were further performed to assess the prediction scaled sensitivities for hypothetical contaminant concentrations using the combined groundwater flow and solute transport models. Including all pertinent parameters into the sensitivity analysis identified the hydraulic conductivity of the upper aquifer as the most sensitive parameter with regard to the prediction of contaminant concentrations.  相似文献   

16.
《Sedimentology》2018,65(4):1246-1276
Submarine mass‐transport deposits represent important stratigraphic heterogeneities within slope and basinal sedimentary successions. A poor understanding of how their distribution and internal architecture affect the fluid flow migration pathway may lead to unexpected compartmentalization issues in reservoir analysis. Studies of modern carbonate mass‐transport deposits mainly focus on large seismic‐scale slope failures; however, the near‐platform basinal depositional environment often hosts mass‐transport deposits of various dimensions. The small‐scale and meso‐scale (metres to several tens of metres) carbonate mass‐transport deposits play a considerable role in distribution of sediment and therefore have an impact on the heterogeneity of the succession. In order to further constrain the geometry and internal architecture of mass‐transport deposits developed in near‐slope basinal carbonates, a structural and sedimentological analysis of sub‐seismic‐scale mass‐transport deposits has been undertaken on the eastern margin of the Apulian carbonate platform in the Gargano Promontory, south‐east Italy. These mass‐transport deposits, that locally comprise a large proportion (50 to 60%) of the base of slope to basinal sediments of the Cretaceous Maiolica Formation, typically display a vertically bipartite character, including debrites and slump deposits of varying volume ratios. A range of brittle and ductile deformation styles developed within distinct bed packages, together with the presence of both chert clasts, folded chert layers and spherical chert nodules, suggest that sediments were at different stages of lithification prior to downslope movement associated with mass‐transport deposits. This study helps elucidate the emplacement processes, frequency and character of subseismic‐scale mass‐transport deposits within the basinal carbonate environment, and thereby reduces the uncertainties in the characterization of subsurface carbonate geofluid reservoirs.  相似文献   

17.
阿舍勒-冲乎尔火山沉积盆地是一个以火山岩并伴有陆源和内源沉积物的堆积场所,即该区泥盆纪地层是以事件作用的产物为其主要成分特点,显然,对其古地理的恢复,相的划分采用传统的相序分析方法是不合适的。作者认为要阐明该区泥盆纪的古地理面貌,必须在区域地层研究基础上,综合分析各种标志和沉积物的特点,寻求它们的物理成因和地质成因解。  相似文献   

18.
Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.  相似文献   

19.
Multiple-point statistics (MPS) provides a flexible grid-based approach for simulating complex geologic patterns that contain high-order statistical information represented by a conceptual prior geologic model known as a training image (TI). While MPS is quite powerful for describing complex geologic facies connectivity, conditioning the simulation results on flow measurements that have a nonlinear and complex relation with the facies distribution is quite challenging. Here, an adaptive flow-conditioning method is proposed that uses a flow-data feedback mechanism to simulate facies models from a prior TI. The adaptive conditioning is implemented as a stochastic optimization algorithm that involves an initial exploration stage to find the promising regions of the search space, followed by a more focused search of the identified regions in the second stage. To guide the search strategy, a facies probability map that summarizes the common features of the accepted models in previous iterations is constructed to provide conditioning information about facies occurrence in each grid block. The constructed facies probability map is then incorporated as soft data into the single normal equation simulation (snesim) algorithm to generate a new candidate solution for the next iteration. As the optimization iterations progress, the initial facies probability map is gradually updated using the most recently accepted iterate. This conditioning process can be interpreted as a stochastic optimization algorithm with memory where the new models are proposed based on the history of the successful past iterations. The application of this adaptive conditioning approach is extended to the case where multiple training images are proposed as alternative geologic scenarios. The advantages and limitations of the proposed adaptive conditioning scheme are discussed and numerical experiments from fluvial channel formations are used to compare its performance with non-adaptive conditioning techniques.  相似文献   

20.
随着全球页岩油气勘探开发的不断深入,海相泥页岩作为非常规油气沉积学的重要研究对象之一,其沉积特征及沉积过程日益引起学者们的关注,并在岩相分类、沉积过程解释及层序地层划分方面取得了显著进展。海相泥页岩的岩相划分已达到厘米—毫米尺度,其命名原则主要考虑结构(粒度)、层理及矿物成分等特征,结构方面着重参考砂级、粗泥级和细泥级的占比,层理描述注重纹层的连续性、形态和几何关系,矿物成分则比较黏土矿物、石英和碳酸盐岩矿物的相对含量。海相泥页岩的沉积过程多样,浅海陆棚环境下细粒沉积物(岩)的主要搬运与沉积过程包括风成输入、异轻流、重力流(前三角洲浊流、异重流和浪控沉积物重力流)及风暴流,而在平坦而广泛的陆表海则以季风或潮汐驱动下的远岸底流搬运沉积为主。重力流及风暴流沉积一般会形成特征相似但各具特色的粒序层理,而在远岸底流作用下则会形成纹层状泥页岩。沉积特征的精细描述推动了泥页岩层序地层格架的建立,三级层序界面的识别在于寻找可对比的主要侵蚀面,而层序内准层序的划分除了要寻找次要侵蚀面外,还需依靠对精细描述的沉积特征进行定量化统计与分析。上述研究表明,沉积特征的精细化描述与定量分析,并结合水槽模拟实验,是今后海相泥页岩沉积过程研究的主体发展方向。  相似文献   

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

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