基于DEM的太阳辐射模拟 .doc

基于数字高程模型的太阳辐射模拟 汤国安1，2，杨 昕2，张 茜2 ，邓凤东3 （1.西北大学 教育部大陆动力学重点开放实验室，2.西北大学 城市与资源学系，陕西 西安 710069； 3.陕西省农业遥感信息中心，陕西 西安， 710015） 摘 要：针对地形变异对太阳辐射在地面分布的影响，为最大限度地利用太阳能，提出了利用数字高程模型为基本信息源，提取的数字坡度模型和数字坡向模型，并采用多层面复合分析的方法，通过地形要素与其他气候要素的复合分析，建立了区域地面的太阳辐射模型。其中，每一个栅格单元瞬时太阳辐射能量的计算是其中最基础的工作，通过对瞬时太阳辐射量的时段累计，可获得日、月、年的太阳辐射总量。该结果可望进一步转换成对农作物生长有重要影响的积温、无霜期长度等重要的小气候指标，为精细农业的实施提供重要的前提条件。 关 键 词：数字高程模型； 太阳辐射； 模拟 中图分类号：P28 文献标识码：A 文章编号：1000-274X(2003)000-08 太阳辐射量大小是影响区域气候最重要的因子。一个地区由于所处一定的地理纬度，已经决定了其太阳辐射的天文总量。然而，由于地面起伏变化造成局部地面接收阳光的状况存在很大的差异，太阳辐射量在地面还存在一个重新分配的过程，不同的坡度、坡向、坡形、坡位，所受太阳辐射量的不同，造成局部的小气候特征的差异，同时，也直接影响到诸如土壤水分、地面无霜期以及作物生长适宜性程度等多项重要的农业生产指标。因此，对于微观地块尺度，如何利用数字高程模型（DEM）数据与GIS空间分析方法，实现对地面所获得太阳辐射量的数值模拟，是需要加紧解决的问题。国外同类研究首推Dozer（1990）的工作[8，9]，他提出了利用数字高程模型模拟太阳辐射的方法，Gilabert (1990), Civco(1993) 也在这方面探讨了任意地形条件下的辐照入射度模拟[7，10,11]。我国对于任意地形条件下太阳辐射的开创性研究是由傅抱璞[1，2]做出的，他给出的计算坡地临界时角的公式和日照时断的判定方法简化了山区太阳辐射的计算。李占清、翁笃鸣、朱志辉等发展了这一方法[5]，并引入了实际山地中的散射辐射和总辐射的模拟[3，4]。1999年，李新等利用DEM进行了太阳辐射的模拟实验[6，7]，但其所用的DEM分辨率较低（2.4km2×2.4km2），仍不能明显地体现太阳辐射在局地的再分配与小气候的差异性。本研究利用国家地形数据库高精度的DEM（1：10 000比例尺，5 m分辨率），选择黄土丘陵区为实验样区，研究在相同太阳辐射天文值的条件下，由于地形起伏造成辐射能量在空间的分布规律。即：通过模拟某一瞬时太阳辐射强度，获得太阳辐照度分布特征；通过对日照时段的积分，获得一个地区太阳辐射日总量分布，进而获得任一地块尺度月、年的太阳辐射数据。该研究结果对于区域小气候研究、农业区划、太阳能有效利用乃至精细农业的实施，有着重要的参考意义。 1 太阳直接辐射的计算 1.1 基础参数 1)地球轨道修正因子 式中称日角，即=。 其中：t = N - 1，N为积日，即日期在年内的顺序号。 2）赤 纬 式中的含义同式（1）。 3）太阳时角 =) ， 其中：时和分表示真太阳时,太阳时角经过年度订正、经度订正、时刻订正，以太阳正午时刻为0点的，顺时针方向（下午）为正，反之为负。 4）太阳高度角 ， 太阳方位角 ， 其中为地理纬度。 5）平地日出日没时角的计算 其中：为日出时角；为日落时角。 6）天文辐射瞬时值模型 , 其中：为坡度；为坡向；为太阳常数1 367 W/m2 。 7）太阳天文日辐射量的模型 式中：T为日长；n为方位划分数，取经验值n = 36；为微分时段内的日出和日没时角。 1.2 DEM及在太阳辐射模拟中的作用 数字高程模型是地面连续高程变化的数字表示, 又以栅格数字矩阵为多。图2(a)为栅格化地面的三维模型，图2(b)显示所建立的栅格单元，具有一定的地面坡度与坡向，对于任一太阳入射光，必然具有相应的瞬时所接受的辐射量。因此，DEM的建立为栅格单元太阳辐射的计算创造了重要的条件。其重要性主要表现在：能够提取模拟所用到的坡度、坡向值；能够利用光线追踪算法计算地形遮蔽信息，以及模拟结果的区域可视化表达等。 值得注意的是，数字高程模型的地形描述精度在很大程度上受其信息源精度与栅格分辨率的影响，了解研究的空间尺度条件，