清华大学姜启源jiangqy@tsinghuaeducn.ppt
清华大学 姜启源 jiangqytsinghua.edu.cn,储油罐的变位识别与罐容表标定,全国大学生数学建模竞赛 2010年A题,题目 题目的背景与分析 问题（1）的通常做法 问题（2）的通常做法 一种建立储油量模型的新方法 对学生论文的评述 题目的启示,通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。,题 目,请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 （2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度 ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。,附件1 实验数据,附件2 检测数据,来自于加油站设备研究与生产企业的一个实际课题.,（2）根据实际检测数据，正确识别罐体是如何变位的，具体变了多少？同时要给出罐容表的修正标定方法和结果，属于“反问题”。,题目背景,问题分析,问题由两部分组成： （1）为了观察检验罐体变位对罐容表的影响，在已知变位参数的情况下，检测出油位高度和油量的对应数值，建模分析罐容表的变化规律，并给出修正的罐容表，属于“正问题”。,问题（1）的通常做法,根据附件1 的实验数据画出罐体不变位（=0）和罐体变位（=4.10）时油位高度h和储油量V的曲线.,实验数据表明，在同样的油位高度h下，罐体变位时储油量V减少.,z=H下的面积,油位高度h时坐标y处z=H,油位高度h时坐标y处的截面面积,罐体变位时油位高度h和储油量(体积)V的数学模型,化重积分为定积分,高度h与储油量V的模型,高度h与储油量V的模型 V(h,),=4.10的罐容表(部分),按照模型讨论变位对罐容表的影响,与实验数据表示的影响一致.,变位后储油量平均约小200L.,模型结果与实验数据的比较,=4.10,=0,模型结果与实验数据基本吻合，而在同样的油位高度h下，储油量v的实验数据均比模型结果小.,可以有各种解释.,是否要修正、怎样修正模型？,问题（2）的通常做法,1. 罐体变位,时油位高度h和储油量V的数学模型,整体考虑,写出左右球面、圆柱面及油位高度h的油平面方程.,将罐体分为3个区域：D1,D2,D3,在3个区域中分别写出体积的重积分，并化为（对z的）定积分.,将罐体分为圆柱体和球缺，圆柱体按照问题（1）考虑,z=H下的面积,截面面积S (y,h,),圆柱体体积V1 (h,),将罐体分为圆柱体和球缺，球缺部分单独考虑,精确计算球缺部分体积V2,将倾斜角的油平面近似为水平面（=0）计算球缺部分体积V2.,直接将罐体的球缺部分按照体积不变条件折合成圆柱，一起计算.,用各种方法得到罐体体积V (h,) ,一般是积分表达式, 给定, 可以数值计算油位高度h时的储油量V.,问题（2）的通常做法,2. 根据储油量的数学模型V (h,)和实测数据（附件2）辨识罐体变位参数,Vi =Vi (hi, 0,0)不能用！,应该用数据ui , hi与模型V (h,) 辨识参数,确定,应使储油量的改变Vi 与出油量ui 相吻合！,辨识准则,确定,使Q最小,搜索法：按照, 的可能范围划分区间，逐步搜索.,直接利用MATLAB的非线性曲线拟合程序 lsqcurvefit,辨识方法,辨识结果,ui , hi取进油前的300组数据,=2.110, =4.310的修正罐容表(部分),罐体变位后显示储油量与实际储油量的比较,模型及结果检验,取进油后的300组数据检验辨识结果, , 的灵敏性检验,的灵敏度,的灵敏度,纵向倾斜角比横向偏转角 对罐容量的影响显著得多！,一种建立模型V (h,) 的新方法,绕x轴转动,绕z轴转动,T2T1= T1T2 ?,油浮子坐标 (0, h-r, L/2-d),过油浮子的油平面方程,坐标变换,油位高度h的储油区域为D,油位高度h时的储油量V,定义示性函数,利用MATLAB的三重积分程序计算储油量,triplequad(x,y,z)(-x*cos(a)*sin(b)+y*cos(a)* cos(b)+z*sin(a)=-s-sqrt(R2-x.2-y.2).*(z=s+sqrt(R2-x.2-y.2),-r,r,-r, r,-(s+R),s+R),triplequad (function, x1,x2, y1,y2, z1,z2),被积函数 function 在长方体中的三重积分:,用于function 是示性函数 I(x,y,z)的情况:,计算时间较长；计算精度受积分域大小的影响.,在上述坐标旋转下罐体纵向倾斜角为负(正负无关).,1）罐体变位后储油量的计算（重积分）： 两端有油和一端有油区域的划分； 积分限的确定； 积分顺序的选择。,对学生论文缺陷的评述,2）用所谓“等效高度”：得出不变位时高度与体积的关系V(h)；变位后根据(主)截面面积不变得到等效高度与原高度的关系h' =f(h)；将它代入原式V(h) 得到变位后高度与体积的关系V' (h') 。,3）辨识变位参数所用的数据和准则不对：用附件2中的 “显示油量”与模型在“显示油高”下的计算油量作拟合。,如果模型正确，用这种错误方法得到的变位参数应该是=0, =0 。,反映对该问题的理解根本错误,4）含混地叙述：“用数据和最小二乘得到 = , = ”,“显示油量”是罐体未变位时的计量值.,5）对问题（1）用实验数据修正模型投入的精力过多（解释误差、用多项式拟合误差、用比例系数修正模型等）.,6）对问题（1）的要求“罐体变位后对罐容表的影响”未作讨论.,7）对问题（2）的参数辨识不做检验和灵敏性分析 。,8）问题（1）,（2）未按要求给出罐容表.,9）从文献中引用工程算法、近似公式的问题.,10）用随机模拟(Monte-Carlo)方法计算重积分的问题.,题目的启示,微积分在工程实际问题中是有用的！,熟练地使用软件对有效、快速地求解模型是重要的！,仔细审题，准确、完整地理解题目的含义和要求是做好论文的首要的和关键的一步！,
