储油罐的变位识别和罐容表标定的方法分析

承诺书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A

我们的参赛报名号为（如果赛区设置报名号的话）：

所属学校（请填写完整的全名）：中国计量学院

参赛队员(打印并签名) ：1. 章苹

2. 文一章

3. 余舒婷

指导教师或指导教师组负责人(打印并签名)：数模组

日期： 2010 年 9 月 13 日赛区评阅编号（由赛区组委会评阅前进行编号）：

编号专用页

赛区评阅编号（由赛区组委会评阅前进行编号）：

全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：

储油罐的变位识别和罐容表标定的方法分析

摘要

本文研究了储油罐的变位识别和罐容表标定的问题。储油罐在使用一段时间后会发生变位，研究储油罐的变位识别和罐容表标定是加油站在管理中要进行的必要工作。

对于问题一，我们利用积分计算体积的方法，建立了罐内油量随油位高度和纵向倾斜角度变化的函数关系式。在对罐容表进行间隔为1cm 的标定过程中，由上述关系式计算得到的值中有很少的几个为虚数，我们对根式内的表达式加上1010-，使得所有计算值均为实数，并且原来的实数值在精度为410-时不产生变化，在这种情况下得到了无变位和倾斜变位时的理论标定值。接着，我们在分析理论值与实验值误差的基础上，对理论标定值进行了三次修正，依次为进出油量实验值与理论值误差的修正，罐内油量实验值与理论值误差的修正，以及倾斜无变位标定值中异常值的处理。修正后可以得到无变位和倾斜变位的最终罐容表标定值。通过分析储油罐纵向倾斜变位的示意图，我们得到了罐内油量一定时，纵向倾斜角对变位时油位高度与无变位时油位高度差值的影响：当纵向倾斜角为顺时针旋转时，角度增大，变位时油位高度与无变位时油位高度的差值会增大，从而相同高度下倾斜变位的罐容表标定值会减少。考虑到罐体变位对罐容表造成的影响，我们对装置进行了改进，将油位探针的位置移至储油罐中心轴处，消除了罐体变位对罐容表的影响。

对于问题二，首先我们用积分计算球冠形储油罐变位的体积，建立了罐内油量与油位高度和变位参数一般关系的模型二。再利用题给附表2中的数据来计算变位参数，计算每次进油后，油位高度变化所对应的罐内油量的变化，我们利用网格搜索法求解实际出油量与计算得到的出油量的差值平方和最小时的变位参数，逐步提高精度，让变位参数值精确到小数点后两位,用一次性进油前的数据来计算变位参数得到

1.96, 5.08αβ??==，用一次性进油后的数据来计算变位参数得到 1.96, 5.04αβ??

==，对这两个变位参数分别求取平均值得到最终的变位参数为 1.96, 5.06αβ??== 。然后将求得的变位参数和油位高度代入模型二，得到高度间隔为10cm 的罐容表标定值。接着我们在储油罐无变位的情况下，将题给附表2中的显示油高带入模型二，求解得到的罐内油量与题给的显示油量完全吻合，由此我们可以验证模型的正确性。再将

1.96, 5.08αβ??

==代入模型二中，

通过题给附表二中一次性进油后每次出油显示油高的不同得到罐内油量的变化，用它与一次性进油后的实际出油量做相对误差，得到所有相对误差的平均值为0.58%，用同样的方法，将 1.96, 5.04αβ??==代入模型二中，得到所有相对误差的平均值为0.52%，它们的值都小于1%，误差很小，从而可以验证我们的方法是可靠的。

最后我们对模型的优缺点进行了评价，提出了求解棱台体积来简化储油罐体积的积分运算量和用非线性回归的方法求解变位参数的方法，并对模型进行了简单的推广。

关键词：罐容表标定；误差修正；变位识别；网格搜索法；

一．问题的提出与重述

1.1问题的提出

加油站的地下储油罐有一个与之配套的“油位计量管理系统”，通过预先标定罐容表可以得到罐内储油量与油位高度之间的关系。储油罐在长时间的使用中，其位置会因为地基变形等原因发生改变。如何识别储油罐的变位情况及对罐容表进行重新标定是加油站关心的问题。

1.2问题的重述

