储油罐的变位识别与罐容表标定数模优秀论文.doc

储油罐的变位识别与罐容表标定摘要通常加油站的地下储油罐都是通过预先标定的罐容表进行实时计算，以便得到罐内油位高度与储油量的变化情况，然而，储油罐在使用一定时间后，罐体会因地基变形等原因发生变位，因此需定期对罐容表进行重新标定。本文针对这一情况，建立了储油罐体积积分模型，综合运用立体几何、微积分、数据拟合、matlab编程、excel表格等知识，研究解决了储油罐的变位识别与罐容表标定问题。针对问题一，我们用了数学中的积分原理和几何知识建立了图像处理模型。在对图像处理模型改进的基础上建立了体积积分模型、理想模型。对模型进行了合理的理论证明和推导，得出了理想状态下变位前后罐内油量和油位高度的函数关系式，储油罐变位前的函数关系为储油罐变位后的函数关系为：然后借助于matlab软件，对附件中所提供的数据进行处理，并分别绘制变位前后实测值、理论值的图像，编程得出罐体变位后油位高度间隔为1cm的罐容表标定值。通过比较图变位前后实测值、理论值的图像，可以得出变位后的关系图像相对于变位前的曲线变化缓慢，即变位后据预先标定的罐容表所测得的油位高度高于实际油位高度，等同于罐容表的刻度标记值偏大，应在保持分度值不变的情况下加大刻度间距或是在保持间距不变的情况下减小分度值。针对问题二，针对问题2我们同样运用了数学中的积分原理和几何知识建立了图像处理模型。在对图像处理模型改进的基础上建立了体积积分模型、理想状态模型。对模型进行了合理的理论证明和推导，得出了理想状态下变位后罐内油量和油位高度的函数关系式，然后借助于matlab软件，对附件中所提供的数据进行处理。关键词：储油罐、罐容表、体积积分模型、理想模型、液面高度19一 问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为a=4.1的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。二 模型假设1 题中所给储油罐的尺寸均为内壁尺寸；2 外界因素如温度、湿度、外力等不会造成储油罐几何形状的改变；3 忽略储油罐内其它设备占用的体积；4 题中所给数据均为精确值；5 储油罐密封性良好；6 所给数据均在储油罐内液面稳定的状态下测定；7 忽略每次进油和出油时滞留在进油管和出油管内的油量；8 储油罐无变位（倾斜变位）进油后罐内油位高度终值即为无变位（倾斜变位）出油前罐内油位高度的初始值；三 符号及变量说明表1 符号说明符号说明单位a两端平头的椭圆柱体截面椭圆的长半轴分米b两端平头的椭圆柱体截面椭圆的短半轴分米v小椭圆型储油罐储油罐内油的体积升l小椭圆型储油罐柱体的长度分米h储油罐内油位的实际高度分米so小椭圆油罐截面示意图阴影部分的面积分米探针与储油罐内壁接触点距水平面的垂直距离分米小椭圆储油罐向上倾斜一端椭圆面离水平面最近的一点与水平面的垂直距离分米小椭圆储油罐向下倾斜一端椭圆面离水平面最远的一点与水平面的垂直距离分米小椭圆储油罐向上倾斜一端椭圆面离水平面最远的一点与水平面的垂直距离分米纵向倾斜角度小椭圆油罐变位后第一部分的体积立方分米小椭圆油罐变位后第二部分的体积立方分米小椭圆油罐变位后第三部分的体积立方分米h据油位探测装置所得的油位高度分米实际储油罐向上倾斜一端椭圆面离水平面最近的一点与水平面的垂直距离分米实际储油罐向下倾斜一端椭圆面离水平面最远的一点与水平面的垂直距离分米实际储油罐向上倾斜一端椭圆面离水平面最远的一点与水平面的垂直距离分米球体的半径分米球冠体最大横截面的半径分米r任意横截圆的半径分米四 问题分析问题一本问中需研究储油罐罐体变位后对罐容表的影响并得到罐体变位后油位高度间隔为1cm的罐容表标定值。求解本问的步骤为：第一步：分别建立变位前后的积分模型导出储油量与油位高度的函数关系式；第二步：针对变位前、变位后两种情况，使用matlab绘制理论所得函数关系图像，同时，据附件1数据绘制储油量与油位高度的关系图；第三步：针对变位前后的函数关系图象进行比较，由此得罐体变位后对罐容表的影响；第四步：由相应的模型可得出罐体变位后油位高度间隔为1cm的罐容表标定值。问题二此问解决的主要是对参数的标定，首先解决的是求出相应罐内储油量与油位高度及变位参数之间的关系，然后使用使用曲线拟合求出相应的参数，根据得出的函数关系可求出并给出罐体变位后油位高度间隔为10cm的罐容表标定值。主要问题在于积分的求解，由符号积分不能得出。五 模型建立及求解问题一的模型建立与求解1 变位前的函数关系模型两端平头的椭圆柱体截面示意图如图1，正面示意图如图2。yxso图1 小椭圆油罐截面示意图lv图2 小椭圆油罐正面示意图建立坐标系如图1所示。储油罐内油品体积是油品液面高度的函数：v=f(h)由椭圆的标准方程得： 用极限法求体积，则 通过matlab软件绘制几何解析的方法所得的储油罐内油量与油位高度系的图形以及附件1所给实际测量值的储油罐内油量与油位高度关系的图形。图3 变位前罐内油量与油位高度的关系图2 变位后的函数关系模型建立空间直角坐标系，绘出小椭圆型储油罐三维立体图：2b=1.2m2a=1.78m油位探针h1lozxh3h2yh0图4 变位后罐体示意图yxso图5 椭圆阴影部分示意图如图4所示，将罐体高度分为三段来求出油量与相应油位高度的函数关系，由积分及相应几何知识可得 由图5所示阴影部分面积为 则可得函数如下：对以上使用matlab编程可以很方便的得出相应图如下：图6 变位后罐内油量与油位高度的关系示意图3 分析得出罐体变位后对罐容表的影响由于储油罐的纵向倾斜，使得储油罐与水平面之间形成一定的倾斜角，导致据预先标定的罐容表所得的油位值与实际油位高度有一定偏差，但该偏差在允许的范围内，则据几何知识可得实际油位高度与通过变位前预先标定的罐容表所得油位高度的关系。 观察图5、图6可知，理论值与实测值存在一定偏差，但均在误差允许的范围内；通过比较图5、图6两条曲线，图6所示曲线相对于图5曲线变化缓慢，即变位后据预先标定的罐容表所测得的油位高度高于实际油位高度，等同于罐容表的刻度标记值偏大，应在保持分度值不变的情况下加大刻度间距或是在保持间距不变的情况下减小分度值。4 罐体变位后油位高度间隔为1cm的罐容表标定值综合运用matlab、excel软件编程得出罐体变位后与刻度相对应的油量值，继而得出罐体变位后油位高度间隔为1cm的罐容表标定值。表2 罐体变位后油位高度间隔为1cm的罐容表标定值刻度/dm油量/l01.67440.13.5310.26.26350.39.97480.414.7560.520.6910.627.8540.736.3160.846.1420.957.39410.93747.2113776.611.13805.311.2383311.33859.811.43885.611.53910.311.63933.911.73956.111.83976.711.93995.5124012.7问题二的模型建立与求解图7 实际储油罐变位后示意图1mr03m图8实际储油罐部分示意图问题2中将题中所示储油罐正面示意图理想化为图7中所示立体，将由倾斜引起球冠体切线的变化暂不考虑。该立体的横向转动对相应油面高度的油量的求解没有影响，但会导致油位探针有所偏转而引起预先标定的罐容表有所偏差。 图9实际储油罐截面示意图阴影部分面积为 图10 实际储油罐变位后部分示意图由此图10可知当纵向偏转角大于时，所建横轴与圆弦相切；当纵向偏转角大于时，所建横轴与圆弦相交，故分为与两种情况来考虑。由已知可得当时：图11 纵向偏转角小于时示意图参考零水平以下高度为相应函数第二部分 油罐两头的总体积为一个球冠体的体积 第三部分 当时：图12 纵向偏转角大于时示意图相应得函数 首部分 第二部分 油罐两头的总体积为一个球冠体的体积第三部分 由以上各式可得实际油位的高度与油量的函数，因考虑为较理想的状况，图13 实际储油罐的正视图由图可得，六 模型的应用与推广本文采用的储油罐体积积分模型有效的模拟出储油罐油量与油位高度以及变位参数的关系，具有较强的现实意义，可在现实生活中广泛应用于计算各种油罐的储油量，同时，该积分模型运用的微积分思想可推广到求解各种不规则形状物体的面积、体积。七 模型的评价与改进1模型的优点该模型综合运用微积分原理、数据拟合、绘图比较，基于matlab、excel等软件，较为精确地模拟出储油罐油量与油位高度以及变位参数的关系。此外，该模型建立在储油罐为理想状态下，即忽略了罐体内各种设备占用的体积以及各种可能导致储油罐几何形体发生变化的因素，使得模型简化、易于计算。2模型的缺点由于该模型将储油罐考虑为理想状态，忽略了罐体内各种设备占用的体积、各种可能导致储油罐几何形体发生变化的因素以及储油罐罐壁厚度，使得通过积分模型得到的储油量与油位高度的关系与真实值存在一定误差。3 模型的改进针对模型存在的缺点，在使用题中所给储油罐尺寸进行计算时应在所给尺寸的基础上减去储油罐罐壁相应的厚度；在使用积分方法求解储油罐内储油量时应除去在液面下的油位探针、进油管以及出油管占用的体积，加上进油管、出油管内残留的油品体积。参考文献1 （美）肖温格特（schowengerdt,r.a.），微波成像技术国家重点实验室(译)，遥感图像处理模型与方法（第三版），出版社：电子工业出版社，2010。2 同济大学应用数学系，高等数学（第五版 下册），出版社：高等教育出版社，1978。3 罗柏森，浅述“物理学中的理性模型”，/coursefile/rexue_20090509/?action=content&m=6514&a=8055&todo=show，访问时间：2010年9月12日。4 管冀年，赵海，卧式储油罐罐内油品体积标定的实用方法，20045 谭永基，数学模型，出版社：复旦大学出版社，1997，1：1110。6 宗容，施继红，尉洪，李海燕，数学实验与数学建模，出版社：云南大学出版社，2009， 18。7 王研玲，李明，椭圆形式封头卧式贮罐液与容积对应关系的建立，齐齐哈尔大学学报j，20028 张志涌，杨祖樱，matlab教程，出版社：北京航空航天大学出版社，2006,177191。9周建明,项忠权,储液罐翘离影响的近似计算a, 北京:地震出版社,1990 , 140 161附录附录一源程序代码小椭圆储油罐变位前理论值与实际值图像syms h ya = 17.8/2;b = 12/2;l = 24.5;v = 2*int(sqrt(a2-a2/b2*(y-b)2),0,h)*l; v1 = zeros(1,78);h1 = zeros(1,78);v1 = xlsread(data1.xls,c2:c79);v1 = v1 + 262;h1 = xlsread(data1.xls,d2:d79); h = 0:0.01:12;vv = double(subs(v, h);plot(h, vv)hold onplot(h1/100, v1, *)hold offxlabel(h),ylabel(v)legend(,)小椭圆储油罐变位后理论值与实际值图像syms y h;a = 1.78/2*10;b = 1.2/2*10;l = 24.5;h1 = l*sin(4.1*pi/180);h2 = 2*b*cos(4.1*pi/180);h3 = h1 + h2; s0 = 2*int(sqrt(a2-a2/b2*(y-b)2),y,0,y);s1 = subs(s0,y,h/cos(4.1*pi/180);v1 = int(s1/cos(pi/2-4.1*pi/180),h,0,h); %0=hh1hh1=0:0.01:h1;v11 = double(subs(v1,hh1); s2 = subs(s0,y,h/cos(4.1*pi/180) - subs(s0,y,h/cos(4.1*pi/180) - l*tan(4.1*pi/180);v2 = int(s2/cos(pi/2-4.1*pi/180),h,h1,h); %h1=hh2hh2 = h1:0.01:h2;v22 = double(subs(v2,hh2) + max(v11); s3 = pi*a*b - subs(s0,y,(h-h1)/cos(4.1*pi/180);v3 = int(s3/cos(pi/2-4.1*pi/180),h,h2,h); hh3 = h2:0.01:h3;v33 = double(subs(v3,hh3) + max(v22); v1 = zeros(1,53);h1 = zeros(1,53);v1 = xlsread(data1.xls,c2:c54);v1 = v1 + 215;h1 = xlsread(data1.xls,d2:d54);plot(h1/100,v1,*)hold onplot(hh1,v11),plot(hh2,v22),plot(hh3,v33) legend(,)xlabel(h),ylabel(v(h)小椭圆储油罐变位后罐体变位后油位高度间隔为1cm的罐容表标定值%syms y h;a = 1.78/2*10;b = 1.2/2*10;l = 24.5;ph = 4.1*pi/180;h1 = l*sin(ph);h2 = 2*b*cos(ph);h3 = h1 + h2; s0 = 2*int(sqrt(a2-a2/b2*(y-b)2),y,0,y);s1 = subs(s0,y,h/cos(ph);v1 = int(s1/cos(pi/2-ph),h,0,h); %0=hh1v10 = double(subs(v1,h,h1); s2 = subs(s0,y,h/cos(ph) - subs(s0,y,h/cos(ph) - l*tan(ph);v2 = int(s2/cos(pi/2-ph),h,h1,h); %h1=hh2v20 = double(subs(v2,h,h2); s3 = pi*a*b - subs(s0,y,(h-h1)/cos(ph);v3 = int(s3/cos(pi/2-ph),h,h2,h); %h2=h=0 & hh(i)=h1 & hh(i)=h2 & hh(i)=h3 vv(i) = double(subs(v3,h,hh(i) + v10 + v20; endend 附录二刻度表刻度/dm油量/l0167440.13.5310.26.26350.39.97480.414.7560.520.6910.627.8540.736.3160.846.1420.957.394170.1271.184.3971.2100.251.3117.751.4136.921.5157.821.6180.261.72041.8228.911.9254.882281.862.1309.762.2338.542.3368.142.4398.532.5429.662.6461.492.74942.8527.142.9560.93595.253.1630.153.2665.583.3701.533.4737.963