CN1125329C - 一种测定渣油组分含量的方法 - Google Patents
一种测定渣油组分含量的方法 Download PDFInfo
- Publication number
- CN1125329C CN1125329C CN 99109677 CN99109677A CN1125329C CN 1125329 C CN1125329 C CN 1125329C CN 99109677 CN99109677 CN 99109677 CN 99109677 A CN99109677 A CN 99109677A CN 1125329 C CN1125329 C CN 1125329C
- Authority
- CN
- China
- Prior art keywords
- residual oil
- visible absorption
- absorption spectra
- principal component
- accordance
- 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 - Lifetime
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种渣油各组分含量的测定方法,是将具有代表性的渣油样品经浓度归一化后的紫外-可见吸收光谱组成光谱矩阵,以主成分分析的得分为特征变量进行模糊K-均值聚类;将聚类后的每种类型渣油组成的校正集的紫外-可见吸收光谱与相应洗脱色谱法测定的基础数据进行回归分析,建立校正模型;对未知渣油样品的紫外-可见吸收光谱进行主成分分析,识别渣油类型并根据相应类型渣油的校正模型计算出渣油各组分含量。本方法具有工业适用性强,准确、快速的特点,特别适合于成批相同类型渣油的常规分析。
Description
技术领域
本发明是关于渣油组分的测定方法,具体地说是采用紫外-可见吸收光谱结合化学计量学测定渣油组分含量的方法。
背景技术
渣油的充分利用是提高石油炼厂经济效益的重要手段之一。渣油的组成和性质对选择合适的渣油加工条件,充分开发和利用重油资源是十分重要的。但因为渣油中烃类数目十分繁多,结构又十分复杂,详细烃族的分离和鉴定都存在困难。渣油一般由饱和烃(Saturates)、芳烃(Aromatics)、胶质(Resins)和沥青质(Asphaltenes)四组分组成,所以一般所说的渣油组分的测定,即是对渣油中这四个组分含量的测定,也称SARA四组分测定法。
常规的渣油四组分测定方法是洗脱色谱法(EC),它的实验设备都是实验室中容量分析使用的常用设备,在简易条件下就可以开展工作,是在科研和生产中被广泛采用的最常规的方法,但是其测定过程中要用到大量对人体有害的溶剂,分析时间比较长,单个试样平行测定两次约需12个小时。
化学计量学是以数学、统计学和计算机学为手段,设计或选择最优化学量测方法,并通过解析化学量测数据,获取有关物质系统的化学及其它相关信息;使用光谱仪器分析有操作简单,速度快的特点。
Fuel,1990,Vol 69,p1381~1385中报道了一种将紫外吸收光谱与化学计量学相结合测定渣油族组分的方法,是将渣油在250~450nm的紫外光谱分段进行数学处理,再将处理后得到的6个参数,通过数学变换得到饱和烃含量,但是测得的饱和烃含量与洗脱色谱法测得的结果有9个百分点的绝对偏差。
Fuel,1993,Vol 72,p505~509中则进一步报道了在上述方法基础上改进了的饱和烃含量的测定方法,是将上述方法中得到的6个参数中的4个用偏最小二乘法校正得到一个参数,再通过数学变换得到饱和烃含量,但是其测定结果与洗脱色谱法方法测得的结果仍有6.9个百分点的绝对偏差。
Fuel,1992,Vol 71,p559~563中报道了建立渣油中沥青质含量校正模型的方法,但由于使用的是部分紫外吸收光谱信息,使沥青质测定结果与洗脱色谱法的相对偏差高达23%。
上述关于渣油组分测定方法的文献均使用的是部分紫外吸收光谱的信息,不但测定的饱和烃或沥青质的含量偏差较大,且未能对渣油中芳烃和胶质的含量进行测定。
渣油中芳烃、胶质和沥青质由于其分子结构不同,致使它们在190~700nm的紫外-可见波长区域内有不同的吸收光谱;与前述三组分不同,饱和烃在此波长区域内没有吸收,其含量对渣油的吸收有负贡献,即饱和烃在渣油中含量越小,吸收越弱。因此,190~700nm的波长区域内的紫外-可见吸收光谱实际上包含了渣油四组分的信息。
在一篇题为“高效液相色谱及紫外光谱法快速分析重油馏份的族组成”的博士论文(石油化工科学研究院,1996)中报道了渣油的紫外吸收光谱与高效液相色谱通过偏最小二乘法(PLS)相关联测定渣油四组分的方法,校正模型是使用包括减压渣油和二次加工渣油在内的渣油样品的混合模型,基础数据是由高效液相色谱(HPLC)方法测定的。这种方法存在如下缺点:(1)基础数据由高效液相色谱测得,与工厂实际操作中采用的洗脱色谱法对渣油组分的切割点不同,且高效液相色谱仪的操作比较复杂,费用大,所以对工厂实际操作的适用性差。(2)由于各种渣油的组成差异较大,异常样品对建立的校正模型的准确性有很大影响,从而影响测定结果的准确性。
发明内容
本发明的目的是在现有技术的基础上,提供一种较现有技术更适合于工业应用,快速且更准确的测定渣油四组分的方法。
本发明的发明点在于:将因子分析和模糊聚类识别引入到紫外-可见吸收光谱结合化学计量学测定渣油组分的过程中。这是因为渣油组分的繁杂性,用一个渣油组分校正模型难于得到准确的结果。
因子分析中较完善的一种分析方法是主成分分析,它通过将特征压缩,把原有的各个特征利用线性变换得到一批新的特征,每个特征都是原有特征的函数,但新特征总数少于原有特征数,这样新特征既保留了原有特征的主要信息,滤除噪音,又减少了特征个数。将主成分分析得到的新特征变量引入模糊K-均值聚类算法中,可使渣油样品聚集分类,为针对性建立渣油组分含量的校正模型建立基础。
渣油组分的基础数据采用工厂中普遍使用的洗脱色谱法获得,是本发明的另一个发明点。
本发明提供的方法包括以下步骤:
1、将具有代表性的渣油样品经浓度归一化后的紫外-可见吸收光谱组成光谱矩阵,进行主成分分析,并以主成分得分为特征变量进行模糊K-均值聚类。
2、将模糊聚类后区分开的各类型渣油组成的校正集的紫外-可见吸收光谱与相应洗脱色谱法测定的基础数据用数学方法进行回归分析,建立各类型渣油的校正模型。
3、对未知渣油样品的紫外-可见吸收光谱进行主成分分析,根据主成分得分识别渣油类型并由步骤2中的相应类型渣油的校正模型计算各组分含量。
本发明所说的渣油紫外-可见吸收光谱图的扫描范围是190~700nm。
所说的具有代表性的渣油样品可包括常压渣油、减压渣油、二次加工油,二次加工油可以包括加氢渣油、回炼油、油浆和焦化蜡油。
本发明所说的对渣油的紫外-可见吸收光谱的主成分分析,是将光谱的量测矩阵Y进行分解,分解为三个矩阵的乘积,即Y=USVt,其中,S为对角矩阵,它收集了Y矩阵的特征值,U和Vt分别收集了特征值对应的列特征矢量和行特征矢量,即主成分载荷矩阵和主成分得分矩阵。
所说的模糊K-均值聚类算法,目标函数Jm(U,V)为 式中
为样本Xi={xki}与聚类中心Vj={vkj}之间的欧氏距离,其中:i=1,2,...n,n为样品数;j=1,2,...c,c为分类数;k=1,2,...s,s为特征数;uij为样本xi对第j类的隶属度;m为加权指数。
目标函数Jm(U,V)表示样本xi与各个聚类中心Vj的带权距离平方和,其权重为样本xi隶属于类Cj的隶属度uij的m次方,而最佳聚类是使目标函数Jm(U,V)最小。因此,要得到最佳的聚类结果,就要求得适当的隶属度uij和聚类中心Vj。,当m>1,Xi≠Vj时,可用式(2)、(3)迭代计算出隶属度uij和聚类中心Vj。 具体算法如下:(1).固定分类数c,加权指数m;收敛门限ε;选取初始隶属度矩阵U(0),其元素uij满足:
0≤uij≤1,i,j
(2).根据式(3)和U(q)求聚类中心Vj (q),q为迭代次数。
(3).根据式(2)和求得的Vj (q),求U(q+1)。
(4).若max{|U(q)-U(q+1)|}≤ε,则停止迭代,U(q+1)及相应的Vj (q)为所求.否则返回步骤(2),继续迭代。
(5).在得到的隶属度矩阵U中,令每列中最大元素为1,其余为0,得到一个普通分类矩阵,即为分类结果。
本发明所说的获取渣油四组分基础数据的洗脱色谱法,已披露于《石油化工分析方法》(科学出版社,1990)第31页上记载的RIPP10-90方法中。
本发明所说的洗脱色谱法测得的四组分基础数据,与相应的紫外-可见吸收光谱数据进行回归分析的数学方法是指多元校正方法,可以是经典最小二乘法(CLS)、逆最小二乘法(ILS)、多元线性回归(MLR)、主成分回归(PCR)、偏最小二乘法(PLS)、稳健偏最小二乘法(RPLS)或人工神经网络(ANN),其中优选偏最小二乘法(PLS)。
采用多元校正方法将紫外-可见吸收光谱与相应基础数据进行回归分析,即可获得相应组分的校正模型。
当采用偏最小二乘法(PLS)时,其根据是比耳-朗白定律:Y=XK+E,其中:
Y:m个样品,n个波长点数组成的(m×n)吸光度矩阵;
X:m个样品,1个组分含量组成的(m×1)浓度矢量;
K:1个组分,n个波长点数组成的(1×n)灵敏度矢量;
E:m×n吸光度残差矩阵。
其校正模型的建立过程是:
(1).吸光度矩阵Y和浓度矢量X的各元素进行均值中心化处理,即所有数据减去其对应数据组(列向量)的平均值。
(2).对归一后的吸光度矩阵Y和浓度矢量X进行主成分分解,
rk=bktk (3)其中:tk(m×1)为吸光度矩阵的因子得分,
vk(1×n)为吸光度矩阵的因子载荷,
rk(m×1)为浓度矢量的因子得分,
qk(1×1)是个数,浓度矢量的因子载荷,
bk(1×1)为rk和tk的回归系数,
p为主因子数。
(3).求解特征矢量及主因子数p:
下面是对特征矢量及主因子数p求解的一般过程:
当忽略残差阵E时,根据式(1)及(2),p=1时有:
Y=tvt 左乘tt得:vt=ttY/ttt
右乘v得:t=Yv/vtv
X=rq 左乘rt得:q=rtX/rtr,两边同除以q得:r=X/q
[1]求吸光度矩阵的权重向量w,
取浓度阵X的某一列作r的起始迭代值,以r代替t,计算w
方程:Y=rwt 解:t=Yw/wtw
[2]归一化w:wt 归一=wt/‖wt‖
[3]求吸光度矩阵的因子得分t,由归一化后w计算t,
方程:Y=twt 解:t=Yw/wtw
[4]求浓度矢量的权重u值,以t代替r计算u,
方程:X=tu 解:u=ttX/ttt
[5]求浓度矢量的因子得分r,由u计算r,
方程:X=ru 解:r=X/u
再以此r代替t返回第[1]步计算wt,由wt计算t新,如此反复迭代,如t已收敛(‖t新-t旧‖≤10-6‖t新‖),继续下步运算,否则返回步骤[1]。
[6]由收敛后的t求吸光度矩阵的载荷向量v,
方程:Y=tvt 解:vt=ttY/ttt
[7]由r求浓度矢量的载荷q值,
方程:X=rq 解:q=rtX/rtr
由此求出对应于第一个主因子的r1,q1,t1,v1 t
代入公式(3)求出b1:b1=r1t1/t1 tt1
[8]由公式(1)和(2)计算残差E,
EX,1=X-r1q1=X-b1t1q1 (4)
EY,1=Y-t1v (5)
[9]以Ex,1代替X,EY,1代替Y,返回步骤[1]计算下一个组分的
r2,q2,t2,v2 t,b2
[10]按式(4)及(5)计算EX,2、EY,2,以此类推,求出X、Y的诸
主因子。
用交互检验法确定主因子数p。
经过以上过程,校正模型建立完成。
由未知试样的紫外-可见吸收光谱y(未知)及校正模型建立过程中存储的vk计算tk(未知):tk(未知)=y(未知)vk/vk tvk;由求出的tk(未知)及校正过程中存储的bk计算rk(未知):rk(未知)=bktk(未知);再由求出的rk(未知)及校正模型建立过程中存储的qk由公式(2)
计算组分含量x(未知)。
本发明采用的紫外-可见吸收光谱与化学计量学相结合的渣油组分含量的测定方法,具有如下优点:
1、因为校正集的基础数据由工业上普遍采用的洗脱色谱法(EC)获得,所以此方法的工业适用性强。
2、由于首先对渣油样品进行了类型识别,然后再选用对应的校正模型进行组分含量测定,提高了测定结果的准确性。通过对30个未知渣油样品的含量测定结果可以看出,本发明方法与洗脱色谱法的测定结果相比,绝对偏差均小于2个百分点,满足洗脱色谱法对渣油组分含量的测量要求。
3、本方法对单个渣油样品平行测定两次的时间为45分钟,明显少于洗脱色谱法的12小时。特别是它更适合于成批相同类型渣油的常规分析,当对成批相同类型渣油进行组分测定时,为每小时15个样品。
具体实施方式
以下的实施例将对本发明作进一步的说明,但本发明的保护范围并不受实施例的限制。
在实施例中,使用PE Lambda 16紫外光谱仪和奔腾586微机;试剂为二氯甲烷(分析纯)和正庚烷(分析纯)。
实例
这些实例说明采用本发明方法对30个未知渣油样品的四组分含量的测定过程,其结果表明本发明方法与常规的洗脱色谱法相比,各组分含量的绝对偏差均小于2个百分点,可满足对渣油组分含量的测量要求。
1、渣油紫外-可见吸收光谱的测定:
称取渣油样品16.0-17.0mg,先用0.5mL二氯甲烷溶解,然后用正庚烷稀释到250mL,得到浓度约为0.06mg/mL的二氯甲烷与正庚烷混合溶剂的样品溶液。
以2‰(体积比)的二氯甲烷与正庚烷的混合溶剂为参比样品,在波长范围190~700nm进行紫外吸收光谱测定,比色皿厚度1.0cm。渣油原始紫外-可见吸收光谱除以配样的浓度,得到归一化的渣油紫外-可见吸收光谱。
2、选取75个渣油样品组成聚类集,由步骤1获得每个样品的归一化紫外-可见吸收光谱并对其进行主成分分析和模糊K-均值聚类。
此聚类集渣油的前三个主成分得分累积贡献率达到99.9%,即只用前三个主成分已几乎包含了原紫外-可见吸收光谱的全部信息。其中第一主成分得分约占90%,信息强,对渣油紫外吸收光谱的贡献最大,主要反映的是渣油中芳烃的信息。若以第一主成分得分作为聚类特征变量,实际上是以渣油芳烃含量作为聚类特征变量的主体信息,会掩盖其它组分的信息,不能得到满意的聚类结果。而选取第二主成分和第三主成分得分为特征变量进行模糊K-均值聚类,可得到满意的聚类结果,将聚类集中渣油样品可聚类为三种类型,编号为A、B、C。聚类集中每种渣油样品的类型见表1。
3、将步骤2中25个A类型渣油样品组成校正集,采用《(石油化工分析方法》(科学出版社,1990)第31页上记载的RIPP10-90方法测定四组分含量基础数据,测定结果列于表2中。采用偏最小二乘法将紫外-可见吸收光谱与相应基础数据进行回归分析,获得A类渣油四组分含量的校正模型。
B类和C类渣油的校正模型的建立同A类。
4、对30个未知渣油样品测定紫外吸收光谱,进行主成分分析后选取第二主成分和第三主成分得分识别渣油类型,再由未知渣油的紫外-可见吸收光谱,调用步骤3中相应渣油校正模型计算各组分含量。测定结果见表2,其中:EC为采用洗脱色谱法测定的组分含量;UV为采用本发明方法测定的组分含量;DE为本发明方法与洗脱色谱法测定的组分含量的差值,即EC-UV。
表1聚类集样品序号、隶属度及聚类结果
序号NO. uiA uiB uiC 聚类结果
1 0.9958 0.0039 0.0004 A
2 0.9982 0.0005 0.0013 A
3 0.9930 0.0050 0.0020 A
4 0.9990 0.0008 0.0001 A
5 0.9688 0.0095 0.0217 A
6 0.9998 0.0002 0.0000 A
7 0.9988 0.0007 0.0005 A
8 0.9943 0.0013 0.0044 A
9 0.9983 0.0014 0.0003 A
10 0.9997 0.0002 0.0001 A
11 1.0000 0.0000 0.0000 A
12 0.9993 0.0005 0.0002 A
13 0.9999 0.0001 0.0000 A
14 0.9964 0.0025 0.0011 A
15 1.0000 0.0000 0.0000 A
16 1.0000 0.0000 0.0000 A
17 0.9986 0.0006 0.0008 A
18 0.9992 0.0007 0.0001 A
19 0.9944 0.0027 0.0029 A
20 0.9990 0.0007 0.0003 A
21 1.0000 0.0000 0.0000 A
22 0.0000 1.0000 0.0000 B
23 0.0320 0.9624 0.0056 B
24 0.0003 0.9996 0.0001 B
25 0.0044 0.9949 0.0007 B
26 0.0179 0.9808 0.0013 B
27 0.0012 0.9987 0.0001 B
28 0.0029 0.9969 0.0002 B
29 0.0011 0.9987 0.0002 B
30 0.0002 0.9998 0.0000 B
31 0.0682 0.9302 0.0016 B
32 0.0099 0.9886 0.0015 B
33 0.0023 0.9971 0.0006 B
34 0.0018 0.9978 0.0004 B
35 0.0000 1.0000 0.0000 B
36 0.0344 0.9645 0.0011 B
37 0.0060 0.9929 0.0011 B
38 0.0112 0.9879 0.0009 B
39 0.0184 0.9811 0.0005 B
40 0.0016 0.9983 0.0001 B
41 0.0191 0.9800 0.0009 B
42 0.2212 0.7764 0.0024 B
43 0.0020 0.9976 0.0004 B
44 0.1966 0.7438 0.0596 B
45 0.0004 0.9996 0.0000 B
46 0.0016 0.9983 0.0001 B
47 0.0005 0.9994 0.0001 B
48 0.0372 0.9621 0.0007 B
49 0.0000 0.0000 1.0000 C
50 0.0000 0.0000 1.0000 C
51 0.0006 0.0002 0.9992 C
52 0.0001 0.0000 0.9999 C
53 0.0000 0.0000 1.0000 C
54 0.0127 0.0007 0.9866 C
55 0.0003 0.0000 0.9997 C
56 0.0057 0.0005 0.9938 C
57 0.0009 0.0001 0.9990 C
58 0.0000 0.0000 1.0000 C
59 0.0099 0.0011 0.9890 C
60 0.0015 0.0002 0.9983 C
61 0.0007 0.0001 0.9992 C
62 0.0007 0.0002 0.9991 C
63 0.0000 0.0000 1.0000 C
64 0.0002 0.0000 0.9998 C
65 0.0001 0.0000 0.9999 C
66 0.0000 0.0000 1.0000 C
67 0.0024 0.0003 0.9973 C
68 0.0002 0.0001 0.9997 C
69 0.0000 0.0000 1.0000 C
70 0.0003 0.0000 0.9997 C
71 0.0000 0.0000 1.0000 C
72 0.0000 0.0000 1.0000 C
73 0.0000 0.0000 1.0000 C
74 0.0002 0.0000 0.9998 C
75 0.0008 0.0002 0.9990 C
表2 30个未知渣油样品四组分含量的测定结果
未知样品编号 | 类型识别结果 | 饱和烃 | 芳烃 | 胶质 | 沥青质 | ||||||||
EC | UV | DE | EC | UV | DE | EC | UV | DE | EC | UV | DE | ||
1# | A | 54.33 | 54.53 | 0.10 | 22.74 | 23.30 | 0.56 | 20.93 | 20.70 | -0.23 | 2.0 | 2.30 | 0.30 |
2# | A | 39.26 | 39.00 | -0.26 | 37.36 | 36.69 | -0.67 | 19.78 | 20.44 | 0.66 | 3.6 | 3.85 | 0.25 |
3# | A | 52.10 | 52.27 | 0.17 | 21.78 | 21.11 | -0.67 | 24.72 | 24.46 | -0.26 | 1.4 | 1.79 | 0.39 |
4# | A | 37.17 | 36.89 | -0.28 | 39.09 | 39.16 | 0.07 | 18.34 | 19.74 | 1.40 | 5.4 | 4.85 | 0.55 |
5# | A | 38.93 | 37.91 | -1.02 | 38.44 | 38.79 | 0.35 | 18.33 | 17.76 | -0.57 | 4.3 | 3.85 | 0.45 |
6# | A | 49.70 | 49.32 | -0.38 | 22.19 | 23.16 | 0.97 | 26.41 | 26.00 | -0.41 | 1.7 | 2.07 | 0.37 |
7# | A | 43.01 | 43.09 | 0.08 | 34.60 | 34.56 | -0.04 | 18.79 | 18.12 | -0.67 | 3.6 | 3.16 | 0.44 |
8# | A | 39.51 | 39.16 | -0.35 | 44.21 | 44.07 | -0.14 | 13.18 | 14.00 | 0.82 | 3.1 | 3.21 | 0.11 |
9# | A | 35.38 | 36.49 | 1.11 | 45.96 | 44.97 | 0.99 | 15.56 | 14.68 | -0.88 | 3.1 | 2.87 | -0.23 |
10# | A | 34.86 | 34.17 | -0.69 | 44.69 | 45.94 | 1.25 | 17.15 | 15.97 | -1.18 | 3.3 | 3.35 | 0.05 |
11# | B | 15.46 | 13.67 | -1.79 | 32.32 | 33.42 | 1.10 | 49.02 | 48.18 | -0.84 | 3.2 | 3.54 | 0.34 |
12# | B | 32.41 | 33.17 | 0.76 | 31.05 | 31.72 | 0.67 | 30.34 | 32.32 | 1.98 | 6.2 | 7.05 | 0.85 |
13# | B | 13.59 | 14.52 | 0.93 | 28.24 | 29.99 | 1.75 | 50.07 | 49.64 | 0.42 | 8.1 | 8.41 | 0.31 |
14# | B | 15.95 | 17.65 | 1.70 | 33.01 | 31.31 | -1.70 | 44.14 | 45.70 | 1.56 | 6.9 | 7.08 | 0.18 |
15# | B | 11.79 | 10.00 | -1.78 | 44.75 | 43.42 | -1.33 | 33.86 | 33.02 | -0.84 | 9.6 | 9.35 | -0.25 |
16# | B | 13.36 | 11.48 | -1.88 | 46.82 | 46.70 | -0.12 | 28.92 | 30.34 | 1.42 | 10.9 | 10.65 | -0.25 |
17# | B | 15.32 | 15.01 | -0.31 | 37.81 | 37.27 | -0.54 | 41.17 | 39.94 | -1.23 | 5.7 | 5.48 | -0.22 |
18# | B | 23.93 | 24.39 | 0.46 | 33.31 | 34.53 | 1.22 | 39.86 | 41.21 | 1.35 | 2.9 | 3.03 | 0.14 |
19# | B | 33.23 | 33.20 | -0.03 | 38.58 | 38.65 | 0.07 | 27.09 | 27.02 | -0.07 | 1.1 | 1.72 | 0.62 |
20# | B | 8.15 | 6.77 | -1.38 | 34.14 | 35.64 | 1.50 | 48.41 | 48.41 | 0.00 | 9.3 | 9.46 | 0.16 |
21# | C | 56.38 | 57.02 | 0.64 | 25.26 | 25.91 | 0.65 | 17.36 | 18.27 | 0.91 | 1.0 | 0.89 | -0.11 |
22# | C | 57.53 | 58.61 | 1.08 | 24.02 | 24.04 | 0.02 | 17.25 | 17.14 | -0.11 | 1.2 | 1.22 | 0.02 |
23# | C | 45.31 | 45.82 | 0.51 | 39.66 | 38.88 | -0.78 | 13.43 | 14.90 | 1.47 | 1.6 | 1.06 | -0.54 |
24# | C | 56.88 | 56.75 | -0.13 | 27.41 | 29.13 | 1.72 | 12.61 | 12.97 | 0.36 | 3.10 | 2.56 | -0.54 |
25# | C | 54.24 | 54.20 | -0.04 | 33.23 | 33.86 | 0.63 | 11.63 | 11.76 | 0.13 | 0.9 | 0.45 | 0.45 |
26# | C | 51.14 | 51.13 | -0.01 | 24.88 | 23.09 | -1.79 | 21.38 | 22.28 | 0.90 | 2.6 | 2.59 | -0.02 |
27# | C | 50.28 | 50.53 | 0.25 | 24.08 | 23.04 | -1.04 | 22.84 | 24.08 | 1.24 | 2.8 | 2.79 | -0.01 |
28# | C | 51.07 | 50.97 | -0.10 | 23.47 | 23.35 | -0.12 | 22.86 | 22.51 | -0.35 | 2.6 | 2.87 | 0.27 |
29# | C | 49.24 | 48.53 | -0.71 | 24.42 | 24.83 | 0.41 | 24.84 | 23.20 | -1.64 | 1.5 | 2.06 | 0.56 |
30# | C | 50.58 | 50.34 | -0.24 | 35.71 | 35.60 | -0.11 | 12.41 | 12.78 | 0.37 | 1.3 | 0.71 | -0.59 |
Claims (6)
1、一种紫外-可见吸收光谱结合化学计量学测定渣油各组分含量的方法,其特征是:
(1)将具有代表性的渣油样品经浓度归一化后的紫外-可见吸收光谱组成光谱矩阵,进行主成分分析,并以主成分得分为特征变量进行模糊K-均值聚类;
(2)将模糊聚类后区分开的各类型渣油组成的校正集的紫外-可见吸收光谱与相应洗脱色谱法测得的基础数据用数学方法进行回归分析,建立各类型渣油的校正模型;
(3)对未知渣油样品的紫外-可见吸收光谱进行主成分分析,根据主成分得分识别渣油类型并由步骤(2)中建立的相应类型渣油的校正模型计算组分含量。
2、按照权利要求1所述的方法,其中所说的紫外-可见吸收光谱的波长范围是190~700nm。
3、按照权利要求1所述的方法,其中所说的具有代表性的渣油样品是常压渣油、减压渣油和二次加工油。
4、按照权利要求3所述的方法,其中所说的二次加工油是加氢渣油、回炼油、油浆和焦化蜡油。
5、按照权利要求1所述的方法,其中步骤(2)所说的数学方法是多元校正方法。
6、按照权利要求5所述的方法,其中所说的多元校正方法是经典最小二乘法(CLS)、逆最小二乘法(ILS)、多元线性回归(MLR)、主成分回归(PCR)、偏最小二乘法(PLS)、稳健偏最小二乘法(RPLS)或人工神经网络(ANN)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 99109677 CN1125329C (zh) | 1999-07-06 | 1999-07-06 | 一种测定渣油组分含量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 99109677 CN1125329C (zh) | 1999-07-06 | 1999-07-06 | 一种测定渣油组分含量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1283790A CN1283790A (zh) | 2001-02-14 |
CN1125329C true CN1125329C (zh) | 2003-10-22 |
Family
ID=5274074
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 99109677 Expired - Lifetime CN1125329C (zh) | 1999-07-06 | 1999-07-06 | 一种测定渣油组分含量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1125329C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100460858C (zh) * | 2005-12-13 | 2009-02-11 | 上海神开石油化工装备股份有限公司 | 钻井液中含油气在线光谱测定方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100419409C (zh) * | 2003-11-26 | 2008-09-17 | 广东省农业科学院作物研究所 | 一种非破坏性测定花生单粒种子含油量的方法 |
CN102103079B (zh) * | 2010-12-31 | 2013-12-11 | 聚光科技(杭州)股份有限公司 | 一种光谱分析方法 |
CN104165821A (zh) * | 2014-08-01 | 2014-11-26 | 湖北富邦科技股份有限公司 | 一种蜡油中胶质的测定方法 |
CN107709983B (zh) * | 2015-06-10 | 2020-05-22 | 真空紫外线分析学股份有限公司 | 使用真空紫外光谱和气相色谱对复杂样品进行详细的批量分类分析的方法 |
CN106290676A (zh) * | 2016-08-25 | 2017-01-04 | 何恺源 | 一种结合分析化学技术与模型延展获取减压渣油分子组成的方法 |
CN109668856B (zh) * | 2017-10-17 | 2021-06-11 | 中国石油化工股份有限公司 | 预测lco加氢原料与产物的烃族组成的方法和装置 |
CN109668854B (zh) * | 2017-10-17 | 2021-06-11 | 中国石油化工股份有限公司 | 预测lco加氢原料与产物的烃族组成的方法和装置 |
CN108387673B (zh) * | 2018-03-07 | 2020-08-14 | 河南理工大学 | 一种混合物成分闪速定性识别方法 |
CN109115710A (zh) * | 2018-07-27 | 2019-01-01 | 重庆中烟工业有限责任公司 | 一种烟用爆珠内液质量监测方法 |
-
1999
- 1999-07-06 CN CN 99109677 patent/CN1125329C/zh not_active Expired - Lifetime
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100460858C (zh) * | 2005-12-13 | 2009-02-11 | 上海神开石油化工装备股份有限公司 | 钻井液中含油气在线光谱测定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1283790A (zh) | 2001-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1125331C (zh) | 一种测定渣油组分含量的方法 | |
CN1125329C (zh) | 一种测定渣油组分含量的方法 | |
CN1838107A (zh) | 处理系列波谱尤其是nmr波谱的方法 | |
CN1423749A (zh) | 分析选定的多组分样品的方法 | |
EP0871936B1 (en) | Method for predicting chemical or physical properties of complex mixtures | |
Lupoi et al. | Characterization of woody and herbaceous biomasses lignin composition with 1064 nm dispersive multichannel Raman spectroscopy | |
CN1690713A (zh) | 对样本进行分析以提供表征数据的方法、系统和指令程序 | |
CN1251128C (zh) | 文字列匹配装置和文字列匹配方法 | |
Titaley et al. | Automating data analysis for two-dimensional gas chromatography/time-of-flight mass spectrometry non‐targeted analysis of comparative samples | |
CN1143188A (zh) | 在光谱测定法中稳定光谱的方法 | |
CN1796979A (zh) | 近红外光谱测定汽油中二烯烃含量的方法 | |
CN117309838A (zh) | 一种基于三维荧光特征数据的工业园区水体污染溯源方法 | |
CN111650151A (zh) | 一种利用近红外光谱法评估湿地松林分松脂价值的方法 | |
CN100443883C (zh) | 用近红外光谱测定加氢尾油链烷烃族组成的方法 | |
CN1796980A (zh) | 红外光谱测定c10~c13正构烃中二烯烃含量的方法 | |
CN1125330C (zh) | 沥青中蜡含量的测定方法 | |
CN115616204A (zh) | 一种聚对苯二甲酸乙二醇酯再生料的鉴别方法及系统 | |
CN112834448B (zh) | 一种用于水污染物分析的光谱分解方法 | |
CN1270179C (zh) | 锶同位素质谱分析的方法 | |
CN1662858A (zh) | 设计试验的方法和计算机 | |
CN113138181B (zh) | 一种对清香型原酒品质分级的方法 | |
CN1010614B (zh) | 局部分子结构分析仪 | |
CN1683923A (zh) | 氟代烃浓度的测定方法、测定装置和全氟化碳的制造方法 | |
Guo et al. | Turning Metabolomics Data Processing from a “Black Box” to a “White Box” | |
Berteau | Authors’ answer to a letter to the editor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CX01 | Expiry of patent term | ||
CX01 | Expiry of patent term |
Granted publication date: 20031022 |