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

《储油罐的变位识别与罐容表标定》由会员分享，可在线阅读，更多相关《储油罐的变位识别与罐容表标定（9页珍藏版）》请在金锄头文库上搜索。

1、储油罐的变位识别与罐容表标定薛申芳邢台学院数学系， 邢台：054001摘要该文就储油罐的变位识别与罐容表标定问题1，用四次多项式灌内油体离散变率进行回归得到连续变化率，建立了连续型微分方程模型，通过该微分方程初值问题，得到倾斜灌内油体体积随高度的变化规律，然后把无变位与倾斜情况的灌内油体随测量高度的变化引起的灌内油量差异，以及当灌内油体相同时，而油标高度测量值差差异来判断分析影响。该文就问题2的倾斜旋转油罐，通过坐标系旋转方法，建立了计算灌内油体与油位关系的精确积分模型；由于灌内油体不是油位、倾斜角以及旋转角的初等函数，利用所给测量数据和所建立的数学模型进行变位参数确实定时采用了最小二乘方法，

2、且对局部被积函数的二元二次多项式展开近似处理，以及对倾斜角和旋转角的离散化把倾斜角和旋转角的范围化为一些离散小区间搜索，计算结果显示倾斜角为，旋转角为；最后利用所给出的模型和所求出的倾斜角和旋转角给出灌容标定值。关键词数学建模；罐容表标定；多项式回归；坐标旋转变换；积分1 问题 通常加油站都有假设干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表即罐内油位高度与储油量的对应关系进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜

3、和横向偏转等变化以下称为变位，从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。现用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 1对小椭圆型储油罐两端平头的椭圆柱体，分别对罐体无变位和倾斜角为a=4.10的纵向变位两种情况做了实验(实验数据参看CMCM2021A附件1),去建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。2对中间局部为柱面，两端为球冠面的实际储油罐，建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数纵向倾斜角度a和横向偏转角度b之间的一般关系、利用罐体变位后在进/出油过程中的实际检测

4、数据参看CMCM2021A附件2确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。2 问题12.1 问题1分析图2 出油变化率图1 进油变化率从附表1中，可以得到油位以及油体的改变量，从而可以建立微分方程模型。图1和图2分别给出了无变位时进油和出油情况的罐内油体增加量出油时罐内油体增加量取负与油位增量之比油体变化率随油位变的离散图形。在图1中有几个不正常的点，可能由于测量误差造成的；而从图2中的点比拟正常，所以在无变位情况建模时采用附表1中的出油数据。对倾斜情况，通过对附件1中的数据分析，采用了附件1中进油测量数据，且在计算变化率时把一组点：累计进油量为0.0017l该值由积

5、分计算得到时，油位高度为0。上述两种情况，通过对灌内油体的离散变化率进行拟合，得到连续变化率，进而可以求解连续变量的微分方程去得到灌内油体体积随高度的变化规律。2.2 问题1变量说明记：罐内油体的体积；: 罐内油位；：油位增量；: 罐内油体的增加量；：，即为油体关于的变化率的拟合函数；其中，为1时表示无变位情况，为2时表示倾斜情况。3问题1数学模型及其及解下面通过建立模型，得到无变位和倾斜两种情况下罐内油体的体积与油位高度的关系。3.1 无变位情况由2.1的分析以及2.2定义的变量知，现选用出油数据对随着油位高度变化的数据，用四次多项式进行回归得1：1这里显著性水平取为，其决定系数为,，可知回

6、归方程非常显著，回归模型成立。通过解微分方程初值问题： 2得到，其解为3上式即为无变位情况下罐内油体的体积与油位高度的关系。3.2 倾斜情况为了让结果更符合实际，利用附件1中倾斜进油变位数据，对随着油位高度变化的数据用四次多项式进行回归得： 4显著性水平也取为，决定系数为,，回归模型成立经验证比用出油数据拟合效果好。通过解微分方程初值问题： 5 得到，其解为 6无变位与倾斜情况的灌内油体随测量高度的变化以及灌内油量差参看图3。当灌内油体相同时，而油标高度测量值差差异的局部数据参看表1。倾斜灌容标定值参看附表1。图3 无变位与倾斜比拟油量差表1 相同油体对应的油位V(m3)0.00170.501

7、71.00171.50172.00172.50173.00173.5017h1( m)0.00150.21900.35940.48430.60390.72350.84790.9870h2 (m)00.18540.32700.45820.58220.70.81440.9333注：上表中为无变位油位，为倾斜油位。4问题2下面建立中间局部为柱面，两端为球冠面的实际储油罐变位后标定罐容表的数学模型。4.1问题2分析油罐变位包括倾斜和旋转两种变位方式，这里先对倾斜进行讨论，然后再加上旋转。首先建立适当的坐标系，对于倾斜变位，相当于实行坐标系旋转，在新的坐标系下，计算罐体中的油体与坐标高度的关系；对于旋转

