PYROLYSIS KINETICS OF HULLESS BARLEY STRAW USING THE DISTRIBUTED ACTIVATION ENERGY MODEL ( DAEM ) BY THE TG / DTA TECHNIQUE AND SEM / XRD CHARACTERIZATIONS FOR HULLESS BARLEY STRAW DERIVED BIOCHAR

The pyrolysis kinetics of hulless barley straw at different heating rates of 5, 10, 15, 20 and 30 oC/min were investigated via thermogravimetry, and the activation energy distribution E and pre-exponential factor k0 were calculated using the Distributed Activation Energy Model (DAEM) from thermogravimetric analysis (TGA) curves, and characterizations of pyrolysis product of biochar were analyzed by techniques of scanning electron microscopy (SEM) and X-ray diffraction (XRD). The pyrolysis process consisted mainly of the dehydration stage (50-150 oC), the active pyrolysis stage (200-400 oC) and the passive pyrolysis stage (400-800 oC). The E ranged from 73.45 to 214.11 kJ/mol within the conversion rate range of 0.10-0.55, and changed from 214.11 to 141.55 kJ/mol within the conversion rate range of 0.55-0.90, and the average value of E was 172.23 kJ/mol. The values of k0 changed greatly with E values at different mass conversion. The wide E and k0 distributions obtained from the kinetic analysis are attributed to the complex chemical reactions of pyrolysis. The structure of biochar was degraded or ruptured due to the increase in temperature. The XRD analysis confirmed that the biochar was amorphous, dominated by disordered graphitic crystallites.


