中深层地热热泵单孔同轴套管式换热器传热数值分析方法.docxVIP

附录B单孔同轴套管式换热器传热数值分析方法 B.1假设： 1） 将中深层地热换热器周围的岩土层看作为几个均匀介质的水平地层，并忽略可能的地下 水流动的影响，将岩土中传热视为单纯的导热问题。 2） 数值模拟区域选择在距离埋管中心足够远的地方作为径向边界，该处的温度分布被认为 不受中深层地热换热器的影响。 3）忽略空气以及岩土层表面的全年温度波动的影响，地面的空气温度被认为是不变的。 ）通过整个介质的大地热流被认为是恒定的。 ）管内循环液的热对流是钻孔内轴向传热的主要途径，忽略钻孔内的轴向导热 B.2控制方程 ????????每个岩土层的导热方程可写为： ???? ???? 1 ???? 1 ?? ????????? ?????????? ??上式中， ?? 上式中，？?=云?为岩土层的热扩散率，单位 m2/s; t为岩土层温度，单位 C； t为时间，单 位s; r和z为径向长度与竖直方向长度，单位 m。 在径向采用不同的步长，采用坐标变换。设定一个新的坐标: ????= ?? ??= ????( ?? 在新坐标下的导热方程写作: 在6坐标系当中均匀的步长1 ???? 1 ???? ?????????? ?? 在6坐标系当中均匀的步长 ???+1_△ b相当于在r ???+1_ ?? =??????? = ?? ?? 上式中B为径向步长放大倍数，即松弛因子。 套管内流体温度的控制方程应考虑流体的流动方向。流动方向不同时，方程中对流项 前的正负号也随之改变。 以下方程中设定 z坐标的方向向下流动方式为内进外出（流体从内管流入套管，从两 管的夹层中向上流出）时， 内管流体的能量方程为： ???? ???1- ???2 ???? ?? = - ?? ???? ?? …???? 上式中tf1与tf2分别为内管流体温度和外管流体温度，单位 C ； R2为内管管壁热阻及内管内 侧对流热阻与内管外侧对流热阻之和，单位 K/W ； C2为单位长度的内管中流体以及内管管 壁的热容量，单位 kJ/（Km）; C=Mc指的是循环液的热容流量 kJ/（sK）。 流动方式为内进外出时外管（两管的夹层）中流体的能量方程为 単？??^+ ??^1+ ?彈 ???? ?? ?? …???? R1为外管壁以及钻孔内回填料热阻，单位 K/W； C1单位长度的外管中流体以及外管管壁的 热容量与填充层热容量，单位 kJ/（Km）; tb为填充层外壁温度，单位 C。 流动方式为外进内出时外管流体的能量方程为 ： ??+1 ??+1 其中: ???? ???2- ???1 ?? = + ???? ?? 流动方式为外进内出时内管流体的能量方程为 ???? ???1_ ?? ???- ???1 ?? ???? ■ ???? ???2 ???? ??— + ?????? …???? n 2 2 Cl = 4(dii - d2。) n 2C2= 4d2i pc 1 ???逊)2???% …(???? ?? = + —1— ?????°) 2 ?????2 2????2 ??? kg、kp1和kp2代表回填料、外管和内管的导热系数，单位 W/(m 2 K); rb代表钻孔半径，单位m;di和do管道的内径和外径，单位m ；1和h2代表套管外管和内管的对流换热系数， 2 - B.3初始条件和边界条件 岩土中的初始温度分布在径向是均匀的，在深度方向有地温梯度。因为认为大地热流是均匀的，由于不同地质层导热系数会有不同， 所以在不同的水平地层中有不同的地温梯度。 在任一深度处的地层中的初始地温可以表达为： 1 1 ?????1 1 1 ??? + ????p) 2????…??° 1 + ???? ?1 W/(m 2 K); W/(m 2 K)。 ??_1 ??? ??? ??? ?(??) = ???+ + 刀 (??7 ???_1) + (??_ ???_1 ) ??? ??=1 ??? ??? T= 0, Hm_lW z Hm , rbWr 锋bnd 其中Hj是第j层地层底部的坐标；qg是大地热流，单位 w/m2； ta是地表温度，单位°C ； ha 是空气与地表的对流换热系数，单位 W/(m 2 K); km 对于岩土中的导热问题，圆柱坐标的径向边界设在 在这个径向边界上设定第一类边界条件， 化，即卩 ??? ?(????????)= ???+ — + ??? W/(m 2 K)。 匕壤的导热系数，单位 r= rbnd处，距离套管的中心足够远。 即认为边界上的温度在研究所涉及的时间内没有变 ??_1 ??? ??? 刀 鬼(？加???_1) + 看(？?_ ?加_1) ??=1 ?… ?… ??= W ?* ???,?? 0 在地表的边界上设定第三类边界条件，即假