Method and apparatus for evaluating catalytic cracking reaction product yield

Abstract

The application relates to a method and device for evaluating the yield of catalytic cracking reaction products, which belongs to the field of catalytic cracking, and comprises the following steps: firstly, performing lumping division; secondly, establishing a catalytic cracking reaction network; thirdly, establishing a catalytic cracking kinetics model according to the established reaction network; and finally, predicting the content of each raw material lumping through raw material properties, and substituting the content into the catalytic cracking kinetics model to predict the yield of each product lumping. Through optimization of lumping division, the raw material oil of catalytic cracking is divided into three layers of residual oil layer, wax oil layer and diesel oil layer according to the distillation range, and is finely divided into 11 raw material lumpings, and the products are divided into 6 lumpings, so that the catalytic cracking reaction law is fully reflected, the yield of each product lumping is predicted by relying on the catalytic cracking kinetics model, the prediction accuracy is greatly improved, and the product distribution and product properties can be accurately predicted for different catalytic cracking reaction processes.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of catalytic cracking, in particular to a method and device for evaluating the yield of catalytic cracking reaction products. BACKGROUND

[0002] Catalytic cracking is one of the important heavy oil lightening processes, and plays a very important role in the production of light oil products such as gasoline and diesel. The catalytic cracking reaction process model can be used to guide the design of the catalytic cracking device and optimize the operation of the device, and has an important role in improving the economic benefit of the device. The catalytic cracking reaction process is complex, and there are a large number of parallel sequential reactions. When modeling the reaction process, the lumped kinetic modeling method is generally used.

[0003] An ideal lumped kinetic model generally has the following characteristics: (1) applicable to different types of feedstocks; (2) high prediction accuracy for product distribution and product properties; (3) low requirements for analysis means, most refineries have analysis means. The first requirement is that the model divides the feedstock into lumped components, which can fully reflect the properties of the feedstock such as density, distillation range, viscosity, carbon residue, molecular weight, elemental composition, structural group composition, and metal content. The second requirement is that the model divides the products into lumped components, which can not only predict the product distribution, but also predict the product properties. The third requirement is that the model cannot be too complex, i.e. the number of lumps cannot be too large, and it can fully reflect the catalytic cracking reaction rules.

[0004] In the prior art, the following several representative lumped kinetic models are available:

[0005] 1) Chen Yang, Liu Jinze, Yu Zhaoshen, et al. Industrial simulation of 12-lump kinetic model for heavy oil catalytic cracking, published in East China Journal of Technology (Natural Science Edition), 2021, No. 47(3), pp. 262-271, which divides heavy oil catalytic cracking into 12 lumps and then performs kinetic modeling;

[0006] 2) Sun Shiyuan, Meng Fandong, Yan Hongfei, et al. Research on 14-lump kinetic model for heavy oil catalytic cracking, published in Henan Chemical Industry, 2017, No. 34(7), pp. 29-37, which divides heavy oil catalytic cracking into 14 lumps and then performs kinetic modeling;

[0007] 3) Hou Kaijun, Ouyang Fusheng, Gao Yongfu, et al. Research on 9-lump kinetic model for heavy oil catalytic cracking for producing more propylene, published in Petrochemical Technology and Application, 2017, No. 35(5), pp. 339-343, which only divides into 9 lumps;

[0008] Because of the difference in lumping division, the above lumping kinetic models are difficult to completely simulate all reactions in the catalytic cracking process, thereby leading to inaccurate prediction of product yield. SUMMARY

[0009] In order to solve the problem in the prior art that the established lumping kinetic model is not accurate enough in predicting the product yield of catalytic cracking reaction due to the inaccuracy of lumping division, the present application provides an evaluation method and device for the product yield of catalytic cracking reaction, which first divides the feedstock oil for catalytic cracking into three layers of residual oil layer, wax oil layer and diesel oil layer according to the distillation range, and finely divides them into 11 feedstock lumps, and at the same time divides the products into 6 lumps, then establishes a reaction network and a catalytic cracking kinetic model according to the divided lumps, and then predicts the content of each feedstock lump through the properties of the feedstock, and relies on the catalytic cracking kinetic model to predict the yield of each product lump, thereby greatly improving the accuracy of the prediction.

[0010] The technical scheme adopted by the present application to solve the above technical problem is as follows: an evaluation method for the product yield of catalytic cracking reaction, comprising the following steps:

[0011] (1) according to the principle of lumping kinetics and the actual needs of industry, the feedstock and products are divided into lumps;

[0012] (2) according to the divided lumps and the mechanism of catalytic cracking reaction, a catalytic cracking reaction network is established;

[0013] (3) according to the established reaction network, a catalytic cracking kinetic model is established;

[0014] (4) the content of each feedstock lump is predicted through the properties of the feedstock, and is substituted into the catalytic cracking kinetic model of step (3) to predict the yield of each product lump;

[0015] The lumping division of step (1) first divides the feedstock oil for catalytic cracking into three layers of residual oil layer, wax oil layer and diesel oil layer according to the distillation range, and finely divides them into 11 feedstock lumps, which are respectively the paraffin lump of residual oil layer, the naphthene lump of residual oil layer, the aromatic lump of residual oil layer, the aromatic lump in gum and asphalt, the paraffin lump of wax oil layer, the naphthene lump of wax oil layer, the aromatic lump of wax oil layer, the paraffin lump of diesel, the naphthene lump of diesel, the single-ring aromatic lump of diesel layer and the double-ring aromatic lump above the diesel layer;

[0016] The reaction products of catalytic cracking are divided into 6 lumps, which are respectively the paraffin lump of gasoline, the olefin lump of gasoline, the naphthene lump of gasoline, the aromatic lump of gasoline, the dry gas lump and the coke lump.

[0017] As an optimization of the above-mentioned method for evaluating the yield of catalytic cracking reaction products, in step (4), the specific operation for predicting the collective content of each feedstock by feedstock property is as follows:

[0018] 4.1) Obtain the distillation range distribution of the feedstock oil by using the existing analysis method, and calculate the mass content of the residue oil layer, the wax oil layer and the diesel oil layer respectively by using the following formula:

[0019]

[0020]

[0021] Y H =100-Y M -Y L (9)

[0022] In the formula, X represents the distillation range of the residue oil layer, the wax oil layer or the diesel oil layer, and the distillation range of the residue oil layer, the wax oil layer or the diesel oil layer is set with a sampling point every 10% from the initial boiling point to the final boiling point, together with the initial boiling point and the final boiling point of the distillation range, a total of 11 sampling points, and i and j represent the i-th and j-th points respectively;

[0023] 4.2) Calculate the average number of carbon atoms and hydrogen atoms in the molecules of the feedstock oil respectively

[0024] The average number of C atoms in the molecules of the feedstock oil is:

[0025]

[0026] In the formula, MW represents the relative molecular mass of one molecule in the feedstock oil, with the unit of kg / kmol, y C represents the mass content of carbon in the feedstock oil, with the unit of w%, and MW and y C are obtained by instrument analysis;

[0027] The average number of H atoms in the molecules of the feedstock oil is:

