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

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