一种面向变动回水区鱼类产卵场生境重塑的生态调度方法

1.本发明属于水库管理技术领域，涉及鱼类产卵场生态调度，尤其涉及一种面向变动回水区鱼类产卵场生境重塑的生态调度方法。


背景技术：

2.水库群在发挥其防洪、发电、灌溉等兴利功能的同时，也改变了河流的天然水文情势，使得河流的地形、水动力、水温等环境状况发生了变化，进而破坏了水生生物的原有繁殖生长条件。水库调节性能越强，对河流的天然水文情势影响越大，以鱼类为代表的水生生物的生活习性也将因此改变。为了降低水库的这种不利影响，需要进行生态调度研究。但目前研究水库建设对鱼类产卵场的影响大都集中聚焦在坝址下游河段，对水库变动回水区关注甚少，使得水库变动回水区的产卵场被弃用。由于水库的可调度的特殊性，水库变动回水区的产卵场并不会完全消失，最具代表性的是长江三峡水库的变动回水区鱼类产卵场生境研究，在进行水库野外观测时，在原有天然河道产卵场附近发现了鱼卵，这表明鱼类仍然会在适宜其产卵的变动回水区进行繁殖。
3.目前对变动回水区鱼类产卵场生境的研究仅针对平原地区水库的四大家鱼，但四大家鱼对水动力条件的适宜性范围较大，且主要是分析大坝建立前后的鱼类有效栖息地面积(weighted usable area)(以下简称为wua)变化，定性评价水库运行管理要兼顾考虑变动回水区产卵场。这些研究方法很难满足山区河道型水库变动回水区以及例如齐口裂腹鱼之类的底栖性冷水鱼类的产卵场定量生态调度。


技术实现要素：

