A prediction method and device for molecular catalytic cracking reaction products based on pressure changes

Abstract

An embodiment of the present application provides a method and device for predicting molecular-level catalytic cracking reaction products based on pressure changes. Among them, the method for predicting molecular-level catalytic cracking reaction products based on pressure changes includes: dividing the cracking reaction zone of the riser reactor into multiple differential units; predicting the products of the first catalytic cracking reaction of the molecular component material in the first differential unit according to the first pressure value by using the molecular dynamics reaction equation; obtaining the second pressure value at the inlet of the second differential unit according to the molar content and molecular weight of each product molecule generated by the first catalytic cracking reaction; predicting the products of the second catalytic cracking reaction of the molecular component material in the second differential unit according to the second pressure value by using the molecular dynamics reaction equation until the products of the cracking reaction in each differential unit among the multiple differential units are predicted.

Description

Technical Field

[0001] The present application relates to the technical field of oil refining and processing, and particularly relates to a method and device for predicting the products of molecular catalytic cracking reactions based on pressure changes. Background Art

[0002] The core idea of the structure oriented lumping (SOL) method is that all complex hydrocarbon molecules in oil products can be disassembled into molecular fragments or molecular structure groups, and this molecular structure group is called a structure vector. In China, in recent years, many studies on the SOL model have been carried out, such as predicting product yields and properties, optimizing the configuration of raw materials, and adjusting processing plans, etc., and the application effects are relatively good.

[0003] In the prior art, due to the complexity of the SOL model itself and the excessive computational amount of the reaction model at the molecular level, restricted by the computing power of the computer in the early stage, the reaction process model was simplified, without considering the pressure gradient distribution caused by the changes in the density and molecular weight of the product mixture during the catalytic cracking reaction in the riser reactor, and it was considered that the pressure in the reactor remained constant throughout the reaction process. However, in the actual catalytic cracking reaction process, as the reaction progresses, the density and molecular weight of the reaction products change continuously, and the pressure in the riser reactor will change accordingly. And the prediction of the reaction products of molecular cracking is related to the pressure. Assuming that the pressure in the riser reactor remains constant during the reaction process will result in a large error in the calculation result of the reaction rate, thereby reducing the accuracy of the prediction result of the reaction process model.

[0004] Therefore, how to more accurately predict the products of molecular component materials using the reaction process model based on the pressure change in the riser reactor during the catalytic cracking reaction process is a technical problem to be solved urgently. Summary of the Invention

[0005] The embodiments of the present application provide a method and device for predicting the products of molecular catalytic cracking reactions based on pressure changes, which are used to more accurately predict the products of molecular component materials using the reaction process model based on the pressure change in the riser reactor during the catalytic cracking reaction process.

[0006] One embodiment of the present application provides a method for predicting the products of a molecular-level catalytic cracking reaction based on pressure changes. The method includes: dividing the cracking reaction zone of the riser reactor into multiple differential units; predicting the products of the first catalytic cracking reaction of the molecular component material in the first differential unit according to the first pressure value using the molecular dynamics reaction equation; wherein, the first pressure value is the pressure value at the inlet of the cracking reaction zone, and the first differential unit is the differential unit located at the inlet of the cracking reaction zone; obtaining the second pressure value at the inlet of the second differential unit according to the molecular structure and molar content of each product molecule generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction; predicting the products of the second catalytic cracking reaction of the molecular component material in the second differential unit according to the second pressure value using the molecular dynamics reaction equation until the prediction of the products of the cracking reaction in each differential unit among the multiple differential units is completed; wherein, the second differential unit is adjacent to the first differential unit.

[0007] In some embodiments, the obtaining the second pressure value at the inlet of the second differential unit according to the molecular structure and molar content of each product molecule generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction includes: calculating the molecular weight of each product molecule generated by the first catalytic cracking reaction according to the molecular structure of each product molecule generated by the first catalytic cracking reaction and the molar mass of the structural groups represented by the structure-oriented lumping method; obtaining the average molecular weight of the products generated by the first catalytic cracking reaction according to the molar content and molecular weight of each product molecule generated by the first catalytic cracking reaction; obtaining the mixed density of the products generated by the first catalytic cracking reaction according to the average molecular weight of the products generated by the first catalytic cracking reaction; calculating the second pressure value according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction.

