999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于組合賦權(quán)-二維云模型的鐵路線路方案優(yōu)選研究

2024-03-13 01:53:08趙長石
關(guān)鍵詞:評(píng)價(jià)模型

趙長石

(中國鐵路設(shè)計(jì)集團(tuán)有限公司線路站場樞紐設(shè)計(jì)研究院,天津 300308)

引言

作為路網(wǎng)骨架,線路走向優(yōu)劣情況關(guān)系到整體路網(wǎng)功能的發(fā)揮和線路自身的作用效果。由于鐵路方案評(píng)價(jià)是一個(gè)多因素相互作用的復(fù)雜系統(tǒng),各因素本身的特性和其他因素的相互影響,使得系統(tǒng)中的部分評(píng)價(jià)指標(biāo)具有不確定性,評(píng)價(jià)指標(biāo)值的確定具有隨機(jī)性[1-3]。因此,建立一種能夠較好處理系統(tǒng)中不確定性、隨機(jī)性和模糊性的評(píng)價(jià)模型,是使評(píng)估方法更接近工程實(shí)際情況的重要保障。

近年來,關(guān)于鐵路線路方案比選的研究已開展了大量工作,國內(nèi)外眾多學(xué)者運(yùn)用層次分析法[4]、模糊綜合評(píng)價(jià)[5]、灰色關(guān)聯(lián)分析[6-7]、聚類分析[8],以及模糊數(shù)學(xué)理論中包含的其他評(píng)價(jià)方法,反映評(píng)價(jià)線路的優(yōu)劣情況,雖上述文獻(xiàn)已建立了相對(duì)完善的線路優(yōu)選模型,但在決策過程中線路方案多從工程建設(shè)結(jié)果考慮,缺乏對(duì)指標(biāo)體系工后恢復(fù)情況等的多維性思考。同時(shí),評(píng)價(jià)指標(biāo)選取也缺乏對(duì)具體環(huán)境的針對(duì)性,忽略了不同地區(qū)影響線路走向因素的重點(diǎn)差異,導(dǎo)致評(píng)價(jià)結(jié)果的準(zhǔn)確性仍有待進(jìn)一步提高。多維云模型是由楊朝暉[9]提出的一種衍生正態(tài)云,其在決策過程中能從不同維度考慮指標(biāo)間的交互作用,使評(píng)價(jià)結(jié)果更貼合工程實(shí)際,其已在風(fēng)險(xiǎn)評(píng)估[10-11]、質(zhì)量安全[12]等領(lǐng)域得到廣泛應(yīng)用,但在線路方案評(píng)價(jià)中還鮮有研究。

鑒于此,本文將二維云模型引入線路方案評(píng)價(jià)優(yōu)選模型,以沈白鐵路永陵至通化段線路為研究背景,選取工程建設(shè)結(jié)果和工后恢復(fù)成本作為二維云模型的基礎(chǔ)數(shù)字變量,再結(jié)合組合賦權(quán)法提高主客觀權(quán)重的精度與可信度,通過MATLAB編程輸出不同線路方案的二維評(píng)價(jià)云圖,直觀反映線路方案的評(píng)價(jià)結(jié)果。以期為沈白鐵路永陵至通化段線路方案比選提供理論依據(jù),也對(duì)類似工程提供借鑒參考。

1 影響線路走向因素分析

1.1 線路概況

沈白鐵路位于中國東北地區(qū)遼寧省東部和吉林省東南部,是連接遼寧省沈陽市與吉林省延邊朝鮮族自治州的高速鐵路,線路全長430 km,設(shè)計(jì)時(shí)速350 km,途經(jīng)沈陽、撫順,通化、延邊等市、縣及自治州,是《中長期鐵路網(wǎng)規(guī)劃》中“八縱八橫”高速鐵路主通道的重要區(qū)域連接線,其建設(shè)對(duì)于完善東北區(qū)域快速客運(yùn)網(wǎng)絡(luò)布局有重要意義。

新建沈白鐵路永陵至通化段線路西起撫平市永陵縣,東至通化市通化縣西部,工程區(qū)域內(nèi)地形地貌以低山重丘為主,沿線存在多處經(jīng)濟(jì)據(jù)點(diǎn)和生態(tài)文化保護(hù)區(qū)。線路建設(shè)前期,通過對(duì)工程布設(shè)地區(qū)的生態(tài)環(huán)境現(xiàn)場調(diào)查及地質(zhì)勘測,在掌握其環(huán)境特征的情況下,研究提出經(jīng)新賓、通港方案,經(jīng)新賓、通化縣西方案和經(jīng)桓仁北、通化縣西方案三條線路走向方案,確定方案比較范圍為CK125+000~CK225+160。具體線路方案及沿線環(huán)境如圖1所示。