INTRODUCTION
It is well documented that biomass can be converted to various chemicals, biofuels, and solid biochar using modern processes.Biomass is currently increasingly perceived as a renewable resource rather than as an organic solid waste (Liu et al., 2017).Hulless barley is one of the most extensively planted crops in the Qinghai Tibetan Plateau of China.In China's Tibet alone, the yield of Hulless barley reached 680542 t in 2014 (Statistical Bureau of Tibet Autonomous Region, 2015).Subsequently, harvesting and processing of hulless barley result in considerable quantities of crop straw, hulless barley straw, which are usually burned (incineration) in the open, emitting CO 2 , NO x , heavy metals, fine particles, and persistent organic pollutants, or by landfilling, producing leachate with various contaminants such as heavy metals, toxic organic compounds, NH 4 + , and NO 3 − leading to serious air, soil and groundwater pollutions.However, when properly handled, hulless barley straw becomes an abundantly available resource.Therefore, the proper and efficient treatment of crop straw has important significance to maintain environmental quality and for sustainable development (Liu et al., 2017).
Among the methods (microbially mediated transformations and thermochemical route) for the conversion of waste biomass into biofuels or chemical feedstocks, the thermochemical route (pyrolysis, gasification and hydrothermal liquefaction) is one of the prospective approaches.Pyrolysis, as an environmentally friendly and cost-effective technology for biomass conversion, has several advantages over the methods of incineration and landfilling because of its low energy consumption (only approximately 10% of the energy content of biomass is consumed for the pyrolysis itself), and the harmful gas emission in biomass pyrolysis is remarkably weaker than that in incineration.Therefore, in the past few years, pyrolysis has attracted growing interest as a promising versatile platform to convert biomass into valuable resources, and the number of scientific journal papers regarding biomass pyrolysis published per year increases rapidly, and has been thoroughly reviewed (Sharma et al., 2015;Anca-Couce, 2016;Dilks et al., 2016;Hassan et al., 2016;Kan et al., 2016;Patel et al., 2016;Tripathi et al., 2016;Zhang et al., 2016;Chiaramonti et al., 2017;Dhyani & Bhaskar, 2017;Li & Jiang, 2017;Liu et al., 2017;Murillo et al., 2017;Wang et al., 2017;Zeng et al., 2017).
Carbonaceous feedstock pyrolysis is the thermal decomposition of raw materials in the absence of oxygen into gaseous products (mainly CO 2 , CO, H 2 , CH 4 , C 2 H 2 , C 2 H 4 , C 2 H 6 , etc.), liquid products (tar, highmolecular-weight hydrocarbons and water) and solid products (biochar).The pyrolysis process, based on the heating rate and reaction time, can be categorized into three types: slow, fast, and flash pyrolyses (Babu, 2008).The yields and properties of the bio-oil, biochar, and gas are influenced by many parameters, including vapor residence time, temperature, and heating rate.In general, with a low reaction temperature, a slow heating rate, and a long vapor residence time, the pyrolysis produces mainly biochar, while high temperature and heating rate and long vapor residence time increase gas yields.Moderate temperature, high heating rate, and short vapor residence time are beneficial for higher pyrolysis oil yields (Demirbas, 2004;Liu et al., 2017).In slow pyrolysis, lower heating rate (approximately 5−7 K/min) and longer vapor residence times (more than 1 h) provide a suitable ambience and sufficient time for the secondary reactions to complete.Moreover, longer vapor residence time allows the vapors produced during the secondary reactions to be removed.This ultimately results in the formation of solid carbonaceous biochar.
Biochar is a solid carbon-rich by-product (typically more that 60% carbon) of thermal stabilization of biomass or any other organic matter.Biochar is mainly composed of aromatic hydrocarbons, simple carbon substances, or carbon with graphene structure, exhibiting better biological stability and thermostability than the precursor biomass.The microporous structure of the biochar and high carbon content makes it useful for several industrial applications.In agriculture it is used to upgrade the soil quality.Addition of biochar to the soil will increase the rate of carbon sequestration in soil, slowing down the rate of decomposition of nutrients from soil and, hence, finally enhancing the soil quality (Lehmann and Joseph, 2009;Atkinson et al., 2010;Uchimiya et al., 2011;Abdel-Fattah et al., 2015;Hansen et al., 2015;Qambrani et al., 2017).In the purification industry, biochar is used to remove heavy metals such as Cr, Cd, Ni, Hg and Pb.etc., and also for textile dye removal.In power generation biochar can be used as a fuel because biochar contains a high carbon percentage (Yip et al., 2010;Tripathi et al., 2016;Qambrani et al., 2017), and more importantly, the easily tuned surface functionality and porosity make biochar a promising platform to synthesize many carbon-based functional materials (Liu et al., 2015).
Thermogravimetry (TG) is a convenient technique for studying the degradation mechanism and for the determination of the kinetic parameters of pyrolysis and other thermochemical conversion processes (Carrier et al., 2011;White et al., 2011;Anca-Couce, 2016;Dhyani & Bhaskar, 2017;Wang et al., 2017).Kinetic analysis is very important for designing and establishing efficient, safe and reasonable processes.Kinetic parameters obtained are necessary for the prediction of reaction behaviour, and the optimization of the pyrolysis process, so that chemicals or bio-energy can be produced.However, the multiscale complexity of the biomass structures and reactions involved in the pyrolysis process makes the investigation of pyrolysis mechanisms very difficult.
Generally, the chemical mechanism for biomass pyrolysis can be divided into the mechanisms of the three components of cellulose, hemicellulose, and lignin, which mainly depend on the vapor residence time, heating rate, temperature, reactor type, and particle size (Liu et al., 2017).Macroscopic kinetic modeling is often used to simulate the thermodegradation rates of feedstocks or to predict the formation rates of products during biomass pyrolysis, providing valuable information for reactor design and process optimization for biomass pyrolysis.With the development of kinetic theories and computer technologies, several macroscopic kinetic modeling methods with comprehensive mechanism schemes, such as the Distributed Activation Energy Model (DAEM), isoconversional method, lumped kinetic model involving the detailed evolution pathways of products, chemical percolation devolatilization model and kinetic Monte Carlo model, have been developed (Wang et al., 2017).Compared with common modelfitting methods DAEM can simulate mass loss across the overall conversion process with a much higher fitting precision.DAEM assumes that pyrolysis occurs through numerous independent parallel reactions with different activation energies, and the activation energy distribution of these reactions can be described by a continuous function.DAEM does not require a priori assumptions and mathematical model fitting for obtaining the kinetic parameters, only requiring a set of experimental data obtained with three different heating profiles, and has been widely used due to its simplicity and accuracy (Miura,1995;Miura et al.,1998;Cai et al., 2014;Soria-Verdugo et al., 2015).
Although, there are a number of studies on the pyrolysis of agricultural by-products, the significant physicochemical processes associated with the thermal degradation of biomass and their mechanisms are not understood in many important features up to now.To date, there is no reported work regarding the pyrolysis of hulless barley straw using DAME.
The main objective of the present study is to investigate the pyrolysis of hulless barley straw by thermogravimetry (TG) at different heating rates of 5, 10, 15, 20 and 30 ºC/min from ambient temperature to 800 ºC under a nitrogen atmosphere.Based on the TG curves, the activation energy and pre-exponential factor were calculated using the Distributed Activation Energy Model (DAEM) in order to elucidate the mechanism of pyrolysis of hulless barley straw, providing necessary information to fill the present technology gap on the pyrolysis of hulless barley straw.The characterizations of the pyrolysis product biochar were performed by the techniques of scanning electron microscopy (SEM) and X-ray diffraction (XRD).

