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

    A New Safety Assessment Method Based on Belief Rule Base With Attribute Reliability

    2021-10-23 02:21:38ZhichaoFengWeiHeZhijieZhouXiaojunBanChanghuaHuandXiaoxiaHan
    IEEE/CAA Journal of Automatica Sinica 2021年11期

    Zhichao Feng,Wei He,Zhijie Zhou,Xiaojun Ban,Changhua Hu,and Xiaoxia Han

    Abstract—Safety assessment is one of important aspects in health management.In safety assessment for practical systems,three problems exist:lack of observation information,high system complexity and environment interference.Belief rule base with attribute reliability (BRB-r)is an expert system that provides a useful way for dealing with these three problems.In BRB-r,once the input information is unreliable,the reliability of belief rule is influenced,which further influences the accuracy of its output belief degree.On the other hand,when many system characteristics exist,the belief rule combination will explode in BRB-r,and the BRB-r based safety assessment model becomes too complicated to be applied.Thus,in this paper,to balance the complexity and accuracy of the safety assessment model,a new safety assessment model based on BRB-r with considering belief rule reliability is developed for the first time.In the developed model,a new calculation method of the belief rule reliability is proposed with considering both attribute reliability and global ignorance.Moreover,to reduce the influence of uncertainty of expert knowledge,an optimization model for the developed safety assessment model is constructed.A case study of safety assessment of liquefied natural gas (LNG)storage tank is conducted to illustrate the effectiveness of the new developed model.

    I.INTRODUCTION

    WITH the increasing complexity of practical system,the difficulty of health management increases.Safety assessment is one of important means to ensure the reliability of system and it is applied to prevent failure [1]–[4].

    The safety assessment model applies multi-source observation information to estimate the safety state of complex system [5],[6].In the safety assessment for complex system,many excellent works have been conducted [7]–[11].For example,Ahmadiet al.[8] developed a safety assessment model based on data-driven model and a theorem of safety assessment is proposed,Syd Aliet al.[9] constructed a safety assessment framework based on fault tree analysis for automatic teller machine (ATM)system,Zhouet al.[10]proposed a hidden behavior prediction model for complex system based on belief rule base (BRB).

    In safety assessment for practical systems,three problems exist.Firstly,due to the high cost of complex system,large amounts of observation data cannot be gathered,especially fault data [12],[13].Moreover,a complex system contains many subsystems and each subsystem is interrelated.Experts can provide some knowledge about the complex system [14].However,the expert knowledge is not enough to build accurate safety assessment model and the accuracy of the knowledge-based safety assessment is limited in practical system [12].Thirdly,influenced by disturbance in engineering practice,the gathered system information may not be fully reliable and noise exists in the information [15].When the unreliable system information is used in the safety assessment model,its modeling accuracy is decreased.The above three problems in practical systems pose challenges for safety assessment.Hence,the gathered unreliable observation data and the expert knowledge should both be applied.The BRB model is developed based on fuzzy theory,Dempster-Shafer theory and IF-THEN rule,which can handle uncertain,vague and incomplete information [14],[16]–[18].Then,the BRB model with attribute reliability (BRB-r)is proposed and it is a more general BRB model,where the attribute reliability is further considered [15].In the BRB-r model,the environment noise is addressed by attribute reliability.So,in this way,the BRB-r model provides an excellent framework to address the above three problems.Thus,in this paper,a safety assessment model is constructed based on BRB-r.

    In the BRB-r model,the belief rule is constructed based on attributes by experts [19],[20].When the unreliable attributes are used in the BRB-r model,the reliability of belief rule is decreased and its output reliability is influenced.Moreover,due to the complex environment and the limitation of the experts themselves,the expert knowledge applied in the belief rule may not be fully reliable.Consequently,influenced by these two factors,the belief rule in the constructed model is unreliable and the reliability of its output is also decreased.In BRB-r,the belief rules are combined by the evidential reasoning (ER)algorithm to generate safety state of practical system [21],[22].It should be noted that in existing works of BRB-r,only the influence of unreliable attribute has been considered in the reliability of belief rule [15].The influence of unreliable expert knowledge cannot be addressed.On the other hand,the BRB-r model is an expert system that complexity is determined by the amount of belief rule.When the number of belief rules of BRB-r is too small,its modeling ability for complex system is influenced.On the contrary,once the number of belief rules is too large,the BRB-r model complexity increases and the constructed safety assessment model cannot be trained by the small amount of observation data [23]–[25].Thus,for safety assessment of complex system,some belief rules in BRB-r should be reduced and its structure should be adaptively adjusted according to the practical system.In existing studies of BRB,some researches have been conducted to optimize the model structure based on the effectiveness of the belief rule [19],and these methods are essentially data-driven structural adjustment method.In these methods,the reliability of attribute and expert knowledge is assumed to be fully reliable.Once the input characteristics and expert knowledge are unreliable,the belief rule is unreliable and its effectiveness has no meaning,and the assessment accuracy is decreased.Thus,to improve the accuracy of safety assessment model and reduce the model complexity adaptively,the unreliable belief rule should be reduced according to the influence of complex working environment.To solve the above two problems in safety assessment for complex system,a new safety assessment model based on BRB-r is developed in this paper,where a calculation method for the belief rule reliability is proposed.The belief rule in BRB-r is constructed based on the IF-THEN rule [26]–[28].Its structure is determined by experts and the unreliable attributes are used as inputs.The reliability of attributes and expert knowledge both influence the reliability of the belief rule.Once the expert knowledge is unreliable and uncertain,the constructed belief rule is unreasonable that influences its discriminate ability for safety states.The output framework of belief rule in BRB-r consists of singleton propositions and universal set.The residual belief degree assigned to the universal set is the global ignorance that is an important way to address uncertainty.The greater the global ignorance is,the less the belief rule will recognize the output state and the less effect a rule has on the model output.The global ignorance is determined by the uncertainty of attributes and the way the belief rule is constructed [20],[29].So,in this way,the global ignorance represents the influence of uncertainty of attributes and the expert knowledge applied to construct the belief rule.It should be noted that the influence of expert knowledge comes from its uncertainty and unreliability.The uncertainty and unreliability denote two different aspects of the expert knowledge.The reliability of expert knowledge is its objective aspect determined by the ability of experts while the uncertainty of expert knowledge is its subjective aspect determined by the personal judgement of experts.Moreover,uncertainty exists in the natural language of the experts that is influenced by the interference factors such as the working environment and the mood of the experts.Therefore,the uncertainty and unreliability of expert knowledge should both be handled in the constructed safety assessment model.Based on the above analysis,in this paper,a calculation method of the belief rule reliability is proposed with considering both attribute reliability and global ignorance for the first time.Then,the belief rule reliability is applied to adjust the structure of the developed model.The reduction threshold of belief rule reliability is determined according to the requirement of the practical system.It should be noted that the reduction threshold of the belief rule reliability should consider both the model complexity and modeling accuracy.In this paper,the reduction threshold of belief rule reliability is determined by experts.To deal with the influence of uncertainty expert knowledge,a new optimization model is constructed based on the projection covariance matrix adaption evolution strategy (P-CMA-ES)algorithm [15],[30].In the optimization model,the belief rule is adaptively reduced based on its reliability while the model parameters are trained.It should be noted that the parameters in the reduced belief rule are no longer used as training parameters.That is to say,after the training part,the number of parameters in the safety assessment model will be reduced,and its modelling accuracy can be guaranteed.

    The reminder of this paper is organized as follows.In Section II,the problems in the safety assessment for complex system are formulated and then a new safety assessment model based on BRB-r is developed.The calculation method of belief rule reliability is proposed in Section III with considering both attribute reliability and global ignorance.Then,in Section IV,the inference process of the developed model is presented.A case study is conducted in Section V.This paper is concluded in Section VI.

    II.PROBLEM FORMULATION

    In safety assessment,the amount and reliability of system observation information will directly affect the estimation accuracy of the safety assessment model.In this section,the problems in safety assessment for practical system are formulated and then a new safety assessment model based on BRB-r is developed.

    A.Problem Formulation of Safety Assessment

    There are three problems in safety assessment that can be concluded as follows.

    Problem 1:In engineering practice,limited by the monitoring costs and monitoring methods,a large number of high-value observation data cannot be gathered.In addition,for practical system,although some expert knowledge can be provided,the experts cannot construct its accurate mathematical model.Thus,the first problem in safety assessment for practical system is how to aggregate the gathered multi-source information to construct the following safety assessment model [31],[32].

    whereS(t)is the estimated safety state of practical system.x1(t),x2(t),...,xM(t)are theMcharacteristics.Eis the vector of model parameters andRdenotes the expert knowledge applied in the constructed model.

    Problem 2:In engineering practice,complex environment interference factors influence the reliability of the observation information and further influence the estimation accuracy of the safety assessment model.Moreover,when the safety assessment model is too complicated,the update time of model structure and parameters is long,and its modeling timeliness is limited,and a lot of system information is needed to construct an accurate safety assessment model.Then,how to adjust the safety assessment model with considering both model accuracy and complexity is the second problem.

    Problem 3:For practical system,experts cannot provide accurate information to construct the safety assessment model,and uncertainty exists in the expert knowledge.So,in this way,the initial safety assessment model may not be able to accurately assess the safety state of the practical system.Moreover,the unreliable gathered information also influences the model accuracy.Therefore,how to adjust the model structure and train the model parameters simultaneously is the third problem that needs to be solved.

    B.Construction of the New Safety Assessment Model

    To solve the above three problems in engineering practice,a new BRB-r based safety assessment model is proposed in this subsection.

    In the new safety assessment model,there are multiple belief rules and thekth rule is profiled as [15]

    wherex1,x2,...,xMare theMcharacteristics of the practical system.denote theMreference points of thekth belief rule of the system characteristics and they are determined by experts.D1,D2,...,DNare the safety state levels and β1,k,β2,k,...,βN,kare their corresponding belief degree.θkis the weight of thekth belief rule.δ1,δ2,...,δMare the weight of theMcharacteristics andr1,r2,...,rMrepresent their radiabilities.Lis the number of the belief rule in the constructed safety assessment model based on BRB-r.Rkdenotes the reliability of thekth belief rule and it is calculated based on the characteristic reliability and global ignorance.The belief rule reliability can quantitatively measure the influence of unreliable characteristics and unreliable expert knowledge.The modeling process of the developed safety assessment model is shown in Fig.1.

    III.CALCULATION METHOD OF BELIEF RULE RELIABILITY

    The reliability of belief rule in BRB-r represents the influence degree of the unreliable attributes and unreliable expert knowledge to the belief rule.In the constructed safety assessment model,the belief rule reliability is applied as the key characteristic for model reduction.Thus,in this section,with considering both attribute reliability and expert knowledge reliability,a calculation method of belief rule reliability is proposed.

    Fig.1.The safety assessment model with both considering model complexity and accuracy.

    In the BRB-r based safety assessment model,the reliability of belief rule consists of two parts:the attribute reliability and the expert knowledge reliability.The reliability of belief rule is represented by

    The authors in [15] proposed a calculation method of attribute reliability based on the statistical method,where the disturbance factors are assumed to be static and the attribute reliability is a constant.In this paper,assume that the disturbance factors are static.Hence,the attribute reliability is calculated by

    whereTis the number of observation data.yirepresents the sum of unreliable observation data that can be eliminated by the fluctuation range.τ is the adjustment coefficient of the fluctuation range and represents the requirement of the complex system.When the observation data is reliable and little noise exists,the observation data can represent the system information andyi,j=1,j=1,2,...,T;otherwise,yi,j=0.Then,.

    In the output of belief degree,the global ignorance is the residual belief degree that cannot be distinguished which safety state should be assigned.It is generated because of the influence of expert knowledge and attribute.When the expert knowledge is unreliable,the ability of the constructed belief rule to discriminate the safety states of complex system is decreased.In other words,limited by the ability of experts,the recognition ability of the constructed belief rule is constrained.Thus,the global ignorance in the output of the belief rule can be used to calculate the reliability of expert knowledge [15],[33].Then,the reliability of belief rule can be calculated by

    wherer1,r2,...,rMare the comprehensive radiabilities of attributes.(1?βD,k/κ)is the reliability of expert knowledge in thekth belief rule,and κ is the scale factor that represents the influence ratio of attribute uncertainty and unreliability expert knowledge.As shown in the above formula,once the global ignorance βD,kis 0,the unreliability expert knowledge does not influences the belief rule reliability and the reliability of the belief rule is only influenced by the reliability of attributes.

    Remark 1:There is need to explain the effectiveness of attribute uncertainty,attribute reliability,expert knowledge uncertainty and expert knowledge reliability.In the BRB-r based safety assessment model,the uncertainty in attribute exists in qualitative knowledge,such as natural language,and it is mainly caused by incomplete information,inconsistent source,inaccurate expression.For the attribute reliability,it represents the influence of the disturbance factors,such as environment noise,and it is the objective aspect of the characteristic.The uncertainty of expert knowledge is the incompleteness,approximation and randomness in the natural language,and it may be influenced by the environment and other factors.In the BRB-r based safety assessment model,the uncertainty of expert knowledge can be addressed by the optimization model based on the gathered observation data.The reliability of expert knowledge represents the expert’s own ability and has nothing to do with external factors.

    IV.THE SAFETY ASSESSMENT MODEL BASED ON BRB-R

    In this section,a new safety assessment model based on BRB-r with considering belief rule reliability is developed.In the new safety assessment model,its structure is adaptively adjusted based on the belief rule reliability while its estimation accuracy is the optimization goal.So,in this way,the model complexity and accuracy can be both considered.In this section,the modeling process of the constructed safety assessment model based on BRB-r is presented in Section IVA,and then an optimization model is constructed based on the P-CMA-ES algorithm in Section IV-B.In Section IV-C,the modeling process of the developed model for complex system is presented.

    A.Inference Process of the New Safety Assessment Model

    The reliability of belief rule is determined by the attribute reliability and the expert knowledge reliability,and its calculation method is shown in Section III.Influenced by the unreliability belief rule,the accuracy of the estimated safety state generated by the ER algorithm is decreased.In addition,the BRB-r model is an expert system that complexity is determined by the number of belief rules.Thus,the unreliable belief rule should be eliminated.In general,in the BRB-r based safety assessment model,the unreliable belief rule will not only affect the accuracy of the model output,but also increase the complexity of the model.Thus,the structure and accuracy of the BRB-r based safety assessment model should both be considered and the inference process is presented in this section.

    When the observation data is available,due to the different formats of observation data for different characteristics,the observation data of characteristics cannot be aggregated directly.Thus,the observation data should be converted into a uniform format by the following formula [34],[35].

    wherexi(k+1)andxikare the reference values of theith characteristic in the (k+1)th and thekth belief rule,respectively.denotes the input of theith characteristic.Lis the numbers of the belief rules in the constructed model.

    Based on the frame of the belief rule as shown in (2),the total matching degree of the multiple characteristics can be obtained by the following formulas.

    where αkis the matching degree of theMinput characteristics to thekth belief rule.Cidenotes the new weight with considering the weight and reliability of theith attribute.represents the relative weight of theith attribute [15].Note that 0≤≤1.

    Once the observation information of system is available,it will activate the belief rules in the safety assessment model and the activation weights can be calculated by

    wherewkis the activation weight of thekth belief rule,which is determined by the input of system characteristics.The activation weight represents the effectiveness of the model inputs to the belief rule.θkdenotes the weight of thekth belief rule andLis the number of the belief rule in the BRB-r model.

    When the belief rule is unreliable,the estimated output is influenced and it increases the complexity of the safety assessment model.Thus,these belief rules should be eliminated according to the practical system.To adaptively adjust the structure of safety assessment model in the optimization model,the belief rule is reduced based on its reliability by

    whereRthresholdis the reduction threshold of the belief rule,and it can be determined by experts according to the model complexity and accuracy requirement.Once the reliability of belief rule cannot meet the requirement,the effectiveness of the belief rule is set to 0 andwk=0.So,in this way,the structure of the constructed safety assessment model can be adjusted.

    Then,after the activation weights of the belief rules are calculated,they can be combined by the ER algorithm and its analytic algorithm is shown as

    where βnis the belief degree of thenth safety state generated by the developed safety assessment model.Note that+βD=1.βDis the residual belief degree unassigned to any safety state,and it is influenced by the uncertainty of expert knowledge and unreliable expert knowledge.In the new safety assessment model,βDis used to calculate the belief rule reliability.L′is the number of the belief rule that satisfies the reliability requirement.

    After the belief degree of each safety state has been obtained,the final safety state of complex system can be calculated by

    whereu(Di)is the utility of theith safety levelDi,which is determined by experts according to the safety state of the practical system.denotes the final estimated safety state of the complex system.

    Remark 2:In the inference process of the developed safety assessment model,the attribute reliability and the residual belief degree are used to calculate the belief rule reliability.Once the belief rule reliability is lower than the thresholdRthreshold,the weight of the belief rule is set to 0 and it will not influence the final output of the safety assessment model.

    B.Optimization Model for the New Safety Assessment Model

    The initial safety assessment model is constructed based on the expert knowledge.Due to the complexity of practical system,the experts cannot provide accurate system information and uncertainty exists in the expert knowledge.Moreover,limited by the ability of experts,the expert knowledge may not be fully reliable.So,in this way,the safety assessment model constructed by the expert knowledge cannot estimate the safety state of complex system accurately.Thus,the gathered observation data should be applied to adaptively adjust the structure and parameters according to the practical system,and an optimization model for the new safety assessment model is developed in this subsection.

    In the BRB-r based safety assessment model,the optimization parameters consist of attribute weights,belief rule weights and belief degree in output of belief rules.There are two parts in the optimization model:the structure adjustment and parameter optimization.For the structure adjustment of the safety assessment model,once the reliability of belief rule cannot meet the requirement,its weight is set to 0 and its output belief degree is not trained in the optimization model.In addition,the BRB-r model is an expert system whose parameters must have physical meaning.Therefore,it parameters should be trained in constraints.Based on the original physical meaning,the constraints of optimization parameters can be listed as

    The mean square error (MSE)between the estimated safety state and actual safety state can represent the accuracy of the safety assessment model,and it can be calculated by the following formula [30].

    whereoutputactualis the actual safety state of the complex system,which is determined by experts according to the system mechanism.outputestimateddenotes the estimation output of the safety assessment model and it is calculated by(14).

    The optimization goal is to minimize the error between the actual safety state and the estimated safety state by the constructed model.Hence,the optimization model of the safety state assessment model can be profiled as

    Remark 3:In the optimization model,the parameters need to be trained and the model structure is adjusted adaptively by the reliability of belief rule simultaneously.It should be noted that when the belief rule reliability is lower than the threshold,the rule weight is set to 0 and its output belief degree is not trained in the training part.That is to say,the belief rule does not work in the constructed safety assessment model.So,in this way,the model structure can be adjusted according to the practical system,and the constructed safety assessment model can address the influence of the unreliable characteristics and the unreliable expert knowledge.Thus,by the optimization process,the accuracy of the optimization model can be improved while its model complexity is decreased.

    C.Modeling Process of the New Safety Assessment Model

    The modeling process of the new safety assessment model can be concluded by the following steps.

    Step 1:Construct the safety assessment model based on BRB-r.Note that the structure and parameters of the initial safety assessment model are both determined by experts.

    Step 2:Calculate the reliability of characteristics based on the calculation method presented in Section III,and the reliability reduction threshold of belief rule is determined by experts.

    Step 3:Reduce the belief rule based on the attribute reliability and global ignorance in belief rule output.At the same time,the parameters in the residual belief rules are optimized by the training data.

    Step 4:Test the optimized safety assessment model by the testing data.Note that in the testing part,the structure of the safety assessment model is not needed to be adjusted adaptively while its parameters of the safety assessment model are not changed.

    V.CASE STUDY

    To illustrate the effectiveness of the developed model,a case study of safety assessment for liquefied natural gas(LNG)storage tank in Hainan province,China,is conducted.

    A.Problem Formulation in Safety Assessment Model for LNG Storage Tank

    In this section,the LNG storage tank in Hainan province,China,is selected as the experimental object.The LNG is applied as one of the important energy sources that is increasingly used in China.In Hainan province,two storage tanks have been built by reclaiming land from the sea as shown in Fig.2.In this way,the foundation of the storage tank is unstable,and influenced by the excessive pressure of the storage tank,uneven settlement may occur in the LNG storage tank.Once the storage tank is unevenly settled,it will cause cracks in the storage tank,which will further lead to the LNG leakage and safety accidents.Thus,it is necessary to estimate the safety state accurately [10],[36].

    Fig.2.LNG storage tank in Hainan province,China.

    In the safety assessment of the LNG storage tank,due to its high cost,a large number of experiments cannot be conducted.Therefore,the number of observation data that can be gathered is small,and an accurate safety assessment model based on the quantitative data cannot be established.Moreover,due to the influence of the disturbance factors in practical environment,noise exists in the gathered observation information,and its reliability is decreased.On the other hand,the LNG storage tank is a typical complex system,where many subsystems exist and they are related to each other.The experts cannot provide accurate information for the storage tank and an analytical model cannot be obtained.So,it is necessary to aggregate the expert knowledge and the gathered observation data to construct the safety assessment model for the LNG storage tank,and the unreliability of the gathered information needs to be solved.When the expert knowledge and the characteristics are unreliable,the reliability of the belief rule is influenced,which further influences the accuracy of the safety assessment model.In addition,the number of belief rule determines the complexity of the safety assessment model.Consequently,to improve the estimation accuracy and decrease the model complexity,in the training part,the belief rule in the rule base should be adjusted by its reliability adaptively according to the actual working environment of the LNG storage tank.Thus,in the experiment,a new safety assessment model for the LNG storage tank is constructed based on the developed model in this paper.

    In the experiment,12 monitoring points are installed in the storage tank,and 42 sets of observation data are gathered from 2012 to 2019.Based on the analysis of the structure and location of the storage tank,the first monitoring point (FP)and the seventh monitoring point (SP)are selected by experts as two attributes of the BRB-r model.The observation data of these two attributes are shown in Fig.3.It can be seen from Fig.3 that the LNG storage tank’s settlement is irregular over time and has a large change range.Based on the building deformation measurement specifications (JGJ8-2007),the greater the gap between the two monitoring points is,the lower the safety state of the LNG storage tank is.Hence,the actual safety state of the LNG storage tank is determined by experts at each monitoring point.It should be noted that the actual safety state of the LNG storage tank should also consider its location.

    Fig.3.The observation data of FP and SP.

    Remark 4:In engineering practice,there are several factors that lead to the LNG storage tank failure,such as the uneven settlement,uneven filling of the insulation material,the failure of the vacuum and the residual gas inside the tank,etc.Based on the investigation in Hainan LNG Company,the uneven settlement should be paid more attention to in the face of other failure factors by the construction mode of the storage tank.So,in this paper,the uneven settlement is selected as the major concern and assume that the other aspects of the storage tank are normal.

    B.Construction of the Adaptive Safety Assessment Model for LNG Storage Tank

    To solve the above analyzed problems and improve the safety assessment accuracy for the LNG storage tank,in this subsection,an adaptive safety assessment model is constructed based on the developed model.

    In the constructed safety assessment model,the first monitoring point and the seventh monitoring point are selected as two attributes,represented byFPandSP,respectively.The safety state of the LNG storage tank is divided into three levels:the high safety state (H),the medium safety state (M)and the low safety state (L).Based on (2),the belief rule in the BRB-r based safety assessment model can be profiled as

    whereyis the estimated safety state of the LNG storage tank.In the training process,once the reliability of thekth belief ruleRkdoes not satisfy the requirement,its rule weight θkis set to 0 and it has no effect in the safety assessment model.

    The reference point indicates the degree of division of the characteristics.It should be noted that the number of the reference points determines the complexity and accuracy of the safety assessment model.Thus,in this experiment,five reference points ofFPandSPare selected:low (L),slightly low (SL),medium (M),slightly high (SH),high (H).For the reference values ofFPandSP,they should be determined by experts in accordance with the three principles of uniform distribution,covering the variation range of two characteristics and focusing on the data-intensive area.Thus,combined with the characteristics of the storage tank under different working states,the reference values ofFPandSPare given by experts as shown in Tables I and II.Moreover,three reference points are selected for safety state and their reference values are shown in Table III.

    Based on the reference value of both attributes and (21),the parameters in the initial safety assessment model can be given by experts.As shown in Table IV,all the belief rules are assumed to be fully important and two characteristics are assumed to be fully important.In the initial safety assessment model,the expert knowledge is assumed to be fully reliable and the global ignorance of each belief rule is set to 0.

    TABLE I REFERENCE POINT AND VALUE OF FP

    TABLE II REFERENCE POINT AND VALUE OF SP

    TABLE III REFERENCE POINT AND VALUE OF THREE SAFETY STATES

    C.Training and Testing of the Constructed Safety Assessment Model

    In the initial safety assessment model,due to the uncertainty of expert knowledge,its estimation accuracy is influenced.Moreover,as shown in Table IV,there are 25 belief rules in the safety assessment model and 127 parameters exist.However,in the experiment,only 42 sets of observation data have been gathered,which cannot adequately train the model parameters.On the other hand,due to the disturbance factors in actual working environment and the limitation of the expert’s own capabilities,the characteristics and the expert knowledge are not fully reliable.As profiled in (21),once the characteristics and the expert knowledge are unreliable,the reliability of the belief rule in the assessment model is influenced.These unreliable belief rules not only influence the accuracy of the safety assessment model but also increase the model complexity.Thus,these unreliable belief rules should be eliminated.Based on the optimization model developed in Section IV-B,in the training process,the structure of the safety assessment model is adjusted while the parameters are optimized.

    In this section,the reduction threshold of the reliability of belief rule is determined by experts andRk=0.3,it is determined by experts according to the estimation accuracy and model complexity.On the basis of the training data and the constructed optimization model,the safety assessment can be optimized as shown in Table V.The estimated output of the new safety assessment model is shown in Fig.4.As shown in Fig.4,the developed model can estimate the safety state of the LNG storage tank accurately,and the MSE of the developed safety assessment model is 0.0509.There are 22 belief rules and 112 parameters in the optimized safety assessment model.The number of the parameters of the optimized model reduces 11.8% and the number of belief rule reduces 12% compared with the initial safety assessment model.Although some observation points cannot be accurately assessed,the error does not exceed half of the actual safety state of the LNG storage tank,which can still provide a reference for decision-making in engineering practice.

    Fig.4.Estimated safety state of the developed model.

    TABLE IV INITIAL BRB-R BASED SAFETY ASSESSMENT MODEL

    TABLE V OpTIMIzED BRB-R BASED SAFETY ASSESSMENT MODEL

    Fig.5.Relationship between reduction threshold of belief rule reliability and number of belief rule.

    Fig.6.Relationship between MSE of the developed model and reduction threshold of belief rule reliability.

    Fig.7.Relationship between error of the developed model and reduction threshold of belief rule reliability.

    The complexity of the developed safety assessment model is determined by the number of the belief rules,and there is a direct relationship between the reduction threshold of belief rule reliability and the complexity of the safety assessment model.As shown in Fig.5,when the reduction threshold of the belief rule reliability changes from 0.1 to 1,the amount of belief rule gradually decreases.The experiment is conducted for 40 times.WhenRk>0.3,the number of belief rules quickly decreases,and once the reduction threshold is greater than 0.6,the number of belief rules becomes 0.On the other hand,as the reduction threshold of the belief rule reliability changes,the accuracy of the safety assessment model also changes.As shown in Figs.6 and 7,when the reduction threshold of belief rule reliability changes from 0.1 to 1,the MSE and error of the constructed safety assessment model first decrease and then increase,and whenRk>0.6,the model MSE and error between actual safety state and estimated output become infinite.That is to say,as the reduction threshold increases from 0.1 to 1,the estimation accuracy of the developed model first increases and then decreases,and the best accuracy of the safety assessment model isRk=0.3.Moreover,as the reduction thresholdRkincreases,the distribution range of MSE first decreases and then increases.Thus,the robustness of the safety assessment model first increases and then decreases.Consequently,in this way,we can obtain a conclusion that the reduction threshold of belief rule reliability can influence the complexity and accuracy of the developed safety assessment model.For practical systems,the reduction threshold of the belief rule reliability cannot be greater than 0.5,and 0

    D.Comparative Studies

    To illustrate the effectiveness of the developed safety assessment model in this paper,comparative studies are conducted in this subsection.

    First,the safety assessment model developed in this paper is constructed based on BRB-r while its structure is adjusted according to the belief rule reliability.In the safety assessment model based on the traditional BRB-r model,its structure is not changed and only the model parameters are trained.Compared with the traditional BRB-r model,the complexity of the proposed safety assessment model in this paper decreases.Moreover,the BRB-r model is one of the expert systems,where the observation data and expert knowledge are aggregated simultaneously.For the LNG storage tank,there are only 42 monitoring data,which is a typical small sample problem,and it is also a complex system.Consequently,in this subsection,comparative studies with the traditional BRBr model,BP neural network and fuzzy inference method are conducted.

    The MSEs of the developed safety assessment model in this paper,the safety assessment models developed by the traditional BRB-r model,the BP neural network and the fuzzy inference method are shown in Table VI.The MSE of the fuzzy inference method is 0.1160,where the initial rule base is determined by experts and the inputs are the observation data of the two characteristics of the LNG storage tank [12].Moreover,the MSE of the BP neural network is 0.0837,where the training data and the testing data are same as the developed model in this paper,and the number of iterations and learning rate are set to 500 and 0.1,respectively.Therefore,compared with the BP neural network and the fuzzy inference method,the new developed safety model based on BRB-r can aggregate the observation data and expert knowledge,and its estimation accuracy is improved by 56.1%and 39.1%,respectively.The MSE of the safety assessment model based on traditional BRB-r is 0.0722,where the model structure is not changed.Assume that the expert knowledge applied in the two BRB-r models is same.That is to say,the reliability of expert knowledge in the two models is same that can be handled by the new developed model.Compared with the traditional BRB-r model,the accuracy of the proposed model improves 29.5%.Moreover,compared with the safety assessment model developed in this paper,there are 25 belief rules and 127 parameters in the safety assessment model developed by the traditional BRB-r model.The complexity of the new proposed model is decreased and its accuracy is higher than that of the traditional BRB-r model.Thus,the developed model can handle the unreliable observation data and the unreliable expert knowledge better,which is more suitable for modeling practical systems.

    TABLE VI COMpARATIvE STUDIES

    VI.CONCLUSION

    In this paper,a new safety assessment model with model structure adaptive adjustment is developed to solve three problems in engineering practice:lack of observation data,high system complexity and environment interference factors.To deal with the influence of the unreliable belief rule to the model accuracy and complexity,a new calculation method of belief rule reliability is proposed,where the reliability of system characteristics and the reliability of expert knowledge are both considered.In the optimization model of the developed safety assessment model,its structure is adjusted adaptively while the model parameters are optimized.In this way,the new belief rule base with attribute reliability (BRB-r)based safety assessment model can be adjusted according to the practical system.

    There are three innovations in this paper.Firstly,for practical system,to solve the lack of observation data,high complexity of system and complex environment interference factors,a new safety assessment model based on BRB-r is developed.Secondly,limited by the ability of the experts and influenced by the unreliable observation data,the belief rules in the constructed safety assessment model may not be fully reliable.These unreliable belief rules cannot estimate the safety state of complex system,but increase the complexity of the safety assessment model.Therefore,a new calculation method of belief rule is developed to reduce the model structure.Lastly,to consider both the model accuracy and complexity,an optimization model is constructed to adaptively adjust the model structure based on belief rule reliability and optimize the model parameters.

    In the case study for the LNG storage tank in Hainan province,China,only the uneven settlement is considered and we assume that the storage tank is not affected by other factors,such as uneven filling of the insulation material.In practical system,these failure causes should be considered simultaneously.Moreover,in this paper,the developed BRB-r based safety assessment model is a general method.Some specific problems in other areas,such as traffic control,need to be further studied and the BRB-r based safety assessment model should also be adjusted according to the subject.These challenges need to be solved in future works.

    久久久久国产网址| 欧美高清性xxxxhd video| 免费播放大片免费观看视频在线观看 | 热99re8久久精品国产| 国产麻豆成人av免费视频| 别揉我奶头 嗯啊视频| 中文字幕av在线有码专区| 成人亚洲欧美一区二区av| 亚洲精品,欧美精品| 99在线人妻在线中文字幕| 亚洲av免费在线观看| 欧美精品一区二区大全| 一级黄色大片毛片| 岛国毛片在线播放| 成人亚洲欧美一区二区av| 男的添女的下面高潮视频| 在线观看av片永久免费下载| 精品一区二区免费观看| 丰满少妇做爰视频| 精品久久久久久成人av| 欧美激情久久久久久爽电影| 国内揄拍国产精品人妻在线| 好男人视频免费观看在线| 国产精品蜜桃在线观看| 欧美激情久久久久久爽电影| 亚洲美女视频黄频| 最后的刺客免费高清国语| 亚洲欧美精品综合久久99| 久久亚洲国产成人精品v| 国产乱人视频| 国产黄a三级三级三级人| 日本欧美国产在线视频| 欧美潮喷喷水| 国产av在哪里看| 丰满乱子伦码专区| 桃色一区二区三区在线观看| 亚洲性久久影院| 欧美又色又爽又黄视频| 成人一区二区视频在线观看| 天堂中文最新版在线下载 | 国产精品.久久久| eeuss影院久久| 日本免费a在线| 国产乱来视频区| 国产午夜精品一二区理论片| 久久99蜜桃精品久久| 久久久精品欧美日韩精品| 午夜激情福利司机影院| 熟妇人妻久久中文字幕3abv| 成人国产麻豆网| 色播亚洲综合网| 99久国产av精品| 不卡视频在线观看欧美| 亚洲怡红院男人天堂| 最近中文字幕高清免费大全6| 超碰97精品在线观看| av在线亚洲专区| 国产极品天堂在线| 狠狠狠狠99中文字幕| 国产一级毛片在线| 三级男女做爰猛烈吃奶摸视频| 亚洲怡红院男人天堂| 亚洲五月天丁香| 亚洲图色成人| 中文亚洲av片在线观看爽| 国产精品一及| 在线免费十八禁| 亚洲丝袜综合中文字幕| 2021天堂中文幕一二区在线观| 欧美性感艳星| 欧美成人午夜免费资源| 欧美三级亚洲精品| 欧美性猛交黑人性爽| 最近的中文字幕免费完整| 国产又色又爽无遮挡免| 国产三级在线视频| 国产高清国产精品国产三级 | 国产成人一区二区在线| 日韩一区二区视频免费看| 日本熟妇午夜| 成人毛片a级毛片在线播放| 日韩中字成人| 床上黄色一级片| 午夜福利在线观看吧| 久久草成人影院| 国内精品一区二区在线观看| 国产一区亚洲一区在线观看| 亚洲国产色片| 一个人看视频在线观看www免费| 美女大奶头视频| av视频在线观看入口| 亚洲久久久久久中文字幕| 夫妻性生交免费视频一级片| 久久精品综合一区二区三区| 国产亚洲一区二区精品| 亚洲成色77777| 国产精品久久久久久精品电影| 一区二区三区乱码不卡18| 国产色爽女视频免费观看| 麻豆久久精品国产亚洲av| 变态另类丝袜制服| 亚洲成人久久爱视频| 国产精品人妻久久久久久| 2021天堂中文幕一二区在线观| 搡老妇女老女人老熟妇| 国产精品麻豆人妻色哟哟久久 | 日韩欧美三级三区| 1024手机看黄色片| 内射极品少妇av片p| 男人舔女人下体高潮全视频| 男插女下体视频免费在线播放| 精品久久久噜噜| 天天躁夜夜躁狠狠久久av| 亚洲av电影在线观看一区二区三区 | 国产亚洲最大av| 麻豆乱淫一区二区| 成年版毛片免费区| 国产欧美另类精品又又久久亚洲欧美| 午夜激情福利司机影院| 国内精品一区二区在线观看| 天堂√8在线中文| 99热6这里只有精品| 伊人久久精品亚洲午夜| 黑人高潮一二区| 成人av在线播放网站| 欧美日本视频| av免费在线看不卡| 欧美丝袜亚洲另类| 欧美日韩国产亚洲二区| 午夜福利在线观看免费完整高清在| 全区人妻精品视频| 夜夜爽夜夜爽视频| 亚洲性久久影院| 你懂的网址亚洲精品在线观看 | 亚洲精品影视一区二区三区av| 麻豆av噜噜一区二区三区| 国产欧美日韩精品一区二区| 国产精品人妻久久久久久| 婷婷色麻豆天堂久久 | 午夜久久久久精精品| 国产成人一区二区在线| 精品国产露脸久久av麻豆 | 身体一侧抽搐| av国产免费在线观看| 嫩草影院精品99| 熟女人妻精品中文字幕| 在线播放国产精品三级| 亚洲三级黄色毛片| 一级二级三级毛片免费看| 国产老妇伦熟女老妇高清| 蜜桃亚洲精品一区二区三区| 国产伦精品一区二区三区四那| 三级国产精品欧美在线观看| 色吧在线观看| 亚洲欧美清纯卡通| 亚洲精品日韩在线中文字幕| 亚洲美女搞黄在线观看| 国产大屁股一区二区在线视频| 可以在线观看毛片的网站| 你懂的网址亚洲精品在线观看 | 日韩亚洲欧美综合| 级片在线观看| 免费观看在线日韩| 一个人免费在线观看电影| 亚洲精品一区蜜桃| 精品国产露脸久久av麻豆 | 看免费成人av毛片| 亚洲欧美精品专区久久| 超碰av人人做人人爽久久| 久久人人爽人人片av| 亚洲三级黄色毛片| 亚洲色图av天堂| 欧美日韩在线观看h| 精品久久久久久久久av| 日本一二三区视频观看| 精品久久久久久久久久久久久| 身体一侧抽搐| 嫩草影院入口| 亚洲国产日韩欧美精品在线观看| 99久久中文字幕三级久久日本| 18禁裸乳无遮挡免费网站照片| 亚洲丝袜综合中文字幕| 99热网站在线观看| 99久国产av精品| 国产成人精品久久久久久| 亚洲最大成人av| 婷婷色av中文字幕| 中文在线观看免费www的网站| 插阴视频在线观看视频| 一级黄片播放器| 中文字幕制服av| 亚洲在久久综合| 国产精品福利在线免费观看| 国产精品爽爽va在线观看网站| 国产成人aa在线观看| 大又大粗又爽又黄少妇毛片口| 日本五十路高清| 一区二区三区乱码不卡18| 国产在线男女| 黑人高潮一二区| 女人被狂操c到高潮| 插逼视频在线观看| 中国国产av一级| 99久久人妻综合| 22中文网久久字幕| 日韩av在线大香蕉| 在线播放无遮挡| 色综合亚洲欧美另类图片| 精品人妻熟女av久视频| 精品国产一区二区三区久久久樱花 | 激情 狠狠 欧美| 亚洲国产成人一精品久久久| 大话2 男鬼变身卡| av在线播放精品| 国产免费又黄又爽又色| 国产不卡一卡二| 亚洲国产色片| 国产精品久久久久久久电影| 免费搜索国产男女视频| 国产伦一二天堂av在线观看| 国产黄a三级三级三级人| 亚洲av不卡在线观看| 日韩欧美 国产精品| 两个人视频免费观看高清| 精品国内亚洲2022精品成人| 国产成人91sexporn| 国产乱来视频区| 亚洲电影在线观看av| 国产精品女同一区二区软件| 国产 一区 欧美 日韩| 国产精品国产三级专区第一集| 18禁裸乳无遮挡免费网站照片| 尾随美女入室| 99在线人妻在线中文字幕| 老司机福利观看| 午夜福利网站1000一区二区三区| 亚洲av熟女| 观看免费一级毛片| 天天一区二区日本电影三级| 男女那种视频在线观看| 国产激情偷乱视频一区二区| 成人亚洲欧美一区二区av| 亚洲精品日韩av片在线观看| 麻豆国产97在线/欧美| 18禁裸乳无遮挡免费网站照片| 美女高潮的动态| 成人午夜精彩视频在线观看| 久久久久久久久久久免费av| 国产不卡一卡二| 国产亚洲最大av| 欧美xxxx性猛交bbbb| 女的被弄到高潮叫床怎么办| 国产成人一区二区在线| 成人三级黄色视频| 亚洲aⅴ乱码一区二区在线播放| 老司机福利观看| 亚洲欧美成人精品一区二区| 波多野结衣高清无吗| 热99在线观看视频| 日韩av在线大香蕉| 观看免费一级毛片| 中文字幕av成人在线电影| 亚洲av电影在线观看一区二区三区 | 尤物成人国产欧美一区二区三区| 色5月婷婷丁香| 神马国产精品三级电影在线观看| 噜噜噜噜噜久久久久久91| 久久精品久久久久久久性| 赤兔流量卡办理| 久久99精品国语久久久| 午夜亚洲福利在线播放| 国产精品一区www在线观看| 欧美性感艳星| 韩国av在线不卡| 国产亚洲精品av在线| 只有这里有精品99| 观看美女的网站| 精品久久久久久成人av| 高清午夜精品一区二区三区| 久久精品综合一区二区三区| 欧美潮喷喷水| 长腿黑丝高跟| 中文字幕久久专区| 日韩欧美三级三区| 最近最新中文字幕大全电影3| 26uuu在线亚洲综合色| 99热这里只有是精品在线观看| 纵有疾风起免费观看全集完整版 | 91久久精品国产一区二区成人| 国产又色又爽无遮挡免| 成年女人永久免费观看视频| 国产探花在线观看一区二区| 国产成人精品婷婷| 成人综合一区亚洲| 国产亚洲91精品色在线| 亚洲人成网站高清观看| 爱豆传媒免费全集在线观看| 晚上一个人看的免费电影| 日本黄色视频三级网站网址| 亚洲18禁久久av| 在线免费观看不下载黄p国产| 内射极品少妇av片p| 亚洲伊人久久精品综合 | 免费一级毛片在线播放高清视频| 日韩av不卡免费在线播放| 偷拍熟女少妇极品色| 成年免费大片在线观看| 看十八女毛片水多多多| 国产精品久久久久久精品电影| 亚洲五月天丁香| 99热精品在线国产| 国产一区二区三区av在线| 毛片一级片免费看久久久久| 十八禁国产超污无遮挡网站| 神马国产精品三级电影在线观看| 51国产日韩欧美| 寂寞人妻少妇视频99o| 一级毛片我不卡| 人妻制服诱惑在线中文字幕| eeuss影院久久| 中国美白少妇内射xxxbb| 一级黄色大片毛片| 成人特级av手机在线观看| 婷婷色av中文字幕| 97超碰精品成人国产| 久久久亚洲精品成人影院| 国产精品爽爽va在线观看网站| 免费黄色在线免费观看| 小蜜桃在线观看免费完整版高清| 亚洲精品,欧美精品| 久久99热这里只频精品6学生 | 成人三级黄色视频| 亚洲精品成人久久久久久| 亚洲欧美精品专区久久| 亚洲内射少妇av| 成人鲁丝片一二三区免费| 国产成人freesex在线| 国产精品久久久久久精品电影| 亚洲成人久久爱视频| 国产伦一二天堂av在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产视频内射| 成年版毛片免费区| 国产成人福利小说| 精品99又大又爽又粗少妇毛片| 国产女主播在线喷水免费视频网站 | 国产精品1区2区在线观看.| 精品午夜福利在线看| 国产v大片淫在线免费观看| 两个人的视频大全免费| 久久久久九九精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 天堂√8在线中文| 亚洲最大成人av| 免费在线观看成人毛片| 网址你懂的国产日韩在线| 99在线人妻在线中文字幕| 亚洲人成网站在线播| 国产极品精品免费视频能看的| 激情 狠狠 欧美| videos熟女内射| 日本与韩国留学比较| 一区二区三区乱码不卡18| 91狼人影院| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久精品熟女亚洲av麻豆精品 | 国产免费福利视频在线观看| 国产黄色小视频在线观看| 十八禁国产超污无遮挡网站| 在线免费观看的www视频| 免费观看在线日韩| 久久久久久久久久久丰满| ponron亚洲| 性插视频无遮挡在线免费观看| 18+在线观看网站| 我的女老师完整版在线观看| 免费看av在线观看网站| 午夜福利在线观看免费完整高清在| 久久综合国产亚洲精品| 国产综合懂色| 欧美日韩在线观看h| 午夜福利视频1000在线观看| 国产探花极品一区二区| 亚洲av福利一区| 亚洲人成网站在线播| 亚洲精品亚洲一区二区| 日本一二三区视频观看| 搞女人的毛片| 亚洲aⅴ乱码一区二区在线播放| 看十八女毛片水多多多| 一区二区三区免费毛片| 亚洲中文字幕日韩| 乱系列少妇在线播放| 亚洲精品乱码久久久久久按摩| 欧美xxxx性猛交bbbb| 午夜福利在线在线| 久热久热在线精品观看| 99久久精品热视频| 亚洲人成网站高清观看| 97在线视频观看| 97超碰精品成人国产| 青青草视频在线视频观看| 国产精品99久久久久久久久| 男女视频在线观看网站免费| 插逼视频在线观看| 能在线免费观看的黄片| 亚洲乱码一区二区免费版| 女人被狂操c到高潮| 亚洲成人中文字幕在线播放| 高清毛片免费看| 七月丁香在线播放| 美女被艹到高潮喷水动态| a级毛片免费高清观看在线播放| 国产成人午夜福利电影在线观看| 69人妻影院| 女人久久www免费人成看片 | 十八禁国产超污无遮挡网站| 欧美日韩精品成人综合77777| 欧美色视频一区免费| 久久久久国产网址| 天堂√8在线中文| 尤物成人国产欧美一区二区三区| 亚洲精品日韩av片在线观看| 午夜激情欧美在线| 日本免费在线观看一区| 精品免费久久久久久久清纯| 中文字幕av在线有码专区| 亚洲av福利一区| 国产在视频线在精品| 欧美日韩在线观看h| 午夜激情福利司机影院| 你懂的网址亚洲精品在线观看 | 久久久久免费精品人妻一区二区| 国产亚洲91精品色在线| 天堂网av新在线| 变态另类丝袜制服| 亚洲,欧美,日韩| 丝袜美腿在线中文| 看十八女毛片水多多多| 一区二区三区四区激情视频| 国产伦理片在线播放av一区| 日韩强制内射视频| 麻豆一二三区av精品| 不卡视频在线观看欧美| 国产成人免费观看mmmm| 久久精品久久久久久久性| av在线观看视频网站免费| 在线免费观看不下载黄p国产| 熟女电影av网| 嫩草影院入口| 久久鲁丝午夜福利片| 一级毛片电影观看 | 国产精品综合久久久久久久免费| 性插视频无遮挡在线免费观看| 久久久久国产网址| 国产精品爽爽va在线观看网站| 国产大屁股一区二区在线视频| 日本wwww免费看| 内地一区二区视频在线| 亚洲国产精品sss在线观看| 一本一本综合久久| 久久亚洲国产成人精品v| 久久99蜜桃精品久久| 亚洲一区高清亚洲精品| 国产三级中文精品| 久久久久网色| 三级经典国产精品| 淫秽高清视频在线观看| 久久久欧美国产精品| 九九久久精品国产亚洲av麻豆| 日韩一区二区三区影片| 日韩人妻高清精品专区| 国产亚洲av嫩草精品影院| 国产精品1区2区在线观看.| 国产免费视频播放在线视频 | 亚洲无线观看免费| 欧美一区二区国产精品久久精品| 亚洲在线自拍视频| 久久精品影院6| 国产成人a区在线观看| 少妇丰满av| 一级av片app| 国产真实伦视频高清在线观看| 欧美xxxx性猛交bbbb| 亚洲av中文av极速乱| 日本一二三区视频观看| 人人妻人人澡人人爽人人夜夜 | 我要看日韩黄色一级片| 尾随美女入室| 久久久久久久午夜电影| 日韩欧美在线乱码| 国产精品日韩av在线免费观看| 免费一级毛片在线播放高清视频| 日韩三级伦理在线观看| 天堂影院成人在线观看| 人人妻人人澡欧美一区二区| 非洲黑人性xxxx精品又粗又长| 婷婷色麻豆天堂久久 | 麻豆久久精品国产亚洲av| 欧美色视频一区免费| 国产女主播在线喷水免费视频网站 | 激情 狠狠 欧美| 看黄色毛片网站| 麻豆国产97在线/欧美| 国产激情偷乱视频一区二区| 日韩强制内射视频| 91狼人影院| 免费看美女性在线毛片视频| 国产成人午夜福利电影在线观看| 久久精品综合一区二区三区| 亚洲av中文字字幕乱码综合| 久久久久久久午夜电影| 国产激情偷乱视频一区二区| 国内精品一区二区在线观看| 我要搜黄色片| 国产成人一区二区在线| 少妇裸体淫交视频免费看高清| 欧美日本视频| 免费看美女性在线毛片视频| 欧美+日韩+精品| 18禁在线播放成人免费| 两性午夜刺激爽爽歪歪视频在线观看| 综合色丁香网| 黄色配什么色好看| 久久午夜福利片| 99视频精品全部免费 在线| 欧美xxxx黑人xx丫x性爽| 色综合色国产| 久久午夜福利片| 看十八女毛片水多多多| 97热精品久久久久久| 看十八女毛片水多多多| 日韩高清综合在线| 国产黄片美女视频| 亚洲婷婷狠狠爱综合网| 精品不卡国产一区二区三区| 18禁在线无遮挡免费观看视频| 黄色配什么色好看| 国产黄片视频在线免费观看| 亚洲,欧美,日韩| 久久6这里有精品| 国产高清不卡午夜福利| 18禁在线无遮挡免费观看视频| 三级男女做爰猛烈吃奶摸视频| 丰满少妇做爰视频| 日本免费在线观看一区| 日本欧美国产在线视频| 日韩av不卡免费在线播放| 久久精品久久久久久久性| 日韩精品青青久久久久久| 亚洲欧美中文字幕日韩二区| 国产免费一级a男人的天堂| 精品一区二区三区视频在线| 日产精品乱码卡一卡2卡三| 桃色一区二区三区在线观看| 亚洲精品成人久久久久久| 亚洲久久久久久中文字幕| 七月丁香在线播放| 欧美97在线视频| 天堂√8在线中文| 日韩欧美在线乱码| 久久精品久久精品一区二区三区| 欧美不卡视频在线免费观看| 久久6这里有精品| 黄色日韩在线| 国产精华一区二区三区| 在线观看av片永久免费下载| 国产乱来视频区| 精品久久久久久久久亚洲| 别揉我奶头 嗯啊视频| 高清日韩中文字幕在线| 亚洲人成网站高清观看| a级毛片免费高清观看在线播放| 亚洲精华国产精华液的使用体验| 国产黄色小视频在线观看| av福利片在线观看| 国产精品熟女久久久久浪| 一区二区三区高清视频在线| 狂野欧美白嫩少妇大欣赏| 国产精品无大码| 亚洲天堂国产精品一区在线| 精品久久久久久久久av| 亚洲无线观看免费| 男女视频在线观看网站免费| 亚洲电影在线观看av| 国产精品一区www在线观看| 又粗又爽又猛毛片免费看| 男人狂女人下面高潮的视频| 草草在线视频免费看| 在线播放国产精品三级| 亚州av有码| 18+在线观看网站| 日日摸夜夜添夜夜爱| 色5月婷婷丁香| 日日干狠狠操夜夜爽| 高清午夜精品一区二区三区| 日本欧美国产在线视频| 久久精品人妻少妇| 色视频www国产| 精品一区二区免费观看| 黄色欧美视频在线观看| 久久久精品欧美日韩精品| 亚洲国产精品合色在线| 最近手机中文字幕大全| 亚洲av中文av极速乱| 亚洲欧美精品综合久久99| 国产精品av视频在线免费观看| 久久久精品大字幕| 长腿黑丝高跟| 尾随美女入室| 久久6这里有精品| 全区人妻精品视频| 国产黄色视频一区二区在线观看 | 老司机影院成人| 人人妻人人澡欧美一区二区|