[0008] In some embodiments, the obtaining the average molecular weight of the products generated by the first catalytic cracking reaction according to the molar content and molecular weight of each product molecule generated by the first catalytic cracking reaction includes:

[0009] Obtaining the average molecular weight of the products generated by the first catalytic cracking reaction using the following formula:

[0010]

[0011] Wherein, is the molar content of the i-th product molecule generated by the first catalytic cracking reaction, is the molecular weight of the i-th product molecule generated by the first catalytic cracking reaction, is the average molecular weight of the products generated by the first catalytic cracking reaction.

[0012] In some embodiments, obtaining the mixed density of the products generated by the first catalytic cracking reaction according to the average molecular weight of the products generated by the first catalytic cracking reaction includes:

[0013] Obtaining the mixed density of the products generated by the first catalytic cracking reaction through the following formula :

[0014]

[0015] where P is the system pressure in the cracking reaction zone of the riser reactor, is the average molecular weight of the products generated by the first catalytic cracking reaction, T is the temperature in the cracking reaction zone of the riser reactor, and R is the gas constant.

[0016] In some embodiments, calculating the second pressure value according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction includes: obtaining a first pressure drop value according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by it, using a preset pressure drop calculation method; obtaining the second pressure value according to the first pressure drop value.

[0017] In some embodiments, the preset pressure drop calculation method is calculated based on the pressure drop calculation formula of the vertical pipe dilute phase conveying system. Obtaining the first pressure drop value according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction using the preset pressure drop calculation method includes: obtaining the pressure drop calculation term of the products generated by the first catalytic cracking reaction in the pressure drop calculation formula of the vertical pipe dilute phase conveying system according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction; obtaining the pressure drop calculation term of the solid catalyst carried by the products generated by the first catalytic cracking reaction in the pressure drop calculation formula of the vertical pipe dilute phase conveying system according to the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction; obtaining the first pressure drop value according to the pressure drop calculation term of the products generated by the first catalytic cracking reaction and the pressure drop calculation term of the solid catalyst carried by the products generated by the first catalytic cracking reaction.

[0018] In some embodiments, the calculation formula of the second pressure value is as follows:

[0019]

[0020] Among them, is the first pressure value, is the first pressure drop value, is the second pressure value.

[0021] One embodiment of the present application provides a prediction device for molecular catalytic cracking reaction products based on pressure changes. The device includes: a differential unit division module for dividing the cracking reaction zone of the riser reactor into multiple differential units; a first prediction module for predicting the products of the first catalytic cracking reaction of the molecular component material in the first differential unit according to the first pressure value by using the molecular dynamics reaction equation; wherein, the first pressure value is the pressure value at the inlet of the cracking reaction zone, and the first differential unit is the differential unit located at the inlet of the cracking reaction zone; a second pressure value obtaining module for obtaining the second pressure value at the inlet of the second differential unit according to the molar content and molecular weight of each product molecule generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction; a second prediction module for predicting the products of the second catalytic cracking reaction of the molecular component material in the second differential unit according to the second pressure value by using the molecular dynamics reaction equation until the prediction of the products of the cracking reaction in each of the multiple differential units is completed; wherein, the second differential unit is adjacent to the first differential unit.

[0022] The above technical solutions provided by the embodiments of the present application have at least the following advantages compared with the prior art:

[0023] In the embodiments provided by the present application, the cracking reaction zone of the riser reactor is divided into multiple differential units; according to the first pressure value, the products of the first catalytic cracking reaction of the molecular component material in the first differential unit are predicted by using the molecular dynamics reaction equation; according to the molar content and molecular weight of each product molecule generated by the first catalytic cracking reaction, the second pressure value at the inlet of the second differential unit is obtained; according to the second pressure value, the products of the second catalytic cracking reaction of the molecular component material in the second differential unit are predicted by using the molecular dynamics reaction equation until the prediction of the products of the cracking reaction in each of the multiple differential units is completed. Since the change in gas pressure caused by the change in the concentration and molecular weight of the reaction products is considered during the prediction process, more accurate prediction results can be obtained. Description of the Drawings

[0024] This application will be further described by way of exemplary embodiments, which will be described in detail with reference to the accompanying drawings. These embodiments are not restrictive, and in these embodiments, the same numbers represent the same structures, where:

[0025] Figure 1 is a schematic diagram of an application scenario of a method for predicting molecular-level catalytic cracking reaction products based on pressure changes according to some embodiments of the present application;

