CN101759236A - 一种工业园区点源污染物的调控方法 - Google Patents
一种工业园区点源污染物的调控方法 Download PDFInfo
- Publication number
- CN101759236A CN101759236A CN200910264222A CN200910264222A CN101759236A CN 101759236 A CN101759236 A CN 101759236A CN 200910264222 A CN200910264222 A CN 200910264222A CN 200910264222 A CN200910264222 A CN 200910264222A CN 101759236 A CN101759236 A CN 101759236A
- Authority
- CN
- China
- Prior art keywords
- pollution factor
- regulation
- control
- sequence
- industrial park
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 12
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 12
- 230000001105 regulatory effect Effects 0.000 title abstract description 5
- 230000001276 controlling effect Effects 0.000 title abstract 2
- 238000009826 distribution Methods 0.000 claims abstract description 27
- 230000033228 biological regulation Effects 0.000 claims abstract description 23
- 239000002351 wastewater Substances 0.000 claims abstract description 20
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 238000007689 inspection Methods 0.000 claims abstract description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims abstract description 4
- 239000010786 composite waste Substances 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000007704 transition Effects 0.000 claims description 9
- 230000001186 cumulative effect Effects 0.000 claims description 7
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 5
- 238000007599 discharging Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000002844 melting Methods 0.000 claims description 3
- 230000008018 melting Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000009827 uniform distribution Methods 0.000 claims description 2
- 239000010865 sewage Substances 0.000 abstract description 8
- 238000012795 verification Methods 0.000 abstract 1
- 239000010842 industrial wastewater Substances 0.000 description 5
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 description 4
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 108010014173 Factor X Proteins 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000004043 dyeing Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 230000007096 poisonous effect Effects 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 238000011112 process operation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001988 toxicity Effects 0.000 description 1
- 231100000419 toxicity Toxicity 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000003827 upregulation Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Landscapes
- Sewage (AREA)
- Activated Sludge Processes (AREA)
Abstract
本发明公开了一种工业园区点源污染物的调控方法,属于污染控制领域。其步骤为:概率分布的确定;进行马氏链模拟,包括污染因子序列状态空间的统计、马氏链轨道的模拟和模拟序列的确定;预测模型的检验规则,采用双检验规则确定污染因子序列实测序列和模拟序列的相关系数;确定混合废水的调控方案,根据混合废水模型公式计算工业园区待调控污染因子的混合浓度,根据待调控污染因子的指标阈值判断该浓度是否超过调控要求。本发明将随机理论和废水混合原理应用到工业园区点源污染物的调控中,经过验证本方法具有较高的吻合度,通过建立的混合废水的调控方案进行调控,可为工业园区污水处理厂长效稳定建设运行、设施智能化控制和接管收费等提供指导。
Description
技术领域
本发明属于环境水污染物控制领域,具体的说是一种工业园区点源污染物的调控方法。
背景技术
工业园区包含了若干类不同性质的工业企业,这些企业占有一定面积的土地,且共同拥有一个对进入工业园区的企业提供必要的基础设施、服务、管理等的行政主管单位或公司。目前,我国已经将印染、化工、焦化、制药等对水环境有较大污染的企业集中到工业园区。这类园区内企业点源的工业废水,需经预处理达到接管标准后纳入工业园区污水处理厂集中处理。由于高浓度废水或在生产过程中使用了有毒有害化学品,点源废水水质水量不均衡、波动大、毒性高、盐度大,进入园区污水处理厂需要进行调控以稳定运行污水处理工艺。本发明可为工业园区集中污水处理厂长效稳定建设运行、设施智能化控制和接管收费等提供指导和借鉴。
发明内容
1.发明解决的技术问题
由于工业废水水质水量变化大,存在着较大的不确定性,容易导致集中污水处理工艺运行的不稳定。本发明的目的是提供一种工业园区点源工业废水污染物的调控方法。针对其特点,本发明从随机理论出发建立马氏链预测模型,预测工业园区点源工业废水主要污染因子的不确定性变化,建立工业园区混合工业废水的调控方案,从而达到控制点源污染物成分的目的。
2.技术方案
本发明的原理如下:
本发明基于随机原理,将点源工业废水污染因子作为随机变量处理,运用马尔可夫链蒙特卡洛的理论和方法研究预测其不确定性变化,基于废水混合原理建立调控方案,以达到发明的目的。
2.1马氏链预测模型的原理
(1)蒙特卡洛方法原理:
当所求解问题是某种随机事件出现的概率,通过随机数发生器重复多次地模拟该随机变量,以这种事件出现的频率估计这一随机事件的概率。
(2)马氏链原理:
如果时间参数集为离散集合(一般为正整数序列)且取值集合(状态空间)是离散集,称马尔可夫过程为马尔可夫链,简称马氏链,记为{Xt=X(t),t=0,1,2,K,N-1}。实际监测的污染因子数据即为一条马氏链。
马氏链预测将来的唯一有用的信息就是过程当前的状态,而与以前的状态无关。设污染因子Xt可以取m个离散状态,状态空间记为I={s1,s2,K,sm}。马氏链的性质完全由它的转移概率(转移核)来决定,用P(i,j)表示从状态si到状态sj的一步转移概率,意味着矩阵每一行元素的和为1,则一步状态转移概率矩阵为
Xt=si的概率记为状态概率Ei(t)。用E(t)表示在t时刻状态空间概率的行向量,Et=[E1(t),E2(t),K,Em(t)]。污染因子的初始向量用E(0)表示,通常E(0)中只有一个分量为1,其余全部是0,意味着马氏链从一个特定的状态开始,随着时间的变化,概率值扩散到整个状态空间。则基本方程可表示为:
E(t+1)=E(t)·P (1)
2.2废水混合模型
工业废水进入调蓄水池混合,可认为废水水质指标皆瞬时、完全混合,则有方程:
式中:Ck表示调蓄水池混合废水第k种污染因子浓度;Q1,...,Qm表示园区m家企业排放的废水流量;C1k,...,Cmk表示园区m家企业排放的第k种污染因子浓度。
本发明的技术方案如下:
一种工业园区点源污染物的调控方法,其步骤为:
(1)概率分布的确定:
将工业园区点源污染因子作为随机变量统计分析,采用频率直方图统计污染因子序列的分布规律,假设与其相似的概率分布,估计概率分布的参数,对于概率分布不明显的污染因子序列,可作为均匀分布处理;
(2)马氏链模拟
该部分是模拟污染因子不确定性变化的重要步骤,包括污染因子序列状态空间的统计、马氏链轨道的模拟和模拟序列的确定。
状态空间的统计:采用平行于时间轴且等间距直线的方法将污染因子序列划分为若干状态,状态的数目应根据污染因子序列的最大值、最小值及精度要求而定;根据此状态区间的划分,统计污染因子序列的初始概率,转移概率矩阵和累积概率矩阵;
初始概率分布是指污染因子序列在各个状态空间中的概率,即各状态空间的样本数量与所有样本数量的比值;转移概率是指污染因子序列从某一状态到另一状态的概率,计算方法为Pij=Mij/Mi。其中Mij为状态Ei转移到Ej的离散时刻个数;Mi为所有状态转移到Ei状态的离散个数。累积概率是指发生在某一状态区间的可能性有多大,计算这个状态区间内所有可能取值的概率之和。
(3)预测模型的检验规则
预测模型的检验规则决定着污染因子不确定性变化的模拟效果。本发明采用双检验规则(关联度检验和误差检验)确定污染因子序列实测序列和模拟序列的相关系数。
关联度检验:采用灰色关联分析作为关联度检验方法。灰色关联分析的基本思想是根据系统内部因素之间发展态势的相似、相异程度来衡量因素之间关联程度的一种方法。灰色相关度(grey relational degree,GRG)被定义为:
式中:Δ0i(k)是两个比较序列的绝对值,即Δ0i(k)=|x0(k)-xi(k)|;
(4)混合废水的调控方案
根据混合废水模型公式(2)计算工业园区待调控污染因子的混合浓度,根据待调控污染因子的指标阈值判断该浓度是否超过调控要求。
超过阈值的污染因子调控方案如下:
(1)确定超标的污染因子;
(2)计算超标最大的污染因子;
(3)查找废水混合前,排放该污染因子浓度最大的企业;
(4)根据混合废水模型,计算园区其他企业的混合废水浓度;
(5)若仍有超标的污染因子,返回(1);相反,则调控完毕。根据该调控模型得到的是待控制的排污企业,此时应将该企业排放的废水单独引入事故池,处理后再排入蓄水池。
3.有益效果
本发明提供了一种工业园区点源污染物的调控方法。将随机理论和废水混合原理应用到工业园区点源污染物的调控中,通过建立马氏链预测估算点源工业废水污染因子的不确定性变化,经过验证本方法具有较高的吻合度,通过建立的混合废水的调控方案进行调控以稳定运行污水处理工艺。本发明可为工业园区污水处理厂长效稳定建设运行、设施智能化控制和接管收费等提供指导。
附图说明
图1为本发明流程图;
图2为随机模拟步骤图
图3为废水COD指标的频率直方图
图4为污染因子原始序列与模拟序列对比图
具体实施方式
以下结合附图进一步说明本发明。
以某工业园区为研究对象,对园区企业点源废水水质进行监测采样。以园区3家企业、调控污染因子为COD指标为例说明本发明的应用过程。
1、污染因子的预测
(1)确定污染因子的概率分布
统计园区企业废水的污染因子的历史数据,其频率直方图见图3。
由图3可以看出污染因子的频率分布符合正态分布的特征,故假设其符合正态分布的规律,正态分布的分布函数为:
采用极大似然估计法估算正态分布的参数平均值为1121.6,方差为68490。由于柯尔莫哥洛夫检验法能够利用实测样本经验分布和推论总体分布之间的最大差异作为检验统计量。故本发明采用柯尔莫哥洛夫检验法,给定显著水平α=0.01,经检验DCOD<D31(0.01),故接受原假设,即确定该原始序列得总体分布为正态分布。
(2)马氏链的模拟
采用平行于时间轴且等间距直线的方法将污染因子序列分成5个状态。统计各状态的初始概率、转移概率矩阵、累积概率矩阵分别为:
a=|0.087,0.279,0.317,0.25,0.067|
根据其累积概率采用蒙特卡洛法随机抽样获得马氏链的一条轨道,再根据污染因子的概率分布获得每个状态序列的数值,重复多次即可获得模拟序列的集合。
(3)模拟序列的确定
设定关联度阈值为0.75,MAPE阈值为0.2,即模拟序列与原始序列的关联度系数大于0.75且MAPE小于0.2判定为有效。确定一条符合检验规则的模拟序列作为最终的模拟结果。污染因子序列与模拟序列的对比结果如图4。
由图4污染因子原始序列与模拟序列对比结果可知,本发明获得的模拟序列能够较好的反映点源工业废水污染因子的变化规律,为优化调控提供数据基础。同理可模拟园区其他企业污染因子的模拟序列。
2、混合废水的调控方案
通过马氏链预测模型可获得园区企业的废水流量和待调控污染因子的浓度,通过混合废水模型可计算得到污染因子的混合浓度计算结果见表1。
表1园区企业废水
设定污染因子的调控阈值上限为4500mg/l,从表1可知T=5和T=11时刻超过调控阈值,根据调控方案可得到该两个时刻企业C排放废水的污染因子浓度最高,将企业C于T=5和T=11时刻废水引入单独的事故池单独处理。调控后的污染因子浓度见表1最后一列。
Claims (2)
1.一种工业园区点源污染物的调控方法,其步骤为:
(1)概率分布的确定:
将工业园区点源污染因子作为随机变量统计分析,采用频率直方图统计污染因子序列的分布规律,假设与其相似的概率分布,估计概率分布的参数,对于概率分布不明显的污染因子序列,可作为均匀分布处理;
(2)进行马氏链模拟:
包括污染因子序列状态空间的统计、马氏链轨道的模拟和模拟序列的确定;
状态空间的统计:采用平行于时间轴且等间距直线的方法将污染因子序列划分为若干状态,状态的数目应根据污染因子序列的最大值、最小值及精度要求而定;根据此状态区间的划分,统计污染因子序列的初始概率,转移概率矩阵和累积概率矩阵;
马氏链轨道的模拟:运用蒙特卡洛方法产生均匀分布的一个随机数,根据污染因子序列状态空间的初始概率分布和累积概率分布矩阵来判定该数值所对应的状态,生成一个马氏链的轨道;
确定模拟序列:根据生成的马氏链轨道以及污染因子序列的概率分布,运用蒙特卡洛方法随机抽样获得每个状态空间的模拟值,从而可获得多个模拟序列;
(3)预测模型的检验规则:
采用双检验规则确定污染因子序列实测序列和模拟序列的相关系数;
(4)混合废水的调控方案
根据混合废水模型公式计算工业园区待调控污染因子的混合浓度,根据待调控污染因子的指标阈值判断该浓度是否超过调控要求。
2.根据权利要求1所述的工业园区点源污染物的调控方法,其特征在于对于步骤(4)中超过阈值的污染因子调控方案如下:
1)确定超标的污染因子;
2)计算超标最大的污染因子;
3)查找废水混合前,排放该污染因子浓度最大的企业;
4)根据混合废水模型,计算园区其他企业的混合废水浓度;
5)若仍有超标的污染因子,返回1);相反,则调控完毕。
根据该调控模型得到的是待控制的排污企业,此时应将该企业排放的废水单独引入事故池,处理后再排入蓄水池。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102642223A CN101759236B (zh) | 2009-12-31 | 2009-12-31 | 一种工业园区点源污染物的调控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102642223A CN101759236B (zh) | 2009-12-31 | 2009-12-31 | 一种工业园区点源污染物的调控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101759236A true CN101759236A (zh) | 2010-06-30 |
CN101759236B CN101759236B (zh) | 2012-06-27 |
Family
ID=42490643
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102642223A Expired - Fee Related CN101759236B (zh) | 2009-12-31 | 2009-12-31 | 一种工业园区点源污染物的调控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101759236B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254081A (zh) * | 2011-04-06 | 2011-11-23 | 北方工业大学 | 基于随机概率模型的突发河流水污染的统计定位方法 |
CN102306235A (zh) * | 2011-07-01 | 2012-01-04 | 北京航空航天大学 | 基于单边u检测的突发污染感知方法 |
CN103632304A (zh) * | 2013-10-27 | 2014-03-12 | 沈阳建筑大学 | 一种印染行业污水处理调控方法 |
CN103810537A (zh) * | 2014-02-12 | 2014-05-21 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
CN106548442A (zh) * | 2016-11-10 | 2017-03-29 | 湖南省环境保护科学研究院 | 一种大气污染物排放源清单精细化动态管理方法 |
CN107585890A (zh) * | 2017-09-22 | 2018-01-16 | 天津大学 | 基于蒙特卡罗模拟的新兴污染物污水优化处理方法 |
CN109190838A (zh) * | 2018-09-20 | 2019-01-11 | 信息产业电子第十设计研究院科技工程股份有限公司 | 基于工业园区废水集中处理预测系统 |
CN117819706A (zh) * | 2024-03-01 | 2024-04-05 | 广东新泰隆环保集团有限公司 | 一种用于农村污水处理的生物转盘污水处理系统 |
-
2009
- 2009-12-31 CN CN2009102642223A patent/CN101759236B/zh not_active Expired - Fee Related
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102254081A (zh) * | 2011-04-06 | 2011-11-23 | 北方工业大学 | 基于随机概率模型的突发河流水污染的统计定位方法 |
CN102254081B (zh) * | 2011-04-06 | 2014-11-05 | 北方工业大学 | 基于随机概率模型的突发河流水污染的统计定位方法 |
CN102306235A (zh) * | 2011-07-01 | 2012-01-04 | 北京航空航天大学 | 基于单边u检测的突发污染感知方法 |
CN103632304A (zh) * | 2013-10-27 | 2014-03-12 | 沈阳建筑大学 | 一种印染行业污水处理调控方法 |
CN103810537A (zh) * | 2014-02-12 | 2014-05-21 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
CN103810537B (zh) * | 2014-02-12 | 2016-09-28 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
CN106548442A (zh) * | 2016-11-10 | 2017-03-29 | 湖南省环境保护科学研究院 | 一种大气污染物排放源清单精细化动态管理方法 |
CN107585890A (zh) * | 2017-09-22 | 2018-01-16 | 天津大学 | 基于蒙特卡罗模拟的新兴污染物污水优化处理方法 |
CN109190838A (zh) * | 2018-09-20 | 2019-01-11 | 信息产业电子第十设计研究院科技工程股份有限公司 | 基于工业园区废水集中处理预测系统 |
CN117819706A (zh) * | 2024-03-01 | 2024-04-05 | 广东新泰隆环保集团有限公司 | 一种用于农村污水处理的生物转盘污水处理系统 |
CN117819706B (zh) * | 2024-03-01 | 2024-05-24 | 广东新泰隆环保集团有限公司 | 一种用于农村污水处理的生物转盘污水处理系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101759236B (zh) | 2012-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101759236B (zh) | 一种工业园区点源污染物的调控方法 | |
Zeng et al. | Optimization of wastewater treatment alternative selection by hierarchy grey relational analysis | |
Gillot et al. | Optimization of wastewater treatment plant design and operation using simulation and cost analysis | |
CN106295800A (zh) | 一种基于递归自组织rbf神经网络的出水总氮tn智能检测方法 | |
CN106600509B (zh) | 一种基于基础数据分析判断企业兑水排污行为的方法 | |
CN105574326A (zh) | 一种基于自组织模糊神经网络的出水氨氮浓度软测量方法 | |
CN110232256B (zh) | 一种基于kpls和rwfcm的污水处理过程监测方法 | |
Li et al. | Integration of life cycle assessment and statistical analysis to understand the influence of rainfall on WWTPs with combined sewer systems | |
CN103632032A (zh) | 一种城市污水处理过程出水指标在线软测量预测方法 | |
CN103886369A (zh) | 一种基于模糊神经网络的出水总磷tp预测方法 | |
CN106203756A (zh) | 一种基于特征数据库的工业污水污染物溯源分析方法 | |
Su et al. | Water–energy–carbon nexus: greenhouse gas emissions from integrated urban drainage systems in China | |
Zhao et al. | Seasonal artificial neural network model for water quality prediction via a clustering analysis method in a wastewater treatment plant of China | |
CN108088974A (zh) | 一种厌氧同时反硝化产甲烷过程出水硝氮的软测量方法 | |
Xu et al. | Prediction and evaluation of indirect carbon emission from electrical consumption in multiple full-scale wastewater treatment plants via automated machine learning-based analysis | |
Meng et al. | Regulatory implications of integrated real-time control technology under environmental uncertainty | |
CN102778548B (zh) | 一种污水处理过程污泥膨胀指数预测方法 | |
CN108549740A (zh) | 一种基于混合智能算法的厌氧系统出水氨氮软测量方法 | |
Canale et al. | Stochastic phosphorous model for Onondaga Lake | |
Yang et al. | Cloud-model-based feature engineering to analyze the energy–water nexus of a full-scale wastewater treatment plant | |
Rafati et al. | Determine the most effective process control parameters on activated sludge based on particle swarm optimisation algorithm (Case Study: South wastewater treatment plant of Tehran) | |
Crabtree et al. | A comparison of two river quality models | |
CN201330211Y (zh) | 污水处理厂运行参数自寻优模拟系统 | |
CN102787582B (zh) | 通过构建水资源模型来预警的方法 | |
Poch et al. | Where are we in wastewater treatment plants data management? A review and a proposal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120627 Termination date: 20211231 |
|
CF01 | Termination of patent right due to non-payment of annual fee |