CN111524035B - 一种面向平原城市河网水环境改善的水动力调控阈值确定方法 - Google Patents
一种面向平原城市河网水环境改善的水动力调控阈值确定方法 Download PDFInfo
- Publication number
- CN111524035B CN111524035B CN202010402856.7A CN202010402856A CN111524035B CN 111524035 B CN111524035 B CN 111524035B CN 202010402856 A CN202010402856 A CN 202010402856A CN 111524035 B CN111524035 B CN 111524035B
- Authority
- CN
- China
- Prior art keywords
- water
- river
- flow velocity
- water environment
- concentration
- 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 228
- 230000033228 biological regulation Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000006872 improvement Effects 0.000 title claims abstract description 17
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 21
- 231100000719 pollutant Toxicity 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 239000000428 dust Substances 0.000 claims abstract description 9
- 238000000746 purification Methods 0.000 claims abstract description 9
- 238000010521 absorption reaction Methods 0.000 claims abstract description 5
- 239000013049 sediment Substances 0.000 claims description 38
- 238000012360 testing method Methods 0.000 claims description 38
- 230000008859 change Effects 0.000 claims description 29
- 241000196324 Embryophyta Species 0.000 claims description 11
- 238000001179 sorption measurement Methods 0.000 claims description 7
- 239000000356 contaminant Substances 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 5
- 230000015556 catabolic process Effects 0.000 claims description 4
- 238000006731 degradation reaction Methods 0.000 claims description 4
- 238000005286 illumination Methods 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 claims description 2
- 239000000843 powder Substances 0.000 claims description 2
- 230000001737 promoting effect Effects 0.000 claims description 2
- 238000011144 upstream manufacturing Methods 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims 3
- 235000012149 noodles Nutrition 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 10
- 230000009471 action Effects 0.000 abstract description 2
- 230000002401 inhibitory effect Effects 0.000 abstract description 2
- 230000007246 mechanism Effects 0.000 abstract description 2
- 239000002245 particle Substances 0.000 description 11
- 230000008569 process Effects 0.000 description 7
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 6
- 229910052698 phosphorus Inorganic materials 0.000 description 6
- 239000011574 phosphorus Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 239000008187 granular material Substances 0.000 description 5
- 229910052739 hydrogen Inorganic materials 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 229910052717 sulfur Inorganic materials 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 244000005700 microbiome Species 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 239000010865 sewage Substances 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 238000012351 Integrated analysis Methods 0.000 description 1
- VVQNEPGJFQJSBK-UHFFFAOYSA-N Methyl methacrylate Chemical compound COC(=O)C(C)=C VVQNEPGJFQJSBK-UHFFFAOYSA-N 0.000 description 1
- 229920005372 Plexiglas® Polymers 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 239000007900 aqueous suspension Substances 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000033558 biomineral tissue development Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000011888 foil Substances 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 230000036284 oxygen consumption Effects 0.000 description 1
- 239000002096 quantum dot Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000006228 supernatant Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000009466 transformation Effects 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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/152—Water filtration
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Activated Sludge Processes (AREA)
- Separation Of Suspended Particles By Flocculating Agents (AREA)
Abstract
本发明涉及一种面向平原城市河网水环境改善的水动力调控阈值确定方法，基于泥水界面的临底流速与切应力作用机理，建立扰动强度与实际河道表面流速的关系，提出以抑制平原城市河道底泥快速释放为准则的水动力调控上限阈值范围；首次采用最小水环境容量理论，考虑河道沿程点源污染、面源污染、直接入河的粉尘污染、底泥污染物释放、河道水体自净、水生植物吸收等多种影响河道水质的因素，利用总体达标方法推导了改善城市水环境的水动力调控下限阈值计算公式。根据提出的平原城市水动力调控阈值确定方法，可以判定平原城市实际水动力调控的合理流速范围，为改善水环境的水动力调控流速合理取值提供依据，提升水资源利用效率、提高水环境改善效果。
Description
技术领域
本发明涉及水利工程领域中的水动力调控理论，特别涉及一种面向平原城市河网水环境改善的水动力调控阈值确定方法。
背景技术
平原城市经济发达、人口密集、城镇化率高，地势低平、河网密布、水动力条件弱。近年来，我国正处于社会经济高速发展的时期，入河污染负荷持续增强，城市河道被挤占、河网分割、水系畅通性差，控源截污不到位、水环境保护措施建设滞后，导致平原城市河网水环境污染问题日趋严重。根据近几年《水环境质量公报》，我国城市中心区大部分河道属于Ⅴ类和劣Ⅴ类，水功能区达标率低，城市河网水环境承载能力不足已经成为制约经济社会可持续发展和宜居环境的瓶颈之一，改善城市河网水环境迫在眉睫。
城市水环境治理是一个综合系统的工程，包括控源截污、清淤疏浚、河道整治、水系连通、动力调控、强化净化、生态修复等城市水环境综合治理措施。其中，水动力调控是改善平原城市河网水环境的“灵魂”，能够调节河道水体置换时间、改善水动力条件、提升河流溶解氧水平、提高水体自净能力，科学有效地增加城市水资源量，提高水环境承载力，使有限的水资源发挥更大的效益，目前，水动力调控已成为河道水质改善和长效维持不可或缺的重要措施之一。
近几年，在平原城市河网区，特别是太湖流域多个城市实施了水动力调控工程，大幅提升了区域河道水体流动性，改善了部分城市的水环境，但是，由于目前水动力调控改善水环境的基础理论不够完善，缺乏指导水环境改善的水动力调控阈值，目前大多数平原城市选用水动力调控方法改善水环境时，将河道流速调节到多大合理、如何调控河道流速，往往都采取试探方式，或者根据经验判断，经常导致河道流速忽高忽低，流向反复变化，水资源利用率不高，水环境改善效果不好。因此，为确定平原河网水动力调控中合理的水动力条件，为更科学高效地利用有限的水资源改善平原城市水环境，急需开展水动力-水质调控阈值确定方法研究，本研究成果将对平原城市水动力调控的实际应用具有重要的指导意义。
发明内容
发明目的：本发明针对平原城市水动力调控改善水环境理论研究的不足，考虑在动力调控中易引起的底泥内源污染释放这一重要因素，并从保障城市水环境的角度，提出改善平原城市水环境的水动力调控阈值确定方法。
技术方案：一种面向平原城市河网水环境改善的水动力调控阈值确定方法，包括以下步骤：
步骤1、在城市河道内采集新鲜底泥和河水水样；
步骤2、建立数学模型；
步骤3、将该切应力作为中间变量，建立室内试验中的扰动转速与野外天然河道表面流速之间的相关关系；
步骤4、获取实际河道表面流速与上覆水水质指标浓度间的响应关系；
步骤5、按照水环境容量计算的理论，根据水质本底值与目标水质，确定水动力调控改善水环境中的流速下限阈值；
步骤6、在步骤5的基础上代入河道进口断面水流流速值、进口和出口断面的水深值，计算得出能够保证河道水质的流速下限阈值。
在进一步的实施例中，步骤1进一步包括：将采样得到的样本避光密封保存并带回实验室进行试验，模拟分析不同水动力条件下的底泥释放与吸附规律，根据试验结果，得到促使底泥快速释放的扰动转速。
在进一步的实施例中，步骤2进一步包括：计算不同扰动强度条件下圆筒试验装置中的流场，提取泥水界面的流速数据，采用切应力计算公式，计算得到装置中泥-水界面的切应力：
式中，τ为切应力，N/m2；u为流速，m/s；y为水深，m；μ为水的动力粘度。
在进一步的实施例中，步骤3进一步包括：采用对数流速分布公式计算实际不同水深河道的表面流速：
式中，u为流速，m/s；u*为摩阻流速，m/s；y为水深，m；Δ为绝对粗糙度，mm；ρ为水的密度；τ为河道底部切应力，N/m2。
在进一步的实施例中，步骤4进一步包括：基于扰动转速与野外天然河道表面流速之间的相关关系和试验中上覆水水质浓度变化，得到实际河道表面流速与上覆水水质指标浓度间的响应关系，以水质变化拐点处的流速作为平原城市河网水动力调控的流速上限阈值。
在进一步的实施例中，步骤5进一步包括：选择总体达标法计算水环境容量，总体达标计算法采用零维模型进行水质计算，考虑点源污染、面源污染、直接入河的粉尘、底泥污染物的释放、河道水体的自净、水中植物对污染物的吸收等多种影响河道水质的因素，由污染物的质量守恒，得到以下公式：
式中，Q上、C上分别为河道上游来水的流量(m3/s)以及水质浓度(mg/L)；W为水环境容量(mg/s)；q点、C点分别为入河点源污染的流量(m3/s)和污染物浓度(mg/L)；q面、C面分别为入河面源污染的流量(m3/s)和污染物浓度(mg/L)；f粉为直接入河粉尘所含污染物的质量函数；rs为底泥污染物的释放速率(mg/(m2*d))；L为河道长度(m)；B为河宽(m)；Cs为河道出口断面的目标水质浓度(mg/L)；K为污染物降解系数(1/d)；V为河段内的水体体积(m3)；f植为水生植物吸收的污染物的质量函数，与光照、温度、植物种类、种植密度以及水深等因素有关。
在进一步的实施例中，步骤6进一步包括：设河道进口断面水流流速为u，进口和出口断面的水深分别为h1和h2，代入步骤5的公式中得：
上式化简得：
式中，若水环境容量W＝0，即该河道水体不再能够承受污染物的排放，此时的流速u为能够保证河道水质的流速下限阈值u小，即：
上式中，各符号含义同上。
有益效果：本发明根据发明提出的平原城市水动力调控阈值确定方法，可以判定平原城市实际水动力调控的合理流速范围，为改善水环境的水动力调控流速合理取值提供依据，提升水资源利用效率、提高水环境改善效果。
附图说明
图1为本发明的研究思路图。
图2为本发明水环境容量计算示意图。
图3为本发明试验装置示意图。
图4为本发明计算模型图。
图5为本发明模型计算网格图。
图6为本发明不同转速条件纵剖面流速分布图。
图7为本发明距轴线12cm位置沿水深方向流速分布图。
图8为本发明转速-流速相关关系图。
图9为本发明水质-水动力相关关系图。
图10.1为上覆水DO浓度变化过程图。
图10.2为上覆水浊度变化图。
图10.3为试验第18天M点上覆水样对比照片。
图10.4为TP浓度和释放速率变化图。
图10.5为NH3-N浓度及释放速率变化图。
图10.6为TN浓度和释放速率变化图。
具体实施方式
申请人认为，目前大多数平原城市选用水动力调控方法改善水环境时，将河道流速调节到多大合理、如何调控河道流速，往往都采取试探方式，或者根据经验判断，经常导致河道流速忽高忽低，流向反复变化，水资源利用率不高，水环境改善效果不好。
为此，本发明基于泥水界面的临底流速与切应力作用机理，建立了扰动强度与实际河道表面流速的关系，提出了以抑制平原城市河道底泥快速释放为准则的水动力调控上限阈值范围；首次采用最小水环境容量理论，考虑河道沿程点源污染、面源污染、直接入河的粉尘污染、底泥污染物释放、河道水体自净、水生植物吸收等多种影响河道水质的因素，利用总体达标方法推导了改善城市水环境的水动力调控下限阈值计算公式。
下面通过具体实例，并结合附图，对本发明的技术方案做进一步具体说明。
一种面向平原城市河网水环境改善的水动力调控阈值确定方法，分为以下步骤：
步骤1，在城市河道内采集新鲜底泥和河水水样，避光密封保存并带回实验室进行室内试验。
本实例在典型平原城市苏州和常州分别采集底泥样四处(H、S、M及X点)，选用有机玻璃圆筒作为反应装置，反应器直径27cm，高50cm，有效容积28.6L。距离反应器上边缘5cm以下部分的侧面贴铝箔纸，模拟河道侧面的避光环境，使光线仅从上方照射，采用电动恒速搅拌器，型号JB60-SH，通过设置不同转速模拟不同的水动力条件，试验装置如图3所示。共设置六种工况，分别为转速100rpm、200rpm、300rpm、400rpm、500rpm以及静置组，为减少试验误差，每种工况设3个平行组，共18套装置同时试验。
试验分析了上覆水体浊度、DO、TP、TN、NH3-N等水质指标的变化规律以及底泥对氮磷营养盐的吸附和释放规律，主要分析成果如下：
1)上覆水DO变化规律
试验期间，H点、S点、M点和X点四种泥样的上覆水溶解氧浓度变化如图10.1所示，对比H点、S点、M点和X点四种泥样的上覆水DO浓度，可以看出不同泥样DO浓度随时间的变化规律类似，扰动工况与静置工况均存在一定差异。扰动工况条件，四种泥样的上覆水DO浓度均随试验的进行而增大，并在扰动刚开始时，DO浓度迅速增加，响应时间短，随着试验的进行，除S点的DO浓度稍有下降外，其余点位均缓慢增加；静置工况条件，试验初期，DO浓度均有不同程度的下降，后期缓慢上升并趋于稳定。上覆水DO浓度主要与微生物活动有关，水体中的氧气为微生物的活动提供条件，影响微生物的有机物降解、矿化速率，污染较轻的河道水体耗氧量小，H点所在河道为景观河道，因此其静置工况的DO浓度降幅较小。
另外，在整个试验期间，H点、S点、M点和X点四种泥样的上覆水DO浓度均随扰动强度的增大而增加，统计结果也显示，四个底泥样本在五种扰动工况下的DO浓度均与静置工况之间差异性极显著(P<0.01)，由此可见，扰动能够增加水体中的DO浓度，流动水体中的DO浓度明显高于静置水体，但扰动工况之间，扰动强度相差较大的工况之间具有显著性差异(P<0.05)，说明增大扰动强度对于提升水体DO浓度是有益的。
2)上覆水浊度变化规律
试验期间，不同工况条件下四个采样点的上覆水浊度变化情况如图10.2所示。对比H点、S点、M点和X点上覆水浊度变化，可以发现，四种浊度随时间的变化趋势不完全相同，H点底泥样本在试验初期受扰动影响较大，颗粒迅速向上覆水悬浮，然后又逐渐下降，且该底泥颗粒粒径相对较大，容易碰撞凝聚成更大颗粒而沉降，而S点、M点和X点三种底泥在试验初期受突加扰动作用的影响小，在试验中后期由于其粒径相对小，当r>400rpm时，颗粒大量起悬，且小颗粒沉速较慢，因此，后期稳定时，其浊度比H点上覆水的浊度大，而且，H点、S点、M点和X点四种底泥的粒径大小关系为：H点>S点>X点>M点，粒径越小的颗粒越容易起悬，因此试验后期稳定时，浊度的大小关系应为：M点>X点>S点>H点，而本文的试验结果也显示出此种关系，说明本文的试验结果合理。另外，根据前文含水率监测成果，M点>X点>S点>H点，因此M点更容易悬浮，X点、S点次之，H点最不容易悬浮。
对比不同工况之间的浊度数值，可以看出，四种底泥试验的上覆水浊度随扰动强度的变化趋势相同，均随扰动强度的增大，浊度逐渐增大，且当转速r＝0～400rpm时，随着扰动强度的增大，四个点位上覆水浊度数值略有上升，增幅不明显，但当转速r＝500rpm时，浊度大幅增加，水体感官明显变差，H点、S点、M点和X点四个点位的500rmp工况时的上覆水浊度最大值比其他扰动工况分别增大约140NTU、270NTU、390NTU和200NTU，由此可见，引起上述四种底泥颗粒大量悬浮的临界扰动强度均介于400～500rpm之间。选取M点在试验后期各扰动强度工况的一个平行样，其上覆水感官状况如图10.3所示，可以看出，静置、100rpm、200rpm和300rpm工况的上覆水几乎清澈见底，400rpm工况有部分底泥颗粒悬浮，500rpm工况则有大量底泥颗粒悬浮，特别是M点，其500rpm工况的上覆水非常浑浊。
另外，统计结果显示，H点、S点、M点和X点四种底泥几乎每个扰动工况的浊度之间均具有显著性差异(P<0.05)，而且，大部分工况之间的差异性极显著(P<0.001)。由此可见，浊度是随扰动条件变化的敏感指标，随扰动强度的增大，浊度逐渐增大。
3)底泥对TP的吸附与释放规律
H点、S点、M点和X点上覆水TP浓度及TP释放速率随时间的变化如图10.4所示。对比各点不同工况之间的TP浓度和TP释放速率可以发现，静置和100～400rpm工况基本一致，由于底泥颗粒未大量起悬，此时底泥中可溶性磷是底泥释放磷的主要途径，TP浓度数值和变化趋势均无较大差异，但500rpm试验工况的TP浓度则明显高于其他组次，这是因为高速扰动的底泥颗粒再悬浮，导致了TP浓度的显著增加，国内外较多学者也研究表明，扰动悬浮引起的总磷增大，是因为悬浮颗粒态磷的增加。由此可见，与浊度分析成果类似，400～500rpm是本文选取的四种底泥上覆水TP变化的临界扰动强度。
H点、S点、M点和X点各点不同工况之间TP浓度的差异性分析统计结果表明，各点500rmp工况的上覆水TP浓度与其他工况间具有显著性差异(P<0.05)，但0～400rmp工况间的上覆水TP浓度不具有显著性差异。因此，TP是随扰动强度变化而改变的敏感指标，TP浓度随扰动强度的增大呈上升趋势，但较小的扰动条件对TP影响不显著，当扰动强度较大时，水体中的TP浓度则会显著增加。
4)底泥对NH3-N的吸附与释放规律
不同扰动工况上覆水NH3-N浓度及释放速率随时间的变化情况如图10.5所示。综合分析不同地点的NH3-N变化规律，可以发现，NH3-N浓度随时间呈先上升后下降或逐渐下降的趋势，不同采样点位由于初始浓度的不同，其变化规律也不尽相同，但对比同一时刻不同工况的NH3-N浓度变化规律，发现不同点位均相似，即随扰动强度增大，NH3-N浓度呈下降趋势，这是因为扰动强度越大，DO浓度越高，硝化作用更迅速，因此，NH3-N下降速度越快，这也说明了增大扰动强度可以控制底泥中NH3-N污染物质的释放。另外，统计学结果显示，各点的静置工况或低扰动工况与高扰动工况的NH3-N浓度差异性极显著(P<0.01)，由此可见NH3-N是随扰动强度改变而变化的敏感指标，但只要稍增大扰动强度，增加受污染水体中的DO浓度，达到降低水体的NH3-N指标浓度的效果，实现水质的改善。
5)底泥对TN的吸附与释放规律
H点、S点、M点和X点四种泥样的上覆水TN浓度和释放速率变化如图10.6所示。TN在底泥沉积物和上覆水中的迁移转化是一个动态过程，底泥间隙水和上覆水之间的TN浓度差有一定关系，当底泥间隙水的TN浓度高于上覆水时，表现出底泥向上覆水的释放过程，且浓度差越大，TN的释放量越大，但当上覆水中的TN浓度逐渐上升至一定数值时，又会发生底泥对上覆水中TN的吸附过程，这与以往的相关研究成果一致的。由于H点、S点、M点和X点四个点位的底泥和河道原水的TN浓度不同，其TN浓度和释放速率随时间变化规律不一致，但随扰动强度的变化幅度均不大，统计分析发现，各扰动工况之间的TN释放强度不具有显著性差异(P>0.05)，因此，TN随扰动强度变化不敏感。
步骤2，为了确定不同转速下试验装置中的流场分布特征，采用FLUENT软件进行模拟计算。根据圆筒中的水深和直径构建圆柱形三维模型(直径27cm，高40cm)，距离圆柱顶部以下5cm距离的中心位置布置十字形旋浆，旋浆的尺寸根据恒速电动搅拌器搅拌棒和底部十字形旋叶的实际尺寸设置，构建的模型如图4所示。采用非结构化网格划分方法对模型进行网格划分，划分的网格数量为287269，网格分布如图5所示。
本实例的研究对象是复杂的高速动边界，在桨叶周边区域存在高雷诺数的复杂湍流，为此选择标准k-ε模型作为湍流模型多参考坐标系处理桨叶动边界，根据相应的搅拌器转速，在桨叶边界及周边区域设置旋转参考坐标系，采用VOF模型处理自由表面边界。将搅拌器转速r分别设置为100rpm、200rpm、300rpm、400rpm、500rpm五种工况分别进行数值模拟，各工况从旋桨下缘至圆筒底部的纵剖面流速分布如图6所示。计算结果显示，装置中心轴线位置流速较小，靠近壁面处流速较大，从装置轴线至壁面水流速度逐渐增大；由于装置中水流的螺旋运动，靠近底部中心轴线位置一定水深的水流流速接近于0；随着转速的增大，装置中水流的最大流速逐渐增大。
步骤3，试验装置中搅拌器促使水体旋转会产生离心力，易造成水体震荡，导致同一水平面的切应力不等，这种现象的形成也很显然，为此，为尽可能地规避旋转水体带来的切应力不均的问题，本文选用同一水平面的切应力平均值作为每种扰动工况最终的切应力计算结果，具体计算方法如下。
分别统计分析不同转速条件，距试验装置轴线0、1cm、2cm～13cm位置的速度垂向分布，例如，距轴线12cm位置的速度分布如图7所示。流体的切应力代表速度梯度的大小，由图7可见，不同转速条件，均为距泥-水交界面深度1cm以下的速度变化最大，1cm以上的速度基本恒定，因此，取距泥-水交界面深度1cm水深范围内、距试验装置中轴线不同位置的切应力的平均值作为每个扰动工况的切应力，切应力计算公式如下：
式中，τ为切应力，N/m2；u为流速，m/s；y为水深，m；μ为水的动力粘度，试验过程中的水温基本维持在15℃左右，因此这里取15℃时水的动力粘度1.197*10-3Pa·s。
不同转速的切应力如表1所示。
表1不同转速泥水界面切应力统计表
步骤3，将该切应力作为中间变量，采用对数流速分布公式计算实际不同水深河道的表面流速，建立室内试验中的扰动转速与野外天然河道表面流速之间的相关关系。
H点、S点、M点和X点四处底泥采样点位所在河道的河道平均水深分别约为2.5m、2.5m、2.0m和2.0m，利用对数流速分布公式即可计算水深分别为2.0m、2.5m的河道表面流速，如图8所示。
步骤4，基于图8中扰动转速与野外天然河道表面流速之间的相关关系和试验中上覆水水质浓度变化，得到H点、S点、M点和X点的上覆水浊度、DO、TP和NH3-N四项敏感水质指标随流速的变化情况，如图9所示。
在上述四项水质指标中，DO浓度随水流流速的增大而增加，NH3-N释放强度随水流流速的增大而减小，该两项指标的变化情况说明，增大流速对河道水质改善有益，一定程度上解释了“流水不腐”的科学道理，但浊度和TP释放强度的变化趋势则不同，均随水流流速的增大而增加，因此，从浊度和TP指标来说，增大流速会促进底泥悬浮，并向上覆水释放TP，对河道水环境带来不利影响，且H点、X点、M点和S点四个点位分析结果均发现u＝0.13～0.15m/s为临界流速范围，当u<0.13m/s时，浊度和TP释放的变化量不大，增幅不明显，可以忽略，但是当u>0.13～0.15m/s时，底泥快速悬浮，浊度显著增大，总磷释放量显著增加。由此可见，水动力调控中抑制底泥快速释放的流速上限阈值u大的范围为0.13～0.15m/s。
步骤5，按照水环境容量计算的理论，根据水质本底值与目标水质，确定水动力调控改善水环境中的流速下限阈值。
以苏州市平江河、临顿河为例，以NH3-N浓度达到IV类标准(Cs＝1.5mg/L)为目标，计算分析其流速下限阈值。
如图2，选择总体达标法计算水环境容量，考虑点源污染、面源污染、直接入河的粉尘、底泥污染物的释放、河道水体的自净、水中植物对污染物的吸收等多种影响河道水质的因素。由污染物的质量守恒，得到以下公式：
步骤6，假设河道进口断面水流流速为u，进口和出口断面的水深分别为h1和h2，则：
化简得：
上式中，若水环境容量W＝0，即该河道水体不再能够承受污染物的排放，此时的流速u为能够保证河道水质的流速下限阈值u小，即：
平江河河道长度L＝2850m，平均河宽B＝9m，进出口断面水深取平江河的平均水深2.5m；引水浓度取西塘河2014-2017年平均NH3-N浓度(西塘河为苏州古城区的主要引水水源)，即C上＝0.85mg/L；平江河的点源污染由苏州市水利局提供，经调查，平江河沿线排放的NH3-N总浓度C点＝7.68mg/L，污水排放流量q点＝0.1m3/s；rs取前文四组底泥释放试验中NH3-N释放率最大值的平均值，为150mg/(m2*d)；K取文献中太湖流域NH3-N降解系数率定结果的平均值0.085d-1；忽略粉尘和面源污染，且无植物净化作用。
将上述相关数据代入u小计算公式得，平江河的流速下限阈值u小＝0.04m/s，即在该条件下，当平江河流速高于0.04m/s时，才能保证该河道的NH3-N浓度达到IV类标准。
2)临顿河
临顿河河段长度L＝1600m，河宽B＝8m，水深h＝2.5m，调查污水浓度C点＝5.28mg/L，引水浓度C上＝0.85mg/L，rs取值与前文相同，150mg/(m2*d)，同样忽略粉尘和面源污染，且无植物净化作用。
根据步骤6中的公式得，临顿河的流速下限阈值u小＝0.03m/s。该条件下，当临顿河流速高于0.03m/s时，可保证该河道的NH3-N浓度达到IV类标准。
如上所述，尽管参照特定的优选实施例已经表示和表述了本发明，但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下，可对其在形式上和细节上做出各种变化。
Claims (6)
1.一种面向平原城市河网水环境改善的水动力调控阈值确定方法，其特征是包括以下步骤：
步骤1、在城市河道内采集新鲜底泥和河水水样；
步骤2、建立数学模型；
步骤3、将切应力作为中间变量，建立室内试验中的扰动转速与野外天然河道表面流速之间的相关关系；
步骤4、获取实际河道表面流速与上覆水水质指标浓度间的响应关系：
基于扰动转速与野外天然河道表面流速之间的相关关系和试验中上覆水水质浓度变化，得到实际河道表面流速与上覆水水质指标浓度间的响应关系，以水质变化拐点处的流速作为平原城市河网水动力调控的流速上限阈值；
步骤5、按照水环境容量计算的理论，根据水质本底值与目标水质，确定水动力调控改善水环境中的流速下限阈值；
步骤6、在步骤5的基础上代入河道进口断面水流流速值、进口和出口断面的水深值，计算得出能够保证河道水质的流速下限阈值。
2.根据权利要求1所述的一种面向平原城市河网水环境改善的水动力调控阈值确定方法，其特征在于，步骤1进一步包括：
将采样得到的样本避光密封保存并带回实验室进行试验，模拟分析不同水动力条件下的底泥释放与吸附规律，根据试验结果，得到促使底泥快速释放的扰动转速。
3.根据权利要求1所述的一种面向平原城市河网水环境改善的水动力调控阈值确定方法，其特征在于，步骤2进一步包括：
计算不同扰动强度条件下圆筒试验装置中的流场，提取泥水界面的流速数据，采用切应力计算公式，计算得到装置中泥-水界面的切应力：
式中，
为试验装置中泥-水界面的切应力，N/m2；
为试验装置中的水流流速，m/s；y为水深，m；μ为水的动力粘度。
4.根据权利要求1所述的一种面向平原城市河网水环境改善的水动力调控阈值确定方法，其特征在于，步骤3进一步包括：
采用对数流速分布公式计算实际不同水深河道的表面流速：
式中，
为河道的表面流速，m/s；u *为摩阻流速，m/s；y为水深，m；Δ为绝对粗糙度，mm；ρ为水的密度；
为河道底部切应力，N/m2。
5.根据权利要求1所述的一种面向平原城市河网水环境改善的水动力调控阈值确定方法，其特征在于，步骤5进一步包括：
选择总体达标法计算水环境容量，总体达标计算法采用零维模型进行水质计算，考虑点源污染、面源污染、直接入河的粉尘、底泥污染物的释放、河道水体的自净、水中植物对污染物的吸收的因素，由污染物的质量守恒，得到以下公式：
式中，Q 上、C 上分别为河道上游来水的流量m3/s以及水质浓度mg/L；W为水环境容量mg/s；q 点、C 点分别为入河点源污染的流量m3/s和污染物浓度mg/L；q 面、C 面分别为入河面源污染的流量m3/s和污染物浓度mg/L；f 粉为直接入河粉尘所含污染物的质量函数；r s为底泥污染物的释放速率mg/(m2*d)；L为河道长度m；B为河宽m；C s为河道出口断面的目标水质浓度mg/L；K为污染物降解系数1/d；V为河段内的水体体积m3；f 植为水生植物吸收的污染物的质量函数，与光照、温度、植物种类、种植密度以及水深有关。
6.根据权利要求5所述的一种面向平原城市河网水环境改善的水动力调控阈值确定方法，其特征在于，步骤6进一步包括：
设河道进口断面水流流速为
，进口和出口断面的水深分别为h 1和h 2，代入步骤5的公式中得：
上式化简得：
式中，若水环境容量W=0，即该河道水体不再能够承受污染物的排放，此时的河道进口断面水流流速
为能够保证河道水质的流速下限阈值u 小，即：
上式中，各符号含义同上。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010402856.7A CN111524035B (zh) | 2020-05-13 | 2020-05-13 | 一种面向平原城市河网水环境改善的水动力调控阈值确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010402856.7A CN111524035B (zh) | 2020-05-13 | 2020-05-13 | 一种面向平原城市河网水环境改善的水动力调控阈值确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111524035A CN111524035A (zh) | 2020-08-11 |
| CN111524035B true CN111524035B (zh) | 2021-05-14 |
Family
ID=71912267
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010402856.7A Active CN111524035B (zh) | 2020-05-13 | 2020-05-13 | 一种面向平原城市河网水环境改善的水动力调控阈值确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111524035B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112632865B (zh) * | 2020-12-16 | 2022-12-09 | 中建三局第一建设工程有限责任公司 | 一种内河水系最优流速的调控方法 |
| CN115014702B (zh) * | 2022-05-17 | 2024-06-14 | 长江水利委员会长江科学院 | 一种量化水生生物运动行为对水质生态修复效果扰动规律的装置及方法 |
| CN117236674B (zh) * | 2023-11-16 | 2024-02-02 | 水利部交通运输部国家能源局南京水利科学研究院 | 城市河网水动力精准调控和水环境提升方法及系统 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102663223B (zh) * | 2011-12-27 | 2016-05-25 | 中国科学院生态环境研究中心 | 复杂河网突发污染事件动态污染场计算方法 |
| CN106202163A (zh) * | 2016-06-24 | 2016-12-07 | 中国环境科学研究院 | 通江湖泊生态监测信息管理及预警系统 |
| CN107084765A (zh) * | 2017-03-30 | 2017-08-22 | 中国农业大学 | 一种渠道流量计算方法及装置 |
| CN107103188B (zh) * | 2017-04-11 | 2020-05-08 | 山东大学 | 一种城市河道生态流量阈值的确定方法 |
| CN107045050B (zh) * | 2017-04-17 | 2019-03-12 | 华北电力大学 | 一种突发污染事故下利用流速阈值的分级预警方法 |
| GB201717116D0 (en) * | 2017-10-18 | 2017-11-29 | Severn Trent Water Ltd | Water distribution network |
| CN107808054A (zh) * | 2017-10-27 | 2018-03-16 | 河海大学 | 一种基于野外观测的沉积物再悬浮模型构建方法 |
| CN108320261A (zh) * | 2018-01-18 | 2018-07-24 | 北控水务(中国)投资有限公司 | 一种城市河道断面动态水质目标确定方法 |
-
2020
- 2020-05-13 CN CN202010402856.7A patent/CN111524035B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111524035A (zh) | 2020-08-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111524035B (zh) | 一种面向平原城市河网水环境改善的水动力调控阈值确定方法 | |
| Shum et al. | Organic matter processing in continental shelf sediments—the subtidal pump revisited | |
| Zhou et al. | Quantitative relationship between biofilms components and emitter clogging under reclaimed water drip irrigation | |
| Narnoli et al. | Sludge blanket of UASB reactor: mathematical simulation | |
| CN104075872B (zh) | 一种模拟往复流作用下沉积物再悬浮的循环直水槽装置 | |
| Ding et al. | Characteristics of sediment resuspension in Lake Taihu, China: a wave flume study | |
| Dou et al. | Modeling the interaction of light and nutrients as factors driving lake eutrophication | |
| Drury et al. | Transport of 1‐μm latex particles in Pseudomonas aeruginosa biofilms | |
| Logan et al. | Potential for increased nutrient uptake by flocculating diatoms | |
| Liang et al. | Study on nutrient distribution and interaction with sediments in a macro-tidal estuary | |
| Wang et al. | Response of bacterial communities to variation in water quality and physicochemical conditions in a river-reservoir system | |
| Fang et al. | Biofilm growth on cohesive sediment deposits: Laboratory experiment and model validation | |
| Yu et al. | Modelling the impact of hydrodynamic turbulence on the competition between Microcystis and Chlorella for light | |
| Sun et al. | Quantifying the effects of submerged aquatic vegetation on internal loading in lake: A modeling study of the largest shallow lake in North China | |
| CN103043784B (zh) | 一种活性污泥污水处理的多维、多相、多过程耦合模拟方法 | |
| Cheng et al. | Transport mechanisms of contaminants released from fine sediment in rivers | |
| Dou et al. | Eutrophication model driven by light and nutrients condition change in sluice-controlled river reaches | |
| Yang et al. | Microorganisms and their metabolic activities affect seepage through porous media in groundwater artificial recharge systems: A review | |
| Gantzer et al. | Effect of long-term water velocity changes on streambed biofilm activity | |
| Dai et al. | Numerical investigation on Re-oxygenation efficiency of stepped overflow weir in urban stream | |
| CN106483254A (zh) | 河流中微污染物自净处理系统 | |
| Hansen et al. | Uptake of dissolved nickel by Elodea canadensis and epiphytes influenced by fluid flow conditions | |
| Khezri et al. | Determination of the effect of wind velocity and direction changes on turbidity removal in rectangular sedimentation tanks | |
| CN117313303A (zh) | 一种考虑水体不同运动状态的微塑料迁移轨迹预测方法 | |
| Wan et al. | Critical velocity in phosphorus exchange processes across the sediment-water interface |
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 |