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

    Research on Condenser Deterioration Evolution Trend Based on ANP-EWM Fusion Health Degree

    2024-01-20 13:01:56HongQianHaixinWangGuangjiWangandQingyunYan

    Hong Qian ,Haixin Wang ,Guangji Wang and Qingyun Yan

    1College of Automation Engineering,Shanghai University of Electric Power,Shanghai,200090,China

    2Shanghai Key Laboratory of Power Station Automation Technology,Shanghai,200072,China

    3Empower Information Technology(Shanghai)Co.,Ltd.,Shanghai,201508,China

    4China Nuclear Power Engineering Design Co.,Ltd.,Shenzhen,518052,China

    ABSTRACT This study presents a proposed method for assessing the condition and predicting the future status of condensers operating in seawater over an extended period.The aim is to address the problems of scaling and corrosion,which lead to increased loss of cold resources.The method involves utilising a set of multivariate feature parameters associated with the condenser as input for evaluation and trend prediction.This methodology offers a precise means of determining the optimal timing for condenser cleaning,with the ultimate goal of improving its overall performance.The proposed approach involves the integration of the analytic network process with subjective expert experience and the entropy weight method with objective big data analysis to develop a fusion health degree model.The mathematical model is constructed quantitatively using the improved Mahalanobis distance.Furthermore,a comprehensive prediction model is developed by integrating the improved Informer model and Markov error correction.This model takes into account the health status of the equipment and several influencing factors,including multivariate feature characteristics.This model facilitates the objective examination and prediction of the progression of equipment deterioration trends.The present study involves the computation and verification of the field time series data,which serves to demonstrate the accuracy of the condenser health-related models proposed in this research.These models effectively depict the real condition and temporal variations of the equipment,thus offering a valuable method for determining the precise cleaning time required for the condenser.

    KEYWORDS Condenser;health degree;improved Mahalanobis distance;GC-Informer model;Markov error correction

    1 Introduction

    In the preceding thirty years,the secondary loop systems of nuclear power plants have encountered diverse levels of degradation and deterioration in their crucial equipment.Hence,conducting research on the health status of nuclear power unit equipment carries significant importance[1].The condenser,being a crucial heat exchange apparatus within the conventional island of nuclear power plants,has a direct influence on the safety of the turbine unit by ensuring its proper and steady operation.The extended exposure of condensers to seawater might result in the occurrence of scaling and corrosion.Presently,the prevailing method for cleaning condensers involves the utilisation of conventional online rubber ball cleaning equipment in order to uphold cleanliness and optimise heat transfer efficiency.Nevertheless,the utilisation of rubber balls for cleaning purposes incurs a substantial financial burden,while inadequate cleaning intervals,whether too prolonged or excessively brief,can lead to unnecessary energy consumption [2,3].The determination of condenser cleaning time holds considerable importance in enhancing the operational dependability of the condensing system and the economic efficiency of the turbine unit.Therefore,it is important to conduct monitoring and analysis of the condenser’s status in order to determine the appropriate timing for cleaning the condenser equipment.This can be achieved by quantitatively assessing the degree of equipment health.

    For the purpose of successfully cleaning the condenser,it is essential to accurately evaluate its overall condition.Currently,the predominant approach employed to assess the condition of the equipment is based on the notion of health degree.Nevertheless,the concept of a “healthy degree”lacks quantifiability and necessitates the identification and analysis of pertinent aspects that contribute to its attainment.The establishment of a scientifically valid assessment approach for determining the status of equipment is deemed required [4,5].The assessment methods can be broadly categorised into subjective weighting methods and objective weighting methods.Subjective weighing methods are predominantly exemplified by two prominent approaches,namely the Analytic Hierarchy Process(AHP) and the Fuzzy Comprehensive Evaluation method (FCE).The literature employs the use of AHP[6,7]in order to establish an evaluation model that is based on scientific principles.This approach presents the benefits of a straightforward procedure and yielding intuitive outcomes.Nevertheless,the susceptibility of this phenomenon to subjective factors,such as expert knowledge and experience,should be acknowledged.As a result,it is possible that this could result in inaccuracies when evaluating the outcomes.The FCE[8]methodology employs fuzzy mathematics in order to thoroughly evaluate a wide range of criteria that pose challenges in terms of quantitative assessment.In their study,Hassan et al.[9]utilised a Fuzzy Analytic Hierarchy Process(FAHP)methodology to evaluate the level of satisfaction among borrowers.Nevertheless,their methodology exclusively facilitates a thorough fuzzy assessment,without immediately offering a quantitative depiction of the level of health of the item.The Principal Component Analysis (PCA) and the Entropy Weighting Method (EWM) are the primary objective weighting approaches.In their study,Giorgi et al.[10] employed an enhanced Principal Component Analysis (PCA) technique to extract the feature parameters that carry the most relevant information for the purpose of monitoring the health status of engines.Nevertheless,the use of PCA involves the conversion of correlated indicators into principle components,hence resulting in the potential loss of certain information.This loss of information can subsequently lead to the wrong determination of weights.In their study,Erkhembaatar et al.[11] employed the entropy weighting approach to ascertain the significance coefficients of indicators.In their study,Zavadskas et al.[12] utilised the Entropy Weighting Method to evaluate alternative solutions.Their objective was to identify the most optimal alternative solution by identifying the most suitable criteria.In their study,Gorgij et al.[13]employed the Exponential Weighted Moving Average(EWM)technique to construct a model for evaluating groundwater quality.The researchers evaluated the weights of pertinent factors in order to assess the water quality level of 21 groundwater samples.This model is well-suited for scenarios that involve several indications,since it has the capability to objectively reflect the significance of these indicators based on big data analysis.Nevertheless,it is worth noting that this issue can be delicate,potentially resulting in the inefficiency of parameters.Hence,there is a need for a more precise approach to ascertain the weights of feature parameters in order to quantify the degree of equipment health.By integrating the benefits of subjective weighting methods and objective weighting methods,a more comprehensive depiction of equipment health status can be attained through mutually reinforcing approaches.

    The assessment of the operational condition and estimation of the cleaning duration of the condenser are crucial tasks in evaluating the health status of the equipment,as per the healthy degree model.In their study,Ma [14] focused on the determination of the equipment health index and its status after a one-year period.This was achieved by establishing a functional link between the equipment health index and the failure rate.Furthermore,Alinejad et al.[15] proposed a complete approach to determine the equipment health index by employing the Monte Carlo method.Nevertheless,the aforementioned approaches solely consider the past health condition,neglecting the influence of previous running parameters on the future condition of the equipment.On the other hand,deep learning techniques facilitate the implementation of multilayer nonlinear transformations,hence enhancing the ability to extract and discern correlations among input parameters.The utilisation of layer-by-layer coding networks and Deep Neural Networks (DNN) for condition monitoring of wind turbine main bearings and gearboxes has been proposed in references[16,17].The operational data of condensers comprises multidimensional time series,wherein there exists temporal dependence between the current and historical moment data.The modelling methods mentioned above have limitations due to their reliance on single time point input and their inability to account for time series dependencies.In recent times,there has been an increase in the utilisation of LSTM techniques by researchers for the purpose of forecasting the operational condition of nuclear power plants.In their study,Zhang et al.[18]employed a Long Short-Term Memory(LSTM)neural network to forecast the pressure of steam generators by utilising a fusion technique for multi-sensor signals.The experimental outcomes provided confirmation of the LSTM model’s effectiveness in predicting the operational condition of nuclear power plants.Babu et al.[19]employed a Long Short-Term Memory(LSTM)-based neural network to estimate the state of water health.Similarly,Tuerxun et al.[20]utilised LSTM to predict renewable energy generation in wind farms.However,it should be noted that the training pace of this model is relatively slow while achieving comparable prediction accuracy.In comparison to alternative deep learning approaches,the Informer model is employed for the purpose of extracting time connections between input sequences,resulting in superior accuracy in temporal prediction.In their study,Bommidi et al.[21] put out the utilisation of the Informer model for the prediction of wind speed indicators characterised by intricate uncertainties.Li et al.[22] introduced a novel approach for predicting the Remaining Useful Life(RUL)of rolling bearings.This method combines multiscale degradation metrics with the Informer model,addressing the limitations of traditional prediction models in terms of slow runtime and limited reliance on long time series data.However,it is important to note that the generalisation of this approach still requires improvement,particularly when handling feature parameters with varying probability distributions.In response to the variability in condenser performance caused by external factors,certain scientists have proposed integrating error correction techniques with predictive models as a means to improve the resilience and adaptability of these models.The researchers in the study conducted by Theocharides et al.[23] utilised linear regression as a method to address the prediction bias associated with solar radiation angles.This approach resulted in improved accuracy of day-ahead forecasts for photovoltaic power generation.The aforementioned model possesses the merits of being straightforward and exhibiting a notable degree of computational efficacy.Nevertheless,the performance of this method is influenced by the distribution of data and the underlying assumptions of the model,limiting its ability to capture nonlinear correlations.Zou et al.[24] employed a Markov chain methodology to rectify the prediction outcomes pertaining to China’s Gross Domestic Product (GDP).The present model incorporates the empirical data collected throughout a specific historical timeframe,effectively accounting for the cumulative effects of inaccuracies.This study aimed to examine the progression of condenser health status.This paper presents a novel approach by integrating the Informer model with the MC error correcting method to estimate the health of condensers,taking into account the natural state transition patterns in health evolution.This methodology enables enhanced precision in forecasting the condition of equipment,particularly within the framework of extensive data settings.

    In order to effectively schedule the cleaning time of equipment and optimise the monitoring of condenser health status,as well as accurately predict the equipment’s deterioration trend,this paper proposes two models:a fusion health degree model utilising high-dimensional Mahalanobis distance,and a combined prediction model that incorporates this health degree.The initial step involves utilising a collection of multivariate feature parameters associated with health degree as the input.This input is then utilised to construct a fusion health degree model,employing high-dimensional Mahalanobis distance.The fusion model integrates both subjective Analytic Network Process(ANP)and objective Entropy Weight Method (EWM),and is commonly referred to as ANP-EWM.In this study,we propose the utilisation of a refined Informer model in conjunction with a Markov error correction technique referred to as Generalised Correntropy-Informer-Markov Error Correction(GC-Informer-MEC).The Informer model utilises the generalized correlation entropy as its loss function to enhance the correlation between feature data and the sequential correlation of time-series data.This approach addresses the issue of low prediction accuracy observed in Informer models when handling non-Gaussian probability distribution parameters.Furthermore,the issue of unpredictability in the process of predicting the level of health degree is mitigated through the implementation of enhanced error correcting techniques.This research shows that the condenser’s calculated health degree using the proposed approach is congruent with its real health condition based on the examination of field timeseries data.The provided analysis delivers a quantitative assessment of the equipment’s health status based on its running time.This analysis effectively forecasts the future trend of the condenser and provides a dependable reference for selecting the optimal cleaning time for the condenser.

    2 Building of Weight Fusion Health Degree Model Based on High-Dimensional Mahalanobis Distance

    LetX=[x1,x2,x3,...,xm] be a multidimensional vector composed ofmfeatures,then the Mahalanobis distance between theith step feature vectorXiand the health featureX0can be expressed as[25]:

    whereXi=[xi1,···,xik,···,xim],xikdenotes thekeigenparameter value in theith step of monitoring;X0=[x1,···,xk,···,xm],xkrepresents thekeigenparameter value when the equipment is in full health;Σis the covariance matrix of the multidimensional variablesXiandX0.

    Taking into account the varying degrees of influence that each characteristic parameter within the multivariate characteristic parameters has on the equipment’s health status,we introduce weights to enhance their significance in determining the health degree.Thus,the fusion healthiness model can be formulated as follows:

    whereH(Xi,X0)represents the health degree of stepi,and the health degree range is defined as[0,1];The weight vector,denoted asW,is derived by integrating the analytic network process incorporating subjective experience and the entropy weight method based on objective big data analysis.

    2.1 Calculation of W1 Based on Analytic Network Process

    The Analytic Network Process (ANP) is a hierarchical multi-criteria decision-making method that builds upon the AHP method.ANP takes into account the multi-level structure and the interdependencies among numerous indicators,resulting in an enhanced approach [26].This study focuses on the condenser as the subject of investigation,employing a thermodynamic mechanism analysis.The primary objective is to develop a health evaluation system comprising four distinct characteristic factors.The comprehensive evaluation of the condenser equipment’s performance can be achieved to a certain extent by considering the overall heat transfer coefficient and cleaning coefficient.Condensate subcooling degreeXΔtgl,circulating water temperature riseXΔtxhs,circulating water pump motor currentXIsband end differenceXΔtdcall have certain influence on heat transfer coefficient and cleaning coefficient.Analyzing the working principle of the condenser,it can be obtained that there is a coupling relationship between the four characteristic parameters[27].Establishing a health assessment framework for the condenser requires analyzing the extent to which its characteristic parameters affect overall performance.This involves calculating the weights of different characteristic parameters in the evaluation of equipment health.The use of the analytic network process allows for the comprehensive consideration of the interrelationships among these characteristic parameters.It facilitates the development of a multi-criteria weight calculation method,enabling accurate assessment of the impact of different characteristic parameters on equipment health degree.The diagram in Fig.1 illustrates the hierarchical structure of the algorithm applied in the analytic network method.

    Figure 1:The hierarchical structure of analytic network process

    Based on the above hierarchical structure,with the control layer as the main criterion and the four characteristic parameters in the network layer as the secondary criterion,the elements are compared in pairs to build an initial unweighted supermatrixWs:

    where the column vector ofWijrepresents the indirect dominance degree comparison of the influence of the elements in element groupion the elements in element groupjwith the elements in element groupjas sub-criteria:

    Taking the overall heat transfer coefficient and the overall cleaning coefficient as the secondary criterion,a judgment matrix is established by pairwise comparison,and then the weight supermatrixAsis obtained:

    On the basis of formulas(3) and (5),the weighted supermatrixWasshown in formula(6) is constructed:

    Finally,the weight is obtained by calculating the limit weighted supermatrix,and the priority of the feature variable is determined:

    whereW1is the weight vector calculated by the ANP.

    2.2 Calculation of W2 Based on Entropy Weight Method

    According to the above four characteristic parameters and the data set ofngroups of samples,an evaluation matrixM=is built.In instances when the initial evaluation matrix comprises values that are equivalent to zero,it indicates that these characteristic parameters do not contribute to any loss of information in the health assessment.As a result,the entropy values associated with these parameters reach a state of zero,hence causing a proportional augmentation in their respective weights.This phenomenon has the potential to induce bias into the outcomes of health assessments.In real-world scenarios,feature parameters may face instances of zero values as a result of measurement mistakes and various other factors.In order to mitigate any bias in the weights,this study proposes the incorporation of a little perturbation when employing the Entropy Weight Method (EWM) to calculate weights,particularly in cases when zero values are present.The purpose of this perturbation is to reduce the potential influence of zero values on the process of determining weights.Subsequently,mijis standardized to obtain a standard evaluation matrixV=[28].

    The entropy valueejof thejth feature parameter can be expressed as:

    wherehij=i=1,2,···,n;j=1,2,···,4 denotes the weight of theith sample value under thejth characteristic parameter.

    Further,the weight vector calculated by the entropy weighting method is obtained as[29]:

    whereW2is the weight vector calculated by entropy weight method,ωj==1,2,···,4.

    2.3 Health Degree Model of Weight Fusion

    To enhance the precision of assessing the health degree,the health degree model is constructed by incorporating the high-dimensional Mahalanobis distance.Additionally,a fusion weighting method is employed,which combines both subjective and objective weighting methods,to determine the weights assigned to the characteristic parameters that represent the health degree.This methodology circumvents the inherent subjectivity associated with the ANP and addresses the computational bias resulting from the limited data available in the EWM.Wis a linear combination ofW1andW2.Due to the subjective sources of representing health degrees using the ANP and the objective sources of representing health degrees using the EWM,this study considers both subjective and objective sources,each contributing equally to the composite weight,to characterize equipment health degree.The proposed model for determining the fusion weight of health degree is as follows:

    whereαandβare the proportions of the weights calculated by the ANP and the EWM in the comprehensive weight,respectively.

    3 Combined Prediction Model Based on GC-Informer and Markov Error Correction

    Various aspects in the operational environment of nuclear power plants have an impact on the operating parameters of the condenser.The equipment parameter data in this context generally follows a non-Gaussian distribution.The health degree time-series data is subject to the influence of the equipment operating environment and many characteristics,resulting in a non-Gaussian distribution.In order to further understand the deterioration trend of the condenser and improve the stability and accuracy of the health degree prediction model,aiming at the problem that the Informer model has low prediction accuracy for non-Gaussian distribution parameters,an improved combination prediction model of Informer and MEC is established.By changing model loss function and increasing the error correction amount to achieve the purpose of accurately predicting the health degree of the condenser.This model comprehensively takes into account the impact of historical running data on health status.It utilizes various time series as input variables in order to forecast the health degree of the condenser.

    3.1 GC-Informer Model

    The GC-Informer prediction model obtains the temporal dependence of the five-dimensional feature parameter set including the historical health degree of the condenser through the encoder,and generates the health prediction result through the decoder.Considering the significance of time in health degree prediction,before inputting the model,time stamp encoding is used to encode the year,month,day,hour,minute and other information in the data to improve the predictive ability of the model in long-term series.The operation is as follows:

    3.1.1Encoder

    The encoder is mainly composed of multi-head probsparse self-attention mechanism modules and distillation modules.It takes the historical running data Xt=of the condenser as input and uses the multi-head sparse self-attention mechanism to extract long-term dependencies in the temporal data.The attention distillation mechanism is then applied to reduce the temporal dimension by half for each individual layer’s features.The feature dimension is further compressed and essential information is extracted using a 1D convolutional neural network[30].The distillation process fromjlayer toj+1 layer is illustrated by Eq.(14),enhancing the model’s capability to handle long time series data.

    Figure 2:Encoder structure

    3.1.2Multi-HeadProbsparseSelf-AttentionMechanism

    The self-attention mechanism assigns weight coefficients by calculating the similarity between the Query and Key[31]:

    whereQ∈RLQ×d,K∈RLK×d,V∈RLν×dis the feature matrix obtained by linear transformation of the input feature variables,dis the input dimension,and Softmax(·)is the activation function.To further discuss the self-attention mechanism,let theith row ofQ,K,Vbeqi,ki,vi,respectively.The attention coefficient of theith row of Query is defined as a kernel smoother in the probabilistic form:

    where the first term is the Log-Sum-Exp(LSE)ofqion all Keys,and the second term is their arithmetic mean.If theith row of Query obtains a largerMthe dot product has a larger contribution to attention.The probsparse self-attention mechanism is implemented by allowing each Key to focus only on the Query with higher importance,thus reducing the model computational complexity.

    3.1.3Decoder

    The input to the decoder is a five-dimensional data sequence,including the historical health index of the condenser,along with a segment of zeros that is equal to the prediction step length.Zeros are used as placeholders for the predicted values to prevent the model from accessing future sequence information during training.The input vector is represented as shown in Eq.(19).

    Within the decoder,the input time series data is first calculated by the masked multi-head probsparse self-attention mechanism module,and then the sparse self-attention operation is performed with the feature map output by the encoder.The resulting computation is then passed through a fully connected layer to obtain the prediction resultY=[yt+1,yt+2,···,yt+m,···,yt+n].The prediction results are compared with the actual values to calculate the loss function,which is used to continuously optimize the model based on the errors.In order to address the probability distribution characteristics of the running parameters in the subject of this study,the generalized correntropy function is applied as the loss function for training the Informer prediction model.The mathematical formula for the GC is shown in Eq.(20):

    whereyr+mis the predicted value,Hr+msamples the true value,Nis the number of data,α>0 is the shape parameter of the generalized Gaussian distribution,β=σis the Gaussian kernel parameter,andΓ(·) is the gamma function.Input the test data set into the trained GC-Informer model,and output the health degree prediction result.Its overall structure is shown in Fig.3.

    Figure 3:Informer model structure

    3.2 Markov Error Correction

    Considering the presence of randomness in the prediction process and aiming to enhance the prediction accuracy of condenser health degree,a Markov chain is applied to partition the development pattern and potential outcomes of condenser health into different states.The transition rules between different states are studied to modify the predicted results of condenser health degree.Characterized by the transition matrix,in the process of solving the state probability,the states are divided intontypes(namelyS1,S2,···,Sn)through the analysis of the condenser health degree prediction error sequence,and the state division is based on the following[32]:

    where the EVEN(·)function is the nearest even number taken along the direction of increasing absolute value,εis the sequence of residuals,anddis the difference between the maximum and minimum residuals.Next,based on the statistics of health prediction errors,the transition probabilities between states are calculated.Specifically,the probability of relative error transitioning from stateSito stateSjoverktime steps is given by:

    whereMij(k)represents the number of times stateSitransitions to stateSjafterktime steps,andMirepresents the number of occurrences of relative error in stateSi.The matrix form is:

    Assuming that the current state isSi,in theith row of the transition probability matrix,Pijis the maximum value in this row,then the possibility of the stateSjappearing at the next moment is the greatest,and it is taken as the corresponding prediction result.

    Finally,an error correction valueet+mis taken for each state to modify it,resulting in the predicted health degree.The correction basis is as follows:

    whereΔmaxandΔminare the upper and lower limits of the relative error state interval,respectively,andyt+mis the predicted health degree of GC-Informer model.

    3.3 Establishment of GC-Informer-MEC Combined Prediction Model

    Based on the above theories,this paper combines GC-Informer and MEC to establish a combined predicting model that improves Informer and Markov error correction.The input of the prediction model is Xt=The inputs of the model are presented in Table 1.The historical time-series data of the condenser’s characteristic parameters include the condensate subcooling degreeXΔtgl,circulating water temperature risecirculating water pump motor currentand end differenceAdditionally,the historical time-series data of the condenser’s health degree is represented by the health degree.One instance of variable data at a particular moment(row)can be observed as[0.4977,7.1874,0.5846,6.7075,0.9571].The time-series data,which consists of five dimensions,is divided into segments for the composite prediction model.The segmentation is achieved by utilising a sliding time-series window with a length of 200.In every iteration of the model,a matrix with dimensions of 200*5 is employed as input data,which is segmented using the sliding window technique.

    Table 1: The inputs of the GC-Informer-MEC combined prediction model

    In order to enhance the reliability of predicting outcomes,the data utilised for the model is divided into three distinct sets:training datasets,validation datasets,and testing datasets.These sets are allocated in a ratio of 6:2:2,respectively.The training datasetand validation datasetare used to train the model,and the testing datasetis fed into the trained GC-Informer model for health degree prediction.The first part of the proposed combined prediction model is the GC-Informer prediction model.This model compares the Informer method’s predicted results against actual values using the training and validation datasets,calculating the GC loss function.Continuous optimization of the model is achieved through error minimization,resulting in a well-trained GC-Informer model.Using the testing dataset,the trained GC-Informer model generates predictions.The latter half of the combined prediction model applys the Markov Error Correction algorithm to calculate errors between the GC-Informer predicted results and actual values.It subdivides the error state space and calculates the state transition matrix,obtaining error correction values.The final prediction of the model is a combination of GC-Informer predictions and MC correction values.The overall process is illustrated in Fig.4.

    4 Example Verification

    This paper utilizes an example dataset comprising 12,960 sets of field operation data from January to March 2022 of a domestic 1000 MW nuclear power plant unit for validation and analysis.Using four characteristic parameters,the health degreeH(X,X0)of the condenser is determined.In addition,the computed condenser health degree is utilised in conjunction with these four distinctive factors to forecast the trajectory of changes in the health degree of the condenser.In order to improve the quality of the data,several procedures were adopted during the processes of data collecting and conversion.The missing data were resolved by employing the mean fill approach,and the outliers were managed with the plural fill method.

    Figure 4:GC-Informer-MEC model

    4.1 Calculation of Condenser Health Degree

    Through the hierarchical analysis of the fusion health degree model,the overall heat transfer coefficient and the overall cleanliness coefficient are taken as the determining factors,and the initial unweighted supermatrix for the control layer health degree indicator is obtained as follows:

    Using the overall heat transfer coefficient being superior to the overall cleanliness coefficient as a sub-criterion,the weight supermatrix is obtained as:

    The weighted supermatrixWasis obtained according to Eq.(6):

    Ultimately,the weight vector of the feature parameters based on the analytic network process according to Eq.(7)is:

    The weight vectorW2of the feature parameters calculated based on the entropy weight method is:

    Weighting was performed using the fusion weighting method based on Eq.(11).The subjective AHP and objective EWM based weighting methods contributed to proportionsαandβ,respectively,in the comprehensive weight,with both being 0.5.The final weights of each feature parameter are presented in Table 2.The weights of the health feature parameters were substituted into the fusion health degree model Eq.(11),resulting in the health degree values of the condenser at each time step.

    4.2 Validation of Condenser Health Degree Trend Prediction

    On the basis of the fusion health degree model,the health degree of the condenser is predicted in multiple steps with 12960 sets of data in the example data set,and then the deterioration evolution trend of the condenser can be judged.Next,based on Eqs.(21) and (22),the training set data were used to partition the error state space.By considering the relative error,the error sequence was divided into four state intervals,which are defined as follows:S1(-0.059,-0.005),S2(-0.005,-0.001),S3(-0.001,0.005),S4(0.005,0.047).The state transition matrixPis obtained by counting the state transition as follows:

    Table 2: Feature weight coefficient

    In order to demonstrate the improvement in the accuracy of condenser health degree prediction achieved by the multidimensional inputs GC-Informer-MEC combination prediction model proposed in this paper,predictions were made using only one-dimensional time-series data consisting of historical health information as the model input.These health prediction results were then compared with the approach presented in this paper,as depicted in Fig.5.The image depicts three curves:green,blue,and orange.These curves correspond to the actual health degree,health degree forecasts derived from five-dimensional time-series data as input,and health degree predictions derived from onedimensional time-series data as input,respectively.The graphic illustrates three vertical dashed lines,where the black line represents the commencement of the prediction,and the two red lines represent the anticipated condenser cleaning times for the multi-input and single-input scenarios,respectively.

    Figure 5:GC-Informer-MEC model prediction results

    Upon careful examination of Fig.5,it becomes apparent that the utilisation of both historical running parameters and historical health degree as inputs to the model regularly results in projected health degrees that closely coincide with the genuine degree.This alignment demonstrates a notable level of accuracy and few discrepancies.In contrast,employing a solitary metric of previous health degree for the purpose of predicting future health outcomes results in diminished levels of accuracy.In this study,the determination of the optimal time point for condenser cleaning is based on the application of a health degree threshold of 0.8.The single-input prediction model exhibits a delay of 508 sample points,each representing a 10-minute interval,compared to the multi-input prediction model in predicting the cleaning time point.As a result,the single-input prediction model indicates a cleaning time point that is approximately 84 h later than that of the multi-input prediction model.The delay has the potential to have a negative impact on the operational quality of the equipment,resulting in extended operation of the condenser at a reduced efficiency level.This,in turn,can lead to a decline in the economic efficiency of the secondary loop.When the condenser operates below a health degree of 0.8,performing timely cleaning can increase its efficiency by 1%.Based on the thermal calculation of a 1000 MW unit with a coal consumption of 320 grams per kWh,it can be determined that the impact of this timely condenser cleaning on coal consumption is:

    When calculated at 460 RMB per metric ton of standard coal,utilizing the multi-input prediction model for timely cleaning of the condenser equipment can result in a cost saving of 123,648 RMB.

    The prediction errors are presented in Table 3.The analysis table illustrates that incorporating numerous data inputs significantly improves the predictive accuracy of the model.Furthermore,it indicates that the selected historical running parameters of the condenser have a certain impact on the forecast of the health degree.When solely utilising health degree time series data as the primary input for the prediction model,there is a notable discrepancy in the projected health outcomes,despite the approximate alignment of the predicted trend with the actual health progression.This approach provides an impartial reflection of the anticipated future evolution trend of condenser equipment.The predictive model demonstrates strong performance in accurately projecting future trends for a single time-series.In situations when exactitude is not the primary concern,a solitary indicator may be employed to forecast the forthcoming condition of the apparatus.In scenarios where the precise prediction of equipment condition is of utmost importance,the construction of multi-dimensional time series allows for more accurate forecasts of the health level.This,in turn,aids operators in effectively monitoring the health status of the condenser equipment and facilitates the development of a well-informed cleaning schedule for the condenser.

    Table 3: Predictive errors under multiple indicators

    In order to assess the enhanced efficacy of the GC-Informer-MEC prediction model in comparison to utilising a solitary model,three commonly applied prediction models,specifically Informer,Long Short-Term Memory (LSTM),and Deep Neural Network (DNN),were chosen as control experimental groups for the purpose of comparing their respective prediction performances.The sample data for this study consists of 12,960 sets of field operating data obtained from a 1000 MW unit of a nuclear power plant.To create a training set,the historical data from the prior 7,800 sets at the projected time were picked.The comparative analysis of the projected outcomes derived from the utilisation of multi-source time-series data is presented in Fig.6.

    Figure 6:Comparison of prediction results of different models

    Upon examination of Fig.6,it becomes apparent that the chosen prediction models,namely Informer,LSTM,and DNN,demonstrate differing levels of accuracy in forecasting the future trend of the health degree.These models closely align with the progression of the actual health degree.Nevertheless,the efficacy of forecasting future health outcomes varies based on the selected predictive model.The present study highlights the effectiveness and accuracy of the suggested GC-Informer-MEC combined prediction model,as evidenced by the results presented in Fig.6.Table 4 displays the discrepancies observed in the forecast outcomes across several prediction models.

    Table 4: Error results from multiple predictive models

    5 Conclusion

    Aiming at the problem of determining the cleaning time for condensers,this paper proposes a fusion health degree model based on improved Mahalanobis distance.On the basis of the fusion health degree model,a combined prediction model based on GC-Informer-MEC is established to predict the health degree of the condenser,and then the state trend of the deterioration and evolution of the condenser can be judged.Applying the above method to an example,the following conclusions are obtained:

    The improvement of the Mahalanobis distance is achieved by the utilisation of a weighting method that combines the analytic network process and the entropy weight method.This approach allows for the integration of expert knowledge and objective observations,enabling a quantitative assessment of the actual status of the condenser’s health.The GC-Informer-MC combination prediction model is constructed based on the health degree,utilising past running parameters and historical health degree data.In contrast to a single-input model,the combined model has the capability to anticipate condenser deterioration durations at 508 sample locations,hence facilitating prompt equipment cleaning and resulting in an estimated cost savings of roughly 123,648 RMB.In contrast to a singular prediction model,the utilisation of this approach yields a more refined forecast and study of the degradation pattern exhibited by the condenser.Consequently,it furnishes a more precise benchmark for ascertaining the optimal timing for equipment cleaning.

    Health indicators,which are derived from the measurement of health status,are extensively utilised in various industrial processes.However,further enhancements are needed in order to get a more precise depiction of equipment status through the utilisation of health degree indicators.In our future research,we aim to improve the health degree model by employing optimisation methods.Additionally,we plan to implement online weight updates for the subjective and objective source weights,denoted asαandβrespectively,which are used to quantify the health degree.This approach is expected to enhance the accuracy of equipment status quantification.

    Acknowledgement:The author greatly appreciates the editor’s encouragement and the anonymous reviewers’suggestions to improve our paper.

    Funding Statement:This work was supported by the National Natural Science Foundation of China(51906133).

    Author Contributions:The authors confirm contribution to the paper as follows: study conception and design:Hong Qian;data collection:Guangji Wang,Qingyun Yan;analysis and interpretation of results:Haixin Wang;draft manuscript preparation:Haixin Wang,Hong Qian.All authors reviewed the results and approved the final version of the manuscript.

    Availability of Data and Materials:The data used this study was obtained from the actual operation records of the condenser at Fangchenggang Nuclear Power Plant from January to March 2022.Due to safety and confidentiality requirements of the nuclear power plant,this data cannot be released.

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

    日日干狠狠操夜夜爽| 久久精品国产自在天天线| 亚洲欧美精品专区久久| 床上黄色一级片| 免费看美女性在线毛片视频| 毛片女人毛片| 十八禁国产超污无遮挡网站| 在线天堂最新版资源| 国产亚洲5aaaaa淫片| 国产黄a三级三级三级人| 在线观看免费高清a一片| 51国产日韩欧美| 国产大屁股一区二区在线视频| 人人妻人人看人人澡| 国产老妇女一区| 26uuu在线亚洲综合色| 欧美一区二区亚洲| 99热全是精品| 色5月婷婷丁香| 国产永久视频网站| av一本久久久久| 在线 av 中文字幕| 久久精品熟女亚洲av麻豆精品 | av国产免费在线观看| 一级a做视频免费观看| 国语对白做爰xxxⅹ性视频网站| 久热久热在线精品观看| 人妻制服诱惑在线中文字幕| 日韩中字成人| 免费观看无遮挡的男女| 国产久久久一区二区三区| 亚洲国产日韩欧美精品在线观看| 成人毛片a级毛片在线播放| 99久久精品热视频| 青春草亚洲视频在线观看| 最近2019中文字幕mv第一页| 亚洲欧美成人精品一区二区| 日本色播在线视频| 欧美成人a在线观看| 熟妇人妻不卡中文字幕| 国产有黄有色有爽视频| 国产国拍精品亚洲av在线观看| 免费观看性生交大片5| 日本黄色片子视频| 国产成人精品福利久久| 国产成人91sexporn| 美女内射精品一级片tv| 久久久久久久亚洲中文字幕| 国产精品综合久久久久久久免费| 天堂√8在线中文| 国产亚洲av嫩草精品影院| 欧美日韩亚洲高清精品| 岛国毛片在线播放| 亚洲精品日韩av片在线观看| 亚洲国产精品sss在线观看| 国精品久久久久久国模美| 国产精品一区二区三区四区久久| 有码 亚洲区| 69人妻影院| 一本久久精品| 青春草国产在线视频| 久久精品国产鲁丝片午夜精品| 欧美一级a爱片免费观看看| 高清视频免费观看一区二区 | 最近最新中文字幕免费大全7| 日本猛色少妇xxxxx猛交久久| 国产午夜精品久久久久久一区二区三区| 欧美+日韩+精品| 日韩一区二区视频免费看| 成人亚洲欧美一区二区av| 久久久久久久久大av| 成人无遮挡网站| 亚洲天堂国产精品一区在线| 国产成人freesex在线| 国产精品无大码| 天美传媒精品一区二区| 午夜福利成人在线免费观看| 国内精品一区二区在线观看| 2022亚洲国产成人精品| 国精品久久久久久国模美| 晚上一个人看的免费电影| 久久鲁丝午夜福利片| 国产不卡一卡二| 69av精品久久久久久| 精品久久久久久成人av| 精品不卡国产一区二区三区| 国产成人免费观看mmmm| 久久国内精品自在自线图片| 欧美变态另类bdsm刘玥| 亚洲欧美精品专区久久| 丰满少妇做爰视频| 国产一区有黄有色的免费视频 | 校园人妻丝袜中文字幕| 别揉我奶头 嗯啊视频| 三级男女做爰猛烈吃奶摸视频| 免费播放大片免费观看视频在线观看| 国产av国产精品国产| 内地一区二区视频在线| 日本-黄色视频高清免费观看| av免费在线看不卡| 午夜爱爱视频在线播放| 欧美激情国产日韩精品一区| 丰满人妻一区二区三区视频av| 午夜福利在线观看吧| 国产大屁股一区二区在线视频| av福利片在线观看| 欧美精品国产亚洲| 51国产日韩欧美| 又爽又黄a免费视频| 国产一区二区在线观看日韩| 欧美激情国产日韩精品一区| 欧美日本视频| 精品不卡国产一区二区三区| 丰满人妻一区二区三区视频av| av黄色大香蕉| 亚洲精品自拍成人| 偷拍熟女少妇极品色| 蜜桃亚洲精品一区二区三区| 你懂的网址亚洲精品在线观看| 26uuu在线亚洲综合色| 能在线免费看毛片的网站| 亚洲av成人精品一区久久| 精品一区二区三区人妻视频| 亚洲av一区综合| 特大巨黑吊av在线直播| 国产伦在线观看视频一区| videos熟女内射| 精品少妇黑人巨大在线播放| 97超碰精品成人国产| 97超碰精品成人国产| 成人高潮视频无遮挡免费网站| 边亲边吃奶的免费视频| 国精品久久久久久国模美| videos熟女内射| 国产一级毛片七仙女欲春2| 我的女老师完整版在线观看| 日韩 亚洲 欧美在线| kizo精华| 国内揄拍国产精品人妻在线| 18禁裸乳无遮挡免费网站照片| 直男gayav资源| 夜夜看夜夜爽夜夜摸| www.色视频.com| 日韩亚洲欧美综合| 久久久欧美国产精品| 亚洲精品,欧美精品| 免费观看性生交大片5| 日韩欧美国产在线观看| 一区二区三区高清视频在线| 天美传媒精品一区二区| 中文字幕久久专区| 亚洲真实伦在线观看| 超碰av人人做人人爽久久| 高清毛片免费看| 亚洲第一区二区三区不卡| 国产亚洲午夜精品一区二区久久 | 又爽又黄a免费视频| 国产精品综合久久久久久久免费| 亚洲av一区综合| 尤物成人国产欧美一区二区三区| 成人亚洲欧美一区二区av| 日本猛色少妇xxxxx猛交久久| 国产黄色视频一区二区在线观看| 99re6热这里在线精品视频| 国产成人91sexporn| 成人性生交大片免费视频hd| 国产精品一区二区三区四区免费观看| 亚洲国产欧美在线一区| 欧美丝袜亚洲另类| av福利片在线观看| 国产精品爽爽va在线观看网站| 青春草国产在线视频| 免费av不卡在线播放| 国产av不卡久久| 亚洲最大成人手机在线| www.av在线官网国产| 丰满人妻一区二区三区视频av| 国产精品嫩草影院av在线观看| 天堂av国产一区二区熟女人妻| 久久99精品国语久久久| 美女cb高潮喷水在线观看| 80岁老熟妇乱子伦牲交| 日韩制服骚丝袜av| 51国产日韩欧美| 国产精品.久久久| 五月玫瑰六月丁香| 女人被狂操c到高潮| 三级男女做爰猛烈吃奶摸视频| 国产精品一区二区三区四区久久| 欧美丝袜亚洲另类| 免费观看av网站的网址| 中国国产av一级| 亚洲av在线观看美女高潮| 国产高潮美女av| 最近中文字幕高清免费大全6| 国产v大片淫在线免费观看| 网址你懂的国产日韩在线| 一级片'在线观看视频| 国产一级毛片在线| 久久久亚洲精品成人影院| av在线播放精品| 国产精品麻豆人妻色哟哟久久 | 亚洲av不卡在线观看| 精品一区二区三区人妻视频| 女人十人毛片免费观看3o分钟| 高清毛片免费看| 男女边摸边吃奶| 午夜福利视频精品| 尾随美女入室| 午夜激情久久久久久久| 欧美高清性xxxxhd video| 精品人妻偷拍中文字幕| 国产精品久久视频播放| 亚洲国产日韩欧美精品在线观看| 国产成人91sexporn| 久久久久性生活片| 熟女电影av网| 成人漫画全彩无遮挡| 高清在线视频一区二区三区| 免费看a级黄色片| 国模一区二区三区四区视频| 一个人看视频在线观看www免费| 99久久中文字幕三级久久日本| 人妻少妇偷人精品九色| 成年女人看的毛片在线观看| 午夜福利视频精品| 不卡视频在线观看欧美| 国产精品爽爽va在线观看网站| 蜜桃久久精品国产亚洲av| 国产一级毛片七仙女欲春2| 国产精品精品国产色婷婷| 亚洲美女搞黄在线观看| 2022亚洲国产成人精品| 男人和女人高潮做爰伦理| 国产伦理片在线播放av一区| 国产欧美另类精品又又久久亚洲欧美| 国产中年淑女户外野战色| av福利片在线观看| 欧美日韩视频高清一区二区三区二| 久热久热在线精品观看| 国产亚洲5aaaaa淫片| 一本久久精品| 国产精品一区二区在线观看99 | 丝袜喷水一区| 搡女人真爽免费视频火全软件| 2021少妇久久久久久久久久久| 欧美高清成人免费视频www| 91av网一区二区| 69av精品久久久久久| 欧美性猛交╳xxx乱大交人| 久久久久久久久久久免费av| 久久精品人妻少妇| 五月玫瑰六月丁香| 高清av免费在线| 天堂√8在线中文| 久久99精品国语久久久| 最近视频中文字幕2019在线8| 成人性生交大片免费视频hd| 久久久欧美国产精品| 国产精品无大码| 午夜福利视频精品| 中文在线观看免费www的网站| 成年版毛片免费区| 亚洲美女视频黄频| 国精品久久久久久国模美| 久久精品国产鲁丝片午夜精品| 色网站视频免费| 亚洲精品色激情综合| 欧美zozozo另类| 国产精品蜜桃在线观看| 男人舔奶头视频| 97热精品久久久久久| av在线亚洲专区| 免费电影在线观看免费观看| 免费看av在线观看网站| 久久热精品热| 五月天丁香电影| 国产男女超爽视频在线观看| 日本av手机在线免费观看| 波多野结衣巨乳人妻| 久久精品久久久久久噜噜老黄| 五月天丁香电影| 九草在线视频观看| 亚洲18禁久久av| 亚洲内射少妇av| 国产又色又爽无遮挡免| 成年版毛片免费区| 五月伊人婷婷丁香| 国产精品女同一区二区软件| 蜜桃亚洲精品一区二区三区| 中文资源天堂在线| 中文天堂在线官网| 国产成年人精品一区二区| 99热6这里只有精品| 久久久久国产网址| 国产男人的电影天堂91| 欧美性感艳星| 国产麻豆成人av免费视频| 国产乱人偷精品视频| 日本-黄色视频高清免费观看| 日韩欧美三级三区| 午夜爱爱视频在线播放| 色5月婷婷丁香| 亚洲精品成人久久久久久| 国产男人的电影天堂91| 国产精品久久久久久久久免| 别揉我奶头 嗯啊视频| 日韩人妻高清精品专区| 日韩av在线大香蕉| 乱系列少妇在线播放| 免费黄频网站在线观看国产| 国产黄片美女视频| 亚洲欧美日韩无卡精品| 一级毛片 在线播放| 久久国内精品自在自线图片| 国产精品人妻久久久久久| 婷婷色综合大香蕉| 亚洲国产精品成人综合色| 别揉我奶头 嗯啊视频| 国产高潮美女av| 日日摸夜夜添夜夜爱| 免费少妇av软件| 神马国产精品三级电影在线观看| 国产精品人妻久久久影院| 精品一区二区三区视频在线| 国产伦理片在线播放av一区| 免费播放大片免费观看视频在线观看| 免费大片18禁| 十八禁国产超污无遮挡网站| 26uuu在线亚洲综合色| 日韩 亚洲 欧美在线| 18+在线观看网站| 少妇丰满av| 日本猛色少妇xxxxx猛交久久| 欧美激情久久久久久爽电影| 国产一级毛片七仙女欲春2| 亚洲欧美日韩无卡精品| 亚洲欧洲日产国产| 晚上一个人看的免费电影| 尾随美女入室| 蜜臀久久99精品久久宅男| 亚洲国产精品成人综合色| 在线免费观看不下载黄p国产| 黄色配什么色好看| 又粗又硬又长又爽又黄的视频| 国产老妇伦熟女老妇高清| 亚洲精华国产精华液的使用体验| 国产精品无大码| 日本欧美国产在线视频| 久久久久免费精品人妻一区二区| 精品人妻视频免费看| 久久久久九九精品影院| 午夜免费观看性视频| 精品久久久久久久久亚洲| 亚洲av免费高清在线观看| 99re6热这里在线精品视频| 亚州av有码| 欧美丝袜亚洲另类| 久久久久久久久久人人人人人人| 黑人高潮一二区| 啦啦啦啦在线视频资源| 国产有黄有色有爽视频| 大香蕉97超碰在线| 成人毛片60女人毛片免费| 亚洲av免费高清在线观看| 国产成人精品久久久久久| 亚洲av免费高清在线观看| 欧美xxⅹ黑人| 久久久久网色| 22中文网久久字幕| 在线免费观看不下载黄p国产| 两个人的视频大全免费| 男人爽女人下面视频在线观看| 全区人妻精品视频| 最后的刺客免费高清国语| 国产熟女欧美一区二区| 亚洲精品日本国产第一区| 国产女主播在线喷水免费视频网站 | 啦啦啦中文免费视频观看日本| kizo精华| 久久99热这里只频精品6学生| 免费看不卡的av| 嫩草影院精品99| 国产精品麻豆人妻色哟哟久久 | 少妇的逼水好多| 97在线视频观看| 欧美高清成人免费视频www| 国产国拍精品亚洲av在线观看| 亚洲人成网站高清观看| 美女国产视频在线观看| 黄色欧美视频在线观看| 麻豆国产97在线/欧美| 午夜视频国产福利| 在线免费观看的www视频| 99久国产av精品| 97精品久久久久久久久久精品| 中文乱码字字幕精品一区二区三区 | 亚洲人成网站在线播| 亚洲乱码一区二区免费版| 国产成人一区二区在线| 天堂√8在线中文| 亚洲av电影不卡..在线观看| 久久热精品热| 久久国内精品自在自线图片| 亚洲色图av天堂| 国产av不卡久久| 国产亚洲一区二区精品| 亚洲欧美一区二区三区黑人 | 一边亲一边摸免费视频| 久久精品国产亚洲av天美| 亚洲真实伦在线观看| 一个人观看的视频www高清免费观看| 可以在线观看毛片的网站| 26uuu在线亚洲综合色| 免费在线观看成人毛片| 人体艺术视频欧美日本| 免费观看无遮挡的男女| .国产精品久久| 亚洲欧美清纯卡通| 亚洲精品aⅴ在线观看| 国产精品一区www在线观看| 国产精品综合久久久久久久免费| 国产精品嫩草影院av在线观看| 大话2 男鬼变身卡| 久久久久性生活片| 亚洲精品日韩av片在线观看| 亚洲国产精品国产精品| 2018国产大陆天天弄谢| 亚洲精华国产精华液的使用体验| 青青草视频在线视频观看| 国产毛片a区久久久久| 亚洲国产最新在线播放| 国产免费福利视频在线观看| 99热6这里只有精品| 国产av在哪里看| 韩国高清视频一区二区三区| eeuss影院久久| 国产久久久一区二区三区| www.色视频.com| 午夜精品一区二区三区免费看| 一级片'在线观看视频| 国产又色又爽无遮挡免| 免费观看无遮挡的男女| 亚洲第一区二区三区不卡| 精品人妻熟女av久视频| 欧美激情在线99| 最近最新中文字幕大全电影3| 日韩欧美三级三区| 午夜免费男女啪啪视频观看| 成人午夜高清在线视频| 国产白丝娇喘喷水9色精品| 午夜老司机福利剧场| 国产一级毛片七仙女欲春2| 亚洲av免费高清在线观看| 男女啪啪激烈高潮av片| 成年免费大片在线观看| 街头女战士在线观看网站| 99久国产av精品| 99热网站在线观看| 午夜视频国产福利| 久久6这里有精品| 亚洲成人av在线免费| 日韩亚洲欧美综合| 黄片无遮挡物在线观看| 亚洲美女搞黄在线观看| 2018国产大陆天天弄谢| 久久久久免费精品人妻一区二区| 校园人妻丝袜中文字幕| 国产精品久久久久久久电影| 99热6这里只有精品| 亚洲真实伦在线观看| 久热久热在线精品观看| 一级毛片 在线播放| a级一级毛片免费在线观看| 青青草视频在线视频观看| 人人妻人人澡人人爽人人夜夜 | 欧美不卡视频在线免费观看| 久久久久九九精品影院| 天天躁夜夜躁狠狠久久av| 国产一级毛片在线| 亚州av有码| 亚洲欧美成人精品一区二区| 六月丁香七月| 中文字幕av成人在线电影| 国产单亲对白刺激| 七月丁香在线播放| 亚洲欧美成人精品一区二区| 精品久久久久久久末码| 两个人的视频大全免费| 国产亚洲一区二区精品| 国产精品一及| 午夜免费观看性视频| 午夜亚洲福利在线播放| 亚洲第一区二区三区不卡| 国产乱来视频区| 日本wwww免费看| 亚洲av成人精品一区久久| 五月玫瑰六月丁香| 免费看不卡的av| 男女视频在线观看网站免费| 国内少妇人妻偷人精品xxx网站| 少妇熟女欧美另类| 菩萨蛮人人尽说江南好唐韦庄| 欧美一级a爱片免费观看看| 欧美一区二区亚洲| 一级毛片黄色毛片免费观看视频| 51国产日韩欧美| 男女下面进入的视频免费午夜| 高清毛片免费看| 精品久久国产蜜桃| 白带黄色成豆腐渣| 好男人在线观看高清免费视频| 国产精品一区二区性色av| 国产视频内射| 菩萨蛮人人尽说江南好唐韦庄| 精品一区二区三卡| 两个人视频免费观看高清| 一级毛片黄色毛片免费观看视频| 午夜激情欧美在线| 国产精品久久久久久精品电影| 成年女人在线观看亚洲视频 | 精品人妻熟女av久视频| 日韩精品有码人妻一区| a级一级毛片免费在线观看| 一二三四中文在线观看免费高清| 亚洲欧美精品自产自拍| 亚洲精品456在线播放app| 国产成年人精品一区二区| 日韩在线高清观看一区二区三区| av卡一久久| 久久精品人妻少妇| 亚洲人成网站在线播| 有码 亚洲区| 免费无遮挡裸体视频| 国产亚洲av嫩草精品影院| 亚洲精品久久久久久婷婷小说| 国产精品蜜桃在线观看| 一边亲一边摸免费视频| 伦理电影大哥的女人| 国产片特级美女逼逼视频| 国产一区二区在线观看日韩| 日韩一区二区视频免费看| 在线观看免费高清a一片| 中文字幕av成人在线电影| 别揉我奶头 嗯啊视频| 国产亚洲最大av| 人人妻人人澡人人爽人人夜夜 | 免费观看av网站的网址| 亚洲精品456在线播放app| 又粗又硬又长又爽又黄的视频| 伦精品一区二区三区| 亚洲最大成人手机在线| 免费观看a级毛片全部| 成人综合一区亚洲| 国产高清有码在线观看视频| 永久网站在线| 国产一级毛片在线| 一级a做视频免费观看| 国产精品1区2区在线观看.| 水蜜桃什么品种好| 亚洲美女视频黄频| 黄色欧美视频在线观看| 亚洲国产精品成人综合色| 纵有疾风起免费观看全集完整版 | 身体一侧抽搐| 精品99又大又爽又粗少妇毛片| 五月玫瑰六月丁香| 亚洲在久久综合| 亚洲国产色片| 久久这里只有精品中国| 欧美变态另类bdsm刘玥| 久久久久九九精品影院| 自拍偷自拍亚洲精品老妇| 真实男女啪啪啪动态图| 久久97久久精品| 成人毛片a级毛片在线播放| 国产在线一区二区三区精| 欧美性感艳星| 国产免费一级a男人的天堂| 蜜桃久久精品国产亚洲av| 国产亚洲午夜精品一区二区久久 | 亚洲国产精品专区欧美| 欧美+日韩+精品| 久久久久精品性色| 亚洲国产精品专区欧美| 欧美激情国产日韩精品一区| 简卡轻食公司| 七月丁香在线播放| 成人高潮视频无遮挡免费网站| 亚洲精品日本国产第一区| 久久这里只有精品中国| 久久久久久久久中文| 成人无遮挡网站| 久久国产乱子免费精品| 亚洲av中文字字幕乱码综合| 卡戴珊不雅视频在线播放| 国产精品国产三级专区第一集| 色尼玛亚洲综合影院| 午夜福利高清视频| 超碰97精品在线观看| 麻豆国产97在线/欧美| 99视频精品全部免费 在线| 99久国产av精品| 91久久精品国产一区二区三区| h日本视频在线播放| 人妻系列 视频| 久久人人爽人人爽人人片va| 成年版毛片免费区| 天堂俺去俺来也www色官网 | 免费看av在线观看网站| 男女边吃奶边做爰视频| 国产有黄有色有爽视频| 午夜精品国产一区二区电影 | 观看美女的网站| 人妻夜夜爽99麻豆av|