数学建模A题参考答案.docx

数学建模A题参考答案.docx

《数学建模A题参考答案.docx》由会员分享，可在线阅读，更多相关《数学建模A题参考答案.docx（26页珍藏版）》请在冰豆网上搜索。

数学建模A题参考答案

题目储油罐的变位识别与罐容表标定

摘要

本文分别建立了小椭圆型储油罐及实际储油罐的变位识别模型。

针对小椭圆型储油罐的变位识别问题，采用积分方法，给出无变位时储油量与油位高度的计算公式并得到正常的罐容表标定。

对于小椭圆型储油罐纵向倾斜变位问题，讨论了其截面是三角形和梯形两种情况，利用积分法给出了纵向倾斜变位问题的计算公式，给出了修正后的罐容表标定值，并与正常标定值进行比较。

针对实际大储油罐的变位识别问题，给出无变位时储油量与油位高度的计算公式，根据计算公式得到正常罐容表标定值。

对于倾斜变位问题，用积分方法在不同油高下分别计算出球冠部分和中间圆柱体部分的油量，并求和给出大储油罐纵向倾斜变位后的修正公式。

然后对储油罐横向偏转角度进行分析，给出横向偏转后实际油面高度与正常时油面高度的关系式。

最后结合纵向倾斜角度及横向偏转角度参数公式推导得到罐内储油量与油位高度及两个变位参数间的函数式。

结合附件二中所给数据，利用非线性最小二乘法通过遍历搜索算法求出纵向倾斜角度及横向偏转角度值，最后利用附件二中的数据对模型的可靠性进行了检验，检验结果表明模型较为合理。

关键词：

积分，数值积分，复化梯度法，非线性最小二乘法，罐容表，标定

一、问题的重述

通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，我们可以采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。

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

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

我们采用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题，并解决以下两个问题。

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

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

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

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

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

二、问题的假设

（1）向罐内注入的油量数都是通过流量计来完成，是准确的；

（2）罐内的储油量只有通过加油机加油流出，并且加油机的计量误差在允许的范围内；

（3）不计外部环境的变化对内部油量所产生的影响。

（4）浮标的大小相对于溶剂来说可以忽略。

（5）储油罐中油的密度是均匀的，不考虑水汽、重油等因素的影响。

（6）储油罐的厚度可以忽略。

（7）设注油期间油量无损耗。

（8忽略出油管、检查口、注油口以及油位探针对油位高度的影响。

部分假设在题中给予说明

三、问题的分析

问题一是利用小椭圆型储油罐模型研究变位对罐容表的影响。

在无变位的情况下，储油罐的储油量就是对小椭圆型储油罐进行积分；在变为后，要分三步计算油的体积，第一步，在油平面未到达右端底部时，可以沿垂直于油面和地面的方向截得三角形切面，以油面到椭圆原点的距离为L，求出三角形面积，然后再积分；第二步，当油平面到达右端底部之后，可以用先前的结果减去虚拟部分的结果；第三步，当油平面上升到左端的上沿后，储油量为总体积减去上部空余部分的体积，空余部分的体积和第一步的算法相同。

问题二是一个求实际储油罐变位参数的问题，由于平位时储油罐内液体的体积是一个比较规则的立体图形，因此可以用三重积分的方法求出平位时不同高度时液体体积的理论值，即罐容表的理论值，然后再利用积分的方法求出罐内液体体积与纵向偏移角度、横向偏移角度的关系,建立一个体积与变位参数的关系模型，用这个关系模型求出的相关数据和题中给出的数据进行对比，利用最小二乘法实际的变位参数。

四、符号说明

符号表示

含义

单位

油位高度

油位高度为的两端冠球体储油量总和

油位高度为时圆柱体的储油量

油位高度为时的储油总量

L

贮油罐中间圆柱形的长度

剩余的符号在解题的过程中说明

五、模型的建立与求解

第

（1）题

（一）首先建立罐体未变位时罐容表标定值模型

设油位高度为h，截面作对应的面积为S，对应的罐容表的标定值为

图1-1储油罐横截面坐标系图1-2整个储存罐的坐标表示

正常时高度是已知的，只需求出截面的储油面积：

带入得到体积的公式：

根据此函数可以得到理论值，与数据中的值在同一图中用MATLAB进行拟合，可以得到图1-3所示图形：

