CN111241467A - 一种多沙河流水库排沙比的计算方法 - Google Patents
一种多沙河流水库排沙比的计算方法 Download PDFInfo
- Publication number
- CN111241467A CN111241467A CN202010013687.8A CN202010013687A CN111241467A CN 111241467 A CN111241467 A CN 111241467A CN 202010013687 A CN202010013687 A CN 202010013687A CN 111241467 A CN111241467 A CN 111241467A
- Authority
- CN
- China
- Prior art keywords
- sand
- section
- inlet
- outlet
- calculating
- 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
Images
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
-
- 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
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
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)
- Barrages (AREA)
Abstract
本发明公开了一种多沙河流水库排沙比的计算方法,属于水利工程领域中的水库调度研究技术领域,设定水库坝址上游进口为0‑0断面,下游排沙处出口为1‑1断面,计算进口0‑0断面和出口1‑1断面的水流挟沙力;建立出口断面1‑1的含沙量垂线分布公式;根步骤6中的含沙量垂线分布公式,计算排沙孔处的出库含沙量S2;将出库沙量除以入库沙量即可得到水库的排沙比η。本发明考虑到实际情况中水库的排沙底孔皆偏底部,平均含沙量大于断面平均含沙量,故而用传统排沙比计算的断面平均含沙量导出的排沙比公式计算值偏小,应该根据实际排沙口的位置计算排沙比;本发明可以依据不同的水库特征修改排沙比公式中的参数,以达到本发明中的排沙比公式适用于不同的水库。
Description
技术领域
本发明属于水利工程领域中的水库调度研究技术领域,具体涉及一种多沙河流水库排沙比的计算方法。
背景技术
在多沙河流上的水库都会出现因泥沙淤积而降低有效库容,减少水库的使用寿命。我国针对多沙河流上的水库,保持其库容的方法多为蓄清排浊。在实际水库中,为了高效合理的排水排沙,通常利用排沙比来反映水库的排沙效果,排沙比为出库沙量与入库沙量的比值。研究排沙比的计算方法可以为水库合理的调度泥沙提供有力的依据。
目前,国内外对排沙比的计算方法开展了大量的研究。有学者采用成因分析法,根据水库排沙的特点找出影响排沙比的7个因素,按照线性多元回归计算方法进行数学处理,根据恒山浑水水库的实测资料,建立排沙比与无量纲因素的数学表达式。
但是该经验公式对于其他水库应用时,需要验证和修正参数,并且此方法在入库洪水挟带泥沙颗粒较细的地区适宜使用,对于入库泥沙颗粒较粗条件,计算误差较大。也有学者将小浪底水库实测数据代入已有的排沙比公式,分析这些公式的适用性,找出与排沙比有关的因素,最后采取洪水滞留系数、进出库平均流量变化系数、进库平均含沙量三个因素,采用回归分析法得到小浪底水库排沙比的公式。
此外,有学者考虑区间来沙的情况,将三峡水库区间分成若干区域,在满足水文比拟法原则的基础上,计算出三峡水库排沙比的计算公式。由于水库排沙是通过排沙低孔或泄洪洞排沙,排沙的位置靠近水库底部。由含沙量的垂线分布公式可知,在垂线上,越靠近底部,含沙量越高,越靠近水面,含沙量越低。目前,排沙比的计算多采用垂线平均的含沙量代替排沙孔的含沙量,在理论上偏小。
发明内容
发明目的:本发明的目的在于提供一种多沙河流水库排沙比的计算方法,实际情况中大坝的排沙底孔皆偏底部,平均含沙量大于断面平均含沙量,故而用断面平均含沙量导出的排沙比公式计算值偏小;该方法根据均匀沙不平衡输沙基本方程和二维恒定流中,平衡情况下含沙量沿垂线分布的基本微分方程推导出排沙口的含沙量;排沙比即为排沙口出口断面含沙量和进口断面含沙量的比值。
技术方案:为实现上述目的,本发明提供如下技术方案:
一种多沙河流水库排沙比的计算方法,包括如下步骤:
步骤1,设定水库坝址上游进口为0-0断面,下游排沙处出口为1-1断面,测量进出口断面间的距离L,测量出口排沙孔处中心点距河床的距离为y0;
步骤2,统计进口河道单宽流量q0、进口含沙量s0和进出口断面的距离L;
步骤3,测量入库悬移质泥沙的中值粒径D50,计算该粒径的沉速ω;
步骤4,测量进口断面和出口断面的平均水深H0和H1,计算进口0-0断面和出口1-1断面的水流挟沙力,分别记为S0*和S1*;
步骤5,采用下式计算1-1断面的平均含沙量:
其中,s1为出口断面的含沙量、s1*为出口断面的水流挟沙力,单位为kg/m3;s0为进口断面的平均含沙量;s0*为进口断面水流挟沙力,单位为kg/m3;ω为悬移质泥沙沉速,单位为m/s;L为进出口断面间的距离,单位为m;q为进口河道单宽流量,单位为m2/s;
步骤6,建立出口断面1-1的含沙量垂线分布公式;
步骤7,根步骤6中的含沙量垂线分布公式,计算排沙孔处的出库含沙量S2;
步骤8,将出库沙量除以入库沙量即可得到水库的排沙比η。
进一步的,所述的步骤3泥沙沉速的计算方法为:
其中,ω为悬移质泥沙沉速,单位为m/s;D50为悬沙泥沙的粒径,单位为m;ν为运动粘性系数,取0.000001,单位为m2/s。
进一步的,所述的步骤4进口为0-0断面和1-1断面的水流挟沙力S*计算式为:
其中,q为进口河道单宽流量,单位为m/s;g为重力加速度,取9.8m/s2;H为断面的平均水深,进口断面和出口断面的平均水深分别取H0和H1,单位为m;ω为悬移质泥沙沉速,单位为m/s。
进一步的,所述的步骤6进口河道单宽流量q=q0;步骤6出口断面含沙量沿垂线分布公式为:
其中,s2为距河床底部高y处的含沙量,单位为kg/m3;H1为出口断面的平均水深;ω为悬移质沉速,单位为m/s;κ为挟沙水流卡门常数,取κ=0.21;u*为摩阻流速。
进一步的,所述的步骤6摩阻流速u*的计算方法为:
进一步的,所述的步骤7排沙底孔处的含沙量计算方法为:
进一步的,所述的步骤8排沙比η的计算公式即为:
其中:s1为出口断面平均含沙量,单位为kg/m3;ω为悬移质泥沙沉速,单位为m/s;L为进出口断面间的距离,单位为m;s0为进口断面的平均含沙量;κ为挟沙水流卡门常数,取k=0.21;u*为摩阻流速。
有益效果:与现有技术相比,本发明的一种多沙河流水库排沙比的计算方法,根据水库实际情况,考虑到实际情况中水库的排沙底孔皆偏底部,平均含沙量大于断面平均含沙量,故而用传统排沙比计算的断面平均含沙量导出的排沙比公式计算值偏小,应该根据实际排沙口的位置计算排沙比;本发明可以依据不同的水库特征修改排沙比公式中的参数,以达到本发明中的排沙比公式适用于不同的水库。
附图说明
图1是本发明多沙河流水库排沙比的计算方法的流程图。
具体实施方式
下面结合附图和具体实施例对本发明作更进一步的说明。
一种多沙河流水库排沙比的计算方法,包括如下步骤:
步骤1,设定水库坝址上游进口为0-0断面,下游排沙处出口为1-1断面,测量进出口断面间的距离L,测量出口排沙孔处中心点距河床的距离为y0;
步骤2,统计进口河道单宽流量q0、进口含沙量s0和进出口断面的距离L;
步骤3,测量入库悬移质泥沙的中值粒径D50,计算该粒径的沉速ω;
步骤4,测量进口断面和出口断面的平均水深H0和H1,计算进口0-0断面和出口1-1断面的水流挟沙力,分别记为S0*和S1*;
步骤5,采用下式计算1-1断面的平均含沙量:
其中,s1为出口断面的含沙量、s1*为出口断面的水流挟沙力,单位为kg/m3;s0为进口断面的平均含沙量;s0*为进口断面水流挟沙力,单位为kg/m3;ω为悬移质泥沙沉速,单位为m/s;L为进出口断面间的距离,单位为m;q为进口河道单宽流量,单位为m2/s;
步骤6,建立出口断面1-1的含沙量垂线分布公式;
步骤7,根步骤6中的含沙量垂线分布公式,计算排沙孔处的出库含沙量S2;
步骤8,将出库沙量除以入库沙量即可得到水库的排沙比η。
其中,步骤3泥沙沉速的计算方法为:
其中,ω为悬移质泥沙沉速,单位为m/s;D50为悬沙泥沙的粒径,单位为m;ν为运动粘性系数,取0.000001,单位为m2/s。
其中,步骤4进口为0-0断面和1-1断面的水流挟沙力S*计算式为:
其中,q为进口河道单宽流量,单位为m/s;g为重力加速度,取9.8m/s2;H为断面的平均水深,进口断面和出口断面的平均水深分别取H0和H1,单位为m;ω为悬移质泥沙沉速,单位为m/s。
其中,步骤6进口河道单宽流量q=q0;步骤6出口断面含沙量沿垂线分布公式为:
其中,s2为距河床底部高y处的含沙量,单位为kg/m3;H1为出口断面的平均水深;ω为悬移质沉速,单位为m/s;κ为挟沙水流卡门常数,取κ=0.21;u*为摩阻流速。
其中,步骤6摩阻流速u*的计算方法为:
其中,步骤7排沙底孔处的含沙量计算方法为:
其中,步骤8排沙比η的计算公式即为:
其中:s1为出口断面平均含沙量,单位为kg/m3;ω为悬移质泥沙沉速,单位为m/s;L为进出口断面间的距离,单位为m;s0为进口断面的平均含沙量;κ为挟沙水流卡门常数,取k=0.21;u*为摩阻流速。
实施例
如图1所示,本发明多沙河流水库排沙比的计算方法的流程如下:
步骤1,设定水库坝址上游进口为0-0断面,下游排沙处出口为1-1断面,进出口断面间的距离40000m,测量出口排沙孔处中心点距河床的距离为10m;
步骤2,统计进口单宽流量q0为20m3/s、进口含沙量s0为10kg/m3;
步骤3,测量入库泥沙的中值粒径D50=0.00003m,计算该粒径的沉速ω;
其中,ω为悬移质泥沙沉速,单位为m/s;D50为泥沙的粒径,单位为m;ν为运动粘性系数,可取0.000001,单位为m2/s。计算得到泥沙沉速为:0.00058m/s。
步骤4,进口断面平均水深10m,出口断面平均水深50m,采用下式计算进口0-0断面的水流挟沙力S0*和出口1-1断面的水流挟沙力s1*;
其中,U为断面的平均流速,单位为m/s;g为重力加速度,取9.8m/s2;H为断面的平均水深,单位为m。
步骤5,采用下式计算1-1断面的平均含沙量:
其中,s1为出口断面的含沙量、s1*为出口断面的水流挟沙力,单位为kg/m3;s0为进口断面的平均含沙量;s0*为进口断面水流挟沙力,单位为kg/m3;ω为悬移质泥沙沉速,单位为m/s;L为进口断面到出口断面的长度,单位为m;q为河道单宽流量,单位为m2/s。计算得到s1=8.94kg/m3。
步骤6,排沙底孔处的含沙量:
其中,s1为出口断面的含沙量,单位为kg/m3;ω为悬移质泥沙沉速,单位为m/s;L为进口断面到出口断面的长度,单位为m;κ为挟沙水流卡门常数,取k=0.21;H1为出口断面的平均水深,H0为进口断面平均水深,y1为排沙孔至河床的距离。
步骤7,计算水库的排沙比:
其中:s1为出口断面平均含沙量,单位为kg/m3;ω为悬移质泥沙沉速,单位为m/s;L为两断面总长度,单位为m;s0为进口断面的平均含沙量;κ为挟沙水流卡门常数,取k=0.21;u*为摩阻流速。最后综合得到水库排沙比为:0.93。
Claims (7)
1.一种多沙河流水库排沙比的计算方法,其特征在于:包括如下步骤:
步骤1,设定水库坝址上游进口为0-0断面,下游排沙处出口为1-1断面,测量进出口断面间的距离L,测量出口排沙孔处中心点距河床的距离为y0;
步骤2,统计进口河道单宽流量q0、进口含沙量s0和进出口断面的距离L;
步骤3,测量入库悬移质泥沙的中值粒径D50,计算该粒径的沉速ω;
步骤4,测量进口断面和出口断面的平均水深H0和H1,计算进口0-0断面和出口1-1断面的水流挟沙力,分别记为S0*和S1*;
步骤5,采用下式计算1-1断面的平均含沙量:
其中,s1为出口断面的含沙量、s1*为出口断面的水流挟沙力,单位为kg/m3;s0为进口断面的平均含沙量;s0*为进口断面水流挟沙力,单位为kg/m3;ω为悬移质泥沙沉速,单位为m/s;L为进出口断面间的距离,单位为m;q为进口河道单宽流量,单位为m2/s;
步骤6,建立出口断面1-1的含沙量垂线分布公式;
步骤7,根步骤6中的含沙量垂线分布公式,计算排沙孔处的出库含沙量S2;
步骤8,将出库沙量除以入库沙量即可得到水库的排沙比η。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010013687.8A CN111241467B (zh) | 2020-01-07 | 2020-01-07 | 一种多沙河流水库排沙比的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010013687.8A CN111241467B (zh) | 2020-01-07 | 2020-01-07 | 一种多沙河流水库排沙比的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111241467A true CN111241467A (zh) | 2020-06-05 |
CN111241467B CN111241467B (zh) | 2022-06-03 |
Family
ID=70872924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010013687.8A Active CN111241467B (zh) | 2020-01-07 | 2020-01-07 | 一种多沙河流水库排沙比的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111241467B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112149314A (zh) * | 2020-10-16 | 2020-12-29 | 黄河勘测规划设计研究院有限公司 | 一种基于虚拟库容修正的多沙水库库容冲淤模拟方法 |
CN117556184A (zh) * | 2024-01-12 | 2024-02-13 | 长江水利委员会水文局 | 多沙区域水库下游河道输沙量的还原方法、系统及介质 |
CN118013895A (zh) * | 2024-04-10 | 2024-05-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种水库深孔前淤积厚度的计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITRM20070134A1 (it) * | 2007-03-15 | 2008-09-16 | Univ Palermo | Metodo di misura della portata di piena in un corso d'acqua e relativo sistema |
CN106320256A (zh) * | 2016-08-22 | 2017-01-11 | 河海大学 | 一种多沙河流抽水蓄能电站过机泥沙通量的确定方法 |
CN110046469A (zh) * | 2019-05-13 | 2019-07-23 | 水利部交通运输部国家能源局南京水利科学研究院 | 多约束条件下水电站坝前河床冲淤变形的计算方法 |
-
2020
- 2020-01-07 CN CN202010013687.8A patent/CN111241467B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITRM20070134A1 (it) * | 2007-03-15 | 2008-09-16 | Univ Palermo | Metodo di misura della portata di piena in un corso d'acqua e relativo sistema |
CN106320256A (zh) * | 2016-08-22 | 2017-01-11 | 河海大学 | 一种多沙河流抽水蓄能电站过机泥沙通量的确定方法 |
CN110046469A (zh) * | 2019-05-13 | 2019-07-23 | 水利部交通运输部国家能源局南京水利科学研究院 | 多约束条件下水电站坝前河床冲淤变形的计算方法 |
Non-Patent Citations (3)
Title |
---|
JIE QIN ET.AL: "Sediment exchange between groin fields and main-stream", 《ADVANCES IN WATER RESOURCES 》 * |
吴腾 等: "多沙河流水库自校正排沙运用研究", 《水利水电技术》 * |
徐国宾: "水库坝前及泄水底孔内泥沙淤积形态", 《天津大学学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112149314A (zh) * | 2020-10-16 | 2020-12-29 | 黄河勘测规划设计研究院有限公司 | 一种基于虚拟库容修正的多沙水库库容冲淤模拟方法 |
CN112149314B (zh) * | 2020-10-16 | 2021-03-12 | 黄河勘测规划设计研究院有限公司 | 一种基于虚拟库容修正的多沙水库库容冲淤模拟方法 |
CN117556184A (zh) * | 2024-01-12 | 2024-02-13 | 长江水利委员会水文局 | 多沙区域水库下游河道输沙量的还原方法、系统及介质 |
CN117556184B (zh) * | 2024-01-12 | 2024-03-26 | 长江水利委员会水文局 | 多沙区域水库下游河道输沙量的还原方法、系统及介质 |
CN118013895A (zh) * | 2024-04-10 | 2024-05-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种水库深孔前淤积厚度的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111241467B (zh) | 2022-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111241467B (zh) | 一种多沙河流水库排沙比的计算方法 | |
CN110287571B (zh) | 一种河流险工冲刷安全分析与岸坡稳定性判定方法 | |
Eshev et al. | Non-eroding speed of water flow of channels running in cohesive soils | |
Patra et al. | Flow interaction of meandering river with floodplains | |
CN107423484B (zh) | 一种泥石流在拦挡坝后的流量计算方法及其应用 | |
CN103898863B (zh) | 一种研究非淹没刚性植物条件下泥沙起动的方法 | |
CN103761583A (zh) | 一种水库排沙发电多目标优化调度方法 | |
CN113673179A (zh) | 一种长距离泥浆管道输送动态模型及应用 | |
Bazarov et al. | Improving methods of increasing reliability without dam water intake | |
CN110472367B (zh) | 一种多沙河流干支流水沙全交互模拟方法及系统 | |
CN103469772A (zh) | 水库竖井式溢洪道 | |
CN113737710B (zh) | 天然河道推移质输沙率的估算方法 | |
CN104195979A (zh) | 河道交汇水流停滞区楔锥体及其构建方法与应用 | |
CN112989565A (zh) | 一种基于水流挟沙能力的造床流量计算方法 | |
CN111560912B (zh) | 一种封闭式人工沙滩前的潜堤高程确定方法 | |
CN105547638B (zh) | 一种基于压力传感器的薄层水流滚波测量系统与方法 | |
CN113806851B (zh) | 一种疏浚挖槽水动力变化引起的航槽淤积量预测方法 | |
Bishwakarma | Settling basin design criteria and trap efficiency computation methods | |
Ahmad et al. | Estimation of trapped sediment load into a trench weir | |
Lei et al. | A Dynamic Regulation and Control Technology of Sediment Based on Hydraulic Power and Strong Artificial Measure--Pneumatic Sluicing | |
WANG et al. | Theory and practice of hyperconcentrated sediment-laden flow in China | |
Duarte | Hopper system design using CFD–Warsak HPS Water-intake (Pakistan) | |
CN111931276B (zh) | 一种快速优选跌坎消力池体型的方法 | |
CN109934506B (zh) | 一种基于定量化判断航道整治工程类型的选取方法 | |
Khidirov et al. | Justification of parameters of sump in machine channel |
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 |