CN106124448A - 一种原油切换状态下的常减压装置进料性质预测方法 - Google Patents

一种原油切换状态下的常减压装置进料性质预测方法 Download PDF

Info

Publication number
CN106124448A
CN106124448A CN201610409301.9A CN201610409301A CN106124448A CN 106124448 A CN106124448 A CN 106124448A CN 201610409301 A CN201610409301 A CN 201610409301A CN 106124448 A CN106124448 A CN 106124448A
Authority
CN
China
Prior art keywords
crude oil
crude
atmospheric
surge tank
oil
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.)
Granted
Application number
CN201610409301.9A
Other languages
English (en)
Other versions
CN106124448B (zh
Inventor
陈夕松
罗凡
方鑫
梅彬
王斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd
Original Assignee
NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd filed Critical NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd
Priority to CN201610409301.9A priority Critical patent/CN106124448B/zh
Publication of CN106124448A publication Critical patent/CN106124448A/zh
Application granted granted Critical
Publication of CN106124448B publication Critical patent/CN106124448B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3577Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
  • Feedback Control In General (AREA)
  • Vaporization, Distillation, Condensation, Sublimation, And Cold Traps (AREA)

Abstract

本发明提出一种原油切换状态下的常减压装置进料性质预测方法,首先获取当前加工原油和待切换原油的性质数据,然后以原油切换操作开始时间作为基准,基于实际工艺参数,每隔一定周期,对常减压装置进料性质进行预测。装置进料性质包括具有线性调合关系与非线性调合关系的原油性质。通过对常减压装置进料性质的准确预测,有利于控制系统和优化系统在原油切换操作期间及早应对进料性质的变化带来的扰动,从而提高装置的控制和优化性能,稳定产品质量。

Description

一种原油切换状态下的常减压装置进料性质预测方法
技术领域
本发明涉及一种原油性质的预测方法,尤其是在原油切换状态下的常减压装置原油进料预测方法。
背景技术
在炼化企业的生产过程中,同一套常减压装置往往需要通过原油切换操作,对不同的来料原油进行生产加工。但在实际生产中,原油切换操作较为频繁,切换操作期间,装置瞬时波动较大,控制性能降低,直接影响企业的产品质量。如果能够对原油切换操作期间常减压装置的进料性质进行准确预测,并将其作为前馈信号及时传递给控制系统,对装置操作参数提前进行合理调整,就能够有效应对此期间内进料性质的变化带来的扰动。
此外,当前加工原油和待切换原油性质的准确性也尤为关键。过去,大部分炼化企业采用传统的历史原油性质数据对常减压装置进行模拟和优化,历史原油数据本身往往不能准确表征原油的实时性质,因而利用历史原油数据计算的结果不够准确。目前,近红外原油快速评价技术已日趋成熟,利用该技术可以准确获取原油的性质数据,从而提高预测结果的精度。
发明内容
本发明提出一种原油切换状态下的常减压装置进料性质预测方法。首先利用原油近红外快速评价技术获取当前加工原油和待切换原油的性质数据,然后以原油切换操作开始时间作为基准,基于实际工艺参数,每隔一定周期,计算当前缓冲罐内的混合原油性质,并由此预测常减压装置进料性质。该方法包括以下步骤:
1)获取当前加工原油和待切换原油的性质数据;
2)初始化参数,包括原油切换操作开始时间、预测周期;
3)计算常减压装置前缓冲罐内当前的混合原油性质,即常减压装置进料性质;
4)判断混合原油性质是否接近待切换原油性质,是则停止,否则等待一个预测周期后,返回步骤3)。
上述过程中,利用近红外快速评价技术得到当前加工原油的性质数据的内容,详见发明专利《一种基于复合预测技术的原油性质快速检测方法》,专利号:CN 103364364A。
优选的,上述对常减压装置进料性质的预测包括具有线性调合关系与非线性调合关系的原油性质,其中具有线性调合关系的原油性质主要包括酸值、硫含量和API密度,具有非线性调合关系的原油性质主要包括原油密度和倾点。
优选的,上述计算混合原油的线性调合性质,依据线性调合规则,如式(1)所示:
f(k)=X1(k)×f+X2(k)×f(k-1) (1)
其中,f(k)为k时刻缓冲罐内混合原油的线性调合性质,f(k-1)为k-1时刻缓冲罐内混合原油的线性调合性质,f为待切换原油的线性调合性质,X1(k)为k时刻待切换原油的质量调合比,X2(k)为k时刻混合前缓冲罐内原油的质量调合比。
优选的,上述计算混合原油的密度,如式(2)所示:
ρ ( k ) = 1 X 1 ( k ) ρ + X 2 ( k ) ρ ( k - 1 ) - - - ( 2 )
其中,ρ(k)为k时刻缓冲罐内混合原油的密度,ρ(k-1)为k-1时刻缓冲罐内混合原油的密度,ρ为待切换原油的密度,X1(k)为k时刻待切换原油的质量调合比,X2(k)为k时刻混合前缓冲罐内原油的质量调合比。
优选的,上述混合原油的倾点采用下式计算:
P ( k ) = ( l n P I ( k ) - 1.85 ) / 0.042 - 32 1.85 - - - ( 3 )
其中,P(k)为k时刻缓冲罐内混合原油的倾点,PI(k)为k时刻缓冲罐内混合原油的倾点指数。
优选的,上述混合原油的倾点指数采用下式计算:
PI(k)=X1(k)e1.85+0.042×(P×1.8+32)+X2(k)e1.85+0.042×(P(k-1)×1.8+32) (4)
其中,PI(k)为k时刻缓冲罐内混合原油的倾点指数,P(k-1)为k-1时刻缓冲罐内混合原油的倾点,P为待切换原油的倾点,X1(k)为k时刻待切换原油的质量调合比,X2(k)为k时刻混合前缓冲罐内原油的质量调合比。
优选的,上述计算质量调合比,如式(5)、式(6)所示:
X 1 ( k ) = F T ρ A l ρ ( k - 1 ) - F T ρ ( k - 1 ) + F T ρ - - - ( 5 )
X 2 ( k ) = A l ρ ( k - 1 ) - F T ρ ( k - 1 ) A l ρ ( k - 1 ) - F T ρ ( k - 1 ) + F T ρ - - - ( 6 )
其中,F为原油流量,A为缓冲罐的截面积,l为缓冲罐内原油液位高度,ρ为待切换原油的密度,ρ(k-1)为k-1时刻缓冲罐内混合原油的密度,T为预测周期。
优选的,缓冲罐的横截面积A为已知数据,原油流量F、缓冲罐内原油液位高度l、预测周期T为初始化参数,待切换原油的密度ρ、k-1时刻缓冲罐内混合原油的密度ρ(k-1)通过近红外快速评价技术得到。
有益效果:
本发明提出一种原油切换状态下的常减压装置进料性质预测方法,通过利用原油近红外快速评价技术获取当前加工原油和待切换原油的性质数据,对常减压装置的进料性质进行准确预测,有利于控制系统和优化系统在原油切换操作期间及早应对进料性质变化带来的扰动,从而提高装置的控制和优化性能,稳定产品质量。
附图说明
图1原油加工工艺示意图。
图2原油切换状态下的常减压装置进料性质预测流程图。
具体实施案例
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
以某炼化企业为例,该企业常减压装置具有炼化企业典型工艺,首先通过机泵将组分原油输送至原油静态混合器进行调合,调合后的原油然后进入原油缓冲罐,最后输送至常减压装置进行加工。2016年4月15日调合组分油为杰诺、卡宾达和达尔混合油,在11时25分开始进行原油切换操作,调合组分油改为沙特重油、巴士拉申和南帕斯。
图1为原油加工工艺示意图。其中缓冲罐的截面积A=400m2
下面以酸值、原油密度和倾点为例,给出原油切换状态下的常减压装置进料性质预测的详细过程。
1)当前加工原油由杰诺、卡宾达和达尔混合油构成,三种组分油经过一定配比调合得到当前加工原油,其中酸值C0=1.0mgKOH/g、原油密度ρ0=872.69kg/m3、倾点P0=-7.4℃。
2)待切换原油为沙特重油、巴士拉申和南帕斯的混合油,三种组分油经过一定配比调合得到待切换原油,其中酸值C1=0.89mgKOH/g、原油密度ρ1=873.39kg/m3、倾点P1=-3.5℃。
3)原油切换开始时间为11时25分,此时循环迭代标识符k=0,原油流量F=744.18m3/h,液位稳定在l=12.4m,预测周期T=0.2h。
4)计算当前缓冲罐内的混合原油的酸值C(k)、密度ρ(k)以及倾点P(k):
k=0时
当前缓冲罐内的混合原油的酸值C(0)=C0=1.0mgKOH/g,密度ρ(0)=ρ0=872.69kg/m3,倾点P(0)=P0=-7.4℃。
k>0时
根据式(5)得到k时刻待切换原油的质量调合比
X 1 ( k ) = FTρ 1 A l ρ ( k - 1 ) - F T ρ ( k - 1 ) + FTρ 1 = 26.796 ρ ( k - 1 ) + 26.796
根据式(6)得到k时刻混合前缓冲罐内原油的质量调合比
X 2 ( k ) = A l ρ ( k - 1 ) - F T ρ ( k - 1 ) A l ρ ( k - 1 ) - F T ρ ( k - 1 ) + FTρ 1 = ρ ( k - 1 ) ρ ( k - 1 ) + 26.796
根据式(1)可得当前缓冲罐内的混合原油的酸值
C ( k ) = X 1 ( k ) × C 1 + X 2 ( k ) × C ( k - 1 ) = 23.848 + ρ ( k - 1 ) × C ( k - 1 ) ρ ( k - 1 ) + 26.796
根据式(2)可得当前缓冲罐内的混合原油的密度
ρ ( k ) = 1 X 1 ( k ) ρ 1 + X 2 ( k ) ρ ( k - 1 ) = 0.9702 ρ ( k - 1 ) + 25.998
根据式(4)可得当前缓冲罐内的混合原油的倾点指数
P I ( k ) = X 1 ( k ) e 1.85 + 0.042 × ( P 1 × 1.8 + 32 ) + X 2 ( k ) e 1.85 + 0.042 × ( P ( k - 1 ) × 1.8 + 32 ) = ρ ( k - 1 ) e 1.85 + 0.042 × ( P ( k - 1 ) × 1.8 + 32 ) + 501.52 ρ ( k - 1 ) + 26.796
根据式(3)可得当前缓冲罐内的混合原油的倾点
P ( k ) = ( l n P I ( k ) - 1.85 ) / 0.042 - 32 1.85
5)判断混合原油性质是否接近待切换原油性质,是则停止,否则k=k+1,返回步骤4)。
经过计算,在2016/4/15 16:45时刻得到的缓冲罐内的混合原油性质已接近切换原油性质,计算结束。各时刻缓冲罐内的混合原油的酸值、密度以及倾点数据如表1所示。
表1缓冲罐内的混合原油的酸值、密度、倾点数据
从表1可以看出,在原油切换状态下,常减压装置的进料性质由切换前原油性质逐步过渡到切换原油的性质,通过该方法可以对常减压装置的进料性质进行准确预测,有利于控制系统和优化系统在原油切换操作期间及早应对进料性质的变化带来的扰动,从而提高装置的控制和优化性能,稳定产品质量。

