汽轮机全机型热力计算技能（精准版）

适配：非再热/一次再热/背压/抽凝/抽背全机型，含IF97物性表、完整公式、单位说明、多工况模板

一、效率参数（全局固定）

> 注意：本版修正发电机效率为0.985（原技能0.99），相应ηmg=0.9653。后续所有计算以此为准。

二、机组类型判定（执行计算前必须先判定）

输入：主汽压力P₀(MPa)、有无再热、抽汽口数、是否有背压排汽
↓
【第1层】有无中间再热？
  ├─ 无再热 → 第2层
  └─ 有再热 → 再热机组流程（下文第六节）

【第2层】是否纯凝汽轮机（无抽汽、无背压）？
  ├─ 是 → 纯凝机组：单段焓降计算
  └─ 否 → 第3层

【第3层】是背压排汽还是凝汽排汽？
  ├─ 背压排汽（Pk<0.1MPa且机组标为背压/抽背）→ 背压机组
  └─ 凝汽排汽（Pk≥0.1MPa表压，或标为凝汽）→ 第4层

【第4层】是否有工业抽汽？
  ├─ 无抽汽 → 纯凝机组（单段）
  ├─ 有抽汽且抽汽口≥1 → 抽凝机组（分段膨胀）
  └─ 排汽侧有供热抽汽（抽背）→ 抽背机组（分段膨胀）

三、压力换算（强制规则）

所有蒸汽物性查询只允许绝对压力，无例外。

> 核心规则（2026-05-18修正）：
> - 除非用户明确标注"绝压"或给出真空度/负压值，一律默认表压，+0.1 MPa转绝对压力
> - 背压机组的背压（如0.98 MPa）同样按表压处理 → 绝压 = 0.98 + 0.1 = 1.08 MPa(a)
> - 只有凝汽机组排汽（7~15 kPa）因数值极小，直接作为绝对压力使用
> - 示例：8.83 MPa表压 = 8.93 MPa(a)；3.03 MPa表压 = 3.13 MPa(a)

四、内效率选用规则

4.1 非再热机组（整机等效内效率，不分缸）

> 修正说明：本版将抽凝/供热/背压/抽背合并一档（原技能背压列的0.81/0.83/0.84），并与纯凝分开。原技能将"背压/抽背/供热"三档与纯凝混列，逻辑不清晰，本版修正。

4.2 一次再热三缸机组（分缸取效率）

五、IAPWS-IF97 蒸汽物性参考表

> 使用规则：所有焓/熵值通过双线性插值查取。表中相邻节点间线性插值；超出表范围则外推（注明风险）。

5.1 主蒸汽区（高温高压，0.1~16 MPa，250~560 °C）

> 格式说明：数值 = h (kJ/kg) / s (kJ/kg·K)，两值以"/"分隔
> 例：8.0MPa/500°C → h=3399.5, s=6.7266

5.2 低中压区（0.7~4.0 MPa，200~450 °C）

5.3 抽汽常用压力级（1.5~4.0 MPa，精确覆盖）

5.4 IAPWS-IF97 焓值计算（iapws库，2026-05-18修正）

> 核心原则：禁止hardcode任何固定焓值字典。所有焓值必须通过IAPWS-IF97标准函数根据实际压力+温度实时计算。
>
> 工具：pip install iapws → from iapws.iapws97 import IAPWS97_PT
>
> Python计算示例：
> python > from iapws.iapws97 import IAPWS97_PT > > def get_h(P_MPa_abs, t_C): > """P: 绝对压力(MPa)，t: 温度(°C)，返回焓h(kJ/kg)""" > return IAPWS97_PT(P_MPa_abs, t_C + 273.15).h > > # 用法示例（8.83表压→8.93绝压，535°C） > h0 = get_h(P0_a, t0) # 主汽 > h1 = get_h(P1_a, t1) # 第1抽汽 > h2 = get_h(P2_a, t2) # 第2抽汽 > hk = get_h(Pk_a, tk) # 背压排汽 >
>
> 工程验证（iapws vs 图纸参考，偏差均在±3 kJ/kg内）：
>
> | 节点 | P(a) | T(°C) | h_IF97 | 图纸参考 | 偏差 |
> |------|------|--------|--------|---------|------|
> | 主汽 h0 | 8.93 | 535.0 | 3475.59 | 3476.6 | -1.01 |
> | 第1抽汽 h1 | 3.13 | 401.3 | 3232.34 | 3234.1 | -1.76 |
> | 第2抽汽 h2 | 1.82 | 315.4 | 3064.30 | 3066.9 | -2.60 |
> | 背压 hk | 1.08 | 300.0 | 3049.59 | 3049.56 | +0.03 |
>
> 注意：IF97是标准函数值，厂家图纸是修正值，偏差属正常范围。不再使用固定焓值字典。

