盆地地下水系统的数值模拟研究——以山西灵石静升盆地为例

以山西灵石静升盆地为例,选择三维流数学模型,合理刻画了模型边界,利用CFC定年技术对模型进行识别校正,提高了模型仿真程度,构建了适合盆地地下水系统的数值模型.最后对不同水平年盆地水动力特征进行了输出分析.     

CFCs法在鄂尔多斯白垩系地下水盆地浅层地下水年龄研究中的应用

鄂尔多斯白垩系地下水盆地是中国重要的能源基地,利用CFCS方法对鄂尔多斯白垩系地下水盆地年轻地下水年龄进行测定,取得了比较准确的结果。结果表明,盆地地下水系统可以分为局部地下水系统、中间地下水系统和区域地下水系统。地下水盆地内局部地下水系统地下水循环更替快,年龄为20 a左右,地下水年龄随深度增加而增加;中间地下水系统和区域地下水系统地下水循环更替相对较慢,年龄大于70 a,可以寻求14C等其他同位素方法确定准确年龄。     

单井混合抽水试验分层计算含水层渗透系数测试研究

本文根据施普德混合抽水试验的基本计算原则,在既无观测孔、也无法分层止水的情况下,对新疆焉耆盆地南部山前5号水井,开展了单井混合抽水试验分层计算含水层渗透系数的测试研究。首先采用图解法原理,计算出混合水井下部承压含水层的静止水位为47.232 m,并绘制了该含水层的Q-S曲线。在此基础上,分层计算得到了承压含水层的渗透系数Kb=6.968m d-1,潜水含水层的渗透系数Ka=4.979m d-1。通过与实际抽水试验计算结果的对比分析,验证了该方法的有效性和准确性,可为单井混合抽水试验分层计算含水层渗透系数提供参考。     

敦煌盆地早更新世沉积物粒度分析、36Cl定年及其构造隆升意义

青藏高原早更新世末期的快速隆升对全球气候变化、中国西部盆山地貌形成与荒漠化的发展有着重大影响。高原东北缘敦煌盆地早更新世河湖相沉积物粒度分析和红色泥岩³⁶Cl 断代法定年结果表明：早更新世晚期（1.164 ~1.087 Ma B.P.）,青藏高原北部党金山快速隆升,砾石层和粗砂层发育,沉积速率加快;1.087~ 0.809 Ma B. P.,山脉隆升速度变缓,沉积速率降低;0.809 Ma B. P.以来,受山脉隆升的影响,敦煌盆地气候逐渐变得干旱,沙漠形成。区域分析表明,沉积速率所反映的中国中西部山脉隆升最新的年代存在差别,不同时期高原的隆升产生不同的气候效应。因此,构造隆升和气候变化的细节尚待进一步研究。     

下扬子天目山盆地火山岩锆石LA-ICP-MS定年及地质意义

天目山盆地是下扬子江南隆起带保存较完整的中生代火山盆地,中生代火山岩系岩性自下而上主要为流纹岩-英安岩-安山岩。对盆地内黄尖组下段流纹岩和英安岩分别进行了锆石 LA-ICP MS定年,分别获得了133.6±1.5 Ma(MSWD=0.73)和135.0±2.1 Ma(MSWD=0.78)的锆石U-Pb年龄,指示天目山盆地黄尖组火山岩时代为早白垩世。天目山盆地火山活动起始时间和长江中下游地区晚中生代火山活动基本一致,说明江南隆起带和长江中下游地区在早白垩世均处于强烈拉张环境。     

鄂尔多斯盆地东北部差异隆升过程裂变径迹分析