[0028]

[0029] In the formula, MW represents the relative molecular mass of one molecule in the feedstock oil, with the unit of kg / kmol, y H represents the mass content of hydrogen in the feedstock oil, with the unit of w%, and MW and y H are obtained by instrument analysis;

[0030] 4.3) Calculate the naphthenic carbon content and aromatic carbon content in the feedstock oil respectively

[0031] The naphthenic carbon content C N and the aromatic carbon content C A in the feedstock oil, with the unit of w%, are calculated by the following specific formula respectively:

[0032]

[0033]

[0034] wherein, denotes the density of the feed oil at 20°C in kg / m 3 , γ1 100 denotes the viscosity of the feed oil at 100°C in mm 2 / s, both parameters being obtained by instrumental analysis;

[0035] 4.4) The mass content of each feed set is calculated separately

[0036] In the residue layer, the mass content of the residue layer paraffin set, the residue layer naphthene set, the residue layer aromatic set and the aromatic set in the gum asphaltenes is respectively:

[0037]

[0038]

[0039]

[0040] Y HP = Y H - Y FA - Y HN - Y HA (17)

[0041] In the gas oil layer, the mass content of the gas oil layer paraffin set, the gas oil layer naphthene set and the gas oil layer aromatic set is respectively:

[0042]

[0043]

[0044] Y MP = Y M - Y MN - Y MA (20)

[0045] In the diesel layer, the mass content of the diesel paraffin set, the diesel naphthene set, the diesel layer monocyclic aromatic set and the diesel layer di- and polycyclic aromatic set is respectively:

[0046]

[0047]

[0048]

[0049] Y LP = Y L - Y LN - Y LA - Y LHA (24)

[0050] In the above formulae, CCR is the mass content of residual carbon in the feedstock oil, and this parameter is obtained by instrumental analysis;

[0051] 4.5) The mass content of each feedstock aggregate calculated in step 4.4) is substituted into the catalytic cracking kinetic model of step (3) to predict the aggregate yield of each product.

[0052] As another optimization scheme of the above-mentioned evaluation method of the product yield of catalytic cracking reaction, the specific operation of predicting the aggregate yield of each product in step 4.5) is as follows:

[0053] The initial content of each feedstock aggregate calculated according to formula (14) to formula (24) is then used to calculate the aggregate yield of each product by formula (5);

[0054]

[0055] In the formulae:

[0056] Y = [y HP , y HN , y HA , y FA , y MP , y MN , y MA , y LP , y LN , y LA , y LHA , y GP , y GO , y GN , y GA , y GAS , y CK ] T is the mass fraction vector of each aggregate component;

[0057] P represents the system pressure, R is the gas constant, 8.314 J / (mol·K); T represents the system temperature, X = x / H represents the dimensionless relative distance at the x cross section in the bed, H represents the total length of the catalyst bed, S WH represents the real weight hourly space velocity, K is the reaction rate constant matrix;

[0058] y i represents the mass content of the i-th aggregate, M iMw(i) represents the average relative molecular weight of the ith group; and a(i) represents the concentration of the ith group i The calculation formula is:

[0059]

[0060] In formula (4), a i The unit is moles / g.

[0061] The K is a reaction rate constant matrix, which is composed of reaction rate constants k j of all reactions in the catalytic cracking reaction network.

[0062]

[0063] In formula (6), A j represents the pre-exponential factor of the jth reaction in the catalytic cracking reaction network, E j represents the activation energy of the jth reaction, k j is the reaction rate constant of the jth reaction in the constructed catalytic cracking reaction network.

[0064] As another optimization scheme of the above-mentioned evaluation method of the catalytic cracking reaction product yield, the calculation method of the content of each component gas in the dry gas group is as follows: according to the structural group composition of the raw oil and the reaction temperature, reaction pressure, reaction time and catalyst-oil ratio of the catalytic cracking reaction, the yields of hydrogen, methane, ethane, ethylene, propane, propylene, n-butane, isobutane, isobutene, n-butene, cis-butene and trans-butene are calculated.

[0065] Reaction correlation formula:

[0066] y i = (C(i, 1) · T + C(i, 2) · P + C(i, 3) · t + C(i, 4) · CO)

[0067] · (y CPH+CPM · C(i, 5) + y CNH+CNM · C(i, 6) + y CAH+CAM+HA · C(i, 7)) / 100 (25)

[0068] In the formula

[0069] C(11,7) = [0.04, 0.10, 0.09, 0.37, 0.11, 0.02, 0.10, 0.01, 0.60, 0.13, 0.02; 0.06, 0.38, 0.11, 0.36, 0.01, 0.31, 0.37, -0.06, 0.13, 0.05, -0.96; -0.43, 0.08, -0.23, -0.02, 0.42, 0.20, 1.00, 0.30, 0.01, -0.15, -0.49; 0.01, -0.39, 0.10, -0.86, -0.03, -0.23, -0.12, -0.09, 0.66, 0.03, 0.23; 0.15, 0.02, -0.23, -0.48, 0.20, -0.30, 0.03, -0.33, -0.14, 0.08, -0.39; 0.87, 0.20, 0.47, 0.15, 0.17, -0.13, -0.27, -0.44, -0.40, -0.13, -0.27; -0.29, 0.22, -0.09, 0.12, 0.56, 0.34, 1.14, 0.44, 0.15, -0.01, -0.35];

[0070] T represents the reaction temperature in °C; P represents the reaction pressure in MPa; t represents the reaction time in s; CO represents the catalyst to oil ratio in kg catalyst / kg feedstock; y i represents the yield of the i-th gas in w%.

[0071] As another optimization of the method for evaluating the yield of the catalytic cracking reaction product, the catalytic cracking reaction network established in the step (2) is as follows:

[0072] The residue layer paraffin lump is respectively reacted with the wax oil layer paraffin lump, the diesel paraffin lump, the gasoline paraffin lump, the gasoline aromatic lump, the gasoline olefin lump, the gasoline naphthene lump, the dry gas lump and the coke lump;

[0073] The residue layer naphthene lump is respectively reacted with the diesel naphthene lump, the wax oil layer naphthene lump, the gasoline paraffin lump, the gasoline olefin lump, the dry gas lump, the gasoline naphthene lump, the gasoline aromatic lump and the coke lump;

[0074] The residue layer aromatic lump is respectively reacted with the gasoline paraffin lump, the gasoline olefin lump, the wax oil layer aromatic lump, the diesel layer monocyclic aromatic lump, the gasoline naphthene lump, the coke lump, the gasoline aromatic lump and the diesel layer polycyclic aromatic lump;

[0075] The aromatic lump in the gum asphalt is respectively reacted with the wax oil layer aromatic lump, the diesel layer polycyclic aromatic lump and the coke lump;

[0076] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0077] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0078] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0079] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0080] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0081] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0082] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0083] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0084] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0085] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0086] The wax oil layer paraffin aggregate is respectively reacted with a gasoline paraffin aggregate, a dry gas aggregate, a gasoline olefin aggregate, a gasoline aromatic aggregate, a gasoline naphthene aggregate, and a coke aggregate;

