井筒摩阻计算原理与方法

1、井筒摩阻计算第一节水头损失及其分类实际流体具有粘性，在通道内流动时，流体内部流层之间存在相对运动和流动阻力。流动阻力和水头损失的规律，因流体的流动状态和流动的边界条件而异一、水头损失分类流体在流动的过程中，在流动的方向、壁面的粗糙程度、过流断面的形状和尺寸均不变的均匀流段上产生的流动阻力称之为沿程阻力，或称为摩擦阻力。沿程阻力的影响造成流体流动过程中能量的损失或水头损失（习惯上用单位重量流体的损失表示）。沿程阻力均匀地分布在整个均匀流段上，与管段的长度成正比，一般用hf表示。另一类阻力是发生在流动边界有急变的流场中，能量的损失主要集中在该流场及附近流场，这种集中发生的能量损失或阻力称为局部阻力

2、或局部损失，由局部阻力造成的水头损失称为局部水头损失。二、水头损失分类1.沿程阻力损失对于圆管:hfhf入1v24R2glv2d2g式中：1管长；R水力半径；d管径；v断面平均流速；g一重力加速度；沿程阻力系数，也称达西系数。一般由实验确定。式中的无量纲系数不是一个常数，它与流体的性质、管道的粗糙程度以及流速和流态有关，在大多数工程问题中，hf确实与v2成正比。此外，这样做可以把阻力损失和流速水头合并在一起，便于计算。2.局部阻力损失hj2g式中：一一局部阻力系数，一般由实验确定。整个管道的阻力损失，应该等于各管段的沿程损失和所有局部损失的总和。第二节粘性流体流动流态一、粘性流体流动流态当流速

3、较小时，沿程损失与流速一次方成正比，当流速较大时，沿程损失几乎与流速的平方成正比，如图所示，并且在这两个区域之间有一个不稳定区域。当阀门B慢慢打开，并打开颜色水阀门D,此时管中的水流流速较小，可以看到玻璃管中一条线状的颜色水。它与水流不相混合，如图63(b)所示。从这一现象可以看出，在管中流速较小时，管中水流沿管轴方向呈层状流动，各层质点互不掺混，这种流动状态称为层流。当阀门B逐渐开大，管中的水流流速也相应增大。此时会发现，在流速增加到某一数值时，颜色水原直线的运动轨迹开始波动，线条逐渐变粗，如图63(c)所示。继续增加流速，则颜色水迅速与周围的清水混合，6-3(d)所示。这表明液体质点的运动

4、轨迹不规则，各层液体相互剧烈混合，产生随机的脉动，这种流动称为紊流。水流流速从小变大。沿程阻力曲线的走线为A-B-C一D。如图62所示。若实验时流速由大变小。则上述观察到的流动现象以相反的程序重演，但有紊流转变为,层流的流速Vc(下临界流速)要小于由层流转变为紊流的流速Vc(上临界流速)。如图6-2所示。沿径阻力曲线的走线为D-C-A。如图6-2所示。实验进一步表明，同一实验装置的临界流速是不固定的，随着流动的起始条件和实验条件不同，外界干扰程度不同，其上临界流速差异很大，但是，其下临流流速却基本不变。在实际工程中，扰动是普遍存在的，上临界流速没有实际意义，一般指的临界流速即指下临界流速。二、

5、流态的判别准则流态不仅与断面平均流速V有关系，而且与管径d、液体粘性、密度有关。即流态既反映管道中流体的特性，同时又反映管道的特性。将上述四个参数合成一无量纲数(无具体单位)，称为雷诺数，用Re表示。vdpvdp式中：Q雷t若数,V流速，m/s,p流体密度，kg/m3,流体粘度，Pa.s;运动粘度，m2/s注，粘度单位1Pa.s=1N*s/m2=10P=103cp=1Kcps=1kg/(m*s)1N=1kg*m/s2对应于临界流速的雷诺数，称为临界雷诺数，通常用表示。大量实验表明，在不同的管道、不同的液体以及不同的外界条件下临界雷诺数不同。通常情况下，临界雷诺数总在2300附近，Rec=230

