一种基于底栖动物多样性的河床稳定性阈值确定方法

本发明属于河流生态修复与水文监测，具体涉及一种基于底栖动物多样性的河床稳定性阈值确定方法。

背景技术：

1、河流是陆地生态系统中最活跃的物质和能量输送通道，其物理结构的完整性是维持水生生物多样性的基础。然而，在气候变化与强人类活动的双重影响下，全球范围内河流生境退化问题日益突出，河床稳定性下降成为威胁河流生态健康的重要因素之一。

2、我国西南山区受青藏高原抬升作用影响，河流坡降大、河床下切剧烈，崩塌、泥石流等地质灾害频发，河床物质组成松散、抗冲性弱，导致区域内多数河流河床稳定性极差。快速的河床冲淤变化不仅严重威胁沿岸居民的生命财产安全，同时破坏底栖动物的栖息环境，进而引发河流生态系统结构失衡、生物多样性下降。底栖动物作为河流生态系统的核心组成部分，其群落结构与多样性能够直接反映河流生态健康状况，且与河床稳定性存在密切的耦合关系。研究表明，稳定的河床可为底栖动物提供相对稳定的生存空间和充足的食物来源，而频繁的冲刷淤积则会直接导致底栖动物死亡、栖息地破坏及食物短缺，从而制约其种群繁衍与维持。

3、目前，河床稳定性评估主要存在以下三方面技术瓶颈。其一，评估指标与生态响应脱节。现有评估方法多依赖于水文动力参数（如流速、坡降）或底质参数（如中值粒径），这些物理指标虽然能够反映河床抵抗水流作用的力学特性，但未充分考虑生态系统的响应特征，难以准确评估河床稳定性变化对河流生态系统产生的实际生态效应。其二，物理监测与生态评价体系相互割裂。河床冲淤动态监测属于河流动力学范畴，底栖动物群落调查属于生态学范畴，两套体系长期独立运行，缺乏定量耦合的技术手段，导致工程治理中“重物理稳定、轻生态恢复”的倾向普遍存在。其三，生态阈值判定缺乏标准化方法。现有河流治理工程因缺乏科学的稳定性阈值指导，普遍存在治理效果不佳、生态修复效率低等问题。

4、为此，亟需建立一种基于底栖动物多样性的河床稳定性阈值确定方法，将生态响应与河床稳定性进行量化耦合，从而为西南山区河流治理提供科学依据，具有重要的理论与实践意义。

技术实现思路

1、本发明的目的在于提供一种基于底栖动物多样性的河床稳定性阈值确定方法，旨在通过构建河床稳定性指数与底栖动物多样性的关联模型，精准界定维持河流生态健康的河床稳定性阈值，从而为河流的稳定性评估与生态修复工程提供科学、可量化的决策依据。

2、为实现上述目的，本发明采用的技术方案是：

3、一种基于底栖动物多样性的河床稳定性阈值确定方法，包括以下步骤：

4、s1、在目标流域内筛选不少于四条分别覆盖稳定、较稳定、较不稳定、不稳定四种河床稳定类型的河流断面，设定不少于一个完整水文周期的监测周期；

5、s2、采用全站仪，在所述监测周期的起始时刻和结束时刻，分别测量各河流断面的水平距离与高程数据，同步采集河床宽度、水深、流量基础水文参数，获得断面形态变化信息；

6、s3、根据所述断面形态变化信息，分别确定各河流断面在所述监测周期内的河床下切面积和淤积面积；

7、s4、基于所述下切面积、淤积面积、河床宽度和监测周期时长计算各河流断面的河床稳定性指数；

8、s5、在各河流断面采集底栖动物样本，鉴定并统计底栖动物的物种丰度和香农-维纳多样性指数；

9、s6、分别构建物种丰度、香农-维纳多样性指数与河床稳定性指数的非线性回归关联模型，通过拐点检测确定所述物种丰度和所述香农-维纳多样性指数发生突变对应的河床稳定性指数临界值，该临界值即为河床稳定性阈值。

10、进一步地，步骤s1中所述监测周期涵盖丰水期、平水期和枯水期，且监测周期时长为1年。

