CN115759841B - 一种基于统计年鉴面板数据的区域下水道甲烷计算方法 - Google Patents
一种基于统计年鉴面板数据的区域下水道甲烷计算方法 Download PDFInfo
- Publication number
- CN115759841B CN115759841B CN202211458691.0A CN202211458691A CN115759841B CN 115759841 B CN115759841 B CN 115759841B CN 202211458691 A CN202211458691 A CN 202211458691A CN 115759841 B CN115759841 B CN 115759841B
- Authority
- CN
- China
- Prior art keywords
- sewer
- methane
- bod
- total
- sewage treatment
- 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.)
- Active
Links
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种基于统计年鉴面板数据的区域下水道甲烷计算方法,包括以下步骤:获取待测区域的相关信息,所述相关信息包括城市人口信息、城市污水处理和污水排放信息;获取待测区域的产污系数,基于所述相关信息和所述产污系数获取污染物产生量;获取待测区域的下水道甲烷排放因子,基于所述污染物产生量和所述下水道甲烷排放因子,采用排放因子法获取待测区域的下水道甲烷排放总量。本发明提供了一种基于统计年鉴面板数据的区域下水道甲烷计算方法,采用排放因子法为基础从宏观角度计算下水道甲烷排放量,方法简单,可操作性强。
Description
技术领域
本发明属于市政环境管理技术领域,尤其涉及一种基于统计年鉴面板数据的区域下水道甲烷计算方法。
背景技术
IPCC的第六份评估报告强调了甲烷带来的气候挑战,甲烷占目前人为温室气体驱动的变暖的0.5℃,约三分之一,仅次于二氧化碳。在各种甲烷排放源中,微生物源估计会释放300-400百万吨每年的甲烷排放量,或超过甲烷总排放量的50%。
废水含有高有机含量,使其成为甲烷的重要微生物来源,2015年污水处理厂排放了全球非二氧化碳温室气体排放的9.6%。除了污水处理厂外,向污水处理厂输送污水的下水道是另一个关键的甲烷来源,但量化污水收集***的甲烷排放的研究有限。Yuan等人对澳大利亚两种不同的上升主下水道进行了监测,证明下水道***产生了大量的甲烷。一项对城市街道空气的分析也证实了下水道管道释放的生物甲烷。Shah等人试图通过直接测量来自排放人孔的甲烷气体通量来量化下水道***的甲烷排放。目前大部分专利申请以下水道甲烷监测、收集设备为主,如专利公开号CN111912935A发明了一种可监测下水道沼气含量的智能井盖,结构设计简单、可在线对下水道甲烷浓度进行实时监测,一旦发生甲烷浓度超过预设值,将实时向服务器上传报警信号,专利CN209778852U则公开了一种城市道路下水道甲烷分离收集装置,其采用向下排空气法将城市道路下水道中空气中的甲烷收集出来,主要通过上过滤体和下过滤体过滤空气杂质和水分,处理后的水通过出水阀排出,实现了甲烷资源的回收利用。
然而,下水道中产生的甲烷将被释放到重力污水管道、泵站或污水处理厂入口的大气中,因为在下水道中没有发现显著的甲烷汇或甲烷营养活动。因此,虽然在某一点位测量甲烷浓度具有较高的精度,但很难将甲烷测量结果归因于究竟是废水的生物生成、破损下水道的泄漏抑或是重力管道中的气体扩散,这对于整个城市下水道***中甲烷排放量的宏观估算并不适用。
发明内容
本发明的目的在于提出一种基于统计年鉴面板数据的区域下水道甲烷计算方法,采用排放因子法为基础从宏观角度计算下水道甲烷排放量,方法简单,可操作性强。
为实现上述目的,本发明提供了一种基于统计年鉴面板数据的区域下水道甲烷计算方法,包括以下步骤:
获取待测区域的相关信息,所述相关信息包括城市人口信息、城市污水处理和污水排放信息;
获取待测区域的产污系数,基于所述相关信息和所述产污系数获取污染物产生量;
获取待测区域的下水道甲烷排放因子,基于所述污染物产生量和所述下水道甲烷排放因子,采用排放因子法获取待测区域的下水道甲烷排放总量。
可选的,所述相关信息的来源包括:
所述城市人口信息从《城市建设统计年鉴》获取,所述城市人口信息包括城镇常住人口和暂住人口,镇常住人口和暂住人口相加为计算参数所需城镇总人口数;
所述城市污水处理和污水排放信息包括污水排放量、污水处理量、污水处理厂进水口BOD年均浓度;所述污水排放量和所述污水处理量均从《城市建设统计年鉴》获取,所述污水处理厂进水口BOD年均浓度由目标城市当地所有污水处理厂的一年连续监测数据经过加权平均获得。
可选的,获取待测区域的产污系数的方法为:
α=k·β(1)
其中,α为产污系数,克BOD/天;k为BOD与化学需氧量(COD)的比值,k取值参考《省级温室气体清单编制指南》;β为排放生活废水中COD的日人均负荷,克COD/天。
可选的,基于所述相关信息和所述产污系数获取污染物产生量,所述污染物产生量包括:待测区域产生的生活污水中的BOD总量、污水处理的BOD总量、直接排放至地表水的BOD总量以及被下水道降解的BOD总量。
可选的,获取污染物产生量的方法为:
M0=365·α·P/100 (2)
M1=0.01·Q1·c (3)
M2=0.01· (Q0-Q1)·c (4)
M3=M0-M1-M2 (5)
其中,M0为生活污水中BOD总量,单位为吨;M1为污水处理的BOD总量,单位为吨;M2为直接排放至地表水的BOD总量,单位为吨;M3为被下水道降解的BOD总量,单位为吨;α为产污系数,P为参数所需城镇总人口数,Q0为污水排放量,单位为万立方米;Q1为污水处理量,单位为万立方米;c为污水处理厂进水口BOD年均浓度。
可选的,获取待测区域的下水道甲烷排放因子具体包括:收集所述待测区域的年均温,根据所述待测区域的年均温获取所述待测区域的温度范围,根据所述温度范围获取对应的下水道甲烷排放因子。
可选的,获取待测区域的下水道甲烷排放因子的方法为:
EFj=B0·MCFj (6)
MCFj=-0.0009t2+0.0422t-0.0019 (7)
其中,EFj为下水道甲烷排放因子,B0为甲烷最大生产能力,单位为吨甲烷/吨BOD;MCFj为修正因子;t为温度,单位为摄氏度。
可选的,采用排放因子法获取待测区域的下水道甲烷排放总量的方法为:
G2=M3·EFj=M3·B0·MCFj (8)
其中,G为下水道甲烷排放总量,单位为吨;M3为被下水道降解的BOD总量,单位为吨;EFj为下水道甲烷排放因子,B0为甲烷最大生产能力,MCFj为修正因子。
本发明技术效果:本发明公开了一种基于统计年鉴面板数据的区域下水道甲烷计算方法,采用排放因子法为基础从宏观角度计算下水道甲烷排放量,方法简单,可操作性强,且从下水道甲烷的产生机理计算,避免了其他专利和研究只能进行小范围测量的局限性和由测量带来的误差,进一步推动了城市下水道甲烷的量化评估工作。
附图说明
构成本申请的一部分的附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为本发明实施例基于统计年鉴面板数据的区域下水道甲烷计算方法的流程示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机***中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
如图1所示,本实施例中提供一种基于统计年鉴面板数据的区域下水道甲烷计算方法,包括以下步骤:
获取待测区域的相关信息,所述相关信息包括城市人口信息、城市污水处理和污水排放信息;
获取待测区域的产污系数,基于所述相关信息和所述产污系数获取污染物产生量;
获取待测区域的下水道甲烷排放因子,基于所述污染物产生量和所述下水道甲烷排放因子,采用排放因子法获取待测区域的下水道甲烷排放总量。
进一步优化方案,所述相关信息的来源包括:
所述城市人口信息从《城市建设统计年鉴》获取,所述城市人口信息包括城镇常住人口和暂住人口,镇常住人口和暂住人口相加为计算参数所需城镇总人口数;
所述城市污水处理和污水排放信息包括污水排放量、污水处理量、污水处理厂进水口BOD年均浓度;所述污水排放量和所述污水处理量均从《城市建设统计年鉴》获取,所述污水处理厂进水口BOD年均浓度由目标城市当地所有污水处理厂的一年连续监测数据经过加权平均获得。
进一步优化方案,获取待测区域的产污系数的方法为:
α=k·β(1)
其中,α为产污系数,克BOD/天;k为BOD与化学需氧量(COD)的比值,k取值参考《省级温室气体清单编制指南》;β为排放生活废水中COD的日人均负荷,克COD/天。
进一步优化方案,基于所述相关信息和所述产污系数获取污染物产生量,所述污染物产生量包括:待测区域产生的生活污水中的BOD总量、污水处理的BOD总量、直接排放至地表水的BOD总量以及被下水道降解的BOD总量。
进一步优化方案,获取污染物产生量的方法为:
M0=365·α·P/100 (2)
M1=0.01·Q1·c (3)
M2=0.01· (Q0-Q1)·c (4)
M3=M0-M1-M2 (5)
其中,M0为生活污水中BOD总量,单位为吨;M1为污水处理的BOD总量,单位为吨;M2为直接排放至地表水的BOD总量,单位为吨;M3为被下水道降解的BOD总量,单位为吨;α为产污系数,P为参数所需城镇总人口数,Q0为污水排放量,单位为万立方米;Q1为污水处理量,单位为万立方米;c为污水处理厂进水口BOD年均浓度。
进一步优化方案,获取待测区域的下水道甲烷排放因子具体包括:收集所述待测区域的年均温,根据所述待测区域的年均温获取所述待测区域的温度范围,根据所述温度范围获取对应的下水道甲烷排放因子。
进一步优化方案,获取待测区域的下水道甲烷排放因子的方法为:
EFj=B0·MCFj (6)
MCFj=-0.0009t2+0.0422t-0.0019 (7)
其中,EFj为下水道甲烷排放因子,B0为甲烷最大生产能力,单位为吨甲烷/吨BOD;MCFj为修正因子;t为温度,单位为摄氏度。
进一步优化方案,采用排放因子法获取待测区域的下水道甲烷排放总量的方法为:
G2=M3·EFj=M3·B0·MCFj (8)
其中,G为下水道甲烷排放总量,单位为吨;M3为被下水道降解的BOD总量,单位为吨;EFj为下水道甲烷排放因子,B0为甲烷最大生产能力,MCFj为修正因子。
以北京市2019年数据为例。
Step1.1:北京市2019年总计城市人口1865万人,污水排放量199335.69万立方米,污水处理量197953.87万立方米;
Step1.2:北京市污水处理厂进水口BOD年均浓度c为97.9毫克/升;
Step2:北京市属于华北地区,k取为0.45,β查表为88克COD/天,则北京市:
α=k·β=0.45×88=39.6(克BOD/天)
Step3.1:北京市产生的生活污水中BOD总量M0(单位:吨):
M0=365·α·P/100=365×39.6×1865÷100=269567.1
Step3.2:
污水处理的BOD总量为M1(单位:吨):
M1=0.01·Q1·c=0.01×197953.87×97.9=193796.8
Step3.3:
直接排放至地表水的BOD总量为M2(单位:吨):
M2=0.01·(Q0-Q1)·c=0.01×(199335.69-197953.87)×97.9=1352.8
Step3.3:
被下水道降解的BOD总量为M3(单位:吨):
M3=M0-M1-M2=269567.1-193796.8-1352.8=74417.5
Step4.1:北京市下水道甲烷排放因子EFj(单位:吨甲烷/吨BOD)
EFj=B0·MCFj
其中,B0为0.6吨甲烷/吨BOD;北京市年均气温为13.8摄氏度
MCFj=-0.0009×(13.8)2+0.0422×(13.8)-0.0019=0.409
具体地,Step5的具体实施步骤如下:
下水道甲烷产量为G(单位:吨):
G2=M3·EFj=M3·B0·MCFj=74417.5×0.6×0.409=18262。
本发明公开了一种基于统计年鉴面板数据的区域下水道甲烷计算方法,采用排放因子法为基础从宏观角度计算下水道甲烷排放量,方法简单,可操作性强,且从下水道甲烷的产生机理计算,避免了其他专利和研究只能进行小范围测量的局限性和由测量带来的误差,进一步推动了城市下水道甲烷的量化评估工作。
以上所述,仅为本申请较佳的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应该以权利要求的保护范围为准。
Claims (4)
1.一种基于统计年鉴面板数据的区域下水道甲烷计算方法,其特征在于,包括以下步骤:
获取待测区域的相关信息,所述相关信息包括城市人口信息、城市污水处理和污水排放信息;
获取待测区域的产污系数,α=k·β (1)
其中,α为产污系数,克BOD/天;k为BOD与化学需氧量的比值,k取值参考《省级温室气体清单编制指南》;β为排放生活废水中COD的日人均负荷,克COD/天;
基于所述相关信息和所述产污系数获取污染物产生量,
M0=365·α·P/100 (2)
M1=0.01·Q1·c (3)
M2=0.01·(Q0-Q1)·c (4)
M3=M0-M1-M2 (5)
其中,M0为生活污水中BOD总量,单位为吨;M1为污水处理的BOD总量,单位为吨;M2为直接排放至地表水的BOD总量,单位为吨;M3为被下水道降解的BOD总量,单位为吨;α为产污系数,P为参数所需城镇总人口数,Q0为污水排放量,单位为万立方米;Q1为污水处理量,单位为万立方米;c为污水处理厂进水口BOD年均浓度;获取待测区域的下水道甲烷排放因子,
EFj=B0·MCFj (6)
MCFj=-0.0009t2+0.0422t-0.0019 (7)
其中,EFj为下水道甲烷排放因子,B0为甲烷最大生产能力,单位为吨甲烷/吨BOD;MCFj为修正因子;t为温度,单位为摄氏度;
基于所述污染物产生量和所述下水道甲烷排放因子,采用排放因子法获取待测区域的下水道甲烷排放总量,G2=M3·EFj=M3·B0·MCFj(8),其中,G2为下水道甲烷排放总量,单位为吨;M3为被下水道降解的BOD总量,单位为吨;EFj为下水道甲烷排放因子,B0为甲烷最大生产能力,MCFj为修正因子。
2.如权利要求1所述的基于统计年鉴面板数据的区域下水道甲烷计算方法,其特征在于,所述相关信息的来源包括:
所述城市人口信息从《城市建设统计年鉴》获取,所述城市人口信息包括城镇常住人口和暂住人口,城镇常住人口和暂住人口相加为计算参数所需城镇总人口数;
所述城市污水处理和污水排放信息包括污水排放量、污水处理量、污水处理厂进水口BOD年均浓度;所述污水排放量和所述污水处理量均从《城市建设统计年鉴》获取,所述污水处理厂进水口BOD年均浓度由目标城市当地所有污水处理厂的一年连续监测数据经过加权平均获得。
3.如权利要求1所述的基于统计年鉴面板数据的区域下水道甲烷计算方法,其特征在于,基于所述相关信息和所述产污系数获取污染物产生量,所述污染物产生量包括:待测区域产生的生活污水中的BOD总量、污水处理的BOD总量、直接排放至地表水的BOD总量以及被下水道降解的BOD总量。
4.如权利要求1所述的基于统计年鉴面板数据的区域下水道甲烷计算方法,其特征在于,获取待测区域的下水道甲烷排放因子具体包括:收集所述待测区域的年均温,根据所述待测区域的年均温获取所述待测区域的温度范围,根据所述温度范围获取对应的下水道甲烷排放因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211458691.0A CN115759841B (zh) | 2022-11-17 | 2022-11-17 | 一种基于统计年鉴面板数据的区域下水道甲烷计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211458691.0A CN115759841B (zh) | 2022-11-17 | 2022-11-17 | 一种基于统计年鉴面板数据的区域下水道甲烷计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115759841A CN115759841A (zh) | 2023-03-07 |
CN115759841B true CN115759841B (zh) | 2023-08-18 |
Family
ID=85333947
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211458691.0A Active CN115759841B (zh) | 2022-11-17 | 2022-11-17 | 一种基于统计年鉴面板数据的区域下水道甲烷计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115759841B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005081182A (ja) * | 2003-09-04 | 2005-03-31 | Univ Of Tsukuba | メタン発生方法及びそれに用いる二相式メタン発生装置 |
CN102161553A (zh) * | 2011-03-04 | 2011-08-24 | 广东工业大学 | 一种餐厨废物制沼气废水处理方法 |
WO2018134560A1 (en) * | 2017-01-23 | 2018-07-26 | Bactest Limited | Wastewater treatment plant control systems |
CN110555545A (zh) * | 2019-07-25 | 2019-12-10 | 南京大学 | 一种基于时空特征的生活污水产污负荷量测算方法 |
JP2020003274A (ja) * | 2018-06-26 | 2020-01-09 | 株式会社アイカムス・ラボ | 吐出量測定システム及び吐出量測定方法 |
CN112084649A (zh) * | 2020-09-04 | 2020-12-15 | 中国市政工程华北设计研究总院有限公司 | 一种基于城镇污水处理全流程环节的碳排放模型计算方法 |
CN113391045A (zh) * | 2021-03-31 | 2021-09-14 | 南开大学 | 一种测试厌氧污泥活性的方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2982953B1 (fr) * | 2011-11-23 | 2019-05-24 | Ams France | Procede, kit et systeme de mesure directe de multiples biodegradabilites. |
US11378563B2 (en) * | 2020-05-30 | 2022-07-05 | Ali Reza Abedini | Method for quantifying fugitive methane emissions rate using surface methane concentration |
US20220324737A1 (en) * | 2021-04-03 | 2022-10-13 | Power and Water Utility Company for Jubail and Yanbu (MARAFIQ) | Methods for monitoring marine waterbody sustainability at discharge sites |
-
2022
- 2022-11-17 CN CN202211458691.0A patent/CN115759841B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005081182A (ja) * | 2003-09-04 | 2005-03-31 | Univ Of Tsukuba | メタン発生方法及びそれに用いる二相式メタン発生装置 |
CN102161553A (zh) * | 2011-03-04 | 2011-08-24 | 广东工业大学 | 一种餐厨废物制沼气废水处理方法 |
WO2018134560A1 (en) * | 2017-01-23 | 2018-07-26 | Bactest Limited | Wastewater treatment plant control systems |
JP2020003274A (ja) * | 2018-06-26 | 2020-01-09 | 株式会社アイカムス・ラボ | 吐出量測定システム及び吐出量測定方法 |
CN110555545A (zh) * | 2019-07-25 | 2019-12-10 | 南京大学 | 一种基于时空特征的生活污水产污负荷量测算方法 |
CN112084649A (zh) * | 2020-09-04 | 2020-12-15 | 中国市政工程华北设计研究总院有限公司 | 一种基于城镇污水处理全流程环节的碳排放模型计算方法 |
CN113391045A (zh) * | 2021-03-31 | 2021-09-14 | 南开大学 | 一种测试厌氧污泥活性的方法 |
Non-Patent Citations (1)
Title |
---|
天津市生活污水甲烷排放量估算研究;孙祥超;陈伟;贺姝峒;张宁;;资源节约与环保(04);133-134,137 * |
Also Published As
Publication number | Publication date |
---|---|
CN115759841A (zh) | 2023-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Removal of heavy metals from urban stormwater runoff using bioretention media mix | |
Yu et al. | Water availability for shale gas development in Sichuan Basin, China | |
Crawford et al. | The significance of embodied energy in certified passive houses | |
Antunes et al. | Potential for potable water savings in buildings by using stormwater harvested from porous pavements | |
US10807037B2 (en) | Process for controlling air purification systems | |
Ma et al. | Temporal and spatial distribution changing characteristics of exogenous pollution load into Dianchi Lake, Southwest of China | |
Li et al. | Assessing the effectiveness and cost efficiency of green infrastructure practices on surface runoff reduction at an urban watershed in China | |
CN103268424B (zh) | 地表饮用水源集水区生态监测方法 | |
Noyola et al. | Methane correction factors for estimating emissions from aerobic wastewater treatment facilities based on field data in Mexico and on literature review | |
CN204675969U (zh) | 一种多级跌水充氧式沟渠连通生态塘农村污水处理*** | |
Negi et al. | Analysing water-energy-GHG nexus in a wastewater treatment plant of Mumbai Metropolitan Region, India | |
Song et al. | Analysis of potential risks associated with urban stormwater quality for managed aquifer recharge | |
CN111810849A (zh) | 一种雨污管网零直排溯源监测方法 | |
CN115409377A (zh) | 一种新建污水处理厂可持续性评价方法及装置 | |
Zhao et al. | Role of Low-Impact Development in generation and control of urban diffuse pollution in a pilot sponge city: A paired-catchment study | |
CN115392757A (zh) | 一种城镇污水处理设施碳排放量化评价方法及装置 | |
CN115759841B (zh) | 一种基于统计年鉴面板数据的区域下水道甲烷计算方法 | |
Neo et al. | Evaluation of active, beautiful, clean waters design features in tropical urban cities: A case study in Singapore | |
Li et al. | Field study of the road stormwater runoff bioretention system with combined soil filter media and soil moisture conservation ropes in North China | |
Zhang et al. | Impacts of rainstorm characteristics on runoff quantity and quality control performance considering integrated green infrastructures | |
Zúñiga-Estrada et al. | Preliminary model-based evaluation of water conservation strategies in a semi-arid urban zone | |
CN114897233A (zh) | 基于空间聚类改进算法的农村生活污水治理模式优化方法 | |
Giwa et al. | An extensive analysis of the engineering design of underground sewage plants in China | |
Indriyani et al. | Impacts of Land Use and Land Use Change In River Basin to Water Quality of Cirarab River, Indonesia | |
Mermillod-Blondin et al. | Stormwater infiltration practices in urban area: influences on groundwater ecosystems |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Zhu Hongtao Inventor after: Wang Jian Inventor before: Zhu Hongtao Inventor before: Wang Jian |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |