储油罐的变位识别与罐容表标定.docxVIP

储油罐的变位识别与罐容表标定 1 储油罐变位识别与罐容表标定 油藏地下储水池通常配备“油位计量系统”。通过计量和油位计算，测量输入和输出油位的油位和容器中的油位高度，并实时计算预测值的容器表（即容器中的油位高度与容器中的相应关系），以获得容器中的油位和存储量的变化。因为很多储存油罐会使用一段时间，由于地基变形等原因，储存油罐的位置会发生纵向倾斜和横向位移等变化（以下简称变异点）。根据相关规定，需要定期重新测量罐容表。 文献提出了研究解决储油罐的变位识别与罐容表标定的问题: (1)为了掌握罐体变位后对罐容表的影响,利用小椭圆型储油罐(两端平头的椭圆柱体),分别对罐体无变位和倾斜角为α=4.1°的纵向变位两种情况做了实验.请建立数学模型研究罐体变位后对罐容表的影响,并给出罐体变位后油位高度间隔为1cm的罐容表标定值. (2)对于实际储油罐,试建立罐体变位后标定罐容表的数学模型,即罐内储油量与油位高度及变位参数(纵向倾斜角度α和横向偏转角β)之间的一般关系.请利用罐体变位后在进/出油过程中的实际检测数据,根据所建立的数学模型确定变位参数,并给出罐体变位后油位高度间隔为10cm的罐容表标定值.进一步利用实际检测数据来分析检验模型的正确性与方法的可靠性. 2 储油罐的划分假设 (1)假设储油罐壁的厚度很小,在测量油位高度H时无需考虑壁的影响; (2)假设外界的物质不会渗入到储油罐里面; (3)不考虑储油罐自身的腐蚀和变形对储油量产生的影响; (4)不考虑储存油量的体积随着季节的变化而变化,例如热胀冷缩; (5)不考虑出/进油管以及油位探测装置的体积. 3 模型的构建和求解 3.1 罐内储油量与罐内储油量的函数关系 在问题一中我们只考虑两端平头的小椭圆型储油罐,利用定积分法求出了罐体在无变位与纵向变位两种情况下分别关于油位高度及纵向倾斜角α的体积表达式,然后利用文献附件1中的无变位数据对模型进行了检验,利用纵向变位体积V2与倾斜角α、油位高度H1之间的函数关系,分析了α的变化对罐容表的影响,并求出了罐体纵向变位后油位高度间隔1cm的罐容表标定值. 设小椭圆型储油罐无变位时储油量体积为V1,纵向变位时储油量体积为V2,S(y)为在y点时截面的面积,H1为罐体内油位高度,H2为无变位时油位高度,L为小椭圆罐体的长度,l1为油位探测装置到油罐两底面中较近一面的距离. 第一步 只考虑罐体无变位情况下体积与油位高度的函数关系.椭圆截面面积如图1所示,设截面面积为S(H2),则 S(Η2)=2abS(H2)=2ab∫-b+Η2-b√b-y2dy∫?b+H2?bb?y2?????√dy =12πab+F(m),=12πab+F(m), 其中, m=Η2b-1?F(m)=ab(m√1-m2+arcsinm). 利用截面积已知,筒体体积的计算方法,得到无变位时储油罐的储油量V1与罐内油位高度H2的关系如下: V1=∫L0S(H2)dz=[12πab+F(m)]L.(1) 第二步 考虑罐体纵向变位时储油量.如图2所示,设罐体变位时任意一点z处的高度是h(z),有: h(z)=H1+l1tanα-ztanα. 由上式求得任意椭圆截面横截面积: S(z,Η1)=12πab+F(m), 其中, m1=Η1+l1tanα-ztanαb-1=p-qz, p=m+tanαbl1,m=Η1b-1,q=tanαb. 利用截面积已知时,立体体积的计算公式,有 V2(Η1,α)={∫l10S(z,Η1)dz,Η1=0;∫z0S(z,Η1)dz,0Η1≤(L-l1)tanα;∫L0S(z,Η1)dz,(L-l1)tanαΗ1≤2b-l1tanα;πabl+∫Ll1-1S(z,Η1)dz,2b-l1tanαΗ12b;πabl1+∫Ll1S(z,Η1)dz,Η1=2b.(2) 其中l=l1-(2b-H1)cotα. 把文献附件1中无变位进油情况、无变位出油情况、纵向变位进油情况以及纵向变位出油情况的油位高度分别代入对应的体积表达式,求出理论上的储油量,并将它们分别与所对应的实际出油进行误差分析,得到各自的平均误差(见表1),因此我们认为变位对罐容表存在一定影响.于是我们根据模型V2=f(H1,α),对罐容表进行了重新标定(见表2). 3.2 罐体横向偏转运动体积的非线性检验 为了掌握罐体变位后对罐容表的影响,需求出罐体在无任何变位情况下的体积公式V3,再求出罐内储油量与油位高度及变位参数(纵向倾斜角α和横向倾斜角β),利用非线性最小二乘法求解出两个参数α,β,并通过罐内储油量与油位高度一般关系式对罐容表重新标定,进一步利用文献中附件2的数据来检验模型的准确性和方法的可靠性. 第一步 求出罐体无变位时的体积.为了求出球冠储油罐在水平面上无偏转时的体积V3,建立如下模型: 设球冠的球面