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

1、卧式储油罐的变位识别与罐容 表标定卧式储油罐的变位识别与罐容表标定摘要本文分析卧式储油罐无变位和有变位时的罐内油位高度与储油量的对应关 系，建立了储油量与油位高度的关系模型。通过变位时罐内储油体积与无变位 时储油体积的变化差来分析罐体变位后对罐容表的影响。在问题一中：我们针对两端平头的小椭圆型储油罐，首先建立罐体无变位的储油体积关于油位高度的积分模型，通过对横截面积的微分ds，求出横截面积s ds，对截面积S在Y轴上积分，得出罐内储油体积模型： V y2sdy。y1 J 然后同样采用数学积分建立罐体变位后储油体积关于油位高度及纵向倾斜角 的关系模型：V yfs( )dy。在分析罐体变位对罐容表

2、的影响时，我们通过建 立罐体变位时容积v2关于变位角度的变化量v()的关系：v() 獎，来体d现罐体变位对罐容表的影响。然后在精确模型可以求得无变位时的容积 v(0)和变 位角度 的变化量v()的前提下，得到变位后容积的一般模型： V2( )v(0) v()。在计算罐体变位后罐内油面高度间隔为1cm时的罐容表时，我们分别采用精确模型和近似模型计算出理论标定值和近似标定值，同时采用分段二次拉格 朗日插值算法根据实际数据计算罐体变位后高度每隔0.01m时对应的标定值，通过比较和误差分析验证了所建立的近似模型的正确性与可行性，从而我们可 采用所建立的近似模型来给出高度间隔为 0.01m时的罐容表。在

3、问题二中：我们先考虑实际储油罐无变位时的罐内储油量与油位高度的 积分模型，将储油罐体积分成圆柱体 V1和两个球缺V2来计算，得到储油罐体积模型为：V(h) V1(h) 2V2(h)。在罐体有变位的储油体积计算时，我们先考虑横向偏转的影响，得到垂直 于油面的高度：h0 (h r) cos r，从而本问可同样近似采用问题一中所建立 的近似计算模型来计算变位后容积。根据附录所给数据，采用MATLAB编程求 解，得到纵向倾斜角和横向偏转角度的值为：3.45,5.45。最后，通过对近似模型计算得出的理论值与插值法根据实际数据计算出的数值进行相对 比较及误差分析检验模型的准确性与合理性，证实了模型是准确可

4、行的，并利 用近似模型给出了高度间隔为 0.1m的罐容表。关键词：回归分析插值法罐容表近似计算一、问题重述3油6m油 位岀油油浮注检 油查通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进 /出油量与罐内油位高度 等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行 实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生 纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按 照有关规定，需要定期对罐容表进行重新标定。图 1是一种典型的储油罐尺

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

6、进/出油过程中的实际检测数据（附件2）， 根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为 10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型 的正确性与方法的可靠性。二、模型假设与符号说明2.1模型的假设假设1:不考虑温度、压强对燃油体积的影响假设2:油位探针是始终垂直指向罐底的假设3:油罐体中燃油的进油和出油都是单独进行的假设4:罐体变位后对进/出油流量无影响假设5:罐体内附件体积对罐容量影响可忽略不计假设6:题中所给实验数据真实可靠油位探针ffl地平线油位探测装置2.2符号说明h由油位探针测出的油位高度ho有变位油罐体内的实际油咼y燃油覆盖所达

7、大最左端占的 Y 轴坐标值y i燃乂油覆盖所达最左端点 口 JT 轴坐标、值y2燃油覆盖所达最右端点的Y轴坐标值V储油罐内的油料体积a椭圆的长半轴长度b椭圆的短半轴长度s无变位储油罐内燃油没过横向截面的面积s平均误差rs平均相对误差ms最大误差mrs最大相对误差Zh无变位储油罐内油位高度为h时椭圆截面的浮油高度m小椭圆型油罐体内原点到左边椭圆之间的距离n小椭圆型油罐体内原点到右边椭圆之间的距离Sy有变位储油罐在丫轴坐标点为y时的横截面被燃油没过的面积zsh当油位高度为h时截面Sy上的燃油高度a储油罐的纵向倾斜角B储油罐的横向倾斜角V(0)罐体无变位时的储油体积V2罐体变位时的储油体积u表示ta

