CN113240171B - 河口叶脉状河网干流最优引水流量确定方法 - Google Patents
河口叶脉状河网干流最优引水流量确定方法 Download PDFInfo
- Publication number
- CN113240171B CN113240171B CN202110506833.5A CN202110506833A CN113240171B CN 113240171 B CN113240171 B CN 113240171B CN 202110506833 A CN202110506833 A CN 202110506833A CN 113240171 B CN113240171 B CN 113240171B
- Authority
- CN
- China
- Prior art keywords
- flow
- maximum
- river
- main flow
- ammonia nitrogen
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 21
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims abstract description 34
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 20
- 231100000719 pollutant Toxicity 0.000 claims abstract description 20
- 238000011010 flushing procedure Methods 0.000 claims abstract description 12
- 239000010865 sewage Substances 0.000 claims abstract description 11
- 238000000611 regression analysis Methods 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 7
- 238000003911 water pollution Methods 0.000 abstract description 2
- 230000003313 weakening effect Effects 0.000 abstract 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 239000001301 oxygen Substances 0.000 description 4
- 229910052760 oxygen Inorganic materials 0.000 description 4
- 230000000630 rising effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000000746 purification Methods 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 208000001848 dysentery Diseases 0.000 description 1
- 238000012851 eutrophication Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000009991 scouring Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- 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/10—Services
- G06Q50/26—Government or public services
-
- 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
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/40—Protecting water resources
- Y02A20/402—River restoration
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Civil Engineering (AREA)
- Health & Medical Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Educational Administration (AREA)
- Ocean & Marine Engineering (AREA)
- Mechanical Engineering (AREA)
- Primary Health Care (AREA)
- Structural Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种河口叶脉状河网干流最优引水流量确定方法。本发明分四部分，一是调查收集河网平面分布及河流几何特征，确定河网干流加权坡降、加权河宽；二是确定口门断面关键水力要素，包括潮差与最大涨潮和落潮单宽流量关系、口门最大水深、最小水深、最大涨潮流速和最大落潮流速；三是确定干流污染源氨氮最大平均排放浓度和排放流量沿程分布规律及峰值；四是确定干流最优引水冲污流量及相应最大平均氨氮浓度。河口叶脉状河网干流最优引水流量确定方法，给出在潮流与径流复杂博弈条件下，削弱河道污染物氨氮浓度最佳效果的引水流量求解方法，解决了调水引流作为改善平原河网水污染问题最为快捷和经济可行方法中的核心技术难题。
Description
技术领域
本发明涉及水力学及河流动力学领域，具体的说是河口叶脉状河网干流最优引水流量确定方法。
背景技术
我国地域辽阔，水系发达，河网密布，伴随着社会经济的快速发展和城镇化步伐的加快，河网的污染物治理问题日益突出，据不完全统计，截止2016年，全国295个地级以上城市中有216座城市查出黑臭水体1811个。河网污染物治理问题已成为关系到社会文明、城市发展、人民健康的重大生态课题。控源截污是改善河网水环境的最根本措施，但实施难度大，调水冲污是解决平原河网水污染问题最为快捷和经济可行的技术方法。20世纪60年代，日本东京母亲河隅田川BOD质量浓度曾一度达到40mg/L，成为臭气熏天的黑臭河流，为此，日本从利根川引水50m3/s，经武藏野水路入荒川，最后经新河岸川引入隅田川，1990年隅田川生态需氧量BOD质量浓度下降了90％，收到明显成效。此外，美国密西西比河、德国鲁尔河、俄罗斯莫斯科河引水工程，我国的福州市城区内河引水冲污工程、杭州市钱塘江引水入西湖工程、成都市府南河综合整治工程、上海市苏州河综合调水工程等均利用引水改善河道水动力条件，增加水体复氧量，提高了水体自净能力，对修复水环境起了重要作用。然而，究竟引入多大流量对污染物治理效果最好，这一直是该工程措施的主要技术瓶颈。引入流量太大，不仅浪费宝贵水资源，增加工程投资，而且对防洪和天然河道冲刷造成不利影响，引入流量太小，对水环境的改善效果不明显。通过对河流自然属性、污染物排放沿程分布、潮流特性以及引水流量对水动力条件和污染物运移的影响研究发现，某一引水流量下，污染物浓度能够降到最低，再增加引水流量，并不能继续降低污染物浓度，换言之，调水冲污工程存在一个最优引水流量。
从包括中国专利在内的有关资料检索表明，目前国内外尚无河口叶脉状河网干流最优引水流量确定方法的相关报道。
发明内容
(1)要解决的技术问题
平原地区河网，特别是河口人口密集的城市河网，因大量污水排入河道，且水动力不足，河网流速缓慢，水流流动性差，导致污染物常滞留在河道中无法顺利排出，造成水体黑臭，加上潮流的逆向作用和纵横交错的复杂边界条件，使得水环境治理成为突出的技术难题。本发明要解决的主要技术问题是：根据河网特点、潮径流条件及污染物分布，确定调水冲污的最优流量，以获取河网污染物治理的最佳效果。
(2)技术方案
本发明解决其技术问题所采用的技术方案是：河口叶脉状河网干流最优引水流量确定方法，其特征是：其特征是：河口叶脉状河网最优引水流量确定方法分五部分，一是调查收集河网平面分布及河流几何特征，确定河网干流加权坡降、加权河宽；二是确定口门断面关键水力要素，包括潮差与最大涨潮和落潮单宽流量关系、口门最大水深、最小水深、最大涨潮流速和最大落潮流速；三是确定干流污染源氨氮最大平均排放浓度和排放流量沿程分布规律及峰值；五是确定干流最优引水冲污流量及相应最大平均氨氮浓度。
河网干流加权坡降、加权河宽由下式确定：
干流加权坡降
干流加权河宽
其中Lms为干流总计算河长，lms(i)为干流第i段河长(m)，rsms(i)为干流第i段坡降(单位：百分之一)，rbms(i)为干流第i段河宽，nms为干流总分段数。
口门断面关键水力要素由以下各式确定口门最大潮差：TRMAXms0＝SLmax-SLmin (3)其中SLmax为口门最高潮位，单位m，SLmin为口门最低潮位，单位m涨潮口门系数：
B0为口门净宽，单位m。
口门最大涨潮单宽流量：
qrmax单位m3/s·m，g为重力加速度，取9.81m/s2。
口门最大涨潮流量：Qrmax＝qrmax×B0(6)Qrmax单位m3/s。
落潮口门系数：
口门最大落潮单宽流量：
qemax单位m3/s·m。
口门最大落潮流量：Qemax＝qemax×B0(9)
Qemax单位m3/s。
口门最大水深：
Hmax0＝SLmax-SLb0 (10)
其中SLb0为口门河床平均底高程，单位m。
口门最小水深：
Hmin0＝SLmin-SLb0 (11)
口门最大涨潮流速：
vrmax单位m/s。
口门最大落潮流速：
vemax单位m/s。
潮区界断面涨潮最大水深：
Hmax1＝η×[SLmax-(SLb0+LTRMAXms×RSms)] (14)
LTRMAXms为干流潮区界距干流口门距离m，η为沿程阻力系数，取0.6532。
干流污染源氨氮最大平均排放浓度沿程分布方程：
Canms＝Ams1×Li 3+Ams2×Li 2+Ams3×Li+Ams4 (15)
式中Li为干流各污染源距口门距离，单位m。系数Ams1、Ams2、Ams3、Ams4通过原型观测取得的污染物排放分布实测值，经一元三次方程回归分析确定。
干流污染源氨氮最大平均排放浓度出现的位置：
干流污染源氨氮最大平均排放浓度：
CMAXanms＝Ams1×LMAXCms 3+Ams2×LMAXCms 2+Ams3×LMAXCms+Ams4 (17)干流污染源氨氮排放流量沿程分布方程：
Qanms＝Bms1×Li 3+Bms2×Li 2+Bms3×Li+Bms4 (18)
式中Li为干流各污染源距口门距离，单位m。系数Bms1、Bms2、Bms3、Bms4通过原型观测取得的污染物排放流量分布实测值，经一元三次方程回归分析确定。
干流污染源氨氮沿程最大排放流量出现的位置：
干流污染源氨氮沿程最大排放流量：
QMAXanms＝Bms1×LMAXQms 3+Bms2×LMAXQms 2+Bms3×LMAXQms+Bms4 (20)最优引水冲污流量由下式确定
其中系数C1optimum由下式确定
λ为潮力与污染物相互作用影响因子，取0.1756
系数C2optimum由下式确定：
其中C3optimum为径流等于0时，干流氨氮最大平均浓度，初设时可取枯季干流氨氮最大平均浓度，
为河流特性与污染物相互作用影响因子，取92.271。
最优引水流量下干流氨氮最大平均浓度：
ANCmax＝C1optimum×Qad 2+C2optimum×Qad+C3optimum (24)
以上公式适用条件：河流天然坡降0～1.25‰，河宽6～52.9m，潮差2.3～3.0m。
(3)本发明的有益效果
河网水体黑臭和富营养化所造成的严重污染，不仅影响社会文明进步和生态环境，而且降低生活质量，甚至直接威胁人民身体健康，阻碍社会经济环境高质量发展。对河网水体进行治理修复，是社会经济发展、城市景观品质和生态环境建设的迫切需要。引进外部清洁水源来改善河道水质，是一种投资少、成本低、见效快的处理工程。调水冲污结合河网和潮流的自然流动规律和污染源沿程分布排放特点，提出最优引水流量，可创造河网优势水动力条件，激发水体运动能力，加快污染物运移，保证水体上下水层之间的持续频繁交流，强化河流自净能力，对安全生态水系建设和健康河流具有重要意义。
附图说明
以下结合附图和实施例对本发明进一步详细说明。
图1是本发明河道平面布置图
图中干流南大塘排涝渠1，支流后山溪2，支流古溪溪3，支流后岗溪4
具体实施方式
实施例1：
第一步：调查收集河网平面分布及河流几何特征，确定河网干流加权坡降、加权河宽
根据干流南大塘排涝渠25个实测断面河流几何特征，见表1。
表1干流南大塘排涝渠实测河流几何特征
| 序号 | 断面名称 | 断面桩号 | 平均河宽 | 河床平均底高程 |
| 1 | ndt1 | 0-303 | 16 | -1.07 |
| 2 | 干4源 | 0-219.38 | 46.4 | -1.04 |
| 3 | ndt2 | 0-137 | 68.1 | -1.17 |
| 4 | ndt3 | 0+000 | 27.7 | -1.04 |
| 5 | ndt4 | 0+368 | 27.9 | -0.86 |
| 6 | ndt5 | 0+639 | 18.2 | -0.76 |
| 7 | 干3源 | 0+660.67 | 18 | -0.77 |
| 8 | ndt6 | 0+849 | 18 | -0.68 |
| 9 | ndt7 | 1+284 | 16 | -0.55 |
| 10 | ndt8 | 1+574 | 16 | -0.51 |
| 11 | 干2源 | 1+664.69 | 16 | -0.51 |
| 12 | ndt9 | 1+819 | 16 | -0.49 |
| 13 | ndt10 | 2+155 | 16 | -0.46 |
| 14 | ndt11 | 2+346 | 14 | -0.43 |
| 15 | 干1源 | 2+608.24 | 14 | -0.19 |
| 16 | ndt12 | 2+686 | 14 | -0.1 |
| 17 | ndt13 | 3+018 | 14 | 0.21 |
| 18 | ndt14 | 3+058 | 6 | 0.26 |
| 19 | ndt15 | 3+190 | 6 | 0.32 |
| 20 | ndt16 | 3+305 | 6 | 0.45 |
| 21 | ndt17 | 3+436.5 | 6 | 0.55 |
| 22 | ndt18 | 3+565 | 6 | 0.66 |
| 23 | 支1源 | 3+673.37 | 6 | 0.74 |
| 24 | ndt19 | 3+690 | 6 | 0.75 |
| 25 | ndt20 | 3+798 | 6 | 0.85 |
根据(1)式计算干流南大塘排涝渠加权坡降RSms＝0.046818，根据(2)式计算干流加权河宽RBms＝17.701m。
第二步：确定口门断面关键水力要素
根据式(3)式计算口门最大潮位差
TRMAXms0＝SLmax-SLmin＝2.877-(-0.118)＝2.995(m)
根据(4)式计算涨潮口门系数
根据(5)式计算口门最大涨潮单宽流量：
根据(6)式计算口门最大涨潮流量：
Qrmax＝qrmax×B0＝4.432×11.2＝49.638(m3/s)
根据(7)式计算落潮口门系数：
根据(8)式计算口门最大落潮单宽流量：
根据(9)式计算口门最大落潮流量：
Qemax＝qemax×B0＝2.411×11.2＝27.003(m3/s)
根据(10)式计算口门最大水深为
Hmax0＝SLmax-SLb0＝2.877-(-1.07)＝3.947(m)
根据(11)式计算口门最小水深为
Hmin0＝SLmin-SLb0＝-0.118-(-1.07)＝0.952(m)
根据(12)式计算口门最大涨潮流速
根据(13)式计算口门最大落潮流速
根据(14)式计算潮区界断面涨潮最大水深
Hmax1＝η×[SLmax-(SLb0+LTRMAXms×RSms)]
＝0.6532×[2.877-(-1.07+5691×0.046818/100)]
＝0.6532×[2.877-1.594]＝0.6532×1.283＝0.838(m)
根据原型观测资料表2中技术数据，通过一元三次方程回归分析，确定式(15)中干流污染源氨氮最大平均排放浓度方程各系数，即Ams1
＝2.5448，Ams2＝-17.159，Ams3＝30.119，Ams4＝20.337
表2干流南大塘排涝渠沿程污染源排放浓度分布
| 断面桩号 | 排放浓度(mg/l) | 污染源距始点距离(km) | 排放浓度沿程密度(mg/l·km) |
| 0-219.38 | 20 | 3.89275 | 27.22295875 |
| 0+660.67 | 20 | 3.01225 | 26.5582206 |
| 1+664.69 | 20 | 2.00868 | 29.87036263 |
| 2+608.24 | 20 | 1.06513 | 37.55410138 |
| 3+673.37 | 20 | 0 | 20 |
得到干流南大塘排涝渠污染源氨氮最大平均排放浓度沿程变化规律如下：
Canms＝Ams1×Li 3+Ams2×Li 2+Ams3×Li+Ams4
＝2.5448×Li 3-17.159×Li 2+30.119×Li+20.337
根据式(16)确定干流南大塘排涝渠污染源氨氮最大平均排放出现位置为
CMAXanms＝Ams1×Lms 3+Ams2×Lms 2+Ams3×Lms+Ams4
＝2.5448×LMAXCms 3-17.159×LMAXCms 2+30.119×LMAXCms+20.337
＝2.5448×1.1963-17.159×1.1962+30.119×1.196+20.337＝36.168(mg/l·km)
干流峰值位于上游首个污染源排放点下游1.196km处，即距口门2.792km，干流污染源最大排放密度CMAXanms＝36.168mg/l·km。
根据原型观测资料表3中技术数据，通过一元三次方程回归分析，确定式(17)中干流污染源氨氮排放流量方程各系数，
表3干流南大塘排涝渠沿程氨氮最大平均浓度污染源排放流量分布
| 断面桩号 | 排放流量(m<sup>3</sup>/s) | 污染源距始点距离(km) | 排放流量沿程密度(m<sup>3</sup>/s·km) |
| 0-219.38 | 0.05 | 3.89275 | 0.068057397 |
| 0+660.67 | 0.05 | 3.01225 | 0.066395551 |
| 1+664.69 | 0.05 | 2.00868 | 0.074675907 |
| 2+608.24 | 0.05 | 1.06513 | 0.093885253 |
| 3+673.37 | 0.05 | 0 | 0.05 |
即Bms1＝0.0064，Bms2＝-0.0429，Bms3＝0.0753，Bms4＝0.0508，得到干流南大塘排涝渠污染源氨氮沿程排放最大流量变化规律如下：
Qanms＝0.0064×Li 3-0.0429×Li 2+0.0753×Li+0.0508
干流南大塘排涝渠污染源氨氮沿程排放最大流量峰值位置：
氨氮沿程排放最大流量：
QMAXanms＝Bms1×LMAXQms 3+Bms2×LMAXQms 2+Bms3×LMAXQms+Bms4
＝0.0064×1.23-0.0429×1.22+0.0753×1.2+0.0508
＝0.09044(m3/s·km)
系数C1optimum：
已知枯季干流南大塘排涝渠最大平均氨氮浓度C3optimum为19.059(mg/l)
系数C2optimum：
干流最优引水冲污流量关系由下式确定：
则在最优引水流量条件下，干流南大塘排涝渠相应最大平均氨氮浓度
ANCmax＝C1optimum×Qoptimum 2+C2optimum×Qoptimum+C3optimum
＝3.102×2.3732-14.725×2.373+19.059＝1.584(mg/l)
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
Claims (2)
1.一种河口叶脉状河网干流最优引水流量确定方法，其特征是：河口叶脉状河网干流最优引水流量确定方法分四部分，一是调查收集河网平面分布及河流几何特征，确定河网干流加权坡降、加权河宽；二是确定口门断面关键水力要素，包括潮差与最大涨潮和落潮单宽流量关系、口门最大水深、最小水深、最大涨潮流速和最大落潮流速；三是确定干流污染源氨氮最大平均排放浓度和排放流量沿程分布规律及峰值；四是确定干流最优引水冲污流量及相应最大平均氨氮浓度；
河网干流加权坡降、加权河宽由下式确定：
干流加权坡降
干流加权河宽
其中Lms为干流总计算河长，lms(i)为干流第i段河长(m)，rsms(i)为干流第i段坡降(单位：百分之一)，rbms(i)为干流第i段河宽，nms为干流总分段数；
干流污染源氨氮最大平均排放浓度和排放流量沿程分布规律及峰值位置由以下各式确定：
干流污染源氨氮最大平均排放浓度沿程分布方程：
Canms＝Ams1×Li 3+Ams2×Li 2+Ams3×Li+Ams4 (15)
式中Li为干流各污染源距口门距离，单位m，系数Ams1、Ams2、Ams3、Ams4通过原型观测取得的污染物排放分布实测值，经一元三次方程回归分析确定；
干流污染源氨氮最大平均排放浓度出现的位置：
干流污染源氨氮最大平均排放浓度：
CMAXanms＝Ams1×LMAXCms 3+Ams2×LMAXCms 2+Ams3×LMAXCms+Ams4 (17)
干流污染源氨氮排放流量沿程分布方程：
Qanms＝Bms1×Li 3+Bms2×Li 2+Bms3×Li+Bms4 (18)
式中Li为干流各污染源距口门距离，单位m，系数Bms1、Bms2、Bms3、Bms4通过原型观测取得的污染物排放流量分布实测值，经一元三次方程回归分析确定；
干流污染源氨氮沿程最大排放流量出现的位置：
干流污染源氨氮沿程最大排放流量：
QMAXanms＝Bms1×LMAXQms 3+Bms2×LMAXQms 2+Bms3×LMAXQms+Bms4 (20)
最优引水冲污流量由下式确定：
其中系数C1optimum由下式确定：
λ为潮力与污染物相互作用影响因子，取0.1756；Qrmax为口门最大涨潮流量；Qemax为口门最大落潮流量；
系数C2optimum由下式确定：
其中C3optimum为径流等于0时，干流氨氮最大平均浓度，初设时取枯季干流氨氮最大平均浓度，
为河流特性与污染物相互作用影响因子，取92.271；Hmax0为口门最大水深，Hmax1为潮区界断面涨潮最大水深，vrmax为口门最大涨潮流速，vemax为口门最大落潮流速；
最优引水冲污流量下干流污染源氨氮最大平均浓度：
ANCmax＝C1optimum×Qoptimum 2+C2optimum×Qoptimum+C3optimum (24)
以上公式适用条件：河流天然坡降0～1.25‰，河宽6～52.9m，潮差2.3～3.0m。
2.根据权利要求1所述的河口叶脉状河网干流最优引水流量确定方法，其特征是：口门断面关键水力要素由以下各式确定：
口门最大潮差：TRMAXms0＝SLmax-SLmin (3)
其中SLmax为口门最高潮位，单位m，SLmin为口门最低潮位，单位m；
涨潮口门系数：
B0为口门净宽，单位m；
口门最大涨潮单宽流量：
qrmax单位m3/s·m，g为重力加速度，取9.81m/s2；
口门最大涨潮流量：Qrmax＝qrmax×B0(6)
Qrmax单位m3/s；
落潮口门系数：
口门最大落潮单宽流量：
qemax单位m3/s·m；
口门最大落潮流量：Qemax＝qemax×B0 (9)
Qemax单位m3/s；
口门最大水深：
Hmax0＝SLmax-SLb0 (10)
其中SLb0为口门河床平均底高程，单位m；
口门最小水深：
Hmin0＝SLmin-SLb0 (11)
口门最大涨潮流速：
vrmax单位m/s；
口门最大落潮流速：
vemax单位m/s；
潮区界断面涨潮最大水深：
Hmax1＝η×[SLmax-(SLb0+LTRMAXms×RSms)] (14)
LTRMAXms为干流潮区界距干流口门距离m，η为沿程阻力系数，取0.6532。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110506833.5A CN113240171B (zh) | 2021-05-10 | 2021-05-10 | 河口叶脉状河网干流最优引水流量确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110506833.5A CN113240171B (zh) | 2021-05-10 | 2021-05-10 | 河口叶脉状河网干流最优引水流量确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113240171A CN113240171A (zh) | 2021-08-10 |
| CN113240171B true CN113240171B (zh) | 2022-06-03 |
Family
ID=77132967
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110506833.5A Active CN113240171B (zh) | 2021-05-10 | 2021-05-10 | 河口叶脉状河网干流最优引水流量确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113240171B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106777959A (zh) * | 2016-12-12 | 2017-05-31 | 郑州大学 | 人工干扰无水文资料地区河流环境流量分区界定计算方法 |
| CN107679021A (zh) * | 2017-09-27 | 2018-02-09 | 中国科学院南京地理与湖泊研究所 | 一种河流入湖库流量计算方法 |
| CN108570963A (zh) * | 2018-05-15 | 2018-09-25 | 中国水利水电科学研究院 | 平原河网区非汛期水动力条件改善智能控制方法 |
| CN110458359A (zh) * | 2019-08-13 | 2019-11-15 | 河海大学 | 一种基于未完全截污条件下的平原河网水量水质优化调控方法 |
| CN111027264A (zh) * | 2019-12-03 | 2020-04-17 | 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） | 基于生态恢复目标的平原区城市河网水循环调控方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111581579B (zh) * | 2020-03-06 | 2023-03-24 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种堰塞湖泄流调控方法 |
-
2021
- 2021-05-10 CN CN202110506833.5A patent/CN113240171B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106777959A (zh) * | 2016-12-12 | 2017-05-31 | 郑州大学 | 人工干扰无水文资料地区河流环境流量分区界定计算方法 |
| CN107679021A (zh) * | 2017-09-27 | 2018-02-09 | 中国科学院南京地理与湖泊研究所 | 一种河流入湖库流量计算方法 |
| CN108570963A (zh) * | 2018-05-15 | 2018-09-25 | 中国水利水电科学研究院 | 平原河网区非汛期水动力条件改善智能控制方法 |
| CN110458359A (zh) * | 2019-08-13 | 2019-11-15 | 河海大学 | 一种基于未完全截污条件下的平原河网水量水质优化调控方法 |
| CN111027264A (zh) * | 2019-12-03 | 2020-04-17 | 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） | 基于生态恢复目标的平原区城市河网水循环调控方法 |
Non-Patent Citations (3)
| Title |
|---|
| 一维河网模型在平原感潮河网地区的应用;刘军英 等;《广东化工》;20150430;第42卷(第08期);第135-137页 * |
| 平原河网区突发性水污染事件应急调水数值模拟分析;马小雪 等;《中国农村水利水电》;20150430(第04期);第47-50,59页 * |
| 慈溪市城区河道改善水环境引水配水方案研究;詹平;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅰ辑 (月刊)》;20160515(第05期);第1-63页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113240171A (zh) | 2021-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Qi et al. | Water environmental degradation of the Heihe River Basin in arid northwestern China | |
| CN110188959B (zh) | 一种城市内河水量水质优化调控方法 | |
| CN111807645A (zh) | 一种河道旁侧生态净化系统 | |
| CN113402033B (zh) | 一种黑臭河流梯级治理净化系统及方法 | |
| CN113240171B (zh) | 河口叶脉状河网干流最优引水流量确定方法 | |
| CN205023952U (zh) | 一种垂直流潜流式人工湿地 | |
| CN109912036B (zh) | 一种河道漂浮垃圾收集及水质净化系统 | |
| CN204569552U (zh) | 一种无动力下行流河岸人工湿地水质净化系统 | |
| CN207143844U (zh) | 一种基质更替式生态滤坝 | |
| CN102583740B (zh) | 感潮河涌的水处理装置 | |
| CN102583760A (zh) | 跌跃式薄水层人工湿地复氧槽 | |
| CN202465389U (zh) | 感潮河涌的水处理装置 | |
| CN113688522A (zh) | 流域生态水文及生态补水调控方法 | |
| CN211770572U (zh) | 河流砾间湿地结构 | |
| CN209740826U (zh) | 一种生态岸带净化结构及带有该结构的河道净化系统 | |
| CN203383184U (zh) | 适用于城市河道的阶梯型生态景观拦水坝 | |
| CN103321182B (zh) | 圆弧形生态拦水坝及其修建方法 | |
| CN103195019A (zh) | 一种河流净化系统 | |
| CN208603975U (zh) | 一种城市河道增氧生态净化及生境改善集成系统 | |
| Yujing et al. | River Channel Improvement and River Ecological Restoration. | |
| CN206328220U (zh) | 一种生态浮床 | |
| CN219239407U (zh) | 一种河口旁侧近自然湿地净化系统 | |
| CN214087908U (zh) | 一种生态平衡湿地处理系统 | |
| CN109516569A (zh) | 一种黑臭河综合治理方法 | |
| CN203144181U (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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: No.158, Dongda Road, Fuzhou, Fujian 350000 Applicant after: Fujian Water Resources and Hydropower Survey, design and Research Institute Co.,Ltd. Address before: 350000 No. 158, Dongda Road, Fuzhou, Fujian Applicant before: FUJIAN PROVINCIAL INVESTIGATION DESIGN & Research Institute OF WATER CONSERVANCY AND HYDROPOWER |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |