凌河扇地下垫面对地表径流形成与运动的影响分析.docxVIP

凌河扇地下垫面对地表径流形成与运动的影响分析 1 地表单因素模型的建立 在传统的地下水-地表耦合模拟过程中，地表径流被视为光滑表面的平坦运动。地下水和地表水之间的水流量差只考虑两者的水差以及介质的渗透性质。降水和地表径流是一个复杂的过程。降水通过植物表面植被的渗透和蒸发反应而进入气体带，然后产生地表径流。地表径流的形成与运动过程中降水和空气带的性质有关，同时也受到各种因素的影响。 由于实际的地表面是一个极其不规则的曲面,不能通过数学模型对其进行描述,而传统的处理方法是假设地表是平缓光滑的,认为地表的洼地是由于DEM数据的采集和插值过程中产生的“错误”,在建立模型之前首先进行DEM的填洼处理,直至近似平滑的曲面。但真实的下垫面条件是坑洼不平的表面,地表的洼地往往可以截留大量的降水,减缓流域的产汇流速度,增加地表径流的滞留时间及对地下水的补给,对流域内地下水与地表水的形成及分布产生影响,传统的处理方法显然不能够真实反映和模拟降水产汇流及对地下水的补给过程。为了准确的模拟流域内水文循环过程,为当地的水资源评价及防洪工程提供可靠的依据,应当将洼地作为一个重要的影响因素进行考虑。 本文在地表径流运动方程中加入地表孔隙度项(Surface Porosity),推导出地表孔隙度和径流深度与下垫面的关系式,并以辽宁省凌海市大、小凌河扇形地为例,分析并讨论下垫面的洼地对当地水文循环的影响。 2 地表径流剖面位置及边界条件 当研究区比较平缓,地表水流动缓慢时,地表径流的水平运动尺度远大于垂直运动,忽略垂向流速。为了考虑下垫面对地表径流的影响,采用如下修正形式的二维圣维南(Saint Venant)扩散波方程及定解条件来描述地表水的运动。 式中:Φ0为地表孔隙度,无量纲;ho为水面高程(m)(ho=do+z);do为地表径流深度(m);z为地表高程(m);Kox,Koy为x,y方向上的地表传导系数(Surface Conductivity)(ms-1);Qo为地表水的源汇项(ms-1);Do为地表水范围。 h(x,y,t)为给定水深的第一类边界条件;qo1(x,y,t)为给定通量的第二类边界条件;Kn为边界法向的地表传导系数张量(ms-1);qo2(x,y,t)为临界深度边界条件(Critical Depth Boundary Condition),用于强制边界上的水深为临界深度;qo3(x,y,t)为零深梯度边界条件(Zero-depth Gradient Boundary Condition),用于强制边界上的水深梯度为地表坡度;So为地表坡度,无量纲;h(x,y,0)为初始条件。 其中Kox,Koy可由下式求出: 式中:nx,ny分别为x和y方向上的曼宁粗糙系数(m-1/3s)。 3 影响研究与实例分析 3.1 单次填洼处理下的单次消减储水量 在基于上述数学模型对地表径流进行模拟时,由于实际的地表面是粗糙不平的,地表的洼地、沟渠、植被及建筑物都会对地表径流的形成与运动产生影响。下面就在理论上推导下垫面条件与地表径流形成与运动的关系。 下垫面对地表径流的影响可以通过洼地储水量和消减储水量两个概念来进行评价: 洼地储水量(Depression Storage/Rill Storage):由于地表的坑洼、沟槽等洼地引起地表水体储存量的滞留,即地表降水产流之前必须先填满的储存量,以洼地储水高度Hds(m)来表示。 消减储水量(Storage Exclusion/Obstruction Storage):由于地表植被、地形、建筑物等障碍物引起储存量的消减,同样还会引起摩擦阻力增加与能量的耗散,以消减储存高度Hos(m)来表示。 当地表水高程ho在LS~LS+Hds+Hos(LS为地表高程)间变化时,由地表水所覆盖的地表面积Ω在0~Do间变化,地表孔隙度在0~1之间变化,以体积高度VH(m)(Volumetric Height)来代替地表径流深度do。VH~do关系如图1所示。 当地表径流深度do≥Hds+Hos时,VH=do,Φ0=1;当do∈(0,Hds+Hos)时,具有如下关系: 对上式进行Hermite插值求解,构建出VH~do函数: 即有如下关系: 因此可以将式(4)和式(5)表示公式(1)中的ho,式(6)表示,得出考虑下垫面影响的地表径流运动公式。 洼地储水高度可以通过“洼地填充”的方法获取。利用Arcgis软件中的Spatial Analysis模块和Hydrology工具箱对研究区的原始DEM数据进行填洼处理,比较填洼处理前后DEM的体积,将填洼处理前后的DEM相减,即得到洼地储水高度的分布情况;将体积变化量除以研究区的面积即得到平均洼地储水高度;消减储水高度可以通过对地表的地形特征进行统计分析获取,本文不做深入探讨。 3.2 地表洼