Claims (10)

1.一种原油切换状态下的常减压装置进料性质预测方法,其特征在于该方法包括以下步骤:
1)获取当前加工原油和待切换原油的性质数据;
2)初始化参数,包括原油切换操作开始时间、预测周期;
3)计算常减压装置前缓冲罐内当前的混合原油性质,即常减压装置进料性质;
4)判断混合原油性质是否接近待切换原油性质,是则停止,否则等待一个预测周期后,返回步骤3)。
2.根据权利要求1所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于当前加工原油和待切换原油的性质数据通过近红外快速评价技术得到。
3.根据权利要求1所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于预测的性质包括具有线性调合关系与非线性调合关系的原油性质。
4.根据权利要求3所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于所述具有线性调合关系的原油性质主要包括酸值、硫含量和API密度,具有非线性调合关系的原油性质主要包括原油密度和倾点。
5.根据权利要求4所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于混合原油的线性调合性质采用下式计算:
f(k)=X1(k)×f+X2(k)×f(k-1)
其中,f(k)为k时刻缓冲罐内混合原油的线性调合性质,f(k-1)为k-1时刻缓冲罐内混合原油的线性调合性质,f为待切换原油的线性调合性质,X1(k)为k时刻待切换原油的质量调合比,X2(k)为k时刻混合前缓冲罐内原油的质量调合比。
6.根据权利要求4所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于混合原油的密度采用下式计算:
ρ ( k ) = 1 X 1 ( k ) ρ + X 2 ( k ) ρ ( k - 1 )
其中,ρ(k)为k时刻缓冲罐内混合原油的密度,ρ(k-1)为k-1时刻缓冲罐内混合原油的密度,ρ为待切换原油的密度,X1(k)为k时刻待切换原油的质量调合比,X2(k)为k时刻混合前缓冲罐内原油的质量调合比。
7.根据权利要求4所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于混合原油的倾点采用下式计算:
P ( k ) = ( ln P I ( k ) - 1.85 ) / 0.042 - 32 1.85
其中,P(k)为k时刻缓冲罐内混合原油的倾点,PI(k)为k时刻缓冲罐内混合原油的倾点指数。
8.根据权利要求7所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于混合原油的倾点指数采用下式计算:
PI(k)=X1(k)e1.85+0.042×(P×1.8+32)+X2(k)e1.85+0.042×(P(k-1)×1.8+32)
其中,PI(k)为k时刻缓冲罐内混合原油的倾点指数,P(k-1)为k-1时刻缓冲罐内混合原油的倾点,P为待切换原油的倾点,X1(k)为k时刻待切换原油的质量调合比,X2(k)为k时刻混合前缓冲罐内原油的质量调合比。
9.根据权利要求5或6或8中任一项所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于质量调合比采用下式计算:
X 1 ( k ) = F T ρ A l ρ ( k - 1 ) - F T ρ ( k - 1 ) + F T ρ
X 2 ( k ) = A l ρ ( k - 1 ) - F T ρ ( k - 1 ) A l ρ ( k - 1 ) - F T ρ ( k - 1 ) + F T ρ
其中,F为原油流量,A为缓冲罐的横截面积,l为缓冲罐内原油液位高度,ρ为待切换原油的密度,ρ(k-1)为k-1时刻缓冲罐内混合原油的密度,T为预测周期。
10.根据权利要求9所述的一种原油切换状态下的常减压装置进料性质预测方法,其特征在于缓冲罐的横截面积A为已知数据,原油流量F、缓冲罐内原油液位高度l、预测周期T为初始化参数,待切换原油的密度ρ、k-1时刻缓冲罐内混合原油的密度ρ(k-1)通过近红外快速评价技术得到。
CN201610409301.9A 2016-06-12 2016-06-12 一种原油切换状态下的常减压装置进料性质预测方法 Active CN106124448B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610409301.9A CN106124448B (zh) 2016-06-12 2016-06-12 一种原油切换状态下的常减压装置进料性质预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610409301.9A CN106124448B (zh) 2016-06-12 2016-06-12 一种原油切换状态下的常减压装置进料性质预测方法

