1986—2020年山东省地表水时空变化特征分析

李艺，范俊甫，张志锟，左吉伟，时宗闻，高宇

(山东理工大学 建筑工程与空间信息学院，山东 淄博 255049)

地表水是维持区域生态平衡、促进经济社会发展的关键因素[1],其作为日益短缺的关键性基础资源,近几年在山东省的总量降幅明显,季节变化显著,供需矛盾日益突出[2]。作为中国典型的沿海省份,目前学者们对该区域的地表水开发利用[3]、水资源管理[4]和水保护[5]已做了大量的研究,开展了一系列基于水文气象站点、雨量数据等资料的定量分析[6],但由于数据的来源渠道多样、标准不规范,且缺乏连续性,因此较难反映研究区地表水的整体时空演化特征。卫星遥感数据能大区域、长时段、高效率监测地表水状态[7],为实现对区域地表水的长时间跨度连续观测提供高可靠的数据源;然而目前研究多局限于局部典型水域或较短时间序列,缺乏对山东省整体地表水长期连续动态变化的研究,且未对地表水的季节性变化特征进行分析比较。从导致地表水面积变化的驱动力方面来看,山东省地表水季节性动态变化受到大气环境影响和社会经济影响等诸多因子的综合影响比较明显,因此单个季度的地表水监测并无法全面反映地表水变化。

今天，当我们提到清江保护立法的时候，这似乎成了一个轻松的话题。而对于生态系统的长远来看，平稳长效的举措更加重要。

获取地面大范围、大尺度的地面信息主要是通过遥感影像,国内外学者在遥感影像提取水体方面已做了大量研究。目前Landsat影像提取方法可大致分为两大类:一类是基于影像波段信息,通过原始波段处理或波段间相组合所构成的指数模型,如McFeeters[8]提出的归一化水体指数(NDWI),虽能提取水体信息但也掺杂大量背景噪音,徐涵秋[9]通过波段变换的方式提出改进后的归一化水体指数(MNDWI),在抑制植被等无效信息方面取得良好效果;另一类是基于图像显示的特征,通过特征分类别判定的分类器法,如Kalke等[10]通过支持向量机的方法对河流水体进行提取,但其中核函数参数的选取有一定的难度,薛源等[11]利用水体指数和决策树结合DEM河网实现对山区河流的自动提取。目前深度学习算法也开始被应用在水体的分类提取当中,但其需要大量准确样本支持的特性是提取工作的痛点[12],因此,高效准确的指数模型法因其普适性依然是当今大范围数据提取工作的首选[13]。

本文基于1986—2020年长时序的Landsat遥感影像数据,通过水体指数模型提取水体以分析山东省地表水的动态度、变化率,揭示了多年山东省地表水的时空演变特征,结合自然资源数据和社会经济数据,对影响地表水面积变化的影响因子进行了相关性分析。该研究旨在揭示山东省地表水时空演变规律及其与气候变化和人类活动的关系,为山东正在实施的区域协调发展战略[14]的顺利实施提供科学的参考依据,为中国沿海地区地表水资源的保护及合理开发利用提供科学参考。

在智慧社区的创建过程中，5G为信息传输提供了高效、可靠的网络基础。加强社区综合服务平台的建设，可以有效地将智慧家庭、物业管理、周边商圈、公共服务统一规划，依靠手机APP、家庭平台等智能终端，形成智慧化、一体化、数字化的社区服务网。

1 研究区与数据来源

1.1 研究区概况

山东省地处中国东部沿海,位于世界公认的黄金纬度的海岸线上(34°25′N～38°23′N, 114°36′E～122°43′E)。地处温带季风型气候,雨热同期,降水集中且季节分配不均,易发生严重春旱和夏涝;加之利用效率低、过度开发等人为因素综合作用,地表水资源短缺逐渐成为常态,给各行业发展和人民的生产生活带来显著影响,成为严重制约工业和农业持续健康发展的瓶颈[15]。

1.2 数据来源与预处理

Landsat系列卫星遥感影像数据来源于美国地质调查局官网;人口和GDP数据来自中国国家统计局;降雨量和用水量数据来自山东省水利厅;耕地数据来自山东省统计局和公开文献[16];气温数据来自欧洲中期数值预报中心的气候再分析数据。

