2010年数学建模全国一等奖论文--储油罐的变位识别与罐容表标定.doc
《2010年数学建模全国一等奖论文--储油罐的变位识别与罐容表标定.doc》由会员分享,可在线阅读,更多相关《2010年数学建模全国一等奖论文--储油罐的变位识别与罐容表标定.doc(32页珍藏版)》请在三一文库上搜索。
1、1 储油罐的变位识别与罐容表标定 摘要 本文对储油罐的变位识别和罐容表的标定问题进行了深入探讨,建立了储油量和 油位高度以及变位参数之间关系的数学模型,主要应用了 mtalab 进行求解。 针对问题一,我们利用积分的方法推导出小椭圆储油罐在无变位和发生纵向倾斜 变位时的一般公式。讨论了在储油罐发生纵向倾斜变位后对罐容表的影响,定义了平 均影响率 (变位前后储油量之差绝对值的平均值占总罐体容积的比例)作为评价罐体 变位对罐容表的影响程度的大小的指标,求出 。并分别给出了小椭圆储油罐4.87% 在无变位和在纵向倾斜变位角取 的罐容表。4.1 表 1 小椭圆储油罐罐容表(纵向变位 ).1 油位高度
2、/hm0.00 0.01 0.02 0.03 0.04 0.05 储油体积 3V0 0.00170.0035 0.0063 0.0100 0.0148 0.0207 1.15 1.16 1.17 1.1713 1.18 1.19 1.20 3.9103 3.9338 3.9560 3.9588 3.9766 3.9955 4.1017 4.1101 针对问题二,将储油罐分为5个区域分别进行讨论,考虑到在球冠处的体积表达式 过于复杂,我们省略了球冠处的一小部分体积,进行了近似求解,得出了罐内储油量 与油位高度以及变位参数之间的一般关系的数学模型。 在利用储油罐的实际测量值估计变位参数时,我们建立
3、了最小二乘拟合模型,得 到了最佳的变位参数为:纵向倾斜变位 ,横向偏转变位 。并据此对2.164.50 储油罐的罐容表进行了标定(见表3)。 在模型验证中,我们又采用蒙特卡洛模拟的方法对在问题二的模型中忽略的部分 球冠体积进行了模拟计算。得到用问题二模型中求出的总储油量与模拟得出的总储油 量一致度达到了99%,误差非常小,验证了我们所建立的模型的合理性和准确性。 关键词 平均影响率 最小二乘参数估计法 蒙特卡洛模拟 2 一 问题重述 通常加油站都有若干个储存燃油的地下储油罐,并且一般都有与之配套的“油位 计量管理系统” ,采用流量计和油位计来测量进/出油量与罐内油位高度等数据,通过 预先标定的
4、罐容表(即罐内油位高度与储油量的对应关系)进行实时计算,以得到罐 内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后,由于地基变形等原因,使罐体的位置会发生纵向 倾斜和横向偏转等变化(以下称为变位) ,从而导致罐容表发生改变。按照有关规定, 需要定期对罐容表进行重新标定。 问题一 为了掌握罐体变位后对罐容表的影响,利用小椭圆型储油罐(两端平头的椭圆柱 体) ,分别对罐体无变位和倾斜角为=4.1的纵向变位两种情况做了实验。请建立数学 模型研究罐体变位后对罐容表的影响,并给出罐体变位后油位高度间隔为1cm的罐容表 标定值。 问题二 对于实际的储油罐,试建立罐体变位后罐内储油量与油位高度及变
5、位参数(纵向 倾斜角度和横向偏转角度 )之间的一般关系。利用罐体变位后在进/出油过程中的 实际检测数据,根据你们所建立的数学模型确定变位参数,并给出罐体变位后油位高 度间隔为10cm的罐容表标定值。然后进一步用实际检测数据来分析检验模型的正确性 与方法的可靠性。 二 模型假设 1、假设题中所给数据均为储油罐内壁测量值; 2、不考虑由于温度、压强变化等原因而引起储油罐的体积变化; 3、油位探针被固定在储油罐上,其上油浮子能够准确测量油位高度; 三 符号说明 椭圆的半长轴长a 椭圆的半短轴长b 储油罐的总长L 油位探针到油罐底部左侧的距离l 储油罐的总体积TV 储油罐的纵向倾斜角度 储油罐的横向偏
6、转角度 油位高度h 球冠体的半径R 圆柱体的底面半径0 注:未说明符号在文中用到时注明 四 问题一的解答 小椭圆储油罐罐体变位前后都可以应用积分的方法求出罐体的储油量和油位高度 3 之间的关系。对于纵向倾斜的小椭圆储油罐,考虑分段求出其储油量和油位高度之间 的关系,从而得到重新标定后的罐容表。 4.1 小椭圆储油罐无变位时的模型 由于此时的椭圆无变位,考虑先对二维椭圆进行积分。为方便表示油位高度,建 立如图所示的坐标系,椭圆的半长轴长为 ,半短轴长为 ,则椭圆方程为ab22()1 0xyb 图1 对椭圆的积分示意图 在 方向上取椭圆面中的一微元 积分得到油的侧面积ydy2() 2100()21
7、ybhahD ybsxxad 储油罐内油的体积为 20()21hybVsLad 查积分表得到 (1)2()arcsin()2bahb 利用matlab计算得到 152232()arcsiLhVb 经验证两种方法得到的体积公式完全等价,(1)式即为小椭圆储油罐无变位时的储油 量和油位高度关系的模型。 根据此模型,我们可以求出小椭圆储油罐无变位时罐容表标定值(油位高度间隔取 ,结果见附录一) 。1cm 4 4.2 小椭圆储油罐纵向倾斜变位时的模型 储油罐纵向倾斜之后,油位计在油位过高或者过低时将不起作用(如图2所示的 和 区域),考虑到倾斜角 变化一般不会很大,所以我们可以将储油罐按液面高低1v5
8、 分成五个部分 ,来求其储油量和油位高度之间的关系。我们讨论的是小椭圆储15v 油罐纵向倾斜变位为逆时针旋转,如图2。对于储油罐顺时针旋转变位(即 为负值)时 的情况与此非常类似,在此不再详细讨论。 图2 储油罐分区示意图 4.2.1 对区域 的讨论 1v 在区域 ,其油位低于油位探针的油浮子,所以油位计量系统中显示油位高度为 零。 当油位计刚开始有示数时,计算其储油体积。将区域 放大得到图31v 图3 区域 的放大图1v 图中,从原点纸面向里为 轴,利用三重积分可以得到x 5 2tan11tan002yblblDVdxyzdxdz 其中 为油位探针到储油罐左侧的距离 l 积分得到 (2) 2
9、1tanta21blylyVdb 4.2.2 对区域 的讨论2v 由区域 很容易得到区域 的储油量和油位高度的变化关系,直接给出结论:12v2tan12(tan)00yhlbbbhlDVdxyzdxdz 所以 (3) 22(tan)ta1bhlylbyd 4.2.3 对区域 的讨论3v 图4 区域 示意图3v 在小椭圆储油罐无变位模型中我们已经求出了 的计算公式,同区域 中的积分1v 原理可以计算出 ,我们就可以得到此时的油量体积为abv和 (4)3abv 其中 6 2()tantan1002ybhylhlavdxdz2()arcsi()2bLb 2()1tan()tan0ybhLb hylL
10、lvdxdz 4.2.4 对区域 的讨论4 由区域4和区域2的相似性,将(3)式中的 换为 ,将 换为 ,并用总h(2)bl()Ll 体积减去 即为区域4的储油体积和油位高度的变化关系。V 21.()tan11.2()tan00yhlbbbTThLlDdxyzVdxdz 其中 为小椭圆储油罐的总体积 化简并积分可得 (5) 24()tan()tan21tbThLlybhLlyV db 4.2.5 对区域 的讨论5v 在此区域中油浮子到达油位探针顶点,无法进一步测量油位高度。无法测量的总 体积为: (6)5 2()tan()tan21bLlVylydb 4.2.6 综合各区域的罐容表标定的数学模
11、型 综上所述,我们得到了储油量 和油位高度 、纵向倾斜角 之间的分段函数关h 系式: 表2 分段函数关系(,)V 区 域 油位高度 /hm储油量 3/Vm 1 0 21tan ta0, 21blylyVdb 7 2 (0, -)tanhLl 22(tan) ta1bhlylbyVd 3 (-)tan, 2hLlb22()tant13002()1tan()tan0 rcsi()2ybhyllhybhLhylLldxzbdxdz 4 (2-tan, hbl 24()tan()t21tbThLlylV yb 5 5 2()tan( ,()antT TbLl Vld 根据储油量和油位高度的分段函数关系
12、我们可以得到罐体纵向倾斜变位( )4.1 后油位高度间隔为 的罐容表标定值(见附录一)。1cm 当 时,各区域油位高度及体积变化范围为:4. 表3 各区域油位高度和储油量变化范围 区域 油位高度 /h储油量 3/Vm 1 0 10,.7 2 (,.147 2(.,.53 3 .,.3 h 30.1,.98V 4 (1.7,.2 4(.5,.27 5 . h.01,. 4.3 罐体变位后对罐容表的影响 为了能更加准确地刻画罐体的纵向倾斜变位对罐容表的影响,我们分别对罐体变 位前后的理论值和测量值进行多方面的比较。 4.3.1 罐体变位前理论值与测量值比较 根据附件一中所给数据,我们计算出在附件中
13、所给的油位高度下理论值和测量值, 8 并画出其曲线。 0 0.2 0.4 0.6 0.8 1 1.2 1.40 0.5 1 1.5 2 2.5 3 3.5 4 4.5 位 位 位 位 h/m 位位 位V /m3 位 位 位 位 位 位 图5 罐体变位前的 曲线对比V 通过对比我们发现,对于任意 ,储油量的理论值和实际值始终成如 (7)式的比例h 关系。 (7)0.963测 量理 论 4.3.2 罐体变位后理论值与测量值比较 同样根据附件一中所给数据,我们计算出在附件中所给的油位高度下理论值和测 量值,并画出其曲线,如图6。 9 0.4 0.5 0.6 0.7 0.8 0.9 1 1.10.5
14、1 1.5 2 2.5 3 3.5 4 位位位位 h/m 位位 位/ m3 位位位 位位位 图6 罐体变位后的 曲线对比V 从图6中可以看出测量值仍然始终小于理论值,进一步求得理论值与测量值之差的 变化范围为0.0454,0.0910,测量值的相对误差范围为 1.56%,5.18%。 4.3.3 罐体纵向倾斜变位 前后理论值比较4.1 0 0.2 0.4 0.6 0.8 1 1.2 1.40 0.5 1 1.5 2 2.5 3 3.5 4 4.5 位位位位h/m 位位 位位 位位 位V /m3 位位位 位位位 图7 罐体变位前后的 曲线对比V 10 图8 同一高度下储油量的理论值与测量值之差变
15、化关系 图9 储油量的测量值的相对误差随油位高度的变化关系 由以上各图可以清晰地看出纵向倾斜变位后,使得在同一个油位高度下,变位后 比变位前的储油量减小。但是这样仍不够直观,我们需要找到一个指标来定量刻画罐 体变位后对罐容表的影响。从图9中可以看出,当油位高度 较小时(0.1 附近),变位hm 11 后相对于变位前的相对误差几乎达到了60%以上,但是此时的储油量的差别并不大, 鉴于此,我们定义平均影响率: 0niV:( 变 位 前 ) ( 变 位 后 )( 总 ) 来刻画罐体变位后对罐容表的影响。可以求出在纵向倾斜变位 时,4.1 4.87%。 五 问题二的解答 如图 8 实际的储油罐示意图,
16、其主体为圆柱体,两端为球冠体,在储油罐无变位 时我们计算其各部分体积。 图 10 实际储油罐示意图 圆柱体积计算公式 2VLR圆 柱 带入数据得到:主体圆柱体积为 56.5487 3m 一端球缺体积计算公式 20H球 缺 带入数据得到:两端球缺总体积为 ,则储油罐的总体积为 64.664538.15m3m 5.1 考虑不发生变位时储油量和油位高度的关系 圆柱内的油体积随油位高度变化关系: (8) 220hLRVdh圆 柱 一端球缺内油体积随油位高度变化关系: 12 (9) 200 22coshhSrRHVdarRHdh 圆 000圆球 缺 圆 其中 22rRh0圆 为圆柱体底面半径 为球缺对应
17、的半径0 为球内小圆半径r圆 总的储油罐内的油量对油位高度的变化关系为: 2VhVh球 缺 圆 柱 220220cos hLRRHrarRHddh 000圆 圆 圆 5.2 只考虑横向偏转变位时的储油量和油位高度的关系 只考虑储油罐横向偏转为 时:对实际的油位高度没有影响,但此时的油位探针 已经随储油罐发生偏转(如下图) 图 11 只考虑横向偏转示意图 由油浮子测量得到的油位高度与实际油位高度的关系为: (10)cos Rh 所以 对罐容表影响转换公式:cosRhVh 13 5.3 只考虑纵向倾斜变位时的储油量和油位高度的关系 只考虑纵向倾斜变位时,我们利用问题一中的思想,将储油罐分成 5 个
18、区域,分 别计算储油量和油位高度的关系(如图 12 所示)。 图 12 实际储油罐分区示意图 5.3.1 对区域 1 的讨论 油位探针测得的油位高度 始终为零,分成三部分来计算这部分的体积,如区域h 1 的放大图(13) 图 13 实际储油罐分区示意图 区域 1 的总体积应为: 1acbV 对圆柱体部分进行三重积分得到: 2tan 2tan00tan2tanc yRlRRRl lyldyxdzRydV 由球缺部分的体积随油位高度的计算公式(9)得: tbVl球 缺 由于 部分的体积精确计算非常复杂,而且储油罐的纵向倾斜角度一般不会高a 于 所我们考虑将这部分体积省略,进行近似计算。由于倾斜角
19、较小,所以区域 35 14 占据了储油罐的绝大部分,而在区域 3 中这种近似计算的误差将由于左右两个球冠的 省略体积一正一负而有所减小,所以 0cV 由此得到区域 1 的总体积公式: (11) 2tan2tantanRlylVRydVl 球 缺 5.3.2 对区域 2 的讨论 在区域 2 中,油位高度 的变化范围 。h0tanLl, 各部分储油体积和油位高度的变化关系: 2tantan002t t2yRhlRRhllcVdyxdzly , bl球 缺 aV 总储油量和油位高度变化关系: (12)2tan2tantanRhlyhlRydV hl 球 缺 5.3.3 对区域 3 的讨论 在区域 3
20、 中,油位高度 的变化范围 。tan 2taLll, 圆柱体部分的储油量: 12cVhV圆 柱 其中 22tantan()100()2 tan()tanhylhlRyLhylhLlVdxdz 球缺部分的近似储油量: ,t tanbVlVl球 缺 球 缺 0aV 总储油量和油位高度变化关系: (13) 312t tnhhlhLl圆 柱 球 缺 球 缺 5.3.4 对区域 4 的讨论 15 在区域 4 中,油位高度 的范围 。h2tanRl, 圆柱体部分的储油量: 2()tan()ttRchLlyhLlVRyd 圆 柱 球缺部分的总体积: ,b V球 缺球 缺 0a 总储油量和油位高度变化关系:
21、(14)24()tan()tan2 tnRhLlyhlVRydhLlV 圆 柱 球 缺 球 缺 5.3.5 对区域 5 的讨论 在区域 4 中,油位高度 始终等于 。 各部分体积分别为: 2()tan()tan2RcLlylRVyd圆 柱 ,Vb球 缺球 缺 0a 总储油量和油位高度变化关系: 25()tan()ta2 2tanRLlylVRydRLlV 圆 柱 球 缺 球 缺 ()ta , t tTRllVl 圆 柱 球 缺 球 缺 5.4 综合考虑储油罐纵向倾斜和横向偏转 根据如上讨论,我们可以得出结论,可以直接把 5.3 中各区域的 公式用式Vh (10)进行变换即可得到综合考虑了储油罐
22、纵向倾斜和横向偏转的 一般关系(, ) 式。 cosRhVh 我们依然考虑将储油罐分成五个区域分别求解。 5.4.1 对区域 1 的讨论 当储油罐未发生偏转时,油位高度为 ,而当储油罐发生横向偏转后,就可能使h 得油位探针测得的示数变为零。 16 图 14 对区域一的讨论示意图 即当 coshR 时,储油罐发生横向偏转 后,油位探针测得的油位高度恰好为零。本属于第二区域 的部分横向偏转变位后转为了变位后的第一区域。在积分时只需要将 替 costanRl 换在第一区域的积分即可得到 的一般关系式。(, )V 在此区域 依然恒等于零, 的(因此区域 ,所以只是求其最大体积)一般h 0h 关系式为
23、2(1cos)tant0(1cos)tan2 1costanyRlRldxdzVVRl 球 缺 (15) 2tlyl球 缺 5.4.2 对区域 2 的讨论 由式(22) cosRh 可以得到 ()h 在第二区域须满足条件 ()tanLl 计算可得 t1cos0lRh 再由(12)式: 17 (16) 2tan2tan00() tanyRhlRRhlVdyxdzVhl 球 缺 然后只需要做下式的变换即可得到 一般关系式:(, )V22cosRhh 5.4.3 对区域 3 的讨论 同区域 2 中讨论, 应该满足h()tan1cosLlRhcsh2tanRl 得到: ()tan1cos1costan
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2010 数学 建模 全国 一等奖 论文 储油罐 变位 识别 罐容表 标定
链接地址:https://www.31doc.com/p-25886.html