本发明涉及水体中新兴污染物的，涉及一种基于蒙特卡罗模拟的含抗生素污水优化处理方法。
背景技术：
随着经济的快速发展，水体中新型污染物的得到了大量的使用。新兴污染物主要包括抗生素、消炎药、中枢兴奋药等常见药用化合物以及护肤品、化妆品、洗涤剂等日常护理用品。通常，新型污染物会随着污水一起流入城市生活污水处理厂。抗生素作为一种报道次数多、检出率高的一种新兴污染物，常常难以被常规污水处理工艺(如活性污泥法)去除。研究表明，深度处理工艺可以比较有效地实现对污染物的去除，而实际情况中城市污水处理厂常采用常规工艺与深度处理工艺的组合。因此，如何为抗生素选择最优的工艺组合来实现最大的去除率成为了关键问题。
由于通常情况下污水中抗生素的浓度极低，常使用高效液相色谱-质谱联用技术测定污水中的抗生素浓度。这种方法存在操作繁琐、费用较高等问题。近年来，一种应用蒙特卡罗模拟来研究环境中新兴污染物归趋的数学方法得到了广泛的应用。蒙特卡罗模拟是一种利用随机数进行数值模拟的方法，主要依据概率统计理论中的中心极限定理和大数定律，并在此基础上通过随机抽样的方式反映所研究问题发生的某种规律。
技术实现要素：
本发明的目的在于克服传统方法中操作繁琐、费用较高的问题，提供一种基于蒙特卡罗模拟的新兴污染物污水优化处理方法，该方法提高了污水处理工艺组合对抗生素的处理率，同时降低了城市污水处理厂最终出水中抗生素的环境风险。
本发明的技术方案:基于蒙特卡罗模拟的新兴污染物污水优化处理方法,包括以下步骤：
(1)确定自变量，包括以下内容：
a)、城市污水处理厂进水中抗生素的浓度；
b)、二级处理工艺对抗生素的去除率范围；
c)、三级处理工艺对抗生素的去除率范围；
(2)描述自变量的概率分布，采用matlab软件，使用混合高斯模型对自变量分布进行拟合；
(3)在拟合得到的抗生素进水浓度分布函数中，进行随机抽样；
(4)选择二级和三级污水处理工艺组合；
(5)根据工艺组合，在去除率分布函数中分别进行随机抽样；
(6)计算总去除率及出水中抗生素的浓度；
(7)随机抽样得出抗生素在任何一种工艺组合下的去除率范围及均值；
(8)根据抗生素的类型选择确定最优的工艺组合。
与现有技术相比，本发明具有的优点：
(1)本发明得到的去除率是在查阅大量文献调研的基础上，并结合污水处理厂实际的去除率，进行十万次蒙特卡罗随机抽样过程的结果，具有一定的数学意义和物理意义，是可信的；
(2)采用蒙特卡罗模拟法，避免了传统方法中操作繁琐、费用较高的问题，可以根据抗生素自身的物理化学性质选择最优的污水处理工艺组合，从而提高污水处理厂的处理效率。
附图说明
图1是基于蒙特卡罗模拟的新兴污染物污水优化处理方法的流程示意图。
具体实施方式
下面通过具体实施例和附图对本发明作进一步的说明。本发明的实施例是为了更好地使本领域的技术人员更好地理解本发明，并不对本发明作任何的限制。
本发明选取天津市某污水处理厂作为研究对象，该厂进水中典型抗生素包括磺胺甲恶唑(smx)、甲氧苄啶(tmp)、红霉素(rox)、罗红霉素(rox)、诺氟沙星(nor)、环丙沙星(cip)。
自变量抗生素进水浓度c进水根据实际情况确定，二级去除率r1、三级去除率r2，通过调研大量文献确定，如表1所示：
表1抗生素的二级去除率r1、三级去除率r2范围
表中：cas-传统活性污泥法；a/o-厌氧-好氧法；a/a/o-厌氧-缺氧-好氧法；od-氧化沟；ro-反渗透；pac-活性炭吸附；o3-臭氧氧化。
根据已知的抗生素进水浓度c进水，依托matlab软件，选择混合高斯模型对二级去除率r1、三级去除率r2的分布进行拟合。
混合高斯模型的概率分布密度函数可以通过m个单高斯模型的加权函数表示：
其中，
表示第j个单高斯模型。
公式中，xj表示数据值，在本发明中分别为抗生素进水浓度、二级工艺去除率和三级工艺去除率；αj为不同单高斯模型在混合高斯模型中的比例；c为协方差矩阵，μ为期望值。
在拟合得到的抗生素进水浓度分布函数f(c进水)中，进行10万次随机抽取一个浓度值。
选择cas、a/o、a/a/o、od为四种二级工艺，ro、o3、pac为三种三级工艺进行工艺组合，每种抗生素在二级和三级工艺组合下模拟10万次，每一次随机过程得到一个最终的去除率；
按照表2中的二级、三级工艺组合方式，选择其中一种待研究的工艺组合。
根据上述选定的工艺组合，在其去除率分布函数f(r1)、f(r2)中分别进行10万次随机抽样。
根据随机抽取的平均抗生素进水浓度c进水、二级去除率r1、三级去除率r2，得出总去除率r总，如表2所示，根据相应进水浓度可以计算得出抗生素出水浓度c出水。
表2蒙特卡罗模拟后抗生素总去除率(％)
随机抽样过程重复10万次，可得出抗生素在某一种工艺组合下的去除率范围及均值。
最终可以根据抗生素的类型选择最优的工艺组合，实现最佳的去除效果。
应当理解的是，这里所讨论的实施方案及实例只是为了说明，对本领域技术人员来说，可以加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。