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 =TT1 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.