施工导流联合泄流计算模型及导流方案风险分析

施工导流联合泄流计算模型及导流方案风险分析甄红锋1柴换成2（1.2.山东水利职业学院山东日照276826）摘要：通过分别研究明渠导流、隧洞导流以及坝体缺口导流的水力计算，提出切实可行的联合泄流水力计算模型。在施工导流系统联合泄流能力的计算基础上，采用Monte-Carlo方法，模拟施工期洪水入库过程和导流建筑物的泄流工况，用系统仿真方法进行施工洪水调洪演算得到堰前水位的随机样本数据，通过聚类分析确定导流系统动态风险。实例验证说明，施工导流联合泄流水力计算以及风险分析计算模型是可靠的、适用的。关键词：施工导流；联合泄流；风险分析；Monte-Carlo方法；聚类分析分类号:TV6ComputationalModelforJointDischargeandRiskAnalysisofConstructionDiversionSchemesHongfengZhen1，HuanchengChai（1.2.ShandongWaterPolytechnic,Rizhao,276826）Abstract:Throughtheseparateanalysisofchanneldiversion、tunneldiversionanddamgapdischarge,thispaperestablishesthecomputationalmodelforjointdischargeduringconstructionbyapplyingthehydrauliccomputationtheory.Onthejointdischargebasisofconstructiondiversionsystem,usingtheMonte-carlomethodtosimulatethefloodanddischargeprocess,theriskmodelanalysesthedistributionofthefloodlevelinfrontofthecofferdamanddynamicriskunderdifferentdiversionstandardbyfloodroutingmodelandtheclusteranalysismodel.Thecasestudyinstitutesthatthecomputationalmodelforjointdischargeandriskanalysismodelofconstructiondiversionschemesarereasonableandeffective.Keywords:constructiondiversion;jointdischarge;riskanalysis;Monte-Carlomethod;clusteranalysis1.引言施工导流工程是风险工程，在进行导流方案确定、导流方式选择、设计优化时，工程的风险应作为一个指标参与决策，所以如何确定施工导流工程的综合风险就显得十分重要。目前，施工导流工程中所存在的风险问题己经引起广泛的重视，并将风险与导流标淮[1]的选择联系起来进行研究。而进行施工导流工程风险分析计算的前提是首先针对导流工程不同导流时段其不同导流方式和导流建筑物的组合进行联合泄流水力计算[2]。本文阐述了一种广泛适用的联合泄流计算原理，提出了切实可行的联合泄流水力计算算法，提高水力计算的速度和精度。在联合泄流水力计算的基础上，采用Monte-Carlo方法对土石围堰挡水导流风险进行分析，利用Monte-Carlo方法模拟施工洪水过程和导流建筑物泄流工况，通过系统仿真进行施工洪水调洪演算，对得到的上游堰前水位分布进行聚类分析，从而确定在不同的导流标准条件下围堰运行的动态风险。2.施工导流联合泄流水力计算模型1国内外水利水电导流工程实践表明，无论是土石坝枢纽，还是混凝土坝枢纽，施工导流过程中，一般都是几个导流建筑物联合泄流。混凝土坝施工后期的导流渡汛，往往采用坝体预留缺口和底孔泄流方式；在峡谷地区进行水电工程施工导流时，常需要采用隧洞群的导流方式。导流洞还会遇到多层、多条群洞联合泄流，同一洞不同时期或不同洞同一时期可能同时存在有压、半有压、无压的复杂局面。2.1明渠导流水力计算模型明渠泄流时，水流可能呈均匀流或非均匀流[34]。明渠均匀流的水流保持匀速直线运动，可采用谢才公作者简介：甄红锋（1973－），男，硕士在读，山东巨野人，从事水利工程施工与水力学的教学研究工作。柴换成（1980－），男，学士，甘肃天水人，从事水利工程施工与水利工程管理的教学研究工作。式RiACQ计算。人工渠道的水工流绝大部分是明渠恒定非均匀流。对于渐变流水面线的计算采用分段求和法。2.2隧洞导流水力计算模型隧洞泄水能力曲线时是按已知流量Q推求上游水位H。（1）隧洞无压流时上游水位：UPsZmBQH324272.4（1）式中，s为淹没系数，当CHH25.1时取1.0，当CHH25.1时，由HHC确定，CH为隧洞进口内的收缩水深；B为临界水深以下断面平均宽度，kkHAB；流量系数m由进口体形确定。（2）隧洞半有压流时上游水位计算半有压短洞上游水位：UPDDZHAQH205102.0（2）半有压长洞上游水位：DOWNDDZHAQH205102.0（3）（3）隧洞有压流时上游水位计算有压淹没流：DOWNTDZHAQH205102.0（4）有压自由出流：DOWNDDZHAQH205102.0（5）2.3缺口泄流水力计算模型缺口泄流，当堰顶长度l与水头H的关系在区间(HH20,5.2)时，按如下宽顶堰公式计算：2/302HgmBQ（6）式中，B为过水宽度；0H为缺口底槛以上上游水头；为侧收缩系数；m为流量系数；为淹没系数，当自由出流时，取1.0，当淹没出流时，可查宽顶堰淹没系数表得到。2.4联合泄流水力计算模型[5]无论是隧洞群联合导流还是缺口与底孔联合导流，在已知总泄流量和下游水位的情况下，任何泄水建筑物联合泄流必须同时满足以下两个条件：nniiZZZZniQQ211),,2,1(（7）式中，iQ为同一上游水位下，一个泄流建筑物的泄流量；Z为上游水位。N个泄水建筑物联合导流泄流能力计算步骤：(1)根据原河床坝址处天然水位～流量关系，先拟定a个下游水位downH，用插值法找到相对应的a个流量),2,1(aiQi。(2)求出各建筑物对应每一级流量iQ的上游水位的最小值2Z，并找出当0iQ时对应各建筑物的上游水位的最小值1Z，然后在区间(21,ZZ)上取初值upZ，用单条隧洞水力计算公式分别计算对应的泄流量)N21(,,jQij。(3)将1iQ、2iQ…iNQ相加得到NjijQ1，如果EPSQQiNjij1(误差允许值)，则，假设的upZ即为下泄流量iQ的上游水位，否则，根据逐步逼近方法，用(ZZup)继续试算，直到EPSQQiNjij1.(4)将计算所得的a组(iQ，upZ)值绘成曲线，进行光滑处理，即为联合泄流的上游水位upZ～泄流能力iQ的关系曲线。(5)如果围堰是过水围堰，则当iQPQ(围堰挡水设计流量)，n个泄水建筑物泄流且围堰参与过水。同理给定下游水位downH，按台形堰溢流计算公式计算过水围堰泄流曲线，最后，把计算所得曲线和n个导流建筑物的联合泄水曲线依照上述(1)～(4)方法迭加生成上游水位和泄流能力upZ～iQ的关系曲线。3.施工导流风险分析模型施工导流风险度定义为：在导流施工过程中发生超过上游围堰高程的洪水的频率。水电站在施工截流成功并实现戗堤闭气后，为了保证围堰填筑进度满足第一年洪水的渡汛要求，必须确定分月围堰的填筑高度及相应的风险率[67]。因此，围堰填筑过程中分月堰前水位超过堰顶高程的可能性为：上堰高上HtZPR)(（8）其中：Z上(t)为第t月上游围堰堰前水位；H上堰高为第t月上游围堰堰顶高程。由施工导流风险率表达式可见，导流风险主要由河道来流洪峰流量的不确定性、洪水过程的不确定性及导流设施的泄流能力不确定性[8]所引起的，只有综合反映这三方面因素，才能全面合理地确定导流风险率，因此，对导流中的水文随机性及水力不确定性进行分析。3.1洪峰流量的不确定性分析洪峰流量的不确定性对导流风险的影响最大，我国最大洪峰流量系列采用P-Ⅲ型分布，其概率密度[9]为)(1e)()()(xxxf(δx∞)（9）式中：x为洪峰流量；)21(svCCx；24sC；svCCx2；x为年最大洪峰流量系列的均值；Cv为变差系数；Cs为偏态系数。Cv、Cs均有水文统计资料来确定。3.2洪水过程的不确定性分析洪水过程的不确定性主要表现在洪量的不确定和洪水历时的不确定性。经过大量实测资料统计分析，认为洪量的分布符合P-Ⅲ型分布，其概率密度如式上式。洪水过程历时可以认为服从正态分布，其密度函数为：22221)(uxexf（10）式中u为历时的数学期望，2为历时的方差。3.3泄流能力的不确定分析计算实例表明，导流洞的泄流能力较为接近三角形分布，其密度函数[1011]为：其他0x))(()(2x))(()(2)(cbbcacxcbaacabaxxf（11）式中a为泄流能力下限；b为平均泄流能力；c为泄流能力上限；a、b、c参数通过模拟泄流能力的统计资料来确定。3.4调洪演算模型根据水量平衡方程式：2121212/)(2/)(VVtqqtQQ（13）其中，21,QQ分别为计算时段t始、末入库流量；21,qq分别为计算时段t始、末出库流量；21,VV分别为计算时段t始、末的水库蓄水量。当已知入库洪水过程线时，21,QQ为已知，V1、q1为计算的初始条件，V2、q2未知。又由于，根据导流建筑物的泄流能力曲线H~q和水库的水位库容曲线H~V，可确定q与V的函数关系，即q=f（V）。采用常用的二分法求解方程，推求水库下泄过程线，从而可求得最大下泄流量、防洪所需库容以及相应的最高洪水位。3.5导流动态风险计算确定导流围堰上游水位风险，要系统考虑河道来流、导流建筑物泄流以及枢纽的其他特征，通过系统模拟的方法来实现。由水文随机参数和分布得到一个随机的上游洪水Fin；同时由水力随机参数和分布得到对应的泄流建筑物的泄流能力Qout。通过反复的抽样、计算可以得到一个任意长的围堰上游水位模拟历史系列，上游水位的概率分布只要对这个历时系列进行聚类分析就可以得到了，而设计水位的风险率，实际上就是这个模拟历史系列的密度函数的一些分位点的值。输入水文、水力原始数据、参数拟合洪水过程线统计堰前水位分布产生洪水过程随机数产生泄流过程随机数拟合初期导流系统的泄流过程线系统模型模拟计算上游围堰堰前水位分析在不同围堰高度条件下导流风险率R和保证率P图1导流风险分析计算流程4.实例分析及其计算结果4.1联合泄流能力计算现以金沙江中游河段上某大型水电工程折坝线混合坝初期导流设计方案为例，说明本文所提出的联合泄流水力计算模型及施工导流风险分析模型的应用。混合坝右岸厂房全年导流方案共布置了4条导流隧洞，隧洞净空断面为修正马蹄形，净空最大宽度为15m，底宽11m，净空高度15m，净空断面积191.5587m2。进口高程定在1020.00m，出口高程1012.00m。混合坝河中厂房方案全年导流共布置了7条导流底孔，导流底孔断面为矩形，单孔宽度8m，高度14m（底孔出口段孔口高度压低为12.5m）。底孔前段8m×14m净空断面积112.00m2，后段8m×12.5m净空断面积100.00m2，底孔进口高程定在1020m，出口高程1020m，底孔纵坡i=0。孔口泄流控制断面在出口，8m×12.5m（宽×高）。应用本文所介绍联合泄流水力计算模型计算，联合泄流曲线如图2：050001000015000101010501100Q(m/s)3高程（m）115011801～4导流隧洞＃