本发明为一种由近红外光谱预测混合原油组成的方法，具体地说，是一种由混合原油的近红外光谱预测其所含各单种原油含量的方法。
背景技术：
随着市场对油品需求的不断增长和国内原油资源的不足，炼油厂正在加工越来越多的进口原油和海洋原油，且加工的品种不断增多，由于油种多，进厂批量少，因此，原油掺炼的情况普遍存在。原油的合理掺炼，不仅可以提高拔出率，而且，当一种原油的性质不适应设备的时候，可以通过原油的掺炼来改善，以提高处理量和改善产品质量等。因此，在进蒸馏装置前需要及时准确地测定混兑原油的比例。另外，原油管道输送过程中，也会出现几种原油混合的情况，因此，在输送末站也需要及时测定原油混合比例，以便得到所接收原油的性质状况。
近红外光(nir)的波长范围为780～2526nm，是由于分子振动的非谐振性使分子振动从基态向高能级跃迁时产生的，反映的是含氢基团x－h(x＝c、n、o)振动的倍频和合频吸收，因此，非常适合用于油品的物化参数测量。早在上世纪90年代前期，近红外光谱技术就被用于原油的快速分析和在线分析，可以在1min内测定原油和混兑原油的关键物化性质数据，如实沸点蒸馏曲线(tbp)、残碳、硫含量和api度等，在原油输送、原油贸易、原油调合等领域得到了应用，取得了显著的经济和社会效益。由于不同原油的近红外光谱存在一定差异，因此，有可能通过光谱方法快速测定原油的混兑比例。
cn200910157459.1“由近红外光谱预测混合原油中单种原油含量的方法”提出了一种测定混兑原油比例的方法，该方法只需知道混兑前各单种原油的近红外光谱，通过采用混兑模拟光谱建模的方式，便可准确快速测定混兑原油的比例。这种方法能够准确快速测定混兑原油的比例，从而可快速获取混兑原油的评价数据，用于原油输送和原油调合等领域。
但在实际应用中，由于该方法需要基于混兑模拟光谱利用化学计量学方法建立多元校正模型，建模过程较为繁琐，需优化选择较多建模参数，例如偏最小二乘主因子数等，因此，在很大程度上限制了这种方法的应用。
技术实现要素：
本发明的目的是提供一种由近红外光谱预测混合原油组成的方法，该法计算过程简便，预测速度快，实用性强。
本发明提供的由近红外光谱原油预测混合原油组成的方法，包括如下步骤：
(1)分别测定单种原油和混合原油的近红外光谱，并进行二阶微分处理，计算每个单种原油和混合原油在w1＝4000～5000cm-1和w2＝5000～6500cm-1两个谱区吸光度的极坐标投影得分，
(2)将所有单种原油的极坐标投影得分组成矩阵x，将混合原油的极坐标投影得分按与矩阵x同样的排列顺序组成向量y，将其代入式③，
y＝bx③
用非负约束最小二乘法求出向量b，记为
本发明方法将原油在特征谱区的吸光度分解为极座标得分，由待测混合原油的近红外光谱的极座标得分和参与混合的单种原油的近红外光谱的极座标得分预测混合原油的组成，方法简单，计算量少，预测快速，实用性强。
具体实施方式
本发明将混兑前各单种原油的近红外光谱和混合原油的近红外光谱进行极坐标投影变换，得到混兑前各单种原油的极坐标投影得分和混合原油的极坐标投影得分，再通过非负最小二乘法计算出混合原油中单种原油的含量。所用方法无需建立复杂的预测模型，预测速度快。
本发明方法中，分别按式①、②计算原油样品特征谱区的吸光度在x轴、y轴的极坐标投影得分，
式①、②中，pcxwp表示原油样品特征光谱区wp的吸光度在x轴的极坐标投影得分，pcywp表示原油样品特征光谱区wp的吸光度在y轴的极坐标投影得分，n为原油样品在特征光谱区wp采集吸光度的波长点数，ii+1为原油样品在特征光谱区wp第i+1个波长点处的吸光度，i＝0,1,2,…,n-1，p为特征谱区的序号，p为1或2。
对于每个原油样品，分别按式①、②逐一计算两个特征谱区吸光度在x、y轴的极坐标投影得分，计算时p取1或2，分别代表w1、w2特征谱区，当p＝1时，计算4000～5000cm-1谱区的吸光度在x、y轴的极坐标投影得分，p＝2时，计算5000～6500cm-1谱区的吸光度在x、y轴的极坐标投影得分。每个原油样品有四个极坐标投影得分，按一定顺序排列后，组成原油样品的得分向量，将每个单种原油样品的得分向量组成矩阵x，混合原油的得分向量为y。
测定混合原油样品和单种原油样品的近红外光谱的方法应保持一致。构成矩阵x和向量y的极坐标投影得分的排列可以是任意的，但矩阵x和向量y中的极坐标投影得分排列顺序相同。
本发明方法计算原油样品吸光度的极坐标投影得分所述的波长点数(n)为在该特征谱区按波长点采集吸光度的数目。
在4000～5000cm-1谱区的波长点数n优选220～280，在5000～6500cm-1谱区的波长点数n优选350～400。
本发明方法中，每个原油样品可得到四个极坐标投影得分：
将所有参与混合的单种原油的极坐标投影得分组成矩阵x，将混合原油的极坐标投影得分组成向量y，将其代入式③的方程，解出其中的系数向量b，向量b即为混合原油中各单种原油的含量。
对于两种单种原油(a、b)参与混兑的情况，矩阵x可为：
矩阵x的维数为2×4，
混合原油(f)的向量y为：
向量y的维数为1×4。
将矩阵x和向量y代入式③，求出向量b，为待求的混合原油中各单种原油含量，维数为1×2。
对于三种单种原油(a、b、c)参与混兑的情况，矩阵x可为：
矩阵x的维数为3×4，
将矩阵x和向量y代入式③，求出向量b，为待求的混合原油中各单种原油含量，维数为1×3。
式③中，向量b采用非负约束最小二乘法求出，需满足式④的目标函数，
式④中，m为参与混合的单种原油数目，
本发明方法适用于掺混原油，如原油掺炼、混输原油中各单种原油含量的快速确定。可通过单种原油性质及含量预测混合原油的性质。
下面通过实例进一步详细说明本发明，但本发明并不限于此。
实例中，原油的近红外光谱采用thermoantarisii傅立叶变换近红外光谱仪测定，光谱范围3800～10000cm-1，分辨率8cm-1，累积扫描次数64次，透射方式测量。
实例1
石蜡基原油a和中间基原油b两种原油混兑比例的测定。
将石蜡基原油a和中间基原油b按35％和65％(质量比)进行混兑，得到混合原油f。
(1)测定单种原油和混合原油的近红外光谱
采集原油a、原油b和混合原油f的近红外光谱，分别对原油a、原油b和混合原油f的近红外光谱进行二阶微分处理，取4000～5000cm-1和5000～6500cm-1两个区间的吸光度。
(2)计算光谱的极坐标投影得分
按式①、②，分别计算原油a、原油b和混合原油f在4000～5000cm-1和5000～6500cm-1两个谱区间的极坐标投影得分。
单种原油a得到的四个极坐标投影得分
单种原油b得到的四个极坐标投影得分
混合原油f得到的四个极坐标投影得分
(3)计算混合原油中单种原油的含量
将单种原油a和单种原油b近红外光谱在上述两个谱区吸光度的极坐标投影得分组成矩阵x，混合原油f在上述两个谱区吸光度的极坐标投影得分组成向量y，如下所列，
y＝[-0.01270.0430-0.07410.2505]
将矩阵x和向量y代入式③，采用非负约束最小二乘法计算出的向量
实例2
中间基原油a和中间基b原油两种原油混兑比例的测定。
将中间基原油b和中间基g按35％和65％(质量比)进行混兑，得到混合原油f。
按实例1(1)、(2)步的方法测定单种原油和混合原油的近红外光谱，进行二阶微分处理后，取4000～5000cm-1和5000～6500cm-1两个谱区间的吸光度，计算其极坐标投影得分。
单种原油a得到的四个极坐标投影得分
单种原油b得到的四个极坐标投影得分
混合原油f得到的四个极坐标投影得分
(3)计算混合原油中单种原油的含量
将单种原油a和单种原油b近红外光谱在上述两个谱区吸光度的极坐标投影得分组成矩阵x，混合原油f在上述两个谱区吸光度的极坐标投影得分组成向量y，如下所列，
y＝[-0.01150.0394-0.06190.2326]
将矩阵x和向量y代入式③，采用非负约束最小二乘法计算出的向量
实例3
中间基原油a、中间基原油b和石蜡基原油c三种原油混兑比例的测定。
将中间基原油a、中间基原油b和石蜡基原油c按35％、40％和25％(质量比)进行混兑，得到混合原油f。
按实例1(1)、(2)步的方法测定单种原油和混合的近红外光谱，进行二阶微分处理后，取4000～5000cm-1和5000～6500cm-1两个谱区间的吸光度，计算其极坐标投影得分。
单种原油a得到的四个极坐标投影得分
单种原油b得到的四个极坐标投影得分
单种原油c得到的四个极坐标投影得分
混合原油f得到的四个极坐标投影得分
(3)计算混合原油中单种原油的含量
将单种原油a、单种原油b和单种原油c的近红外光谱在上述两个谱区吸光度的极坐标投影得分组成矩阵x，混合原油f在上述两个谱区吸光度的极坐标投影得分组成向量y，如下所列，
y＝[-0.01070.0361-0.05560.1878]
将矩阵x和向量y代入方程③，采用非负约束最小二乘法计算出的向量
由上述结果可知，本方法计算得到的混合原油中单种原油含量与实际的混兑含量较为相符，且可预测不同类型的原油混兑后单种原油的含量，本发明方法可以用于实际预测混合原油中单种原油的含量。