CN103955778A - 一种高性能汽油调合非线性优化方法 - Google Patents

Abstract

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

Description

一种高性能汽油调合非线性优化方法

技术领域

本发明涉及炼油企业汽油调合领域，尤其是汽油调合的优化方法，具体是一种高性能汽油调合非线性优化方法。

背景技术

随着我国环保要求的日益提高，要求汽油满足更加严格的排放标准。而与此同时，随着世界油价上涨，炼油企业为提高效益，必须实现质量卡边控制。上述问题要求汽油调合必须实现更高精度的控制。

为方便实施，现有的汽油调合技术主要采用线性的调合规则来近似，但调合组分之间存在复杂的调合效应，尤其是辛烷值，采用线性调合规则很难准确预测调合结果。部分汽油调合技术虽考虑了辛烷值的非线性问题，但还有其它性质，如雷德蒸汽压，也具有非线性特性。此外，部分性质有多个调合规则，如辛烷值有相互作用法、斯图尔特(Stewart)法等非线性调合规则，调合系统如何灵活地根据炼油企业自身情况选择不同调合规则是炼油企业现阶段亟需解决的问题。

发明内容

为解决现有技术存在的问题，本发明提出了一种高性能汽油调合非线性优化方法，该方法能够在满足设备工艺等约束条件的前提下，采用线性和非线性调合规则，支持多种优化方式，以汽油性质含量目标偏差最小且经济最优为目标进行优化调合，保证汽油质量指标的同时，提高企业经济效益。

一种高性能汽油调合非线性优化方法，通过求解目标函数来获得各组分油占比，该目标函数既考虑汽油参与优化的不同性质的调合值在优化范围内与优化目标值偏差最小，同时考虑生产成本最小化，具体包括以下步骤：

a、进行本次调合参数初始化：设定调合总量，调合总流速的上限和下限，各掺炼线流速的上限和下限，参与调合的组分油种类，各组分油量的上限，各组分油占比的上限和下限；一条掺炼线上输入一种组分油，组分油可通过调合组分罐输入掺炼线也可以从生产装置直接输入掺炼线参与调合；

b、选择汽油参与优化的性质种类、调合规则、各性质的优化范围、优化目标值、优化方式和权重，组分油的价格目标和权重，不同性质种类选择不同的性质规则；

c、获取各掺炼线中组分油的当前性质；

d、通过目标函数计算各组分油占比；

e、占比下发到控制系统执行；

f、判断本次调合是否完成：完成则结束，未完成则返回步骤c。

本发明中，权重的选择依据人工经验进行设定。

本发明中，所述目标函数为下式：

$\min \mathrm{\Φ}\left(X\right)=\underset{i=1}{\overset{k}{\mathrm{\Σ}}}{\left[{\mathrm{\λ}}_k(f{\left(X\right)}_k-{\mathrm{goal}}_k)\right]}^2+{\left[{\mathrm{\λ}}_{\mathrm{Price}}(\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{price}}_i{x}_i-{\mathrm{Price}}_{\min })\right]}^2$

s.t.∑X＝1

X_min≤X≤X_max

FT＝F₁/x₁＝F₂/x₂＝…＝F_m/x_m      (1)

FT_min≤FT≤FT_max

F_min≤FT·X≤F_max

t·FT·X≤M

rangeL_k≤f(X)_k≤rangeU_k

