999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

Assessing the Efficacy of Improved Learning in Hourly Global Irradiance Prediction

2023-12-15 03:58:12AbdennasserDahmaniYaminaAmmiNadjemBailekAlbanKuriqiNadhirAlAnsariSalahHaniniIlhamiColakLaithAbualigah9andElSayedElkenawy5
Computers Materials&Continua 2023年11期

Abdennasser Dahmani,Yamina Ammi,Nadjem Bailek,Alban Kuriqi,Nadhir Al-Ansari,Salah Hanini,Ilhami Colak,Laith Abualigah9,0,,,3,4 and El-Sayed M.El-kenawy5

1Department of Mechanical Engineering,Faculty of Science and Technology,GIDD Industrial Engineering and Sustainable Development Laboratory,University of Relizane,Bourmadia,Relizane,48000,Algeria

2Laboratory of Biomaterials and Transport Phenomena(LBMPT),University of Medea,Urban Pole,Medea,26000,Algeria 3Energies and Materials Research Laboratory,Faculty of Sciences and Technology,University of Tamanghasset,Tamanrasset,10034,Algeria

4Sustainable Development and Computer Science Laboratory,Faculty of Sciences and Technology,Ahmed Draia University of Adrar,Adrar,Algeria

5CERIS,Instituto Superior Técnico,Universidade de Lisboa,Lisbon,1049-001,Portugal

6Civil Engineering Department,University for Business and Technology,Pristina,10000,Kosovo

7Department of Civil,Environmental and Natural Resources Engineering,Lulea University of Technology,Lulea,97187,Sweden 8Engineering and Architecture Faculty,Nisantasi University,Istanbul,Turkey

9Computer Science Department,Prince Hussein Bin Abdullah Faculty for Information Technology,Al al-Bayt University,Mafraq,25113,Jordan

10Hourani Center for Applied Scientific Research,Al-Ahliyya Amman University,Amman,19328,Jordan

11MEU Research Unit,Middle East University,Amman,11831,Jordan

12Applied Science Research Center,Applied Science Private University,Amman,11931,Jordan

13School of Computer Sciences,Universiti Sains Malaysia,Pulau Pinang,11800,Malaysia

14Department of Computing and Information Systems,School of Engineering and Technology,Sunway University Malaysia,Petaling Jaya,27500,Malaysia

15Faculty of Artificial Intelligence,Delta University for Science and Technology,Mansoura,35712,Egypt

ABSTRACT Increasing global energy consumption has become an urgent problem as natural energy sources such as oil,gas,and uranium are rapidly running out.Research into renewable energy sources such as solar energy is being pursued to counter this.Solar energy is one of the most promising renewable energy sources,as it has the potential to meet the world’s energy needs indefinitely.This study aims to develop and evaluate artificial intelligence (AI)models for predicting hourly global irradiation.The hyperparameters were optimized using the Broyden-Fletcher-Goldfarb-Shanno(BFGS)quasi-Newton training algorithm and STATISTICA software.Data from two stations in Algeria with different climatic zones were used to develop the model.Various error measurements were used to determine the accuracy of the prediction models,including the correlation coefficient,the mean absolute error,and the root mean square error(RMSE).The optimal support vector machine(SVM)model showed exceptional efficiency during the training phase,with a high correlation coefficient(R=0.99)and a low mean absolute error(MAE=26.5741 Wh/m2),as well as an RMSE of 38.7045 Wh/m2 across all phases.Overall,this study highlights the importance of accurate prediction models in the renewable energy,which can contribute to better energy management and planning.

KEYWORDS Renewable energy;energy prediction;global irradiation;artificial intelligence;BFGS quasi-Newton training algorithm

1 Introduction

Increasing concern about the effects of climate change and the need to diversify energy sources has led to a significant increase in the development of renewable energy sources.Among these,solar energy has emerged as a viable option due to its abundance and potential to reduce carbon emissions.Moreover,as oil and gas resources become less available,developing renewable energy sources is becoming increasingly important for the country’s long-term energy security.According to the recent the intergovermental panel on climate change(IPCC)report,solar energy has the potential to meet a significant portion of the world’s energy needs,and Algeria is no exception[1].

Accurately estimating the amount of solar radiation hitting the Earth’s surface is critical for various applications such as photovoltaic systems,heating,medical research,agriculture,and architecture.This is usually done with solar measurement devices such as solarimeters or pyranometers.However,it is difficult to measure solar radiation in many places in Algeria because the meters are too expensive,and the systems are very complex.Even if there are several meteorological stations in different locations in Algeria,measurements may not always be available due to power outages or limitations on the number of variables that can be recorded[2].

