《多元线性回归模型探析.docx》由会员分享,可在线阅读,更多相关《多元线性回归模型探析.docx(8页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、多元线性回归模型探析摘要:基于液态烃脱硫装置生产操作运行数据,针对实际液态烃脱除H2S情况构建多元线性回归模型,分析相关操作参数与H2S含量的之间影响因素,并通过拟合优度检验、方程显著性检验以及变量显著性检验对模型进行了验证,得到了较为可靠的多元线性回归模型,最后基于回归模型提出降低液态烃中H2S含量解决方案.关键词:脱硫;H2S含量;相关分析;多元线性回归;模型检验脱硫装置原料液态烃来自上游炼厂重油催化装置及焦化装置,原料液态烃进入抽提塔底部,与塔上部注入的溶剂逆流接触、液液抽提,经抽提脱除绝大部分硫化氢.脱硫后液态烃H2S含量常出现异常偏高的现象,影响后续碱洗脱硫醇流程.本文通过分析液态烃
2、进料量、脱硫前液态烃硫化氢含量、脱硫剂循环量、脱硫剂醇胺含量、抽提塔上部温度等因素对脱硫后液态烃H2S含量的影响,并逐一进行验证,建立液态烃H2S含量与各影响因素之间的多元线性回归模型.1数据预处理1.1自变量的选取根据实际操作经历,选取了正常生产中5个可能的影响因素数据,即液态烃进料量1a、脱硫前液态烃硫化氢含量2a、脱硫剂循环量3a、脱硫剂醇胺含量4a、抽提塔上部温度5a.考虑到液态烃脱硫遭到液态烃处理、溶剂循环量、原料硫化氢含量、溶剂醇胺含量的综合影响,在正常操作条件下,对数据进行比值处理,新增液态烃进料量与循环溶剂量比163aaa(料剂比)、液态烃总硫与循环溶剂量比1273aaaa(总
3、硫剂比)、液态烃总硫与溶剂循环总醇胺比值12834aaaaa(总硫醇胺比).将液态烃H2S含量作为因变量b,其它影响因素作为自变量a,H2S含量与挑选后的因素。1.2源数据无量纲化为H2S含量与各因素散点图.H2S含量与其影响因素之间并不是呈简单的线性函数,而是其它类函数,如指数函数.那么对源数据取对数就将其关系转换成线性函数,并进行线性回归.2相关性分析2.1相关概念解释2.1.1Pearson相关系数Pearson相关系数能够相当准确地反响两个变量之间的线性相关程度,判定这些自变量与MSBE含量的线性相关强度.其中相关系数用r表示,描绘的是两个变量间线性相关强弱的程度以及相关性的方向,n为
4、样本量,iX和iY为变量的观测值,X和Y为观测值的均值.r的绝对值越大,表明相关性越强.r为负数,表明因变量与自变量呈负相关关系;r为正数,表明因变量与自变量呈正相关关系.一般而言,r的绝对值大于0.95为显著性相关,大于等于0.8为高度相关,在0.50.8之间为中度相关,在0.30.5之间为低度相关,小于0.3则为不相关.2.1.2P值P值是一个概率,一种在原假设为真的前提下出现观察样本以及更极端情况的概率.在因变量与自变量的相关性分析中,原假设为因变量与自变量无相关性,P值表示对原假设的支持程度2.通过SPSS计算出P值后,将给定的显著性水平(通常取0.05)与P值比拟,就可作出检验的结论
5、:假如p值,则在显著性水平下拒绝原假设,即因变量与自变量相关;假如p值,则在显著性水平下接受原假设,即因变量与自变量不相关.2.2双变量相关性分析由于所选取的自变量较多,为了简化分析,先通过两两变量之间的相关性分析对自变量进行挑选,排除不相关的因素.将取对数后的数据导入SPSS中进行双变量相关性分析,显著性水平取值为0.05。在各自变量(1,2,8)iXi与因变量Y的相关性分析中,排除Pearson相关系数r大于0.3且P值大于0.05的自变量,即排除和溶剂循环量(3X)和抽提塔上部温度(5X)的因素.其余自变量因素P值小于0.05,拒绝原假设,讲明二者相关,其中1X、2X、6X、7X、8X自
6、变量与与脱硫后液态烃H2S含量呈中度正相关,且7X总硫料比Pearson相关系数r到达了0.774.4X自变量与与脱硫后液态烃H2S含量呈中度负相关,即醇胺含量越高液态烃H2S含量越低.从双变量相关性分析中,能够初步得出单因素自变量与H2S含量相关性较低,综合因素自变量十分是液态烃总硫与溶剂循环量比值7X能够显著性地反映出与脱硫后液态烃H2S含量的相关性.3多元线性回归模型的建立3.1模型的建立从相关性分析中,我们确定了脱硫后液态烃H2S含量的影响因素,但却难以区分各影响因素的主次,以及重复影响因子,此时需要运用多元回归模型进行区分和预测.在SPSS软件中对因变量(液态烃H2S含量的含量)和自
7、变量(液态烃进料量、脱硫前液态烃硫化氢含量、脱硫剂醇胺含量、料剂比、总硫剂比、总醇胺比值总硫醇胺比)作多元线性逐步回归.在多元线性逐步回归经过中,其基本思想是:对全部的自变量X按它们对Y奉献的大小进行比拟,并通过F检验法,选择偏回归平方和显著的变量进入回归方程,每一步只引入一个变量,同时建立一个偏回归方程.当一个变量被引入后,对原已引入回归方程的变量,逐个检验他们的偏回归平方和.假如由于引入新的变量而使得已进入方程的变量变为不显著时,则及时从偏回归方程中剔除.在引入了两个自变量以后,便开场考虑能否有需要剔除的变量.只要当回归方程中的所有自变量对Y都有显著影响而不需要剔除时,再考虑从未选入方程的
8、自变量中,挑选对Y有显著影响的新的变量进入方程.不管引入还是剔除一个变量都称为一步.不断重复这一经过,直至无法剔除已引入的变量,也无法再引入新的自变量时,逐步回归经过结束.在多元线性逐步回归经过中,经建模分析确定多元线性回归的自变量X7,排除了其它自变量的组合.并根据回归标准化残差直方图,能够判定出本模型的残差数据服从正态分布.确定的数学模型为3.2模型的检验模型检验是由统计理论决定的,它是运用数理统计的方法,对方程进行检验、对模型参数估计值的可靠性进行检验.主要包括拟合优度检验、方程显著性检验、变量显著性检验3等.3.2.1拟合优度检验(R2检验)为了检验描绘MSBE含量的多元线性回归方程与
9、观测值的拟合优度,我们采取拟合优度检验的方法来讲明.拟合优度检验方法是构造一个能够表征拟合程度的指标R2.其中总变差平方和S总是各个观察值与样本均值之差的平方和,反映了全部数据之间的差异;回归平方和S回是总变差平方和中被回归方程解释的部分,是回归方程中所选取的解释变量Xi对被解释变量Y所造成的影响;S残是总变差平方和中未被回归方程解释的部分.R2的取值范围是,R2的值越接近1,讲明回归直线对观测值的拟合程度越好;反之,R2的值越接近0,讲明回归直线对观测值的拟合程度越差.通过SPSS软件实现R2检验.模型(1)的R2检验结果为0.599,调整R2为0.589拟合优度非常好.3.2.2方程显著性
10、检验(F检验)解释自变量与液态烃H2S含量的含量之间的线性关系,我们计算统计量F的值来讲明总体上被解释变量与所有解释变量之间的线性关系能否显著.通过公式FFINV(,回归自由度,误差项自由度),计算在特定显著度0.05,F临界值=FINV(0.05,1,42)的条件下F检验的临界值.当统计量F的值大于临界值时,即可以为在总体上,自变量与因变量呈显著的线性关系.计算结果见表8,通过比拟,模型(1)F值明显大于F临界值,通过了F检验,且呈显著的线性关系.3.2.3变量显著性检验(t检验)通过t检验我们能够看出每个自变量对MSBE含量的影响,假如某个自变量对H2S含量的影响并不显著,就从回归方程中将
11、它剔除.从以上检验结果中能够看出,模型(1)中的自变量X7和截距T检验的P值都小于0.05,即模型(1)通过T检验.3.2.4模型检验结果通过对模型进行拟合优度检验、方程显著性检验以及变量显著性检验,得出模型一回归效果都特别显著,拟合优度较佳.同时,结合实际生产结果及理论研究,液态烃总硫含量和溶剂循环量在实际操作中对H2S脱除都有重要的影响,即本模型为多元线性逐步回归的最优结果.通过上述分析,我们得出影响脱硫液态烃H2S含量与综合影响因素的多元线性回归方程7Y1.632X15.146.由于作多元线性回归的数据是对原数据取对数后的数据,现将其复原为1231.6321n(/)15.146e.aaa
12、b(2)结果表明液态烃进料量和液态烃原料H2S含量对脱硫液态烃H2S含量的影响较大,在正常温度和醇胺控制范围内,通过提高溶剂循环量,增大脱硫溶剂与液态烃的液液接触面积,可明显降低脱硫液态烃H2S含量.同时发现脱硫溶剂有效组分醇胺含量在较低浓度120g/L期间仍能够知足液态烃脱H2S的生产要求.本文模型建立在正常醇胺浓度120g/L250g/L范围,因而本模型中醇胺对脱硫H2S影响较小,经检验,建模及回归方程符合实际生产运行情况.通过实时的原料液态烃分析,当进料量和原料液态烃H2S发生明显变化时,可根据脱硫建模方程的数据计算,适当提高或降低溶剂循环量;溶剂醇胺含量低于120g/L,建议参加新鲜MDEA脱硫溶剂,保持适当醇胺组分浓度,提高脱硫效果.