Industrial Investigation of the Combined Action of Vacuum Residue Hydrocracking and Vacuum Gas Oil Catalytic Cracking While Processing Different Feeds and Operating under Distinct Conditions

: Ebullated bed vacuum residue hydrocracking and ﬂuid catalytic cracking (FCC) are among the most proﬁtable processes in modern reﬁning. Their optimal performance is vital for petroleum reﬁning proﬁtability. That is why a better understanding of their combined action and the interrelations between these two heavy oil conversion processes in a real-world reﬁnery could provide valuable information for further performance optimization. Nine distinct petroleum crudes belonging to the extra light, light, and medium petroleum crude types were processed in the LUKOIL Neftohim Burgas reﬁnery to study the combined performance of two processes: FCC of vacuum gas oil and ebullated bed vacuum residue H-Oil hydrocracking. The operating conditions along with the characterization data of the feeds and products of both processes were evaluated through the employment of intercriteria analysis to deﬁne the variables with statistically signiﬁcant relationships. Maple 2023 Academic Edition mathematics software was used to develop models to predict the vacuum residue conversion level under different operating conditions. The plug ﬂow reactor model with an activation energy of 215 kJ/mol and a reaction order of 1.59 was found to provide the highest accuracy of vacuum residue conversion, with an average absolute deviation of 2.2%. H-Oil yields were found to correlate with the vacuum residue conversion level and the content of FCC slurry oil (SLO), the recycling of partially blended fuel oil, a material boiling point below 360 ◦ C, and the vacuum gas oil (VGO) in the H-Oil feed. FCC conversion was found to depend on the H-Oil VGO content in the FCC feed and the content of FCC SLO in the H-Oil feed.


Introduction
World climate change necessitates the development of measures to lessen carbon emissions-from traditional fossil fuels to renewable and sustainable clean energy.The truth is, however, that an immediate turn to renewables or clean energy is very difficult because of underdeveloped technologies, inadequate resources, and an unsteady provision of energy [1][2][3].Over the next few decades, fossil fuels are expected to remain the main source of energy for the transport sector.Therefore, instead of neglecting fossil fuels altogether and looking for impractical technologies to reduce carbon emissions, improvement in the performance of energy-intensive refining plants can deliver the desired reduction in greenhouse gas emissions [1][2][3].This study was devoted to the thermo-catalytic industrial processes of heavy petroleum fractions-fluid catalytic cracking and hydrocracking-and searching for possibilities to improve their performance.
Heavy oil petroleum fractions are those that cannot be removed by the atmospheric distillation.Their boiling point is usually above 350 • C, and they have a specific gravity higher than 0.935 at 15.6 • C [4].Due to their significantly lower demand (about 15% versus 49% for middle distillates and 36% for light distillates), their price is considerably lower than that of the light petroleum fractions [5].This is why the transformation of heavy oil fractions into transportation fuels and feeds for petrochemistry, which takes place in conversion units, makes them a key factor in modern refining [6][7][8][9].
Heavy oils are separated in vacuum gas oils and vacuum residues during vacuum distillation [10].Vacuum gas oils are typically processed in fluid catalytic cracking (FCC) or hydrocracking, while vacuum residue can be converted by several processes classified in two groups: carbon rejection and hydrogen addition [11][12][13][14][15][16][17][18].The hydrocracking of vacuum gas oils is commonly used in the industry in two fixed bed reactors, the first designed for hydrotreating VGO and the second for hydrocracking the hydrotreated VGO [19,20].This processing scheme, unfortunately, is vulnerable to the treatment of secondary VGOs coming from residue coking and hydroprocessing due to their high content of nitrogen, carbon residue, and metals, which are hydrocracking catalyst poisons.FCC seems to be the better option for the conversion of secondary VGOs due to the higher tolerance of FCC catalysts.The highest conversion level of vacuum residues was reported to be achieved by two commercially proven hydroprocessing technologies-slurry hydrocracking [21,22] (95% vacuum residue conversion) and H-Oil ebullated bed hydrocracking (93% vacuum residue conversion) [23].The use of FCC and H-Oil in the LUKOIL Neftohim Burgas (LNB) refinery enabled a light petroleum product yield enhancement from 64 to 81%.The synergy between FCC and ebullated bed hydrocracking was thoroughly investigated in a study by Stratiev at al. [23].The lower demand for naphtha produced by H-Oil hydrocrackers in 2023 necessitated a reduction in process severity.While the operation of vacuum residue hydrocracking at a high severity has been extensively investigated and discussed [21][22][23][24][25][26][27], it has not been studied or discussed at a low severity to the best of our knowledge.H-Oil at a lower severity along with the refinement of various types of crude oils and their relationship with FCC unit performance have not been investigated or discussed yet, especially on an industrial scale.In the published research, knowledge about the industrial performance of both ebullated bed vacuum residue operating at low severity and FCC, which processes the hydrocracked VGO obtained at a low severity, is completely missing.The severity of vacuum residue hydrocracking depends on the reaction time or liquid hourly space velocity (LHSV) and reaction temperature.The LHSV design of the LUKOIL Neftohim Burgas (LNB) H-Oil hydrocracker is 0.25 h −1 , and the design reaction temperature is 423 • C, providing a vacuum residue conversion of 70 wt.%when Urals vacuum residue is processed [27].During the COVID 19 pandemic, the LNB H-Oil hydrocracker operated at an LHSV as low as 0.10 h −1 , and the reaction temperature was as high as 436 • C, achieving a vacuum residue conversion of 93 wt.% [23].Starting from July 2023, the H-Oil hydrocracker began to operate at a reaction temperature of 405 • C and an LHSV of about 0.13 h −1 .The influence of the low severity of the ebullated bed hydrocracking H-Oil unit operation on FCC unit performance has also not been investigated yet.This is the reason we conducted this study, in which H-Oil conversion was varied between 51.8 and 91.7 wt.% and FCC conversion varied between 71.5 and 78.6%.The impact of the feed quality and operating conditions in both processes, H-Oil and FCC, on vacuum residue, vacuum gas oil conversion levels, hydrodeasphaltization, hydrodemetalation, and the quality of the products was investigated using intercriteria analysis.
The purpose of this research was to quantitatively determine the relationship between the conversion, product yields, extent of metals, sulfur removal, and product quality of the H-Oil and FCC processes and the operation conditions, crude slate processing, and processing of recycled and FCC decant oil in the H-Oil hydrocracker.

