[发明专利]一种高性能汽油调合非线性优化方法

摘要

本发明公开了一种高性能汽油调合非线性优化方法，采用线性调合并结合非线性调合规则，提高汽油加工企业成品油调合精度。该方法支持调合后汽油性质含量与目标偏差最小和经济最优相结合的综合指标，且支持汽油性质含量在范围内定值、向下限优化、向上限优化及范围内优化四种优化方式，提高了生产的灵活性，系统自动生成汽油调合配方，下发到控制系统执行。该方法可为炼厂汽油调合提供满足国标要求的汽油调合配方，避免质量过剩，提高了生产效率和经济效益，特别适用于炼油企业的汽油加工。

主权项

一种高性能汽油调合非线性优化方法，通过求解目标函数来获得各组分油占比，其特征在于该目标函数既考虑汽油参与优化的不同性质的调合值在优化范围内与优化目标值偏差最小，同时考虑生产成本最小化，具体包括以下步骤：a、进行本次调合参数初始化：设定调合总量，调合总流速的上限和下限，各掺炼线流速的上限和下限，参与调合的组分油种类，各组分油量的上限，各组分油占比的上限和下限；b、选择汽油参与优化的性质种类、调合规则、各性质的优化范围、优化方式、优化目标值和权重，组分油的价格目标和权重，不同性质种类选择不同的调合规则；汽油参与优化的性质在硫含量、氧含量、硫醇硫含量、烯烃含量、芳烃含量、苯含量、研究法辛烷值、马达法辛烷值、抗暴指数和雷德蒸汽压中选择一到数种，采用不同的调合规则：硫含量、氧含量和硫醇硫含量，采用线性质量调合规则如式a=Σi=1maixi---(2)]]>式(2)中，a为汽油相应性质的调合值，ai为组分油i的相应性质的参数，m为参与调合的组分油的个数，xi为组分油i的质量占比；烯烃含量、芳烃含量和苯含量采用线性体积调合规则如式a=Σi=1maivi---(3)]]>式(3)中，a为汽油相应性质的调合值，ai为组分油i的相应性质的参数，m为参与调合的组分油的个数，vi为组分油i的体积占比；研究法辛烷值、马达法辛烷值和抗暴指数采用相互作用法调合规则如式a=Σi=1maivi+Σi=1m-1viΣj=i+1mci,jvj---(4)]]>式(4)中，a为汽油相应性质的调合值，ai为组分油i的相应性质的参数，m为参与调合的组分油的个数，ci,j为两种组分油之间的调合效应系数，vi为组分油i的体积占比；研究法辛烷值、马达法辛烷值和抗暴指数或者采用Stewart法调合规则如式a=Σi=1m(viDi(ai+0.13Pi))/Σi=1mviDi---(5)]]>式(5)中，a为汽油相应性质的调合值，ai为组分油i的相应性质的参数，m为参与调合的组分油的个数，Pi为调合汽油中组分油i体积含量与不饱和烃的体积含量差，Di为根据Pi对应的权重指数，vi为组分油i的体积占比；其中，不饱和烃为烯烃和芳烃的集合；雷德蒸汽压采用Chevron法调合规则如式a=(Σi=1mai1.25vi)0.8---(6)]]>式(6)中，a为汽油相应性质的调合值，ai为组分油i的相应性质的参数，m为参与调合的组分油的个数，vi为组分油i的体积占比；式(3)‑(6)中，vi与xi的关系满足：vi=xiρiΣi=1mxiρi---(7)]]>其中，ρi为组分油i的密度；c、获取各掺炼线中组分油的当前性质；d、通过目标函数计算各组分油占比；e、占比下发到控制系统执行；f、判断本次调合是否完成：完成则结束，未完成则返回步骤c，所述目标函数为下式：minΦ(X)=Σj=1k[λj(f(X)j-goalj)]2+[λPrice(Σi=1mpriceixi-Pricemin)]2s.t.ΣX=1Xmin≤X≤XmaxFT=F1/x1=F2/x2=...=Fm/xmFTmin≤FT≤FTmaxFmin≤FT·X≤Fmaxt·FT·X≤MrangeLj≤f(X)j≤rangeUj---(1)]]>式(1)中，Φ(X)为目标函数；k表示汽油参与优化的性质种类个数；f(X)j为汽油参与优化的第j种性质所采用的调合规则，不同种性质采用不同的调合规则；m表示参与调合的掺炼线个数，即组分油个数；i表示第i条掺炼线，即第i个组分油；X为各组分油的质量占比，X＝[x1,x2…xm]，X≥0；xi表示第i个组分油的质量占比；Xmin和Xmax分别为各组分油占比的上限和下限，Xmin≥0，Xmin＝[xmin,1,xmin,2…xmin,m]，Xmax≥0，Xmax＝[xmax,1,xmax,2…xmax,m]；λj为汽油参与优化的第j种性质的权重，λj≥0；goalj为汽油参与优化的第j种性质的优化目标值，goalj≥0；Price为各组分油的价格，Price≥0，Price＝[price1,price2,…pricem]≥0；pricei表示第i个组分油的价格；Pricemin为各组分油的价格最低值，λPrice为价格权重，λPrice≥0；FT为当前调合总流速，FT≥0，FTmax为调合总流速的上限，FTmin为调合总流速的下限；Fmax为各掺炼线流速的上限，Fmax≥0，Fmax＝[Fmax,1,Fmax,2…Fmax,m]，Fmin为各掺炼线流速的下限，Fmin≥0，Fmin＝[Fmin,1,Fmin,2…Fmin,m]；t为本次调合预计执行时间，t≥0；M为各组分油可供调合总量，即各组分油量的上限，M＝[m1,m2…mm]，M≥0；F为各组分油的流速，F＝[F1,F2…Fm]；rangeLj为汽油参与优化的第j种性质的目标下限，rangeUj为汽油参与优化的第j种性质的目标上限。