储油罐的变位识别与罐容表标定(数学建模论文)

摘要

加油站都有若干个储存燃油的地下储油罐，需要采用流量计和油位计来测

量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。但是许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化，从而导致罐容表发生改变。需要定期对罐容表进行重新标定。在求解过程中，我们对于罐体无变位、罐体产生纵向变位、罐体在水平和纵向都产生变位三种情况，利用解析几何的方式计算出体积与变位参数之间的关系，同时应用契比雪夫多项式对体积值进行近似多项式展开用以对标高和出油量的关系进行拟合表示，得到较为满意的效果。

第一问、(1)针对无变位情况，我们计算得到椭圆油罐容积表达式为：

hbhb2hb'V椭1()arcsinabl，利用契比雪夫多项式方法在提高拟合精

vvv2度的前提下用5阶多项式拟合处标高和容量之间的函数关系；(2)对于纵向变位的情况，当椭圆型罐体发生变位纵向变位角度4.1时,我们利用体积等效思想，讲上述罐内不规则油量容积的计算转为(1)中规则油容进行计算，利用附件

(1)中数据利用最小二乘拟合方法算出油位高度的真实值，继而利用拟合多项式：

5432 V椭变0.0320H -1.1361H 13.2498H-52.5322H395.7748 H-408.5976进行间隔为1cm的此罐容表进行标定，得出的表标定值如下： 1cm 2cm 3cm 4cm  118cm 119cm 120cm  0L 0L 0L 0L 4017.26L 4050.08L 4082.80L 第二问、(1)利用第一问中等体积的思想，我们同样可以对纵向倾斜角度和横向倾斜角度时进行数学模型的建立。(2)在模型的建立过程中得到一个关于浮游子高度H和偏转角、以及等效高度h之间的一个表达式，从而利用

最小二乘拟合确定变位参数、。(3)利用已给数据求得表达式： 继而再次利用拟合拟合多项式得出间隔为10cm值： z0Rtanh12tan，

10cm 1352.17 20cm 2223.85 30cm 3082.13   280cm 60303.88 290cm 60690.98 300cm 61807.31 利用附表（2）中的数据进而进行模型正确性与可靠性的检验。

关键词：储油罐 罐容表 变位 契比雪夫公式 等效高度 最小二乘拟合

1

一、问题的重述

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

许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。储油罐的主体为圆柱体，两端为球冠体。

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

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

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

二、问题的分析

综述，该题要解决的问题是储油罐的变位识别和罐容表的标定。显然这是几何求解的应用问题。可以通过所学过知识求解。

对于问题一，题目中给出一个小椭圆型的平头油罐的模型和纵向的倾斜角的值，并且也给出了对这个罐体无变位和有纵向变位两种情况的实验数据，所以我们可以把这个问题分成两部分来解决，首先解决当罐体无变位时罐中储油量与油液面高度的关系，其次解决当罐体发生倾斜时罐体中油的储油量随油液面高度的变化。然后对这两种情况求出的结果进行比较，判断出油罐变位后对罐容表的影响。最后重新确定变位后罐容表的标定值。分析发现当油罐无变位时，问题比较好解决，我们可以利用积分的方法求出罐中的储油量和油液面高度的关系。当油罐变位时，首先要选好坐标系，确定倾斜后液面高度是怎样变化的。也是通过积分的方法进行求解，但是我们也可以寻找一点，可以把它等效成无变位的情况去解决。

对于问题二，给的是一个实际储油罐，其主体为圆柱体，两端为球冠体。该问题主要是要解决罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度 ）之间的一般关系。对这个问题可以看成两部分求解。首先也是要选好坐标系，确定好油面高度和变为参数的关系，分别求解主体圆柱体和两端球冠体中的储油量和油位高度的关系。建立出一定的数学模型，然后通过题目给出的数据，确定数据的可用性求出模型中的变位参数。最后通过误差分析检验模型的正确性和方法的可靠性。

2

三、问题的假设

1、忽略油罐中温度，压强对储油量的影响； 2、假设油罐的形状是规则的；