安康市地处秦巴山区，经济发展落后，近30年来农民收入的主要来源是外出务工，一部分女青年因婚姻、工作而流失，导致当地农村男青年的失婚现象更为严重。仅对2010年安康市人口数据进行分析发现，现阶段安康农村20～44岁的男青年比女青年多出近4万人，理论上也就意味着近4万男性可能会沦为“光棍”，这不仅影响到当地人口再生产和社会经济活动的正常进行，还会增加未来政府的负担。

选择1986年、1995年、2004年、2013年、2020年的旱季、雨季共10期Landsat卫星遥感数据作为数据源。旱季选取一月份至三月份、雨季选取八月份中旬至九月下旬,且云量小于5%的数据。为消除差异影响,根据影像特征确定合适的参数和大气模型,对遥感影像进行辐射定标和大气校正等处理后得到可满足研究要求的数据集,相关操作在ENVI 5.3软件中完成。

1.3 研究方法

1.3.1 提取方法

将文献[17]中现有的22种水体指数模型分别应用在研究区地表水提取中,对实验结果进行对比分析发现,MNDWI是分离度最高、提取效果最好的水体指数模型,这与屈慧慧等[18]研究结果一致,其计算原理见式(1)[9]。

近几年全球海平面逐年上升[20],加上研究区北部大面积的养殖池扩建和盐田开发导致海岸线边界不断扩展,海岸侵蚀导致滩涂下边界向内陆推进[21],故本研究将海水从GMNDWI指数提取的研究区地表水的实验结果中进行剔除处理。参照山东省海洋局发布的《海岸线调查技术规范》(山东省地方标准)[22],结合行政区划矢量数据目视解译2020年Landsat遥感影像的海岸线边界,将海水剔除。

(1)

式中:MNDWI为水体指数模型,Green为绿光波段的像元值,MIR1为中红外波段的像元值。该方法能准确提取绝大部分大面积水体,但部分细小河流未被成功提取。采用调整最佳阈值的方法以解决该问题,将初始阈值减小0.01后,水体未被提取的情况得到改善,但部分建筑物、植被等干扰物也被归为水体。针对该问题,本研究引入提取植被效果较好的植被指数模型NDVI[19](式(2))和提取建筑物较为完善的建筑物指数模型NDBI[19](式(3)),分别二值化后构建新型的水体指数GMNDWI,其计算公式见式(4)。

在建设工程项目管理过程中，要将成本计划、成本监督、成本跟踪、成本诊断等工作落实到位，形成成本控制的细节化管理[5]。成本计划是对整个工程的设计、采购、施工等提出科学合理的方案，然后按照方案对施工工作各方面进行管理。在各项工作进行期间，要确保成本监督工作落实到位。及时审核各项费用支出是否合理，符合合同和施工的要求。确保各工程款是否到位，工程是否完成进行具体的施工阶段等等。成本跟踪则是对具体的施工阶段发生的不确定因素如项目变动等进行统计和报告。最后，对整个工程的成本进行诊断分析。确保工程成本是否符合预定计划，如不符合分析其中的具体影响，方便做好下一阶段的战略部署。

(2)

式中:NDVI为植被指数模型,NIR为近红外波段的像元值,RED为红光波段的像元值。

(3)

GMNDWI=MNDWI-NDVI-NDBI。

(4)

1.3.2 时空分布特征分析

时空分布特征分析利用空间网格化处理和分区统计来完成。空间网格化处理就是用规定大小的格网将一定范围的平面进行分割,从而得到单元数据的过程。通过ArcGIS软件,采用1 km×1 km的格网将研究区内旱季水和雨季水进行分割,经过分区统计得到水体面积占各网格的位置及比例,从而分析讨论研究区内旱季水和雨季水的空间分布和时间变化特征。

不同程序中当事人所提供的现有技术证据不尽相同，导致裁决者认识到的专利申请实际上的技术贡献以及对技术细节的理解也不尽相同，对技术方案的最终解释也会产生不可预知的影响。

3.2.1 地表水面积变化率

通过动态度计算来反映地表水面积变化的剧烈程度,分析地表水面积变化的趋势,其计算方法为研究区内研究期始末地表水的面积变化量占初始期地表水面积的比重,并与研究时段的比值。动态度的绝对值越大,表明在一段时间内地表水的面积变化越剧烈。此外,为直观展现研究区地表水面积变化在空间上的分布,采用格网法将研究区分割成1 km×1 km的空间格网,用研究期始末地表水面积占该格网面积的比重变化反映该地区地表水面积变化的趋势,若其值为正值,表明呈增加趋势,反之则为减少趋势。

1.3.4 相关性分析

为探究各影响因子对研究区内旱季水和雨季水的影响程度,分析地表水区域差异特征,运用皮尔森相关系数法将每年地表旱季水和雨季水面积与各影响因子进行相关性分析,通过该系数反映地表水面积与各影响因子的紧密程度。

油茶苗木的种植，对于促进我国农业技术水平的提升具有重要作用，对于改善我国部分地区的农民收入和农村面貌做出了重要贡献。因此，对于油茶的苗木培育和栽培管理技术进行深入研究，不断实践，总结与探讨，促进我国油茶种植技术的不断完善，更好更快地形成有效的油茶产业。

2 精度评价与结果

2.1 提取结果

基于IEET认证的高职专业培养目标的改革——以广东轻工职业技术学院通信技术专业为例………………陈岗,成超，张晋勇，董兵(56)

草莓育苗可在露天的环境下进行，选择相应的匍匐茎，在形成秧苗后与母株分开，形成新的草莓苗株。在母株栽培的过程中，应保证每株之间的距离约为70厘米，及时去除生长出花蕾，保证匍匐茎的快速生长。母株生长的过程中，根据土地的实际情况，针对性地调整施肥和浇水，保证土壤处于湿润的状态，促进匍匐茎的根部生长。由于草莓比较矮小，易受到杂草的影响，在栽培的过程中，除草工作非常重要，通常情况下，每年8～9月为定植期。当草莓的匍匐茎长出幼叶时，为促进匍匐茎生根，幼叶的前端用土覆盖，并喷洒生根剂，在定植期，每株草莓的匍匐茎保留5个作用最佳，若匍匐茎的数量较多，应该及时的去除，确保留下的匍匐茎可以生长成壮苗。

2.2 精度评价

3.2.2 地表水面积动态度

some一般用于肯定句，any用于疑问句、否定句和条件句中。但在疑问句中，当表示说话人希望得到肯定回答或表达请求、建议时应用some。

表1 精度验证结果 单位:%

整体来看,1986—2020年间研究区旱季水面积和雨季水面积的变化趋势大体一致,均为先增加后减少。雨季水面积呈现出波动增加的趋势,面积增加14.2%;旱季水面积呈现出波动减少的趋势,面积减少7.1%;雨季水在1986年时,面积最小为3 238 km2;旱季水在2020年面积最小为2 223 km2。在2004年旱季水面积和雨季水面积均达到顶峰,分别为4 474 km2和2 847 km2。

3 研究区地表水时空演变特征

3.1 研究区域地表水空间区划

为更直观地反应研究区内地表水面积的时空变化,参考水利部水利水电规划设计总院发布的《中国水功能区划》[23]及山东省水利部门发布的《山东省水资源公报》[24],将研究区划分为四部分分别为徒骇马颊河区、花园口以下区、沂沭泗河区和山东半岛诸河区,如图1所示。

图1 地表水分区

3.1.1 地表水空间变化特征

结合图2,从地表水面积的空间占比来看,研究区内永久水(即全年维持水域状态的地表水)和季节水(即随季节变动的地表水)在空间分布上的疏密程度不均,永久水在花园口以下区和山东半岛诸河区的面积占比高于季节水,而徒骇马頬河区和沂沭泗河区内的永久水和季节水占比均等。1986—2020年间永久水的占比有所下降,但仍高于季节水占比。

(a)1986年

3.1.2 地表水时间变化特征

经精度验证发现,本研究所提出方法的精度评价、Kappa系数、总体精度、用户精度均比JRC数据集精度有提升,充分证明了这种水体指数方法精度的提高。

