DB44T749-2010饮用水水源保护区划分技术指引.pdf 18页

'ICS13.060.20DB44广东省地方标准DB44/T749-2010饮用水水源保护区划分技术指引TechnicalGuidelineforDelineatingSourceWaterProtectionAreas(发布稿)2010-05-14发布2010-07-14实施广东省环境保护厅发布广东省质量技术监督局 目次前言……………………………………………………………………………………………………II1适用范围………………………………………………………………………………………………12规范性引用文件……………………………………………………………………………………………13术语和定义………………………………………………………………………………………………14总则……………………………………………………………………………………………24.1饮用水水源保护区的设置与划分……………………………………………………………………24.2划分的一般技术原则………………………………………………………………………………35技术指标……………………………………………………………………………………………45.1保护区分级………………………………………………………………………………………45.2水质要求……………………………………………………………………………………………46河流型饮用水水源保护区的划分方法………………………………………………………………46.1一级保护区………………………………………………………………………………………46.2二级保护区………………………………………………………………………………………56.3准保护区……………………………………………………………………………………………67水库型饮用水水源保护区的划分方法……………………………………………………………67.1一级保护区………………………………………………………………………………………67.2二级保护区………………………………………………………………………………………77.3准保护区……………………………………………………………………………………………78地下水饮用水水源保护区的划分方法………………………………………………………………89其他……………………………………………………………………………………………………810饮用水水源保护区的最终定界……………………………………………………………………811监督实施……………………………………………………………………………………………8附录Ａ编写技术文件的基本要求（规范性附录）……………………………………………………9附录Ｂ水质模型基本方程及解析解（资料性附录）……………………………………………10I DB44/T749-2010前言为贯彻落实《中华人民共和国水污染防治法》、《广东省饮用水源水质保护条例》等相关的法律法规，加强饮用水水源地环境保护、防治饮用水水源地污染，改善环境质量，按照《标准化工作导则》（GB/T1.1-2009）制定本标准。本标准在《饮用水水源保护区划分技术规范》（HJ/T338－2007）技术要求的基础上，结合广东省实际情况，进一步细化了地表水饮用水水源保护区划分的基本方法及饮用水水源保护区划分技术文件的编制要求。本标准为首次发布。本标准由广东省环境保护厅提出。本标准起草单位：广东省环境监测中心。本标准起草人：易雯岑世柏刘乙敏陈小文陈松明温丽蓉本标准广东省环境保护厅和广东省质量技术监督局2010年5月14日批准。本标准自2010年7月14日起实施。本标准由广东省环境保护厅解释。II DB44/T749-2010饮用水水源保护区划分技术指引1适用范围本标准适用于集中式地表水、地下水饮用水水源保护区（包括现用、备用和规划水源地）的划分。2规范性引用文件下列文件对于本文件的应用是必不可少的。凡是注明日期的引用文件，仅注日期的版本适用于本文件。凡是不注明日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB3838-2002地表水环境质量标准GB5749-2006生活饮用水卫生标准GB15618-1997土壤环境质量标准GB/T14848-1993地下水环境质量标准3术语和定义下列术语和定义适用于本标准。3.1饮用水水源地Drinkingwatersource指提供居民生活及公共服务用水取水工程的水源地域，包括河流、湖泊、水库、地下水等一定区域。3.2集中式饮用水水源地Centralizeddrinkingwatersource指集中取水进入输水管网送到用户和具备一定规模(供水人口一般大于1000人或日供水3量大于200m)的现用、备用和规划饮用水水源地。3.3界河Sharedjurisdictionriver指流经两个互不隶属行政区之间、边界线所经过的河流，也称共河。3.4饮用水水源保护区Drinkingwatersourceprotectionraeas指国家为防治饮用水水源地污染、保证水源地环境质量而划定，并要求加以特殊保护的一定面积的水域和陆域。3.5山区型水库Mountainousreservoir1 DB44/T749-2010指用拦河坝横断河谷，拦截河川径流，抬高水位形成的水库。3.6平原型水库Plainreseroir指在平原地区，利用天然湖泊、洼淀、河道，通过修筑围堤和控制闸等建筑物形成的水库。3.7山脊线Ridgeline山脊是由两个坡向相反、坡度不一的斜坡相遇，组合而成条形脊状延伸的凸形地貌形态。山脊最高点的连线即两个斜坡的交线，叫做山脊线。3.8分水岭Watershed分隔相邻两个流域的地带，可以是山地、高原或是微有起伏的山丘、平原。3.9中泓线Midstreamofchannel河道各横断面表面最大流速点的连线。3.10库容系数Storagecapacityindex水库兴利库容与多年平均来水量的比值。4总则4.1饮用水水源保护区的设置与划分4.1.1饮用水水源保护区分为地表水饮用水水源保护区和地下水饮用水水源保护区，地表水饮用水水源保护区包括一定面积的水域和陆域，地下水饮用水水源保护区指地下水饮用水水源地的地表区域。4.1.2集中式饮用水水源地都应设置饮用水水源保护区；饮用水水源保护区一般划分为一级保护区和二级保护区。饮用水水源地二级保护区外围一定距离内分布有对水环境可能造成重大影响的建设项目或污染源密集区，存在较大污染风险时，宜设置准保护区，并将其纳入准保护区进行管理，其划分范围可参照二级保护区的划分方法确定，但不宜采用类比经验方法确定。4.1.3饮用水水源保护区的设置应纳入当地社会经济发展规划和水污染防治规划；跨地区的饮用水水源保护区（含界河）的设置应纳入有关流域、区域、城市社会经济发展规划和水污染防治规划。4.1.4在水环境功能区划和水功能区划分中，应将饮用水水源保护区的设置和划分放在最优先位置；跨地区的（含界河）河流、水库、输水渠道，不得影响上、下游（或相邻）地区饮用水水源保护区对水质的要求，并应保证下游有合理水资源量；对同一流域上下游相II DB44/T749-2010距10公里以内且分处不同城市的两个水源地，两城市共同协调划定保护区，或由两城市共同的上级政府划定保护区。4.1.5跨地区（含界河）和向异地供水的水源地水源保护区的划分按照有关法律法规的规定，由用水城市商水源地所在城市提出划分方案，由地级以上市人民政府报省人民政府批准；存在异议时，可由省人民政府环境保护行政主管部门会同有关行政主管部门和相关市、县人民政府提出划定方案。4.1.6应对现有集中式饮用水水源地进行评价和筛选；对于因受污染而达不到饮用水源水质要求且经技术、经济论证证明饮用水功能难以恢复的水源地，应采取措施调整饮用水水源地。4.1.7饮用水水源保护区的水环境监测与污染源监督应作为监督管理工作的重点，纳入地方环境管理体系中，若不能满足保护区规定水质的要求，则应及时调整保护区范围。4.2划分的一般技术原则4.2.1水源地分类4.2.1.1地表水水源地4.2.1.1.1河流型水源地非潮汐河段水源地：不受潮汐影响。潮汐河段水源地：受潮汐影响。拦河闸(坝)型水源地：用拦河闸(坝)横断河流，抬高水位形成的河流型水源地。4.2.1.1.2水库型水源地依据水库规模大小，将水库型饮用水水源地进行分类。3大型水库：水库总库容≥1亿m3中型水库：0.1亿m3≤水库总库容＜1亿m3小型水库：水库总库容＜0.1亿m4.2.1.2地下水水源地地下水按含水层介质类型的不同分为孔隙水、基岩裂隙水和岩溶水三类；按地下水埋藏条件分为潜水和承压水两类。地下水饮用水水源地按开采规模分为中小型水源地（日开采量小于5万立方米）和大型水源地（日开采量大于等于5万立方米）。4.2.2饮用水水源保护区划分的技术指标，应考虑以下因素：当地的地理位置、水文、气象、集水汇流特征、地质特征、土地利用、水动力特性、水域污染类型、污染特征、污染源分布、排水区分布、水源地规模、水量需求等。地表水饮用水水源保护区范围，应按照不同水域特点进行水质预测，并考虑当地具体条件，保证在规划设计的水文条件、污染负荷和供水量时，保护区的水质能满足相应的标准。地下水饮用水水源保护区范围，应根据当地的水文地质条件、供水量、开采方式和污染源分布确定，并保证开采规划水量时能达到所要求的水质标准。4.2.3划定的水源保护区范围，应防止水源地附近人类活动对水源的直接污染；应足以使所选定的主要污染物（含影响当地人群健康的特征污染物）在向取水点（或开采井、井群）输移（或运移）过程中，衰减到所期望的浓度水平；在正常情况下保证取水水质达到3 DB44/T749-2010规定要求；一旦出现污染水源的突发事件，有采取紧急补救措施的充足时间和缓冲地带。4.2.4对高强度开发区域，为加强区域环境管理，水源保护区划分范围宜适当扩大。4.2.5划分的技术方法可采用水质模型或类比经验方法。但由于各地自然条件不同，为了保证计算的科学合理性，各地可根据当地的水文地质特征选用合适的模型进行计算。同时开展跟踪验证监测。5技术指标5.1保护区分级5.1.1一级保护区：在取水口附近划定一定范围的水域和陆域，以保证取水口水质安全。5.1.2二级保护区：在一级保护区外围划定一定范围的水域和陆域，以保证在正常情况下满足水质要求，同时在出现污染饮用水源的突发情况下，保证有足够的采取紧急措施的时间和缓冲地带。5.1.3准保护区：在二级保护区外围划分一定范围的水域和陆域，以防范二级保护区上游高危项目或污染源较集中区域的污染风险，有效控制污染，并保证有足够的采取紧急措施的时间和缓冲地带。5.2水质要求5.2.1地表水饮用水水源保护区水质要求5.2.1.1地表水饮用水源一级保护区的水质必须满足GB3838中的Ⅱ类标准要求，且补充项目和特定项目应满足该标准规定的限值要求。5.2.1.2地表水饮用水源二级保护区的水质必须满足GB3838中的Ⅲ类标准要求，并保证流入一级保护区的水质满足一级保护区水质标准的要求。5.2.1.3地表水饮用水源准保护区的水质应保证流入二级保护区的水质满足二级保护区水质标准的要求。5.2.2地下水饮用水源保护区水质要求地下水饮用水水源保护区（包括一级、二级和准保护区）水质各项指标不得低于GB/T14848中的Ⅲ类标准。6河流型饮用水水源保护区的划分方法6.1一级保护区6.1.1水域范围6.1.1.1通过分析计算，确定一级保护区水域长度。6.1.1.1.1非潮汐河段水源地，应用二维水质模型计算得到一级保护区范围，一级保护区水域长度范围内应满足GB3838Ⅱ类水质标准的要求。二维水质模型及其解析解参见附录Ｂ，大型、边界条件复杂的水体采用数值解方法，对小型、边界条件简单的水体可采用解析解方法进行模拟计算。6.1.1.1.2潮汐河段水源地，运用非稳态水动力－水质模型模拟，计算可能影响水源地II DB44/T749-2010水质的最大范围。非稳态水动力－水质模型参见附录Ｂ。16.1.1.1.3一级保护区上、下游范围不得小于卫生部门规定的饮用水源卫生防护带范围。6.1.1.2在水文、源强、边界条件等参数无法确定的情况下，可采用类比经验方法确定一级保护区水域范围，同时开展跟踪监测。若发现划分结果不合理，应及时予以调整。河流型饮用水水源地一级保护区水域长度划分见表1。表1河流型饮用水水源地一级保护区水域长度划分(类比经验方法)河流型污染源取水口上游2000米距取水口上游距取水口下游水期水源地类型类型河段多年平均流速水域长度水域长度＜1.5米/秒≥1500米面源丰水期≥1.5米/秒，＜2.5米/秒≥2000米≥2.5米/秒≥2500米非潮汐河段≥100米＜0.8米/秒≥1500米点源枯水期≥0.8米/秒，＜1.5米/秒≥2000米≥1.5米/秒≥2500米潮汐河段≥1500米≥1500米6.1.1.3一级保护区水域宽度没有堤防的河流型水源地，为5年一遇洪水所能淹没的区域；有堤防的为至堤防内侧的水域范围；通航河道在多年平均水位时，水面宽度达500米以上的河流，以河道中泓线为界，保留一定宽度的航道外，规定的航道边界线到取水口范围即为一级保护区范围。6.1.2陆域范围一级保护区陆域范围的确定，以确保一级保护区水域水质为目标，采用以下分析比较确定。6.1.2.1陆域沿岸长度不小于相应的一级保护区水域长度。6.1.2.2河流两岸为浅滩、平原、小山丘的水源地其陆域沿岸纵深与河岸的水平距离不小于50米，若有防洪堤则为至堤外侧的距离；同时，一级保护区陆域沿岸纵深不得小于饮用水水源卫生防护范围。6.1.2.3河流两岸为陡峭山峰的水源地其一级保护区陆域范围为沿岸侧纵深至第一重山山脊线。6.2二级保护区6.2.1水域范围6.2.1.1通过分析计算，确定二级保护区水域范围6.2.1.1.1二级保护区水域范围应用二维水质模型计算得到，二级保护区上游侧边界到一级保护区上游边界的距离应大于污染物从GB3838Ⅲ类水质标准浓度水平衰减到GB3838Ⅱ类水质标准浓度所需的距离。二维水质模型及其解析解参见附录Ｂ，大型、边界条件复杂的水体采用数值解方法，对小型、边界条件简单的水体可采用解析解方法进行模拟计算。6.2.1.1.2潮汐河段水源地，二级保护区采用模型计算方法；按照排污对取水口影响的程度,计算确定二级保护区上、下游边界位置。1卫监发[2001]161号文生活饮用水集中式供水单位卫生规范5 DB44/T749-20106.2.1.2在水文、源强、边界条件等参数无法确定的情况下，可采用类比经验方法确定二级保护区水域范围。河流型饮用水水源地二级保护区水域长度划分见表2。表2河流型饮用水水源地二级保护区水域长度划分(类比经验方法)河流型污染源距一级保护区下游侧外边界距一水期该河段多年平均流速水源地类型类型上游边界级保护区边界向上游延伸＜1.5米/秒≥2500米向上游延伸面源丰水期≥1.5米/秒，＜2.5米/秒≥3000米向上游延伸≥2.5米/秒非潮汐≥3500米≥200米河段向上游延伸＜0.8米/秒≥2500米向上游延伸点源枯水期≥0.8米/秒，＜1.5米/秒≥3000米向上游延伸≥1.5米/秒≥3500米潮汐河段不宜采用类比经验方法确定6.2.1.3二级保护区水域宽度：一级保护区水域向外10年一遇洪水所能淹没的区域，有防洪堤的河段二级保护区的水域宽度为防洪堤内的水域宽度。6.2.2陆域范围二级保护区陆域范围的确定，以确保水源保护区水体水质为目标，采用以下分析比较确定。6.2.2.1二级保护区陆域沿岸长度不小于一级和二级水域保护区河长。6.2.2.2二级保护区沿岸纵深范围自一级保护区陆域和二级保护区水域沿岸向外不小于1000米，具体可依据自然地理、环境特征和环境管理需要，通过分析地形、植被、土地利用、地面径流的集水汇流特性、集水域范围等确定；两岸是陡峭山峰的河流型水源地，其第一重山山脊线高于50年一遇的洪水线时可不划二级保护区；有物理隔离区的、封闭输水河（渠）的水源地可视具体情况不划或适当划分二级陆域保护区；对于流域面积小于100平方公里的小型流域，二级保护区可以是整个集水范围。6.3准保护区准保护区划分可参照二级保护区的划分方法确定范围。7水库型饮用水水源保护区的划分方法7.1一级保护区7.1.1水域范围7.1.1.1小型水库正常水位线以下的全部水域面积划为一级保护区。7.1.1.2大中型水库采用模型分析计算方法确定一级保护区范围。7.1.1.2.1当大、中型水库部分水域面积划定为一级保护区时，应对水体进行水动力（流动、扩散）特性和水质状况的分析、二维水质模型模拟计算，确定水源保护区水域面积，即一级保护区范围内主要污染物浓度满足GB3838Ⅱ类水质标准的要求。具体方法参见附录Ｂ,宜采用数值计算方法。II DB44/T749-20107.1.1.2.2一级保护区范围不得小于卫生部门规定的饮用水源卫生防护范围。7.1.1.3在水文、源强、边界条件等参数无法确定的情况下，采用类比经验方法确定一级保护区水域范围，同时开展跟踪验证监测。若发现划分结果不合理，应及时予以调整。7.1.1.3.1中型水库水域范围为不小于取水口半径300米范围内的区域。7.1.1.3.2大型水库为不小于取水口半径500米范围内的区域。7.1.2陆域范围水库沿岸陆域一级保护区范围，以确保水源保护区水体水质为目标，采用以下分析比较确定：7.1.2.1小型水库至流域分水岭，大、中水库为一级水域保护区沿岸正常水位线以上200米范围内的陆域或至流域分水岭。7.1.2.2一级保护区陆域沿岸纵深范围不得小于饮用水水源卫生防护范围。7.2二级保护区7.2.1水域范围7.2.1.1通过模型分析计算方法，确定二级保护区范围。二级保护区边界至一级保护区的径向距离大于所选定的主要污染物或水质指标从GB3838Ⅲ类水质标准浓度水平衰减到GB3838Ⅱ类水质标准浓度所需的距离，具体方法参见附录Ｂ，宜采用数值计算方法。7.2.1.2在技术条件有限的情况下，采用类比经验方法确定二级保护区水域范围。7.2.1.2.1中型水库一级保护区边界外的水域面积设定为二级保护区。7.2.1.2.2大型水库以一级保护区外径向距离不小于2000米区域为二级保护区水域面积，但不超过水面范围。7.2.2陆域范围二级保护区陆域范围确定，应重点考虑面污染源，要依据流域内自然地理、环境特征和环境管理的需要，通过分析地形、植被、土地利用、森林开发、地面径流的集水汇流特性、集水域范围等确定。二级保护区陆域边界不超过相应的流域分水岭范围。7.2.2.1小型水库可将上游整个流域（一级保护区陆域外区域）设定为二级保护区。山区型小型水库可不设二级保护区7.2.2.2平原型中型水库的二级保护区范围是正常水位线以上（一级保护区以外），水平距离2000米区域，山区型中型水库二级保护区的范围为水库周边第一重山山脊线以内（一级保护区以外）及入库河流上溯3000米的汇水区域。7.2.2.3大型水库可以划定一级保护区外不小于3000米的区域为二级保护区范围。7.3准保护区根据水库流域范围、高危项目污染源分布，在二级保护区以外的汇水区域设定准保护区。8地下水饮用水水源保护区的划分方法地下水饮用水水源保护区的划分，应在收集相关的水文地质勘查、长期动态观测、水源地开采现状、规划及周边污染源等资料的基础上综合确定。广东省地下水饮用水水源保护区划分方法采用《饮用水水源保护区划分技术规范》7 DB44/T749-2010HJ/T338-2007中地下水饮用水水源保护区的划分方法。9其他类型饮用水水源保护区的划分方法9.1如果饮用水源一级保护区或二级保护区内有支流汇入，应从支流汇入口向上游延伸一定距离，作为相应的一级保护区或二级保护区，划分方法可参照上述河流型保护区划分方法划定。根据支流汇入口所在的保护区级别高低和距取水口距离的远近，其范围可适当减小。9.2非压力式饮用水输水渠道，应划一级保护区，其范围可比河流型保护区适当减小。有自然径流汇入的饮用水输水渠道的水源保护区划分，可参照河流型的划分方法。9.3以水库为主要补给源的河流型水源地，其饮用水水源保护区范围应包括水库一定范围内的水域和陆域，保护级别按具体情况参照水库型水源地的划分办法确定。9.4对非主要利用潮汐动力进行水体交换的河口拦河闸（坝）型水源地，应同时满足水库型饮用水水源保护区和潮汐河段型饮用水水源保护区的划分要求。对主要利用潮汐动力进行水体交换的河口拦河闸（坝）型水源地参考潮汐河段型饮用水水源保护区的划分。9.5对于库容系数大于10%的具有调节功能的拦河闸（坝）型水源地水源保护区的划分参照水库型饮用水水源保护区的划分方法，其他类型的拦河闸(坝)型水源地水源保护区的划分参照河流型保护区的划分，闸坝下游可不划保护区。9.6入库河流（引水渠）的保护区水域和陆域范围的确定，以确保水库饮用水水源保护区水质为目标，参照河流型饮用水水源保护区的划分方法确定一、二级保护区的范围。10饮用水水源保护区的最终定界10.1为便于开展日常环境管理工作，依据保护区划分的分析、计算结果，结合水源保护区的地形、地标、地物特点，最终确定各级保护区的界线。10.2充分利用具有永久性的明显标志如水分线、行政区界线、公路、铁路、桥梁、大型建筑物、水库大坝、水工建筑物、河流汊口、输电线、通讯线等标示保护区界线，并应设置专门标志。10.3最终确定的各级保护区坐标红线图、表，作为政府部门审批的依据，也作为规划国土、环保部门土地开发审批的依据。II DB44/T749-2010附录Ａ编写技术文件的基本要求（规范性附录）A．1“XXXX饮用水水源保护区划分技术报告”格式划分饮用水水源保护区，应编写正式的“XXXX饮用水水源保护区划分技术报告”技术文件。技术文件的基本内容应包括以下几个部分：A．1.1划分依据A．1.1.1相关法律法规；A．1.1.2相关已经批准实施的规划。A．1.2保护区背景分析1.2.1饮用水水源保护区所在区域或流域的自然状况；1.2.2饮用水水源保护区所在区域或流域的社会经济状况；1.2.3饮用水水源地的资源、环境质量评价。评价的基本内容包括水量、水质状况及发展趋势，可能对水源地产生污染影响的主要污染源、污染物及污染影响途径，作为饮用水源开采的前景；与相邻水体的关系，包括饮用水水源保护区上、下游或相邻水域（或相邻区域）的水体功能、保护区的水量和水质是否受本行政区外的影响；若受到其影响，列出影响途径、影响程度（水量、水质、生态、经济、人体健康等）等实测数据、定量计算和定性分析结果；1.2.4划分饮用水水源保护区的必要性和充分性说明。A．1.3技术方法与计算结果1.3.1根据各级保护区的划分方法，说明选用的技术指标、数值计算方法；1.3.2计算结果及分析，各级保护区定界的技术说明；1.3.3用图表示各级保护区的范围，并用表格确定红线坐标，保护区内污染源、集水区、排水区分布特性等。A．1.4饮用水水源保护区的监督与管理措施饮用水水源保护区内的水质监测网站的布置，水质项目的监测，陆源污染的监督等；若水质尚未达标，应确定水质达标期限和相应的管理与控制措施。A．1.5饮用水水源保护区划分方案、图件及有关说明饮用水水源保护区划分方案的说明，表明保护区详细情况（包括监测点的位置等）的图集、饮用水水源保护区登记表、保护区详细情况的文字说明，准保护区划分的必要性及意义等。A.2饮用水水源保护区调整的可行性研究技术报告可参照A1。9 DB44/T749-2010附录Ｂ水动力和水质模型基本方程及解析解（资料性附录）二维水质模型的基本方程为：22∂C∂C∂C∂C∂C=D+D−u−u−KC…………………（附1）x2y2xy∂t∂x∂y∂x∂y∂C在稳态条件下，=0，上式可变形为：∂t22∂C∂C∂C∂CD+D−u−u−KC=0………………………(附2)x2y2xy∂x∂y∂x∂y对于应用于水质模拟的二维模型，会涉及到有无边界影响两类情况。B.1无边界水体边界点源的稳态排放在均匀流场中，当强度为M的点源排放到无限宽的水体中，见图B.1：Ox在边界条件为：∂C=0时，式（附2）的解析解为：∂yy=0图B.1宽度无限水体中的点源排放M⎛(y−uyxux)⎞⎛x⎞C()x,y=exp⎜−⎟exp⎜−K⎟………（附3）4πh()xu2DD⎜4Dxu⎟⎜u⎟xxy⎝yx⎠⎝x⎠式中：uy——y方向的流速分量；Dy——y方向的扩散系数；H——平均水深；3K——污染物的降解速率（m/s）；虚源B源B实源OxO宽度无限宽度有限图B.2污染物的边界排放II DB44/T749-2010如果是顺直河道，在水深变化不大的情况下横向流速很小，近似为零；纵向扩散项远小于推流的影响，即可以忽略uy和Dx项，则式（附2）可简化为：2∂C∂CD−u−KC=0……………………………………………（附4）y2x∂y∂x相应的解析解为：M⎛uy2⎞⎛x⎞C()x,y=exp⎜−x⎟exp⎜−K⎟……………（附5）uh4πDx/u⎜4Dx⎟⎜u⎟xyx⎝y⎠⎝x⎠B.2有边界水体连续点源的稳态排放在有边界的情况下，污染物的扩散会因受到边界的阻碍而产生反射，这种反射可以通过设立虚源来模拟，即设想边界为一面镜子，镜子后面有一个与实际源强度相同，距离相同的虚拟反射源。当有两个边界时，反射会成为连锁式的。当污染源在边界上，对于宽度无限大的环境，有：⎛2⎞2M⎜uxy⎟⎛x⎞C()x,y=exp−exp⎜−K⎟………………（附6）uh4πDx/u⎜4Dx⎟⎜u⎟xyx⎝y⎠⎝x⎠可以看出，对于全反射的边界（不考虑扩散物质被边界吸附），污染物的浓度是没有反射时的两倍。对于宽度为B的环境，则：2M⎧⎪+∞⎛u()2nB−y⎞⎫⎪⎛x⎞C()x,y=⎨∑exp⎜−x⎟⎬exp⎜−K⎟………（附7）uh4πDx/u⎪⎩n=−∞⎜4Dyx⎟⎪⎭⎜⎝ux⎟⎠xyx⎝⎠当污染源在两个边界的中间时，有：M⎧⎪+∞⎛u()nB−y⎞⎫⎪⎛x⎞C()x,y=⎨∑exp⎜−x⎟⎬exp⎜−K⎟………（附8）uh4πDx/u⎪⎩n=−∞⎜4Dyx⎟⎪⎭⎜⎝ux⎟⎠xyx⎝⎠y虚源边界的反射的影响随着距离的增加(n的增大)而衰减很B/2快，当n>4以后，计算结果基本趋于稳定，计算时取n=4～5就足够了。B实源x如果污染源的位置既不位于边界，也不位于河流正中央，而是位于距岸y0（0≤y0≤B）的位置，即可以表达为：OB/2虚源图B.3双边界的中心排放M⎧⎪+∞⎛u(y−(2nB±y))2⎞⎫⎪⎛x⎞C()x,y=⎨∑exp⎜−x0⎟⎬exp⎜−K⎟……（附9）uh4πDx/u⎪⎩n=−∞⎜4Dyx⎟⎪⎭⎜⎝ux⎟⎠xyx⎝⎠瞬时点源排放时，无边界阻碍的情况下，边界条件为：∂Cy=±∞,=0∂y11 DB44/T749-2010时，其解析解为：22M⎛⎜()x−uxt(y−uyt)⎞⎟⎛x⎞C()x,y,t=exp−−exp⎜−K⎟………（附10）4uhDDt2⎜4Dt4Dt⎟⎜u⎟xxy⎝xy⎠⎝x⎠有边界阻碍时，可将上式修正为：MC()x,y,t=24uhDDtxxy⎧2222⎫⎪⎛⎜()x−uxt()y−uyt⎞⎟⎛⎜()x−uxt()2b+y−uyt⎞⎟⎪⎛x⎞⎨exp−−+exp−−⎬exp⎜⎜−K⎟⎟⎜4Dt4Dt⎟⎜4Dt4Dt⎟u⎪⎩⎝xy⎠⎝xy⎠⎪⎭⎝x⎠………（附11）式中，b－污染源到边界的距离当为岸边排放时，即b=0时，上式可变为：222M⎡()x−uxt(y−uyt)⎤C()x,y,t=exp⎢−−−⎥exp(−Kt)………（附12）4uhDDt2⎢4Dt4Dt⎥xxy⎣xy⎦B.3感潮河段保护区范围计算由于河道的横向宽度与河口地区相较窄，水质要素（主要指污染物浓度）在河道横向上的差异并不显著，可以认为各水力要素在河段断面上是一致的。因此选择一维水动力模型和一维对流扩散方程建立感潮河网一维数学模型。一维河网水动力模型基本方程表达如下：连续性方程：1∂Q∂Z+=qL………………………………………………(附13)B∂x∂t动力学方程：∂u∂u∂Zuu+u+g+g=0…………………………………(附14)2∂t∂x∂xCRSQ式中，Z表示断面水位；Q表示流量；u=表示平均流速；A为河道过水断面面积；gA为重力加速度；B表示不同水位下的水面宽度；qL为旁侧入流流量；R是水力半径；CS是谢才(Chezy)系数；x、t是位置和时间座标。在汊点处的衔接条件目前使用的最广泛的是水流连续和能量守恒两个条件。其表述如下：（1）流量连接条件进出每一节点的流量必须与该节点内实际水量的增减率相平衡，即：II DB44/T749-2010j+1j∂WZ−Z∑Qi==A………………………………………(附15)∂tΔt式中，Qi表示通过i断面进入节点的流量；W为节点的蓄水量；A为调蓄节点的蓄水面积j+1j(汇合区面积)；Z、Z为调蓄节点j+1时刻和j时刻的水位，若调蓄节点面积很小，则：∑Qi=0………………………………………………………(附16)（2）动力连接条件如果节点可以概化成一个几何点，出入各节点的水位平缓，不存在水位突变情况，则各节点相连汊道的水位应相等，等于该点的平均水位，即：Z=Z=L=Z…………………………………………………(附17)120如果各断面的过水面积相差悬殊，流速有较明显的差别，当略去汊点的局部损耗时，得伯努力(Bernouli)方程如下：22uu12Z+=Z+………………………………………………(附18)122g2g它表明各支汊河道上在汊点处流体微段的动力学特性趋于一致。采用四点隐式差分格式，差分格点的布置如图B4所示n+1（1-θ）△tq△tθ△t△xinii+1图B4一维河网水动力模型使用的四点偏心格式n+1n+11(−θ)nnθnnq=θ(q+q)+(q+q)−(q+q)………………………(附19)ii+1i+1ii+1i22∂qθn+1n1(−θ)nn=(q−q)+(q−q)………………………………………(附20)i+1ii+1i∂xΔxΔxii∂q1n+1n+1nn=(q+q−q−q)………………………………………………(附21)i+1ii+1i∂t2Δt∂qθnn+1nn+11(−θ)nnnnθnnnnq=(qq−qq)+(q+q)(q−q)−(q+q)(q−q)i+1i+1iii+1ii+1ii+1ii+1i∂xΔx2ΔxΔxiii………………………………………………(附22)13 DB44/T749-2010式中，q为变量，可代表流量、水位等，将四点隐式差分格式代入连续性方程(附13)和运动学方程(附14)，经整理后得差分方程为：H−AQ+H+AQ=A…………………………………(附23)i2ii+12i+15−BH+BQ+BH+BQ=B…………………………(附24)1i2i3i+14i+15其中：Δt11A=θ(+)2ΔxBBiii+1nΔt11nnnnA=H−1(−θ)(+)(Q−Q)+H+q5ii+1ii+1LΔxBBiii+12uuB=gθΔt{+()}122iΔxCRi⎡Δt2u⎤1B2=⎢1−2θui+gθΔt(2)i⎥ΔxCRS⎣i⎦⎡2uu⎤B3=gθΔt⎢−(22)i+1⎥ΔxCR⎣i⎦⎡Δt2u⎤1B4=⎢1+2θui+1+gθΔt(2)i+1⎥ΔxCRS⎣i⎦nnΔtnnnn2ΔtnnB=u+u−1(−θ)(u+u)(u−u)−g1(−θ)(H−H)−5ii+1i+1ii+1ii+1iΔxΔxii⎡uunuun⎤uug1(−θ)Δt⎢(2)i+(2)i+1⎥−gθΔt(22)iHi−⎣CRCR⎦CRuuΔt22uuuugθΔt()H−θ(u−u)+gθΔt()+gθΔt()22i+1i+1ii+12i2i+1CRΔxCRCRi…………………………(附25)式中：S为过水面积。对河网非恒定流基本方程组采用一种按隐式计算的三级联合解法，即将参加计算的方程分成微段、河段、汊点三级，逐级处理，再联合运算，求得河网中各微段断面的水位、流量等值。河网水质模型基本方程表达如下：∂(AC)∂(QC)∂∂C+−(AE)+S−S=0………………（附26）XC∂t∂X∂X∂XNLdz∑(QC)l,j=(CΩ)J()j………………………………（附27）l=1dt以上式（附26）是河道方程，式（附27）是河道交汊点方程。式中Z，Q是水位及流量；II DB44/T749-2010A是河道断面面积；EX是纵向分散系数；C是水流输送的物质浓度；Ω是河道交叉点—汊点的水面面积；j是汊点编号；l是与汊点j相联接的河道编号；SC是与输送物质浓度有关的衰减项，在这里写成SC=KdACC；Kd是衰减因子；S是外部的源汇项。河网水质模型求解骤(1)根据河道的流态，利用式（附28）～(附31)，建立每一条河道上各断面的递推方aC+bC+cC=z程组ii−1iiii+1i(i=1,…,n)（附28）式中ai，bi，ci是系数；Ci是i断面时段末的浓度；n是某一河道的断面数。对于一般断面（i=1,…,n）有⎧ai=−[]rc1D11+rd1D21+Fc1Δt/V⎪⎪bi=[]rc1D11+rc2D22+rd1D21+rd2D32+Fc2−Fd2Δt/V+[]rc1Kd,i−1+rd2Kd,iΔt+.1⎨c=−[]rD+rD−FΔt/V⎪ic222d232d3⎪k[]z=αC+rSΔx+rSΔxΔt/V⎩iiic1i−1i−1d2ii…………………………(附29)对于首、末断面（i=1或n）有⎧a1=0⎪⎪b1=[]rd2D32−Fd2Δt/V2+rd2Kd1,Δt+.1⎨…………………………(附30)c=−[]rD−FΔt/V⎪1d232d32⎪kz=αC+rSΔxΔt/V⎩111d2112以及⎧an=−[]rc1D11+Fc1Δt/V1⎪⎪bn=[]rc1D11+Fc2Δt/V1+rc1Kd,n−1Δt+.1⎨………………………（附31）c=0⎪n⎪kz=αC+rSΔxΔt/V⎩nnnc1n−1n−11(2)建立汊点浓度方程组(3)根据汊点浓度方程组，求得河网中每个汊点的浓度值(4)将汊点浓度值回代给与各汊点相连的河段首、末断面浓度未知量。15'