基于 modis 时间序列数据的洞庭湖区洪水淹没频率研究-许超.pdf

《基于 modis 时间序列数据的洞庭湖区洪水淹没频率研究-许超.pdf》由会员分享，可在线阅读，更多相关《基于 modis 时间序列数据的洞庭湖区洪水淹没频率研究-许超.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第 卷第期 年月灾害学 许超，蒋卫国，万立冬，等基于 时间序列数据的洞庭湖区洪水淹没频率研究灾害学， ， （）： ， ， ， ， ， （）： 基于 时间序列数据的洞庭湖区洪水淹没频率研究许超，蒋卫国，万立冬，王文杰，白雪（ 北京师范大学减灾与应急管理研究院，北京 ； 中国环境科学研究院，北京 ； 北京市石景山区环境监测站，北京 ； 环境保护部机关服务中心，北京 ）摘要：选用 卫星的 合成地表反射率数据（ ）和 合成 数据（ ），构建了 年的时间序列数据集，结合：万数字高程数据，通过遥感信息提取方法对洞庭湖区 年水体淹没情况进行了研究，形成 期水体分布数据，并进行了洞庭湖区淹没频率的计算和变化分

2、析。结果表明：利用遥感技术和 语言实现长时间序列影像数据信息提取是可行的；洞庭湖区水体呈明显季节性变化， 月面积较大， 月次年月面积较小，近 年来水体面积总体呈减小趋势；洞庭湖区约 的区域是永久水面或永久陆地，其余 的区域是水陆交替区，是洪水防控的重点区域；近 年来洞庭湖区水陆交替区域年内水面消退时间提前，滩地出露天数增加，枯水期水量减少明显。关键词： 数据；洞庭湖区；水域；面积变化；淹没频率中图分类号： 文献标志码： 文章编号： （ ） ： 洪水灾害是影响人类生存和发展的最主要、最严重的自然灾害之一，而中国是洪水灾害发生最频繁的国家之一，从 年共发生 次水灾，每年经济损失占国民生产总值的 ，

3、占全国主要自然灾害总损失的 以上， 以上的大城市、半数以上的人口和 以上的工农业产值分布在洪水威胁地区。长江流域是我国第一大流域，是洪水灾害最严重的地区，仅 世纪以来就先后爆发了 年、 年、 年和 年等几次特大洪水灾害，对长江中下游地区造成了严重的经济和社会损失。湖泊是流域内重要的天然水利资源，具有供水、灌溉、调节气候等多种重要功能，同时在洪水调蓄、防洪滞洪方面发挥着重要作用，湖泊水体的分布特征和变化规律直接影响着洪水灾害的加剧和减缓。在气候变化、人类活动等多重因素影响下，湖泊水体的分布特征和变化规律发生着明显变化，从而导致了洪水灾害风险的变化。淹没频率是反映某一位置被水体淹没时间长短的重要指

4、标，淹没频率的变化能够直接体现该位置被水体淹没时间长短的变化，综合分析研究区内每个位置的淹没频率和变化规律，就可以获取研究区湖泊水体分布特征和变化规律，因此加强对湖泊水体的长时间、高频次的监测，分析区域淹没频率及其变化，对洪水灾害风险研究和防灾减灾具有重要意义。随着遥感技术的快速发展，遥感技术在湖泊水体监测中得到广泛应用，由于水体几乎可以全部吸收近红外和中红外波段的能量，而土壤、植被等对这两个波段能量的吸收非常少，因此水体反射率很小，植被和土壤反射率很高，目前通过遥感技术提取水体大多都是基于这个特性实现的，通常有两种方法：单波段阈值法，通过设定阈值提取水体 ；水体指数法，通过差值、比值等方式构

5、建指数，增大水体与其他地表辐射的差异，从而提取水体，主要的指数有归一化植被指数（ ）、归一化差异水体指数（ ）、改进的归一化差异水体指数（ ）、混合水体指数（ ）、开放水体似然性指数（ ）等 ，其中归一化植被指数（ ）应用的最为广泛。收稿日期： 修回日期： 基金项目：国家自然科学基金（ ）；中央高校基本科研业务费专项资金；国家重大科技专项（ ）作者简介：许超（ ），男，河北石家庄人，硕士研究生，主要从事生态遥感监测和洪水灾害风险研究 ： 通讯作者：蒋卫国（ ），男，湖南衡阳人，博士，副教授，主要从事遥感信息识别和洪水灾害与湿地生态研究 ： 期许超，等：基于 时间序列数据的洞庭湖区洪水淹没频率研

