一种管输原油的混油界面检测方法与流程

1.本发明涉及原油管道输送领域，具体为一种管输原油的混油界面快速检测方法。


背景技术：

2.在原油管道输送领域，通常采用顺序输送的方法，即在一条管道中按照一定批量和次序，连续输送不同油种的一种输送方式。顺序输送可以提高管道利用率，运输量大，经济性好，自动控制方便。
3.在管道顺序输送原油的过程中，两种原油交替处会形成混油段，因此接收原油的收油站需要实时检测混油界面，以便及时、准确地进行原油切割收油。目前，工业中混油界面检测多采用密度测量法，该方法操作简单，对于密度相差较大的原油具有较好的分辨能力，应用广泛。然而，这种方法对于密度相近的不同种类原油却难以区分，进而降低了管道输送的灵活性。
4.近红外光谱分析技术近年来已开始运用于原油性质快速分析领域，相较于手工化验，具有测量速度快、精度高等优点。目前已有将基于近红外光谱的分析技术应用于混油界面检测的研究，但采用的方法仍以主元统计分析等机器学习方法为主，存在建模较为复杂、计算量偏大、实时性偏弱等问题。


技术实现要素：

5.本发明针对背景技术中存在的问题，利用近红外光谱分析及导数加和技术，对混油界面进行实时检测与跟踪。该方法具有以下步骤：
6.1)实时采集管道顺序输送原油的近红外光谱，选取5200cm
‑1～6200cm
‑1波数范围内的光谱作为特征区光谱，并对特征区内的光谱进行矢量归一化预处理；
7.2)以当前时刻k为起点，向前连续读取p个时刻的光谱吸光度a
ij
，并针对每一个波数点i，计算连续p次光谱的吸光度标准偏差s
ik
：
[0008][0009][0010]
式中，i＝1,2,...,m，m为特征区内光谱的波数点总数，a
ij
为时刻j采集的光谱序列中第i个波数点的光谱吸光，j＝k
‑
p+1,k
‑
p+2,...,k；p＝0.025πd2lt
‑1q
‑1，其中l为混油段长度，d为管道内直径，q为油品在管道中的体积流量，t为光谱采集周期，p向下取整。
[0011]
3)采用savitzky
‑
golay卷积求导方法，计算吸光度标准偏差s
ik
的一阶导数s
ik
′
：
[0012]
[0013]
式中，g
t
为权重因子，平滑窗口宽度为2ω+1，ω为整数，
△
k为差分宽度；平滑窗口宽度设置为5，以5点二次平滑方法求取一阶导数s
ik
′
：
[0014][0015]
4)计算光谱特征区内各波数点一阶导数加和值：
[0016]
5)依据加和值s
k
′
的大小，按照以下步骤判断管道输送原油的混油阶段：
[0017]5‑
1)初始化混油标志位tag＝0；
[0018]5‑
2)实时采集管道顺序输送原油的近红外光谱；
[0019]5‑
3)计算当前时刻特征区内光谱各波数点一阶导数加和值s
k
′
；
[0020]5‑
4)判断s
k
′
是否小于等于0.005：若是则转步骤5
‑
5)；否则判断为混油开始阶段，置tag＝1，转步骤5
‑
7)；
[0021]5‑
5)判断s
k
′
是否大于等于
‑
0.005：若是则转步骤5
‑
6)；否则判断为混油结束阶段，置tag＝0，转步骤5
‑
7)；
[0022]5‑
6)判断tag是否为0：若是则判断为单一原油阶段；否则判断为混油中间阶段；
[0023]5‑
7)输出原油混油界面的判断结果；
[0024]5‑
8)转步骤5
‑
2)。
[0025]
6)输出判断结果至原油管输调度控制系统，及时对管输原油进行切割收油。
[0026]
有益效果：
[0027]
本发明公开了一种管输原油的混油界面检测方法，该方法采用近红外光谱分析及导数加和技术，实时对混合油品进行检测与跟踪。该方法即使面对密度相近的原油，也能够实时分析原油特征光谱变化情况，可准确判断原油混油界面开始、中间和结束阶段，且建模简单，运算速度快，可及时指导收油站进行原油切割收油。这对于提高原油管输调度过程的实时性，提升管输灵活性等具有重要价值。
附图说明
[0028]
图1为本发明一种管输原油的混油界面检测方法流程图；
[0029]
图2为本发明实施例中某管线原油光谱样本的吸光度标准偏差计算结果；
[0030]
图3为本发明实施例中某管线原油吸光度标准偏差的一阶导数加和值计算结果；
[0031]
图4为本发明实施例中某管线分批分次输送不同种类原油的混油界面判断结果。
具体实施方式
[0032]
下面结合附图和具体实施例对本发明作进一步说明。由具体的操作流程说明本方法在管道输送原油在线检测混油界面过程中的实施效果。本实施案例在以本发明技术方案为前提下进行实施，但本发明的保护范围不限于下述的实施例。
[0033]
某长输管线顺序输送卡斯原油及荣卡多原油，管线上安装ftpa2000
‑
260在线光谱仪，采用2mm光程流通池，40℃恒温后送至近红外光谱仪扫描。为准确判断原油混油界面的开始、中间和结束阶段，决定采用本发明技术进行分析检测。具体过程如下：
[0034]
1)选取波数段为5200cm
‑1～6200cm
‑1的近红外光谱并对其进行矢量归一化预处
理。其中，表1为预处理后管输原油近红外光谱的部分数据。
[0035]
表1矢量归一化处理后的管输原油近红外光谱
[0036][0037][0038]
2)在每一个波数点下，计算连续p次的光谱样本吸光度标准偏差。
[0039]
经工艺理论即前期实测历史数据推算，管输原油的混油段长度约为3.0km，流速约为3.0ms，光谱采集周期设定为10s，则p＝0.025πd2lt
‑1q
‑1＝10.2，取p＝10。
[0040]
如表1所示，针对波数点5207.42cm
‑1，以光谱标号01为起点，计算01～10共计p＝10个的吸光度a
1j
标准偏差s
11
：
[0041][0042][0043]
同理，还是针对波数点5207.42cm
‑1，以02位起点，光谱标号02～11的吸光度标准偏
差s
12
的值为0.0001；类似的，其它波数点及不同时间序列计算得到的吸光度标准偏差部分数据如表2所示：
[0044]
表2连续序列光谱的吸光度标准偏差
[0045][0046][0047]
对其样本吸光度标准偏差数据进行可视化处理，其结果如图2所示。
[0048]
3)采用savitzky
‑
golay卷积求导方法，计算吸光度标准偏差的一阶导数。
[0049]
如表2所示，以波数点5207.42cm
‑1为例，时间序列16～20的吸光度标准偏差s
1k
分别为[0.0146,0.0135,0.0113,0.0085,0.0061]，以5点二次平滑方法，计算吸光度标准偏差s
1(20)
的一阶导数s
1(20)
′
：
[0050][0051]
同理，其它计算得到的一阶导数值的部分数据如表3所示：
[0052]
表3卷积平滑处理的吸光度标准偏差一阶导数
[0053][0054][0055]
4)计算光谱特征区内各波数点一阶导数加和值。
[0056]
以表3中时间序列15为例，其在光谱特征区内各波数点一阶导数s
i(15)
′
为[0.0001,0.0001,...,0.0007,0.0007]，则时间序列15在光谱特征区内各波数点一阶导数的加和值s
15
′
按下式计算：
[0057][0058]
光谱特征区内各波数点一阶导数加和值的计算结果如图3所示。
[0059]
5)依据加和值s
k
′
的大小，按照如下顺序判断管道输送原油的混油阶段：
[0060]5‑
1)初始化混油标志位tag＝0；
[0061]5‑
2)实时采集管道顺序输送原油的近红外光谱；
[0062]5‑
3)计算当前时刻特征区光谱内各波数点一阶导数加和值s
k
′
；
[0063]5‑
4)判断s
k
′
是否小于等于0.005：若是则转步骤5
‑
5)；否则判断为混油开始阶段，置tag＝1，转步骤5
‑
7)；
[0064]5‑
5)判断s
k
′
是否大于等于
‑
0.005：若是则转步骤5
‑
6)；否则判断为混油结束阶段，置tag＝0，转步骤5
‑
7)；
[0065]5‑
6)判断tag是否为0：若是则判断为单一原油阶段；否则判断为混油中间阶段；
[0066]5‑
7)输出原油混油界面的判断结果；
[0067]5‑
8)转步骤5
‑
2)。
[0068]
导数值加和值从0至最大值，再由最大值至最小值，又由最小值至0的过程，记作1个混油段的混合油周期。图3清晰展现了1个混合油周期的变化过程，可见，该方法可以准确判断混油界面的开始、中间和结束阶段，且由于该方法无需复杂的建模过程，运算速度快，实时性好。
[0069]
6)将结果输出至原油管输调度控制系统，下游收油站可以对油品进行及时切割收油。
[0070]
为进一步论证本方法的可靠性，针对管道中多次输送不同种类的原油，存在多个混油界面的情况进行检测。与前述过程类似，连续计算光谱特征区内各波数点吸光度标准偏差一阶导数加和值，对多次输送不同种类原油的混油界面的判断结果如图4所示。图中清晰展现了3段混油周期曲线，显然，该方法均可实时、准确地对混合油品进行检测与跟踪。