第十一讲 估计水箱水流量模型

第十一讲估计水箱水流量模型

一、问题的提出

随着社会和经济的不断发展，环境和资源问题日益突出，水便是其中的主要问题之一。1997年联合国水资源会议曾郑重向全世界发出警告：“水，不久将成为继石油危机之后的下一个社会危机”。我国是一个缺水的国家，人均水资源拥有量仅为2150m3/a（按13亿人计），不到世界人均水平的四分之一，排在世界第109位。特别是“三北”（东北、华北和西北）地区和经济发达的沿海地区，水的供需矛盾已十分突出。有关资料表明，我国每年因缺水而

影响工业产值已达2300多亿元。预计到本世纪末，全国年总需水量将达到700亿m3，而缺水量也将达到70亿m3，水资源短缺已成为制约我国经济和社会发展的重要因素。

某些地区的用水管理机构为了达到节约用水的目的，需估计公众的用水速度（单位是G/h）和每天总用水量的数据。现在许多地方没有测量流入或流出水箱流量的设备,而只能测量水箱中的水位（误差不超过5%）。当水箱水位低于某最低水位L时，水泵抽水，灌入水箱内直至水位达到最高水位H为止，但是也无法测量水泵的流量，因此在水泵启动时不易建立水箱中水位和水泵工作时用水量之间关系。水泵一天灌水1～2次，每次约2h。试估计在任意

法测量水泵的流量，因此在水泵启动时不易建立水箱中水位和水泵工作时用水量之间关系。水泵一天灌水1～2次，每次约2h。试估计在任意时刻(包括水泵灌水期间)流出水箱的流量，并估计一天的总用水量。

表1给出了某镇中某一天的真实用水数据，表中测量时间以秒为单位，水位以E为单位。例如3316s以后，水箱中的水深降至31.10E时，水泵自动启动把水输入水箱；而当水位回升至35.5E时，水泵停止工作。

本问题中使用的长度单位为E(=30.24cm）;容积单位为G(=3.785L(升))。水箱为圆柱体，其直径为57E.

二、问题的分析

1. 数据变换：水位与时间的关系转化为水箱中水的体积与时间的关系.

2. 探求水箱中水的流量与时间的函数关。

定水位对流速没有影响。类似地，还假设大气情况，温度变化等对水流速均无直接影响。

3．水泵的灌水速度为常数，不随时间变化也不是已灌水量的函数，因此假设水泵大约在水位27E时开始灌水，在水位35.5E时停止灌水。同时假设水泵不会损坏或不需要维护。

4．从水箱中流出的最大流速小于水泵的灌水速度。为了满足公众的用水需求不让水箱中的水流尽。

5．每天的用水量分布都是相似的。因为公众对水的消耗量是以全天的活动（诸如洗澡、做饭、洗衣服等）为基础的，所以每日用水类型是相似的。

6．水箱的流水速度可用光滑曲线来近似。每个用户的用水需求量与整个区的用水需求量相比微不足道，而且它与整个社区需求量的增减情况是不相似的。

四、分析与建模

引入如下记号：

V、V i —水的容积，时刻t i水的容积G）；

t i—时刻（h）；

p —水泵的灌水速度（G/h）;

T0 —初始数据的当天测量时间;

T —当天的时间（以24小时制）;

f(t) —流出水箱的流速是时间的函数（G/h）;

从此表的数据可知,在水泵的第二次抽水时间区间内,从水箱中流出的平均流速必定大于第一次抽水时的流速。因为第二次灌水时间长得多。但是，第二次水泵停止灌水的确切时刻是未知的，在水泵第二次灌水后所测得的水位与水泵实际停止运行的水位是不同的,因此,水泵停止运行的时刻一定在图2右边阴影部分两次测量时间之间.我们需要估计水泵开始和停止灌水的时间.

水泵第一次约在8.968h后开始运行的,因为这时水箱中水的体积约为514800G,水泵停止灌水时间在10.926h～10.954h之间的0.028长的区间内。

我们可以确定水泵第二次开始运行时间为20.839h,因为这时水的体积又为14800G ,也可确定水泵第二次停止运行时间为22.958h, 这是由于在紧接着的测量时刻23.88h 测得水的体积为663400G,与水泵停止时的水位67700G相比流掉了14300G。

计算出在相邻时间区间的中点及在时间区间内水箱中流出水的平均速度，并将其列成表3,并划出其散点图。

下面就对表3的数据进行拟合,最后对其积分计算出每天的用水总量。借助于软件包（Mathematica ）用9次多项式进行拟合，得

23456

576889()17.362311.52277.96127 2.7085 +0.4846130.04868160.0028516 9.6757810 1.7630110 1.3347510f t t t t t t t t t

t −−−=−+−−+−×+×+×拟合的复相关系数为0.980767,t 的定义由下式给出.

000024()0.4606,()0.4606,

T T T T t T T T T +−+<⎧=⎨−+≥⎩当时.当时.

五、模型求解

由上面的“分析余建模”中得到水的流速函数f (t )在t = 0.46 (h)和t =24.46(h)时的值分别近似为14170 G/h 和13730 G/h ,相差仅3%，从而可认为得到了f (t )的连续性。于是，一天里的用水总量近似地等于曲线f (t )在24小时周期内的积分

这个数字近似对应于水箱总体积的1.5倍。因为按常规每1000人的用水量为105000(G/d),其中d为天，因此估计这个区大约有3200人。240.46060.4606()338000(G)

f t dt +=∫

(3) 任意时刻从水箱中流出水的流速都可通过多项式模型计算出来。

(4) 用9次多项式曲线拟合所给的数据其复相关系数为0.981, 且回归值与原始数据点没有很大的波动。

(5) 在24小时周期的端点,模型的取值非常接近,可以推测几天的流速。

(6) 人们自然会将用水量与用电量联系起来,特别是对家庭不用天然气的情况。例如，烧饭需要大量的水来洗碗。同时烧饭和照明等也要用电,洗澡的耗水量也很大,并伴随着电力的消耗(热水器、电吹风等)。通过调查比较得到,用水分布类型是及其相似的。