大连大魏家水源地海水入侵过程中水文地球化学作用分析及定量模拟

大连大魏家水源地位于中国北方典型滨海岩溶地区。近年来,地下淡水不合理开采造成的地下水位降落漏斗引发了严重的海水入侵。本文以大魏家水源地海水入侵为研究对象,通过大量的水文地质调查和水化学样品测试分析,探讨了海水入侵形成的水动力条件,通过绘制Piper图、以Cl^-为天然示踪剂计算了海水和地下淡水的混合比例、实际浓度和理论浓度的对比、水化学组成特征分析等方法,探讨了该区岩溶含水层在海水入侵过程中所发生的主要水文地球化学作用,并对其进行了定量模拟,进而阐明了岩溶含水层中的海水入侵机理。研究结果表明,研究区微成水和成水并不是地下水淡水和海水简单混合而成。利用反向水文地球化学模拟揭示了控制滨海岩溶含水层中水化学演化的主要水文地球化学反应有方解石、蒙脱石和石膏的溶解作用,伊利石的沉淀作用以及Ca—Na离子交换作用,伴随着CO2的释放,石膏的溶解并没引起白云石化作用。     

环境同位素特征对滨海岩溶地区海水入侵过程的指示意义

大连大魏家水源地位于中国北方典型滨海岩溶地区。近30年来,地下淡水的不合理开采造成的地下水位降落漏斗引发了严重的海水入侵。以大魏家水源地为研究对象,通过大量的水文地质调查和水化学及同位素采样测试分析,探讨海水入侵形成的水动力条件,通过分析滨海岩溶含水层中地下水主要水化学和多种同位素(δ2H-δ18O,δ34S,δ13C)组成特征,识别了海水入侵过程中发生的主要水文地球化学作用,并对其进行了定量模拟,从而阐明了岩溶含水层中的海水入侵机理。研究结果表明:大连大魏家海水入侵主要通道为大魏家地区存在的导水断裂、岩溶裂隙以及第四系松散地层。对δ2H-δ18O同位素的组成分析表明,研究区地下水主要来自大气降水补给,结合Cl-浓度分布,认为除海水入侵淡水含水层后增加了地下水中的盐分外,浅层地下水的蒸发也对地下水中盐分的累积起到了重要作用。根据不同水体中δ34SSO4,δ13CHCO3等同位素特征,结合水化学成分(如SO2-4,Cl-)分析认为,研究区微咸水和咸水并不是地下水淡水和海水简单混合而成。利用反向水文地球化学模拟揭示了控制滨海岩溶含水层中水化学演化的主要水文地球化学反应有方解石、蒙脱石和石膏的溶解作用,伊利石的沉淀作用以及Ca-Na离子交换作用,伴随着CO2的释放。     

中国沿海地区海水入侵机理及防治措施研究

早在1964年,中国大连市就发生了过量开采地下水引发的海水入侵现象。目前有十几座沿海城市发生了不同程度的海水入侵。海水入侵不仅恶化了滨海地区的生态环境,而且阻碍了沿海地区的经济发展。中国沿海地区海水入侵面积经历了20世纪80年代的快速增长后,近些年来呈缓慢增长态势,甚至有些地区海水入侵面积不再增加。为了深入了解海水入侵的机制,从20世纪70年代后期开始,有关学者相继开展了一系列研究,如通过水-岩间阳离子交换实验来了解海水入侵的地球化学作用,建立三维有限元变密度的过渡带数值模型模拟与预报海水入侵过程及发展趋势。为了预防和治理海水入侵,除了控制地下水开采量及开采布局外,沿海地区还园地制宜采取了相应的防治措施,包括人工补给、地下防渗帷幕、节水灌溉以及远距离调水等工程措施。     

海水入侵研究现状与展望

文中讨论了海水入侵的研究进展：海岸带海水入侵界面或过渡带的模拟，海岛地下淡水透镜体的动力学研究，井孔在界面或过渡带上部抽水所引起的升锥问题以及岸边地下水位动态研究及海平面变化对海水入侵的影响。海水入侵研究经历了从理想假定到合理概化，从室内实验模型、理想模型到数值模型这一过程。数值方法已成为模拟和求解海水入侵问题的最有力工具，海水入侵的过渡带模型成为主要的研究方向。分析了研究中存在的困难，展望了海水入侵研究的未来趋势：边界条件确定，水文地质参数选取，含水层结构估计，海水入侵输运与化学作用的耦合以及海水入侵预报。     

