污染物在河流中的混合PPT文档.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,河流的水质好坏对直接工农业生产和人民的生活。,自环境问题出现以来，人们对河流污染的预测和防治进行了大量的研究，已取得了很多成果。,第五章 污染物在河流中的混合,混合：是指污水进入环境水体之后的混掺和输移的过程,本章将对河流在稳态和动态情况下初始段和远区的浓度计算问题进行介绍，其中对污染带的计算将给出较详细的分析和论述。,1,污水排入河流之后的混合过程，可以将其划分为三个阶段：,图,5-1,污水与河流的混合过程,第一节 河流中的混合过程,2,（,1,）第一阶段,(,垂直混合阶段，也称为初始稀释阶段,),是从排污口开始到污水在垂直方向完全混合为止。该阶段实际上是一种三维混合过程，只是由于水深（垂向）的尺度比其他两个方向的尺度要小的多，所以首先完成垂向混合。,混合情况与污水排出时的初始动量和浮力以及排污的位置等有关：,如果污水排出的流速大于河水流速，则低流速的河水会卷吸到高流速的污水之中，从而加强了污水的初始稀释。如果排出的流速较小，可以不考虑这种卷吸作用。,如果污水的密度比河水要小，就有浮力作用，例如热电厂的冷却水（水温较河水高）要考虑浮力作用。如果污水的密度比河水大，就有下沉作用。如果两者密度相差很小，就不考虑浮力和下沉的影响。,第一节 河流中的混合过程,3,排污口的位置有表面排放与淹没排放两类。如果污水在水下较深处排放，则可利用较大的水深使污水在河流中达到较好的初始稀释。,在初始稀释过程中，射流的动量和浮力的作用也将随之减弱，在第二阶段就不考虑其影响。,对第一阶段的计算需要用浮力射流理论等有关知识。,第一节 河流中的混合过程,4,（,2,）第二阶段,(,横向混合段或初始段,),从污水在垂直均匀混合之后算起至河流横向（在断面上）均匀混合为止。,在本阶段中，初始动量和浮力已经消失，混合取决于河流中的二次环流和横向紊动的作用。在此过程中，横向的污染范围逐渐变宽。如果污水的出流是恒定的时间连续源（即稳态情形）在本段将形成一条稳定的污染带。横向混合的结果导致达到全断面的均匀混合。,为了简化分析，可以对该阶段的流速，浓度和横向混合系数都各自沿水深平均，只研究垂线上平均值的纵向和横向变化，按水平二维的混合过程处理，可以应用二维紊流扩散方程作为控制方程和进行计算。,第一节 河流中的混合过程,5,（,3,）第三阶段,(,纵向分散段,),从河流横断面均匀混合以后起算的阶段。,在本阶段中，在横断面上的浓度分布是均匀的，服从一维纵向分散方程，同时必须考虑污染物质的非保守性。,图,5-2,污水与河流的混合过程,第一节 河流中的混合过程,6,第一阶段在排污口附近，称为近区。一般是三维问题，需要浮力射流理论。,第二、第三阶段发生在离排污口较远的区域，称为远区。,顺直河流断面完全混合时的距离（河长）：,中心排污：,L=0.1VB,2,/M,y,岸边排污：,L=0.4VB,2,/M,y,M,y,为河流横向混合扩散系数。,第一节 河流中的混合过程,7,对三个阶段的划分也不是严格的，因为各个方向的混合并不是截然分开的，这样的划分只是反映了混合过程各个时期的主要特征，况且在实际问题中，也不一定都按三个阶段进行处理：,（,1,）污染物质的非射流排放，第一阶段的距离可能很短，可,忽略不计；,（,2,）如果河流的宽度比深度大的多，垂向混合与横向混合相,比可认为是瞬时完成，也可忽略第一阶段；,（,3,）由于河流不太宽，而射流排放的初始动量很大，垂向混,合段很长，污染物质在第一阶段就扩展至全断面，所以,不存在第二阶段；,（,4,）如果进行水质规划，从大范围来看河流的混合，相对说,来，第一和第二阶段很短，第三阶段才是主要的，此时,也可以忽略第一和第二阶段。,第一节 河流中的混合过程,8,第二节 矩形河道均匀流污染带的计算,就计算方法而言，有确定性方法和随机方法两类：,（,1,）确定性方法以紊流扩散为控制方程，对浓度等问题进行求解；,（,2,）随机方法从扩散位移是随机过程的观点出发，采用概率论的数学方法处理。,本章主要介绍在矩形河道均匀流中和在不规则河道渐变流中污染带计算的确定性方法。,9,在大多数河流中，河宽远大于水深。例如：有一河流宽,W,=30m,、深,h,=1m,，初步用下式估算垂向紊动扩散系数,E,z,和横向混合系数,M,y,：,通常认为垂向混合相对于横向混合来说是瞬时完成的,由量纲分析可知，混合时间,t,L,2,/E,i,（,L,表示某一特征长度），故有：,M,y,10,E,z,第二节 矩形河道均匀流污染带的计算,10,在实际应用中，为了研究方便，可以假定污染物质在开始时就是沿垂线均匀混合的（即忽略第一阶段），或者说开始时就可以作为一条垂直均匀混合线源来分析其水平二维扩散问题。如果第二阶段的距离不是太长，此时可以忽略污染物质的非守恒性，作为示踪物质考虑。,第二节 矩形河道均匀流污染带的计算,11,一、污染带的浓度,当河槽近似为矩形棱柱体，水流近似为均匀流，水深和断面平均流速分别为,h,和,V,，假设：,断面上所有点流速,u,V,u,=,w,=0;,污染源为时间连续沿水深的线源,单位时间内沿水深方向上注入的污染物质质量为,(,量纲为,MT,-1,L,-1,);,不考虑岸边对横向扩散的反射作用和污染物质的非守恒性；,浓度和横向混合系数沿水深平均，研究浓度的垂线平均值在纵向和横向上的变化。,第二节 矩形河道均匀流污染带的计算,12,（,5-2-1,）,式中：,c,和,M,y,均为沿水深的平均值,，在不致引起混乱的情况下,省去在字母两侧的代表沿水深取平均的两条竖线。,在上述假设下，求污染带的浓度问题便简化为求一维纵向随流一维横向紊动扩散的稳态解，其控制方程：,随流紊动扩散方程,第二节 矩形河道均匀流污染带的计算,13,当坐标原点与污染原点重合时,参照,连续无限长恒定线源一维随流一维横向紊动扩散的稳态情形的解的形式：,（,5-2-2a,）,可得式,(5,-2-1,),的解为,:,第二节 矩形河道均匀流污染带的计算,14,将式（,5-2-2a,）改写为：,（,5-2-2a,）,令 ，便有,（,5-2-2b,）,式中：,c,m,的意义为污水与河水完全（均匀）混合后的浓度；,Q,为河流流量；,Q,d,和,c,d,分别为从排污口注入河流的污水流量和浓度。,第二节 矩形河道均匀流污染带的计算,15,图,5-3,污染源点位于,y,0,或,（,5-2-3a,）,（,5-2-3b,）,第二节 矩形河道均匀流污染带的计算,16,图,5-3,河流污染带的起始计算断面和坐标,进一步设河槽的宽度为,W,，,坐标原点取在左岸水边，污染源位于,x,=0,、,y,=,y,0,处。假设两岸边界为完全反射，则可在解式（,5-2-2a,）的基础上用像源法解决。,第二节 矩形河道均匀流污染带的计算,17,对真源,y,=,y,0,，在,y,=0,和,y,=,W,均有完全反射壁：,2W+y,0,W,y,0,0,-y,0,-2W+y,0,2W-y,0,图,5-4,两岸反射的像源法,第二节 矩形河道均匀流污染带的计算,18,考虑到两岸的反射，利用像源法便得污染带的浓度解：,无边壁反射：,2W+y,0,W,y,0,0,-y,0,-2W+y,0,2W-y,0,图,5-4,两岸反射的像源法,第二节 矩形河道均匀流污染带的计算,19,（,5-2-4a,）,（,5-2-4b,）,或,式中：,n,取整数,在实际应用中，一般只取,n,=0,，,+1,，,-1,计算就足够准确了。,第二节 矩形河道均匀流污染带的计算,20,二、污染带的长度（带长）,定义：从污染源的断面（或从垂向扩散完成的断面）开始至完全混合,(,c=c,m,),的断面为止的一段纵向距离。,断面上各点的浓度,c,均满足,|(,c-c,m,),/,c,m,|,5%,，可近似认为此时已达到完全混合。,第二节 矩形河道均匀流污染带的计算,图,5-1,中心排放,21,设污染源点位于河中心线上（,y,0,=,1/2,），据式,（,5-2-4b,）,:,分别算出沿中心线（,y,=,1/2,）和沿岸边线（,y,=0,或,y,=1,）的相对浓度值,c/c,m,。,第二节 矩形河道均匀流污染带的计算,从上式得不到带长的显式解析解，费希尔提出了用数值解来求带长的方法。,22,图,5-5,中心排放时沿中线和岸边的浓度曲线,通过计算，当,x,0.1,，,断面上各点的浓度,c,均满足,:,|(,c-c,m,)/,c,m,|5%,，可近似认为此时已达到完全混合。,当中心排放时，由,x,=,0.1,（忽略第一阶段的长度），有带长,L,p,的近似式：,（,5-2-5,）,第二节 矩形河道均匀流污染带的计算,23,如果水流条件和边界条件不变，但将中心排放改为在岸边一侧排放。此时，岸边排放的污染带形状与中心排放的污染带的一半是相似的。也就是说，岸边排放具有的横向扩散宽度是中心排放的一侧宽度的两倍。,图,5-1,中心排放,图,5-2,岸边排放,第二节 矩形河道均匀流污染带的计算,24,上述两个带长公式并,没有得到实验的支持，主要是因为浓度沿纵向变化很慢，对完全混合的标准也难以掌握，以致对带长的量测有很大的不确定性。目前，对带长公式的实验验证仍然有困难。,（,5-2-6,）,岸边排放时的带长近似式：,以,2,W,代替,W,第二节 矩形河道均匀流污染带的计算,25,（,5-2-7,）,式中：,K,为带长系数，视点源或线源以及源的位置而定。,对点源情形，有：,应用最大熵原理，从理论上导出满足带长定义的带长公式：,（,5-2-8,）,在特殊情形下：当中心排放,y,0,/,W,=1/2,时有,K,=1/24,；,当岸边排放,y,0,/,W,=0,或,1,时有,K,=1/6,。,点源位置,第二节 矩形河道均匀流污染带的计算,26,对横向线源情形，近似有：,式中,:,y,01,和,y,02,分别为线源的始点和终点的横坐标，且规定,y,01,和,y,02,必须同在河中心线的一侧，即,y,01,、,y,02,W,/2,。,（,5-2-9,）,对点源情形,y,01,=,y,02,：,线源,y,01,y,02,第二节 矩形河道均匀流污染带的计算,27,三、污染带的宽度（带宽）,带宽是指污染带的横向宽度,图,5-1,中心排放,图,5-2,岸边排放,从理论上说，只要根据带边的浓度应为零的要求就可以求带边的,y,值,从而求出带宽。,由于恒定线源一维随流一维横向紊动扩散的稳态解是正态型，当,y,，才有,c,0,。况且，也不能将稍有一点污染物就说成受污染了。,必须对污染带的污染含义给出数值性的规定。,第二节 矩形河道均匀流污染带的计算,28,三、污染带的宽度（带宽）,常用的带宽定义有两种,第一种：认为带边的浓度为同断面上最大浓度的,5%,。据此规定，便可求出相应于各个,x,值的带边,W,p,值。,当中心排放时，各断面的中心点就是出现最大浓度的点；,当岸边排放时，各断面的最大浓度点均位于排放岸边处。,第二节 矩形河道均匀流污染带的计算,29,第二种：由于浓度在横向上服从正态分布，宽度为,4,y,的正态分布曲线下的面积占总面积的,95.4%,，习惯上取,4,y,的宽度来代表正态曲线的宽度。,对岸边排放时的带宽,:,（,5-2-10,）,（,5-2-11,）,对中心排放，如果污染物质尚未扩展到岸边，此时取带宽为,4,y,，便得中心排放时的带宽,:,图,5-6,一维随流一维横向紊动扩散的稳态解,第二节 矩形河道均匀流污染带的计算,30,四、时间连续横向线源的污染带,如果排入河流的污水具有足够的初始动量和浮力，则垂向混合阶段会较长，可能不宜忽略第一阶段。在第一阶段的末断面上，已有一部分达到横向完全混合，设其浓度为,c,i,(,y,),，现在将它看作是已知的时间连续横向线源来计算下游的污染带。,线源,y,01,y,02,第二节 矩形河道均匀流污染带的计算,31,（,5-2-12,）,时间连续点源的污染带浓度解式：,式中：,y,01,和,y,02,分别为横向线源两端点的横向无量纲坐标,。,对上式按横向线源积分，得时间连续横向线源的污染带浓度,:,第二节 矩形河道均匀流污染带的计算,32,如果横向线源的浓度是均匀的,亦即,c,i,(,y,),=,c,0,(,常数,),有：,令,上式变为,:,第二节 矩形河道均匀流污染带的计算,33,（,5-2-13a,）,最后积分得,:,关于带宽和带长,可据上式并参照点源的做法及式,(5-2-9),或,（,5-2-13b,）,第二节 矩形河道均匀流污染带的计算,34,（,5-2-14,）,式,（,5-2-13,）,也可用于污染源是来自一条横置于水流中的污水扩散器（该扩散器是一条水管，在管壁沿管轴连续装有很多喷嘴使污水射入河水中），其中假设射入河中的污水瞬时完全垂直混合，混合后的横向线源浓度为：,式中,：,(,y,01,-,y,02,),为扩散器的无量纲长度,(,y,01,=,y,01,/,W,y,02,=,y,02,/,W,),；,c,d,和,Q,d,分别为污水的浓度和流量；,Q,为河流流量。,第二节 矩形河道均匀流污染带的计算,35,例：某一河流，有微弯，边岸和断面的变化不大，流量为,141m,3,/s,，近似为均匀流，河宽为,124m,，水深为,1.86m,，河床糙率为,0.025,。有工业污水排入河中，污水流量为,0.132m,3,/s,，含有害的守恒物质，浓度为,200mg/L,。,（,1,）设排污口位于河中心，求下游,1km,断面上的最大浓度及该处的带宽,;,（,2,）设排污口位于左岸边，求下游,1km,断面上的最大浓度及分别离左岸,31m,和,62m(,河中心,),处的浓度，并求该断面的带宽,;,（,3,）分别求上述中心排放和岸边排放的带长。,第二节 矩形河道均匀流污染带的计算,36,解,:,断面平均流速：,剪切流速：,按题给的河流情况,横向混合系数取,:,单位时间内在单位水深上注入的污染物质质量为,:,第二节 矩形河道均匀流污染带的计算,37,（,1,）求某指定点的浓度,当中心排放和岸边排放时，断面上的最大浓度分别出现在该断面上的河中心和排放岸的岸边上。为便于计算，令,第二节 矩形河道均匀流污染带的计算,38,表 污染带浓度计算,污染源位置,C,k,/mgL,-1,y,0,/m,x/m,y/m,n,E,1,E,2,c,/mgL,-1,中,心,排,放,0.74,62,1000,62,(,河中心,),0,1,0,0.74,+1,0,0,-1,0,0,岸,边,排,放,0.74,0,1000,0,(,左岸边,),0,1,1,1.48,+1,0,0,-1,0,0,0.74,0,1000,31,0,0.047,0.047,0.07,+1,0,0,-1,0,0,0.74,0,1000,62,(,河中心,),0,0,0,0,+1,0,0,-1,0,0,第二节 矩形河道均匀流污染带的计算,39,（,2,）求带宽,方法,1,：根据带边的浓度为同断面上最大浓度的,5%,的要求。,中心排放：,y,0,=62m,岸边排放：,y,0,=0,第二节 矩形河道均匀流污染带的计算,40,（,2,）求带宽,用试算法求解上式中,y,的值。,第一种方法：根据带边的浓度为同断面上最大浓度的,5%,的要求，有：,（,5-2-15,）,第二节 矩形河道均匀流污染带的计算,41,对岸边排放：,仍保留坐标轴原点在起始断面的左岸边，便有,y,0,=,0,。,对,x,=1000m,，用试算法得：当,y,=34m,，,f,(,y,)=0.051,，于是有带宽,W,p,=,y,34m,。,在同断面上，岸边排放的带宽比中心排放的带宽要小,50%,。,对中心排放：,为使计算简化，将坐标原点改置于起始断面的河中心，便有,y,0,=,0,。,对,x,=1000m,，用试算法：分别设,y,=40,、,35,、,34m,，依次得,f,（,y,）,=0.012,、,0.041,、,0.051,，故认为,y,=34m,，于是有带宽,W,p,=2,y,68m,。,第二节 矩形河道均匀流污染带的计算,42,第二种方法,对岸边排放，由式（,5-2-11,）计算，得带宽,:,对中心排放，由式（,5-2-10,）计算，得带宽：,可见，第二种方法比第一种方法算得的带宽要小。,第二节 矩形河道均匀流污染带的计算,43,（,3,）求带长,对岸边排放，带长,:,对中心排放，由式 计算，得带长：,对只有,124m,宽的河流来说，这两个那么长的带长会令人不敢相信，改用式（,5-2-7,）和式（,5-2-8,）,计算：,对岸边排放,对中心排放,第二节 矩形河道均匀流污染带的计算,44,例：两条河流汇合时的混合问题在实际中是常见的。例如某城市的水源，一部分来自上游的一条河，另一部分来自本地区的集雨，两者水质不同，在进行水处理之前，先将两者混合。设每个水源流量均为,1.42m,3,/s,，现设计一条宽为,6.10m,的矩形渠道，纵坡,S,为,0.001,，糙率,n,为,0.030,，以便两个水源在该渠汇合后达到完全混合（见图,5-7,），试问该渠道要多长？,图,5-7,两个水源（或两条支流）的汇合,第二节 矩形河道均匀流污染带的计算,45,解 先求出该设计渠道的水深,h,和断面平均流速,V,：,已知连续方程,Q,=,hWV,，式中的断面平均流速由谢才公式和满宁糙率确定，于是有,（,5-2-16,）,式中：,Q,为两个水源渠合后的流量,，为,2.84m,3,/s,。,用试算法解得上式的,h,=0.670,（,m,）,由,V,=,Q,/(,hW,),得,V,=0.695,（,m/s,）,由 得剪切流速,u,=8.11,10,-2,（,m/s,）,采用费希公式,M,y,=0.15,hu,得：,横向混合系数,M,y,=8.16,10,-3,（,m,2,/s,）,第二节 矩形河道均匀流污染带的计算,46,带长系数按式（,5-2-9,）计算，得：,代入,带长计算公式得带长：,如果,采用带末端断面上各点浓度均达到,|(,c-c,m,)/,c,m,|5%,的要求来进行计算,由数值解可得,x,=,0.35,便得带长,L,p,=1109m,，该值比,264m,大,3.2,倍,。,在渠的上游断面,(,x,=0),有横向线源，线源起端为,y,01,=0,，末端,y,02,=,W,/2,，设线源浓度为,c,0,，则：,在渠下游达到完全混合时的浓度为,c,m,=,c,0,(,Q,/2)/,Q,=,c,0,/2,第二节 矩形河道均匀流污染带的计算,47,第三节 不规则河道非均匀流污染带的计算,累积流量法,不规则河道,主要是指河道断面形状沿纵向变化大，河道具有较大的弯曲。,不规则河道水流的非均匀性显著，水深,h,是,x,、,y,的函数，流速,u,是,x,、,y,、,z,的函数。,分析水平二维扩散时，将流速和浓度按水深平均考虑，虽然可以认为垂向流速,w,=0,，并假设纵向流速,u,沿水深不变，只是,x,、,y,的函数,但不能忽略横向流速,v,，它沿水深的平均也是,x,、,y,的函数，它在不规则河道的流动和横向随流扩散中起重要作用。,u(x,y),v(x,y),图 不规则河道示意图,48,自然坐标系比直角坐标系更适合于二维水质模型的解析解法以及河道基本均直、水文水力学条件沿河流纵向和横向分布均匀变化的情况,What,？,自然坐标系,49,自然坐标系,如果质点做,平面曲线运动,，且被,约束在已知的轨道,上，则可采用“自然坐标系”。,所谓“自然”，意即“顺其自然”，把轨道当作坐标的“轴”。,质点的坐标是代表路程。,质点的运动方向规定为轨道切线的正方向,。,50,自然坐标系的建立,将此轨道曲线作为一维坐标的轴线，在其上,任意选一点,O,作为坐标原点,。,质点在轨道上的位置可以用从原点,O,算起的,弧长度,s,来表示,，,s,称为,弧坐标,。,运动方程：,自然坐标系是建立在物体运动的轨迹上的。,切向坐标,沿运动轨迹的切线并指向质点运动的方向；,法向坐标,n,沿运动轨迹的法线方向并指向曲线凹侧。,在质点上建立两个的坐标轴：,切向坐标,和,法向坐标,。,51,强调：,自然坐标系是建立在运动质点上的，它随质点一起运动在轨道曲线上。轨道上各点的自然坐标系的二个坐标轴的方位是不断变化的。,为单位矢量，,大小,不变（模为,1,不变），但,方向,随时间改变。,52,建立一个平面的自然坐标系：,x,坐标与流线重合，,y,坐标与,x,坐标（流线）垂直。作为,x,轴的那条流线把河流流量分为一半，沿,x,轴的各分段长度,x,彼此相等，沿,y,轴的各分段长度,y,也彼此相等。纵向坐标线都是流线，横向坐标线都是过水面线，它们处处与纵向坐标线垂直。,图 不规则河道的自然坐标系,第三节 不规则河道非均与流污染带的计算,53,设垂向流速 ，由于纵向紊动扩散项比同方向的随流项小得多，可将纵向紊动扩散项忽略，则上式又简化为：,将上式沿水深取积分，并以符号,代表，即,（,5-3-2,）,=,于是有,（,5-3-3,）,第三节 不规则河道非均与流污染带的计算,一、基本方程,对稳态情形（），三维随流紊动扩散方程变为：,（,5-3-1,）,54,式中：各个沿水深平均值为,（,5-3-4,）,令,（,5-3-5,）,和 分别为沿水深的浓度偏离系数和沿水深的横向流,速偏离系数。,第三节 不规则河道非均与流污染带的计算,55,显然有,（,5-3-6,）,利用式 和式 有,（,5-3-7,）,第三节 不规则河道非均与流污染带的计算,56,（,5-3-8,）,式中：,（,5-3-9,）,第三节 不规则河道非均与流污染带的计算,57,（,5-3-8,）,式中：,如果在河底和水面没有污染物质的源和汇，则：,（,5-3-9,）,即：,式（,5-3-8,）中还有,:,第三节 不规则河道非均与流污染带的计算,58,根据横向流速沿水深分布的非均匀性而产生分散的概念，定义,式中：,M,y,为横向混合系数，主要反映横向二次流的分散作用将上式代入式（,5-3-10,），有,（,5-3-10,）,（,5-3-11,）,（,5-3-12,）,第三节 不规则河道非均与流污染带的计算,59,将横向紊动扩散项取消,第三节 不规则河道非均与流污染带的计算,60,u,和,v,沿水深平均后的二维水流连续方程为：,（,5-3-15,）,则式 变为：,（,5-3-16,）,是作为累积流量法基础的方程式,第三节 不规则河道非均与流污染带的计算,为了书写简化，取消式中表示沿水深平均的双竖线，得,（,5-3-13,）,（,5-3-14,）,61,式中：,m,x,和,m,y,分别为坐标系的度量系数。,m,x,等于沿纵向坐标线量度的距离与在,x,轴上量度的距离之比,m,x,等于沿横向坐标线量度的距离与在,y,轴上量度的距离之比,x,轴线上的,m,x,=1,，,y,轴上的,m,y,=1,由于各条纵向坐标线上的,x,一般都分别与,x,轴上相应的,x,的长度不相等，各条横向坐标线上的,y,一般也分别与,y,轴上相应的,y,的长度不相等，故在自然坐标系情况下，式（,5-3-16,）应写为,（,5-3-17,）,第三节 不规则河道非均与流污染带的计算,62,同样，水流连续方程应改写为,（,5-3-18,）,只有当流线剧烈弯曲，才有必要引入度量系数。,在一般情形下，从实用上说是不需要的。,在,自然坐标系,情况下，二维随流一维横向扩散方程应写为,（,5-3-17,）,第三节 不规则河道非均与流污染带的计算,63,已知单宽流量可表为,q,=,hu,，引入累积流量,q,c,坐标，建立,q,c,与,y,的关系：,式中,m,y,只是,y,的函数，而与,x,无关。,（,5-3-19,）,流管内的流量沿程不变（,q,c,/,x,=0,）,由式,得：,0,第三节 不规则河道非均与流污染带的计算,64,q,c,/,y,=,m,y,hu,为方便计算，有时将上式改写为,进一步将累积流量坐标,q,c,加以无量纲化，令,式中：,h,(,y,),称为无量纲累积流量坐标，,Q,为河流流量。,（,5-3-21a,）,（,5-3-21b,）,（,5-3-22,）,将式 通过复合函数表示为,（,5-3-20,）,第三节 不规则河道非均与流污染带的计算,65,令,称,D,y,(x,h,),为横向扩散因素，量纲为,L,5,T,-2,。,（,5-3-23,）,（,5-3-24,）,D,y,综合地反映了水流的流动，包括纵向流速、紊动、二次环流以及水深、河床地形、河道弯曲等因素对横向扩散输移的影响。,累积流量坐标下计算污染带的基本方程,第三节 不规则河道非均与流污染带的计算,66,累积流量法的精髓：将上两式比较可知，由于累积流量的坐标线与流线重合，使得,h,u,=0,，在累积流量坐标下计算污染带的基本方程中便不会出现横向流速,u,。这样，既通过使用累积流量坐标对,u,加以考虑，而同时又在计算中避开了当,u,出现时的麻烦，这样就将河槽矩形化了，使计算大为简化。,第三节 不规则河道非均与流污染带的计算,67,q,c,(,y,),值可由实测流速、流量、水深和断面形状等资料求出。,如果没有流速资料，而只有流量、水深和断面形状资料，则可用,Sium,提出的求单宽流量的经验公式,（,5-3-25,）,式中：为平均单宽流量，,=,Q,/,W,W,为水面宽；,为断面平均水深，,=,A,/,W,A,为过水断面面积,;,b,o,和,b,1,为经验参数，当,b,o,=,1,、,b,1,=,5/3,便与曼宁公式一致,对顺直河道：当,50,W,/,70,b,o,=,0.92,，,b,1,=,7/4,二、无量纲累积流量坐标的建立,第三节 不规则河道非均与流污染带的计算,68,式中：,y,=,y,1,+,y,2,+,y,j,。,为了使式（,5-3-26,）满足,(,y=W,)=1,，在计算中可能要对经验参数值作一些修正。,有,:,（,5-3-26,）,将式,代入式,:,h,i,h,i+1,j+1,j,j,第三节 不规则河道非均与流污染带的计算,69,例：有一不规则的弯曲河段，在断面,A,上，水面宽,W,为,300m,，过水断面面积为,378.75m,2,，,y,轴原点在左岸水边，取,y,均为,25m,，每个,y,的中点水深见表。初取,b,o,=1,，,b,1,=5/3,，试求无量纲累积流量坐标,与,y,的数值关系。,解：,第三节 不规则河道非均与流污染带的计算,断面平均水深,70,1.00,0.993,0.968,0.921,0.864,0.780,0.679,0.551,0.388,0.224,0.090,0.016,0,7,488,483,471,448,720,380,330,268,189,109,43.7,7.87,0,6,0.217,0.488,0.918,1.102,1.629,1.966,2.493,3.175,3.175,2.614,1.431,0.015,0,5,0.40,0.6,0.95,1.06,1.34,1.50,1.73,2.00,2.00,1.78,1.24,0.50,0,h,i,/m,4,300,275,250,225,200,175,150,125,100,75,50,25,0,y/m,3,25,25,25,25,25,25,25,25,25,25,25,25,0,2,12,11,10,9,8,7,6,5,4,3,2,1,0,序号,i,1,表,断面,A,的,h,(,y,),计算,该结果不等于,1,不合理的,y,i,/m,第三节 不规则河道非均与流污染带的计算,71,现在修改,b,o,的取值，以使,h,(,y,=,W,),=1,，由,最后由式,得表中第,7,行,h,(,y,）的,值。,得：,第三节 不规则河道非均与流污染带的计算,72,1.00,0.993,0.968,0.921,0.864,0.780,0.679,0.551,0.388,0.224,0.090,0.016,0,7,488,483,471,448,720,380,330,268,189,109,43.7,7.87,0,6,0.217,0.488,0.918,1.102,1.629,1.966,2.493,3.175,3.175,2.614,1.431,0.015,0,5,0.40,0.6,0.95,1.06,1.34,1.50,1.73,2.00,2.00,1.78,1.24,0.50,0,h,i,/m,4,300,275,250,225,200,175,150,125,100,75,50,25,0,y/m,3,25,25,25,25,25,25,25,25,25,25,25,25,0,2,12,11,10,9,8,7,6,5,4,3,2,1,0,序号,i,1,表 断面,A,的,h,(,y,),计算,y,i,/m,第三节 不规则河道非均与流污染带的计算,73,M,y,反映了二次环流和紊动对横向扩散的影响。,如果河段中断面形状沿程变化较大，则可将该河段再分为若干段，每段都有它的代表断面，此时对,M,y,值取分段平均。,如果断面形状沿程变化不大，此时对,M,y,值可取全长平均。（,m,x,h,2,u,）反映了河道弯曲、水深和纵向流速对横向扩散的影响，（,m,x,h,2,u,）随点位置（,x,，,y,）而变，可视情况取断面平均或取全长平均。,Lau,等人认为，对,M,y,取分段平均，但保持（,m,x,h,2,u,）是,x,和,y,的函数，这样计算效果较好。,三、横向扩散因素的处理及最简化的基本方程,第三节 不规则河道非均与流污染带的计算,74,有,令,式中：,D,也称为横向扩散因素，量纲为,L,-1,。,对较简单的边界条件，上式可能有解析解。,（,5-3-28,）,（,5-3-29,）,如果对,M,y,和（,m,x,h,2,u,）都取全长平均，则,D,y,是常数：,（,5-3-27,）,第三节 不规则河道非均与流污染带的计算,75,M,y,/,(,VW,2,),D,第四节 累积流量法的污染带方程的解析解,稳态情况下一维纵向随流一维横向紊动扩散的控制方程可改写为：,只要,D,为常数，就可以将,稳态情况下一维纵向随流一维横向紊动扩散的方程,的定解和由该式引导出来的各种定解，套用到不规则河道非均匀流污染带的相应定解问题中来。,稳态情况累积流量坐标下污染带计算的控制方程为：,y,(,y,/,W,),76,1,、时间连续点源，未受岸壁反射,参照（,5-2-2b,），并注意 ，有解,2,、时间连续点源，考虑两岸反射,参照式（,5-2-4a,），有解,（,5-4-1,）,（,5-4-2,）,第三节 不规则河道非均与流污染带的计算,77,3,、时间连续线源，考虑两岸反射,参照式（,5-2-13a,），有解,（,5-4-3,）,以上各式中：,c,o,为已知常数，对扩散器则：,第三节 不规则河道非均与流污染带的计算,78,例：考虑前面例子给出的河段，在该例中已对水面宽为,300m,的断面,A,，求出,(,y,),了。现补充下列数据：,河水流量为,162m,3,/s,，污水在断面,A,离左岸,25m,处有一排污口，污水流量为,0.20m,3,/s,，浓度为,140mg/L,；该河段的平均水深为,1.22m,，平均水面宽为,312m,，断面平均流速为,0.42m/s,，河床糙率为,0.025,；在断面下游,500m,处的断面,B,上，水面宽为,312m,，取,y=,26m,，每个,y,的中点水深见表；平均纵向度量系数均为,1.0,，取无量纲横向混合系数,M,y,/(,u,*,)=0.6,。试求断面,B,上左岸水边和离左岸,25m,处的浓度。,y=25m,x=500m,W,a,=300m,Q=162m,3,/s,A,B,W,b,=312m,第三节 不规则河道非均与流污染带的计算,79,1,序号,i,0,1,2,3,4,5,6,7,8,9,10,11,12,2,0,26,26,26,26,26,26,26,26,26,26,26,26,3,0,26,52,78,104,130,156,182,208,234,260,286,312,4,0,0.38,0.63,0.93,1.00,1.30,1.45,1.68,1.95,1.95,1.70,1.20,0.45,5,0,0.199,0.463,0.866,1.00,1.55,1.86,2.37,3.04,3.04,2.42,1.36,0.264,6,0,5.18,17.2,40.3,66.3,1.7,155,217,296,375,438,473,480,7,0,0.011,0.036,0.084,0.138,0.222,0.323,0.452,0.617,0.783,0.914,0.988,1.00,表 断面,B,的,（,y,）计算,解：在断面,A,上，,y,=25,相应的源点位置为,o,=0.016,。,建立断面,B,的,(,y,),：暂取,b,o,=1,、,b,1,=5/3,，计算结果见表,其中,b,o,修正为：,在断面,B,上有：,y,=0,和,y,=25,分别相应有,=0,和,=0.0106,（由内插得到）,第三节 不规则河道非均与流污染带的计算,80,剪切流速,u,计算：,继有：,第三节 不规则河道非均与流污染带的计算,81,据时间连续点源、考虑两岸反射的解析解式：,求断面,B,上,=0,和,=0.0106,处的浓度。,表 污染带的累积流量法计算,x,(m),h,h,o,(mg/L),4,Dx,n,E,1,E,2,C,(mg/L),500,0,(,左岸水边,),0.016,（,A,断面离岸,25m),2.93,0.00111,0,0.794,0.794,4.65,+1,0,0,-1,0,0,500,0.0106,(,y,=25m),2.93,0.00111,0,0.974,0.529,4.40,+1,0,0,-1,0,0,其中：,第三节 不规则河道非均与流污染带的计算,82,