6、究对湖泊水体的遥感监测，一方面要求具有较高的空间分辨率，以便能够开展水体的遥感识别，另一方面要求具有较高的时间分辨率，以便获取较高频次的水体数据用于分析湖泊水体变化特征，在众多卫星遥感器中， 具有 的最大空间分辨率，能够满足中大型湖泊水体遥感识别要求，每 可获取一次全球观测数据，时间分辨率较高，能够满足多频次的监测需求，同时其数据成本低廉，社会公开度高，是湖泊水体监测较为理想的数据源。洞庭湖是我国第二大淡水湖，是长江中下游仅存的两个通江湖泊之一，也是历史上多次大型洪水灾害的主要受灾区域，洞庭湖南接湘江、资江、澧水、沅江四水，北通长江，承担着调蓄长江和湖南四水的重要任务，是长江流域极为重要的调蓄

7、滞洪区 。本文以洞庭湖区（地 县区）为研究区，基于 数据，利用遥感技术提取 年洞庭湖区洪水淹没范围的时间序列数据，分析区域水体淹没频率和变化特征，为该地区洪水风险管理、防洪减灾提供参考。 研究区和数据洞庭湖位于湖南省北部，长江南岸，洞庭湖区主要涉及岳阳、益阳、常德个地市，包括岳阳市的云溪区、岳阳楼区、君山区、岳阳县、临湘县、华容县、汩罗市和湘阴县，益阳市的资阳区、赫山区、沅江市和南县，常德市的武陵区、鼎城区、汉寿县、安乡县、津市市、临澧县和澧县，区域面积 。该地区属热带季风湿润气候区，年均气温 ，年均降水约 ，地形中心低四周高，东南西三面环山，最高海拔约 。本文所使用的数据主要有类：遥感影像数

8、据，选用的 卫星遥感数据主要有两种，一是 合成的 分辨率 数据（ ），时相为 年月 年月，共计 幅，另一种是 合成的 分辨率地表反射数据（ ），时相为 年月 年月，共计 幅，这些数据来源于美国航空航天局（ ： ）；行政区划数据，选取岳阳市、益阳市、常德市的 个县级行政区作为研究区，行政区划数据来源于国土资源部第二次全国土地调查数据；数字高程数据，选用的数字高程数据是：万高精度数字高程数据（ ），空间分辨率 ，来源于国家基础地理信息中心。 研究方法和技术路线 数据处理方法在长时间序列数据中，不可避免地存在部分质量较差、不能满足研究要求的数据，针对这种情况可以采用多期数据最值合成的方法对数据进行处

9、理，降低数据质量影响 。本文将每期 数据（ ）作为一组，通过最大值合成法，形成期 数据；将每期地表反射数据（ ）作为组，通过最小值合成法，形成期地表反射数据。 水体提取方法本文采用单波段阈值法和归一化植被指数（ ）法相结合的方法，设置以下条规则用于水体提取，并将两条规则下提取的水体进行叠和，减少自动提取中的漏提现象。（） 值满足 月和 月的 ， 月及 月的 ， 月的 ，则提取为水体。通常情况下 为正值的地表为植被，因此以 为条件提取水体，但受到水生植物、土壤含水及泥沙等因素的影响，水体的 值常常稍高于，杜涛等人通过反复测试并修改阈值，最终得到用 提取水体的合理阈值，即 月和 月的 ， 月及 月

10、的 ， 月的 。（）地表反射数据（ ）第二波段的值小于 ，则提取为水体。水体在地表反射数据第二波段的近红外区（ ）反射率较低，而其他地物在该波段具有较高的反射率 ，因此可以利用此特性区分水体与其他地物；通过多次变化阈值，并与同时期卫星影像进行目视对比，最终确定阈值为 。在以上基础上，为了消除云的影响，以 为阈值分别对地表反射数据的两个波段进行二值化，两个波段 值均高于 的像元定义为云，进行剔除；另外为消除山体阴影影响，通过：万 数据生成坡度数据，将坡度大于 的像元定义为山体阴影，予以剔除 。以上各步骤由于数据序列较长，均通过 语言编写批处理程序实现（图）。图 水体提取技术路线图 水体提取精度评