6、0【例61】有一直径d25mm的室内上水管，如管中流速v1.0m/s水温t10c。(1) .试判别管中水的流态；(2) .试求管内保持层流状态的最大流速为多少？解：(1)10c时，水的运动粘性系数v1.31106m2/s,此时，管内雷诺数rvd100.025191002300,故管中水流为紊流。e1.31106(2)保持层流的最大流速就是临界流速，Re空2300所以vc23001.311060.0250.12ms第三节沿程水头损失与切应力的关系一、均匀流动方程式沿程阻力（均匀流内部流层间的切应力）是造成沿程水头损失的直接原因。建立沿程水头损失与切应力的关系式，再找出切应力的变化规律，就能解决沿

7、程水头损失的计算问题。在圆管恒定流均匀流段上设11和22断面，如图所示。作用于流段上的外力：压力、壁面切应力重力相平衡。即：p1Ap2AAlcosw10QLm-ra一工1JQ式中w壁面切应力，湿周，a圆管倾斜角，A圆管断面截面积，1管段长度。由几何关系得:1cosz1z2,除以A整理得：R7卫2w1zlz2AA并由断面1和断面2的能量方程得:7R7压z1Z2hf,故:hfw1RhfR1式中：R水力半径，RJ水力坡度，J、圆管过流段面上切应力分布0即圆管均匀过流断面上切应力呈直线分布，管轴处0,管壁处切应力达最大值三、壁剪切速度下面在均匀流动方程式的基础上，推导沿程摩阻系数和壁面切应力的关系。.

8、lV2以务工、工出口d2g将J代入均匀流动方程式，整理得：Jw_具有速度的量纲，称为壁剪切速度擦速度）。则:上式是沿程摩阻系数和壁面切应力的关系式，该式在紊流的研究中广为引用。第四节、圆管中的层流运动层流常见于很细的管道流动，或者低速、高粘流体的管道流动，如阻尼管、润滑油管、原油输油管道内的流动。圆管中层流运动的流动特征如前述，层流各流层质点互不掺混，对于圆管来说，各层质点沿平行管轴线方向运动。与管壁接触的一层速度为零，管轴线上速度最大，整个管流如同无数薄壁圆筒一个套着一个滑动。各流层间切应力服从牛顿内摩擦定律，即满足式dudrdu由yr。r圆管层流的断面流动分布因讨论圆管层流运动，所以可用牛

9、顿内摩擦定律来表达液层间的切应力:dududydr式中为动力粘性，u为离管轴距离r处的切应力（即离管壁距离y处）的流速。对于均匀管流而言，在半径等于r处的切应力应为：于均匀管流而言，根据式（621）,在半径等于r处的切应力应为:durdr2积分得:Jr2C4利用管壁上的边界条件，确定上式中的积分常数当r0时u0,得：CJ.2J/22、uL0r)66所示。过流断面上上式表明，圆管中均匀层流的流速分布是一个旋转抛物面，如图流速呈抛物面分布，这是圆管层流的重要特征之一。将r0代入上式，得到管轴处最大流速为UmaxJr20平均流速为:AudAA：r2rdr0r210Lro2rJ(r。2r2)2rdr三

10、、圆管层流的沿程阻力损失)2将直径d代替式（634）中的2r032进而可得水力坡度32以Jhf/l代入上式,可得沿程阻力损失为:这就从理论上证明了圆管的均匀层流中.沿程阻力损失hf与平均流速v的一次方成正比，这与雷诺实验的结果相符。上式还可以进一步改写成达西公式的形式hf32d_264lvvdd2g26419Red2gitd该式为达西和魏斯巴哈提出的著名公式，雷诺数的函数，与管壁粗糙情况无关。例题62设有一恒定有压均匀管流.已知管径v0.12m/s,水温t10c时水的运动粘度由上式可得:64Re此公式表明圆管层流中的沿程阻力系数只是d20mm,管长l20mm,管中水流流速，62v1.30610

