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

    A Hybrid Neural Network Model for Marine Dissolved Oxygen Concentrations Time-Series Forecasting Based on Multi-Factor Analysis and a Multi-Model Ensemble

    2021-04-22 11:35:02HuiLiuRuiYangZhuDuanHaipingWu
    Engineering 2021年12期

    Hui Liu*, Rui Yang, Zhu Duan, Haiping Wu

    Institute of Artificial Intelligence and Robotics(IAIR),Key Laboratory of Traffic Safety on Track of Ministry of Education,School of Traffic and Transportation Engineering,Central South University, Changsha 410075, China

    Keywords:Dissolved oxygen concentrations forecasting Time-series multi-step forecasting Multi-factor analysis Empirical wavelet transform decomposition Multi-model optimization ensemble

    ABSTRACT Dissolved oxygen (DO) is an important indicator of aquaculture, and its accurate forecasting can effectively improve the quality of aquatic products. In this paper, a new DO hybrid forecasting model is proposed that includes three stages: multi-factor analysis, adaptive decomposition, and an optimizationbased ensemble. First, considering the complex factors affecting DO, the grey relational (GR) degree method is used to screen out the environmental factors most closely related to DO. The consideration of multiple factors makes model fusion more effective.Second,the series of DO,water temperature,salinity,and oxygen saturation are decomposed adaptively into sub-series by means of the empirical wavelet transform (EWT) method. Then, five benchmark models are utilized to forecast the sub-series of EWT decomposition. The ensemble weights of these five sub-forecasting models are calculated by particle swarm optimization and gravitational search algorithm(PSOGSA).Finally,a multi-factor ensemble model for DO is obtained by weighted allocation. The performance of the proposed model is verified by timeseries data collected by the pacific islands ocean observing system (PacIOOS) from the WQB04 station at Hilo. The evaluation indicators involved in the experiment include the Nash–Sutcliffe efficiency(NSE), Kling–Gupta efficiency (KGE), mean absolute percent error (MAPE), standard deviation of error(SDE),and coefficient of determination(R2).Example analysis demonstrates that:①The proposed model can obtain excellent DO forecasting results;②the proposed model is superior to other comparison models; and ③the forecasting model can be used to analyze the trend of DO and enable managers to make better management decisions.

    1. Introduction

    Dissolved oxygen (DO) plays an important role in the sustainable development of aquaculture and marine ecology in coastal cities.Nowadays,with the rapid development of industry,large volumes of greenhouse gas emissions are increasingly intensifying global warming[1].As a result,the seawater temperature is rising and DO is escaping from seawater [2]. If the DO concentrations in seawater are too low, marine organisms such as fish and shrimp will be deprived of oxygen,leading to slow growth and even suffocation[3].For offshore aquaculture farmers,the massive scale of the death of aquatic products will bring huge economic losses [4]. More seriously,it is difficult for DO to diffuse into the depths of the ocean,and a lack of oxygen will force aquatic organisms to change their habitats[5]. The marine ecosystem could change significantly as a result.Once the virtuous cycle of the ecosystem is broken,the survival of marine organisms will be greatly threatened [6]. Accurate DO content forecasting is very necessary for marine ecological balance and for the aquaculture development of coastal cities. It can provide a trend analysis of future changes in DO and facilitate the management of decision-makers by enabling them to take effective measures in advance[7,8].Influenced by factors such as nutrients,climate, ecological environment, and the life activities of aquatic organisms, the DO content in seawater is nonlinear and has a considerable time delay. This characteristic makes DO difficult to predict.In addition,it is easy for data obtained by sensors to be lost during transmission, and such data contains uncertain outliers,posing more challenges to the accuracy of forecasting[9].

    1.1. Related works

    In recent years, the rise of artificial intelligence(AI)technology has spawned a wave of AI models. On this basis, DO forecasting models can be divided into two types: physics-based models and data-driven models. Data-driven models include AI-based models and statistics-based models. To comprehensively evaluate stateof-the-art efforts in DO models, we summarize the literature on statistics-based models and AI-based models.

    Statistics-based models have the advantage of easy calculation in application. In addition, a statistics-based model can predict an unknown situation by mining the potential relationship in the DO series. In practical applications, DO series are nonlinear and time-delayed. Faruk [10] used the autoregressive integrated moving average (ARIMA) model for DO forecasting to deal with the instability and nonlinearity of time series. In the experiments of Li et al. [11], the grey model was proposed to predict the trend term of DO.Huan et al.[12]obtained the optimal parameter values of the model through a Bayesian evidence framework. Khan et al.[13] compared the new Bayesian regression model with the new autoregressive modified fuzzy linear regression method.The comparison results show that the method based on fuzzy numbers can better capture DO changes in an urban river environment. Khan et al. [14] predicted changes in DO concentration by constructing multiple linear regression(MLR) models, in which the uncertainty of DO change was effectively characterized and propagated.In later research by Kisi and Parmar [15], the MLR model was used as an experimental comparison to explore the performance of other prediction models. In general, these statistics-based models have obtained the ideal forecasting effect for DO. However, forecasting accuracy remains to be improved.

    Compared with statistics-based models, AI-component models have much better forecasting performance.The core of this section focuses on describing models with AI components.Commonly used AI algorithms include the multi-layer perceptron(MLP)[16],radial basis neural network (RBNN) [17], back propagation neural network (BPNN), extreme learning machine (ELM) [18], least square support vector machine (LSSVM) [19], and fuzzy neural network(FNN) [20,21], among others [22]. Since the selection of different parameters has a huge impact on the performance of an AI model,it is necessary to determine the optimal parameters by means of an optimization algorithm.Ren et al.[23]used a genetic algorithm(GA)to optimize the parameters of FNN.In addition,particle swarm optimization (PSO) [24], cauchy particle swarm optimization(CPSO)[25],the firefly algorithm(FFA)[26],and other optimization algorithms are used for parameter selection for intelligent models[27]. With its rapid development, deep learning has been widely studied for its excellent adaptive ability. Ma et al. [28] solved the problem of a sparse matrix by deep matrix factorization. Ren et al. [29] used a deep brief network (DBN)’s powerful feature extraction and functional representation capabilities to deal with highly complex nonlinear DO time-series data.Nevertheless,given the complex laws of DO changes and many other factors, simple machine learning basic predictors and parameter optimization methods are insufficient to allow accurate DO prediction.To further improve the forecasting performance of AI models, an increasing number of scholars are using feature selection, decomposition,and ensemble methods for the state-of-the-art research of DO forecasting.Further details about such research include the following:

    (1) The change in DO concentration in the ocean and in freshwater environments is not isolated,as it is affected by many factors including temperature,turbidity,pH,chlorophyll,and specific conductance [30]. Changes in these factors have seemingly invisible effects on DO concentration. To effectively use these factors, the specific effects need to be studied. Shi et al. [22] used a clustering method to segment water quality time series,and finally improved the forecasting accuracy of DO.Ren et al.[23]arranged the factors that are positively and negatively correlated with DO through correlation analysis.It was found that the model showed a better forecasting performance when factors strongly related to DO were entered into the model [31]. From the perspective of the generation and consumption of DO, these influencing factors should be organically incorporated rather than directly and crudely input into the prediction model. This is not a comprehensive consideration. The direct superposition of a multi-factor forecasting series may prove to reduce the DO forecasting accuracy of a model. The weight allocation of each influencing factor needs to be analyzed,and the optimal combination method needs to be identified [32].Therefore, exploring methods involving the multi-factor analysis of DO forecasting has great research value [33].

    (2) The nonstationary properties of the original DO series are not favorable for forecasting. The extreme volatility of the DO series makes it difficult to predict. The birth of the decomposition method largely solved this problem. A decomposed sub-series has independent oscillation components and is more stable, making it easier to predict [34,35]. To apply this theoretical effect to practical problems,many data preprocessing decomposition methods for AI models have been proposed.Widely used decomposition algorithms include ensemble empirical mode decomposition(EEMD) [12], discrete wavelet transform (DWT), and variational mode decomposition (VMD) [36]. DO series are more predictable after decomposition. However, the current mainstream decomposition methods require users to determine the number of decomposition layers based on experience. This practice inevitably adds human error [37]. To improve the forecasting model performance,a decomposition algorithm that can adaptively determine the decomposition level should be investigated.

    (3) After years of research and development, some ensemble models can be used for the intelligent forecasting of DO.An ensemble model can combine the advantages of multiple forecasting models to achieve better forecasting performance [35]. By means of a layered ensemble,Zhu et al.[38]successfully solved the problem of it being impossible to simultaneously achieve the accuracy of high and low concentration forecasting.Iterated stepwise multiple linear regression (ISMLR) is used to preprocess the data. An artificial neural network (ANN) and MLP are used to run layered forecasting.Finally,compromise programming(CP)is used to evaluate and select the best results. Kisi et al. [39] proposed a new bayesian model averaging (BMA) method as an ensemble model.Through comparison experiments with ELM, ANNs, an adaptive neuro-fuzzy inference system (ANFIS), a classification and regression tree (CART), and MLR, it was found that the BMA ensemble method is very effective for DO prediction. The essence of an ensemble model is to improve the forecasting accuracy and scope of model application by making use of the complementary advantages of multiple models. Unfortunately, the ensemble methods described above only use multiple models for hierarchical prediction, and then simply superimpose the results. A direct ensemble of multiple models will greatly increase the complexity of the ensemble model, which can easily lead to overfitting. At the same time, the shortcomings of a single model cannot be compensated for by other models.Therefore,it is very important to select a reasonable ensemble of models.Different machine learning prediction models have individual characteristics and can cope with time series in different states. If they can be combined organically, the ensemble model will have complementary advantages. Appropriate weight allocation can cover up the prediction defects of benchmark models and enhance the advantages of benchmark models.Therefore,it is necessary to explore the scientific ensemble method of multiple models.

    The DO forecasting models discussed above have been demonstrated to be effective. Still, there are gaps in scientific research that need to be filled in order to further improve DO forecasting performance. Table 1 [10–19,21–30,36,38–40] summarizes the reviewed state-of-the-art works using data-driven methods in DO forecasting.

    1.2. The novelty of this study

    To summarize the above literature,it can be seen that there are few comparative studies on multi-factor analysis,adaptive decomposition analysis, and scientific ensemble analysis in DO forecasting research. In view of the above limitations, a novel forecasting model is herein proposed for DO forecasting: the MF-RNNs-EWTBEGOE model. MF represents the proposed multi-factor analysis method. RNNs represents the replicator neural networks outlier detection method, EWT stands for empirical wavelet transform,BEGOE stands for BFGS-ENN-GRNN-ORELM-ELM. This model can be divided into the three stages of multi-factor analysis, adaptive decomposition, and optimization-based ensemble. Each sub-layer under the four factors uses the same type of benchmark model to predict;finally,the forecasting results of the model are obtained through a combination of the results.In the optimization ensemble phase, the ensemble weights for the Broyden–Fletcher–Goldfarb–Shanno(BFGS)model,elman neural network(ENN)model,general regression neural network (GRNN) model, outlier-robust extreme learning machine (ORELM) model, and ELM model are optimized by particle swarm optimization and gravitational search algorithm(PSOGSA) to obtain the BEGOE model. Section 2 provides the details of the proposed model.To fully understand the role of each part of the model, we performed numerous comparative experiments.

    Table 1 A summary of the reviewed state-of-the-art data-driven DO forecasting methods.

    The innovations and contributions of this study are described below.

    (1) A hybrid three-stage ocean DO forecasting ensemble model is designed. The model considers multi-factor analysis, adaptive decomposition analysis, and scientific ensemble analysis, thereby filling three of the research gaps described in Section 1.1. The multi-factor analysis reasonably considers the effects of temperature, salinity, turbidity, chlorophyll, and oxygen saturation on the concentration of DO in the ocean. EWT decomposition is utilized to adaptively improve the predictability of time series. The nonstationarity of DO series is reduced. The PSOGSA optimization method is used to optimize the weights of multiple benchmark models. The complementation of the benchmark model further improves the performance of the hybrid model. These three modules work together to achieve accurate forecasting of DO.

    (2) Several factors that affect the concentration of DO in the ocean are reasonably considered.The multi-factor analysis method not only considers the influence of single and multiple factors on model accuracy, but also considers correlations among multiple factors.In this way,a hybrid model with better fusion performance can be obtained, avoiding the unilateral influence and local disadvantage of a single factor.This method also makes up for the shortcomings of direct multi-factor input, as such input may have a negative effect. This approach can take into account the reasons for a change in DO concentration from various aspects.

    (3)The adaptive data preprocessing method has achieved ideal results in regard to the decomposition effect. Compared with the mainstream decomposition method, adaptive decomposition avoids the error caused by the artificial selection of decomposition layers.In this way,the decomposition layer can be predicted more easily. In addition, the EWT decomposition method discards the residual signal of the decomposition mode to extract meaningful information from the ocean DO time series. The EWT method also compensates for the sensitivity of some decomposition algorithms to noise and sampling to a certain extent. The support of mathematical theory allows the machine learning algorithm performance to develop effectively.

    (4)The organic combination of multiple benchmark models fills the research gap for the ensemble model. A meta-heuristic optimization algorithm is used to assign weights to BFGS,ENN,GRNN,ORELM, and ELM. The ensemble approach can achieve the advantages of multiple benchmark models while reducing the adverse effects of a single model. The ensemble model solves the poor robustness problem of some forecasting models and fully demonstrates the superiority of the complementary performance. Therefore,a multi-model complementary ensemble framework based on meta-heuristic optimization algorithm has great value.

    2. The proposed multi-factor forecasting model

    As mentioned earlier,the modeling process of the proposed MFRNNs-EWT-BEGOE model can be divided into three stages: multifactor feature extraction, EWT decomposition, and multi-model optimization ensemble optimization ensemble. Fig. 1 shows the modeling flow. It is worth mentioning that the proposed model uses the multi-input multi-output strategy (MIMOS), which is commonly used in the field of time-series forecasting [41]. Compared with recursive strategy(RS),the MIMOS has a smaller cumulative error[42].Fig.2 shows the schematic diagram of the RS and MIMOS.

    2.1. Stage 1: The multi-factor analysis method

    2.1.1. Substage 1.1: Replicator neural networks outlier detection

    Outlier rejection is an important process for DO data. Preprocessing makes the data look more orderly and makes it more easily learned by the forecasting model, without paying attention to the sources of the outliers.It is worth noting that only the training sets are corrected, while the test sets remain the same. This is a more convincing way to verify the validity of the model. RNNs is used to detect outliers of the original ocean DO series. This is a multilayer feed-forward neural network [43].

    Fig.1. The modeling flow of the proposed three-stage model.GR:grey relational;GWO:grey wolf optimizer;BA:bat algorithm;MVO:multi-verse optimization algorithm;WOA: whale optimization algorithm; NBA: new bat algorithm.

    Fig. 2. Schematic diagram of multi-step forecasting strategies.

    Assuming that the output of the ith neuron in the k layer of the RNNs is Sk(Iki), Ikirepresents the input of the ith neuron in the k layer,and Skrepresents the activation function used in the k layer.The input of the neuron can be expressed as follows:

    where N represents the number of steps and a3represents the rate of promotion to the next level.

    In general,continuous data entering the third layer is converted into a batch of discrete values through the S3stepwise activation function of the RNNs.This is equivalent to mapping the series sample to N clusters.Finally,the RNNs can calculate a single outlier and a small cluster of outliers. This completes the outlier detection of marine data at the WQB04 site. Algorithm S1 in the Appendix A shows the pseudo-code of the RNNs outlier detection method.

    2.1.2. Substage 1.2: Multi-factor analysis of the grey relational (GR)degree method

    To screen out the most favorable factors for DO forecasting,the GR degree method is used to select the features of the original multi-factor series. In this example, the GR degree method is used to calculate the correlation of environmental factors such as temperature, salinity, turbidity, chlorophyll, and oxygen saturation.The greater the degree of grey correlation between time series,the closer the distance between the two series. Feature selection

    2.2. Stage 2: The EWT-decomposition method

    2.3. Stage 3: The multi-model optimization ensemble method

    To further improve the forecasting performance,the BFGS,ENN,GRNN,ORELM,and ELM are combined to obtain the BEGOE model.The ensemble weights are optimized by the PSOGSA. Fig. 4 shows the weighted optimization scheme. The optimization objective function of PSOGSA is minimizing the mean absolute percent error(MAPE).The MAPE refers to the average absolute percentage error between different samples’ actual and forecasting values. The objective function is calculated in the validation data. The specific objective function is shown as follows:

    Fig. 3. EWT-decomposed and layer-by-layer forecasting.

    Fig. 4. The PSOGSA optimization ensemble process under multi-factor analysis.

    Considering that the DO, water temperature, salinity, and oxygen saturation series are involved in this example, the input of the BEGOE model are four series, and the output is a DO series.By minimizing the objective function MAPE,appropriate ensemble weights can be obtained.The effects of water temperature,salinity,and oxygen saturation are added to the DO series.Finally,the forecasting results of the BEGOE can be obtained. Algorithm S3 of the Appendix A shows the pseudo-code of the multi-model optimization ensemble method.

    3. Experiment and analysis

    In this section, multiple comparative experiments are carried out to verify the effectiveness of the proposed model. In Section 3.1,six real-time series sets collected by the ocean observing system are introduced. In Section 3.2, the performance indices for evaluating the forecasting results are demonstrated. In Section 3.3, the implementation details and the contributions of each stage are explained. In Section 3.4, two experiments are set up.Experiment I compares the contributions of the model components horizontally, and Experiment II compares several state-of-the-art models to demonstrate the superiority of the proposed model.

    3.1. Data description

    In this section, time-series data such as DO, temperature,salinity, turbidity, chlorophyll, and oxygen saturation collected by the Pacific Islands Ocean Observing System from the WQB04 station at Hilo are used for experimental research. The time interval of the six time-series data is 15 min. To effectively build the subsequent models, the dataset is divided into three parts: a training set, validation set, and test set. The data corresponding to the six factors were all selected from the time period of 1–30 December 2016. Among them, the data for 16–30 December is regarded as Dataset #1, while the data for 1–15 December is regarded as Dataset #2. The total length of each series corresponding to each factor is 1440. The training set contains the 1st–1248th data, the validation set contains the 1249th–1344th data, and the test set contains the 1345th–1440th data. Fig. S1 and Table S1 in the Appendix A present the dataset information after RNN outlier processing.

    3.2. Performance evaluation

    An index set is applied to evaluate the deterministic forecasting performance,including the coefficient of determination(R2),Nash–Sutcliffe Efficiency (NSE) [47] index, Kling–Gupta fficiency (KGE)[48],standard deviation of error(SDE),and MAPE.The expressions of these five indices are shown as follows:

    3.3. Modeling analysis

    3.3.1. Multi-factor analysis

    The DO in the ocean is affected by many factors.In this paper,a variety of factors affecting the increase and decrease of DO concentrations are analyzed. Fig. 5 shows the system dynamics model.

    As shown in Fig.5,many factors affect DO concentrations.Usually, atmospheric diffusion, water exchange, mechanical aerobics,aquatic photosynthesis, and the like are the main sources of DO.Under natural conditions, the DO content in the ocean is closely related to the partial pressure of oxygen in the air and the temperature of the water.When the rate of oxygen dissolved into the water is equal to the rate of oxygen escaping from the water,the dissolution reaches a dynamic equilibrium. The equilibrium of DO in the ocean is also affected by salinity, turbidity, pH, and other factors.In addition, chemical oxidation reactions, organic decomposition,breathing by phytoplankton at night, and aquatic respiration are the main factors of marine DO consumption. In short, there are many factors affecting changes in DO content, and the law of distribution variation is difficult to explain by a simple mechanism.

    In this experiment, the relationship between temperature,salinity, turbidity, chlorophyll, oxygen saturation, and DO concentrations is calculated by using the GR degree method. Table 2 shows the specific results.According to Table 2,temperature,salinity,and oxygen saturation have the highest correlation degree with DO. The correlation indicators reach 0.61, 0.68, and 0.87, respectively,and all exceed the set threshold of 0.5. It is worth mentioning that salinity contains nutrients, and nutrients play an important role in DO concentration.Considering that high concentrations of turbidity will influence the photosynthesis of plants in water and affect oxygen production, turbidity was also input into the model as a DO-related factor. The experiment adopts the form of multi-factor variable input and single-factor variable output. In an actual comparison of repeated experiments, it was found that the forecasting accuracy is the best when temperature, salinity,and oxygen saturation are input into the model as the influential factors of DO. The experimental results show that if lowcorrelation factors such as turbidity,chlorophyll, and pH are input into the model, the forecasting accuracy decreases accordingly. In this case, the index of chlorophyll always maintains a value of-8.8 × 10-7. Therefore, it does not show a correlation with the change in DO, so the correlation coefficient between chlorophyll and DO is not a number.

    3.3.2. Decomposition model analysis

    In this section,the EWT is used to reduce the non-stationarity of the DO series,the water temperature series,the salinity series,and the oxygen saturation series, all of which are adaptively decomposed. It is worth mentioning that only the training set data is decomposed.After the original DO series and the highly correlated series are decomposed, the sub-series having more stable characteristics are easier to predict.The primary predictor simply concentrates on a certain frequency band within these sub-series.Eventually,the final result after data preprocessing can be obtained by superimposing the prediction result of the sub-predictor. This section compares the performance changes of the forecasting models by decomposing the input series.

    Fig. 5. System dynamics model diagram of dissolved oxygen change.

    Table 2 Correlation coefficient between environmental factors and dissolved oxygen in the WQB04 dataset.

    3.3.3. Optimization method analysis

    In this section, the PSOGSA is compared with multiple alternative methods, including the grey wolf optimizer (GWO), bat algorithm (BA), multi-verse optimization algorithm (MVO), PSO,whale optimization algorithm (WOA), and new bat algorithm(NBA).Most of these optimization algorithms are based on heuristic bionic algorithms, and it is difficult to clearly explain which is more suitable from a theoretical point of view. Therefore, experiments are needed to verify the DO prediction effect of each optimization algorithm in this situation.

    3.4. Comparison study

    3.4.1. Experiment I: Comparisons with each component of the proposed model

    In Experiment I,five baseline models based on the components of the proposed model are constructed and used to evaluate the proposed model. Table 3 lists their specific expressions. To compare the effects of environmental parameters related to DO in actual calculations, we established a multi-factor analysis experiment. Models 1–3 are trained without multi-factor analysis. In the same situation, Models 4 and 5 and the proposed model are added with multi-factor analysis. In other words, Models 4 and 5 and the proposed model are trained with environmental factorsincluding DO series,temperature series,salinity series,and oxygen saturation series. To compare the effects of the optimization ensemble, a multi-model ensemble comparison experiment is established. Compared with Models 1 and 4, Models 2 and 5 are optimized and combined with five benchmark models. To verify the superiority of the EWT decomposition method,a set of experimental comparisons with or without data preprocessing are also implemented. More specifically, Model 3 and the proposed model add a module for EWT data preprocessing to reduce the instability of the original series. Tables 4 and 5 list the DO testing dataset forecasting results. In addition, in order to reflect the impact of optimization algorithms, the seven optimization algorithms described in Section 3.3.3 are used in the multi-model optimization ensemble module. Fig. 6 shows testing dataset forecasting results in 3-step of the PSOGSA and other optimization algorithms.

    Table 3 Experimental setup of baseline models and the proposed model.

    The conclusion can be summarized from Tables 4 and 5, and Fig. 6 as follows:

    (1) Model 2 (RNNs-BEGOE) is significantly better than Model 1(RNNs-BFGS/ENN/GRNN/ORELM/ELM) regardless of whether in Dataset #1 or Dataset #2, which shows that the multi-model ensemble method is effective. This method combines the advantages of five benchmark models.The multiple iterations of the optimization algorithm complete the combination of each model through the optimal weight allocation.The ensemble model incorporates the advantages of each benchmark model while discarding their negative performance.Taking the 1-step forecasting results of Dataset #1 as an example, the MAPE of RNNs-BFGS, RNNs-ENN,RNNs-GRNN, RNNs-ORELM, RNNs-ELM, and RNNs-BEGOE are 1.76%, 1.51%, 2.74%, 1.49%, 1.61%, and 1.32%, respectively. Other evaluation indicators such as NSE, KGE, SDE, and R2all show that Model 2 is better than Model 1 in forecasting performance.

    (2) Model 3 (RNNs-EWT-BEGOE) is significantly better than Model 2 (RNNs-BEGOE), which shows that the decompositionbased data preprocessing method can have a positive effect on improving model prediction performance. The EWT adaptively decomposes the original series into multiple sub-series. The decomposed series has independent oscillation components and strong stability. The EWT has strong mathematical support whileeliminating the interaction between different sub-series.The combination of these advantages makes EWT an excellent candidate for data preprocessing.Taking the 1-step forecasting results of Dataset#1 as an example, the NSE, KGE, MAPE, SDE, and R2are 0.70, 0.75,1.32%, 0.10 mg·L-1, and 0.83 for Model 2, and 0.97, 0.96, 0.38%,0.04 mg·L-1, and 0.97 for Model 3, respectively.

    Table 4 Evaluation indices of the compared models in Experiment I (Dataset #1).

    (3) Models 4 (MF-RNNs-BFGS/ENN/GRNN/ORELM/ELM) and 5(MF-RNNs-BEGOE) and the proposed model with multi-factor analysis performed significantly better in comparison with the corresponding Models 1–3 without multi-factor analysis. This work shows that environmental factors such as temperature, salinity,and oxygen saturation have an impact on DO concentration.There is a complex nonlinear relationship between them,and the consideration of a single DO factor is not comprehensive. Taking the 1-step forecasting results of Dataset #1 as an example, the NSE,KGE, MAPE, SDE, and R2are 0.70, 0.75, 1.32%, 0.10 mg·L-1, and 0.83 for Model 2, and 0.91, 0.90, 0.67%, 0.07 mg·L-1, and 0.91 for Model 5, respectively. The experimental results show that the multi-factor consideration is correct. When multiple factors are reasonably input into the forecasting model, a more accurate and more scientific forecasting result can be obtained.

    (4)Considering all the experimental datasets,forecasting steps,and different model evaluation indicators, the proposed hybrid model has the most outstanding forecasting performance. The excellent robustness and accuracy of the proposed model make its prediction result very close to the actual DO value. Taking the 1-step forecasting results of Dataset #1 as an example, the NSE,KGE, MAPE, SDE, and R2of the proposed model are 1.00, 1.00,0.11%, 0.01 mg·L-1, and 1.00, respectively. Taking the 1-step forecasting results of Dataset #2 as an example, the NSE, KGE, MAPE,SDE, and R2of the proposed model are 1.00, 0.98, 0.18%,0.01 mg·L-1,and 1.00,respectively.Both sets of experiments show that the proposed model has good durability.

    In summary, multi-factor analysis, adaptive decomposition analysis,and scientific ensemble model analysis are indispensable for improving the performance of the proposed hybrid model. The framework of the hybrid model has strong interpretability, and it provides a direction for the improvement of DO forecasting.

    3.4.2. Experiment II: Comparison with existing models

    In recent years, many AI-based models have been proposed for DO prediction.In Experiment II,we reproduced three state-of-theart models that were published in 2018–2020 to compare them with our proposed hybrid model in order to verify the DO forecasting performance. These models are relatively complex, and basic machine learning models can also achieve satisfactory results. To complete a reasonable and scientific model evaluation, we added the support vector machines (SVM) model based on ANN, and the DBN model based on deep learning for comparison.The specific model introduction is elaborated as follows:

    Table 5 Evaluation indices of the compared models in Experiment I (Dataset #2).

    Li et al. [31] used a recurrent neural network, long short-term memory network (LSTM), and gate recurrent unit (GRU) to construct three DO forecasting models,respectively,in order to determine the most suitable one. In the experiment, the correlation coefficients between pH, turbidity, temperature, NH, and DO are calculated,and then these series are uniformly input into the forecasting model. Compared with Li’s direct inputting of all parameters, we set thresholds and experimental screening methods to finally determine the factors that can have a positive effect on DO forecasting. Li’s approach risks reducing accuracy. In fact, the correlation coefficient of the impact factors can only be used as a reference indicator, and the true relationship between these environmental factors and DO is very complicated.

    Huan et al.[12]combined the EEMD data preprocessing method and LSSVM to improve the model’s performance for DO forecasting.EEMD is used to decompose the DO series,and then LSSVM is used to predict these intrinsic mode function (IMF) components separately. The synthesized prediction result is obtained by superposition.Before the final DO forecasting results are generated,BPNN is used to reconstruct the forecast results, thereby eliminating the wrong analysis. It is a pity that EEMD lacks sufficient theoretical support in terms of mathematical definition, which results in decomposition flaws. We use a novel signal processing tool EWT to overcome this shortcoming. At the same time, the adaptive decomposition characteristics of EWT can determine the optimal number of decomposition layers to maximize the decomposition effect.

    Ren et al. [23]used GA to optimize the center and width of the center layer of the FNN to determine the best combination and improve DO prediction performance. In the proposed model, we apply the optimization algorithm to the ensemble of multiple benchmark models. The advantage of this approach is to reasonably combine the superior performance of multiple AI models while discarding bad performance. Compared with a single parameter-optimized predictor, the complementarity of multiple models can make the model more mature and allow it to obtain more adequate training. Among the models used, the reasonable weight allocation obtained by multiple iterations plays a key role.

    More specifically,Tables 6 and 7,and Figs.7 and 8 show the DO testing dataset forecasting results of the proposed model and stateof-the-art models.

    The conclusion can be summarized from Figs. 7 and 8, and Tables 6 and 7, as follows:

    Fig.6. Testing dataset forecasting results in 3-step of the PSOGSA and other optimization algorithms(the above-mentioned different optimization methods are all based on the multi-factor-EWT-BEGOE ensemble model to conduct control variable comparison experiments)(Dataset#1).Time(15 min)indicates that the sampling interval of data is 15 minutes.

    Table 6 The DO testing dataset forecasting results of the proposed model and state-of-the-art models (Dataset #1).

    Considering all datasets, forecasting steps, and different timeseries evaluation indicators,the proposed model has the best forecasting performance in comparison with the other state-of-the-art models. Taking the 1-step forecasting results of Dataset #1 as an example, the MAPE of SVM, DBN, Li’s model, Huan’s model, Ren’s model, and the proposed model are 1.41%, 1.44%, 0.52%, 0.56%,0.51%, and 0.11%, respectively. The proposed model performs well in forecasting robustness and can still perform satisfactorily even in higher prediction steps. Taking the 3-step forecasting results of Dataset #1 as an example, the MAPE of SVM, DBN, Li’s model,Huan’s model, Ren’s model, and the proposed model are 1.99%,2.06%, 1.05%,0.89%, 1.21%,and 0.17%, respectively. Based on other persistence indices,the performance of the proposed model is very outstanding.Taking the 3-step forecasting results of Dataset#1 as an example, the NSE, KGE, MAPE, SDE, and R2of the proposed model are 0.98, 0.98, 0.17%, 0.04 mg·L-1, and 0.99, respectively.Taking the 3-step forecasting results of Dataset #2 as an example,the NSE, KGE, MAPE, SDE, and R2of the proposed model are 0.96,0.94, 0.58%, 0.04 mg·L-1, and 0.98, respectively.

    To sum up, the forecasting results of different datasets show that the proposed hybrid forecasting model has a reasonable framework and achieves a satisfactory forecasting performance.This is mainly because the proposed model hybridizes multifactor analysis,adaptive decomposition analysis,and the scientific ensemble model analysis method.This combination is very helpful in improving the forecasting accuracy and robustness of the proposed model.

    Table 7 DO testing dataset forecasting results of the proposed model and state-of-the-art models (Dataset #2).

    Fig. 7. DO testing dataset forecasting results of the proposed model (MF-RNNs-EWT-BEGOE) and other models (Dataset #1).

    4. Application potential

    The level of DO concentration in a marine environment is caused by two effects: the oxygen consumption effect, which decreases DO concentration;and the oxygen recovery effect,which increases DO concentration.The life activities of aquatic organisms consume a certain amount of oxygen, thereby reducing the DO concentration. The photosynthesis of aquatic plants and the diffusion of oxygen in the air will increase the DO concentration.Fig. 9(a) shows the changes of actual DO time series data and 1-step, 2-step, and 3-step forecasting results. The following facts are illustrated in Fig. 9: ①As the temperature increases, the DO concentration decreases, and as the temperature decreases, the DO concentration increases. This is because higher temperatures will force DO in the water to escape to the surface (Fig. 9(b)).② Salinity and DO are negatively correlated, due to oxygen consumption. Sufficient nutrients will promote the life activities of aquatic organisms, resulting in increased breathing and oxygen consumption. An increase in salinity will also promote DO to a level closer to saturation (Fig. 9(c)). ③High turbidity will reduce the photosynthesis of aquatic plants and affect the production of DO, so turbidity has a certain inhibitory relationship with the increase of DO(Fig.9(d)).④The oxygen content dissolved in water is also known as the oxygen saturation. When the oxygen exchange between water and the atmosphere is in equilibrium,the concentration of DO in the water is most similar to the fluctuation of oxygen saturation(Fig.9(e)).Mathematical modeling combined with AI technology can make a relatively scientific prediction of the changes in DO concentration. This prediction can provide useful guidance to people in related industries.

    Fig. 8. Scatter plots of actual observations and corresponding forecasting results generated by the proposed model and other models.

    This study proposes a novel hybrid model that can generate reliable DO forecasting results. The model takes temperature,salinity, turbidity, oxygen saturation, and many other factors closely related to DO as input variables,and finally outputs a DO forecasting result.Fig.9(f)shows the division of DO data and influence factors data. The specific multiple-input–single-output forecasting results are shown in Fig. 9. Potential applications for users of the proposed model are presented as follows:

    · Aquaculture farmers can utilize the proposed DO forecasting model to determine future DO concentration changes in ponds and offshore aquaculture areas.They can then take corresponding emergency measures under guidance to reduce their economic losses.

    · Fishery department managers can fully understand the situation of and trends in DO by means of accurate forecasting.This can provide a scientific decision-making basis for improving the living environment of marine life, which is conducive to the restoration and regulation of ecological water quality.

    5. Conclusions

    The development of DO forecasting technology has brought more possibilities to seafood farmers and related industries.Accurate real-time forecasting of DO can provide managers with scientific guidance for decision-making.Unfortunately,the instability of the marine DO series and various other potential factors present challenges to DO forecasting. In our study, a novel hybrid model is proposed for DO forecasting: the MF-RNNs-EWT-BEGOE. The main characteristics of this model include multi-factor analysis,adaptive decomposition, and a multi-model ensemble. After the analysis of the experimental results in the case study,the following conclusions can be drawn:

    (1)Considering that many factors affect DO concentration in the ocean, a multi-factor analysis method is proposed to incorporate the effects of these factors on DO concentration. The influence of abundant physical, chemical, and biological systems makes the DO fluctuation mechanism more complicated. Reasonable multifactor consideration will make model forecasting more accurate.

    (2) An adaptive data preprocessing decomposition method is used to reduce the instability of marine DO series, making subseries easier to predict.The experimental results demonstrate that this method of forecasting after decomposing the original series is effective.

    (3)A novel multi-model ensemble method is proposed to complete the combination of multiple models, which can absorb the excellent quality of each benchmark model to make the ensemble model more mature. A comparison with the experimental results of multiple benchmark models shows that the complementary advantages of multiple models can indeed cause the performance of the ensemble model to develop in a better direction.

    (4) Each component of the proposed hybrid model can significantly improve the model’s forecasting performance for DO. By detecting and eliminating outliers through RNNs, the negative impact of extreme outliers on the model can be reduced.The environmental factors that affect the change of DO can be considered more comprehensively through the multi-factor correlation analysis by GR. EWT adaptive decomposition effectively improves the forecasting accuracy of the model and reduces the impact of series non-stationarity.The BEGOE ensemble model combines the advantages of multiple benchmark models to maximize the model’s forecasting performance. Thus, the proposed hybrid model is significantly better than the other models used for comparison.

    Fig. 9. DO vs input data change trend chart of the multiple–input–single–output model.

    Given the potential for application described in Section 4, it is clear that accurate DO prediction can provide scientific guidance enabling aquaculture farmers to reduce their economic losses.Accurate DO forecasting can also provide a scientific basis for fishery department managers to help them improve the living environment of ponds and marine life.It can also help decision-makers to regulate and restore ecological water quality.The proposed hybrid model has excellent multi-step forecasting performance. The MIMOS mechanism provides sufficient information about future DO fluctuations for relevant industry personnel, thereby allowing effective management.

    Acknowledgments

    The study was fully supported by the National Natural Science Foundation of China(61873283),the Changsha Science&Technology Project(KQ1707017),and the innovation-driven project of the Central South University (2019CX005).

    Compliance with ethics guidelines

    Hui Liu,Rui Yang,Zhu Duan, and Haiping Wu declare that they have no conflict of interest or financial conflicts to disclose.

    Appendix A. Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.eng.2020.10.023.

    xxxhd国产人妻xxx| 久久中文字幕人妻熟女| 亚洲美女黄片视频| 飞空精品影院首页| 午夜激情av网站| 国产亚洲欧美精品永久| 99热网站在线观看| tocl精华| 成年人黄色毛片网站| 国产熟女午夜一区二区三区| 青草久久国产| 日韩一卡2卡3卡4卡2021年| 波多野结衣一区麻豆| 国产色视频综合| a在线观看视频网站| 最近最新免费中文字幕在线| 青草久久国产| 黄色片一级片一级黄色片| 国产麻豆69| 超色免费av| 视频在线观看一区二区三区| 亚洲国产精品sss在线观看 | xxx96com| 免费在线观看视频国产中文字幕亚洲| 欧美成人午夜精品| tube8黄色片| 日本精品一区二区三区蜜桃| 久久精品国产亚洲av高清一级| 如日韩欧美国产精品一区二区三区| 18禁裸乳无遮挡动漫免费视频| 777米奇影视久久| 黄频高清免费视频| 在线观看日韩欧美| av片东京热男人的天堂| 欧美中文综合在线视频| 人人澡人人妻人| 国产精品久久久久成人av| 成人影院久久| svipshipincom国产片| 黄色毛片三级朝国网站| 国产高清国产精品国产三级| 国产精品乱码一区二三区的特点 | 亚洲五月天丁香| 99国产综合亚洲精品| 亚洲精品自拍成人| 国产黄色免费在线视频| 亚洲免费av在线视频| 高清av免费在线| 在线播放国产精品三级| 日日夜夜操网爽| 亚洲午夜理论影院| 手机成人av网站| 高清欧美精品videossex| 亚洲成人免费电影在线观看| 老熟妇乱子伦视频在线观看| 欧美日韩视频精品一区| 麻豆av在线久日| 国产免费现黄频在线看| 国产一卡二卡三卡精品| 在线观看免费视频网站a站| 国产男靠女视频免费网站| 新久久久久国产一级毛片| 天天躁夜夜躁狠狠躁躁| 天天影视国产精品| 亚洲视频免费观看视频| 捣出白浆h1v1| 国产有黄有色有爽视频| 最近最新免费中文字幕在线| 欧美日韩一级在线毛片| 老司机福利观看| 国产99久久九九免费精品| 不卡av一区二区三区| 一进一出抽搐gif免费好疼 | av不卡在线播放| 飞空精品影院首页| 久久久精品免费免费高清| 视频区欧美日本亚洲| 啪啪无遮挡十八禁网站| 制服诱惑二区| 欧美日韩视频精品一区| 精品亚洲成a人片在线观看| 黑人巨大精品欧美一区二区蜜桃| 丁香欧美五月| 久热这里只有精品99| 欧美日韩国产mv在线观看视频| 精品一区二区三区四区五区乱码| 大码成人一级视频| 国产成人av激情在线播放| 好看av亚洲va欧美ⅴa在| 国产乱人伦免费视频| 成人永久免费在线观看视频| 亚洲色图av天堂| 99国产精品一区二区三区| 国产精品 欧美亚洲| 精品无人区乱码1区二区| 捣出白浆h1v1| a级片在线免费高清观看视频| 精品亚洲成a人片在线观看| 老汉色∧v一级毛片| 欧美日韩精品网址| 色综合婷婷激情| 两个人看的免费小视频| 欧美成人免费av一区二区三区 | 男人操女人黄网站| 国产精品久久久久久人妻精品电影| 女人久久www免费人成看片| 国产欧美日韩一区二区三| 91成人精品电影| 成人av一区二区三区在线看| 亚洲视频免费观看视频| 香蕉丝袜av| 中亚洲国语对白在线视频| 国产精品久久视频播放| 国产蜜桃级精品一区二区三区 | 国产精品久久久久成人av| 一本综合久久免费| 亚洲专区国产一区二区| 91成人精品电影| 国产欧美日韩精品亚洲av| 啪啪无遮挡十八禁网站| 欧美黄色淫秽网站| 久久精品亚洲精品国产色婷小说| 亚洲视频免费观看视频| 人人妻人人澡人人看| 女人爽到高潮嗷嗷叫在线视频| 无遮挡黄片免费观看| cao死你这个sao货| 精品国产乱子伦一区二区三区| 日韩制服丝袜自拍偷拍| 国产av精品麻豆| e午夜精品久久久久久久| 欧美日韩瑟瑟在线播放| 欧美激情 高清一区二区三区| 一区在线观看完整版| 精品高清国产在线一区| 久久精品熟女亚洲av麻豆精品| 在线观看午夜福利视频| 看免费av毛片| 国产真人三级小视频在线观看| 日韩大码丰满熟妇| 在线国产一区二区在线| 80岁老熟妇乱子伦牲交| 极品人妻少妇av视频| 老司机亚洲免费影院| 在线观看免费视频网站a站| 久久精品亚洲av国产电影网| 国产精品久久久久成人av| 热99久久久久精品小说推荐| 人人妻人人澡人人看| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产超薄肉色丝袜足j| 国产亚洲欧美在线一区二区| 国产精品 欧美亚洲| 亚洲精品国产区一区二| 久久中文字幕人妻熟女| 人妻久久中文字幕网| 一区二区三区激情视频| 国产有黄有色有爽视频| 黄色a级毛片大全视频| 老熟女久久久| 日韩人妻精品一区2区三区| 人人澡人人妻人| 亚洲伊人色综图| 亚洲av日韩在线播放| av天堂久久9| 亚洲全国av大片| 亚洲国产精品一区二区三区在线| 亚洲 欧美一区二区三区| 精品亚洲成a人片在线观看| 最近最新免费中文字幕在线| 香蕉久久夜色| 12—13女人毛片做爰片一| 1024香蕉在线观看| 国产伦人伦偷精品视频| 桃红色精品国产亚洲av| 欧美日韩福利视频一区二区| 国产高清国产精品国产三级| 激情在线观看视频在线高清 | 日本wwww免费看| 久久人人爽av亚洲精品天堂| 国产高清视频在线播放一区| 亚洲av熟女| 9色porny在线观看| 啦啦啦免费观看视频1| 1024视频免费在线观看| 咕卡用的链子| 久久精品国产清高在天天线| 久久久久久亚洲精品国产蜜桃av| 精品国产超薄肉色丝袜足j| 大陆偷拍与自拍| 国产精品永久免费网站| 男女午夜视频在线观看| tocl精华| 亚洲专区国产一区二区| 国产一区二区三区在线臀色熟女 | 亚洲精品成人av观看孕妇| 五月开心婷婷网| 777米奇影视久久| 999精品在线视频| 一二三四社区在线视频社区8| 亚洲精品一二三| 免费在线观看完整版高清| 亚洲av第一区精品v没综合| 国产午夜精品久久久久久| bbb黄色大片| 一本大道久久a久久精品| 99re在线观看精品视频| 中国美女看黄片| 免费久久久久久久精品成人欧美视频| 国产乱人伦免费视频| av欧美777| 午夜福利在线观看吧| 欧美一级毛片孕妇| 十分钟在线观看高清视频www| 国产单亲对白刺激| 老司机亚洲免费影院| 国产真人三级小视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 成年女人毛片免费观看观看9 | 国产在视频线精品| 午夜视频精品福利| 亚洲av成人一区二区三| 久久国产乱子伦精品免费另类| av超薄肉色丝袜交足视频| av天堂久久9| 999久久久国产精品视频| 欧美日本中文国产一区发布| 亚洲成a人片在线一区二区| 久久精品国产综合久久久| 亚洲精品久久成人aⅴ小说| 香蕉国产在线看| 美女福利国产在线| 国产又色又爽无遮挡免费看| 满18在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 精品少妇一区二区三区视频日本电影| 久久亚洲精品不卡| 少妇裸体淫交视频免费看高清 | 午夜福利视频在线观看免费| 俄罗斯特黄特色一大片| 视频区欧美日本亚洲| 精品国产一区二区三区四区第35| 午夜福利欧美成人| 亚洲久久久国产精品| 大陆偷拍与自拍| 国产成人免费观看mmmm| 人人妻人人添人人爽欧美一区卜| 777久久人妻少妇嫩草av网站| 波多野结衣av一区二区av| 午夜亚洲福利在线播放| 三级毛片av免费| 国产成人免费无遮挡视频| 色婷婷久久久亚洲欧美| 精品一品国产午夜福利视频| 亚洲国产精品一区二区三区在线| 成人国语在线视频| 欧美乱码精品一区二区三区| 成年版毛片免费区| 黄色女人牲交| 在线看a的网站| 午夜老司机福利片| 黄色毛片三级朝国网站| 十八禁高潮呻吟视频| 黄色视频,在线免费观看| 777米奇影视久久| 国产一区二区三区视频了| 国产亚洲欧美在线一区二区| 亚洲精品乱久久久久久| 国产在视频线精品| 日韩 欧美 亚洲 中文字幕| 成在线人永久免费视频| 欧美亚洲 丝袜 人妻 在线| 两个人看的免费小视频| 久久人人97超碰香蕉20202| 中文字幕另类日韩欧美亚洲嫩草| 国产极品粉嫩免费观看在线| 亚洲精品av麻豆狂野| 国产在线观看jvid| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区视频了| 捣出白浆h1v1| 精品无人区乱码1区二区| 手机成人av网站| 在线天堂中文资源库| 国产精品亚洲av一区麻豆| 欧美日韩精品网址| 国产一区在线观看成人免费| 黑人操中国人逼视频| 新久久久久国产一级毛片| 日韩人妻精品一区2区三区| 成人18禁高潮啪啪吃奶动态图| 身体一侧抽搐| 欧美午夜高清在线| 亚洲成a人片在线一区二区| 亚洲成人免费电影在线观看| 久久国产精品大桥未久av| 正在播放国产对白刺激| 精品亚洲成国产av| 欧美成狂野欧美在线观看| 国产在视频线精品| 国产精品二区激情视频| 久久久久国产一级毛片高清牌| 一边摸一边做爽爽视频免费| 大香蕉久久网| 亚洲第一欧美日韩一区二区三区| av免费在线观看网站| 19禁男女啪啪无遮挡网站| 1024视频免费在线观看| 亚洲av欧美aⅴ国产| 人妻丰满熟妇av一区二区三区 | 丰满的人妻完整版| 在线观看免费午夜福利视频| 超碰成人久久| 99精品欧美一区二区三区四区| 亚洲精品自拍成人| 国产三级黄色录像| 久久草成人影院| 久久久久国产精品人妻aⅴ院 | 成人18禁在线播放| 伊人久久大香线蕉亚洲五| 免费人成视频x8x8入口观看| 亚洲国产欧美网| 精品乱码久久久久久99久播| 美女 人体艺术 gogo| 老司机深夜福利视频在线观看| 免费黄频网站在线观看国产| 在线观看午夜福利视频| 久久国产精品人妻蜜桃| 日韩欧美国产一区二区入口| 亚洲七黄色美女视频| √禁漫天堂资源中文www| 纯流量卡能插随身wifi吗| 免费观看人在逋| 国产精品欧美亚洲77777| 日韩欧美国产一区二区入口| 天天操日日干夜夜撸| 50天的宝宝边吃奶边哭怎么回事| 免费黄频网站在线观看国产| 这个男人来自地球电影免费观看| 淫妇啪啪啪对白视频| 天天影视国产精品| 国产精品美女特级片免费视频播放器 | 黄色毛片三级朝国网站| 中文字幕色久视频| 岛国在线观看网站| 欧美黄色片欧美黄色片| 国产成人一区二区三区免费视频网站| 欧美日韩国产mv在线观看视频| 国产xxxxx性猛交| 国产精品美女特级片免费视频播放器 | 欧美在线黄色| 午夜精品久久久久久毛片777| 久热爱精品视频在线9| 成人亚洲精品一区在线观看| 国产激情久久老熟女| 男人舔女人的私密视频| 国产精品 国内视频| 日韩人妻精品一区2区三区| 高潮久久久久久久久久久不卡| 日韩成人在线观看一区二区三区| 欧美人与性动交α欧美精品济南到| 久久精品亚洲熟妇少妇任你| 极品人妻少妇av视频| 狠狠狠狠99中文字幕| bbb黄色大片| 九色亚洲精品在线播放| 国产av精品麻豆| 一级毛片女人18水好多| 天天躁狠狠躁夜夜躁狠狠躁| 电影成人av| 久久国产亚洲av麻豆专区| 在线免费观看的www视频| 校园春色视频在线观看| 自线自在国产av| 国产av又大| 女性被躁到高潮视频| 一进一出好大好爽视频| 久久草成人影院| 亚洲伊人色综图| 波多野结衣一区麻豆| 欧美精品亚洲一区二区| 久久国产精品大桥未久av| 中文字幕最新亚洲高清| 国产一区有黄有色的免费视频| 制服诱惑二区| 黄色毛片三级朝国网站| 女人精品久久久久毛片| 国产精品一区二区精品视频观看| 少妇裸体淫交视频免费看高清 | 国产男女超爽视频在线观看| 91大片在线观看| 久久国产乱子伦精品免费另类| 黄色片一级片一级黄色片| 飞空精品影院首页| 欧美精品人与动牲交sv欧美| 亚洲午夜理论影院| 亚洲成人国产一区在线观看| 国产亚洲欧美98| 欧美亚洲日本最大视频资源| 日本a在线网址| 看片在线看免费视频| videos熟女内射| 午夜精品国产一区二区电影| 国产91精品成人一区二区三区| 91麻豆av在线| 男女午夜视频在线观看| 国产av又大| 色老头精品视频在线观看| 国产午夜精品久久久久久| 老司机午夜福利在线观看视频| 国产无遮挡羞羞视频在线观看| 性色av乱码一区二区三区2| av国产精品久久久久影院| 国产又色又爽无遮挡免费看| 久99久视频精品免费| 国产免费男女视频| 老司机影院毛片| 91精品国产国语对白视频| 久久久水蜜桃国产精品网| 久久国产精品人妻蜜桃| 亚洲成人国产一区在线观看| 国产高清国产精品国产三级| av欧美777| 日本wwww免费看| 欧美 日韩 精品 国产| 制服诱惑二区| 不卡一级毛片| 最新的欧美精品一区二区| 亚洲国产精品合色在线| 日韩欧美免费精品| av电影中文网址| 欧美日韩中文字幕国产精品一区二区三区 | 9色porny在线观看| 91麻豆精品激情在线观看国产 | 久久影院123| 午夜福利视频在线观看免费| 在线观看日韩欧美| 欧美成人午夜精品| 19禁男女啪啪无遮挡网站| 一a级毛片在线观看| 久久久国产成人免费| 亚洲专区国产一区二区| 性色av乱码一区二区三区2| 91在线观看av| 国产不卡av网站在线观看| 极品人妻少妇av视频| 国产高清视频在线播放一区| 午夜福利在线观看吧| 久久亚洲真实| 成人影院久久| avwww免费| 桃红色精品国产亚洲av| 欧美日韩乱码在线| 国产精品 欧美亚洲| 精品一区二区三卡| 日韩一卡2卡3卡4卡2021年| 美女视频免费永久观看网站| 精品福利永久在线观看| 美国免费a级毛片| 日韩熟女老妇一区二区性免费视频| 国产片内射在线| 伊人久久大香线蕉亚洲五| 欧美乱妇无乱码| 麻豆成人av在线观看| 久久中文字幕人妻熟女| 国产男女超爽视频在线观看| 国产97色在线日韩免费| 在线观看www视频免费| 男女下面插进去视频免费观看| 老司机深夜福利视频在线观看| 久久精品aⅴ一区二区三区四区| 亚洲国产中文字幕在线视频| 久久精品国产亚洲av香蕉五月 | 极品少妇高潮喷水抽搐| 久久久精品国产亚洲av高清涩受| 亚洲专区中文字幕在线| 免费在线观看亚洲国产| 欧美日韩亚洲高清精品| 99精品久久久久人妻精品| 亚洲专区中文字幕在线| 亚洲五月婷婷丁香| 亚洲精品国产区一区二| 99香蕉大伊视频| 韩国精品一区二区三区| 99久久精品国产亚洲精品| 免费黄频网站在线观看国产| 欧美一级毛片孕妇| 日本欧美视频一区| 婷婷丁香在线五月| 午夜影院日韩av| 久久这里只有精品19| 激情视频va一区二区三区| 最近最新中文字幕大全免费视频| 香蕉久久夜色| 叶爱在线成人免费视频播放| 亚洲av成人一区二区三| 人人妻人人澡人人爽人人夜夜| 一级毛片高清免费大全| 国产欧美亚洲国产| 亚洲 欧美一区二区三区| 中文字幕人妻丝袜一区二区| 亚洲在线自拍视频| 人人妻人人爽人人添夜夜欢视频| 国产高清国产精品国产三级| av视频免费观看在线观看| 国产精品久久电影中文字幕 | 国产视频一区二区在线看| 久久精品91无色码中文字幕| 国产精品亚洲一级av第二区| 91国产中文字幕| 亚洲精品国产区一区二| 免费在线观看黄色视频的| 中文字幕制服av| 免费在线观看黄色视频的| 国产成人欧美| 视频区图区小说| 国产男靠女视频免费网站| 精品久久蜜臀av无| avwww免费| 制服人妻中文乱码| 亚洲久久久国产精品| 日日爽夜夜爽网站| 美女高潮到喷水免费观看| 亚洲 欧美一区二区三区| 黄色视频不卡| 人人妻人人添人人爽欧美一区卜| 人人妻人人爽人人添夜夜欢视频| 91大片在线观看| 咕卡用的链子| 亚洲欧美激情在线| 丝袜人妻中文字幕| 国产视频一区二区在线看| av超薄肉色丝袜交足视频| 侵犯人妻中文字幕一二三四区| 男女下面插进去视频免费观看| 下体分泌物呈黄色| 夜夜爽天天搞| 欧美精品亚洲一区二区| 一边摸一边抽搐一进一小说 | 91精品三级在线观看| а√天堂www在线а√下载 | 麻豆成人av在线观看| 久久精品91无色码中文字幕| 一边摸一边做爽爽视频免费| 国产精品成人在线| 嫩草影视91久久| 免费在线观看亚洲国产| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 午夜影院日韩av| 欧美成人免费av一区二区三区 | 欧美乱妇无乱码| 中文字幕制服av| 欧美丝袜亚洲另类 | 一二三四社区在线视频社区8| 99精品在免费线老司机午夜| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲色图av天堂| 桃红色精品国产亚洲av| 亚洲精品国产区一区二| 99热国产这里只有精品6| 黄色成人免费大全| 亚洲三区欧美一区| av网站免费在线观看视频| 精品卡一卡二卡四卡免费| 亚洲,欧美精品.| 成人18禁高潮啪啪吃奶动态图| 在线免费观看的www视频| 老司机福利观看| 十八禁高潮呻吟视频| 国产单亲对白刺激| 啦啦啦免费观看视频1| 多毛熟女@视频| 久久精品国产清高在天天线| 亚洲色图综合在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 日韩有码中文字幕| 飞空精品影院首页| 欧美日韩福利视频一区二区| av线在线观看网站| 女人精品久久久久毛片| 在线视频色国产色| 国产激情久久老熟女| 亚洲第一青青草原| 啦啦啦在线免费观看视频4| 一个人免费在线观看的高清视频| av电影中文网址| 啪啪无遮挡十八禁网站| 日本黄色视频三级网站网址 | 男人舔女人的私密视频| 亚洲精品乱久久久久久| 人妻 亚洲 视频| 国产99久久九九免费精品| 久久国产精品大桥未久av| 美女国产高潮福利片在线看| 丝袜人妻中文字幕| 久久久久久人人人人人| 18禁裸乳无遮挡动漫免费视频| 国产淫语在线视频| 夜夜躁狠狠躁天天躁| 中文字幕另类日韩欧美亚洲嫩草| 变态另类成人亚洲欧美熟女 | 欧美亚洲日本最大视频资源| 精品电影一区二区在线| tocl精华| 亚洲精品在线观看二区| 国产亚洲欧美精品永久| 久久精品国产清高在天天线| 丁香六月欧美| а√天堂www在线а√下载 | 亚洲av日韩精品久久久久久密| 欧美日韩瑟瑟在线播放| 午夜精品在线福利| 午夜免费成人在线视频|