Experimental and computational validation and verification of the Stokes–Darcy and continuum pipe flow models for karst aquifers with dual porosity structure

In our previous study, we developed the Stokes–Darcy (SD) model was developed for flow in a karst aquifer with a conduit bedded in matrix, and the Beavers–Joseph (BJ) condition was used to describe the matrix–conduit interface. We also studied the mathematical well‐posedness of a coupled continuum pipe flow (CCPF) model as well as convergence rates of its finite element approximation. In this study, to compare the SD model with the CCPF model, we used numerical analyses to validate finite element discretisation methods for the two models. Using computational experiments, simulation codes implementing the finite element discretisations are then verified. Further model validation studies are based on the results of laboratory experiments. Comparing the results of computer simulations and experiments, we concluded that the SD model with the Beavers–Joseph interface condition is a valid model for conduit–matrix systems. On the other hand, the CCPF model with the value of the exchange parameter chosen within the range suggested in the literature perhaps does not result in good agreement with experimental observations. We then examined the sensitivity of the CCPF model with respect to the exchange parameter, concluding that, as has previously been noted, the model is highly sensitive for small values of the exchange parameter. However, for larger values, the model becomes less sensitive and, more important, also produces results that are in better agreement with experimental observations. This suggests that the CCPF model may also produce accurate simulation results, if one chooses larger values of the exchange parameter than those suggested in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

A new solution for a partially penetrating constant‐rate pumping well with a finite‐thickness skin

A mathematical model describing the constant pumping is developed for a partially penetrating well in a heterogeneous aquifer system. The Laplace‐domain solution for the model is derived by applying the Laplace transforms with respect to time and the finite Fourier cosine transforms with respect to vertical co‐ordinates. This solution is used to produce the curves of dimensionless drawdown versus dimensionless time to investigate the influences of the patch zone and well partial penetration on the drawdown distributions. The results show that the dimensionless drawdown depends on the hydraulic properties of the patch and formation zones. The effect of a partially penetrating well on the drawdown with a negative patch zone is larger than that with a positive patch zone. For a single‐zone aquifer case, neglecting the effect of a well radius will give significant error in estimating dimensionless drawdown, especially when dimensionless distance is small. The dimensionless drawdown curves for cases with and without considering the well radius approach the Hantush equation (Advances in Hydroscience. Academic Press: New York, 1964) at large time and/or large distance away from a test well. Copyright © 2007 John Wiley & Sons, Ltd.     

Heterogeneous aquitard properties in sedimentary successions in the Apennine chain: case studies in southern Italy

The aim of this study was to analyse the hydrogeologic role of siliciclastic deposits and marly‐clayey‐calcareous successions within the carbonate Apennine chain (southern Italy). The study was carried out along the northern part of the Matese carbonate massif through (1) the hydraulic characterization of siliciclastic rocks in a test site, by means of Lugeon tests, and (2) the identification of the groundwater flow system discharging at an important spring located within a marly‐clayey‐calcareous succession in a second test site, by means of isotopic investigations. The results showed that the investigated siliciclastic deposits and marly‐clayey‐calcareous successions may allow significant groundwater discharge from carbonate aquifers. Thus, they do not everywhere behave as aquitard, contrary to the previous model. Instead, groundwater flows through the upper part of these successions, where stress release fracturing enhanced rock permeability in the near‐surface bedrock. Thus, these successions may locally be a new groundwater source within the southern Apennine chain. Copyright © 2009 John Wiley & Sons, Ltd.     

南海西部盐水层CO2埋存潜力评估

南海西部东方1-1、崖城13-1和乐东等气田具有CO2含量高的特点,每年从生产出的天然气中分离出大量的CO2。为了减少CO2在大气中的排放,考虑在莺歌海地区选择合适的盐水层埋存体,拟进行CO2地质埋存示范工程。根据CO2在盐水层中的各种埋存机理,并考虑盐水层构造特征对CO2运移分布的影响,提出了一种较为准确的CO2埋存潜力评估方法,并利用此方法对筛选出的5个备选盐水层进行了评估。评估结果表明,各盐水层埋存潜力巨大,都远远大于示范工程期限内预计的CO2埋存总量,并可将这些盐水层作为将来海南省及临近广东省人为CO2的埋存场所。南海西部东方1-1、崖城13-1和乐东等气田具有CO2含量高的特点,每年从生产出的天然气中分离出大量的CO2。为了减少CO2在大气中的排放,考虑在莺歌海地区选择合适的盐水层埋存体,拟进行CO2地质埋存示范工程。根据CO2在盐水层中的各种埋存机理,并考虑盐水层构造特征对CO2运移分布的影响,提出了一种较为准确的CO2埋存潜力评估方法,并利用此方法对筛选出的5个备选盐水层进行了评估。评估结果表明,各盐水层埋存潜力巨大,都远远大于示范工程期限内预计的CO2埋存总量,并可将这些盐水层作为将来海南省及临近广东省人为CO2的埋存场所。     

