CN104361536B - 一种渠道突发水污染事故限制扩散的调控方法 - Google Patents
一种渠道突发水污染事故限制扩散的调控方法 Download PDFInfo
- Publication number
- CN104361536B CN104361536B CN201410704332.8A CN201410704332A CN104361536B CN 104361536 B CN104361536 B CN 104361536B CN 201410704332 A CN201410704332 A CN 201410704332A CN 104361536 B CN104361536 B CN 104361536B
- Authority
- CN
- China
- Prior art keywords
- formula
- gate
- time
- close
- check gate
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000003911 water pollution Methods 0.000 title claims abstract description 19
- 238000009792 diffusion process Methods 0.000 title claims abstract description 11
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 28
- 231100000719 pollutant Toxicity 0.000 claims abstract description 28
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 230000004888 barrier function Effects 0.000 claims description 10
- 238000000205 computational method Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000011109 contamination Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 5
- 230000000007 visual effect Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000001925 catabolic effect Effects 0.000 abstract description 2
- 238000011065 in-situ storage Methods 0.000 abstract description 2
- 108090000623 proteins and genes Proteins 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 3
- 230000006378 damage Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035755 proliferation Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 230000009182 swimming Effects 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
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
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)
- Barrages (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明公开了一种渠道突发水污染事故限制扩散的调控方法,在污染事故发生后可通过关闭下游节制闸,将污染物控制在闸内,通过计算下游各节制闸临界关闸时间,可判断关闭下游节制闸的具体位置,解决了无法确定关闭下游具体节制闸及何时关闭的问题,具有物理概念清晰、计算量小、结果直观等优点,可为渠道污染事故原位处置技术提供技术支撑。
Description
技术领域
本发明涉及突发水污染事件应急处置领域,具体地指一种渠道突发水污染事故限制扩散的调控方法。
背景技术
引水渠道跨越长度大,干渠交叉建筑物众多,易发生突发水污染事件。突发水污染事故不同于一般的水污染,它没有固定的排放方式和排放途径,是突然发生,瞬时或短时间内大量地排放污染物质,对环境造成严重污染和破坏。其特点主要表现在事件发生的不确定性强、危害程度高、管控难度大、处置时间紧等。突发水污染事故处置不当或不及时,将会造成环境污染、生态破坏等危害,甚至发生停水、人员伤亡等。
目前所见关于突发污染事故的专利多为处置方法、模型预测、预警系统等方面的专利,未见对渠道内污染物限制扩散方面的专利,如一种基于云技术的突发性水环境风险预测系统及方法(CN104050388A)、基于物联网的大型水域突发水污染应急处置系统和方法(CN103236020A)、基于水动力约束原理的突发水污染生物预警装置及方法(CN103149336A)、一种廊道式水质污染预警装置及预警方法(CN103105398A)、突发水污染事件监控预警系统(CN102855526A)等。
实际上突发水污染事故发生后,为了方便快速处置、减少污染影响的范围和程度,首先应限制污染物影响的范围,尽量减少受污染水体总量。对渠道而言,可采用限制扩散的调控技术,即在一定时间内关闭下游节制闸将污染物限制在闸内,然后再进行集中处置。
发明内容
本发明的目的就是要解决上述背景技术的不足,提供一种渠道突发水污染事故限制扩散的调控方法,以使在突发水污染事故发生后,可快速处置、减少污染影响范围和程度。
本发明的技术方案为:一种渠道突发水污染事故限制扩散的调控方法,它包括以下步骤:
1)计算污染物泄漏点下游节制闸最晚关闸时间,该最晚关闸时间即临界关闸时间,在所述临界关闸时间内关闭下游节制闸以将污染物限制在下游节制闸内,所述临界关闸时间的计算方法为:
1.1)计算污染物泄漏至开始闭闸时间间隔t(min)内污染团峰值输移距离Dv:
Dv=60vst 式1)
式1中,vs为时间间隔t(min)内污染团峰值输移速度;
1.2)开始闭闸至水体静止污染团峰值输移距离DR包扩两部分,即节制闸完全关闭的时间tclose tclose(min)内水体输移距离DM DM和水流往复污染团输移距离DF DF;
DM=0.5×60×vstclose 式2)
DF=Catclose+Cb 式3)
DR=DM+DF=0.5×60×vstclose+Catclose+Cb 式4)
式中:vs为关闸前水流速度(m/s);
Ca和Cb分别为拟合参数;
1.3)在T(T=t+tclose)时间内污染团纵向长度W表示为:
式5)
式6)
式中:v为污染团纵向拉伸速度,Cc为纵向拉伸系数;
1.4)计算临界关闸时间:
Dv+DM+DF+0.5W<L 式7)
式7)中:L为污染物泄漏点至下游节制闸的距离;
将式1)、式2)、式3)、式6)代入式7)即可得到一个关于t的非线性方程:
式8)
利用式8)即可求得污染物泄漏点至开始闭闸的时间间隔t(min),即为临界关闸时间;
2)利用步骤1)计算污染物泄漏点下游第一个节制闸临界关闸时间,若第一个节制闸临界关闸时间小于实际操作时间,则关闭该节制闸将无法限制污染物在该节制闸内,需计算下游第二个节制闸临界关闸时间;
3)以下节制闸依次类推,直至临界关闸时间大于实际时间,该节制闸即为需要关闭的节制闸。
其中:式8)中污染物泄漏点至开始闭闸的时间间隔(t)使用迭代法进行求解的方法为:
令函数f(t)为:
则函数的导数f'(t)为:
当|tn+1-tn|小于某一微小变量,tn+1即为临界关闸时间。
本发明的一种渠道突发水污染事故限制扩散的调控方法,解决了无法确定关闭下游具体节制闸及何时关闭的问题,具有物理概念清晰、计算量小、结果直观等优点,可为渠道污染事故原位处置技术提供技术支撑。
具体实施方式
一种渠道突发水污染事故限制扩散的调控方法在intel fortran11.1上开发,运行环境为Microsoft visual studio 2012,详细流程如下,
1、首先输入基本参数,包括节制闸完全关闭时间tclose、纵向拉伸系数Cc等参数;
2、计算污染物泄漏点下游第一个节制闸临界关闸时间,计算方法为:
2.1)计算污染物泄漏至第一个节制闸开始关闭时的时间间隔t(min)内污染团峰值输移距离Dv:
Dv=60vst 式1)
式1中,vs为时间间隔t(min)内污染团峰值输移速度,为已知量;
2.2)第一个节制闸开始关闭至水体静止污染团峰值输移距离DR包扩两部分,即第一个节制闸完全关闭的时间tclose(min)内水体输移距离DM和水流往复污染团输移距离DF;
DM=0.5×60×vstclose 式2)
DF=Catclose+Cb 式3)
DR=DM+DF=0.5×60×vstclose+Catclose+Cb 式4)
式中:vs为关闸前水流速度(m/s),为已知量;
Ca和Cb分别为拟合参数,可根据水动力模拟结果拟合求出或者类比相关工程给出;
tclose为第一个节制闸完全关闭时间,可根据节制闸的特性进行设置;
1.3)在T(T=t+tclose)时间内污染团纵向长度W表示为:
式5)
式6)
式中:v为污染团纵向拉伸速度,可参考《南水北调中线明渠段事故污染特性模拟方法研究》;
Cc为纵向拉伸系数,可根据文献或者类比相关工程给出;
1.4)计算临界关闸时间:
Dv+DM+DF+0.5W<L 式7)
式7)中:L即为污染物泄漏点至第一个节制闸的距离,可实际测量得出;
将式1)、式2)、式3)、式6)代入式7)即可得到一个关于t的非线性方程:
式8)
利用式8)即可求得污染物泄漏点至第一个节制闸开始关闭的时间间隔t,该时间间隔t即为污染物泄漏点下游第一个节制闸临界关闸时间;
3、第一个节制闸临界关闸时间小于实际操作时间,则关闭该节制闸将无法限制污染物在第一个节制闸内,需计算下游第二个节制闸临界关闸时间;
4、以下节制闸依次类推,直至临界关闸时间大于实际时间,该节制闸即为关闭的节制闸;
5、最后输出闸址信息及其临界关闸时间。
式8)中污染物泄漏点至开始闭闸的时间间隔(t)可使用迭代法进行求解,其求解方法为:
令函数f(t)为:
则函数的导数f'(t)f'(t)为:
当|tn+1-tn|小于某一微小变量,tn+1即为临界关闸时间。
临界关闸时间控制方程f(t)为非线性方程,方程含有多解,使用迭代法进行求解可能会出现不符合物理现象的解,此时初始值的选取尤为重要。为确保解的合理性,临界关闸时间可表示为污染物泄漏点至下游节制闸距离L与水体流速vs的比值。
本发明在污染事故发生后可通过关闭下游节制闸,将污染物控制在闸内,通过计算下游各节制闸临界关闸时间,可判断关闭下游节制闸的具体位置,解决了无法确定关闭下游具体节制闸及何时关闭的问题,具有物理概念清晰、计算量小、结果直观等优点,可为渠道污染事故原位处置技术提供技术支撑。
Claims (2)
1.一种渠道突发水污染事故限制扩散的调控方法,其特征在于,它包括以下步骤:
1)计算污染物泄漏点下游节制闸最晚关闸时间,该最晚关闸时间即临界关闸时间,在所述临界关闸时间内关闭下游节制闸以将污染物限制在下游节制闸内,所述临界关闸时间的计算方法为:
1.1)计算污染物泄漏至开始闭闸时间间隔t(min)内污染团峰值输移距离Dv:
Dv=60vst 式1)
式1中,vs为时间间隔t(min)内污染团峰值输移速度;
1.2)开始闭闸至水体静止污染团峰值输移距离DRDR包括 两部分,即节制闸完全关闭的时间tclose(min)内水体输移距离DM和水流往复污染团输移距离DF;
DM=0.5×60×vstclose 式2)
DF=Catclose+Cb 式3)
DR=DM+DF=0.5×60×vstclose+Catclose+Cb 式4)
式中:vs为关闸前水流速度(m/s);
Ca和Cb分别为拟合参数;
1.3)在T(t=t+tclose)时间内污染团纵向长度W W表示为:
式5)
式6)
式中:v为污染团纵向拉伸速度,Cc为纵向拉伸系数;
1.4)计算临界关闸时间:
Dv+DM+DF+0.5W<L 式7)
式7)中:L为污染物泄漏点至下游节制闸的距离;
将式1)、式2)、式3)、式6)代入式7)即可得到一个关于t的非线性方程:
式8)
利用式8)即可求得污染物泄漏点至开始闭闸的时间间隔(t),即为临界关闸时间;
2)利用步骤1)计算污染物泄漏点下游第一个节制闸临界关闸时间,若第一个节制闸临界关闸时间小于实际操作时间,则关闭该节制闸将无法限制污染物在该节制闸内,需计算下游第二个节制闸临界关闸时间;
3)以下节制闸依次类推,直至临界关闸时间大于实际时间,该节制闸即为需要关闭的节制闸。
2.根据权利要求1所述的渠道突发水污染事故限制扩散的调控方法,其特征在于,式8)中污染物泄漏点至开始闭闸的时间间隔(t)使用迭代法进行求解的方法为:
令函数f(t)为:
则函数的导数f'(t)为:
当|tn+1-tn|小于某一微小变量,tn+1即为临界关闸时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410704332.8A CN104361536B (zh) | 2014-11-28 | 2014-11-28 | 一种渠道突发水污染事故限制扩散的调控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410704332.8A CN104361536B (zh) | 2014-11-28 | 2014-11-28 | 一种渠道突发水污染事故限制扩散的调控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104361536A CN104361536A (zh) | 2015-02-18 |
CN104361536B true CN104361536B (zh) | 2017-11-03 |
Family
ID=52528794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410704332.8A Active CN104361536B (zh) | 2014-11-28 | 2014-11-28 | 一种渠道突发水污染事故限制扩散的调控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104361536B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105243202A (zh) * | 2015-09-25 | 2016-01-13 | 天津大学 | 长距离输水工程突发水污染快速调控决策方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2804978A1 (fr) * | 2000-01-27 | 2001-08-17 | Francois Faraud | Barrage antipollution par confinement des nappes de produits polluants flottant a la surface de l'eau |
CN101853436A (zh) * | 2010-04-30 | 2010-10-06 | 中国环境科学研究院 | 一种水污染事故风险源识别方法 |
CN102129614A (zh) * | 2011-02-28 | 2011-07-20 | 中国环境科学研究院 | 一种饮用水源地突发污染事故风险诊断及分级方法 |
-
2014
- 2014-11-28 CN CN201410704332.8A patent/CN104361536B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2804978A1 (fr) * | 2000-01-27 | 2001-08-17 | Francois Faraud | Barrage antipollution par confinement des nappes de produits polluants flottant a la surface de l'eau |
CN101853436A (zh) * | 2010-04-30 | 2010-10-06 | 中国环境科学研究院 | 一种水污染事故风险源识别方法 |
CN102129614A (zh) * | 2011-02-28 | 2011-07-20 | 中国环境科学研究院 | 一种饮用水源地突发污染事故风险诊断及分级方法 |
Non-Patent Citations (1)
Title |
---|
突发水污染事故下明渠输水工程应急调控研究;刘婵玉;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20130515(第5期);正文第1-52页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104361536A (zh) | 2015-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105468844A (zh) | 管道内水-气耦合瞬变流的模拟方法 | |
Chosie et al. | Experimental and numerical investigation on the motion of discrete air pockets in pressurized water flows | |
Ardalan et al. | CFD and experimental study of 45° inclined thermal-saline reversible buoyant jets in stationary ambient | |
CN112036043B (zh) | 长距离重力流大管径供水管道破裂模型计算的方法 | |
Lowe | Sewer ventilation: factors affecting airflow and modeling approaches | |
CN104361536B (zh) | 一种渠道突发水污染事故限制扩散的调控方法 | |
Hamza et al. | Simulation of pollutant dispersion of a free surface flow in coastal water | |
Jiang et al. | Modelling the oil spill transport in inland waterways based on experimental study | |
Shibata et al. | Model evaluation of faecal contamination in coastal areas affected by urban rivers receiving combined sewer overflows | |
CN204941461U (zh) | 模拟解堵剂解堵的装置 | |
Elíasson et al. | Large hazardous floods as translatory waves | |
CN106777770A (zh) | 基于有限体积法的输水管道中空穴流的模拟方法 | |
Zolghadr et al. | Assessment of MIKE21 model in dam and dike-break simulation | |
Kho et al. | Quantitative dam break analysis on a reservoir earth dam | |
León et al. | Discussion of “Numerical oscillations in pipe-filling bore predictions by shock-capturing models” by JG Vasconcelos, SJ Wright, and PL Roe | |
Razad et al. | Simulation of breach outflow for earthfill dam | |
Ginting et al. | Numerical simulation of dam break using finite volume method case study of Situ Gintung | |
Kumar et al. | Effect of polymer additives on dynamics of water level in an open channel | |
Zhang et al. | Control optimization for liquid pipelines under leak condition | |
Dorin et al. | Mathematical modelling and numerical analysis of hydraulic system behaviour. A case study with application in HEC-RAS | |
Pokharel et al. | Inundation analysis of the effectiveness of an estuarygate in Hori River | |
da Silva et al. | LYSIS OF MULTIPHASE FLOW UNDER PARTIALLY OPEN FLOODGATES IN A SITUATION OF MAXIMUM FLOOD USING A NUMERICAL MODEL | |
Yao et al. | Estimation and Prediction Method of Water Level, Flow and Bursting Time of Baige Barrier Lake in Jinsha River | |
Panici et al. | A multi-phase modelling example of the formation and burst of landslide dams | |
Payesteh et al. | Interdependence of Hydraulic Parameters in Transient Induced Contaminant Intrusion in a Pipeline |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |