一种基于蒙特卡洛法的化工过程参数敏感性确定方法与流程

本发明属于化工过程参数敏感性确定领域，具体涉及一种基于蒙特卡洛法的化工过程参数敏感性确定方法。

背景技术：

化工生产过程绝大多数为多参数非线性系统，具有两个显著的特征：一是参数众多，二是参数之间存在着强烈的交互效应，由于化学反应系统的复杂性，尤其是系统状态对操作条件的变化相当敏感，影响化学反应系统过程的设计参数和操作参数是多方面的，一个参数微小变动便可能使系统的输出结果产生很大的变化。对化工过程参数敏感性的研究有助于理论上揭示出反应系统安全操作范围，使可能引起的不良后果能在设计和操作中得以避免。

参数敏感性分析有助于确定模型的关键输入值(参数和初始条件)，同时定量得到输入值的不确定性如何影响模型结果，有助于定量得到一个模型的可靠程度，从而能够确定合适的模型选择，同时能够帮助我们识别显著影响结果的主要参数及其影响大小。前人关于安全临界判据的预测大多采用局部敏感度分析方法甚至是简单的数据直接拟合，往往脱离了反应的本征动力学，在实际复杂工况下难以精确预测反应结果。单因素敏感度分析的方法较为简单，但其不足在于忽略了诸多因素之间的相关性。实际上，一个因素的变动往往也伴随着其他因素的变动，多因素敏感度分析，即全局参数敏感性分析，考虑了这种相关性，能反映几个因素同时变动产生的综合影响，弥补了单因素分析的局限性，全面地解释了事物的本质。全局参数敏感性分析在辨识复杂模型的不确定性方面有明显的优势，其重要目的之一是确定模型的关键输入值(参数和初始条件)，同时定量得到输入值的不确定性如何影响模型结果。

蒙特卡洛(montecarlo)方法作为一种有代表性的数值求解方法，在重要性度量、模型正交分解等方面具有独特的优势，在计算物理学、化学反应工程、生物医学、宏观经济学等领域应用广泛。特别是近年来随着非线性科学的发展，基于蒙特卡洛法的全局参数敏感性的思路被认为能够有效地解决高维复杂模型建立的难题。

技术实现要素：

针对现有技术中存在的上述技术问题，本发明提出了一种基于蒙特卡洛法的化工过程参数敏感性确定方法，设计合理，克服了现有技术的不足，具有良好的效果。

为了实现上述目的，本发明采用如下技术方案：

一种基于蒙特卡洛法的化工过程参数敏感性确定方法，按照如下步骤依次进行：

步骤1：确定输入变量与分布

全面分析待考察的反应系统，识别相关变量，并确定变量波动范围；

步骤2：建立反应系统模型

根据反应系统特点，建立一维/二维和拟均相/非均相模型，建立物料平衡与能量平衡方程，结合反应动力学与热力学参数，建立反应器温度、冷却介质温度、转化率之间的代数关系：y＝f(x)为函数表达式，x＝(x1,x2,x3…,xn)为n维输入变量，每一个变量均有一个概率密度函数；

步骤3：数据分解

将待测值分解成信号值(s)与误差(e)，以输出参数为纵坐标，输入参数为横坐标，将输入参数值和对应的输出参数结果表示在散点图上，定性分析输入参数的不确定性对输出参数的影响，观察输入参数与输出参数之间的关系，输入参数与输出参数之间是否存在非线性关系等，然后利用全局敏感性分析方法做进一步的深入分析；

步骤4：随机抽样模拟

利用蒙特卡洛方法模拟生成所需的抽样样本，加入输入参数的不确定性表征，代入反应系统模型中计算相应的输出与误差偏差，定量分析输入参数对输出结果的影响，考虑参数间的耦合作用，分别计算每个输入变量对输出结果的总方差的贡献；

步骤5：偏差计算与计算结果分析

由各输入参数及各参数间耦合作用共同得到模型的总方差，进而将模型输出方差分解，计算每个参数直接贡献比重和通过参数之间耦合作用间接对模型输出总方差的贡献比重之和，实现误差溯源，得到各个参数的总敏感度指数。

本发明所带来的有益技术效果：

本发明通过蒙特卡洛法计算得到各操作参数之间的内在联系及其对系统安全关键参数的影响规律，以及不同参数变化对系统安全性的高阶影响规律，避免了热失控等危险结果的发生，进而有助于确定多参数共同影响下反应系统的安全边界条件，克服了传统的局部敏感性分析方法脱离反应本征动力学、忽略变量间的相互作用等缺陷，能够为化工过程安全稳定运行提供科学指导与保障。