[0026] Figure 2 is an exemplary flowchart of a method for predicting molecular-level catalytic cracking reaction products based on pressure changes according to some embodiments of the present application;

[0027] Figure 3 is an exemplary flowchart of multiple differential units according to some embodiments of the present application;

[0028] Figure 4 is an exemplary schematic diagram of a device for predicting molecular-level catalytic cracking reaction products based on pressure changes according to some embodiments of the present application;

[0029] Figure 5 is an exemplary schematic diagram of 24 groups included in the structure-oriented lumping method according to some embodiments of the present application. Detailed Embodiments

[0030] To more clearly illustrate the technical solutions of the embodiments of the present application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some examples or embodiments of the present application. For those of ordinary skill in the art, without creative efforts, the present application can also be applied to other similar scenarios based on these drawings. Unless obvious from the language context or otherwise stated, the same reference numerals in the figures represent the same structures or operations.

[0031] It should be understood that the "system", "device", "unit" and / or "module" used herein is a way to distinguish different components, elements, parts, portions or assemblies at different levels. However, if other words can achieve the same purpose, they can be replaced by other expressions.

[0032] As shown in the present application and the claims, unless the context clearly indicates an exception, words such as "a", "an", "one" and / or "the" are not specifically singular and may also include plural. Generally speaking, the terms "comprising" and "including" only indicate the inclusion of the clearly identified steps and elements, and these steps and elements do not constitute an exclusive list. The method or device may also include other steps or elements.

[0033] In this application, flowcharts are used to illustrate the operations performed by the system according to the embodiments of this application. It should be understood that the preceding or subsequent operations do not necessarily need to be executed precisely in sequence. On the contrary, the steps can be processed in reverse order or simultaneously. At the same time, other operations can also be added to these processes, or one or more steps can be removed from these processes.

[0034] Figure 1 It is a schematic diagram of the application scenario of the prediction method for the molecular catalytic cracking reaction product based on pressure change shown in some embodiments of this application.

[0035] As Figure 1 shown, the application scenario can include a server 110, a terminal 120, and a network 130.

[0036] In some embodiments, data or information can be exchanged between the server 110 and the terminal 120 through the network 130. For example, the server 110 can obtain information and / or data in the terminal 120 through the network 130, or can send information and / or data to the terminal 120 through the network 130.

[0037] The terminal 120 is an electronic device used by the user to predict the products of the molecular component material in the catalytic cracking reaction process. In some embodiments, the terminal 120 can predict the products of the molecular component material in the catalytic cracking reaction process according to the method provided in the embodiments of this application. When the computing resources of the terminal 120 are limited, the server 110 can predict the products of the molecular component material in the catalytic cracking reaction process according to the method provided in the embodiments of this application, and return the prediction result to the terminal 120, so that the terminal 120 can display the prediction result to the user. The terminal 120 can be one or any combination of devices with input and / or output functions such as mobile devices, tablet computers, etc.

[0038] The server 110 can be a single server or a server group. The server group can be centralized or distributed (for example, the server 110 can be a distributed system), can be dedicated, or can also provide services by other devices or systems simultaneously. In some embodiments, the server 110 can be regional or remote. In some embodiments, the server 110 can be implemented on a cloud platform or provided in a virtual manner. Only by way of example, the cloud platform can include a private cloud, a public cloud, a hybrid cloud, a community cloud, a distributed cloud, an internal cloud, a multi-layer cloud, etc. or any combination thereof.

[0039] In some embodiments, network 130 may be any one or more of a wired network or a wireless network. For example, network 130 may include a local area network (LAN), a wide area network (WAN), a wireless local area network (WLAN), a metropolitan area network (MAN), etc., or any combination thereof.

[0040] For ease of understanding, the technical solutions of this application are introduced below in conjunction with the accompanying drawings and embodiments.

[0041] Figure 2 is an exemplary flowchart of a method for predicting the products of a molecular-level catalytic cracking reaction based on pressure changes as shown in some embodiments of this application. As Figure 2 shown, the method for predicting the products of a molecular-level catalytic cracking reaction based on pressure changes includes the following steps:

[0042] Step S210: Divide the cracking reaction zone of the riser reactor into multiple differential units.

