储油罐的变位识别与罐容表标定数学建模竞赛.doc
《储油罐的变位识别与罐容表标定数学建模竞赛.doc》由会员分享,可在线阅读,更多相关《储油罐的变位识别与罐容表标定数学建模竞赛.doc(31页珍藏版)》请在三一办公上搜索。
1、全 国 大 学 生数 学 建 模 竞 赛储油罐的变位识别与罐容表标定参赛学校:重庆工商大学2010高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性如有违反竞赛规则的行为,我们将受到严肃处理我们参赛选择的
2、题号是(从A/B/C/D中选择一项填写): A 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名): 重庆工商大学 参赛队员 (打印并签名) :1 王文姣 2 白 洋 3 吴 静 指导教师或指导教师组负责人 (打印并签名): 袁德美 日期: 2010 年 9 月 13日赛区评阅编号(由赛区组委会评阅前进行编号):2010高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编):储油罐的变位识别与罐容表
3、标定摘要油品的数量管理在油品的经营过程中占有很重要的地位,其中储油罐罐容表的标定是加油站中油品管理的关键.但由于储油罐的长时间使用会导致地基变形,罐体的位置会发生纵向倾斜和横向偏转等变化(以下称为变位), 从而需要定期对罐容表进行重新标定.因此能够正确地解决好罐容表的标定问题,将会给现实生活中加油站等储油行业的操作带来方便.本文主要解决储油罐的变位识别及罐容表的标定问题.我们根据积分“无限细分,无限求和”的思想,通过建立积分模型,将储油罐划分为无数个连续的椭圆形截面.在进行储油量的计算时,由于油液面将这无数个椭圆截成了无数个弓形,故计算储油量的过程即转化为了对这无数个弓形在一定范围内求积分的问
4、题.问题一,在准确的模型假设的前提下,根据油位高度与各弓形面积的关系和弓形面积与油罐体体积的关系,分别对罐体无变位和变位的情况建立积分模型,然后利用附件的实测数据,对模型进行误差分析与拟合修正,最后给出罐体变位后油位高度间隔为1cm的罐容表标定值(结果请见表1).问题二,在问题一的基础上,首先我们同样采用积分的思想求得罐内储油量与油位高度及变位参数(纵向倾斜角度a和横向偏转角度b )之间的一般关系.然后根据对问题二的模型所求得的数据与附件所给的实际检测数据进行运算可以得到理想的a、b值,我们求解得出,.进而利用,得到油位高度间隔为10cm的罐容表标定值(结果请见表2). 另外在去掉温度对储油量
5、不会产生影响的假设条件下,我们对模型进行了进一步的改进. 为了消除温度的影响,我们考虑了油品的体积随温度变化的关系.利用经验公式.将油品体积全部转化为固定温度下的数据,然后再进行比较分析.关键词:优化处理;拟合;罐容表标定;微积分模型;最小二乘法.一、问题重述通常加油站都有若干个储存燃油的地下储油罐,并且一般都有与之配套的油位计量管理系统,采用流量计和油位计来测量进/出油量与罐内油位高度等数据,通过预先标定的罐容表(即罐内油位高度与储油量的对应关系)进行实时计算,以得到罐内油位高度和储油量的变化情况许多储油罐在使用一段时间后,由于地基变形等原因,使罐体的位置会发生纵向倾斜和横向偏转等变化(以下
6、称为变位),从而导致罐容表发生改变按照有关规定,需要定期对罐容表进行重新标定在不考虑外界环境的影响下,现解决如下问题:1为了掌握罐体变位后对罐容表的影响,利用小椭圆型储油罐(两端平头的椭圆柱体),分别对罐体无变位和倾斜角为的纵向变位两种情况做了实验,得出实验数据并在所得数据的基础上建立数学模型,研究罐体变位后对罐容表的影响,并算出罐体变位后油位高度间隔为1cm的罐容表标定值2在实际情况下,罐体变位后标定罐容表的标定值与理论上是有偏差的,但也存在着一定的联系,因此问题二需要找出罐内储油量与油位高度及变位参数(纵向倾斜角度a和横向偏转角度b )之间的一般关系在对实际情况下罐体变位后进/出油过程中的
7、实际检测数据进行分析与运算后,我们建立一数学模型,并通过其确定变位参数,同时求得罐体变位后油位高度间隔为10cm的罐容表标定二、问题分析储油罐罐体的变位识别是油位计量管理系统中的重要环节之一,而油品的数量管理是加油站等经营部门的基础工作,同时它又在其经营过程中占有重要地位目前,由于地基变形等原因,出现了一些不规范的问题故对罐体变位识别是确定一个规范的、科学的、精确的油位计量管理系统的必要前提问题一要解决的是小椭圆形罐体纵向倾斜变位后对罐容表的影响问题对于此类问题,我们通常利用高等数学中的定积分方法来求解其一般思想为“求和、取极限” 我们根据附件1所给出的小椭圆形罐体在无变位和变位时的进/出油量
8、与油位高度的实验数据最后来修正模型综上所述,先讨论小椭圆形罐体无变位时,储油量与油位高度之间的关系,建立积分模型一并且根据模型求出无变位时的罐容表然后再讨论当储油罐发生纵向倾斜后的情况,建立积分模型二模型二涉及二重积分的知识对模型二分盲区和非盲区两种情况进行讨论其中盲区包含两个部分:一、油面刚好接触油位探测装置底部,此油位探针的读数为0但实际油量不为0;二、油位探针刚好接触储油罐顶部,油位探针的读数为12,但此时储油罐并没有装满对于非盲区情况也需要进行分类讨论最后将模型数据和实测数据通过MATLAB软件进行拟合,我们可以得出两种情况下模型数据与实测数据间的关系,通过该关系进一步对原来的模型进行
9、修正最后确定变位后的罐容表,并给出罐体变位后油位高度间隔为1cm的罐容表标定值问题二要解决罐内储油量与油位高度及变位参数(纵向倾斜角度和横向偏转角度)之间的一般关系,且与未知,通过对题意的理解和对图形的分析,我们决定在问题一的基础上运用积分的知识建立数学模型三首先,我们将油罐体横向分为五个部分,并依次求得各部分截面面积;其次,我们又将油罐体纵向分为三个部分,依据之前求得的截面面积,纵向依次对其进行积分运算,从而得到各部分的体积,而油量的总体积即为各部分体积之和,该和式即为罐内油量与油位高度及变位参数与的关系式根据附件2所给出的数据确定与,然后通过对模型数据与实测数据之差(即离差)的平方和求出离
10、差最小时,与的取值,进而确定罐体在变位后油位高度间隔为罐容表标定值最后,再用附件2给定的数据,利用最小二乘法对我们所建立的“罐体纵、横向变位后模型”进行检验.下面为该问题的解法流程图:不变位变位问题一盲区非盲区变位问题二横向纵向特殊一般三、模型假设1 累计进/出油量与罐内油位高度为连续型变量;2 空气对油品的氧化情况不存在,注入油料时没有气泡的存在;3 地下储油罐的外界环境适宜如气压为常压,温度在19200,考虑到数据为8月份的数据,设温度为固定温度30;4 忽略储油罐壁厚和油浮子所占用的体积和罐底污泥厚度;5 系统稳定,不存在信号、噪声等外界因素带来的随机误差,也不考虑观测误差、连续问题离散
11、化所产生的误差,附录所给的数据真实、准确、可靠;6 该储油罐为两端平头且为椭圆的柱体;7 忽略温度对储油罐储油量的影响,储油罐储油量不随温度的变化而变化;8 储油罐密封性好,没有泄露和蒸发损失的情况;9 不考虑液体静压力对罐壁的作用而对油罐容积产生的影响;10储油罐罐壁平滑,不存在变形;11当高度达到1.2时,不再向储油罐内注油, 这是从单位经济效益方面考虑的12忽略油罐内部气体压强对注油这一过程的影响四、符号说明:储油罐截面圆圆心, ; :变位与无变位罐容表标定值的相似度;:储油罐纵向倾斜的角度,单位为度;:储油罐横向偏转角度,单位为度;:建立三维坐标轴,单位为;:建立三维坐标轴,单位为;:
12、小椭圆型油罐椭圆截面长半轴长,单位为;:小椭圆型油罐椭圆截面短半轴长,单位为;:小椭圆型油罐连续椭圆截面到储油罐罐底的距离,单位为;:以椭圆截面的中心为坐标原点,建立的横坐标,单位为; :第种情况下油位探针测得储油器的油位的高度, ,单位为;:在第问中第种情况下油罐在点处弓形截面高度,单位为;:第问中第种情况下油罐在阶段形成弓形截面面积与的关系, ,单位为;:在第问中第种情况下储油罐在第部分内的储油量关于的函数, ,单位为;:第种情况下求得的储油量,单位为;:第种情况下给出的储油量,单位为;:第种情况下求得的储油量的绝对误差,单位为;:第种情况下误差调节函数,单位为;:替换变量,单位为m;:储
13、油罐截面圆的半径,单位为;:球冠体球心到的距离,单位为;:储油器的油位的实际高度,单位为;:包含球冠体的球体的半径,单位为;:点纵坐标,单位为;:点的纵坐标,单位为;:储油罐各分段截面的面积,单位为;,:图形上相应的点;:图中相应圆的周长.五、问题一模型的建立与求解5.1 模型一的建立5.1.1 油罐无变位时模型的建立小椭圆型油罐无变位时,油位探针所测得的油位高度与椭圆截面的弓形高度始终是相等的,即此时,小椭圆型平头油罐椭圆截面的弓形面积如图1-1-1中阴影所示: 图1-1-1该椭圆的方程为: ,对阴影部分积分得弓形面积:,由图中弓形所形成的体的体积为: 5.1.2 罐体无变位时模型的求解利用
14、牛顿莱布尼茨公式求解得: (1.1)将给定的无变位时进油量的实验采集的数据和题中已知的数据代入式(11)中,用MATLAB编程求出模型一的结果,将其与给定的数据进行比较分析(程序见附录一)可得误差结果(见附录表1-1)5.1.3 误差分析及修正从附录表1-1中可以看出,绝对误差值随着储油量的增大而增大经分析产生误差的因素有:1油品中的气泡当油品中混有气泡时,由于气泡具有体积,从而使油位探针的读数比实际的读数大,且随着油量的增大气泡的所占的体积也增大;2油品储油罐罐壁的厚度由于储油罐罐壁包括内壁和外壁,我们计算的体积包括壁的厚度所占的体积所以随着油容量的增加,壁厚所占的体积就增大,我们所测量的体
15、积与实际油量的容积差就增大3储油罐的变形储油罐的变形是指罐体壁的凹凸变形,无论是凹还是凸都会使油位探针的读数与实际值不符,当罐壁凹进去时,实际容量比油位探针的读数小;当罐壁凸出来时,实际容量比油位探针的读数大在本题中,由于误差随储油量的增大而增大,因此可以猜测为罐壁凸时的情况;4外界温度油品的性质与外界温度有必然的联系,当外界的温度越高时,油的体积就相对越大为较正误差,我们在MATLAB软件中对附录表1-1中所得出的绝对误差值与油量高度进行了拟合(程序见附录二),得出了校正误差的调节函数关系式如下:,所以得到较正后的函数为:.下图为对理论数据调节前、后的曲线与实际曲线的拟合图,图1-1-2所示
16、:图1-1-2从图中可以看出,修正后的理论数据与实测数据能很好的吻合用MATLAB编程(程序见附录三)求出无变位情况下油位高度间隔为1cm时罐容表标定值(见表一)5.2.1 罐体变位后模型的建立在上面模型的基础上小椭圆型油罐在地基变形的情况下,发生了纵向倾斜角的倾斜,我们建立三维坐标系以油罐身长的延长线作轴,以油罐左底面的纵向对称轴为轴,以垂直于平面过点作轴,如图1-1-4所示:图1-1-31考虑盲部分:由于储油罐发生纵向倾斜,导致储油罐存在有部分油料体积无法准确测得的情况这就是所谓的盲区情况进一步说:所谓的盲区是指由于液位计的选型和安装位置不同形成的无法测量的区域出现盲区的情况又分为两种:(
17、1)第一种盲区情况如图1-1-5(盲区一)所示: 此时,由不变位时模型中椭圆截面弓形面积公式易得: ,积分得阴影部分体积(即盲区一的体积)得:,(2)同理,可得如图1-1-6中阴影部分体积(盲区二):综上所述:即当满足或者时,测量出油位的高度是有误差的,为了减小误差我们有必要将盲区考虑到模型中去2接下来研究非盲区情况:根据图1-1-4进行分析,可以将非盲区在分为三个部分,这三个部分在图中之间(1)当时,即在之间的区域内:此时的椭圆截面弓形面积公式为: ,求得储油量的公式为 : ,(2)当时,即在之间的区域内:此时的椭圆截面弓形面积公式为: ,求得储油量的公式为:(3)当时,即在之间的区域内:总
18、之,综合盲区和非盲区情况,可以将整个储油罐的储油量分为五个阶段,得到如下结果: 5.2.2 罐体变位后模型的求解(1)盲区两种情况储油量的计算,利用MATLAB编程求解(程序见附录四),得到结果,.此模型的求解利用MATLAB编程(程序见附录五)将附件1中的变位进油量的实验采集的数据导入,将得出的结果与实际结果进行比较分析可得误差结果(见附录表1-2)MATLAB编程进行误差拟合,得到此模型的误差拟合曲线,即调节函数:对此模型同样用MATLAB编程(程序见附录三)求出变位情况下油位高度间隔为1cm的罐容表标定值如下表一表一:无变位和变位时罐容表标定值高度h/cm无变位时的标定结果(L)倾斜变位
19、时的标定结果(L)高度h/cm无变位时的标定结果(L)倾斜变位时的标定结果(L)高度h/cm无变位时的标定结果(L)倾斜变位时的标定结果(L)01.7111.854412027.8961843.298822891.5472746.86616.4085.358422070.0111886.612832930.6222788.600215.4258.112432112.1021929.996842969.4072830.084327.20511.786442154.1581973.411853007.8872871.388441.18416.610452196.1682016.835863046.
20、0462912.412557.03222.515462238.1202060.289873083.8692953.186674.52529.699472280.0012103.743883121.3402993.671793.49438.123482321.8022147.167893158.4413033.8758113.80847.937492363.5082190.591903195.1563073.7599135.36059.211502405.1092233.946913231.4663113.34310158.06171.925512446.5922277.260923267.35
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 储油罐 变位 识别 罐容表 标定 数学 建模 竞赛
链接地址:https://www.31ppt.com/p-4062383.html