基于模拟退火算法的储油罐变位识别模型.pdf

第 卷第 期计算机应用与软件 年 月 基于模拟退火算法的储油罐变位识别模型张万菊姚凯学何勇席雷鹏（贵州大学计算机科学与技术学院贵州 贵阳 ）收稿日期：。张万菊，硕士生，主研领域：计算机测控，物联网技术。姚凯学，教授。何勇，副教授。席雷鹏，硕士生。摘要储油罐在长期使用的过程中因地基变形会导致发生变位（横向倾斜或纵向倾斜），从而造成出油量与油位高度变化不匹配。针对这样的情况，用定积分的方式建立储油罐油位高度变化、发生横向及纵向倾斜变位角和出油量之间的数学关系；通过非线性最小二乘法建立变位参数辨识模型的目标函数，采用模拟退火算法和 求解变位角的最优解。利用 年全国大学生数学建模竞赛 题的实际测量数据及文献数据对比验证，结果表明了该模型的优越性和实用性。关键词储油罐变位识别模拟退火算法最优解中图分类号 文献标志码 ：（，），引言目前我国储油罐数量众多且存放地理位置较为分散 ，尤其是近几年来，为了满足社会经济发展的需要，我国公路、桥梁、高速路等施工建设行业得到大力发展，这些修建工地地理位置通常都比较分散且偏僻。工地上使用的挖掘机、铲车、大卡车等大型工程车需消耗大量柴油，若将工程车直接开到加油站，不仅要浪费油料和时间，而且也有诸多不便，因此，通常在工地现场安装储油罐，使工程车可以方便加油，提高工作效率和经济效益。偏僻的地理位置使得工地上的储油罐常存在被偷油、漏油无人察觉、人工记录数据混乱、油料补添不及时的情况，因此，需要对储油罐油位高度变化与出油量情况及油罐中剩余油量进行实时无人在线远程监测。影响储油罐监测精确度的一个重要因素是储油罐在使用一段时间后，很多会由于地基变形导致变位书书书第 期张万菊，等：基于模拟退火算法的储油罐变位识别模型 （即油罐横向倾斜 角或纵向倾斜 角），从而出现储油罐罐容表不准确的情况 ，同时导致油位变化与出油量不匹配。针对罐容表不准确的情况，许多学者采用定积分的方式建立倾斜储油罐油位高度和储油量之间的数学模型，并采用误差补偿方法 、非线性最小二乘法 、蒙特卡洛算法 、插值法 、神经网络 、遗传 网络 等方法求解所建立模型中的最优变位角度，修正模型后对罐容表做标定，取得了相应的研究成果。本文采用定积分的方式，建立了发生变位后的储油罐油位高度变化量、变位角和出油量之间的数学模型，采用非线性最小二乘法建立求解该模型最优变位角的目标函数，并利用模拟退火算法和 寻找模型中的最优倾斜角。从 年全国大学生数学建模竞赛 题中随机选择 组数据求解出最优变位角，用剩下的 数据对矫正后的数学模型做检验，证明了该模型参数辨识方法的正确性。通过与实测数据、文献 和文献 中所使用方法的实验结果对比，证明了本文建立的数学模型可以更好地反映油位变化和出油量的关系，对于工地上防止偷油、漏油及监测油量变化是否正常的情况有很好的作用。定积分建立变位数学模型目前工地上使用最为广泛的储油罐是圆柱形和椭圆柱形卧式储油罐。本文通过对使用较为广泛的卧式圆柱形储油罐的分析，根据发生纵向和横向倾斜变位的油罐几何形状及油面，采用定积分的方式建立储油罐油位高度变化、变位角和出油量之间的数学模型。纵向倾斜储油罐柱体半径为 ，油位探测针到左侧球冠体距离为，到右侧球冠体距离为。根据储油罐发生纵向倾斜 角后的几何形状及不同油位高度下油面的形状，对储油罐做分段处理。油罐分段和坐标示意图如图 所示。图 纵向倾斜油罐分段示意图油位高度从（油位高度为从罐底向 轴正方向的高度）下降到 的体积变化即为储油罐的出油量 。左冠 柱体 右冠（）油位高度为 ，油平面 （）满足关系：（，）（）油位高度变化处于区域 和 内时，油位高度超出油位探测针的测量范围，出油量与油位高度变化关系无效，且油面几乎不会处于这两个区域内（正常使用下），因此忽略这两个区域的油位变化和出油量的关系。（）柱体部分出油量计算。储油罐柱体部分、区域内油位高度从 下降到 的出油量分别为：当 时，即油位在 内。（，）（，）槡 （，）槡（）当 时，即油位在 内。（，）（，）槡（）当 时，即油位在 内。（）（）（，）槡（）（）球冠体部分出油量计算。储油罐两端的球冠体高为 ，半径为 ，当油罐倾斜 角时，左侧部分球冠体油位变化对应油面及坐标示意图如图 所示。图 左球冠体油位对应油面由图 可知，左侧球冠体油位高度从 下降到 的理论出油量 左为：左 （）其中：（，）（，）槡 （）同理，当油位高度从 下降到 时，右侧球冠体的理论出油量为：计算机应用与软件 年 右（，）（，）槡 （）式（）、式（）中的 为：槡()横向倾斜根据卧式圆柱形储油罐的几何形状可知，油罐发生横向倾斜时，液面形状及油位高度均无变化 ，倾斜前后的油位测量值代表的出油量不同。横向倾斜后油罐内油位探测针的位置如图 所示。图 储油罐横向倾斜前后油位探测针发生横向偏转后，油位探测针所测得实际油位高度 与未偏转油位 的关系为：（）（）基于模拟退火算法的油罐变位识别模拟退火算法 （）最早是由 等提出来的一种基于蒙特卡罗思想的随机全局寻优算法，算法来源于物理学中固体物质的退火过程与组合优化间的相似性。在某一初温下，伴随着温度参数的不断下降，利用具有概率突跳性的 准则在解空间中随机搜索，在局部最优解时以一定概率接受新解并跳出局部最优解，最终趋于全局最优 ，计算量相比于蒙特卡罗法、最小二乘法显著减少。具体而言：在某一温度 下，在当前状态 （）基础下根据某一指标（梯度下降）产生新状态 （），两个状态对应的能量分别为 和 ，若 ，则将新状态 （）接受为当前状态；否则系统能量变大，此算法并不会直接抛弃新状态 （）（容易陷入局部最优），而是产生一个 ，的随机数 ，若满足 （为接受概率），则接受新状态，否则保持原状态。接受概率 为：（）（）模拟退火算法流程如图 所示。图 模拟退火算法流程储油罐第 次出油前的油位高度为，出油后的油位高度为 ，做横向倾斜矫正后的油位高度为 ，实际出油量为 ，理论出油量为 （，）。采用非线性最小二乘法建立模拟退火算法中的目标函数为：（）（，）（）模拟退火算法中新解的产生流程描述如下：（）随机产生两个 （，）的随机数 和；（）计算温度采用率：槡（）（）根据解 产生新解：（）（）新解有效性检查，若新解无效，则在解空间内随机产生新解。求解储油罐的最优变位角，即是求解使目标函数（）函数值最小的 角和 角。第 期张万菊，等：基于模拟退火算法的储油罐变位识别模型 实验结果分析设置模拟退火算法的初始温度 为 ，退火速率 等于 ，外层循环结束条件为循环迭代 次，内层循环抽样稳定准则为循环迭代 次，解空间为 且 。随机选择文献 附件 中的 组实测出油量及油位高度变化的数据对模型做变位识别。退火适应度进化曲线如图 所示。图 退火进化曲线由图 可知，模拟退火算法输出的模型的最优变位角为 、。将该变位角代入储油罐油位变化、倾斜角与出油量关系模型矫正模型，随机从剩余实测数据中选择 组数据代入矫正后的出油量数学模型，计算变位识别出油量绝对误差。同时，计算出未做变位识别的出油量数学模型的绝对误差，对比结果如图 所示。图 变位识别和无变位识别与实测出油量对照变位识别与未变位识别模型理论出油量与实测出油量对照图如图 所示。可以看出，所建立数学模型做变位识别矫正后，理论出油量更接近实际出油量。将模拟退火算法求解出的最优变位角代入倾斜储图 变位识别与无变位识别对照 图油罐出油量数学模型，计算油位高度从 下降到 的理论出油量，并与实际测量数据和文献 中实验结果比较，计算模型出油量与实测加油机出油量相对误差，可以得到如表 所示的结果。表 实测出油量与模型出油量对比表对比项数值油位高度 出油量 加油机 实测罐容表 遗传算法 文献 网络 本文模型 相对误差实测罐容表 遗传算法 文献 网络 本文模型 油位从 下降到 时，罐容表显示容积变化量即为罐容表显示出油量。根据表 可知，文献 中实测罐容表出油量相对误差约为 ，文献 中采用遗传算法计算出的出油量相对误差约为 ，网络的相对误差约为 ，文献 模型的相对误差约为 ，本文出油量模型相对误差约为 。本文所建油位高度变化与出油量的数学模型经模拟退火算法变位识别后，可以很好地反映加油机出油量与油位高度变化的关系。通过出油量与油位高度关系建立的数学模型在相同油位变化下的油量容积变化的误差较大，其中可能存在累积误差、模型误差及出油机误差等。计算机应用与软件 年 结语本文采用定积分的方式建立储油罐油位高度变化量、倾斜角与出油量的数学模型，通过 、非线性最小二乘法和模拟退火算法求解模型最优变位角。对发生变位后的出油量数学模型做矫正，并利用实测数据检验、文献数据对比分析证明了本文模型和变位识别方式的正确性及实用性和建模方式减少累积误差的优势。储油罐在实际使用过程中，影响出油量与油位变化不匹配的因素还有很多，比如温度、液位计精度等，还需要做更进一步的研究。参考文献赵庆平 储油罐远程监管系统设计与实现 唐山学院学报，（）：，贾跃，刘远征，车秉正，等 储油罐的变位识别与罐容表标定的建模求解 电子世界，（）：邵树琴 储油罐的变位识别与罐容表标定的模型研究 现代盐化工，（）：窦霁虹，梅钰，陈振勋，等 储油罐的变位识别与罐容表标定模型 纯粹数学与应用数学，（）：，：，：李宁靨，邱煜淳，李建蕊 基于 插值法的储油罐变位识别与罐容表标定 南通职业大学学报，（）：，李彬，周斌，李国平 基于神经网络的变位储油罐罐容表标定 中原工学院学报，（）：孙娓娓，盛兴平，黄映雪 基于遗传 网络的变位储油罐识别模型 计算机工程与应用，（）：，廖嘉，谢中华，刘寅立 储油罐的变位识别与罐容表的标定 天津科技大学学报，（）：聂柏松，孙继斌，张丽 年全国大学生数学建模竞赛题的研究 （）：王凌 智能优化算法及其应用 北京：清华大学出版社，：孙丽萍，李佳琪，张希萌，等 基于 算法的优化排料研究 计算机应用与软件，（）：袁汪凰，游晓明，刘升，等 求解 问题的自适应模拟退火蚁群算法 计算机应用与软件，（）：岳秀，张超峰，张伟，等 基于 和改进模拟退火算法的航迹规划 控制工程，（）：唐煜，曹小鹏，张莹 基于模拟退火和历史存档融合的风驱动优化算法 计算机应用与软件，（）：，中国工业与应用数学学会 年全国大学生数学建模竞赛试题 （）：（上接第 页），（）：，：，（）：，：，（）：，（）：，：，（）：，（）：，：，（）：，