[0087] As another optimization of the above-mentioned evaluation method of the catalytic cracking reaction product yield, the basic equation for establishing the catalytic cracking kinetics model in step (3) is:

[0088]

[0089] In the formula, a iLet represent the concentration of the i-th ensemble, in molesi / g; P represent the system pressure; R is the gas constant, 8.314 J / (mol·K); T represents the system temperature; X = x / H represents the dimensionless relative distance at section x in the bed; H represents the total length of the catalyst bed; S WH Represents the actual weight-time space velocity, K is the reaction rate constant matrix, and a = [a1, ..., a2]. ni ] T This is the lumped component concentration vector.

[0090] As another optimization scheme for the evaluation method of the above-mentioned catalytic cracking reaction product yield, the basic equation of the catalytic cracking kinetic model is derived from the continuity equation and the reaction rate equation:

[0091] The continuity equation is:

[0092]

[0093] The reaction rate equation is:

[0094]

[0095] In the formula, subscript i represents the lumped component, subscript j represents the j-th reaction, and ρ represents the density of the oil-gas mixture in g / cm³. 3 , ai represents the concentration of the i-th set, in molesi / g, t represents the reaction time, and G v This represents the mass flow rate of oil and gas across a cross section, in g / (cm). 2 ·h), where x represents the distance from the riser inlet into the reactor, R i Let n represent the reaction rate of the i-th ensemble. r v represents the number of reactions. i,j R represents the stoichiometric coefficient of i lumped in reaction j. j k represents the reaction rate of reaction j. j ρ represents the rate constant of reaction j. c Catalyst density relative to reactor volume, in g / cm³ 3 ε represents the porosity, P represents the system pressure, R is the gas constant, and T represents the system temperature. This indicates the effect of catalyst coking on activity, and assumes that the catalyst coking rate is only a function of the catalyst residence time, t. c η is the catalyst residence time, and β is the catalyst deactivation constant.

[0096] A model-building apparatus for evaluating the yield of catalytic cracking reaction products, the apparatus comprising a modeling module for modeling using the aforementioned evaluation method for the yield of catalytic cracking reaction products.

[0097] An electronic device includes a processor and a memory, the memory storing computer-readable instructions that, when executed by the processor, enable the performance evaluation method described above.

[0098] A readable storage medium having a computer program stored thereon, which, when executed by a processor, performs the steps of the evaluation method described above.

[0099] Step (1) of this invention is the lumped component division, which divides the catalytic cracking reaction system into several lumped components according to the kinetic characteristics of each hydrocarbon molecule. In the kinetic study, each lumped component is considered as a virtual single molecular component. This invention divides the feedstock into three layers according to the distillation range: diesel layer, wax oil layer, and residue oil layer. Each layer is divided into alkyl carbon, cycloalkane carbon, and aromatic carbon according to the structural group composition. At the same time, considering that the aromatic carbon in gum and asphaltenes is mostly fused ring aromatic carbon, and its cracking performance is different from that of the aromatic carbon in the residue oil layer, it is treated as a separate lumped component. In addition, considering that the aromatic carbon of monocyclic aromatic hydrocarbons and the aromatic carbon of bicyclic or higher aromatic hydrocarbons in the diesel layer have significantly different cracking performance, they are treated as two separate lumped components. Thus, the feedstock can be divided into 11 lumped components. Among the products, gasoline is divided into alkanes, alkenes, cycloalkanes, and aromatics according to the PONA composition, and the remaining products are divided into gases and coke. Thus, the products are divided into 6 lumped components.

[0100] The establishment of the catalytic cracking kinetic reaction network in step (2) generally follows these principles: the catalytic cracking reaction process is a parallel-sequential reaction, gasoline is both a reactant and a product, and gas and coke are the final products; the catalytic cracking reaction process is mainly based on the reaction of reorganization and splitting to generate light components; in the reactions of cracking the residue layer to generate the wax oil layer, the diesel layer, and the wax oil layer to generate the diesel layer, alkyl carbons can only crack to generate alkyl carbons, cycloalkane carbons can only crack to generate cycloalkane carbons, and aromatic carbons can only crack to generate aromatic carbons; for the residue layer, wax oil layer, and diesel layer, there is no interaction between the lumped layers in the same layer; the reaction of olefins and cycloalkanes in the gasoline layer to generate alkanes and aromatics through hydrogen transfer is considered. Based on this, the reaction network established in this invention includes 17 lumped layers totaling 77 reactions, including alkane cracking reaction, condensation reaction, cycloalkane cracking reaction, aromatization reaction, condensation reaction, hydrogen transfer reaction, olefin cracking reaction, condensation reaction, hydrogen transfer reaction, and aromatic cracking reaction and condensation reaction.

[0101] Compared with the prior art, the present invention has the following beneficial effects:

[0102] 1) This invention optimizes the lumped layer division, first dividing the feedstock oil of catalytic cracking into three layers according to the distillation range: residue oil layer, wax oil layer, and diesel oil layer, and further subdividing them into 11 feedstock lumped layers. At the same time, the products are divided into 6 lumped layers. Then, a reaction network and catalytic cracking kinetic model are established according to the divided lumped layers. The content of each feedstock lumped layer is predicted by the feedstock properties, and the yield of each product lumped layer is predicted by the catalytic cracking kinetic model, which greatly improves the accuracy of the prediction.

[0103] 2) The 17-lumped kinetic model proposed in this invention divides the feedstock into 11 lumped groups according to its distillation range and structural group composition. It can fully reflect the influence of feedstock density, distillation range, viscosity, carbon residue, molecular weight, elemental composition, structural group composition, metal content, and other properties, and has strong adaptability to different types of feedstocks. This invention can calculate the content of the 11 lumped groups of feedstocks based on their conventional properties, and has low requirements for analytical methods. The 17-lumped kinetic model proposed in this invention includes 77 reactions, covering all reaction types in catalytic cracking processes such as catalytic cracking, thermal cracking, condensation, and hydrogen transfer. It fully reflects the reaction law of catalytic cracking and can accurately predict product distribution and product properties for different catalytic cracking reaction processes. Attached Figure Description

[0104] Figure 1 This is a schematic diagram of the reaction network of the 17 lumped reaction kinetic models for catalytic cracking in this invention;

[0105] Figure labels: HP, aggregate of alkanes in residual oil layer; HN, aggregate of cycloalkanes in residual oil layer; HA, aggregate of aromatics in residual oil layer; FA, aggregate of aromatics in asphaltene; MP, aggregate of alkanes in wax oil layer; MN, aggregate of cycloalkanes in wax oil layer; MA, aggregate of aromatics in wax oil layer; LP, aggregate of alkanes in diesel fuel layer; LN, aggregate of cycloalkanes in diesel fuel layer; LA, monocyclic aromatics in diesel fuel layer; LHA, bicyclic or more aromatics in diesel fuel layer; GP, aggregate of alkanes in gasoline layer; GO, aggregate of olefins in gasoline layer; GN, aggregate of cycloalkanes in gasoline layer; GA, aggregate of aromatics in gasoline layer; GS, aggregate of dry gas layer; CK, aggregate of coke layer. Detailed Implementation

