CN113408194B - 常减压装置的通盘优化方法 - Google Patents

常减压装置的通盘优化方法 Download PDF

Info

Publication number
CN113408194B
CN113408194B CN202110647013.8A CN202110647013A CN113408194B CN 113408194 B CN113408194 B CN 113408194B CN 202110647013 A CN202110647013 A CN 202110647013A CN 113408194 B CN113408194 B CN 113408194B
Authority
CN
China
Prior art keywords
optimization
variables
model
variable
atmospheric
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.)
Active
Application number
CN202110647013.8A
Other languages
English (en)
Other versions
CN113408194A (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.)
Beijing Scienco Technology Co ltd
Beijing Yineng Gaoke Technology Co ltd
Original Assignee
Beijing Scienco Technology Co ltd
Beijing Yineng Gaoke Technology 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 Beijing Scienco Technology Co ltd, Beijing Yineng Gaoke Technology Co ltd filed Critical Beijing Scienco Technology Co ltd
Priority to CN202110647013.8A priority Critical patent/CN113408194B/zh
Publication of CN113408194A publication Critical patent/CN113408194A/zh
Application granted granted Critical
Publication of CN113408194B publication Critical patent/CN113408194B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Mathematical Physics (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明属于石油炼化领域,公开了一种常减压装置的通盘优化方法,该方法包括:确定常减压装置模拟的范围和涉及的变量;建立常减压装置的优化模型,包括优化目标方程和约束条件;通过物理逻辑重建方法对选取的变量进行优化求解,得到通盘优化所需的初值;在优化平台根据编程规则编写优化程序,给各变量代入所述初值,选用非线性求解器和制定的优化策略对优化模型进行求解,使优化值跳出局部最优,最终获得目标方程全局最优解及对应的优化变量的调整方向和幅度。本发明同时考虑多个现场操作变量的变化对优化目标和复杂现场约束的影响,避免了过分挖掘短期优化效益而造成长周期运行的重大操作风险和经济损失,对装置的长期稳定运行具有积极的效果。

Description

常减压装置的通盘优化方法
技术领域
本发明涉及石油炼化技术领域,尤其是涉及一种常减压装置的通盘优化方法。
背景技术
炼化行业生产决策方式,从以往的依靠经验发展到目前的工艺专家建立机理模型测算,经历了两个阶段的变化。
凭经验决策阶段:缺乏测算工具和模型,计划人员依据统计数据和经验来为生产计划分配物流走向;操作员根据自己的操作经验和操作手册,结合现场的生产情况给出操作参数调整方向,并通过化验分析结果一步步小幅度调整。
工艺专家用机理模型测算指导阶段:通过工况标定,建立和校正机理模型,再由工艺专家根据专业知识提出的优化思路,利用机理模型模拟测算或者设置简单的优化策略进行测算来验证,最后得到生产计划和操作方案用于指导生产。流程模拟软件20世纪50年代问世,经过几十年的发展,在业界得到了广泛的应用,并培养了一批优秀的建模和工艺优化专家。这些工艺专家利用机理模型和专业知识测算装置加工方案和生产计划,为生产提供指导,取得了一定的经济效益。
但是,常减压装置的生产工况常常处于动态调整过程中,目前依赖专家和模拟软件进行工艺优化,这远远跟不上变化的需求。企业需要常减压装置优化的整体解决方案,以进一步提升企业原油加工产业链的决策管理水平,将企业的经营目标转化成生产操作中的操作目标。
常规优化思路通常一次仅考虑单个操作变量产生的影响,很难考虑多个操作变量同时调整对整个装置各个方面造成的影响,如何放眼到整个装置的运行周期从生产或经济效益的角度出发,快速获得一套对装置当前状态经过全局考虑的可行且准确的优化方案,是目前亟待解决的技术问题。
发明内容
针对现有技术中存在的上述问题,本发明的目的是提供一种常减压装置的通盘优化方法,其同时考虑多个操作变量的调整对系统或设备的影响,兼顾装置当前状态和装置长周期稳定运行需求,提出经过全局考虑的优化方案,进而使装置达到节能增效和长周期运行的目的。
为解决上述技术问题,本发明采用的技术方案如下:
一种常减压装置的通盘优化方法,包括:
确定常减压装置囊括的范围和涉及的变量;
建立常减压装置的优化模型,包括优化目标方程和约束条件;
通过物理逻辑重建方法获得选取的变量在所需工况下的数据,得到通盘优化所需的初值、原油性质和装置特性参数;
在优化平台根据编程规则编写优化程序,给各变量代入所述初值,选用非线性求解器和制定的优化策略对优化模型进行求解,使优化值跳出局部最优,最终获得目标方程全局最优解及对应的优化变量的调整方向和幅度。
进一步地,所述的优化目标方程的建立包括如下过程:
确定优化目标,其中,所述的优化目标包括经济效益最大化、产品产量最大/最小化、公用工程消耗最小化,或者它们的均衡组合;
建立优化目标方程。
进一步地,所述的约束条件的建立包括如下过程:
通过降维方法,选择模型的输入输出变量;
采用机器学习方法,训练某个输出变量与对应的输入变量间的关系模型;采用机器学习方法和经验机理混合算法,建立具有多个神经元网络群的常减压装置的AI模型,用于常减压装置优化模型中的等式约束条件;
采用统计方法对历史数据进行统计,获得变量允许的变化范围,从而获得常减压装置优化模型中不等式约束条件,结合上一步形成的等式约束条件,形成了常减压装置优化模型中的约束条件。
进一步地,在每一步的优化方向搜索中,对选取的变量对整个装置优化目标和潜在瓶颈产生的影响进行统一测算,其中,所述的变量包括各个塔的处理量、常压炉及减压炉出口温度、各处汽提蒸汽流量、各处回流量和产品产量。
进一步地,在每一步的优化方向搜索中,基于设定的现场约束对优化模型同时进行测算是否满足生产要求,其中,所述的现场约束包括物料平衡、能量平衡、相平衡、热量传递和设备性能五大现场约束条件。
进一步地,所述的约束条件为:
yn=fn(F,T,P,yp) n,p∈N,n≠p
Propm=fm(F,T,P,yn)
yn,L≤yn≤yn,U
ΔFmin≤ΔF≤ΔFm,ax
ΔTmin≤ΔT≤ΔTmax
ΔPmin≤ΔP≤ΔPmax
Propm,minspec+Propm,accuracy≤Propm≤Propm,maxspec-Propm,accuracy
其中,F、T、P分别表示的是在优化中的流量、温度和压力可调自变量,yn表示除了侧线产品中工艺指标要求的物性和可调自变量以外的变量,yp表示除了侧线产品中工艺指标要求的物性、可调自变量和yn以外的变量;f代表了物料平衡、能量平衡、相平衡、压力平衡、热量传递多个物理模型;N表示1到N的集合,即1,2,3……N;min表示最小值,max表示最大值;Propm表示的是侧线产品中工艺指标要求的物性,Propm,minspec,Propm,maxspec分别表示指标要求的上下限,Propm,accuracy表示这个物性的模型预测精度;ΔF,ΔT,ΔP是变量幅度,相关不等式是对优化幅度的约束。
进一步地,对优化变量按照属性进行调整频率分类,其中,将变量的频率分为高频次、中频次和低频次三类,调整级别分为三级,3级调整包括所有变量,2级调整包括中、低频次变量,1级调整包括低频次变量,根据调整级别涉及的变量类别确定优化变量,不在调整级别内的变量则固定为常量。
进一步地,采用广义既约梯度算法优化求解,在每一步优化方向搜索中都基于高精度ANN模型同时进行测算是否满足设备性能,包括:分馏塔塔板水力学性能、换热器传热能力、输送设备可操作范围、加热炉负荷瓶颈、设备高温位结焦倾向。
进一步地,在使用所述优化模型前,还判断当前工况是否超出了模型范围,若超出范围,会在当前工况附近范围进行工况增强,新增部分样本,将新增的所述部分样本信息汇总到之前的样本中对模型进行校正训练。
与现有技术相比,本发明所述的一种常减压装置的通盘优化方法,具有如下技术效果:
1、本发明采用通盘优化方法,以物理逻辑重建的结果为基础,综合考虑现场的实际可行性约束,在可控的优化策略下,以提高系统经济效益为主要目标,利用数学规划算法,优化求解得到操作条件变化的方向和大小,从而为下一过程的调整步骤生成提供目标,进而起到节能增效的目的。
2、优化变量全面,综合考虑现场所有变量对装置的影响,从而更充分的挖掘节能增效空间。
3、从实际出发,根据现场操作方法制定优化策略,使优化结果更易实施,避免与现场脱钩。
4.建立模型的数据量大,且均来自现场数据或严格机理模拟数据,模型可自学习,不断学习装置特性,从而提高模型表征装置的准确度。
5、通过物理逻辑重建获得原油性质,从而确保优化的准确度。
6、通过物理逻辑重建获得优化的初值,将优化时间减少至5分钟以内。
7、通过模型来评估设备运行现状,避免过分挖掘短期优化效益而造成长周期运行的重大操作风险和经济损失,从而满足装置长期稳定运行需求。
附图说明
图1为本发明实施例所述的常减压装置的通盘优化方法的流程示意图。
具体实施方式
以下结合附图对本发明作进一步详细说明,但不作为对本发明的限定。
参照图1所示,本发明实施例所公开的一种常减压装置的通盘优化方法,包括如下步骤:
步骤S1:确定装置的覆盖范围和涉及的变量。常减压装置主要包括三个精馏塔(初馏塔、常压塔和减压塔)、两个加热炉(常压炉和减压炉)、一个电脱盐和一个换热网络,变量包括各个塔的处理量、常压炉及减压炉出口温度、各处汽提蒸汽流量、各中段回流量、产品产量和换热器旁路流量等。
步骤S2:建立常减压装置优化模型,包括优化目标方程和约束条件。
其中,优化目标方程建立包括如下步骤:
步骤S21:确定优化目标,包括经济效益最大化、产品产量最大/最小化、公用工程消耗最小化等;
步骤S22:建立优化目标方程。
通盘优化的目标是根据项目需求实现所需的输出,比如经济效益最大化、产能最大化、公用工程最小化等,或者它们的均衡组合,或者设定的其它任一可行的目标方案,以在可操作区域内最大化经济效益为例,目标函数实质为最大收益函数,经济效益中一般包含三部分内容:各侧线产品收益、原料成本(在固定原料结构和处理量的情况下,原料成本一般固定)、以及公用工程成本,故最大收益函数实质是求解出各侧线产品收益减去原料成本及公用工程成本后的结果值。因此影响经济效益的变量包括原料加工量、产品的产量和公用工程的消耗量,那么该目标的方程可表示为:
Figure BDA0003110310320000061
其中,p是价格,Fi为产品i的产量,Fj为原料j加工量,Fk为公用工程k用量。
其中,常减压装置优化约束条件包括等式约束和不等式约束。等式约束采用人工神经元网络(ANN)训练和经验机理模拟结合获得。ANN模型模拟人脑神经元网络结构以及对信息的处理方式,采用大量的神经元之间相互联接构成的多层方程结构实现对输入信息的计算输出,是一种结构复杂但计算速度快速的模型,可模拟非常复杂并且参数众多的各种过程。另外ANN模型具备自学习功能,因此非常适用于本身工况多变的常减压装置。当工况出现较大变化时,随着现场新数据的输入,ANN模型具备自我校正功能,可逐步对模型参数进行调整以适应新工况,模拟精度也能得到逐步改善。
建立常减压装置优化约束条件的具体步骤如下:
步骤S23:通过降维方法,选择模型的输入输出变量;
步骤S24:采用机器学习方法,训练某个输出与其输入变量间的关系模型,通过改变机器学习中的参数,使训练获得的模型既能较好拟合训练数据又能较好预测测试数据;
步骤S25:若模型精度达不到要求,则通过反复试验方法,改变该输出对应的输入变量,继续采用步骤S24训练模型,直至精度达到要求;
步骤S26:采用机器学习方法和经验机理混合算法,建立具有多个神经元网络群的常减压装置的AI模型。具体获得产品性质软仪表预测模型、原油/中间流股性质软仪表预测模型、操作条件软仪表预测模型、流股单位焓值预测模型、塔板气液相负荷软仪表预测模型、减压炉和转油线相关温度和压力预测模型和换热网络模型,结合质量模型、能量平衡模型和减压塔全塔压降相关的经验模型,形成了常减压装置AI模型,用于常减压装置优化模型中的等式约束条件。
步骤S27:采用统计方法对历史数据进行统计,获得变量允许的变化范围,从而获得常减压装置优化模型中不等式约束条件,结合上一步形成的等式约束条件,形成了常减压装置优化模型中的约束条件。
举例来说,首先,根据某常减压装置的特点选出286个产品性质、93个原油/中间流股性质、42个操作条件、124个塔板气液相负荷、46个流股单位、减压炉和转油线相关温度和压力和98个换热温度作为目标变量,在数十万套运行数据基础上,通过降维方法,选出影响目标变量的输入变量。例如,产品性质会作为模型的目标变量。将其它变量包括原料组成、设备特性、操作条件等用降维方法分析对目标变量的影响,以此确定对目标变量影响较大的一些变量。再采用机器学习方法对目标变量和影响目标变量的变量样本数据进行训练,获得目标变量的AI模型。
最终根据需求训练了286个产品性质软仪表预测模型、93个原油/中间流股性质软仪表预测模型、42个操作条件软仪表预测模型、124个塔板气液相负荷软仪表预测模型、98个带有详细换热器结构的换热网络模型、46个流股单位焓值预测模型、2个减压炉和转油线相关温度和压力预测模型。结合上述模型和2个减压塔全塔压降相关模型、6个全塔物料平衡模型、3个全塔能量平衡模型形成了整个常减压装置优化模型的约束条件。
优化约束方程示意如下:
yn=fn(F,T,P,yp) n,p∈N,n≠p
Propm=fm(F,T,P,yn)
yn,L≤yn≤yn,U
ΔFmin≤ΔF≤ΔFmax
ΔTmin≤ΔT≤ΔTmax
ΔPmin≤ΔP≤ΔPmax
Propm,minspec+Propm,accuracy≤Propm≤Propm,maxspec-Propm,accuracy
其中,F、T、P分别表示的是在优化中的流量、温度和压力可调自变量,yn表示除了侧线产品中工艺指标要求的物性和可调自变量以外的变量,yp表示除了侧线产品中工艺指标要求的物性、可调自变量和yn以外的变量;f代表了物料平衡、能量平衡、相平衡、压力平衡、热量传递多个物理模型;N表示1到N的集合,即1,2,3,…,N;min表示最小值,max表示最大值;Propm表示的是侧线产品中工艺指标要求的物性;Propm,minspec;Propm,maxspec分别表示指标要求的上下限,Propm,accuracy表示这个物性的模型预测精度;ΔF,ΔT,ΔP是变量幅度,相关不等式是对优化幅度的约束。
步骤S3:采用非线性求解器对物理逻辑重建模型进行优化求解,获得表征所需工况的所有参数,一方面获取通盘优化所需的初值,另一方面获取优化时采用的原油性质和装置性能参数。
具体来说,以常减压装置AI模型为约束,以处理后的稳态数据为输入,采用数学规划法求解各变量预测值和测量值差值最小时各变量的优化值,用此方法获得的变量数据实现了历史工况数据和实时工况数据的校正,以及获取变量的软测量,包括原油性质、装置特性参数,从而还原现场工况,此方法获得的变量数据满足优化模型中约束条件,用该套数据作为通盘优化的初值可将优化时间缩减至5分钟以内。
物理逻辑重建模型示意如下:
Figure BDA0003110310320000091
yj=fj(xcal,i,yj′) j,j′∈J,j≠j′
xcal,i=fi(xcal,i′,yj) j,j′∈I,i≠i′
yj,L≤yj≤yj,U
xi,L≤xi,cal≤xi,U
其中,xmsd,i为经过数据处理之后的现场测量仪表数据;xsyserr,i为仪表i系统误差;xcal,i为最后计算输出的数据;xcal,i’为除了xcal,i之外的最后计算输出的数据;wi为仪表i的权重;Maxi为仪表i最大值;Mini为仪表i最小值,yj为其他变量;包括除了有现场测量之外的输出变量和中间变量;yj,为其他变量,包括除了有现场测量和yi之外的输出变量和中间变量;函数f代表所有x和y的关系,主要反映物料平衡、能量平衡、相平衡、压力平衡、热量传递等物理模型;J表示1到J的集合,即1,2,3,…,J;I表示1到I的集合,即1,2,3,…,I;yj,L和yj,U分别为yj的下限和上限,xi,L和xi,U分别为xi的下限和上限。
步骤S4:把优化变量根据调整的频次分为低、中、高频3类,如表1所示,调整级别分为1、2、3级,根据调整级别选择要调整的变量类别,如表2所示。本发明根据变量属性,设置不同的调节频率:对于不易进行反复调整和对系统影响较大的参数,尽量减少调节,比如加热炉出口温度变动对整个系统的变动影响较大,故设置低频次调节;对于一些需要经常进行调整或是有自动化手段辅助的操作参数,可高频次调节,如侧线抽出温度、常压塔塔底吹汽等;也有一些属于中等频次调节的,如减压塔塔底吹汽、中段回流流量的操作。对于1级调整,仅可以调整高频变量,中频变量和低频变量尽量不调整,对于2级调整,高频变量和中频变量均可以调整,低频变量尽量不调整,对于3级调整,高、中、低频均可以调整,其中,针对需要调整的参数,每一次的调整幅度控制在20%以内较佳。针对不同变量,进行不同频次的调整以及不同的调整策略,进而实现系统的长周期稳定运行。
表1变量优化类别
Figure BDA0003110310320000111
表2优化变量调整策略
Figure BDA0003110310320000112
步骤S5:根据优化平台编程规则依次定义前述步骤中涉及的变量、优化模型,并给变量赋前述计算出的初值,选用非线性求解器和制定的优化策略对模型进行优化求解。其中,优化平台可选用GAMS建模平台,求解器基于广义既约梯度(Generalized ReducedGradient,GRG)算法,在每一步操作变量调整时都综合评估所有变量的变化对系统所有潜在瓶颈的影响,确保优化方向的可行性。在每一步优化方向搜索中,多个常减压装置的变量,包括各个塔的处理量、常压炉、减压炉出口温度、各处汽提蒸汽流量、各处回流量等对整个系统优化目标和潜在瓶颈产生的影响都统一进行测算,也即对变量同时进行调整,而该调整也需要满足相互间的约束关系,进而评估产生满足装置操作可行条件下的最佳优化方向和幅度。同时,为了确保优化方案可行,五大类现场约束在每一步优化方向搜索中都基于高精度ANN模型同时进行测算是否满足系统的性能要求,包括分馏塔塔板水力学性能、换热器传热能力、输送设备可操作范围、加热炉负荷瓶颈、设备高温位结焦倾向。
通过灵敏度分析选取影响装置性能较大的关键变量,通过对关键变量分段将求解空间划分为几个部分,在每个空间进行优化求解,比较所有空间的结果,选择最优结果,以此避免优化结果陷入局部最优。
步骤S6:模型的自学习和自校正。现场的设备性能和加工的原料发生较大变化时,可能会导致工况超出上述AI模型覆盖的范围,模型对范围以外的预测效果比较差。因此在使用该模型之前,本发明还判断当前工况是否超出了模型范围,若超出范围,会在当前工况附近范围进行工况增强,新增部分样本,将新增的所述部分样本信息汇总到之前的样本中对模型进行校正训练,实现模型的自学习和自校正。另外通过物理逻辑重建计算结果获得装置当前特性如塔板板效和转油线压降等,更新模型中装置特性参数,使优化模型自适应装置特性变化。
本发明实施例中的优化约束主要是考虑到物料平衡、能量平衡、相平衡、热量传递和设备性能五大因素;变量的一些调整幅度的限制,如产品指标约束范围,需要考虑模型精度所引起的冗余量、关键塔板内回流变化量、三个塔处理量的变化幅度、去下游装置的侧线流量变化幅度、分支流量的变化幅度、减压塔各填料段的喷淋密度范围。
本实施例中,通盘优化后的结果主要包括:优化前后的总效益变化、优化前后产品效益变化、优化前后能耗变化、优化前后操作参数变化、优化前后各侧线产品性质变化、优化前后塔板水力学变化等。
本实施例中,不仅能针对装置当前状态提出经过全局考虑的优化方案,还能够放眼整个装置运行周期,对装置的长周期稳定运行进行优化。例如:短期来看,减压塔洗涤油量减少有利于提高减压塔拔出率。然而,长期来看过小的洗涤油量容易造成减压塔下部结焦,降低减压塔设备性能,从而导致更大的经济损失。通盘优化技术通过ANN模型来评估设备高温位结焦倾向,避免过分挖掘短期优化效益而影响长周期运行的重大操作风险。
本发明利用机理模型产生大量数据,以弥补生产数据范围窄而无法全面模拟生产装置的不足,该方法建立的模型能和机理模型一样表达出生产装置的特性,在计算方面比机理模型计算速度快;且自主学习装置具有实时的特性,根据装置实时特性修正模型,从而提高模型的准确度。
本发明的通盘优化是以物理逻辑重建的结果为基础,物理逻辑重建的结果是表征当前工况的所有参数数据,该结果将作为通盘优化的初始值,并给通盘优化提供原油性质和装置性能参数数据,从而确保快速准确地获得优化结果。
本发明根据还原的现场工况确定优化方向和幅度,避免了与现场工况脱钩,同时考虑调整多个变量对效益的影响,符合常减压真实的运行逻辑,能够及时地搜索出最佳的参数优化方向,使得经济效益最大化,不仅能针对装置当前状态提出经过全局考虑的优化方案,还能够放眼整个装置运行周期,对装置的长周期稳定运行进行优化,对生产具有较高的实用性。
上述说明示出并描述了本发明的若干优选实施例,但如前所述,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。

Claims (9)

1.一种常减压装置的通盘优化方法,其特征在于,该方法包括:
确定常减压装置囊括的范围和涉及的变量;
建立常减压装置的优化模型,包括优化目标方程和约束条件;
通过物理逻辑重建方法获得选取的变量在所需工况下的数据,得到通盘优化所需的初值、原油性质和装置性能参数;
在优化平台根据编程规则编写优化程序,给各变量代入解出的所述初值,选用非线性求解器和制定的优化策略对优化模型进行求解,使优化值跳出局部最优,最终获得目标方程全局最优解及对应的优化变量的调整方向和幅度;
其中,所述物理逻辑重建方法为采用非线性求解器对物理逻辑重建模型进行优化求解;
所述物理逻辑重建模型表达如下:
Figure FDA0003812931040000011
yj=fj(xcal,i,yj′) j,j′∈J,j≠j′
xcal,i=fi(xcal,i′,yj) i,i′∈I,i≠i′
yj,L≤yj≤yj,U
xi,L≤xi,cal≤xi,U
其中,xmsd,i为经过数据处理之后的现场测量仪表数据;xsyserr,i为仪表i系统误差;xcal,i为最后计算输出的数据;xcal,i’为除了xcal,i之外的最后计算输出的数据;wi为仪表i的权重;
Figure FDA0003812931040000012
为仪表i最大值;Mini为仪表i最小值,yj为其他变量;包括除了有现场测量之外的输出变量和中间变量;yj’为其他变量,包括除了有现场测量和yj之外的输出变量和中间变量;函数f代表所有x和y的关系,主要反映物料平衡、能量平衡、相平衡、压力平衡、热量传递的物理模型;J表示1到J的集合,即1,2,3,…,J;I表示1到I的集合,即1,2,3,…,I;yj,L和yj,U分别为yj的下限和上限,xi,L和xi,U分别为xi的下限和上限。
2.根据权利要求1所述的通盘优化方法,其特征在于,所述的优化目标方程的建立包括如下过程:
确定优化目标,其中,所述的优化目标包括产品产量最大/最小化、公用工程消耗最小化,或者它们的均衡组合;
建立优化目标方程。
3.根据权利要求2所述的通盘优化方法,其特征在于,所述的约束条件的建立包括如下过程:
通过降维方法,选择模型的输入输出变量;
采用机器学习方法,训练某个输出变量与对应的输入变量间的关系模型;
采用机器学习方法和经验机理混合算法,建立具有多个神经元网络群的常减压装置的AI模型,用于常减压装置优化模型中的等式约束条件;
采用统计方法对历史数据进行统计,获得变量允许的变化范围,从而获得常减压装置优化模型中不等式约束条件,结合上一步形成的等式约束条件,形成了常减压装置优化模型中的约束条件。
4.根据权利要求1所述的通盘优化方法,其特征在于,在每一步的优化方向搜索中,对选取的变量对整个装置优化目标和潜在瓶颈产生的影响进行统一测算。
5.根据权利要求4所述的通盘优化方法,其特征在于,在每一步的优化方向搜索中,基于设定的现场约束对优化模型同时进行测算是否满足生产要求,其中,所述的现场约束包括物料平衡、能量平衡、相平衡、热量传递和设备性能五大现场约束条件。
6.根据权利要求5所述的通盘优化方法,其特征在于,所述的约束条件为:
yn=fn(F,T,P,yp) n,p∈N,n≠p
Propm=fm(F,T,P,yn)
yn,L≤yn≤yn,U
ΔFmin≤ΔF≤ΔFmax
ΔTmin≤ΔT≤ΔTmax
ΔPmin≤ΔP≤ΔPmax
Propm,minspec+Propm,accuracy≤Propm≤Propm,maxspec-Propm,accuracy
其中,F、T、P分别表示的是在优化中的流量、温度和压力可调自变量,yn表示除了侧线产品中工艺指标要求的物性和可调自变量以外的变量,yp表示除了侧线产品中工艺指标要求的物性、可调自变量和yn以外的变量;f代表了物料平衡、能量平衡、相平衡、压力平衡、热量传递多个物理模型;N表示1到N的集合,即1,2,3……N;min表示最小值,max表示最大值;Propm表示的是侧线产品中工艺指标要求的物性,Propm,minspec,Propm,maxspec分别表示指标要求的上下限,Propm,accuracy表示这个物性的模型预测精度;ΔF,ΔT,ΔP是变量幅度,相关不等式是对优化幅度的约束。
7.根据权利要求3~6任一所述的通盘优化方法,其特征在于,对变量的优化按照属性进行调整频率分类,其中,对变量的频率分为高频次、中频次和低频次三类,调整级别分为三级,3级调整包括所有变量,2级调整包括中、低频次变量,1级调整包括低频次变量,根据调整级别涉及的变量类别确定优化变量,不在调整级别内的变量则固定为常量。
8.根据权利要求5或6所述的通盘优化方法,其特征在于,采用广义既约梯度算法优化求解,在每一步优化方向搜索中都基于高精度ANN模型同时进行测算是否满足设备性能,包括:分馏塔塔板水力学性能、换热器传热能力、输送设备可操作范围、加热炉负荷瓶颈、设备高温位结焦倾向。
9.根据权利要求1所述的通盘优化方法,其特征在于,在使用所述优化模型前,还判断当前工况是否超出了模型范围,若超出范围,会在当前工况附近范围进行工况增强,新增部分样本,将所述新增部分样本的信息汇总到之前的样本中对模型进行校正训练,实现模型自学习。
CN202110647013.8A 2021-06-10 2021-06-10 常减压装置的通盘优化方法 Active CN113408194B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110647013.8A CN113408194B (zh) 2021-06-10 2021-06-10 常减压装置的通盘优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110647013.8A CN113408194B (zh) 2021-06-10 2021-06-10 常减压装置的通盘优化方法

Publications (2)

Publication Number Publication Date
CN113408194A CN113408194A (zh) 2021-09-17
CN113408194B true CN113408194B (zh) 2022-12-02

Family

ID=77683408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110647013.8A Active CN113408194B (zh) 2021-06-10 2021-06-10 常减压装置的通盘优化方法

Country Status (1)

Country Link
CN (1) CN113408194B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106444428A (zh) * 2016-10-25 2017-02-22 中国石油化工股份有限公司 基于流程模拟软件的常减压装置优化操作系统及方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9014858B2 (en) * 2008-09-30 2015-04-21 Rockwell Automation Technologies, Inc. Energy optimizer for dehydrating biofuels through distillation towers and molecular sieves
CN108536108A (zh) * 2018-05-17 2018-09-14 杭州辛孚能源科技有限公司 一种常减压装置的实时优化方法
US11205028B2 (en) * 2018-09-06 2021-12-21 Terrafuse, Inc. Estimating physical parameters of a physical system based on a spatial-temporal emulator
CN109783916B (zh) * 2019-01-02 2021-06-18 大连理工大学 一种基于仿真技术的空气压缩机群组优化调度决策方法
CN111914382B (zh) * 2019-05-07 2023-04-21 宁波大学 一种基于代理模型的常减压装置的约束进化优化方法
CN111914381B (zh) * 2019-05-07 2023-04-25 宁波大学 一种基于kplsr模型的常减压装置操作优化方法
CN111860938B (zh) * 2020-06-01 2023-05-30 浙江中控技术股份有限公司 用于原油储运系统的全局调合调度优化方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106444428A (zh) * 2016-10-25 2017-02-22 中国石油化工股份有限公司 基于流程模拟软件的常减压装置优化操作系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于流程模拟的常减压装置过程操作与生产计划集成优化;董晓杨等;《化工学报》;20150115;第66卷(第01期);第237-243页 *

Also Published As

Publication number Publication date
CN113408194A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
CN111353656B (zh) 一种基于生产计划的钢铁企业氧气负荷预测方法
WO2016101182A1 (zh) 基于贝叶斯网络和极限学习机的区间型指标预报方法
CN107430398A (zh) 用于调谐工艺模型的系统和方法
CN102540879A (zh) 基于群决策检索策略的多目标评价优化方法
CN107464025A (zh) 基于大数据的供应链需求预测与智能决策方法
CN113420500B (zh) 智能常减压系统
CN103631140A (zh) 基于性能评估的焦炉加热燃烧过程火道温度自动调节方法
Wu et al. Integrated soft sensing of coke-oven temperature
CN112085156B (zh) 一种分压线损率预测及降损潜力评估方法及系统
CN107464017A (zh) 基于带时间差分贝叶斯网络的自适应软测量预测方法
CN103699947A (zh) 一种基于元学习的电力系统时变非线性负荷组合预测方法
CN117389236B (zh) 一种环氧丙烷生产过程优化方法及系统
CN103778466B (zh) 一种基于矢量误差的转炉炼钢温度建模预测方法及系统
CN104765347A (zh) 一种渣油延迟焦化过程中收率实时预测方法
CN114036452B (zh) 一种应用于离散型生产线的产能评估方法
Subramanian et al. A prediction-optimization framework for site-wide process optimization
CN111832215A (zh) 一种在线预测钢板组织性能的方法
Wang et al. A novel self-adaptive fractional grey Euler model with dynamic accumulation order and its application in energy production prediction of China
CN113408194B (zh) 常减压装置的通盘优化方法
CN111680452A (zh) 基于全要素数据挖掘的电网工程精准投资决策仿真方法
CN114936742A (zh) 一种供水系统调度代理决策方法
Al Jlibawi et al. Optimization of distribution control system in oil refinery by applying hybrid machine learning techniques
CN110766234B (zh) 基于信息融合的水泥冷却过程篦下压力预测方法
CN116662925A (zh) 一种基于加权稀疏神经网络工业过程软测量方法
CN110728031B (zh) 一种基于ann建模平衡复杂石化过程产能的多目标优化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20211202

Address after: 6 / F, block AB, ocean Guanghua international, No. 10, Jintong West Road, Chaoyang District, Beijing 100020

Applicant after: Beijing Yineng Gaoke Technology Co.,Ltd.

Applicant after: BEIJING SCIENCO TECHNOLOGY CO.,LTD.

Address before: 6 / F, block AB, ocean Guanghua international, No. 10, Jintong West Road, Chaoyang District, Beijing 100020

Applicant before: Beijing Yineng Gaoke Technology Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant