CN109583791B - 分湖控藻绩效评估法 - Google Patents

分湖控藻绩效评估法 Download PDF

Info

Publication number
CN109583791B
CN109583791B CN201811590581.3A CN201811590581A CN109583791B CN 109583791 B CN109583791 B CN 109583791B CN 201811590581 A CN201811590581 A CN 201811590581A CN 109583791 B CN109583791 B CN 109583791B
Authority
CN
China
Prior art keywords
point
value
time point
trajectory
filtering
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
Application number
CN201811590581.3A
Other languages
English (en)
Other versions
CN109583791A (zh
Inventor
赵健
富国
郝晨林
徐香勤
韩雪梅
刘庆庆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chinese Research Academy of Environmental Sciences
Original Assignee
Chinese Research Academy of Environmental Sciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chinese Research Academy of Environmental Sciences filed Critical Chinese Research Academy of Environmental Sciences
Publication of CN109583791A publication Critical patent/CN109583791A/zh
Application granted granted Critical
Publication of CN109583791B publication Critical patent/CN109583791B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152—Water filtration

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Cultivation Of Seaweed (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于滤波轨线的分湖控藻绩效评估法，包含以下步骤：均值滤波轨线建立步骤：采用滤波方法建立CHLa值与限制营养物浓度的均值滤波轨线；等效率轨线建立步骤：在均值滤波轨线上寻找等效率轨线，并根据等效率轨线建立等产藻效率值的等效率轨线；极值射线建立步骤：根据等效率轨线在均值滤波轨线上寻找拐点建立至少一条过零点的极值射线；效益值获得步骤：在均值滤波轨线上定义起点时间点及终点时间点，并通过起点时间点或终点时间点的营养物浓度值与等效率轨线及极值射线获得评估参考点，根据起点时间点、终点时间点及评估参考点获得视在效益值、实在效益值及潜在效益值。

Description

分湖控藻绩效评估法
技术领域
本发明属于环境保护与管理领域，具体地说，尤其涉及基于滤波轨线的一种分湖控藻绩效评估法。
背景技术
目前的湖泊富营养化评价大致有达标评价、营养状态评估、藻华灾害程度评估三个层面：
1)地表水达标评价
以TP为主、以TN为辅的一主一辅的氮磷双控考核，认为其可以反映藻华(CHLa)的变化。
2)富营养化状态评估
以CHLa为主，以TP、TN、高锰酸盐指数和透明度为辅的一主四辅的加权平均(综合营养指数TLIΣ)评价。
3)藻华灾害程度评估
目前太湖等湖泊有藻华的灾害等级评价，强调范围及强度，以及年度评价中增加的频次。即考虑藻华密度(CHLa)的年内峰值、分布面积等。
但在实际实施中，发现存在着以下缺陷：
1)地表水达标评价
缺点是：(1)全国一标，必然会存在较大范围的过保护及欠保护情况，5个分类虽降低过保护及欠保护情况。(2)不考核CHLa，在背景环境因子出现趋势变化时，对很多湖泊的富营养化状态产生误判。如水资源枯竭、气候变暖等趋势下，欠保护程度及范围会逐步加大。例如，“依法被污染”的情况-氮磷达标，藻华严重。
2)富营养化状态评估
缺点是：(1)全国一标，且采用的是回归线控制，约50％的保证率，背景条件平稳时过保护及欠保护大致相当，严格的讲欠保护范围可达40-50％；背景条件恶化时过欠保护增加，欠保护范围超过40-50％；。(2)由于辅助指标过多，任何一个辅助指标不是限制因子时，均对富营养化水平及程度产生一定影响，而多个辅助指标不是限制因子时，影响加大。例如，“被依法污染”的情况-依CHLa评价持续贫营养状态，TN持续劣V类重富营养，综合评价可能是轻、中富营养化，但实际贫营养化被认为超载。(3)忽视单指标的TLI的定位功能，如TLI(CHLa)指标最大时，实际的富营养化危害具有决定性以包含辅助指标的作用及其他作用，计算中以四个辅助指标的影响权重降低富营养化危害程度的评估是不合理的，当且仅当主控指标低才具有一定的合理性。(4)TLI方法仅适用于CHLa与营养物存在协调关系(正相关)时段可用，而目前许多湖泊处在非典型的非协调段，并不适用于改方法。
3)藻华灾害程度评估
缺点是：目前标准中对TP、TN并无对应的时间保证率的明确要求要求及对应藻类密度的控制要求，导致即使TP、TN100％达标，也可能对应60-90％的藻华爆发。
因此急需开发一种克服上述缺陷的分湖控藻绩效评估法。
发明内容
为了克服上述现有技术存在的问题，本发明的目的在于提供一种基于滤波轨线的分湖控藻绩效评估法，其中，包含以下步骤：
均值滤波轨线建立步骤：采用滤波方法建立CHLa值与限制营养物浓度的均值滤波轨线；
等效射线建立步骤：在所述均值滤波轨线上寻找等效率轨线，并根据所述等效率轨线建立等产藻效率值的等效射线；
极值射线建立步骤：根据所述等效射线在所述均值滤波轨线上寻找拐点建立至少一条过零点的极值射线；
效益值获得步骤：在所述均值滤波轨线上定义起点时间点及终点时间点，并通过所述起点时间点或所述终点时间点的营养物浓度值与所述等效射线及所述极值射线获得评估参考点，根据所述起点时间点、所述终点时间点及所述评估参考点获得视在效益值、实在效益值及潜在效益值。
上述的分湖控藻绩效评估法，其中，于所述效益值获得步骤中包括：
定义时间点步骤：在所述均值滤波轨线上定义起点时间点及终点时间点；
评估参考点获得步骤：通过所述起点时间点或所述终点时间点的营养物浓度值与所述等效射线及所述极值射线的最近或最适交点获得评估参考点；
效益值获得步骤：根据所述起点时间点、所述终点时间点及所述评估参考点的CHLa值获得视在效益值、实在效益值及潜在效益值。
上述的分湖控藻绩效评估法，其中，于所述均值滤波轨线建立步骤中所述均值滤波轨线为年度均值滤波轨线，所述起点时间点为起点年份，所述终点时间点为终点年份。
上述的分湖控藻绩效评估法，其中，于所述等效射线建立步骤中所述等效率轨线为与过零点的所述极值射线重合较多的轨线。
上述的分湖控藻绩效评估法，其中，于所述等效射线建立步骤中所述等效射线为多条。
上述的分湖控藻绩效评估法，其中，于所述极值射线建立步骤中，所述拐点为产藻效率最高或最低的拐点。
本发明还提供一种基于滤波轨线的分湖控藻绩效评估法，其中，包含以下步骤：
均值滤波轨线建立步骤：采用滤波方法建立CHLa值与限制营养物浓度的均值滤波轨线；
趋势外推线建立步骤：在所述均值滤波轨线上寻找拐点，并根据所述拐点建立趋势外推线；
效益值获得步骤：在所述均值滤波轨线上定义起点时间点及终点时间点，并通过所述起点时间点或所述终点时间点的营养物浓度值与所述趋势外推线获得评估参考点，根据所述起点时间点、所述终点时间点及所述评估参考点获得视在效益值、实在效益值及潜在效益值。
上述的分湖控藻绩效评估法，其中，于所述效益值获得步骤中包括：
定义时间点步骤：在所述均值滤波轨线上定义起点时间点及终点时间点；
评估参考点获得步骤：通过所述起点时间点或所述终点时间点的营养物浓度值与所述趋势外推线的最近或最适交点获得评估参考点；
效益值获得步骤：根据所述起点时间点、所述终点时间点及所述评估参考点的CHLa值获得视在效益值、实在效益值及潜在效益值。
上述的分湖控藻绩效评估法，其中，于所述趋势外推线建立步骤中包括：
分段步骤：根据所述拐点对所述均值滤波轨线进行分段；
趋势外推线获得步骤：在所有分段选择CHLa-限制性营养物协调分段建立所述趋势外推线。
上述的分湖控藻绩效评估法，其中，于所述均值滤波轨线建立步骤中所述均值滤波轨线为年度均值滤波轨线，所述起点时间点为起点年份，所述终点时间点为终点年份。
针对于现有技术本发明的功效在于：基于滤波轨线的分湖控藻绩效评估法，可评估治理项目和背景条件对控藻过程的贡献，采用视在效率(治理过程前后叶绿素a浓度的减量)、实在效率(治理项目措施减量)、潜在效率(背景条件影响减量)的解析方法，可以分析背景条件发生趋势变化时，如何从视在效率中看到减污工程的实在效率，可以识别实在效率与潜在效率的相对变化，从而不会因视在效率低而否认减污过程的贡献；也不会因视在效率高而夸大减污过程的贡献，从而避免误判投资效益；而潜在效率的持续为负，往往会促进治理投资向湖泊生境的恢复及保护倾斜。合理的控藻绩效评估可为富营养化湖泊恢复措施的调整提供技术支持，以利于实现精准的水质目标过程管理。
附图说明
图1为本发明分湖控藻绩效评估法第一实施例的流程图；
图2为图1的分步骤流程图；
图3为基于图1的控藻绩效评估图；
图4为本发明分湖控藻绩效评估法第二实施例的流程图；
图5-6为图4的分步骤流程图；
图7为基于图4的第二实施例的控藻绩效评估图；
图8为基于图4的第三实施例的控藻绩效评估图；
图9为案例一和二的控藻绩效评估图；
图10为案例三的控藻绩效评估图。
具体实施方式
下面结合附图与具体实施例对本发明作进一步详细描述：本实施例在以本发明技术方案为前提下进行实施，给出了实施方式和操作过程，但本发明的保护范围不限于下述的实施例。
为更好地阐述本发明，以下将涉及本发明的相关名词进行相应地解释。
视在效益BS：治理项目前后CHLa浓度的减量。直接体现的直观或表观效果。综合影响可由湖内的监控藻类初级生产力指标CHLa和营养物浓度指标TP(TN)趋势变化表达出来。
实在效益BR：减污项目对CHLa浓度的减量贡献。治理措施产生的效益，即项目实施后的自身达到的效果(不是指设计效益。这里不讨论设计效益与实施效益之差及原因)。背景条件无趋势性变化时，有无治理项目的CHLa浓度的减量改变。
潜在效益BP：背景条件对CHLa的减量贡献。为视在效益与实在效益之差,体现背景条件的贡献或治理措施之外的其他贡献。当治理期间背景条件过程平稳，则视在效益与实在效益等同。无治理项目时，背景条件趋势性变化导致CHLa浓度的减量改变。
与效益对应的效率也可同样定义，将效益除以初值，可得视在效率ES、实在效率ER和潜在效率EP。
CHLa为湖水中叶绿素a含量(mg/m3)；
TP为湖水中总磷浓度(mg/m3)。
以下请参照图1-3，图1为本发明分湖控藻绩效评估法第一实施例的流程图；图2为图1的分步骤流程图；图3为基于图1的控藻绩效评估图。如图1-3所示，本发明分湖控藻绩效评估法包括：
均值滤波轨线建立步骤S11：采用滤波方法建立CHLa值与限制营养物浓度的均值滤波轨线L0，其中所述均值滤波轨线为年度均值滤波轨线，所述起点时间点为起点年份，所述终点时间点为终点年份，但本发明并不以此为限；
等效射线建立步骤S12：在所述均值滤波轨线上寻找等效率轨线，并根据等效率轨线等产藻效率值的等效射线，其中所述等效射线为与过零点的极值射线重合较多的轨线，所述等效射线可以为多条；
极值射线建立步骤S13：根据所述等效射线在所述均值滤波轨线上寻找拐点建立至少一条过零点的极值射线，所述拐点为产藻效率最高或最低的拐点；
效益值获得步骤S14：在所述均值滤波轨线上定义起点时间点及终点时间点，并通过所述起点时间点或所述终点时间点的营养物浓度值与所述等效射线及所述极值射线获得评估参考点，根据所述起点时间点、所述终点时间点及所述评估参考点获得视在效益值、实在效益值及潜在效益值。
进一步地，于所述效益值获得步骤S14中包括：
定义时间点步骤S141：在所述均值滤波轨线上定义起点时间点及终点时间点；
评估参考点获得步骤S142：通过所述起点时间点或所述终点时间点的营养物浓度值与所述等效射线及所述极值射线的最近或最适交点获得评估参考点；
效益值获得步骤S143：根据所述起点时间点、所述终点时间点及所述评估参考点的CHLa值获得视在效益值、实在效益值及潜在效益值，进而分析影响三效比例的原因，评估治理过程的得失，提出改进的对策。
以下结合图3进一步说明本发明分湖控藻绩效评估法。如图3所示，本实施例基于零点射线法进行三效分解，由于CHLa-TP总体非协调，加三条过零点的极值射线L1、L2、L3作为降藻效率评估的辅助线。根据设计的起点时间点S和评价终点时间点E，利用等值TP得到起点评估参考点Rs和终点评估参考点Re。表1和表2分别为根据终点估参考点和起点评估参考点估算的三效估值。可以看到治理的实在效益为正，但潜在效益超过实在效益导致视在效益为负。潜在效益的来源可以认为包括气候变暖(各种背景气象因子变化产藻效率提升)、运行水位抬高(强化藻型稳态不利大型水生植物恢复)等多种因素。本例假设终点时间点延伸到终点评估参考点产藻效率基本稳定，采用终点评估参考点评估可体现不治污大致可以达到CHLa高点(重富营养化)；假设起点时间点尚未进入稳固的藻型稳态，采用起点评估参考点可体现背景条件维持早期水平大致可以达到的治理后CHLa低点(轻富营养)。表3为两者效益的平均结果，从定性准确的角度看更合理一些。
表1控藻治理效率评估—参考终点
基础数据 标识点 年份 CHLa(μg/L)
起点 S 1992 22.6
终点 E 2011 41.5
参照终点 Re 2011 78.8
效率指标 CHLa减量算式 CHLa降效(％) CHLa(μg/L)减量
视在效率-效益 CHLa(S)-CHLa(E) -84％ -18.9
实在效率-效益 CHLa(Re)-CHLa(E) 165％ 37.3
潜在效率-效益 CHLa(S)-CHLa(Re) -249％ -56.2
表2控藻治理效率评估—参考起点
Figure BDA0001920153140000071
表3控藻治理效率评估平均结果
效率指标 CHLa降效(％) CHLa(μg/L)减量
视在效率-效益 -84％ -18.9
实在效率-效益 106％ 24.0
潜在效率-效益 -190％ -42.9
以下请参照图4-6，图4为本发明分湖控藻绩效评估法第二实施例的流程图；图5-6为图4的分步骤流程图。如图4-6所示，本发明分湖控藻绩效评估法包括：
均值滤波轨线建立步骤S21：采用滤波方法建立CHLa值与限制营养物浓度的均值滤波轨线L0，其中所述均值滤波轨线为年度均值滤波轨线，所述起点时间点为起点年份，所述终点时间点为终点年份，但本发明并不以此为限；
趋势外推线建立步骤S22：在所述均值滤波轨线上寻找拐点，并根据所述拐点建立趋势外推线L4；
效益值获得步骤S23：在所述均值滤波轨线上定义起点时间点及终点时间点，并通过所述起点时间点或所述终点时间点的营养物浓度值与所述趋势外推线获得评估参考点，根据所述起点时间点、所述终点时间点及所述评估参考点获得视在效益值、实在效益值及潜在效益值。
进一步地，于所述趋势外推线建立步骤S22中包括：
分段步骤S221：根据所述拐点对所述均值滤波轨线进行分段；
趋势外推线获得步骤S222：在所有分段选择CHLa-限制性营养物协调分段建立所述趋势外推线。
更进一步地，于所述效益值获得步骤S23中包括：
定义时间点步骤S231：在所述均值滤波轨线上定义起点时间点及终点时间点；
评估参考点获得步骤S232：通过所述起点时间点或所述终点时间点的营养物浓度值与所述趋势外推线的最近或最适交点获得评估参考点；
效益值获得步骤S233：根据所述起点时间点、所述终点时间点及所述评估参考点的CHLa值获得视在效益值、实在效益值及潜在效益值，进而分析影响三效比例的原因，评估治理过程的得失，提出改进的对策。
以下结合图7进一步说明本发明分湖控藻绩效评估法，图7为基于图4的第二实施例的控藻绩效评估图。如图7所示，本实施例基于单调轨线延伸法，即采用实际分段滤波轨线得到参考点的方法，并利用单调轨线确定参考点的方法。若评价点不在轨线数据区间内，用实际分段滤波轨线向外延伸得到参考点。单调轨线延伸法假设：在滤波轨线的拐点之间一般为单调曲线，如果该线是协调的，假定背景条件平稳，则该线的外推是合理的，因此，通过该线上的点做回归线作为终点线的势轨线，从而根据起点确定参考点。单调轨线延伸法则不适合非协调段。
由于CHLa-TP总体协调，四段轨线及其延伸线可作为降藻效率评估的趋势外推线。根据设计的起点时间点S和终点时间点E，利用等值TP得到起点评估参考点Rs和终点评估参考点Re。在本实施例中，为根据终点评估参考点Re估算三效估值的算例，其算法与前述实施例大致相同，因此在此就不再赘述了。可以看到治理的实在效益为正，但潜在效益超过实在效益导致视在效益为负。潜在效益的来源可以认为包括气候变暖(各种背景气象因子变化产藻效率提升)、运行水位抬高(强化藻型稳态不利大型水生植物恢复)、入湖水量减少(产藻效率提升)等多种因素。从近期数据来看，近期TP总量减排及补水工程对上述负效应有抵消作用，保持入湖水量对于降藻更为有效。
表4控藻治理效率评估—参考终点
Figure BDA0001920153140000091
图8为基于图4的第三实施例的控藻绩效评估图。如图8所示，在本实施例中，为根据起点评估参考点Rs估算三效估值的算例，由于CHLa-TP总体协调，四段轨线及其延伸线可作为降藻效率评估的辅助线。本实施例根据设计的评价起点S和评价终点E，利用等值TP得到参考起点Rs。
表5控藻治理效率评估—参考起点
基础数据 标识点 滤波中点年份 CHLa(μg/L)
起点 S 1991 20
终点 E 1996 60
参照起点 Rs 1996 30
效率指标 CHLa减量算式 CHLa降效(％) CHLa(μg/L)减量
视在效率-效益 CHLa(S)-CHLa(E) -200％ -40
实在效率-效益 CHLa(S)-CHLa(Rs) -150％ -30
潜在效率-效益 CHLa(Rs)-CHLa(E) -50％ -10
表5为根据起点评估参考点Rs估算三效估值算例。这是污染早期的算例，可以看到三效均为负值，其中增污的实在效率为-50％，但潜在效益为-150％，为实在效益的3倍，即非污染损害远超直接污染损害。可以认为潜在效益的包括背景条件变化气候变暖影响以及与直接污染损害间接导致的草藻迁跃及沉水植物消亡，磷的产藻效率的提高，污染加1倍，损害翻两番。所以，浅水湖泊的保护应绝对避免“先污染后治理”的方式，否则将进入“越治越污”的极难恢复的轨道；而滞留时间长的浅水湖泊更应该避让为先，不让超背景的增量氮磷负荷入湖。
请参照图9及图10，图9为案例一和二的控藻绩效评估图；图10为案例三的控藻绩效评估图；以下结合图9-10具体说明三效分解方法及关系：
案例一：减污藻增升轨案例
起点时间点S、终点时间点E和终点评估参考点Re位置见图6的标记，并依据表7进行计算设治理项目前后背景条件由低势轨线L5-中势轨线L6-高势轨线L7。
视在效益：BS＝CHLa(S)-CHLa(E)＝40-55＝-15μg/L
实在效益：BR＝CHLa(R)-CHLa(E)＝90-55＝35μg/L
潜在效益：BP＝CHLa(S)-CHLa(R)＝40-90＝-50μg/L
可以看到视在效率为负，但减污效率为正。减污的成效被背景条件的变化淹没。
表7减污升轨降藻效率(效益)计算表
基础数据 标识点 年份 CHLa(μg/L)
起点 S 1999 40
终点 E 2007 55
参照终点 R 2007 90
效率指标 CHLa减量算式 CHLa降效(％) CHLa(μg/L)减量
视在效率-效益 CHLa(S)-CHLa(E) -38％ -15
实在效率-效益 CHLa(R)-CHLa(E) 88％ 35
潜在效率-效益 CHLa(S)-CHLa(R) -125％ -50
案例二：增污藻减降轨案例
起点时间点S、终点时间点E和终点评估参考点Re位置见图6的标记，设增污前后背景条件由高势轨线L7-中势轨线L6-低势轨线L5。该例的三效评估结果见表8。可以看到视在效率为正，但减污效率为负。增污负效益被背景条件的正效益淹没。
表8增污升轨降藻效率(效益)计算表
基础数据 标识点 年份 CHLa(μg/L)
起点 S 1999 55
终点 E 2007 40
参照终点 R 2007 30
效率指标 CHLa减量算式 CHLa降效(％) CHLa(μg/L)减量
视在效率-效益 CHLa(S)-CHLa(E) 27％ 15
实在效率-效益 CHLa(R)-CHLa(E) -18％ -10
潜在效率-效益 CHLa(S)-CHLa(R) 45％ 25
案例三：减污藻减升轨案例
起点时间点S、终点时间点E和终点评估参考点Re位置见图7的标记，设治理项目前后背景条件由低势轨线L5-中势轨线L6-高势轨线L7。该例的三效评估结果见表9。可以看到视在效率较小，减污效率很大，为视在效率的6倍。若不实施减污，CHLa会增加2倍；实施减污，CHLa降低25％。
表9增污升轨降藻效率(效益)计算表
基础数据 标识点 年份 CHLa(μg/L)
起点 S 2000 110
终点 E 2010 100
参照终点 R 2010 155
效率指标 CHLa减量算式 CHLa降效(％) CHLa(μg/L)减量
视在效率-效益 CHLa(S)-CHLa(E) 9％ 10
实在效率-效益 CHLa(R)-CHLa(E) 50％ 55
潜在效率-效益 CHLa(S)-CHLa(R) -41％ -45
从以上三例可以看出，在背景条件发生趋势变化时，可从视在效率中看到控藻工程项目的实在效率，能有效识别实在效率与潜在效率的相对变化，减少效率评估的误判。
综上所述，本发明采用分湖年尺度的CHLa和TP数据，基于低通滤波轨线识别背景条件的趋势变化引起的CHLa-TP响应关系变化，进行控藻绩效三效分解，同时首创零点射线法、单调轨线延伸法应用于控藻绩效评估。因此基于滤波轨线的分湖控藻绩效评估法，可评估治理项目和背景条件对控藻过程的贡献，采用视在效率(治理过程前后叶绿素a浓度的减量)、实在效率(治理项目措施减量)、潜在效率(背景条件影响减量)的解析方法，可以分析背景条件发生趋势变化时，如何从视在效率中看到减污工程的实在效率，可以识别实在效率与潜在效率的相对变化，从而不会因视在效率低而否认减污过程的贡献；也不会因视在效率高而夸大减污过程的贡献，从而避免误判投资效益；而潜在效率的持续为负，往往会促进治理投资向湖泊生境的恢复及保护倾斜。合理的控藻绩效评估可为富营养化湖泊恢复措施的调整提供技术支持，以利于实现精准的水质目标过程管理。
需要说明的是：以上实施例仅仅用以说明本发明，而并非限制本发明所描述的技术方案；同时，尽管本说明书参照上述实施例对本发明进行了详细的说明，但是，本领域的普通技术人员应当理解，仍然可以对本发明进行修改或等同替换；因此，一切不脱离本发明的精神和范围的技术方案及其改进，均应涵盖在本发明所附权利要求的保护范围之内。

