(11.17修改)基于数值积分的变位储油罐罐容表标定改进算法讲述.doc

基于数值积分的变位储油罐罐容表标定的改进算法※ 侯国亮1，李希瑞2，孙敏2，林泳坤2 （1.长春师范大学数学学院，吉林长春130032；2.长春师范大学工程学院，吉林长春130032） [摘要] 本文采用细化积分区域的方法给出了标定变位储油罐罐容表的数值积分新算法. 该算法克服了文献[1]中所提算法在油面高度较低时计算误差大的缺点. 数值实验表明，本文所提算法正确有效. [关键词] 细化积分区域；数值积分；罐容表；变位储油罐 中图分类号 O151；O24 一 问题背景介绍 通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之相配套的“油位计量管理系统”，采用流量计和油位计来测量进、出油量与罐内油位高度等数据，通过预先标定好的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况. 许多储油罐在使用一段时间后，由于地基变形等原因，就可能会使罐体的位置发生纵向倾斜和横向偏转等变化（以下称之为变位），从而导致罐容表发生改变. 按照有关规定，需要定期对罐容表进行重新标定. 图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体. 图2是其罐体发生纵向倾斜变位的示意图，图3是罐体发生横向偏转变位的截面示意图. 现阶段，变位储油罐罐容表的标定工作均是由人工或半人工操作实现的. 需要在加油站停工条件下进行，而且操作程序繁琐复杂，技术要求严格，时间周期长，这样做就既影响加油站的正常运转，也需要较高的经济成本，同时还会有一定的系统误差和测量误差. 为此，运用记录的储油罐进、出油量的实际数据，进行建模分析，正确识别罐体变位参数，实时修正罐容表，将是一项非常有社会现实意义的研究课题. 针对该问题，文献[1]运用坐标旋转变换、示性函数和Matlab中计算重积分的程序来计算任意变位后的罐内油量体积与油面高度的相互关系，大大简化了繁琐的数学推导过程(，是解决该问题的一种实用算法，值得推广. 但该方法的一个致命缺点是当变位储油罐内油面高度较低时计算的罐内油量与油面高度相对应的数据误差较大，因此如何克服此缺点就成了此方法是否能得到广泛推广应用的一个关键性问题. 本文在深入分析该法的计算程序基础之上，采用细化积分区域的方法提出了此法的改进算法程序，成功的克服了该致命缺点. 数值实验表明，新算法正确有效，为该类方法的实际推广应用奠定了扎实的基础. 二 算法的理论分析 1. 坐标系的建立 建立如图1、3(a)所示的固连于油罐、并随油罐一起转动的空间直角坐标系，以油罐中心为坐标原点，中央横截面为平面，轴、轴的正方向如图所示，轴正方向由右手定则确定，即沿中轴线方向. 规定，如无特别说明，下文中涉及到的区域及其约束条件的数学表达式均是在此坐标系下得到的. 2. 积分区域的确定 依据多重积分的几何意义可知，对罐内油量所占据的空间区域进行积分即可得到油的体积. 该问题中所涉及的积分区域受两类条件约束，一是油罐外壳的约束，二是油面所在平面的约束. 其中，油罐外壳约束可表述为 ，， 其中，表示油罐两端的球冠所在的球面半径，表示坐标原点到球冠所在的球面球心的距离. 下面讨论储油罐发生变位后油面所在平面的数学表达式. 为了叙述的方便，规定储油罐发生纵、横向变位时的两个转角按逆时针转动时为正，如图4所示. 根据相对运动的思想，当储油罐绕轴转动角、绕轴转动角时，我们可以设想成储油罐未动，而是罐内实际油面所在的平面以油浮子所在位置为定点，绕轴转动了角，绕轴转动了角，其中表示储油罐“油位计量管理系统”显示的油位高度. 这种等效转化的思想能在不影响体积计算的前提下，使得关于积分区域约束条件的讨论大幅简化. 接下来，本文就运用这种思想给出油面所在的平面以油浮子所在位置为定点，绕轴转动角、轴转动角后所得对应平面的数学表达式. 首先，由坐标变换相关知识可求出绕轴转动角，绕轴转动角的变换矩阵. 设变换前点的坐标为，变换后点的坐标为，则绕轴转动角的变换可表述为 ，， 绕轴转动角的变换为 ，. 一般来说，绕定点转动物体的有限角位移具有不可交换性，但由于该问题中所涉及的转动角，均很小，以致于变位次序对变位后的结果几乎没有影响，所以，旋转前的点与旋转后的点之间的关系可表示为 ， 或均可. 按展开得 ， 代入即可得由坐标平面绕轴旋转角，并再绕轴旋转角后所得对应平面的数学表达式为 . 其次，若把此时罐内油面所在的平面按照上述要求旋转后的平面记为平面，则平面与平面平行，区别在于平面经过油浮子所在的位置，也即油浮子所表示的点在平面上，所以平面的方程可表述为 ， 其中为一常数. 若将油浮子坐标代入，即可得 ， 所以按要求转动后油面所在的平面方程为 . 而油