4.针对缺乏山区河道型水库变动回水区生态调度，尤其是缺乏针对底栖性冷水鱼类产卵场定量生态调度现状，本发明旨在提供一种基于水库变动回水区鱼类产卵场确定水库生态调度的方法，能够对不同工况下鱼类有效栖息地面积进行预测，进而得到水库变动回水区鱼类有效栖息地面积(wua)和水库调度的关系，从而为水库的运行调度提供技术支持，有利于水库管理和变动回水区鱼类产卵场保护。
5.本发明的发明思路为，首先建立水库模型，模拟天然河道的wua，进而确定目标wua值，然后再模拟建库以后不同工况下的wua，进而明确能够使变动回水区产卵场满足目标wua值的工况，根据筛选出的工况进行水库调度。
6.基于上述发明思路，本发明提供了一种面向变动回水区鱼类产卵场生境重塑的生态调度方法，其包括以下步骤：
7.步骤一，建立包含水库变动回水区的水库深度二维平面数值模型，记作二维水动力学模型；
8.步骤二，利用二维水动力学模型计算天然河道多年平均流量不同百分比下变动回水区的底栖性冷水鱼类有效栖息地面积wua，以确定目标wua值；
9.步骤三，依据历史水文数据，确定底栖性冷水鱼类产卵期时出现过的流量和水位
工况；
10.步骤四，利用二维水动力学模型计算步骤三中不同流量和水位工况下变动回水区的水动力条件，并根据鱼类流速水深适宜性曲线计算变动回水区的wua值；
11.步骤五，根据满足目标wua值以上的流量和水位工况，确定水库水位调度方案。
12.上述步骤一，收集水库地形资料、水文资料和鱼类流速水深适宜性曲线，水文资料包括入库流量、出库流量、坝前水位，流速水深适宜性曲线从文献中查得，包括流速适宜性曲线和水深适宜性曲线。收集的水文资料用于模型验证和统计产卵期水库出现过的流量水位情况。统计产卵期出现的流量水位工况将可为后续明确计算工况，所以应尽可能多的收集水文资料以全面反映水库变动回水区鱼类wua对水文资料的响应。
13.利用收集的地形资料建立水库深度二维平面数值模型；利用收集的水文资料对建立的数值模型进行参数率定，确保模型模拟的水动力条件是可靠的。
14.具体操作为：利用记事本将地形资料和边界资料转换成xyz格式文件，将xyz格式文件导入mike zero中的mesh generator模块生成计算网格和地形插值，生成mesh文件并导入。在mike21中选择flow mode模块，导入mesh文件，设置模拟时间、模拟步长、干湿边界、密度、涡粘系数、底床糙率、风场、降雨量、蒸发量、流量和水位数据，运行生成m21fm格式的文件，完成水动力学模型构建；生成的二维水动力学模型中变动回水区网格数量占整个水库网格数量的70～98％。根据收集的水文资料对水动力学模型进行率定，模拟值和实测值的误差控制在10％以内。
15.上述步骤二，模拟水库变动回水区天然河道情况下多年平均流量不同百分比下变动回水区的wua值，得到流量-wua曲线，根据曲线确定最大wua值，将wua最大值的60％作为目标wua值，具体步骤如下：
16.(1)利用二维水动力学模型模拟包含水库变动回水区的水库天然河道情况下多年平均流量不同百分比下的水动力条件；所述水动力条件包括流速和水深；
17.(2)提取各计算单元不同百分比情况下的流速和水深数据；
18.(3)从流速水深适宜性曲线查得每个流速水深数据对应的适宜度，并用几何平均法计算每个计算单元的综合栖息地适宜性指数hsi；hsi计算公式为：
[0019][0020]
其中，hsii为计算单元的综合栖息地适宜性指数，由si
vi
，si
di
和si
ci
(即流速、水深、河道适宜度(si
ci
综合考虑了底质和覆盖物状况))组合而成，这三者是河流物理栖息地最具有代表性的变量，三者的取值范围均为0～1；对于已经明确可以作为鱼类产卵场的水库变动回水区，计算中河道适宜指数si
ci
均为1；
[0021]
(4)计算变动回水区的鱼类有效栖息地面积wua；wua是指鱼类可在河流中生存的有效栖息地面积，定义为计算单元的栖息地适宜性指数和该单元面积的乘积，计算公式为：
[0022][0023]
其中，n为水库变动回水区的计算单元个数；ai为计算单元的面积；由此计算出各个工况下变动回水区的鱼类有效栖息地面积wua；
[0024]
(5)依据多年平均流量不同百分比对应流量及相应的鱼类有效栖息地面积wua确
定流量-wua曲线，并根据曲线得到最大wua值，将wua最大值的60％作为目标wua值。
[0025]
本发明所针对的底栖性冷水鱼类是指高原地区特产、常生活于急流中的裂腹鱼，例如齐口裂腹鱼、重口裂腹鱼等。
[0026]
上述步骤三，根据收集到的水文资料，统计底栖性冷水鱼类产卵期出现的来流流量和水库水位情况，包括流量阈值范围和水位阈值范围。对于流量而言，从最低流量开始，每间隔20-100m3/s设置一个流量工况，直至最高流量；对水位而言，从最低水位开始，每间隔2-5米设置一个水位工况，直至最高水位；利用流量工况和水位工况构建若干流量水位组合。
[0027]
上述步骤四，利用二维水动力学模型计算步骤三中不同来流流量和水库水位工况下变动回水区的水动力条件，并确定水库变动回水区的鱼类有效栖息地面积wua值，具体步骤参见步骤二中的(2)-(4)。
[0028]
上述步骤五，结合目标wua值以上工况，筛选出满足目标wua以上的流量水位组合情况；根据鱼类产卵期来流情况(即水库来流流量)，确定水库水位范围，进而调度水库水位以使变动回水区的wua在目标值以上。
[0029]
与现有技术相比，本发明提供的面向变动回水区鱼类产卵场生境重塑的生态调度方法，具有以下有益效果：
[0030]
本发明通过整合二维水动力学模型和鱼类栖息地模型，建立水库流量-水位-鱼类有效栖息地面积三者的关系，明确水库生态调度措施；通过本发明提供的生态调度方法可有效改善河流生态，为维护重建变动回水区产卵场的水库提供有效参考。
附图说明
[0031]
图1为面向变动回水区鱼类产卵场生境重塑的生态调度方法的流程示意图。
[0032]
图2为本发明实施例选取的岷江齐口裂腹鱼产卵期流速水深适宜性曲线。
[0033]
图3为本发明实施例不同工况下变动回水区wua图；其中，(a)为不同流量、水位情况下的wua变化图；(b)为满足目标wua以上的流量-水位关系图。
[0034]
本发明中，水位是指水面高于特征海平面(常用黄海基面)的高程；水深是指水面到河底的距离。
具体实施方式
[0035]
下面将结合附图对本发明进行进一步说明。
[0036]
实施例1
[0037]
岷江是长江中上游的一条二级支流，其与中国西南山区众多河流具有相同的水流湍急、坡降大等特点，岷江干流也同样进行了梯级水电开发，河流的水生生境遭到破坏，库区产卵场被淹没。紫坪铺水库为岷江干流中最大的水利工程，水库的正常蓄水位为877m，死水位为817m，水位变幅达60m。紫坪铺水库的变动回水区长12km，据调查，水库蓄水前的库区古溪沟村附近有齐口裂腹鱼(属于底栖性冷水鱼)产卵场存在，然而建库后的回水淹没及水位变幅影响造成急流生境消失，使得原有的齐口裂腹鱼等急流性鱼类的存活机会及存活空间大幅减小。齐口裂腹鱼是西南山区河流中的典型性优势鱼类，也是四川省省级保护动物，其喜好在水温较低、溶解氧较高的急、缓流交界处生活，产卵季节有短距洄游需求，上述习
性与库区环境并不相符使得人们少有关注其生活动态。2022年5月，对紫坪铺水库开展野外观测时发现原有产卵场位置出现鱼类产卵繁殖迹象，这表明了变动回水区内仍存在适宜鱼类产卵的生境，故本实施例以紫坪铺水库为例，通过本发明提供的生态调度方法重塑水库变动回水区齐口裂腹鱼的产卵生境。
[0038]
本实施例提供的面向变动回水区鱼类产卵场生境重塑的生态调度方法，如图1所示，其包括以下步骤：
[0039]
步骤一，建立包含水库变动回水区的水库深度二维平面数值模型，记作二维水动力学模型。
[0040]
本步骤中，收集水库地形资料、水文资料和鱼类流速水深适宜性曲线。水文资料包括入库流量、出库流量、坝前水位，流速水深适宜性曲线从文献中查得。
[0041]
利用收集的地形资料建立水库深度二维平面数值模型；利用收集的水文资料对建立的数值模型进行参数率定，确保模型模拟的水动力条件是可靠的。具体操作为：
[0042]
(1)利用记事本将地形资料和边界资料转换成xyz格式文件，将xyz格式文件导入mike zero中的mesh generator模块生成计算网格和地形插值，生成mesh文件并导入；在mike21中选择flow mode模块，导入mesh文件，设置模拟时间、模拟步长、干湿边界、密度、涡粘系数、底床糙率、风场、降雨量、蒸发量、流量和水位数据，运行生成m21fm格式的文件，完成水动力学模型构建；
[0043]
(2)根据收集的水文资料对水动力学模型进行率定，模拟值和实测值的误差控制在10％以内。
[0044]
本实施例中，紫坪铺水库是西南山区典型的河道型水库，变动回水区坡降大，根据具体实施方式步骤一从紫坪铺水库管理处收集到2016年地形资料、水文资料和流速水深适宜性曲线等基础数据，水文资料包括50年(1971-2020年)水库坝址处的逐旬流量和10年(2011-2020年)间紫坪铺水库齐口裂腹鱼产卵期(3月～5月)实测的逐日流量及水位情况。岷江齐口裂腹鱼产卵期流速水深适宜性曲线采用陈明千在岷江研究获得的成果，岷江中齐口裂腹鱼产卵的流速适宜区间为0.5～2.5m/s，其中1.4～1.6m/s为最优适宜区间，适宜度为1；对于水深而言，适宜区间为0.5～1.5m，详细流速与水深适宜度见图2所示。
[0045]
按照步骤一构建紫坪铺水库水动力mike21数值模型(即二维水动力学模型)并进行率定。首先将收集的2016年紫坪铺水库地形资料转化为平面坐标和高程，导入mike21中生成23452个非结构网格，其中在变动回水区内构建了面积不大于700m2的22391个加密网格(生成的二维水动力学模型中变动回水区网格数量约占整个水库网格数量的95％)，建立的紫坪铺水库深度平均二维mike21数值模型(即二维水动力学模型)用于计算评估水动力条件。利用2016年3月主库岷江入库流量、支库寿溪河入库流量和水库出库流量模拟每日坝前水位，比较模型计算得到的坝前水位模拟值和坝前水位实测值，模型模拟值与实测值的最大绝对误差为0.32m，满足误差要求。
[0046]
步骤二，利用二维水动力学模型计算天然河道多年平均流量不同百分比下变动回水区的底栖性冷水鱼类有效栖息地面积wua，以确定目标wua值。
[0047]
本步骤具体包括以下分步骤：
[0048]
(1)利用二维水动力学模型模拟包含水库变动回水区的水库天然河道情况下多年平均流量不同百分比下的水动力条件；所述水动力条件包括流速和水深。
[0049]
紫坪铺水库坝址处多年平均流量为469m3/s，百分比从15％逐渐增加至100％的过程中设置了10种工况进行计算，分别是15％、30％、40％、50％、55％、60％、65％、70％、80％和100％。
[0050]
基于设置的10种不同流量工况，利用二维水动力学模型进行模拟，获取相应工况下的水动力条件。
[0051]
(2)提取各计算单元不同百分比情况下的流速和水深数据。
[0052]
根据步骤(1)模拟结果，导出每一种流量工况下的流速和水深数据。
[0053]
(3)从流速水深适宜性曲线查得每个流速水深数据对应的适宜度，并用几何平均法计算每个计算单元的综合栖息地适宜性指数hsi；hsi计算公式为：
[0054][0055]
其中，hsii为计算单元的综合栖息地适宜性指数，由si
vi
，si
di
和si
ci
(即流速、水深、河道适宜度，si
ci
综合考虑了底质和覆盖物状况)组合而成，这三者是河流物理栖息地最具有代表性的变量，三者的取值范围均为0～1；对于已经明确可以作为鱼类产卵场的水库变动回水区，计算中河道适宜指数si
ci
均为1。
[0056]
根据图2查得每个计算单元的流速水深对应的适宜度，然后按照上述公式计算每个计算单元的综合栖息地适宜性指数hsi。
[0057]
(4)计算变动回水区的鱼类有效栖息地面积wua；wua是指鱼类可在河流中生存的有效栖息地面积，定义为计算单元的栖息地适宜性指数和该单元面积的乘积，计算公式为：
[0058][0059]
其中，n为水库变动回水区的计算单元个数；ai为计算单元的面积；由此计算出各个工况下变动回水区的鱼类有效栖息地面积wua。
[0060]
根据wua公式计算出每一种流量工况下变动回水区的wua值，进而得到流量-wua曲线。
[0061]
(5)依据多年平均流量不同百分比对应流量及相应的鱼类有效栖息地面积wua确定流量-wua曲线，并根据曲线得到最大wua值，将wua最大值的60％作为目标wua值。
[0062]
根据流量-wua曲线，百分比从15％逐渐增加至100％的过程中，wua先增大后减少。65％的多年平均流量(即304.85m3/s)时，wua达到峰值634889m2，故而紫坪铺水库变动回水区目标wua即为wua峰值的60％，也就是380933m2。
[0063]
步骤三，依据历史水文数据，确定底栖性冷水鱼类产卵期时出现过的流量和水位工况，包括流量阈值范围和水位阈值范围。对于流量而言，从最低流量开始，每间隔20-100m3/s设置一个流量工况，直至最高流量；对水位而言，从最低水位开始，每间隔3米设置一个水位工况，直至最高水位；利用流量工况和水位工况构建若干流量水位组合。
[0064]
统计收集到的2011-2020年的10年间紫坪铺水库齐口裂腹鱼产卵期实测的逐日流量及水位情况，依据统计结果流量设计变化范围为50-900m3/s，当流量高于900m3/s时，由于流速过大，变动回水区wua值较低，明显不适宜营造产卵场生境；水位设计变化范围为817-856m，故确定13种流量梯度和14种水位梯度如表1所示。13种流量梯度和14种水位梯度共组合成182组工况开展模拟计算变动回水区的水动力条件。
[0065]
表1流量和水位工况
[0066][0067]
步骤四，利用二维水动力学模型计算步骤三中不同流量和水位工况下变动回水区的水动力条件，并根据鱼类流速水深适宜性曲线计算变动回水区的wua值。
[0068]
本步骤具体包括以下分步骤：
[0069]
(1)利用二维水动力学模型模拟包含水库变动回水区的水库天然河道步骤三中不同流量和水位工况下的水动力条件；所述水动力条件包括流速和水深。
[0070]
基于设置的182组不同流量和水位工况，利用二维水动力学模型进行模拟，获取相应工况下的水动力条件。
[0071]
(2)提取各计算单元不同百分比情况下的流速和水深数据。
[0072]
根据步骤(1)模拟结果，导出每一种流量、水位工况下的流速。
[0073]
(3)从流速水深适宜性曲线查得每个流速水深数据对应的适宜度，并用几何平均法计算每个计算单元的综合栖息地适宜性指数hsi；hsi计算公式为：
[0074][0075]
其中，hsii为计算单元的综合栖息地适宜性指数，由si
vi
，si
di
和si
ci
(即流速、水深、河道适宜度(si
ci
综合考虑了底质和覆盖物状况))组合而成，这三者是河流物理栖息地最具有代表性的变量，三者的取值范围均为0～1；对于已经明确可以作为鱼类产卵场的水库变动回水区，计算中河道适宜指数si
ci
均为1。
[0076]
根据图2查得每个计算单元的流速水深对应的适宜度，然后按照上述公式计算每个计算单元的综合栖息地适宜性指数hsi。
[0077]
(4)按照以下公式计算变动回水区的鱼类有效栖息地面积wua：
[0078][0079]
其中，n为水库变动回水区的计算单元个数；ai为计算单元的面积；由此计算出各个工况下变动回水区的鱼类有效栖息地面积wua。
[0080]
以紫坪铺水库入库流量为300m3/s，坝前水位为841m为例，模型部分计算单元流速水深对应适应性如表2所示。
[0081]
表2计算单元流速水深对应适应性
[0082]
[0083][0084]
对于已经明确可以作为鱼类产卵场的水库变动回水区，计算中河道适宜指数均为1。根据hsi公式计算入库流量300m3/s、坝前水位841m时的各计算单元的hsii，结果见表3。然后根据wua公式水库的有效栖息地面积(即wua值)。
[0085]
表3计算单元hsi
[0086]
[0087][0088]
最终计算得到紫坪铺水库入库流量300m3/s、坝前水位841m时变动回水区的wua为414667.5185m2。
[0089]
其他工况按上述工况计算，最终得到不同工况下紫坪铺变动回水区的wua结果见图3。
[0090]
步骤五，根据满足目标wua值以上的流量和水位工况，确定水库水位调度方案。
[0091]
本步骤的具体实施过程如下：
[0092]
(1)整理不同流量水位组合情况下水库变动回水区的wua变化，筛选出满足目标
wua以上的流量水位组合情况；
[0093]
(2)根据水库来流流量，确定水库水位范围，调度水库水位使水库变动回水区的wua满足目标wua。
[0094]
本实施例中，由图3(a)图可见，当水库水位固定时，随着来流流量从50m3/s增大至900m3/s，水库变动回水区的wua先增大到峰值后逐渐减小。不同水位工况下均显示水库变动回水区的wua在流量300m3/s时达到峰值，并在300m3/s后稳步下降，有着极高的一致性。而当来流流量固定时，随着水位从817m升至856m，水库变动回水区的wua呈现显著单调性，水位越高，wua越小。图中白色平面为目标wua值380933m2，白色平面以上即满足目标wua，水库水位在844m以上时，wua值均在白色平面以下，不能达到目标wua。所以对紫坪铺水库而言，当水位高于844m时，不管来流情况如何，均不能营造目标值以上的wua，故齐口裂腹鱼产卵期应尽可能使水库水位在844m以下运行，但需注意不应低于紫坪铺水库设计的死水位817m。
[0095]
提取满足目标wua值的流量水位工况(图3(b)图)。图中颜色越深则表明水库变动回水区的wua越高，反之颜色越浅则表明水库变动回水区的wua越低，斜线区域为目标wua值以上。可见水库变动回水区的wua高值集中在流量300m3/s附近，且水位越低，wua越大。来流流量低于114.72m3/s或高于652.64m3/s时，无法通过调度水库水位满足目标wua值。来流流量介于[260m3/s，370m3/s]时，水库水位在844m以下均可满足目标wua，但需注意不应低于水库设计的死水位817m。
[0096]
紫坪铺水库在齐口裂腹鱼产卵期可根据来流流量调度水位使wua值落于图3斜线区域，以达到在变动回水区营造足够鱼类wua的目的。
[0097]
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