[0043] As Figure 3 shown, the riser reactor is sequentially provided with a pre-lift section, a feed section, and a cracking reaction zone from bottom to top. Figure 3 Only as an example, the actual length of the cracking reaction zone is much larger than the diameter of the riser reactor. The riser reactor is for vertical pneumatic conveying. After the molecular component material enters the cracking reaction zone from the feed section, it undergoes a catalytic cracking reaction during the rising process.

[0044] In the specific implementation process, as Figure 3 shown, the cracking reaction zone can be divided into multiple continuously arranged differential units. The length of each differential unit can be the same or different, and is not limited by the description in this specification. The actual catalytic cracking reaction is a continuous process. In the embodiments of this application, in order to obtain the pressure value in the riser reactor during the reaction process, according to the differential units, the catalytic cracking reaction process is sequentially divided into the first catalytic cracking reaction, the second catalytic cracking reaction... As the molecular component material rises along the cracking reaction zone, the cracking reaction depth gradually increases, and the density and molecular weight of the reaction products also change continuously, and the pressure in the cracking reaction zone also changes accordingly.

[0045] Step S220: Predict the products of the first catalytic cracking reaction of the molecular component material in the first differential unit by using the molecular dynamics reaction equation according to the first pressure value.

[0046] The first pressure value is the pressure value at the inlet of the cracking reaction zone. In the specific implementation process, the first pressure value can be obtained by instrument measurement. As Figure 3As shown, the first differential unit is the differential unit located at the inlet of the cracking reaction zone. The catalytic cracking reaction mainly involves cracking and is a reaction with an increase in the number of moles. As the molecular component material ascends along the riser reactor, the reaction depth gradually increases, the concentration of small molecule products gradually increases, the system expands, and the pressure drop generated gradually increases. Therefore, after the first catalytic cracking reaction, the pressure value at the outlet of the first differential unit changes compared with the pressure value at the inlet.

[0047] Step S230: Obtain the second pressure value at the inlet of the second differential unit according to the molecular structure and molar content of each product molecule generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the product generated by the first catalytic cracking reaction.

[0048] In some embodiments, according to the molecular structure of each product molecule generated by the first catalytic cracking reaction, the molecular weight of each product molecule generated by the first catalytic cracking reaction can be calculated using the molar mass of the structural groups represented by the structure-oriented lumping method.

[0049] The 24 groups of the structure-oriented lumping method are as Figure 5 shown, where A6 is a benzene ring; A4 is the four-carbon aromatic ring increment attached to other aromatic rings; A2 is the two-carbon aromatic ring increment; N6 and N5 are aliphatic rings with 6 carbons and 5 carbons respectively; N4, N3, N2, and N1 are the aliphatic ring increments representing 4 carbons, 3 carbons, 2 carbons, and 1 carbon connected to aromatic rings or cycloalkane rings respectively; R is the total number of carbons except for the carbons on the ring; me refers to the number of methyl groups attached to the molecular aromatic ring or aliphatic ring; br is the number of alkyl substituents connected to hydrocarbons, alkenes, or alkyl branches; AA represents the bridge bond between two rings; IH is the hydrogen increment used to specify the degree of molecular unsaturation (except for the unsaturation on the aromatic ring); NS, NN, and NO are sulfur, nitrogen, and oxygen atoms connecting two carbon atoms; RS, RN, and RO represent sulfur, nitrogen, and oxygen atoms between carbon and hydrogen respectively; AN represents the nitrogen atom on the aromatic ring; KO represents the oxygen atom of the carbonyl or aldehyde group; Ni and V represent metal nickel and vanadium atoms.

[0050] In some embodiments, the average molecular weight of the product generated by the first catalytic cracking reaction can be obtained according to the molar content and molecular weight of each product molecule generated by the first catalytic cracking reaction; the mixed density of the product generated by the first catalytic cracking reaction can be obtained according to the average molecular weight of the product generated by the first catalytic cracking reaction.

[0051] In the specific implementation process, the average molecular weight of the product generated by the first catalytic cracking reaction can be obtained using the following formula:

[0052] (1)

[0053] Among them, is the molar content of the i-th product molecule generated by the first catalytic cracking reaction, is the molecular weight of the i-th product molecule generated by the first catalytic cracking reaction, is the average molecular weight of the products generated by the first catalytic cracking reaction.

[0054] In some embodiments, the mixed density of the products generated by the first catalytic cracking reaction can be obtained through the following formula :

[0055] (2)

