2010年我校参加全国大学生数学建模竞赛获全国二等奖论文A题

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

【摘要】

本文针对油罐的变位识别和罐容表标定的具体问题，通过建立模型将储油分成无数个小微元，研究小微元的规律从而利用微积分的方法推导出，当纵向倾斜通过建立数学模型将储油罐中剩余油的体积与?角度和横向偏移?角度一定时，

油位探针显示高度的函数关系表示出来，并制定罐容表标定值。

针对问题一：研究罐体变位对灌容表的影响，首先利用了逼近法计算油罐未变位前的储油体积，在MATLAB中对实际与理论数据进行拟合利用误差分析公式

E?理论值?真实值真实值，求得其误差E范围为0.014866?E?0.034917，求储油罐

中储油量时将油罐分成多部分考虑，利用微元思想和积分方法求得其储油量的体积与油位探针的读数h及变位角?，?的函数关系V?f?h,?,??在此问中

??0,??4.1?时得出V?f?h,4.1?,0?，求出其一定h时的V。用模型求出的理论

值与题目附表1中的实际值相比较，得出其误差1.42%?E?3.37%。并标出变位后间距为1cm罐容表。

针对问题二：对实际储油罐建立罐体变位罐容表的标定模型。在问题一的理论和方法的基础上加上球冠中的储油体积即可得到实际储油罐储油体积，并采用最小二乘法推导出所求变位参数?及?。并得出当?及?一定时油位探针的高度与剩余储油体积的关系V?f?h?，进而制定间距为10cm罐容表标定值。 本文充分运用了数学分析、高等数学等知识对储油体积积分，并通过MATLAB软件模拟的方法对理论数据进行了误差分析，以及运用最小二乘法估计其变位参数值?，?。最后对模型的优缺点进行了评价，并给出了改进方向。

关键词：MATLAB数据拟合、微元法、最小二乘法、罐容表标定

1

1 问题重述

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

通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。

许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示图3是罐体横向偏转变位的截面示意图。

请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。

（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为??4.1°的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。

（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度?和横向偏转角度? ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。

2 模型假设

[1] 假设油浮子不会超过油位探针颈口，在要到油位探针劲口时就停止注油； [2] 假设温度和压强对油的体积不影响；

[3] 假设油浮子，油位探针和出入油导管的体积不计；

[4] 假设问题一中的储油罐为标准的椭圆柱体，没有凹凸地方；

[5]假设储油罐只沿着题目中所示的纵向倾斜，当其沿反方向倾斜时其处理方法

类似，在此题中我们不予考虑； [6] 假设纵向倾斜变位的角度??9.5?。

3 符号说明

? 油罐纵向倾斜的角度

? 油罐横向偏转的角度

2

h 油位探针的显示高度 L 油罐的长度 V 储油的体积 a 椭圆的长半轴 b 椭圆的短半轴

4 问题的分析

本文主要研究油罐变位对罐容和罐容表的影响，首先我们利用微元思想与体积积分的算法建立油面高度与体积关系的模型。然后对问题一利用附件1的油位高度所对应的实际油量的体积与模型求出的理论油量的体积进行数据拟合。比较拟合曲线得到相应的误差范围。

在问题一里面由于椭圆柱的纵向偏离4.1°,罐容和罐容表都将受到影响。我们将倾斜的椭圆柱体分成三部分，并对每一部分通过微元和积分的思想和方法对其进行计算。从而推导纵向变位（倾斜）对罐容的影响，并标定出以1cm为间隔的罐容表。

在问题二里面我们利用第一题的思想方法，纵向倾斜为例，将球冠和圆柱都沿其轴线切截，从而油面形成弓形，利用微积分，分别对两边的球冠和倾斜的圆柱积分，从而得出储油的体积。

5模型的建立与求解

本节主要研究变位对椭圆柱油罐罐容表影响，分别对它在变位前后油罐中油体积的计算，从而推导倾斜对罐容表的影响，标定间隔为1cm的罐容表。

由于油罐是不规则容器，因此油罐中储油量的体积可用微元法的思想来求解。 微元法求几何体体积的基本原理

如图1所示的一个空间立体物体，设它在x 处S(x) 截面面积为S(x)，利用极坐标下推导面积公式的思想求出它的体积？

如果纵向切把它切成薄片，则每个薄片可近似看作直柱体，其体积等于底面积乘高，所有薄片体积加在一起就近似等于该立体的体积。

a 在[a, b]上任取一点x，并且任给x的一个⊿x 增量?x，这样就得到一个非常薄的薄片，这个小

薄片我们可以近似地把它看成柱体，于是这个微小的柱体体积为：

b dv?s?x??x?A?x?dx 图1 把这些小体积加起来，就是我们要求的体积即V??S?x?dx

ab5.0油罐中理论油量与实际油量的数据拟合

3

X2Y2如图2所示建立坐标系，设截面椭圆的方程为2?2?1,椭圆弓形的高

ab为h，图2（储油截面图）中带阴影部分为储油截面，用定积分求椭圆的体积。

设弓形的面积为S(h)则

2??h?b2ah?b2h?bh?b???2?ab (1) S?h??b?ydy???1????arcsin??bbbb??2?b???h?b?由0?h?2b,可知?1?x?1?油罐的长为L，储油的体积为V?x?，可得设x?b???V?x????x1?x2?arcsinx?abL (2)

?2?

图2 储油截面图

利用公式（2）所求的结果和附件1的油面高度得到对应的储油量体积，将这些点拟合得到理论的储油量变化曲线，并与题目给出的实际的储油量的拟合曲线比较，如图3所示。

450040003500300025002000150010005000-500 0200400600800100012001400理论数据理论拟合曲线实验数据实验拟合曲线

图3 无变位时的数据拟合图

4