Nothing Special   »   [go: up one dir, main page]

CN106295135A - 分段回归突变点确定浅水湖泊水质管理目标方法 - Google Patents

分段回归突变点确定浅水湖泊水质管理目标方法 Download PDF

Info

Publication number
CN106295135A
CN106295135A CN201610615551.8A CN201610615551A CN106295135A CN 106295135 A CN106295135 A CN 106295135A CN 201610615551 A CN201610615551 A CN 201610615551A CN 106295135 A CN106295135 A CN 106295135A
Authority
CN
China
Prior art keywords
logm
water quality
catastrophe point
concentration
total phosphorus
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
Application number
CN201610615551.8A
Other languages
English (en)
Other versions
CN106295135B (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201610615551.8A priority Critical patent/CN106295135B/zh
Publication of CN106295135A publication Critical patent/CN106295135A/zh
Application granted granted Critical
Publication of CN106295135B publication Critical patent/CN106295135B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Water Treatment By Sorption (AREA)

Abstract

本发明公开了一种分段回归突变点确定浅水湖泊水质管理目标方法,包括以下步骤:(1)搜集需要确定水质管理目标的浅水湖泊n组叶绿素a的浓度Zi,总磷浓度Xi以及总氮浓度Yi值,获得对应的值logZi、logXi及logYi;(2)初步确定总磷突变点;(3)精确确定总磷的突变点;(4)用总氮浓度代替总磷浓度,重复步骤(2)和(3)精确确定总氮的突变点;(5)确定浅水湖泊候选营养盐水质管理目标;(6)确定浅水湖泊水体用途所要求水质并与候选营养盐水质管理目标比较,最终确定具体湖泊营养盐水质管理目标。本发明完全利用了客观观测数据中自然展现的统计突变点,精确的确定水质等级区间,易于在各种受污染浅水湖泊水质管理目标的确定中推广使用。

Description

分段回归突变点确定浅水湖泊水质管理目标方法
技术领域
本发明涉及分段回归突变点确定浅水湖泊水质管理目标方法,属于湖泊领域。
背景技术
我国大部分淡水湖泊都是浅水湖泊且富营养化问题较为严重,水质趋于恶化;而氮和磷是湖泊水体中最主要的营养盐,控制湖泊特别是浅水湖泊中氮和磷的含量对于降低湖泊富营养化程度,保护湖泊水体质量具有决定性的作用。制定恰当的氮磷营养盐水质管理目标是有效保护浅水湖泊水环境质量的前提条件,对于保护浅水湖泊水环境具有十分重要的意义。目前,浅水湖泊营养盐水质目标的制定多数依靠专家和领导的经验作出决策,其人为主观性较强;这导致制定的湖泊营养盐水质管理目标时常脱离湖泊水质的实际情况。事实上,若浅水湖泊水环境管理目标过于宽松则不利于湖泊水环境保护,导致浅水湖泊进一步富营养化,水质更加恶化;若管理目标过于严格则既无必要又可能妨碍当地经济社会发展。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种分段回归突变点确定浅水湖泊水质管理目标方法,完全利用了客观观测数据中自然展现的统计突变点,精确的确定水质等级区间,易于在各种受污染浅水湖泊水质管理目标的确定中推广使用。
技术方案:为解决上述技术问题,本发明的一种分段回归突变点确定浅水湖泊水质管理目标方法,包括以下步骤:
(1)搜集需要确定水质管理目标的浅水湖泊n组叶绿素a的浓度、总磷浓度以及总氮浓度,n为大于50的整数,Zi、Xi和Yi分别表示第i组叶绿素a的浓度、总磷浓度和总氮浓度的值,统一叶绿素a的浓度,总磷浓度和总氮浓度的单位;并将叶绿素a的浓度、总磷浓度以及总氮浓度分别取自然对数e为底的对数,获得对应的值logZi、logXi及logYi
(2)以叶绿素a的浓度对数值为因变量,以总磷浓度对数值为自变量做非参数回归,并初步确定总磷突变点;
(3)以叶绿素a浓度对数值为因变量,以总磷浓度对数值为自变量做分段线性回归,用步骤(2)中确定的总磷突变点作为求解初始值,精确确定总磷突变点;
(4)用总氮浓度代替总磷浓度,重复步骤(2)和(3)精确确定总氮的突变点;
(5)根据总氮突变点和总磷突变点及湖泊所在区域的社会经济发展情况和国家水质标准,确定浅水湖泊候选营养盐水质管理目标;
(6)确定浅水湖泊水体用途所要求水质并与候选营养盐水质管理目标比较,最终确定具体湖泊营养盐水质管理目标。
作为优选,所述步骤(2)具体包括以下步骤:
a.根据湖泊总磷浓度和叶绿素a的浓度的观测值数据数量n,用h=n-1/5计算最佳带宽值h;
b.在总磷浓度对数值的取值范围内均匀取1000个数值作为自变量logM的值logMp,即logM1、logM2、logM3...logM1000
c.从logM1开始计算得到所有满足条件logXi≤logM1+3h及logXi≥logM1-3h的总磷对数值logXi,对每一个满足条件的logXi及logM1计算相应的核函数K(logXi-logM1),作为logXi在点logM1上的权重,不满足条件logXi≤logM1+3h及logXi≥logM1-3h的logXi权重为0;其中:
K ( log X i - log M 1 ) = 1 2 π h e - 1 2 h 2 ( logX i - logM 1 ) 2
d.利用已有数值求解法,对logM1依次求解两组实数a(logM1),b(logM1);具体即要求a(logM1),b(logM1)使下式最小,
Σ i = 1 n { log Z i - a ( log M 1 ) - b ( log M 1 ) log X i } 2 K ( log X i - log M 1 )
e.重复步骤c和d,依次将logM2、logM3...logM1000替换logM1计算出两组实数a(logMp),b(logMp);
f.画出回归方程logZ=a(logMp)+b(logMp)logM的曲线,其中logZ代表回归方程的因变量,logM是回归方程的自变量;再次注意a(logMp),b(logMp)是回归方程的系数,是随着自变量logM的取值logMp的不同而不同的两组实数,观察方程logZ=a(logMp)+b(logMp)logM的曲线中的转折点,并记录该变量的突变点logγ对应的值作为分段回归求解迭代的初始值,如有多个总磷突变点,取最小的总磷突变点。
作为优选,所述步骤(3)的求解方法是求解四个实数β0,β1,δ和使下式最小,求解方法可见参考文献(Ulm,1991):
其中是示性函数,也即满足下式:
上面方程可以用已有标准非线性回归求解程序迭代计算,即为所需求解的总磷精确的突变点,迭代计算初始值的好坏直接影响求解的效率和结果的准确性;将logγ代入方程作为的迭代初始值能很好的提高计算效率和结果的准确性。的初始值为logγ。
作为优选,所述步骤(5)步骤具体包括以下步骤:
a.根据湖泊所在区域的经济发展指标,确定该区域经济发展等级,我国各地区经济发展等级包括:
A.发达地区:区域人均国民生产总值超过全国人均国民生产总值的2倍,即区域人均国民生产总值大大超过全国人均国民生产总值,人民生活富足;
B.较发达地区:区域人均国民生产总值是全国人均国民生产总值的1.3倍至2倍,即区域人均国民生产总值超过全国人均国民生产总值,人民生活比较富足;
C.中等地区:区域人均国民生产总值是全国人均国民生产总值的0.8至1.3倍,即区域人均国民生产总值与全国人均国民生产总值差不多,人民生活水平中等;
D.落后地区:区域人均国民生产总值低于全国人均国民生产总值0.8倍,即区域人均国民生产总值大大低于全国人均国民生产总值,人民生活较贫困。
b.将总氮、总磷和叶绿素a的突变点值与国家《地表水环境质量标准(GB3838-2002)》比较,确定突变点所在水质等级区间;
c.发达地区应制定更严格的浅水湖泊营养盐管理目标,因此在突变点值的基础上,升高一个等级作为营养盐水质管理目标的候选;较发达地区取突变点所在水质等级区间对应的高等级作为营养盐水质管理目标的候选;中等地区取突变点所在水质等级区间对应的低等级作为营养盐水质管理目标的候选;落后地区在突变点基础上,降低一个等级作为营养盐水质管理目标的候选。
作为优选,所述步骤(6)步骤具体包括以下步骤:
a.确定湖泊水体用途及所要求水质:按照国家水质标准水体用途可分为自然保护区、饮用水源地、渔业、工业和农业五大类;依次对应于《地表水环境质量标准(GB3838-2002)》中的I、II、III、IV和V类水;
b.将候选营养盐水质管理目标与该湖泊水体用途所要求水质进行比较,最终取两者中水质较为严格的结果作为具体湖泊营养盐水质管理的目标。
在本发明中,环境学的研究指出,浅水湖泊富营养化的主要指标叶绿素a浓度与营养盐总氮、总磷浓度之间存在着密切关系;且这些环境因子彼此关系中存在突变点。本发明依据以上观点,利用分段线性回归方法精确确定总氮、总磷和叶绿素a关系中的突变点;结合国家地表水环境质量标准、浅水湖泊所在区域的社会经济发展情况及湖泊水体的具体用途客观精确确定浅水湖泊营养盐水质管理目标的方法。本发明结合湖泊当地经济社会发展水平和湖泊水体具体用途,客观精确确定浅水湖泊水质管理目标的方法。该方法能帮助环境管理部门决定正确客观的管理目标,达到既能充分利用湖泊水资源,为经济社会发展提供必要的条件又能有效保护和改善湖泊水质,降低浅水湖泊富营养化程度。
有益效果:本发明的分段回归突变点确定浅水湖泊水质管理目标方法,完全利用了客观观测数据中自然展现的统计学突变点和国家权威的标准,几乎没有人为主观干预的步骤,克服了已有方法主观性和人为干预较强的缺点,可精确的确定湖泊区域的水质标准,且综合考虑了经济社会发展,湖泊水体用途因素,易于在各种受污染浅水湖泊水质管理目标的确定中推广使用。
附图说明
图1为本发明的流程图;
图2为总磷-叶绿素a局部线性回归图;
图3为总磷-叶绿素a分段线性回归图;
图4为总氮-叶绿素a局部线性回归图;
图5为总氮-叶绿素a分段线性回归图。
具体实施方式
下面结合附图和我国东部某浅水湖泊实际观测数据对本发明作更进一步的说明。
(1)按照图1所示流程图,首先从文献及实测中获得该浅水湖泊总磷、总氮和叶绿素a的浓度值共151组;总氮和总磷浓度的单位均为毫克/升,叶绿素a浓度单位一般为微克/升,将其除以1000后统一与总氮、总磷浓度单位一致的数值。
(2)将总磷、总氮和叶绿素a的浓度取对数值,同时将n=151,代入h=n-1/5得到理论最佳带宽h=0.366,本次计算中取h=0.37。
设定在总磷对数值的取值范围内要计算的自变量点的个数,为了最后显示的局部回归线较为光滑,本次在总磷对数值的取值范围-4.6至-0.5内均匀取1000个数值作为logMp,即为logM1、logM2、logM3…logM1000,对于每一个logMp对比总磷观测值的对数logXi,所有满足条件logXi≤logMp+3h及logXi≥logMp-3h的logXi,分别计算核函数K(logXi-logMp),其中:
K ( log X i - log M p ) = 1 2 π h e - 1 2 h 2 ( logX i - logM p ) 2
求解a(logMp),b(logMp)使下式取值最小:
Σ i = 1 n { log Z i - a ( log M p ) - b ( log M p ) log X i } 2 K ( log X i - log M p )
这个求解用已有的Newton切线法完成,对于每一个logMp,在满足上式和最小的条件下得到一组a(logMp),b(logMp),如图2所示画出方程logZ=a(logMp)+b(logMp)logM的曲线,这条曲线的横纵坐标分别代表总磷浓度的对数值和对应的叶绿素a的对数值,图上有较明显的多个突变点,出于保护湖泊环境的目的,其目标浓度应该较小,所以应该取总磷值最小的突变点,在这个例子中突变点对应的总磷值约为0.045毫克/升,叶绿素a的值约为0.0037毫克/升。
(3)将上一步求得的总磷0.045毫克/升,作为突变点迭代的初始值代入中,求解β0,β1,δ和使这个和最小,具体求解方法参考Ulm,K.,1991.A statistical method for assessing a thresholdin epidemiological studies.Statistics in medicine 10,341-349;其中是示性函数,也即下式:
对应的分段线性回归方程为:
如图3,在本例中,最终求出各参数的值:β0=-6.12;β1=-0.32;δ=0.50;最重要的是精确估计的总磷突变点为0.061毫克/升。
(4)如图4和图5所示,将研究对象换为总氮-叶绿素a重复步骤(2)-(3)得到总氮突变点对应浓度为0.71毫克/升。
(5)对照《地表水环境质量标准(GB3838-2002)》可以发现总磷突变点0.061毫克/升,位于III类水和IV类水标准之间;而总氮突变点0.71毫克/升位于II类水和III类水标准之间。该浅水湖泊位于江苏南部的苏州和无锡两市辖区内,该区域人均国民生产总值远高于全国人均国民生产总值(区域人均国民生产总值超过全国人均国民生产总值的2倍,人民生活富足),人民生活富裕,是我国经济最发达的区域之一。因此,升高一个等级作为营养盐水质管理目标的候选,即总氮的候选水质管理目标定为I类水水质标准,总磷的水质管理目标定为II类水。
(6)经调查该湖泊主要的用途为农业和渔业用水,其对应的水质标准为III类水;因此,经过与候选的总氮和总磷水质管理目标比较,最终该湖泊的总氮水质管理目标是I类水水质标准,总磷的水质管理目标定为II类水水质标准。
本发明完全利用了客观观测数据中自然展现的统计学突变点和国家权威的标准,几乎没有人为主观干预的步骤,克服了已有方法主观性和人为干预较强的缺点,能够精确的确定目标区域的水质标准,便于目标区域的水质管理。该方法对所用数据无任何特殊要求,且综合考虑了经济社会发展,湖泊水体用途因素,易于在各种受污染浅水湖泊水质管理目标的确定中推广使用。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (5)

1.一种分段回归突变点确定浅水湖泊水质管理目标方法,其特征在于,包括以下步骤:
(1)搜集需要确定水质管理目标的浅水湖泊n组叶绿素a的浓度、总磷浓度以及总氮浓度,n为大于50的整数,Zi、Xi和Yi分别表示第i组叶绿素a的浓度、总磷浓度和总氮浓度的值,统一叶绿素a的浓度、总磷浓度和总氮浓度的单位;并将叶绿素a的浓度、总磷浓度以及总氮浓度分别取自然对数e为底的对数,获得对应的值logZi、logXi及logYi
(2)以叶绿素a的浓度对数值为因变量,以总磷浓度对数值为自变量做非参数回归,并初步确定总磷突变点;
(3)以叶绿素a的浓度对数值为因变量,以总磷浓度对数值为自变量做分段线性回归,用步骤(2)中确定的总磷浓度突变点作为求解初始值,精确确定总磷突变点;
(4)用总氮浓度代替总磷浓度,重复步骤(2)和(3)精确确定总氮的突变点;
(5)根据总氮突变点和总磷突变点及湖泊所在区域的社会经济发展情况和国家水质标准,确定浅水湖泊候选营养盐水质管理目标;
(6)确定浅水湖泊水体用途所要求水质并与候选营养盐水质管理目标比较,最终确定具体湖泊营养盐水质管理目标。
2.根据权利要求1所述分段回归突变点确定浅水湖泊水质管理目标方法,其特征在于:所述步骤(2)具体包括以下步骤:
a.根据湖泊总磷浓度和叶绿素a的浓度的观测值数据数量n,用h=n-1/5计算最佳带宽值h;
b.在总磷浓度对数值的取值范围内均匀取1000个数值作为自变量logM的值logMp,即logM1、logM2、logM3...logM1000
c.从logM1开始计算得到所有满足条件logXi≤logM1+3h及logXi≥logM1-3h的总磷对数值logXi,对每一个满足条件的logXi及logM1计算相应的核函数K(logXi-logM1),作为logXi在点logM1上的权重,不满足条件logXi≤logM1+3h及logXi≥logM1-3h的logXi权重为0;其中:
K ( log X i - log M 1 ) = 1 2 π h e - 1 2 h 2 ( logX i - logM 1 ) 2
d.利用已有数值求解法,对logM1依次求解两组实数a(logM1),b(logM1);具体即要求a(logM1),b(logM1)使下式最小,
Σ i = 1 n { logZ i - a ( log M 1 ) - b ( log M 1 ) log X i } 2 K ( log X i - log M 1 )
e.重复步骤c和d,依次将logM2、logM3...logM1000替换logM1计算出两组实数a(logMp),b(logMp);
f.画出回归方程logZ=a(logMp)+b(logMp)logM的曲线,其中logZ代表回归方程的因变量,logM是回归方程的自变量;再次注意a(logMp),b(logMp)是回归方程的系数,是随着自变量logM的取值logMp的不同而不同的两组实数,观察方程logZ=a(logMp)+b(logMp)logM的曲线中的转折点,并记录该变量的突变点logγ对应的值作为分段回归求解迭代的初始值,如有多个总磷突变点,取最小的总磷浓度突变点。
3.根据权利要求2所述分段回归突变点确定浅水湖泊水质管理目标方法,其特征在于:所述步骤(3)的求解方法是求解四个实数β0,β1,δ和使下式最小,
其中是示性函数,也即满足下式:
β0,β1,δ的初始值为0,的初始值为logγ。
4.根据权利要求1所述分段回归突变点确定浅水湖泊水质管理目标方法,其特征在于:所述步骤(5)具体包括以下步骤:
a.根据湖泊所在区域的经济发展指标,确定该区域经济发展等级,我国各地区经济发展等级包括:
A.发达地区:区域人均国民生产总值超过全国人均国民生产总值的2倍,人民生活富足;
B.较发达地区:区域人均国民生产总值是全国人均国民生产总值的1.3至2倍,人民生活比较富足;
C.中等地区:区域人均国民生产总值是全国人均国民生产总值的0.8至1.3倍,人民生活水平中等;
D.落后地区:区域人均国民生产总值低于全国人均国民生产总值0.8倍,人民生活较贫困。
b.将总氮、总磷和叶绿素a的突变点值与国家《地表水环境质量标准(GB3838-2002)》比较,确定突变点所在水质等级区间;
c.发达地区应制定更严格的浅水湖泊营养盐管理目标,升高一个等级作为营养盐水质管理目标的候选;较发达地区取突变点所在水质等级区间对应的高等级作为营养盐水质管理目标的候选;中等地区取突变点所在水质等级区间对应的低等级作为营养盐水质管理目标的候选;落后地区在突变点基础上,降低一个等级作为营养盐水质管理目标的候选。
5.根据权利要求1所述分段回归突变点确定浅水湖泊水质管理目标方法,其特征在于:所述步骤(6)具体包括以下步骤:
a.确定湖泊水体用途及所要求水质:按照国家水质标准水体用途可分为自然保护区、饮用水源地、渔业、工业和农业五大类;依次对应于《地表水环境质量标准(GB3838-2002)》中的I、II、III、IV和V类水;
b.将候选营养盐水质管理目标与该湖泊水体用途所要求水质进行比较,最终取两者中水质较为严格的结果作为具体湖泊营养盐水质管理的目标。
CN201610615551.8A 2016-07-29 2016-07-29 分段回归突变点确定浅水湖泊水质管理目标方法 Expired - Fee Related CN106295135B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610615551.8A CN106295135B (zh) 2016-07-29 2016-07-29 分段回归突变点确定浅水湖泊水质管理目标方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610615551.8A CN106295135B (zh) 2016-07-29 2016-07-29 分段回归突变点确定浅水湖泊水质管理目标方法

Publications (2)

Publication Number Publication Date
CN106295135A true CN106295135A (zh) 2017-01-04
CN106295135B CN106295135B (zh) 2017-08-15

Family

ID=57663546

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610615551.8A Expired - Fee Related CN106295135B (zh) 2016-07-29 2016-07-29 分段回归突变点确定浅水湖泊水质管理目标方法

Country Status (1)

Country Link
CN (1) CN106295135B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108492042A (zh) * 2018-03-29 2018-09-04 河海大学 确定湖泊水质管理目标可接受范围及其可达性风险方法
CN109346134A (zh) * 2018-11-09 2019-02-15 扬州大学 一种黑臭水体关键控制变量及其目标值的确定方法
CN109460919A (zh) * 2018-11-09 2019-03-12 扬州大学 一种高精度确定湖泊营养物参照状态浓度的方法
CN117347327A (zh) * 2022-11-22 2024-01-05 中国环境科学研究院 利用叶绿素a和沉水植物评价浅水湖泊稳态的系统及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006101527A1 (en) * 2005-03-21 2006-09-28 Sensis Corporation Methods and apparatuses for evaluating water samples
CN102053140A (zh) * 2009-10-27 2011-05-11 华中科技大学 一种富营养化浅水湖泊生态修复后水质的诊断和评价方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006101527A1 (en) * 2005-03-21 2006-09-28 Sensis Corporation Methods and apparatuses for evaluating water samples
CN102053140A (zh) * 2009-10-27 2011-05-11 华中科技大学 一种富营养化浅水湖泊生态修复后水质的诊断和评价方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
汪靓等: "《中国东部浅水湖泊参照状态确定的非参数方法》", 《水科学进展》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108492042A (zh) * 2018-03-29 2018-09-04 河海大学 确定湖泊水质管理目标可接受范围及其可达性风险方法
CN108492042B (zh) * 2018-03-29 2022-04-08 河海大学 确定湖泊水质管理目标可接受范围及其可达性风险方法
CN109346134A (zh) * 2018-11-09 2019-02-15 扬州大学 一种黑臭水体关键控制变量及其目标值的确定方法
CN109460919A (zh) * 2018-11-09 2019-03-12 扬州大学 一种高精度确定湖泊营养物参照状态浓度的方法
CN109346134B (zh) * 2018-11-09 2022-05-10 扬州大学 一种黑臭水体关键控制变量及其目标值的确定方法
CN117347327A (zh) * 2022-11-22 2024-01-05 中国环境科学研究院 利用叶绿素a和沉水植物评价浅水湖泊稳态的系统及方法

Also Published As

Publication number Publication date
CN106295135B (zh) 2017-08-15

Similar Documents

Publication Publication Date Title
Luo et al. Evaluating the coordinated development of social economy, water, and ecology in a heavily disturbed basin based on the distributed hydrology model and the harmony theory
Kazbekov et al. Evaluating planning and delivery performance of water user associations (WUAs) in Osh Province, Kyrgyzstan
CN106295135B (zh) 分段回归突变点确定浅水湖泊水质管理目标方法
Longyang Assessing the effects of climate change on water quality of plateau deep-water lake-A study case of Hongfeng Lake
He et al. Reduction pathways identification of agricultural water pollution in Hubei Province, China
Xie et al. Exploring the dynamic correlation of landscape composition and habitat fragmentation with surface water quality in the Shenzhen river and deep bay cross-border watershed, China
Zhang et al. Effects of ecological protection and restoration on phytoplankton diversity in impounded lakes along the eastern route of China's South-to-North Water Diversion Project
CN112651548B (zh) 一种高原湖滨生态景观修复规划的评估识别方法
Yao et al. Evaluation of arable land suitability based on the suitability function-A case study of the Qinghai-Tibet Plateau
Ghezzo et al. Modeling the inter-annual variability of salinity in the lagoon of Venice in relation to the water framework directive typologies
CN115953061B (zh) 一种水资源短缺度估算方法、装置及计算机可读存储介质
Xu et al. Quantifying and optimizing agroecosystem services in China's Taihu Lake Basin
Lermontov et al. A fuzzy water quality index for watershed quality analysis and management
CN108629519A (zh) 服务于主体功能区的农业发展水平综合评估方法
Gikas Water quantity and hydrochemical quality monitoring of Laspias River, North Greece
Reder et al. European scenario studies on future in-stream nutrient concentrations
Zuo et al. A new method for water quality assessment: By harmony degree equation
CN115713448A (zh) 基于汇水单元的流域水污染常规因子快速溯源方法
Liao et al. A novel irrigation canal scheduling model adaptable to the spatial-temporal variability of water conveyance loss
Li et al. Assessment of water pollution in the Tibetan Plateau with contributions from agricultural and economic sectors: a case study of Lhasa River Basin
CN116562698A (zh) 蓄滞洪区协调发展评估方法、系统、设备和存储介质
Fan et al. The impact of clean water on infant mortality: evidence from China
CN116797095A (zh) 流域水生态安全评估方法、电子设备及存储介质
Kroiss et al. daNUbs: lessons learned for nutrient management in the Danube Basin and its relation to Black Sea euthrophication
CN107169874A (zh) 一种山区生态养殖的综合效益评价指标系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20170815

Termination date: 20200729