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

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

加油站的储油罐通常配备“油位计量管理系统”。其原理是用流量计测量进给油量和罐内油量之间的油量和高度，并实时计算预测量的罐容表（即罐内油量与储油量之适当关系），以保持罐内油量和储油量的变化。然而,许多储油罐在使用一段时间后,由于地基变形等原因,使罐体的位置发生纵向倾斜和横向偏转等变化(以下称为变位),从而导致罐容表发生改变。按照有关规定,需要定期对罐容表重新进行标定。针对该问题,作如下假设:储油罐变位时未发生形变;储油罐在进出油过程中形状保持不变;室外温度的变化对油体积变化的影响忽略不计;油罐壁光滑均匀、忽略壁厚。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.89m,b=0.6m,L=2.45m,联立式(2)、(3)、(6)、(7),代入相应的数值,可得V(H)与H之间的关系为其中0<H≤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)当-b<h0≤z1tanα-b时,罐容表的作用区域为图3中的区域Ⅰ。3)当z1tanα-b<h0<b-z2tanα时,罐容表的作用区域为图3中的区域Ⅱ。4)当b-z2tanα≤h0<b时,罐容表的作用范围为图3中的区域Ⅲ。5)当h0≥b时,H=2b,由于罐体纵向倾斜,储油罐不能继续进油,罐容表的测量值达到最大值,此时储油罐并未装满,罐容表出现第2个测量盲区(即图3中区域Ⅲ上方的测量盲区),罐中的储油量为综上所述,当储油罐罐体纵向倾斜后,储油罐内的油位高度H与罐内储油量V(H)之间的关系为其中H=h0+b,即V(H)反映了储油罐罐体变位后对罐容表的影响。1.3油位高度计算利用罐体纵向倾斜后油位高度标定罐容表的数学模型,即式(9),可以给出油位高度间隔为1cm的罐容表标定值(油位高度从1～120cm),如表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)进行误差分析,此时各组数据的相对误差η≤1.5%,说明此时模型效果较好。对表1得出的罐容表标定值进行修正,得到的修正值如表2所示。从修正结果来看,修正效果明显的是储油罐油位高度在25～113cm内的标定值。油罐出油时的模型检验及修正与油罐入油时的模型检验及修正方法一致,在此不再赘述。2实际储油位的识别与容器体积表的测定2.1油罐截面坐标系的建立当实际储油罐水平卧置时,参照问题1中同一背景下类似问题的思想与解法,可建立实际储油罐水平卧置时V(H)与H之间的关系模型,此时油罐的截面为圆面,建立如图4所示的坐标系。图4中,hz=4−4h20+10.5625√+438hz=4-4h02+10.5625+438。已知2R=3,L=8,c=1,其中c为球冠体的高,设Va(H)、Vc(H)为油罐左右两端球冠体的储油量,Vb(H)为圆柱体的储油量,V(H)为油罐储油量,油罐水平卧置时有:1rh00时20r时2.2罐体倾斜后的总储油量的计算罐体纵向倾斜变位后,以罐体上油浮子与油位探针所在直线为y轴,罐体的中心轴线为z轴,2条轴线的交点为原点O,过原点O垂直于罐体yOz平面的直线为x轴,建立如图5所示的坐标系。当储油罐处于纵向倾斜变位时,可将储油罐罐体划分为7个区域(如图5所示),其中每一个区域都是能反映罐容表对罐内储油量测量的作用。通过对式(3)及图5的分析可知,在坐标系中无论h0取何值,H与h0之间仍然满足关系式:因而,可用h0与V(H)关系来代替讨论H与V(H)之间的关系。现对倾斜后油罐上半部分的球冠体进行单独分析,球冠体在yOz平面的投影如图6所示。根据图形几何关系可求得球冠体的半径、球冠体方程。通过对图5进行进一步分析可知,若对图中所示的各个区域分别积分,则运算麻烦,且不易实施。故可以将倾斜后的罐体分为3个部分,即右侧球冠体、圆柱体、左侧球冠体。此时,对罐体的3个部分分别进行分析与讨论,分别求出罐体3个部分的体积为Va(H)、Vb(H)、Vc(H)。现对罐体3个部分进行分析,得:1)左侧球冠体①当-R<h0≤-(1+z2)tanα时,②当-(1+z2)tanα<h0≤R-z2tanα时,③当R-z2tanα<h0≤R时,2)中间部分(圆柱体)①当-R<h0≤z1tanα-R时,②当z1tanα-R<h0≤R-z2tanα时,③当R-z2tanα≤h0≤R时,3)右侧球冠体①当z1tanα-R<h0≤(1+z1)tanα时,②当(1+z1)tanα<h0≤R时,根据上述分析,可知罐体倾斜后的总储油量为现对图5中的7个区域逐一进行分析,结果如下:1)当h0≤-R时,罐容表无法测量出罐内的储油量,罐容表的工作区域为测量盲区Ⅰ;2)当-R<h0≤z1tanα-R时,V(H)=Vb(H)+Va(H),罐容表的工作区域为区域Ⅱ;3)当z1tanα-R<h0≤-(1+z2)tanα时,V(H)=Va(H)+Vb(H)+Vc(H),罐容表的工作区域为区域Ⅲ;4)当-(1+z2)tanα<h0≤(1+z1)tanα时,V(H)=Va(H)+Vb(H)+Vc(H),罐容表的工作区域为区域Ⅳ;5)当(1+z1)tanα<h0≤R-z2tanα时,V(H)=Va(H)+Vb(H)+Vc(H),罐容表的工作区域为区域Ⅴ;6)当R-z2tanα<h0<R时,V(H)=Va(H)+Vb(H)+Vc(H),罐容表的工作区域为区域Ⅵ;7)当h0≥R时,罐内储油量无法测量,罐容表的工作区域为测量盲区Ⅶ。2.3储油罐内油位与水平偏转角间的关系罐体横向偏转后,罐体截面示意图如图7所示。分析图7可得,储油罐水平偏转后油面坐标h0与油位计的示数h′以及水平偏转角之间的关系为又有H=h0+R,可得储油罐内油位高度H与水平偏转角β之间的关系为将式(10)代入本文2.2节中所建模型,即可得出H、V(H)与偏转角α、β之间的关系。利用H、V(H)与α、β之间的关系式,即可确定油罐在不同油位高度H时所对应的储油量V(H)。2.4储油罐变位参数由本文2.2节与2.3节的分析结果可知,H、α、β共同作用于V(H),综合罐体内油位高度H以及油罐纵向倾斜变位与横向偏转变位对油罐储油量的影响,可得油罐的储油量为V(H,α,β)。由于利用积分的方法求解变位参数α、β运算复杂且不易实施,故通过设定参数(α,β)的取值范围,借助遍历搜索算法,最终可确定出储油罐的变位参数α、β,具体作法如下:利用相关数据信息,当对储油罐进行加油或者放油操作时,相邻的2个操作流水号之间油罐的实际进油量或出油量为ΔVT(n),即油位高度变化ΔH时的储油量变化ΔV(H)。例如,当储油罐在流水号为201与202之间进行操作时,油罐内油量的变化为ΔVT(202)=149L,此时将2个操作流水号所对应的油位高度H201与H202代入V(H,α,β)的表达式中,不断对α、β赋值,就可以得到不同变位情况下储油罐的进出油量计算值ΔVC(202)。采用相同的方法不断对α、β赋值,分别求出不同流水号时的储油罐内油量变化的实际值ΔVT(n)与理论值ΔVC(n),当α、β取不同值时ΔVC(n)与ΔVT(n)之间总的相对误差最小,所对应的α、β即为油罐的变位参数,即目标函数为由于ΔVT(n)=VT(n)-VT(n-1),不同组数据叠加后将会前后抵消,因而使得目标函数中的大部分数据相互抵消,使结果的准确性大大降低。为提高数据利用率以及结果的准确性,选取参考数据中的前400组数据,并将其等分为2组。令ΔVT(n)=VT(m+200)-VT(m),1≤m≤200,则目标函数变为在0≤α≤4°,0≤β≤5°内,以步长0.1进行遍历搜索,即可求出α、β的近似最优解,最终确定储油罐的变位参数为α≈3.4°,β≈2.8°。2.5罐体横向倾斜与横向偏转后油位高度l将求解出的α、β值分别代入本文2.2节中7个区域的H与V(H)关系式,即可求得油罐在各个区域的储油量V(H,α,β),即V(H),进而可以确定罐体纵向倾斜与横向偏转后油位高度标定罐容表的数学模型,并可以给出油位高度间隔为10cm的罐容表标定值(油位高度10～300cm),如表3所示。2.6罐容表的标定与模型的建立对现有数据进行分析可知,H与V(H)的参考值绝大多数落在本文2.2节中所建模型7个区域中的Ⅲ、Ⅳ、Ⅴ区域。为保证所得罐容表标定值具有普遍适用性,故在Ⅲ、Ⅳ、Ⅴ区域中均匀选取4个不同的油位高度H,并利用油罐变位后罐容表的标定模型求出相应的储油量V(H),如表4所示。将罐容表的标定值与题中给出的油罐出入油量进行比较,并用逐差法进行误差分析,具体做法如下:将计算所得12个数据依次编号为(1)～(12)并进行分组,其中(1)～(6)为第1组,(7)～(12)为第2组。用第1组中的数据减去第2组中的数据求得不同ΔH下的ΔV(H),并与所给的