储油罐的变位识别与罐容表标定的优化算法(获国家奖).docVIP

储油罐的变位识别与罐容表标定的优化算法 摘要: 本文针对不同卧式储油罐. 关键词：卧式储油罐优化算法 参数估计 （1） 再由柱体体积与面积之间关系，得罐身中油体积计算公式如下： （2） 2.2小椭圆储油罐有纵向倾斜时体积的计算模型 角度后，可将总体体积分成若干个截面椭圆中的面积在求微分和，油面高度分为以下三种情形： （3） 对应三种情形对应的示意图如下所示（图1），其中为取任意位置处垂直于油罐底面的垂直油面高度： 情况1 情况2 情况3 图1 不同油面高度示意图 当时（图1 情况2）， （4） 此椭圆截面上对应面积可似公式（1）得到，进而体积计算公式为： = （5） 令，又由（4）式，可得： （6） 则有： = （7） 其中，； 2）当时，得： （8） 3）当时， 为左边椭圆柱体体积，为右边纵截面为梯形时油的体积，其中： （9） 仍然为（7）式，只是其积分下界值变为1。 2.3 实际储油罐有横、纵向变位时体积的计算模型 角度后，做如下处理：以罐体中截面与油面的交线为旋转轴，将油面逆时针旋转角度，得到切割平面，并如图2所示进行割补，此时球冠体中油的体积即转化为水平位置的情形，罐身仍然采用2.2类似方法进行求解。 首先确定横向变位后实际液面高度与浮标显示数值的关系，则需要求出横向变位前油面高， 图2 实际油罐纵向变位示意图 图3横向变位示意图 由如图所示，可得： （10） 2.3.1罐身体积的推导 ，对于罐身体积的计算可类似于2.2中的推导，其三种情形分别如下： 1） 当时， 令：，得： = (11) 其中，； 2) 当时，体积计算式仍为（11）式，只是. 3）当时，，其中切割部分可用（11）式计算，只是；且左边圆柱体体积：. 2.3.2实际油罐罐头油量的计算 图4 水平位置罐头油品容积示意图 因球冠体为球体的一部分，故其体积计算公式推导如下: 如图4所示建立空间直角坐标系，设球心坐标为,则其球面方程为： 由几何关系得球面半径为，，则 再根据体积计算积分公式： 其中积分区域D为： 由近似处理的特点，此时.经计算得： （12） 当出现类似于2.2中情况1和3时，仍然采用割补的思想，公式（12）仍然适用，就不再进行推导. 从而，由（11）、（12）可得到实际储油罐体积计算. 但因(12)式中含有定积分，不易进行一般性的推广，因此采用泰勒展式和曲面面积近似代替的方法并利用MATLAB6.5软件做近似求解. 3.模型的求解 3.1无变位情形 .2纵向变位情形 时，得到油面高度区间分为[0,430.2]，(430.2,2856.6)，[2856.6,3000]三段，将变位后进油数据代入相应体积公式，模拟累加进油量，给出实际与模拟出的累加进油量每个误差值在[1.59%,5.47%]范围内波动，平均误差为3.782%；出油时每个误差值在[1.59%,5.47%]范围内波动，平均误差为2.155%. 同时运用MATLAB6.5软件，给出油位高度对应的罐容表标定值（但由于数据量较大，因此仅给出间隔为0.1m的12个标定值），结果见下表. 表1：油位高度间隔为0.1m的罐容表标定值 按0.1m渐变油位高度 体积（应变量）/L 按0. 1m渐变油位高度 体积（应变量）/L 0.1 619.48 0.7 2232.5 0.2 281.86 0.8 2661.4 0.3 595.25 0.9 3072.4 0.4 965.66 1.0 3450.7 0.5 1371.9 1.1 3776.6 0.6 1798.5 1.2 4353.4 3.3变位参数估计 ，然后以理论出油量与实际出油量偏差绝对值的平均值最小，即： 为标准，得到参数的估计值。 分别采用梯形面积分割法和泰勒公式法对公式（12）中定积分部分做近似计算，并运用MATLAB6.5软