一种月尺度蒸散发量的定量估算方法技术

本发明专利技术公开了一种月尺度蒸散发量的定量估算方法，步骤为，收集流域内的水文气象数据资料，将所有数据调整到逐月尺度；构建陆地水储量变化数据的多元线性方程，对GRACE重力卫星反演的陆地水储量变化数据进行空间降尺度和滤波处理；获得各子流域的实际蒸散发量月值序列，分析典型流域内水文气候季节值的变化发展趋势和时空演变特征；优选Budyko水热耦合平衡方程，利用陆地水储量变化，构建月尺度流域实际蒸散发定量估算模型；率定月尺度流域实际蒸散发定量估算模型参数，根据率定的模型参数计算全球蒸散发月值序列。本发明专利技术丰富了流域蒸散发估算方法，发布空间分辨率较高的蒸散发月值序列，为科学研究以及防洪抗旱等生产工作提供科学依据。

【技术实现步骤摘要】
一种月尺度蒸散发量的定量估算方法
本专利技术属于水文水资源应用
，特别涉及了一种月尺度蒸散发量的定量估算方法。
技术介绍
全球气候变化加剧了水文极端事件和水文过程的时空变异性，进而引起区域乃至全球性的水资源分布不均问题。另外，高强度人类活动也通过改变流域下垫面特征，影响水文循环过程，导致水资源时空变异性加剧，严重威胁社会经济的可持续发展。蒸散发作为陆面水文过程中最关键的变量之一，其产生过程及演变与气候变化和下垫面覆被演变联系紧密。蒸散发不仅在全球和区域气候中起着至关重要的作用，而且在水资源中也影响着流域地下水储量和地表径流。因此，关于蒸散发估算方法的研究已成为当今水文科学领域的热点问题。由于蒸散发及其时空分布与气象状况、土壤水分和植被覆盖等因素彼此关联又相互制约，难以准确获取，因此，有关蒸散发测定方法和估算模型的研究一直是水文学中的热点问题之一。目前蒸散发估算/模拟方法众多，然而由于时空尺度扩展、物理机制不明、参数较多和不确定性等问题，区域蒸散发的估算和模拟研究仍然处在不断探索和完善中。在逐年和多年平均尺度，基于Budyko假设的水量平衡方法在减少不确定性方面具有优势，如何减少年内时间尺度的不确定性取决于对区域水储量变化的合理估算上。GRACE时变重力场卫星提供地球物质迁移变化的信息，反演的陆地水储量变化为区域实际蒸散发估算提供了有效的数据支撑和技术支持。当前针对陆地水储量及其变化的研究多集中在全球或面积较大的流域，在小流域上由于受其较粗的空间分辨率的限制，造成GRACE反演结果与水文模型结果相差较大。
技术实现思路
为了解决上述
技术介绍
提出的技术问题，本专利技术旨在提供一种月尺度蒸散发量的定量估算方法，结合实测水文气象数据和多源卫星遥感产品，基于水热耦合平衡理论构建月尺度蒸散发的定量估算模型。为了实现上述技术目的，本专利技术的技术方案为：一种月尺度蒸散发量的定量估算方法，包括以下步骤：(1)选择典型研究流域，收集流域内的水文气象数据资料，下载多种卫星遥感产品，包括GRACE重力卫星反演的陆地水储量变化数据，研究区域内及附近监测井观测的地下水提取数据，全球陆面数据同化系统GLDAS的降雨量、土壤含水量、冠层含水量、雪水当量、潜在蒸散发量以及中国自动站与CMORPH融合的逐时降水量0.1°网格数据，对这些数据进行预处理，将所有数据调整到逐月尺度；(2)基于流域实测水文气象数据资料和多种卫星遥感产品，构建陆地水储量变化数据的多元线性方程，将GRACE重力卫星反演的陆地水储量变化数据的空间尺度调整与全球陆面数据同化系统GLDAS数据的空间尺度一致，对陆地水储量变化数据进行空间滤波处理以消除误差，获取高精度陆地水储量变化数据；(3)根据水量平衡原理获得各子流域的实际蒸散发量月值序列，采用多种数理统计方法分析典型流域内水文气候季节值的变化发展趋势和时空演变特征；(4)在典型流域评价多种Budyko水热耦合平衡方程的适用性，优选出最能表达流域水热耦合平衡状态的Budyko水热耦合平衡方程，利用陆地水储量变化，构建月尺度流域实际蒸散发定量估算模型；(5)基于流域实测水文气象数据资料率定月尺度流域实际蒸散发定量估算模型参数，根据率定的模型参数计算全球蒸散发月值序列。进一步地，在步骤(1)中，GRACE重力卫星反演的陆地水储量变化数据的空间尺度为1°×1°，全球陆面数据同化系统GLDAS数据的空间尺度为0.25°×0.25°；在步骤(2)中，将GRACE重力卫星反演的陆地水储量变化数据降尺度至0.25°×0.25°。进一步地，在步骤(2)中，所述陆地水储量变化数据的多元线性方程如下：ΔS＝ΔSM+ΔSWE+ΔGW+ΔSC其中，ΔS为陆地水储量变化数据，ΔSM为土壤水变化量，ΔSWE为雪水当量，ΔGW为地下水变化量，ΔSC为冠层含水量变化量。进一步地，在步骤(3)中，采用的数理统计方法包括线性趋势分析法、Mann-Kendall趋势显著性诊断法、小波周期性分析法和GIS空间分析法。进一步地，在步骤(3)中，根据水量平衡原理获得各子流域的实际蒸散发量月值序列如下式所示：E＝P-ΔTWS-R其中，E为实际蒸散发量，P为降水量，R为径流深，ΔTWS为流域储水量变化量。进一步地，在步骤(4)中，优选出最能表达流域水热耦合平衡状态的Budyko水热耦合平衡方程：其中，E0为潜在蒸散发量，参数ω反映流域下垫面特征。进一步地，在步骤(4)中，构建的月尺度流域实际蒸散发定量估算模型如下：其中，Pe为有效降水，Pe＝P-ΔTWS，参数ω′反映流域下垫面特征。进一步地，在步骤(5)中，采用粒子群优化算法对参数ω′进行优化，得到率定的参数ω′值。采用上述技术方案带来的有益效果：本专利技术结合区域内的历史实测水文气候要素数据和多源卫星遥感产品，提出GRACE重力卫星数据的空间降尺度与误差消除技术，基于水热耦合平衡理论构建月尺度蒸散发的定量估算模型，丰富流域蒸散发估算方法，发布空间分辨率较高的蒸散发月值序列，为科学研究以及防洪抗旱等生产工作提供科学依据。附图说明图1是本专利技术的方法流程图；图2是典型流域多年平均尺度Budyko方程参数与多年平均月尺度Budyko方程参数对比图；图3是本专利技术构建的Budyko月尺度蒸散发模型模拟的典型流域蒸散发与通过月水量平衡获取的蒸散发对比的模拟结果图，包括(a)(b)(c)3幅子图，分别代表黄河、长江、澜沧江源区三个流域。具体实施方式以下将结合附图，对本专利技术的技术方案进行详细说明。一种月尺度蒸散发量的定量估算方法，如图1所示，具体步骤如下。步骤1：选择典型研究流域，收集流域内的水文气象数据资料，下载多种卫星遥感产品，包括GRACE重力卫星反演的陆地水储量变化数据，研究区域内及附近监测井观测的地下水提取数据，全球陆面数据同化系统GLDAS的降雨量、土壤含水量、冠层含水量、雪水当量、潜在蒸散发量以及中国自动站与CMORPH融合的逐时降水量0.1°网格数据，对这些数据进行预处理，将所有数据调整到逐月尺度。在GRACE重力卫星反演的陆地水储量变化数据中，采用美国德克斯大学空间研究中心(CSR,CenterforSpaceResearch)、波茨坦地学研究中心(GFZ,GeoForschungsZentrumPotsdam)和美国宇航局(NASA)喷气动力实验室(JPL,JetPropulsionLaboratory)发布的RL05数据，最大阶数60阶，时间跨度为2003年1月至2015年12月。由于RL05数据产品在处理过程中已经去除了各种潮汐影响(海潮、极潮、固体潮等)以及非潮汐的大气和海洋影响，所以在陆地区域主要反映为水储量变化。针对推求过程中的位模型截断阶数的选择、高斯平滑半径等关键性问题，拟采用以下方法进行。由质量变化引起的球谐系数变化的表达式为：其中，a为地球长半轴；γ为质点距微分单元的距离；θ,λ为地心余纬和地心经度；l,m为球谐展开的阶和次；为完全规格化球谐系数；为完全规格化Legendre缔和函数；ρa＝5517kgm-3为地球平均密度；Δρ(γ,θ,λ)为物质的体密度变化，并假定其表现为地球表面厚度为H的一薄层。由于在地球重力场模型高阶项上，GRACE数据的误差占主要成分，因此需要对重力场模型进行截断l本文档来自技高网...