根据地表水区域划分,本文分别统计了各时期各分区内旱雨季水的面积,其变化趋势如图3所示,由图可知,徒骇马頬河区旱季水面积稳步增加,雨季水面积增加后趋于稳定,花园口以下区的旱雨季水面积均呈波动减少的态势,沂沭泗河区与山东半岛诸河区旱雨季水的面积均呈先增加后减少的趋势。

图3 分区地表水面积

3.2 地表水面积测度动态变化

1.3.3 变化趋势分析

分别计算旱雨季水在四个研究时段的变化率及1986—2020年均变化率,以此来分析研究区内地表水面积的变化程度。由图4可知,旱季水增加面积的变化率呈现先减少后增加的趋势,特别是在2004—2013年时间段变化率最小为10%;而旱季水减少面积的变化率在各研究时段较为平稳,保持在15%左右。由此可见,旱季水总体变化率与增加面积的变化率基本保持一致。另外,雨季水增加面积的变化率呈现先减少后增加的趋势,尤其是在2004—2013年时间段变化率最小为6.7%,成为在各研究期变化最为平稳的地表水类型。雨季水减少面积的变化率在前三个研究时间段保持了平稳变化,而在2013—2020年时间段雨季面积减少的幅度高达39.87%,成为在各研究期变化最剧烈的地表水类型。

图4 分区地表水的面积变化率

从地表水面积变化动态度的空间分布(图5)来看,1986—2020年旱雨季水动态度均整体表现为负,其中旱季水动态度在徒骇马頬河区的中部、沂沭泗河区的东部及微山湖水域西部和南部、山东半岛诸河区东北部均表现为负,而在动态度表现为正的区域较为集中如徒骇马頬河区的西北部、山东半岛诸河区的峡山水库和淮河水域、黄河流域花园口以下区的大汶河流域以及沂沭泗河区的沂河流域和微山湖水域。雨季水动态度表现为负的地区包括黄河三角洲地区、沂沭泗河区的微山湖水域西部、山东半岛诸河区的青岛胶州湾等地,而微山湖水域、黄河流域、大汶河流域、潍河流域等地的雨季水动态度则表现为正。整体来看,内陆地表水总体面积呈减少趋势。

表2 研究区地表水面积变化率 单位:%

先以JRC数据集为参考,验证本文水体指数模型方法的精度,利用ArcGIS软件在提取的矢量水体数据内选取随机点,使用eCognition软件对水体提取结果进行基于对象样本的混淆矩阵法精度评价。选取研究区境内的黄河流域对比高分辨率的遥感影像进行人工目视解译,作为精度评价的验证样本,以50 m为最小允许距离,随机选取653个样本点对水体提取结果进行Kappa系数精度、总体精度和用户精度验证,验证结果见表1。

总体来看,在1986—1995年、1995—2004年、2004—2013年三个研究期时段内,旱季水的总体变化率高于雨季水的总体变化率,说明在该时段内旱季水变化较为剧烈,而在2013—2020年时间段内雨季水的变化剧烈程度高于旱季水,详见表2。

(a)旱季水

3.3 地表水变化的驱动力分析

3.3.1 影响因子分析

1)气候变化。山东省地区属温带季风气候,雨热同期,因此选取气温和降雨量2个重要的气候因子分析气候变化对地表水的影响,如图6所示。

图6 气候影响因子年际变化

2)人类活动。综合前人研究,选取人口、GDP、耕地面积、用水量4个社会因子作为指标,用于分析人类活动对山东地表水面积的影响程度。从各指标的变化趋势来看,人口、耕地面积和GDP呈现增长趋势;用水量趋势不明显,在220～260亿m3之间波动,如图7所示。

(a)人口和GDP

3.3.2 相关性分析

由表3可知,旱雨季地表水的面积变化受降雨量、GDP、人口及耕地面积等因素的综合影响作用。从雨季水的影响因素看,降雨量对其面积变化的影响效果最为显著,相关系数达到0.895,造成该现象的主要原因是研究区内雨季降水较为集中。旱季水面积变化的影响因素中,耕地面积的影响程度最大,其相关系数为-0.997,其原因是受到研究区农业结构的影响,冬季大量种植高需水量的小麦、油菜等作物。由于研究区范围较大、地形复杂、人口分布不均匀,其余因素与地表水面积变化的相关性较低,因此其余单一影响因素对地表水面积变化的影响并非是简单的线性关系,其影响效果并不显著。

