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

黄河小浪底调水调沙问题
一、问题的提出
2004年6月至7月黄河进行了第三次调水调沙实验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙实验获得成功，整个试验期为20多天，小浪底从6月19日开始预泄放水，直到7月13日结束并恢复正常供水。小浪底水利工程按设计拦沙量为75.5亿m3，在这之前，小浪底共积泥沙达14.15亿t。这次调水调沙试验一个重要目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙，在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于6月29日先后到达小浪底，7月3日达到最大流量2700m3/s，使小浪底水库的排沙量也不断增加。表1是由小浪底观测站从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现在，根据试验数据建立数学模型研究下面的问题：
给出估计任意时刻的排沙量及总排沙量的方法；
确定排沙量与水流量的关系。
二、模型的建立与求解
2.1问题一的模型
1、观测时间（时刻）的确定
以6月29日0时开始计时，各观测时刻（离开始计时的时间）分别为：，
其中，计时单位s。
2、排沙量的确定
记第次观测时水流量为，含沙量为，则第次观测时的排沙量。其数据如下表2。
表2时刻对应的排沙量排沙量单位：102kg
节点1234567889101112时刻t1t2t3t4t5t6t7t8t9t10t11t12排沙量57611401575187020702352250026522862302431283074节点131415161718192021222324时刻t13t14T15t16t17t18t19t20t21t22t23t24排沙量30683000271423101600111091054045530080453、模型建立
在上述已经知道24对数据的基础上，建立任意时刻的排沙量的函数，可以通过插值或拟合的方法来实现。考虑到实际中的排沙量应该是时间的连续函数，顾采用三次样条函数进行插值。
在求出三次样条函数的基础上，通过积分可以得到总的排沙量为：
。
4、程序
wv=[180019002100220023002400250026002650270027202650...
26002500230022002000185018201800175015001000900];
sth=[326075859098100102108112115116118120118105...
80605030262085];
i=1:24;
t=(12*i-4)*3600;
y=wv.*sth;
t1=t(1);
t2=t(end);
pp=csape(t,y');%或采用三次B样条插值：pp=spapi(4,t,y');两种结果一样
xsh=pp.coefs%求得插值多项式的系数矩阵，每一行是一个区间上的多项式系数
TL=quadl(@(tt)ppval(pp,tt),t1,t2)
5、结果：
xsh=[

-6.389389

6.622972162971060e-011-7.270705780568926e-0065.330865266189391e-0012.352000000000000e+005

-4.092849612657905e-0112.663045325228995e-0064.474499496724854e-0012.652000000000000e+005

-5.873884115581812e-0111.860563831575575e-0073.423238799069683e-0013.024000000000000e+005

-5.690704384073027e-011-1.342536338439519e-0066.792363910504279e-0033.068000000000000e+005

-1.77396058190