To address these challenges,researchers have developed models that use readily available meteorological data to predict global solar radiation (GSR) more accurately.These predictive models are becoming more advanced daily,but results vary by location.Therefore,it is important to use sophisticated GSR prediction techniques to improve solar energy potential prediction accuracy in Algeria[3].

Many research efforts have been made to predict solar radiation (SR) in different areas of the world using various techniques such as artificial intelligence and empirical methods.One popular method is Multilayer Perceptron (ANN-MLP) artificial neural network technology.However,other methods,such as decision tree models,support vector machines (SVMs),and feed-forward radial basis functions(FF-RBFs),have also been used to estimate solar radiation.Researchers such as[4-10] have used SVMs,and others have also used FF-RBF.Loutfi et al.[11] presented three different Feed-Forward Neural Network (FFNN) model topologies for generating global,direct and diffuse hourly solar radiation in Fez,Morocco.In order to perform a comparative study,different models were implemented,including the decision tree model,random forest model,generalized linear models,artificial neural network,linear regression model,and adaptive fuzzy neural interference system model[12].Bamisile et al.[13]developed and compared eight artificial intelligence models for solar radiation prediction at different time intervals (hourly,every minute,and daily average) using datasets from 6 African countries.They found that different AI models suited different solar radiation estimation tasks.Extreme gradient boost algorithm (XG-Boost) boost was the best model for 10 of the 13 case studies considered in this work.They concluded that hourly solar irradiance prediction is more accurate for the models than the daily average and minute time step.

Solar energy is a focus of Algeria’s ambitious energy policy,which allocates significant resources to solar thermal and photovoltaic resources.Projections indicate that solar energy will account for more than 37%of the country’s electricity generation by 2030[14].The annual sunshine duration in Algeria is more than 3,900 h in the Sahara and 3,000 h on the plateaus.The daily energy gain on a horizontal surface of 1 m2averages 5 KWh [10],in most regions of the country [15].Thus,solar energy is a good basis to help the country meet its energy needs.However,accurate solar radiation estimates are needed to exploit this potential fully.In this study,we use state-of-the-art machine learning techniques to improve the predictability of solar energy potential in Algeria.

The main objective of this research is to develop a method to optimize the hyperparameters of traditional machine learning using the multilayer perceptron (MLP) and support vector machines,thus increasing the reliability of hourly predictions of global irradiance.We used the FNN-MLP and SVM models to produce a reliable forecast of global solar irradiance at one-hour intervals at stations with different climates in Algeria.The following overview provides the framework for this research work: the materials and processes are discussed in Section 2,while the construction of the model is covered in Section 3,and the results and discussion are presented in Section 4.The paper concluded with a conclusion.

2 Materials and Methods

2.1 Studied Region and Database Collection

In this investigation,two radiometric stations were used to compile the database.The first station,“Shems”,located in Bouzareah in Algeria,recorded experimental data using Kipp and Zonen pyranometers to measure the global horizontal irradiation (GHI).The second station,located in Tamanrasset in the Sahara Desert in southern Algeria,is equipped with an Eppley PSP pyranometer and has the highest solar energy resources in an arid desert environment.Table 1 presents comprehensive information on the two stations used for training and testing purposes,including station ID,station name,latitude,longitude,elevation,climate zone,and data periods.The table indicates that the training station,Bouzareah(BOU),is located in the Mediterranean climate zone and covers the period between January 01,2014,and December 31,2014.In contrast,the testing station,Tamanrasset(TAM),is in the hot desert climate zone and covers the period from July 01,2019,to December 17,2020.The information in the table is crucial for understanding the data used in the investigation and interpreting the results obtained.

Table 1: Geographical region and period covered by stations in this investigation

Table 2 provides statistical analysis results for the input and output variables of solar radiation prediction at the BOU and TAM stations for temperature (TMP),humidity (HUM),wind speed(WSP),and global horizontal irradiance at the BOU and TAM stations.The BOU station has a maximum temperature of 44.11°C,a minimum humidity of 8.28%,a mean wind speed of 4.62 m/s,and a mean GHI of 517.88 Wh/m2.The TAM station has a maximum temperature of 38.5°C,a minimum humidity of 2%,a mean wind speed of 5.19 m/s,and a mean GHI of 678.76 Wh/m2.The correlation between these variables and solar radiation is further explored in Table 3.The results show that temperature positively correlates with solar radiation at both stations,with humidity having a negative correlation at BOU and a weak negative correlation at TAM.Wind speed has a weak negative correlation at BOU and a positive correlation at TAM.These correlations can be useful for developing accurate solar irradiance prediction models and designing effective solar energy systems.

