• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于Erlang分布蓄水容量曲线的流域产流模型

闫宝伟 李正坤 段美壮 江慧宁 刘昱

闫宝伟, 李正坤, 段美壮, 江慧宁, 刘昱. 基于Erlang分布蓄水容量曲线的流域产流模型[J]. 水科学进展, 2021, 32(1): 120-126. doi: 10.14042/j.cnki.32.1309.2021.01.012
引用本文: 闫宝伟, 李正坤, 段美壮, 江慧宁, 刘昱. 基于Erlang分布蓄水容量曲线的流域产流模型[J]. 水科学进展, 2021, 32(1): 120-126. doi: 10.14042/j.cnki.32.1309.2021.01.012
YAN Baowei, LI Zhengkun, DUAN Meizhuang, JIANG Huining, LIU Yu. A watershed runoff yield model based on Erlang distribution water storage capacity curve[J]. Advances in Water Science, 2021, 32(1): 120-126. doi: 10.14042/j.cnki.32.1309.2021.01.012
Citation: YAN Baowei, LI Zhengkun, DUAN Meizhuang, JIANG Huining, LIU Yu. A watershed runoff yield model based on Erlang distribution water storage capacity curve[J]. Advances in Water Science, 2021, 32(1): 120-126. doi: 10.14042/j.cnki.32.1309.2021.01.012

基于Erlang分布蓄水容量曲线的流域产流模型

doi: 10.14042/j.cnki.32.1309.2021.01.012
基金项目: 

国家重点研发计划资助项目 2016YFC0402708

中央高校基本科研业务费资助项目 2017KFYXJJ195

详细信息
    作者简介:

    闫宝伟(1981-), 男, 山东滨州人, 副教授, 主要从事水文基础理论研究。E-mail:bwyan@hust.edu.cn

  • 中图分类号: TV11

A watershed runoff yield model based on Erlang distribution water storage capacity curve

Funds: 

the National Key R & D Program of China 2016YFC0402708

