玛纳斯河流域山前平原地下水资源分析及合理开发利用研究

以玛纳斯河流域山前平原为例 ,论述我国西北内陆盆地地下水资源特点和组成 ,即地表水 -地下水联系密切、转化频繁。在已建立的研究区地下水流数值模拟模型的基础上 ,进行地下水均衡分析 ,得出研究区多年地下水补给资源量为 5 9816× 10 4 m3/a、可开采资源量为 396 90× 10 4 m3/a、地下水可利用量为 44 6 90× 10 4 m3/a。给出了不同水文地质分带的地表水 -地下水转化关系和转化量 ,预测需水量以每年 1.5 6 %～ 2 .17%的速度增长情况下 ,到 2 0 2 0洪积扇的中上部地下水位下降幅 2～ 5m ,而在溢出带以下下降小于 2m ,含水层的贮存量累计减少了6 85 0 0 .88× 10 4 m3。提出了地下水的合理开发利用模式 ,并将研究区分为 4个水资源开发利用分区 :冲洪积扇上部地表水利用区、中部地下水集中开采区、冲积扇缘井泉灌区和冲积平原井渠并灌区。     

高放废物深地质处置地下水数值模拟应用综述

开展地下水数值模拟研究是高放废物处置场地安全评价的重要组成部分,然而深地质处置介质类型的复杂性、基岩深部资料的相对匮乏性导致模拟结果存在不确定性,如何刻画深部地下水动力场并评估可能引起的风险已成为高放废物处置安全评价中重点关注的问题。在大量文献调研的基础上,综述了世界典型国家高放废物深地质处置场地的地下水数值模拟与不确定性分析应用,并归纳总结该领域研究经验,得到以下认识:（1）深地质处置场深部构造、裂隙的发育与展布决定了地下水循环条件,探究适用于基岩裂隙地区新的水文地质试验方法是提高地下水数值模型仿真性的基础;（2）不同尺度模型融合是解决深地质处置地下水模拟的有效技术方法,区域尺度多采用等效连续介质法,场地尺度使用等效连续多孔介质和离散裂隙网络耦合模型,处置库尺度使用离散裂隙网络方法,其次需重点关注未来大时间尺度下放射性核素在地质体中的迁移转化规律,模拟预测场址区域地下水环境长期循环演变对核素迁移的潜在影响;（3）考虑到不同的处置层主岩岩性以及在多介质中发生的THMC（温度场—渗流场—应力场—化学场）过程,目前国内外常用的地下水模拟软件有:Porflow、Modflow、GMS及MT3DMS等用于模拟孔隙或等效连续介质,Connectflow、Feflow及FracMan等用于模拟地下水和核素在结晶岩、花岗岩等裂隙中的迁移,TOUGH系列软件主要应用于双重介质的水流、溶质及热运移模拟;（4）指导开展有针对性的模型和参数的不确定性分析工作,减少投入工作量,提高模型精度,并可针对处置库长期演变、废物罐失效、极端降雨等多情景预测模拟,为处置库安全评价及设计提供基础数据支撑;（5）针对我国深地质处置地下水数值模拟研究现状,下一步应加强区域地质、水文地质、裂隙测量以及现场试验等相关的调查及监测工作,多介质耦合、多场耦合模拟及不确定性分析研究将会是未来的研究重点。     

基于数值模拟的岩溶地下水源保护区划分技术研究

以山东省邹城市某岩溶水源地为例,研究数值模拟方法在岩溶地下水水源地保护区划分中的应用,以及水文地质参数对于水源地保护区划分结果的影响。基于GMS软件构建了研究区三维非稳定流地下水数值模型,采用质点追踪技术,通过反向追踪示踪粒子在运移100 d、1 000 d以及25 a后的位置及迁移轨迹,对该岩溶水源地进行保护区划分;改变数值模型中的主要水文地质参数,以变化后的保护区面积作为衡量参数敏感度的指标,并计算各参数的敏感度系数,研究水文地质参数对于保护区划分范围的影响。结果表明:当参数改变幅度在20%范围内时,一层垂向渗透系数(VK1)的敏感度系数最大可达到2.63×10-3,二层垂向渗透系数(VK2)的敏感度系数最大可达到3.64×10-3;垂向渗透系数的敏感度明显高于其他参数,说明垂向渗透系数对保护区划分范围的影响最大。因此,在应用数值模拟法对岩溶地下水源地进行保护区划分时,应尤其注意模型中各含水层垂向渗透系数取值的准确性和合理性。      

基于多源数据的弱透水层水文地质参数反演研究

