CN108427654A - 一种中型以上淤地坝已淤积库容快速演算方法 - Google Patents
一种中型以上淤地坝已淤积库容快速演算方法 Download PDFInfo
- Publication number
- CN108427654A CN108427654A CN201810079277.6A CN201810079277A CN108427654A CN 108427654 A CN108427654 A CN 108427654A CN 201810079277 A CN201810079277 A CN 201810079277A CN 108427654 A CN108427654 A CN 108427654A
- Authority
- CN
- China
- Prior art keywords
- storage capacity
- deposited
- silt arrester
- sized
- medium
- 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.)
- Granted
Links
- 238000003860 storage Methods 0.000 title claims abstract description 76
- 238000004364 calculation method Methods 0.000 title claims abstract description 23
- 238000004062 sedimentation Methods 0.000 claims abstract description 42
- 238000013461 design Methods 0.000 claims abstract description 14
- 238000011160 research Methods 0.000 claims abstract description 6
- 230000002596 correlated effect Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 description 11
- 239000013049 sediment Substances 0.000 description 6
- 238000005259 measurement Methods 0.000 description 4
- 239000002689 soil Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004162 soil erosion Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000006424 Flood reaction Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种中型以上淤地坝已淤积库容快速演算方法，具体按照以下步骤实施：步骤1、收集研究区若干座中型以上淤地坝的设计文件，拟合已淤积高程～已淤积库容的关系；步骤2、针对步骤1得到的样本关系建立函数关系；步骤3、针对步骤2得到的函数关系拟合数学关系，进而得到已淤积库容，本发明解决了现有技术中存在的因淤地坝设计文件缺失而难以计算已淤积库容且估算工作量大、精度低的问题。
Description
技术领域
本发明属于地理学及水土保持监测技术领域，具体涉及一种中型以上淤地坝已淤积库容快速演算方法。
背景技术
淤地坝是黄土高原地区主要的水土保持工程措施之一，在减少黄河洪水泥沙问题、确保黄河安澜方面发挥了重要作用。明确淤地坝的拦沙量，对于评价淤地坝建成后的拦泥效益和拦泥潜力，分析各下垫面因素的减沙贡献率，预测黄河未来的来沙情况等具有重要意义。
国内外目前推算已淤积库容的方法主要包括以下几种：(1)坝高～库容曲线法，通过实测已淤积高程，查找设计文件中的坝高～库容曲线确定已淤积库容；(2)水文站实测泥沙资料推算法，基于流域范围内水文站多年泥沙观测资料，按水文站控制面积和淤地坝控制面积比推算；(3)泥沙淤积体概化法：根据坝区泥沙淤积体的形状，将横断面概化为规则的锥体，通过测算特征要素，计算泥沙淤积体体积估算；(4)泥沙淤积年限粗估法：根据设计淤积年限计算公式来粗估淤积库容。
这4种方法在估算已淤积库容时均存在一定的局限性：(1)受多种因素影响，大量淤地坝的设计文件缺失或不完整。尤其是早期淤地坝，大多由群众自发修建，无设计文件；(2)部分水文站长系列实测资料不完整，无法很好地提供泥沙历史资料。水文站控制面积远远大于淤地坝控制面积，坡面产生的泥沙到达水文站还存在输移、沉积等多个环节，与淤地坝的实际淤积泥沙存在较大差异。淤地坝一般修建在水土流失较为严重的沟道上，按控制面积比推算默认的前提假设是水文站控制范围内下垫面均一，必然导致估算结果误差较大；(3)黄土高原沟道的地形地貌复杂，泥沙淤积体通常很不规则，往往难以概化。即使粗略估算，需要测量的要素特征很多，工作量较大。因适用性较差，实际工作中很少采用；(4)近年来黄土高原坡面植被明显好转，实施了大量水土流失治理措施，坡面、沟道产沙明显减少，实际淤积年限已远远大于早期的设计淤积年限，该方法已不适用于现阶段淤地坝已淤积库容的估算。
发明内容
本发明的目的是提供一种中型以上淤地坝已淤积库容快速演算方法，解决了现有技术中存在的因淤地坝设计文件缺失而难以计算已淤积库容且估算工作量大、精度低的问题。
本发明所采用的技术方案是，一种中型以上淤地坝已淤积库容快速演算方法，具体按照以下步骤实施：
步骤1、收集研究区若干座中型以上淤地坝的设计文件，拟合已淤积高程～已淤积库容的关系；
步骤2、针对步骤1得到的样本关系建立函数关系；
步骤3、针对步骤2得到的函数关系拟合数学关系，进而得到已淤积库容。
本发明的特点还在于，
步骤1中，各种不同坝高的淤地坝，其已淤积库容与已淤积高程均呈幂函数关系：
y＝axb (1)
式中，y为已淤积库容，单位为万m3；x为已淤积高程，单位为m；a、b为与总坝高和总库容的相关系数。
步骤2具体按照以下步骤实施：
步骤(2.1)、对任一淤地坝而言，在淤满时，必定满足下式：
V＝aHb (2)
式中，H为坝高，V为总库容，当H、V已知时，根据式(2)得a、b关系：
对步骤1中用于建模的若干座淤地坝样本，分别建立函数关系，计算出总坝高和总库容的相关系数a、b，建立a、b关系散点图；
步骤(2.2)、由步骤(2.1)得到的a、b关系散点图可知，a、b呈负相关的对数函数关系，通过非线性函数拟合得到二者的数学关系如下：
b＝1.57-0.16ln(a) (4)
由式(3)和式(4)，可知：
步骤3具体为：
结合式2、式4和式5，可得：
式中，x为已淤积高程。
步骤1中收集研究区30～40座中型以上淤地坝的设计文件。
本发明的有益效果是，一种中型以上淤地坝已淤积库容快速演算方法，只需提供坝高、总库容和已淤积高程，就能较为准确的快速演算出已淤积库容，由评价结果可知，本方法计算的淤地坝已淤库容与实际已淤库容的误差平均值为4.08％，标准差为15.19％，单坝最大误差在30％内。
附图说明
图1(a)是本发明一种中型以上淤地坝已淤积库容快速演算方法中序号为1～5座淤地坝已淤积高程与已淤积库容的关系曲线；
图1(b)是本发明一种中型以上淤地坝已淤积库容快速演算方法中序号为6～10座淤地坝已淤积高程与已淤积库容的关系曲线；
图1(c)是本发明一种中型以上淤地坝已淤积库容快速演算方法中序号为11～15座淤地坝已淤积高程与已淤积库容的关系曲线；
图1(d)是本发明一种中型以上淤地坝已淤积库容快速演算方法中序号为16～20座淤地坝已淤积高程与已淤积库容的关系曲线；
图1(e)是本发明一种中型以上淤地坝已淤积库容快速演算方法中序号为21～25座淤地坝已淤积高程与已淤积库容的关系曲线；
图1(f)是本发明一种中型以上淤地坝已淤积库容快速演算方法中序号为26～30座淤地坝已淤积高程与已淤积库容的关系曲线；
图1(g)是本发明一种中型以上淤地坝已淤积库容快速演算方法中序号为31～35座淤地坝已淤积高程与已淤积库容的关系曲线；
图2是本发明一种中型以上淤地坝已淤积库容快速演算方法中a与b的散点图及关系曲线图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种中型以上淤地坝已淤积库容快速演算方法，具体按照以下步骤实施：
步骤1、收集研究区30～40座中型以上淤地坝的设计文件，拟合已淤积高程～已淤积库容的关系，其中，各种不同坝高的淤地坝，其已淤积库容与已淤积高程均呈幂函数关系：
y＝axb (1)
式中，y为已淤积库容，单位为万m3；x为已淤积高程，单位为m；a、b为与总坝高和总库容的相关系数；
步骤2、针对步骤1得到的样本关系建立函数关系，具体按照以下步骤实施：
步骤(2.1)、对任一淤地坝而言，在淤满时，必定满足下式：
V＝aHb (2)
式中，H为坝高，V为总库容，当H、V已知时，根据式(2)得a、b关系：
对步骤1中用于建模的30～40座淤地坝样本，分别建立函数关系，计算出总坝高和总库容的相关系数a、b，建立a、b关系散点图；
步骤(2.2)、由步骤(2.1)得到的a、b关系散点图可知，a、b呈负相关的对数函数关系，通过非线性函数拟合得到二者的数学关系如下：
b＝1.57-0.16ln(a) (4)
由式(3)和式(4)，可知：
步骤3、针对步骤2得到的函数关系拟合数学关系，进而得到已淤积库容，具体为：
结合式2、式4和式5，可得：
式中，x为已淤积高程。
实施例
不同坝高淤地坝已淤积高程与已淤积库容函数关系列表如表1所示：
表1不同坝高淤地坝已淤积高程与已淤积库容函数关系
根据表1制图1(a)～图1(g)，得到不同坝高淤地坝已淤积高程与已淤积库容的关系曲线，对图1中建模的35座淤地坝样本，分别建立上述式1的函数关系，计算出a、b，建立a、b关系散点图，如图2所示。
本实施例中，选取位于山西省的5座淤地坝，已知每座淤地坝的总坝高H和总库容V，已淤积高程通过实测计算确定，详见表2。总坝高和总库容来源于统计数据，已淤积高程的实测方法参照《水土保持监测理论与方法》(郭索彦，2010)中的加权平均淤积高程法确定，在需要测量淤地坝的淤积面上均匀选取3个断面，每个断面上均匀选取6个测点，实测各断面各测点淤积高程，取6个测点的平均淤积高程计算每个断面的淤积高程，取3个断面的平均淤积高程计算淤地坝的平均淤积高程H已淤。
表2典型淤地坝基本信息表
将每座淤地坝的H、V和H已淤代入式(4)、式(5)，计算得到参数a、b，见表3。将计算得到的a、b代入式(6)，便可计算得出已淤积库容。y1,2…5分别为编号为1,2…5的淤地坝的已淤积库容。
表3典型淤地坝计算参数表
精度分析
为了验证本计算方法的合理性，在淤地坝分布较为集中的山西省西部、陕西省榆林、延安地区，选择了20座淤地坝进行精度验证分析。
通过统计资料，获得每座淤地坝的总坝高H和总库容V，已淤积高程通过实测确定。已淤积高程的实测方法参照《水土保持监测理论与方法》中的加权平均淤积高程法确定，在需要测量淤地坝的淤积面上均匀选取3个断面，每个断面上均匀选取6个测点，实测各断面各测点淤积高程，取6个测点的平均淤积高程计算每个断面的淤积高程，取3个断面的平均淤积高程计算淤地坝的平均淤积高程H已淤。基于搜集的淤地坝设计文件中坝高～库容曲线，按H已淤查找确定实际淤积库容。将H、V和H已淤代入式(6)，得到计算已淤积库容。误差＝(计算已淤库容—实际已淤库容)/实际已淤库容。由评价结果可知，本方法计算的淤地坝已淤库容与实际已淤库容的误差平均值为4.08％，标准差为15.19％，单坝最大误差在30％内，详见表4，能够满足淤地坝淤积量估算的精度需要。
表4验证淤地坝基本信息表
Claims (5)
1.一种中型以上淤地坝已淤积库容快速演算方法，其特征在于，具体按照以下步骤实施：
步骤1、收集研究区若干座中型以上淤地坝的设计文件，拟合已淤积高程～已淤积库容的关系；
步骤2、针对步骤1得到的样本关系建立函数关系；
步骤3、针对步骤2得到的函数关系拟合数学关系，进而得到已淤积库容。
2.根据权利要求1所述的一种中型以上淤地坝已淤积库容快速演算方法，其特征在于，所述步骤1中，各种不同坝高的淤地坝，其已淤积库容与已淤积高程均呈幂函数关系：
y＝atb (1)
式中，y为已淤积库容，单位为万m3；x为已淤积高程，单位为m；a、b为与总坝高和总库容的相关系数。
3.根据权利要求2所述的一种中型以上淤地坝已淤积库容快速演算方法，其特征在于，所述步骤2具体按照以下步骤实施：
步骤(2.1)、对任一淤地坝而言，在淤满时，必定满足下式：
V＝aHb (2)
式中，H为坝高，V为总库容，当H、V已知时，根据式(2)得a、b关系：
对步骤1中用于建模的若干座淤地坝样本，分别建立函数关系，计算出总坝高和总库容的相关系数a、b，建立a、b关系散点图；
步骤(2.2)、由步骤(2.1)得到的a、b关系散点图可知，a、b呈负相关的对数函数关系，通过非线性函数拟合得到二者的数学关系如下：
b＝1.57-0.16ln(a) (4)
由式(3)和式(4)，可知：
4.根据权利要求3所述的一种中型以上淤地坝已淤积库容快速演算方法，其特征在于，所述步骤3具体为：
结合式2、式4和式5，可得：
式中，x为已淤积高程。
5.根据权利要求1所述的一种中型以上淤地坝已淤积库容快速演算方法，其特征在于，所述步骤1中收集研究区30～40座中型以上淤地坝的设计文件。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810079277.6A CN108427654B (zh) | 2018-01-26 | 2018-01-26 | 一种中型以上淤地坝已淤积库容快速演算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810079277.6A CN108427654B (zh) | 2018-01-26 | 2018-01-26 | 一种中型以上淤地坝已淤积库容快速演算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108427654A true CN108427654A (zh) | 2018-08-21 |
| CN108427654B CN108427654B (zh) | 2021-10-22 |
Family
ID=63156209
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810079277.6A Expired - Fee Related CN108427654B (zh) | 2018-01-26 | 2018-01-26 | 一种中型以上淤地坝已淤积库容快速演算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108427654B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109685312A (zh) * | 2018-11-19 | 2019-04-26 | 西安理工大学 | 一种流域次降雨事件下淤地坝系失效风险评价方法 |
| CN114299136A (zh) * | 2021-12-31 | 2022-04-08 | 西北农林科技大学 | 淤地坝泥沙淤积量的测量方法、装置、计算机及存储介质 |
| CN117495939A (zh) * | 2023-09-28 | 2024-02-02 | 黄河流域水土保持生态环境监测中心 | 基于dsm数据快速计算淤地坝库容及高程信息的方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101205710A (zh) * | 2007-12-20 | 2008-06-25 | 马向前 | 一种治理水土流失的流域堰系及其构筑方法 |
| CN102852114A (zh) * | 2012-09-13 | 2013-01-02 | 中国水电顾问集团北京勘测设计研究院 | 一种水库泥沙淤积的计算方法 |
| US20130166262A1 (en) * | 2008-12-15 | 2013-06-27 | Chevron U.S.A. Inc. | System And Method For Evaluating Dynamic Heterogeneity In Earth Models |
| CN105654488A (zh) * | 2015-12-31 | 2016-06-08 | 西北农林科技大学 | 一种基于gis的淤地坝拦淤量估算方法 |
| CN106815450A (zh) * | 2017-04-12 | 2017-06-09 | 西安理工大学 | 一种淤地坝系减沙量的计算方法 |
-
2018
- 2018-01-26 CN CN201810079277.6A patent/CN108427654B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101205710A (zh) * | 2007-12-20 | 2008-06-25 | 马向前 | 一种治理水土流失的流域堰系及其构筑方法 |
| US20130166262A1 (en) * | 2008-12-15 | 2013-06-27 | Chevron U.S.A. Inc. | System And Method For Evaluating Dynamic Heterogeneity In Earth Models |
| CN102852114A (zh) * | 2012-09-13 | 2013-01-02 | 中国水电顾问集团北京勘测设计研究院 | 一种水库泥沙淤积的计算方法 |
| CN105654488A (zh) * | 2015-12-31 | 2016-06-08 | 西北农林科技大学 | 一种基于gis的淤地坝拦淤量估算方法 |
| CN106815450A (zh) * | 2017-04-12 | 2017-06-09 | 西安理工大学 | 一种淤地坝系减沙量的计算方法 |
Non-Patent Citations (6)
| Title |
|---|
| LIU XIAOYAN: "Response of sediment yield to vegetation restoration at a large spatial scale in the Loess Plateau", 《SCIENCE CHINA TECHNOLOGICAL SCIENCES》 * |
| 刘晓燕 等: "黄土高原仍有拦沙能力的淤地坝数量及分布", 《人民黄河》 * |
| 周庆东: "根据水库总淤积量修改库容曲线的方法与应用", 《水利管理技术》 * |
| 尚桢 等: "水土保持淤地坝特性曲线绘制软件研发", 《甘肃水利水电技术》 * |
| 王雄世 等: "小水库库容及淤积体的数学模型与简单测算", 《人民黄河》 * |
| 高洪波 等: "基于非线性拟合的水库泥沙预报模型研究——以丹江口水库为例", 《企业经济》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109685312A (zh) * | 2018-11-19 | 2019-04-26 | 西安理工大学 | 一种流域次降雨事件下淤地坝系失效风险评价方法 |
| CN109685312B (zh) * | 2018-11-19 | 2023-09-29 | 西安理工大学 | 一种流域次降雨事件下淤地坝系失效风险评价方法 |
| CN114299136A (zh) * | 2021-12-31 | 2022-04-08 | 西北农林科技大学 | 淤地坝泥沙淤积量的测量方法、装置、计算机及存储介质 |
| CN114299136B (zh) * | 2021-12-31 | 2023-03-24 | 西北农林科技大学 | 淤地坝泥沙淤积量的测量方法、装置、计算机及存储介质 |
| CN117495939A (zh) * | 2023-09-28 | 2024-02-02 | 黄河流域水土保持生态环境监测中心 | 基于dsm数据快速计算淤地坝库容及高程信息的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108427654B (zh) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Di Stefano et al. | Flow resistance equation for rills | |
| Guo et al. | An experimental study on the effects of grass root density on gully headcut erosion in the gully region of China's Loess Plateau | |
| Bruno et al. | Field investigation on rilling in the experimental Sparacia area, South Italy | |
| Ferro et al. | Applying hypothesis of self-similarity for flow-resistance law in Calabrian gravel-bed rivers | |
| Verstraeten et al. | Variability of dry sediment bulk density between and within retention ponds and its impact on the calculation of sediment yields | |
| Rejman et al. | Rill characteristics and sediment transport as a function of slope length during a storm event on loess soil | |
| CN108427654A (zh) | 一种中型以上淤地坝已淤积库容快速演算方法 | |
| Pizzuto | An empirical model of event scale cohesive bank profile evolution | |
| Nearing et al. | Determining erodibility parameters from rangeland field data for a process-based erosion model | |
| Todisco et al. | Practical thresholds to distinguish erosive and rill rainfall events | |
| CN113297753B (zh) | 山区河流透过性坝沿程水深预测方法 | |
| Zwolenik et al. | Effect of water surface slope and friction slope on the value of the estimated Manning’s roughness coefficient in gravel-bed streams | |
| García‐Serrana et al. | Description of soil micro‐topography and fractional wetted area under runoff using fractal dimensions | |
| Cheng et al. | Adjustment of flood discharge capacity with varying boundary conditions in a braided reach of the Lower Yellow River | |
| CN108960680A (zh) | 一种泥石流松散物源活跃性判定方法 | |
| Damangir et al. | System identification and subsequent discharge estimation based on level data alone—Gradually varied flow condition | |
| Sidorchuk et al. | Variability of shallow overland flow velocity and soil aggregate transport observed with digital videography | |
| Bulygina et al. | DWEPP: a dynamic soil erosion model based on WEPP source terms | |
| Choo et al. | Estimation of river discharge using mean velocity equation | |
| CN111666667A (zh) | 一种游荡性河道造床流量确定方法 | |
| Schügerl et al. | Effect of aquatic vegetation on Manning roughness coefficient value-case study at the Šúrsky channel | |
| Ebissa | Estimation of Open Channel Roughness by using Gradual Varied Flow Profiles | |
| Radoslav | Field study for determine Mannings roughness coefficient with different flow conditions | |
| Khalifa et al. | Experimental Study to Determine Flow Parameters over Roughed Crump Weir Models | |
| Fard et al. | Study factors influencing the hydraulic roughness coefficient of the Karun river (Iran) |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211022 |