共查询到19条相似文献，搜索用时 250 毫秒
1.
双重介质模型在岩溶地下水流动系统模拟中的应用 总被引：1，自引：0，他引：1
文章采用不稳定层状裂隙水流模型,应用有限元-卷积结合法,模拟了岩溶地下水流动系统。相对于等效多孔介质模型,双重介质模型模拟的水位下降速度更小,达到稳定的时间更长。距抽水井较近的观测孔处,双重介质模型模拟的水位变化过程阶段性较明显,水位变化过程类似于承压含水层-弱透水层的释水过程,裂隙和孔隙基质分别相当于承压含水层和弱透水层。然后讨论了影响裂隙和孔隙基质间水流交换项的因素,分析了孔组抽水后的渗流场特征。 相似文献
2.
块状岩溶介质可视为双重孔隙介质, 其含水层为承压含水层, 水文地质模型概化为各向同性非均质。根据块状岩溶承压含水层地下水运动规律, 建立二维承压含水层非稳定流数学模型, 采用有限体积法对其求解, 并以深圳葵涌区为例对模型进行识别和验证。结果表明: 各时段的计算水位值与实测水位值非常接近, 所建立的数学模型正确可应用。 相似文献
3.
4.
本文将北方岩溶介质划分为两种基本类型:溶隙网络—溶孔型和溶隙网络—微裂隙型。概化出裂隙岩溶介质的双重介质水流模型。解决了下面几个问题:（1）认识北方岩溶介质所具有的双重介质性;（2）分别对孔隙和微裂隙所构成的岩块域和溶隙网络域建立数学模型,并用岩块和裂隙的水量交换项将两者耦合起来。模型中岩块并未被概化成某一特殊几何形态如球体或平行板,而是按其实际产出形态处理;（3）在岩块中也应用了二维水流模型。以往作者多用一维流方程,目的在于求其解析解,这有许多不足之处;（4）在裂隙和岩块域均采用有限单元法,将两个域的模型联立求数值解,克服了解析法或解析—数值法联合的局限性。 在给定参数情况下,本文进行了典型水文地质条件下的计算,证明了模型是能够较好地刻画溶隙网络—溶孔或微裂隙型介质的特征和水流特点的。同时,本文成功地应用了济南西南郊岩溶含水层的实际资料进行了模拟计算,对双重介质方法的实用进行了听的尝试。 相似文献
5.
焦作矿区双重介质裂隙网络渗流与渗流模型研究 总被引：3，自引：1，他引：2
焦作矿区矿井水害严重,位于煤层底板下的石炭系薄层灰岩中的第八、第二层,以及奥系灰岩岩溶含水层是矿井主要充水水源。裂隙是矿区岩溶水渗流的主要通道,裂隙介质的不均匀性、各向异性、不连续性等特征非常明显,且不同规模的裂隙对岩溶水渗流的贡献不同。根据裂隙介质对岩溶水渗流所起的作用,对那些控制矿区岩溶水渗流分布且起主导渗透作用的巨型和大型裂隙建立非连续介质裂隙网络三维渗流模型;对大型和巨型裂隙之间的裂隙岩块建立连续介质渗透张量模型,两模型耦合在一起组成双重介质裂隙网络三维渗流模型。该模型能较好地刻画矿区岩溶水渗流和突水特征 相似文献
6.
江夏法泗岩溶塌陷区冲孔桩施工引发岩溶塌陷的成因机理 总被引：2，自引：2，他引：0
依托武汉市岩溶地面塌陷区监测系统,开展江夏法泗冲孔桩施工过程裂隙岩溶含水层的水位、水压力和影响半径等关键参数的动态监测,研究施工过程岩溶含水层水压力的动态响应,探讨其诱发覆盖层突然塌陷的成因机理,结果表明:在第四系孔隙承压水含水层施工时,对该含水层的影响较大,但水位变幅相对较为平缓;在碳酸盐岩裂隙岩溶水含水层施工时,由于含水介质的不均一性,未揭穿岩溶通道时,对岩溶水的影响较小;当揭穿岩溶通道时,岩溶水水位可能产生突增,此时,诱发岩溶地面塌陷的风险加大。 相似文献
7.
张之淦 《水文地质工程地质》1979,(6)
岩溶水运动的性质远比孔隙水、裂隙水复杂。当前研究岩溶水运动的重要课题之一是通过研究岩溶含水介质的特征和岩溶水运动的性质建立起岩溶水运动的基本水文地质模型。近些年来许多区域性研究逐渐揭示出岩溶含水层具有宏观的双重介质的性质。通常岩溶水的水位变化十分剧烈,一般不能忽略它们的非稳定流特征。因此在一些天然条件下,岩溶水运动可以近似地看做宏观双重介质内的非稳定流运动。本文打算用老庄子泉动态资料表述一下宏观双重介质非稳定流水文地质模型的一些特征,并从这一模型出发进而讨论一下老庄子泉动态成因的一些问题。 相似文献
8.
山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究 总被引：11，自引：0，他引：11
在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上 ,对裂隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时 ,通常用等价多孔介质模型 ,但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙网络 ,导水系数大 ,地下水的实际平均流速比孔隙水大得多 ,但给水度和贮水系数小。当用等价多孔介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟 ,要同时求解水头方程和对流弥散方程 ,可采用MODFLOW和MT3D软件进行模拟。研究区裂隙岩溶水水头的数值计算表明 ,等效多孔介质模型水头的拟合误差能满足国标GB/T144 97- 93的要求。各时段地下水水量均衡计算的精度也满足要求。对流弥散方程的数值计算 ,由于Peclet数高达 95 .6 7,对流占绝对优势 ,可能存在数值弥散和数值振荡 ,因而采用多种方法进行了比较。对于同一问题 ,同时采用上游有限差分法 (UFDM) ,混合的欧拉拉格朗日方法 (特征线法MOC、改进特征线法MMOC和混合特征线法HMOC) ,总变异消减法(TVD)进行计算 ,并比较其结果。结果表明 ,混合特征线法 (HMOC)和总变异消减法 (TVD)比较适合于对流占优势的运移问题计算。由于渗透系数K和有效孔隙度θ对溶质运移结果的影响很大 ,? 相似文献
9.
溶蚀作用下古岩溶盆地系统中介质场演化模拟 总被引：2，自引：1，他引：1
为了研究古岩溶盆地在溶蚀作用下碳酸盐岩油气储层介质场的演化规律,基于地下水渗流理论和碳酸盐溶蚀动力理论,用数值模型模拟了均质盆地和非均质盆地2种岩溶盆地含水系统发育过程。岩溶盆地含水系统演化过程中总伴随着裂隙差异溶蚀,不管均质盆地还是非均质盆地,都越来越非均质化,强溶蚀带集中在潜水位及优势裂隙附近,在侵蚀基准面和构造裂隙及层理等大裂隙处会形成高孔隙率、渗透率的良好储层。受介质场非均质化反作用,岩溶高地区潜水位随系统演化不断下降,潜水面处的强烈溶蚀随水位不断下切而使高地区改造成竖直裂隙发育的厚储层;坡地区水平径流活跃,易于形成水平裂隙发育的储层;平原区水位相对稳定,溶蚀作用主要发生大裂隙和侵蚀基准面处,但在大裂隙网内经溶蚀可形成局部高孔隙率储层。对比两子模型发现,在岩溶盆地含水系统中,主导渗流场都要经历由局部流场向二级流场、二级流场向全局流场的转变,前者发生在坡地区内部,后者发生在坡地区与平原区之间;流畅的渗流场转换更有利于介质场发育,受大裂隙网导水作用,非均质盆地比均质盆地更快捷、更顺利地实现主导流场升级,溶蚀作用更强烈,3 ka后非均质盆地比均质盆地的孔隙率增幅大34%。 相似文献
10.
试验场地内现场调查发现固结致密砂岩层中存在一个倾角47°的倾斜裂隙含水层。为了研究倾角所产生的影响,通过建立倾斜裂隙微水试验的双重孔隙介质理论模型,利用理论和现场试验数据分析方法,讨论倾角对水文地质参数推估的影响。结果表明:低渗透条件下,倾角越大非振荡水位恢复速度越快;高渗透条件下,倾角越大振荡水位振幅越大。裂隙渗透系数越小,母岩贮水率越大,双重孔隙介质特征参数越大时,倾角上限越小,倾角影响越明显。根据该结论确定出模型适用于倾角大于30°的裂隙含水层。忽略倾角会导致裂隙渗透系数估值偏高,母岩贮水率估值偏低,双重孔隙介质特征参数估值偏低。 相似文献
11.
We derive a macroscopic model for single-phase, incompressible, viscous fluid flow in a porous medium with small cavities called vugs. We model the vuggy medium on the microscopic scale using Stokes equations within the vugular inclusions, Darcy's law within the porous rock, and a Beavers–Joseph–Saffman boundary condition on the interface between the two regions. We assume periodicity of the medium and obtain uniform energy estimates independent of the period. Through a two-scale homogenization limit as the period tends to zero, we obtain a macroscopic Darcy's law governing the medium on larger scales. We also develop some needed generalizations of the two-scale convergence theory needed for our bimodal medium, including a two-scale convergence result on the Darcy–Stokes interface. The macroscopic Darcy permeability is computable from the solution of a cell problem. An analytic solution to this problem in a simple geometry suggests that: (1) flow along vug channels is primarily Poiseuille with a small perturbation related to the Beavers–Joseph slip, and (2) flow that alternates from vug to matrix behaves as if the vugs have infinite permeability. 相似文献
12.
Many pieces of evidence have confirmed that the seepage in fine-grained soils can deviate from traditional Darcy's law that contributes to the consolidation theory. In this paper, a numerical model, referred to as consolidation with non-Darcian flow 2 (CNDF2), is developed for the 1D large strain consolidation of a saturated porous medium with the non-Darcian flow. The algorithm accounts for vertical strain, general constitutive relationships, the relative velocity of the fluid and solid phases, variable compressibility and permeability relation during consolidation, time-dependent loading, and unload/reload effects. Compared with the CS2 model proposed in previous studies, some parts of the CNDF2 model are modified in order to adapt to the non-Darcian flow, especially the equivalent series hydraulic conductivity and fluid flow. The verification examples of the CNDF2 model demonstrate excellent accuracy for both small strain and large strain consolidation. According to the applicable conditions of non-Darcian flow law, the CNDF2 model is most suitable for the fine-grained soil. The development of CNDF2 is first presented, followed by examples to analyze the influences of the non-Darcian flow on the consolidation behavior. 相似文献
13.
In this paper we are concerned with modeling single phase flow in a medium with known fractures. In particular we are interested
in the case in which the flow rate in the fractures is large enough to make it appropriate to use Forchheimer’s law for modeling
the flow in the fractures even though the flow in the surrounding domain is such that Darcy’s law is adequate. We describe
a model in which the fractures are treated as interfaces. We also consider the case of intersecting fractures and the case
of nonconforming meshes. 相似文献
14.
以页岩油气为代表的低品位油气资源勘探与开发不断取得重大突破,已经成为中国重要的接替性资源。开展页岩气渗流机理和数值模拟模型研究有助于实现页岩气藏开发的动态变化过程,为认识页岩气渗流规律、优化数学模型、产能评价和预测奠定技术基础。围绕页岩气多尺度多流态多重介质下的运移机理,系统地阐述了页岩气数值模拟模型的研究进展,页岩气数值模拟模型可分为等效连续介质模型、离散裂缝网络模型和混合模型,总结了这3类数值模拟模型的优缺点。等效连续介质模型原理简单,追求宏观尺度的等效,忽略了储层内部真实流态,适用于裂缝发育程度低的均质页岩气藏;离散裂缝网络模型准确反映复杂裂缝网络的渗流特征,可以描述高度离散裂缝的形态规律,适用于勘探程度高且裂缝高度发育的页岩气藏;混合模拟模型结合两者的优点,能够准确反映复杂裂缝网络和流体运移规律,在满足计算精度的同时又节约了大量的计算资源,随着计算处理能力的增强,混合模拟模型是今后的发展趋势。最后分析了页岩气藏数值模拟模型中存在的问题,并指出了发展方向。 相似文献
15.
In this paper, 3D steady‐state fluid flow in a porous medium with a large number of intersecting fractures is derived numerically by using collocation method. Fluid flow in the matrix and fractures is described by Darcy's law and Poiseuille's law, respectively. The recent theoretical development presented a general potential solution to model the steady‐state flow in fractured porous media under a far‐field condition. This solution is a hypersingular integral equation with pressure field in the fracture surfaces as the main unknown. The numerical procedure can resolve the problem for any form of fractures and also takes into account the interactions and the intersection between fractures. Once the pressure field and then the flux field in fractures have been determined, the pressure field in the porous matrix is computed completely. The basic problem of a single fracture is investigated, and a semi‐analytical solution is presented. Using the solution obtained for a single fracture, Mori‐Tanaka and self‐consistent schemes are employed for upscaling the effective permeability of 3D fractured porous media. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
16.
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. 相似文献
17.
18.
We consider a stationary flow of an incompressible non-Newtonian flow through a porous medium, induced by an injection velocity when inertial effects are negligible. At the pore scale, the governing equations are based on a nonlinear relation between the stress and the rate of deformation. In such a situation, the limit problem obtained when the pore size tends to zero, is called the homogenized problem that leads to the filtration law. This filtration law is given by a non-linear system coupling a local problem on a typical cell of the porous medium to a global problem at the scale of the whole porous medium. We propose, in this work, a numerical method to solve this homogenized problem and apply this method when the velocity dependent viscosity is given by the power law. Finally, we propose some numerical experiments to illustrate our approach. 相似文献
19.
高放废物地质处置过程中涉及的核素在围岩裂隙地下水中的迁移问题已引起广泛关注,数值模拟是研究核素粒子运移的重要方法。目前裂隙介质中渗流模型主要是等效连续介质模型、双重介质模型和离散裂隙网络模型。对于岩体尺度裂隙地下水的流动,离散裂隙网络模型能充分表现裂隙介质的各向异性、不连续性等特征。因此,针对裂隙介质准确概化及核素迁移模拟等难点,文章结合Monte Carlo随机生成裂隙方法、裂隙渗流有限元算法和高放射性核素衰变方程等方法,依据花岗岩深钻孔裂隙统计数据,采用离散裂隙网络模型对内蒙古阿拉善高放废物地质处置预选区展开了核素粒子迁移数值模拟研究,并讨论了实例预测分析结果。结果显示:针对设定的地质模型,核素粒子从中心运移到边界的迁移路径长度平均为1293.35 m,粒子运移到边界耗费的时间平均为1.70E+11 d。 相似文献