CN118261090A - 中小型水库群对下游交叉建筑物的冲击风险评价方法、系统 - Google Patents

中小型水库群对下游交叉建筑物的冲击风险评价方法、系统 Download PDF

Info

Publication number
CN118261090A
CN118261090A CN202410677418.XA CN202410677418A CN118261090A CN 118261090 A CN118261090 A CN 118261090A CN 202410677418 A CN202410677418 A CN 202410677418A CN 118261090 A CN118261090 A CN 118261090A
Authority
CN
China
Prior art keywords
reservoir
risk
reservoirs
dam
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.)
Pending
Application number
CN202410677418.XA
Other languages
English (en)
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.)
Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
Original Assignee
Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
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 Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources filed Critical Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
Priority to CN202410677418.XA priority Critical patent/CN118261090A/zh
Publication of CN118261090A publication Critical patent/CN118261090A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0635—Risk analysis of enterprise or organisation activities
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/08—Construction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Geometry (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Primary Health Care (AREA)
  • Mathematical Optimization (AREA)
  • General Health & Medical Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Computational Mathematics (AREA)
  • Quality & Reliability (AREA)
  • Civil Engineering (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Architecture (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Structural Engineering (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种中小型水库群对下游交叉建筑物的冲击风险评价方法、系统，针对位于水库下游的重要基础设施可能受到的洪水冲击风险，综合考虑水库规模、工程安全情况及分布模式等多因素的影响，快速有效地对上游流域水库（群）进行风险评估。首先通过上游流域水库基础信息排查明确影响范围内水库数量及分布模式，从溃坝可能性、后果严重性及对下游建筑物冲击程度三个方面进行风险评估。若上游存在多座水库则在风险评估前需要进行风险评价单元划分及库群风险排序，根据水库群不同分布模式选取风险评价方法。本发明从保障下游交叉建筑物等重要基础设施工程安全的角度出发，实现中小型水库（群）冲击风险的快速评估。

Description

中小型水库群对下游交叉建筑物的冲击风险评价方法、系统
技术领域
本发明属于水库大坝工程风险分析领域，具体涉及一种中小型水库（群）对下游交叉建筑物的冲击风险评价方法。
背景技术
截至目前，我国共有各类水库九万余座，其中，中小型水库占比达99%以上。然而，受全球气候变化影响，我国气候情势复杂多变，水旱灾害的突发性、异常性、不确定性突出，多发重发态势加强，严重威胁着水库大坝的运行安全，中小型水库群出险导致下游重大基础设施的风险问题突出。
当前有关水库风险评价的研究工作均是针对单一水库风险或是梯级水库连溃风险展开，尚未考虑流域上水库（群）对某一特定基础设施的冲击风险，尤其对于引调水工程、地铁、高速公路等线性工程，由于其线路较长，沿线周边可能存在众多水库，目前仍缺少一种能够快速评价水库冲击风险的方法。
发明内容
本发明所要解决的技术问题：提出一种中小型水库（群）对下游交叉建筑物的冲击风险评价方法，综合考虑了中小型水库（群）与下游交叉建筑物的相对位置及水库群的分布模式，采用经验公式估算溃坝洪水影响并通过隐患排查确定溃坝概率，为水库对下游重要基础设施的冲击风险评估提供了理论和技术支撑。
本发明为解决以上技术问题而采用以下技术方案：
排查交叉建筑物上游流域，初步确定会对下游交叉建筑物产生影响的流域及流域上中小型水库数量；收集流域及水库相关基础资料，根据同一流域上水库情况确定不同的风险评价方法；其中流域及水库相关基础资料包括水文气象资料、工程特性资料；水文气象资料包括：上游流域暴雨、洪水特征、水库设计洪水及其过程；工程特性资料包括：水库地理位置、工程等级、防洪标准、特征水位与库容、水工建筑物基本情况。
若上游流域仅有一座中小型水库，采用单库冲击风险评价LSD法计算单库溃决对下游交叉建筑物的冲击风险；若上游流域同时存在两座或两座以上中小型水库时，则在单库冲击风险评价LSD法的基础上进一步考虑水库连溃的风险叠加效应，采用库群冲击风险评价方法。
具体技术方案如下：
S1、若上游流域仅有一座中小型水库，采用单库冲击风险评价LSD法计算单库溃决对下游交叉建筑物的冲击风险，计算公式为：
（1）
式中：R代表风险值；L代表溃坝可能性系数；S代表溃坝后果严重性系数；D代表冲击系数，与距离、河段断面、比降、糙率因素相关；
S2、若上游流域同时存在两座或两座以上中小型水库时，则采用库群冲击风险评价方法；包括以下子步骤：
S201、划分并明确水库分布模式：将包括关键交叉建筑物及其所跨流域范围内的所有水库作为一个风险评价单元；根据流域范围内水库相对位置关系，将水库群分为串联库群、并联库群、混联库群三种分布模式；
S202、根据不同风险特点划分控制水库、薄弱水库及一般水库，通过控制水库风险特征值计算单元风险R值；
根据单库冲击风险评价LSD法分别计算单元内各单库风险的L、S值大小；L值最大，即单库溃坝概率最高的为薄弱水库；S值最大即为控制水库，S值较大者在单元内截断风险传导的能力较强；
S203、通过确定不同分布模式下库群整体的溃坝可能性系数L、溃坝后果严重性系数S及冲击系数D，计算风险值R；具体如下：
对于串联库群，L值为考虑上游溃坝风险叠加时控制水库的溃坝概率，S值为溃坝水库后果严重性的叠加值，D值为洪水叠加后库群内临近交叉建筑物的水库对下游的冲击系数；
对于并联库群，首先应计算并联水库同时溃坝时风险值R，其中L值为并联水库同时溃坝的概率，S值为溃坝水库后果严重性的叠加值，D值为并联水库同时溃坝时组合洪水对下游的冲击系数，再分别计算单座水库溃坝风险值，通过比较并联多座水库同时溃坝与单座水库溃坝风险值的大小，取其大值代表库群风险值；
对于混联库群，通过由多个串、并联结构组合而成，其风险值R通过串、并联风险评价方法组合计算而成。
进一步的，本发明还提出一种计算机系统，包括存储器、处理器及存储在存储器上的计算机程序，所述处理器执行所述计算机程序以实现本发明前面提出方法的步骤。
进一步的，本发明还提出一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现本发明前面所提出方法的步骤。
最后，本发明还提出一种计算机程序产品，包括计算机程序或指令，其特征在于，该计算机程序或指令被处理器执行时实现本发明所提出方法的步骤。
本发明采用以上技术方案，与现有技术相比，其显著技术效果如下：
（1）本发明方法充分考虑了水库群分布模式在大坝风险评价中的重要性，通过规定不同的串、并、混联模式对应的溃坝概率估算及溃坝后果计算方法，能够快速识别库群中的控制水库、薄弱水库及一般水库，为快速锚定流域防洪的薄弱环节提供了技术支撑。
（2）本发明综合考虑了中小型水库（群）与下游交叉建筑物的相对位置及水库群的分布模式在流域水库风险评价中的重要性，说明水库上下游之间的水力连接造成了溃坝风险在水库群之间传递的特点，解决了同一流域水库群风险动态复杂难以快速评价的难题。
附图说明
图1是本发明的整体步骤流程图。
图2是本发明的不同水库分布模式的水库群示意图。
图3是本发明实施例的三座水库与下游倒虹吸A工程的地理位置示意图。
具体实施方式
下面将结合说明书附图和具体实施方式对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。
本发明提出的一种中小型水库（群）对下游交叉建筑物的冲击风险评价方法，包括：
排查交叉建筑物上游流域，初步确定可能对下游交叉建筑物产生影响的流域及流域上中小型水库数量。全面收集流域及水库相关基础资料，包括但不限于：上游流域暴雨、洪水特征、水库设计洪水及其过程等水文气象资料，水库地理位置、工程等级、防洪标准、特征水位与库容、水工建筑物基本情况等工程特性资料。
根据同一流域上水库情况确定不同的风险评价方法：
S1、若上游流域仅有一座中小型水库，则考虑单库溃决对下游交叉建筑物的冲击风险，采用单库冲击风险评价LSD法；
S2、若上游流域同时存在两座或两座以上中小型水库时，则在单库冲击风险评价LSD方法的基础上进一步考虑水库连溃的风险叠加效应，采用库群冲击风险评价方法。
本发明所提出的中小型水库（群）对下游交叉建筑物的冲击风险评价方法，所述步骤S1中具体步骤中，单库冲击风险评价LSD法是指综合考虑单座水库溃坝可能性、后果严重性及对下游建筑物冲击程度确定风险值R的大小，可按式（1）进行计算：
（1）
式中：R代表风险值；L代表溃坝可能性系数；S代表溃坝后果严重性系数；D代表冲击系数，与距离、河段断面、比降、糙率因素相关。
对于单座水库来说，L值可通过隐患排查水库病险程度定性确定溃坝可能性并赋值。隐患排查主要包括工程隐患及管理隐患两个方面，按隐患严重程度可分为一般、较重、严重三个等级，具体取值如下：
S值主要根据水库坝高及库容信息确定，基于结构系统可靠度理论的安全性等级划分公式，S值的确定可按式（2）进行计算：
（2）
式中：H表示最大坝高，单位为m；V表示库容，单位为万m3。
为了便于快速评价水库对下游交叉建筑物的冲击风险，采用冲击系数D定性描述溃坝洪水演进至交叉建筑物断面时最大流量Qxm与交叉建筑物校核标准流量Q校核、设计标准流量Q设计之间的相对关系：若Qxm≥Q校核，则D取1；若Q校核＞Qxm≥Q设计，则D取0.8；若Qxm＜Q设计，则D取0.6。
溃坝洪水Qxm计算可采用《水利水电工程设计洪水计算手册》中推荐的公式，具体公式见式（3）。以土石坝逐步溃决过程为例，坝址最大流量Qm采用圣维南公式进行估算，具体公式见式（4）。
（3）
式中：Qxm表示距上游坝址X的溃坝最大流量，单位是m3/s；Qm表示上游坝址溃坝最大流量，有溢洪道的水库需加上溢洪道最大泄洪流量，单位是m3/s；W表示下泄总水量，单位是m3；X表示下游断面距上游坝址距离，单位是m；K表示综合反映河道断面、比降、糙率因素的系数。
坝址最大流量Qm采用圣维南公式进行估算，即：
（4）
式中：b—垮坝长度（m）；B—坝址断面的平均宽度（m）；h0—溃坝时坝前水深（m）；g—重力加速度（m/s2）。
本发明所提出的中小型水库（群）对下游交叉建筑物的冲击风险评价方法，所述步骤S2中具体步骤中，包括以下子步骤：
S201、评价单元划分并明确水库分布模式：将包括关键交叉建筑物及其所跨流域范围内的所有水库作为一个风险评价单元。根据流域内水库相对位置关系，可将水库群分为“串联”“并联”“混联”三种基本类型，示意图见图2。
S202、明确评价单元内控制水库、薄弱水库及一般水库：为保证评价方法易于操作，从偏保守的角度对评价过程适当简化，根据不同风险特点划分控制水库、薄弱水库及一般水库，通过控制水库风险特征值计算单元风险R值。
根据单库冲击风险评价方法（LSD法）分别计算单元内各单库风险的L、S值大小。L值最大，即单库溃坝概率最高的为薄弱水库；S值最大即为控制水库，由于计算S值时综合考虑了坝高、库容两方面因素，一般来说，S值较大者在单元内截断风险传导的能力较强。
S203、根据不同分布模式，分为串联库群风险评价、并联库群风险评价、混联库群风险评价。
本发明所提出的中小型水库（群）对下游交叉建筑物的冲击风险评价方法，所述步骤S203中具体步骤中，串联库群风险评价主要指：
若A、B两座水库经判定为串联（图2中的（a）），单元风险可用式（5）计算。上游水库风险叠加时B水库的溃坝概率可用式（6）进行计算，
（5）
式中：B为控制水库，SA、SB表示A、B水库分别溃坝对下游造成危害的严重程度，D表示洪水叠加后交叉建筑物上游临近水库的溃坝洪水对下游的冲击系数，此处需要将多座水库洪水叠加；表示上游水库风险叠加时B水库的溃坝概率，用下式计算：
（6）
式中：LA、LB分别代表两水库各自溃坝概率，其中LAB表示A水库溃坝情况下B水库发生溃坝概率，此条件概率通过库容差值及上下游出入库流量比较得出：假设无任何预警及应急处置措施，当上库A的总库容＜下库B的防洪库容时，认为下库能够容纳上库溃坝洪水，故LAB=0；当上库A的总库容≥下库B的防洪库容时，若下库最大泄洪流量≥上库溃坝到达时最大组合流量，认为下库能够一定程度减小上库溃坝风险，取LAB=0.5；若下库最大泄洪流量小于上库溃坝到达时最大组合流量，认为上库洪水将不断流入下库最终导致漫坝，故LAB=1。
A、B、C（C为控制水库）三座水库串联时，溃坝概率，以此类推。
本发明所提出的中小型水库（群）对下游交叉建筑物的冲击风险评价方法，所述步骤S203中具体步骤中，并联库群风险评价主要指：
由于并联各水库间并无直接影响，因此A、B两座（图2中的（b））并联水库的风险可按式（7）计算。
（7）
式中：R为该单元风险值；LA、LB分别代表两水库各自溃坝概率；SA、SB分别代表两水库溃坝造成危害的严重程度；D表示上游邻近两水库同时溃坝时组合洪水对下游的冲击系数；
需要注意的是，由于两水库同时溃坝概率远小于单一水库溃坝概率，因此，对于并联库群来说，应当分别计算单座水库溃坝风险值RA、RB，比较R、RA、RB风险值大小后，以最大值代表该单元风险水平。
本发明所提出的中小型水库（群）对下游交叉建筑物的冲击风险评价方法，所述步骤S203中具体步骤中，混联库群风险评价主要指：
混联库群可看作多个串、并联结构组合而成，对于最简单的A、B、C三水库混联分布模式（图2中的（c）），控制水库C在考虑风险叠加效应后可代表该单元整体风险水平，RC可按（8）计算。表示考虑A、B水库溃坝风险叠加时C水库的溃坝概率，可用式（9）计算，其他情况可由该方法进行类推。
（8）
式中：SA、SB、SC表示A、B、C水库溃坝造成危害的严重程度，D表示洪水叠加后交叉建筑物上游临近水库的溃坝洪水对下游的冲击系数；表示考虑A、B水库溃坝风险叠加时C水库的溃坝概率，用下式计算：
（9）
式中：LA、LB、LC分别代表三座水库各自溃坝概率，LAC、LBC、LAB,C分别代表A、B水库分别溃坝及同时溃坝时下游C水库的溃坝概率。
本发明的一个实施例是对河北邯郸三座水库及下游一倒虹吸工程进行风险评价分析。
ZZ水库为小（2）型，均质土坝，总库容87.14万m3，最大坝高10.9m，洪水标准为10年一遇洪水设计，50年一遇洪水校核。SQ水库为小（2）型水库，总库容38.42万m3，大坝为均质土坝，最大坝高10.11m，坝顶长136m。水库防洪标准为10年一遇洪水设计，50年一遇洪水校核，最大下泄流量42.8 m³/s。LJH水库与SQ水库并联，为小（2）型水库，总库容68.7万m3，大坝为均质土坝，最大坝高12.6m，洪水标准为20年一遇洪水设计，200年一遇洪水校核。
ZZ水库下游2.19km处为一排水倒虹吸A工程，洪水标准为50年一遇洪水设计，设计流量为319.6m3/s，200年一遇洪水校核，校核流量为490.9m3/s。三座水库与倒虹吸工程位置示意图见图3，三座水库的工程基础信息如下表所示，其中溃前水深为假定值，溃口平均宽度根据（针对库容小于100万m3，K2根据坝体施工和管理质量确定，此处均取9.1）。
由于缺少现场隐患排查的信息，此处假定三座水库的工程隐患及管理隐患排查结果均为一般，具体赋值如下：
根据溃坝洪水估算结果，通过对比与下游倒虹吸工程设计及校核流量可确定D值的大小，由于三座水库溃坝洪水演进至倒虹吸A工程处最大流量均超过该工程校核流量490.9m3/s，因此D值均取为1.0。三座水库分别溃坝对下游倒虹吸工程的风险R值及L、S、D值如下表，其中ZZ为控制水库，由于L值均假定为0.2，不存在相对薄弱水库。由单库冲击风险评价结果可以看出，ZZ水库对下游的冲击风险相对较大，其次是SQ水库，最后是ZZ水库。
进一步考虑三座水库之间的关联对下游倒虹吸A工程的冲击风险，采用式（5）~（9）计算可得出以下结论：
1）假设SQ和LJH水库串联溃坝，其风险值R为0.17；
2）假设LJH水库溃坝导致下游ZZ水库连溃，其风险值R为1.84；
3）假设SQ水库溃坝导致下游ZZ水库连溃，其风险值R为1.29；
4）假设LJH和SQ水库同时溃坝导致下游ZZ水库溃坝，其风险值为1.584。
结合以上单库溃坝及水库群溃坝风险评价结果可知，LJH水库溃坝导致ZZ水库连溃情况下对下游倒虹吸A工程的冲击风险最大，应当予以重视。
需要说明的是，本申请实施例系统的描述，与上述方法实施例的描述是类似的，具有同方法实施例相似的有益效果，因此不做赘述。
作为本发明的另一个实施例，本发明还提出一种计算机系统，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现本发明方法的步骤。
最后，本发明提出一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现前述本发明方法的步骤。
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、场可编程门阵列（FPGA）、专用集成电路（ASIC）、专用标准产品（ASSP）、芯片上系统的系统（SOC）、负载可编程逻辑设备（CPLD）、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
用于实施本申请的方法的程序代码可以采用一个或多个编程语言的任何组合来编写。这些程序代码可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器或控制器，使得程序代码当由处理器或控制器执行时使流程图和/或框图中所规定的功能/操作被实施。程序代码可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本申请的上下文中，机器可读介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的程序。机器可读介质可以是机器可读信号介质或机器可读储存介质。机器可读介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器（RAM）、只读存储器（ROM）、可擦除可编程只读存储器（EPROM或快闪存储器）、光纤、便捷式紧凑盘只读存储器（CD-ROM）、光学储存设备、磁储存设备、或上述内容的任何合适组合。
以上所述仅是本发明的部分实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种中小型水库群对下游交叉建筑物的冲击风险评价方法，其特征在于，包括：
S1、若上游流域仅有一座中小型水库，采用单库冲击风险评价LSD法计算单库溃决对下游交叉建筑物的冲击风险，计算公式为：
（1）
式中：R代表风险值；L代表溃坝可能性系数；S代表溃坝后果严重性系数；D代表冲击系数，与距离、河段断面、比降、糙率因素相关；
S2、若上游流域同时存在两座或两座以上中小型水库时，则采用库群冲击风险评价方法；包括以下子步骤：
S201、划分并明确水库分布模式：将包括关键交叉建筑物及其所跨流域范围内的所有水库作为一个风险评价单元；根据流域范围内水库相对位置关系，将水库群分为串联库群、并联库群、混联库群三种分布模式；
S202、根据不同风险特点划分控制水库、薄弱水库及一般水库，通过控制水库风险特征值计算单元风险R值；
根据单库冲击风险评价LSD法分别计算单元内各单库风险的L、S值大小；L值最大，即单库溃坝概率最高的为薄弱水库；S值最大即为控制水库，S值较大者在单元内截断风险传导的能力较强；
S203、通过确定不同分布模式下库群整体的溃坝可能性系数L、溃坝后果严重性系数S及冲击系数D，计算风险值R；具体如下：
对于串联库群，L值为考虑上游溃坝风险叠加时控制水库的溃坝概率，S值为溃坝水库后果严重性的叠加值，D值为洪水叠加后库群内临近交叉建筑物的水库对下游的冲击系数；
对于并联库群，首先应计算并联水库同时溃坝时风险值R，其中L值为并联水库同时溃坝的概率，S值为溃坝水库后果严重性的叠加值，D值为并联水库同时溃坝时组合洪水对下游的冲击系数，再分别计算单座水库溃坝风险值，通过比较并联多座水库同时溃坝与单座水库溃坝风险值的大小，取其大值代表库群风险值；
对于混联库群，通过由多个串、并联结构组合而成，其风险值R通过串、并联风险评价方法组合计算而成。
2.根据权利要求1所述的方法，其特征在于，对于单座水库，L值是通过隐患排查水库病险程度定性确定溃坝可能性并赋值；隐患排查包括工程隐患及管理隐患，按隐患严重程度分为三级，由低到高分别赋值0.1、0.3、0.5。
3.根据权利要求1所述的方法，其特征在于，S值是根据水库坝高及库容信息确定，基于结构系统可靠度理论的安全性等级划分公式，S值根据下式确定：
（2）
式中：H表示最大坝高，单位为m；V表示库容，单位为万m3。
4.根据权利要求1所述的方法，其特征在于，采用冲击系数D定性描述溃坝洪水演进至交叉建筑物断面时最大流量Qxm与交叉建筑物校核标准流量Q校核、设计标准流量Q设计之间的相对关系：若Qxm≥Q校核，则D取1；若Q校核＞Qxm≥Q设计，则D取0.8；若Qxm＜Q设计，则D取0.6；
（3）
式中：Qxm表示距上游坝址X的溃坝最大流量，单位是m3/s；Qm表示上游坝址溃坝最大流量，有溢洪道的水库需加上溢洪道最大泄洪流量，单位是m3/s；W表示下泄总水量，单位是m3；X表示下游断面距上游坝址距离，单位是m；K表示综合反映河道断面、比降、糙率因素的系数。
5.根据权利要求1所述的方法，其特征在于，串联库群风险评价指：
若A、B两座水库经判定为串联，单元风险用下式计算：
（5）
式中：B为控制水库，SA、SB表示A、B水库分别溃坝对下游造成危害的严重程度，D表示洪水叠加后交叉建筑物上游临近水库的溃坝洪水对下游的冲击系数，此处需要将多座水库洪水叠加；表示上游水库风险叠加时B水库的溃坝概率，用下式计算：
（6）
式中：LA、LB分别代表两水库各自溃坝概率，LAB表示A水库溃坝情况下B水库发生溃坝概率，此条件概率通过库容差值及上下游出入库流量比较得出：
假设无任何预警及应急处置措施，当上库A的总库容＜下库B的防洪库容时，认为下库能够容纳上库溃坝洪水，故LAB=0；当上库A的总库容≥下库B的防洪库容时，若下库最大泄洪流量≥上库溃坝到达时最大组合流量，认为下库能够一定程度减小上库溃坝风险，取LAB=0.5；若下库最大泄洪流量小于上库溃坝到达时最大组合流量，认为上库洪水将不断流入下库最终导致漫坝，故LAB=1；
设C为控制水库，A、B、C三座水库串联时，溃坝概率，以此类推。
6.根据权利要求1所述的方法，其特征在于，对于并联库群风险评价，首先计算A、B两座并联水库同时溃坝时的风险值，按下式计算：
（7）
式中：R为该单元风险值；LA、LB分别代表两水库各自溃坝概率；SA、SB分别代表两水库溃坝造成危害的严重程度；D表示上游邻近两水库同时溃坝时组合洪水对下游的冲击系数；
其次，分别计算单座水库溃坝风险值RA、RB，比较R、RA、RB风险值大小后，以最大值代表该单元风险水平。
7.根据权利要求1所述的方法，其特征在于，混联库群风险评价指：
将混联库群看作多个串、并联结构组合而成，对于A、B、C三水库混联分布模式，控制水库C在考虑风险叠加效应后代表该单元整体风险水平，RC按下式计算：
（8）
式中：SA、SB、SC表示A、B、C水库溃坝造成危害的严重程度，D表示洪水叠加后交叉建筑物上游临近水库的溃坝洪水对下游的冲击系数；表示考虑A、B水库溃坝风险叠加时C水库的溃坝概率，用下式计算：
（9）
式中：LA、LB、LC分别代表三座水库各自溃坝概率，LAC、LBC、LAB,C分别代表A、B水库分别溃坝及同时溃坝时下游C水库的溃坝概率。
8.一种计算机系统，包括存储器、处理器及存储在存储器上的计算机程序，其特征在于，所述处理器执行所述计算机程序以实现权利要求1-7任一所述方法的步骤。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1-7任一所述方法的步骤。
10.一种计算机程序产品，包括计算机程序或指令，其特征在于，该计算机程序或指令被处理器执行时实现权利要求1-7任一所述方法的步骤。
CN202410677418.XA 2024-05-29 2024-05-29 中小型水库群对下游交叉建筑物的冲击风险评价方法、系统 Pending CN118261090A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410677418.XA CN118261090A (zh) 2024-05-29 2024-05-29 中小型水库群对下游交叉建筑物的冲击风险评价方法、系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410677418.XA CN118261090A (zh) 2024-05-29 2024-05-29 中小型水库群对下游交叉建筑物的冲击风险评价方法、系统

Publications (1)

Publication Number Publication Date
CN118261090A true CN118261090A (zh) 2024-06-28

Family

ID=91613569

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410677418.XA Pending CN118261090A (zh) 2024-05-29 2024-05-29 中小型水库群对下游交叉建筑物的冲击风险评价方法、系统

Country Status (1)

Country Link
CN (1) CN118261090A (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101341437B1 (ko) * 2012-12-27 2013-12-13 성균관대학교산학협력단 다기준 홍수위험도 평가방법
CN107044105A (zh) * 2017-05-27 2017-08-15 中国三峡建设管理有限公司 模拟上游大坝溃坝对下游大坝产生冲击影响的实验装置
CN111260213A (zh) * 2020-01-15 2020-06-09 扬州大学 针对梯级水库群的多风险源的水资源风险评估方法
CN113971362A (zh) * 2021-10-26 2022-01-25 中国水利水电科学研究院 基于能量准则的混联水库群防洪优化调度方案生成方法
CN114565211A (zh) * 2021-12-29 2022-05-31 郑州大学 一种风险传递及叠加作用下梯级水库群溃坝风险后果评估方法
CN117610272A (zh) * 2023-11-27 2024-02-27 北京工业大学 一种梯级水电枢纽群大坝连溃风险计算方法和系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101341437B1 (ko) * 2012-12-27 2013-12-13 성균관대학교산학협력단 다기준 홍수위험도 평가방법
CN107044105A (zh) * 2017-05-27 2017-08-15 中国三峡建设管理有限公司 模拟上游大坝溃坝对下游大坝产生冲击影响的实验装置
CN111260213A (zh) * 2020-01-15 2020-06-09 扬州大学 针对梯级水库群的多风险源的水资源风险评估方法
CN113971362A (zh) * 2021-10-26 2022-01-25 中国水利水电科学研究院 基于能量准则的混联水库群防洪优化调度方案生成方法
CN114565211A (zh) * 2021-12-29 2022-05-31 郑州大学 一种风险传递及叠加作用下梯级水库群溃坝风险后果评估方法
CN117610272A (zh) * 2023-11-27 2024-02-27 北京工业大学 一种梯级水电枢纽群大坝连溃风险计算方法和系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
FANG WANG 等: "A Bayesian method for dam failure risk analysis using causal loop diagrams and word frequency analysis", NATURAL HAZARDS, 23 September 2023 (2023-09-23) *
李同春 等: "基于改进SPH模型的溃坝洪水演进模拟方法", 水科学进, 10 October 2023 (2023-10-10) *
胡亚: "长江支流流域梯级水库溃坝风险分析", 河南水利与南水北调, no. 08, 30 August 2018 (2018-08-30) *
谢玉葳 等: "中小流域梯级水库群风险评估及防范对策研究", 人民黄河, vol. 45, no. 8, 31 August 2023 (2023-08-31) *

Similar Documents

Publication Publication Date Title
CN104898183B (zh) 暴雨城市内涝建模评估方法
CN112749475B (zh) 一种确定梯级水库群连续溃坝风险分析方法
Barbaro Master Plan of solutions to mitigate the risk of coastal erosion in Calabria (Italy), a case study
CN104111091A (zh) 一种泥石流力学参数监测系统及泥石流预警系统
Osti et al. Method to improve the mitigative effectiveness of a series of check dams against debris flows
CN104915891A (zh) 一种核电厂地震叠加外部水淹风险评估方法
CN115203993A (zh) 一种高烈度区内河架空直立式框架码头抗震韧性评估方法、系统、设备以及介质
CN115861012B (zh) 一种多源数据融合预警发布方法、装置、电子设备及介质
CN104535346B (zh) 一种混凝土重力危坝的抗滑稳定性检测方法
CN113743658A (zh) 基于临界雨量的中小流域地质灾害及洪水预警方法
CN104482991A (zh) 一种确定危坝坝体安全库水位的测定参数与预警方法
CN106504158A (zh) 一种山丘区沿河村落山洪灾害成灾水位及危险区划分方法
CN114565211A (zh) 一种风险传递及叠加作用下梯级水库群溃坝风险后果评估方法
CN117610272A (zh) 一种梯级水电枢纽群大坝连溃风险计算方法和系统
CN118261090A (zh) 中小型水库群对下游交叉建筑物的冲击风险评价方法、系统
Sato et al. Study on a simplified method of tsunami risk assessment
CN113221440A (zh) 一种排水系统监测点优化布置与实时全局反演方法
CN104596735B (zh) 一种桥梁群的优化布置方法
KR101272854B1 (ko) 단위 지진 해일에 의한 지진해일 데이터베이스 구축 및 지진해일 계산 방법
Waskito et al. Simulation of Saguling dam break using the HEC-RAS software
CN103952998B (zh) 水库下游设计标准下泄流量的测算方法
Lai et al. A numerical modeling study of sediment bypass tunnels at shihmen reservoir, Taiwan
Chongxun et al. Risk analysis for earth dam overtopping
KR101272855B1 (ko) 단위 지진 해일에 의한 지진해일 데이터베이스 구축 및 지진해일 계산 시스템
Sakamoto Case studies in loss of life risk assessment for dam facilities using a simulation model

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