CN111797131B - 一种基于遥感降水产品的极端降水区域频率分析方法 - Google Patents
一种基于遥感降水产品的极端降水区域频率分析方法 Download PDFInfo
- Publication number
- CN111797131B CN111797131B CN202010515991.2A CN202010515991A CN111797131B CN 111797131 B CN111797131 B CN 111797131B CN 202010515991 A CN202010515991 A CN 202010515991A CN 111797131 B CN111797131 B CN 111797131B
- Authority
- CN
- China
- Prior art keywords
- precipitation
- region
- grid
- area
- frequency analysis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001556 precipitation Methods 0.000 title claims abstract description 133
- 238000004458 analytical method Methods 0.000 title claims abstract description 59
- 238000009826 distribution Methods 0.000 claims abstract description 19
- 238000011160 research Methods 0.000 claims abstract description 10
- 238000012360 testing method Methods 0.000 claims abstract description 10
- 238000005315 distribution function Methods 0.000 claims abstract description 8
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 6
- 238000010586 diagram Methods 0.000 claims abstract description 6
- 239000000872 buffer Substances 0.000 claims abstract description 4
- 238000007621 cluster analysis Methods 0.000 claims abstract description 4
- 238000013461 design Methods 0.000 claims abstract description 4
- 238000012216 screening Methods 0.000 claims abstract description 4
- 238000012795 verification Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 39
- 238000004088 simulation Methods 0.000 claims description 11
- 238000005259 measurement Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000009191 jumping Effects 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 239000006185 dispersion Substances 0.000 claims description 2
- 238000011835 investigation Methods 0.000 claims description 2
- 238000013508 migration Methods 0.000 claims description 2
- 230000005012 migration Effects 0.000 claims description 2
- 238000013112 stability test Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 11
- 238000012937 correction Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 125000005842 heteroatom Chemical group 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2462—Approximate or statistical queries
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Probability & Statistics with Applications (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Fuzzy Systems (AREA)
- Computational Linguistics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Operations Research (AREA)
- Artificial Intelligence (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了基于遥感降水产品的极端降水区域频率分析方法，包括以下步骤：S1，首先在研究区或缓冲区范围内的站点气象环境数据基础上，将遥感降水产品根据其适用性和精度进行订正；S2，针对极端降水指标利用订正后的遥感降水产品数据生成降水序列，并筛选可用的网格；S3，根据地理特征和统计特征对研究区内所有的网格进行聚类分析，并在聚类结果基础上进行降水一致区识别；S4，采用拟合优度检验为降水一致区选择最佳的频率分布函数；S5，利用最佳频率分布函数对降水一致区中的网格进行区域频率分析，并采用考虑了站点间相关性的蒙特卡洛模拟对区域频率分析结果进行精度验证；S6，描绘研究区内网格在各种重现期下的极端降水设计值和空间分布图。
Description
技术领域
本发明涉及一种区域频率分析方法，特别涉及一种基于遥感降水产品的极端降水区域频率分析方法，属于水文频率分析领域。
背景技术
在全球变暖背景下，探究极端降水的强度、频率等时空特征对预测及应对洪涝灾害至关重要[1]。传统的极端降水时空特征主要应用基于站点数据的区域频率分析进行研究[2-4]，主要包括四个过程：数据筛选、降水一致区识别、频率分布拟合、区域频率分析。
然而地面气象站点受到地形等客观条件影响[5]，分布稀疏，得到的数据完整度不够，难以全面反映地区内降水的空间分布机制。特别是稀疏的站点降水数据用于区域频率分析时，得到的降水一致区的边界不够准确[6]，另外有些无站点地区无法应用区域频率分析方法。因此，基于卫星的网格降水数据、遥感降水数据或者再分析数据，如热带降雨测量计划(Tropical Rainfall Measurement Mission,TRMM)[7]，全球降水气候学项目计划(Global Precipitation Climatology Project,GPCP)[8]，人工神经网络融合GPCP得到的PERSIANN-CDR产品[9]，融合卫星产品与地面观测的全球陆面同化系统(Global Land DataAssimilation System,GLDAS)[10]等，由于其更高的空间分辨率和空间覆盖度，可以当作站点数据的替代。
虽然降水产品可以在一定程度上弥补站点数据的不足，但它在小空间尺度(如单个网格尺度)和小时间尺度(如日尺度)上的精度不能够满足实际的模拟需求[11-12]。特别是对降水极端事件的敏感性较低，降水量与实际值存在较大偏差[13]。在遥感降水产品在投入水文模型或预测模型中使用前，有必要对其进行适用性分析和订正。
区域频率分析的降水输入从站点数据变为遥感降水产品的过程中，不仅改变了输入数据的类型，也为频率分析带来了一些挑战。在区域频率分析研究的一致区识别过程中，调整区域时的人为操作耗时费力[14]，存在较强的主观性和经验性[15]。人为调整在已有的研究中没有得到明显地改进，特别是随着降水输入的量级增长调整时间会成倍增长。此外在频率分析过程中，关于数据独立的基本假设常常被忽视[16-17]，而遥感降水产品是空间连续的，站点独立的要求在实际中难以满足，忽略站点间相关性会使得分析结果过于理想化。
综上所述，利用遥感降水产品替代站点数据进行区域频率分析的方法仍存在一些需要改进的方面，目前在无站点地区缺乏高效可行的极端降水频率分析方法。
文中涉及如下参考文献：
[1]Kendon E.J.,2018.Why extreme rains are gaining strength as theclimate warms[J]. Nature,563(7732),458-460.
[2]Hosking,J.R.M.,Wallis,J.R.,Wood,E.F.,1985.Estimation of theGeneralized Extreme-Value Distribution by the Method of Probability-WeightedMoments.Technometrics,27, 251–261.
[3]Alvarez,M.,Puertas,J.,Pena,E.,2016.Regional frequency analysis ofextremes precipitations in Northern of Mozambique.Ing.Del Agua,20,29–42.
[4]Liang,Y.Y.,Liu,S.G.,Guo,Y.P.,Hua,H.,2017.L-Moment-Based RegionalFrequency Analysis of Annual Extreme Precipitation and its UncertaintyAnalysis.Water Resour.Manag.,31, 3899–3919.
[5]Worqlul,A.W.,Ayana,E.K.,Maathuis,B.H.,MacAlister,C.,Philpot,W.D.,Leyton,J.M.O.,Steenhuis,T.S.,2018.Performance of bias corrected MPEG rainfallestimate for rainfall-runoff simulation in the upper Blue Nile Basin,Ethiopia.Journal of Hydrology,556,1182-1191.
[6]Yang,T.,Shao,Q.X.,Hao,Z.C.,Chen,X.,Zhang,Z.X.,Xu,C.Y.,Sun,L.,2010.Regional frequency analysis and spatio-temporal pattern characterizationof rainfall extremes in the Pearl River Basin,China.J.Hydrol.,380,386–405.
[7]Zhang,Y.,Hong,Y.,Wang,X.G.,Gourley,J.J.,Xue,X.W..,Saharia,M.,Ni,G.H.,Wang,G.L.,Huang,Y.,Chen,S.,Tang,G.Q.,2015.Hydrometeorological Analysisand Remote Sensing of Extremes:Was the July 2012Beijing Flood EventDetectable and Predictable by Global Satellite Observing and Global WeatherModeling Systems？J.Hydrometeorol.,16,381–395.
[8]Basheer,M.,Elagib,N.A.,2019.Performance of satellite-based andGPCC 7.0rainfall products in an extremely data-scarce country in the NileBasin.Atmos.Res.,215,128–140.
[9]Gao,F.,Zhang,Y.H.,Chen,Q.H.,Wang,P.,Yang,H.R.,Yao,Y.J.,Cai,W.Y.,2018.Comparison of two long-term and high-resolution satellite precipitationdatasets in Xinjiang,China. Atmos.Res.,212,150–157.
[10]Bai,P.,Liu,X.M.,2018.Intercomparison and evaluation of threeglobal high-resolution evapotranspiration products across China.J.Hydrol.,566,743–755.
[11]Seyyedi,H.,Anagnostou,E.N.,Beighley,E.,McCollum,J.,2015.Hydrologic evaluation of satellite and reanalysis precipitation datasetsover a mid-latitude basin.Atmos.Res.,164-165, 37-48.
[12]Tang,G.,Clark,M.P.,Papalexiou,S.M.,Ma,Z.,Hong,Y,2020.Havesatellite precipitation products improved over last two decades？Acomprehensive comparison of GPM IMERG with nine satellite and reanalysisdatasets.Remote Sens.Environ.,240,111697.
[13]Muhammad,E.,Muhammad,W.,Ahmad,I.,Muhammad Khan,N.,Chen,S.,2020.Satellite precipitation product:Applicability and accuracy evaluation indiverse region.Sci.China Technol.Sc.,doi:10.1007/s11431-019-1457-3
[14]Hosking,J.R.M.,Wallis,J.R.,1997.Regional FrequencyAnalysis.Cambridge,United Kingdom.
[15]Sun,H.M.,Wang,G.J.,Li,X.C.,Chen,J.,Su,B.D.,Jiang,T.,2017.Regionalfrequency analysis of observed sub-daily rainfall maxima over easternChina.Adv.Atmos.Sci.,34,209–225.
[16]Yin,Y.,Chen,H.,Xu,C.Y.,Xu,W.,Chen,C.,Sun,S.,2016.Spatio-temporalcharacteristics of the extreme precipitation by L-moment-based index-floodmethod in the Yangtze River Delta region,China.Theor.Appl.Climatol.,124,1005–1022.
[17]Chen,Y.D.,Zhang,Q.,Xiao,M.,Singh,V.P.,Leung,Y.,Jiang,L.,2014.Precipitation extremes in the Yangtze River Basin,China:Regionalfrequency and spatial-temporal patterns. Theor.Appl.Climatol.116,447–461.
发明内容
针对以上站点降水数据不完整、遥感降水产品精度较低以及区域频率分析的应用条件不等问题，本发明提出了一种基于遥感降水产品的极端降水区域频率分析方法，从数据输入的角度，考虑无站点地区的降水数据来源及精度，建立空间连续的高分辨率遥感降水产品与布设稀疏的高精度地面气象站点降水数据之间的关联性，得到基于遥感降水产品的高精度降水输入。此外，在区域调整过程中，围绕统计特征与地理条件等基本因素，将调整过程自动化；在区域频率分析时严格控制基本假设的应用条件，结合物理实际对极端降水进行分析。因此结合遥感降水产品的高精度高分辨率输入以及更加自动化、实际性的分析特点，提出了一种实用性更广泛、客观性更强、可信度更高的极端降水区域频率分析方法。与传统方法相比，本发明首先利用订正后的高分辨率的遥感降水产品有效解决了地面气象站点降水数据不完整的问题，其次提出自动主观调整的思想明显改善了人为调整的性能，最后考虑独立性假设。本发明的核心是尽可能精确地为频率分析提供高精度高分辨率降水输入，高效地为频率分析划分边界明显、空间连续的降水一致区，真实地分析极端降水的频率特征。
本发明提出的技术方案为一种基于遥感降水产品的极端降水区域频率分析方法，包括以下步骤：
步骤1，利用研究区或其缓冲区范围内影响降水的站点气象环境数据与站点对应位置的网格降水值构建K近邻回归模型，预测订正后的遥感降水值；
步骤2，基于订正后的遥感降水数据，采用极端降水指标生成极端降水序列，并进行平稳性检验和独立性检验，筛选能够进行区域频率分析的网格；
步骤3，根据地理特征和统计特征对研究区内所有网格进行聚类分析，并在聚类结果基础上进行降水一致区识别；
步骤4，针对得到的降水一致区，按照拟合优度检验结果选择最佳的频率分布函数；
步骤5，结合降水一致区内的极端降水序列和频率分布函数进行区域频率分析；
步骤6，对所有网格的频率分析结果进行插值，描绘出研究区内网格在重现期下的极端降水设计值和空间分布图。
进一步的，步骤1所述的K近邻回归模型构建方法如下，
以站点气象环境因子，包括经度、纬度、高程、坡度、植被指数、空气温度、相对湿度、待订正的遥感降水以及站点实测降水为样本，按照距离测度计算距离实例i最近的K个实例，将实例i附近的K个实例的输出值的均值等于该实例的预测输出值，按照留一法交叉验证结果，得到最佳K值，从而构建日降水数据订正的KNN回归模型。
进一步的，步骤3中进行一致区识别的具体方法如下：
步骤3.1，对于每个聚类结果，即初始类簇，每个类簇中有N个网格，记录第i个网格的长度为ni，网格i处的样本线性变差系数为t(i)、样本线性偏态系数为t3 (i)、样本线性峰度系数为t4 (i)，计算每个类簇的异质性测度H值，对于H值大于k的类簇记为异质区域r，异质性测度公式如下；
式中，初始类簇的离散度V计算公式如下；
其中tR为区域平均线性变差系数，μV和σV分别代表V的均值和标准差；
步骤3.2，对于异质区域r，计算每个网格的不一致性测度D值并按降序排列，D值越大说明网格在整个区域中越不一致，由此得到需要被移动的网格g；
其中，ui＝[t(i) t3 (i) t4 (i)]T，为异质区域r中的样本线性矩比系数的中心，A定义为/>
步骤3.3，记异质区域r的异质性测度为H，检测网格g的D值是否超过了阈值，如果是，直接删除异质区域r内的网格g，再次计算删除网格g后的区域的异质性测度H′d，如果满足均质性要求，即异质性测度<k，则区域删除操作完成；若H′d＞k&H′d＜H：重复该步骤；若H′d＞k&H′d＞H：撤销所有删除操作，输出提示“该区域是异质的！”，异质区域r的调整结束；如果不是，跳到步骤3.5；
步骤3.4，检查删除的网格是否位于所属异质区域r的边缘：如果是，跳到步骤3.6；否则直接删除网格g；
步骤3.5，检测调整后的区域是否存在网格数少于一定阈值或者区域相互交错的情况：如果是，撤销所有删除操作，对初始类簇执行合并或分割操作；计算异质性测度H′m，是否小于 k，若H′m＜k：合并或分割操作完成，异质区域r的调整结束；记下得到的一个或多个均质区域rs和删除的网格gs；若H′m＞k&H′m＜H：重复步骤3.4；若H′m＞H&H′m＞H：撤销合并或分割操作，输出提示“这些区域都是异质的！”，异质区域r的调整结束；
步骤3.6，移动网格g；计算移出区域的异质性测度Hout和移入区域的异质性测度Hin，若 Hout＜k&Hin＜k：移动操作完成，区域r的调整结束，记下移出区域rout和移入区域rin；若 Hout＜k||Hin＜k：重复步骤3.4；若Hout＞k&Hin＞k：撤销移动操作；
步骤3.7，重复步骤1～6，直到每一个初始类簇都被识别为一个均质区或异质区，记下最终的一致区Homs和最终的删除网格Gs，至此调整操作结束。
进一步的，步骤3.3中检测网格g的D值是否超过了阈值，即D>3或D＞(N-1)/3。
进一步的，步骤3.5中网格数少于一定阈值，即网格数低于研究区域内网格总量的1％。
进一步的，k的取值为2。
进一步的，步骤5中还包括，采用蒙特卡洛模拟对极端降水频率分析结果进行精度评估，具体实现方法方式如下：
步骤5.1，计算网格间的相关性：假设Pi是网格i在第k年的降水量，rij代表网格i和j之间在第k年的样本相关性，为网格i的多年平均降水量，则区域样本相关性ρ根据rij的平均值估计，其中：
步骤5.2，蒙特卡洛模拟设定：模拟区域具有与实际区域相同的条件，包括站点数量、每个点的记录长度、区域平均线性矩比和区域站点间相关性；在第m次模拟中，网格i处的具体不超越概率F处的分位数估计值为估计分位数的相对均方根误差(relativeRMSE) 是在所有模拟的基础上开方并取平均值来估计；Pi(F)为实际区域中网格i处的分位数；对于大于M次的模拟次数，近似的相对RMSE为：
估计分位数的平均相对RMSE反映了该区域所有网格的估计分位数精度。
进一步的，步骤3中所述地理特征包括经度、纬度、高程，所述统计特征包括年平均降水、年降水的标准差。
本发明具有以下优点和积极效果：
(1)适用于气象站点布设少、无站点地区。不同于直接利用地面气象站点降水数据进行分析，本方法采用了订正后的高分辨率、高覆盖率以及高精度的遥感降水产品作为降水输入，规避了站点稀疏地区甚至是无站点地区的降水资料短缺的问题。
(2)区域频率分析的效率提高。本方法考虑了常见的几种调整思想，采用自动化的形式降低了由于降水数据输入的量级的增大导致的调整时间和迭代次数的增长，减少了调整过程的主观性和经验性，提高了降水一致区识别的速度和精度。
(3)区域频率分析结果更加符合物理实际。本方法在处理区域频率分析的基本假设时，充分尊重物理实际，严格控制分析、模拟过程中的站点独立的设定，提高了频率分析结果的可靠性。
(4)区域频率分析所用的数据、方法、特征和知识具有可扩展性。本发明的重点是尽可能准确地反映无站点地区的极端降水的设计值及空间分布，原则上对这一步骤使用到的数据、方法、特征和知识没有限制。
附图说明
图1是本发明实例的研究区；
图2是本发明实例的总体技术流程图；
图3是本发明实例的遥感降水订正前后精度对比图；
图4是本发明实例的区域调整过程图；
图5是本发明实例的极端降水频率分析结果及精度验证结果图，a)-f)表示实例中研究区划分得到的7个降水一致区(记为I区-VII区)；
图6是本发明实例的结果图。
具体实施方式
具体实施时，本发明所提供的技术方法和流程可通过计算机软件技术进行实现，为了使本发明技术更容易理解掌握，下面以具体实施结合附图与实例对本发明作进一步说明：
实施例选取长江流域为研究区域，如图1所示。使用的降水数据包括长江流域50km缓冲区范围内的177个气象站点数据，数据集来源于国家气象科学数据中心。以及GLDAS日降水数据(空间分辨率为0.25°×0.25°)，数据来源于GES DISC网站。实施例选取的数据资料时间范围为1970-2014年。
实施例提出的基于遥感降水产品的极端降水区域频率分析方法采用如下的流程，总体方法如图2所示：
步骤1，首先对GLDAS日降水数据集在长江流域的适用性进行分析，在此基础上对GLDAS日降水数据进行订正，订正模型为由177个地面站点的气象环境数据学习得到的KNN回归模型，实施例选取的气象环境数据为经度、纬度、高程、坡度、植被指数、空气温度、相对湿度、待订正的遥感降水。然后根据每日的KNN模型预测订正后的GLDAS日降水数据，订正前后GLDAS日降水在年最大1日降水事件上的精度如图3所示。
步骤2，基于订正后的GLDAS日降水数据，生成年最大1日降水序列。对该序列进行平稳性检验和独立性检验。本实施例选择改进的Mann Kendall检验检验序列的平稳性，利用滞后数为1的自相关系数检验序列独立性，筛选出95％置信区间下，能够进行区域频率分析的网格。
步骤3，根据网格的地理特征和统计特征进行聚类分析。本实施例选取的地理特征为经度、纬度、高程，统计特征为年平均降水、年降水的标准差，对长江流域所有网格进行K均值聚类，得到7个初始类簇。对初始类簇逐一进行调整，得到长江流域7个降水一致区的过程如图4所示。从图4中的7个一致区的一致区和删除网格的分布可以看出，经过调整的迭代，长江流域被划分为边界明显、空间连续的一致区。较之人为调整方法，自动调整结果的速度快，可重复性高。
步骤4，为长江流域的7个降水一致区选择最佳的频率分布函数，本实施例选择了5种三参数分布(广义Logistic分布GLO、广义极值分布GEV、广义正态分布GNO、皮尔逊III型分布PE3、广义帕鲁托分布GPA)和1种五参数分布(Wakeby分布WAK)作为候选分布，拟合优度检验可得到每个降水一致区的最佳分布。
步骤5，根据长江流域7个降水一致区和区内最佳频率分布函数进行区域频率分析，结合每个一致区内的平均相关系数进行蒙特卡洛模拟。本实施例以区域I为例，说明在考虑相关性(R＝0.301)和不考虑相关性的情况下，均质区域I在重现期为2-100年下的区域增长曲线、相对均方根误差和99％误差界，结果如表1所示。
表1.长江流域区域I的增长曲线值和精度评估结果表
从表1中可以看出，考虑了降水一致区中实际存在的站点间相关性后，相对均方根误差 (RMSE)变大，99％误差界也更宽，但是区域频率分析结果仍在误差界以内，说明频率分析结果在重现期为2-100年情况下都是可靠的；模拟区的平均异质性测度H与实际区的异质性测度相等，说明考虑了相关性结果也更加符合实际情况。对全部降水一致区进行区域频率分析，得到长江流域所有网格的极端降水的区域频率分析结果。
步骤6，通过ArcMap的地统计功能对长江流域所有网格的频率分析结果进行插值，得到重现期为100年(即百年一遇)下的极端降水空间分布图如图6所示。
以上所述仅为本发明的优选实施例而已，不能认定本发明的具体实施只限于此实例，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，都应当视为属于本发明的保护范围。
Claims (7)
1.一种基于遥感降水产品的极端降水区域频率分析方法，其特征在于，包括以下步骤：
步骤1，利用研究区或其缓冲区范围内影响降水的站点气象环境数据与站点对应位置的网格降水值构建K近邻回归模型，预测订正后的遥感降水值；
步骤2，基于订正后的遥感降水数据，采用极端降水指标生成极端降水序列，并进行平稳性检验和独立性检验，筛选能够进行区域频率分析的网格；
步骤3，根据地理特征和统计特征对研究区内所有网格进行聚类分析，并在聚类结果基础上进行降水一致区识别；
步骤3中进行一致区识别的具体方法如下：
步骤3.1，对于每个聚类结果，即初始类簇，每个类簇中有N个网格，记录第i个网格的长度为ni，网格i处的样本线性变差系数为t(i)、样本线性偏态系数为t3 (i)、样本线性峰度系数为t4 (i)，计算每个类簇的异质性测度H值，对于H值大于k的类簇记为异质区域r，异质性测度公式如下；
式中，初始类簇的离散度V计算公式如下；
其中tR为区域平均线性变差系数，μV和σV分别代表V的均值和标准差；
步骤3.2，对于异质区域r，计算每个网格的不一致性测度D值并按降序排列，D值越大说明网格在整个区域中越不一致，由此得到需要被移动的网格g；
其中，ui＝[t(i)t3 (i)t4 (i)]T，为异质区域r中的样本线性矩比系数的中心，A定义为
步骤3.3，记异质区域r的异质性测度为H，检测网格g的D值是否超过了阈值，如果是，直接删除异质区域r内的网格g，再次计算删除网格g后的区域的异质性测度H′ d，如果满足均质性要求，即异质性测度<k，则区域删除操作完成；若H′ d>k且H′ d<H：重复该步骤；若H′d>k且H′d>H：撤销所有删除操作，输出提示“该区域是异质的！”，异质区域r的调整结束；如果不是，跳到步骤3.5；
步骤3.4，检查删除的网格是否位于所属异质区域r的边缘：如果是，跳到步骤3.6；否则直接删除网格g；
步骤3.5，检测调整后的区域是否存在网格数少于一定阈值或者区域相互交错的情况：如果是，撤销所有删除操作，对初始类簇执行合并或分割操作；计算异质性测度H′m，是否小于k，若H′m<k：合并或分割操作完成，异质区域r的调整结束；记下得到的一个或多个均质区域rs和删除的网格gs；若H′m>k且H′m<H：重复步骤3.4；若H′m>k且H′m>H：撤销合并或分割操作，输出提示“这些区域都是异质的！”，异质区域r的调整结束；
步骤3.6，移动网格g；计算移出区域的异质性测度Hout和移入区域的异质性测度Hin，若Hout<k且Hin<k：移动操作完成，区域r的调整结束，记下移出区域rout和移入区域rin；若Hout<k或Hin<k：重复步骤3.4；若Hout>k且Hin>k：撤销移动操作；
步骤3.7，重复步骤1～6，直到每一个初始类簇都被识别为一个均质区或异质区，记下最终的一致区Homs和最终的删除网格Gs，至此调整操作结束；
步骤4，针对得到的降水一致区，按照拟合优度检验结果选择最佳的频率分布函数；
步骤5，结合降水一致区内的极端降水序列和频率分布函数进行区域频率分析；
步骤6，对所有网格的频率分析结果进行插值，描绘出研究区内网格在重现期下的极端降水设计值和空间分布图。
2.如权利要求1所述的一种基于遥感降水产品的极端降水区域频率分析方法，其特征在于：步骤1所述的K近邻回归模型构建方法如下，
以站点气象环境因子，包括经度、纬度、高程、坡度、植被指数、空气温度、相对湿度、待订正的遥感降水以及站点实测降水为样本，按照距离测度计算距离实例i最近的K个实例，将实例i附近的K个实例的输出值的均值等于该实例的预测输出值，按照留一法交叉验证结果，得到最佳K值，从而构建日降水数据订正的KNN回归模型。
3.如权利要求1所述的一种基于遥感降水产品的极端降水区域频率分析方法，其特征在于：步骤3.3中检测网格g的D值是否超过了阈值，即D>3或D>(N-1)/3。
4.如权利要求1所述的一种基于遥感降水产品的极端降水区域频率分析方法，其特征在于：步骤3.5中网格数少于一定阈值，即网格数低于研究区域内网格总量的1％。
5.如权利要求1所述的一种基于遥感降水产品的极端降水区域频率分析方法，其特征在于：k的取值为2。
6.如权利要求1所述的一种基于遥感降水产品的极端降水区域频率分析方法，其特征在于：步骤5中还包括，采用蒙特卡洛模拟对极端降水频率分析结果进行精度评估，具体实现方法方式如下：
步骤5.1，计算网格间的相关性：假设Pi是网格i在第k年的降水量，rij代表网格i和j之间在第k年的样本相关性，为网格i的多年平均降水量，则区域样本相关性ρ根据rij的平均值估计，其中：
步骤5.2，蒙特卡洛模拟设定：模拟区域具有与实际区域相同的条件，包括站点数量、每个点的记录长度、区域平均线性矩比和区域站点间相关性；在第m次模拟中，网格i处的具体不超越概率F处的分位数估计值为估计分位数的相对均方根误差是在所有模拟的基础上开方并取平均值来估计；Pi(F)为实际区域中网格i处的分位数；对于大于M次的模拟次数，近似的相对RMSE为：
估计分位数的平均相对RMSE反映了该区域所有网格的估计分位数精度。
7.如权利要求1所述的一种基于遥感降水产品的极端降水区域频率分析方法，其特征在于：步骤3中所述地理特征包括经度、纬度、高程，所述统计特征包括年平均降水、年降水的标准差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010515991.2A CN111797131B (zh) | 2020-06-09 | 2020-06-09 | 一种基于遥感降水产品的极端降水区域频率分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010515991.2A CN111797131B (zh) | 2020-06-09 | 2020-06-09 | 一种基于遥感降水产品的极端降水区域频率分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111797131A CN111797131A (zh) | 2020-10-20 |
| CN111797131B true CN111797131B (zh) | 2024-02-20 |
Family
ID=72802944
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010515991.2A Active CN111797131B (zh) | 2020-06-09 | 2020-06-09 | 一种基于遥感降水产品的极端降水区域频率分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111797131B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112507052B (zh) * | 2020-12-03 | 2024-08-23 | 清华大学 | 一种大气水循环分析方法、装置、电子设备和存储介质 |
| CN113570104A (zh) * | 2021-03-08 | 2021-10-29 | 中电建华东勘测设计研究院(郑州)有限公司 | 一种流域水生态功能分区方法 |
| CN113407596B (zh) * | 2021-06-23 | 2022-07-15 | 北京华云星地通科技有限公司 | 基于交叉辐射定标的气象遥感数据精度检验分析系统 |
| CN113408211B (zh) * | 2021-07-12 | 2022-02-11 | 中国水利水电科学研究院 | 基于机器学习的城市内涝积水点积水过程预测方法及系统 |
| CN114398592B (zh) * | 2021-12-27 | 2022-09-02 | 中国人民武装警察部队警官学院 | 一种基于异质格网高程线性分解模型的高程值解算方法 |
| CN114677059B (zh) * | 2022-05-26 | 2022-08-23 | 水利部交通运输部国家能源局南京水利科学研究院 | 集成时空指标的反演降水产品精度综合评估方法及系统 |
| WO2023240509A1 (zh) * | 2022-06-15 | 2023-12-21 | 中山大学 | 一种基于降水预报与遥相关对应关系的空间概率分析方法及系统 |
| CN116050647B (zh) * | 2023-02-03 | 2023-08-22 | 中国水利水电科学研究院 | 一种基于遥感的流域洪涝受灾强度重现期推求方法 |
| CN117668442B (zh) * | 2023-12-11 | 2024-06-21 | 重庆市气候中心 | 一种无有效气象观测资料拟算点推算重现期降水量的方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103033856A (zh) * | 2012-12-06 | 2013-04-10 | 南京信息工程大学 | 基于水文气象一致性地区分析的雨量频率估计方法 |
| CN104732092A (zh) * | 2015-03-25 | 2015-06-24 | 河海大学 | 一种基于聚类的水文降雨一致区分析方法 |
| KR20170134830A (ko) * | 2016-05-26 | 2017-12-07 | 경상대학교산학협력단 | 혼합분포형을 이용한 지구 기후 모델의 극치 강수자료의 오차 보정방법 |
| CN110059745A (zh) * | 2019-04-17 | 2019-07-26 | 武汉大学 | 一种基于星地融合的流域降雨产品订正方法及系统 |
| CN110555557A (zh) * | 2019-08-30 | 2019-12-10 | 北京师范大学 | 区域性小时级暴雨灾害预警方法、装置与电子设备 |
-
2020
- 2020-06-09 CN CN202010515991.2A patent/CN111797131B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103033856A (zh) * | 2012-12-06 | 2013-04-10 | 南京信息工程大学 | 基于水文气象一致性地区分析的雨量频率估计方法 |
| CN104732092A (zh) * | 2015-03-25 | 2015-06-24 | 河海大学 | 一种基于聚类的水文降雨一致区分析方法 |
| KR20170134830A (ko) * | 2016-05-26 | 2017-12-07 | 경상대학교산학협력단 | 혼합분포형을 이용한 지구 기후 모델의 극치 강수자료의 오차 보정방법 |
| CN110059745A (zh) * | 2019-04-17 | 2019-07-26 | 武汉大学 | 一种基于星地融合的流域降雨产品订正方法及系统 |
| CN110555557A (zh) * | 2019-08-30 | 2019-12-10 | 北京师范大学 | 区域性小时级暴雨灾害预警方法、装置与电子设备 |
Non-Patent Citations (1)
| Title |
|---|
| 基于GPD 分布的黑河流域极端降水频率特征分析;李占玲;王武;李占杰;地理研究;第33卷(第11期);2169-2179 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111797131A (zh) | 2020-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111797131B (zh) | 一种基于遥感降水产品的极端降水区域频率分析方法 | |
| CN105160192B (zh) | 基于M5‑LocalR的TRMM卫星降雨数据降尺度方法 | |
| CN110751094B (zh) | 一种基于gee综合遥感影像和深度学习方法的作物估产方法 | |
| CN110222911B (zh) | 一种卫星遥感与地面数据协同的雨量站网优化布局方法 | |
| WO2020043027A1 (zh) | 一种大气污染预测中对区域边界优化的方法 | |
| CN111401602A (zh) | 基于神经网络的卫星以及地面降水测量值同化方法 | |
| Sridhar et al. | Implications of climate-driven variability and trends for the hydrologic assessment of the Reynolds Creek Experimental Watershed, Idaho | |
| CN105760970A (zh) | 空气质量指数的预测方法 | |
| CN105974495B (zh) | 利用分类拟合法预判目标区域未来平均云量的方法 | |
| Tian et al. | Evaluation of six precipitation products in the Mekong River Basin | |
| CN105550423A (zh) | 基于Fuzzy-OLS的CMORPH卫星降雨数据降尺度方法 | |
| CN103529189A (zh) | 一种基于定性和定量辅助变量的土壤有机质空间分布预测方法 | |
| WO2024152865A1 (zh) | 一种基于气象站的生态系统碳水通量计算方法及系统 | |
| CN108981616B (zh) | 一种由无人机激光雷达反演人工林有效叶面积指数的方法 | |
| CN113108918B (zh) | 一种极轨气象卫星热红外遥感数据反演气温方法 | |
| CN109472283B (zh) | 一种基于多重增量回归树模型的危险天气预测方法和装置 | |
| CN111160680A (zh) | 一种基于信息同化融合的农业干旱评估方法 | |
| CN117763970B (zh) | 基于图像深度学习的稀缺资料地区水文模型参数重建方法 | |
| CN113836808A (zh) | 一种基于重污染特征约束的pm2.5深度学习预测方法 | |
| CN111860974B (zh) | 一种基于状态空间和联合分布的干旱多级预测方法 | |
| CN115203643A (zh) | 融合水文与生态要素的水源涵养功能定量诊断方法及系统 | |
| CN117114176A (zh) | 基于数据分析和机器学习的土地利用变化预测方法及系统 | |
| Amin et al. | Evaluation of the performance of SWAT model to simulate stream flow of Mojo river watershed: in the upper Awash River basin, in Ethiopia | |
| CN109543911A (zh) | 一种日照辐射预测方法及系统 | |
| Scuderi | Quantification of long-term erosion rates from root exposure/tree age relationships in an alpine meadow catchment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |