储油罐的变位识别与罐容表-同济大学数学系.docVIP

储油罐的变位识别与罐容表标定 摘要 加油站的地下储油罐中油量的储量和损耗都直接关系到加油站的经济利益，所以必须的对储油罐的储油容量进行精确的标定，从而使加油站进行完善的进销存控制。但目前加油站的储油罐多为一次性埋放并长期使用的，油罐会因地基变形等各种原因而发生变位，从而导致罐容表发生改变，因此需要定期对罐容表进行重新标定。 对于罐容表的重新精确标定，即研究解决储油罐的变位识别与罐容表标定的问题。对于这样一个问题，我们采用高等数学中的定积分知识，加上对于数据的理解，并用MATLAB编程求解，从而求出罐体变位后对罐容表的影响。 对于简化的小椭圆型储油罐，我们用MATLAB求出罐体变位前后油位高度与理论储油量之间的函数曲线，并结合附件中的实验数据做出实际函数曲线。将无变位进出油和倾斜变位进出油时我们的计算值和实际值进行对比，并进行了误差分析。在此基础上利用曲线拟合对模型进行了改进，并给出了罐体变位后油位高度间隔为１ｃｍ的罐容表标定值。 对于实际储油罐，我们再次采用积分的方法，将实际储油罐分为三部分即右侧球冠体、中间圆柱体和左侧球冠体分别进行积分，得到其体积与油位高度、纵向倾斜角度α和横向偏转角度β的积分表达式。在此基础上，我们适当简化了附件2中所给的数据，并用MATLAB画出了罐体储油量对油位高度的导数与油位高度的关系曲线。针对所得曲线的图形特征，我们用二次函数进行了拟合，从而得到其函数表达式。这样做可以使数据得到了充分和有效的利用，又避免了处理数据时的庞大计算量；并进一步利用待定系数法多次枚举α和β的值，积分求出其上述导数的理论值，与其拟合函数的实际值相比较，从而不断缩小纵向倾斜角度α和横向偏转角度β的范围，使两者的值尽可能的接近，最终确定变位参数：α=2.1°，β=4.3°；并给出了罐体变位后油位间隔高度为10cm的罐容表标定值，并直接带进附件数据，求其储油量的理论值与实际值的误差的标准差，最终得到的标准差仅为6.0068L，误差不到1‰。 关键字： MATLAB，曲线拟合，枚举法，标定，储油罐 问题重述 作为加油站最重要的产品——燃料油，其储量、损耗将直接关系到加油站的经济利益，而储油罐容积表的精确与否又直接决定了油站是否可以进行完善的进销存控制。通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表进行实时计算，以得到罐内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。 原题中的附件1为利用小椭圆型储油罐（两端平头的椭圆柱体）分别对罐体无变位和倾斜角为(=4.1°的纵向变位两种情况做实验而得的实验数据。根据这些数据，研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 对于实际储油罐（主体为圆柱体，两端为球冠体），建立罐内储油量与油位高度及变位参数（纵向倾斜角度(和横向偏转角度( ）之间的一般关系。根据附件2中罐体变位后在进/出油过程中的实际检测数据，利用所建模型确定变位参数，给出罐体变位后油位高度间隔为10cm的罐容表标定值，并进一步分析检验所建模型的正确性与方法的可靠性。 问题分析 小椭圆型储油罐在纵向变位后两端的截面均为椭圆圆弧与平行于该椭圆长轴的直线所围成的特殊弓形，小椭球型储油罐可以看成由无数个弓形拼接而成。每个弓形的高度在沿着油面的方向上成线性分布，根据椭圆弧与平行于其长轴的直线所形成的弓形面积与，以椭圆中心点为圆心，椭圆短半轴为半径的圆弧与相应直线所围成的弓形面积成比例的关系，可以推算出截面椭圆的面积，然后再沿着油面所在方向积分算出储油罐的储油容积，从而确定储油罐变位对罐容表的影响。 对于椭圆，有；而对于圆，有。所以有，则，即如图椭圆下部分阴影部分的面积可用圆的同高度的弓形面积的倍代替。 图2 ：椭圆油面横截面示意图 根据附表1提供的数据，我们发现在进出油油位大致相等时，无变位和倾斜变位时储油罐的油量初值、累加进油量和累加出油量之和为一常数，并且总小于储油罐的理论总容积4108.06升。例如，在无变位的情况下，罐内油量初值为262L，进油使油位高度为450.40mm累加进油量为1100升，出油使油位高度为451.43mm时，累加出油量为2602.72升，所得和为3964.72升；进油使油位高度为678.54mm和出油使油位高度为678.63mm时的数据之和亦为3964.72升。在倾斜变位的情况下，同理比较得出在进出油油位大致相等时三者之和为常数3512.73升。