草地水量问题的数学模型2.doc

草地水量问题的数学模型 刘海景1 林锡波2 王燕芹3 1．韶关学院2004级数学系应用数学班 2．韶关学院2006机电系06自动化 3．韶关学院2006级数学系数学教育班 摘 要 本文针对网球场草地的降水后情况,首先分析草地干燥的过程,定义了水高度(水面与草地地面的高度差,也称水头差)的概念,巧妙地把问题转换成对水高度的分析的问题,建立了描述草地干燥过程的微分方程模型;其次对微分方程中的未知量进行理论分析，并且在假设描叙下雨速度的数学函数遵循某一曲线的情况下分析了模型,分析了雨过后草地水量的减少过程.得到水位高最高值是在下雨速度由高到低这一过程取得,时间为下雨速度等于蒸发、蒸腾和渗透速度之和这一时刻.并根据所建立的模型预测雨停后经过时间时草地变干燥. 关键词:水头差； 渗透系数；水压梯度；蒸发率 1 问题的提出 草地网球比赛常因下雨而被迫中断,由于防水层不一定有效,往往需要经过一段时间使草地的最上层充分干后,才能继续进行比赛,雨停之后,部分雨水直接渗入地下,部分蒸发到空气中去.虽然有一些机械装置可用来加速干燥过程,但为避免损伤草皮,最好让草地自然地变干.试建立一个描述这种干燥过程的数学模型,并说明雨停止以后草地里的水量减少过程. 2 模型的假设 2.1 不考虑除自然变干外的其他因素(如有排水道,一些加速干燥的机器等); 2.2 所研究的对象除了下雨带来的降水外和外界水体系没有联系; 2.3 草地的草分布是均匀的并且每一棵小草的蒸腾的速率一样的; 2.4 在所研究的面积内草地自由水面的蒸发速率处处相等. 3 符号的约定 3.1 f(t) t时刻下雨的速度，(毫米/秒); 3.2 h(t) t时刻草地上面的水高度(也称水头差)， (毫米); 3.3 v(t) t时刻渗入地面的速度(随着水高度的变化而变化) ，(厘米/秒); 3.4 (t) t时刻草地的蒸腾速度(毫米/秒); 3.5 (t) t时刻草地蒸发的速度(毫米/秒); K 渗透系数; 3.7 H 土柱上水头差（厘米）即静水压力； 3.8 I 水压梯度; 3.9 下雨时间(秒); 3.10 L 发生水分渗透作用的土层的厚度（厘米）即渗透路程。 4 问题的分析 此题的目的是要我们分析从下雨开始时到草地最上层干燥后这一过程草地水量变化的问题.草地水量变化的问题可以对草地水的质量变化的分析,也可以是对水体积大小的分析,这里我们采用草地水高度变化来分析草地水量变化的问题,当水位最高时表示此时草地最为潮湿,当水位高度为零时就说草地土壤表面干燥了. 首先，采用微分方程模型对草地水量进行分析，得到相应的微分方程。经过分析当水高度为零时候说明草地土壤表面干燥了。 其次，对微分方程一些未知量进行讨论。通过参照相关书籍和网上相关论文，给出微分方程未知量相应的推导理论依据。 最后，对模型的讨论和推广。 5模型的建立和求解 5.1 下面采用微分方程模型对此问题进行分析，如图： 下 雨 蒸腾 水面蒸发 水高度（水头差） 地面渗透作用 通过分析显然有以下的结论和式子: (1) 式中：表示时刻草地上的水高度；表示时刻草地蒸腾的速度；表示时刻草地地面渗透的速度；表示时刻水面蒸发的速度。 5.2式(1)中各符号的讨论: 5.1.1 的分析。显然是一个变化的量,每一场雨它的下雨速度的数学函数都不同,但是雨速的变化都有两个过程:小大,大小.这里我们把下雨速度达到最大值的时间定在，并且下雨过程中的雨速函数的轨迹遵循某一抛物线.经分析这抛物线有以下几点性质:开口向下、经过坐标原点.所以我们可以写出如下的函数方程. (2) 其中是下雨时间,为下雨的最大速度.经过分析其最大速度可以是反映这一类型雨的雨速的最大值 (比如是中雨的f(/2)，就用中雨和大雨的分界雨速来表示). 5.2.2 渗透系数的测量理论分析。在饱和水分土壤中，渗透性按照达西公式计算如下： 式中：为渗透速度，每秒钟通过1平方厘米土壤断面的水的流量，以立方厘米表示；为水压梯度，即渗透层中单位距离内的水压降；为渗透系数，在单位水压梯度（I=1）下，单位时间内通过单位截面积的流量（毫升/分或小时）；为土柱上水头差（厘米）即静水压力；为发生水分渗透作用的土层的厚度（厘米）即渗透路程。 在时间t内渗透过一定截面积A（平方厘米）的水量Q，可以用下列的方程式来表示： 因此渗透系数 （毫米/厘米2/分或小时） 土壤渗透性的测定有室外法（渗透简法）及室内法（环刀法）。 5.2.3 v(t)的分析。