泉水的合理开发利用论文.doc

泉水的合理开发与利用 摘要 本文首先对问题进行合理的分析假设，为建模建立了良好的环境；其次，通过更深入的分析问题，建立适当的模型，然后在模型中分离出题目中提出的几个问题，并采用分离变量、线性拟合等方法逐个求解，进行相关预测；最后，对建模求解的结论进行简要的综合论述，并提出几点对合理开发与利用泉水的建议。 正文 问题的分析与假设 泉是地下水的天然集中地表出露，含水层或含水通道与地面相交处产生地下水涌出地表的现象。多分布于山谷和山麓，是地下水的一种重要排泄方式。它是在一定的地形、地质和的结合下产生的。适宜的地形、地质条件下，潜水和承压水集中排出地面成泉。泉水流量主要与泉水补给区的面积和降水量的大小有关补给大、降水越多，则泉水流量越大。水流在重力作用下呈下降运动，季节性变化:泉流量，单位 :降雨量，单位 :地下水开采量，单位 :泉水补给区面积，单位 ：外来水的总补给量，单位 ：排除降雨的作用时外来水对泉水的总补给量，单位 ：年份，单位 年 ：常数系数，在后面的式子中详细给出，=1,2,3… 1.2分析 泉及其相关因素组成的系统可近似等效为如图1所示的连通器。 如图1中，B液面可看作泉的水面，A液面比B液面高出h的部分来自降雨和地下水的对泉的补给。由于A、B存在液面差，B液面会有水溢出，其水流量等效为泉流量。 实际中，某一段时间泉水存储总量远大于外来水的补给量，外来水的总补给量正比于这里的h，B处的水流量可近似看作与h成正比，因此有如下关系式： 泉流量 假设没有降雨和地下水开采时，外来总补给量Z是一个与地下水状况随时间变化有关的随机函数，此时Z是关于时间的函数，设为 为带有很大随机性的的函数，需通过对已知数据的拟合以确定其形式。 当加入降雨量时，外来水的总补给量为 其中S为泉水补给区面积，为对应常数。 当地下水被开采时，外来水的总补给量可近似看作随开采量线性下降，则有 所以泉流量与降雨量和地下水开采量的关系为 可简化为 1.3建立模型并求解 如图2将该泉看做一个黑箱，这个黑箱使“输入”降雨量J、地下水开采量C、随机函数F（t）和“输出”泉流量L满足关系式： 利用Matlab7.0软件通过已知数据分别绘制出L-J、L-C、J-t的坐标图，并进行统计分析，找出、 、的确切表示式。 既表征地下水状况的变化，也直接表达了泉流量与随时间变化的地下水状况的关系，地下水状况的变化受很多因素影响，是自然性的随机变化，根据概率知识：[3]如果一个量是由大量的相互独立的随机因素影响所造成的，而每一个因素在总影响中所起的作用不很大，则这个量通常近似服从正态分布。故取值近似服从正态分布。 ① 对L-J(1983年以前)坐标图分析 如图3为利用1983年以前的数据绘制的L-J坐标曲线图，从图中可看出曲线的总体趋势为：未开采地下水以前，随降雨量的增加，泉流量也曾上升趋势。利用Matlab7.0软件可拟合出该曲线对应的线性曲线： 又1983年以前，地下水开采量为零，所以泉流量为 故。 再利用式子可计算出的一组对应数据，从中抽样十组数据得： 4.354 3.445 4.590 5.335 3.158 2.547 2.818 2.343 2.960 ② 对L-C（1984年以后）坐标图分析 如图4为利用1984年以后的数据绘制的L-C坐标曲线图，从图中可看出曲线的总体趋势为：随地下水开采量的增加，泉流量呈下降的趋势，降雨量对泉流量的影响不再明显。此时泉流量符合公式： 利用Matlab7.0可绘制出以为纵轴，为横 轴的坐标内的点图，如图5 所示；其中直线为通过点拟合的一条近似线性直线： 结合的取值服从正态分布和公式： 可确定。 再用导出式 可求出对应于不同的、值时的取值，下面对抽样十组数据： 4.008 3.669 2.818 2.570 1.541 2.720 2.257 3.078 2.377 2.007 分析前后对抽样的二十组数据，其最大值为5.335，最小值为1.514，并可算出其样本均值为3.1055，可近似看作的期望值，样本方差为0.934，可近似看作的方差值，则服从正态分布，的取值在这一范围内的概率为91.49%，同时从的取值看，后一组数据总体上比前一组数据小，由此可见的取值同时具有两个特性：⑴在某一段时间内的取值近似服从正态分布，⑵的取值整体上呈减小的趋势。 综上可得泉流量的关系式： ③ 对J-t坐标图分析 如图6为降雨量随年份变化的坐标曲线图，从图中分析可知：1979年以前降雨量呈周期性的起伏，分别在1959年、1965年、1