1、摘 要本文通过对储油罐中油位高度及变位参数之间的不同情形的储油量进行分析并建立相应的数学模型，在该过程中先利用投影法、截面法及微元法得出储油量与油位高度及变位参数的函数关系。再由 Matlab 编程可知各高度储油量的理论数据，最后分析误差及评价模型的合理性。对于问题一的任一种情形，我们均建立笛卡尔坐标系，当储油罐无变位时，利用微元法得到体积关于 的公式，当储油罐发生变位时，根据储油罐中油量h的多少分成三种情形，就每一类利用微元法得到体积关于 的公式。代人附件h1 实验数据中的高度得到储油罐中的理论油量 。根据理论油量及实际油量得V出误差，判断误差所服从的分布，再利用相对误差进行误差分析并评价模
2、型的合理性。由上述得到储油罐发生变位时体积关于 的公式我们给出了罐体变位后油位高度间隔为 1cm 的罐容表标定值（即进/出油量与罐内油位高度的表格） 。对于问题二中的储油罐，我们先将问题进行简化考虑，得出了储油罐水平卧放时油量与浮油子高度的函数关系；再考虑储油量与油位高度及变位参数（纵向倾斜角度 和横向偏转角度 ）的一般情况，在该过程中，我们进行近似处理，利用投影法和截面法得出了储油量关于油位高度及变位参数的函数关系；并在固定的横向偏转角度 条件下，就纵向倾斜角度 的变化进行分成三类讨论，这三类又可以分成八种情形，得到了每一种情形下实际储油罐罐内储油量与油位高度的函数关系。在模型的改进中，我们
3、就问题二储油量与油位高度及变位参数的一般情况进行了仔细的考虑，将含油部分的体积分成四个部分，每一个部分将上述所提到的积分方法相结合，得到了各个部分的储油量与油位高度及变位参数的函数关系，从而可得总储油量与油位高度及变位参数的函数关系；并据此利用Matlab 编程和实际测量的数据求得 和 值；与此同时我们可以得出在固定 、值时各高度下的理论储油量；根据理论油量及实际油量得出误差，判断误差所服从的分布再利用相对误差进行误差分析并评价模型的合理性。由上述得到储油罐发生变位时体积关于 的公式我们给出了罐体变位后油位高度间隔为h10cm 的罐容表标定值。【关键词】投影法 截面法 微元法 Matlab 编
4、程11 问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统” ，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位） ，从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。（1）为了掌握罐体变位后对罐容表的影响，利用如图 4 的小椭圆型储油罐（两端平头的椭圆柱体
5、） ，分别对罐体无变位和倾斜角为 的纵向变位两.1种情况做了实验，实验数据如附件 1 所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为 1cm 的罐容表标定值。（2）对于图 1 所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度 和横向偏转角度 ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件 2） ，根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm 的罐容表标定值。进一步利用附件 2 中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。 （以上涉及的图
6、14 均在附录中）2 模型的假设与符号的约定2.1 模型的假设与说明（1）在储油罐倾斜的情况下，忽略油浮子高度为 0 时油所占的体积；（2）在储油罐倾斜的情况下，假设当油浮子高度达到最大后不再进油；（3）油的挥发速度很慢，忽略因油的挥发而造成储油量的减少；（4）储油罐的材料为钢体，忽略因渗出油而造成储油量的减少；（5）储油罐管理妥当，不会因特殊情况而造成储油量的变化。2.2 符号的约定与说明V表 示 储 油 罐 中 油 的 体 积 ;L表 示 储 油 罐 圆 柱 体 部 分 的 长 度a表 示 任 一 椭 球 截 面 的 长 半 轴 ;2b表 示 任 一 椭 球 截 面 的 短 半 轴 ;1a
7、表 示 油 浮 子 在 圆 柱 体 高 方 向 上 投 影 至 两 端 的 较 小 值 ;h表 示 油 浮 子 到 圆 柱 体 高 方 向 的 距 离 ;1表 示 储 油 罐 接 地 一 端 油 面 到 地 面 得 距 离表 示 纵 向 倾 斜 角 度 ;表 示 横 向 倾 斜 角 度0;h表 示 球 冠 高b表 示 球 冠 底 半 径3 问题的分析3.1 问题一的分析当储油罐无变位时，储油罐圆柱体的接地一端为原点，以圆柱体高方向为轴，建立笛卡尔坐标系，利用微元法得到体积关于 的公式，代人附件 1 实z h验数据中的高度得到储油罐中的理论油量 。利用附件 1 实验数据中得到储油V罐中的实际油量
8、，根据理论油量及实际油量就可以得出误差，判断误差所服从的分布，利用相对误差进行误差分析。当储油罐发生变位时，以储油罐圆柱体的接地一端为原点，圆柱体高方向为 轴，建立笛卡尔坐标系。根据储油罐中油量的多少分成三类，然后就每一z类利用微元法得到体积关于 的公式，代人附件 1 实验数据中的高度得到储油h罐中的理论油量 。利用附件 1 实验数据中得到储油罐中的实际油量，根据理V论油量及实际油量就可以得出误差，判断误差所服从的分布，利用相对误差进行误差分析。由上述得到储油罐发生变位时体积关于 的公式可以给出罐体变位后油位h高度间隔为 1cm 的罐容表标定值（即进/出油量与罐内油位高度的表格） 。3.2 问
9、题二的分析对于实际储油罐，我们首先将问题进行简化考虑，得出了当实际储油罐水平卧放时实际储油罐中油量与浮油子高度的函数关系；然后我们先考虑实际储油罐罐内储油量与油位高度及变位参数（纵向倾斜角度 和横向偏转角度 ）的一般情况，在该过程中，我们进行近似处理，利用投影法和截面法得出了储油量关于油位高度及变位参数的函数关系，再在固定的横向偏转角度 条件下，就纵向倾斜角度 的变化进行分类讨论，一共有三种情形，得到了每一种情形下实际储油罐罐内储油量与油位高度的函数关系。3最后我们先利用附件 2 中的少量实际数据得出了附件 2 所处状态下的纵向倾斜角度 和横向偏转角度 ，再利用附件 2 中给定各高度进行代人，
10、得到实际储油罐理论的储油量，与实际储油量进行比较，求出误差及相对误差。由上述得到储油罐发生变位时体积关于 的公式可以给出罐体变位后油位h高度间隔为 10cm 的罐容表标定值（即进/出油量与罐内油位高度的表格） 。4 模型的建立与求解4.1 问题一4.1.1 模型一当储油罐体无变位时，储油罐圆柱体的接地一端为原点，以圆柱体高方向为 轴，建立笛卡尔坐标系，如图 4-1 所示，在高度为 时，利用微元法过垂直z h轴的方向做截面 ，对 关于 进行积分，得到体积关于 的公式。()Sz()z图 4- 12 21.xyaxbyab由 2022() 1b)arcsinhSzdahb ()dVSz402 2()
11、 1barcsinLVSzdahbh 2 2()iLb利用 Matlab 中的命令 subs 代人附件 1 实验数据中的各高度得到储油罐中的理论油量 。V由附件 1 实验数据中进油量、出油量及储油罐罐内油量初值可以得到储油罐中的实际储油量，根据理论油量及实际油量就可以得出误差。由附录中的程序 youliang1，我们得到了理论储油量，误差及相对误差。进油后理论储油量与实际储油量随高度的变化规律如图 4-2 所示：图 4- 2出油后理论储油量与实际储油量随高度的变化规律如图 4-3 所示：5图 4- 3无变位进油和无变位出油的储油理论量和储油实际量及误差和相对误差的数据如表 4-1 所示：表 4
12、- 1无变位进油 无变位出油储油量实际值储油量理论值 (误差)i相对误差储油量实际值储油量理论值 (误差)i相对误差312 322.9 10.88 3.37% 3916.2 4052.8 136.61 3.37%362 374.6 12.63 3.37% 3866.2 4001.0 134.86 3.37%412 426.4 14.36 3.37% 3816.2 3949.3 133.13 3.37%462 478.1 16.13 3.37% 3766.2 3897.6 131.38 3.37%512 529.9 17.85 3.37% 3716.2 3845.8 129.65 3.37%56
13、2 581.6 19.61 3.37% 3666.2 3794.1 127.88 3.37%612 633.4 21.35 3.37% 3616.2 3742.3 126.16 3.37%662 685.1 23.08 3.37% 3566.2 3690.6 124.39 3.37%712 736.8 24.85 3.37% 3516.2 3638.8 122.65 3.37%762 788.6 26.58 3.37% 3466.2 3587.1 120.92 3.37%812 840.3 28.33 3.37% 3416.2 3535.4 119.17 3.37%862 892.1 30.0
14、6 3.37% 3366.2 3483.6 117.43 3.37%912 943.8 31.80 3.37% 3316.2 3431.9 115.70 3.37%962 995.5 33.54 3.37% 3266.2 3380.1 113.94 3.37%1012 1047.3 35.30 3.37% 3216.2 3328.4 112.18 3.37%1062 1099.1 37.05 3.37% 3166.2 3276.6 110.45 3.37%1112 1150.8 38.81 3.37% 3116.2 3224.9 108.70 3.37%1162 1202.6 40.55 3.
15、37% 3066.2 3173.2 106.96 3.37%1212 1254.3 42.29 3.37% 3016.2 3121.4 105.22 3.37%1262 1306.0 44.03 3.37% 2966.2 3069.7 103.46 3.37%61312 1357.8 45.77 3.37% 2916.2 3017.9 101.73 3.37%1362 1409.5 47.49 3.37% 2866.2 2966.2 100.00 3.37%1412 1461.2 49.24 3.37% 2816.2 2914.4 98.23 3.37%1462 1513.0 50.98 3.
16、37% 2766.2 2862.7 96.50 3.37%1512 1564.7 52.74 3.37% 2716.2 2811.0 94.76 3.37%1562 1616.5 54.49 3.37% 2666.2 2759.2 93.02 3.37%1612 1668.2 56.24 3.37% 2616.2 2707.5 91.26 3.37%1662 1720.0 57.98 3.37% 2566.2 2655.7 89.51 3.37%1712 1771.7 59.73 3.37% 2516.2 2603.9 87.76 3.37%1762 1823.5 61.46 3.37% 24
17、66.2 2552.2 86.03 3.37%1812 1875.2 63.19 3.37% 2416.2 2500.5 84.30 3.37%1862 1927.0 64.95 3.37% 2366.2 2448.7 82.53 3.37%1912 1978.7 66.68 3.37% 2316.2 2397.0 80.80 3.37%1962 2030.4 68.43 3.37% 2266.2 2345.3 79.07 3.37%2012 2082.2 70.20 3.37% 2216.2 2293.5 77.32 3.37%2062 2134.0 71.95 3.37% 2166.2 2
18、241.7 75.55 3.37%2112 2185.7 73.67 3.37% 2116.2 2190.0 73.84 3.37%2162 2237.4 75.43 3.37% 2066.2 2138.3 72.08 3.37%2212 2289.2 77.16 3.37% 2016.2 2086.5 70.32 3.37%2262 2340.9 78.89 3.37% 1966.2 2034.8 68.60 3.37%2312 2392.7 80.67 3.37% 1916.2 1983.0 66.85 3.37%2315.8 2396.6 80.78 3.37% 1866.2 1931.
19、3 65.12 3.37%2365.8 2448.4 82.54 3.37% 1816.2 1879.5 63.35 3.37%2367.1 2449.6 82.56 3.37% 1766.2 1827.8 61.61 3.37%2417.1 2501.4 84.34 3.37% 1716.2 1776.1 59.88 3.37%2467.1 2553.1 86.05 3.37% 1666.2 1724.3 58.12 3.37%2517.1 2604.9 87.82 3.37% 1616.2 1672.6 56.37 3.37%2567.1 2656.6 89.53 3.37% 1566.2
20、 1620.8 54.64 3.37%2617.1 2708.3 91.28 3.37% 1516.2 1569.1 52.87 3.37%2667 2760.0 93.03 3.37% 1466.2 1517.3 51.14 3.37%2668.8 2761.9 93.11 3.37% 1416.2 1465.6 49.42 3.37%2718.8 2813.7 94.83 3.37% 1366.2 1413.8 47.65 3.37%2768.8 2865.4 96.59 3.37% 1316.2 1362.1 45.91 3.37%2818.8 2917.2 98.34 3.37% 12
21、66.2 1310.4 44.18 3.37%2868.8 2968.9 100.09 3.37% 1216.2 1258.6 42.42 3.37%2918.8 3020.7 101.84 3.37% 1166.2 1206.9 40.68 3.37%2968.8 3072.4 103.58 3.37% 1116.2 1155.1 38.94 3.37%3018.8 3124.1 105.31 3.37% 1066.2 1103.4 37.19 3.37%3068.8 3175.9 107.06 3.37% 1016.2 1051.6 35.43 3.37%3118.8 3227.6 108
22、.80 3.37% 966.2 999.9 33.71 3.37%3168.8 3279.4 110.55 3.37% 916.2 948.2 31.96 3.37%3168.9 3279.5 110.55 3.37% 866.2 896.4 30.20 3.37%3218.9 3331.2 112.27 3.37% 816.2 844.6 28.45 3.37%3268.9 3382.9 114.03 3.37% 766.2 792.9 26.71 3.37%73318.9 3434.7 115.76 3.37% 716.2 741.2 24.99 3.37%3368.9 3486.4 11
23、7.52 3.37% 666.2 689.4 23.25 3.37%3418.9 3538.2 119.26 3.37% 616.2 637.7 21.51 3.37%3468.9 3589.9 121.01 3.37% 566.2 585.9 19.74 3.37%3518.9 3641.7 122.76 3.37% 516.2 534.2 18.02 3.37%3568.9 3693.4 124.51 3.37% 466.2 482.5 16.28 3.37%3618.9 3745.1 126.23 3.37% 416.2 430.7 14.50 3.37%3668.9 3796.9 12
24、7.98 3.37% 366.2 379.0 12.77 3.37%3718.9 3848.6 129.74 3.37% 316.2 327.2 11.02 3.37%3768.9 3900.4 131.48 3.37% 266.2 275.5 9.28 3.37%3818.9 3952.1 133.23 3.37%3868.9 4003.9 134.95 3.37%3918.9 4055.6 136.70 3.37%3968.9 4107.4 138.45 3.37%由上述的表格可以得出相对误差稳定，不会随高度发生变化。当储油罐体发生变位时，我们就油面及油浮子的位置与倾斜角度的变化情况分成将该问题三类，然后就每一类进行分析，找出了储油罐中油的体积与油浮子位置的函数关系，其立体图如图 4-4 所示：图 4- 4（1）当油浮子的高度在 时，如图 4-5 所示，截面图如图10tanhz4-6 所示8图 4- 5 图 4- 6 截面图截面面积为： 1tan220 21112 21()()t)tan()tan()t()arcsinhzSbydhzhzbzb 所以油罐体积为：（1.1） 111tan0ta 2110tan2 210()()t2(tan()tan()ta)rcsihhVfSzdbhzhzbzb 将 代入公式（1.1）1tah，则公式（1.1）nzbt令
