CN102135531A - 大型浅水湖泊72小时蓝藻水华预测方法 - Google Patents
大型浅水湖泊72小时蓝藻水华预测方法 Download PDFInfo
- Publication number
- CN102135531A CN102135531A CN2010106042071A CN201010604207A CN102135531A CN 102135531 A CN102135531 A CN 102135531A CN 2010106042071 A CN2010106042071 A CN 2010106042071A CN 201010604207 A CN201010604207 A CN 201010604207A CN 102135531 A CN102135531 A CN 102135531A
- Authority
- CN
- China
- Prior art keywords
- blue
- chlorophyll
- bloom
- concentration
- algae
- 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.)
- Pending
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种大型浅水湖泊72小时蓝藻水华预测方法,根据蓝藻原位生长率和漂移速率参数,构建水体叶绿素a含量预测模型和水华发生概率预报模型;分析监测水域即时的叶绿素a含量和分布并采集监测期内气象信息数据,输入构建的模型后,输出监测水域未来72小时叶绿素a浓度的预测值和蓝藻水华出现的水域及概率。本发明方法操作性强、用于预测的参数易获取,对蓝藻水华的暴发能做出迅速而准确的预判,时效性强。可用于水源地或重点景观水域蓝藻水华的预报。
Description
技术领域
本发明涉及一种湖泊、水库等水域的短时间尺度内蓝藻水华的监测预报方法,特别是涉及一种大型浅水湖泊72小时蓝藻水华预测方法。
背景技术
蓝藻水华是目前全球许多国家所面临的环境问题,在很多富营养的湖泊或水库中,蓝藻在局部湖区或库区堆积,并在高温下分解,形成恶臭;特别是当蓝藻在水源地取水口附近大量集聚时,就有可能引起水源地的水质恶化,危及供水安全。
目前还缺少有效的治理手段彻底地消除这一环境问题,蓝藻水华仍会大量出现于许多湖泊和水库中。因此,在研究蓝藻水华治理方法的同时,发展短期的蓝藻水华预测、预报技术,提前预知敏感湖区、尤其是水源地和重点景观湖区水华的发生几率,可提高环境管理部门的决策能力,并有利于相关部门及时采取应急措施应对水华污染。
目前国内外还没有关于短期蓝藻水华预报方法的报道。有关蓝藻水华监测预报的报道主要是采用逐步回归的统计学方法对浮游植物在长时间尺度上进行预测,如连续观测水温,透明度,pH,总氮(TN), 总磷(TP),溶解无机磷(DIP)等指标,经过连续观察,建立叶绿素与环境因子的关系,并在此基础上利用监测的数据来预测藻类密度。蓝藻水华是一定密度的藻类在特定气象、水文条件下上浮并大量聚集于水面的现象,因此,即时的藻类密度、气象、水文条件与水华的发生密切相关,如果能阐明其中的相互关系,并结合大量的野外调查数据,并根据相关参数建立模型,则能够据以对蓝藻水华进行预报。目前已有水体细胞色素快速监测的仪器,但是这些仪器在应用于藻类密度很大的水体时往往会有很大的误差,用于进行蓝藻监测时给预测结果带来较大的误差。卫星遥感图片可以从整体上识别水体中的蓝藻水华分布情况,可用于蓝藻水华预报工作准确性的评价。
发明内容
本发明提供了一种大型浅水湖泊72小时蓝藻水华的预测方法,基于水体中即时的蓝藻状况和未来的气象信息,利用构建的预测模型对蓝藻水华的发生进行短时间尺度的预报。
本发明方法的技术内容如下:
一种大型浅水湖泊72小时蓝藻水华预测方法,其特征在于:根据蓝藻原位生长率和漂移速率参数,构建水体叶绿素a含量预测模型和水华发生概率预报模型;分析监测水域即时的叶绿素a含量和分布状况;将监测水域即时的叶绿素a含量、分布状况,以及监测期内气象信息数据输入模型,输出监测水域预测期内叶绿素a浓度的预测值和蓝藻水华出现的水域及其概率,发布水华监测预报。
本发明方法包括以下步骤:
1)在监测水域设置采样点,采集分析监测水域即时的水质、蓝藻生物量和分布状况;
可根据监测水域的特点和水域面积,结合历史数据设置采样点,通过人工野外巡测,或使用离线或在线监测仪采集并分析获得监测水域即时叶绿素a浓度值和分布状况,以获得本发明方法进行水华预测所需的监测水域即时水质、蓝藻生物量和分布状况数据。
预报期间每周至少采样两次,采样后可利用快速水质监测仪器测定采样点各层水体的平均叶绿素浓度。
快速水质监测仪器在测量野外群体蓝藻叶绿素a含量时存在一定的误差,本发明方法建立了水体平均叶绿素浓度与快速水质监测仪器测定结果的校正关系式,即:
水体平均叶绿素a浓度=Chla YSI×7.76 (Chl YSI≤40μg/L)
水体平均叶绿素a浓度=Chla YSI×2.2 (Chl YSI>40μg/L)
式中Chla YSI为快速水质监测仪器测定结果。
2)采集预测期内的包括风速、风向、降雨在内的气象信息数据
根据气象部门的预报,采集监测水域预测期内未来天气情况(晴、阴以及降雨)以及风速、风向、温度等气象信息。
本发明选取并采集特定的监测数据,并据以建立短时间尺度的水华预报模型。本发明方法利用生态学原理,根据不同气候条件下蓝藻原位生长率,以及不同风向和风速条件下水华蓝藻的漂移速率参数,构建水体叶绿素a含量预报模型和水华发生概率预报模型。
蓝藻水华发生与很多因素相关,如藻类密度,风速,风向,降雨等,建立蓝藻水华的预报模型通常用如下公式进行表达:
F (水华概率)= f (藻类密度,风速,风向,降雨,有效光辐射,湖流,总氮,总磷,
水温,电导率,溶解无机氮,溶解无机磷,pH,溶氧,…) 公式1
为本发明目的,根据历史数据可将营养盐等在长时间尺度上影响蓝藻生物量的环境因素忽略,将公式1简化为如下公式:
F=f(N t Vw, R) 公式2
式中,F为水华的发生频率,Nt为对应时刻水体中的藻类密度,Vw为该时刻的风速,R为降雨情况。风速和降雨情况可以通过气象预报获得,而水体中藻类的数量则可根据生态学原理由如下公式计算:
N t = N0+(Bt-Dt)+(It-Et) 公式3
式中,Nt为所预报时刻(t)的水体中藻类密度,N0为监测时刻(0)藻类密度,Bt、Dt、It、Et分别是0—t时段该湖区藻类的生长、死亡、迁入、迁出的数量。
本发明根据上述生态学原理,分别建立用于短时间内蓝藻水华的水体叶绿素a浓度预测模型和水华发生概率预报模型:
3)构建水体叶绿素a浓度预测模型,输入监测水域即时的叶绿素a浓度值、分布状况以及气象信息数据,输出监测水域预测期内叶绿素a浓度的预测值;
为定量描述水域叶绿素a浓度的分布及变化过程,将监测水体网格化,水体网格单元大小为am×bm的矩形,通常a、b可根据预测要求在200---500m范围内选择。
如上文所述,水体中藻类的数量则可根据公式3计算:
N t = N0+(Bt-Dt)+(It-Et)
蓝藻由于生长、死亡引起的变化量,可用以下公式计算:
式中,q为叶绿素a浓度的日变化率,各气象条件下藻类日变化率用表1的参数来进行计算,N0可以通过步骤1)现场监测获得。
对于水体网格单元内迁入(或迁出)的蓝藻量,风速是决定蓝藻迁移的主要因素,迁入(或迁出)的蓝藻量可用以下公式计算:
E=a×V×t×d×P×C 公式5
式中,a为网格的边长,E为迁入(或迁出)的藻类量,V为风速,t为时间(1h,3600s),d为水深,P为发生漂移的藻类占总藻量的百分比,C为水体中叶绿素a的浓度,相关参数根据观测数据,按下表2选取。以此计算出所有漂移出(移入)的藻类。
表1 不同气象条件下蓝藻由于原位生长死亡引起的日变化率(q)
天气、风速 | 1-2级 | 2-3级 | 3-4级 | 4-5级 | 5级以上 |
晴 | 0.6 | 0.4 | 0.3 | 0.08 | 0 |
晴到多云 | 0.5 | 0.3 | 0.2 | 0.06 | 0 |
多云 | 0.4 | 0.2 | 0.07 | 0.05 | 0 |
多云到阴 | 0.2 | 0.1 | 0.05 | 0.05 | 0 |
阴,降雨 | 0 | 0 | 0 | 0 | 0 |
表2 不同风力条件下漂移的藻类比例和对应的漂移速度
风力(风速m/s) | 漂移藻类占总藻的百分比(%) | 漂移速度(m/s) |
1-2级(0.3-3.3) | 50 | 0.07 |
3级(3.4-5.4) | 40 | 0.04 |
4级(5.5-7.9) | 25 | 0.03 |
5级(8.0-10.7) | 10 | 0.01 |
5级以上(>10.8) | 5 | 0.005 |
4)构建水华概率预测模型,输入监测水域叶绿素a浓度值、分布状况的预测值以及气象信息数据,输出蓝藻水华出现的水域及其概率;
为实现本发明目的,本发明方法将叶绿素a浓度,降雨和风速作为水华形成的共同的触发因子,预测模型各因子之间的关系如下:
F=f1 (Nt)×f2(V)×f3(R) 公式6
其中f1 (Nt)为t时刻由藻类密度引起的水华发生概率;f2(V)为由风速条件引起的水华发生概率;f3(R)为降雨条件引起的水华发生概率。
根据实际监测数据,将叶绿素a含量划分为5个等级,并将相应气象条件下的水华发生概率赋值(表3)。
表3蓝藻水华预报中叶绿素a含量、风速和降雨与水华发生概率的对应表
叶绿素a(μg/L) | f1(Nt) | 风力(风速m/s) | f2(V) | 降雨 | f3(R) |
60以上 | 1 | 1-2级(0.3-3.3) | 1 | 晴、多云 | 1 |
50 | 0.9 | 3级(3.4-5.4) | 0.9 | 阴、小雨 | 0.9 |
40 | 0.8 | 4级(5.5-7.9) | 0.8 | 阵雨、雷阵雨 | 0.8 |
30 | 0.6 | 5级(8.0-10.7) | 0.7 | 中雨 | 0.7 |
20 | 0.4 | 5级以上(>10.8) | 0.5 | 中、大、暴雨 | 0 |
将监测水域蓝藻生物量、分布状况以及气象信息数据输入预报模型,得到预测期内蓝藻水华在特定水域出现的概率及其强度。
5)水华预报发布
根据模型输出数据,发布湖区未来三天的叶绿素a浓度等值线图和未来水华趋势的预报。水华的预报以文字阐述并结合图表对未来三天特定水域蓝藻水华的概率和程度进行描述,同时提供相关的建议。
本发明方法还可以利用遥感影像数据对水华预报进行修正,并对本方法的预报方法进行评价。
附图说明
图1为本发明的72小时蓝藻水华预测方法的工作流程图(T0表示预测开始的时间,T0+4h表示四小时以后,余类推);
图2A-2C 为本发明方法预测的太湖水域未来72小时叶绿素a浓度的等值线图;
图3A-3C 为本发明方法预测的太湖水域未来72小时叶绿素a的浓度分布图。
具体实施方式
按照本发明方法在2007,2008,2009年水华发生期间对太湖蓝藻水华进行预报。工作流程图参见图1。
72小时蓝藻水华预测方法,包括以下步骤:
1)资料的采集和应用
包括分析太湖水域即时的水质、蓝藻生物量和分布状况和气象信息数据的采集。
人工巡测:预报期间每周采样两次,根据历史数据,北太湖的三个湖湾(梅梁湾、贡湖湾和竺山湾)是蓝藻水华发生比较频繁的区域,而梅梁湾、贡湖湾又是水源地和景观水域,因此,在太湖区域对蓝藻水华的预报主要关注梅梁湾、贡湖湾的水华状况。在这两个湖区的每个水厂取水口附近和景观水域设立采样点,采样分析以获得关注区域的叶绿素a浓度,其它水域则选择有代表性的位点设立采样点,以获得该水域整体的叶绿素a浓度。
到达采样点后,用2m的PVC管采集整水柱,置于测量桶中摇匀,将多功能水质参数仪(YSI6600)放入测量桶中进行测定,以获得采样点各层水体的平均叶绿素浓度。
数据的校正:由于YSI6600在测量野外群体蓝藻叶绿素a含量时存在一定的误差,根据YSI6600测定结果与室内化学法测定的比较分析,按下式根据YSI6600测定结果对水体平均叶绿素浓度a进行校正:
水体平均叶绿素a浓度=Chla YSI×7.76 (Chl YSI≤40μg/L)
水体平均叶绿素a浓度=Chla YSI×2.2 (Chl YSI>40μg/L)
式中Chla YSI为YSI6600测定结果。
遥感数据:除非湖面上有云层覆盖,每天10点左右自太湖上空获得MODIS遥感影像。
气象数据:来自气象部门的关于太湖区域的气象预报,应用的主要参数有天气情况(晴、阴以及降雨)、风速、风向、温度。
2007年 7 月 26 日采集的数据如下表4所示:
表4 2007年7月26日太湖各巡测点叶绿素a含量以及未来72h气象信息
人工巡测数据和气象数据用于输入模型进行叶绿素a含量和水华概率预测,遥感数据能获得湖面蓝藻水华大致的分布状况,可用于人工巡测数据的校正和预报结果的评价。
2)水华预报
1. 水体中叶绿素含量的预测
为定量描述太湖叶绿素a浓度的变化过程,本发明将太湖水体网格化,水体网格单元大小为500m×500m。建立生态学模型,以人工巡测获得初始时刻的藻类密度(N0)和采集的气象数据为输入量,分别预测蓝藻由于原位生长、死亡的变化量,和蓝藻由于迁入、迁出的变化量,从而预测每个水体单元格中在预测期内的不同时刻叶绿素a浓度。
对于所构建的模型,采用特定的计算机语言(如python)进行编程,确立太湖的边界条件,并根据划分形成的特定大小的网格,在windows条件下输入起始条件(叶绿素a浓度、气象预报等数据),即可进行运算,得到每个网格内t时段后叶绿素a的浓度,进而得到各个湖区的藻类密度情况。
根据已有的研究结果(Tsjimura,2000;吴晓东等,2008),蓝藻生长速率多数在0.2-0.3,裂解速率在0.02-0.2之间,藻类的生长和死亡速率均与气象条件相关。用Kriging空间插值法以500m×500m的大小获得各湖区网格单元,具体网格内的藻类总量按该网格内插值所得到的叶绿素a值计算总量,以如下公式计算t时刻后每个网格内藻类原位生长、死亡后的总量:
式中,q为叶绿素a浓度的日变化率,各气象条件下藻类日变化率用表1的参数来进行计算。该公式中原位生长的藻类以原有藻类的基数为准,不考虑漂移的变化量。
风速是决定蓝藻上浮的主要因素,根据实际情况和相关研究结果,将不同的风力条件下蓝藻的漂移量和漂移速度按如表2计算。
藻类漂移与风速风向相关。由于风是矢量,结合相关部门提供的气象预报信息,将东、东南、南、西南、西、西北、北和东北的8个方位角依次赋值为1、2、3、4、5、6、7、8,并以此确定不同风向条件下模型的计算网格,其中2、4、6、8风向条件下将网格的漂移方向按水平面的两个方向计算,依据漂移的距离来确定漂移的藻类进入相邻网格的量。
例如,在风速为1(东)方向的条件下,确定靠近东岸的网格为起始运算网格,其每个网格发生漂移而出的藻类量为:
E1=500×V1×t1×d1×P1×C1 公式6
式中,500为网格的边长,E1为漂移出的藻类量,V1为风速,t1为时间(1h,3600s),d1为水深,P1为该网格内发生漂移的藻类占总藻量百分比,C1为该网格内水体中叶绿素a的浓度。以此计算出所有网格内漂移出的藻类,至于最西岸的网格,由于其没有迁出,将这些网格内的E1设为0。
同理,该网格内由上一网格迁移而入的藻类量,即为上一网格漂移出的藻类量,计算公式为:
I1=E1’=500×V1×t1×d1×P1×C1’ 公式7
式中各参数同以上公式,C1’为东边网格中的藻类密度。
同理,计算出风向为3,5,7情况下各个网格内藻类漂移迁出和迁入情况。
在风向为2(东南)的条件下,将藻类漂移速度分解为向东和向南两个方向的矢量,则每个方向的漂移速度为0.707V。运算开始时,以最靠近东岸和南岸的所有网格为计算的初始网格,不过藻类的漂移总是在与其相邻的三个网格内发生(即北部、西部和西北部的网格),由于风速相同,风向视为理想条件下的45度,则漂移进入该三个网格内藻类总量相等:
EN=EW=ENW= 0.707 V1×0.707 V1××t1×d1×P1×C1 公式8
式中,EN, EW, ENW分别为漂移进入北部、西部和西北部网格内的藻类总量,其他参数同公式,总漂移迁出量为(EN+EW+ENW)。
漂移进该网格的数量则由东边、南边和东南边的三个网格内藻类而来,记为IE, IS和ISE,其计算公式同公式,只是藻类密度分别为对应网格内的藻类密度,即:
I=IE+ IS+ I
计算时,将最靠近东边和南边的网格中将漂移进入的藻类设为0;其中最靠近北边的网格向北漂移的量设为0,最靠近西边的网格内向西的漂移量设为0。
同理,计算出4,6,8风向条件下各个网格内的漂移藻类密度;
根据以上公式,分别计算出(B t -D t )、I t 和E t ,即可以计算出不同气象条件下不同时间内的各个网格内的藻类密度。
2. 水华概率预测
公式2确定了主要的模型参数,预报中将降雨和风速视为水华形成的共同的触发因子,从而预测模型各因子之间的关系如下:
F=f1 (Nt)×f2(V)×f3(R) 公式6
其中f1 (Nt)为t时刻由藻类密度引起的水华发生概率;f2(V)为由风速条件引起的水华发生概率;f3(R)为降雨条件引起的水华发生概率。
蓝藻水华发生的基础是叶绿素a含量达到某一阈值。通常情况下,形成水华的水体的叶绿素a浓度高于10μg/L,但由于不同的水体中叶绿素浓度和水华的关系之间存在差异,而且不同的水体人们的关注程度不同,因此对水华的准确界定也存在差异,目前还没有统一的标准来确定水华发生时绿素浓度的阈值。太湖水体中夏季叶绿素含量高,根据长期太湖的野外观测,发现当水体叶绿素a含量达20μg/L时,即有可能发生水华,且水华概率随着叶绿素a含量的增加而增加,因此,将叶绿素a含量划分为5个等级,并将相应的水华发生概率赋值(表3)。
降雨情况本身是一个二态数据(即降雨有降雨或没有降雨),仅在雨量方面有大小之分。可是当前常规的气象预报都无法获得对应时刻的Vw(风速)和R(降雨),因此也需要对天气预报中的参数进行赋值。在预报中,根据太湖实际监测的情况,将f2(V)和f3(R)分别赋值如表3所示。
3) 水华预报的内容和形式
水华的预报包括各湖区监测点即时叶绿素a含量、未来3日气象数据、未来三天特定水域蓝藻水华的概率和程度预测、叶绿素a的浓度等值线图、水华现状及未来水华趋势的文字描述。叶绿素a含量可用surfer软件制等高线图输出,如将上文表4数据输入模型后,可输出图2A-2C所示的2007年7 月26日预测的未来72小时叶绿素a的等高线图,并做出7月27日至29日三天内水华预报如下:
梅梁湾蓝藻水华发生概率90%,其中马山水厂周围比较严重;贡湖湾蓝藻水华发生概率90%,主要区域在北部沿岸,特别是南泉水厂附近。
又如,在太湖梅梁湾、贡湖湾2009年4月15日监测或采集如表5所示的数据,根据本发明方法得到如图3A-3C所示的叶绿素a浓度分布图,并预测4月16日至 4月18日三天内梅梁湾水面(直湖港、湾口西、湾东北)灾害性蓝藻水华发生概率30%,贡湖湾水面(贡湖与大太湖的交界处)灾害性蓝藻水华发生概率为10%,其中叶绿素a含量在Arcgis条件支持下编程后直接输出。
表5 2009年4月15日太湖各巡测点叶绿素a含量以及未来72h气象信息
4) 水华预报的评价
采用了MODIS遥感影像作为参照,得到实际的蓝藻水华发生情况,根据不同的湖区分别对水华预报概率准确率和地点准确率进行评价。
当预报的概率>50%(或<50%)时,对应湖区出现(或未出现)蓝藻水华,评价预报为正确,反之为错误;小面积湖区(湖湾)蓝藻预报有水华湖区和实际发生水华的湖区在70%以上水域相符时,则预报地点为正确;水域面积大的湖区,将每次预报的水华水域与实际水华水域面积对比后加以评价。
采用本发明方法在2009年发布了44份蓝藻水华预报,每次预报三天,共预报132天。梅梁湾内水华概率预报正确率为70%;共30天发生水华,水华发生地点预报正确26天,在概率预报正确的情况下,水华地点的正确率为64%。贡湖湾内水华概率预报正确率为87%;共15天发生水华,在概率预报正确前提下,水华发生地点预报正确10天,正确率为79%。
预报结果评价详见表6,其中预报概率的正确性在80%左右,水华地点的预报正确性也在50%以上。
年太湖蓝藻水华预报的回顾性评价
湖区 | 评价次数(天) | 概率预报正确率 | 地点预报正确率 |
梅梁湾 | 63 | 70% | 65% |
贡湖 | 63 | 87% | 79% |
其他湖区 | 63 | 95% | 81%(基本吻合为56%,预测面积为实际发生面积一半水域为29%) |
从更长时间来看,2007年、2008年和2009年分别发布了32、54和44份蓝藻水华预报。其中,蓝藻水华发生概率预测的准确性在大太湖为88%-95%,在贡湖和梅梁湾为60%-87%之间;对蓝藻水华发生的具体地点的预测,贡湖湾和梅梁湾准确性介于60%-84%。
Claims (9)
1. 一种大型浅水湖泊72小时蓝藻水华预测方法,其特征在于:根据蓝藻原位生长率和漂移速率参数,构建水体叶绿素a含量预测模型和水华发生概率预报模型;分析监测水域即时的叶绿素a含量和分布;将监测水域即时的叶绿素a含量、分布,以及监测期内气象信息数据输入模型,输出监测水域预测期内叶绿素a浓度的预测值和蓝藻水华出现的水域及概率,发布水华监测预报。
2. 根据权利要求1所述的蓝藻水华预测方法,其特征在于所述方法包括以下步骤:
1)在监测水域设置采样点,采集分析监测水域即时叶绿素a浓度值和分布;
2)采集预测期内的包括风速、风向、降雨在内的气象信息数据;
3)构建水体叶绿素a浓度预测模型,输入监测水域即时的叶绿素a浓度值、分布以及气象信息数据,输出监测水域预测期内叶绿素a浓度的预测值;
4)构建水华发生概率预测模型,输入监测水域叶绿素a浓度值、分布的预测值以及气象信息数据,预测蓝藻水华出现的水域及其概率;
5)输出数据,发布水华监测预报。
3. 根据权利要求1或2所述的蓝藻水华预测方法,其特征在于所述的水体中叶绿素a浓度采用快速水质监测仪器测定并根据下列公式校正:
水体平均叶绿素a浓度=Chla YSI×7.76 (Chl YSI≤40μg/L)
水体平均叶绿素a浓度=Chla YSI×2.2 (Chl YSI>40μg/L)
式中Chla YSI为快速水质监测仪器测定结果。
4. 根据权利要求1或2所述的蓝藻水华预测方法,其特征在于所述的叶绿素a浓度预测模型将监测水域划分为网格单元,并按下式得到水体网格单元内叶绿素a浓度的预测值:
N t = N0+(Bt-Dt)+(It-Et)
式中,Nt为所预报时刻(t)的水体中藻类密度,N0为监测时刻(0)藻类密度,Bt、Dt、It、Et分别是0—t时段该湖区藻类的生长、死亡、迁入、迁出的数量。
6. 根据权利要求4所述的蓝藻水华预测方法,其特征在于水体网格单元内迁入或迁出的蓝藻量用以下公式计算:
Et=a×V×t×d×P×C
式中,a为网格的边长,Et为迁入或迁出的藻类量,V为风速,t为时间,d为水深,P为发生漂移的藻类占总藻量百分比,C为水体中叶绿素a的浓度。
7.根据权利要求1或2所述的蓝藻水华预测方法,其特征在于所述的水华发生概率预测模型,按下式计算水华发生概率:
F=f1 (Nt)×f2(V)×f3(R)
其中f1 (Nt)为t时刻由藻类数量引起的水华发生概率;f2(V)为由风速条件引起的水华发生概率;f3(R)为降雨条件引起的水华发生概率;根据输入的叶绿素a和气象数据,按下表1对水华发生概率赋值;
表1蓝藻水华预报中叶绿素a含量、风速和降雨与水华发生概率的对应表
。
8. 根据权利要求1或2所述的蓝藻水华预测方法,其特征在于所述的输出数据包括监测水域叶绿素a浓度预测值的等值线图或浓度分布图、特定水域蓝藻水华发生的概率。
9. 根据权利要求1或2所述的蓝藻水华预测方法,其特征在于利用遥感影像数据对水华预报进行修正,并对预报方法进行评价。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010106042071A CN102135531A (zh) | 2010-12-24 | 2010-12-24 | 大型浅水湖泊72小时蓝藻水华预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010106042071A CN102135531A (zh) | 2010-12-24 | 2010-12-24 | 大型浅水湖泊72小时蓝藻水华预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102135531A true CN102135531A (zh) | 2011-07-27 |
Family
ID=44295376
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010106042071A Pending CN102135531A (zh) | 2010-12-24 | 2010-12-24 | 大型浅水湖泊72小时蓝藻水华预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102135531A (zh) |
Cited By (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102539398A (zh) * | 2011-12-08 | 2012-07-04 | 上海市环境科学研究院 | 浮游藻类叶绿素a的现场实时测定方法 |
CN102735671A (zh) * | 2012-07-02 | 2012-10-17 | 中国科学院南京地理与湖泊研究所 | 大型浅水湖泊越冬蓝藻种源疏浚位点的确定方法 |
CN103063202A (zh) * | 2012-12-30 | 2013-04-24 | 同济大学 | 基于遥感影像的蓝藻生物量时空变化监测与可视化方法 |
CN103743700A (zh) * | 2014-01-17 | 2014-04-23 | 中国科学院南京地理与湖泊研究所 | 一种大型浅水湖泊蓝藻水华modis卫星高精度监测方法 |
CN103913422A (zh) * | 2014-03-03 | 2014-07-09 | 中国科学院南京地理与湖泊研究所 | 一种基于hj-ccd影像的浅水湖泊水生植物快速监测方法 |
CN104091072A (zh) * | 2014-07-11 | 2014-10-08 | 中国科学院南京地理与湖泊研究所 | 一种保障水源地安全的蓝藻密度阈值的建立方法 |
CN104569306A (zh) * | 2014-12-24 | 2015-04-29 | 中国科学院重庆绿色智能技术研究院 | 一种预判藻类种群演替和水华的方法 |
CN104615885A (zh) * | 2015-02-04 | 2015-05-13 | 中国科学院南京地理与湖泊研究所 | 一种富营养化浅水湖泊藻源性湖泛的短期预报方法 |
CN104699979A (zh) * | 2015-03-24 | 2015-06-10 | 北京工商大学 | 基于复杂网络的城市湖库藻类水华混沌时间序列预测方法 |
CN104899653A (zh) * | 2015-06-02 | 2015-09-09 | 北京工商大学 | 基于专家***及蓝藻生长机理时序模型的湖库蓝藻水华预测方法 |
CN105403680A (zh) * | 2015-12-02 | 2016-03-16 | 浙江省海洋监测预报中心 | 一种基于浮标在线监测数据的赤潮短期预警方法及用途 |
CN105550493A (zh) * | 2015-12-02 | 2016-05-04 | 江西理工大学 | 基于云模型和声搜索的水体叶绿素a浓度预测方法 |
CN105606075A (zh) * | 2015-12-18 | 2016-05-25 | 中国科学院南京地理与湖泊研究所 | 一种基于局部风速的大型浅水湖泊水华微囊藻垂向分布模式的判别方法 |
CN103886218B (zh) * | 2014-04-08 | 2016-08-24 | 北京工商大学 | 基于多元非平稳时间序列分析与神经网络及支持向量机补偿的湖库藻类水华预测方法 |
CN106990216A (zh) * | 2017-04-01 | 2017-07-28 | 河海大学 | 一种浅水湖泊水华风险分析预警***及其分析预警方法 |
CN107340365A (zh) * | 2017-06-19 | 2017-11-10 | 中国科学院南京地理与湖泊研究所 | 一种面向湖泊蓝藻灾害的立体监控及数据挖掘***和方法 |
WO2018098736A1 (zh) * | 2016-11-30 | 2018-06-07 | 深圳益强信息科技有限公司 | 一种基于遥感成像的蓝藻监测方法及装置 |
CN108226437A (zh) * | 2018-04-09 | 2018-06-29 | 山东建筑大学 | 一种应用于湖库水体的修正水质监测结果的方法 |
CN108416460A (zh) * | 2018-01-19 | 2018-08-17 | 北京工商大学 | 基于多因素时序—随机深度置信网络模型的蓝藻水华预测方法 |
CN108893519A (zh) * | 2018-07-24 | 2018-11-27 | 成都蔚瑞生物科技有限公司 | 一种基于藻春化太湖铜绿微囊藻暴发的藻密度预测方法 |
CN109029366A (zh) * | 2018-07-02 | 2018-12-18 | 中国科学院南京地理与湖泊研究所 | 一种通过飞行器拍摄测定藻斑水平漂移速率的方法 |
CN109598439A (zh) * | 2018-12-03 | 2019-04-09 | 河海大学 | 一种利用风速辅助浅水湖泊水环境改善的方法 |
CN109657200A (zh) * | 2018-12-05 | 2019-04-19 | 北京师范大学 | 一种确定湖泊水库蓝藻水华爆发几率的方法 |
CN109711066A (zh) * | 2018-12-28 | 2019-05-03 | 南开大学 | 一种浅水型小型湖库水华预测方法及预测模型 |
CN109959621A (zh) * | 2017-12-22 | 2019-07-02 | 北京航天泰坦科技股份有限公司 | 一种基于中高分辨率时序遥感数据的水环境分析方法 |
CN110836683A (zh) * | 2018-08-15 | 2020-02-25 | 光大水务科技发展(南京)有限公司 | 河道管理方法、河道管理装置及终端 |
CN112989281A (zh) * | 2021-02-20 | 2021-06-18 | 中国科学院南京地理与湖泊研究所 | 一种基于遥感藻总量的藻华预测方法 |
CN113158762A (zh) * | 2021-02-23 | 2021-07-23 | 中国科学院南京地理与湖泊研究所 | 大型湖库滨岸带蓝藻水华堆积风险评估方法 |
CN113711087A (zh) * | 2019-04-15 | 2021-11-26 | 国立研究开发法人理化学研究所 | 环境因子预测装置、方法、程序、学习完毕模型以及存储介质 |
CN114169646A (zh) * | 2021-12-31 | 2022-03-11 | 中科三清科技有限公司 | 水华预测方法、装置、电子设备及计算机可读存储介质 |
CN114324796A (zh) * | 2021-12-28 | 2022-04-12 | 中国科学院南京地理与湖泊研究所 | 一种水草上附着丝状藻类水华的预测方法 |
CN114485781A (zh) * | 2021-12-30 | 2022-05-13 | 常州市环境科学研究院 | 一种用于浅水湖泊的沉浮式蓝藻水华监测***及预报方法 |
CN114564883A (zh) * | 2022-02-11 | 2022-05-31 | 水利部交通运输部国家能源局南京水利科学研究院 | 融合机理与ML的湖泊叶绿素a短期集合预报方法及*** |
CN115508287A (zh) * | 2022-10-28 | 2022-12-23 | 生态环境部卫星环境应用中心 | 湖库蓝藻水华提取方法和装置 |
CN115719174A (zh) * | 2022-10-26 | 2023-02-28 | 生态环境部卫星环境应用中心 | 土地利用类型与蓝藻水华风险定量关系确定方法和装置 |
CN116485045A (zh) * | 2023-06-21 | 2023-07-25 | 中科三清科技有限公司 | 一种建立水华预测模型的方法、水华预测方法及装置 |
CN118212970A (zh) * | 2024-05-16 | 2024-06-18 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于时序和空间分布的蓝藻预测方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004294380A (ja) * | 2003-03-28 | 2004-10-21 | Japan Science & Technology Agency | 水中分光照度からのクロロフィルa鉛直分布の推定方法 |
CN101900687A (zh) * | 2010-07-06 | 2010-12-01 | 重庆大学 | 一种基于图像处理的小型水域水华监测预警方法 |
-
2010
- 2010-12-24 CN CN2010106042071A patent/CN102135531A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004294380A (ja) * | 2003-03-28 | 2004-10-21 | Japan Science & Technology Agency | 水中分光照度からのクロロフィルa鉛直分布の推定方法 |
CN101900687A (zh) * | 2010-07-06 | 2010-12-01 | 重庆大学 | 一种基于图像处理的小型水域水华监测预警方法 |
Non-Patent Citations (3)
Title |
---|
刘苑等: "YSI(多参数水质检测仪)测定叶绿素a浓度的准确性及误差探讨", 《湖泊科学》 * |
孔繁翔等: "太湖蓝藻水华的预防、预测和预警的理论与实践", 《湖泊科学》 * |
黄佳聪等: "《蓝藻水华预报模型及基于遗传算法的参数优化》", 《生态学报》 * |
Cited By (59)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102539398A (zh) * | 2011-12-08 | 2012-07-04 | 上海市环境科学研究院 | 浮游藻类叶绿素a的现场实时测定方法 |
CN102735671A (zh) * | 2012-07-02 | 2012-10-17 | 中国科学院南京地理与湖泊研究所 | 大型浅水湖泊越冬蓝藻种源疏浚位点的确定方法 |
CN102735671B (zh) * | 2012-07-02 | 2014-07-02 | 中国科学院南京地理与湖泊研究所 | 大型浅水湖泊越冬蓝藻种源疏浚位点的确定方法 |
CN103063202A (zh) * | 2012-12-30 | 2013-04-24 | 同济大学 | 基于遥感影像的蓝藻生物量时空变化监测与可视化方法 |
CN103063202B (zh) * | 2012-12-30 | 2015-04-15 | 同济大学 | 基于遥感影像的蓝藻生物量时空变化监测与可视化方法 |
CN103743700B (zh) * | 2014-01-17 | 2016-07-06 | 中国科学院南京地理与湖泊研究所 | 一种大型浅水湖泊蓝藻水华modis卫星高精度监测方法 |
CN103743700A (zh) * | 2014-01-17 | 2014-04-23 | 中国科学院南京地理与湖泊研究所 | 一种大型浅水湖泊蓝藻水华modis卫星高精度监测方法 |
CN103913422A (zh) * | 2014-03-03 | 2014-07-09 | 中国科学院南京地理与湖泊研究所 | 一种基于hj-ccd影像的浅水湖泊水生植物快速监测方法 |
CN103886218B (zh) * | 2014-04-08 | 2016-08-24 | 北京工商大学 | 基于多元非平稳时间序列分析与神经网络及支持向量机补偿的湖库藻类水华预测方法 |
CN104091072A (zh) * | 2014-07-11 | 2014-10-08 | 中国科学院南京地理与湖泊研究所 | 一种保障水源地安全的蓝藻密度阈值的建立方法 |
CN104091072B (zh) * | 2014-07-11 | 2017-05-24 | 中国科学院南京地理与湖泊研究所 | 一种保障水源地安全的蓝藻密度阈值的建立方法 |
CN104569306A (zh) * | 2014-12-24 | 2015-04-29 | 中国科学院重庆绿色智能技术研究院 | 一种预判藻类种群演替和水华的方法 |
CN104615885A (zh) * | 2015-02-04 | 2015-05-13 | 中国科学院南京地理与湖泊研究所 | 一种富营养化浅水湖泊藻源性湖泛的短期预报方法 |
CN104699979A (zh) * | 2015-03-24 | 2015-06-10 | 北京工商大学 | 基于复杂网络的城市湖库藻类水华混沌时间序列预测方法 |
CN104899653B (zh) * | 2015-06-02 | 2018-01-12 | 北京工商大学 | 基于专家***及蓝藻生长机理时序模型的湖库蓝藻水华预测方法 |
CN104899653A (zh) * | 2015-06-02 | 2015-09-09 | 北京工商大学 | 基于专家***及蓝藻生长机理时序模型的湖库蓝藻水华预测方法 |
CN105550493B (zh) * | 2015-12-02 | 2017-12-15 | 江西理工大学 | 基于云模型和声搜索的水体叶绿素a浓度预测方法 |
CN105403680B (zh) * | 2015-12-02 | 2017-09-08 | 浙江省海洋监测预报中心 | 一种基于浮标在线监测数据的赤潮短期预警方法及用途 |
CN105403680A (zh) * | 2015-12-02 | 2016-03-16 | 浙江省海洋监测预报中心 | 一种基于浮标在线监测数据的赤潮短期预警方法及用途 |
CN105550493A (zh) * | 2015-12-02 | 2016-05-04 | 江西理工大学 | 基于云模型和声搜索的水体叶绿素a浓度预测方法 |
CN105606075A (zh) * | 2015-12-18 | 2016-05-25 | 中国科学院南京地理与湖泊研究所 | 一种基于局部风速的大型浅水湖泊水华微囊藻垂向分布模式的判别方法 |
CN105606075B (zh) * | 2015-12-18 | 2017-11-28 | 中国科学院南京地理与湖泊研究所 | 一种基于局部风速的大型浅水湖泊水华微囊藻垂向分布模式的判别方法 |
WO2018098736A1 (zh) * | 2016-11-30 | 2018-06-07 | 深圳益强信息科技有限公司 | 一种基于遥感成像的蓝藻监测方法及装置 |
CN106990216B (zh) * | 2017-04-01 | 2019-03-29 | 河海大学 | 一种浅水湖泊水华风险分析预警***及其分析预警方法 |
CN106990216A (zh) * | 2017-04-01 | 2017-07-28 | 河海大学 | 一种浅水湖泊水华风险分析预警***及其分析预警方法 |
CN107340365A (zh) * | 2017-06-19 | 2017-11-10 | 中国科学院南京地理与湖泊研究所 | 一种面向湖泊蓝藻灾害的立体监控及数据挖掘***和方法 |
CN107340365B (zh) * | 2017-06-19 | 2019-04-26 | 中国科学院南京地理与湖泊研究所 | 一种面向湖泊蓝藻灾害的立体监控及数据挖掘***和方法 |
CN109959621B (zh) * | 2017-12-22 | 2021-07-20 | 北京航天泰坦科技股份有限公司 | 一种基于中高分辨率时序遥感数据的水环境分析方法 |
CN109959621A (zh) * | 2017-12-22 | 2019-07-02 | 北京航天泰坦科技股份有限公司 | 一种基于中高分辨率时序遥感数据的水环境分析方法 |
CN108416460A (zh) * | 2018-01-19 | 2018-08-17 | 北京工商大学 | 基于多因素时序—随机深度置信网络模型的蓝藻水华预测方法 |
CN108416460B (zh) * | 2018-01-19 | 2022-01-28 | 北京工商大学 | 基于多因素时序—随机深度置信网络模型的蓝藻水华预测方法 |
CN108226437A (zh) * | 2018-04-09 | 2018-06-29 | 山东建筑大学 | 一种应用于湖库水体的修正水质监测结果的方法 |
CN109029366A (zh) * | 2018-07-02 | 2018-12-18 | 中国科学院南京地理与湖泊研究所 | 一种通过飞行器拍摄测定藻斑水平漂移速率的方法 |
CN108893519A (zh) * | 2018-07-24 | 2018-11-27 | 成都蔚瑞生物科技有限公司 | 一种基于藻春化太湖铜绿微囊藻暴发的藻密度预测方法 |
CN110836683A (zh) * | 2018-08-15 | 2020-02-25 | 光大水务科技发展(南京)有限公司 | 河道管理方法、河道管理装置及终端 |
CN109598439A (zh) * | 2018-12-03 | 2019-04-09 | 河海大学 | 一种利用风速辅助浅水湖泊水环境改善的方法 |
CN109598439B (zh) * | 2018-12-03 | 2021-06-04 | 河海大学 | 一种利用风速辅助浅水湖泊水环境改善的方法 |
CN109657200A (zh) * | 2018-12-05 | 2019-04-19 | 北京师范大学 | 一种确定湖泊水库蓝藻水华爆发几率的方法 |
CN109657200B (zh) * | 2018-12-05 | 2020-10-23 | 北京师范大学 | 一种确定湖泊水库蓝藻水华爆发几率的方法 |
CN109711066A (zh) * | 2018-12-28 | 2019-05-03 | 南开大学 | 一种浅水型小型湖库水华预测方法及预测模型 |
CN109711066B (zh) * | 2018-12-28 | 2023-01-31 | 南开大学 | 一种浅水型小型湖库水华预测方法及预测模型 |
CN113711087A (zh) * | 2019-04-15 | 2021-11-26 | 国立研究开发法人理化学研究所 | 环境因子预测装置、方法、程序、学习完毕模型以及存储介质 |
CN113711087B (zh) * | 2019-04-15 | 2023-11-14 | 国立研究开发法人理化学研究所 | 环境因子预测装置、方法、程序、学习模型及存储介质 |
CN112989281A (zh) * | 2021-02-20 | 2021-06-18 | 中国科学院南京地理与湖泊研究所 | 一种基于遥感藻总量的藻华预测方法 |
CN112989281B (zh) * | 2021-02-20 | 2023-09-12 | 中国科学院南京地理与湖泊研究所 | 一种基于遥感藻总量的藻华预测方法 |
CN113158762A (zh) * | 2021-02-23 | 2021-07-23 | 中国科学院南京地理与湖泊研究所 | 大型湖库滨岸带蓝藻水华堆积风险评估方法 |
CN113158762B (zh) * | 2021-02-23 | 2024-05-28 | 中国科学院南京地理与湖泊研究所 | 大型湖库滨岸带蓝藻水华堆积风险评估方法 |
CN114324796A (zh) * | 2021-12-28 | 2022-04-12 | 中国科学院南京地理与湖泊研究所 | 一种水草上附着丝状藻类水华的预测方法 |
CN114324796B (zh) * | 2021-12-28 | 2023-03-10 | 中国科学院南京地理与湖泊研究所 | 一种水草上附着丝状藻类水华的预测方法 |
CN114485781A (zh) * | 2021-12-30 | 2022-05-13 | 常州市环境科学研究院 | 一种用于浅水湖泊的沉浮式蓝藻水华监测***及预报方法 |
CN114169646A (zh) * | 2021-12-31 | 2022-03-11 | 中科三清科技有限公司 | 水华预测方法、装置、电子设备及计算机可读存储介质 |
CN114564883B (zh) * | 2022-02-11 | 2023-05-19 | 水利部交通运输部国家能源局南京水利科学研究院 | 融合机理与ML的湖泊叶绿素a短期集合预报方法及*** |
CN114564883A (zh) * | 2022-02-11 | 2022-05-31 | 水利部交通运输部国家能源局南京水利科学研究院 | 融合机理与ML的湖泊叶绿素a短期集合预报方法及*** |
CN115719174A (zh) * | 2022-10-26 | 2023-02-28 | 生态环境部卫星环境应用中心 | 土地利用类型与蓝藻水华风险定量关系确定方法和装置 |
CN115719174B (zh) * | 2022-10-26 | 2023-05-02 | 生态环境部卫星环境应用中心 | 土地利用类型与蓝藻水华风险定量关系确定方法和装置 |
CN115508287B (zh) * | 2022-10-28 | 2023-03-24 | 生态环境部卫星环境应用中心 | 湖库蓝藻水华提取方法和装置 |
CN115508287A (zh) * | 2022-10-28 | 2022-12-23 | 生态环境部卫星环境应用中心 | 湖库蓝藻水华提取方法和装置 |
CN116485045A (zh) * | 2023-06-21 | 2023-07-25 | 中科三清科技有限公司 | 一种建立水华预测模型的方法、水华预测方法及装置 |
CN118212970A (zh) * | 2024-05-16 | 2024-06-18 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于时序和空间分布的蓝藻预测方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102135531A (zh) | 大型浅水湖泊72小时蓝藻水华预测方法 | |
Rovere et al. | The analysis of Last Interglacial (MIS 5e) relative sea-level indicators: Reconstructing sea-level in a warmer world | |
Winstral et al. | Spatial snow modeling of wind-redistributed snow using terrain-based parameters | |
JP4682178B2 (ja) | 全国合成レーダ雨量を用いた分布型流出予測システム | |
JP2007011582A (ja) | 洪水予測システム | |
CN109871637B (zh) | 一种云天条件下近地面气温估算方法 | |
CN113821925B (zh) | 一种基于水生土三要素的湿地动态边界确定方法 | |
Ibrakhimov | Spatial and temporal dynamics of groundwater table and salinity in Khorezm (Aral Sea Basin), Uzbekistan | |
CN110516856A (zh) | 基于卷积神经网络的估算海洋次表层温度的方法 | |
Bliefernicht et al. | The WASCAL hydrometeorological observatory in the Sudan Savanna of Burkina Faso and Ghana | |
CN113591288A (zh) | 一种基于克里格插值的土壤湿度数据预测方法及装置 | |
Pachri et al. | Development of water management modeling by using GIS in Chirchik river basin, Uzbekistan | |
CN112926468A (zh) | 一种潮滩高程自动提取方法 | |
Wang et al. | Temporal continuous monitoring of cyanobacterial blooms in Lake Taihu at an hourly scale using machine learning | |
CN115293262B (zh) | 一种岸线功能与结构长时序检测方法 | |
Adhikary et al. | Geostatistical analysis of groundwater level fluctuations in the shallow aquifer of northwestern Bangladesh | |
Polat | An investigation of ancient water collection and storage systems near the Karahantepe neolithic site using UAV and GIS | |
CN108875292A (zh) | 基于遥感的流域水文的仿真***及方法 | |
Rasel et al. | Geospatial analysis of groundwater level variations using Kriging Method | |
Liu et al. | Influences of different surveying and mapping methods on fractal characteristics of gully-head shoulder lines | |
Peng | Application of marine remote sensing technology in the development of fishery economy | |
Rejith et al. | A geo-processing modeling of deltaic suspended sediment variability | |
Khemiri et al. | DEM-Based GIS Algorithms and 3D Geospatial Mapping for Creation of Hydrogeological Models Data in Foussana Basin (Central Tunisia) | |
CN114528672B (zh) | 基于3s技术的城市水文站网布设方法及*** | |
Hastings | Comparison of digital terrain and field-based channel derivation methods in a subalpine catchment, Front Range, Colorado |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20110727 |