11、定方法选取美国陆地资源卫星影像数据，通过目视解译和监督分类方法对研究区陆地与水体进行区灾害学 卷图 水体提取精度检验对比图分，并以此为准确的水陆划分结果，通过二者水体提取的总面积、重叠面积等数据的叠加对比，分析基于 的水体提取结果的精度。 淹没频率计算方法淹没频率是指在整个研究期限内，某一位置被洪水淹没的时间与研究总时间的比值，原则上讲应该获取研究区每天的水体淹没空间分布数据，但受到数据时间分辨率、数据质量、天气状况等因素的影响，获取逐日数据非常困难。本文用每期水体淹没范围代替本期与下期数据之间几天的淹没范围，认为在固定的时间间隔内，其淹没范围均与前一期水体淹没范围提取结果相同，并以此构建淹没

12、频率的计算方法如下： （ ） 。（）式中：（ ）表示淹没频率；表示水体数据的期数；（ ）是第期的水体分布数据；（ ）是数据的时间间隔；是研究期限内总天数。 结果分析 水体提取结果通过上述方法，共完成了水体提取数据 期，其中受数据本身质量影响，剔除幅，最终获取水体分布数据 期。 精度分析选取 年月底到月初的 和 年月底到月初的 两期影像数据，通过人工目视解译结合监督分类的方法进行水体和陆地划分，并与对应时相 提取的水体进行叠加对比分析；从叠加分析图上看，误差主要有两个方面： 影像的空间分辨率为 ，远低于 影像，大量面积较小、宽度较窄的河流水系在 影像中很难被识别出来，出现了漏提的情况，我们认为这

13、方面的差异是由数据本身特性决定的，因此不列入精度对比范畴内；受空间分辨率及混合像元的影响，通过 提取水体的过程中将水陆交替区域的部分陆地错误的归入了水体（图）。通过两种影像的水体提取面积对比，两期数据的提取精度分别为 和 （表），考虑到 数据的空间分辨率以及 数据云量的影响，我们认为这样的精度是可以接受的。表 水体提取精度计算表年份 提取水体面积 二者一致的水体面积 陆地面积 精度 水体演变特征从年内变化上看，取 年中各月水体面积平均值，洞庭湖水体淹没面积总体呈波状演替， 月为增长期，水体淹没面积迅速增加， 月为丰水期，水体淹没面积最大，形成波峰；月次年月为减少期，水体淹没面积迅速减少， 月为

14、枯水期，水体淹没面积最小，形成波谷；月水体淹没面积最小，约为 ，月最大，约为 ，是月均值的近倍；一般年份湖南省 月降水量可达到汛期（ 月）降水量的 ，少数年份能达到 ，因此月容易形成淹没面积高值，长江主汛期一般在每年月到月上旬，因此月底、月淹没面积也会达到高值，水体提取结果显示，丰水期月面积最大，月稍有回落，月又有所增加，与汛期淹没面积规律一致（图、图）。图 平均水体面积年内变化特征图 期许超，等：基于 时间序列数据的洞庭湖区洪水淹没频率研究图 年各月湖面水体面积变化图从年际变化上看，以三峡水库蓄水为时间节点，对 年前（ 年）和 年后（ 年）的淹没面积进行比较发现：（）月均淹没面积减少， 年之

15、前洞庭湖区月平均水体淹没面积约 ， 年之后到 年初，月平均水体淹没面积约为 ，减少了约 ；（）丰水期淹没面积减少， 年之前 月的淹没面积明显高于 年之后，其中、月平均水体淹没面积均在 以上， 年后、月平均水体淹没面积明显下降，均未超过 ；（）丰水期和枯水期水体淹没面积差减少， 年前丰水期月平均最大淹没面积约 ，月平均最小淹没面积约为 ，相差约 ， 年后月平均淹没面积最大值与最小值相差约 ，减少约 （图、图）。图 各月平均水体面积年际变化特征图图 年 年月湖面水体变化图表 年洞庭湖区淹没频率统计表淹没频率 面积 比例 灾害学 卷 淹没频率计算结果根据淹没频率计算方法，对研究区范围内各像元进行淹没