3、假设储容表的读数是无读数上的技术误差； 4、假设油罐变位的角度不会太大； 5、假设罐中的油分布均匀；

6、假设油罐不会因为出油的变化而发生位置变化；

四、模型的建立与求解

Ⅰ问题一：

对于本问题，首先可以把它分解成两个问题来考虑，第一个问题是确定小椭圆型储油罐没有变位时储油量与油面高度的关系，利用解析几何的方法确定标高和油量之间的关系。第二个问题确定油体容量和标高及变位角之间的关系，利用等效体积替换的思想把它转化为上一问中未变位的问题进行求解。同时我们分析发现变位数据的油高不是真实值，所以考虑利用已有读数寻求真值。依据上述思想分别对这两个小问题建立合理的模型。

㈠ 无变位时：

小椭圆型储油罐的横截面很明显是一个标准的椭圆，首先建立合适的坐标系，并且椭圆的图形如图1:

2b h 0 xx a 图（1） 1、符号声明： a---------------------椭圆罐横截椭圆面的半长轴

b---------------------椭圆罐横截椭圆面的半长轴 h---------------------无变位椭圆油罐中的油位高度

S椭无变位时椭圆的横切面面积

l---------------------小椭圆油罐的母线长

R2-----------------------方差和

3

Vi根据模型得出的数据

ˆ实测数据 Vi2、模型建立： 椭圆的方程为：

x2(yb)2122ab ------------------------------ -- (1)

式中： a-----------------椭圆罐横截椭圆面的半长轴且为0.米

b-----------------椭圆罐横截椭圆面的半长轴且为0.6米

油面高度为h时，储油罐横截面的面积可求得为：

'S椭2ah22bybdy0 ---------------------（2） b 式中：

h---------------无变位椭圆油罐中的油位高度； 可推导得：

2hbhbhb'1 s椭 abarcsinbbb（3）2------------

我们可以设thb，由0h2b可以得：1t1 b 可求得小椭圆油罐未变位时储油量随液面高度的变化为：

' V椭2 t1tarcsintabl -----------------------（4）

2 式中：

l---------------------小椭圆油罐的母线长且为2.45米；

在求解时我们可以利用最佳平方逼近多项式（即契比雪夫多项式,），求出储油量的近似值。根据公式（1）我们可以准确的计算出储油量，但是公式里有反三角函数没有多项式计算方便，所以我们用多项式逼近方法解决此问题。（具体理论证明见附录二）

求得 ：

2 V椭 t1tarcsintabl------------------------（5）

2'1 a0T0tanTntabl

n122 4

当k=1时， 求得V椭近似是一个2阶的多项式，当k=2时，求得的V椭近似

''是一个五阶的多项式。由于拟合时考虑的误差，选择五阶多项式拟合效果较好，所以我们用当k=2时的五阶契比雪夫多项式逼近V椭

'' V椭16615x80x348x5abl---------------（6）

21575 由于储油罐的形状以及温度压强的关系这里的系数不太准确，我们可以根据

题目（附件1中变位进油量数据）给出的累加进油量和油位高度的数据进行拟合。

用MATLAB软件拟合后得到的结果是:

（h）-0.0209h 50.6284h 4 - 9.2692h 376.1928h 299.4693h-3.9321 V 椭

对该模型进行误差分析（误差分析的数据源：附件1位变位出油量数据）：将数据代入到上述5阶拟合多项式中与已知油量数据对比验证： 验证数据方差和R2为：

'1nˆViVi Sn1i1221.0517 ---------------（7）

R12VˆVnii2Vi1i1niVi1 ----------------------（8）

2由以上可知拟合的效果比较好。

㈡、有变位时：

椭圆油罐有纵向倾斜且倾斜角4.1，选定的坐标系及有关参数如图2所示。分析发现变位后油高数据由于变位后产生影响，而油量数据我们认为是真实的，因此首先利用几何关系求出油量容积的解析表达式，利用最小二乘逼近的方法求出真实反映油量的油高读数。然后利用拟合的方法建立油高和油量的模型，从而解决标高问题。

5

y

d

c

L2: Q 

0

图（2） 1、 符号说明：

S2标准差

x 水平线 -------------纵向倾角

