静止液体中射流和羽流污染混合区的理论分析
来源:华拓科技网
第31卷第4期 V0L 31 No.4 2010 青岛理工大学学报 Journal of Qingdao Technological University 静止液体中射流和羽流污染混合区的理论分析 武周虎 (青岛理工大学环境与市政工程学院,青岛266033) 摘要:以无限空间静止液体中自由紊动射流和密度差产生羽流的浓度分布理论解为基础,从水环境功能区 管理的角度出发对不同条件下射流和羽流的污染混合区进行了理论探讨和求解.分别给出了平面和圆形射 流以及羽流的污染混合区长度、最大宽度或直径、面积、体积的理式以及外边界等浓度标准曲线或曲面方 程.结果表明,射流污染混合区长度与喷口起始浓度和水环境功能区所执行的标准浓度之比G/G、喷口处断 面几何尺度有关;羽流污染混合区长度与比值 /Cs、起始流量和起始比浮力通量有关;平面和圆形射流、平 面和圆形羽流的污染混合区宽长比依次为0.1317,0.1095,0.1564,0.1310,污染混合区为近似椭圆或近似椭 球形状.可为生产实践中射流和羽流污染混合区的计算提供理论依据. 关键词:静止液体;射流;羽流;污染混合区;标准曲线;解析方法 中图分类号:X143;X52 文献标志码:A 文章编号:1673—46O2(2O1O)04—0012一O7 Analytical Solution to Pollutant Mixing Zone of Jets and Plumes in Static Liquid WU Zhou—hu (School of Environmental and Municipal Engineering,Qingdao Technological University,Qingdao 266033,China) Abstract:Based on the theoretical solution of the concentration distribution of simple turbu— lent jets and plumes produced by density deficiency in infinite space static liquid,the paper gives the theoretical study and solution of the pollutant mixing zone of jets and plumes in dif- ferent conditions from the angle of the management of the water environmental function are— a.The theoretical formula of the length,maximum width or diameter,area,volume of the pollutant mixing zone of plane and round j ets,plane and round plumes,as well as its outer boundary equal concentration standard curve or surface equation are given respectively.The result indicates that the length of the pollutant mixing zone of the j ets is related with the ge— ometry scale in the nozzle cross—section and Co/Ca which is the ratio of initia1 concentration of the nozzle and standard concentration implemented by the water environment function zones, and the length of the pollutant mixing zone of the plumes is related with the C0/C。,initial discharge and initial buoyancy ratio flux.The ratio of width to length of the pollutant miing xzone of plane and round jets,plane and round plumes are 0.1317,0.1095,0.1564,0.1310 in order.The shape of the pollutant mixing zone is approximately elliptical or ellipsoidal, 收稿日期:2010一Ol一06 基金项目:国家自然科学基金资助项目(50979036,50839001) 编辑部约稿 作者简介:武周虎(1959一),男,陕西岐山人.教授,主要从事环境水力学与水污染防治研究.E-mail:wuzh2008@yahoo.corr1.on. 第4期 武周虎:静止液体中射流和羽流污染混合区的理论分析 13 which can provide theoretical basis for the calculation of the pollutant mixing zone of j ets and plumes in production practice. Key words:static liquid;turbulent j et;plume;pollution mixing zone;standard curve;analytic method 工业或城镇污水经处理达到相应的排放标准泄入自然水体后,其近区流动大多属于射流性质.环境质 量标准一般都比排放标准严格,在排污口附近水域会形成污染混合区[1].前人关于无限空间静止液体中等 密度自由紊动射流和密度差产生羽流的研究成果l_2 ],分别给出了射流、羽流的流场和浓度场的理论解,各 断面流速和浓度均服从高斯正态分布,浓度分布曲线比流速分布曲线要平坦些,各断面浓度混合层厚度是 射流边界混合层厚度的入倍.然而,各断面浓度混合层厚度不是依据水环境功能区标准值确定,该厚度是 随射流轴线纵坐标增加线性增大延伸的区域,无法直接评价浓度混合层中的水质状况.在人类生产生活对 水质要求较高的区域,对排污口附近水体污染混合区允许范围的确定,是实施总量控制、确保水环境功能 目标实现的关键所在 J. 国内外对射流和羽流的研究成果,大多集中在浓度场的测量方法及浓度分布特性的试验和数学模型 研究,很少涉及表征水环境功能的污染混合区理式.Birch A D等[6]利用激光腊曼测试系统对自由射 流中的浓度场进行了测量,对射流出口下游的浓度场进行了平均计算并对浓度概率密度的分布进行了统 计;Yoda M等[ 利用LIF技术获得了不同流速下反向流动环境中浓度射流长度的变化规律;黄真理等l8] 采用平面激光诱导荧光技术测量了横流中射流的浓度场;王超E9]采用虚拟扰动源的处理方法,给出了扩散 器排污下游流速和浓度分布的计算公式;樊靖郁等[1 采用RNG湍流模型对浅水环境中含污染物横向射 流近区的三维流场和浓度场进行了数值计算;吴航等_11]应用全场数学模型模拟了深圳市政污水排海工程 的排污混合区范围.武周虎等E12q5 对排污口下游远区的移流扩散问题,采用解析方法推导了简化二维和简 化三维移流扩散的污染混合区最大长度、最大宽度与相应纵坐标、面积以及后者的最大深度等理式; 陈永灿等[1 对三峡水库采用二、三维水质数学模拟方法研究了交汇河段和岸边污染混合区范围.目前,对 污染混合区范围和最大允许污染负荷的计算主要是采用水质模型来实现,虽然水质模型具有地形适应性 强等特点,但不能直观地给出污染混合区范围和最大允许污染负荷与各参数之间的函数关系,给水环境功 能区管理带来不便El 71 8].张丽娜等l】 系统地综述和分析了射流理论在水环境保护中的应用,并指出应从 理论、试验与实际应用等各方面加强对射流的研究. 笔者分别从射流、羽流浓度分布的理论解出发,对静止液体中射流和羽流污染混合区进行理论探讨和 求解,推导污染混合区长度、最大宽度或直径及相应纵坐标、面积、体积的理式以及污染混合区外边界 等浓度标准曲线或曲面方程,进行污染混合区理式的分析与讨论.其研究成果可为相应问题的实际应 用提供简便、快捷的计算公式,也可以对水质模型计算和数据的分析归纳起到理论指导作用. 1污染混合区的理论分析 1.1平面自由紊动射流 从狭长的矩形孔口或缝隙喷出的射流可按平面二维问题分析,当射流出口雷诺数Re一 >30时 可认为射流为紊流.当射流中含有的污染物质对其密度和流场不产生影响,则射流中污染物的浓度分布与 流速分布可以分开计算.在射流主体段的断面上保守物质的浓度分布存在相似性,在背景浓度为零的 无限空问静止液体等密度自由紊动射流中,平面射流的污染物浓度分布为_2 Cm—e一xpp[_l一(一I ) JJ ] ㈩ 其中,射流轴线上浓度C 为 Cn1—2.34√ co (2) 青岛理工大学学报 第31卷 式(2)中的系数2.34是经试验修正的积分常数,下同.在有些新出的环境水力学类书籍中将其作为 理论积分常数,未考虑试验修正因素,笔者根据文献[2—4]原著加注. 射流断面浓度混合层半厚度b ,即射流断面上浓度等于Cm/e点的Y坐标值为 b =2b 一 (3) 式中:41:为沿射流轴线的纵坐标; 为垂直于射流轴线的横坐标;坐标原点取在射流主体段上、下边界延长 线与轴线的相交点O;b。为射流孔口的半高度; 。为射流喷口处起始流速;Co为射流喷口处污染物起始浓 度;6 为射流边界混合层半厚度(即6 ===£z);平面射流边界混合层半厚度扩展系数e=0.154;平面射流浓 度混合层与射流边界混合层的厚度比 —1.41. 污染混合区是从水环境功能区管理的角度出发,针对排污产生的污染物超标区域提出来的概念,是指 由于排污而引起污染物超标水域的影响范围.根据式(1)污染物断面浓度的正态分布和式(2)射流轴线上 浓度Cn 随 -n 减小的特征,令排污引起的允许浓度升高值等于水环境功能区所执行的浓度标准值Cs,则 该等浓度曲线所包围区域为污染混合区.由式(1)一式(3)可得污染混合区外边界等浓度曲线方程为 Cs=2.34√ coexp[一( ) ] 在式(4)中,令 —O可以得到污染混合区长度的计算公式为 ㈤ (5) (6) Ls 34 ×26。(爱) _5.48X2bo(星) 染混合区最大半宽度为 b 一 L 2√e 式(4)两边对z求导,令dy/dx=0,可以求得污染混合区最大半宽度和相应纵坐标的计算公式.则污 污染混合区最大宽度相应纵坐标为 L 一L /√e (7) 将式(5)、式(6)代人式(4)化简得到污染混合区外边界等浓度标准曲线方程为 ( ) 一e(芒) ・n( ) 或 一 8 6 麦√ln( ) (9) 图1给出了静止液体中平面自由紊动射流浓度分布和污染混合区外边界等浓度标准曲线分布.由图 1可以看出,污染混合区外边界等浓度标准曲线近似于椭圆形区域,在靠近排放口一一端出现稍尖形状,在 污染混合区下游边界出现钝头,污染混合区最大宽度出现在纵坐标L 为最大长度L 的1 ≈0.6065处. 污染混合区外边界是采用水环境功能区等浓度曲线确定的,该污染混合区范围反映了污染物扩散稀释后 的缩小过程,在污染混合区内为污染物超过某一水质标准的影响区域,不像浓度混合层厚度随射流轴线纵 坐标增加线性增大无限延伸,浓度单调衰减. 对式(9)在z[O,L ]上求面积积分,可以求得污染混合区的面积计算公式,推导如下: 5===2 √ n(芒)~do:= ̄/eLsb √ n(芒) d(Lx:)。S一√eL b I ̄/ln( )d叩 (10 (11) 进行变量替换,令(x/L ) 一刀,则有 由积分表查得I、//ln( )dr]一√兀/2,代入上式得到类似于椭圆面积的污染混合区面积计算公式 s一÷ ̄/ eL (2b )一0.7306L (2b。) ‘(12) 第4期 武周虎:静止液体中射流和羽流污染混合区的理论分析 乓 童 鼍辞 蟋 纵坐标x/m 图1 自由素动射流浓度分布和污染混合区范围 …一’射流轴线;…一浓度混合层边界线i——|_一射流轴线浓度分布; —污染混合区范围2 一射流横向浓度分佑;— ‘_-一污染混合区范围1;—1.2 圆形自由紊动射流 圆形断面喷口的射流在牛产实践中最为常见,具 轴对称性质,可采用柱坐标分析.在背景浓度为零 的无限空间静止液体等密度自由紊动射流中,圆形射流的污染物浓度分布为 。 I 594Coexp.i…(\ )^£ j’} ~]j ’(13) 式中:r为垂直于射流轴线的径向坐标; 为射流孔口的直径;e一0.114; 一1.12;其他符号同前. 由式(13)可得污染混合区外边界等浓度曲面方程为 _5.59 c。Z exp[L一( \ 门^e _/ J 式(14)为旋转体方程,令r一0可以得到污染混合区长度的计算公式为 L 一5.59d L, ) (15) 式(14)两边对 求导,令dr/'dx一0,可以求得污染混合区最大半径和相应纵坐标的计算公式.则污染 混合区最大半径为 r === L。 (16) 、/2e 污染混合区最大半径相应纵坐标的计算公式与式(7)相同. 将式(15)、式(16)代入式(14)化简得到污染混合区外边界等浓度标准曲面方程为 ( ) =--e( ) n( ) 或 ’, ,一 一 √ ^、/m /n(f) (18) (19) 对式(18)在xEo,L。]上求面积积分,可以求得污染混合区的最大投影面积计算公式,推导同式(12). s一 1 L (2r )一0.7306L (2 ) 对式(18)在z[O,L。]上求体积积分,可以求得污染混合区的体积计算公式,推导如下. 一 e L (芒)。 n( )d(芒) c2。, 16 青岛理工大学学报 第31卷 进行变量替换,令x/L 一 ,进而再令 一 ,则有 V—nelLs j。 In( )d 一一号ne Ls J。ln (21) 由积分表查得1 1n,/dy一一1,代人上式得到污染混合区的体积计算公式 V一 7re L。 (22) 1.3 平面羽流 羽流是由于射流与周围环境流体之间存在密度差,使从喷口出来的射流受到浮力的作用,并且射流的 起始动量很小可以忽略,羽流主要在浮力作用下继续流动和扩散.在背景浓度为零的无限空间静止液体 屯,平面羽流的污染物浓度分布形式与式(1)相同,羽流轴线上浓度Cm为[2 Cm一2.40 qoCo (23) 式中:B。一 gqo为起始比浮力通量,量纲为[L。T-。];q0为起始单宽流量;Pa为环境流体密度; 为环 境流体与射流起始密度之差;z为沿羽流轴线的纵坐标,一般取垂直向上;其他符号同前. 由式(1)、式(3)和式(23)可得污染混合区外边界等浓度曲线方程为 Cs—z.40艉e [一(恚) ] 式中:£一0.147; —1.24. 在式(24)中,令Y一0可以得到污染混合区长度的计算公式为 L。_2.40畿 (25) 式(24)两边对z求导,令dy/dx一0,可以求得污染混合区最大半宽度和相应纵坐标的计算公式.则 污染混合区最大半宽度为 ‘ b。一 L (26) √Ze 污染混合区最大宽度相应纵坐标、污染混合区外边界等浓度标准曲线方程以及污染混合区面积的计算公 式分别与式(7)、式(8)或式(9)和式(12)相同. 1.4 圆形羽流 在背景浓度为零的无限空间静止液体中,圆形羽流的污染物浓度分布为[2 c一11. 7 exp[一 ) ] (27) 式中:B。一a ̄ogQo为起始比浮力通量,量纲为[L r。];Q0为起始流量;e一0.102 ̄2—1.16;其他符号同 前. 由式(27)可得污染混合区外边界等浓度曲面方程为 c 一 7 QoCo exp[一(志)。] (28) 式(28)为旋转体方程,令r一0可以得到污染混合区长度的计算公式为 4.25\OcoC。o)。 (29) 式(28)两边对X求导,令dr/dx一0,可以求得污染混合区最大半径和相应纵坐标的计算公式.则污染混合 区最大半径为 √丧 (30) 污染混合区最大半径相应纵坐标、污染混合区外边界等浓度标准曲面方程、污染混合区面积和体积的计算 第4期 武周虎:静止液体中射流和羽流污染混合区的理论分析 公式分别与式(7)、式(17)或式(18)、式(19)和式(22)相同. 2 分析与讨论 现将静止液体中不同条件下射流和羽流污染混合区主要特性参数列于表I,以便对照分析其变化规 律. 表I 射流和羽流污染混合区特性参数 由表1可以看出,平面和圆形射流、平面和圆形羽流的污染混合区长度分别与喷口处污染物起始浓度 和水环境功能区所执行标准浓度之比Co/c 的2,1,1,0.6次方成正比,平面射流比圆形射流的污染混合 区长度增长快,平面羽流比圆形羽流的污染混合区长度增长快,射流比羽流的污染混合区长度增长快.另 外,射流的污染混合区长度还与喷口处断面几何尺度成正比;而平面羽流的污染混合区长度还与起始单宽 流量成正比、与起始比浮力通量的1/3次方成反比;圆形羽流的污染混合区长度还与起始流量的0.6次方 成正比、与起始比浮力通量的1/5次方成反比. 射流和羽流污染混合区最大半宽度或最大半径与长度、射流边界混合层半厚度扩展系数、射流浓度混 合层与射流边界混合层的厚度比均成正比,污染混合区最大宽度或直径相应的纵坐标均为L。--L。/√e.平 面和圆形射流、平面和圆形羽流的污染混合区最大宽度或直径与长度的比值依次为0.1317,0.1095, 0.1564,0.1310,说明平面羽流的污染混合区相对较宽,显得圆胖一些,而圆形射流的污染混合区相对较 窄,显得细长一些.有趣的是,射流和羽流污染混合区外边界等浓度标准曲线方程式(8)或曲面方程式(17) 具有相类似的形式——近似椭圆或近似椭球.对平面问题为近似椭圆,污染混合区面积的计算公式为式 (12);对轴对称问题为近似椭球,污染混合区的最大投影面积计算公式为式(19),污染混合区的体积计算 公式为式(22). 将式(7)代人式(3)得到污染混合区最大宽度相应断面的浓度混合层半厚度为 b,o一 L。 √e (31) 将射流和羽流相同条件下的污染混合区半宽度6。与相应断面的浓度混合层半厚度6 相除得到6。/6 分别 列于表1.由表l可以看出,该比值在平面射流时为0.5最小,在圆形羽流时为0.913最大.说明浓度混合 层半厚度与污染混合区最大半宽度或最大半径没有一致性关系,前者反映占射流断面最大浓度I/e的混 合层边界半厚度,是一个相对的概念;而后者反映污染物浓度超过某一水质标准区域的最大半宽度或最大 半径.因此,射流和羽流污染混合区长度、最大宽度或直径及相应纵坐标、面积、体积的理式以及污染 混合区外边界等浓度标准曲线或曲面方程,对水环境功能区管理十分有益. 3 结论 笔者分别以无限空间静止液体中等密度自由紊动射流和密度差产生羽流浓度分布的理论解为基础, 18 青岛理工大学学报 第31卷 通过理论求解分别给出了平面和圆形射流、平面和圆形羽流的污染混合区长度、最大宽度或直径及相应纵 坐标、面积、体积的理式以及外边界等浓度标准曲线或曲面方程.结果表明,在静止液体中射流污染混 合区长度与喷V1起始浓度和水环境功能区所执行的标准浓度之比Co/G、喷口处断面几何尺度有关;羽流 污染混合区长度与比值C0/C 、起始流量和起始比浮力通量有关;平面和圆形射流、平面和圆形羽流的污 染混合区宽长比依次为0.1317,0.1095,0.1564,0.1310;最大宽度或直径相应的纵坐标均为L。一 0.6065L ;污染混合区为近似椭圆或近似椭球形状,面积或轴对称问题的最大投影面积为0.7306L。(2b。) 或0.7306L。(2r ),其值相当于椭圆面积的93.O2 ;轴对称问题的体积为1.898,詈L ,其值相当于椭球体 积的90.61 9/5. 参考文献(References): [1]武周虎.倾斜岸水库污染混合区的理论分析及简化条件LJ].水动力学研究与进展:A辑,2009,24(3):296—304. wU Zhou-hu.Theoretical Analysis and Simplified Condition for Pollutant Mixing Zone in Reservoir of Gradient Bank[J].Chin ̄e Journal of Hydrodynamics:Set.A,2009,24(3):296—304. [2]Fischer H B,Imberger J,List E J,et a1.Mixing in Inland and Coastal Watersl'M].New York:Academic Press,1979:315—340. [3]赵文谦.环境水力学[M].成都:成都科技大学出版社,1986:84—108. ZHAO Wen-qian.Environmental Hydraulics[M].Chengdu:Chengdu University of Science and Technology Press,1986:84—108. [4]余常昭.环境流体力学导论[M].北京:清华大学出版社,1992:187—238. YU Chang-zhao.An Introduction to Environmental Hydrodynamics[M].Beijing:Tsinghua University Press,1992:187—238. [5]周丰,刘永,黄凯,等.流域水环境功能区划及其技术关键[J].水科学进展,2007,18(2):216—222. ZHOU Feng,L1V Yong,HUANG Kai,et a1.Water Environmental Function Zoning at Watershed Scale and Its Key Problems[J]. Advances in Water ciSence,2007,18(2):216—222. [6]Birch A D,Brown D R,Dodson M G,et a1.The Turbulent Concentration Field of a Methane Jet[J].J.Huid Mech,1978,88(3): 431—449. [7]Yoda M,Fiedler E H.The Round Jet in a Uniform Counterflow Visualization and Mean Concentration Measurements[J].Experi— ments in Fluids,1996,21:427—436. E8]黄真理,李玉梁,余常昭.平面激光诱导荧光技术测量横流中射流浓度场的研究[J].水利学报,1994(n):1-7. HUANG Zhenqi,H Yu-liang,YU Chang-zhao.Planar Laser-Induced Fluorescence Technique and Its Measurements of oncentCration Field of Turbulent Jet in Cross-Flow[J].Journal of Hydraulic Engineering,1994(11):1f7. [9]王超.多孔射流扩散器在宽浅型江河中排放特性研究[J].水动力学研究与进展:A辑,1993,8(4):426—437. WANG Chao.A Study on Multiport Diffuser Discharge Characteristics in Shallow Open Channels[J].Chinese Journal of Hy& ̄y— namias:Set.A,1993,8(4):426 437. [1O]樊靖郁,张燕,王道增.浅水环境中含污染物横向射流三维浓度分布特性研究l,c]//第十八届全国水动力学研讨会文集.北京;海 洋出版社,2004:590-596. FAN Jing—yu,ZHANG Yan,WANG Dao—zeng.Investigation on 3D Concentration Distribution Features of a Pollutant-Contained Transverse Jet in Shallow Waters[C]//Proceedings of the 18th National Conference on Hydrohynamics.Beijing..China Ocean Press, 2004:590—596. ,1l1]吴航,王泽良,黄剑,等.深圳市政污水排海工程排污混合区范围的确定厂J].环境监测管理与技术,2002,14(6):37—40. WU Hang。WANG Ze-liang,HUANG Jian,et a1.Determination of Pollutant Discharge Region’s Area on Municipal Wastewater Release to Sea Projeet of Shenzhen口].The Administration and Technique of Environmental Monitoring,2002,14(6):37—40. -112]武周虎.明渠移流扩散中无量纲数与相似准则及其应用I-J].青岛理工大学学报,2008,29(5):17—22. WU Zhou-hta The Concept and Its Application of the Dimensionless Numbers&Conform Rules by Adveetion-Diffusion in the Open Channel Flowl'J].Journal of Qingdao Technological University,2008,29(5):17—22. ,1la]武周虎,贾洪玉.河流污染混合区的解析计算方法[J].水科学进展,2009,20(4):544—548. wU Zhou—hu,JIA Hong-yu.Analytic Method for Pollutant Mixing Zone in Riverl'J].Advances in Water Science,2009,20(4):544— 548。 [14]武周虎.水库倾斜岸坡地形污染混合区的三维解析计算方法 ].科技导报,2008,26(18):30—34. wu Zhou_hu.3-D nalAytic Compulational Method of Pollutant M[ixing Zone for Reservoir with Incline nk[J].Science&Technol— ogy Review,2008,26(18):30 34. (下转第4l页) 第4期 申建红,等:基于HHT和wT识别高层建筑模态参数的比较研究 41 [4]Huang N E,Shen Z,Long S R,et a1.The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Non-Stationary Time Series Analysis[C]//Proceeding of Royal Society.London,1998:903—995. Es]Yag nJ N,Lei Y,Pan S W.et a1.System Identification of Linear Structures Based on Hilbert—Huang Spectral Analysis[C]//Part 1: Normal Modes,Earthquake Engineerig&Structuraln Dynamics,2003:1443—1467. [6]陈隽,徐幼麟.HHT方法在结构模态参数中的应用[J].振动工程学报,2003,16(3):383—388. CHEN Juan,XU You—lin.Application of HHT for Modal Parameter Identiifcation to Civil Structures[J].Journal of Vibration Engi— neering,2003,16(3):383—388. ET]Xu Y L,Chen S W,Zhang R C Modal Identification of Di Wang Building Under Typhoon York Using the Hilbert-Huang Transform Method[J].The Structural Design of Tall and Special Buildings,2003,12(1):21—47. [8]Ruzzene M,Fasana A,Garibaldi L,et a1.Natural Frequency and Damping Identification Usig Wavelnet Transfom:Application to Real DataEJ].Mechanical Systems and Signal Processing Journal,1997,11:207—218. -I9]Slavic J,Simonovski I,Boltezar^,L Damping Identifiaticon Using A Continuous Wavelet Transform:Application to Real Data[-J]. Journal of Sound and Vibration,2003,262:291—307. [1o]Le T P,Argoul P.Continuous Wavelet Transform for Modal Identification Using Free ecDay Response[J].Journal of Sound and Vi— bration,2004,277:73—100. [】1] Kijewski T,Kareem八Wavelet Transforms for System Identification in Ciilv Engineering[J].ComputePAided Civil and Infrastruc— ture Engineering,2003,18:339—355. E12]Kijewski T,Kareem Efficacy of Hilbert and Wavelet Transforms for Time-Frequency Analysis[J].Journal of Engineering Me— chanics,2006,132(1O):1037~1049. [13]Peng Z K,Peter W,Chu Tse F L A omparCison Study of Improved Hilbert-Huang Transform and Wavelet Transform:Application to Fault Diagnosis for Rolling Bearing[J].Mechanical Systems and Signal Processing,2005,19:974-988. [14]Li Q S,Xiao Y Q,Fu J Y,et a1.Full Scale Measurements of Wind Effects on the Jin Mao Building[J].Journal of Wind Engineerig nand Industrial Aerodynamics,2007,95:445—446. (上接第18页) [15]武周虎.水库铅垂岸地形污染混合区的三维解析计算方法[J].西安理工大学学报,2009,25(4):436-440. wU Zhou hu.3-D Analytic Computational Method of Pollutant Mixing Zone for Reservoir with Vertical Bank[J].Journal of Xi’an University of Technology,2009,25(4):436—440. [16]陈永灿,刘昭伟,李闯.三峡库区岸边污染混合区数值模拟与分析[C]//中国水力学2000.成都:四川大学出版社,2000:501—508. CHEN Yong-can,LIU Zhao ̄wei,LI Chuang.Numerical Modelling and Analysis of Bank Pollutant Mixing Characteristics in Three Gorge Reservoir[-C]//China Hydraulics 2000.Chengdu:Sichuan University Press,2000:501—508. [17]黄真理,李玉探,陈永灿,等.三峡水库水质预测和环境容量计算[M].北京:中国水利水电出版社,2006:433—438. HUANG Zhen-li,LI Yu-liang,CHEN Yong-can,et a1.Water Quality Prediction and Water Environmental Carr ̄ng Capacity Calcu— lation for Three Gorges Reservoir[M].Beijing:China Water Resources and Hydropower Press,2006:433—438. [183廖文根,李锦秀,彭静.水体纳污能力量化问题探讨[J].中国水利水电科学研究院学报,2003(3):211—215. LIAO Wen gen,LI Jin-xiu,PENG Jig.Dinscussion on Quantification of Pollutant Carrying Capacity[-J].Journal of China Institute of Water Resources and Hydropower Research,2003(3):211-215. [19]张丽娜,王垣.射流理论在水环境保护中的应用[J].水资源保护,2006,22(6):71—75,80. ZHANG Li—na,WANG Xuan.Application of Jet Theory in Water Environment Protection[J].Water Resources Protection,2006,22 (6):7卜75,8O.