• <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ǔ)處理的常用措施
    国产精品一区二区精品视频观看| 91精品国产国语对白视频| 老司机深夜福利视频在线观看| 18禁观看日本| 一区二区日韩欧美中文字幕| 久久中文字幕人妻熟女| 99香蕉大伊视频| 亚洲欧洲精品一区二区精品久久久| 亚洲一区二区三区不卡视频| 神马国产精品三级电影在线观看 | 视频区图区小说| 欧美日韩黄片免| 亚洲一区二区三区色噜噜 | 首页视频小说图片口味搜索| 在线国产一区二区在线| 丰满饥渴人妻一区二区三| 国产成人系列免费观看| 久久人妻福利社区极品人妻图片| 日本wwww免费看| 亚洲精品久久午夜乱码| 动漫黄色视频在线观看| 国产精品免费一区二区三区在线| 国产黄色免费在线视频| 男女高潮啪啪啪动态图| 久热爱精品视频在线9| 精品一区二区三卡| 最近最新中文字幕大全免费视频| 性欧美人与动物交配| 欧美日韩国产mv在线观看视频| 亚洲全国av大片| 亚洲av成人av| tocl精华| 黄色怎么调成土黄色| 丰满饥渴人妻一区二区三| 亚洲精品一二三| 久久天躁狠狠躁夜夜2o2o| 成熟少妇高潮喷水视频| 老司机亚洲免费影院| 国产深夜福利视频在线观看| 90打野战视频偷拍视频| 色婷婷av一区二区三区视频| av电影中文网址| 亚洲伊人色综图| 久久久久久久午夜电影 | 国产xxxxx性猛交| 一夜夜www| 久久久国产成人精品二区 | 日韩高清综合在线| 最近最新中文字幕大全电影3 | 午夜成年电影在线免费观看| 国产精品99久久99久久久不卡| 日韩 欧美 亚洲 中文字幕| 狠狠狠狠99中文字幕| 成人国语在线视频| 日韩国内少妇激情av| 亚洲自拍偷在线| 亚洲九九香蕉| www国产在线视频色| 成人永久免费在线观看视频| 亚洲国产精品sss在线观看 | 亚洲av成人一区二区三| av天堂在线播放| 亚洲狠狠婷婷综合久久图片| 国产极品粉嫩免费观看在线| 麻豆一二三区av精品| 999久久久国产精品视频| 国产午夜精品久久久久久| 在线观看日韩欧美| 真人做人爱边吃奶动态| 久久精品人人爽人人爽视色| 黄频高清免费视频| 老司机福利观看| 亚洲av五月六月丁香网| 国产精品久久久人人做人人爽| 老鸭窝网址在线观看| 午夜久久久在线观看| 欧美不卡视频在线免费观看 | 久久人人爽av亚洲精品天堂| 国产一区二区三区视频了| 一二三四社区在线视频社区8| 欧美日韩视频精品一区| 人人妻人人添人人爽欧美一区卜| 一进一出抽搐gif免费好疼 | 在线观看66精品国产| 他把我摸到了高潮在线观看| 黄片大片在线免费观看| 人妻丰满熟妇av一区二区三区| 99久久久亚洲精品蜜臀av| 国产精品九九99| 99久久国产精品久久久| 在线国产一区二区在线| 国产在线观看jvid| 日韩av在线大香蕉| 亚洲国产中文字幕在线视频| 黑人巨大精品欧美一区二区mp4| 91麻豆精品激情在线观看国产 | 婷婷精品国产亚洲av在线| 丝袜美腿诱惑在线| 精品欧美一区二区三区在线| 亚洲欧美精品综合一区二区三区| 在线免费观看的www视频| 制服诱惑二区| 国产成人精品久久二区二区免费| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产欧美一区二区综合| 成人18禁高潮啪啪吃奶动态图| 精品国产超薄肉色丝袜足j| 日本黄色日本黄色录像| 免费在线观看亚洲国产| 一本综合久久免费| 一区福利在线观看| 青草久久国产| 可以在线观看毛片的网站| 婷婷精品国产亚洲av在线| 免费女性裸体啪啪无遮挡网站| 一级片'在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 18禁国产床啪视频网站| 欧美日韩一级在线毛片| 一本综合久久免费| 黄色女人牲交| 麻豆国产av国片精品| 美女福利国产在线| 国产又爽黄色视频| 高清毛片免费观看视频网站 | 久久精品国产清高在天天线| 麻豆国产av国片精品| 亚洲国产欧美网| 香蕉久久夜色| 国产高清激情床上av| a在线观看视频网站| 国产99久久九九免费精品| 日韩欧美三级三区| www.www免费av| xxxhd国产人妻xxx| 男女午夜视频在线观看| 亚洲国产精品一区二区三区在线| 国产人伦9x9x在线观看| 999久久久国产精品视频| 亚洲精品国产精品久久久不卡| 老司机午夜十八禁免费视频| 欧美日韩瑟瑟在线播放| 国产精品免费一区二区三区在线| 无人区码免费观看不卡| 色老头精品视频在线观看| 精品乱码久久久久久99久播| 精品无人区乱码1区二区| 日韩国内少妇激情av| 午夜免费观看网址| 很黄的视频免费| 老司机亚洲免费影院| 亚洲精品一二三| 日本一区二区免费在线视频| 真人做人爱边吃奶动态| 精品国产亚洲在线| 午夜免费观看网址| 超色免费av| 国产日韩一区二区三区精品不卡| 天堂俺去俺来也www色官网| 1024香蕉在线观看| 精品国产乱码久久久久久男人| 在线av久久热| 91在线观看av| 精品人妻1区二区| 免费高清在线观看日韩| 视频在线观看一区二区三区| 国产人伦9x9x在线观看| 一级黄色大片毛片| 精品一区二区三卡| av有码第一页| 日韩人妻精品一区2区三区| 狂野欧美激情性xxxx| 怎么达到女性高潮| 不卡av一区二区三区| 97碰自拍视频| 琪琪午夜伦伦电影理论片6080| 99在线视频只有这里精品首页| 可以免费在线观看a视频的电影网站| 日本 av在线| 亚洲性夜色夜夜综合| 日本wwww免费看| 国产日韩一区二区三区精品不卡| 亚洲欧美激情在线| av欧美777| 大型av网站在线播放| 操出白浆在线播放| 黑人操中国人逼视频| 亚洲伊人色综图| 91成年电影在线观看| 大陆偷拍与自拍| 国产成人一区二区三区免费视频网站| 成人影院久久| 免费在线观看视频国产中文字幕亚洲| 日韩中文字幕欧美一区二区| 国产在线观看jvid| videosex国产| 如日韩欧美国产精品一区二区三区| 欧美大码av| 久久久久久免费高清国产稀缺| 久久天堂一区二区三区四区| 亚洲精品中文字幕一二三四区| 亚洲精品中文字幕在线视频| 国产又爽黄色视频| 美女高潮喷水抽搐中文字幕| 大香蕉久久成人网| 丰满人妻熟妇乱又伦精品不卡| 欧美精品啪啪一区二区三区| 国产精品影院久久| 国产欧美日韩精品亚洲av| 嫁个100分男人电影在线观看| 国产人伦9x9x在线观看| 美国免费a级毛片| 一本大道久久a久久精品| 国产av在哪里看| 一个人观看的视频www高清免费观看 | 美国免费a级毛片| 精品一品国产午夜福利视频| 成人手机av| 久久精品亚洲熟妇少妇任你| 美女高潮喷水抽搐中文字幕| av网站在线播放免费| 每晚都被弄得嗷嗷叫到高潮| 老鸭窝网址在线观看| 看黄色毛片网站| 一区二区日韩欧美中文字幕| 桃色一区二区三区在线观看| 亚洲精品国产区一区二| 精品免费久久久久久久清纯| 久久 成人 亚洲| 国产高清激情床上av| 精品国产一区二区久久| 久久香蕉精品热| 欧美日韩乱码在线| 99精品在免费线老司机午夜| 他把我摸到了高潮在线观看| 老司机午夜福利在线观看视频| 老汉色∧v一级毛片| 51午夜福利影视在线观看| 久久青草综合色| av欧美777| 免费在线观看亚洲国产| 中国美女看黄片| 中文字幕高清在线视频| 久久中文看片网| 国产成人欧美在线观看| 男人舔女人的私密视频| 1024视频免费在线观看| 一本综合久久免费| 五月开心婷婷网| 正在播放国产对白刺激| 搡老岳熟女国产| 国产精品美女特级片免费视频播放器 | 久久久久久免费高清国产稀缺| 亚洲一码二码三码区别大吗| 丝袜美足系列| 亚洲成av片中文字幕在线观看| 国产精品亚洲av一区麻豆| 亚洲专区国产一区二区| 天堂中文最新版在线下载| 久久 成人 亚洲| 亚洲激情在线av| 首页视频小说图片口味搜索| 男人舔女人的私密视频| tocl精华| 国产精品亚洲一级av第二区| 满18在线观看网站| 精品国产国语对白av| 国产精品一区二区三区四区久久 | 麻豆av在线久日| 搡老熟女国产l中国老女人| www.999成人在线观看| 校园春色视频在线观看| 男女午夜视频在线观看| 国产欧美日韩综合在线一区二区| 一区二区三区激情视频| 夜夜爽天天搞| 一级毛片精品| 在线视频色国产色| 18禁裸乳无遮挡免费网站照片 | 桃红色精品国产亚洲av| svipshipincom国产片| 国产精品免费视频内射| 久久中文看片网| 精品熟女少妇八av免费久了| 免费久久久久久久精品成人欧美视频| 国产乱人伦免费视频| 岛国视频午夜一区免费看| 色精品久久人妻99蜜桃| 精品久久久久久久久久免费视频 | 啦啦啦 在线观看视频| 青草久久国产| 美女福利国产在线| 啦啦啦免费观看视频1| 午夜福利在线观看吧| 丁香六月欧美| 欧美亚洲日本最大视频资源| 啦啦啦免费观看视频1| 亚洲中文字幕日韩| 少妇被粗大的猛进出69影院| 一级,二级,三级黄色视频| 夜夜躁狠狠躁天天躁| 中文亚洲av片在线观看爽| 88av欧美| 久久精品亚洲av国产电影网| 欧洲精品卡2卡3卡4卡5卡区| 欧美+亚洲+日韩+国产| 一边摸一边抽搐一进一小说| 淫妇啪啪啪对白视频| 亚洲精品一区av在线观看| 丝袜美腿诱惑在线| 国产单亲对白刺激| 欧美日韩瑟瑟在线播放| 88av欧美| 一进一出抽搐动态| 人人澡人人妻人| 日韩一卡2卡3卡4卡2021年| 免费日韩欧美在线观看| 亚洲激情在线av| 久久久精品国产亚洲av高清涩受| xxx96com| 在线观看一区二区三区激情| 欧美日韩视频精品一区| 怎么达到女性高潮| 黑人巨大精品欧美一区二区mp4| 亚洲男人天堂网一区| 成人黄色视频免费在线看| 看片在线看免费视频| av中文乱码字幕在线| 久久精品国产亚洲av高清一级| 亚洲国产中文字幕在线视频| 天堂俺去俺来也www色官网| 视频区欧美日本亚洲| 亚洲第一欧美日韩一区二区三区| 亚洲午夜精品一区,二区,三区| 亚洲精品久久成人aⅴ小说| 欧美精品一区二区免费开放| 国产成人精品久久二区二区免费| 久久国产精品男人的天堂亚洲| 午夜免费激情av| 欧美日韩国产mv在线观看视频| 在线观看免费视频网站a站| 日韩人妻精品一区2区三区| 亚洲国产欧美网| 久久精品aⅴ一区二区三区四区| 99热国产这里只有精品6| 无限看片的www在线观看| 90打野战视频偷拍视频| 91在线观看av| 国内毛片毛片毛片毛片毛片| 欧美乱码精品一区二区三区| 亚洲九九香蕉| 99国产精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 女人爽到高潮嗷嗷叫在线视频| 露出奶头的视频| 热re99久久国产66热| 黄色视频不卡| 精品高清国产在线一区| 成在线人永久免费视频| 搡老熟女国产l中国老女人| 黑丝袜美女国产一区| 亚洲av成人一区二区三| 国内毛片毛片毛片毛片毛片| 久久精品亚洲熟妇少妇任你| 精品福利永久在线观看| 久久国产精品影院| 一区二区三区精品91| 一边摸一边做爽爽视频免费| 男人的好看免费观看在线视频 | av网站免费在线观看视频| 日韩免费高清中文字幕av| 久久人妻福利社区极品人妻图片| 99国产精品一区二区三区| 51午夜福利影视在线观看| 欧美国产精品va在线观看不卡| 香蕉丝袜av| 好看av亚洲va欧美ⅴa在| 极品教师在线免费播放| 亚洲自偷自拍图片 自拍| 日本免费a在线| 超碰成人久久| cao死你这个sao货| 国产视频一区二区在线看| 桃红色精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 久久这里只有精品19| 日韩精品青青久久久久久| 久久这里只有精品19| 女人高潮潮喷娇喘18禁视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产极品粉嫩免费观看在线| 亚洲国产欧美日韩在线播放| 99久久国产精品久久久| 久久久国产欧美日韩av| 国产精品久久电影中文字幕| 一区二区日韩欧美中文字幕| 欧美久久黑人一区二区| 在线天堂中文资源库| 热re99久久国产66热| 极品教师在线免费播放| 曰老女人黄片| 嫁个100分男人电影在线观看| 国产亚洲av高清不卡| 中文欧美无线码| 黄色 视频免费看| av网站免费在线观看视频| 9191精品国产免费久久| а√天堂www在线а√下载| 国产熟女xx| 十八禁网站免费在线| 国产高清国产精品国产三级| 长腿黑丝高跟| 99久久99久久久精品蜜桃| 91av网站免费观看| 国产乱人伦免费视频| 国产黄色免费在线视频| 成人av一区二区三区在线看| 精品无人区乱码1区二区| 久久人人精品亚洲av| 色综合欧美亚洲国产小说| 午夜福利欧美成人| 国产精品一区二区在线不卡| 久久久水蜜桃国产精品网| 色综合婷婷激情| av在线天堂中文字幕 | 久久久久久久午夜电影 | 麻豆av在线久日| 精品一区二区三区视频在线观看免费 | 制服人妻中文乱码| 亚洲av美国av| 色播在线永久视频| 亚洲精品中文字幕一二三四区| 精品福利观看| 成人黄色视频免费在线看| www.自偷自拍.com| 亚洲欧美精品综合久久99| 久热这里只有精品99| 免费av中文字幕在线| 欧美在线黄色| 国产97色在线日韩免费| 日日摸夜夜添夜夜添小说| 美女国产高潮福利片在线看| 岛国视频午夜一区免费看| 欧美中文日本在线观看视频| 纯流量卡能插随身wifi吗| 十八禁人妻一区二区| 好看av亚洲va欧美ⅴa在| 五月开心婷婷网| 高潮久久久久久久久久久不卡| 丰满的人妻完整版| 亚洲国产中文字幕在线视频| 亚洲成人精品中文字幕电影 | 天堂影院成人在线观看| 精品国产一区二区三区四区第35| 日韩一卡2卡3卡4卡2021年| av中文乱码字幕在线| 91国产中文字幕| 男人舔女人的私密视频| 亚洲人成伊人成综合网2020| av欧美777| 黄片播放在线免费| av在线天堂中文字幕 | 精品国产乱码久久久久久男人| 手机成人av网站| 国产精华一区二区三区| 波多野结衣高清无吗| 国产单亲对白刺激| 欧美老熟妇乱子伦牲交| 精品高清国产在线一区| 精品福利永久在线观看| 99久久久亚洲精品蜜臀av| 国产亚洲欧美98| 国产精品香港三级国产av潘金莲| 亚洲情色 制服丝袜| 国产精品日韩av在线免费观看 | 久久久久国产精品人妻aⅴ院| 麻豆av在线久日| 丝袜在线中文字幕| 欧美丝袜亚洲另类 | 看片在线看免费视频| 少妇被粗大的猛进出69影院| 我的亚洲天堂| 琪琪午夜伦伦电影理论片6080| 午夜91福利影院| 亚洲av熟女| 麻豆成人av在线观看| 久久午夜亚洲精品久久| 午夜免费成人在线视频| 老司机福利观看| 久热这里只有精品99| 欧美久久黑人一区二区| 亚洲av日韩精品久久久久久密| 亚洲欧美一区二区三区久久| 高清黄色对白视频在线免费看| 亚洲精品成人av观看孕妇| 制服诱惑二区| 午夜免费成人在线视频| 多毛熟女@视频| 欧美黄色片欧美黄色片| 真人做人爱边吃奶动态| 51午夜福利影视在线观看| 国产精品1区2区在线观看.| 成人黄色视频免费在线看| 黑人巨大精品欧美一区二区蜜桃| 99久久综合精品五月天人人| 亚洲激情在线av| 亚洲色图综合在线观看| 国产av精品麻豆| 亚洲一区中文字幕在线| 久久精品成人免费网站| 三上悠亚av全集在线观看| 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看| 久久青草综合色| 免费观看人在逋| 老司机在亚洲福利影院| 如日韩欧美国产精品一区二区三区| 国产深夜福利视频在线观看| 精品欧美一区二区三区在线| 欧美老熟妇乱子伦牲交| 久久精品国产综合久久久| 亚洲av熟女| 一区二区三区激情视频| 国产国语露脸激情在线看| 一级a爱片免费观看的视频| 成人18禁在线播放| 啦啦啦 在线观看视频| 一本综合久久免费| 在线观看午夜福利视频| 久久久久久人人人人人| 欧美丝袜亚洲另类 | 757午夜福利合集在线观看| 久久精品91无色码中文字幕| 成人国产一区最新在线观看| 人妻久久中文字幕网| 亚洲第一青青草原| 亚洲午夜理论影院| 电影成人av| 欧美日韩亚洲高清精品| 在线av久久热| 99热只有精品国产| 咕卡用的链子| 久久天躁狠狠躁夜夜2o2o| 欧美日韩亚洲高清精品| 久久久久国产一级毛片高清牌| 欧美午夜高清在线| 窝窝影院91人妻| 日本免费a在线| 欧美成人午夜精品| 99热国产这里只有精品6| 国产又爽黄色视频| 国产成人av激情在线播放| 国产一区二区三区在线臀色熟女 | 午夜免费激情av| 免费看a级黄色片| 夫妻午夜视频| 亚洲精品国产精品久久久不卡| 国产精品乱码一区二三区的特点 | 最新在线观看一区二区三区| 欧美老熟妇乱子伦牲交| 亚洲欧美一区二区三区久久| 中亚洲国语对白在线视频| 国产成人欧美| 嫩草影院精品99| 水蜜桃什么品种好| 一级作爱视频免费观看| 免费少妇av软件| 国产日韩一区二区三区精品不卡| 亚洲九九香蕉| 黄网站色视频无遮挡免费观看| 三级毛片av免费| 久久狼人影院| 窝窝影院91人妻| 亚洲一码二码三码区别大吗| 国产成人欧美| 少妇的丰满在线观看| 久久精品国产亚洲av香蕉五月| 一级毛片女人18水好多| 亚洲av第一区精品v没综合| 日本黄色视频三级网站网址| 女人精品久久久久毛片| 宅男免费午夜| 高清欧美精品videossex| 真人一进一出gif抽搐免费| 色在线成人网| 国产高清videossex| 亚洲av电影在线进入| 国产精品1区2区在线观看.| 精品久久久久久电影网| 日韩欧美一区视频在线观看| 我的亚洲天堂| 日韩欧美免费精品| 热re99久久精品国产66热6| 久久草成人影院| 深夜精品福利| 欧美 亚洲 国产 日韩一| 久久精品亚洲av国产电影网| 在线观看一区二区三区激情| 9色porny在线观看| 日本免费a在线| 国产一区二区三区视频了| 中文字幕色久视频| 欧美黑人欧美精品刺激| 欧美国产精品va在线观看不卡| 国产欧美日韩一区二区精品| 国产一区在线观看成人免费| 看免费av毛片| 夜夜爽天天搞| 天天躁夜夜躁狠狠躁躁| 精品一区二区三区四区五区乱码| 色在线成人网| 91麻豆av在线| 欧美精品啪啪一区二区三区| 高清欧美精品videossex| 交换朋友夫妻互换小说|