表3 各影响因子与地表水面积相关性

4 总结与讨论

通过分析新指数方法对研究区地表水面积的提取结果,得出结论如下:

1) 从地表水时空分布来看,研究区地表水空间上时空分布不均衡,大致表现出“东密西疏”的特征。研究区每期旱季水和雨季水面积之比均大致稳定在3∶5。

2) 从地表水测度动态变化来看,研究区1986—2020年地表水面积整体呈先增加再减少的趋势。除徒骇马頬河区西北部沿海地区地表水面积增加外,其他地区地表水面积明显减少。地表水总体面积从2004年开始呈减少趋势。

3) 地表水面积与降雨量和耕地面积的相关系数明显高于其他因素,表明二者的共同作用是导致旱雨季地表水面积变化的关键因素。

总体来说,本研究提出的指数(GMNDWI)对于大多数水体提取效果都较为准确,能基于Landsat遥感影像对研究区地表水进行精细提取,尤其对于面积较大的水体,提取边界极其吻合,准确度高;但影像重访周期较长,无法及时捕捉短期内暴发的洪水,极端天气也一定程度地影响着地表水提取精度。

公式（3）为基于污染物滞留时间的河岸植被缓冲带宽度计算模型，式中，Lp是河岸植被缓冲带除污效率比达到p时，规划河岸植被缓冲带的宽度。

2.2 浓厚的攻防技击意识 “峨眉派拳术带着深厚的攻防技击观念（即带假想敌观念），做到与假想敌进行模拟搏斗，精神还需要保持高度精神集中，技击方法要求瞬间的变换，追求有形打形，无形打影。”[7]通过武术习练者长期的实践，峨眉派拳术技艺一直在探求其本质技击性，同时也是拳师们的大脑这种高级组织的特殊物质的机能长期演变而来的，习练峨眉派拳术技艺十分强调气、力、功法的结合，讲究内外呼吸引导、硬、软、轻、毒的融合，更讲究其功法与门派的结合，塑造峨眉武术的拳种风格无定势，随行多变。

通过本文所用相关性分析指数的结果发现,雨季水增加主要是因为降雨量,而旱季水减少与耕地面积相关性极高,考虑主要是因为农业用水导致。在丰水期,降雨量显著增加,同时由于雨季的灌溉,农业用水量减少,进一步增加了地表水资源量,加上河流湖泊汛期到来,此时易发生洪涝灾害使得人民的生命财产受到损失。在枯水期,由于降水量减少和耕地开发,使得农业用水量显著增加,此时由于水资源的短缺,造成农产品的价格波动,同样会对社会经济的发展造成不利影响。

根据本文研究结果中研究区冬末初春降雨稀少、河流干旱,而夏末雨量充沛、洪涝灾害频发的特点,应分周期种植适宜生长的农作物,在夏季重点做好黄河流域、小清河流域和微山湖周边的水资源调节工作以及旱涝灾害的应急预案。根据雨季丰水期地表水时空变化规律和分布规律,建议着重在黄河流域加固堤坝以保障财产安全。在枯水期做好防旱工作,适当增加人工降雨。山东省地域辽阔,由沿海向内陆延伸的过程中气候及地域差异明显,加上各地区经济水平发展差异较大,因此地表水存在很大的空间差异,基于此本文提出地表水分区分级的治理方案,在一级分区中按流域划分为四个区,实现研究区内地表水的整体调控和方案部署。在此条件下可按行政区划、经济发展情况、流域分布进行二级划分,在该级别中做好水资源的调配,做好灾情的预警、灾时的救援以及灾后的重建。对易发生灾情的湖泊、河流实行三级划分,由当地政府施行专案措施,采取改进传统的灌溉方式、保护和修复水生态等因地制宜的措施,以促进当地经济的可持续发展。

此外,本文侧重于山东省地表水的宏观研究,在小区域地表水体提取方面存在诸多不足,今后将进一步加强联合分析,开展地方地表水综合研究和区域发展规划,促进地表水资源的合理开发利用。