椭球形储油罐油量计算.doc

椭球形储油罐油量计算
菏泽一中高二36班杨阳指导教师孙国林
摘要:椭球形储油罐的油量计算较为复杂,本文采用随机模拟的方法,在试验次数较多的情况下,使得频率可近似代替概率,计算出油量,。最后,本方法具有一般性,对于能方便在三维直角坐标系的不规则图形,都能较快计算出体积。
关键字:椭球形油量计算随机模拟
一、问题提出及分析
在一次参观某小型油田时,了解到油田开采出来的石油是通过储油罐运送到炼油厂,储油罐里的油量一般用流量计来估计。而实际上由于石油的特殊性,油量估计结果往往误差较大,于是希望能直接通过测量储油罐内油的高度直接估计油量。
据调查了解到,运输的储油罐一般多为球形和椭球形。学****了课本知识,教科书中一般给出规则体积的计算公式,比如球体积计算公式,圆台计算公式,但是对于稍微不太“好看”的体积,比如椭球的一部分(如上右图)计算就比较麻烦。因此,本文以较为复杂的椭球形储油罐为例,对于储油罐的油量进行计算。
二、模型假设
,即储油罐是一个标准的椭球形;
、压强变化等原因而引起储油罐的体积变化;
,储油罐罐壁可以忽略;
模型建立

我们假想有个正方体能够将椭球形储油罐内的石油都包括进去,然后向这个假想的正方体
随机投点,那么这些点一部分落在液体区域,一部分落在液体外的区域,落在液体区域的点的个数除以总投的点的个数就是正方体的点落在液体区域的概率。用正方体的体积乘以这个概率就是液体的体积。当投射的点足够多时就能保证所求的油量非常接近真实油量。

建立下图所示的坐标系,油罐中心在坐标系原点,其中
阴影代表油面;
表示方向半轴长;
表示方向半轴长;
表示方向半轴长;


Y
Y
i=i+1
i<NUM

输出V
开始
输入,,,
取试验个数和变量初值
调用随机函数产生随机点
点在液体空间区域
add=add+1
N
N
x
z
y
a
b
c

,油液面的方程为,椭球形油罐的方程表示为,则储油罐液体任意一点可由下面方程表示

,其中为,,中的最大值,立方体的中心位于坐标原点,三个坐标轴分别垂直过立方体的面。利用计算机编程可以产生若干组随机数,每组随机数包含三个数,三个数的范围都是,这三个数分别代表一个点的坐标的,,值,然后判断这点是否在液体所在的空间区域内。若干次试验后,计算出所有组的点落在液体区域内的频率。试验次数足够多,频率代替概率,这个概率乘以正方体的体积就是所求的液体的体积。
(如右图)
模型求解
对于实际问题,当,,时,计算储油罐内的油量。
首先,我们简单对模型进行验证。假设,即在坐标系中上平面为,此时为液体体积刚好是半个椭球。查找相关公式,半椭球的体积为


我们取不同量级的投点个数,即不同的NUM值,用MATLAB编程,结果如下表:
投点个数对结果的影响
投点NUM(10^)
2
3
4
5
6
7
体积V
140
134