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

    基于集對分析的土石壩施工期風(fēng)險評價

    2016-12-20 07:24:42李宗坤楊朝霞
    水力發(fā)電 2016年9期
    關(guān)鍵詞:施工期石壩置信度

    李宗坤,李 奇,楊朝霞

    (鄭州大學(xué)水利與環(huán)境學(xué)院,河南鄭州450001)

    ?

    基于集對分析的土石壩施工期風(fēng)險評價

    李宗坤,李 奇,楊朝霞

    (鄭州大學(xué)水利與環(huán)境學(xué)院,河南鄭州450001)

    針對土石壩施工期安全穩(wěn)定性差,風(fēng)險事件多發(fā)的特點,引入集對分析理論,建立聯(lián)系度評價矩陣,從同、異、反3方面全面刻畫評價指標(biāo)與各風(fēng)險等級的同一和對立程度。同時考慮傳統(tǒng)AHP方法存在的專家打分偏于絕對、專家打分的可靠性及專家間動態(tài)差異未被充分考慮的不足,運用區(qū)間數(shù)學(xué)理論建立不確定型判斷矩陣,計算專家置信度并確定風(fēng)險評價指標(biāo)的動態(tài)權(quán)重,進而建立基于集對分析的土石壩施工期風(fēng)險評價模型。將此方法應(yīng)用于工程實例,結(jié)果表明,所建立的土石壩施工期風(fēng)險評價的集對分析方法可行而有效。

    土石壩;風(fēng)險評價;集對分析;不確定型AHP;動態(tài)權(quán)重

    0 引 言

    土石壩施工期一般具有施工周期長、作業(yè)環(huán)境差、受不確定因素影響大的特點,是各類風(fēng)險事件發(fā)生的重要時期。風(fēng)險事件的發(fā)生不僅會延長施工周期,增加工程投資,還會導(dǎo)致各種質(zhì)量事故,為土石壩正常運行埋下隱患。據(jù)統(tǒng)計[1],28%的土石壩失事事故發(fā)生在施工期。因此,對土石壩施工期進行風(fēng)險評價具有重要意義。

    目前,土石壩風(fēng)險分析主要有層次分析法、故障樹法、蒙特卡羅法、模糊綜合評價法等[2]。針對土石壩施工期風(fēng)險評價的研究也取得了一定進展,葛巍等[3]構(gòu)建了WBS-RBS判斷矩陣進行風(fēng)險辨識并運用AHP方法計算土石壩施工期風(fēng)險;李宗坤等[4]采用S形曲線分布作為擬合函數(shù)對常規(guī)突變評價法進行改進,并對土石壩施工期風(fēng)險進行評價;姜樹海等[5]運用Bayes模型,實現(xiàn)對土石壩安全等級動態(tài)概率的評定。然而土石壩施工期的不確定因素復(fù)雜繁多,對土石壩施工期的風(fēng)險評價方法有待進一步研究完善。另外,風(fēng)險指標(biāo)權(quán)重的確定是土石壩施工期風(fēng)險評價的重要組成部分?,F(xiàn)有研究多采用傳統(tǒng)AHP方法確定權(quán)重,但該方法未充分考慮專家自身置信度的差異且用確定分值對兩兩風(fēng)險指標(biāo)重要程度進行比較,使得計算得到的權(quán)重較難反映評價指標(biāo)重要性的客觀程度。

    集對分析作為處理不確定性問題的系統(tǒng)理論方法,明確了評價指標(biāo)在各個評價級別都具有同異反聯(lián)系的客觀事實,在水質(zhì)評價、信息系統(tǒng)、水資源預(yù)測等領(lǐng)域得到了廣泛應(yīng)用[6- 9]。本文基于集對分析理論,建立聯(lián)系度評價矩陣,從同、異、反三方面全面分析各風(fēng)險指標(biāo)值與各風(fēng)險等級的統(tǒng)一和對立程度。同時,引入?yún)^(qū)間數(shù)學(xué)理論,考慮專家對風(fēng)險指標(biāo)認(rèn)知的局限性和專家之間置信度的動態(tài)差異,提出基于不確定型AHP的風(fēng)險指標(biāo)權(quán)重確定方法,從而建立土石壩施工期風(fēng)險評價的集對分析模型,并開展應(yīng)用研究。

    1 指標(biāo)綜合權(quán)重的確定

    風(fēng)險指標(biāo)權(quán)重的確定不僅是風(fēng)險評價的重要步驟,而且對于評價結(jié)果的客觀性具有重要意義。由于同一專家對不同風(fēng)險指標(biāo)的認(rèn)知具有局限性,且不同專家的知識儲備及掌握的工程信息等也具有差異性[10],故本文引入?yún)^(qū)間數(shù)學(xué)理論,首先用區(qū)間數(shù)代替具體值來表示專家判斷,建立不確定型判斷矩陣;其次,通過計算不確定型判斷矩陣的相似度和差異度來確定專家置信度,進而求得風(fēng)險評價指標(biāo)的綜合權(quán)重,使風(fēng)險評價指標(biāo)的權(quán)重更加客觀。

    1.1 專家置信度的確定

    首先,計算判斷矩陣的相似度。如式(1)~(3)所示:

    (1)

    (2)

    (3)

    式中,As,At表示任兩個判斷矩陣;s,t取值(1,2,…,q);γst=cosαst指vec(As)與vec(At)的夾角,反映As與At的相似度;γt指At與其他判斷矩陣的相似度之和;λt指對γt歸一化處理后得到的At與其他判斷矩陣的相似度。

    其次,計算不確定型判斷矩陣的差異度。如式(4)~(5)所示:

    (4)

    (5)

    式中,atj表示不確定型判斷矩陣的元素;σt表示At與其他判斷矩陣的差異度總和;δt表示歸一化處理后,At與其他判斷矩陣的差異度。

    最后,計算專家置信度。即

    (6)

    式中,rt表示第t位專家的專家置信度;λt和δt分別表示歸一化處理后,判斷矩陣At與其他判斷矩陣的相似度和差異度。

    1.2 權(quán)重區(qū)間的確定

    (7)

    (8)

    令A(yù)-=(a-ij)n×n,A+=(a+ij)n×n則可求得:

    (9)

    則權(quán)重向量區(qū)間為ω=[lω-,mω+]=[u1,u2]。

    1.3 指標(biāo)權(quán)重的確定

    (1)指標(biāo)主觀權(quán)重為

    (10)

    (2)指標(biāo)客觀權(quán)重為

    (11)

    式中,ωi2表示第i項評價指標(biāo)客觀權(quán)重;gi反映了權(quán)重ωi1的可靠度,其值越小說明不確定型判斷矩陣的相對可靠性就越高,反之則越低;bdi為計算中間變量。

    (3)指標(biāo)綜合權(quán)重為

    (12)

    2 土石壩施工期風(fēng)險集對分析模型

    2.1 集對分析基本原理

    集對分析(Set Pair Analysis,簡稱SPA)由我國學(xué)者趙克勤[11]于1989年首次提出。該方法將具有確定性和不確定性信息的兩個集合A和B構(gòu)成集對H=(A,B),并用“同一”和“對立”來描述確定性信息,用“差異”來描述不確定性信息。通過聯(lián)系度μ對集對中兩集合的特性從同、異、反3方面進行定量刻畫,其中聯(lián)系度μ的表達式為

    (13)

    式中,μ表示同、異、反聯(lián)系度;a、b、c分別表示集合A和B的同一度、差異度和對立度,a,b,c∈[0,1],且a+b+c=1;N表示集合的總特性數(shù),P,S分別為兩集合的對立特性個數(shù)和共有特性個數(shù),F(xiàn)=N-S-P;i為差異度系數(shù),且i∈[-1,1];j為對立度系數(shù),通常取j=-1。

    2.2 土石壩施工期風(fēng)險評價指標(biāo)體系

    建立基于集對分析的土石壩施工期風(fēng)險評價模型,首先必須對土石壩施工期風(fēng)險進行辨識,并建立土石壩施工期風(fēng)險評價指標(biāo)體系,如圖1所示。該體系包含3個一級風(fēng)險指標(biāo)Ri(i=1,2,3)和11個二級風(fēng)險指標(biāo)Rij(j=1,2,3,4,5,6,7或1,2,3),該評價體系可較好的反映土石壩施工期的風(fēng)險因素。

    圖1 土石壩施工期風(fēng)險評價指標(biāo)體系

    在上述評價體系基礎(chǔ)上,還需建立合理的風(fēng)險分級標(biāo)準(zhǔn)。本文根據(jù)目前風(fēng)險分析中較為常用的五級分類法[12- 13],建立土石壩施工期風(fēng)險等級的劃分標(biāo)準(zhǔn)及相應(yīng)接受準(zhǔn)則,并根據(jù)現(xiàn)有研究的量化標(biāo)準(zhǔn)[14],將風(fēng)險等級進行量化處理,如表1所示。

    表1 風(fēng)險分析評價標(biāo)準(zhǔn)

    風(fēng)險等級定性描述定量描述Ⅰ風(fēng)險極低,可忽略,不需要采取風(fēng)險處理措施和檢測(90,100]Ⅱ風(fēng)險較低,可容許,不需要采取風(fēng)險處理措施,但需引起注意,進行常規(guī)管理(80,90]Ⅲ中等風(fēng)險,可接受,但需要引起重視,設(shè)計上必須考慮,施工時要制定詳細(xì)的管理計劃(70,80]Ⅳ高風(fēng)險,部分可接受,風(fēng)險事故一旦發(fā)生對整個工程造成較為嚴(yán)重的影響,必須采取風(fēng)險處理措施盡可能將風(fēng)險降低至三級,還需要加強監(jiān)測并建立預(yù)警措施(60,70]Ⅴ風(fēng)險極高,不可接受,必須高度重視并采取措施進行規(guī)避和轉(zhuǎn)移,否則要不惜代價將風(fēng)險至少降低至四級[50,60]

    2.3 集對分析計算模型

    建立基于集對分析的土石壩施工期風(fēng)險評價模型。從同、異、反三方面全面分析各風(fēng)險指標(biāo)與各風(fēng)險等級的同一和對立程度,進而確定施工期整體風(fēng)險等級,具體步驟為

    (1)集對模型的構(gòu)建。設(shè)專家對各二級風(fēng)險評價指標(biāo)的評分集合為Qij,風(fēng)險評價標(biāo)準(zhǔn)集合為Xkj,根據(jù)集對分析理論,建立兩集合的集對H(Qij,Xkj),其中:

    (14)

    (2)二級評價矩陣的確定。①計算各二級風(fēng)險指標(biāo)對5個風(fēng)險等級的聯(lián)系度μⅠ、μⅡ、μⅢ、μⅣ、μⅤ,并確定各二級風(fēng)險指標(biāo)的聯(lián)系度向量(μijⅠ,μijⅡ,μijⅢ,μijⅣ,μijⅤ),即

    (15)

    (16)

    (17)

    (18)

    (19)

    ②確定二級評價矩陣。設(shè)Ui(i=1,2,3)表示3個一級風(fēng)險指標(biāo)Ri分別對應(yīng)的二級風(fēng)險指標(biāo)Rij(j=1,2,3,4,5,6,7或1,2)的聯(lián)系度向量所構(gòu)成的二級評價矩陣,即

    (20)

    例如,U1表示一級風(fēng)險指標(biāo)R1所包含的7個二級風(fēng)險指標(biāo)的聯(lián)系度向量所構(gòu)成的7×5階矩陣。

    (3)風(fēng)險等級的確定。根據(jù)以上計算內(nèi)容,建立土石壩施工期風(fēng)險評價計算模型,即

    (21)

    式中,Ui表示二級評價矩陣;Wi表示各級風(fēng)險指標(biāo)權(quán)向量;Ni表示二級評價結(jié)果向量,即3個一級風(fēng)險指標(biāo)所對應(yīng)的二級風(fēng)險指標(biāo)的綜合評價結(jié)果向量;N0表示3個一級風(fēng)險指標(biāo)綜合評價結(jié)果向量,即土石壩施工期風(fēng)險等級的聯(lián)系度向量;“o”表示模糊乘運算算子。

    3 實例分析

    表2 二級風(fēng)險指標(biāo)評分均值

    二級風(fēng)險指標(biāo)專家評分(qij)方案一方案二方案三二級專家方案一方案二方案三R11899196R21848487R12848484R22868591R13968491R23898989R14939393R31898689R15969696R32878484R16918988R33919196

    表3 一級風(fēng)險指標(biāo)權(quán)重

    專家專家置信度一級風(fēng)險指標(biāo)權(quán)重區(qū)間R1R2R3A01975[03035,03132][03432,03575][03348,03478]B02802[04065,04258][02913,03085][02769,02946]C02452[03852,03985][03040,03157][02934,03031]D02771[03671,03773][02925,03007][03250,03374]

    注:專家置信度下R1、R2、R3的綜合權(quán)重分別為0.382、0.310、0.308。

    表4 安全風(fēng)險二級指標(biāo)權(quán)重

    專家專家置信度安全風(fēng)險因素權(quán)重區(qū)間R11R12R13R14R15R16A02547[01030,01086][00840,00873][01472,01573][02140,02229][02104,02113][02240,02429]B02462[01056,01097][00715,00732][01506,01556][02148,02260][02136,02209][02348,02460]C02488[00778,00809][01210,01265][01295,01369][02160,02298][02079,02135][02360,02498]D02503[01010,01064][00947,00986][01407,01505][02079,02145][02109,02154][02179,02345]

    注:專家置信度下R11、R12、R13、R14、R15、R16的綜合權(quán)重分別為0.101、0.095、0.145、0.215、0.211、0.233。

    表5 工期風(fēng)險二級指標(biāo)權(quán)重

    專家專家置信度工期風(fēng)險因素權(quán)重區(qū)間R21R22R23A02515[04803,04979][03023,03149][01974,02073]B02320[04020,04313][03342,03524][02232,02468]C02617[04803,04979][03223,03349][01774,01873]D02636[04193,04502][03189,03369][02273,02473]

    注:專家置信度下R21、R22、R23的綜合權(quán)重分別為0.450、0.328、0.222。

    表6 成本風(fēng)險二級指標(biāo)權(quán)重

    專家專家置信度成本風(fēng)險因素權(quán)重區(qū)間R31R32R33A02683[03705,03922][03527,03698][02538,02609]B02673[03715,03932][03517,03688][02538,02609]C02401[03630,03938][03660,03915][02381,02477]D02244[03912,04112][03204,03425][02632,02714]

    注:專家置信度下R31、R32、R33的綜合權(quán)重分別為0.385、0.256、0.359。

    表7 權(quán)向量Wi計算結(jié)果

    權(quán)向量權(quán)向量值W0[0382,0310,0308]W1[0101,0095,0145,0215,0211,0233]W2[0450,0328,0222]W3[0385,0256,0359]

    本文僅以方案三為例進行計算,具體過程如下:

    (1)根據(jù)式(10)~(13)計算各風(fēng)險評價指標(biāo)權(quán)重,計算結(jié)果如表3~ 7所示。

    (2)根據(jù)式(1)~(8)計算各二級風(fēng)險指標(biāo)與5個風(fēng)險等級的聯(lián)系度并確定二級評價矩陣Ui,其中

    (3)由式(9)計算二級評價結(jié)果向量:N1=[0.6409+0.3401i+0019j,0.3496+0.5880i+0.0624j,0.0095+0.1787i+0.8118j,0+0i+1j,0+0i+1j];N2=[0.3756+0.6244i+0j,0.6244+0.3756i+0j,0+0.3144i+0.6856j,0+0i+1j,0+0i+1j];N3=[0.5130+0.4358i+0.0512j,0.4614+0.4668i+0.0718j,0.0256+0.3074i+0.6670j,0+0i+1j,0+0i+1j]。進而得到一級評價結(jié)果向量N0=[0.5193+0.4577i+0.0230j,0.4692+0.4848i+0.0460j,0.0115+0.2604i+0.7281j,0+0i+1j,0+0i+1j]。

    實例分析結(jié)果表明,方案三的風(fēng)險等級比方案一和方案二低,應(yīng)優(yōu)先選擇方案三。燕山水庫實際建設(shè)中也正是選擇了方案三,在按期完成工程建設(shè)任務(wù)的同時,有效保障了工程安全并合理控制了工程投資,獲得了良好的經(jīng)濟和社會效益。由此可知,根據(jù)本文方法得到的評價結(jié)果與實際相符,集對分析方法運用于土石壩施工期的風(fēng)險評價是合理可行的。

    5 結(jié) 語

    本文從安全、工期、成本3個方面建立土石壩施工期風(fēng)險評價體系,較好地反映了土石壩施工期的風(fēng)險特性。采用不確定型AHP法確定指標(biāo)權(quán)重,使得權(quán)重計算更具合理性和可操作性。建立基于集對分析的土石壩施工期風(fēng)險評價模型,從同、異、反三方面刻畫風(fēng)險指標(biāo)與各風(fēng)險等級的同一和對立程度,進而確定土石壩施工期風(fēng)險等級。將該方法應(yīng)用于燕山水庫得到的評價結(jié)果與實際相符,說明該方法有效可行,為土石壩施工期風(fēng)險評價提供了一種新的思路。

    [1]KRIMSKY S. The role of theory in risk studies[M]∥KRIMSKY S, GOLDING D, ed. Social Theories of Risk, Westport, Praeger, 1992.

    [2]麻榮永. 土石壩風(fēng)險分析方法及應(yīng)用[M]. 北京: 科學(xué)出版社, 2004.

    [3]葛巍, 李宗坤, 王文姣, 等. 基于WBS-RBS和AHP的土石壩施工期風(fēng)險評估[J]. 人民黃河, 2013, 35(6): 121- 123.

    [4]李宗坤, 葛巍, 王娟, 等. 改進的突變評價法在土石壩施工期風(fēng)險評價中的應(yīng)用[J]. 水利學(xué)報, 2014, 45(10): 1256- 1260.

    [5]姜樹海, 范子武. 應(yīng)用Bayes方法對土石壩安全等級的動態(tài)概率評定[J]. 水利學(xué)報, 2008, 39(8): 922- 926.

    [6]YUE Wencong, CAI Yangpeng, RONG Qiangqiang, et al. A hybrid life-cycle and fuzzy-set-pair analyses approach for comprehensively evaluating impacts of industrial wastewater under uncertainty[J]. Journal of Cleaner Production, 2014, 80(1): 57- 68.

    [7]李文君, 邱林, 陳曉楠, 等. 基于集對分析與可變模糊集的河流生態(tài)健康評價模型[J]. 水利學(xué)報, 2011, 42(7): 775- 782.

    [8]WU F F, WANG X. Eutrophication Evaluation Based on Set Pair Analysis of Baiyangdian Lake, North China[J]. Procedia Environmental Sciences, 2012(13): 1030- 1036.

    [9]LIU Yang. Study on Environment Evaluation and Protection Based on Set Pair Analysis—A Case Study of Chongqing [J]. Energy Procedia, 2012(14): 14- 19.

    [10]翟友成, 曹文貴, 王江營, 等. 基于不確定型層次分析法的邊坡穩(wěn)定模糊評價辦法 [J]. 巖土力學(xué), 2011, 32(2): 539- 543.

    [11]趙克勤. 集對分析及其初步應(yīng)用[M]. 杭州: 浙江科學(xué)技術(shù)出版社, 2000.

    [12]曹文貴, 翟友成, 王江營, 等. 山嶺隧道塌方風(fēng)險的集對分析方法 [J]. 中國公路學(xué)報, 2012, 25(2): 90- 98.

    [13]張國棟, 李雷, 彭雪輝. 基于大壩安全鑒定和專家經(jīng)驗的病險程度評價技術(shù)[J]. 中國安全科學(xué)學(xué)報, 2008, 18(9): 158- 166.

    [14]曹文貴, 張永杰. 地下結(jié)構(gòu)巖體質(zhì)量分類的變權(quán)重二級模糊綜合評判方法研究[J]. 巖土力學(xué)與工程學(xué)報, 2006, 25(8): 1612- 1618.

    [15]王建武, 柳鋒波, 等. 水工程施工期技術(shù)研究與應(yīng)用[M]. 北京: 中國水利水電出版社, 2011.

    (責(zé)任編輯 王 琪)

    Risk Evaluation of Earth-rock Dam during Construction Period Based on Set Pair Analysis

    LI Zongkun, LI Qi, YANG Chaoxia

    (School of Water Conservancy and Environment, Zhengzhou University, Zhengzhou 450001, Henan, China)

    Because many risk events are usually occurred in the construction period of earth-rock dam, the Set Pair Analysis theory is introduced to evaluate the risks. Firstly, the connection degree evaluation matrixes are established, which can comprehensively describe the relationships between evaluation indexes and risk ranks from three aspects of identity, difference and antagonism. And then, an uncertain type of AHP method based on expert’s credibility is proposed herein to compute experts’ confidences and determine the dynamic weights of risk factors after considering the shortages of expert scoring to be absolute, less considering the credibility of expert scoring and the dynamic difference between experts in traditional AHP evaluation. Finally, the risk evaluation model of earth-rock dam during construction period based on Set Pair Analysis method is developed. The application of model shows that the Set Pair Analysis method is feasible and effective.

    earth-rock dam; risk analysis; Set Pair Analysis; uncertain type of AHP; dynamic weight

    2016- 04- 05

    國家自然科學(xué)基金資助項目(51379192)

    李宗坤(1961—),男,河南鄧州人,教授,博士生導(dǎo)師,博士,研究方向為水利水電工程健康診斷及安全性評價研究.

    TV641

    A

    0559- 9342(2016)09- 0054- 06

    猜你喜歡
    施工期石壩置信度
    硼鋁復(fù)合材料硼含量置信度臨界安全分析研究
    正負(fù)關(guān)聯(lián)規(guī)則兩級置信度閾值設(shè)置方法
    水電站施工期超標(biāo)洪水應(yīng)對措施研究
    歡迎訂閱《碾壓式土石壩設(shè)計》
    珠三角世界級城市群建設(shè)進入施工期
    新型加筋土技術(shù)在土石壩除險加固中的應(yīng)用
    置信度條件下軸承壽命的可靠度分析
    軸承(2015年2期)2015-07-25 03:51:04
    施工期高速公路改擴建作業(yè)區(qū)分流組織研究
    淺談土石壩輸水涵管安全隱患及處理措施
    土石壩基礎(chǔ)處理的常用措施
    精品电影一区二区在线| 日本撒尿小便嘘嘘汇集6| 免费搜索国产男女视频| 久久九九热精品免费| 国产精品美女特级片免费视频播放器 | 久久人妻av系列| 国产亚洲精品久久久久久毛片| 午夜老司机福利片| 欧美日本中文国产一区发布| 亚洲伊人色综图| 亚洲欧美日韩高清在线视频| 满18在线观看网站| 日本 av在线| 神马国产精品三级电影在线观看 | 啪啪无遮挡十八禁网站| 中文字幕人妻丝袜一区二区| 午夜福利欧美成人| 欧美日韩亚洲高清精品| 国产免费av片在线观看野外av| 91麻豆精品激情在线观看国产 | 国产蜜桃级精品一区二区三区| 国产精品一区二区精品视频观看| 曰老女人黄片| 最新美女视频免费是黄的| 老汉色av国产亚洲站长工具| 黄网站色视频无遮挡免费观看| 男男h啪啪无遮挡| 国产成人av教育| 18禁裸乳无遮挡免费网站照片 | a级毛片黄视频| 国产精品一区二区免费欧美| 电影成人av| 日本a在线网址| 成人三级黄色视频| 成人av一区二区三区在线看| 国产乱人伦免费视频| 老司机福利观看| 一级,二级,三级黄色视频| 国产一卡二卡三卡精品| 久久狼人影院| 999久久久国产精品视频| 中文亚洲av片在线观看爽| 欧美激情久久久久久爽电影 | 久久性视频一级片| cao死你这个sao货| 亚洲aⅴ乱码一区二区在线播放 | 日韩av在线大香蕉| 两人在一起打扑克的视频| 国产精品一区二区在线不卡| 中文字幕高清在线视频| 日韩精品中文字幕看吧| 午夜视频精品福利| 欧美黄色片欧美黄色片| 日韩av在线大香蕉| 一进一出抽搐动态| 久久人人爽av亚洲精品天堂| 岛国在线观看网站| 午夜精品国产一区二区电影| 久久精品aⅴ一区二区三区四区| 99国产精品免费福利视频| 久久中文字幕一级| 一边摸一边抽搐一进一小说| 最新在线观看一区二区三区| 麻豆一二三区av精品| 欧美精品亚洲一区二区| 身体一侧抽搐| 99热只有精品国产| 人成视频在线观看免费观看| 久久这里只有精品19| 男女高潮啪啪啪动态图| 免费在线观看视频国产中文字幕亚洲| 国产人伦9x9x在线观看| 亚洲av成人一区二区三| 日韩免费av在线播放| 婷婷六月久久综合丁香| 免费高清在线观看日韩| 欧美日韩黄片免| 久久久久久久精品吃奶| 亚洲中文av在线| 亚洲熟女毛片儿| 国产精品99久久99久久久不卡| 搡老熟女国产l中国老女人| 午夜成年电影在线免费观看| 9色porny在线观看| 99热国产这里只有精品6| 中文字幕人妻熟女乱码| 国产欧美日韩一区二区精品| 高清毛片免费观看视频网站 | 国产无遮挡羞羞视频在线观看| 精品国产乱子伦一区二区三区| 黑人猛操日本美女一级片| 淫妇啪啪啪对白视频| 99精品在免费线老司机午夜| 水蜜桃什么品种好| 亚洲国产欧美日韩在线播放| 午夜视频精品福利| 免费看a级黄色片| 久久狼人影院| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日韩中文字幕国产精品一区二区三区 | 麻豆国产av国片精品| 国产亚洲精品久久久久久毛片| 免费观看精品视频网站| xxx96com| 国产熟女xx| 热re99久久精品国产66热6| 亚洲欧美日韩高清在线视频| 久久精品人人爽人人爽视色| www.999成人在线观看| 淫秽高清视频在线观看| 在线国产一区二区在线| 天天躁狠狠躁夜夜躁狠狠躁| 18禁观看日本| 婷婷精品国产亚洲av在线| 亚洲片人在线观看| 欧美日韩亚洲综合一区二区三区_| 久久久久久久久久久久大奶| 欧美黄色淫秽网站| 自拍欧美九色日韩亚洲蝌蚪91| 久久国产精品人妻蜜桃| 亚洲性夜色夜夜综合| 国产精品国产av在线观看| 亚洲精品一区av在线观看| 国产精品香港三级国产av潘金莲| 国产又爽黄色视频| 两个人免费观看高清视频| 免费av毛片视频| 视频区欧美日本亚洲| 老鸭窝网址在线观看| 免费看十八禁软件| 免费在线观看亚洲国产| 亚洲欧洲精品一区二区精品久久久| 51午夜福利影视在线观看| a在线观看视频网站| 美女 人体艺术 gogo| √禁漫天堂资源中文www| 色综合欧美亚洲国产小说| 国产精品一区二区在线不卡| 一级毛片高清免费大全| 超碰成人久久| 丝袜在线中文字幕| 美女大奶头视频| 国产亚洲欧美在线一区二区| 99国产精品一区二区蜜桃av| av网站免费在线观看视频| 国产精华一区二区三区| aaaaa片日本免费| 50天的宝宝边吃奶边哭怎么回事| 成人三级做爰电影| 免费不卡黄色视频| 国产在线精品亚洲第一网站| 91九色精品人成在线观看| 久久精品91蜜桃| 久久久久亚洲av毛片大全| 一进一出好大好爽视频| 丝袜美腿诱惑在线| 国产精品美女特级片免费视频播放器 | 精品国产国语对白av| 女性被躁到高潮视频| 国产一区二区三区视频了| 国产精品久久久av美女十八| 国产一区二区三区在线臀色熟女 | 国产成人精品久久二区二区免费| 女同久久另类99精品国产91| 欧美日韩亚洲高清精品| 欧美在线黄色| 日本精品一区二区三区蜜桃| 丝袜美腿诱惑在线| 久久草成人影院| 国内久久婷婷六月综合欲色啪| 亚洲在线自拍视频| av电影中文网址| 老司机在亚洲福利影院| 精品人妻1区二区| 人人妻人人爽人人添夜夜欢视频| 亚洲中文av在线| 久久久精品欧美日韩精品| 99国产精品一区二区蜜桃av| 国产精华一区二区三区| 久久热在线av| 亚洲av五月六月丁香网| 国产精品影院久久| 美女高潮到喷水免费观看| 啪啪无遮挡十八禁网站| 国产精品秋霞免费鲁丝片| tocl精华| 夫妻午夜视频| 国产亚洲精品一区二区www| av天堂在线播放| 天堂中文最新版在线下载| 久久久国产成人免费| 美女 人体艺术 gogo| 成人18禁在线播放| 热re99久久国产66热| 女同久久另类99精品国产91| 黄片小视频在线播放| 久久久久久久精品吃奶| 三级毛片av免费| 窝窝影院91人妻| 精品国产亚洲在线| 国产精品久久久av美女十八| 精品国产乱码久久久久久男人| tocl精华| 黄色丝袜av网址大全| 亚洲自偷自拍图片 自拍| 日本五十路高清| 国产一卡二卡三卡精品| 国产一区二区激情短视频| 成年人黄色毛片网站| 久久久久久亚洲精品国产蜜桃av| 亚洲自拍偷在线| 一本综合久久免费| 悠悠久久av| 精品国产一区二区三区四区第35| 日日爽夜夜爽网站| 淫秽高清视频在线观看| 亚洲av日韩精品久久久久久密| 一级a爱视频在线免费观看| www.自偷自拍.com| 国产精品二区激情视频| 久久久精品欧美日韩精品| 女同久久另类99精品国产91| 午夜精品国产一区二区电影| 久久人人97超碰香蕉20202| 免费在线观看完整版高清| 亚洲一卡2卡3卡4卡5卡精品中文| 1024视频免费在线观看| avwww免费| 欧美日韩瑟瑟在线播放| 麻豆久久精品国产亚洲av | 啦啦啦 在线观看视频| 精品国产一区二区久久| 两个人看的免费小视频| 亚洲欧美日韩另类电影网站| 丰满的人妻完整版| 十八禁人妻一区二区| 免费观看精品视频网站| 亚洲成人免费av在线播放| 亚洲黑人精品在线| 搡老乐熟女国产| 香蕉久久夜色| 91字幕亚洲| 国产一区二区在线av高清观看| 欧美日韩亚洲综合一区二区三区_| 无人区码免费观看不卡| 国产aⅴ精品一区二区三区波| 欧美最黄视频在线播放免费 | 亚洲人成电影观看| 国产精品日韩av在线免费观看 | 极品教师在线免费播放| 真人一进一出gif抽搐免费| 欧美中文综合在线视频| 免费在线观看视频国产中文字幕亚洲| 一边摸一边抽搐一进一出视频| 91成人精品电影| 女性生殖器流出的白浆| 成人手机av| 欧美日韩亚洲国产一区二区在线观看| 黑人操中国人逼视频| 韩国精品一区二区三区| 一二三四社区在线视频社区8| 亚洲欧美一区二区三区黑人| 中文字幕色久视频| 欧美乱色亚洲激情| 久久精品国产亚洲av香蕉五月| 日日干狠狠操夜夜爽| av片东京热男人的天堂| 久久久久国产精品人妻aⅴ院| 午夜a级毛片| www.999成人在线观看| 最近最新中文字幕大全免费视频| 一级,二级,三级黄色视频| 欧美日韩一级在线毛片| 看黄色毛片网站| 免费看a级黄色片| 电影成人av| 怎么达到女性高潮| 少妇粗大呻吟视频| 久久中文字幕一级| 最近最新免费中文字幕在线| 欧美激情 高清一区二区三区| 91麻豆精品激情在线观看国产 | 国产麻豆69| 亚洲中文日韩欧美视频| 91成人精品电影| 国产av一区在线观看免费| 色播在线永久视频| 男女午夜视频在线观看| 日本免费a在线| aaaaa片日本免费| 成年版毛片免费区| 亚洲久久久国产精品| 电影成人av| 久久久久久久久免费视频了| 亚洲伊人色综图| cao死你这个sao货| 淫妇啪啪啪对白视频| 日韩精品青青久久久久久| 69av精品久久久久久| 黑人欧美特级aaaaaa片| 12—13女人毛片做爰片一| 香蕉丝袜av| 成熟少妇高潮喷水视频| 久久久久久免费高清国产稀缺| 欧美黑人欧美精品刺激| 色婷婷av一区二区三区视频| 女性被躁到高潮视频| 黄色女人牲交| 女警被强在线播放| 国产乱人伦免费视频| 久久精品国产亚洲av香蕉五月| 国产亚洲精品久久久久5区| av网站在线播放免费| 亚洲人成网站在线播放欧美日韩| 亚洲av成人不卡在线观看播放网| 午夜福利,免费看| 两性夫妻黄色片| 亚洲 欧美一区二区三区| 三上悠亚av全集在线观看| 国产不卡一卡二| 91国产中文字幕| 亚洲欧美精品综合久久99| 亚洲中文字幕日韩| 亚洲久久久国产精品| 十八禁人妻一区二区| 国产激情欧美一区二区| 欧美在线黄色| 精品午夜福利视频在线观看一区| 久久久国产欧美日韩av| 三级毛片av免费| 一进一出抽搐gif免费好疼 | 777久久人妻少妇嫩草av网站| 一个人免费在线观看的高清视频| 久久青草综合色| 免费av毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 黄色女人牲交| 制服人妻中文乱码| 一区二区日韩欧美中文字幕| 美女扒开内裤让男人捅视频| 久久精品aⅴ一区二区三区四区| 国产精品电影一区二区三区| 在线永久观看黄色视频| 日韩欧美在线二视频| 一区福利在线观看| av网站在线播放免费| 一级a爱片免费观看的视频| 成人国语在线视频| 女人高潮潮喷娇喘18禁视频| 日本撒尿小便嘘嘘汇集6| 日本免费一区二区三区高清不卡 | 国产精品久久电影中文字幕| 久热这里只有精品99| 日本wwww免费看| www.自偷自拍.com| 中文字幕人妻熟女乱码| 欧美一区二区精品小视频在线| 亚洲黑人精品在线| 侵犯人妻中文字幕一二三四区| 亚洲九九香蕉| 国产精品久久久久成人av| 午夜福利,免费看| 欧美中文综合在线视频| 亚洲视频免费观看视频| 亚洲情色 制服丝袜| 黄色女人牲交| 国产高清国产精品国产三级| 777久久人妻少妇嫩草av网站| 99国产综合亚洲精品| 精品第一国产精品| 老汉色av国产亚洲站长工具| 在线观看免费高清a一片| 亚洲精品粉嫩美女一区| 亚洲av成人一区二区三| 日韩欧美免费精品| 日韩免费高清中文字幕av| 免费在线观看完整版高清| 成人影院久久| 久久久国产欧美日韩av| 操美女的视频在线观看| 一级a爱片免费观看的视频| 亚洲欧美精品综合久久99| 亚洲国产看品久久| 亚洲精品国产色婷婷电影| 老熟妇仑乱视频hdxx| 国产熟女午夜一区二区三区| 91字幕亚洲| 很黄的视频免费| 亚洲伊人色综图| 亚洲免费av在线视频| 亚洲国产毛片av蜜桃av| 桃色一区二区三区在线观看| 亚洲国产精品999在线| 老司机午夜福利在线观看视频| 日韩免费av在线播放| 久99久视频精品免费| www.自偷自拍.com| 少妇的丰满在线观看| 久久国产精品男人的天堂亚洲| 两人在一起打扑克的视频| 日韩成人在线观看一区二区三区| 老司机午夜十八禁免费视频| 神马国产精品三级电影在线观看 | 神马国产精品三级电影在线观看 | 精品国产超薄肉色丝袜足j| 成年人黄色毛片网站| 成人av一区二区三区在线看| 精品熟女少妇八av免费久了| 大型黄色视频在线免费观看| 国产亚洲欧美98| 无遮挡黄片免费观看| 婷婷丁香在线五月| 91麻豆精品激情在线观看国产 | 国产一区二区三区综合在线观看| 一级毛片精品| 黑丝袜美女国产一区| 真人做人爱边吃奶动态| 午夜亚洲福利在线播放| 亚洲专区国产一区二区| 波多野结衣一区麻豆| 日本免费一区二区三区高清不卡 | 丝袜美足系列| 俄罗斯特黄特色一大片| 757午夜福利合集在线观看| 免费搜索国产男女视频| 久久久久久免费高清国产稀缺| 久久午夜亚洲精品久久| 亚洲 国产 在线| av片东京热男人的天堂| 久久国产乱子伦精品免费另类| 夜夜夜夜夜久久久久| 国产在线观看jvid| 天堂影院成人在线观看| 国产一区二区在线av高清观看| 亚洲 欧美一区二区三区| 午夜福利在线免费观看网站| www日本在线高清视频| 一级,二级,三级黄色视频| 女同久久另类99精品国产91| tocl精华| 久久久久亚洲av毛片大全| 国产精品亚洲av一区麻豆| 午夜a级毛片| 久久草成人影院| 乱人伦中国视频| 日本欧美视频一区| 男女午夜视频在线观看| 51午夜福利影视在线观看| 看免费av毛片| 天堂动漫精品| 最新美女视频免费是黄的| 国产三级在线视频| aaaaa片日本免费| 在线观看免费视频日本深夜| 狠狠狠狠99中文字幕| 欧美日韩国产mv在线观看视频| 国产高清激情床上av| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av成人不卡在线观看播放网| 亚洲av第一区精品v没综合| 91在线观看av| 国产三级黄色录像| 黄色怎么调成土黄色| 天天躁狠狠躁夜夜躁狠狠躁| 老司机福利观看| 亚洲第一青青草原| 中文字幕高清在线视频| 高清黄色对白视频在线免费看| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲精品久久久久5区| 中文字幕精品免费在线观看视频| 男女下面插进去视频免费观看| 欧美不卡视频在线免费观看 | 成人影院久久| 在线观看一区二区三区| 人妻丰满熟妇av一区二区三区| 美女高潮到喷水免费观看| 日韩国内少妇激情av| 免费看十八禁软件| 韩国精品一区二区三区| www.精华液| 一二三四社区在线视频社区8| 一级毛片高清免费大全| 日韩av在线大香蕉| 老司机靠b影院| 亚洲成国产人片在线观看| 另类亚洲欧美激情| 国产精品电影一区二区三区| а√天堂www在线а√下载| 亚洲中文字幕日韩| 操美女的视频在线观看| 午夜免费观看网址| 99精品在免费线老司机午夜| 亚洲七黄色美女视频| 一级黄色大片毛片| 久久久精品国产亚洲av高清涩受| 97碰自拍视频| 国产亚洲欧美精品永久| 岛国视频午夜一区免费看| 9191精品国产免费久久| 电影成人av| 久久久久久人人人人人| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区日韩欧美中文字幕| 黄色毛片三级朝国网站| 51午夜福利影视在线观看| 久热爱精品视频在线9| 中文字幕人妻熟女乱码| av在线播放免费不卡| 精品久久久久久久久久免费视频 | xxx96com| 精品国产乱码久久久久久男人| 成年版毛片免费区| 99热国产这里只有精品6| 国产三级在线视频| 精品第一国产精品| 两个人免费观看高清视频| 精品一区二区三区av网在线观看| 露出奶头的视频| 黄色成人免费大全| 99热国产这里只有精品6| 无遮挡黄片免费观看| 成人18禁高潮啪啪吃奶动态图| 琪琪午夜伦伦电影理论片6080| 无限看片的www在线观看| 嫁个100分男人电影在线观看| 大型黄色视频在线免费观看| √禁漫天堂资源中文www| 长腿黑丝高跟| 一边摸一边抽搐一进一出视频| 青草久久国产| 99国产极品粉嫩在线观看| 精品午夜福利视频在线观看一区| 午夜福利,免费看| 99久久99久久久精品蜜桃| 久久天躁狠狠躁夜夜2o2o| 男人舔女人的私密视频| 国产成人啪精品午夜网站| 琪琪午夜伦伦电影理论片6080| 日韩国内少妇激情av| 后天国语完整版免费观看| 亚洲成人精品中文字幕电影 | 国产精品影院久久| 欧美大码av| 精品国产乱码久久久久久男人| 亚洲五月色婷婷综合| 女警被强在线播放| 少妇被粗大的猛进出69影院| 午夜日韩欧美国产| www.www免费av| 亚洲美女黄片视频| 黄色怎么调成土黄色| 91av网站免费观看| 精品少妇一区二区三区视频日本电影| 国产成人av教育| 一进一出好大好爽视频| 午夜福利影视在线免费观看| 1024视频免费在线观看| 最好的美女福利视频网| 国产成人av激情在线播放| 一夜夜www| 丝袜在线中文字幕| 精品欧美一区二区三区在线| 日韩三级视频一区二区三区| 777久久人妻少妇嫩草av网站| 精品熟女少妇八av免费久了| 男女做爰动态图高潮gif福利片 | 日韩精品中文字幕看吧| 18禁黄网站禁片午夜丰满| 在线永久观看黄色视频| 亚洲一区二区三区不卡视频| 久久性视频一级片| 精品国产一区二区三区四区第35| 看免费av毛片| a级毛片在线看网站| 伊人久久大香线蕉亚洲五| 国产高清videossex| 亚洲性夜色夜夜综合| 热re99久久精品国产66热6| e午夜精品久久久久久久| 亚洲成av片中文字幕在线观看| 欧美乱色亚洲激情| 麻豆国产av国片精品| 一区福利在线观看| 欧美日韩国产mv在线观看视频| 国产成人影院久久av| 丁香欧美五月| 国产精品野战在线观看 | 亚洲精品在线美女| 亚洲第一欧美日韩一区二区三区| 80岁老熟妇乱子伦牲交| 黑人操中国人逼视频| 亚洲色图 男人天堂 中文字幕| 欧美日本亚洲视频在线播放| 看免费av毛片| 视频在线观看一区二区三区| 亚洲av成人不卡在线观看播放网| 亚洲精品国产色婷婷电影| 俄罗斯特黄特色一大片| 每晚都被弄得嗷嗷叫到高潮| www.熟女人妻精品国产| 久久中文字幕人妻熟女| 亚洲欧美激情在线| 人妻久久中文字幕网| 成人免费观看视频高清| 又紧又爽又黄一区二区| 老司机午夜福利在线观看视频| 日本黄色日本黄色录像| 色婷婷久久久亚洲欧美| 一级作爱视频免费观看| 国产精品久久久人人做人人爽| 午夜91福利影院| 成人手机av| 国产又色又爽无遮挡免费看| 看黄色毛片网站|