天然气设计计算公式大全

天然气设计计算公式大全说明:天然气设计计算涉及的主要内容可划分为以下几大领域：管道水力计算（威莫斯、潘汉德尔A、潘汉德尔B、AGA、IGT等公式）是输气管网和城市燃气管网设计的核心技术。这些公式分别适用于不同的管径、压力和流态范围。采用科尔布鲁克-怀特方程及其显式近似（如海兰德公式）可精确计算管道摩阻系数。管道局部阻力通过当量长度法或局部阻力系数法折算。物性计算（压缩因子、密度、粘度、热值、沃泊指数等）贯穿于所有工程设计环节。AGA8号报告提供的高精度计算方法适用于工艺设计，而CNGA简化公式适用于工程快速估算。设备设计计算涵盖了压缩机（压头、功率计算）、分离器（液滴沉降速度）、脱水塔（塔径、填料层高度）、调压阀（Kv值、阻塞流）等关键工艺设备。传热计算和温度分布计算对埋地管道的设计至关重要。苏霍夫温度分布公式和焦耳-汤姆逊效应温降公式在埋地管道温度场分析中应用广泛。液化工艺和气质分析则涉及更复杂的相平衡计算和热物性计算，是LNG工厂设计和高精度天然气贸易计量（孔板流量计、音速喷嘴）的基础。

