数学建模全国一等奖论文储油罐的变位识别与罐容表标定.doc

《数学建模全国一等奖论文储油罐的变位识别与罐容表标定.doc》由会员分享，可在线阅读，更多相关《数学建模全国一等奖论文储油罐的变位识别与罐容表标定.doc（32页珍藏版）》请在优知文库上搜索。

1、储油罐的变位识别与罐容表标定摘要本文对储油罐的变位识别和罐容表的标定问题进行了深入探讨，建立了储油量和油位高度以及变位参数之间关系的数学模型，主要应用了mtalab进行求解。针对问题一，我们利用积分的方法推导出小椭圆储油罐在无变位和发生纵向倾斜变位时的一般公式。讨论了在储油罐发生纵向倾斜变位后对罐容表的影响，定义了平均影响率(变位前后储油量之差绝对值的平均值占总罐体容积的比例)作为评价罐体变位对罐容表的影响程度的大小的指标，求出。并分别给出了小椭圆储油罐在无变位和在纵向倾斜变位角取的罐容表。表1 小椭圆储油罐罐容表(纵向变位)油位高度0.000.010.020.030.040.05储油体积00

2、.00170.00350.00630.01000.01480.02071.151.161.171.17131.181.191.203.91033.93383.95603.95883.97663.99554.10174.1101针对问题二，将储油罐分为5个区域分别进行讨论，考虑到在球冠处的体积表达式过于复杂，我们省略了球冠处的一小部分体积，进行了近似求解，得出了罐内储油量与油位高度以及变位参数之间的一般关系的数学模型。在利用储油罐的实际测量值估计变位参数时，我们建立了最小二乘拟合模型，得到了最佳的变位参数为：纵向倾斜变位，横向偏转变位。并据此对储油罐的罐容表进行了标定(见表3)。在模型验证中，我

3、们又采用蒙特卡洛模拟的方法对在问题二的模型中忽略的部分球冠体积进行了模拟计算。得到用问题二模型中求出的总储油量与模拟得出的总储油量一致度达到了99%，误差非常小，验证了我们所建立的模型的合理性和准确性。关键词 平均影响率 最小二乘参数估计法 蒙特卡洛模拟一 问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转

4、等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。问题一为了掌握罐体变位后对罐容表的影响，利用小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为a=4.1的纵向变位两种情况做了实验。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。问题二对于实际的储油罐，试建立罐体变位后罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。利用罐体变位后在进/出油过程中的实际检测数据，根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。然后进

5、一步用实际检测数据来分析检验模型的正确性与方法的可靠性。二 模型假设1、假设题中所给数据均为储油罐内壁测量值；2、不考虑由于温度、压强变化等原因而引起储油罐的体积变化；3、油位探针被固定在储油罐上，其上油浮子能够准确测量油位高度；三 符号说明 椭圆的半长轴长 椭圆的半短轴长 储油罐的总长 油位探针到油罐底部左侧的距离 储油罐的总体积 储油罐的纵向倾斜角度 储油罐的横向偏转角度 油位高度 球冠体的半径圆柱体的底面半径 注：未说明符号在文中用到时注明四 问题一的解答小椭圆储油罐罐体变位前后都可以应用积分的方法求出罐体的储油量和油位高度之间的关系。对于纵向倾斜的小椭圆储油罐，考虑分段求出其储油量和油

6、位高度之间的关系，从而得到重新标定后的罐容表。4.1 小椭圆储油罐无变位时的模型由于此时的椭圆无变位，考虑先对二维椭圆进行积分。为方便表示油位高度，建立如图所示的坐标系，椭圆的半长轴长为，半短轴长为，则椭圆方程为图1 对椭圆的积分示意图在方向上取椭圆面中的一微元积分得到油的侧面积储油罐内油的体积为查积分表得到 (1)利用matlab计算得到经验证两种方法得到的体积公式完全等价，(1)式即为小椭圆储油罐无变位时的储油量和油位高度关系的模型。根据此模型，我们可以求出小椭圆储油罐无变位时罐容表标定值(油位高度间隔取，结果见附录一)。4.2 小椭圆储油罐纵向倾斜变位时的模型储油罐纵向倾斜之后，油位计在

7、油位过高或者过低时将不起作用(如图2所示的和区域)，考虑到倾斜角变化一般不会很大，所以我们可以将储油罐按液面高低分成五个部分，来求其储油量和油位高度之间的关系。我们讨论的是小椭圆储油罐纵向倾斜变位为逆时针旋转，如图2。对于储油罐顺时针旋转变位(即为负值)时的情况与此非常类似，在此不再详细讨论。图2 储油罐分区示意图4.2.1 对区域的讨论 在区域，其油位低于油位探针的油浮子，所以油位计量系统中显示油位高度为零。当油位计刚开始有示数时，计算其储油体积。将区域放大得到图3图3 区域的放大图图中，从原点纸面向里为轴，利用三重积分可以得到其中为油位探针到储油罐左侧的距离 积分得到 (2)4.2.2 对