式(1)中，Φ(X)为目标函数；k表示汽油参与优化的性质种类个数，f(X)_k为汽油参与优化的第k种性质所采用的调合规则，不同种性质采用不同的调合规则；m表示参与调合的掺炼线个数，即组分油个数；X为各组分油的质量占比，X＝[x₁,x₂…x_m]，X≥0；X_min和X_max分别为各组分油占比的上限和下限，X_min＝[x_min,1,x_min,2…x_min,m]，X_max≥0，X_max＝[x_max,1,x_max,2…x_max,m]；λ_k为汽油参与优化的第k种性质的权重，λ_k≥0；goal_k为汽油参与优化的第k种性质的优化目标值，goal_k≥0；Price为各组分油的价格，Price≥0，Price＝[price₁,price₂,…price_m]≥0；Price_min为各组分油的价格最低值，λ_Price为价格权重，λ_Price≥0；FT为当前调合总流速，FT≥0，FT_max为调合总流速的上限，FT_min为调合总流速的下限；F_max为各掺炼线流速的上限，F_max≥0，F_max＝[F_max,1,F_max,2…F_max,m]，F_min为各掺炼线流速的下限，F_min≥0，F_min＝[F_min,1,F_min,2…F_min,m]；t为本次调合预计执行时间，t≥0；M为各组分油可供调合总量，即各组分油量的上限，M＝[m₁,m₂…m_m]，M≥0；F为各组分油的流速，F＝[F₁,F₂…F_m]；rangeL_k为汽油参与优化的第k种性质的目标下限，rangeU_k为汽油参与优化的第k种性质的目标上限，rangeL_k≤f(X)_k≤rangeU_k即汽油参与优化的第k种性质的优化范围。

本发明中，优化方式选择：

向下限优化，goal_k＝rangeL_k；

向上限优化，goal_k＝rangeU_k；

范围内定值优化，goal_k＝rangeF_k，rangeF_k为设定的汽油参与优化的第k种性质的优化目标值，rangeF_k∈(rangeL_k,rangeU_k)；

范围内无目标优化，需满足f(X)_k∈[rangeL_k,rangeU_k]，此时f(X)_k-goal_k＝0。

本发明中汽油参与优化的性质在硫含量、氧含量、硫醇硫含量、烯烃含量、芳烃含量、苯含量、研究法辛烷值、马达法辛烷值、抗暴指数和雷德蒸汽压中选择一到数种，采用不同的调合规则：

硫含量、氧含量和硫醇硫含量的调合规律符合质量叠加关系，采用线性质量调合规则如式

$a=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{a}_i{x}_i---\left(2\right)$

式(2)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，x_i为组分油i的质量占比；

烯烃含量、芳烃含量和苯含量的调合规律符合体积叠加关系，采用线性体积调合规则如式

$a=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{a}_i{v}_i---\left(2\right)$

式(3)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，v_i为组分油i的体积占比；

研究法辛烷值、马达法辛烷值和抗暴指数的调合规律并不满足简单的叠加关系，在调合过程中组分油之间存在着调合效应，采用相互作用法调合规则如式

$a=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{a}_i{v}_i+\underset{i=1}{\overset{m-1}{\mathrm{\Σ}}}{v}_i\underset{j=i+1}{\overset{m}{\mathrm{\Σ}}}{c}_{i,j}{v}_j---\left(4\right)$

式(4)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，c_i,j为两种组分油之间的调合效应系数，v_i为组分油i的体积占比；

研究法辛烷值、马达法辛烷值和抗暴指数或者采用Stewart法调合规则如式

$a=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}\left(v_i{D}_i(a_i+0.13P_i)\right)/\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{v}_i{D}_i---\left(5\right)$

式(5)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，P_i为调合汽油中组分油i体积含量与不饱和烃的体积含量差，D_i为根据P_i对应的权重指数，v_i为组分油i的体积占比；其中，不饱和烃为烯烃和芳烃的集合，权重指数的选择依据人工经验进行设定；

雷德蒸汽压采用Chevron法调合规则如式

$a={\left(\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{a}_i^{1.25}{v}_i\right)}^{0.8}---\left(6\right)$

式(6)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，v_i为组分油i的体积占比；

式(3)-(6)中，v_i与x_i的关系满足：

$v_i=\frac{\frac{x_i}{{\mathrm{\ρ}}_i}}{\underset{i=1}{\overset{m}{\mathrm{\Σ}}}\frac{x_i}{{\mathrm{\ρ}}_i}}---\left(7\right)$

其中，ρ_i为组分油i的密度。

作为优选的实施方式，步骤c中通过实验室数据管理系统或设置采样点或二者的结合来获取各掺炼线中组分油的当前性质。

有益效果：

