一种基于遥感影像与降雨监测数据的雨水去向不确定性分析方法及系统与流程

本发明涉及雨水去向分析领域，尤其涉及一种基于遥感影像与降雨监测数据的雨水去向不确定性分析方法及系统。
背景技术：
：目前，我国正处于快速城市化的发展时期，城市区域自然覆被减少，硬化面积和建筑物密度增加，人口和产业高度集中，导致热量的超常排放和下垫面热力性质的改变，这些都影响和改变着城市的气候气象和水循环过程，使得城市化地区降雨量明显多于非城市化地区。与此同时下垫面硬化面积和硬化率快速增加，城市区域不透水面积增加，径流系数增大，雨水量增加，汇流时间变短，在汛期短时间内较易造成地面集水和局部洪灾。当雨水作为一种用来满足人类生产、生活和生态需要的物质资料时，我们就称之为雨水资源，作为城市防洪和缓解水资源危机的一种有效措施，城市雨水资源化利用显得尤为重要。要进行雨水资源化利用，可通过雨水集蓄入渗系统渗透截留，也可通过不同类型源头设施控制利用，也可通过城市排水管网控制利用，而这些必须建立在掌握降雨—径流—排放以及对雨水去向的分析过程的基础上。雨水降落到地面后先通过城市下垫面自然消化吸收，再经过源头设施控制利用，在年径流总量控制率高的城市，小降雨事件一般不再产生外排雨水，逢大降雨或持续降雨事件，则雨水再通过排水管网等方式外排。收集降雨和城市下垫面实时监测数据，则可进行外排雨量估算。但由于城市各类下垫面的径流系数具有不确定性，当进行大范围径流损失估算时，会产生一定的误差，需对综合径流损失量的不确定性进行定量分析评估，从而为海绵城市的雨水资源化利用或内涝防治工作提供支撑。技术实现要素：针对上述问题，本发明提供一种基于遥感影像与降雨监测数据的雨水去向不确定性分析方法及系统，可以实现运用不确定性理论对雨水外排量的估算，为雨水资源化利用或内涝防控提供支撑。本发明的技术方案如下：一种基于遥感影像与降雨监测数据的雨水去向不确定性分析方法，包括以下步骤：1)基于遥感影像，依据数字高程模型中数据点的高程差确定水流方向及上游集水区，并利用上游集水区内部和上游集水区之间的高程数据按照设定阈值提取水系的高程数据点，建立河网空间拓扑关系及编码，进行汇水面分割得到各子汇水面；2)按照各子汇水面的地面种类加权平均计算得出综合径流系数，并依据区域总降雨量以及综合径流系数得到综合径流损失量；3)对不同下垫面的径流系数进行抽样计算，并依据抽样计算结果定量计算各类下垫面的径流系数不确定所造成的综合径流损失量误差，得到校正后综合径流损失量；4)根据区域总降雨量、校正后综合径流损失量和源头设施的控制降雨量，得到外排雨水量，所述源头设施的控制降雨量为年径流总量控制率乘以区域总降雨量。进一步地，采用arcgis的水文分析方法，通过数字高程模型(dem)提取地表水径流模型信息，进行所述汇水面分割。进一步地，依据《gb50014—2006室外排水设计规范》，计算所述综合径流系数其中si为单一下垫面种类的面积，为单一下垫面种类的径流系数，s为红线内场地总面积，n为下垫面数量。进一步地，采用envi软件创建regionofinterest(roi)训练样区方法，得到不同所述下垫面的具体分布及其面积。进一步地，采用的是拉丁超立方抽样方法(lhs)，定量各类下垫面的径流系数不确定所造成的误差，包括以下步骤：1)设定对n种下垫面的径流系数区间范围进行区间划分，将每一维拆分成互斥的m个区间，使得每个区间有相同的概率；2)分别针对每种下垫面的每个区间中，均随机抽取一个值，组成一个m×n的径流系数抽样结果矩阵；3)对该矩阵的每一列中的元素，在其所在列内进行随机排序，形成一个新的m×n的矩阵，其中每一行即为一个径流系数组合。进一步地，根据《海绵城市建设技术指南》，得到所述年径流总量控制率。进一步地，所述年径流总量控制率的计算方法为选取多年的日降雨资料，将降雨量日值按雨量大小分类，过小的降雨由于不产生径流而忽略不计，统计小于某一降雨量的降雨事件在总降雨记录中的比例，该比例即为年径流总量控制率。进一步地，所述多年的日降雨资料为选取至少30年的、反映长期的降雨规律和近年气候的变化的且不包括降雪的日降雨资料。进一步地，区域总降雨量减去校正后综合径流损失量小于源头设施的控制降雨量，则外排雨水量为零。一种基于遥感影像与降雨监测数据的雨水去向不确定性分析系统，包括：1)基于遥感影像的综合径流损失估算模块，用以处理遥感影像、基于遥感影像分割汇水面、计算汇水面综合径流系数以及综合径流损失量；2)不确定性分析模块，用以对不同下垫面的径流系数进行抽样计算，并依据抽样计算结果定量各类下垫面的径流系数不确定所造成的误差，估算校正后雨水综合径流损失量；3)源头设施的控制降雨量估算模块，用以依据年径流总量控制率和区域总降雨量估算源头设施的控制降雨量；4)外排雨水量估算模块，用以计算总雨量与校正后雨水综合径流损失量及源头设施的控制降雨量的差值，估算外排雨水量。本方法运用不确定性理论，按照径流损失量、受到海绵城市的源头控制设施控制或雨水外排等不同的去向，计算雨水的量，并给出相应的不确定性评估。从而为海绵城市的雨水资源化利用或内涝防控提供支撑。附图说明图1北京地区遥感合成影像图。图2北京地区不同下垫面监督分类结果图。图3雨水去向分析结果的可利用柱状图。具体实施方式为了使本发明的目的、技术方案及优点更加清晰，结合以下北京地区实施例，对本发明进行进一步详细阐述。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。一、基于遥感影像的综合径流损失量估算综合径流损失的雨水，主要是由于雨水通过地表下渗或截留作用造成的，在计算中主要通过综合径流系数来反映损失水量的占比。由于汇水单元在水文统计时是一个相对独立的区域，为了能够较为精确地分析径流情况，首先需要进行汇水面分割；第二，不同土地利用类型的径流系数不同，需通过遥感影像识别出不同土地利用类型，再按照子汇水面的地面种类加权平均计算得出综合径流系数，以此估算出雨水综合径流损失量。1、遥感影像处理方法原理：遥感影像是运用人造地球卫星、航空等平台上的传感器或遥感器在远离目标和非接触目标物体条件下，对目标地物电磁波的辐射、反射特性进行探测，利用不同波段对物体不同的穿透性，达到不同颜色的成像，以此来识别出不同的地物的方法。本方法所运用到的美国航天局陆地卫星landsat8携带有oli陆地成像仪和tirs热红外传感器。某些特殊的光谱波段可排除大气干扰，看到一些特殊的地物特征，比如在近红外波段(nir)通道下，植被的反射率非常高，可利用该波段有效监测植被；短波红外波段(swir)能反应出裸土表面的湿度情况，则可利用该波段有效监测裸土。在雨水径流估算中，无需精度很高的地物识别，主要为了识别城市的不透水地面和绿地等。可通过landsat8传感器数据的不同波段组合进行识别，具体运用到的波段参数说明见下表：landsat8数据部分波段参数波段波长范围(μm)空间分辨率(m)蓝波段0.450～0.51530绿波段0.525～0.60030红波段0.630～0.68030近红外波段0.845～0.88530短波红外11.560～1.66030短波红外22.100～2.30030请参考附图1，关于土地利用识别方面，以北京为例，使用landsat8影像绿光波段、红光波段、红外光波段的北京区域进行合成。合成影像能够直接表现北京地区的下垫面情况。影像中深色表示植被与裸地、浅色表示不透水面(道路、建筑等)、黑色表示水体。请参考附图2，利用遥感图像处理平台envi创建regionofinterest(roi)训练样区，对不同类型下垫面进行监督分类，可得到不同下垫面的具体分布及其面积。2、汇水面分割汇水面是指流经其中的水流和其它物质从一个公共的出水口排出从而形成一个集中的排水区域。在水文分析中，研究区经常要基于这种更小的流域单元进行分析，因而需要对研究区进行汇水面分割，确定小级别流域的出水口的位置，从而精确分析该子流域单元内的径流情况。运用arcgis软件的水文分析模块，通过数字高程模型dem(digitalelevationmodel)提取地表水径流模型的水流方向、汇流量积累、水流长度、河流网络(包括河流网络的分级等)，最终对北京地区流域进行分割。原理：从dem中自动提取流域自然水系是基于水是沿斜坡最陡方向流动的原理，依据dem中数据点之间的高程差来确定水流方向，然后根据水流方向计算每一个数据点的上游集水区；再利用集水区内部和集水区之间的高程数据，通过设置阈值提取所属水系的高程数据点；最后基于水流方向数据，从源头追溯出整个水系，同时划分子流域，建立河网空间拓扑关系及编码。具体操作：首先对原始dem数据进行洼地填充，利用水文分析(surfacehydrologicanalysis)模块中填洼(fillsinks)功能，把洼地点的高程升高，与洼地边缘的最低点的高程相同，得到无洼地的dem；继续应用surfacehydrologicanalysis模块下的流向功能(flowdirection)，生成格网流向图，每个格网值表示它流向相邻像元的方向；汇流网络提取依次执行surfacehydrologicanalysis模块中汇流累积量(flowaccumulation—stream)，河网结构定义(definition—stream)，河流弧段提取(segmentation)一系列功能，提取该流域的一级汇流网络；子汇流区域划分依次执行surfacehydrologicanalysis模块中矢量汇流区生成(catchmentgriddelineation—catchmentpolygonprocessing)功能，划分流域的子汇流区域，并提取子汇流区域边线，生成矢量图。3、汇水面综合径流系数计算方法径流系数说明了流域内自然地理要素将降水量转化为径流量的能力，运用公式则指代任意时段内径流深度d与同时段内降水深度p之比，用符号α表示，即α＝d/p，式中：α为径流系数；d为径流深度，mm；p为降水深度，mm。径流损失量包括下渗、截留雨量，其总量为总雨量×α根据gb50014—2006室外排水设计规范,各种地面种类的径流系数如下表:径流系数汇水面积通常是由各种性质的地面覆盖组成的，随着它们占有的面积比例变化，α值也各异，因此整个汇水面积的径流系数应采用按地面种类加权平均计算得出的综合径流系数。雨量综合径流系数计算公式：式中：为综合径流系数；si为单一下垫面种类的面积，m2；为单一下垫面种类的径流系数；s为红线内场地总面积，m2。以下是常见综合径流系数：综合径流系数4、综合径流损失量计算方法雨水通过地表径流系统下渗、截留部分称为综合径流损失量，通过总雨量和综合径流系数求得。计算公式如下：式中：l为综合径流损失量，mm；r为总雨量，mm；为综合径流系数。二、不确定性分析由于各类下垫面的径流系数的具体数值不能完全确定，因此造成的综合径流损失量具有一定的不确定性，在计算的过程中难免会产生误差。为了能够定量识别这种误差的大小，需要进行不确定性分析。进行不确定性分析采用的是基于拉丁超立方抽样方法(lhs)，它提供了一种有效的从多元参数分布中抽样的方法。它的目的是要保证所有的采样区域都能够被采样点覆盖，也是为了从多维分布中通过抽样得到一个可信的参数值的分布。本发明运用该方法基于取值范围对不同下垫面的径流系数进行抽样计算，具体步骤如下：(1)设定对每种下垫面(假设共有n种下垫面)的径流系数区间范围进行区间划分，将每一维拆分成互斥的m个区间，使得每个区间有相同的概率；(2)分别针对每种下垫面的每个区间中，均随机抽取一个值，组成一个m×n的径流系数抽样结果矩阵；(3)对该矩阵的每一列中的元素，在其所在列内进行随机排序。最终形成一个新的m×n的矩阵，其中每一行即为一个径流系数组合。三、源头设施的控制降雨量估算雨水源头控制是指通过源头分散的小型控制设施，维持和保护场地自然水文功能，有效缓解不透水面积增加造成的洪峰流量增加、径流系数增大、面源污染负荷加重的城市雨水管理理念。根据《海绵城市建设技术指南》对我国各城市年径流总量控制率提出的要求，年径流总量控制率一般在60％～90％之间，可依据不同城市要求判断对应的设计降雨量。源头设施的设计降雨量与年径流总量控制率之间的关系，可参照《海绵城市建设技术指南》，该实施例中，北京地区的年径流总量控制率为85％，设计降雨量为33.6mm。具体方法：选取近30年的日降雨资料，将降雨量日值按雨量大小分类，过小的降雨由于不产生径流而忽略不计，统计小于某一降雨量的降雨事件在总降雨记录中的比例，该比例即为年径流总量控制率，对应降雨量为设计降雨量，源头设施的控制降雨量就等于年径流总量控制率乘以区域总降雨量，公式如下：h＝α×f式中：h为源头设施的控制降雨量，mm；α为年径流总量控制率，％；f为总降雨量，mm。四、外排雨水量估算外排雨水量等于总雨量减去前两项(综合径流损失量和源头设施的控制降雨量)之和。若减去综合径流损失量后的雨量小于源头设施的控制降雨量，则外排雨水量为零；式中：w为外排雨水量，mm；r为总雨量，mm；为综合径流系数；h为源头设施的控制降雨量，mm。请参考附图3，以北京三个汇水区(下垫面情况各不相同)的一场降雨量为100mm的降雨为例，使用本方法进行计算，所得到的带有不确定性的雨水去向分析结果可利用柱状图表示，综合径流损失量、源头设施的控制降雨量和外排雨水量三者在数值上此消彼长，但维持总量不变。以上实施例仅用以说明本发明的技术方案而非对其进行限制，本领域的普通技术人员可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明的精神和范围，本发明的保护范围应以权利要求书所述为准。当前第1页1 2 3