储油罐变位识别罐容表标定08538.docVIP

储油罐的变位识别与罐容表标定 薛申芳 （邢台学院数学系， 邢台：054001） 摘要储油罐的变位识别与罐容表标定用四次多项式进行回归，旋转角为；最后利用所给出的模型和所求出的倾斜角和旋转角给出灌容标定值。 关键词 数学建模；罐容表标定；坐标旋转变换；积分1 问题 通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。现用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）对小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为(=4.10的纵向变位两种情况做了实验(实验数据参看CMCM2010A附件1),去建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 （2）对中间部分为柱面，两端为球冠面的实际储油罐，建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度(和横向偏转角度( ）之间的一般关系、利用罐体变位后在进/出油过程中的实际检测数据（参看CMCM2010A附件2）确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。 2 问题1 2.1 问题1分析 从附表1中，可以得到油位以及油体的改变量，从而可以建立微分方程模型。图1和图2分别给出了无变位时进油和出油情况的罐内油体增加量（出油时罐内油体增加量取负）与油位增量之比随油位变的离散图形。 在图1中有几个不正常的点，可能由于测量误差造成的；而从图2中的点比较正常，所以在无变位情况建模时采用附表1中的出油数据。 对倾斜情况，通过对附件1中的数据分析，采用了附件1中进油测量数据，且在计算变化率时把一组点：累计进油量为0.0017l（该值由积分计算得到）时，油位高度为0。 上述两种情况，通过对灌内油体的离散变化率进行拟合，得到连续变化率，进而可以求解连续变量的微分方程去得到灌内油体体积随高度的变化规律。 2.2 问题1变量说明 记 ：罐内油体的体积（）； : 罐内油位（）； ：油位增量（）； : 罐内油体的增加量（）； ：，即为油体关于的变化率，为1时表示无变位情况，为2时表示倾斜情况3 问题1数学模型 下面通过建立模型，得到无变位和倾斜两种情况下罐内油体的体积与油位高度的关系。 3.1 无变位情况 由2.1的分析以及2.2变量知，选用出油数据随着油位高度变化用四次多项式进行回归得： （） 这里显著性水平取为，决定系数为,,，可知回归方程非常显著，回归模型成立。微分方程初值问题： （） ，其解为 （） 上式即为无变位情况下罐内油体的体积与油位高度的关系。 3.2 倾斜情况 为了让结果更符合实际，利用附件1中倾斜进油变位数据，随着油位高度变化用四次多项式进行回归得： （4） 显著性水平取为，决定系数为,,回归模型成立。微分方程初值问题： （） ，其解为 （） 参看图3。当灌内油体相同时，而油标高度测量值差差异的部分数据参看表1。倾斜灌容标定值参看附 表1。 表1 相同油体对应的油位 V(m3) 0.0017 0.5017 1.0017 1.5017 2.0017 2.5017 3.0017 3.5017 h1( m) 0.0015 0.2190 0.3594 0.4843 0.6039 0.7235 0.8479 0.9870 h2 (m) 0 0.1854 0.3270 0.4582 0.5822 0.7 0.8144 0.9333 注：上表中为无变位油位，为倾斜油位。 4问题2中间部分为柱面，两端为球冠面的实际储油罐变位后标定罐容表的数学模型问题分析油罐变位和。 4.2问题：油罐纵向倾斜角(弧度),假设不超过10o ×π/180 )，且假定油罐右端上斜； ：油罐旋转角（弧度，假设不超过10o ×π/180 ）； ：油罐柱面半径（1.5）； ：油罐柱体的长度（8）； ：油罐油面铅垂坐标()； ：变位后罐内油浮子的高度,即浸泡在罐油中的浮杆长度（m）； V：罐油体积()； ：球冠半径(1.625)。