Claims (9)

1.一种基于滤波轨线的分湖控藻绩效评估法，其特征在于，包含以下步骤：
均值滤波轨线建立步骤：采用滤波方法建立CHLa值与限制营养物浓度的均值滤波轨线；
等效射线建立步骤：在所述均值滤波轨线上寻找等效率轨线，并根据所述等效率轨线建立等产藻效率值的等效射线；
极值射线建立步骤：根据所述等效射线在所述均值滤波轨线上寻找拐点建立至少一条过零点的极值射线；
所述等效率轨线为与过零点的所述极值射线重合较多的轨线；
效益值获得步骤：在所述均值滤波轨线上定义起点时间点及终点时间点，并通过所述起点时间点或所述终点时间点的营养物浓度值与所述等效射线及所述极值射线获得评估参考点，根据所述起点时间点、所述终点时间点及所述评估参考点获得视在效益值、实在效益值及潜在效益值；
所述评估参考点获得步骤：通过所述起点时间点或所述终点时间点的营养物浓度值与所述等效射线及所述极值射线的最近点获得评估参考点；
所述视在效益值为治理项目前后CHLa浓度的减量；
所述实在效益值为减污项目对CHLa浓度的减量贡献；
所述潜在效益值为背景条件对CHLa浓度的减量贡献。
2.如权利要求1所述的分湖控藻绩效评估法，其特征在于，于所述效益值获得步骤中包括：
定义时间点步骤：在所述均值滤波轨线上定义起点时间点及终点时间点；
效益值获得步骤：根据所述起点时间点、所述终点时间点及所述评估参考点的CHLa值获得视在效益值、实在效益值及潜在效益值。
3.如权利要求1所述的分湖控藻绩效评估法，其特征在于，于所述均值滤波轨线建立步骤中所述均值滤波轨线为年度均值滤波轨线，所述起点时间点为起点年份，所述终点时间点为终点年份。
4.如权利要求3所述的分湖控藻绩效评估法，其特征在于，于所述等效射线建立步骤中所述等效射线为多条。
5.如权利要求1所述的分湖控藻绩效评估法，其特征在于，于所述极值射线建立步骤中，所述拐点为产藻效率最高或最低的拐点。
6.一种基于滤波轨线的分湖控藻绩效评估法，其特征在于，包含以下步骤：
均值滤波轨线建立步骤：采用滤波方法建立CHLa值与限制营养物浓度的均值滤波轨线；
趋势外推线建立步骤：在所述均值滤波轨线上寻找拐点，并根据所述拐点建立趋势外推线；
效益值获得步骤：在所述均值滤波轨线上定义起点时间点及终点时间点，并通过所述起点时间点或所述终点时间点的营养物浓度值与所述趋势外推线获得评估参考点，根据所述起点时间点、所述终点时间点及所述评估参考点获得视在效益值、实在效益值及潜在效益值；
所述评估参考点获得步骤：通过所述起点时间点或所述终点时间点的营养物浓度值与所述趋势外推线的最近或最适交点获得评估参考点；
所述视在效益值为治理项目前后CHLa浓度的减量；
所述实在效益值为减污项目对CHLa浓度的减量贡献；
所述潜在效益值为背景条件对CHLa浓度的减量贡献。
7.如权利要求6所述的分湖控藻绩效评估法，其特征在于，于所述效益值获得步骤中包括：
定义时间点步骤：在所述均值滤波轨线上定义起点时间点及终点时间点；
效益值获得步骤：根据所述起点时间点、所述终点时间点及所述评估参考点的CHLa值获得视在效益值、实在效益值及潜在效益值。
8.如权利要求6所述的分湖控藻绩效评估法，其特征在于，于所述趋势外推线建立步骤中包括：
分段步骤：根据所述拐点对所述均值滤波轨线进行分段；
趋势外推线获得步骤：在所有分段选择CHLa-限制性营养物协调分段建立所述趋势外推线。
9.如权利要求6所述的分湖控藻绩效评估法，其特征在于，于所述均值滤波轨线建立步骤中所述均值滤波轨线为年度均值滤波轨线，所述起点时间点为起点年份，所述终点时间点为终点年份。
CN201811590581.3A 2018-05-22 2018-12-20 分湖控藻绩效评估法 Active CN109583791B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2018104971046 2018-05-22
CN201810497104 2018-05-22

