实用新型内容
本实用新型的目的是提供一种水中重金属污染风险在线预警系统,对水中重金属污染进行在线监测和风险预警,实现水污染源监控的智能化和自动化。
为实现上述目的,本实用新型采取的技术方案是:
一种水中重金属污染风险在线预警系统,其包括在线检测装置、数据处理装置和终端设备,所述在线检测装置、数据处理装置和终端设备之间通过有线或无线网络连接;
在线检测装置,用于在线获取水质中的重金属浓度数据并传送至数据处理装置,其包括水样采集模块、在线消解模块、重金属富集模块以及光谱检测模块,水样采集模块按一定时间间隔连续自动采集水样,所述在线消解模块对水样进行消化处理,所述重金属富集模块对消化处理后的水样进行富集和洗脱,所述光谱检测模块对洗脱出来的水样进行重金属浓度检测,并将重金属浓度数据传送至数据处理装置;
数据处理装置,用于根据重金属浓度数据获得水质平均污染指数和偏离系数,并向终端设备发送预警信号,具体地,数据处理装置上设有水质评价模块和预警模块,水质评价模块根据重金属浓度数据获得平均污染指数,预警模块利用指数平滑法获得重金属浓度的预测值,并将实际重金属浓度与预测值相对比获得偏离系数,如果平均污染指数和偏离系数超出正常阈值,则向终端设备发送预警信号;
终端设备,用于接收和显示来自数据处理装置的数据和信号。
本实用新型通过自动采集水样的监测数据,根据预置的水质评价模块和预警模块判断水质状况,能减少人的工作量,节约时间,并通过平均指数模型可避免某个异常值达到阈值而导致的频发预警,在风险评价上更接近于真实水样;而且本实用新型能及时自动向用户发出预警信号,可实现无人值守监控。
作为本实用新型的一种改进,还包括数据库服务器装置,所述的数据库服务器装置用于存取水质实时监测数据及历史数据,仪器运行日志,包括已确认的每日、周、月的水质检测数据,以及水质状况和污染预警事件信息。
进一步地,所述在线检测装置为基于原子发射光谱技术的重金属在线监测设备。
进一步地,所述重金属包括Cd、Cu、Cr、Co、Fe、Mn、Ni、Ti、V、Pb、Zn和Bi。
一种水中重金属污染风险在线预警方法,其包括以下步骤;
通过在线检测获得时间t时水质中的重金属浓度数据;
通过下式计算单项污染指数以及平均污染指数,计算公式如下:
式中:K—平均污染指数;Ci—某金属浓度的实际测定值;C0—某金属浓度的水质标准限值;N—金属的种类数;
通过下式(1)计算出在时间t时某重金属的浓度平均值,
式中:Xt为在时间t时某重金属的浓度平均值,Xi为在时间t时第i个某重金属的浓度数据,n为某重金属浓度数据个数;
根据某重金属的浓度平均值,利用指数平滑法,计算时间t的一次指数平滑值,通过公式(2)计算:
St (1)=aXt+(1-a)St-1 (1) (2)
式中:St (1)为在时间t时某重金属浓度的一次指数平滑值;a为常数,取值范围为[0,1];
St-1 (1)为时间t-1时某重金属浓度的一次指数平滑值;
利用指数平滑法,计算时间t的二次指数平滑值,通过公式(3)计算:
St (2)=aSt (1)+(1-a)St-1 (2) (3)
St (2)为在时间t时某重金属浓度的二次指数平滑值,St-1 (2)为在时间t-1时某重金属浓度的二次指数平滑值;
利用指数平滑法,计算时间t的三次指数平滑值,通过公式(4)计算:
St (3)=aSt (2)+(1-a)St-1 (3) (4)
St (3)为在时间t时某重金属浓度的三次指数平滑值,St-1 (3)为在时间t-1时某重金属浓度的三次指数平滑值;
利用三次指数平滑法计算某重金属在时间t+1时的浓度预测值Pt+1,具体地:
Pt+1=at+bt+ct (5)
at=3St (1)-3St (2)+St (3) (6)
根据某重金属在时间t+1时的浓度预测值Pt+1,通过与在时间t+1时某重金属的实际值P相对比,计算出偏离系数b:
如果平均污染指数K和偏离系数均b超过设定阈值,则向用户发送预警信号;否则,不发送预警信号。
所述水质的扩散模型为一维水质扩散模型。
所述预警信号设有5个档次,按偏离系数b的大小分级,分别为一级预警、二级预警、三级预警、四级预警和五级预警。
与现有技术相比,本实用新型具有以下优点:
本实用新型通过自动采集水样的监测数据,根据预置的水质评价模块和预警模块判断水质中多种重金属参数状况,能减少人的工作量,节约时间,并通过平均指数模型可避免某个异常值达到阈值而导致的频发预警,在风险评价上更接近于真实水样;而且本实用新型能及时自动向用户发出预警信号,可实现无人值守监控。
实施例
请参考图1,一种水中重金属污染风险在线预警系统,其包括在线检测装置10、数据处理装置20和终端设备30,所述在线检测装置10、数据处理装置20和终端设备30之间通过有线网络连接;
在线检测装置10,用于在线获取水质中的重金属浓度数据并传送至数据处理装置20,其包括水样采集模块11、在线消解模块12、重金属富集模块13以及光谱检测模块14,水样采集模块11按一定时间间隔连续自动采集水样,所述在线消解模块12对水样进行消化处理,所述重金属富集模块13对消化处理后的水样进行富集和洗脱,所述光谱检测模块14对洗脱出来的水样进行重金属浓度检测,并将重金属浓度数据传送至数据处理装置20;
数据处理装置20,用于根据重金属浓度数据获得水质平均污染指数和偏离系数,并向终端设备发送预警信号,具体地,数据处理装置上设有水质评价模块21和预警模块22,水质评价模块21根据重金属浓度数据获得平均污染指数,预警模块22利用指数平滑法获得重金属浓度的预测值,并将实际重金属浓度与预测值相对比获得偏离系数,如果平均污染指数和偏离系数超出正常阈值,则向终端设备发送预警信号;
终端设备30,用于接收和显示来自数据处理装置的数据和信号。
在线检测装置的工作过程为:水样采集模块按一定时间间隔连续自动采集水样,经适当的沉降和过滤后,进入在线消解模块对水样进行消化处理,使其达到光谱检测的要求;将消化处理后的溶液泵入到填充有活性炭的重金属富集模块中进行富集,通过一定倍数的富集后,再用相应的洗脱液进行洗脱使得水中痕量的重金属达到仪器的检测限;利用光谱检测模块对洗脱出来的溶液同时进行重金属检测,并收集检测数据。
本实用新型通过自动采集水样的监测数据,根据预置的水质评价模块和预警模块判断水质状况,能减少人的工作量,节约时间,并通过平均指数模型可避免某个异常值达到阈值而导致的频发预警,在风险评价上更接近于真实水样;而且本实用新型能及时自动向用户发出预警信号,可实现无人值守监控。
在本实施例中,本在线预警系统还包括数据库服务器装置40,所述的数据库服务器装置40用于存取水质实时监测数据及历史数据,仪器运行日志,包括已确认的每日、周、月的水质检测数据,以及水质状况和污染预警事件信息。
在本实施例中,所述在线检测装置为基于原子发射光谱技术的重金属在线监测设备。基于原子发射光谱技术的重金属在线监测设备可以对水中多种重金属实现同时在线监测,而且该设备使用空气或者氮气作为工作气体,在安全性能提高的同时也降低了仪器的运行成本。
在本实施例中,所述重金属包括Cd、Cu、Cr、Co、Fe、Mn、Ni、Ti、V、Pb、Zn和Bi。也即是本实用新型可以同时对多种重金属(包括Cd、Cu、Cr、Co、Fe、Mn、Ni、Ti、V、Pb、Zn、Bi)进行检测,并收集检测数据。
请参考图2,一种水中重金属污染风险在线预警方法,其包括以下步骤;
通过在线检测获得时间t时水质中的重金属浓度数据;
通过下式计算单项污染指数以及平均污染指数,计算公式如下:
式中:K—平均污染指数;Ci—某金属浓度的实际测定值;C0—某金属浓度的水质标准限值;N—金属的种类数;
通过下式(1)计算出在时间t时某重金属的浓度平均值,
式中:Xt为在时间t时某重金属的浓度平均值,Xi为在时间t时第i个某重金属的浓度数据,n为某重金属浓度数据个数;
根据某重金属的浓度平均值,利用指数平滑法,计算时间t的一次指数平滑值,通过公式(2)计算:
St (1)=aXt+(1-a)St-1 (1) (2)
式中:St (1)为在时间t时某重金属浓度的一次指数平滑值;a为常数,取值范围为[0,1];
St-1 (1)为时间t-1时某重金属浓度的一次指数平滑值;
利用指数平滑法,计算时间t的二次指数平滑值,通过公式(3)计算:
St (2)=aSt (1)+(1-a)St-1 (2) (3)
St (2)为在时间t时某重金属浓度的二次指数平滑值,St-1 (2)为在时间t-1时某重金属浓度的二次指数平滑值;
利用指数平滑法,计算时间t的三次指数平滑值,通过公式(4)计算:
St (3)=aSt (2)+(1-a)St-1 (3) (4)
St (3)为在时间t时某重金属浓度的三次指数平滑值,St-1 (3)为在时间t-1时某重金属浓度的三次指数平滑值;
利用三次指数平滑法计算某重金属在时间t+1时的浓度预测值Pt+1,具体地:
Pt+1=at+bt+ct (5)
at=3St (1)-3St (2)+St (3) (6)
根据某重金属在时间t+1时的浓度预测值Pt+1,通过与在时间t+1时某重金属的实际值P相对比,计算出偏离系数b:
如果平均污染指数K和偏离系数均b超过设定阈值,则向用户发送预警信号;否则,不发送预警信号。
所述水质的扩散模型为一维水质扩散模型。一维水质扩散模型的假定条件为横向和垂直混合速度极快,被监测的断面污染物浓度均匀,污染物的浓度仅在河流的纵向上发生变化,因此适用于较为宽浅、重金属能够在短时间混合均匀的河流。
其中,所述预警信号设有5个档次,按偏离系数b的大小分级,分别为一级预警、二级预警、三级预警、四级预警和五级预警。
以下通过具体例子对利用上述在线预警系统进行风险预警的讲解:
目前某水厂已将本实用新型所述的设备投入到了实际使用中,用于监测水源地的重金属浓度,因此,本实用新型实施例选择其水源地作为对象,应用本实用新型的预警系统进行水质评价以及风险等级评价。某水厂水源地的取水口位于经度113°26’12”、纬度22°35’50”,水厂的服务人口为30万人,年取水量为8090万t,为典型的河流型水源地,污染物的平流输送远比扩散显著,适用于一维水质扩散模型。
利用上述在线预警系统对水源地进行水质评价,具体过程为:
首先根据水质标准确定待监测重金属的标准值限值C0。
标准值限值C0的会直接影响单项污染指数以及平均污染指数的大小,C0过高或过均会严重影响水质的分类以及预警等级的判断。目前我国国内规定了水中重金属浓度限值的标准主要有环保部办公厅在2012年印发的《全国集中式生活饮用水水源地水质监测实施方案》、《地表水环境质量标准》(GB3838-2002)以及《地下水质量标准》(GB/T14848-2007),本实用新型实施例主要是参考对比了以上几项标准,确定了各类重金属的标准值C0。
表1各类重金属监测浓度标准值
重金属 |
Zn |
Cd |
Pb |
Cu |
Cr |
Fe |
Ni |
Mn |
C0≤ |
1.0 |
0.005 |
0.05 |
1.0 |
0.05 |
0.3 |
0.02 |
0.1 |
第一步:通过在线检测获得时间t时水质中的重金属浓度数据。
下表为使用本实用新型所述的重金属在线监测设备得到的不同时刻的各类重金属浓度。
表2重金属监测数据(单位:mg/L)
第二步:通过下式计算单项污染指数以及平均污染指数,计算公式如下:
式中:K—平均污染指数;Ci—某金属浓度的实际测定值;C0—某金属浓度的标准值;N—金属的种类数。
第三步:根据计算所得平均污染指数通过参照污染指数级别分类的阈值从而判断水质状况,具体的水质分类阈值级别分类如下:
表3平均污染指数级别阈值分类
水质级别 |
Ⅰ |
Ⅱ |
Ⅲ |
Ⅳ |
Ⅴ |
K |
≤0.3 |
0.3~0.4 |
0.4~0.7 |
0.7~1.0 |
≥2.0 |
具体的计算结果如下表:
表4预警级别及污染程度评价结果
从上述结果可以看出,此河流支流的水质状况较为良好,在0:00~8:00的监测时刻内均属于Ⅰ、Ⅱ类水质,在12:00的时刻水质情况较差,属于Ⅲ类水质,其中Fe的Pi值达到了2.36,污染程度较为严重。
利用上述在线预警系统进行风险在线评价方法,具体过程为:
所述的在线预警系统的预警模块是根据重金属实际浓度与预测浓度的偏离大小评价风险程度的,这种方法具有计算简便、实用性较强的优点。本实施例以某河流支流在1~8月份的在8.00监测时刻的历史监测数据为例,按照下列步骤具体说明,步骤如下:
步骤S1:通过在线检测获得1~8月份重金属在8.00时刻时水质中的重金属浓度数据;
步骤S2:通过下式(1)计算出在时间t时某重金属的浓度平均值,
式中:Xt为在时间t时某重金属的浓度平均值,Xi为在时间t时第i个某重金属的浓度数据,n为某重金属浓度数据个数;
本实用新型实施例以实际的监测数据为例,计算得到的1~8月份重金属在8.00时刻的监测浓度平均值,具体数据如下:
表5某河流支流1~8月份在8:00时刻的监测浓度平均值单位(μg/L)
月份 |
Zn |
Cd |
Pb |
Cu |
Cr |
Fe |
Ni |
Mn |
1 |
10.6 |
2.06 |
4.6 |
9.45 |
2.01 |
82.1 |
9.1 |
56 |
2 |
8.05 |
1.05 |
4.3 |
9.2 |
3.45 |
84.4 |
9.5 |
47.2 |
3 |
14 |
3.15 |
2.35 |
5.1 |
4.08 |
170 |
7.04 |
12 |
4 |
15.5 |
4.22 |
3.5 |
4.99 |
2.5 |
260 |
8.03 |
44 |
5 |
10.6 |
3.12 |
5.4 |
8.21 |
3.21 |
78.5 |
9.1 |
56 |
6 |
12.5 |
4.22 |
3.5 |
5.02 |
2.5 |
150 |
8.46 |
44 |
7 |
14 |
3.12 |
2.35 |
4.78 |
3.78 |
170 |
8.04 |
12 |
8 |
5.05 |
1.99 |
2.13 |
8.65 |
3.45 |
84.4 |
9.5 |
47.2 |
步骤S3:根据某重金属的浓度平均值,利用指数平滑法,计算时间t的一次指数平滑值,通过公式(2)计算:
St (1)=aXt+(1-a)St-1 (1) (2)
式中:St (1)为在时间t时某重金属浓度的一次指数平滑值;a为常数,取值范围为[0,1];
St-1 (1)为时间t-1时某重金属浓度的一次指数平滑值;
具体地,初始值即St-1 (1)是根据前三个月的平均值作为起始浓度预测值;
表6某河流支流1~8月份在8:00时刻的监测浓度的一次指数平滑值
步骤S4:根据步骤3所得的一次指数平滑值,利用指数平滑法,计算时间t的二次指数平滑值,通过公式(3)计算:
St (2)=aSt (1)+(1-a)St-1 (2) (3)
St (2)为在时间t时某重金属浓度的二次指数平滑值,St-1 (2)为在时间t-1时某重金属浓度的二次指数平滑值;
具体地,所得的二次指数平滑值计算结果如下:
表7某河流支流1~8月份在8:00时刻的监测浓度的二次指数平滑值
步骤S5:根据步骤S4所得的二次指数平滑值,利用指数平滑法,计算时间t的三次指数平滑值,通过公式(4)计算:
St (3)=aSt (2)+(1-a)St-1 (3) (4)
St (3)为在时间t时某重金属浓度的三次指数平滑值,St-1 (3)为在时间t-1时某重金属浓度的三次指数平滑值;
表8某河流支流1~8月份在8:00时刻的监测浓度的三次指数平滑值
月份 |
Zn |
Cd |
Pb |
Cu |
Cr |
Fe |
Ni |
Mn |
起始值 |
10.88 |
2.09 |
3.75 |
7.92 |
3.18 |
112.17 |
8.55 |
38.40 |
1 |
10.85 |
2.08 |
3.86 |
8.11 |
3.03 |
108.41 |
8.62 |
40.60 |
2 |
10.48 |
1.95 |
3.98 |
8.36 |
2.99 |
103.06 |
8.77 |
42.80 |
3 |
10.69 |
2.02 |
3.84 |
8.09 |
3.12 |
108.55 |
8.64 |
40.05 |
4 |
11.47 |
2.35 |
3.69 |
7.50 |
3.13 |
131.59 |
8.46 |
38.55 |
5 |
11.82 |
2.65 |
3.83 |
7.26 |
3.13 |
138.66 |
8.45 |
40.14 |
6 |
12.03 |
2.99 |
3.90 |
6.90 |
3.05 |
141.62 |
8.46 |
41.80 |
7 |
12.36 |
3.18 |
3.73 |
6.44 |
3.09 |
146.13 |
8.42 |
38.92 |
8 |
11.63 |
3.11 |
3.41 |
6.47 |
3.17 |
140.87 |
8.53 |
37.94 |
步骤S6:利用三次指数平滑法计算某重金属在时间t+1时的浓度预测值Pt+1,具体地:
Pt+1=at+bt+ct (5)
at=3St(1)-3St(2)+St(3) (6)
表9某河流支流1~8月份在8:00时刻的监测浓度的浓度预测值
月份 |
Zn |
Cd |
Pb |
Cu |
Cr |
Fe |
Ni |
Mn |
1 |
10.46 |
2.05 |
5.03 |
10.22 |
1.43 |
67.07 |
9.38 |
64.80 |
2 |
6.63 |
0.53 |
4.58 |
9.84 |
3.59 |
70.52 |
9.98 |
51.60 |
3 |
15.63 |
3.69 |
1.44 |
3.31 |
4.82 |
206.43 |
6.15 |
0.00 |
4 |
18.57 |
5.55 |
3.08 |
2.92 |
2.31 |
346.50 |
7.43 |
41.30 |
5 |
10.24 |
3.57 |
6.39 |
8.68 |
3.06 |
55.23 |
9.52 |
68.10 |
6 |
11.84 |
4.65 |
3.62 |
4.66 |
2.17 |
124.66 |
8.72 |
48.98 |
7 |
14.56 |
2.91 |
1.40 |
3.91 |
4.11 |
175.06 |
7.86 |
0.00 |
8 |
1.25 |
0.98 |
1.15 |
10.07 |
3.77 |
51.35 |
10.00 |
47.41 |
步骤S7:根据某重金属在时间t+1时的浓度预测值Pt+1,通过与在时间t+1时某重金属的实际值P相对比,计算出偏离系数b:
综上所述,如果平均污染指数K和偏离系数均b超过设定阈值,则向用户发送预警信号;否则,不发送预警信号。
具体地,本实用新型实施例根据9月份某日在8.00时间监测时刻的数据进行计算,计算结果如下表所示。
表10某河流支流9月份在8:00时刻的某次实际监测浓度与偏离系数
根据所得的平均污染指数K和偏离系数b,判断是否向用户发出预警。
具体地,由于当日8:00时刻的重金属浓度检测数据的单项污染指数属于正常范围,因此按照预警判断树的判断结果,当日8:00时刻的水质为正常水质,不进行预警。
上述实施例仅用以说明本专利而并非限制本专利所描述的技术方案;因此,尽管本说明书参照上述的各个实施例对本专利已进行了详细的说明,但是,本领域的普通技术人员应当理解,仍然可以对本专利进行修改或者等同替换;而一切不脱离本专利的精神和范围的技术方案及其改进,其均应涵盖在本专利的权利要求范围当中。