H-------------变位后液面在Y坐标轴上的分量

图中：

c=0.4m, d=2.05m;

设浮标的距离罐底的距离为hf，可以知道浮标的坐标为(c,hf） 直线AB的斜率为tan ，所以可求得AB的直线方程为： L2: ytanxhfctan ---------------------（9）

显然油液面的高度随横坐标x的变化为： Htanxhfctan -----------------------（10）

2abH 此时储油罐横截面的面积为： S椭0b2(yb)2dy---------------------------------（11）

2HbHbHbab 1 arcsin2bbb 所以变位时小椭油罐中油容量与油位高度的关系是：

V椭l0S椭dx --------------------------------------------（12）

根据公式（2）我们可以准确的计算出不同油位高的储油量，但是，公式中

还有反三角函数和开根号，我们计算起来很不方便，所以我们可以寻找一点做体积上的等效转换仿效上一问的方法。我们设那一点为Q，那么它的坐标可设为

6

,H，见图(2)。

'并且该点Q也在直线l1上，根据l1的直线方程可以求得H'为： H'tanhfctan ------------------（13）

对于一个固定外壳的容器，其中装定量体积的液体，无论如何对容器发生偏

转，其液体的体积始终是一个固定值。根据等效的原理可以肯定点Q可以把楔状油形体积转化为标准的椭圆柱的体积，并且此时的油位高度为H'。所以这样就可以转化成没变位的情况去解决问题。得到转化后的椭圆面积为： S'椭12aHb0'b2(yb)2dy --------------------（14）

2H'b'H'bHbab 1arcsin =bb2b 可以求得体积为：

2'''HbHbHb'abl---（15） V椭1 1arcsinbb2b

H'b 我们可以设t 得到：

b'''2' Vt1tarcsintabl -----------------（16）

2'椭1得到的体积公式的形式和未变位时得到的公式形式相似，所以我们同样可以 利用契比雪夫多项式来近似地求解储容量与油位高度的关系。 求得：

16'V615t'80t'348t'5 椭21575（17） abl -------------

题目中附录的高度是变位后的读数，为了求H'，我们可以用最小二乘法

' minV椭Vii1n2

算出的值。

7

讲求出来的带入H'中可以得到真实等价的油位高度H'，我们用这个真实等价的油位高度H'和题目给出的累加进油量的数据进行5阶多项式拟合，求得的V椭变近似的为：

5432 V椭变0.0320H -1.1361H 13.2498H-52.5322H395.7748 H-408.5976对这个模型进行误差分析根据上面的公式求得： 标准差为： S22.1602 R检验 R21

可见拟合的效果很好，符合要求。

然后根据这个模型就可以算出罐体变位后油位高度间隔为1cm的罐容表标定值如下(详细结果见附录)：

1cm 2cm 3cm 4cm  118cm 119cm 120cm  0L 0L 0L 0L 4017.26L 4050.08L 4082.80L 利用附件1中的有变位的出油数据分别带入上述两个模型中，做出的图如下： 从图中可以看出，当油位高度一样时变位后的油罐中储油量的值比没变位的储油量低。说明罐容表实际的读数比理论的要偏大。

8

图（3）

问题二

该问题是问题一的变形，是研究实际储油罐变位后储容量与油位高度之间的关系，储油罐的变位是横向和纵向都有一定的偏角，并且储油罐两端不是平头而是球冠体，首先建立合适的坐标系，并且要设定一定的参数，如图（4）: z R y h x 水平线 图（4）

（1）符号声明：

h hf ------------------------实际测量的油位高度  --------------------------储油罐横向偏转角度

R ----------------------------------------圆柱截面的半径 r --------------------------球缺的半径

b(1m) ---------------------- 球缺的顶高

V柱---------------------------柱体的体积

储油罐的中油的体积为：

VV柱V冠 ------------------------- （18） （一）、对储油罐圆柱体部分的计算

如图在XOZ坐标面内，Z轴方向的水深为

9

z h1  x 图（5） h1(hfR)cosR ---------------------------（19） 式中 ： hf ------------------------- 实际测量的油位高度  ------------------------- 储油罐横向偏转角度 R ----------------------------------------- 圆柱截面的半径 浮游子在ZOY面投影的点的坐标为（2，hf-R）；

ZOY坐标面与油液面的交线l2的斜率为tan，将浮游子在ZOY坐标面投影的点的坐标带入l2的方程得

zh1Rtany2 --------------------------(20) 即求得l2的方程得：

ztanyh1R2tan

对于这个计算我们也可以寻找一个点P,z0，使得z0+R变为没有变位时油面的高度，即

z0tanh1R2tan 在YOZ平面内截圆柱的方程为 z2y2R2 它的截面积为 S2z0RR2z2dz ----------------------------(21)

10

z0z0z02R 1arctan 22RRR

z0z0z021arctan V柱lS2 Rl------------------（22）2RRR' 由第一问可知它可以用5次多项式近似代换。

（二）、对储油罐两段球缺部分的计算

由球的关系有：

R rrb --------------------------------(23)

2222式中 ：

r ----------------- 球缺的半径 b(1m) ----------------- 球缺的顶高

解得 r=1.625m

在XOZ截面内用平面z=z0截到一个平面 面积为： Sz2

另p 有Szr2y2

prbr2z2rbr2z2y2dy - ------------------(24)

p2y2dy

p2p2rb1 =arcsinrbp2rb

242p2 V球z0-1.5Szdz

p2p2rb12 arcsinrbprbdz --------(25) 1.5242p2z0 这样储油罐的体积就可以表示为： VV柱V球 -------------------(26)

V也可以用多项式拟合，不过不能直接用附录里的数据，哪里的体积是不准

确的，它是没有考虑变位时对体积的影响，我们就用第一个体积，把它作为储油罐的初始体积，分别用当前行的体积减去下一行的出油量，这样我们认为此时的

11

体积即为实际体积，同第一问中问题二的解决方法一样通过最小二乘法计算出

和： 求得： 2.7 3.6

然后反向应用拟合出5次多项式，即：

V实a0h5a1h4a2h3a3h2a4ha5 ---------------(27)

5432 0.1206h-6.3084h  178.2758h 479.2440h 800.8409  -0.0015h

对这个模型进行误差分析求得： 标准差为： S20.8848 R检验 R21

符合要求，然后根据这个模型就可以算出罐体变位后油位高度间隔为1cm的罐容表标定值如下：

显示油量容积显示油量容积显示油高/cm 显示油高/cm

/L /L 10 1352.17 160 31568.85 20 2223.85 170 33378.42 30 3082.13 180 36168.99 40 3195.81 190 36925.2 50 4936.02 200 40631.3 60 6876.07 210 44270.94 70 91.21 220 46826.99 80 11258.52 230 49281.36 90 13656.67 240 50614.79 100 16165.76 250 52806.71 110 17767.13 260 55835.04 120 20443.18 270 59675.99 130 23177.17 280 60303.88 140 25953.09 290 60690.98 150 28755.38 300 61807.31

五、模型评价

12

地下储油罐的使用在现实中是比较常见的一种应用方式，而且出现罐体的偏转也是正常现象，所以我们建立的模型在一定程度上是可以解决一些社会现实问题的，具有一定的实用性。

本模型在适当的假设下，找到适当的角度，运用数据拟合的思想，建立数学模型，使用问题得以求解。在模型的建立过程中，我们用Matlab语言对所给的数据进行拟合，得出偏转角测量高度H、实际高度h、和与罐体中的油量的体积关系。从而求解出题目中的各个问题的答案。

以上我们所建立的数学模型，能比较好的解决题目中的要求，但我们的模型是建立在一定的条件下的，对于实际问题具有一定的局限性，不太能直接应用于现实。

参考文献

[1]、李岳生,黄友谦. 数值逼近,北京:人民教育出版社:281～283，1978 [2]、樊映川等编. 高等数学讲义下册,北京:高等教育出版社,1993 :62. [3]、袁亚湘，非线性优化计算方法，北京，科学出版社：100～120，2008 [4]、张志涌，等编著，Matlab 教程，北京航天航空大学出版社,2004.

13