Table 2:Statistical analysis of input and output variables for solar irradiance prediction at two stations

Table 3: Climatic-output correlations for solar irradiance prediction at two stations

Fig.1 displays the frequency counts of GHI between the two stations,BOU and TAM,at different bin center intervals.The plot reveals that TAM has higher frequency counts of GHI than BOU at all bin center intervals.The highest frequency counts for TAM are at the bin centers of 250,350,and 450 Wh/m2,with 508,254,and 192 counts,respectively.In contrast,BOU has its highest frequency counts at the bin centers of 250,450,and 550 Wh/m2,with 508,490,and 458 counts,respectively.These results indicate that TAM receives more GHI than BOU due to their different geographic locations and climatic conditions.Such analysis is valuable for understanding the variability of GHI and designing solar energy systems.

Figure 1:Graphical depiction of the GHI as a function of the relative frequency of each station

2.2 Feedforward Neural Networks Multi-Layer Perceptron(FNN-MLP)

Feedforward Neural Networks’multilayer perceptron,also known as FNN-MLP,is modeled on the human brain’s information processing.Their known ability to learn from their environment makes them ideal for nonlinear modeling systems that are difficult to characterize analytically.Even though the architecture allows for arbitrarily small approximation errors with related weight values,there is still an obstacle to their efficiency in some form of training.The multilayer perceptron,whose architecture defines multiple layers of neurons,is today’s most widely used supervised neural network for approximation problems.

The“FNN-MLP”consists of three layers:an input layer,a hidden layer,and an output layer.Each of these layers hides information from the other two layers.Synaptic weights,Wij.connect the neurons in a layer to the neurons below it.These weights determine the relative importance of each input to the output of each neuron.An activation function ensures that each neuron picks up all the information transmitted by the neurons that preceded it in the layered structure.After this step,an output signal is generated,ready for transmission to the neurons in the subsequent layer[11,16].An“FNN-MLP”with three layers is shown in Fig.2 and is used to predict GHI.

wI,andwh:connection weights(input-hidden,hidden-output)

bh,andbo:columns vector of neuron bias hidden and output,respectively.

The following Eqs.(4) and (5) represent the assimilation of the GHI in an accurate model that includes all inputs xi.

The instance outputsZjof the hidden layer:

The following mathematical formula,which accounts for all inputs and represents the global solar radiation,is produced when Eqs.(1)and(2)are combined.

The “FNN-MLP” framework was optimized for GHI prediction using MATLAB 2020b.The methods,database distribution,layer depth,neuron count,and activation functions are all included.Table 4 displays the optimized FNN-MLP model’s structure.

Figure 2:Three-layer feed-forward neural network multilayer perceptron

Table 4: Structure of the improved“FNN-MLP”model

2.3 Support Vector Machines(SVM)

The Support Vector Machine,commonly known as SVM,is a supervised learning method that has gained significant popularity recently for its ability to predict meteorological data such as temperature[17]and wind speed[18].The ease of use and adaptability of the SVM method makes it suitable for a wide range of classification and regression problems across various sectors,including mechanical engineering,energy,finance [19],and more.Despite its potential use in studies with small sample sets,the SVM method has been shown to provide balanced predicted performance due to its unique characteristics [20].In an SVM model,the regression function can model nonlinear relationships between input and output.The output of an SVM model can be determined by solving a specific equation[21]:

f(xi):The predicted data.

?(xi):The implicitly constructed nonlinear function.

ω:The SVM model’s weight vector.

b:The SVM model’s bias.

The dataset has a D-dimensional input vector xi∈RDand a scalar output yi∈R.

The following equations provide the SVM optimization model(for the training set):

C:The factor that balances model complexity with empirical risk‖w‖2

The problem above can be solved in the same manner as a standard nonlinear restricted optimization problem by utilizing the concepts of Lagrange multipliers to generate a dual optimization problem:

2.4 Model Development

