Quest for the Co-Pyrolysis Behavior of Rice Husk and Cresol Distillation Residue: Interaction, Gas Evolution and Kinetics

: With the tremendous prosperity of industry, more and more hazardous waste is discharged from industrial production processes. Cresol distillation residue is a typical industrial hazardous waste that causes severe pollution without proper treatment. Herein, the co-pyrolysis of rice husk and cresol distillation residue was studied using thermogravimetry–mass spectrometry and kinetic studies. The Coats and Redfern method was employed to calculate the activation energy. The results indicated that the pyrolysis process of cresol distillation residue and RH/CDR (Rice Husk and Cresol Distillation Residue) blends can be divided into four stages and three stages for RH. The introduction of RH not only improved the thermo-stability of cresol distillation residue at a low temperature but also reduced the activation energy of the blends. The activation energy was the lowest when the proportion of rice husk in the blend was 60%. The main gaseous pyrolysis products included CH 4 , H 2 O, C 2 H 2 , CO 2 , C 3 H 6 and H 2 . There existed an unusual combination of synergistic and inhibitive interactions between RH and cresol distillation residue, respectively, within different temperature ranges. The synergistic interaction decreased the reaction’s activation energy, whereas the inhibitive interaction reduced the emission of main gaseous products, such as CH 4 and CO 2 . It was concluded that the addition of RH was conducive to improving the pyrolytic performance of cresol distillation residue and the resource utilization of cresol distillation residue.


Introduction
Distillation processes dominate 60% of separation in the chemical industry [1], but there are ca. 2.5 million tons of distillation residues produced in China every year. Distillation residues have been included in the national hazardous waste lists of different countries [2]. In general, distillation residue is mainly treated by landfill or incineration approaches [3][4][5]. The landfill of distillation residues generates large amounts of leachate that severely contaminate the soil and underground water [6][7][8][9]. The incineration of distillation residues is a high-energy consumption process that also produces severe secondary pollution [10,11]. The resource utilization of distillation residues for the production of value-added products seems to be an environment-benign approach [12], but only a handful of papers exist regarding the conversion of distillation residues to diesel and lubricating oil [13].
At present, varieties of biomass with clean and renewable characteristics have been used for the resource utilization of solid waste, including sewage sludge, food waste and municipal solid waste, through co-pyrolysis technology [14][15][16][17]. There is an inhibitive or synergistic interaction between biomass and solid waste during thermal treatment. We note that only one paper reported an inhibitive interaction between solid waste and biomass, in which the undecomposed lignite particles prevented the release of volatile matters in solid waste derived from refining and chemical wastewater at a lower temperature [18]. In contrast, synergistic interactions have been extensively studied. Several researchers have reported a synergistic interaction between sewage sludge and biomass that reduced the release of gaseous sulfur substances and NOx [19][20][21]. The structure of pyrolysis products can be optimized through synergistic interactions, as exemplified by the improved surface area of combustion ashes from textile dyeing sludge [22]. In addition, the synergistic interactions of blends have also resulted in higher reactivity and better combustion performance. For example, the synergistic interactions between textile dyeing sludge and microalgae improve the combustion performance of textile dyeing sludge because the density of its blends is larger than single microalgae [23]. We note that the co-pyrolysis of industrial distillation residues with biomass remains relatively underexplored. There are fewer than five studies focusing on co-pyrolysis of biomass with the distillation residue from lab-scale bio-oil production [24][25][26][27]. The interaction between industrial distillation residues and biomass remains relatively unexplored.
In this contribution, we reported an unexpected synergistic effect that combined the high-temperature inhibitive and low-temperature synergistic processes in a sequential manner during the co-pyrolysis of cresol distillation residue and rice husk. Cresol distillation residue is a typical industrial waste from the production of p-cresol, which has been widely used for the synthesis of pharmaceuticals, herbicides, antioxidants and dyes [28]. It is estimated that 10-12 tons of cresol residue can be produced for every 100 tons of pcresol [29]. This work aims to investigate the interactions and product characteristics during co-pyrolysis of cresol distillation residue (CDR) and rice husk (RH) at various mixing ratios through thermogravimetric analysis. The interactions between cresol distillation residue and rice husk were investigated using the deviation of weight loss TG (∆W) between the calculated and experimental values in detail. TG coupled with mass spectrometry (TG-MS) enables the tracing of thermal reactions and the characterization of evolved gases. The kinetic parameters and apparent activation energy during thermal decomposition were calculated using the Coats and Redfern model. The elucidation of interactions between RH and CDR during the co-pyrolysis process is likely to provide scientific support for the effective utilization of CDR and to reduce related environmental hazards.

