基于HHLP理论的直接优化方法在换料方案搜索中应用.pdfVIP

第十三届反应堆数值计算与粒子输运学术会议 基于HHLP 理论的直接优化方法 在换料方案搜索中的应用 王文聪，蒋校丰，赵荣安 上海交通大学核科学与工程学院 上海市东川路800 号,200240 摘要：本文采用HHLP 理论（Hybrid Harmonics and Linear Perturbation ，HHLP ）和混合整数 线性规划（Mixed Integer Linear Programming ，MILP ）方法相结合的直接优化方法来求解换料优 化问题。基于HHLP 理论，可快速、准确地评价换料方案，而将它和混合整数规划（Mixed Integer Programming MIP ） 方法相结合，可将换料方案的搜索过程和评价过程整合到一个数学规划模 型中。由于基于 HHLP 理论建立的仍然是特征值问题，因此直接建立的是一个混合整数非线性 规划（Mixed Integer Nonlinear Programming ，MINLP ）模型。目前商用的数学规划软件不能高效 地求解这样的MINLP 问题，因此本文基于一阶微扰理论建立了MILP 模型。该MILP 模型需要 迭代求解，尽管不能得到全局最优解，但迭代过程是非常高效的。本文应用商用的数学规划软 件IBM CPLEX 来求解该MILP 问题。数值试验表明，该优化方法可获得高质量的优化解。 关键词：换料方案；直接优化；MILP ；HHLP ；CPLEX 中图分类号：TL323 文献标识码：A 由于换料方案的搜索空间非常庞大且 扩散方程是非线性的，这将导致数学规划问 评价过程非常耗时，至今换料方案的优化问 题无法高效求解。为此，需要将中子扩散方 题仍然是极具挑战性的。对于典型压水堆， 程线性化以建立线性规划模型。但从某个参 即使在1/8 堆芯对称装料的情况下，仍然有 考换料方案出发所建立的线性化模型不可 26 能对搜索空间内的所有换料方案都适用，因 10 以上个换料方案。尽管换料优化方法已 经研究多年，但真正工程应用的方法是非常 此必须要将整个搜索空间划分为很多个子 少的。换料优化的方法通常有两类，一类是 空间，而线性规划模型则需要在这些子空间 随机优化算法，如遗传算法[1] （genetic 内独立求解。 algorithm ，GA ）和模拟退火算法 HHLP （Hybrid Harmonics and Linear [2](simulated annealing algorithm ，SAA)，另 Perturbation)理论是有上海交通大学物理组 一类是确定论算法，如基于专家守则的方法 开发的。利用该理论，即使对扰动非常大的 和数学规划方法[3，4] 。对于随机优化算法， 换料方案，仍然可以快速、精确地评价[5，6] 。 换料方案的搜索过程和评价过程是分离的。 将HHLP 理论和数学规划方法相结合，可以 换料方案的搜索过程具有随机性，缺乏有效 将换料方案的搜索过程和评价过程整合到 的物理机制去引导换料方案的搜索，而换料 一个数学规划模型中。由于基于HHLP 理论 方案的评价过程则需要求解中子扩散方程， 建立的仍然是特征值问题，因此以上所建立 非常耗时。而基于专家守则的确定论方法依 的是一个混合整数非线性规划（Mixed 赖于以往的经验。与随机性的优化算法不 Integer Nonlinear Programming ，MINLP ）模 同，在数学规划的方法中，方案的评价过程 型。目前商用的数学规划软件不能高效地求 和搜索过程同时进行，是一种直接优化方 解这样的MINLP