[0106] The technical solution of the present invention will be further described in detail below with reference to specific embodiments. Any parts not explicitly described in the following embodiments should be understood as prior art known or should be known by those skilled in the art.

[0107] Example 1

[0108] A method for evaluating the yield of catalytic cracking reaction products includes the following steps:

[0109] (1) According to the principle of lumped dynamics and the actual needs of industry, raw materials and products are lumped together;

[0110] The lumped layer division in this step is based on the distillation range. The feedstock oil for catalytic cracking is first divided into three layers: residue oil layer H, wax oil layer M, and diesel oil layer L. These layers are further subdivided into 11 feedstock lumps, namely: residue oil layer alkane lump HP, residue oil layer cycloalkane lump HN, residue oil layer aromatic lump HA, aromatic lump FA in asphaltenes, wax oil layer alkane lump MP, wax oil layer cycloalkane lump MN, wax oil layer aromatic lump MA, diesel oil alkane lump LP, diesel oil cycloalkane lump LN, diesel oil layer monocyclic aromatic lump LA, and diesel oil layer bicyclic or higher aromatic lump LHA.

[0111] The reaction products of catalytic cracking are divided into 6 aggregates: gasoline alkane aggregate GP, gasoline olefin aggregate GO, gasoline cycloalkane aggregate GN, gasoline aromatics aggregate GA, dry gas aggregate GS, and coke aggregate CK.

[0112] (2) Based on the lumped sum catalytic cracking reaction mechanism, establish a catalytic cracking reaction network;

[0113] like Figure 1 As shown, the residue oil layer alkane aggregate HP was reacted with the wax oil layer alkane aggregate MP, the diesel oil alkane aggregate LP, the gasoline alkane aggregate GP, the gasoline aromatics aggregate GA, the gasoline olefins aggregate GO, the gasoline cycloalkanes aggregate GN, the dry gas aggregate GS, and the coke aggregate CK, respectively.

[0114] The naphthenic aggregate HN in the residue oil layer was reacted with the naphthenic aggregate LN in the diesel oil layer, the naphthenic aggregate MN in the wax oil layer, the naphthenic aggregate GP in the gasoline layer, the olefin aggregate GO in the gasoline layer, the dry gas aggregate GS, the naphthenic aggregate GN in the gasoline layer, the aromatics aggregate GA in the gasoline layer, and the coke aggregate CK, respectively.

[0115] The aromatic hydrocarbon aggregate HA in the residual oil layer was reacted with the gasoline alkane aggregate GP, the gasoline olefin aggregate GO, the wax oil aromatic hydrocarbon aggregate MA, the diesel monocyclic aromatic hydrocarbon aggregate LA, the gasoline cycloalkane aggregate GN, the coke aggregate CK, the gasoline aromatic hydrocarbon aggregate GA, and the diesel bicyclic or higher aromatic hydrocarbon aggregate LHA, respectively.

[0116] The aromatic hydrocarbon aggregate FA in the asphaltene was reacted with the aromatic hydrocarbon aggregate MA in the wax oil layer, the bicyclic or higher aromatic hydrocarbon aggregate LHA in the diesel layer, and the coke aggregate CK, respectively.

[0117] The wax oil layer alkane aggregate MP was reacted with dry gas aggregate GS, diesel alkane aggregate LP, gasoline olefin aggregate GO, gasoline aromatics aggregate GA, gasoline cycloalkanes aggregate GN and coke aggregate CK, respectively.

[0118] The wax oil layer cycloalkanes MN was reacted with gasoline alkane aggregates GP, dry gas aggregates GS, gasoline olefin aggregates GO, diesel cycloalkanes LN, gasoline cycloalkanes GN, gasoline aromatics aggregates GA, and coke aggregates CK, respectively.

[0119] The aromatic hydrocarbon aggregate MA of the wax oil layer is reacted with the dry gas aggregate GS, the monocyclic aromatic hydrocarbon aggregate LA of the diesel layer, the bicyclic or higher aromatic hydrocarbon aggregate LHA of the diesel layer, the gasoline aromatic hydrocarbon aggregate GA, and the coke aggregate CK, respectively.

[0120] The diesel alkane aggregate LP was reacted with the dry gas aggregate GS, the gasoline alkane aggregate GP, the gasoline olefin aggregate GO, the gasoline cycloalkane aggregate GN, the gasoline aromatics aggregate GA, and the coke aggregate CK, respectively.

[0121] The diesel cycloalkane aggregate LN was reacted with the gasoline alkane aggregate GP, dry gas aggregate GS, gasoline olefin aggregate GO, gasoline cycloalkane aggregate GN, gasoline aromatics aggregate GA, and coke aggregate CK, respectively.

[0122] The diesel layer monocyclic aromatic hydrocarbon cluster LA was reacted with the dry gas cluster GS, the gasoline aromatic hydrocarbon cluster GA, and the coke cluster CK, respectively.

[0123] The diesel layer bicyclic or higher aromatic hydrocarbon cluster LHA was reacted with the dry gas cluster GS and the coke cluster CK, respectively.

[0124] The gasoline alkane aggregate GP was reacted with the dry gas aggregate GS and the coke aggregate CK, respectively.

[0125] The gasoline olefin aggregate GO was reacted with the gasoline alkane aggregate GP, the dry gas aggregate GS, and the coke aggregate CK, respectively.

[0126] The gasoline cycloalkane aggregate GN was reacted with the dry gas aggregate GS, the gasoline aromatics aggregate GA, and the coke aggregate CK, respectively.

[0127] The gasoline aromatic hydrocarbon aggregate GA was reacted with the dry gas aggregate GS and the coke aggregate CK, respectively.

[0128] (3) Based on the established reaction network, a catalytic cracking kinetic model was established;

[0129] The basic equations for establishing the catalytic cracking kinetic model in this step are:

[0130]

[0131] In the formula, a iLet represent the concentration of the i-th ensemble, in molesi / g; P represent the system pressure; R is the gas constant, 8.314 J / (mol·K); T represents the system temperature; X = x / H represents the dimensionless relative distance at section x in the bed; H represents the total length of the catalyst bed; S WH Represents the actual weight-time space velocity, K is the reaction rate constant matrix, and a = [a1, ..., a2]. ni ] T This represents the lumped component concentration vector;

[0132] The basic equations of the catalytic cracking kinetic model are derived from the continuity equation and the reaction rate equation:

[0133] The continuity equation is:

[0134]

[0135] The reaction rate equation is:

[0136]

