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

    Catalytic Hydrothermal Liquefaction of Algae for Production of Bio-oil with Solid Superacid CatalystSO42?/ZrO2

    2022-10-24 05:44:52LIUMengfanYUXinYUXiaofanZHAOYongnianFENGLijuanLIXianguoandYAOShuo
    Journal of Ocean University of China 2022年5期

    LIU Mengfan, YU Xin, YU Xiaofan, ZHAO Yongnian, FENG Lijuan, LI Xianguo, and YAO Shuo

    Catalytic Hydrothermal Liquefaction of Algae for Production of Bio-oil with Solid Superacid CatalystSO42?/ZrO2

    LIU Mengfan, YU Xin, YU Xiaofan, ZHAO Yongnian, FENG Lijuan, LI Xianguo, and YAO Shuo*

    Key Laboratory of Marine Chemistry Theory and Technology, Ministry of Education, College of Chemistry and Chemical Engineering, Ocean University of China, Qingdao 266100, China

    Solid superacid SO42?/ZrO2as heterogeneous catalyst was prepared to upgrade the bio-oil in the progress of hydrothermal liquefaction (HTL) for the represented algaeofandThe solid superacid catalyst could obviously adjust the composition of the bio-oil and improve the higher heating values (HHVs). The catalytic performance could be regulated by adjusting the acid amount and acid strength of SO42?/ZrO2. Furthermore, it was explored the catalytic effects of SO42?/ZrO2by the HTL for algae major model components, including polysaccharides, proteins, lipids, binary mixture and ternary mixture. The results showed that the introducing of SO42?/ZrO2catalyst could increase the yields of bio-oil from proteins and lipids, and avoid the Maillard reaction between polysaccharides and proteins. Moreover, a possible reaction pathway and mechanisms has proposed for the formation of bio-oils from HTL of algae catalyzed by SO42?/ZrO2based on the systematic research of the producing bio-oil from major model components.

    bio-oil; heterogeneous catalyst; solid superacid catalyst; hydrothermal liquefaction; algae major model components

    1 Introduction

    Bio-oil from renewable biomass sources is the most promising replacement for the exhaustibility fossil fuels. Algae biomass, which own the characteristics of fast growth rate, low requirements on growth environment, inhibition of water pollution and greenhouse effect, are considered as sustainable development and environmental-friendly feedstock for bio-oil production (Biswas, 2017; Wang, 2017; Sharma, 2021).The mainly approaches to produce bio-oil from algae are tran- sesterification conversion, biochemical conversion and thermochemical method (Ambaye, 2021). Compared with above methods, hydrothermal liquefaction (HTL) of algae belonging to the thermochemical method features priority for the following aspects: i) it is an economical process especially for the wet algae feedstock since it avoids the dewatering and drying steps (Saber, 2016; Sheng, 2018; Marzbali, 2021); ii) the components of algae can be fully utilized; iii) the oil products generally possess higher energy content and better stability properties than pyrolysis oil (López Barreiro, 2013; Hu, 2017; Xu, 2020);iv) the decomposition of biochemical compositions into aque-ous-phase, heavy oil and light oil products requires quite low activation energy (Kim, 2013; Vo, 2017a, 2017b; SundarRajan., 2021).

    In order to increase the yield and quality of bio-oil, the solvent water in HTL is heated and pressurized to a subcritical state, which has excellent solubility to biological macromolecules such as proteins and polysaccharides, and can promote the progress of acid-catalyzed reaction (Zhu, 2017; Xu, 2020). It has been found that the ionic product constant (w) of water increases from 10?14to 10?12when the temperature increases from 25 to 350℃ (Toor, 2011). The resulting high levels of H+and OH?can promote the acid or base-catalyzed of HTL process (Cheng, 2017; Galadima and Muraza, 2018). In addition, the introducing of catalysts is an effective strategy to increase bio-oil yield and improve the quality. The application of various homogeneous and heterogeneous catalysts in algal HTL has attracted extensive attention (Gai, 2015;Saber, 2016). Yang(2017) introduced sulfuric acid and acetic acid as the homogeneous catalysts in HTL and the bio-oil yield has been obviously increased. However, heterogeneous catalyst has received extensive attention for the advantages with environment-friendly, convenient recovery and cyclic utilization,.(Carrero, 2015; Navalon, 2016). Solid acid catalyst is a typical heterogeneous catalyst, which is widely used in alkylation, esterification, isomerization, HTL,. In order to achieve prominent catalytic effect, it should focus on the acidity and stability of the solid acid catalyst. Relevant research efforts have confirmed that increasing the acidic of catalyst can improve the fluidity and reduce the boiling point of bio-oil (Yang, 2017). Accordingly, solid superacid which provide strong acidity can be regarded as a recommendable kind of heterogeneous catalyst. SO42?/ZrO2with super acidity have been seen as a representative solid superacid for the advantages of easy preparation, high hydrothermal stability and strong activity. Therefore, the use of solid superacid catalyst SO42?/ZrO2in HTL of algae exhibits broad application prospects.

    The biochemical composition of algae is also an important factor to influence the product yield and quality. Biller(2011) preliminarily reported the related effects using soya oil and microalgae. However, it is still worth further exploring other model compounds such as polysaccharides and proteins. In addition, few studies have an insight into the effects of catalysts on the distribution and properties of HTL product for algae major model components (Peng, 2018). Considering comprehensively,with high polysaccharideandwith high proteins and lipids were selected as representative algae to produce bio-oil by HTL.

    Herein, a series of solid superacid SO42?/ZrO2were prepared by calcining at different temperature to explore the effect of the catalyst acidity on its catalytic performance. Subsequently, it investigated the optimum temperature and dosage of the prepared catalyst on bio-oil production for the selected algae from HTL. Furthermore, it has been conducted the same performance for main components of algae including polysaccharides, proteins and lipids. The analysis of bio-oil product by means of higher heating values (HHVs), thermal gravimetric analysis (TGA) and gas chromatography-mass spectrometry (GC- MS) was conductive to approve the effects of each biochemical composition for bio-oil and the formation pathways.The recycling of the SO42?/ZrO2catalyst used in the HTL confirmed the catalyst still exhibits good catalytic activity after 4 cycles. In this study, it firstly proposed the formation pathways of bio-oil from algae and model com- pounds with the solid superacid SO42?/ZrO2as heterogeneous catalyst. The results could provide a fundamental basis for the screening of catalyst in the HTL of algae.

    2 Msterials and Methods

    2.1 Materials

    Crude polysaccharides, proteins anddry powder were purchased from Shanxi Pioneer Biotech Co., Ltd. Oleic acid and glycerol were mixed in a molar ratio of 1:3 representing a model lipid material since lipids are easily hydrolyzed to fatty acids and glycerol under hydrothermal conditions., which was collected from Stone old man area in Qingdao, China, was used as algae liquefaction raw materials. All chemicals and reagents are analytical grade without treatment. The proximate analysis of moisture, ash, volatile matter, fixed carbon and ultimate analysis of C, H, N of the same batch of algae have been previously reported by our group (Zhao, 2017).

    2.2 Preparation of Catalyst

    SO42?/ZrO2catalyst was prepared by precipitation-im- pregnation. Zirconium oxychloride solution at a concentration of 0.125gmL?1was prepared and the added ammonium hydroxide drop by drop to adjust the pH value to 9 with rapidly stirring. After aging for 12h at room temperature, the precipitate was filtered and washed with deionized water until chlorine-free. The filter cake was dried at 110℃ for 24h and then grind. The resulting powder was immersed in 1 mol L?1dilute sulfuric acid for 12h and filtered. The resulting solid was dried at 110℃ for 2h and then calcined for 4h.

    2.3 Characterization of Catalyst

    The method of back titration is selected to determine the total acid amount of SO42?/ZrO2catalyst (Yang, 2015). 0.1g calcined SO42?/ZrO2catalyst was added to 20 mL 0.1molL?1sodium hydroxide standard solution and reacted for 1h. The liquid was obtained by centrifugation and the remaining sodium hydroxide solution was titrated with hydrochloric acid standard solution calibrated by NaHCO3. The formula for quantity of total acid amount of SO42?/ZrO2catalyst is as follows:

    [H+] (mmolg?1) = (NaOH×NaOH?HCl×HCl)/. (1)

    SO42?/ZrO2catalystwas characterized by D8-Advance X-ray diffractometer (XRD) from Bruker Company, Ger- many. Test conditions: anode target Cu target, diffraction wavelength 1.5406?, operating voltage 40kV, current 40mA, NaI scintillation counter, step length 0.05?step?1, scanning speed 3?min?1, scan angle 2=10?–70?. The sur- face morphology and elements of SO42?/ZrO2were observed by scanning electron microscopy (SEM, Zeiss Su- pra 55) and energy X-ray dispersive spectrometer (EDS, X-MAX-20), respectively.

    2.4 Apparatus and Experimental Procedure

    The HTL process was performed in a stainless steel autoclave (CJK-(0.2)), equipped with an electrically heated furnace, a magnetic stirrer and a temperature controller. In a typical run, the autoclave was charged with suitable polysaccharides, protein powder, and oleic acid- glycerol mixtures or algaeThe ratio of material to liquid is 0.1gmL?1. The autoclave was sealed and purged with nitrogen for 5min before heated from 200 to 300℃ with 20℃ intervals and held for 10min at each temperature. Then autoclave was cooled down to room temperature. The gas products produced by HTL were released from the reactor and then the reactor was opened. Liquid products were separated by a series of filtration and extraction procedures.

    Dichloromethane (CH2Cl2) was added after filtering out the residue to extract crude oil phase. The CH2Cl2phase was evaporated at 35℃ under reduced pressure to remove the solvent, and the remaining liquid product is defined as bio-oil. The aqueous phase was vacuum evaporated at 60℃ to remove water and the residue was defined as water-soluble organics (WSOs). The residue was dried at 110℃ for 12h and then weighed. Product yield was calculated by the following equations:

    y=m/m×100%, (2)

    y=m/m×100%, (3)

    y=m/m×100%, (4)

    y=100%?y?y?y, (5)

    wherey,y,yandycorrespond to the yield of bio-oil, residue, water-soluble organics and gas, respectively.m,m,mandmrefer to the mass of bio-oil, residue, water-soluble organics and feed, respectively.

    2.5 Analysis of Bio-oil

    Higher heating values (HHVs) of bio-oils were analyzed by using the ZR-3R combustion heat experimental device produced by Nanjing Dosuke Technology Development Co., Ltd. with 0.001℃ temperature difference re- solution, and the temperature resolution is 0.01℃.

    Gas chromatography-mass spectrometry (GC-MS) ana- lyses of bio-oils were carried out using an Agilent 7890A/5975C (Santa Clara, CA, US) with a DB-5 column (30m×250μm×0.25μm). The carrier gas was helium with a flow rate of 1.5mLmin?1. 1mL of CH3OH solution of the sample (0.1–0.2gmL?1) was injected into the column. The injector was set at split mode with an inlet temperature of 280℃ and the diversion flow was 20mLmin?1. The GC-MS oven temperature was held at 40℃ for 3min and then raised to 300℃ with a heating rate of 10℃min?1. Compounds were identified by means of the National Institute of Standards and Technology (NIST) library of mass spectra.

    Thermal gravimetric analysis (TGA) of bio-oils was performed in nitrogen atmosphere. Samples were purged to a constant weight at room temperature and then heated to 900℃ with a heating rate of 15℃min?1. The gas flow rate was 30mLmin?1.

    2.6 Reusability of Catalyst

    The catalyst was used for successive batches. During each cycle, the residue obtained by HTL of.were collected and white powder was obtained after calcined at 550℃ for 4h. Then the white powder was added into the HTL of.for the next cycle. The yields and HHVs of bio-oil were measured according to the assay described in the previous paragraph.

    3 Results and Discussion

    3.1 Characterization of Solid Superacid Catalyst SO42?/ZrO2

    Acid amount is an important property of acid catalyst, while the influence of total acid amount on the catalytic performance of SO42?/ZrO2catalyst is less studied. The calcination temperature has a particularly prominent effect on its acid amount. Therefore, it was measured the total acid amount of SO42?/ZrO2catalyst calcined at different temperature (Fig.1a). The total acid amount of SO42?/ZrO2catalyst fluctuates slightly but maintains above 2.5mmolg?1when the calcination temperature ranges from 200 to 400℃. With the further increase of calcination temperature, the total acid amount of SO42?/ ZrO2catalyst shows a decreasing trend until it is close to zero at 800℃. The decrease of acid amount can be ascribed to the reduction of total hydroxyl content on the surface, which is influenced by the surface hydroxyl density and surface area of solid acid catalyst with the calcination temperature increases (Qin, 2004).

    Fig.1 Total acid amount (a) and XRD patterns (b) of SO42?/ZrO2 catalyst calcined at different temperatures.

    Acid strength is another fundamental property of solid acids. The acid strength of SO42?/ZrO2catalyst calcined at 400, 550 and 600℃ had been studied by our group (Zhao, 2017). The proportion of strong acidity of SO42?/ ZrO2catalyst calcined at 550℃ was more than 60% which is higher than 400 and 600℃. Herein, we further explore the essential reasons by XRD, which indicates the difference of acid strength can be ascribed to the different crystalline state of SO42?/ZrO2catalyst at different calcination temperature. Since sulfate trends to combine with materials with tetragonal crystal system, the SO42?/ZrO2catalyst with tetragonal crystal system possesses strong acidity and shows better catalytic performance than other types. As the XRD patterns in Fig.1b, the broad peak indicates the solid superacid catalyst SO42?/ZrO2calcined at 400℃ was amorphous. Therefore, the lower calcination temperature is not conducive to the formation of surface acid sites for SO42?/ZrO2catalyst. The peaks of tetragonal crystals appeared with the calcination temperature rising to 550℃, indicating that more tetragonal crystals appeared in the SO42?/ZrO2catalyst and the surface acidity enhanced. In SO42?/ZrO2solid superacid, sulfate acts in a double coordination form with metal ions on the surface of metal oxide. The formation of SO42?/ZrO2solid superacid center is mainly due to the coordination adsorption of SO42?on the surface of metal oxide. However, the intensity of the peaks of tetragonal crystals decreased when the calcination temperature further increased to 600℃, which indicates that an excessively high calcination temperature would have an adverse effect on the formation of tetragonal crystals and reduced the acidity of SO42?/ZrO2catalyst. Overall, it can further conclude that 550℃ is the optimum calcination temperature for the preparation of SO42?/ZrO2catalyst. Fig.2 shows the surface morphology and elements of SO42?/ZrO2catalyst, which has a smooth surface, indicating that the material supported on the surface has not been shed by high temperature (Wang, 2009). Sulfur was detected in the EDS spectrum, suggesting that the sulfate ions had been successfully supported on ZrO2.

    Fig.2 SEM images and EDS spectrum of SO42?/ZrO2 catalyst.

    3.2 Effect of Temperature and SO42?/ZrO2 Catalyst Dosage on Bio-oil Yield from HTL of Algae

    The HTL of algae is an extremely complex process that involves interaction among polysaccharides, proteins and lipids. In order to understand the effect of different biochemical composition of algae on the yield of bio-oil obtained from HTL,with high polysaccharide andwithhigh proteins and lipids are selected as representative algae.

    In order to explore the optimum temperature for the HTL ofand, the experiment without catalysts were conducted firstly at different liquefaction temperature and the product distribution results were shown in Fig.3a. The optimum liquefaction temperatures forandare 280 and 300℃ and the corresponding maximum yield of bio-oil are 12.63% and 23.39%, respectively (Fig.3b). The yields of residue fromHTL are higher than those ofat different temperatures.The optimum temperature and maximum bio-oil yield forare both lower than those for, which indicated that algae with higher proteins and lipids content required higher liquefaction temperature and higher bio-oil yield and lower residue yield can be obtained.

    Fig.3 Product yields from HTL of Enteromorpha prolifera (E) and Chlorella vulgaris (C) versus liquefaction temperature (a) and the dosage of SO42?/ZrO2 (b).

    Keeping the optimum reaction temperature and adding SO42?/ZrO2catalyst, the bio-oil yield produced byfluctuates slightly compared to which was catalytic with ZrO2carrier or no catalyst, and maintains at about 10% and the yield of residue is almost constant. Nevertheless, the yield of bio-oil forincreases firstly and then decreases with the increase of SO42?/ZrO2catalyst dosage compared to which was catalytic with ZrO2carrier or no catalyst. It indicates that excessive use of SO42?/ZrO2catalyst would inhibit the formation of bio-oil. The residue yields ofare higher than those ofat different SO42?/ZrO2dosage. As a result, SO42?/ ZrO2catalyst is more effective for the HTL of algae with higher proteins and lipids content, which can significantly increase the bio-oil yield and decrease the residue yield.

    3.3 Effect of Temperature and SO42?/ZrO2 Catalyst Dosage on Product Yields from HTL of Model Compounds of Algal Main Components

    To further explore the influence of main composition of algae including polysaccharides, proteins and lipids on the yield and composition of bio-oil and understand the bio-oil formation pathway, the similar experimental were performed to study the effect of temperature and addition of SO42?/ZrO2catalyst on product distribution obtained from HTL of model compounds of algal main components.

    HTL experiments were then conducted for three major model components at different temperatures (Fig.4a). The yield of bio-oil from lipids is significant higher than that of polysaccharides and proteins, since polysaccharides are freely soluble in water at room temperature and easily de- composed into water-soluble polar small molecules under the subcritical hydrothermal environment. Hence, polysaccharides had little contribution to the formation of bio- oil. In virtue of the corresponding optimal yield, 200, 300 and 280℃ were respectively selected as the optimal liquefaction temperature for polysaccharides, proteins and lipids.

    In addition, it is also investigated the effect of SO42?/ ZrO2catalyst dosage on the HTL for three major model components at the corresponding optimal liquefaction temperature (Fig.4b). It can conclude that WSOs and gas products are the main products for polysaccharides in HTL. The bio-oil yield does not affect by the amount of SO42?/ZrO2catalyst, while the residue yield increases with SO42?/ZrO2catalyst dosage. This is due to the hydrolysis of polysaccharides in a subcritical water environment to produce 5-Hydroxymethylfurfural (5-HMF) which is prone to polymerize in an acidic environment to generate a coke-like substance. For proteins and lipids, only bio-oil yields are the subject of discussion. The bio-oil yields of proteins and lipids reach the highest when the amount of SO42?/ZrO2catalyst is 5% and 9%, respectively. By comparing the performance of three major model components in HTL, it is further verified that the contribution of main components of algae to bio-oil is lipids>proteins>polysaccharides.

    Fig.4 Product yields from HTL of polysaccharides (Ps), proteins (Pr) and lipids (L) versus liquefaction temperature (a) and the dosage of SO42?/ZrO2 (b).

    3.4 Effect of Temperature and SO42?/ZrO2 Catalyst Dosage on Bio-oil Yield from HTL of Binary Mixtures

    The effects of reaction between polysaccharides and proteins on the HTL of algae can be roughly judged by the bio-oil yield distribution and bio-oil composition. A series of different proportions of polysaccharide and protein mixtures are hydrothermally liquefied at 200, 250 and 300℃, respectively (Fig.5a). When liquefaction tem- perature is 200℃, the increase of polysaccharides content contributes to the formation of bio-oil. However, the yield of bio-oil shows a decreasing trend in varying degrees with the increase of polysaccharides content at liquefaction temperature of 250 and 300℃. This is consistent with the higher bio-oil yield of proteins HTL at high temperatures and the higher bio-oil yield of polysaccharides HTL at low temperatures.

    The influence of SO42?/ZrO2catalyst dosage with different proportion of polysaccharides-proteins binary mix- tures on bio-oil yield is evaluated by the experimental bio-oil yield (solid line) and theoretical bio-oil yield (dash line) at the optimum liquefaction temperature of 300℃. As seen in Fig.5b, the bio-oil yield slightly increases with the addition of SO42?/ZrO2catalyst for the binary mixtures with 20% and 40% polysaccharide, while the yield of bio-oil decreases slightly with polysaccharides content of 80%. In addition, the experimental bio-oil yields with polysaccharides content of 80% are closer to the theoretical bio-oil yields comparing with the polysaccharide content of 20% and 40%.

    HTL experiments are conducted for equal mass proportion of binary mixtures for polysaccharides-lipids and proteins-lipids, and the yields of bio-oil varying with liquefaction temperature are shown in Fig.5c. The yield of bio-oil from polysaccharides-lipids decreases with liquefaction temperature increasing which is opposite of the result of lipids HTL alone. It is speculated that the increase of liquefaction temperature is not conducive to polysaccharide HTL and even cause side reactions of polysaccharides and lipids, which led to the decrease of bio-oil yield. The similar results have also been found by Yang’s group (Yang., 2015, 2019). On the contrary, the yield of bio-oil from proteins-lipids increases with the increase of liquefaction temperature and reaches the highest at 300℃.

    Fig.5 Bio-oil yields from HTL of polysaccharides-proteins versus polysaccharides content (a) and the dosage of SO42?/ZrO2 (b). Bio-oil yields from HTL of polysaccharides-lipids and proteins-lipids binary mixtures versus liquefaction temperature (c) and the dosage of SO42?/ZrO2 (d). Solid line, test bio-oil yield; dash line, theoretical bio-oil yield.

    HTL experiments with varying amount of SO42?/ZrO2catalyst are conducted on equal mass proportion of binary mixtures for polysaccharides-lipids and proteins-lipids (Fig.5d). The yield of bio-oil from proteins-lipids HTL reached the highest with 5% SO42?/ZrO2catalyst addition. On the contrary, the yield of bio-oil decreases for polysaccharides-lipids HTL with the addition of SO42?/ZrO2catalyst.

    3.5 Effect of Temperature and SO42?/ZrO2 Catalyst Dosage on Bio-oil Yield from HTL of Ternary Mixtures Simulated Algae

    In order to make use of HTL results of model compounds more accurately to infer the actual liquefaction process of algae, three major components were mixed to simulate algae. Simulatingandwere performed by mixing polysaccharides, proteins and lipids at a mass fraction of 80%, 18% and 2%, 10%, 60% and 30%, respectively. HTL experiments are conducted at different liquefaction temperatures and different addition of SO42?/ZrO2catalyst. The optimal bio-oil yield of simulatedis at 220℃ shown in Fig.6a. The trend of bio-oil yield varying with liquefaction temperature of simulatedis different with the real one for the complexity of algae structure. With respect to simulating, the bio-oil yield is significantly higher than that of simulatingat the optimal liquefaction temperature of 300℃, which is similar to the HTL results of real algae. This indicates that higher liquefaction temperature can promote HTL of algae with higher protein and lipids content and increase the yield of bio-oil.

    The results of bio-oil distribution from simulated and real algae varying with the dosage of SO42?/ZrO2catalyst are shown in Fig.6b. For the simulation of, the yield of bio-oil hardly changes with the amount of SO42?/ZrO2catalyst. The possible reason is that polysaccharides occupy a larger proportion in the ternary mixtures and the products distribution results are similar to that of polysaccharides HTL. Besides, as for the simulated, the yield of bio-oil first increases and then decreases with the increase of SO42?/ZrO2catalyst dosage.

    Fig.6 Bio-oil yields from HTL of ternary mixtures and real algae versus liquefaction temperature (a) and the dosage of SO42?/ZrO2 (b).

    3.6 Recycling of SO42?/ZrO2 Catalyst

    The reusability of the SO42?/ZrO2catalyst was studied by successive batches. As shown in Fig.7, the SO42?/ZrO2catalyst still exhibits 92.06% of its initial conversion after 4 cycles. From the XRD patterns of SO42?/ZrO2catalyst after each circle (Fig.8), the tetragonal peak of SO42?/ ZrO2catalyst is substantially the same as original SO42?/ ZrO2catalyst after calcination, indicating the good cyclic regeneration performance. Therefore, the solid superacid SO42?/ZrO2with high hydrothermal stability is considered to be a kind of good heterogeneous catalyst for regeneration recycling.

    Fig.7 Bio-oil yields and HHVs from HTL of Chlorella vul- garis varying with cycling times of SO42?/ZrO2 catalyst.

    Fig.8 XRD pattern versus reuse times of SO42?/ZrO2 catalyst: without reuse (a); reused once (b); reused twice (c); reused three time (d).

    3.7 Elemental Analysis of Bio-oils

    Elemental analysis results with HHVs of bio-oils produced from three major model compounds are shown in Table 1. It shows that hydrogen content for catalytic conditions of bio-oil increases significantly and carbon content decreases obviously, which resulted in an increase in the H/C ratio and HHVs of the bio-oil obtained in the presence of SO42?/ZrO2catalyst. High H/C ratios cause low molecular weights, which associates with TGA analysis, indicating that the light component proportion increases greatly when SO42?/ZrO2catalyst are employed. By using GC-MS, it can further conclude that the addition of SO42?/ZrO2catalyst actually promotes the decarboxylation of amino acids and causes higher H/C ratios, which is an efficient way to increase HHVs.

    Table 1 Elemental analysis of bio-oil from HTL of algal major components

    Note:?calculated by difference.

    3.8 TGA of Bio-oils

    Although thermal decomposition may occur during the heating process, TGA provides an estimate of the boiling point range of bio-oil (Yuan, 2010;Nazari, 2015; Yang, 2017). In this study, it carried out TGA of bio-oils obtained from catalytic and non-catalytic HTL of three major model compounds (Fig.9). The estimated boiling point distributions of bio-oils are shown in Table 2, which can be roughly divided into three temperature intervals. The results show that the addition of SO42?/ ZrO2catalystslightly increases the content of light component for bio-oil obtained from polysaccharides and lipids. It also proves that SO42?/ZrO2catalyst possesses strongly acidity and the acid condition is beneficial to the cracking reaction to get smaller molecular weight organics. Among three model compounds, the highest content of light components in bio-oil is obtained from lipids, which reaches 73.68% with SO42?/ZrO2catalyst.

    For bio-oil from HTL of polysaccharides, the medium component proportions decrease greatly and mainly convert into light components in the presence of SO42?/ZrO2catalyst. It indicates that the SO42?/ZrO2catalyst promotes algae to split into more small organic molecules and significantly narrow the boiling point range of bio-oil of polysaccharides.

    Fig.9 TGA curves of bio-oils obtained from HTL of polysaccharides (a), proteins (b) and lipids (c).

    Table 2 Boiling point distributions of bio-oil estimated by TGA

    3.9 GC-MS Analysis of Bio-oils and General Reaction Network

    Compared with bio-oil, the composition of gas and aqueous phase products from HTL of algal major components is relatively simple. The gas phase products are mainly composed of CO2, H2and hydrocarbons. The largest proportion of gas is CO2from HTL of polysaccharides and proteins analyzed by GC. The main components in the aqueous phase formed by HTL of polysaccharides are acids and alcohols, while that from HTL of proteins is Nitrogen-containing compounds such as pyrazine, pipe- razine and pyrrole (Yang, 2017). Herein, we emphatically analyze the composition of bio-oil from HTL of algal major components.

    During the HTL, the biochemical composition of bio- oil is mainly influenced by catalyst. Therefore, the chemi- cal composition of bio-oil obtained with and without SO42?/ZrO2catalyst was analyzed at the same time by GC-MS. The total ion chromatogram (TIC) of bio-oil obtained is shown in Fig.10. Although SO42?/ZrO2catalysthas little effect on the yield of bio-oil, it significantly changed the composition of bio-oil from HTL of polysaccharides. Organic products of polysaccharides liquefaction are mainly divided into five categories: cyclic ketones, furfural, furans, hydrocarbons and phenols (Table 3). Among these organic products, 5-HMF has the largest proportion without the catalyst, while the content of 5-HMF decreased and the content of furan increased significantly in the presence of SO42?/ZrO2catalyst. It is speculated that the addition of SO42?/ZrO2catalyst promotes the degradation of 5-HMF into furans in the HTL process. To some extent, furans are inclined to polymerize to undesired coke-like substances, which led to an increase in the yield of residue with the addition of SO42?/ZrO2catalystimproves. Parsa(2017) reported that furfurals or furans are more unstable under alkali condition or neutral conditions. They tended to decompose by dehydration and partial degradation of furfural or furans to produce phenolic compounds. In this study, due to the strongly acidity of SO42?/ZrO2catalyst, the hydrolysis of furfural or furans to phenols is inhibited, resulting in a decrease in phenolic substances. Cyclic ketones can be formed by dehydration, isomerization, and cyclization of sugars and reducing sugars. The formation of cyclic ketones is inhibitd by SO42?/ZrO2catalyst. The way of bio-oil deoxidation during HTL is generally achi- eved by removing CO2or water molecules. The removal of oxygen in the form of water is inhibited, but the formation of CO2is promoted, and the H/C in bio-oil is improved, so the HHVs of bio-oil are promoted in the presence of SO42?/ZrO2catalyst.

    Fig.10 TIC spectraof bio-oil obtained from HTL of polysaccharides (a), proteins (b), lipids (c) and polysaccharides-pro- teins (d).

    Table 3 Major chemical compositions of bio-oils obtained from HTL of polysaccharides

    Major bio-oils components from HTL of proteins can be divided into ketones, phenols, amines, amides, alcohol and various N-containing heterocyclic compounds (Table 4). The ketones in bio-oil are formed by the process of intermolecular condensation, hydrolysis and isomerization of amino acids produced by proteins. The results show that N-containing heterocyclic compounds posses- ses a larger proportion without catalyst and the total content decreases with the addition of catalyst, which proves that SO42?/ZrO2catalyst promotes the denitrification reaction. In addition, the reduction of amine compounds indicates that SO42?/ZrO2catalyst may inhibit the decarboxylation reaction of amino acid, resulting in the decrease of released CO2, H/C and HHVs of bio-oil. The substantial increase of alcohols indicates the addition of SO42?/ZrO2catalyst promoted the deamination reaction of amino acids, while the main product carboxylic acids further reduced to alcohols. The hydrocarbon substance increases significantly, which is possibly due to SO42?/ ZrO2catalyst promote the decarboxylation of amino acids to generate alkanes and alkenes.

    Lipids are the mixtures of glycerol and oleic acid in which glycerol is generally dissolved in the aqueous phase without being converted to oil phase products. The oleic acid is mainly reconstituted by HTL to form four kinds of fatty acids (Table 5). The addition of SO42?/ZrO2catalyst promoted the formation of linoleic acid, palmitic acid and-Hexadecenoic acid but inhibited the formation of 6-Octadecenoic acid. In addition, a small amount of hydrocarbons was detected, which mainly derived from polysaccharides and proteins. The total content of fatty acid decreased, which is possibly due to the further decarboxylation of fatty acid to form hydrocarbon in the presence of SO42?/ZrO2catalyst.

    Table 4 Major chemical compositions of bio-oils obtained from HTL of proteins

    Table 5 Major chemical compositions of bio-oils obtained from HTL of lipids

    The resulting bio-oil from HTL of polysaccharides- proteins binary mixtures is mainly composed of N-con- taining heterocyclic compounds, cyclic ketones, amides, phenols, alcohols and other substances (Table 6). Among them, N-containing heterocyclic compounds such as pyra- zines and pyrroles have a large proportion. Yang(2017) concluded that these N-containing heterocyclic compounds were formed by Maillard reactions between polysaccharides and proteins degradation products. Gai(2015) reported that a part of N-containing heterocyclic compounds were formed by the reaction of protein- produced ammonia and polysaccharide-produced cyclic oxygen compounds. By introducing SO42?/ZrO2catalyst, the pyrazines content in the bio-oil decreases which indicates that the catalyst might inhibit the above two reactions to increase the bio-oil yield. In addition, the content of cyclic ketones in the bio-oil also decreases with the addition of SO42?/ZrO2catalyst, since the cyclic ketones are the product of Maillard reaction in which are inhibited by the catalyst. However, the contents of phenols and al- cohols increases obviously, indicating the SO42?/ZrO2catalyst promoted the deamination of amino acids and the generation of alcohols reduced from carboxylic acids. The resulting alcohols can reduce the viscosity of bio-oil to make it exhibits good fluidity.

    Table 6 Major chemical compositions of bio-oils from HTL of polysaccharides-proteins

    ()

    ()

    Category CompoundArea (%) No catalystWith catalyst Pyrroles3-Ethyl-2,4,5-trimethyl-1H-pyrrole3.084.01 N-Ethylpyrrole2.231.97 2-Methyl-3,4-diethyl-1H-pyrrole1.532.77 2-Ethyl-3,4,5-trimethyl-1H-pyrrole4.220.84 N-Acetylpyrrolidine–3.13 2-Ethyl-3,4,5-trimethyl-1H-pyrrole–1.26 Phenols2,3,5,6-Tetramethylhydroquinone2.333.38 3-Hydroxymethyl-5-methoxyphenol1.432.64 3-Ethoxy-4-methoxybenzaldehyde3.163.01 5-Methoxy-2-(methoxymethyl)phenol2.854.69 Alcohols2-Ethylhexanol–10.47 3-Methoxybenzyl alcohol1.73– AmidesN-Methylsuccinimide4.632.41 Hexadecanamide6.880.41

    Fig.11 Formation pathways of bio-oils from HTL of algae three major model components in the presence of SO42?/ZrO2 catalyst. Red and blue represents the promoted and inhibited process by the solid acid catalyst.

    Overall, the addition of SO42?/ZrO2catalyst can significantly regulate the composition and content of bio-oil and affect the degradation and interaction pathway of polysaccharides, proteins and lipids (Fig.11). During the HTL process, the reaction pathways can be divided into the following steps: 1) at low temperature stage, the main algae components polysaccharides, proteins and lipids hydrolyzed into small molecules; 2) with the increasing of temperature, the small molecules were further degrad- ed, dehydrated, decarboxylated and deaminated to decompose into active intermediates; 3) as temperature continues to rise, Maillard reaction, amidation reaction, esterification reaction were tended to occur between the active intermediates (Galadima and Muraza, 2018; Wang, 2018; Marzbali, 2021;).The coordination of sulfate radical with metal atoms on the oxide surface provides activity for the SO42?/ ZrO2solid acid (Li, 2018). The addition of solid acids has an important influence on macromolecular hydrolysis during the whole process of producing bio-oil, and controls the further reaction rate of bio-oil formation process, which is beneficial to accelerate the formation of bio-oil from cellular components (Yang., 2017).

    4 Conclusions

    Solid superacid SO42?/ZrO2has been prepared by precipitation-impregnation method and the acid amount showed a decreasing trend as the calcination temperature increased. The prepared SO42?/ZrO2can been seen as excellent catalyst in the HTL ofand, which promote the increasing of light component proportion and the HHVs, while inhibited the Maillard reaction between polysaccharides and proteins to upgrade the bio-oil producing. Based on the analysis of the properties of bio-oil produced by the major model components, it has put forward the formation path- ways of bio-oils from HTL of real algae in the presence of solid superacid catalyst SO42?/ZrO2.

    Acknowledgements

    This research was supported by the Shandong Provincial Natural Science Foundation, China (No. ZR2019BB 033) and the Fundamental Research Funds for the Central Universities of Ocean University of China (No. 201813 031).

    Ambaye, T. G., Vaccari, M., Bonilla-Petriciolet, A., Prasad, S., van Hullebusch, E. D., and Rtimi, S., 2021. Emerging tech- nologies for biofuel production: A critical review on recent progress, challenges and perspectives.,290: 112627, DOI:10.1016/j.jenvman. 2021.112627.

    Biller, P., Riley, R., and Ross, A. B., 2011. Catalytic hydrothermal processing of microalgae: Decomposition and upgrading of lipids., 102 (7): 4841-4848, DOI: 10.1016/j.biortech.2010.12.113.

    Biswas, B., Kumar, A. A., Bisht, Y., Singh, R., Kumar, J., and Bhaskar, T., 2017. Effects of temperature and solvent on hydrothermal liquefaction ofalgae., 242: 344-350, DOI: 10.1016/j. biortech.2017.03.045.

    Carrero, A., Vicente, G., Rodríguez, R., Peso, G. L. D., and Santos, C., 2015. Synthesis of fatty acids methyl esters (FAMEs) frommicroalga using heterogeneous acid catalysts., 97: 119- 124, DOI: 10.1016/j.bej.2015.02.003.

    Cheng, F., Cui, Z., Chen, L., Jarvis, J., Paz, N., Schaub, T.,., 2017. Hydrothermal liquefaction of high- and low-lipid algae: Bio-crude oil chemistry., 206: 278-292, DOI: 10.1016/j.apenergy.2017.08.105.

    Gai, C., Zhang, Y. H., Chen, W. T., Zhang, P., and Dong, Y. P., 2015. An investigation of reaction pathways of hydrothermal liquefaction usingand, 96: 330-339, DOI: 10.1016/j.enconman.2015.02.056.

    Galadima, A., and Muraza, O., 2018. Hydrothermal liquefaction of algae and bio-oil upgrading into liquid fuels: Role of heterogeneous catalysts., 81: 1037-1048, DOI: 10.1016/j.rser.2017.07.034.

    Hu, Y. L., Feng, S. H., Yuan, Z. S., and Xu, C. B., 2017. Investigation of aqueous phase recycling for improving bio-crude oil yield in hydrothermal liquefaction of algae., 239: 151-159, DOI: 10.1016/j.biortech.2017.05. 033.

    Kim, S. S., Ly, H. V., Kim, J., and Choi, J. H., 2013. Thermogravimetric characteristics and pyrolysis kinetics of algasp. biomass., 139: 242-248, DOI: 10.1016/j.biortech.2013.03.192.

    Li, L., Yan, B., Li, H. X., Yu, S. T., Liu, S. W., Yu, H. L.,., 2018. SO42?/ZrO2as catalyst for upgrading of pyrolysis oil by esterification., 226: 190-194, DOI: 10.1016/j.fuel.2018. 04.006.

    López Barreiro, D., Zamalloa, C., Boon, N., Vyverman, W., Ronsse, F., Brilman, W.,., 2013. Influence of strain-specific parameters on hydrothermal liquefaction of microalgae., 146: 463-471, DOI: 10.1016/j.bior- tech.2013.07.123.

    Marzbali, M. H., Kundu, S., Halder, P., Patel, S., Hakeem, I. G., Paz-Ferreiro, J.,., 2021. Wet organic waste treatmenthydrothermal processing: A critical review., 279: 130557, DOI: 10.1016/j.chemosphere.2021.130557.

    Navalon, S., Dhakshinamoorthy, A., Alvaro, M., and Garcia, H., 2016. Metal nanoparticles supported on two-dimensional graphenes as heterogeneous catalysts., 312: 99-148, DOI: 10.1016/j.ccr.2015.12.005.

    Nazari, L., Yuan, Z., Souzanchi, S., Ray, M. B., and Xu, C. B., 2015. Hydrothermal liquefaction of woody biomass in hot- compressed water: Catalyst screening and comprehensive characterization of bio-crude oils., 162: 74-83, DOI: 10. 1016/j.fuel.2015.08.055.

    Parsa, M., Jalilzadeh, H., Pazoki, M., Ghasemzadeh, R., and Abduli M. A., 2017. Hydrothermal liquefaction ofandmacro-algae for biocrude production., 56: 234-244, DOI: 10. 1016/j.biortech.2017.10.059.

    Peng, Y. B., Xu, Y. F., Dearn, K. D., Geng, J., and Hu, X. G., 2018. Noveltribo-catalysis for improved tribological properties of bio-oil model compound., 212: 546-553, DOI: 10.1016/j.fuel.2017.10.080.

    Qin, J. L., Yu, D. H., and Chen, J. D., 2004. Effect of calcination temperature on surface acidity of TiO2/SiO2., 18 (6): 472-474, DOI: 10.16084/j.cnki.issn 1001-3555.2004.06.014.

    Saber, M., Golzary, A., Hosseinpour, M., Takahashi, F., and Yoshikawa, K., 2016. Catalytic hydrothermal liquefaction of microalgae using nanocatalyst., 183: 566-576, DOI: 10.1016/j.apenergy.2016.09.017.

    Saber, M., Nakhshiniev, B., and Yoshikawa, K., 2016. A review of production and upgrading of algal bio-oil., 58: 918-930, DOI: 10.1016/j. rser.2015.12.342.

    Sharma, N., Jaiswal, K. K., Kumar, V., Vlaskin, M. S., Nanda, M., Rautela, I.,., 2021. Effect of catalyst and temperature on the quality and productivity of HTL bio-oil from micro- algae: A review., 174: 810-822, DOI:10. 1016/j.renene.2021.04.147.

    Sheng, L. L., Wang, X., and Yang, X. Y., 2018. Prediction model of biocrude yield and nitrogen heterocyclic compounds analysis by hydrothermal liquefaction of microalgae with model compounds., 247: 14-20, DOI: 10.1016/j.biortech.2017.08.011.

    SundarRajan, P., Gopinath, K. P., Arun, J., GracePavithra, K., Adithya Joseph, A., and Manasa, S., 2021.Insights into valuing the aqueous phase derived from hydrothermal liquefaction., 144: 111019,DOI:10.1016/j.rser.2021.111019.

    Toor, S. S., Rosendahl, L., and Rudolf A., 2011. Hydrothermal liquefaction of biomass: A review of subcritical water technologies., 36 (5): 2328-2342, DOI: 10.1016/j.energy. 2011.03.013.

    Vo, T. K., Kim, S., Ly, H. V., Lee, E. Y., Lee, C. G., and Kim, J., 2017a. A general reaction network and kinetic model of the hydrothermal liquefaction of microalgaesp., 241: 610-619, DOI: 10.1016/j.biortech. 2017.05.186.

    Vo, T. K., Ly, H. V., Lee, O. K., Lee, E. Y., Kim, C. H., Seo, J. E.,., 2017b. Pyrolysis characteristics and kinetics of microalgalsp. KRS101., 118: 369-376, DOI: 10.1016/j.energy.2016.12.040.

    Wang, B., Zhu, J. P., and Ma, H. Z., 2009. Desulfurization from thiophene by SO42?/ZrO2catalytic oxidation at room temperature and atmospheric pressure., 164 (1): 256-264, DOI: 10.1016/j.jhazmat.2008. 08.003.

    Wang, W. J., Yu, Q., Meng, H., Han, W., Li, J., and Zhang, J. L., 2018. Catalytic liquefaction of municipal sewage sludge over transition metal catalysts in ethanol-water co-solvent., 249: 361-367, DOI: 10.1016/j.biortech. 2017.09.205.

    Wang, X., Tang, X. H., and Yang, X. Y., 2017. Pyrolysis mechanism of microalgaesp. based on model compounds and their interaction., 140: 203-210, DOI: 10.1016/j.enconman.2017. 02.058.

    Xu, Y. F., Hu, Y. H., Peng, Y. B., Yao, L. L., Dong, Y. H., Yang, B. X.,., 2020. Catalytic pyrolysis and liquefaction behavior of microalgae for bio-oil production., 300: 122665, DOI: 10.1016/j.biortech.2019.122 665.

    Xu, Y. F., Liu, K., Hu, Y. H., Dong, Y. H., and Yao, L. L., 2020. Experimental investigation and comparison of bio-oil from hybrid microalgaesuper/subcritical liquefaction., 279: 118412, DOI: 10.1016/j.fuel.2020.118412.

    Yang, W. C., Li, X. G., Li, Z. H., Huang, F. B., Sun, X. M., and Feng, L. J., 2015. Application of solid acid SO42?/Fe2O3in thermal liquefaction of., 31 (5): 1075-1081, DOI: 10.3969/j.issn.1001-8719.2015.05.007.

    Yang, W. C., Li, X. G., Li, Z. H., Tong, C. H., and Feng, L. J., 2015. Understanding low-lipid algae hydrothermal liquefaction characteristics and pathways through hydrothermal liquefaction of algal major components: Crude polysaccharides, crude proteins and their binary mixtures., 196: 99-108, DOI: 10.1016/j.biortech.2015.07.020.

    Yang, W. C., Li, X. G., Zhang, D. H., and Feng, L. J., 2017. Catalytic upgrading of bio-oil in hydrothermal liquefaction of algae major model components over liquid acids., 154: 336-343, DOI: 10.1016/j. enconman.2017.11.018.

    Yang, W. C., Wang, Z. W., Han, J. B., Song, S., Zhang, Y., and Gong, W. M., 2019. The role of polysaccharides and proteins in bio-oil production during the hydrothermal liquefaction of algae species., 9: 41962-41969, DOI: 10.1039/ c9ra07150d.

    Yuan, Z., Cheng, S., Leitch, M., and Xu, C. B., 2010. Hydrolytic degradation of alkaline lignin in hot-compressed water and ethanol., 101 (23): 9308-9313, DOI: 10.1016/j.biortech.2010.06.140.

    Zhao, Y. N., Feng, L. J., Liu, S. S., Li, D., and Yang, W. C., 2017. Application of solid acid SO42?/ZrO2in ethanol thermal liquefaction ofand., 33 (3): 447-455, DOI: 10.3969/j.issn.1001- 8719.2017.03.008.

    Zhu, Z. B., Si, B. C., Lu, J. W., Watson, J., Zhang, Y. H., and Liu, Z. D., 2017. Elemental migration and characterization of products during hydrothermal liquefaction of cornstalk., 24: 9-16, DOI: 10.1016/j.biortech.2017. 06.085.

    June 30, 2021;

    April 11, 2022;

    April 28, 2022

    ? Ocean University of China, Science Press and Springer-Verlag GmbH Germany 2022

    . E-mail: yaoshuo@ouc.edu.cn

    (Edited by Ji Dechun)

    日产精品乱码卡一卡2卡三| 99视频精品全部免费 在线| av视频免费观看在线观看| 亚洲国产毛片av蜜桃av| 国产日韩欧美在线精品| 自线自在国产av| 久久久a久久爽久久v久久| tube8黄色片| 一区二区日韩欧美中文字幕 | 91精品伊人久久大香线蕉| 久久午夜福利片| 能在线免费看毛片的网站| 丝袜美足系列| 久久国产亚洲av麻豆专区| 国产免费视频播放在线视频| 亚洲精品久久午夜乱码| 国产精品久久久久久av不卡| 婷婷色综合www| 免费观看无遮挡的男女| 国产免费视频播放在线视频| 国产精品熟女久久久久浪| 国产日韩欧美视频二区| 美女视频免费永久观看网站| 丁香六月天网| 日韩电影二区| 久久精品久久久久久久性| 夜夜骑夜夜射夜夜干| 五月开心婷婷网| 久久国产精品大桥未久av| 大话2 男鬼变身卡| 母亲3免费完整高清在线观看 | 天堂8中文在线网| 久久国产亚洲av麻豆专区| 精品久久久精品久久久| 国产在视频线精品| 九色亚洲精品在线播放| 性色av一级| 丝袜在线中文字幕| 自线自在国产av| 精品亚洲乱码少妇综合久久| 日日啪夜夜爽| 久久久a久久爽久久v久久| 男人添女人高潮全过程视频| 免费看不卡的av| 亚洲在久久综合| www.av在线官网国产| 老司机影院毛片| 人妻夜夜爽99麻豆av| 蜜臀久久99精品久久宅男| 日韩中字成人| 一级毛片 在线播放| 国产亚洲av片在线观看秒播厂| 国产欧美日韩一区二区三区在线 | 美女主播在线视频| 九九久久精品国产亚洲av麻豆| 美女主播在线视频| av免费在线看不卡| 国产视频首页在线观看| 亚洲色图综合在线观看| 欧美日韩国产mv在线观看视频| 成年人午夜在线观看视频| 超碰97精品在线观看| 丝袜喷水一区| 日韩一区二区视频免费看| a级片在线免费高清观看视频| 国产在线免费精品| 亚洲欧美精品自产自拍| 久久久午夜欧美精品| 精品久久久久久久久亚洲| 亚洲精品色激情综合| 欧美日韩视频高清一区二区三区二| 国模一区二区三区四区视频| 人妻少妇偷人精品九色| 高清毛片免费看| 青春草视频在线免费观看| 午夜福利,免费看| 黑丝袜美女国产一区| 日韩一本色道免费dvd| 搡女人真爽免费视频火全软件| 中文精品一卡2卡3卡4更新| 日日摸夜夜添夜夜爱| 欧美精品亚洲一区二区| 亚洲国产成人一精品久久久| 日韩视频在线欧美| 国产精品不卡视频一区二区| 色94色欧美一区二区| 一级,二级,三级黄色视频| 最新中文字幕久久久久| 色吧在线观看| 美女xxoo啪啪120秒动态图| 好男人视频免费观看在线| 在线 av 中文字幕| 建设人人有责人人尽责人人享有的| 精品亚洲乱码少妇综合久久| 色5月婷婷丁香| 在线天堂最新版资源| av在线观看视频网站免费| 爱豆传媒免费全集在线观看| 久久97久久精品| 国产精品熟女久久久久浪| 美女视频免费永久观看网站| tube8黄色片| 99re6热这里在线精品视频| 插逼视频在线观看| 汤姆久久久久久久影院中文字幕| 成人国产麻豆网| 中文精品一卡2卡3卡4更新| 人人妻人人添人人爽欧美一区卜| 18禁在线无遮挡免费观看视频| 麻豆乱淫一区二区| 涩涩av久久男人的天堂| 亚洲精华国产精华液的使用体验| 五月开心婷婷网| 亚洲av欧美aⅴ国产| 精品国产乱码久久久久久小说| 亚洲综合色网址| 色视频在线一区二区三区| 亚洲av免费高清在线观看| 欧美亚洲日本最大视频资源| 欧美精品高潮呻吟av久久| 久久久久久久久久成人| 插逼视频在线观看| 超色免费av| 国产欧美另类精品又又久久亚洲欧美| av在线播放精品| 久久久久精品性色| 国产欧美日韩一区二区三区在线 | xxx大片免费视频| 最近中文字幕高清免费大全6| 在线 av 中文字幕| 熟妇人妻不卡中文字幕| kizo精华| 夜夜看夜夜爽夜夜摸| 黄片无遮挡物在线观看| 久久久久久久久大av| 亚洲av电影在线观看一区二区三区| 成人无遮挡网站| 蜜桃国产av成人99| 日本欧美视频一区| 久久精品国产亚洲av涩爱| 亚洲av日韩在线播放| av免费观看日本| 男女边吃奶边做爰视频| 美女内射精品一级片tv| 热re99久久精品国产66热6| 国产免费又黄又爽又色| 国产精品99久久久久久久久| 精品国产一区二区三区久久久樱花| 91国产中文字幕| 一边亲一边摸免费视频| 久久青草综合色| 免费观看的影片在线观看| 婷婷色av中文字幕| 久久婷婷青草| 大陆偷拍与自拍| 国产黄频视频在线观看| 国产一区亚洲一区在线观看| 啦啦啦视频在线资源免费观看| 国产 精品1| 中文字幕人妻熟人妻熟丝袜美| 精品人妻熟女毛片av久久网站| 人成视频在线观看免费观看| 一区二区av电影网| 激情五月婷婷亚洲| 老司机影院毛片| av免费观看日本| 男女国产视频网站| 2021少妇久久久久久久久久久| 美女视频免费永久观看网站| 80岁老熟妇乱子伦牲交| 久久久久久久久久人人人人人人| 国产极品天堂在线| 九色成人免费人妻av| 性色avwww在线观看| 久热久热在线精品观看| 成人午夜精彩视频在线观看| 久久女婷五月综合色啪小说| 2022亚洲国产成人精品| 一级片'在线观看视频| videossex国产| 久久久精品区二区三区| 久久婷婷青草| 中文字幕人妻丝袜制服| 观看av在线不卡| 91久久精品国产一区二区三区| 婷婷色综合www| 亚洲美女搞黄在线观看| 欧美日本中文国产一区发布| freevideosex欧美| 九色成人免费人妻av| 最后的刺客免费高清国语| 国产亚洲一区二区精品| 欧美一级a爱片免费观看看| 99视频精品全部免费 在线| 七月丁香在线播放| 免费观看的影片在线观看| 黑人巨大精品欧美一区二区蜜桃 | 欧美+日韩+精品| 亚洲精品亚洲一区二区| 日本猛色少妇xxxxx猛交久久| 亚洲精品日本国产第一区| 高清在线视频一区二区三区| 伦理电影大哥的女人| 午夜av观看不卡| 亚洲情色 制服丝袜| 亚洲五月色婷婷综合| 亚洲精华国产精华液的使用体验| 亚洲精品自拍成人| 夫妻午夜视频| 老司机亚洲免费影院| 日本91视频免费播放| 夫妻性生交免费视频一级片| 国产 一区精品| 免费大片18禁| 少妇熟女欧美另类| 成年女人在线观看亚洲视频| 免费观看av网站的网址| 日本vs欧美在线观看视频| 热re99久久国产66热| 18禁观看日本| 国产精品99久久久久久久久| 亚洲精品,欧美精品| 欧美国产精品一级二级三级| 久久久久网色| 久久久a久久爽久久v久久| 久久国产精品大桥未久av| 日韩精品有码人妻一区| 极品人妻少妇av视频| 国产综合精华液| 最近2019中文字幕mv第一页| 成人综合一区亚洲| 精品一品国产午夜福利视频| 成人亚洲精品一区在线观看| 久久人妻熟女aⅴ| 观看av在线不卡| 色视频在线一区二区三区| 性色av一级| a级片在线免费高清观看视频| 国产精品 国内视频| 丰满乱子伦码专区| 男人爽女人下面视频在线观看| 亚洲国产毛片av蜜桃av| 亚洲第一区二区三区不卡| 99热网站在线观看| 亚洲精品av麻豆狂野| 国产精品久久久久久精品古装| 久久久精品免费免费高清| 两个人的视频大全免费| 啦啦啦在线观看免费高清www| 国产精品熟女久久久久浪| 国产免费又黄又爽又色| 高清午夜精品一区二区三区| 一级毛片 在线播放| av不卡在线播放| 老女人水多毛片| 午夜福利在线观看免费完整高清在| 国产精品人妻久久久影院| 免费日韩欧美在线观看| 男人添女人高潮全过程视频| 两个人的视频大全免费| 国产亚洲欧美精品永久| 国产成人aa在线观看| 少妇的逼水好多| 久久免费观看电影| 免费av中文字幕在线| 一边摸一边做爽爽视频免费| a级片在线免费高清观看视频| 亚洲av国产av综合av卡| 国国产精品蜜臀av免费| 日韩一区二区三区影片| 国产在线免费精品| 成人手机av| 一级二级三级毛片免费看| 亚洲av综合色区一区| 国产一区二区三区综合在线观看 | 69精品国产乱码久久久| 日本爱情动作片www.在线观看| 狠狠精品人妻久久久久久综合| 欧美最新免费一区二区三区| 亚洲国产精品专区欧美| 亚洲成人手机| 精品亚洲成国产av| 欧美xxⅹ黑人| 精品久久久久久久久av| 久久久久国产精品人妻一区二区| 久久久欧美国产精品| 国产精品国产三级专区第一集| 丰满迷人的少妇在线观看| 在线观看www视频免费| av在线观看视频网站免费| 有码 亚洲区| a级毛色黄片| 免费看光身美女| 少妇被粗大猛烈的视频| 亚洲,一卡二卡三卡| 亚洲av在线观看美女高潮| 国产男女超爽视频在线观看| 桃花免费在线播放| 热99国产精品久久久久久7| 亚洲第一av免费看| 亚洲av欧美aⅴ国产| 国产视频首页在线观看| 超碰97精品在线观看| 国产成人精品在线电影| 一级毛片 在线播放| 午夜激情av网站| a级片在线免费高清观看视频| 亚洲天堂av无毛| 一二三四中文在线观看免费高清| 亚洲av在线观看美女高潮| 亚洲av欧美aⅴ国产| 青春草亚洲视频在线观看| 国产视频内射| 欧美日本中文国产一区发布| 少妇人妻 视频| 2021少妇久久久久久久久久久| 两个人免费观看高清视频| 亚洲色图 男人天堂 中文字幕 | 桃花免费在线播放| 日韩欧美精品免费久久| 天天影视国产精品| 亚洲精品亚洲一区二区| 青春草亚洲视频在线观看| 中文字幕免费在线视频6| 性高湖久久久久久久久免费观看| av免费在线看不卡| 欧美日韩av久久| 高清视频免费观看一区二区| 国产精品秋霞免费鲁丝片| av女优亚洲男人天堂| 久久久久久久久久成人| 又大又黄又爽视频免费| 日本欧美国产在线视频| h视频一区二区三区| 如日韩欧美国产精品一区二区三区 | 亚洲成人av在线免费| 大片电影免费在线观看免费| 久久久精品区二区三区| 国产一区二区三区综合在线观看 | 国产男人的电影天堂91| 麻豆乱淫一区二区| 性高湖久久久久久久久免费观看| 国产精品人妻久久久影院| 麻豆乱淫一区二区| 日韩熟女老妇一区二区性免费视频| 最近最新中文字幕免费大全7| 在现免费观看毛片| 国产精品国产av在线观看| 国产一区二区在线观看日韩| 一级毛片我不卡| 80岁老熟妇乱子伦牲交| 日韩中字成人| 亚洲av免费高清在线观看| 欧美人与善性xxx| 亚洲欧美成人综合另类久久久| 日韩免费高清中文字幕av| 少妇被粗大的猛进出69影院 | 精品国产一区二区三区久久久樱花| 免费久久久久久久精品成人欧美视频 | 亚洲色图 男人天堂 中文字幕 | xxx大片免费视频| 丰满饥渴人妻一区二区三| 岛国毛片在线播放| 最近最新中文字幕免费大全7| 亚洲欧美一区二区三区国产| 国产精品99久久久久久久久| 精品国产一区二区久久| 下体分泌物呈黄色| 日韩,欧美,国产一区二区三区| 亚洲精品色激情综合| 亚洲精品中文字幕在线视频| 99精国产麻豆久久婷婷| av一本久久久久| 99久久人妻综合| 伊人亚洲综合成人网| 少妇猛男粗大的猛烈进出视频| 午夜免费观看性视频| 日韩中字成人| 老司机亚洲免费影院| 久久久久久久久久成人| 日日撸夜夜添| 亚洲精品第二区| 免费av不卡在线播放| 黑人猛操日本美女一级片| 特大巨黑吊av在线直播| 寂寞人妻少妇视频99o| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜添av毛片| 欧美3d第一页| 中国美白少妇内射xxxbb| 99久久中文字幕三级久久日本| 国产黄色视频一区二区在线观看| 伦精品一区二区三区| 午夜激情久久久久久久| 天堂8中文在线网| 久久久久久久亚洲中文字幕| av播播在线观看一区| 免费日韩欧美在线观看| 国产在线免费精品| 夫妻午夜视频| 乱人伦中国视频| 婷婷色麻豆天堂久久| 一级毛片我不卡| 99国产精品免费福利视频| 又粗又硬又长又爽又黄的视频| 嘟嘟电影网在线观看| 高清黄色对白视频在线免费看| 91精品国产国语对白视频| 80岁老熟妇乱子伦牲交| 久久久久久久久大av| 久久久精品区二区三区| 亚洲第一区二区三区不卡| 最近中文字幕高清免费大全6| 97在线人人人人妻| 97在线视频观看| 边亲边吃奶的免费视频| 18禁在线播放成人免费| 在线观看人妻少妇| 只有这里有精品99| 女性被躁到高潮视频| 日日摸夜夜添夜夜添av毛片| 亚洲综合色网址| 亚洲精品乱码久久久v下载方式| 国产黄色免费在线视频| 亚洲欧美清纯卡通| 精品少妇黑人巨大在线播放| 国产精品免费大片| 久久久久久久久久久丰满| 国产精品成人在线| 热99国产精品久久久久久7| 伦理电影大哥的女人| 男女国产视频网站| 精品亚洲乱码少妇综合久久| 久久99精品国语久久久| 五月天丁香电影| 国产高清不卡午夜福利| 久久久久久久精品精品| 久久久精品免费免费高清| 十分钟在线观看高清视频www| 亚洲国产欧美日韩在线播放| 欧美精品人与动牲交sv欧美| 亚洲高清免费不卡视频| 日日摸夜夜添夜夜添av毛片| 久久久国产一区二区| 一边摸一边做爽爽视频免费| 国产毛片在线视频| 如日韩欧美国产精品一区二区三区 | 亚洲av成人精品一区久久| 永久免费av网站大全| 一边亲一边摸免费视频| 狂野欧美白嫩少妇大欣赏| 在线播放无遮挡| 欧美日韩视频精品一区| 久久综合国产亚洲精品| 特大巨黑吊av在线直播| 女性生殖器流出的白浆| 一本—道久久a久久精品蜜桃钙片| 狂野欧美激情性bbbbbb| 免费日韩欧美在线观看| 国产精品人妻久久久久久| 国产精品国产三级国产av玫瑰| 亚洲国产成人一精品久久久| 汤姆久久久久久久影院中文字幕| 99久久精品国产国产毛片| 亚洲av二区三区四区| 色视频在线一区二区三区| 久久久精品免费免费高清| 中文精品一卡2卡3卡4更新| 在线观看人妻少妇| 亚洲美女视频黄频| 免费看av在线观看网站| 亚洲精品aⅴ在线观看| 好男人视频免费观看在线| 久久久久国产网址| 伦理电影免费视频| 男女无遮挡免费网站观看| 桃花免费在线播放| 国产欧美日韩一区二区三区在线 | 国产黄色视频一区二区在线观看| 97超视频在线观看视频| 久久久久久人妻| 日本午夜av视频| 亚洲成人一二三区av| 男女边摸边吃奶| 亚洲国产精品成人久久小说| 久久久久久久久久人人人人人人| 高清不卡的av网站| 亚洲国产精品999| √禁漫天堂资源中文www| 亚洲av成人精品一区久久| 女人精品久久久久毛片| 国产精品蜜桃在线观看| 久久99蜜桃精品久久| 亚洲五月色婷婷综合| 亚洲av欧美aⅴ国产| 热99国产精品久久久久久7| 大陆偷拍与自拍| 国精品久久久久久国模美| 亚洲国产精品999| 久久午夜综合久久蜜桃| 赤兔流量卡办理| 精品久久久久久久久av| 精品视频人人做人人爽| xxx大片免费视频| 十八禁高潮呻吟视频| 欧美人与性动交α欧美精品济南到 | 久久久久久人妻| 欧美bdsm另类| 晚上一个人看的免费电影| 亚洲色图 男人天堂 中文字幕 | 王馨瑶露胸无遮挡在线观看| 五月开心婷婷网| 免费观看无遮挡的男女| 欧美+日韩+精品| 视频中文字幕在线观看| 亚洲性久久影院| 青春草亚洲视频在线观看| 熟女av电影| 国产高清国产精品国产三级| 久久久精品免费免费高清| 国产精品免费大片| 久久狼人影院| 国产亚洲欧美精品永久| 久久久久久久亚洲中文字幕| 亚洲一级一片aⅴ在线观看| 91成人精品电影| 久久久久久伊人网av| 久久精品国产a三级三级三级| 亚洲美女视频黄频| 亚洲av不卡在线观看| 日韩中字成人| 午夜精品国产一区二区电影| 美女主播在线视频| 欧美成人精品欧美一级黄| 国产熟女午夜一区二区三区 | 中文欧美无线码| 婷婷成人精品国产| av免费观看日本| 蜜桃久久精品国产亚洲av| 美女主播在线视频| 纯流量卡能插随身wifi吗| 国产欧美日韩综合在线一区二区| 自线自在国产av| 99久久中文字幕三级久久日本| 高清在线视频一区二区三区| 亚洲欧美成人综合另类久久久| 欧美成人午夜免费资源| 欧美人与性动交α欧美精品济南到 | 日本欧美视频一区| 国产成人精品婷婷| 秋霞伦理黄片| 男女高潮啪啪啪动态图| 国产毛片在线视频| 高清不卡的av网站| av有码第一页| 热re99久久国产66热| 久久国产精品大桥未久av| 最近手机中文字幕大全| 欧美三级亚洲精品| 99久久精品国产国产毛片| 纵有疾风起免费观看全集完整版| 曰老女人黄片| 亚洲激情五月婷婷啪啪| 51国产日韩欧美| 久久久午夜欧美精品| 久久99热6这里只有精品| 成人黄色视频免费在线看| 特大巨黑吊av在线直播| 国产精品秋霞免费鲁丝片| 十八禁高潮呻吟视频| 在线观看三级黄色| 成年女人在线观看亚洲视频| 欧美少妇被猛烈插入视频| 欧美+日韩+精品| 亚洲成人手机| 三级国产精品欧美在线观看| 51国产日韩欧美| 制服诱惑二区| 久久久久久人妻| 亚洲一级一片aⅴ在线观看| 男女无遮挡免费网站观看| 中文字幕最新亚洲高清| 国产在线免费精品| 99九九在线精品视频| 全区人妻精品视频| 美女国产高潮福利片在线看| 国产精品国产三级专区第一集| 青春草国产在线视频| 看非洲黑人一级黄片| 在线观看人妻少妇| 国产精品三级大全| 欧美激情国产日韩精品一区| 美女福利国产在线| 大陆偷拍与自拍| 我要看黄色一级片免费的| 中文天堂在线官网| 欧美3d第一页| 九九久久精品国产亚洲av麻豆| av在线老鸭窝| 国产精品三级大全| 人人妻人人爽人人添夜夜欢视频| 啦啦啦中文免费视频观看日本| 国产成人精品在线电影| 一本一本综合久久| 亚洲av国产av综合av卡| 国产精品久久久久久久久免| 极品人妻少妇av视频| 亚洲三级黄色毛片| 成人18禁高潮啪啪吃奶动态图 | 欧美精品一区二区免费开放| 高清毛片免费看| 午夜激情久久久久久久| 在线播放无遮挡| 91精品伊人久久大香线蕉| 成人毛片60女人毛片免费|