毕业论文：储油罐的变位识别与罐容表标定大学生数学建模一等奖.doc

储油罐的变位识别与罐容表标定 摘要 本文对储油罐的变位识别和罐容表的标定问题进行了深入探讨，建立了储油量和油位高度以及变位参数之间关系的数学模型，主要应用了mtalab进行求解。 针对问题一，我们利用积分的方法推导出小椭圆储油罐在无变位和发生纵向倾斜变位时的一般公式。讨论了在储油罐发生纵向倾斜变位后对罐容表的影响，定义了平均影响率(变位前后储油量之差绝对值的平均值占总罐体容积的比例)作为评价罐体变位对罐容表的影响程度的大小的指标，求出。并分别给出了小椭圆储油罐在无变位和在纵向倾斜变位角取的罐容表。 表1 小椭圆储油罐罐容表(纵向变位) 油位高度 0.00 0.01 0.02 0.03 0.04 0.05 … 储油体积 00.0017 0.0035 0.0063 0.0100 0.0148 0.0207 … … 1.15 1.16 1.17 1.1713 1.18 1.19 1.20 … 3.9103 3.9338 3.9560 3.9588 3.9766 3.9955 4.10174.1101 针对问题二，将储油罐分为5个区域分别进行讨论，考虑到在球冠处的体积表达式过于复杂，我们省略了球冠处的一小部分体积，进行了近似求解，得出了罐内储油量与油位高度以及变位参数之间的一般关系的数学模型。 在利用储油罐的实际测量值估计变位参数时，我们建立了最小二乘拟合模型，得到了最佳的变位参数为：纵向倾斜变位，横向偏转变位。并据此对储油罐的罐容表进行了标定(见表3)。 在模型验证中，我们又采用蒙特卡洛模拟的方法对在问题二的模型中忽略的部分球冠体积进行了模拟计算。得到用问题二模型中求出的总储油量与模拟得出的总储油量一致度达到了99%，误差非常小，验证了我们所建立的模型的合理性和准确性。 关键词 平均影响率 最小二乘参数估计法 蒙特卡洛模拟 一 问题重述 通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。 问题一 为了掌握罐体变位后对罐容表的影响，利用小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为(=4.1的纵向变位两种情况做了实验。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 问题二 对于实际的储油罐，试建立罐体变位后罐内储油量与油位高度及变位参数（纵向倾斜角度(和横向偏转角度( ）之间的一般关系。利用罐体变位后在进/出油过程中的实际检测数据，根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。然后进一步用实际检测数据来分析检验模型的正确性与方法的可靠性。 二 模型假设 1、假设题中所给数据均为储油罐内壁测量值； 2、不考虑由于温度、压强变化等原因而引起储油罐的体积变化； 3、油位探针被固定在储油罐上，其上油浮子能够准确测量油位高度； 三 符号说明 椭圆的半长轴长 椭圆的半短轴长 储油罐的总长 油位探针到油罐底部左侧的距离 储油罐的总体积 储油罐的纵向倾斜角度 储油罐的横向偏转角度 油位高度 球冠体的半径 圆柱体的底面半径 注：未说明符号在文中用到时注明 四 问题一的解答 小椭圆储油罐罐体变位前后都可以应用积分的方法求出罐体的储油量和油位高度之间的关系。对于纵向倾斜的小椭圆储油罐，考虑分段求出其储油量和油位高度之间的关系，从而得到重新标定后的罐容表。 4.1 小椭圆储油罐无变位时的模型 由于此时的椭圆无变位，考虑先对二维椭圆进行积分。为方便表示油位高度，建立如图所示的坐标系，椭圆的半长轴长为，半短轴长为，则椭圆方程为 图1 对椭圆的积分示意图 在方向上取椭圆面中的一微元积分得到油的侧面积 储油罐内油的体积为 查积分表得到 (1) 利用matlab计算得到 经验证两种方法得到的体积公式完全等价，(1)式即为小椭圆储油罐无变位时的储油量和油位高度关系的模型。 根据此模型，我们可以求出小椭圆储油罐无变位时罐容表标定值(油位高度间隔取，结果见附录一)。 4.2 小椭圆储油罐纵向倾斜变位时的模型 储油罐纵向倾斜之后，油位计在油位过高或者过低时将不起作用(如图2所示的和区域)，考虑到倾斜角变化一般不会很大，所以我们可以将储油罐