CN107291975A - A kind of method and system of catalytic cracking reaction product hard measurement - Google Patents

A kind of method and system of catalytic cracking reaction product hard measurement Download PDF

Info

Publication number
CN107291975A
CN107291975A CN201710305383.7A CN201710305383A CN107291975A CN 107291975 A CN107291975 A CN 107291975A CN 201710305383 A CN201710305383 A CN 201710305383A CN 107291975 A CN107291975 A CN 107291975A
Authority
CN
China
Prior art keywords
catalytic cracking
cracking reaction
model
optimized
parameter
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.)
Pending
Application number
CN201710305383.7A
Other languages
Chinese (zh)
Inventor
吴迎亚
曹道帆
蓝兴英
高金森
徐春明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201710305383.7A priority Critical patent/CN107291975A/en
Publication of CN107291975A publication Critical patent/CN107291975A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The embodiments of the invention provide a kind of method and system of catalytic cracking reaction product hard measurement, this method includes:Gather the actual measurement Product distribution data and duty parameter of the different operating modes of catalytic cracking reaction;According to duty parameter, numerical simulation, the corresponding calculating Product distribution data of the different operating modes of acquisition are carried out to catalytic cracking reaction kinetic model using Fluid Mechanics Computation method;If the deviation summation for calculating Product distribution data and actual measurement Product distribution data is more than predetermined threshold value, the model parameter to catalytic cracking reaction kinetic model is optimized, until deviation summation is less than or equal to predetermined threshold value, obtains Optimized model parameter;According to Optimized model parameter, using Fluid Mechanics Computation method, the corresponding optimized product distributed data of the different operating modes of acquisition completes catalytic cracking reaction product hard measurement.The system is used to perform this method.The catalytic cracking production hard measurement of different operating modes is realized in the embodiment of the present invention, the precision and adaptability of catalytic cracking production hard measurement is improved.

Description