[0137] In the formula, subscript i represents the lumped component, subscript j represents the j-th reaction, and ρ represents the density of the oil-gas mixture in g / cm³. 3 , ai represents the concentration of the i-th set, in molesi / g, t represents the reaction time, and G v This represents the mass flow rate of oil and gas across a cross section, in g / (cm). 2 ·h), where x represents the distance from the riser inlet into the reactor, R i Let n represent the reaction rate of the i-th ensemble. r v represents the number of reactions. i,j R represents the stoichiometric coefficient of i lumped in reaction j. j k represents the reaction rate of reaction j. j ρ represents the rate constant of reaction j. c Catalyst density relative to reactor volume, in g / cm³ 3 ε represents the porosity, P represents the system pressure, R is the gas constant, and T represents the system temperature. This indicates the effect of catalyst coking on activity, and assumes that the catalyst coking rate is only a function of the catalyst residence time, t. c Where β is the catalyst residence time, and β is the catalyst deactivation constant;

[0138] (4) Predict the total content of each raw material by the properties of the raw materials, and substitute it into the catalytic cracking kinetic model in step (3) to predict the yield of each product.

[0139] The specific steps are as follows:

[0140] 4.1) Using existing analytical methods, obtain the distillation range distribution of the feedstock oil, and calculate the mass content of the residue layer H, wax oil layer M, and diesel oil layer L using the following formula:

[0141]

[0142]

[0143] Y H =100-Y M -Y L (9)

[0144] In the formula, X represents the distillation range of the residue oil layer H, wax oil layer M, or diesel oil layer L. The distillation range of the residue oil layer H, wax oil layer M, or diesel oil layer L is set at 10% intervals from the initial boiling point to the final boiling point, with a total of 11 sampling points including the initial boiling point and the final boiling point of the distillation range. i and j represent the i-th and j-th points, respectively.

[0145] 4.2) Calculate the average number of carbon and hydrogen atoms in the crude oil molecules respectively.

[0146] The average number of carbon atoms in the crude oil molecule:

[0147]

[0148] In the formula, MW represents the relative molecular mass of one molecule in the feedstock, with units of kg / kmol.

[0149] y C This indicates the mass content of carbon in the crude oil, expressed in w%, MW, and y. C All were obtained through instrumental analysis;

[0150] The average number of H atoms in the crude oil molecule:

[0151]

[0152] In the formula, MW represents the relative molecular mass of one molecule in the feedstock, with units of kg / kmol.

[0153] y H This indicates the mass content of hydrogen in the crude oil, expressed in w%, MW, and y. H All were obtained through instrumental analysis;

[0154] 4.3) Calculate the naphthenic carbon content and aromatic carbon content in the feedstock oil respectively.

[0155] The naphthenic and aromatic carbon contents in the feedstock, expressed in w%, are calculated using the following formulas:

[0156]

[0157]

[0158] In the formula, This indicates the density of the crude oil at 20℃, expressed in kg / m³. 3 γ 100 This indicates the viscosity of the feedstock oil at 100℃, in mm. 2 / s, these two parameters are obtained through instrumental analysis;

[0159] 4.4) Calculate the mass content of each raw material aggregate separately.

[0160] In the residual oil layer H, the mass contents of the residual oil layer alkane aggregate HP, residual oil layer cycloalkane aggregate HN, residual oil layer aromatic hydrocarbon aggregate HA, and aromatic hydrocarbon aggregate FA in the asphaltene are respectively:

[0161]

[0162]

[0163]

[0164] Y HP =Y H -Y FA -Y HN -Y HA (17)

[0165] In wax oil layer M, the mass contents of the wax oil layer alkane aggregate MP, the wax oil layer cycloalkane aggregate MN, and the wax oil layer aromatic aggregate MA are respectively:

[0166]

[0167]

[0168] Y MP =Y M -Y MN -Y MA (20)

[0169] In diesel layer L, the mass contents of diesel alkane aggregate LP, diesel cycloalkane aggregate LN, diesel monocyclic aromatic hydrocarbon aggregate LA, and diesel bicyclic or higher aromatic hydrocarbon aggregate LHA are respectively:

[0170]

[0171]

[0172]

[0173] Y LP =Y L -Y LN -Y LA -Y LHA (twenty four)

[0174] In the above formulas, CCR is the mass content of residual carbon in the feedstock, and this parameter is obtained through instrumental analysis.

[0175] 4.5) Substitute the mass content of each raw material aggregate calculated in step 4.4) into the catalytic cracking kinetic model in step (3) to predict the yield of each product aggregate;

[0176] The specific steps for predicting the aggregate yield of each product in step 4.5 are as follows:

[0177] The initial content of each raw material aggregate is calculated according to formula (14)-formula (24), and then the yield of each product aggregate is calculated by formula (5).

[0178]

[0179] In the formula:

[0180] Y = [y HP ,y HN ,y HA ,y FA ,y MP ,y MN ,y MA ,y LP ,y LN ,y LA ,y LHA ,y GP ,y GO ,y GN ,y GA ,y GAS ,y CK ] T This is the mass fraction vector of each ensemble component;

[0181] P represents the system pressure, R is the gas constant, 8.314 J / (mol·K); T represents the system temperature, X = x / H represents the dimensionless relative distance at section x in the bed, H represents the total length of the catalyst bed, and S... WH Represents the actual gravity-time space velocity, and K is the reaction rate constant matrix;

[0182] y i M represents the total mass content of the i-th set. i The average relative molecular mass of the i-th set; the concentration a of the i-th set. i The calculation formula is:

[0183]

[0184] In equation (4), a i The unit is molesi / g;

[0185] K is the reaction rate constant matrix, which is composed of the reaction rate constants k of all reactions in the catalytic cracking reaction network. j constitute:

[0186]

[0187] In equation (6), A j E represents the pre-exponential factor of the j-th reaction in the catalytic cracking reaction network. j k represents the activation energy of the j-th reaction. j Let be the reaction rate constant of the j-th reaction in the constructed catalytic cracking reaction network.

[0188] Example 2

[0189] This embodiment is an improvement on Embodiment 1. The main method is the same as that of Embodiment 1. The improvement is that the dry gas aggregate GS is composed of: hydrogen, methane, ethane, ethylene, propane, propylene, n-butane, isobutane, isobutene, n-butene, cis-butene, and trans-butene. This embodiment was designed to calculate the content of each component gas in the dry gas aggregate GS.

[0190] Based on the structural composition of the feedstock and the reaction temperature, pressure, time and catalyst-to-oil ratio of the catalytic cracking reaction, the yields of hydrogen, methane, ethane, ethylene, propane, propylene, n-butane, isobutane, isobutene, n-butene, cis-butene and trans-butene are calculated using the following reaction correlation.

[0191] Reaction correlation:

[0192] y i =(C(i,1)·T+C(i,2)·P+C(i,3)·t+C(i,4)·CO)·(y CPH+CPM ·C(i,5)+y CNH+CNM ·C(i,6)+y CAH+CAM+HA ·C(i,7)) / 100(25)

[0193] In the formula

