[发明专利]一种POT洪水自动提取方法及系统

权利要求书

1.一种POT洪水自动提取方法，其特征在于，所述POT洪水自动提取方法包括：

1）准备好日观测流量序列文件；

2）给定欲提取POT样本的起始年Y_begin、终止年Y_end，年均POT洪水发生次数、流域汇流时间或两个洪水峰间隔时间T_interval；

3）根据给定欲提取POT样本的起始年、终止年，年均POT洪水发生次数，计算最终要得到的POT洪水次数， N_POTs= Y_POT*( Y_end- Y_begin+1)；

4）计算起始年、终止年逐日流量序列的最大值；

5）随机给定一个大流量值作为初始阈值 Threshold0，要求该值大于流量序列中的最大值Flow-max；

6）逐日判断流量序列内日流量是否大于初始阈值Threshold0且大于前后2天流量值，计算满足条件的流量出现次数N次，同时记录满足条件的流量值；

7）如果满足条件的流量出现次数N次大于最终要得到的POT洪水次数，则执行8）；否则，以Threshold0＝Threshold0 *0.1作为新的阈值返回6）重新计算；

8）对6）挑选出的流量值进行升值排序，最大流量值作为第1个洪水样本，洪水样本数为1；从第2最大值开始，按流量从大到小依次计算各流量与已选洪水样本的最小时间间隔，所述最小时间间隔大于给定的T_interval，则将流量值作为新的洪水样本；然后从这个流量开始继续判断剩余流量与前面已选洪水样本的最小间隔时间，如果大于给定间隔时间T_interval，则在洪水样本中加入该流量，洪水样本加1；如此循环一直到最小流量，记录下所有入选的洪水样本；

9）如果全部入选洪水样本数小于最终要得到的洪水样本N_POTs，计算Threshold0=Threshold0*0.1返回6）；如果全部入选洪水样本数大于等于最终要得到的洪水次数N_POTs，执行10）；

10）输出所挑选洪水样本数，如果该所挑选洪水样本数远大于最终要得到的洪水样本数，返回5）重新给一个新的初始阈值Threshold0依次执行各步，直到所挑选洪水样本数等于或略大于N_POTs为止；

11）输出所挑选POT样本前N_POTs个，最后一次所用阈值作为最终阈值。

2.一种信息数据处理终端，其特征在于，用于实现权利要求1所述POT洪水自动提取方法。

3.一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如权利要求1所述的POT洪水自动提取方法。

4.一种如权利要求1所述的POT洪水自动提取方法的POT洪水自动提取系统，其特征在于，所述POT洪水自动提取系统包括：

日观测流量序列准备模块，用于准备好日观测流量序列文件，文件内容包含年、月、日、观测流量；

阈值设定模块，用于给定欲提取POT样本的起始年、终止年，年均POT洪水发生次数、流域汇流时间或两个洪水峰峰间隔时间；

自动试错计算模块，用于运行程序通过试错方法计算不断调整阈值，确定最终阈值使POT洪水满足年均POT洪水发生次数、峰间间隔两个条件；

结果输出模块，用于根据确定阈值及年均POT洪水发生次数，筛选出起止年内POT洪水样本及各样本对应日期。