【技术保护点】
1.一种月尺度蒸散发量的定量估算方法，其特征在于，包括以下步骤：(1)选择典型研究流域，收集流域内的水文气象数据资料，下载多种卫星遥感产品，包括GRACE重力卫星反演的陆地水储量变化数据，研究区域内及附近监测井观测的地下水提取数据，全球陆面数据同化系统GLDAS的降雨量、土壤含水量、冠层含水量、雪水当量、潜在蒸散发量以及中国自动站与CMORPH融合的逐时降水量0.1°网格数据，对这些数据进行预处理，将所有数据调整到逐月尺度；(2)基于流域实测水文气象数据资料和多种卫星遥感产品，构建陆地水储量变化数据的多元线性方程，将GRACE重力卫星反演的陆地水储量变化数据的空间尺度调整与全球陆面数据同化系统GLDAS数据的空间尺度一致，对陆地水储量变化数据进行空间滤波处理以消除误差，获取高精度陆地水储量变化数据；(3)根据水量平衡原理获得各子流域的实际蒸散发量月值序列，采用多种数理统计方法分析典型流域内水文气候季节值的变化发展趋势和时空演变特征；(4)在典型流域评价多种Budyko水热耦合平衡方程的适用性，优选出最能表达流域水热耦合平衡状态的Budyko水热耦合平衡方程，利用陆地水储量变化，构建月尺度流域实际蒸散发定量估算模型；(5)基于流域实测水文气象数据资料率定月尺度流域实际蒸散发定量估算模型参数，根据率定的模型参数计算全球蒸散发月值序列。...

【技术特征摘要】
1.一种月尺度蒸散发量的定量估算方法，其特征在于，包括以下步骤：(1)选择典型研究流域，收集流域内的水文气象数据资料，下载多种卫星遥感产品，包括GRACE重力卫星反演的陆地水储量变化数据，研究区域内及附近监测井观测的地下水提取数据，全球陆面数据同化系统GLDAS的降雨量、土壤含水量、冠层含水量、雪水当量、潜在蒸散发量以及中国自动站与CMORPH融合的逐时降水量0.1°网格数据，对这些数据进行预处理，将所有数据调整到逐月尺度；(2)基于流域实测水文气象数据资料和多种卫星遥感产品，构建陆地水储量变化数据的多元线性方程，将GRACE重力卫星反演的陆地水储量变化数据的空间尺度调整与全球陆面数据同化系统GLDAS数据的空间尺度一致，对陆地水储量变化数据进行空间滤波处理以消除误差，获取高精度陆地水储量变化数据；(3)根据水量平衡原理获得各子流域的实际蒸散发量月值序列，采用多种数理统计方法分析典型流域内水文气候季节值的变化发展趋势和时空演变特征；(4)在典型流域评价多种Budyko水热耦合平衡方程的适用性，优选出最能表达流域水热耦合平衡状态的Budyko水热耦合平衡方程，利用陆地水储量变化，构建月尺度流域实际蒸散发定量估算模型；(5)基于流域实测水文气象数据资料率定月尺度流域实际蒸散发定量估算模型参数，根据率定的模型参数计算全球蒸散发月值序列。2.根据权利要求1所述月尺度蒸散发量的定量估算方法，其特征在于，在步骤(1)中，GRACE重力卫星反演的陆地水储量变化数据的空间尺度为1°×1°，全球陆面数据同化系统GLDAS数据的空间尺度为0.25°×...

【专利技术属性】
技术研发人员：邢万秋，王卫光，董青，丁一民，郑佳重，傅健宇，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32