Materials
The cresol distillation residue was supplied by a local company in Nantong, and rice husks were purchased from an agricultural product store. RH was mixed with CDR in a ratio of 20-80% to examine the effect of biomass additives on thermal decomposition during the co-pyrolysis process. Residue-biomass mixtures were prepared with a biomass content of 20%, 40%, 60% and 80%. All samples were oven-dried at 105 ± 5 • C for 24 h to remove water before pyrolysis.
The elemental content and chemical features of raw materials are exhibited in Table 1. The ash content, volatile matter and moisture were measured following the Chinese coal industry method (Chinese standard methods, GB/T 212-2008). The elemental analysis was performed using an elemental analyzer (Euro Vector EA3000, NETZSCH, Italy). The element ratio of H/C was adopted to characterize the amount of CO 2 released during pyrolysis [30]: the higher the H/C, the lower the level of CO 2 per unit of energy produced [31]. In this work, it was demonstrated that in comparison to cresol distillation residue, rice husk has potentially lower CO 2 emissions. The elemental content of Ca, Fe, Al, K, P, Cl and Si were determined using X-ray fluorescence ( Table 2).  CDR was mixed with RH ground into a particle size of less than 100 mesh via vigorous stirring, with an RH weight percentage of 20-80 wt.%. Pyrolysis of all samples was performed on STA 449-QMS 403 TG-MS analyzers. About 10 ± 0.5 mg samples with particle sizes less than 100 mesh were heated from room temperature to 1000 • C with a heating rate of 20 • C/min and a N 2 flow rate of 20 mL/min. The generated gaseous products were monitored, including H 2 (m/z = 2), CH 4 (m/z = 16), H 2 O (m/z = 18), C 2 H 2 (m/z = 26), C 3 H 6 (m/z = 42) and CO 2 (m/z = 44), according to the database of the National Institute of Standards and Technology (NIST). All the pyrolysis test conditions were repeated two or three times, and the average data were taken to ensure repeatability.

Kinetic Analysis
We selected the Coats and Redfern model to calculate kinetic parameters based on TG analysis [32,33].
The conversion fraction is the mass fraction of a decomposed sample [34,35], which is described by: where x is the conversion extent of samples, W 0 = initial mass, W t = instantaneous mass and W ∞ = final mass. Isothermal reaction rate: A constant heating rate: Thus, Equation (2) is also equal to the following equation: The integral of Equation (4): First, perform partial integration on the right side of Equation (5); then, ignore the higher-order terms to obtain the following expression: Equations (6) and (7) can be written as a logarithm: Because 2RT/E 1, we simplify Equations (8) and (9): It is assumed that the pyrolysis reactions follow the first-order reaction kinetics at n = 1, which gives: For easy calculating, Equation (12) can be changed to the straight-line formula: x = instantaneous conversion ratio, T = absolute temperature (K), β = heating rate ( • C/min), R = Gas constant (J·mol −1 ·K −1 ), A 0 = pre-exponential factor (min −1 ) and E = activation energy (K·J·mol −1 ).

