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

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

水力屏障和截渗墙在海水入侵防治中的数值模拟研究

基于SEAWAT-2000程序构建室内二维砂箱试验中咸水入侵数值模型,利用该模型分析了针对不同补给井井位、补给井流量、截渗墙位置及贯穿深度等多种情景下的咸淡水界面运移规律。二维砂箱实验模拟结果表明,当注水井位于盐水楔前锋附近,距咸水边界40 cm、砂箱顶部边界5 cm处时,注水井工程措施能达到最佳海水入侵驱退效果,回退系数达21.5%。当截渗墙布设于距咸水边界10 cm处,贯穿深度为35 cm时,截渗墙工程措施能达到最佳海水入侵驱退效果,回退系数达81.1%。在此基础上,结合实际场地条件,构建山东龙口地区滨海含水层中某典型二维剖面的海水入侵数值模型,探讨了不同截渗墙布设情景模式下海水入侵状况。模拟结果表明,当截渗墙布设于距海岸线600 m处,贯穿深度为18 m时,截渗墙工程措施达到最佳海水入侵驱退效果,回退系数达28.4%。研究结果揭示了补给井井位、补给井流量、截渗墙位置及贯穿深度等因素对咸淡水界面运移规律的影响,可为场地条件下滨海含水层海水入侵防治中的工程管理措施优化提供参考依据。     

内陆单斜构造内咸水入侵淡水含水层三维数值模拟：——以山西柳林？…

本文建立了内陆单斜构造内咸水入侵淡水含水层的三维可混溶溶质运移模型。同时,为了尽可能提高模拟精度,还建立了区域二维水流模型和水质模型,为三维模型提供边界条件。模型应用于描述柳林电厂水源地咸－淡水界面运移特征,模拟效果良好,还据此预报了咸水入侵的可能性并探讨了淡水流的强弱对咸－淡水界面运移的影响。     

内陆单斜构造内咸水入侵淡水含水层三维数值模拟--以山西柳林泉区柳林电厂水源地为例

本文建立了内陆单斜构造内咸水入侵淡水含水层的三维可混溶溶质运移模型.同时,为了尽可能提高模拟精度,还建立了区域二维水流模型和水质模型,为三维模型提供边界条件.模型应用于描述柳林电厂水源地咸-淡水界面运移特征,模拟效果良好,还据此预报了咸水入侵的可能性并探讨了淡水流的强弱对咸-淡水界面运移的影响.     

电阻率法在防治莱州湾地区海水入侵中的应用

阐述电阻率法在莱州湾地区海水入侵防治研究中应用的理论依据、监测指标.对海、咸水入侵区地层电阻率的变化规律,咸、淡水界面的确定,圈定海水入侵通道及界面运移规律做了分析.     

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

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

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

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

人工回灌条件下滨海砂质含水层渗透性的时空演变特征

人工回灌补给滨海含水层是海水入侵修复经常采用的工程措施,然而人工回灌条件下滨海砂质含水层经常出现渗透性显著降低的现象,其含水介质渗透性的时空演变规律仍不清楚。本文以青岛市大沽河下游咸水入侵区含水层砂样为研究对象,通过利用不同尺度的室内砂柱淡咸水驱替试验,对定水头和定流速条件下人工回灌咸淡水驱替过程中砂质含水层渗透性的时空演变特征进行了研究。结果表明：当淡水驱替咸水时含水介质的渗透性会发生显著变化,整个砂柱的渗透性先降低后回升,原因是回灌过程中砂柱中的黏土矿物发生释放、迁移、沉积;定水头条件下,在砂柱的前半段存在一个淘空区,渗透系数是初始渗透系数的1.6~2.0倍,砂柱其它部分渗透性是先降低再轻微回弹;定流速条件下,砂柱的渗透系数随时间的变化均呈现先降低后升高趋势,且砂柱各段离入水口的距离越远其渗透系数的值越小;砂柱黏粒含量的变化规律与渗透性的演化特征相吻合。研究成果可为人工回灌治理滨海含水层海水入侵提供一定的科学依据。     

海水入侵研究现状与展望

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

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

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

Modeling the dynamics of the freshwater-saltwater interface in response to construction activities at a coastal site

