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

    礦山尾礦庫(kù)風(fēng)險(xiǎn)預(yù)測(cè)時(shí)空模型的研究

    2016-06-29 05:28:12顏世航洑佳紅
    關(guān)鍵詞:潰壩尾礦庫(kù)矩陣

    顏世航,洑佳紅

    (浙江工業(yè)大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,浙江 杭州 310023)

    礦山尾礦庫(kù)風(fēng)險(xiǎn)預(yù)測(cè)時(shí)空模型的研究

    顏世航,洑佳紅

    (浙江工業(yè)大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,浙江 杭州 310023)

    摘要:尾礦庫(kù)風(fēng)險(xiǎn)預(yù)測(cè)是增強(qiáng)尾礦庫(kù)安全管理水平和提高突發(fā)應(yīng)對(duì)能力的重要手段之一.針對(duì)尾礦庫(kù)內(nèi)部結(jié)構(gòu)復(fù)雜,風(fēng)險(xiǎn)預(yù)測(cè)指標(biāo)沒(méi)有明確界定,穩(wěn)定性評(píng)價(jià)及預(yù)警手段尚未科學(xué)量化等問(wèn)題,從優(yōu)化尾礦庫(kù)安全監(jiān)測(cè)指標(biāo),研究監(jiān)測(cè)數(shù)據(jù)變化規(guī)律入手,建立多層次數(shù)據(jù)分析模型,進(jìn)而定量評(píng)估尾礦庫(kù)風(fēng)險(xiǎn)狀況和實(shí)時(shí)預(yù)測(cè)尾礦庫(kù)安全發(fā)展趨勢(shì),最后通過(guò)實(shí)例驗(yàn)證了模型的有效性.

    關(guān)鍵詞:地質(zhì)災(zāi)害;風(fēng)險(xiǎn)評(píng)估;尾礦庫(kù);時(shí)空信息

    截至2014年底,全國(guó)共有尾礦庫(kù)11 359座,2014年全國(guó)尾礦庫(kù)共發(fā)生事故8起,其中4起造成環(huán)境污染,4起造成5人死亡[1].我國(guó)尾礦庫(kù)普遍存在建設(shè)基礎(chǔ)薄弱,管理水平低下等問(wèn)題,在安全性研究以及監(jiān)測(cè)預(yù)警方面有許多不足,主要原因是尾礦庫(kù)自身結(jié)構(gòu)的復(fù)雜性;尾礦庫(kù)筑壩方式、材料特性影響傳統(tǒng)力學(xué)模型的建立以及行業(yè)技術(shù)標(biāo)準(zhǔn)對(duì)尾礦庫(kù)安全性研究未給出明確監(jiān)測(cè)指標(biāo).當(dāng)前學(xué)者們?cè)谖驳V庫(kù)風(fēng)險(xiǎn)預(yù)測(cè)方面開(kāi)展了大量的研究工作,其中比較典型的成果有:1) 潰壩因素研究,其中Rico M等[2]在分析全球147座尾礦庫(kù)潰壩原因后提出尾礦庫(kù)潰壩事故與壩高有關(guān),Chakraborty D等[3]研究了地震對(duì)尾礦庫(kù)的變形狀況,Mahmood K等[4]認(rèn)為地震將造成尾礦液化導(dǎo)致潰壩事故,研究發(fā)現(xiàn)僅分析單一因素對(duì)尾礦風(fēng)險(xiǎn)影響還遠(yuǎn)遠(yuǎn)不夠.2) 安全監(jiān)測(cè)技術(shù),李忠奎等[5]將浸潤(rùn)線(xiàn)、壩體位移和降雨量作為尾礦庫(kù)主要安全指標(biāo)設(shè)計(jì)潰壩預(yù)警系統(tǒng),謝旭陽(yáng)等[6-7]利用支持向量機(jī)等機(jī)器學(xué)習(xí)方法建立尾礦庫(kù)區(qū)域預(yù)警模型,趙志軍等[8]將ZigBee無(wú)線(xiàn)傳感網(wǎng)絡(luò)技術(shù)應(yīng)用于尾礦大壩安全監(jiān)測(cè)系統(tǒng),利用新型技術(shù)監(jiān)測(cè)各類(lèi)指標(biāo)來(lái)提高尾礦監(jiān)測(cè)水平,但在數(shù)據(jù)后續(xù)處理上略顯不足.3) 壩體穩(wěn)定性評(píng)價(jià),李全明等[9]提出了層次結(jié)構(gòu)模型得到各指標(biāo)權(quán)重,羅飛飛等[10]提出模糊灰色綜合評(píng)價(jià)法解決專(zhuān)家評(píng)價(jià)在較小分?jǐn)?shù)差異對(duì)象上誤差較大的問(wèn)題,鄭瑞等[11]提出可用信息熵客觀計(jì)算權(quán)重從而建立模糊評(píng)價(jià)模型.

    綜合以上分析,從現(xiàn)有的尾礦庫(kù)安全研究狀況來(lái)看,要完善尾礦庫(kù)風(fēng)險(xiǎn)預(yù)測(cè)還需進(jìn)行以下工作:1)尾礦庫(kù)諸多致災(zāi)因素對(duì)尾礦庫(kù)潰壩影響程度需科學(xué)量化;2) 建立多維多層次的評(píng)價(jià)模型,并在模型中增加由于致災(zāi)因素的數(shù)據(jù)采集節(jié)點(diǎn)空間位置不同對(duì)尾礦庫(kù)穩(wěn)定性的影響;3) 能實(shí)時(shí)依據(jù)尾礦庫(kù)安全現(xiàn)狀對(duì)未來(lái)趨勢(shì)做出預(yù)警判斷.采用在尾礦庫(kù)分析過(guò)程中能全面反映尾礦庫(kù)各指標(biāo)空間信息及安全狀況的層次分析法,建立空間層次評(píng)價(jià)模型并結(jié)合能將指標(biāo)定性的安全狀況轉(zhuǎn)化成定量評(píng)估分值的模糊隸屬度函數(shù)對(duì)目標(biāo)尾礦庫(kù)的安全現(xiàn)狀作出準(zhǔn)確快速評(píng)析,科學(xué)量化當(dāng)前尾礦庫(kù)安全狀況并實(shí)時(shí)作出風(fēng)險(xiǎn)預(yù)測(cè).

    1空間層次評(píng)價(jià)模型

    基于應(yīng)用網(wǎng)絡(luò)系統(tǒng)理論和多目標(biāo)綜合評(píng)價(jià)方法,美國(guó)運(yùn)籌學(xué)家薩蒂提出了層次分析法,該方法將影響決策的各類(lèi)元素分解成目標(biāo)、標(biāo)準(zhǔn)和方案等層次,再進(jìn)一步展開(kāi)定性和定量分析[12].空間層次評(píng)價(jià)模型的評(píng)價(jià)原則如下:首先將決策問(wèn)題分級(jí)系列化,根據(jù)決策目標(biāo)自身性質(zhì)和決策要求,將決策問(wèn)題分解歸納為各類(lèi)影響因素,確立因素間的互相影響關(guān)系和隸屬關(guān)系將其分級(jí)聚類(lèi),在層次評(píng)價(jià)模型的基礎(chǔ)上增加空間影響因素,將空間因素置于層次評(píng)價(jià)模型最底層,形成一個(gè)遞階、有序的層次結(jié)構(gòu)模型,然后判斷空間層次模型中每一層因素以及各層因素之間的相互影響力,確定該層全部因素對(duì)上一層因素相對(duì)重要性的權(quán)值,最后從最底層開(kāi)始通過(guò)綜合計(jì)算,逐層往上,得到最終目標(biāo)問(wèn)題的評(píng)價(jià)結(jié)果[13].該層次空間模型涉及到對(duì)目標(biāo)決策的模型構(gòu)造和空間因素預(yù)測(cè),基于空間數(shù)據(jù)分析和空間統(tǒng)計(jì)分析進(jìn)行區(qū)位預(yù)測(cè)和優(yōu)化.

    1.1空間數(shù)據(jù)分析

    尾礦庫(kù)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)的選取是影響尾礦庫(kù)風(fēng)險(xiǎn)評(píng)估預(yù)測(cè)準(zhǔn)信度的重要因素.各類(lèi)風(fēng)險(xiǎn)指標(biāo)中,尾礦庫(kù)的空間評(píng)價(jià)指標(biāo)除了在空間上有其定位特性,指標(biāo)間的空間關(guān)系還具備復(fù)雜的連接屬性.其中的作用關(guān)系體現(xiàn)在尾礦庫(kù)空間指標(biāo)自相關(guān)特點(diǎn)和與之相伴隨的可變區(qū)域單位問(wèn)題、尺度和邊緣效應(yīng).通過(guò)空間指標(biāo)構(gòu)建空間模型后,通過(guò)Geary系數(shù)進(jìn)行空間自相關(guān)檢驗(yàn),以此保證各個(gè)空間指標(biāo)間相互獨(dú)立性.各指標(biāo)間的相關(guān)系數(shù)在[0,2]之間取值且用γ表示,γ<1表示正相關(guān),γ>1表示負(fù)相關(guān),γ=1表示不相關(guān).各指標(biāo)間進(jìn)行相關(guān)性計(jì)算,認(rèn)定相關(guān)性較大的風(fēng)險(xiǎn)指標(biāo)為同一因素而舍去,最終建立尾礦庫(kù)風(fēng)險(xiǎn)評(píng)估指標(biāo)體系如圖1所示的多級(jí)遞階層次結(jié)構(gòu)模型.該模型層次結(jié)構(gòu)分為P,U,R,K四級(jí).其中P級(jí)為目標(biāo)尾礦庫(kù);U級(jí)為第二級(jí),選取地下位移,庫(kù)水位等影響尾礦庫(kù)安全的六個(gè)主要因素;R級(jí)分別選取影響地下位移和地下傾斜角的三個(gè)不同空間位置的因素;K級(jí)為第四級(jí),將代表不同空間位置的傳感器作為第四級(jí)因素,增加了空間信息對(duì)尾礦庫(kù)的影響.

    圖1 尾礦庫(kù)風(fēng)險(xiǎn)評(píng)估指標(biāo)Fig.1 Tailings risk assessment index system

    1.2模型建立

    設(shè)有尾礦庫(kù)有A1,A2,…,An共n個(gè)安全指標(biāo),它們的重要程度可表示成重要向 量w=(a1,a2,…,an)T,若將它們兩兩比較重量,其比值可構(gòu)成n×n的矩陣A,即

    為了因素間兩兩比較得到科學(xué)量化的判斷矩陣,根據(jù)人們區(qū)分信息等級(jí)的極限能力采用1~9比例標(biāo)度法(表1)對(duì)尾礦庫(kù)指標(biāo)兩兩比較建立評(píng)判矩陣,用aij表示因素ai與因素aj相對(duì)隸屬目標(biāo)的影響程度之比,通常用數(shù)字1~9及其倒數(shù)作為程度比較的標(biāo)度,標(biāo)度aij={2,4,6,8,1/2,1/4,1/6,1/8}表示重要性等級(jí)在aij={1,3,5,7,9,1/3,1/5,1/7,1/9}之間[14].

    表1 判斷矩陣標(biāo)度及含義

    尾礦庫(kù)與二級(jí)指標(biāo)構(gòu)成的模糊評(píng)判矩陣為

    二級(jí)指標(biāo)地下位移與三級(jí)指標(biāo)地下位移上部、中部及下部構(gòu)成的模糊評(píng)判矩陣為

    二級(jí)指標(biāo)地下傾斜角與三級(jí)指標(biāo)地下傾斜角上部、中部及下部構(gòu)成的模糊評(píng)判矩陣為

    四級(jí)空間指標(biāo)與上一級(jí)指標(biāo)間構(gòu)成的模糊評(píng)判矩陣為

    運(yùn)用方根法近似計(jì)算求解得到矩陣特征向量W,即

    (1)

    (2)

    其中 aij為n維評(píng)判矩陣A的第i行第j列元素.

    當(dāng)各層的諸因素的相對(duì)權(quán)重都得到后,進(jìn)行組合權(quán)重計(jì)算.設(shè)有由P,U,R,K構(gòu)成的四層層次模型.目標(biāo)尾礦庫(kù)對(duì)U層的相對(duì)權(quán)重為

    U層的各指標(biāo)Ui對(duì)R層的n個(gè)指標(biāo)間的相對(duì)權(quán)重為

    1.3一致性檢驗(yàn)

    依據(jù)矩陣?yán)碚?,利用模糊評(píng)判矩陣進(jìn)行決策時(shí),在保持元素間兩兩比較序的傳遞性又要保持整體偏差的可接受性,因此定義了計(jì)算滿(mǎn)意一致性矩陣的概念[15],將矩陣一致性檢驗(yàn)公式定義為

    (3)

    式中:RI為平均隨機(jī)一致性指標(biāo),判斷矩陣的維數(shù)越大,判斷矩陣的一致性將越差,矩陣的隨機(jī)一致性指標(biāo)可以通過(guò)表2得到;(AW)i為向量的第i個(gè)元素;A為判斷矩陣;CR為隨機(jī)一致性比率.一般,當(dāng)CR小于0.10時(shí),認(rèn)為A的不一致程度在容許范圍內(nèi),有滿(mǎn)意的一致性,可用其歸一化特征向量作為權(quán)重向量,否則要重新構(gòu)造判斷矩陣,對(duì)aij加以調(diào)整[16].

    表2 1~9階矩陣隨機(jī)一致性指標(biāo)

    運(yùn)用馬氏距離衡量某一尾礦庫(kù)中所有同級(jí)評(píng)價(jià)權(quán)重向量之間的差異,以此削減專(zhuān)家確定判斷矩陣造成的主觀性,對(duì)于l個(gè)專(zhuān)家給出的判斷矩陣,根據(jù)式(1,2)得到的特征向量W1~Wl,確立最合理向量Wk作為最終的評(píng)價(jià)權(quán)重.協(xié)方差矩陣記為S,向量Wi與Wj之間的馬氏距離定義為

    (4)

    (5)

    協(xié)方差矩陣中每個(gè)元素是各個(gè)矢量元素之間的協(xié)方差表示為Cov(X,Y),E為數(shù)學(xué)期望,Cov(X,Y)=E{[X-E(X)][Y-E(Y)]}.

    2模糊綜合評(píng)價(jià)

    模糊數(shù)學(xué)是由美國(guó)自動(dòng)控制專(zhuān)家L.A.Zadeh提出的一種分析復(fù)雜系統(tǒng)的新方法,模糊綜合評(píng)價(jià)是模糊數(shù)學(xué)的一種具體應(yīng)用方法[17].用模糊數(shù)學(xué)對(duì)受到多種因素制約的事物或?qū)ο笞龀鲆粋€(gè)總體的評(píng)價(jià),能較好地解決模糊的、難以量化的問(wèn)題,適合各種非確定性問(wèn)題的解決[18].為了減少人為因素對(duì)尾礦庫(kù)風(fēng)險(xiǎn)評(píng)估指標(biāo)的主觀評(píng)定,運(yùn)用模糊變換原理和最大隸屬度原則,考慮因素集 U={u1,u2,…,un}對(duì)尾礦庫(kù)總體風(fēng)險(xiǎn)狀況相關(guān)程度來(lái)對(duì)其作出判斷集 V={v1,v2,…,vn}.對(duì)單個(gè)因素 ui(i=1,2,…,n)的評(píng)判,得到V上的模糊集 (ri1,ri2,…,rin),所以它是從 U到 V的一個(gè)模糊映射.

    2.1隸屬度函數(shù)

    若對(duì)論域 U中的任一元素x,都有一個(gè)數(shù)F(x)∈[0,1]與之對(duì)應(yīng),則稱(chēng)F為論域上的模糊集,F(xiàn)(x)稱(chēng)為x對(duì)F的隸屬度[19].用取值于區(qū)間[0,1]的隸屬函數(shù)F(x)表征x屬于F的程度高低,以此通過(guò)建立隸屬度函數(shù)取代專(zhuān)家打分來(lái)確立尾礦庫(kù)指標(biāo)分值.尾礦庫(kù)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)數(shù)量較多且自身的模糊特性也不相同,庫(kù)水位、地表位移、地下位移、地下傾斜角和地下水位等五個(gè)指標(biāo)的評(píng)價(jià)值應(yīng)隨屬性值的減小而增大,同時(shí)指標(biāo)干灘長(zhǎng)度的評(píng)價(jià)值應(yīng)隨著屬性值的減小而減小.尾礦庫(kù)指標(biāo)一年內(nèi)的數(shù)據(jù)值集合含有大量的數(shù)值且傳感器采樣頻率較高和指標(biāo)特性使得集合內(nèi)的值等同于連續(xù),因此采用計(jì)算簡(jiǎn)單、高效的一般分段的線(xiàn)性隸屬度函數(shù),即

    (6)

    式中:X(t)定義為t時(shí)刻傳感器采集的數(shù)據(jù);S(t)為t時(shí)刻指標(biāo)評(píng)價(jià)分值;a, b分別代表尾礦庫(kù)采集時(shí)間范圍內(nèi)該指標(biāo)數(shù)據(jù)量的最小、最大值,a=minX(t),b=maxX(t).

    2.2綜合評(píng)價(jià)

    在復(fù)雜的系統(tǒng)中,由于因素較多且各因素之間還有層次之分,為了更好地比較系統(tǒng)中事務(wù)間的優(yōu)劣次序,得到有意義的評(píng)判結(jié)果,還需要進(jìn)行多級(jí)綜合評(píng)判.尾礦庫(kù)評(píng)價(jià)標(biāo)準(zhǔn)采用五級(jí)評(píng)語(yǔ),將尾礦庫(kù)安全風(fēng)險(xiǎn)預(yù)測(cè)的分值與表3對(duì)應(yīng),實(shí)時(shí)給出安全風(fēng)險(xiǎn)等級(jí).每一級(jí)綜合評(píng)價(jià)式為

    (7)

    表3 安全風(fēng)險(xiǎn)等級(jí)表

    3實(shí)例分析

    采用浙江省淳安縣某尾礦庫(kù)2013年11月至2014年7月間10 類(lèi)風(fēng)險(xiǎn)指標(biāo)傳感器采集的400 條數(shù)據(jù)作為測(cè)試樣本,并運(yùn)用數(shù)據(jù)預(yù)處理技術(shù)對(duì)少量缺失數(shù)據(jù)填寫(xiě)空缺值,平滑噪聲數(shù)據(jù).

    3.1建立判斷矩陣

    確定尾礦庫(kù)各級(jí)因素之間的相對(duì)重要性,得到U,R,K級(jí)九級(jí)標(biāo)度法,10 位專(zhuān)家對(duì)每一層的致災(zāi)因素進(jìn)行衡量評(píng)判并寫(xiě)出所有評(píng)判矩陣式,即

    式中:l= 1,2,…, 10;當(dāng)i= j時(shí),uij,mij,nij表示單一因素對(duì)自身的重要性,依據(jù)判斷矩陣標(biāo)度其值為1.根據(jù)式(1,2)計(jì)算特征向量并利用馬氏距離求出最合適的特征向量作為最終的權(quán)重向量.

    3.2計(jì)算權(quán)重

    利用馬氏距離選擇最合適的特征向量作為權(quán)重向量,得

    3.3矩陣一致性檢驗(yàn)

    通過(guò)查表2獲得相應(yīng)階數(shù)矩陣的隨機(jī)一致性指標(biāo),將最終的權(quán)重向量所對(duì)應(yīng)的評(píng)判矩陣代入式(3)計(jì)算,分別得到隨機(jī)一致性比率CRu=0.043,CRRi=0.025,CRRj=0.017,CRk=0.019,其值都小于0.1,均滿(mǎn)足矩陣一致性要求.

    3.4尾礦庫(kù)時(shí)空信息

    依據(jù)當(dāng)前尾礦庫(kù)風(fēng)險(xiǎn)情況的判斷預(yù)測(cè)未來(lái)尾礦庫(kù)發(fā)展趨勢(shì).選取時(shí)間、空間位置及地表位移評(píng)價(jià)分值建立尾礦庫(kù)指標(biāo)評(píng)價(jià)三維圖(圖2),選取時(shí)間,地表位移綜合評(píng)價(jià)分值建立地表位移歷史安全狀況圖(圖3).從圖2中觀察單一傳感器的安全狀況隨時(shí)間變化情況發(fā)現(xiàn),6個(gè)空間傳感器分別反映地表位移指標(biāo)的安全狀況隨時(shí)間波動(dòng)較大且彼此間的狀況也很不一致,說(shuō)明尾礦庫(kù)在不同地表位置移動(dòng)情況不同,不能根據(jù)單一指標(biāo)的快速移動(dòng)或是單一指標(biāo)的相對(duì)穩(wěn)定而判斷出尾礦庫(kù)的安全狀況.而通過(guò)圖3觀察發(fā)現(xiàn)從總體趨勢(shì)上來(lái)說(shuō)地表位移總體態(tài)勢(shì)良好,且安全狀況趨于上升.在2013年12月1日這一時(shí)間點(diǎn)開(kāi)始,根據(jù)淳安縣12月未來(lái)15天的天氣狀況判斷尾礦庫(kù)將會(huì)由于大面積降雪而導(dǎo)致地表勢(shì)能增加而擴(kuò)大潰壩風(fēng)險(xiǎn).在查閱淳安縣2013年歷史天氣狀況后發(fā)現(xiàn)淳安縣在2013年12月中旬有大面積降雪降雨,地表位移安全狀況明顯下降,造成了圖3中2013年12月地表位移綜合評(píng)價(jià)分值的明顯降低.

    圖2 地表位移時(shí)空信息表Fig.2 Information about surface displacement condition

    圖3 地表位移歷史狀況Fig.3 Surface displacement safety history

    圖4反映了尾礦庫(kù)安全狀況時(shí)空綜合評(píng)價(jià),發(fā)現(xiàn)2013年12月及2014年夏季該尾礦庫(kù)安全狀況出現(xiàn)了兩次明顯的下降,在對(duì)所有指標(biāo)進(jìn)行時(shí)空分析后得出:庫(kù)水位的升高是造成安全趨勢(shì)下降的主要原因,圖5所示7號(hào)傳感器所處的尾礦庫(kù)位置由于水位上升增加了潰壩風(fēng)險(xiǎn)傳感器.查閱尾礦庫(kù)當(dāng)?shù)靥鞖鉅顩r發(fā)現(xiàn)在12月及夏季,淳安縣降雪、降雨量同比大幅度快速增長(zhǎng).而12月由于氣溫較低,地表濕冷且積雪覆蓋堆積給尾礦庫(kù)潰壩造成較大風(fēng)險(xiǎn);相比之下,夏季淳安縣雖然雨水豐沛,但由于地表溫度較高導(dǎo)致的水分蒸發(fā)迅速使得庫(kù)水位維持在安全水平,地下傾斜角偏移不大,對(duì)尾礦庫(kù)的安全狀況威脅較小.

    圖4 基于時(shí)間的尾礦庫(kù)風(fēng)險(xiǎn)評(píng)估結(jié)果Fig.4 Comprehensive evaluation of tailings security situation

    圖5 庫(kù)水位歷史狀況Fig.5 Reservoir water level safety history

    比較圖4,6可以發(fā)現(xiàn):地下傾斜角上中下三個(gè)空間位置所反映的尾礦庫(kù)安全狀況在特定時(shí)刻對(duì)尾礦庫(kù)總體安全狀況的影響.將傾斜角不同位置的安全狀況分值與尾礦庫(kù)總體分值對(duì)比發(fā)現(xiàn)由于各種不穩(wěn)定因素而產(chǎn)生形變的尾礦庫(kù)區(qū)域,便于尾礦庫(kù)管理人員重點(diǎn)監(jiān)測(cè),將潰壩事故防患于未然.

    圖6 地下傾斜角反映的尾礦庫(kù)風(fēng)險(xiǎn)評(píng)估結(jié)果Fig.6 Underground tilt angle reflects the results of the risk assessment

    4結(jié)論

    基于所建立的模型,在原有建立判斷矩陣基礎(chǔ)上,為代表指標(biāo)不同空間信息的各傳感器也建立判斷矩陣,運(yùn)用方根法驗(yàn)證矩陣一致性并利用馬氏距離保證判斷矩陣的嚴(yán)謹(jǐn)性與科學(xué)性,由下到上逐層評(píng)價(jià)直至完成對(duì)尾礦庫(kù)安全現(xiàn)狀的判斷,并結(jié)合空間信息分析尾礦庫(kù)潰壩隱患并對(duì)未來(lái)風(fēng)險(xiǎn)作出實(shí)時(shí)預(yù)測(cè),最后以淳安縣尾礦庫(kù)為例證明了方法的有效性.

    參考文獻(xiàn):

    [1]國(guó)家安全監(jiān)管總局等七部門(mén).全國(guó)尾礦庫(kù)綜合治理行動(dòng)2014年工作總結(jié)和2015年重點(diǎn)工作安排的通知[EB/OL].[2015-07-17].http://www.chinasafety.gov.cn/newpage/Contents/Channel_6288/2015/0717/254838/content_254838.html.

    [2]RICOM,BENITOG,SALGUEIROAR.Reportedtailingsdamfailures:areviewoftheEuropeanincidentsintheworldwidecontext[J].Journalofhazardousmaterials,2008,152(2):846-852.

    [3]CHAKRABORTYD,CHOUDHURYD.Investigationofthebehavioroftailingsearthendamunderseismicconditions[J].Americanjournalofengineeringandappliedsciences,2009,2(3):559-564.

    [4]PETERMB,MAHMOODSK.Seismicstabilityofimpoundments[C]//VancouverGeotechnicalSociety.Proceedingsofthe17thVGSonGeotechnicalEngineeringforGeoenvironmentalApplications.BritishColumbia:Bi-TechPublishersLtd,2003:77-84.

    [5]李忠奎,廖國(guó)禮.尾礦庫(kù)潰壩監(jiān)測(cè)預(yù)警系統(tǒng)設(shè)計(jì)研究[J].有色金屬(礦山部分),2008,60(6):33-35.

    [6]謝旭陽(yáng),江田漢,王云海,等.基于支持向量機(jī)的尾礦庫(kù)災(zāi)害區(qū)域預(yù)警[J].中國(guó)安全生產(chǎn)科學(xué)技術(shù),2008,4(4):17-21.

    [7]陳國(guó)定,姚景新.基于堆棧式自編碼器的尾礦庫(kù)安全評(píng)價(jià)[J].浙江工業(yè)大學(xué)學(xué)報(bào),2015,43(3):326-331.

    [8]趙志軍,閻高偉,謝克明.尾礦大壩安全監(jiān)測(cè)系統(tǒng)研究[J].現(xiàn)代電子技術(shù),2010(5):197-199.

    [9]李全明,陳仙,王云海,等.基于模糊理論的尾礦庫(kù)潰壩風(fēng)險(xiǎn)評(píng)價(jià)模型研究[J].中國(guó)安全生產(chǎn)科學(xué)技術(shù),2008,4(6):57-61.

    [10]羅飛飛,李慶軍.基于灰色模糊理論的尾礦庫(kù)安全評(píng)價(jià)方法研究[J].工業(yè)安全與環(huán)保,2009,35(8):46-48.

    [11]鄭銳,楊振宏,潘成林.基于熵技術(shù)的尾礦庫(kù)安全模糊綜合評(píng)價(jià)體系研究[J].中國(guó)安全生產(chǎn)科學(xué)技術(shù),2011,7(6):107-111.

    [12]徐曉敏.層次分析法的運(yùn)用[J].統(tǒng)計(jì)與決策,2008(1):156-158.

    [13]ZHANGYue,GAOXu,LIFeng,etal.StudyonsafetyevaluationfortailingreservoirsbasedontheIAHP-fuzzycomprehensiveevaluation[C]//IEEE.IEEEconferenceanthology.Washington:IEEEComputerSociety,2013:1-6.

    [14]劉俊萍.基于AHP和GIS的海塘工程安全評(píng)價(jià)研究[J].浙江工業(yè)大學(xué)學(xué)報(bào),2014,42(4):383-387.

    [15]吳力平,馮楊,李剛.基于層次分析法的施工現(xiàn)場(chǎng)平面布置方案評(píng)估[J].浙江工業(yè)大學(xué)學(xué)報(bào),2010,38(1):111-114.

    [16]戴建華,李軍,申文明,等.模糊判斷矩陣的滿(mǎn)意一致性檢驗(yàn)[J].系統(tǒng)工程與電子技術(shù),2006(8):1174-1177.

    [17]ZADEHLA.模糊集合、語(yǔ)言變量及模糊邏輯[M].陳國(guó)權(quán),譯.北京:科學(xué)出版社,1982:199.

    [18]叢卓紅,吳正鶯,鄭南翔,等.高速公路養(yǎng)護(hù)預(yù)算定額測(cè)定對(duì)象的優(yōu)選方法[J].重慶交通大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,30(4):823-827.

    [19]湯建國(guó),佘堃,祝峰.一種新的覆蓋粗糙模糊集模型[J].控制與決策,2012,27(11):1653-1662.

    (責(zé)任編輯:陳石平)

    Study on temporal and spatial risk prediction model of tailing reservoir

    YAN Shihang, FU Jiahong

    (College of Computer Science and Technology, Zhejiang University of Technology, Hangzhou 310023, China)

    Abstract:The risk prediction of tailings reservoir is one of the important means to strengthen the safety management level of tailings dam and improving the emergency response capability. Aimed at the tailings complex internal structure, the risk forecast index are not clearly defined, stability evaluation and warning method has not yet been scientific quantified, the tailings dam safety monitoring indicators are optimized, researching the change law of monitoring data, multi-level data analysis model was established to quantitative evaluate the risk situation of tailing pond and real-time predict of tailing pond safety development trend. Finally, the effectiveness of this method is verified through the example.

    Keywords:geologic hazard; risk assessment; tailings; spatio-temporal information

    收稿日期:2015-11-05

    作者簡(jiǎn)介:顏世航(1979—),男,黑龍江哈爾濱人,講師,研究方向?yàn)槲锫?lián)網(wǎng)、大數(shù)據(jù)分析和云計(jì)算平臺(tái)等,E-mail:shyan@zjut.edu.cn.

    中圖分類(lèi)號(hào):TP391

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1006-4303(2016)03-0254-06

    猜你喜歡
    潰壩尾礦庫(kù)矩陣
    尾礦庫(kù)空間信息提取與形變監(jiān)測(cè)應(yīng)用
    尾礦庫(kù)的環(huán)保防滲設(shè)計(jì)分析
    云南化工(2021年5期)2021-12-21 07:41:42
    筑牢尾礦庫(kù)安全防線(xiàn)
    徐家河尾礦庫(kù)潰壩分析
    初等行變換與初等列變換并用求逆矩陣
    潰壩涌浪及其對(duì)重力壩影響的數(shù)值模擬
    潰壩波對(duì)單橋墩作用水力特性研究
    基于改進(jìn)控制方程的土石壩潰壩洪水演進(jìn)數(shù)值模擬
    矩陣
    南都周刊(2015年4期)2015-09-10 07:22:44
    矩陣
    南都周刊(2015年3期)2015-09-10 07:22:44
    日韩不卡一区二区三区视频在线| 中文字幕免费在线视频6| 成人午夜精彩视频在线观看| 制服丝袜香蕉在线| eeuss影院久久| 97在线视频观看| 老司机影院成人| 熟妇人妻不卡中文字幕| 久久女婷五月综合色啪小说 | 欧美xxxx性猛交bbbb| 蜜桃久久精品国产亚洲av| 国产乱人偷精品视频| 国产69精品久久久久777片| 好男人在线观看高清免费视频| 天堂俺去俺来也www色官网| 亚洲av中文字字幕乱码综合| 久久精品人妻少妇| av在线天堂中文字幕| 老司机影院成人| 国产精品麻豆人妻色哟哟久久| av福利片在线观看| 精品一区二区免费观看| 国产精品久久久久久av不卡| 97在线人人人人妻| 国产精品久久久久久久电影| 99久久中文字幕三级久久日本| 欧美bdsm另类| 欧美 日韩 精品 国产| 国产精品国产三级国产专区5o| 黄色怎么调成土黄色| 国产成人精品婷婷| 国产一区亚洲一区在线观看| 国产精品人妻久久久影院| 黑人高潮一二区| 久久精品综合一区二区三区| 国产黄片视频在线免费观看| 好男人视频免费观看在线| 一级毛片我不卡| 久久久久九九精品影院| 五月天丁香电影| 丝袜脚勾引网站| 亚洲成人一二三区av| 久久精品国产鲁丝片午夜精品| 中文字幕久久专区| 你懂的网址亚洲精品在线观看| 草草在线视频免费看| 丝袜美腿在线中文| 国产综合懂色| 亚洲av日韩在线播放| 亚洲精品乱久久久久久| 国产 精品1| 成人欧美大片| 国产精品国产三级专区第一集| 久久精品综合一区二区三区| 国产成人aa在线观看| 国产永久视频网站| 2018国产大陆天天弄谢| 少妇熟女欧美另类| 看免费成人av毛片| 国产成人91sexporn| 亚洲婷婷狠狠爱综合网| 国产日韩欧美亚洲二区| 午夜福利在线观看免费完整高清在| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产成人一精品久久久| 国产高清不卡午夜福利| 在线观看三级黄色| 国产精品嫩草影院av在线观看| 国产一区有黄有色的免费视频| 一级毛片黄色毛片免费观看视频| 高清av免费在线| 男人爽女人下面视频在线观看| 国产淫语在线视频| 久久99精品国语久久久| 69av精品久久久久久| 99热全是精品| 成人欧美大片| 在线播放无遮挡| 国产精品精品国产色婷婷| 亚洲精品国产成人久久av| 国产高清有码在线观看视频| 亚洲高清免费不卡视频| 色婷婷久久久亚洲欧美| 在线 av 中文字幕| 国产精品久久久久久av不卡| 性色avwww在线观看| 下体分泌物呈黄色| 蜜桃亚洲精品一区二区三区| 成年av动漫网址| av女优亚洲男人天堂| 色视频www国产| 九九爱精品视频在线观看| 午夜福利视频1000在线观看| 人妻少妇偷人精品九色| 春色校园在线视频观看| 在线免费十八禁| 久久久久久久久久久丰满| 老司机影院毛片| 国产精品久久久久久精品电影| 亚洲欧美清纯卡通| 少妇的逼水好多| 亚洲精品久久久久久婷婷小说| 狂野欧美激情性xxxx在线观看| 男插女下体视频免费在线播放| 女的被弄到高潮叫床怎么办| 日韩欧美一区视频在线观看 | 亚洲av福利一区| av在线播放精品| 国产淫片久久久久久久久| tube8黄色片| 国产欧美亚洲国产| 看非洲黑人一级黄片| 最新中文字幕久久久久| 国产男女超爽视频在线观看| av在线观看视频网站免费| 国产国拍精品亚洲av在线观看| 久久午夜福利片| 少妇 在线观看| 日韩一区二区视频免费看| 三级国产精品欧美在线观看| 免费看a级黄色片| 亚洲av免费在线观看| 少妇裸体淫交视频免费看高清| 国产 精品1| 免费少妇av软件| 成人高潮视频无遮挡免费网站| 男的添女的下面高潮视频| 少妇被粗大猛烈的视频| 国产高清有码在线观看视频| 女的被弄到高潮叫床怎么办| 少妇高潮的动态图| 熟女电影av网| 深夜a级毛片| 18禁在线播放成人免费| av女优亚洲男人天堂| 欧美丝袜亚洲另类| 欧美成人精品欧美一级黄| 精品视频人人做人人爽| 午夜福利网站1000一区二区三区| av又黄又爽大尺度在线免费看| 一边亲一边摸免费视频| 午夜免费鲁丝| 免费看日本二区| 欧美日本视频| 国产高清不卡午夜福利| 亚洲自偷自拍三级| 一级二级三级毛片免费看| 日日摸夜夜添夜夜爱| 在现免费观看毛片| 久热久热在线精品观看| 男女那种视频在线观看| kizo精华| 国产男女超爽视频在线观看| 国产黄频视频在线观看| 黑人高潮一二区| 最近中文字幕高清免费大全6| 国产精品嫩草影院av在线观看| 高清毛片免费看| 大香蕉久久网| 黄片wwwwww| 国产国拍精品亚洲av在线观看| 久久精品综合一区二区三区| 最近的中文字幕免费完整| av黄色大香蕉| 熟女av电影| 日本-黄色视频高清免费观看| 日韩伦理黄色片| 校园人妻丝袜中文字幕| 亚洲欧美日韩另类电影网站 | 日本猛色少妇xxxxx猛交久久| 免费大片18禁| 男人和女人高潮做爰伦理| 欧美 日韩 精品 国产| 小蜜桃在线观看免费完整版高清| 网址你懂的国产日韩在线| 免费观看在线日韩| 汤姆久久久久久久影院中文字幕| 日韩不卡一区二区三区视频在线| 国产高清国产精品国产三级 | 中文字幕亚洲精品专区| 久久热精品热| 一区二区三区免费毛片| 人妻一区二区av| 精品一区二区免费观看| 在线免费观看不下载黄p国产| 日韩国内少妇激情av| 18禁裸乳无遮挡动漫免费视频 | 欧美bdsm另类| 亚洲欧美日韩另类电影网站 | 欧美精品人与动牲交sv欧美| 欧美成人a在线观看| 国产色婷婷99| 国内揄拍国产精品人妻在线| 在线免费观看不下载黄p国产| 九草在线视频观看| av一本久久久久| 99热这里只有精品一区| 亚洲精品日韩在线中文字幕| 亚洲精品影视一区二区三区av| 日韩视频在线欧美| 日韩人妻高清精品专区| 免费大片18禁| 亚洲av中文字字幕乱码综合| 免费电影在线观看免费观看| 女的被弄到高潮叫床怎么办| 亚洲精品久久久久久婷婷小说| 好男人视频免费观看在线| 看黄色毛片网站| 亚洲国产av新网站| 成人鲁丝片一二三区免费| 国产av码专区亚洲av| 噜噜噜噜噜久久久久久91| 亚洲不卡免费看| 午夜福利高清视频| 亚洲国产最新在线播放| 亚洲欧美日韩东京热| 大片电影免费在线观看免费| 国产乱人偷精品视频| 亚洲精品久久久久久婷婷小说| 丰满少妇做爰视频| 成人高潮视频无遮挡免费网站| 大片电影免费在线观看免费| 中文在线观看免费www的网站| 婷婷色麻豆天堂久久| 久久午夜福利片| 亚洲精品一二三| 免费电影在线观看免费观看| 内地一区二区视频在线| 中文欧美无线码| 禁无遮挡网站| 网址你懂的国产日韩在线| 激情五月婷婷亚洲| 乱码一卡2卡4卡精品| 免费看a级黄色片| 高清欧美精品videossex| 内射极品少妇av片p| 久久人人爽人人爽人人片va| 亚洲精品国产av蜜桃| 日韩欧美 国产精品| 亚洲欧美日韩卡通动漫| 久久亚洲国产成人精品v| 亚洲av免费在线观看| 晚上一个人看的免费电影| 成人毛片a级毛片在线播放| 丝袜美腿在线中文| 搞女人的毛片| 舔av片在线| 在线观看一区二区三区激情| 99热全是精品| 日本黄大片高清| 国产综合精华液| 日本三级黄在线观看| 欧美激情国产日韩精品一区| 久久精品国产亚洲网站| 成人综合一区亚洲| 日韩一本色道免费dvd| 黑人高潮一二区| 国产亚洲精品久久久com| 丝瓜视频免费看黄片| 国产探花极品一区二区| freevideosex欧美| 一级毛片电影观看| 简卡轻食公司| videos熟女内射| 亚洲av二区三区四区| 日韩欧美精品v在线| 爱豆传媒免费全集在线观看| 婷婷色综合www| 亚洲精品一区蜜桃| 欧美变态另类bdsm刘玥| 国产极品天堂在线| 亚洲色图综合在线观看| 国产伦理片在线播放av一区| 综合色av麻豆| 成人午夜精彩视频在线观看| 少妇的逼好多水| 新久久久久国产一级毛片| 亚洲熟女精品中文字幕| 国产av国产精品国产| 男插女下体视频免费在线播放| 国产日韩欧美亚洲二区| 国产精品一区二区在线观看99| 欧美日韩精品成人综合77777| 超碰av人人做人人爽久久| 婷婷色综合www| 国产 一区 欧美 日韩| 九九在线视频观看精品| 亚洲av.av天堂| 欧美激情在线99| 熟妇人妻不卡中文字幕| 久久久久久久久久久免费av| 欧美日韩精品成人综合77777| 在线观看国产h片| 一级av片app| 人妻制服诱惑在线中文字幕| 国内揄拍国产精品人妻在线| 最近手机中文字幕大全| 亚洲成色77777| 亚洲精品久久午夜乱码| 内射极品少妇av片p| 亚洲精品国产色婷婷电影| 成人综合一区亚洲| 久久人人爽人人爽人人片va| 国产免费又黄又爽又色| 亚洲av福利一区| 又爽又黄a免费视频| 日本午夜av视频| a级毛色黄片| 亚洲精品aⅴ在线观看| 最近手机中文字幕大全| 国产精品久久久久久久电影| 日韩免费高清中文字幕av| 欧美激情久久久久久爽电影| 亚洲色图综合在线观看| 黄色一级大片看看| 欧美zozozo另类| 搞女人的毛片| 日本爱情动作片www.在线观看| 免费av观看视频| 精品久久久久久久久av| 联通29元200g的流量卡| 国产亚洲91精品色在线| 精品久久久噜噜| 国产精品国产三级国产av玫瑰| 亚洲精品乱码久久久v下载方式| 精华霜和精华液先用哪个| 国产69精品久久久久777片| 亚洲精品日韩在线中文字幕| 麻豆乱淫一区二区| freevideosex欧美| 麻豆乱淫一区二区| 国产精品福利在线免费观看| 超碰97精品在线观看| 免费高清在线观看视频在线观看| 麻豆乱淫一区二区| av播播在线观看一区| 久久精品国产亚洲网站| 国产精品成人在线| 狂野欧美白嫩少妇大欣赏| 1000部很黄的大片| av女优亚洲男人天堂| 亚洲欧美一区二区三区黑人 | 国产黄a三级三级三级人| 天堂网av新在线| 中国三级夫妇交换| 一级爰片在线观看| 777米奇影视久久| 久久久久久久亚洲中文字幕| 日韩三级伦理在线观看| 最近最新中文字幕大全电影3| 久久久成人免费电影| 丝瓜视频免费看黄片| 国内揄拍国产精品人妻在线| 精品亚洲乱码少妇综合久久| 成人国产av品久久久| 97超视频在线观看视频| 最近最新中文字幕大全电影3| 国产精品国产三级专区第一集| 日本av手机在线免费观看| 亚洲自偷自拍三级| 少妇人妻精品综合一区二区| 在线精品无人区一区二区三 | 久久久久国产网址| 日本午夜av视频| 黄色怎么调成土黄色| 久久久久久久大尺度免费视频| av在线亚洲专区| 22中文网久久字幕| 免费大片18禁| 欧美成人精品欧美一级黄| 插逼视频在线观看| 一级毛片电影观看| tube8黄色片| 亚洲怡红院男人天堂| 欧美成人午夜免费资源| 亚洲av男天堂| 亚洲在线观看片| 777米奇影视久久| 黄色欧美视频在线观看| 国产av不卡久久| 男女无遮挡免费网站观看| 中文精品一卡2卡3卡4更新| 亚洲av成人精品一区久久| 亚洲国产欧美在线一区| 一级毛片aaaaaa免费看小| 亚洲内射少妇av| 午夜免费观看性视频| 欧美另类一区| 亚洲无线观看免费| 国产男人的电影天堂91| 成人亚洲精品av一区二区| 夜夜看夜夜爽夜夜摸| 五月玫瑰六月丁香| 国产老妇伦熟女老妇高清| 成人欧美大片| 亚洲欧美精品专区久久| 亚洲欧美日韩东京热| 日韩欧美精品v在线| 熟女人妻精品中文字幕| 美女脱内裤让男人舔精品视频| 青青草视频在线视频观看| 在线亚洲精品国产二区图片欧美 | a级一级毛片免费在线观看| 麻豆久久精品国产亚洲av| 听说在线观看完整版免费高清| 蜜桃久久精品国产亚洲av| 欧美xxⅹ黑人| 久久久国产一区二区| 乱系列少妇在线播放| av在线亚洲专区| 精品国产三级普通话版| 最近2019中文字幕mv第一页| 成年女人看的毛片在线观看| 六月丁香七月| 麻豆成人av视频| 五月玫瑰六月丁香| 一级黄片播放器| 亚洲性久久影院| 久久99精品国语久久久| 一级a做视频免费观看| 秋霞在线观看毛片| 亚洲欧美成人精品一区二区| 中文欧美无线码| 亚洲最大成人手机在线| 国产免费一区二区三区四区乱码| 青春草视频在线免费观看| 一区二区av电影网| 99热全是精品| 天堂中文最新版在线下载 | 国产黄色视频一区二区在线观看| 白带黄色成豆腐渣| 国国产精品蜜臀av免费| 欧美成人午夜免费资源| 深夜a级毛片| 水蜜桃什么品种好| 国产有黄有色有爽视频| 国产精品女同一区二区软件| 国产精品麻豆人妻色哟哟久久| 国产精品偷伦视频观看了| 校园人妻丝袜中文字幕| 成年av动漫网址| 国产免费一级a男人的天堂| 青青草视频在线视频观看| 亚洲天堂国产精品一区在线| 亚洲精品视频女| 午夜老司机福利剧场| 成人欧美大片| 国产精品av视频在线免费观看| 色婷婷久久久亚洲欧美| 久久久久国产网址| 夫妻午夜视频| 亚洲人成网站在线播| 91久久精品电影网| 国产免费视频播放在线视频| 80岁老熟妇乱子伦牲交| 日韩中字成人| 热re99久久精品国产66热6| 亚洲精品亚洲一区二区| 国产精品久久久久久精品电影小说 | 国产亚洲av嫩草精品影院| 亚洲欧美一区二区三区黑人 | 久久国内精品自在自线图片| 日韩av不卡免费在线播放| 亚洲精品中文字幕在线视频 | 色吧在线观看| 性插视频无遮挡在线免费观看| 国产伦精品一区二区三区四那| 菩萨蛮人人尽说江南好唐韦庄| 国产精品一二三区在线看| 大码成人一级视频| 麻豆精品久久久久久蜜桃| 成人亚洲欧美一区二区av| 国产av国产精品国产| 国产探花极品一区二区| 嫩草影院新地址| 黄色一级大片看看| 大片免费播放器 马上看| 在线播放无遮挡| 九色成人免费人妻av| 乱码一卡2卡4卡精品| 久久久久精品久久久久真实原创| 香蕉精品网在线| 婷婷色综合大香蕉| 女的被弄到高潮叫床怎么办| 97人妻精品一区二区三区麻豆| 亚洲成人久久爱视频| 中文天堂在线官网| 老司机影院成人| 午夜福利在线观看免费完整高清在| 好男人视频免费观看在线| 久久久久久久久久久免费av| 69av精品久久久久久| 岛国毛片在线播放| 亚洲色图综合在线观看| 在线精品无人区一区二区三 | 久久精品久久精品一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | av在线老鸭窝| 我要看日韩黄色一级片| 国产男女超爽视频在线观看| 水蜜桃什么品种好| 男女下面进入的视频免费午夜| 国产成人a区在线观看| 成人漫画全彩无遮挡| 夫妻性生交免费视频一级片| 狂野欧美白嫩少妇大欣赏| 男女下面进入的视频免费午夜| 精品一区二区三卡| 啦啦啦中文免费视频观看日本| 一级毛片黄色毛片免费观看视频| 建设人人有责人人尽责人人享有的 | 亚洲av国产av综合av卡| 免费看日本二区| 哪个播放器可以免费观看大片| 日本黄色片子视频| 97在线人人人人妻| 91aial.com中文字幕在线观看| 最近手机中文字幕大全| 亚洲欧美日韩另类电影网站 | 高清在线视频一区二区三区| 国内揄拍国产精品人妻在线| 国产精品久久久久久久久免| 亚洲成人精品中文字幕电影| 99热6这里只有精品| 亚洲精品,欧美精品| 成人二区视频| 精品午夜福利在线看| 在线精品无人区一区二区三 | 国产淫语在线视频| 亚洲欧美日韩卡通动漫| 22中文网久久字幕| 亚洲久久久久久中文字幕| 国产精品三级大全| 国产白丝娇喘喷水9色精品| 丝袜脚勾引网站| 又大又黄又爽视频免费| 午夜免费观看性视频| 午夜福利视频精品| 91午夜精品亚洲一区二区三区| 国产久久久一区二区三区| 国产亚洲午夜精品一区二区久久 | 亚洲国产精品成人久久小说| 在线观看一区二区三区| 成人午夜精彩视频在线观看| 国产免费一区二区三区四区乱码| 午夜福利在线在线| 免费黄频网站在线观看国产| 黄色怎么调成土黄色| 三级国产精品欧美在线观看| 午夜精品国产一区二区电影 | 一级爰片在线观看| 69人妻影院| 成人毛片60女人毛片免费| 成人亚洲精品av一区二区| 国产精品三级大全| 老司机影院毛片| 99热这里只有是精品在线观看| 日韩大片免费观看网站| 99re6热这里在线精品视频| 色婷婷久久久亚洲欧美| 久久精品久久久久久噜噜老黄| 美女被艹到高潮喷水动态| 一二三四中文在线观看免费高清| 日韩一本色道免费dvd| 精品亚洲乱码少妇综合久久| 18禁在线播放成人免费| 亚洲av日韩在线播放| 亚洲精品国产av蜜桃| 久久久精品94久久精品| 色网站视频免费| 欧美3d第一页| 国产国拍精品亚洲av在线观看| 成年版毛片免费区| 亚洲欧美日韩另类电影网站 | 国产精品精品国产色婷婷| 成年女人看的毛片在线观看| 在线亚洲精品国产二区图片欧美 | .国产精品久久| 99热这里只有是精品在线观看| 七月丁香在线播放| 国产伦理片在线播放av一区| 久久久色成人| 最新中文字幕久久久久| 欧美高清成人免费视频www| 最近手机中文字幕大全| 五月伊人婷婷丁香| 另类亚洲欧美激情| 能在线免费看毛片的网站| 亚洲熟女精品中文字幕| 久久这里有精品视频免费| 能在线免费看毛片的网站| 国产 一区精品| 成人一区二区视频在线观看| 一区二区三区免费毛片| 国产真实伦视频高清在线观看| 夜夜爽夜夜爽视频| 亚洲精品久久午夜乱码| 插逼视频在线观看| 亚洲人与动物交配视频| 91狼人影院| 欧美日韩精品成人综合77777| 久久精品国产亚洲av天美| 亚洲av男天堂| 欧美变态另类bdsm刘玥| 女人十人毛片免费观看3o分钟| 天堂中文最新版在线下载 | 国产乱人偷精品视频| 国产一区二区亚洲精品在线观看| 日本熟妇午夜| 免费观看的影片在线观看| 蜜臀久久99精品久久宅男| 亚洲精品乱码久久久v下载方式| 爱豆传媒免费全集在线观看| 亚洲人成网站在线观看播放| 18禁在线播放成人免费|