共查询到20条相似文献，搜索用时 31 毫秒
1.
随机水文地质学方法,较传统的确定性方法而言,是解决非均质含水层中水流和溶质运移问题的一种更为合理的手段。据以往研究,假设渗透系数场遵循对数正态分布,利用直接傅立叶变换方法来生成渗透系数随机场。应用基于随机理论的蒙特卡罗方法,来研究渗透系数的空间变异性对污染物运移结果的影响。实例研究表明,污染物在含水层中运移过程中污染羽的展布范围(二阶矩)随着渗透系数空间变异方差的增大而扩大,而污染羽在空间上的质心位置(一阶矩)基本不受方差的影响,仅取决于渗透系数随机场的均值大小。另外还分析了污染羽在各点的浓度变化方差和变异系数分别随渗透系数变异方差的变化状况。 相似文献
2.
Wang Mingyu The University of Arizona Tucson Arizona USA Department of Water Resources Environmental Engineering China University of Geosciences Beijing Chen Jinsong Wan Li Department of Water Resources Environmental Engineering 《中国地质大学学报(英文版)》2001,12(3)
ＩＮＴＲＯＤＵＣＴＩＯＮＧｒｏｕｎｄｗａｔｅｒｏｒｆｌｕｉｄｆｌｏｗｍｏｄｅｌｉｎｇｉｎｆｒａｃｔｕｒｅｄｒｏｃｋｓｉｓａｃｏｍｐｌｉｃａｔｅｄｔｈｅｏｒｅｔｉｃａｌａｎｄａｐｐｌｉｅｄｔｏｐｉｃ.Ｂｏｔｈｔｈｅｏｒｅｔｉｃａｌｌｙａｎｄｏｐｅｒａｔｉｏｎａｌｌｙ ,ｉｔｉｓｉｍｐｏｒｔａｎｔｉｎｍａｎｙｆｉｅｌｄｓｓｕｃｈａｓｇｅｏｌｏｇｉｃａｌａｎｄｈｙｄｒｏｇｅｏｌｏｇｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ ,ｅｎｖｉｒｏｎｍｅｎｔａｌｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｅｔｒｏｌｅｕｍｅｎｇｉｎｅｅｒｉｎ… 相似文献
3.
在地下水污染模拟预报中,弥散参数是很难确定的一个模型参数。因实验室小尺度弥散规律一般不能用于大尺度弥散过程,而野外示踪试验却耗资大、周期长,限制了其实用性。文中利用随机数值模拟手段、基于随机理论的蒙特卡罗方法及序贯高斯模拟技术来生成渗透系数随机场,并研究渗透系数对数场的方差、相关长度以及变异函数类型在不同尺度上对纵向弥散度的影响,进而建立纵向弥散度与随机分布渗透系数场的方差和相关长度的统计定量关系,并与Gelhar理论计算结果进行比较。数值模拟结果表明,经过一定迁移距离后纵向弥散度与随机分布渗透系数对数场的方差和相关长度具有良好的线性统计关系,与Gelhar理论公式表达的关系类型类似。但对于较大的方差,纵向弥散度模拟结果明显大于Gelhar理论计算值,而对于较大相关长度在迁移距离不很大时,纵向弥散度模拟结果明显小于Gelhar理论计算值。本研究可为野外大尺度地下水污染预报模型中水动力弥散参数的确定提供方法借鉴。 相似文献
4.
Numerical evaluation of multicomponent cation exchange reactive transport in physically and geochemically heterogeneous porous media 总被引：1，自引：0，他引：1
Experimental evidence and stochastic studies strongly show that the transport of reactive solutes in porous media is significantly
influenced by heterogeneities in hydraulic conductivity, porosity, and sorption parameters. In this paper, we present Monte
Carlo numerical simulations of multicomponent reactive transport involving competitive cation exchange reactions in a two-dimensional
vertical physically and geochemically heterogeneous medium. Log hydraulic conductivity, log K, and log cation exchange capacity (log CEC) are assumed to be random Gaussian functions with spherical semivariograms. Random
realizations of log K and log CEC are used as input data for the numerical simulation of multicomponent reactive transport with CORE2D, a general purpose reactive transport code. Longitudinal features of the fronts of reactive and conservative species are
computed from the temporal and spatial moments of depth-averaged concentrations. Monte Carlo simulations show that: (1) the
displacement of reactive fronts increases with increasing variance of log K, while it decreases with the variance of log CEC; (2) second-order spatial moments increase with increasing variances of
log K and log CEC; (3) uncertainties in the mean arrival time are largest (smallest) for negatively (positively) correlated log
K and Log CEC; (4) cations undergoing competitive cation exchange exhibit different apparent velocities and retardation factors
due to both physical and geochemical heterogeneities; and (5) the correlation between log K and log CEC affects significantly apparent cation retardation factors in heterogeneous aquifers. 相似文献
5.
This work presents the application of a Monte Carlo simulation method to perform an statistical analysis of transient variably saturated flow in an hypothetical random porous media. For each realization of the stochastic soil parameters entering as coefficients in Richards' flow equation, the pressure head and the flow field are computed using a mixed finite element procedure for the spatial discretization combined with a backward Euler and a modified Picard iteration in time. The hybridization of the mixed method provides a novel way for evaluating hydraulic conductivity on interelement boundaries. The proposed methodology can handle both large variability and fractal structure in the hydraulic parameters. The saturated conductivity K
s and the shape parameter vg in the van Genuchten model are treated as stochastic fractal functions known as fractional Brownian motion (fBm) or fractional Gaussian noise (fGn). The statistical moments of the pressure head, water content, and flow components are obtained by averaging realizations of the fractal parameters in Monte Carlo fashion. A numerical example showing the application of the proposed methodology to characterize groundwater flow in highly heterogeneous soils is presented. 相似文献
6.
潜水水流的动态随机模拟是一个复杂而难解决的问题.通过建立二维潜水非稳定流模拟的摄动随机有限元模型, 把控制方程的主要参数渗透系数和给水度随机变量、及源汇项和边界条件看作随机变量.在充分考虑4种随机因素的条件下, 推导出求解潜水二维非稳定流均值和方差的9个方程; 重点介绍了不同方程数值离散的特殊处理方法.通过设定理想例子对模拟结果进行了分析, 表明随机变量中边界条件值方差、渗透系数方差变化对水头方差变化的影响很小, 给水度方差的变化对水头方差的变化影响很大.本模型考虑因素全面, 对一般的潜水非稳定流随机模拟都可应用.本研究给出了边界、渗透系数、给水度的随机因素对潜水动态模拟的影响, 丰富和补充了地下水运动的随机理论. 相似文献
7.
离散裂隙渗流方法与裂隙化渗透介质建模 总被引：4，自引：1，他引：4
流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。 相似文献
8.
A new perturbation technique and finite-element method, which incorporates an undetermined-coefficients approach, is proposed to conduct stochastic simulation for problems of groundwater flow. Formulas of the mean value and variance of groundwater head are derived to exclude the process of computing the first and second-order derivatives of groundwater head with respect to stochastic variables (coefficient of transmissibility, storage coefficient, etc.) in the calculation process of common finite element methods. For ordinary stochastic groundwater problems, the statistical properties of groundwater head can be easily obtained by using the proposed method, which is especially efficient for the stochastic problems that have fewer input stochastic variables. A stochastic numerical simulation of a two-dimensional confined-groundwater flow problem has been conducted to validate effectiveness and determine the limitations of the proposed method. The results illustrate that the accuracy and efficiency of the method proposed are satisfactory, and that the computing time is shorter when there is a small number of input stochastic variables. 相似文献
9.
In this paper, a new enrichment scheme is proposed to model fractures and other conduits in porous media flow problems. Inserting
this scheme into a partition of unity based method results in a new numerical method that does not require the mesh to honor
the specific geometry of these subsurface features. The new scheme involves a specially designed integration procedure and
enrichment functions, which can capture effects of local heterogeneity introduced by subsurface features on the pressure solution.
The new method is also capable of modeling fractures with low as well as high conductivity. Another feature of the proposed
scheme is that, even though two enrichment functions are used to model the permeability change at the two rock/fracture interfaces
of a fracture, only one element partition is made for numerical integration. To demonstrate the accuracy and effectiveness
of the proposed approach, production problems for wells that were stimulated or completed by longitudinal fracture, transverse
fractures, and perforations are studied. 相似文献
10.
The problem of multiphase phase flow in heterogeneous subsurface porous media is one involving many uncertainties. In particular, the permeability of the medium is an important aspect of the model that is inherently uncertain. Properly quantifying these uncertainties is essential in order to make reliable probabilistic-based predictions and future decisions. In this work, a measure-theoretic framework is employed to quantify uncertainties in a two-phase subsurface flow model in high-contrast media. Given uncertain saturation data from observation wells, the stochastic inverse problem is solved numerically in order to obtain a probability measure on the space of unknown permeability parameters characterizing the two-phase flow. As solving the stochastic inverse problem requires a number of forward model solves, we also incorporate the use of a conservative version of the generalized multiscale finite element method for added efficiency. The parameter-space probability measure is used in order to make predictions of saturation values where measurements are not available, and to validate the effectiveness of the proposed approach in the context of fine and coarse model solves. A number of numerical examples are offered to illustrate the measure-theoretic methodology for solving the stochastic inverse problem using both fine and coarse solution schemes. 相似文献
11.
多孔介质渗透系数的空间尺度效应研究进展 总被引：6，自引：1，他引：5
多孔介质渗透系数的空间尺度问题是一个与地下流体运动和溶质运移的数值模拟密切相关的应用性课题,广泛的应用需求和新的计算方法使其成为近年的热门课题之一。它涉及到相互联系的两个方面：①非均质介质场渗透系数空间尺度行为的分析与模拟;②将局部测量尺度下的试验参数转化为数值模拟网格尺度下的参数输入值的升尺度（ｕｐｓｃａｌｉｎｇ）方法和计算模型。首先介绍了该课题在概念上的拓展及其物理含义,进而以方法为主线,对这一领域具有代表性的研究成果进行了分类和评述,讨论了该课题的研究对地下水流和溶质运移的模拟分析乃至整个多孔介质流体运动研究的意义。 相似文献
12.
基于择多算子的随机搜索法建立数字岩心的新技术 总被引：3，自引：0，他引：3
介绍了用于刻画岩心孔隙空间特征的两点概率函数、线性路径函数,提出了建立数字岩心的新方法--基于择多算子的随机搜索法,系统阐述了该方法的建模理论及基于择多算子的孔隙点、骨架点的随机搜索策略,指出搜索孔隙点、骨架点的择多算子的适宜取值区间分别为[9,12]和[4,7]。研究表明,该法进行数字岩心建模不仅可行而且效果良好。它不受岩性地限制;所建数字岩心的孔隙形态及空间分布特征与真实岩心孔隙相似;尽管岩心含有一定数量的孤立孔隙和岩石颗粒,但其规模很小;岩心具有良好的孔隙连通性、各向同性和渗透性,故可用来研究多孔介质中流体的分布和渗流等问题;由于该法建立的数字岩心的连通性随孔隙度的减小而变差,故该法更适于建立大孔隙度岩心。 相似文献
13.
14.
矿物热分解动力学的研究方法探讨 总被引：5，自引：0，他引：5
:矿物分解动力学的研究可以大大提高我们对地质作用过程中矿物生长、元素和流体迁移等过程的了解。本文在介绍固相反应动力学原理的基础上 ,讨论了常压下研究矿物分解动力学的热重分析法和电导率法 ,指出精确确定高压下的反应分数 ,是研究高压下矿物热分解动力学的关键。提出了高压下确定反应分数的方法。 相似文献
15.
孔隙液体的可压缩性对非饱和孔隙介质弹塑性分析的影响 总被引：2，自引：0，他引：2
首先在多相孔隙介质的混合物理论框架内探讨了如何利用固-液两相孔隙介质模型描述介质的非饱和问题。其次,讨论了一个实用的描述孔隙液体可压缩性的方法,并建立了小变形条件下固-液两相孔隙介质的弹塑性本构方程。最后,采用俞茂宏统一强度理论,针对一个具体的结构例子,分析了孔隙液体的可压缩性对结构弹塑性分析的影响。 相似文献
16.
地下水污染事件难以被直接观测,其发生环境通常也具有多种不确定性,风险评价可以量化地下水污染事件的危害性。利用随机配点模型和多项式抽样技术,建立了高效的风险评价模型;考虑了渗透系数、孔隙率和弥散度等多种不确定性因子,探讨了复杂条件下的风险函数的分布特征。研究结果表明:基于随机配点法的风险评价模型避免了多次重复求解对流弥散方程,通过计算成本低廉的随机配点技术得到浓度随机场的拉格朗日多项式,进行多项式抽样获取浓度样本并得到风险函数;与传统的解析算法相比,该方法无需对输入参数和浓度的分布形态做出假设;与传统的蒙特卡罗(Monte Carlo)算法相比,该模型具有明显的效率优势和优越的收敛速度;输入参数的分布类型对风险分布产生显著影响。 相似文献
17.
溶蚀岩体随机结构模型建立及其在岩体渗漏评价中的应用 总被引：5，自引：1，他引：4
针对中小规模裂隙和溶蚀洞穴造成的溶蚀岩体渗漏 ,本文提出了溶蚀岩体渗漏评价的一种新思路 :首先对现场开挖所揭露的有限的裂隙和溶蚀洞穴观测资料进行统计分析 ,利用统计分析结果建立溶蚀岩体随机结构模型 ;然后根据随机结构模型所反映出的裂隙和溶蚀洞穴有关参数计算渗透参数 ;最后结合区域水力学条件计算岩体渗漏量。多次模拟后可以得到渗漏量的分布规律。其中 ,建立合理的随机结构模型是关系本方法成败的关键 ,文中给出了建立溶蚀岩体随机结构模型的步骤 ,提供了随机模拟的编程思路和程序框图 ,并结合具体的工程实例对该方法的应用前景进行了讨论 相似文献
18.
19.
This study experimentally investigates the effect of particle size, particle concentration and flow velocity on the migration of suspended particles of size 1.02–47 μm in porous media. The results show that at the same flow velocity, the peak values of the breakthrough curves decrease and corresponding pore volumes increase slightly with increasing particles size. The migration velocity of smaller suspended particles is even greater than water flow velocity, which is attributed to the size exclusion effect. With increase of the injected particle concentration, the deposition coefficients of small single particles increase at first and then tend to a steady state or even decrease slightly, explained by the maximum retention concentration. The dispersivity of small particles decreases with increasing velocity. However, at a high flow velocity, the hydrodynamic dispersivity becomes increasingly dominant with the increase of particle size. The deposition coefficients for large-sized particles are higher than those for small-sized particles, which is attributed to considerable mass removal due to straining. An analytical solution, considering the release effect of sorbed particles, is developed to account for the one-dimensional flow and dispersive effect using a source function method, and then three transport parameters—dispersivity, deposition coefficient and release coefficient—are fitted using the experimental results. Finally, suspended-particle migration is predicted by the proposed model for short-time constant-concentration injection and repeated three-pulse injection. Overall, particle size has a significant effect on the seepage migration parameters of suspended particles in porous media such as the particle velocity, dispersivity and deposition coefficient. 相似文献
20.
In this paper, a stochastic collocation-based Kalman filter (SCKF) is developed to estimate the hydraulic conductivity from
direct and indirect measurements. It combines the advantages of the ensemble Kalman filter (EnKF) for dynamic data assimilation
and the polynomial chaos expansion (PCE) for efficient uncertainty quantification. In this approach, the random log hydraulic
conductivity field is first parameterized by the Karhunen–Loeve (KL) expansion and the hydraulic pressure is expressed by
the PCE. The coefficients of PCE are solved with a collocation technique. Realizations are constructed by choosing collocation
point sets in the random space. The stochastic collocation method is non-intrusive in that such realizations are solved forward
in time via an existing deterministic solver independently as in the Monte Carlo method. The needed entries of the state covariance
matrix are approximated with the coefficients of PCE, which can be recovered from the collocation results. The system states
are updated by updating the PCE coefficients. A 2D heterogeneous flow example is used to demonstrate the applicability of
the SCKF with respect to different factors, such as initial guess, variance, correlation length, and the number of observations.
The results are compared with those from the EnKF method. It is shown that the SCKF is computationally more efficient than
the EnKF under certain conditions. Each approach has its own advantages and limitations. The performance of the SCKF decreases
with larger variance, smaller correlation ratio, and fewer observations. Hence, the choice between the two methods is problem
dependent. As a non-intrusive method, the SCKF can be easily extended to multiphase flow problems. 相似文献