Two models were used to develop an accurate prediction of hourly global solar irradiance: The FNN-MLP and the SVM.The process used to evaluate and improve the structure of the FNN-MLP and SVM models is shown in detail in Fig.3.The datasets were divided into different subsets for each model.The FNN-MLP model’s datasets were divided into three subsets:The training,validation,and testing phases.For the SVM model,on the other hand,the datasets were divided into two subsets,the training phase,and the testing phase.Both subsets were created from the entire data set.

Figure 3:Flow diagram for the two models(FNN-MLP)and(SVM)development

Several techniques were used to determine the most effective FNN-MLP model.The BFGS quasi-Newton training algorithm(trainbfg)was used,four activation functions(log sigmoid,tangent sigmoid,exponential,and sin)were used in the hidden layer,and a single transfer function(identity)was used.Experiments were performed with different sizes for the hidden layers(from 3 to 25 neurons)to obtain the most accurate model possible.An iterative testing process was performed to determine the FNN-MLP model with the best performance.An optimal SVM model was developed using the support vector machine learning strategy for the SVM technique.The selection of appropriate kernel features is critical to the success of the SVM model.The STATISTICA software provides a wide range of kernel functions for SVM models.The penalty term for the Gaussian radial basis function parameters was set to nu=1.0000,C=10.0000,and Gamma=13.93.This process determined the optimal values for the target parameters of the SVM model.

2.5 Evaluation Criteria

In this study,various error measures were employed to determine the level of accuracy of the prediction models.These error measures include the Correlation Coefficient (R),mean percentage error(MPE),and Root Mean Squared Error(RMSE).These measures are mathematically represented by Eqs.(8)-(10)as described in[22-32].These error measures allow for a comprehensive evaluation of the performance of the prediction models,providing a clear understanding of their strengths and weaknesses.

nis the number of data points;Yi,expandYi,calare the experimental and calculated data points of global solar radiation,respectively;andis the mean experimental data.

3 Results and Discussion

This subsection presents the results of the models developed in the study to predict hourly global irradiation.Initially,data collected from the Bouzareah station generated two models: The FNNMLP and SVM.The performance of these models was evaluated using three different data divisions for training,validation,and testing.The results were visualized in Fig.4,which depicts the correlation coefficient(R)error values obtained for each division.It can be observed that division 3 outperforms the other two divisions in terms of R-values for the testing phase,yielding R=0.9567 for the FNN-MLP model and R=0.9715 for the SVM model.For the FNN-MLP model,division 3 consisted of 60%of the data for training,20%for validation,and 20%for testing.For the SVM model,division 3 had 60%of the data for training and 40%for testing.The results suggest that division 3 provides the most accurate predictions,making it the optimal choice for testing the FNN-MLP and SVM models.

Figure 4:Effect of the division of the database in term coefficient correlation(R)for the testing phase

The performance statistics of the optimal FNN-MLP model for the training,validation,testing,and overall phases in terms of R,MPE,and RMSE are presented in Table 5.The coefficient R-values of 0.9543 or higher during the training phase indicate a strong agreement between the predicted and experimental values.The testing phase’s correlation coefficient “R”measures the model’s ability to interpolate the experimental data accurately.The correlation value of 0.9567 for the testing phase demonstrates a high level of consistency between the experimental and predicted global solar radiation values(see Fig.5),indicating the effectiveness of the FNN-MLP model.

Figure 5:Comparison of predicted and experimental data in the testing phase

Table 5: Statistical evaluation of FNN-MLP model performance on solar irradiation prediction

Table 6 evaluates several support vector machine(SVM)models utilizing various kernels,including linear,polynomial,radial basis function (RBF),and sigmoid kernels.The performance of each model is evaluated using the root mean square error (RMSE) and the correlation coefficient (R) in different phases,including the training phase,testing phase,and total phase.

Table 6: Performance evaluation of SVM models with various kernels

The linear SVM model with C=10 and E=0.1 achieved an RMSE of 193.312 Wh/m2and R of 0.625 in the training phase,an RMSE of 190.394 Wh/m2and R of 0.628 in the testing phase,and an overall RMSE of 192.150 Wh/m2and R of 0.626.The polynomial SVM model with C=10,nu=1,Degree=3,and Gamma=0.125 achieved the lowest RMSE of 118.351 Wh/m2and highest R of 0.884 in the training phase,an RMSE of 120.717 Wh/m2and R of 0.873 in the testing phase,and an overall RMSE of 119.304 Wh/m2and R of 0.880.

