一种基于环境风险规避的水源地优化选址方法

摘要

本发明公开了一种基于环境风险规避的水源地优化选址方法，属于环境科学与环境风险技术领域，包括以下三个步骤：步骤1：计算饮用水源地周边各个风险源大小S

说明书

技术领域

本发明属于环境科学与环境风险技术领域，具体涉及一种基于环境风险规避的水源地优化选址方法。

背景技术

水是生命之源，生产之要，生态之基。饮用水安全直接关系到公众健康，关系到经济社会的可持续发展，保障和提供安全饮用水已成为社会各界及各级政府关注的焦点。湖泊、河流是我国主要饮用水源地，维护河湖水质安全和水生态健康是保障饮用水安全的基础。

近年来，随着社会经济快速发展和人类活动强度加剧，水源地周边污染源密度逐渐上升，部分高危点源的污染负荷排放导致水源地达标率出现严重下滑。部分水源地岸线开发利用布局不合理，上下游岸线建有码头、化工企业、高新工业开发区，通过排污口或雨污混排口排放废水形成污染带，直接威胁水源地的供水安全，同时制约了局部优良岸线的开发利用。

目前，针对高强度岸线开发背景下饮用水源地的优化布局研究尚处于起步阶段，国内外对水源地的风险评估大多集中于计算污染源对水源地水质影响的直接风险，相关技术手段有待进一步提升。

因此，亟需提出一种针对高强度岸线利用下饮用水源地环境风险评估和规避的技术方法，为高风险水源地规避风险和优化选址提供技术依据，对于降低水源地环境风险和支撑水源地供水安全具有重要意义。

发明内容

发明目的：本发明的目的是提供一种基于环境风险规避的水源地优化选址方法，避免主观判断的不确定性，为保障水源地的供水安全提供技术支持。

技术方案：为了实现上述发明目的，本发明采用如下技术方案：

一种基于环境风险规避的水源地优化选址方法，包括以下三个步骤：

步骤1：计算饮用水源地周边各个风险源大小S_i、事故发生概率P_i，得到风险源风险值R_si；

步骤2：基于风险源风险值，结合事故危害程度，利用软件ArcGIS>

步骤3：划分饮用水源地风险等级，按照风险规避原则提出水源地优化选址方案或环境风险源调整方案。

进一步的，步骤1中，所述的风险源大小S_i获取方式包括如下步骤：

(1)收集水源地取水口及水源地保护区的地理信息；

(2)收集水源地保护区上下游10km范围内的排污口地理信息、行业类别、年均排污量、废水达标排放率；

(3)收集水源地保护区上下游10km范围内港口码头的地理信息、行业类别、年均吞吐量、岸线利用长度、危化品比例；

(4)收集水源地保护区上下游10km范围内入河支流的地理信息、年均流量、年均污染物浓度、入河污水排放达标率；

(5)将获取的数据进行标准化处理，采用层次分析法计算各风险源的大小S_i，如式(I)所示：

式(I)中，根据层次分析法确定权重ω_j，并根据各指标大小将各指标划分为若干级别，并对各级别进行赋值得到S_ij，进而得到风险源大小S_i(数据1)，如表5所示。式(I)中，n为饮用水源地保护区上下游10km范围内环境风险源的个数。

进一步的，步骤1中，所述的风险源事故发生概率P_i获取方式如式(II)所示：

式(II)中，P_ave为行业平均事故概率，P_s为不同类型风险源事故概率，P_a为不同受体事故概率，P_m为不同管理水平事故概率。

进一步的，步骤1中，各风险源风险值R_si获取方式如式(III)：

R_si＝S_i×P_i>

根据各个风险源大小与事故发生概率，得到该风险源的风险值。

进一步的，步骤2中，包括以下步骤：

(1)将水源地相关信息输入至软件ArcGIS 10.2中；

(2)将各个风险源地理信息及其风险值输入至软件ArcGIS 10.2中；

(3)采用反距离权重插值法，计算风险等值线分布图。

进一步的，步骤3中，包括以下步骤：

(1)根据风险等值线分布图，得到各个水源地的最大风险值；

(2)求得各个水源地风险值R_i的平均值μ与标准差s，如式(IV)和式(V)所示：

式(IV)和式(V)中，n为研究区域中饮用水源地的个数；

(3)判断水源地风险等级：当R_i≥μ+1.15s时，该水源地处于高风险；当μ+1.15s>R_i>μ－0.05s时，水源地处于中风险；当R_i≤μ－0.05s时，水源地处于低风险；

(4)水源地风险规避

基于经济性、有效性、合理性、可操作性并重的原则，针对高风险水源地，将饮用水源地从高风险区迁移或合并到低风险区，或通过关闭和迁移水源地周边环境风险源的方式降低饮用水源地的环境风险值，从而达到对高风险水源地进行风险规避的目的。

