储油罐的变位识别与罐容表标定

《储油罐的变位识别与罐容表标定》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、2010高教社杯全国大学生数学建模竞赛题目A题储油罐的变位识别与罐容表标定通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况.许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变.按照有关规定，需要定期对罐容表进行重新标定.图1是一种典型的储油罐尺寸及

2、形状示意图，其主体为圆柱体，两端为球冠体.图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图.请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题.（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为a=4.1°的纵向变位两种情况做了实验，实验数据如附件1所示.请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值.（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储

3、油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b）之间的一般关系.请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值.进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性.附件1：小椭圆储油罐的实验数据附件2：实际储油罐的检测数据油油浮子出油管油位探测装置注油口检查口地平线2m6m1m1m3m油位高度图1储油罐正面示意图油位探针9油位探针α地平线图2储油罐纵向倾斜变位后示意图油油浮子出油管油位探测装

4、置注油口检查口水平线图3储油罐截面示意图（b）横向偏转倾斜后正截面图β地平线垂直线油位探针（a）无偏转倾斜的正截面图油位探针油位探测装置3m(b)小椭圆油罐截面示意图α油油浮子出油管油位探针注油口水平线2.05mcm0.4m1.2m1.2m1.78m(a)小椭圆油罐正面示意图图4小椭圆型油罐形状及尺寸示意图92010年全国大学生数学建模竞赛A题解法评述该问题是来自于加油站设备研究与生产企业的一个实际课题，问题由两大部分组成：（1）为了观察检验罐体变位对罐容表的影响，在已知变位参数的情况下，检测出油位高度和油量的对应数值，建

5、模分析罐容表的变化规律，并给出修正的罐容表，属于“正问题”.（2）根据实际检测数据，正确识别罐体是如何变位的，具体变了多少？同时要给出罐容表的修正标定方法和结果，属于“反问题”.具体解法如下：第一部分：小椭圆型实验罐的有关问题按如下方式引入空间直角坐标系OXYZ（左手系）：以油位探针所在直线与罐体底部的交点为原点O；以罐体经过原点O的高为X轴，X轴正向指向罐体右端；以油位探针为Z轴，h增加的方向为Z轴正方向.符号说明：a、b分别为罐体截面椭圆的长、短半轴（a=0.89m,b=0.6m）；L为罐体长度，L1、L2分别为罐体左

6、、右截面到油位探针的距离（L1=0.4m,L2=2.05m,L=L1+L2=2.45m）；为罐内的油位高度；a为罐体纵向倾斜角度；为罐体纵向倾斜a角时、油位高度为h时的实际罐容量的理论值；为对应的实际值；为罐体纵向倾斜a角、油位高度为h时，截面X=x处油面到罐底的距离；简记为，简记为.于是有（1）无变位情况下，，，容易求得：代入参数的具体值，可计算得到无变位情况下的罐容表（理论值，见表1）.表1：无变位情况下小椭圆罐的罐容表部分结果油位高度/cm油量/L油位高度/cm油量/L油位高度/cm油量/L油位高度/cm油量/L10

7、163.59401199.31702489.151003659.8820450.27501621.00802910.841103946.5530803.54602055.07903306.611204110.15计算表1的MATLAB脚本：9a=8.9;b=6;L=24.5;V0=@(h)L.*(pi/2*a*b+a/b.*(h-b).*sqrt(2*b.*h-h.*h)+a*b.*asin(h./b-1));H=(0.1:.1:12)';D0=V0(H);评注1：对附件一在无变位情况下进油（或出油）的油位高度与实际罐容量

8、数据作3次样条插值，并与理论值比较，可以得到修正公式，但是这一步评阅要点没有要求。（2）罐体纵向倾斜a角时，仅考虑a>0.容易看出（）：现在，，.可以分三种情况讨论：情况一，当时，即，有此时，有(1)情况二，当时，即，有此时，有(2)情况三，当时，即，有此时，有(3)记函数9用MATLAB的数值积分函数