基于生态系统完整性的河流水生态健康评价方法及其应用与流程

本发明涉及河流生态系统完整性评价的技术研究领域，涉及一种基于鱼类生物完整性且包括水体物理完整性、化学完整性、生境完整性在内的水生态健康评价方法及其应用。
背景技术：
：随着社会的发展，人类对河流水资源的利用率越来越高，河流生态系统稳定性遭到严重破坏。水质污染、栖息地环境破坏及三生物多样性丧失等诸多问题频繁出现。针对这一现象，越来越多的国内外学者开始关注河流健康并开展对其进行健康评价的研究。河流评价始于20世纪中期，目前关于河流生态系统健康评价的方法有很多，最初主要是对水体生物需氧量、化学需氧量以及大肠菌群数量等的测定，但这种方法较片面，不能准确反映河流的总体状态。后来，karr等于1981年提出了生物完整性的定义，建立了鱼类生物完整性指数(fishindexofbiologicalintegrity,f-ibi)，应用生物信息对河流生态系统完整性进行评估。1999年，ladson提出溪流状态指数（indexofstreamcondition,isc），从水文、物理、河岸带状况、水质、水生生物特征5个方面来评价澳大利亚河流的健康状况。wrightjf等于2010年提出o/e指数，运用河流无脊椎动物预测与分类系统（rivpacs）模型计算物种观测值(o)和期望值(e)的比值，从而反映所调查点位的健康状况。2010年stevensonrj等以硅藻为研究对象提出单因子评价，选择5个评价指标分别对河流健康进行评价。由于河流的生物完整性并不等同于河流的生态系统健康，生物信息只能反映河流生物类群的状况，而无法准确指示水生态环境质量的优劣，使得目前普遍使用的应用生物信息指示河流健康状况的方法存在一定的缺陷。比如，鱼类是最先作为指示生物评价河流的健康状况，其生物等级高，在食物链中的地位和对水质的灵敏程度高，能够综合表征河流水生态健康状况，是评价河流健康较适宜的指标。但近年来，中度干扰假说愈发受到人类的关注，在中等程度干扰的情况下，鱼类生物量反而是较多的，因此，仅以生物群落作为河流健康状况的评价依据，结果往往有失偏颇。技术实现要素：本发明要解决的技术问题，是针对目前应用生物完整性评价不能准确指示河流健康状况的现状，提供一种基于生态系统完整性的河流水生态健康评价方法及其应用。本发明的目的通过以下技术方案予以实现：本发明以淡水生态系统中的鱼类生物完整性指标体系及其水环境物理完整性指标体系、水环境化学完整性指标体系及水环境生境完整性指标体系为监测对象，根据各指标的变化，表征人为干扰作用下的水质污染和生境退化对河流生态系统的影响情况，为水生态系统健康评价体系提供一套有效的技术体系。所述方法包括三个关键步骤：s1.基于生态系统结构特征，整合生态系统全部生物要素和环境要素，通过遥感-实测-人工解译-智能分析模型，实现河流生态系统健康评价的全面性、准确性及可靠性；s2.四大体系权重比；s3.评价标准等级化；其中，步骤s1所述基于生态系统结构特征，整合生态系统全部要素（生物和环境），通过遥感-实测-人工解译-智能分析模型，实现河流生态系统健康评价的全面性、准确性及可靠性；其中，所述基于河流生态完整性，整合生态系统全部要素是综合鱼类生物完整性、水环境物理完整性、水环境化学完整性及水环境生境完整性四部分指标体系，反应河流生态系统生物、物理、化学及生境的不同层面。优选地，所述四大体系指标类型具体包括50个候选指标，所述50个候选指标为：鱼类完整性指标体系包括雅罗鱼亚科个体百分比（%）、鳑鲏亚科个体百分比（%）、鮈亚科个体百分比（%）、鲤亚科个体百分比（%）、鳅科鱼类个体百分比（%）、鰕虎鱼科鱼类百分比（%）、本地特有鱼种（葛氏鲈塘鳢）百分比（%）、经济鱼类个体百分比（%）、肉食性鱼类数量比例（%）、植食性鱼类数量比例（%）、杂食性鱼类的数量比例（%）、敏感性物种百分比（%）、耐污物种百分比（%）、中上层鱼类个体百分比（%）、中下层鱼类个体百分比（%）、底层鱼类个体百分比（%）、冷水鱼百分比（%）、有护卵行为鱼类个体百分比（%）共20个候选指标；物理完整性指标体系包括ph、电导率（cond）、盐度（s）、水温（t）、溶解氧（do）、总溶解固体（tds）及悬浮物（ses）共7个候选指标；化学完整性指标体系包括氯离子（cl-）、氨氮（nh4+-n）、总氮（tn）、总磷（tp）、活性磷（po43-）、硬度（har）、高锰酸盐指数（codmn）、亚硝酸态氮（no2-n）、硝酸态氮（no3-n）、总有机碳（toc），总无机碳（tic）及总溶解碳（tdc）共12个候选指标；生境完整性指标体系包括流速（cv）、流量（flo）、水深（wd）、河宽（w）、底质含沙量（sc）及土地利用方式含缓冲区流域面积（tot）、森林面积（fa）、草地面积（ga）、建设用地面积（cla）、水田面积（pfa）和旱田面积（fie）共11个候选指标；本发明确定的指标体系能够综合反映鱼类生物群落、水体理化性质及生境的变化对河流生态系统的影响，可以确保有的效监测水生态健康功能状况。其中，通过遥感-实测-人工解译-智能分析模型，实现河流生态系统健康评价的全面性、准确性及可靠性；优选地，所述遥感-实测-人工解译-智能分析模型具体是：土地利用方式的获取是基于遥感技术（rs），结合上游区域的土地利用特征，以左右岸800米为界与土地利用数据叠加分析获得；生物、物理及生境完整性指标体系是通过现场实测与实验室实测相结合的方法获取；经过人工解译初步筛选指标，结合软件canoco4.5进行主成分分析（pca）以及软件spss2.0进行pearson相关性分析筛选核心指标，利用标准化方法并进行数据智能分析。步骤s1所述智能分析方法包括主成分分析（pca）、pearson相关性分析以及标准化方法；其中，所述主成分分析（pca）是选取轴1或轴2的贡献率p>0.6的指标；优选地，所述主成分分析（pca）是保留轴1或轴2特征贡献率p>0.6的指标，即保留对评价结果贡献率较大的指标，删除轴1和轴2的特征贡献率均小于0.6的指标；其中，pearson相关性分析是依照相关性显著性水平判断生物间的信息重叠程度，如果指标间的信息重叠程度较高，则删除其中一个指标；优选地，依据相关性显著性水平，判断信息间的重叠程度；生物完整性的相关系数为|r|<0.9，物理完整性、化学完整性、生境完整性的相关系数为|r|<0.7，如果超过r值，则保留其中一个指标即可。其中，标准化方法分为两种方法，目的是为了统一量纲；优选地，所述标准化方法分为两种：对于随干扰增强而数值减小的指标，以其95th分位数为期望值记分；对于随干扰增强而数值增加的指标，以其5th分位数为期望值记分。所述计算指标分值的计算公式分别为：fij=xij/xi95×100fij=（xmax-xij）/（xmax-xi5）×100式中，fij为原始观测值；xij为j个采样点的第i个指标的标准化指数；xi95为全部采样点中第i个指标95th分位数的原始观测值；xmax为全部采样点中的最大值；xi5为全部采样点中第i个指标5th分位数的原始观测值。步骤s2所述四大体系权重比运用层次分析法对四部分指标体系进行权重，赋予鱼类生物完整性、物理完整性、化学完整性、生境完整性分别为40%、20%、20%、20%的权重，并进行加和平均得出生态系统完整性的河流水生态健康评价的得分。步骤s3所述评价标准等级化是以大于90th分位数作为健康的标准；将小于90th分位数的部分进行三等分，即60th~90th为较好的标准；30th~60th为一般的标准；小于30th为较差的标准，最终得出水生态健康评价的结果。本发明提供的所述方法在淡水生态系统评估方面的应用，应用本发明方法，可以对河流、湖泊为代表的内陆淡水生态系统进行科学、综合的水体质量监测和生态环境评价。本发明的有益效果在于：生态系统是一个完整的体系，包括生物因素和非生物因素，本发明以河流鱼类群落完整性为基础，结合水环境物理完整性、水环境化学完整性及水环境生境完整性对河流健康状况进行评价，因此该方法是一种更为有效和可靠的健康评价方法。本发明建立基于鱼类生物完整性且包括水体物理、化学、生境完整性的水生态健康评价方法既包括与周围环境关系密切的、受干扰后反应敏感的多个生物指标又包括对环境变化响应较快的理化指标，即从多方位、多角度对河流健康状况进行了检测。本发明评价体系有效弥补现有传统的水体理化指标只能反映瞬时水质状况的缺点和技术不足；同时对生物完整性健康评价进行了补充，能够体现人为干扰对生境破坏的累积效果，预测多重污染胁迫下的有害生物效应，是一种经济、定量、智能、多指标体系的综合手段，为评价各个监测点位及整个流域的生态健康状况提供有力的技术支撑。本发明方法具有指标全面、适用面广、操作简单等优点，在水生态评价
技术领域：
具有重要推广应用价值。附图说明图1是本发明方法流程示意图；图2是鱼类生物完整性主成分分析图；图3是物理完整性主成分分析图；图4是化学完整性主成分分析图；图5是生境完整性主成分分析图；图6是水生态健康评价结果分布图。具体实施方式下面结合附图和具体实施例进一步说明本发明方法。下述实施例和附图仅用于示例性说明，不能理解为对本发明的限制。除非特别说明，下述实施例中选用的评价河流为随机选取，并不因此限定本发明应用范围和对象。本发明方法中使用的试剂原料为常规或商业途径获得的生试剂原料，除非特别说明，下述实施例中使用的方法和设备为本领域常规使用的方法和设备。实施例1对浑太河流域生态系统进行水生态健康评价按照本发明方法，选定浑太河流域生态系统进行健康评价本发明所述基于生态系统完整性包括鱼类生物完整性、物理完整性、化学完整性及生境完整性的水生态健康评价方法，方法流程（步骤）如图1所示，其体系构建过程及健康评价标准如下：s1.基于河流生态完整性，整合生态系统全部要素，通过遥感-实测-人工解译-智能分析模型，实现河流生态系统健康评价的全面性、准确性及可靠性；（1）对浑太河流域共设置65个采样点，太子河41个点位（t1-t36，t39-t43），浑河25个点位(h1-h16，h18-h25)。对鱼类生物完整性指标体系进行采集，分类与鉴定，本次调查共采集鱼类11920尾，分为10科33属40种，其中，鲤科最多，占物种数的47.5%，其次是鳅科、鰕虎鱼科，它们分别占物种数的12.5%和7.5%。浑太河流域的常见种是北方条鳅(nemachilusnudusbleeker)、洛氏鱥（phoxinuslagowskiidybowski）、棒花鱼(gobiogobiorivuloidesnichols)。ph、电导率（cond）、盐度（s）、水温（t）、溶解氧（do）、总溶解固体（tds）采用水质分析仪（ysi85）现场获取数据；流速（cv）、流量（flo）、水深（wd）采用流速仪(fp111)获取数据；河宽（w）使用测距仪(leupoldｒx-iv)获取数据；其他指标根据gb3838—2002《地表水环境质量标准》于实验室内测定。土地利用方式的获取基于遥感技术（rs），根据研究区的实际情况以及参考2007年《土地利用现状分类》，结合上游区域的土地利用特征，以左右岸800米为界与土地利用数据叠加分析获得。选用鱼类群落候选指标20个、物理候选指标7个、化学候选指标12个、生境候选指标11个，以上候选指标的变化能够对河流生态系统的健康状况产生影响，可以确保有效的监测水生态健康功能状况。（2）对候选指标进行主成分分析（pca），选取轴1或轴2的贡献率p>0.6的指标，删除轴1和轴2的特征贡献率p＜0.6的指标。在鱼类生物完整性的pca分析中（图2），鱼类总分类单元数和总渔获量两个指标的贡献率p>0.6，其余十八个指标的贡献率p＜0.6，均不符合条件；在物理完整性的pca分析中（图3），电导率、盐度及总溶解固体三个指标的贡献率p>0.6，其余四个指标的贡献率p＜0.6；在化学完整性的pca分析中（图4），氨氮、硝酸态氮、总氮及总溶解碳四个指标的贡献率p>0.6，其余指标均不符合条件；在生境完整性的pca分析中（图5），缓冲区流域面积、森林面积、建设用地面积、水田面积、旱田面积、水深、流量及底质含沙量八个指标的贡献率p>0.6。将符合条件的指标进行下一步分析。（3）pearson相关性分析对通过主成分分析（pca）的候选指标进行pearson相关性分析，结果表明：鱼类生物完整性核心指标筛选中（表1），鱼类总分类单元数与总渔获量的相关系数为0.054.，两个指标间相关系数较低，均被保留下来；物理完整性核心指标筛选中（表2），由于电导率与盐度、总溶解固体的相关系数较高，所以，仅保留电导率一个指标；化学完整性核心指标筛选中（表3），由于氨氮与总氮的相关系数较高，信息重叠程度高，故保留氨氮，剔除总氮，余下氨氮、硝酸态氮、总溶解碳三个指标作为核心指标；生境完整性核心指标筛选中（表4），由于缓冲区流域面积与森林面积、建设用地面积、水田面积相关系数高，旱田面积与建设用地面积相关系数高，流量与水田面积、底质含沙量与水深相关性高，因此，仅保留森林面积、流量及底质含沙量三个核心指标。表1鱼类生物完整性候选指标相关性分析m1m2m11m20.0541*.在0.05显著性相关;**.在0.01显著性相关注：m1：鱼类总分类单元数；m2：总渔获量；表2物理完整性候选指标相关性分析p1p2p3p11p20.968**1p30.995**0.970**1*.在0.05显著性相关;**.在0.01显著性相关注：p1：电导率；p2：盐度；p3：总溶解固体；表3化学完整性候选指标相关性分析c1c2c3c4c11c20.488**1c30.822**0.897**1c40.422**-0.1040.1431*.在0.05显著性相关;**.在0.01显著性相关注：c1：氨氮；c2：硝酸态氮；c3：总氮；c4：总溶解碳；表4生境完整性候选指标相关性分析s1s2s3s4s5s6s7s8s11s20.311*1s30.813**0.0581s40.768**-0.0520.545**1s50.807*-0.0210.769**0.548**1s60.351**-0.374**0.297*0.493**0.451**1s70.583**-0.0960.445**0.702**0.645**0.605**1s80.400**-0.373**0.371**0.438**0.547**0.747**0.348**1*.在0.05显著性相关;**.在0.01显著性相关注：s1：缓冲区流域面积；s2：森林面积；s3：建设用地面积；s4：水田面积；s5：旱田面积；s6：水深；s7：流量；s8：底质含沙量；（4）标准化方法本发明对核心指标标准化统一量纲，其分为两种：对于随干扰增强而数值减小的指标，以其95th分位数为期望值记分；对于随干扰增强而数值增加的指标，以其5th分位数为期望值记分。所述计算指标分值的计算公式分别为：fij=xij/xi95×100fij=（xmax-xij）/（xmax-xi5）×100式中，fij为原始观测值；xij为j个采样点的第i个指标的标准化指数；xi95为全部采样点中第i个指标95th分位数的原始观测值；xmax为全部采样点中的最大值；xi5为全部采样点中第i个指标5th分位数的原始观测值；s2.四大体系权重比；步骤s2所述四大体系权重比是运用层次分析法对四部分指标体系进行权重，赋予鱼类生物完整性、物理完整性、化学完整性、生境完整性分别为40%、20%、20%、20%的权重，并进行加和平均得出生态系统完整性的河流水生态健康评价的得分。s3.评价标准等级化；步骤s3所述确定评价标准是以大于90th分位数作为健康的标准；将小于90th分位数的部分进行三等分，即60th~90th为较好的标准；30th~60th为一般的标准；小于30th为较差的标准，最终得出水生态健康评价的结果，如表5所示。表5生态系统完整性水生态健康评价标准健康较好一般较差生态完整性>76.7065.52~76.7059.68~65.52<59.68结果显示（图6）：浑太河流域健康点位4个，占总点位的6.15%；较好点位26个，占总点位的40%；一般点位14个，占总点位的21.55%；较差点位21个，占总点位的32.30%。从整体来看，浑太河流域一般以下点位超过一半，整体状况一般偏下。在下达河、二道河、细河、小汤河西支等点位呈现较好的水生态状况，在海城河、杨柳河、蒲河等点位呈现较差的状态。尤其是在海城河及杨柳河周边水系（t31-t40），属城镇区，人类干扰强度较大，水生态状况较差。当前第1页12