弱透水层是含水层系统的重要组成部分,其水文地质参数的准确获取一直以来都是研究难点,传统室内试验难以克服应力状态及环境变化对土层参数的影响,因此如何在原位状态下评价其水文地质参数有着重要意义。选取呼和浩特盆地某淤泥层为研究对象,基于抽水试验和溶质运移试验获得的多源数据,采用基于控制体积有限差分法的非结构化网格（Unstructured Grid）方法,建立三维地下水流-溶质运移耦合模型,对弱透水层水文地质参数进行反演。结果显示:模型计算的抽水井中水位和溶质浓度变化规律与实际观测数据能够较好拟合,模型识别出的抽水含水层渗透参数（4.8 m/d）与通过解析法计算得出的结果（4.17 m/d）相近,建立的数值模型符合实际水文地质条件;反演得到淤泥质弱透水层垂向渗透系数为1.2×10?4 m/d、储水率为1.0×10?5 m?1。基于多源数据的参数反演方法可为弱透水层参数研究提供借鉴。     

不同作物覆盖对农业区地下水入渗补给的影响分析

灌溉入渗是卫宁平原地下水的主要补给来源,研究不同农作物对地下水的灌溉入渗补给对于准确评价地下水资源量十分重要。本研究基于卫宁平原灌区2个包气带水分运移原位试验点系统观测数据,运用Hydrus-1D软件建立了试验点的包气带水分运移数值模型;设置了单次灌溉量、生长期天数、最大根系埋深和叶面积指数四种影响因子,应用模型分析了其对地下水入渗补给量的影响,计算了不同作物种植期内地下水垂向入渗补给量。结果表明:试验点地下水入渗补给量以灌溉入渗补给量为主。单次灌溉量的大小对地下水垂向入渗补给量的影响最为显著,其次是生长期天数和最大根系埋深,叶面积指数对地下水垂向入渗补给量的影响最小。随着地下水位埋深的增加,农作物种植因子对地下水入渗补给的影响也会增大。不同时期的降雨入渗系数为0.02～0.25;受次降雨量和降雨频率影响差异较大。灌溉入渗系数大小与作物种类关系密切:玉米种植期内的灌溉入渗系数为0.78,茄子种植期内的灌溉入渗系数为0.51,枸杞种植期内的灌溉入渗系数为0.6～0.63。综合考虑研究区作物类型和地下水位埋深(117～267 cm),给出了研究区农田区域在作物单次灌溉量为50～150 mm情况下,对应的灌溉入渗系数参考值。     

天津市地下水流-地面沉降耦合模型

天津市平原区地面沉降主要由地下水大量开采引起,影响范围广、危害大,已成为天津市主要的环境地质问题。分析了研究区的水文地质条件,结合地下水开发利用状况,将研究区概化为6个含水层组,地下水流考虑三维非稳定流,地面沉降选用一维固结压缩模型,运用地下水流模型Modflow 2005和地面沉降模拟模块 Sub,建立了天津市平原区地下水流-地面沉降数值耦合模型,模型面积为1.1×10⁴ km²,利用1998-2008年地下水位等值线、过程线、地面沉降过程线等资料对模型进行了识别。模拟期的地下水均衡分析表明,在多年开采条件下,越流补给、压缩释水、侧向边界流入分别占深层含水层补给量的41.84％、32.15％和24.17％。将调试后的模型应用于南水北调实施后地下水控采条件下的地面沉降趋势预测,显示出停采或减少地下水的开采,有利于减缓地面沉降下降速度,且表现出开采层位越往下,地面沉降恢复难度越大的变化趋势。     

鲁西南地区高氟水分布规律与成因分析

地方性氟中毒是我国北方地区最为典型的地方病之一,查明高氟地下水的空间分布及其成因是除氟改水、防治氟害的前提。通过对鲁西南地区不同层位地下水水质分析结果及水文地质、地质条件等多个环境影响因素的综合分析,查明了鲁西南地区高氟水的空间分布规律,并分析了影响浅层和深层高氟水形成的环境因素。浅层高氟水呈片状分布于洼地、缓平坡地等地势较低的区域,占鲁西南地区总面积的47%,大部分地区高氟水氟离子含量为1.2～2.0 mg/L,局部大于4.0 mg/L,其形成受气候、地质环境、地形地貌特征及水化学环境等多个因素的影响,成因类型为溶滤-蒸发浓缩型。深层高氟水具有水平分带性,占鲁西南地区总面积的65%,大部分高氟水氟离子含量为2.0～4.0 mg/L,氟离子含量分布与晚更新世沉积相带呈现很好的相似性,推测其为地质历史时期形成的古地下水。     

基于MODFLOW-CFP的贵州大井流域岩溶地下水数值模拟

贵州大井流域岩溶分布广泛，岩溶水是当地人民生产和生活的主要来源。由于对岩溶水资源的不合理开发利用，水资源短缺现象经常发生。大井流域水文地质条件复杂，管道-多孔介质双重介质特征明显。文章采用MODFLOW-CFP耦合模型对大井流域展开数值模拟，进而掌握大井流域地下水运动规律、准确评价岩溶水资源，促进其合理开发利用。结果表明：大井流域管道与多孔介质交换量为6 719.1 m³·a^-1，主要集中在上游和中游；总补给水量为10 977.3×10⁴ m³·a^-1，补给模数为133.495m³·km^-2·a^-1，其中降雨汇入量和降雨入渗量占总补给量的81.35%，而总排泄量为10 813.47×10⁴ m³·a^-1，主要在地下河出口排泄。