一种耦合生境适宜模型的高效高精度栖息地模拟方法

1.本发明属于生物栖息地模拟方法技术领域，具体涉及一种耦合生境适宜模型的高效高精度栖息地模拟方法。


背景技术：

2.一直以来，人们采取一系列水库生态运行和河流生态恢复的措施来缓解水利工程对河流生态系统的压力，为此提出许多模拟方法来解决此类问题，其中生境模型成为世界上应用最广泛的方法，如phabsim模型和river2d 模型。而生境因子的准确性对于生境模型更准确地模拟河流生态环境和评价水生生境的质量至关重要。近些年，为提高生境变量的时空分辨率，引入越来越多的水力模型与生境适宜性模型相结合，如delft3d、mike21、ssiim 模型等二维或三维流体动力学模型。但对于高分辨率的地形数据，以及大型且具有复杂流态的河流而言，这些模型仍然存在模拟精度不足，模拟时间长等缺点。因此，亟需发展一种具有能够模拟高精度二维水动力过程的生境模型，来评价单种或多种水生生物的生境质量。基于此，构建一种耦合生境适宜模型的高效高精度栖息地模拟方法，以准确模拟河流水力特性和生物栖息地，使其成为评价生物生境质量的有力工具。


技术实现要素：

3.本发明的目的是提供一种耦合生境适宜模型的高效高精度栖息地模拟方法，该方法可提高生境变量的时空分辨率；对于高分辨率地形数据及大型且具有复杂流态的河流而言，还可弥补当前很多模型存在的模拟精度不足，模拟时间长等缺点，能够更准确地模拟河流生态环境，评价单种或多种水生生物的生境质量。
4.本发明所采用的技术方案是，一种耦合生境适宜模型的高效高精度栖息地模拟方法，具体按照以下步骤实施：
5.步骤1、利用arcgis软件生成高分辨率dem地形数据表征真实地形，以ascii格式输出；
6.步骤2、调研收集研究河段目标物种以及目标物种的生境因子，利用偏好适宜曲线法确定出目标物种的生境因子适宜曲线，包括水深、流速和底质；
7.步骤3、将高分辨率dem地形数据离散成矩形网格，根据已知边界条件对相关参数进行赋值；
8.步骤4、将得出的dem地形数据及各生境因子适宜条件在gpu中分配内存，并复制到gpu变量中；
9.步骤5、基于godunov格式的中心有限体积法空间离散二维浅水方程，求解通量项、底坡源项和摩阻源项；
10.步骤6、由muscl重组法和龙格库塔法进一步实现时间和空间的二阶精度，时间和空间步长计算稳定性由cfl条件限制；
11.步骤7、基于ifim理论，由得到的目标物种高精度生境因子及目标物种的生境适宜
条件，确定各生境因子的适宜性指数hsi计算出综合适宜性指数 csi；
12.步骤8、通过cudamemcpy函数将gpu计算得到的各生境因子数据及综合适宜指数csi复制到cpu中，计算目标物种有效栖息地加权可利用面积wua；
13.步骤9、根据计算得到的水深、流速和底质数据，以及水深、流速和底质三种生境因子的综合适宜指数csi数据，利用matlab软件输出水深、流速、底质和综合适宜性指数csi分布图；至此，耦合生境适宜模型的高效高精度栖息地模拟方法构建完成。
14.本发明的特点还在于，
15.步骤5中二维浅水方程表示如下：
[0016][0017][0018][0019]
式中：t为时间步长，单位为s；x、y分别为水平坐标和纵向坐标； q为流量变量矢量；h表示水深，单位为m；u、v分别为x、y方向上的平均流速，单位为m/s；f、g分别为x、y方向上的通量矢量；s为源项矢量；其中，sb为底坡源项；sf为摩阻源项；zb为河床高程，m；cf为床面糙率系数，cf＝gn2/h
1/3
，n为曼宁系数，s/m
1/3
。g为重力加速度，取 9.81，单位为m/s2。
[0020]
步骤5具体如下：
[0021]
步骤5.1、第i个控制单元网格内，采用中心有限体积法对swes控制方程组积分，基于高斯散度定理通量项面积线积分形式为：
[0022][0023]
通量f(q)
·
n线积分为：
[0024][0025]
式中：ω为第i个控制单元网格的体积，单位为m3；г为第i个控制单元网格的边界；n为边界г所对应的外法线方向的单位向量；sb表示底坡源项通量；sf表示摩阻源项通量；k为控制单元网格边的编号；lk为第i个控制单元网格第k个边的边长；
[0026]
步骤5.2、在godunov方法中实现hllc黎曼求解器步骤为：
[0027]
左波速s
l
、中间波速s
*
及右波速sr计算公式为，
[0028]
[0029][0030][0031]
其中，中间变量h
*
和根据式计算，
[0032][0033][0034]
从而计算通量f(q)
·
n的hllc黎曼求解器计算公式为，
[0035][0036]
其中，控制体单元网格界面中间通量f
*
为，
[0037][0038]
控制体单元网格界面通量中间波两侧，f
*l
和f
*r
为，
[0039][0040]
式中，f
l
和fr分别为单元网格界面左、右侧的通量；f
*
、f
*l
和f
*r
分别为控制体单元网格界面通量中间波两侧的左波通量、中间波通量和右波通量；u
⊥
为沿交界面的切向速度，u
⊥
＝-uny+vn
x
，u
⊥
l
和u
⊥r分别为控制体单元网格交界面的左侧、右侧切向速度，单位为m3/s；n
x
和ny分别为沿x 和y的指示方向；u和v分别为沿x和y的单宽流量，单位为m/s；变量 q
⊥
＝[h,q
xnx
+qyny]；q
x
和qy分别为沿x和y的单宽流量，单位为m2/s；s
l
、 sm、sr分别为左波速、中间波速、右波速，单位为m/s；u
⊥
为通过网格单元界面垂直的速度，单位为m/s；u
⊥
＝un
x
+vny；h
*
和u
*
⊥
为中间变量；g 为重力加速度，取9.81，单位为m/s2；
[0041]
步骤5.3、底坡源项采用底坡通量法求解，将单元内底坡源项转化为单元边界上的通量。
[0042]
采用底坡通量法得到某一边f处的底坡通量的矢量形式为，
[0043][0044]
式中，f
sf
(q)
·
为某一边f处的底坡通量；nf为面f单位外法向量；n
fx
和n
fy
是nf与(面
f单位外法向量)的组成部分，h
l
和z
bl
分别为网格单元中心处的水深和底部高程，与z
ml
网格边上的水深和底部高程，g为重力加速度，单位为m/s2；
[0045]
步骤5.4、摩阻源项采用显隐式摩阻处理法求解，将隐式转化为显式，采用显隐式摩阻处理法求解摩阻源项计算公式为，
[0046][0047]
其中，单宽流量更新为时间步长n+1时刻计算公式为，
[0048][0049][0050]
式中，sf为摩阻源项通量，s
fx
和s
fy
分别为x和y方向的摩阻源项通量；cf为床面糙率系数，cf＝gn2/h
1/3
，n为曼宁系数，s/m
1/3
；h为水深，单位为m；q
x
、qy分别为x、y方向上平均流速，单位为m2/s；q
xn+1
、q
yn+1
分别为时间步长n+1时刻在x、y方向上的单宽流量，单位为m2/s；q
xn
、q
yn
分别为时间步长n时刻在x、y方向上的单宽流量，单位为m2/s；δt为时间步长，单位为s；δq
x
为时间步长n时刻和n+1时刻的x方向上单宽流量之差，即δq
x
＝q
xn+1-q
xn
，单位为m2/s；为了避免分母为零，θ取1.0e-12
；
[0051]
步骤5.5、引入静水重构法求解干湿边界问题，即将界面处重构的负水深值强制为零来维持非负水深，并修改相应地水位、底坡高程和界面两侧的流量，并将重构的新值代入近似黎曼求解器中来计算质量和动量的通量，保证干湿界面处的平衡状态和质量守恒。
[0052]
步骤6具体如下：
[0053]
步骤6.1、求解出通量项、底坡源项和摩阻源项之后，采用muscl重构法在第i个控制单元网格通过同一限制变量梯度插值求解出第i个控制单元网格界面上的变量值；
[0054]
步骤6.2、由两步runge-kutta方法获得二阶时间精度，在新时间步长中第i个控制单元网格的值为：
[0055][0056]
中间变量值为
[0057][0058]
[0059]
式中：i为控制单元网格数；ω为第i个控制单元网格的体积，单位为 m3；q
n+1
为时间步长n+1时刻的单宽流量，单位为m2/s，qn为时间步长n 时刻的单宽流量，单位为m2/s；δt为时间步长，单位为s；s(qn)为时间步长n时刻的源项；
[0060]
步骤6.3、在空间步长δx下，时间步长δt的取值的计算公式如下：
[0061][0062]
式中：δx
min
表示网格单元中心到对应界面的最小距离，单位为m；u为流速，单位为m/s；g为重力加速度，取值9.81m/s2；h为水深，单位为m。
[0063]
步骤7具体如下：
[0064]
步骤7.1、根据目标物种高精度生境因子及目标物种适宜条件，采用偏好适宜曲线法，确定出目标物种各生境因子适宜指数hsi，具体为，根据第 i个控制单元网格内水深veli、流速depi和底质subi三种生境因子的模拟数据，采用偏好适宜曲线法，计算出目标物种水深、流速和底质的适宜指数种水深、流速和底质的适宜指数和规定生境适宜指数范围为0～1，最适宜生境适宜指数为1，最不适宜生境适宜指数为0；
[0065]
步骤7.2、第i个控制单元网格内综合适宜指数csi由下式得出：
[0066][0067][0068][0069]
式中：i为控制单元网格数；csi为第i个控制单元网格内的综合适宜指数；hsi为水深、流速和底质等每个生境因子在第i个控制单元网格内的生境适宜指数；depi,veli,subi分别为第i个控制单元网格内水深、流速、底质或覆盖物，分别为第i个控制单元网格内水深、流速、底质或覆盖物的适宜指数；式(4)将所选择的生境因子适宜指数相乘，表示影响因子综合作用结果；式(5)考虑当某一生境因子较为不利时，组成栖息地生境因子之间的补偿影响；式(6)将最不适于鱼种生存的生境因子适宜指数作为组合适宜指数。
[0070]
步骤8具体如下：
[0071]
加权可利用面积wua计算式；
[0072][0073]
式中：wua为栖息地加权可利用面积，单位为m2；i为控制单元网格数；ai为第i个控制单元网格的面积，单位为m2。
[0074]
本发明的有益效果是，一种耦合生境适宜模型的高效高精度栖息地模拟方法，是针对现有数值模型在处理高分辨率的地形数据，以及大型且具有复杂流态的河流方面，存在模拟精度不足，模拟时间长等问题，提出一种具有高精度高效的生境适宜模型，以提高生境模型中生境因子精确性，为准确模拟河流水力特性和生物栖息地提供有效的技术支撑，使其成为评价生物生境质量的有力工具。
附图说明
[0075]
图1是研究河段的高分辨率地形数据；
[0076]
图2是花斑裸鲤产卵期生境因子适宜性曲线，其中，图2(a)是水深适宜性曲线图，图2(b)是流速适宜性曲线图，图2(c)是底质适宜性曲线图；
[0077]
图3是基于godunov格式的二维浅水方程中心有限体积法求解过程；
[0078]
图4是基于gpu加速的生境适宜模型与高精度水动力过程耦合过程；
[0079]
图5是gast和river2d模型下不同流量与wua关系模拟结果比较图；
[0080]
图6是流量74m3/s条件下水深分布图；
[0081]
图7是流量74m3/s条件下流速分布图；
[0082]
图8是流量74m3/s条件下综合适宜指数csi分布图。
具体实施方式
[0083]
下面结合附图和具体实施方式对本发明进行详细说明。
[0084]
本发明一种耦合生境适宜模型的高效高精度栖息地模拟方法，具体按照以下步骤实施：
[0085]
步骤1、利用arcgis软件生成高分辨率dem地形数据表征真实地形，以ascii格式输出；
[0086]
步骤2、调研收集研究河段目标物种以及目标物种的生境因子，利用偏好适宜曲线法确定出目标物种的生境因子适宜曲线，包括水深、流速和底质；
[0087]
步骤3、将高分辨率dem地形数据离散成矩形网格，根据已知边界条件对相关参数进行赋值；
[0088]
步骤4、引入gpu加速技术以提高模拟效率，将得出的dem地形数据及各生境因子适宜条件在gpu(graphics processing unit,即图形处理器,又称为显示核心)中分配内存，并复制到gpu变量中；
[0089]
步骤5、基于godunov格式的中心有限体积法空间离散二维浅水方程，求解通量项、底坡源项和摩阻源项，获取更高空间离散精度，以提高激波捕捉的分辨率；
[0090]
步骤5中二维浅水方程表示如下：
[0091][0092][0093][0094]
式中：t为时间步长，单位为s；x、y分别为水平坐标和纵向坐标； q为流量变量矢量；h表示水深，单位为m；u、v分别为x、y方向上的平均流速，单位为m/s；f、g分别为x、y方向
上的通量矢量；s为源项矢量；其中，sb为底坡源项；sf为摩阻源项；zb为河床高程，m；cf为床面糙率系数，cf＝gn2/h
1/3
，n为曼宁系数，s/m
1/3
。g为重力加速度，取 9.81，单位为m/s2。
[0095]
步骤5具体如下：
[0096]
步骤5.1、第i个控制单元网格内，采用中心有限体积法对swes控制方程组积分，基于高斯散度定理通量项面积线积分形式为：
[0097][0098]
通量f(q)
·
n线积分为：
[0099][0100]
式中：ω为d第i个控制单元网格的体积，单位为m3；г为第i个控制单元网格的边界；n为边界г所对应的外法线方向的单位向量；sb表示底坡源项通量；sf表示摩阻源项通量；k为控制单元网格边的编号；lk为第i个控制单元网格第k个边的边长；
[0101]
步骤5.2、通量项计算中，godunov格式考虑迎风格式(upwind scheme)的信息源特性；对于求解中出现一个任意间断的阶梯函数(黎曼问题)，通过hllc黎曼求解器，可直接得到控制体间数值通量的近似值，并可自动满足 godunov格式。在godunov方法中实现hllc黎曼求解器步骤为：
[0102]
考虑到存在网格单元无水的情况，左波速s
l
、中间波速s
*
及右波速sr计算公式为，
[0103][0104][0105][0106]
其中，中间变量h
*
和根据式计算，
[0107][0108][0109]
从而计算通量f(q)
·
n的hllc黎曼求解器计算公式为，
[0110][0111]
其中，控制体单元网格界面中间通量f
*
为，
[0112][0113]
控制体单元网格界面通量中间波两侧，f
*l
和f
*r
为，
[0114][0115]
式中，f
l
和fr分别为单元网格界面左、右侧的通量；f
*
、f
*l
和f
*r
分别为控制体单元网格界面通量中间波两侧的左波通量、中间波通量和右波通量；u
||
为沿交界面的切向速度，u
⊥
＝-uny+vn
x
，u
⊥
l
和u
⊥r分别为控制体单元网格交界面的左侧、右侧切向速度，单位为m3/s；n
x
和ny分别为沿x和y的指示方向；u和v分别为沿x和y的单宽流量，单位为m/s；变量q
⊥
＝[h,q
xnx
+qyny]；q
x
和qy分别为沿x和y的单宽流量，单位为m2/s； s
l
、sm、sr分别为左波速、中间波速、右波速，单位为m/s；u
⊥
为通过网格单元界面垂直的速度，单位为m/s；u
⊥
＝un
x
+vny；h
*
和为中间变量；g为重力加速度，取9.81，单位为m/s2；
[0116]
步骤5.3、底坡源项采用底坡通量法求解，将单元内底坡源项转化为单元边界上的通量，稳定性更佳。
[0117]
采用底坡通量法得到某一边f处的底坡通量的矢量形式为，
[0118][0119]
式中，f
sf
(q)
·
为某一边f处的底坡通量；nf为面f单位外法向量；n
fx
和n
fy
是nf与(面f单位外法向量)的组成部分，h
l
和z
bl
分别为网格单元中心处的水深和底部高程，与z
ml
网格边上的水深和底部高程，g为重力加速度，单位为m/s2；
[0120]
步骤5.4、摩阻源项采用显隐式摩阻处理法求解，将隐式转化为显式，消除常规隐式的冗余迭代，以兼顾计算精度和效率。采用显隐式摩阻处理法求解摩阻源项计算公式为，
[0121][0122]
其中，单宽流量更新为时间步长n+1时刻计算公式为，
[0123]
[0124][0125]
式中，sf为摩阻源项通量，s
fx
和s
fy
分别为x和y方向的摩阻源项通量；cf为床面糙率系数，cf＝gn2/h
1/3
，n为曼宁系数，s/m
1/3
；h为水深，单位为m；q
x
、qy分别为x、y方向上平均流速，单位为m2/s；q
xn+1
、q
yn+1
分别为时间步长n+1时刻在x、y方向上的单宽流量，单位为m2/s；q
xn
、q
yn
分别为时间步长n时刻在x、y方向上的单宽流量，单位为m2/s；δt为时间步长，单位为s；δq
x
为时间步长n时刻和n+1时刻的x方向上单宽流量之差，即δq
x
＝q
xn+1-q
xn
，单位为m2/s；为了避免分母为零，θ取1.0e-12
；
[0126]
步骤5.5、引入静水重构法求解干湿边界问题，即将界面处重构的负水深值强制为零来维持非负水深，并修改相应地水位、底坡高程和界面两侧的流量，并将重构的新值代入近似黎曼求解器中来计算质量和动量的通量，保证干湿界面处的平衡状态和质量守恒。以保持干湿界面处的平衡状态和质量守恒，达到满足全稳条件。
[0127]
步骤6、由muscl重组法和龙格库塔法(runge-kutta)进一步实现时间和空间的二阶精度，时间和空间步长计算稳定性由cfl(courant friedrichs lewy)条件限制；
[0128]
步骤6具体如下：
[0129]
步骤6.1、求解出通量项、底坡源项和摩阻源项之后，采用muscl重构法在第i个控制单元网格通过同一限制变量梯度插值求解出第i个控制单元网格界面上的变量值，实现空间二阶精度，以平衡数值扩散和振荡问题；
[0130]
步骤6.2、由两步runge-kutta方法获得二阶时间精度，在新时间步长中第i个控制单元网格的值为：
[0131][0132]
中间变量值为
[0133][0134][0135]
式中：i为控制单元网格数；ω为第i个控制单元网格的体积，单位为m3；q
n+1
为时间步长n+1时刻的单宽流量，单位为m2/s，qn为时间步长n 时刻的单宽流量，单位为m2/s；δt为时间步长，单位为s；s(qn)为时间步长n时刻的源项；
[0136]
步骤6.3、在空间步长δx下，时间步长δt的取值的计算公式如下：
[0137][0138]
式中：δx
min
表示网格单元中心到对应界面的最小距离，单位为m；u为流速，单位为m/s；g为重力加速度，取值9.81m/s2；h为水深，单位为m。
[0139]
步骤7具体如下：
[0140]
步骤7.1、根据目标物种高精度生境因子及目标物种适宜条件，采用偏好适宜曲线法，确定出目标物种各生境因子适宜指数hsi，具体为，根据第 i个控制单元网格内水深depi、流速veli和底质subi三种生境因子的模拟数据，采用偏好适宜曲线法，计算出目标物种水深、流速和底质的适宜指数种水深、流速和底质的适宜指数和规定生境适宜指数范围为0～1，最适宜生境适宜指数为1，最不适宜生境适宜指数为0；
[0141]
步骤7.2、结合实际情况，选择合适的计算方法，计算出综合适宜指数 csi(combined suitability index)；第i个控制单元网格内综合适宜指数csi 由下式得出：
[0142][0143][0144][0145]
式中：i为控制单元网格数；csi为第i个控制单元网格内的综合适宜指数；hsi为水深、流速和底质等每个生境因子在第i个控制单元网格内的生境适宜指数；depi,veli,subi分别为第i个控制单元网格内水深、流速、底质或覆盖物，分别为第i个控制单元网格内水深、流速、底质或覆盖物的适宜指数；式(4)将所选择的生境因子适宜指数相乘，表示影响因子综合作用结果；式(5)考虑当某一生境因子较为不利时，组成栖息地生境因子之间的补偿影响；式(6)将最不适于鱼种生存的生境因子适宜指数作为组合适宜指数。
[0146]
步骤7、基于ifim理论，由得到的目标物种高精度生境因子及目标物种的生境适宜条件，确定各生境因子的适宜性指数hsi(habitat suitability index)，结合实际情况，选择合适的计算方法，计算出综合适宜性指数csi (combined suitability index)；
[0147]
步骤8、通过cudamemcpy函数将gpu计算得到的各生境因子数据及综合适宜指数csi复制到cpu中，计算目标物种有效栖息地加权可利用面积 wua(weighted usable area)；
[0148]
步骤8具体如下：
[0149]
根据综合适宜性指数csi计算出目标物种适宜有效栖息地加权可利用面积wua(weighted usable area)。与river2d模拟结果进行比较，以验证该方法的准确性；
[0150]
加权可利用面积wua计算式：
[0151][0152]
式中：wua为栖息地加权可利用面积，单位为m2；i为控制单元网格数； ai为第i个控制单元网格的面积，单位为m2。
[0153]
步骤9、根据计算得到的水深、流速和底质数据，以及水深、流速和底质三种生境因子的综合适宜指数csi数据，利用matlab软件输出水深、流速、底质和综合适宜性指数csi分布图；至此，耦合生境适宜模型的高效高精度栖息地模拟方法构建完成。
[0154]
以黄河上游某河段1～7km为研究河段，如图1所示，河段高分辨率地形空间离散成404,544个矩形网格，单元网格边长为5m，河道曼宁糙率系数取值为0.05s/m
1/3
。选择花斑裸鲤为目标鱼种，水深、流速和底质为花斑裸鲤的生境因子。由水生专家调查研究知，花斑裸
鲤产卵期适宜水深为0.1～ 2.0m，最适宜水深为0.4～1.2m；适宜流速为0.25～1.5m/s，最适宜流速为 0.5～1.0m/s；适宜基质条件为砂砾底质。根据花斑裸鲤产卵期水深、流速和底质三种生境因子条件，采用偏好适宜法绘制生境适宜性曲线，生境适宜指数hsi范围为0～1。其中最适宜生境适宜指数为1，最不适宜生境适宜指数为0。从而得水深、流速和底质三种生境因子的生境适宜曲线。
[0155]
图2(a)是水深适宜性曲线图。该图根据上述调查产卵期适宜水深条件；利用适宜性指数的规定原则，最适宜的产卵水深适宜性指数为1.0，最不适宜产卵水深适宜性指数为0.0，之间为0～1.0.从得到的水深适宜性曲线图。图 2(b)是流速适宜性曲线图。该图根据上述调查产卵期适宜流速条件，利用适宜性指数的规定原则，最适宜的产卵流速适宜性指数为1.0，最不适宜产卵流速适宜性指数为0.0，之间为0～1.0，从而得到流速适宜性曲线图。图2 (c)是底质适宜性曲线图。该图根据河段基质条件多为砂砾底质，又知花斑裸鲤产卵适宜基质条件为砂砾底质，利用适宜性指数的规定原则，底质的适宜性指数为1.0。
[0156]
通过将构建的生境适宜模型与高精度水动力过程相耦合，得到一种耦合生境适宜模型的高效高精度栖息地模拟方法，如图3和图4所示。为验证耦合生境适宜模型的高效高精度栖息地模拟方法的计算结果准确性，选取流量 26.5m3/s、53m3/s、74m3/s、79.5m3/s、106m3/s、159m3/s、212m3/s、318m3/s、 434m3/s、530m3/s、689m3/s、848m3/s、1060m3/s，分别采用gast模型和 river2d模型计算花斑裸鲤有效栖息地面积(以加权可用面积表示)。采用两种模型模拟不同流量条件下花斑裸鲤有效栖息地面积，结果如图5所示。由图5可知，在两种模型模拟下，花斑裸鲤产卵的有效栖息地面积均表现为：随着流量增加，有效栖息地面积不断增大，达到某一流量时，有效栖息地面积出现峰值。之后，随着流量增加，有效栖息地面积不断减少。根据对比图两种模型模拟结果几乎相同，规律一致，计算得到有效栖息地面积误差最大仅为5.366％。表明gast模型结果具有较好的准确性。
[0157]
此外，选择流量74m3/s，得出该流量下水深和流速分布图，如图6和图7所示。结合生境适宜模型，考虑花斑裸鲤产卵期水深、流速和底质三种生境因子，计算得花斑裸鲤产卵期的综合适宜指数csi，csi的分布如图8所示。发现在流量74m3/s时，花斑裸鲤产卵期有效栖息地面积达到最大，适宜产卵繁殖的有效栖息地分布区域向河道中间聚拢。且在流量74m3/s时，采用gast模型和river2d模型模拟得到的加权可利用面积wua分别为 202492.08m2和198775.00m2。发现gast模型得到的有效产卵栖息地面积略大，这是由于gast模型利用高分辨率地形条件模拟得到的水深和速度更为准确。