Materials
The hulless barley straw biomass was obtained from the Napahai Plateau in Southwest China.The hulless barley straw was washed in deionizer water and dried in an oven at 110 ºC, which was then smashed in a crushing machine (Yili QE-300, China) and passed through a 40 mesh sieve.The samples were pyrolysed in a stainless steel box with a lid under a nitrogen atmosphere at different target temperatures (300, 500 and 700 ºC) for 3 h.The treated biochars were designated as QK-300, QK-500, QK-700.
The proximate analysis was performed according to ASTM standards D 1102, E871 and E872 for ash, moisture content and volatile matter, respectively.The fixed carbon is calculated as a difference.The ultimate analysis was performed using a Vario MICRO Cube (Elementar, Germany) elemental analyzer to determine the weight fractions of non-mineral major elements such as carbon (C), hydrogen (H), nitrogen (N), and oxygen (O), using 200 mg oven dried samples in three replicates.For each analysis, the standard samples (i.e., sulfanilic acid and benzoic acid) were first analyzed to check the experimental error within ± 0.1%.The proximate and ultimate analyses of hulless barley straw biomass sample are listed in Table 1.Table 2 lists the ultimate analyses of hulless barley straw derived biochars.

Thermogravimetric analysis
The pyrolysis of hulless straw were performed on a TG 209 F3 Tarsus thermogravimetric analyzer (Netzsch, Germany).A thin layer of the sample (ca. 4 mg) was distributed evenly in a ceramic crucible to limit the thermal gradients in the conditions used.The instrument was previously calibrated for Brazilian Journal of Chemical Engineering assuming infinite number of irreversible first-order parallel reactions with different rate parameters occur simultaneously, where the difference in activation energies is represented by a distribution function f(E), has been extensively used for analyzing the complex reactions occurring during the pyrolysis of different biomass fuels for understanding the reaction kinetics of the pyrolysis process (Sonobe and Worasuwannarak, 2008;Navarro et al., 2009;Shen et al., 2011;Soria-Verdugo et al., 2013;2015;Cai et al., 2014;Chen et al., 2016;Bhavanam et al., 2015;Huang et al., 2017).Miura (1995) and Miura and Maki (1998) proposed two methods to estimate f(E) and k 0 (E), a differential method and an integral method, from a set of three TG experiments at different heating rates, without assuming k 0 values and a functional form for f(E).Although both methods are the same in principle, the integral method is more accurate and does not need a tedious differentiation step.All the reaction activation energies have the same k 0 at the same conversion rate and the activation energy has a continuous distribution.
When the model is used to analyze biomass pyrolysis, the change in the total volatiles, V, with time, t, is given by (Miura and Maki, 1998): (1) where V/V* is the reacted fraction, R the universal constant, T the temperature for each devolatilization, V * is the effective volatile content, V is the volatile content at temperature T, k 0 is the frequency factor corresponding to the E value, and f(E) is the normalized distribution curve of the activation energy that represents the difference in the activation energies of the many first-order irreversible reactions.Usually, the pre-exponential factor is considered to be a constant over the temperature range studied.Then Eq. ( 1) can be simplified to Eq. ( 2) (Miura, 1995;Miura and Maki, 1998).The Ф(E, T) function varies steeply with E at a given temperature T for the heating rate β. (2) (3) where, (4) In this simplified model, the Arrhenius equation can be described as follows:  the temperature under each thermal program.The pyrolysis experiments were performed at heating rates of 5, 10, 15, 20 and 30 ºC/min from ambient temperature to 800 ºC.During heating, mass and temperature of the sample were simultaneously and continuously recorded.In the pyrolysis experiments, N 2 of ultra-high purity was passed through the TGA at a constant flow rate of 60 mL/min to provide the inert gas environment and to sweep away the volatiles released.Every test was replicated at least twice.No significant differences were observed upon comparing the results of replicated tests.

