• <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.

    国产伦精品一区二区三区视频9 | 国产真人三级小视频在线观看| 99热这里只有是精品50| 久久久国产成人精品二区| 亚洲五月婷婷丁香| 黄色视频,在线免费观看| 成人特级黄色片久久久久久久| 在线播放国产精品三级| tocl精华| 亚洲熟妇中文字幕五十中出| 此物有八面人人有两片| 天天一区二区日本电影三级| 国产成人啪精品午夜网站| 高清毛片免费观看视频网站| 久久久水蜜桃国产精品网| 午夜免费成人在线视频| 成人国产综合亚洲| 观看免费一级毛片| 欧美中文日本在线观看视频| 亚洲欧美日韩卡通动漫| 少妇丰满av| 亚洲成人久久性| av片东京热男人的天堂| 久久精品综合一区二区三区| 天天一区二区日本电影三级| 亚洲欧洲精品一区二区精品久久久| 日韩欧美 国产精品| 韩国av一区二区三区四区| 色av中文字幕| 色视频www国产| 国产精品一区二区三区四区久久| 亚洲 国产 在线| 亚洲精品在线美女| 欧美一级毛片孕妇| 99在线视频只有这里精品首页| 香蕉av资源在线| 嫩草影视91久久| 欧美精品啪啪一区二区三区| 黄色片一级片一级黄色片| 亚洲欧美一区二区三区黑人| 日本 av在线| 1000部很黄的大片| 国产精品99久久久久久久久| 久久久国产成人精品二区| 国产成人影院久久av| 噜噜噜噜噜久久久久久91| 久久久久亚洲av毛片大全| 欧美激情久久久久久爽电影| 中文字幕最新亚洲高清| bbb黄色大片| 岛国在线免费视频观看| 男女床上黄色一级片免费看| 中文字幕人成人乱码亚洲影| 成人av一区二区三区在线看| 亚洲成人免费电影在线观看| 亚洲欧美日韩高清在线视频| 亚洲人与动物交配视频| 免费搜索国产男女视频| 亚洲精品美女久久av网站| 日韩成人在线观看一区二区三区| 日韩中文字幕欧美一区二区| 欧美日本视频| 国产人伦9x9x在线观看| 最新在线观看一区二区三区| 免费观看的影片在线观看| 国产伦人伦偷精品视频| 天堂动漫精品| 国产不卡一卡二| 国产爱豆传媒在线观看| 国产乱人伦免费视频| 一区二区三区国产精品乱码| www.熟女人妻精品国产| 欧美色欧美亚洲另类二区| 国产真实乱freesex| 麻豆av在线久日| 波多野结衣高清作品| 国产免费av片在线观看野外av| 亚洲国产精品999在线| 欧美激情久久久久久爽电影| 啦啦啦韩国在线观看视频| 中文字幕最新亚洲高清| 在线看三级毛片| 国模一区二区三区四区视频 | 欧美乱妇无乱码| 午夜亚洲福利在线播放| 亚洲,欧美精品.| 97人妻精品一区二区三区麻豆| 日韩欧美 国产精品| 国产伦一二天堂av在线观看| 久久久久久久久免费视频了| www国产在线视频色| 亚洲人成电影免费在线| 色精品久久人妻99蜜桃| 午夜亚洲福利在线播放| 波多野结衣高清作品| 久久精品国产亚洲av香蕉五月| 三级国产精品欧美在线观看 | 国产一区二区在线av高清观看| 黄色女人牲交| 亚洲 欧美 日韩 在线 免费| 窝窝影院91人妻| 亚洲第一电影网av| 美女cb高潮喷水在线观看 | 日本黄色片子视频| 十八禁人妻一区二区| 精品国产乱码久久久久久男人| 99在线人妻在线中文字幕| 国产高清激情床上av| 99国产极品粉嫩在线观看| 黄色 视频免费看| 久久香蕉精品热| 日日干狠狠操夜夜爽| avwww免费| 色在线成人网| 91字幕亚洲| 国内精品美女久久久久久| 黄色成人免费大全| 少妇熟女aⅴ在线视频| 午夜日韩欧美国产| 搡老熟女国产l中国老女人| 国产单亲对白刺激| 国产真实乱freesex| 这个男人来自地球电影免费观看| 亚洲七黄色美女视频| 国产精品自产拍在线观看55亚洲| 亚洲熟妇中文字幕五十中出| 啦啦啦免费观看视频1| 久久精品夜夜夜夜夜久久蜜豆| 男女下面进入的视频免费午夜| 国内精品一区二区在线观看| 午夜日韩欧美国产| 亚洲 欧美一区二区三区| 91麻豆av在线| 国内精品久久久久久久电影| 久久九九热精品免费| 亚洲精品久久国产高清桃花| 国内毛片毛片毛片毛片毛片| 亚洲精品粉嫩美女一区| 国产又色又爽无遮挡免费看| 一a级毛片在线观看| 国产又色又爽无遮挡免费看| 巨乳人妻的诱惑在线观看| 岛国在线观看网站| 网址你懂的国产日韩在线| 国产伦一二天堂av在线观看| 亚洲无线观看免费| 老司机深夜福利视频在线观看| 国产成+人综合+亚洲专区| www.熟女人妻精品国产| 国产成人啪精品午夜网站| 亚洲avbb在线观看| 午夜福利成人在线免费观看| 国产单亲对白刺激| 成年女人毛片免费观看观看9| 在线看三级毛片| 老司机午夜十八禁免费视频| 午夜激情欧美在线| aaaaa片日本免费| 免费在线观看成人毛片| 非洲黑人性xxxx精品又粗又长| 午夜激情福利司机影院| 三级毛片av免费| 18禁黄网站禁片午夜丰满| 一进一出抽搐动态| 国产麻豆成人av免费视频| 无人区码免费观看不卡| 可以在线观看的亚洲视频| 精品一区二区三区视频在线 | 巨乳人妻的诱惑在线观看| 精品久久久久久,| 国产蜜桃级精品一区二区三区| 国产一区二区在线观看日韩 | 亚洲第一电影网av| 欧美国产日韩亚洲一区| 夜夜夜夜夜久久久久| 久久久久性生活片| 国产精品女同一区二区软件 | 国产精品av视频在线免费观看| 草草在线视频免费看| 色综合亚洲欧美另类图片| 婷婷六月久久综合丁香| 国内揄拍国产精品人妻在线| 亚洲美女视频黄频| 一进一出好大好爽视频| 亚洲成人中文字幕在线播放| 18禁黄网站禁片免费观看直播| 婷婷精品国产亚洲av在线| 国产一区在线观看成人免费| 99久久久亚洲精品蜜臀av| 两个人的视频大全免费| 男女午夜视频在线观看| 村上凉子中文字幕在线| 午夜影院日韩av| www日本黄色视频网| 九九久久精品国产亚洲av麻豆 | 亚洲性夜色夜夜综合| 观看美女的网站| 夜夜爽天天搞| 久久午夜亚洲精品久久| 久久九九热精品免费| 久久午夜亚洲精品久久| 国产淫片久久久久久久久 | 色综合欧美亚洲国产小说| svipshipincom国产片| 国产精品亚洲av一区麻豆| 91麻豆av在线| 夜夜躁狠狠躁天天躁| 亚洲av中文字字幕乱码综合| 国产 一区 欧美 日韩| 国产人伦9x9x在线观看| 国产成人啪精品午夜网站| 夜夜躁狠狠躁天天躁| 国产精品野战在线观看| 亚洲中文日韩欧美视频| 超碰成人久久| 欧美不卡视频在线免费观看| 久久国产精品影院| 国内精品美女久久久久久| 伊人久久大香线蕉亚洲五| 亚洲欧美日韩无卡精品| 少妇丰满av| 免费高清视频大片| 90打野战视频偷拍视频| 免费搜索国产男女视频| 国产亚洲av嫩草精品影院| 亚洲成人久久爱视频| 18禁黄网站禁片午夜丰满| 免费在线观看亚洲国产| 欧美日韩亚洲国产一区二区在线观看| 亚洲欧洲精品一区二区精品久久久| 一级毛片女人18水好多| 国产伦精品一区二区三区四那| 久久精品亚洲精品国产色婷小说| 动漫黄色视频在线观看| 国产精品香港三级国产av潘金莲| 亚洲美女视频黄频| 中文字幕人妻丝袜一区二区| 精品久久久久久久毛片微露脸| 欧美成人性av电影在线观看| 欧美乱妇无乱码| 最新在线观看一区二区三区| www.熟女人妻精品国产| 成人欧美大片| 又黄又爽又免费观看的视频| 国产视频一区二区在线看| 久久精品人妻少妇| 美女免费视频网站| 国产精品,欧美在线| 免费大片18禁| 毛片女人毛片| 国产97色在线日韩免费| 一卡2卡三卡四卡精品乱码亚洲| 亚洲 欧美一区二区三区| 国产精品一区二区精品视频观看| 1000部很黄的大片| 欧美一区二区精品小视频在线| 亚洲真实伦在线观看| 婷婷六月久久综合丁香| 国产成年人精品一区二区| 亚洲av日韩精品久久久久久密| 亚洲七黄色美女视频| www.www免费av| 美女 人体艺术 gogo| 色av中文字幕| 伦理电影免费视频| 国产精品野战在线观看| 日本成人三级电影网站| 伊人久久大香线蕉亚洲五| 欧美日本亚洲视频在线播放| 国产主播在线观看一区二区| www.999成人在线观看| 亚洲无线在线观看| 欧美xxxx黑人xx丫x性爽| 三级毛片av免费| 色尼玛亚洲综合影院| 色综合婷婷激情| 国产亚洲av高清不卡| 成人国产一区最新在线观看| 国产不卡一卡二| 国产一区二区激情短视频| 啪啪无遮挡十八禁网站| 在线观看免费午夜福利视频| 国产激情久久老熟女| 亚洲av中文字字幕乱码综合| 老熟妇乱子伦视频在线观看| 亚洲成人久久性| 久久这里只有精品中国| 一个人免费在线观看的高清视频| 一本一本综合久久| 亚洲成a人片在线一区二区| 九色国产91popny在线| 国内精品久久久久精免费| 黄色日韩在线| ponron亚洲| 久久精品夜夜夜夜夜久久蜜豆| 国产精品免费一区二区三区在线| 亚洲七黄色美女视频| 免费观看的影片在线观看| 免费看日本二区| 天堂av国产一区二区熟女人妻| 在线观看66精品国产| 国产av在哪里看| 亚洲精品美女久久久久99蜜臀| 日韩av在线大香蕉| 美女 人体艺术 gogo| 久久人妻av系列| 日韩国内少妇激情av| 日韩人妻高清精品专区| 免费电影在线观看免费观看| 一夜夜www| 99热6这里只有精品| 日本一本二区三区精品| 在线播放国产精品三级| 熟女少妇亚洲综合色aaa.| 欧美日韩国产亚洲二区| 国产男靠女视频免费网站| 美女高潮喷水抽搐中文字幕| 欧美成人性av电影在线观看| 亚洲人成伊人成综合网2020| 天堂网av新在线| 国产视频一区二区在线看| 草草在线视频免费看| 久久久久久久精品吃奶| 三级毛片av免费| 伦理电影免费视频| 成年版毛片免费区| 午夜福利18| 国产91精品成人一区二区三区| 欧美成人性av电影在线观看| 一进一出抽搐动态| 国产高清三级在线| 手机成人av网站| 欧美av亚洲av综合av国产av| 久久香蕉国产精品| 给我免费播放毛片高清在线观看| 欧美成人免费av一区二区三区| 小蜜桃在线观看免费完整版高清| 不卡一级毛片| 国产私拍福利视频在线观看| 亚洲国产精品成人综合色| 成年女人毛片免费观看观看9| 亚洲成av人片在线播放无| 欧美一区二区精品小视频在线| 哪里可以看免费的av片| 99在线视频只有这里精品首页| 久久伊人香网站| 黑人操中国人逼视频| 午夜福利高清视频| 亚洲av免费在线观看| 麻豆国产av国片精品| 国产熟女xx| 亚洲欧美日韩无卡精品| 日韩高清综合在线| 丁香六月欧美| avwww免费| 真人一进一出gif抽搐免费| 成人18禁在线播放| 老熟妇乱子伦视频在线观看| 搡老岳熟女国产| 制服丝袜大香蕉在线| 97超级碰碰碰精品色视频在线观看| 国产1区2区3区精品| 亚洲成人久久性| 长腿黑丝高跟| 香蕉国产在线看| 男女那种视频在线观看| 又紧又爽又黄一区二区| 在线观看一区二区三区| 最近最新中文字幕大全免费视频| 狂野欧美激情性xxxx| 国产精品一区二区精品视频观看| 亚洲精品粉嫩美女一区| 亚洲男人的天堂狠狠| 午夜免费成人在线视频| 久9热在线精品视频| 手机成人av网站| 露出奶头的视频| 在线播放国产精品三级| 午夜免费成人在线视频| 久久久精品欧美日韩精品| 国内精品久久久久久久电影| 人人妻人人看人人澡| 人人妻,人人澡人人爽秒播| 精品久久久久久久毛片微露脸| 成人鲁丝片一二三区免费| 亚洲欧美日韩东京热| 99re在线观看精品视频| 色播亚洲综合网| 午夜久久久久精精品| 日本黄大片高清| 亚洲精品美女久久久久99蜜臀| 久久午夜亚洲精品久久| 久久午夜综合久久蜜桃| 亚洲成人免费电影在线观看| 五月玫瑰六月丁香| 好男人在线观看高清免费视频| 欧美最黄视频在线播放免费| 麻豆久久精品国产亚洲av| 久久精品影院6| 亚洲成人久久爱视频| 亚洲av成人精品一区久久| 久久精品91蜜桃| 色综合站精品国产| 国产成人av教育| 国产亚洲精品综合一区在线观看| av中文乱码字幕在线| av视频在线观看入口| 丰满人妻熟妇乱又伦精品不卡| 一个人看视频在线观看www免费 | 免费在线观看视频国产中文字幕亚洲| 国产成人系列免费观看| 欧美在线一区亚洲| 精品久久久久久久久久久久久| 久久精品91蜜桃| 亚洲精品国产精品久久久不卡| 久久这里只有精品19| 无限看片的www在线观看| 色噜噜av男人的天堂激情| 中国美女看黄片| 欧美一级a爱片免费观看看| 日韩免费av在线播放| 极品教师在线免费播放| 特大巨黑吊av在线直播| 黄色日韩在线| 亚洲国产欧洲综合997久久,| 国产精品美女特级片免费视频播放器 | 最新在线观看一区二区三区| 老司机深夜福利视频在线观看| 国产一区二区在线观看日韩 | 波多野结衣高清无吗| 一夜夜www| av黄色大香蕉| 国产高潮美女av| 99国产综合亚洲精品| 欧美一级毛片孕妇| 亚洲中文字幕一区二区三区有码在线看 | 亚洲欧美日韩卡通动漫| 日韩欧美在线二视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲欧美日韩卡通动漫| 日本撒尿小便嘘嘘汇集6| 精品久久久久久久人妻蜜臀av| 日本与韩国留学比较| 18禁美女被吸乳视频| 亚洲国产中文字幕在线视频| 亚洲成人免费电影在线观看| 久久香蕉国产精品| 成年女人毛片免费观看观看9| 久久国产精品影院| 日韩三级视频一区二区三区| 国产黄a三级三级三级人| 热99re8久久精品国产| www日本在线高清视频| 国产高清视频在线观看网站| 亚洲一区高清亚洲精品| 亚洲专区国产一区二区| 久久香蕉国产精品| 老鸭窝网址在线观看| 国产精品 欧美亚洲| 老熟妇乱子伦视频在线观看| 国产爱豆传媒在线观看| 91麻豆精品激情在线观看国产| 啪啪无遮挡十八禁网站| 真实男女啪啪啪动态图| 午夜福利在线观看免费完整高清在 | 丰满人妻熟妇乱又伦精品不卡| 18禁观看日本| 性色av乱码一区二区三区2| 久久草成人影院| 亚洲中文字幕一区二区三区有码在线看 | 国产精品女同一区二区软件 | 免费观看精品视频网站| 九九久久精品国产亚洲av麻豆 | 国产精品98久久久久久宅男小说| 午夜亚洲福利在线播放| 人人妻人人看人人澡| 久久香蕉国产精品| 国产又黄又爽又无遮挡在线| 成人av一区二区三区在线看| 亚洲av成人av| 美女高潮的动态| 精品不卡国产一区二区三区| 淫秽高清视频在线观看| 亚洲欧美精品综合久久99| 麻豆久久精品国产亚洲av| 叶爱在线成人免费视频播放| 午夜久久久久精精品| 好男人电影高清在线观看| 18美女黄网站色大片免费观看| 国产探花在线观看一区二区| 此物有八面人人有两片| 国产精品久久久久久人妻精品电影| 亚洲在线自拍视频| 午夜福利18| 国模一区二区三区四区视频 | 日韩欧美三级三区| 国产精品亚洲av一区麻豆| 一级黄色大片毛片| 国产欧美日韩精品亚洲av| 精品久久蜜臀av无| 欧美一区二区精品小视频在线| 色播亚洲综合网| 9191精品国产免费久久| 欧美午夜高清在线| 首页视频小说图片口味搜索| 身体一侧抽搐| 久久亚洲真实| 亚洲无线观看免费| 长腿黑丝高跟| svipshipincom国产片| 人人妻人人澡欧美一区二区| 午夜精品在线福利| 网址你懂的国产日韩在线| 国产精品自产拍在线观看55亚洲| 99热这里只有是精品50| 成人三级黄色视频| 很黄的视频免费| 色哟哟哟哟哟哟| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品一区二区www| 女生性感内裤真人,穿戴方法视频| 伦理电影免费视频| 国产又色又爽无遮挡免费看| 首页视频小说图片口味搜索| 老汉色∧v一级毛片| 亚洲精华国产精华精| 又黄又粗又硬又大视频| 丰满人妻一区二区三区视频av | 变态另类丝袜制服| 黄色丝袜av网址大全| 啦啦啦免费观看视频1| 在线观看免费午夜福利视频| 麻豆成人av在线观看| 精品免费久久久久久久清纯| 国产1区2区3区精品| 岛国视频午夜一区免费看| 久久伊人香网站| 美女大奶头视频| 老司机午夜十八禁免费视频| 每晚都被弄得嗷嗷叫到高潮| 又大又爽又粗| 国语自产精品视频在线第100页| 香蕉丝袜av| 国产亚洲精品综合一区在线观看| 欧美黄色片欧美黄色片| 熟女少妇亚洲综合色aaa.| 国产野战对白在线观看| 热99re8久久精品国产| 我的老师免费观看完整版| 又大又爽又粗| 免费无遮挡裸体视频| 亚洲精品中文字幕一二三四区| 99久久99久久久精品蜜桃| 亚洲精品粉嫩美女一区| 淫妇啪啪啪对白视频| 久久精品国产清高在天天线| 国产成人啪精品午夜网站| 搡老岳熟女国产| 最新在线观看一区二区三区| 日韩精品青青久久久久久| 欧美色视频一区免费| av福利片在线观看| а√天堂www在线а√下载| 久久久久免费精品人妻一区二区| 女生性感内裤真人,穿戴方法视频| 欧美国产日韩亚洲一区| 亚洲美女黄片视频| 亚洲五月天丁香| 色老头精品视频在线观看| 亚洲九九香蕉| 久9热在线精品视频| 亚洲av五月六月丁香网| 久久这里只有精品中国| 成人亚洲精品av一区二区| 日韩欧美精品v在线| 久久久久性生活片| 欧美三级亚洲精品| 午夜福利在线观看吧| 久久精品国产综合久久久| 一个人免费在线观看的高清视频| 日日夜夜操网爽| 草草在线视频免费看| www.精华液| av欧美777| xxxwww97欧美| 女人高潮潮喷娇喘18禁视频| 男人的好看免费观看在线视频| 久久国产精品影院| 午夜福利在线观看吧| 国产成人精品久久二区二区免费| 欧美一区二区国产精品久久精品| 窝窝影院91人妻| 亚洲性夜色夜夜综合| 一a级毛片在线观看| 一进一出抽搐动态| 日韩精品青青久久久久久| 欧美成人一区二区免费高清观看 | 国产伦一二天堂av在线观看| 黄色视频,在线免费观看| 国产人伦9x9x在线观看| 久久精品国产亚洲av香蕉五月| 在线观看日韩欧美| 国语自产精品视频在线第100页| 亚洲av熟女| 国产伦在线观看视频一区| 国产伦人伦偷精品视频| 黄色日韩在线| 久久精品91无色码中文字幕| 国产一级毛片七仙女欲春2| 国产淫片久久久久久久久 | 久久国产精品影院| 亚洲,欧美精品.| 淫妇啪啪啪对白视频| 十八禁网站免费在线| 久久午夜综合久久蜜桃|