A kind of method and system of catalytic cracking reaction product hard measurement
Technical field
The present embodiments relate to technical field of petrochemical industry, and in particular to a kind of catalytic cracking reaction product hard measurement Method and system.
Background technology
Catalytic cracking unit (FCCU, Fluid Catalytic Cracking Unit) is the key device of Petrochemical Enterprises, It is heavy oil deep processing, the important means of lighting, aggravates with heaviness, the in poor quality of crude oil, catalytic cracking unit Effect is more protruded, and very important status is occupied in petroleum refining industry.Heavy crude by catalytic cracking unit by lighting, Dry gas (C1, C2, H2, H2S), liquefied gas (C3, C4), gasoline (C5~C11), diesel oil (C10~C20), slurry oil are generated (with condensed ring Based on aromatic hydrocarbons), coke (condensation product) etc..It is limited to the factors such as plant bulk, e measurement technology and catalytic cracking reaction Product yield in the feature that the Complex Flows reaction system of itself is highly coupled, catalytic cracking reaction product is can not to use instrument Table direct measurement is obtained, therefore using online measurable auxiliary variable, such as temperature, pressure, flow are estimated to be difficult to examine in real time The soft-measuring technique of the product yield of survey is widely applied in catalytic cracking unit.
In the prior art, catalytic cracking soft-measuring technique uses lumped reaction kinetics mostly, passes through lumping kinetics mould Type, obtains the product yield of catalytic cracking reaction.But this traditional lumped watershed hydrologic model model is in modeling process, only Focus on cracking reaction dynamic process, it is assumed that the state of reacting gas is isothermal piston flow, and be have ignored in catalytic cracking system Flowing, heat transfer, the influence of mass transfer and turbulence pulsation to reaction so that numerical simulation has certain deviation, to catalytic cracking The precision of reaction hard measurement result affects.Moreover, the model parameter adaptability obtained under conventional method is poor, it is past Toward being result for being obtained under special reaction condition and operating mode, therefore when the factors such as reaction condition, catalyst property change, Need to carry out model correction and compensation by collection in worksite data.
Therefore, how to propose a kind of scheme, it is possible to increase the adaptability of catalytic cracking reaction flexible measurement method, further carry The precision of high measurement result, as urgent problem to be solved.
The content of the invention
For defect of the prior art, the embodiments of the invention provide and a kind of catalytic cracking reaction product hard measurement Method and system.
On the one hand, the embodiments of the invention provide a kind of choosing method of catalytic cracking reaction hard measurement auxiliary variable, bag Include:
The corresponding floor data of the different operating modes of catalytic cracking reaction is gathered, the floor data includes catalytic cracking reaction Survey Product distribution data and duty parameter;
According to the duty parameter, the catalytic cracking reaction kinetic model set up using Fluid Mechanics Computation method to pre-selection Numerical simulation is carried out, the corresponding calculating Product distribution data of the different operating modes of catalytic cracking reaction is obtained;
If judging, the deviation summation for knowing the calculating Product distribution data and the actual measurement Product distribution data is more than in advance If threshold value, then the model parameter to the catalytic cracking reaction kinetic model is optimized, until the deviation summation is less than Or equal to the predetermined threshold value, obtain Optimized model parameter;
According to the Optimized model parameter and the catalytic cracking reaction kinetic model, the Fluid Mechanics Computation is utilized Method, obtains the corresponding optimized product distributed data of the different operating mode of catalytic cracking reaction, using the optimized product distributed data as The result of the catalytic cracking reaction product hard measurement.
Further, it is described that the catalytic cracking reaction kinetic model that pre-selection is set up is carried out using Fluid Mechanics Computation method Numerical simulation, obtains the corresponding calculating Product distribution data of the different operating modes of catalytic cracking reaction, including:
If the numerical simulation is first simulation, model parameter initial value is set, and according to the corresponding institute of different operating modes Duty parameter is stated, parameter of a flow model is set;
According to the model parameter initial value and the parameter of a flow model, pre-selection is set up using Fluid Mechanics Computation method Catalytic cracking reaction kinetic model carry out numerical simulation, obtain that different operating modes described in catalytic cracking reaction are corresponding to calculate production Thing distributed data.
Further, if it is described judgement know it is described calculating Product distribution data and the actual measurement Product distribution data it is inclined Poor summation is more than predetermined threshold value, then the model parameter to the catalytic cracking reaction kinetic model is optimized, including:
If judgement knows that the deviation summation is more than predetermined threshold value, the model parameter is carried out using genetic algorithm excellent Change.
Further, it is described according to the Optimized model parameter and the catalytic cracking reaction kinetic model, utilize institute Fluid Mechanics Computation method is stated, the optimized product distributed data of catalytic cracking reaction is obtained, including:
The Optimized model parameter is substituted into the catalytic cracking reaction kinetic model, the Fluid Mechanics Computation is utilized Method, numerical simulation is carried out to the catalytic cracking reaction kinetic model, obtains the optimized product distributed data.
Further, the catalytic cracking reaction kinetic model is ten four-lump kinetics models.
Further, the duty parameter includes:Temperature data, pressure data and data on flows.
Further, the actual measurement Product distribution data includes:Dry gas yield, yield of liquefied gas, yield of gasoline, diesel oil are received Rate and coking yield.
On the other hand, the embodiment of the present invention provides a kind of system of catalytic cracking reaction product hard measurement, including:Data are adopted Collect unit, for gathering the corresponding floor data of the different operating modes of catalytic cracking reaction, the floor data includes Catalytic Cracking Unit of Measure The actual measurement Product distribution data and duty parameter answered;
Product distribution unit is calculated, for according to the duty parameter, being set up using Fluid Mechanics Computation method to pre-selection Catalytic cracking reaction kinetic model carries out numerical simulation, obtains the corresponding product that calculates of the different operating modes of catalytic cracking reaction and is distributed Data;
Model Parameter Optimization unit, if for judging to know the calculating Product distribution data and actual measurement product distribution The deviation summation of data is more than predetermined threshold value, then the model parameter to the catalytic cracking reaction kinetic model is optimized, Until the deviation summation is less than or equal to the predetermined threshold value, Optimized model parameter is obtained;
Hard measurement unit, for according to the Optimized model parameter and the catalytic cracking reaction kinetic model, utilizing The Fluid Mechanics Computation method, obtains the corresponding optimized product distributed data of the different operating modes of catalytic cracking reaction, by the optimization Product distribution data as the catalytic cracking reaction product hard measurement result.
Another further aspect, the embodiment of the present invention provides a kind of electronics chosen for catalytic cracking reaction hard measurement auxiliary variable Equipment, including:
At least one processor;And
At least one memory being connected with the processor communication, wherein:
The memory storage, which has, to call described program to refer to by the programmed instruction of the computing device, the processor Order is able to carry out above-mentioned method.
Another aspect, the embodiment of the present invention provides a kind of non-transient computer readable storage medium storing program for executing, the non-transient calculating Machine readable storage medium storing program for executing stores computer instruction, and the computer instruction makes the computer perform the above method.
The method and system of the catalytic cracking reaction product hard measurement that example is provided in real time of the invention, utilizes Fluid Mechanics Computation Method CFD carries out the numerical simulation of catalytic cracking reaction, and the model parameter of catalytic cracking reaction kinetic model is optimized, The model parameter of the catalytic cracking reaction kinetic model for the condition that meets is obtained, Fluid Mechanics Computation method CFD is re-used and carries out The numerical simulation of catalytic cracking reaction, obtains the optimized product distributed data of final acquisition catalytic cracking reaction, realizes catalysis Reaction product hard measurement.Catalytic cracking reaction kinetic model is combined with Fluid Mechanics Computation method CFD, by optimizing power Model parameter is learned, catalytic cracking reaction product hard measurement is realized, improves the precision of catalytic cracking reaction product hard measurement result. In addition, the embodiment of the present invention is when catalytic cracking operating mode changes, it is only necessary to corresponding adjustment catalytic cracking reaction dynamics The model parameter of model, you can the corresponding calculating Product distribution data of the different operating modes of acquisition, realizes the catalytic cracking of different operating modes Reaction product hard measurement, improves the adaptability of catalytic cracking reaction flexible measurement method.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are this hairs Some bright embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, can be with root Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the schematic flow sheet of catalytic cracking reaction product flexible measurement method in the embodiment of the present invention;
Fig. 2 is the reaction network schematic diagram of ten four-lump kinetics models in the embodiment of the present invention;
Fig. 3 is the schematic flow sheet of another catalytic cracking reaction product flexible measurement method in the embodiment of the present invention;
Fig. 4 shows for the flow for carrying out catalytic cracking reaction numerical simulation using Fluid Mechanics Computation method in the embodiment of the present invention It is intended to;
Fig. 5 is the schematic flow sheet of the method for another catalytic cracking reaction product hard measurement in the embodiment of the present invention;
Fig. 6 be the embodiment of the present invention in carry out Model Parameter Optimization when fitness function curve synoptic diagram;
Fig. 7 is the structural representation of catalytic cracking reaction product hard measurement system in the embodiment of the present invention;
Fig. 8 be the embodiment of the present invention in be used for catalytic cracking reaction product hard measurement electronic equipment structural representation.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is A part of embodiment of the present invention, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art The every other embodiment obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
Fig. 1 is the schematic flow sheet of catalytic cracking reaction product flexible measurement method in the embodiment of the present invention, as shown in figure 1, The method of catalytic cracking reaction product hard measurement provided in an embodiment of the present invention includes:
S1, the corresponding floor data of the different operating modes of collection catalytic cracking reaction, the floor data include Catalytic Cracking Unit of Measure The actual measurement Product distribution data and duty parameter answered;
Specifically, the corresponding floor data of the different operating modes of catalytic cracking reaction is gathered as needed, and floor data includes urging Change the actual measurement Product distribution data and duty parameter of cracking reaction.Wherein duty parameter includes:Temperature data, pressure data and stream Data are measured, actual measurement Product distribution data includes:Dry gas yield, yield of liquefied gas, yield of gasoline, diesel yield and coking yield. Duty parameter can gather the number for obtaining real-time working condition parameter or being preserved from plant historian by industrial site According to middle acquisition history duty parameter.The duty parameter collected is designated as data matrix and is designated as C, data matrix by the embodiment of the present invention In include but is not limited to:1) temperature data:Gas outlet temperature, fresh feed temperature, freshening oil temperature, reprocessed oil slurry temperature, Atomizing steam temperature, pre-lift vapor (steam) temperature and catalyst temperature;2) pressure data:Reactor outlet pressure, reactor pressure Drop;3) data on flows:Fresh feed flow, recycle oil flow, reprocessed oil slurry flow, the mass flow of catalyst, pre-lift are steamed Vapour amount, atomizing steam amount.It is the real-time actual measurement Product distribution data or history number that Site Detection is obtained to survey Product distribution data Product distribution data is surveyed according to the history stored in storehouse.Actual measurement product point in real time can be obtained by industrial site collection History actual measurement Product distribution data is obtained in cloth data or the data preserved from plant historian, the embodiment of the present invention can P is designated as so that actual measurement Product distribution data is designated as into data matrix.
Wherein floor data can be obtained by the monitoring system or catalytic cracking DCS system of catalytic cracking unit, or logical Other modes acquisition is crossed, the embodiment of the present invention is not especially limited.
S2, according to the duty parameter, the catalytic cracking reaction dynamics for utilizing Fluid Mechanics Computation method to set up pre-selection Model carries out numerical simulation, obtains the corresponding calculating Product distribution data of the different operating modes of catalytic cracking reaction;
Specifically, get after the floor data of catalytic reaction, Fluid Mechanics Computation method is utilized according to these floor datas CFD (Computational fluid dynamics) carries out Numerical-Mode to the catalytic cracking reaction kinetic model that pre-selection is set up Intend, obtain the corresponding calculating Product distribution data of the different operating modes of catalytic cracking reaction.I.e. according to the corresponding operating mode number of different operating modes According to the corresponding calculating Product distribution data of acquisition.
If S3, judgement know that the deviation summation of the calculating Product distribution data and the actual measurement Product distribution data is more than Predetermined threshold value, then the model parameter to the catalytic cracking reaction kinetic model is optimized, until the deviation summation is small In or equal to the predetermined threshold value, Optimized model parameter is obtained;
Specifically, the calculating Product distribution data got and the actual measurement Product distribution data collected in advance are made poor, Judge whether calculate Product distribution data and the deviation summation of actual measurement Product distribution data is more than predetermined threshold value, to urging if being more than The model parameter for changing cracking reaction kinetic model is optimized, until the deviation summation is less than or equal to the default threshold Value, obtains Optimized model parameter;If the deviation summation for calculating Product distribution data and actual measurement Product distribution data is less than or equal to Predetermined threshold value, then regard the model parameter of catalytic cracking reaction kinetic model now as Optimized model parameter.
S4, according to the Optimized model parameter and the catalytic cracking reaction kinetic model, utilize the calculating fluid Mechanical Method, obtains the corresponding optimized product distributed data of the different operating modes of catalytic cracking reaction, by the optimized product distributed data It is used as the result of the catalytic cracking reaction product hard measurement.
Specifically, according to Optimized model parameter and catalytic cracking reaction kinetic model, Fluid Mechanics Computation method is utilized CFD, obtains the corresponding optimized product distributed data of the different operating modes of catalytic cracking reaction, regard optimized product distributed data as catalysis The result of cracked reaction product hard measurement, realizes the product hard measurement under the different operating modes of catalytic cracking reaction.
It should be noted that when the operating mode of catalytic cracking reaction changes, floor data is adjusted accordingly, and urge Change the model parameter of cracking reaction kinetic model, re-execute step S2-S4, obtain corresponding catalysis under new working condition The optimized product distributed data of cracking reaction, realizes the catalytic cracking reaction product hard measurement of different operating modes.
The method of the catalytic cracking reaction product hard measurement that example is provided in real time of the invention, utilizes Fluid Mechanics Computation method CFD The numerical simulation of catalytic cracking reaction is carried out, and the model parameter of catalytic cracking reaction kinetic model is optimized, is obtained The model parameter of the catalytic cracking reaction kinetic model of condition is met, Fluid Mechanics Computation method CFD is re-used and is catalyzed The numerical simulation of cracking reaction, obtains the optimized product distributed data of final acquisition catalytic cracking reaction, realizes catalytic reaction Product hard measurement.Catalytic cracking reaction kinetic model is combined with Fluid Mechanics Computation method CFD, by optimizing kinetic simulation Shape parameter, realizes catalytic cracking reaction product hard measurement, improves the precision of catalytic cracking reaction product hard measurement result.This Outside, the embodiment of the present invention is when catalytic cracking operating mode changes, it is only necessary to corresponding adjustment catalytic cracking reaction kinetic simulation The model parameter of type, you can the corresponding calculating Product distribution data of the different operating modes of acquisition, realizes the Catalytic Cracking Unit of Measure of different operating modes Product hard measurement is answered, the adaptability of catalytic cracking reaction flexible measurement method is improved.
On the basis of above-described embodiment, the utilization Fluid Mechanics Computation method is moved to the catalytic cracking reaction that pre-selection is set up Mechanical model carries out numerical simulation, obtains the corresponding calculating Product distribution data of the different operating modes of catalytic cracking reaction, including:
If the numerical simulation is first simulation, model parameter initial value is set, and according to the corresponding institute of different operating modes Duty parameter is stated, parameter of a flow model is set;
According to the model parameter initial value and the parameter of a flow model, pre-selection is set up using Fluid Mechanics Computation method Catalytic cracking reaction kinetic model carry out numerical simulation, obtain that different operating modes described in catalytic cracking reaction are corresponding to calculate production Thing distributed data.
Specifically, when carrying out catalytic cracking reaction numerical simulation for the first time, it is necessary to set catalytic cracking kinetic model Model parameter initial value, specifically can rule of thumb be configured, or be set using other method, and the embodiment of the present invention is not made to have Body is limited.Set after the model parameter initial value of catalytic cracking kinetic model, according to the corresponding duty parameter of different operating modes, Set parameter of a flow model.According to model parameter initial value and parameter of a flow model, using Fluid Mechanics Computation method CFD to pre-selection The catalytic cracking reaction kinetic model of foundation carries out numerical simulation, obtains corresponding calculate of the different operating modes of catalytic cracking reaction and produces Thing distributed data.
Wherein, catalytic cracking reaction kinetic model, which belongs in reaction model, actual industrial reactor, also has complexity Gas-particle Flows, heat transfer, mass transfer, description Gas-particle Flows, heat transfer equation be called flow model (because heat transfer can be regarded as heat Amount flowing, mass transfer are considered as component flowing, and flow model or TRANSFER MODEL are typically referred to as in chemical industry).Flowing can shadow Reaction is rung, reaction can also influence flowing in turn, and flow model and reaction model are coupled together, that is, utilized by the embodiment of the present invention Fluid Mechanics Computation method CFD carries out numerical simulation to the catalytic cracking reaction kinetic model that pre-selection is set up, so that what is obtained urges Change the corresponding Product distribution data that calculates of the different operating modes of cracking reaction more accurate, further improve catalytic cracking reaction product soft The accuracy of the measurement result of measurement.
Wherein, according to model parameter initial value and parameter of a flow model, pre-selection is set up using Fluid Mechanics Computation method CFD Catalytic cracking reaction kinetic model carry out numerical simulation process it is as follows:
1) 3-D geometric model of catalytic cracking riser is set up;
2) mesh generation is carried out to 3-D geometric model in 1);
3) each entry condition, boundary condition and flow field are set using the duty parameter of the catalytic cracking reaction of pre-selection collection Primary condition, sets parameter of a flow model;
4) according to the model parameter initial value of catalytic cracking reaction kinetic model, reaction model is set;
5) catalytic cracking riser is solved using solver, obtains calculating Product distribution data.
On the basis of above-described embodiment, if the calculating Product distribution data and the actual measurement product are known in the judgement The deviation summation of distributed data is more than predetermined threshold value, then the model parameter progress to the catalytic cracking reaction kinetic model is excellent Change, including:
If judgement knows that the deviation summation is more than predetermined threshold value, the model parameter is carried out using genetic algorithm excellent Change.
Specifically, get after the calculating Product distribution data of catalytic cracking reaction, judge calculate Product distribution data and Whether the deviation summation for surveying Product distribution data is more than predetermined threshold value.If being more than, using genetic algorithm to catalysis now The model parameter of cracking reaction kinetic model is optimized, until calculating Product distribution data and surveying Product distribution data Deviation summation is less than or equal to predetermined threshold value, obtains Optimized model parameter.The Optimized model parameter is substituted into the catalysis to split Change reaction Kinetics Model, using the Fluid Mechanics Computation method, numerical value is carried out to the catalytic cracking reaction kinetic model Simulation, obtains the optimized product distributed data.The Optimized model model parameter obtained according to optimization, resets catalysis and splits Change the reaction model of reaction, numerical simulation is carried out to catalytic cracking reaction kinetic model using Fluid Mechanics Computation method CFD, obtained Obtain the optimized product distributed data of catalytic cracking reaction.
On the basis of above-described embodiment, the catalytic cracking reaction kinetic model is ten four-lump kinetics models.
Specifically, the catalytic cracking reaction kinetic model that the embodiment of the present invention is used for ten four-lump kinetics models, It is a kind of reaction Kinetics Model simplified with lumping strategy, and reaction system is divided into 14 lumps, specific partiting step For:Feedstock oil is divided into diesel oil layer (l), recycle oil-wax oil layer (m) and decompression residuum-slurry oil layer according to its flow scope first (h) three layers, each layer is constituted according to structural group is divided into alkyl carbon (P), cycloalkanes carbon (N) and aromatic carbon (A) again.Due to decompression residuum With there is polycyclic aromatic hydrocarbon in slurry oil, further collected by increasing new polycyclic aromatic hydrocarbon (Fah) in decompression residuum and slurry oil Always.Finally, it is considered to coke (CK), gasoline (GO), natural gas (LPG) and dry gas (DG), 14 lumps for obtaining feedstock oil are drawn Point.
Fig. 2 is the reaction network schematic diagram of ten four-lump kinetics models in the embodiment of the present invention, as shown in Fig. 2 14 Lumped reaction kinetics are related to 49 reactions altogether.Ten four-lump kinetics models are related to 49 kinetic parameters undetermined, i.e. frequency Factor K, K is the vector that dimension is 49, the chemical reaction rate expression of ten each gaseous components of four-lump kinetics model (lump) Formula is equation below (1):
In formula:A0--- represent catalyst activity correction factor;
--- represent the effective active of catalyst, KhAttach most importance to aromatic ring absorption constant;
Φ(Cc)=(1+0.51Cc)-2.78--- represent the influence factor of catalyst coking;
Cc--- represent catalyst coking yield percentage;
kj--- lump j reaction rate constant is represented,k0Represent pre-exponential factor;
E --- represent activation energy;
εg--- represent voidage;
ραj--- represent lump j molar concentration.
Fig. 3 is the schematic flow sheet of another catalytic cracking reaction product flexible measurement method in the embodiment of the present invention, such as Fig. 3 institutes Show, with reference to certain refinery catalytic cracking device riser reactor actual production data, specifically introduce the embodiment of the present invention Technical scheme:
The catalytically cracked material property of the operating mode 1 of table 101
Project Unit Numerical value
Oil density [kg/m3] 914.8
Viscosity under different temperatures mm2/s 33.57(50℃)/6.493(100℃)
Basic n content [ppmwt] 728
Total nitrogen content [ppmwt] 2776
Sulfur content [%] 0.5412
Conradson carbon residue [%] 3.86
Content of vanadium [ppmwt] 5
Nickel content [ppmwt] 15
Sodium content [ppmwt] 10
Iron content [ppmwt] 10
Copper content [ppmwt] 0.3
The equilibrium catalyst property of the operating mode 1 of table 102
Project Data
Activity 63
Compare surface 118
Pore volume 0.3
Heap ratio 0.85
Carbon 0.056
0-20μm 0
0-40μm 8.9
0-80μm 54.8
0-105μm 75.3
0-149μm 92.7
APS 75.4
Iron 5410
Nickel 3760
Vanadium 1800
Nickel vanadium ratio 2.09
Sodium 3800
Antimony 1960
Antimony nickel ratio 0.52
Calcium 6150
The riser operation condition of the operating mode 1 of table 103
Project Unit Numerical value
Feedstock oil treating capacity [t/h] 152
Atomizing steam mass flowrate [t/h] 9.6
Lift medium volume/mass flowrate [NM3/h] Dry gas/4574
Rich gas mass flowrate [NM3/h] 41472
Flue gas cooling water speed [kg/h] Smokeless cooling water
Temperature of charge [C] 220
Atomizing steam temperature [C] 290
Outlet temperature of riser [C] 520±2
Lift medium temperature [C] 37-38
Stripping gas temperature [C] 400
Rich gas temperature [C] 35.5
Air-blaster delivery temperature [C] 170
Environment temperature [C] 20
Material pressure [kPa] 800
Atomizing steam pressure [kPa] 1000
Lift pressure medium [kPa] 1000
Strip atmospheric pressure [kPa] 1000
Rich gas pressure [kPa] 1090
Environmental pressure [kPa] 101
Reactor pressure [kPa] 170
Regenerator pressure [kPa] 200
Dense bed bulk density [kg/m3] 480
Material mixing temperature [C] 551.4
Material preheating temperature [C] 220
Riser mixing temperature [C] 550
Outlet temperature of riser [C] 520
Regenerator dense bed temperature [C] 680
Regenerator rotation point temperature [C] 680、705
Flue-gas temperature [C] 647
Catalyst circulation rate [t/h] 1156
Coke content in regenerated catalyst [%] 0.08
Oil ratio 7-8
Catalyst inventory [kg] 180000
The supplement speed of fresh catalyst [t/h] 0.3
Recycle ratio 18~20%
T1, gathered data.The creation data that is collected from catalytic cracking unit production scene, laboratory assay data, i.e., Floor data is gathered, mainly includes duty parameter C and actual measurement Product distribution data P, 1 time parameter of certain operating mode such as table 101~103 It is shown.
T2, the ten four-lump kinetics models for setting up catalytic cracking reaction.The specific foundation of ten four-lump kinetics models Method is consistent with above-described embodiment, and here is omitted.
T3, acquisition calculate Product distribution data.The Three-dimensional simulation model of catalytic cracking reaction is set up, by step T2 The ten four-lump kinetics models set up substitute into three-dimensional numerical value numerical simulator, are carried out using computational fluid dynamics Calculate, obtain calculating Product distribution data Y.According to Three-dimensional simulation model, counted using computational fluid dynamics The method for calculating Product distribution data Y is as follows:
1) 3-D geometric model of catalytic cracking riser is set up;
2) mesh generation is carried out to the 3-D geometric model of catalytic cracking riser;
3) each entry condition, boundary condition and flow field primary condition are set using the duty parameter C gathered in step T1, Set parameter of a flow model;
4) using the ten four-lump kinetics models and model parameter K in step T2, reaction model is set.For the model Parameter K, if first calculate, is calculated using model parameter initial value, if not first calculate, is then used after optimizing in step T4 Optimized model parameter calculated;
5) riser is solved using solver, obtains calculating Product distribution data Y.
T4, setting fitness function.By the actual measurement Product distribution data P's calculated in Product distribution data Y and step T1 Deviation summation is defined as fitness function, if more than predetermined threshold value TOL, being optimized using genetic algorithm to model parameter K, Obtain Optimized model parameter.
T5, model parameter optimization.The Three-dimensional simulation mould that Optimized model parameter after optimization is substituted into step T3 Type, repeats T3-T4, until fitness function is less than or equal to predetermined threshold value TOL.Solved for computational fluid dynamics Detailed process not in the range of the narration of the embodiment of the present invention, with business software, open source software or oneself can write code Calculated.Fig. 4 is the flow in the embodiment of the present invention using Fluid Mechanics Computation method progress catalytic cracking reaction numerical simulation Fortran language can be used to write code calculating in schematic diagram, the embodiment of the present invention.
The contrast (wt%) of riser reactor numerical simulation calculation result and commercial plant data under the original operating mode of table 2
Composition Dry gas Liquefied gas Gasoline Diesel oil Coke
Industry value 3.75 16.36 41.65 29.03 7.53
The analogue value 3.60 16.13 40.52 28.22 7.55
Frequency factor pair is answered in the reaction network of table 3 pre-exponential factor and activation energy
Result of calculation and analog result after the operating mode of table 4 changes
Product Dry gas Liquefied gas Gasoline Diesel oil Coke
Industrial actual value 4.01 18.42 37.75 30.23 7.78
The GA-CFD analogues value 3.78 18.43 37.74 30.23 7.30
Relative error 6.08% 0.05% 0.03% 0.00% 6.50%
T6, hard measurement result are obtained.When operating mode does not change, calculated and produced with step T3 Three-dimensional simulation model Thing distributed data Y is used as hard measurement result.When operating mode changes, invocation step T4 to step T5 is carried out to model parameter Adapt to adjustment.Fig. 5 is the schematic flow sheet of the method for another catalytic cracking reaction product hard measurement in the embodiment of the present invention, is such as schemed Shown in 5, when working conditions change, Fig. 1 dotted portions will be had an effect, and model parameter is adjusted automatically.When operation operating mode hair During changing, the model after being corrected in a short time adapts to operation working conditions change.The embodiment of the present invention is counted by coupling Ten four-lump kinetics of fluid operator mechanics CFD approach and the three-dimensional flow-reaction coupled model of genetic algorithm to catalytic cracking reaction The model parameter of model is corrected, and hard measurement result is had good precision and adaptability.
Under operating mode 1, as shown in table 2, when working conditions change, model parameter is adjusted result of calculation automatically.Fig. 6 is this hair The curve synoptic diagram of fitness function during Model Parameter Optimization is carried out in bright embodiment, as shown in fig. 6, after operating mode changes, utilizing Genetic algorithm carries out the adjustment and optimization of model parameter, and after 400 iteration, fitness curve tends to be steady, and now 49 Individual model parameter is optimal solution, that is, obtains Optimized model parameter.
By in the flow-reaction coupled model for optimizing the Optimized model parameter substitution catalytic cracking riser reactor obtained, pass through Calculate and obtain product yield at leg outlet, and by calculated value and industrial actual comparison, such as table 4.
The method of the catalytic cracking reaction product hard measurement that example is provided in real time of the invention, utilizes Fluid Mechanics Computation method CFD The numerical simulation of catalytic cracking reaction is carried out, and the model parameter of catalytic cracking reaction kinetic model is optimized, is obtained The model parameter of the catalytic cracking reaction kinetic model of condition is met, Fluid Mechanics Computation method CFD is re-used and is catalyzed The numerical simulation of cracking reaction, obtains the optimized product distributed data of final acquisition catalytic cracking reaction, realizes catalytic reaction Product hard measurement.Catalytic cracking reaction kinetic model is combined with Fluid Mechanics Computation method CFD, by optimizing kinetic simulation Shape parameter, realizes catalytic cracking reaction product hard measurement, improves the precision of catalytic cracking reaction product hard measurement result.This Outside, the embodiment of the present invention is when catalytic cracking operating mode changes, it is only necessary to corresponding adjustment catalytic cracking reaction kinetic simulation The model parameter of type, you can the corresponding calculating Product distribution data of the different operating modes of acquisition, realizes the Catalytic Cracking Unit of Measure of different operating modes Product hard measurement is answered, the adaptability of catalytic cracking reaction flexible measurement method is improved.
Fig. 7 is the structural representation of catalytic cracking reaction product hard measurement system in the embodiment of the present invention, as shown in fig. 7, Catalytic cracking reaction product hard measurement system provided in an embodiment of the present invention includes data acquisition unit 71, calculates product distribution list Member 72, Model Parameter Optimization unit 73 and hard measurement unit 74, wherein:
Data acquisition unit 71 is used to gather the corresponding floor data of the different operating modes of catalytic cracking reaction, the floor data Actual measurement Product distribution data and duty parameter including catalytic cracking reaction;Calculating product distribution unit 72 is used for according to the work Condition parameter, carries out numerical simulation to the catalytic cracking reaction kinetic model that pre-selection is set up using Fluid Mechanics Computation method, obtains The corresponding calculating Product distribution data of catalytic cracking reaction difference operating mode;If Model Parameter Optimization unit 73 is used to judge to know institute The deviation summation for calculating Product distribution data and the actual measurement Product distribution data is stated more than predetermined threshold value, then the catalysis is split The model parameter for changing reaction Kinetics Model is optimized, until the deviation summation is less than or equal to the predetermined threshold value, is obtained Obtain Optimized model parameter;Hard measurement unit 74 is used for according to the Optimized model parameter and the catalytic cracking reaction kinetic simulation Type, using the Fluid Mechanics Computation method, obtains the optimized product distributed data of catalytic cracking reaction, by the optimized product point Cloth data as the catalytic cracking reaction product hard measurement result.
Specifically, data acquisition unit 71 gathers the corresponding floor data of the different operating modes of catalytic cracking reaction as needed, Floor data includes the actual measurement Product distribution data and duty parameter of catalytic cracking reaction.Wherein duty parameter includes:Temperature number According to, pressure data and data on flows, actual measurement Product distribution data includes:Dry gas yield, yield of liquefied gas, yield of gasoline, diesel oil Yield and coking yield.Duty parameter can obtain real-time working condition parameter or from factory's history number by industrial site collection According to acquisition history duty parameter in the data preserved in storehouse.Actual measurement product point in real time can be obtained by industrial site collection History actual measurement Product distribution data is obtained in cloth data or the data preserved from plant historian.Get catalytic reaction Floor data after, calculate product distribution unit 72 and utilize Fluid Mechanics Computation method CFD according to these floor datas (Computational fluid dynamics) carries out numerical simulation to the catalytic cracking reaction kinetic model that pre-selection is set up, Obtain the corresponding calculating Product distribution data of the different operating modes of catalytic cracking reaction.Obtained according to the corresponding floor data of different operating modes Take corresponding calculating Product distribution data.Model Parameter Optimization unit 73 is by the calculating Product distribution data got and adopts in advance It is poor that the actual measurement Product distribution data collected is made, and judges that calculating Product distribution data and the deviation summation of actual measurement Product distribution data is It is no to be more than predetermined threshold value, optimized if the model parameter to catalytic cracking reaction kinetic model more than if, until described inclined Poor summation is less than or equal to the predetermined threshold value, obtains Optimized model parameter;If calculating Product distribution data and actual measurement product point The deviation summation of cloth data is less than or equal to predetermined threshold value, then by the model parameter of catalytic cracking reaction kinetic model now It is used as Optimized model parameter.Hard measurement unit 74 utilizes meter according to Optimized model parameter and catalytic cracking reaction kinetic model Fluid operator Mechanical Method CFD, obtains the corresponding optimized product distributed data of the different operating modes of catalytic cracking reaction, optimized product is distributed Data realize the product hard measurement under the different operating modes of catalytic cracking reaction as the result of catalytic cracking reaction product hard measurement.
It should be noted that when the operating mode of catalytic cracking reaction changes, floor data is adjusted accordingly, and urge Change the model parameter of cracking reaction kinetic model, by the new corresponding floor data of work condition environment and catalytic cracking reaction power The model parameter input for learning model calculates product distribution unit 72, Model Parameter Optimization unit 73 and hard measurement unit 74, again The hard measurement of catalytic cracking reaction product is carried out, the optimized product point of corresponding catalytic cracking reaction under new working condition is obtained Cloth data, realize the catalytic cracking reaction product hard measurement of different operating modes.The embodiment of the present invention is by coupling Fluid Mechanics Computation The model of CFD approach and the three-dimensional flow-reaction coupled model of genetic algorithm to ten four-lump kinetics models of catalytic cracking reaction Parameter is corrected, and hard measurement result is had good precision and adaptability.
Wherein, the catalytic cracking reaction kinetic model that the embodiment of the present invention is used for ten four-lump kinetics models, its It is a kind of reaction Kinetics Model simplified with lumping strategy, the method for building up of the model is consistent with above-described embodiment, herein Repeat no more.
System provided in an embodiment of the present invention is used to perform the above method, and its embodiment is with above-described embodiment one Cause, here is omitted.
The method and system of the catalytic cracking reaction product hard measurement that example is provided in real time of the invention, utilizes Fluid Mechanics Computation Method CFD carries out the numerical simulation of catalytic cracking reaction, and the model parameter of catalytic cracking reaction kinetic model is optimized, The model parameter of the catalytic cracking reaction kinetic model for the condition that meets is obtained, Fluid Mechanics Computation method CFD is re-used and carries out The numerical simulation of catalytic cracking reaction, obtains the optimized product distributed data of final acquisition catalytic cracking reaction, realizes catalysis Reaction product hard measurement.Catalytic cracking reaction kinetic model is combined with Fluid Mechanics Computation method CFD, by optimizing power Model parameter is learned, catalytic cracking reaction product hard measurement is realized, improves the precision of catalytic cracking reaction product hard measurement result. In addition, the embodiment of the present invention is when catalytic cracking operating mode changes, it is only necessary to corresponding adjustment catalytic cracking reaction dynamics The model parameter of model, you can the corresponding calculating Product distribution data of the different operating modes of acquisition, realizes the catalytic cracking of different operating modes Reaction product hard measurement, improves the adaptability of catalytic cracking reaction flexible measurement method.
Fig. 8 be the embodiment of the present invention in be used for catalytic cracking reaction product hard measurement electronic equipment structural representation, As shown in figure 8, described device can include:Processor (processor) 81, memory (memory) 82 and communication bus 83, Wherein, processor 81, memory 82 completes mutual communication by communication bus 83.Processor 81 can call memory 82 In logical order, to perform following method:Gather the corresponding floor data of the different operating modes of catalytic cracking reaction, the operating mode number According to actual measurement Product distribution data and duty parameter including catalytic cracking reaction;According to the duty parameter, using calculating fluid Mechanical Method carries out numerical simulation to the catalytic cracking reaction kinetic model that pre-selection is set up, and obtains the different operating modes of catalytic cracking reaction Corresponding calculating Product distribution data;If judging to know the calculating Product distribution data and the actual measurement Product distribution data Deviation summation is more than predetermined threshold value, then the model parameter to the catalytic cracking reaction kinetic model is optimized, until institute Deviation summation is stated less than or equal to the predetermined threshold value, Optimized model parameter is obtained;According to the Optimized model parameter and described Catalytic cracking reaction kinetic model, using the Fluid Mechanics Computation method, obtains the different operating modes of catalytic cracking reaction corresponding Optimized product distributed data, using the optimized product distributed data as the catalytic cracking reaction product hard measurement result.
In addition, the logical order in above-mentioned memory 82 can be realized by the form of SFU software functional unit and is used as solely Vertical production marketing in use, can be stored in a computer read/write memory medium.Understood based on such, this hair The part or the part of the technical scheme that bright technical scheme substantially contributes to prior art in other words can be with soft The form of part product is embodied, and the computer software product is stored in a storage medium, including some instructions are to make Obtain a computer equipment (can be personal computer, server, or network equipment etc.) and perform each embodiment of the invention The all or part of step of methods described.And foregoing storage medium includes:USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. it is various Can be with the medium of store program codes.
The embodiment of the present invention provides a kind of non-transient computer readable storage medium storing program for executing, the non-transient computer readable storage Medium storing computer is instructed, and the computer instruction makes the computer perform the side that above-mentioned each method embodiment is provided Method, for example including:The corresponding floor data of the different operating modes of catalytic cracking reaction is gathered, the floor data includes Catalytic Cracking Unit of Measure The actual measurement Product distribution data and duty parameter answered;According to the duty parameter, pre-selection is set up using Fluid Mechanics Computation method Catalytic cracking reaction kinetic model carry out numerical simulation, obtain that the different operating mode of catalytic cracking reaction is corresponding to calculate product point Cloth data;If judging the deviation summation for knowing the calculating Product distribution data and the actual measurement Product distribution data more than default Threshold value, then the model parameter to the catalytic cracking reaction kinetic model optimize, until the deviation summation be less than or Equal to the predetermined threshold value, Optimized model parameter is obtained;According to the Optimized model parameter and the catalytic cracking reaction power Model is learned, using the Fluid Mechanics Computation method, the corresponding optimized product distributed data of the different operating modes of catalytic cracking reaction is obtained, Using the optimized product distributed data as the catalytic cracking reaction product hard measurement result.
Above example is merely to illustrate technical scheme, rather than its limitations;Although with reference to the foregoing embodiments The present invention is described in detail, it will be understood by those within the art that:It still can be to foregoing each implementation Technical scheme described in example is modified, or carries out equivalent substitution to which part technical characteristic;And these are changed or replaced Change, the essence of appropriate technical solution is departed from the spirit and scope of various embodiments of the present invention technical scheme.

