CN112201309A - P-graph-based ethylene full-process superstructure model modeling method - Google Patents
P-graph-based ethylene full-process superstructure model modeling method Download PDFInfo
- Publication number
- CN112201309A CN112201309A CN202010557272.7A CN202010557272A CN112201309A CN 112201309 A CN112201309 A CN 112201309A CN 202010557272 A CN202010557272 A CN 202010557272A CN 112201309 A CN112201309 A CN 112201309A
- Authority
- CN
- China
- Prior art keywords
- ethylene
- graph
- superstructure
- ethane
- full
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 96
- 239000005977 Ethylene Substances 0.000 title claims abstract description 56
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 title claims abstract description 48
- 230000008569 process Effects 0.000 claims abstract description 65
- 238000005336 cracking Methods 0.000 claims abstract description 35
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 claims abstract description 24
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 238000007142 ring opening reaction Methods 0.000 claims abstract description 14
- 238000010791 quenching Methods 0.000 claims abstract description 8
- 230000000171 quenching effect Effects 0.000 claims abstract description 8
- 150000001335 aliphatic alkanes Chemical class 0.000 claims abstract description 4
- 150000001336 alkenes Chemical class 0.000 claims abstract description 4
- JRZJOMJEPLMPRA-UHFFFAOYSA-N olefin Natural products CCCCCCCC=C JRZJOMJEPLMPRA-UHFFFAOYSA-N 0.000 claims abstract description 4
- 239000002994 raw material Substances 0.000 claims description 31
- 239000000446 fuel Substances 0.000 claims description 9
- 230000010354 integration Effects 0.000 claims description 9
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical group [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 6
- 239000007858 starting material Substances 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 229910052799 carbon Inorganic materials 0.000 claims description 4
- 239000002253 acid Substances 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000003197 catalytic effect Effects 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract description 16
- 230000008901 benefit Effects 0.000 abstract description 15
- 238000012216 screening Methods 0.000 abstract description 5
- 238000005457 optimization Methods 0.000 description 13
- 238000007670 refining Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 238000000926 separation method Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 230000007613 environmental effect Effects 0.000 description 6
- 239000007789 gas Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000003915 liquefied petroleum gas Substances 0.000 description 2
- 125000004805 propylene group Chemical group [H]C([H])([H])C([H])([*:1])C([H])([H])[*:2] 0.000 description 2
- 238000004064 recycling Methods 0.000 description 2
- 238000002910 structure generation Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 208000037656 Respiratory Sounds Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- -1 and meanwhile Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003776 cleavage reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003502 gasoline Substances 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013386 optimize process Methods 0.000 description 1
- QQONPFPTGQHPMA-UHFFFAOYSA-N propylene Natural products CC=C QQONPFPTGQHPMA-UHFFFAOYSA-N 0.000 description 1
- 238000000197 pyrolysis Methods 0.000 description 1
- 206010037833 rales Diseases 0.000 description 1
- 230000007017 scission Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 239000013585 weight reducing agent Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
The invention discloses a P-graph-based modeling method for a full-flow superstructure of ethylene, which is used for establishing a model by taking a process of converting straight-chain alkane into olefin as a target object, and comprises four parts of cracking, quenching, compressing and separating, wherein the quenching and compressing processes are simplified into an operation unit type node, ethane circulation existing in the process is subjected to ring opening, and an optimal solution set and a near suboptimal solution set obtained on the basis of a P-graph algorithm are analyzed. The invention provides an ethylene full-process superstructure model based on a P-graph framework, and has the advantage of being capable of obtaining a series of suboptimal solutions. In the modeling process, a method for screening related equipment of a complex industrial process is provided according to the characteristics of 'three-pass-one-reverse' and a problem to be researched.
Description
Technical Field
The invention belongs to the technical field of ethylene cracking, and particularly relates to a P-graph-based modeling method for a full-flow superstructure model of ethylene, which is called PECMA (P-graph based ethylene cracking and analysis method for short).
Background
The ethylene production in China faces the problems of heavier cracking raw materials, higher energy consumption and long-term shortage of the cracking raw materials, and meanwhile, oil refining is in a state of excess capacity for a long time due to the adjustment of the economic structure in China. Therefore, the ethylene production enterprises and the oil refining production enterprises together convert the surplus oil refining capacity into the high-quality ethylene cracking light raw material, and become the first choice of the current Chinese ethylene production enterprises. In the process, the excess refinery gas in the refinery is pretreated and converted into high-quality light raw materials to improve the benefit of an ethylene production device, and meanwhile, the improvement of the yield of high value-added light products and the insufficient treatment capacity of subsequent equipment of a cracking furnace in the ethylene plant form a new contradiction, thereby greatly limiting the further improvement of economic benefit and environmental benefit. Since the analysis resolves this conflict involves individual units including cleavage, separation, etc. Therefore, it is necessary to establish an ethylene full-flow model and search the optimum degree of raw material lightening caused by refining integration under the current equipment condition based on the model.
In addition, the solution of the contradiction relates to the synchronous optimization of the structural working condition and the parameters. The method is characterized in that models of equipment possibly involved in all processes are placed in a network, logistics are represented by edges of the network, and the established model capable of synchronously performing parameter and structure optimization is called a superstructure model and is one of the best methods for synchronously solving structure and parameter optimization. As a superstructure modeling method, a series of algorithms (including a maximum structure generation algorithm (MSG), a solution structure generation algorithm (SSG) and an acceleration branch definition Algorithm (ABB)) in a P-graph method enable a superstructure model to obtain an optimal solution and a large number of near-to-suboptimal solutions through one-time optimization calculation, the model constructed based on the method has better visualization performance, and a high-quality suboptimal solution set is greatly beneficial to guiding optimization decision of a process or analyzing bottleneck equipment. A basic P-graph diagram is shown in figure 1.
In the past, superstructure modeling based on P-graph does not relate to ethylene or similar large-scale industrial devices, and a P-graph method is directly applied to superstructure modeling of a whole process of a large-scale ethylene industrial device, and some technical problems need to be solved, for example:
1) large industrial devices such as ethylene often contain a large amount of equipment, and in order to simplify the complexity of superstructure modeling, it is necessary to study how to screen equipment to be modeled for the problem to be studied and merge other equipment in the process;
2) aiming at the contradiction between the logistics circulation of large industrial devices such as ethylene and the like and the stability of a P-graph algorithm, how to perform ring-opening on the established maximum structure of the superstructure and ensure the minimum influence of the ring-opening on the process is required to be researched.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a P-graph-based ethylene full-flow superstructure model modeling method, model complexity is simplified by combining problem screening equipment and associated process equipment in a model, and contradiction between logistics circulation and P-graph algorithm stability is solved by realizing maximum structural ring-opening of the superstructure and ensuring minimum influence of ring-opening on the process.
The technical scheme is as follows:
a P-graph-based ethylene full-process superstructure model modeling method comprises the following steps:
1) establishing a model by taking a process of converting straight-chain alkane into olefin as a target object, wherein the process comprises four parts of cracking, quenching, compressing and separating, and the quenching and compressing processes are simplified into an operation unit type node;
2) the ethane circulation existing in the process is subjected to ring-opening, a ring-opening superstructure is constructed, and M is used2And M1Representing ethane, the product of the ethane cyclic cracking process being M3And M4The corresponding yields are then as follows:
in the formula, rnThe yield of recycled ethane for the ratio corresponding to the superstructure edge n is shown in equation (3):
each time the cycle number is ncThen pass through nc=NcAfter a second cycle, product M3And M4The yields of (a) are expressed by the following equations (4) and (5), respectively:
over a number of iterations, the percentage of remaining ethane is shown in equation (6):
if eta2< 0.0005, the remaining amount of ethane is approximately neglected, r'1-r1Then, the corrected variation value is:
3) analyzing the optimal and near suboptimal solution sets obtained based on the P-graph algorithm:
generating an objective function and a constraint according to the maximum structure and parameters of the PNS problem, and giving an optimal solution and a large number of suboptimal solutions based on an accelerated branch definition algorithm, wherein the established model comprises three types of constraints: (1) upper feed supply and lower product yield (ethylene) constraints; (2) upper and lower limits of the processing capacity of the equipment are restricted; (3) material balance of each node;
calculating an objective function as shown in equation (9) according to the raw material price and the product price:
the optimal parameters under each feasible structure are solved according to the objective function (9), and sequencing is carried out according to the objective function maxprofit so as to give optimal and suboptimal solutions; firstly, an index for measuring the weight of the raw material is calculated based on the optimal solution and the secondary solution, and the weight index adopts the weight ratio of the light raw material as the formula (10)LTo show that:
in the formula, FiIs a raw material RiOf (2) a traffic, coefficienti(l) Is reacted with a starting material RiCoefficient of property when starting material RiWhen the average carbon chain length l of (2) is not more than 6, coefficienti(l) 1, otherwise coefficienti(l)=0;
The influence of the refinery-ethylene integration process on the environment is characterized by the index SPI as shown in formula (11):
andis the carbonaceous percentage of the feedstock and product, FRiIs the mass flow rate, A is the carbon emission fixed value which can be treated by the forest in unit area in one year, and 2t/km is taken in the modeling process2。
Further, the fuel consumed by the cracking process is synchronously considered due to sustainability indicatorsThe fuel consumption of the cracked section is estimated based on historical data fit based on the flow rate of the feedstock and the properties of the feedstock.
Further, the catalytically dry gas in the feed needs to be removed from the acid gases that are not involved in the cracking process.
The technical effects obtained by the invention are as follows:
1. an ethylene full-Process superstructure modeling algorithm constructed based on a P-graph (Process graph) algorithm is used for carrying out simulation on an ethylene full-Process under a refining and chemical integration background, equipment required by modeling is analyzed, screened and deleted by applying domain knowledge, the problem of expression of an ethane circulation phenomenon in a superstructure model is solved based on a ring-opening technology, an optimal solution scheme and a near-suboptimal solution scheme are calculated, and the optimal light raw material ratio suitable for a current device is obtained. The optimal solution and a large number of near suboptimal solutions of the process can be obtained based on the P-graph algorithm, and decision and analysis are facilitated.
2. The established superstructure modeling comprises various working conditions of the existing equipment under the condition of load change, and the process can be synchronously subjected to structural optimization and parameter optimization based on the superstructure modeling, so that the search domain of the process optimization problem is enlarged.
3. Based on the limit concept, the ethane circulation process which is easy to generate dead circulation for superstructure calculation is subjected to ring opening by adjusting parameters, and the modeling precision of a superstructure model is ensured, and meanwhile, the complexity of the model and optimization calculation is reduced.
4. The degree of process-refining integration is described by the mass ratio of the light feedstock mass to the total cracking feedstock mass.
Drawings
FIG. 1 is a schematic diagram of a P-graph based modeling method for a full-flow superstructure model of ethylene.
Fig. 2 is a loop-opening flow diagram of the superstructure.
FIG. 3 is a diagram of a P-graph-based maximum structure model for ethylene cracking full-flow modeling.
FIG. 4 is a flow chart of a P-graph based modeling method for a full-flow superstructure model of ethylene.
Detailed Description
In order to make those skilled in the art better understand the technical solution of the present invention, the P-graph-based ethylene full-flow superstructure modeling method provided by the present invention is described in detail below with reference to the following examples. The following examples are intended to illustrate the invention only and are not intended to limit the scope of the invention.
The invention aims at carrying out modeling research on a specific ethylene plant under the integrated background of oil refining and ethylene. In the actual research process, an optimization model is established by taking economic benefit as a target, material flow and equipment working conditions of an ethylene process as decision variables and actual conditions of equipment and materials as constraints, and an optimal scheme is comprehensively determined based on economic indexes and carbon emission corresponding to suboptimal solutions so as to guide actual production. In different schemes, the optimization of the analysis process by taking the ratio of the light raw material to the total raw material (light hydrocarbon rate) as an independent variable and taking economic and environmental benefits as a dependent variable is used for exploring the optimal degree of refining integration for a specific device.
Equipment for screening and deleting models required to be modeled by using process knowledge
The goal of building an ethylene full-process superstructure model is to optimize process paths and operating parameters, and for an ethylene production process comprising a large number of devices, the devices that affect the above goal should be focused, and it is not necessary to connect all the devices into the network of the superstructure, so at the very beginning of the modeling process, the devices that are closely related to the analysis index should be screened out for modeling. On the other hand, a plurality of series or parallel devices should be simplified into one operation unit as much as possible to control the complexity of the P-graph modeling and improve the modeling efficiency.
The set analysis index mainly focuses on the influence of the economic benefit and the process of cracking on the environment, while the economic benefit index mainly focuses on the increment generated in the process of converting low-value hydrocarbon (mainly straight-chain alkane) into high-value-added olefin, so that a model of equipment (cracking furnace group) related to chemical reaction needs to be established.
The main flow of the cracking process comprises four parts of cracking, quenching, compressing and separating. The rapid cooling and compression processes are mainly physical processes, relate to energy consumption and phase change and are the key points for researching energy efficiency, but in the process of researching economy and environment, a large number of equipment models related to the two parts can be simplified into an O-type node (using OMC)iRepresentation). Under normal production conditions (i.e. conditions involved in the plant), the separation part has little influence on the economic benefits and environmental indicators (emissions) of the whole process, but because the upgrading of the co-produced raw material is not matched with the processing capacity of the existing separation part, and the mismatching can greatly limit the upgrading degree of the feasible cracked raw material, the separation part is also one of the key points of modeling. In addition, since the research involves environmental influences, such as carbon dioxide emissions, fuel systems not included in the main process should also be added to the model, taking into account in the modeling and optimization process simultaneously.
Second, the ethane circulation existing in the process is subjected to ring-opening based on the limit thought
Because in the actual cracking process, ethane generated by cracking heavy raw materials needs to be recycled as cracking raw materials, and the process containing recycling does not meet the non-recycling process axiom of 5 axioms in P-graph. For this purpose, a superstructure of the solution ring needs to be constructed. For complex P-graph networks, the ring is often very difficult to separate, and deleting an edge may mean modifying parameters of a large number of edges and nodes, but the P-graph structure of other products generated by ethane circulation in the ethylene cracking process is relatively simple, and only the P-graph structure is usedOnly the weight of the edge needs to be modified. As shown in FIG. 2, M2And M1All represent ethane, O1Represents the corresponding operation unit of the cracking process (actually, the whole cracking process, which is simplified to be one operation unit for convenience). Assuming that the product of the ethane recycle cracking process is M3And M4Their yields can be expressed by equation (1) and equation (2), respectively:
the yield of recycled ethane is shown in equation (3):
each time the cycle number is ncThen pass through nc-NcAfter a second cycle, product M3And M4The yield of (c) can be expressed by the following formula (4) and formula (5), respectively:
the percentage of remaining ethane is shown in equation (6):
if eta2< 0.0005, the amount of ethane remaining can be approximately neglected, usually r'1=r1Then, the corrected variation value is:
thirdly, analyzing the optimal and near suboptimal solution sets obtained based on the P-graph algorithm
And automatically generating an objective function and a constraint according to the maximum structure and the parameters of the PNS problem, and giving an optimal solution and a large number of suboptimal solutions based on an accelerated branch-and-bound (ABB) algorithm. The established model contains three types of constraints: (1) the upper feed supply and lower product yield (ethylene) constraints, which are determined by the feed capacity of the refinery light feedstock; (2) the upper and lower limits of the processing capacity of the equipment are usually determined by design values or data of actual operation of the plant; (3) material balance of each node.
An objective function maxprofit shown in formula (9) is automatically calculated according to the provided raw material price and product price:
and solving the optimal parameters under each feasible structure according to the objective function, and sequencing according to the objective function maxprofit to give optimal and suboptimal solutions. Firstly, an index for measuring the light weight of the raw material is calculated based on the optimal solution and the secondary solution, and the light weight index adopts the light raw material proportion rale as the formula (10)nTo show that:
in the formula, FiIs a raw material RiOf (2) a traffic, coefficienti(l) Is reacted with a starting material RiIs of a natureCoefficient of relevance when starting material RiWhen the average carbon chain length l of (2) is not more than 6, coefficienti(l) -1, otherwise coefficienti(l)-0。
The influence of the refinery-ethylene integration process on the environment is characterized by the index SPI as shown in formula (11):
andis the carbonaceous percentage of the feedstock and product, FRiIs the mass flow rate, A is the carbon emission amount which can be treated by the forest in unit area in one year, is a fixed value, and is taken as 2t/km in the modeling process2. Because sustainability index needs to synchronously consider fuel consumed by cracking process, a raw material is addedRepresenting the fuel. The fuel consumption of the cracked section is estimated based on historical data fit based on the flow rate of the feedstock and the properties of the feedstock. In addition, the catalytically dry gas in the feed requires removal of the pretreated portion, i.e., the separated acid gas, because this portion of the feed does not participate in the cracking process.
And then, carrying out multi-index analysis on the optimal and suboptimal solution sets obtained by modeling, and analyzing the limit of light weight and the influence of the light weight on the economy and emission of the whole ethylene flow by combining the visual hypergraph chart corresponding to each solution so as to provide an operation guide with operability.
Example 1
To verify the effectiveness of the process, there were 11 cracking furnaces (N) in a given refinery-ethylene co-production ethylene plant before co-productionF11), one of which is ready for use. The ethylene plant uses naphtha (NAP, R)1),Hydrogenated tail oil (HVGO, R)2) Liquefied petroleum gas (LPG, R)3) Recycle of ethane (R)4) And C5 (R)5) As cracking feedstock, hydrogen (H2, P) is produced1) Methane (CH4, P)2) Ethane (C2H6, P)3) Ethylene (C2H4, P)4) Propylene (C3H6, P)5),C4(P6) And pyrolysis gasoline and diesel (C5+, (P)7) Etc.).
The distribution of the raw materials in the whole plant before co-production is shown in Table 1, in which case the weight reduction of the raw materials is 4.16% and the SPI is 3.405X 10 by applying the formulas 9, 10 and 114year·km2The economic benefit is 51582 $/h.
TABLE 1 Pre-Cogeneration Condition
According to the calibration test of the whole ethylene flow, the working conditions of different separation systems are changed when the load is greatly changed as shown in table 2. The calibration test is carried out on the working state of the whole process under the condition of load change.
TABLE 2 variation of operating conditions of the separation system with load
In the refining and chemical integration process, light raw materials are increased. The supply of light feedstocks to refiners upstream of the ethylene production facility is shown in table 3. This supply is the maximum degree of process lightening.
TABLE 3 feedstocks for refinery-ethylene integrated post-ethylene plants
Based on the data and the equipment screening method, the main equipment for modeling is determined to be a cracking furnace group and analysis equipment, the compression part and the quenching part are combined for modeling, and material balance modeling needs to be carried out on the refinery gas pretreatment and fuel pipe network system. And then, performing ring-opening calculation based on the formulas 1 to 8 to establish a superstructure model. Some of the parameters used for the superstructure model optimization are shown in table 4.
TABLE 4 parameters associated with index calculation
The superstructure created is shown in fig. 3. The calculation of multiple indexes is performed on the obtained optimal solution and the near suboptimal solution based on the formulas 9 to 11, and the results are as follows:
the maximum profit appears as three peaks in all sub-optimal solutions. The highest peak in economic indicator corresponds to a maximum profit of $ 126677/hour, an increase of 21.79% compared to the performance before refinery-ethylene integration ($ 100722/hour). The SPI value for this solution was 47.77(km2 years), although better than that before co-production (48.02km2 years), but it was not the optimal solution under current environmental policy tightening. The second highest peak in the economic indicator corresponds to an SPI of 26.07km2 x years, which represents an economic benefit of $ 122336/hour, and the economic benefit of this solution is slightly reduced (3.43%) but the environmental protection effect is greatly improved (SPI is reduced by 45.43%) compared to the most profitable solution. In conclusion, the optimized proposal that the light hydrocarbon ratio is 11.91 percent, the economic benefit is improved by 21.47 percent, and the negative impact (SPI) on the environment is reduced by 40.27 percent is recommended.
The full-process ethylene superstructure model based on the P-graph framework has the advantage of being capable of obtaining a series of suboptimal solutions. In the modeling process, a method for screening related equipment of a complex industrial process is provided according to the characteristics of 'three-pass-one-reverse' and a problem to be researched. Based on the limit principle, a ring-opening strategy of a specific situation in the logistics circulation of the complex industrial process is given. Practical cases verify the effectiveness of the modeling method, and reveal the advantages and limitations of the light ethylene cracking raw material under the integrated oil refining-ethylene background: for the plants involved in the case, the light feedstocks are temporarily not completely converted into feedstocks for ethylene cracking, limited by the processing capacity of the separation apparatus. Under the optimal condition, the light degree (light raw material ratio) of the raw materials can reach 11.91 percent, the economic efficiency of the ethylene process can be improved by 21.47 percent, and the emission can be reduced by 40.27 percent.
The present invention is not limited to the above-described examples, and various changes can be made without departing from the spirit and scope of the present invention within the knowledge of those skilled in the art.
Claims (3)
1. A P-graph-based ethylene full-process superstructure model modeling method is characterized by comprising the following steps:
1) establishing a model by taking a process of converting straight-chain alkane into olefin as a target object, wherein the process comprises four parts of cracking, quenching, compressing and separating, and the quenching and compressing processes are simplified into an operation unit type node;
2) the ethane circulation existing in the process is subjected to ring-opening, a ring-opening superstructure is constructed, and M is used2And M1Representing ethane, the product of the ethane cyclic cracking process being M3And M4The corresponding yields are then as follows:
in the formula, rnThe yield of recycled ethane for the ratio corresponding to the superstructure edge n is shown in equation (3):
each time of the loop is numbered asncThen pass through nc=NCAfter a second cycle, product M3And M4The yields of (a) are expressed by the following equations (4) and (5), respectively:
over a number of iterations, the percentage of remaining ethane is shown in equation (6):
if eta2< 0.0005, the remaining amount of ethane is approximately neglected, r'1=T1Then, the corrected variation value is:
3) analyzing the optimal and near suboptimal solution sets obtained based on the P-graph algorithm:
generating an objective function and a constraint according to the maximum structure and parameters of the PNS problem, and giving an optimal solution and a large number of suboptimal solutions based on an accelerated branch definition algorithm, wherein the established model comprises three types of constraints: (1) upper feed supply and lower product yield (ethylene) constraints; (2) upper and lower limits of the processing capacity of the equipment are restricted; (3) material balance of each node;
calculating an objective function as shown in equation (9) according to the raw material price and the product price:
the optimal parameters under each feasible structure are solved according to the objective function (9), and sequencing is carried out according to the objective function maxprofit so as to give optimal and suboptimal solutions; firstly, an index for measuring the weight of the raw material is calculated based on the optimal solution and the secondary solution, and the weight index adopts the weight ratio of the light raw material as the formula (10)LTo show that:
in the formula, FiIs a raw material RiOf (2) a traffic, coefficienti(l) Is reacted with a starting material RiCoefficient of property when starting material RiWhen the average carbon chain length l of (2) is not more than 6, coefficienti(l) 1, otherwise coefficienti(l)=0;
The influence of the refinery-ethylene integration process on the environment is characterized by the index SPI as shown in formula (11):
2. The P-graph based ethylene of claim 1The modeling method of the full-flow superstructure model is characterized in that the fuel consumed by the cracking process is synchronously considered due to the sustainability indexThe fuel consumption of the cracked section is estimated based on historical data fit based on the flow rate of the feedstock and the properties of the feedstock.
3. The modeling method for the full-flow superstructure model of ethylene based on P-graph according to claim 1, characterized in that the catalytic dry gas in the feed needs to remove the acid gas not participating in the cracking process.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010557272.7A CN112201309B (en) | 2020-06-18 | 2020-06-18 | Ethylene full-flow super-structure model modeling method based on P-graph |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010557272.7A CN112201309B (en) | 2020-06-18 | 2020-06-18 | Ethylene full-flow super-structure model modeling method based on P-graph |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112201309A true CN112201309A (en) | 2021-01-08 |
CN112201309B CN112201309B (en) | 2024-02-02 |
Family
ID=74006114
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010557272.7A Active CN112201309B (en) | 2020-06-18 | 2020-06-18 | Ethylene full-flow super-structure model modeling method based on P-graph |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112201309B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040117040A1 (en) * | 2002-12-09 | 2004-06-17 | Bijan Sayyarrodsari | System and method of adaptive control of processes with varying dynamics |
CN108073784A (en) * | 2017-12-11 | 2018-05-25 | 北京谱豪胜世纪工程技术有限公司 | A kind of method for solving of the optimization of ethylene feedstock scheduling scheme of Kernel-based methods figure |
CN108535990A (en) * | 2017-03-03 | 2018-09-14 | 夏普株式会社 | Image forming apparatus and image forming method |
CN109902925A (en) * | 2019-01-22 | 2019-06-18 | 北京化工大学 | Flow scheduling modeling method based on carbon emission folder point analysis and procedure chart |
-
2020
- 2020-06-18 CN CN202010557272.7A patent/CN112201309B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040117040A1 (en) * | 2002-12-09 | 2004-06-17 | Bijan Sayyarrodsari | System and method of adaptive control of processes with varying dynamics |
CN108535990A (en) * | 2017-03-03 | 2018-09-14 | 夏普株式会社 | Image forming apparatus and image forming method |
CN108073784A (en) * | 2017-12-11 | 2018-05-25 | 北京谱豪胜世纪工程技术有限公司 | A kind of method for solving of the optimization of ethylene feedstock scheduling scheme of Kernel-based methods figure |
CN109902925A (en) * | 2019-01-22 | 2019-06-18 | 北京化工大学 | Flow scheduling modeling method based on carbon emission folder point analysis and procedure chart |
Non-Patent Citations (3)
Title |
---|
牟鹏;顾祥柏;朱群雄: "基于P-graph的乙烯裂解原料调度建模与优化", 化工学报, no. 002 * |
耿志强;朱群雄;顾祥柏;林晓勇;: "基于多群竞争PSO-RBFNN的乙烯裂解深度智能优化控制", 化工学报, no. 08 * |
郝泽;张毅;姜晓滨;贺高红;肖武;: "P-图理论优化炼厂气轻烃回收过程工艺参数", 计算机与应用化学, no. 01 * |
Also Published As
Publication number | Publication date |
---|---|
CN112201309B (en) | 2024-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Investigation and optimization analysis on deployment of China coal chemical industry under carbon emission constraints | |
Wang et al. | Sustainable design and synthesis of hydrocarbon biorefinery via gasification pathway: Integrated life cycle assessment and technoeconomic analysis with multiobjective superstructure optimization | |
Zhou et al. | Produce petrochemicals directly from crude oil catalytic cracking, a techno-economic analysis and life cycle society-environment assessment | |
CN109698014B (en) | Method for predicting catalytic reformate | |
CN106444672A (en) | Molecular-level real time optimization (RTO) method for oil refining and petrochemical device | |
CN104789256A (en) | Real-time yield predicting method for catalytic cracking device | |
Jia et al. | Paraffin-based crude oil refining process unit-level energy consumption and CO2 emissions in China | |
CN103524284A (en) | Forecasting and optimizing method for ethylene cracking material configuration | |
CN103087750B (en) | Naphtha industrial cracking furnace value maximization model construction method | |
Wu et al. | Design and optimization of bio-oil co-processing with vacuum gas oil in a refinery | |
Zhou et al. | Crude oil hierarchical catalytic cracking for maximizing chemicals production: Pilot-scale test, process optimization strategy, techno-economic-society-environment assessment | |
CN111475957B (en) | Oil refining process production plan optimization method based on device mechanism | |
CN102768702A (en) | Oil refining production process schedule optimization modeling method on basis of integrated control optimization | |
Bi et al. | Knowledge expression, numerical modeling and optimization application of ethylene thermal cracking: From the perspective of intelligent manufacturing | |
CN104765347A (en) | Yield real-time prediction method in residual oil delayed coking process | |
Han et al. | Energy consumption hierarchical analysis based on interpretative structural model for ethylene production | |
Zhou et al. | Opportunities for utilizing waste cooking oil in crude to petrochemical process: Novel process design, optimal strategy, techno-economic analysis and life cycle society-environment assessment | |
Zhang et al. | Prediction of research octane number loss and sulfur content in gasoline refining using machine learning | |
Gong et al. | Multi-level and multi-granularity energy efficiency diagnosis scheme for ethylene production process | |
CN106597853A (en) | Active dynamic regulation method in hydrocracking process | |
CN109902925A (en) | Flow scheduling modeling method based on carbon emission folder point analysis and procedure chart | |
CN112201309A (en) | P-graph-based ethylene full-process superstructure model modeling method | |
CN103087753B (en) | Heavy raw material industrial cracking furnace value maximization model construction method | |
WO2023143106A1 (en) | Multi-objective optimization method and apparatus for fluid catalytic cracking process | |
Mahmoud et al. | Systematic methodology for optimal enterprise network design between bio-refinery and petroleum refinery for the production of transportation fuels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |