法律状态公告日
法律状态信息
法律状态
2017-12-05
授权
授权
2015-06-10
实质审查的生效 IPC(主分类):G06F19/00 申请日:20150119
实质审查的生效
2015-05-13
公开
公开
技术领域
本发明属于湖泊水库水质污染控制领域,具体涉及一种等温层曝气器曝气室表观水流速 度的预测方法。
背景技术
水温分层现象广泛存在于水深较大的湖泊水库,是引起此类水体水质问题的重要诱因.。 一旦湖泊水库水温分层,跃温层会阻碍表层变温层和底部等温层水体之间物质和能量的交换, 同时底部等温层水体由于受各种生物或化学反应而消耗溶解氧,使等温层水体溶解氧浓度逐 渐降低,甚至为零。在缺氧或厌氧条件下,湖库底泥中的有机与无机污染物会向底部水体释 放和溶解,导致湖库水质的内源污染。对于作为饮用水水源的湖库,湖库底泥释放的有机物 质、氮、磷为藻类的生长提供营养盐,在气温较高的季节引起藻类的大量繁殖,增加饮用水 处理难度和成本。
等温层曝气是目前解决分层湖库内源污染的主导技术,主要是在不破坏水温分层的前提 下提高等温层水体溶解氧浓度,抑制湖库底泥中内源污染物的释放,同时给冷水生物提高良 好的栖息场所。其结构可基本简化为一种通用型式(图1),它主要包括进水口、曝气室、顶 部气液分离区、下降区和出水口等。气体通过气体释放器进入曝气室后,气体中的氧与水发 生传质,增加出水的溶解氧浓度,未能溶解的残余气体会缓慢进入等温层之上水体。
曝气室内表观水流速度与曝气量、曝气孔径、曝气室结构、流体特性相关,并直接影响 曝气室内的充氧效率和等温层水体的溶解氧改善效果;但由于对曝气室内气液两相流的复杂 行为缺乏准确合理的数学表征,曝气室内表观水流速度的预测精度欠佳,如Ashley的经验模 型以及Little的能量平衡模型,以及气泡涡流产生的能量损失的表达式不能很准确表达该能 量损失,从而导致绝大多数等温层曝气器设计不合理,不同程度地存在使用效果不理想的实 际问题。有基于此,非常有必要建立一种能准确预测表观水流速度的新模型及其求解方法, 以指导等温层曝气器的设计和运行优化。
发明内容
本发明的目的在于针对上述等温层曝气器运行效果不佳和曝气室表观水流速度预测精度 低的问题和不足,提供了一种等温层曝气器曝气室表观水流速度的预测方法。
为实现上述目的,本发明采用如下的技术方案予以实现:
一种等温层曝气器曝气室表观水流速度的预测方法,包括以下步骤:
1)分析等温层曝气器曝气室内水流运动的净输入能量和总损失能量的组成,确定基于气 体膨胀做功的净输入能量和曝气室和下降区边壁、进出口、顶部、涡流所产生的能量损失的 数学表达式;
2)根据能量守恒原理,建立等温层曝气器曝气室内水流运动的净输入能量和总损失能量 的能量平衡方程;
3)根据气含率的概念,建立等温层曝气器曝气室内气含率与表观气速、表观水流速度和 气液相对滑移速度的方程;
4)联立等温层曝气器曝气室内能量平衡方程和气含率的方程,建立基于MATLAB真域 算法的求解方法,计算可得不同曝气室结构及曝气条件下的等温层曝气器曝气室表观水流速 度。
本发明进一步的改进在于,步骤1)中,具体分为如下步骤:
1-1)根据曝气室内气体膨胀做功产生的总能量等效于在曝气室内外净水头差产生的总势 能的原理,曝气室内水流运动的净输入能量PI用公式(4)计算如下:
PI=PIr-PId={ρl-(ρlεlr+ρgεgr)}gHr-{ρl-(ρlεld+ρgεgd)}gHd (4)
式中:PIr为曝气室气体等温膨胀输入能量,pa;PId为下降区气体等温膨胀输入能量,pa; ρl为水的密度,kg/m3;εlr为曝气室水的体积分率;εld为下降筒水的体积分率;ρg为气体的 密度,kg/m3;εgr为曝气室气体的体积分率;εgd为下降筒气体的体积分率;g为重力加速度, m/s2;Hr为曝气室内水深,m;Hd为下降筒内水深,m;
1-2)曝气室顶部产生的能量损失PT用公式(13)计算;根据等温层曝气器的实际运行数 据,确定了公式(13)中顶部能量损失系数KT的无量纲表达式,并用公式(15)表示;
式中:KT为曝气室顶部能量损失系数;Vlr为曝气室内水流速度,m/s;
式中:Ar为曝气室横截面积,m2;
式中:Ad为下降筒横截面积,m2;
式中:Rer为曝气室雷诺数;
式中:υl为水的粘度,N.s/m2;
式中:Qgr为曝气室曝气量,m3/s;
1-3)曝气室内涡流产生的能量损失PW用公式(14)计算;根据等温层曝气器的实际运 行数据,确定了公式(14)中涡流能量损失系数KW的无量纲表达式,并用公式(17)表示;
式中:KW为曝气室涡流能量损失系数;Dr为曝气室相对轴心的直径,m;D0为中心管 直径,m;Fr为弗劳德数;
式中:g为为重力加速度,m/s2;
1-4)曝气室和下降区边壁的能量损失PFr和PFd分别采用公式(7)和(12)计算:
式中:CFr为曝气室边壁摩擦阻力系数;Ulr为曝气室内表观水流速度,m/s;Ugr为曝气 室内表观气速,m/s;
式中:为下降筒边壁摩擦阻力系数;Uld为下降筒内表观水流速度,m/s;
1.5)曝气室进出口产生的能量损失PE用公式(19)计算:
式中:KEn为曝气室底部入口处局部水头损失系数;KEx为下降筒底部出口处局部水头损 失系数;UEx为下降筒底部出口处表观水流速度,m/s;
式中:AEx为下降筒底部出口处横截面积,m2。
本发明进一步的改进在于:步骤2)中,等温层曝气器曝气室内水流运动的净输入能量 和总损失能量的能量平衡方程如下式所示:
PI=PT+PE+PW+PF (1)
式中:PI表示曝气室内水流运动的净输入能量,PT、PE、PW、PF分别表示曝气室顶部、 进出口、气泡涡流和曝气室边壁产生的能量损失。
本发明进一步的改进在于:步骤3)中,等温层曝气器曝气室内气含率与表观水流速度 之间关系表示如下:
气液相对滑移速度Vb与气泡直径rb的关系如下:
当rb<7×10-4m时,Vb=4474rb1.357;
当7×10-4≤rb≤5.1×10-3m,Vb=0.23; (22)
当rb>5.1×10-3m,Vb=4.202rb0.547。
与现有技术相比,本发明具有如下的技术效果:
本发明通过分析等温层曝气器内气水两相流运动的能量平衡及气液漂移速度,建立无量 纲的曝气室顶部能量损失系数KT和涡流能量损失系数KW,准确表征曝气室内较大的顶部能 量损失和涡流能量损失,建立一维等温层曝气水动力学模型及其基于MATLAB的解析方法, 并成功应用于美国Prince湖的等温层曝气器曝气室水流速度的预测。
具体来说,本发明系统分析曝气室内气液两相流运动的能量平衡关系,确定了曝气室内 顶部能量损失和涡流能量损失的数学表达式以及对应的无量纲能量损失系数,并将气含率与 表观水流速度进行关联,得到由能量平衡方程和气含率方程组成的曝气室一维水动力学模型, 利用MATLAB软件的真域算法求解该水动力学模型,预测曝气室内表观水流速度。
本发明曝气室一维水动力学模型及其求解方法,针对美国Prince湖实际应用的等温层曝 气器结构和曝气量条件,成功预测了不同曝气量条件下曝气室表观水流速度,预测误差低于 8%,远优于目前采用其他方法的预测误差20%。
并且,本发明曝气室一维水动力学模型及其求解方法,可用于预测不同曝气室结构及曝 气量、气泡直径条件下的曝气室表观水流速度。
附图说明
图1为等温层曝气器的通用构造示意图;
图2为采用本发明对曝气室表观水流速度的预测值与实测值对比图;
图3为忽略涡流能量损失时曝气室表观水流速度的预测值与实测值对比图;
图4采用其他方法对曝气室表观水流速度的预测值与实测值对比图;
图5为采用本发明对不同气泡直径条件下曝气室表观水流速度的预测值图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
参见图1,本发明一种等温层曝气器曝气室表观水流速度的预测方法,具体包括如下内 容:
水动力学模型的构建方法
系统分析气液两相流运动的能量平衡,建立等温层曝气的一维水动力学模型及主要模型 参数的无量纲表达式,根据曝气室内气、液运动的驱动能量来自气体释放器释放出的压缩气 体的膨胀做功的原理,曝气室内净输入能量等效于曝气室内、外水体密度差产生的有效作用 能量,并和曝气室内气、水两相流运动的总能量损失(包括曝气室顶部、边壁、进出口及气 泡涡流产生的能量损失)相平衡,由此建立曝气室内水流的能量平衡如方程(1)所示:
PI=PT+PE+PW+PF (1)
式中:PI表示曝气室内水流运动的净输入能量,PT、PE、PW、PF分别表示曝气室顶部、 进出口、气泡涡流和曝气室边壁产生的能量损失。
1、曝气室气体等温膨胀的净输入能量:曝气室、下降区的气体等温膨胀推动力为气、水 之间的密度差,上升段流体为气水二相流,由于气体不符合理想状态条件,故不适宜采用理 想气体等温膨胀做功公式计算输入能量,气体等温膨胀输入能量可用下述公式表示:
PIr=ρlgHr-(ρlεlr+ρgεgr)gHr (2)
式中:PIr为曝气室气体等温膨胀输入能量,pa;ρl为水的密度,kg/m3;εlr为曝气室水 的体积分率;ρg为气体的密度,kg/m3;εgr为曝气室气体的体积分率;g为重力加速度,m/s2; Hr为曝气室内水深,m;
PId=ρlGHd-(ρlεld+ρgεgd)GHd (3)
式中:PId下降区气体等温膨胀输入能量,pa;εld为下降筒水的体积分率;εgd为下降筒 气体的体积分率;Hd为下降筒内水深,m;
曝气室内水流运动的净输入能量PI为:
PI=PIr-PId={ρl-(ρlεlr+ρgεgr)}GHr-{ρl-(ρlεld+ρgεgd)}GHd (4)
液体含率和气体含率之和为1,即:
εgr+εlr=1 (5)
εld+εgd=1 (6)
2、曝气室和下降区边壁的能量损失:曝气室边壁的能量损失PFr可用下式表示:
式中:CFr为曝气室边壁摩擦阻力系数;Ulr为曝气室内表观水流速度,m/s;Ugr为曝气 室内表观气速,m/s;Dr为曝气室相对轴心的直径,m;D0为中心管的直径,m;
式中:为曝气室雷诺数;
式中:υl为为水的粘度,N.s/m2;
式中:Qgr为曝气室曝气量,m3/s;Ar为曝气室横截面积,m2;
下降区边壁的能量损失PFd同样可表示如下:
式中:为曝气室边壁摩擦阻力系数;Uld为下降筒内表观水流速度,m/s;
3、曝气室顶部和气泡涡流产生的能量损失占总能量损失的比例很大(可达20%-50%以 上),而目前缺乏相关准确的顶部和涡流能量损失计算方法,本发明提出一种新的顶部和涡流 能量损失PT和PW表达式,分别用公式14、15表示:
式中:KT为曝气室顶部能量损失系数;Vlr为曝气室内水流速度,m/s;
式中:KW为曝气室涡流能量损失系数;
4、目前国内外对等温层曝气室顶部能量损失系数KT的研究相对较少,本发明利用等温 层曝气器的实际运行数据(表2),拟合一个含有Ad/Ar及Rer的无量纲表达式(公式16)来 表征KT,应用公式16计算得到不同工况下KT的值与根据表2实测能量损失计算的值吻合良 好(表3);涡流损失系数KW与曝气室尺寸结构Hr/(Dr-D0)、气含率εgr、弗劳德数Fr等因 子有关,本发明拟合一个包含以上无量纲参数的无量纲表达式(公式17)来表征KW,应用 公式17计算得到不同工况下KW的值与根据表2实测能量损失计算的值吻合良好(表4)。
式中:Ad为下降筒横截面积,m2;
式中:g为重力加速度,m/s2;
5、曝气室进出口产生的能量损失分别用公式19表示,其中出水口和曝气室表观水流速 度的关系用公式(20)表示:
式中:KEn为曝气室底部入口处局部水头损失系数;KEx为下降筒底部出口处局部水头损 失系数;UEx为下降筒底部出口处表观水流速度,m/s;
式中:AEx为下降筒底部出口处横截面积,m2。
6、下降区的气含率对曝气室内的水流速度影响不大,故在本解析模型中忽略其影响,即 εgd=0。曝气室内的气含率与表观水流速度之间关系表示如下:
式中:Vb为气液相对滑移速度;
气液相对滑移速度Vb与气泡直径rb的关系如下:
当rb<7×10-4m时,Vb=4474rb1.357;
当7×10-4≤rb≤5.1×10-3m,Vb=0.23; (22)
当rb>5.1×10-3m,Vb=4.202rb0.547;
7、联立式(1)-(20),并利用和化简得:
联合公式(19)和(23),得到曝气室内水流运动的一维水动力学模型,利用MATLAB 软件的真域算法求解该水动力学模型,可预测曝气室表观水流速度。
本发明一种等温层曝气器曝气室表观水流速度的预测方法,其核心思想是根据曝气室内 气、液运动的净输入能量来自气体释放器释放出的压缩气体的膨胀做功的原理,曝气室内净 输入能量等效于曝气室内、外水体密度差产生的有效作用能量(公式4),并和曝气室内气、 水两相流运动的总能量损失(包括曝气室顶部、边壁、进出口及气泡涡流产生的能量损失) 相平衡,由此建立曝气室内水流的能量平衡如方程(公式1);曝气室顶部及涡流产生的能量损 失占总能量损失的比重较大,重点建立了曝气室顶部及涡流产生的能量损失表达式,并利用 等温层曝气器的实际运行数据,确立了曝气室顶部及涡流能量损失系数的无量纲表达式,此 两项能量损失采用公式(14-18)进行分项计算;曝气室边壁的能量损失用公式(7-13)计算, 曝气室进出口的能量损失采用公式(19-20)计算;曝气室内的气含率与表观水流速度之间关 系用公式(21)表示,其中气液相对滑移速度与气泡直径的关系用公式(22)表示;联立公 式(1-20)得到公式(23),联合公式(21)和(23)得到曝气室水流运动一维水动力学方程, 建立基于MATLAB真域算法的模型求解方法,可以预测不同曝气室结构和曝气量、气泡直径 条件下的曝气室表观水流速度。
对照图1,美国Prince湖实际的等温层曝气器曝气室结构尺寸如表1所示;实际的曝气 量、实测表观水流速度和气含率等运行数据如表2所示,具体数据来源参见其他文献(Burris V L,Mcginnis D F,Little J C.Predicting oxygen transfer and water velocity in airlift aerators[J].Water Research,2002,36(18):4605–4615)。
表1:实际的等温层曝气器结构尺寸表。
该等温层曝气器中气泡直径db在2.1—3.1mm之间,选取气泡上升滑移速度为0.23m/s, 即取Vb=0.23m/s。曝气器进出口能量损失系数KEn、KEx参见Chisti的指导手册,即取KEn=1.4、 KEx=1.5。
根据表1和表2所示数据,拟合出无量纲的曝气室顶部和涡流能量损失系数(公式16-17); 结合表1和表2所示数据,分别应用公式16和17预测了不同工况条件下该等温层曝气器曝 气室顶部和涡流能量损失系数,分别如表3和表4所示,计算值与实算值的相关系数(R) 分别为0.99和1.00,说明拟合公式的预测精度非常良好。
表2:实际的等温层曝气运行数据表。
表3:KT的拟合值与实际计算值的对比表。
表4:Kw的拟合值与实际计算值的对比表。
应用本发明所提出的方法,对该等温层曝气器曝气室表观水流速度的预测结果如图2所 示,预测值与实测值的吻合非常理想,表观水流速度的预测误差低于8%。如果忽略曝气室内 涡流能量损失,曝气室表观水流速度的预测误差将提高至46%,如图3所示)。如果采用其 他方法,如气体膨胀做功量采用理想气体的等温膨胀做功量来计算、曝气室内涡流能量损失 采用计算,则曝气室表观水流速度的预测误差将提高至25%,如图4所示。
应用本发明所提出的方法,针对该等温层曝气器的曝气室结构尺寸和曝气量,分别采用 气泡直径2mm、0.2mm和0.02mm,曝气室表观水流速度的预测结果如图5所示。当气泡直 径为2mm,当曝气量从0.018m3.s-1逐渐增加到0.063m3.s-1时,预测的表观水流速度相应从 0.44m.s-1增加到1.13m.s-1;在曝气量一定的条件下,当气泡直径从2mm减小10倍时,表观 水流速度与气泡直径负相关;当曝气量从0.018m3.s-1逐渐增加到0.063m3.s-1时,表观水流速 度的增加幅度从34%逐渐减小到9%。当气泡直径从0.2mm进一步减小10倍时,表观水流速 度与气泡直径基本无关,当曝气量从0.018m3.s-1逐渐增加到0.063m3.s-1时,表观水流速度的 增加幅度几乎为0。此预测结果与理论预期一致,气泡直径对表观水流速度的影响结果也基 本与国内外其他研究成果类似。
本发明可用于预测不同等温层曝气器曝气室结构和曝气条件下曝气室表观水流速度,为 等温层曝气器的优化设计和运行提供科学可靠的指导。
机译: 一种淬灭入射流的动能的方法,一种水流的能量耗散池(选件),一种用于通过其横截面区分液体速度的装置,一种差能淬灭器,一种淬灭水蒸气的方法。液体的动能,保护河流免受溢洪道侵蚀的方法,保护海岸免受侵蚀的沿海地区,保护海岸免受侵蚀的设备,冲洗因污染造成的海岸污染的方法沿海地区海浪的能量,一种由于nergii沿海地区海浪冲刷海洋的设备
机译: 一种水处理装置,其通过从过滤层上部返还一些原水,并通过搅拌叶片使水流通过过滤层,通过搅拌叶片搅拌过滤介质,从而在过滤层上方形成涡流。
机译: 连续可变传输的输入旋转速度预测装置,输入惯性计算装置的输入惯性扭矩输入,控制装置中装有任意一种的输入装置,输入旋转速度预测方法和输入惯量