CRU稳定汽油初馏点软测量建模.docx

1、题：CRU稳定汽油初馅点软测量建模摘要初馅点是稳定汽油生产过程中的一个重要指标，它表示轻组分含量的多少。在现在的生产中，初俺点基本是在化验室由人工化验得到。由于化验的时间比较长，不能用于实时控制。随着计算机技术的发展和控制媒介的提高，软测量技术应运而生。软测量技术是利用统计、辨识、人工智能等方法建立数学模型，将不可测的量作为主导变量，通过测量与主导变量相关的量（辅助变量）来推测主导变量。软测量技术可以连续不断的预测不可测量，解决了过程生产中重要指标靠化验时间滞后大的问题，达到了以软件代替硬件的目的。关键词：初僧点;软测量;主导变量;辅助变量生产过程中，试图确定各种辅助变量的取值范围内，选择记录

2、点。但同时应注重价值的数据点应合理分布范围，即这些数据应该代表:样品应均匀分布特征附近。每个特征点的数据应合理分配近，甚至，应该不会耗费太多的样品和特征点的边缘点，少索取样品，所以极易造成模型泛化能力削弱甚至建模失败。所以点的选择记录应注意合理分布，而且对样本点“缩小”，以确保数据的代表性。1.4.2数据选择与处理软测量技术是通过计算过程数据血实现的，为了保证软测量的额准确性和有效性，对采集来的数据进行筛选时，要注意数据的有效性和精简性，均匀合理的分配采样点，最大限度的拓宽数据的范围，减少信息的重复，避免信息的冗余。除此之外，由于仪表测量误差和环境噪声的影响，对输入数据进行预处理是软测量技术必

3、不可少的一步。输入数据的准确性直接影响软仪表的精度。对输入数据预处理主要在以下两方面进行：数据变换和误差处理。数据的变换对模型的精度，非线性的映射能力有很大的影响。数值变换包括标度、转换和全函数三部分。标度用于解决测量量数据太大的问题，从而改变算法的精度和稳定性。转换用于解决非线性的问题，通过直接转换或者寻找新变量代替原变量的方法来降低非线性。权函数主要用于对变量动态特性的补偿。对数据进行误差处理是十分重要的步骤，误差处理保证了数据的准确性，误差的存在会导致数据的失效，任何数据的失效都会降低模型的精度，甚至是整个软测量失败。误差可以分为随机误差和粗大误差两种。随机误差是由随机因素引起的，可能是

4、操作过程中的微小扰动环境噪声等等，随机误差是不可避免的，但是随机误差符合统计规律，可以用滤波的方法减小，如算数平均滤波，阻尼滤波和小波滤波等。粗大误差的处理包括利用3。法、随机搜索法、神经网络等。1.4.3软测量模型的辨识和在线校正软测量技术的关键是模型辨识，辨识得到的模型的精度直接影响软测量的结果。随着软测量技术的发展，因为软测量采用的理论工具和所针对的实际对象的不同，现在有多种软件测量可以建模。这其中应用最多的是统计回归方法，如模式识别建模，机理模型建模，神经网络技术等等。由于软测量对象随时间变化，非线性模型本身不完整性等特点，该模型在生产过程的软测量建模在线校正是不可缺少的一部分。在线修

5、正的模型结构，包括校准和模型参数标定。因为调整模型结构需要大量的生产数据和长期的训练时间，所以通常意义上在线校正主要是指对参数修正。目前，软测量中使用最广泛的是与主导变量的动态特性相近，关系密切的可测量参数，例如精馄塔的温度。软测量技术在控制领域的研究和实践中取得了很大的成功，但是目前没有形成系统的理论。第1.5节文章的研究内容和结构安排1.5.1文章研究的主要内容本文主要研究了工业过程软测量建模方法，非线性，根据工业过程参数的多重相关性的特点，基于多元逐步回归和偏最小二乘法的方法来解决此类问题的软测量模型。同时熟悉掌握了催化重整炼油工艺，并将此方法应用在催化重整的炼油生产过程中稳定汽油初馅点

6、软测量模型，验证了该方法的有效性。本论文的主要完成的工作如下：（1）掌握重整装置的生产原理和工艺流程，从机理上认识个变量对初儒点的影响。（2）应用多元逐步回归方法和机理分析选择的辅助变量建立软测量模型。（3）采用偏最小最小二乘法选取主元，建立软测量模型。（4）利用主元分析的方法选取主元，再利用偏最小二乘法建立软测量模型。（5）把建立的数学模型应用到催化重整装置中稳定汽油初馅点软测量建模，并检验模型的有效性。1.5.2文章的结构安排本文的结构安排如下：一. 绪论：阐述了论文研究的理论依据和行业背景，明确的课题的研究意义。总结了软测量建模研究的现状并有一个综合评价方法。最后，总结论文的结构安排。二

