城镇扩张背景下湘江中游水系变化及其响应研究_宣润泽.pdf

第4 1卷第5期2023年5月水 电 能 源 科 学W a t e r R e s o u r c e s a n d P o w e rV o l.4 1 N o.5M a y 2 0 2 3D O I:1 0.2 0 0 4 0/j.c n k i.1 0 0 0-7 7 0 9.2 0 2 3.2 0 2 2 1 1 5 9城镇扩张背景下湘江中游水系变化及其响应研究宣润泽1,杨 柳1,2,3,冯 畅1,2,3,詹 琴4,庾艳花1,吕 凌5(1.衡阳师范学院地理与旅游学院,湖南 衡阳 4 2 1 0 0 2;2.湖南省环境教育与可持续发展研究基地,湖南 衡阳 4 2 1 0 0 2;3.古村古镇文化遗产数字化传承协同创新中心,湖南 衡阳 4 2 1 0 0 2;4.吉林师范大学旅游与地理科学学院,吉林 四平 1 3 6 0 0 0;5.湖南省水文水资源勘测局衡阳分局,湖南 衡阳 4 2 1 0 9 9)摘要:鉴于城镇扩张使大量水域被挤占压缩甚至消失,导致洪涝风险剧增。以湘江中游典型城市衡阳为例,基于2 0 0 12 0 1 6年研究区水系特征数据与城镇化数据,构建了水系变化与多维城镇发展评价体系,剖析了城镇扩张背景下水系变化响应。结果表明,湘江中游地区城镇发展可分为三阶段,城镇扩张呈现速度先快后慢、强度相对较低的特征;研究期内水系数量总体呈衰减趋势,河道形态整体呈萎缩状态,水系数量少、结构单一的流域衰减剧烈;城镇扩张与水系变化关系紧密,其对水系数量变化影响大于结构变化,且与水系数量少、结构单一的流域关联度高。由此可见,城镇建设与扩张等人类活动要尽可能减少对此类流域的不良影响。关键词:城镇扩张;水系变化;响应研究;湘江流域中游中图分类号:F 2 9 9.2 3;T V 8 2 文献标志码:A 文章编号:1 0 0 0-7 7 0 9(2 0 2 3)0 5-0 0 3 0-0 4收稿日期:2 0 2 2-0 5-3 1,修回日期:2 0 2 2-0 8-1 4基金项目:国家自然科学基金项目(4 1 9 0 1 0 2 6,4 2 0 0 1 0 2 4);湖南省自然科学基金项目(2 0 2 2 J J 4 0 0 1 5,2 0 2 1 J J 4 0 0 1 1);湖南省教育厅科学研究项目青年项目(2 1 B 0 6 2 5,2 1 B 0 6 4 6);教育部大学生创新创业训练计划项目(S 2 0 2 0 1 0 5 4 6 0 2 5)作者简介:宣润泽(1 9 9 8-),女,硕士研究生,研究方向为城镇化下水系适应性,E-m a i l:x r z 2 6 31 6 3.c o m通讯作者:杨柳(1 9 8 6-),女,博士、副教授、硕导,研究方向为水文与水资源环境、防灾减灾等,E-m a i l:p a l a p a l a y l 2 0 0 61 2 6.c o m1 概况衡阳作为湘江中游流域典型区,其位于湖南省 中 南 部(2 6.1 2 2 7.4 6 N、1 1 0.5 3 1 1 3.2 8 E),面积1 5 3 1 0 k m2。湘江在境内长2 2 6 k m,占总里程的3 9.7%1,呈S型贯穿其中。区内主要支流有洣水、蒸水、耒水、舂陵水,水系形态呈向心型树枝状辐聚式2。区内降水集中于春夏季,多年平均降水量1 6 9 8.6 mm3,属亚热带季风气候。以盆地地形为主,地势西南高东北低。衡阳为环长株潭都市圈的重要城市,城镇快速扩张,城镇化率从2 0 0 0年的2 9.6 0%增长到2 0 2 1年的5 5.2 3%4。城镇扩张大量侵占林地、水域等蓄水面积,极端强降水事件频发,使区内洪涝风险剧增。现有研究表明5,湘江中游为洪水灾害高风险区,其危险度与易损度均较高,因此开展城镇扩张背景下湘江中游水系变化及其响应研究具有现实意义和迫切性。为此,本文基于2 0 0 12 0 1 6年水系特征数据与城镇化数据,构建了水系变化与多维城镇发展评价体系,剖析了湘江中游城镇扩张与水系变化的联动关系,旨在为推进“绿色城镇化”和防洪减灾提供科学支撑。2 数据来源、处理及研究方法2.1 数据来源和处理水系 数 据 源 于 地 理 空 间 数 据 云(h t t p s:/www.g s c l o u d.c n),2 0 0 1、2 0 0 6、2 0 1 1年L a n d-s a t 4-5 TM影像分辨率为3 0 m,2 0 1 6年L a n d s a t 8 O L I-T I R S影像分辨率为1 5 m。城镇发展数据来自 湖南统计年鉴。城镇新增面积主要指实际城乡、工矿、居民用地面积,来自中国科学院资源环境科学与数据中心(h t t p s:/www.r e s d c.c n)土地利用现状遥感监测数据。为获取科学准确、可对比的水系数据,对所获取的影像进行大气校正、辐射定标、镶嵌裁剪等预处理。运用最邻近分配法6将2 0 1 6年1 5 m分辨率影像重采样为3 0 m,统一遥感影像分辨率。采用监督分类最大似然法7提取水系数据。第4 1卷第5期宣润泽等:城镇扩张背景下湘江中游水系变化及其响应研究但拓扑检查和实地调查发现该结果存在河流不连续、边界移位等错误。同时较多湖泊因云层、建筑等阴影及水体质量产生颜色差异被误判、漏判。故以11 0 0 0同时期G o o g l e E a r t h高清影像为底图,在A r c M a p中调整河道边界、连接断裂处,手动增添可见河流与湖泊等水体,删减非湖泊斑块。最终得到经反复校准后的高精度水系数据。并用3 0 m D EM数据提取了湘江中游及四大子流域边界线。2.2 研究方法为多维度展现水系变化特征,在湘江中游衡阳段全域基础上进一步细分子流域(洣水、蒸水、耒水、舂陵水)。城镇发展水平选择人口城镇化、经济城镇化和空间城镇化三个维度表征,并用熵值法确定各指标权重8,探究三者的贡献率及其相互作用关系。着重分析城镇新增面积、城镇扩张速度和强度三个指标,进一步评价城镇空间扩张水平。表1为水系与城镇化定量表达的指标体系。表1 水系与城镇化定量表达的指标体系T a b.1 I n d e x s y s t e m f o r q u a n t i t a t i v e e x p r e s s i o n o f r i v e r n e t w o r k s y s t e m a n d u r b a n i z a t i o n目标层准则层指标与计算方法含义水系变化1水系数量水面率Wp=Wp 1+Wp 2,%各河流面积Sa,k m2变化特征河面率Wp 1=(Sa/A)1 0 0%,%区域面积A,k m2湖面率Wp 2=(Aw/A)1 0 0%,%湖泊总面积Aw,k m2河网密度Dr=L/A,k m/k m2河流总长度L,k m2水系结构河网复杂度CR=N L/LmN河道等级数 河流主干道长度Lm,k m变化特征支流发育系数K=Ln/Lm各级支流总长度Ln,k m河网结构稳定度SRt=(Lt/RAt)/(Lt-n/RAt-n)(n 0,tn)SRt:第t年河网过去n年的稳定度Lt、RAt、Lt-n、RAt-n分别为第t年、第t-n年河流总长度(k m)、总面积(k m2)城镇发展1人口城镇化人口结构城镇化率/%城镇人口占总人口的比重定量表达生活质量:城镇居民人均可支配收入,元居民全部现金收入中用于安排家庭日常生活的那部分收入2经济城镇化产业结构:第三产业占G D P比重,%第三产业占G D P比重定量表达经济状况:国内生产总值,亿元一定时期内,生产出的全部最终产品和劳务的价值3空间城镇化Ui:城镇新增面积,k m2实际城乡、工矿、居民用地所达到的范围定量表达城镇扩张速度指数Mi=1 0 0 Ui/(t UL Ai)Ui:时段内扩张数量UL Ai:时段初期面积城镇扩张强度指数Li=1 0 0 Ui/(t TL Ai)Ui:时段内扩张数量TL Ai:区域土地总面积 灰色关联度法可用于分析两种动态事物发展的相关性,探索系统内部发展的主要影响因素。本文选择2 0 0 1、2 0 0 6、2 0 1 1、2 0 1 6年的衡阳段全域与各子流域的河面率(X1)、河网密度(X2)、河网复杂度(X3)、支流发育系数(X4)、湖面率(X5)和城镇扩张指数(Y)组成数列,计算关联度,分析城镇化与水系的联动关系。3 结果与讨论3.1 城镇扩张与发展水平分析图1为多维城镇发展水平测度(图例后所标数字为权重)。图1中,湘江中游城镇发展分为三个阶段:2 0 0 6年前,多维城镇发展水平较低,经济城镇化水平高于空间、人口城镇化,且人口与空间城镇化发展趋势拟合度高。2 0 0 62 0 1 4年,多维城镇发展水平提高,空间城镇化增速快,处领先地位,经济城镇化落后。2 0 1 4年后,多维城镇发展水平高且速度快,经济城镇化反超。为一步评价湘江中游空间城镇化水平在中部地区的地位,参照文献9 研究增加2 0 0 62 0 1 65.55.04.54.03.53.02.52.01.51.00.53430262218141062多维城镇化发展水平 无量纲()各指标累计值 无量纲()2001200620112016年份图1 多维城镇发展水平测度F i g.1 T h e l e v e l o f u r b a n i z a t i o n年数据(表2)进行比较分析,发现湘江中游城镇扩张面积、速度、强度呈正相关;空间城镇化先快后慢的发展趋势与南昌一致。在中部六省中,湘表2 湘江中游城镇的空间变化T a b.2 T h e s p a t i a l e x p a n s i o n o f u r b a n o f X i a n g j i a n g R i v e m i d s t r e a m年份城镇新增面积/k m2扩张速度指数扩张强度指数年份城镇新增面积/k m2扩张速度指数扩张强度指数2 0 0 1 2 0 0 6 1 7.2 6 6 1.8 30.0 2 2 0 0 6 2 0 1 6 2 3 4.9 3 01 0.3 60.1 42 0 0 6 2 0 1 1 1 2 6.0 9 81 2.2 30.1 6 2 0 0 1 2 0 1 6 2 5 2.1 9 68.9 00.3 32 0 1 1 2 0 1 6 1 0 8.8 3 2 6.5 50.1 413江中游城镇扩张面积仅小于郑州、合肥,城镇扩张速度仅小于合肥,扩张强度略弱于武汉、长沙、太原,仅强于南昌。可见,湘江中游的城镇扩张速度快但强度低。3.2 水系变化特征规律研究研究区内蒸水、耒水、舂陵水、洣水子流域面积分别占总面积的1 8.9 4%、1 8.5 3%、8.3 3%、7.0 5%。图2为湘江中游水系数量、结构指标变化率。由图2可知,2 0 0 12 0 1 6年湘江中游水系变化显著,不同时段各流域差异大。(a)WP1(b)Dr(c)CR(d)K3.7200120062006201120112016200120162006年河流流域线9.2200120062006201120112016200120162001年河流流域线8.1200120062006201120112016200120162011年河流流域线15200120062006201120112016200120162016年河流流域线Wp1CRDrK图2 湘江中游水系数量、结构指标变化率F i g.2 T h e c h a n g i n g r a t e o f r i v e r n e t w o r k s y s t e m s i n m i d s t r e a m o f X i a n g j i a n g R i v e r从数 量 变 化 看,2 0 0 12 0 1 6年 湘 江 中 游Wp 1、Wp 2、Dr总体呈衰减趋势,分别减少0.7 8%、2 3.2 8%、3.1 1%,与长沙相比,变化趋势相同但衰减