圖1 沈白鐵路永陵至通化段線路方案

經(jīng)新賓、通港方案(方案Ⅰ):線路自清永陵南側(cè)向東出發(fā),上跨撫通高速公路,于新賓縣南側(cè)設(shè)新賓站,再向東北延伸,在老米溝村北側(cè)設(shè)置橋梁穿越喇蛄河國家濕地公園合理利用區(qū),隨后線位繼續(xù)向東,經(jīng)鶴大高速公路、G303國道后至方案比較終點(diǎn)。

經(jīng)新賓、通化縣西方案(方案Ⅱ):線路自清永陵南側(cè)向東出發(fā),上跨撫通高速公路后,在新賓縣南側(cè)建設(shè)新賓站,再向東南繞行,沿線經(jīng)喇蛄河國家濕地公園和水源保護(hù)區(qū)北側(cè),并在距通化縣西北側(cè)約10 km處金斗鄉(xiāng)設(shè)通化縣西站,出站后線路向東北方向行進(jìn),上跨鶴大高速、下穿通化市一級(jí)路后至通化縣方案比較終點(diǎn)。

經(jīng)桓仁北、通化縣西方案(方案Ⅲ):線路自方案比較起點(diǎn)向東南行進(jìn),上跨永恒高速公路,經(jīng)高麗溝、火眼嶺進(jìn)入本溪市境內(nèi),在恒仁縣北部設(shè)恒仁北站,出站后折向東北沿鶴大高速公路行進(jìn),經(jīng)大廟溝、小倉庫溝上跨撫通高速公路、蝲蛄河國家濕地公園后在通化縣西北側(cè)設(shè)通化西站,出站后線路繼續(xù)向東北方向行進(jìn),至通化縣方案比較終點(diǎn)。

1.2 影響線路走向因素

線路方案優(yōu)選的基礎(chǔ)是構(gòu)建科學(xué)合理的線路方案評(píng)價(jià)指標(biāo)體系,工程中各線路方案均有其自身特點(diǎn),目前關(guān)于線路走向的優(yōu)劣評(píng)價(jià)并沒有統(tǒng)一規(guī)范標(biāo)準(zhǔn)約束,且受地理環(huán)境及施工技術(shù)的影響,不同工程環(huán)境下影響線路走向的因素也并不完全相同[13]。因此,評(píng)價(jià)指標(biāo)選取時(shí)應(yīng)針對(duì)具體工程環(huán)境進(jìn)行構(gòu)建。

根據(jù)沈白鐵路永陵至通化段線路實(shí)際情況可知,線路布設(shè)地區(qū)周圍環(huán)境復(fù)雜,人文及生態(tài)自然保護(hù)區(qū)較多,且地形環(huán)境以山區(qū)為主,鐵路在不斷建設(shè)時(shí)易對(duì)沿線地區(qū)的生態(tài)系統(tǒng)完整性產(chǎn)生破壞,引起水土流失、泥石流等惡性地質(zhì)災(zāi)害,對(duì)交通線路的正常運(yùn)營構(gòu)成了極大威脅。本文針對(duì)線路所處的工程環(huán)境,通過深入了解線路所處地理環(huán)境和施工特征,以技術(shù)可行和生態(tài)平衡兩方面為導(dǎo)向,并結(jié)合前人已有研究成果[14-16]與專家咨詢,從線路設(shè)計(jì)、自然人文景觀破壞、水土資源污染、施工難度及安全四個(gè)方面構(gòu)建線路走向的評(píng)價(jià)指標(biāo)層次模型,并對(duì)不同類型影響因素的建設(shè)結(jié)果進(jìn)行等級(jí)劃分,如表1所示。

表1 線路方案優(yōu)選指標(biāo)層次模型

2 組合賦權(quán)二維云評(píng)價(jià)模型

2.1 評(píng)價(jià)指標(biāo)體系組合賦權(quán)

2.1.1 IGAHP法主觀賦權(quán)

