《1储油罐的变位识别与罐容表标定王力工指导解析》由会员分享，可在线阅读，更多相关《1储油罐的变位识别与罐容表标定王力工指导解析（26页珍藏版）》请在金锄头文库上搜索。
1、储油罐的变位识别与罐容表标定摘要加油站都有地下储油罐，一般是通过查看罐容表得知储油罐的剩余油量。由于地基变形等原因，储油罐会发生变位，从而导致罐容表的改变。因此研究储油罐的变位识别与罐容表重新标定具有重要意义。问题一，首先，把求平底椭圆柱储油罐的储油量的问题转化为求体积积分的 问题，做三重积分时，利用了平行于椭圆柱母线的截面微元，并利用Maple软件， 求出储油量关于油位高度和倾斜角:的关系表达式；然后通过分析理论值和实验 值的相对误差，利用=0的实验数据对该表达式进行了误差系数补偿，得到了 储油量与罐内油位高度及倾斜角的关系的 V（h=）数学模型，利用:=4.1的实验 数据对补偿后的模型进行
2、了检验，平均相对误差由补偿前的5%变为不到2%；在分析变位后对罐容表的影响时取定不同油位高度，研究储油量关于倾斜角的变 化关系，得到了当油位高度较高、较低及在一定中间范围时不同的变化规律。然 后基于此模型得到了变位后油位高度间隔为 1 cm的罐容表标定值（见14、15页）；问题二，研究了主体为圆柱、两端为球冠的储油罐变位后的罐容问题。首先 将该问题转化为球冠和圆柱所含油的体积积分问题。圆柱所含油的体积可利用问题一中的模型求解（其中油位高度要经过一定转化）；球冠部分利用三重积分直接运算很难计算；而后我们通过分析球冠采用了近似积分算法，得到了储油量关于显示油位高度、两个变为参数（横向偏转角及纵向倾
3、斜角）的一般关系式；为 求-1，基于此一般关系式V（h,，J，建立了目标为理论出油量和实际出油量 之差的平方和最小的优化模型，利用附表二中的出油量的前半部分数据，并利用逐步减小区间的搜索算法， -:同时逐次以1、0.1、0.01为步长，用Matlab进 行了三次搜索，求得=2.13J =4.44；;然后利用附表二中显示油高和显示储油 量两组数据，与、1取为0时显示油高对应的理论出油量进行比较，得到的相 对误差的数量级为10*（见24页图）；并利用附表二中出油量后半部分数据，与 :=2.13/ =4.44时所得出的理论出油量进行比较，得到平均相对误差为 0.57%，从而检验了模型的正确性与方法的
4、可靠性。最后利用此模型得到了变位 后油位高度间隔为10cm的储容表标定值（见23页）。关键词：储油罐 油位高度 储油量 纵向倾斜角 横向偏转角#一、问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进 /出油量与罐内油位高度 等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系） 进行实 时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生 纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照 有关规定，需要定期对罐容表进行重新标定。图 1
5、是一种典型的储油罐尺寸及形 状示意图，其主体为圆柱体，两端为球冠体。图 2是其罐体纵向倾斜变位的示意 图，图3是罐体横向偏转变位的截面示意图。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两 端平头的椭圆柱体），分别对罐体无变位和倾斜角为：=4.10的纵向变位两种情况 做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的 影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角
6、度 和横向偏转角度一：）之间 的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2）,根 据你们所建立的数学模型确定变位参数， 并给出罐体变位后油位高度间隔为10cm 的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正 确性与方法的可靠性。二、基本假设1、注油管、出油管及浮杆所占的体积忽略不计2、问题一中所给的是油桶内径的尺寸（即忽略壁厚对结果的影响）3、问题二中储油罐在没有变位时的罐容表准确无误4、假设在所研究的时间段内纵向倾斜角度：和横向偏转角度保持不变三、符号说明符号符号说明h:测得油浮子高度a:纵向倾斜角p横向倾斜角c油浮子距离油罐左端距离d油浮子距离
7、油罐右端距离a椭圆截面半长轴b椭圆截面半短轴3:柱状罐体左端面油面咼度S2:柱状罐体右端面油面咼度L油罐柱状体长度n相对误差v储油量R球冠头所在球面的半径n:油罐柱状部分截面半径四、问题的分析4.1概论该题主要是解决解决储油罐的变位识别与罐容表标定两个问题，而核心问 题是求储油量关于显示油位高度、 储油罐两个方向的转角的关系表达式。 求此表 达式的思路是先求油所占的体积关于参量s（参量s为和油高相关的参量）及两角度的关系式，进而利用参量s和显示油高的关系式求得最终表达式；由于问题二 中的圆筒具有旋转对称性（旋转时所要积分的体积形状不改变），故问题一和问题二中储油量关于和油高相关的参量 s的表达
8、式的求法（实际为三重积分）在本质 上是一样的，而区别仅在于参量s和显示油高的的关系表达式。4.2问题一对于问题一，可通过三重积分和Maple软件先求解储油量关于参量s和倾斜 角的一般表达式，然后找到参量 s和显示油高h的关系，代入表达式可得到储油 量关于显示油高和倾斜角的表达式； 然后将理论值与实验值对比，参阅资料得出 误差修正方法，利用倾角为0时的数据进行修正，并利用倾角为4.1的数据对修 正后的表达式进行检验，得到了较理想的模型。然后基于此模型求出当倾斜角变 化时对罐容表的影响，并得到4.1：时的罐容表标定值。4.3问题二对于问题二，同问题一一样，也需找到储油量关于参量 s和两方向倾角的一
9、 般表达式，由于积分过程很复杂，有必要对此积分进行近似求解，得到一般表达 式后，再找到参量s和显示油高h的关系式，代入即可得到储油量关于参量 s和 两倾角-1的一般关系。然后利用附表二中出油量和显示油高，基于一般关系 式，确立优化目标为理论出油量和实验出油量的误差平方和最小，变量为两偏角的优化模型，利用逐步逼近的搜索算法，得到了两偏角的值。变位后的罐容表很 自然的可利用一般表达式得到。 最后利用附表二中的显示油高和显示油量 （实际 为无偏角时的油高关于油量的精确对应关系），对模型的准确性验证；利用实验 出油量和理论出油量的对比，对模型的可靠性进行验证。五、模型的建立与求解5.1问题一-两端平头
10、的椭圆柱体5.1.1数据预处理根据附件1,在无变位的情况下，我们根据进油过程和出油过程分别作出 储油量V随h变化的图，il呂亡顷5图1无变位进出油过程储油量与高度的关系对比由上图可以看出，在无变位的情况下，用进油数据求得的V进（h）变化曲线几乎与用出油数据求得的V出（h）变化曲线相同。此外，在有倾角的时候，同理作出下图。1UU0科闻甯叶出沖廿号占油星耳注蛊哥建莫系iio 血 EJO JDD an 90C 1MD 1100 沌吃咼度血2aJU-社融归扫匸程林舉油左高度关系410 SO EOO 70J 3B K 11R HCD 幻忙専誉伽L - _ :； 30印DOWJDO过00山臣 - & rL
11、 =3 r- & r- 5 3 3 ? 3 1 1T軒ft毎图2有倾角进出油过程储油量与高度的关系对比同样可以得到上述结论。这就意味着，进油过程和出油过程地位是相等的（基本是相同的），如果我们可以对所有的进油量和出油量求其对应的储油量V，即可增加样本容量，减小误差。因此，在后文的计算过程中，我们会将所有的进油量值和出油量值当做求 储油量的样本点。5.1.2倾斜角为：（，0）时计算储油量V的一般模型1）状态分析讨论先过哪个点：当（0）一定时，桶内油量不同，对应的储油量计算公式不同。根据储 油量的多少，可以将此问题分为 5个状态。如下图图3油罐液面5种状态图状态序号状态名称I油量特少n油量较少川油
12、量中等IV油量较多V油量特多表1另一方面，不同的倾斜角所对应的5种状态有所不同，液面直线可能先过。2 , 或者先过04点。临界条件如下图所示:_j I 二一已,_图4临界条件示意图根据上图，我们可以求出此临界状态下的倾斜角:二arctan丄2、26.1，即2.45只要小于这个角，当储油量增加时，液面直线必定先过 02点。而在此题中， 储油罐的倾角是由于地基变形等原因引起的， 角度不可能太大。在此问题中，我 们认为，不会大于临界值，即液面直线必定先过 02点。因此，图4所表示的状态即为桶内油量的所有状态讨论先求那种状态：通过分析可知，在油罐倾斜的时候，如果运用积分的方法，不同的油位高度 h，储油
13、量的求解公式也不同。故求出的储油量 V,的公式应当为关于h的分段函 数。2）储油量V的计算-油量特少和油量特多在油量特少的情况下，油面高度h一直保持为零。即无法通过h的变化，求 得储油量V。即V在此情况下无法计算。同理，对于油量特多的情况，V,无法计算。3）储油量V,的计算-油量较少根据已知条件，我们以O点为坐标原点，建立空间直角坐标系，示意如下：图5油量较少立体示意图在图5的基础上，我们设定相关参数，得下图图6油罐yoz面截面图在上图中，截取在区间z,z+dz内的体积元dV。然后将dV在xoz面的截面 图，如下:图7油罐xoz面的截面图参考图5后，体积元dV为长方体。若以xoy面为底面，其面
14、积为 图6中的各个参数，可以求得此矩形底面的长和宽。长：在图6中，有如下等量关系:tan :OC =e tan :ACOD 二 OC tan :二 h etan：据几何关系得,EF宽:底面面积:体积元体积:tan :tan:-! =2a 1(z - b)2 b2(z-b)2h etan：zb2tan。(z-b)2h etan：-z1b2tan。S底-21-dz宽=PQdV 二 s底面 dz = 2a . 1DE = OD - OE = h +ctano - zDE h ctan ：- -zS底，根据(1)(3)(5)(7)(8)积分区间：求储油量须从O到D点对z积分，积分区间为0,OD，其中 由(9)式得 OD = h ctan：。积分模型：储油量y(h)满足Mh) = 0 n eg2a、. 1 _ (z - b)2 h Ctan zeta n化简结果：用Maple求解(15)式，化简后得到b2tan：dz(9)a 2b s J-22V3b ( b S”(一 2 a rcskn - . )s 2b (S b) 2b s(10)6b t am Ib-J其中，s = s(h,
