碳酸盐岩单裂隙渗流-溶蚀耦合模型及其参数敏感性分析

岩溶地下水系统是由碳酸盐岩裂隙含水介质演化形成的,系统初始的裂隙网络介质特征及边界条件决定了其演化过程。为揭示岩溶系统演化过程中裂隙介质特征和边界条件的影响程度,建立了裂隙溶蚀扩展的渗流-溶蚀耦合模型,并对不同边界条件下不同隙宽的单裂隙溶蚀扩展特征进行了模拟分析。结果表明:裂隙溶蚀扩展受水的侵蚀性(CO2分压)、水动力条件(水力梯度)、裂隙介质特征(裂隙初始隙宽)等综合作用影响,Ca2+的平衡浓度、水力梯度以及裂隙初始隙宽等参数的增加均能促进裂隙的快速扩展。在这些参数中,初始隙宽B0对岩溶发育的影响最为敏感,水力梯度J和Ca2+平衡浓度Ceq对岩溶发育具有相同的敏感性;此外,随着各参数值的不断增大,参数变化对岩溶发育的敏感程度越来越低。     

层流—紊流共存流场中岩溶裂隙网络演化过程的数值模拟方法

研究岩溶水系统的演化过程对许多资源与环境问题有着重要意义。岩溶地区水资源、油气资源的预测与开发、水土流失成因及防治等课题都与此密切相关。岩溶管道是地下水对裂隙的逐渐溶蚀扩展形成的,岩溶演化初期,所有裂隙宽度不大,流场整体成层流状态。随着溶蚀的进行,一部分裂隙优先扩大使其中的水流进入紊流状态,而另一部分裂隙中水流仍呈层流状态。文章提出一种数值方法,能够模拟层流-紊流共存流场的岩溶演化过程。利用蒙特卡洛(Monte-Carlo)方法模拟初始裂隙网络,通过非连续介质方法模拟裂隙网络中的渗流。裂隙扩展的速度通过岩石表面溶蚀速度经验公式计算,使用迭代方法求解层流-紊流共存条件下的水头非线性方程。构建了能够利用解析法求解的模型,把数值解和解析解的结果进行对比,验证了本研究的数值法及软件的可行性。      

高山峡谷区岩溶水示踪试验──以川西锦屏地区为例

示踪结果显示,大理岩是以溶蚀裂隙为主的含水介质,其岩溶水流迭较溶洞管道型含水层低一个数量级。在高山峡谷区,包气带深厚,示踪水流的“分流”现象十分普遍;隔水边界易于遭受破坏,不同单元之间越流水力联系较为广泛;且深部裂隙系统开启程度好,岩溶发育深度超过２０００～３０００ｍ,由于补给区与排泄区高差大,为深部区域水流系统的发育提供了介质和水动力条件。      

河间隐伏型岩溶裂隙含水系统演化的数值模拟

为揭示隐伏型岩溶含水系统在不同边界条件与介质特征下的发育状况,基于地下水渗流和碳酸盐岩溶蚀理论构建耦合裂隙网络渗流与碳酸盐岩溶蚀的岩溶裂隙网络溶蚀扩展模型,模拟再现了河间隐伏型岩溶含水系统平面上发育演化过程。结果显示:岩溶系统发育过程中流量变化呈现三个阶段:第Ⅰ阶段为缓慢层流阶段,流量增速为2.2×10 -6 mL·(s·a) -1,持续时间与初始水头差呈对数关系;第Ⅱ阶段为极快速紊流阶段,流量增速为5.5×10 2 mL·(s·a) -1,持续时间与初始水头差呈线性关系;第Ⅲ阶段为快速紊流阶段,流量增速为26 mL·(s·a) -1。优势裂隙的存在使得岩溶裂隙网络溶蚀演化加快,非均质系统进入极快速紊流阶段所用时间较均质系统缩短了53%,裂隙溶蚀扩展的速度加快了12.5%。      

兴文县溶洞特征、成因及演化研究

兴文县位于四川省东南部山地,县内拥有大量溶洞资源,结合兴文县溶洞资源调查与保护规划项目野外调查,对区内溶洞数量、规模、分布、形态、洞穴沉积物、地质背景、发育方向、分布高程等进行了统计和研究。该地区溶洞数量多,洞穴规模中等;溶洞空间分布与发育规模受地层与岩性控制;溶洞发育方向受构造影响明显,具体每个溶洞的发育方向主要受构造裂隙与岩石层面裂隙控制;区内中小型溶洞发育演化主要集中在沿节理裂隙的岩溶早期阶段和顺层岩溶作用的岩溶中期阶段,大型溶洞则表现为水平溶蚀与垂向溶蚀交替出现的岩溶中晚期阶段。     

裂隙岩溶含水系统溢流泉演化过程的数值模拟

本文利用非连续裂隙网络介质模型,耦合水动力条件和碳酸盐溶蚀动力条件,建立岩溶含水系统演化的数学模型,并用数值法进行模拟分析。在已建数值模型基础上,构造了两组实例模型(降水入渗补给的裂隙岩溶含水系统和河流补给的裂隙岩溶含水系统)来研究裂隙岩溶系统中溢流泉早期演化过程。模拟发现,裂隙岩溶含水系统演化过程中,裂隙不断被溶蚀拓宽,致使系统中潜水位不断下降,许多小泉被疏干;同时不同裂隙之间存在溶蚀差异,在潜水位及优势通道附近的裂隙溶蚀速度快,而其它位置则相对慢得多,差异溶蚀促使系统形成优势溶管,优势溶管逐渐袭夺水流,形成大泉。另外对比两模型发现,降水入渗为主要补给源的裂隙岩溶含水系统,泉管道从汇向源发育;而河流入渗为主要补给源的裂隙岩溶含水系统,泉管道则由源向汇发育。      

岩溶水系统演化中河间地块分水岭消失过程的数值模拟分析

本文通过数值模拟的方法研究了岩溶系统发展演化过程,并对其中分水岭消失现象进行了探讨。岩溶系统发育模拟中,假设裂隙段为光滑平板裂隙,裂隙中的水流处于平流状态,可以应用立方率描述水流在裂隙中的流动规律;碳酸盐岩裂隙的溶蚀过程可以通过溶蚀动力学的经验公式表示,即每个时间步长内,水流都处于稳定流状态。研究结果表明,岩溶系统溶蚀具有差异性,岩溶最为发育的部位位于潜水面附近,而且垂向的裂隙段比水平向的裂隙段发育更为强烈,深部岩溶系统发育缓慢,其裂隙张开度扩展率与潜水面附近的相差8200多倍。差异性溶蚀还造成潜水面附近形成优先水流通道,至分水岭从系统消失时,其强烈的汇水作用使几乎90%的流出量从此优先水流通道排出。分水岭消失时间受河间地块的范围和降雨量的影响,当河间地块的宽度增大至原来的5倍时,分水岭的存在时间几乎延长至原来的10倍,河间地块的范围愈大,分水岭越难以消失;降雨量增加对分水岭消失时间的影响表现出先增加后减小的趋势。     

山东峄城盆地岩溶发育的基本规律有与水源地形成的关系

山东省枣庄市峄城水源地是个典型的北方埋藏型岩溶水水源地。峄城盆地出露的岩溶地貌有干谷、溶洞、岩溶准平原、溶蚀裂隙、溶沟、溶槽、溶痕等，地下岩溶形态主要为溶洞、溶蚀裂隙和溶孔，岩溶的发育具有层控性、水平分带和垂直分带等特点和规律。该盆地岩溶水系统的含水介质为寒武系、奥陶系的灰岩、白云岩，以溶隙、溶孔为主加少量小型溶洞构成了一个统一的互相连通的地下水空隙网络系统，为地下水的运移、储存提供了通道和空间。岩溶水流在此空隙网络系统内具有渗流性质，流态以层流为主，形成一个基本统一的单斜扇形水动力网，并具有统一的水化学场，该盆地分间接补给、直接补给、汇集排泄和汇集四人输导功能区。并在汇集区形成了具有供水意义的大型岩溶水水源地。     

基于CFP的岩溶管道流溶质运移数值模拟研究

