基于最大龙口流速指标的立堵截流施工风险随机模拟研.docVIP

　　基于最大龙口流速指标的立堵截流施工风险随机模拟研 摘要：截流系统受到水文水力等随机因素的影响，因此，设计的截流方案的风险如何是需要全面考虑的。本文从截流系统的风险认识入手，确定出能够反映截流系统风险的系统量，在理论分析的基础上，提出采用随机模拟计算方法来计算既考虑水文不确定性又考虑水力不确定性的截流系统风险计算方法，实例验证该方法的可行性。从而可以为截流系统的方案设计与比较起到辅助决策作用。 关键词：截流系统 风险分析 随机模拟 在施工组织设计阶段要考虑截流成功的把握性有多大，就此要对截流的风险有一个明确的认识。立堵截流过程中有两类风险存在，一类是戗堤堤头发生结构坍塌的风险；另一类是截流材料能否抵抗水流的冲刷而保持稳定的风险。前一类主要发生在低流速大水头状况；后一类发生在大流速低水头状况。对于在枯水期进行截流施工的工程，后一类风险是主要的；因此 ，在这里针对此类状况的施工截流系统风险度进行分析和研究。施工截流系统存在风险是因为系统受到很多不确定性因素的影响，归纳这些因素并分类如下： 1 风险率模型的建立 1.1 风险变量的确定 从施工水力学单个材料稳定的概念来看，立堵风险的本质是要看抛投的材料能否抵抗得住龙口水流的冲刷，如果能则截流基本上能顺利完成；如果不能则系统必然存在风险。抛投材料在龙口的稳定是由材料本身的物理特性和龙口水力状况确定的。在伊兹巴什稳定公式(1)中，可以反映这层关系： (1) 当设计确定的抛投材料粒径为ds，此时对应的龙口最大流速为vms；如果实际龙口截流过程中抛投材料的粒径为d，由公式(1)知其所能抵抗的龙口最大流速为vm，若d大于ds，则龙口的最大流速vm大于vms，此时，抛投材料在龙口不能稳定，系统出现了失效情况。这说明龙口流速能够反映风险发生。考虑天然来流量的随机性和分流建筑物泄流量的随机性，龙口泄流量也具有随机性，而龙口流量的变化是引起龙口流速变化的主要原因。从上面分析看出，vm是一个反映了众多不确定因素的综合随机变量；因此，选择它作为风险率模型的具体风险变量。 1.2 风险率模型 设功能函数g=vm-vs；式中：vm为考虑了水文、水力不确定 6 7 8 9 因素的截流过程中龙口最大流速，它最具有某种概率分布的随机变量；vs为设计截流龙口最大流速。利用功能函数可以反映以下关系： (2) 设系统的风险率为R，则R可用下式计算： (3) 式中：f(vm)为vm的概率密度函数。 2 对f(vm)的分析 截流过程中龙口出现的最大流速vm是由龙口泄流量决定的，龙口泄流量又与上游来流量、分流建筑物分流量、戗堤渗流量、龙口体形等随机量相关。认为在截流过程中某时刻龙口的体形是一个确定值。因此，f(vm)只与来流量随机性、分流量随机性和渗流量随机性因素有关。 2.1 来流量水文随机因素 水文随机性主要体现在来流量Q上，对Q的描述，工程界一般采用P-Ⅲ型分布，也有采用对数正态分布和极值Ⅰ型分布等。PⅢ型分布的概率密度函数为： (4) 其分布函数为 (5) 上两式中：b=mQ(1-2Cv／Cs)；α＝4／C2s；β=2／(mQCsCv)，其中：mQ为年洪峰流量系列均值，;Cv变差系数，Cs偏态系数。 2.2 分流建筑分流量水力随机因素 截流过程中，分流建筑物分流量的随机性产生的本质分流建筑中某些水力参数存在不确定性。施工截流一般安排在河道流量较小的枯水期，分流建筑物中的水流形态一般为明流。不失一般性，采用图2作为过流的计算截面，分流量计算公式为： (6) 式中：A为过水断面面积；χ为湿周；S为底坡；n为糙率。 对于上面各水力因子xi(参数S、n、b、m)，可以认为是服从三角形分布的随机变量，相应的有： (7) 式中：ai,bi,ci分别为xi的最小值，最可能值和最大值；各水力因素xi的平均值用xi表示。在Qd处按泰勒级数展开Qd,分别取均值和方差，同时根据许瓦兹不等式，略去高阶小量，可以得到： (8) (9) 有关研究表明，概率分布的形式对保证率的影响并不敏感。因此，为简化计算，假设分流建筑物的泄流能力服从位置参数μQd=Qd和尺度参数σQd=CQdQd的正态分布。其概率密度函数为： (10) 2.3 其它随机因素 其它不确定因素还有渗流量不确定性、库容水位流量关系、河道水位流量关系等不确定性因素。对于渗流量的确定有模型试验方法和水力计算方法等，由于影响渗流量的因素比较多，且难于考虑；同时，许多的工程实测资料显示，截流过程中最大渗流量占到截流流量的3%～10%内，随机因素对渗流量产生的变化对总的截流过程中龙口水力参数影响不太大，因此，对渗流量的不确定性暂不考虑，戗堤渗透流量参照《水利水电施工组织设计与系统分析》介绍的方法进行计算。公式如下： (11) 式中：Qsmax为最大渗透流量