A numerical model was developed to evaluate the response of groundwater flow and the fresh-saltwater interface in relation to the construction of a particle accelerator at the coastal plain of Tokaimura, Japan. Undisturbed conditions were initially simulated and validated against field observations as a prerequisite for the analysis of predictive scenarios. Groundwater heads and the shape of the saltwater interface were appropriately described by the model, although it tended to underestimate salinity concentrations. Saltwater penetrated up to 250 m inland during predevelopment conditions, reaching more than 400 m at the dewatering phase. Flushing of entrapped saline groundwater might occur in addition to seawater intrusion. In depth, multiple saltwater fronts develop in response to the hydraulic properties of the sediments. Groundwater discharges offshore through the sandy aquifers, but salinity fronts prevail in the relatively impermeable layers. Routes for freshwater outflow turned into pathways of seawater intrusion during the pumping phase. The equilibrium would be reestablished within 2 years from the end of the stress, with no evidence of a permanent deterioration of neighbor residential wells. Nonetheless, after construction the accelerator forms a barrier that leads to a sharp rise in piezometric levels and creates a new and long-term disequilibrium in the saltwater wedge. Despite further work is still necessary to test many of the ideas proposed, the present study makes a new contribution to enhance the understanding of the processes occurring in coastal aquifers subjected to anthropogenic influence.     

Rapid and Intense Phosphate Desorption Kinetics When Saltwater Intrudes into Carbonate Rock

It is important to understand how phosphate sorption dynamics of coastal carbonate aquifers are affected by seawater intrusion, because many coastal aquifers are composed of carbonate rocks and subject to an increase in saltwater intrusion during relative sea-level rise. Twelve carbonate rock and unconsolidated sediment specimens were acquired from a test corehole spanning the full thickness of the Biscayne aquifer in southeastern Florida. All 12 samples exhibit low phosphorus content but variable contents of iron. Column leaching experiments were conducted with two carbonate aquifer samples, alternating between freshwater and saltwater flow. With the first influx of saltwater, phosphate concentration in leachate increased rapidly from a freshwater value of approximately 0.2 μM to peaks of between 0.8 and 1.6 μM. The phosphate concentration began to diminish as saltwater continued to flow, but sustained desorption continued for over 2 h. Overall, seawater drove sorption behavior much more than chemical composition for the aquifer rocks and sediment from the seven rock samples for which we did isotherm sorption experiments. Our results indicate that an immediate and intense pulse of phosphate desorption from carbonate rock and sediment with low phosphorus content occurs in response to an influx of seawater and that the duration of desorption will vary by layer within a single aquifer.     

Finite element methods for variable density flow and solute transport

Saltwater intrusion into coastal freshwater aquifers is an ongoing problem that will continue to impact coastal freshwater resources as coastal populations increase. To effectively model saltwater intrusion, the impacts of increased salt content on fluid density must be accounted for to properly model saltwater/freshwater transition zones and sharp interfaces. We present a model for variable density fluid flow and solute transport where a conforming finite element method discretization with a locally conservative velocity post-processing method is used for the flow model and the transport equation is discretized using a variational multiscale stabilized conforming finite element method. This formulation provides a consistent velocity and performs well even in advection-dominated problems that can occur in saltwater intrusion modeling. The physical model is presented as well as the formulation of the numerical model and solution methods. The model is tested against several 2-D and 3-D numerical and experimental benchmark problems, and the results are presented to verify the code.     

Saltwater intrusion in coastal regions of North America

Saltwater has intruded into many of the coastal aquifers of the United States, Mexico, and Canada, but the extent of saltwater intrusion varies widely among localities and hydrogeologic settings. In many instances, the area contaminated by saltwater is limited to small parts of an aquifer and to specific wells and has had little or no effect on overall groundwater supplies; in other instances, saltwater contamination is of regional extent and has resulted in the closure of many groundwater supply wells. The variability of hydrogeologic settings, three-dimensional distribution of saline water, and history of groundwater withdrawals and freshwater drainage has resulted in a variety of modes of saltwater intrusion into coastal aquifers. These include lateral intrusion from the ocean; upward intrusion from deeper, more saline zones of a groundwater system; and downward intrusion from coastal waters. Saltwater contamination also has occurred along open boreholes and within abandoned, improperly constructed, or corroded wells that provide pathways for vertical migration across interconnected aquifers. Communities within the coastal regions of North America are taking actions to manage and prevent saltwater intrusion to ensure a sustainable source of groundwater for the future. These actions can be grouped broadly into scientific monitoring and assessment, engineering techniques, and regulatory approaches.     

