当前位置: 首页 > 专利查询>祁超祯专利>正文

污染物衰减系数估算方法技术

技术编号:19009707 阅读:28 留言:0更新日期:2018-09-22 09:30
本发明专利技术公开了一种污染物衰减系数估算方法,包括确定稳态一维忽略污染物的沉降和弥散作用水质数学模型的基本方程;解析基本方程;检测河流数据和入河排污口数据;确定目标函数;目标函数变形;利用牛顿法迭代目标函数;判断是否终止迭代,如终止迭代,最后一次迭代值就是所求的最优值,如不终止迭代,继续迭代,直到出现最优值。本发明专利技术的有益效果为:提出了利用牛顿(切线)求极值法和计算机技术,解决了有排污口存在的河段污染物衰减系数的估算难题,为实际的水资源保护和河流管理起到指导作用。

Estimation method of pollutant attenuation coefficient

The invention discloses a method for estimating the attenuation coefficient of pollutants, including determining the basic equations of a steady-state one-dimensional water quality mathematical model neglecting the sedimentation and dispersion of pollutants; analysing the basic equations; detecting the river data and the sewage outlet data; determining the objective function; deforming the objective function; and iterating the objective function by Newton method. Number; determine whether to terminate the iteration, such as termination of the iteration, the last iteration value is the optimal value, if not termination of the iteration, continue to iterate until the optimal value. The invention has the beneficial effect that the extremum method of Newton (tangent) and the computer technology are used to solve the difficult problem of estimating the attenuation coefficient of pollutants in the river section with sewage outlets, and plays a guiding role in the actual water resources protection and river management.

【技术实现步骤摘要】
污染物衰减系数估算方法
本专利技术涉及水环境保护
,具体涉及一种污染物衰减系数估算方法。
技术介绍
污染物进入河流在输移过程中通过物理、化学及生物的作用发生浓度衰减,其衰减系数反映了污染物在水体作用下降解速度的快慢。河流中污染物的衰减系数是研究河流水质污染变化、计算水环境容量以及纳污能力的重要参数;在区域排污总量控制计划的制定、总量负荷指标的科学分配、控制计划执行过程的管理等工作中发挥了重要作用;为实现污染物排放浓度管理制度向总量管理制度的转变奠定了坚实的基础。同时,在水资源保护工作中,对计算水体纳污能力、预测污染物浓度,制定污染物控制方案以及制定河流水质管理规划方案等整个过程的准确性和合理性的影响显著。本文系统地概述了河流中污染物衰减系数的国内外研究现状,并对衰减系数的影响因素加以分析,以期为我国河流污染物衰减系数的研究工作提供科学依据。由于水体中污染物的种类繁多、性质复杂,水质标准中常用综合指标COD、BOD、氨氮等确定污染物含量。因此,水体中污染物衰减系数的研究重点主要是COD、BOD、氨氮的衰减。目前,无有效的污染物衰减系数确定方法,无法对水资源保护和河流管理起指导作用。
技术实现思路
本专利技术的目的在于提供一种污染物衰减系数估算方法,用以解决现有的不能确定污染物衰减系数从而导致的无法指导水资源保护和河流管理工作的问题。为实现上述目的,本专利技术的技术方案为一种污染物衰减系数估算方法,包括:S1、确定稳态一维忽略污染物的沉降和弥散作用水质数学模型的基本方程S2、稳定一维均匀河流水质模型的解析解为:c=c0e-kt(2)S3、检测河流数据和入河排污口数据;S4、确定目标函数S5、设y=e-k,则c=c0e-kt,变为c=c0yt,则目标函数变为:S6、利用牛顿法迭代公式,y2=y1-f′(y)/f″(y),式中f′(y)、f″(y)分别为:S7、当y2-y1≤0.0001时,并且根据给出的初值或迭代值和m次实测水质、水量资料,计算出f′(y)和f″(y),判别步骤S3的河流数据和入河排污口数据是否全部累加完,如果全部累加完,退出循环,计算完毕,最后一次迭代值就是所求的最优值;S8、把以上所求的最优值y,通过公式k=-lny转化成k;S9、当y2-y1>0.0001时,或者没有全部累加完,重复步骤S6和S7。进一步的,所述河流数据包括第一水文数据和水质数据,所述第一水文数据包括流量、流速、风向和风力数据,所述水质数据包括水温、溶解氧、化学需氧量、高锰酸盐指数和氨氮数据。进一步的,所述入河排污口数据包括排污口数据、第二水文数据和水质数据。进一步的,所述排污口数据为入河排污口位于河流的里程数。进一步的,所述第二水文数据为取样位置的流量。进一步的,所述第二水文数据包括流量、流速、风向和风力数据,所述水质数据包括水温、溶解氧、化学需氧量、高锰酸盐指数和氨氮数据。本专利技术具有如下优点:提出了利用牛顿(切线)求极值法和计算机技术,解决了有排污口存在的河段污染物衰减系数的估算难题,为实际的水资源保护和河流管理起到指导作用。具体实施方式以下实施例用于说明本专利技术,但不用来限制本专利技术的范围。本专利技术公开了一种污染物衰减系数估算方法,其包括:S1、确定稳态一维忽略污染物的沉降和弥散作用水质数学模型的基本方程S2、稳定一维均匀河流水质模型的解析解为:c=c0e-kt(2)S3、检测河流数据和入河排污口数据;S4、确定目标函数S5、设y=e-k,则c=c0e-kt,变为c=c0yt,则目标函数变为:S6、利用牛顿法迭代公式,y2=y1-f′(y)/f″(y),式中f′(y)、f″(y)分别为:S7、当y2-y1≤0.0001时,并且根据给出的初值或迭代值和m次实测水质、水量资料,计算出f′(y)和f″(y),判别步骤S3的河流数据和入河排污口数据是否全部累加完,如果全部累加完,退出循环,计算完毕,最后一次迭代值就是所求的最优值;S8、把以上所求的最优值y,通过公式k=-lny转化成k;S9、当y2-y1>0.0001时,或者没有全部累加完,重复步骤S6和S7。具体地,河流数据包括第一水文数据和水质数据,第一水文数据包括流量、流速、风向和风力数据,水质数据包括水温、溶解氧、化学需氧量、高锰酸盐指数和氨氮数据。入河排污口数据包括排污口数据、第二水文数据和水质数据。排污口数据为入河排污口位于河流的里程数。第二水文数据为取样位置的流量。第二水文数据包括流量、流速、风向和风力数据,水质数据包括水温、溶解氧、化学需氧量、高锰酸盐指数和氨氮数据。上述方法的确定过程如下:1水质模型的确定水环境容量计算的设计条件是水质最不利的情况。发生这种条件的河流流量,一般为一定保证率的低流量,这时的水流状态一般为稳定状态均匀流。污染物在中小河流迁移过程中,可忽略污染物的沉降和弥散作用。中小河流水环境容量计算,一般采用的水质模型是稳定状态下一维忽略污染物的沉降和弥散作用数学模型。稳态一维忽略污染物的沉降和弥散作用水质数学模型的基本方程:式中:x-上下断面距离(km),u-河段平均流速(km/d),c-污染物浓度(mg/l),k-污染物衰减系数(l/d)。2目标函数的建立2.1基本假设监测河段的水流水质状态和入河排污口水量水质变化稳定;污染物质在河流横断面分布均匀,符合一级反应动力学;忽略污染物的沉降、弥散作用;稳定一维均匀河流水质模型的解析解为:c=c0e-kt(2)式中:c-河段下断面污染物浓度(mg/l),c0-河段上断面污染物浓度(mg/l),k--污染物衰减系数(l/d),t-污染物自上断面到下断面的移动时间,d。2.2资料收集按照以上基本假设条件,布设水文水质监测断面,监测、收集资料。水质模型参数估算所需要资料如表1。表1水质模型参数K值估算所需要资料2.3目标函数的建立由实测值与模型计算值之差的平方,构成目标函数。即目标函数:式中:f(k)—目标函数;k—污染物衰减系数,1/d;i—水质指标实测数;n—排污口个数;Q上i---上断面第i次实测流量(m3/s);Q下i---下断面第i次实测流量(m3/s);C上i---上断面第i次污染浓度(mg/L);C下i---下断面第i次污染浓度(mg/L);A1i—第1个排污口第i次污染浓度(mg/L);q1i—第1个排污口第i次实测流量(m3/s);t1i—第1个排污口污染物移动至下断面的时间(d);Ani—第n个排污口第i次污染浓度(mg/L);qni—第n个排污口第i次实测流量(m3/s);ti—污染物从上断面移动至下断面的时间(d);tni—第n个排污口污染物质移动至下断面的时间(d);m—各水质指标的实测次数。3模型参数K值的估算3.1确定模型参数K值的估算公式在目标函数中,K是所需估算的参数,这是一个单变量函数求极值的问题。在无约束条件下,用牛顿(切线)法,求出目标函数中最小值时的参数。为了方便,设y=e-k,则c=c0e-kt,变为c=c0yt,则目标函数变为:式中符号含义同上。判断标准,一般二次迭代值之差小于或等于0.0001,迭代值接近于最优值。牛顿(切线)法迭代公式:y2=y1-f′(y)/f″(y),式中f′(y)、f″(y)分别为:式中:y1--迭代计算初值,y2—迭代计算终值。本文档来自技高网...

【技术保护点】
1.一种污染物衰减系数估算方法,其特征在于,包括:S1、确定稳态一维忽略污染物的沉降和弥散作用水质数学模型的基本方程

【技术特征摘要】
1.一种污染物衰减系数估算方法,其特征在于,包括:S1、确定稳态一维忽略污染物的沉降和弥散作用水质数学模型的基本方程S2、稳定一维均匀河流水质模型的解析解为:c=c0e-kt(2)S3、检测河流数据和入河排污口数据;S4、确定目标函数S5、设y=e-k,则c=c0e-kt,变为c=c0yt,则目标函数变为:S6、利用牛顿法迭代公式,y2=y1-f′(y)/f″(y),式中f′(y)、f″(y)分别为:S7、当y2-y1≤0.0001时,并且根据给出的初值或迭代值和m次实测水质、水量资料,计算出f′(y)和f″(y),判别步骤S3的河流数据和入河排污口数据是否全部累加完,如果全部累加完,退出循环,计算完毕,最后一次迭代值就是所求的最优值;S8、把以上所求的最优值y,通过公式k=-lny转化成k;S9、当y2-y1>0.0001时,或...

【专利技术属性】
技术研发人员:祁超祯张含明
申请(专利权)人:祁超祯
类型:发明
国别省市:山东,37

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1