储油罐的变位识别与罐容表标定数模优秀论文_毕业论文

2010高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。

如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：云南大学参赛队员(打印并签名) ：1. 洪建武2. 张艳3. 刘继萍指导教师或指导教师组负责人(打印并签名)：李海燕日期： 2010 年 9 月 10 日2010高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：储油罐的变位识别与罐容表标定摘要通常加油站的地下储油罐都是通过预先标定的罐容表进行实时计算，以便得到罐内油位高度与储油量的变化情况，然而，储油罐在使用一定时间后，罐体会因地基变形等原因发生变位，因此需定期对罐容表进行重新标定。

本文针对这一情况，建立了储油罐体积积分模型，综合运用立体几何、微积分、数据拟合、MATLAB 编程、EXCEL 表格等知识，研究解决了储油罐的变位识别与罐容表标定问题。

针对问题一，我们用了数学中的积分原理和几何知识建立了图像处理模型】【1。

在对图像处理模型改进的基础上建立了体积积分模型]2[、理想模型]3[。

对模型进行了合理的理论证明和推导，得出了理想状态下变位前后罐内油量和油位高度的函数关系式，储油罐变位前的函数关系为()()()⎰--==hdyb y ba aL h S h v 022220L 2储油罐变位后的函数关系为：()()()()⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎨⎧≤<+-⎪⎭⎫⎝⎛--=≤<+-⎪⎭⎫⎝⎛--⎪⎭⎫ ⎝⎛<≤-⎪⎭⎫⎝⎛=⎰⎰⎰3220103211000100212/cos cos )(2/cos tan cos S cos )0()2/cos(cos S h h h v dh h h S ab v h h h v dh L h h S h h dh h h v h h h h h απαπαπααααπα然后借助于MATLAB 软件，对附件中所提供的数据进行处理，并分别绘制变位前后实测值、理论值的图像，编程得出罐体变位后油位高度间隔为1cm 的罐容表标定值。

通过比较图变位前后实测值、理论值的图像，可以得出变位后的关系图像相对于变位前的曲线变化缓慢，即变位后据预先标定的罐容表所测得的油位高度高于实际油位高度，等同于罐容表的刻度标记值偏大，应在保持分度值不变的情况下加大刻度间距或是在保持间距不变的情况下减小分度值。

针对问题二，针对问题2我们同样运用了数学中的积分原理和几何知识建立了图像处理模型。

在对图像处理模型改进的基础上建立了体积积分模型、理想状态模型。

对模型进行了合理的理论证明和推导，得出了理想状态下变位后罐内油量和油位高度的函数关系式，然后借助于MATLAB 软件，对附件中所提供的数据进行处理。

关键词：储油罐、罐容表、体积积分模型、理想模型、液面高度一问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。

许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。

按照有关规定，需要定期对罐容表进行重新标定。

图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。

图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。

请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。

（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为α=4.1°的纵向变位两种情况做了实验，实验数据如附件1所示。

请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。

（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度α和横向偏转角度β）之间的一般关系。

请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm 的罐容表标定值。

进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。

二模型假设1 题中所给储油罐的尺寸均为内壁尺寸；2 外界因素如温度、湿度、外力等不会造成储油罐几何形状的改变；3 忽略储油罐内其它设备占用的体积；4 题中所给数据均为精确值；5 储油罐密封性良好；6 所给数据均在储油罐内液面稳定的状态下测定；7 忽略每次进油和出油时滞留在进油管和出油管内的油量；8 储油罐无变位（倾斜变位）进油后罐内油位高度终值即为无变位（倾斜变位）出油前罐内油位高度的初始值；三符号及变量说明四问题分析问题一本问中需研究储油罐罐体变位后对罐容表的影响并得到罐体变位后油位高度间隔为1cm的罐容表标定值。

求解本问的步骤为：第一步：分别建立变位前后的积分模型导出储油量与油位高度的函数关系式；第二步：针对变位前、变位后两种情况，使用MATLAB 绘制理论所得函数关系图像，同时，据附件1数据绘制储油量与油位高度的关系图；第三步：针对变位前后的函数关系图象进行比较，由此得罐体变位后对罐容表的影响；第四步：由相应的模型可得出罐体变位后油位高度间隔为1cm 的罐容表标定值。

问题二此问解决的主要是对参数的标定，首先解决的是求出相应罐内储油量与油位高度及变位参数之间的关系，然后使用使用曲线拟合求出相应的参数，根据得出的函数关系可求出并给出罐体变位后油位高度间隔为10cm 的罐容表标定值。

主要问题在于积分的求解，由符号积分不能得出。

五 模型建立及求解问题一的模型建立与求解 1 变位前的函数关系模型两端平头的椭圆柱体截面示意图如图1，正面示意图如图2。

图1 小椭圆油罐截面示意图图2 小椭圆油罐正面示意图建立坐标系如图1所示。

储油罐内油品体积是油品液面高度的函数：V=f(h) 由椭圆的标准方程得：()()12222b y ba a x --=用极限法求体积，则()()()()2L 2022220⎰--==hdyb y ba a L h S h v通过MATLAB 软件绘制几何解析的方法所得的储油罐内油量与油位高度系的图形以及附件1所给实际测量值的储油罐内油量与油位高度关系的图形。

图3 变位前罐内油量与油位高度的关系图2 变位后的函数关系模型建立空间直角坐标系，绘出小椭圆型储油罐三维立体图：图4 变位后罐体示意图图5 椭圆阴影部分示意图如图4所示，将罐体高度分为三段来求出油量与相应油位高度的函数关系，由积分及相应几何知识可得()()()5;4cos 23sin 21321h h h b h L h +===；；αα由图5所示阴影部分面积为()()()()()()()8)0()(7)2/cos(cos S 62S 101101022220h h dh h S h h h S dyb y ba a y hyv <≤=-⎪⎭⎫ ⎝⎛=--=⎰⎰则体积为第一部分απα()()()()()()()11)()(10)(92/cos tan cos S cos 10221220021v h v h v h h h dhh S h v L h h S h S hh +=≤<=-⎪⎭⎫⎝⎛--⎪⎭⎫ ⎝⎛=⎰απααα第二部分()()()()()()()14132/cos 12cos 2033133132v h v h v h h h dhS v h h S ab S h h +=≤<-=⎪⎭⎫⎝⎛--=⎰απαπ第三部分则可得函数()h v 如下：()()()()()⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎨⎧≤<+-⎪⎭⎫⎝⎛--=≤<+-⎪⎭⎫ ⎝⎛--⎪⎭⎫ ⎝⎛<≤-⎪⎭⎫⎝⎛=⎰⎰⎰3220103211000100212/cos cos 15)(2/cos tan cos S cos )0()2/cos(cos S h h h v dh h h S ab v h h h v dh L h h S h h dh h h v h h h h h απαπαπααααπα对以上使用MATLAB 编程可以很方便的得出相应图如下：图6 变位后罐内油量与油位高度的关系示意图3 分析得出罐体变位后对罐容表的影响由于储油罐的纵向倾斜，使得储油罐与水平面之间形成一定的倾斜角，导致据预先标定的罐容表所得的油位值与实际油位高度有一定偏差，但该偏差在允许的范围内，则据几何知识可得实际油位高度与通过变位前预先标定的罐容表所得油位高度的关系。

()()()()()17cos /16cos cos )2tan 4(cos 000αααααh v H V h H h b h h h h H =+=∴+≤≤-=观察图5、图6可知，理论值与实测值存在一定偏差，但均在误差允许的范围内；通过比较图5、图6两条曲线，图6所示曲线相对于图5曲线变化缓慢，即变位后据预先标定的罐容表所测得的油位高度高于实际油位高度，等同于罐容表的刻度标记值偏大，应在保持分度值不变的情况下加大刻度间距或是在保持间距不变的情况下减小分度值。

4 罐体变位后油位高度间隔为1cm 的罐容表标定值综合运用MATLAB 、EXCEL 软件编程得出罐体变位后与刻度相对应的油量值，继而得出罐体变位后油位高度间隔为1cm 的罐容表标定值。

表2 罐体变位后油位高度间隔为1cm的罐容表标定值问题二的模型建立与求解图7 实际储油罐变位后示意图图8实际储油罐部分示意图问题2中将题中所示储油罐正面示意图理想化为图7中所示立体，将由倾斜引起球冠体切线的变化暂不考虑。

该立体的横向转动对相应油面高度的油量的求解没有影响，但会导致油位探针有所偏转而引起预先标定的罐容表有所偏差。

图9实际储油罐截面示意图阴影部分面积为()()18arcsin,r d rdr d r S ⋅-=图10 实际储油罐变位后部分示意图由此图10可知当纵向偏转角大于θ时，所建横轴与圆弦相切；当纵向偏转角大于θ时，所建横轴与圆弦相交，故分为θα<与θα≥两种情况来考虑。