SEM and XRD analysis
To observe the surface morphology pattern of hulless barley straw after pyrolysis, the biochar samples were imaged using scanning electron microscopy (SEM) (Hitachi S-4800, Japan).
XRD analysis was carried out to identify crystallographic structure in the biochar samples using a computer-controlled X-ray diffractometer (X′Pert Pro MPD, Dutch) equipped with a Cu Kα radiation source operated at 65 kV and 50 mA.

Kinetic analysis using DAEM
The so-called distributed activation energy model (DAEM), originally developed by Vand (1943), (5) Equation ( 5) develops a linear relationship between ln(β/T 2 ) and (1/T) with the slope of (-E/R).Miura and Maki (1998) proposed a procedure following the steps: Measure and plot the devolatilization rate, V/V * , as a function of the temperature, T, for the different heating rates, β Plot ln(β/T 2 ) versus 1/T for the different heating rates at each devolatilization rate (Arrhenius plot).
For each V/V * in the Arrhenius plot, linearize the data of the different heating rates and obtain the activation energy E and the pre-exponential factor k 0 from the slope and the intercept, respectively.

Thermogravimetric analysis
The interpretation of derivative thermogravimetric (DTG) curves is often required as a prerequisite to extract preliminary and important thermophysical characteristics of materials.The main DTG peak is usually attributed to cellulose degradation, accompanied by a shoulder at lower temperature, which is related to hemicellulose degradation, and a tail at high temperature corresponding to lignin degradation (Carrier et al., 2011).
Figure 1 shows the thermogravimetric (TG) and derivative thermogravimetric (DTG) curves of hulless barley straw pyrolysis at different heating rates of 5, 10, 15, 20 and 30 ºC/min from ambient temperature to 800ºC.The shape and area of the peaks (Figure 1) follow classical trends.It was found that with increasing heating rate, the TG curve shifted to higher temperature and the peak of the DTG curve became bigger.Below 150 ºC, there was a small weight loss that represents the evaporation of moisture, corresponding to a peak in the DTG curves.The devolatilization of hulless barley straw starts at above 200 ºC and it is completed by 600 ºC.There were three distinct stages of weight loss at the three different heating rates: a dehydration stage (50-150 ºC), a second active pyrolysis stage (200-400 ºC), and a passive pyrolysis stage (400-800 ºC).
During stage I, the small change in loss of the mass is attributed to the loss of water and light volatile compounds in the hulless barley straw samples.Stage II is characterized by a major mass loss, mainly from the decomposition of hemicellulose, cellulose and lignin, which is in the agreement with other biomass pyrolysis (Liu et al., 2017).TG curves (Figure 1a) of the hulless barley straw samples exhibited two mass loss steps.A peaks was observed in the DTG curves (Figure 1b), with the peak temperature at 340 ºC, and the rate of mass loss reached its maximum.In this stage, most of the organic materials were gradually released, resulting in a large loss of mass (more than 70 % of total volatiles) and the formation of the main pyrolytic products.In stage III, the carbonaceous matter decomposed at a slow rate over a broad temperature range (Carrier et al., 2011).
The peak temperature, maximum rate of degradation and the residual weight after 800 ºC for all the biomasses at different heating rates are presented in Table 3.The TG curves showed that they shifted toward the right as the heating rate increased.It generally follows a nonisothermal pyrolysis process.It can be seen from Figure 1b that the maximum rate of decomposition tends to increase at higher heating rate because there is more thermal energy to facilitate better heat transfer between the surroundings and the insides of the samples; the higher heating rates provided less time for reactions to form more volatile products.The temperatures corresponding to the maximum mass loss all increased, which were 314, 329, 331, 333 and Brazilian Journal of Chemical Engineering 336 ºC at heating rates of 5, 10, 15, 20 and 30ºC/min, respectively.

