摘要：截流系统受到水文水力等随机因素的影响，因此，设计的截流方案的风险如何是需要全面考虑的。本文从截流系统的风险认识入手，确定出能够反映截流系统风险的系统量，在理论分析的基础上，提出采用随机模拟计算方法来计算既考虑水文不确定性又考虑水力不确定性的截流系统风险计算方法，实例验证该方法的可行性。从而可以为截流系统的方案设计与比较起到辅助决策作用。

　　在施工组织设计阶段要考虑截流成功的把握性有多大，就此要对截流的风险有一个明确的认识。立堵截流过程中有两类风险存在，一类是戗堤堤头发生结构坍塌的风险;另一类是截流材料能否抵抗水流的冲刷而保持稳定的风险。前一类主要发生在低流速大水头状况;后一类发生在大流速低水头状况。对于在枯水期进行截流施工的工程，后一类风险是主要的;因此 ，在这里针对此类状况的施工截流系统风险度进行分析和研究。施工截流系统存在风险是因为系统受到很多不确定性因素的影响，归纳这些因素并分类如下：

　　1 风险率模型的建立

　　1.1 风险变量的确定 从施工水力学单个材料稳定的概念来看，立堵风险的本质是要看抛投的材料能否抵抗得住龙口水流的冲刷，如果能则截流基本上能顺利完成;如果不能则系统必然存在风险。抛投材料在龙口的稳定是由材料本身的物理特性和龙口水力状况确定的。在伊兹巴什稳定公式(1)中，可以反映这层关系：

　　(1)当设计确定的抛投材料粒径为ds，此时对应的龙口最大流速为vms;如果实际龙口截流过程中抛投材料的粒径为d，由公式(1)知其所能抵抗的龙口最大流速为vm，若d大于ds，则龙口的最大流速vm大于vms，此时，抛投材料在龙口不能稳定，系统出现了失效情况。这说明龙口流速能够反映风险发生。考虑天然来流量的随机性和分流建筑物泄流量的随机性，龙口泄流量也具有随机性，而龙口流量的变化是引起龙口流速变化的主要原因。从上面分析看出，vm是一个反映了众多不确定因素的综合随机变量;因此，选择它作为风险率模型的具体风险变量。

　　1.2 风险率模型 设功能函数g=vm-vs;式中：vm为考虑了水文、水力不确定因素的截流过程中龙口最大流速，它最具有某种概率分布的随机变量;vs为设计截流龙口最大流速。利用功能函数可以反映以下关系：

　　(2)设系统的风险率为R，则R可用下式计算：

　　(3)式中：f(vm)为vm的概率密度函数。

　　2 对f(vm)的分析

　　截流过程中龙口出现的最大流速vm是由龙口泄流量决定的，龙口泄流量又与上游来流量、分流建筑物分流量、戗堤渗流量、龙口体形等随机量相关。认为在截流过程中某时刻龙口的体形是一个确定值。因此，f(vm)只与来流量随机性、分流量随机性和渗流量随机性因素有关。

　　2.1 来流量水文随机因素 水文随机性主要体现在来流量Q上，对Q的描述，工程界一般采用P-Ⅲ型分布，也有采用对数正态分布和极值Ⅰ型分布等。PⅢ型分布的概率密度函数为：

　　2.2 分流建筑分流量水力随机因素 截流过程中，分流建筑物分流量的随机性产生的本质分流建筑中某些水力参数存在不确定性。施工截流一般安排在河道流量较小的枯水期，分流建筑物中的水流形态一般为明流。不失一般性，采用图2作为过流的计算截面。

　　式中：Qsmax为最大渗透流量，按5%×Qr计算，Qr为截流设计流量;Bos为龙口起始平均宽度;B为截流过程中的龙口平均宽度;Hsmax为龙口合龙但未闭气时的上游水头;Zsmax龙口合龙但未闭气时的绝对落差;H为截流过程中的上游水头;Z为截流过程中的上下游水位差(绝对落差)。在进行截流设计过程中对河道的调蓄作用都简化不予考虑，一方面是使设计偏安全，另一方面是调蓄作用对截流设计的影响不大。河道的调蓄作用主要取决与库容流量关系，在这里，对它的不确定性不予考虑。河道水位流量关系的不确定性也不考虑。

　　2.4 f(vm)的讨论 vm是龙口截流过程中出现的最大流速值，它的获得需要经过截流水力计算，

　　(12)式中：t是反映截流过程的参数;Q为随机来流量;Qdi是泄流建筑物的随机泄流量，它的随机性反映了泄流建筑物中的水力参数的随机性，而在某一截流过程中，它又是随截流时间变化的量;Qst是t时刻戗堤的渗流量;At是截流某时刻的龙口过流面积。

　　分析式(12)，通过解析方法同按P-Ⅲ型分布的f(Q)和按正态分布的f(Qd)获得f(vm)，可以利用“当量正态化”等近似处理方法;但是在处理max{

[1] [2] 下一页