CN107273667A - 一种工业废水中污染物排放量驱动因子识别方法 - Google Patents
一种工业废水中污染物排放量驱动因子识别方法 Download PDFInfo
- Publication number
- CN107273667A CN107273667A CN201710387636.XA CN201710387636A CN107273667A CN 107273667 A CN107273667 A CN 107273667A CN 201710387636 A CN201710387636 A CN 201710387636A CN 107273667 A CN107273667 A CN 107273667A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- delta
- mtr
- mtd
- 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
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2219/00—Indexing scheme relating to application aspects of data processing equipment or methods
- G06F2219/10—Environmental application, e.g. waste reduction, pollution control, compliance with environmental legislation
Landscapes
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Treatment Of Water By Ion Exchange (AREA)
Abstract
本发明公开了一种考虑节水减排和污染物治理作用的工业废水中污染物排放量驱动因子识别方法，其步骤包括：收集不同年份可比价的工业分行业增加值、新鲜用水量、废水排放量和废水中污染物排放数据；从污染物产生、治理与排放全过程出发，考虑节水减排和节水减污作用，构造了工业废水中污染物排放量恒等式，将工业废水中污染物排放量变化分解成工业增加值、排污结构、用水效率、废水排放系数、污染物产生强度和污染物治理6个驱动因子；利用LMDI方法，推导贡献值测算公式，测算驱动因子贡献与节水减排、节水减污作用，无残差、定量解释污染物排放量变化贡献，并通过贡献值大小，识别关键因子，揭示驱动机制。
Description
技术领域
本发明涉及一种考虑节水减排和污染物治理作用的工业废水中污染物排放量驱动因子识别方法。
背景技术
工业废水及其污染物排放是造成我国水体污染严重因素之一。依据1997～2015年《中国环境统计年报》中点源废水和主要污染物数据发现，废水排放量逐年增长，由1997年的415.7亿吨增加到2015年的734.7亿吨；其中工业废水排放量明显下降，由1997年的226.7亿吨降至2015年199.5亿吨。而点源化学需氧量(COD)排放量和氨氮排放量逐年降低。COD排放量由1997年的1757.0万吨降低到2015年的1140.4万吨，其中工业废水排放量明显下降，由1997年的226.7亿吨降至2015年199.5亿吨。虽然工业废水及其污染物排放形势有所好转，但当前，我国一些地区水环境质量差、水生态受损重、环境隐患多等问题十分突出，影响和损害群众健康，不利于经济社会持续发展。为切实加大水污染防治力度，保障国家水安全，国务院关于印发水污染防治行动计划，简称“水十条”，提出全面控制污染物排放，狠抓工业污染防治。识别影响工业废水与污染物排放驱动因子，揭示驱动机制，对于掌握未来工业废水与污染物排放变化趋势，提出合理可行的污染物控制和污染防治策略，具有重要意义。
现有研究没有考虑节水减污作用，无法测算节水对于减污贡献。本发明从污染物产生、治理与排放全过程出发，考虑节水减污作用，构造了工业废水中污染物排放量恒等式，将工业废水中污染物排放量变化分解成工业增加值、排污结构、用水效率、废水排放系数、污染物产生强度和污染物治理6个驱动因子；利用对数平均迪氏指数法(logarithmicmean divisia index，LMDI)方法，推导贡献值测算公式，定量测算驱动因子贡献与节水减排、节水减污作用，无残差、定量解释污染物排放量变化贡献；通过贡献值大小，识别关键因子，揭示驱动机制。
发明内容
本发明的目的在于提供一种考虑节水减排和污染物治理作用的工业废水中污染物排放量驱动因子识别方法。
本发明的目的通过以下技术方案实现，其步骤包括：
1、一种考虑节水减排和污染物治理作用的工业废水中污染物排放量驱动因子识别方法，其步骤包括：
(1)收集不同年份工业分行业可比价增加值Gt i、新鲜用水量Wt i、废水排放量Pt i和工业废水中污染物产生量At i、污染物排放量Et i数据，计算得到不同年份工业总体的增加值Gt、新鲜用水量Wt、废水排放量Pt和工业废水中污染物产生量At、污染物排放量Et数据，其中
(2)建立工业废水中污染物排放量驱动因子识别模型：利用如下污染物排放量恒等式：
式中：G、W、P、A和E分别表示工业的增加值、新鲜用水量、废水排放量、污染物产生量和排放量，下标i和t分别表示工业行业和年份；将工业废水中污染物排放量变化分解成工业增加值EGa、排污结构EGs、用水效率EWe、废水排放系数EPe、污染物产生强度Eg和污染物治理Er共6个驱动因子：
式中：则t～t+1年期间工业废水排放量ΔE变化可表示成这6个驱动因子贡献值：
(3)基于LMDI方法，推导出如式(4)所示的驱动因子贡献值测算公式，利用这一公式测算每个行业6个驱动因子贡献值以及工业总体的6个驱动因子贡献值：
式中：ΔEGai、ΔEGsi、ΔEWei、ΔEPei、ΔEgi和ΔEri分别为第i个行业的排污规模、增加值比重、用水效率、废水排放系数、污染物产生强度和污染物治理因子贡献值，EGat＝Gt，
(4)按照式(5)计算工业总体的6个驱动因子的贡献率
式中：pEGa、pEGs、pEWe、pEPe、pEg和pEr分别为排污规模Ga、工业结构Gs、行业用水效率We、废水排放系数Pe、污染物产生强度Eg和污染物治理Er6个驱动因子的贡献率；
(5)按照式(6)分别计算每个行业6个驱动因子的贡献率：
式中：pEGai、pEGsi、pEWei、pEPei、pEgi和pEri分别为排污规模Ga、工业结构Gs、行业用水效率We、废水排放系数Pe、污染物产生强度Eg和污染物治理Er6个驱动因子的贡献率。
工业行业分类的依据是《国民经济行业分类与代码》(GB/4754-2011)，分析过程中具体采用工业行业分类及数目，是在收集到的增加值、废水排放、污染物产生、污染物排放等数据中行业分类交集的基础上，通过行业合并确定。
在测算出工业废水中污染物排放量6个驱动因子贡献值后，通过驱动因子贡献值和贡献率大小与正负，定量地解释任意两年间污染物排放量变化，分析驱动机制。当因子贡献值和贡献率为负值时，表示它是抑制污染物排放量增长的因子，简称抑制因子，反之则为推动污染物排放量增长的因素，简称推动因子。贡献值和贡献率的绝对值大小表示抑制和推动作用的强弱。贡献值和贡献率绝对值较大的因子是主要驱动因子。
本发明方法考虑了节水减污作用，从污染物产生、治理与排放全过程出发，建立了工业废水中污染物排放量驱动因子集；利用对数平均迪氏指数法(logarithmic meandivisia index，LMDI)方法，推导出了驱动因子贡献值测算公式，定量测算驱动因子贡献，完整解释灌溉用水量变化；通过贡献值大小，识别关键因子，揭示驱动机制。
附图说明
图1为山东省2002～2012年工业废水中COD排放量变化驱动因子识别结果。
具体实施方式
为了更好地理解本发明，下面以山东省工业废水中COD排放量为例来详细说明本发明的技术方案，但是本发明并不局限于此。
(1)收集山东省2002年、2007年和2012年24个工业行业增加值新鲜取用水量废水排放量COD产生量和COD排放量并计算山东省2002年、2007年和2012年工业总体的增加值Gt、新鲜用水量Wt、废水排放量Pt和工业废水中污染物产生量At、污染物排放量Et数据；
(2)建立工业废水中污染物排放量驱动因子识别模型：利用如下污染物排放量恒等式：
将工业废水中污染物排放量变化分解成工业增加值EGa、排污结构EGs、用水效率EWe、废水排放系数EPe、污染物产生强度Eg和污染物治理Er 6个驱动因子：
其中：于是t～t+1年期间工业废水中污染物排放量ΔE变化可表示成6个驱动因子贡献值：
(3)基于LMDI方法，推导出如式(4)所示的驱动因子贡献值测算公式，利用这一公式测算每个行业6个驱动因子贡献值以及工业总体的6个驱动因子贡献值：
其中ΔEGai、ΔEGsi、ΔEWei、ΔEPei、ΔEgi和ΔEri分别为第i个行业的排污规模、增加值比重、用水效率、废水排放系数、污染物产生强度和污染物治理因子贡献值，计算出了驱动因子贡献值。表1和图1中给出了山东省2002～2007年、2007～2012年以及2002～2012年三个时期工业总体的6个驱动因子贡献值。表2给出了山东省2002～2012年期间24个行业的6个驱动因子贡献值。
(4)按照式(5)计算6个驱动因子的贡献率
式中：pEGa、pEGs、pEWe、pEPe、pEg和pEr分别为排污规模Ga、工业结构Gs、行业用水效率We、废水排放系数Pe、污染物产生强度Eg和污染物治理E r6个驱动因子的贡献率。表3中给出了山东省2002～2007年、2007～2012年以及2002～2012年三个时期工业总体的6个驱动因子贡献率。
(5)按照式(6)分别计算每个行业6个驱动因子的贡献率：
式中：pEGai、pEGsi、pEWei、pEPei、pEgi和pEri分别为排污规模Ga、工业结构Gs、行业用水效率We、废水排放系数Pe、污染物产生强度Eg和污染物治理Er6个驱动因子的贡献率。表4给出了山东省2002～2012年期间24个行业的6个驱动因子贡献率。
表1山东省2002～2012年工业废水COD排放量变化驱动因子识别结果单位万吨
表2山东省2002～2012年工业废水COD排放量驱动因子分行业贡献值单位万吨
表3山东省2002～2012年工业废水中COD排放量变化驱动因子贡献率测算结果单位：％
表4山东省2002～2012年工业废水COD排放量驱动因子分行业贡献率单位％
在测算出工业废水中污染物排放量6个驱动因子贡献值后，通过驱动因子贡献值和贡献率大小与正负，定量地解释任意两年间污染物排放量变化，分析驱动机制。当因子贡献值和贡献率为负值时，表示它是抑制污染物排放量增长的因子，简称抑制因子，反之则为推动污染物排放量增长的因素，简称推动因子。贡献值和贡献率的绝对值大小表示抑制和推动作用的强弱。贡献值和贡献率绝对值较大的因子是主要驱动因子。
由表1和表3可知，工业行业用水效率提高、污染物产生强度降低和污染物治理力度加大，抑制了COD排放量增加，是COD排放的抑制因子；在2002～2012年期间，这三个抑制因子分别使COD排放量下降了40.55万吨、20.60万吨和19.29万吨。而排污规模扩张、废水排放率增大、工业结构不合理发展，则使得COD排放量增加，是COD排放的促进因子；在2002～2012年期间，分别使COD排放量增加了34.15万吨、20.00万吨和0.75万吨。这两类因子共同作用，使得2002年到2012年期间工业废水中COD排放量下降了25.55万吨，由39.49万吨降至13.94万吨，降幅达64.7％。从因子贡献值绝对值看，用水效率、污染物产生强度和污染物治理力度是主要的抑制因子，而排污规模和废水排放率则是主要的促进因子。工业行业内部结构调整对COD排放影响很小，贡献值仅为0.75万吨，总体偏向于高污染。山东省下阶段应加大工业结构调整力度。
由表2和表4可知，绝大部分行业的排污规模、用水效率、废水排放系数、污染物产生强度、污染物治理力度5个因子以及COD排放量的贡献值和贡献率率均为正值。这说明绝大部分行业处于排污规模(行业增加值)扩张、用水效率提高、耗水数(耗水量与新鲜用水量比值)降低、废水排放系数(废水排放量与新鲜用水量比值)提高、污染物产生强度降低、污染物治理力度加大，COD排放量降低的状态。对这5个因子24个行业贡献率从大到小排序，发现造纸印刷和文教体育用品、食品和烟草、化学产品和纺织四个高污染行业的贡献率均排在前列，累积贡献率均超过了80％，而且这四个行业COD排放量贡献率也位于前列，累积贡献率超过了85％。这说明，在2002年至2012年期间，山东省对这四个高污染行业加大了管理力度，减排效果显著。
本专利实施例中以COD排放量来表征污染物的排放，在实际应用中，也可用氨氮、总磷、有机物等排放量来表征污染物的排放，获得工业废水中氨氮、总磷或有机物排放量驱动因子贡献值和贡献率。
Claims (1)
1.一种考虑节水减排和污染物治理作用的工业废水中污染物排放量驱动因子识别方法，其步骤包括：
(1)收集不同年份工业分行业可比价增加值Gt i、新鲜用水量Wt i、废水排放量Pt i和工业废水中污染物产生量At i、污染物排放量Et i数据，计算得到不同年份工业总体的增加值Gt、新鲜用水量Wt、废水排放量Pt和工业废水中污染物产生量At、污染物排放量Et数据，其中
(2)建立工业废水中污染物排放量驱动因子识别模型：利用如下污染物排放量恒等式：
<mrow>
<msup>
<mi>E</mi>
<mi>t</mi>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mrow>
<mo>(</mo>
<msup>
<mi>G</mi>
<mi>t</mi>
</msup>
<mo>*</mo>
<mfrac>
<msubsup>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msup>
<mi>G</mi>
<mi>t</mi>
</msup>
</mfrac>
<mo>*</mo>
<mfrac>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mfrac>
<mo>*</mo>
<mfrac>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>W</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mfrac>
<mo>*</mo>
<mfrac>
<msubsup>
<mi>A</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mfrac>
<mo>*</mo>
<mfrac>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>A</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：G、W、P、A和E分别表示工业的增加值、新鲜用水量、废水排放量、污染物产生量和排放量，下标i和t分别表示工业行业和年份；将工业废水中污染物排放量变化分解成工业增加值EGa、排污结构EGs、用水效率EWe、废水排放系数EPe、污染物产生强度Eg和污染物治理Er共6个驱动因子：
<mrow>
<msup>
<mi>E</mi>
<mi>t</mi>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mrow>
<mo>(</mo>
<msup>
<mi>EGa</mi>
<mi>t</mi>
</msup>
<mo>*</mo>
<msubsup>
<mi>EGs</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>EWe</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>EPe</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>Eg</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>Er</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：EGat＝Gt，则t～t+1年期间工业废水排放量ΔE变化可表示成这6个驱动因子贡献值：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>E</mi>
<mo>=</mo>
<msup>
<mi>E</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>E</mi>
<mi>t</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>EGa</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>EGs</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>EWe</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>EPe</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>Eg</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>Er</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>EGa</mi>
<mi>t</mi>
</msup>
<mo>*</mo>
<msubsup>
<mi>EGs</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>EWe</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>EPe</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>Eg</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>*</mo>
<msubsup>
<mi>Er</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>G</mi>
<mi>a</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>G</mi>
<mi>s</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>W</mi>
<mi>e</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>P</mi>
<mi>e</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>g</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>r</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
(3)基于LMDI方法，推导出如式(4)所示的驱动因子贡献值测算公式，利用这一公式测算每个行业6个驱动因子贡献值以及工业总体的6个驱动因子贡献值：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>G</mi>
<mi>a</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msub>
<mi>&Delta;EGa</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>L</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>EGa</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>/</mo>
<msup>
<mi>EGa</mi>
<mi>t</mi>
</msup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>G</mi>
<mi>s</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msub>
<mi>&Delta;EGs</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>L</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>EGs</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>EGs</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>W</mi>
<mi>e</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msub>
<mi>&Delta;EWe</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>L</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>EWe</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>EWe</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>P</mi>
<mi>e</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msub>
<mi>&Delta;EPe</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>L</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>EPe</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>EPe</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>g</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msub>
<mi>&Delta;Eg</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mi>L</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>L</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>Eg</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>Eg</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>r</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msub>
<mi>&Delta;Er</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mi>L</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>L</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>Er</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>Er</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
式中：ΔEGai、ΔEGsi、ΔEWei、ΔEPei、ΔEgi和ΔEri分别为第i个行业的排污规模、增加值比重、用水效率、废水排放系数、污染物产生强度和污染物治理因子贡献值，EGat＝Gt，
(4)按照式(5)计算工业总体的6个驱动因子的贡献率
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>p</mi>
<mi>E</mi>
<mi>G</mi>
<mi>a</mi>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>G</mi>
<mi>a</mi>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>p</mi>
<mi>E</mi>
<mi>G</mi>
<mi>s</mi>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>G</mi>
<mi>s</mi>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>p</mi>
<mi>E</mi>
<mi>W</mi>
<mi>e</mi>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>W</mi>
<mi>e</mi>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>p</mi>
<mi>E</mi>
<mi>P</mi>
<mi>e</mi>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>P</mi>
<mi>e</mi>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>p</mi>
<mi>E</mi>
<mi>g</mi>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>g</mi>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>p</mi>
<mi>E</mi>
<mi>r</mi>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>r</mi>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：pEGa、pEGs、pEWe、pEPe、pEg和pEr分别为排污规模Ga、工业结构Gs、行业用水效率We、废水排放系数Pe、污染物产生强度Eg和污染物治理Er6个驱动因子的贡献率；
(5)按照式(6)分别计算每个行业6个驱动因子的贡献率：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>pEGa</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;EGa</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>G</mi>
<mi>a</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>pEGs</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;EGs</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>G</mi>
<mi>s</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>pEWe</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;EWe</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>W</mi>
<mi>e</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>pEPe</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;EPe</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>P</mi>
<mi>e</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>pEg</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;Eg</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>g</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>pEr</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;Er</mi>
<mi>i</mi>
</msub>
<mo>/</mo>
<mi>&Delta;</mi>
<mi>E</mi>
<mi>r</mi>
<mo>*</mo>
<mn>100.0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：pEGai、pEGsi、pEWei、pEPei、pEgi和pEri分别为排污规模Ga、工业结构Gs、行业用水效率We、废水排放系数Pe、污染物产生强度Eg和污染物治理Er6个驱动因子的贡献率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710387636.XA CN107273667B (zh) | 2017-05-27 | 2017-05-27 | 一种工业废水中污染物排放量驱动因子识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710387636.XA CN107273667B (zh) | 2017-05-27 | 2017-05-27 | 一种工业废水中污染物排放量驱动因子识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107273667A true CN107273667A (zh) | 2017-10-20 |
| CN107273667B CN107273667B (zh) | 2019-07-30 |
Family
ID=60065518
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710387636.XA Active CN107273667B (zh) | 2017-05-27 | 2017-05-27 | 一种工业废水中污染物排放量驱动因子识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107273667B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109191329A (zh) * | 2018-09-14 | 2019-01-11 | 云南电网有限责任公司 | 一种基于多层lmdi的电力消费量影响因素分解方法 |
| CN112287280A (zh) * | 2020-10-10 | 2021-01-29 | 鞍钢股份有限公司 | 一种循环水系统的污染因子富集规律的计算及应用方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101866397A (zh) * | 2010-06-21 | 2010-10-20 | 南京大学 | 一种确定影响污染物排放主导因素的方法 |
| CN104978475A (zh) * | 2014-04-10 | 2015-10-14 | 鞍钢股份有限公司 | 一种钢铁联合企业废水排放源解析方法 |
| CN106600027A (zh) * | 2016-10-31 | 2017-04-26 | 上海市政工程设计研究总院（集团）有限公司 | 城市交通碳排放量测算系统和测算方法 |
-
2017
- 2017-05-27 CN CN201710387636.XA patent/CN107273667B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101866397A (zh) * | 2010-06-21 | 2010-10-20 | 南京大学 | 一种确定影响污染物排放主导因素的方法 |
| CN104978475A (zh) * | 2014-04-10 | 2015-10-14 | 鞍钢股份有限公司 | 一种钢铁联合企业废水排放源解析方法 |
| CN106600027A (zh) * | 2016-10-31 | 2017-04-26 | 上海市政工程设计研究总院（集团）有限公司 | 城市交通碳排放量测算系统和测算方法 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109191329A (zh) * | 2018-09-14 | 2019-01-11 | 云南电网有限责任公司 | 一种基于多层lmdi的电力消费量影响因素分解方法 |
| CN112287280A (zh) * | 2020-10-10 | 2021-01-29 | 鞍钢股份有限公司 | 一种循环水系统的污染因子富集规律的计算及应用方法 |
| CN112287280B (zh) * | 2020-10-10 | 2024-01-09 | 鞍钢股份有限公司 | 一种循环水系统的污染因子富集规律的计算及应用方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107273667B (zh) | 2019-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | Operational energy performance assessment system of municipal wastewater treatment plants | |
| CN101258108B (zh) | 在分级循环序列式生物反应器中处理富氮排放物的方法和设备 | |
| Delanka-Pedige et al. | Evaluating wastewater treatment infrastructure systems based on UN Sustainable Development Goals and targets | |
| Hu et al. | Greenhouse gas emission by centralized wastewater treatment plants in Chinese industrial parks: Inventory and mitigation measures | |
| CN107273667A (zh) | 一种工业废水中污染物排放量驱动因子识别方法 | |
| Remy et al. | Identifying energy and carbon footprint optimization potentials of a sludge treatment line with life cycle assessment | |
| Li et al. | Assessment of energy use and environmental impacts of wastewater treatment plants in the entire life cycle: A system meta-analysis | |
| CN112084649A (zh) | 一种基于城镇污水处理全流程环节的碳排放模型计算方法 | |
| CN110526405A (zh) | 复合人工湿地就地去除稀土矿区氨氮污染的方法及系统 | |
| CN114037583A (zh) | 一种污水处理厂尾水湿地碳排放计算方法 | |
| Zhou et al. | Study on the impact of electric power and thermal power industry of Beijing–Tianjin–Hebei region on industrial sulfur dioxide emissions—From the perspective of green technology innovation | |
| Chen et al. | The GHG mitigation opportunity of sludge management in China | |
| Teow et al. | Comparative environmental impact evaluation using life cycle assessment approach: a case study of integrated membrane-filtration system for the treatment of aerobically-digested palm oil mill effluent | |
| Whitehead et al. | The Wheal Jane wetlands model for bioremediation of acid mine drainage | |
| CN106169101A (zh) | 一种基于结构分解的能源需求预测方法 | |
| Sobańtka et al. | Extended statistical entropy analysis (eSEA) for improving the evaluation of Austrian wastewater treatment plants | |
| CN105152467B (zh) | 生活污水集约式深度去污及脱氮处理系统与方法 | |
| CN113673117A (zh) | 基于蒙特卡罗法模拟构建海油污染物减排运行成本的方法 | |
| CN105129982A (zh) | 一种处理难降解有机废水的多级潮汐流湿地系统及方法 | |
| CN102249407A (zh) | 一种高效厌氧氨氧化反应器 | |
| Garbowski et al. | Analysis and assesment of the Wastewater Treatment Plant Operation in the City of Kłodzko | |
| Quispe et al. | Optimization of the operation of a municipal wastewater treatment plant with hydrocotyle ranunculoides | |
| CN107274076B (zh) | 考虑节水减排作用的工业废水排放量驱动因子识别方法 | |
| CN104386809A (zh) | 一种好氧颗粒污泥处理养猪废水的优化方法 | |
| CN105549388A (zh) | 一种基于自适应回归核函数的污水处理过程能耗预测方法 |
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 |