毕业论文-基于DEM的南水北调中线水源区降水量空间插值方法研究

第第页基于DEM的南水北调中线水源区降水量空间插值方法研究作者：指导老师：摘要：本研究的目的是比较不同空间插值方法的插值效果，为分析南水北调中线水源区降水空间分布特征奠定基础。采用南水北调中线水源区29个水文及雨量站点20a（1991－2010年）的降水量观测资料以及高程数据，运用GIS技术，进行空间插值运算，建立20年年平均降水量与站点坐标、高程之间的回归方程。选取普通克里格法、反距离加权法、协同克里格法等空间插值方法，对南水北调中线水源区29个站点的多年平均降水量进行空间插值计算，采用5个站点的降水量数据对插值结果进行验证，并进行对比误差分析，根据插值结果，进一步分析南水北调中线水源区降水空间分布特征。结果表明，考虑了地形因素的协同克里格方法插值效果最好，且较为真实地反映了南水北调中线水源区降水空间分布特征，普通格里格方法次之，反距离加权法最差。插值结果表明，南水北调中线水源区降水呈现明显的纬度地带性分布特征，且降水受地形、地貌的因素影响较大。得出的结论是运用空间插值技术，应综合考虑地形因素，这样可以提高研究区域降水量空间插值结果的精度，并可为分析该区域降水量时空结构特征提供更为精细的数据，为南水北调中线工程的调配和运行提供技术支持。关键词：南水北调中线水源区；降水量；空间插值0引言降水数据是研究区域水文、水资源的基础资料[1]，是农林、生态、地理、全球变化等学科研究的重要数据来源，也是水文分析模型的重要输入因子[2-5]。准确地获取降雨量的空间和时间分布，是进行水资源管理，特别洪水预报、水质管理、径流模拟的重要环节。然而，由于经济和人力、自然环境等原因，降水量观测点的数量极其有限，且分布极为不均，采用有限的站点观测数据来研究区域降水的空间整体分布具有重要意义。利用GIS技术，对地面离散的降水观测点进行空间插值，来得到其他未知区域的降水空间分布结果，已经成为国内外研究的重点和热点[6]。利用有限的观测点来估算未知点的降雨量，常见方法有协同克里格算法(CoKriging)、反距离加权算法（InverseDistanceWeighted，IDW）、张力样条函数（TheTensionSplineFunction）等[7]。由于研究区域和时间尺度的差异，不同插值方法的适应性不同，插值结果和精度也存在较大差异。研究区域越大，地形结构越复杂，涉及的气象类型也就越多，在研究其降水空间插值时，还应该充分考虑插值过程中对降水有重要影响的因子（如地形、风速等）[8]。庄立伟等利用东北日降水数据进行空间插值方法比较，认为IDW插值方法优于克里金插值法[9]。周锁铨等在对长江流域降水数据进行插值研究中，运用了逐步插值法（StepwiseInterpolationAssociation，SIA），并结合GIS技术和多元逐步回归方法，提高了降水空间插值结果的精度[10]。刘志勇等基于DEM的榆林市降水量空间插值方法分析，运用GIS空间插值技术，并综合考虑地形因素，可以提高研究区降水空间插值结果的精度，并为分析区域降水量时空结构特征提供更为精细的数据[11]。本研究区选在南水北调中线工程的水源地。南水北调中线工程主要向输水沿线的河南、河北、北京、天津等四省市20余座城市提供生活和工业用水。南水北调工程是缓解我国北方水资源严重短缺局面的重大战略性工程，对华北地区的人民的生产和生活有着重大影响。而水源地的水质、水量，决定了供水的水质和水量。因此，研究分析南水北调中线水源区降水量的空间分布，有利于进一步掌握供水能力，为政府制定决策提供技术和数据方面的参考。1研究区域概况及数据来源 1.1区域概况南水北调中线水源区地域地跨湖北、河南、陕西三省，位于东经105°-113°，北纬30°-36°，地势由西到东逐渐降低，海拔高度在72.2m至2212.5m之间。地貌有南阳盆地、秦岭山脉、伏牛山系、江汉平原等主要类型，平均海拔为696m。该地区地跨暖温带季风气候和北亚热带两大自然单元，具有明显的过渡特征：冬季寒冷少雨雪，春短干旱多风沙，夏季炎热，秋季晴朗，历年平均降水量为804mm。1.2降水量资料及地形数据收集整理南水北调中线水源区29个降雨量观测点20年（1990—2010年）的年平均降水量以及高程数值，通过对各个站点年降水量资料原始数据的统计处理，并对其准确性进行验证,数据的高程值作为参考因素。如表1所示：表1降水站点基本资料站名所在省份经度纬度海拔（m）西峡河南111.466733.2833224南阳河南112.516732.4333127安康陕西109.016732.7000249镇坪陕西109.516731.88331536三门峡河南111.200034.7667424卢氏河南111.033334.0500869栾川河南111.600033.78331163太白陕西107.316734.05002140永寿陕西108.133334.6833969武功陕西108.200034.2500450西安陕西108.933334.3333392华山陕西110.133334.5500503秦都陕西108.700034.3167393华县陕西109.766734.5000397略阳陕西106.150033.31671434留坝陕西106.916733.61671209汉中陕西106.016733.0667512佛坪陕西107.983333.51671008商县陕西110.016733.8167939镇安陕西109.150033.4167634商南陕西110.866733.5167866宁强陕西106.250032.8167954石泉陕西108.233333.0333380镇巴陕西107.883332.51671283宜昌湖北111.283330.683396老河口湖北111.666732.3833114十堰湖北110.783332.6167297郧县湖北110.800032.8333334郧西湖北110.416732.98332992研究方法2.1方法简介地统计分析是借助相关变异函数，以区域化变量为基础，分析既具有随机性又具有结构性，或空间相关性和依赖性的自然现象的一门学科。其前提理论基础包括前提假设、区域量变化、半变异分析和空间估值。空间插值分为确定性插值、地统计插值。确定性插值分为全局性插值和局部性插值，局部性插值可以分为反距离加权插值、径向基插值、局部多项式插值等。地统计插值又可以分为：普通克里格、简单克里格、泛克里格、协同克里格等。这些方法的体系和关系如图1所示。图1空间插值方法分类图2.2数据探查（预处理）数据分析工具可以让人们更加全面的了解所使用的数据，以便选取合适的参数及方法。如，数据是否服从正态分布，是否存在某种趋势等。在地统计分析里，首先对29个降水观测点的降水量数据（图2）进行数据探查,同时加载湖北省、河南省、陕西省三省的行政县界图，导入观测点。经对比发现，观测点分布大致呈正态分布。如下图3所示：图2南水北调水源区各地历年降水平均值（单位：mm）图3数据探查2.3反距离加权插值反距离加权插值算法（IDW）是基于相似相近的原理--即两个物体相离越近，它们的性质就越相似；反之，离得越远，相似性越小。它以插值点与样本点间的距离为权重进行加权平均，离插值点越近的样本点赋予的权重越大。反距离加权算法的一般公式如下：。（1）其中，Z（s0）为s0处的预测值；N为预测计算过程中要使用的预测点周围样点的数量；为预测计算过程中使用的各个样点的权重，该值随着样点与预测点之间的距离增加而减少；Z(si)是在Si处获得的测量值。2.4普通克里格插值普通克里格（OrdinaryKriging）是区域化变量的线性估计，它假设数据变化成正态分布，认为区域化变量Z的期望值是未知的。插值过程类似于加权滑动平均，权重值的确定来自于空间数据分析。普通克里格插值法以空间统计学作为理论基础，克服了内插中误差难以分析的问题，能对空间分布数据进行最优、线性、无偏内插估计，缺点是计算量大且过程复杂[12]。其计算表达式为：。（2）式中：式中：i为参与插值的站点对估算点降水要素的权重，Xi表示降水站的位置。在克里金插值中，权重不仅建立在已知点和预测点距离的基础上，而且还要依据已知点的位置和已知点值的整体空间分布和排列特征。因此，运用普通克里金插值方法时，权重取决于已知点与预测点之间的距离及空间关系[13]。2.5协同克里格插值协同克里格算法（CoKriging）是统计插值方法中的一种，该方法用一个或多个次要变量对所感兴趣的变量进行插值估算，次要变量和主要变量都有相关关系。并且假设变量之间的相关关系能用于提高变量预测值的精度。用协同克里格方法对降水进行插值时，如果将高程作为一个协变量会获得较好的估值结果。协同克里格插值方法比普通克里格插值方法有2个方面的改进：一是理论上，假设两个随机变量是相关的，通过估算两种估算方法所得到的方差可知，前者在估算精度上比后者有明显提高；二是在采样方面，协同克里格插值方法可以通过较为容易获得的变量对难以得到的变量进行估算，此方法可以提高精度和运算效率。3研究结果与分析3.1三种插值方法插值结果在南水北调中线水源区进行空间插值研究采用ArcGIS软件中的IDW法、普通克里格插值法、协同克里格插值法等3种不同的插值方法，对29个降水量观测站点20年（1991-2010年）的多年平均降水量进行插值，结果见下图4-6。图4基于IDW方法的南水北调中线水源区年平均降水量插值结果（IDW法）IDW法以插值点与样本点间的距离为权重进行加权平均，认为与插值点距离越近的若干样本点的权重越大，其贡献率与距离成反比[14]。该方法简单易行，但是容易出现因极值产生的“牛眼”现象[6]。所以效果并不太理想。图5基于普通克里格方法的南水北调中线水源区年平均降水量插值结果图6基于协同克里格方法的南水北调中线水源区降水量平均值（CoKriging法）3.2三种插值方法插值结果的检验本研究选用空间分布均匀的5个站点的多年平均降水量数据（站点名称如表2所示），对3种插值结果进行准确度检验。从表2和图7可以看出，总体而言，3种插值方法都有较高的插值精度。由于IDW方法是基于贡献率与距离成反比这一原理，因此在样本点处得到的值和原值相差最小。除此之外，从整体上说协同克里格插值法的插值结果与实测值差值的绝对值较小，普通克里格插值法次之，这主要是由于协同克里格插值法在插值过程中充分考虑了地形因素。表2三种插值方法结果及其与原值的误差站名原始值IDW插值结果IDW绝对误差普通克里格插值结果普通克里格绝对误差协同克里格插值结果协同克里格绝对误差秦都520.5580.8660.36603.9183.41624.41103.91华县684.95680.134.82647.0737.88650.0234.93汉中764.14912.52148.38935.90171.76940.24176.10商南1035.25748.02287.23752.73282.52757.95277.30郧县732.3853.63121.33850.73118.43846.38114.08图7普通克里格和协同克里格插值方法结果绝对误差比较图3.3分析从反距离加权算法、协同克里格算法和普通克里格算法的结果图上可以看出，南水北调水源地的空间降水分布，大致呈现西南地区较多，南部的部分地区的年降水量可以达到984mm，西北地区较少，年降水量为657mm。同时，也呈现由山地向平原地区逐渐减少的趋势。反距离加权算法得到的结果在离观测点近的地方，精度较好，但是局部出现极值现象；在采样点较少的地方，精度较低，效果不明显。普通克里格插值算法和协同克里格插值算法得到的结果整体上分布合理，符合自然现象，没有出现反距离加权的局部异常值的现象，是研究本区域较为理想的插值方法。4结论（1）通过数据统计发现，南水北调中线水源区采用的29个站点的年降水量平均值是804mm，年平均降水量最大的是镇巴县的1491mm，年平均降水量最小的是秦都县的520mm，区域内最大降水差约为971mm，分析认为是由于纬度和地形因素影响所致。（2）在利用反距离加权（IDW）算法时出现明显的牛眼现象，如在区域西南角的宁强县，区域中下部的镇巴县，区域东边的商南县和区域右下角的宜昌市等，这些地区的年降水量明显的比较多，就形成了以地区为中心的同心圆状图形。（3）在本区域的插值方法结果研究中发现普通克里格和协克里格两种插值方法的结果很相似，但是，利用5个站点检验发现协同克里格插值法的插值结果与实测值差值的绝对值较小，普通克里格插值法次之，这主要是由于协同克里格插值法在插值过程中充分考虑了地形因素。（4）通过利用29个降水量观测点的20年（1991—2010）降水量平均数据，生成南水北调水源区的降水空间分，再通过对数据的预处理、探查、变换等操作，生成基础数据。利用地统计分析，采用不同的插值方法，得到不同的空间插值图。通过对比分析，发现在这个研究区域里采用协同克里格插值方法得到的插值效果较好，分布基本符合水平地带规律。参考文献[1]陈贺,李原园,杨志峰，等.地形因素对降水分布影响的研究[J].水土保持研究，2007,14(1):119-122.[2]赵传燕,冯兆东,南忠仁.陇西祖厉河流域降水插值方法的对比分析[J].高原气象，2008，27(1):208-214．[3]储少林，周兆叶，袁雷，等．降水空间插值方法应用研究：以甘肃省为例[J]．草业科学，2008,25（6）：19-23.[4]何䶮,傅德平,赵志敏,等．基于GIS的新疆降水空间插值方法分析[J].水土保持研究，2008,15(6):35-37.[5]陆忠艳，李长青，袁子鹏，等．基于GIS的辽宁气温和降水空间扩展方法[J]．中国农业气象，2008，29(1):90-93，14.[6]蔡福,于贵瑞,祝青林,等.气象要素空间化方法精度的比较研究：以平均气温为例[J].资源科学，2005，27(50:173-179.[7]谭继强,丁明柱.空间数据插值方法的评价[J].测绘与空间地理信息，2004，27(4):11-13.[8]甄计国,陈全功,韩涛.甘肃省各流域降水量的GIS模块插值估计与改进[J].气象科学，2009,29(4):467-474.[9]庄立伟,王石立.东北地区逐日气象要素的空间插值方法应用研究[J].应用气象学报，2003，14(5)：605-616.[10]周锁铨,孙琪,肖桐松,等.长江中上游区基于GIS的不同时间尺度降水插值方法探讨[J].高原气象，2008,27(5):1021-1034.[11]刘智勇，张鑫，方睿红.基于DEM的榆林市降水空间插值方法分析[J].西北农林科技大学学报（自然科学版），2010,38（7）:227-234.[12]冯锦明,赵天保,张英娟.基于台站降水资料对不同空间内插方法的比较[J].气候与环境研究，2004,9(2):261-277.[13]邵晓梅,严昌荣,魏红兵.基于Kriging插值的黄河流域降水时空分布格局[J].中国农业气象，2006,27(2):65-69.[14]何红艳,郭志华,肖文发.降水空间插值技术的研究进展[J].生态学杂志，2005,24(10):1187-1191.BasedontheDEMofthesouth-to-northwatertransferprojectcatchmentareasPrecipitationspatialinterpolationmethodsZHANGHai-taoAbstract：Theaimofthisstudywastocomparedifferentspatialinterpolationmethodsofinterpolationeffect,fortheanalysisofsouth-to-northwatertransfermainforestprecipitationspatialdistributioncharacteristicsoffoundation.USESthesouth-to-northwatertransferprojectmainforesthydrologyandrainfallsite20a29(1991-2010)ofrainfallobservationdataandtheelevationdata,usingGIStechnology,spatialinterpolationoperation,build20annualaveragerainfallandtheregressionequationbetweensitecoordinatesandelevation.Chooseordinarykrigingmethod,theinversedistanceweightingmethod,collaborativekrigingmethodforspatialinterpolationmethods,suchas,thesouth-to-northwatertransferprojectin29mainforestsiteyearsofaveragerainfallspatialinterpolationcalculation,usingprecipitationdataoffivesitesoninterpolationresultsforvalidation,erroranalysisandcomparison,accordingtotheinterpolationresults,andfurtheranalysisofsouth-to-northwatertransferprojectmainforestprecipitationspatialdistributioncharacteristics.Resultsshowthatconsideringtheterrainfactorscollaborativekrigingmethodofinterpolationeffectisbest,andtrulyreflectsthesouth-to-northwatertransferprojectma