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

黄河小浪底调水调沙问题
摘要：本文利用插值拟合的方法通过Matlab工具模拟出了排沙量与时间、排沙量与水流量的函数关系，并且求出了总排沙量为1.704亿吨。整个模型简单且方便计算，其中排沙量与水流量的函数关系为分段函数。

关键词：调水调沙Matlab插值拟合

一、问题重述

2004年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功．整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日恢复正常供水结束．小浪底水利工程按设计拦沙量为75.5亿立方米，在这之前，小浪底共积泥沙达14.15亿吨．这次调水调试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙．在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700立方米/每秒，使小浪底水库的排沙量也不断地增加．下面是由小浪底观测站从6月29日到7月10日检测到的试验数据：
表1:试验观测数据单位：水流为立方米/秒，含沙量为公斤/立方米
日期6.296.307.17.27.37.4时间8:0020:008:0020:008:0020:008:0020:008:0020:008:0020:00水流量180019002100220023002400250026002650270027202650含沙量326075859098100102108112115116日期7.57.67.77.87.97.10时间8:0020:008:0020:008:0020:008:0020:008:0020:008:0020:00水流量26002500230022002000185018201800175015001000900含沙量11812011810580605030262085现在，根据试验数据建立数学模型研究下面的问题：
(1)给出估算任意时刻的排沙量及总排沙量的方法；
(2)确定排沙量与水流量的变化关系。

二、模型假设

假设所给数据客观准确的反应了现实情况
假设所给数据遵循一定规律变化，即是连续的
假设模型中不需要考虑一些外在因素
假设可将时间化为等分的时间点进行计算

三、符号说明

t:时间或时间点
v:水流量
S:含沙量
V:排沙量

四、问题分析

假设水流量和含沙量都是连续的，那么某一时刻的排沙量V=v(t)S(t)，其中v(t)为t时刻的水流量，而S(t)为t时刻的含沙量。通过观察数据，这些数据是每个12小时采集一次，所以我们可以将时间设为时间点t，依次为1，2，3，……，24，单位时间为12h。为了找到排沙量与时间的关系，我们就要先找到水流量和含沙量与时间的关系，一但找到水流量和含沙量与时间的关系，那么所要求的问题也就不难解决了。

五、模型的建立与求解

通过分析，我们假设水流量和含沙量都是连续的，那么我们开始对问题“(1)给出估算任意时刻的排沙量及总排沙量的方法”进行求解。
我们通过Matlab工具将所知道的数据显示为直观的图像，如下所示，具体程序见附录的%tuxing.m。

通过观察图像，我们可以看出其变化并不光滑，而且也没有特定的表现出服从某种分布的趋势。
但是为了得到具体的计算函数，我们就必须对数据进行拟合，所以通过Matlab先利用spline方法对数据进行插值，从而提高精确度，使图像变得光滑，然后利用多项式进行拟合，当多项式次数越高拟合也越准确，但是由于数据受到的影响较多，所以这里的数据也不是准确值，因此我们可以只取三次进行拟合，也方便了后续的计算。
于是我们分别对含沙量和水流量进行插值拟合，便可以得到下面图像和结果，具体程序见附录%hansha.m和%liuliang.m。

所得到的拟合函数为：y=0.014*x^{3}-1.3*x^{2}+21*x+16
即含沙量与时间的关系式为：S=0.014*t^3-1.3*t^2+21*t+16

所得到的拟合函数为：y=0.13*x^{3}-14*x^{2}+2.4e+002*x+1.5e+003
即水流量与时间的关系式为：v=0.13*t^3-14*t^2+2.4e+002*t+1.5e+003
因为某一时刻的排沙量V=v(t)S(t)，所以我们可以将所拟合出来的多项式带入上式，通过Matlab进行计算可以得到下面答案，程序见附录%jisuan.m。
ans=91/50000*t^6-73/200*t^5+2429/100*t^4-14573/25*t^3+2866*t^2+35340*t+24000
即排沙量与时间的关系为：
V=0.0018*t^6-0.365*