多重岩溶含水介质的复杂性导致岩溶地下水流动及溶质运移的数学模拟成为地下水研究难点之一。为了探讨岩溶多重含水介质中地下水流溶质运移特征,文章构建了管道流CFP水流模型和MT3DMS溶质运移三维耦合数值模型。在阐述管道流CFP和MT3DMS基本原理的基础上,通过建立水文地质概念模型算例（1个落水洞、4个直管道）,探讨岩溶管道水流及溶质运移规律,分析讨论不同水文地质参数对浓度穿透曲线的影响。研究结果表明:管道流CFP模型能够刻画岩溶管道与基岩裂隙水流交换特征,MT3DMS模型能够模拟穿透曲线的拖尾现象,符合实际岩溶区特征。随着水力梯度、管道直径及管道渗透系数增大,孔隙度减小,浓度曲线峰值越大,峰值到达时间越快,浓度穿透曲线越对称。得出结论:耦合CFP水流模型和MT3DMS溶质运移模型能够刻画岩溶管道流溶质运移规律,为研究岩溶复杂介质污染物运移特征提供一种思路和途径。     

基于CFP的岩溶管道流数值模拟研究——以桂林寨底地下河子系统为例

为了描述岩溶多重含水介质存在非达西流问题,研究管道流CFP模型在岩溶区数值模拟中的应用。在分析管道流CFP基本原理的基础上,通过建立概念模型算例(5个管道,上下临界雷诺数分别为2 000和4 000),探讨不同管道水文地质参数对出口流量及水流状态的影响,并对参数敏感性进行分析,最后将CFPM1模型应用于广西寨底岩溶地下河系统,探讨在实际应用中的管道流模型适用性。概念模型结果表明,参数敏感度从大到小依次为:管道直径、水力梯度、管道渗透系数、弯曲度及粗糙度,流量与雷诺数随管道参数增大而增大,水流状态从层流变为紊流。寨底岩溶管道流模型中3个观测井及地下河出口实测流量表明CFPM1能较好地模拟出岩溶地下水水位变化趋势。得出结论:CFPM1管道流模型允许岩溶管道与基岩含水层进行水流交换,能够较好地刻画岩溶区管道介质非达西流特征,但管道参数获取较为困难,且不能刻画岩溶管道形态变化。该方法具有一定通用性,可为实际应用提供指导。     

Simulating flow in karst aquifers at laboratory and sub-regional scales using MODFLOW-CFP

Groundwater flow in a well-developed karst aquifer dominantly occurs through bedding planes, fractures, conduits, and caves created by and/or enlarged by dissolution. Conventional groundwater modeling methods assume that groundwater flow is described by Darcian principles where primary porosity (i.e. matrix porosity) and laminar flow are dominant. However, in well-developed karst aquifers, the assumption of Darcian flow can be questionable. While Darcian flow generally occurs in the matrix portion of the karst aquifer, flow through conduits can be non-laminar where the relation between specific discharge and hydraulic gradient is non-linear. MODFLOW-CFP is a relatively new modeling program that accounts for non-laminar and laminar flow in pipes, like karst caves, within an aquifer. In this study, results from MODFLOW-CFP are compared to those from MODFLOW-2000/2005, a numerical code based on Darcy’s law, to evaluate the accuracy that CFP can achieve when modeling flows in karst aquifers at laboratory and sub-regional (Woodville Karst Plain, Florida, USA) scales. In comparison with laboratory experiments, simulation results by MODFLOW-CFP are more accurate than MODFLOW 2005. At the sub-regional scale, MODFLOW-CFP was more accurate than MODFLOW-2000 for simulating field measurements of peak flow at one spring and total discharges at two springs for an observed storm event.     

西南某岩溶区地下水系统示踪试验与解析

以西南某典型岩溶区为例,解析示踪试验在岩溶管道连通性以及获取水文地质参数中的应用。选择落水洞为投放点,分别从落水洞西侧和东侧寻找地下水出露点作为接收点,判别落水洞地下径流的实际去向以及落水洞与接收点之间的水力联系。结合Qtracer2软件对示踪试验成果进行定量解析,确定示踪剂回收率、地下水平均流速、最快流速,估算出岩溶管道结构特征和水文地质参数。结果表明:落水洞与接收点JS01、JS03之间不存在直接水力联系;落水洞与接收点JS02存在水力联系且岩溶管道极为发育,含水介质不均匀,地下水运移路径较为通畅,为典型的紊流流态;落水洞地下径流的主要方向是由西向东,但在丰水期雨量较大期间,接收点JS04能够接收到荧光素钠,说明丰水期水位上涨后两者间会有水力联系,导致部分水量向落水洞西侧排泄。     