[0056] Among them, P is the system pressure in the cracking reaction zone of the riser reactor, is the average molecular weight of the products generated by the first catalytic cracking reaction, T is the temperature in the cracking reaction zone of the riser reactor, and R is the gas constant.

[0057] In some embodiments, according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, as well as the density and flow rate of the solid catalyst carried by them, using a preset pressure drop calculation method, a first pressure drop value can be obtained. Only as an example, the preset pressure drop calculation method can be the pressure drop calculation method for the dilute-phase transportation system of a vertical pipe. The formula for the pressure drop calculation method of the dilute-phase transportation system of a vertical pipe is as follows:

[0058] (3)

[0059] In formula (3), is the total pressure drop. The pressure drop calculation terms for the products generated by the first catalytic cracking reaction include: the pressure drop of the accelerating gas , the pressure drop of the gas friction with the pipe wall , the static head of the gas . The pressure drop calculation terms for the solid catalyst carried by the products generated by the first catalytic cracking reaction include: the pressure drop of the accelerating solid , the pressure drop of the solid friction with the pipe wall , the static head of the solid .

[0060] Among them, the pressure drop of the accelerating gas is calculated through the following formula:

[0061] (4)

[0062] Among them, the pressure drop of the accelerating solid is calculated through the following formula:

[0063] (5)

[0064] Among them, the frictional pressure drop between the gas and the pipe wall is calculated by the following formula:

[0065] (6)

[0066] Among them, the frictional pressure drop between the solid and the pipe wall is calculated by the following formula:

[0067] (7)

[0068] Among them, the static head of the solid is calculated by the following formula:

[0069] (8)

[0070] Among them, the static head of the gas is calculated by the following formula:

[0071] (9)

[0072] Among them, is the linear velocity of the gas mixture, is the density of the gas mixture, is the gas friction coefficient, is the slip velocity, is the true velocity of the gas flow, is the true velocity of the catalyst solid particles, is the solid mass flow rate, is the acceleration due to gravity, L is the length of the riser section, is the diameter of the riser section;

[0073] Among them, is calculated by the following formula:

[0074] (10)

[0075] Among them, is calculated by the following formula:

[0076] (11)

[0077] In some embodiments, the pressure drop calculation term of the product generated by the first catalytic cracking reaction in the pressure drop calculation formula of the vertical pipe dilute phase conveying system can be obtained according to the average molecular weight and mixed density of the product generated by the first catalytic cracking reaction; the pressure drop calculation term of the solid catalyst carried by the product generated by the first catalytic cracking reaction in the pressure drop calculation formula of the vertical pipe dilute phase conveying system can be obtained according to the density and flow rate of the solid catalyst carried in the product. For example, the following method can be used to obtain each pressure drop calculation term: the molar flow rate of the product generated by the first catalytic cracking reaction is obtained according to the average molecular weight of the product generated by the first catalytic cracking reaction, and the pressure drop calculation term of the product generated by the first catalytic cracking reaction is calculated according to the mixed density and flow rate of the product generated by the first catalytic cracking reaction; the pressure drop calculation term of the solid is calculated according to the density and flow rate of the solid catalyst carried in the product generated by the first catalytic cracking reaction.

[0078] After obtaining the pressure drop calculation term of the product generated by the first catalytic cracking reaction and the pressure drop calculation term of the solid catalyst carried by the product generated by the first catalytic cracking reaction, formula (3) can be used to calculate the first pressure drop value (corresponding to the total pressure drop ).

[0079] After calculating the first pressure drop value, the second pressure value can be calculated according to the following formula:

[0080] (12)

[0081] In formula (12), is the first pressure value, is the first pressure drop value, is the second pressure value.

[0082] Step S240, according to the second pressure value, using the molecular dynamics reaction equation, predict the product of the second catalytic cracking reaction of the molecular component material in the second differential unit until the prediction of the product of the cracking reaction in each of the multiple differential units is completed.

[0083] The second differential unit is adjacent to the first differential unit. In a specific implementation process, the method described in steps S210 to S240 can be repeatedly executed to obtain the pressure value at the inlet of each differential unit, and the product of the cracking reaction occurring in each differential unit can be predicted according to the pressure value at the inlet.

[0084] Figure 4 is an exemplary schematic diagram of a prediction device for molecular-level catalytic cracking reaction products based on pressure changes according to some embodiments of the present application.