the Fundamental Research Funds for the Central Universities, China 2017KFYXJJ195

  • 摘要: 蓄水容量曲线是反映流域缺水量空间分布不均匀性的特征曲线,对流域产流计算有直接影响。通过对多个典型流域蓄水容量空间分布的分析,发现Erlang分布能更好地拟合流域的蓄水容量曲线,进一步基于Erlang分布进行流域产流的推导,提出了基于Erlang分布蓄水容量曲线的流域产流模型。应用结果表明,基于Erlang分布的流域产流模型,增加了模型的适应性,模拟结果更接近流域实际的产流过程,比新安江模型能取得更高的模拟精度。此外,该产流模型的参数可由流域地形和土壤类型数据估算,为无资料地区的产流计算提供了一种可行的途径。
  • 图  1  典型流域蓄水容量的空间分布

    Figure  1.  Spatial distribution of watershed storage capacity in the selected basins

    图  2  典型流域蓄水容量分布的拟合图

    Figure  2.  Fitness of distribution of the storage capacity in the selected basins

    图  3  流域蓄水容量分布曲线与降雨径流关系

    Figure  3.  Watershed storage capacity distribution curve and its rainfall runoff relationship

    图  4  栾川流域检验期洪水过程模拟结果

    Figure  4.  Simulated flood hydrographs in the validation period in Luanchuan station

    表  1  所选典型流域属性表

    Table  1.   Feature attributes of the selected typical basins

    流域名称 流域面积/km2 所在地区 所属流域 气候区
    松柏 184 湖北 长江流域 湿润区
    恒口 924 陕西 长江流域 半湿润区
    栾川 555 河南 黄河流域 半干旱区
    青狮潭 747 广西 珠江流域 湿润区
    下载: 导出CSV

    表  2  模型参数模拟结果

    Table  2.   Simulation results of model parameters

    计算模块 参数 参数含义 松柏 恒口 栾川 青狮潭
    蒸散发计算 K 蒸发折减系数 0.11 1.08 1.06 0.19
    UM 上层张力水蓄水容量/mm 10.00 13.31 11.64 12.42
    LM 下层张力水蓄水容量/mm 20.01 23.20 23.12 21.79
    DM 深层张力水蓄水容量/mm 20.00 23.40 36.34 30.41
    C 深层蒸发系数 0.05 0.09 0.05 0.08
    产流计算 n 蓄水容量曲线参数 3 2 4 9
    λ 蓄水容量曲线参数 8.35 6.15 9.17 5.99
    IM 不透水面积比例 0.01 0.01 0.03 0.01
    水源划分 SM 自由水蓄水容量/mm 29.79 33.81 32.23 34.26
    EX 自由水蓄水容量曲线方次 1.11 1.16 0.50 1.45
    KI 壤中流出流系数 0.37 0.35 0.14 0.10
    KG 地下水出流系数 0.44 0.30 0.39 0.24
    汇流计算 CS 地表水消退系数 0.80 0.80 0.72 0.80
    CI 壤中流消退系数 0.81 0.87 0.91 0.56
    CG 地下水消退系数 0.97 0.99 0.99 0.98
    下载: 导出CSV

    表  3  洪水模拟结果平均值统计及对比

    Table  3.   Mean statistics and comparison of flood simulation results

    时期 流域名称 场次 洪峰相对误差 峰现时间误差/h 洪水总量相对误差 确定性系数
    XAJ E-XAJ XAJ E-XAJ XAJ E-XAJ XAJ E-XAJ
    率定期 松柏 12 0.23 0.22 1.25 1.25 0.18 0.18 0.84 0.86
    恒口 14 0.26 0.26 1.36 1.36 0.24 0.24 0.80 0.80
    栾川 11 0.29 0.27 0.73 0.73 0.37 0.31 0.73 0.80
    青狮潭 11 0.20 0.19 0.64 0.55 0.09 0.07 0.87 0.88
    检验期 松柏 4 0.26 0.26 1.00 0.75 0.25 0.24 0.79 0.80
    恒口 6 0.26 0.26 1.36 1.36 0.24 0.24 0.80 0.80
    栾川 3 0.26 0.21 1.33 1.00 0.14 0.11 0.80 0.83
    青狮潭 3 0.13 0.09 0.67 0.67 0.16 0.15 0.89 0.89
      注:XAJ表示新安江模型;E-XAJ表示基于Erlang分布蓄水容量曲线的新安江模型。
    下载: 导出CSV
  • [1] 周买春, JAYAWARDENA W A.利用双抛物线型土壤蓄水容量曲线对新安江产流模型的改进[J].水利学报, 2002, 33(12):38-43. doi:  10.3321/j.issn:0559-9350.2002.12.008

    ZHOU M C, JAYAQARDENA W A. Modified Xin'anjiang model for runoff calculation[J]. Journal of Hydraulic Engineering, 2002, 33(12):38-43. (in Chinese) doi:  10.3321/j.issn:0559-9350.2002.12.008
    [2] GAO H K, CAI H Y, DUAN Z. Understanding the impacts of catchment characteristics on the shape of the storage capacity curve and its influence on flood flows[J]. Hydrology Research, 2017, 49(1):90-106.
    [3] 郭方, 刘新仁, 任立良.以地形为基础的流域水文模型:TOPMODEL及其拓宽应用[J].水科学进展, 2000, 11(3):298-301. http://skxjz.nhri.cn/article/id/1602

    GUO F, LIU X R, REN L L. A topography based hydrological model:TOPMODEL and its widened application[J]. Advances in Water Science, 2000, 11(3):298-301. (in Chinese) http://skxjz.nhri.cn/article/id/1602
    [4] 石朋, 芮孝芳, 瞿思敏, 等.一种通过地形指数计算流域蓄水容量的方法[J].水科学进展, 2008, 19(2):264-267. doi:  10.3321/j.issn:1001-6791.2008.02.017

    SHI P, RUI X F, QU S M, et al. Calculating storage capacity with topographic index[J]. Advances in Water Science, 2008, 19(2):264-267. (in Chinese) doi:  10.3321/j.issn:1001-6791.2008.02.017
    [5] 孔凡哲, 宋晓猛.一种新安江模型中土壤水蓄水容量空间分布的计算方法[J].水文, 2011, 31(5):1-5. https://www.cnki.com.cn/Article/CJFDTOTAL-SWZZ201105003.htm

    KONG F Z, SONG X M. A computation method for spatial distribution of soil moisture storage capacity in Xin'anjiang model[J]. Journal of China Hydrology, 2011, 31(5):1-5. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SWZZ201105003.htm
    [6] 向小华, 宋琪峰, 陈喜, 等.融合地形和土壤特征的流域蓄水容量模型[J].水科学进展, 2013, 24(5):651-657. http://skxjz.nhri.cn/article/id/2336

    XIANG X H, SONG Q F, CHEN X, et al. A storage capacity model integrating terrain and soil characteristics[J]. Advances in Water Science, 2013, 24(5):651-657. (in Chinese) http://skxjz.nhri.cn/article/id/2336
    [7] WANG D B. A new probability density function for spatial distribution of soil water storage capacity leads to the SCS curve number method[J]. Hydrology and Earth System Sciences, 2018, 22(12):6567-6578. doi:  10.5194/hess-22-6567-2018
    [8] GAO H, BIRKEL C, HRACHOWITZ M, et al. A simple topography-driven and calibration-free runoff generation module[J]. Hydrology and Earth System Sciences, 2019, 23(2):787-809. doi:  10.5194/hess-23-787-2019
    [9] TROCH P. Conceptual basin-scale runoff process models for humid catchments: analysis, synthesis and applications[D]. Ghent: Ghent University, 1992: 62-67.
    [10] GAO H, HRACHOWITZ M, SCHYMANSKI S J, et al. Climate controls how ecosystems size the root zone storage capacity at catchment scale[J]. Geophysical Research Letters, 2015, 41(22):7916-7923.
    [11] WANG-ERLANDSSON L, BASTIAANSSEN W G M, GAO H K, et al.Global root zone storage capacity from satellite-based evaporation[J].Hydrology and Earth System Sciences, 2016, 20(4):1459-1481. doi:  10.5194/hess-20-1459-2016
    [12] 曹凤云, 张建伟.机场调度的队列模型及延误研究[J].交通运输工程与信息学报, 2010, 8(4):77-80. doi:  10.3969/j.issn.1672-4747.2010.04.013

    CAO F Y, ZHANG J W. Research on the queue model and delays of airport scheduling[J]. Journal of Transportation Engineering and Information, 2010, 8(4):77-80. (in Chinese) doi:  10.3969/j.issn.1672-4747.2010.04.013
    [13] 闫宝伟, 于雨, 郭生练, 等.离散广义Nash汇流模型及其在河道洪水演算中的应用[J].水利学报, 2017, 48(7):773-778. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201707003.htm

    YAN B W, YU Y, GUO S L, et al. The discrete generalized Nash Model and its application in the river flood routing[J]. Journal of Hydraulic Engineering, 2017, 48(7):773-778. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201707003.htm
  • [1] 赵建飞, 梁忠民, 刘金涛, 李彬权, 段雅楠.  山丘区变动产流层分布式水文模型 . 水科学进展, 2022, 33(3): 429-441. doi: 10.14042/j.cnki.32.1309.2022.03.007
    [2] 童冰星, 李致家, 姚成.  基于SoilGrids的栅格新安江模型参数空间分布估算 . 水科学进展, 2022, 33(2): 219-226. doi: 10.14042/j.cnki.32.1309.2022.02.006
    [3] 田富强, 徐冉, 南熠, 李琨彪, 贺志华.  基于分布式水文模型的雅鲁藏布江径流水源组成解析 . 水科学进展, 2020, 31(3): 324-336. doi: 10.14042/j.cnki.32.1309.2020.03.002
    [4] 陈华, 霍苒, 曾强, 杨无双, 陈杰, 郭生练, 许崇育.  雨量站网布设对水文模型不确定性影响的比较 . 水科学进展, 2019, 30(1): 34-44. doi: 10.14042/j.cnki.32.1309.2019.01.004
    [5] 王卫光, 李进兴, 魏建德, 邵全喜, 邓超, 余钟波.  基于蒸散发数据同化的径流过程模拟 . 水科学进展, 2018, 29(2): 159-168. doi: 10.14042/j.cnki.32.1309.2018.02.002
    [6] 李致家, 梁珂, 阚光远, 李巧玲, 姚成, 张珂.  新安江模型中河网汇流参数Cs的一种计算方法 . 水科学进展, 2016, 27(5): 652-661. doi: 10.14042/j.cnki.32.1309.2016.05.002
    [7] 李妍敏, 安翼, 刘青泉.  植被分布对小流域产流影响的数值实验 . 水科学进展, 2015, 26(2): 187-195. doi: 10.14042/j.cnki.32.1309.2015.02.005
    [8] 洪安宇, 王正中, 赵延风.  层流区Shields曲线分布规律试验 . 水科学进展, 2015, 26(6): 837-845. doi: 10.14042/j.cnki.32.1309.2015.06.010
    [9] 李致家, 黄鹏年, 姚成, 李巧玲, 赵丽霞, 余钟波.  灵活架构水文模型在不同产流区的应用 . 水科学进展, 2014, 25(1): 28-35.
    [10] 张小琴, 刘可新, 包为民, 李佳佳, 赖善证.  产流误差比例系数的系统响应修正方法 . 水科学进展, 2014, 25(6): 789-796.
    [11] 向小华, 宋琪峰, 陈喜, 吴晓玲, 王船海.  融合地形和土壤特征的流域蓄水容量模型 . 水科学进展, 2013, 24(5): 651-657.
    [12] 司伟, 包为民, 瞿思敏.  洪水预报产流误差的动态系统响应曲线修正方法 . 水科学进展, 2013, 24(4): 497-503.
    [13] 宋晓猛, 孔凡哲, 占车生, 韩继伟.  基于统计理论方法的水文模型参数敏感性分析 . 水科学进展, 2012, 23(5): 642-649. doi: CNKI: 32.1309.P.20120614.2158.003
    [14] 阚光远, 刘志雨, 李致家, 姚成, 周赛.  新安江产流模型与改进的BP汇流模型耦合应用 . 水科学进展, 2012, 23(1): 21-28. doi: CNKI:32.1309.P.20120104.2012.004
    [15] 刘金涛, 梁忠民.  坡地径流入渗机制对水文模拟的影响分析 . 水科学进展, 2009, 20(4): 467-472.
    [16] 梁忠民, 施晔, 李彬权, 余钟波.  基于统计理论的产流模型 . 水科学进展, 2009, 20(6): 789-793.
    [17] 张金存, 芮孝芳.  分布式水文模型构建理论与方法述评 . 水科学进展, 2007, 18(2): 286-292.
    [18] 高龙华, 张文海.  非点源污染空间连续分布模型研究 . 水科学进展, 2007, 18(3): 439-443.
    [19] 张小峰, 许全喜, 裴莹.  流域产流产沙BP网络预报模型的初步研究 . 水科学进展, 2001, 12(1): 17-22.
    [20] 颜开, 陈树娥, 陈信华.  消退曲线超渗产流模型 . 水科学进展, 1997, 8(1): 90-93.
  • 加载中