本发明提出了一种高性能汽油调合非线性优化方法，针对辛烷值、雷德蒸汽压等性质，采用非线性的调合规则，并针对辛烷值提供多种调合规则供企业选择。该发明支持汽油性质含量在范围内定值、向下限优化、向上限优化及范围内优化四种优化方式，提高了生产的灵活性。以汽油性质含量目标偏差最小且经济最优为目标进行优化调合，既符合更为严格的质量排放标准，又满足生产效益和利益最大化的需求。

附图说明

图1是汽油调合系统总体框图。

图2是汽油调合优化子系统工作流程图。

具体实施方式

下面结合附图和实施案例对本发明作进一步的说明。

下面以本发明在某炼化企业实际实施情况并结合一具体算例，给出详细的计算过程和具体的操作流程。该企业有5路调合组分，其中第1路精致汽油(1#掺炼线)和第2路混合芳烃(2#掺炼线)为直调组分即直接从生产装置出来输入掺炼线参与调合，另外3路为别为己烷轻石(3#掺炼线)、重整抽余油(4#掺炼线)和MTBE(5#掺炼线)为非直调组分即从调合组分罐出来输入掺炼线参与调合。调合按批次进行，每批次调合一定的总量，本批次调合完成后再进行下一批次调合。本实施例在以本发明技术方案为前提下进行实施，但本发明的保护范围不限于下述的实施例。

如图1所示，本汽油调合系统包括以下部分：汽油调合优化计算机1、调合控制器2、实验室数据管理系统5和现场的设备及仪表。设备及仪表具体包括近红外在线性质检测仪3、硫含量在线检测仪4、调合组分罐6、产品罐7、调节阀8、流量计9、调合泵10和静态混合器11。调合控制器2接受来自汽油调合优化计算机1的参与调合各组分油流量占比指令，控制现场的设备及仪表完成流量控制。

如图2所示，本方法的工作流程主要包括以下步骤：

步骤a：本批次调合参数初始化。

结合一具体算例，假设本批次调合总量m_T＝1000t，调合总流速上限FT_max＝200t，调合总流速下限FT_min＝50t/h。

掺炼管线流速上限设置。假设1#掺炼线最大流速F_max，1＝57.5t/h，2#掺炼线最大流速F_max，2＝7t/h，3#掺炼线最大流速F_max，3＝20t/h，4#掺炼线最大流速F_max，4＝20t/h，5#掺炼线最大流速F_max，5＝22t/h，则F_max＝[57.5 7 20 20 22]；1#掺炼线最小流速F_min，1＝57.5t/h，2#掺炼线最小流速F_min，2＝7t/h，3#掺炼线最小流速F_min，3＝17t/h，4#掺炼线最小流速F_min，4＝3t/h，5#掺炼线最小流速F_min，5＝5t/h，则F_min＝[57.5 7 17 3 5]。

直调组分为第1路精致汽油和第2路混合芳烃，这两路组分油无调合组分罐6限制，直接从生产装置中获得组分油，故可参与的调合量不受限制，因此可给一个较大的数值，在本发明中默认等于总调合量m₁＝1000t，m₂＝1000t；本实施例中，通过采集其余三个组分的调合组分罐6的数据，获取各组分油的调和组分罐6可参与的调合量，采集得出：ZF01的可调量m₃＝1000t；ZF02的可调量m₄＝800t；ZFO3的可调量m₅＝1500t；则M＝[1000 1000 1000 800 1500]。

该优化系统允许人工根据经验预先设定各组分油的质量占比下限和上限。本实施例中，1#掺炼线所属组分油的质量占比下限为x_1min＝40％，上限为x_1max＝100％；2#掺炼线所属组分油的质量占比下限为x_2min＝10％，上限为x_2max＝40％；3#掺炼线所属组分油的质量占比下限为x_3min＝5％，上限为x_3max＝30％；4#掺炼线所属组分油的质量占比下限为x_4min＝0％，上限为x_4max＝20％；5#掺炼线所属组分油的质量占比下限为x_5min＝0％，上限为x_5max＝20％。从而：

X_min＝[40％ 10％ 5％ 0％ 0％]，X_max＝[100％ 40％ 30％ 20％ 20％]

步骤b、选择汽油参与优化的性质种类、并设定各性质的优化目标值、优化方式和权重。

本实施例中，汽油参与优化的性质种类选为研究法辛烷值(RON)、硫含量、雷德蒸汽压、苯含量，优化目标值和优化方向分别为：RON＝[92.292.5]，向下限优化，权重＝0.0001；硫含量/ppm＝[010]，范围内优化，权重＝0.1；雷德蒸汽压/KPa＝[4587]，范围内优化，权重＝1；苯量/v％＝[01]，范围内优化，权重＝1。从而：

rangeL＝[rangel₁ rangel₂ rangel₃ rangel₄]＝[92.2 0 45 0]

rangeU＝[rangeu₁ rangeu₂ rangeu₃ rangeu₄]＝[92.5 10 87 1]

goal₁＝92.5，goal₂∈[0，10]，goal₃∈[45，87]，goal₄∈[0，1]，即只有RON有明确的优化目标值。

Λ＝[λ₁ λ₂ λ₃ λ₄]＝[0.0001 0.1 1 1]。

优化时，考虑价格因素，本实施例中，各组分油的价格如下：

Price＝[price₁ price₂ price₃ price₄ price₅]＝[5602 4644 4955 6050 6401]，价格权重选择λ_Price＝1，价格优化选为向下限优化，故Price_min＝[4644]。

各性质采用的调合规则为：研究法辛烷值采用相互作用法；硫含量采用线性质量调合规则；雷德蒸汽压采用Chevron法；苯含量采用线性体积调合规则。

其中，研究法辛烷值各组分间的相互作用系数如表1所示：

表1各组分间的相互作用系数

组分油	精致汽油	混合芳烃	己烷轻石	重整抽余油	MTBE
						精致汽油	0	-2.67	20.5	0.13	1.36
混合芳烃	-2.67	0	-0.38	116.7	28.9
						己烷轻石	20.5	-0.38	0	40.29	3.25
重整抽余油	0.13	116.7	40.29	0	-49.26
						MTBE	1.36	28.9	3.25	-49.26	0

步骤c：获取掺炼线组分油的当前性质。

由于重整抽余油、MTBE性质相对稳定，因此在4#掺炼线和5#掺炼线上不设置近红外采样点；而精制汽油、混合芳烃和己烷轻石性质不稳定，因此在1#掺炼线、2#掺炼线上和3#掺炼线上均设置近红外采样点，通过近红外在线性质检测仪检测组分油各性质；同时精制汽油的硫含量波动大，需实时监控其实际硫含量，故在精制汽油1#掺炼线上设置硫含量采样点。

各掺炼线组分油性质值不在近红外在线性质检测仪3和硫含量在线检测仪4处获取的，均通过实验室数据管理系统5获取，本实施例中，当前掺炼汽油性质数据本实施例中，采集的具体数据如表2所示：

表2当前掺炼组分油性质数据

性质	1#	2#	3#	4#	5#
						RON	92	106	81.2	65	118
硫含量/wt％	2	1	1.7	0	0
						雷德蒸汽压/KPa	60	47	75.2	4.9	55
苯含量/wt％	1	0	0.4	0	0
						密度/(kg/m3)	732	869.3	656.8	870	740

则

$A=\left[\begin{array}{ccccc}{a}_{\mathrm{1,1}}& a_{\mathrm{1,2}}& a_{\mathrm{1,3}}& a_{\mathrm{1,4}}& a_{\mathrm{1,5}}\\ a_{\mathrm{2,1}}& a_{\mathrm{2,2}}& a_{\mathrm{2,3}}& a_{\mathrm{2,4}}& a_{\mathrm{2,5}}\\ a_{\mathrm{3,1}}& a_{\mathrm{3,2}}& a_{\mathrm{3,3}}& a_{\mathrm{3,4}}& a_{\mathrm{3,5}}\\ a_{\mathrm{4,1}}& a_{\mathrm{4,2}}& a_{\mathrm{4,3}}& a_{\mathrm{4,4}}& a_{\mathrm{4,5}}\end{array}\right]=\left[\begin{array}{ccccc}92& 106& 81.2& 65& 118\\ 2& 1& 1.7& 0& 0\\ 60& 47& 75.2& 4.9& 55\\ 1& 0& 0.4& 0& 0\end{array}\right]$

步骤d：计算掺炼线组分油配方。

利用式(1)自动计算五条掺炼原油质量调合占比X。式(1)的求解计算方法很多，本实施例中通过调用MATLAB优化求解函数fgoalattain，得到调合配方：x₁＝53.91％,x₂＝6.56％,x₃＝18.75％，x₄＝16.09％，x₅＝4.69％，即X＝[53.91％6.56％18.75％16.09％4.69％]。

根据以上优化出的配方计算得到调合后汽油各性质值分别为：研究法辛烷值＝92.2，硫含量/wt％＝1.46，雷德蒸汽压/KPa＝56.04，苯含量/wt％＝0.63，满足调合优化目标要求。各掺炼线调合流速分别为：1#掺炼线流速＝57.5t/h，2#掺炼线流速＝7t/h，3#掺炼线流速＝20t/h，4#掺炼线流速＝17.16t/h，#掺炼线流速＝5t/h，调合总流速＝106.66t/h，满足设备约束条件。

步骤e：占比下发到控制系统执行。

掺炼汽油配方经通讯接口下发到调合控制系统执行。

步骤f：判断本次调合是否完成。

如果本批次调合完成，则结束。如果未完成，则返回步骤c。

本发明未涉及方法均与现有技术相同或可采用现有技术加以实现。

Claims (7)

1.一种高性能汽油调合非线性优化方法，通过求解目标函数来获得各组分油占比，其特征在于该目标函数既考虑汽油参与优化的不同性质的调合值在优化范围内与优化目标值偏差最小，同时考虑生产成本最小化，具体包括以下步骤：

a、进行本次调合参数初始化：设定调合总量，调合总流速的上限和下限，各掺炼线流速的上限和下限，参与调合的组分油种类，各组分油量的上限，各组分油占比的上限和下限；

b、选择汽油参与优化的性质种类、调合规则、各性质的优化范围、优化方式、优化目标值和权重，组分油的价格目标和权重，不同性质种类选择不同的调合规则；

c、获取各掺炼线中组分油的当前性质；

d、通过目标函数计算各组分油占比；

e、占比下发到控制系统执行；

f、判断本次调合是否完成：完成则结束，未完成则返回步骤c。

2.根据权利要求1所述的一种高性能汽油调合非线性优化方法，其特征在于所述目标函数为下式：

$\min \mathrm{\Φ}\left(X\right)=\underset{i=1}{\overset{k}{\mathrm{\Σ}}}{\left[{\mathrm{\λ}}_k(f{\left(X\right)}_k-{\mathrm{goal}}_k)\right]}^2+{\left[{\mathrm{\λ}}_{\mathrm{Price}}(\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{price}}_i{x}_i-{\mathrm{Price}}_{\min })\right]}^2$

s.t.∑X＝1

X_min≤X≤X_max

FT＝F₁/x₁＝F₂/x₂＝…＝F_m/x_m     (1)

FT_min≤FT≤FT_max

F_min≤FT·X≤F_max

t·FT·X≤M

rangeL_k≤f(X)_k≤rangeU_k

式(1)中，Φ(X)为目标函数；k表示汽油参与优化的性质种类个数，f(X)_k为汽油参与优化的第k种性质所采用的调合规则，不同种性质采用不同的调合规则；m表示参与调合的掺炼线个数，即组分油个数；X为各组分油的质量占比，X＝[x₁,x₂…x_m]，X≥0；X_min和X_max分别为各组分油占比的上限和下限，X_min＝[x_min,1,x_min,2…x_min,m]，X_max≥0，X_max＝[x_max,1,x_max,2…x_max,m]；λ_k为汽油参与优化的第k种性质的权重，λ_k≥0；goal_k为汽油参与优化的第k种性质的优化目标值，goal_k≥0；Price为各组分油的价格，Price≥0，Price＝[price₁,price₂,…price_m]≥0；Price_min为各组分油的价格最低值，λ_Price为价格权重，λ_Price≥0；FT为当前调合总流速，FT≥0，FT_max为调合总流速的上限，FT_min为调合总流速的下限；F_max为各掺炼线流速的上限，F_max≥0，F_max＝[F_max,1,F_max,2…F_max,m]，F_min为各掺炼线流速的下限，F_min≥0，F_min＝[F_min,1,F_min,2…F_min,m]；t为本次调合预计执行时间，t≥0；M为各组分油可供调合总量，即各组分油量的上限，M＝[m₁,m₂…m_m]，M≥0；F为各组分油的流速，F＝[F₁,F₂…F_m]；rangeL_k为汽油参与优化的第k种性质的目标下限，rangeU_k为汽油参与优化的第k种性质的目标上限。