Publications (2)

Publication Number Publication Date
CN106124448A true CN106124448A (zh) 2016-11-16
CN106124448B CN106124448B (zh) 2019-02-01

Family

ID=57270246

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610409301.9A Active CN106124448B (zh) 2016-06-12 2016-06-12 一种原油切换状态下的常减压装置进料性质预测方法

Country Status (1)

Country Link
CN (1) CN106124448B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109634105A (zh) * 2018-12-11 2019-04-16 南京富岛信息工程有限公司 一种面向加工方案切换的原油动态调合方法
CN110726830A (zh) * 2019-10-23 2020-01-24 汉谷云智(武汉)科技有限公司 一种多进多出混合油性质预测方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101286065A (zh) * 2008-05-23 2008-10-15 清华大学 一种原油掺炼多周期优化的原油调合调度方法
CN102393636A (zh) * 2011-09-22 2012-03-28 任季明 一种对石油炼制过程的控制系统及方法
CN103365205A (zh) * 2013-07-25 2013-10-23 南京富岛信息工程有限公司 一种原油调合双调合头在线协调控制方法
CN105389648A (zh) * 2015-10-21 2016-03-09 南京富岛信息工程有限公司 一种常减压装置稳态工况的判别方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101286065A (zh) * 2008-05-23 2008-10-15 清华大学 一种原油掺炼多周期优化的原油调合调度方法
CN102393636A (zh) * 2011-09-22 2012-03-28 任季明 一种对石油炼制过程的控制系统及方法
CN103365205A (zh) * 2013-07-25 2013-10-23 南京富岛信息工程有限公司 一种原油调合双调合头在线协调控制方法
CN105389648A (zh) * 2015-10-21 2016-03-09 南京富岛信息工程有限公司 一种常减压装置稳态工况的判别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张栋栋等: "浅谈原油在线调合技术", 《内蒙古石油化工》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109634105A (zh) * 2018-12-11 2019-04-16 南京富岛信息工程有限公司 一种面向加工方案切换的原油动态调合方法
CN109634105B (zh) * 2018-12-11 2022-04-08 南京富岛信息工程有限公司 一种面向加工方案切换的原油动态调合方法
CN110726830A (zh) * 2019-10-23 2020-01-24 汉谷云智(武汉)科技有限公司 一种多进多出混合油性质预测方法及装置
CN110726830B (zh) * 2019-10-23 2022-03-04 汉谷云智(武汉)科技有限公司 一种多进多出混合油性质预测方法及装置