[0194] C(11,7)=[0.04,0.10,0.09,0.37,0.11,0.02,0.10,0.01,0.60,0.13,0.02;0.06,0.38,0.11,0.36,0.01,0.31,0.37,-0.06,0 .13,0.05,-0.96;-0.43,0.08,-0.23,-0.02,0.42,0.20,1.00,0.30,0.01,-0.15,-0.49;0.01,-0.39,0.10,-0.86,-0.03,-0. 23, -0.12, -0.09, 0.66, 0.03, 0.23; 0.15, 0.02, -0.23, -0.48, 0.20, -0.30, 0.03, -0.33, -0.14, 0.08, -0.39; 0.87, 0.20, 0.47, 0.15, 0.17, -0.13, -0.27, -0.44, -0.40, -0.13, -0.27; -0.29, 0.22, -0.09, 0.12, 0.56, 0.34, 1.14, 0.44, 0.15, -0.01, -0.35];

[0195] T represents the reaction temperature in °C; P represents the reaction pressure in MPa; t represents the reaction time in seconds; CO represents the catalyst-to-oil ratio in kg catalyst / kg feedstock; y i This represents the yield of the i-th gas, expressed in w%.

[0196] Example 3

[0197] A model-building apparatus for evaluating the yield of catalytic cracking reaction products, the apparatus comprising a modeling module for modeling using the aforementioned evaluation method for the yield of catalytic cracking reaction products.

[0198] Example 4

[0199] An electronic device includes a processor and a memory, the memory storing computer-readable instructions that, when executed by the processor, enable the performance evaluation method described above.

[0200] In this embodiment, the electronic device includes a modeling device, a processor, a memory, a storage controller, a peripheral interface, an input / output unit, an audio unit, and a display unit, etc.

[0201] Specifically, the memory, storage controller, processor, peripheral interface, input / output unit, audio unit, and display unit are electrically connected directly or indirectly to each other to achieve data transmission or interaction. The modeling device includes at least one software functional module that can be stored in the memory or embedded in the operating system (OS) of the modeling device in the form of software or firmware. The processor is used to execute the executable module stored in the memory, including the software functional module or computer program.

[0202] The memory can be, but is not limited to, random access memory (RAM), read-only memory (ROM), programmable read-only memory (PROM), erasable read-only memory (EPROM), electrically erasable read-only memory (EEPROM), etc. The memory stores programs, and after receiving execution instructions, the processor executes the corresponding program. The method executed by the server defined in the flow process of this application can be applied to the processor or implemented by the processor.

[0203] A processor can be an integrated circuit chip with signal processing capabilities. The processors mentioned above can be general-purpose processors, including central processing units (CPUs), network processors (NPs), etc.; they can also be digital signal processors (DSPs), application-specific integrated circuits (ASICs), off-the-shelf programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0204] Peripheral interfaces couple various input / output devices to the processor and memory. Peripheral interfaces, processors, and memory controllers can be implemented in a single chip, or they can be implemented separately by independent chips.

[0205] The input / output unit, audio unit, and display unit are all existing technologies. For example, the input / output unit is used to provide users with input data to enable interaction between the user and the server (or local terminal), and can be a mouse, keyboard, etc.; the audio unit provides users with an audio interface, which may include one or more microphones, one or more speakers, and audio circuitry; the display unit provides an interactive interface (e.g., a user operation interface) between the electronic device and the user or is used to display image data for the user's reference.

[0206] Example 5

[0207] A readable storage medium having a computer program stored thereon, which, when executed by a processor, performs the steps of the evaluation method described above.

[0208] To verify the effectiveness of the present invention, the following experimental examples were conducted:

[0209] Experimental Example

[0210] Following the method in Example 1, a lumped kinetic model for catalytic cracking was established, and the parameters in the model were then fitted based on experimental data. The properties of the experimental feedstock are shown in Table 1, and the product distribution is shown in Table 2. The comparison between the fitted values ​​calculated in Example 1 and the experimental values ​​is shown in Table 3. The model was calibrated based on industrial plant data to obtain the plant calibration factor. The comparison between the actual data from the industrial plant and the fitted data is shown in Table 4.

[0211] Comparative Example 1

[0212] Comparative Example 1 is based on the lumped kinetic model of catalytic cracking 12 proposed by Chen Yang, Liu Jinze, Yu Zhaochen et al. in 2021, which establishes the lumped kinetic model of catalytic cracking 12. The parameters in the model are fitted according to the same experimental data as in the experimental example. The comparison between the fitted values ​​and the experimental values ​​is shown in Table 3 below.

[0213] Comparative Example 2

[0214] Comparative Example 2 is based on the lumped kinetic model of catalytic cracking 14 proposed by Sun Shiyuan, Meng Fandong, Yan Hongfei and others in 2017. The lumped kinetic model of catalytic cracking 14 was established. The parameters in the model were fitted according to the same experimental data as in the experimental example. The comparison between the fitted values ​​and the experimental values ​​is shown in Table 3 below.

[0215] Comparative Example 3

[0216] Comparative Example 3 is based on the 9-lumped kinetic model of catalytic cracking proposed by Hou Kaijun, Ouyang Fusheng, Gao Yongfu and others in 2017, which established the 9-lumped kinetic model of catalytic cracking. The parameters in the model were fitted with the same experimental data as in the experimental example. The comparison between the fitted values ​​and the experimental values ​​is shown in Table 3 below.

[0217] As can be seen from the data comparison in Table 3, the relative error of the product distribution in the experimental example is significantly smaller than that in comparative examples 1, 2 and 3, indicating that the modeling method proposed in this invention has higher prediction accuracy.

[0218] Table 1 Properties of Crude Oil

[0219] Item Feed oil 1 Feed oil 2 Feed oil 3 Density, (20°C kg / m3 3) ]]> 924.5 925.5 914.6 viscosity, (100°C, mm 2 / s) 15.2 12.7 5.4 Residual carbon, w% 3.54 2.24 0.79 Molecular weight, kg / kmol 516 508 498 Hydrogen to carbon ratio 7.14 7.22 7.12 Saturated hydrocarbon content, w% 71.6 60.4 42.6 Aromatics + resin + asphaltenes content, w% 28.4 39.6 57.4 500°C content, w% 4.0 8.8 27.2

[0220] Table 2 Product Distribution

[0221]

[0222]

[0223] Table 3. Average relative errors of experimental and fitted values

[0224] Item Experimental example Comparative example 1 Comparative example 2 Comparative example 3 Product distribution / weight % Dry gas 2.36 3.68 4.16 4.85 Liquefied gas 2.07 4.18 3.88 4.13 Gasoline 0.13 1.09 1.66 1.93 Diesel 1.08 2.11 2.54 3.66 Slurry oil 1.94 4.81 4.89 5.44 Coke 1.77 3.66 3.87 4.32

[0225] Table 4. Average Relative Error of Actual Data and Fitted Values ​​for Industrial Equipment

[0226] Item Experimental example Product distribution / weight % Dry gas 0.03 Liquefied gas 0.03 Gasoline 0.01 Diesel 0.01 Slurry oil 0.02 Coke 0.02

Claims