3.根据权利要求2所述的一种高性能汽油调合非线性优化方法，其特征在于优化方式选择：

向下限优化，goal_k＝rangeL_k；

向上限优化，goal_k＝rangeU_k；

范围内定值优化，goal_k＝rangeF_k，rangeF_k为设定的汽油参与优化的第k种性质的优化目标值，rangeF_k∈(rangeL_k,rangeU_k)；

范围内无目标优化，需满足f(X)_k∈[rangeL_k,rangeU_k]，此时f(X)_k-goal_k≡0。

4.根据权利要求2所述的一种高性能汽油调合非线性优化方法，其特征在于汽油参与优化的性质在硫含量、氧含量、硫醇硫含量、烯烃含量、芳烃含量、苯含量、研究法辛烷值、马达法辛烷值、抗暴指数和雷德蒸汽压中选择一到数种，采用不同的调合规则：

硫含量、氧含量和硫醇硫含量，采用线性质量调合规则如式

$a=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{a}_i{x}_i---\left(2\right)$

式(2)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，x_i为组分油i的质量占比；

烯烃含量、芳烃含量和苯含量采用线性体积调合规则如式

$a=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{a}_i{v}_i---\left(3\right)$

式(3)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，v_i为组分油i的体积占比；

研究法辛烷值、马达法辛烷值和抗暴指数采用相互作用法调合规则如式