题给两种储油罐：第一种为典型油罐，其主体为圆柱体，两端为球冠体；第二种为小椭圆形储油罐，其主体为椭圆柱体，两端为平头。在长期的使用过程中，储油罐可能发生横向偏转倾斜和纵向偏转倾斜。横向偏转倾斜即油位探针所在直线偏离地平线垂直线一个角度β，纵向偏转倾斜即油罐底部所在直线偏离水平线一个角度α。为了获得罐内准确的油位高度和储油量情况，我们需要定期对罐容表进行标定，罐容表即为罐内油位高度与储油量的对应关系。

已知这两种罐体的尺寸，请解决如下关于储油罐的变位识别与罐容表标定的问题：（1）我们对第二种罐体分别对其无变位和纵向变位倾角为0

α=的两种情况做了

4.1

实验，得到了两种不同位置放置的罐体分别在进油和出油两种情况下的4组实验数据。请建立数学模型研究这种罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。

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

二．问题分析

对于问题一，我们通过三重积分寻找小椭圆形储油罐罐内油量与测量油位高度和纵向倾斜角度的理论关系，代入间隔为1cm的不同的高度，通过此关系计算出罐容表理论标定值，接着我们通过题给4组进出油数据对罐容表理论标定值进行修正，其中对进出油量和罐内容量的理论值与实际值的差值进行修正，从而可以得到贴近实际情况的无变位罐容表标定值与倾斜变位罐容表标定值。通过分析罐内油量一定值时，寻找倾斜变位后测得的油位高度与无变位时测得的油位高度的关系，从而我们可以分析得到罐体变位后对罐容表的影响，接着我们根据罐体变位对罐容表的影响对装置进行了改进。

对于问题二，我们寻找球冠形储油罐罐内油量与油位高度和变位参数的一般关系，通过题给附表2中一次性进油前的进油量和油位高度数据计算变位参数。然后通过求解出的变位参数与一般关系的表达式得到油位高度间隔为10cm的罐容表标定值。接着用题给附表2中显示罐内油量与油位高度在无变位的情况下对模型进行检验，并通过一次性进油后的数据对方法的可靠性进行检验。

三．模型假设

（1）每次实验前，进油管与出油管中均无油量残余；

（2）忽略储油罐在制造过程中产生的缝隙、突出物等对容积造成的误差；

（3）忽略注油管、出油管在罐体内所占体积；

（4）假设每次测量在相同环境下进行；

四．变量与符号说明

h ：表示油浮子测得的油位高度（即油浮子与油罐底部之间的距离）

，单位米； 'h ：表示发生横向变位时，液面距油罐底部距离，单位米； 0

h ：表示无变位时测得的油位高度，单位米；

a ：表示小椭圆油罐横截面的长半轴长，单位米； b

：表示小椭圆油罐横截面的短半轴长，单位米；

c ：表示油位探针距左端的距离，单位米；

d

：表示油位探针距右端的距离，单位米；

S ：表示油罐横截面的面积，单位平方米； T ：表示小椭圆油罐的长度，单位米；

i V ：表示储油罐中第i 部分油量体积，单位立方米； 'V :表示第一次修正后的标定值，单位立方米； ''V :表示第二次修正后的标定值，单位立方米； σ：表示第一次修正值，单位立方米；

1ξ：表示无变位进油的第二次修正值，单位立方米；

2ξ：表示倾斜变位进油的第二次修正值，单位立方米；

其余变量均在后文中说明

五．模型建立与求解

5.1问题一的解答 5.1.1模型一的建立

只需将储油罐中油量的体积V 用油位探针测得的油高h 和油罐纵向倾斜的倾角α表示出来，即可得到储油罐发生纵向变位时罐中油量体积与油位高度间的关系，从而可以得到理论的罐容表标定值。

我们以油罐左侧面椭圆的圆心为坐标原点，平行于油罐底部的轴为z 轴，平行于椭圆长轴的轴为x 轴，平行于椭圆短轴的轴为y 轴建立坐标系。沿z 轴将储油罐分为若干个体积微元，则每一个微元的体积可看成椭圆面积的一部分与微元长度dz 的乘积，即：

dV Sdz = （1）

对于每一个体积微元，当储油罐发生倾角为α的纵向倾斜时，其对应椭圆的面积可表示为：

H b b

S dy --=

?

?

（2）

式中，H 表示dx 对应液面距离油罐底部的高度，可以通过油位探针测得的油位高度为h 计算得到，具体计算式为：

()tan H h c z α=+-? （3）