高密度电阻率探针原位监测海水入侵过程试验研究

为研究砂质含水层中海水入侵问题,确定海水入侵过程中含水介质内海水入侵楔形体的发展演化过程,本文采用砂槽实验模拟海水入侵,利用高密度电阻率探针观测系统远程实时监测海水入侵过程中土体的垂向电阻率变化,据电阻率变化曲线分析海水入侵发生发展过程。研究结果表明:高密度电阻率探针对海水入侵界面变化反应灵敏,能够监测到咸淡水突变界面形成、过渡带演化及咸淡水界面变迁等一系列现象,较好地实时监测海水入侵发生发展过程及海水入侵状态并进行预警。     

基于FEFLOW的海水入侵数值模拟

在系统分析珠江口东岸地区地质及水文地质条件的基础上,建立了研究区海水入侵三维溶质模型,利用基于有限元原理的FEFLOW软件对建立的模型进行求解,用地下水位和浓度的动态观测资料对模型进行了识别和校验。运用识别后的模型预测了在不同开采条件下,研究区海水入侵的趋势。研究结果表明减少地下水开采量,海水入侵面积可逐步减少,入侵可得到有效控制;在减少开采量的情况下,相较于基岩裂隙含水层,第四系含水层海水入侵面积收缩较慢。     

含水层中海水入侵的数学模型

提出了考虑过渡带的适合潜水含水层和承压含水层的海水入侵数学模型,导出了相应的方程和定解条件.经实例验证,模型是合理和可靠的.     

控制海水入侵的地下水多目标模拟优化管理模型

为实现滨海含水层地下水开采-回灌方案优化、控制海水入侵面积和降低海水入侵损失等多重管理目标,建立了海水入侵条件下地下水多目标模拟优化管理模型SWT-NPTSGA。模拟模型采用基于变密度流的数值模拟程序SEAWAT来模拟海水入侵过程。优化模型采用小生境Pareto禁忌遗传混合算法NPTSGA来求解,该算法在保证多目标权衡解的收敛性和计算效率的前提下,能维护整个进化群体的全局多样性。将SWT-NPTSGA程序应用于一个理想滨海含水层地下水开采方案和人工回灌控制海水入侵的优化设计中,结果表明该管理模型能够同时处理最大化总抽水流量、最小化人工回灌总量和最小化海水入侵范围等3个目标函数之间的权衡关系。通过采用人工回灌海水入侵区的减灾策略,既能增加滨海地区的供水量,又可减少海水入侵的范围,由此进一步验证了模型的有效性和可靠性。     

西班牙巴塞罗那Llobregat三角洲含水层中的咸水

位于西班牙巴塞罗那西南的Ilobregat三角洲含水层中含有重要的地下水资源,其详细描述见其它文章(Custodio,1981年)。由于承压地层中强烈开采地下水,出现了海水污染问题,近期入侵的海水不是咸水唯一的可能来源。要建立有效管理准则就必须充分了解海水入侵,要有效地防止和制止海水入侵,很大程度上取决于该准则。上覆不透水层为一非常新的地层,它含有原生源的咸水。几年前就很好地确定出了氯离子剖面图推出的垂直渗透率资料,该剖面易于符合数值模拟模型结果。为更详细了解深层含水层和弱含水层咸水特征,对深层含水层和弱含水层孔隙水重     

基于GIS的海水入侵可视化数值模拟系统

将变密度三维海水入侵模型(SUTRA模型)与地理信息系统(GIS)进行有机集成开发出可视化海水入侵数值模拟系统.该系统充分利用了GIS的空间分析和数据管理功能及其可视化功能,实现了SUTRA模型的可视化构建1、含水层空间和模拟时间的可视化离散、各种参数的可视化赋值、模型的可视化模拟及其模拟结果的可视化表达等功能,为海水入侵研究提供了有效的模拟和分析手段.     

Seawater intrusion and associated processes in a small coastal complex aquifer (Castell de Ferro,Spain)

An attempt to identify the hydrogeochemical processes that accompany current or past intrusion of seawater in the Castell de Ferro coastal aquifer (SE Spain) was made using ionic deltas, ionic ratios and saturation indices. Alluvial and coastal Quaternary sediments and Triassic marbles constitute this complex aquifer. Cation exchange, dolomitization, and calcite dissolution were identified as operating in the aquifer. The cation exchange processes in the detrital materials are different from those in the marbles. In addition, it was determined that while the processes of marine intrusion associated with periods of low groundwater levels still persist in the detrital aquifer, in the carbonate sector the flushing process characteristic of wetter conditions has begun.     