$a=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{a}_i{v}_i+\underset{i=1}{\overset{m-1}{\mathrm{\Σ}}}{v}_i\underset{j=i+1}{\overset{m}{\mathrm{\Σ}}}{c}_{i,j}{v}_j---\left(4\right)$

式(4)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，c_i,j为两种组分油之间的调合效应系数，v_i为组分油i的体积占比；

研究法辛烷值、马达法辛烷值和抗暴指数或者采用Stewart法调合规则如式

$a=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}\left(v_i{D}_i(a_i+0.13P_i)\right)/\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{v}_i{D}_i---\left(5\right)$

式(5)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，P_i为调合汽油中组分油i体积含量与不饱和烃的体积含量差，D_i为根据P_i对应的权重指数，v_i为组分油i的体积占比；其中，不饱和烃为烯烃和芳烃的集合；

雷德蒸汽压采用Chevron法调合规则如式

$a={\left(\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{a}_i^{1.25}{v}_i\right)}^{0.8}---\left(6\right)$

式(6)中，a为汽油相应性质的调合值，a_i为组分油i的相应性质的参数，m为参与调合的组分油的个数，v_i为组分油i的体积占比；

式(3)-(6)中，v_i与x_i的关系满足：

$v_i=\frac{\frac{x_i}{{\mathrm{\ρ}}_i}}{\underset{i=1}{\overset{m}{\mathrm{\Σ}}}\frac{x_i}{{\mathrm{\ρ}}_i}}---\left(7\right)$

其中，ρ_i为组分油i的密度。

5.根据权利要求1所述的一种高性能汽油调合非线性优化方法，其特征在于步骤c中通过实验室数据管理系统来获取各掺炼线中组分油的当前性质。

6.根据权利要求1所述的一种高性能汽油调合非线性优化方法，其特征在于步骤c中设置采样点来获取各掺炼线中组分油的当前性质。

7.根据权利要求1所述的一种高性能汽油调合非线性优化方法，其特征在于步骤c中结合实验室数据管理系统和设置采样点来获取各掺炼线中组分油的当前性质。