当(0,tan )h d α∈?时，储油罐中油量体积与测得油位高度之间的关系可用下式表示：

cot ()tan 0

2h c h c z b b

V dz dy αα?++-?--=??

?

?

（4）

当[tan ,2tan ]h d b c αα∈?-?时，储油罐中油量体积与测得油位高度之间的关系可用下

式表示：

()tan 00

2T

h c z b b

V dz dy α+-?-=???

?

（5）

当[tan ,2tan ]h d b c αα∈?-?时，储油罐中油量体积与测得油位高度之间的关系可用下式表示：

()tan 0

2T

h c z b m b

V dz dy abm

απ+-?--=?+??

?

（6）

其中(2)cot m c b h α=--?。

5.1.2计算理论标定值

对于小椭圆储油罐，将T =2.45，a =0.89，b =0.6，c =0.4，d =2.05，l =0代入（4）、（5）、（6）式，我们可以计算油位高度在0.01m 与1.2m 之间所对应的罐内油量体积的计算值。

我们得到的计算值时，有很少的几值会产生虚数，我们发现虚数的虚部系数是一个

根式内的值在y 变化时会产生负数，我们

对根号内的表达式进行修正，修正方法为对b 加上一个极小的正值，在本问题中修正值取1010-，这样会使得计算得到的理论罐内油量为实数，且使得原来的实数值在精度为410-时不产生变化。由此我们可以得到无变位理论标定值(如表1所示)。

对 4.1α?=时用同样的方法得到倾斜变位理论标定值(如表2所示)。

当罐内油量低于 1.6744L 时，随着油量的增加，油浮子始终未浮起，油位高度为0m 时表明储油量为1.6744L 以下，当储油量高于4110.1L 时，随着油量的增加，油浮子不再上浮，油位高度为1.2m 时表明储油量为4110.1L 以上，因此，我们不对这两个极限位置进行标定。

5.1.3误差分析

我们用题给附录1及表1、表2中的数据，得到了无变位及倾斜变位后的实际罐内油量与理论罐内油量的曲线，从图1中可以看出，实际罐内油量的曲线总是在理论罐内油量曲线的上方。

罐内油量(L)

油位高度（m m ）

进油时测量值与理论值的比较

图1 进油时测量值与理论值的比较

从上面四组数值比较中不难看出，罐容表的标定还要受到诸多外部因素如压力、温度、蠕变等的影响。

综上所述，我们考虑利用实验得到的数据对无变位与倾斜变位两种情况下的理论罐容表标定值分别作出修正，以消除上述众多因素带来的误差，使得罐容表标定的数据更接近罐内实际油量。 5.1.4对标定值的修正