5.5 饱和水与饱和蒸汽参数（常用压力级）

六、全节点焓值取用优先级

> 特别说明：实际工程中抽汽口温度往往由现场实测或设计给出，与等熵膨胀终点存在偏差（节流损失、散热损失等），此时以给定值为准，在功率方程中强制使用给定抽汽焓。

七、计算公式（完整+单位明确）

7.1 焓降计算

有实测焓值：Δh_实测 = h_上游 - h_下游
等熵膨胀：Δhs = h_上游 - h_s（等熵）
实际膨胀：Δh = ηᵢ × Δhs
实际出口焓：h_下游 = h_上游 - Δh

7.2 多级抽汽功率方程（单位：kg/h，MW）

原始方程：

D₀×Δh₁ + (D₀-D₁)×Δh₂ + (D₀-D₁-D₂)×Δh₃ + ... = Pi(kW) × 3600

最终形式（推荐使用）：

Pi(MW) = [D₀×Δh₁ + (D₀-D₁×1000)×Δh₂ + ...] / 3600   （D单位：kg/h）
Pe(MW) = Pi × ηmg

> 单位说明：D₀等为kg/h，Δh为kJ/kg，/3600后得到MW。若D以t/h输入，则分子需乘1000：
> Pi(MW) = [D₀(t/h)×1000×Δh₁ + (D₀-D₁)×1000×Δh₂ + ...] / 3600

7.3 已知Pe反求D₀（一元一次方程）

设 Pi = Pe × 1000 / ηmg   （kW）
代入功率方程：
Pi × 3600 = D₀×Δh₁ + (D₀-D₁)×Δh₂ + ...
展开：Pi×3600 = D₀×(Δh₁+Δh₂+...) - D₁×Δh₂ - D₂×Δh₃ - ...
整理：D₀ = [Pi×3600 + Σ(各抽汽流量×对应后续焓降×1000)] / Σ(各级焓降)

即：D₀ = [Pi×3600 + D₁×1000×(Δh₂+Δh₃+...) + D₂×1000×(Δh₃+...) + ...] / (Δh₁+Δh₂+Δh₃+...)

单抽汽特例：D₀ = [Pi×3600 + D₁×1000×Δh₂] / (Δh₁+Δh₂)   [kg/h]

7.4 汽耗率

d = D₀(kg/h) / [Pe(MW)×1000] × (1/1000) = D₀(t/h) / Pe(MW)   [t/(MW·h)]

7.5 热耗率

非再热机组：

R = d(t/MWh) × (h₀ - h_f@凝汽压力)   [kJ/(kW·h)]
其中 h_f@凝汽压力 = 凝汽绝对压力对应的饱和凝结水焓（查饱和水表）

一次再热机组：

总吸热量 = D₀×1000×(h₀ - h_f@P_k) + D₀×1000×(h_rh - h_HPT_out)
         = D₀×1000×[(h₀-h_f@P_k) + (h_rh - h_HPT_out)]

R = 总吸热量 / Pe(MW)   [kJ/(kW·h)]

参数说明：
  h₀       = 主汽焓（kJ/kg）
  h_f@P_k  = 凝汽压力对应饱和凝结水焓（kJ/kg）
  h_rh     = 再热后蒸汽焓（kJ/kg，直接给定或查表）
  h_HPT_out= 高压缸排汽焓（kJ/kg，需计算或给定）

抽凝/抽背机组再热热耗率（★修正公式）：
R = [D₀×1000×(h₀-h_f) + D₀×1000×(h_rh-h_HPT_out)] / Pe
  = D₀×1000/Pe × [(h₀-h_f) + (h_rh-h_HPT_out)]

> 抽凝机组特别注意：若有工业抽汽从再热前抽出，则再热蒸汽流量≠D₀，应扣除该段抽汽量。此时：
> 总吸热量 = D₀×1000×(h₀-h_f) + (D₀-Dₓ)×1000×(h_rh-h_HPT_out)

7.6 给水焓（禁止350 kJ/kg固定值）

步骤1：确认凝汽器实际排汽绝对压力 P_k
步骤2：查饱和水表得该压力下饱和液态焓 h_f
步骤3：以此h_f作为给水进水焓参与热耗计算
注意：极低真空凝汽器（P_k≈0.008~0.01 MPa），查饱和水表h_f≈173~192 kJ/kg

八、多缸分段计算规则

8.1 非再热多缸机组

统一用整机等效内效率，不拆分缸体，各自等熵膨胀过程串接。

8.2 一次再热三缸机组（★重要）

主汽 → 高压缸 → 高压缸排汽 → [再热器加热] → 再热后蒸汽 → 中压缸 → 低压缸 → 凝汽器
                  ↑
              高压抽汽D_HP         高压抽汽在 高压缸 内取，用ηHIV
                                   中压抽汽在 再热后/中压缸 内取，用ηIP
                                   低压抽汽在 低压缸 内取，用ηLP