8、 ne表示cosR实际油罐体横切面所得圆的半径r实际油罐体中球缺所在球的半径三、问题分析本文要求用数学建模的方法来研究解决储油罐的变位识别与罐容表标定的 问题。针对问题一：此问要求我们利用一个小椭圆型储油罐来研究分析罐体变位后对罐容表的影响，并给出该罐体变位后油位高度间隔为1cm的罐容表标定值。为分析变位对罐容表的影响，我们希望通过探究有变位时与变位角度a相关的 体积变化量来研究变位对罐容表的影响。为此我们首先用积分法建立模型来求 解无变位时的储油体积和有变位时的储油体积，然后通过题中所给实验数据来 检验该体积求解模型，最后将有变位时的体积求解模型进行对a的微分得出与 变位角a相关的体积变化量

9、，此变化量即为罐体变位对罐容表的影响。为得出 该罐体变位后油位高度间隔为 1cm的罐容表，我们先用精确模型和近似模型进 行计算得出理论标定值和近似标定值，再运用插值算法根据附录中的数据进行 计算求解，以此来对比验证两类模型的准确性。针对问题二：此问首先要求我们建立一个罐体变位后标定罐容表的数学模 型，对于该模型的建立我们可沿用第一问中的体积变化量思想，用变位产生的 体积变化量与对应的无变位时储油体积之和来求解有变位时的储油体积。然后 针对该问用附件2中的部分实验数据求解确定模型的变位参数，并进一步利用 其余实际检测数据来检验模型的正确性。最后，我们用两种方法标定罐容表， 一种是采用本文所提出的

10、模型，另一种是以实际检测数据用插值法，比较了两 种方法给出的罐容表的一致性。四、模型的建立与求解问题一：研究变位对平头椭圆储油罐罐容表的影响，并给出变位后油位高度间 隔为1cm的罐容标定值4.1.1模型的建立为了分析两端平头的椭圆柱体罐体变位对罐容量的影响，我们先考虑水平放置时油罐容油量与油位高度之间的关系模型，然后考虑纵向倾斜时储油罐容量 与油位高度之间的关系模型。(1)罐体无变位的模型1为便于分析计算，我们以油位探针与罐底的交点为原点0,以油位探针为Z轴，罐底水平线为丫轴，以垂直于Y0Z平面且过0点的直线为X轴，建立三维 坐标系，如图1中(a)和(b)所示：油傥探計(a)在小椭圆油罐内建立

11、坐标系(b)横向椭圆截面的坐标示意图图2无变位小椭圆型罐体的坐标系示意图为建立数学模型，首先可设立如下变量：h:由油位探针测出的油位高度；y1 :燃油覆盖所达最左端点的丫轴坐标值；y2 :燃油覆盖所达最右端点的丫轴坐标值；V :储油罐内的燃油体积； a :椭圆的长半轴长度； b :椭圆的短半轴长度； S:无变位储油罐内燃油没过横向截面上椭圆的面积；:无变(1)Sdyzh :无变位储油罐内油位高度为h时椭圆截面的浮油高度。(位时h=)椭圆柱体横向截面椭圆方程为：x2鼻学2 1a2b2当燃油高度为h时罐内的燃油横截面的面积：s oh27/b2 (Z b)2dz(2)b对截面积S在丫轴上积分，被积表

12、达式：积分区间为：y1，y2得到体积的积分模型：v :sdy(4)将s表达式代入，即：V 护 A 2 ab2 (Z b)2dzdya(h b) ,b2 (h b)2 abarcsinU 岂 y2 ybb 221（2）罐体有变位的模型1与罐体无变位时一样建立坐标系，如图3中（&）和（b）所示建立三维坐标系：（a）在小椭圆内建立坐标系（b）横向椭圆截面的坐标示意图图3有变位小椭圆型罐体的坐标系示 意图由于此问中只分析了无变位和纵向变位两种情况，所以这里的变位模型只建立有纵向变位的模型。首先可设立如下变量：h:由油位探针测出的油位高度；:罐体的纵向倾斜角度；m原点到左边椭圆之间的距离；n原点到右边椭

13、圆之间的距离；y,:燃油覆盖所达最左端点的 丫轴坐标值；y2 :燃油覆盖所达最右端点的 丫轴坐标值；V :储油罐内的储油体积；a：椭圆的长半轴长度；b:椭圆的短半轴长度；Sy :储油罐在丫轴坐标点为y时的横截面被燃油没过的面积；zsh:当油位高度为h时截面Sy上的燃油高度。其中由平面几何知识可得zsh与各截面在丫轴上坐标的关系为：shh yta n与无变位时一样，我们对每个横截面被燃油没过的面积求解，有:Sy:b)2dZ(7)b对截面积Sy在Y轴上积分得到体积的计算模型：y2VSydy(8)yi y将Sy的求解公式代入体积模型可得：V 弋 Zsh2j7b2 (Z b)2dzdyb (9)弋 y

