石油和天然气工业中的热力学过程：液位、体积和压力变化的计算

在石油和天然气工业中，与温度变化及其对流体和结构的影响相关的精确计算起着关键作用。本文讨论了三个实际任务，演示温度波动如何影响:

油箱中的油位,
汽油冷却时的体积,
温度变化时管路中的压力。

这些任务由一个共同的概念统一起来：物质和材料的热膨胀（或压缩），以及其对工程系统的影响。在每种情况下，热膨胀系数，弹性模量和热力学定律被用来预测液体和结构的行为。

让我们用一个脚本开始我们的分析，用于计算温度波动期间油箱油位的变化。

# 输入数据
t1 = 20.0   # 白天温度，°C
t2 = 10.0   # 夜间温度，°C
H = 6.0     # 油箱高度，m
ρ20 = 850.0 # 20°C时的油密度，kg/m3
β = 0.0007  # 油的热膨胀系数，1/°C

# 计算电平变化
Δt = t1 - t2
h = H * β * Δt

# 用公式输出结果
println("任务1:")
println("公式：h=H*β*(t1-t2)")
println("取代：h=△H*△β*(△t1-△t2)")
println("计算：h=△H*△β*△(Δt)= ", H*β*Δt)
println("\油位变化h= ", round(h, digits=4), " m")
println("液位变化不依赖于罐的直径，如")
println("体积膨胀仅影响液柱的高度。")

Задача 1:
Формула: h = H * β * (t1 - t2)
Подстановка: h = 6.0 * 0.0007 * (20.0 - 10.0)
Расчет: h = 6.0 * 0.0007 * 10.0 = 0.041999999999999996

Изменение уровня нефти h = 0.042 м
Изменение уровня не зависит от диаметра резервуара, так как
объемное расширение влияет только на высоту столба жидкости.

物理解释

-油加热时膨胀，冷却时收缩。
-系数 β = 0.0007 1/°C 这意味着当在1°C下加热时，油的体积会增加 0.07%.
-在这种情况下，温差（Δt = 10°C)导致4.2cm的水平变化。

该计算演示了液体的简单热力学属性（热膨胀）如何影响储罐的操作参数。为了进行准确的库存监测，必须考虑到这些变化，特别是在每日温度波动较大的地区。

现在考虑下面的例子：冷却过程中汽油体积的变化。该脚本计算汽油体积的变化，因为它从运输温度冷却到地面温度。使用两步计算:

调整音量至标准条件(20°C),
体积到地面温度的转换。

# 输入数据
t = 25.0    # 运输温度，°C
t1 = 15.0   # 地面温度，°C
V = 10.0    # 运输过程中的体积，m3
ρ20 = 730.0 # 20°C时的汽油密度，kg/m3
β = 0.0012  # 汽油热膨胀系数，1/°C

# 20°C下的体积计算
V20 = V / (1 + β * (t - 20))

# 地面温度下体积的计算
V1 = V20 * (1 + β * (t1 - 20))

# 用公式输出结果
println("\n任务2:")
println("1. 找到20°C的体积:")
println("   V20 = V / (1 + β * (t - 20))")
println("   V20 = $V / (1 + $β * ($t - 20)) = ", round(V20, digits=4), " m3")

println("2. 我们找到地面温度下的体积:")
println("   V1 = V20 * (1 + β * (t1 - 20))")
println("   V1 = $(round(V20, digits=4)) * (1 + $β * ($t1 - 20))")

println("\冷却后汽油的最小量: ", round(V1, digits=2), " m3")

Задача 2:
1. Находим объем при 20°C:
   V20 = V / (1 + β * (t - 20))
   V20 = 10.0 / (1 + 0.0012 * (25.0 - 20)) = 9.9404 м³

2. Находим объем при температуре грунта:
   V1 = V20 * (1 + β * (t1 - 20))
   V1 = 9.9404 * (1 + 0.0012 * (15.0 - 20))

Итоговый объем бензина после остывания: 9.88 м³

冷却后，汽油体积从10.0m3降至9.88m3。

物理解释

汽油具有显着的膨胀系数（β=0.0012 1/°C），比第一任务的油高1.7倍。
在10°C（从25°c到15°C）下冷却可使体积减少1.2％。
使用20°C的参考温度进行计算可确保精度，因为密度p20是针对这些条件精确设置的。

实际结论

对于商业会计，有必要考虑到温度校正–1.2％的差异对于大量而言是显着的。
在设计储罐时，重要的是提供体积的"热储备"。
降低到20°C的计算方法是炼油行业的标准做法。

特点：与第一项任务不同，在这里，数量的变化直接影响到可销售产品的数量，这对经济具有直接的重要性。

现在让我们来看看如何计算燃料冷却过程中管道中的压力变化。

