储油罐的变位识别与罐容表标定（定稿）

储油罐的变位识别与罐容表标定摘要通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，通过预先标定好的罐容表，可得到罐内油位高度与储油量的变化关系。但许多储油罐使用一段时间以后，由于地基变形等原因，使罐体的位置发生纵向倾斜和横向偏转，从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。因而建立储油罐变位后储油量与油高及变位参数（纵向倾斜和横向偏转）之间的一般关系，对罐体储油量的真实计算及加油站的经营管理具有重要意义。对于问题一，本文先建立没有变位时的罐体储油量和油位高度的关系，将计算值与实际值进行比较，进行图形仿真和误差分析，从而检验模型的可靠性和准确性。对于发生纵向倾斜后的椭圆型储油罐，在油液面低于柱体右端最低点和高于左端最高点，及两者之间，储油量与油位高度有不同的关系式，因而我们分了三段积分处理，得出储油量与油位高度的函数关系式。用建立好的函数关系式计算出给定油位探针监测高度的储油量，和实际储油量进行图像曲线对比，并进行误差分析，从而验证建立的函数关系式的准确性。在用建立好的模型对变位和未变位的两种情况的储油量随探针监测油位高度变化的曲线进行对比并列表分析，从而得出罐体变位后同一监测高度，变位后罐容体的实际储油量比原先罐容表上标定的值小，并计算出罐体变位后油位高度间隔1CM的罐容表标定值。对于问题二，本文利用几何关系，将横向偏转修正，以消除其对储油量的影响，将问题归结为只需要计算纵向偏转对储油量的影响，将储油量的计算分成三部分圆柱体和左右球冠体，圆柱体可直接积分得到，球冠体通过柱面坐标变换，将二重积分转换为定积分，然后利用微分中值定理近似计算该定积分。三者相加得到整个储油量体积，且和问题一一样分为油液面低于圆柱体部分右端最低点和高于左端最高点，及两者之间三段，再整合为一个函数关系式。得出的计算值与实际数据比较，进行误差分析，从而用线性拟合的方法对函数关系式进行修正使其与实际值的误差更小。最后利用循环迭代并结合矩形套定理，逐步缩小范围，以确定偏转角，和以使误差在一定精度范围内符合实际值，最后将得到的偏转角，218代入建立的函数关系式，用以模拟检验，得出结果与实际相符。之后我2504们给出了油位高度间隔为10CM的罐容表标定值。最后，本文对模型进行了进一步的讨论和改进，对问题二建议制定出不同的对应储油量体积增长的拐点的表，只要根据实际数据利用二阶差分近似求得拐点和位置，只需查表即可得到。和关键字罐容体储油量分段积分微分中值定理线性拟合循环迭代一、问题的背景通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。二、问题的提出与重述由于地基变形等原因，使罐体的位置发生变位，从而导致罐容表不能显示实际的储油量，罐容表误差过大而不能正常使用，造成加油站油品虚假盈亏。这样，就给加油站经营管理带来一些问题。如造成加油站虚假盈亏，无法对油品数量进行正确的监控和管理，以及年底盘底或新旧站长变更时，无法进行正常的油品库存交接。因而需要对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为的纵向变位两种情况做了实验，41实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1CM的罐容表标定值。（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10CM的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。图1储油罐正面示意图图2储油罐纵向倾斜变位后示意图油图3储油罐截面示意图（B）横向偏转倾斜后正截面图地平线地平线垂直线油位探针（A）无偏转倾斜的正截面图油位探针油位探测装置地平线油3M油图3储油罐界面示意图图4小椭圆油罐截面示意图三、基本假设1假设油浮子始终处于水平状态，并且油浮子的体积不记，视为质点；2假设油位探针是固定的，不发生任何转动；3忽略外界因素对储油罐内部的影响并排除储油罐的机械故障；4假设油位不受温度、压力等因素的影响；5假设油位探针检测液位控制灵敏，罐容表标定无误；6假设储油罐容器壁光滑平整，没有凹凸现象；7忽略储油罐内各器件所占的体积。四、模型的主要符号变量说明问题一的主要符号说明V无变位时椭圆形储油罐的储油量；V1椭圆型储油罐变位后，底部部分覆盖时的储油量；V2椭圆型储油罐变位后，底部覆盖顶部未覆盖时的储油量；V3椭圆型储油罐变位后，底部覆盖顶部部分覆盖时的储油量；椭圆型储油罐的总体积；总探针监测到的油位高度；椭圆型储油罐变位后，底部部分覆盖时，油位与底部的相交线离左端的距离；0Z椭圆型储油罐变位后，顶部部分覆盖时，油位与顶部的相交线离左端的距离；1问题二的主要符号说明储油罐总体积V油罐圆柱体部分的体积1左端球冠体体积2右端球冠体体积3油位探针监测到的高度0H横向偏转修正后的高度五、问题的分析题目中的第一问要求我们建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1CM的罐容表标定值。我们首先建立了无变位时小椭圆形储油量与油位高度的一般关系函数式，并用附件1中的数据检验模型的正确性，由此得到无变位情况下理论值与实际值的相对误差A。由于小椭圆型储油罐纵向变形后，在油量少到低于油位探针最底端时和油量多于探针与椭圆柱体顶部的交点时，不能写出储油量与油位高度对应函数关系式，对这两种情况不做出具体的对应关系式。纵向变位后，在油液面低于柱体右端最低点和高于左端最高点，及两者之间，储油量与油位高度有不同的关系式，因而我们分了三段处理，得出储油量与油位高度的函数关系式。将油量理论值和实际值比较而得到相对误差B，并和误差A进行比较，检验建立的函数关系式的正确性，并用所得的理论计算公式对变位前后储油量同一高度储油量进行比较，从而得出罐体变位后对罐容表的影响以及油位高度间隔为1CM的罐容表标定值。第二问要求我们对图1所示的罐体建立变位后标定罐容表的数学模型，得出罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度）之间的一般关系式，并根据附件2的实际检测数据，用所建立的模型确定变位参数。由于横向偏转不会引起油液面的变化，只会影响油位探针测得的油位高度，所以我们对纵横向变位后的标定高度转换成只有纵向变位时的油位标定高度，设为发生横向偏转后的油位探0H针测得的高度，为转换成只有纵向偏转的油位探针测得的高度，由于探针必经过探H针所在圆柱横截面的圆心，如图5，有，00COSCOSRRH所以把研究储油量与油位高度及参数，的函数关系转变为只研究储油量与油位高度及参数的关系。则研究方法和问题一类似，也要分为三段，油液面低于圆柱体部分右端最低点，油液面高于圆柱体部分左端最高点和介于两者之间的三段储油量和油位高度的关系式。其中对圆柱体部分和左右球冠体分别积分求油量体积，三部分油量体积相加得出三段储油量和油位高度的关系式，综合得出储油量与油位高度及变位参数的关系式，得出计算0,VH值，并与实际值比较，进行误差分析。然后运用线性拟合的方法对V进行修正，再利用附件2中的数据，用二分法原理，对，划一个比较宽的范围，结合闭矩形套定理编程，求出附件2的数据所对应的参数，从而得出关系式。由此可求出储油量的理论值，0H结合实际数据进行误差分析，验证函数关系式的可靠性，从而给出罐体变位后油位高度间隔为10CM的罐容表标定值。六、问题一的模型建立与求解由于对于特定的椭圆型储油罐，当其所处状态变位或未变位确定时，对于进油和出油的研究都一样，所以本问题只用进油这一情况进行分析。研究对于图4的小椭圆型储油罐，我们首先建立无变位时罐内储油量与油位高度的函数关系式。建立如图6所示的坐标系，设椭圆的长半轴为，短半轴为B,得A椭圆方程，则，设储油量的体积为V,椭21XYAB2XBY圆柱体的长度为L，油位高度为H，则22ARCSIN1VHB其中A178/2089M，B12/206M，L245M。把附件1工作表“无变位进油”中的油位高度一栏的高度值代入式，计算出储油量的理论值，计算MATLAB程序见附件一（PART1），实际值与理论计算值随油位高度变化的图像见图7。XXYXYAABB图6图5从图7可以看出实际值和理论计算值的曲线吻合的比较好。又将对应高度的实际值与计算值列入EXCEL表中并计算理论计算的储油量与实际储油量的差值，并算出理论计算的储油量对实际储油量的相对误差。其计算结果见附件二（SHEET1），现截取其中十行见下表一。从附件二SHEET1中的计算结果可以看出，对于同一油面高度，理论计算的储油量对实际储油量的相对误差，其相对误差的最大值理论计算值实际值实际值为34917316，最小值为3486559，总体平均误差为34883831,近似为3488，说明无变位情况下计算值和实际值的相差比例可以看成常数。图和表的结果说明了理论公式的科学性，同时也说明了积分求理论公式这种方法的合理性。从而得出了未变位时小椭圆罐的罐容表每隔1CM的标定值，程序见附件一（PART2），结果见附件二SHEET2，表二给出了罐容表的部分理论标定值。而计算所得的相对误差结果可以为椭圆型储油罐变位后罐容表的重新标定提供参考依据。00204060811214050010001500200025003000350040004500位位位位位位表一油位高度/MM理论计算储油量/L实际储油量/L理论值实际值/L（理论值实际值）/实际值159023228826312108826003488012817614374633362126330034897791925942636494121436490034866262208504781318462161318003491731622393529851951217851900348669922389758160585621960580034885765253666333526122135200348888892680468508166223081003486555928216736846871224846800348971912960378857777622657770034878871表二图7油位高度/M理论计算储油量/L油位高度/M理论计算储油量/L油位高度/M理论计算储油量/L油位高度/M理论计算储油量/L01163590411993107248915136598802450270516210008291084113946550380354062055070933066112411015椭球型储油罐发生倾角为410纵向变位后，在油液面低于柱体右端最低点和高于左端最高点，及两者之间，储油量与油位高度有不同的关系式，因而我们分了三段处理，得出储油量与油位高度的函数关系式（坐标系见图8）由以上数据可得，2178YX1在油液面低于柱体右端最低点时，设Z0为柱体底部部分覆盖时，油位与底部的相交线离柱体左端的距离，由题中图4可知Z004M时，油位探针不能检测到油，所以此时油位探针检测到的油位高度H值为H0。又柱体长度L245M，所以对于0404M时，油位探针不能监测到储油量的值，此处不予考虑，所以Z104M。因而有（为油位离柱体顶部的距离），则11TAN04TAN04YZ1Y2H对于Z处的油截面（为Z处油位离柱体顶部的距离），所以211TANTANYZ2Y可得Y的积分下限为21106T06768Z下限由上可得此时储油量的表达式1124506322450621781ZYZYVDZXDYDY下限下限总总式中为椭圆柱体的总体积，即V总308962451MABL总所以11224506378ZY036ZVDZD（）此处0Z104,H12,即11713MH12MTAN以上将椭圆柱体倾斜纵向变位时储油量与油位高度的关系分为三段进行考虑,得出了储油量与测得的油位高度的关系式、。而题目中附件1的工作表“倾斜变位进油”油位高度一栏油位高度的范围为，工作表4129MH0356“倾斜变位出油”油位高度一栏油位高度的范围为，所以72M对附件1中的数据只需用函数关系式（01469M24504TANTAN0620624504TANTAN060617821HZHZVDZXDYDYH11713M）。首先建立对应于式的体积积分函数VTUO（H,），程序见附件三。把附件一工作表“倾斜变位进油”油位高度一栏的数据代入公式计算，得出储油量的理论计算值，其计算的MATLAB程序见附件四。在加油过程中，理论值和实际值随高度变化的图像见图9。从图9可以看出理论计算值和实际值的吻合效果很好。同时，将对应高度的实际值与计算值列入EXCEL表中并并计算理论计算的储油量与实际储油量的差值，并算出理论计算的储油量对实际储油量的相对误差。其计算结果见附件五，现截取其中一部分，见下表三。从下表三和附件五知相对误差有一个波动范围，相对误差的最大值为4587221，最小值为12719502，平均误差值为32303292。此处的平均误差与椭圆型储油罐无变位时的平均误差值34883831很接近，这说明从整体上看，理论计算值与实际值产生的误差是一样的，即产生误差的原因是一样的。纵向变位后，从附件五的误差数值和油位高度可以看出误差数值随油位高度呈现出先升后降的趋势，这可能是由于椭圆型储油罐发生纵向变位后，沿Z轴各处液位深度不一致及油罐壁厚各处不完全均匀引起的。0405060708091115001000150020002500300035004000位位位位位位图9表三油位高度/MM理论计算储油量/L实际储油量/L理论计算值实际值（理论计算值实际值）/实际值4112999916499628636304900377052744234510474491012863458900341498334383311071651062864430480041684512450541156651111286437907003934969446391211263116286484033004162435789292303333296273706004002382950990434307876301273660302002191706591734312994930627367219100219474469299317884311273661103002123868894142322316531627360435300191085879546327323332127360503200188323339680933237243262736099420018694222980143368141331273554110016726688从以上对结果的分析可以知道，以上建立的函数关系式可以用来对罐体变位后储油量的计算。可用以上公式、对罐容体进行重新标定。在时，150,6H对变位和没有变位时同一高度的储油量比较的图10，MATLAB程见附件六。由图10可以看出变后油位探针监测同一油位高度时，变位前的储油量大于变位后的储油量。所以需要对罐容表进行重新标定。首先建立对应于式的体积积分函数V1TUO1（H,）（见附件七）和式的体积积分函数V3TUO3（H,）（见附件八），再分别编写程序算出对应的V10,1469H值，01469MH11713M对应的V2值（V2和附件三中的V等价），11713MH12M对应的V3的值，程序见附件九。从而得出罐体变位后油位高度间隔为1CM的罐容表标定值，标定值见表四。综合图形和表格可以看出，在误差允许的范围内，该模型计算得到的结果具有可行性，进而可以说明该模型是正确可靠的。从而问题一得到解决。表四油位高度/M储油量/L油位高度/M储油量/L油位高度/M储油量/L油位高度/M储油量/L0001674403059524520601798524090307242700135310316301462061184179709131120010026263503266558080621885131092315123400399748033701525606319285130933190110041475630347379584064197193109432286120052069080357748577065201537209532667220062785420368122030662058824096330442100736316303784997470672102275097334169100846142403888815370682145713098337851100957393503992672170692189126099341486201070127040965660807022325100345072011843968041100495407122758241013486064012100254104210445840722319086102352087013117747504310845350732362273103355511401413692304411247910742405372104358876901515781840451165336075244837210536218080161802591046120615507624912591063654201720399940471247234077253402107368591501822890660481288557078257664310837169180192548849049133011107926191151093747171020281857705013718810802661423110377663602130976080511413854081270355211138052660223385387052145601508227454911123833013023368142605314983520832787225113385981902439852850541540851084282874114388561802542965670551583499085287002211539103320264614906056162628308629110571163933859027493996705716691908729518311739560560285271438058171220808829923261183976510295609024059175532308930325311193995392120401259900204060811214050010001500200025003000350040004500位位位位位位位七、问题二的模型建立与求解1体积公式推导由以上问题分析式可得发生纵向偏转和横向偏转时的油位探针监测高度H0图10转换为只有纵向变位时的油位探针监测高度H，有00COSCOSHRHRH则以下只需建立只有纵向变位时储油量与油位高度的数学模型。纵向偏转下储油量与油位高度的关系式建立如下图11如图11（储油罐发生变位后的等价示意图）所示，表示油位探针与圆柱体左端1L面的距离，表示油位探针与圆柱体右端面的距离，H1，H2分别表示液面与圆柱体左端2L面、右端面的交线到圆柱体底部的距离，则有62TANL由图1和题目中的数据可得，R为球冠体的半径，R为圆柱体横21R截面的半径，由R15求得R1625M。图11中为油液面与左球冠体的交线和油液1XH面与圆柱体左端面的交线的高度差，设D为球冠体球心到对应圆柱体端面的距离，则DR10625M，由勾股定理，所以可以得到和22111RTANXX1XH的关系式。同理，可得到和的关系式2XH22TXXHDHRR我们把实际储油罐体积分为三部分计算，在问题的分析中已有说明，即为中间圆柱体、左球冠体、右球冠体体积的计算。（1）中间圆柱体体积的计算由图12所示的油面和圆柱体左右端面都相交的情况简单容易计算体积，所以先计算这种情况下的体积。如下所示。其他较复杂的情况油面低于圆柱体右端面最低点和油液面高于圆柱体左端面最高点时的情况，在计算总体积时计算。图12图13油面法向量为,且油面过点,由几何关系（如图13）有0,SINCO0Z,故油面方程为1202TALZHR00SICOSTANYZY如图13，分别为左、右端面的Y轴坐标值，又，则可得圆12,TY2215XZR柱体部分储油量为，记为。201TA251TTYZVDD2,TTVY（2）左球冠体、右球冠体体积的计算。为了方便计算左球冠体、右球冠体储油的体积，如图14，我们把这每一侧的球冠体的体积用一平行于平面的截面切割成两部分，且两截面与圆柱体底部的距离分0XY别为H1，H2。这两部分体积我们分别形象地称之为、。阴平首先计算下图14左球冠体阴影部分的体积的公式V阴，左图14图15由几何关系有（图15）0625D11TANTANXXZHRZHRD又斜面过点,斜面的单位法向量为，故斜面方程为10625,HRSIN,0CO1SINCOS0XZ经柱面坐标变换,代入方程得SINXYZ1COT0625SZHR在Z坐标确定时,由图16中几何关系得，为柱面坐标下的2065ARCOS1XZDZZ积分上限；，为的积分下限10COTHR0利用对称性,21016502XZZHRVDD阴，左1220XZ1221COT065652TAN2XHRZZHRDZA由于纵向偏角不会太大,故X很小,为了便于计算，利用微分中值定理做近似计算有12211265TACOT065XXZZHZRVHZHA阴，左上述结果记为1,XVH阴，左阴，左同理可得右球冠体的阴，右2222COT065165TANXHRZZHRZVDZAA阴，右近似为22222TXHZZZR阴，右（）（记为，且等价于，所以程序里用,XH阴，右阴，右V阴，右3,XH阴，左代替，以便于编程。23,X阴，左阴，右然后计算下图17的阴影部分体积,由文献容器内存留液体体积与液位高度函平数关系3有本题中为1M,236VHCDH平CD为3M，C为球冠体的厚度，D为圆柱体的底面直径，H油罐体水平时的油面高度。最后根据上述两半部分公式,对储油罐两端球冠体进行平行罐体分割后再计算储油罐两端球状体的体积左球冠体体积（图18），记为因为也由211,XVHVH阴，左平21H1X1H图16图19唯一确定右球冠体体积（图19），记为322,XVHH平阴，右32VH图17图18从上面可以看出中间段、左球冠状、右球冠状储油体积计算的表达式都可以实现。（3）总体积的计算V经过问题分析，可以明显地得到油面和实际储油罐的对应函数关系分成如下图20所示的三种情况情形I（图21）由图中几何关系有积分下限，；积分上限，,且，1V132TANHY24TY1TTLY即为油液面与罐体顶部交线所L在位置到圆柱体左端面的距离。即2112,TYLR23公式不变化,3V32VH情形II所有,同原公式1带入情形III（图22）由几何关系,积分上限为1V2TANTHY图21图20公式不变化,2V21VH为03综合以上三种情况可得,212132,TYLRMIN,ATTHY2,TTTY1X0TTLY122IN3,AX,HH2体积函数的修正利用以上公式用MATLAB计算体积，由于油罐开始标记时并未发生偏转（见题目的附件2），因此可认为显示高度与显示体积即为两种偏转角皆为0时的对应值。将油位高度代入理论公式，并与实际值进行比较，计算程序见附件十，得出理论储油量与实际储油量的比较图，见图2318192212223242526272842444648552545658662X104图23中红线为计算值，蓝线为实际值。从图中可以看出两者有误差，误差来源可能是在计算体积中应用微分中值定理产生的，设这一误差为。因此对计算函数进0OH行平移修正为（单位M3）212132102,TVYVHLRK其中为待定系数。为了求得该系数，12,K12132,TYVHLR图22图23则原公式变为，则Y为线性函数。在无变位情况下，将代入各部分102YKH0H储油量体积计算公式，从而得出Y值，再利用MATLAB中的POLYFIT函数对Y与拟0合，从而得到，程序见附件十四。再将油位高度代1249,01385K入修正后的理论公式，并与实际值进行比较，得出理论储油量与实际储油量的比较图，见图2418192212223242526272842444648552545658662X104图24从图24可以看出修正后的函数与实际数据吻合得很好，因而修正后的计算公式可以用来确定变位罐容体的变位参数和罐容表值的重新标定。3偏转角的计算简记前述体积计算公式为。利用二分法的原理，对平面0,VH,划一范围不断四等份搜索偏转角。具体算法如下（1）取题目中的附件2中的一组油高，记为。对应累计出油为。又在理IHIV论上，从到的累计出油为。则总误IH1I1,IIVVHA差为；2,IIEVA（2）取定区域为，四等份，取每一小区域的中点11,，计算。,2,34IP,34IEP（3）比较的大小，把总误差最小的点所在的小区域的四个顶点赋给IE，返回（2）继续计算直到总误差达到足够的精度停止。121,在以上算法中，记步骤（3）第J次重复得到的区域为，必1212,JJJJ有，边长以的速度收敛到11121222,JJJJJJJJ120，由闭矩形套定理必定能收敛到一个点。,利用附件十六的搜索数据（由原题提供的附表二计算得到）及附件十一提供的程序，循环10次，计算得，2182504将，代入修正后的体积计算公式与附表二提供的数据比218504较如图25、图26，程序见附件十五。所有具体的数值与相对误差见附件十二（一次性补充进油前的数据）与附件十三（一次性补充进油后的数据），部分数据如表五。从图25、图26可知理论计算值和实际数据几乎重合。而从附件十二中可得相对误差的最大值为405,最小值为0，平均值为0533附件十三可得相对误差的最大值为461,最小值为0，平均值为0591。从图和这些数据可得求得的两偏转,角精度很高，可靠性和准确性很好，可以用来对题中附件2所对应的变位罐容体的罐容表的重新标定。因而，将，这组偏角代入校正后的罐容体储油量计算公2182504式，即可得出高度与体积的关系，程序见附件十七，罐体变位后油位高度间隔为10CM的罐容表标定值见表六。0051152253050100150200250300350图25（一次性补充进油前的数据。红点为计算数据，蓝叉为实际数据）005115225050100150200250300350图26一次性补充进油后的数据,红点为计算数据，蓝叉为实际数据表五显示高度M实际数据L计算数据L相对误差26243149091492500010731772620768456874600043243242610319927198040006172532606670057072300096074232599613636135650005206806258762327423389000494113625821079710919001129943525796492448981000525995125754806581825001456912625695120291190400103915542564110824106860012749446255988346862260033141625254852299322994434915E052539618171805200064942212528238522393800036055682521613179132430004856211251022383323835839173E052508242924330400089468782500117134170960002217813表六（变为后罐容表重新标定值）高度M体积L高度M体积L高度M体积L018033116704244122014501311193042146759021125112219772249313032258613247062351768043730714274772454105055460315302732556300067401816330812658330079521917358842760166081179418386692861767091419519414202963078363976同时，为便于比较计算结果，这里给出无变位时储油罐的罐容表每隔10CM的储油量标定值，见表七，程序只需将附件十七中的都换成0即可。,表七（无变位时罐容表理论计算值）高度M体积L高度M体积L高度M体积L01332811858524617901693151121253214876902177412239832251265033192613267582353647044876714295622455894056780315323822557983068869516352012659887071111717380062761571081349918407812862990091599419435102964071364630综上图形和表格可以看出，在同一油位高度下计算得到储油罐的油量容积和实际储油罐的油量容积吻合的很好，因此在误差允许的范围内，该模型计算得到的结果具有一定的可行性，进而可以说明该模型是正确可靠的。八、模型的进一步讨论和改进对问题一的改进椭圆型储油罐纵向变位后，从附件五的误差数值和油位高度可以看出误差数值随油位高度呈现出先升后降的趋势。虽然计算值和实际值的误差很小，但为了计算更准确，更接近实际，我们对误差随监测油位高度的变化进行了多项式拟合，发现三次拟合效果最好。拟合图形和散点图见图27（散点为附件五中的（理论计算值实际值）/实际值，蓝线为拟合曲线），从图中可以看出拟合效果很好。因而设误差随油位高度变化的函数关系式为。有理论计算函数关系式为32134FXAXA，又，为使矫正后的结果更接近实际值，因而,VH理论计算V理论实际实际相对误差可近似得到（V理论V矫正）/V矫正FH,导出。由以上公式重新标1VFH理论矫正定罐容表的值将更准确。400500600700800900100011000010015002002500300350040045005图27对问题二的改进在前述模型建立的的函数后，对于问题（2）中求偏转角度的解答是通过不断迭V代尝试搜索，只能逐步逼近，方法上不够简洁，且需要计算机帮助。由于偏转角的不同，导致罐体位置的变化，会影响油量进出与高度升降的速度，即关于0,VH有拐点位置。同时又有时的即为拐点。因此，我们可以改进偏转0H20,H0H角度的方法，以简化寻找的难度。首先在偏转角的一定范围内计算不同高度,的油的体积，在进行一次数值差分，再重复进行二次插分IIV11IIVDH，用以代替二阶导数，最小的点即为拐点。以建立偏转角与122IIIDH2I,拐点的近似对应关系表。只要对实