(一）第一次修正----对进出油量实验值与理论值误差的修正

对于题给附表1中的四组实验数据，我们可以得到各组实验的每次实际进出油量。通过油位高度变化，可以得到理论的进出油量，从而可以得到实验值与理论值的误差，我们修正此种误差，就是我们要进行的第一次修正。 （1）对无变位情况的修正：

我们首先分析无变位进出油的两组数据，整理每次实际进出油量为50L 的数据，同时通过理论计算得到各数据油位高度差所对应的理论进出油量()V h （具体数值见附表1和附表2）。

分析数据得到每次理论进出油量稳定在一个值左右，我们计算所有理论进油量的均值分别51.7437L 、理论出油量的均值为51.7442L 。从而可以得到实际进出油50L ，理论进出油量比实际进出油量要多1.744L ，故我们得到了第一次修正的修正值的表达式：

1.744

50

V σ?=

（7）

式中，σ表示第一次修正的修正值。

通过理论罐容表中油位高度所对应的油量不同，对罐容表标定值进行修正。

V V σ'=+ （8） V '即为第一次修正后的标定值。

为了表明第一次修正后标定值、实验数据、理论标定值间的关系，我们作出了第一次修正后的比较图像（如图2）：

050010001500

20002500

300035004000

罐内容量(L)

油位高度(m m )

第一次修正后的标定值

图2 第一次修正的标定值与实验数据、理论标定值的关系图

（2）对倾斜变位情况的修正：

我们分析倾斜变位进出油的两组数据，整理每次实际进出油量为50L 的数据，同时通过理论计算得到各数据油位高度差所对应的理论进出油量（见附表3和附表4）：

分析这些数据，它们的值在50L 左右波动较大，且理论出油均值为50.2262L ，理论进油均值为50.1845L 。由于波动较大修正意义不大，且均值与50L 接近，所以无需进行此种修正。

（二）第二次修正----对罐内油量实验值与理论值误差的修正

对于题给的两组进油实验数据，我们可以得到每次进油后的罐内油量。通过油位高度，可以得到理论的罐内油量，从而可以得到实验值与理论值的误差，我们修正此种误差，就是我们要进行的第二次修正。 （1）对无变位情况的修正：

我们通过无变位进油的数据中的累计进油量和罐内油量初值得到实际罐内油量，再将无变位进油数据中的0α?=及油位高度h 带入（4）、（5）、（6）式，得到理论罐内油量，用（8）式对其进行第一次修正， 再将其与实际罐内油量比较，并做差值得到一组样本容量为78的修正值数据。

对这组数据求均值，得到1ξ=-2.62563m 。对第一次修正后的标定值再进行平移修正，修正后的标定值为：

1V V ξ'''=+ （9）

从而，我们可以分别得到无变位（表3）时储油罐修正后的标定。

25 596.09 65 2191 105 3678 26 630.52 66 2232.9 106 3705.4 27 665.44 67 2274.7 107 3732 28 700.81 68 2316.5 108 3757.7 29 736.64 69 2358.2 109 3782.5 30 772.88 70 2399.7 110 3806.3 31 809.53 71 2441.2 111 3829 32 846.57 72 2482.5 112 3850.6 33 883.97 73 2523.6 113 3871 34 921.73 74 2564.6 114 3890 35 959.83 75 2605.5 115 3907.6 36 998.27 76 2646.1 116 3923.6 37 1037 77 2686.6 117 3937.8 38 1076 78 2726.9 118 3949.8 39 1115.3 79 2766.9 119 3959.1 40 1154.9 80 2806.7 120 3964.2

为了更直观地表现修正后的效果，我们做出了对无变位的罐体最终的标定值与实验数据、理论标定值的关系图，如图3：

图3 旋转修正的标定值与实验数据、理论标定值的关系图

（2）对倾斜变位情况的修正：

我们通过倾斜变位进油数据中的累计进油量和罐内油量初值得到实际罐内油量，再将倾斜变位进油数据中的 4.1α?=及油位高度h 带入（4）、（5）、（6）式，得到理论罐内油量，将其与实际罐内油量比较，并做差值得到一组样本容量为53的修正值数据。

对这组数据求均值，得到平移修正值2ξ=0.07663m ，修正后标定值为：

2V V ξ'''=- （10）

（三）第三次修正----倾斜无变位标定值中异常值的处理

由于倾斜变位进油实验给出的数据较少，倾斜变位理论标定值经过第二次修正后得到的曲线在油量适中时能较好地接近我们实验得到的曲线。然而，在油位高度较小时，会标定出负数值的油量，在油位高度较大时，会与实际情况不符。

我们将油位高度为0.01m和1.19m时的倾斜变位理论值作为最终标定值，然后将第二次修正得到的0.07663m的修正值均匀的加到0.01m到0.35m和0.96m到1.19m 所对应的标定值。

从而我们可以得到倾斜变位时的最终标定值 (见表4)

为了更直观地表现修正后的效果，我们做出了对倾斜变位罐体最终标定值与实验数据、理论标定值的关系图，如图4：

05001000

150020002500

300035004000

200

400

600

800

1000

1200

罐内容量(L)

油位高度(m m )

倾斜变位罐容表标定值

图4 倾斜变位最终标定值与实验数据、理论标定值的关系图

5.1.5变位对罐容表的影响

当罐内油量为定值，油罐只发生倾斜角为α的纵向倾斜时（为顺时针倾斜），油位探针测得的高度与无变位时通过油位探针测得的高度有如下关系：

0tan h h l α=+? （11） 其中l 为油位探针到中心轴的垂直距离, 0h 表示无变位时测得的油位高度，h 表示倾斜变位时测得的油位高度。

随着纵向倾斜角α的增大，油位探针测得的油位高度与无变位时测得的油位高度间的差值会增大。因此变位的角度增大，相同高度下的倾斜变位的罐容表标定值会减少。以 4.1α?=为例，得到的无变位罐容表标定值与倾斜变位标定值见表3与表4，可以清楚的看到相同高度下罐容表标定值会减少。 5.1.6装置的改进

分析小椭圆油罐正面示意图，以及根据（11）式，当油位探针位于油罐中心轴时(即l =0时)，无论油罐的纵向倾斜角为多少，通过油浮子测得的油位高度可等同于油罐未

发生倾斜时测得的实际油位高度。

由以上的分析，如果我们将油位探针移动到油罐中心轴时，发生纵向倾斜时就不需要重新进行罐容表标定，以先前的罐容表得到的得到罐内容量是准确的。

5.2问题二的求解 5.2.1模型的建立

我们将油罐内储油量的体积分为圆柱体内油量和两个球罐体内油量三部分来求取。 在求取圆柱体油量体积的过程中，我们以储油罐右端为坐标原点，过圆柱体圆心平行于储油罐罐底的直线为z 轴，平行于圆柱体与冠体交线的直线为y 轴，垂直于的直线为zO y 平面的直线为x 轴建立坐标系。由于罐内油液面始终与水平线相平，则当储油罐发生倾斜角为α的纵向倾斜时，罐内不同位置的液面距底部距离可用测量液面高度来表示：

1'tan tan H h d c αα=-?+? （12） 其中'cos()h h β= （1） 当(0,tan )h d α∈?时，

储油罐内油量体积由左端冠体与中间圆柱体内油量的体积组成，左端冠体内油量的计算公式为：

cot 2

(1)

cot 1

2(1)h r

r

V dy R dx αα?+--=?-?

?

（13）

中部圆柱体内油量的体积为：

(2)

1

cot 0

2Z

H b d h b

V dz dy α

---=?

?

?

（14）

式中，tan tan H h d z αα=-?+?。 则储油罐内油量的总体积为

(1)

(2)

11

1

V V V =+ （15）

（2）当(tan ,(2)tan )h d r d αα∈??-?时：

圆柱体内油量与测量液面高度间的函数关系为：

18

(1)

2

00

2H r r

V dt dy --=??

（16）

右端冠体内油量体积：

2(2)

2

2(1))H r r