Publications (2)

Publication Number Publication Date
CN109583791A CN109583791A (zh) 2019-04-05
CN109583791B true CN109583791B (zh) 2023-04-14

Family

ID=65932484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811590581.3A Active CN109583791B (zh) 2018-05-22 2018-12-20 分湖控藻绩效评估法

Country Status (1)

Country Link
CN (1) CN109583791B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102053140A (zh) * 2009-10-27 2011-05-11 华中科技大学 一种富营养化浅水湖泊生态修复后水质的诊断和评价方法
CN102628856A (zh) * 2012-02-24 2012-08-08 中国环境科学研究院 一种湖泊营养物基准参照状态确定的组合方法
WO2015109957A1 (zh) * 2014-01-23 2015-07-30 深圳市城市规划设计研究院有限公司 一种规划人工水体富营养化预警分析方法
CN106501196A (zh) * 2016-10-13 2017-03-15 中国科学院南京地理与湖泊研究所 一种基于水体吸收系数的湖泊富营养化评价方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7720615B2 (en) * 2006-03-20 2010-05-18 Sensis Corporation System for detection and prediction of water quality events

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102053140A (zh) * 2009-10-27 2011-05-11 华中科技大学 一种富营养化浅水湖泊生态修复后水质的诊断和评价方法
CN102628856A (zh) * 2012-02-24 2012-08-08 中国环境科学研究院 一种湖泊营养物基准参照状态确定的组合方法
WO2015109957A1 (zh) * 2014-01-23 2015-07-30 深圳市城市规划设计研究院有限公司 一种规划人工水体富营养化预警分析方法
CN106501196A (zh) * 2016-10-13 2017-03-15 中国科学院南京地理与湖泊研究所 一种基于水体吸收系数的湖泊富营养化评价方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
变权综合营养状态指数(TLI_(CW))评价方法及其在三峡库区支流龙河的应用;韩雪梅 等;《环境工程技术学报》;20181120(第06期);全文 *

