储油罐的变位识别与罐容表标定.doc

储油罐的变位识别与罐容表标定摘要本文是关于储油罐的变位识别与罐容表标定的文章。加油站通常都有若干个储存燃油的地下储油罐，许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化，从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。本文利用微积分的知识，建立相关模型，通过matlab编程求解，给出了重新标定罐容表的具体方法。对于第一题，在无变位的情况下，建立合适的空间直角坐标系，算出横截面积与油位高度的关系，然后乘上桶的长度，即可得到储油量与油位高度的关系，这样就得到了无变位情况下的理想罐容表。在倾斜变位的情况下，同样地，建立空间直角坐标系，以桶的底边母线为Y轴，分三段算出横截面积与油位高度和Y轴坐标的关系，然后利用matlab通过积分的方法编程计算倾斜变位情况下的理想罐容表。通过计算可以得到无变位和倾斜变位的情况下理论值的平均误差分别为3.3%和3.9%，可以直接将无变位和倾斜变位的情况下的理想罐容表比较。从图表可以看出倾斜变位的情况下理想曲线与实际存在误差，由理想值总是大于实际值可以得出误差为系统误差，通过matlab实现最小二乘拟合，得到更加精确的罐体变位后油位高度间隔为1cm的罐容表标定值。对于第二题，将实际储油罐分成两个部分进行求解，即两端的球罐体和中间的圆柱体。对于球罐体中的液体体积计算可以假设油液面在球冠体内是平齐的，高即为圆柱体和球冠体连接处的液面高，由于少算了左端球罐体内的油量而多算了右端球冠体内的油量，左右两个球罐体内的计算误差基本可以抵消，这样就可以很好的简化模型的计算，而中间的圆柱体内油量的计算可以参照第一小题的计算方法，这里不再赘诉。将、分别从1取到10，通过matlab软件编程计算出最优的、值，确定最符合实际值的参数以后即得到罐容表函数，然后以10厘米为间隔，给出罐容表标定值，再与进出油的实际值比较误差。最后利用附件二中的数据对模型的可靠性进行了检验，检验结果表明模型较为合理。 关键词 微积分 matlab 拟合 最小二乘法一、 问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。题目给出了一种典型的储油罐，主体为圆柱体，两端为球冠体。还分别给出罐体纵向倾斜变位的示意图和横向偏转变位的截面示意图。请用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用小椭圆型储油罐（如图，两端平头的椭圆柱体），分别对罐体无变位和倾斜角为a=4.1度的纵向变位两种情况做了实验，实验数据由附件1给出。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。(b) 小椭圆油罐截面示意图 油油浮子出油管油位探针注油口水平线2.05mcm0.4m1.2m1.2m1.78m(a) 小椭圆油罐正面示意图图4 小椭圆型油罐形状及尺寸示意图（2）对于实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。并且进一步利用附件2中的实际检测数据来分析检验模型的正确性与方法的可靠性。二、 问题分析第一题：要研究罐体变位后对罐容表的影响，需要将无变位和变位情况下的罐容表作比较，因此需要作出两张罐容表。在无变位的情况下，以椭圆桶左边椭圆的最低点为原点，母线方向为Y轴方向，与母线垂直向上为Z轴方向，垂直纸面向外为X轴方向建立空间直角坐标系。利用积分算出横截面积与油位高度的关系，然后乘上桶的长度，即可得到储油量与油位高度的函数关系，这样就可以得到无变位情况下的理想罐容表。在倾斜变位的情况下，以同样的方法建立空间直角坐标系，如图：此时的横截面积与油位高度和Y坐标有关，利用双重积分计算出储油量与油位高度的函数关系（积分时要油位高度要分三种情况），得到倾斜变位情况下的理想罐容表。通过matlab将两条曲线画在同一张图表上进行比较即可得到变位对罐容表的影响。 通过比较计算得出的变位理想罐容表和附件中给出的实际数值，发现存在一定的误差，可以通过matlab实现最小二乘拟合，给出更加符合实际的罐容表曲线，从而得到更加精确的变位后油位高度间隔为1cm的罐容表标定值。第二题： 要得到实际储油罐变位后的罐容表，难点在于同时考虑纵向倾斜和横向偏转，还有就是储油罐两端的球冠体。首先用和第一题一样的方法建立空间直角坐标系，然后将纵向倾斜和横向偏转分开考虑，横向偏转对计算造成的影响相对不大，只需通过一定的函数将有横向偏转的油位高度等效转化为无横向偏转的油位高度；对于纵向倾斜造成的影响，计算时我们将油的体积分为两部分，中间圆柱体内为一部分，两端球冠体内为另一部分。圆柱体内的一部分可以利用第一题的结论，将椭圆方程改为圆方程即可求解。两端球冠体内的部分，由于油液面在球冠体内是倾斜的，这里运用近似的方法，假设油液面在球冠体内是平齐的，高即为圆柱体和球冠体连接处的液面高，这样虽然产生了一定误差，但是由于少算了左端球罐体内的油量而多算了右端球冠体内的油量，总体上不会产生大的偏差，对结果影响不大。要得到最优化的、的解决方法是使确定参数的罐容表与实际值的误差最小。令、分别从1取到10，则共有100种组合，根据附件中给出的出油实际值，在每一种组合中算出301对油位高度和储油量的数据，分别求每一组的301对理论值与实际值的差的绝对值之和，得出偏差最小的一种组合，然后再在该组合的周围利用同样的方法精确计算最优解，精确到0.1。确定最符合实际值的参数以后即得到罐容表函数，然后以10厘米为间隔，给出罐容表标定值，再与进出油的实际值比较误差。三、 模型假设1.假设油浮子相对于油液面理想为质点；2.假设对油罐进行进出油测数据时管道内没有残留；3.假设所有数据都是在相同的外界条件下测量；4.假设油浮子到达最高点不再进油，到达最低点不再出油；5.假设储油罐的倾斜变位总是左低右高；6.假设实际储油罐的纵向倾斜角不超过10；四、 符号说明h 油浮子显示的油位高度(dm)v1 小椭圆型储油罐无变位情况下的储油量v2 小椭圆型储油罐倾斜变位情况下的储油量 实际储油罐变位的纵向倾斜角度 实际储油罐变位的横向偏转角度h 有横向偏转的油位高度等效转化后的无横向偏转的油位高度五、 模型的建立与求解第一题：在无变位的情况下，储油量可以表示为：在倾斜变位情况下，储油量的表达式为： 计算得到无变位和倾斜变位的情况下理论值的平均误差分别约为3.3%和3.9%，误差很小，可以直接将无变位和倾斜变位的情况下的理想罐容表比较，在matlab上作出罐容表曲线：通过比较两条曲线可以得到变位后对罐容表的影响：油液面很少，到达浮游子之前，油位高度显示始终为0；变位后的第一阶段内（），曲线斜率小于变位前，储油量变化较慢；变位后的第二阶段内（），曲线的斜率和无变位情况下曲线的斜率差不多，储油量增长幅度相同，两者的差值基本不变；第三阶段（），曲线斜率趋缓，当油浮子高度为12分米时，油箱还未加满，但根据假设，此时不再进油。将变位情况下的理论罐容表曲线和实际值在图上做出（如下图），可以看出两者有很小的误差，并且总是实际值小于理论值，可以解释为罐内一些管道占据了很小一部分体积，这是系统误差，可以通过计算将其缩小。 通过matlab实现最小二乘拟合，可以将系统误差减小，由下图给出的变位情况下拟合后罐容表曲线与实际值可以看出在第二阶段内，实际值与罐容表曲线符合的很好。另外由于当油浮子高度很小时，仍然使用拟合的数据会产生很大的误差。不符合实际，所以h在0到0.5分米的范围内，罐容表使用理论值。变位后油位高度间隔为1cm的罐容表标定值如下表：油位高度(dm)储油量(L)油位高度(dm)储油量(L)油位高度(dm)储油量(L)油位高度(dm)储油量(L)01.67443.1577.50176.21802.9139.33115.6450.13.5313.2611.39976.31845.9639.43155.0190.26.26353.3645.84686.41889.0869.53194.0390.39.97483.4680.82066.51932.2739.63232.6880.414.75633.5716.29986.61975.5099.73270.9460.520.69083.6752.26396.72018.7839.83308.7930.626.248783.7788.69326.82062.0839.93346.2090.732.203253.8825.56866.92105.39610.03383.1720.839.560553.9862.87197.02148.70910.13419.660.948.381674.0900.58517.12192.01210.23455.6491.058.723954.1938.6917.22235.29110.33491.1141.170.641464.2977.1737.32278.53410.43526.0291.284.185354.31016.0157.42321.72810.53560.3661.399.404134.41055.27.52364.86110.63594.0961.4116.34394.51094.7137.62407.9210.73627.1871.5135.04234.61134.547.72450.89310.83659.6051.6155.32494.71174.6657.82493.76610.93691.3131.7176.9464.81215.0737.92536.52711.03722.2711.8199.77284.91255.758.02579.16211.13752.4331.9223.70955.01296.6838.12621.65911.23781.7512.0248.67965.11337.8578.22664.00311.33810.1662.1274.61885.21379.2598.32706.18111.43837.6142.2301.47175.31420.8758.42748.17911.53864.0142.3329.18945.41462.6918.52789.98311.63889.2672.4357.72795.51504.6968.62831.57911.73913.232.5387.04765.61546.8768.72872.95211.83935.6332.6417.11195.71589.2178.82914.08711.93956.3582.7447.88725.81631.7088.92954.96812.03975.4472.8479.34225.91674.3369.02995.5812.9511.44776.01717.0879.13035.9093.0544.17626.11759.9519.23075.936第二题：首先将有横向偏转的油位高度等效转化为无横向偏转的油位高度：中间圆柱体内的储油量va与h的函数关系： 左端球罐体内vb与h的函数关系： 右端球罐体内vc与h的函数关系： 将va,vb,vc相加，即可得到带有参数、的罐容表函数关系：求、的最优解：目标函数：通过matlab使用遍历搜索法求的最优解：以10厘米为间隔的罐容表标定值：油位高度（dm）储油量（L）油位高度（dm）储油量（L）油位高度（dm）储油量（L）1581.60561122316.742153877.0521389.1561225382.712256828.1132694.6741328514.532359669.4644415.0431431694.582462380.1256417.8531534905.672564937.0968655.6771638130.882667314.46711093.261741353.422769481.71813701.811844556.472871400.62916456.451947723.052973017.381019334.932050835.853068863.53六、 模型评价与推广优点： 本文计算储油量均采用了微