非牛顿液体管道输送流变修正设计计算书

总则1.1编制目的本计算书针对非牛顿流体在长距离管道输送中因剪切变稀（或屈服）特性引起的流动阻力畸变问题，建立以幂律流变模型为核心、兼顾宾汉塑性与Herschel-Bulkley模型的流变修正计算体系。通过对输送介质流变参数的测定与回归、圆管层流流速分布与流量方程的推导、广义雷诺数与摩阻系数的修正，系统给出沿程摩阻、局部阻力、位差压头、泵扬程与功率，并进一步开展热力—流变耦合修正与多维度安全校核，使设计成果既满足现行标准要求，又真实反映浆体浆料的非牛顿流动行为。常规按牛顿流体假设的设计方法，仅以单一表观黏度参与计算，无法体现黏度随剪切速率的变化，在高浓度精矿浆体、煤浆、油水乳状液、聚合物溶液等工况下将产生显著偏差。本计算书的编制目的，即是消除这一偏差，为管道选型、泵组配置、运行控制与安全管理提供可靠依据。1.2适用范围本计算书适用于以水或其他液相为载体的高浓度固液两相似均质流浆体、以及油水乳化液、聚合物溶液等符合幂律、宾汉塑性或Herschel-Bulkley流变规律的非牛顿液体的长距离管道输送工程设计。设计内容包括输送工艺参数确定、流变修正水力计算、热力修正、强度与稳定性校核及运行控制建议。对于含粗颗粒、存在明显滑移或时变性（触变、震凝）的浆体，应按本计算书第3章的方法重新测定并确认模型适用性。1.3设计依据本计算书编制过程中执行及参考的主要标准、规范与文献如下：《输油管道工程设计规范》GB50253—2014（含2020年版局部修订）；《输气管道工程设计规范》GB50251—2015；《油气输送管道穿越工程设计规范》GB50423—2013；《浆体长距离管道输送工程设计标准》T/CECS98—2019；《石油天然气工业管线输送系统用钢管》GB/T9711—2017；《钢质管道外腐蚀控制规范》GB/T21447—2018；《油气输送用钢制感应加热弯管》SY/T5257—2016；《流体输送用不锈钢无缝钢管》GB/T14976及相关管件标准；《化工原理》《工程流变学》中非牛顿流体圆管层流与广义雷诺数的计算理论；介质流变特性室内试验报告及现场试运行数据。1.4设计原则（1）安全可靠：管道强度、刚度与稳定性按承载能力极限状态和正常使用极限状态双重校核，工作压力、水击压力与温度应力均留足安全裕度。（2）流变为本：以实测流变参数为基础选择本构模型，严禁以固定黏度替代非牛顿流体的真实流动特性。（3）经济合理：在满足不淤流速、不超压和不沉降的前提下优化管径与泵功率，降低全生命周期能耗。（4）可监测可调：设置流量、压力、温度与流变在线监测点，为运行工况的动态修正提供数据支撑。2主要符号与流变学基础2.1主要符号说明本计算书主要符号及其含义、单位如下表所示。未列出的符号在首次出现处予以说明。表2-1主要符号及计量单位符号含义单位τ剪切应力Paτ屈服应力Paτ管壁剪切应力Paμ动力黏度Pa·sμ塑性黏度Pa·sμ表观黏度Pa·sμ载体（水）动力黏度Pa·sK幂律稠度系数Pa·sⁿK参考温度下稠度系数Pa·sⁿn幂律流性指数—γ剪切速率s⁻¹γ管壁剪切速率s⁻¹ρ浆体密度kg/m³ρ固相真密度kg/m³ρ载体密度kg/m³C固相质量浓度%C固相体积浓度%Q体积流量m³/sV断面平均流速m/sD管道外径mmd管道内径mmL管道输送长度mδ管壁厚度mmΔp管道压降MPaΔz首末站高差mH泵所需扬程mf达西摩阻系数—Re雷诺数（含广义雷诺数）—g重力加速度m/s²σ管道环向应力MPa[σ]钢管许用应力MPaσ钢管最低屈服强度MPaϕ焊缝系数—w颗粒自由沉降末速m/sd颗粒中值粒径ma水击波速m/sT温度℃T环境（地温）温度℃T管道起点温度℃c浆体比热容J/(kg·K)G浆体质量流量kg/sβ稠度系数温度修正系数1/℃η泵机组效率—T阀门关闭时间sΔv流速变化量m/sζ局部阻力系数—K液体体积模量PaE钢材弹性模量Pan安全系数—2.2牛顿流体与非牛顿流体凡剪切应力与剪切速率之间呈线性关系、且斜率（黏度）不随剪切速率变化的流体称为牛顿流体，其内摩擦规律服从牛顿黏性定律： τ=μ d式中，τ为剪切应力（Pa）；μ为动力黏度（Pa·s）；du/dy为速度梯度即剪切速率（s⁻¹）。当剪切应力与剪切速率呈非线性关系，或黏度随剪切速率、剪切持续时间而变化时，该流体即为非牛顿流体。非牛顿流体的表观黏度不再是常数，而是剪切速率的函数，这正是流变修正计算的出发点。工程实践表明，高浓度精矿浆体、煤浆、钻井泥浆、油水乳状液、聚合物驱油剂溶液及食品浆料等，普遍呈现出剪切变稀、屈服或剪切增稠等非牛顿特性。若沿用牛顿流体假设，其水力计算与泵组选型将出现系统性偏差。2.3非牛顿流体分类按黏度函数与剪切持续时间是否相关，非牛顿流体可分为时变性流体与非时变性流体两大类。非时变性流体包括宾汉塑性流体、假塑性流体（剪切变稀）和胀塑性流体（剪切增稠）；时变性流体包括触变性、震凝性及黏弹性流体。本工程输送的精矿浆体在稳定输送温度下呈假塑性（剪切变稀）特征，且基本无时变效应，故按非时变性流体处理。从流变方程形式看，工程上常用的非时变性非牛顿流体本构模型包括幂律模型、宾汉塑性模型、Herschel-Bulkley模型和卡森模型，其选取方法见第3章。2.4流变修正的技术路线本计算书采用“参数测定—模型回归—流态判别—摩阻修正—热力耦合—安全校核”的技术路线：先由流变试验建立剪切应力—剪切速率关系曲线并回归本构参数；继而由圆管层流平衡方程推导幂律流体的流速分布与流量公式；再以广义雷诺数（Metzner-Reed）修正流态判别与摩阻系数；随后计入沿程温降对流变参数的影响完成耦合修正；最后从强度、刚度、沉降、水击与泵组工况五个维度开展安全核验。3流变本构模型与参数确定3.1常用流变本构模型幂律模型适用于中高剪切速率区间的假塑性与胀塑性流体，形式简单、参数少、外推性好，是浆体管道工程应用最广的模型： τ=K γ式中，K为稠度系数（Pa·sⁿ），反映流体稠度；n为流性指数（无因次），n<1表示剪切变稀，n>1表示剪切增稠，n=1退化为牛顿流体。与之对应的表观黏度为： μa宾汉塑性模型适用于存在明显屈服应力、且超过屈服值后呈线性流动的流体（如高浓度矿浆、水泥浆）： τ=τ式中，τ₀为屈服应力（Pa）；μₚ为塑性黏度（Pa·s）。当流体同时具备屈服应力和非线性流动特征时，应采用Herschel-Bulkley模型： τ=τ对于高剪切速率下表现出渐近黏度的流体，可采用卡森模型描述： τ=以上模型中，幂律模型（式2）、Herschel-Bulkley模型（式5）与宾汉塑性模型（式4）分别为本工程的预选模型。3.2模型比选与适用性判别模型比选以回归决定系数R²、参数显著性及物理合理性为依据：优先选择R²高、参数置信区间不含零、且外推至设计剪切速率区间仍稳定的模型。对同时可被多个模型拟合的数据，应根据本工程实际剪切速率范围（本工程管壁剪切速率约50～60s⁻¹）择取在该区间内拟合精度最高者。经对流变试验数据回归，宾汉塑性模型的屈服应力τ₀置信区间下限接近零，表明本浆体的屈服特征不显著；Herschel-Bulkley模型虽R²略高，但增加参数后对设计区间的预测改善有限；幂律模型在本工程剪切速率区间内R²达0.998以上，参数少、稳定性好。综合判断，本工程输送浆体为剪切变稀型假塑性流体，采用幂律模型作为流变修正计算的标准本构模型。3.3流变参数测定与回归流变参数采用旋转黏度计在恒温条件下测定，剪切速率取10～500s⁻¹区间，覆盖管道输送的实际剪切工况。试验数据经幂律模型双对数回归，得稠度系数K=1.85Pa·sⁿ、流性指数n=0.72，R²=0.9986。试验回归结果见表3-1，回归得到的本构曲线可作为后续水力计算的依据。表3-1浆体流变试验及幂律回归结果剪切速率γ̇/s⁻¹剪切应力τ/Pa表观黏度μₐ/Pa·s109.710.9712016.050.8035030.940.61910050.950.51020083.900.420500163.90.328由表3-1可见，随剪切速率增大，剪切应力非线性增长而表观黏度持续下降，呈现典型的剪切变稀特征。当剪切速率由10s⁻¹增至500s⁻¹时，表观黏度由0.971Pa·s降至0.328Pa·s，降幅达66.2%。这说明若以单一黏度值近似非牛顿流体，必然产生显著偏差。3.4流变参数的温度修正稠度系数对温度敏感，温度升高时流体稠度下降。工程上可用指数型温度修正关系描述： K(T)=K式中，β为稠度系数温度修正系数（本工程经试验取0.018/℃）；T为实际温度；T₀为参考温度。该关系将用于第6章的热力—流变耦合修正，以量化冬、夏工况下的摩阻差异。4设计基础条件4.1输送介质特性本工程输送介质为某铜精矿浆体，属似均质流。介质基本物性参数见下表。浆体密度由固相真密度与体积浓度按下式确定： ρm质量浓度与体积浓度的换算关系为： Cw经计算，当质量浓度Cw=60%、固相真密度取4200kg/m³、载体密度取1000kg/m³时，体积浓度Cv=26.3%，浆体密度ρₘ=1840kg/m³，与实测值吻合。介质物性汇总见表4-1。表4-1输送介质物性参数项目数值输送介质铜精矿浆体（似均质流）固相真密度ρₛ/(kg/m³)4200载体（水）密度ρ_w/(kg/m³)1000浆体密度ρₘ/(kg/m³)1840固相质量浓度Cw/%60固相体积浓度Cv/%26.3中值粒径d₅₀/μm45幂律稠度系数K/(Pa·sⁿ)1.85幂律流性指数n0.72浆体比热容c/(J/(kg·K))35004.2管道与设备条件管道线路全长8.0km，采用DN250无缝钢管，规格φ273.1×8.0mm，内径d=257.1mm；管材为L360（X52），最低屈服强度σs=359MPa，焊缝系数取1.0。管道全线为埋地敷设，土壤温度夏季按15℃、冬季按−5℃考虑。首站设置2台同型号离心泵机组（1用1备），配套电机按泵轴功率并留取裕度选型。管道设计内压取6.3MPa，设计温度范围−20～60℃。4.3设计参数汇总表4-2主要设计参数汇总参数取值单位管道外径D273.1mm管道内径d257.1mm壁厚δ8.0mm管材L360（X52）—最低屈服强度σs359MPa设计内压P6.3MPa输送长度L8000m设计流量Q0.085（306m³/h）m³/s首末站高差Δz30m总传热系数K_h1.2W/(m²·K)环境温度（夏/冬）15/−5℃浆体质量流量G156.4kg/s5流变修正水力计算5.1管径与平均流速管道流通面积与断面平均流速为：A=πd²/4=π×0.2571²/4=0.05192m²；V=Q/A=0.085/0.05192=1.637m/s。计算得平均流速V=1.637m/s，处于精矿浆体长距离管道输送的合理流速区间（1.2～2.5m/s），同时高于后文校核得到的临界不淤流速。5.2流态判别非牛顿流体的流态判别不能直接沿用牛顿流体雷诺数，需采用广义雷诺数（Metzner-Reed）： Re相应的幂律流体层流—湍流临界雷诺数采用Ryan-Johnson关系式： Re将ρ=1840kg/m³、V=1.637m/s、D=0.2571m、K=1.85Pa·sⁿ、n=0.72代入式16，得广义雷诺数Re_MR=1177。由式17计算得临界雷诺数Re_c=2269。由于Re_MR=1177<Re_c=2269，判定管道内浆体处于层流状态。层流工况有利于减小摩阻和磨蚀，是本工程期望的流态。5.3壁面剪切速率与剪切应力幂律流体在圆管中流动时，管壁处的剪切速率按下式计算： γw代入数据得γ̇_w=[(3×0.72+1)/(4×0.72)]×(8×1.637)/(0.2571)=1.0972×50.94=55.9s⁻¹。管壁剪切应力由圆管层流平衡方程确定，管壁剪应力与压降的关系为： τw由式2反向计算，管壁剪应力τw=K·γ̇wⁿ=1.85×55.9^0.72=33.5Pa。该值等于后文由压降反算的管壁剪应力，二者相互印证，表明计算自洽。表观黏度μₐ=τw/γ̇w=0.600Pa·s。5.4沿程摩阻损失（流变修正）幂律流体在圆管层流状态下的流量方程由层流流速分布积分得到，其流速分布为： u=n对流速分布积分得到体积流量与压降的关系： Q=n进一步可得断面平均流速表达式： V=Q层流状态下达西摩阻系数与广义雷诺数成反比： f=64沿程水头损失采用计入流变修正的达西—魏斯巴赫公式： hf将幂律流体参数代入并将流量方程化为压降显式形式，可直接计算沿程压降： Δp=4KL计算得：8V/D=50.94s⁻¹；f=64/1177=0.0544；h_f=0.0544×(8000/0.2571)×[1.637²/(2×9.81)]=0.0544×31116×0.1366=231.1m。沿程压降Δp=4KL/D×[(3n+1)/4n]ⁿ×(8V/D)ⁿ=4.17MPa。经式11反算，管壁剪应力为33.5Pa，与5.3节结果一致。5.5局部阻力损失局部阻力损失由弯管、阀门、进出站管件等引起，按局部阻力系数法计算： hj长距离浆体管道局部损失通常较小，为便于安全取用，本工程按沿程水头损失的5%计取，即h_j=5%×231.1=11.6m。5.6位差与总扬程首末站高差Δz=30m（末站高于首站），位差所需压头为： Δp泵站所需总扬程为沿程、局部与位差压头之和： H=h计算得H=231.1+11.6+30=272.6m。泵出口压力为： P2代入得P₂=1840×9.81×272.6=4.92MPa，低于设计内压6.3MPa，满足要求。5.7泵功率泵机组轴功率按下式计算： P=ρ取泵机组效率η=0.75，代入得P=1840×9.81×0.085×272.6/(1000×0.75)=558kW。考虑1.1～1.15的电动机容量裕度，选用配套电机额定功率630kW；首站设2台同型号机组（1用1备），单机轴功率558kW。6热力—流变耦合修正计算6.1沿程温降计算浆体在输送过程中的散热采用舒霍夫温降公式计算： TL式中，T_L为终点温度；T_a为环境温度；T₀为起点温度；K_h为总传热系数；G为质量流量；c为比热容。夏季工况取T_a=15℃、T₀=25℃，代入得T_L=15+(25−15)×exp(−0.0150)=24.85℃；冬季工况取T_a=−5℃、T₀=15℃，代入得T_L=−5+(15+5)×exp(−0.0150)=14.70℃。由于管输量较大、沿程温降很小，全年温度基本稳定在设计区间内。6.2温度对流变参数的影响温度降低使浆体稠度升高，需按式26对幂律稠度系数进行修正。以夏季平均温度24.93℃为参考，冬季平均温度为14.85℃，温差10.08℃，则冬季稠度系数K_冬=1.85×exp(0.018×10.08)=2.22Pa·sⁿ，较夏季升高20.0%。流性指数n在5～25℃区间内的变化不显著，按常数0.72取值。6.3耦合修正结果与冬季工况以冬季修正后的稠度系数重新计算水力工况：广义雷诺数Re_MR=1177×(1.85/2.22)=982，摩阻系数f=64/982=0.0652，沿程水头损失h_f=277.1m，较夏季工况（231.1m）增加19.9%，相应沿程压降由4.17MPa增至5.00MPa。据此，运行中应按冬季工况校核泵组能力，并预留不小于20%的扬程裕度，避免因温度降低导致的输送能力不足与超压风险。7多维度校核计算7.1管道强度校核管道在设计与水击工况下的环向应力按下式计算： σθ钢管许用应力按GB50253—2014第5.2.1条确定，管线一般地段设计系数K_d取0.72，焊缝系数φ取1.0，则[σ]=0.72×1.0×359=258.5MPa。设计内压下环向应力σθ=6.3×257.1/(2×8.0)=101.2MPa，安全系数为： ns计算得安全系数n_s=258.5/101.2=2.55，管壁应力水平仅为许用应力的39.2%，强度裕度充足。考虑水击压力（详见7.4节）叠加后，最大环向应力为91.9MPa，仍在许用应力范围内。7.2壁厚与刚度校核按GB50253—2014第5.4.1条，直管段计算壁厚为： δ=P D代入设计内压P=6.3MPa、外径D=273.1mm、许用应力[σ]=258.5MPa，得计算壁厚δ=3.33mm。计入浆体磨蚀裕量1.5mm，所需最小壁厚为4.83mm。考虑长距离矿浆输送的磨蚀、第三方破坏防护及埋地施工刚度要求，选用壁厚8.0mm。管道径厚比D/δ=34.1，远小于100的刚度限值，抗椭圆化与抗外压屈曲能力良好。7.3临界流速与颗粒沉降校核浆体能否维持似均质流、不产生沉积，取决于颗粒自由沉降末速与管流紊动强度的对比。颗粒自由沉降末速按斯托克斯公式计算： w0代入d₅₀=45μm、ρₛ=4200kg/m³、ρ_w=1000kg/m³、μ_w=1.0×10⁻³Pa·s，得w₀=(45×10⁻⁶)²×3200×9.81/(18×1.0×10⁻³)=3.53×10⁻³m/s，即3.53mm/s。沉降末速极小，说明颗粒悬浮性好、浆体为似均质流，其临界不淤流速可按经验下限取V_c=1.0m/s。设计流速1.637m/s大于临界流速，裕度63.7%，可有效避免管道底部沉积。7.4水击压力校核突然停泵或阀门快速关闭将产生水击。先计算浆体中的水击波速： a=K取液体体积模量K_w=2.2×10⁹Pa、浆体密度ρₘ=1840kg/m³、管径D=273.1mm、钢弹性模量E=2.06×10¹¹Pa、壁厚δ=8.0mm，代入得波速a=936m/s。水击相长t_c=2L/a=2×8000/936=17.1s。将阀门关闭时间控制为T_s=60s（大于相长），属间接水击，水击压力为： Δp计算得水击压力Δpₐ=1840×1.637×(2×8000/60)=0.803MPa。则工作压力与水击压力叠加后的最大压力为4.92+0.80=5.72MPa，低于设计内压6.3MPa，水击工况安全。工程上应设置缓闭阀门、空气阀与水击泄放装置，确保关阀时间不小于60s。7.5泵机组工况校核夏季设计工况泵扬程272.6m、轴功率558kW；冬季工况扬程升至321m、轴功率约658kW。为此，配套电机额定功率应按冬季工况选型并留取裕度，或将泵组设计为变频调速，通过降低转速适应冬季摩阻上升，避免电机过载与管道超压。泵组工作点应落在高效区（效率不低于75%），并避开喘振与汽蚀区间。首站采用1用1备配置，可保证在单泵检修时维持不低于90%的输送能力。8安全核验8.1安全系数与设计裕度本工程从压力、强度、流速、水击与温降五个方面设置安全裕度，各项核验结果汇总见表8-1。各项指标均满足现行标准及工程实践要求，设计裕度合理。表8-1多维度安全核验结果汇总核验项目计算值控制值结论裕度/安全系数压力校核工作压力4.92MPa设计内压6.3MPa通过1.28强度校核环向应力101.2MPa许用应力258.5MPa通过2.55刚度校核径厚比D/δ=34.1限值100通过2.93沉降校核设计流速1.637m/s临界流速1.0m/s通过1.64水击校核最大压力5.72MPa设计内压6.3MPa通过1.10温降校核冬季摩阻277.1m夏季摩阻231.1m通过≥1.208.2极限工况组合分析考虑最不利的极限工况组合，即“冬季最低地温+最高浓度+最大流量+阀门误动作”。该组合下浆体稠度最大、摩阻最高，瞬时压力最高。经核算，冬季最高摩阻对应压力为5.00MPa，叠加最大水击压力0.80MPa后为5.80MPa，仍低于设计内压6.3MPa。针对该工况，运行规程已设置浓度上限报警、流量联锁与缓闭阀保护，可将风险控制在可接受水平。8.3风险识别与防护措施（1）沉积堵管风险：低流量或长时间停输时颗粒可能沉降，应在停输前用清水顶替，并设置低流速报警；再启动时按低速启动、逐步升速。（2）超压风险：设置压力高报警、泵组联锁停泵与泄放阀，水击防护采用缓闭阀与空气阀组合。（3）磨蚀风险：弯头、三通等易磨部位加厚管壁或采用耐磨衬里，定期测厚。（4）温度风险：对冬季工况进行专项核算，必要时对埋地管道加强保温或增设加热措施。（5）流变漂移风险：定期取样复测流变参数，当K、n偏离设计值超过10%时重新核算水力工况。9结果分析与敏感性分析9.1计算结果汇总表9-1流变修正水力计算结果汇总计算项目结果断面平均流速V/(m/s)1.637广义雷诺数Re_MR1177流态判别层流管壁剪切速率γ̇_w/(s⁻¹)55.9管壁剪应力τ_w/Pa33.5表观黏度μₐ/(Pa·s)0.600沿程水头损失h_f/m231.1局部水头损失h_j/m11.6位差压头Δz/m30.0泵所需扬程H/m272.6沿程压降Δp/MPa4.17泵出口压力P₂/MPa4.92泵轴功率P/kW5589.2敏感性分析为量化流变参数与输送工况对水力计算结果的敏感性，以设计工况为基准，分别对稠度系数K、流性指数n、流量Q及温度进行±10%扰动，结果见表9-2。表9-2流变与工况参数敏感性分析参数基准值变化沿程压降/MPa相对偏差稠度系数K1.85+10%4.59+10.0%稠度系数K1.85−10%3.75−10.0%流性指数n0.72+10%5.45+30.6%流性指数n0.72−10%3.19−23.4%流量Q0.085+10%4.47+7.2%温度（冬季）24.9℃−10℃5.00+19.9%由表9-2可见：流量与稠度系数的影响近似线性，而流性指数n的敏感性最高。当n增大10%时，沿程压降抬升30.6%；减小10%时，压降下降23.4%。这说明流性指数是决定非牛顿流体水力特性的关键参数，其测定精度直接决定设计成果的可靠性。工程上必须保证流变参数在规定剪切速率区间内的高精度测定，任何以经验值替代实测值的做法都可能带来重大偏差。9.3与牛顿流体假设的对比误差分析为评估传统牛顿流体假设的误差，分别以两种方式近似：一是取管壁表观黏度μₐ=0.600Pa·s，二是取低剪切速率（10s⁻¹）下的表观黏度0.971Pa·s。计算对比如下：（1）按管壁表观黏度近似：Re=1291，f=0.0496，h_f=210.8m，Δp=3.80MPa，较流变修正值4.17MPa偏低8.9%。（2）按低剪切速率表观黏度近似：Re=798，f=0.0802，h_f=341.1m，Δp=6.16MPa，较流变修正值偏高47.6%。上述结果表明，用单一黏度值替代幂律流变特性，会因剪切速率取值不同而产生−8.9%～+47.6%的偏差，且偏差方向与取值位置相关，具有很强的不确定性。这正是本计算书坚持采用流变修正计算的工程价值所在。10工程指导意义与实施建议10.1设计建议（1）管径优化应以“流速不低于临界不淤流速、摩阻与能耗最小”为目标，通过不同管径方案的技术经济比选确定；本工程DN250方案兼顾了输送能力、流速裕度与磨蚀控制。（2）泵组应按冬季最不利工况选型或采用变频调节，预留不小于20%的扬程裕度，避免因温度降低导致输送能力不足。（3）管材选择应兼顾强度、刚度与耐磨蚀性能，弯头、三通等局部部位应加厚或采用耐磨衬里。（4）水击防护应采用缓闭阀、空气阀与泄放装置组合措施，确保关阀时间不小于60s。10.2运行控制建议（1）启动采用低速启动、逐步升速方式，避免瞬时剪切与压力冲击；停输前应用清水顶替管道内浆体，防止沉积堵管。（2）建立以流量、压力、温度为输入