[0085] As Figure 4As shown in the figure, the prediction device for the molecular catalytic cracking reaction products based on pressure change includes: a differential unit division module 410, a first prediction module 420, a second pressure value acquisition module 430, and a second prediction module 440.

[0086] The differential unit division module 410 is configured to divide the cracking reaction zone of the riser reactor into a plurality of differential units.

[0087] The first prediction module 420 is configured to predict the products of the first catalytic cracking reaction of the molecular component material in the first differential unit according to the first pressure value by using the molecular dynamics reaction equation; wherein, the first pressure value is the pressure value at the inlet of the cracking reaction zone, and the first differential unit is the differential unit located at the inlet of the cracking reaction zone.

[0088] The second pressure value acquisition module 430 is configured to obtain the second pressure value at the inlet of the second differential unit according to the molecular structure and molar content of each product molecule generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction.

[0089] The second prediction module 440 is configured to predict the products of the second catalytic cracking reaction of the molecular component material in the second differential unit according to the second pressure value by using the molecular dynamics reaction equation until the prediction of the products of the cracking reaction in each of the plurality of differential units is completed; wherein, the second differential unit is adjacent to the first differential unit.

[0090] In some embodiments, the obtaining the second pressure value at the inlet of the second differential unit according to the molecular structure and molar content of each product molecule generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction includes: calculating the molecular weight of each product molecule generated by the first catalytic cracking reaction according to the molecular structure of each product molecule generated by the first catalytic cracking reaction and the molar mass of the structural groups represented by the structure-oriented lumping method; obtaining the average molecular weight of the products generated by the first catalytic cracking reaction according to the molar content and molecular weight of each product molecule generated by the first catalytic cracking reaction; obtaining the mixed density of the products generated by the first catalytic cracking reaction according to the average molecular weight of the products generated by the first catalytic cracking reaction; and calculating the second pressure value according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction.

[0091] In some embodiments, obtaining the average molecular weight of the products formed by the first catalytic cracking reaction based on the molar content and molecular weight of each product molecule formed by the first catalytic cracking reaction includes:

[0092] Obtaining the average molecular weight of the products formed by the first catalytic cracking reaction using the following formula:

[0093]

[0094] Wherein, is the molar content of the i-th product molecule formed by the first catalytic cracking reaction, is the molecular weight of the i-th product molecule formed by the first catalytic cracking reaction, is the average molecular weight of the products formed by the first catalytic cracking reaction.

[0095] In some embodiments, obtaining the mixed density of the products formed by the first catalytic cracking reaction based on the average molecular weight of the products formed by the first catalytic cracking reaction includes:

[0096] Obtaining the mixed density of the products formed by the first catalytic cracking reaction through the following formula :

[0097]

[0098] Wherein, P is the system pressure in the cracking reaction zone of the riser reactor, is the average molecular weight of the products formed by the first catalytic cracking reaction, T is the temperature in the cracking reaction zone of the riser reactor, R is the gas constant.

[0099] In some embodiments, calculating the second pressure value based on the average molecular weight and mixed density of the products formed by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products formed by the first catalytic cracking reaction includes: based on the average molecular weight and mixed density of the products formed by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products formed by the first catalytic cracking reaction, using a preset pressure drop calculation method to obtain a first pressure drop value; obtaining the second pressure value based on the first pressure drop value.

[0100] In the embodiments of the above prediction device for molecular-level catalytic cracking reaction products based on pressure changes, the specific processing of each module and the technical effects brought thereby can be respectively referred to the relevant descriptions in the corresponding method embodiments, which will not be elaborated herein.

[0101] The basic concepts have been described above. Obviously, for those skilled in the art, the above detailed disclosure is merely an example and does not constitute a limitation to this application. Although not explicitly stated here, those skilled in the art may make various modifications, improvements, and corrections to this application. Such modifications, improvements, and corrections are suggested in this application, so such modifications, improvements, and corrections still fall within the spirit and scope of the exemplary embodiments of this application.

[0102] Meanwhile, this application uses specific terms to describe the embodiments of this application. Such as "one embodiment", "an embodiment", and / or "some embodiments" mean a certain feature, structure, or characteristic related to at least one embodiment of this application. Therefore, it should be emphasized and noted that the "one embodiment" or "an embodiment" or "an alternative embodiment" mentioned twice or more at different positions in this application does not necessarily refer to the same embodiment. In addition, certain features, structures, or characteristics in one or more embodiments of this application can be appropriately combined.