1. A method for evaluating the yield of catalytic cracking reaction products, comprising the following steps: (1) According to the principle of lumped dynamics and the actual needs of industry, raw materials and products are lumped together; (2) Based on the lumped sum catalytic cracking reaction mechanism, establish a catalytic cracking reaction network; (3) Based on the established reaction network, a catalytic cracking kinetic model was established; (4) Predict the content of each feedstock aggregate by means of the properties of the feedstock, and substitute it into the catalytic cracking kinetic model in step (3) to predict the yield of each product aggregate; characterized in that: the aggregate division in step (1) is based on the distillation range, the feedstock oil of catalytic cracking is first divided into three layers: residue oil layer (H), wax oil layer (M) and diesel oil layer (L), and further subdivided into 11 feedstock aggregates, namely, residue oil layer alkane aggregate (HP), residue oil layer cycloalkane aggregate (HN), residue oil layer aromatics aggregate (HA), aromatics aggregate in asphaltene (FA), wax oil layer alkane aggregate (MP), wax oil layer cycloalkane aggregate (MN), wax oil layer aromatics aggregate (MA), diesel oil alkane aggregate (LP), diesel oil cycloalkane aggregate (LN), diesel oil layer monocyclic aromatics aggregate (LA) and diesel oil layer bicyclic and above aromatics aggregate (LHA); The reaction products of catalytic cracking are divided into six clusters: gasoline alkane cluster (GP), gasoline olefin cluster (GO), gasoline cycloalkane cluster (GN), gasoline aromatics cluster (GA), dry gas cluster (GS), and coke cluster (CK).

2. The method for evaluating the yield of catalytic cracking reaction products according to claim 1, characterized in that, In step (4), the specific operation of predicting the total content of each raw material based on its properties is as follows: 4.1) Using existing analytical methods, obtain the distillation range distribution of the feedstock oil, and calculate the mass content of the residue layer (H), wax oil layer (M), and diesel oil layer (L) using the following formula: AND H =100-Y M -AND L (9) In the formula, X represents the distillation range of the residue oil layer (H), wax oil layer (M), or diesel oil layer (L), and a sampling point is set at every 10% interval from the initial boiling point to the final boiling point of the residue oil layer (H), wax oil layer (M), or diesel oil layer (L), for a total of 11 sampling points including the initial boiling point and the final boiling point of the distillation range; i and j represent the i-th and j-th points, respectively. 4.2) Calculate the average number of carbon and hydrogen atoms in the crude oil molecule. The average number of C atoms in the crude oil molecule: In the formula, MW represents the relative molecular mass of one molecule in the feedstock, with units of kg / kmol, and y C This indicates the mass content of carbon in the crude oil, expressed in w%, MW, and y. C All were obtained through instrumental analysis; The average number of H atoms in the crude oil molecule: In the formula, MW represents the relative molecular mass of one molecule in the feedstock, with units of kg / kmol, and y H This indicates the mass content of hydrogen in the crude oil, expressed in w%, MW, and y. H All were obtained through instrumental analysis; 4.3) Calculate the naphthenic carbon content and aromatic carbon content in the feedstock oil respectively. (The text abruptly ends here, likely due to an incomplete sentence or missing information.) N and aromatic carbon content C A The unit is w%, and the specific formulas are as follows: In the formula, This indicates the density of the crude oil at 20℃, expressed in kg / m³. 3 γ 100 This indicates the viscosity of the feedstock oil at 100℃, in mm. 2 / s, these two parameters are obtained through instrumental analysis; 4.4) Calculate the mass content of each raw material aggregate separately. In the residue layer (H), the mass contents of the total alkane (HP), total cycloalkanes (HN), total aromatics (HA), and total aromatics (FA) in the asphaltene are as follows: AND HP =Y H -AND FA -AND HN -AND HA (17) In the wax oil layer (M), the mass contents of the wax oil layer alkane aggregate (MP), wax oil layer cycloalkanes aggregate (MN), and wax oil layer aromatics aggregate (MA) are as follows: AND MP =Y M -AND MN -AND MA (20) In the diesel fuel layer (L), the mass contents of the diesel fuel alkane aggregate (LP), diesel fuel cycloalkane aggregate (LN), diesel fuel monocyclic aromatic hydrocarbon aggregate (LA), and diesel fuel bicyclic or higher aromatic hydrocarbon aggregate (LHA) are as follows: AND LP And L -AND LN -AND LA -AND LHA (24) In the above formulas, CCR is the mass content of residual carbon in the feedstock, and this parameter is obtained through instrumental analysis. 4.5) Substitute the mass content of each raw material aggregate calculated in step 4.4) into the catalytic cracking kinetic model in step (3) to predict the yield of each product aggregate.

3. The method for evaluating the yield of catalytic cracking reaction products according to claim 2, characterized in that, The specific operation for predicting the aggregate yield of each product in step 4.5) is as follows: The initial content of each raw material aggregate is calculated according to formula (14)-formula (24), and then the yield of each product aggregate is calculated by formula (5). In the formula: Y = [y HP y HN y HA y FA y MP y MN y MA y LP y LN y LA y LHA y GP y GO y GN y GA y GAS y CK ] T This is the mass fraction vector of each ensemble component; P represents the system pressure, R is the gas constant, 8.314 J / (mol·K); T represents the system temperature, X = x / H represents the dimensionless relative distance at section x in the bed, H represents the total length of the catalyst bed, and S... WH Represents the actual gravity-time space velocity, and K is the reaction rate constant matrix; y i M represents the total mass content of the i-th set. i The average relative molecular mass of the i-th set; the concentration a of the i-th set. i The calculation formula is: In equation (4), a i The unit is molesi / g; K is the reaction rate constant matrix, which is composed of the reaction rate constants k of all reactions in the catalytic cracking reaction network. j constitute, In equation (6), A j E represents the pre-exponential factor of the j-th reaction in the catalytic cracking reaction network. j k represents the activation energy of the j-th reaction. j Let be the reaction rate constant of the j-th reaction in the constructed catalytic cracking reaction network.

4. The method for evaluating the yield of catalytic cracking reaction products according to claim 1, characterized in that, The calculation method for the content of each component gas in the dry gas aggregate (GS) is as follows: based on the structural group composition of the feedstock oil and the reaction temperature, reaction pressure, reaction time and catalyst-to-oil ratio of the catalytic cracking reaction, calculate the yield of hydrogen, methane, ethane, ethylene, propane, propylene, n-butane, isobutane, isobutene, n-butene, cis-butene and trans-butene. Reaction correlation: y i =(C(i,1)·T+C(i,2)·P+C(i,3)·t+C(i,4)·CO)·(y CPH+CPM ·C(i,5)+y CNH+CNM ·C(i,6)+y CAH+CAM+HA C(i,7)) / 100 (25) In the formula C(11，7)＝[0.04，0.10，0.09，0.37，0.11，0.02，0.10，0.01，0.60，0.13，0.02；0.06，0.38，0.11，0.36，0.01，0.31，0.37，-0.06，0.13，0.05，-0.96；-0.43，0.08，-0.23，-0.02，0.42，0.20，1.00，0.30，0.01，-0.15，-0.49；0.01，-0.39，0.10，-0.86，-0.03，-0.23，-0.12，-0.09，0.66，0.03，0.23；0.15，0.02，-0.23，-0.48，0.20，-0.30，0.03，-0.33，-0.14，0.08，-0.39；0.87，0.20，0.47，0.15，0.17，-0.13，-0.27，-0.44，-0.40，-0.13，-0.27；-0.29，0.22，-0.09，0.12，0.56，0.34，1.14，0.44，0.15，-0.01，-0.35]； T represents the reaction temperature in °C; P represents the reaction pressure in MPa; t represents the reaction time in seconds; CO represents the catalyst-to-oil ratio in kg catalyst / kg feedstock; y i This represents the yield of the i-th gas, expressed in w%.