图(4) / 表 (3)
计量
  • 文章访问数:  326
  • HTML全文浏览量:  60
  • PDF下载量:  58
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-08-01
  • 网络出版日期:  2020-09-09
  • 刊出日期:  2021-01-30

基于Erlang分布蓄水容量曲线的流域产流模型

doi: 10.14042/j.cnki.32.1309.2021.01.012
    基金项目:

    国家重点研发计划资助项目 2016YFC0402708

    中央高校基本科研业务费资助项目 2017KFYXJJ195

    作者简介:

    闫宝伟(1981-), 男, 山东滨州人, 副教授, 主要从事水文基础理论研究。E-mail:bwyan@hust.edu.cn

  • 中图分类号: TV11

摘要: 蓄水容量曲线是反映流域缺水量空间分布不均匀性的特征曲线,对流域产流计算有直接影响。通过对多个典型流域蓄水容量空间分布的分析,发现Erlang分布能更好地拟合流域的蓄水容量曲线,进一步基于Erlang分布进行流域产流的推导,提出了基于Erlang分布蓄水容量曲线的流域产流模型。应用结果表明,基于Erlang分布的流域产流模型,增加了模型的适应性,模拟结果更接近流域实际的产流过程,比新安江模型能取得更高的模拟精度。此外,该产流模型的参数可由流域地形和土壤类型数据估算,为无资料地区的产流计算提供了一种可行的途径。