11、m/s。求沿程阻力损失。解：Re64Revd0.120.021.30610618382300为层流6418380.0352一一2hf0.026mH2Olv20(0.12)0.035d2g0.0229.8第五节紊流运动分析一、粘性底层在紊流运动中，并不是整个流场都是紊流。由于流体具有粘滞性，紧贴管壁或槽壁的流体质点将贴附在固体边界上，无相对滑移，流速为零，继而它们又影响到邻近的流体速度也随之变小，从而在紧靠近面体边界的流层里有显著的流速梯度，粘滞切应力很大，但紊动则趋于零。各层质点不产生混掺，也就是说，在取近面体边界表面有厚度极薄的层流层存在，称它为粘性底层或层流底层，如图68所示。在层流底层之

12、外，还有一层很簿的过渡层。在此之外才是紊层，称为紊流核心区。在粘性底层中，速度按线性分布，在壁面上速度为零。粘性底层虽然很薄，但它对紊流的流速分布和流速阻力却有重大的影响。紊流紊流过流断面上流速成对数曲线分布，同层流过流断面上流速成抛物线分布相比，的流速分布均匀很多。第六节沿程水头损失系数的变化规律是计算沿程损失的关键。但由于紊流的复杂性，直到目前还不能像层流那样严格地从理论上推导出适合紊流的值来，所以值的确定，现有的方法仍然只有经验和半经验方法。一、阻力系数的影响因素在紊流中，除与反映流动状态的雷诺数有关之外，还因为突入紊流核心的粗突起会直接影响流动的紊动程度，因而壁面粗糙度是影响阻力系数的

13、另一个重要因素。可以采用一个指标即检验突起高度ks来表示壁面的粗糙程度，ks称为绝对粗糙度。绝对粗糙度具有长度量纲，所以仍感到有所不便，因而引入了量纲的相对粗糙度，即ks与直径(或半径)之比ks/d(或ks/r0),它是一个能够在不同直径的管流中用来反映管壁粗糙度的量，由以上分析可知，影响紊流沿程阻力系数的因素是雷诺数和相对粗糙度ks/d,写成函数关系式为f(Re,ks.d)二、尼古拉兹实验尼古拉兹在1933在人工粗糙管中系统地进行了沿程阻力系数和断面流速v的测定。他的实验涉及的参数范围比较大，相对粗糙度范围为ks/d130-11014；雷诺数范围为Re500106纵坐标为lg(100),横坐

14、标为lgRe,再算出和Re,取对数点绘在坐标纸上，就得到f(Re,ks/d)曲线，即尼古拉兹曲线图。管道的流动可分为五个区域。第一个区域是层流区，对应的雷诺数Re2300(lgRe3.36)，试验点均落在直线ab上。表明与相对粗糙ks/d无关，只是Re的函数，并符合64/Reo还可知，沿程阻力损失hf与断面平均流速v成正比，这也与雷诺试验的结果一致。第二个区域为层流与紊流之间的过渡区，Re23004000(lgRe3.363.6)试验点落在bc附近，表明与相对粗糙ks/d无关，只是Re的函数。此区是层流向紊流过渡，这个区的范围很窄，实用意义不大，不予讨论。第三个区域为紊流光滑区，Re4000(

15、lgRe3.6)不同的相对粗糙管的试验点都先后落在同一条cd线上。表明与相对粗糙ks/d无关，是Re的函数。随着Re的增大，k$/d大的管道，实验点在Re较低时便离开此线，而ks/d较小的管道，在Re较大时才离开。第四个区域是紊流过渡区，不同的相对粗糙管实验点分别落在不同的曲线上。表明既与Re有关,又与ks/d有关。第五个区域是紊流粗糙区，不同的相对粗糙管实验点分别落在不同水平直线上，表明与ks/d有关，与Re无关。在这个阻力区里，对于一定的管道(ks/d一定)，是常数。沿程水头损失与流速的平方成正比，故有称为阻力平方区。三、速度分布1 .紊流光滑区5.75lg”5.5v2 .紊流粗糙区uy半