Materials and Methods
The characteristics of the nine petroleum crudes refined in the LNB are presented in Table 1.Table 2 summarizes characteristics of the vacuum residues derived from the nine crude oils that were processed in the H-Oil unit.
Table 3 shows the characteristics of the vacuum gas oils derived from the nine crude oils.
A technological scheme of the H-Oil vacuum residue hydrocracker and its detailed explanation are shown in Figure 1.A detailed explanation of the H-Oil hydrocracker and its application in the LNB refinery the reader can be found in our recent research [23].The catalyst used in the H-Oil unit (Figure 1) was a low-sediment Ni-Mo catalyst.Notes: * This is the kinematic viscosity of mixtures of 70%VR/30% FCC HCO at 80 • C. ** The molecular weight was calculated using the correlation established in our recent research [28].
The process scheme of the studied vacuum gas oil fluid catalytic cracker is displayed in Figure 2. The catalyst used in the FCC unit (Figure 2) during the study was an octane barreltype catalyst, and its properties were summarized in our recent research [26] (Catalyst H from Table 2 in [29]).oils.A technological scheme of the H-Oil vacuum residue hydrocracker and its detailed explanation are shown in Figure 1.A detailed explanation of the H-Oil hydrocracker and its application in the LNB refinery the reader can be found in our recent research [23].The catalyst used in the H-Oil unit (Figure 1) was a low-sediment Ni-Mo catalyst.The process scheme of the studied vacuum gas oil fluid catalytic cracker is displayed in Figure 2. The catalyst used in the FCC unit (Figure 2) during the study was an octane barrel-type catalyst, and its properties were summarized in our recent research [26] (Catalyst H from Table 2 in [29]).A simplified layout of the heavy oil technological chain in the LUKOIL Neftohim Burgas refinery that includes the H-Oil and FCC processes is presented in Figure 3.
Figure 3 shows that the H-Oil and FCC processes in the present study are parallelsequential with an intricate relationship between them, in which the H-Oil VGO product feeds the FCC unit along with the hydrotreated VGO (HTVGO) and FCC HCO, and the SLO products are feed components for the H-Oil unit together with the straight-run vacuum residue.
The SARA composition of the primary and secondary (hydrocracked) vacuum residues, the H-Oil feed, and the vacuum gas oils were analyzed by the liquid chromatography inhouse method described in [30].
A net 540 • C+ vacuum residue conversion, which involves the transformation of material boiling above 540 • C from the feed into material boiling below 540 • C (percentage of feed), asphaltene conversion (HDAs), and hydrodemetallization (HDM) extent were computed following Equations ( 1)-( 3), as reported in our recent investigation [31].A simplified layout of the heavy oil technological chain in the LUKOIL Neftohim Burgas refinery that includes the H-Oil and FCC processes is presented in Figure 3.A simplified layout of the heavy oil technological chain in the LUKOIL Neftohim Burgas refinery that includes the H-Oil and FCC processes is presented in Figure 3.The sediment content in the H-Oil atmospheric tower bottom product (ATB) was determined following the procedure IP 375 (total sediment existent (TSE)).By contrast, the sediment content after thermal aging in the partially blended fuel oil (PBFO), which was prepared from about 60% H-Oil VTB (vacuum tower bottom product), about 30% fluid catalytic cracking heavy-cycle oil (FCC HCO), and about 10% FCC light cycle oil (LCO), and slurry oil (SLO) was determined following the procedure IP 390 (total sediment potential (TSP)).The composition of PBFO in the study is shown in Table S1.
The conversion of vacuum gas oil (VGO) in the FCC unit was computed as reported in our recent research [29].
Figure 4 shows how the contents of individual petroleum crudes in processed petroleum mixtures varied during the study.
A net 540 °C+ vacuum residue conversion, which involves the transformation of material boiling above 540 °C from the feed into material boiling below 540 °C (percentage of feed), asphaltene conversion (HDAs), and hydrodemetallization (HDM) extent were computed following Equations ( 1)-( 3), as reported in our recent investigation [31].
The sediment content in the H-Oil atmospheric tower bottom product (ATB) was determined following the procedure IP 375 (total sediment existent (TSE)).By contrast, the sediment content after thermal aging in the partially blended fuel oil (PBFO), which was prepared from about 60% H-Oil VTB (vacuum tower bottom product), about 30% fluid catalytic cracking heavy-cycle oil (FCC HCO), and about 10% FCC light cycle oil (LCO), and slurry oil (SLO) was determined following the procedure IP 390 (total sediment potential (TSP)).The composition of PBFO in the study is shown in Table S1.
The conversion of vacuum gas oil (VGO) in the FCC unit was computed as reported in our recent research [29].
Figure 4 shows how the contents of individual petroleum crudes in processed petroleum mixtures varied during the study.During the study (10 April 2022-24 July 2023), the LNB refinery processed the mixtures of the petroleum crudes Arab Light, Arab Medium, Basrah Heavy, Johan Sverdrup, Kirkuk, Light Siberian, Urals, and Western Desert.The contents of individual petroleum crudes in processed mixtures are shown in Figure 4.One can see from the data in Figure 4 that the main crude oils processed were Urals and Light Siberian crude oils.Their mutual content in the crude oil blend varied between 74.5 and 100%.During the study (10 April 2022-24 July 2023), the LNB refinery processed the mixtures of the petroleum crudes Arab Light, Arab Medium, Basrah Heavy, Johan Sverdrup, Kirkuk, Light Siberian, Urals, and Western Desert.The contents of individual petroleum crudes in processed mixtures are shown in Figure 4.One can see from the data in Figure 4 that the main crude oils processed were Urals and Light Siberian crude oils.Their mutual content in the crude oil blend varied between 74.5 and 100%.
An intercriteria analysis (ICrA) assessment of the data for both the H-Oil and FCC processes was carried out to search for the statistically significant relationships between the studied variables.A µ-value = 0.70 ÷ 1.00; υ = 0 ÷ 0.30 indicates a statistically significant positive relationship, in which a strong positive association exhibits values of µ = 0.90 ÷ 1.00; υ = 0 ÷ 0.1 and a weak positive association demonstrates values of µ = 0.70 ÷ 0.80; υ = 0.20 ÷ 0.30.Accordingly, the values of negative association with µ = 0 ÷ 0.30; υ = 0.70 ÷ 1.00 allude to a statistically significant negative relationship, in which a strong negative association displays values of µ = 0 ÷ 0.1; υ = 0.90 ÷ 1.00 and a weak negative association shows values of µ = 0.20 ÷ 0.30; υ = 0.70 ÷ 0.80.All other cases are considered dissonance.A detailed explanation of ICrA application in petroleum chemistry and processing can be found in [32].
Two models were developed to relate vacuum residue conversion to reaction temperature and reaction time.The first model was a regression similar to that developed in our recent research [31].The second used the kinetic expression employed for a continuous stirred tank reactor (CSTR) because ebullated bed reactors are assumed to behave as such a reactor [33].The mass balance of the CSTR is given as Equation (1) [34].
where: τ = reaction time = 1 LHSV , hour; Conv = 540 • C+ net conversion (calculated by Equation ( 1)); k 0 = collision factor; h −1 ×frac.−n ; E A = activation energy, kJ/mol; R = universal gas constant; (R = 8.31446261815324 joules/mol.K); T = reaction temperature, K; n = reaction order.After the rearrangement of Equation ( 1), the expression for conversion takes the form of Equation ( 2): The objective function (the sum of the square of the difference between the experimental and calculated conversions) was minimized, and the best set of kinetic parameters (k 0 , E A , and n) was obtained, as shown in Equation (3): The bounded nonlinear optimization problem in Equation ( 3) is iteratively minimized using a differential evolution algorithm (diffevol) in Maple's Global Optimization Package.
Schweitzer and Kressmann, who studied ebullated bed reactor modeling, assumed that the ebullated bed reactor is modelled by the plug flow model [35].A very simplified mass balance equation of the plug flow reactor is shown in Equation (4): where: X = conversion of the vacuum residue (540 • C+), frac.; t = reaction time, h; t = 1/LHSV; n = reaction order.