Saline fresh water interface structure in Mahanadi delta region, Orissa, India

 Saline/fresh water interface structure is one of the most important and basic hydrogeological parameter that needs to be estimated for studies related to coastal zone management, well-field design and understanding saline water intrusion mechanism/processes. The success and stability of a groundwater structure in a coastal region depend upon an accurate estimate of interface structure between saline and fresh water zones, aquifer-aquiclude boundaries and their lateral continuities and the interstitial water qualities of aquifers. Self-potential and resistivity logs provide a reasonably good basis for such estimates and for sustainable development of fresh groundwater resources. The interface depth structure for the Mahanadi delta region, as obtained and interpreted through self-potential and resistivity logs, provides a fairly clear picture of the regional extensions and boundaries of aquifers, aquicludes and interstitial water quality patterns. Aquifers in the northern sector of the basin and within the framework of Birupa and Mahanadi are characterized by an interface depth range that varies between 40 and 280 m below ground level (bgl) with brackish water on the top underlain by freshwater aquifers. The aquifers in the southern sector within the framework of Khatjori/Devi and Koyakhai/Daya/Kushbhadra/Bhargavi are characterized by an interface depth range that varies from 10 to 120 m with freshwater aquifers near the surface underlain by saline, brackish water aquifers. The inversion of these major fluid systems appears to have taken place over a narrow zone between Mahanadi and Khatjori tributaries, possibly over a wide subsurface ridge with separate basin characteristics. Received: 29 November 1999 · Accepted: 2 May 2000     

咸淡水界面位置确定的综合方法(TEcG)及其应用

准确确定咸淡水界面位置是评价咸水入侵范围的前提。对于咸淡水界面位置的确定,传统研究以水化学法和物探法为主。为克服单一使用传统方法造成的人力、物力和财力的大量浪费,以莱州湾西南岸广饶县小清河以南咸水入侵区为研究对象,在综合分析研究区地下水开采现状、水化学监测结果的基础上,选择3个典型断面,采用水化学法中的野外现场电导率法,快速判断咸淡水界面大致位置（某两眼监测井之间）;再在一咸一淡的两眼监测井之间,采用高密度电阻率法快速、准确地确定咸淡水界面的空间分布。结果表明：该地区的地下水电导率若大于1.61 mS/cm,即可认为此处受到咸水入侵;咸淡水界面位置的视电阻率特征值为11~13 Ω·m,咸淡水界面附近咸水体呈舌状入侵并主要发生在地表以下13 m内的浅层地下水中。     

Geochemical processes in the saltwater–freshwater transition zone: comparing results of a sand tank experiment with field data

Geochemical processes, occurring in a stable transition zone between saltwater and freshwater, were simulated in a 2D, multi-layer flow chamber experiment. Mixing, calcite dissolution, and oxidative degradation of organic matter were identified as the main controlling factors. The results of the chamber experiment were compared to field data and verified by thermodynamic modeling. Similarity in most ion distributions suggests the general applicability of the experimental method. Differences in the redox conditions between field and experiment were reflected by the oxidants involved in the mineralization of organic carbon; while field data show evidence of sulfate reduction, the presence of oxygen in the laboratory experiment resulted in the reoxidation of sulfides.     

Dynamic of the freshwater–saltwater interface in a karstic aquifer under extraordinary recharge action: the Merida Yucatan case study

A hydraulic analysis of the interface between freshwater–saltwater behavior was done in the Merida Yucatan zone, two machines that constantly register the groundwater levels were installed, and three electric conductivity logs were taken from wells. When comparing the measured results with the ones obtained using theoretical equations developed to calculate the freshwater–saltwater interface position, it was proved that in some cases these equations can be applied, and in others not. Two effects that rule the behavior of karst aquifers in extraordinary conditions were found.     

Marine electrical resistivity imaging of submarine groundwater discharge: sensitivity analysis and application in Waquoit Bay, Massachusetts, USA

Electrical resistivity imaging has been used in coastal settings to characterize fresh submarine groundwater discharge and the position of the freshwater/salt-water interface because of the relation of bulk electrical conductivity to pore-fluid conductivity, which in turn is a function of salinity. Interpretation of tomograms for hydrologic processes is complicated by inversion artifacts, uncertainty associated with survey geometry limitations, measurement errors, and choice of regularization method. Variation of seawater over tidal cycles poses unique challenges for inversion. The capabilities and limitations of resistivity imaging are presented for characterizing the distribution of freshwater and saltwater beneath a beach. The experimental results provide new insight into fresh submarine groundwater discharge at Waquoit Bay National Estuarine Research Reserve, East Falmouth, Massachusetts (USA). Tomograms from the experimental data indicate that fresh submarine groundwater discharge may shut down at high tide, whereas temperature data indicate that the discharge continues throughout the tidal cycle. Sensitivity analysis and synthetic modeling provide insight into resolving power in the presence of a time-varying saline water layer. In general, vertical electrodes and cross-hole measurements improve the inversion results regardless of the tidal level, whereas the resolution of surface arrays is more sensitive to time-varying saline water layer.