大学生数学建模竞赛a题获奖论文-储油罐的变位识别与罐容表标定-本科论文.doc

摘要本文通过对储油罐中油位高度及变位参数之间的不同情形的储油量进行分析并建立相应的数学模型，在该过程中先利用投影法、截面法及微元法得出储油量与油位高度及变位参数的函数关系。再由Matlab编程可知各高度储油量的理论数据，最后分析误差及评价模型的合理性。对于问题一的任一种情形，我们均建立笛卡尔坐标系，当储油罐无变位时，利用微元法得到体积关于的公式，当储油罐发生变位时，根据储油罐中油量的多少分成三种情形，就每一类利用微元法得到体积关于的公式。代人附件1实验数据中的高度得到储油罐中的理论油量。根据理论油量及实际油量得出误差，判断误差所服从的分布，再利用相对误差进行误差分析并评价模型的合理性。由上述得到储油罐发生变位时体积关于的公式我们给出了罐体变位后油位高度间隔为1cm的罐容表标定值（即进/出油量与罐内油位高度的表格）。对于问题二中的储油罐，我们先将问题进行简化考虑，得出了储油罐水平卧放时油量与浮油子高度的函数关系；再考虑储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度）的一般情况，在该过程中，我们进行近似处理，利用投影法和截面法得出了储油量关于油位高度及变位参数的函数关系；并在固定的横向偏转角度条件下，就纵向倾斜角度的变化进行分成三类讨论，这三类又可以分成八种情形，得到了每一种情形下实际储油罐罐内储油量与油位高度的函数关系。在模型的改进中，我们就问题二储油量与油位高度及变位参数的一般情况进行了仔细的考虑，将含油部分的体积分成四个部分，每一个部分将上述所提到的积分方法相结合，得到了各个部分的储油量与油位高度及变位参数的函数关系，从而可得总储油量与油位高度及变位参数的函数关系；并据此利用Matlab编程和实际测量的数据求得和值；与此同时我们可以得出在固定、值时各高度下的理论储油量；根据理论油量及实际油量得出误差，判断误差所服从的分布再利用相对误差进行误差分析并评价模型的合理性。由上述得到储油罐发生变位时体积关于的公式我们给出了罐体变位后油位高度间隔为10cm的罐容表标定值。【关键词】投影法 截面法 微元法 Matlab编程1问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。（以上涉及的图14均在附录中）2模型的假设与符号的约定2.1模型的假设与说明（1）在储油罐倾斜的情况下，忽略油浮子高度为0时油所占的体积；（2）在储油罐倾斜的情况下，假设当油浮子高度达到最大后不再进油；（3）油的挥发速度很慢，忽略因油的挥发而造成储油量的减少；（4）储油罐的材料为钢体，忽略因渗出油而造成储油量的减少；（5）储油罐管理妥当，不会因特殊情况而造成储油量的变化。2.2符号的约定与说明3问题的分析3.1问题一的分析当储油罐无变位时，储油罐圆柱体的接地一端为原点，以圆柱体高方向为轴，建立笛卡尔坐标系，利用微元法得到体积关于的公式，代人附件1实验数据中的高度得到储油罐中的理论油量。利用附件1实验数据中得到储油罐中的实际油量，根据理论油量及实际油量就可以得出误差，判断误差所服从的分布，利用相对误差进行误差分析。当储油罐发生变位时，以储油罐圆柱体的接地一端为原点，圆柱体高方向为轴，建立笛卡尔坐标系。根据储油罐中油量的多少分成三类，然后就每一类利用微元法得到体积关于的公式，代人附件1实验数据中的高度得到储油罐中的理论油量。利用附件1实验数据中得到储油罐中的实际油量，根据理论油量及实际油量就可以得出误差，判断误差所服从的分布，利用相对误差进行误差分析。由上述得到储油罐发生变位时体积关于的公式可以给出罐体变位后油位高度间隔为1cm的罐容表标定值（即进/出油量与罐内油位高度的表格）。3.2问题二的分析对于实际储油罐，我们首先将问题进行简化考虑，得出了当实际储油罐水平卧放时实际储油罐中油量与浮油子高度的函数关系；然后我们先考虑实际储油罐罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度）的一般情况，在该过程中，我们进行近似处理，利用投影法和截面法得出了储油量关于油位高度及变位参数的函数关系，再在固定的横向偏转角度条件下，就纵向倾斜角度的变化进行分类讨论，一共有三种情形，得到了每一种情形下实际储油罐罐内储油量与油位高度的函数关系。最后我们先利用附件2中的少量实际数据得出了附件2所处状态下的纵向倾斜角度和横向偏转角度，再利用附件2中给定各高度进行代人，得到实际储油罐理论的储油量，与实际储油量进行比较，求出误差及相对误差。由上述得到储油罐发生变位时体积关于的公式可以给出罐体变位后油位高度间隔为10cm的罐容表标定值（即进/出油量与罐内油位高度的表格）。4模型的建立与求解4.1问题一4.1.1模型一当储油罐体无变位时，储油罐圆柱体的接地一端为原点，以圆柱体高方向为轴，建立笛卡尔坐标系，如图4-1所示，在高度为时，利用微元法过垂直轴的方向做截面，对关于进行积分，得到体积关于的公式。图4- 1利用Matlab中的命令subs代人附件1实验数据中的各高度得到储油罐中的理论油量。由附件1实验数据中进油量、出油量及储油罐罐内油量初值可以得到储油罐中的实际储油量，根据理论油量及实际油量就可以得出误差。由附录中的程序youliang1，我们得到了理论储油量，误差及相对误差。进油后理论储油量与实际储油量随高度的变化规律如图4-2所示：图4- 2出油后理论储油量与实际储油量随高度的变化规律如图4-3所示：图4- 3无变位进油和无变位出油的储油理论量和储油实际量及误差和相对误差的数据如表4-1所示：表4- 1无变位进油无变位出油储油量实际值储油量理论值(误差)相对误差储油量实际值储油量理论值(误差)相对误差312322.910.883.37%3916.24052.8136.613.37%362374.612.633.37%3866.24001.0134.863.37%412426.414.363.37%3816.23949.3133.133.37%462478.116.133.37%3766.23897.6131.383.37%512529.917.853.37%3716.23845.8129.653.37%562581.619.613.37%3666.23794.1127.883.37%612633.421.353.37%3616.23742.3126.163.37%662685.123.083.37%3566.23690.6124.393.37%712736.824.853.37%3516.23638.8122.653.37%762788.626.583.37%3466.23587.1120.923.37%812840.328.333.37%3416.23535.4119.173.37%862892.130.063.37%3366.23483.6117.433.37%912943.831.803.37%3316.23431.9115.703.37%962995.533.543.37%3266.23380.1113.943.37%10121047.335.303.37%3216.23328.4112.183.37%10621099.137.053.37%3166.23276.6110.453.37%11121150.838.813.37%3116.23224.9108.703.37%11621202.640.553.37%3066.23173.2106.963.37%12121254.342.293.37%3016.23121.4105.223.37%12621306.044.033.37%2966.23069.7103.463.37%13121357.845.773.37%2916.23017.9101.733.37%13621409.547.493.37%2866.22966.2100.003.37%14121461.249.243.37%2816.22914.498.233.37%14621513.050.983.37%2766.22862.796.503.37%15121564.752.743.37%2716.22811.094.763.37%15621616.554.493.37%2666.22759.293.023.37%16121668.256.243.37%2616.22707.591.263.37%16621720.057.983.37%2566.22655.789.513.37%17121771.759.733.37%2516.22603.987.763.37%17621823.561.463.37%2466.22552.286.033.37%18121875.263.193.37%2416.22500.584.303.37%18621927.064.953.37%2366.22448.782.533.37%19121978.766.683.37%2316.22397.080.803.37%19622030.468.433.37%2266.22345.379.073.37%20122082.270.203.37%2216.22293.577.323.37%20622134.071.953.37%2166.22241.775.553.37%21122185.773.673.37%2116.22190.073.843.37%21622237.475.433.37%2066.22138.372.083.37%22122289.277.163.37%2016.22086.570.323.37%22622340.978.893.37%1966.22034.868.603.37%23122392.780.673.37%1916.21983.066.853.37%2315.82396.680.783.37%1866.21931.365.123.37%2365.82448.482.543.37%1816.21879.563.353.37%2367.12449.682.563.37%1766.21827.861.613.37%2417.12501.484.343.37%1716.21776.159.883.37%2467.12553.186.053.37%1666.21724.358.123.37%2517.12604.987.823.37%1616.21672.656.373.37%2567.12656.689.533.37%1566.21620.854.643.37%2617.12708.391.283.37%1516.21569.152.873.37%26672760.093.033.37%1466.21517.351.143.37%2668.82761.993.113.37%1416.21465.649.423.37%2718.82813.794.833.37%1366.21413.847.653.37%2768.82865.496.593.37%1316.21362.145.913.37%2818.82917.298.343.37%1266.21310.444.183.37%2868.82968.9100.093.37%1216.21258.642.423.37%2918.83020.7101.843.37%1166.21206.940.683.37%2968.83072.4103.583.37%1116.21155.138.943.37%3018.83124.1105.313.37%1066.21103.437.193.37%3068.83175.9107.063.37%1016.21051.635.433.37%3118.83227.6108.803.37%966.2999.933.713.37%3168.83279.4110.553.37%916.2948.231.963.37%3168.93279.5110.553.37%866.2896.430.203.37%3218.93331.2112.273.37%816.2844.628.453.37%3268.93382.9114.033.37%766.2792.926.713.37%3318.93434.7115.763.37%716.2741.224.993.37%3368.93486.4117.523.37%666.2689.423.253.37%3418.93538.2119.263.37%616.2637.721.513.37%3468.93589.9121.013.37%566.2585.919.743.37%3518.93641.7122.763.37%516.2534.218.023.37%3568.93693.4124.513.37%466.2482.516.283.37%3618.93745.1126.233.37%416.2430.714.503.37%3668.93796.9127.983.37%366.2379.012.773.37%3718.93848.6129.743.37%316.2327.211.023.37%3768.93900.4131.483.37%266.2275.59.283.37%3818.93952.1133.233.37%3868.94003.9134.953.37%3918.94055.6136.703.37%3968.94107.4138.453.37%由上述的表格可以得出相对误差稳定，不会随高度发生变化。当储油罐体发生变位时，我们就油面及油浮子的位置与倾斜角度的变化情况分成将该问题三类，然后就每一类进行分析，找出了储油罐中油的体积与油浮子位置的函数关系，其立体图如图4-4所示：图4- 4（1）当油浮子的高度在时，如图4-5所示，截面图如图4-6所示图4- 5图4- 6截面图截面面积为：所以油罐体积为：（1.1）将代入公式（1.1），则公式（1.1）（2）当油浮子的高度在，如图4-7所示， 图4- 7截面面积为：(0.1)所以油罐装油的体积为；(0.2)将公式(0.3)，：，则（3）当油浮子的高度时，如图4-8所示，图4- 8其中，此时截面面积为：所以油罐装油的体积为则上式公式变为：利用Matlab中的命令subs代人附件1实验数据中的各高度得到储油罐中的理论油量。由附件1实验数据中进油量、出油量及储油罐罐内油量初值可以得到储油罐中的实际油量，根据理论油量及实际油量就可以得出误差。由附录中的程序youliang2，我们得到了理论油量，误差及相对误差。倾斜时进油后理论储油量与实际储油量随高度的变化规律如图4-9所示：图4- 9倾斜时出油后理论储油量与实际储油量随高度的变化规律如图4-10所示：图4- 10将倾斜进油和倾斜变位出油的储油理论值和储油实际值及误差和相对误差的数据如表4-2所示：表4- 2倾斜变位进油倾斜变位出油储油量实际量储油量理论值(误差)相对误差储油量实际值储油量理论值(误差)相对误差962.91010.07.194.67%3464.73523.158.371.66%1012.91058.345.474.30%3414.73478.163.351.82%1062.91118.055.194.94%3364.73430.465.671.91%1112.91167.554.674.68%3314.73382.067.281.99%1162.91222.159.294.85%3264.73330.966.191.99%1212.91279.266.345.19%3214.73289.474.692.27%1262.91327.564.634.87%3164.73234.569.772.16%1312.81382.669.825.05%3114.73188.974.172.33%1362.81433.670.814.94%3064.73137.372.572.31%1412.71485.072.234.86%3014.73088.974.112.40%1462.71536.173.354.78%2964.73040.175.402.48%1512.71591.979.144.97%2914.72989.274.482.49%1562.71645.182.425.01%2864.72947.782.982.82%1612.71697.484.674.99%2814.72891.176.362.64%1662.71749.987.164.98%2764.72845.480.622.83%1712.71796.683.894.67%2714.72794.079.312.84%1762.71848.886.084.66%2664.72747.582.723.01%1812.71900.087.324.60%2614.72695.881.073.01%1862.71952.790.014.61%2564.72644.579.783.02%1912.71999.286.484.33%2514.72598.183.333.21%1962.72052.789.974.38%2464.72547.182.383.23%2012.72103.290.504.30%2414.72500.285.473.42%2062.72148.485.723.99%2364.72451.987.113.55%2112.72202.389.544.07%2314.72402.988.143.67%2162.72252.790.014.00%2264.72349.784.973.62%2212.72303.791.013.95%2214.72298.884.113.66%2262.72352.289.443.80%2164.72256.091.294.05%2312.72402.790.023.75%2114.72204.489.664.07%2362.72452.289.513.65%2064.72156.691.834.26%2412.72497.985.173.41%2014.72106.791.974.37%2462.72548.785.933.37%1964.72053.188.354.30%2512.72599.686.833.34%1914.72005.690.814.53%2562.72644.781.963.10%1864.71953.588.784.54%2612.72696.383.553.10%1814.71905.190.394.74%2662.72748.886.113.13%1764.71852.888.024.75%2712.72798.986.143.08%1714.71801.586.774.82%2762.72847.284.492.97%1664.71752.787.995.02%2812.72895.983.132.87%1614.71699.384.554.98%2862.72942.179.342.70%1564.71645.981.184.93%2912.72992.679.842.67%1514.71602.087.275.45%2962.73044.281.482.68%1464.71544.179.355.14%3012.73089.676.912.49%1414.71492.677.895.22%3062.73140.878.102.49%1364.71444.179.415.50%3112.73189.776.992.41%1314.71392.377.565.57%3162.73234.071.322.21%1264.71333.468.625.15%3212.73284.171.392.17%1214.71280.565.785.14%3262.73334.671.882.16%1164.71230.665.915.36%3312.73379.066.291.96%1114.71175.160.375.14%3362.73423.560.821.78%1064.71124.759.975.33%3412.73473.260.461.74%1014.71067.853.104.97%3462.73517.754.931.56%964.71011.847.054.65%3512.73569.556.731.59%3514.73573.258.491.64%由上述的程序还可以描绘出倾斜时储油量的误差随罐内油位高度的变化情况，如图4-11所示：图4- 11由上述图我们可以知道倾斜时不管是进油后储油量的误差还是出油后储油量的误差均随罐内油位高度的变化呈正态分布，说明该模型建立得贴合实际。由上述得到储油罐发生变位时体积关于的公式我们可以给出罐体变位后油位高度间隔为1cm的罐容表标定值（即进/出油量与罐内油位高度的表格），如表4-3所示：表4- 3罐体变位后油位高度间隔为1cm的罐容表标定值高度（h）cm体积（v）L高度（h）cm体积（v）L高度（h）cm体积（v）L高度（h）cm体积（v）L13.53126461.5511413.9762491.326.26352749452145677253439.974828527.1531498.4782576.6414.756329560.9541540.9792619.1520.690830595.2551583.5802661.4627.854231630.1561626.3812703.6736.316332665.6571669.2822745.5846.142433701.5581712.2832787.2957.393534738591755.3842828.71070.12735774.9601798.58528701184.396836812.2611841.8862911.112100.254137850621885.1872951.813117.747538888.2631928.5882992.314136.92339926.7641971.9893032.515157.824640965.7652015.4903072.416180.493941100.5662058.891311217204.971421044.6672102.3923151.218228.9431084.5682145.7933190.119254.94411248692189.1943228.620281.9451165.3702232.5953266.721309.8461206.2712275.8963304.422338.5471247.2722319.1973341.723368.1481.2886732362.3983378.524398.5491330.1742405.4993414.925429.7501371.9752448.41003450.74.2问题二我们一开始考虑当储油罐水平放置时油量与浮油子高度的函数关系。当储油罐水平放置无横向与纵向倾斜时，设油高为h，所建立的空间直角建立坐标系如图4-12所示：图4- 12储油罐中储油的体积可分为两部分，圆柱体中储存油的体积与两端的球罐体中油的体积，根据对称性，两端的球罐体中油的体积相等。用垂直于z轴 的平面截储油罐，根据截面法求体积公式可求出储油的体积。第一种情况：1、对于圆柱体部分：2、，如图4-13所示：图4- 13图4- 14由图4-13可看出，根据几何关系：用平面截球罐体，设截得小圆半径为，所图4-14所示，根据几何关系得：截面圆的油高为，积分解得：第二种情况： 对于圆柱体的部分同第一种情况对于球罐体体积截面圆的油高为这种情况与第一种情况相同，因此与第一种情况相同。模型的建立：V与h的函数关系为：我们就油所占的体积进行初步讨论,建立笛卡尔坐标系如图4-15,可以得出球冠体的球面方程为:图4- 15如图4-15所示，油罐体中所含油的部分分成两份，分别计算它们的体积。因为右半部分的体积利用投影法可以由以下约束公式算出由上述第一个公式可以得出曲顶柱体的盖面方程为左半部分的体积近似利用截面法可以由以下公式算出，由图4-15可以得出点的坐标为；过垂直于轴做截面，设过的截面面积为；由图4-15可知过的截面面积为近似可表示如下：有上述的可以得到储油罐中油量与油位高度及变位参数的函数关系：以下为在固定的横向偏转角度条件下，就纵向倾斜角度的变化进行分类讨论。令情形当时，（1）该情况的图4-16所示：图4- 16由图4-16可知 和式代入式，得：由此可得：（2）该情况如图4-17所示：图4- 17同上可得（3）该情况如图4-18所示：图4- 18,设实际储油罐总体积为，情形当时，（1）该情况的图示如下：图4- 19（2）该情况如图4-20所示：图4- 20（3）该情况如图4-21所示：图4- 21,情形当时，（1）该情况如图4-22所示：图4- 22（2）该情况如图4-23所示：图4- 235模型的进一步讨论就第二问我们在上述讨论一般情况分析的基础上对左半部分进行了详细的分析。以下方程用来求：首先我们求出截面和球冠体的交线：6模型的优缺点6.1优点（1）对于问题一，从储油罐无变位与变位两个方面建立了模型，求得的储油量理论值与储油量实际值吻合得较好，相对误差小；特别是水平时相对误差稳定，符合要求，因此能准确的预测出罐体变位后油位高度间隔为1cm的罐容表标定值。（2）对于问题二，先就一般情况考虑，再细分成8种情形，思维缜密，考虑周到。利用Matlab软件求出 值，求得各高度下的储油量理论值与储油量实际值吻合得较好，相对误差小，符合要求，因此能准确的预测出罐体变位后油位高度间隔为10cm的罐容表标定值。（3）模型实用性好，具有很强的现实应用指导意义，数学推导严谨，理论性强。6.2缺点模型的假设不够充分，对误差的检验不够详细，没有从多方面进行误差诊断。7参考文献1韩中庚.数学建模竞赛-获奖论文精选与点评M.北京：科学出社,2007.2王文波.数学建模及其基础知识详解M.武汉：武汉大学出版社，2006.3王庚，王敏生.现代数学建模方法M.北京:科学出版社,2006.4赵静,但琦.数学建模与数学实验M.北京:高等教育出版社,2007.5姜启源,谢金星,叶俊.数学模型M.北京:高等教育出版社,2003.6冯杰,黄力伟,王勤,尹成义.数学建模原理与案例M.北京:科学出版社,2007.7附录7.1附件表7-1无变位进油流水号油罐号累加进油量/L油位高度/mm采集时间11150159.022010-08-18 10:32:18121100176.142010-08-18 10:33:18131150192.592010-08-18 10:34:18141200208.52010-08-18 10:35:18151250223.932010-08-18