Investigation of H-Oil Vacuum Residue Hydrocracking
The main problems in the ebullated bed vacuum residue hydrocracking operation are sediment formation and the related fouling [23,27,[36][37][38][39][40][41][42][43].The increase in reaction temperature typically leads to the exponential enhancement of the sediment formation rate and fouling [42,43].Figure 5 presents a graph of the H-Oil atmospheric tower bottom (ATB) sediment content variation (TSE) and that of the partially blended fuel oil (PBFO) after thermal aging (TSP) along with the weighted average bed temperature (WABT) of both H-Oil hydrocracking reactors.

Investigation of H-Oil Vacuum Residue Hydrocracking
The main problems in the ebullated bed vacuum residue hydrocracking operation are sediment formation and the related fouling [23,27,[36][37][38][39][40][41][42][43].The increase in reaction temperature typically leads to the exponential enhancement of the sediment formation rate and fouling [42,43].Figure 5 presents a graph of the H-Oil atmospheric tower bottom (ATB) sediment content variation (TSE) and that of the partially blended fuel oil (PBFO) after thermal aging (TSP) along with the weighted average bed temperature (WABT) of both H-Oil hydrocracking reactors.The data in Figure 5 indicate that irrespective of the WABT enhancement from 414.5 to 430 °C, the ATB TSE did not show any trend toward augmentation.Instead, the ATB TSE remained under the safe limit of about 0.3 wt.%.The reduction in WABT from 430 to 405 °C also showed no effect on the ATB TSE at the end of the study period.This supports the correctness of the decision to replace the cascade mode of fresh catalyst addition with the parallel mode and the optimization of the fresh catalyst addition rate, as discussed in our recent research [31].The data shown in Figure 6 confirm the very low fouling rate registered during this study.The data in Figure 5 indicate that irrespective of the WABT enhancement from 414.5 to 430 • C, the ATB TSE did not show any trend toward augmentation.Instead, the ATB TSE remained under the safe limit of about 0.3 wt.%.The reduction in WABT from 430 to 405 • C also showed no effect on the ATB TSE at the end of the study period.This supports the correctness of the decision to replace the cascade mode of fresh catalyst addition with the parallel mode and the optimization of the fresh catalyst addition rate, as discussed in our recent research [31].The data shown in Figure 6 confirm the very low fouling rate registered during this study.
The unusually high PBFO TSP observed on 18 April 2023, 24 April 2023, and 28 April 2023 can be explained by the shutdown of the FCC unit for repairs during the period of 14 April 2023-30 April 2023.The lack of highly aromatic FCC gas oils (LCO, HCO, and SLO) due to the FCC unit shutdown negatively affected the sediment content after the thermal aging of the PBFO (Figure 5).This was reported in our earlier studies [23,27].
Figure 7 shows the graphs of variation of TSE and the TSP of blends of H-Oil VTB obtained at very low severities (WABT of 405 • C); high aromatic FCC LCO, HCO, and SLO; and low aromatic FCCPT diesel with the alteration of diluent concentration.The properties of these four gas oils are shown in Table S2.The unusually high PBFO TSP observed on 18 April 2023, 24 April 2023, and 28 April 2023 can be explained by the shutdown of the FCC unit for repairs during the period of 14 April 2023-30 April 2023.The lack of highly aromatic FCC gas oils (LCO, HCO, and SLO) due to the FCC unit shutdown negatively affected the sediment content after the thermal aging of the PBFO (Figure 5).This was reported in our earlier studies [23,27].
Figure 7 shows the graphs of variation of TSE and the TSP of blends of H-Oil VTB obtained at very low severities (WABT of 405 °C); high aromatic FCC LCO, HCO, and SLO; and low aromatic FCCPT diesel with the alteration of diluent concentration.The properties of these four gas oils are shown in Table S2.The unusually high PBFO TSP observed on 18 April 2023, 24 April 2023, and 28 April 2023 can be explained by the shutdown of the FCC unit for repairs during the period of 14 April 2023-30 April 2023.The lack of highly aromatic FCC gas oils (LCO, HCO, and SLO) due to the FCC unit shutdown negatively affected the sediment content after the thermal aging of the PBFO (Figure 5).This was reported in our earlier studies [23,27].
Figure 7 shows the graphs of variation of TSE and the TSP of blends of H-Oil VTB obtained at very low severities (WABT of 405 °C); high aromatic FCC LCO, HCO, and SLO; and low aromatic FCCPT diesel with the alteration of diluent concentration.The properties of these four gas oils are shown in Table S2.It is obvious from the data in Figure 7 that the sediment content of the blend of H-Oil VTB-diluent is specific in relationship with the diluent, a finding that was also observed in our recent research [44].The FCC HCO turned out to be the most effective diluent related to sediment content reduction, while the low aromatic FCCPT diesel exhibited the highest TSP, confirming that the shutdown of the FCC unit results in PBFO TSP enlarge- It is obvious from the data in Figure 7 that the sediment content of the blend of H-Oil VTB-diluent is specific in relationship with the diluent, a finding that was also observed in our recent research [44].The FCC HCO turned out to be the most effective diluent related to sediment content reduction, while the low aromatic FCCPT diesel exhibited the highest TSP, confirming that the shutdown of the FCC unit results in PBFO TSP enlargement.It deserves mentioning here that at a higher reaction severity (WABT of 430 • C), the TSE of the blend of H-Oil VTB-FCCPT diesel was about 1.0 wt.% [44], while at a WABT of 405 • C, as obvious from the data in Figure 7d, it was about 0.2 wt.%.This suggests that at higher a reaction severity, the H-Oil VTB asphaltenes become less soluble.
The operating variables affecting vacuum residue conversion, hydrodeasphaltization (HDAs) extent, and hydrodemetallization (HDM), such as the reaction temperature and reaction time ( 1LHSV ), are summarized in Table S3.Using the data from Table S3 and considering that the conversion depends on the reaction temperature (WABT, • C) and reaction time ( 1LHSV) ), two models, as mentioned in Section 2 (Materials and Methods), were developed to relate the vacuum residue conversion to the reaction temperature and reaction time.
The regression equation that predicts the conversion takes the form of Equation ( 9): where: LHSV = liquid hourly space velocity = Feed f low rate EBR volume ; h −1 ; feed flow rate in m 3 /h; EBR = ebullated bed reactor volume in m 3 ; WABT = weighted average bed temperature, • C. For the dataset in Table S3, the following kinetic parameters for the CSTR model (Equation (2)) were found: k 0 = 1.4369 × 10 11 ; E A = 152.67kJ/mol; n = 1.82.
For the dataset in Table S3, the following kinetic parameters for the plug flow model (Equation (8)) were found: k 0 = 4.42875 × 10 15 ; E A = 215 kJ/mol; n = 1.59, following the iteratively minimization of Equation ( 8) using the differential evolution algorithm (diffevol) in Maple's Global Optimization Package.
In Figure 8, the training set is related to development of Equations ( 2), ( 8) and ( 9) with the data in Table S3.The test set was used to verify Equations ( 2), ( 8) and ( 9) with data not included in the process of equation development.
The statistical parameters used to evaluate the precision of conversion prediction are presented in Table 4.A total of 121 datasets were used to calculate them.Equation ( 8) (the plug flow model) seemed to outperform the regression (Equation ( 9)) and Equation (2) (CSTR model).2), ( 8) and ( 9) to predict vacuum residue (540  In Figure 8, the training set is related to development of Equations ( 2), ( 8) and ( 9) with the data in Table S3.The test set was used to verify Equations ( 2), ( 8) and ( 9) with data not included in the process of equation development.
The statistical parameters used to evaluate the precision of conversion prediction are presented in Table 4.A total of 121 datasets were used to calculate them.Equation ( 8) (the plug flow model) seemed to outperform the regression (Equation ( 9)) and Equation (2) (CSTR model).The data in Table S3 along with the data in Tables S4-S9, which encompass the H-Oil operating conditions, conversions, yields as well as the characteristics of the H-Oil feeds and H-Oil VTBs were evaluated by the use of ICrA.The results of ICrA in terms of the µand υ-values with statistically meaningful relationships between the variables (HDAs, HDM, conversion, reactor temperatures, ATB TSE, PBFO TSP, product yields, and H-Oil VTB properties) are summarized in Tables S10 and S11.
The data in Tables S10 and S11 indicate that both de-asphaltization degrees (HDAs (C 5 ) and C 7 ), hydrodemetallization (HDM), and 540 • C+ conversion had a statistically significant positive association.This implies that the factors contributing to 540 • C+ conversion enhancement contribute to HDAs and HDM enlargement.The data in Tables S10 and S11 also indicate that the higher the molecular weight of the H-Oil feed, the lower the HDAs and HDM because they have a statistically meaningful negative association.This suggests that higher-molecular-weight H-Oil feeds may impede HDAs and HDM.It is evident from the data in Tables S10 and S11 that the 540 • C+ conversion and VTB density and the VTB C 5 and C 7 asphaltene contents had a statistically significant positive association, along with a statistically meaningful negative association with VTB molecular weight and T 50% .This means that the magnification of conversion is accompanied by the enhancement of VTB density and VTB C 5 and C 7 asphaltene contents as well as a reduction in VTB molecular weight.This can be explained by the cracking of the H-Oil 540 • C+ feed material, which leads to the production of lower-molecular-weight products like gas, naphtha, diesel, and VGO and leaves a lower-molecular-weight, higher aromaticity unconverted 540 • C+ material, denoted as VTB. Figure 9 illustrates the relationship between the VTB molecular weight and density and the 540 • C+ conversion.
suggests that higher-molecular-weight H-Oil feeds may impede HDAs and HDM.It is evident from the data in Tables S10 and S11 that the 540 °C+ conversion and VTB density and the VTB C5 and C7 asphaltene contents had a statistically significant positive association, along with a statistically meaningful negative association with VTB molecular weight and T50%.This means that the magnification of conversion is accompanied by the enhancement of VTB density and VTB C5 and C7 asphaltene contents as well as a reduction in VTB molecular weight.This can be explained by the cracking of the H-Oil 540 °C+ feed material, which leads to the production of lower-molecular-weight products like gas, naphtha, diesel, and VGO and leaves a lower-molecular-weight, higher aromaticity unconverted 540 °C+ material, denoted as VTB. Figure 9 illustrates the relationship between the VTB molecular weight and density and the 540 °C+ conversion.It is worth mentioning that along with the VTB density augmentation as a consequence of the increase in H-Oil conversion, the gas oils from the H-Oil, which make up the FCC feed, also followed the trend of VTB density enhancement, as illustrated in Figure 10.The data in Figure 10 indicate the same slope of the variation of the density of HVGO and LVGO with the VTB density alteration, whereas the HAGO density slope was substantially lower, suggesting that this H-Oil gas oil is less vulnerable to secondary cracking than HVGO and LVGO.Thus, it can be expected that the HAGO conversion observed in the FCC plant is less dependent on the severity of the H-Oil hydrocracker reaction than LVGO and HVGO.It is worth mentioning that along with the VTB density augmentation as a consequence of the increase in H-Oil conversion, the gas oils from the H-Oil, which make up the FCC feed, also followed the trend of VTB density enhancement, as illustrated in Figure 10.The data in Figure 10 indicate the same slope of the variation of the density of HVGO and LVGO with the VTB density alteration, whereas the HAGO density slope was substantially lower, suggesting that this H-Oil gas oil is less vulnerable to secondary cracking than HVGO and LVGO.Thus, it can be expected that the HAGO conversion observed in the FCC plant is less dependent on the severity of the H-Oil hydrocracker reaction than LVGO and HVGO.Another interesting fact observed in the data in Tables S10 and S11 is the relatively strong positive association (μ = 0.89; υ = 0.00) between the VTB sulfur content and the sulfur content of the H-Oil feed.The regression that relates the VTB sulfur content to the H-Oil feed sulfur content for this data set is given in Equation ( 10): = 0.5713 ×    − 0.3516 R 2 = 0.8704 (10) where: HOil VTB Sul = Sulfur content in the H-Oil VTB, wt.%; HOil Feed Sul = Sulfur content in the H-Oil feed, wt.%.
A statistically meaningful weak association between the VTB vanadium content and the VTB sulfur content can be seen in the data in Tables S10 and S11, suggesting that the  Another interesting fact observed in the data in Tables S10 and S11 is the relatively strong positive association (µ = 0.89; υ = 0.00) between the VTB sulfur content and the sulfur content of the H-Oil feed.The regression that relates the VTB sulfur content to the H-Oil feed sulfur content for this data set is given in Equation ( 10): HOil VTB Sul = 0.5713 × HOil Feed Sul − 0.3516 R 2 = 0.8704 (10) where: HOil VTB Sul = Sulfur content in the H-Oil VTB, wt.%; HOil Feed Sul = Sulfur content in the H-Oil feed, wt.%.A statistically meaningful weak association between the VTB vanadium content and the VTB sulfur content can be seen in the data in Tables S10 and S11, suggesting that the higher sulfur content of the H-Oil feeds are more difficult to demetallize.
Figure 11 presents the graphs of the function of the yields of H-Oil gas, naphtha, diesel, VGO, and VTB yields of conversion.It can be observed in the data in Figure 11 that the yields of gas, naphtha, and diesel, which are primary and secondary stable products [10], increased, while the VTB yield decreased as conversion increased.The yield of VGO, which is a primary unstable product, should exhibit a decrease as conversion increases beyond 70 wt.%,as shown in our recent investigation [23].The presence of a trend toward augmentation with conversion enhancement can be ascribed to the processing of VGO material in the H-Oil feed during the study.Using the data in Table 4 and Tables S4-S9 and employing statistical analyses, the following regression equations for the yields of H-Oil gas and diesel products were established:  Using the data in Tables 4 and S4-S9 and employing statistical analyses, the following regression equations for the yields of H-Oil gas and diesel products were established:  11b,e.
It should be noted here that none of the H-Oil catalyst properties not shown in this work due to confidential reasons exhibited any statistically meaningful relationships with sediment formation or fouling rates or with the conversion level or product yields.

Investigation of Fluid Catalytic Cracking (FCC)
The operating conditions and calculated parameters from the heat balance of the FCC unit for the cases studied are summarized in Table S12, and the conversion and product yields obtained fromt the FCC unit for cases studied are presented in Table S13.
Figure 12 plots the variation of the FCC product yields as conversion changed.Figure 12a indicates that yield of dry gas did not change with the conversion variation, supporting the thermal cracking nature of this product [45].The PPF, BBF, and gasoline yields show an increasing trend with increasing conversion, while the LCO and SLO yields exhibited a decreasing trend with conversion augmentation.The HCO yield demonstrated a difficult-to-distinguish trend of reduction as conversion increased.The reason for the relatively large dispersion in the data of the yields of FCC gas oils may be attributed to the various separation efficiencies of the FCC main fractionator, the FCC gasoline stabilizer, and the section of absorption and gas fractionation in this study.This is related to the production of different grades of finished automotive gasoline and diesel fuels.
Tables S14 and S15 present the µ and υ values of ICrA evaluation of the relationships between the FCC performance variables and the H-Oil performance variables.The data in these tables indicate that the H-Oil conversion level showed a statistically meaningful positive association with the density of H-Oil VGO that was processed in the FCC unit; with coke yield, ∆ coke; and with FCC HCO yield, as well as a statistically meaningful negative association with the FCC LCO yield.
This implies that increasing the H-Oil conversion is associated with an enhancement of the coke yield, ∆ coke, and the FCC HCO yield, along with a decrease in the FCC LCO yield.
The TSP of PBFO showed a statistically meaningful negative association with the FCC HCO yield, meaning that a decrease in the HCO yield is accompanied by an increase of PBFO TSP.
The content of FCC SLO in the H-Oil feed had a statistically meaningful negative association with FCC conversion and a statistically meaningful positive association with the FCC SLO yield.This suggests that an increase in FCC SLO content in the H-Oil feed is accompanied by a decrease in the FCC conversion and an increase inthe FCC SLO yield.This may also imply that the FCC SLO does not convert in the H-Oil unit, and the recycling of this material between the H-Oil and FCC unit occurs.This finding supports Marques et al.'s [46] observation that the FCC SLO material was refractory to hydroconversion under ebullated bed hydroconversion conditions.
work due to confidential reasons exhibited any statistically meaningful relationships with sediment formation or fouling rates or with the conversion level or product yields.

Investigation of Fluid Catalytic Cracking (FCC)
The operating conditions and calculated parameters from the heat balance of the FCC unit for the cases studied are summarized in Table S12, and the conversion and product yields obtained fromt the FCC unit for cases studied are presented in Table S13.
Figure 12 plots the variation of the FCC product yields as conversion changed.Figure 12a indicates that yield of dry gas did not change with the conversion variation, supporting the thermal cracking nature of this product [45].The PPF, BBF, and gasoline yields show an increasing trend with increasing conversion, while the LCO and SLO yields exhibited a decreasing trend with conversion augmentation.The HCO yield demonstrated a difficult-to-distinguish trend of reduction as conversion increased.The reason for the relatively large dispersion in the data of the yields of FCC gas oils may be attributed to the various separation efficiencies of the FCC main fractionator, the FCC gasoline stabilizer, and the section of absorption and gas fractionation in this study.This is related to the production of different grades of finished automotive gasoline and diesel fuels.The content of the H-Oil VGO in the FCC feed had a statistically meaningful negative association with FCC conversion and FCC gasoline and a statistically meaningful positive association with the FCC SLO yield.This means that the enhancement of the H-Oil VGO share in the FCC feed is accompanied by FCC conversion reduction, a drop in the gasoline yield, and an increase in the FCC SLO yield.
The FCC conversion was found to depend on only two variables for the studied cases, and these were the content of H-Oil VGO in the FCC feed and the content of FCC SLO in the H-Oil feed.The dependence of FCC conversion on content of H-Oil VGO in the FCC feed and content of FCC SLO in the H-Oil feed is shown in Equation (14).

Discussion of the Results of the Investigation of H-Oil Vacuum Residue Hydrocracking
Using the data shown in Tables 1 and 2 and in Figure 4, one can easily calculate the range of the variation of properties of the primary vacuum residual oil blend treated in the H-Oil hydrocracker.This, along with reactivity index, which presents the ratio of the vacuum residue hydrocracking reactivity to the reactivity of the Urals vacuum residue, is shown in Table 5.The reactivity index (RI) was calculated using the data shown in Table 2 and following the procedure described in our recent study [23].
Table 5. Range of variation of properties of the primary vacuum residual oil blend processed in the H-Oil unit during the study, along with the reactivity index.It is obvious from data in Table 5 that the reactivity of the primary vacuum residual oil blend treated in the H-Oil unit in the study was the same.Using the data in Table S4, the range of the variation of the RI of the H-Oil feed for the 40 cases explored was estimated, and it was established that the RI of the H-Oil feeds varied between 0.97 and 0.99, which means that practically no variation occurred in the H-Oil feed RI compared with that of the primary vacuum residual oil mixture.An interesting observation made in our earlier research [23] was that the vacuum residual oils with the same RI were expected to have the same sediment formation rate.Thus, we can expect that the processed H-Oil feeds for the 40 cases studied have the same sediment formation affinity.Therefore, we may deduce that the processed H-Oil feeds for the 40 cases studied should have the same sediment formation affinity.Indeed, the data in Figure 5 confirm that the ATB sediment content was almost the same, with some small variations.The curiosity in this case is that the ATB TSE was almost the same but at a different reaction temperature.Moreover, it is known that sedimentation increases exponentially with reaction temperature magnification [42,43].Such an increase in the sediment formation rate with an increase in reaction temperature cannot be seen in the data in Figure 5.This may be attributed to the suitable catalyst addition rate strategy and the application of a parallel mode of fresh catalyst addition, as discussed in our recent research [31].Sediment formation affinity during vacuum residue hydrocracking has been the subject of numerous investigations [36][37][38][39][40][41][42][43].Alonso et al. [47] classified residual oil hydrocracking feeds as high-, medium-, and low-sediment formation feeds.They reported that the high-sediment formation feeds originated from heavy and extra-heavy crude oils.The medium sediment formation feeds originated from medium crude oils, and the low sediment formation feeds originated from light crude oils [47].In our study, the crude oils processed were extra light, light, and medium, as apparent in the data in Table 1.Thus, light and medium sediment formation affinity may be expected during the processing of vacuum residual oils originating from the nine crude oils investigated in this study according to Alonso et al.'s [47] classification.However, Alonso et al. [47] also showed that the sediment formation rate increased exponentially as conversion was enhanced.The conversion level variation, as observed in data in Table S3, was between 51.8 and 91.7 wt.%, while the sediment content in the H-Oil ATB product (TSE), as seen in data in Figure 5, varied in a relatively narrow range between 0.02 and 0.46 wt.%.The colloidal instability index of the H-Oil VTB (CI I = Saturates+Asphaltenes Aromatics+Resins ) used to evaluate the colloidal stability and hence the sediment formation affinity [48], as indicated in data in Figure 13, did not show any trend twoard enhancement as conversion increased, supporting the observed lack of sediment augmentation with the increase in conversion.

Density at 15
uum residue hydrocracking has been the subject of numerous investigations [36][37][38][39][40][41][42][43].Alonso et al. [47] classified residual oil hydrocracking feeds as high-, medium-, and lowsediment formation feeds.They reported that the high-sediment formation feeds originated from heavy and extra-heavy crude oils.The medium sediment formation feeds originated from medium crude oils, and the low sediment formation feeds originated from light crude oils [47].In our study, the crude oils processed were extra light, light, and medium, as apparent in the data in Table 1.Thus, light and medium sediment formation affinity may be expected during the processing of vacuum residual oils originating from the nine crude oils investigated in this study according to Alonso et al.'s [47] classification.However, Alonso et al. [47] also showed that the sediment formation rate increased exponentially as conversion was enhanced.The conversion level variation, as observed in data in Table S3, was between 51.8 and 91.7 wt.%, while the sediment content in the H-Oil ATB product (TSE), as seen in data in Figure 5, varied in a relatively narrow range between 0.02 and 0.46 wt.%.The colloidal instability index of the H-Oil VTB ( = +ℎ + ) used to evaluate the colloidal stability and hence the sediment formation affinity [48], as indicated in data in Figure 13, did not show any trend twoard enhancement as conversion increased, supporting the observed lack of sediment augmentation with the increase in conversion.The increase in conversion, as shown in the data in Figure 9b, was accompanied by VTB density enhancement.As the density of VTB was correlated with the density of asphaltenes, as reported in our recent research [23], the density of asphaltenes is expected to grow with conversion magnification.Their higher density is related to their lower solubility, as discussed in our earlier study [49], and when blended with more paraffinic oil, they may facilitate the process of sediment formation.This hypothesis was verified with the experimental data on blending 50% of four H-Oil VTBs obtained at different reaction severities during the processing four different vacuum residual oil blends with 50% higher The increase in conversion, as shown in the data in Figure 9b, was accompanied by VTB density enhancement.As the density of VTB was correlated with the density of asphaltenes, as reported in our recent research [23], the density of asphaltenes is expected to grow with conversion magnification.Their higher density is related to their lower solubility, as discussed in our earlier study [49], and when blended with more paraffinic oil, they may facilitate the process of sediment formation.This hypothesis was verified with the experimental data on blending 50% of four H-Oil VTBs obtained at different reaction severities during the processing four different vacuum residual oil blends with 50% higher paraffinic FCCPT diesel.The data in Figure 14 clearly indicate that C 5 asphaltene solubility, judged by its density, was correlated with the sediment content of the blend of 50% H-Oil VTB/50% FCCPT diesel.
Therefore, one may conclude that during conversion increment the asphaltenes become denser and less soluble, but at lower saturate concentration in the VTB, they do not have the driving force to flocculate and precipitate.The conversion enlargement is associated with saturate reduction due to their cracking, and although the asphaltenes concentrate in the VTB, the colloidal instability index does not grow.This might be a plausible explanation for the observed lack of an increase in the sediment formation rate as conversion was enhanced.Actually, the stability of asphaltenes and its relationship to sediment formation is much more complex than the oversimplified explanation given above, and it has been shown to depend on asphaltene structure [50][51][52][53][54][55][56][57][58][59][60], the presence of impurities [50,58], free radicals [59], the content of metals [60,61], the solubility of asphaltenes, and the solubility power of the maltene fraction [62][63][64][65][66][67].An investigation by Dreillard et al. [24] also shoed that the hydrocracked residue sediment content at a high conversion depended not on asphaltene solubility, which was low, but on its concentration.The intricate matter of sediment formation is very well illustrated by the data shown in Figure 7.The data in Figure 7 indicate that the best sediment formation suppressor of the VTB obtained at a low severity (WABT of 405 • C) is the FCC HCO.This observation coincides with the findings of our resent research [18], indicating that FCC HCO was also the best sediment formation suppressor of the VTB obtained at a high severity (WABT of 430 • C) when Urals and Sib Light-derived vacuum residues were processed in H-Oil hydrocracker.However, when Basrah Heavy vacuum residue at a concentration of about 14% was processed along with Urals and Sib Light-derived vacuum residues, the FCC LCO became the best sediment formation suppressor of the VTB obtained at a high severity (WABT of 430 • C).Our first investigation dedicated to sediment formation in H-Oil hydrocracking demonstrated that among the high aromatic FCC gas oil diluents, there was no difference in their sediment formation inhibition when they were added to the H-Oil VTB [17].Marques et al. [46] announced that poly-condensed tri-, tetra-, and penta-aromatics containing in FCC SLO retard the sediment formation rate during Urals vacuum residue hydrocracking.If we take a look at the data on the diluents used in this study (Table S2), we can see that the aromatic ring index (ARI)-a correlation developed by Abutaqiya et al. [68] to relate the refractive index and molecular weight of nonpolar hydrocarbons to the average number of condensed aromatic ring structures-of FCC SLO is 3.6, which coincides with the one obtained by data reported by Marques et al. [62].The results reported in Figure 7 and Table S2, however, show that the FCC HCO with an ARI of 2.4 was much more efficient in the process of sediment formation retardation than the FCC SLO or the FCC LCO with an ARI of 1.5.Therefore, in this case, the polycondensed di-and tri-aromatics contained in the FCC HCO seem to be more efficient in sediment formation inhibition than the mono-and di-aromatics contained in the FCC LCO and the tri-, tetra-, and penta-aromatics contained in the FCC SLO.In a study by Tirado and Ancheyta [69], it was shown that mainly di-aromatics containing in the FCC LCO were the best suppressors of sediment formation during the hydroprocessing of a heavy crude oil.The results of Tirado and Ancheyta [69] indicated that the addition of 10 wt.% of FCC SLO led to more sediment formation than the hydrocracking of pure vacuum residue, which is a very convincing indicator of the very complex matter of sediment formation during vacuum residue hydrocracking.Understandably, such a process is very difficult to model and fully comprehend.
Another interesting observation was made during this study, and it was related to the content of metals in the hydrocracked vacuum residue (VTB) obtained at a low H-Oil reaction severity (WABT of 405 • C).The balance of aluminum over the H-Oil unit showed a higher amount of aluminum leaving the H-Oil unit with products than coming in with the vacuum residue feed.We hypothesized because of the lower viscosity of the oil in the H-Oil reactors, catalyst fines could be entrained with the hydrocracked vacuum residue (VTB) product.To verify this hypothesis, we measured the metal level in the VTB before and after filtration in the apparatus for the hot filtration test, and the results shown in Table S8 indicated that the metal level in the H-Oil VTB after filtration was considerably lower.The same procedure was applied to a sample of H-Oil VTB obtained at high reaction severity (430 • C), and the results showed no differences in metal level between the nonfiltrated and filtrated samples, suggesting that operation at a low severity is conducive to the entrainment of catalyst fines, with the hydrocracked vacuum residue increasing its metal content.
While the sediment formation in vacuum residue hydrocracking depends on many factors, conversion enhancement maily depends on two process variables: reaction time (throughput) and reaction temperature.The feedstock quality also affects conversion level, as discussed in [23].However, as mentioned at the beginning of this section, the quality of the mixed feed processing during this study did not vary substantially.Therefore, the process variables (reaction time and reaction temperature) were the only factors controlling the vacuum residue conversion.These two factors participate in the three models developed in this work to predict the conversion (see Equations ( 2), ( 8) and ( 9)).The data in Figure 8 and Table 4 suggest that the plug flow (Equation ( 8)) model had the highest accuracy of vacuum residue conversion in the LNB commercial ebullated bed H-Oil hydrocracker in comparison with the CSTR model (Equation ( 2)) and the linear regression (Equation ( 9)).The data in Figure 8 show some points with relatively a high dispersion for all three models.One may suppose that this high dispersion data concern the processing of feeds with different levels of quality.In fact, the main H-Oil feedstock was vacuum residue from Urals or a blend of Urals with Siberian Light, and the maximum content of other vacuum residues originating from different crudes was not higher than the 30% from the commissioning of the LNB H-Oil hydrocracker in 2015 [23].For this reason, the effect of feedstock on the relatively high deviation of observed conversion from the predicted one should not be considered.Of 121 datasets, 5 (4% of total data set) exhibited high deviation between the observed and predicted conversion.These data are related to the initial exploitation of the H-Oil hydrocracker, and therefore, a suspicious incorrect measurement of the quantity of feed and products might be the reason for this high deviation during the first stages of H-Oil exploitation.paraffinic FCCPT diesel.The data in Figure 14 clearly indicate that C5 asphaltene solubility, judged by its density, was correlated with the sediment content of the blend of 50% H-Oil VTB/50% FCCPT diesel.Therefore, one may conclude that during conversion increment the asphaltenes become denser and less soluble, but at lower saturate concentration in the VTB, they do not have the driving force to flocculate and precipitate.The conversion enlargement is associated with saturate reduction due to their cracking, and although the asphaltenes concentrate in the VTB, the colloidal instability index does not grow.This might be a plausible explanation for the observed lack of an increase in the sediment formation rate as conversion was enhanced.Actually, the stability of asphaltenes and its relationship to sediment formation is much more complex than the oversimplified explanation given above, and it has been shown to depend on asphaltene structure [50][51][52][53][54][55][56][57][58][59][60], the presence of impurities [50,58], free radicals [59], the content of metals [60,61], the solubility of asphaltenes, and the solubility power of the maltene fraction [62][63][64][65][66][67].An investigation by Dreillard et al. [24] also shoed that the hydrocracked residue sediment content at a high conversion depended not on asphaltene solubility, which was low, but on its concentration.The intricate matter of sediment formation is very well illustrated by the data shown in Figure 7.The data in Figure 7 indicate that the best sediment formation suppressor of the VTB obtained at a low severity (WABT of 405 °C) is the FCC HCO.This observation coincides with the findings of our resent research [18], indicating that FCC HCO was also the best sediment formation suppressor of the VTB obtained at a high severity (WABT of 430 °C) when Urals and Sib Light-derived vacuum residues were processed in H-Oil hydrocracker.However, Figure 9 indicates that molecular weight of hydrocracked vacuum residue fell and the density grew as the conversion increased, which is easily explained by the cracking of the aliphatic (paraffinic and naphthenic) moieties, leaving a lower-molecular-weight, higher-aromatic unconverted hydrocracked vacuum residue.This observation is in line with the findings reported by Dreillard et al. [24].
The data in Figure 10 suggest not only that the hydrocracked vacuum residue cracked in the H-Oil process, making it denser and more aromatic as reaction proceeded, but also that the three gas oil fractions: HVGO, LVGO, and HAGO were subject to secondary cracking, making them denser and more aromatic.Surprisingly, for the studied period, no statistically meaningful relationships between H-Oil VGO density and FCC conversion were found, irrespective of the wider range of H-Oil conversion variation and the consequent alteration of H-Oil gas oil densities.The scope of the alteration of the density of H-Oil HAGO, LVGO, and HVGO for the studied period estimated from the VTB density data shown in Table S6 and the regression equations embedded in Figure 10 is presented in Table 6.The H-Oil product yields variation during the study shown in Figure 11 indicate that they depend mainly on conversion level but also on the content of recycled material, FCC SLO, material boiling at up to 360 • C, and VGO in the H-Oil feed.Thus, by varying the content of these H-Oil feed components to different extents, a distinct product yield slate can be obtained.

Discussion the Results from the Investigation of the Fluid Catalytic Cracking (FCC)
The data in Table 3 exhibit a relatively wide variation in the quality of the FCC feedstock of SRVGOs originating from the nine crude oils processed in the LNB refinery in the study.Irrespective of the wide range of the variation of density (between 0.884 and 0.949 g/cm 3 ) and the related content of saturates plus mononuclear aromatics (between 61 and 83.1 wt.%, see Table 3), which determines the FCC feedstock crackability [70], the density of the blended SRVGOs estimated on the basis of the data in Tables 1 and 3 and the data in Figure 1 shows a very narrow variation between 0.913 and 0.922 g/cm 3 .This range of the variation of density corresponded to variation in the contents of saturates plus mononuclear aromatics between 70.3 and 73.3 wt.%, which is within the uncertainty of its determination.Therefore, the variation of the quality of the SRVGO as an FCC feed should be considered negligible.
The dispersion of the data in Figure 12 showing the relationship between the yields of dry gas, the propane-propylene fraction (PPF), the butane-butylene fraction (BBF), gasoline, LCO, HCO, and SLO and the FCC conversion is considerably higher than that observed in the H-Oil investigation (for comparison, see Figure 11).The reason for this high dispersion is the variable efficiency of the separation columns in the FCC unit, which was affected by the high natural gas prices, requiring maximum extraction of dry gas from the FCC unit while penalizing the quality of separation between dry gas and PPF.The production of different grades of gasolines with different specifications also affected the efficiency of the separation between gasoline and BBF and between gasoline and LCO.
In our previous research [23], we found that severity of the H-Oil unit affected FCC conversion and the regenerator temperature, whereas in this study, the only relationship between the quality of H-Oil gas oil, expressed by its density, and the FCC unit performance was the regenerator temperature (see Figure 15).In our previous research [23], we found that severity of the H-Oil unit affected FCC conversion and the regenerator temperature, whereas in this study, the only relationship between the quality of H-Oil gas oil, expressed by its density, and the FCC unit performance was the regenerator temperature (see Figure 15).Probably, the reconstruction of the FCC unit implemented in late 2021 and discussed in our recent research [26] affected the performance of the FCC unit in a way that was less dependent on the severity of the H-Oil operation.Now, as regression Equation (13) shows, the FCC conversion for the operation condition window displayed in Table S13   Probably, the reconstruction of the FCC unit implemented in late 2021 and discussed in our recent research [26] affected the performance of the FCC unit in a way that was less dependent on the severity of the H-Oil operation.Now, as regression Equation (13) shows, the FCC conversion for the operation condition window displayed in Table S13 is correlated with the amount of H-Oil VGO processed in the FCC unit and the amount of FCC SLO processed in the H-Oil unit.The optimization of the quantity of the FCC SLO processed in the H-Oil unit, considering its inhibiting impact on the sediment formation rate, can deliver higher FCC conversion while keeping the sediment formation rate and fouling in H-Oil hydrocracker within reasonable limits.

Conclusions
Nine distinct crude oils belonging to the extra light, light, and medium petroleum crude types were refined as diverse mixtures in the LNB refinery to study the combined performance of the two most profitable processes: FCC of VGO and VR ebullated bed H-Oil hydrocracking.Although the qualities of the individual vacuum residues and vacuum gas oils originating from the nine processed crude oils significantly differed, those of the blended SRVRs and SRVGOs were almost the same.The difference in the quality of the H-Oil feed came from the distinct amount of co-processed FCC SLO, recycled PBFO, and VGO with SRVR.The variation in the quality of the FCC feed came from the diverse amount and quality of co-processed H-Oil VGO with the SRVGO.For the studied range of operating conditions and feed quality, the H-Oil conversion level was found to depend on only two process variables: feed throughput and WABT.Three models were developed using CAS Maple 2023 Academic Edition Software: regression, a CSTR model, and a plug flow reactor model to predict the H-Oil vacuum residue conversion.The conversion level was found to be modeled best by the plug flow reactor model with fitted parameters of k 0 = 4.42875 × 10 15 ; E A = 215 kJ/mol; n = 1.59.The H-Oil product yields were found to depend on the conversion level and content of FCC SLO, recycled PBFO, the material boiling below 360 • C, and VGO in the H-Oil feed.The sediment formation rate in the H-Oil unit did not exhibit the typical exponential increase as conversion increased irrespective of the asphaltene density enhancement and its solubility reduction.It was found that the asphaltene density was correlated with the sediment content of the blend of 50% H-Oil VTB/50% high paraffinic FCCPT diesel.The molecular weight of the hydrocracked vacuum residue decreased while its density increased as the vacuum residue conversion was enhanced in the H-Oil hydrocracker.The density of the H-Oil gas oils HAGO, LVGO, HVGO followed the trend of the H-Oil VTB density augmentation with conversion magnification.It was found that during the low severity mode of H-Oil operation, the entrainment of catalyst fines with the hydrocracked vacuum residue product occurred, leading toa higher metal level in this product.
Among all the variables of the studied window of operating conditions of the FCC unit, only the H-Oil VGO content in the FCC feed and the FCC SLO share in the H-Oil feed were found to control the vacuum gas oil conversion in the FCC process.The empirical models developed in this research are intended to be used for performance improvement of both FCC and H-Oil units, including the amelioration of energy consumption and, subsequently, CO 2 emissions reduction.and H-Oil VTB properties; Table S11: υ-Value of the ICrA of relationships of HDAs, HDM, conversion, reactor temperatures, ATB TSE, PBFO TSP, product yields, and H-Oil VTB properties; Table S12: Operating conditions and calculated parameters from the heat balance of the FCC unit for the studied cases; Table S13: Conversion and product yields obtained from the FCC unit for the studied cases; Table S14: µ-Value of the ICrA of relationships of H-Oil performance variables and FCC performance variables; Table S15: υ-Value of the ICrA evaluation of relationships of H-Oil performance variables and FCC performance variables.

Conflicts of Interest:
The authors declare no conflict of interest.

Figure 2 .
Figure 2. Process scheme of LNB FCC feed hydrotreating and FCC reactor-regenerator sections.

Figure 3 .
Figure 3.A simplified layout of the heavy oil technological chain in LUKOIL Neftohim Burgas refinery.

Figure 3 .
Figure 3.A simplified layout of the heavy oil technological chain in LUKOIL Neftohim Burgas refinery.

Figure 3 .
Figure 3.A simplified layout of the heavy oil technological chain in LUKOIL Neftohim Burgas refinery.

Figure 4 .
Figure 4. Variation of individual crude oil contents in the refined petroleum mixture in the LNB refinery during the study.

Figure 4 .
Figure 4. Variation of individual crude oil contents in the refined petroleum mixture in the LNB refinery during the study.

Figure 5 .
Figure 5. Variation of ATB TSE, and PBFO TSP along with the WABT of both H-Oil hydrocracking reactors in the study.

Figure 5 .
Figure 5. Variation of ATB TSE, and PBFO TSP along with the WABT of both H-Oil hydrocracking reactors in the study.

Figure 6 .
Figure 6.Variation of atmospheric tower skin temperature and the difference between the ATB product temperature and the atmospheric tower bottom skin temperatures measured in diverse radial positions (denoted as DT1, DT2, DT3, and DT4).

Figure 6 .
Figure 6.Variation of atmospheric tower skin temperature and the difference between the ATB product temperature and the atmospheric tower bottom skin temperatures measured in diverse radial positions (denoted as DT1, DT2, DT3, and DT4).

Processes 2023 , 27 Figure 6 .
Figure 6.Variation of atmospheric tower skin temperature and the difference between the ATB product temperature and the atmospheric tower bottom skin temperatures measured in diverse radial positions (denoted as DT1, DT2, DT3, and DT4).

FCC 27 Figure 7 .
Figure 7. TSE and TSP of blends of H-Oil VTB obtained at very low severity (405 °C) with FCC HCO (a), FCC LCO (b), FCC SLO (c), and FCCPT diesel (d) as a function of diluent content in the blend.

FCCFigure 7 .
Figure 7. TSE and TSP of blends of H-Oil VTB obtained at very low severity (405 • C) with FCC HCO (a), FCC LCO (b), FCC SLO (c), and FCCPT diesel (d) as a function of diluent content in the blend.

Figure 9 .
Figure 9. Variation of H-Oil VTB molecular weight (a) and density (b) with conversion alteration.

Figure 9 .
Figure 9. Variation of H-Oil VTB molecular weight (a) and density (b) with conversion alteration.

Figure 10 .
Figure 10.Relationship between H-Oil VTB density and the densities of the three H-Oil gas oils (HVGO, LVGO, and HAGO).

Figure 10 .
Figure 10.Relationship between H-Oil VTB density and the densities of the three H-Oil gas oils (HVGO, LVGO, and HAGO).

Figure 15 .
Figure 15.Relationship between the H-Oil VGO density and the FCC regenerator dense bed temperature.

Figure 15 .
Figure 15.Relationship between the H-Oil VGO density and the FCC regenerator dense bed temperature.

Author Contributions:
Conceptualization, D.S. and I.S.; methodology, V.T.; software, D.D.S., S.N., and V.B.; validation, D.P.; formal analysis, S.V.; investigation, V.T. and D.S.; resources, K.A.; data curation, S.V.; writing-original draft preparation, D.S.; writing-review and editing, D.S.All authors have read and agreed to the published version of the manuscript.Funding: This was funded by the European Union's NextGenerationEU through the National Recovery and Resilience Plan of the Republic of Bulgaria project № BG-RRP-2.004-0002,"BiOrgaMCT".Data Availability Statement: Data are contained within the article.

Table 1 .
Characteristics of investigated crude oils.

Table 2 .
Characteristics of vacuum residual (VR) oils distilled from the studied crude oils by laboratory atmospheric and vacuum distillations.

Table 3 .
Characteristics of vacuum gas oil (VGO) fractions distilled from the investigated crude oils.

Table 3 .
Characteristics of vacuum gas oil (VGO) fractions distilled from the investigated crude oils.

Table 4 .
Statistical analysis of Equations ( • C+) conversion in the LNB H-Oil hydrocracker.

Table 6 .
Range of the variation of the density of the H-Oil gas oils in the vacuum residue conversion range of 51.8-91.7 wt.%.