技术特征：
1.一种面向变动回水区鱼类产卵场生境重塑的生态调度方法，其特征在于，包括以下步骤：步骤一，建立包含水库变动回水区的水库深度二维平面数值模型，记作二维水动力学模型；步骤二，利用二维水动力学模型计算天然河道多年平均流量不同百分比下变动回水区的底栖性冷水鱼类有效栖息地面积wua，以确定目标wua值；步骤三，依据历史水文数据，确定底栖性冷水鱼类产卵期时出现过的流量和水位工况；步骤四，利用二维水动力学模型计算步骤三中不同流量和水位工况下变动回水区的水动力条件，并根据鱼类流速水深适宜性曲线计算变动回水区的wua值；步骤五，根据满足目标wua值以上的流量和水位工况，确定水库水位调度方案。2.根据权利要求1所述的面向变动回水区鱼类产卵场生境重塑的生态调度方法，其特征在于，生成的二维水动力学模型中变动回水区网格数量占整个水库网格数量的70-98％。3.根据权利要求1或2所述的面向变动回水区鱼类产卵场生境重塑的生态调度方法，其特征在于，步骤二中，模拟水库变动回水区天然河道情况下多年平均流量不同百分比下变动回水区的wua值，得到流量-wua曲线，根据曲线确定最大wua值，将wua最大值的60％作为目标wua值，具体步骤如下：(1)利用二维水动力学模型模拟包含水库变动回水区的水库天然河道情况下多年平均流量不同百分比下的水动力条件；所述水动力条件包括流速和水深；(2)提取各计算单元不同百分比情况下的流速和水深数据；(3)从流速水深适宜性曲线查得每个流速水深数据对应的适宜度，并用几何平均法计算每个计算单元的综合栖息地适宜性指数his；hsi计算公式为：其中，hsi
i
为计算单元的综合栖息地适宜性指数，si
vi
，si
di
和si
ci
分别表示流速、水深、河道适宜度；(4)计算变动回水区的鱼类有效栖息地面积wua；wua是指鱼类在河流中生存的有效栖息地面积，定义为计算单元的栖息地适宜性指数和该单元面积的乘积，计算公式为：其中，n为水库变动回水区的计算单元个数；a
i
为计算单元的面积；由此计算出各个工况下变动回水区的鱼类有效栖息地面积wua；(5)依据多年平均流量不同百分比对应流量及相应的鱼类有效栖息地面积wua确定流量-wua曲线，并根据曲线得到最大wua值，将wua最大值的60％作为目标wua值。4.根据权利要求3所述的面向变动回水区鱼类产卵场生境重塑的生态调度方法，其特征在于，步骤三中，对于流量而言，从最低流量开始，每间隔20-100m3/s设置一个流量工况，直至最高流量；对水位而言，从最低水位开始，每间隔2-5米设置一个水位工况，直至最高水位；利用流量工况和水位工况构建若干流量水位组合。5.根据权利要求3所述的面向变动回水区鱼类产卵场生境重塑的生态调度方法，其特征在于，步骤五中，结合目标wua值以上工况，筛选出满足目标wua以上的流量水位组合情
况；根据鱼类产卵期来流情况，确定水库水位范围，进而调度水库水位以使变动回水区wua在目标值以上。

技术总结
本发明公开了一种面向变动回水区鱼类产卵场生境重塑的生态调度方法，建立水库变动回水区深度二维平面数值模型，记作二维水动力学模型；模拟变动回水区天然河道时多年平均流量下的鱼类有效栖息地面积(WUA)，进而确定目标WUA值，然后再模拟建库以后不同流量、水位工况下的WUA，进而明确能够使变动回水区产卵场满足目标WUA值的工况，针对来流情况选定满足营造需求的水位区间进行水位调度。本发明通过整合二维水动力学模型和鱼类栖息地模型，建立水库流量-水位-鱼类有效栖息地面积三者的关系，明确水库生态调度措施。通过本发明提供的生态调度方法可有效改善河流生态，为维护重建变动回水区产卵场的水库提供有效参考。回水区产卵场的水库提供有效参考。回水区产卵场的水库提供有效参考。


技术研发人员：张陵蕾 李福霜 陈旻 李嘉 廖宁 安瑞冬 邓云 脱友才 方思倩 李东海
受保护的技术使用者：四川大学
技术研发日：2023.04.28
技术公布日：2023/8/4