CN110322140A - 一种基于mike21的感潮河段水环境容量计算方法 - Google Patents
一种基于mike21的感潮河段水环境容量计算方法 Download PDFInfo
- Publication number
- CN110322140A CN110322140A CN201910576346.9A CN201910576346A CN110322140A CN 110322140 A CN110322140 A CN 110322140A CN 201910576346 A CN201910576346 A CN 201910576346A CN 110322140 A CN110322140 A CN 110322140A
- Authority
- CN
- China
- Prior art keywords
- water
- sewage draining
- draining exit
- tidal reach
- environmental
- 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.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 73
- 230000007613 environmental effect Effects 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000004364 calculation method Methods 0.000 title claims abstract description 16
- 239000010865 sewage Substances 0.000 claims abstract description 32
- 238000003908 quality control method Methods 0.000 claims description 19
- 238000012544 monitoring process Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000012876 topography Methods 0.000 claims description 3
- 238000007599 discharging Methods 0.000 claims 1
- 230000002457 bidirectional effect Effects 0.000 abstract description 4
- 238000013178 mathematical model Methods 0.000 abstract description 3
- 238000007726 management method Methods 0.000 abstract description 2
- 239000013535 sea water Substances 0.000 abstract 1
- 239000003344 environmental pollutant Substances 0.000 description 2
- 231100000719 pollutant Toxicity 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于MIKE21的感潮河段水环境容量计算方法,通过MIKE21水环境数学模型的前期计算,得出每个排污口对每个控制断面的贡献度系数,从而将水环境容量计算演变成线性规划问题的求解,采用单纯形法即可计算出感潮河段的水环境容量。以解决现有方法计算量大,且计算结果不唯一,受人为主观因素影响较大等问题。涉及感潮河段水环境容量及限制排污总量控制等环境管理领域,可广泛应用于受潮汐作用的双向流水体环境容量的计算,甚至可应用于海域水环境容量计算。
Description
技术领域
本发明涉及一种基于MIKE21的感潮河段水环境容量计算方法,属于感潮河段水环境容量及限制排污总量控制等环境管理领域。
背景技术
单向河流一般是采用一维水质模型计算其水环境容量,将天然河道概化成顺直河道,将非稳态水流简化为稳态水流,认为支流汇入主流即达到完全混合状态,认为污染物浓度在断面横向方向上变化不大,对横向和垂向的污染物浓度梯度予以忽略,将沿河排污口进行概化,在确定模型参数取值后,由河道一维水质模型计算出水环境容量。双向流的流态与单向流有很大不同,其一维水质计算模型不适用于双向流水体,现有感潮河段水环境容量计算方法多为试算法,即通过建立和调校能反映河段“污染源~水质动态响应关系”的水质数学模型,然后根据环境功能分区设定一系列控制点的水质目标,对水质模型进行反解,从而可得出满足水质目标的允许排放量。该方法计算量大,且计算结果不唯一,受人为主观因素影响较大。
发明内容
本发明提供一种基于MIKE21的感潮河段水环境容量计算方法,以解决现有方法计算量大,且计算结果不唯一,受人为主观因素影响较大等问题。
为解决上述问题,拟采用这样一种基于MIKE21的感潮河段水环境容量计算方法,具体如下:
①前期基础资料收集及现场调研,对感潮河段周边建设项目进行全面调查,了解现状河段取退水情况,了解所有支流的水质情况,了解感潮河段排污口水质、排放量、排放方式情况,水闸按排污口进行考虑,另外需收集水下地形数据,及相关水位站和水文站的水位、流量数据;
②对感潮河段进行控制断面选取,以水质目标为主要控制参数来划分不同控制单元;
③制定现场水文水质采样方案,在感潮河段选取合适断面对水量、水位、流速、流向、水质参数进行逐时同步现状监测;
④对监测数据进行分析和整理,并对MIKE21水动力、水质参数进行率定与验证;
⑤选取合适水动力、水质边界进行感潮河段背景浓度计算,同时计算各排污口对各控制断面的贡献度系数;
⑥采样线性规划的单纯形法计算出感潮河段水环境容量。
上述步骤⑤的计算方法为
目标函数:
约束条件:
xj≥O(j=1,…,n) (3)
式中,i为水质控制点编号;m为水质控制点数目;j为排污口编号;n为排污口数目;x为负荷量;L为总负荷量,即环境容量;aij为第j个排污口的单位负荷量对第i个水质控制点的污染贡献度系数;cb为水质控制点的污染背景浓度;为水质控制点处的环境标准值。
贡献度系数aij(i=1,…,m;j=1,…,n)的求解是以线性规划法求解环境容量的关键,也是工作量最大的部分,根据定义,假想在一个排污口给1个单位负荷量,即xj=1,而其它排污口无负荷量排出,即xk=0(k=1,…,n;k≠j),然后用二维水质模型计算出在这种情况下的浓度分布,确定出m个水质控制点的浓度值,即为aij(i=1,…,m),改变排污口,重复以上步骤,就可以求出每个排污口的aij,在以上所有计算中,均不考虑背景浓度。
与现有技术相比,本发明基于MIKE21水环境数学模型可较好地模拟出双向流水体中的水动力、水质因子随潮汐涨退的变化情况,由MIKE21模型计算出各排污口对各水质控制断面的污染贡献度系数,及水体背景浓度,从而将水环境容量的计算转化为线性规划矩阵的求解,采用单纯形法计算出该矩阵的最优解,即为感潮河段最大水环境容量。此方法避免人为主观因素对模型输入条件的干扰,同时避免了模型反演试算法的大尺度工作量。可较好的适用于感潮河段、潮汐河口、甚至是海湾水环境容量的计算。
附图说明
图1是本发明的技术线路图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细说明,应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
实施例
参照图1,本实施例提供一种基于MIKE21的感潮河段水环境容量计算方法,具体如下:
①前期基础资料收集及现场调研,对感潮河段周边建设项目进行全面调查,了解现状河段取退水情况,了解所有支流的水质情况,了解感潮河段排污口水质、排放量、排放方式情况,水闸按排污口进行考虑,另外需收集水下地形数据,及相关水位站和水文站的水位、流量数据;
②对感潮河段进行控制断面选取,以水质目标为主要控制参数来划分不同控制单元;
③制定现场水文水质采样方案,在感潮河段选取合适断面对水量、水位、流速、流向、水质参数进行逐时同步现状监测;
④对监测数据进行分析和整理,并对MIKE21水动力、水质参数进行率定与验证;
⑤选取合适水动力、水质边界进行感潮河段背景浓度计算,同时计算各排污口对各控制断面的贡献度系数;
⑥采样线性规划的单纯形法计算出感潮河段水环境容量。
步骤⑤的计算方法为
目标函数:
约束条件:
xj≥0(j=1,…,n) (3)
式中,i为水质控制点编号;m为水质控制点数目;j为排污口编号;n为排污口数目;x为负荷量;L为总负荷量,即环境容量;aij为第j个排污口的单位负荷量对第i个水质控制点的污染贡献度系数;cb为水质控制点的污染背景浓度;为水质控制点处的环境标准值。
贡献度系数aij(i=1,…,m;j=1,…,n)的求解是以线性规划法求解环境容量的关键,也是工作量最大的部分,根据定义,假想在一个排污口给1个单位负荷量,即xj=1,而其它排污口无负荷量排出,即xk=0(k=1,…,n;k≠j),然后用二维水质模型计算出在这种情况下的浓度分布,确定出m个水质控制点的浓度值,即为aij(i=1,…,m),改变排污口,重复以上步骤,就可以求出每个排污口的aij,在以上所有计算中,均不考虑背景浓度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于MIKE21的感潮河段水环境容量计算方法,其特征在于,具体如下:
①前期基础资料收集及现场调研,对感潮河段周边建设项目进行全面调查,了解现状河段取退水情况,了解所有支流的水质情况,了解感潮河段排污口水质、排放量、排放方式情况,水闸按排污口进行考虑,另外需收集水下地形数据,及相关水位站和水文站的水位、流量数据;
②对感潮河段进行控制断面选取,以水质目标为主要控制参数来划分不同控制单元;
③制定现场水文水质采样方案,在感潮河段选取合适断面对水量、水位、流速、流向、水质参数进行逐时同步现状监测;
④对监测数据进行分析和整理,并对MIKE21水动力、水质参数进行率定与验证;
⑤选取合适水动力、水质边界进行感潮河段背景浓度计算,同时计算各排污口对各控制断面的贡献度系数;
⑥采样线性规划的单纯形法计算出感潮河段水环境容量。
2.根据权利要求1所述一种基于MIKE21的感潮河段水环境容量计算方法,其特征在于,上述步骤⑤的计算方法为:
目标函数:
约束条件:
xj≥0(j=1,…,n) (3)
式中,i为水质控制点编号;m为水质控制点数目;j为排污口编号;n为排污口数目;x为负荷量;L为总负荷量,即环境容量;aij为第j个排污口的单位负荷量对第i个水质控制点的污染贡献度系数;cb为水质控制点的污染背景浓度;为水质控制点处的环境标准值。
根据定义,假想在一个排污口给1个单位负荷量,即xj=1,而其它排污口无负荷量排出,即xk=0(k=1,…,n;k≠j),然后用二维水质模型计算出在这种情况下的浓度分布,确定出m个水质控制点的浓度值,即为aij(i=1,…,m),改变排污口,重复以上步骤,求出每个排污口的aij。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910576346.9A CN110322140A (zh) | 2019-06-28 | 2019-06-28 | 一种基于mike21的感潮河段水环境容量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910576346.9A CN110322140A (zh) | 2019-06-28 | 2019-06-28 | 一种基于mike21的感潮河段水环境容量计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110322140A true CN110322140A (zh) | 2019-10-11 |
Family
ID=68121315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910576346.9A Pending CN110322140A (zh) | 2019-06-28 | 2019-06-28 | 一种基于mike21的感潮河段水环境容量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110322140A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539158A (zh) * | 2020-03-31 | 2020-08-14 | 安徽国祯环保节能科技股份有限公司 | 一种水体的最大污染负荷分配方法及系统 |
CN111814411A (zh) * | 2020-06-12 | 2020-10-23 | 长江水利委员会长江科学院 | 一种基于mike21和盲数理论的雨源型河流水环境容量计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218485A (zh) * | 2013-04-03 | 2013-07-24 | 郑州大学 | 一种gis技术支持下的小流域环境容量计算方法 |
CN108549744A (zh) * | 2018-03-14 | 2018-09-18 | 西安理工大学 | 一种河流水环境容量分配方法 |
CN109326332A (zh) * | 2018-10-30 | 2019-02-12 | 江苏省环境科学研究院 | 潮汐水网区水环境容量计算方法 |
CN109657418A (zh) * | 2019-01-31 | 2019-04-19 | 湖北省水利水电规划勘测设计院 | 一种基于mike21的湖泊水环境容量计算方法 |
-
2019
- 2019-06-28 CN CN201910576346.9A patent/CN110322140A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103218485A (zh) * | 2013-04-03 | 2013-07-24 | 郑州大学 | 一种gis技术支持下的小流域环境容量计算方法 |
CN108549744A (zh) * | 2018-03-14 | 2018-09-18 | 西安理工大学 | 一种河流水环境容量分配方法 |
CN109326332A (zh) * | 2018-10-30 | 2019-02-12 | 江苏省环境科学研究院 | 潮汐水网区水环境容量计算方法 |
CN109657418A (zh) * | 2019-01-31 | 2019-04-19 | 湖北省水利水电规划勘测设计院 | 一种基于mike21的湖泊水环境容量计算方法 |
Non-Patent Citations (1)
Title |
---|
万军明: "感潮河段水环境容量总量控制方案的研究", 《环境科学与管理》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539158A (zh) * | 2020-03-31 | 2020-08-14 | 安徽国祯环保节能科技股份有限公司 | 一种水体的最大污染负荷分配方法及系统 |
CN111539158B (zh) * | 2020-03-31 | 2023-06-16 | 中节能国祯环保科技股份有限公司 | 一种水体的最大污染负荷分配方法及系统 |
CN111814411A (zh) * | 2020-06-12 | 2020-10-23 | 长江水利委员会长江科学院 | 一种基于mike21和盲数理论的雨源型河流水环境容量计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109657418B (zh) | 一种基于mike21的湖泊水环境容量计算方法 | |
CN109815608A (zh) | 一种浅水湖泊群水质水量水生态耦合调度分析方法 | |
CN108256177A (zh) | 一种河流水沙模型的参数优化方法和系统 | |
CN107657912A (zh) | 一种基于Mike11的河流水环境容量确定方法 | |
CN111598754B (zh) | 基于水质模型的河网区水质断面污染贡献率核算方法 | |
CN113626923B (zh) | 一种平原河网区域水资源调度精准率计算与提升方法 | |
CN106777987B (zh) | 基于水动力、水质和som神经网络的河流水量聚类及水质评价方法 | |
Hayder et al. | Implementation of machine learning methods for monitoring and predicting water quality parameters | |
CN102662039A (zh) | 一种基于bp神经网络的水体溶解氧饱和度预测方法 | |
CN110570031A (zh) | 一种预测平原河网两级分汊河道分流比的方法 | |
CN109598428A (zh) | 一种基于行政单元和水系的污染物削减分配方法 | |
CN110322140A (zh) | 一种基于mike21的感潮河段水环境容量计算方法 | |
Fovet et al. | Modelling periphyton in irrigation canals | |
CN109726891A (zh) | 一种河流灰水足迹评估方法及其应用 | |
Samandar | A model of adaptive neural-based fuzzy inference system (ANFIS) for prediction of friction coefficient in open channel flow | |
Yao | Application of the ecohydraulic model on hydraulic and water resources engineering | |
Wang et al. | Integrated simulation and assessment of water quantity and quality for a river under changing environmental conditions | |
CN115422850B (zh) | 基于efdc和边缘计算判断河流入河排污口污染贡献率的方法 | |
Li et al. | Effects of the Three Gorges Project on the environment of Poyang Lake | |
CN116432439A (zh) | 一种基于数值模拟的城市河道纳污能力规划方法及系统 | |
Zhang et al. | Water exchange and pollutant diffusion law in Gangnan reservoir | |
CN100397405C (zh) | 基于变分算法的陆源排污量反演方法 | |
Yu et al. | Experimental study of the impact of the floating-vegetation island on mean and turbulence structure | |
CN113343601A (zh) | 一种复杂水系湖泊水位和污染物迁移的动态模拟方法 | |
Liu et al. | Parameter identification of river water quality models using a genetic algorithm |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191011 |
|
RJ01 | Rejection of invention patent application after publication |