渗透回流模式白云石化作用对碳酸盐岩储层的影响

白云石化作用对含油气的碳酸盐岩储层有明显的改造作用,而前人的研究大部分是针对白云石的成因模式,未对其形成过程给出具体描述,更未对白云石化作用对储层的改造程度做出定量分析。针对白云石化中的“渗透回流模式”,建立了白云石化作用的一维模型,利用多组分、多相流溶质运移模拟软件TOUGHREACT,研究渗透回流模式对储层物性影响的深度范围。结果表明,当高Mg/Ca比的浓缩海水进入含水层后,导致方解石溶解,白云石沉淀,储层孔隙度与渗透率均比初始值大,模拟5万年时,依靠重力作用,渗透回流白云石化作用能够影响到50 m深的储层范围,产生一定规模的白云岩储层。     

Saltwater intrusion and nitrate pollution in the coastal aquifer of Dar es Salaam, Tanzania

Dar es Salaam Quaternary coastal aquifer is a major source of water supply in Dar es Salaam City used for domestic, agricultural, and industrial uses. However, groundwater overdraft and contamination are the major problems affecting the aquifer system. This study aims to define the principal hydrogeochemical processes controlling groundwater quality in the coastal strip of Dar es Salaam and to investigate whether the threats of seawater intrusion and pollution are influencing groundwater quality. Major cations and anions analysed in 134 groundwater samples reveal that groundwater is mainly affected by four factors: dissolution of calcite and dolomite, weathering of silicate minerals, seawater intrusion due to aquifer overexploitation, and nitrate pollution mainly caused by the use of pit latrines and septic tanks. High enrichment of Na⁺ and Cl^? near the coast gives an indication of seawater intrusion into the aquifer as also supported from the Na–Cl signature on the Piper diagram. The boreholes close to the coast have much higher Na/Cl molar ratios than the boreholes located further inland. The dissolution of calcite and dolomite in recharge areas results in Ca–HCO₃ and Ca–Mg–HCO₃ groundwater types. Further along flow paths, Ca²⁺ and Na⁺ ion exchange causes groundwater evolution to Na–HCO₃ type. From the PHREEQC simulation model, it appears that groundwater is undersaturated to slightly oversaturated with respect to the calcite and dolomite minerals. The results of this study provide important information required for the protection of the aquifer system.     

Identification of hydrogeochemical process linked to marine intrusion induced by pumping of a semiconfined mediterranean coastal aquifer

This article analyses the hydrogeochemical processes, linked to the freshwater–seawater mixing zone, which can be caused by continuous pumping from a detrital coastal aquifer. It was necessary to construct an experimental plot, drilling three boreholes along a line perpendicular to the coast. A complete physico-chemical analysis was done of all water samples taken. The percentage of seawater, calculated from Chloride and ¹⁸Oxygen concentrations, varied between 55 and 90 %. The ionic deltas (Δ) calculated, and the saturation indices (SI) of mineral phases susceptible to precipitation or dissolution, allowed a series of hydrogeochemical processes to be identified that occur as a consequence of the advance of marine intrusion into the coastal band, and of aquifer flushing. Based on the major elements, the fraction of exchange (βI) was calculated for samples ranging from seawater to freshwater, and this revealed that differences in βI could explain the hydrochemistry of the mixing zone. The main processes recognised include precipitation of dolomite, dissolution of gypsum, fixation of sulphur salts and cation exchange. Most of the ion exchange took place between Na and Calcium + Magnesium ions. The process of fixation or liberation of these ions is probably determined by the advance or recession of the saline wedge, and/or by recharge during rainy periods. The behaviour of Magnesium is more sensitive to small variations in salinity, whilst Calcium behaves more homogeneously. The high percentage of seawater in the samples studied favours the speed and magnitude of processes such as ion exchange, and the intervention of magnesium is also a key.     

The integrated impacts of natural processes and human activities on groundwater salinization in the coastal aquifers of Beihai,southern China

Salinization in coastal aquifers is usually related to both seawater intrusion and water–rock interaction. The results of chemical and isotopic methods were combined to identify the origin and processes of groundwater salinization in Daguansha area of Beihai, southern China. The concentrations of the major ions that dominate in seawater (Cl^?, Na⁺, Ca²⁺, Mg²⁺ and SO ₄ ^2– ), as well as the isotopic content and ratios (²H, ¹⁸O, ⁸⁷Sr/⁸⁶Sr and ¹³C), suggest that the salinization occurring in the aquifer of the coastal plain is related to seawater and that the prevailing hydrochemical processes are evaporation, mixing, dissolution and ion exchange. For the unconfined aquifer, groundwater salinization has occurred in an area that is significantly influenced by land-based sea farming. The integrated impacts of seawater intrusion from the Beibuwan Gulf and infiltration of seawater from the culture ponds are identified in the shallowest confined aquifer (I) in the middle of the area (site BBW2). Leakage from this polluted confined aquifer causes the salinization of groundwater in the underlying confined aquifer (II). At the coastal monitoring site (BBW3), confined aquifer I and lower confined aquifer II are heavily contaminated by seawater intrusion. The weak connectivity between the upper aquifers, and the seaward movement of freshwater, prevents saltwater from encroaching the deepest confined aquifer (III). A conceptual model is presented. Above all, understanding of the origin and processes of groundwater salinization will provide essential information for the planning and sustainable management of groundwater resources in this region.     

Hydrogeochemical Analysis of Processes Through Modeling of Seawater Intrusion Impacts in Biscayne Aquifer Water Quality,USA

Hydrogeochemical processes that accompany seawater intrusion in coastal aquifers can alter the resulting water quality and are important ingredients in coastal aquifer management. The presence of dissolution–precipitation reactions and ion exchange in the mixing zone of the Biscayne aquifer (FL, USA) are suggested based on changes in major ion concentrations and mineral saturation indices (SI). Major ion concentrations from 11 groundwater samples are compared with theoretical mixing between freshwater and seawater. PHREEQC code was used to calculate saturation indices of the samples with respect to common phases in the Biscayne aquifer. High Ca²⁺ and HCO₃ ^? content of the samples is typical of waters in contact with carbonate aquifers. Water quality of the samples is mainly attributed to mixing and precipitation–dissolution reactions with calcite and dolomite. The samples were saturated with calcite (SI ~ 0) and undersaturated for dolomite (SI < 0), while a few samples showed dolomite saturation. Because gypsum and halite SI could be predicted by theoretical mixing, reactions with those minerals, if present, are thought to be insignificant. In the active intrusion areas, cation exchange also appears to modify water quality leading to excess Ca²⁺, but depleted Na⁺, Mg²⁺ and K⁺ concentrations. On the other hand, samples from previous intrusion areas plotted very close to the theoretical mixing line and approached equilibrium with the seawater.     

Hydrogeochemistry and geochemical simulations to assess water–rock interactions in complex carbonate aquifers: The case of Aguadulce (SE Spain)

The hydrogeological unit of Aguadulce (Campo de Dalías aquifers, SE Spain) has a complex geometry. This fact, together with a continuous rise in water demand due to intensive agriculture and tourism create problems for groundwater quantity and quality. In this paper classic geochemical tools managed by means of GIS software and geochemical simulations are combined to delineate, identify and locate the possible physicochemical processes acting in the Aguadulce groundwater. Two main aquifers can be distinguished: the carbonate or lower aquifer of Triassic age, and the calcodetritic or upper aquifer of Plio-Quaternary age. Groundwaters from the latter are more saline and, assuming all chlorinity originates from seawater intrusion, the seawater contribution to their composition would be up to 7%. Nevertheless the carbonate aquifer appears not to be homogeneous: it is compartmentalised into 4 zones where different processes explain the different groundwaters compositions. Zone 4 samples (E margin of the carbonate aquifer) resemble those of the Plio-Quaternary aquifer, where calcite precipitation, dolomite and gypsum dissolution and some cation exchange (water–rock interaction) together with seawater–freshwater mixing occur. In contrast, water–rock interaction predominates in zones 1 and 3 of the carbonate aquifer. Moreover, zone 2 samples, located between zones 1 and 3, are explained by water–rock interaction in addition to mixing with Plio-Quaternary aquifer waters. The combination of geochemical simulations with GIS and hydrogeochemical analyses has proven to be effective in identifying and locating the different physicochemical processes in the aquifer areas, thus improving understanding of hydrogeochemistry in complex aquifers.     

川东三叠系飞仙关组白云岩锶含量、锶同位素组成与白云石化流体

川东飞仙关组白云岩包括微晶白云岩、具原始结构的粒屑白云岩和结晶白云岩等3种主要的结构类型,其中结晶白云岩是重要的天然气储集岩。根据白云岩的锶舍量和锶同位素组成,结合相关岩石学和地球化学资料,本文研究了白云石化过程中锶的迁移和相应的同位素变化。三种类型白云岩具有类似的锶含量,但其(87)~Sr/(86)~Sr 比值均不同程度地高于同期海水值,锶在白云石和方解石(或文石)中分配系数的差别是影响白云石化过程中白云石锶含量最为重要的因素。不同流体中锶的混合方程计算表明:如果白云石化流体是大气淡水和海水的混合流体,则结晶白云岩的白云石化流体相当于94%的大气淡水和6%的海水混合,具原始结构的粒屑白云岩的白云石化流体相当于81%的大气淡水和19%的海水混合,微晶白云岩的白云石化流体相当于92%的大气淡水和8%的海水混合。但极低的锰含量、相对较高的锶含量和包裹体均一化温度等资料均不支持有94%的大气淡水参与了结晶白云岩的白云石化作用,其白云石化流体应是埋藏成岩过程中囚禁的近同期、但是非同期的海源流体,其来源可能与埋藏成岩过程中该时间段蒸发盐的溶解有关。锰含量、锶含量、(87)~Sr/(86)~Sr 比值都表明具原始结构的粒屑白云岩和微晶白云岩的白云石化过程都存在大气淡水的介入,微晶白云岩的白云石化可能与潮坪环境的蒸发泵机理有关;具原始结构的粒屑白云岩的白云石化可能与混合水白云石化作用有关,这类白云岩分布于向上浅滩化旋回滩体的顶部,这是埋藏前或浅埋藏阶段大气淡水与海水混合白云石化作用最容易发生地带。     

四川盆地川中地区震旦系灯影组碳酸盐沉积物成岩作用与孔隙流体演化

以四川盆地川中高石梯-磨溪地区钻井岩芯和水样资料为主,盆地周边地面剖面露头资料为辅,采用微量元素、阴极发光、包裹体、同位素等分析手段进行岩石学和地球化学分析,在此基础上讨论了灯影组碳酸盐岩的成岩作用变化和成岩作用过程中孔隙流体的演化。研究结果表明:灯影组碳酸盐岩沉积物在沉积作用和成岩作用过程中广泛地受到孔隙流体和成岩作用变化的影响,发生过海水的胶结作用、早期大气水成岩作用（溶解作用和沉淀作用）、表生喀斯特化作用、构造控制的热液白云石化和密西西比谷型矿化作用、碳氢产生时的溶解作用、油热裂解时焦沥青的充填孔洞作用等,其中喀斯特化和热液白云石化对储层的发育起着重要的作用。然而,灯影组孔隙流体广泛地被白云石化作用和喀斯特化作用的修饰。现今地层水资料分析表明这些成岩作用明显地影响蒸发海水的地层水流体。地层水分析资料在Ca-Br及Mg-Br关系图上显示出蒸发海水的特征,Ca的浓度比蒸发海水分布曲线要高,Mg的浓度比蒸发海水曲线要低,为灯影组沉积物在成岩作用过程中发生过白云石化作用的结果。     

Regional studies of dolomites and their included fluids: recognizing multiple chemically distinct fluids during the complex diagenetic history of Lower Carboniferous (Mississippian) rocks of the Irish Zn-Pb ore field

This study provides a regional framework within which studies of ore-related dolomite and dolomite cements may be placed. Fluid inclusion data indicate the presence of three distinct fluids following early dolomitization: 1) a ubiquitous low-temperature, higher salinity fluid found in saddle dolomite; 2) a low-temperature, lower salinity fluid limited to sub-Waulsortian and Waulsortian carbonates; and 3) a higher temperature, lower salinity fluid found in Waulsortian and supra-Waulsortian rocks. Similar fluids have been reported in ore-associated minerals and postmineralization dolomite (Type 1) and in ore-stage carbonates and sphalerite (Types 2 and 3). The halogen geochemistry of included fluids indicates genesis from evaporated seawater. Type 1 fluids are enriched in chloride relative to bromide, reflecting a component of salinity derived via dissolution of halite or from dehydration of seawater. These data suggest that dolomitization and mineralization of Mississippian rocks in the Midlands requires both regionally extensive and localized fluid flow.