14、ta门 2Jb2 (Z b)2dzdy这是一个关于h的分段函数。我们分段考虑y1和y2的具体取值。我们知道油罐的纵向变位有两种情况，即向左端倾斜和向右端倾斜。由于 储油罐是倾斜放置的，因此当燃油油位超过油浮子所测最大范围后，油罐体仍 可进油，燃油的储量是不可测的，当油位高度显示为0时，油罐内不一定无油，燃油的储量同样是未知不可测的。若以油位接近油位探针底端为临界点，则罐 体向不同的方向倾斜会有不同的最大未知储油值。可知向左倾斜和向右倾斜对 油罐罐容标定值有一定的影响。下面我们分别讨论向左、向右倾斜的情况。罐体向左倾斜由图4可得只有当燃油表面介于L1L4之间时模型才可求解： 当罐体内油位垂直高度

15、低于直线L1，即hr时的图形(b)高度hr时的图形图19横向偏转后的高度关系图图中的油位探针高度h=AB实际高度ho = MN横向偏转角度为，横切面圆 的半径为r，由图易得：(a)当高度 h r 时：ho=OM ONOAcosOB(h r) cosr。(b)当高度h r时：hoON OMOBOA cosr (r h)cos 。因为在两种咼度下得出的h0表达式(hr) cosrr (r h)cos ，所以无论探针油位咼度为多少都有:ho (hr)cosr o(36)此时的ho为竖向截面内垂直高度，因此可先考虑横向偏转的影响，在求得ho的基础上转而只求竖向变位的影响。4.2.4有变位的储油体积计算

16、模型由问题一求得的的实验小椭圆型储油罐体积计算模型可知：变位时体积的变化量是一个与高度h和纵向倾斜角度有关的函数关系式，此时只要求得体积 关于的变化量，我们便可从罐体无变位时的容积出发求得罐体有变位时的容 积。由问题一中求得的变位后的容积 V2(u)关于U和h的一般模型为：2 2V2(u)v(0) (y2 yj(kih kzhu ksu kq)在本问题中，虽然条件有所改变，但v2(u)与U和h的相关关系仍然存在，因此我们同样可延用问题一中的模型来近似计算变位后的容积V2(u);不过此时用k1 k2 k3 k4和u作为参数来求得v2(u)与u和h的关系。令cos e,可以得到：V2(u)V(0)

17、 u(h r)er)2(k2ukJks(hr)e r)gu2k4)(37)4.2.5模型求解据附表二前半部分出油的数据，对式(34)进行回归分析，用 MATLAB程求 得各项系数的值为(程序见附录八)：k10.8239 k2-4.8307 k3-17.0932 k4=6.3798 u=0.0603 e 0.9955由 e cos 0.9955 得5.45。由 u tan 0.0603 得3.45。因此可建立V2(u)关于h的回归模型为：v2(u)v(0) 2.45(0.8239(h 1.5)0.9955 1.5)24.8307(h 1.5)0.9955 1.5)u 0.9955u2 6.379

18、8)(38)进一步可利用此模型，通过 MATLABS程(程序见附录八)，据附表二倾斜变 位的出油数据，可绘出相同高度下关于油体积的理论值和实际值图形如图20:65.55V 4.543.5315001600170018001900200021002200230024002500h/mm图20由图形可知：理论值与实际值相当吻合，从而说明了模型的可行性与准确性。426给出罐体变位后油位高度间隔为 10c m的罐容表标定值我们根据附录中给出的罐体变位后油面高度及油料容积值数据，先利用求得的近似模型计算出油位高度每隔 0.1m对应的油料容积，再应用三点拉格朗日插 值法计算出油位高度每隔0.1m对应的油料

19、容积，采用 MATLAB解，将计算出 的数据绘制成曲线，其结果如图 21和图22所示。(数据见附录二，程序见附录 九和附录十)3 m近似模型的预测有变位岀油55504540352h/m图21图22由式(18)，(19)，(20)，(21)可得:平均误差：s 21.58(L)平均相对误差：rs 0.0048最大误差:ms 22.35(L)最大相对误差：mrs 0.0066由图形和误差分析可知，插值算法与近似模型算得的结果近似一致，两者差值保持在一个很小的范围内，从而再一次验证了近似模型的准确性。4.2.7给出罐容表在已验证近似模型较为准确的情况下，我们利用近似模型来给出罐容表 通过MATLAB程