深层水漏斗区开采量组成变化特征与机制

为阐明强烈开采条件下华北平原深层水的开采水量来源及组成问题,以德州地区深层水漏斗分布区为重点研究区,通过基于大量钻孔资料对含水层组及其粘性土夹层或透镜体的分析和野外地下水采样、动态监测及地面沉降测量资料等分析,采用分层、分区计算方法对漏斗中心至外围的各分区弹性释水量、粘性土压密释水量、侧向流入和越流补给量进行均衡计算。结果表明:在强烈开采条件下,深层水开采量中弹性释水量和压密释水量所占比率不是常量,随着深层水位下降的幅度和深度增大、超采释压持续时间愈长、开采层中粘性土夹层或透镜体愈多,则区内弹性释水和压密释水的强度及累积地面沉降量愈大;愈靠近漏斗中心,弹性和压密释水强度愈大,相同降深下平均单位降深的地面沉降量愈小。在深层水开采量中粘性土压密释水占35.5%~61.9%,与地面沉降密切相关。     

综合物探技术在煤系上覆薄硬岩矿区的应用初探——以山西潞安集团司马矿区为例

山西潞安集团司马矿区东高西低,西部延至长治盆地,矿区水系发育。主采3#煤层埋藏较浅,层位稳定。新生界不整合于下伏基岩之上,厚度变化较大,属典型的薄硬岩煤矿区。为准确提供矿区煤炭开采上限的地质依据,采用了综合勘探技术,其中利用三维地震勘探技术追踪新生界底面反射波及第四系"天窗"位置,利用直流电法与瞬变电磁法提供电性分层,并对其解释结果进行相互印证。勘探实例表明,三维地震勘探、电法勘探与钻探的综合运用不仅可解决薄硬岩煤矿区新生界厚度、含水层富水性及第四系缺失等问题,还可提供新生界含水层与基岩含水层之间的水力联系、制约因素等资料。     

山西新阳矿太原组下组煤层底板岩溶水突水安全性评价

新阳矿9-10-11号煤层底板低于奥灰岩溶水位,承受水压较高,存在带压开采问题。依据矿区水文地质勘查资料,通过对9-10-11号煤层底板隔水岩层隔水性能、奥灰岩溶水富水程度、底板断裂构造发育情况、煤层底板承压大小以及采煤扰动底板破坏深度等突水因素分析,表明井田西南部断裂构造发育地段为突水危险区。分别采用突水系数法和突水危险度法对全井田9-10-11号煤层底板突水危险性进行评价,评价结果显示井田南部及断裂构造发育地段为突水危险区。理论分析与定量评价结果基本一致。针对9-10-11号煤层底板突水安全评价结果,提出了相应的防突水对策。     

56m深TRD工法搅拌墙在深厚承压含水层中的成墙试验研究

上海国际金融中心项目基坑面积约为48 860m2,开挖深度为26.5²7.9m,周边环境复杂。为控制抽降承压水对周边环境的影响,经方案比选,基坑周边设置厚700mm、深56m等厚度水泥土搅拌墙(TRD)作为承压水悬挂隔水帷幕。在上海地区施工如此深TRD墙体尚属首例,为此现场开展了试成墙试验,试成墙监测表明,墙身在深厚承压含水层中水泥土强度达到0.8427.9m,周边环境复杂。为控制抽降承压水对周边环境的影响,经方案比选,基坑周边设置厚700mm、深56m等厚度水泥土搅拌墙(TRD)作为承压水悬挂隔水帷幕。在上海地区施工如此深TRD墙体尚属首例,为此现场开展了试成墙试验,试成墙监测表明,墙身在深厚承压含水层中水泥土强度达到0.84¹.38 MPa。室内渗透性试验表明,渗透系数由10-3cm/s提高到10-7cm/s,满足隔水帷幕设计要求;墙体施工期间,地表最大沉降约8mm,主要影响范围约5m;土体侧向位移主要产生在距离墙体5m的范围内,TRD墙体施工对周边环境影响很小。试验墙体的顺利实施为后续正式墙体的施工提供了依据,也为类似工程提供了重要参考。     

矿山开发中的含水层破坏评估方法初探——以白山市恒基煤矿为例

本文以白山市恒基煤矿为例,在阐述矿区地质环境基础上,探讨了矿山开采中的含水层结构破坏评估方法,含水层水位、水量变化、和采矿活动对含水层水质造成的污染。最后提出了防治措施,为矿山地质环境恢复治理提供科学依据,也为其他矿山开采提供参考。     

陕北榆横矿区小纪汗井田水文地质特征及矿床充水条件分析

从用顶板充水、底板突水两个角度对矿坑充水强度进行了分析,为小纪汗井田及周边开采提供了参考。