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

    基于復(fù)雜網(wǎng)絡(luò)和灰色關(guān)聯(lián)分析的引水隧洞坍塌事故致因分析

    2024-10-26 00:00:00張勇張琪劉少乾王沛豐
    項(xiàng)目管理技術(shù) 2024年10期

    摘要:引水隧洞建設(shè)中坍塌事故多發(fā)且后果惡劣,為有效預(yù)防坍塌事故,提出基于復(fù)雜網(wǎng)絡(luò)和灰色關(guān)聯(lián)分析的引水隧洞坍塌事故關(guān)鍵致因識(shí)別方法。首先,深入分析歷史事故,從4個(gè)維度系統(tǒng)地總結(jié)風(fēng)險(xiǎn)因素集;其次,依據(jù)風(fēng)險(xiǎn)因素及其影響關(guān)系構(gòu)建坍塌事故致因網(wǎng)絡(luò)模型,并多指標(biāo)、多維度地對(duì)節(jié)點(diǎn)進(jìn)行重要度分析;最后,鑒于單一特性可能無(wú)法全面呈現(xiàn)節(jié)點(diǎn)在網(wǎng)絡(luò)中的影響力,引入灰色關(guān)聯(lián)分析可有效融合節(jié)點(diǎn)的多個(gè)特征指標(biāo),進(jìn)而確定各風(fēng)險(xiǎn)因素的重要程度,并從工程實(shí)踐角度提出建議。

    關(guān)鍵詞:引水隧洞;坍塌事故;灰色關(guān)聯(lián)分析;復(fù)雜網(wǎng)絡(luò);致因分析;風(fēng)險(xiǎn)管理

    0 引言

    20世紀(jì)以來(lái),我國(guó)對(duì)水利工程建設(shè)的資金投入不斷增加,水利工程事業(yè)快速發(fā)展,先后興建了滇中引水、引漢濟(jì)渭、引江濟(jì)淮等工程,實(shí)現(xiàn)了水資源的調(diào)配,促進(jìn)了缺水地區(qū)的經(jīng)濟(jì)發(fā)展[1]。引水隧洞

    作為引水工程的關(guān)鍵設(shè)施,具有建設(shè)周期長(zhǎng)、建設(shè)難度大、施工風(fēng)險(xiǎn)高等特點(diǎn),對(duì)施工效率和施工安全構(gòu)成了巨大挑戰(zhàn)。由于地質(zhì)條件、施工技術(shù)等多重限制,引水隧洞施工建設(shè)中事故頻發(fā),不僅加大了施工難度和建設(shè)成本,還可能引發(fā)次生危害,威脅社會(huì)穩(wěn)定,其中坍塌事故占比最大且人員傷亡最嚴(yán)重[2]。因此,如何有效預(yù)防引水隧洞坍塌事故,保障施工安全,是一項(xiàng)重要課題。

    以往關(guān)于引水隧洞坍塌事故的研究聚焦于兩個(gè)方面。一是基于實(shí)際數(shù)據(jù)和特征參數(shù)分析的預(yù)測(cè)預(yù)警研究。例如,Chen等[3]以引松引水工程的18組坍塌數(shù)據(jù)為基礎(chǔ),引入特征參數(shù)TPI(巖石扭剪掘進(jìn)指標(biāo))對(duì)塌方段進(jìn)行了準(zhǔn)確預(yù)測(cè);劉詩(shī)洋等[4]將FPI(現(xiàn)場(chǎng)貫入指標(biāo))和TPI(巖石扭剪掘進(jìn)指標(biāo))作為機(jī)器學(xué)習(xí)的訓(xùn)練對(duì)象,結(jié)合時(shí)間序列預(yù)測(cè)法分析了工程塌方概率;Hou等[5]以松花江輸水工程TBM運(yùn)行數(shù)據(jù)和地質(zhì)條件數(shù)據(jù)為基礎(chǔ),基于LSTM算法構(gòu)建了隧道坍塌預(yù)警指標(biāo)體系;裴成元等[2]基于掘進(jìn)數(shù)據(jù)和工程塌方數(shù)據(jù)計(jì)算出圍巖特征參數(shù),并通過(guò)理論、實(shí)驗(yàn)等方式驗(yàn)證,以此作為塌方的輔助判斷依據(jù)。二是集中于引水隧洞坍塌風(fēng)險(xiǎn)評(píng)估研究,學(xué)者們多采用模糊數(shù)學(xué)與專(zhuān)家系統(tǒng)相結(jié)合的方法。例如,Shin等[6]基于相互作用矩陣和神經(jīng)網(wǎng)絡(luò),提出了KICT隧道坍塌危險(xiǎn)指數(shù)(KTH-Index)用于評(píng)估隧道坍塌危險(xiǎn)等級(jí),并將評(píng)估程序系統(tǒng)化,為隧道坍塌風(fēng)險(xiǎn)研究奠定了基礎(chǔ)。后來(lái),學(xué)者們致力于運(yùn)用組合模型以提升風(fēng)險(xiǎn)評(píng)估的精度。例如,譚立偉等[7]選用改進(jìn)G2法和熵權(quán)法對(duì)風(fēng)險(xiǎn)因素賦權(quán),構(gòu)建G2-EWM-TOPSIS模型對(duì)滇中工程引水隧洞塌方風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià),并得出相應(yīng)的風(fēng)險(xiǎn)等級(jí);李燕等[8]考慮了風(fēng)險(xiǎn)的不確定性及評(píng)價(jià)的模糊性,建立了基于熵理論-可拓云模型的引水隧洞塌方風(fēng)險(xiǎn)評(píng)價(jià)模型,并驗(yàn)證了該方法的有效性;呂擎峰等[9]考慮到社會(huì)、經(jīng)濟(jì)等影響因素,提出一種基于FAHP法和后果當(dāng)量估計(jì)法的隧道塌方風(fēng)險(xiǎn)評(píng)估模型;Li等[10]提出了一種基于改進(jìn)量化理論III和客觀熵權(quán)重法-層次分析法(EW-AHP)耦合權(quán)重的隧道坍塌風(fēng)險(xiǎn)評(píng)估模型,并通過(guò)實(shí)際工程驗(yàn)證了該方法的可靠性;Ou等[11]基于D-S證據(jù)理論和超前地質(zhì)預(yù)報(bào),建立了隧道坍塌綜合風(fēng)險(xiǎn)評(píng)估模型,豐富和完善了隧道坍塌風(fēng)險(xiǎn)評(píng)估方法;Chen等[12]提出了一種基于案例推理、粗糙集理論和未確知測(cè)度集對(duì)分析理論(UM-SPA)的隧道塌方風(fēng)險(xiǎn)綜合評(píng)價(jià)方法,提高了風(fēng)險(xiǎn)評(píng)估的可靠性。

    以往的引水隧洞坍塌預(yù)警預(yù)測(cè)和風(fēng)險(xiǎn)評(píng)估為預(yù)防坍塌事故提供了有力指導(dǎo)。但在評(píng)估風(fēng)險(xiǎn)致險(xiǎn)能力過(guò)程中,以往的研究未能考慮風(fēng)險(xiǎn)傳遞對(duì)事故的催化作用。同時(shí),以往研究更重視地質(zhì)因素和施工因素對(duì)隧洞坍塌事故的“貢獻(xiàn)”,對(duì)人為因素和管理因素的考量較為欠缺。因此,本研究以實(shí)際案例為研究基礎(chǔ),首先,結(jié)合相關(guān)文獻(xiàn),從工程地質(zhì)風(fēng)險(xiǎn)、自然風(fēng)險(xiǎn)、設(shè)計(jì)施工風(fēng)險(xiǎn)、組織管理風(fēng)險(xiǎn)4個(gè)方面構(gòu)建引水隧洞坍塌事故風(fēng)險(xiǎn)因素集;其次,構(gòu)建引水隧洞坍塌事故致因網(wǎng)絡(luò)模型并分析節(jié)點(diǎn)的特征指標(biāo);最后,為消除獨(dú)立的指標(biāo)分析對(duì)于節(jié)點(diǎn)重要程度評(píng)價(jià)的偏差,本研究利用灰色關(guān)聯(lián)分析定量計(jì)算各個(gè)指標(biāo)的單因素關(guān)聯(lián)度,融合節(jié)點(diǎn)的多個(gè)評(píng)價(jià)指標(biāo),更加全面地衡量節(jié)點(diǎn)重要性。

    1 事故風(fēng)險(xiǎn)因素識(shí)別

    1.1 事故統(tǒng)計(jì)及信息提取

    為識(shí)別影響引水隧洞坍塌的風(fēng)險(xiǎn)因素,對(duì)引水隧洞坍塌事故案例進(jìn)行梳理。通過(guò)各地區(qū)、各級(jí)應(yīng)急管理廳網(wǎng)站,結(jié)合相關(guān)文獻(xiàn)、公開(kāi)報(bào)道等信息渠道,搜集整理出近年來(lái)發(fā)生的引水隧洞坍塌事故共41起。為保證研究的可靠性和準(zhǔn)確性,篩選并排除數(shù)據(jù)資料不完整或記錄不詳實(shí)的案例,最終得到32份典型的引水隧洞坍塌事故報(bào)告,部分引水隧洞坍塌事故報(bào)告見(jiàn)表1。

    1.2 事故致因風(fēng)險(xiǎn)因素集

    以事故報(bào)告為基礎(chǔ),結(jié)合相關(guān)文獻(xiàn)和引水隧洞施工難點(diǎn),初步整理致因因素。參考《水工隧洞設(shè)計(jì)規(guī)范》(SL279—2016)和《水工隧洞安全監(jiān)測(cè)技術(shù)規(guī)范》(SL764—2018)等規(guī)范標(biāo)準(zhǔn),以系統(tǒng)思維為原則,進(jìn)一步篩選和修正風(fēng)險(xiǎn)因素。最終,從工程地質(zhì)風(fēng)險(xiǎn)、自然風(fēng)險(xiǎn)、設(shè)計(jì)施工風(fēng)險(xiǎn)、組織管理風(fēng)險(xiǎn)4個(gè)方面總結(jié)出引水隧洞坍塌事故風(fēng)險(xiǎn)因素集,共包含27個(gè)二級(jí)風(fēng)險(xiǎn)因素,引水隧洞坍塌事故風(fēng)險(xiǎn)因素集見(jiàn)表2。

    2 基于復(fù)雜網(wǎng)絡(luò)和灰色關(guān)聯(lián)分析的事故致因分析方法

    2.1 復(fù)雜網(wǎng)絡(luò)

    復(fù)雜網(wǎng)絡(luò)將具有高復(fù)雜性的實(shí)際問(wèn)題抽象為網(wǎng)絡(luò)圖的形式,能夠?qū)崿F(xiàn)系統(tǒng)內(nèi)部要素間的關(guān)系、幾何性質(zhì)和演化規(guī)律的可視化。復(fù)雜網(wǎng)絡(luò)中,節(jié)點(diǎn)的各項(xiàng)特征指標(biāo)從不同角度體現(xiàn)了節(jié)點(diǎn)在網(wǎng)絡(luò)中的影響程度,常見(jiàn)的特征指標(biāo)主要包括以下4點(diǎn):

    (1)節(jié)點(diǎn)的度。

    有向網(wǎng)絡(luò)中,節(jié)點(diǎn)的度分為入度值和出度值,分別指目標(biāo)節(jié)點(diǎn)連接到其他節(jié)點(diǎn)的邊的數(shù)目和其他節(jié)點(diǎn)指向目標(biāo)節(jié)點(diǎn)的邊的數(shù)目??偠戎礙i公式如下

    式中,aij為節(jié)點(diǎn)i指向其他節(jié)點(diǎn)j的邊數(shù);aji為其他節(jié)點(diǎn)j指向節(jié)點(diǎn)i的邊數(shù);n為節(jié)點(diǎn)個(gè)數(shù)??偠戎翟酱螅瑒t與該節(jié)點(diǎn)相連的節(jié)點(diǎn)越多,在局部范圍內(nèi)對(duì)周邊節(jié)點(diǎn)的影響程度越高。

    (2)接近中心度。

    接近中心度被用來(lái)描述網(wǎng)絡(luò)中節(jié)點(diǎn)分布的緊湊程度,指目標(biāo)節(jié)點(diǎn)與網(wǎng)絡(luò)中其他節(jié)點(diǎn)的平均最短路徑。接近中心度CC(i)-1公式如下

    式中,dij為節(jié)點(diǎn)i到節(jié)點(diǎn)j的網(wǎng)絡(luò)距離;n為節(jié)點(diǎn)個(gè)數(shù)。接近中心度越大,則該節(jié)點(diǎn)到達(dá)其他節(jié)點(diǎn)的距離越短,信息在網(wǎng)絡(luò)中的傳播速率越快。

    (3)介數(shù)中心度。

    介數(shù)中心度反映單個(gè)節(jié)點(diǎn)作為媒介傳遞風(fēng)險(xiǎn)的能力,指網(wǎng)絡(luò)中經(jīng)過(guò)目標(biāo)節(jié)點(diǎn)的最短路徑數(shù)目與所有最短路徑數(shù)目的比值。介數(shù)中心度CB(e)公式如下

    式中,gij為節(jié)點(diǎn)i到節(jié)點(diǎn)j的最短路徑數(shù)量;gij(e)為網(wǎng)絡(luò)中g(shù)ij中經(jīng)過(guò)節(jié)點(diǎn)e的最短路徑數(shù)量;n為節(jié)點(diǎn)個(gè)數(shù)。介數(shù)中心度越大,說(shuō)明節(jié)點(diǎn)在整個(gè)網(wǎng)絡(luò)中的作用和影響力越大。

    (4)特征向量中心度。

    特征向量中心度認(rèn)為節(jié)點(diǎn)的重要性與其鄰接節(jié)點(diǎn)的數(shù)量及其影響力有關(guān),特征向量中心度ECi公式如下

    式中,aij為節(jié)點(diǎn)i指向其他節(jié)點(diǎn)的邊數(shù);xj為一個(gè)相關(guān)特征向量中心度;n為節(jié)點(diǎn)個(gè)數(shù);λ為常數(shù)。

    2.2 灰色關(guān)聯(lián)分析

    灰色關(guān)聯(lián)分析能夠量化系統(tǒng)中各個(gè)因素的相關(guān)性,其計(jì)算步驟如下:

    (1)確定特征指標(biāo)矩陣。假設(shè)網(wǎng)絡(luò)中有m個(gè)節(jié)點(diǎn),每個(gè)節(jié)點(diǎn)有n個(gè)特征指標(biāo)。aij表示第i個(gè)節(jié)點(diǎn)的第j個(gè)特征指標(biāo)值,其中i=1,2,…,m;j=1,2,…,n。則網(wǎng)絡(luò)的特征指標(biāo)矩陣A為

    (2)歸一化處理。由于各指標(biāo)對(duì)節(jié)點(diǎn)的評(píng)價(jià)維度不同,為了減少數(shù)據(jù)波動(dòng),研究采用均值法對(duì)初始數(shù)據(jù)進(jìn)行處理,公式如下

    式中,xij為處理后的值。歸一化處理后,得到矩陣X為

    (3)確定比較序列和參考序列。矩陣X中的每一列都為比較序列,選取各指標(biāo)標(biāo)準(zhǔn)化后的最大值作為參考序列C=(C1,C2,…,Cn),Ci=maxxij。

    M2ttQyHOr16E+gwlRdm+sQ==

    (4)計(jì)算灰色關(guān)聯(lián)系數(shù)ζj(k)。公式如下

    式中,ζj(k)為矩陣X中xij與參考序列C的灰色關(guān)聯(lián)系數(shù);Cj-xkj為差序列,即各數(shù)據(jù)與參考序列的差距;β為分辨系數(shù),通常取β=0.5,k=1,2,…,m。

    (5)計(jì)算灰色關(guān)聯(lián)度R(C,xj)。公式如下

    關(guān)聯(lián)度R(C,xi)越大,說(shuō)明比較序列xi與參考數(shù)列C越相似,即節(jié)點(diǎn)i越接近理想對(duì)象,其重要程度也越高。

    (6)計(jì)算各指標(biāo)的權(quán)重rj。公式如下

    式中,Rn為節(jié)點(diǎn)第n項(xiàng)指標(biāo)的關(guān)聯(lián)度。

    (7)根據(jù)各節(jié)點(diǎn)的特征指標(biāo)初始值和指標(biāo)權(quán)重,計(jì)算節(jié)點(diǎn)的重要度Zk,公式如下

    式中,Zk為第k個(gè)節(jié)點(diǎn)的重要度;rj為節(jié)點(diǎn)第j個(gè)特征指標(biāo)的權(quán)重;xkj為第k個(gè)節(jié)點(diǎn)的第j項(xiàng)指標(biāo)標(biāo)準(zhǔn)化后的樣本值。

    2.3 基于復(fù)雜網(wǎng)絡(luò)和灰色關(guān)聯(lián)分析的事故致因分析方法

    早期研究[13-14]通過(guò)復(fù)雜網(wǎng)絡(luò)模型進(jìn)行節(jié)點(diǎn)重要性評(píng)價(jià)時(shí)僅基于單一屬性,評(píng)價(jià)結(jié)果并不準(zhǔn)確。后期[15-16]考慮了節(jié)點(diǎn)的多項(xiàng)特征指標(biāo),再通過(guò)各項(xiàng)指標(biāo)下的節(jié)點(diǎn)排序、分析和總結(jié)來(lái)識(shí)別關(guān)鍵因素,但未能定量整合節(jié)點(diǎn)各方面的表現(xiàn)力。復(fù)雜網(wǎng)絡(luò)中節(jié)點(diǎn)的各項(xiàng)指標(biāo)對(duì)于評(píng)價(jià)節(jié)點(diǎn)重要性各有側(cè)重,同時(shí)也各有欠缺。

    節(jié)點(diǎn)的度值作為最基礎(chǔ)的指標(biāo),直接反映節(jié)點(diǎn)的連接數(shù)量,但忽略了連接“質(zhì)量”;接近中心度通過(guò)節(jié)點(diǎn)間的最短距離反映節(jié)點(diǎn)傳遞信息的快慢,卻忽略了網(wǎng)絡(luò)局部結(jié)構(gòu);特征向量中心度考慮了節(jié)點(diǎn)的連接“質(zhì)量”,從全局視角評(píng)估節(jié)點(diǎn)的重要性;介數(shù)中心度能夠體現(xiàn)節(jié)點(diǎn)在網(wǎng)絡(luò)中的聯(lián)結(jié)能力和控制能力,卻忽略了節(jié)點(diǎn)局部影響力。單獨(dú)分析某項(xiàng)指標(biāo)并不能全面體現(xiàn)節(jié)點(diǎn)的影響力,鑒于此,本研究提出一種基于復(fù)雜網(wǎng)絡(luò)和灰色關(guān)聯(lián)分析的綜合方法,用于引水隧洞塌方事故致因研究。在復(fù)雜網(wǎng)絡(luò)事故致因模型分析的基礎(chǔ)上,綜合各指標(biāo)的考量角度,通過(guò)灰色關(guān)聯(lián)分析融合節(jié)點(diǎn)的度值、接近中心度、介數(shù)中心度和特征向量中心度,計(jì)算出系統(tǒng)內(nèi)各節(jié)點(diǎn)的關(guān)聯(lián)度,實(shí)現(xiàn)多角度的全面評(píng)價(jià)節(jié)點(diǎn)重要性。

    3 基于復(fù)雜網(wǎng)絡(luò)和灰色關(guān)聯(lián)分析的事故致因分析

    3.1 致因網(wǎng)絡(luò)構(gòu)建

    在引水隧洞施工安全系統(tǒng)中,將風(fēng)險(xiǎn)因素抽象為網(wǎng)絡(luò)節(jié)點(diǎn),將要素間的依賴(lài)關(guān)系抽象為邊,基于復(fù)雜網(wǎng)絡(luò)構(gòu)建引水隧洞坍塌事故致因網(wǎng)絡(luò)。通過(guò)Gephi軟件將網(wǎng)絡(luò)可視化,最終構(gòu)建的致因模型包括27個(gè)節(jié)點(diǎn),96條連邊,引水隧洞坍塌事故致因網(wǎng)絡(luò)如圖1所示。

    3.2 網(wǎng)絡(luò)節(jié)點(diǎn)拓?fù)涮匦苑治?/p>

    描述網(wǎng)絡(luò)結(jié)構(gòu)和節(jié)點(diǎn)特征的拓?fù)鋮?shù)是定量評(píng)估復(fù)雜網(wǎng)絡(luò)模型的主要手段。本研究選取節(jié)點(diǎn)的度、接近中心度、中介中心度和特征向量中心度,用來(lái)描述節(jié)點(diǎn)在系統(tǒng)中的特征。

    3.2.1 節(jié)點(diǎn)的度

    節(jié)點(diǎn)的度分布如圖2所示。C8是出度最大的節(jié)點(diǎn),表明該因素直接誘發(fā)其他風(fēng)險(xiǎn)的能力最強(qiáng)。出度值較高的節(jié)點(diǎn)C8、D2、D4、D1、D7大多屬于組織管理風(fēng)險(xiǎn),表明它們?cè)陲L(fēng)險(xiǎn)鏈中處于上游,可以看作引水隧洞坍塌事故的根源。A3是入度最大的節(jié)點(diǎn),說(shuō)明圍巖變形失穩(wěn)在系統(tǒng)中最易受到其他因素的影響。同時(shí),A3也是總度值最大的節(jié)點(diǎn),說(shuō)明該節(jié)點(diǎn)在網(wǎng)絡(luò)中的重要性較突顯,易與其他因素互相影響。C1、D4、C8、A2、D2的總度值也較大,與其他風(fēng)險(xiǎn)因素聯(lián)系緊密。針對(duì)性管控這些度值較大的節(jié)點(diǎn)能夠切斷風(fēng)險(xiǎn)傳播路徑,有效降低風(fēng)險(xiǎn)鏈的連通性。

    3.2.2 接近中心度

    節(jié)點(diǎn)的接近中心度分布如圖3所示。A3具有最高的接近中心度,表示該節(jié)點(diǎn)與其他節(jié)點(diǎn)的平均距離較小。由于更接近網(wǎng)絡(luò)的中心位置,節(jié)點(diǎn)A3在網(wǎng)絡(luò)中能夠快速將風(fēng)險(xiǎn)擴(kuò)散至其他節(jié)點(diǎn)。C1、D4、C8、D7節(jié)點(diǎn)也表現(xiàn)出較大的接近中心度,傳遞風(fēng)險(xiǎn)至其他節(jié)點(diǎn)的速度也相對(duì)較快。高接近中心度的節(jié)點(diǎn)是風(fēng)險(xiǎn)傳播的加速器,因此,管控此類(lèi)節(jié)點(diǎn)能夠有效降低風(fēng)險(xiǎn)的傳播效率保障風(fēng)險(xiǎn)網(wǎng)絡(luò)的穩(wěn)定性。

    3.2.3 介數(shù)中心度

    介數(shù)中心度高的節(jié)點(diǎn)在網(wǎng)絡(luò)中也具有關(guān)鍵風(fēng)險(xiǎn)傳遞作用,節(jié)點(diǎn)的介數(shù)中心度分布如圖4所示。該致因網(wǎng)絡(luò)中存在9個(gè)介數(shù)中心度為0的節(jié)點(diǎn),如A1、A3、B3、C8、D2等,這些節(jié)點(diǎn)未發(fā)揮媒介效應(yīng)。介數(shù)中心度最大的節(jié)點(diǎn)為C1,說(shuō)明支撐支護(hù)工程不合理在該風(fēng)險(xiǎn)網(wǎng)絡(luò)中的調(diào)節(jié)能力和控制能力較強(qiáng)。該節(jié)點(diǎn)的微小變化極大概率引起其他節(jié)點(diǎn)的狀態(tài)改變,進(jìn)而引起整體風(fēng)險(xiǎn)波動(dòng)。實(shí)際施工過(guò)程中,若配置了無(wú)效或不及時(shí)的支撐支護(hù)工程,后續(xù)作業(yè)中極大可能會(huì)引發(fā)圍巖結(jié)構(gòu)失穩(wěn)、應(yīng)力集中等危險(xiǎn)情況,進(jìn)而導(dǎo)致局部甚至整體坍塌。B1、A4、D4、D6等節(jié)點(diǎn)的介數(shù)中心度也較大,對(duì)這些節(jié)點(diǎn)進(jìn)行控制,可以降低風(fēng)險(xiǎn)點(diǎn)之間的傳播效率,減少網(wǎng)絡(luò)中的連鎖反應(yīng)。

    3.2.4 特征向量中心度

    特征向量中心度認(rèn)為,如果網(wǎng)絡(luò)中有許多重要節(jié)點(diǎn)與某節(jié)點(diǎn)相連,那么該節(jié)點(diǎn)在網(wǎng)絡(luò)中的重要程度也對(duì)應(yīng)較高,節(jié)點(diǎn)的特征向量中心度分布如圖5所示。特征向量中心度最高的風(fēng)險(xiǎn)因素為A3,為0.446 59,表明圍巖變形失穩(wěn)與多個(gè)關(guān)鍵節(jié)點(diǎn)相連,其穩(wěn)定性對(duì)整個(gè)工程安全至關(guān)重要。C1、A2、C5、D7等節(jié)點(diǎn)的特征向量中心度也較高,表明支護(hù)工程的設(shè)計(jì)與實(shí)施、結(jié)構(gòu)面黏聚力強(qiáng)度、施工方法、安全防護(hù)措施等在網(wǎng)絡(luò)中的影響力也較大。

    3.3 致因因素重要度分析

    由上文可知,節(jié)點(diǎn)A3和C8的總度值和接近中心度較高,通常意味著它們?cè)诰W(wǎng)絡(luò)中具有較高的連接度和接近性,理論上應(yīng)該在風(fēng)險(xiǎn)傳遞中起著重要作用。然而,它們的出度或入度為0,即它們既不向其他節(jié)點(diǎn)發(fā)送信息,也不從其他節(jié)點(diǎn)接收信息,即它們?cè)陲L(fēng)險(xiǎn)傳遞過(guò)程中并未發(fā)揮作用,所以A3和C8的中介中心度為0。這一現(xiàn)象表明,通過(guò)單一指標(biāo)衡量節(jié)點(diǎn)的重要性有失偏頗,可能造成關(guān)鍵節(jié)點(diǎn)判別失誤,影響研究結(jié)果的合理性和科學(xué)性。因此,本研究以節(jié)點(diǎn)的多項(xiàng)特征指標(biāo)為基礎(chǔ),進(jìn)行灰色關(guān)聯(lián)度分析。根據(jù)初始的特征指標(biāo)矩陣,通過(guò)式(5)對(duì)所有值進(jìn)行歸一化處理,確定參考序列和比較序列后,根據(jù)式(6~7)計(jì)算得出指標(biāo)的灰色關(guān)聯(lián)系數(shù)和灰色關(guān)聯(lián)度。根據(jù)式(8~9)得出指標(biāo)權(quán)重和節(jié)點(diǎn)的綜合重要度,引水隧洞坍塌事故致因因素重要度排序見(jiàn)表3。

    3.4 結(jié)果分析

    (1)結(jié)果表明,節(jié)點(diǎn)C1的綜合重要度最大,即支撐支護(hù)工程不合理是對(duì)引水隧洞坍塌事故影響程度最深的因素,這說(shuō)明強(qiáng)支護(hù)工程是確保施工穩(wěn)定性的關(guān)鍵措施。實(shí)際工程中,支護(hù)工程強(qiáng)度不夠、支撐不及時(shí)等情況較為常見(jiàn)。因此,施工前應(yīng)深入地質(zhì)勘探,依據(jù)地層條件、巖土性質(zhì)等合理設(shè)計(jì)和選擇支護(hù)類(lèi)型、支護(hù)結(jié)構(gòu)和材料;施工中應(yīng)嚴(yán)格遵守設(shè)計(jì)方案,同時(shí)高度監(jiān)測(cè)支護(hù)結(jié)構(gòu)的應(yīng)力、位移等重要參數(shù),出現(xiàn)異常監(jiān)測(cè)結(jié)果及時(shí)調(diào)整支護(hù)設(shè)計(jì),保證隧洞的承載能力和結(jié)構(gòu)完整性。

    (2)D4、A3、B1、D7的影響程度也較高,這些節(jié)點(diǎn)也是引水隧洞坍塌防控工作需重點(diǎn)考慮的因素。因此,應(yīng)定期進(jìn)行風(fēng)險(xiǎn)管理培訓(xùn)和工程風(fēng)險(xiǎn)評(píng)估,并確保相關(guān)人員及時(shí)了解風(fēng)險(xiǎn)狀況;加強(qiáng)圍巖變形監(jiān)測(cè)工作,及時(shí)干預(yù)異常情況,動(dòng)態(tài)調(diào)整施工方案以適應(yīng)環(huán)境變化;設(shè)計(jì)階段考慮防水措施,加強(qiáng)對(duì)地下水情的監(jiān)測(cè)并評(píng)估其影響,科學(xué)制訂隧洞排水計(jì)劃;明確安全責(zé)任,加強(qiáng)安全監(jiān)督和安全檢查,建立健全安全管理體系。

    (3)整體上看,重要度較高的節(jié)點(diǎn)中,工程地質(zhì)風(fēng)險(xiǎn)和組織管理風(fēng)險(xiǎn)兩類(lèi)風(fēng)險(xiǎn)占比達(dá)到72.73%,包括D4、A3、D7、A2、A7、A4、D6、D2等,說(shuō)明這兩類(lèi)風(fēng)險(xiǎn)是造成引水隧洞坍塌事故發(fā)生的主要原因。在勘探設(shè)計(jì)階段,應(yīng)深入調(diào)查工程地質(zhì)情況并進(jìn)行全面的風(fēng)險(xiǎn)評(píng)估,合理地選擇和設(shè)計(jì)隧洞結(jié)構(gòu)、施工技術(shù)、施工方法等。施工階段加強(qiáng)監(jiān)測(cè)、及時(shí)干預(yù)并嚴(yán)格執(zhí)行應(yīng)急預(yù)案,能極大降低工程地質(zhì)風(fēng)險(xiǎn)對(duì)塌方事故的威脅。另外,應(yīng)優(yōu)化組織的安全管理流程,同時(shí)加強(qiáng)安全教育管理與培訓(xùn),強(qiáng)化團(tuán)隊(duì)的專(zhuān)業(yè)能力,提高其風(fēng)險(xiǎn)識(shí)別和風(fēng)險(xiǎn)應(yīng)對(duì)能力。

    4 結(jié)語(yǔ)

    (1)以官方出具的事故報(bào)告為基礎(chǔ),構(gòu)建了由4個(gè)一級(jí)風(fēng)險(xiǎn)因素和27個(gè)二級(jí)風(fēng)險(xiǎn)因素組成的引水隧洞坍塌事故致因系統(tǒng),相比于其他文獻(xiàn)總結(jié)等方式,基于歷史案例的風(fēng)險(xiǎn)因素分析更具有現(xiàn)實(shí)性。此外,在充分考慮風(fēng)險(xiǎn)因素間的依賴(lài)關(guān)系的基礎(chǔ)上構(gòu)建事故致因網(wǎng)絡(luò)模型,并從多角度衡量各風(fēng)險(xiǎn)因素在系統(tǒng)中的影響力。

    (2)針對(duì)單一指標(biāo)衡量節(jié)點(diǎn)重要程度的片面性和主觀性,本研究引入灰色關(guān)聯(lián)分析有效融合節(jié)點(diǎn)的多項(xiàng)指標(biāo),從節(jié)點(diǎn)的連接次數(shù)、連接質(zhì)量、路徑距離等多角度綜合考量節(jié)點(diǎn)重要性,更準(zhǔn)確地識(shí)別引水隧洞坍塌事故關(guān)鍵致因,為施工安全管理提供理論支持和實(shí)踐指導(dǎo)。

    (3)由于事故部分資料的缺乏,本研究未能考慮時(shí)間效應(yīng)對(duì)風(fēng)險(xiǎn)狀態(tài)的影響。后續(xù)研究中可基于不同時(shí)間節(jié)點(diǎn),探究風(fēng)險(xiǎn)狀態(tài)及風(fēng)險(xiǎn)發(fā)生概率對(duì)于各風(fēng)險(xiǎn)因素致險(xiǎn)能力的影響。另外,本研究在分析事故報(bào)告時(shí)仍存在主觀性,后續(xù)研究將考慮利用文本挖掘技術(shù)使得致因因素體系構(gòu)建更客觀。

    參考文獻(xiàn)

    [1]貢力,路瑞琴,靳春玲,等.基于博弈-改進(jìn)可拓理論的寒冷地區(qū)長(zhǎng)距離明渠冬季運(yùn)行安全評(píng)價(jià)[J].自然災(zāi)害學(xué)報(bào),2019,28(6):81-92.

    [2]裴成元,張?jiān)旗?,劉軍生,?基于TBM掘進(jìn)大數(shù)據(jù)和特征參數(shù)的引水隧洞塌方分析[J].隧道建設(shè)(中英文),2024,44(5):952-963.

    [3]CHEN Z, ZHANG Y, LI J, et al. Diagnosing tunnelcollapse sections based on TBM tunneling big data and deep learning: a case study on the Yinsong project, China[J]. Tunnelling and Underground Space Technology, 2021(108): 103700.

    [4]劉詩(shī)洋,陳祖煜,張?jiān)旗?,?基于卷積神經(jīng)網(wǎng)絡(luò)對(duì)TBM塌方段的反演分析[J].固體力學(xué)學(xué)報(bào),2021,42(3):287-301.

    [5]HOU S, LIU Y. Early warning of tunnel collapse based on Adam-optimised long short-term memory network and TBM operation parameters[J]. Engineering Applications of Artificial Intelligence, 2022 (112): 104842.

    [6]SHIN H S, KWON Y C, JUNG Y S, et al. Methodology for quantitative hazard assessment for tunnel collapses based on case histories in Korea[J]. International Journal of Rock Mechanics and Mining Sciences, 2009(46): 1072-1087.

    [7]譚立偉,溫少卿,李開(kāi)華,等.基于改進(jìn)G2-EWM-TOPSIS模型的引水隧洞塌方風(fēng)險(xiǎn)評(píng)價(jià)[J].科學(xué)技術(shù)與工程,2023,23(30):13105-13111.

    [8]李燕,靳春玲,貢力,等.基于熵理論-可拓云模型的引水隧洞塌方風(fēng)險(xiǎn)評(píng)價(jià)[J].長(zhǎng)江科學(xué)院報(bào),2022,39(2):70-76,81.

    [9]呂擎峰,霍振升,趙本海,等.基于模糊層次和后果當(dāng)量法的隧道塌方風(fēng)險(xiǎn)評(píng)估[J].隧道建設(shè)(中英文),2018,38(S2):31-38.

    [10]LI L, NI B, ZHANG S, et al. Tunnel collapse risk assessment based on improved quantitative theory III and EW-AHP couplingweight[J]. Scientific Reports, 2022(12): 1600eKVGLh4rbgaz1zeVHb12is2C6jHIq+Fbkh5rvL+JYM=54.

    [11]OU G, JIAO Y, ZHANG G, et al. Collapse risk assessment of deep-buried tunnel during construction and its application[J]. Tunnelling and Underground Space Technology, 2021(115): 104019.

    [12]CHEN W, ZHANG G, JIAO Y, etal. Unascertained measure-set pair analysis model of collapse risk evaluation in mountain tunnels and its engineering application[J]. KSCE Journal of Civil Engineering, 2021(25):451-467.

    [13]BRIN S and PAGE L. The anatomy of a large-scalehypertextual Web search engine[J]. Computer Networks and ISDN Systems, 1998(30):107-117.

    [14]KLEINBERG J M. Authoritative sources in a hyperlinkedenvironment[J]. Journal of the ACM, 1998(46): 604-632.

    [15]王海英,楚林,許建.融合N-K模型的復(fù)雜網(wǎng)絡(luò)架橋機(jī)施工安全風(fēng)險(xiǎn)因素分析[J].中國(guó)安全科學(xué)學(xué)報(bào),2023,33(9):10-17.

    [16]張兵,蔡天昊,熊婷蜓.基于復(fù)雜網(wǎng)絡(luò)的基坑事故致因分析[J].項(xiàng)目管理技術(shù),2023,21(7):89-94.

    收稿日期:2024-08-15

    作者簡(jiǎn)介:

    張勇(1965—),男,博士,教授級(jí)高級(jí)工程師,碩士研究生導(dǎo)師,研究方向:土木建造與安全管理、項(xiàng)目管理。

    張琪(通信作者)(2001—),女,研究方向:水利工程安全風(fēng)險(xiǎn)管理。

    劉少乾(2000—),男,研究方向:項(xiàng)目管理。

    王沛豐(1999—),男,研究方向:項(xiàng)目管理。

    91在线精品国自产拍蜜月| 内射极品少妇av片p| 欧美三级亚洲精品| 久久精品久久久久久噜噜老黄 | 最近视频中文字幕2019在线8| 成人无遮挡网站| 精品久久久久久成人av| 欧美+日韩+精品| 日韩欧美 国产精品| 国产精品,欧美在线| 嫩草影院入口| 亚洲第一欧美日韩一区二区三区| 欧美日本视频| 男人舔奶头视频| 精品久久国产蜜桃| 久久久久九九精品影院| 久久伊人香网站| 欧美绝顶高潮抽搐喷水| 久久欧美精品欧美久久欧美| 免费高清视频大片| 久久久久九九精品影院| 99久久成人亚洲精品观看| 看十八女毛片水多多多| 精品久久久久久久久亚洲 | 最新在线观看一区二区三区| 少妇的逼水好多| 精品国产亚洲在线| 我的女老师完整版在线观看| 18美女黄网站色大片免费观看| 成人亚洲精品av一区二区| 国产精品亚洲美女久久久| 人人妻,人人澡人人爽秒播| 国产亚洲欧美在线一区二区| av女优亚洲男人天堂| 99热这里只有是精品50| 日韩中文字幕欧美一区二区| 亚洲综合色惰| 97热精品久久久久久| 久久伊人香网站| 网址你懂的国产日韩在线| 日韩中字成人| 两个人视频免费观看高清| 久久国产乱子免费精品| www日本黄色视频网| 国产一区二区在线av高清观看| 老鸭窝网址在线观看| 老司机午夜十八禁免费视频| 日韩欧美国产在线观看| netflix在线观看网站| 日韩中文字幕欧美一区二区| 久久精品久久久久久噜噜老黄 | 草草在线视频免费看| 亚洲精品色激情综合| 嫩草影院精品99| 精品一区二区三区视频在线| 人妻丰满熟妇av一区二区三区| 不卡一级毛片| 天堂动漫精品| 日韩亚洲欧美综合| 高潮久久久久久久久久久不卡| 欧美日韩黄片免| 亚洲av免费在线观看| 禁无遮挡网站| 淫妇啪啪啪对白视频| 999久久久精品免费观看国产| 88av欧美| 国产欧美日韩一区二区精品| 免费人成在线观看视频色| 一区二区三区高清视频在线| 变态另类丝袜制服| 成人av在线播放网站| 窝窝影院91人妻| 国产爱豆传媒在线观看| 国产精品日韩av在线免费观看| 可以在线观看的亚洲视频| 免费一级毛片在线播放高清视频| 国内少妇人妻偷人精品xxx网站| 99在线人妻在线中文字幕| 国产淫片久久久久久久久 | 国产91精品成人一区二区三区| 高清毛片免费观看视频网站| 国产高清三级在线| 国产一区二区三区视频了| 一二三四社区在线视频社区8| 欧美在线一区亚洲| 18禁裸乳无遮挡免费网站照片| 日韩中文字幕欧美一区二区| 日本免费一区二区三区高清不卡| 欧美在线黄色| 最近视频中文字幕2019在线8| 亚洲黑人精品在线| 男插女下体视频免费在线播放| av女优亚洲男人天堂| 久久久国产成人精品二区| 国产成人啪精品午夜网站| 日韩大尺度精品在线看网址| 亚洲va日本ⅴa欧美va伊人久久| 国产老妇女一区| 国产极品精品免费视频能看的| 此物有八面人人有两片| 国产精品av视频在线免费观看| 少妇的逼好多水| 亚洲七黄色美女视频| 国产美女午夜福利| 午夜久久久久精精品| 成年女人看的毛片在线观看| 九色成人免费人妻av| 狂野欧美白嫩少妇大欣赏| 韩国av一区二区三区四区| 欧美又色又爽又黄视频| 天美传媒精品一区二区| 日韩欧美在线二视频| 午夜两性在线视频| 一二三四社区在线视频社区8| 狂野欧美白嫩少妇大欣赏| 搡老熟女国产l中国老女人| 成人三级黄色视频| 色播亚洲综合网| 又爽又黄a免费视频| 午夜福利在线观看免费完整高清在 | 亚洲aⅴ乱码一区二区在线播放| 日本三级黄在线观看| 夜夜看夜夜爽夜夜摸| 久久人人爽人人爽人人片va | 午夜福利在线观看免费完整高清在 | 高清在线国产一区| 三级男女做爰猛烈吃奶摸视频| 丰满人妻一区二区三区视频av| 国产真实乱freesex| 国产亚洲精品av在线| 精品久久久久久久久亚洲 | av在线老鸭窝| 日本免费一区二区三区高清不卡| 日韩欧美免费精品| 亚洲精品乱码久久久v下载方式| 夜夜看夜夜爽夜夜摸| 亚洲精品粉嫩美女一区| 97超级碰碰碰精品色视频在线观看| 欧美日韩瑟瑟在线播放| 久久精品国产清高在天天线| 久久久国产成人免费| 人人妻人人看人人澡| 18禁黄网站禁片免费观看直播| 亚洲精华国产精华精| 99热这里只有是精品50| 亚洲精品在线观看二区| 91久久精品国产一区二区成人| 国产高清三级在线| 神马国产精品三级电影在线观看| 亚洲七黄色美女视频| 中文字幕人妻熟人妻熟丝袜美| 每晚都被弄得嗷嗷叫到高潮| 免费看美女性在线毛片视频| 丰满人妻熟妇乱又伦精品不卡| 精品一区二区三区人妻视频| 18禁在线播放成人免费| 亚洲久久久久久中文字幕| 精品国产亚洲在线| 我要看日韩黄色一级片| 99热这里只有精品一区| 亚洲第一区二区三区不卡| 亚洲av电影在线进入| 特大巨黑吊av在线直播| 国产午夜福利久久久久久| 舔av片在线| 亚洲人成网站在线播| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人精品一区久久| 国产精品久久久久久亚洲av鲁大| 啪啪无遮挡十八禁网站| 国产色婷婷99| 亚洲国产精品999在线| 18禁裸乳无遮挡免费网站照片| 麻豆一二三区av精品| 婷婷色综合大香蕉| 亚洲av免费高清在线观看| 国产亚洲精品av在线| 国产免费av片在线观看野外av| 一级作爱视频免费观看| 亚洲欧美日韩东京热| 91麻豆精品激情在线观看国产| 欧美日韩福利视频一区二区| 天天一区二区日本电影三级| 欧美黑人欧美精品刺激| 欧美丝袜亚洲另类 | 日韩欧美一区二区三区在线观看| 亚洲成人久久爱视频| 欧美日韩瑟瑟在线播放| 精品午夜福利在线看| 日本黄色片子视频| 国产精品久久久久久久电影| 听说在线观看完整版免费高清| 亚洲第一电影网av| 又粗又爽又猛毛片免费看| 成人午夜高清在线视频| 国产精品电影一区二区三区| 婷婷色综合大香蕉| 欧美色视频一区免费| 国产免费男女视频| 国产高清视频在线观看网站| 免费看日本二区| 亚洲电影在线观看av| 极品教师在线免费播放| 国产aⅴ精品一区二区三区波| 波多野结衣高清无吗| 亚洲五月天丁香| 国产激情偷乱视频一区二区| av欧美777| 精品人妻熟女av久视频| 欧美不卡视频在线免费观看| 久久久久久国产a免费观看| 99热这里只有精品一区| 国内少妇人妻偷人精品xxx网站| 欧美在线一区亚洲| 真人做人爱边吃奶动态| 久久亚洲精品不卡| 国产精品美女特级片免费视频播放器| 看十八女毛片水多多多| 国产精品爽爽va在线观看网站| 久久伊人香网站| 神马国产精品三级电影在线观看| 一边摸一边抽搐一进一小说| 日日摸夜夜添夜夜添av毛片 | 欧美乱妇无乱码| 内射极品少妇av片p| 天堂动漫精品| 亚洲久久久久久中文字幕| 91字幕亚洲| 久久午夜福利片| 欧美高清成人免费视频www| 精品一区二区免费观看| 麻豆av噜噜一区二区三区| 国产极品精品免费视频能看的| 51午夜福利影视在线观看| 啦啦啦观看免费观看视频高清| 97人妻精品一区二区三区麻豆| 久久午夜亚洲精品久久| 国产毛片a区久久久久| 欧美一区二区国产精品久久精品| 国产av一区在线观看免费| 日本五十路高清| 一本久久中文字幕| 日韩国内少妇激情av| 精品久久久久久久久av| 黄色丝袜av网址大全| 精品人妻1区二区| 禁无遮挡网站| 婷婷亚洲欧美| 国产一区二区三区在线臀色熟女| 日本一本二区三区精品| 亚洲精品成人久久久久久| 久99久视频精品免费| 亚洲色图av天堂| 久久久精品欧美日韩精品| 国产又黄又爽又无遮挡在线| 国产真实伦视频高清在线观看 | 日日摸夜夜添夜夜添小说| 国产精品av视频在线免费观看| xxxwww97欧美| 三级国产精品欧美在线观看| 中亚洲国语对白在线视频| 中文资源天堂在线| 哪里可以看免费的av片| 直男gayav资源| 午夜福利18| 1024手机看黄色片| 亚洲专区国产一区二区| 亚洲av第一区精品v没综合| 99久久精品一区二区三区| 久久精品国产清高在天天线| 精品午夜福利视频在线观看一区| 午夜福利在线在线| 国产淫片久久久久久久久 | 两个人视频免费观看高清| 欧美乱妇无乱码| 伦理电影大哥的女人| 亚洲性夜色夜夜综合| 亚洲在线观看片| av欧美777| 麻豆av噜噜一区二区三区| 日韩欧美 国产精品| 国产成人影院久久av| 在线十欧美十亚洲十日本专区| 亚洲av日韩精品久久久久久密| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 成人av在线播放网站| 小说图片视频综合网站| 丝袜美腿在线中文| 一本综合久久免费| 午夜福利免费观看在线| 欧美最黄视频在线播放免费| 九九久久精品国产亚洲av麻豆| 又爽又黄无遮挡网站| 又粗又爽又猛毛片免费看| 成年人黄色毛片网站| 波野结衣二区三区在线| 91在线精品国自产拍蜜月| 亚洲国产日韩欧美精品在线观看| 精品久久久久久久久av| 亚洲av成人精品一区久久| 久久99热这里只有精品18| 久久草成人影院| 国产精品久久久久久亚洲av鲁大| 一级作爱视频免费观看| 欧美三级亚洲精品| 亚洲精品久久国产高清桃花| 国产乱人伦免费视频| 51午夜福利影视在线观看| 在线播放无遮挡| 国产一区二区激情短视频| 久久久久久久亚洲中文字幕 | 成人精品一区二区免费| 国产欧美日韩精品亚洲av| 乱码一卡2卡4卡精品| 搡女人真爽免费视频火全软件 | 老鸭窝网址在线观看| 免费看光身美女| 俄罗斯特黄特色一大片| 久久久久久久久中文| 床上黄色一级片| 天堂网av新在线| 国产私拍福利视频在线观看| www.熟女人妻精品国产| 久久这里只有精品中国| 91麻豆av在线| 久久精品夜夜夜夜夜久久蜜豆| 桃色一区二区三区在线观看| 欧美+亚洲+日韩+国产| 真人一进一出gif抽搐免费| 九九在线视频观看精品| 欧美成人一区二区免费高清观看| 一个人看的www免费观看视频| 日本a在线网址| 国产午夜精品论理片| 九九热线精品视视频播放| 99精品在免费线老司机午夜| 波多野结衣高清作品| 久久人妻av系列| 免费av毛片视频| 亚洲美女搞黄在线观看 | 国产av在哪里看| 无人区码免费观看不卡| 在线观看av片永久免费下载| 国产精品美女特级片免费视频播放器| 91久久精品电影网| 在现免费观看毛片| 身体一侧抽搐| 在现免费观看毛片| 精品一区二区三区人妻视频| 精品欧美国产一区二区三| 午夜两性在线视频| 男插女下体视频免费在线播放| 日本成人三级电影网站| www.www免费av| 日本成人三级电影网站| 午夜老司机福利剧场| 麻豆成人av在线观看| 国产精品久久久久久人妻精品电影| 一个人看视频在线观看www免费| 国产精品自产拍在线观看55亚洲| 国产精品不卡视频一区二区 | 亚洲av二区三区四区| 99久久无色码亚洲精品果冻| 国产白丝娇喘喷水9色精品| 免费人成视频x8x8入口观看| 搞女人的毛片| 51午夜福利影视在线观看| 久久精品国产亚洲av涩爱 | 全区人妻精品视频| 久久性视频一级片| 精品一区二区三区人妻视频| 99精品久久久久人妻精品| 国产精品综合久久久久久久免费| 国产人妻一区二区三区在| 日本免费一区二区三区高清不卡| 三级毛片av免费| 精品一区二区三区人妻视频| 成人av在线播放网站| 久久久久久九九精品二区国产| 成人性生交大片免费视频hd| 91字幕亚洲| 国产中年淑女户外野战色| 亚洲精品在线观看二区| 欧美黄色淫秽网站| 搡老岳熟女国产| 99久国产av精品| 婷婷色综合大香蕉| 国产成人啪精品午夜网站| 久久婷婷人人爽人人干人人爱| 别揉我奶头 嗯啊视频| 丰满的人妻完整版| 乱码一卡2卡4卡精品| 3wmmmm亚洲av在线观看| 成人性生交大片免费视频hd| av天堂中文字幕网| 我要搜黄色片| 黄色配什么色好看| 久久午夜福利片| 亚洲成人免费电影在线观看| 日韩亚洲欧美综合| 又粗又爽又猛毛片免费看| 一级黄片播放器| 2021天堂中文幕一二区在线观| 性欧美人与动物交配| 亚洲自拍偷在线| 婷婷丁香在线五月| 久久99热这里只有精品18| 国产aⅴ精品一区二区三区波| 88av欧美| 精品久久久久久久人妻蜜臀av| 最近视频中文字幕2019在线8| 欧美乱色亚洲激情| 欧美在线一区亚洲| 中文亚洲av片在线观看爽| 青草久久国产| 在线观看av片永久免费下载| 夜夜夜夜夜久久久久| 日韩成人在线观看一区二区三区| 免费黄网站久久成人精品 | 亚洲av中文字字幕乱码综合| 俄罗斯特黄特色一大片| 日韩精品中文字幕看吧| 久久天躁狠狠躁夜夜2o2o| 成人特级黄色片久久久久久久| 欧美乱妇无乱码| 日本免费一区二区三区高清不卡| 日本黄色片子视频| 男人舔女人下体高潮全视频| 国产 一区 欧美 日韩| АⅤ资源中文在线天堂| 亚洲成av人片在线播放无| 午夜亚洲福利在线播放| 欧美最新免费一区二区三区 | 欧美最新免费一区二区三区 | 一级黄色大片毛片| 欧美成狂野欧美在线观看| 亚洲专区中文字幕在线| 又爽又黄无遮挡网站| 久久精品夜夜夜夜夜久久蜜豆| 级片在线观看| 国产黄a三级三级三级人| 亚洲午夜理论影院| 狠狠狠狠99中文字幕| 国产精品1区2区在线观看.| 一本久久中文字幕| 久久精品人妻少妇| 久久久国产成人精品二区| 亚洲精品色激情综合| 久久精品国产亚洲av天美| 免费电影在线观看免费观看| 波多野结衣巨乳人妻| 日本免费a在线| 精品福利观看| 久久久国产成人免费| av中文乱码字幕在线| 色综合欧美亚洲国产小说| 九九久久精品国产亚洲av麻豆| 国产成人a区在线观看| 国产精品1区2区在线观看.| 青草久久国产| 日本五十路高清| 欧美另类亚洲清纯唯美| 99热只有精品国产| 亚洲午夜理论影院| 老鸭窝网址在线观看| 日本黄大片高清| x7x7x7水蜜桃| 少妇丰满av| av在线老鸭窝| 国产精品精品国产色婷婷| 国产精品国产高清国产av| 禁无遮挡网站| 亚洲精品一区av在线观看| 亚洲人成网站在线播| 欧美黑人欧美精品刺激| 好看av亚洲va欧美ⅴa在| 性插视频无遮挡在线免费观看| av在线蜜桃| 99久久精品国产亚洲精品| 亚洲片人在线观看| 免费观看精品视频网站| 嫩草影院新地址| 九九在线视频观看精品| 非洲黑人性xxxx精品又粗又长| 欧美bdsm另类| 日韩大尺度精品在线看网址| 国产伦人伦偷精品视频| 在线天堂最新版资源| 狂野欧美白嫩少妇大欣赏| 中文字幕人成人乱码亚洲影| 国产色婷婷99| 亚洲七黄色美女视频| 一区二区三区激情视频| 很黄的视频免费| 午夜福利欧美成人| 看免费av毛片| 日日干狠狠操夜夜爽| 日韩中文字幕欧美一区二区| 99热这里只有是精品在线观看 | 精品国内亚洲2022精品成人| 午夜亚洲福利在线播放| 欧美三级亚洲精品| 欧美黄色淫秽网站| 日韩欧美三级三区| 国产伦一二天堂av在线观看| 亚洲欧美清纯卡通| 国产三级在线视频| 中文在线观看免费www的网站| 九九热线精品视视频播放| 免费av观看视频| 69人妻影院| 永久网站在线| 久久草成人影院| 欧美黑人欧美精品刺激| 搞女人的毛片| 在现免费观看毛片| 亚洲成人中文字幕在线播放| 夜夜躁狠狠躁天天躁| 一级毛片久久久久久久久女| 动漫黄色视频在线观看| 又爽又黄无遮挡网站| 亚州av有码| 搡老熟女国产l中国老女人| 亚洲精品久久国产高清桃花| av中文乱码字幕在线| 日韩欧美国产在线观看| 欧美性猛交╳xxx乱大交人| 欧美区成人在线视频| 很黄的视频免费| 亚洲综合色惰| 亚洲无线观看免费| 久久性视频一级片| 美女大奶头视频| 国产色婷婷99| 久久久久久久精品吃奶| 日韩免费av在线播放| 蜜桃久久精品国产亚洲av| 国产一区二区在线av高清观看| 婷婷亚洲欧美| 热99re8久久精品国产| 国产伦精品一区二区三区视频9| 赤兔流量卡办理| 国产精品久久视频播放| 日本五十路高清| 国产aⅴ精品一区二区三区波| 俄罗斯特黄特色一大片| 精品国产亚洲在线| 亚洲第一欧美日韩一区二区三区| 哪里可以看免费的av片| 国产白丝娇喘喷水9色精品| 亚洲人与动物交配视频| 他把我摸到了高潮在线观看| 日日摸夜夜添夜夜添av毛片 | 亚洲欧美日韩高清在线视频| 亚洲七黄色美女视频| 91久久精品电影网| 午夜福利成人在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲经典国产精华液单 | 久久人人精品亚洲av| 老司机福利观看| 他把我摸到了高潮在线观看| 久久久久免费精品人妻一区二区| 久久久久九九精品影院| 欧美一区二区精品小视频在线| 欧美色欧美亚洲另类二区| 国内揄拍国产精品人妻在线| 免费av观看视频| 免费观看的影片在线观看| 少妇人妻一区二区三区视频| 成人一区二区视频在线观看| 亚洲精华国产精华精| 亚洲真实伦在线观看| 99热这里只有是精品50| 女人十人毛片免费观看3o分钟| 精品人妻偷拍中文字幕| 别揉我奶头 嗯啊视频| 嫩草影院精品99| 精品一区二区三区视频在线观看免费| 蜜桃亚洲精品一区二区三区| 有码 亚洲区| 久久国产精品影院| 久久久国产成人精品二区| 亚洲内射少妇av| 变态另类成人亚洲欧美熟女| 麻豆av噜噜一区二区三区| 长腿黑丝高跟| 亚洲av二区三区四区| 色视频www国产| 天堂√8在线中文| 国产精品1区2区在线观看.| 亚洲成a人片在线一区二区| 波多野结衣高清作品| 国产黄片美女视频| 麻豆av噜噜一区二区三区| 哪里可以看免费的av片| 波野结衣二区三区在线| 国产精品亚洲av一区麻豆| 国产一区二区在线av高清观看| 国产老妇女一区| 亚洲第一欧美日韩一区二区三区| 久久精品综合一区二区三区| 国产亚洲精品久久久久久毛片| 又爽又黄无遮挡网站| 国产三级黄色录像| 亚洲中文字幕日韩| 欧美极品一区二区三区四区| 九九热线精品视视频播放| 动漫黄色视频在线观看| 免费黄网站久久成人精品 | 国产亚洲精品久久久久久毛片| 极品教师在线视频| 丰满乱子伦码专区| 成年女人看的毛片在线观看| 国产精品久久电影中文字幕|