共查询到20条相似文献，搜索用时 15 毫秒
1.
Widely used numerical models of solute transport processes in subsurface aquifers are limited to nonlocally refined rectangular, or logically rectangular, structured grids. This presents an unsuitable option to efficient numerical simulations maintaining an acceptable level of accuracy. Optimal selection of locally refined cells for efficient solute transport models is challenging to the current generation of numerical models. We present a novel and relatively simple to implement algorithm addressing these shortcomings. This method operates in four steps involving travel times simulations, a grid coarsening stage followed by a selective local grid refinement based on a cell-wise indicator, and a final postprocessing step. The refinement index is the sum of weighted logarithmic distributions of scaled forward and backward travel times. We calculate representative flow and transport properties at the two scales of the composite grid with a flow-based upscaling technique. We present two test problems to demonstrate the performances of this new gridding algorithm. We obtain the most important speedups for composite grids generated with the highest indicator thresholds. When hydrodynamic dispersion effects increase, we obtain less important speedups. An important outcome of this work is that grid design depends on nature and strength of the underlying flow and solute transport processes. Therefore, we suggest developing solute transport workflows integrating this grid generation algorithm as an integral component to build comprehensive and efficient groundwater models. 相似文献
2.
We present a method to determine equivalent permeability of fractured porous media. Inspired by the previous flow-based upscaling methods, we use a multi-boundary integration approach to compute flow rates within fractures. We apply a recently developed multi-point flux approximation Finite Volume method for discrete fracture model simulation. The method is verified by upscaling an arbitrarily oriented fracture which is crossing a Cartesian grid. We demonstrate the method by applying it to a long fracture, a fracture network and the fracture network with different matrix permeabilities. The equivalent permeability tensors of a long fracture crossing Cartesian grids are symmetric, and have identical values. The application to the fracture network case with increasing matrix permeabilities shows that the matrix permeability influences more the diagonal terms of the equivalent permeability tensor than the off-diagonal terms, but the off-diagonal terms remain important to correctly assess the flow field. 相似文献
3.
We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells. 相似文献
4.
《Advances in water resources》2003,26(10):1041-1060
A new technique for generating coarse scale models of highly heterogeneous subsurface formations is developed and applied. The method uses generic global coarse scale simulations to determine the boundary conditions for the local calculation of upscaled properties (permeability or transmissibility). An iteration procedure assures consistency between the local and global calculations. Transport processes are simulated using a subgrid velocity reconstruction technique applied in conjunction with the local–global upscaling procedure. For highly heterogeneous (e.g., channelized) systems, the new method is shown to provide considerably more accurate coarse scale models for flow and transport, relative to reference fine scale results, than do existing local (and extended local) upscaling techniques. The applicability of the upscaled models for different global boundary conditions is also considered. 相似文献
5.
《Advances in water resources》1999,22(7):697-710
Flownets are useful tools for the visualization of groundwater flow fields. Using orthogonal flownets as grids for transport modeling is an effective way to control numerical dispersion, especially transverse to the direction of flow. Therefore tools for automatic generation of flownets may be seen both as postprocessors for groundwater flow simulations and preprocessors for contaminant transport models. Existing methods to generate streamline-oriented grids suffer from drawbacks such as the inability to include sources in the interior of the grid. In this paper, we introduce a new method for the generation of streamline-oriented grids which handles wells in the grid interior, and which produces orthogonal grids for anisotropic systems. Streamlines are generated from an accurate velocity field obtained from the solution of the mixed-hybrid finite element method for flow, while pseudopotentials, which are orthogonal to the streamlines, are obtained by a standard finite element solution of the pseudopotential equation. A comprehensive methodology for the generation of orthogonal grids, including the location of stagnation points and dividing streamlines, is introduced. The effectiveness of the method is illustrated by means of examples. A related paper presents a compatible formulation of the solution for reactive transport, while a second related paper gives a detailed quantitative assessment of the various forms of modelled mixing and their effect on the accuracy of simulations of the biodegradation of groundwater contaminants. 相似文献
6.
This article demonstrates that permeability upscaling, which can require complex techniques, is not necessary to significantly decrease the CPU time in reactive transport modeling. CPU time depends more on the geochemistry than the flow calculation. Flow rate upscaling is proposed as an alternate method to permeability upscaling, which is more suited to time-consuming flow resolution. To apply this method, a finite volume approach is most convenient.Considering the equality of flow as the equivalence criterion, when the coarse grid overlays the fine grid, flow rate upscaling leads, by construction, to the exact results, whereas the accuracy of permeability upscaling methods often depends on specific conditions. Some focus is put on the limitations of a common permeability upscaling technique, the simplified renormalization. In stationary flow, the gain in CPU time is the same for both flow rate upscaling and permeability upscaling. In transient flow, flow rate upscaling is slightly less time-efficient but the ratio between both CPU times decreases when the geochemistry is more complex.Working with an accurate flow rate field in the upscaled case reveals that porosity upscaling is a surprisingly tricky issue. Solution mixing is induced and residence times can be significantly affected. These changes have potentially important consequences on reactive transport modeling. They are not specific to the flow rate upscaling method; they are a general issue. Some simplified cases, assuming a homogeneous mineralogy, are examined. At this stage, a simple heuristic method is proposed, which yields reliable results under particular conditions (high heterogeneity). Porosity upscaling remains an open research field. 相似文献
7.
Subsurface formations are characterized by heterogeneity over multiple length scales, which can have a strong impact on flow and transport. In this paper, we present a new upscaling approach, based on time-of-flight (TOF), to generate upscaled two-phase flow functions. The method focuses on more accurate representations of local saturation boundary conditions, which are found to have a dominant impact (in comparison to the pressure boundary conditions) on the upscaled two-phase flow models. The TOF-based upscaling approach effectively incorporates single-phase flow and transport information into local upscaling calculations, accounting for the global flow effects on saturation, as well as the local variations due to subgrid heterogeneity. The method can be categorized into quasi-global upscaling techniques, as the global single-phase flow and transport information is incorporated in the local boundary conditions. The TOF-based two-phase upscaling can be readily integrated into any existing local two-phase upscaling framework, thus more flexible than local–global two-phase upscaling approaches developed recently. The method was applied to permeability fields with different correlation lengths and various fluid-mobility ratios. It was shown that the new method consistently outperforms existing local two-phase upscaling techniques, including recently developed methods with improved local boundary conditions (such as effective flux boundary conditions), and provides accurate coarse-scale models for both flow and transport. 相似文献
8.
In this work we propose upscaling method for nonlinear Forchheimer flow in heterogeneous porous media. The generalized Forchheimer law is considered for incompressible and slightly-compressible single-phase flows. We use recently developed analytical results (Aulisa et al., 2009) [1] and formulate the resulting system in terms of a degenerate nonlinear flow equation for the pressure with the nonlinearity depending on the pressure gradient. The coarse scale parameters for the steady state problem are determined so that the volumetric average of velocity of the flow in the domain on fine scale and on coarse scale are close. A flow-based coarsening approach is used, where the equivalent permeability tensor is first evaluated following streamline methods for linear cases, and modified in order to take into account the nonlinear effects. Compared to previous works (Garibotti and Peszynska, 2009) [2], (Durlofsky and Karimi-Fard) [3], this approach can be combined with rigorous mathematical upscaling theory for monotone operators, (Efendiev et al., 2004) [4], using our recent theoretical results (Aulisa et al., 2009) [1]. The developed upscaling algorithm for nonlinear steady state problems is effectively used for variety of heterogeneities in the domain of computation. Direct numerical computations for average velocity and productivity index justify the usage of the coarse scale parameters obtained for the special steady state case in the fully transient problem. For nonlinear case analytical upscaling formulas in stratified domain are obtained. Numerical results were compared to these analytical formulas and proved to be highly accurate. 相似文献
9.
Use of upscaled elevation and surface roughness data in two-dimensional surface water models 总被引：1，自引：0，他引：1
In this paper, we present an approach that uses a combination of cell-block- and cell-face-averaging of high-resolution cell elevation and roughness data to upscale hydraulic parameters and accurately simulate surface water flow in relatively low-resolution numerical models. The method developed allows channelized features that preferentially connect large-scale grid cells at cell interfaces to be represented in models where these features are significantly smaller than the selected grid size. The developed upscaling approach has been implemented in a two-dimensional finite difference model that solves a diffusive wave approximation of the depth-integrated shallow surface water equations using preconditioned Newton-Krylov methods. Computational results are presented to show the effectiveness of the mixed cell-block and cell-face averaging upscaling approach in maintaining model accuracy, reducing model run-times, and how decreased grid resolution affects errors. Application examples demonstrate that sub-grid roughness coefficient variations have a larger effect on simulated error than sub-grid elevation variations. 相似文献
10.
It has been shown that convective mixing in porous media flow is important for applications such as saltwater intrusion and geological storage of carbon dioxide. In the latter case, dissolution from the injected phase to the resident brine is assisted by convective mixing, which leads to enhanced storage security through reduced buoyancy. Here, we focus on the effect of horizontal barriers on the efficiency of convective mixing. Previous investigations of the effect of heterogeneity on mixing efficiency have focused on random permeability fields or barriers of small extent compared to the intrinsic finger wavelength. The effect of horizontal barriers of larger extent, such as mudstone inclusions or thin shale deposits, has not been given sufficient attention. We perform detailed numerical investigations to represent the continuous solution of this problem in semi-infinite domains with barriers arranged in a periodic manner. The results show that mass flux into the domain, which is a measure of the efficiency of redistribution of the solute, is inversely proportional to the barrier length and proportional to the horizontal and vertical aperture between the barriers, for the cases studied. The flow structure is complex, and it depends not only on the total area of barriers but also largely on the distribution of barriers. Therefore, neither simple analytical models nor simple upscaling methods that lack information about the flow paths, can be used to predict the behavior. However, we compute the effective vertical permeability by flow-based upscaling and show that it can be used to directly obtain a first-order approximation to the mass flux into the domain. 相似文献
11.
本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性. 相似文献
12.
On the generation of coastline-following grids for ocean models—trade-off between orthogonality and alignment to coastlines 总被引：1，自引：1，他引：0
Regional ocean models usually utilize orthogonal curvilinear grids that are fit to the coastline of the modeled regions. While the orthogonality of the grid is required from the perspective of the numerical algorithms, the alignment to the irregular coastlines improves the characterization of the land-sea distribution and the ocean simulation. In this article, we carry out fractal analysis of two representative coastal regions and discuss the trade-offs between the orthogonality and coastline alignment during the grid generation of these regions. A new grid generation method based on Schwarz-Christoffel conformal mappings is proposed, with automatic coastal boundary retrieval algorithm that generates resolution dependent boundary for grid generation and alleviates the human efforts involved in traditional methods. We show that for the southeastern Pacific region, the coastline is smooth with low fractal dimension and there exists effective trade-off with a coastline boundary that adjusts to the desired grid resolution. On the contrary, there is no effective trade-off for southeast China seas where the coastline is of higher fractal dimension, and a coarser coastline boundary is recommended for better orthogonality with little loss in coastline alignment. Further numerical study of coastal trapped Kelvin waves for the typical regions demonstrate that the new coastline-fitting grids achieve smaller error in numerical dispersion and higher accuracy. Through analysis, we conclude that for grid generation for regional ocean modeling, modelers should bring into consideration of the multi-scale fractal characteristics of the coastline. 相似文献
13.
Giovanni Battista Chirico Andrew W. Western Rodger B. Grayson Günter Blschl 《水文研究》2005,19(13):2539-2556
Specific catchment area (SCA) patterns are commonly computed on grids using flow direction algorithms that treat the flow as coming from a point source at the pixel centre. These algorithms are all ambiguous in the definition of the flow width to be associated with a pixel when computing the SCA. Different methods for computing the flow width have been suggested, without giving an objective reason. In the few cases where this issue has been specifically discussed, the flow width is derived from subjective analysis and incorrect conceptualizations. This paper evaluates alternative approaches for defining the flow width when computing SCA patterns using the D∞ and D8 algorithms, by comparing theoretical and computed SCA patterns on sloping planes, inward and outward cones. Two new methods of defining the flow width are also analysed for both the D∞ and D8 algorithms. The performances of the different methods are discussed in relation to two dimensionless parameters: (1) the global resolution, defined as the ratio of a characteristic length of the study area to the grid size and (2) the upslope area resolution, defined as the ratio of the theoretical SCA to the grid size. The optimal methods are identified by specific threshold values of these dimensionless parameters. We conclude that assuming the flow width invariant and equal to the grid size is generally the best approach in most practical circumstances, both for the D∞ and D8 algorithm. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
14.
The Fourier pseudospectral method has been widely accepted for seismic forward modelling because of its high accuracy compared to other numerical techniques. Conventionally, the modelling is performed on Cartesian grids. This means that curved interfaces are represented in a ‘staircase fashion‘causing spurious diffractions. It is the aim of this work to eliminate these non-physical diffractions by using curved grids that generally follow the interfaces. A further advantage of using curved grids is that the local grid density can be adjusted according to the velocity of the individual layers, i.e. the overall grid density is not restricted by the lowest velocity in the subsurface. This means that considerable savings in computer storage can be obtained and thus larger computational models can be handled. One of the major problems in using the curved grid approach has been the generation of a suitable grid that fits all the interfaces. However, as a new approach, we adopt techniques originally developed for computational fluid dynamics (CFD) applications. This allows us to put the curved grid technique into a general framework, enabling the grid to follow all interfaces. In principle, a separate grid is generated for each geological layer, patching the grid lines across the interfaces to obtain a globally continuous grid (the so-called multiblock strategy). The curved grid is taken to constitute a generalised curvilinear coordinate system, where each grid line corresponds to a constant value of one of the curvilinear coordinates. That means that the forward modelling equations have to be written in curvilinear coordinates, resulting in additional terms in the equations. However, the subsurface geometry is much simpler in the curvilinear space. The advantages of the curved grid technique are demonstrated for the 2D acoustic wave equation. This includes a verification of the method against an analytic reference solution for wedge diffraction and a comparison with the pseudospectral method on Cartesian grids. The results demonstrate that high accuracies are obtained with few grid points and without extra computational costs as compared with Cartesian methods. 相似文献
15.
16.
《Advances in water resources》2007,30(4):1027-1045
Streamline methods have shown to be effective for reservoir simulation. For a regular grid, it is common to use the semi-analytical Pollock’s method to obtain streamlines and time-of-flight coordinates (TOF). The usual way of handling irregular grids is by trilinear transformation of each grid cell to a unit cube together with a linear flux interpolation scaled by the Jacobian. The flux interpolation allows for fast integration of streamlines, but is inaccurate even for uniform flow. To improve the tracing accuracy, we introduce a new interpolation method, which we call corner-velocity interpolation. Instead of interpolating the velocity field based on discrete fluxes at cell edges, the new method interpolates directly from reconstructed point velocities given at the corner points in the grid. This allows for reproduction of uniform flow, and eliminates the influence of cell geometries on the velocity field. Using several numerical examples, we demonstrate that the new method is more accurate than the standard tracing methods. 相似文献
17.
Three-dimensional grids representing a heterogeneous, ground water system are generated at 10 different resolutions in support of a site-scale flow and transport modeling effort. These grids represent hydrostratigraphy near Yucca Mountain, Nevada, consisting of 18 stratigraphic units with contrasting fluid flow and transport properties. The grid generation method allows the stratigraphy to be modeled by numerical grids of different resolution so that comparison studies can be performed to test for grid quality and determine the resolution required to resolve geologic structure and physical processes such as fluid flow and solute transport. The process of generating numerical grids with appropriate property distributions from geologic conceptual models is automated, thus making the entire process easy to implement with fewer user-induced errors. The series of grids of various resolutions are used to assess the level at which increasing resolution no longer influences the flow and solute transport results. Grid resolution is found to be a critical issue for ground water flow and solute transport. The resolution required in a particular instance is a function of the feature size of the model, the intrinsic properties of materials, the specific physics of the problem, and boundary conditions. The asymptotic nature of results related to flow and transport indicate that for a hydrologic model of the heterogeneous hydrostratigraphy under Yucca Mountain, a horizontal grid spacing of 600 m and vertical grid spacing of 40 m resolve the hydrostratigraphic model with sufficient precision to accurately model the hypothetical flow and solute transport to within 5% of the value that would be obtained with much higher resolution. 相似文献
18.
在大地电磁二维反演中,如果网格不均匀或者测点在网格中的位置不恰当,则将出现测点偏离现象,影响反演的精度.此外,由于正演精度和稳定性的需要,反演中网格间距不能差别太大,需要具有一定的光滑性.本文提出了几种自动构建测点中心网格的方法来消除测点偏离现象,并以正则化约束条件来保证网格的光滑程度.本文还提出了一套网格质量检验指标.算例对比研究表明,对于不同的测点集,最佳网格构建的方法可能不相同,但其中“最平模型约束动态配点法”、“最光模型约束动态配点法”拥有最强的适应能力,它们在各种情况下所构建的网格一般同时具备较光滑的网格间距变化和较小的测点偏离程度.此外,采用高质量（光滑性好、测点偏离程度低）网格在一定程度上不仅可提高反演结果的可靠性,还能提高反演的计算效率.本文的方法为获得高质量的大地电磁二维反演网格提供了简单有效的途径. 相似文献
19.
Eulerian–Lagrangian localized adjoint methods (ELLAMs) provide a general approach to the solution of advection-dominated advection–diffusion equations allowing large time steps while maintaining good accuracy. Moreover, the methods can treat systematically any type of boundary condition and are mass conservative. However, all ELLAMs developed so far suffer from non-physical oscillations and are usually implemented on structured grids. In this paper, we propose a finite volume ELLAM which incorporates a novel correction step rendering the method monotone while maintaining conservation of mass. The method has been implemented on fully unstructured meshes in two space dimensions. Numerical results demonstrate the applicability of the method for problems with highly non-uniform flow fields arising from heterogeneous porous media. 相似文献
20.
《Advances in water resources》2004,27(9):899-912
This paper describes work that extends to three dimensions the two-dimensional local-grid refinement method for block-centered finite-difference groundwater models of Mehl and Hill [Development and evaluation of a local grid refinement method for block-centered finite-difference groundwater models using shared nodes. Adv Water Resour 2002;25(5):497–511]. In this approach, the (parent) finite-difference grid is discretized more finely within a (child) sub-region. The grid refinement method sequentially solves each grid and uses specified flux (parent) and specified head (child) boundary conditions to couple the grids. Iteration achieves convergence between heads and fluxes of both grids. Of most concern is how to interpolate heads onto the boundary of the child grid such that the physics of the parent-grid flow is retained in three dimensions. We develop a new two-step, “cage-shell” interpolation method based on the solution of the flow equation on the boundary of the child between nodes shared with the parent grid. Error analysis using a test case indicates that the shared-node local grid refinement method with cage-shell boundary head interpolation is accurate and robust, and the resulting code is used to investigate three-dimensional local grid refinement of stream-aquifer interactions. Results reveal that (1) the parent and child grids interact to shift the true head and flux solution to a different solution where the heads and fluxes of both grids are in equilibrium, (2) the locally refined model provided a solution for both heads and fluxes in the region of the refinement that was more accurate than a model without refinement only if iterations are performed so that both heads and fluxes are in equilibrium, and (3) the accuracy of the coupling is limited by the parent-grid size—a coarse parent grid limits correct representation of the hydraulics in the feedback from the child grid. 相似文献