调水工程受水区水生态系统风险评估方法
【技术领域】
[0001] 本发明涉及生态风险评估领域，尤其涉及基于健康指数的调水工程受水区水生态 风险评估方法。
【背景技术】
[0002] 生态风险评价始于风险评价，始于20世纪70年代末80年代初的美国（Cairns， 1978)，最初的风险评价主要用于单一化学污染物对环境和人类健康影响的毒理研究。20世 纪80年代末，生态风险评价的工具和方法在一些研究中开始出现，但内容仍然侧重生物生 态毒理研究，尺度一般限于单一种群或者群落。20世纪80年代后，风险评价研究的内容开 始逐渐从毒理风险、人体健康风险向生态风险转变，尺度也从种群、群落向生态系统扩展。 很多学者开始把研究尺度扩展到了区域、景观和流域尺度，风险来源也由单一的化学污染 物延伸至化学、物理、生物等多领域的负荷风险源及可能造成生态风险的事件，并开始考虑 人类活动对生态系统的干扰（孙洪波、杨桂山、苏伟等，生态风险评价研究进展[J]。生态学 杂志，2009, 28 (2) :335-341)。调水工程作为人为干扰的重要表现形式，涵盖物理、化学、生 物、能量通量等多种风险源类型，由此带来诸多潜在水生态风险。但目前，风险评估工作集 中于工程风险和环境风险（刘恒，南水北调运行风险管理关键技术问题研究[M]，科学出版 社，2011 ;王发廷.，大型调水工程风险管理与保险[M]，黄河水利出版社，2011)，对受水区 水生态系统风险评估研究尚不多见，缺乏有效风险评估方法。
[0003] 当前制约调水工程受水区水生态系统风险评估方法发展的主要原因包括：1)水 生态系统风险的最终体现形式难以量化表征；2)水生态系统模拟预测工具的缺乏，使得调 水后水生态系统各组分的变化无法精确掌握，研究结论不足以支撑风险管理和决策。因此 在现有成果基础上，开发一种科学简洁的水生态系统风险评估方法，对调水工程受水区风 险评估和风险管理具有非常现实的意义。
【发明内容】
[0004] 本发明针对目前水生态系统风险评估方法欠缺问题，在充分考虑调水工程对受水 区水生态影响机制和风险特点的基础上，提出了基于健康指数的水生态风险评估方法。具 体流程见附图1。
[0005] 为实现上述目的，本发明包括如下技术方案：
[0006] 一种调水工程受水区水生态系统风险评估方法，包括如下步骤：
[0007] I.建立受水区水生态健康评价指标体系，该评价指标包括物理化学指标和生态指 标；
[0008] II.根据评估要求确定评价对象，该评价对象为时间段、空间方位或其结合；评价 对象的数量为一个或多个；
[0009] III.采集每个评价对象对应的评价指标现状数据；
[0010] IV.采用主成分分析法确定指标权重
[0011] (I)分别对各评价对象的每个评价指标现状数据进行归一化，将每个评价指标现 状数据中相对最佳的值归一化为1 ;其中针对越大越满意的指标采用公式（1)计算，针对越 小越满意指标采用公式（2)计算；
[0014] 其中，X1为某个评价对象中第i个评价指标的指标值；X为该指标的最佳值；
[0015] (II)计算各评价指标的权重值
[0016] 采用主成分分析法，应用SPSS软件来实现，输入全部评价对象的各项评价指标现 状数据，得到各评价指标对总体方差的贡献矩阵A = (X，y，z)，同时得到各评价指标的贡献 矩阵，即载荷矩阵：
[0018] 其中，η表示第η个评价指标，m表示第m个评价对象；
[0019] 求各评价指标的标准权重：W = AXBt= (f i，f2,…，fn);
[0020] 对标准权重进行归一化，得到各评价指标的权重W1;
[0021] V.计算现状受水区水生态健康指数
[0022] 根据公式（4)计算各评价对象的现状健康指数；
[0024] 其中，。表示该评价对象的健康指数；I i为该评价对象中第i种评价指标的归一 化值J1表示第i种评价指标的权重；
[0025] VI.受水区水生态模型建立与模拟
[0026] 应用MIKE软件模拟物理化学指标，应用AQUAT0X软件模拟水生态指标；对评价对 象进行模拟，具体步骤如下：
[0027] (I)MIKE 模型模拟
[0028] ①构建研究区地形文件；
[0029] ②输入模拟条件，包括边界条件和初始条件；确定模拟参数，在数据缺乏时模型参 数采用默认值；
[0030] ③根据评价时段和研究需要确定模拟时间和时间步长；
[0031] ④模型模拟及输出指标，输出指标包括水位、水深、流速以及调水后物理化学指 标；
[0032] (2) AQUAT0X 模型模拟
[0033] ①输入模拟条件，包括边界条件和初始条件，初始条件为步骤（I)MIKE模型输出 的水位、水深和流速指标；确定模拟参数，在数据缺乏时模型参数采用默认值；
[0034] ②根据评价时段和研究需要确定模拟时间和时间步长，步长应不小于MIKE模拟 步长；
[0035] ③1?拟及结果输出，1?拟结果为调水后生态指标；
[0036] VII.风险评估
[0037] 根据公式（1)、（2)计算调水后受水区各评价指标的归一化值，根据公式（4)计算 调水后受水区各评价对象的健康指数，根据公式（5)计算调水前后受水区健康指数变化程 度R;
[0039] 其中，Iatl为调水后健康指数；I ?为调水前健康指数。
[0040] VIII.确定风险等级
[0041] 根据R值确定风险等级：R多0.5,极尚；0.4<R<0.5，尚；0.2<R<0.4，中； 0.1 < R < 0. 2,低；0 < R < 0. 1，极低；R < 0,表明调水后水生态健康情况不变或有所改 善，不存在风险。
[0042] 如上所述的调水工程受水区水生态系统风险评估方法，优选地，所述物理化学指 标包括：溶解氧、高锰酸盐指数、总氮和总磷；所述生态指标包括：叶绿素 a、浮游植物生物 量、浮游动物生物量和浮游动值物生物量之比。
[0043] 如上所述的调水工程受水区水生态系统风险评估方法，优选地，所述步骤（I)MIKE 模型模拟过程中，边界条件包括：来水量和来水水质；初始条件包括：水位、流速及水质情 况；模拟参数包括：横向扩散系数、糙率、热传导系数、风速摩阻系数和水质降解系数。
[0044] 如上所述的调水工程受水区水生态系统风险评估方法，优选地，所述步骤（I)MIKE 模型模拟过程中，模拟条件包括：边界条件、外部环境条件和初始条件；其中边界条件包 括：来水量、来水水质和出水量；外部环境条件包括：风速、气温、降雨和蒸发；初始条件包 括：水位、流速及水质情况，现状水质情况包含评价指标现状数据。
[0045] 如上所述的调水工程受水区水生态系统风险评估方法，优选地，所述步骤（2) AQUAT0X模型模拟过程中，边界条件包括：来水量和来水水质；初始条件为步骤（I)MIKE模 型输出的水位、水深和流速指标；模拟参数包括藻类最大生长系数、最大死亡率和半饱和常 数。
[0046] 如上所述的调水工程受水区水生态系统风险评估方法，优选地，所述步骤（2) AQUAT0X模型模拟过程中，模拟条件包括：边界条件、外部及内部环境条件和初始条件；其 中边界条件包括：来水量、来水水质和出水量；外部及内部环境条件包括：风速、气温、降 雨、蒸发和水质；初始条件为步骤（I)MIKE模型输出的水位、水深和流速指标；模拟参数包 括：藻类最大生长系数、最大死亡率和半饱和常数。
[0047] 如上所述的调水工程受水区水生态系统风险评估方法，优选地，所述步骤（2) AQUAT0X模型模拟过程输出参数包括叶绿素 a、浮游植物生物量和浮游动物生物量；然后根 据浮游动植物生物量计算二者之比。
[0048] 表 1
[0049]
[0053] 本发明在对受水区特点分析的基础上，筛选能够反应水生态系统结构、功能和系 统特征的指标（如表1所示），并采用主成分分析法确定指标权重，建立健康评价体系，计算 现状健康指数值；耦合水动力和水生态模型，借助耦合的MIKE和AQUAT0X软件，对调水后受 水区水生态情况进行模拟，计算调水后健康指数值；通过指健康指数值变化程度确定风险 等级（如表2所示）。本发明的有益效果在于以下几个方面：
[0054] (1)科学：健康指数反映水生态系统整体健康情况已被科研界广泛接受，利用健 康指数的变化反映风险程度具有