16、经验公式：5.75lg8.48Vksd)0u指数公式：umax四、的半经验公式1 .光滑区沿程摩阻系数尼古拉兹光滑管公式2lgRe12.512 .粗糙区沿程摩阻系数尼古拉兹粗糙管公式2lg3.7dks五、工业管道实验曲线忌浙巴a.OLC：A-fe-tfjnnMH一同，；歌mHR娶里必塞=1：J!：=aEE；r_Tikni；就加&例5我M21C.N弧015Q.8C.dlCl3fled0g1优2门2门产)C.CD(J31WB9M|Hnii|2(107)d001fl.。，g丸OOC4Q.OJOCOlt.CMC!吟1在尼古拉兹试验中，紊流有明显的光滑区。因为人工粗糙砂粒的直径是一致的。只要粘性底层的厚

17、度大于砂粒直径，流动就处于光滑区。而工业管道、出于工业加工的缘故，不可能制造出粗糙度完全一致的管道。壁面的粗糙部分，从微观上讲，高低不一。因此没有明显的光滑区，或者光滑区的跨越范围很窄，无法进行对比。人工粗糙区，无沦是人工管道，还是工业管道，由于粗糙面完全暴露在紊流中，其水头损失的变化规律也是一致的。因此，在相同的情况下。可用人工管道的相对粗糙度来表示工业管道的相对粗糙度，即当量粗糙度。当量粗糙度是用直径相同，在紊流粗糙区相同的人工管道的粗糙度ks,来定义该工业1939管道的粗糙度，表列出了常用工业管道的当量粗糙度。由表中数据可知，工业管道的计算方法与人工管道的计算方法一样。但尼古拉兹阻力系数

18、公式在紊流过渡区是不适用的。年，柯列勃洛克和怀特给出了工业管道紊流区中的计算公式:21lgks3.7d2.51Reks工业管道的当量粗糙度。常用工业管道的当量粗糙管道材料ksmm管道材料ksmm管道材料ksmm新氯乙烯00.002钢管0.046新铸铁管0.150.5铅铜玻璃0.01涂沥青铸铁管0.12旧出铁管11.5镀锌钢管0.15混凝土管0.33.0为了将式(661)图形化，1944年，美国工程师穆迪以该公式为基础，以当量粗糙度为参数，用对数坐标绘制出工业管道摩阻损失系数曲线图，即穆迪图，见图611。六、沿程摩阻系数的经验公式布拉修斯(Blasius)公式1931年德国水力学家布拉修斯在总结

19、前人实验的基础上总结并提出了紊流光滑区经验公式(662)0.3164C0.25Re该式形式简单，11算方便。在Re105范围内，有极高的精度，得到广泛的应用。希弗林松公式_.,ks、0.250.11()(663)d希弗林松粗糙区公式，该式形式简单，计算方便，工程界经常采用。谢才公式和谢才系数22ghfhf.将达西公式(62)变换形式vd，以d4R,J,代入上式，整理得:llv8g.RJC.RJ(664)式中：v断面平均流速；R水力半径;J水力坡度；c谢才系数。上式最初是1769年法国工程师谢才直接根据渠道和塞纳河的实测资料提出的，是水力学最古老的公式之一，称为谢才公式。(665)式（665）给

20、出了谢才系数C和沿程摩阻系数的关系，谢才系数含有阻力的因素。流动阻力越大，谢才系数越小，反之亦然。1895年，爱尔兰工程师曼宁提出了计算谢才系数的经验公式：(666)C1R16n式中：n反映壁面粗糙性质并与流动性质无关的系数，称为粗糙系数。Weymouth方程：0.0094073d为Darcy-Weisban摩阻系数，无量纲；D为管道内径，mWeymouth方程取管壁粗糙度k=0.0508mm（美国取k=0.02mm,苏联取k=0.03mm）,该方程适用于管径小、输量不大、净化程度较差的矿区集气管网。适用于管径范围254-508mm,阻力平方区Zegarola方程:114991.884*log

21、10-Re*.该方程适用于3.2*10A4Re3.5*10A7Colebrook-White方程:在计算摩阻系数的方程中，CW方程最常用，是一个经验方程1k/D2log10(3T2.51)Re*适用于各个雷诺数，相对粗糙度0.00001-0.03之间，水力光滑区，混合摩擦区，阻力平方区。Gerg方程:k/D)n3.711.4992，一10g10Kn为Darcy-Weisban摩阻系数，无量纲；D为管道内径，m;f为阻力因子（等同于输送效率），无量纲；n为哥指数，描述从光滑管到粗糙管转变的变化剧烈程度；Re为雷诺数，无量纲。方程的特点：（1）对于完全粗糙流态（例如非常大的雷诺数），与Colebr

