《储油罐的变位识别与罐容表标定(数学建模论文)》由会员分享，可在线阅读，更多相关《储油罐的变位识别与罐容表标定(数学建模论文)（13页珍藏版）》请在金锄头文库上搜索。
1、1储油罐的变位识别与罐容表标定摘要加油站都有若干个储存燃油的地下储油罐，需要采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。但是许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化，从而导致罐容表发生改变。需要定期对罐容表进行重新标定。在求解过程中，我们对于罐体无变位、罐体产生纵向变位、罐体在水平和纵向都产生变位三种情况，利用解析几何的方式计算出体积与变位参数之间的关系，同时应用契比雪夫多项式对体积值进行近似多项式展开用以对标高和出油量的关系
2、进行拟合表示，得到较为满意的效果。第一问、(1)针对无变位情况，我们计算得到椭圆油罐容积表达式为：，利用契比雪夫多项式方法在提高拟合精ablvhvbhVrcsin)(122椭度的前提下用 5 阶多项式拟合处标高和容量之间的函数关系；(2)对于纵向变位的情况，当椭圆型罐体发生变位纵向变位角度 时,我们利用体积等效思1.4想，讲上述罐内不规则油量容积的计算转为(1)中规则油容进行计算，利用附件(1)中数据利用最小二乘拟合方法算出油位高度的真实值，继而利用拟合多项式：408.5976-H 395.752.-13.498H .6- 0.32H235 椭 变V进行间隔为 1cm 的此罐容表进行标定，得出
3、的表标定值如下：1cm 2cm 3cm 4cm 118cm 119cm 120cm0L 0L 0L 0L 4017.26L 4050.08L 4082.80L第二问、(1)利用第一问中等体积的思想，我们同样可以对纵向倾斜角度和横向倾斜角度 时进行数学模型的建立。(2)在模型的建立过程中得到一个关于浮游子高度 和偏转角 、 以及等效高度 之间的一个表达式，从而利Hh用最小二乘拟合确定变位参数 、 。(3)利用已给数据求得表达式： ，继而再次利用拟合拟合多项式得出间隔为 10cmtan2ta10hRz值：10cm 20cm 30cm 280cm 290cm 300cm1352.17 2223.85
4、 3082.13 60303.88 60690.98 61807.31利用附表（2）中的数据进而进行模型正确性与可靠性的检验。2关键词：储油罐 罐容表 变位 契比雪夫公式 等效高度 最小二乘拟合一、问题的重述加油站一般都有若干个储存燃油的地下储油罐，并且也都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定
5、，需要定期对罐容表进行重新标定。储油罐的主体为圆柱体，两端为球冠体。 用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为 的纵向变位两种情况o1.4做了实验，实验数据已给出。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为 1cm 的罐容表标定值。（2）对于实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度 和横向偏转角度 ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（已给出），根据所
6、建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为 10cm 的罐容表标定值。进一步利用已给出的实际检测数据来分析检验模型的正确性与方法的可靠性。二、问题的分析综述，该题要解决的问题是储油罐的变位识别和罐容表的标定。显然这是几何求解的应用问题。可以通过所学过知识求解。对于问题一，题目中给出一个小椭圆型的平头油罐的模型和纵向的倾斜角的值，并且也给出了对这个罐体无变位和有纵向变位两种情况的实验数据，所以我们可以把这个问题分成两部分来解决，首先解决当罐体无变位时罐中储油量与油液面高度的关系，其次解决当罐体发生倾斜时罐体中油的储油量随油液面高度的变化。然后对这两种情况求出的结果进行比较，判断出
7、油罐变位后对罐容表的影响。最后重新确定变位后罐容表的标定值。分析发现当油罐无变位时，问题比较好解决，我们可以利用积分的方法求出罐中的储油量和油液面高度的关系。当油罐变位时，首先要选好坐标系，确定倾斜后液面高度是怎样变化的。也是通过积分的方法进行求解，但是我们也可以寻找一点，可以把它等效成无变位的情况去解决。对于问题二，给的是一个实际储油罐，其主体为圆柱体，两端为球冠体。该问题主要是要解决罐内储油量与油位高度及变位参数（纵向倾斜角度 和横向偏转角度 ）之间的一般关系。对这个问题可以看成两部分求解。首先也是要选好坐标系，确定好油面高度和变为参数的关系，分别求解主体圆柱体和两端球冠体中的储油量和油位
8、高度的关系。建立出一定的数学模型，然后通过题目给出的数据，确定数据的可用性求出模型中的变位参数。最后通过误差分析检验模型的正确性和方法的可靠性。3三、问题的假设1、忽略油罐中温度，压强对储油量的影响；2、假设油罐的形状是规则的；3、假设储容表的读数是无读数上的技术误差；4、假设油罐变位的角度不会太大；5、假设罐中的油分布均匀；6、假设油罐不会因为出油的变化而发生位置变化；四、模型的建立与求解问题一：对于本问题，首先可以把它分解成两个问题来考虑，第一个问题是确定小椭圆型储油罐没有变位时储油量与油面高度的关系，利用解析几何的方法确定标高和油量之间的关系。第二个问题确定油体容量和标高及变位角之间的关
9、系，利用等效体积替换的思想把它转化为上一问中未变位的问题进行求解。同时我们分析发现变位数据的油高不是真实值，所以考虑利用已有读数寻求真值。依据上述思想分别对这两个小问题建立合理的模型。1 无变位时：小椭圆型储油罐的横截面很明显是一个标准的椭圆，首先建立合适的坐标系，并且椭圆的图形如图 1:图（1） 1、符号声明： a-椭圆罐横截椭圆面的半长轴 b-椭圆罐横截椭圆面的半长轴h-无变位椭圆油罐中的油位高度面 积无 变 位 时 椭 圆 的 横 切 面椭 S-小椭圆油罐的母线长lxx0h2ba4-方差和2R根 据 模 型 得 出 的 数 据iV实 测 数 据i2、模型建立：椭圆的方程为： 1)(22b
10、yax- - (1) 式中： a-椭圆罐横截椭圆面的半长轴且为 0.89 米b-椭圆罐横截椭圆面的半长轴且为 0.6 米油面高度为 h 时，储油罐横截面的面积可求得为：-dybaSh022椭（2）式中：h-无变位椭圆油罐中的油位高度；可推导得：-abhbhs rcsin122椭（3）我们可以设 ，由 可以得：bht01t可求得小椭圆油罐未变位时储油量随液面高度的变化为：-ablttVrcsin122椭（4） 式中：-小椭圆油罐的母线长且为 2.45 米；l在求解时我们可以利用最佳平方逼近多项式（即契比雪夫多项式, ） ， 求出储油量的近似值。根据公式（1）我们可以准确的计算出储油量，但是公式里
11、有反三角函数没有多项式计算方便，所以我们用多项式逼近方法解决此问题。 （具体理论证明见附录二）求得 ：5-ablttVrcsin122椭（5）ltTtann10当 k=1 时， 求得 近似是一个 2 阶的多项式，当 k=2 时，求得的 近似椭V椭V是一个五阶的多项式。由于拟合时考虑的误差，选择五阶多项式拟合效果较好，所以我们用当 k=2 时的五阶契比雪夫多项式逼近 椭V-椭 ablxx5348061572（6）由于储油罐的形状以及温度压强的关系这里的系数不太准确，我们可以根据题目（附件 1 中变位进油量数据）给出的累加进油量和油位高度的数据进行拟合。用 MATLAB 软件拟合后得到的结果是:3
12、.921-.46h76.1928 9.2h -0.684-.29h V3 5 ）（椭对该模型进行误差分析（误差分析的数据源：附件 1 位变位出油量数据）：将数据代入到上述 5 阶拟合多项式中与已知油量数据对比验证：验证数据方差和 为：2R-1.05712niiiVS（7） -1122niiiiiiVR（8）由以上可知拟合的效果比较好。、有变位时：椭圆油罐有纵向倾斜且倾斜角 ，选定的坐标系及有关参数如图 2 所o1.4示。分析发现变位后油高数据由于变位后产生影响，而油量数据我们认为是真实的，因此首先利用几何关系求出油量容积的解析表达式，利用最小二乘逼近6的方法求出真实反映油量的油高读数。然后利用
13、拟合的方法建立油高和油量的模型，从而解决标高问题。图（2） 1、 符号说明：标准差2S-纵向倾角-变位后液面在 Y 坐标轴上的分量H图中：c=0.4m, d=2.05m;设浮标的距离罐底的距离为 ，可以知道浮标的坐标为(c, ） 直线 ABfhfh的斜率为 ，所以可求得 AB 的直线方程为：tanL2: -tantcxyf（9）显然油液面的高度随横坐标 x 的变化为：-ttanhHf（10）yL2:x0 水平线Qcd7此时储油罐横截面的面积为：-dybbSH022)(a椭（11）abHb arcsin122所以变位时小椭油罐中油容量与油位高度的关系是：-ldxSV0椭椭（12） 根据公式（2）
14、我们可以准确的计算出不同油位高的储油量，但是，公式中还有反三角函数和开根号，我们计算起来很不方便，所以我们可以寻找一点做体积上的等效转换仿效上一问的方法。我们设那一点为 Q，那么它的坐标可设为 ，见图 (2)。,H并且该点 Q 也在直线 l1 上，根据 l1 的直线方程可以求得 为：H-tantan chf（13）对于一个固定外壳的容器，其中装定量体积的液体，无论如何对容器发生偏转，其液体的体积始终是一个固定值。根据等效的原理可以肯定点 Q 可以把楔状油形体积转化为标准的椭圆柱的体积，并且此时的油位高度为 。所以这样H就可以转化成没变位的情况去解决问题。得到转化后的椭圆面积为：-dybbaSH0221 )(椭