8、区域的讨论由区域很容易得到区域的储油量和油位高度的变化关系，直接给出结论：所以 (3)4.2.3 对区域的讨论图4 区域示意图在小椭圆储油罐无变位模型中我们已经求出了的计算公式，同区域中的积分原理可以计算出，我们就可以得到此时的油量体积为 (4)其中4.2.4 对区域的讨论由区域4和区域2的相似性，将(3)式中的换为，将换为，并用总体积减去即为区域4的储油体积和油位高度的变化关系。其中为小椭圆储油罐的总体积化简并积分可得 (5)4.2.5 对区域的讨论在此区域中油浮子到达油位探针顶点，无法进一步测量油位高度。无法测量的总体积为： (6)4.2.6 综合各区域的罐容表标定的数学模型综上所述，我们

9、得到了储油量和油位高度、纵向倾斜角之间的分段函数关系式：表2 分段函数关系区域油位高度储油量12345根据储油量和油位高度的分段函数关系我们可以得到罐体纵向倾斜变位()后油位高度间隔为的罐容表标定值(见附录一)。当时，各区域油位高度及体积变化范围为：表3 各区域油位高度和储油量变化范围区域油位高度储油量123454.3 罐体变位后对罐容表的影响为了能更加准确地刻画罐体的纵向倾斜变位对罐容表的影响，我们分别对罐体变位前后的理论值和测量值进行多方面的比较。4.3.1 罐体变位前理论值与测量值比较根据附件一中所给数据，我们计算出在附件中所给的油位高度下理论值和测量值，并画出其曲线。图5 罐体变位前的

10、曲线对比通过对比我们发现，对于任意，储油量的理论值和实际值始终成如(7)式的比例关系。 (7)4.3.2 罐体变位后理论值与测量值比较同样根据附件一中所给数据，我们计算出在附件中所给的油位高度下理论值和测量值，并画出其曲线，如图6。图6 罐体变位后的曲线对比从图6中可以看出测量值仍然始终小于理论值，进一步求得理论值与测量值之差的变化范围为0.0454,0.0910，测量值的相对误差范围为1.56%，5.18%。4.3.3 罐体纵向倾斜变位前后理论值比较图7 罐体变位前后的曲线对比图8 同一高度下储油量的理论值与测量值之差变化关系图9 储油量的测量值的相对误差随油位高度的变化关系由以上各图可以清

11、晰地看出纵向倾斜变位后，使得在同一个油位高度下，变位后比变位前的储油量减小。但是这样仍不够直观，我们需要找到一个指标来定量刻画罐体变位后对罐容表的影响。从图9中可以看出，当油位高度较小时(0.1附近)，变位后相对于变位前的相对误差几乎达到了60%以上，但是此时的储油量的差别并不大，鉴于此，我们定义平均影响率：来刻画罐体变位后对罐容表的影响。可以求出在纵向倾斜变位时，4.87%。五 问题二的解答如图8实际的储油罐示意图，其主体为圆柱体，两端为球冠体，在储油罐无变位时我们计算其各部分体积。图10 实际储油罐示意图圆柱体积计算公式 带入数据得到：主体圆柱体积为56.5487一端球缺体积计算公式带入数

12、据得到：两端球缺总体积为，则储油罐的总体积为64.66455.1考虑不发生变位时储油量和油位高度的关系圆柱内的油体积随油位高度变化关系： (8)一端球缺内油体积随油位高度变化关系： (9)其中为圆柱体底面半径为球缺对应的半径为球内小圆半径总的储油罐内的油量对油位高度的变化关系为：5.2只考虑横向偏转变位时的储油量和油位高度的关系只考虑储油罐横向偏转为时：对实际的油位高度没有影响，但此时的油位探针已经随储油罐发生偏转(如下图)图11 只考虑横向偏转示意图由油浮子测量得到的油位高度与实际油位高度的关系为： (10)所以对罐容表影响转换公式：5.3只考虑纵向倾斜变位时的储油量和油位高度的关系只考虑纵

13、向倾斜变位时，我们利用问题一中的思想，将储油罐分成5个区域，分别计算储油量和油位高度的关系(如图12所示)。图12 实际储油罐分区示意图5.3.1对区域1的讨论油位探针测得的油位高度始终为零，分成三部分来计算这部分的体积，如区域1的放大图(13)图13 实际储油罐分区示意图区域1的总体积应为：对圆柱体部分进行三重积分得到：由球缺部分的体积随油位高度的计算公式(9)得：由于部分的体积精确计算非常复杂，而且储油罐的纵向倾斜角度一般不会高于所我们考虑将这部分体积省略，进行近似计算。由于倾斜角较小，所以区域3占据了储油罐的绝大部分，而在区域3中这种近似计算的误差将由于左右两个球冠的省略体积一正一负而有

14、所减小，所以由此得到区域1的总体积公式： (11)5.3.2 对区域2的讨论在区域2中，油位高度的变化范围。各部分储油体积和油位高度的变化关系:, 总储油量和油位高度变化关系： (12)5.3.3 对区域3的讨论在区域3中，油位高度的变化范围。圆柱体部分的储油量：其中球缺部分的近似储油量：，总储油量和油位高度变化关系： (13)5.3.4 对区域4的讨论在区域4中，油位高度的范围。圆柱体部分的储油量：球缺部分的总体积：，总储油量和油位高度变化关系： (14)5.3.5 对区域5的讨论在区域4中，油位高度始终等于。各部分体积分别为：，总储油量和油位高度变化关系：5.4综合考虑储油罐纵向倾斜和横向