Individual Samples
Thermogravimetric (TG) and differential thermal gravity (DTG) curves showed a difference in thermal behaviors between RH and CDR ( Figure 1). RH consists of (hemi)cellulose, lignin and other organic minorities [36,37] so that the weight loss of RH can be divided into three main stages ( Figure 1a). The release of adsorbed gases and water vapor occurred from 50 to 218 • C, followed by the thermal decomposition of volatiles from (hemi)cellulose with the range of 218 to 409 • C (weight loss = 57.3%). The final stage is the slow decomposition of lignin and carbonaceous residue with a weight loss of 14.9%. In comparison, the pyrolysis of CDR was divided into four individual stages ( Figure 1b). Water evaporation occurred before 129 • C, after which the pyrolysis of (hemi)celluloses, phenols and ethers occurred until 510 • C. The weight loss from 510 • C to 750 • C corresponded to the charring reaction of lignin and the remaining hydrocarbon, followed by the thermal decomposition of inorganic substances at a very low rate [38]. As the DTG curves show, the maximum weight loss rate of RH (19.16 wt.%/min) is much higher than CDR (9. 4 wt.%/min) (Table 3), indicating the higher reactivity of RH than CDR during the co-pyrolysis process. Compared with RH, the peak temperature of maximum weight loss of CDR was lower than that of RH, indicating that CDR was more unstable during the co-pyrolysis process. Therefore, RH exhibited better pyrolysis performance than CDR, and it is estimated that the blends of RH with CDR enhance the pyrolysis process of CDR.
Water evaporation occurred before 129 °C, after which the pyrolysis of (hemi)celluloses, phenols and ethers occurred until 510 °C. The weight loss from 510 °C to 750 °C corresponded to the charring reaction of lignin and the remaining hydrocarbon, followed by the thermal decomposition of inorganic substances at a very low rate [38]. As the DTG curves show, the maximum weight loss rate of RH (19.16 wt.%/min) is much higher than CDR (9. 4 wt.%/min) ( Table 3), indicating the higher reactivity of RH than CDR during the co-pyrolysis process. Compared with RH, the peak temperature of maximum weight loss of CDR was lower than that of RH, indicating that CDR was more unstable during the copyrolysis process. Therefore, RH exhibited better pyrolysis performance than CDR, and it is estimated that the blends of RH with CDR enhance the pyrolysis process of CDR.  Ti is the initial decomposition temperature; Tf is the final decomposition temperature; Tmax is the first peak temperature; DTGmax is the maximum weight loss rate; Mf is the residue mass.
We calculate the error with Mf as the reference:   T i is the initial decomposition temperature; T f is the final decomposition temperature; T max is the first peak temperature; DTG max is the maximum weight loss rate; M f is the residue mass.
The uncertainty values u 1 , u 2 , u 3 , u 4 , u 5 and u 6 are introduced by the measurement error of RH, 80RH, 60RH, 40RH, 20RH and CDR [39]. u is written as: where δ is error, and k = √ 3. We calculate the error with M f as the reference: The larger the uncertainty value is, the more unstable the sample is. The most unstable sample is CDR, which is consistent with the above analysis.

Blend Samples
The TG/DTG curves of the blends of RH and CDR can be roughly divided into four stages, similar to CDR, but the beginning temperature (T i ), maximum decomposition rate (T max ) and ending temperature (T f ) of the blends are totally different (Figure 2a,b, Table 3). The addition of RH to CDR reduced the heat transfer from the surface to the core of the blend samples, giving rise to the shift of both blends' T i and T f to higher temperatures [40]. The addition of RH significantly reduced the maximum rate of weight loss (DTG max ) and residual mass (M f ) of the blend samples from 20RH to 60RH. We noted that the DTG max peak gradually moved toward the high-temperature zone with the increase in RH content, indicating the addition of RH reduced the reactivity of the mixtures. However, when the proportion of RH is 80%, the DTG curve of 80RH/CDR is similar to pure RH. In Figure 2c, the conversion of all samples occurred from 100 to 900 • C. The conversion curve of the mixture almost coincides with the RH curve from 372 to 472 • C. The pyrolysis of all samples was delayed when the temperature rose above 472 • C, corresponding to the highly reduced value of the second DTG peak in Figure 2b. stages, similar to CDR, but the beginning temperature (Ti), maximum decomposition rate (Tmax) and ending temperature (Tf) of the blends are totally different (Figure 2a,b, Table 3). The addition of RH to CDR reduced the heat transfer from the surface to the core of the blend samples, giving rise to the shift of both blends' Ti and Tf to higher temperatures [40]. The addition of RH significantly reduced the maximum rate of weight loss (DTGmax) and residual mass (Mf) of the blend samples from 20RH to 60RH. We noted that the DTGmax peak gradually moved toward the high-temperature zone with the increase in RH content, indicating the addition of RH reduced the reactivity of the mixtures. However, when the proportion of RH is 80%, the DTG curve of 80RH/CDR is similar to pure RH. In Figure 2c, the conversion of all samples occurred from 100 to 900 °C. The conversion curve of the mixture almost coincides with the RH curve from 372 to 472 °C. The pyrolysis of all samples was delayed when the temperature rose above 472 °C, corresponding to the highly reduced value of the second DTG peak in Figure 2b.