Claims (10)

1. a kind of method of catalytic cracking reaction product hard measurement, it is characterised in that including:
The corresponding floor data of the different operating modes of catalytic cracking reaction is gathered, the floor data includes the actual measurement of catalytic cracking reaction Product distribution data and duty parameter;
According to the duty parameter, the catalytic cracking reaction kinetic model that pre-selection is set up is carried out using Fluid Mechanics Computation method Numerical simulation, obtains the corresponding calculating Product distribution data of the different operating modes of catalytic cracking reaction;
If judging, the deviation summation for knowing the calculating Product distribution data and the actual measurement Product distribution data is more than default threshold Value, then the model parameter to the catalytic cracking reaction kinetic model is optimized, until the deviation summation is less than or waited In the predetermined threshold value, Optimized model parameter is obtained;
According to the Optimized model parameter and the catalytic cracking reaction kinetic model, using the Fluid Mechanics Computation method, The corresponding optimized product distributed data of the different operating mode of catalytic cracking reaction is obtained, using the optimized product distributed data as described The result of catalytic cracking reaction product hard measurement.
2. according to the method described in claim 1, it is characterised in that the utilization Fluid Mechanics Computation method is urged pre-selection foundation Change cracking reaction kinetic model and carry out numerical simulation, obtain the corresponding calculating product distribution number of the different operating modes of catalytic cracking reaction According to, including:
If the numerical simulation is first simulation, model parameter initial value is set, and according to the corresponding work of different operating modes Condition parameter, sets parameter of a flow model;
According to the model parameter initial value and the parameter of a flow model, pre-selection foundation is urged using Fluid Mechanics Computation method Change cracking reaction kinetic model and carry out numerical simulation, obtain the corresponding calculating product point of difference operating modes described in catalytic cracking reaction Cloth data.
3. according to the method described in claim 1, it is characterised in that if the judgement know the calculating Product distribution data and The deviation summation of the actual measurement Product distribution data is more than predetermined threshold value, then to the mould of the catalytic cracking reaction kinetic model Shape parameter is optimized, including:
If judgement knows that the deviation summation is more than predetermined threshold value, the model parameter is optimized using genetic algorithm.
4. according to the method described in claim 1, it is characterised in that described to be split according to the Optimized model parameter and the catalysis Change reaction Kinetics Model, using the Fluid Mechanics Computation method, obtain the optimized product distributed data of catalytic cracking reaction, bag Include:
The Optimized model parameter is substituted into the catalytic cracking reaction kinetic model, using the Fluid Mechanics Computation method, Numerical simulation is carried out to the catalytic cracking reaction kinetic model, the optimized product distributed data is obtained.
5. according to the method described in claim 1, it is characterised in that the catalytic cracking reaction kinetic model is 14 lumps Kinetic model.
6. according to the method described in claim 1, it is characterised in that the duty parameter includes:Temperature data, pressure data and Data on flows.
7. according to the method described in claim 1, it is characterised in that the actual measurement Product distribution data includes:Dry gas yield, liquid Change gas yield, yield of gasoline, diesel yield and coking yield.
8. a kind of system of catalytic cracking reaction product hard measurement, it is characterised in that including:
Data acquisition unit, for gathering the corresponding floor data of the different operating modes of catalytic cracking reaction, the floor data includes The actual measurement Product distribution data and duty parameter of catalytic cracking reaction;
Product distribution unit is calculated, for according to the duty parameter, using Fluid Mechanics Computation method to preselecting the catalysis set up Cracking reaction kinetic model carries out numerical simulation, obtains the corresponding calculating product distribution number of the different operating modes of catalytic cracking reaction According to;
Model Parameter Optimization unit, if for judging to know the calculating Product distribution data and the actual measurement Product distribution data Deviation summation be more than predetermined threshold value, then the model parameter to the catalytic cracking reaction kinetic model optimize, until The deviation summation is less than or equal to the predetermined threshold value, obtains Optimized model parameter;
Hard measurement unit, for according to the Optimized model parameter and the catalytic cracking reaction kinetic model, using described Fluid Mechanics Computation method, obtains the corresponding optimized product distributed data of the different operating modes of catalytic cracking reaction, by the optimized product Distributed data as the catalytic cracking reaction product hard measurement result.
9. a kind of electronic equipment for catalytic cracking reaction product hard measurement, it is characterised in that including:
At least one processor;And
At least one memory being connected with the processor communication, wherein:
The memory storage has can be by the programmed instruction of the computing device, and the processor calls described program to instruct energy Enough perform the method as described in any one of claim 1 to 7.
10. a kind of non-transient computer readable storage medium storing program for executing, it is characterised in that the non-transient computer readable storage medium storing program for executing is deposited Computer instruction is stored up, the computer instruction makes the computer perform the method as described in any one of claim 1 to 7.
CN201710305383.7A 2017-05-03 2017-05-03 A kind of method and system of catalytic cracking reaction product hard measurement Pending CN107291975A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710305383.7A CN107291975A (en) 2017-05-03 2017-05-03 A kind of method and system of catalytic cracking reaction product hard measurement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710305383.7A CN107291975A (en) 2017-05-03 2017-05-03 A kind of method and system of catalytic cracking reaction product hard measurement

