毕业论文-沿海油库海水淡化后浓海水处置技术的研发排海扩散器的设计

毕 业 论 文（设计） 题 目 ： 沿海油库海水淡化后浓海水处置技术的研发 排海扩散器的设计 学 院： 石油化工学院 学生姓名： 虞邦炜 专 业： 油气储运工程 班 级： A06储运 指导教师： 竺柏康 起 止 日期： 2011.2.21-2011.4.29 2011 年 4 月 25 日 目录 中文摘要 . I 英文摘要 . 1 前言 . 1 1.1 背景 . 1 1.2 浓海水处理方法 . 1 1.3 浓海水处理排海工程 . 1 1.4 本文的工作 . 2 2 海水淡化的尾液影响 . 3 2.1 浓海水对海洋生态系统的影响 . 3 2.2 海水淡化厂污染物排放 . 3 2.3 排放水物理性质变化对海洋生态系统的影响 . 3 3 模型理论 . 4 3.1 近区模型研究现状 . 4 3.2 远区模型研究现状 . 6 3.3 水动力学参数 . 7 3.4 海洋水动力学状况评价方法 . 7 4 排海扩散器的设计 . 11 4.1 扩散器结构与要求 . 11 4.2 扩散器结构概念设计的原理和流程 . 13 4.3 设计中主要影响因素 . 14 4.4 孔口出流系数法 . 14 4.5 动压力水头法及其优化设计程序 . 18 4.6 六横岛台门海水淡化工程 . 22 5 总结与建议 . 28 5.1 总结 . 28 5.2 建议 . 28 6 参考文献 . 29 I 沿海油库海水淡化后浓海水处置技术的研发 排海扩散器的设计 虞邦炜 (浙江海洋学院石油化工学院，浙江 舟山 316004) 摘要 浓 海 水处理排海工程就是将海水经过淡化处理的浓海水通过海洋放流管输送到离海岸一定距离、一定深度的 强流区域，由海洋放流管尾部的扩散器排放，即充分利用海洋的扩散、降解和自净能力，达到浓海水处理的 目 的。 海水经淡化之后，产生浓盐水尾液，如果它不加处理，直接排放，势必对环境造成严重影响。利用扩散器排放浓海水，能很好的利用自然海域的特点，海域本身的自净溶解能力使浓海水及时得到扩散稀释，达到净化目的。这样，既解决了环保问题，又使处理浓海水尾液成本得以降低。本文从我国水资源的严重缺乏角度，阐述了海水淡化的重要性，进而引出海水淡化之后所产生的浓海水尾液的影响。通过对近些年排海工程理论模型的了解，再对海洋水动力学进行研究 ，提出一些模拟方法。本文着重阐述如何进行排海扩散器的设计，从结构、要求、设计原理、流程和主要影响因素等多角度考虑，并介绍两种扩散器水力设计方法 孔口出流系数法和 动压力水头法。 在 六横岛台门海水淡化工程 中，采用 T 型走向扩散器 ， 通过数值模拟，对 轴线流速、浓度及稀释度的沿程变化 进行 比较 ，得出结论。 关键词 扩散器；浓海水尾液；孔口出流系数法；动压力水头法； T型走向扩散器 II The coastal depot sea water desalination with the r&d disposal technology,the emissions sea diffuser design Yu bangwei (Institute of petrochemical technology, Zhejiang Ocean University, Zhoushan 316004) Abstract Concentrated seawater discharging project is that seawater through desalinaton by marine discharging pipe deliver to thet coast which is distance from a certain depth of the high current region, discharge by the diffuser of sea discharging pipe rear, or take by the full advantage of the ocean spread, degradation and self-purification ability, so as to achieve the purpose of concentrated seawater. After the water desalinated,then produced the concentrated seawater, if it is not addressed, discharges directly, is bound to have serious impact on the environment. Concentrated sea water discharges by using the diffuser, the use of well characteristics of natural waters, waters dissolve their self-purification ability to dilute the concentration of water diffusion in time to achieve purpose of purification. In this way, not only to solve the environmental problem, but also to deal with the cost of concentrated sea water solution.From the perspective of a serious shortage of water resources in China, explaining the importance of desalination, and then leads the generated after the end of solution of concentrated sea water. Through the recent discharge engineering of the theoretical model to understand, and to study the dynamics of ocean water, we give out some simulation. From the structure, requirements, design principles, processes, and the main factors and many other considerations, this article focuses on how to design the diffuser, and describes two methods of hydraulic design of the diffuser-orifice the flow coefficient and the dynamic pressure head method. The desalination projects of island of Taimen Liuheng, through numerical simulation with T diffuser, compared with the axis velocity, concentration and dilution,then make the conclusions. Key words diffuser； thick liquid water； the flow coefficient method of the orifice；dynamic pressure head method； T toward of the diffuser 1 1 前言 1.1 背景 淡水资源短缺是全球目前面临的主要社会问题之一，解决淡水短缺问题除了传统的节约用水、废水利用、远途调水等方法外，利用现代技术大规模开辟新的水源则首推海水淡化技术。 我国是一个水资源严重短缺的国家，人均淡水占有量仅为世界平均水平的 1/4。沿海地区经济发达、人口稠密，淡 水供需矛盾更加突出。在陆地淡水资源日益短缺的严峻形势下，人们把目光投向了海洋。我国海水淡化技术的研究起步较早， 1967 1969年全国组织海水淡化会战，同时开展电渗析（ ED）、反渗透（ RO）和蒸馏多种海水淡化方法的研究。 1981年建成西沙 200m3/d电渗析海水淡化装置； 1997年，浙江省重大科技攻关项目“ 500m3/d反渗透海水淡化示范工程”在浙江省嵊泗县嵊山岛建成投产 1； 2000年，在国家科技部重点科技攻关项目“日产千吨级反渗透海水淡化系统及工程技术开发”的支持下，先后在山东长岛、浙江嵊泗建成了 1000m3/d反渗透海水淡化示范工程； 2003年，国家发改委高技术产业化项目“山东荣成日产10000吨级反渗透海水淡化示范工程”一期 5000m3/d机组在荣成市石岛建成投产； 2004年，国家科技部科技攻关项目“低温多效海水淡化示范工程”， 3000m3/d低温多效海水淡化装置在青岛市黄岛电厂建成 2。经过近 40年的研发和示范，我国海水淡化技术已日趋成熟，为大规模应用打下了良好基础。我国已成为世界上少数几个掌握海水淡化先进技术的国家之一 3。目前已建成运行的海水淡化水产量约为 12万 m3/d， 在建和待建的工程规 模为 38万 m3/d。根据国家海水利用专项规划，我国海水淡化能力 2010年将达到 80万 100万 m3/d4， 2020年将达到 250万 300万 m3/d。 1.2 浓海水处理方法 世界上常用的浓海水处理方法可分为两大类。一类是直接排放，如排入海洋、地表水、污水处理系统等；第二类将浓海水进行再利用，如地表灌溉、制盐、提取化工原料等 5。其中，直接排放投资少、效益高，但其处理不当会对环境造成严重污染，故一般采用扩散器来加速浓海水的稀释，达到保护环境的作用。 1.3 浓海水处理排海工程 浓 海 水处理排海工程就是将海水经过淡化 处理的浓海水通过海洋放流管输送到离海岸一定距离、一定深度的强流区域，由海洋放流管尾部的扩散器排放，即充分利用海洋的扩散、降解和自净能力，达到浓海水处理的 目 的。排海工程的规划、设计、施工、运行，在国外已有几十年的历史。早期建造的排海工程一般只是一条简单的放流管，末端开口，不带扩散器，排放的污水也不经过任何预处理。至上个世纪 20年代，排海工程开始在放流管的末端增加了一段带有多孔的扩散器，而且也逐步开始了排海前的预处理。例如，建于 1925年的英国 Heng istbury污水排海工程安装了一个有 6个喷孔的扩散器。利 用多孔扩散器将污水分散排放更有利于污水与海水的 混合 ，它通过其喷嘴将浓 海 水射入环境水体，是一种用来增强污水与环境水体掺混稀释能力的工程措施。 合理 地设计扩散器可使射出的浓 海 水 在较小范围内获得高倍数稀释，并可以提高稀释效率，降低局部海区的污染程度，避免形成稳定的浓海水场，便于 2 浓海水的进一步输移扩散。这样既充分利用纳浓海水体的环境容量和自净能力，显著减少浓海水处理费用，又保证环境目标的实现 6。 1.4 本文的工作 本文应用数值模拟的方法对海水淡化浓海水尾液排放过程进行研究，预测污染物对环境可能造成的影响范围和程度，针对 浓海水尾液，设计出合理的扩散器形式，并提供可行的减轻不良环境影响的措施，为管理部门提供决策依据，为开发工程的环境保护设计提供科学依据，为作业者实施完善的环境管理措施提供可操作依据，对沿海社会、经济的可持续发展也有着较为重要的现实意义。 3 2 海水淡化的尾液影响 海水淡化可为内陆地区节省更多可以利用的淡水资源，这对于长远解决我国水资源短缺问题具有战略意义。但在大力发展海水淡化事业的同时，也不可避免地带来海洋生态环境问题。海水淡化厂排放的浓海水及其所含的污染物 (重金属、化学添加剂等 )以及物理性质（如温度、密度 )的变化，如未经适当的处理而直接排放入海，将对海洋生态环境造成相当的冲击。 2.1 浓海水对海洋生态系统的影响 海水淡化厂排放的浓海水的盐度一般是取用海水的 2倍 7。若这些浓海水排放方式不当，将导致排放海域盐度的升高。以胶州湾为例，按海水淡化水产量 20万 m3/d、胶州湾海水交换周期为 60天计算，若产生的浓海水全部排入胶州湾，则胶州湾的平均盐度将每年上升约 0.3个盐度单位， 30年后胶州湾的平均盐度将超过 40，与死海的盐度相当。 盐度的升高会改变海洋生物本身体液与其生活 环境海水中渗透压的平衡，从而降低海洋生物的繁殖力 (主要是幼虫和幼仔 )，甚至使其灭绝。研究发现许多类海洋生物的呼吸及排泄能力，都与其周遭环境的盐度有密切的关系。有些海洋生物被称为狭盐性，因为它们仅能在一个狭窄的盐度变化范围内保持其体液与周围环境间渗透压的平衡。而能忍受环境中较大盐度变化的海洋生物，则被称为广盐性。有些海洋生物虽然在盐度增高至某一程度时仍能生存，但其细胞的增殖能力却已大为减低。此外，由于底栖生物无足够的移栖能力，因此浓海水排放对排水口附近的底栖生物的影响尤其严重。 2.2 海水淡化厂污染物排放 海水淡化 厂排放水中污染物来源主要有两类：一类是化学添加剂，如生物杀灭剂 (通常为氯气或次氯酸钠 )、抑垢剂 (通常为聚磷酸盐 )、防沫剂、防蚀剂、酸洗剂等；另一类是由管路腐蚀产生的毒性重金属，如 Cu、 Ni、 Mo、 Cr、 Zn 等。研究表明，以上这些污染物都会对海洋生态系统产生危害。以生物杀灭剂 氯气为例，多级闪急蒸馏法排放水中游离氯含量一般在 (0.2 0.5)10-6，若按淡水 /排水 (含冷却水 )比 1 9 来估算，一个 10 万 m3/d 的淡化厂排放游离氯量为 180450kg/d。虽然稀释作用和降解作用会降低接受水体的游离氯浓度 ，但即使很低的游离氯浓度也会对海洋生物产生毒害作用。首先，游离氯是一种高效的生物杀灭剂 (这也是作为海水淡化生物杀灭剂的原因 )；其次，游离氯与海水中的有机物发生化学反应，产生若干有致癌作用或毒害作用的卤化物。美国国家环保局规定海水中游离氯不得超过 0.01310-6(短期 )和 0.0075 10-6(长期 )。另外，在浓海水中浓度最大的重金属是铜，其含量比自然海水高 12 个数量级 8。 2.3 排放水物理性质变化对海洋生态系统的影响 排放水物理性质的改变主要有两点：温度升高和密度增大。排放水密度的增大主要影响接受水体的 物理性质，由于浓海水的密度大于自然海水，其入海后易于沉降在水底，阻碍了海水的垂直混合，并在排水口附近形成高盐沙漠。而排放水温度的升高 (尤其是蒸馏法 )则直接影响海洋生物的生长和繁殖。大部分的海洋生物都是在一定的温度范围内生长和繁殖，温度的改变会影响海洋生物的生理机能，并影响其产卵、生长及幼虫孵化能力。此外，排放水温 4 度的升高将导致接受水体溶解氧含量的降低，而间接对海洋生物和水质产生不利影响。 3 模型理论 3.1 近区模型研究现状 国内外对城市污水海洋处置工程近区的稀释扩散规律的研究已有较长的历史，与污染近区有关的射流理论 也由规则边界中静止环境内的平面与单孔射流向复杂流动中的复杂射流发展，如横流、分层流、浅水域射流，潮汐流中的多孔射流、表面射流、旋动射流等。 射流是指一般流体从各种排放口或喷口流入周围环境流体，并同其发生强烈混合的流动状态。根据射流的形成可将射流分为三种：喷口处初始动量对流动起支配作用的射流为纯射流或动量射流，由于喷口处流体与环境流体的密度差产生的浮力作用形成的射流为羽流和卷流，而浮射流则是即受动量作用又受浮力作用，同时具有射流和羽流特性的一种射流。射流的基本特征就是由于射流与周围环境流体之间速度差形成的卷吸 和掺混，使射流断面不断扩大，射流流速沿程减小，浓度或温度降低。目前，射流理论被广泛的应用在水利、航空航天、环保、冶金、化工、交通、矿山等多个领域。 浮射流是环境水力学和流体力学近年来的重要研究对象。从上世纪 20 年代开始研究无限空间同类流体中等密度自由紊动射流理论以来，国内外学者进行了大量有关浮射流的研究工作，使浮射流理论有了较大的发展。浮射流特性研究的目的主要是在于确定它的稀释度、射流轴线的轨迹、射流扩展的范围和射流中流速的分布等。对于变密度、非等温和挟带有污染物质的射流还要确定密度分布、温度分布、和挟带物 浓度分布。研究方法主要有物理模型和数学模型两种。 3.1.1 物理模型 物理模型即根据水力相似性原理，采用将原型缩小到一定比尺后进行污染物排放物理实验，通过实测数据分析污水经排放口出流后浮射流的行为特征，从而得到描述其运动规律的经验半经验公式。 量纲分析是物理模型常采用的方法，通过分析影响射流运动的诸多因素，如环境水体因素、扩散器几何要素、浮射流自身排放要素等，得到描述浮射流运动特性的关系式。在浮射流初始研究阶段，量纲分析方法在分析整理试验资料和数据上发挥了巨大的作用，获得了有关浮射流的许多基本特征和基本参 数，为人们认识和研究浮射流现象提供了强有力的手段，即使在当今的江河湖海尾水排放的物理模型实验中，该方法仍不失为有效的手段。 List 和 Papanicolaious（ 1988）和 Fischer（ 1979）在分析了大量前人的研究过后分别确定了射流区和羽流区内的速度和浓度的高斯分布函数中的指数系数，以及断面扩展半宽度等特征参数。他们还对各流速沿轴线的衰减规律进行了研究，并在采用量纲分析方法对实验数据整理分析的基础上进一步证实了时均速度和浓度断面符合高斯分布形式。 国内的许多专家和学者利用该方法对不同形式和排放条件 的污水排放口进行了模型试验研究，取得了大量的宝贵资料和工程经验，对排放口的设计和工程运行起到了很好的指导作用，同时也为其它类似污水排放口的建设提供了宝贵的经验，推动了我国污水江河湖海处置技术的发展。采用该方法，韦鹤平（ 1989）分别对上海星火工业园区污水排海工程排放口模型进行了实验研究，得出了适用于该工程的排放口参数 9；徐高田（ 1997）结合上海河流污水 5 治理二期工程白龙岗排放口水利模型试验，采用量纲分析法对近区初始稀释度的变化规律进行了研究，得出了描述该排放口条件下初始稀释度变化规律的经验公式；严忠民、 蒋传丰等人（ 1991）还依据无量纲分析结果，对污水排放管的物理模型试验准则和实验方法进行了有益的探索。 3.1.2 数学模型 数学模型主要包括积分模型和场模型，以及长度比模型。 1）积分模型 数学模型中最经典的是射流积分模型，它始于 20 世纪 60 年代后期， 70 年代有较快的发展，提出的模型有 20 多种。射流积分模型的基本思路是：首先假定沿射流各个断面上的流速分布、浓度（温度）分布相似性，通常将射流近区分为初始段和主体段，初始段是射流流态形成区，沿程侧向断面上各因变量的变化多用椭圆分布近似，主体段是射流流态发展区， 沿程侧向断面上因变量的分布多用高斯分布逼近。其次必须对射流的厚度作线性扩展假定或对浮射流从侧边卷吸流体的流量作出一个卷吸假定，常用卷吸系数又可近似是常量。对有环境来流的情况还应确定绕流阻力系数。射流积分模型是上述假定的基础上，根据质量、动量和物质守恒的原理，对与射流轴线垂直的横断面积分获得的。然而，积分模型由于采用了许多对称假设，使之不能很好的应用于边界条件较为复杂的情况，只能适用于简单的流动计算。Morton， Taylor 和 Turner（ 1973）提出了半无限水域、水体静止且密度均匀的理想环境中单孔或长孔 射流研究的积分模型， List 和 Imberger 以及 Lee 等应用该方法对不同环境水体特性条件下的浮射流运动规律进行了研究，取得了大量成果，丰富了积分模型理论。我国学者韩保新（ 1996）等人研究了一种污水排海的近远区计算与控制方法，并应用于大亚湾水容量计算和污水排海。槐文心和李炜（ 1993）应用此方法研究了浮射流卷吸模式本身的特性。徐田高和韦鹤平（ 2000）研究了积分控制模型，并应用于嘉兴污水海洋处置工程中。 从文献可以看出，有关排污射流的研究其环境流体为恒定横流的较多，并取得了丰富的成果，但实际上多数排污工程 的受纳域是受潮汐作用强烈的河流和近海水域，因此对非恒定横流作用下射流性质的研究更具有实际意义，但就目前的科学研究的状况来说，认识和描述脉动量及非恒定问题存在很大的困难，对潮流中的排放问题的实验研究也较少，因此也只能在恒定横流研究的基础上近似模拟非恒定横流中污染物排放的射流问题。 2）场模型 在浮射流的近区，水流紊动剧烈，有强剪切层存在，一般的时均化模型较难于应用，特别是在动水条件下，采用紊流模型才能较准确的把握。紊流模型较上述积分法的应用范围要广，并可以考虑复杂的流动，考虑潮流的动态变化对浓度场的影响，但是 紊流模型相对比较复杂，其发展尚不成熟。紊流方程一般分为：零方程（混合长度方程），单方程（紊动动能方程 K 方程），双方程（ K f， K方程等）。还有很多学者对紊流模型进行了改进： Rodi 首先采用了代数应力 /通量模型对热电厂的温排放浮羽流特性进行了模拟； Chen 和 Li 采用此模型及改进的单尺度、多尺度模型，对不同分层环境射流进行了较好的模拟；此外 Gu Jie 等采用模型结合尾水排放近区的模拟； Kim 等应用垂向坐标变换对多孔扩散器近远区进行了数值模拟。在“七五”期间，我国学者深入研究了潮汐流动中底部排放污染混和区近区 K模型，岸边排放污染混合区完全深度平均模型及实用化计算方法等课题，取得了一定的成果。 3）长度比模型 自 80 年代以来，一种以量纲分析为基础来对流动进行分区的长度比模型得以开发和应用。 6 它容易处理较为复杂的环境条件和各种不同的排放方式。这种模型中引入了大量以实验为依据的经验参数，使其计算结果以较强的实践背景为依托，而计算工作量较场模型大大减少，为工程问题混合区流场的分析提供了便利。 3.1.2 智能化软件 伴随着污水海洋处置工程近区模型的日益完善，智能化的软件模型也得到开发和应用，已知的各种智能化软件以美国 开发的为多，其中多半是由美国国家环保局（ USEPA）组织开发的。 1985 年 USEPA 推荐了 5 个污水排海稀释度计算模型（ UPLUME， UOUTPLM， UM，UDKHDEN 和 ULINE）。最早的 UPLUME 模型针对在静水中的圆形单孔排放，基本上只是讲Abraham 的浮射流理论计算机化，它适合于环境水体为任意分层，射流的倾角可以是 -5o90o。UOUTPLM 模型则是基于 Frick 冷却塔热羽流模型，适用于环境流体为均匀流动中的单孔排放，环境流体的密度分层可以是任意的，射流的倾角可以是 -5o90o。 UM 模型则比 UPLUME 模型大进一步，其可以用于单孔圆形浮射流，也可用于多孔排放圆形浮射流，并允许环境流体有流速分布。 UDKHDEN 模型运用欧拉积分方法计算浮羽流的轨迹、半宽、浓度和温度方程，适用于浮羽流的单孔或者多孔排放。 ULINE 模型是针对在流动水体中的线状扩散器的，是对Roberts 在密度不分层的横流中的线源浮力羽流实验数值外推，环境水体密度分层和流速分布可以是任意的。后来 Roberts 又完成了在密度分层横流中的线源浮力羽流实验，将此结果合并入 ULINE 模型形成了一个新的模型，称之为 RSB 模型。 1995 年 USEPA 在颁布这些模型时，又把 UM 和 RSB 模型并入含有远区稀释度计算的 PLUMES 软件。 2001 年， USEPA 改进了PLUMES 模型，发布了 Visual Plume 模型软件，包括 UM3（ UM 的升级版本）、 DKHW（ UDKHDEN 的升级版本）、 PDSW、 NRFIELD（ RSB）、 PDS、 DOS PLUMES 等模型。 此外还有 DKHEN 模型及基于拉格朗日方法的 OUTPLM 模型。 1993 年康奈尔大学更开发了一个专家系统 CORMIX，它把水流的各种水动力条件，如射流与水体表面和水底的动力相互作用都考虑在内，采用特征长度综合 了大量的实验资料，因而该专家系统可以扫描各种可能的初始混合情景。 Lee 开发的 Jetlag 模型 10对环境水流成任意夹角的圆形浮射流进行全场模拟，并可应用于在同一喷头上有多个喷孔的情形。槐文信等也采 用数学模型对射流进行了一系列研究。 3.2 远区模型研究现状 近区研究侧重根据污水的初始动量、浮力、环境水体的分层情况以及水流的作用，远区研究的重点是污水进入海洋水体，污染物经紊动混合与输运扩散后时空分布特性。近岸海域是陆地、海洋和大气之间各种过程最活跃的界面，该水域的环境和生态系统受到来自陆地和海洋的双重作用的影响， 对大范围内各种自然过程变化所引起的波动和人类活动的影响十分敏感，生态系统相当脆弱。因此，近海水域环境污染物的迁移转化机理一直是环境科学的热门课题。污染物经过海流不断的紊动混合与输运扩散过程逐渐形成光滑变化的浓度场，称其为平衡浓度场，平衡浓度场是污染物时空分布特性表征，其理论基础就是污染物的扩散输运模型。近海水域的水质预测、评价和管理中，最为实用的流场模型为 2D 浅水环流模型和二维输运模型。这些模型的应用中，由于控制方程和边界条件的复杂性，其求解过程仍需借助于数值方法，应用最多的如 FDM， FEM， FVM， FAM 等。 Georgetal（ 2003）则将二维水动力模型与生态模型偶合，进行了复合模型研究； Tsanisetal 及张行南等（ 2001， 2004）分别将二维 7 水质模型与 GIS 技术相结合；刘成等（ 2003）应用二维模型分别对尾水排入长江口工程的污染物运动轨迹进行了模拟；娄安刚、王学昌等（ 1993 2006）分别对渤海海峡、胶州湾等海域应用二维数值模型进行水质分析。 近年来三维水质模型得到了应用，如沈永明（ 2000 2004）就将三维模型应用于香港维多利亚港、日本九州岛博多湾等地污染预报中；陈祖军等（ 2003）应用三维 模型对长江口水质模拟研究；林卫强、逄勇等（ 2003）在珠江口海域三维水质模型的应用；闫菊等（ 2001）将三维模型在胶州湾水质模拟中的应用。美国普林斯顿大学三维海洋模式（ POM）及河口陆架海洋模式（ ECOMSI）因国际海洋界的认可而得到广泛应用。国内经过吸收改进，成功将这两个模式用在我国的河口海湾的数值模拟中。另外，还有不少三维模型与生态或富营养模型等混合模型。近年还出现有超标概率场模型、概率分析法模型等模型的应用。 目前文献中常见的三维水质模型系统还有 WASP， CE-QUAL-ICM， EFDC/HEM3D，MIKE3， Delft3D 和 RMA10 等，可实现河流、湖泊、水库、河口和沿海水域等一系列水质问题的模拟 11。 3.3 水动力学参数 排入海域后浓海水浓度的分布除受海洋生物、化学过程影响外，主要受海洋物理过程的控制，其中海洋水动力学是支配海洋中浓海水浓度分布的最主要因子。因此，掌握浓海水排放海域的水动力学状况是了解浓海水浓度分布和变化的关键。 1） 潮汐 评价海域内潮汐类型及潮汐特征值，包括平均海平面、深度基准面、最高高潮高、最低低潮高、平均高潮高、平均低潮高、最大潮差、平均潮差、平均涨潮时间、平均落潮时间。 2） 潮流 日 平均流速、最大流速、最小流速。 3） 波浪 常波向、强浪向、最大波高、平均波高、波浪平均周期。 4） 海水密度 海水平均密度及排放口附近海水密度范围和季节变化。 5） 非周期海流 非周期海流的基本特征和季节变化，包括表层流速、底层流速、海水交换量、近岸区和离岸区的紊流流速和流向。 6） 表面漂流 表面漂流的流速及流向。 3.4 海洋水动力学状况评价方法 为了预测污水在海流驱动下的输运状况，需要更深入地了解排海区域的海水动力学现状。通过建立适合排污海区的潮波模型进行数值模拟，并在此基础上建立水质模型讨论排污后对周围海域的影响，同时还应建立拉 格朗日质点追踪模型，模拟污水标识质点输运的轨迹，从而选择出较佳的污水排放海域。 3.4.1 潮波数值模拟 8 根据排放海域的实际情况，选择合适的潮波数值模型模拟该海域的潮流场，比较模拟和实测的潮位潮流结果，在计算潮流场的基础上，阐明欧拉余流场的特征。潮流场的计算结果可以作为进一步模拟污染物浓度和轨迹的依据。 由于一般拟选海域的水平尺度远大于垂直尺度，海水混合比较强烈，盐度、温度等要素的垂直分布较均匀，可采用非线性二维潮波动力学模型进行海区的潮波数值模拟，其基本方程为： 22221 ( ) ( )abXXU U U A U UU V f V gt X Y X H X Y （ 1） 22221 ( ) ( )abYYV V V A V VU V f V gt X Y X H X Y （ 2） ( ) ( ) 0H U H Ut X Y （ 3） 式中： UV、 分别为 x、 y 方向的平均流速分量； Hh水体的深度， h 为从静水面起算的水深； 自静止水面起算的水位高度； 海水密度； t 时间坐标； f 柯氏参量； aaxY、 海面风应力，如不考虑风影响，则： ax 0， aY 0； A 侧向湍粘性系数； bbxY、 海底摩擦力的 X、 Y 方向分量。采用二次律公式： 1 / 2222( ) ( , )bbxY UVg U VC 、 （ 4） 式中： c chezy 系数，主要取决于海底粗糙度和 水深，一般采用美国水利学界普遍采用的 Manning 公式： 1/ 64.64CHn （ 5） 式中 ： n 表征海底粗糙度的 manning 系数。 3.4.2 浓度预测 根据二维浓度模型公式： ( ) ( ) ( ) ( ) ( )H P H P U H P V P PH D x H D y ft x y x x y y （ 6） 式中： P 垂直平均后的浓海水浓度； Dx Dy、 污染扩散系数在 xy、 方向的分量。按 Elder 的公式： 1 / 2 1( ) 5 . 9 3 ( )gD x D y H C U、 、 V （ 7） 式中： f 污染物质的排放速度或者由于生物、化学、物理等因子引起的衰减速率。 确定定解条件和离散格式，利用不同时刻的流速作 为平流速度场，计算拟建排污口在不 9 同时刻 (包括落潮中间时、低潮时、涨潮中间时、高潮时 )排放的 COD 浓度增量。 3.4.3 浓海水运移轨迹 浓海水通过扩散器进入海水后，污染物大多处于溶解态或悬浮态，它们在海洋中的输移与海水质点一样，其输运途径和去向可以用拉格朗日标识质点轨迹来描述。 假定：1tt时刻在计算海域里投下一标识海水微元，其坐标位置为1()jXt， j 表示标识微元编 号。当21t t t t 时刻，该微元的位置为2()jXt，可表示为： 2121( ) ( ) ( ) ,tj j L jtX t X t U X t t d t （ 8） 式中： ( ),LjU X t t为微元的拉格朗日速度，它是微元的位置和时间的函数。由于潮流场变化的连续性，可以认为 t 足够小时，2t时刻标识微元的拉格朗日流速可用该微元1()jXt位置上展开的欧拉流速的泰勒级数逼近： 2 2 1 2 1 2( ) , ( ) , ( ) ,L j L j L jU X t t U X t t X H U X t t LL（ 9） 若 略 去 2X 以上的项并将 X 近似地取为： 211( ) ,tLjtX U X t t d t （ 10） 则该微元从1nt到nt的位置改变可用下式确定： 111 1 1 1( ) , ( ) , ( ) , ( ) , ( ) , nnttj n j n L j n L j n L j nX t t X t t U X t t U X t t d t H U X t t d t （ 11） 由 此得到一个潮周期的拉格朗日漂移和拉格朗日漂移速度 (或拉格朗日 余流 )； 00( ) ( ) ( )L j j m jX t X t X t（ 12） 00 ( ) ( )() j m jLj X t X tUt T （ 13） 式中： T 一个潮周期， 0mtt、 一个周期的终了和开始的时间。 根据模型计其结果，可以预测浓海水经过一个潮周期运动以后 的实际输送方向和速度，并以此结果划分出海域水交换的活跃区和滞缓区，评价海域的稀释扩散能力。 3.4.4 工程海域稀释扩散模拟实验 为了保证数值模拟结果的可靠性，使排海口选择更科学合理，一般都需要在现场进行示踪剂扩散试验，用现场实测数据校核理论计算结果，可更直观地了解该海区的稀释扩散能力和规律。 目前国外在水中用作示踪剂的荧光染料有 8 9 种之多，国内一般选用国产罗丹明 B 作为示踪剂。结合工况方案，在选择好的投放点以点派形式投放到海水表面，间隔一定时间跟踪取样。通过实验室分析，确定罗丹明 B 浓度与时 间的函数关系，测定罗丹明 B 在不同潮时的运移轨迹，计算稀释因子和海水扩散系数。 1）稀释因子 污染物质进入水体之后，随时间过程被逐渐稀释的 倍数，即称为稀释因子 (diffusion factor，DF)。 mg mgDF 原 始 污 染 物 质 的 浓 度 （ /L ）不 同 时 间 海 水 中 污 染 物 浓 度 （ /L ） （ 19） 10 DF 值愈大，表示污染物浓度愈低。在相同时间内， DF 值愈大，表示稀释愈快，反之亦然。 2）扩散系数 K 对海区的表层平流扩散系数可采用下列公式计算： 2 2e x p44X U t YMKt c t（ 20） 式中： M 投放罗丹明 B 的用量 (g)； t 时间 (s)； c 测定的罗丹明 B 的浓度 (mg/L)。 根据对污染物浓度增量的预测和扩散实验结果的分析，可以初步确定排污口离岸的距离范围，结合进一步的工程设计参数的研究，可以确定出更适宜的距离。 11 4 排海扩散器的设计 浓海水排海工程是浓海水处理和处置的一种常用方法。在浓海水排海中，需要合理利用海洋的稀释和净化容量，否则会造成海洋污染和海洋的生态破坏。扩散器的结构直接影响浓海水排放后在近区的浓海水场的形成，不同的环境保护需求所要求的水质不同，对扩散器结构的要求可能也不一样。以前的工程多采用线源模 型指导扩散器的结构设计，这与实际工程中的导流管型扩散器不相吻合。近年来海洋排放近区模型研究取得很大进展，出现了一些三维模型。运用这些模型来模拟污水排放后的近区浓海水场状况 (如初始稀释度 )，将浓海水排放近区的水质影响同扩散器结构紧密地联系起来。尤其在浓海水海洋处置的环境影响研究或可行性研 究阶段，在扩散器结构设计还没有进行的情况下，采用本论文提出的扩散器结构 概念设计的方法不仅可以评价浓海水排海对环境的水质影响，还可以帮助优化扩散器结构，对实际的设计具有指导作用。 4.1 扩散器结构与要求 4.1.1 扩散器结构的参数 扩散器设计的主要目的是使浓海水水得到足够初始稀释度，同时又投资较省，更充分地利用天然水体的容量。其设计参数包括： 1） 扩散器离岸距离 (放流管长度 ) 放流管占工程造价的很大一部分，放流管越长，造价越高，其决定与环境因子水深有直接关系也是工程选址所需考虑的决定因素。 2） 扩散器长度 3） 扩散器形式及走向 扩散器形式是为了充分地利用海流，同时使潮周内最小稀释度或扩散器的有效长度尽可能的大，保证在潮周内初始稀释度的振荡幅度最小，浓度的变化更均匀，其实质仍是个扩散器长度问题。初始稀释度预测公式及实验均表明，除了静水及流速很小的 情况下，扩散器的走向对初始稀释度相当敏感，直观的判断也很清楚，对于往复流而言，参加稀释的水体总是在扩散器的范围之内，因而，如果环境流向与扩散器正交，则通过扩散器的流量为最大，若为零，则通过扩散器的流量也就最小，一般来说，长的扩散器，首尾两个喷嘴出来的羽流半宽相对于扩散器长度是微不足道。因此，就设计的角度而言，扩散器走向与流向夹角为 90是最佳的选择。但是在环境流速流向不断变化的复杂流动的情况下，扩散器的设计就需要考虑其他类型，而不是总是考虑一种形式。 扩散器形式可分为（见图 1）： I型 即扩散器是放流管的直 线延伸，多适用于具有往复流特性的沿岸流地区。 L型 即扩散器在与放流管末端连接时有一定的夹角，其用意在于满足地形要求及流向要求，使总的管道长度最小，投资最低，而又能使扩散器与环境流向正交。 T型 扩散器的中部与放流管连接呈 90 角，浓海水从扩散器中部分两个方向流入扩散器。这样设计扩散器管径可减少变化。这种 T 型结构主要运用于向岸流及离岸的情况，使扩散器 12 能与流向正交。 Y型 与 T 型相似，扩散器在放流管末端形成两个分支，但与放流管轴线的夹角不是 90 ，一般为 45 ，这种结构运用于旋转流的地区可使扩散器的有效长度 能保持在一定的水平上，而不至于流向的变动，使扩散器有效长度降至或接近于零，影响污水的有效稀释。 图 1 扩散器形式 Fig.1. Diffuser form 4） 喷口高程 5） 喷口角度 6） 喷口直径 其中，前三项的作用最为关键，是主控参数，后几个参数相对比较次要，为辅控参数，一般对韧始稀释度的影响较小。在负浮力射流中，喷射倾角成为主控参数之一。 早期的扩散器只是一根末段开口的管子，而现在的扩散器则沿管道方向在末段设置多个导流管，并在导流管末端设置多个喷口，称为导流管型扩散器。扩散器结构设计的内容在于确定扩散器各部分的 尺寸，将扩散器结构进行简化后 ， 包括：扩散器的管径、导流管根数、每根导流管上的喷口数、喷口直径和扩散器长度。喷口的喷射流量、流速和角度是影响近区初始稀释度和近区污水场的关键因素。这些因素反映到扩散器的结构上与导流管根数、每根导流管上的喷口数和喷口直径有关。扩散器的管径由水力学计算得到，在明确导流管根数和导流管间距后，扩散器的长度也就得到了。因此，考虑影响近区浓海水场水质的扩散器结构主要集中在导流管根数、每根导流管上的喷口数和喷口直径这 3 个因素上。扩散器结构概念设计将讨论在设计的排放流量下，初始稀释度与导流管根 数、每根导流管上的喷口数和喷口直径的关系。 4.1.2 扩散器结构设计的要求 现代排海工程中，对于扩散器的结构有下面几点共识： 1） 为保证扩散器最末端仍能保持一定的喷射流速，同时能量损耗不太大，喷口的平均喷射流速应选择在一定范围内，如 2.03.0m/s12，也有建议选 2.53.5m/s。喷口平均喷射流速的计算公式如下： 24jjQVN n d(21) 式中：jV 喷口的喷射流速， m/s； 13 Q 设计的浓海水排放流量， m3/s； N 导流管根数； n 每根导流管上的喷口数； jd 喷口直径 ， m； 2） 扩散器在运行期间应避免发生海水倒灌现象 13。试验结果表明，当喷口弗劳德 (Froude)数 大于 1.0 时即满足要求。为更安全，有人建议采用更严格的要求，如要求弗劳德数大于 2或 2。弗劳德数的计算公式如下： irjVFgd （ 22） 式中：rF 喷口弗劳德数； jd 喷口直径， m； g 有效重力加速度， 2/ms， 00agg ，a、0分别为排放口处海水和原污水的密度， g 为重力加速度； iV 喷口的喷射流速， m/s。 3） 为防止喷口被阻塞，扩散器 喷口直径应大于 0.055m；也有建议不小于 0.08m14，如果浓海水只进行了预处理，建议的喷口直径不小于 0.15m15。 4） 每根导流管上的喷口数通常设置为偶数，这样便于均匀布置，方便施工。喷口数增加有利于分散出口的污水，但喷口数也不能太多。国外的实验室研究表明，喷口数超过 8 个将由于密集喷射的污水形成环状反而不利于稀释 16。 4.2 扩散器结构概念设计的原理和流程 近区模型用于扩散器结构概念设计的原理是：利用近区模型广泛模拟一定范围内的扩散器结构 (不同导流管数、每根导流管上的喷口数和喷口直径范围的组合 )在排放浓海水时的初始稀释度和近区边缘水质，依据近区边缘水质目标的要求，寻求满足水质要求的合适的扩散器结构范围。近区水质的要求仅仅是限制扩散器结构的一个因素，此外还需要考虑水力学要求 (满足一定的喷射流速的要求 )和工程要求 (防止海水倒灌 )。当然，如果还有其它要求或限制因素，相应也可以考虑进去。 扩散器结构的概念设计流程见图 2。 14 图 2 扩散器结构概念设计流程 Fig.2.The structure concept design process of diffuser 4.3 设计中主要影响因素 影响扩散器设计的因素主要有： 1）排放流量 排放流量是扩散器水力设计中最主要的影响因素，排放流量的变化能带来： 管道内流速的变化，在排放流量达不到设计流量时，会造成悬浮物和颗粒物的沉积； 各喷口出流不均匀。 2）喷口面积与主管面积之比 下游扩散器所有喷口面积与相应主管截面积之比超过 0.5 时，喷口出流量开始变得不均匀。 3）管道坡度及上升管高度 管道坡度对扩散器的影响有两个： 对喷口流 量分布的影响，坡度变大，喷口流量将不均匀； 对海水入侵的影响，在坡度变大时，所需阻止海水入侵的最小流量增加。 4） 管道摩擦系数 摩擦系数增大，扩散器喷口出流量将不均匀。在实际工程中，管道摩擦系数的增大往 往是由于管道内颗粒和悬浮物的沉积造成，由于管道内的沉积，流量变得不均匀。 4.4 孔口出流系数法 4.4.1 孔口出流系数法 随着排海工程技术的发展，扩散器的水力设计方法 (即扩散器水力计算 )已趋成熟主要是采用水力学原理，应用能量方程和流量连续性方程，确定扩散器各部尺寸。 设计海水排放流量 环境水文物理条件 导流管 根数 每根 导流管 上的根数 喷口直径 近区模型 预测初始 稀释度 浓海水水质及近区水质标准 最小出口密度弗劳德数 喷口喷 射流速 满足前面共同要求的合适的扩散器结构 调整扩散器结构尺寸 否 否 是否 防止海水倒灌 是 合适的扩散器结构 之三 满足喷射流速要求 合适的扩散器结构 之二 是否满足水质要求 合适的扩散器结构之一 15 图 3 扩 散器简图 Fig.3.The diagram of diffuser 1)喷口流量方程 在第 j 喷口，污水出流的能量方程为： 2 2 212 2 2j j j a j jjjV P U P Uh h kg g g (23) 式中：1k 局部阻力系数； g 重力加速度。 从上式可以看出：11 / ( 1 ) 2jjU k g E（ 24） 式中： 22j a j jjP P VEg （ 25） 令：j j j jq C a u（ 26） 式中：jC是收缩系数，合并两式得： 同时令：11 / (1 )jjC C k （ 27） 2j j j jq C a g E （ 28） 对于 圆 形喷口的jC经验公式 (Brooks)为： 2 3 / 8 3 / 80 . 9 7 5 1 / 2 0 . 9 7 5 1 j j jC V g E j （ 29） 此公式仅应用于喷口 (或喷嘴 )的收缩系数为 4： 1 或更大及 0.10jjdD的条件下，而且由于jV、jE的不同，jC在各喷口是变化的。 在以上计算中，为获得各喷口的流量值，通常是从扩散器远端开始计算，并假定最远端的条件是已知的，即可假定流量1jQ、jq、jQ；流速1jV、jU和jV，压力1jP和jP及管道高程jh和1jh是已知的。 2）管段的能量方程 对第 J 管段可得如下能量方程： 16 1 jjj j ljPPh h h （ 30） 其中： 22jlj jjVlhfDg（ 31） 式中：jf 达西摩擦系数 在上式中1jh、jh、jP、 是已知的，并且ljh可根据方程（ 31）算出，因此 P 可以求出。 在计算中管道压水头是恒定的，可以得出 ： 2211j j j jV P V Pgg （ 32） 公式（ 32）右边是已知的，将公式（ 25）的右边与公式（ 32）的左边比较，当假设1ajP是已知的，1jE可以得出，但由于 21 /2jVg不知道，式（ 29）中也需要知道1jV来计算1jC。在最初的计算中， 2 /2jVg可以代替式（ 29）中 21 /2jVg，来第一次计算1jC， 则第一次的1jq可以用公式（ 28）计算得到。 利用公式：11j j jQ Q q（ 33） 可以得到第一次1jV的值，1jV的计算可利用流量连续方程： 21 1 1/ ( )4j j jV Q P （ 34） 这样可以第二次得到 21 /2jVg，以上过程可重复计算，直到两次计算获得1jV的值较接近为止。 在以上计算中，可以利用计算机编成程序 方便地完成，直到 Q 等于NQ为止。 4.4.2 孔口出流收缩系数的计算 上已述及，在常用水力设计方法中，对于管壁出流圆形喷口的收缩系数计算经验公式为： 2 3 / 80 . 9 7 5 1 / 2 j j jC V g E （ 35） 除此之外，其它经验公式还有： 1） 管壁出流尖角喷口收缩系数经验计算公式： 20 . 6 3 0 . 5 8 ( / 2 )j j jC V g E （ 36） 2） 立管 单喷口收缩系数近似计算公式： 2 2 44( / 2 ( 1 ( / 2 ) ) )j j j j jjjY V g E V g E rCr （ 37） 4121( ) ( )prj e n r er r cDLfD D C (38) 17 图 4 立管 单喷口示意图 Fig.4. The stand pipe-Single spout schematic 式中：rL 立管长度； rD 立管直径； pD 喷口口径最小处直径； cC 射流收缩系数，等于 /pjUU； r 等于 /prUU； jU 出流速度； en 立管进口局部阻力系数； 1 肘形弯管局部阻力系数； e 收缩段局部阻力系数； rf 立管流程阻力系数。 （ 37）、（ 38）两式 参见图 4。收缩系数cC可以根据自由流线理论导出，它是 /prDD和 的函数， 由下式确定： ( ) / 2 r p ca r c t g D D L （ 39） 式中：cL 收缩段长度。 相应于不同的 /prDD比值和 值的cC值列于表 1。 表 1 射流收缩系数 Cc Table.1.Jet contraction coefficient of Cc Dp/Dr 45 90 135 180 0.0 0.746 0.611 0.537 0.500 0.1 0.747 0.612 0.546 0.513 0.2 0.747 0.616 0.555 0.528 0.3 0.748 0.622 0.566 0.544 0.4 0.749 0.631 0.580 0.564 18 0.5 0.752 0.644 0.599 0.586 0.6 0.758 0.662 0.620 0.613 0.7 0.768 0.687 0.652 0.646 0.8 0.789 0.722 0.698 0.691 0.9 0.829 0.781 0.761 0.760 1.0 1.000 1.000 1.000 1.000 当 / 1 / 4prDD时，en可由下式确定： 2211( / )e n n rB r V U （ 40） 式中： B 为常数，对于锐角进口，取 0.406B ，对于圆角 进口，取 0.1 0.2B ；1 /rr D D；rU为立管内流速。 1取决于肘形弯管的曲率半径。对于图 4 所示的 90 肘形弯管，1 0.80 。 21( 1)e k （ 41） 20 . 0 4 30 . 5 71 . 1 ( / )prDD （ 42） 式中： k 为与 2 有关的系数，可由表 2 选取。 表 2 k 与 2的系数 Table.2. The coefficient of k and 2 2 10 20 40 60 80 100 140 k 0.40 0.25 0.20 0.20 0.30 0.40 0.60 3） 立管 多喷口的jC对于立管 多喷口的出流系数，尚未有可用的近似的方法，只能依靠实验来确定。 4.4.3 孔口出流系数法的局限性及缺陷 从上述分析可看出： 1)孔口出流系数法的关键在于喷口出流收缩系数jC的计算，而且目前只有上述两种 (管壁、立管 单喷口 )计算出流系数的计算经验公式，对 于其它形式的喷口还无法使用，因此，这就限制了该方法应用的范围，使该方法有一定的局限性。 2)在孔口出流收缩系数的计算中，大多为经验公式不能利用现有计算阻力损失的方法，使该方法的应用有一定的缺陷。 4.5 动压力水头法及其优化设计程序 4.5.1 动压力水头法及其优化设计程序 针对孔口出流系数法的局限性及缺陷，目前扩散器水力设计方法主要采用动压水头法。以下将对动压水头法计算公式进行推导，并根据排海管道的特点设计计算优化程序。 19 图 5 排海管道计算示意图 Fig.5.discharge seawater pipe calculating schematic 设1q离岸最远端第一个上升管的浓海水排放量，则对第一个上升管和第一段扩散管，根据能量方程及连续方程可得如下方程： 2 2 1 1 1 11 2 1() ()2 wP V V HZ Z hg （ 43） 2 2 21 1 1 1 1 1 1()4 4 4q V d U d V D （ 44） 2 2 2 21 1 1 1 1 1 11 1 1 1 111()2 2 2 2wl V U H U VhD g g d g g （ 45） 计算1P、1q。 对第 2 个上升管和第 2 段扩散管，根据能量方程及连续方程可得如下方程： 222 1 1 2 22 3 2 1()()2 wP P V UZ Z hg （ 46） 2 2 2 1 / 2122 2 2 2 1 1 2 2 2 2 ( ) (1 ) ( ) 4PHq d g H V （ 47） 212222 1 22 2 2 2 1 221 ()2 2 2wl V VhUD g g g （ 48） 计算2P、2q。 20 式中： 2 2 2( / ) ； 22 2 2( / )D； 21 2 1( / )D D V ； 22 2 2 2 2 2( / )BH 。 类似，对第 i 个上升管和第 i 段扩散管，根据能量方程及连续方程可得如下方程： 22111 ( 1 )()()2i i i i ii i w i iP P V UZ Z hg （ 49） 2 2 2 1 / 212 ( 1 ) 1 2 ( ) (1 ) ( ) 4 iii i i i i i i iPHq d g H V （ 50） ( 1 )222 1( 1 )1 ()2 2 2w i i i i ii i i i i iil V VhUD g g g （ 51） 由此计算iP、iq。 式中： ( / )i i i ； 2( / )i i iD； 211( / )i i iD D V ； 2( / )i i i i i iBH 。 对于高位井或出水泵站，其水位满足： 22010 1 11()n n n T n nT n nnD P V Z Z D VZZZ g D g （ 52） 根据 以上各方程，可以得到计算的排放总量iQq计，若计算结果满足 QQ计 实（允许误差），则计算结束；否则，重新设定1q，重复上述计算。 4.5.2 扩散器水力设计计算优化程序见图 21 图 6 扩散器设计计算优化程序图 Fig.6.The program diagram of the diffuser design calculating optimize 4.5.3 扩散器水力设计计算参数的选取 1）管道沿程阻力 沿程水头损失反映水流为克服摩擦阻力做功而消耗的能量，因而其数值的大小主要取决于管径和管材。根据 谢才公式可以得 到 阻力系数 28/gC ，其中谢才系数 C 可由曼宁公式计算得到，即 1/ 61 ()4DC n， D 为管径， n 为管壁粗糙系数，则单位管道长度的水头损失22ViLDg 。 2） 局部阻力系数 局部 水头损失反映由于局部边界急剧变化，水流在急剧调整过程中所消耗的能量，因而局部阻力系数 的大小主要由管道的局部构造所决定，一般视情况作如下处理： 放流管沿程各转弯处均按急转弯考虑，相应局部阻力系数根据管线沿程转弯角度按经验值选取。 扩散段变径管道按突缩考虑。 上升管与扩散管之间按三通管考虑，其阻力系数按塔里叶夫公式进行计算： 21 1 1 1 1 1 11 ( / ) 2 ( / ) c o s m m m mQ A Q A Q A Q A （ 53） QQ计 实 22 22 2 2 2 2 2 21 ( / ) 2 ( / ) c o s m m m mQ A Q A Q A Q A （ 54） 上升管各喷口可简化为一等效喷口，在计算时按 90 急转弯加突缩情况考虑，即211 . 0 0 . 5 (1 / )AA ，2A、1A分别为喷口截面积和上升管截面积，而对于其他形 式的喷口，可以按等阻力进行校核，对于精确的计算可以从物理模型试验获得。 4.5.4 程序的优点 对比孔口出流系数法及程序本身设计特点，可以得出该程序设计的如下优点： 1) 该程序不仅能确定扩散器各工艺参数，而且可以比较多方案水头损失的大小。 2） 该程序能反映各扩散器设计参数对喷口出流量的影响，如摩擦阻力系数、局部阻力系数、喷口开孔大小、上升管直径、上升管间距等。 3） 该程序在喷口阻力的处理中，是按 90o 弯头逐渐收缩，因此，可以利用现有公式进行计算，计算结果较准确，而对于其它形式的喷口，可以按其喷口的特点，对程序本 身进行很方便的改写，也可按实验或计算得到的阻力相等的原则进行换算，而 不对计算结果产生影响。 4.6 六横岛台门海水淡化工程 某大型海岛油库海水淡化厂位置位于舟山六横岛台门海域，日产淡化海水 1.25 万吨，主要为油库码头渔船供水、油库生产用水、生活用水和消防冷却用水，并向周围制冰和渔业加工等企业供水。经淡化后的浓海水排放海域处于二岛屿之间的狭长水道中，为保护海洋生态环境，要求设计一个高效率大规模浓海水排海扩散器。 4.6.1 选取设计方案 由于 浓海水排放海域处于二岛屿之间的狭长水道中 ，排海扩散器走向较适合采用 T 型 17。因为从水道的利用功能来说，大中型江河多数为航运水道，船只往返频繁，水生生物回游范围要求较大，不允许扩散器布置在水道中央，并且，从施工及运行管理来看，中心布置施工困难，不便于维修及反冲洗，严重影响通航。 T 型对江河的航运、生物回游等影响较小，施工与管理都较为方便，工程投资少，值得推广使用 18。 T 型扩散器就是将扩散器布置成顺水流方向，与输水管道成 T 型，如图 7 所示的排放近区物理形态模型。浓海水射入环境水体的掺混稀释特性及流态的变化规律，取决于 T 型扩散器的几何尺寸、喷嘴射流参数及河道水流条件 19。因此， 排放近区特征参数可表达为 0 1 0 0 0( , , , , , , , , , , , , , )s D a a a af a L h L u C h u C 23 图 7 T 型扩散器排放近区物理形态模型 Fig.7.Flow model of T diffuser discharge in nearfield region 式中： 为扩散器喷嘴排放角； 0a为喷嘴孔口面积； sL为喷嘴孔口间距； 1h为孔口距河床高度； ah为河道水深； 为环境水体运动粘滞系数； 为扩散器与河道水体流动方向的夹角， T 型扩散器 0o 4.6.2 T 型扩散器排放近区的数值模拟 为了预报 T 型扩散器排放近区的特性，运用本文所取的扩散器在宽浅江河中排放的流 场和浓度数值模型来模拟计算 T 型扩散器排放近区。要进行计算，首先将 T 型扩散器排放口由点排放概化到所在网格单位时间单位体积的质量源，其表达式为0 0 0 / ( )sau a L h y ，则喷嘴所在网格的连续方程为： 0 0 0 / ( )sah h u h u h u a L h yt x y （ 55） 式 中 ： y 为网格尺度 ； 水深ahh。 喷嘴所在网 格的动量源 0xm ， 20 0 0 / ( )y s am a u L h y， 则得到 X 方向动量方程 ： () 2 ( ) ( ) b b xe e xg h Z hh u h u h u u u u hu v v h v h mt x y x x x y