Relative importance and chemical effects of diffuse and focused recharge in an eogenetic karst aquifer: an example from the unconfined upper Floridan aquifer, USA

Karst aquifer studies often focus on allogenic water inputs and large conduit flow. However, diffuse recharge can be significant, particularly in unconfined eogenetic karst aquifers that retain high matrix permeability. This study examines an unconfined region of the upper Floridan aquifer (USA) that hosts a sinking stream, its resurgence, and a large conduit system. Daily diffuse recharge was approximated using a water-budget method and ranged from 17% of precipitation during a low precipitation year to >53% during the highest precipitation year, illustrating the highly variable nature of diffuse recharge in this region. The total allogenic input via the sinking stream over the 5 years of the study was significantly larger than the volume of diffuse recharge. However, only about 2% of the allogenic recharge flows from the conduit into the surrounding aquifer. That flow is restricted to storm events when hydraulic heads in the conduits exceed those in the surrounding aquifer. The estimated volume of dissolution is similar for allogenic recharge and diffuse recharge to the unconfined region surrounding the conduits, but dissolution from the diffuse recharge is distributed over a larger area than dissolution from allogenic recharge. These results exemplify how recharge type impacts flow and water–rock interactions in eogenetic karst aquifers.     

Mise-a-la-masse method,tracer tests and cave detection analysis of underground-river karst-conduit distribution and structure: a case study of the Dafengdong conduits in Guizhou Province,China

In karst areas, groundwater movement is dominated by conduit flow; thus, understanding the distribution and structure of karst conduits has great significance for water research, groundwater protection and engineering construction. With the Dafengdong (Guizhou Province, China) underground river karst as an example, a mise-a-la-masse method, tracer tests and cave detection were conducted to study the distribution and structure of karst conduits. Combined with information on the geological and hydrogeological conditions, the geological factors that form karst conduits were determined. It can be concluded that: (1) Under the influence of faults, karst conduits usually develop towards or along faults; (2) multiple karst conduits form easily on both sides of tectonic fracture zones; (3) both lithology and geological structure affect the formation of blue holes, and when the directions of karst conduits change, blue holes easily form in weak parts of the structure; (4) at springs where two aquifers intersect, with strong chemical dissolution of the lower karst aquifer and mechanical collapse of the upper aquifer, blue holes also form easily; and (5) integrated mise-a-la-masse method, tracer tests and cave detection can accurately discern the distribution and structure of karst conduits. Geological factors can be used to preliminarily delineate the distribution and structure of karst conduits in similar areas based on hydrogeological conditions. Such methods hold great significance for groundwater extraction and protection and engineering construction in karst areas.

     

Application of MODFLOW for simulating groundwater flow in the Trifilia karst aquifer, Greece

Karst aquifers can have a complex flow as a result of the formation of large conduits from dissolution features. As a result, a three-dimensional finite-difference groundwater flow model (equivalent porous media) may not apply as the dual porosity nature of karst features and the effects of turbulent flow cannot be directly simulated. Statistical analysis of karst hydrographs of the Trifilia aquifer in Greece showed the existence of a slightly karstified mass with high primary porosity that regulates the flow. An equivalent porous media model was developed to simulate the Trifilia karst aquifer using MODFLOW. Steady state and transient state calibration gave encouraging results for the equivalent porous media approach, which does not consider pipe flow or turbulence. Detailed hydrogeological research conducted in the area helped define the aquifer hydraulic conductivity zones and extent; and flux to/from the aquifer. Only hydraulic conductivity and specific yield were adjusted during calibration, as the flux to/from the system was considered known and applied as boundary conditions. Small mean absolute and RMS piezometric head error of the model under both steady and transient state conditions were achieved.     

Applying a modified conduit flow process to understand conduit-matrix exchange of a karst aquifer