Kinetic analysis
The thermal decomposition reactions of biomass involve a large number of reactions in parallel and in series.The kinetic analysis includes recording of weight loss curves at different heating rates in order to deduce the dependence of the kinetic parameters on conversion level.In this work a DAEM is evaluated from the set of three experimental runs at each different heating rate (5, 10, 15, 20 and 30 ºC/min).The V/V * values with respect to the temperatures (T) are calculated from the TGA experimental data.From the data obtained from the V/V * vs. T relationships at different heating rates, the values of β/T 2 at selected V/V * values (0.1-0.9) are calculated.These values are plotted in the form of ln (β/T 2 ) vs. 1/T plot at the selected values of V/V * .The E and k 0 values are determined from the slope and the intercept of the Arrhenius plots at different V/V * values using the relationship in Eq. ( 5).The activation energies calculated by the DAEM method are listed in Table 3. Figure 2 shows the relationships of V/V * vs. T for hulless barley straw measured at heating rates β of 5, 10, 15, 20 and 30 ºC/min, these being representative plots for the main stage of mass loss.Figure 3 shows the Arrhenius plots of ln β/T 2 vs. 1/T at selected V/V* values for hulless barley straw, showing that all plots had fairly high linear correlation coefficients greater than 0.90 (Table 4), which meant that the linear fitting equation was in accord with the rule of weight losing during the hulless barley straw pyrolysis, indicating that for different V/V*, the volatiles conversion during the thermal decomposition of the hulless barley straw can be described by a set of parallel first-order reactions.The kinetic parameters of the thermal reactions significantly affect their reactivity and sensitivity.The activation energy E is defined as the minimum amount of energy required for the reagents to react; therefore, higher E values mean slower overall reaction rates.The relationship between V/V * vs. E could be obtained by plotting the V/V * value against the corresponding E values, as shown in Figure 4.The activation energy revealed fluctuation at different conversion rates V/V * .Table 4 lists the activation energies and preexponential factors obtained at different conversion rates.The E ranged from 73.45 to 214.11 kJ/mol in the conversion rate range of 0.10-0.55,and changed from 214.11 to 141.55 kJ/mol in the conversion rate range of 0.55-0.90, and the average value of E was 172.23 kJ/mol.The values of k 0 changed greatly with E values at different mass conversions.It is obvious that the activation energy of hulless barley straw pyrolysis is widely distributed, and each conversion rate has a corresponding individual activation energy during the hulless barley straw pyrolysis, and the k 0 value cannot be assumed to be a constant for the pyrolysis of hulless barley straw.Chen et al. (2016) used DAEM to study the kinetics of coffee industrial residue pyrolysis, concluding that the E values changed from 163.42 to 239.49 kJ/mol -1 (average 187.86 kJ/mol -1 ), Sonobe and Worasuwannarak (2008) used DAEM to study several agricultural residue (rice straw, rice husk, corncob and cellulose) pyrolyses and found an E spread from 120 to 250 kJ/mol -1 .That is to say, the energy requirement for hulless barley straw pyrolysis is close to those of other agricultural residue pyrolyses.These results can be attributed to the energy requirements of the complex chemical reactions that occur during thermal conversion.Furthermore, they highlight the effect of conversion on the activation energy.
Figure 5 shows the relationship between k 0 and V/V*, it can be found from Figure 5 and Table 4 that k 0 increased with V/V* within the range of 0.1-0.6, and then decreased with V/V* at the range of 0.6-0.9.
An important measure used to judge the validity of the kinetic model is to examine whether the model with the kinetic parameters obtained can predict the conversion degree of the biomass pyrolysis.Utilizing the f(E) and k 0 (E) obtained by DAEM, the calculated TG curves could be obtained without any mathematical fitting techniques, as shown in Fig. 6.Both the experimental curves and the predicted curves for hulless barley straw at five different heating rates correspond well to each other.Overall, the values calculated for V/V* by the method are in good agreement with the experimental values, demonstrating that the model has the ability to predict the pyrolysis process of the hulless barley straw biomass.