基于磷灰石裂变径迹(AFT)的分析方法,探讨鄂尔多斯盆地东北缘差异隆升过程及其隆升强度,为鄂尔多斯盆地东北缘(晋西挠摺带府谷—吴堡区段)构造演化历史及其与多种能源矿产耦合关系提供新的认识。不同构造单元及其不同层系样品的AFT分析表明:研究区北段府谷—兴县地区构造抬升相对较早,且经历了白垩纪晚期(86~56Ma)和古近纪(44~37 Ma)两次隆升过程,平均隆升速率分别为24.5 m/Ma和41.8 m/Ma;研究区中段紫金山地区抬升相对较晚,主控构造事件发生在晚白垩世末期—古近纪早期(68~56 Ma)和古近纪中晚期(35 Ma),平均隆升速率分别为48.8 m/Ma和49.2 m/Ma;研究区南段临县—吴堡地区抬升最晚(35~21 Ma),平均隆升速率为73.9 m/Ma。因此,鄂尔多斯盆地东北缘晚白垩世以来的差异隆升过程具有北段抬升早、中段抬升相对较晚和南段抬升更晚的特点,南北区段统一的强烈构造抬升活动主要发生在古近纪以来的晚近时期,且构造隆升强度由南向北逐渐减弱。结合已有的成矿(藏)年代学资料分析表明,鄂尔多斯盆地东北缘关键构造事件及其隆升强度与多种矿产耦合成矿(藏)事件关系密切,构造事件与成藏(矿)事件呈现出显著的协同耦合特点。     

黑龙江省鸡西盆地基性岩40Ar/39Ar同位素定年及其地质意义

黑龙江省鸡西盆地张新地区辉绿玢岩呈岩床状侵入晚中生代城子河组地层中;辉绿玢岩较纯净,无捕虏体,呈斑状结构,斑晶由普通辉石、斜长石组成,斑晶体积分数为35%～40%;基质由微晶斜长石、辉石、磁铁矿及玻璃质组成。辉绿玢岩⁴⁰Ar/³⁹Ar同位素定年结果表明,其等时线年龄与坪年龄在误差范围内较为接近,该辉绿玢岩的形成时代介于96～101 Ma,即形成于早白垩世晚期-晚白垩世早期的Albian-Cenomanian阶。结合区域上该期岩浆事件的特点,暗示早白垩世晚期-晚白垩世时期鸡西盆地处于古太平洋板块(Izanagi-Kula 板块 )向东亚大陆边缘斜向俯冲背景之下的弧后伸展环境。     

安徽繁昌盆地牛山岩体锆石U-Pb定年及其意义

牛山岩体位于长江中下游地区繁昌火山盆地内,为一钻孔揭露的隐伏花岗斑岩岩体。近年来,勘探发现了与该岩体密切相关的铜、金、铅、锌多金属矿化。应用LA-ICP-MS锆石微区分析技术对该岩体进行了精确定年。岩体含有大量岩浆成因锆石,以高的Th/U比和振荡结晶环带为特征,17个测点的206Pb/238U年龄加权平均值为(121.8±1.9)Ma,属于长江中下游地区白垩纪岩浆晚期活动的产物。该岩体的精确定年和成矿发现,为本地区火山盆地内寻找铜、金矿开辟了新的前景。     

运城盆地地下水同位素年龄特征及其演化过程和可持续利用