有益效果：与现有技术相比，本发明的一种基于环境风险规避的水源地优化选址方法，综合考虑了风险源大小、事故发生概率与危害程度，构建了水源地环境风险评价与风险规避方法，有利于避免主观判断的不确定性，为保障水源地的供水安全提供技术支持。

附图说明

图1为饮用水源地优化选址和风险规避方法示意图；

图2为长江南京段县级以上集中式饮用水源地和取水口GIS示意图；

图3为长江南京段县级以上集中式饮用水源地、取水口和环境风险源GIS示意图；

图4为环境风险源的GIS反距离权重插值示意图；

图5为长江南京段县级以上集中式饮用水源地的环境风险场等值线分布示意图；

图6为基于环境风险计算的饮用水源地优化选址技术流程示意图。

具体实施方式

下面结合附图和具体实施方式对本发明作进一步的说明。

目前国内外对水源地的风险评估大多集中于直接计算污染源对水源地风险，无法指导水源地选址风险规避，关于高强度岸线周边水源地周边选址的风险规避尚有待进一步完善。因此，本发明针对高强度岸线利用周边水源地风险的影响因素数量众多、类型复杂、差异明显的特点，提出了一种水源地选址的风险规避方法。本发明高强度岸线周边污染源与水源地出发，根据风险三要素(风险源大小、事故发生概率、事故危害程度)，构建了水源地环境风险评价与风险规避方法。

图6为基于环境风险计算的饮用水源地优化选址技术流程示意图。一种基于环境风险规避的水源地优化选址方法，包括以下三个步骤：

步骤1：计算饮用水源地周边各个风险源大小S_i(数据1)、事故发生概率P_i(数据2)，得到风险源风险值R_si(数据3)；

步骤2：基于风险源风险值R_si，考虑事故危害程度，利用软件ArcGIS>

步骤3：划分饮用水源地风险等级，按照风险规避原则提出水源地优化选址方案或环境风险源调整方案。

以下分别对上述三个步骤进行详细描述。

步骤1中风险源大小S_i(数据1)获取方式如下：

(1)收集水源地取水口及水源地保护区的地理信息，如表1所示；

(2)收集水源地保护区上下游10km范围内的排污口地理信息、行业类别、年均排污量、废水达标排放率，如表2所示；

(3)收集水源地保护区上下游10km范围内港口码头的地理信息、行业类别、年均吞吐量、岸线利用长度、危化品比例如表3所示；

(4)收集水源地保护区上下游10km范围内入河支流的地理信息、年均流量、年均污染物浓度、入河污水排放达标率如表4所示；

(5)将表1～表4中获取的数据进行标准化处理，采用层次分析法计算各风险源的大小S_i，如式(I)所示：

式(I)中，根据层次分析法确定权重ω_j，并根据各指标大小将各指标划分为若干级别，并对各级别进行赋值得到S_ij，进而得到风险源大小S_i(数据1)，如表5所示。式(I)中，n为饮用水源地保护区上下游10km范围内环境风险源的个数。

表1水源地信息

表2排污口及直排企业信息

编号名称经度纬度行业类别年均排污量废水排放达标率…………………

表3港口码头信息

表4入河支流信息

编号名称经度纬度年均流量(亿/m³)水质达标率………………

表5各风险源大小S_i

编号名称经度纬度风险源大小S_i……………

步骤1中风险源事故发生概率P_i(数据2)获取方式如式(II)所示：

式(II)中，P_ave为行业平均事故概率，P_s为不同类型风险源事故概率，P_a为不同受体事故概率，P_m为不同管理水平事故概率。上述各项概率的取值如表6所示。

表6风险源各项概率取值

步骤1中各风险源风险值R_si(数据3)获取方式如式(III)：

R_si＝S_i×P_i>

根据各个风险源大小与事故发生概率，可以得到该风险源的风险值，如表7所示。

表7各风险源风险值

编号名称经度纬度风险源大小事故发生概率风险源风险值…………………

所述步骤2主要包括以下内容：

(1)将水源地相关信息输入至软件ArcGIS 10.2中。

(2)将各个风险源地理信息及其风险值输入至软件ArcGIS 10.2中。

(3)采用反距离权重插值法，得到风险等值线分布图。

所述步骤3主要包括以下内容：

(1)根据风险等值线分布图，得到各个水源地的最大风险值。

(2)求得各个水源地风险值R_i的平均值(μ)与标准差(s)，如式(IV)和式(V)所示：

式(IV)和式(V)中，n为研究区域中饮用水源地的个数。

(3)判断水源地风险等级。

当R_i≥μ+1.15s时，该水源地处于高风险；当μ+1.15s>R_i>μ－0.05s时，水源地处于中风险；当R_i≤μ－0.05s时，水源地处于低风险。

(4)水源地风险规避。

基于经济性、有效性、合理性、可操作性并重的原则，针对高风险水源地，将饮用水源地从高风险区迁移或合并到低风险区，或通过关闭和迁移水源地周边环境风险源的方式降低饮用水源地的环境风险值，从而达到对高风险水源地进行风险规避的目的，如图1所示。

