《储油罐的变位识别与罐容表标定数模竞赛论文-毕业论文.doc》由会员分享，可在线阅读，更多相关《储油罐的变位识别与罐容表标定数模竞赛论文-毕业论文.doc（51页珍藏版）》请在金锄头文库上搜索。
1、储油罐的变位识别与罐容表标定2010高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参
2、赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 2010 年 9 月 13 日赛区评阅编号（由赛区组委会评阅前进行编号）：2008高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号储油罐的变位识别与罐容表标定摘要本文是与储油罐的变位识别与罐容表标定相关的数学问题。问题一中，我们运用高等数学中的微元的思
3、想，分别建立了储油罐在变前位、变位后的数学模型，并对其求解，得出变位前与变位后，在相同油高下，间隔为1cm的罐容表标定值（详见39页）。油体高度变位后的体积变位前的体积121.6016.1718222.2024.46211093658.103786.66201103687.703809.9516通过对变位前、变位后标定值的对比，在时得出倾斜对罐容表的影响为：倾斜之后的油罐，油体体积的标定值在任意高度下都要比倾斜之前标定体积要小。当储油罐内装油过多或过少时，倾斜对油体体积的标定值影响较小；当储油罐内油面深度接近罐体中线时影响最大。问题二中，储油罐发生纵向变位的同时还发生横向变位。通过几何分析找出
4、倾角、偏角与油体高度间的函数关系。按照问题一中的分析方式，利用微元法得到体积与深度、纵向倾角和横向偏角的关系式，但方程比较复杂，不能够得到解析解。因此我们采用近似算法，找出相邻油层间的体积，利用最优化理论得到、的最优解，再将代入上述关系式中求解，即得到罐体变位后油位高度间隔为的罐容表标定值。(详见41页)高度标定值高度标定值403.68183626058.4239505.41510627060.2533607.35832928061.851我们用和高度算出体积，计算出理论出油体积，与实际出油体积作对比，算出相对误差均值在0.1%左右，证明我们模型的准确性较高。最后，我们用蒙特卡洛法对模型中的高
5、度进行扰动时，发现对高度添加一个服从正态分布随机扰动值时产生的相对误差与的值产生的相对误差在一个数量级。从而证明了我们所建立的模型具有较好稳定性。一、问题重述（略）二、问题分析加油站的地下储油罐通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高和储油量的变化情况。罐体发生变位后会发生纵向倾角和横向偏角，影响油浮子的位置，致使读出的油位高度与实际高度不符，从而导致罐容表发生改变。通过对附表1数据的分析，可以看出进油曲线和出油曲线完全重合，即没有漏油存在，且流量计读数完全正确。我们不考虑温度、压力及不计浮盘上下浮动对罐内油量体积的影响，但要考虑到油罐内有注油管道
6、和出油管道（未达油罐底部）以及油位探测管道（焊接至油罐底部），会占一部分体积。考虑到在实际情况中，压力及温度对油罐变形的影响，为确保安全储油罐不可能完全装满，也不可能在油量很少的情况下供油，且储油罐的变位不可能很大，纵向倾角和横向偏角都在小范围内。简化的储油罐如下图，我们只考虑在油量在阴影部位的情况。 简化的储油罐2.1对问题一的分析：我们首先找出油罐倾斜前、倾斜后体积与倾角，高度间的函数关系，然后通过互相对比得出倾角对罐容表的影响程度。未倾斜的油罐：如图：对于横截面为椭圆形的储油罐，封闭的罐体内油面为水平面，在统一的坐标系下，油罐内处处油体深度相等，且处处横截面面积相等。只需计算出任一截面面
7、积，并对其长度方向积分即可求出体积。通过得到的体积与高度间的关系式，即可得到任意高度时储油罐中的理论储油量。但实际中储油罐内含有输油、进油管道等其他影响油体体积的因素。其标定值可能不等于理论值。我们通过修正得到它们之间的差值，从而找到正确的标定值。倾斜的油罐：罐体纵向倾斜时，油深不再是处处相等，各截面面积也此时倾角、油体高度对体积都有影响。我们打算通过微元法建立他们的函数关系，对其进行积分求解。同样找出修正值，从而得出油罐倾斜时的正确标定值。2.2对问题二的分析：附件2中给出的数据是在倾角和偏角同时影响下的实际检测数据。同时考虑倾角和偏角对罐容表标定的影响，为了简化数学建模的求解，我们分为两步
8、求解，为此我们首先考虑横向偏转的情况，找出偏转之后油体的实际高度与读数高度、偏转角之间的函数关系。将读数高度用偏转角和油体实际高度的函数表达，即建立：这样既可将偏转后实际油体高度转为读数高度，再考虑纵向变位时，就可将问题转化为问题一的情况。从而得到的关系式。但问题二中的储油罐的主体是圆柱体，两端为球冠体。这就使得建立积分方程比较复杂，无法得出解析解。我们采用两种方法来解决这个问题。1. 对积分方程进行数值解；2. 采用近似法求得两油层间体积的近似解，步骤如下：1）首先找出液面所在的截面面积与高度、变位角度之间的函数关系；2）取上下两层距离为的两油层，相邻高度的体积可以用一个近似柱体来代替，其计
9、算公式为：表示上层油面面积，表示下层油面面积，表示两油层间的间距 3）由题中所给数据，得到相邻高度间的采集出油量，用最小二乘法思想建立优化模型，可得到的估计值。把高度h、代入积分函数，便可以算出任意高度体积。从而给出油位高度间隔为的罐容表标定值。相比之下，第二种方法比第一种方法更易编程实现，因此我们采用第二种方式进行建模。三、模型假设1、流量计所测数据准确；2、储油罐为刚体，即油罐不会变形3、忽略因温度变化而对油体容积的影响；4、不计油品的静压力作用而引起得油体的膨胀和收缩的容积；5、油浮子与标尺为点接触，不会因为油浮子影响油体体积；6、罐体发生变位时纵向倾角的取值范围为，横向偏角的取值范围为
10、四、符号说明-出油罐内油体的读数高度-油体的实际高度-椭圆的长半轴-椭圆的短半轴-储油罐的长度-按照读数高度得出的油罐内油体体积-纵向倾角-横向偏角-椭圆横截面面积-两相邻为的油层间的体积-上层油面面积-下层油面面积-两油层间的间距-真实体积部分变量在使用时说明 五、模型的建立与求解5.1问题（一）微分方程模型5.1.1油罐未变位时的模型根据对前面对问题一的分析，首先建立了未变位时，储油罐内油量与油位高度的数学程模型1。建立如图5.1的坐标系，设油罐内液面高度为，为恒定值；椭圆的长半轴为,短半轴为，xyh简化小椭圆储油罐AAAA图5.1椭圆轨迹方程为： (5.1)对半弓横截面（阴影部分）进行积
11、分得出面积： (5.2)由（1）式得出为： (5.3)将(3)式代入（2）式得： (5.4)积分后可得： (5.5)面积再乘以长度即为体积： (5.6)其中，半弓椭圆横截面面积，表示油罐长度，表示油罐内油体体积。代入已知值以及变值，用求解得出结果部分见表1，详见附录1。累加进油量/L实际累加体积/L油位高度/mm油高/m计算体积/m3实际体积/m3差值50312159.020.159020.322882590.3120.01088259100362176.140.176140.3746330340.3620.012633034365739191152.361.152364.0556137683
12、.918910.136703768370739691193.491.193494.1073621283.968910.138452128累加出油量/L实际累加体积/L油位高度/mm油高/m计算体积/ m3实际体积/ m3差值52.7239161150.721.150724.0527979523.916190.136607952102.738661123.991.123994.0010492213.866190.1348592213653316.2160.480.160480.327208490.316190.011018493703266.2142.620.142620.2754733340.266190.009283334表1注：实际累加体积是在罐内油量初值的基础上累次加上每次进油量得出，计算体积是按照上述函数式（5）代值求解的。不考虑其他任何影响因素的情况下，油深应该与高度的关系为（6）式所得积分值。则计算体积应该与实际体积相等，没有差值产生。但由结果可见计算体积与实际体积之间有差值。对题中所给附件1中