改進(jìn)的AHP法能夠有效量化主觀判斷的不確定性,克服評(píng)價(jià)思維與判斷矩陣的矛盾,提高權(quán)重計(jì)算的精確性和便捷性[17]。其具體步驟如下。

(1)

(3)計(jì)算各評(píng)價(jià)指標(biāo)權(quán)重wi(w1,w2,…,wm),即各評(píng)價(jià)指標(biāo)權(quán)重系數(shù)

(2)

2.1.2 熵權(quán)法客觀賦權(quán)

熵表示一個(gè)系統(tǒng)評(píng)判指標(biāo)的離散程度,屬于客觀賦權(quán)方法。評(píng)價(jià)體系中某個(gè)指標(biāo)的信息熵值越大,其對(duì)評(píng)價(jià)結(jié)果的影響程度就越弱,賦予權(quán)重就越小;反之,某指標(biāo)的信息熵值越小,其對(duì)評(píng)價(jià)結(jié)果的影響程度就越強(qiáng),賦予權(quán)重越大[18]。熵權(quán)法的具體計(jì)算方程如下。

(3)

式中,m為待評(píng)價(jià)項(xiàng)目數(shù);Qij為第j個(gè)指標(biāo)下第i個(gè)項(xiàng)目的評(píng)價(jià)值;n為評(píng)價(jià)指標(biāo)個(gè)數(shù);Hj和Wj分別為第j個(gè)指標(biāo)的熵值、熵權(quán)。

2.1.3 線性加權(quán)法組合賦權(quán)

線性加權(quán)法即按照各目標(biāo)的重要程度給予其對(duì)應(yīng)的權(quán)系數(shù),再對(duì)其進(jìn)行線性組合得到計(jì)算結(jié)果。線性賦權(quán)法的主要計(jì)算公式如下

W=αW′+βW″

(4)

式中,W′和W″分別為IGAHP法和熵權(quán)法計(jì)算的主、客觀權(quán)重向量;α、β為向量的加權(quán)系數(shù),且α、β∈(0,1)。

(5)

式中,n為評(píng)價(jià)指標(biāo)個(gè)數(shù);Pi(i=1,2,3,…,n)為主觀權(quán)重向量W′中各評(píng)價(jià)指標(biāo)由小到大排列后的權(quán)重值。

2.2 二維云模型

在一維云模型的基礎(chǔ)上,引入二維云模型概念來描述2個(gè)因素協(xié)同作用下的復(fù)雜概念。設(shè)F為服從正態(tài)分布的二維隨機(jī)函數(shù),其數(shù)字特征由兩組期望(Ex)、熵(En)和超熵(He)組成,反映2個(gè)維度的定性概念與定量數(shù)值進(jìn)行轉(zhuǎn)化時(shí)綜合評(píng)價(jià)不確定問題的優(yōu)劣情況[19]。

(6)

稱符合式(6)的云滴drop(xi,yi,ui)構(gòu)建的云模型為二維云模型,其中xi和yi為云滴坐標(biāo),ui為隸屬度。

2.2.1 評(píng)價(jià)云模型

鐵路修建時(shí)的工程建設(shè)結(jié)果和建成后恢復(fù)成本共同反映線路方案的優(yōu)劣情況,因此以評(píng)價(jià)指標(biāo)為度量,選取建設(shè)成本與建后恢復(fù)作為線路方案評(píng)價(jià)的二維基礎(chǔ)變量。通過邀請(qǐng)工程建設(shè)時(shí)的評(píng)估專家對(duì)所選底層指標(biāo)進(jìn)行打分分析,則每個(gè)評(píng)價(jià)指標(biāo)衡量的建設(shè)結(jié)果和恢復(fù)成本分別形成一個(gè)云滴,組成各指標(biāo)的結(jié)果云和成本云,二者統(tǒng)稱為二維綜合評(píng)價(jià)云。按照式(7)運(yùn)用MATLAB編程計(jì)算建設(shè)結(jié)果和恢復(fù)成本云的各自數(shù)字特征。

(7)

式中,q為樣本數(shù)量;xk為專家對(duì)樣本打分值;Ex為期望;En為熵,數(shù)學(xué)期望和熵反映樣本信息的隨機(jī)性;He為超熵,在評(píng)價(jià)云圖中反映云滴凝聚的松緊程度。

