The whole reaction time t of delayed coking process was divided
into n intervals. Each interval was Δti (i =1-n ). In the reaction time
Δt1 , the molecular composition matrix of feed oil
and reaction temperature T were input, and the reaction network
was established based on the reaction rules of the SOL method. Combining
with the reaction rate constants, the reaction kinetics differential
equations were constructed, and then solved using the improvedRunge-Kutta method to get the product molecular composition
matrix of Δt1 , which was the feed molecular
composition matrix of the reaction time Δt2 . At
the same time, the total reaction heat in the reaction time
Δt1 was calculated and converted into the
temperature difference ΔT1 to obtain the reaction
temperature
Δt2 =T +ΔT1 for the
next reaction interval. Then the coking process of the reaction time
Δt2 was calculated. The calculation was repeated
until the reaction time reached t . The calculation diagram of the
delayed coking model was shown in Figure 2.