The RBF-SVM model with C=10,nu=1,and Gamma=13.93 achieved the lowest RMSE of 32.414 Wh/m2and highest R of 0.991 in the training phase but had a higher RMSE of 57.326 Wh/m2and lower R of 0.972 in the testing phase,resulting in an overall RMSE of 38.706 Wh/m2and R of 0.988.The sigmoid SVM model with C=10,nu=0.1,and Gamma=0.125 had the highest RMSE of 229.925 Wh/m2and lowest R of 0.414 in the training phase,an RMSE of 223.432 Wh/m2and R of 0.434 in the testing phase,and an overall RMSE of 227.348 Wh/m2and R of 0.422.

Overall,the RBF-SVM model with C=10,nu=1,and Gamma=13.93 outperformed the other evaluated models,achieving the lowest RMSE and the highest R-value in the training phase,as well as the second-lowest RMSE and the second-highest R-value in the testing phase.This resulted in the lowest overall RMSE and the highest overall R-value.In comparison to the linear and sigmoid SVM models,the RBF-SVM model demonstrated a substantial improvement in both RMSE and R-values during the training and testing phases,indicating its efficacy in predicting hourly global horizontal irradiation.Furthermore,Fig.6 reveals a robust alignment between the predicted and actual values gathered from the solar irradiation measurement station,further substantiating the models’dependability and precision in predicting global irradiation using the chosen input features.

Figure 6:Comparison of predicted and actual hourly global solar irradiation

Following a comprehensive verification process,we present the performance evaluation of the FNN-MLP and RBF-SVM models for solar radiation prediction at BOU and TAM stations in Table 7.The evaluation criteria include the R,MPE,and RMSE across the training,testing,and total phases.

Table 7: Performance evaluation of developed models on solar irradiation prediction for BOU and TAM stations

Both models demonstrated reasonable accuracy in predicting solar radiation.However,the SVMRBF model exhibited superior performance compared to the FNN-MLP model in terms of R and RMSE values for both stations,particularly during the testing phase.During the training phase,the FNN-MLP model achieved R values of 0.9544 and 0.9322 for the BOU and TAM stations,respectively.In contrast,the SVM-RBF model achieved notably higher R-values of 0.9914 and 0.9204 for the same stations.

Moreover,incorporating additional inputs such as MON,PRE,and WID did not yield improvements in the model’s performance.These findings indicate that the SVM-RBF model provides a more accurate GHI prediction than the FNN-MLP model.

The performance of the present work model in predicting hourly global solar radiation was compared with different techniques used in previous studies.Table 8 compares different machine learning models used in literature studies,providing insights into the effectiveness of the present work model compared to other techniques.The RBF-SVM model used in the present work achieved the highest prediction accuracy (R=0.9876),outperforming other models such as random forest,artificial multi-neural,and adaptive approach,which also achieved high prediction accuracies ranging from R=0.95 to R=0.96.The K-means clustering-NAR model had the lowest prediction accuracy(R=0.93).

Table 8:Comparison of the present results with the literature studies in predicting hourly global solar irradiance

The comparison results suggest that machine learning models have great potential in predicting hourly global solar radiation.However,the performance of these models can vary based on various factors,such as the quality and quantity of input data,feature selection,and the specific algorithm used.Therefore,it is important to carefully consider and test different models to achieve the best results for a particular application.The high accuracy achieved by the present work using RBF-SVM indicates that it could be a useful model for future predictions in this area.

4 Conclusion

This study aims to improve the accuracy of hourly global horizental irradiation prediction by using advanced machine-learning techniques.The main objective is to develop a method that optimizes the hyperparameters of conventional machine learning models,specifically multilayer Perceptron Feedforward Neural Networks(FNN-MLP)and Support Vector Machines(SVM).To achieve this,two models were used:the FNN-MLP and the SVM.

To create the most effective model FNN-MLP,the BFGS quasi-Newton method was used as the training algorithm,four activation functions were tested in the hidden layer,and a single transfer function was used.The dimensions of the hidden layers were also varied to obtain the most accurate model possible.Regarding power and performance,the SVM model with the radial basis function(RBF) kernel function gives significantly better results than the SVM models with other functions.The RBF kernel function also shows a superior capacity in characterizing the SVM model’s hourly global solar irradiance forecast.