綜合評(píng)價(jià)云是由一級(jí)工程建設(shè)結(jié)果云和工后恢復(fù)成本云數(shù)字特征矩陣與對(duì)應(yīng)的組合賦權(quán)權(quán)重矩陣合成得到,同樣,一級(jí)評(píng)價(jià)云也由二級(jí)評(píng)價(jià)云與其對(duì)應(yīng)權(quán)重合成而來。綜合評(píng)價(jià)云和一級(jí)評(píng)價(jià)云數(shù)字特征的具體合成方程如下

(Ex′,En′,He′)

(8)

式中,Ex′,En′,He′為上一等級(jí)風(fēng)險(xiǎn)云的數(shù)字特征。

2.2.2 標(biāo)準(zhǔn)云模型

(9)

式中,k為常數(shù),其數(shù)值范圍在0.01~1之間[20],本文根據(jù)變量的模糊閾值,k值取0.01。

建設(shè)結(jié)果評(píng)價(jià)等級(jí)已按表1進(jìn)行劃分,恢復(fù)成本的評(píng)價(jià)等級(jí)在查閱資料和借鑒前人研究成果[21]后進(jìn)行劃分,標(biāo)準(zhǔn)云模型兩方面基礎(chǔ)變量的區(qū)間分值、數(shù)字特征相同,具體等級(jí)描述見表2。

表2 評(píng)價(jià)等級(jí)描述

2.2.3 云模型相近度

為精確分辨評(píng)價(jià)云圖的隸屬區(qū)間和等級(jí),引入相近度概念,通過云模型數(shù)字特征計(jì)算最貼近標(biāo)準(zhǔn)云,具體公式如下

(10)

2.3 正向云發(fā)生器

通過MATLAB編程運(yùn)用正向云發(fā)生器輸出二維評(píng)價(jià)云圖,其關(guān)鍵流程如圖2所示。

圖2 二維正向云發(fā)生器

3 分析討論

3.1 評(píng)價(jià)過程

為滿足工程實(shí)際情況需要,根據(jù)前文理論建立基于組合賦權(quán)二維云的線路評(píng)價(jià)模型。邀請(qǐng)6位參與該項(xiàng)目設(shè)計(jì)專家和2位選線領(lǐng)域的高校學(xué)者,按表1所劃分等級(jí)情況,分別對(duì)不同線路方案修建后的工程建設(shè)結(jié)果和工后恢復(fù)成本進(jìn)行評(píng)分,分?jǐn)?shù)記為x1和x2,以經(jīng)新賓、通港方案為例,詳細(xì)評(píng)價(jià)結(jié)果如表3所示。根據(jù)改進(jìn)AHP法計(jì)算影響因素主觀權(quán)重,式(3)熵權(quán)法計(jì)算影響因素客觀權(quán)重,式(4)、式(5)進(jìn)行組合賦權(quán)求得各級(jí)指標(biāo)最終權(quán)重,計(jì)算結(jié)果如表4所示。

表3 建設(shè)結(jié)果和修復(fù)成本評(píng)分值

表4 各級(jí)指標(biāo)加權(quán)后權(quán)重

3.2 結(jié)果分析

為判斷線路方案在不同方面的優(yōu)劣情況,選取一級(jí)評(píng)價(jià)指標(biāo)進(jìn)行對(duì)比,將二級(jí)指標(biāo)專家評(píng)分值代入式(7)計(jì)算,得到二級(jí)指標(biāo)建設(shè)結(jié)果云和恢復(fù)成本云的數(shù)字特征,應(yīng)用式(8)將二級(jí)評(píng)價(jià)云數(shù)字特征與其對(duì)應(yīng)權(quán)重矩陣進(jìn)行合成運(yùn)算,可得一級(jí)指標(biāo)評(píng)價(jià)云數(shù)字特征。通過MATLAB正向云發(fā)生器輸入不同線路方案風(fēng)險(xiǎn)云數(shù)字特征,得到不同線路方案一級(jí)指標(biāo)的評(píng)價(jià)云與標(biāo)準(zhǔn)云,如圖3所示。

圖3 不同方案一級(jí)指標(biāo)評(píng)價(jià)云對(duì)比(圖中綠、紅、藍(lán)分別代表方案Ⅰ、Ⅱ、Ⅲ級(jí)指標(biāo)的二維評(píng)價(jià)云)

