基于微元法的储油罐变位识别与罐容表标定.docxVIP

基于微元法的储油罐变位识别与罐容表标定 加油站的储油罐通常配备“油位计量管理系统”。其原理是用流量计测量进给油量和罐内油量之间的油量和高度，并实时计算预测量的罐容表（即罐内油量与储油量之适当关系），以保持罐内油量和储油量的变化。然而, 许多储油罐在使用一段时间后, 由于地基变形等原因, 使罐体的位置发生纵向倾斜和横向偏转等变化 (以下称为变位) , 从而导致罐容表发生改变。按照有关规定, 需要定期对罐容表重新进行标定。针对该问题, 作如下假设:储油罐变位时未发生形变;储油罐在进出油过程中形状保持不变;室外温度的变化对油体积变化的影响忽略不计;油罐壁光滑均匀、忽略壁厚。 1 小椭圆形储油罐的变位识别与容器体积表的测定 1.1 储油罐内油位计算 小椭圆型储油罐 (两端平头的椭圆柱体) 的罐体可以被看成是一段卧置的椭圆柱体, 其正面、截面示意图如图1所示。 以小椭圆的长轴所在直线为x轴, 罐体水平卧置时油浮子所在直线为y轴, 罐体的中心轴线为z轴, 3条轴线的交点为原点O, 建立如图2所示坐标系。 设储油罐截面的椭圆长轴长为2a, 短轴长为2b, 罐内油位的高度为H, 此时相应的罐内储油量为V(H) , 则有截面的椭圆方程: 本问题归结为求储油罐内油品体积与油品液面高度的函数, 即: 其中:H=0,V(H) =0;H=2b,V(H)max=f(2b) 。 设储油罐内油位高为H时油面与椭圆截面的交点坐标为 (x,y) , 储油罐罐体长为L,Sxy表示罐内油位高为H时, 油面沿y轴方向纵截面的面积,h0表示油位高为H时, 油面在坐标系内的y轴坐标, 即: 根据微元法的思想及定积分的概念, 有面积元素 由式 (1) 可得 联立式 (1) ～ (5) 可解得 储油罐内油位高度为H时的储油量为 现可知,a=0.89 m,b=0.6 m,L=2.45 m, 联立式 (2) 、 (3) 、 (6) 、 (7) , 代入相应的数值, 可得V(H) 与H之间的关系为 其中0H≤1.2。 1.2 罐容表对于储油罐的作用 当倾斜角α存在时, 储油罐发生纵向变位, 此时可建立如图3所示坐标系。 图3中的z1, -z2分别表示油罐两端在z轴上投影点的z坐标, z0表示储油罐倾斜后任意纵截面与z轴交点的z坐标, 此时该截面上的油位高度为hy。 采用与建立储油罐无变位时模型相同的思想, 即可求出当储油罐纵向倾斜时, 储油罐纵截面的面积表达式Sxy以及罐内储油量表达式V (H) 与油位计测量出的罐内油位的高度H之间的关系。 当储油罐处于纵向倾斜变位时, 可将储油罐罐体划分为5个区域 (如图3所示) , 其中每一个区域都是能反映罐容表对罐内储油量测量的作用区域。 由式 (3) 以及对图3的综合分析可知, 在坐标系中无论h0取何值, H与h0之间的关系始终为 因而, 可用h0与V (H) 关系来代替讨论H与V (H) 之间的关系。 现对图3中的5个区域逐一进行分析: 1) 当h0≤-b时, 由于油位计不能测量出罐内油面的实际高度, 使罐容表不能测量出罐内的储油量, 出现罐容表的第1个测量盲区 (图3中区域Ⅰ下方的测量盲区) , 此时V (H) 与h0无关。 2) 当-bh0≤z1tanα-b时, 罐容表的作用区域为图3中的区域Ⅰ。 3) 当z1tanα-bh0b-z2tanα时, 罐容表的作用区域为图3中的区域Ⅱ。 4) 当b-z2tanα≤h0b时, 罐容表的作用范围为图3中的区域Ⅲ。 5) 当h0≥b时, H=2b, 由于罐体纵向倾斜, 储油罐不能继续进油, 罐容表的测量值达到最大值, 此时储油罐并未装满, 罐容表出现第2个测量盲区 (即图3中区域Ⅲ上方的测量盲区) , 罐中的储油量为 综上所述, 当储油罐罐体纵向倾斜后, 储油罐内的油位高度H与罐内储油量V (H) 之间的关系为 其中H=h0+b, 即V(H) 反映了储油罐罐体变位后对罐容表的影响。 1.3 油位高度计算 利用罐体纵向倾斜后油位高度标定罐容表的数学模型, 即式 (9) , 可以给出油位高度间隔为1 cm的罐容表标定值 (油位高度从1～120 cm ) , 如表1所示。 1.4 储油罐油位高度的误差分析 根据储油罐倾斜变位进油的相关数据, 对所建立的罐体纵向变位后罐容表标定模型进行效果检验及误差分析, 分别计算出罐内油位高度的标定值与实际测量值之间的绝对误差δ与相对误差η, 经计算可知45δ92, 1.5%η5.5%, 误差较大。对不同的油位高度H与绝对误差δ两项进行分析, 可以求得H与δ之间的关系为 利用式 (9) 、 (10) 可求出不同的H所对应的V(H) 与δ, 根据V′ (H) =V(H) -δ, 即可得到对罐容表标定出来的V(H) 修正后的修正值V′ (H) 。再对修正后的V′ (H)