以长江江苏段27个集中式饮用水源地为实施案例，对本发明的技术方案进行阐述。

步骤1：计算饮用水源地周边各个风险源大小S_i(数据1)、事故发生概率P_i(数据2)，得到风险源风险值R_si(数据3)。

(1)收集水源地取水口和水源地二级保护区的地理信息，如表8。

表8水源地信息

(2)收集水源地周边排污口及直排企业的地理信息、行业类别、年均排污量、废水达标排放率，如表9。

表9排污口及直排企业信息

(3)收集水源地周边港口码头的地理信息、行业类别、年均吞吐量、岸线利用长度、危化品比例如表10。

表10港口码头信息

(4)收集水源地周边入河支流的地理信息、年均流量、年均污染物浓度、周边城镇入河污水排放达标率如表11。

表11入河支流信息

(5)用层次分析法计算各风险源的大小S_i。

根据层次分析法确定权重ω_j(如表12)，并根据各指标大小将各指标划分为若干级别，并对各级别进行赋值得到S_ij(如表13)，进而得到S_i数据1，如表14所示。

表12各指标权重ω_j

表13各风险源指标分级及其赋值S_ij

表14各风险源大小S_i

步骤1中风险源事故发生概率P_i(数据2)获取方式如下：

式中，P_ave为行业平均事故概率，P_s为不同类型风险源事故概率，P_a为不同受体事故概率，P_m为不同管理水平事故概率。各项概率的取值如表15所示。

表15风险源各项概率取值

步骤1中各风险源风险值R_si(数据3)获取方式如下：

R_si＝S_i×P_i>

根据各个风险源大小与事故发生概率，可以得到该风险源的风险值，如表16所示。

表16各风险源风险值R_si

步骤2：基于数据3，考虑事故危害程度，利用软件ArcGIS 10.2，采用反距离权重插值法，构建饮用水源地风险等值分布图。

(1)将水源地相关信息输入至软件ArcGIS 10.2中，如图2所示。

(2)将各个风险源地理信息及其风险值输入至软件ArcGIS 10.2中，如图3所示。

(3)采用反距离权重插值法，得到风险等值线分布图。

对风险值进行反距离权重插值，如图4所示。

进而得到风险等值线图，如图5所示。

步骤3：划分饮用水源地风险等级，按照风险规避原则提出水源地优化选址方案或环境风险源调整方案。

(1)根据风险等值线分布图，计算各个水源地的最大风险值R_i，如表17所示。

表17长江江苏段各水源地风险值R_i

水源地R_i水源地R_i水源地R_i子汇洲水源地1.253黄港水源地1.004长青沙水源地1.930夹江水源地1.133三江营水源地0.629张家港水源地1.385江浦水源地1.102二墩港水源地1.214狼山水源地0.975燕子矶水源地1.949永安洲水源地1.175洪港水源地0.954八卦洲水源地1.115魏村水源地0.918常熟水源地1.126龙潭水源地2.205西石桥水源地2.446海门水源地1.075仪征水源地0.918小湾水源地0.631浪港水源地0.959瓜洲水源地0.881肖山水源地0.636浏河水源地1.166征润州水源地0.993蟛蜞港水源地1.064

(2)根据式(IV)和式(V)，求得水源地风险值的平均值(μ)和标准差(s)分别为1.186和0.457。

(3)判断水源地风险等级。

当R_i≥1.71时，水源地为高风险；当1.71>R_i>1.16，水源地为中风险；当R_i≤1.16时，水源地为低风险。根据环境风险等级划分标准，燕子矶水源地、龙潭水源地、西石桥水源地、长青沙水源地为高风险；子汇洲水源地、二墩港水源地、永安洲水源地、张家港水源地和浏河水源地为中风险，其余16个水源地均为低风险。

(4)水源地风险规避。

燕子矶水源地、龙潭水源地、西石桥水源地、长青沙水源地为高风险饮用水源地。

通过迁移或合并水源地、迁移或关闭风险源等措施对高风险水源地进行风险规避，如表18所示。

表18高风险水源地风险规避对策

发明目的和意义：本发明描述了一种高强度岸线利用背景下基于风险规避的饮用水源地优化选址方法，全面考虑了污染源的风险大小与水源地的脆弱性，构建了由污染源大小、事故发生概率和危害程度三个要素组成的饮用水源地环境风险评价模型，提出了基于风险规避的水源地优化选址方法。本发明有利于实现对饮用水源地环境风险水平的全面评估，为降低饮用水源地环境风险、保障饮用水源地供水安全提供技术指导，具有重要的应用价值和指导意义。

本发明的描述和应用实例是说明性的，发明使用范围并非限制在上述实施例中。本发明专利使用人员及本领域技术人员须知晓并执行的是，如对本发明技术方案的细节或形式进行的合理修改或替换，均落入本发明的专利保护范围内。