一种内河航道整治工程生态影响模拟预测方法与流程

本发明涉及一种内河航道整治领域，特别涉及一种内河航道整治工程生态影响模拟预测方法。

背景技术：

目前，河流水文情势和动力条件等水文机制的改变是河流生态系统退化的主要驱动力。内河航道整治工程(疏浚、清礁、筑坝、护滩、护岸等)改变了航道水位、流态、流速在三维空间上的显著变化及河床变形与冲淤分布特征，进而导致河流水生生物栖息地的时空分布发生不同程度的变化。目前的水生态模拟理论构筑起了水力学计算与生态学特征表述间的理论联系，能够依据特有鱼类的栖息、产卵习性通过对水文过程的调节，人工提供出适宜的栖息环境。水生生物，特别是指示鱼类及其关键生境，是航道整治工程生态影响预测和后评估的关注焦点，当前基于河道内流量增加方法开展了内河航道整治工程对物理栖息地影响的模拟和评估研究，仍多采用一、二维模型，因物理栖息地三维空间特性对生物多样性有重要影响，上述模拟方法仍有一定的局限性，因此合理科学预测内河航道整治工程对物理栖息地的影响是十分必要的。

技术实现要素：

针对内河航道整治工程建设对河流栖息地条件的影响预测难题，本发明基于不同通航条件下的物理栖息地指示性物种的特征与变化趋势，建立三维水流泥沙数学模型水力参量与栖息地特征量间的响应过程，为科学评价内河航道工程方案生态影响提供一种更加科学有效的模拟预测方法。本发明为解决公知技术中存在的技术问题而提供一种三维的内河航道整治工程生态影响模拟预测方法。

本发明为解决公知技术中存在的技术问题所采取的技术方案是：一种内河航道整治工程生态影响模拟预测方法，调查统计整治河段指示水生生物在流速、水深、底质和含沙量的某个区间范围的出现频次，采用多元线性回归建立整治河段指示性水生生物响应与栖息地生境因子变量之间的关系，建立该水生生物栖息地的适宜度模型；建立整治河段的三维水沙数学模型，由三维水沙数学模型，预测整治河段在不同航道整治措施参数下水深、流速、流态及含沙量的三维时空分布数据；将预测得到的水深、流速、流态及含沙量三维时空分布数据，输入适宜度模型中，预测整治河段在不同航道整治措施参数下的该水生生物栖息地的适宜度的变化，通过指示生物栖息地生境因子适宜度的分布范围和适宜生境面积的大小，优化航道整治措施参数和平面布局。

进一步地，根据整治河段水生生物状况的历史资料，以及通过渔获物和鱼探仪相结合的生态调查与监测方法，确定整治河段的指示水生生物以及构成该水生生物栖息地生境因子及生境因子的时空分布。

进一步地，指示水生生物为指示鱼类。

进一步地，该水生生物的生境因子包括栖息地流速、流态、水深、含沙量。

进一步地，由三维水沙控制方程，建立整治河段的三维水沙数学模型。

进一步地，三维水沙控制方程包括流体连续性方程、流体动量方程、标准k-ε紊流方程以及泥沙输移控制方程。

进一步地，三维水沙数学模型，采用有限体积法进行离散求解,三维计算域在平面上采用三角形网格离散、沿水深方向进行分层离散。

本发明具有的优点和积极效果是：本发明着重考虑了内河航道整治工程对河流生态系统中指示性生物关键生境要素的三维分布影响，采用生态学、水力学及河流动力学、计算水力学等多种方法科学合理确定内河航道整治工程对指示性生物栖息地条件的三维分布影响，在开发河流通航功能的同时，兼顾河流生态功能的保护与修复。

本发明建立的内河航道整治工程生态影响模拟预测方法对内河航道生态化建设具有通用性，通过本发明可准确模拟预测分析不同整治参数条件下鱼类栖息环境三维分布变化，在实现河流通航功能前提下，为水生生物提供多样性的栖息条件，实现河流通航--生态功能的融合，预测内河航道整治工程措施参数条件下整治位点水文变化和河道物理结构改变对指示性物种生境的影响。内河航道整治布置方案时，统筹考虑航道整治效果和河流生态功能的保护，既注重航道尺度的满足，同时又关注对河流洪水边滩、河湾、汊道、江心洲、沙洲等浅水生境的保留与营造。

附图说明

图1为本发明步骤流程图；

图2为模型计算网格剖分示意图；

图3为本发明实施例中一种指示水生生物的流速适宜度曲线；

图4为本发明实施例中指示水生生物的水深适宜度曲线；

图5为本发明实施例中指示水生生物的含沙量适宜度曲线；

图6为本发明实施例中整治河段中指示水生生物对应流速的适宜性指数分布图；

