共查询到20条相似文献，搜索用时 437 毫秒
1.
Manuel Valera Zhengyang Guo Priscilla Kelly Sean Matz Vito Adrian Cantu Allon G. Percus Jeffrey D. Hyman Gowri Srinivasan Hari S. Viswanathan 《Computational Geosciences》2018,22(3):695-710
Structural and topological information play a key role in modeling flow and transport through fractured rock in the subsurface. Discrete fracture network (DFN) computational suites such as dfnWorks (Hyman et al. Comput. Geosci. 84, 10–19 2015) are designed to simulate flow and transport in such porous media. Flow and transport calculations reveal that a small backbone of fractures exists, where most flow and transport occurs. Restricting the flowing fracture network to this backbone provides a significant reduction in the network’s effective size. However, the particle-tracking simulations needed to determine this reduction are computationally intensive. Such methods may be impractical for large systems or for robust uncertainty quantification of fracture networks, where thousands of forward simulations are needed to bound system behavior. In this paper, we develop an alternative network reduction approach to characterizing transport in DFNs, by combining graph theoretical and machine learning methods. We consider a graph representation where nodes signify fractures and edges denote their intersections. Using random forest and support vector machines, we rapidly identify a subnetwork that captures the flow patterns of the full DFN, based primarily on node centrality features in the graph. Our supervised learning techniques train on particle-tracking backbone paths found by dfnWorks, but run in negligible time compared to those simulations. We find that our predictions can reduce the network to approximately 20% of its original size, while still generating breakthrough curves consistent with those of the original network. 相似文献
2.
3.
基于离散裂隙网络模型的裂隙水渗流计算 总被引：1，自引：1，他引：0
离散裂隙网络模型(Discrete Fracture Network（DFN）)是研究裂隙水渗流最为有效的手段之一。文章根据裂隙几何参数和水力参数的统计分布,利用Monte Carlo随机模拟技术生成二维裂隙网络,基于图论无向图的邻接矩阵判断裂隙网络的连通,利用递归算法提取出裂隙网络的主干网或优势流路径。基于立方定律和渗流连续性方程,利用数值解析法建立了二维裂隙网络渗流模型,分析不同边界条件下裂隙网络中的流体流动。结果表明,该方法可以模拟区域宏观水力梯度和边界条件下,裂隙网络水力梯度方向总的流量,以及节点的水位、节点间的流量和流动方向的变化特征,为区域岩溶裂隙水渗流计算提供了一种实用、可行的方法。 相似文献
4.
Sung-Hoon Ji Kyung Woo Park Doo-Hyun Lim Chunsoo Kim Kyung Su Kim William Dershowitz 《Hydrogeology Journal》2012,20(7):1341-1353
The development and implementation of a hybrid discrete fracture network/equivalent porous medium (DFN/EPM) approach to groundwater flow at the Gyeong-Ju low- and intermediate-level radioactive waste (LILW) disposal site in the Republic of Korea is reported. The geometrical and hydrogeological properties of fractured zones, background fractures and rock matrix were derived from site characterization data and implemented as a DFN. Several DFN realizations, including the deterministic fractured zones and the stochastic background fractures, whose statistical properties were verified by comparison with in-situ fracture and hydraulic test data, were suggested, and they were then upscaled to continuums using a fracture tensor approach for site-scale flow simulations. The upscaled models were evaluated by comparison to in-situ pressure monitoring data, and then used to simulate post-closure hydrogeology for the LILW facility. Simulation results demonstrate the importance of careful characterization and implementation of fractured zones. The study highlighted the importance of reducing uncertainty regarding the properties and variability of natural background fractures, particularly in the immediate vicinity of repository emplacement. 相似文献
5.
G. Lecca C. Lai F. Murgia R. Biddau L. Fanfani P. Maggi 《Earth Science Informatics》2009,2(1-2):83-95
AQUAGRID is the subsurface hydrology computational service of the Sardinian GRIDA3 infrastructure, designed to deliver complex environmental applications via a user-friendly Web portal. The service aims to provide to water professionals integrated modeling tools to solve water resources management problems and aid decision making for contaminated soil and groundwater. In this paper, the AQUAGRID application concept and enabling technologies are illustrated. At the heart of the service are the computational models to simulate complex and large groundwater flow and contaminant transport problems and geochemical speciation. AQUAGRID is built on top of compute-Grid technologies by means of the EnginFrame Grid framework. Distributed data management is provided by the Storage Resource Broker data-Grid middleware. The resulting environment allows end-users to perform groundwater simulations and to visualize and interact with their results, using graphs, 3D images and annotated maps. The problem solving capability of the platform is demonstrated using the results of two case studies deployed. 相似文献
6.
Fracture flow modelling based on satellite images of the Wajid Sandstone, Saudi Arabia 总被引：1，自引：0，他引：1
Conny Zeeb Daniel Göckus Paul Bons Hussain Al Ajmi Randolf Rausch Philipp Blum 《Hydrogeology Journal》2010,18(7):1699-1712
Large-scale geological features have been identified by satellite imagery and global positioning system data in the Wajid Sandstone in Saudi Arabia. The main objective is to evaluate the importance of fractures for the overall flow behaviour in this fractured rock aquifer and to estimate in-situ hydraulic apertures. Data on fractures and lineaments were available for three outcrops. By applying a “cut-out” routine on the fracture endpoint data of these fracture trace windows, three deterministic discrete fracture networks (DFN), with an area of 100 m?×?100 m, could be generated. These were used to simulate the fracture flow and to determine the hydraulic conductivity tensors. Using additional data on hydraulic pumping tests and matrix conductivities, in-situ hydraulic apertures could be determined. Average in-situ hydraulic apertures range from 1,300 to 1,700 µm. Observations from the field support these results. In addition, a hydraulic conductivity ratio between the matrix and fracture system was used to identify the contribution of the DFN to the overall fluid transport. A ratio of 10.4 was determined, which indicates that the effective flow behaviour in the Wajid Sandstone aquifer is not entirely dominated by the fracture system, though evidently strongly controlled by it. 相似文献
7.
Groundwater flow in fractured rocks is modeled using a coupled model based on the domain decomposition method. In the model, the fractured porous medium is divided into two non-overlapping sub-domains. One is the rock matrix, in which the medium is described using a continuum model. The other consists of deep fractures and fissure zones, where the medium is described using a discrete fracture network (DFN) model. The two models are coupled through the continuity of the hydraulic heads and fluxes on the common boundaries. The coupled model is used to simulate groundwater flow in a hydropower station. The results show that the model simulates groundwater levels that are in agreement with the measured groundwater levels. Furthermore, the model’s parameters relating to deep fractures and fissure zones are verified by comparing three different models (the continuum model, coupled model, and DFN model). The results show that the coupled model can capture and duplicate the hydrogeological conditions in the study domain, whereas the continuum model overestimates and the DFN model underestimates the measured hydraulic heads. A sensitivity analysis shows that fracture aperture has a considerable effect on the groundwater level. So, when the fracture aperture is large, the coupled model or DFN model is more appropriate than the continuum model in the fracture domain. 相似文献
8.
对大空间尺度和长时间跨度的地下水流动及污染物质运移进行预测的需求, 使水文地质研究面临异乎寻常的挑战.这些需求来自于对核废料地质储放方法的安全性评价、地下水污染状况评价及其治理方案的选择.流动系统的非均质性是地下水流动及物质运移模拟中最主要的困难之一, 这种困难来自对非均质系统进行特征描述(通过原位观测实现)、概念化及模拟.评述了非均质介质中流动运移模拟的一些重要问题与挑战, 讨论了解决的途径.讨论的主题包括: 动力流动的沟道化, 示踪剂穿透曲线, 裂隙岩石中流体流动的多尺度, 观测的不同尺度, 模拟、预测与非均质性以及系统特征描述和预测性模拟的分析. 相似文献
9.
In fractured-rock aquifers, the geometric and hydraulic properties of the fractures commonly have a dominant influence on transport. Tracer tests are often used to estimate directly the gross transport properties of a fractured rock mass. The prospects for understanding characteristics of the heterogeneities in a fractured porous medium were explored from evidence provided by tracer experiments. The approach was to simulate flow and transport on a large set of prescribed fracture networks in a two-dimensional homogeneous permeable medium, thus generating synthetic tracer test data. The fracture orientation, aperture, spacing and network geometry were systematically altered from one case to the next. A classification scheme was devised for the tracer breakthrough curves using principal component analysis and this classification was linked to the fracture pattern properties. Even under highly simplified and controlled conditions, quite different fracture patterns can produce very similar breakthrough curves. The classification scheme thus demonstrates that a single breakthrough curve cannot reveal the fracture geometry with any precision. However, the scheme provided a methodology for rejecting geometric properties that do not belong to the fracture pattern under investigation, thus reducing the uncertainty in fracture geometry. 相似文献
10.
We present a discussion of the state-of-the-art on the use of discrete fracture networks (DFNs) for modelling geometrical characteristics, geomechanical evolution and hydromechanical (HM) behaviour of natural fracture networks in rock. The DFN models considered include those based on geological mapping, stochastic generation and geomechanical simulation. Different types of continuum, discontinuum and hybrid geomechanical models that integrate DFN information are summarised. Numerical studies aiming at investigating geomechanical effects on fluid flow in DFNs are reviewed. The paper finally provides recommendations for advancing the modelling of coupled HM processes in fractured rocks through more physically-based DFN generation and geomechanical simulation. 相似文献
11.
Hui-Hai Liu Sumit Mukhopadhyay Nicolas Spycher Burton M. Kennedy 《Hydrogeology Journal》2011,19(6):1151-1160
Precipitation-dissolution reactions are important for a number of applications such as isotopic tracer transport in the subsurface. Analytical solutions have been developed for tracer transport in both single-fracture and multiple-fracture systems associated with these reactions under transient and steady-state transport conditions. These solutions also take into account advective transport in fractures and molecular diffusion in the rock matrix. For studying distributions of disturbed tracer concentration (the difference between actual concentration and its equilibrium value), effects of precipitation-dissolution reactions are mathematically equivalent to a “decay” process with a decay constant proportional to the corresponding bulk reaction rate. This important feature significantly simplifies the derivation procedure by taking advantage of the existence of analytical solutions for tracer transport associated with radioactive decay in fractured rock. It is also useful for interpreting tracer breakthrough curves, because the impact of a decay process is relatively easy to analyze. Several illustrative examples are presented, which show that the results are sensitive to fracture spacing, matrix diffusion coefficient (fracture surface area), and bulk reaction rate (or “decay” constant), indicating that the relevant flow and transport parameters may be estimated by analyzing tracer signals. 相似文献
12.
Multi-borehole pumping and tracer tests on the 10 to 100-m scale were conducted in a fractured chalk aquitard in the Negev
Desert, Israel. Outcrop and core fracture surveys, as well as slug tests in packed-off intervals, were carried out at this
site to obtain the parameters needed for construction of a stochastic discrete fracture network (DFN). Calibration of stochastic
DFNs directly to the multiple borehole test data was inadequate. Instead, two equivalent deterministic DFN flow models were
used: the vertical-fractures (VF) model, consisting of only vertical fractures, and the fractures’ intersections (INT) model,
consisting of vertical and horizontal fractures with enhanced transmissivity at their intersections. Both models were calibrated
against the multi-borehole response of one pumping test and their predictions were tested against three other independent
pumping tests. The average accuracies of all transient drawdown predictions of the VF and INT models were 65 and 66%, respectively.
In contrast to this equality in average drawdown predictions of both models, the INT model predicted better important breakthrough
curve features (e.g., first and peak arrival times), than the VF model. This result is in line with previously assumed channeled
flow, derived from analytical analysis of these pumping and tracer tests.
Ronit Nativ, deceased, may her memory be blessed. 相似文献
13.
Analysis of contaminant transport through fractured crystalline rocks has received considerable attention, particularly with regard to subsurface nuclear waste repositories. Most of the studies have employed the dual continuum approach, with the fractures and the rock matrix as the two continuums, assuming that fractures control the overall conductivity of the rock and the porous matrix just provides storage. However, field observations of rock fractures have shown that the real situation can be very complex. Based on some recent investigations, it has been reported that the portion of the rock matrix adjacent to many open fractures is physically and chemically altered. These alterations, referred to as the fracture skin, can have different sorption and diffusion properties compared to those of the undisturbed rock matrix and this may influence the transport of solutes through such formations. In the present study, a numerical model is developed to simulate conservative solute transport in a fractured crystalline rock formation using the triple continuum approach ?? with the fracture, fracture skin and the rock matrix as the three continuums. The model is solved using a fully implicit finite difference scheme. Contaminant migration in the fractured formation with and without skin has been simulated. It is observed that contaminant penetration along the fracture is enhanced at large flow velocities. The effect of flow velocity on conservative solute transport is investigated for different fracture apertures and fracture skin thicknesses. The influence of flow velocity on contaminant transport is demonstrated to be more with change in fracture aperture than with change in skin thickness. 相似文献
14.
Simulation of subsurface heterogeneity is important for modeling subsurface flow and transport processes. Previous studies have indicated that subsurface property variations can often be characterized by fractional Brownian motion (fBm) or (truncated) fractional Levy motion (fLm). Because Levy-stable distributions have many novel and often unfamiliar properties, studies on generating fLm distributions are rare in the literature. In this study, we generalize a relatively simple and computationally efficient successive random additions (SRA) algorithm, originally developed for generating Gaussian fractals, to simulate fLm distributions. We also propose an additional important step in response to continued observations that the traditional SRA algorithm often generates fractal distributions having poor scaling and correlation properties. Finally, the generalized and modified SRA algorithm is validated through numerical tests. 相似文献
15.
Artificial tracer experiments were conducted in the mature karst system of Jeita (Lebanon) under various flow conditions using surface and subsurface tracer injection points, to determine the variation of transport parameters (attenuation of peak concentration, velocity, transit times, dispersivity, and proportion of immobile and mobile regions) along fast and slow flow pathways. Tracer breakthrough curves (TBCs) observed at the karst spring were interpreted using a two-region nonequilibrium approach (2RNEM) to account for the skewness in the TBCs’ long tailings. The conduit test results revealed a discharge threshold in the system dynamics, beyond which the transport parameters vary significantly. The polynomial relationship between transport velocity and discharge can be related to the variation of the conduit’s cross-sectional area. Longitudinal dispersivity in the conduit system is not a constant value (α?=?7–10 m) and decreases linearly with increasing flow rate because of dilution effects. Additionally, the proportion of immobile regions (arising from conduit irregularities) increases with decreasing water level in the conduit system. From tracer tests with injection at the surface, longitudinal dispersivity values are found to be large (8–27 m). The tailing observed in some TBCs is generated in the unsaturated zone before the tracer actually arrives at the major subsurface conduit draining the system. This work allows the estimation and prediction of the key transport parameters in karst aquifers. It shows that these parameters vary with time and flow dynamics, and they reflect the geometry of the flow pathway and the origin of infiltrating (potentially contaminated) recharge. 相似文献
16.
Estimating the hydraulic properties of fractured aquifers is challenging due to the complexity of structural discontinuities that can generally be measured at a small scale, either in core or in outcrop, but influence groundwater flow over a range of scales. This modeling study uses fracture scanline data obtained from surface bedrock exposures to derive estimates of permeability that can be used to represent the fractured rock matrix within regional scale flow models. The model is developed using PETREL, which traditionally benefits from high resolution data sets obtained during oil and gas exploration, including for example seismic data, and borehole logging data (both lithological and geophysical). The technique consists of interpreting scanline fracture data, and using these data to generate representative Discrete Fracture Network (DFN) models for each field set. The DFN models are then upscaled to provide an effective hydraulic conductivity tensor that represents the fractured rock matrix. For each field site, the upscaled hydraulic conductivities are compared with estimates derived from pumping tests to validate the model. A hydraulic conductivity field is generated for the study region that captures the spatial variability of fracture networks in pseudo-three dimensions from scanline data. Hydraulic conductivities estimated using this approach compare well with those estimated from pumping test data. The study results suggest that such an approach may be feasible for taking small scale fracture data and upscaling these to represent the aquifer matrix hydraulic properties needed for regional groundwater modeling. 相似文献
17.
Tomás Morales-Juberías Martín Olazar Jose María Arandes Pedro Zafra Iñaki Antigüedad Félix Basauri 《Environmental Geology》1997,30(3-4):143-151
A model based on numerical solutions, which allows for solving the dispersion equation under variable recharge and velocity conditions, is developed to simulate solute transport in conduit flow aquifers during flow recession periods. As an example, the evolution of a tracer in the little known karst conduit that links the sinking stream of Oma valley to the Olalde spring is investigated in the karstic region of Santa Eufemia-Ereñozar (Basque Country, Spain). The model, with different hypothetical structures, allows for obtaining series of tracer breakthrough curves, which are fitted to experimental data using an optimization algorithm. These results, although they can be used to simulate the tracer evolution between the two points considered, do not allow for determining the internal structure and spatial disposition of contributions in the aquifer. 相似文献
18.
19.
Victor Ginting Felipe Pereira Michael Presho Shaochang Wo 《Computational Geosciences》2011,15(4):691-707
In this paper, we develop a procedure for subsurface characterization of a fractured porous medium. The characterization involves
sampling from a representation of a fracture’s permeability that has been suitably adjusted to the dynamic tracer cut measurement
data. We propose to use a type of dual-porosity, dual-permeability model for tracer flow. This model is built into the Markov
chain Monte Carlo (MCMC) method in which the permeability is sampled. The Bayesian statistical framework is used to set the
acceptance criteria of these samples and is enforced through sampling from the posterior distribution of the permeability
fields conditioned to dynamic tracer cut data. In order to get a sample from the distribution, we must solve a series of problems
which requires a fine-scale solution of the dual model. As direct MCMC is a costly method with the possibility of a low acceptance
rate, we introduce a two-stage MCMC alternative which requires a suitable coarse-scale solution method of the dual model.
With this filtering process, we are able to decrease our computational time as well as increase the proposal acceptance rate.
A number of numerical examples are presented to illustrate the performance of the method. 相似文献
20.
We present a fully implicit formulation of coupled flow and geomechanics for fractured three-dimensional subsurface formations. The Reservoir Characterization Model (RCM) consists of a computational grid, in which the fractures are represented explicitly. The Discrete Fracture Model (DFM) has been widely used to model the flow and transport in natural geological porous formations. Here, we extend the DFM approach to model deformation. The flow equations are discretized using a finite-volume method, and the poroelasticity equations are discretized using a Galerkin finite-element approximation. The two discretizations—flow and mechanics—share the same three-dimensional unstructured grid. The mechanical behavior of the fractures is modeled as a contact problem between two computational planes. The set of fully coupled nonlinear equations is solved implicitly. The implementation is validated for two problems with analytical solutions. The methodology is then applied to a shale-gas production scenario where a synthetic reservoir with 100 natural fractures is produced using a hydraulically fractured horizontal well. 相似文献