V dy R dx

--=-?

?

（17）

其中2223

(1)()2

R R -+= ， 2'tan H h d α=-? （R 为圆球冠体的所在球的半径）

左端冠体内油量体积：

3(3)

2

2(1))H r r

V dy R dx

--=-?

（18）

其中，3'tan H h c α=+? 则储油罐内油量的总体积为：

(1)

(2)

(3)

22

2

2

V V V V =++ （19）

当油罐中油量极少或超过一定高度时，我们需要重新考虑油位高度与储油体积间的对应关系。

（3）当(2tan ,2)h r d r α∈-?时：

储油罐内油量体积由左端冠体体积以及中间圆柱体、右端冠体内体积组成，左端冠体内

油量的计算公式为：

3

(1)

3

2(1))r r

V dy R dx

--=-?

（20）

中部圆柱体内油量的体积为：

6(2)cot (6)tan (2)

2

3

2[2(2)cot ]

b h h z r r

V dz dy r b h ααπα+---?--=+--?

?

（21）

右端冠体内油量的体积为：

tan

(3)

3

2(1)]h d r r

V dy R dx

α-?--=-?

?

（22）

则罐内油量总体积为：

(1)

(2)

(3)

33

3

3

V V V V =++ (23)

5.2.2数据的拟合

为了得到实验所用储油罐的纵向偏转角α与横向偏转角β，我们用最小二乘法对题中一次性进油之前的数据进行了拟合。

将题给附表2中每一个高度值i h 代入（8）式，可以得到理论储油量i V 。则每次输出油量的理论值为：

11,1,2,...,301i i i V V V i ++?=-= （24）

由题给附表2中数据，我们可以直接得到每次加入油量的实际值1'i V +?，

1,2,...,302i =。

我们将偏转角α与β的求解转化为最优化问题，即建立实际油量和理论油量差值的平方组成的目标函数，通过极小化该目标函数得到α与β的解，如式（24）。

301

2

1

11

