• <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板加固混凝土梁的抗彎承載力研究
    耐火鋼圓鋼管混凝土柱耐火極限和承載力
    潛艇極限承載力計算與分析
    国产精品一区二区在线不卡| 激情视频va一区二区三区| 成年av动漫网址| 欧美变态另类bdsm刘玥| 夫妻性生交免费视频一级片| 极品人妻少妇av视频| xxxhd国产人妻xxx| 婷婷成人精品国产| 久久精品国产亚洲av涩爱| 国产免费福利视频在线观看| 国产一区亚洲一区在线观看| 最后的刺客免费高清国语| 久久久久久久国产电影| 国产精品.久久久| 国产 一区精品| 两个人看的免费小视频| 欧美 日韩 精品 国产| 日韩成人伦理影院| 制服诱惑二区| 黄色毛片三级朝国网站| 久久精品国产亚洲av天美| 亚洲国产日韩一区二区| 日韩制服骚丝袜av| 国产精品久久久久久av不卡| 乱码一卡2卡4卡精品| 亚洲经典国产精华液单| av片东京热男人的天堂| 只有这里有精品99| 丰满迷人的少妇在线观看| 亚洲欧美一区二区三区国产| 国产不卡av网站在线观看| 日本-黄色视频高清免费观看| av片东京热男人的天堂| 中文乱码字字幕精品一区二区三区| av在线app专区| 国产白丝娇喘喷水9色精品| 观看av在线不卡| 久久ye,这里只有精品| 欧美精品亚洲一区二区| 精品一品国产午夜福利视频| 91午夜精品亚洲一区二区三区| 亚洲综合色网址| 99久久精品国产国产毛片| 亚洲少妇的诱惑av| 午夜福利在线观看免费完整高清在| 日本午夜av视频| 搡女人真爽免费视频火全软件| 桃花免费在线播放| 在线观看国产h片| 日韩成人av中文字幕在线观看| 久久97久久精品| 丝袜美足系列| 亚洲五月色婷婷综合| 建设人人有责人人尽责人人享有的| 999精品在线视频| 欧美成人午夜免费资源| 免费看不卡的av| 欧美人与善性xxx| 欧美激情极品国产一区二区三区 | av又黄又爽大尺度在线免费看| 国产男人的电影天堂91| 亚洲精品久久午夜乱码| 国产精品嫩草影院av在线观看| 欧美精品亚洲一区二区| 国产69精品久久久久777片| 欧美成人精品欧美一级黄| 久久久久久人人人人人| 国产日韩欧美亚洲二区| 欧美激情 高清一区二区三区| 成年人免费黄色播放视频| 成人亚洲欧美一区二区av| 制服丝袜香蕉在线| 国产黄色免费在线视频| 制服人妻中文乱码| 秋霞在线观看毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男的添女的下面高潮视频| 自拍欧美九色日韩亚洲蝌蚪91| 春色校园在线视频观看| 一区二区日韩欧美中文字幕 | av卡一久久| 99久久人妻综合| 欧美最新免费一区二区三区| 国产精品人妻久久久影院| 欧美人与性动交α欧美软件 | 狂野欧美激情性xxxx在线观看| 一边摸一边做爽爽视频免费| 国产精品.久久久| 精品国产一区二区久久| 日本欧美国产在线视频| 18+在线观看网站| 亚洲精品久久久久久婷婷小说| av天堂久久9| 日韩一区二区三区影片| 久久久久久久久久久免费av| 国产一区亚洲一区在线观看| 成年人午夜在线观看视频| 国产淫语在线视频| 久久99蜜桃精品久久| 免费大片18禁| 亚洲伊人色综图| 99精国产麻豆久久婷婷| 日本猛色少妇xxxxx猛交久久| 国产精品99久久99久久久不卡 | 老熟女久久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 在线观看美女被高潮喷水网站| 天堂8中文在线网| 最近手机中文字幕大全| a级毛色黄片| 国产色婷婷99| av线在线观看网站| 国产精品麻豆人妻色哟哟久久| 亚洲成色77777| 美女国产视频在线观看| 国产 精品1| a级片在线免费高清观看视频| 日韩av不卡免费在线播放| 青春草亚洲视频在线观看| 韩国高清视频一区二区三区| 大陆偷拍与自拍| 91午夜精品亚洲一区二区三区| av免费观看日本| 亚洲精品乱久久久久久| 国精品久久久久久国模美| 女性被躁到高潮视频| 热re99久久精品国产66热6| 亚洲精品一区蜜桃| 欧美日本中文国产一区发布| 亚洲精品久久久久久婷婷小说| 日韩制服骚丝袜av| 一二三四中文在线观看免费高清| 51国产日韩欧美| 久久久久久人妻| 久久婷婷青草| 肉色欧美久久久久久久蜜桃| 免费观看a级毛片全部| 你懂的网址亚洲精品在线观看| 久久影院123| 少妇的逼好多水| 女性被躁到高潮视频| 美女视频免费永久观看网站| 亚洲精品中文字幕在线视频| 国产精品人妻久久久久久| 亚洲人成77777在线视频| 免费在线观看黄色视频的| 久久青草综合色| 国产成人aa在线观看| 黄色 视频免费看| 男人添女人高潮全过程视频| 大香蕉久久网| 欧美激情国产日韩精品一区| 高清黄色对白视频在线免费看| 国产极品天堂在线| 成人毛片60女人毛片免费| xxx大片免费视频| 精品少妇内射三级| 亚洲av免费高清在线观看| 丝瓜视频免费看黄片| 女性生殖器流出的白浆| 如日韩欧美国产精品一区二区三区| 午夜久久久在线观看| xxx大片免费视频| 久久热在线av| 一级黄片播放器| 十八禁网站网址无遮挡| 大码成人一级视频| 国产成人欧美| 男女下面插进去视频免费观看 | 热99久久久久精品小说推荐| 我的女老师完整版在线观看| 激情视频va一区二区三区| 岛国毛片在线播放| 中文字幕人妻丝袜制服| 三级国产精品片| 宅男免费午夜| 亚洲精品,欧美精品| 最新中文字幕久久久久| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品成人久久小说| 男人舔女人的私密视频| 99国产精品免费福利视频| 超色免费av| 欧美国产精品一级二级三级| 精品国产国语对白av| 亚洲精品国产av蜜桃| 成人免费观看视频高清| av播播在线观看一区| 日韩伦理黄色片| 综合色丁香网| 97人妻天天添夜夜摸| 飞空精品影院首页| 国产精品麻豆人妻色哟哟久久| 岛国毛片在线播放| 国产精品一国产av| 美女xxoo啪啪120秒动态图| 国产免费一区二区三区四区乱码| 国产精品蜜桃在线观看| 尾随美女入室| 亚洲一区二区三区欧美精品| 欧美另类一区| 国产精品99久久99久久久不卡 | 亚洲欧洲精品一区二区精品久久久 | 狠狠婷婷综合久久久久久88av| 欧美另类一区| 国产精品久久久av美女十八| 久久精品久久久久久噜噜老黄| 下体分泌物呈黄色| 国产白丝娇喘喷水9色精品| 日韩av免费高清视频| 欧美 日韩 精品 国产| 最近最新中文字幕大全免费视频 | 国产有黄有色有爽视频| 免费观看性生交大片5| 成人综合一区亚洲| 国产一区二区在线观看av| 黄色视频在线播放观看不卡| 老司机影院毛片| 久久久精品94久久精品| 欧美精品高潮呻吟av久久| 黄片无遮挡物在线观看| 国产精品久久久久成人av| 久久av网站| xxx大片免费视频| 女的被弄到高潮叫床怎么办| 亚洲av国产av综合av卡| 日本欧美国产在线视频| 成人亚洲欧美一区二区av| 成人18禁高潮啪啪吃奶动态图| 免费黄频网站在线观看国产| 精品国产乱码久久久久久小说| 国产精品久久久久久精品电影小说| 亚洲成av片中文字幕在线观看 | 欧美日韩视频高清一区二区三区二| 久久久国产欧美日韩av| 免费在线观看黄色视频的| 51国产日韩欧美| 国产av精品麻豆| 久久人人爽av亚洲精品天堂| 欧美变态另类bdsm刘玥| 亚洲综合色网址| 自线自在国产av| 男女边摸边吃奶| 99九九在线精品视频| 人人妻人人爽人人添夜夜欢视频| 精品久久国产蜜桃| 国产xxxxx性猛交| av福利片在线| 激情视频va一区二区三区| 久久99一区二区三区| av在线播放精品| 日韩精品免费视频一区二区三区 | 黄色 视频免费看| 交换朋友夫妻互换小说| 日韩免费高清中文字幕av| 人人妻人人爽人人添夜夜欢视频| 伦理电影免费视频| 日本欧美视频一区| 国产国语露脸激情在线看| 欧美xxxx性猛交bbbb| 在线天堂最新版资源| 国产av国产精品国产| 午夜免费男女啪啪视频观看| 中文精品一卡2卡3卡4更新| 中文字幕制服av| 日韩av不卡免费在线播放| 夜夜骑夜夜射夜夜干| 岛国毛片在线播放| 日本爱情动作片www.在线观看| 国产色婷婷99| av又黄又爽大尺度在线免费看| 熟妇人妻不卡中文字幕| 国产一区二区在线观看日韩| 久热久热在线精品观看| 日韩av在线免费看完整版不卡| 午夜老司机福利剧场| 99国产精品免费福利视频| 国产精品欧美亚洲77777| 久久久久久久国产电影| 国产av国产精品国产| 午夜福利乱码中文字幕| 777米奇影视久久| 国产在线一区二区三区精| 制服诱惑二区| 欧美丝袜亚洲另类| 亚洲精品,欧美精品| av女优亚洲男人天堂| 欧美最新免费一区二区三区| 久久国产亚洲av麻豆专区| 国产不卡av网站在线观看| 国产精品秋霞免费鲁丝片| 欧美成人午夜精品| 亚洲av国产av综合av卡| 少妇 在线观看| 亚洲一码二码三码区别大吗| 春色校园在线视频观看| 国产精品久久久久久久电影| 美女内射精品一级片tv| 黄色配什么色好看| 日日爽夜夜爽网站| 欧美激情极品国产一区二区三区 | 亚洲成人av在线免费| 亚洲 欧美一区二区三区| 尾随美女入室| 黄片播放在线免费| 久久久久久久久久成人| 一级a做视频免费观看| 国产精品麻豆人妻色哟哟久久| 搡女人真爽免费视频火全软件| 国产成人欧美| 国产熟女午夜一区二区三区| 国产69精品久久久久777片| 日韩av不卡免费在线播放| www日本在线高清视频| 久久久久久久久久人人人人人人| 久久99热这里只频精品6学生| 热99久久久久精品小说推荐| 狂野欧美激情性xxxx在线观看| 成年动漫av网址| 大码成人一级视频| 日韩视频在线欧美| 欧美精品国产亚洲| 日韩电影二区| 少妇的逼好多水| 在线观看一区二区三区激情| 一本色道久久久久久精品综合| 免费黄网站久久成人精品| av播播在线观看一区| 一级黄片播放器| 久久免费观看电影| 亚洲五月色婷婷综合| 久久久久久久久久久久大奶| 亚洲人与动物交配视频| 亚洲av国产av综合av卡| 一本久久精品| 亚洲精品一二三| 91aial.com中文字幕在线观看| 精品国产国语对白av| 国产毛片在线视频| 国产精品秋霞免费鲁丝片| 人妻人人澡人人爽人人| 亚洲欧美成人精品一区二区| 在线观看免费日韩欧美大片| av不卡在线播放| 欧美精品高潮呻吟av久久| 婷婷色麻豆天堂久久| av国产精品久久久久影院| 免费黄频网站在线观看国产| av在线观看视频网站免费| 男男h啪啪无遮挡| 秋霞伦理黄片| 一区在线观看完整版| 国内精品宾馆在线| 欧美激情 高清一区二区三区| 亚洲国产精品专区欧美| 免费黄色在线免费观看| 国产成人精品在线电影| 日本猛色少妇xxxxx猛交久久| 中文字幕亚洲精品专区| 寂寞人妻少妇视频99o| videossex国产| 久久国产精品大桥未久av| 亚洲婷婷狠狠爱综合网| 日本免费在线观看一区| 五月天丁香电影| 永久网站在线| 美女国产高潮福利片在线看| 欧美精品人与动牲交sv欧美| 亚洲欧美一区二区三区国产| 欧美亚洲 丝袜 人妻 在线| 亚洲精品第二区| 18在线观看网站| 51国产日韩欧美| 亚洲精品乱码久久久久久按摩| 国产黄色视频一区二区在线观看| 视频在线观看一区二区三区| 人妻系列 视频| 蜜臀久久99精品久久宅男| 麻豆精品久久久久久蜜桃| 亚洲精品久久久久久婷婷小说| 国产精品国产三级国产专区5o| 欧美亚洲 丝袜 人妻 在线| 天堂俺去俺来也www色官网| 国产精品国产三级专区第一集| 青春草国产在线视频| 久久99蜜桃精品久久| 亚洲av综合色区一区| 成人黄色视频免费在线看| 飞空精品影院首页| 亚洲国产av影院在线观看| 极品少妇高潮喷水抽搐| 欧美精品国产亚洲| 国产又色又爽无遮挡免| 久久久久久久精品精品| 国产精品一区二区在线观看99| 久久精品国产综合久久久 | 丝袜脚勾引网站| 日韩av在线免费看完整版不卡| 少妇熟女欧美另类| 如日韩欧美国产精品一区二区三区| 亚洲精品456在线播放app| 欧美日韩视频精品一区| 国产深夜福利视频在线观看| 黑人高潮一二区| 日韩精品免费视频一区二区三区 | 国产日韩欧美在线精品| 国产精品欧美亚洲77777| 国产又色又爽无遮挡免| 看免费av毛片| 久久女婷五月综合色啪小说| 日韩免费高清中文字幕av| 日本黄大片高清| 中文乱码字字幕精品一区二区三区| 伊人亚洲综合成人网| 最近2019中文字幕mv第一页| av在线观看视频网站免费| 亚洲精品456在线播放app| 国产欧美日韩综合在线一区二区| 91成人精品电影| 亚洲综合色网址| 赤兔流量卡办理| 久久久久人妻精品一区果冻| 免费大片黄手机在线观看| 久久韩国三级中文字幕| 国产有黄有色有爽视频| a 毛片基地| 亚洲精品视频女| 建设人人有责人人尽责人人享有的| 成人手机av| 亚洲国产欧美日韩在线播放| 七月丁香在线播放| 一区二区日韩欧美中文字幕 | 新久久久久国产一级毛片| 精品久久蜜臀av无| 午夜视频国产福利| 亚洲 欧美一区二区三区| 夫妻性生交免费视频一级片| 亚洲国产色片| 两性夫妻黄色片 | 欧美精品一区二区免费开放| 久久久久久久国产电影| 亚洲精品色激情综合| 考比视频在线观看| 一区二区三区四区激情视频| 黑丝袜美女国产一区| 国产老妇伦熟女老妇高清| 18在线观看网站| 女人被躁到高潮嗷嗷叫费观| 亚洲成国产人片在线观看| 两个人免费观看高清视频| 青春草亚洲视频在线观看| 看免费成人av毛片| 一本久久精品| 这个男人来自地球电影免费观看 | 少妇 在线观看| 一区在线观看完整版| 十八禁高潮呻吟视频| 国产国拍精品亚洲av在线观看| 欧美 亚洲 国产 日韩一| 国产色婷婷99| 一边亲一边摸免费视频| 高清视频免费观看一区二区| 热re99久久国产66热| 免费观看在线日韩| 国产乱来视频区| 日韩伦理黄色片| 老女人水多毛片| 日韩av免费高清视频| 色5月婷婷丁香| 国产免费视频播放在线视频| 久久人妻熟女aⅴ| 欧美xxxx性猛交bbbb| 免费观看a级毛片全部| videos熟女内射| 国产男女内射视频| 久久精品国产亚洲av涩爱| 午夜91福利影院| 国产永久视频网站| 国产欧美日韩一区二区三区在线| 桃花免费在线播放| 成人午夜精彩视频在线观看| 日韩成人av中文字幕在线观看| 国产一区二区三区综合在线观看 | 日韩人妻精品一区2区三区| 国产69精品久久久久777片| 日韩制服骚丝袜av| 国产av一区二区精品久久| 综合色丁香网| 狠狠婷婷综合久久久久久88av| 日韩 亚洲 欧美在线| 只有这里有精品99| 涩涩av久久男人的天堂| 搡老乐熟女国产| 不卡视频在线观看欧美| 精品少妇内射三级| 母亲3免费完整高清在线观看 | 婷婷成人精品国产| 免费不卡的大黄色大毛片视频在线观看| 午夜免费男女啪啪视频观看| 丝袜美足系列| 男的添女的下面高潮视频| 飞空精品影院首页| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利,免费看| 99久久综合免费| 欧美精品一区二区大全| 国产欧美日韩一区二区三区在线| 久久久国产一区二区| 久久精品国产亚洲av涩爱| 春色校园在线视频观看| 黄片无遮挡物在线观看| av在线老鸭窝| 视频中文字幕在线观看| 国产精品久久久久久久电影| 性色av一级| 最近最新中文字幕免费大全7| 亚洲欧美成人精品一区二区| 午夜老司机福利剧场| freevideosex欧美| 亚洲少妇的诱惑av| 成人漫画全彩无遮挡| 国产精品一区二区在线观看99| 日韩 亚洲 欧美在线| 久久久久精品久久久久真实原创| h视频一区二区三区| 亚洲精品av麻豆狂野| 纵有疾风起免费观看全集完整版| 国产精品久久久久久久久免| 高清av免费在线| 黄色一级大片看看| 97人妻天天添夜夜摸| 视频区图区小说| 日韩精品免费视频一区二区三区 | 国产成人一区二区在线| 国语对白做爰xxxⅹ性视频网站| 日产精品乱码卡一卡2卡三| 精品人妻一区二区三区麻豆| 丰满迷人的少妇在线观看| 插逼视频在线观看| 国产亚洲一区二区精品| 一级毛片 在线播放| 国产毛片在线视频| 日本欧美国产在线视频| 欧美精品一区二区大全| 熟女人妻精品中文字幕| 精品福利永久在线观看| 色婷婷久久久亚洲欧美| 国产免费一级a男人的天堂| 美国免费a级毛片| 国产有黄有色有爽视频| av女优亚洲男人天堂| 成人手机av| 国产乱人偷精品视频| 亚洲成国产人片在线观看| 日本欧美视频一区| 精品午夜福利在线看| 日本欧美视频一区| 国产片内射在线| 日韩欧美一区视频在线观看| 2018国产大陆天天弄谢| av视频免费观看在线观看| 男女高潮啪啪啪动态图| 久久av网站| 久久 成人 亚洲| 在线天堂最新版资源| 曰老女人黄片| 这个男人来自地球电影免费观看 | 大片电影免费在线观看免费| 秋霞伦理黄片| 精品第一国产精品| 国产在线视频一区二区| 18禁在线无遮挡免费观看视频| 九九在线视频观看精品| 男人爽女人下面视频在线观看| 亚洲激情五月婷婷啪啪| 免费看av在线观看网站| 国产欧美日韩综合在线一区二区| videos熟女内射| 男女边吃奶边做爰视频| 亚洲国产精品999| 国产精品久久久久久av不卡| 国产极品天堂在线| 亚洲第一av免费看| 久久韩国三级中文字幕| 蜜臀久久99精品久久宅男| 在线观看一区二区三区激情| 国产成人免费无遮挡视频| 久久免费观看电影| 老女人水多毛片| 亚洲,一卡二卡三卡| 伊人亚洲综合成人网| av在线app专区| 亚洲精品久久午夜乱码| 久久久亚洲精品成人影院| 国产亚洲av片在线观看秒播厂| 久久青草综合色| 视频在线观看一区二区三区| 国产一区有黄有色的免费视频| 日日爽夜夜爽网站| 蜜臀久久99精品久久宅男| 国产xxxxx性猛交| 好男人视频免费观看在线| 91aial.com中文字幕在线观看| 日韩电影二区| 午夜老司机福利剧场| 美女主播在线视频| 日韩在线高清观看一区二区三区| 欧美国产精品va在线观看不卡| 搡女人真爽免费视频火全软件| 在线观看人妻少妇| 中国国产av一级| 人妻 亚洲 视频| 国产毛片在线视频| 青青草视频在线视频观看| 日本黄色日本黄色录像|