CN101169387B - 一种在线确定常压塔顶石脑油质量指标的软测量方法 - Google Patents
一种在线确定常压塔顶石脑油质量指标的软测量方法 Download PDFInfo
- Publication number
- CN101169387B CN101169387B CN2007101711161A CN200710171116A CN101169387B CN 101169387 B CN101169387 B CN 101169387B CN 2007101711161 A CN2007101711161 A CN 2007101711161A CN 200710171116 A CN200710171116 A CN 200710171116A CN 101169387 B CN101169387 B CN 101169387B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- atmospheric
- naphtha
- tower
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 10
- 230000008016 vaporization Effects 0.000 claims abstract description 6
- 230000001537 neural effect Effects 0.000 claims abstract 2
- 238000010992 reflux Methods 0.000 claims description 28
- 238000004519 manufacturing process Methods 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 11
- 239000013598 vector Substances 0.000 claims description 11
- 238000010606 normalization Methods 0.000 claims description 10
- 230000001419 dependent effect Effects 0.000 claims description 9
- 238000012821 model calculation Methods 0.000 claims description 9
- 238000004821 distillation Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 5
- 238000003062 neural network model Methods 0.000 claims description 5
- 238000009834 vaporization Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 8
- 230000004907 flux Effects 0.000 abstract 5
- 238000012706 support-vector machine Methods 0.000 abstract 1
- 239000003921 oil Substances 0.000 description 27
- 239000010779 crude oil Substances 0.000 description 13
- 239000000047 product Substances 0.000 description 11
- 239000003208 petroleum Substances 0.000 description 6
- 238000009835 boiling Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 239000012071 phase Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 239000003502 gasoline Substances 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000009897 systematic effect Effects 0.000 description 3
- 238000005292 vacuum distillation Methods 0.000 description 3
- 238000010420 art technique Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000002283 diesel fuel Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005504 petroleum refining Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000007670 refining Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 239000006227 byproduct Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000005194 fractionation Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000010946 mechanistic model Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 238000001577 simple distillation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000004148 unit process Methods 0.000 description 1
Images
Landscapes
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种在线确定常压塔顶石脑油质量(干点)指标的软仪表技术:选取常压炉出口温度、常压塔处理量、塔顶温度、塔顶压力、顶回流单位处理量温差、常压塔顶石脑油流量、常一线流量、常二线流量、常三线流量、常压塔顶循单位处理量温差、常一中单位处理量温差、常二中单位处理量温差、汽化段温度、汽提蒸汽流量、以及前一时刻常压塔顶石脑油干点人工分析值,作为常压塔顶石脑油干点软仪表的输入变量;采集工业装置样本数据,通过自适应谐振神经网络-支持向量机回归组合算法建立常压塔顶石脑油干点模型,实现在线计算;通过模型系统偏差的分析,对模型计算值进行加权修正,在线确定常压塔顶石脑油干点。
Description
技术领域
本发明属于石油炼制与过程控制交叉领域,涉及石油炼制常减压生产装置中常压塔顶产品——石脑油质量(干点)指标的在线确定方法。
背景技术
由于世界能源问题的日趋严重,原油市场和成品油市场变化莫测,原油的品质和价格成为世界炼油工业选择炼制原油的重要因素。常减压装置是整个炼油工业的首要生产环节,其中初馏塔、常压塔和减压塔作为常减压装置的关键设备,负责完成原油一次加工的主要过程。常减压装置直接处理原油,将原油切割成各种不同馏分的产品并消耗大量的能量。这些馏分产品或作为后续工段的进料,或经调和后作为成品油出售,它们有一定的质量要求。其中常压塔是常减压装置的核心精馏塔。
某一典型的常压塔流程如图1所示,初馏塔底油经常压炉加热后进入常压塔第57层塔盘上。常压塔顶油气换热后进入到常压塔顶回流罐,液相经常压塔顶回流泵一部分打至常压塔顶第一层塔盘上作塔顶回流,剩余部分进重石脑油出装置。未凝油气经常压塔顶空冷器二级冷却后进入常压塔顶产品罐,不凝气可去一级减压炉作燃料,紧急情况也可放火炬。常压塔顶油(即石脑油)经常压塔顶产品油泵可送去脱丁烷塔进料缓冲罐,或者去初馏塔顶回流罐,也可以直接出装置。常压塔设置有三个侧线:常一线自第28层塔盘流出,进入常压汽提塔上段,油气返回常压塔第28层塔盘上方,常一线油由常一线泵抽出,经冷却后出装置;常二线自常压塔第40层塔盘流出,进入常压汽提塔的中段,油气返回到常压塔第40层塔盘的上方,常二线油由常二线泵抽出,经冷却后出装置;常三线自常压塔第50层塔盘上流出,进入常压汽提塔的下段,油气返回到常压塔第50层塔盘的上方,常三线油由常三线泵抽出,经冷却后出装置。常压塔设置3个循环回流:顶循环回流油经常压塔顶循回流泵从第6层塔盘抽出,换热后返回到第3层塔盘上;常一中回流由常一中泵从第32层塔盘抽出,经换热冷却后,与初侧线合并后返回第29层塔盘;常二中回流由常二中回流泵从第44层塔盘抽出,经换热冷却后返回至第41层塔盘。常压塔底油则由常底泵抽出,分六路进入一级减压炉,经加热后进入一级减压塔。
石油馏分是复杂的混合物,一般常用沸点范围来表征其蒸发及汽化性能。沸点范围又称沸程。石油馏分的沸程测定是一种在标准设备中,按照GB6536-86规定的方法进行的简单蒸馏。国外将此类方法成为ASTM(American Society forTesting Material,美国材料试验学会)蒸馏或恩氏(Engler)蒸馏。其测定过程是,将100mL(20℃下)油品放入标准的蒸馏瓶中按规定的速度进行加热,其馏出第一滴冷凝液时的气相温度称为初馏点。随后,其温度逐步升高而不断地馏出,依次记下馏出液达10mL、20mL直至90mL时的气相温度,称为10%、20%、......、90%馏出温度。当气相温度升到一定数值后,它就不再上升反而回落,这个最高的气相温度称为干点(或终馏点)。
常压塔顶石脑油干点是常压塔重要的控制指标,其控制的好坏,直接影响轻质油收率的提高和后续工序的操作稳定。由于在常减压装置生产过程没有合适的在线仪表能对常压塔顶石脑油干点进行实时分析,一般操作人员仅能根据每4或8个小时分析样的分析结果对生产进行调整,而这些分析样需要2小时的分析时间,这些因素都不可避免的带来了分析结果的太滞后和生产操作条件不能及时调整,从而给产品质量的控制带来困难。
近年来软仪表技术在工业中获得了大量成功的应用,解决了许多不可测关键控制指标的“测量”问题。软仪表技术是通过可测的辅助变量建立关于不可测主导变量的模型,从而在线估计出实时连续的质量指标估计值。软仪表技术常用建立数学模型的方法,主要采用机理建模和统计建模。石油馏分干点的在线预测机理模型方法如:中国专利(ZL95101183.9)公开了一种用于石油分馏塔塔顶汽油干点和中部侧线产品柴油95%点的在线确定方法。利用温度、油汽分压和油品组分分率,以及分馏塔及其相关设备的动态模型,在线实时计算分馏塔内物流流量、内回流流量和油品组分分率,进而随时确定汽油干点和柴油95%点。其他一些机理模型方法,主要采用实测的分馏塔温度和压力,按稳态关系估计汽油或柴油产品质量,如:Applied Automation INC.(A.A.公司),兰州炼油化工总厂自动化所,上海高桥石化公司炼油厂等等。同时,大量公开的文献也报道各种统计方法建立馏分油质量指标模型,主要有:各种全局多元线性回归方法,它们具有简洁明确的解析表达形式,但难以描述石油分馏这样高度非线性系统。有多种非线性建模方法(如神经网络方法、支持向量基回归方法、关联向量基回归方法、核PCA回归方法等等),它们可以处理高度非线性体系,建模时有较高的拟合精度,但模型预测性能依赖于输入变量的选择以及建模样本的质量。
现有的各种在线确定馏分油产品质量指标方法存在如下的一些缺点:
(1)石油分馏塔处理的油品组成复杂、经常变换且油品性质通常难以及时获得,操作变量多、存在交互作用且对产品质量影响复杂,准确的机理模型难以及时建立。
(2)影响馏分油质量的操作变量众多且存在交互作用,统计模型的输入变量通常或存在大量冗余变量,使建立的模型稳定性差;或过于精简,难以描述馏分油质量指标的变化规律。
(3)由于炼制原油频繁变换且原油性质无法及时获得,因此大量统计模型方法未考虑炼制原油性质的变化,造成模型精度下降、甚至严重偏离采样分析值。
同时,常减压装置通常根据生产负荷和产品质量需求的改变、以及炼制原油性质的变化等进行操作条件的调整,从而积累大量具有一定类别特征的生产数据。针对常压塔呈高度非线性、且具有一定类别特性等特征,本发明采用了基于自适应谐振神经网络(Adaptive Resonance Theory,ART)将样本模式分割成若干模式特性相近的子空间;然后,对各子空间分别采用支持向量回归(Support VectorRegression,SVR)建立各自模型,实现基于样本模式空间分割的“分段”建模,降低各子空间模型的非线性程度,提高模型预测性能。其中ART网络具有自组织、自适应分类功能;SVR是一种基于结构风险最小化原则下的机器学习方法,具有强非线性表达能力,且建立模型具有很好的泛化能力,在过程建模中得到了广泛的应用。
发明内容
本发明目的是提供一种在线确定常压塔顶石脑油干点的软仪表技术。选取常压炉出口温度(x1,℃)、常压塔处理量(x2,T/Hr)、塔顶温度(x3,℃)、塔顶压力(x4,MPa)、顶回流单位处理量温差(x5,℃)、常压塔顶石脑油流量(x6,T/Hr)、常一线流量(x7,T/Hr)、常二线流量(x8,T/Hr)、常三线流量(x9,T/Hr)、常压塔顶循单位处理量温差(x10,℃)、常一中单位处理量温差(x11,℃)、常二中单位处理量温差(x12,℃)、汽化段温度(x13,℃)、汽提蒸汽流量(x14,Kg/Hr)、以及前一时刻常压塔顶石脑油干点人工分析值(x15,℃),作为常压塔顶石脑油干点软仪表的输入变量;利用常压塔相关的测量仪表值,直接测量或间接计算获得x1~x14的实时测量值,以及通过生产装置集散控制系统的人机界面获得前一时刻常压塔顶石脑油干点人工分析值x15;基于软仪表,在线确定常压塔顶石脑油干点。
1,输入变量的选取(或构造)
在实际生产过程中,影响常压塔顶石脑油干点的主要因素:常压炉出口温度、常压塔处理量、塔顶温度、塔顶压力、各中段循环以及顶回流带出能量、常压塔顶石脑油流量以及各侧线流量、汽化段温度、汽提蒸汽流量以及炼制原油的性质。其中各中段循环以及顶回流带出能量的计算,涉及流股物性数据,因此无法在线获得;炼制原油的性质经常无法及时获得;其他影响因素都有测量仪表直接获得或间接计算获得。
本发明的特点是:
(1)构造了顶回流单位处理量温差、常压塔顶循单位处理量温差、常一中单位处理量温差、常二中单位处理量温差的关系式,分别表征了常压塔在单位处理量下顶回流带出能量、常压塔顶循带出能量、常一中带出能量、以及常二中带出能量。
(2)将前一时刻常压塔顶石脑油干点人工分析值作为炼制原油性质的间接表征。
(3)通过ART网络与SVR建模方法相结合,实现基于建模空间分割的分段建模。
为此,常压塔顶石脑油干点软仪表的输入变量选取如下:
(1)常压炉出口温度(x1,℃)
(2)常压塔处理量(x2,T/Hr)
(3)塔顶温度(x3,℃)
(4)塔顶压力(x4,MPa)
(5)顶回流单位处理量温差(x5,℃)
(6)常压塔顶石脑油流量(x6,T/Hr)
(7)常一线流量(x7,T/Hr)
(8)常二线流量(x8,T/Hr)
(9)常三线流量(x9,T/Hr)
(10)常压塔顶循单位处理量温差(x10,℃)
(11)常一中单位处理量温差(x11,℃)
(12)常二中单位处理量温差(x12,℃)
(1 3)汽化段温度(x13,℃)
(14)汽提蒸汽流量(x14,Kg/Hr)
(15)前一时刻常压塔顶石脑油干点人工分析值(x15,℃)
其中,顶回流单位处理量温差、常压塔顶循单位处理量温差、常一中单位处理量温差、常二中单位处理量温差定义(或构造)如下:
式中t1:塔顶回流温度(℃);m0:常压塔顶回流量(T/Hr);m1:常压塔顶循回流量(T/Hr);m2:常一中回流量(T/Hr);m3:常二中回流量(T/Hr);Δt1:常压塔顶循抽出与返塔温差(℃);Δt2:常一中抽出与返塔温差(℃);Δt1:常二中抽出与返塔温差(℃)。
x1、x2、x3、x4、x6、x7、x8、x9、x13、x14、t1、m0、m1、m2、m3、Δt1、Δt2、Δt1均可以由常压塔的相关测量仪表直接获得。x5、x10、x11、x12可以由公式(1)实时间接计算获得。前一时刻常压塔顶石脑油干点人工分析值x15,通过生产装置集散控制系统的人机界面,由操作人员输入获得。因此,常压塔顶石脑油干点软仪表的所有输入变量都可以实时、在线获取。
2.建模样本的预处理
为了消除量纲的影响,对采集的样本数据进行归一化预处理。输入变量利用下式(2)进行归一化处理:
输出变量利用下式(3)进行归一化处理:
(3)式中,sy表示因变量归一化处理后的值,[ymin,ymax]表示因变量的变化范围,归一化处理后因变量的变化范围为[a,b]。
采集到n组代表性的工业装置数据,其中每组数据包含输入变量(x1、x2、x3、x4、x5、x6、x7、x8、.......、x15)与对应的常压塔顶石脑油干点(y),经(2)、(3)式归一化处理后为[sx1,sx2,…,sx15,sy],形成建模样本。
3.基于ART-SVR的常压塔顶石脑油干点模型
设对于样本容量为n的建模样本,采用ART技术,自适应地将其分割成m组模式特性相近的子空间(这样分割可以降低各子空间的模型的非线性程度,提高软仪表模型预测精度与稳定性),子空间建模样本容量分别为n1,n2,…,nm 与子空间相对应的ART网络权向量分别为w1,w2,…,wm。
设与第i个子空间相对应的ART网络权向量为:
wi=[w1,w2,…,w15]i。
设第i个子空间的建模样本为:
[sx1,sx2,…,sx15,sy]1 i、[sx1,sx2,…,sx15,sy]2 i、......、[sx1,sx2,…,sx15,sy]ni i。
基于第i个子空间的建模样本,采用SVR技术建立第i个子空间的模型。设为:
sy=fi(sx1,sx2,…,sx15)
设sx=[sx1,sx2,…,sx5],
则,
其中, 是RBF径向基核函数;σ2是RBF径向基核函数的宽度;αj≥0,j=1,…,ni,为拉格朗日乘子;bi是偏置值;fi(·)为第i个子空间[sx1,sx2,…,sx15]与sy之间的SVR模型。
则,所有m个子空间的关联模型f1(·),f2(·),…,fm(·),形成整个常压塔顶石脑油干点模型,且模型可以表示为:
sy=f(sx1,sx2,…,sx15)=fi(sx1,sx2,…,sx15)假如
其中,
sx=[sx1,sx2,…,sx15], 表示输入变量组成的向量;sx=[sx1,sx2,…,sx15]与第i个子空间相对应的ART网络权向量最相近。
4.模型计算值系统偏差的加权修正
该模型可以投运在生产装置的计算机集散控制系统(Distributed controlsystem,DCS)中,通过以上15个输入变量,在线计算常压塔顶石脑油干点。其采样时刻的模型计算值可以与该时刻的人工分析值进行对比,并对计算系统偏差进行加权修正。
则,依次分析到前m时刻的计算偏差为:
5.常压塔顶石脑油干点软仪表的在线计算
附图说明
图1常压塔工艺流程图;
图2常压塔顶石脑油干点软仪表的在线计算流程
具体实施方式
通过以下实施例的说明将有助于理解本发明,但并不限制本发明的内容。
实施例1
1.预处理样本
采集214组样本数据,利用(2)式,对上述各自变量进行归一化处理:x1的变化范围[325,336],x2的变化范围[344,496],x3的变化范围[123,139],x4的变化范围[0.02,0.06],x5的变化范围[0.49528,6.6955],x6的变化范围[6,67],x7的变化范围[16,27],x8的变化范围[24,49],x9的变化范围[21,42],x10的变化范围[3.3714,15.932],x11的变化范围[7.288,14.962],x12的变化范围[13.227,41.771],x13的变化范围[323.5,329],x14的变化范围[0.97,1.785],x15的变化范围[156.5,183],[a,b]的取值通常为[0,1]、[0.1,0.9]、[0.2,0.8]等值,此处取a=0,b=1,进行归一化计算:
利用(3)式,对上述因变量进行归一化处理:y的变化范围[160,180],取a=0,b=1,进行归一化计算:
2.基于ART-SVR的常压塔顶石脑油干点模型
采用ART技术,通过自组织聚类将输入样本分成2个子空间,在每个子空间内的输入样本个数分别为144、70。其中ART网络权向量分别为:
w1=[0.1119 0.2443 0.3261 0.4446 0.5129 0.3335 0.4575 0.4068 0.3558 0.3217 0.3506 0.273 0.4419 0.4076 0.4773]T
w2=0.0935 0.5934 0.3839 0.3127 0.3576 0.3952 0.7675 0.6695 0.7917 0.511 0.45634 0.2849 0.4974 0.2842 0.6809]T
基于第1个子空间的建模样本,采用SVR技术建立第1个子空间的模型,其具体的模型参数如下:
(1)RBF径向基核函数的宽度σ2=29.932
(2)拉格朗日乘子α
α1~α144分别为:
1.4573,-1.6002,1.8153,0.1408,-0.0007,-0.9855,0.4812,2.5634,1.3435,-0.867,-0.5571,1.983,-2.2542,-0.9366,-3.7789,1.2029,0.6732,-1.1335,1.438,2.0002,1.6584,0.6511,-4.043,1.718,1.8181,-2.5,1.4361,0.7836,0.4084,1.4349,3.1544,0.9104,-2.2101,0.1726,0.8672,0.1032,-0.0998,0.9106,-1.2496,-1.4905,1.1103,-2.8416,0.489,-0.1718,1.4444,2.6551,-3.9308,0.0841,1.2677,-0.5811,0.1061,2.0973,0.8195,-0.2483,0.1867,-2.2516,-0.5896,0.93,1.3695,-0.4883,0.3457,0.387,-1.6526,-0.7869,-1.4065,0.0274,3.2889,-0.9519,3.134,0.2523,-2.8262,-0.6062,0.5301,2.4521,-0.8114,-1.5919,1.0099,0.043,-2.1753,-0.3385,1.125,-1.1376,1.9628,-0.0297,-1.7062,0.8358,0.4288,2.1005,-0.0123,0.6933,-0.1881,-0.1368,-0.2693,-0.9212,-2.7911,-2.2705,-0.1995,-1.3494,-2.2154,-1.2452,2.9256,3.2517,0.3455,1.4421,0.3933,-1.0607,1.8145,-2.0639,-2.2807,-2.4791,2.0502,1.294,3.2381,-3.6608,1.3081,3.0046,-3.6628,1.0179,-1.7468,0.044,0.5745,0.4211,-0.4587,-2.6081,-0.9536,-1.8756,1.7219,-0.2548,-0.3253,0.9397,-1.3933,0.2054,3.1652,0.2519,-3.0052,0.118,-0.7773,0.5745,1.0878,-2.046,-1.2502,0.6056,-0.8753,1.1105
(3)偏置值b1=-0.151
所得模型为:
基于第2个子空间的建模样本,采用SVR技术建立第2个子空间的模型,其具体的模型参数如下:
(1)RBF径向基核函数的宽度σ2=72.235
(2)拉格朗日乘子α
α1~α70分别为:
-3.7343,4.0046,-0.1908,3.3423,2.8114,2.9569,1.4128,-2.9388,1.7654,-0.0902,-4.046,-0.6466,4.5443,-1.3812,5.0608,-1.2393,4.5572,1.6977,-4.9586,1.4616,-3.6372,-2.6511,-1.2519,-0.0001,1.3424,-2.4017,-1.9785,2.5433,3.191,1.4485,1.5652,-2.1666,1.1395,2.886,0.2472,-0.9538,3.104,0.6747,-2.3058,1.5871,-0.4164,0.7716,-4.0417,-3.4048,-3.2497,0.1065,-5.9586,0.6434,1.0818,3.3086,-2.3385,0.8748,-0.8998,1.3952,0.8945,-3.1362,-1.461,1.4207,0.0064,-2.9581,3.0328,-2.6056,-4.5223,1.1028,0.144,4.3972,-0.3444,2.9789,-1.9382,-1.6554
(3)偏置值b2=0.4387
所得模型为:
则,对于输入变量组成矢量sx=[sx1,sx2,…,sx15],常压塔顶石脑油干点神经网络模型为:
上述描述了,如何基于工业装置生产过程数据建立常压塔顶石脑油干点神经网络模型。以下将描述,如何对常压塔顶石脑油干点神经网络模型的计算值进行系统偏差的加权修正。
设当前时刻为t,15个输入变量的采集值为:
经归一化计算得到:
判断 属于哪一个子空间的模型:
设,依次分析到前m=6时刻的计算偏差;在实际投运时,设t时刻它们的偏差为:
[Δyt-1Δyt-2…Δyt-6]=[0.3-0.23-0.45 0.61-0.73-0.82]
设,权值取为:
[w1w2…w6]=[0.3 0.2 0.2 0.1 0.1 0.1]
在实际应用中,常压塔顶石脑油干点软仪表的在线精度可达到相对误差5%左右。
Claims (2)
1.一种在线确定常压塔顶石脑油质量指标的软测量方法,其特征在于,包括如下步骤:
a、选取常压炉出口温度(x1,℃)、常压塔处理量(x2,T/Hr)、塔顶温度(x3,℃)、塔顶压力(x4,MPa)、顶回流单位处理量温差(x5,℃)、常压塔顶石脑油流量(x6,T/Hr)、常一线流量(x7,T/Hr)、常二线流量(x8,T/Hr)、常三线流量(x9,T/Hr)、常压塔顶单位处理量温差(x10,℃)、常一中单位处理量温差(x11,℃)、常二中单位处理量温差(x12,℃)、汽化段温度(x13,℃)、汽提蒸汽流量(x14,Kg/Hr)以及前一时刻常压塔顶石脑油干点人工分析值(x15,℃),作为常压塔顶石脑油干点软仪表的输入变量;
b、利用与常压塔相关的测量仪表值和生产装置集散控制系统的人机界面,直接测量或间接计算获得步骤a所述的x1~x15的实时测量值;
c、根据步骤b得到的所述实时测量值,基于软仪表在线计算因变量常压塔顶石脑油干点(y,℃);
所述x1、x2、x3、x4、x6、x7、x8、x9、x13、x14均可以由初馏塔的相关测量仪表直接获得;所述x15通过生产装置集散控制系统的人机界面,由操作人员输入获得;
所述顶回流单位处理量温差、常压塔顶单位处理量温差、常一中单位处理量温差、常二中单位处理量温差由下式计算后间接获得:
式中t1:塔顶回流温度(℃);m0:常压塔顶回流量(T/Hr);m1:常压塔顶循回流量(T/Hr);m2:常一中回流量(T/Hr);m3:常二中回流量(T/Hr);Δt1:常压塔顶循抽出与返塔温差(℃);Δt2:常一中抽出与返塔温差(℃);Δt1:常二中抽出与返塔温差(℃);
输入变量利用下式进行归一化处理:
常压塔顶石脑油干点(y,℃),即因变量,利用下式进行归一化处理:
式中,sy表示所述因变量归一化处理后的值,[ymin,ymax]表示所述因变量的变化范围,归一化处理后的变化范围为[a,b];
采集工业装置数据,经归一化处理后,采用自适应谐振神经网络-支持向量回归建立[sx1,sx2,…,sx15]与sy的关联模型:
sx=[sx1,sx2,…,sx15],表示所述输入变量组成的向量;sx=[sx1,sx2,…,sx15]与第i个子空间相对应的ART网络权向量最相近;则通过对所述sy的反归一化,就可以求得常压塔顶石脑油干点的模型计算值即
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101711161A CN101169387B (zh) | 2007-11-28 | 2007-11-28 | 一种在线确定常压塔顶石脑油质量指标的软测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101711161A CN101169387B (zh) | 2007-11-28 | 2007-11-28 | 一种在线确定常压塔顶石脑油质量指标的软测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101169387A CN101169387A (zh) | 2008-04-30 |
CN101169387B true CN101169387B (zh) | 2012-02-08 |
Family
ID=39390080
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101711161A Expired - Fee Related CN101169387B (zh) | 2007-11-28 | 2007-11-28 | 一种在线确定常压塔顶石脑油质量指标的软测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101169387B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102116773B (zh) * | 2010-01-06 | 2013-06-05 | 北京化工大学 | 一种乙烯工业能效指标提取分析方法 |
CN101880543B (zh) * | 2010-07-02 | 2013-07-17 | 沈阳石蜡化工有限公司 | 一种常压塔塔盘结盐的在线清洗方法 |
CN102880905B (zh) * | 2011-07-15 | 2016-01-20 | 中国石油天然气股份有限公司 | 一种常顶油干点在线软测量方法 |
CN103472722B (zh) * | 2013-09-03 | 2016-05-11 | 邵光震 | 用于石化行业的控制方法及系统 |
CN107220252B (zh) * | 2016-03-22 | 2020-06-02 | 中国石油化工股份有限公司 | 常减压装置常压塔常一线闪点预测方法 |
CN107220705B (zh) * | 2016-03-22 | 2020-10-09 | 中国石油化工股份有限公司 | 常减压装置常压塔顶干点预测方法 |
CN107220393B (zh) * | 2016-03-22 | 2020-07-28 | 中国石油化工股份有限公司 | 常减压装置常压塔常一线干点预测方法 |
CN111523779B (zh) * | 2020-04-13 | 2024-05-31 | 湖州同润汇海科技有限公司 | 一种在线确定常减压直馏柴油馏程曲线的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1114419A (zh) * | 1995-01-19 | 1996-01-03 | 中国石油化工总公司 | 汽油与柴油质量指标的在线确定方法 |
US5572030A (en) * | 1994-04-22 | 1996-11-05 | Intevep, S.A. | Method for determining parameter of hydrocarbon |
-
2007
- 2007-11-28 CN CN2007101711161A patent/CN101169387B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5572030A (en) * | 1994-04-22 | 1996-11-05 | Intevep, S.A. | Method for determining parameter of hydrocarbon |
CN1114419A (zh) * | 1995-01-19 | 1996-01-03 | 中国石油化工总公司 | 汽油与柴油质量指标的在线确定方法 |
Non-Patent Citations (2)
Title |
---|
张笑天等.基于多神经网络模型的石脑油干点软测量.《控制工程》.2004,第11卷52-54,64. * |
颜学峰等.基于自适应偏最小二乘回归的初顶石脑油干点软测量.《化工学报》.2005,第56卷(第8期),1511-1515. * |
Also Published As
Publication number | Publication date |
---|---|
CN101169387A (zh) | 2008-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101169387B (zh) | 一种在线确定常压塔顶石脑油质量指标的软测量方法 | |
CN104965967B (zh) | 一种常减压蒸馏装置的收率实时预测方法 | |
CN101201331B (zh) | 一种在线确定初馏塔顶石脑油质量指标的软测量方法 | |
CN106845796B (zh) | 一种加氢裂化流程产品质量在线预测方法 | |
CN104765346B (zh) | 一种炼油过程全流程建模方法 | |
Mahalec et al. | Inferential monitoring and optimization of crude separation units via hybrid models | |
CA2397555A1 (en) | Petroleum distillation method and system | |
CN109324014B (zh) | 一种自适应的原油性质近红外快速预测方法 | |
CN102213949B (zh) | 一种乙烯装置价值优化方法 | |
CN107976419B (zh) | 一种由油品近红外光谱预测其性质的方法 | |
CN111209967A (zh) | 一种基于支持向量机的轧制过程板凸度预测方法 | |
CN108760789A (zh) | 一种原油快评方法 | |
CN106709654A (zh) | 一种加氢裂化流程全局运行状况评估与质量追溯方法 | |
CN116644655A (zh) | 一种基于加权目标特征回归神经网络工业过程软测量方法 | |
Yuan et al. | A comparative study of adaptive soft sensors for quality prediction in an industrial refining hydrocracking process | |
CN102998276B (zh) | 由红外光谱测定原油实沸点蒸馏曲线的方法 | |
CN109978059A (zh) | 一种原油蒸馏过程初馏塔冲塔故障预警方法 | |
CN109817287A (zh) | 一种延迟焦化模型集成方法 | |
CN111523779B (zh) | 一种在线确定常减压直馏柴油馏程曲线的方法 | |
EP3861321B1 (en) | Method of determining octane number of naphtha and of determining cetane number of diesel fuel or jet fuel using infrared spectroscopy | |
CN110763649B (zh) | 一种由近红外光谱及性质选择目标原油调合配方的方法 | |
CN111914214B (zh) | 一种基于缩减kpls模型的pta生产过程软测量方法 | |
CN113793652A (zh) | 一种基于分段智能优选的谱图化学计量学解析方法 | |
CN106709169A (zh) | 一种原油加工过程的性质估计方法 | |
CN111598305A (zh) | 一种轻烃分离装置运行状态优化及预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120208 Termination date: 20191128 |