(,)(')m in i i i f V

V αβ++==

?-?→∑ （25）

在用网格搜索法寻找最优αβ、参数的过程中，为了兼顾搜索速度和搜索精度，我们采用逐步细化的方法，保持搜索步长不变，减少每次搜索范围的长度，在搜索到较优αβ、的周围设置更为密集的穷举点，进行第二次甚至更多次搜索，直到达到所要求的精度位置。首先，我们对α取00.1为间隔，对β取01为间隔进行粗略搜索，得到了(,)f αβ的最小值大致在00(1.9,2.1)α∈与00(4,6)β∈ 之间，将搜索步长缩小10倍，得到了一组较优的αβ、组合。为了使精度更加准确，再次缩小搜索步长，可以得到(,)f αβ最小时对应的αβ、及其周围的点，如表5所示：

表5 三次搜索后最小(,)f αβ值及邻近点值

从表5中可以看出，min (,)246.19f αβ=， 1.96, 5.08αβ??==。 接着，我们用一次性进油后的数据，用同样的方法进行了三次搜索，得到了最小的(,)f αβ值及对应的参数值，最后搜索结果如表6所示：

表6 三次搜索后最小(,)f αβ值及邻近点值

由表6 1.96, 5.04αβ??==。对前面得到的得到的变位参数求均值，得到最终的变位参数：

1.96, 5.06αβ??

==

5.2.3罐容表的标定

对于实际储油罐，36,2,,82

d c r Z ===

=。代入题给数据，我们可以计算出油位高

度在0.1米与3.0米之间以0.1米为间隔的罐容表标定理论值，如表5所示：

（1）模型正确性的检验：

由于题给附表2中的油量容积是在通过无变位时得到的罐容表标定的，故我们将显示的油高带入（19）式中，得到理论的油量容积，与题给数据中的显示油量容积完全吻合，从而可以得知我们建立的模型二是正确的。 （2）方法的可靠性的检验：

我们通过附表2中未进油前的所有出油数据和显示高度(前302个数据)进行了数据拟合得到了变位参数，我们通过进一次油后的其它出油数据和显示高度对方法的可靠性进行研究。我们将显示高度与得到的变位参数带入（19）式中求得理论罐内油量值，再计算差值就得到了每次出油量1P ，将它与实际出油数据2P 进行比较可以计算一次实验的误差：

12

2

P P e P -=

（26）

我们用一次性进油前的数据计算变位参数，用一次性进油后的数据对得到的变位参数进行检验，得到所有误差的平均值：10.58%1%e =<

同样的，我们用一次性进油后的数据计算变位参数，用一次性进油前的数据对得到

的变位参数进行检验，得到所有误差的平均值

20.52%1%

e=<

从而可以检验我们得到的变位参数是很准确的，用网格搜索法求解最小出油差值平方和的方法是可靠的。

六．模型的评价与改进

6.1模型的评价

6.1.1优点

（1）我们对装置进行了改进，在只有纵向倾斜角度的情况下，如果将油位探针移至储油罐中心轴位置，可以消除纵向变位对罐容表的影响。

（2）通过精确的公式得到理论的罐容表标定值，并通过题给数据进行了三次修正，得到了较为准确，贴近实际的标定值。

（3）用两组数据求解了两次变位参数，并且对,αβ的值的准确性进行了检验；

6.1.2缺点：

由于在计算,αβ的时候我们只取了精度为0.01，由于计算方法的特殊性没有进行进一步的取值。

6.2模型的改进：

（1）我们可以寻找更好的积分方法得到储油罐的体积与油位高度的关系，从而简化计算量；

（2）计算变位参数时，取更高精度求解，并分析精度的提高对误差的影响；

（3）用求解棱台体积来简化储油罐体积的积分运算量；

（4）用非线性回归的方法求解变位参数的方法。

七．模型的应用与推广

我们的模型不仅适用于储油罐的变位识别和罐容表标定，还适用于模具制造及各种设计计算容量的场合。

参考文献：

[1]廉育婴，梁昭燕.立式罐底边界条件对静压力效应修正值的影响[J].现代计量测试，1996年第4期：27-29

[2]程继元.浅谈影响油罐标定与计量的因素及其修正方法[J].石油商技，第22卷第2期：33-35

附录

附录1：理论标定值计算的matlab程序：

（1）%纵向倾斜4.1度时的理论标定值

syms y z h

a=0.89

b=0.6

j=1;

for h=0.01:0.01:0.14

f1=a*sqrt(1-y^2/b^2);

f2=int(f1,'y',-b,h+(0.4-z)*tand(4.1)-b);

f3=int(f2,'z',0,(h*cotd(4.1)+0.4));

f4(j)=2*double(f3);

j=j+1;

end

j=15;

for h=0.15:0.01:1.17

f1=a*sqrt(1-y^2/b^2);

f2=int(f1,'y',-b,h+(0.4-z)*tand(4.1)-b);

f3=int(f2,'z',0,2.45);

f4(j)=2*double(f3)

j=j+1;

end

j=118;

for h=1.18:0.01:1.2

f1=a*sqrt(1-y^2/b^2);

f2=int(f1,'y',-b,h+(0.4-z)*tand(4.1)-b);

f3=int(f2,'z',0.4-(2*b-h)*cotd(4.1),2.45);

f4(j)=2*double(f3)+pi*a*b*(0.4-(2*b-h)*cotd(4.1)); j=j+1;

end（2）%无变位时的理论标定值

syms y z

a=0.89

b=0.6

j=1;

f1=a*sqrt(1-y^2/b^2);

for h=0.01:0.01:1.2

f2=int(f1,'y',-b,h-b);

f3=2*2.45*f2;

f4(j)=double(f3);

j=j+1;

end

附录2.修正值计算的matlab程序：

（1）%无变位修正及作图

t0=load('G:\无变位进油.txt');

x0=t0(:,1);

y0=t0(:,2);

syms y z

a=0.89

b=0.6

g1=a*sqrt(1-y^2/b^2);

for h=1:length(y0)

g2=int(g1,'y',-b,y0(h)/1000-b);

g3=2*2.45*g2;

g4(h)=double(g3);

end

for i=1:length(y0)

g5(i)=1000*(g4(i)-g4(i)*1.744/50);

wucha(i)=g5(i)-(x0(i)+262);

end%计算平移修正值

load('G:\无变位标准值.mat');

k=linspace(10,1200,120)

for i=1:120

f5(i)=f4(i)-f4(i)*1.744/50

end%旋转修正

hold on

%plot(1000*f5,k,'r-')

hold on

t0=load('G:\无变位进油.txt');

x0=t0(:,1);

y0=t0(:,2);

plot(x0+262,y0)

xuanzhuan=load('G:\无变位旋转修正值.txt');

k=linspace(10,1200,120)

xzz=xuanzhuan+sum(wucha)/78;%平移修正

plot(xzz,k,'g--')

axis([0 4000 0 1200])

hold on

plot(1000*f4,k,'k-')%画修正图（2）%倾斜变位修正值计算syms y z

a=0.89

b=0.6

chuyouliang=0.215;

f1=a*sqrt(1-y^2/b^2);

t0=load('D:\倾斜变位进油.txt')

x0=t0(:,1)/1000;

y0=t0(:,2)/1000;

for h=1:53

f2=int(f1,'y',-b,y0(h)+(0.4-z)*tand(4.1)-b); f3=int(f2,'z',0,2.45);

f4(h)=2*double(f3);

end

xiuzhengzhi=f4-x0'-chuyouliang;

junzhi=sum(xiuzhengzhi)/length(xiuzhengzhi);

附录3.理论值修正后的校验图像matlab程序：（1）load('D:\倾斜变位标准值.mat');

load('D:\变位进油高度理论值.mat');

plot(1000*(f4-0.0766),1000*k,'r-');

t0=load('D:\倾斜变位进油.txt')

x0=t0(:,1);

y0=t0(:,2);

hold on

plot(x0+215,y0,'b-');

（2）load('D:\无变位标准制.mat');

load('D:\变位进油高度理论值.mat');

plot(1000*(f4-0.0755),1000*k,'g-');

t0=load('D:\倾斜变位进油.txt')

x0=t0(:,1);

y0=t0(:,2);

hold on

plot(x0+262,y0,'k-');

附录4.对参数进行求解的程序

clear

syms x y z h

r=1.5;

f1=sqrt(r^2-y^2);

g1=sqrt((13/8)^2-x^2-y^2)-5/8;

g2=int(g1,'x',0,sqrt(r^2-y^2));

for w=2:0.1:3

for s=4.5:0.1:5.5

c=w*pi/180;

p=s*pi/180;

chazhi=0;

m=1;

f2=int(f1,'y',-r,h*cos(p)+(z-6)*tan(c)-r); f3=2*int(f2,'z',0,8);

g3=2*int(g2,'y',-r,h*cos(p)-6*tan(c)-r);

k3=2*int(g2,'y',-r,h*cos(p)+2*tan(c)-r);

kk=f3+g3+k3;