7、. 基于统计回归技术的软测量方法:着重研究了统计分析的软测量建模方法，对多元逐步回归及偏最小二乘法做了详细的阐述和研究，通过数值试验分析验证了其局限性。三. 掌握重整装置的生产原理和工艺流程，从机理上认识个变量对初饲点的影响。四. 应用多元逐步同归方法和机理分析选择的辅助变量建立软测量模型。基于主元分析选择主元，利用偏最小二乘法建立数学模型，将模型应用到软测量建模中，并检验其有效性。五. 把建立的数学模型应用到催化重整装置中稳定汽油初僧点软测量建模，并检验模型的有效性。六. 结论与展望:对全文内容进行归纳和总结，并对今后的研究工作进行展望和规划。第2章基于统计学的回归建模方法统计回归技术主要是

8、通过统计理论建立输入和输出数据的数学模型。因为它不是完全依赖于生产过程的机理和工艺分析，为复杂工业过程的统计回归技术可以绕过复杂的过程和机理分析，通过统计学习理论，对生产过程的数学模型。统计回归主要可分为线性回归和非线性回归:线性回归分析要求自变量和因变量之间的近似线性关系见面，所以只适用于简单的线性拟合关系;非线性回归是指自变量和因变量是之间的非线性关系。线性回归和非线性回归方程形式是已知的，对模型参数进行了估计，当然，常用的方法:最小二乘法拟合方法，多项式回归方法，回归方法，利用偏最小二乘法等。但是在复杂工业过程中，多种因素的影响主要变量，它很难给出一个方程和机制，因此如何之机制，建立一个

9、完全基于样本模型成为众多学者的研究。目前比较广泛应用的“黑箱”的建模方法主要有神经网络方法等。本章讨论统计回归技术中的多元逐步回归方法和偏最小二乘方法。第2.1节多元逐步回归方法A. 思想多元回归分析方法(基础最小二乘法)；特点是双向筛选；B. 计算步骤：1) 确定F检验值在进行逐步回归计算前要确定检验每个变量是否显若的检验水平，以作为引人或剔除变量的标准。检验水平要根据具体问题的实际情况来定。一般地，为使最终的同归方程中包含较多的变量，水平不宜取得过高，即显著水平a不宜太小，水平还与自由度有关2) 逐步计算(a) 计算全部自变量的贡献(偏回归平方和)Vi。(b) 在已引入的自变量中，检查是否

10、有需要剔除的不显著变量。这就要在已引入的变量中选取具有最小值的一个并计算其F值，如果FVFi,表示该变量不显著，应将其从回归方程中剔除，计算转至(c)。如F>Fi则不需要剔除变量，这时则考虑从未引入的变量中选出具有最大值的一个并计算值，如果F>F2，则表示该变量显著，应将其引入回归方程，计算转至(c)。如果FVF2,表示已无变量可选入方程，则逐步计算阶段结束，计算转入3)。（c）剔除或引人一个变量后，相关系数矩阵进行消去变换，计算结束。其后重复（a）（c）再进行下步计算。3）其他计算主要是计算回归方程入选变量的系数、复相关系数及残差等统计量。逐步回归的具体计算步骤：（1）假设有n个