Publications (1)

Publication Number Publication Date
CN107291975A true CN107291975A (en) 2017-10-24

Family

ID=60094494

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710305383.7A Pending CN107291975A (en) 2017-05-03 2017-05-03 A kind of method and system of catalytic cracking reaction product hard measurement

Country Status (1)

Country Link
CN (1) CN107291975A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114446397A (en) * 2020-10-19 2022-05-06 中国石油化工股份有限公司 Memory, refining process flow segmentation simulation and optimization method, device and equipment
CN114495686A (en) * 2022-01-26 2022-05-13 华东理工大学 Real-time simulation method and system for industrial catalytic cracking device
CN114638030A (en) * 2021-11-30 2022-06-17 国家电投集团科学技术研究院有限公司 Thermal hydraulic model correction method and device based on expert knowledge base
CN115841050A (en) * 2023-02-20 2023-03-24 新疆独山子石油化工有限公司 Method for calculating premixing temperature of catalytic cracking regenerated catalyst and molecular composition raw oil
CN116821695A (en) * 2023-08-30 2023-09-29 中国石油大学(华东) Semi-supervised neural network soft measurement modeling method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1563983A (en) * 2004-04-15 2005-01-12 石油大学(北京) Soft measuring method for yield rate of catalytic cracking reaction
CN101520453A (en) * 2009-02-10 2009-09-02 浙江大学 Soft measuring instrument of dry point of catalytic crude gasoline
CN101980230A (en) * 2010-11-03 2011-02-23 天津大学 Catalytic cracking reaction system process simulation optimization model and solution method thereof
CN102004846A (en) * 2010-11-03 2011-04-06 天津大学 Catalytic cracking reaction kinetics model of distillate
CN102663235A (en) * 2012-03-20 2012-09-12 浙江大学 Modeling method for catalytic cracking main fractionator with varying-population-size DNA genetic algorithm
CN105631107A (en) * 2015-12-23 2016-06-01 浙江中控软件技术有限公司 Prediction method and apparatus for fluidized catalytic cracking reaction

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1563983A (en) * 2004-04-15 2005-01-12 石油大学(北京) Soft measuring method for yield rate of catalytic cracking reaction
CN101520453A (en) * 2009-02-10 2009-09-02 浙江大学 Soft measuring instrument of dry point of catalytic crude gasoline
CN101980230A (en) * 2010-11-03 2011-02-23 天津大学 Catalytic cracking reaction system process simulation optimization model and solution method thereof
CN102004846A (en) * 2010-11-03 2011-04-06 天津大学 Catalytic cracking reaction kinetics model of distillate
CN102663235A (en) * 2012-03-20 2012-09-12 浙江大学 Modeling method for catalytic cracking main fractionator with varying-population-size DNA genetic algorithm
CN105631107A (en) * 2015-12-23 2016-06-01 浙江中控软件技术有限公司 Prediction method and apparatus for fluidized catalytic cracking reaction

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
刘方: "聚酯反应动力学研究及反应器内流场的模拟与优化", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *
洪建旭: "工业催化裂化装置反再系统建模研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *
王敏 等: "FCC提升管反应器中终止剂注入对裂化反应的影响", 《化工学报》 *
袁慎芳: "《结构健康监控》", 30 April 2007 *
高成 等: "软测量技术", 《传感器与检测技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114446397A (en) * 2020-10-19 2022-05-06 中国石油化工股份有限公司 Memory, refining process flow segmentation simulation and optimization method, device and equipment
CN114638030A (en) * 2021-11-30 2022-06-17 国家电投集团科学技术研究院有限公司 Thermal hydraulic model correction method and device based on expert knowledge base
CN114495686A (en) * 2022-01-26 2022-05-13 华东理工大学 Real-time simulation method and system for industrial catalytic cracking device
CN115841050A (en) * 2023-02-20 2023-03-24 新疆独山子石油化工有限公司 Method for calculating premixing temperature of catalytic cracking regenerated catalyst and molecular composition raw oil
CN116821695A (en) * 2023-08-30 2023-09-29 中国石油大学(华东) Semi-supervised neural network soft measurement modeling method
CN116821695B (en) * 2023-08-30 2023-11-03 中国石油大学(华东) Semi-supervised neural network soft measurement modeling method