16、频率计算，并按照淹没频率（）划分为个等级，即 、 、 、 、 、 、 ，计算结果如表所示。根据计算结果可以判断， 年初 年时间里，研究区平均淹没频率为 ，永久水面（ ）主要分布在东洞庭湖、南洞庭湖的核心区，以及中部大通湖，面积约为 ，占研究区总面积的 ；永久陆地（ ）面积约 ，占研究区总面积的 ，这两部分属于较为稳定的区域，不经历水陆更替，基本没有洪水风险；其余 的区域有可能经历水陆更替，存在一定的洪水风险，应该是洪水风险防控的重点区域，其中约 （总面积的 ）的区域平均每年被淹没的时间在 （ ）以上，这些区域基本都是湖泊中心区和主要河流干流区域；约 （总面积的 ）的区域平均每年被淹没的时间在

17、（ ）之间，这些区域主要是湖泊和主干河流周边湿地，丰水期被水体淹没，枯水期暴露在外，水陆更替较频繁，是最容易受到洪水威胁的区域，同时也是洪水调蓄、防洪滞洪的最重要区域（图）。图 年水体淹没频率分布图年际变化方面依然以三峡水库蓄水为时间节点，对比分析 年以前和 年以后洞庭湖区淹没频率发现：（）永久陆地增大，永久水面和水陆交替区域减小。 年以后永久陆地（ ）面积增加了 ，而永久水面（ ）面积减少了 ，水陆交替（ ）区域面积减少了 （表）。表 年前后淹没频率对比统计表频率 年 年变化面积 比例 面积 比例 面积 比例 （）绝大多数水陆交替区域和少量永久水面的淹没频率降低。水陆交替区域和永久水面区域内

18、共有 的区域出现淹没频率降低的情况，占二者总面积的 ，其中降低超过 个百分点（ ）的面积约 ，主要分布在东洞庭湖、南洞庭湖及西洞庭湖，降低 个百分点（ ）的面积达到约 。淹没频率是指某一位置淹没时间与总时间的比值，能够反映该位置水体淹没的持续时间，三峡工程建成蓄水后，汛末蓄水必然导致下游河流湖泊水情发生变化，根据 年前后淹没频率的对比分析结果可知（图）， 年后水陆交替的绝大部分区域出现了淹没频率下降的情况，这就意味着这些区域的水面消退时间提前，滩地的出露天数增加，与赖锡军 、邹邵林 等人通过水动力模型模拟和水位数据统计的结果相一致；根据对比计算结果，约有 区域滩地露出天数增加超过 ，约 区域滩

19、地出露天数增加 ，约 区域滩地出露天数增加 ，约 区域滩地出露天数增加 。图 年前后淹没频率变化图 结论与讨论本文应用 年 的 数据和地表反射数据提取水体信息，获取了长时间序列（ 期）洞庭湖区水体分布数据，分析了洞庭湖水体年内和年际变化规律，并开展了洞庭湖区水体淹没频率的计算，及三峡工程蓄水前后淹没频率的对比分析。（）利用遥感手段，将 数据用于湖泊水体动态监测具有可行性，应用 和地表反射数据相结合的方式，设定合理阈值水体提取的结果具有较高的准确性和可靠性；在针对长时间序列、多时相遥感数据，通过与 相适应的 语言编写批处理程序，可以实现快速、准确的遥感数据批量处理，具有快速、准确、高效等优越性。

20、（）从年内变化上看，洞庭湖区水体面积呈现明显的季节变化规律，水体面积从月份开始增长， 月是丰水期， 月份开始减少， 月 期许超，等：基于 时间序列数据的洞庭湖区洪水淹没频率研究为枯水期；年际变化方面，年均水体面积和丰水期水体面积均呈减少趋势，丰水期与枯水期之间的面积差呈减小趋势，可以说近 年来洞庭湖水体面积总体呈减小趋势。（）近 年来洞庭湖区约 （ ）的区域是永久水面，约 （ ）的区域是永久陆地，其余 （ ）是水陆交替区，存在洪水风险，是洪水风险防控的重点区域；年际变化方面， 年前后洞庭湖区呈现永久陆地增大趋势，绝大多数水陆交替区的淹没频率降低，表明水陆交替区域年内水面消退时间提前，滩地的出露