由圖3可知,3個(gè)線路方案的優(yōu)缺點(diǎn)都十分鮮明,在線路設(shè)計(jì)方面,經(jīng)新賓、通港方案(方案Ⅰ)和經(jīng)新賓、通化縣西方案(方案Ⅱ)采用航空直線原則,線路均較為順直,其中方案Ⅰ線路正線里程最短,總長約91.37 km,橋隧比也相對(duì)其他線路方案最小(79.8%),表明該線路對(duì)地形環(huán)境具有較好的適應(yīng)性,方案I線路走向由永陵南側(cè)直接向西展線,途中穿越蝲蛄河濕地公園合理緩沖區(qū),雖對(duì)原有生態(tài)環(huán)境產(chǎn)生了一定改變,但并不破壞現(xiàn)有動(dòng)植物的正常生存,且后期修復(fù)成本屬于可接受范圍,而方案Ⅱ受設(shè)站位置影響,正線長度較方案Ⅰ超出約3 km,沿線工程區(qū)域內(nèi)的生態(tài)保護(hù)區(qū)較密集,后半段線路與蝲蛄河水資源保護(hù)區(qū)較近,工程建設(shè)時(shí)易對(duì)原有的自然景觀和水土環(huán)境產(chǎn)生破壞,后期恢復(fù)成本極大。經(jīng)桓仁北、通化縣西方案(方案Ⅲ)直接向西南繞行,穿越龍崗山省級(jí)自然保護(hù)區(qū)的實(shí)驗(yàn)區(qū),減小了鐵路建設(shè)對(duì)生態(tài)環(huán)境的影響,但線路里程長度較方案Ⅰ超出約30.6 km,橋隧長度占正線總長的81.2%,受山區(qū)復(fù)雜地形環(huán)境的制約,較差的線形與地形環(huán)境的適應(yīng)能力,極大提高了線路建設(shè)的工程成本和維護(hù)成本。通過上述說明解釋了不同方案的鐵路修建時(shí)直接產(chǎn)生的Z1、Z2、Z3評(píng)價(jià)云相差的原因。另外,在施工安全及難度方面,經(jīng)新賓、通港方案(方案Ⅰ)沿線工程地質(zhì)環(huán)境良好,壓覆礦產(chǎn)資源較少,且施工區(qū)距新賓、通化縣城較近,在材料獲取便捷性和交通運(yùn)輸條件上都相對(duì)高效,而另外兩線路方案由于線路的走向問題,施工材料的運(yùn)輸環(huán)境較差,建設(shè)過程將不可避免破壞現(xiàn)有植被環(huán)境,且經(jīng)桓仁北、通化縣西方案(方案Ⅲ)沿線壓覆礦產(chǎn)資源多,工程地質(zhì)屬性較差,鐵路建設(shè)的難度及后期修復(fù)成本也較高。

為進(jìn)一步確定不同線路方案的評(píng)價(jià)結(jié)果等級(jí),將綜合評(píng)價(jià)云和標(biāo)準(zhǔn)云數(shù)字特征輸入MATLAB正向云發(fā)生器,得到綜合評(píng)價(jià)云與標(biāo)準(zhǔn)云的對(duì)比圖,見圖4。

圖4 不同線路方案綜合評(píng)價(jià)云

圖5為線路方案的綜合評(píng)價(jià)云俯視圖,可判斷線路評(píng)價(jià)結(jié)果所處的標(biāo)準(zhǔn)云等級(jí)范圍。

圖5 不同線路方案綜合評(píng)價(jià)云俯視

由綜合評(píng)價(jià)云圖可以看出,方案Ⅰ評(píng)價(jià)等級(jí)介于良好和優(yōu)秀之間,更貼近良好;方案Ⅱ、方案Ⅲ評(píng)價(jià)等級(jí)介于一般和良好之間,更貼近一般;為更準(zhǔn)確地確定不同線路的評(píng)價(jià)等級(jí),應(yīng)用式(10)進(jìn)行計(jì)算綜合風(fēng)險(xiǎn)云與標(biāo)準(zhǔn)云之間的貼近度,結(jié)果表明方案Ⅰ的綜合評(píng)價(jià)等級(jí)為“良好”,方案Ⅱ、Ⅲ的綜合評(píng)價(jià)等級(jí)均為“一般”。因此,推薦經(jīng)新賓、通港方案(方案Ⅰ)為線路建設(shè)方案,這也與該項(xiàng)目專家所推薦方案相一致。

4 結(jié)論