Interactive Effects Analysis
For the purpose of investigating the interactions between RH and CDR during copyrolysis, the weight loss deviations between experimental and calculated values within the whole temperature range were calculated according to the following equation [41,42]: where W RH and W CDR represented the weight loss (TG) of each sample. x 1 and x 2 were the proportion of RH and CDR in the blend, respectively. W calculated is the sum of the component weight based upon its fraction at certain temperatures [43][44][45]. ∆W refers to the deviation of weight loss of the blend according to the TG curves, which could be used as an indicator of the interaction degree. There is no interaction if the value of ∆W is 0 [46]. Three stages of co-pyrolysis could be identified from Figure 3. The negative ∆W value reflects the synergistic effect between RH and CDR at temperatures below 374 • C due to the catalytic effect of alkali and alkaline earth metals in the rice husk (Table 2) [47,48]. In comparison, the calculated TG curves were below the experimental TG curves above 455 • C, demonstrating an inhibitive effect between RH and CDR during pyrolysis. The inhibitive effect is ascribed to the adherence of CDR pyrolysis products to the blend's surface, which, in turn, prevented further volatilization and attenuated the heat/mass transfers [49]. Moreover, it is the first time a perfect match between the calculated TG curves and the experimental TG curves between 374 and 455 • C has been reported, indicating the complete degradation of all volatiles generated from the initial pyrolysis of blends in this temperature range.
indicator of the interaction degree. There is no interaction if the value of ΔW is 0 [46].
Three stages of co-pyrolysis could be identified from Figure 3. The negative ΔW value reflects the synergistic effect between RH and CDR at temperatures below 374 °C due to the catalytic effect of alkali and alkaline earth metals in the rice husk (Table 2) [47,48]. In comparison, the calculated TG curves were below the experimental TG curves above 455 °C, demonstrating an inhibitive effect between RH and CDR during pyrolysis. The inhibitive effect is ascribed to the adherence of CDR pyrolysis products to the blend's surface, which, in turn, prevented further volatilization and attenuated the heat/mass transfers [49]. Moreover, it is the first time a perfect match between the calculated TG curves and the experimental TG curves between 374 and 455 °C has been reported, indicating the complete degradation of all volatiles generated from the initial pyrolysis of blends in this temperature range.

TG-MS Analysis
The ion fragments of the main gaseous products generated during pyrolysis were monitored by TG-MS analysis ( Figure 4). CH4, H2O and CO2 were major products during the pyrolysis process. CH4 has three ionic strength peaks in the range of 100-309 °C, 309-478 °C and 478-900 °C, respectively. The synergistic stage from 100 °C to 309 °C is associated with the cracking of side chains of aliphatic hydrocarbons. The conversion of longchained aromatic groups and alkyl groups to methane then occurred from 309 °C to 478 °C [50]. In the inhibitive stage, CH4 is mainly derived from the conversion of methoxyl groups in lignin after 478 °C [51,52]. The CH4 peak intensity in the synergistic zone was lower than that in the inhibitory zone, which is mainly attributed to the speed of weight loss. However, the amount of CH4 released in both zones decreased after adding rice husk.

