• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Thermogravimetric analysis on the characteristics of oxy-fuel co-combustion of sub-bituminous coal and semi-coke

    2022-09-09 02:23:00LIZhaoyangNIUShengliHANKuihuaLIYingjieWANGYongzhengLUChunmei
    燃料化學學報 2022年8期

    LI Zhao-yang ,NIU Sheng-li ,HAN Kui-hua ,LI Ying-jie ,WANG Yong-zheng ,LU Chun-mei

    (School of Energy and Power Engineering, Shandong University (Shandong 250061,China))

    Abstract: The co-combustion of the low-rank coal with coal derived semi-coke is of great significance to solve the urgent problem of excessively produced semi-coke in China.In this research, the oxy-fuel co-combustion characteristics of Zhundong sub-bituminous coal with bituminous coal derived semi-coke are systematically investigated using thermogravimetric analysis.Compared with air combustion, oxy-fuel atmosphere increased the ignition and burnout temperature by 10 and 40 °C, respectively.Increasing the oxygen concentration to 30% strongly compensated for the slight reduction of the combustion parameters under oxy-fuel condition and much better co-combustion performance was obtained.Three iso-conversional methods, namely, Flynn-Wall-Ozawa (FWO), Kissinger-Akahira-Sunose (KAS) and Starink, were applied to estimate the activation energy, which can be divided into two stages during the co-combustion process.The average activation energy of sub-bituminous coal, the blend and semi-coke were 49.31, 50.82 and 59.00 kg/mol, respectively.Further,the pre-exponential factor and thermodynamic parameters of the enthalpy change, Gibbs free energy change and entropy change were calculated.Interaction indices were innovatively used for both kinetic-thermodynamic parameters and DTG values.An obvious interaction can be observed during the co-combustion process.The kinetic and thermodynamic results demonstrated that the 30% semi-coke ratio was beneficial to co-combustion.Meanwhile, X-ray fluorescence (XRF) and ash fusion analyses proved that the slagging tendency of sub-bituminous coal ash reduced by blending of semi-coke.

    Key words: oxy-fuel co-combustion;coal semi-coke;thermogravimetric analysis;kinetics;thermodynamics;ash sintering

    Cascade utilization of coal is a reliable technology to improve the utilization efficiency and reduce the environmental pollution in recent years[1].This technology extracts high-value oil and gas from lowrank coals through pyrolysis or gasification without direct burning process, and thus realizes the clean and efficient usage.A by-product left by cracking raw coal is called semi-coke, which is commonly used in smelting, calcium carbide and other fields[2,3].However,the supply of the over-produced coal semi-coke has exceeded the demand in recent years and how to deal with the excessive semi-coke becomes an urgent problem[4].

    Semi-coke is a carbon-based material with low volatile and high calorific value.In theory, combustion of semi-coke can provide a lot of thermal energy to realize its reasonable application.However, due to the extensively low content of volatile substances,hydrogen and oxygen, semi-coke is widely recognized with the difficulty in ignition, unstable flame and low combustion efficiency.Thus, it is practically hard to be applied as a single fuel in the pulverized coal furnaces[5].Experiments about semi-coke combustion in circulating fluidized bed boilers (CFBB) have been reported, where the high-temperature preheating combustion approach (>1000 °C), fuel-rich conditions and optimization of the positions of secondary and tertiary air ports can apparently reduce NOxgeneration during the semi-coke combustion[6-8].However, the high cost introduced by the auxiliary heating equipment obviously cannot meet the requirements of the power plants.

    Among the combustion strategies, the cocombustion appears to be proper in improving the poor combustion performance of semi-coke, which attracts many attentions spontaneously.The combustion of blends is believed to achieve a better combustion behavior than pure substances because of synergetic effect between the two fuels.Researchers have carried out lots of studies on co-combustion, such as the binary blending of coals[9-11], coal and biomass blending[12,13],coal and solid waste blending[14].In some of the semi-coke co-combustion studies, synergetic effect occurs during the oil shale semi-coke and biomass cocombustion process[15-17], while other studies observe no significant interaction[18,19].But it should be noted that huge differences exist between the oil shale semi-coke and coal derived semi-coke.For example, the ash content of oil shale semi-coke is extremely high (more than 80%), and its calorific value is almost one order lower than that of the raw coal[20].On the contrary, the calorific value of the coal semi-coke is higher than that of the raw coal[2].By contrasting with the cocombustion studies of oil shale semi-coke and other fuels, kinetic characteristics investigation into the cocombustion of coal derived semi-coke and low-rank coal is rarely employed.From a small-scale plant test,the co-combustion of semi-coke with coal was proved feasible in pulverized coal injection (PCI) operation[21].Further, industrial-scale turbulent chemical simulations showed the semi-coke and bituminous coal can achieve stable co-combustion under certain conditions[22].

    Co-combustion of the high volatile raw coal with the coal semi-coke is a theoretically feasible solution to consume the excessively produced semi-coke, as the addition of semi-coke may strengthen the slagging resistance of the low-rank coals.In previous studies,the ignition and combustion stability of oil shale semicoke could be improved through the high volatile lowrank coal addition[23].Wang et al.[24]used drop tube furnace (DTF) to simulate the co-combustion process in a pulverized coal boiler and concluded that the NOxand unburned carbon ratio generated by the semi-coke combustion were higher.Moreover, the blending ratio significantly determines the NOxformation[4,25].Zhu et al.[26]studied the co-combustion of Shenmu coal with char and found that the ignition performance and combustion stability of the semi-coke were significantly improved.However, the existing experiments did not pay much attention to the byproduct (coal semi-coke) from low-rank coal cascade utilization and may underestimated the great value of coal derived semi-coke by confusing it with oil shale semi-coke.

    Meanwhile, the carbon capture and storage (CCS)is a reliable way to reduce greenhouse gas emission caused by the coal combustion[27].In particular, with the peak of carbon emissions by 2030 and carbon neutral by 2060 plan proposed by China, many coal-fired power plants turn to explore the viable CCS methods[28,29].The oxy-fuel combustion technology can achieve a carbon enrichment higher than 90%.Thus, is considered as a promising commercial approach for CCS[30,31].Moreover, the oxy-fuel combustion was chosen as the only CCS technology introduced in coal power industry, which however has not been put into commercial application in China[32].To achieve effective carbon emission reduction, unconventional combustion technologies like oxy-fuel co-combustion must undergo basic experimental verification to provide theoretical support for further large-scale industrial applications.Meanwhile, the notable slagging problem of low-rank coal may also be relieved due to the blending of coal semi-coke.

    To demonstrate the feasibility of co-combustion of low-rank coal and coal derived semi-coke (not oil shale semi-coke) in thermal power plants, the thermogravimetric analysis was used to study the kinetic and thermodynamic perspectives for cocombustion of low-rank coal and coal derived semi-coke in this paper.Detailed intentions are to (1) investigate the possibility of substituting oxy-fuel combustion for air combustion to achieve convenient carbon capture; (2) evaluate the specific influences of cocombustion conditions on the co-combustion behavior;(3) innovatively assess the possible interaction during co-combustion processes employing synergistic index for both kinetic and thermodynamic parameters;(4) determine the slagging risk of industrial cocombustion application by X-ray fluorescence (XRF)and ash fusion analyses.

    1 Materials and methods

    1.1 Sample preparation

    To obtain the coal semi-coke under real condition,bituminous coal was put into a tube furnace and heated from room temperature to the medium temperature(600 °C) at the heating rate of 5 °C/min and kept for 60 min before the natural cooling process.The whole heating process was conducted under pure N2atmosphere.The raw coal used in this study is Zhundong sub-bituminous coal, which is a well-known low-rank coal in China[33].Both the two samples were ground and sieved to a particle size of less than 75 μm and dried in air at 110 °C for 24 h to ensure the mass kept constant.Table 1 gives the proximate and ultimate analyses of the coal semi-coke and sub-bituminous coal.Compared with sub-bituminous coal, semi-coke contains extremely low content of volatile matter(8.78%) but relatively high content of ash (32.27%).Moreover, the fixed-carbon content is 58.45%, which is much higher than oil-shale semi-coke and reveals the excellent fuel property of coal derived semi-coke.After cracking, the atomic ratio O/C and H/C of coal semicoke diminished to 0.029 and 0.38, respectively,indicating that the bituminous coal rank could be improved and upgraded as another clean fuel as coal semi-coke.The atomic ratio O/C and H/C of Zhundong sub-bituminous coal are 0.22 and 0.61, respectively.To some extent, the blending of semi-coke with subbituminous coal may compensate for the poor combustion performance of the pure semi-coke.

    Table 1 Proximate and ultimate analyses of sub-bituminous coal and the bituminous coal derived semi-coke

    1.2 TGA experiment

    The combustion processes were tested in a thermogravimetric analyzer (NETZSCH STA 2500 Regulus, Germany) with a sensitivity of 10-3mg.Each sample (10±0.5) mg was placed in an Al2O3crucible in the atmosphere composed of O2and CO2to analyze the effect of the oxy-fuel combustion and heated from 30 to 950 °C to observe the whole combustion process.To eliminate the influence of the diffusion effect on combustion, the total flow rate of the carrier gas was controlled at 100 mL/min.In this work, three different heating rates, namely, 20, 25 and 30 °C/min, were chosen to obtain the data for the iso-conversional methods.

    1.3 Combustion characteristics

    Several combustion parameters were analyzed from the TG-DTG curves.The ignition temperature (Ti)is obtained through extrapolation of the TG curve at the DTG peak temperature.The burnout temperature (Tf)was defined as the temperature at which the conversion degree reached 98%.The comprehensive combustion index (CCI) , as described in Eq.(1), was extensively used to evaluate the combustion characteristics, where a higherCCIvalue means a better combustion performance[34].Besides, the combustion stability index(Rw) described as Eq.(2) was applied as well to investigate the combustion stability[19].

    where,DTGmaxandDTGaveare the maximum and average mass loss rates (data fromTitoTf),respectively.Tmrepresents the DTG peak temperature.

    1.4 Kinetic parameters calculation

    The iso-conversional approaches based on the multiple heating rates are considered as the most reliable methods to calculate activation energy for solid fuel combustion.According to International Confederation for Thermal Analysis and Calorimetry(ICTAC) Kinetics Committee recommendations, the iso-conversion methods are also categorized as differentiation method and integration method.The integration method for kinetics is considered more accurate than the differentiation method[35].So, the integration methods are used to calculate the kinetic parameters.

    The kinetic equation for the burning processes is expressed as Eq.(3)[36]:

    wheretis time,f(α)is the mechanism function related with the combustion fraction.The mass conversion,α,can be represented as Eq.(4):

    wherem0,mtandmfrefer to the mass of a sample at initial stage, timetand the end, respectively.

    k(T) is an absolute temperature dependent rate constant that is described by the Arrhenius law in Eq.(5):

    whereEa,k0,Rare the activation energy, preexponential factor and universal gas constant,respectively.

    Taking Eq.(5) and the heating rate ofβ=dT/dtinto Eq.(3), Eq.(6) is then obtained:

    The integration form of Eq.(6) is expressed as Eq.(7)[37]:

    Equation (7) can be solved by the approximation methods or numerical methods because it has no exact solution.At the same time, different iso-conversional methods use different approximation approaches to calculate the kinetic parameters.In these isoconversional methods, Flynn-Wall-Ozawa method(FWO) and Kissinger-Akahira-Sunose method (KAS)are the most classic and most used integration methods for solid fuel combustion.In addition, a more accurate method called Starink based on these two methods was also used to calculate and compare the kinetic parameters.

    1.4.1 Flynn-Wall-Ozawa method (FWO)

    Doyle’s approximation was applied in the integration form of Eq.(7) and then the equation can be written as Eq.(8)[38]:

    At every same conversion degree for the different heating rates, different corresponding temperaturesTcan be obtained.Thus, the regression lines of ln(β)versus 1/Tare fitted and the activation energy is calculated from the slope.

    1.4.2 Kissinger-Akahira-Sunose method (KAS)

    KAS method[39]is described as Eq.(9):

    At every heating rate, a chain of temperaturesTare recorded and the regression lines of ln(β/T2) versus 1/Tcan be obtained at each conversion degree.Then,the activation energy is calculated from the slope.

    1.4.3 Starink method

    Starink method is modified from FWO and KAS and more accurate, as described in Eq.(10)[40]:

    whereCsis a constant.Like KAS method, plot ln(β/T1.92) versus 1/Tat each selectedαand the activation energy is obtained from the slope.

    1.4.4 Pre-exponential factor

    Based on the activation energy calculated, the preexponential factor for the combustion is suggested by the ICTAC guidelines for model free methods through Eq.(11)[34]:

    1.5 Thermodynamic parameters calculation

    The present study also calculated three thermodynamic parameters as a supplement of the combustion behavior, which are the enthalpy change(ΔH), Gibbs free energy change (ΔG) and entropy change (ΔS), as described from Eq.(12) to Eq.(14)[40,41]:

    whereTαrepresents the temperature at the correspondingα, KBis the Boltzmann constant of 1.381×10-23J/K, andhis the Plank constant of 6.626×10-34J/s.

    1.6 Interaction indices

    In the co-combustion process, the combustion behavior differs from both the sub-bituminous coal and semi-coke.In order to illustrate the interaction effect between the coal and semi-coke during the cocombustion process, the theoretical value of the sample blends was calculated from the pure individual sample using the average weight method[42].

    whereYcal,YliandYscrefer to the theoretical value of the blend, pure sub-bituminous coal and pure semi-coke,respectively.XliandXscare the mass percentage of subbituminous coal and semi-coke in the blend,respectively.

    TheMR(Mean Absolute Error) andRMS(Root Mean Square Error) indices[43], as the mostly widely applied evaluation methods, were conducted to evaluate the interaction effect:

    where,andrepresent each experimental value, each theoretical value and the averagetheoretical value.nis the total number of experimental points.A non-zeroMRvalue means interaction effect exists during the co-combustion process andRMSdemonstrates the intensity of the interaction, where a higher value refers to a stronger interaction behavior.

    1.7 Ash composition analysis and ash fusion temperature measurement

    A major factor restricting the utilization of lowrank coal is the risk of ash slagging caused by the high content of alkali and alkaline earth metals.To obtain ash samples, the semi-coke combustion was conducted in a muffle furnace in air at (815 ± 10) °C.As for subbituminous coal, the temperature was set as (500 ±10) °C to avoid the alkali metals transformation during combustion.Then the major elements of the ash samples were analyzed using an X-ray fluorescence(XRF) spectrometer (SHIMADAZU EDX-7200,Japan).Besides, the four characteristic temperatures of the deformation temperatureTD, softening temperatureTS, hemispheric temperatureTHand flowing temperatureTFof the ash samples were tested by an ash fusion point apparatus (ZDHR-3, China).

    2 Results and discussion

    2.1 Combustion characteristics in air and oxy-fuel atmospheres

    The TG-DTG curves of sub-bituminous coal and semi-coke at the heating rate of 20 K/min under 20%O2/80% CO2and 20% O2/80% N2are shown in Figure 1.According to the TG curves, the semi-coke shows almost no mass loss in the temperature range lower than 480 °C, while the sub-bituminous coal loses nearly 5% of its mass, which corresponds to the bound water and small amount light volatiles released from sub-bituminous coal.During the semi-coke preparation process (600 °C cracking), the free water and bound water as well as the most volatiles were released, while sub-bituminous coal only released free water with the air-dried treatment.The main weight loss for subbituminous coal and semi-coke occurs at 400–600 °C and 480–700 °C, respectively.At the accomplishment of the combustion, the mass remained was about 10%for sub-bituminous coal and nearly 40% for semi-coke,which can be ascribed to the large proportion of ash after the removal of volatiles from semi-coke.Compared with the ash content in Table 1, both of the two fuel residuals were higher than ash weight (subbituminous coal ash 3.5%, semi-coke ash 32.27%),indicating that part of the fixed carbon was not fully burned out.As can be seen from Figure 1(b), oxy-fuel atmosphere decreases the maximum mass loss rate by about 15% when compared with air atmosphere.For example, the DTG peak of sub-bituminous coal is about 9.5%/min in air condition while about 8%/min in oxy-fuel condition.The higher specific heat and density of CO2and lower diffusivity of O2in CO2than N2weaken the transport of O2to the fuel particle surface[43].This results in the suppressed volatiles combustion and reactivity and thus, lower the maximum mass loss rate.

    Figure 1 Sub-bituminous coal and semi-coke combustion TG and DTG curves in different atmosphere

    Overall, the two fuels displayed highly similar mass loss characteristics in both oxy-fuel and air atmosphere.But the ignition and burnout temperatures of sub-bituminous coal are 364 and 602 °C in air atmosphere and about 374 and 641 °C in oxy-fuel atmosphere.It is indicated that the overall combustion characteristics change little due to the replacement of N2with CO2.However, the ignition and burnout temperatures increase by about 10 and 40 °C,respectively, due to the higher specific heat of CO2thanN2to postpone the combustion process.As CO2suppresses the maximum mass loss rate and the burnout temperature is delayed under the same condition, the peak values in the oxy-fuel atmosphere are smaller than those in the air atmosphere.Notably, a small peak at about 900 °C in the oxy-fuel atmosphere is observed and it may be caused by char-CO2gasification reaction in the high temperature zone[23].Although the ignition temperature and burnout temperature have increased in oxy-fuel atmosphere, the relatively mild combustion environment is conducive to the smooth operation of the fuels in the thermal plant.

    2.2 Combustion performances of the blends

    The blending ratio is an important parameter to implement the coal co-combustion, where a suitable blending ratio is beneficial to improve the combustion performance and keep stable operation of the thermal plants.Thus, this section focuses on the influence of the blending ratio of semi-coke on co-combustion behavior.

    The TG-DTG curves of the pure sub-bituminous coal, pure semi-coke and semi-coke ratios of 10%, 30%,50%, and 70% in the oxy-fuel atmosphere with 20% O2concentration at the heating rate of 20 K/min are shown in Figure 2.In Figure 2, the curves of the blends fall between those of the pure sub-bituminous coal and pure semi-coke sequentially.With the increased proportion of semi-coke in blends, the ash content and the burnout temperature increase correspondingly.On the other hand, DTG curves show significant differences among the blends.As shown in Figure 2(b),the DTG curve of sub-bituminous coal has a small peak(~100 °C) and a large peak (~440 °C).The former peak corresponds to the loss of the free and bound water,and the latter peak corresponds to the combustion of volatiles and fixed-carbon.In biomass and coal cocombustion process, there are generally more than two large peaks, which are caused by the combustion of volatiles and fixed-carbon, respectively[44-46].In this study, the first half of the sub-bituminous coal combustion peak is steep, and the second half is relatively smooth, indicating that the combustion of the fixed-carbon and volatiles do not separate completely but overlap in partial, which is quite different from the biomass combustion.Further, the combustion of the earlier released light volatile matter may promote the ignition of fixed-carbon and therefore realize more continuous and stabler combustion than biomass.The DTG curve of semi-coke presents a large peak and two small peaks.The large peak is almost ascribed to the fixed-carbon combustion as the volatile content in semi-coke is particularly low.The small peak near 100 °C is caused by the release of the bound water and the peak near the end of combustion is for the char-CO2gasification.Besides, the DTG peaks of the sample changes are minor when the blending ratio of semicoke exceeds 30%, which means this blending ratio begins to limit the volatile combustion in low-rank coal and the combustion rate is dominated by coal semicoke.Yao et al.[23]also found similar phenomenon.This may be due to the coating effect of the porous structure of semi-coke that preventing the contact of low-rank coal volatiles with O2.

    Figure 2 TG and DTG curves of the blends under oxy-fuel atmosphere

    The combustion parameters for each sample of the TG-DTG curves are listed in Table 2, which provide a quantitative analysis to further understand the influence of the blending ratio on the oxy-fuel co-combustion characteristics.As the increase of semi-coke ratio,Ti,TmandTfrange from 363.6–510, 439.8–620.6 and 640.6–750.2 °C, respectively.These characteristic temperatures do not increase monotonically with theincrease ratio of coal semi-coke.This finding is in good agreement with previous study by Zheng et al.[2].In this study, the ignition temperatures of sub-bituminous coal and semi-coke are 363.6 and 510 °C, respectively.The low ignition temperature of sub-bituminous coal is consistent with light volatiles combustion, and it is the main reason for the spontaneous combustion of this kind of low-rank coal.On the contrary, the small amount of heavy volatile residuals and relatively high content of fixed-carbon in coal semi-coke result in the highTiand delay the values ofTi,Tm,Tfas well,therefore cut down the reactivity of semi-coke and its blends.As a result, as the blending ratio of the semicoke increases, the characteristic temperatures of the blended fuels all approach pure semi-coke.

    Table 2 Combustion characteristic parameters of the blends

    According to the trend ofCCI, the value decreases from 5.62×10-7to 2.63×10-7with the blending ratio of semi-coke increase from 0 to 30%, indicating the comprehensive combustion performance of the blended fuel becomes relatively unsatisfactory with a semi-coke blending ratio of 30%.Rwcharacterizes the combustion stability and can be divided into three levels: 4.34×103,2.63×103-3.24×103and 1.63×103-2.16×103.Only the former two levels are relatively suitable for actual combustion.It is inferred that the combustion stability of the pure semi-coke is much lower than that of the pure sub-bituminous coal.With a higher blending ratio of semi-coke, the combustion stability turns worse,where a low combustion stability is not conducive to the smooth operation of the coal-fired boilers.In this study, the suggested proportion of coal semi-coke is 30% to ensure a safe combustion process.

    2.3 Effect of O2 concentration

    As discussed above, the combustion performance will be slightly reduced when air combustion atmosphere switches to oxy-fuel for the different thermo-physical properties between CO2and N2.To enhance the fuel combustion characteristics, increasing O2concentration is a reliable and feasible method.According to the blending ratio results, the 30% semicoke sample is chosen in this section.The sample under the O2concentrations of 10%, 15%, 20%, 25%and 30%, respectively, were adopted at the heating rate of 20 K/min in oxy-fuel atmosphere.

    As shown in Figure 3, the TG curves overlap as the sample heated from room temperature to about 300 °C, where the O2concentration shows little influence on the dehydration and devolatilization process but temperature is the main driving force.In the temperature range of 400–900 °C of the main combustion zone, the TG and DTG curves shifted to the lower temperature sharply with the increased O2concentration and show significant differences in fuel combustion performance.In the main combustion zone,the TG curve gets steep with the increased O2concentration, which also leads to the remarkable changes in DTG curves.The DTG peak at the O2concentration of 10% is only 4%/min and remarkably reaches about 9%/min after increasing the O2concentration to 30%.At the same time, the DTG peak temperature moves from 515 to 465 °C and the peak width reduced by 50% when the O2concentration increases to 30%, which demonstrates the increase of O2concentration causes the intensive combustion of the flammable ingredients and therefore provides higher combustion efficiency of the blended fuel.After the combustion process, the burnout rates are stable at 20% when O2concentration ranges from 10% to 20%,while increasing the O2concentration to 30% will improve the burnout rate to about 15%.The possible reason is that higher O2concentration increase the mass flux rate of O2to the particle surfaces and pores as well as compensates for low flame propagation speed andstability in oxy-fuel atmosphere[43,47].This also provides extra heat feedback to the particle, enhancing devolatilization and combustion[48].

    Figure 3 TG and DTG curves of the 30% semi-coke blend under different O2 concentration

    Several co-combustion characteristics are listed in Table 3 to validate the influences of O2concentration on the co-combustion performance.Notably,TiandTfsignificantly decrease from 454.1 to 330.9 °C and 869.3 to 609.6 °C with the increased O2concentration.The main reason can be attributed to the large porosity and specific surface area of coal semi-coke to be beneficial to O2diffusion.Moreover, the higher O2concentration can oxidize the fuel at lower temperature.At the same time, the maximum mass loss rate is twice that when O2concentration increases from 10% to 30%.As for the other parameters,CCIandRwincrease from 0.74×10-7to 7.21×10-7and 1.56×103to 4.84×103,respectively, to indicate that the influence of O2concentration on comprehensive combustion performance and burning stability is remarkable.Compared with Table 2, an increased O2concentration can achieve even better combustion performance than pure sub-bituminous coal and offset the instability of co-combustion caused by low reactivity of coal semicoke.Also, the increased O2concentration can make up for the relatively negative effects brought by replacing N2with CO2.Based on the above analyses, an increased O2concentration coupled with oxy-fuel combustion is a favorable choice to achieve proper combustion performance and carbon capture at the same time.

    Table 3 Combustion parameters of 30% semi-coke blend in different O2 concentration

    2.4 Kinetic analysis

    Compared with the model-fitting methods, the isoconversional approaches, also named as the model free method, can calculate activation energy without assuming specific models for the reaction and allow the activation energy to be a function of the conversion degree[49].In this part, the pure sub-bituminous coal,pure semi-coke and the 30% semi-coke blend sample were chosen to evaluate the activation energy in the 30%O2concentration oxy-fuel atmosphere, where the temperature heating rates were set as 20, 25 and 30 K/min.And the iso-conversional methods of FWO,KAS and Starink were applied.Further, the activation energies derived from the Starink method were used to calculate the pre-exponential factork0.

    The fitting lines to calculate the activation energy of semi-coke are shown in Figure 4 and theEaresults of the three samples are listed in Table 4.The variations of the slopes confirm thatEais highly correlated with the conversion degreeαthough the contours of the slopes do not show much variation with differentα.Take Starink results of sub-bituminous coal for an example,Eacan divide the combustion into two stages:αfrom 10% to 20% (stage 1) and 20% to 95%(stage 2).At stage 1,Eaincreases dramatically from 35.59 to 115.4 kJ/mol.At stage 2,Eagradually declines from 115.4 to 23.61 kJ/mol.Meanwhile, the regression lines derived from the KAS and Starink methods are similar to each other due to the similarity of the approximate equations used by the two methods.However, the slopes derived from the FWO method look larger than the other two methods.But according to the high value of the correlation coefficient (R), the fitting lines of the three methods all give satisfactory results.

    Figure 4 Regression lines of different iso-conversional methods for pure semi-coke activation energy calculation: (a)FWO, (b) KAS, (c) Starink

    From Table 4, it is observed that the values ofEacalculated by KAS and Starink are close to each other but about 10 kJ/mol lower than FWO method.This phenomenon can be attributed to the high similarity of the two methods, as shown in Eq.(9) and Eq.(10).In fact, both the FWO and KAS methods include some assumptions and approximations.For instance, in FWO method the analytic solution for the temperature integral is not available, and the Doyle’s approximation is used.The iso-conversional methods require the use of different heating rates to calculate the activation energy for a specific mass conversion.While the heating rates at that mass conversion correspond to different temperatures.Hence, it is more practical to integrate temperature and heating rate, like ln(β/T2) in KAS method.Since the temperature termTis in the denominator, which can alleviate the variation of thelnterm caused by the change of the heating rate.Thus, the calculated slope is smaller than that of the FWO method and further the activation energy is also smaller.It is very noteworthy that sub-bituminous coal,which is easier to ignite (see Table 2), has an activation energy larger than coal semi-coke and their blend in the initial stage.At the mass conversion of 20%, the activation energy of sub-bituminous coal, the blend,and semi-coke are 115.40, 94.21 and 76.46 kJ/mol,respectively.The activation energy refers to the energy required for a molecule to transform from the normal state to the activated state, where the chemical reactions are prone to occur.From the results in Table 4, it is inferred that although sub-bituminous coal is easy to ignite (low ignition temperature), the molecular state of sub-bituminous coal is not so active at the beginning of combustion due to its physicalchemical properties.Thus, the energy performance of sub-bituminous coal at stage 1 is not better than that of semi-coke as expected.

    Table 4 Activation energies (Ea) for the samples using different iso-conversional approaches

    Meanwhile,Eaof the blend is between the two pure fuels to indicate the blending of semi-coke in this co-combustion stage improves the energy state of the blend.At stage 2,Eaof both sub-bituminous coal and semi-coke decrease while showing a slight increase within the range ofαfrom 80% to 95% as the 30%semi-coke blend sample.Similar phenomenon was found by Zheng et al.[2].The ash generated from subbituminous coal covers the micropores of the semicoke gradually, which results in the increase ofEaat this stage.Notably,Eaof sub-bituminous coal is lower than that of semi-coke and the value of the blend is also between the two pure fuels.As calculated by Starink,the averageEafor sub-bituminous coal, the blend sample and semi-coke are 49.31, 50.82 and 59 kJ/mol,respectively.Thus, sub-bituminous coal improves the combustion performance of semi-coke in a wide combustion range.Moreover, the pre-exponential factork0increases quickly at stage 1 and decreases sharply at stage 2, indicating that a higher rate of molecular collisions and reaction energy is needed in the beginning surface reaction process[40].The variation trends ofk0are similar withEa, which could be explained by the kinetic compensation effect[38].

    2.5 Thermodynamic analysis

    According to the discussion above, The Starink method is obtained by improving the experimental data on the basis of FWO and KAS methods and its calculation accuracy is higher.Thus, the activation energy derived from Strink method was selected for the calculation of thermodynamic parameters to evaluate the three thermodynamic parameters of enthalpy change ΔH, Gibbs free energy change ΔGand entropy change ΔS, as listed in Table 5.ΔHillustrates the energy difference between the reagent and the activated complex[49].As can be seen in eq.(12), ΔHhas a strong correlation withEa, so it displaysasimilar change tendency withαlikeEadoes.At stage 1 (α< 20%), the ΔHof the blend increases from 53.04 to 88.52 kJ/mol,while decreasing from 88.52 to 35.51 kJ/mol at stage 2(α> 20%).It indicates that more heat energies are required for the initial volatile matters to dissociate the bonds, which agrees with the characteristics ofEa.

    The value of ΔGrepresents the total energy required by the system from the initial state to the activated state[50].As listed in Table 5, the ΔGrange of sub-bituminous coal, the blend and semi-coke is 191.19–199.91 kJ/mol, 202.92–208.97 kJ/mol and 256.64–260.56 kJ/mol, respectively.The ΔGfor cocombustion of sub-bituminous coal with semi-coke is close to pure sub-bituminous coal but much smaller than pure semi-coke.As stated in section 3.2, the 30%ratio of semi-coke cannot significantly reduce the combustion performance of the blending fuel.At the same time, the small increase in ΔGwill not have much negative effect on the stable operation of the boiler.

    The value of ΔSis a state function of the reaction system, representing the disorder degree of the system.Negative value of ΔSstates that disorder degree of the combustion products is lower than reagents and bigger value of ΔSmeans higher reactivity of the sample[51].As shown in Table 5, the values of ΔSare all negative.At stage 1, the reactivity of the three samples all decline but rise up with increasingαat stage 2, which has a same trend asEachanges.Considering the whole process, the blend gains higher values of the apparent entropy change than semi-coke but lower than subbituminous coal to demonstrate the improved reactivity of the pure semi-coke through blending with subbituminous coal.

    Table 5 Thermodynamic parameters of the samples

    2.6 Interaction evaluation

    Two interaction indices are applied to assess the probability of the interaction existence during the cocombustion process of the blend in this study.TheMRindex represents the difference between the actual value and the theoretical value.For the DTG value, a positiveMRindicates the actual mass loss rate of the blend is higher than expected, which means the existence of the positive interaction.But for kineticthermodynamic parameters, a negativeMRillustrates the actual energy required for combustion is smaller than expected, thus, a positive interaction occurs.The judgement methods forMRof DTG and kineticthermodynamic parameters are exactly opposite.Besides, theRMSindex characterizes the interaction degree.

    As shown in Figure 5(a), the DTGMRfactor varies under different semi-coke blending ratios.As the ratio of semi-coke increases,MRdecreases monotonically.When the semi-coke proportion is 30%,MRis positive with a biggestRMSvalue to indicate strong positive interaction occurs during cocombustion.When the semi-coke ratio is more than 30%, theMRbecomes negative with relatively largeRMSvalues, revealing that notable negative interaction occurs and the semi-coke ratio should not increase any more.The possible reason is that the combustion of the sub-bituminous coal volatiles promotes the semi-coke ignition and char combustion, but the high ash content of coal semi-coke can cover up the particle surface and suppress the burning of sub-bituminous coal once the semi-coke ratio exceeds the threshold.In general, the semi-coke blending ratio should be no more than 30%to obtain a good combustion performance under current conditions in this study.

    Figure 5(b) shows the kinetic-thermodynamic interaction behaviors of the blend of 30% semi-coke,which is another vital tool to determine interaction from the energy perspective.TheMRvalues of the four parameters are all negative, indicating that the subbituminous coal and semi-coke affects each other and achieve mutual promotion during the co-combustion.Moreover, the values ofRMSindicate that the blended fuel is easier to activate, which is consistent with the above activation energy analyses.In general, the subbituminous coal and semi-coke in this study show obvious interaction during the co-combustion process and similar discussions can be found in previous studies[52,53].While some other studies show weak interaction during the co-combustion of semi-coke and coals/biomass[19,54], where the variations in different research conclusions can be attributed to the differences in the structure and physical-chemical properties of the fuels used.

    Figure 5 Interaction indices for: (a) DTG of the blends, (b)kinetic-thermodynamic parameters of the 30% semi-coke blend

    2.7 Ash analysis

    The ash compositions and four characteristic ash sintering and fusion temperatures are presented in Figure 6(a) and Figure 6(b), respectively.Besides the small amounts of alkali and alkali earth metals, the semi-coke ash is mainly composed of clay minerals with the content of SiO2and Al2O3more than 70%.Chen et al.[55]found that the SiO2in ash could form stable network and the Al2O3could support the skeleton structure of ash to increase the melting point.On the other hand, the contents of the alkali metals in the sub-bituminous coal ash are 43.3%, 6.3% and 10.6% for CaO, Na2O and MgO, respectively.Li et al.[56]studied the ash fusion characteristics during Zhundong sub-bituminous coal combustion and found that calcium in the compound could form eutectics with aluminum and iron that cause an easier fusion of subbituminous coal ash.

    Figure 6 Ash analysis of semi-coke and sub-bituminous coal:(a) XRF, (b) ash fusion temperature

    Moreover, Figure 6(b) shows four characteristic temperatures of sub-bituminous coal ash which are all lower than that of semi-coke ash.Due to the presence of alkali and alkaline earth metals, sub-bituminous coal ash deforms at about 1130 °C while semi-coke ash deforms only at 1350 °C.During the co-combustion process, SiO2in semi-coke may react with CaO in subbituminous coal to form the calcium compounds and increase the fusion temperature of the blends.The analyses suggest that blending sub-bituminous coal with semi-coke can relieve the risk of slagging and contribute to the cascade utilization of low-rank coals.

    In addition to the influence on slagging, the mineral composition also interferes with the cocombustion process of coal and semi-coke.As shown in Figure 6(a), the combustion of sub-bituminous coal produced abundant calcium oxide and other metal oxides.These metal oxides are natural catalysts for coal combustion.For example, the Na+, K+, Ca2+, Mn2+and Fe3+can evidently reduce the ignition point of coal and promote the reactivity of coal[57].In the oxidation reaction, these metal cations enter the molecular structure of coal and semi-coke by ion exchange and further increase the reactivity of volatile substances.The Ca2+has high fluidity and chemical activity, and the formed CaO has better oxidation catalytic performance on coal combustion as well.Zhao et al.[58]also demonstrated that the alkali and alkaline-earth metal (AAEM) species as important parts of mineral matter could be catalytically active during the coal/char combustion.The Zhundong sub-bituminous coal obtains higher combustion activity due to the high content of AAEM.These minerals not only facilitate the combustion of the low-rank coal itself, but also improve the reactivity of the semi-coke.This phenomenon is consistent with the above discussion in section 2.2.

    3 Conclusions

    The co-combustion performance of Zhundong sub-bituminous coal and bituminous coal derived semicoke is systematically studied by kinetic and thermodynamic analyses.The main conclusions are:

    Compared with air atmosphere, the oxy-fuel atmosphere achieved similar co-combustion behaviors but slightly increase the ignition and burnout temperature by about 10 and 40 °C.An increased O2concentration to 30% enhanced the co-combustion performances greatly.

    The kinetic parameters were compared by different iso-conversional methods.The averageEacalculated by FWO method is about 10 kJ/mol higher than both KAS and Starink method.The averageEaof 30% semi-coke blend sample (50.82 kJ/mol) is close to sub-bituminous coal (49.31 kJ/mol) but much lower than coal semi-coke (59.00 kJ/mol).

    Co-combustion with sub-bituminous coal improved the reactivity of semi-coke.In addition,obvious synergistic effected exists within the blending ratio of semi-coke less than 30%.

    The blending of the semi-coke reduced the risk of slagging.This conclusion can be extended to the cocombustion application of coal semi-coke and other low-rank coals.

    国产精品一区二区性色av| 国产精品一及| 人人妻人人澡人人爽人人夜夜 | 伊人久久精品亚洲午夜| 中文乱码字字幕精品一区二区三区 | 91精品伊人久久大香线蕉| 国产一区亚洲一区在线观看| 久久久久国产网址| 久久久久久久久久久丰满| 亚洲五月天丁香| av又黄又爽大尺度在线免费看 | 爱豆传媒免费全集在线观看| 日本五十路高清| 亚洲最大成人中文| 精品欧美国产一区二区三| 偷拍熟女少妇极品色| 亚洲自拍偷在线| 高清在线视频一区二区三区 | 欧美高清性xxxxhd video| 成人性生交大片免费视频hd| 久久久久久大精品| 26uuu在线亚洲综合色| 国产精品.久久久| 免费人成在线观看视频色| 中文字幕亚洲精品专区| 一个人观看的视频www高清免费观看| 久久久久久大精品| av又黄又爽大尺度在线免费看 | 我要看日韩黄色一级片| 国产不卡一卡二| 亚洲一级一片aⅴ在线观看| 久久久久久久午夜电影| videossex国产| 少妇猛男粗大的猛烈进出视频 | 精品国产一区二区三区久久久樱花 | 国产精品乱码一区二三区的特点| 少妇的逼水好多| av黄色大香蕉| 非洲黑人性xxxx精品又粗又长| 精品免费久久久久久久清纯| 久久精品夜夜夜夜夜久久蜜豆| 久久热精品热| 国产精品一区二区在线观看99 | 狠狠狠狠99中文字幕| 亚洲av成人av| 五月伊人婷婷丁香| 日本黄大片高清| 国产精品99久久久久久久久| 高清在线视频一区二区三区 | 99热精品在线国产| 亚洲国产成人一精品久久久| 一边摸一边抽搐一进一小说| 超碰av人人做人人爽久久| 99热网站在线观看| 国产精品久久久久久精品电影| 永久免费av网站大全| 国产精品国产高清国产av| av专区在线播放| 精品99又大又爽又粗少妇毛片| 午夜激情欧美在线| 午夜福利视频1000在线观看| 久久6这里有精品| 一级毛片久久久久久久久女| 成人国产麻豆网| 国产不卡一卡二| 亚洲成人中文字幕在线播放| 亚洲成人av在线免费| 亚洲av中文字字幕乱码综合| 欧美一级a爱片免费观看看| 精品人妻一区二区三区麻豆| 男插女下体视频免费在线播放| 综合色丁香网| 色哟哟·www| 国产成人91sexporn| 亚洲av福利一区| 国产亚洲5aaaaa淫片| 网址你懂的国产日韩在线| 中文字幕av在线有码专区| 别揉我奶头 嗯啊视频| 国产中年淑女户外野战色| videossex国产| 一个人免费在线观看电影| 精品久久久久久久久av| 精品一区二区免费观看| 天堂中文最新版在线下载 | 国内少妇人妻偷人精品xxx网站| 国产成人精品久久久久久| 大香蕉久久网| 国产精品一区二区在线观看99 | 校园人妻丝袜中文字幕| 少妇的逼好多水| 国国产精品蜜臀av免费| 国产成人福利小说| 精品久久久久久电影网 | 99久久精品国产国产毛片| 日日摸夜夜添夜夜爱| 国产黄片视频在线免费观看| 三级国产精品欧美在线观看| 亚洲无线观看免费| 又黄又爽又刺激的免费视频.| 免费黄色在线免费观看| 国产精品熟女久久久久浪| 亚洲欧美一区二区三区国产| 春色校园在线视频观看| 国产精华一区二区三区| 91久久精品国产一区二区成人| 91久久精品国产一区二区成人| 国产黄色小视频在线观看| 美女xxoo啪啪120秒动态图| 亚洲国产欧洲综合997久久,| www.色视频.com| 久久久久精品久久久久真实原创| 少妇被粗大猛烈的视频| 成人三级黄色视频| 亚洲av中文av极速乱| 国产三级在线视频| 少妇人妻一区二区三区视频| 日本与韩国留学比较| 国产一区二区三区av在线| 人人妻人人澡欧美一区二区| 神马国产精品三级电影在线观看| 国产av在哪里看| 亚洲精品成人久久久久久| 国内揄拍国产精品人妻在线| 波多野结衣高清无吗| 成人毛片60女人毛片免费| 国产老妇女一区| 三级男女做爰猛烈吃奶摸视频| 国产精品乱码一区二三区的特点| 尾随美女入室| 日日摸夜夜添夜夜添av毛片| 亚洲三级黄色毛片| 黄色欧美视频在线观看| 欧美变态另类bdsm刘玥| 中文字幕久久专区| 成人午夜精彩视频在线观看| 久久亚洲国产成人精品v| 国产伦精品一区二区三区四那| 尤物成人国产欧美一区二区三区| 国产色爽女视频免费观看| 久久久久久久久久久丰满| 国产色婷婷99| 日韩在线高清观看一区二区三区| 欧美日本视频| 国产成人freesex在线| 麻豆成人av视频| 伦精品一区二区三区| 网址你懂的国产日韩在线| 九九在线视频观看精品| av.在线天堂| 日本黄色视频三级网站网址| 日本免费在线观看一区| 在线天堂最新版资源| 全区人妻精品视频| 欧美日韩一区二区视频在线观看视频在线 | 级片在线观看| 亚洲图色成人| 国产在线一区二区三区精 | 日韩国内少妇激情av| 特级一级黄色大片| 亚洲国产精品合色在线| 亚洲真实伦在线观看| 男人舔女人下体高潮全视频| 国产伦精品一区二区三区视频9| 久久精品国产自在天天线| 亚洲精品成人久久久久久| 少妇丰满av| 久久这里只有精品中国| 精品久久久久久电影网 | 中文字幕亚洲精品专区| 午夜福利成人在线免费观看| 国产91av在线免费观看| 亚洲久久久久久中文字幕| 精品人妻一区二区三区麻豆| 欧美一级a爱片免费观看看| 中国美白少妇内射xxxbb| 精品不卡国产一区二区三区| 欧美最新免费一区二区三区| 在线a可以看的网站| 最近中文字幕2019免费版| 我的老师免费观看完整版| 男人舔奶头视频| 五月伊人婷婷丁香| 丰满人妻一区二区三区视频av| 久久人人爽人人爽人人片va| 国产精品精品国产色婷婷| 欧美性猛交╳xxx乱大交人| 毛片一级片免费看久久久久| 国产老妇伦熟女老妇高清| 国产69精品久久久久777片| 午夜a级毛片| 久久99精品国语久久久| 亚洲综合色惰| 日本猛色少妇xxxxx猛交久久| 舔av片在线| 日本爱情动作片www.在线观看| 亚洲高清免费不卡视频| 亚洲欧美日韩高清专用| 又爽又黄无遮挡网站| 又爽又黄a免费视频| 波多野结衣巨乳人妻| 久久久久久久久久久免费av| 91午夜精品亚洲一区二区三区| 黄片无遮挡物在线观看| 国产精品av视频在线免费观看| 免费电影在线观看免费观看| 听说在线观看完整版免费高清| 亚洲高清免费不卡视频| 欧美xxxx性猛交bbbb| 欧美极品一区二区三区四区| av黄色大香蕉| 成人午夜精彩视频在线观看| 三级经典国产精品| 91久久精品电影网| 男人舔女人下体高潮全视频| 高清av免费在线| 听说在线观看完整版免费高清| 国产免费视频播放在线视频 | 亚洲精品影视一区二区三区av| 麻豆精品久久久久久蜜桃| 97超碰精品成人国产| 婷婷色综合大香蕉| 六月丁香七月| 别揉我奶头 嗯啊视频| 人人妻人人看人人澡| 岛国毛片在线播放| 日本-黄色视频高清免费观看| 亚洲精品亚洲一区二区| www日本黄色视频网| 国产高清视频在线观看网站| 最近视频中文字幕2019在线8| 国产免费视频播放在线视频 | 精品久久久久久久末码| 亚洲欧美中文字幕日韩二区| 久热久热在线精品观看| 观看免费一级毛片| 久久亚洲国产成人精品v| .国产精品久久| 久久精品人妻少妇| 男女下面进入的视频免费午夜| 精品午夜福利在线看| 舔av片在线| 午夜激情福利司机影院| 中文在线观看免费www的网站| 国产伦在线观看视频一区| 一区二区三区高清视频在线| 欧美bdsm另类| 亚洲av一区综合| 日本一二三区视频观看| 中文乱码字字幕精品一区二区三区 | 久久婷婷人人爽人人干人人爱| 国产一区二区在线av高清观看| 69av精品久久久久久| 亚洲欧美一区二区三区国产| 久久久久久久久久久丰满| 亚洲欧美成人精品一区二区| 免费观看在线日韩| 亚洲欧洲国产日韩| 国产成人精品婷婷| 成人无遮挡网站| 蜜桃亚洲精品一区二区三区| 可以在线观看毛片的网站| 亚洲av二区三区四区| 国产av在哪里看| АⅤ资源中文在线天堂| 性插视频无遮挡在线免费观看| 一级黄片播放器| av在线老鸭窝| 纵有疾风起免费观看全集完整版 | 久久久久久久亚洲中文字幕| 搞女人的毛片| 国产免费福利视频在线观看| 黄片wwwwww| 丰满乱子伦码专区| 插阴视频在线观看视频| 亚洲一区高清亚洲精品| 亚洲在线自拍视频| 亚洲最大成人中文| 久久久亚洲精品成人影院| 亚洲av男天堂| 天美传媒精品一区二区| 水蜜桃什么品种好| 尤物成人国产欧美一区二区三区| 欧美日本视频| 身体一侧抽搐| 国内揄拍国产精品人妻在线| 舔av片在线| 国产老妇伦熟女老妇高清| 欧美另类亚洲清纯唯美| 天天一区二区日本电影三级| 成人毛片60女人毛片免费| 国产精品,欧美在线| 免费不卡的大黄色大毛片视频在线观看 | 水蜜桃什么品种好| 国产在视频线在精品| 网址你懂的国产日韩在线| 国产 一区精品| 18禁在线播放成人免费| 嫩草影院精品99| 亚洲av免费高清在线观看| 亚洲图色成人| www.av在线官网国产| 亚洲成人精品中文字幕电影| 麻豆国产97在线/欧美| 国产亚洲av嫩草精品影院| 国产av一区在线观看免费| 国产探花极品一区二区| 国产精品福利在线免费观看| 九色成人免费人妻av| 日韩高清综合在线| 在线免费十八禁| 在线播放国产精品三级| 搞女人的毛片| 国产亚洲5aaaaa淫片| 99久久中文字幕三级久久日本| 岛国在线免费视频观看| 久久精品国产亚洲网站| 久久精品熟女亚洲av麻豆精品 | 三级国产精品片| 亚洲丝袜综合中文字幕| 午夜福利高清视频| 在线免费观看的www视频| 日本一二三区视频观看| 真实男女啪啪啪动态图| 99视频精品全部免费 在线| 我的女老师完整版在线观看| 一级毛片我不卡| 久久久久久久久久久丰满| 欧美不卡视频在线免费观看| 亚洲自拍偷在线| 91久久精品国产一区二区三区| 男人的好看免费观看在线视频| 精品人妻一区二区三区麻豆| .国产精品久久| 久久婷婷人人爽人人干人人爱| 国产免费视频播放在线视频 | 国产老妇伦熟女老妇高清| 在线播放国产精品三级| 成人亚洲欧美一区二区av| 在线观看一区二区三区| 两个人视频免费观看高清| 99热这里只有精品一区| 国产精品乱码一区二三区的特点| 成人三级黄色视频| 69人妻影院| 精品99又大又爽又粗少妇毛片| 国产免费男女视频| 男人舔奶头视频| 97在线视频观看| 精品人妻熟女av久视频| 国产淫语在线视频| 1000部很黄的大片| 中文精品一卡2卡3卡4更新| 亚洲精品国产av成人精品| 尤物成人国产欧美一区二区三区| 国产黄片美女视频| 搡女人真爽免费视频火全软件| 亚洲aⅴ乱码一区二区在线播放| 晚上一个人看的免费电影| 两性午夜刺激爽爽歪歪视频在线观看| 精品人妻视频免费看| 在线观看美女被高潮喷水网站| 国产视频内射| 国产伦理片在线播放av一区| 只有这里有精品99| av卡一久久| 波多野结衣巨乳人妻| 少妇被粗大猛烈的视频| 成人午夜高清在线视频| 看黄色毛片网站| 久久久久久久久久黄片| 女的被弄到高潮叫床怎么办| 成人毛片a级毛片在线播放| 国产高清视频在线观看网站| 午夜爱爱视频在线播放| 国产成人a区在线观看| 婷婷色麻豆天堂久久 | 亚洲久久久久久中文字幕| 精品不卡国产一区二区三区| av又黄又爽大尺度在线免费看 | 亚洲av电影不卡..在线观看| 黄色日韩在线| 日韩制服骚丝袜av| 国产精华一区二区三区| 成人一区二区视频在线观看| 亚洲精品日韩在线中文字幕| 最后的刺客免费高清国语| 免费一级毛片在线播放高清视频| 国产 一区精品| 在线观看美女被高潮喷水网站| 嘟嘟电影网在线观看| 国内精品一区二区在线观看| 九九久久精品国产亚洲av麻豆| 一区二区三区乱码不卡18| 十八禁国产超污无遮挡网站| 欧美一区二区精品小视频在线| 欧美一区二区国产精品久久精品| 亚洲一区高清亚洲精品| 少妇熟女aⅴ在线视频| 久久精品久久久久久噜噜老黄 | 亚洲欧美成人精品一区二区| 国产一区二区三区av在线| 91久久精品国产一区二区成人| 亚洲美女搞黄在线观看| 人体艺术视频欧美日本| 麻豆国产97在线/欧美| 免费看av在线观看网站| 免费无遮挡裸体视频| 亚洲自偷自拍三级| 一级黄片播放器| 国产精品国产高清国产av| 亚洲av成人av| 建设人人有责人人尽责人人享有的 | 五月玫瑰六月丁香| 一级黄片播放器| www.av在线官网国产| 亚洲av免费在线观看| 少妇熟女欧美另类| 网址你懂的国产日韩在线| 国产高清视频在线观看网站| 国产中年淑女户外野战色| 精品酒店卫生间| 国产免费视频播放在线视频 | 国产黄色视频一区二区在线观看 | 深爱激情五月婷婷| 女的被弄到高潮叫床怎么办| 插逼视频在线观看| 日韩强制内射视频| 午夜福利成人在线免费观看| 亚洲第一区二区三区不卡| 国产视频首页在线观看| 国产91av在线免费观看| 国产精品久久久久久久电影| 亚洲国产欧美人成| 国产精品伦人一区二区| 一级黄色大片毛片| 精品熟女少妇av免费看| 国产黄色视频一区二区在线观看 | 日韩av不卡免费在线播放| 久久99热这里只频精品6学生 | 少妇熟女aⅴ在线视频| 日韩欧美三级三区| 秋霞在线观看毛片| 国产白丝娇喘喷水9色精品| 国内揄拍国产精品人妻在线| 国产片特级美女逼逼视频| 久久久久久久久久成人| 国产精品一及| 少妇熟女欧美另类| 日韩在线高清观看一区二区三区| 亚洲精品乱码久久久久久按摩| 狠狠狠狠99中文字幕| av线在线观看网站| 永久免费av网站大全| 国产精品精品国产色婷婷| 久久精品国产鲁丝片午夜精品| 日韩欧美 国产精品| 国产黄色小视频在线观看| 高清在线视频一区二区三区 | 亚洲av男天堂| 免费黄色在线免费观看| 亚洲经典国产精华液单| 国产成人福利小说| 老司机福利观看| 国产午夜福利久久久久久| 我要看日韩黄色一级片| 亚洲国产精品国产精品| 国产三级中文精品| 好男人视频免费观看在线| 亚洲av中文字字幕乱码综合| 国产不卡一卡二| 简卡轻食公司| 欧美不卡视频在线免费观看| 亚洲内射少妇av| 欧美成人免费av一区二区三区| 亚洲精品亚洲一区二区| 丰满少妇做爰视频| 国产成人午夜福利电影在线观看| 亚洲精品乱久久久久久| 久久久国产成人精品二区| 国产精品一区二区性色av| 成人午夜高清在线视频| 在线播放无遮挡| 国国产精品蜜臀av免费| 丝袜喷水一区| 久久久成人免费电影| 桃色一区二区三区在线观看| 小说图片视频综合网站| 视频中文字幕在线观看| 国产精品电影一区二区三区| 九九久久精品国产亚洲av麻豆| 久久精品夜色国产| 国产高清不卡午夜福利| 免费大片18禁| 可以在线观看毛片的网站| 美女国产视频在线观看| 干丝袜人妻中文字幕| 日本免费一区二区三区高清不卡| 能在线免费看毛片的网站| 搞女人的毛片| 国产视频首页在线观看| 久久人妻av系列| 久久99热这里只有精品18| 国产乱来视频区| 麻豆久久精品国产亚洲av| 国产精品久久视频播放| 欧美+日韩+精品| 久久精品综合一区二区三区| 直男gayav资源| av国产免费在线观看| 一级黄片播放器| 三级国产精品欧美在线观看| 一个人观看的视频www高清免费观看| 日韩在线高清观看一区二区三区| av又黄又爽大尺度在线免费看 | 欧美潮喷喷水| 99热网站在线观看| 男人狂女人下面高潮的视频| 久久热精品热| 欧美激情在线99| 国产亚洲精品av在线| 亚洲av成人av| 级片在线观看| 3wmmmm亚洲av在线观看| 在线免费观看不下载黄p国产| 欧美三级亚洲精品| 国产成人a∨麻豆精品| 中文资源天堂在线| 99在线人妻在线中文字幕| 精品人妻一区二区三区麻豆| 黄色日韩在线| 18禁在线播放成人免费| 免费人成在线观看视频色| 国产精品一及| 亚洲不卡免费看| 欧美日韩在线观看h| 国产成人aa在线观看| 国产极品精品免费视频能看的| 国产精品久久久久久av不卡| 中文字幕av在线有码专区| 国产精品国产三级国产av玫瑰| 色尼玛亚洲综合影院| 热99在线观看视频| 久久人人爽人人片av| 一区二区三区高清视频在线| 亚洲av免费高清在线观看| 国产精品一区二区在线观看99 | 久久精品人妻少妇| 久久久色成人| 最近中文字幕2019免费版| 男女视频在线观看网站免费| 一个人看视频在线观看www免费| 日韩强制内射视频| 国产淫语在线视频| 亚洲av成人av| 国产成人a区在线观看| 美女国产视频在线观看| 日本三级黄在线观看| 国产不卡一卡二| 免费看日本二区| 人人妻人人澡欧美一区二区| 国产精品久久久久久精品电影小说 | 亚洲精品影视一区二区三区av| 国产精品99久久久久久久久| 狂野欧美激情性xxxx在线观看| 秋霞伦理黄片| 国产精品久久久久久精品电影| 午夜精品国产一区二区电影 | 欧美一区二区亚洲| 麻豆一二三区av精品| 少妇的逼好多水| 日韩视频在线欧美| 五月玫瑰六月丁香| 男女那种视频在线观看| 免费电影在线观看免费观看| 精品久久久久久久久久久久久| 午夜a级毛片| 精品久久久久久成人av| 久久精品国产自在天天线| 久久久久久久国产电影| 我的女老师完整版在线观看| 久久人妻av系列| 床上黄色一级片| 看片在线看免费视频| 免费观看人在逋| 夫妻性生交免费视频一级片| 国产精品熟女久久久久浪| 国产极品天堂在线| 精品久久久久久成人av| 亚洲自偷自拍三级| 波多野结衣巨乳人妻| 一区二区三区四区激情视频| 精品国产一区二区三区久久久樱花 | 国产伦在线观看视频一区| 偷拍熟女少妇极品色| 亚洲,欧美,日韩| 日本av手机在线免费观看| 日本熟妇午夜| av视频在线观看入口| 亚洲色图av天堂| 亚洲婷婷狠狠爱综合网| 最近中文字幕2019免费版| 最近手机中文字幕大全| 免费观看在线日韩| 精品欧美国产一区二区三| 永久网站在线| 国产人妻一区二区三区在| 国产爱豆传媒在线观看| 51国产日韩欧美| av在线观看视频网站免费| 免费看av在线观看网站| 少妇被粗大猛烈的视频| 日本三级黄在线观看| 国产私拍福利视频在线观看| 中文在线观看免费www的网站|