CN107585890A - 基于蒙特卡罗模拟的新兴污染物污水优化处理方法 - Google Patents
基于蒙特卡罗模拟的新兴污染物污水优化处理方法 Download PDFInfo
- Publication number
- CN107585890A CN107585890A CN201710864299.9A CN201710864299A CN107585890A CN 107585890 A CN107585890 A CN 107585890A CN 201710864299 A CN201710864299 A CN 201710864299A CN 107585890 A CN107585890 A CN 107585890A
- Authority
- CN
- China
- Prior art keywords
- antibiotic
- clearance
- random sampling
- combination
- monte carlo
- 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
Landscapes
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
Abstract
本发明涉及水体中新兴污染物的,涉及一种基于蒙特卡罗模拟的含抗生素污水优化处理方法,主要包括:确定自变量;采用Matlab软件,使用混合高斯模型对自变量分布进行拟合;在拟合得到的抗生素进水浓度分布函数中,进行随机抽样;选择二级三级污水处理工艺组合;根据工艺组合,在去除率分布函数中分别进行随机抽样;计算总去除率及出水中抗生素的浓度;随机抽样得出抗生素在某一种工艺组合下的去除率范围及均值;根据抗生素的类型选择确定最优的工艺组合。采用蒙特卡罗模拟法,避免了传统方法中操作繁琐、费用较高的问题。
Description
技术领域
本发明涉及水体中新兴污染物的,涉及一种基于蒙特卡罗模拟的含抗生素污水优化处理方法。
背景技术
随着经济的快速发展,水体中新型污染物的得到了大量的使用。新兴污染物主要包括抗生素、消炎药、中枢兴奋药等常见药用化合物以及护肤品、化妆品、洗涤剂等日常护理用品。通常,新型污染物会随着污水一起流入城市生活污水处理厂。抗生素作为一种报道次数多、检出率高的一种新兴污染物,常常难以被常规污水处理工艺(如活性污泥法)去除。研究表明,深度处理工艺可以比较有效地实现对污染物的去除,而实际情况中城市污水处理厂常采用常规工艺与深度处理工艺的组合。因此,如何为抗生素选择最优的工艺组合来实现最大的去除率成为了关键问题。
由于通常情况下污水中抗生素的浓度极低,常使用高效液相色谱-质谱联用技术测定污水中的抗生素浓度。这种方法存在操作繁琐、费用较高等问题。近年来,一种应用蒙特卡罗模拟来研究环境中新兴污染物归趋的数学方法得到了广泛的应用。蒙特卡罗模拟是一种利用随机数进行数值模拟的方法,主要依据概率统计理论中的中心极限定理和大数定律,并在此基础上通过随机抽样的方式反映所研究问题发生的某种规律。
发明内容
本发明的目的在于克服传统方法中操作繁琐、费用较高的问题,提供一种基于蒙特卡罗模拟的新兴污染物污水优化处理方法,该方法提高了污水处理工艺组合对抗生素的处理率,同时降低了城市污水处理厂最终出水中抗生素的环境风险。
本发明的技术方案:基于蒙特卡罗模拟的新兴污染物污水优化处理方法,包括以下步骤:
(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万次,可得出抗生素在某一种工艺组合下的去除率范围及均值。
最终可以根据抗生素的类型选择最优的工艺组合,实现最佳的去除效果。
应当理解的是,这里所讨论的实施方案及实例只是为了说明,对本领域技术人员来说,可以加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (1)
1.基于蒙特卡罗模拟的新兴污染物污水优化处理方法,包括以下步骤:
(1)确定自变量,包括以下内容:
a)、城市污水处理厂进水中抗生素的浓度;
b)、二级处理工艺对抗生素的去除率范围;
c)、三级处理工艺对抗生素的去除率范围;
(2)描述自变量的概率分布,采用Matlab软件,使用混合高斯模型对自变量分布进行拟合;
(3)在拟合得到的抗生素进水浓度分布函数中,进行随机抽样;
(4)选择二级和三级污水处理工艺组合;
(5)根据工艺组合,在去除率分布函数中分别进行随机抽样;
(6)计算总去除率及出水中抗生素的浓度;
(7)随机抽样得出抗生素在任何一种工艺组合下的去除率范围及均值;
(8)根据抗生素的类型选择确定最优的工艺组合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710864299.9A CN107585890B (zh) | 2017-09-22 | 2017-09-22 | 基于蒙特卡罗模拟的新兴污染物污水优化处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710864299.9A CN107585890B (zh) | 2017-09-22 | 2017-09-22 | 基于蒙特卡罗模拟的新兴污染物污水优化处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107585890A true CN107585890A (zh) | 2018-01-16 |
CN107585890B CN107585890B (zh) | 2020-07-24 |
Family
ID=61048312
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710864299.9A Expired - Fee Related CN107585890B (zh) | 2017-09-22 | 2017-09-22 | 基于蒙特卡罗模拟的新兴污染物污水优化处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107585890B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111652445A (zh) * | 2020-06-11 | 2020-09-11 | 广东科创工程技术有限公司 | 基于高斯分布的污水设备优化运行控制方法 |
CN113673117A (zh) * | 2021-09-03 | 2021-11-19 | 天津大学 | 基于蒙特卡罗法模拟构建海油污染物减排运行成本的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100094202A1 (en) * | 2005-06-17 | 2010-04-15 | Bayer Technology Services Gmbh | Device for the Time-controlled Intravenous Administering of the Anesthetic Propofol |
CN101759236A (zh) * | 2009-12-31 | 2010-06-30 | 南京大学 | 一种工业园区点源污染物的调控方法 |
CN103810537A (zh) * | 2014-02-12 | 2014-05-21 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
CN105095997A (zh) * | 2015-07-30 | 2015-11-25 | 浙江大学 | 基于蒙特卡洛和层次分析法的突发性水污染事故预警方法 |
-
2017
- 2017-09-22 CN CN201710864299.9A patent/CN107585890B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100094202A1 (en) * | 2005-06-17 | 2010-04-15 | Bayer Technology Services Gmbh | Device for the Time-controlled Intravenous Administering of the Anesthetic Propofol |
CN101759236A (zh) * | 2009-12-31 | 2010-06-30 | 南京大学 | 一种工业园区点源污染物的调控方法 |
CN103810537A (zh) * | 2014-02-12 | 2014-05-21 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
CN105095997A (zh) * | 2015-07-30 | 2015-11-25 | 浙江大学 | 基于蒙特卡洛和层次分析法的突发性水污染事故预警方法 |
Non-Patent Citations (1)
Title |
---|
王林媛: "工业园区水污染防治技术路线的优化方法研究", 《万方硕士学位论文》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111652445A (zh) * | 2020-06-11 | 2020-09-11 | 广东科创工程技术有限公司 | 基于高斯分布的污水设备优化运行控制方法 |
CN111652445B (zh) * | 2020-06-11 | 2024-03-22 | 广东科创智水科技有限公司 | 基于高斯分布的污水设备优化运行控制方法 |
CN113673117A (zh) * | 2021-09-03 | 2021-11-19 | 天津大学 | 基于蒙特卡罗法模拟构建海油污染物减排运行成本的方法 |
CN113673117B (zh) * | 2021-09-03 | 2024-01-26 | 天津大学 | 基于蒙特卡罗法模拟构建海油污染物减排运行成本的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107585890B (zh) | 2020-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Haynes et al. | Diatom-derived carbohydrates as factors affecting bacterial community composition in estuarine sediments | |
Mohd Hanafiah et al. | Performance of wild-Serbian Ganoderma lucidum mycelium in treating synthetic sewage loading using batch bioreactor | |
Liu et al. | Shifting nutrient-mediated interactions between algae and bacteria in a microcosm: evidence from alkaline phosphatase assay | |
Ortiz-Vera et al. | Influence of water quality on diversity and composition of fungal communities in a tropical river | |
CN107585890A (zh) | 基于蒙特卡罗模拟的新兴污染物污水优化处理方法 | |
Mu et al. | Impacts of water flow on epiphytic microbes and nutrients removal in constructed wetlands dominated by Vallisneria natans with decreasing temperature | |
Chatterjee et al. | Mycoremediation of 17 β‐estradiol using trichoderma citrinoviride strain AJAC3 along with enzyme studies | |
Zong et al. | Taxonomy and phylogeny of four new species in Absidia (Cunninghamellaceae, Mucorales) from China | |
Del Rio et al. | First record of an Icacinaceae Miers fossil flower from Le Quesnoy (Ypresian, France) amber | |
Addo-Bankas et al. | Challenges of urban artificial landscape water bodies: Treatment techniques and restoration strategies towards ecosystem services enhancement | |
Khemmuk et al. | Wongia gen. nov.(Papulosaceae, Sordariomycetes), a new generic name for two root-infecting fungi from Australia | |
Mehrabi et al. | Description of Allocanariomyces and Parachaetomium, two new genera, and Achaetomium aegilopis sp. nov. in the Chaetomiaceae | |
Tokuda et al. | First real-time observation of transverse division in azooxanthellate scleractinian corals | |
Froufe et al. | The evolutionary history of sharp-and blunt-snouted lenok (Brachymystax lenok (Pallas, 1773)) and its implications for the paleo-hydrological history of Siberia | |
Tseplik et al. | Achnanthidiumgladius sp. nov.(Bacillariophyceae)–a new monoraphid diatom species from Indonesia | |
CN102329046B (zh) | 一种光催化还原和沉水植物联合处理沉积物磷的方法 | |
Xiong et al. | Comparative chloroplast genomics reveals the phylogeny and the adaptive evolution of Begonia in China | |
Ma et al. | Impacts of River Damming on Sediment Methylmercury Dynamics | |
CN205473215U (zh) | 生活污水处理装置 | |
Zhang et al. | Removal of Nitrate Nitrogen from Municipal Wastewater Using Autotrophic Denitrification Based on Magnetic Pyrite | |
Hadei et al. | Determination of biokinetic coefficients for an adsorption/bio-oxidation process on municipal wastewater in pilot-scale | |
Schmidt et al. | Groundwater ecosystems in changing times | |
Jamshidi et al. | Performance evaluation and comparison of facultative ponds in series and parallel for wastewater treatment (Case study: Delijan WWTP) | |
CN108946998A (zh) | 一种污水的深度处理工艺 | |
Vairavel et al. | Residence time distribution studies and modeling of rotating biological contactor reactor for decolorization of Congo red from synthetic dye wastewater |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200724 Termination date: 20200922 |