11、自变量Xio,Xii,Xii,因变量为y,有k个观测点。根据最小二乘原理，y的预测值为：/=bi（）Xio+.+biiXii+bn（2-1）其中，0<i0<il.<il<n-l且各xit（t=l,2,3.1）是从n个自变量中按一定显著性水平筛选出来的检验为显著的因子。筛选过程如下：首先做出（n+1）*（n+l）格式的系统初始相关矩阵rrOOrrOOR=rOy(2-2)ryOryOry,y矩阵中的元素是:(2-3)Z(呵-印(号-寻)*=1£（心-吊）£（勺-弓）2jt=iV&=1其中第n个变量为因变量y(2-4)(2-4)习=一£易

12、(，=1,2,?+1)i=（2）计算偏回归平方和Vi=誓(i=o,l,2,Vi=誓(i=o,l,2,n-1)(2-5)（3）若Vv（）,则对应的冶是已被选入的因子。若所有Vi<0的Vi中选出的Vmin=minVI,其对应的因子为xmin，然后检验因了Xmin的显著性。若按一定显著性水平选取F出:(2-6)（5二噤）（-，-1）<Flliyy则剔除因子Xmin,并对系统相关矩阵R进行该因子的消元。转（2）（4）若Vi>0,则对应的Xi是未被入的因子若所有Vi>0的Vi中选出的Vmax=maxVi，其中对应的因子为Xmax，然后检验因子Xmax的显著性，若按一定显著性水平选

13、取F讥：(2-7)任+|）.（一/一2）尸脾）_'句>f_一成gyyFi.niax则因子Xmax被选入，并对系统相关矩阵R进行该因子的增元改变，转（2）直到上述过程没有因子可以剔除为止。进入回归方程各因子的系数:(2-8)常数项：如=亍-2岫(2-9)iel回归值和偏差值："二如+人土，偏差值£=叫一丸（2-10）is/第2.2节偏最小二乘方法2.2.1偏最小二乘法的特点1. PLS是一种能处理多个因变量更多论述建模方法。尤其是当变量在设定更高度的相关性，使用PLS回归建模分析和比较了因变量，是一个多元回归更有效、更可靠整体性强的建模方法。2. PLS能解决许

14、多以前使用常见的多元回归分析方法不能解决这个问题。诸如独立变量之间的多重相关性问题和样本点不宜太小等问题。3. PLS能实现多种数据分析方法的综合应用。它可以设置多元线性回归法、主成分分析和典型相关分析为一体的基本技能。在一个PLS计算中，不仅可以得到更多的多变量的因变量回归模型，可以分析两组变量之间的关系,同时观察样本点相似结构。这使得数据系统分析的内容更加丰富，同时回归模型为更多的细节的实际的解释。4. PLS允许在最终模型包含所有原来的自变量，最大限度地利用数据信息，在相同的数据信息环境下PLS比普通的多元回归模型具有较高的效率。5. 在建模方案的数据结构简化，可以在二维平而的多维数据特

15、征的观察，强大的图形功能。2.2.2偏最小二乘法的建模步骤1. 将X与Y进行标准化处理后，得到自变量矩阵E<)=(Eoi,Eo2,.,EoP)nxp和因变量矩阵Fo=(FOi,Fo2,.,Foq)nxq。标准化处理是为了公式表达上的方便和减少运算误差。&=(气)"(2-11)其=或=4：s，Sy(2-12)/=j=1,2,.,pk=l,2,.,q2. 记ti是E()的第1个成,ti=Eoa)i,wi是E()的第1个轴，为一个单位向量，既|wi|=lo记ui是Fo的第1个成分，ui=Focioci是Fo的第1个轴，并且|ci|=lo如果要ti和山能分别很好地代表X与Y中的

16、数据变异信息，根据主成分分析原理，应该有：Var(ti)-maxVar(ui)今max(2-13)另外,因为回归建模的需要，要求h对山有很大的解释能力，根据典型相关分析的思想山与ui的相关度应达到最大值，即：r(ti,ui)3max(2-14)综上所述，在偏最小二乘法回归中，我们要求ti与ui的协方差达到最大，即：Cov(ti,ui)=r(ti,ui)>max(2-15)正规的数学表述如F：colwl=1clcl=1因此,将在|CO1|2=1和|C1=1的约束条件下，去求coi'Eo'Foci的最大值。采用拉格朗日算法（过程略）得出：是矩阵Eo'FoFo'

17、Eo的特征向量，对应的特征值为Oi2,0,是目标函数值，它要求取最大值,所以，是对应于Eo'FoFo'Eo矩阵最大特征值的单位特征向量。而另一方面,5是对应于矩阵E<）ToF（）o最大特征值0/的单位特征向量。求得和ci后，即可得到成分：tl=EOcolul=FOcl(2-17)然后,分别求Eo和Fo对h的回归方程:EO=tlpl'+ElFO=tlrl*+Fl(2-18)回归系数向量为:pl'=(tl'tl)Ttl'EOrl'=(tl'tl)Ttl'FO(2-19)3.用残差矩阵Ei和Fi取代Eo和Foo然后，求第2

18、个轴皿和c2以及第2个成分12和U2,有:t2=Elco2u2=Flc2.92=<t2,u2>co2ElFlc2(2-20)同理，有：<02是对应于矩阵EiTiFiTi最大特征值02的特征向量，C2是对应于矩阵最大特征值的特征向量。所以得到回归方程：El=t2p2+E2Fl=t2r2'+F2(2-21)回归系数向量为:CEl't2P2=rf2=些I|t2|2(2-22)4.如此计算下去，如果X的秩是A,则会有:（E0=tlpl+t2p2+tApA'（223）F0=tlrl'+t2r2，+.tArA由于ti,t2,.,tA均可以表示成Eoi,Eo

19、2,.,Eop的线性组合，因此上式可还原成yk*=Fok关于yk*=Eoj的回归方程形式，即：yk*=aklxkl*+ak2xk2*+akpxkp*+FAK（2-24）Fak是残差矩阵Fa的第K列5.确定抽取成分的个数交义有效性下面要讨论在现有的数据表下,怎么确定最好的回归方程。在许多情况下，偏最小二乘回归方程并不需要选用全部的成分h,t2,.,tA进行回归建模，而是可以像在主成分分析时一样，采用截尾的方式选择前m个成分（m<A,A=秩（X）,仅用这m个后续的成分就可以得到一个预测性能较好的模型。在多元回归分析中,经常采用抽样测试法来确定回归模型是否适于预测应用。该方法是把观测到的样本点

20、分成2部分:第1部分数据用于建立回归方程，求出回归系数估计量bB，拟合值yB以及残差均方和於；再用第2部分数据作为试验点，代入所求得的回归方程，由此求出0T和名。一般地,若有弟=扉，则回归方程会有更好的预测效果;若弟2谛，则回归方程不宜用于预测。在PLS建模中，究竟该选取多少个成分为宜，这可通过考察增加一个新的成分后，能否对模型的预测功能有明显改进来考虑。采用类似于抽样测试法的工作方式，把所有n个样本点分成2部分:第1部分除去某个样本点i的所有样本点集合（共含n-1个样本点），用这部分样本点并使用h个成分拟合一个回归方程;第二部分是把刚才被排除的样本点i代入前面拟合的回归方程，得到力在样本点i

21、上的拟合值yhj（-i）o第2.3节主元分析回归建模2.3.1主成分分析基本思想主成分分析是数学上的数据降维的一种方法。基本的想法是尽量许多原来的有一些相关的指标（如p指标），新组合的一套新的无关的综合指标来代替原来的指标。通常数学上处理是将原p做一个线性组合，作为一种新的综合指数。在所有的线性组合的F1的选择应该是最大的方差，故称F1是第一的重要组成部分。如果第一主成分不足以代表原单指标的信息，来考虑选择第二个选择F2线性组合。为了有效地反映出原始信息,F1现有信息将不会出现在F2,用数学语言表达要求Cov（Fl、F2）=0,F2即为第二主成分。以此类推可以构造出第三、第四、第p个重要的组成

22、部分（主元）。2.3.2主成分分析的数学模型设有p项指标（变量）:X,.,X,每个指标n个观测值，p得到原始数据资料阵如下：乩xipX=工21工22*2p=(X|,X2，.,Xp)(2-25)Xtn.Xttn其中Xj=(xjj,X27，xnp(i=1,2,.»p)用数据矩阵X的p个列向量(即p个指标向量)X,x2,Xp作线性组合，得综合指标向量(特征向量)：鸟=1X+2iX?+.+apiXpF、=q°X+.+%,Xp''一一(2-26)Fp=ApXi+%X2+.+cippX记为：Fj=ajjXj+32jX2+.+apjXp(i=1,2,p)每个主成分的系数平

23、方和为1。即："*+福+“；=Li=l,2,.,p(2-27)ABSTRACTInitialboilingpointisoneofstabilizedgasolineproductionfigure,whichsayshowmuchlightcomponentcontent.Inthemodernproduction,initialboilingpointisgetbyartificialinlaboratorytests.Duetothetesttimeistoolong,whichcannotbeusedinreal-timecontrol.Withthedevelopmento

24、fcomputertechnologyandcontrolmediaimprovement,thesoftmeasurementtechnologyarisesatthehistoricmoment.Thesoftmeasurementtechnologyusesstatistics,identification,artificialintelligencemethodtoestablishthemathematicalmodel,regardingthemeasurablevariableasprimaryvariable,throughmeasuringrelatedsecondaryva

25、riablepredicttheprimaryvariable.Thesoftmeasurementtechnologycanbecontinuousforecasttheunmeasurablevariable.Solvingtheproductionprocessofimportantfigureontesttimedelay,reachingthepurposeofreplacinghardwaretosoftware.Keywords：Initialboilingpoint,Softsensor,primaryvariable,Secondaryvariable(2-28)(2-29)

26、(2-30)aj为单位向量：aj'aj=1,且由下列原则决定1) Fj与Fj(详j,i,j=1,p)互不相关,主元之间无相关信息，即Cov(FpFj)=0,并有Var(Fi)=af£ai,其中£为X的协方差阵。2) F是X,X2,Xp的一切线性组合(系数满足上述要求)中方差最大的，P即Var()=maxVar(ciX.),其中c=(C|,c?,c)'oC°r=1F2与F是互不相关的X,X2，.»Xp-切线性组合中方差最大的，Fp是与F,F2,，Fp都不相关的X,X2,.»Xp的一切线性组合中方差最大的。满足上述要求的综合指标向量

27、F,F2,Fp就是主成分，这p个主成分从原始指标所提供的信息总量中所提取的信息量依次递减，每一个主成分所提取的信息量用方差来度量，主成分方差的贡献就等于原指标相关系数矩阵相应的特征值每一个主成分的组合系数如下：af=(，ii，a2iJJapj)就是特征值人i所对应的单位特征向量方差的贡献率如下:ai=&/£Ai=l%越大，说明相应的主成分反映综合信息的能力越强。2.3.3主成分分析的步骤计算协方差矩阵计算变量数据的协方差矩阵：E=(Sij)pxp，其中£(与_可(杨_弓)i，j=1,2,.,pn-1*=i求出£的特征值X及相应的特征向量a求出协方差矩阵&#

28、163;的特征值入2入22>°及相应的正交化单位特征向量:则X的第i个主成分为Fi=afX(i=1,2,p)o(1) 选择主成分在已确定的全部p个主成分中合理选择m个来确定主元。一般用方差贡献率作为评价指标。解释主成分与所含信息量的大小，主元个数的确定以累计贡献率G(z)=Z4/£&.(2-32)r=l*=1达到足够大(一般在85%以上)为原则。(2) 计算主成分得分向量计算n个样品在m个主成分上的得分：Fj=%X+6z2/.X2+.+apiXp(2,33)标准化实际应用时，指标的量纲往往不同，所以在主成分计算之前应先消除量纲的影响。消除数据的量纲有很多方法，

29、常用方法是将原始数据标准化，即做如下数据变换：=X(iXj，=1,2,.,；/=1,2,.,(2.34)Sj其中'广板厂习尸，j=1,2,po标准化后的数据阵记为X*,其中每个列向量(标准化变量)的均值为0,标准差为1,数据无量纲。标准化后变量的协方差矩阵(CovarianceMatrix)£=(Sjj)pxp，即原变量的相关系数矩阵(CorrelationMatrix)R=(r：；)nvn:IIxx-xZ(%T)-习)J7=1Jt=V/=!V/=ln-7?-li,j=1,2,.»p此时n个样品在m个主成分上的得分应为：Fj=ajX+a2jX2*+.+apjXp*j

30、=1,2,m(2-35)2.3.4主成分分析的应用因为主成分分析方法具有良好特性，主成分分析方法在建模及需要数据处理的许多领域发挥了很大的作用。主成分析方法的主要应用方向为：利用主成分分析方法来减少建模所需要的变量，降低数据空间的维数。通过分析主成分与原变量之间的系数关系，弄清原变量之间相关关系。1）多维数据的主成分分析的方法，处理后，选择前两个主成分，根据主成分所对应的矩阵特征值的大小，得出n个样品的平面分布,可直接由图形分析每个样本的主要成分的位置,也可以被分类，来样加工，也可发现的图形与大多数的样本点的离群点。2）利用主成分分析的方法对构造回归模型，是利用最主要的组成部分，取代原来的独立

31、变量作为一种新型的独立变量做回归分析。3）运用主成分分析方法进行筛选的回归变量。返回变量的选择具有重要的实际意义,为了使模型本身很容易做到结构的分析、控制和预测、构成最佳优变量子集，去选择最好的变量，使优化变量量子设置。与主成分分析，筛选变量，可以用较小的计算量选择变量，去选择最好的改变影响量子集。这也是主成分分析中的应用，最重要的是,但本文利用主成分分析法，选择最优的辅助变量集作为偏最小二乘法的输入。4）利用主成分分析方法建立综合评价函数，这种方法的主要思想是通过主成分分析的一个重要的组成部分，选择m个主成分的方差贡献作为加权函数、建立综合评价函数来进行评估。2.3.5主成分分析在本文的应用

32、利用主成分分析的方法进行辅助变量的选择。辅助变量的选择具有重要的实际意义,为了使模型本身很容易做到结构的分析、控制和预测,最好从原来的辅助变量集合中选举辅助变量形成的一个子集，构成最优的辅助变量集合。利用主成分分析筛选变量，可以用较小的计算量选择变量，去选择最好的改变影响量子集。本文在建模中辅助变量选择、也是偏最小二乘法的输入，利用主成分分析法的原理，选取最少并旦所含信息量最大的辅助变量集，作为片最小二乘的输入，以汽油初僧点为输出建立CRU主分馅塔粗汽油初馅点点软测量模型。第3章多元逐步回归在建立CRU稳定汽油初馅点软测量模型中的应用第3.1节装置简介装置设备为比利时雪夫隆炼油厂引进的二手设备

33、，原投产于一九七零年，设计能力为40万吨/年，我厂引进后由独山子石化公司设计院重新设计为32万吨/年。2000年装置进行了45万吨/年配套项目的部分改造，现加工能力达40万吨/年。装置主要分为两部分，即预加氢部分和重整部分，预加氢部分包括预加氢反应系统，精制油分僧脱水系统；重整部分包括重整反应系统、分馄稳定系统；辅助系统包括拔头油脱硫系统、分子筛干燥系统及重整加热炉发汽系统。第3.2节工艺原理3.2.1预加氢预加氢原料油为直馅僧分，杂质(硫、氮、烯炷)含量少，其反应条件较为缓和，反应在一定的氢分压和催化剂存在的情况下进行，所用催化剂为FDSTA,重整原料油在预加氢反应中主要诳行脱硫、脱氮、脱氧

34、、脱金属、烯炷饱和等反应。a.脱硫反应：为了保证重整催化剂有较好的选择性和稳定性，重整原料中应有较低的硫含量,对于PRT-C、PRT-D催化剂要求重整原料硫含量少于0.5ppmo脱硫反应速度较快，存在于原料中的硫化物很快生成烧类与硫化氢。化学反应方程式：(1) 硫醇RSH+FhRH+FhS(2) 硫酰RSR+2H2>2RH+2H2S(3) 二硫化物(RS)2+3H2>2RH+2H2Sb脱氮反应重整原料中氮化物所含的氮，经加成反应转化成氨和相应的炷，但脱氮比脱硫困难的多，多数直儒重整原料油僧分中氮化物含量少，但是由二次加工汽油作重整原料时(如焦化汽油)，其氮化物较直儒原料要高的多，需

35、经较高深度的加氢才能较彻底的脱除。对于PRT-C、PRT-D催化剂要求原料中氮含量小于0.5ppmo化学反应方程式：(1)烷基胺R-CH2-NH3+H3R-CHs+NH(2)毗嚏C5H6N+5H2C5H12+NH3(3)日引噪C8H7N+6H2C8Hi6+NH3(4)毗呢C4H5N+4H2>C4H10+NH3(5)哇I琳C11H7N+6H2C9H18+NH3c. 烯炷饱合烯炷饱合生成烷炷，其加氢反应速度比脱硫反应略慢。原料中由于烯炷的存在,会增加催化剂上的积碳，缩短生产周期。化学反应方程式：(1) 单烯炷R-CH-CH2+H2R-CH2-CH3(2) 双烯炷R-CH=CH-CH=CH+2

36、H2R-CH2CH2CH2CH3d. 脱氧反应：原料油中含氧化合物在加氢条件下与氢反应生成水及相应的炷分子，通常很容易脱除。原料中的含氧化合物主要是环烷酸，在二次加工产品中也有酚类，如不除去,进入重整反应系统后加氢生成水会使系统中存水过多，从而使催化剂减活。化学反应方程式：(1) 酚类C6H6O+H2C6H8+H2。(2) 环烷酸C7H12O2+3H2C7H14+2H2Oe. 脱卤素反应有机卤素物经加氢可以生成相应的炷类及卤化氢，有机卤化物的脱卤比脱硫更困难一些。化学反应方程式：RCL+H2一RH+HCLf. 脱金屈反应原料中含砰、铜、铅的化合物，在加氢条件下分解生成金属，然后吸附在加氢催化剂

37、上，如不除去进入重整反应系统会使重整催化剂减活或失去活性。3.2.2重整催化重整以一定沸点范围的石油馅份（属于燃料型为80-170°C或180C,化工型装置为70-140°C）,经重整催化剂的作用，并在氢气保护和一定的反应条件下，发生芳构化反应和其它反应，使原来含少量芳炷的原料发生分子结构的重排，而成为富含芳炷和异构烷炷的生成物。这种产物经进一步加工处理，可获得高纯度的芳烷或高辛烷值的汽油。催化重整主要进行下述反应a.六员环烷炷脱氢：在所有重整反应中，六员环烷炷脱氢反应速度最快，而且能充分转化成芳炷，是重整的最基本反应。C7H|4=C7H8+3田AHARON(3-1)(RO

38、N74.8)(RON120)+205.8KJ/mol45.2b.五员环烷炷异构化成六员环烷炷：C7H14=C7H14AHARON(3-2)(RON80.6)(RON74.8)-18.9KJ/moI-5.8c.加氢裂化反应:C=C=C=C=C=C=C+H_C=C=C+C=C-(C)2AHRON(3-5)(RONO)(RON102)-58KJ/mol102d.脱甲基反应：C6H14+H2C5H12AHARON(3-6)(RON24.8)(RON61.8)-67.6KJ/mol37.9C7H8+H2C&H6AHARON(3-7)(RON120)(RON100)-54.7KJ/mol-20e.

39、芳炷脱烷基反应:CgHio+H?>C7H8+CH4AHARON(38)(RON120)(RON120)-56.8KJ/mol0f.积炭反应：炷类的深度脱氢，生成烯炷和二烯炷，烯炷进一步聚合及环化，形成稠环芳炷,吸附在催化剂上，最终转化成积炭，而使催化剂失活。上述反应中芳构化反应和异构化反应是所希望的，而加氢裂解反应、脱甲基反应和积炭反应是不希望的，应尽量减小或避免。第3.3节工艺流程说明3.3.1预加氢部分经脱神处理后的重整原料（砰V150ppb,馅程：干点V175°C）自原料罐区来，经原料泵P101升压后与预加氢循环压缩机C-101来的氢气混合后进入立换E-101的管程、卧式

40、换热器E-101A/B的壳程与反应产物换热到200C左右进入原料加热炉F-101加热，加热到280C左右的物料依次通过加氢反应器R101、脱氯器R-102进行反应，其中在R-101中进行原料的脱硫、脱氮、脱氧、烯炷饱合以及脱金属反应；在R-102中进行脱氯反应（预加氢反应如图3.3.1所示）。反应后的产物进入E101A/B的管程、立换E-101的壳程与预加氢原料换热降温到100°C左右进入空冷E102、水冷E-103冷却至小于40C后进入预加氢高分罐V-101进行油气分离。V101压力控制在2.0MPa,含氢气体去柴油加氢装置或去瓦斯系统；如氢油比不足一部分氢气可返回压缩机入口分液罐

41、V104进行循环。液相产物E104管程与脱水塔T-101底精制油换热后进T-101进行脱水；T-101还有一部分进料来自馅分油加氢裂化的重石脑油（进料图如图3.3.2所示），由P-201C/D升压经T-101垫油线送至E-104管程入口，同V-101减油一起进入T-101（脱水塔如图3.3.3所示）。V-101切水包中的含硫污水进入下水。T101塔底油一部分经P-102升压后进再沸炉F102加热到215°C左右返塔，为T101提供热量；另一部分经E104壳程与V-10I减油换热到80°C左右后去重整作重整原料。塔底部分物料也可经E-104壳程、E301壳程冷却，作为精制油储

42、存于906#、907#罐。塔顶产物经E105,E106冷却至40°C后进入回流罐V103进行汽液分离，罐内液体经回流泵P-103升压，一部分做回流打入T-101的顶部，一部分经脱硫罐脱硫后作为乙烯原料外送或去加氢T-403脱硫做石脑油送原料车间，顶部气体进入瓦斯系统。a>i5町视剧U11201：顼加氢反应器R101-103P11142.101MPaP11142.101MPaTH08304rPHWOOMMPiT114<R1010420RIOtSftYJ10150R101EHYJ101T1149TI126C184rT1143»72rTI144287.7T?TI144

43、287.7T?T1145rT1145PD10400.014MPamie2.018E101ST101ItC202AC202BI喊合通4RC201IV301R401T101ItC202AC202BI喊合通4RC201IV301R401图31预加氢反应器鸟15Kia£iU11201:F-101U11201:F-101P»W1O1T51RXL2CT2GM.O%FC104AGMPaPB32C144ITI5O2CarTW03G650PT15O4CS25TTI5O5G0E101E*12.095MPaTH0402M4rFC101GC0213MPaTC101OA28OTE102ABCMPaA

44、痂4"M02GT101P101CE101AACWTI129G264,rFS122G”1228PH220OAMMR10Y81011664FO1140331406Nm30Nm3hStKS1013A5R102ORPD104O0013TH13。PIW3G1177PaFI10TG7821Nm3hFcioecvW101GfFC151RWlMPD18G0.420MPaloicjE10-WMAT101KFI10TGA1Mt4Nm3hE10241.1CPt2242.758MPaF1105V2CMT11050420VFC102CA42184k°hA581rTI1MTill®T1151

45、A79Qr1680cF1012124r1M理1P1Q,T1151BiHCfiOIGM50J%283.1T?nisic2849rnisiDni5iEMl7rT1181FF0191G0：14M190Nm34T8Nma*hLX101C202A1C202B1整合通It4RC»11V301R401T402图32预加氢进料图32预加氢进料immamimmamU11201:T-101E1O0ABCFC123GFC10TCV10诔WIWG81MPa36GF102AF102BP103BW17O200.9CT1117AOEQT?。物MPaFO12CO0169601Wkf2044kghbimbo210.4

46、TFC1W0M50235kghPC102A0.8&8MPaFQ11BOl01TM2eNm323Hm3hU111OMSS1X102.WO2F01120_0531517Na3100Nm3hF001PCI100A0225MP，FO113O0283X14151Nm3HWW3(jF1123H1WBC2138Pmor种0PI1090图3-3脱水塔T-101TC1ZSA123.8TE1MFC1WOM50108kflhFC110CMA6413kghTilFC1110M0182«kaliT1132C2M4T：r102MBW104GJ1T胶|St&SiUTC11SCHC102GMO.OS

47、OPC108CTI134AC210TY14RT201C202AC201V301U合谓F301R401怦噩叶IFI111G0341MP1NB*心Qiao05075378kg4382S53.3.2重整部分预加氢来的精制油经重整原料泵P-201A/B升压后与重整循环压缩机C-201来的氢气混合，然后进入原料换热器E-201管程与重整反应产物换热到420°C左右，依次经过F201,R201,F202,R202,F203,R203,F204,R204完成重整反应。因为反应过程主要是吸热反应，反应中会有温降产生，其中R-20I的温降约为75°C,R202的反应温降约为40°C

48、,R203的反应温降约为25°C,R204的反应温降约在10-20°C（重整反应器如图3.3.4）o反应后的产物进入E201壳程与重整原料换热到120°C左右，然后依次通过空冷E203、水冷E204壳程，冷却到30-40°C后进入气液分离罐容V201oV201控制压力在1.25MPa,顶部氢气进C-201进行压缩循环；C-201出口氢气一部分作重整循环氢用，一部分经水冷器E211壳程冷却后进V104作为预加氢的新氢（气液分离如图3.3.5所示）：V201压控阀出来的气体进燃料系统；V201底的液体经E205A/B/C壳程与稳定塔T201底稳定汽油换热到1

49、20C左右进T201；塔底稳定汽油经E205A/B/C管程与V-201减油换热到90°C左右，再经E209、E210冷却至小于40°C后外送。塔底设有以4.0MPa蒸汽为热源的再沸器E208为T201提供热量。T201顶出来的产物经空冷E206A/B,水冷换E207冷却到小于40°C后进入回流罐V202,气相经压控进燃料系统，液相经回流泵P203升压，一部分打回流返塔，一部分作为液态炷外送出装置。为更好的保护重整催化剂少受硫的污染，在蒸发脱水塔旁设有精制油脱硫罐V-102,若预加氢反应正常，V-102也可不投用。为防止铉盐的沉积，在E101管程出口上设有自加氢凝结

50、水泵P405来的凝结水注水点。为提高T101的精饲效果，C201的出口线上设有条到T101底的氢气汽提线。设置有注氯泵P-303,抽二氯乙烷注入一反和三反入口，以保持催化剂的水氯平衡。为抑制新鲜或再生催化剂的裂解活性，在重整各个反应器入口均设有预硫化设施。用高压氮气将注硫罐V-310/1,2,3,4内的硫化物压入各个反应器入口进行硫化。为脱除和降低系统存水，重整临氢系统设有分子筛干燥系统及分子筛再生系统。为更好的再生重整催化剂,C201出口及R202出口均设有工业风注入点（稳定汽油蒸僧如图3.3.6所示）。前言1第1章绪论1第1.1节基于传统方法的软测量建模2第1.2节基于统计回归分析的软测量

51、建模3第1.3节人工智能建模方法4第1.4节软测量方法在工业过程中的实现61.4.1软测量辅助变量的选取71.4.2软测量的数据选择与处理81.4.3软测量模型的辨识和在线校正8第1.5节文章的研究内容和结构安排91.5.1文章研究的主要内容91.5.2文章的结构安排10第2章基于统计学的回归建模方法10第2.1节多元逐步回归方法11第2.2节偏最小二乘方法132.2.1偏最小二乘法的特点132.2.2偏最小二乘法的建模步骤14第2.3节主元分析建模162.3.1主成分分析基本思想162.3.2主成分分析的数学模型172.3.3主成分分析的步骤182.3.4主成分分析的应用202.3.5主成分

52、分析在本文的应用20第3章多元逐步回归在建立CRU稳定汽油初僧点软测量模型中的应用21第3.1节装置简介21H24GLC2O7GFC2O2CAFO258O14244Nr0Nm3hU11201：V-201.V-204E204ASV101|C202V202LSM2«>GBLSK2030C5*晦纹空虹登的气吁芽故LC2O2GA50.0、C201E201/2PI1O4G2.754MP>S213LC20tGA502FI107CTVHNm»hP101FI107GA10558Nm3hPC20SOXidOiMFTC202T2011C202A1F101岐含姻道ITWIC201V3

53、01R401图3-5气液分离罐U11201:T-201TK21SBC2205rFC203OABIM232142hLC2CCCAO2、FC204GBM232S8kahTlZiaAO221-3VFC2O4GC42SB8k.hRXL1Q28.6%AI1023-M%Tt205C"rT1216O<44rPCEGA1.000MP1-rtrafji研一T1245蚯汽*«mi7.G22S4FCZWGAT1257O124.3rTH”682T?HC2O4CM0.090MPaLoraIA(W.1、TC222GA221.4T：C202APS217CPS217CAPttlTG0*2MPa210

54、.7VTI247210rIAC2O3G(M04GF020SGO3M63t4Nm3160Naai)FO222Q0:224472.1t449th图3-6稳定汽油分馅塔第3.4节MSR建模3.4.1辅助变量的选择根据机理对数据进行筛选：汽油精德发生在塔T201,故与初僧点最相关的因素主要为分馅塔的操作参数，主要有塔顶温度、塔底温度、液位、压力、进料量及进料温度、回流量及回流温度等、由于分馅塔的压力在实际操作中，保持某一特定的值，且一般不作为调节手段，所以压力与初馋点的关系没有操作数据作为参考依据。由此剔除以下数据：2预加氢进料量，64T101进料温度,65T101塔顶温度，68T101回流流量，73

55、T101塔底温度,164段混氢流量，165二段混氢流量，176R201入口温度，189R202入口温度，204R203入曰温度,216R204入曰温度，289V201温度,290V201压力，295V201液位，317V205压力，368V202进料温度,369V202液位。剩下F205出口温度X20,T201塔顶回流X22。（声明一点这两个变量是和现场操作工人讨论后得出的）作为辅助变量建模。3.4.2数据处理在一组条件完全相同的重复试验中，个别的测量值可能会出现异常。如果测量值过大或过小，这些测量数据不不正常的，或称为可疑的。对于这些可以数据应该使用树立统计的方法辨别其真伪，并决定取舍。常用的方法有拉一达法、小维那特法、格拉布斯法。当试验次数较多时，可简单地用3倍标准偏差（3S）作为确定可疑数据取舍的标准。当某一测量数据（xi）与其测量结果的算术平均值（又）之差大于3倍标准偏差时，用公式表示为：Ixi一又I>3S（3-9）则该测量数据应舍弃。这是美国混凝土标准中所采用的方法，由于该方法是以3倍标准偏差作为判别标准，所以亦称3倍标准偏差法，简称3S法。取3S的理由是:根据随机变量的正态分布规律，在多次试验中，测量值落在X-，一3S与X十3S之间