没有合适的资源?快使用搜索试试~ 我知道了~
全国数学建模储油罐.pdf
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 137 浏览量
2022-11-15
19:09:22
上传
评论
收藏 1.83MB PDF 举报
温馨提示
试读
26页
。。。
资源推荐
资源详情
资源评论
储油罐的变位识别与罐容表标定
摘要
本文研究的是小椭圆形储油罐与实际储油罐在发生纵向倾斜与横向偏转倾斜等变
化后,储油量与实测油高的关系,从而对变位后的储油罐的罐容表进行重新标定。
本文采用的是微积分知识中分割求和取极限以及等效转化的思想,求得储油量与实
测油高的关系。
问题一,首先将变位后的椭圆储油罐分割成三局部并建立坐标系,分别求得每一局
部水平截面面积与坐标 y 的关系,用 MATLAB 对其进行求积分,得到新的罐容表。运用
给出的倾斜变位储油量和油位高度数据与新罐容表进行比对求误差,得其平均相对误差
为 5%。将题目所给数据与模型得到的数据进行比对,并对误差进行多项式拟合,利用
拟合结果改良罐容表,最终平均误差为 2%。
问题二,分别从数值解与解析解两个角度建立模型。既形象又精确的表现储油量与
实测油高的关系。
首先将储油罐分割为三局部并建立坐标系,参考问题一中微积分的方法得到储油罐
三局部的横截面关于坐标 y 的解析式进行计算。但由于其为超越函数,实际应用中较为
复杂,于是采用微积分中精密分割、求和的思想及坐标旋转变换的关系式,利用 MATLAB
进行数值积分,得到实测油高与实际储油量的关系,即得到标定后的罐容表。运用附件
二中出油量与显示油高的数据进行无限逼近的方法使得实验数据与理论数据的平均误
差和标准差之和最小的方式求解得到了角度α = 3.3750°,β = 4.5000°。
模型二将储油罐中封头局部假设为椭球体,利用其在无变位条件下局部体积随高度
变化的函数较为简单的优点,通过寻找等效液面将实测油位高度转化为无变位条件下的
油位高度,再代入原函数式中得到较为精确的解析解。并最终得到与模型一相似的结果。
对于问题二中的两个模型进行验证,通过题目所给显示油高与显示油量容积的关系
和模型得到的数据进行误差比对;以及通过出油量与显示油高和模型已得到其变位参数
的条件下进行比对,都得到了误差。数据说明模型一较为精确,模型二的误差在允许范
围之内,模型具有较好的正确性与可靠性。
关键词
微积分 多项式拟合 几何学 坐标变换 体积等效
.
一、 问题重述
石油被称为“工业的血液〞,是最重要的战略能源。同时,它也是正在加速枯竭的
非可再生资源。在国际油价居高不下的今天,准确的石油计量既满足了可持续开展的要
求,也保证了相关企业的经济利益。
加油站是常见的供油单位,通常在地下有假设干储存燃油的储油罐,并有与之配套
的“油位计量管理系统〞。油罐在安装后即会进行容积标定,制作储油量与油位高度一
一对应的罐容表。以后就通过测量油位高度,再查表来获取储油量信息。
储油罐在使用一定时间后,由于地基变形等原因,可能会发生纵向与横向的倾斜。
此时,之前的以油罐水平放置为前提的罐容表就不再适用了,需要制定新的罐容表。
此题目分两问,第一问给出了一个平头椭圆罐的尺寸与纵向变位参量,要求结合实
验数据,研究变位对罐容表造成的影响。并给出变位后新的罐容表。第二问给出了某实
际储油罐发生变位后,进行进出油实验所得的数据。要求建立罐内储量与油位高度及变
位参数之间的一般关系。并利用给出的数据,根据模型,确定变位参数,并给出新的罐
容表。最后仍用数据来分析模型的正确性与方法的可靠性。
二、 问题假设
1)
2)
3)
4)
假设不考虑储油罐的形状误差;
假设不考虑温度、气压等因素造成的影响;
假设储油罐的封头可看作球缺状与椭球状两种形式;
假设储油罐倾斜变位参数的范围在0°~10°之间。
三、 符号说明
V
h
S
D
R
L
b
d
a
储油罐的体积〔L〕
储油罐的实测高度〔mm〕
储油罐的水平截面面积〔mm
2
〕
储油罐筒体底面直径〔3m〕
储油罐筒体底面半径〔〕
储油罐筒体的高〔8m〕
储油罐封头的高度〔1m〕
储油罐油位探针到筒体左底面的距离〔2m〕
小椭圆柱形储油罐底面椭圆半长轴〔〕
.
四、 问题分析
根据标准 GB/T19779-2005《石油和液体石油产品油量计算 静态计量》 规定:当
储油罐的纵向倾斜度在 1% ~ 8% 之间时,其罐容表标定值要求采用实际测量所得数据进
行标定;而当其纵向倾斜度在 0 ~ 1%之间时,其罐容表标定值要求套用标准的罐容表,
并注明液高修正值。其中罐容表反映了储油罐中任意高度下的容积,即从油罐底部基准
点起任一垂直高度下该储油罐的有效容积。
由题意知罐体的倾斜角α = 4.1°,倾斜度在 1% ~ 8%之间,因此需要采用实际测量的
数据对其进行罐容表的标定。首先将罐体分为三局部,对于每一局部分别求得高度 h对
应小椭圆型油罐水平截面的面积表达式,再利用积分的方法即可得到罐体变位后油位高
度间隔为 1cm 的罐容表标定值。
对于问题二,针对问题一的微积分方法建立模型,假设较为复杂,那么考虑采用解
析解与数值解的两个角度对模型进行分析改良,从而更为形象准确的得出罐内储油量与
油位高度及变位参数〔纵向倾斜角度α和横向偏转角度β 〕之间的一般关系。
[
1
]
五、 模型建立与求解
问题一
模型的建立
如图 1 所示,按此方法将纵向倾斜变位后的小椭圆型油罐分为三局部:
图 1 纵向倾斜变位后的小椭圆油罐分为三局部示意图
对于第一局部,按图 2 对面积进行积分,
.
图 2 面积积分方法
从而求得高度h对应小椭圆型油罐水平截面的面积S的表达式:
ℎ∙𝑡𝑎𝑛𝛼
𝑆=∫
0
(
𝑏−𝑦𝑡𝑎𝑛𝛼
)
2
2𝑎∙
√
1− 𝑑𝑦
𝑏
2
经计算求得:
对于第二局部,同理按图 3 所示对面积进行积分,
𝑏
ℎ−
𝑡𝑎𝑛𝛼
𝑎𝑡𝑎𝑛𝛼 𝑏 𝑏
2
𝑏
2
𝑏
2
𝑆= ∙[(ℎ− )∙
√
( ) −(ℎ− ) +( ) ∙𝑠𝑖𝑛
−1
]
𝑏
𝑏 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼
𝑡𝑎𝑛𝛼
图 3 面积积分方法
从而求得高度h对应小椭圆型油罐水平截面的面积S的表达式:
𝑆=∫
经计算求得:
8
𝑐𝑜𝑠𝛼
0
{
𝑏−
[
ℎ+
(
𝑙−𝑦
)
𝑡𝑎𝑛𝛼
]}
2
2𝑎∙
√
1− 𝑑𝑦
𝑏
2
.
𝑎𝑡𝑎𝑛𝛼 𝑏−ℎ+2.45𝑡𝑎𝑛𝛼 𝑏
2
𝑏−ℎ+2.45𝑡𝑎𝑛𝛼
2
𝑏
2
S= ∙[𝑥 − ∙
√
−𝑥 − +
𝑏 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼
∙𝑠𝑖𝑛
−1
𝑥−
𝑏−ℎ+2.45𝑡𝑎𝑛𝛼
𝑡𝑎𝑛𝛼
]
𝑏
𝑡𝑎𝑛𝛼
由此亦得第三局部:
𝑏
2 2 2
ℎ
′
−
𝑡𝑎𝑛𝛼
𝑎𝑡𝑎𝑛𝛼 𝑏 𝑏 𝑏 𝑏
𝑆= ∙[ℎ
′
− ∙
√
−ℎ
′
− + ∙𝑠𝑖𝑛
−1
]
𝑏
𝑏 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼 𝑡𝑎𝑛𝛼
𝑡𝑎𝑛𝛼
ℎ
′
=8𝑠𝑖𝑛𝛼+3𝑐𝑜𝑠𝛼−ℎ
在高度h对应小椭圆型油罐水平截面的面积 S的情况下,运用 matlab 求得相隔的
13000 个h所对应的平面面积,并运用积分函数对其积分,从而得到罐体变位后油位高
度间隔为 1cm 的罐容表标定值,如表 1 所示。
表 1 罐体变位后每隔 1cm 高度的罐容表标定值
高度〔cm〕 容积〔L〕 高度〔cm〕 容积〔L〕 高度〔cm〕 容积〔L〕 高度〔cm〕 容积〔L〕
0≤1.664 120≥4012.626
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
.
剩余25页未读,继续阅读
资源评论
春哥111
- 粉丝: 1w+
- 资源: 5万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 堆排序(Heap Sort)是一种基于比较的排序算法
- ebatis 是一个简单方便上手的声明式 Elasticsearch ORM 框架
- 威纶通触摸屏编程软件Easy builder pro V6.09.02安装包(2024.06).txt
- ES查询客户端,elasticsearch可视化工具 elasticsearch查询客户端
- html css js网页制作实例 dldtdd实现列表功能
- 用python制作的tts语音小工具
- 三菱PLC编程参考手册
- 吃豆人代码源码全套.cpp
- 快速了解学习「编译原理」都需要掌握哪些基础知识.pdf
- Verilog示例代码,以SMIC 12nm工艺库为例给出Tessent TCL脚本示例
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功