11、进一步地，步骤s1中筛选河流断面时，兼顾流域地质背景、植被覆盖度及人类活动干扰程度的一致性，排除非河床稳定性因素对底栖动物多样性的干扰。

12、进一步地，步骤s2中，全站仪测量时，测点沿断面左岸至右岸均匀布设，常规断面测量点间距不大于2m，各测点测量误差不超过±0.01m；河床地形复杂断面的测点间距调整为1~1.5m。

13、进一步地，步骤s3中根据两个时刻的断面水平距离与高程数据，绘制各河流断面的初始剖面图和末期剖面图，分别核算各河流断面初始剖面图与末期剖面图的面积，通过对比确定监测周期内的下切面积和淤积面积。

14、进一步地，步骤s4中所述河床稳定性指数的计算公式为：

15、

16、其中，为河床宽度，为监测周期时长，为监测周期内河床下切面积，为监测周期内河床淤积面积。

17、进一步地，所述河床稳定性指数的量纲为年/米，表征单位河宽条件下，河床横断面发生单位面积冲淤变化所需的平均时间；值越大，表示河床抵抗冲淤变化的能力越强，底栖动物栖息地的物理稳定性越高。

18、进一步地，步骤s5中在丰水期结束后1~2周内采集底栖动物样本，在每个河流断面均匀布设3~5个采样子点，采用彼得森采泥器进行采集，每个子点采集面积为0.0625，采集深度为河床表层0~20cm；所述香农-维纳多样性指数的计算公式为：

19、

20、其中，为第个物种的个体数量，为该采样点底栖动物总个体数量，为物种丰度。

21、进一步地，步骤s6中拐点检测具体流程：对散点数据进行loess平滑处理消除局部波动；对平滑后的曲线计算数值二阶导数，设定二阶导数变号阈值筛选拐点候选点；以拐点候选点为分界将数据分为两段分别进行线性回归，当两段回归的决定系数均大于0.85且两段斜率经t检验差异显著时，确认该点为有效拐点。

22、进一步地，分别获取物种丰度与河床稳定性指数关联模型的拐点对应的值、香农-维纳多样性指数与河床稳定性指数关联模型的拐点对应的值，取两个值的算术平均值作为河床稳定性阈值；若两个拐点值的相对差异超过20%，则以香农-维纳多样性指数与河床稳定性指数关联模型对应的拐点值作为最终河床稳定性阈值。

23、上述技术方案的有益效果为：

24、第一，本发明提出了一种新的河床稳定性评价指标指数。该指数以单位河宽河床发生单位面积冲淤变化所需的平均时间作为稳定性度量，具有明确的物理含义。相较于传统依赖流速、坡降或中值粒径等单一参数的评估方法，指数综合反映了监测周期内河床的冲刷与淤积总强度，能够更全面地刻画河床的动态变化特征。将该指数与底栖动物多样性建立关联，实现了从生态视角评价河床稳定性的技术突破，弥补了现有评估方法缺乏生态响应的不足。

25、第二，本发明确定了河床稳定性影响底栖动物多样性的临界阈值。通过构建指数与物种丰度s、香农-维纳多样性指数h'的非线性回归关联模型，结合二阶导数拐点检测和分段线性回归显著性检验，能够精准识别底栖动物多样性发生突变的临界值。实施例验证结果表明，s-模型和-模型的决定系数分别达到0.92和0.90，拟合效果好，所确定的阈值能够有效区分底栖动物群落处于适宜状态还是退化状态。

26、第三，本发明为西南山区河流治理工程提供了明确的量化目标。治理工程可将指数提升至阈值以上作为设计依据和验收标准，解决了以往工程中缺乏科学指导的问题。应用该阈值指导不稳定河段构建阶梯深潭系统后，后期跟踪监测显示底栖动物多样性显著提升，验证了本方法确定的河床稳定性阈值具有实际工程指导价值。

27、第四，本发明方法规范性强，可操作性好。整个操作流程涵盖断面筛选与监测、全站仪高精度测量、断面剖面图绘制、稳定性指数计算、拐点检测阈值确定等环节，各步骤均采用成熟的测量和绘图工具，操作步骤清晰明确，便于行业内推广应用。该方法有助于提升治理工程的科学性和针对性，改善河流生态环境，保障区域生命财产安全与生态安全，具有重要的工程应用价值。