Also Published As

Publication number Publication date
CN109583791A (zh) 2019-04-05

Similar Documents

Publication Publication Date Title
CN110134919B (zh) 风电机组异常数据清洗方法
CN106204491A (zh) 一种基于暗通道先验的自适应图像去雾方法
CN112434962B (zh) 基于电力负荷数据的企业用户状态评估方法和系统
CN103576060A (zh) 基于小波自适应阈值的局部放电信号去噪方法
CN110009145B (zh) 基于时频变换特征提取及自回归趋势预测技术的电力负荷预测方法及装置
CN107067375A (zh) 一种基于暗通道先验和边缘信息的图像去雾方法
CN112730938A (zh) 一种基于用电采集大数据的窃电用户判断方法
CN106683405B (zh) 基于道路车流密度的道路交通状态评价方法
CN109583791B (zh) 分湖控藻绩效评估法
CN110212592B (zh) 基于分段线性表达的火电机组负荷调节最大速率估计方法及系统
CN114756817A (zh) 一种基于Copula函数的复合洪涝灾害联合概率分析方法
CN116796403A (zh) 一种基于商业建筑综合能耗预测的建筑节能方法
CN115968088A (zh) 隧道智慧调光方法、系统及计算机储存介质
CN111898097B (zh) 结合概率密度和保证率的生态流量确定方法
CN105414774A (zh) 一种能够自主切割的激光切割装置
CN108628960A (zh) 一种专变电量数据异常快速诊断流程方法
CN105677770B (zh) 近海海洋环境数据监测的自适应采样方法
CN114897918B (zh) 基于人工智能的光伏清洁机器人毛刷功率调整方法及系统
CN107203837B (zh) 基于负荷率和同时率的区域整体设备经济寿命测算方法
CN103761525A (zh) 基于显微ct数字图像处理铝电解用阴极炭块组分检测方法
CN108983301A (zh) 一种k指数现报方法
CN111985563B (zh) 一种基于谐波发射水平的工业用户电力设备分类方法
CN111242892B (zh) 一种太阳能光伏电池缺陷的检测方法
CN114912664A (zh) 基于水质-水位二元响应关系的适宜生态水位预测方法
CN110117798B (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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zhao Jian

Inventor after: Fu Guo

Inventor after: Hao Chenlin

Inventor after: Xu Xiangqin

Inventor after: Han Xuemei

Inventor after: Liu Qingqing

Inventor before: Fu Guo

Inventor before: Hao Chenlin

Inventor before: Zhao Jian

Inventor before: Xu Xiangqin

Inventor before: Han Xuemei

Inventor before: Liu Qingqing

GR01 Patent grant
GR01 Patent grant