(1)通過深入了解沈白鐵路通化段所處地理環(huán)境和施工特征,從技術(shù)可行性和生態(tài)平衡兩方面選取影響線路走向因素,并對(duì)其結(jié)果進(jìn)行等級(jí)劃分,建立了針對(duì)具體線路方案優(yōu)選的評(píng)價(jià)指標(biāo)層次模型。

(2)運(yùn)用組合賦權(quán)法確定各指標(biāo)權(quán)重,優(yōu)化單一賦權(quán)過度依賴主觀判斷或樣本信息的局限性,并結(jié)合二維云理論,構(gòu)建主客觀方法相結(jié)合的二維云線路方案評(píng)價(jià)模型。

(3)通過MATLAB正向云發(fā)生器,直觀反映各線路方案在工程建設(shè)結(jié)果和工后恢復(fù)成本兩方面的隸屬評(píng)價(jià)等級(jí),推薦線路與工程實(shí)際選擇相一致,驗(yàn)證了模型的可行性和有效性,也為沈白鐵路通化段線路方案優(yōu)選提供了理論依據(jù)。

猜你喜歡
評(píng)價(jià)模型
一半模型
SBR改性瀝青的穩(wěn)定性評(píng)價(jià)
石油瀝青(2021年4期)2021-10-14 08:50:44
中藥治療室性早搏系統(tǒng)評(píng)價(jià)再評(píng)價(jià)
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
基于Moodle的學(xué)習(xí)評(píng)價(jià)
關(guān)于項(xiàng)目后評(píng)價(jià)中“專項(xiàng)”后評(píng)價(jià)的探討
保加利亞轉(zhuǎn)軌20年評(píng)價(jià)
主站蜘蛛池模板: 波多野结衣视频网站| 久热中文字幕在线| 九色视频一区| 亚洲欧美成人在线视频| 国产sm重味一区二区三区| 欧洲熟妇精品视频| 久久精品国产精品青草app| av在线5g无码天天| 亚洲全网成人资源在线观看| 欧美精品啪啪一区二区三区| 国产成人亚洲综合A∨在线播放| 精品国产三级在线观看| 最新痴汉在线无码AV| 欧美三级不卡在线观看视频| 亚洲国产欧洲精品路线久久| 成人午夜视频免费看欧美| 亚洲国产高清精品线久久| 欧美午夜理伦三级在线观看| 久久久久久久久久国产精品| 丁香五月亚洲综合在线| 亚洲精品麻豆| 四虎精品免费久久| 中文字幕无码制服中字| 精品久久久久久成人AV| 曰韩人妻一区二区三区| 91在线播放国产| 欧美成人影院亚洲综合图| 国产欧美精品专区一区二区| 国产精品美女在线| 亚洲高清免费在线观看| 在线精品视频成人网| 狠狠色噜噜狠狠狠狠色综合久| 亚洲精品国偷自产在线91正片| a级高清毛片| 国产女人综合久久精品视| a级高清毛片| 中国国产高清免费AV片| 国产精品爆乳99久久| 亚洲人成网站在线观看播放不卡| 欧美成a人片在线观看| 热久久综合这里只有精品电影| 伊人色综合久久天天| 欧美精品亚洲二区| 狠狠操夜夜爽| 日本尹人综合香蕉在线观看| 免费无码AV片在线观看中文| 婷婷综合色| 91福利免费| 黄色网址手机国内免费在线观看| 精品伊人久久大香线蕉网站| 99er精品视频| 高清无码手机在线观看| 精品人妻无码中字系列| 久久国产精品波多野结衣| 午夜国产小视频| 男女男免费视频网站国产| 国产日韩欧美中文| 亚洲手机在线| 18禁影院亚洲专区| 久久综合丝袜日本网| 波多野结衣二区| WWW丫丫国产成人精品| 欧美狠狠干| 国产精品第| 国产精品九九视频| 久久综合色天堂av| 久久精品这里只有国产中文精品| 视频一区视频二区中文精品| 人人艹人人爽| 日韩欧美国产成人| 国产91色| 欧美激情一区二区三区成人| 国产偷国产偷在线高清| 91网在线| 日韩无码视频播放| 中文字幕亚洲无线码一区女同| 国产麻豆91网在线看| 丰满人妻一区二区三区视频| 在线毛片网站| 亚洲第一成人在线| 色老头综合网| 激情综合激情|