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

    基于AHP-OWA-EM-FUZZY模型的西藏高原典型小流域水資源承載力研究

    2022-08-29 05:09:24文梅君彭明勇劉國勇張位首
    水利規(guī)劃與設計 2022年9期
    關鍵詞:典型承載力流域

    文梅君,余 江,彭明勇,劉國勇,張位首

    (1.西藏自治區(qū)水利電力規(guī)劃勘測設計研究院,西藏 拉薩 850000;2.四川大學,四川 成都 610000)

    目前,我國水資源短缺問題突出,主要分為資源型缺水、工程性缺水和水質性缺水3大類[1]。水資源承載力是對水資源安全的一個基本度量,研究水資源承載力對于認識和建設水資源安全保障體系具有重要意義。近年來,國內外學者開展了各種單一模型、復合模型的水資源承載力相關研究[2- 9],但是相關專家、學者對西藏高原流域的水資源承載力研究較少?,F(xiàn)階段關于西藏的雅魯藏布江[10]、拉薩河[11]、瀾滄江[12]、金沙江[13]等大江大河研究比較多,而西藏有大量小流域,且小流域中有許多人口、牛羊、耕地、林地、草地分布,但是針對小流域的研究卻比較少。2014年以來,西藏先后啟動了多個流域,包括中小流域的綜合規(guī)劃編制工作[14],對流域進行統(tǒng)籌管理,迫切需要對西藏高原流域進行水資源承載力深入研究。

    1 西藏高原典型流域概況

    西藏河流眾多,流域面積100km2及以上河流3361條;流域面積1000km2及以上河流331條;流域面積10000km2及以上河流28條[15]。本文選取的6條典型流域包括西藏東部的金沙江右岸一級支流嘎托河流域、西藏中部的雅魯藏布江左岸一級支流湘河流域、西藏西部阿里地區(qū)KHL、西藏北部的怒江二級支流索曲左岸一級支流連曲(巴青曲)、西藏東部怒江一級支流冷曲、西藏東部怒江右岸一級支流昌西曲(舍曲河)。6條流域分布在那曲市、昌都市、林芝市、阿里地區(qū)、日喀則市,占西藏七個地市中的5個地市,地市涉及率達到70%以上。6條流域遍布了西藏東部、中部、西部,因此選取的6條流域具有典型性。

    本文數(shù)據(jù)基準年為2015年,采用《2016年西藏統(tǒng)計年鑒》《2016年昌都市統(tǒng)計年鑒》《2016年日喀則市統(tǒng)計年鑒》《2016年阿里地區(qū)統(tǒng)計年鑒》《2016年林芝市統(tǒng)計年鑒》、流域綜合規(guī)劃和《2015年西藏水資源公報》等數(shù)據(jù),統(tǒng)計6條流域的自然資源、社會經(jīng)濟和水資源狀況。

    2 西藏高原典型流域水資源承載力評價指標體系的建立

    2.1 總體思路

    決定水資源承載力狀況的是一個復雜的耦合的系統(tǒng),這個耦合的系統(tǒng)主要包括資源稟賦、社會經(jīng)濟、生態(tài)環(huán)境等[16]。對流域水資源承載力進行評價就需要定量分析水資源承載力狀況,全面系統(tǒng)分析資源稟賦、社會經(jīng)濟、生態(tài)環(huán)境的相互關系,初步篩選出各個子系統(tǒng)對應的量化指標。結合西藏高原流域特點,篩選出反映西藏高原流域水資源承載力的評價指標。評價指標數(shù)量不宜太多,一般每個系統(tǒng)評價指標建議不超過10個,7—9個最好[17]。對評價指標體系進行量化分級,確定其權重,從而科學合理構建水資源承載力指標體系。

    2.2 指標體系的選取分析

    根據(jù)以上思路,本研究把西藏高原流域水資源承載力指標體系分為3個層次5個大類28項具體指標。評價指標為正向時,其值越大,承載力越強;評價指標為逆向時,其值越大,承載力越弱[18]。3個層次分別為目標層(A)、準則層(B)、指標層(C)。目標層(A)為西藏高原典型流域水資源承載力。準則層的建立本著從宏觀全面反映與水資源承載力相關的資源、社會、經(jīng)濟、生態(tài)等情況,同時結合用水總量控制、用水效率控制、水功能區(qū)限制納污控制這“三條紅線”的要求和最嚴格水資源管理制度的內容,把水功能區(qū)限制納污這條紅線對應指標歸入生態(tài)環(huán)境準則層中,因此準則層分為5個大類,分別為資源稟賦(B1)、社會經(jīng)濟(B2)、生態(tài)環(huán)境(B3)、用水總量(B4)、用水效率(B5)。指標層共計28項,詳見表1。為便于評價指標的分級,反映各指標相對大小,利于橫向比較各流域中相應指標,總體來說大部分指標基本是一個比值。

    2.3 指標體系的分級

    為研究西藏高原典型河流的水資源承載力,需對評價指標建立評價等級,評價等級建立得是否科學合理將影響水資源承載力評價結果的可靠性。

    由于西藏的水資源稟賦、社會經(jīng)濟、生態(tài)環(huán)境、生產(chǎn)生活方式和全國相比具有相當大的差異性,參考現(xiàn)有的水資源承載力指標體系分級研究[19],結合西藏自治區(qū)省級行政區(qū)劃28個評價指標的平均狀況,建立評價指標分級體系,見表2。

    本文假設西藏自治區(qū)省級行政區(qū)的28個評價指標基本處于V2等級,將西藏高原典型流域水資源承載力28個評價指標劃分為V1、V2、V3三個等級,詳見表3。V1表明水資源承載力好,V2表明水資源承載力一般,V3表明水資源承載力差?!?”說明該指標為正向性指標,指標數(shù)值越大,水資源承載狀況越好。“-”說明該指標為負向性指標,指標數(shù)值越大,水資源承載狀況越差。

    表1 西藏高原典型流域水資源承載力評價指標體系

    3 AHP-OWA-EM-FUZZY模型構建

    首先通過層次分析模型(AHP)及單因素方差分析(One Way ANOVA)建立與西藏高原典型流域水資源承載力評價指標體系相應的主觀權重;其次,根據(jù)西藏東部、中部、西部的典型流域的現(xiàn)狀,用熵值模型(Entropy Method)建立指標體系的客觀權重。綜合權重為層次單因素方差模型(AHP-OWA)和熵值模型(EM)的算術平均值。

    把求得的綜合權重作為模糊綜合評價模型(FUZZY)的輸入項,通過構造隸屬度函數(shù)和模糊評價矩陣,計算綜合評價指數(shù)。

    3.1 層次單因素方差分析模型(AHP-OWA)

    層次分析法(AHP)能綜合定性和定量2方面的指標,將人的主觀判斷用數(shù)量表達出來,應用廣泛。由于人的主觀判斷得出的權重可能有很大的差異性,為了避免主觀判斷的差異性,用單因素方差分析法對專家主觀判斷的權重進行分析,以保證主觀權重數(shù)據(jù)的合理性和有效性。

    采用美國運籌學家T.L.Saaty教授于20世紀70年代初期提出的層次分析法(AHP)[20],把建立的研究目標的指標體系表1在yaahp軟件上建立層次結構模型,生成AHP專家調查表,如圖1、圖2所示等,第3層要素調查表在yaahp軟件上。

    咨詢專家填寫AHP模型相對重要性調查表。本研究共咨詢了10位西藏水資源權威專家,分析了水資源承載力并對指標層和準則層的重要性進行判斷。咨詢的專家中,副教授2位,高級工程師2位,工程師6位。

    建立判斷矩陣并求特征值和權重,對判斷矩陣進行一致性檢驗。首先根據(jù)每位專家填寫的AHP模型相對重要性調查表,計算每位專家指標層中指標對A目標層的權重排序。

    表3 西藏高原典型流域水資源承載力評價指標分級

    使用yaahp軟件的群決策對專家數(shù)據(jù)進行分析。設每5個專家為一個群決策組,每個群決策組中有1位副教授、1位高級工程師、3位工程師,因此把10個專家分為群決策組J1、J2。把2組群決策結果的權重進行EXCEL“數(shù)據(jù)分析”中的單因素方差分析(One Way ANOVA),判斷兩組數(shù)據(jù)的差異顯著性水平。如2組數(shù)據(jù)差異性不顯著,則專家主觀判斷取得的權重可信。

    每組群決策計算結果詳見表4。Wj1′為群決策組J1對指標層C的權重,Wj2′為群決策組J2對指標層C的權重。對Wj1′、Wj2′2組數(shù)據(jù)采用算術平均法計算層次模型權重得到Wj′。Wj1′、Wj2′的數(shù)據(jù)趨勢分布如圖3所示。通過分析,群決策組J1、J2由層次分析法得出的權重較大的指標集中在B1、B3、B5準則層,且2組數(shù)據(jù)權重趨勢具有一致性。

    表4 C指標層中指標對A目標層的權重(Wj1′、Wj2′、Wj′)

    由于專家在填寫調查問卷的主觀性強,現(xiàn)用單因素方差分析法對專家主觀判斷的權重進行分析,從而判斷層次分析法2組群決策專家組調查得出的權重是否有效。

    采用EXCEL單因素方差分析,以α=0.05及95%的可靠性推斷,得出方差分析結果,見表5。得出F

    3.2 熵值模型[21]

    熵值法是計算權值的一種客觀方法。在賦值過程中沒有人為主觀因素加入,依靠實測和統(tǒng)計資料通過對數(shù)據(jù)計算確定評價指標權重,能相對真實地反映出指標之間的關系。確定權重的步驟如下:

    表5 Wj1′,Wj2′單因素方差分析結果表

    圖1 西藏高原典型流域水資源承載力調查表

    圖2 西藏高原典型流域水資源承載力調查表

    圖3 Wj1′,Wj2′趨勢分布圖

    (1)建立有m個樣本、n個評價指標實測或統(tǒng)計的數(shù)據(jù)的矩陣:

    R=(xij)m×n(i=1,2…,n;j=1,2,…m)

    (1)

    (2)把樣本數(shù)據(jù)矩陣進行非負化處理,同時為了避免求熵值得出無意義的數(shù)據(jù),同時平移數(shù)據(jù),非負化公式如下:

    (2)

    式中,指標bij越大越優(yōu)。

    (3)

    式中,指標bij越小越優(yōu)。

    (3)由B=(bij)m×n計算j項指標下第i個樣本占該指標的比重Pij:

    (4)

    (4)第j個評價指標的熵為:

    (5)

    (6)

    (5)計算評價指標的權重:

    (7)

    各指標的權重滿足指標的熵越大,其熵權越小,該指標越不重要,而且滿足:

    (8)

    3.3 綜合權重及計算結果

    AHP模型得出的各評價指標的權重為專家、學者的主觀判斷。熵值模型計算的各評價指標為典型流域的客觀權重。為使權重的確立更加科學有效,因此用綜合權重計算評價指標。綜合權重為層次單因素方差模型(AHP-OWA)和熵值模型(EM)的算術平均值。

    計算公式如下:

    (9)

    評價指標的綜合權重結果如下:

    B1、B2、B3、B4、B5對目標A的權重為WA=(0.2026,0.1555,0.2209,0.1576,0.2635);

    在資源稟賦B1中,C1、C2、C3、C4、C5、C6對目標B1的權重為WB1=(0.2195,0.1767,0.2602,0.1418,0.1758,0.0260);

    在社會經(jīng)濟B2中,C7、C8、C9、C10、C11、C12對目標B2的權重為WB2=(0.1203,0.1571,0.1781,0.2281,0.1404,0.1759);

    在生態(tài)環(huán)境B3中,C13、C14、C15、C16、C17對目標B3的權重為WB3=(0.1577,0.3147,0.1408,0.1674,0.2194);

    在用水總量B4中,C18、C19、C20、C21、C22對目標B4的權重為WB4=(0.2352,0.2064,0.2151,0.1877,0.1556);

    在用水效率B5中,C23、C24、C25、C26、C27、C28對目標B5的權重為WB5=(0.1384,0.1473,0.1828,0.1951,0.1702,0.1662);

    評價指標層C對目標層A的綜合權重WAC=(0.0445,0.0358,0.0527,0.0287,0.0356,0.0053,0.0187,0.0244,0.0277,0.0355,.0218,0.0274,0.0348,0.0695,0.0311,0.0370,0.0485,0.0371,0.0325,0.0339,0.0296,0.0245,0.0365,0.0388,0.0481,0.0514,0.0448,0.0438)。

    其中綜合權重大于0.04以上的從高到低依次為:C14,C3,C26,C17,C25,C27,C1,C28,這些指標屬于資源稟賦、生態(tài)環(huán)境和用水效率準則層,說明在西藏高原典型流域水資源承載力影響因素中,B1、B3、B5的作用更為突出。

    3.4 模糊綜合評價模型

    1965年,美國L.A.Zadeh教授提出了模糊集合論[22]。模糊綜合評價法是以模糊數(shù)學隸屬度的基本理論為基礎,為了某種目標,考慮多種因素的綜合影響,對復雜問題做出決定的綜合評價方法??偟膩碚f,模糊綜合評價法是通過評價因素集、對應評語集通過隸屬度函數(shù)構造模糊矩陣,從而綜合了評價對象的評價因素信息,得出評價結果,避免了單指標評價信息流失[23]。

    3.4.1設立評價因素集

    建立西藏高原典型流域水資源承載力評價因素集。因素集是評價對象的各評價因素的集合。設C={C1,C2,…Cn},Ci(i=1,2,…n)為評價因素。本文根據(jù)西藏高原典型流域的特點,設置28個評價指標(因素)。其中A代表西藏高原典型流域水資源承載力這個目標。設A=(B1,B2,B3,B4,B5),其中Bi為準則層構建的評價因素集。

    設B1=(C1,C2,C3,C4,C5,C6);

    B2=(C7,C8,C9,C10,C11,C12);

    B3=(C13,C14,C15,C16,C17);

    (10)

    B4=(C18,C19,C20,C21,C22);

    B5=(C23,C24,C25,C26,C27,C28)。

    3.4.2設定評價因素的評判集

    評價對象為西藏高原典型流域,包括昌西曲、嘎托河、冷曲、連曲、湘河、KHL。對評價對象根據(jù)評價因素屬于評價分級標準的相應等級做出評語的集合為評判集,如下:

    V=(V1,V2,…Vm)

    (11)

    式中,V—C中各因素對應評價等級的集合。

    3.4.3建立評價因素的權重

    評價因素的權重為層次單因素方差模型(AHP-OWA)和熵值模型(EM)的算術平均值求得的綜合權重:

    (12)

    3.4.4構造等級之間能平滑過渡的隸屬度函數(shù)[24]

    隸屬度函數(shù)能使各評價指標相差不大的數(shù)值或者相差很大的數(shù)值通過模糊化處理后,各指標能在等級間平滑過渡[25]。本文評判集劃分為3個等級V=(V1,V2,V3)。

    因此設評價因素對評判集的隸屬度函數(shù)為[26]:

    (13)

    (14)

    (15)

    設k1為V1與V2的臨界值,k3為V2與V3的臨界值,k2為V2水平區(qū)間中點值[27]。當評價因素越小,水資源承載力越好時,采用式子(13)至(15)計算隸屬度函數(shù)。如評價因素越大,水資源承載力越好時,需將上式右邊的所有區(qū)間號改變方向,即大于改成小于,小于等于改成大于等于,再用來計算評價因素的隸屬度函數(shù)。

    3.4.5根據(jù)隸屬函數(shù)構造評價因素與評判集的隸屬度矩陣

    根據(jù)隸屬度函數(shù),定義評價因素集中B中的第i個元素Ci(i=1,2,…,n)對應評判集中的第j個元素(j=1,2,…,m)的隸屬度為rij,把評價因素集為維度構造模糊矩陣R如下:

    (16)

    3.4.6模糊綜合評價計算

    由加權模型,模糊綜合評價可以表示為:

    Z=WR,Z=(z1,z2,…zj…zn),zj∈[0,1]

    (17)

    式中,W—下一層指標相對于上一層指標的權重;zj—等級Vj對Z的隸屬度,表示模糊綜合評價結果。

    再賦予評判集中各等級分值aj,設V1等級中a1=0.95,V2等級中a2=0.50,V3等級中a3=0.05。用下式計算水資源承載力的綜合評分值及水資源承載力準則層的綜合評分值a。a值越高,水資源承載力的潛力越大。根據(jù)AHP-OWA-EM-FUZZY模型及評價指標分級標準設計,V2等級約為西藏水資源承載力平均水平。

    (18)

    4 評價結果及分析

    4.1 6條典型流域水資源承載力的綜合評價結果

    根據(jù)AHP-OWA-EM-FUZZY模型,由公式(17)可以求得5個準則層的模糊綜合評價結果,見表6—11。

    表6 嘎托河流域水資源承載力模糊綜合評價結果表

    表7 湘河流域水資源承載力模糊綜合評價結果表

    表8 KHL流域水資源承載力模糊綜合評價結果表

    表9 昌西曲流域水資源承載力模糊綜合評價結果表

    表10 連曲流域水資源承載力模糊綜合評價結果表

    表11 冷曲流域水資源承載力模糊綜合評價結果表

    4.2 6條典型流域水資源承載力綜合評價結果分析

    4.2.1從流域屬于的等級分析

    根據(jù)表6—11,Z值對V1等級隸屬度最大的有昌西曲、嘎托河、連曲,說明昌西曲、嘎托河、連曲水資源承載力處于較好的狀態(tài);Z值對V2等級隸屬度最大的流域有冷曲、KHL,說明冷曲、KHL水資源承載力處于一般的狀態(tài);Z值對V3等級隸屬度最大的流域有湘河,說明湘河水資源承載力處于較差的狀態(tài)。

    4.2.2從流域水資源承載力模糊綜合評價得分分析

    西藏高原6條典型流域水資源承載力得分詳見表12。得分越高,水資源承載力越強。西藏高原典型流域水資源承載力模糊綜合評價得分從高到低依次為昌西曲、連曲、嘎托河、KHL、冷曲、湘河,說明昌西曲流域水資源承載力最好,湘河水資源承載力相對比較差。

    表12 西藏高原典型流域水資源承載力模糊綜合評價得分表

    4.2.3從流域所處的地理位置分析

    西藏地域面積廣闊,資源稟賦、社會經(jīng)濟、生態(tài)環(huán)境狀況從東到西、從北到南差異性極大。本文把西藏高原6條典型流域從林芝市所在地的經(jīng)度分為東面3條流域、西面3條流域分別分析其水資源承載力。林芝市東面3條流域所處地理位置為藏東,從南到北分別為昌西曲、嘎托河、冷曲。林芝市西面3條流域所處地理位置為藏中和藏西北部,從南到北分別為湘河、KHL、連曲。在藏東從南向北,昌西曲、嘎托河、冷曲的水資源承載力逐漸降低。在藏中和藏西北部從南向北,湘河、KHL、連曲的水資源承載力逐漸升高。如圖4所示。

    圖4 西藏高原典型流域水資源承載力趨勢圖

    4.2.4從準則層分析

    從準則層分析,社會經(jīng)濟、用水總量得分相對較高。資源稟賦、生態(tài)環(huán)境、用水效率得分相對較低,如需提高西藏高原流域水資源承載力的安全保障能力,需從這3個方面著力。

    由于水資源承載力是由多系統(tǒng)綜合作用決定的,難以簡單歸納說明某一方面的因素決定水資源承載力高低。但是本文研究發(fā)現(xiàn),西藏高原6條典型流域水資源承載力高低和用水效率水平高低呈現(xiàn)基本一致性。

    5 結論

    本文構建適合西藏高原不同流域的水資源承載力評價模型(AHP-OWA-EM-FUZZY),通過建立評價指標量化評價研究西藏高原典型小流域的水資源承載力,對西藏高原小流域的水資源優(yōu)化配置、提高高原流域的水資源承載能力、更有效地提高農(nóng)牧民生活用水、一產(chǎn)二產(chǎn)三產(chǎn)用水、生態(tài)用水的保證率具有很好的借鑒作用,為我國水資源承載能力評估指南建立提供基礎資料。

    由于西藏高原流域基礎數(shù)據(jù)的缺乏本文僅選取了西藏高原6條典型小流域進行空間研究。在空間區(qū)域上,未來可以增加西藏高原典型流域的數(shù)量,以建立更完善的適合西藏高原典型流域水資源承載力體系,研究各流域水資源承載力的差異性及共性,更好地為流域綜合管理與規(guī)劃服務。

    猜你喜歡
    典型承載力流域
    用最典型的事寫最有特點的人
    壓油溝小流域
    多項式求值題的典型解法
    堡子溝流域綜合治理
    羅堰小流域
    典型胰島素瘤1例報道
    打造智慧流域的思路及構想——以討賴河流域為例
    CFRP-PCP板加固混凝土梁的抗彎承載力研究
    耐火鋼圓鋼管混凝土柱耐火極限和承載力
    潛艇極限承載力計算與分析
    免费在线观看视频国产中文字幕亚洲| 最新的欧美精品一区二区| 日本撒尿小便嘘嘘汇集6| 搡老乐熟女国产| 欧美乱色亚洲激情| 欧美+亚洲+日韩+国产| 亚洲精品一二三| 国产av精品麻豆| 欧美日韩视频精品一区| 国产在线精品亚洲第一网站| 天堂俺去俺来也www色官网| 国产不卡一卡二| 俄罗斯特黄特色一大片| 亚洲成a人片在线一区二区| 日日夜夜操网爽| 变态另类成人亚洲欧美熟女 | 一级毛片精品| 亚洲专区中文字幕在线| 人人妻人人澡人人爽人人夜夜| 亚洲人成电影观看| 欧美久久黑人一区二区| 婷婷成人精品国产| 深夜精品福利| 最新的欧美精品一区二区| 黑人猛操日本美女一级片| 久久国产亚洲av麻豆专区| 老司机深夜福利视频在线观看| 国产亚洲精品久久久久久毛片 | 黄色怎么调成土黄色| 国产精品亚洲一级av第二区| 亚洲精品美女久久久久99蜜臀| 又大又爽又粗| 人人妻人人添人人爽欧美一区卜| 中文字幕制服av| 日本精品一区二区三区蜜桃| 欧美精品一区二区免费开放| 免费在线观看日本一区| www.自偷自拍.com| 中文字幕人妻丝袜制服| 久久久国产一区二区| 免费在线观看视频国产中文字幕亚洲| 日韩精品免费视频一区二区三区| cao死你这个sao货| 欧美精品av麻豆av| 久久国产精品影院| 99精品欧美一区二区三区四区| 十八禁人妻一区二区| 丰满人妻熟妇乱又伦精品不卡| 国产精品一区二区精品视频观看| 精品久久久久久,| 一个人免费在线观看的高清视频| 大码成人一级视频| 亚洲精品一卡2卡三卡4卡5卡| 国产真人三级小视频在线观看| 91av网站免费观看| av网站在线播放免费| 欧美人与性动交α欧美软件| 日韩精品免费视频一区二区三区| 又大又爽又粗| 国产99久久九九免费精品| 老司机靠b影院| 十八禁网站免费在线| 波多野结衣一区麻豆| 国产成人av激情在线播放| 丰满饥渴人妻一区二区三| av视频免费观看在线观看| 国产成人精品无人区| 亚洲成国产人片在线观看| 国产精品久久久久成人av| 国产野战对白在线观看| 满18在线观看网站| 国产日韩一区二区三区精品不卡| ponron亚洲| 99在线人妻在线中文字幕 | 天天躁狠狠躁夜夜躁狠狠躁| 夜夜爽天天搞| 欧美久久黑人一区二区| 成人国产一区最新在线观看| 免费女性裸体啪啪无遮挡网站| 国产亚洲欧美98| 女人爽到高潮嗷嗷叫在线视频| 午夜影院日韩av| 精品国产一区二区三区久久久樱花| 日韩一卡2卡3卡4卡2021年| 少妇猛男粗大的猛烈进出视频| 大型黄色视频在线免费观看| 中文字幕人妻丝袜制服| 正在播放国产对白刺激| 亚洲欧美精品综合一区二区三区| 国产欧美日韩精品亚洲av| 成人黄色视频免费在线看| 中文字幕最新亚洲高清| 欧美 亚洲 国产 日韩一| 中文字幕制服av| 美女午夜性视频免费| 亚洲aⅴ乱码一区二区在线播放 | 麻豆乱淫一区二区| 亚洲精品粉嫩美女一区| 欧美黄色淫秽网站| 午夜激情av网站| 精品亚洲成a人片在线观看| 亚洲av成人不卡在线观看播放网| 亚洲av成人不卡在线观看播放网| 午夜免费成人在线视频| 国产主播在线观看一区二区| 久久中文字幕人妻熟女| 亚洲人成电影免费在线| 日韩欧美在线二视频 | 99精品欧美一区二区三区四区| 欧美激情高清一区二区三区| 国产国语露脸激情在线看| 欧洲精品卡2卡3卡4卡5卡区| 99国产精品一区二区三区| 久久国产乱子伦精品免费另类| 女人久久www免费人成看片| 国产欧美日韩一区二区精品| 成人18禁在线播放| 日日爽夜夜爽网站| 亚洲午夜精品一区,二区,三区| 亚洲精品中文字幕在线视频| 午夜影院日韩av| 99热只有精品国产| 国产精品一区二区在线观看99| 精品福利永久在线观看| 亚洲熟女精品中文字幕| 成人国语在线视频| 两个人看的免费小视频| 无限看片的www在线观看| 老司机在亚洲福利影院| 丝袜在线中文字幕| 少妇粗大呻吟视频| 国产成人系列免费观看| 久久精品aⅴ一区二区三区四区| 欧美黄色片欧美黄色片| 久久精品成人免费网站| 国内久久婷婷六月综合欲色啪| av不卡在线播放| 韩国精品一区二区三区| 亚洲成av片中文字幕在线观看| 午夜福利在线观看吧| 老司机午夜十八禁免费视频| 欧美精品av麻豆av| 亚洲欧美一区二区三区久久| 动漫黄色视频在线观看| 精品人妻熟女毛片av久久网站| 国产成人av教育| 欧美乱色亚洲激情| www.999成人在线观看| 久久这里只有精品19| 国产免费男女视频| 99热只有精品国产| 精品久久蜜臀av无| 热re99久久精品国产66热6| 欧美日韩中文字幕国产精品一区二区三区 | 91麻豆精品激情在线观看国产 | 看免费av毛片| 国产在线观看jvid| 两性夫妻黄色片| 9191精品国产免费久久| 一区二区日韩欧美中文字幕| 99热国产这里只有精品6| 国产成+人综合+亚洲专区| 精品一品国产午夜福利视频| 免费观看a级毛片全部| 欧美日韩国产mv在线观看视频| 日韩精品免费视频一区二区三区| 国产成人av教育| 老熟妇乱子伦视频在线观看| 欧美激情极品国产一区二区三区| 成人三级做爰电影| 国产深夜福利视频在线观看| a在线观看视频网站| 黄色毛片三级朝国网站| 欧美午夜高清在线| 亚洲成a人片在线一区二区| 97人妻天天添夜夜摸| 十八禁高潮呻吟视频| 久久99一区二区三区| 最近最新中文字幕大全电影3 | 一区在线观看完整版| 欧美乱码精品一区二区三区| 99re6热这里在线精品视频| 久久久久久亚洲精品国产蜜桃av| 国产无遮挡羞羞视频在线观看| 亚洲人成77777在线视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区中文字幕在线| 高清黄色对白视频在线免费看| 欧美久久黑人一区二区| 国产又爽黄色视频| 日本wwww免费看| 手机成人av网站| 国产欧美亚洲国产| 久久久国产一区二区| 婷婷成人精品国产| 一夜夜www| 国产在线精品亚洲第一网站| 久久精品熟女亚洲av麻豆精品| 别揉我奶头~嗯~啊~动态视频| 成人永久免费在线观看视频| 少妇被粗大的猛进出69影院| 中文字幕色久视频| 午夜精品国产一区二区电影| 中国美女看黄片| 精品卡一卡二卡四卡免费| 久久亚洲精品不卡| 欧美黄色淫秽网站| 国产亚洲一区二区精品| 久久久久久久精品吃奶| 丝袜美足系列| 欧美黑人精品巨大| 亚洲一区高清亚洲精品| 亚洲专区中文字幕在线| 母亲3免费完整高清在线观看| 久久 成人 亚洲| 午夜福利一区二区在线看| 免费观看a级毛片全部| 欧美日韩瑟瑟在线播放| 午夜福利免费观看在线| 怎么达到女性高潮| 1024香蕉在线观看| 国产男女内射视频| 极品少妇高潮喷水抽搐| 99热国产这里只有精品6| 亚洲色图av天堂| 国产三级黄色录像| 黄片播放在线免费| 久久精品aⅴ一区二区三区四区| 国产成人一区二区三区免费视频网站| 亚洲欧美精品综合一区二区三区| 51午夜福利影视在线观看| 亚洲情色 制服丝袜| 久久久国产欧美日韩av| 日本wwww免费看| 久久香蕉精品热| 国产亚洲精品久久久久5区| 亚洲成国产人片在线观看| 成人特级黄色片久久久久久久| 俄罗斯特黄特色一大片| 国产成人av激情在线播放| 啦啦啦免费观看视频1| 99香蕉大伊视频| av免费在线观看网站| 色94色欧美一区二区| av一本久久久久| 新久久久久国产一级毛片| 成人手机av| 超碰成人久久| 久久久国产精品麻豆| 成人18禁高潮啪啪吃奶动态图| 国产日韩欧美亚洲二区| 99香蕉大伊视频| www日本在线高清视频| 亚洲精品美女久久av网站| 国产成人啪精品午夜网站| 午夜免费成人在线视频| 18禁黄网站禁片午夜丰满| 黄色成人免费大全| 日本a在线网址| 国产1区2区3区精品| 午夜福利一区二区在线看| 在线av久久热| 老熟女久久久| 99热国产这里只有精品6| 国产视频一区二区在线看| 51午夜福利影视在线观看| 啦啦啦视频在线资源免费观看| 国产蜜桃级精品一区二区三区 | 国产欧美日韩一区二区三| 国产成人免费观看mmmm| 999久久久国产精品视频| 午夜福利在线观看吧| 国产成+人综合+亚洲专区| 亚洲精品粉嫩美女一区| 国产在线精品亚洲第一网站| 午夜福利乱码中文字幕| 大码成人一级视频| 国产一区有黄有色的免费视频| 国产精品 国内视频| 免费av中文字幕在线| 一区二区三区激情视频| 一级,二级,三级黄色视频| 欧美人与性动交α欧美软件| 老汉色av国产亚洲站长工具| 午夜老司机福利片| 欧美精品啪啪一区二区三区| 热99久久久久精品小说推荐| 欧美成人免费av一区二区三区 | 国产精品成人在线| 亚洲精品成人av观看孕妇| 黄频高清免费视频| 精品卡一卡二卡四卡免费| 精品人妻熟女毛片av久久网站| 在线播放国产精品三级| 久久精品aⅴ一区二区三区四区| 亚洲一区高清亚洲精品| 91成人精品电影| 日本一区二区免费在线视频| 午夜精品国产一区二区电影| 99精品久久久久人妻精品| 黄网站色视频无遮挡免费观看| av网站免费在线观看视频| 91九色精品人成在线观看| 99在线人妻在线中文字幕 | 久久国产精品大桥未久av| 亚洲色图综合在线观看| 在线天堂中文资源库| www.熟女人妻精品国产| 三上悠亚av全集在线观看| 久久精品国产亚洲av高清一级| 成人精品一区二区免费| 久久久久久久午夜电影 | 悠悠久久av| 中国美女看黄片| 岛国在线观看网站| 亚洲七黄色美女视频| 亚洲欧美一区二区三区久久| 国产av精品麻豆| 一级毛片高清免费大全| 久久久国产精品麻豆| 99精品在免费线老司机午夜| 国产一区二区三区综合在线观看| 久久久久久久久免费视频了| 亚洲成av片中文字幕在线观看| 国产人伦9x9x在线观看| 另类亚洲欧美激情| 高清在线国产一区| 精品人妻在线不人妻| 亚洲第一欧美日韩一区二区三区| 无遮挡黄片免费观看| 久久久精品国产亚洲av高清涩受| 99国产极品粉嫩在线观看| 精品熟女少妇八av免费久了| 久久久久久人人人人人| 女警被强在线播放| 亚洲全国av大片| 国产欧美日韩一区二区精品| 一进一出抽搐gif免费好疼 | 精品国内亚洲2022精品成人 | 成人18禁在线播放| 国产精品久久久人人做人人爽| 免费久久久久久久精品成人欧美视频| 午夜福利乱码中文字幕| 免费av中文字幕在线| 亚洲精品久久成人aⅴ小说| 大片电影免费在线观看免费| 久热这里只有精品99| 精品午夜福利视频在线观看一区| tube8黄色片| 国产亚洲欧美98| av中文乱码字幕在线| 在线观看免费午夜福利视频| 中文字幕制服av| 在线观看日韩欧美| 大型黄色视频在线免费观看| 老司机午夜十八禁免费视频| 免费看a级黄色片| 久久久精品国产亚洲av高清涩受| 国产日韩欧美亚洲二区| 大码成人一级视频| 亚洲精品久久成人aⅴ小说| av网站在线播放免费| 国产不卡一卡二| 久久久久国产一级毛片高清牌| 午夜免费成人在线视频| 欧美黄色淫秽网站| 国产真人三级小视频在线观看| 欧美成人免费av一区二区三区 | 操美女的视频在线观看| 亚洲欧美色中文字幕在线| 久久草成人影院| 亚洲第一av免费看| 欧美另类亚洲清纯唯美| 久久香蕉激情| 嫩草影视91久久| 日本wwww免费看| av线在线观看网站| 十分钟在线观看高清视频www| 涩涩av久久男人的天堂| 女人被狂操c到高潮| 国产深夜福利视频在线观看| 夜夜夜夜夜久久久久| 久久久国产欧美日韩av| 中出人妻视频一区二区| 午夜福利乱码中文字幕| 久久国产精品人妻蜜桃| www.999成人在线观看| 人妻丰满熟妇av一区二区三区 | 免费久久久久久久精品成人欧美视频| 国产又色又爽无遮挡免费看| 成人18禁在线播放| 日韩欧美三级三区| 脱女人内裤的视频| 美女高潮到喷水免费观看| 久久 成人 亚洲| 91成年电影在线观看| 在线观看舔阴道视频| 精品国产乱码久久久久久男人| 好看av亚洲va欧美ⅴa在| 男女之事视频高清在线观看| 老司机在亚洲福利影院| 一区二区三区精品91| 热re99久久精品国产66热6| 亚洲色图 男人天堂 中文字幕| 老司机靠b影院| 在线观看日韩欧美| av天堂在线播放| 女人被狂操c到高潮| 久久国产精品男人的天堂亚洲| 大片电影免费在线观看免费| 久久天堂一区二区三区四区| 久久国产精品人妻蜜桃| 黄色视频不卡| av国产精品久久久久影院| 麻豆av在线久日| 女同久久另类99精品国产91| 午夜精品在线福利| 99久久人妻综合| 成人国语在线视频| 日韩欧美一区视频在线观看| 久久ye,这里只有精品| 国产精品二区激情视频| 久久青草综合色| 国产欧美日韩精品亚洲av| 亚洲色图综合在线观看| 欧美日韩亚洲综合一区二区三区_| 久久人人97超碰香蕉20202| 国产成人av教育| 国产野战对白在线观看| av网站在线播放免费| 国产精品久久久人人做人人爽| 国产在线精品亚洲第一网站| 午夜福利免费观看在线| 欧美精品啪啪一区二区三区| 超色免费av| 亚洲精品中文字幕在线视频| 日韩一卡2卡3卡4卡2021年| 一级a爱片免费观看的视频| 夫妻午夜视频| 亚洲精品美女久久久久99蜜臀| 亚洲 欧美一区二区三区| 露出奶头的视频| 亚洲专区字幕在线| 狠狠狠狠99中文字幕| 国产精品99久久99久久久不卡| 欧美亚洲 丝袜 人妻 在线| 婷婷成人精品国产| 老司机午夜十八禁免费视频| 精品国产一区二区三区四区第35| 国产精品一区二区在线不卡| 99在线人妻在线中文字幕 | 国产一区二区激情短视频| 淫妇啪啪啪对白视频| 国产91精品成人一区二区三区| 看免费av毛片| 俄罗斯特黄特色一大片| 国产精品98久久久久久宅男小说| 国产熟女午夜一区二区三区| 精品久久久久久,| 美女扒开内裤让男人捅视频| 99热只有精品国产| 王馨瑶露胸无遮挡在线观看| 亚洲黑人精品在线| 动漫黄色视频在线观看| 性少妇av在线| 亚洲片人在线观看| 亚洲在线自拍视频| 老司机影院毛片| 久久精品亚洲熟妇少妇任你| 嫁个100分男人电影在线观看| 日韩免费av在线播放| 亚洲免费av在线视频| 久久精品aⅴ一区二区三区四区| 亚洲在线自拍视频| 99国产综合亚洲精品| 中文字幕精品免费在线观看视频| 色精品久久人妻99蜜桃| 成人av一区二区三区在线看| 久久久久视频综合| 免费在线观看完整版高清| 亚洲中文日韩欧美视频| 大型av网站在线播放| 丰满饥渴人妻一区二区三| 女警被强在线播放| 久久狼人影院| 脱女人内裤的视频| 亚洲欧美一区二区三区黑人| 日韩中文字幕欧美一区二区| 色婷婷av一区二区三区视频| 国产免费av片在线观看野外av| 国产精品电影一区二区三区 | 国产亚洲精品久久久久久毛片 | 午夜福利乱码中文字幕| 国产精品电影一区二区三区 | 久久久久久免费高清国产稀缺| 久久久精品国产亚洲av高清涩受| 亚洲精品久久成人aⅴ小说| 怎么达到女性高潮| 亚洲欧美一区二区三区黑人| 999久久久精品免费观看国产| 免费在线观看亚洲国产| 在线永久观看黄色视频| 国产成人精品久久二区二区免费| 大型av网站在线播放| 婷婷成人精品国产| 精品欧美一区二区三区在线| 99国产精品免费福利视频| 日韩免费高清中文字幕av| 亚洲五月天丁香| 飞空精品影院首页| 人妻久久中文字幕网| 91精品三级在线观看| 日韩免费高清中文字幕av| 亚洲人成电影观看| 国产亚洲精品第一综合不卡| 91大片在线观看| 色94色欧美一区二区| 国产精品99久久99久久久不卡| 人人妻人人澡人人爽人人夜夜| 岛国在线观看网站| 一级作爱视频免费观看| 欧美在线黄色| 人妻丰满熟妇av一区二区三区 | 99久久99久久久精品蜜桃| av超薄肉色丝袜交足视频| 天天躁夜夜躁狠狠躁躁| 国产成人啪精品午夜网站| 国产精品九九99| 女性被躁到高潮视频| 久9热在线精品视频| 国产成人系列免费观看| 高潮久久久久久久久久久不卡| 国产亚洲精品久久久久5区| 激情视频va一区二区三区| 国产精品亚洲av一区麻豆| a级片在线免费高清观看视频| 99久久99久久久精品蜜桃| 免费看十八禁软件| 国产日韩欧美亚洲二区| xxx96com| 日本a在线网址| 亚洲国产欧美日韩在线播放| 一级片'在线观看视频| 十八禁网站免费在线| 亚洲成人手机| 一级,二级,三级黄色视频| av线在线观看网站| av天堂久久9| 国内毛片毛片毛片毛片毛片| 精品国产超薄肉色丝袜足j| 久久久国产成人精品二区 | 欧美日韩福利视频一区二区| 国产精品香港三级国产av潘金莲| 国产1区2区3区精品| 亚洲精华国产精华精| 久久人人97超碰香蕉20202| www.熟女人妻精品国产| 99精国产麻豆久久婷婷| 日韩中文字幕欧美一区二区| 亚洲第一欧美日韩一区二区三区| 美女高潮喷水抽搐中文字幕| 国产免费男女视频| 中文字幕人妻熟女乱码| 80岁老熟妇乱子伦牲交| 天天躁日日躁夜夜躁夜夜| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩黄片免| 久久久精品国产亚洲av高清涩受| 欧美老熟妇乱子伦牲交| 99国产精品一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利乱码中文字幕| 亚洲午夜理论影院| 成人手机av| 午夜福利欧美成人| 99国产精品99久久久久| 18禁裸乳无遮挡免费网站照片 | 免费在线观看亚洲国产| 欧美黄色片欧美黄色片| 久久久久国产精品人妻aⅴ院 | 久久精品熟女亚洲av麻豆精品| 精品高清国产在线一区| 在线观看免费日韩欧美大片| 精品一区二区三卡| 美女福利国产在线| 在线观看免费日韩欧美大片| 日韩 欧美 亚洲 中文字幕| 热re99久久国产66热| 黄片大片在线免费观看| 99re在线观看精品视频| 色播在线永久视频| 成年人黄色毛片网站| av欧美777| 韩国av一区二区三区四区| 后天国语完整版免费观看| 亚洲精品久久成人aⅴ小说| 亚洲欧美一区二区三区久久| 亚洲熟女精品中文字幕| 国产蜜桃级精品一区二区三区 | 嫁个100分男人电影在线观看| 免费少妇av软件| 看黄色毛片网站| 超色免费av| 午夜福利欧美成人| 黑人操中国人逼视频| tocl精华| 亚洲精品久久成人aⅴ小说| 91九色精品人成在线观看| 免费观看a级毛片全部| 国产精品久久久人人做人人爽| 无限看片的www在线观看| 国产欧美日韩综合在线一区二区| 宅男免费午夜| 女性被躁到高潮视频| aaaaa片日本免费|