Compensation effect of the kinetic parameters of pyrolytic process
From the results of kinetic parameters, it can be found that the activation energy (E) and frequency factor (k 0 ) obtained under different heating rates were different.This was due to the effects of heating rates on the shapes of TG curves in the thermogravimetric method.So, it was necessary to introduce a"kinetic compensation effect" (Czajka et al., 2016).
where a and b are constants.
Making plots of lnk 0 versus E at different heating rates of hulless barley straw pyrolysis as illustrated in Fig. 7, we could obtain the expression for the kinetic compensation effect of hulless barley straw pyrolysis (Eq.7), ( 7)   The linear coefficient of correlation (R 2 ) is high (0.977), indicating that the kinetic parameters E and k 0 obtained exhibit the kinetic compensation effect under changing heating rates.

SEM and XRD analysis of hulless barley straw derived biochar
Table 2 lists the ultimate analyses of hulless barley straw derived biochars at different pyrolysis temperatures.It should be noted that the carbon content was slightly above that of its precursor in Table 1, and the contents of C tended to increase slightly with increasing temperature.However, H and O highly significantly decreased with increasing pyrolysis temperature.The contents of O and H were reduced from 19.48 and 4.5 wt%, respectively, at 300ºC to 14.02 and 3.30 wt% at 700 ºC.The decrease of H and O contents should result from removal of water, decomposition of lignocellulose in hulless barley straw, and more CO, CH 4 , CO 2 and chemical compounds containing carbon were produced.
Biochar is a solid carbon rich by-product (typically more that 60% carbon) of thermal stabilization of biomass or any other organic matter, which is usually produced at relatively low temperatures (below 700 ºC) (Lehmann and Joseph, 2009).The pyrolysis temperature has a great influence on the biochar property (Uchimiya et al., 2011), in order to understand the evolved porous structure of hulless barley straw derived biochar obtained at different temperatures of 300, 500 and 500 ºC, the SEM and XRD analyses of hulless barley straw derived biochar were carried out.
The SEM images of the biochar obtained through pyrolysis of hulless barley straw are given in Figure 8, illustrating that the pyrolysis temperature had an effect on the structure of the biochars.With the onset of the pyrolysis process, porous structures were found at the surface of the biochar formed.The porous structure of the biochar is created during pyrolysis by the production of volatile matter.As the matter escapes, pores and cracks begin to appear on the surface of the biochar.The porous structure and cracks of the resulting biochars can be clearly seen in the SEM photographs.Just like TG analysis, when the temperature exceeded 200 ºC, pyrolysis entered the second stage and most of the organic materials were gradually released, resulting in a large of mass loss and formation of the main pyrolytic products.After pyrolysis at 300 ºC (Figure 8a), the tissue and cellulosic components of hulless barley straw biomass has not been fully degraded.Large channel-like structures proliferated on the surface of the biochar, but the remains of the vessel structure from the hulless barley straw indicated at least partial retention of the morphology of the tissue structure of the precursor plants.As the temperature increases, more volatile matter will gradually vaporize, causing the pores of biochar particles to be hollow.The surface of biochars was found to have more irregular porous structures, while at higher pyrolysis temperatures changes of surface morphology were apparent (Figures 8b and c).At 500 ºC, the initiation of degradation and cracking starts to occur in the structure of biochar particles due to the lack of remaining fluid components.This shows that the structure of biochar degrades or ruptures due to the increase in temperature.This results in the formation of irregular structure of biochar, thus reducing the uniformity of biochar pore structure.The biochar formed at high temperature is more brittle compared to biochar formed at lower temperature and cannot withstand pressure due to its fragile structure (Claoston et al., 2014).The porous nature of the biochar structure produced at 500 and 700 ºC was more prevalent, but there was more debris, as shown in Figure 6b and 6c.The biochars were found to have microporous and mesoporous structures that could strongly affect their adsorption properties.Brazilian Journal of Chemical Engineering X-ray diffraction is an applicable method to analyze the crystallinity of cellulose in biomass and the structure of bio-char.Figure 9 shows the XRD patterns of biochars obtained at different temperatures, illustrating that the patterns of the biochars show broad graphite peaks are nearly amorphous.A broad peak in the range of 20-30º may indicate the early presence of a poor crystalline structure and a carbon-rich phase in these samples.These types of chars are often known as turbostratic chars, which are dominated by disordered graphitic crystallites (Uchimiya et al., 2011).There is a decrease in the intensity of this broad peak with the pyrolysis temperature, This result confirms the degradation of cellulose or structural composition of biochar sample; the peak at 20º to 30º, as reported by Keiluweit et al. (2010), of the sample with higher temperature will progressively lose intensity and become broader, indicating a gradual decrease in cellulose.It indicates effective destruction of biomass structure at pyrolysis temperatures higher than 300ºC, which is in agreement with the results of TG analysis.The XRD spectra of the biochar showed several peaks, including the peaks at 27º, 28º and 40º.These results indicated that there were mineral crystals such as SiO 2 , CaCO 3 and Na 4 SiO 4 in the biochar.
revealed that the activation energy E ranged from 73.45 to 214.11 kJ/mol in the conversion range of 0.10-0.55,and changed from 214.11 to 141.55 kJ/mol in the conversion range of 0.55-0.90; the average value of E was 172.23 kJ/mol.The values of k 0 changed greatly with E values at different mass conversion.The wide E and k 0 distributions obtained from the kinetic analysis are attributed to the complex chemical reactions of pyrolysis.The distributed activation energy model accurately predicted the conversion degree of the hulless barley straw pyrolysis.The kinetic parameters exhibited kinetic compensation effects and tend to promote pyrolysis as temperature increased.The porous structure of the biochar is created during pyrolysis by the production of volatile matter.The structure of the biochar degraded or ruptured due to the increase in temperature.XRD analysis confirmed that the biochar was amorphous, dominated by disordered graphitic crystallites.