20、求得罐容表油高 /mm0100200300400500600700800900:容积/L0591016823102478366828768110131339515893油高/mm1000110012001300140015001600170018001900容积/L18488211632390126686295013233235163379794076443501五、模型的评价与推广5.1模型的评价在问题一中我们首先建立无变位时的体积求解模型，再建立有变位时的体积求解模型，而无变位时的体积求解模型即为有变位时a=0的特殊情况，模型建立具有一定的连贯性和综合性，并结合实验数据对模型进行检验验证了

21、我们 体积求解模型的准确性。我们直接用与变位角a相关的体积变化量表示变位对 罐容表的影响，对其影响的反映准确明了。但在模型建立时我们忽略了罐体内 附件体积及温度、压强等对储油体积的影响，而实际中这些因素都是对其有影 响的，因而会产生一定的误差。在问题二中，我们先建立求解实际储油罐无变位时储油体积，然后沿用问 题一中的模型来近似计算变位后的容积，采用拟合，得到变位参数与变位后体 积变化量之间的函数关系。对某一高度，用变位时体积变化量与对应无变位时 储油体积之和来表示有变位时总的储油体积。该近似算法与用积分计算的模型 相比，避免了繁冗复杂的计算，思路简单，可操作性强，并结合第二问中的实 验数据进行

22、参数确定和模型检验，检验结果很好地验证了我们模型的准确性。5.2模型的推广该模型是一个解决储油罐变位后油量体积计算的模型，适用于各种封头形 式的卧式容器无变位和变位后不同液面高度的体积计算。该模型同样适用于土 方等体积的计算。六、参考文献1高纯一，周勇，高等数学，复旦大学出版社，2006。2颜庆津，数值分析M，北京：北京航空航天大学出版社，2000年2004-8-8。3 王郑耀，卧式加油罐剩余油料体积的计算，西安交通大学,附录一（问题一中的罐容表）油高/mm容积/L油高/mm容积/L油高/mm容积/L油高/mm容积/L097.859310623.26201798.29303138.51089.

23、983320654.556301840.994031802086.653330686.486401883.79503221.13086.314340718.976501926.79603261.94088.401350752.016601969.99703302.45092.584360785.596702013.19803342.46098.636370819.676802056.59903382.170106.39380854.26690210010003421.480115.71390889.337002143.510103460.290126.5400924.877102187.210

24、203498.6100138.66410960.877202230.910303536.4110152.11420997.317302274.610403573.7120166.794301034.27402318.410503610.3130182.644401071.47502362.210603646.4140199.64501109.1760240610703681.8150217.634601147.27702449.810803716.5160236.684701185.67802493.610903750.4170256.724801224.47902537.411003783.

25、4180277.74901263.68002581.111103815.6190299.5950013038102624.811203846.8200322.365101342.88202668.411303876.9210345.975201382.98302711.911403905.9220370.415301423.38402755.311503933.5230395.6454014648502798.611603959.7240421.6455015058602841.711703984.2250448.385601546.28702884.711804006.6260475.845

26、701587.68802927.511904026.52705045801629.38902970.212004051.6280532.845901671.29003012.6290562.336001713.49103054.8300592.466101755.79203096.8附录二（问题二匚F的罐容:表）油高 /mm0100200300400500600700800900容积/L0591016823102478366828768110131339515893油高/mm1000110012001300140015001600170018001900容积/L1848821163239012

27、6686295013233235163379794076443501附录三：（图5程序）%无变位的罐容表的制定及其误差分析clearclch=0.159020.296030.413320.522950.629960.714910.814190.17614 0.19259 0.2085 0.22393 0.23897 0.25366 0.26804 0.28216 0.30969 0.32315 0.33644 0.34957 0.36256 0.37542 0.38816 0.40079 0.42576 0.43812 0.4504 0.46262 0.47478 0.48689 0.4989

28、5 0.51097 0.5349 0.54682 0.55872 0.57061 0.58248 0.59435 0.60622 0.61809 0.64185 0.65375 0.66567 0.67763 0.67854 0.69053 0.69082 0.70285 0.72703 0.73919 0.75142 0.7637 0.76416 0.77653 0.78899 0.80154 0.82695 0.83983 0.85284 0.8660.920451.06240.93461 0.94905 0.9638 0.97891 0.99443 1.0104 1.0271.04421.08161.10231.12531.15241.1935;v仁322.88374.63426.36478.13529.85581.61633.35685.08736.85788.58840