黄河小浪底调水调沙问题
更新时间:2023-04-17 00:00:01 阅读量: 实用文档 文档下载
黄河小浪底调水调沙问题
摘要:本文利用插值拟合的方法通过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: 试验观测数据单位:水流为立方米/ 秒,含沙量为公斤/ 立方米
(1) 给出估算任意时刻的排沙量及总排沙量的方法;
(2) 确定排沙量与水流量的变化关系。
二、模型假设
1.假设所给数据客观准确的反应了现实情况
2.假设所给数据遵循一定规律变化,即是连续的
3.假设模型中不需要考虑一些外在因素
4.假设可将时间化为等分的时间点进行计算
三、符号说明
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*t^5+24.29*t^4-582.92*t^3+2866*t^2+35340*t+24000
由于这里的多项式次数过高,不便于计算和传播,所以我们可以对其再进行一次拟合,有下面结果,程序见附录%paisha.m 。
所以拟合后的函数为V= 95*t^3-5.5e+003*t^2+7.7e+004*t-3.2e+004,通过图像可以看出排沙量与时间服从正态分布,所以也可以化成的形式e 的指数形式进行拟合,这里就不再重复计算。
我们得到了拟合函数,下面就可以计算出这几天的总排沙量∫12×60×60×24095?(t 3?5.5e +003?t 2+7.7e +004?t ?3.2e +004)dt ,通过Matlab 编程可以计算出定积分,结果如下,程序详见附录%jisuan.m 。
即总含沙量为1.704亿吨。
下面我们对问题“(2) 确定排沙量与水流量的变化关系。”进行分析计算。以下所有相关程序见附录%paishui.m ,下面就不重复说明。
我们先利用Matlab 将排沙量和水流量的相关数据反映到图像中。
通过观察可以看出,其关系是分段的,所以我们按时间进行分段拟合,拟合原理同问题(1)相同,于是可以得到分段前后的拟合多项式。
y = - 7.5e-005*x^{3} + 0.43*x^{2} - 5.2e+002*x + 3.6e+004
y = 2.3e-005*x^{3} - 0.066*x^{2} + 1.9e+002*x - 1.9e+005
综上,我们可以得到排沙量与水流量的关系式为
- 7.5e-5*v^3+0.43*v^2-5.2e+2*v+3.6e+4 0<=t<9
V=
2.3e-5*v^3-0.066*v^2+1.9e+2*v-1.9e+5 9<=t<=24
六、模型评估
本模型的优点是:建模简单,方便计算,适用度广。但也有最大的缺点为:精确度较低。为了减少误差,我们可以通过增大模型中拟合多项式的次数。当然在日后的模型改进中可以加入误差评估系统,来对模型进行完善。
附录
%tuxing.m
T=1:24;
S=[32 60 75 85 90 98 100 102 108 112 115 116 118 120 118 105 80 60 50 30 26 20 8 5 ];
W=[1800 1900 2100 2200 2300 2400 2500 2600 2650 2700 2720 2650 2600 2500 2300 2200 2000 1850 1820 1800
1750 1500 1000 900];
subplot(2,1,1);
plot(T,S);
hold on;
plot(T,S,'.');
title('时间与含沙量关系');
xlabel('时间t/12h');ylabel('含沙量/公斤每立方米');
subplot(2,1,2);
plot(T,W);
hold on;
plot(T,W,'.');
title('时间与水流量关系');
xlabel('时间t/12h');ylabel('水流量/立方米每秒');
%hansha.m
T=1:24;
S=[32 60 75 85 90 98 100 102 108 112 115 116 118 120 118 105 80 60 50 30 26 20 8 5 ];
x=1:0.1:24;
y=interp1(T,S,x,'spline');
plot(T,S,'.',x,y);
title('时间与含沙量关系拟合图');
xlabel('时间t/12h');ylabel('含沙量/公斤每立方米');
%liuliang.m
T=1:24;
W=[1800 1900 2100 2200 2300 2400 2500 2600 2650 2700 2720 2650 2600 2500 2300 2200 2000 1850 1820 1800 1750 1500 1000 900];
x=1:0.1:24;
y=interp1(T,W,x,'spline');
plot(T,W,'.',x,y);
title('时间与水流量关系拟合图');
xlabel('时间t/12h');ylabel('水流量/立方米每秒');
%jisuan.m
syms t;
S=0.014*t^3-1.3*t^2+21*t+16;
v=0.13*t^3-14*t^2+2.4e+002*t+1.5e+003;
V=v*S;
simple(V);
syms t;
V=95*t^3-5.5e+003*t^2+7.7e+004*t-3.2e+004;
int(12*60*60*V,t,0,24)
%paisha.m
t=1:24;
V=0.0018*t.^6-0.365*t.^5+24.29*t.^4-582.92*t.^3+2866*t.^2+35340*t+24000;
plot(t,V);
title('时间与排沙量关系图')
%paishui.m
t=1:24;
v=0.13*t.^3-14*t.^2+2.4e+002*t+1.5e+003;
V= 95*t.^3-5.5e+003*t.^2+7.7e+004*t-3.2e+004;
plot(v,V,'.');
title('整理图')
figure;
t=1:9;
v=0.13*t.^3-14*t.^2+2.4e+002*t+1.5e+003;
V= 95*t.^3-5.5e+003*t.^2+7.7e+004*t-3.2e+004;
plot(v,V,'.');
title('前半段图')
figure;
t=10:24;
v=0.13*t.^3-14*t.^2+2.4e+002*t+1.5e+003;
V= 95*t.^3-5.5e+003*t.^2+7.7e+004*t-3.2e+004;
plot(v,V,'.');
title('后半段图')
正在阅读:
黄河小浪底调水调沙问题04-17
家庭给我带来了快乐作文700字07-01
论文综述报告05-29
教师队伍建设年工作开展情况03-09
自学考试电子商务(填空)含答案01-25
142页2015年公共基础知识真题演练2000道附答案05-01
组胚试卷10-25
结婚时预订婚宴酒店注意事项04-23
【鼓励病人的话】鼓励病人的话大全02-10
补锌产品有哪些 - 图文09-25
- 教学能力大赛决赛获奖-教学实施报告-(完整图文版)
- 互联网+数据中心行业分析报告
- 2017上海杨浦区高三一模数学试题及答案
- 招商部差旅接待管理制度(4-25)
- 学生游玩安全注意事项
- 学生信息管理系统(文档模板供参考)
- 叉车门架有限元分析及系统设计
- 2014帮助残疾人志愿者服务情况记录
- 叶绿体中色素的提取和分离实验
- 中国食物成分表2020年最新权威完整改进版
- 推动国土资源领域生态文明建设
- 给水管道冲洗和消毒记录
- 计算机软件专业自我评价
- 高中数学必修1-5知识点归纳
- 2018-2022年中国第五代移动通信技术(5G)产业深度分析及发展前景研究报告发展趋势(目录)
- 生产车间巡查制度
- 2018版中国光热发电行业深度研究报告目录
- (通用)2019年中考数学总复习 第一章 第四节 数的开方与二次根式课件
- 2017_2018学年高中语文第二单元第4课说数课件粤教版
- 上市新药Lumateperone(卢美哌隆)合成检索总结报告
- 小浪底
- 调水
- 黄河
- 问题
- 最新xxx-课程微课制作脚本模版word版本
- (完整版)动物疫病检测样品采集、重大动物报告及采样表填写注意事
- 传承井冈山精神坚定理想信念
- 部编七年级历史下册材料七年级历史下册分析题专项材料(含答案)
- 人教版四年级语文下册全册导学案集
- 一年级下册数学第5课时 找规律教案优秀教学设计
- 给妈妈的一封信作文【小学六年级500字】_6
- 小学作文:我的理想是当一名科学家.doc
- 教材解析:必修第一课《生活在人民当家作主的国家》
- 河南省洛阳市东升二中八年级英语上册《Unit 1 Where d
- 青岛市2022年高三统一质量检测数学参考答案
- 幼儿园大班音乐游戏教案《快乐的牛仔》中班音乐游戏教案
- Kl-水冷冷水机组使用手册(盟立)
- 安科瑞多回路配电监控装置在数据中心精密配电柜中的应用
- 提高传接球能力的教学训练讲课稿
- 宜川县职称论文发表-建筑施工绿色施工运用论文选题题目
- (人教版)(全国卷)2022届高考英语一轮复习 Unit 5 First aid课时
- 实现有质量有效益可持续的发展
- 100句经典古诗词名句及解释鉴赏名师优质资料
- 社保卡电子版照片要求及报送流程