TG-MS Analysis
The ion fragments of the main gaseous products generated during pyrolysis were monitored by TG-MS analysis (Figure 4). CH 4 , H 2 O and CO 2 were major products during the pyrolysis process. CH 4 has three ionic strength peaks in the range of 100-309 • C, 309-478 • C and 478-900 • C, respectively. The synergistic stage from 100 • C to 309 • C is associated with the cracking of side chains of aliphatic hydrocarbons. The conversion of long-chained aromatic groups and alkyl groups to methane then occurred from 309 • C to 478 • C [50]. In the inhibitive stage, CH 4 is mainly derived from the conversion of methoxyl groups in lignin after 478 • C [51,52]. The CH 4 peak intensity in the synergistic zone was lower than that in the inhibitory zone, which is mainly attributed to the speed of weight loss. However, the amount of CH 4 released in both zones decreased after adding rice husk. H 2 O evolution can also be divided into three stages. H 2 O peaks in the synergistic region (<317 • C) are attributed to the loss of cellular water and the hydration of the phenolic hydroxyl groups from chemically bonded water and distillation residue [53]. Then, the thermal decomposition of light-volatile components occurred from 317 to 470 • C. The release of H 2 O in the inhibitive region (>470 • C) resulted from the binding of some free radical groups, such as hydroxide radicals, and oxygen ions, as well as the decomposition of O-contained functional groups (especially hydroxyl groups) [54]. During the whole pyrolysis process, the amount of H 2 O released from blends 60RH and 80RH was always higher, as opposed to 20RH and 40RH.
The concentration of CO 2 initially increased in the synergistic region (300-450 • C) because of the breaking of aromatic moieties and carboxyl groups. The synergistic effect led to the rapid release of CO 2 . Then, the second peak of CO 2 was observed between 450 • C and 700 • C, which resulted from the degradation of carbonyl compounds and oxygenated compounds with high thermal stability [55]. When the temperature reached the inhibitive region (>700 • C), the decrease in CO 2 indicated the decomposition of a small amount of CaCO 3 [56]. However, the CO 2 emission of blends decreased compared with CDR. The CO 2 emission was the lowest at 60RH due to the inhibitive interaction.
The H 2 emission occurred continuously over the temperature range of 300-900 • C (single peak) for each sample. The degradation of hydrogen-rich compounds occurred at about 300 • C, and the condensation of (hydro)aromatic compounds or the decomposition of heterocyclic compounds were detected above 600 • C [53,57]. The hydrogen release of CDR is higher than RH owing to the higher hydrogen content in CDR (Table 1). For blends, the peak values of H 2 occurred in the order of 40RH, 20RH, 0RH, 60RH, 80RH and 100RH, which differed from the regular sequence of the RH-addition ratio. Additionally, 20RH and 40RH contributed to the promotion of H 2 production, while other mixtures were just the opposite. H 2 yield was highest when the RH-addition ratio was 40 wt% (40RH).  H2O evolution can also be divided into three stages. H2O peaks in the synergistic region (<317 °C) are attributed to the loss of cellular water and the hydration of the phenolic hydroxyl groups from chemically bonded water and distillation residue [53]. Then, the thermal decomposition of light-volatile components occurred from 317 to 470 °C. The release of H2O in the inhibitive region (>470 °C) resulted from the binding of some free radical groups, such as hydroxide radicals, and oxygen ions, as well as the decomposition of O-contained functional groups (especially hydroxyl groups) [54]. During the whole pyrolysis process, the amount of H2O released from blends 60RH and 80RH was always higher, as opposed to 20RH and 40RH.
The concentration of CO2 initially increased in the synergistic region (300-450 °C) because of the breaking of aromatic moieties and carboxyl groups. The synergistic effect led to the rapid release of CO2. Then, the second peak of CO2 was observed between 450 °C and 700 °C, which resulted from the degradation of carbonyl compounds and oxygenated compounds with high thermal stability [55]. When the temperature reached the in- Aliphatic hydrocarbons C n H m (n ≥ 2) can be generated by two pathways. One is the decomposition of macromolecular components, including long-chain and branched-chains paraffin in CDR, into small molecules. Another is the degradation of hydroaromatic groups, polymethylene and aliphatic bridges (e.g., n-fatty acis). Alkyne (C 2 H 2 ) and alkene (C 3 H 6 ) have strong ionic peaks within the temperature range of 200-700 • C.

Kinetic Analysis
All kinetic parameters were calculated via the Coats and Redfern method within the temperature ranges from 190 to 380 • C and from 260 to 400 • C for CDR and RH, respectively. According to Equation (12), a series of linear-fitting curves were obtained and plotted in  Figure 5. The activation energy E a and the pre-exponential factor A 0 were thereby generated for various blend ratios ( Table 4). The activation energies for RH, CDR and their blends lie in the range of 15-25 kJ/mol. The linear correlation coefficient R 2 of 0.95 indicates the reasonable fitting of the first-order reaction model. As exhibited in Table 4, the E a values of RH and CDR were 21.85 kJ/mol and 24.00 kJ/mol, respectively. The E a values of blends with different proportions were lower than those of RH or CDR, demonstrating that the existence of RH could reduce the energy required for CDR pyrolysis. This result also verified that a higher value of activation energy indicates a slower reaction. As the activation energy represents the critical energy required to initiate a reaction [58], the 60 RH blend, with the lowest activation energy among these blended samples, is recommended. It is consistent with the positive synergistic effect induced by the addition of RH so that the optimal mixing ratio of RH to CDR is 3:2. matic groups, polymethylene and aliphatic bridges (e.g., n-fatty acis). Alkyne (C2H2) and alkene (C3H6) have strong ionic peaks within the temperature range of 200-700 °C.