通过分析运城盆地地下水的碳同位素组成,结合水化学特征,揭示了盆地深层承压地下水的补给期为22～3 ka BP （现代碳百分比(a¹⁴C) 6~38 pmC）。浅层地下水（71~89 pmC）由现代水或现代水和老水混合组成。深层地下水氢氧同位素组成特征（δ¹⁸O~-10‰; δ²H~-70‰）表明地下老水在气候较冷的环境下受到补给,而浅层地下水的氢氧同位素组成（δ¹⁸O~-8‰; δ²H～-51‰）特征与现代西安降水组成相似。浅层地下水NO^-₃平均含量（31mg/L）比深层地下水（1.8 mg/L）高,硝酸盐的δ¹⁵N－δ¹⁸O_NO3 组成 (0‰～5‰）揭示了硝酸盐的主要来源为综合肥料。此外,浅层地下水的TDS由于蒸散发、矿物溶解,可达8.5 g/L（平均2.0 g/L）,深层地下老水TDS可达1.8 g/L（平均1.1g/L）水质相对较好。研究区目前主要开采深层地下水,受断裂带影响,浅层地下水已经侵入中深层地下水并与之发生混合,严重影响了中深层地下水的水质。如果发生大规模的浅层地下水与中深层地下水混合,会造成中深层地下老水的NO^-₃、TDS等含量越来越高。     

吐哈盆地南缘克孜尔卡拉萨依和大南湖花岗质岩基锆石SHRIMP定年及其地质意义

报道了采自吐哈盆地南缘大南湖和克孜尔卡拉萨依花岗质岩基的锆石SHRIMP定年结果。采自大南湖岩基的二长花岗岩中的10粒锆石~(206)Pb/~(238)U年龄为363～400Ma,加权平均为(383±9)Ma,采自克孜尔卡拉萨依岩体花岗闪长岩的11粒锆石加权平均~(206)Pb/~(238)U年龄为(357.3±6.2)Ma,分别为中泥盆世早期和石炭纪初。说明该区大规模的陆缘型中酸性岩浆侵入活动在中泥盆世早期就开始发育,并可能发育古老的陆壳基底。     

Large-scale hydraulic conductivity distribution in an unconfined carbonate aquifer using the ocean tidal propagation

The hydraulic conductivity of an unconfined carbonate aquifer at the uplifted atoll of Minami-Daito, Japan, was evaluated by a combination of cross-spectral analysis, analytical solution, and density-dependent groundwater modeling based on observed groundwater levels in 15 wells and at sea level. The island area was divided into 10 subregions based on island morphology and on inland propagation of ocean tides. The hydraulic conductivity was obtained for each subregion using analytical solutions based on phase lags of M₂ constituents of ocean tides at each well by assuming two aquifer thicknesses (300 and 1,800 m) and two effective porosities (0.1 and 0.3). The density-dependent groundwater model evaluated the hydraulic conductivity of the subregions by reproducing observed groundwater levels. The hydraulic conductivity in the subregions was estimated as 3.46?×?10^?3 to 6.35?×?10^?2 m/s for aquifer thickness of 300 m and effective porosity of 0.1, and as 1.73?×?10^?3 to 3.17?×?10^?2 m/s for aquifer thickness of 1,800 m and the effective porosity of 0.3. It was higher in southern and northern areas, and higher in interior lowland than in the western and eastern areas. Fissures and dolomite distributions on the island control differences of the omnidirectional ocean tidal propagation and cause these differences in hydraulic conductivity. The method used for this study may also be applicable to other small islands that have few or no data for hydraulic conductivity.

     

Modeling water-table fluctuations in a sloping aquifer with random hydraulic conductivity

To prevent environmental problems like water logging and increase in soil salinity which are responsible for the degradation of the top productive soils, an optimum ditch drainage design is required. For this purpose a knowledge of the spatio-temporal variation of the water table is essential. In this study the spatio-temporal variation of the water table in a sloping ditch drainage system has been modeled from a stochastic point of view, incorporating randomness in hydraulic conductivity to get the expression for the mean and the standard deviation of the water-table height. The hydraulic conductivity has been considered to be a realization of a log-normal distribution. Application of these expressions in the prediction of mean water-table variation with the associated error bounds has been demonstrated with the help of a ditch drainage problem of a sloping aquifer. The sensitivity analysis has also been carried out to see the effect of variability in the hydraulic conductivity on the water-table fluctuations. The error bounds quantified on the water-table height will thus help in the decision-making process for proper drainage design.     

咸淡水驱替过程中含水介质渗透性变化的试验研究

在野外水文地质调查的基础上,采集青岛市大沽河下游咸水入侵区砂样,首先对含水介质的粒度和矿物组成进行了分析,然后通过室内砂槽模拟试验,对咸淡水驱替过程中含水介质的渗透性变化规律进行了研究。结果表明,当咸淡水相互驱替时,含水介质的渗透性会发生显著变化,这种变化主要是由于驱替液盐浓度的变化,使得伊利石、高岭石和绿泥石等非膨胀性粘土矿物经释放、迁移、絮凝和沉积,从而重新分布所引起的。含水介质渗透性的变化具有明显的非均质性,即位于距入水处不同水平距离以及不同高度处的含水介质,其渗透性变化规律有所不同。     

含水层渗透系数空间变异性对地下水数值模拟的影响

地下水数值模拟是目前定量研究地下水水量和水质的重要手段。使用基于随机理论的MonteCarlo方法来进行地下水数值模拟。这种方法能较好地考虑水文地质参数的空间变异性。主要将MonteCarlo方法和确定性模型模拟方法的模拟结果在渗透系数场、水头场、速度场和浓度场等方面进行了比较。结果表明:在模拟三维非均质含水层中的溶质运移问题时,充分考虑了含水层渗透系数空间变异性的MonteCarlo法比确定性方法更为有效,模拟精度提高了很多,且对模拟误差及误差来源有合理的数学解释。     

磁共振探测估算含水层渗透系数的原位试验研究

通过一个完整的磁共振探测（MRS）、钻探、抽水试验过程,开展了MRS估算含水层渗透系数的原位试验研究。在分析MRS估算渗透系数准确性的基础上,系统剖析了引起准确性差异的主要因素,并对进一步提高估算效果提出了建议。结果显示：（1）采用现阶段常用的参数设置,与抽水试验计算值相比,MRS估算渗透系数的差值为抽水试验计算值的17.59%。（2）MRS推断的含水层顶板、底板埋深与钻孔揭示的信息相比,差值分别为4.11m、1.03m,表明其对估算渗透系数准确性的影响较小。（3）CP值是影响估算准确性的重要因素,其值为10-9数量级符合大多数地层的特点。另外,通过在已知渗透系数的钻孔附近进行MRS,从而获取CP参考值,应用该值估算的渗透系数准确性高于常用的参数设置。（4）指数a、b设置方面,应用Seevers公式（a=1、b=2）的估算效果优于Kenyon公式（a=4、b=2）。本成果有助于提高MRS估算渗透系数方法在野外条件下的适用性。     

Correlation of hydraulic conductivity with stratigraphy in a fractured-dolomite aquifer, northeastern Wisconsin, USA

This study integrated surface and subsurface stratigraphic data with geophysical logs and hydrogeologic data in order to characterize the hydraulic properties of the Silurian dolomite in northeastern Wisconsin. Silurian stratigraphy consists of predictable alternations of characteristic facies associations. A vertical profile of hydraulic conductivity, obtained from short-interval packer tests in a core hole that penetrates a majority of the Silurian section, indicates that hydraulic conductivity ranges over five orders of magnitude (10^–1 to 10^–6 cm/s). Matrix conductivity is generally low and varies with texture; the finer-grained restricted-marine and transitional facies being less permeable than the coarser-grained open-marine facies. High-conductivity values are generally associated with bedding-plane fractures, and fracture frequency is greater in the restricted-marine facies. High-flow features in 16 wells were identified using fluid-temperature/resistivity and heat-pulse flowmeter logs. Natural-gamma logs were used to identify the stratigraphic position of flow features in each well and to correlate flow features to specific stratigraphic horizons. By combining stratigraphic, geophysical, and hydrogeologic data, 14 high-permeability zones within the Silurian aquifer have been identified and correlated in wells more than 16 km apart. These zones parallel bedding, appear most pronounced at contacts between contrasting lithologies, and are most abundant in restricted-marine lithologies. Electronic Publication     

Estimating hydraulic conductivity using grain-size analyses, aquifer tests, and numerical modeling in a riverside alluvial system in South Korea

Hydraulic conductivity (K) for an alluvial system in a riverbank filtration area in Changwon City, South Korea, has been studied using grain-size distribution, pumping and slug tests, and numerical modeling. The alluvial system is composed of layers: upper fine sand, medium sand, lower fine sand, and a highly conductive sand/gravel layer at the base. The geometric mean of K for the sand/gravel layer (9.89?×?10^?4 m s^?1), as determined by grain-size analyses, was 3.33 times greater than the geometric mean obtained from pumping tests (2.97?×?10^?4 m s^?1). The geometric mean of K estimates obtained from slug tests (3.08?×?10^?6 m s^?1) was one to two orders of magnitude lower than that from pumping tests and grain-size analyses. K estimates derived from a numerical model were compared to those derived from the grain-size methods, slug tests and pumping tests in order to determine the degree of deviation from the numerical model. It is considered that the K estimates determined by the slug tests resemble the uppermost part of the alluvial deposit, whereas the K estimates obtained by grain-size analyses and pumping tests are similar to those from the numerical model for the sand/gravel layer of the riverside alluvial system.     

The hydraulic conductivity field and groundwater flow in the unconfined aquifer system of the Keta Strip,Ghana

This study investigates the hydraulic conductivity field and the groundwater flow pattern as predicted by a calibrated steady state groundwater flow model for the Keta Strip, southeastern Ghana. The hydraulic conductivity field is an important parameter in evaluating aquifer properties in space, and in general basin-wide groundwater resources evaluation and management. This study finds that the general hydraulic conductivity of the unconsolidated unconfined aquifer system of the Keta Strip ranges between 2 m/d and 20 m/d, with an average of 15 m/d. The spatial variation in horizontal hydraulic conductivity appears to take the trend in the variations in the nature of the material in space. Calibrated groundwater recharge suggests that 6.9–34% of annual precipitation recharges the shallow aquifer system. This amount of recharge is significant and suggests high fortunes in terms of groundwater resources development for agriculture and industrial activities in the area. A spatial distribution of groundwater recharge from precipitation is presented in this study. The spatial pattern appears to take the form of the distribution in horizontal hydraulic conductivity, and suggests that the vertical hydraulic conductivity takes the same pattern of spatial variation as the horizontal hydraulic conductivity. This is consistent with observations in other areas. The resulting groundwater flow is dominated by local flow systems as the unconfined system is quite shallow. A general northeast – southwest flow pattern has been observed in the study area.     

Variations in hydraulic conductivity with scale of measurement during aquifer tests in heterogeneous, porous carbonate rocks

 Previous studies have shown that hydraulic conductivity of an aquifer seems to increase as the portion of the aquifer tested increases. To date, such studies have all relied on different methods to determine hydraulic conductivity at each scale of interest, which raises the possibility that the observed increase in hydraulic conductivity is due to the measurement method, not to the scale. This study analyzes hydraulic conductivity with respect to scale during individual aquifer tests in porous, heterogeneous carbonate rocks in southeastern Wisconsin, USA. Results from this study indicate that hydraulic conductivity generally increases during an individual test as the volume of aquifer impacted increases, and the rate of this increase is the same as the rate of increase determined by using different measurement methods. Thus, scale dependence of hydraulic conductivity during single tests does not depend on the method of measurement. This conclusion is supported by 22 of 26 aquifer tests conducted in porous-flow-dominated carbonate units within the aquifer. Instead, scale dependency is probably caused by heterogeneities within the aquifer, a conclusion supported by digital simulation. All of the observed types of hydraulic-conductivity variations with scale during individual aquifer tests can be explained by a conceptual model of a simple heterogeneous aquifer composed of high-conductivity zones within a low-conductivity matrix. Received, January 1997 Revised, August 1997, November 1997 Accepted, November 1997