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

椭球形储油罐油量计算菏泽一中高二 36班杨阳指导教师孙国林摘要: 椭球形储油罐的油量计算较为复杂, 本文采用随机模拟的方法, 在试验次数较多的情况下,使得频率可近似代替概率,计算出油量,并给出每隔 油高度的油量表。最后, 本方法具有一般性,对于能方便在三维直角坐标系的不规则图形,都能较快计算出体积。关键字:椭球形油量计算随机模拟一、问题提出及分析在一次参观某小型油田时, 了解到油田开采出来的石油是通过储油罐运送到炼油厂, 储油罐里的油量一般用流量计来估计。而实际上由于石油的特殊性, 油量估计结果往往误差较大, 于是希望能直接通过测量储油罐内油的高度直接估计油量。据调查了解到,运输的储油罐一般多为球形和椭球形。学****了课本知识,教科书中一般给出规则体积的计算公式, 比如球体积计算公式, 圆台计算公式, 但是对于稍微不太“好看”的体积, 比如椭球的一部分( 如上右图) 计算就比较麻烦。因此, 本文以较为复杂的椭球形储油罐为例,对于储油罐的油量进行计算。二、模型假设 1. 忽略储油罐可能存在的变形情况,即储油罐是一个标准的椭球形; 2. 不考虑由于温度、压强变化等原因而引起储油罐的体积变化; 3. 相对于储油罐体积来说,储油罐罐壁可以忽略; 三、模型建立 模型思想我们假想有个正方体能够将椭球形储油罐内的石油都包括进去, 然后向这个假想的正方体随机投点, 那么这些点一部分落在液体区域, 一部分落在液体外的区域, 落在液体区域的点的个数除以总投的点的个数就是正方体的点落在液体区域的概率。用正方体的体积乘以这个概率就是液体的体积。当投射的点足够多时就能保证所求的油量非常接近真实油量。 坐标系建立建立下图所示的坐标系,油罐中心在坐标系原点,其中阴影代表油面; a 表示 x 方向半轴长; b 表示 y 方向半轴长; c 表示 z 方向半轴长; x zy a b c 具体做法 设油量高度为 h ,油液面的方程为 chz??,椭球形油罐的方程表示为 1 2 22 22 2???c zb ya x , 则储油罐液体任意一点可由下面方程表示???????????????????? 2 22 2 2 22 21 1b ya xcz b ya xcz chz 假想存在一个边长为 d2 的立方体,其中d 为a ,b ,c 中的最大值,立方体的中心位于坐标原点,三个坐标轴分别垂直过立方体的面。利用计算机编程可以产生若干组随机数,每组随机数包含三个数,三个数的范围都是) (dd,?, 这三个数分别代表一个点的坐标的 x ,y ,z 值,然后判断这点) (zyx,, 是否在液体所在的空间区域内。若干次试验后,计算出所有组的点落在液体区域内的频率。试验次数足够多,频率代替概率,这个概率乘以正方体的体积 38d 就是所求的液体的体积。 YY i=i+1 i<NUM NUM add dV?? 38 输出 V 开始输入 a ,b ,c ,d 取试验个数和变量初值调用随机函数产生随机点点在液体空间区域 add=add+1 N N 计算机的算法流程(如右图) 四、模型求解对于实际问题,当 m4?a ,m5?b ,m3?c 时,计算储油罐内的油量。首先, 我们简单对模型进行验证。假设m3?h , 即在坐标系中上平面为 0?z , 此