Similar Documents

Publication Publication Date Title
CN107291975A (en) A kind of method and system of catalytic cracking reaction product hard measurement
CN104789256B (en) A kind of yield real-time predicting method of catalytic cracking unit
CN103914595B (en) Medium temperature coal tar hydrogenation of total effluent cracking lumped reaction kinetics modeling method
Amblard et al. CFD modeling of the coke combustion in an industrial FCC regenerator
CN104804761B (en) A kind of yield real-time predicting method of hydrocracking unit
Han et al. Modeling and optimization of a fluidized catalytic cracking process under full and partial combustion modes
CN104392098A (en) Method for predicting yield of catalytically cracked gasoline
CN104463343A (en) Method for predicting catalytic cracking light oil yield
Sankararao et al. Multi-objective optimization of an industrial fluidized-bed catalytic cracking unit (FCCU) using two jumping gene adaptations of simulated annealing
Bisotti et al. Impact of kinetic models on methanol synthesis reactor predictions: in silico assessment and comparison with industrial data
CN102213949A (en) Method for optimizing value of ethylene device
CN111863145B (en) Low-temperature coal tar full-fraction hydrocracking lumped dynamics model modeling method
CN107301264A (en) The evaluation method and system of a kind of catalytic cracking reaction kinetic parameter
Arandes et al. Modelling FCC units under steady and unsteady state conditions
CN101980230B (en) Catalytic cracking reaction system process simulation optimization model and solution method thereof
Chin et al. Modeling and simulation study of an industrial radial moving bed reactor for propane dehydrogenation process
Gubaydullin et al. Numerical simulation of oxidative regeneration of a spherical catalyst grain
Chen et al. A mass-temperature decoupled discretization strategy for large-scale molecular-level kinetic model
Taşkin et al. FUZZYFCC: Fuzzy logic control of a fluid catalytic cracking unit (FCCU) to improve dynamic performance
Ahmed et al. Modelling and simulation of fluid catalytic cracking unit
CN112770829A (en) Modeling of operational and/or dimensional parameters of a gas treatment device
Li et al. A bilevel data-driven framework for robust optimization under uncertainty–applied to fluid catalytic cracking unit
CN112365071A (en) Real-time optimization method of S Zorb device based on improved differential evolution algorithm
CN103849423B (en) Modified catalytic cracking process device and catalytic cracking process
Jarullah et al. Energy consumption and heat recovery of an industrial fluidized catalytic cracking process based on cost savings

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20171024

RJ01 Rejection of invention patent application after publication