Also Published As

Publication number Publication date
CN106124448B (zh) 2019-02-01

Similar Documents

Publication Publication Date Title
CN106124448A (zh) 一种原油切换状态下的常减压装置进料性质预测方法
CN102375411A (zh) 一种油品在线调合方法及装置
CH703493B1 (de) Systeme und Verfahren zum Überwachen korrosiver Verunreinigungen, die mit flüssigem Brennstoff einhergehen.
Macaskill Computer simulation for mixed-model production lines
CN104950670B (zh) 一种连续搅拌釜式反应器的一体化多模型控制方法
CN107922855A (zh) 用于动态控制烃精制过程的预测分析的系统和方法
Zhong et al. Improvement and application of earned value analysis in coal project management
CA3100786A1 (en) Distillate blending system with online derived and/or indicated cetane number analyzer and related methods
CN104677721A (zh) 铸钢材质力学性能在线预测方法
CN103116315B (zh) 数控系统中实现刀具半径补偿全局干涉的控制方法
CN111257170A (zh) 一种解决混凝土生产过程中坍落度变化的间接观测方法
CN103365208A (zh) 一种具有底油性质偏差修正的原油调合优化方法
CN105400929A (zh) 一种kr终点硫含量的控制方法
CN108131566A (zh) 一种稠油在线输送系统及方法
CN103599713A (zh) 减水剂自动复配排料系统及其配料方法
Rusanescu et al. Mathematical Modelling of the Stress-Strain Curve for 31VMn12 Ecological Steel
CN101770532A (zh) 一种提高电路仿真器收敛性的方法
CN201162541Y (zh) 用于控制向尾砂中补加水的控制设备
CN201156185Y (zh) 用于控制充填料浆浓度的模糊控制设备
CN109589632B (zh) 一种剥离液的回收方法
CN103544549B (zh) 一种水泥沥青砂浆弹性模量的预测方法
US11221406B2 (en) Guided wave radar for consumable particle monitoring
CN102507476A (zh) 塑料样品中镉的测量方法
CN104007658B (zh) 一种具有双重循环优化的汽油在线调合方法
CN106501080A (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
GR01 Patent grant
GR01 Patent grant