Due to the high heterogeneity and complexity of water flow movement for multiple karst water-bearing mediums,the evaluation,effective development,and utilization of karst water resources are significantly limited.Matrix flow is usually laminar,whereas conduit flow is usually turbulent.The driving mechanisms of water exchange that occur between the karst conduit and its adjacent matrix are not well understood.This paper investigates the hydrodynamic characteristics and the mechanism of flow exchange in dual water-bearing mediums(conduit and matrix)of karst aquifers through laboratory experimentation and numerical simulation.A karst aquifer consisting of a matrix network and a conduit was proposed,and the relationship between the water exchange flux and hydraulic head differences generated from the laboratory experiments was analyzed.Two modes of experimental tests were performed with different fixed water level boundaries in the laboratory karst aquifer.The results indicate that the water exchange capacity was proportional to the square root of hydraulic head differences.The linear exchange term in the conduit flow process(CFP)source program was modified according to experimental results.The modified CFP and the original CFP model experimental data results were compared,and it was found that the modified CFP model had better fitting effects.These results showed that the water exchange mechanism between conduit and matrix is very important for solid-liquid interface reaction,water resource evaluation,and understanding of karst hydrodynamic behavior.     

Comparison of a karst groundwater model with and without discrete conduit flow

Karst aquifers exhibit a dual flow system characterized by interacting conduit and matrix domains. This study evaluated the coupled continuum pipe-flow framework for modeling karst groundwater flow in the Madison aquifer of western South Dakota (USA). Coupled conduit and matrix flow was simulated within a regional finite-difference model over a 10-year transient period. An existing equivalent porous medium (EPM) model was modified to include major conduit networks whose locations were constrained by dye-tracing data and environmental tracer analysis. Model calibration data included measured hydraulic heads at observation wells and estimates of discharge at four karst springs. Relative to the EPM model, the match to observation well hydraulic heads was substantially improved with the addition of conduits. The inclusion of conduit flow allowed for a simpler hydraulic conductivity distribution in the matrix continuum. Two of the high-conductivity zones in the EPM model, which were required to indirectly simulate the effects of conduits, were eliminated from the new model. This work demonstrates the utility of the coupled continuum pipe-flow method and illustrates how karst aquifer model parameterization is dependent on the physical processes that are simulated.     

The comparison of hydrodynamic characteristics of karst aquifers: application on two karst formations in Zagros (Asmari and Ilam-Sarvak), southwest Iran

In order to study the function, hydrodynamic behavior, and hydraulic properties of the karst aquifers in Izeh, southwest Iran, time series analysis was applied to the precipitation, spring discharge, and piezometric head data of two representative karst systems of Zagros (Ilam-Sarvak and Asmari Formations). Time series analysis was applied to two karst aquifers, those of Asmari and Ilam-Sarvak Formations. The daily precipitations of anticlines were estimated based on the precipitation–elevation function which was applied on digital elevation model (DEM) of the area. The mean estimated daily precipitations were considered in bivariate time series analysis as input data of each karst system. The total length of time series was about 2.7 years, extending from May 2007 to December 2009. During the research, several one-parameter probe data loggers were installed, which daily measure the water surfaces in karst aquifers. Time series analysis was applied for describing Izeh karst aquifers with a focus on both univariate (autocorrelation and spectral analysis) and bivariate (cross-correlation, gain function, and coherency function) methods. The results show that Asmari karst aquifer in Kamarderaz Anticline has a large storage capacity. Because of lacking a well-organized karst network, in the Asmari karst aquifer, baseflow dominates with low contribution of quick-flow. In the Ilam-Sarvak karst aquifer (Shavish and Tanush Anticlines), the karstification occurred in fractures and small diameter conduits, which caused to quick-flow between dense limestone. The Ilam-Sarvak karst aquifer could be regarded as a transition between two extreme types of karst, e.g., highly karstified system and in the opposite, extremely diffused one. The analysis of well hydrograph in Ilam-Sarvak karst aquifer shows that the karst aquifer has a low storage capacity. Unlike Asmari karst aquifer, the fractures and small diameter conduits in Ilam-Sarvak karst aquifer are more enhanced, producing a better developed karst network. Contrary to the typical karst systems, however, diffuse flow and conduit flow coexist in the Asmari Formation.