The statistical error difference values between the RBF-SVM model and the FNN-MLP model are significant,indicating the higher accuracy of the proposed RBF-SVM model in predicting global solar irradiance compared to the FNN-MLP model.Moreover,all the machine learning methods discussed in this study provide highly accurate predictions of global solar irradiance at different temporal resolutions.However,our results show that the RBF-SVM model performs better than the FNN-MLP-BFGS model in predicting hourly global solar irradiance,with an R-value of 0.99 and an RMSE of 38.70 Wh/m2over all phases.Moreover,this study also investigates the performance of the proposed models in different climatic regions of Algeria,which is crucial for accurately predicting solar radiation at a specific location.In this way,it could help in the design and installation of solar energy systems as well as in the evaluation of thermal conditions in building studies.

In summary,this study provides a promising alternative to the traditional methods currently used in Algeria to predict solar radiation.With its superior accuracy and performance,the RBF-SVM model can be a valuable tool for predicting global solar irradiance at any location,thus supporting the development and implementation of renewable energy sources in the country.In addition,the study opens the possibility of using these techniques in other countries with similar climate and energy needs.

Acknowledgement:The authors would like to thank the University of Relizane and the University of Dr.Yahia Fares-Laboratory Medea’s of Biomaterials and Transport Phenomena for their assistance with this work.

Funding Statement:The authors received no specific funding for this study.

Author Contributions:Study conception and design,Y.Ammi,N.Bailek,S.Hanini,L.Abualigah,E.El-kenawy;analysis and interpretation of results,A.Dahmani,Y.Ammi,S.Hanini,I.Colak,L.Abualigah;draft manuscript preparation,A.Dahmani,Y.Ammi,E.El-kenawy;writing—review and editing: N.Bailek,A.Kuriqi,N.Al-Ansari,I.Colak,L.Abualigah.All authors reviewed the results and approved the final version of the manuscript.

Availability of Data and Materials:The authors confirm that the data supporting the reported findings are all available within the article.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

主站蜘蛛池模板: 欧美另类视频一区二区三区| 亚洲成aⅴ人片在线影院八| 亚洲国产日韩欧美在线| 欧美成人午夜影院| 免费中文字幕一级毛片| 日韩欧美在线观看| 成人免费午夜视频| 色婷婷啪啪| 亚洲高清国产拍精品26u| 欧美国产中文| 国产素人在线| 亚洲精品麻豆| 亚洲av无码人妻| 久久免费成人| 国产高潮流白浆视频| 高清视频一区| 国产一区二区三区夜色| 亚洲系列中文字幕一区二区| 五月天福利视频 | 亚洲一区二区三区国产精品 | 91在线精品免费免费播放| 制服无码网站| 一区二区三区高清视频国产女人| 亚洲性日韩精品一区二区| 国产一二三区在线| 国产AV无码专区亚洲精品网站| 欧美成人午夜视频| 久久一色本道亚洲| 成人精品免费视频| 久久久久久久97| 精品日韩亚洲欧美高清a| 久久综合婷婷| 女人爽到高潮免费视频大全| 欧洲极品无码一区二区三区| 国产丰满大乳无码免费播放 | 2021天堂在线亚洲精品专区| 91无码人妻精品一区| 男女猛烈无遮挡午夜视频| 午夜在线不卡| 91九色国产在线| 国产高潮流白浆视频| 呦视频在线一区二区三区| 亚洲第一福利视频导航| 波多野结衣的av一区二区三区| 思思热精品在线8| 青青国产成人免费精品视频| 视频一区视频二区中文精品| 国产美女无遮挡免费视频| 全免费a级毛片免费看不卡| 亚洲热线99精品视频| 午夜激情婷婷| 国产91熟女高潮一区二区| 国产精品自在线拍国产电影| 成人久久18免费网站| 美女裸体18禁网站| 极品性荡少妇一区二区色欲| 玩两个丰满老熟女久久网| 麻豆精品在线| 中文字幕第4页| 亚洲一区网站| 秋霞午夜国产精品成人片| 制服丝袜在线视频香蕉| 色婷婷成人| 国产美女无遮挡免费视频网站| 国产精品尹人在线观看| 福利在线免费视频| 香蕉eeww99国产精选播放| 四虎综合网| 毛片视频网| 专干老肥熟女视频网站| 一区二区三区国产精品视频| 国产精品2| 亚洲无限乱码| 欧美日本在线一区二区三区| 亚洲第一色网站| 动漫精品中文字幕无码| 在线视频一区二区三区不卡| 亚洲国产第一区二区香蕉| 99久久婷婷国产综合精| 狠狠ⅴ日韩v欧美v天堂| 国产办公室秘书无码精品| 久久精品国产91久久综合麻豆自制|