[0103] In addition, unless clearly stated in the claims, the order of the processing elements and sequences, the use of numerical letters, or the use of other names in this application is not used to limit the order of the processes and methods of this application. Although some currently considered useful invention embodiments are discussed through various examples in the above disclosure, it should be understood that such details only serve the purpose of illustration. The appended claims are not limited to the disclosed embodiments. On the contrary, the claims are intended to cover all modifications and equivalent combinations that conform to the essence and scope of the embodiments of this application. For example, although the system components described above can be implemented by hardware devices, they can also be implemented only through software solutions, such as installing the described system on existing servers or mobile devices.

[0104] Similarly, it should be noted that, in order to simplify the expression of the disclosure of this application and thus help the understanding of one or more invention embodiments, in the previous description of the embodiments of this application, sometimes multiple features are merged into one embodiment, drawing, or description thereof. However, this disclosure method does not mean that the features required by the object of this application are more than those mentioned in the claims. In fact, the features of the embodiments are less than all the features of the individual embodiments disclosed above.

[0105] In some embodiments, numbers are used to describe components and the quantity of attributes. It should be understood that such numbers used in the description of embodiments are modified by the modifiers "about", "approximately" or "substantially" in some examples. Unless otherwise specified, "about", "approximately" or "substantially" indicate that the said numbers allow a variation of ±20%. Accordingly, in some embodiments, the numerical parameters used in the specification and claims are approximate values, which may vary according to the characteristics required by individual embodiments. In some embodiments, the numerical parameters should consider the specified significant digits and adopt the method of retaining general digits. Although the numerical ranges and parameters used in some embodiments of the present application to confirm the breadth of their scope are approximate values, in specific embodiments, such numerical settings are as precise as possible within the feasible range.

[0106] For each patent, patent application, patent application publication, and other materials cited in the present application, such as articles, books, specifications, publications, documents, etc., their entire contents are hereby incorporated into the present application by reference. Except for the application history documents that are inconsistent with or conflict with the content of the present application, and also except for the documents that limit the broadest scope of the claims of the present application (currently or subsequently appended to the present application). It should be noted that if there are any inconsistencies or conflicts between the descriptions, definitions, and / or uses of terms in the supplementary materials of the present application and the content described in the present application, the descriptions, definitions, and / or uses of terms in the present application shall prevail.

[0107] Finally, it should be understood that the embodiments described in the present application are only used to illustrate the principles of the embodiments of the present application. Other variations may also fall within the scope of the present application. Therefore, by way of example and not limitation, alternative configurations of the embodiments of the present application may be considered to be consistent with the teachings of the present application. Accordingly, the embodiments of the present application are not limited to the embodiments explicitly introduced and described in the present application.

Claims

1. A prediction method for molecular catalytic cracking reaction products based on pressure change, characterized in that, The method includes: Dividing the cracking reaction zone of the riser reactor into a plurality of differential units; Predicting the products of the first catalytic cracking reaction of the molecular component material in the first differential unit according to the first pressure value by using the molecular dynamics reaction equation; wherein, the first pressure value is the pressure value at the inlet of the cracking reaction zone, and the first differential unit is the differential unit located at the inlet of the cracking reaction zone; Obtaining the second pressure value at the inlet of the second differential unit according to the molecular structure and molar content of each product molecule generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction, including: Calculating the molecular weight of each product molecule generated by the first catalytic cracking reaction according to the molecular structure of each product molecule generated by the first catalytic cracking reaction and the molar mass of the structural groups represented by the structure-oriented lumping method; Obtaining the average molecular weight of the products generated by the first catalytic cracking reaction according to the molar content and molecular weight of each product molecule generated by the first catalytic cracking reaction, and including: Obtaining the average molecular weight of the products generated by the first catalytic cracking reaction by using the following formula: Wherein, is the molar content of the i-th product molecule generated by the first catalytic cracking reaction, is the molecular weight of the i-th product molecule generated by the first catalytic cracking reaction, is the average molecular weight of the products generated by the first catalytic cracking reaction; Obtaining the mixed density of the products generated by the first catalytic cracking reaction according to the average molecular weight of the products generated by the first catalytic cracking reaction, and including: The mixed density of the product generated by the first catalytic cracking reaction is obtained through the following formula : where P is the system pressure in the cracking reaction zone of the riser reactor, is the average molecular weight of the products generated by the first catalytic cracking reaction, T is the temperature in the cracking reaction zone of the riser reactor, R is the gas constant; Calculating the second pressure value according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction; Predicting the products of the second catalytic cracking reaction of the molecular component material in the second differential unit according to the second pressure value by using the molecular dynamics reaction equation until the prediction of the products of the cracking reaction in each differential unit among the plurality of differential units is completed; wherein, the second differential unit is adjacent to the first differential unit.

