亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于貝葉斯網(wǎng)絡(luò)的電力施工安全承載能力評(píng)估

        2017-06-28 15:09:33來駿陳蕾
        山東科學(xué) 2017年3期
        關(guān)鍵詞:貝葉斯概率要素

        來駿, 陳蕾

        (1.國家電網(wǎng)浙江省電力公司湖州供電公司,浙江 湖州313000;2.國家電網(wǎng)浙江省電力公司,浙江 杭州310007)

        ?

        基于貝葉斯網(wǎng)絡(luò)的電力施工安全承載能力評(píng)估

        來駿1, 陳蕾2

        (1.國家電網(wǎng)浙江省電力公司湖州供電公司,浙江 湖州313000;2.國家電網(wǎng)浙江省電力公司,浙江 杭州310007)

        提出了基于貝葉斯網(wǎng)絡(luò)的動(dòng)態(tài)承載能力評(píng)估模型。通過對承載能力風(fēng)險(xiǎn)要素的識(shí)別和評(píng)估,構(gòu)建了貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)和參數(shù),在此基礎(chǔ)上對電力施工安全承載能力進(jìn)行了風(fēng)險(xiǎn)評(píng)估和診斷。最后,依據(jù)分析得到的結(jié)果進(jìn)行風(fēng)險(xiǎn)控制。該系統(tǒng)可以滿足電力施工安全承載能力估計(jì)的需求。

        貝葉斯網(wǎng)絡(luò);安全承載能力;電力施工;風(fēng)險(xiǎn)評(píng)估

        隨著電網(wǎng)的發(fā)展,供電企業(yè)工作量猛增,生產(chǎn)班組工作量不斷加大,長時(shí)間超負(fù)荷工作,導(dǎo)致不安全作業(yè)時(shí)有發(fā)生,如觸電、火災(zāi)、設(shè)備漏電、設(shè)備短路、高空墜落、人身傷亡等。因此,需要開展生產(chǎn)單位安全承載能力分析。浙江省電力公司近年來在員工安全積分機(jī)制、外協(xié)隊(duì)伍評(píng)價(jià)和安全承載能力分析方面進(jìn)行了很多嘗試,但仍然存在一些問題:

        (1)安全承載能力分析模型不完善,主要以定性分析為主,無法與實(shí)際的安全生產(chǎn)過程緊密結(jié)合,難以真正發(fā)揮安全承載能力分析降低安全風(fēng)險(xiǎn)的作用。

        (2)員工安全積分規(guī)則簡單,主要適用于檢修、操作專業(yè)的一線員工,未涵蓋其他專業(yè)員工和管理人員,無法實(shí)現(xiàn)對員工安全積分的精益化管理。

        (3)外協(xié)隊(duì)伍的準(zhǔn)入門檻較低,評(píng)價(jià)考核手段單一,施工能力評(píng)價(jià)沒有量化標(biāo)準(zhǔn)。

        因此,公司迫切需要對員工安全積分機(jī)制、外協(xié)隊(duì)伍評(píng)價(jià)和安全承載能力分析進(jìn)行深入的研究,建立覆蓋所有專業(yè)和管理人員、符合實(shí)際管理需要的安全承載能力分析模型算法和員工安全積分標(biāo)準(zhǔn)、外協(xié)隊(duì)伍評(píng)價(jià)標(biāo)準(zhǔn),充分利用公司業(yè)務(wù)系統(tǒng)的現(xiàn)有數(shù)據(jù),實(shí)現(xiàn)智能計(jì)算匯總,進(jìn)一步減輕基層班組的負(fù)擔(dān),實(shí)現(xiàn)安全積分、外協(xié)隊(duì)伍評(píng)價(jià)和承載能力分析的智能化計(jì)算和決策分析。

        目前,國內(nèi)外關(guān)于電力企業(yè)安全體系、風(fēng)險(xiǎn)評(píng)估方面的研究有云模型法[1],基于層次分析法(AHP)[2]和故障樹法(FTA)[3]等。傳統(tǒng)的系統(tǒng)安全分析方法在風(fēng)險(xiǎn)評(píng)估中存在一定局限性[4],譬如FTA假設(shè)事件是二值的、相互獨(dú)立的,不能很好地解決復(fù)雜系統(tǒng)的建模問題,計(jì)算精確度低并且浪費(fèi)時(shí)間[5]。此外,基于員工安全積分和外協(xié)隊(duì)伍評(píng)價(jià)的安全承載能力分析的研究報(bào)道非常少。

        貝葉斯網(wǎng)絡(luò)是人工智能、概率理論、圖論、決策分析相結(jié)合的產(chǎn)物,是目前不確定知識(shí)表達(dá)和推理領(lǐng)域中最有效的理論模型之一[6]。近幾年來,貝葉斯方法廣泛應(yīng)用于各個(gè)領(lǐng)域的安全評(píng)估分析[7-9]。用貝葉斯網(wǎng)絡(luò)能很好地彌補(bǔ)傳統(tǒng)安全評(píng)估方法的不足[10]。

        本文將貝葉斯網(wǎng)絡(luò)方法應(yīng)用于基于員工安全積分和外協(xié)隊(duì)伍評(píng)價(jià)的安全承載能力分析,實(shí)現(xiàn)安全承載能力的評(píng)估。

        1 貝葉斯網(wǎng)絡(luò)的基本理論

        1.1 貝葉斯定理

        貝葉斯定理是概率統(tǒng)計(jì)中的應(yīng)用所觀察到的現(xiàn)象對先驗(yàn)概率進(jìn)行修正的標(biāo)準(zhǔn)方法,即當(dāng)分析樣本大到趨近總體數(shù)時(shí),樣本中事件發(fā)生的概率將接近于總體中事件發(fā)生的概率[11]。

        先驗(yàn)概率和后驗(yàn)概率這兩個(gè)概念是相對于某組證據(jù)而言。設(shè)H和E為兩個(gè)隨機(jī)變量,H=h為某一假設(shè),E=e為一組證據(jù),在考慮證據(jù)E=e之前,對事件H=h的概率估計(jì)P(H=h) 稱為先驗(yàn)概率,而在考慮證據(jù)之后,對H=h的概率估計(jì)P(H=h|E=e)稱為后驗(yàn)概率,貝葉斯公式描述了先驗(yàn)概率和后驗(yàn)概率之間的關(guān)系[11]:

        (1)

        1.2 貝葉斯網(wǎng)絡(luò)概述

        圖1 一個(gè)簡單的貝葉斯網(wǎng)絡(luò)示例Fig.1 An example of Bayesian network

        貝葉斯網(wǎng)絡(luò)[12]是一種有向圖模型,在貝葉斯網(wǎng)絡(luò)中,每一個(gè)原因變量和結(jié)果變量都由節(jié)點(diǎn)表示,并用概率表示變量之間的關(guān)系強(qiáng)弱,每個(gè)節(jié)點(diǎn)都有自己的節(jié)點(diǎn)概率分布以表征父節(jié)點(diǎn)對子節(jié)點(diǎn)產(chǎn)生的影響。通過貝葉斯網(wǎng)絡(luò)可以匯總各種數(shù)據(jù)并對這些數(shù)據(jù)進(jìn)行綜合推理。貝葉斯網(wǎng)絡(luò)具有如下特征圖形結(jié)構(gòu):網(wǎng)絡(luò)節(jié)點(diǎn)代表隨機(jī)變量,節(jié)點(diǎn)間的有向連線表示隨機(jī)變量之間的條件依賴關(guān)系,由節(jié)點(diǎn)和連線組成有向無環(huán)圖,節(jié)點(diǎn)條件概率表表征該節(jié)點(diǎn)相對于其父節(jié)點(diǎn)的所有可能的條件概率[13]。如圖1所示的貝葉斯網(wǎng)絡(luò),C1事件與C2事件為C3事件的父節(jié)點(diǎn), C3為C1與C2的子節(jié)點(diǎn),S表示安全狀態(tài),R表示危險(xiǎn)狀態(tài)。根據(jù)貝葉斯規(guī)則可推算C1和C2為R態(tài)的概率分別為54.2%和26.6%,因此C1是使C3處于R態(tài)的主要影響因素[13]。

        2 基于貝葉斯網(wǎng)絡(luò)的安全承載能力評(píng)估

        貝葉斯網(wǎng)絡(luò)通過實(shí)踐積累可對網(wǎng)絡(luò)結(jié)構(gòu)參數(shù)進(jìn)行更新改進(jìn),本文用貝葉斯網(wǎng)絡(luò)實(shí)現(xiàn)電力企業(yè)安全承載能力應(yīng)用過程如圖2所示。

        步驟1:要素分類識(shí)別

        此階段,對電力施工安全承載有影響的風(fēng)險(xiǎn)要素進(jìn)行歸類分級(jí),根據(jù)歷史數(shù)據(jù)、專家經(jīng)驗(yàn),調(diào)查分析等列出重要的風(fēng)險(xiǎn)要素。

        步驟2:風(fēng)險(xiǎn)要素評(píng)估

        此階段,步驟1的每項(xiàng)風(fēng)險(xiǎn)要素進(jìn)行定級(jí),確定風(fēng)險(xiǎn)要素等級(jí)數(shù)據(jù)集便于貝葉斯網(wǎng)絡(luò)的分析。風(fēng)險(xiǎn)要素等級(jí)可由式(2)表示:

        Rrisk=(DL)×(PO) ,

        (2)

        式中,Rrisk為風(fēng)險(xiǎn)要素等級(jí),DL為損失程序,PO為發(fā)生概率。風(fēng)險(xiǎn)要素等級(jí)數(shù)據(jù)集通過圖3的風(fēng)險(xiǎn)要素等級(jí)矩陣進(jìn)行規(guī)范處理后用于貝葉斯網(wǎng)絡(luò)。

        圖2 基于貝葉斯網(wǎng)絡(luò)的安全承載能力評(píng)估過程圖Fig.2 Bayesian network based safety bearing capacity evaluation procedure

        圖3 風(fēng)險(xiǎn)要素等級(jí)矩陣圖Fig.3 Rank matrix for risk factors

        步驟3:貝葉斯網(wǎng)絡(luò)構(gòu)建

        貝葉斯網(wǎng)絡(luò)通過結(jié)構(gòu)性學(xué)習(xí)進(jìn)行構(gòu)建,并用來檢驗(yàn)風(fēng)險(xiǎn)要素之間的關(guān)系。每個(gè)風(fēng)險(xiǎn)要素的可能性可通過參數(shù)學(xué)習(xí)進(jìn)行計(jì)算。設(shè)安全生產(chǎn)指標(biāo)為en,n=1,2,3,…,風(fēng)險(xiǎn)因素H發(fā)生的概率為Pn,則有P(H|e1=P1),P(H|e2=P2),…,P(H|en=Pn),則該風(fēng)險(xiǎn)因素H的發(fā)生概率為P(H|a1,a2,…,an),由于各安全度指標(biāo)之間相互獨(dú)立,則有:

        (3)

        其中,P(H)為H的風(fēng)險(xiǎn)發(fā)生概率。

        通過調(diào)查問卷、歷史數(shù)據(jù)、專家討論,在考慮安全度指標(biāo)體系的情況下,確定貝葉斯網(wǎng)絡(luò)中所有證據(jù)節(jié)點(diǎn)的概率。結(jié)合公式(3),得到安全生產(chǎn)風(fēng)險(xiǎn)概率。

        步驟4:要素控制

        對影響安全生產(chǎn)的風(fēng)險(xiǎn)要素通過貝葉斯網(wǎng)絡(luò)敏感性分析進(jìn)行評(píng)估,對風(fēng)險(xiǎn)較高的要素,采取措施進(jìn)行控制。

        步驟5:要素再評(píng)估

        步驟4中需要控制的風(fēng)險(xiǎn)要素確定后,有些風(fēng)險(xiǎn)要素的改變會(huì)導(dǎo)致電力施工安全承載能力風(fēng)險(xiǎn)的改變,因此需要更新貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)參數(shù),并進(jìn)行風(fēng)險(xiǎn)要素再評(píng)估。

        3 安全承載能力評(píng)估應(yīng)用

        3.1 風(fēng)險(xiǎn)要素分類識(shí)別和評(píng)估

        結(jié)合電力施工的歷史數(shù)據(jù)和參考文獻(xiàn),分析在電力施工過程中存在的風(fēng)險(xiǎn)要素,通過專家討論法確定風(fēng)險(xiǎn)要素和類別,并且明確風(fēng)險(xiǎn)要素發(fā)生的可能性和嚴(yán)重程度,根據(jù)式(2)和圖3中的風(fēng)險(xiǎn)要素等級(jí)矩陣確定主要風(fēng)險(xiǎn)因素和風(fēng)險(xiǎn)等級(jí)(主要風(fēng)險(xiǎn)因素只有子節(jié)點(diǎn)而沒有父節(jié)點(diǎn)),將高風(fēng)險(xiǎn)要素(R3)確定為主要風(fēng)險(xiǎn)要素。得到的電力企業(yè)安全生產(chǎn)的主要風(fēng)險(xiǎn)因素如表1。

        表1 電力施工安全承載的主要風(fēng)險(xiǎn)要素

        3.2 貝葉斯網(wǎng)絡(luò)構(gòu)建與計(jì)算

        貝葉斯網(wǎng)絡(luò)的構(gòu)建通過貝葉斯網(wǎng)絡(luò)分析軟件GeNIe Ver2.0實(shí)現(xiàn),將表1中電力企業(yè)電力施工安全承載的主要風(fēng)險(xiǎn)要素作為貝葉斯網(wǎng)絡(luò)基本節(jié)點(diǎn),各節(jié)點(diǎn)間的因果關(guān)系和條件概率分布通過歷史數(shù)據(jù)、參考文獻(xiàn)和專家討論共同確定。最終的網(wǎng)絡(luò)結(jié)構(gòu)如圖4所示。

        圖4 貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.4 The structure diagram of a Bayesian network

        與此同時(shí),將根據(jù)式(2)和圖3確定的主要風(fēng)險(xiǎn)要素的風(fēng)險(xiǎn)等級(jí)數(shù)據(jù)輸入貝葉斯網(wǎng)絡(luò)模型中,利用貝葉斯網(wǎng)絡(luò)分析軟件參數(shù)學(xué)習(xí)功能對其他節(jié)點(diǎn)(含有父節(jié)點(diǎn))進(jìn)行概率計(jì)算,得出剩余節(jié)點(diǎn)的邊緣概率分布,完成對各風(fēng)險(xiǎn)狀態(tài)的評(píng)估。將歷史數(shù)據(jù)及后續(xù)的實(shí)踐數(shù)據(jù)對模型進(jìn)行訓(xùn)練并反復(fù)修正,最終確定網(wǎng)絡(luò)結(jié)構(gòu)各節(jié)點(diǎn)參數(shù)。對電力施工安全承載能力進(jìn)行評(píng)估,得到電力施工安全承載能力處于危險(xiǎn)狀態(tài)(R)的概率為19%,如圖5所示。

        圖5 貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)分析結(jié)果Fig.5 The analysis result of the Bayesian network for safety bearing capacity evaluation

        3.3 風(fēng)險(xiǎn)再評(píng)估

        貝葉斯網(wǎng)絡(luò)具有自動(dòng)更新的功能,電力施工安全承載能力風(fēng)險(xiǎn)要素發(fā)生變化時(shí),需要及時(shí)更新整個(gè)網(wǎng)絡(luò)結(jié)構(gòu)參數(shù),進(jìn)行風(fēng)險(xiǎn)再評(píng)估,使得網(wǎng)絡(luò)中變量的概率分布發(fā)生相應(yīng)的改變。例如,在電力施工過程中選擇施工經(jīng)驗(yàn)豐富的隊(duì)伍進(jìn)行施工,即A5要素R變?yōu)?,輸入證據(jù)后,得到更新后的變量的邊緣概率,最終電力施工安全承載能力處于危險(xiǎn)狀態(tài)(R)的概率為17%,即換成經(jīng)驗(yàn)足夠豐富的隊(duì)伍進(jìn)行施工,安全承載能力得到了提高,如圖6所示。

        圖6 A5要素R為0時(shí)的葉斯網(wǎng)絡(luò)結(jié)構(gòu)分析結(jié)果Fig.6 The analysis result of the Bayesian network when the factor R of A5 is zero

        4 系統(tǒng)軟件實(shí)現(xiàn)

        圖7 電力施工安全承載能力評(píng)估系統(tǒng)總體架構(gòu)圖Fig.7 The overall architecture of safety bearing scapacity evaluation of power construction

        電力施工安全承載能力評(píng)估系統(tǒng)總體架構(gòu)圖如圖7所示,包括應(yīng)用與服務(wù)層、傳輸層、數(shù)據(jù)層、OS硬件層。安全承載能力評(píng)估的主要功能是在基于貝葉斯網(wǎng)絡(luò)模型的軟件算法對風(fēng)險(xiǎn)源進(jìn)行辨識(shí)的基礎(chǔ)上進(jìn)行風(fēng)險(xiǎn)評(píng)估、衡量風(fēng)險(xiǎn)發(fā)生的概率以及其可能造成的損失,最終確定安全承載能力。

        電力施工安全承載能力評(píng)估系統(tǒng)所用的貝葉斯網(wǎng)絡(luò)數(shù)學(xué)模型采用GeNIe軟件輔助實(shí)現(xiàn),GeNIe 軟件中具備SMILE接口。SMILE庫函數(shù)是關(guān)于決策理論模型和圖形概率、貝葉斯網(wǎng)絡(luò)等影響圖和結(jié)構(gòu)方程模型的接口,允許創(chuàng)建、編輯、保存和加載圖形化模型,并應(yīng)用在概率推理和決策的不確定性中[14]?;赟MILE接口的輸入與輸出流程圖如圖8所示。將參數(shù)值,即第一層次風(fēng)險(xiǎn)因素概率輸入到SMILE接口通過GeNIe軟件中的模型進(jìn)行仿真,將輸出結(jié)果保存為xdsl文件。

        系統(tǒng)軟件界面如圖9所示。在安全承載能力評(píng)估界面,輸入各相關(guān)風(fēng)險(xiǎn)要素級(jí)數(shù)后,點(diǎn)擊計(jì)算時(shí)系統(tǒng)會(huì)調(diào)用事先建立的xdsl文件格式的貝葉斯網(wǎng)絡(luò)模型,并將各文本框中的概率通過SMILE接口輸入到模型中進(jìn)行計(jì)算,將最終的安全承載能力計(jì)算顯示出來。

        圖8 輸入輸出流程圖Fig.8 The input and output flow chart of SMILE interface based GeNIe software

        圖9 系統(tǒng)軟件實(shí)物圖Fig.9 The software interface of evaluation system

        5 結(jié)束語

        本文運(yùn)用貝葉斯網(wǎng)絡(luò)模型實(shí)現(xiàn)了對電力施工安全承載能力的評(píng)估,并可根據(jù)實(shí)際情況對模型參數(shù)進(jìn)行更新。該系統(tǒng)為安全承載能力評(píng)估提供了可量化數(shù)據(jù)支撐和動(dòng)態(tài)自適應(yīng)功能,提高了安全承載能力評(píng)估的科學(xué)性和適應(yīng)性,減輕基層負(fù)擔(dān),提升安全承載能力評(píng)估對基層單位安全生產(chǎn)和現(xiàn)場作業(yè)的指導(dǎo)作用,確?,F(xiàn)場作業(yè)安全。本軟件系統(tǒng)目前已在本單位進(jìn)行了實(shí)際應(yīng)用,反應(yīng)效果良好。

        [1]趙徽.電網(wǎng)企業(yè)安全文化對安全行為影響的評(píng)價(jià)研究[D].華北電力大學(xué)(北京),2010.

        [2]YOU D H, CHEN Q Q, YIN X G. A study of electrical security risk assessment system based on electricity regulation[J]. Energy Policy, 2011, 39(4): 2062-2074.

        [3]YAN F Q, WANG Y, YIN X G. Research on the security risk assessment system for regional complex power grid[C]. 5th International Conference on Electric Utility Deregulation and Restructuring and Power Technologies (DRPT). Changsha, 2015,1-4.

        [4]王開碧,申曉留,李哲,等.電力企業(yè)安全生產(chǎn)風(fēng)險(xiǎn)管理方法的研究[J].黑龍江電力,2012,34 (4):228-291.

        [5]姚錫文,許開立,王貝貝,等.石化裝置火災(zāi)爆炸貝葉斯網(wǎng)絡(luò)與防護(hù)層集成分析術(shù)[J].中國安全生產(chǎn)科學(xué)技術(shù),2014,10(3):95-100.

        [6]車玉龍,蘇宏升.基于貝葉斯網(wǎng)絡(luò)的列控系統(tǒng)安全風(fēng)險(xiǎn)評(píng)估方法[J].計(jì)算機(jī)工程與應(yīng)用,2013,49(24):238-242.

        [7]李曉東,孫悅,2010年我國房屋市政工程生產(chǎn)安全事故分析[J].土木工程學(xué)報(bào),2011,11(S1):225-229.

        [8]李天縱, 王強(qiáng).一種基于貝葉斯網(wǎng)絡(luò)的軟件項(xiàng)目風(fēng)險(xiǎn)管理方法[J].計(jì)算機(jī)系統(tǒng)應(yīng)用,2011,20(2) :226-229.

        [9]任雪利.粗糙集和貝葉斯網(wǎng)絡(luò)在軟件風(fēng)險(xiǎn)評(píng)估中的應(yīng)用[J].計(jì)算機(jī)系統(tǒng)應(yīng)用,2010,19(5) :202-204.

        [10]郭發(fā)蔚,王宏輝.基于Bayesian隧道施工風(fēng)險(xiǎn)模糊綜合評(píng)估方法[J].鐵道科學(xué)與工程學(xué)報(bào),2016,13(2):401-406.

        [11]JONES B, JENKINSON I, YANG Z,et al.The use of Bayesian network modeling for maintenance planning in a manufacturing industry[J]. Reliability Engineering and System Safety,2010,95(3):267-277.

        [12]GREENBERG R, COOK S C, HARRIS D. A civil aviation safety assessment model using a Bayesian belief network (BBN) [J].The Aeronautical Journal, 2005,109(1101):557-568.

        [13]李典慶,鄢麗麗,邵東國.基于貝葉斯網(wǎng)絡(luò)的土石壩可靠性分析[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版),2007,40(6):24-29.

        [14]葛運(yùn)飛.GeNIe的擴(kuò)展研究[D].昆明:云南大學(xué),2012

        Bayesian network based evaluation of the safety bearing
        capacity of power construction

        LAI Jun1, CHEN Lei2

        (1.Huzhou Power Supply Company, Zhejiang Electric Power Co.Ltd, State Grid, Huzhou 313000, China;2.Zhejiang Electric Power Co.Ltd, State Grid, Hangzhou 310007, China)

        ∶The Bayesian network based dynamic evaluating model for safety bearing capacity was put forward.Bayesian network structure and parameters were built based on the identification and assessment of risk factors which may have effects on the safety bearing capacity. Furthermore, the risk assessment and risk diagnosis of the safety bearing capacity of power construction were carried out. Finally, the risk was controlled based on the analytical results. This system can meet the needs of evaluating the safety bearing capacity of power construction.

        ∶Bayesian network; safety bearing capacity; power construction; risk assessment

        10.3976/j.issn.1002-4026.2017.03.019

        2016-08-26

        浙江省電力公司科技計(jì)劃(5211UZ15007V)

        來駿(1981—),男,工學(xué)碩士,研究方向?yàn)殡娏ο到y(tǒng)自動(dòng)化和電力安全管理。E-mail:laijun-huzhou@qq.com

        TM76

        A

        1002-4026(2017)03-0110-07

        猜你喜歡
        貝葉斯概率要素
        第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
        第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
        概率與統(tǒng)計(jì)(一)
        概率與統(tǒng)計(jì)(二)
        掌握這6點(diǎn)要素,讓肥水更高效
        觀賞植物的色彩要素在家居設(shè)計(jì)中的應(yīng)用
        論美術(shù)中“七大要素”的辯證關(guān)系
        貝葉斯公式及其應(yīng)用
        基于貝葉斯估計(jì)的軌道占用識(shí)別方法
        也談做人的要素
        山東青年(2016年2期)2016-02-28 14:25:36
        亚洲成在人线天堂网站| 天天躁日日躁狠狠躁av中文| 欧美一级欧美一级在线播放| avtt一区| 色婷婷一区二区三区久久亚洲| 国产精品天天看天天狠| 亚洲成a v人片在线观看| 亚洲综合色一区二区三区另类| 日本高清一区二区三区视频| 国产一区二区三区日韩在线观看| 无码熟妇人妻av在线影片最多| 亚洲综合在线观看一区二区三区 | 亚洲AV秘 片一区二区三| 岛国熟女一区二区三区| 亚洲中文字幕高清在线视频一区 | 性色欲情网站| 四虎影视4hu4虎成人| 亚洲AV毛片无码成人区httP | 亚洲精品无码久久毛片| av无码天一区二区一三区| 国产目拍亚洲精品二区| 一二三区无线乱码中文在线| 女人喷潮完整视频| 亚洲女同精品一区二区久久| 日本一区不卡在线观看| 精品卡一卡二卡3卡高清乱码| 琪琪的色原网站| 国产精品九九九无码喷水| 东京热东京道日韩av| 亚洲一区二区三区地址| 国产精品久久久久9999赢消| 久久国产精品不只是精品| 国产 无码 日韩| 国产av一区二区日夜精品剧情 | 日本在线视频网站www色下载| 大岛优香中文av在线字幕| 国产成人综合久久久久久 | 国产乱人视频在线看| 国产在线精彩自拍视频| 日韩在线 | 中文| 久激情内射婷内射蜜桃人妖|