22、ook-White公式是一致的；（2）f也成为气流因子，描述完全管流中产生的额外损失；（3）对于水力光滑区（粗糙度k=0）f=1,公式退化为Zagarola方程，低雷诺数粗糙管流同样可以。上=-2b-$,53lU(为+r)玷=/地+一U4显式:_hlO一外_1.499实例:计算数据见表L取GERG方程中参数”二1,广1进行求解s1道的际据当量粗程度上Re管径1)M是否有内涂层防护50.8212k1070.9664LJ由号art计箕结果用对调差Hi（时商1式仲11.079勃？8OT41SS1E-02减值如情0.1umnnue内语为制U7938?B418S1E-020.512出代反法：理re法程计

23、算ZA12JM27W67014221-0234西53斗。3GLOWM73式E1X179387W31891E-023492ME-Q?。叩1714式住1.079337和941国就一。才91M的E750.MV765S武俘1J079789941881-02OjMIkg其他公式:计算公式适用范闱Weymtxith管径范用为25450gmm，阻力平方区PantA）相对粗糙度等f0.0000L雷诺数工乂10*14*IT（水力光滑区）由*阴相衬粗精度V0.00001.宙讲数5*1门阻力平方区）相对粗髓度范圉为0（X）（）05-0.0004之间.水力AGA光滑区相对粗糙曳葩国。0000503之间，阻力平方区Cu

24、kbrook相对粗粗度范围0.0000】-0.03之间，水力比滑区，混合摩擦区,阻力平方区七、当量粗糙度当量粗糙度是将工业管道绝对粗糙度（e或者ks表示）折算成人工粗糙管的绝对粗糙度。具体做法：是在流态处于完全粗糙管范围内，选用与实际管道摩阻系数（用符号f表示）值相等和管径相同的人工粗糙管的绝对粗糙度当做实际管道的当量粗糙度。表1管道当粗嘴睫:后值表管通种费加工及帧用献猊kJmibL1专化战因平均值段瑞.相.幅肯媪的.常卅的莘沐和制的(1OfiLClQl讯口修的管新的.光甫的.祭饰收制的aoois-(iO6ao?无圾诩注L髀的.清勖的.被段良好的1用过几隼后加汉清人的.信谢比的.橙他钙馋的.弧

25、垢ivtrj(LO2-LIJ60.15-0.3ausa2库接蟒管和抑抻记菅1,卜口处坤桂制管|R咨汉向焊装的副U）新的.清沾的翁经清洗后桥惊木星著的旧曾I，轻度精性的in皆土中管锈曲的旧营:卜LJ*寻根呼m辄缱和植健都是焊接的门做缝葬帐椅继脚搭-rt-iJ口;班隆倬梯横建掰掂两排或者百排以上钢灯n03-c.ID1-4L20.2-fll7D.AL5a”lnCi，L2-28105Q150,51OL7I.2I.8段株捆管1钝洋面地冷消萍的新甘2 .彼抨面一股的断管3 .用过几年后的阳营a07-0.i0.Lfl,20,4-fL7aifa5情帙哲1.新管2.涂沙育的舟管L/ifltf的阻管D.Z-ft5a.i*ais0.12-0.503a阳混凝L转及捆的版混t.ff1,无捕藁面层”卜胸模他胞闱n口好*崔健,淞I*木梗检施ra-ft2 .向小版层井拄技尤3 .白M犁面层11）表而用而刷刷过拌西杼加枕北表加阳圳刷刷通但去不拉*:仇工、根。L0L0,2什IS07-1S费工口0.7L20.7L2揩报必轴桂4pg下面平一个魏点为何作具体示能“例,计舞:从井U开鲂第一井段MM)m距离捌点的油管摩机系教,ftI.A,20l(l-+a/Ak7r0.7k.1L”奉敷海符了.tdl_Hf(，市jr|w-a4MP-itT尸2iOIsgp+15.KMPm尸.=1,54i兀.1=LMiZ=(k7C1B.(MJ