黄河小浪底调水调沙问题.pptVIP

黄河小浪底调水调沙问题 主讲人 陈玲华,靳文进,方  内容摘要 为了确定排沙量与时间、排沙量与水流量的函数关系,我 们可以用 Matlab软件做线性回归得到排沙量与时间的函 数关系式,再利用所求函数在区间[0,24上进行积分得到 总排沙量1.93962亿吨。对于排沙量与水流量之间的关系 按时间分为两段进行拟合,最终用 Matlab软件来画出图 像,确定排沙量与排水量之间的函数关系式。  问题的提出 ·在小浪底水库蓄水后,黄河水利委员会进行了多次试验 特别是2004年6月到7月进行的黄河第三次调水调沙试验 具有典型的意义。这次试验首次由小浪底、三门峡和万家 寨三大水库联合调度,进行接力式防洪预泄放水,形成人 造洪峰进行调水调沙试验成功。这次试验的—个重要目的 就是由小浪底上游的三门峡和万家寨水库泄洪,在小浪底 形成人造洪峰,冲刷小浪底的库区的沉积泥沙。在小浪底 开闸泄洪以后,从6月27日开始三门峡水库和万家寨水库 陆续开闸放水,人造洪峰于29日先后到达小浪底,7月3 日达到最大流量2720m^3/s,使小浪底水库的排沙量也 不断的增加。表一是由小浪底观测站从6月29日到7月10 日检测到的试验数据。  表1:试验观测数据单位:水流为立方米/秒,含沙量为公斤/立方米 日期 7.1 7,4 时间8:002008020.008:0020008002008:002000800200 水流量18001900210022002300240025002600265027027202650 含沙量3260 90 102 115116 日期 7.5 7.7 7.8 7.9 7.10 时间8:0020:0080020:008:020:008:0020:008:0020:008:0020:00 水流量2602500230022002001850182018001750150010001900 含沙量118120  向题分析: 1、对于问题一,所给数据中水流量X和含沙量h的乘积即 为该时刻的排沙量y即:y=hx。 2对于问题二,研究排沙量与排水量的关系,从实验数据 中可以看岀,开始排沙量随水量增加而增加,而后随水流 量的增加而减少,显然变化关系并非线性的关系,为此, 把问题分为两部分,从水流量增加到最大值为第一阶段, 从水流量最大值到结束为第二阶段,分别来研究水流量与 排沙量之间的函数关系  模型假设 1、水流量和排沙量都是连续的,不考虑上游泄洪 所带来的含沙量和外界带来的含沙量。 ·2、时间是连续变化的,所取时间点依次为 1,2,3,…,24单位时间为12h  时间与含沙量关系 时间与水汽量关系 时间2h 通过观察图像,我们可以看出其变化并不光滑,而 且也没有特定的表现出服从某种分布的趋势。  但是为了得到具体的计算函数,我们就必须对数 据进行拟合,所以通过 Matlab先利用spne方法对数 据进行插值,从而提高精确度,使图像变得光滑,然 后利用多项式进行拟合,当多项式次数越高拟合也越 准确,但是由于数据受到的影响较多,所以这里的数 据也不是准确值,因此我们可以只取三次进行拟合 也方便了后续的计算。 于是我们分别对含沙量和水流量进行插值拟 便可以得到下面图像和结果  时间与含沙量关系拟合图 12+y=001423-1.3x2+21*+16 如 所得到的拟合函数为 y=0.014*x^{3}-1.3*x^{2}+21*x+16 即含沙量与时间的关系式为 S=0.014*tA3-1.3*t^2+21*t+16