储油罐的变位识别与罐容表标定.docx
《储油罐的变位识别与罐容表标定.docx》由会员分享，可在线阅读，更多相关《储油罐的变位识别与罐容表标定.docx（12页珍藏版）》请在冰豆网上搜索。
储油罐的变位识别与罐容表标定
储油罐的变位识别与罐容表标定
摘要：
本论文主要针对储油罐的变位识别与罐容表标定进行研究和讨论的。
储油罐在我们的现实生活中非常的重要，特别是它的安全存储非常的重要，所以这是一个很好很实际的研究课题。
首先采用了统计方法对数据进行了分析，将数据按无变位和有变位综合分析体积与油品液面度的函数关系,并通过测量油罐实际参数,使该函数具体化,给出了该函数的计算机程序,用maTLAB软件，通过拟合法，得出了油品体积与油品液面高度对照表,提高了油品体积的测量精度,减少了实际测量中的工作量。
模型一：
本模型的建立是研究在无变位情况下，油面高度h与油的体积v之间的函数关系。
通过大量的数据录入，得出十分逼近的拟合图像和函数。
注意的是录入的数据越多，拟合出的图像与真实的图像就越接近，误差就越小。
模型二：
它主要是研究油罐在变位情况下，在模型一的基础上，建立油位高度与油的体积之间的函数关系。
对模型进行了充分的验证和合理的理论推导，所求出的结果都在理论的误差值范围内。
然后借助MATLAB软件的拟合方法进行分析，对附件中的数据进行了筛选，去除异常的数据，对残缺数据进行适当的补充。
结果显示，理论结果跟模拟数据结果基本吻合
模型三：
它是建立在前两个模型基础上的。
通过对前两个模型的研究，来讨论更加复杂的情况。
在有三个变化量的情况下，如何研究呢？
首先我求出在没有变为下，应用积分学的知识求出油面高度与油的体积之间的一般函数关系式（见模型三内的（9））然后再用MATLAB的拟合方法帮助我们拟合出在横向与纵向都变化的情况下，h与v之间的函数关系
关键词：
卧式储油罐 油品体积 变位识别 标定 MATLAB数据的拟合统计学
一、问题重述
现实生活中，许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。
按照有关规定，需要定期对罐容表进行重新标定。
图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。
图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。
（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。
请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。
（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度）之间的一般关系。
请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。
进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。
二、问题的分析
在石油的产地和加工厂，为了储存原油，经常使用大量的储油罐，油罐的外形由三部分组成，中间是一个平放圆柱体，两边是由一个椭球体沿对称轴一分为二的两个半椭球体焊接在圆柱体两端组合成的，上端有住有孔，由于经常油和取油，很难知道注入油和取出油的数量，这给注油和取油数量带来很大麻烦。
因此，要给油桶上安装一个精细的刻度计，工人根据刻度计就能知道每次注入和取出油的量的多少。
如何设计刻度计呢？
这是一个很实际的问题。
下面我将用maTLAB和积分的方法解决这个问题。
三、模型的假定
1．假设所给油罐都是完好的，没有破损；
2.假设题目中所给的数据真实可靠；
3.假设当测量罐体的直径时，罐体的内外经长度相同；（见模型三中的（9））
4.假设油位刻度的变化只与有关的体积和油罐的倾斜度有关，与温度、气压等因素无关；
5.假设每次油量的注入或抽出不会引起油罐随着油量的变化而发生倾斜；
6.假设油罐的焊接处都是圆滑的，圆柱体是标准的，两边的椭球体也是标准的。
四、模型的建立与求解
4.1模型一：
若油罐的高为h时，此时油罐的出油量为v，那么问题就归结为求出函数
h=h（V）
（1）
油罐的高度和油量的函数关系如下：
由MATLAB的拟合法，通过已经给出的数据在MATLAB中进行你和运算
体位无变化时：
x=[159.02,176.14,192.59,208.50,223.93,238.97,253.66,268.04,282.16,296.03,309.69,323.15,336.44,349.57,362.56,375.42,388.16,400.97,413.32,425.76,438.12,450.40,462.62,474.78,486.89,498.95,510.97,522.95,534.90,546.82,558.72,570.61,582.48,594.35,606.22,618.09,629.96,641.85,653.75,665.67,677.63,678.54,690.53,690.82,702.85,714.91,727.03,739.19,751.42,763.70,764.16,776.53,788.99,801.54,814.19,826.95,839.83,852.84,866.00,879.32,892.82,892.84,906.53,920.45,934.61,949.05,963.80,978.91,994.43,1010.43,1026.99,1044.25,1062.37,1081.59,1102.33,1125.32,1152.36,1193.49];
y=[312,362,412,462,512,562,612,662,712,762,812,862,912,962,1012,1062,1112,1162,1212,1262,1312,1362,1412,1462,1512,1562,1612,1662,1712,1762,1812,1862,1912,1962,2012,2062,2112,2162,2212,2262,2312,2315.83,2365.83,2367.06,2417.06,2467.06,2517.06,2567.06,2617.06,2666.98,2668.83,2718.83,2768.83,2818.83,2868.83,2918.83,2968.83,3018.83,3068.83,3118.83,3168.83,3168.91,3218.91,3268.91,3318.91,3368.91,3418.91,3468.91,3518.91,3568.91,3618.91,3668.91,3718.91,3768.91,3818.91,3868.91,3918.91,3968.91];
p=polyfit（x,y,3）
p=
-0.00000240774433 0.00431658990406 1.66987172334547-49.329
得到三次拟合多项式：
-0.00000240774433x
+0.00431658990406
+1.66987172334547x-49.329
xi=159.02:
100:
1193.49;
yi=polyval（p,xi）;
plot（x,y,'.',xi,yi）;
我们得到拟合之后油位高度h与储油量v的函数图：
（关系图像如图1）
（图1）
即h=h（v）=-0.0000024x^3+0.0043x^2+1.6699x-49.1864;
4.2模型二：
倾斜4.1°时，油罐上的刻度h与储油量之间的函数关系
倾斜变位时：
x=[411.29,423.45,438.33,450.54,463.90,477.74,489.37,502.56,514.69,…
526.84,538.88,551.96,564.40,576.56,588.74,599.56,611.62,623.44,…
635.58,646.28,658.59,670.22,680.63,693.03,704.67,716.45,727.66,…
739.39,750.90,761.55,773.43,785.39,796.04,808.27,820.80,832.80,…
844.47,856.29,867.60,880.06,892.92,904.34,917.34,929.90,941.42,…
954.60,968.09,980.14,992.41,1006.34,1019.07,1034.24,1035.36];
y=[962.86,1012.86,1062.86,1112.86,1162.86,1212.86,1262.86,1312.79,…
1362.79,1412.73,1462.73,1512.73,1562.73,1612.73,1662.73,1712.73,…
1762.73,1812.73,1862.73,1912.73,1962.73,2012.73,2062.73,2112.73,…
2162.73,2212.73,2262.73,2312.73,2362.73,2412.73,2462.73,2512.73,…
2562.73,2612.73,2662.73,2712.73,2762.73,2812.73,2862.73,2912.73,…
2962.73,3012.73,3062.73,3112.73,3162.73,3212.73,3262.73,3312.73,…
3362.73,3412.73,3462.73,3512.73,3514.74];
h=polyfit（x,y,3）
h=
-0.00000249721671 0.00539373083922 0.44909628838092 42.10407086898768
得到三次拟合多项式-0.0000025x^3+0.0054x^2+0.449096x+42.104071；
我们得到拟合之后油位高度h与储油量v的函数图：
（图2）
（图2）
将在无变位情况下油位高度与罐内储油量关系与4.1°倾斜变位后比较
无变位与倾斜变位（4.1°）油位高度比较
罐内储油量
实际无变位油位高度
函数计算值（无变位）
计算与给出数据的差值
函数计算值（倾斜变位）
实际倾斜与计算出来结果的差值
312
159.02
157.5046096
1.515390431
217.7106701
60.20606054
362
176.14
174.0726887
2.067311258
233.8631093
59.79042058
412
192.59
190.3297781
2.260221877
249.7732987
59.4435206
462
208.5
206.2856849
2.214315109
265.4477296
59.16204469
512
223.93
221.9502162
1.979783774
280.8928931
58.94267691
562
238.97
237.3331793
1.636820693
296.1152806
58.78210133
612
253.66
252.4443813
1.215618685
311.1213833
58.67700203
662
268.04
267.2936294
0.746370569
325.9176925
58.62406307
712
282.16
281.8907308
0.269269167
340.5106994
58.61996852
762
296.03
296.2454927
-0.215492702
354.9068952
58.66140246
812
309.69
310.3677222
-0.677722218
369.1127712
58.74504895
862
323.15
324.2672266
-1.117226561
383.1348186
58.86759206
912
336.44
337.9538129
-1.513812911
396.9795288
59.02571587
962
349.57
351.4372884
-1.867288447
410.6533929
59.21610444
1012
362.56
364.7274604
-2.167460351
424.1629022
59.43544184
1062
375.42
377.8341358
-2.414135802
437.5145479
59.68041215
1112
388.16
390.767122
-2.607121979
450.7148214
59.94769942
1162
400.79
403.5362261
-2.746226063
463.7702138
60.23398775
1212
413.32
416.1512552
-2.831255234
476.6872164
60.53596118
1262
425.76
428.6220167
-2.862016673
489.4723205
60.8503038
1312
438.12
440.9583176
-2.838317558
502.1320172
61.17369967
1362
450.4
453.1699651
-2.769965069
514.6727979
61.50283286
1412
462.62
465.2667664
-2.646766388
527.1011538
61.83438745
1462
474.78
477.2585287
-2.478528694
539.4235762
62.1650475
1512
486.89
489.1550592
-2.265059166
551.6465562
62.49149708
1562
498.95
500.966165
-2.016164986
563.7765853
62.81042027
1612
510.97
512.7016533
-1.731653332
575.8201545
63.11850113
1662
522.95
524.3713314
-1.421331386
587.7837551
63.41242373
1712
534.9
535.9850063
-1.085006326
599.6738785
63.68887214
1762
546.82
547.5524853
-0.732485333
611.4970158
63.94453044
1812
558.72
559.0835756
-0.363575587
623.2596583
64.17608269
1862
570.61
570.5880843
0.021915732
634.9682972
64.38021297
1912
582.48
582.0758186
0.404181444
646.6294239
64.55360534
1962
594.35
593.5565856
0.79341437
658.2495295
64.69294387
2012
606.22
605.0401927
1.179807328
669.8351053
64.79491263
2062
618.09
616.5364469
1.55355314
681.3926426
64.8561957
2112
629.96
628.0551554
1.904844624
692.9286325
64.87347714
2162
641.85
639.6061254
2.243874602
704.4495664
64.84344102
2212
653.75
651.1991641
2.550835893
715.9619355
64.76277142
2262
665.67
662.8440787
2.825921317
727.4722311
64.6281524
2312
677.63
674.5506763
3.079323694
738.9869443
64.43626803
2315.83
678.54
675.450225
3.089774974
739.8693388
64.41911378
2365.83
690.53
687.2341934
3.295806585
751.3960643
64.16187088
2367.06
690.82
687.5250729
3.294927124
751.6798207
64.15474782
2417.06
702.85
699.3932319
3.456768106
763.2249657
63.83173381
2467.06
714.91
711.353488
3.556512015
774.7946591
63.44117108
2517.06
727.03
723.4156483
3.61435167
786.3953921
62.97974372
2567.06
739.19
735.5895201
3.600479892
798.0336559
62.44413579
2617.06
751.42
747.8849105
3.535089501
809.7159419
61.83103136
2666.98
763.7
760.2916338
3.408366172
821.4299249
61.13829108
2668.83
764.16
760.7540611
3.405938868
821.8650887
61.1110276
2718.83
776.53
773.3270861
3.20291387
833.6570296
60.32994347
2768.83
788.99
786.0513985
2.938601544
845.5126961
59.46129766
2818.83
801.54
798.9368053
2.603194711
857.4385795
58.50177425
2868.83
814.19
811.9931138
2.19688619
869.4411711
57.44805731
2918.83
826.95
825.2301312
1.719868801
881.5269621
56.29683091
2968.83
839.83
838.6576646
1.172335365
893.7024437
55.04477911
3018.83
852.84
852.2855213
0.554478701
905.9741073
53.68858599
3068.83
866
866.1235084
-0.123508371
918.348444
52.22493562
3118.83
879.32
880.181433
-0.86143303
930.8319451
50.65051206
3168.83
892.82
894.4691025
-1.649102457
943.4311019
48.96199939
3168.91
892.84
894.4921521
-1.652152061
943.4513566
48.95920458
3218.91
906.53
909.0197646
-2.489764578
956.1728609
47.15309637
3268.91
920.45
923.7967519
-3.346751913
969.0230138
45.22626187
3318.91
934.61
938.8329212
-4.222921248
982.0083064
43.17538517
3368.91
949.05
954.1380798
-5.088079761
995.1352301
40.99715033
3418.91
963.8
969.7220346
-5.922034634
1008.410276
38.68824142
3468.91
978.91
985.594593
-6.684593046
1021.839936
36.24534251
3518.91
994.43
1001.765562
-7.335562178
1035.4307
33.66513767
3568.91
1010.43
1018.244749
-7.814749208
1049.18906
30.94431097
3618.91
1026.99
1035.041961
-8.051961318
1063.121508
28.07954649
3668.91
1044.25
1052.167006
-7.917005687
1077.234534
25.06752829
3718.91
1062.37
1069.629689
-7.259689495
1091.53463
21.90494044
3768.91
1081.59
1087.43982
-5.849819922
1106.028287
18.58846701
3818.91
1102.33
1105.607204
-3.277204148
1120.721996
15.11479207
3868.91
1125.32
1124.141649
1.178350646
1135.622249
11.48059969
3918.91
1152.36
1143.052963
9.307037282
1150.735537
7.68257395
3968.91
1193.49
1162.350951
31.13904858
1166.06835
3.717398907
平均值
-0.002326173
53.88973544
升级会员 