Kinetic Analysis
All kinetic parameters were calculated via the Coats and Redfern method within the temperature ranges from 190 to 380 °C and from 260 to 400 °C for CDR and RH, respectively. According to Equation (12), a series of linear-fitting curves were obtained and plotted in Figure 5. The activation energy Ea and the pre-exponential factor A0 were thereby generated for various blend ratios ( Table 4). The activation energies for RH, CDR and their blends lie in the range of 15-25 kJ/mol. The linear correlation coefficient R 2 of 0.95 indicates the reasonable fitting of the first-order reaction model. As exhibited in Table 4, the Ea values of RH and CDR were 21.85 kJ/mol and 24.00 kJ/mol, respectively. The Ea values of blends with different proportions were lower than those of RH or CDR, demonstrating that the existence of RH could reduce the energy required for CDR pyrolysis. This result also verified that a higher value of activation energy indicates a slower reaction. As the activation energy represents the critical energy required to initiate a reaction [58], the 60RH blend, with the lowest activation energy among these blended samples, is recommended. It is consistent with the positive synergistic effect induced by the addition of RH so that the optimal mixing ratio of RH to CDR is 3:2.

Discussions
This work will be discussed from the following two points: (1) interaction, (2) Pyrolysis products, and (3) kinetics data.
(a) Interaction: The deviation of weight loss TG (∆W) demonstrated that there were synergistic interaction, no interaction and inhibitive interaction between RH and CDR at 76-374 • C, 374-455 • C and 455-1000 • C, respectively. Low temperatures favor synergistic interaction, which is consistent with some previous studies [26,59]. It is reported that the synergistic mechanism was mainly attributed to the catalytic effects of alkali and alkaline earth metals and the transfer of hydrogen and hydroxy radicals [60]. The phenomenon of no interaction generally occurs at the initial stage of pyrolysis due to the low temperature, at which the sample has not started to degrade yet [61]. It is a new discovery that there is no interaction in the middle temperature. This may be the reason for the temporary pause of the blends as the volatiles decrease. The inhibition mechanism was mainly attributed to the carbonization of biomass at high temperatures [62]. Further decomposition of CDR was hindered by a large number of carbonaceous deposits that covered and blocked the molecule pores of CDR residues. (b) Pyrolysis products: All co-pyrolysis products including CH 4 , H 2 O, CO 2 , H 2 and light hydrocarbon were detected via MS. The addition of rice husk reduced the main gaseous products CH 4 and CO 2 . For CH 4 , RH consistently produced more methane than CDR. This result was mainly attributed to the removal of methoxyl substituents of the lignin, cellulose and hemicellulose and the conversion of the alkyl chain of the lignin [63]. CDR was dominant during the pyrolysis of the blends, which reduced methane production. CDR produces a large amount of CO 2 between 400 • C and 600 • C, indicating that a large number of aliphatic groups in CDR were produced by decarboxylation/decarbonylation reaction [64]. (c) Kinetics data: The activation energy of RH in non-catalytic pyrolysis was 21.85 kJ/mol, which is far lower than the results of other studies. Balasundram et al. [65] revealed that the activation energy of RH under non-catalytic action was 49.78 kJ/mol, lower than 53.10 kJ/mol under catalytic action. The kinetic study of CDR has not been reported before this work. López-González et al. [66] reported activation energy of some biomass samples, such as Nannochloropsis gaditana, Scenedesmus almeriensis and Chlorella vulgaris, during pyrolysis in the range of 135-178 kJ/mol. Zhu et al. [26] reported an activation energy value of 71 kJ/mol for bio-oil distillation residue. Sanchez et al. [67] reported that the activation energy of animal manure, sewage sludge and municipal solid waste are 140, 143 and 173 kJ/mol. All samples studied in this paper have low activation energies, mainly due to the synergistic interaction at low temperatures. The synergistic interaction promoted the reaction process and resulted in a significant decrease in activation energy in corresponding conversion stages [62].
Therefore, the future research direction of distillation residues or other industrial hazardous waste can be started from the perspective of interaction and pyrolysis products. In short, compared with the traditional treatment methods, such as landfilling or incinerating, the co-pyrolysis of industrial hazardous waste and biomass would be a better solution.

Conclusions
In summary, an unexpected interaction existed in the co-pyrolysis of CDR and RH. There is a synergistic interaction between RH and CDR from 76 to 374 • C, which disappears in the medium temperature range. The inhibitive interaction occurs from 500 to 1000 • C. All co-pyrolysis products, including CH 4 , H 2 O, CO 2 , H 2 and light hydrocarbon, were detected via MS. Inhibitive interactions reduced the main gaseous product (CH 4 and CO 2 ), and synergistic interactions decreased the activation energy simultaneously. The optimum blending ratio between RH and CDR based on the lowest activation energy of 15.01 kJ/mol is 3:2. The interaction, gas evolution and kinetic parameters will be helpful to large-scale co-pyrolysis of cresol distillation residue and rice husk, and they also provide a promising solution for other distillation residues.