2. The method according to claim 1, characterized in that, The calculating the second pressure value according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by the products generated by the first catalytic cracking reaction, includes: Obtaining a first pressure drop value by using a preset pressure drop calculation method according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by it; Obtaining the second pressure value according to the first pressure drop value.

3. The method according to claim 2, characterized in that, The preset pressure drop calculation method is calculated based on the pressure drop calculation formula of the vertical pipe dilute phase conveying system. The obtaining a first pressure drop value by using the preset pressure drop calculation method according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction, and the density and flow rate of the solid catalyst carried by it, includes: Obtaining the pressure drop calculation term of the products generated by the first catalytic cracking reaction in the pressure drop calculation formula of the vertical pipe dilute phase conveying system according to the average molecular weight and mixed density of the products generated by the first catalytic cracking reaction; Based on the density and flow rate of the solid catalyst carried by the product generated from the first catalytic cracking reaction, obtain the pressure drop calculation term of the solid catalyst carried by the product generated from the first catalytic cracking reaction in the pressure drop calculation formula of the vertical pipe dilute phase conveying system; Based on the pressure drop calculation term of the product generated from the first catalytic cracking reaction and the pressure drop calculation term of the solid catalyst carried by the product generated from the first catalytic cracking reaction, obtain the first pressure drop value.

4. A prediction device for molecular catalytic cracking reaction products based on pressure change, characterized in that, The device includes: A differential unit division module for dividing the cracking reaction zone of the riser reactor into multiple differential units; A first prediction module for predicting the product of the first catalytic cracking reaction of the molecular component material in the first differential unit according to the first pressure value by using the molecular dynamics reaction equation; wherein, the first pressure value is the pressure value at the inlet of the cracking reaction zone, and the first differential unit is the differential unit located at the inlet of the cracking reaction zone; A second pressure value obtaining module for obtaining the second pressure value at the inlet of the second differential unit according to the molecular structure and molar content of each product molecule generated from the first catalytic cracking reaction, as well as the density and flow rate of the solid catalyst carried by it, and the density and flow rate of the solid catalyst carried by the product generated from the first catalytic cracking reaction, including: According to the molecular structure of each product molecule generated from the first catalytic cracking reaction, calculate the molecular weight of each product molecule generated from the first catalytic cracking reaction by using the molar mass of the structural groups represented by the structure-oriented lumping method; Based on the molar content and molecular weight of each product molecule generated from the first catalytic cracking reaction, obtain the average molecular weight of the product generated from the first catalytic cracking reaction, and include: Use the following formula to obtain the average molecular weight of the product generated from the first catalytic cracking reaction: wherein, is the molar content of the i-th product molecule generated by the first catalytic cracking reaction, is the molecular weight of the i-th product molecule generated by the first catalytic cracking reaction, is the average molecular weight of the products generated by the first catalytic cracking reaction; Based on the average molecular weight of the product generated from the first catalytic cracking reaction, obtain the mixed density of the product generated from the first catalytic cracking reaction, and include: The mixed density of the product formed by the first catalytic cracking reaction is obtained through the following formula :[[]]END]] wherein, P is the system pressure in the cracking reaction zone of the riser reactor, is the average molecular weight of the product formed by the first catalytic cracking reaction, T is the temperature in the cracking reaction zone of the riser reactor, R is the gas constant; Based on the average molecular weight and mixed density of the product generated from the first catalytic cracking reaction, as well as the density and flow rate of the solid catalyst carried by the product generated from the first catalytic cracking reaction, calculate the second pressure value; A second prediction module for predicting the product of the second catalytic cracking reaction of the molecular component material in the second differential unit according to the second pressure value by using the molecular dynamics reaction equation until the prediction of the product of the cracking reaction in each differential unit among the multiple differential units is completed; wherein, the second differential unit is adjacent to the first differential unit.

Images