21、天数增加，这主要是受长江中下游汛期末段三峡水库蓄水，下游河湖水资源补给不足的影响。遥感技术在湖体动态变化研究中具有明显的优势，可以提供准确的水体空间分布信息，甚至可以模拟水体空间演变的过程，这是通过传统水文数据分析难以实现的。同时也存在缺陷，一方面遥感技术兴起时间相对较短，缺乏数据储备，可获取的数据时相要低于传统水文数据；另一方面遥感数据的时间分辨率往往较低，很难像传统水文数据那样精确到天，甚至是小时。另外遥感数据本身的时间分辨率和空间分辨率也存在着矛盾性，空间分辨率较高的数据往往时间分辨率较低，反之亦然，因此在运用较高时间分辨率遥感数据时分析结果的精确度往往会受到影响，而运用较高空间分辨率遥

22、感数据时，分析结果的时间分辨率较低，其动态变化特征的详细程度就会受到影响。参考文献：陈，史培军自然灾害北京：北京师范大学出版社， 张念强基于 的鄱阳湖地区洪水灾害风险评价南昌：南昌大学， 吴赛，张秋文基于 遥感数据的水体提取方法及模型研究计算机与数字工程， ， （）： 郭铌，张杰，梁芸西北地区近年来内陆湖泊变化反映的气候问题冰川冻土， ， （）： 徐萌，李亚春，曾燕，等苏北大型湖泊水域的 遥感监测气象科技， ， （）： 简兴，张振国，邱银国，等基于 数据的水体提取研究进展廊坊师范学院学报：自然科学版， ， （）： 李辉，李长安，张利华，等基于 影像的鄱阳湖湖面积与水位关系研究第四纪研究， ，

23、（）： ， ， ， ， ， （）： 丁莉东，余文华，覃志豪，等基于 的鄱阳湖区水体水灾遥感影像图制作国土资源遥感， ， （）： ， 马明国，宋怡，王雪梅 年新疆若羌湖泊群遥感动态监测研究冰川冻土， ， （）： 孔凡明，蒋卫国，李京，等基于 的 年泰国洪涝受灾信息提取与分析灾害学， ， （）： 刘勇洪，牛铮基于 遥感数据的宏观土地覆盖特征分类方法与精度分析研究遥感技术与应用， ， （）： 徐涵秋利用改进的归一化差异水体指数（ ）提取水体信息的研究遥感学报， ， （）： ， ， ， （ ）： ， ， ， ： ： ， 莫伟华，孙涵，钟仕全，等 水体指数模型（ ）研究及其应用遥感信息， ， （）： ，

24、 黄昌基于遥感和 的河漫滩洪水淹没分析与建模方法研究上海：华东师范大学， 王亚梅，李忠武，曾光明，等洞庭湖土地利用覆被变化及洪涝灾害研究进展四川环境， ， （）： ， 雷璇，杨一鹏，蒋卫国，等基于 数据的洞庭湖水体面积变化分析水资源与水工程学报， ， （）： 杜涛，熊立华，易放辉，等基于 数据的洞庭湖水体面积与多站点水位相关关系研究长江流域资源与环境， ， （）： 吴赛，张秋文基于 遥感数据的水体提取方法及模型研究计算机与数字工程， ， （）： 李景刚，李纪人，黄诗峰，等 时间序列数据在湖泊水域面积动态监测中的应用研究以洞庭湖地区为例自然资源学报， ， （）： 罗伯良，张超湖南主汛期旱涝与大气环流及春季海温的关系高原气象， ， （）： 赖锡军，姜加虎，黄群三峡工程蓄水对洞庭湖水情的影响格局及其作用机制湖泊科学， ， （）： 邹邵林，刘晓清，刘新平，等三峡工程对洞庭湖区滩地出露天数的影响长江流域资源与环境， ，（）： ， ， ， （ ， ， ， ； ， ， ； ， ， ； ， ， ， ） ： （ ） （ ） ， ： ，