CONCLUSIONS
In this study, pyrolysis of hulless barley straw was carried out in a thermogravimetric analyzer.The kinetics of hulless barley straw pyrolysis were analyzed using DAEM, and the pyrolysis product biochars were analyzed by the techniques of SEM and XRD.TG/ DTG analysis indicated that the pyrolysis process consisted mainly of an active pyrolysis stage (200-400 ºC) and a passive pyrolysis stage (400-800 ºC).DAEM kinetics of hulless barley straw using the distributed activation energy model (DAEM) by the TG/DTA technique and SEM/XRD characterizations for hulless barley straw derived biochar

Figure 4 .Figure 5 .
Figure 4. Activation energy E vs. conversion rate V/V* estimated from the Arrhenius plot

Figure 6 .Figure 7 .
Figure 6.Comparison of experimental and predicted values for conversion degrees of hulless barley straw pyrolysis at different heating rates

Figure 9 .
Figure 9. XRD patterns of biochars obtained at different temperatures.

Table 1 .
Proximate and ultimate analysis of hulless barley straw biomass samples.
a Mean of three hulless barley straw samples with standard deviation
a Biochar produced at 300℃, etc. b Dry basis.c Mean of three biochar samples with standard deviation.

Table 3 .
Thermal degradation of hulless barley straw during the pyrolysis process at different heating rates.

Table 4 .
The activation energies and pre-exponential factors obtained at different conversion rates.