CN109615076A - 一种面向鱼类生境保护的河流生态流量过程推求方法 - Google Patents
一种面向鱼类生境保护的河流生态流量过程推求方法 Download PDFInfo
- Publication number
- CN109615076A CN109615076A CN201811523739.5A CN201811523739A CN109615076A CN 109615076 A CN109615076 A CN 109615076A CN 201811523739 A CN201811523739 A CN 201811523739A CN 109615076 A CN109615076 A CN 109615076A
- Authority
- CN
- China
- Prior art keywords
- fish
- habitat
- water
- interest
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 241000251468 Actinopterygii Species 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000008569 process Effects 0.000 title claims abstract description 28
- 238000004364 calculation method Methods 0.000 title claims abstract description 19
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 93
- 230000004044 response Effects 0.000 claims abstract description 33
- 238000002474 experimental method Methods 0.000 claims abstract description 10
- 238000011160 research Methods 0.000 claims abstract description 8
- 230000001932 seasonal effect Effects 0.000 claims abstract description 5
- 208000016117 Idiopathic hemiconvulsion-hemiplegia syndrome Diseases 0.000 claims description 10
- 230000017448 oviposition Effects 0.000 claims description 10
- 238000012216 screening Methods 0.000 claims description 8
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 7
- 239000001301 oxygen Substances 0.000 claims description 7
- 229910052760 oxygen Inorganic materials 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 102000002322 Egg Proteins Human genes 0.000 claims description 3
- 108010000912 Egg Proteins Proteins 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 210000004681 ovum Anatomy 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims 1
- 230000007613 environmental effect Effects 0.000 abstract description 10
- 238000011835 investigation Methods 0.000 abstract description 5
- 239000007787 solid Substances 0.000 abstract description 2
- 241001269238 Data Species 0.000 abstract 1
- 241000894007 species Species 0.000 description 6
- 206010058109 Hangnail Diseases 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000000975 bioactive effect Effects 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 239000003643 water by type Substances 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000009605 growth rhythm Effects 0.000 description 1
- 230000000366 juvenile effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/02—Computing arrangements based on specific mathematical models using fuzzy logic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/04—Inference or reasoning models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Marine Sciences & Fisheries (AREA)
- Biomedical Technology (AREA)
- Automation & Control Theory (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Fuzzy Systems (AREA)
- Mining & Mineral Resources (AREA)
- Computational Linguistics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Revetment (AREA)
Abstract
本发明公开了一种面向鱼类生境保护的河流生态流量过程推求方法，所述方法具体包含：首先对研究区域鱼类资源现状进行调查，确定生态流量推求的目标鱼类；对目标鱼类进行针对性调查，同时进行实验室胁迫实验，筛选关键水环境因子，利用实验室及野外调查获得的数据，建立目标鱼类不同生命阶段与关键环境因子间的定量响应关系；收集研究区域的水文、水质、地形及气象等基础数据，采用模糊数学原理建立鱼类栖息地模型，并与水环境模型耦合，构建鱼类栖息地评价体系，得出栖息地与流量的动态响应关系；确定河流生态流量过程计算的参照体系，综合考虑目标鱼类不同生命阶段的生境需求及天然径流的季节变化，推求河流完整水文年生态流量过程。该发明可为水库面向生态友好调度提供约束条件，并为流域水资源分配及管理提供切实可靠的依据。
Description
技术领域
本发明涉及水利技术领域，尤其涉及一种面向鱼类生境保护的河流生态流量过程推求方法。
背景技术
世界上50％的可利用水资源已被人类开发，到2025年可能会增加至70％。拦河闸坝是人类开发利用河流水资源的重要方式之一。目前，我国过闸流量1m3/s以上的水闸有268476座，已建库容10万m3以上的水库高达98002座。闸坝建设人为改变了河流原有的物质场、能量场、化学场和生物场，直接影响生源要素在河流中的生物地球化学行为(生源要素输送通量、赋存形态、组成比例等)，改变河流生态系统的物种构成、栖息地分布以及相应的生态功能，进而造成河流生态系统退化。随着闸坝建成后生态累积效应日益凸显，如何减缓闸坝运行对河流生态环境的负面效应已经成为水资源开发利用中迫切需要解决的问题。
生态流量是保证河流生态功能的最基本要素，是实现河流生态友好调度及水资源可持续开发利用的基础。但是生态流量作为生态修复的重要手段还处于初步研究阶段，推求方法多采用水文法或水力学法，以河流生态系统需求为基础的生态流量推求方法相对较少，尤其在完整水文年内生态流量过程的推求尚缺乏较为完善的计算方法。
发明内容
发明目的：本发明提供一种面向鱼类生境需求的河流生态流量过程推求方法，为水库面向生态友好调度提供约束条件，并为流域水资源分配及管理提供切实可靠的依据。
技术方案：本发明所述的面向鱼类生境需求的河流生态流量过程推求方法，包括：
(1)对研究区域鱼类资源现状进行调查，确定生态流量推求的目标鱼类；
(2)筛选关键水环境因子，并利用胁迫实验建立目标鱼类不同生命阶段与关键水环境因子间的响应关系；
(3)根据研究区域的水文、水质、地形及气象基础数据建立水动力水环境模型，并结合目标鱼类对关键生境因子的定量响应关系曲线，采用模糊数学原理建立鱼类栖息地模型；
(4)建立鱼类栖息地评价体系，对鱼类栖息地模型的模拟结果进行评价，并构建评价指标与流量的动态响应关系；
(5)确定平原河流生态流量计算的参照体系，并依照该参照体系，综合考虑鱼类不同生命阶段的生境需求及天然径流的季节变化，推求河流完整水文年生态流量过程。
进一步的，所述步骤(1)具体包括：
(1-1)鱼类资源现状调查：依据研究区域鱼类资源现状及历史资料，确定河流生态流量推求备选鱼类，其中，鱼类资源现状调查包括调查鱼类的种群数量，测定鱼类群体的结构状况，掌握渔业利用情况，研究鱼类种群和数量变动规律；
(1-2)目标鱼类筛选：分析鱼类生物学特性，筛选对水环境要求具有代表性的鱼类作为目标鱼类。
进一步的，所述步骤(2)具体包括：
(2-1)关键水环境因子筛选：针对筛选出的目标鱼类，根据鱼类生物学特性选择流速、水深、温度及溶氧为目标鱼类的关键水环境因子；
(2-2)响应关系曲线建立：将鱼类生命阶段分为为产卵期、捕食期及越冬期，利用胁迫实验建立目标鱼类不同生命阶段对关键水环境因子的定量响应关系。
进一步的，所述步骤(3)具体包括：
(3-1)水动力水环境模型建立：收集研究区域的水文、水质、地形及气象基础数据，根据所述基础数据采用开源模型SELFE建立水动力水环境模型，其中，建立的水动力水环境模型的控制方程为：
式中，Qa为流量，单位m3/s；t表示时间，单位为s；H为水位，单位m；h为水深，单位m，u,v为x,y方向上的速度，单位m/s；p为大气压强，单位N/m2；ρ0为水体密度，值为1025kg/m3；ν为水平粘性系数，单位m2/s；f为科氏力系数；τx,τy为底部剪应力，单位N/m；c为物质浓度，单位kg/m3；Dx,Dy为x,y方向扩散系数，单位m2/s；S为源汇项；fR(c,t)为反应项；
(3-2)鱼类栖息地模型建立：将目标鱼类不同生命阶段对关键水环境因子的响应关系和水动力水环境模型做为输入参数，采用模糊数学原理建立鱼类栖息地模型。
进一步的，所述步骤(4)具体包括：
(4-1)鱼类栖息地评价体系建立：采用水力栖息地适应性指数IHHS和栖息地连通性指数HCI作为鱼类栖息地的评价指标；其中：
水力栖息地适应性指数IHHS的计算公式为：
式中，Ai表示对研究区域划分出的网格计算单元i的面积，单位为m2；HSIi为网格计算单元i的栖息地适应性指数，N为网格计算单元数量。
栖息地连通指数HCI的计算公式为：式中，Ai表示对研究区域划分出的网格计算单元i的面积，单位为m2；HSIi为网格计算单元i的栖息地适应性指数，N为网格计算单元数量。
栖息地连通指数HCI的计算公式为：
其式中，M为连通性系数；MDPi为斑块i与其他斑块间最短距离，单位为m；，定义为适应性大于一定阈值(HSI≥HSIcut)且面积大于一定阈值(P≥Pcut)的栖息地斑块为有效斑块pe，单位m2，ne为有效斑块个数，N为网格计算单元数量；
(4-2)将流量设置为不同值，重复步骤(1)-(4)，从而得到目标鱼类栖息地评价指标与流量的动态响应曲线。
(4-2)将流量设置为不同值，重复步骤(1)-(4)，从而得到目标鱼类栖息地评价指标与流量的动态响应曲线。
进一步的，所述步骤(5)具体包括：
(5-1)参照体系及鱼类生境需求确定：构建流量与IHHS和HCI的三维曲面，并选择三维曲面中的最优点相应的栖息地面积作为参照体系，并根据鱼类生物学特性将目标鱼类的产卵期、捕食期及越冬期与水文年的月份对应，且在一个水文年份中，产卵期优先考虑产卵场的需求，其次考虑索饵场；越冬期同时考虑越冬场及索饵场需求；其他时段考率满足鱼类基本生存条件的需求；
(5-2)河流生态流量过程计算：首先计算参照体系下的有效栖息地面积，使用有效栖息地面积恢复比例为指标，确定鱼类栖息地恢复的目标，其次采用建立的栖息地模型反算不同恢复目标对应的流量，并进行栖息地评价，从中选择合适的流量值，最后根据时间尺度和鱼类生活史，确定完整水文年内的流量过程。
有益效果：本发明与现有技术相比，其显著优点是：
(1)针对以往生态流量采用单一流量的不足，提出加入季节性变化及鱼类生活史影响因素，确定河流生态流量过程，比以往单一静态的生态流量更具生态价值。
(2)本发明提供的生态流量过程推求方法具有普适性，可以用于河道内其他生物的生态流量推求，将各典型物种的生态流量过程进行集合运算，生成包络线即可视为河流生物的综合生态流量。
附图说明
图1是本发明的一个实施例的流程示意图；
图2是目标鱼类对关键环境因子的定量响应关系示意图；
图3是某河流鱼类栖息地模拟结果与历史资料验证示意图；
图4是某河流鱼类栖息地评价指数--流量动态响应关系示意图；
图5是某河流流量—连续性—适应性曲面示意图；
图6是某河流面向鱼类生境保护的生态流量过程示意图。
具体实施方式
如图1所示，本实施例的面向河流生境保护的生态流量过程推求方法，包括以下步骤：
(1)对研究区域鱼类资源现状进行调查，确定生态流量推求的目标鱼类。具体包括：
(1-1)鱼类资源现状调查：依据研究区域鱼类资源现状及历史资料，确定河流生态流量推求备选鱼类，其中，鱼类资源现状调查包括调查鱼类的种群数量，测定鱼类群体的结构状况，掌握渔业利用情况，研究鱼类种群和数量变动规律。
(1-2)目标鱼类筛选：分析鱼类生物学特性，筛选对水环境要求具有代表性的鱼类作为目标鱼类。分析组成鱼群个体的年龄、体长、体重、摄食营养、繁殖、洄游分布规律等指标，依据鱼类生物学特性了解鱼类与水环境因子间的响应关系，筛选目标物种。
鱼类资源调查原则上要采集调查水域的全部种类。在每个采集点使用撒网、刺网、地笼等工具捕捞小型鱼类，并联系当地渔民出河捕捞；收集周围农贸市场中当地渔民捕获的野杂鱼类；对于非捕捞水体、非经济鱼类或者稀有和珍贵鱼类，则需要通过专门的采捕而获得。同时收集历史鱼类资源调查数据，最终结合鱼类生物学特性，筛选确定目标鱼类。通过以上方法，本实施例选择光倒刺鲃作为目标鱼类。
(2)筛选关键水环境因子，并利用胁迫实验建立目标鱼类不同生命阶段与关键水环境因子间的响应关系。具体包括：
(2-1)目标鱼类生物学特性研究：了解目标鱼类生命周期及生长规律，调查其产卵场、越冬场及幼鱼育肥场的水域分布位置、面积范围大小、洄游规律等，记录上述场所的水位、水深、温度、透明度、含氧量、流速及底质等环境条件，筛选关键水环境因子，为初步建立环境因子与目标鱼类的响应关系及栖息地模型验证提供数据支持。
确定目标种类后，针对光倒刺鲃的生物习性及产卵场和越冬场进行了调查，光倒刺鲃是一种重要的名优经济鱼类，一般栖息于底质多乱石而水流较湍急的江河中的中下层，其生长环境要求水质清新，最适生长温度为18～28℃。光倒刺鲃性成熟年龄为2.5～3年，每年4～5月间在水流缓慢、水草较多处产粘性卵。调查显示研究区段存在三个传统的光倒刺鲃产卵场。结合目标鱼类的生物特性及调查数据，选取溶氧、水温、流速及水深作为栖息地模型的关键环境因子。
(2-2)响应关系曲线建立：目标鱼类对关键水环境因子的定量响应关系的建立通过实验室水槽实验获得。实验通过观测目标鱼类的出现频率和关键环境因子，建立目标鱼类不同生命阶段对环境因子的定量响应关系。
确定研究物种及生理生态参数大概范围后，为进一步确定物种与水环境因子间的响应关系，在实验室中进行鱼类生长实验及相关环境因子胁迫。实验根据光倒刺鲃生活史分别建立其不同生命阶段(产卵期、捕食期及越冬期)对关键生境因子(溶氧、水温、流速、水深)的定量响应曲线(图2)。
(3)根据研究区域的水文、水质、地形及气象基础数据建立水动力水环境模型，并结合目标鱼类对关键生境因子的定量响应关系曲线，采用模糊数学原理建立鱼类栖息地模型。具体包括：
(3-1)水动力水环境模型建立：收集研究区域的水文、水质、地形及气象基础数据，根据所述基础数据采用开源模型SELFE建立水动力水环境模型，其中，建立的水动力水环境模型的控制方程为：
式中，Qa为流量，单位m3/s；t表示时间，单位为s；H为水位，单位m；h为水深，单位m，u,v为x,y方向上的速度，单位m/s；p为大气压强，单位N/m2；ρ0为水体密度，值为1025kg/m3；ν为水平粘性系数，单位m2/s；f为科氏力系数；τx,τy为底部剪应力，单位N/m；c为物质浓度，单位kg/m3；Dx,Dy为x,y方向扩散系数，单位m2/s；S为源汇项；fR(c,t)为反应项；
(3-2)鱼类栖息地模型建立：将目标鱼类不同生命阶段对关键水环境因子的响应关系和水动力水环境模型做为输入参数，采用模糊数学原理建立鱼类栖息地模型。鱼类栖息地模型基本框架采用模糊数学方法建立，具体步骤包括：输入因子的模糊化；模糊规则库的生成；模糊推理过程；解模糊化。
例如，选取我国西南部某河流典型河段建立水动力水环境模型。计算河段上游入口以日均流量作为边界条件，下游出口以日均水位作为控制条件，溶解氧和温度也以日均浓度作为输入条件，并以实测资料进行验证。基于模糊数学原理，将目标鱼类不同生命阶段(产卵期、捕食期及越冬期)对关键生境因子(流速、水深、温度及溶氧)定量响应关系及水动力水环境模型计算结果输入栖息地模型模拟不同流量下鱼类栖息地状况，并采用调查数据进行验证(图3)。
(4)建立鱼类栖息地评价体系，对鱼类栖息地模型的模拟结果进行评价，并构建评价指标与流量的动态响应关系。具体包括：
(4-1)鱼类栖息地评价体系建立：采用水力栖息地适应性指数IHHS和栖息地连通性指数HCI作为鱼类栖息地的评价指标；其中：
水力栖息地适应性指数IHHS的计算公式为：
式中，Ai表示对研究区域划分出的网格计算单元i的面积，单位为m2；HSIi为网格计算单元i的栖息地适应性指数，N为网格计算单元数量。
栖息地连通指数HCI的计算公式为：式中，Ai表示对研究区域划分出的网格计算单元i的面积，单位为m2；HSIi为网格计算单元i的栖息地适应性指数，N为网格计算单元数量。
栖息地连通指数HCI的计算公式为：
其式中，M为连通性系数；MDPi为斑块i与其他斑块间最短距离，单位为m；，定义为适应性大于一定阈值(HSI≥HSIcut)且面积大于一定阈值(P≥Pcut)的栖息地斑块为有效斑块pe，单位m2，ne为有效斑块个数，N为网格计算单元数量；
(4-2)将流量设置为不同值，重复步骤(1)-(4)，从而得到目标鱼类栖息地评价指标与流量的动态响应曲线。
即根据栖息地模拟结果计算研究区域水力栖息地适应性指数和连通性指数，得到水力栖息地适应性指数--流量关系图和栖息地连通性指数—流量指数关系图(图4)，作为栖息地评价指标与流量的动态响应关系，定量评价径流变化对鱼类栖息地的影响。
(5)确定平原河流生态流量计算的参照体系，并依照该参照体系，综合考虑鱼类不同生命阶段的生境需求及天然径流的季节变化，推求河流完整水文年生态流量过程。具体包括：
(5-1)参照体系及鱼类生境需求确定：构建流量与IHHS和HCI的三维曲面，并选择三维曲面中的最优点相应的栖息地面积作为参照体系，并根据鱼类生物学特性将目标鱼类的产卵期、捕食期及越冬期与水文年的月份对应，且在一个水文年份中，产卵期优先考虑产卵场的需求，其次考虑索饵场；越冬期同时考虑越冬场及索饵场需求；其他时段考率满足鱼类基本生存条件的需求。
如图5所示，利用流量与IHHS和HCI的三维曲面，选择最优点作为参考标准，得出以月为单位的生态流量过程。为确定鱼类生态流量过程，在一个水文年内分时段优先考虑鱼类的不同生长阶段的生境需求。每年的4月到6月优先考虑满足鱼类产卵场要求，7月至10月优先满足鱼类捕食场要求，11月至次年3优先考虑鱼类越冬场要求(表1)。
表1生态流量计算原则
(5-2)河流生态流量过程计算：首先计算参照体系下的有效栖息地面积，使用有效栖息地面积恢复比例为指标，确定鱼类栖息地恢复的目标，其次采用建立的栖息地模型反算不同恢复目标对应的流量(可能非单一)，并进行栖息地(水力栖息地适应性指数、连通指数)评价，从中选择合适的流量值，最后根据时间尺度和鱼类生活史，确定完整水文年内的流量过程。
确定河道鱼类栖息地保护的参照体系，计算参照条件下的有效栖息地面积。根据研究区域不同的栖息地恢复目标，分别计算栖息地为最优情况的50～60％，60～80％及80～100％三种情况的生态流量区间，最终得出研究河段生态流量过程(图6)。
以上所揭露的仅为本发明一种较佳实施例而已，不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。
Claims (6)
1.一种面向鱼类生境保护的河流生态流量过程推求方法，其特征在于该方法包括：
(1)对研究区域鱼类资源现状进行调查，确定生态流量推求的目标鱼类；
(2)筛选关键水环境因子，并利用胁迫实验建立目标鱼类不同生命阶段与关键水环境因子间的响应关系；
(3)根据研究区域的水文、水质、地形及气象基础数据建立水动力水环境模型，并结合目标鱼类对关键生境因子的定量响应关系曲线，采用模糊数学原理建立鱼类栖息地模型；
(4)建立鱼类栖息地评价体系，对鱼类栖息地模型的模拟结果进行评价，并构建评价指标与流量的动态响应关系；
(5)确定平原河流生态流量计算的参照体系，并依照该参照体系，综合考虑鱼类不同生命阶段的生境需求及天然径流的季节变化，推求河流完整水文年生态流量过程。
2.根据权利要求1所述的方法，其特征在于：所述步骤(1)具体包括：
(1-1)鱼类资源现状调查：依据研究区域鱼类资源现状及历史资料，确定河流生态流量推求备选鱼类，其中，鱼类资源现状调查包括调查鱼类的种群数量，测定鱼类群体的结构状况，掌握渔业利用情况，研究鱼类种群和数量变动规律；
(1-2)目标鱼类筛选：分析鱼类生物学特性，筛选对水环境要求具有代表性的鱼类作为目标鱼类。
3.根据权利要求1所述的方法，其特征在于：所述步骤(2)具体包括：
(2-1)关键水环境因子筛选：针对筛选出的目标鱼类，根据鱼类生物学特性选择流速、水深、温度及溶氧为目标鱼类的关键水环境因子；
(2-2)响应关系曲线建立：将鱼类生命阶段分为为产卵期、捕食期及越冬期，利用胁迫实验建立目标鱼类不同生命阶段对关键水环境因子的定量响应关系。
4.根据权利要求1所述的方法，其特征在于：所述步骤(3)具体包括：
(3-1)水动力水环境模型建立：收集研究区域的水文、水质、地形及气象基础数据，根据所述基础数据采用开源模型SELFE建立水动力水环境模型，其中，建立的水动力水环境模型的控制方程为：
式中，Qa为流量，单位m3/s；t表示时间，单位为s；H为水位，单位m；h为水深，单位m，u,v为x,y方向上的速度，单位m/s；p为大气压强，单位N/m2；ρ0为水体密度，值为1025kg/m3；ν为水平粘性系数，单位m2/s；f为科氏力系数；τx,τy为底部剪应力，单位N/m；c为物质浓度，单位kg/m3；Dx,Dy为x,y方向扩散系数，单位m2/s；S为源汇项；fR(c,t)为反应项；
(3-2)鱼类栖息地模型建立：将目标鱼类不同生命阶段对关键水环境因子的响应关系和水动力水环境模型做为输入参数，采用模糊数学原理建立鱼类栖息地模型。
5.根据权利要求1所述的方法，其特征在于：所述步骤(4)具体包括：
(4-1)鱼类栖息地评价体系建立：采用水力栖息地适应性指数IHHS和栖息地连通性指数HCI作为鱼类栖息地的评价指标；其中：
水力栖息地适应性指数IHHS的计算公式为：
式中，Ai表示对研究区域划分出的网格计算单元i的面积，单位为m2；HSIi为网格计算单元i的栖息地适应性指数，N为网格计算单元数量。
栖息地连通指数HCI的计算公式为：
其式中，M为连通性系数；MDPi为斑块i与其他斑块间最短距离，单位为m；，定义为适应性大于一定阈值且面积大于一定阈值的栖息地斑块为有效斑块pe，单位m2，ne为有效斑块个数，N为网格计算单元数量；
(4-2)将流量设置为不同值，重复步骤(1)-(4)，从而得到目标鱼类栖息地评价指标与流量的动态响应关系曲线。
6.根据权利要求5所述的方法，其特征在于：所述步骤(5)具体包括：
(5-1)参照体系及鱼类生境需求确定：构建流量与IHHS和HCI的三维曲面，并选择三维曲面中的最优点相应的栖息地面积作为参照体系，并根据鱼类生物学特性将目标鱼类的产卵期、捕食期及越冬期与水文年的月份对应，且在一个水文年份中，产卵期优先考虑产卵场的需求，其次考虑索饵场；越冬期同时考虑越冬场及索饵场需求；其他时段考率满足鱼类基本生存条件的需求；
(5-2)河流生态流量过程计算：首先计算参照体系下的有效栖息地面积，使用有效栖息地面积恢复比例为指标，确定鱼类栖息地恢复的目标，其次采用建立的栖息地模型反算不同恢复目标对应的流量，并进行栖息地评价，从中选择合适的流量值，最后根据时间尺度和鱼类生活史，确定完整水文年内的流量过程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811523739.5A CN109615076A (zh) | 2018-12-13 | 2018-12-13 | 一种面向鱼类生境保护的河流生态流量过程推求方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811523739.5A CN109615076A (zh) | 2018-12-13 | 2018-12-13 | 一种面向鱼类生境保护的河流生态流量过程推求方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109615076A true CN109615076A (zh) | 2019-04-12 |
Family
ID=66008094
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811523739.5A Pending CN109615076A (zh) | 2018-12-13 | 2018-12-13 | 一种面向鱼类生境保护的河流生态流量过程推求方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109615076A (zh) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111241469A (zh) * | 2020-01-19 | 2020-06-05 | 华北电力大学 | 一种构造河流生态流量隶属度函数的方法 |
| CN111477266A (zh) * | 2020-03-04 | 2020-07-31 | 云南大学 | 一种天然水域鱼类产卵场判识方法 |
| CN112215389A (zh) * | 2019-07-10 | 2021-01-12 | 中国科学院地理科学与资源研究所 | 一种确定河流环境流量过程区间的方法 |
| CN112487640A (zh) * | 2020-11-27 | 2021-03-12 | 交通运输部天津水运工程科学研究所 | 一种内河航道整治工程生态影响模拟预测方法 |
| CN112766594A (zh) * | 2021-01-28 | 2021-05-07 | 中国电建集团北京勘测设计研究院有限公司 | 一种基于生态习性的河流鱼类栖息地适宜面积的计算方法 |
| CN113176328A (zh) * | 2021-04-25 | 2021-07-27 | 中国水产科学研究院珠江水产研究所 | 一种基于耳石微区原位微量元素评估水利枢纽对鱼类生境影响的方法 |
| CN113283743A (zh) * | 2021-05-21 | 2021-08-20 | 中国科学院南京地理与湖泊研究所 | 一种流域中不同生态修复类型生境阈值的判定方法 |
| CN113475431A (zh) * | 2021-07-16 | 2021-10-08 | 四川大学 | 构建与鱼类生命史匹配的支流生境替代适宜性模型的方法 |
| CN113688497A (zh) * | 2021-07-05 | 2021-11-23 | 郑州大学 | 基于生态-生境-流量分组响应机制的生态流量量化方法 |
| CN113780749A (zh) * | 2021-08-16 | 2021-12-10 | 烟台大学 | 一种城市河道鱼类生境评价方法和系统 |
| CN113807017A (zh) * | 2021-09-23 | 2021-12-17 | 西安理工大学 | 一种鱼类偏好栖息地的确定方法及终端设备 |
| CN113837534A (zh) * | 2021-08-17 | 2021-12-24 | 中国海洋大学 | 一种基于栖息地特征的鱼类资源调查与评估方法 |
| CN114722352A (zh) * | 2022-05-05 | 2022-07-08 | 北京师范大学 | 滨海湿地垂向水文连通指数构建方法 |
| CN116432902A (zh) * | 2023-03-31 | 2023-07-14 | 中国水利水电科学研究院 | 一种考虑水环境因子变异的物种长期生存能力评估方法 |
| CN117449243A (zh) * | 2023-11-08 | 2024-01-26 | 中国水利水电科学研究院 | 一种用于栖息地修复的生态丁坝群构建方法及系统 |
| CN117973706A (zh) * | 2024-04-02 | 2024-05-03 | 水利部交通运输部国家能源局南京水利科学研究院 | 融合鱼类生态流量过程水温过程需求的水库生态调度方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107464202A (zh) * | 2017-08-14 | 2017-12-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于鱼类生境需求的多闸坝平原河流生态流量过程推求方法 |
| CN107480884A (zh) * | 2017-08-14 | 2017-12-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种多闸坝平原河流冬季生态流量推求方法 |
-
2018
- 2018-12-13 CN CN201811523739.5A patent/CN109615076A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107464202A (zh) * | 2017-08-14 | 2017-12-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于鱼类生境需求的多闸坝平原河流生态流量过程推求方法 |
| CN107480884A (zh) * | 2017-08-14 | 2017-12-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种多闸坝平原河流冬季生态流量推求方法 |
Non-Patent Citations (4)
| Title |
|---|
| 李若男: "鱼类生境模型及河流生态流量研究", 《万方学位论文数据库》 * |
| 李若男等: "模糊数学方法模拟水库运行影响下鱼类栖息地的变化", 《生态学报》 * |
| 王丽等: "淮河干流基于生态流量的目标鱼类选择研究", 《环境科学学报》 * |
| 陈求稳: "生态水力学及其在水利工程生态环境效应模拟调控中的应用", 《水利学报》 * |
Cited By (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112215389A (zh) * | 2019-07-10 | 2021-01-12 | 中国科学院地理科学与资源研究所 | 一种确定河流环境流量过程区间的方法 |
| CN112215389B (zh) * | 2019-07-10 | 2024-03-08 | 中国科学院地理科学与资源研究所 | 一种确定河流环境流量过程区间的方法 |
| CN111241469A (zh) * | 2020-01-19 | 2020-06-05 | 华北电力大学 | 一种构造河流生态流量隶属度函数的方法 |
| CN111477266A (zh) * | 2020-03-04 | 2020-07-31 | 云南大学 | 一种天然水域鱼类产卵场判识方法 |
| CN111477266B (zh) * | 2020-03-04 | 2021-03-12 | 云南大学 | 一种天然水域鱼类产卵场判识方法 |
| CN112487640B (zh) * | 2020-11-27 | 2023-02-14 | 交通运输部天津水运工程科学研究所 | 一种内河航道整治工程生态影响模拟预测方法 |
| CN112487640A (zh) * | 2020-11-27 | 2021-03-12 | 交通运输部天津水运工程科学研究所 | 一种内河航道整治工程生态影响模拟预测方法 |
| CN112766594A (zh) * | 2021-01-28 | 2021-05-07 | 中国电建集团北京勘测设计研究院有限公司 | 一种基于生态习性的河流鱼类栖息地适宜面积的计算方法 |
| CN113176328A (zh) * | 2021-04-25 | 2021-07-27 | 中国水产科学研究院珠江水产研究所 | 一种基于耳石微区原位微量元素评估水利枢纽对鱼类生境影响的方法 |
| CN113283743B (zh) * | 2021-05-21 | 2023-06-20 | 中国科学院南京地理与湖泊研究所 | 一种流域中不同生态修复类型生境阈值的判定方法 |
| CN113283743A (zh) * | 2021-05-21 | 2021-08-20 | 中国科学院南京地理与湖泊研究所 | 一种流域中不同生态修复类型生境阈值的判定方法 |
| CN113688497A (zh) * | 2021-07-05 | 2021-11-23 | 郑州大学 | 基于生态-生境-流量分组响应机制的生态流量量化方法 |
| CN113688497B (zh) * | 2021-07-05 | 2024-02-02 | 郑州大学 | 基于生态-生境-流量分组响应机制的生态流量量化方法 |
| CN113475431B (zh) * | 2021-07-16 | 2022-05-10 | 四川大学 | 构建与鱼类生命史匹配的支流生境替代适宜性模型的方法 |
| CN113475431A (zh) * | 2021-07-16 | 2021-10-08 | 四川大学 | 构建与鱼类生命史匹配的支流生境替代适宜性模型的方法 |
| CN113780749B (zh) * | 2021-08-16 | 2024-07-05 | 烟台大学 | 一种城市河道鱼类生境评价方法和系统 |
| CN113780749A (zh) * | 2021-08-16 | 2021-12-10 | 烟台大学 | 一种城市河道鱼类生境评价方法和系统 |
| CN113837534B (zh) * | 2021-08-17 | 2023-11-24 | 中国海洋大学 | 一种基于栖息地特征的鱼类资源调查与评估方法 |
| CN113837534A (zh) * | 2021-08-17 | 2021-12-24 | 中国海洋大学 | 一种基于栖息地特征的鱼类资源调查与评估方法 |
| CN113807017B (zh) * | 2021-09-23 | 2022-05-24 | 西安理工大学 | 一种鱼类偏好栖息地的确定方法及终端设备 |
| US11645354B2 (en) | 2021-09-23 | 2023-05-09 | Xi'an University Of Technology | Determination method for preferred habitat of fish and terminal device |
| CN113807017A (zh) * | 2021-09-23 | 2021-12-17 | 西安理工大学 | 一种鱼类偏好栖息地的确定方法及终端设备 |
| CN114722352A (zh) * | 2022-05-05 | 2022-07-08 | 北京师范大学 | 滨海湿地垂向水文连通指数构建方法 |
| CN116432902B (zh) * | 2023-03-31 | 2023-09-12 | 中国水利水电科学研究院 | 一种考虑水环境因子变异的物种长期生存能力评估方法 |
| CN116432902A (zh) * | 2023-03-31 | 2023-07-14 | 中国水利水电科学研究院 | 一种考虑水环境因子变异的物种长期生存能力评估方法 |
| CN117449243A (zh) * | 2023-11-08 | 2024-01-26 | 中国水利水电科学研究院 | 一种用于栖息地修复的生态丁坝群构建方法及系统 |
| CN117973706A (zh) * | 2024-04-02 | 2024-05-03 | 水利部交通运输部国家能源局南京水利科学研究院 | 融合鱼类生态流量过程水温过程需求的水库生态调度方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109615076A (zh) | 一种面向鱼类生境保护的河流生态流量过程推求方法 | |
| CN107464202B (zh) | 基于鱼类生境需求的多闸坝平原河流生态流量过程推求方法 | |
| CN107563610B (zh) | 一种闸坝调控对鱼类栖息地空间特性影响的量化分析方法 | |
| CN107480884B (zh) | 一种多闸坝平原河流冬季生态流量推求方法 | |
| Wang et al. | Modeling oyster growth rate by coupling oyster population and hydrodynamic models for Apalachicola Bay, Florida, USA | |
| Lamberts et al. | Flood pulse alterations and productivity of the Tonle Sap ecosystem: a model for impact assessment | |
| CN106295833B (zh) | 一种太平洋褶柔鱼资源补充量预测方法及其应用 | |
| EP1907321B1 (en) | Method of sequestering carbon dioxide in aqueous environments | |
| Sánchez-Garrido et al. | Decadal-scale variability of sardine and anchovy simulated with an end-to-end coupled model of the Canary Current ecosystem | |
| Sipkay et al. | Trends in research on the possible effects of climate change concerning aquatic ecosystems with special emphasis on the modelling approach | |
| White et al. | Distribution of Nephrops norvegicus (L.) larvae in the western Irish Sea: an example of advective control on recruitment | |
| Izmaniar et al. | The Business Prospect of Climbing Perch Fish Farming with Biofloc Technology at De’Papuyu Farm Banjarbaru | |
| Flynn et al. | Modeling the lateral variation of bottom-attached algae in rivers | |
| Aich et al. | The development of a geographic information system (GIS) to document research in an Everglades physical model | |
| Newell et al. | Farm-scale production models | |
| Mandal et al. | How newly recruited cohorts are formed in the trochid gastropod population (Umbonium moniliferum) on an intertidal sandflat in western Kyushu, Japan | |
| Srisunont et al. | Development of models for sustainable green mussel cultivation under climate change events | |
| Stepputtis | Distribution patterns of Baltic sprat (Sprattus sprattus L.): causes and consequences | |
| Konrad | Approaches for evaluating the effects of bivalve filter feeding on nutrient dynamics in Puget Sound, Washington | |
| Katrasov et al. | Estimation of bivalve mollusk plantation productivity based on simulation results | |
| Hasegawa et al. | Association between the planktonic larval and benthic stages of Manila clam Ruditapes philippinarum in eastern Hokkaido, Japan | |
| Geffen et al. | Lemon sole Microstomus kitt in the northern North Sea: a multidisciplinary approach to the early life‐history dynamics | |
| Junaidi et al. | Identifying the Environmental Factors Affecting Puerulus Settlement of the Spiny Lobster, Panulirus homarus in Awang Bay, Lombok Island | |
| Capinpin Jr | Status of abalone fishery and experiential mariculture as a resource conservation strategy in Carot, Anda, Pangasinan. | |
| CN114223588B (zh) | 一种保障库区粘草性鱼卵孵化的水库生态调度方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Chen Qiuwen Inventor after: Li Ruonan Inventor after: Zhang Jianyun Inventor after: Wang Li Inventor after: Mo Kangle Inventor after: Lin Yuqing Inventor after: Wang Jun Inventor before: Chen Qiuwen Inventor before: Li Ruonan Inventor before: Wang Li Inventor before: Mo Kangle Inventor before: Lin Yuqing Inventor before: Wang Jun |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190412 |