本发明相比于传统的参数敏感性分析方法，考虑了各个变量之间的相互影响，大幅提高了模型的可靠性，研究人员与工程技术人员可利用各参数敏感性系数的结果和对系统主要安全指标的影响规律解决相应的问题，对于确定化工过程安全操作范围、保障化工设备安全运行具有重要的应用推广价值。

本发明方法可得到各种反应模式下的热失控的临界判据，可简便地推广应用于多种危险工艺过程的安全操作域确定，具有运算快捷、结果准确、重复性好的优点。

本发明主要通过蒙特卡洛法进行化工过程全局参数敏感性计算，进而可以确定反应系统安全操作区域，可以运用到间歇式/连续式、平行及串联反应等多种反应系统。

附图说明

图1为一种基于蒙特卡洛法的化工过程参数敏感性确定方法流程图。

图2为某环氧化工艺全局参数敏感性指数示意图。

具体实施方式

下面结合附图以及具体实施方式对本发明作进一步详细说明：

本发明方法主要包括反应模型建立、参数敏感性分析计算两个部分。详细说明如下：

1、反应模型建立：根据不同反应动力学特性建立能量和物料平衡方程，结合设备性质与冷却介质的传热建立反应器模型，基本操作参数包括初始浓度c0、初始进料温度t0、冷却介质温度tc、反应器参数(长度、直径及其他可能存在的参数如搅拌、内构件等)、催化剂参数(直径、密度、孔隙率等)、操作压力p以及其他传递系数等，不同工艺需根据主要的影响变量分别建立模型。

2、全局参数敏感性分析计算：根据反应模型y＝f(x1,x2,x3…,xn)(xi表示模型的第i个参数)，首先构建模拟总体，进而令每个参数在可能的取值范围内进行蒙特卡洛抽样模拟，研究和预测这些参数的变动对模型输出值的影响程度，考察变量在全部输入空间变化的综合影响。针对不同过程非线性的复杂变化规律分析，可以采用基于模型输出结果方差分解的方法或者基于取样的计算方法。

本发明基于蒙特卡洛模拟的化工过程全局参数敏感性确定方法由两大部分五个步骤组成，如图1所示，具体按照如下步骤依次进行：

步骤1：确定输入变量与分布

全面分析待考察的反应系统，识别相关变量，并确定变量波动范围；

步骤2：建立反应系统模型

根据反应系统特点，建立一维/二维和拟均相/非均相模型，建立物料平衡与能量平衡方程，结合反应动力学与热力学参数，建立反应器温度、冷却介质温度、转化率之间的代数关系：y＝f(x)为函数表达式，x＝(x1,x2,x3…,xn)为n维输入变量，每一个变量均有一个概率密度函数；

步骤3：数据分解

将待测值分解成信号值(s)与误差(e)，以输出参数为纵坐标，输入参数为横坐标，将输入参数值和对应的输出参数结果表示在散点图上，定性分析输入参数的不确定性对输出参数的影响，观察输入参数与输出参数之间的关系，输入参数与输出参数之间是否存在非线性关系等，然后利用全局敏感性分析方法做进一步的深入分析；

步骤4：随机抽样模拟

利用蒙特卡洛方法模拟生成所需的抽样样本，加入输入参数的不确定性表征，代入反应系统模型中计算相应的输出与误差偏差，定量分析输入参数对输出结果的影响，考虑参数间的耦合作用，分别计算每个输入变量对输出结果的总方差的贡献；

步骤5：偏差计算与计算结果分析

由各输入参数及各参数间耦合作用共同得到模型的总方差，进而将模型输出方差分解，计算每个参数直接贡献比重和通过参数之间耦合作用间接对模型输出总方差的贡献比重之和，实现误差溯源，得到各个参数的总敏感度指数。

以某固定床环氧化反应工艺为例，按照本发明方法考察了不同参数(包括温度、浓度、压力、空速、冷却介质等)对系统安全参数(包括反应热、尾气氧气含量等)的影响规律，使用拟均相二维模型结合冷却介质的传热建立反应器模型，采用全局敏感性分析方法，结果如图2所示。

以对体系安全风险影响较大的尾气氧气含量为例，由图2可知，变量一阶敏感度排序由大到小依次为：空速，温度，原料ph，醇水比，双氧水浓度，压力，烯水比；而考虑高阶影响之和的大小次序为：温度，空速，原料ph，醇水比，双氧水浓度，烯水比，压力。

当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。