简化算法（工程常用）：
- 中压缸+低压缸合并计算，合并段内效率取加权平均或保守取0.82~0.84
- 不影响各级抽汽做功量、总流量、功率最终结果

再热热量必须全额计入总吸热量：

Q_总 = D₀×(h₀-h_f) + D_rh×(h_rh - h_HPT_out)
其中 D_rh ≈ D₀（若高压抽汽量较小可忽略；若有显著高压抽汽则扣除）

九、工况计算模板（直接套用）

模板A：背压机组（无抽汽）

已知：P₀(表压), t₀, P_k(背压，表压), Pe(MW), ηᵢ(查表)
步骤：
  1. P₀(a)=P₀表+0.1, P_k(a)=P_k表+0.1
  2. 查表得 h₀, s₀（双线性插值）
  3. 查表求等熵排汽焓 h_ks（@P_k(a), s=s₀）
  4. Δhs = h₀ - h_ks; Δh = ηᵢ×Δhs
  5. Pi = Pe×1000/ηmg; D₀ = Pi×3600/Δh  [kg/h]
  6. 反查 h_k=h₀-Δh 对应温度 t_k
  7. 汽耗率 d = D₀(t/h)/Pe(MW)

模板B：抽凝机组（1级抽汽）

已知：P₀(表压), t₀, P₁(a), t₁, D₁(t/h), P_k(表压), Pe(MW), ηᵢ(查表)
步骤：
  1. 换算压力 → P₀(a)=P₀表+0.1, P_k(a)=P_k表+0.1
  2. 查表得 h₀, s₀（主汽节点）
  3. 查表得 h₁（抽汽节点，若给定则直接用）
  4. 【第一段：进汽→抽汽】
     等熵膨胀 @P₁(a), s=s₀ → h₁s → Δhs₁=h₀-h₁s → Δh₁=ηᵢ×Δhs₁
     实际第一段出口焓 = h₀-Δh₁
  5. 【第二段：抽汽→排汽】★关键：用抽汽口实际焓 h₁（给定值）★
     等熵膨胀 @P_k(a), s=s₀ → h_ks → Δhs₂=h₁-h_ks → Δh₂=ηᵢ×Δhs₂
     实际排汽焓 h_k = h₁ - Δh₂
  6. 流量方程：Pi=Pe×1000/ηmg → D₀=[Pi×3600+D₁×1000×Δh₂]/(Δh₁+Δh₂)  [kg/h]
  7. 校验：反算功率，与Pe比较
  8. 反查 t_k

模板C：抽凝机组（2级抽汽）

已知：P₀(表压), t₀, P₁(a), t₁, D₁(t/h), P₂(a), t₂, D₂(t/h), P_k(表压), Pe(MW), ηᵢ
步骤：
  1-4. 同模板B，计算 h₀, s₀, h₁(给定), h₂(给定), h_ks
  5. 第一段：Δh₁ = ηᵢ×(h₀-h₁s)
  6. 第二段：Δh₂ = ηᵢ×(h₁-h₂s)，其中h₂s为@P₂(a), s=s₀
  7. 第三段：Δh₃ = ηᵢ×(h₂-h_ks)
  8. D₀ = [Pi×3600 + D₁×1000×(Δh₂+Δh₃) + D₂×1000×Δh₃] / (Δh₁+Δh₂+Δh₃)
  9. 校验 + 反查排汽温度

模板D：一次再热抽凝机组（简化双段合并）

已知：P₀(表压), t₀, 再热后 t_rh, P₁(a), t₁, D₁(t/h), P_k(表压), Pe(MW)
机组类型→高压缸ηHIV, 合并段ηIC/LP
步骤：
  1. 高压缸膨胀：P₀(a)→P_HP_out(给定或估算≈0.2×P₀)
     h₀, h_HPs_out = 等熵终点 @P_HP_out, s=s₀
     h_HP_out = h₀ - ηHIV×(h₀-h_HPs_out)
     Δh_HPV = h₀ - h_HP_out
  2. 再热：h_rh 查表
  3. 中低压合并段膨胀：h_rh, s_rh → 等熵终点 @P_k, s=s_rh
     Δhs_IC = h_rh - h_ks
     Δh_IC = ηIC×Δhs_IC
     Δh₁_IC = 若有中压抽汽：Δh₁_IC = ηIC×(h_rh-h₁s)
     合并段功率 = (D₀-D₁-...-D抽)×Δh_IC
  4. 总功率 Pi = D₀×Δh_HPV/3600 + (D₀-D₁)×Δh_IC/3600
  5. Pe = Pi×ηmg → 反求D₀

十、硬性禁止条例（8条，本版新增第8条）

十一、典型工况内效率速查表

十二、常见计算错误清单（本版特有）