图7为本发明实施例中整治河段中指示水生生物对应水深的适宜性指数分布图；

图8为本发明实施例中整治河段中指示水生生物综合适宜度因子分布图。

图中，uhsi表示指示水生生物对应流速的适宜性指数；hhsi表示指示水生生物对应水深的适宜性指数；csf表示指示水生生物对应综合适宜度因子指数。

具体实施方式

为能进一步了解本发明的发明内容、特点及功效，兹列举以下实施例，并配合附图详细说明如下：

请参见图1至图8，一种内河航道整治工程生态影响模拟预测方法，调查统计整治河段指示水生生物在流速、水深、底质和含沙量的某个区间范围的出现频次，采用多元线性回归建立整治河段指示性水生生物响应与栖息地生境因子变量之间的关系，建立该水生生物栖息地的适宜度模型；建立整治河段的三维水沙数学模型，由三维水沙数学模型，预测整治河段在不同航道整治措施参数下水深、流速、流态及含沙量的三维时空分布数据；将预测得到的水深、流速、流态及含沙量三维时空分布数据，输入适宜度模型中，预测整治河段在不同航道整治措施参数下的该水生生物栖息地的适宜度的变化，通过指示生物栖息地生境因子适宜度的分布范围和适宜生境面积的大小，优化航道整治措施参数和平面布局。

上述指示水生生物可为指示鱼类。上述水生生物的生境因子可包括栖息地流速、流态、水深、含沙量等。

可根据整治河段水生生物状况的历史资料，以及通过渔获物和鱼探仪相结合的生态调查与监测方法，确定整治河段的指示水生生物以及构成该水生生物栖息地生境因子及生境因子的时空分布。

优选地，采用渔获物、生物学、早期资源和鱼探仪等多种技术手段，探明整治工程建设前后受影响河段的鱼类资源变化情况、产卵场位置、规模大小，尤其是特有鱼类资源现状。同时，结合资源量、受威胁程度及其对环境敏感度，综合确定整治河段的指示性生种，一般可选择鱼类作为指示性生种。鱼类在水生生态系统处于食物链顶端，对其他种群的存在和丰度有着重大影响，且对水环境的变化敏感，如以单一鱼类物种作为指示物种的研究案例达47％，以鱼类种群作指示物种的案例达25％，以底栖动物的达20％。

可通过渔获物、鱼探仪相等监测方法确定指示鱼类出现的位置，采用整治河段的三维水沙数学模型反演指示鱼类流速流态、水深、含沙量及覆盖物生境因子的时空分布特征。

可采用多元线性回归建立整治河段指示性物种响应与栖息地变量水深、流速及含沙量等生境因子之间的关系，获得水生生物对应单个生境因子的适宜性指数。统计整治河段指示鱼类在某一流速、水深、底质和含沙量等区间范围的出现的频次作为该区间的适宜度，即指示水生生物对应生境因子的相对丰度，认为出现频次最高的区间适宜度为1，采用归一化处理即得到该指示生物对应水深、流速以及含沙量等各个生境因子的适宜度曲线，如图3～5。

划分区间采用sturge计算最佳间隔：

式中，i为最佳间隔；r为指标变化范围；n为观察到的指示性生物的数量。

进一步地，根据三维水沙模型计算出整治河段的水深、流速及含沙量分布，结合整治河段指示性物种的栖息地生境因子适宜度曲线，该指示生物栖息地的适宜度模型建立如下：

hsii＝min(vi,di,si)

式中，wua为整治河段区域指示水生生物的加权栖息地面积；ai为第i单元的表面积；hsii为第i单元的栖息地生境因子综合适宜性，综合适宜度因子(combinedsuitabilityfactor，csf)也可采用最小值法确定；vi为第i单元的流速适宜性指数；di为第i单元的水深适宜性指数；si为第i单元的含沙量适宜性指数。

进一步地，采用指示生物栖息地的适宜度模型和三维水沙数学模型，预测整治河段在不同航道整治措施参数下的指示生物栖息地的生境因子适宜度的分布和适宜生境面积的大小和分布范围，如图6和图7分别为本发明实施例中整治河段中指示生物对应水深和流速的适宜性指数分布图，图8为本发明实施例中整治河段中指示生物综合适宜度因子分布图。

优选地，可由三维水沙控制方程，建立整治河段的三维水沙数学模型。三维水沙控制方程可包括流体连续性方程、流体动量方程、标准k-ε紊流方程以及泥沙输移控制方程。

基于navier-stokes建立流体动量方程、标准k-ε紊流方程，navier-stokes，即纳维-斯托克斯方程，纳维-斯托克斯方程描述粘性不可压缩流体动量守恒的运动方程，简称n-s方程。n-s方程反映了粘性流体(又称真实流体)流动的基本力学规律。

