7 结束语
传统的蒸汽供热管网的热耗散计算是以热电厂出来的总热量减去用户得到的总热量之和，在可测量正确的假定下，也无法区分热耗散的性质。这种方法无法对整个管网的供热进行可靠的监测，更无法对整个网络的各个部分工作状态进行完整的评估。
针对蒸汽供热管网中热损耗现象所建立的正常管网散热的热损耗模型具有软仪表的测量模型特质，基于机理分析并依赖于可测信息集的模型。
本文所提出的热损耗模型和递推算法并不要求模型的精度特别高，它主要是用来研究管网能量损耗的，这相对于时下热电厂使用的质量计量有了一个新的突破，使用有用能来计量可以相对准确的分析管网的运行工况。研究编写的软件可对某个因素对管网的影响程度进行分析，通过对大量数据的计算、比对，就可以得到相对能量损耗较小的管线运行参数。
本文研究管网运行模型时考虑到的因素还不是很全，这就有待于进一步研究，例如，模型的结构并不是固定的，不同的模型得到的研究结果自然也是不同的。由于毕业设计时间有限，无法把所有的可能性都考虑周全，希望在今后的学习和研究中进一步改善。
30
参考文献
[1] 卢峰，田贯三，丁国玉，孙永海. 蒸汽供热管网凝结水损失的影响因素及其分析[J]，实用
节能技术，2008，23（10）：50-54
[2] 黄镇海,朱云,钟绍俊,周秀英．蒸汽管网的智能监测系统[J]，中国计量学院学报，2005，
16(2)：127-130
[3] 李瑜,崔国民,涂惟民,姜涛．蒸汽管网的建模和设计优化[J]，节能，2009，10(4)：14-17 [4] 张增刚，李继志，李永安.城市蒸汽管网水力热力耦合系统及其应用[J],微计算机信息，
2008，6（3）：207-209
[5] 姚庆，翁红芳．基于管网模拟技术的蒸汽管网热损分析方法[J]，数字石油和化工，2008，
(2)：47-49
[6] 吴建恒. 热电厂供气管网的节能改造[J]，实用节能技术，2009，（2）：60-63
[7] 刘大山. 不容忽视的节能领域--蒸汽管网系统[J],中国机械工程，2002，13（19）：1634-1636 [8] 李世武. 蒸汽管网运行工况模拟计算方法[J],煤气与热力，2001，21（6）：552-556 [9] 曹燕青，陈志圭．用于蒸汽管网系统的严格在线模拟与智能监测技术[J]，计算机与应用化
学，2006，25(2)：6-8
[10] 赵家军,陈贵军,刘晓华,王旭东. 蒸汽供热管网运行工况的水力模拟计算与软件编制[J]，
节能，2006，9(2)：34-35
[11] 章金华．室外供热管网的水力计算[J]，煤气与热力，2001，35(1)：162-163
[12] 陈志奎，曹雁青. 智能监测技术在蒸汽管网系统中的应用[J]，区域供热，2005，4：1-10 [13] 邓永成. 蒸汽管网实时数据采集与计算机模拟优化系统[J],能源技术，2005，26（1）:14-17 [14] 李祥立，邹平华. 分段积分法蒸汽管网水力计算[J],煤气与热力，2005，25(3):48-51 [15] 刘笑驰，曹瑞忠，吕崇德.大型蒸汽管网的在线仿真研究[J],系统仿真学报，2002,14(3):397-402 [16] 孙玉宝.蒸汽管网水力热力联合计算数学模型及应用方法的研究[J],2006,(10):242-246 [17] J. Okrajni ,K. Mutwil , M. Ciesla .Steam pipelines’effort and durability [J]，
Journal of Achievements in Materials and Manufacturing Engineering，2007，22（2）：63-66
[18] King Mongkut. Drying Parawood with Superheated Steam [J]，Journal of Achievements
in Materials and Manufacturing Engineering，2007，4(4)：215-219
[19] M.Thirumarimurugan,T.Kannadasan and E.Ramasamy. Performance Analysis of Shell and
Tube Heat Exchanger Using Miscible System[J]，American Journal of Applied Sciences，2008，5(5)：548-552
31
附录 源程序
Public P As Double, T1 As Double, D As Double, L As Integer, α As Double, jrchd As Double, ω As Double, hjwd As Double, M As Double, jinya As Double
Public Nu As Double, CP As Double, ρ As Double, V As Double, PRN As Double, RAMD As Double, Re As Double, U As Double
Public Function ckt(ByVal jinya, ByVal T1, ByVal D, ByVal L, ByVal jrchd, ByVal ω, ByVal hjwd, ByVal M, ByVal α) As Double Dim p2(20) As Double Dim A As Double
Dim ls As Double ' 蒸汽流速 Dim dp As Double Dim Tn As Double Dim sanre As Double Dim β As Double Dim Tw As Double Dim Q As Double Dim wj As Double Dim Tf2 As Double Dim N As Integer Dim KF As Double Dim bz As Double Dim rz As Double Dim nb As Double Dim wb As Double
Dim T2(20) As Double
Dim tj(1 To 10) As Integer Do
If D > 600 Or D < 100 Then
MsgBox (\请检查内径数据的正确性\ tj(1) = 0
Exit Function Else: tj(1) = 1 End If
If T1 > 1000 Or T1 < 100 Then
MsgBox (\请检查进口温度输入的正确性\ tj(2) = 0
Exit Function
32
Else: tj(2) = 1 End If
If L > 10000 Or L <= 0 Then
MsgBox (\请检查输入长度的正确性\tj(3) = 0
Exit Function Else: tj(3) = 1 End If
If jinya > 1.5 Or jinya < 0.1 Then MsgBox (\请检查压力输入的正确性\tj(4) = 0
Exit Function Else: tj(4) = 1 End If
If M > 100 Or M < 5 Then
MsgBox (\请检查质量流量输入的正确性\tj(5) = 0
Exit Function Else: tj(5) = 1 End If
If hjwd > 60 Or hjwd <= 0 Then
MsgBox (\请检查环境温度输入的正确性\tj(6) = 0
Exit Function Else: tj(6) = 1 End If
If jrchd > 100 Or jrchd < 0 Then
MsgBox (\请检查绝热层厚度输入的正确性\tj(7) = 0
Exit Function Else: tj(7) = 1 End If
If ω > 1.5 Or ω <= 0 Then
MsgBox (\请检查绝热层导热系数的正确性\tj(8) = 0
Exit Function Else: tj(8) = 1
33
End If
If α > 50 Or α < 0 Then
MsgBox (\请检查空气传热系数的正确性\tj(9) = 0
Exit Function Else: tj(9) = 1 End If
Loop Until tj(1) = 1 Or tj(2) = 1 Or tj(3) = 1 Or tj(4) = 1 Or tj(5) = 1 Or tj(6) = 1 Or tj(7) = 1 Or tj(8) = 1 Or tj(9) = 1
i = 1 N = 1
T2(1) = T1 - 0.04 * L 'jinya = Val(Text4.Text) p2(1) = jinya Do
P = jinya
T = (T1 + T2(N)) / 2
Rem 已知压力(MPa)和温度(℃)，求比容(m^3/kg) Call PT2V97(P, T, V, r)
Rem 已知压力(MPa)和温度(℃)，求定压比热(kJ/(kg.℃)) Call PT2CP97(P, T, CP, r)
Rem 已知压力(MPa)和温度(℃)，求运动粘度(m^2/s) Call PT2U97(P, T, U, r)
Rem 已知压力(MPa)和温度(℃)，求普朗特数 Call PT2PRN97(P, T, PRN, r)
Rem 已知压力(MPa)和温度(℃)，求热传导系数(W/(m.℃))
Re = 4 * M * 1000 * 1000 * V / (3.14159 * D * U * 3600) '计算雷诺数 A = 0.0056 + 0.5 / Re ^ 0.32 ρ = 1 / V
u1 = 4 * M * 1000 / 3600 / 3.14 / D ^ 2 * 1000000 * V
dp = A * L * ρ * u1 ^ 2 / D * 1000 / 1000000 '计算压差 i = i + 1
p2(i) = jinya - dp
Loop Until Abs(p2(i) - p2(i - 1)) <= 0.001 Do
P = (jinya + p2(i)) / 2 T = (T1 + T2(N)) / 2
Rem 已知压力(MPa)和温度(℃)，求比容(m^3/kg) Call DT2V96(P, T, V, r)
34