冷文鹏，陶 亚，孙若涵，夏建新

(1. 中央民族大学 生命与环境科学学院，北京 100081；2. 生态环境部环境规划院，北京 100012)

我国饮用水源地面临水质污染超标、突发水污染事件频发等安全隐患，特别是具有水量丰富、水质比较稳定等特点[1]的地下水水源地被污染后不易修复，因此，要加强对地下水饮用水源地的保护。为保护地下水水质，建立地下水源地保护区是必不可少的措施[2]。但是，地下水流动比较复杂，地下水源地保护区划分存在一定的技术困难。利用数值模拟法确定地下水水源中的污染物运移方向和距离，可为保护区划分提供依据。邓媛媛等[3]利用FEFLOW 软件建立了吴忠市金积傍河型地下水水源地的三维地下水数值模型，将水源地的12 眼水井概化成“一眼、四眼、十二眼”，模拟了3 种地下水开采情景下地下水的运动，简化复杂的边界条件，并利用反向示踪粒子跟踪技术，建立了三级保护区。王涵等[4]在内蒙古呼和浩特市城市水源地运用MODFLOW 建立了该水源地的地下水流场，根据148 口抽水井单独设置了148 个反向示踪粒子点位，计算出反向示踪粒子的运移迹线，确定了水源地各级保护区的界线。肖杰等[5-9]运用Visual MODFLOW 中的modflow 和modpath 模块分别对崇州市地下水源地、渭河傍河水源地、关中盆地水源地进行了数值模拟，按照抽水井位置设置单个示踪粒子，对其进行保护区的划分。目前边界条件的相关研究方向集中在单井概化，缺少对井群概化的研究和应用。边界条件处理方法的差异可能导致模拟结果不同，按照单井概化处理的研究结果可能使保护区范围过小，带来污染风险。……