8、变位，根据罐体的轴对称性，便可以得到油位的坐标高度与测量浮标高度的关系，从而得到灌内油体与测油浮标高度的关系的数学模型。然后利用最小二乘法得到参数和。4.2问题2变量记号记：油罐纵向倾斜角(弧度),假设不超过10o /180)，且假定油罐右端上斜；：油罐旋转角弧度，假设不超过10o /180；：油罐柱面半径1.5；：油罐柱体的长度8；：油罐油面铅垂坐标()；：变位后罐内油浮子的高度,即浸泡在罐油中的浮杆长度m；V：罐油体积()；：球冠半径(1.625)。5问题2模型现在建立问题2的数学模型。5.1无变位油罐的坐标系对无变位油罐而言，以柱面左侧圆的圆心为坐标中心，指向正上方，轴穿过柱面中心轴指向

9、右方，建立右手坐标系参看图4.图4 罐体无变位坐标系 记分别为平面与罐体柱面所交矩形的四个顶点参看图4。那么在坐标系下，油罐柱面方程为：(7)左、右球冠面满足的方程分别为： (8)和 (9)柱面左、右侧所在的平面方程分别为： (10)和 (11)5.2 倾斜油罐的坐标系假设油罐右端向上倾斜，倾斜角度为，这等价于将坐标系绕轴顺时针旋转角度，假设得到的新坐标系为,那么这两个坐标系之关系为2：或(12)其中为旋转矩阵正交矩阵，且 (13)这里旋转矩阵中的取旋转角度的负值绕轴逆时针旋转取正，顺时旋转针取负。于是在坐标系下，油罐柱面方程为： (14)左、右球冠面满足的方程分别为： (15)和 (16)柱

10、面左、右侧所在的平面方程分别为： (17)和 (18)记分别为点在坐标系下的坐标，那么 (19)且在假定最大倾斜角度不会超过10的情况下，易得： (20)5.3倾斜油罐灌容模型下面在坐标系下讨论倾斜油罐内油体与坐标高度和倾斜角的关系： (21)利用微元法知，在方向取微元，用过轴上的点作垂直于轴的平面去截取油罐体，得到的截面面积与有关，记之为,那么微元体的体积为，那么罐油体积可用积分 (22)得到。由于罐面是的分片函数，故对的不同区间范围进行分别积分，即分为三个区间范围分别进行积分3，即 (23)其中为当在不同区间上用过轴上的点作垂直于轴的平面去截取油罐体，得到的截面面积。分别表示用过轴上的点作

11、垂直于轴的平面去截取油罐体左、右球冠局部所得到的截面面积。它们的表达式分别为： (24)(25) (26) (25) (26) 其中将和 、的表达式代入(23)式便得倾斜油罐内油体与坐标高度和倾斜角的关系的数学模型。5.4倾斜旋转油罐灌容模型下面考虑加上旋转发生旋转角度为，参看图9，为罐体倾斜和旋转后的斜，且记：油罐最低点到油面的高度参看图5；：油罐只发生倾斜时，游标测得的高度参看图5和图6；图5 油罐倾斜后坐标图6 油罐旋转断面图由几何关系以及所建的坐标系，得：由于假定油罐右端往上倾斜，取负值，从而得到与的关系为 (27)将(27)式代入(23)式，即得到了罐内油体积与、及的关系的一般模型：

12、 (28)至此，就建立了倾斜旋转油罐灌容的精确数学模型。6 问题2变位参数确定及灌容标定1首先讨论、确实定。利用附件2中E列显示油高和F列显示油量容积的的数据，即由测量值及得： (29)为确定和，理论上讲课采用最小二乘法，令误差为，即通过求极值问题： (30)s.t. (31)去得出和。遇到问题是并不是的初等函数。故采用离散方法和MATLAB软件去解决。方法是把的区间-10/180,0以及的区间0, 10/180均分成一些离散个小区间每个小区间长度适当小，对每个去得到，去求 (32)由于 、的积分表示不是初等函数，故采用近似处理，即把 、积分表达式中的被积函数对、分别在0，-0.6、0，8.3

13、2进行二元二次多项式展开带来的体积误差，然后再进行积分求得 、的近似表示式，由于其解析式繁杂，用软件进行处理，不再写出相应的展开表达式。计算结果为参看程序附录2Vprog2_1。根据所得角度，把局部测量值相应的理论值进行了比拟，绝对误差参看图7。图7 结果比拟2罐容标定。利用上述求得的角度值对罐容进行标定，计算结果参看图8、计算程序参看附录2Vprog2_24。图8 标定结果图7总结问题1用四次多项式灌内油体离散变率进行回归得到连续变化率，建立了连续型微分方程模型，通过该微分方程初值问题，得到倾斜灌内油体体积随高度的变化规律，也可以采用问题2的方法去建立精确地积分模型得到倾斜灌内油体体积随高度

14、的变化规律。问题2的精确地积分模型中，变位油罐内油体体积随高度的变化规律用积分表示的，由于不是油位、倾斜角以及旋转角的初等函数，在变为参数确定是带来困难，计算机编程时只能近似处理，二元被积函数在何处进行多项式展开是导致误差大小的关键问题；为了克服灌内油体体积的积分表达式难以表示成油位、倾斜角以及旋转角的初等函数，采用了把倾斜角和旋转角的范围化为一些离散小区间，在区间各分点得到测量数据与模型数据误差平方和，然后搜索到使得该平方和最小的倾斜角以及旋转角，该方法一会有一定的误差，但倾斜角和旋转角的离散化越细密误差将会越小。参考文献1 茆诗松，程依明，濮晓龙 编。概率论与数理统计教程.北京：高等教育出版社，2004.72 刘林 编。航天器轨道理论。