CN112765852B - 一种用于特小流域的可能最大暴雨计算方法 - Google Patents

一种用于特小流域的可能最大暴雨计算方法 Download PDF

Info

Publication number
CN112765852B
CN112765852B CN202110070383.XA CN202110070383A CN112765852B CN 112765852 B CN112765852 B CN 112765852B CN 202110070383 A CN202110070383 A CN 202110070383A CN 112765852 B CN112765852 B CN 112765852B
Authority
CN
China
Prior art keywords
ultra
rainstorm
small
small watershed
value
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
Application number
CN202110070383.XA
Other languages
English (en)
Other versions
CN112765852A (zh
Inventor
张金良
盖永岗
付健
李超群
王鹏
韦诗涛
蔺冬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yellow River Engineering Consulting Co Ltd
Original Assignee
Yellow River Engineering Consulting Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yellow River Engineering Consulting Co Ltd filed Critical Yellow River Engineering Consulting Co Ltd
Priority to CN202110070383.XA priority Critical patent/CN112765852B/zh
Publication of CN112765852A publication Critical patent/CN112765852A/zh
Application granted granted Critical
Publication of CN112765852B publication Critical patent/CN112765852B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/14Rainfall or precipitation gauges
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90Details of database functions independent of the retrieved data types
    • G06F16/95Retrieval from the web
    • G06F16/953Querying, e.g. by the use of web search engines
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Environmental & Geological Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Hydrology & Water Resources (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Alarm Systems (AREA)

Abstract

本发明提供了一种用于特小流域的可能最大暴雨计算方法,包括:收集设计流域可能最大暴雨等值线图,将位置标定在等值线图中,来计算特小流域的可能最大24h暴雨值;收集典型特大暴雨资料信息;根据露点信息以及可能最大降水估算手册,计算测站相对应Td典型和Td历史的可降水量,对典型特大暴雨实测24h雨量P24h,典测进行放大,重新推求测站经修正的可能最大24h雨量;将特大暴雨发生测站的位置标定在等值线图中,根据测站在等值线图中的位置计算测站在等值线图中的可能最大24h雨量结果,获取上述两者比值作为修正系数;获取特小流域经修正的可能最大24h暴雨值,对获取的经修正的可能最大24h暴雨值进行验证,并将验证结果进行传输显示。便于提高获取的准确性。

Description

一种用于特小流域的可能最大暴雨计算方法
技术领域
本发明涉及水文气象技术领域,特别涉及一种用于特小流域的可能最大暴雨计算方法。
背景技术
目前,可能最大暴雨一般是针对重要的水利水电工程推求可能最大洪水时所需计算的内容,其可能最大暴雨分析计算对象多处于较大流域,分析计算方法已较为成熟。对于核电工程进行防洪安全论证时,需考虑核电厂址周边沟道的可能最大洪水,此时也需计算可能最大暴雨,其分析计算对象一般为小流域或特小流域。目前,计算小流域可能最大暴雨一般包括当地暴雨放大法、暴雨移置法、统计估算法等方法,这些方法的使用要求工程技术人员深入掌握相关专业知识,其中暴雨放大法和暴雨移置法要求具有一定的气象学知识,统计估算法则要求具有一定的统计学知识,对工程技术人员专业知识要求较高,使得可能最大暴雨计算的专业门槛设置过高,不便于广大基层科技工作者掌握和使用。对于一些面广量大的小型工程规划设计,若需计算可能最大暴雨,则需研发便于众多基层科技工作者掌握和使用的简便计算方法。
本发明针对以上情况,以我国在70~80年代组织编制的可能最大暴雨等值线图成果为基础,选用编制等值线图所用资料年份以后发生的典型特大暴雨资料,复核计算典型特大暴雨发生地的可能最大暴雨设计值,与通过可能最大暴雨等值线图计算的值进行比较,以其比值作为对可能最大暴雨等值线图成果的修正系数,进而推求工程点的可能最大暴雨成果。
发明内容
本发明提供一种用于特小流域的可能最大暴雨计算方法,用以通过已知的可能最大暴雨等值线图以及典型特大暴雨资料信息,并以典型特大暴雨资料信息对中国可能最大暴雨等值线图成果进行修正,进而获取工程点的可能最大暴雨结果,提高其获取的准确性。
本发明提供一种用于特小流域的可能最大暴雨计算方法,包括:
步骤1:收集设计流域所在地区的可能最大暴雨等值线图,将设计流域的位置标定在等值线图中,根据设计流域在等值线图中的位置插值计算特小流域的可能最大24h暴雨值,以P24h,小查表示;
步骤2:收集等值线图编制所用资料年份以后设计流域所在区域发生的典型特大暴雨资料信息;
步骤3:根据测站的特大暴雨代表性露点和历史露点信息,分别以Td典型和Td历史表示,根据露点信息以及可能最大降水估算手册,根据Td典型和Td历史的值可查表获取相应的可降水量,分别以Wm和W表示,以K=Wm/W为放大倍比对典型特大暴雨实测24h雨量P24h,典测进行放大,重新推求所述测站经修正的可能最大24h雨量,以P24h,典算表示,即P24h,典算=K*P24h,典测
步骤4:将特大暴雨发生测站的位置标定在等值线图中,根据测站在等值线图中的位置插值计算测站在等值线图中的可能最大24h雨量结果,以P24h,典查表示;
步骤5:以测站的经修正的可能最大24h雨量和在等值线图中计算的可能最大24h雨量成果的比值作为可能最大暴雨等值线图计算结果的修正系数,即修正系数为μ=P24h,典算/P24h,典查
步骤6:获取特小流域经修正的可能最大24h暴雨值P24h,小算=μ*P24h,小查
步骤7:对步骤6获取的经修正的可能最大24h暴雨值进行验证,并将验证结果进行传输显示。
在一种可能实现的方式中,
所述典型特大暴雨资料信息包括:发生特大暴雨的雨量监测站的测站位置与高程、24h实测暴雨量、特大暴雨的代表性露点、雨量监测站的历史最大露点信息。
在一种可能实现的方式中,所述方法还包括:
根据所述等值线图编制所用资料年份以后设计流域所在区域发生的典型特大暴雨资料信息以及设计流域的安全等级确定所述设计流域的洪水发生频率;
根据所述洪水发生频率确定所述设计流域发生洪水的环境参数;
建立特小流域经修正的可能最大24h暴雨值P24h,小算的特征值;
根据所述特小流域经修正的可能最大24h暴雨值P24h,小算的特征值以及设计流域发生洪水的环境参数计算出符合洪水产生要求的目标降雨量;
将所述目标降雨量代入预设双变量联合分布函数中计算出特小区域以特小流域经修正的可能最大24h暴雨值P24h,小算的损失参数;
根据所述等值线图构建特小区域的随机有限元模型;
将所述损失参数作为随机变量输入到所述随机有限元模型中来计算出特小流域经修正的可能最大24h暴雨值P24h,小算在洪水冲击下的随机响应;
根据所述特小流域经修正的可能最大24h暴雨值P24h,小算在洪水冲击下的随机响应创建特小流域环境损坏的极限状态方程;
对所述极限状态方程进行状态模拟,确定特小流域在修正的可能最大24h暴雨值P24h,小算下的多个风险指标的发生概率;
根据每个风险指标的发生概率出计算每个风险指标在所述特小流域环境损坏中所占据的权重值;
根据所述特小流域的地理结构以及修正的可能最大24h暴雨值P24h,小算构建雨水密度分布函数;
根据所述雨水密度分布函数计算出所述特小流域的雨水密度分布情况;
利用每个风险指标所占的权重值和特小流域的雨水密度分布情况计算出所述特小区域的目标损坏度;
将所述目标损坏度与所述特小流域在修正的可能最大24h暴雨值P24h,小算进行关联以进行传输显示。
在一种可能实现的方式中,所述利用每个风险指标所占的权重值和特小流域的雨水密度分布情况计算出所述特小区域的目标损坏度,包括:
根据所述特小流域的雨水密度分布情况和特小流域地理结构的结构参数计算出暴雨在特小流域内的起动流速:
Figure BDA0002905802010000041
其中,V表示为暴雨在特小流域内的起动流速,g表示为重力加速度,m表示为雨滴的平均质量,ρ1表示为特小流域内泥土的密度,ρ2表示为暴雨的密度,S表示为特小流域的面积,h表示为雨滴的粒径,L表示为特小流域内的积水深度,T表示为特小流域内最高点的高度值,T1表示为特小流域内最低点的高度值,t表示为暴雨从开始降落到接触特小流域最高点所需要的时长;
根据所述暴雨在特小流域内的起动流速计算出特小流域的渗透系数:
Figure BDA0002905802010000042
其中,k表示为特小流域的渗透系数,e表示为特小流域的地面孔隙率,Q表示为特小流域的环境安全系数,P表示为单位时间内特小流域的存储的雨水的体积,H表示为雨水的黏滞系数,B表示为特小流域的地表高度,ω表示为特小流域的含沙量影响系数;
根据特小流域的渗透系数和暴雨在特小流域内的起动流速以及每个风险指标所占的权重值计算出特小区域的目标损坏度:
Figure BDA0002905802010000043
其中,M表示为特小区域的目标损坏度,N表示为风险指标的数量,Si表示为第i个风险指标所占的权重值,A表示为特小流域的初始地面结构稳定性指数,λ表示为特小流域内沙粒的流动率;
将所述目标损坏度同预设损坏度进行比较,当所述目标损害度大于等于所述预设损害度时,发出报警提示并显示所述目标损坏度,当所述目标损害度小于所述预设损坏度时,单显示所述目标损坏度。
在一种可能实现的方式中,步骤1中,收集设计流域所在地区的可能最大暴雨等值线图之前,包括:
获取收集与设计流域所在地区的可能最大暴雨等值线图的网址集合;
对所述网址集合中的每个网址进行搜索,筛选界面显示访问成功的第一网址;
同时,对每个所述第一网址进行预处理,根据预处理结果,筛选获得第二网址;
获取所述第二网址对应的显示界面的显示图像,并对每个显示图像进行单独评估,并对两两显示图像进行关联评估,同时,对所述显示界面中的所有图像进行综合评估;
根据单独评估结果、关联评估结果以及综合评估结果,确定所述显示界面的显示内容是否合格,若合格,将所述第二网址进行保留,否则,将所述第二网址进行剔除,获得最终网址。
在一种可能实现的方式中,对每个所述第一网址进行预处理,根据预处理结果,筛选获得第二网址,包括:
记录所述第一网址被访问的次数以及访问过程中,对对应的显示界面的目标图像的获取操作,所述获取操作包括截取、下载操作;
根据所述获取操作,确定操作行为,并建立行为列表,同时,对行为列表中的每个子行为进行分析,获得分析结果;
根据分析结果,判断对对应的显示界面的目标图像的综合获取值;
当所述综合获取值小于预设值时,将所述第一网址进行剔除。
在一种可能实现的方式中,
剔除第一网址后的剩余地址即为第二网址。
在一种可能实现的方式中,
每个合格的第二网址即为最终网址。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1为本发明实施例中一种用于特小流域的可能最大暴雨计算方法的流程图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
本发明提供一种用于特小流域的可能最大暴雨计算方法,如图1所示,包括:
步骤1:收集设计流域所在地区的可能最大暴雨等值线图,将设计流域的位置标定在等值线图中,根据设计流域在等值线图中的位置插值计算特小流域的可能最大24h暴雨值,以P24h,小查表示;
步骤2:收集等值线图编制所用资料年份以后设计流域所在区域发生的典型特大暴雨资料信息;
步骤3:根据测站的特大暴雨代表性露点和历史露点信息,分别以Td典型和Td历史表示,根据露点信息以及可能最大降水估算手册,根据Td典型和Td历史值可查表获取相应的可降水量,分别以Wm和W表示,以K=Wm/W为放大倍比对典型特大暴雨实测24h雨量P24h,典测进行放大,重新推求所述测站经修正的可能最大24h雨量,以P24h,典算表示,即P24h,典算=K*P24h,典测
步骤4:将特大暴雨发生测站的位置标定在等值线图中,根据测站在等值线图中的位置插值计算测站在等值线图中的可能最大24h雨量结果,以P24h,典查表示;
步骤5:以测站的经修正的可能最大24h雨量和在等值线图中计算的可能最大24h雨量成果的比值作为可能最大暴雨等值线图计算结果的修正系数,即修正系数为μ=P24h,典算/P24h,典查
步骤6:获取特小流域经修正的可能最大24h暴雨值P24h,小算=μ*P24h,小查
步骤7:对步骤6获取的经修正的可能最大24h暴雨值进行验证,并将验证结果进行传输显示。
在一种可能实现的方式中,
所述典型特大暴雨资料信息包括:发生特大暴雨的雨量监测站的测站位置与高程、24h实测暴雨量、特大暴雨的代表性露点、雨量监测站的历史最大露点信息。
针对步骤1:收集设计流域所在地区的可能最大暴雨等值线图,将设计流域的位置标定在等值线图中,根据设计流域在等值线图中的位置通过线性插值的方法计算设计流域的可能最大24h暴雨值,以P24h,小查表示。
针对步骤2:可能最大暴雨等值线图编制时所用资料一般截止至70~80年代,所用资料系列长度相对较短,之后的时期内,我国的降雨监测站点数量、监测手段和监测精度不断提高,积累了相当丰富的暴雨资料,应充分运用这一时期的特大暴雨资料,对通过可能最大暴雨等值线图查算的工程点可能最大24h暴雨值进行修正,且应当确保所选用的典型特大暴雨量级是比编制可能最大暴雨等值线图所用的暴雨资料中的最大暴雨还要大。
针对步骤3:进行放大,相当于充分利用了该地区在编制等值线图以后年份的特大暴雨信息来推求该测站经修正的可能最大24h雨量,以P24h,典算表示,即P24h,典算=K*P24h,典测
针对步骤4:由于P24h,典查所依据的暴雨资料年份较短,对大暴雨效率的反应不够充分,因此P24h,典查比要P24h,典算小。
针对步骤5:由于μ=P24h,典算/P24h,典查充分利用了可能最大暴雨等值线图编制所依据的暴雨资料以后发生的典型特大暴雨资料,其特大暴雨效率的反应更加充分,因此,应以μ=P24h,典算/P24h,典查作为该区域可能最大暴雨等值线图查算结果的修正系数;修正系数μ可作为典型测站附近气候一致区的可能最大暴雨等值线图查算结果的修正系数。
上述技术方案的有益效果是:通过已知的可能最大暴雨等值线图以及典型特大暴雨资料信息,并以典型特大暴雨资料信息对中国可能最大暴雨等值线图成果进行修正,进而获取工程点的可能最大暴雨结果,提高其获取的准确性。
本发明提供一种用于特小流域的可能最大暴雨计算方法,所述方法还包括:
根据所述等值线图编制所用资料年份以后设计流域所在区域发生的典型特大暴雨资料信息以及设计流域的安全等级确定所述设计流域的洪水发生频率;
根据所述洪水发生频率确定所述设计流域发生洪水的环境参数;
建立特小流域经修正的可能最大24h暴雨值P24h,小算的特征值;
根据所述特小流域经修正的可能最大24h暴雨值P24h,小算的特征值以及设计流域发生洪水的环境参数计算出符合洪水产生要求的目标降雨量;
将所述目标降雨量代入预设双变量联合分布函数中计算出特小区域以特小流域经修正的可能最大24h暴雨值P24h,小算的损失参数;
根据所述等值线图构建特小区域的随机有限元模型;
将所述损失参数作为随机变量输入到所述随机有限元模型中来计算出特小流域经修正的可能最大24h暴雨值P24h,小算在洪水冲击下的随机响应;
根据所述特小流域经修正的可能最大24h暴雨值P24h,小算在洪水冲击下的随机响应创建特小流域环境损坏的极限状态方程;
对所述极限状态方程进行状态模拟,确定特小流域在修正的可能最大24h暴雨值P24h,小算下的多个风险指标的发生概率;
根据每个风险指标的发生概率出计算每个风险指标在所述特小流域环境损坏中所占据的权重值;
根据所述特小流域的地理结构以及修正的可能最大24h暴雨值P24h,小算构建雨水密度分布函数;
根据所述雨水密度分布函数计算出所述特小流域的雨水密度分布情况;
利用每个风险指标所占的权重值和特小流域的雨水密度分布情况计算出所述特小区域的目标损坏度;
将所述目标损坏度与所述特小流域在修正的可能最大24h暴雨值P24h,小算进行关联以进行传输显示。
上述技术方案的有益效果为:通过计算以特小流域经修正的可能最大24h暴雨值P24h,小算的损失参数可以使得工作人员在了解特小区域的目标降雨量的同时可以同时地知悉在目标降雨量下对特小区域的损坏情况,进而可以提前做好防护措施,提高了安全性,避免出现安全事故的发生,进一步地,通过根据雨水密度分布函数计算出所述特小流域的雨水密度分布情况进而计算出特小区域的目标损坏度可以使用完全数据化的以函数的形式来确定特小流域的损坏度,避免了人工情感等因素的影响,使得计算的数据更加准确和实际,同时,可以使工作人员针对目标损坏度对特小环境内的有价值物品进行特效防护,一定程度上可以避免经济损失,提高了工作人员的使用体验感。
本发明提供一种用于特小流域的可能最大暴雨计算方法,所述利用每个风险指标所占的权重值和特小流域的雨水密度分布情况计算出所述特小区域的目标损坏度,包括:
根据所述特小流域的雨水密度分布情况和特小流域地理结构的结构参数计算出暴雨在特小流域内的起动流速:
Figure BDA0002905802010000101
其中,V表示为暴雨在特小流域内的起动流速,g表示为重力加速度,m表示为雨滴的平均质量,ρ1表示为特小流域内泥土的密度,ρ2表示为暴雨的密度,S表示为特小流域的面积,h表示为雨滴的粒径,L表示为特小流域内的积水深度,T表示为特小流域内最高点的高度值,T1表示为特小流域内最低点的高度值,t表示为暴雨从开始降落到接触特小流域最高点所需要的时长;
根据所述暴雨在特小流域内的起动流速计算出特小流域的渗透系数:
Figure BDA0002905802010000102
其中,k表示为特小流域的渗透系数,e表示为特小流域的地面孔隙率,Q表示为特小流域的环境安全系数,P表示为单位时间内特小流域的存储的雨水的体积,H表示为雨水的黏滞系数,B表示为特小流域的地表高度,ω表示为特小流域的含沙量影响系数;
根据特小流域的渗透系数和暴雨在特小流域内的起动流速以及每个风险指标所占的权重值计算出特小区域的目标损坏度:
Figure BDA0002905802010000103
其中,M表示为特小区域的目标损坏度,N表示为风险指标的数量,Si表示为第i个风险指标所占的权重值,A表示为特小流域的初始地面结构稳定性指数,λ表示为特小流域内沙粒的流动率;
将所述目标损坏度同预设损坏度进行比较,当所述目标损害度大于等于所述预设损害度时,发出报警提示并显示所述目标损坏度,当所述目标损害度小于所述预设损坏度时,单显示所述目标损坏度。
上述技术方案的有益效果为:通过计算暴雨在特小流域内的起动流速可以根据特小流域内的地形复杂情况准确地计算出雨水的流速,无需人工勘测计算,节省了人力,进一步地,通过计算特小流域的渗透系数可以确定在暴雨持续的情况下特小流域内的地面对于雨水的渗透情况,进而根据渗透情况可以准确地评估出暴雨对于特小流域内地表的损害情况,进一步地,通过利用特小流域内的多个环境参数和风险指标参数来计算特小流域的目标损坏度可以在考虑到外在和内在情况的同时来确定暴雨对于特小流域的损坏度,提高了数据的精确性。
在一个实施例中,步骤1中,收集设计流域所在地区的可能最大暴雨等值线图之前,包括:
获取收集与设计流域所在地区的可能最大暴雨等值线图的网址集合;
对所述网址集合中的每个网址进行搜索,筛选界面显示访问成功的第一网址;
同时,对每个所述第一网址进行预处理,根据预处理结果,筛选获得第二网址;
获取所述第二网址对应的显示界面的显示图像,并对每个显示图像进行单独评估,并对两两显示图像进行关联评估,同时,对所述显示界面中的所有图像进行综合评估;
根据单独评估结果、关联评估结果以及综合评估结果,确定所述显示界面的显示内容是否合格,若合格,将所述第二网址进行保留,否则,将所述第二网址进行剔除,获得最终网址。
该实施例中,显示图像与暴雨等值线图有关。
上述技术方案的有益效果是:通过对网址集合进行遍历搜索,便于获取第一网址,且对第一网址进行预处理,便于获取第二网址,通过对显示界面的图像进行单独评估、关联评估以及综合评估,便于确定显示内容是否合格,保证网址的有效性,为获取等值线图提供有效的收集基础。
在一种可能实现的方式中,对每个所述第一网址进行预处理,根据预处理结果,筛选获得第二网址,包括:
记录所述第一网址被访问的次数以及访问过程中,对对应的显示界面的目标图像的获取操作,所述获取操作包括截取、下载操作;
根据所述获取操作,确定操作行为,并建立行为列表,同时,对行为列表中的每个子行为进行分析,获得分析结果;
根据分析结果,判断对对应的显示界面的目标图像的综合获取值;
当所述综合获取值小于预设值时,将所述第一网址进行剔除。
该实施例中,目标图像与暴雨等值线图有关。
上述技术方案的有益效果是:通过对用户的获取操作进行分析,便于有效的分析用户对与该网址的喜爱程度,便于有效的确定,该网址存在的有效信息。
在一种可能实现的方式中,
剔除第一网址后的剩余地址即为第二网址。
在一种可能实现的方式中,
对每个网址进行搜索包括:遍历搜索。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (7)

1.一种用于特小流域的可能最大暴雨计算方法,其特征在于,包括:
步骤1:收集设计流域所在地区的可能最大暴雨等值线图,将设计流域的位置标定在等值线图中,根据设计流域在等值线图中的位置插值计算特小流域的可能最大24h暴雨值,以P24h,小查表示;
步骤2:收集等值线图编制所用资料年份以后设计流域所在区域发生的典型特大暴雨资料信息;
步骤3:根据测站的特大暴雨代表性露点和历史露点信息,分别以Td典型和Td历史表示,根据露点信息以及可能最大降水估算手册,根据Td典型和Td历史值可查表获取相应的可降水量,分别以Wm和W表示,以K=Wm/W为放大倍比对典型特大暴雨实测24h雨量P24h,典测进行放大,重新推求所述测站经修正的可能最大24h雨量,以P24h,典算表示,即P24h,典算=K*P24h,典测
步骤4:将特大暴雨发生测站的位置标定在等值线图中,根据测站在等值线图中的位置插值计算测站在等值线图中的可能最大24h雨量结果,以P24h,典查表示;
步骤5:以测站的经修正的可能最大24h雨量和在等值线图中计算的可能最大24h雨量成果的比值作为可能最大暴雨等值线图计算结果的修正系数,即修正系数为μ=P24h,典算/P24h,典查
步骤6:获取特小流域经修正的可能最大24h暴雨值P24h,小算=μ*P24h,小查
步骤7:对步骤6获取的经修正的可能最大24h暴雨值进行验证,并将验证结果进行传输显示;
所述方法还包括:
根据所述等值线图编制所用资料年份以后设计流域所在区域发生的典型特大暴雨资料信息以及设计流域的安全等级确定所述设计流域的洪水发生频率;
根据所述洪水发生频率确定所述设计流域发生洪水的环境参数;
建立特小流域经修正的可能最大24h暴雨值P24h,小算的特征值;
根据所述特小流域经修正的可能最大24h暴雨值P24h,小算的特征值以及设计流域发生洪水的环境参数计算出符合洪水产生要求的目标降雨量;
将所述目标降雨量代入预设双变量联合分布函数中计算出特小区域以特小流域经修正的可能最大24h暴雨值P24h,小算的损失参数;
根据所述等值线图构建特小区域的随机有限元模型;
将所述损失参数作为随机变量输入到所述随机有限元模型中来计算出特小流域经修正的可能最大24h暴雨值P24h,小算在洪水冲击下的随机响应;
根据所述特小流域经修正的可能最大24h暴雨值P24h,小算在洪水冲击下的随机响应创建特小流域环境损坏的极限状态方程;
对所述极限状态方程进行状态模拟,确定特小流域在修正的可能最大24h暴雨值P24h,小算下的多个风险指标的发生概率;
根据每个风险指标的发生概率计算出每个风险指标在所述特小流域环境损坏中所占据的权重值;
根据所述特小流域的地理结构以及修正的可能最大24h暴雨值P24h,小算构建雨水密度分布函数;
根据所述雨水密度分布函数计算出所述特小流域的雨水密度分布情况;
利用每个风险指标所占的权重值和特小流域的雨水密度分布情况计算出所述特小流域的目标损坏度;
将所述目标损坏度与所述特小流域在修正的可能最大24h暴雨值P24h,小算进行关联以进行传输显示。
2.如权利要求1所述的可能最大暴雨计算方法,其特征在于,
所述典型特大暴雨资料信息包括:发生特大暴雨的雨量监测站的测站位置与高程、24h实测暴雨量、特大暴雨的代表性露点、雨量监测站的历史最大露点信息。
3.如权利要求1所述的可能最大暴雨计算方法,其特征在于,所述利用每个风险指标所占的权重值和特小流域的雨水密度分布情况计算出所述特小区域的目标损坏度,包括:
根据所述特小流域的雨水密度分布情况和特小流域地理结构的结构参数计算出暴雨在特小流域内的起动流速:
Figure FDA0003264136380000031
其中,V表示为暴雨在特小流域内的起动流速,g表示为重力加速度,m表示为雨滴的平均质量,ρ1表示为特小流域内泥土的密度,ρ2表示为暴雨的密度,S表示为特小流域的面积,h表示为雨滴的粒径,L表示为特小流域内的积水深度,T表示为特小流域内最高点的高度值,T1表示为特小流域内最低点的高度值,t表示为暴雨从开始降落到接触特小流域最高点所需要的时长;
根据所述暴雨在特小流域内的起动流速计算出特小流域的渗透系数:
Figure FDA0003264136380000032
其中,k表示为特小流域的渗透系数,e表示为特小流域的地面孔隙率,Q表示为特小流域的环境安全系数,P表示为单位时间内特小流域的存储的雨水的体积,H表示为雨水的黏滞系数,B表示为特小流域的地表高度,ω表示为特小流域的含沙量影响系数;
根据特小流域的渗透系数和暴雨在特小流域内的起动流速以及每个风险指标所占的权重值计算出特小区域的目标损坏度:
Figure FDA0003264136380000033
其中,M表示为特小区域的目标损坏度,N表示为风险指标的数量,Si表示为第i个风险指标所占的权重值,A表示为特小流域的初始地面结构稳定性指数,λ表示为特小流域内沙粒的流动率;
将所述目标损坏度同预设损坏度进行比较,当所述目标损坏 度大于等于所述预设损坏 度时,发出报警提示并显示所述目标损坏度,当所述目标损坏 度小于所述预设损坏度时,单显示所述目标损坏度。
4.如权利要求1所述的可能最大暴雨计算方法,其特征在于,步骤1中,收集设计流域所在地区的可能最大暴雨等值线图之前,包括:
获取收集与设计流域所在地区的可能最大暴雨等值线图的网址集合;
对所述网址集合中的每个网址进行搜索,筛选界面显示访问成功的第一网址;
同时,对每个所述第一网址进行预处理,根据预处理结果,筛选获得第二网址;
获取所述第二网址对应的显示界面的显示图像,并对每个显示图像进行单独评估,并对两两显示图像进行关联评估,同时,对所述显示界面中的所有图像进行综合评估;
根据单独评估结果、关联评估结果以及综合评估结果,确定所述显示界面的显示内容是否合格,若合格,将所述第二网址进行保留,否则,将所述第二网址进行剔除,获得最终网址。
5.如权利要求4所述的可能最大暴雨计算方法,其特征在于,对每个所述第一网址进行预处理,根据预处理结果,筛选获得第二网址,包括:
记录所述第一网址被访问的次数以及访问过程中,对对应的显示界面的目标图像的获取操作,所述获取操作包括截取、下载操作;
根据所述获取操作,确定操作行为,并建立行为列表,同时,对行为列表中的每个子行为进行分析,获得分析结果;
根据分析结果,判断对对应的显示界面的目标图像的综合获取值;
当所述综合获取值小于预设值时,将所述第一网址进行剔除。
6.如权利要求5所述的可能最大暴雨计算方法,其特征在于,
剔除第一网址后的剩余地址即为第二网址。
7.如权利要求4所述的可能最大暴雨计算方法,其特征在于,
对每个网址进行搜索包括遍历搜索。
CN202110070383.XA 2021-01-19 2021-01-19 一种用于特小流域的可能最大暴雨计算方法 Active CN112765852B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110070383.XA CN112765852B (zh) 2021-01-19 2021-01-19 一种用于特小流域的可能最大暴雨计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110070383.XA CN112765852B (zh) 2021-01-19 2021-01-19 一种用于特小流域的可能最大暴雨计算方法

Publications (2)

Publication Number Publication Date
CN112765852A CN112765852A (zh) 2021-05-07
CN112765852B true CN112765852B (zh) 2021-11-02

Family

ID=75703219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110070383.XA Active CN112765852B (zh) 2021-01-19 2021-01-19 一种用于特小流域的可能最大暴雨计算方法

Country Status (1)

Country Link
CN (1) CN112765852B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101094221B1 (ko) * 2011-04-15 2011-12-15 한국수자원공사 수문 및 수리분석 통합시스템
CN106844973A (zh) * 2017-01-25 2017-06-13 云南省水利水电勘测设计研究院 一种推求特小流域设计洪水的方法
CN106951674A (zh) * 2017-02-09 2017-07-14 陕西省水文水资源勘测局 一种无历史资料小流域洪水重现期推求方法
CN106971013A (zh) * 2016-01-13 2017-07-21 辽宁省水文局 山洪灾害小流域设计暴雨洪水计算方法
CN107341586A (zh) * 2017-05-12 2017-11-10 成都理工大学 一种基于降雨的地质灾害发生频率的计算方法
CN111209525A (zh) * 2020-04-24 2020-05-29 中国水利水电科学研究院 一种基于优化推理公式的洪峰流量计算方法及系统
CN111414671A (zh) * 2019-12-20 2020-07-14 福建省水利水电勘测设计研究院 小流域暴雨推求洪水的方法
CN112052602A (zh) * 2020-09-18 2020-12-08 福建省水土保持工作站 一种应用于山区小流域设计暴雨量时程分配计算方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101094221B1 (ko) * 2011-04-15 2011-12-15 한국수자원공사 수문 및 수리분석 통합시스템
CN106971013A (zh) * 2016-01-13 2017-07-21 辽宁省水文局 山洪灾害小流域设计暴雨洪水计算方法
CN106844973A (zh) * 2017-01-25 2017-06-13 云南省水利水电勘测设计研究院 一种推求特小流域设计洪水的方法
CN106951674A (zh) * 2017-02-09 2017-07-14 陕西省水文水资源勘测局 一种无历史资料小流域洪水重现期推求方法
CN107341586A (zh) * 2017-05-12 2017-11-10 成都理工大学 一种基于降雨的地质灾害发生频率的计算方法
CN111414671A (zh) * 2019-12-20 2020-07-14 福建省水利水电勘测设计研究院 小流域暴雨推求洪水的方法
CN111209525A (zh) * 2020-04-24 2020-05-29 中国水利水电科学研究院 一种基于优化推理公式的洪峰流量计算方法及系统
CN112052602A (zh) * 2020-09-18 2020-12-08 福建省水土保持工作站 一种应用于山区小流域设计暴雨量时程分配计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
青衣江某特小流域PMP/PMF的估算;黎恺等;《四川水利》;20140215;第35卷(第01期);第1-3页 *

Also Published As

Publication number Publication date
CN112765852A (zh) 2021-05-07

Similar Documents

Publication Publication Date Title
Kent et al. Evaluation of urban local-scale aerodynamic parameters: implications for the vertical profile of wind speed and for source areas
CN108960599B (zh) 基于反演算法的输电线路暴雨灾害精细化预测方法及系统
CN106991278B (zh) 一种集合降水预报与实时洪水概率预报的耦合方法
CN111651885A (zh) 一种智慧型海绵城市洪涝预报方法
CN110646867A (zh) 一种城市排水监测预警方法和系统
CN111507375B (zh) 一种城市内涝风险快速评估方法及系统
CN114997541B (zh) 一种基于数字孪生技术的城市内涝预测方法及预警平台
Nandalal et al. Event based modeling of a watershed using HEC-HMS
CN111795681A (zh) 一种山洪灾害预警方法、装置、服务器及存储介质
CN111428942B (zh) 基于变网格技术提取微地形因子的线路覆冰厚度预测方法
CN110134907B (zh) 一种降雨缺失数据填补方法、系统及电子设备
CN109472283B (zh) 一种基于多重增量回归树模型的危险天气预测方法和装置
CN112308292A (zh) 一种火灾风险等级分布图的绘制方法
CN112700349B (zh) 测风塔选址方法和装置
CN107807907A (zh) 一种降水分类方法及系统
CN113821895B (zh) 输电线路覆冰厚度预测模型的构建方法、装置及存储介质
CN111080003B (zh) 一种适应复杂山地风电场的测风塔规划选址方法
CN115860268A (zh) 基于cnn-xgboost模型的短时强降水预测方法
CN115100819A (zh) 基于大数据分析的滑坡灾害预警方法、装置及电子设备
CN112613633A (zh) 风电场的气象要素预测方法和装置
Tsakiris et al. Regional drought identification and assessment. Case study in Crete
CN112765852B (zh) 一种用于特小流域的可能最大暴雨计算方法
CN110555538A (zh) 风电场风速预测方法及预测系统
CN111460553A (zh) 一种基于不同纬度评估建筑物安全系数的方法
CN117091675A (zh) 一种大坝水位监测预警方法、系统、电子设备及存储介质

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