供热管网直埋、地沟设置及热损失分析.docxVIP

供热管网直埋、地沟设置及热损失分析 0 保温层厚度计算 通过建设部制定的总体规划超过50%的建筑的加热系统，热网输送效率应超过90%。为了减少保温层的散热损失, 提高保温层的经济性, 应采用经济厚度来计算保温层厚度, 即在考虑年折旧率的情况下, 隔热保温设施的费用和散热量价值之和为最小时的厚度。 但目前在保温层厚度的确定过程中存在不合理的因素。多数省市在制定本省市建筑节能地方标准时并没有进行分析计算, 而是完全照搬文献中给出的保温层厚度的数值。使用现有的简化模型和算法无法满足保温层厚度优化计算的需要。应分析不同地区的气候、经济等条件对管道经济厚度的影响, 以合理确定保温管道经济厚度。 1 理论分析与数值计算 1.1 定解问题的描述 目前普遍采用的直埋管道热力计算公式是20世纪30年代用热源/汇方法得到的, 土壤的温度场分布服从于熟知的舒宾公式, 将管道传热简化为单层域复连通问题。 工程计算方法对于定解问题进行了如下简化: a) 直埋热力管道在初始温度为常数的半无限大均匀介质中为线热源; b) 圆管保温层外缘处温度呈均匀分布, 为常数, 管壁面的热流强度也为常数; c) 与空气接触的大地边界的土壤温度为常数, 等于空气的温度。 这种简化会带来较大的误差, 不能满足后续演算的需要。本文采用边界离散法对具有均匀保温层的地下保温管道的热力状况进行分析, 通过保形映射、分离变量及区域衔接得到温度场的级数解;以分离变量所构造的特殊函数作为基本解, 根据边界条件直接离散赋值来定解。 1.2 温度场解析的解析法 图1a为物理平面上保温管道示意图, 图1b为保形映射后像平面上管道的非正交边界的离散示意图, 图1c为像平面上保温管道及其邻域和大地表面温度场区域的划分示意图, 图1d为像平面上保温管道及其邻域和大地表面上离散温度单元上热流的示意图。 在边界Ci(金属热力管道外表面) 上, 有 Τ=Τi(1) 在边界Co(大地表面) 上, 有 α(Τ-Τf)=-λ2?Τ?y(2) 式中T为物理平面上任意点的温度, K;Ti为金属管道的外表面温度, K;α为大地表面的表面传热系数, W/ (m2·K) ;Tf为大地表面附近的空气温度, K;λ2为土壤的导热系数, W/ (m·K) ;y为物理平面上任意点的纵坐标, m。 采用保形映射的方法, 取ζ(Ζ)=Ζ-(√Η2-r2c)iΖ+(√Η2-r2c)i(i为虚数单位, i2=-1) , 变换后, 把物理平面Z上的曲线Ci,Cc(保温层外表面) 和Co相应地映射到像平面ζ上的曲线Si,Sc和So, 其中Sc和So为同心圆, 且So为单位圆 (其半径ro=1) , 如图1c所示。 将计算域划分为Ⅰ区和Ⅱ区, 用TⅠ和TⅡ分别表示Ⅰ区和Ⅱ区的温度, 对于在极坐标系下任意点 (极径r, 极角θ) 可以得到: 区域Ⅰ?2ΤⅠ?r2+1r?ΤⅠ?r+1r2?ΤⅠ?θ2=0(3) 区域Ⅱ?2ΤⅡ?r2+1r?ΤⅡ?r+1r2?ΤⅡ?θ2=0(4) 在边界Si上, 有 ΤⅠ=Τi 在边界So上, 有 α(ΤⅡ-Τ0)=-λ2?ΤⅡ?r2√Η2-r2cx2+Η2-r2c 式中T0为环境温度, K;H为管道埋深, m;x为坐标, m。 在区域Ⅰ和区域Ⅱ的交界处, 有 ΤⅠ=ΤⅡ 在r=r′c(r′c为Sc的半径) 处, 有 -λ1?ΤⅠ?r=-λ2?ΤⅡ?r 式中λ1为保温材料的导热系数。 在无穷远点即平面ζ上的A(r=ro,θ=0) 点, 有 ΤⅠ=Τ0 对Laplace方程式采用分离变量法求解, 可以得到两个区域 (n=Ⅰ, Ⅱ) 内的一般解, 为 Τn(r,θ)=a(n)0+b(n)0lnr+∞∑m=1(a(n)mrmcosmθ+b(n)mr-mcosmθ+c(n)mrmsinmθ+d(n)mr-msinmθ)(5) 由极点引出2N+1条射线, 相邻两条射线间的夹角为2π/ (2N+1) 。这样就把圆周划分为2N+1个单元, 可将式 (5) 截断、离散成式 (6) (j=1, 2, …, 2N+1) a(Ⅰ)0+b(Ⅰ)0lnrj+Ν∑m=1(a(Ⅰ)mrmjcosmθ+b(Ⅰ)mr-mjcosmθ+c(Ⅰ)mrmjsinmθ+d(Ⅰ)mr-mjsinmθ)=1(6) 式中rj为像平面极坐标系下圆周上离散点j的极径;a(Ⅰ)0,b(Ⅰ)0,a(Ⅰ)m,b(Ⅰ)m,cm(Ⅰ),dm(Ⅰ)(m=1, 2, …,N) 均为待定的常实数;上角 (Ⅰ) 表示区域Ⅰ。 最后将像平面ζ上的解反演到物理平面Z上, 便得出直埋管道保温层内部 (区域Ⅰ) 以及保温层以外 (区域Ⅱ) 的温度场的解析形式的级数解。用于求解温度场一般解的系数矩阵的代数方程组由4个相互耦联的方程组成, 可转化为含有8N+4个未知数的矩阵方程, 使用计算机编写