m第4章 表观热惯量遥感信息模型.docVIP

第4章 表观热惯量遥感信息模型 4．1 热惯量的概念 4．2 表观热惯量(ATI)遥感信息模型 4．3 ATI遥感图像的计算 4．4 小结 物体的热惯量在热力学中是一个不变的物理量。因此，研究地物的特性可从不变的热惯量入手，例如不同的岩石、土壤、水体、植被等等。当热惯量有差别时，则可用热惯量作为识别地物的条件之一。研究热惯量一直是遥感中的难题之一，世界上只有少数几个国家的科学工作者推演了热惯量的遥感信息模型[1]，且至今尚没有发表过推演的过程。推演的假设有几种条件，本章也只能从一种假设出发进行推演。目的是使读者了解推导过程，以培养严谨的思维方法。 4．1 热惯量的概念 (一)热传导方程 物体在不同空间点有不同的温度，且温度还随时间变化，因此可以用温度场T(x,y,z,t)来描述温度的时空分布。温度场中存在热量流动，热量由较热部分向较冷部分流动，且温度差越大，热量流动得也越快。在场中任取一曲面单元dS，且确定法线方向n。根据热传导理论，在时间dt内沿法线n流过dS的热量dQ与dt，dS以及温度沿曲面法线的方向导数成正比，以k(0)为比例系数(称为热传导系数，单位是J/(m·s·K)，即焦/(米·秒·开))。于是，有 这样，当；为负时，即沿法线”的方向T变小时，热量dQ则为正。若以q表示热流向量，由热传导理论有q = -k grad T。于是，即有 dQ = qn dS dt 其中qn为在法线n上的投影，即 在物体内部任取一由闭曲面S所围成的区域V，可以看到：在时间dt内，从区域V流过曲面S向外流出的热量为 这里认为grad T(＝)0，即沿外法线方向的温度是上升的。应用奥氏公式[2]，并将上式变号，可得时间dt内向区域V内部流入的热量为 (4—1) 式中grad表示梯度，div表示散度。 另一方面，设以ρ表示物体的密度(单位是kg/m3)，γH表示物体的比热容(单位是J/(kg·K))，在时间山内，体积元素dV的温度增加 所以在这段时间内，这个体积元素流入热量为 从而在时间dt以内，全区域就需要吸收热量 (4—2) 根据能量平衡原理：流进和流出区域的热量相等，故(4—1)和(4—2)所表达的值相等，即 因为这个等式对物体内任何区域皆成立，从而有如下微分方程 (4—3) 这就是热传导方程。 对于均匀物体来说，k，γH，ρ均为常数，如果展开(4—3)式，则有 所以，可得 (4—4) 式中 (4—5) DH称为热扩散系数，也称为温度传导率(单位m2/s)。(4—4)式是均匀物体的热传导方程。 (二)热惯量 若对方程(4—3)在一维情况下给予一个简单的温度周期函数为边值条件，求解方程，可以得到一个简单的公式(具体推导见4．2节) (4—6) 式中Q为热量；T0为温度变化的振幅；ω为周日角频率；P是物体固有的常量，称为热惯量， 其表达式为 (4—7) 显然，热惯量P的单位是J／(m2·s1/2·K)。 从(4—6)式中可见，当地物吸收或释放的热量相等时，地物温度变化的振幅与地物的热惯量成反比，即：P值大的地物，日温差小；P值小的地物，日温差大。在力学牛顿第二定律F=ma中，当作用力相等时，加速度与质量成反比，人们把物体不变的质量看成是运动的惯性。P与惯性有类同之处，因此称为热惯量。热惯量的实质是地物阻止其温度变化幅度的一种特性。从微观的角度看，热惯量是阻止物体内部分子运动速度变化的阻力。在地物的周民温度变化中，热惯量起着决定性的作用。 热惯量是地物的体内特性，不是表面特性，因此不能用一般的遥感方法直接测量。热惯量P中的热传导系数k、密度ρ、比热容γH都难以直接与各不同波段的遥感信息相对应。因此，必须建立热惯量P的遥感信息模型。 从20世纪40年代以来，有许多科学工作者研究热惯量，直到1976年，Kahle A等人作出了一张加里福尼亚Pisgah火山口熔岩湖的热惯量图[3]。它引起了各国遥感工作者的重视。之后，我国也有人开始研究这一问题；1986年，我们利用航空遥感数据制作出表观热惯量图，建立了表观热惯量遥感信息模型[4]。1990—1992年，又制作出真实热惯量图；建立了真实热惯量遥感信息模型，并在此基础上完成了土壤含水量的遥感信息模型[5]。本章特讨论表观热惯量遥感信息模型，下一章中将讨论真实热惯量和土壤含水量的遥感信息模