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

        ?

        基于地貌單元與灰關(guān)聯(lián)分析的地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)

        2019-01-08 05:51:12李得立李小磊羅德江陳國(guó)東袁名康
        關(guān)鍵詞:風(fēng)險(xiǎn)性樂(lè)山市指標(biāo)值

        李得立,李小磊,羅德江,陳國(guó)東,袁名康

        (1.成都理工大學(xué)管理科學(xué)學(xué)院,成都 610059;2.中國(guó)地質(zhì)環(huán)境監(jiān)測(cè)院,北京 100081)

        中央全面深化改革領(lǐng)導(dǎo)小組第三十五次會(huì)議通過(guò)了《關(guān)于建立資源環(huán)境承載能力監(jiān)測(cè)預(yù)警長(zhǎng)效機(jī)制的若干意見(jiàn)》,強(qiáng)調(diào)“開展承載能力評(píng)價(jià),規(guī)范空間開發(fā)秩序,合理控制開發(fā)強(qiáng)度,促進(jìn)人口、經(jīng)濟(jì)、資源環(huán)境的空間均衡”[1]。地質(zhì)災(zāi)害評(píng)價(jià)作為地質(zhì)環(huán)境承載能力評(píng)價(jià)的主要內(nèi)容之一,其成果對(duì)于區(qū)域空間規(guī)劃的編制、社會(huì)經(jīng)濟(jì)空間布局的優(yōu)化具有重要意義。

        目前地質(zhì)災(zāi)害評(píng)價(jià)的研究主要有3個(gè)方面:一是評(píng)價(jià)單元的劃分,主要有柵格單元、地形地貌單元和行政單元3種[2-5],前者可利用GIS實(shí)現(xiàn)單元的快速剖分,但與地形地貌及地質(zhì)環(huán)境條件缺乏有機(jī)聯(lián)系[5],地形地貌單元?jiǎng)t可有效克服這一不足,行政單元?jiǎng)t利于根據(jù)評(píng)價(jià)結(jié)果制訂相關(guān)措施,便于管理;二是評(píng)價(jià)指標(biāo)體系的構(gòu)建,不同的學(xué)者其評(píng)價(jià)目的不同,所構(gòu)建的評(píng)價(jià)指標(biāo)體系也不一致,但大都包含地層巖性、地質(zhì)構(gòu)造、地形地貌等指標(biāo),部分研究還將人類工程活動(dòng)、植被覆蓋率等納入指標(biāo)體系中[6-9];三是評(píng)價(jià)方法的選取及新方法的引入,應(yīng)用較多的方法有層次分析法[2,6]、線性回歸模型[10-11]、模糊綜合評(píng)價(jià)模型[12]、信息量模型[13-14]等,這些方法模型適用性與特點(diǎn)各不相同,前兩種方法常用于確定地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)指標(biāo)體系的權(quán)重,后兩種方法則根據(jù)所構(gòu)建的評(píng)價(jià)指標(biāo)體系,提取指標(biāo)值,利用模型開展定量評(píng)價(jià)。

        影響地質(zhì)災(zāi)害風(fēng)險(xiǎn)性的因素?cái)?shù)量多且復(fù)雜,在進(jìn)行評(píng)價(jià)時(shí)要考慮所有的因素是不現(xiàn)實(shí)的,因此地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)指標(biāo)體系可以看成是一個(gè)“信息不完全”的系統(tǒng)?!靶畔⒉煌耆毕到y(tǒng)的解決方法主要有模糊數(shù)學(xué)和灰色系統(tǒng)[15]。利用模糊數(shù)學(xué)理論對(duì)地質(zhì)災(zāi)害風(fēng)險(xiǎn)性進(jìn)行評(píng)價(jià)時(shí),關(guān)鍵問(wèn)題不在于評(píng)價(jià)模型的構(gòu)建,而在于各評(píng)價(jià)指標(biāo)分級(jí)閾值的確定。各地區(qū)地質(zhì)構(gòu)造、地形地貌等千差萬(wàn)別,要制定統(tǒng)一的適用于各地區(qū)的閾值是困難和不現(xiàn)實(shí)的。灰關(guān)聯(lián)分析法則無(wú)須為各評(píng)價(jià)指標(biāo)確定分級(jí)標(biāo)準(zhǔn),適用于信息不完全的復(fù)雜系統(tǒng),其評(píng)價(jià)結(jié)果具有客觀性[16]。本文針對(duì)研究區(qū)崩塌、滑坡、泥石流等地質(zhì)災(zāi)害的特性,構(gòu)建地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)的灰關(guān)聯(lián)分析模型并開展評(píng)價(jià)工作,為地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)引入新的方法。

        1 評(píng)價(jià)指標(biāo)與評(píng)價(jià)方法

        1.1 評(píng)價(jià)指標(biāo)體系的構(gòu)建與權(quán)值計(jì)算

        地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)選取地形起伏度C1(m)、地貌類型C2、工程地質(zhì)巖組C3、活動(dòng)斷裂距離C4 (m)、巖溶發(fā)育程度C5、地震動(dòng)參數(shù)C6、年平均降雨量C7 (mm)、礦山分布密度C8(個(gè)/100km2)、土地利用程度C9、地震動(dòng)峰值加速度C10(g)等10個(gè)指標(biāo)[17]。在這些指標(biāo)中,工程地質(zhì)巖組、巖溶發(fā)育程度、地貌類型、土地利用程度等為定性指標(biāo),其定量化見(jiàn)表1。

        表1 定性指標(biāo)的分級(jí)標(biāo)準(zhǔn)

        利用熵值法計(jì)算權(quán)值:設(shè)有m個(gè)地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)與n個(gè)評(píng)價(jià)單元所構(gòu)成的數(shù)據(jù)集為S={x1,x2,…,xn},xi∈Rm,i=1,2,n。

        第一步:采用極差法對(duì)原始數(shù)據(jù)的預(yù)處理(式1),得到極差化處理后的值yij。

        (1)

        第二步:利用比重法(式2)計(jì)算評(píng)價(jià)值pij。

        (2)

        第三步:計(jì)算第j個(gè)地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)指標(biāo)的輸出熵E和差異性系數(shù)G。

        (3)

        Gj=1-Ejj=1,2,…,n

        (4)

        第四步:計(jì)算各個(gè)地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)指標(biāo)的權(quán)重w=(w1,w2,…,wn)′(式5)。

        (5)

        1.2 評(píng)價(jià)方法

        灰關(guān)聯(lián)分析法是灰色理論[18]應(yīng)用最廣泛的方法之一,它通過(guò)計(jì)算灰關(guān)聯(lián)度來(lái)實(shí)現(xiàn)對(duì)地質(zhì)災(zāi)害風(fēng)險(xiǎn)性的評(píng)價(jià),以參與評(píng)價(jià)的各地形地貌單元的指標(biāo)最優(yōu)值構(gòu)成一個(gè)虛擬的“風(fēng)險(xiǎn)性最小”的評(píng)價(jià)單元,以其各指標(biāo)值作為參考數(shù)列,參與評(píng)價(jià)的各地形地貌單元的指標(biāo)值作為比較數(shù)列,計(jì)算各參與評(píng)價(jià)的各地形地貌單元與“風(fēng)險(xiǎn)性最小”評(píng)價(jià)單元之間的關(guān)聯(lián)度,并根據(jù)關(guān)聯(lián)度的大小得出評(píng)價(jià)結(jié)論。具體做法是[16]:

        (1) 根據(jù)評(píng)價(jià)目的,并結(jié)合相關(guān)研究成果構(gòu)建地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)指標(biāo)體系,記為xi(i=1,2,…,m)。

        (2) 提取地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)指標(biāo)體系中各評(píng)價(jià)指標(biāo)值,并以各指標(biāo)的“最優(yōu)值”構(gòu)成參考數(shù)列x*={x1*,x2*,…,xm*},xj*=maxxi(i=1,2,…,m)和比較數(shù)列xk(k=1,2,…,n)。

        (3) 對(duì)原始指標(biāo)值作規(guī)范化處理,將各指標(biāo)值初值化為0~1之間的數(shù),以消除不同量綱對(duì)指標(biāo)值的影響(式1)。記參考數(shù)列和比較數(shù)列構(gòu)造的矩陣G,規(guī)范化后的矩陣Y。

        (4) 確定分辨系數(shù)0<ρ<1,計(jì)算各比較數(shù)列與參考數(shù)列中各點(diǎn)與點(diǎn)之間的關(guān)聯(lián)系數(shù)rki(式6)再以全部關(guān)聯(lián)系數(shù)構(gòu)成關(guān)聯(lián)系數(shù)矩陣R(式7):

        (6)

        (7)

        (5) 利用熵值法計(jì)算出權(quán)重w。

        (6) 由關(guān)聯(lián)系數(shù)矩陣R和權(quán)值向量w求評(píng)價(jià)結(jié)果Z(式8),再根據(jù)Z每一個(gè)分量值的大小對(duì)各評(píng)價(jià)單元地質(zhì)災(zāi)害風(fēng)險(xiǎn)性做出排序和等級(jí)劃分。

        Z=Rw

        (8)

        2 評(píng)價(jià)方法、結(jié)果與分析

        2.1 研究區(qū)概況

        2.1.1 自然地理

        樂(lè)山市位于四川中南部、四川盆地西南部、成都平原至川西南山地過(guò)渡帶。北連眉山市,東鄰自貢市,南接宜賓市和涼山彝族自治州,西界雅安市。介于東經(jīng)102°55′~104°00′、北緯28°25′~29°55′之間。南北長(zhǎng)165 km,東西寬90~120 km,幅員面積12 723 km2。市境地勢(shì)呈西南高、東北低的傾斜狀,地貌以山地為主,丘陵次之,平原較少,分別占全市幅員面積的67%、21%、12%。海拔高程306~4 288 m。氣候?qū)偎疅峤M合優(yōu)越的中亞熱帶季風(fēng)濕潤(rùn)氣候,為我國(guó)內(nèi)地四川盆地西部特殊的“海洋性氣候島”的一部分。地層從前震旦系峨邊群到第四系全新統(tǒng)都有賦存。

        2.1.2 地質(zhì)災(zāi)害現(xiàn)狀

        受自然和人類工程活動(dòng)的綜合作用,樂(lè)山市境內(nèi)地質(zhì)災(zāi)害種類發(fā)育較齊全、數(shù)量多,但分布不均,且具突發(fā)性及不可預(yù)見(jiàn)性。地質(zhì)災(zāi)害種類以滑坡、崩塌(危巖)為主,其次為潛在不穩(wěn)定斜坡、泥石流、地面塌陷(地面裂縫)。其中滑坡約占50%,崩塌約占30%,潛在不穩(wěn)定斜坡、泥石流、地面塌陷、地裂縫等地質(zhì)災(zāi)害約占20%。

        2.2 評(píng)價(jià)單元

        樂(lè)山市全境根據(jù)其地貌成因形態(tài)可分為三大類八亞類,即侵蝕構(gòu)造地形(高山、中山、低山)、構(gòu)造剝蝕地形(包括深切割丘陵、中切割丘陵、淺切割丘陵、)、侵蝕堆積地形(包括沖積扇平原、河流及階地等)。各類地貌亞類及單元(34個(gè))分布面積見(jiàn)表2。

        (1) 高山區(qū)(編號(hào)1-1、1-2)

        主要分布于樂(lè)山市西側(cè)邊緣及西南角,分為2塊,面積444.19 km2。該地域地勢(shì)起伏大,多為山脊、山峰地帶,地理上為樂(lè)山市與雅安市、涼山州的分界地帶。山坡陡峭,一般大于50°,海拔在3 500~4 288 m之間,相對(duì)高差大于1 000 m。 行政隸屬金口河區(qū)、峨邊縣、馬邊縣。

        (2) 中山區(qū)(編號(hào)2-1、2-2、2-3)

        主要分布于樂(lè)山市西南部的大渡河南北兩側(cè)的峨眉山市、金口河區(qū)、沙灣區(qū)西部、峨邊縣、馬邊縣及沐川縣西南邊緣地帶,分為3塊,分布面積大,共計(jì)5 587.23 km2。該地域地勢(shì)起伏很大,河谷切割深,峽谷多為“V”形谷,谷坡陡峭,坡度40°~80°,海拔在1 000~3 500 m之間,相對(duì)高差大,1 000~2 500 m。

        (3) 低山區(qū)(編號(hào)3-1、3-2、3-3、3-4、3-5、3-6)

        主要分布于樂(lè)山市北部、東南部,可分為6塊,分布面積2 535.8 km2。該地域低山突起,地勢(shì)起伏較大,切割較深,可見(jiàn)“V”形谷或不對(duì)稱的“U”形谷,谷坡較陡,一般30°~45°,多為深溝窄谷低山地貌。海拔一般在520~950 m之間,相對(duì)高差大于200 m。

        (4) 深切割丘陵(編號(hào)4-1、4-2、4-3、4-4、4-5、4-6、4-7、4-8、4-9)

        深切割丘陵區(qū)分布于研究區(qū)東北部及東部地區(qū),且分為9塊,分布面積1 222.21 km2。海拔340~560 m,平均海拔415 m,相對(duì)高差100~160 m,谷寬小于100 m,丘陵多呈單斜狀,分水嶺單薄,橫向溝谷發(fā)育,呈“U”形或“V”形,以深丘窄谷地貌為主。

        (5) 中切割丘陵(編號(hào)5-1、5-2、5-3、5-4、5-5)

        分布在研究區(qū)東北部,分割為5塊,面積1 228.93 km2。多為圓頂丘,“U”型谷,谷底較開闊,丘坡較陡。高程400~500 m,切割深度50~60 m。

        (6) 淺切割丘陵(編號(hào)6-1、6-2、6-3、6-4、6-5)

        分布在研究區(qū)東北部,分割為5塊,面積648.01 km2。以臺(tái)狀和趾狀淺丘為主,具有圓丘、塔狀丘、脊?fàn)钋鸬榷喾N形態(tài),丘頂圓緩常呈條帶形或串珠狀,溝谷開闊、平緩或者以殘丘零星點(diǎn)綴,形若平原。海拔在350~442 m之間,平均海拔383 m,相對(duì)高差30~50 m,谷寬大于100 m,以緩丘谷帶小壩地貌相結(jié)合。

        (7) 沖積平原(壩)區(qū)(編號(hào)7-1、7-2、7-3)

        主要分布在研究區(qū)東部的青衣江兩側(cè),涉及峨眉山市和夾江縣,另外在馬邊縣縣城一帶見(jiàn)一平壩。3塊平原(壩)面積494.98 km2。

        (8) 河流、階地、河漫(心)灘區(qū)(編號(hào)8-1)

        主要為岷江、大渡河、青衣江及沿岸,面積為561.66 km2。大部由河水占據(jù),河漫灘及心灘由現(xiàn)代沖洪積物(砂、卵石、礫石)構(gòu)成。

        樂(lè)山市全域可根據(jù)地貌成因形態(tài)分為三大類八亞類,即侵蝕構(gòu)造地形(高山、中山、低山)、構(gòu)造剝蝕地形(包括深切割丘陵、中切割丘陵、淺切割丘陵)、侵蝕堆積地形(包括沖積扇平原、河流及階地等)。各類地貌亞類及單元(34個(gè))分布面積見(jiàn)表2。

        表2 樂(lè)山市地貌類型及地質(zhì)災(zāi)害(隱患)點(diǎn)分布表

        2.3 指標(biāo)值的提取與權(quán)值的確定

        根據(jù)構(gòu)建的地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)指標(biāo)體系,分別提取其指標(biāo)值,對(duì)其中的定性指標(biāo)按表1的規(guī)則進(jìn)行定量化處理。因評(píng)價(jià)指標(biāo)體系中各指標(biāo)的量綱、數(shù)量級(jí)是不同的,在進(jìn)行綜合評(píng)價(jià)之前采用公式(1)對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理,處理后的部分指標(biāo)項(xiàng)如表3所示。

        表3 地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)指標(biāo)數(shù)據(jù)

        地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)的指標(biāo)體系中各指標(biāo)重要性是不同的,合理確定權(quán)重是提高評(píng)價(jià)結(jié)果可靠性的關(guān)鍵點(diǎn)之一。利用熵值法計(jì)算出了各指標(biāo)的權(quán)重w=( 0.103 1,0.104 5,0.105 4,0.103 6,0.104 6,0.105 7,0.105 8,0.105 8,0.056 9,0.104 6)。

        2.4 地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)

        在對(duì)樣本數(shù)據(jù)進(jìn)行處理和計(jì)算出地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)的指標(biāo)體系中各指標(biāo)權(quán)值后,利用構(gòu)建的灰關(guān)聯(lián)分析方法對(duì)樂(lè)山市所劃分出單元進(jìn)行評(píng)價(jià)。首先確定參考數(shù)列,因構(gòu)建的指標(biāo)皆為負(fù)向指標(biāo),故取表3中各列最小值作為最優(yōu)值,x*={25,25,45,25,25,45,65,65,0,25}。利用公式(6)計(jì)算關(guān)聯(lián)系數(shù)并進(jìn)行灰關(guān)聯(lián)評(píng)價(jià)(式8),評(píng)價(jià)的結(jié)果及其各評(píng)價(jià)單元地質(zhì)災(zāi)害風(fēng)險(xiǎn)性排序如表4所示。從表4的評(píng)價(jià)結(jié)果看,灰關(guān)聯(lián)度在0.616 5~0.935 5之間,可劃分為高(III級(jí):≤0.75)、中(II級(jí):0.75

        根據(jù)表4及圖1,高風(fēng)險(xiǎn)性(Ⅲ級(jí))區(qū)主要涉及中山區(qū)(2-1、2-2、2-3單元)、低山區(qū)(3-2、3-3單元)、深切丘陵區(qū)(4-1單元)、淺切丘陵區(qū)(6-1單元)等7個(gè)地貌單元;低風(fēng)險(xiǎn)性(Ⅰ級(jí))區(qū)涉及河谷區(qū)(7-1、7-2、7-3、8-1單元)等4個(gè)地貌單元;其余23個(gè)地貌單元均為中風(fēng)險(xiǎn)性(Ⅱ級(jí))。

        圖1 地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)結(jié)果及分級(jí)圖

        地貌編號(hào)評(píng)價(jià)結(jié)果排序地貌編號(hào)評(píng)價(jià)結(jié)果排序1-10.790684-70.86 11-20.790 674-80.757 922-10.729 7124-90.842 6272-20.710 845-10.799 5212-30.737 335-20.775233-10.774 455-30.803 3113-20.667 8265-40.834 4103-30.616 5165-50.842 8303-40.762 4136-10.743 2243-50.822 6196-20.792 3203-60.816 696-30.788 8144-10.679 4296-40.772 9254-20.756 766-50.825 3184-30.842 8157-10.869 2324-40.775227-20.865314-50.754 4177-30.909334-60.775 5288-10.935 534

        3 結(jié)論與討論

        (1) 針對(duì)地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)單元?jiǎng)澐謫?wèn)題,在對(duì)評(píng)價(jià)區(qū)地貌成因形態(tài)進(jìn)行分析的基礎(chǔ)上,劃分出侵蝕構(gòu)造地形、構(gòu)造剝蝕地形、侵蝕堆積地形等三大類,高山、中山、低山、深切割丘陵、中切割丘陵、淺切割丘陵、沖積扇平原、河流及階地等八亞類共計(jì)34個(gè)單元。這種劃分方法既方便提取各評(píng)價(jià)指標(biāo)值,也利于有針對(duì)性的提出防治措施。

        (2) 文中的評(píng)價(jià)指標(biāo)權(quán)值計(jì)算方法以及灰關(guān)聯(lián)分析法在樂(lè)山市資源環(huán)境承載力評(píng)價(jià)中進(jìn)行了實(shí)際應(yīng)用,并根據(jù)評(píng)價(jià)結(jié)果進(jìn)行了地質(zhì)災(zāi)害風(fēng)險(xiǎn)性的等級(jí)劃分。經(jīng)驗(yàn)證,使用該方法所得到的分區(qū)結(jié)果與實(shí)地調(diào)查結(jié)果吻合較好。

        (3) 文中的地質(zhì)災(zāi)害主要是指崩塌、滑坡、泥石流等。事實(shí)上由于地質(zhì)災(zāi)害的多樣性以及影響因素的復(fù)雜性,使得要建立科學(xué)、實(shí)用的指標(biāo)體系還需要深入的研究;其次要加強(qiáng)方法的理論分析和適用性分析外,還應(yīng)加強(qiáng)多種方法評(píng)價(jià)結(jié)果不一致性問(wèn)題的研究,以提高地質(zhì)災(zāi)害風(fēng)險(xiǎn)性評(píng)價(jià)的可靠性。

        猜你喜歡
        風(fēng)險(xiǎn)性樂(lè)山市指標(biāo)值
        樂(lè)山市在職男護(hù)士健康素養(yǎng)水平調(diào)查
        樂(lè)山市五通橋區(qū)“雙創(chuàng)”五通橋在行動(dòng)
        樂(lè)山市城市管理行政執(zhí)法局
        關(guān)于打造樂(lè)山市特色水上旅游的調(diào)查與思考
        淺談食品中大腸菌群檢測(cè)方法以及指標(biāo)值的對(duì)應(yīng)關(guān)系
        每天超8小時(shí)睡眠 癡呆癥概率增40%
        東西南北(2017年9期)2017-05-24 22:03:10
        維修性定性要求評(píng)價(jià)指標(biāo)融合模型研究
        羅甸縣外來(lái)生物—飛機(jī)草風(fēng)險(xiǎn)性評(píng)估
        綠色科技(2016年24期)2017-03-30 23:48:52
        醫(yī)院藥學(xué)服務(wù)風(fēng)險(xiǎn)性探討
        1995年—2013年地方預(yù)算內(nèi)財(cái)力、中央返還及上解情況
        美女扒开内裤让我捅的视频| 日韩www视频| 国产精品爽爽v在线观看无码| 日本丰满人妻xxxxxhd| 国产午夜精品福利久久| 国产成人自拍视频在线观看网站| 国产三级不卡一区不卡二区在线| 国产亚洲自拍日本亚洲| 国产精品无码一区二区三区在| 国精品无码一区二区三区在线蜜臀 | 国产影片免费一级内射| 丰满少妇高潮惨叫久久久| 女人被爽到呻吟gif动态图视看| 亚洲乱码一区二区三区成人小说| 国产精品美女自在线观看| 亚洲精品国产一二三区| 成片免费观看视频大全| 日批视频免费在线观看| 久久2020精品免费网站| 精品无码久久久久久久久水蜜桃| 久久精品女人天堂av| 精品无码久久久久久久久粉色| 一本色道亚州综合久久精品| 97人妻人人揉人人躁九色| 国产精品va无码一区二区| 亚洲一区二区三区久久不卡| 国产不卡av一区二区三区| 亚洲日韩成人无码| 色偷偷av亚洲男人的天堂| 亚洲av无码av在线播放| 东京热无码人妻中文字幕| 亚洲天堂亚洲天堂亚洲色图| 国产av一区二区三区性入口| 激情偷乱人成视频在线观看| 日本高清一区二区三区水蜜桃| 一区二区三区午夜视频在线观看 | 日本av亚洲中文字幕| 少妇下蹲露大唇无遮挡| 午夜影视啪啪免费体验区入口| 亚洲国产成人av毛片大全| 国产主播一区二区三区蜜桃|