其中：

流体连续性方程为：

流体动量方程为：

标准k-ε紊流方程为：

泥沙输移控制方程为：

上式中：η为水位；u，v，w分别为速度矢量沿三个坐标轴x，y，z的分量；g为重力加速度；ρ为水体密度；pa为大气压；q表示非静压项；ρ0为参考密度；f为科式力系数；vt为涡粘系数，取cμ为常数，取cu＝0.09，ν^h为水平方向的涡粘系数；ν^v为垂直方向的涡粘系数；gk为由于平均速度梯度引起的紊动能生成项，其中xi或xj(i,j＝1,2,3)为x、y、z三个坐标轴分量，x1为x坐标轴分量；x2为y坐标轴分量；x3为z坐标轴分量，ui或uj(i,j＝1,2,3)分别为流速矢量沿x、y、z坐标轴的分量，u1为流速矢量沿x坐标轴的分量，u2为流速矢量沿y坐标轴的分量，u3为流速矢量沿z坐标轴的分量；k为紊动动能；ε为紊动耗散率；xk(k＝1,2,3)为x、y、z三个坐标轴分量，其中x1为x坐标轴分量；x2为y坐标轴分量；x3为z坐标轴分量；常数项：cε1＝1.44，cε2＝1.92，σk＝1.0，σε＝1.3；s为含沙量；p′为床沙质的孔隙率；ws为泥沙沉降速度；kh为水平方向的泥沙扩散系数；kv为垂直方向的泥沙扩散系数；zb为河床总变形。qsx为x方向上的悬移质输沙率；qsy为y方向上的悬移质输沙率；qbx为x方向上的推移质输沙率；qby为y方向上的推移质输沙率。

三维水沙数学模型，可采用有限体积法进行离散求解,三维计算域在平面上可采用三角形网格离散、沿水深方向进行分层离散。

计算在整治河段的水深、流速以及流速的分布，将整治河段的河床地形文件网格化，基于连续性方程、navier-stokes、泥沙输移方程建立研究河段三维水沙模型，采用三角形网格离散平面计算域和水位函数法捕捉自由表面，提高模型适应复杂边界的能力和计算效率。

采用三角形离散平面计算域，沿水深方向采用垂向分层离散，三维计算域被剖分为若干棱柱形单元，如图2，由于这种剖分方式不会产生杂乱的四面体等结构，采用分开来定义变量的水平和垂向的空间布置，水平速度项定义在棱柱形网格的垂向面上，垂向速度定在棱柱形网格的上下面上，水位、非静压项、紊动动能、紊动耗散率等标量定义在棱柱形网格的中心(详见：吕彪.基于非结构化网格的具有自由表面水波流动数值模拟研究[d].大连：大连理工大学，2010)。

采用水位函数法捕捉自由表面，即沿水深方向积分连续方程，并利用自由表面的运动学边界条件和不可穿透底面边界条件，由莱布尼茨公式可得到的沿水深积分的连续方程，称为水位演化方程，求解该方程即可得到自由表面的位置。沿水深方向积分连续性方程有：

上式由leibniz和微积分定理可得：

不可穿透底面边界条件：

自由表面的运动学边界条件：

可得水位演化方程

上式中：η为水位；u，v，w分别为速度矢量沿三个坐标轴x，y，的分量；h为水深；ub为底层流速矢量沿x坐标轴的分量；us为表层流速矢量沿x坐标轴的分量；vb为底层流速矢量沿y坐标轴的分量；vs为表层流速矢量沿y坐标轴的分量；wb为底层流速矢量沿z坐标轴的分量；ws为表层流速矢量沿z坐标轴的分量。

由上述方程求解得到的整治河段的水深、流速、流态及含沙量(详见：吕彪.基于非结构化网格的具有自由表面水波流动数值模拟研究[d].大连：大连理工大学，2010)。

由上述三维水沙数学模型，可以预测整治河段在不同航道整治措施参数下的水深、流速、流态以及含沙量的三维时空分布数据；将预测得到的水深、流速、流态及含沙量三维时空分布数据，与该指示生物栖息地的适宜度模型结合，从而预测整治河段在不同航道整治措施参数下的该指示生物栖息地的适宜度的三维空间变化趋势，结合指示生物栖息地的适宜度的分布范围、适宜生境面积的大小和航道通航目标，优化航道整治措施参数和平面布局，得到最终通航-生态融合的航道整治方案，以达到航道整治目标和水生境保护的有机融合。

以上所述的实施例仅用于说明本发明的技术思想及特点，其目的在于使本领域内的技术人员能够理解本发明的内容并据以实施，不能仅以本实施例来限定本发明的专利范围，即凡本发明所揭示的精神所作的同等变化或修饰，仍落在本发明的专利范围内。