精编版-水电站泄流数值计算管理论文.docx

关注我 实时更新 最新资料 水电站泄流数值计算管理论文 水电站采用不同 的泄流消能方式，对水雾 的产生机理、形态及雾量多寡，存在较大 的差异。对于挑流消能工 程，泄流雾化会造成以下危害：厂房进水、断电、交通中断和边坡失稳。自20世纪70年代以来，国内一些专家学者对挑流雾化问题进行了广泛 的研究，取得了丰硕 的成果；但对于底流消能工 程，其泄流雾化问题很少研究。本文就湾塘水电站底流消能雾化 的数计算作一探讨。1底流消能雾化 的数学模型[1]洪水在下泄和消能过程中，由于水流与空气边界 的相互作用，使得水流自由面失稳和水流紊动加剧[2]，进而部分水体以微小水滴 的形式进入空气中，产生某种形式 的雾源。雾源在自然风和水舌风 的综合作用下，向下游扩散，使水雾分布在下游 的一定空间中。之后，水雾经自动转换过程和碰并过程转变为雨滴，以及水雾和水汽之间发生雾滴 的蒸发或凝结过程。如图1所示，因雨滴数较雾滴少得多，故在本数学模型中不考虑雨滴 的蒸发过程和水汽凝结为雨滴 的过程，在图1中用带虚线箭头来表示。1.1水雾雾源量 的计算根据雾源产生 的机理不同，底流消能雾化 的雾源可分为二个；第一是溢流坝面自掺气而产生雾源；第二是水跃区强迫掺气而产生雾源。理论分析[3]和原型观测[4]都表明，后者为主要雾源，故在本数学模型中仅考虑第二雾源，而不计第一雾源对下游 的影响。如图2所示，高速水流流经水跃区发生强迫掺气，其中跃首处旋涡最强，可以认为掺气点发生在此处，从而形成水气两相流。被旋涡挟持进水中 的空气形成气泡，气泡在水中随着旋涡运动，有 的气泡脱离自由面 的束缚以水滴、水雾 的形式跃出水面，从而形成雾源[5]。根据底流消能[6] 的雾化机理，得到下式：式中：ql为单位长度线源 的水雾雾源量，kg/(s·m)；ρ为水 的密度；Lj为水跃 的长度，Lj=10.8hc(Fr1-1)0.93，hc为跃首处 的水深；vc为跃首处 的流速；q为单宽流量，m2/s·m；u′2为跃首处 的脉动速度均方根；uw为自然风和水舌风 的合成风速。选取ρ、Lj、vc为基本物理量，令qe=ρLjvc，利用量纲分析方法式(2)可得：根据湾塘水电站雾化原型观测 的数据[4}，应用逐步回归分析方法[7]，试建立式(3) 的回归模型：以ql/qe为因变量，以Fr1、Nt1和uw/vc作为可能 的自变量，计算表明，在显著性水平为0.05，Fr1、Nt1和uw／vc这三个量对因变量ql／qe 的影响都不显著；以ql为因变量，以qe作为可能 的自变量，在显著性水平为0.05，qe对因变量ql 的影响显著。故对ql和qe作线性回归，求得：1.2水雾扩散 的数学模型1.2.1基本假设(1)水雾雾源位于跃首 的上方，且为连续线源；(2)水雾扩散满足高斯扩散模式，扩散参数采用布鲁克海汶扩散(BNL)参数系统，时空为小尺度模式；(3)水雾在峡谷内扩散，水雾在下垫面发生沉降和反射；(4)地形采用VALLEY(山谷)修正模式。1.2.2风向与线源垂直时水雾 的扩散[8]图4是一个高架连续线源扩散 的示意图，坐标系oxyz 的y轴与坝轴线平行，x轴为水流方向，z轴为垂直向上，点o位于跃首上方，且高程等于下游水位。设P为下游空间 的任意一点，其坐标分别为x、y、z，其水雾 的浓度为：式中：σy为水雾在y方向 的浓度分布方差；σz为水雾在z方向 的浓度分布方差；h为水雾线源 的高度，h=(05～1)(h″c-hc)，h″c为hc 的共轭水深。y1为水雾线源起点y坐标；y2为水雾线源终点y坐标；φ为下垫面 的反射系数。考虑到峡谷内盛行山谷风，并且其风向变化不大。故扩散参数选用布鲁克海纹扩散(BNL)参数系统(阵风度等级为D)：1.2.3风向与线源成任意角时水雾 的扩散在坐标系oxyz中，假定自然风速为uw1，其风向与x轴正向成β1角；水舌风速为vjw，其风向沿x轴 的正向，则自然风速和水舌风速 的合成速度为uw，其风向与x轴正向成β2角，规定：从x轴 的正向开始，绕点o逆时针转动时，角β为正值；反之，角β为负值，如图5所示。建立风坐标系ox1y1z，使x1轴与uw平行，坐标系oxyz、ox1y1z 的z轴相重合。将线源在y1轴上投影，分别得到虚拟线源在y1轴上 的起点和终点坐标：y01=y1cosβ2；y02=y2cosβ2。这样，合成风速uw与线源成任意角 的情况就转化为合成风速uw垂直流过虚拟线源 的情况。参照式(8)，得到下游任意一点 的水雾浓度分布：1.2.4地形 的修正模式因峡谷内盛行山谷风，并且其风向变化不大，故雾流扩散属于中性或弱不稳定 的情形。选取美国国家环保局(EPA) 的VALLEY(山谷)模式，地形 的修正模式主要体现在修正雾源 的排放高度上。在中性或不稳定 的情