针对线性加热装置在热传导过程中测试点温度随加热时间的变化问题，研究了加热装置中心点处辐射值E与功率P之间的关系式，基于周边传导介质的性质，推导了考虑时间延迟参数的一维热传导模型，并通过加热棒为装置进行了热传导实验，验证了理论模型的适宜性。研究结果表明，受热传导影响而测试点温差逐渐提高时，测试点温度变化斜率逐渐减小并趋于稳定。通过对理论计算的温度值与实测温度值进行了对比表明，实验初期相对误差最大达到20%，但随着加热时间和热源功率的增加，两者的温差逐渐减小并趋于零，说明本文推导的理论计算公式符合实际。 In order to study the relationship between the changes of linear heating device testing point temperature and heating time in the heat transfer process, we study the relevance between radiation values E and the power P at the center of the heating device. Based on the nature of the surrounding conductive medium, we derived one-dimensional heat conduction model which considered the time delay parameters. We verify the suitability of the theoretical model through laboratory experiments by heating bar arrangement. The results show that test point temperature slope is gradually reduced and stabilized when the test point temperature difference gradually rises in the effect of heat conduction. Through the measured temperature values compared with theoretical calculations, it showed that the maximum relative error reached 20% in early experiments, but the temperature difference decreased and tended to zero with the increase of heating time and source of power, showed that theoretical formula derived realistic in this paper.
刘子文*，朴春德#，谢亮，李斌，杨大帅，潘东玥
中国矿业大学资源与地球科学学院，江苏 徐州
收稿日期：2016年3月28日；录用日期：2016年4月18日；发布日期：2016年4月21日
针对线性加热装置在热传导过程中测试点温度随加热时间的变化问题，研究了加热装置中心点处辐射值E与功率P之间的关系式，基于周边传导介质的性质，推导了考虑时间延迟参数的一维热传导模型，并通过加热棒为装置进行了热传导实验，验证了理论模型的适宜性。研究结果表明，受热传导影响而测试点温差逐渐提高时，测试点温度变化斜率逐渐减小并趋于稳定。通过对理论计算的温度值与实测温度值进行了对比表明，实验初期相对误差最大达到20%，但随着加热时间和热源功率的增加，两者的温差逐渐减小并趋于零，说明本文推导的理论计算公式符合实际。
关键词 :一维热传导，温度分布，温度梯度
通过温度场研究渗流场是在国内外对于温度示踪法研究的基础上发展起来的，为研究渗流场提供了新的研究理论和方法 [
综上所述，虽然有关热传导特征和温度识别方法方面取得了不少研究成果，但有关一维热传导条件下与热源不同距离的测试点温度变化理论模型研究方面还需进一步探索。本文基于傅里叶导热定律，推导了考虑时间延迟参数的一维热传导模型，并以加热棒作为热源，通过加热棒在恒温水槽中的温度变化特性实验，得到不同功率下的温度分布特征，验证了理论计算公式的合理性。论文的研究成果对于一维加热环境下岩土体介质等材料的温度传导特性计算及研究方面具有良好的借鉴意义。
对于一维线性加热装置在热传导过程中测试点温度随加热时间的变化问题，假设与热源一定距离的目标质点在单位时间、单位面积内接受的热量值为E (以下简称为辐射值E)，该质点在单位时间、单位面积内散失的能量值为H0，则热传导条件下该质点的辐射值与温度t、某时间段τ的关系为 [
式(1)中，
式(1)中热量从热源散发到达目标质点时，其间一部分热量由目标质点吸收，另一部分则由目标质点与热源之间的填充物所吸收，即目标质点接受的热量值只是热源散发热量的一部分。热量以热传导方式从热源至目标质点之间传播时，存在温度传导的时间差，即时间延迟τ0。因此，对于一维热传导条件下目标质点温度与热源功率及时间的变化关系，应在式(1)的基础上需要考虑目标质点接受的热量值与热源总散发热量的比值关系和时间延迟τ0。
对于长为L，电阻半径为r3，功率为P的加热棒，管壁的表面积为2·π·L·r3，因为辐射值为单位面积散发的能量值，功率为散发的总能量值，则加热棒管壁任一点处其辐射值与功率之间的关系为
式(2)中，n为目标质点接受的热量值与总热源散发热量的比值。当加热功率和周边材料介质的性质及环境温度保持不变时，假定目标质点吸收的能量与热源散发的能量的比值为恒定值。因此，可以把加热装置中心点温度响应特征等效为稳态情况下无限长圆筒壁表面温度的响应特征，根据稳态无限长圆筒壁的温度分布，求出目标质点吸收的能量与热源散发的能量的比值。
通过傅立叶导热定律，可以得到稳态情况下无限长圆筒壁的温度分布方程 [
根据傅立叶导热定律，对于无限长圆筒，距热源中心轴线r处质点温度为
对于某一质点的比热容为C、质量为M，若其温升为Δt，则该质点吸收的能量为C·M·Δt。对式(3)从r1至r2进行积分，得到热源传导过程中其它质点吸收总能量Q1，即
式(4)中C1、M1分别为热源与指定质点之间充填物体的比热容及质量，t0为质点的初始温度。
当质点距中心r2至r3之间时，质点所处位置所吸收的总能量为Q2，即
式(5)中C2、M2分别为圆柱体管壁的比热容及质量。
由式(4)和式(5)可以得到质点吸收的能量占热源散发能量的比值n，即
图1. 无限长圆柱体温度分布图
由式(6)可以得到不同功率下目标质点接受的热量值与热源散发热量的比值n。
当周边介质的温度保持不变时，距热源一定位置的质点其热量损耗值与温度之间有如下关系 [
式(7)中，当实验在水体环境下进行时h为热对流系数，若在岩土介质环境时h为导热系数，t为加热后质点某时刻温度。
在岩土介质环境下，将热量损耗值式(7)及其目标质点辐射值与功率之间的关系式(2)代入式(1)，得到在岩土介质环境下，距热源一定距离的目标质点，其温度的响应特征。即
时间延迟量值τ0，其大小与充填的物体性质有关，可通过实验进行测定。考虑时间延迟影响，将上式(8)时间τ改为τ-τ0，即
实验时将内有电阻加热丝的加热棒放置于恒温水浴箱，并在玻璃管表面中点处黏贴温度传感器，测量石英壁外壁中点处的温度变化。有关恒温水浴箱尺寸为40 cm ×35 cm × 15 cm，如图2所示。加热棒的结构为空心石英玻璃管的中央放置直径为1 cm、长度18 cm的电阻加热丝作为热源，加热丝与玻璃管内部由氧化镁粉末充填。有关空心石英玻璃管的尺寸分别为长度21.3 cm，直径2.0 cm，壁厚0.2 cm。有关加热棒结构示意图，如图3所示。
实验前将功率为200 W的加热棒放置于充满水的水浴箱并进行加热实验，水浴箱中水的初始温度为
图2. 实验装置简图
图3. 加热棒结构示意图
14℃。试验时将电阻丝加热到石英管壁出现温度变化的时间段作为时间延迟τ0，并按照一定温度梯度记录加热时间，当玻璃管表面温度不发生变化时终止实验。按以上步骤，更换水浴箱中的水，放置功率为300 W加热棒进行下一组实验。
通过加热棒外表面中点处安装的温度传感器，得到加热功率分别为200 W和300 W时温度随时间的变化关系。此次实验，式(9)中t为加热棒中点的温度值；t0为其初始温度，取水温值14℃；n为加热棒中点接受的热量值与热源散发热量的比值，此值与热源功率及目标质点位置有关，可由式(6)求得不同功率下的的值；C为加热棒中点(材质为不透明石英玻璃)的比热容891.8 J/(Kg·K)；L为石英玻璃管长21.3 cm。
通过加热棒的参数代入式(9)。得中到加热棒外表面中心位置温度随时间的变化曲线，并与实测的温度值进行了对比，计算了相对误差。有关功率分别为200 W和300 W时加热棒外表面中心温度随时间变化曲线，如图4、图5所示。图中，时间轴的起始值为时间延迟τ0，温度轴的起始值为初始温度。
由图4可知，通过理论计算的温度变化曲线与温度实测值均随加热时间的增加其曲线斜率逐渐减小并趋于稳定。其主要原因为，由于加热初期石英管壁与电阻丝之间的温差较大，石英管的温度上升较快。随着加热时间的增加，石英管壁与电阻丝之间的温差逐渐减小，同时受恒温箱中水温的影响，石英管表面的热量损耗增大，最终管壁接收到的热量与散失热量逐渐相等，石英管壁的温度趋于稳定。通过理论计算的温度值与实测值进行对比表明，加热初期两者相差较大，其相对误差达到12%，但当加热时间增加至50 s后，两者的相对误差将小于3%，说明此时两者的结果基本一致。
由图5可知，功率300 W时因热传导引起的石英管温度变化曲线与功率200 W时大致相同，即随着
图4. 200 W加热棒外表面中心温度随时间变化曲线
图5. 300 W加热棒外表面中心温度随时间变化曲线
加热时间的增加而温度逐渐增高，但其变化率却逐渐减小。通过理论计算得到的温度值与实验值进行对比表明，两者在实验初期的相对误差将达到20%，但随着加热时间的增加两者逐渐趋于一致，特别是当加热时间增加至40 s后，两者的相对误差小于6%。
通过上述分析，得到以下结论：
1) 通过考虑时间延迟参数的一维热传导模型，得到了加热棒表面温度随时间的变化曲线。结果表明，由于加热初期石英管壁与电阻丝之间的温差较大，测试点温度随时间的斜率近似为直线，但随着加热时间的增加，两者的温差减小，温度变化斜率逐渐减小并趋于稳定。
2) 通过不同功率下对理论计算的温度量值与温度实测值进行了对比，其结果两者在实验初期相对误差最大达到20%，但随着加热时间的增加，两者的温度逐渐趋于一致，说明本文推导的理论计算公式符合实际。
国家自然科学基金青年基金(51004103)，江苏省高等学校大学生创新创业训练计划 (201410290007Y)。
刘子文,朴春德,谢亮,李斌,杨大帅,潘东玥. 一维热传导条件下测点温度与加热时间的相关性研究 Correlation between the Measured Point Temperature and the Heating Time under the Conditions of One-Dimensional Heat Conduction[J]. 地球科学前沿, 2016, 06(02): 72-78. http://dx.doi.org/10.12677/AG.2016.62009