5. The method for evaluating the yield of catalytic cracking reaction products according to claim 1, characterized in that, The catalytic cracking reaction network established in step (2) is as follows: The residue oil layer alkane aggregate (HP) was reacted with the wax oil layer alkane aggregate (MP), diesel alkane aggregate (LP), gasoline alkane aggregate (GP), gasoline aromatics aggregate (GA), gasoline olefins aggregate (GO), gasoline cycloalkanes aggregate (GN), dry gas aggregate (GS), and coke aggregate (CK), respectively. The residue oil layer cycloalkane aggregate (HN) was reacted with the diesel cycloalkane aggregate (LN), the wax oil layer cycloalkane aggregate (MN), the gasoline alkane aggregate (GP), the gasoline olefin aggregate (GO), the dry gas aggregate (GS), the gasoline cycloalkane aggregate (GN), the gasoline aromatics aggregate (GA), and the coke aggregate (CK), respectively. The aromatic hydrocarbon cluster (HA) of the residue oil layer is reacted with the gasoline alkane cluster (GP), the gasoline olefin cluster (GO), the aromatic hydrocarbon cluster of the wax oil layer (MA), the monocyclic aromatic hydrocarbon cluster of the diesel layer (LA), the gasoline cycloalkane cluster (GN), the coke cluster (CK), the gasoline aromatic hydrocarbon cluster (GA), and the bicyclic or higher aromatic hydrocarbon cluster of the diesel layer (LHA), respectively. The aromatic hydrocarbon aggregate (FA) in the asphaltene was reacted with the aromatic hydrocarbon aggregate (MA) in the wax oil layer, the bicyclic or higher aromatic hydrocarbon aggregate (LHA) in the diesel layer, and the coke aggregate (CK), respectively. The wax oil layer alkane aggregate (MP) was reacted with dry gas aggregate (GS), diesel alkane aggregate (LP), gasoline olefin aggregate (GO), gasoline aromatics aggregate (GA), gasoline cycloalkane aggregate (GN) and coke aggregate (CK), respectively. The wax oil layer cycloalkane aggregate (MN) was reacted with the gasoline alkane aggregate (GP), dry gas aggregate (GS), gasoline olefin aggregate (GO), diesel cycloalkane aggregate (LN), gasoline cycloalkane aggregate (GN), gasoline aromatics aggregate (GA), and coke aggregate (CK), respectively. The aromatic hydrocarbon clusters (MA) of the wax oil layer are reacted with the dry gas clusters (GS), the monocyclic aromatic hydrocarbon clusters (LA) of the diesel layer, the bicyclic or higher aromatic hydrocarbon clusters (LHA) of the diesel layer, the gasoline aromatic hydrocarbon clusters (GA), and the coke clusters (CK), respectively. Diesel alkane aggregate (LP) was reacted with dry gas aggregate (GS), gasoline alkane aggregate (GP), gasoline olefin aggregate (GO), gasoline cycloalkane aggregate (GN), gasoline aromatics aggregate (GA), and coke aggregate (CK), respectively. Diesel cycloalkane aggregate (LN) was reacted with gasoline alkane aggregate (GP), dry gas aggregate (GS), gasoline olefin aggregate (GO), gasoline cycloalkane aggregate (GN), gasoline aromatics aggregate (GA), and coke aggregate (CK), respectively. The diesel layer monocyclic aromatic hydrocarbon cluster (LA) was reacted with the dry gas cluster (GS), the gasoline aromatic hydrocarbon cluster (GA), and the coke cluster (CK), respectively. The diesel fuel layer bicyclic or higher aromatic hydrocarbon cluster (LHA) was reacted with the dry gas cluster (GS) and the coke cluster (CK), respectively. Gasoline alkane aggregates (GP) are reacted with dry gas aggregates (GS) and coke aggregates (CK), respectively. Gasoline olefin aggregates (GO) are reacted with gasoline alkane aggregates (GP), dry gas aggregates (GS), and coke aggregates (CK), respectively. Gasoline cycloalkane aggregates (GN) are reacted with dry gas aggregates (GS), gasoline aromatics aggregates (GA), and coke aggregates (CK), respectively. Gasoline aromatic hydrocarbon aggregates (GA) are reacted with dry gas aggregates (GS) and coke aggregates (CK), respectively.

6. The method for evaluating the yield of catalytic cracking reaction products according to claim 1, characterized in that, The basic equations for establishing the catalytic cracking kinetic model in step (3) are as follows: In the formula, a i Let represent the concentration of the i-th ensemble, in molesi / g; P represent the system pressure; R is the gas constant, 8.314 J / (mol·K); T represents the system temperature; X = x / H represents the dimensionless relative distance at section x in the bed; H represents the total length of the catalyst bed; S WH Represents the actual weight-time space velocity, K is the reaction rate constant matrix, and a = [a1, ..., a2]. ni ] T This is the lumped component concentration vector.

7. The method for evaluating the yield of catalytic cracking reaction products according to claim 6, characterized in that, The basic equations of the catalytic cracking kinetic model are derived from the continuity equation and the reaction rate equation: The continuity equation is: The reaction rate equation is: In the formula, subscript i represents the lumped component, subscript j represents the j-th reaction, and ρ represents the density of the oil-gas mixture in g / cm³. 3 a i G represents the total concentration of the i-th set, in molesi / g, t represents the reaction time, and G v This represents the mass flow rate of oil and gas across a cross section, in g / (cm). 2 ·h), where x represents the distance from the riser inlet into the reactor, R i Let n represent the reaction rate of the i-th set, and n represent the reaction rate of the set. r v represents the number of reactions. i,j R represents the stoichiometric coefficient of i lumped in reaction j. j k represents the reaction rate of reaction j. j ρ represents the rate constant of reaction j. c Catalyst density relative to reactor volume, in g / cm³ 3 ε represents the porosity, P represents the system pressure, R is the gas constant, and T represents the system temperature. This indicates the effect of catalyst coking on activity, and assumes that the catalyst coking rate is only a function of the catalyst residence time, t. c η is the catalyst residence time, and β is the catalyst deactivation constant.

8. A model-building apparatus for evaluating the yield of catalytic cracking reaction products, characterized in that: The model building apparatus includes a modeling module that uses the evaluation method for the yield of catalytic cracking reaction products according to any one of claims 1-7 for modeling.

9. An electronic device, characterized in that, It includes a processor and a memory, the memory storing computer-readable instructions that, when executed by the processor, enable the evaluation method of any one of claims 1-7 to be performed.

10. A readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it performs the steps of the evaluation method according to any one of claims 1-7.

Images