目录TOC\o"1-1"\h\z\u第一部分：管道水力计算公式 第一部分：管道水力计算公式1.1威莫斯（Weymouth）公式——输气管道流量计算威莫斯公式是小直径高压输气管道（管径约150mm~600mm，压力约7MPa）常用的设计公式。Q=433.5×符号意义：Q—气体流量，标准状态下m³/d433.5—经验常数Tb—基准温度，K（通常为288.7K或293K）Pb—基准压力，kPa（通常为101.325kPa）P1—管道起点绝对压力，kPaP2—管道终点绝对压力，kPaG—气体相对密度（空气=1，无量纲）Tf—管内气体平均流动温度，KL—管道长度，kmZ—平均压力下的气体压缩因子，无量纲D—管道内径，mm应用场景：小直径（≤600mm）、中等压力（≤7MPa）的输气管道设计计算。限制条件：未考虑高程变化，不适合大口径高压管道；经验常数仅适用于工程单位制。工程实例：某天然气集输管线，管径DN300，长度15km，起点压力6.0MPa，终点压力5.2MPa，天然气相对密度0.65，平均温度15℃，计算输气能力约85×10⁴m³/d。1.2潘汉德尔A式（PanhandleAEquation）——输气管道流量计算潘汉德尔A式主要用于中高压输气管道的水力设计。Q=435.87×符号意义：Q—气体流量，标准状态下m³/d435.87—经验常数其余符号同威莫斯公式应用场景：中等管径（300mm\800mm）、中等压力（3.5MPa\7MPa）的输气管道。限制条件：适用于雷诺数在5×10⁶~1.5×10⁷之间的湍流区，不适用于小管径低压管道。1.3潘汉德尔B式（PanhandleBEquation）潘汉德尔B式（修正式）适用于高压大流量输气管道。Q=737×符号意义：同潘汉德尔A式。应用场景：大管径（≥600mm）、高压（5MPa~12MPa）的长距离输气干线。限制条件：粗糙度影响未显式体现，雷诺数适用范围为4×10⁶~4×10⁷。1.4AGA（美国天然气协会）公式AGA公式基于摩阻系数的严格分析，是目前被认为最准确的输气管道水力计算方法之一。首先使用科尔布鲁克-怀特方程（见下文1.6节）确定摩阻系数f，然后代入基本输气方程计算流量：Q=77.54×或P符号意义：Q—气体流量，m³/d77.54—包含单位换算的常数f—摩阻系数，无量纲（Colebrook-White方程求解）D—管道内径，mm其余符号同前应用场景：各类输气管道的精确水力设计。限制条件：需要迭代求解摩阻系数；计算相对复杂，需编程或使用专用软件。1.5IGT（气体技术研究所）公式Q=0.001149×应用场景：城市燃气管网中压、低压管道的水力计算。1.6科尔布鲁克-怀特（Colebrook-White）方程——摩阻系数计算摩阻系数f的隐式方程，是湍流管道摩阻计算的理论基础。1符号意义：f—达西-威斯巴赫摩阻系数，无量纲ε—管道绝对粗糙度，mm（钢管一般取0.046mm~0.2mm）D—管道内径，mmRe—雷诺数，无量纲应用场景：所有湍流管道流动的摩阻系数精确计算。限制条件：这是关于f的隐式方程，需要迭代求解或使用近似显式公式。1.7海兰德（Haaland）显式摩阻系数公式隐式科尔布鲁克-怀特方程的高精度显式近似。1符号意义：同Colebrook-White方程。应用场景：避免迭代求解的管道摩擦系数快速计算。1.8穆迪（Moody）摩阻系数公式——显式近似f=0.0055[1+应用场景：工程快速估算。1.9苏联天然气研究所早期公式Q=k值随管径和流态变化。1.10一般水平管流方程基于达西-威斯巴赫公式的能量方程，可推导出一般水平管流方程：P1.11高山管道水力计算校正——高程校正因子当管道沿线存在高程变化时，需引入高程校正因子。苏联公式法采用如下公式进行高程校正：Q符号意义：Qc—校正后流量Δh—管道起终点高程差C—经验系数1.12达西-威斯巴赫（Darcy-Weisbach）摩擦损失公式Δ符号意义：ΔP—管道摩擦压力降，Paf—达西摩擦系数L—管道长度，mD—管道内径，mρ—气体密度，kg/m³v—气体流速，m/s应用场景：所有流体（气体、液体）管道摩擦损失计算。

第二部分：雷诺数与流态判别公式2.1雷诺数定义Re符号意义：Re—雷诺数，无量纲ρ—流体密度，kg/m³v—平均流速，m/sD—管道内径，mμ—动力粘度，Pa・sν—运动粘度（ν=μ/ρ），m²/s工程实例：天然气管道，D=0.3m，v=10m/s，ρ=40kg/m³（压力约5MPa），μ≈1.2×10⁻⁵Pa・s，Re≈10⁷（湍流）。2.2流态判别标准Re<2000：层流2000≤Re≤4000：过渡流Re>4000：湍流2.3天然气动力粘度计算——Sutherland公式μ=符号意义：μ—温度T下的动力粘度，Pa・sμ₀—参考温度T₀下的动力粘度T—绝对温度，KT₀—参考温度，KC—Sutherland常数（甲烷取164K）2.4天然气运动粘度ν=

第三部分：天然气物性计算公式3.1理想气体状态方程PV=nRT

或

符号意义：P—绝对压力，PaV—体积，m³n—摩尔数，molR—通用气体常数，8.314J/(mol・K)T—绝对温度，Kρ—密度，kg/m³M—摩尔质量，kg/mol限制条件：仅适用于低压（≤1MPa）、高温条件，高压下误差较大。3.2实际气体状态方程——压缩因子修正PV=ZnRT符号意义：Z—气体压缩因子，无量纲其余符号同理想气体3.3CNGA压缩因子公式美国加州天然气协会经验公式，是AGA计算方法的简化版：Z=其中：Pavg为平均压力，kPa；Tf为平均温度，K；G为气体相对密度。适用范围：P≤7MPa，T≥0℃，G≤0.7。3.4AGA8号报告压缩因子计算AGA8-DC92是国际标准ISO12213-2推荐的天然气压缩因子计算方法。详细形式为外延式二维维里方程：Z=1+其中：

B符号意义：Z—压缩因子Bm(T)、Cm(T)—第二、第三维里系数xi、xj、xk—组分摩尔分数ρ—摩尔密度应用场景：天然气压缩因子的高精度计算。3.5彭-鲁宾逊（Peng-Robinson）状态方程——PR方程P=符号意义：P—压力Vm—摩尔体积a、b—物性参数α—与温度和偏心因子有关的函数3.6索阿韦-雷德利希-邝（Soave-Redlich-Kwong）状态方程——SRK方程P=3.7甘内西・戴维斯（Gunn-Yamada）通用状态方程由状态方程计算热力学函数。3.8本尼迪克特-韦布-鲁宾（Benedict-Webb-Rubin）状态方程——BWR方程多常数状态方程，用于高压天然气计算：P=ρRT+(3.9GERG-2008方程GERG-2008是国际标准ISO20765-2推荐的天然气热物性参数计算方法，是当前最精确的天然气状态方程之一。应用场景：LNG工艺设计、高压天然气物性计算、烃露点计算等。3.10气体相对密度G=标准参比条件：0℃，101.325kPa；或15℃，101.325kPa。3.11天然气密度（气体混合物）ρ=符号意义：yi—组分i的摩尔分数Mi—组分i的摩尔质量，g/mol3.12标准状态体积换算V下标N表示标准状态。3.13天然气热值计算（ISO6976:2016）高位发热量：

H低位发热量：

H3.14沃泊指数（华白指数）——WobbeIndex高位沃泊指数：

W低位沃泊指数：

W符号意义：H_G—高位发热量，MJ/m³H_N—低位发热量，MJ/m³G—气体相对密度（空气=1）应用场景：燃气互换性评价——具有相同沃泊指数的不同燃气组分在相同燃烧压力下可释放同等热负荷。

第四部分：管道局部阻力与当量长度公式4.1局部阻力系数法hf=ζv符号意义：hf—局部阻力引起的压头损失，m液柱ζ—局部阻力系数，无量纲v—局部截面处的平均流速，m/sg—重力加速度，9.81m/s²ΔP—局部压力降，Paρ—气体密度，kg/m³4.2当量长度法L符号意义：Leq—当量长度，mζ—局部阻力系数D—管道内径，mλ—管道沿程阻力系数应用场景：将管件的局部阻力折算成相同管径的等长度直管段进行水力计算。4.3总折算长度L4.4弯头局部阻力系数——经验公式90°标准弯头：

ζ=0.51.54.5突然扩大局部阻力系数ζ=4.6突然缩小局部阻力系数ζ=0.5(1-符号意义：A1—上游截面积，m²A2—下游截面积，m²

第五部分：管道传热计算公式5.1天然气管道沿线温度分布——苏霍夫公式T符号意义：Tx—距起点x处气体温度，℃T₀—环境温度（埋地管道为土壤温度），℃Tin—管道起点气体温度，℃D—管道外径，mK—从气体到环境的总传热系数，W/(m²・℃)x—距起点的距离，mm—气体质量流量，kg/sCp—气体定压比热容，J/(kg・℃)应用场景：埋地管道和架空管道沿线温度分布计算。5.2气体平均温度计算T更精确的计算需考虑对数平均温差。5.3总传热系数计算组合传热的热阻关系式——埋地管道包括管内对流、管壁导热、绝缘层导热以及土壤导热等：15.4焦耳-汤姆逊效应——等焓节流温降μ工程简化：Δ应用场景：调压站、节流元件处气体压力降低造成的温度下降。5.5热平衡方程m5.6埋地管道土壤热阻计算R

第六部分：天然气压缩机计算公式6.1压缩机压比ε=符号意义：ε—压缩比（压比），无量纲Pd—排气压力，MPaPs—进气压力，MPa6.2等熵压缩过程绝热指数（κ=cp/cv）：

P6.3等熵压缩功（多变压头）H符号意义：His—等熵压头，J/kgκ—等熵指数（cp/cv）Zs—入口气体压缩因子R—气体常数，J/(kg・K)Ts—进气温度，Kε—压缩比6.4多变压缩过程P6.5多变压缩功（多变头）H符号意义：Hpol—多变压头，J/kgn—多变指数6.6压缩机功率计算等熵功率：

P轴功率：

P符号意义：m—质量流量，kg/sηpol—多变效率（离心式压缩机通常取0.8~0.86）6.7压缩机的总效率关系η6.8压缩机排气温度等熵过程排气温度：

T实际过程排气温度：

T6.9压缩机站的基本方程——特性曲线压缩比ε、多变效率ηn、功率/压头H与流量Q和转速n的关系由特性曲线决定。6.10压缩机相似定律相同气体：Q∝n；H∝n²；P∝n³6.11欧拉方程——透平机械基本方程H=符号意义：H—理论压头u—圆周速度cu—绝对速度的圆周分量6.12连续方程ρ6.13能量方程（稳定流动）h

第七部分：孔板流量计公式7.1孔板流量计基本公式qq符号意义：qm—质量流量，kg/sqv—体积流量，m³/sC—流出系数，无量纲β—直径比β=d/D，无量纲d—节流孔直径，mD—上游管道内径，mρ—流体密度，kg/m³ΔP—孔板前后差压，Pa7.2斯特尔茨（Stolz）流出系数公式（ISO5167）C=0.5961+0.02617.3可膨胀性系数（用于气体）ε=1-(0.351+0.2567.4天然气超压缩因子Fz（SY/T6143标准）天然气超压缩因子Fz是气体偏离理想气体性质的修正参数。7.5实际天然气孔板流量计算q应用场景：天然气计量站、输气管道流量在线计量。7.6节流装置流量计算的通式qα包含流出系数和可膨胀性系数。

第八部分：燃气调压阀计算8.1阀门流量系数Kv值计算K8.2阻塞流判别条件P8.3气体阀门Kv值计算——阻塞流K8.4压力恢复系数FLF8.5气体调节阀Kv通用计算公式K符号意义：Q—气体流量，Nm³/hG—气体相对密度T—气体温度，KFy—流体阻塞流系数

第九部分：音速（临界流）喷嘴——临界流系数9.1临界流条件P应用场景：确定喷嘴是否达到临界流状态。9.2音速喷嘴质量流量公式q9.3临界流系数C*（CriticalFlowFactor）C应用场景：天然气流量标准装置——依据音速喷嘴独特的临界流测量原理，用临界流函数表征喉部流体特性。

第十部分：分离器设计计算公式10.1液滴沉降速度公式——重力沉降W符号意义：Wo—液滴沉降速度，m/sg—重力加速度，9.81m/s²dL—液滴直径，m（通常取60×10⁻⁶~100×10⁻⁶m）ρL—液体密度，kg/m³ρG—气体密度，kg/m³CD—阻力系数，无量纲（Re的函数）应用场景：重力式气液分离器设计，确定液滴能否从气相中分离。10.2小液滴沉降速度——斯托克斯定律W10.3丝网捕雾器最大允许速度V符号意义：Vmax—气体通过丝网的最大允许速度，m/sk—经验系数（0.03~0.13之间）应用场景：分离器出口丝网捕雾段设计。10.4气体通过丝网捕雾器设计速度V10.5立式分离器气体处理能力计算Q=10.6立式分离器设计直径D=其中：

A10.7分离效率计算方法液滴在分离器中的分离效率是液滴直径和气体流速的函数，工程上常用粒级效率模型进行计算。10.8阿基米德数Ar（用于沉降速度计算）Ar=根据Ar数查表求Re和ω（沉降速度）。10.9分离器临界气流速度当气流速度低于临界值Vt时，所需分离直径的液滴即可从气相中分离。

第十一部分：三甘醇（TEG）脱水塔计算11.1吸收塔塔径计算——贝恩-霍根（Bain-Haugen）关联式Δ11.2通用压降关联图法（GPDC）——Sherwood-Eckert关联式Δ11.3塔负荷系数法D=符号意义：D—塔径，mQG—气体体积流量，m³/sC—塔负荷系数，m/s11.4贫甘醇溶液最低浓度确定由气液平衡关系确定欲达到干气平衡露点所必须的贫三甘醇溶液最低浓度。11.5甘醇循环量估算L11.6传质单元数（NTU）方法NTU应用场景：确定脱水塔填料层高度。11.7填料层高度计算Z=其中HTU为传质单元高度。11.8吸收塔塔板数计算——克雷姆塞尔（Kremser）方程y符号意义：A—吸收因子（L/mG）N—理论板数

第十二部分：胺法脱硫（脱CO₂）吸收塔计算12.1填料吸收塔体积总传质系数KGaK12.2化学吸收的传质速率方程N符号意义：NA—组分A的传质通量，mol/(m²・s)kG—气相传质系数kL—液相传质系数12.3双膜理论——气液传质模型1符号意义：KG—总气相传质系数H—亨利常数12.4浮阀塔塔板数确定按净化天然气质量标准和对CO₂吸收率的要求经计算确定吸收塔板数。12.5胺溶液循环量计算L

第十三部分：液化天然气（LNG）工艺计算公式13.1相平衡计算——闪蒸计算给定压力和温度，通过相平衡计算获得汽相和液相组成分布。13.2露点压力/露点温度计算给定温度求泡点压力；给定压力求露点温度：i=1i=1符号意义：Ki—气液平衡常数Ki=yi/xiN—组分总数13.3液化率（闪蒸率）V由闪蒸计算确定。13.4液化冷箱热负荷Q=∑13.5PR方程在LNG工艺中的应用用于计算烃类及混合制冷剂（如丙烷-乙烯）的相平衡特性。13.6丙烷预冷混合制冷剂流程分析需计算原料气压力、组成和预冷温度等参数对能耗的影响。13.7LNG密度计算ρ准确预测LNG密度的状态方程包括PR方程、SRK方程、VT-SRK等。

第十四部分：天然气气质（热物性）计算公式14.1天然气高位发热量计算H14.2天然气低位发热量计算H14.3相对密度计算ρ14.4不同参比条件之间的换算系数——ISO13443用于在不同基准温度压力下换算天然气的体积、密度、相对密度、压缩因子、发热量和沃泊指数。14.5气体混合法则——凯氏法则P14.6道耳顿分压定律——理想混合气体P14.7阿马格分体积定律V14.8比热容计算——混合气体C14.9气体等熵指数κ=应用场景：压缩机设计、音速喷嘴计算、可膨胀性系数计算。14.10水露点与含水量计算根据经验公式从含水量推算天然气在给定压力下的水露点温度。14.11烃露点计算由组成数据通过相平衡计算获得天然气烃露点。要求在输送过程中，烃露点必须比沿线各地段的最低地温低5℃。14.12天然气的临界流系数——AGA8号报告计算由AGA8状态方程计算临界流系数，用于音速喷嘴流量计量。

第十五部分：防火安全与泄放计算15.1管道泄放（吹扫）流量计算根据管道容积和允许时间确定。15.2安全阀泄放面积计算（气体）A=符号意义：A—最小泄放面积，mm²Q—最大泄放流量，kg/hC—气体特性系数Kd—安全阀额定泄放系数M—气体摩尔质量15.3火炬系统