该脚本模拟当所输送燃料的温度降低时封闭管道中的压力变化。联合影响被考虑在内

*燃料的热压缩,

*管道材料的热压缩,

*管液系统的弹性特性。

# 输入数据
d0 = 0.514       # 管道直径，m
δ = 0.008        # 壁厚，m
E = 2e11         # 杨氏模块，Pa
αT = 3.3e-5      # 金属热膨胀系数，1/°C
ρ20 = 840.0      # 燃料密度，kg/m3
K = 1.5e9        # 燃料的弹性模量，Pa
py = 5e3         # 蒸汽弹性，Pa
p1 = 2.5e6       # 初始压力，Pa
t1 = 10.0        # 初始温度，°C
t2 = 7.0         # 最终温度，°C
βT = 0.0007      # 燃料热膨胀系数，1/°C

# 更改的计算
Δt = t2 - t1

println("\n任务3:")
println("1. 改变管道容积:")
println("   ΔVT/V = 3 * αT * Δt")
println("   ΔVT/V = 3 * $αT * ($t2 - $t1) = ", 3*αT*Δt)

# 改变管道容积
ΔVT_V = 3 * αT * Δt

println("2. 燃料体积变化:")
println("   ΔVf/V = βT * Δt")
println("   ΔVf/V = $βT * ($t2 - $t1) = ", βT*Δt)

# 改变燃料量
ΔVf_V = βT * Δt

println("\n3． 总成交量变化:")
println("   ΔV/V = ΔVT/V + ΔVf/V")
println("   ΔV/V = $(ΔVT_V) + $(ΔVf_V) = ", ΔVT_V + ΔVf_V)

# 总成交量变化
ΔV_V = ΔVT_V + ΔVf_V

println("4. 体系的有效弹性模量:")
println("   KT = (E * δ) / d0 = ($E * $δ) / $d0 = ", (E*δ)/d0)
KT = (E * δ) / d0
println("   1/K_eff = 1/K + 1/KT = 1/$K + 1/$KT")
K_eff = 1 / (1/K + 1/KT)
println("   K_eff = ", K_eff)

println("\n5． 压力变化:")
println("   Δp = -K_eff * ΔV/V")
println("   Δp = -$K_eff * $(ΔV_V) = ", -K_eff*ΔV_V)
Δp = -K_eff * ΔV_V

println("6. 最终压力:")
println("   p2 = p1 + Δp = $p1 + $(Δp) = ", p1 + Δp)
p2 = p1 + Δp

# 转换为MPa和atm
p2_MPa = p2 / 1e6
p2_atm = p2 / 101325

println("\结果:")
println("温度＜t2°C时的压力:")
println("  ", round(p2_MPa, digits=2), " MPa,MPa")
println("  ", round(p2_atm, digits=2), " 自动柜员机")

Задача 3:
1. Изменение объема трубопровода:
   ΔVT/V = 3 * αT * Δt
   ΔVT/V = 3 * 3.3e-5 * (7.0 - 10.0) = -0.000297

2. Изменение объема топлива:
   ΔVf/V = βT * Δt
   ΔVf/V = 0.0007 * (7.0 - 10.0) = -0.0021

3. Суммарное изменение объема:
   ΔV/V = ΔVT/V + ΔVf/V
   ΔV/V = -0.000297 + -0.0021 = -0.002397

4. Эффективный модуль упругости системы:
   KT = (E * δ) / d0 = (2.0e11 * 0.008) / 0.514 = 3.1128404669260697e9
   1/K_eff = 1/K + 1/KT = 1/1.5e9 + 1/3.1128404669260697e9
   K_eff = 1.0122311261071277e9

5. Изменение давления:
   Δp = -K_eff * ΔV/V
   Δp = -1.0122311261071277e9 * -0.002397 = 2.426318009278785e6

6. Итоговое давление:
   p2 = p1 + Δp = 2.5e6 + 2.426318009278785e6 = 4.926318009278785e6

Результат:
Давление при температуре 7.0 °C:
  4.93 МПа
  48.62 атм

物理解释

在3°C冷却时:

*管道压缩0.03%,

*燃料压缩0.21%
压力增加初始值的97%。
主要贡献（87％）提供燃料压缩

重要提示：获得的压力（4.92MPa）接近主管道的标准操作压力（5-7MPa），这证实了这种计算的必要性。

结论

所考虑的任务清楚地表明，即使是微小的温度变化也会影响石油和天然气行业的关键参数。

-油箱油位的变化取决于温差和膨胀系数，但不取决于油箱的几何形状。
-在冷却过程中汽油的体积被调整考虑到其密度和温度系数。
-管道中的压力不仅取决于燃料的特性，还取决于管壁的弹性。

这些计算有助于防止事故发生，优化存储和运输，并设计更可靠的系统。