图1-3无变位是理论值与数据值的对比

（二）建立罐体变位时罐容表标定值模型

见如图1-4

图1-4变位时的储存罐的坐标表示

L表示油平面到椭圆的中心o的距离，L可以为负数；

a为长半轴，b为短半轴；

为倾斜角（）；

先表示出投影三角形的面积，然后再对Y轴积分；

图1-5V与外部虚拟部分的关系

图1-6标高h与L的关系

下面分别计算三种情况下的V函数，即V关于h的数学模型模型：

为方便输入，一些复杂的表达式由字母代替

①0

②0.1469

③1.17

其中

把求得的102个值与进油时的数据在同一图中比较，发现吻合的非常好，说明我们的模型是可行的。

如图1-7

图1-7变位后理论值与数据的对比

罐体变位后油位高度间隔为1cm的罐容表标定值如表1。

表1：

变位后油位高度间隔为1cm的罐容表标定值

高度/m

容积/L

高度/m

容积/L

高度/m

容积/L

高度/m

容积/L

0.01

3.5380

0.31

630.1550

0.61

1841.8820

0.91

3112.4600

0.02

6.2710

0.32

665.5850

0.62

1885.1330

0.92

3151.4020

0.03

9.9820

0.33

701.5320

0.63

1928.7240

0.93

3190.3330

0.04

14.7640

0.34

737.9800

0.64

1972.5660

0.94

3229.1200

0.05

20.7080

0.35

774.9620

0.65

2015.4560

0.95

3266.9230

0.06

27.8610

0.36

812.3540

0.66

2058.7600

0.96

3304.6320

0.07

36.3250

0.37

849.9700

0.67

2102.1940

0.97

3342.0310

0.08

46.1540

0.38

888.2460

0.68

2145.7160

0.98

3378.7240

0.09

57.4070

0.39

926.8420

0.69

2190.0650

0.99

3415.0080

0.10

70.1430

0.40

965.6660

0.70

2232.5120

1.00

3451.2140

0.11

84.4350

0.41

1004.9580

0.71

2275.6870

1.01

3486.4320

0.12

100.2660

0.42

1044.7620

0.72

2319.1350

1.02

3521.2250

0.13

117.7600

0.43

1084.8500

0.73

2362.3410

1.03

3555.4440

0.14

136.9320

0.44

1124.9720

0.74

2405.4020

1.04

3589.1290

0.15

157.8320

0.45

1165.4910

0.75

2448.4360

1.05

3622.2550

0.16

180.2680

0.46

1206.3620

0.76

2491.3850

1.06

3654.5060

0.17

204.3410

0.47

1247.2300

0.77

2534.1380

1.07

3686.2140

0.18

228.9180

0.48

1288.5640

0.78

2576.7220

1.08

3717.0120

0.19

254.8880

0.49

1330.1220

0.79

2619.2560

1.09

3747.4820

0.20

281.8650

0.50

1371.8880

0.80

2661.5380

1.10

3776.8250

0.21

309.7630

0.51

1413.6820

0.81

2704.1520

1.11

3805.5080

0.22

338.5420

0.52

1456.2130

0.82

2745.6210

1.12

3833.3290

0.23

368.1490

0.53

1498.1450

0.83

2788.0020

1.13

3860.1200

0.24

398.5630

0.54

1540.6520

0.84

2828.8520

1.14

3885.9250

0.25

429.8720

0.55

1583.6720

0.85

2870.1360

1.15

3910.6660

0.26

461.5910

0.56

1626.3860

0.86

2911.2650

1.16

3934.2850

0.27

494.2300

0.57

1669.2700

0.87

2952.0810

1.17

3956.4690

0.28

527.2400

0.58

1712.3230

0.88

2992.5580

1.18

3977.0180

0.29

560.9870

0.59

1755.4250

0.89

3032.6030

1.19

3996.9020

0.30

595.2520

0.60

1798.6320

0.90

3072.5700

1.20

4013.3450

正常情况下的容积：

高度/m

容积/L

高度/m

容积/L

高度/m

容积/L

高度/m

容积/L

高度/m

容积/L

0.01

5.3

0.25

620.4

0.49

1578.1

0.73

2617.5

0.97

3559.6