English Abstract

闫宝伟, 李正坤, 段美壮, 江慧宁, 刘昱. 基于Erlang分布蓄水容量曲线的流域产流模型[J]. 水科学进展, 2021, 32(1): 120-126. doi: 10.14042/j.cnki.32.1309.2021.01.012
引用本文: 闫宝伟, 李正坤, 段美壮, 江慧宁, 刘昱. 基于Erlang分布蓄水容量曲线的流域产流模型[J]. 水科学进展, 2021, 32(1): 120-126. doi: 10.14042/j.cnki.32.1309.2021.01.012
YAN Baowei, LI Zhengkun, DUAN Meizhuang, JIANG Huining, LIU Yu. A watershed runoff yield model based on Erlang distribution water storage capacity curve[J]. Advances in Water Science, 2021, 32(1): 120-126. doi: 10.14042/j.cnki.32.1309.2021.01.012
Citation: YAN Baowei, LI Zhengkun, DUAN Meizhuang, JIANG Huining, LIU Yu. A watershed runoff yield model based on Erlang distribution water storage capacity curve[J]. Advances in Water Science, 2021, 32(1): 120-126. doi: 10.14042/j.cnki.32.1309.2021.01.012
  • 根据蓄满产流原理, 降水首先补充土壤缺水量, 在流域土壤达到田间持水量以后, 流域产流。蓄水容量曲线用于描述流域土壤缺水量空间分布不均匀的问题, 曲线的形状决定着流域产流面积和产流量的多少, 对于径流的准确模拟具有重要意义。流域蓄水容量受流域下垫面地形和土壤等多种因素影响, 通常难以直接测量, 蓄水容量的空间分布更难以确定, 通常是经验假定。在新安江模型中广泛采用B次抛物线表示蓄水容量曲线, 至于为何采用该类线型, 至今尚无合理解释。为能反映雨季和枯季蓄水容量曲线的灵活性, 周买春和Jayawardena[1]提出了双抛物线型蓄水容量曲线, 以反映不同土壤水分并存状态下的蓄水容量。流域蓄水容量曲线同时也是对流域地形地貌特征的一种反映, Gao等[2]通过对美国404个流域的洪水模拟发现, 流域地形地貌对抛物线型蓄水容量曲线的参数B有较大影响。TOPMODEL中的地形指数同样反映了流域的缺水情况, 其累计频率曲线与蓄水容量曲线实质是一致的, 都是对流域缺水量空间不均匀性的一种描述, 郭方等[3]提出了利用TOPMODEL地形指数的分布推求流域蓄水容量曲线的思路。石朋等[4]通过对流域内栅格地形指数的统计分析, 采用对数韦布尔分布曲线描述流域蓄水容量空间分布的不均匀。孔凡哲和宋晓猛[5]根据各子流域的平均坡度计算蓄水容量的空间分布。向小华等[6]则进一步提出了融合土壤和地形的栅格蓄水容量计算方法, 为研究流域蓄水容量的空间分布提供了一种有效工具。Wang[7]针对SCS模型中的土湿比函数没有考虑土壤蓄水容量空间分布不均匀性的问题, 提出了一种新的分布函数来描述土壤蓄水容量的空间分布。Gao等[8]根据植被根系深度随距离最近河道高程的空间分布规律, 建立了以地形为基础的土壤含水量和变源产流面积间的非线性关系, 提出了一种新的产流模型, 为无资料地区的水文模拟提供了一种可行的工具。

    本文尝试通过栅格蓄水容量的计算, 分析多个典型流域蓄水容量的空间分布型式, 选用Erlang分布曲线进行流域产流的推导, 提出基于Erlang分布蓄水容量曲线的流域产流模型, 一方面可以更加真实地反映流域的产流过程, 另一方面, 也为无资料地区的产流计算提供一种可行途径。

    • 流域蓄水容量的实质是流域下垫面的缺水量, 与流域的地形地貌有较大关系, 对流域的产流有直接影响, 由于蓄水容量难以直接测定, 其空间分布更是难以直接确定。考虑到新安江模型中的蓄水容量与TOPMODEL中缺水量的概念相对应, 而TOPMODEL可以根据地形指数和流域平均缺水量计算出各点的缺水量, 进而可以分析流域缺水量或蓄水容量的空间分布。蓄水容量除了与流域地形有关之外, 还受到土壤类型的影响, 如黏性土壤的土壤颗粒间孔隙较小, 相应的土壤缺水量也较小;而砂质土壤的土壤颗粒间孔隙较大, 相应的土壤缺水量也较大。向小华等[6]以van Genuchten模型表述的土壤水分特征曲线为基础, 推导出了融合地形和土壤特征的流域栅格单元蓄水容量的计算公式:

      $$ {W_i} = \left( {{D_i} - {\psi _{\rm{c}}}} \right)\left( {{\theta _{\rm{f}}} - {\theta _{\rm{r}}}} \right) - \left\{ {\frac{{{D_i}\left( {{\theta _{\rm{s}}} - {\theta _{\rm{r}}}} \right)}}{{{{\left[ {1 + {{\left( {\beta {D_i}} \right)}^m}} \right]}^{\frac{1}{m}}}}} - \frac{{{\psi _{\rm{c}}}\left( {{\theta _{\rm{s}}} - {\theta _{\rm{r}}}} \right)}}{{{{\left[ {1 + {{\left( {\beta {\psi _{\rm{c}}}} \right)}^m}} \right]}^{\frac{1}{m}}}}}} \right\} $$ (1)

      式中: Wi为某点的蓄水容量, mm;Di=Szm[max(tpi)-tpi], 为某点的地下水埋深, mm;Szm为流域内非饱和区的最大蓄水深, mm;tpi为地形指数;θrθfθs分别为土壤凋萎含水量、田间持水量和饱和含水量;ψc为土壤水分达到田间持水量时毛管上升高度, mm;βm为土壤水参数。各种土壤水分特征参数的取值见文献[9], Szm则由流域平均蓄水容量试算而得, 流域平均蓄水容量可通过参数率定得出, 无资料地区可采用Gao等[10]和Wang-Erlandsson等[11]提出的MCT方法进行估算。

      按式(1)计算表 1所示典型流域的栅格蓄水容量, 得到各典型流域蓄水容量的空间分布, 如图 1所示。

      表 1  所选典型流域属性表

      Table 1.  Feature attributes of the selected typical basins

      流域名称 流域面积/km2 所在地区 所属流域 气候区
      松柏 184 湖北 长江流域 湿润区
      恒口 924 陕西 长江流域 半湿润区
      栾川 555 河南 黄河流域 半干旱区
      青狮潭 747 广西 珠江流域 湿润区

      图  1  典型流域蓄水容量的空间分布

      Figure 1.  Spatial distribution of watershed storage capacity in the selected basins

    • Erlang分布是一种连续型概率分布, 该分布与指数分布一样多用来表示独立随机事件发生的时间间隔, 常应用于保险业以及排队论中[12]。Erlang分布作为伽马分布的特例, 计算和应用较为便捷, 且对该分布函数的积分易得出显式表达式, 可以简化后续的产流计算, 由此选用Erlang分布作为新的蓄水容量曲线线型, 其概率密度函数为

      $$ {f_n}\left( x \right) = \frac{1}{{\lambda \Gamma \left( n \right)}}{\left( {\frac{x}{\lambda }} \right)^{n - 1}}\exp \left( { - \frac{x}{\lambda }} \right) $$ (2)

      式中: Γ(·)表示伽马函数;nλ为参数; n值为整数。若式中参数n取值为实数, 则式(2)即为伽马分布的密度函数。Erlang分布是一种伽马分布, 表示独立随机事件发生的时间间隔, 水文中广泛应用的Nash瞬时单位线就是伽马分布函数。对上式积分可得其分布函数:

      $$ {E_n}\left( x \right) = \int_0^x {{f_n}\left( t \right)\;{\rm{d}}t} = 1 - \exp \left( { - \frac{x}{\lambda }} \right)\sum\limits_{i = 0}^{n = 1} {\frac{1}{{i!}}{{\left( {\frac{x}{\lambda }} \right)}^i}} $$ (3)

      式中: En(x)为Erlang分布函数。

      根据流域蓄水容量的空间分布, 进一步统计蓄水容量不同量级所占的面积比例, 绘制流域蓄水容量的面积分布图, 并采用正态分布、对数正态分布、指数分布、Erlang分布等多种不同型式的分布函数以及常用的抛物线型曲线进行拟合, 发现Erlang分布拟合效果更优, 如图 2所示, 图中α表示产流面积比。

      图  2  典型流域蓄水容量分布的拟合图

      Figure 2.  Fitness of distribution of the storage capacity in the selected basins

    • 参考抛物线型蓄水容量曲线定义, 则基于Erlang分布的蓄水容量曲线可以表示为

      $$ \alpha = 1 - {E_n}\left( {{W_{{\rm{mm}}}} - W'} \right) = \exp \left( { - \frac{{{W_{{\rm{mm}}}} - W'}}{\lambda }} \right)\sum\limits_{i = 0}^{n - 1} {\frac{1}{{i!}}\left( {\frac{{{W_{{\rm{mm}}}} - W'}}{\lambda }} \right)} $$ (4)

      式中: W′为流域中某点的蓄水容量;Wmm为流域中最大的点蓄水容量。

      采用Erlang分布蓄水容量曲线计算流域产流, 如图 3所示, 则由蓄满产流原理可知, 流域平均蓄水容量

      $$ {W_{\rm{m}}} = \int_0^{{W_{{\rm{mm}}}}} {\left( {1 - \alpha } \right)} \;{\rm{d}}t = {W_{{\rm{mm}}}} - \lambda \sum\limits_{i = 1}^n {{E_i}\left( {{W_{{\rm{mm}}}}} \right)} = {W_{{\rm{mm}}}} - n\lambda $$ (5)

      图  3  流域蓄水容量分布曲线与降雨径流关系

      Figure 3.  Watershed storage capacity distribution curve and its rainfall runoff relationship

      上式的推导用到了Erlang分布的如下性质[13] :

      $$ \int_0^x {\left[ {1 - {E_n}\left( x \right)} \right]} \;{\rm{d}}x = \lambda \sum\limits_{i = 1}^n {{E_i}} \left( x \right) $$ (6)

      同理可得流域初始蓄水容量

      $$ {W_0}{\rm{ = }}\int_0^A {\left( {1 - \alpha } \right){\rm{d}}W'} = A - \lambda \sum\limits_{i = 1}^n {\left[ {{E_i}\left( {{W_{{\rm{mm}}}}} \right) - {E_i}\left( {{W_{{\rm{mm}}}} - A} \right)} \right]} = A - n\lambda + \lambda \sum\limits_{i = 1}^n {{E_i}\left( {{W_{{\rm{mm}}}} - A} \right)} $$ (7)

      式中: A为初始蓄水容量W0对应的纵坐标。为计算产流量, 需首先确定A的值, 在已知W0的情况下, A可以通过求解式(7)而得, 而该式为A的隐式方程, 很难直接求解, 此处采用牛顿迭代法进行求解。令

      $$ \varphi \left( A \right) = A - n\lambda + \lambda \sum\limits_{i = 1}^n {{E_i}\left( {{W_{{\rm{mm}}}} - A} \right) - {W_0}} $$ (8)

      进一步对上式求导可得

      $$ \varphi '\left( A \right) = 1 - \lambda \sum\limits_{i = 1}^n {{f_i}\left( {{W_{{\rm{mm}}}} - A} \right) = {E_n}\left( {{W_{{\rm{mm}}}} - A} \right)} $$ (9)

      根据牛顿迭代公式可知, 经过i次迭代后, 得

      $$ {A_{i + 1}} = {A_i} - \frac{{\varphi \left( {{A_i}} \right)}}{{\varphi '\left( {{A_i}} \right)}} $$ (10)

      当相邻两次迭代值之差小于允许误差时, 便可求得A的数值解。根据蓄满产流原理, 有效降雨Pe产生的径流为

      $$ \begin{array}{l} R = \int_A^{{P_{\rm{e}}} + A} {\alpha {\rm{d}}W'} = \int_A^{{P_{\rm{e}}} + A} {\left[ {1 - {E_n}\left( {{W_{{\rm{mm}}}} - W'} \right)} \right]} {\rm{d}}W'\\ = \lambda \sum\limits_{i = 1}^n {\left[ {{E_i}\left( {{W_{{\rm{mm}}}} - A} \right) - {E_i}\left( {{W_{{\rm{mm}}}} - A - {P_{\rm{e}}}} \right)} \right]} \end{array} $$ (11)

      综上, 若Pe+A < Wmm, 即局部产流时:

      $$ R = \lambda \sum\limits_{i = 1}^n {\left[ {{E_i}\left( {{W_{{\rm{mm}}}} - A} \right) - {E_i}\left( {{W_{{\rm{mm}}}} - A - {P_{\rm{e}}}} \right)} \right]} $$ (12)

      Pe+AWmm, 即全流域产流时:

      $$ R = {P_{\rm{e}}} - {W_{\rm{m}}} + {W_0} $$ (13)

      在产流计算中, 只要给定流域平均蓄水容量Wm和Erlang分布的参数nλ, 即可实现产流计算。在无资料地区, 也可根据式(1)计算出的栅格蓄水容量, 统计出流域蓄水容量的分布曲线, 采用极大似然估计等参数估计方法估算Erlang分布的参数, 从而可以实现无资料地区的产流计算。在资料丰富的地区, 可以根据降雨径流资料直接进行参数率定。

    • 将新安江模型的产流模块替换成上述基于Erlang分布蓄水容量曲线的流域产流模型, 汇流模块仍然采用原算法, 其中, 壤中流和地下径流汇流采用线性水库法, 河网汇流采用滞后演算法, 由此得到改进的新安江模型。将改进后的模型分别应用于松柏、恒口、栾川和青狮潭4个典型流域的洪水模拟, 以分析模型的应用效果。

      选取松柏流域1983—2012年、恒口流域1980—2017年、栾川流域1998—2010年和青狮潭流域2002—2015年的场次洪水用于模型的率定和检验。其中, 模型的参数采用SCE-UA算法进行优化率定, 计算时段长均为1 h, 模型参数模拟结果如表 2所示。

      表 2  模型参数模拟结果

      Table 2.  Simulation results of model parameters

      计算模块 参数 参数含义 松柏 恒口 栾川 青狮潭
      蒸散发计算 K 蒸发折减系数 0.11 1.08 1.06 0.19
      UM 上层张力水蓄水容量/mm 10.00 13.31 11.64 12.42
      LM 下层张力水蓄水容量/mm 20.01 23.20 23.12 21.79
      DM 深层张力水蓄水容量/mm 20.00 23.40 36.34 30.41
      C 深层蒸发系数 0.05 0.09 0.05 0.08
      产流计算 n 蓄水容量曲线参数 3 2 4 9
      λ 蓄水容量曲线参数 8.35 6.15 9.17 5.99
      IM 不透水面积比例 0.01 0.01 0.03 0.01
      水源划分 SM 自由水蓄水容量/mm 29.79 33.81 32.23 34.26
      EX 自由水蓄水容量曲线方次 1.11 1.16 0.50 1.45
      KI 壤中流出流系数 0.37 0.35 0.14 0.10
      KG 地下水出流系数 0.44 0.30 0.39 0.24
      汇流计算 CS 地表水消退系数 0.80 0.80 0.72 0.80
      CI 壤中流消退系数 0.81 0.87 0.91 0.56
      CG 地下水消退系数 0.97 0.99 0.99 0.98

      选取洪峰相对误差、峰现时间误差和洪水总量相对误差以及确定性系数作为模型精度评价指标, 分别计算每场洪水的评价指标值, 表 3给出了各个流域这些评价指标的平均值以及与新安江模型的结果对比。可以看出, 松柏流域和恒口流域的模拟精度与新安江模型相当, 栾川流域和青狮潭流域的模拟精度有一定提高, 尤其是栾川流域的提高较为明显, 率定期确定性系数由0.73提高到0.80, 其他指标也都有不同程度的提高。图 4进一步展示了栾川流域检验期3场洪水的模拟结果, 可以发现, 改进后的模型无论是洪峰还是洪水过程的模拟都取得了较高的精度, 相比新安江模型也都有提高。

      表 3  洪水模拟结果平均值统计及对比

      Table 3.  Mean statistics and comparison of flood simulation results

      时期 流域名称 场次 洪峰相对误差 峰现时间误差/h 洪水总量相对误差 确定性系数
      XAJ E-XAJ XAJ E-XAJ XAJ E-XAJ XAJ E-XAJ
      率定期 松柏 12 0.23 0.22 1.25 1.25 0.18 0.18 0.84 0.86
      恒口 14 0.26 0.26 1.36 1.36 0.24 0.24 0.80 0.80
      栾川 11 0.29 0.27 0.73 0.73 0.37 0.31 0.73 0.80
      青狮潭 11 0.20 0.19 0.64 0.55 0.09 0.07 0.87 0.88
      检验期 松柏 4 0.26 0.26 1.00 0.75 0.25 0.24 0.79 0.80
      恒口 6 0.26 0.26 1.36 1.36 0.24 0.24 0.80 0.80
      栾川 3 0.26 0.21 1.33 1.00 0.14 0.11 0.80 0.83
      青狮潭 3 0.13 0.09 0.67 0.67 0.16 0.15 0.89 0.89
        注:XAJ表示新安江模型;E-XAJ表示基于Erlang分布蓄水容量曲线的新安江模型。

      图  4  栾川流域检验期洪水过程模拟结果

      Figure 4.  Simulated flood hydrographs in the validation period in Luanchuan station

    • 蓄水容量曲线的形状决定着流域产流面积和产流量的多少, 对于径流的模拟有直接影响。通过对多个典型流域蓄水容量空间分布的分析, 选用Erlang分布曲线进行流域蓄水容量曲线的拟合, 提出了基于Erlang分布蓄水容量曲线的流域产流模型, 主要结论如下:

      (1) 基于Erlang分布推导了一种新的流域产流模型, 相比新安江模型增加了1个产流参数, 模型的适应性增强, 模拟结果更接近实际的产流过程, 一定程度上提高了模型的模拟精度。由于增加了模型参数, 可能会增加其他参数如Wm的不确定性, 需进一步分析。

      (2) 基于Erlang分布蓄水容量曲线产流模型的参数, 可以借助流域的地形和土壤类型数据进行估算, 为无资料地区的产流计算提供了一种可行的途径。

参考文献 (13)

目录

    /

    返回文章
    返回