罐容表 储油罐的变位识别与罐容表标定分析

储油罐的变位识别与罐容表标定分析摘要：为了实现对储油罐的变位识别与罐容表标定的分析，建立了罐内储油量基于油位高度及变位参数的数学模型，过程如下：第一：由题目可知，油位探针与罐体固定，以油位探针与罐底的交点为坐标原点，建立空间直角坐标系o-xyz。由微元法可得无变位时的储油量与油位的函数表达式为：。代入题目所给数据并用Matlab编程进行检验，最大误差为4.42%，平均误差为3.31%；当罐体纵向倾斜变位=时，相当于罐体以o为原点绕z轴旋转，由坐标旋转公式导出小椭圆型储油罐表面方程，油面为水平面，油标杆与油面夹角为，分五种情况利用微元法分别求出了储油量与油位的表达式）。 并依据所建数学模型用Matlab编程得出罐体变位后油位高度间隔1cm的罐容的标定值，具体数据见表3所示。经检验最大误差为6.62%，平均误差为3.47%。第二：油罐若再横向旋转角度，此时旋转轴为油罐中轴线。由于油罐由园柱体和球冠围成，关于中轴线对称，旋转后油罐表面方程保持不表，油面高度保持不变，只是油浮子旋转了角度，由o点变到，故只需在计算中让油罐纵向变化，利用问题一的方法算出储油量与原油标高度的表达式，再由式和进行换算。得出罐体变位后标定罐容表的数学模型和）。不妨设间的理论关系为：，令，用MATLAB导入题目中所给实际储油罐的检测数据，计算得出： ， 。利用所建模型得出罐体变位后油位高度间隔10cm的罐容的标定值，具体数据见表4所示。通过比较标定值和实际采样数据的拟合曲线，得到均方差很小。关键字：空间直角坐标系 微积分 Matlab 误差 标定值一：问题的背景近年来，随着中国国民经济的快速发展、交通基础设施的不断改善和机动车保有量的快速增加，加油站已经成为民众生活中不可或缺的一部分。通过查资料可知[1]，在储油罐的安装过程中和长期使用后，难免会发生倾斜、变形等问题。对于加油站埋地油罐的倾斜、变形等变化在现场技术条件下难以准确测定，这势必造成储存罐容积表的不准确，尤其是近年来随着中石化收购加油站数量的增加，大量系统外的加油站纳入中石化系统管理，这些加油站又多为个体加油站，大部分没有油罐容积表，即使有也很不准确。如何妥善地解决这一问题就成为现在数量管理工作的一大难题。二：问题的重述请结合问题背景中的资料，现知道通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化，从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。请用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐，分别对罐体无变位和倾斜角为?=的纵向变位两种情况做了实验，实验数据已经给出。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。对于实际的储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据，根据所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用实际检测数据来分析检验所建模型的正确性与方法的可靠性。三：基本假设1：假设题目中所给数据真实有效；2：假设地下储油罐一直处于地下恒温的环境中，即汽油不会发生热胀冷缩；3：假设储油罐的罐壁厚度忽略不计；4：假设油位探针、注油口和出油管在罐体内的体积忽略不计；5：假设油位探针进口孔很小，游浮子不会沿着探针浮出；四：模型符号说明罐体纵向变位的倾斜角度罐体横向变位的倾斜角度椭圆型储油罐横截面长半轴椭圆型储油罐横截面短半轴油浮子测得的实际油位高度储油罐的容积储油罐第个区域的容积第个区域，油位所能达到的最大水平高度过油浮子的水平面的水平高度所得切面的面积=时，油标所能达到的最大油位高度储油罐两端球冠体的半径发生横向变位时油罐整体下移高度发生纵向变位时的油位高度发生横向变位后油浮子的水平高度发生横向变位后油浮子所能测得的最大水平高度五：问题的分析1．问题的分析：为了掌握罐体变位后对罐容表的影响，利用小椭圆型储油罐分别对罐体无变位和倾斜角已知的纵向变位做了实验，并给出了实验数据。考虑到油罐的倾斜角已知且油位探针与油罐的的位置固定，以油位探针与罐底的交点为坐标原点，以油罐发生变位前油罐底边为坐标轴，建立空间直角坐标系O-xyz。如下图所示：为得到较准确的结果，需要衡量坐标系建立的合理性，可通过当前建立的坐标求出油罐无变位时油位高度与储油量的理论对应关系，得出理论数据，再结合题目所提供的实验数据相对比，求出相对误差。若误差在允许的范围内，则说明该坐标系的建立可行，可进行后续问题的计算。分析油罐纵向倾斜时，题目倾斜角已知，考虑其是以O为定点，在xOy平面内发生偏转。为得到其储油量与油位高度的关系，必须涉及体积计算问题，考虑用平行于xoz的平面切割储油罐，累加切面面积来计算。由于所取切面的高度不同，所得切面形状也有多种可能，这里可分为五种情况来讨论，具体见模型建立与求解。从而建立起罐体纵向倾斜后储油量基于油位高度计算的数学模型。利用程序计算出罐体变位后油位高度间隔1cm的罐容表标定值，并制excel表将其表示出来。2．问题的分析：实际储油罐的两端为球冠体，可考虑利用问题中所建的空间直角坐标系进行求解。由于油位探针与罐体的位置固定，可先分析罐子横向偏转倾斜对油位高度的单独作用影响，再将两者结合起来，从而得到变位后油罐内储油量与油位高度间的数学模型。由于在问题中得出了小椭圆形储油罐纵向变位后的影响，这里仅仅是形状有些差异，可继续利用此思想。考虑先找出储油罐横向变位前后油位变化关系，然后通过类似问题的解决方法得出油罐未发生横向变位时储油量与油位高度间的关系式，再利用油罐横向变位前后油位变化关系，得出变位后油罐储油量与油位高度及变为参数间的一般关系。在此条件下，可以利用题目所给实际储油罐的检测数据，依靠模型用Matlab编程估计变位参数的取值，得出储油量与油位高度间的具体数学模型，从而可得到罐体变位后油位高度间隔为10cm的罐容表的标定值，并制excel表将其表示出来。再进一步利用所给实际检测数据分析检验模型的正确性与方法的可靠性。可通过计算误差来进行衡量，若误差较大，则需考虑限制条件是否有遗漏，若在允许范围内，则可通过对理想因素的实际化来增加模型的精度，提高可靠性。六：模型的建立及求解1．问题： 根据分析，以油位探针与罐底的交点为坐标原点，建立空间直角坐标系O-xyz。I:当罐体无变位时，罐体在坐标平面内的截面图如下图1所示： SHAPE 参考所查资料[2]并运用所学数学知识[3]得此时小椭圆型储油罐的空间函数为：又题目已知： a=1.78/2=0.89m b=1.2/2=0.6m m=0.4m n=2.05m 用y=切割椭圆柱体，切面是一个关于x轴对称的矩形，积分时只积一半，结果乘以2，其面积为： 将其面积在y轴上积分可得其总体积，即为给定油位时的储油量： = 代入数据，并通过matlab编程可得与h的理论与实际对应值表如表1所示：理论累加进油量实际累加进油量油位高度理论累加出油量实际累加出油量油位高度55.6871159.0252.7253.31681150.7255.6871176.14102.72103.30291123.991196.3619462.622702.722702.9413426.803656.913803.91041152.363652.723652.8032160.483706.913882.32841193.493702.723702.79142.62将题中所给无变位进油和无变位出油的数据代入式进行检验，得出其最大误差为4.42%，平均误差为=3.31%，两误差均较小。 II:当罐体发生纵向倾斜变位时，相当于罐体以o为原点绕z轴逆时针旋转，故罐体在空间平面内的截面图如下图2所示：根据所建坐标系，罐体发生=的位变时，相当于罐体在xoy平面内转动了的角度，设罐体上的一点变位前后的坐标分别为：A,。则有：发生变位后的小椭圆型储油罐在空间坐标系中的空间函数为： 由式可得： 设油位高度为h，过油浮子水平面与y轴的交点为,考虑到油位高度的取值不同，切面的形状不一样，将椭球罐体分为五个区域，具体划分如图2所示，记油位能达到的最大值为:=2b=1.2m，各个区域的取值范围可根据椭圆柱面和两端平面与xoy平面的四条交线的四个交点的y坐标及三角函数关系得到：下面用定积分的元素法来计算罐体内的储油量。用的平面对椭圆型柱体进行切割，把罐体分为五个区域，现对这五个情况进行分别讨论：：当时，，即=0。即罐体内储油量，但油标高度，当=0时，=，此时，此时为区域一的情况，罐体内的储油量为： ：当时，即时，此时为区域二的情况，罐体内的储油量为：。 = 当时，，此时的。：当此时为区域三的情况，罐体内的储油量为：。 = ：当，此时为区域四的情况，罐体内的储油量为：。 = ：当。此时为区域五的情况，油标达到最大值，说明纵向倾斜较多过大，计量系统已经不能准确的确定罐内的储油量，但可以得到它的一个范围为。联立式子-，利用matlab编程可得：=1.625L =33.8L =2730.2L =1222.97L此时罐体总容积在五个区域内的总表达式为： 上式即为建立的研究罐体变位后对罐容表的影响的数学模型。代入数据，通过matlab编程可得发生变位后与h的理论与实际对应值表和根据所建模型可得罐体变位后油位高度间隔为1cm的罐容标定值分别如表2和3所示： 理论累加进油量实际累加进油量油位高度理论累加出油量实际累加出油量油位高度747.86782.6401309411.2953.11225361020.65797.86833.240565423.45103.07285681007.731797.731899.306347670.221102.32383762.093297.733310.5828911034.242501.298477425.833299.743313.9780671035.362551.255482411.731.625123.4694716.15371988.1483433.4341.791606136.0364756.12172051.3433482.2242.292435149.4497797.39052114.8593529.4083.130506163.7415839.95262178.6253574.9134.310846178.9439883.7982242.5663618.6725.840461195.0886928.91362306.6053660.6157.728313212.2068975.28332370.6643700.689.985289230.32921022.8882434.6653738.80412.62415249.48571071.7062498.5253774.92815.6595269.70551121.7112562.1653808.99619.10771291.01671172.8772625.5013840.95422.98686313.4461225.1712688.4513870.75327.31668337.01921278.562750.9293898.34532.11846361.76051333.0052812.8533923.68637.41499387.69261388.4682874.1373946.73843.23044414.83661444.9052934.6983967.46349.59027443.21161502.2692994.4513985.82856.52117472.83491560.5133053.3114001.80464.05092503.72191619.5833111.1964015.36672.20825535.88571679.4273168.0224026.49181.0228569.33711739.9863223.7074035.16390.52495604.08481801.2013278.174041.367100.7457640.13481863.013331.3314045.094111.7166677.49081925.3483383.1114046.336根据误差公式得到此时的最大误差为6.26%，平均误差为3.47%。2．问题：根据第二问的分析，油罐若再横向旋转角度，此时旋转轴为油罐中轴线，由于油罐由柱体和球冠围城关于中轴线对称，旋转后油罐表面方程保持不表，油面高度保持不变，只是油浮子旋转了角度，有o点变到故只需在计算中让油罐只纵向变化，利用问题一的方法算出原油标高度，再进行换算。未发生变位时，以油位探针与油罐的交点O为原心建立空间直角坐标系O-xyz，油罐两端球冠体的半径为R，-3m-2m3/2mm根据空间平面几何关系可得： 未发生变位时的油罐实体方程为： 现分析横向倾斜后油位的变化情况：过油位探针作油罐的横截面图，如图4所示：由于油罐对称，容易观察出经过旋转前后切面油面深度始终都为，旋转后油位刻度为h，可以找出旋转前后油位刻度和h的关系：由图可知： 联立式-得： 又知 下面再单独研究纵向倾斜后储油量和油位高度的关系。｛1｝当时，其截面图如图5所示： 当发生纵向倾斜变位后，设罐体上的一点在变位前后的坐标分别为：，根据矩阵可得前后点的坐标的关系式为： 将式代入可得纵向旋转后，油罐实体的方程为： 记为实际测得的油位高度，为油面的水平高度，依据图示的几何关系可得： ，，， 用的平面对球冠体进行切割，记切面上的左球冠在x轴上的坐标为，右球冠在x轴上的坐标为，通过计算可得： 记切面左球冠所在曲线函数为，切面右球冠所在曲线函数为,柱体切面所在曲线的函数为，通过旋转后的油罐实体方程可知： 同问题，将罐体分为五个区域，现对这五个情况进行分别讨论： I：1、当时，即 当时，,此时 II：当时，即，此时。 III: 。 IV：当。 V：当,此时 此时为区域五的情况，油标达到最大值，计量系统已经不能准确的确定罐内的储油量，但可以得到它的一个范围。联立式子-可得罐内储油量与油位高度及变位参数的关系式为： ｛2｝当时，其截面图如图6所示：3m6m1m2m1mαyxz图6 由图可知： ， 同样分五种情况进行讨论，具体过程如下： 此时得到的罐内储油量与油位高度及变位参数的关系式为式和即为建立的研究罐体变位后对罐容表的影响的数学模型。由以上分析建立了间的理论关系为： 但是在实际的应用过程中，很多理想处理的因素都不能忽略，一般情况下试验得到均不能满式，令 上式可知这是一个关于的二元方程，要使I最小，即满足： 联立方程-，用MATLAB导入题目中所给实际储油罐的检测数据，计算得出 ： ， 将变位参数的数值代入油罐储油量与油位高度的关系模型中，用MATLAB编程得到油罐变位后油位高度间隔为10cm的罐容标定值，具体数据如表4所示：利用Matlab编程，根据所建模型可预测出变位参数的数值，以及罐体变位后油位高度h间隔为10cm的罐容标定值V，具体数据如表4所示： hhh436.2721187.0648759.991969.623926.0451253.953502.9326685.0553620.725036.2629512.6855912.456695.3732322.2357983.368793.4235146.4759987.6911004.5438096.2261957.9613412.8840795.0963908.3215906.5643534.665848.6718487.7846185.5967779.26模型正确性与方法可靠性的检验：通过比较标定值和实际采样数据的拟合曲线，得到均方差很小。拟合曲线如图所示：七：模型的评价及推广1．模型的优点： 利用油位探针与罐体的位置固定，以探针与罐底的交点为坐标零点建立的空间直角坐标系，先计算出无变位时的罐内油位高度与储油量的函数关系，用题中所给数据检验坐标系建立的合理性，误差不大，则保证了模型的可行性。在发生横向偏转倾斜和纵向倾斜变位时，可根据点在空间的坐标变换关系得出变位后球体在空间的表达式，方便积分计算其体积。 根据所建的模型可以很快得出已知变位情况和油位高度条件下的罐内容积表，模型的适应性强、推广性好、应用性广。2．模型的缺点：在模型的建立过程中，对影响罐内容积的微小因素进行了理想的处理，但是在实际应用中应考虑这些因素对模型准确性的影响，如油的热胀冷缩和罐内的设备等对油位的高度有所影响。代入题目所给数据后进行验证时所得的误差虽然在允许的范围内，但可通过将理想处理的影响因素实际化，以减小误差，使得模型的精确度更高。3．模型的推广：该模型不仅适用于加油站的地下储油罐的变位识别和罐容表标定的分析，对根据液面高度计算卧式容器内的物料容积都有很好的指导意义，推广性好。参考文献：[1]：董文辉 冯东，《储油罐容积表自动编制系统的构成与运用》，《石油库与加油站》，第11卷 第4期：第29-30页，2002年8月出版。[2]：王宇 陈志军 苏建伟 郭昕，《椭圆形封头卧式容器容积的计算》，《河南化工》，第27卷 第2期：第75-76页，2010年2月出版。[3]：同济大学应用数学系，《高等数学及其应用》，北京市西城区德外大街4号：高等教育出版社，2004年11月第1版。4：张德丰，《MATLAB数值计算方法》，北京市百万庄大街22号：机械工业出版社，2010年1月第1版。5：John H.Mathews 、Kurtis D.Fink，《数值方法第四版》，北京市海淀区万寿路：电子工业出版社，2005年12月第1版。clear;H=[159.02176.14192.59208.50223.93238.97253.66268.04282.16296.03309.69323.15336.44349.57362.56375.42388.16400.79413.32425.76438.12450.40462.62474.78486.89498.95510.97522.95534.90546.82558.72570.61582.48594.35606.22618.09629.96641.85653.75665.67677.63678.542053.832103.832105.062155.062205.062255.062305.062355.062404.982406.832456.832506.832556.832606.832656.832706.832756.832806.832856.832906.832906.912956.913006.913056.913106.913156.913206.913256.913306.913356.913406.913456.913506.913556.913606.913656.913706.91r=R/1000;A=1.78-2*t;B=1.2-2*t;a=A/2;b=B/2;l=2.45-2*t;h=H/1000;for i=1:1:78 f=-1; c=asin; g=.*h; d=sqrt; v=.*.*d+b.*b.*c+./2)); esp1=abs)./))function [ql,Ak,xk]=guasslegendre% 高斯-勒让德数值积分% 输入参数有效性检验if nargin==1 a=-1;b=1;n=7;tol=1e-8;elseif nargin==3 n=7;tol=1e-8;elseif nargin==4 tol=1e-8;elseif nargin==2|nargin 5 error^,n+1))/);tk=roots; % 求积节点% 计算求积系数Ak=zeros;for i=1:n+1 xkt=tk; xkt=[]; pn=poly; fp=@polyval/polyval); Ak=quadl; % 求积系数% 积分变量代换，将[a,b]变换到[-1,1]xk=/2*tk+/2;% 检验积分函数fun有效性fun=fcnchk*/2;% 计算积分值ql=sum;h=[411.29423.45438.33450.54463.90477.74489.37502.56514.69526.84538.88551.96564.40576.56588.74599.56611.62623.44635.58646.28658.59670.22680.63693.03704.67716.45727.66739.39750.90761.55773.43785.39796.04808.27820.80832.80844.47856.29867.60880.06892.92904.34917.34929.90941.42954.60968.09980.14992.411006.341019.071034.241035.36v1=[747.86797.86847.86897.86947.86997.861047.861097.791147.791197.731247.731297.731347.731397.731447.731497.731547.731597.731647.731697.731747.731797.731847.731897.731947.731997.732047.732097.732147.732197.732247.732297.732347.732397.732447.732497.732547.732597.732647.732697.732747.732797.732847.732897.732947.732997.733047.733097.733147.733197.733247.733297.733299.74v2=/1000;h=h/1000;for i=1:1:78 %v =4.11-1.0.*.*.^+0.534.*asin+0.838).*.^2./.^2+6.00.*./.^2+4.90./.^2-1.0.*.^2./.^2.*-1.0.*./.^2.*.^2+0.333.*.^3./.^2); v=2.05+0.727.*).^.*h-0.436.*).^+1.31.*asin; esp11=abs.^2)./.^2)); esp22=abs)./));vn=.*1000;espp1=sum./78espp2=sum./78function [max_x,max_v]=findmxformat long;eps=1.e-3;max_v=0;max_x=a;fa=subs,findsym),a); %左断点函数fb=subs,findsym),b); %右断点函数step=n/5;tol=1;tmp=0;while tol eps t=a; for j=0: %找出最大值x t=a+j*step*/n; pt=fa+*/; ft=subs,findsym),t); if abs max_v max_v=abs; max_t=t; tol=abs; tmp=max_x; step=step/2;format short;clear;digits;A0=1.78;B0=1.2;a=A0/2;b=B0/2;m=0.4;n=2.05;A=1.2;aerf=*pi;h1=-m*sin;h2=n*sin;h3=a*cos-m*sin;h4=a*cos;bb=b*b;k=2*a/b;syms Sx h x ySx=sym+b-y*cos)^2))xjs1=y/tan;xjx1=)/cos;xjs2=y/tan;xjx2=)/cos;xjs3=)/cos;xjx3=)/cos;xjs4=)/cos;xjx4=-2*b)/sin;digits;v1=k.*simple)v2=k.*simple)v3=k.*simple)v4=k.*simple)yjs1=0;yjx1=h1;yjs2=h2/cos;yjx2=0;yjs3=h3/cos;yjx3=h2/cos;yjs4=h4/cos;yjx4=h3/sin;v1=simplify)v2=simplify)v3=simplify)v4=simplify)function f=lspoly1n=length;b=zeros;f=zeros;for k=1:m+1 f=h .^;a=f *f;b=f *v;c=ab;f=flipud;function[I,step]=combine Traprlif eps=1.0e-4;h=/2;I1=0;I2=,findsym),a)+subs,findsym),b))/h;while abs eps n=n+1; h+/n; I1=I2; I2=0; for i=0;n-1 x=a+h*i; x1=x+h; I2=I2+*,findsym ),x)+subs,findsym),x1) I=I2; step=n;function df=Midpointif nargin==2 h=0.1:else if disp return;y1=subs,findsym ),x0+h);y2=subs,findsym ),x0-h);df=/;

“原创力文档”前称为“文档投稿赚钱网”，本站为“文档C2C交易模式”，即用户上传的文档直接卖给用户，本站只是中间服务平台，本站所有文档下载所得的收益归上传人所有【成交的100%】。原创力文档是网络服务平台方，若您的权利被侵害，侵权客服QQ：3005833200 电话：19940600175 欢迎举报，上传者QQ群:784321556