CN111651895B - 一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法 - Google Patents
一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法 Download PDFInfo
- Publication number
- CN111651895B CN111651895B CN202010541227.2A CN202010541227A CN111651895B CN 111651895 B CN111651895 B CN 111651895B CN 202010541227 A CN202010541227 A CN 202010541227A CN 111651895 B CN111651895 B CN 111651895B
- Authority
- CN
- China
- Prior art keywords
- river
- habitat
- model
- local
- dam
- 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.)
- Active
Links
- 238000010276 construction Methods 0.000 title abstract description 12
- 241000251468 Actinopterygii Species 0.000 claims abstract description 46
- 230000006872 improvement Effects 0.000 claims abstract description 28
- 238000004088 simulation Methods 0.000 claims abstract description 19
- 238000012876 topography Methods 0.000 claims abstract description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 78
- 239000013049 sediment Substances 0.000 claims description 53
- 230000008859 change Effects 0.000 claims description 49
- 238000000034 method Methods 0.000 claims description 31
- 239000004576 sand Substances 0.000 claims description 17
- 238000011144 upstream manufacturing Methods 0.000 claims description 17
- 239000002245 particle Substances 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000003628 erosive effect Effects 0.000 claims description 12
- 239000000758 substrate Substances 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 9
- 238000011010 flushing procedure Methods 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 8
- 238000005859 coupling reaction Methods 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000007613 environmental effect Effects 0.000 claims description 6
- 230000006399 behavior Effects 0.000 claims description 5
- 230000006870 function Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000004062 sedimentation Methods 0.000 claims description 4
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims description 2
- 229920006395 saturated elastomer Polymers 0.000 claims description 2
- 238000013519 translation Methods 0.000 claims description 2
- 239000007787 solid Substances 0.000 claims 1
- 238000013461 design Methods 0.000 description 13
- 101100041681 Takifugu rubripes sand gene Proteins 0.000 description 12
- 239000002023 wood Substances 0.000 description 11
- 241000252185 Cobitidae Species 0.000 description 9
- 125000000484 butyl group Chemical group [H]C([*])([H])C([H])([H])C([H])([H])C([H])([H])[H] 0.000 description 9
- 230000001965 increasing effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 239000010866 blackwater Substances 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 241000894007 species Species 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000001617 migratory effect Effects 0.000 description 3
- 238000009991 scouring Methods 0.000 description 3
- 239000010802 sludge Substances 0.000 description 3
- 241000124879 Grus leucogeranus Species 0.000 description 2
- 238000005336 cracking Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000004083 survival effect Effects 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 241000525154 Micrathena beta Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 235000013601 eggs Nutrition 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000037406 food intake Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000017448 oviposition Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000035755 proliferation Effects 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009418 renovation Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B3/00—Engineering works in connection with control or use of streams, rivers, coasts, or other marine sites; Sealings or joints for engineering works in general
- E02B3/02—Stream regulation, e.g. breaking up subaqueous rock, cleaning the beds of waterways, directing the water flow
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/80—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in fisheries management
- Y02A40/81—Aquaculture, e.g. of fish
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Structural Engineering (AREA)
- Evolutionary Computation (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Ocean & Marine Engineering (AREA)
- Mechanical Engineering (AREA)
- Revetment (AREA)
Abstract
本发明公开了一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法,包括:(1)从目标干流的支流中,选择与目标干流鱼类群落结构特征和生境异质性最相似的支流,作为替代生境支流;(2)构建替代生境支流的一维生态河貌模型,模拟拆坝后替代生境支流的河貌变化和生境质量变化;(3)针对替代生境支流拆坝后生境质量低于预设阈值的局部河段,构建二维生态河貌演变模型,模拟实施不同局部生态改善措施后局部河段的河貌变化及生境质量变化;(4)对选择作为替代生境的支流的水坝进行拆除,并对局部河段按照模拟得到的生境质量改善最大的生态改善措施,对局部河段地貌进行改造。本发明可以构造更适宜鱼类的支流替代生境。
Description
技术领域
本发明涉及水利与环境工程,尤其涉及一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法。
背景技术
水利工程作为人类改造自然、利用自然的重要手段,在防洪、发电、航运以及流域调水等方面为人类社会进步做出了重要的贡献。然而,越来越多拦河坝的兴建也给区域生态环境带来了不同程度的影响,其中又以梯级水电站的建设对生态环境影响最为剧烈。随着全民环保意识的提升,水电开发过程中的生态环境保护问题已经逐渐成为人们关注的焦点,尤其是围绕着水电开发下的鱼类保护问题。鱼类作为水生生态系统中最高级的动物,因其重要性以及对河流系统变化反应的敏感性逐渐成为河流生态系统评价中的关键指示物种之一。大坝的建设阻断了洄游性鱼类通道,致使洄游性鱼类的繁殖规模及数量较建坝前显著下降,同时还影响上下游鱼类个体行为特性及种群间的基因交流,进而导致生物遗传多样性丧失。
为了减缓水电工程对鱼类等重要水生生物资源的影响,目前,人们实施的多种鱼类保护措施已经在不同程度上起到了积极的作用,包括修建过鱼设施、设立人工增殖放流基地、水库多目标生态调度、栖息地综合保护等。受到现阶段科学认识程度、技术水平、管理难度等因素的制约,这些工程措施还难以从根本上解决鱼类因生境破碎、破坏而面临的多样性丧失的问题。近些年来,支流生境替代保护作为水电开发中生态环境保护的一种新思路逐渐受到人们的高度关注,不过目前研究较少,且不能在实施前估量实施后的效果,导致实施困难。
发明内容
发明目的:本发明针对现有技术存在的问题,提供一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法。
技术方案:本发明所述的基于拆坝及局部微地貌人工干预的支流替代生境构造方法包括:
(1)从目标干流的支流中,选择与目标干流鱼类群落结构特征和生境异质性最相似的支流,作为替代生境支流;
(2)构建替代生境支流的一维生态河貌模型,模拟拆坝后替代生境支流的河貌变化和生境质量变化,所述一维生态河貌模型通过耦合一维水动力模型、泥沙输移及河貌变化模型和鱼类栖息地模型得到;
(3)针对替代生境支流拆坝后生境质量低于预设阈值的局部河段,构建二维生态河貌演变模型,模拟实施不同局部生态改善措施后局部河段的河貌变化及生境质量变化;其中,所述二维生态河貌演变模型利用一维生态河貌模型的计算结果作为边界条件,通过耦合二维水动力模型、二维泥沙输移模型与鱼类栖息地模型得到;所述局部生态改善措施具体为:可以调整河床冲淤演变格局、改变局部流场和河道坡降、从而达到增强水域栖息地功能的水流控制结构措施;
(4)对选择的替代生境支流进行拆坝,并对局部河段按照模拟得到的生境质量改善最大的生态改善措施,对局部河段地貌进行改造。
进一步的,步骤(2)具体包括:
(2-1)根据拆坝后替代生境支流当前时刻的河貌,采用一维水动力模型计算替代生境支流当前时刻的水深和流速;
(2-2)根据当前时刻的水深和流速,采用泥沙输移及河貌变化模型计算替代生境支流当前时刻的底质组成和下一时刻的河床糙率和河道地形,并将河床糙率和河道地形作为下一时刻的河貌传递给一维水动力模型;
(2-3)根据当前时刻的水深、流速和底质组成,采用鱼类栖息地模型计算得到替代生境支流当前时刻的生境质量;
(2-4)重复步骤(2-1)到(2-3),从而得到替代生境支流全局河貌变化和生境质量变化过程。
进一步的,步骤(3)中对实施每一局部生态改善措施后局部河段的河貌变化及生境质量变化的模拟方法包括:
(3-1)根据生态改善措施实施后局部河段当前时刻的河貌,采用二维水动力模型计算局部河段当前时刻的水动力条件;
(3-2)根据当前时刻的水动力条件,采用二维泥沙输移模型计算局部河段下一时刻的河床高程、河床糙率和河道地形,并将河床高程、河床糙率和河道地形作为下一时刻的河貌传递给二维水动力模型;
(3-3)根据当前时刻的水动力条件和河貌,采用鱼类栖息地模型计算得到局部河段当前时刻的生境质量;
(3-4)重复步骤(3-1)到(3-3),从而得到实施当前局部生态改善措施后,局部河段的河貌变化及生境质量变化。
进一步的,所述局部生态改善措施包括在局部河段构建丁坝和构建深潭-浅滩结构。
进一步的,所述鱼类栖息地模型具体为适宜度曲线法得到的栖息地模型。
有益效果:本发明与现有技术相比,其显著优点是:本发明可以在实施前模拟得到实施后的效果,且还对局部生境质量不理想河段进行了改善,使得替代支流的生境更适宜鱼类生存。
附图说明
图1是本发明实施例的流程示意图;
图2是中华金沙鳅产卵场生境因子适宜度曲线(A:水深;B:流速;C:底质类型);
图3是一维生态河貌模型基本框架图;
图4是交替方向差分离散格式示意图;
图5是拆坝后河床冲淤计算的边界条件设置(上图:水流边界;下图:悬移质泥沙边界,Year1~Year10分别表示拆坝后不同年份);
图6是老木河水坝拆除后原坝址上下游河道年际冲淤分布预测(图中虚线表示每一年河床冲淤的参考面;虚线上方的实线表示河床淤积;虚线下方的实线表示冲刷);
图7老木河水坝拆除10年以后的河床累积变化预测;
图8是老木河水坝拆除后不同年份河床地貌格局(深潭、浅滩、平滩河段)分布预测;
图9是老木河水坝拆除后中华金沙鳅产卵场生境适宜度分布;
图10不同设计流量下河道水动力计算结果和金沙鳅产卵的水动力适宜度模拟结果,((a)为不同设计流量下的流速适宜度模拟结果,(b)为不同设计流量下的水深适宜度模拟结果,(c)为不同设计流量下的金沙鳅产卵的HSI模拟结果);
图11是配置丁坝后不同年份的冲淤分布和河床底质中值粒径的变化趋势((a)为配置丁坝后不同年份的冲淤分布,(b)为配置丁坝后不同年份的河床底质中值粒径的变化趋势)。
具体实施方式
本实施例提供了一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法,如图1所示,包括以下步骤:
步骤一、从目标干流的支流中,选择与目标干流鱼类群落结构特征和生境异质性最相似的支流,作为替代生境支流。
步骤二、构建替代生境支流的一维生态河貌模型,模拟拆坝后替代生境支流的河貌变化和生境质量变化。
其中,构建的一维生态河貌模型包括一维水动力模型、泥沙输移及河貌变化模型和鱼类栖息地模型,各模型具体介绍如下:
(1)一维水动力模型
一维水力模型的构建基于下述假设:①水体是不可压缩的、均匀的,即可以忽略密度的变化;②水流在垂直方向上加速度忽略不计,并且沿着垂直水流方向的压力变化可用水静力学方法计算。一维水动力数学模型的控制方程包括水流的连续性方程及动量方程,具体如下:
式中:x为沿着水流方向的空间坐标,m;t表示时间坐标,s;A表示过流面积,m2;Q表示流量,m3/s;q表示单位长度的流入或者分出的流量,m2/s;α表示动量校准系数;g表示重力加速度,m/s2;h是水面高程,m;C表示谢才系数,m1/2/s;R表示水力半径,m。β为修正系数,β的取值由Froude数Fr的大小来决定,当Fr≤1时,水流为缓流,可以正常求解方程;当Fr>1时,水流转化为急流,动量方程忽略对流项的影响。
(2)泥沙输移及河貌变化模型
A、泥沙输移模型为:
qt*=udGgr(u/u*)n
式中:Sk为第k组泥沙组分的断面平均含沙量,kg/m3;ρ'为泥沙的干容重,N/m3;A0k为第k组泥沙组分导致的冲淤面积,m2,Z0k为推移质和悬移质引起的总的河床变形量,m;qtk为第k组泥沙组分的单宽输沙率,kg/ms;Db和Eb分别是悬移质泥沙沉降和冲起的数量,kg/m2s。Ls为不平衡输移调整长度,为第k组泥沙组分的饱和单宽平衡输沙率,kg/ms。n是非均匀泥沙组分的总数;d为泥沙组分的中值粒径,mm;u为断面平均流速,m/s;u*表示摩阻流速,m/s;Ggr表示泥沙组分的无量纲输移参数,Ggr的确定依赖于描述水流运动的无量纲运动参数Fgr与描述泥沙运动的无量纲参数Agr以及泥沙粒径的无量纲参数Dgr:s是泥沙颗粒的相对密度,其计算公式可表示为:s=(ρs-ρ)/ρ,ρs和ρ分别是泥沙密度与水的密度,kg/m3;Cgr、n、m分别是和Dgr相关的无量纲系数。且
Dgr=d[g(s-1)/v2]1/3
在Dgr>60的情况下,其取值如下:
n=0
m=1.5
Cgr=0.025
当1<Dgr≤60,其计算的表达式如下:
n=1.0-0.56logDgr
m=9.66Dgr+1.34
logCgr=2.86logDgr-(logDgr)2-3.53
B、河貌变化模型为:
混合层厚度:
河床变形控制方程:
式中,Em为混合层厚度,混合层的上边界为床面,下边界为不受水流波及的原始河床,dl为床沙中不动颗粒的最小粒径,mm;l为对应的颗粒组序数;n为颗粒组数目,ε为泥沙颗粒的孔隙率;Pk为泥沙组分的百分比;Zk为第k组泥沙冲淤造成的河床高程变化,m。
(3)水沙模型边界处理
溢流坝的溢流流量和上下游水位之间的关系可通过堰流公式计算:
淹没出流:
非淹没出流:
式中,m为溢流坝的过流系数;b为溢流坝的宽度,Zs为坝顶高程;Z11和Z12分别是坝址上下游计算节点的水位高程,m。
(4)鱼类栖息地模型
鱼类栖息地通常指某种鱼类群体生存和繁衍的地域或环境类型。本发明选取流速,水深和底质类型(床沙中值粒径)作为决定栖息地适宜度的关键生境因子。鱼类栖息地模型采用了适宜度曲线法:
HSI=(IwIvIs)1/3
式中,HSI表示中华金沙鳅产卵场的栖息地适宜度指数,用于评价生境质量,Iw,Iv和Is分别中华金沙鳅产卵行为对水深,流速和底质的适合度,环境因子的适合度指数在0~1之间,0表示对应条件完全不适合产卵,1表示最适的产卵条件。适合度指数的取值由适宜度曲线确定(图2)。
对于已经构建好的一维生态模型,进行拆坝后替代生境支流的河貌变化和生境质量变化的模拟方法如图3所示,具体包括:
(2-1)根据拆坝后替代生境支流当前时刻的河貌,采用一维水动力模型计算替代生境支流当前时刻的水深和流速;
(2-2)根据当前时刻的水深和流速,采用泥沙输移及河貌变化模型计算替代生境支流当前时刻的底质组成和下一时刻的河床糙率和河道地形,并将河床糙率和河道地形作为下一时刻的河貌传递给一维水动力模型;
(2-3)根据当前时刻的水深、流速和底质组成,采用鱼类栖息地模型计算得到替代生境支流当前时刻的生境质量;
(2-4)重复步骤(2-1)到(2-3),从而得到替代生境支流全局河貌变化和生境质量变化过程。
步骤三、针对替代生境支流拆坝后生境质量低于预设阈值的局部河段,构建二维生态河貌演变模型,模拟实施不同局部生态改善措施后局部河段的河貌变化及生境质量变化。
其中,所述二维生态河貌演变模型利用一维生态河貌模型的计算结果作为边界条件,通过耦合二维水动力模型、二维泥沙输移模型与鱼类栖息地模型得到。具体如下:
(1)二维水动力模型
二维水动力模型是基于不可压缩流体、静水压强假定以及Boussinesq假定条件下连续性方程和Naiver-Stokes方程,采用沿水深的二维平均模式,忽略动量方程中的垂向分量,具体为:
其中:
式中,z为水位,m;h为节点水深,m;u和v分别是在x和y轴方向上沿水深积分的水平流速,m/s;vt为紊动粘性系数(m2/s);ρ0为水体密度,kg/m3;τbx和τby分别为沿着x和y方向的床底切应力,N/m2;f为柯氏力系数,1/s。Mx和My分别是x和y方向的动量源汇项,m2/s;q'是单位面积上由于降水、蒸发、排水和引水等因素引起的水量变化,m/s。
(2)二维泥沙输移模型
二维泥沙输移模型采用全沙模型进行泥沙输移以及河床演变的计算,泥沙输沙率的计算采用与一维模型相同的Ackers-White模型,这里不再赘述。地形更新是基于泥沙冲淤和床面高程变化之间的平衡:
式中,Zb为河床高程,m;qtx和qty分别是沿x方向和y方向的全沙输沙率,kg/(ms)。
由于河床变化的时间尺度通常比水动力变化的时间尺度高1~2个数量级,通过在地貌模型中引入地貌加速因子(Morphological factor,MF)来解决二者在时间尺度上的不均匀性。该方法通过将每个水动力时间步长内的河貌变化量乘以地貌加速因子,使水动力结果被线性扩展到地形更新的时间步长内:
Δtmor=MF×Δthyd
式中,Δtmor和Δthyd分别是水动力计算的时间步长和河貌变化计算的时间步长,s,MF为地貌加速因子。
(3)鱼类栖息地模型与一维中的一致。
(4)一二维模型的耦合
一、二维模型的耦合的原理是:在模型的断面连接处,根据两种模型计算的水位、流量以及输沙率相等来实现的模型的耦合,水流从一维模型流向二维模型的情况下,由一维模型求解出连接断面的物理量作为二维模型的边界条件。具体耦合过程如下:
水位连接条件:z1=z2
式中,z1和z2分别为一、维模型在断面连接处的水位值,m。
一维模型在连接断面的流量值按照连接处的水深值分配到二维模型的入流边界:
Q1=∫uξhξdξ
式中,Q1为一维模型在断面连接处的流量值,m3/s;ξ为二维模型的法向坐标;uξ和hξ分别是断面连接处沿法向分布的流速和水深,m/s,m。
在构建好二维生态河貌演变模型后,步骤三会对实施每一局部生态改善措施后局部河段的河貌变化、水动力条件变化及生境质量的改善进行模拟,所述局部生态改善措施具体为:可以调整河床冲淤演变格局、改变局部流场和河道坡降、从而达到增强水域栖息地功能的水流控制结构措施,包括在局部河段构建丁坝和构建深潭-浅滩结构等。模拟的方法具体包括:
(3-1)根据生态改善措施实施后局部河段当前时刻的河貌,采用二维水动力模型计算局部河段当前时刻的水动力条件;
(3-2)根据当前时刻的水动力条件,采用二维泥沙输移模型计算局部河段下一时刻的河床高程、河床糙率和河道地形,并将河床高程、河床糙率和河道地形作为下一时刻的河貌传递给二维水动力模型;
(3-3)根据当前时刻的水动力条件和河貌,采用鱼类栖息地模型计算得到局部河段当前时刻的生境质量;
(3-4)重复步骤(3-1)到(3-3),从而得到实施当前局部生态改善措施后,局部河段的河貌变化及生境质量变化。
其中,二维模型的数值计算格式采用了交替方向的隐式差分格式(Alternatingdirection implicit,ADI),在水平面上采用交错网格,将对原偏微分方程组的求解转化为求解网格上离散点的变量值,u、v、z、h变量在网格中指定的位置进行计算(图4)。在时间上采用了显、隐格式交替进行计算。具体方法是把每一个时间步长等分成两个时间层,在前半个时间步采用隐式差分格式求解x轴方向上动量方程和连续性方程,得到时间层面上的U值和Z值,然后在y方向上联立动量方程和连续性方程,采用显示差分格式求得/>时间层上的V值;依照前半部的离散方式,在后半个时间步对调求解顺序,首先在y方向上联解动量方程和连续性方程,用隐式差分格式求得(n+1)Δt时间层上的V值和Z值,然后采用显示差分求得(n+1)Δt时间层上的U值。下列所示为ADI差分格式:
其中
通过ADI法对连续性方程和动量方程差分求解可以得到每一个计算时间步长的水动力场分布,并用于泥沙输移场和河床地形的更新,时间步长Δt的选取内,相邻网格在同一方向上的长度比不应超过1.2。此外,时间步长Δt的选取对模型稳定性和精度都有很大影响,时间步长过大可能会使计算不容易收敛,相反时间步长选取过小则会大大增加有效计算时间且放大离散误差,时间步长Δt的确定由Courant数确定:
式中,Cf为Courant数,当其小于时认为时间步长的选取满足ADI算法精度要求;Δx和Δy分别是在x和y轴方向的最小网格尺度,m。
步骤四、对选择的替代生境支流进行拆坝,并对局部河段按照模拟得到的生境质量改善最大的生态改善措施,对局部河段地貌进行改造。
下面采用具体例子来说明。
该例中,选择中国金沙江下游作为目标干流,金沙江下游河流生境具有较高的空间异质性和时间动态变化特征,为不同鱼类的生境选择提供了适宜的场所。在金沙江下游“白鹤滩~乌东德”段各主要支流中,黑水河植被资源变化多样,河道呈“V”字型,曲折蜿蜒,形成了鱼类生境的多样性,白鹤滩水库建成蓄水后,原来急流的环境会逐渐过度到河道型水库环境,急流产卵的土著鱼类将会失去赖以生存的产卵场生境。因此,从鱼类群落结构特征及生境异质性来看,黑水河鱼类群落结构特征和生境异质性与金沙江干流最为相似,选择黑水河作为替代生境支流。选择中华金沙鳅作为目标鱼类。老木河水坝是黑水河下游四级电站中的最下一级,坝址处距离河口41.2km,建于1988年,基本无调蓄能力,因此对老木河水坝拆坝进行模拟。
采用一维生态河貌模型进行模拟,预测原坝址上下游河道在未来10年内的河貌变化以及生境质量(栖息地适宜度)变化情况,其中,第一年模拟的边界条件选用2018年12月~2019年12月的实测水沙序列,其余9年的边界条件选用2009年12月~2018年12月水文系列年的流量和沙量序列(图5)。河床冲淤变化模拟结果如图6所示,可以看出拆坝后河床冲淤变化分布具有显著的差异。在拆坝后第1年(Year1),由于库区淤沙的无控释放,上游河床迅速侵蚀下切,最大冲刷深度出现在坝体附近,并且随着时间推移,侵蚀裂点由坝体逐渐向上游发展;坝址下游河床受冲刷下泄泥沙的影响而出现显著抬高,淤积现象主要集中在坝址下游1.0km以内的河道,而且淤积深度沿水流方向降低。随着库区淤沙的持续冲刷,地形动力结构处于不断调整之中,上游河段出现了溯源侵蚀现象,并且侵蚀裂点向更上游发展的趋势,由此冲刷的沉积物在下游河床以沙垄和沙浪的形式向前推进,并不断拉伸加长。拆坝5年后(Year5),河床形态调整幅度减小,河道地形基本达到稳定。但是至拆坝后第7年(Year7),床面再一次出现了显著的起伏变化,原坝址上游库区河段发生强烈冲刷侵蚀,最大下切深度达到了2.26m,淤积现象主要出现在坝址下游12~13km处。从老木河水坝拆除10年后的河床累积变化来看(图7),坝址上游河床变化主要以侵蚀下切为主,侵蚀现象主要集中在原库区段,库区上游河床未受到较大扰动。长远来看,拆坝后库区河床比降由拆坝前2.2‰增加至3.8‰,坝址上下游河床恢复自然连通,河流形态回归至建坝前的原始状态,提高洄游鱼类和溯河产卵鱼类种群的物种数量;其次,大坝的修建使原有连续的河流生态系统被分割成不联系的片段,导致群体之间出现遗传分化。此外,上游冲刷的沉积物主要淤积在下游3.0~8.0km的河道内,主要是由于河道断面拓宽、水流能量梯度降低导致输沙能力减弱。老木河水坝拆除后不同年份的河床地貌模拟结果表明(图8),拆坝工程显著提高了研究河段的地貌格局多样性,部分河段出现了浅滩地貌,地形地貌的显著变化提高了水生生物栖息地多样性。
拆坝后的栖息地适宜度变化情况如图9所示,该模拟在中华金沙鳅产卵期50%的月平均流量下,基于拆坝后的地形条件对该物种产卵场的生境适宜度进行。由图可以看出,在老木河水坝拆除以前,只有16%的河段(3.0km)产卵条件较为适宜(HSI>0.8),同时有接近42%的河段(7.5km)产卵场适宜度较差(HSI<0.4);拆坝1年后,有45%的河段(8.1km)适合作为产卵场(HSI>0.8),不适合产卵(HSI<0.4)的河段降至26%(4.7km),产卵场质量得到了显著提升;拆坝后第5年至第10年,适合产卵(HSI>0.8)的河段由40%(7.2km)增加至68%(12.2km),同时不适合产卵(HSI<0.4)的河段由13%(2.4km)降至9%(1.7km)。
总体而言,小型水坝导致的河流地貌改变会促进河道水流多样性,进而提高鱼类产卵场的适宜度指数。但是研究区域最下游5km河段生境质量并没有得到显著性改善,其主要原因是上游建坝以后,下游河道因遭受清水冲刷而出现了床质粗化现象,导致原有河床结构退化,而且老木河水坝的水库库容十分有限,拆坝以后侵蚀输送下来的推移质泥沙很难运动至此。因此,该区域河床形态以及水流结构均没有太大改变,需要有针对性地采取人工修复措施(构建丁坝、深潭-浅滩结构)来改造河床结构,丰富水流多样性,从而营造水生生物适宜的栖息环境,并采用二维模型对该区域河形以及鱼类栖息地的变化做进一步研究。
二维生态河貌演变模型的模型选择老木河水坝下游11km处的顺河村至大花地村附近的河段为研究对象,对该河段设计不同设计频率流量下中华金沙鳅产卵场栖息地适合度模拟,模拟结果见图10。根据图10可知,在设计频率为75%、50%和25%的产卵期月平均流量下,该河段为中华金沙鳅的产卵行为提供了良好的水深环境,但流速适合度(Iv)整体较低,几乎所有生境单元的流速适宜度都在0.6以下;在10%的设计流量下,研究河段的流速适宜度有所提高,有接近20%的生境单元达到了较适宜的流速状态(Iv>0.8),但此时部分区域水深过大,超出了适宜范围。对于HIS,在75%的设计流量下,所有生境单元的栖息地适合度指数(HSI)均在0.8以下,表明该河段水力生境较差,不适合中华金沙鳅产卵;在50%的设计流量下,除弯道右岸较短区域外,其他区域基本不适合产卵(HSI<0.4),产卵场整体状况较差;在25%的设计流量下,有接近5%的生境单元产卵条件较为适宜(HSI>0.8),主要集中在下游的弯道处;在10%的设计流量下,产卵场适合度明显提高,有超过34%的生境单元适合产卵,同时,可以发现弯道以上河段水深超过了适宜范围,导致产卵场栖息地质量较差。综上可知,在75%~25%的设计流量时,主要受流速限制,产卵场的整体适合度较差;当流量增加至10%时,产卵场的适合度虽然有所改善,但由于水深超过了适宜范围,适宜中华金沙鳅产卵的区域面积较小,会导致产卵行为相对集中,不利于产卵亲鱼的摄食和鱼卵的发育。因此,有必要采取有针对性的人工修复措施诱导河貌演变,增加水力学条件的多样性,创造多样化的栖息地,从而维持不同流量下产卵场生境质量的稳定性,为溯河到支流产卵的鱼类提供适宜的水流环境。
本例采用丁坝(spur dike)作为局部改善措施。丁坝是从河道岸边伸出、在河流平面上与河岸构成“丁字形”的河道整治结构物,在发挥束水攻沙与河道治理功能的同时也会造成水流方向与泥沙输移的改变。本例将丁坝沿河岸两侧交叉布置形成丁坝群,坝顶顶面高于河流水位或河岸顶面,坝体长度(L)与平滩河宽之比应保持在0.25~0.33之间丁坝顶端至河对岸的距离在70%~80%范围,即坝体长度不应超过河宽的30%,因此本文设计的丁坝长度为24m,为平均河宽的30%。坝体长度(L)和相邻坝体之间的错口距离(D)比值(D/L)不小于3.4,本例选择5,即生态丁坝错口距离为120m,为坝体长度的5倍。模拟结果见图11,根据图可以看出,丁坝有效改善了河床冲淤变化和河床底质组成。综上实验结果,可以看出本发明构造的替代支流生境可以满足干流鱼类生境质量需求。
Claims (5)
1.一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法,其特征在于该方法包括:
(1)从目标干流的支流中,选择与目标干流鱼类群落结构特征和生境异质性最相似的支流,作为替代生境支流;
(2)构建替代生境支流的一维生态河貌模型,模拟拆坝后替代生境支流的河貌变化和生境质量变化,所述一维生态河貌模型通过耦合一维水动力模型、泥沙输移及河貌变化模型和鱼类栖息地模型得到,各模型具体介绍如下:
①、一维水动力模型
包括水流的连续性方程及动量方程,具体如下:
式中:x为沿着水流方向的空间坐标;t表示时间坐标;A表示过流面积;Q表示流量;q表示单位长度的流入或者分出的流量;α表示动量校准系数;g表示重力加速度;h是水面高程;C表示谢才系数;R表示水力半径,β为修正系数,β的取值由Froude数Fr的大小来决定,当Fr≤1时,水流为缓流,可以正常求解方程;当Fr>1时,水流转化为急流,动量方程忽略对流项的影响;
②泥沙输移及河貌变化模型
A、泥沙输移模型为:
qt*=udGgr(u/u*)n
式中:Sk为第k组泥沙组分的断面平均含沙量;ρ'为泥沙的干容重;A0k为第k组泥沙组分导致的冲淤面积,Z0k为推移质和悬移质引起的总的河床变形量;qtk为第k组泥沙组分的单宽输沙率;Db和Eb分别是悬移质泥沙沉降和冲起的数量,Ls为不平衡输移调整长度,为第k组泥沙组分的饱和单宽平衡输沙率,n是非均匀泥沙组分的总数;d为泥沙组分的中值粒径;u为断面平均流速;u*表示摩阻流速;Ggr表示泥沙组分的无量纲输移参数,Ggr的确定依赖于描述水流运动的无量纲运动参数Fgr与描述泥沙运动的无量纲参数Agr以及泥沙粒径的无量纲参数Dgr:且
Dgr=d[g(s-1)/v2]1/3
s是泥沙颗粒的相对密度,其计算公式可表示为:s=(ρs-ρ)/ρ,ρs和ρ分别是泥沙密度与水的密度;Cgr、n、m分别是和Dgr相关的无量纲系数,在Dgr>60的情况下,其取值如下:
n=0
m=1.5
Cgr=0.025
当1<Dgr≤60,其计算的表达式如下:
n=1.0-0.56logDgr
m=9.66Dgr+1.34
logCgr=2.86logDgr-(logDgr)2-3.53
B、河貌变化模型为:
混合层厚度:
河床变形控制方程:
式中,Em为混合层厚度,混合层的上边界为床面,下边界为不受水流波及的原始河床,dl为床沙中不动颗粒的最小粒径;l为对应的颗粒组序数;nk为颗粒组数目,ε为泥沙颗粒的孔隙率;Pk为泥沙组分的百分比;Zk为第k组泥沙冲淤造成的河床高程变化;
③水沙模型边界处理
溢流坝的溢流流量和上下游水位之间的关系可通过堰流公式计算:
淹没出流:
非淹没出流:
式中,mm为溢流坝的过流系数;b为溢流坝的宽度,Zs为坝顶高程;Z11和Z12分别是坝址上下游计算节点的水位高程;
④鱼类栖息地模型
HSI=(IwIvIs)1/3
式中,HSI表示栖息地适宜度指数,用于评价生境质量,Iw、Iv和Is分别表示产卵行为对水深、流速和底质的适合度,环境因子的适合度指数在0~1之间,0表示对应条件完全不适合产卵,1表示最适的产卵条件;
(3)针对替代生境支流拆坝后生境质量低于预设阈值的局部河段,构建二维生态河貌演变模型,模拟实施不同局部生态改善措施后局部河段的河貌变化及生境质量变化;其中,所述二维生态河貌演变模型利用一维生态河貌模型的计算结果作为边界条件,通过耦合二维水动力模型、二维泥沙输移模型与鱼类栖息地模型得到;所述局部生态改善措施具体为:调整河床冲淤演变格局、改变局部流场和河道坡降、从而达到增强水域栖息地功能的水流控制结构措施;
(4)对选择的替代生境支流进行拆坝,并对局部河段按照模拟得到的生境质量改善最大的生态改善措施,对局部河段地貌进行改造。
2.根据权利要求1所述的基于拆坝及局部微地貌人工干预的支流替代生境构造方法,其特征在于:步骤(2)具体包括:
(2-1)根据拆坝后替代生境支流当前时刻的河貌,采用一维水动力模型计算替代生境支流当前时刻的水深和流速;
(2-2)根据当前时刻的水深和流速,采用泥沙输移及河貌变化模型计算替代生境支流当前时刻的底质组成和下一时刻的河床糙率和河道地形,并将河床糙率和河道地形作为下一时刻的河貌传递给一维水动力模型;
(2-3)根据当前时刻的水深、流速和底质组成,采用鱼类栖息地模型计算得到替代生境支流当前时刻的生境质量;
(2-4)重复步骤(2-1)到(2-3),从而得到替代生境支流全局河貌变化和生境质量变化过程。
3.根据权利要求1所述的基于拆坝及局部微地貌人工干预的支流替代生境构造方法,其特征在于:步骤(3)中对实施每一局部生态改善措施后局部河段的河貌变化及生境质量变化的模拟方法为:
(3-1)根据生态改善措施实施后局部河段当前时刻的河貌,采用二维水动力模型计算局部河段当前时刻的水动力条件;
(3-2)根据当前时刻的水动力条件,采用二维泥沙输移模型计算局部河段下一时刻的河床高程、河床糙率和河道地形,并将河床高程、河床糙率和河道地形作为下一时刻的河貌传递给二维水动力模型;
(3-3)根据当前时刻的水动力条件和河貌,采用鱼类栖息地模型计算得到局部河段当前时刻的生境质量;
(3-4)重复步骤(3-1)到(3-3),从而得到实施当前局部生态改善措施后,局部河段的河貌变化及生境质量变化。
4.根据权利要求1所述的基于拆坝及局部微地貌人工干预的支流替代生境构造方法,其特征在于:所述局部生态改善措施包括在局部河段构建丁坝和构建深潭-浅滩结构。
5.根据权利要求1所述的基于拆坝及局部微地貌人工干预的支流替代生境构造方法,其特征在于:所述鱼类栖息地模型具体为适宜度曲线法得到的栖息地模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010541227.2A CN111651895B (zh) | 2020-06-15 | 2020-06-15 | 一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010541227.2A CN111651895B (zh) | 2020-06-15 | 2020-06-15 | 一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111651895A CN111651895A (zh) | 2020-09-11 |
CN111651895B true CN111651895B (zh) | 2023-11-28 |
Family
ID=72347861
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010541227.2A Active CN111651895B (zh) | 2020-06-15 | 2020-06-15 | 一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111651895B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113475431B (zh) * | 2021-07-16 | 2022-05-10 | 四川大学 | 构建与鱼类生命史匹配的支流生境替代适宜性模型的方法 |
CN113886922B (zh) * | 2021-10-08 | 2022-04-15 | 河北省水利水电勘测设计研究院集团有限公司 | 一种基于大数据和现场试验为基础的丁坝群设计方法 |
CN114707371B (zh) * | 2021-12-30 | 2023-02-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种干湿交替带氮迁移转化耦合模拟方法及系统 |
CN117807752B (zh) * | 2023-07-11 | 2024-08-23 | 江苏省水利科学研究院 | 一种沙坝潟湖型潮汐汊道中长期演变预测方法 |
CN117808648A (zh) * | 2024-01-08 | 2024-04-02 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种基于多闸控的河流农业面源污染调控方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563610A (zh) * | 2017-08-14 | 2018-01-09 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种闸坝调控对鱼类栖息地空间特性影响的量化分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140278323A1 (en) * | 2013-03-14 | 2014-09-18 | Wayne S. Wright | System and Method for Ecosystem Habitat Optimization |
-
2020
- 2020-06-15 CN CN202010541227.2A patent/CN111651895B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563610A (zh) * | 2017-08-14 | 2018-01-09 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种闸坝调控对鱼类栖息地空间特性影响的量化分析方法 |
Non-Patent Citations (2)
Title |
---|
邱阳凌."黑水河生态河貌模拟及鱼类替代生境评价".《中国优秀硕士学位论文全文数据库》.2019,第7-11、23-24、39-72页. * |
魏文杰等."丁坝布置形式对中小河流栖息地影响的数值模拟研究".《中国农村水利水电》.2017,第136-139页. * |
Also Published As
Publication number | Publication date |
---|---|
CN111651895A (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111651895B (zh) | 一种基于拆坝及局部微地貌人工干预的支流替代生境构造方法 | |
Zahar et al. | Impacts of large dams on downstream flow conditions of rivers: Aggradation and reduction of the Medjerda channel capacity downstream of the Sidi Salem dam (Tunisia) | |
CN112084671B (zh) | 城市时变增益降雨-径流过程模拟计算方法 | |
CN104695377A (zh) | 一种中小型直线化河道蜿蜒形态的构建方法 | |
CN106012947B (zh) | 一种基于地形地貌条件的鱼类栖息地划分方法 | |
CN112784502B (zh) | 一种水文水力学动态双向耦合的洪水预测方法 | |
CN111382547A (zh) | 一维非稳态河流水动力-水质模拟方法和装置 | |
CN210194532U (zh) | 用于推移质输沙率研究的植被群落河道模型 | |
Zhu et al. | Study of quality maintenance of fish habitats in small-and medium-sized mountain rivers with low flow rate | |
CN112627102B (zh) | 一种天然石块构成型多级阶梯-深潭消能系统设计方法 | |
CN113065689B (zh) | 一种多生境城市生态水系构建系统及方法 | |
CN113888710A (zh) | 一种基于dem的动力学流域产沙量模拟方法 | |
CN111783346B (zh) | 一种考虑水冰沙耦合作用的河冰运动与岸滩侵蚀计算方法 | |
CN111563353B (zh) | 一种模拟泄洪水舌对河道和堆丘形态影响的方法 | |
Yi et al. | The bedform morphology of Chinese sturgeon spawning sites in the Yangtze River | |
Im et al. | Two-dimensional physical habitat modeling of effects of habitat structures on urban stream restoration | |
CN100383343C (zh) | 环保型多用通洪、排水河床 | |
CN1641109A (zh) | 一种生态护岸 | |
CN117494477B (zh) | 一种区域磷污染负荷评估方法及系统 | |
Belli et al. | Physical Model Approach of Erosion and Sedimentation Movement in Front of Sluice Gate in the JSC Main Channel | |
CN115852910B (zh) | 一种模块化鱼坡结构、及其设计方法和施工方法 | |
Liu et al. | Can the regulation of Golden Inland Waterways meet the needs of navigation, flood control, and ecology? A model-based case study | |
CN114215000B (zh) | 一种基于水系重构的湿地生态修复方法 | |
Namu et al. | Effects of continuous flushing on the sediment removal efficiency in settling basins of small scale irrigation projects; A case study of Kiriku-Kiende irrigation project, Embu County, Kenya | |
CN221421858U (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |