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

    基于改進(jìn)TOPSIS法的鹽城市水資源承載力評(píng)價(jià)及障礙因子診斷

    2018-12-28 06:09:48孫雅茹董增川
    中國農(nóng)村水利水電 2018年12期
    關(guān)鍵詞:鹽城市子系統(tǒng)承載力

    孫雅茹,董增川,劉 淼

    (1. 河海大學(xué)水文水資源學(xué)院,南京 210098;2. 江蘇省水文水資源勘測(cè)局,南京 210029)

    水資源是人類賴以生存發(fā)展的最重要的基礎(chǔ)資源之一,但是隨著社會(huì)經(jīng)濟(jì)的迅速發(fā)展,水資源短缺、水環(huán)境惡化、水生態(tài)退化等水問題相繼出現(xiàn)[1],區(qū)域水資源承載力是水資源、社會(huì)、經(jīng)濟(jì)和生態(tài)環(huán)境的綜合反映,是區(qū)域水安全的重要度量[2],因此水資源承載力的研究對(duì)于水安全和區(qū)域的協(xié)調(diào)發(fā)展有著重要意義。國內(nèi)新疆軟課題組最早提出水資源承載力的研究[3],夏軍等運(yùn)用水資源承載力研究水安全[4],趙建世等運(yùn)用水資源承載能力計(jì)算模型分析海河流域的水資源承載能力演變[5],金菊良等建立了不同承載標(biāo)準(zhǔn)下區(qū)域水資源承載力評(píng)價(jià)模型[6]。

    到目前為止,水資源承載力評(píng)價(jià)的方法眾多,例如綜合指標(biāo)法[7]、模糊綜合評(píng)價(jià)法[8]、主成分分析法[9]、投影尋蹤法[10]等。綜合指標(biāo)法雖然簡單方便,但是對(duì)于問題的分析和認(rèn)知 不夠具體細(xì)致;模糊綜合評(píng)價(jià)法應(yīng)用廣泛,但是容易遺失信息而對(duì)評(píng)價(jià)結(jié)果造成影響;主成分分析法可以將多個(gè)指標(biāo)降維處理成幾個(gè)綜合變量,但是有時(shí)會(huì)出現(xiàn)該綜合變量物理意義不明確的情況;投影尋蹤法不需要事先確定各項(xiàng)指標(biāo)的權(quán)重,但是不適用于樣本較少的情況。本文采用改進(jìn)的Topsis法進(jìn)行水資源承載力的評(píng)價(jià),Topsis方法在土地利用評(píng)價(jià)[11]、生態(tài)安全評(píng)價(jià)[12]及電能質(zhì)量評(píng)價(jià)[13]等多方面有所應(yīng)用,該方法在指標(biāo)數(shù)量、樣本數(shù)量和數(shù)據(jù)等方面都沒有嚴(yán)格的限制和要求,原始數(shù)據(jù)利用充分,信息損失較少[11],而且該方法操作簡單方便。鑒于此,本文從水資源、社會(huì)、經(jīng)濟(jì)和生態(tài)環(huán)境4方面構(gòu)建水資源承載力評(píng)價(jià)指標(biāo)體系,以鹽城市為例,基于鹽城2005、2010、2015和2020四個(gè)水平年的指標(biāo)數(shù)據(jù),采用組合賦權(quán)法確定權(quán)重,以改進(jìn)的Topsis法和障礙因子診斷模型為理論基礎(chǔ),進(jìn)行鹽城市水資源承載力評(píng)價(jià)并識(shí)別出影響其的關(guān)鍵因素,為鹽城市水資源規(guī)劃和保護(hù)提供理論依據(jù)。

    1 研究方法

    1.1 指標(biāo)體系構(gòu)建

    水資源承載力是指在一定的社會(huì)背景及技術(shù)條件的支持下,以可預(yù)見期內(nèi)滿足合理的生態(tài)環(huán)境用水和保護(hù)生態(tài)環(huán)境為前提,考慮來水、工程情況和用水需求等因素,水資源對(duì)區(qū)域社會(huì)經(jīng)濟(jì)的最大負(fù)荷[14]。水資源承載力與社會(huì)、經(jīng)濟(jì)和生態(tài)環(huán)境息息相關(guān),彼此相互影響,相互作用,本研究參考相關(guān)文獻(xiàn)[14,15]和相關(guān)研究,從水資源、社會(huì)、經(jīng)濟(jì)和生態(tài)環(huán)境4方面考慮,選擇了15個(gè)指標(biāo)構(gòu)建水資源承載力評(píng)價(jià)指標(biāo)體系(表1),該指標(biāo)體系較為全面描述了水資源承載力的相關(guān)影響因素,從而對(duì)區(qū)域水資源承載力做出合理科學(xué)的評(píng)價(jià)。

    表1 水資源承載力評(píng)價(jià)指標(biāo)體系Tab.1 Evaluation index system for water resources carrying capacity

    1.2 權(quán)重確定

    本文采用基于博弈論的組合賦權(quán)法[16]來確定指標(biāo)體系中各指標(biāo)的權(quán)重,該方法可以綜合主觀與客觀權(quán)重相結(jié)合,提高了結(jié)果的可靠性和科學(xué)性。該方法的基本思想是在不同方法求得的權(quán)重之間需找妥協(xié),使得可能的權(quán)重與各權(quán)重之間的各自偏差最小化。本文采用層次分析法和熵值法計(jì)算各指標(biāo)的權(quán)重,并根據(jù)組合賦權(quán)法的思想,綜合兩種方法的結(jié)果求得最終權(quán)重。

    1.3 改進(jìn)TOPSIS模型

    Topsis法(逼近與理想值的排序方法)可以根據(jù)有限個(gè)數(shù)的評(píng)價(jià)方案與理想化目標(biāo)的接近程度來進(jìn)行排序,是一種適用于多指標(biāo)、多方案進(jìn)行比較選擇的分析方法。相比傳統(tǒng)的Topsis法,改進(jìn)的Topsis法考慮了各評(píng)價(jià)指標(biāo)的相對(duì)重要程度,構(gòu)造了加權(quán)的決策矩陣;傳統(tǒng)Topsis評(píng)價(jià)中運(yùn)用各計(jì)算方案的指標(biāo)數(shù)據(jù)進(jìn)行指標(biāo)標(biāo)準(zhǔn)化,也從標(biāo)準(zhǔn)化矩陣中選取最大值和最小值作為正、負(fù)理想值,這種方法沒有劃分指標(biāo)標(biāo)準(zhǔn),只能反映區(qū)域水資源承載力狀態(tài)變化的趨勢(shì)或者人為劃分最終評(píng)價(jià)等級(jí),在科學(xué)合理地反映區(qū)域的水資源承載力處于何等水平上有一定難度,同時(shí)當(dāng)評(píng)價(jià)對(duì)象自身?xiàng)l件發(fā)生變化時(shí),容易引起正、負(fù)理想解的變化,缺乏一定的穩(wěn)定性[17],因此本研究引入了指標(biāo)分級(jí)標(biāo)準(zhǔn)進(jìn)行改進(jìn),最終求出各方案與理想化目標(biāo)的接近程度,并判斷出各方案中水資源承載力所處的水平。計(jì)算步驟如下。

    (1)數(shù)據(jù)標(biāo)準(zhǔn)化處理。指標(biāo)分為成本型指標(biāo)(指標(biāo)值越大越好)和收益型指標(biāo)(指標(biāo)值越小越好),所選評(píng)價(jià)指標(biāo)因?yàn)榱烤V不同需要進(jìn)行標(biāo)準(zhǔn)化處理。假設(shè)指標(biāo)評(píng)價(jià)體系中一共有k個(gè)指標(biāo),s個(gè)待標(biāo)準(zhǔn)化方案,其中s=n+(m-1),n為待評(píng)價(jià)方案的個(gè)數(shù),m為指標(biāo)標(biāo)準(zhǔn)劃分的等級(jí)數(shù),即有(m-1)個(gè)等級(jí)之間的臨界值組成的方案。

    對(duì)于成本型指標(biāo):

    (1)

    對(duì)于收益型指標(biāo):

    (2)

    式中:xij和pij分別為j個(gè)標(biāo)準(zhǔn)化方案中第i個(gè)指標(biāo)的原始值和標(biāo)準(zhǔn)化值,當(dāng)j=1,…,n時(shí),xij和pij分別為待評(píng)價(jià)方案的指標(biāo)原始值和標(biāo)準(zhǔn)化值,當(dāng)j=n+1,…,n+(m-1)時(shí),xij和pij分別為各指標(biāo)等級(jí)標(biāo)準(zhǔn)劃分的臨界值及其標(biāo)準(zhǔn)化值;max(xij)和min(xij)分別為s個(gè)待標(biāo)準(zhǔn)化方案中第i個(gè)指標(biāo)的最大值和最小值。

    (2)構(gòu)建加權(quán)決策矩陣。由組合賦權(quán)法求得各指標(biāo)的權(quán)重為ωi,結(jié)合指標(biāo)標(biāo)準(zhǔn)化的結(jié)果,可以構(gòu)建加權(quán)決策矩陣V=(vij)k×s:

    vij=ωipij

    (3)

    (4)

    (5)

    (6)

    (7)

    (5)計(jì)算貼近度。貼近度Tj表征各方案與最優(yōu)方案的接近程度,其值越大,表明區(qū)域水資源承載狀態(tài)越好,計(jì)算公式如下:

    (8)

    通過該貼近度的計(jì)算,可以求出指標(biāo)標(biāo)準(zhǔn)(m-1)個(gè)臨界值方案的貼近度,從而可以判斷各評(píng)價(jià)方案的水資源承載力所屬的等級(jí)。

    1.4 障礙因子診斷模型

    為有效提高區(qū)域的水資源承載力,需要找出對(duì)區(qū)域水資源承載力產(chǎn)生負(fù)面影響的主要障礙因素。目前障礙度診斷模型的研究還處于初級(jí)階段,通常方法是采用因子貢獻(xiàn)度、指標(biāo)偏離度和障礙度來進(jìn)行分析診斷[11,12]。因子貢獻(xiàn)度一般指單一指標(biāo)對(duì)總目標(biāo)的貢獻(xiàn)大小,通??梢杂弥笜?biāo)權(quán)重ωi表示;指標(biāo)偏離度指指標(biāo)實(shí)際值與最優(yōu)目標(biāo)值之間的差距,可以用1與各指標(biāo)的標(biāo)準(zhǔn)化值pij之差來表示;障礙度Qij表示第j(j=1,2,…,n)個(gè)評(píng)價(jià)方案中各子系統(tǒng)或者各指標(biāo)對(duì)區(qū)域水資源承載力的影響程度,計(jì)算公式如下:

    (9)

    2 鹽城市水資源承載力評(píng)價(jià)

    鹽城市位于淮河流域的下游,蘇北平原中東部,其多年平均降水量為1 016.6 mm,降雨年內(nèi)分配不均,全年70%的雨量集中在汛期。鹽城市內(nèi)河網(wǎng)密布,溝渠縱橫,多年平均本地水資源量47.89 億m3,但是近年來鹽城市的水污染和水環(huán)境問題比較嚴(yán)重,水質(zhì)型缺水問題突出。鹽城市是全國重點(diǎn)糧、棉、油生產(chǎn)基地,是江蘇省沿海開發(fā)戰(zhàn)略的重點(diǎn)地區(qū),其社會(huì)經(jīng)濟(jì)發(fā)展還有進(jìn)步的空間,為緩解供需水之間的矛盾,營造良好的生態(tài)環(huán)境,需要了解鹽城市的水資源狀況,科學(xué)規(guī)劃和保護(hù)鹽城市水資源,因此有必要對(duì)鹽城市水資源承載力進(jìn)行評(píng)價(jià)及障礙因子診斷。

    2.1 權(quán)重確定

    根據(jù)上文構(gòu)建的鹽城市水資源承載力評(píng)價(jià)指標(biāo)體系和鹽城市的具體情況,先分別采用層次分析法和熵值法計(jì)算得到權(quán)重,然后采用組合賦權(quán)法求得綜合權(quán)重,結(jié)果見表2。

    2.2 改進(jìn)的Topsis水資源承載力評(píng)價(jià)模型

    指標(biāo)標(biāo)準(zhǔn)的等級(jí)劃分對(duì)評(píng)價(jià)結(jié)果有著重要的影響,根據(jù)國家相關(guān)規(guī)范和標(biāo)準(zhǔn),結(jié)合鹽城市的相關(guān)發(fā)展規(guī)劃以及水資源保護(hù)方面的要求,參考相關(guān)文獻(xiàn)和其他地區(qū)的標(biāo)準(zhǔn)[3,18,19],確定了評(píng)價(jià)體系各指標(biāo)的標(biāo)準(zhǔn),并將其劃分為 “好(Ⅰ)”、“較好(Ⅱ)”、“中等(Ⅲ)”、“較差(Ⅳ)”和“差(Ⅴ)”這5個(gè)等級(jí), 5個(gè)等級(jí)由臨界值1(Ⅰ和Ⅱ級(jí)之間)、臨界值2(Ⅱ和Ⅲ級(jí)之間)、臨界值3(Ⅲ和Ⅳ級(jí)之間)和臨界值4(Ⅳ和Ⅴ級(jí)之間)這4個(gè)臨界值劃分,詳見表3。

    表2 指標(biāo)權(quán)重Tab.2 Index weights

    表3 指標(biāo)標(biāo)準(zhǔn)Tab.3 Standards of indexes

    本研究針對(duì)鹽城市2005、2010、2015和2020年這4個(gè)水平年開展水資源承載力評(píng)價(jià),根據(jù)《鹽城市水資源公報(bào)》、《鹽城市統(tǒng)計(jì)年鑒》《鹽城市水資源保護(hù)規(guī)劃》、《鹽城市水利現(xiàn)代化規(guī)劃》、《鹽城市節(jié)水型社會(huì)規(guī)劃》等相關(guān)資料,可以統(tǒng)計(jì)出各指標(biāo)在2005、2010和2015年的數(shù)值,并預(yù)測(cè)出各指標(biāo)在2020水平年的數(shù)值,詳見表4。

    根據(jù)式(1)、(2)將各指標(biāo)的原始值和標(biāo)準(zhǔn)的臨界值進(jìn)行標(biāo)準(zhǔn)化后,結(jié)合權(quán)重計(jì)算生成加權(quán)決策矩陣并根據(jù)公式(4)、(5)確定各指標(biāo)的正、負(fù)理想解,最后根據(jù)式(6)~(8)計(jì)算歐氏距離和貼近度,結(jié)果見表5。

    表4 不同年份指標(biāo)值Tab.4 Index values in different years

    由表5、6可知,鹽城市水資源承載力在2005年處于Ⅳ級(jí)的水平,即鹽城市在2005年水資源承載力處于較差的水平;在2010年和2015年,鹽城市水資源承載力處于Ⅲ級(jí),即中等水平;到2020年,鹽城市水資源承載力將提升至Ⅱ級(jí),即較好的水平。

    根據(jù)在不同水平年的貼近度計(jì)算結(jié)果,2005-2020年,貼近度數(shù)值在不斷上升,說明鹽城市水資源承載力的狀態(tài)在不斷好轉(zhuǎn)提高。從水資源子系統(tǒng)來看,2005-2020年,總體上供水壓力呈現(xiàn)小幅度的增長,由于人口變化并不大,所以人均水資源占有量變化也不大,但是在地下水開采方面,由于《鹽城市地下水壓采方案》等相關(guān)政策的推行,鹽城市的地下水超采情況有明顯的好轉(zhuǎn);從社會(huì)子系統(tǒng)來看,生活方面的用水效率和污水處理率都有了顯著的提升,但是還有一定的提升空間;從經(jīng)濟(jì)子系統(tǒng)來看,隨著人民生活水平的提高,農(nóng)業(yè)和工業(yè)方面的用水效率從較差的水平提高到了較好的水平,使得萬元GDP用水量近年來明顯降低,經(jīng)濟(jì)子系統(tǒng)的狀態(tài)明顯改善;從生態(tài)環(huán)境子系統(tǒng)來看,雖然整體生態(tài)環(huán)境子系統(tǒng)的狀態(tài)在好轉(zhuǎn),但是鹽城市的水環(huán)境和水污染問題一直比較嚴(yán)重,2005-2015年水功能區(qū)水質(zhì)達(dá)標(biāo)率一直低于70%,嚴(yán)重制約了鹽城市的水資源承載力。

    表5 鹽城市水資源承載力評(píng)價(jià)結(jié)果Tab.5 Evaluation results of water resources carrying capacity of yancheng city

    表6 鹽城市水資源承載力評(píng)價(jià)等級(jí)標(biāo)準(zhǔn)Tab.6 Evaluation standards of water resources carrying capacity of yancheng city

    2.3 障礙因子診斷

    根據(jù)式(9)計(jì)算各指標(biāo)以及各子系統(tǒng)的障礙度,基于障礙度排序,找出影響鹽城市水資源承載力的主要障礙因子,本文列出排名前6的指標(biāo)因子,詳見表7,各子系統(tǒng)障礙度變化趨勢(shì)見圖1。

    表7 鹽城市水資源承載力前6位障礙因子Tab.7 Top 6 obstacle factors of water resources carrying capacity of yancheng city

    圖1 各子系統(tǒng)障礙度Fig.1 Obstacle degrees of sub-systems

    從表7來看,2005-2020年,水資源子系統(tǒng)的障礙度隨著時(shí)間明顯在增加,從2005年的19.96%增加到2020年的49.08%,成為影響鹽城市水資源承載力的首要因素;而經(jīng)濟(jì)子系統(tǒng)的障礙度呈現(xiàn)相反的趨勢(shì),由于用水效率和污水處理方面的提高,使得經(jīng)濟(jì)子系統(tǒng)的障礙度明顯降低;社會(huì)子系統(tǒng)的障礙度也呈現(xiàn)輕微下降的趨勢(shì),該子系統(tǒng)2005-2020年障礙度一直小于其他子系統(tǒng);生態(tài)環(huán)境子系統(tǒng)的障礙度在這段時(shí)間內(nèi)變化不大,在2020年位居第二。

    根據(jù)表7,在2005年,排名前6的障礙因子之間障礙度的差距并不大,主要阻礙鹽城水水資源承載力的障礙因子集中在水質(zhì)、用水效率和水量方面,其中水功能區(qū)水質(zhì)達(dá)標(biāo)率(C12)是影響水資源承載力的首要障礙因子,其次是用水效率方面的供水管網(wǎng)扣損率(C5)和萬元GDP用水量(C8)兩個(gè)指標(biāo),接著是水量方面的人均水資源量(C1)和地下水超采面積比(C2)兩個(gè)指標(biāo);到了2010年,排名前6的障礙因子之間障礙度的差距依然不大,由于水質(zhì)情況有一定的改善,水功能區(qū)水質(zhì)達(dá)標(biāo)率(C12)的障礙度降至第二,由于用水效率的提高,相關(guān)指標(biāo)的障礙度排名均有所降低,而人均水資源量(C1)成為首要障礙因子;在2015年,排名前6的障礙因子之間障礙度的逐漸拉大,首要障礙因子和次要障礙因子沒有變化,由于各方面用水效率進(jìn)一步提高,其障礙度進(jìn)一步減少,地下水超采面積比(C2)和城鎮(zhèn)綠化覆蓋率(C15)的障礙度排名上升;預(yù)測(cè)到2020年,排名前6的障礙因子之間障礙度的差距將進(jìn)一步拉大,人均水資源量(C1)和地下水超采面積比(C2)兩個(gè)指標(biāo)依然位居前三,根據(jù)相關(guān)規(guī)劃指標(biāo)水功能區(qū)水質(zhì)達(dá)標(biāo)率(C12)將改善至85%,故在2020年該指標(biāo)障礙度降低,城鎮(zhèn)綠化覆蓋率(C15)的障礙度排名上升,用水效率方面只有灌溉水利用系數(shù)(C9)位于前六,其他用水效率相關(guān)指標(biāo)障礙度依然在不斷減小。

    結(jié)合鹽城市2005、2010、2015和2020年的水資源承載力的評(píng)價(jià)結(jié)果和障礙因子診斷結(jié)果,目前鹽城市的水污染和水環(huán)境問題比較嚴(yán)重,需要積極整治黑臭水體,加強(qiáng)排污監(jiān)管,加大污水處理設(shè)施的建設(shè),嚴(yán)格控制農(nóng)業(yè)面源污染,實(shí)現(xiàn)工業(yè)廢水全部達(dá)標(biāo)排放;同時(shí),一方面需要積極推廣節(jié)水型社會(huì)的建設(shè),尤其鹽城市農(nóng)業(yè)發(fā)達(dá),在農(nóng)業(yè)節(jié)水方面需要加大投入,加快農(nóng)業(yè)節(jié)水工程建設(shè),另一方面需要合理科學(xué)建設(shè)跨區(qū)域調(diào)水工程,緩解鹽城市的供水壓力。

    3 結(jié) 語

    (1)本文采用改進(jìn)的Topsis法構(gòu)建鹽城市水資源承載力評(píng)價(jià)模型,結(jié)合組合賦權(quán)法確定的權(quán)重,對(duì)鹽城市進(jìn)行了水資源承載力評(píng)價(jià),同時(shí)構(gòu)建障礙因子診斷模型,識(shí)別出影響鹽城市水資源承載力的障礙因子。該方法科學(xué)簡便,對(duì)區(qū)域水資源承載力評(píng)價(jià)和診斷有一定的參考價(jià)值。

    (2)根據(jù)鹽城市水資源承載力評(píng)價(jià)結(jié)果,2005-2015年,鹽城市水資源承載狀態(tài)呈現(xiàn)好轉(zhuǎn)的趨勢(shì),從2005年較差的水平提升到2015年中等水平,期間鹽城市大力推廣節(jié)水措施,用水效率迅速提高,但是水環(huán)境問題依然嚴(yán)重,預(yù)計(jì)到2020年,鹽城市水資源承載力將提高到較好的水平。

    (3)根據(jù)障礙因子診斷模型計(jì)算結(jié)果,從準(zhǔn)則層來看,在2005年時(shí),經(jīng)濟(jì)子系統(tǒng)和生態(tài)環(huán)境子系統(tǒng)對(duì)鹽城市水資源承載力影響較大,到2020年時(shí),水資源子系統(tǒng)和生態(tài)環(huán)境子系統(tǒng)的影響較大;從指標(biāo)層來看,用水效率的相關(guān)指標(biāo)的障礙度在不停降低,人均水資源量指標(biāo)由主要障礙因子上升為首要障礙因子,地下水超采面積比和城鎮(zhèn)綠化覆蓋率由一般障礙因子變成主要障礙因子,通過水污染治理等措施,水功能區(qū)水質(zhì)達(dá)標(biāo)率到2020年可由主要障礙因子降為一般障礙因子。

    猜你喜歡
    鹽城市子系統(tǒng)承載力
    鹽城市深化專項(xiàng)整治
    鹽城市扛起使命擔(dān)當(dāng) 鍛造應(yīng)急鐵軍
    不對(duì)中轉(zhuǎn)子系統(tǒng)耦合動(dòng)力學(xué)特性研究
    GSM-R基站子系統(tǒng)同步方案研究
    駝峰測(cè)長設(shè)備在線監(jiān)測(cè)子系統(tǒng)的設(shè)計(jì)與應(yīng)用
    CFRP-PCP板加固混凝土梁的抗彎承載力研究
    耐火鋼圓鋼管混凝土柱耐火極限和承載力
    潛艇極限承載力計(jì)算與分析
    鹽城市
    江蘇年鑒(2014年0期)2014-03-11 17:10:11
    車載ATP子系統(tǒng)緊急制動(dòng)限制速度計(jì)算
    成人亚洲精品一区在线观看| 国产av国产精品国产| 青春草视频在线免费观看| 亚洲欧美中文字幕日韩二区| 免费av中文字幕在线| 日韩制服骚丝袜av| 国产国语露脸激情在线看| 日韩不卡一区二区三区视频在线| 国产高清不卡午夜福利| 少妇被粗大的猛进出69影院| h视频一区二区三区| 国产精品一二三区在线看| 久久久久久久精品精品| 免费av中文字幕在线| 男女边吃奶边做爰视频| 狂野欧美激情性xxxx| 国产免费福利视频在线观看| 无遮挡黄片免费观看| 激情视频va一区二区三区| 在线观看人妻少妇| 日韩欧美精品免费久久| 最近中文字幕2019免费版| bbb黄色大片| av线在线观看网站| 大陆偷拍与自拍| 亚洲情色 制服丝袜| 一区在线观看完整版| 美女高潮到喷水免费观看| 亚洲视频免费观看视频| 在线观看一区二区三区激情| 国产精品一区二区精品视频观看| 国产极品天堂在线| 伦理电影大哥的女人| 国产精品久久久久久精品古装| 人人澡人人妻人| 久久久久久久久免费视频了| 2018国产大陆天天弄谢| 国产男人的电影天堂91| 男女之事视频高清在线观看 | 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩福利视频一区二区| 你懂的网址亚洲精品在线观看| 91成人精品电影| 一级a爱视频在线免费观看| 综合色丁香网| 性少妇av在线| 婷婷色麻豆天堂久久| 亚洲国产日韩一区二区| 久久影院123| 亚洲精品一二三| 亚洲精品日本国产第一区| 欧美日韩成人在线一区二区| 啦啦啦 在线观看视频| 欧美老熟妇乱子伦牲交| 我要看黄色一级片免费的| 黑人欧美特级aaaaaa片| 亚洲国产av影院在线观看| 久久久久国产一级毛片高清牌| 亚洲av在线观看美女高潮| 一区二区三区乱码不卡18| 9色porny在线观看| 久久午夜综合久久蜜桃| 成人影院久久| 免费不卡黄色视频| 1024视频免费在线观看| 51午夜福利影视在线观看| 国产精品 国内视频| 日韩av在线免费看完整版不卡| 人妻 亚洲 视频| 亚洲一区中文字幕在线| 国产激情久久老熟女| 国产成人欧美在线观看 | 精品一区在线观看国产| 亚洲av成人精品一二三区| 亚洲国产av影院在线观看| 女人久久www免费人成看片| 欧美人与善性xxx| 亚洲国产日韩一区二区| 久久国产精品大桥未久av| 18禁观看日本| 97精品久久久久久久久久精品| 一本久久精品| 中文字幕av电影在线播放| 丰满少妇做爰视频| 一二三四在线观看免费中文在| 欧美日韩亚洲高清精品| www.自偷自拍.com| 国产在视频线精品| 国产精品一国产av| 九草在线视频观看| 国产精品av久久久久免费| 亚洲精品视频女| 天天添夜夜摸| 一级爰片在线观看| 婷婷成人精品国产| 黄色视频在线播放观看不卡| 亚洲精品国产av成人精品| 91精品国产国语对白视频| 国产成人av激情在线播放| www.精华液| 观看av在线不卡| 国产成人av激情在线播放| 久久久久人妻精品一区果冻| 中文精品一卡2卡3卡4更新| 成年女人毛片免费观看观看9 | 久久精品久久久久久噜噜老黄| 纵有疾风起免费观看全集完整版| 欧美激情极品国产一区二区三区| av视频免费观看在线观看| 免费黄频网站在线观看国产| 成年av动漫网址| 免费看av在线观看网站| a 毛片基地| 我要看黄色一级片免费的| 亚洲婷婷狠狠爱综合网| 我的亚洲天堂| 国产高清国产精品国产三级| 国产男人的电影天堂91| 国产97色在线日韩免费| 免费在线观看黄色视频的| 国产熟女午夜一区二区三区| 国产探花极品一区二区| 青春草视频在线免费观看| 欧美变态另类bdsm刘玥| 深夜精品福利| 精品国产乱码久久久久久男人| 午夜免费观看性视频| 热re99久久精品国产66热6| 亚洲国产精品成人久久小说| 久热爱精品视频在线9| 天堂俺去俺来也www色官网| 国产黄色视频一区二区在线观看| 曰老女人黄片| 一级黄片播放器| 乱人伦中国视频| 日韩电影二区| www.自偷自拍.com| 免费观看性生交大片5| 美女大奶头黄色视频| 在线天堂最新版资源| 亚洲精品日韩在线中文字幕| 日本午夜av视频| 少妇的丰满在线观看| 香蕉国产在线看| 制服丝袜香蕉在线| 亚洲欧美一区二区三区黑人| 精品福利永久在线观看| 男的添女的下面高潮视频| 国产精品秋霞免费鲁丝片| 久久久精品94久久精品| 中文字幕高清在线视频| 99久久综合免费| 欧美国产精品va在线观看不卡| 秋霞在线观看毛片| 美女中出高潮动态图| 亚洲成人一二三区av| 婷婷成人精品国产| 久久狼人影院| 国产成人精品福利久久| 这个男人来自地球电影免费观看 | 亚洲欧美精品综合一区二区三区| 一本久久精品| 多毛熟女@视频| 人人妻人人添人人爽欧美一区卜| 午夜激情av网站| 伊人久久大香线蕉亚洲五| 欧美黑人精品巨大| 如何舔出高潮| 熟妇人妻不卡中文字幕| 三上悠亚av全集在线观看| 久久人人97超碰香蕉20202| 精品亚洲乱码少妇综合久久| 国产一区有黄有色的免费视频| 亚洲国产看品久久| 亚洲视频免费观看视频| 男女无遮挡免费网站观看| 亚洲第一青青草原| 国产熟女欧美一区二区| 精品国产乱码久久久久久男人| 亚洲成人av在线免费| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕另类日韩欧美亚洲嫩草| 精品免费久久久久久久清纯 | 丝袜美足系列| 国语对白做爰xxxⅹ性视频网站| 2021少妇久久久久久久久久久| 在现免费观看毛片| 国产精品免费视频内射| 最近2019中文字幕mv第一页| 少妇 在线观看| 国产无遮挡羞羞视频在线观看| 欧美亚洲日本最大视频资源| 女的被弄到高潮叫床怎么办| 国产日韩欧美亚洲二区| 亚洲伊人色综图| 久久人妻熟女aⅴ| 这个男人来自地球电影免费观看 | 日韩制服丝袜自拍偷拍| 秋霞伦理黄片| 女性生殖器流出的白浆| 欧美日韩亚洲综合一区二区三区_| 五月开心婷婷网| 男女下面插进去视频免费观看| 只有这里有精品99| 日韩中文字幕视频在线看片| 搡老乐熟女国产| 色婷婷久久久亚洲欧美| 国产野战对白在线观看| 秋霞伦理黄片| 最近中文字幕2019免费版| 亚洲美女黄色视频免费看| 啦啦啦在线观看免费高清www| 综合色丁香网| 老司机在亚洲福利影院| 超碰97精品在线观看| 99热全是精品| 亚洲三区欧美一区| 女人高潮潮喷娇喘18禁视频| kizo精华| 欧美精品一区二区大全| 极品人妻少妇av视频| 国产免费视频播放在线视频| 亚洲成人手机| 美国免费a级毛片| 免费黄色在线免费观看| 亚洲国产精品国产精品| 一本一本久久a久久精品综合妖精| 亚洲欧美精品自产自拍| 国产一区二区在线观看av| 国产99久久九九免费精品| 青春草国产在线视频| 久久毛片免费看一区二区三区| 久久久国产精品麻豆| 亚洲av日韩在线播放| 精品一区二区三区四区五区乱码 | 精品人妻熟女毛片av久久网站| 香蕉丝袜av| bbb黄色大片| 国产av国产精品国产| 亚洲成人一二三区av| 国产免费又黄又爽又色| 国产欧美日韩一区二区三区在线| 女人被躁到高潮嗷嗷叫费观| 久久av网站| 亚洲欧美精品自产自拍| 国产成人91sexporn| 大片免费播放器 马上看| 欧美国产精品va在线观看不卡| 狂野欧美激情性bbbbbb| 热re99久久精品国产66热6| 国产日韩欧美亚洲二区| 欧美老熟妇乱子伦牲交| 中文字幕人妻丝袜制服| 如日韩欧美国产精品一区二区三区| 久久久久视频综合| 亚洲 欧美一区二区三区| 欧美 亚洲 国产 日韩一| 午夜精品国产一区二区电影| 国产毛片在线视频| 大香蕉久久成人网| 一区二区av电影网| 午夜激情久久久久久久| 一区二区三区四区激情视频| 亚洲四区av| 欧美 日韩 精品 国产| 在线观看一区二区三区激情| 中文字幕人妻丝袜一区二区 | 久久热在线av| 国产精品久久久久成人av| 国产精品国产三级国产专区5o| 国产野战对白在线观看| 亚洲男人天堂网一区| 51午夜福利影视在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 七月丁香在线播放| 我的亚洲天堂| 9热在线视频观看99| 亚洲图色成人| 亚洲,欧美精品.| 天堂中文最新版在线下载| 如日韩欧美国产精品一区二区三区| 日韩中文字幕视频在线看片| 亚洲 欧美一区二区三区| 水蜜桃什么品种好| 97在线人人人人妻| 亚洲,欧美,日韩| 51午夜福利影视在线观看| 国产成人精品久久二区二区91 | 国产精品女同一区二区软件| 秋霞在线观看毛片| 午夜福利乱码中文字幕| 精品少妇黑人巨大在线播放| 不卡视频在线观看欧美| 免费高清在线观看日韩| 又大又黄又爽视频免费| xxxhd国产人妻xxx| 婷婷色av中文字幕| 日韩,欧美,国产一区二区三区| 亚洲精品aⅴ在线观看| 日本一区二区免费在线视频| 一区二区三区四区激情视频| 天天躁夜夜躁狠狠躁躁| 色94色欧美一区二区| 极品人妻少妇av视频| 午夜免费鲁丝| 99久国产av精品国产电影| 亚洲四区av| 亚洲图色成人| 亚洲美女视频黄频| 久久久久国产一级毛片高清牌| 亚洲精品国产av成人精品| 国产日韩欧美亚洲二区| 国产高清不卡午夜福利| 女人精品久久久久毛片| 王馨瑶露胸无遮挡在线观看| 如何舔出高潮| 十八禁网站网址无遮挡| 这个男人来自地球电影免费观看 | 老汉色∧v一级毛片| 国产成人午夜福利电影在线观看| 国产成人精品久久二区二区91 | 亚洲一区中文字幕在线| 久久精品aⅴ一区二区三区四区| 亚洲精品自拍成人| 熟女少妇亚洲综合色aaa.| 在线 av 中文字幕| tube8黄色片| 久久久亚洲精品成人影院| 国产亚洲最大av| 久久久久国产精品人妻一区二区| 久久毛片免费看一区二区三区| 亚洲精品视频女| 观看美女的网站| 如何舔出高潮| 大陆偷拍与自拍| 丰满乱子伦码专区| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品在线电影| 90打野战视频偷拍视频| 18在线观看网站| 丰满乱子伦码专区| 国产黄色视频一区二区在线观看| h视频一区二区三区| 国产一区二区三区综合在线观看| 亚洲精品久久久久久婷婷小说| 多毛熟女@视频| 国产无遮挡羞羞视频在线观看| 久久久国产精品麻豆| 少妇被粗大猛烈的视频| 久久97久久精品| 亚洲精品中文字幕在线视频| 国产又色又爽无遮挡免| a级毛片在线看网站| 国产成人啪精品午夜网站| 男人爽女人下面视频在线观看| 搡老岳熟女国产| 欧美日韩综合久久久久久| 欧美日韩一级在线毛片| 亚洲 欧美一区二区三区| 婷婷色av中文字幕| 国产精品二区激情视频| 青草久久国产| 色视频在线一区二区三区| 久久久精品区二区三区| 日本黄色日本黄色录像| 观看av在线不卡| 中文字幕av电影在线播放| 十八禁高潮呻吟视频| 欧美精品一区二区大全| 亚洲av电影在线观看一区二区三区| 国产成人欧美| 国产av精品麻豆| 免费在线观看黄色视频的| 一区二区三区四区激情视频| 国产精品一二三区在线看| 欧美人与善性xxx| 七月丁香在线播放| 精品少妇一区二区三区视频日本电影 | 狂野欧美激情性xxxx| 欧美人与性动交α欧美软件| 精品亚洲成a人片在线观看| 最近最新中文字幕大全免费视频 | av一本久久久久| 亚洲国产精品国产精品| 视频在线观看一区二区三区| 一二三四中文在线观看免费高清| 丝瓜视频免费看黄片| 国产精品av久久久久免费| 哪个播放器可以免费观看大片| 青青草视频在线视频观看| 超碰97精品在线观看| 国产在线视频一区二区| 成人国语在线视频| 久热爱精品视频在线9| www.av在线官网国产| 亚洲av福利一区| 欧美日韩亚洲高清精品| 久久鲁丝午夜福利片| 亚洲av成人精品一二三区| 久久精品国产综合久久久| 久久精品熟女亚洲av麻豆精品| 亚洲精品视频女| 亚洲,一卡二卡三卡| 伊人久久国产一区二区| 美女大奶头黄色视频| 中文字幕av电影在线播放| 精品一区二区三卡| 美女主播在线视频| a 毛片基地| 午夜福利在线免费观看网站| 大码成人一级视频| 欧美日韩av久久| 伦理电影大哥的女人| 夫妻午夜视频| 国产在线一区二区三区精| 亚洲七黄色美女视频| 国产一区二区 视频在线| 18禁观看日本| 国产免费福利视频在线观看| 国产成人av激情在线播放| 日本欧美视频一区| 免费在线观看黄色视频的| 欧美日韩精品网址| 性高湖久久久久久久久免费观看| 色吧在线观看| 午夜激情av网站| 热99久久久久精品小说推荐| 国产极品粉嫩免费观看在线| 国产精品久久久久久精品古装| 国产免费福利视频在线观看| 高清不卡的av网站| 欧美在线黄色| 久久这里只有精品19| 国产日韩欧美视频二区| 国产伦理片在线播放av一区| 中文字幕色久视频| 丝袜在线中文字幕| 国产有黄有色有爽视频| 久久久国产一区二区| 国产一级毛片在线| 亚洲成人av在线免费| 黑人猛操日本美女一级片| 美女视频免费永久观看网站| 狠狠婷婷综合久久久久久88av| av一本久久久久| 韩国精品一区二区三区| 欧美精品人与动牲交sv欧美| 久久99热这里只频精品6学生| 午夜福利视频在线观看免费| 亚洲av国产av综合av卡| 黄色毛片三级朝国网站| 日韩电影二区| 精品一区二区三区四区五区乱码 | 男人添女人高潮全过程视频| 制服丝袜香蕉在线| 久久热在线av| 午夜影院在线不卡| 在线观看免费午夜福利视频| 欧美国产精品va在线观看不卡| 国产一区二区三区综合在线观看| 免费久久久久久久精品成人欧美视频| 国产爽快片一区二区三区| 日韩av在线免费看完整版不卡| 2021少妇久久久久久久久久久| 啦啦啦 在线观看视频| 成人国产av品久久久| 人人妻人人添人人爽欧美一区卜| 美女午夜性视频免费| 亚洲美女视频黄频| 亚洲国产欧美日韩在线播放| av在线观看视频网站免费| 在线 av 中文字幕| 性高湖久久久久久久久免费观看| 欧美国产精品一级二级三级| 精品亚洲成a人片在线观看| 看非洲黑人一级黄片| 日韩中文字幕欧美一区二区 | av网站免费在线观看视频| 一级毛片我不卡| 成人亚洲精品一区在线观看| 亚洲欧洲精品一区二区精品久久久 | 一边摸一边抽搐一进一出视频| 亚洲,欧美,日韩| 大码成人一级视频| www.自偷自拍.com| 欧美日韩视频精品一区| 国产精品一二三区在线看| 成人毛片60女人毛片免费| 色播在线永久视频| 91精品三级在线观看| 国产日韩一区二区三区精品不卡| 国产一区有黄有色的免费视频| av免费观看日本| 亚洲欧美一区二区三区国产| 精品国产露脸久久av麻豆| 人人妻人人添人人爽欧美一区卜| 美女午夜性视频免费| 国产成人免费观看mmmm| 18禁动态无遮挡网站| 天天操日日干夜夜撸| 亚洲欧美一区二区三区久久| 日本黄色日本黄色录像| 精品第一国产精品| 考比视频在线观看| 亚洲国产中文字幕在线视频| 国产精品久久久av美女十八| 少妇被粗大猛烈的视频| 国产亚洲av高清不卡| 成人毛片60女人毛片免费| 亚洲精品国产av成人精品| 久久女婷五月综合色啪小说| 一二三四中文在线观看免费高清| 蜜桃在线观看..| 日韩大片免费观看网站| 黑人巨大精品欧美一区二区蜜桃| 日本午夜av视频| 亚洲精品在线美女| 国产极品天堂在线| 一二三四在线观看免费中文在| 男女下面插进去视频免费观看| 又黄又粗又硬又大视频| 人体艺术视频欧美日本| 91精品国产国语对白视频| 高清视频免费观看一区二区| 欧美另类一区| 国产又爽黄色视频| 国产麻豆69| 嫩草影视91久久| 日韩精品有码人妻一区| 老司机靠b影院| 大香蕉久久成人网| 麻豆乱淫一区二区| 视频在线观看一区二区三区| 国产精品嫩草影院av在线观看| 国产精品一区二区在线观看99| av国产精品久久久久影院| 久久天堂一区二区三区四区| 日本猛色少妇xxxxx猛交久久| 精品亚洲成a人片在线观看| 又大又黄又爽视频免费| 999精品在线视频| 亚洲精品,欧美精品| 少妇人妻精品综合一区二区| 人成视频在线观看免费观看| 欧美国产精品va在线观看不卡| 国产精品av久久久久免费| 热re99久久国产66热| 亚洲精品美女久久av网站| 欧美久久黑人一区二区| 另类亚洲欧美激情| 夫妻午夜视频| 成人午夜精彩视频在线观看| 久久久精品94久久精品| 国产av一区二区精品久久| 男女高潮啪啪啪动态图| 日本欧美视频一区| 七月丁香在线播放| 亚洲专区中文字幕在线 | 日韩免费高清中文字幕av| 人人妻人人澡人人爽人人夜夜| 黄片播放在线免费| 熟女av电影| 国产片内射在线| 国产av一区二区精品久久| 亚洲欧洲日产国产| 搡老岳熟女国产| 一级片'在线观看视频| 精品国产一区二区三区四区第35| 精品亚洲成a人片在线观看| 中文欧美无线码| 老司机影院成人| 欧美xxⅹ黑人| 亚洲欧美中文字幕日韩二区| 又大又爽又粗| 侵犯人妻中文字幕一二三四区| 一级,二级,三级黄色视频| 国产精品久久久久成人av| 国产免费福利视频在线观看| a级毛片黄视频| 亚洲一区二区三区欧美精品| 国产乱来视频区| 中文字幕人妻熟女乱码| 久久久久久久国产电影| 性少妇av在线| 搡老岳熟女国产| 多毛熟女@视频| 黄片小视频在线播放| 黄片播放在线免费| 亚洲在久久综合| 看免费av毛片| 亚洲在久久综合| 在线观看免费高清a一片| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 国产成人精品久久二区二区91 | 精品亚洲成国产av| 一本—道久久a久久精品蜜桃钙片| 男女午夜视频在线观看| 各种免费的搞黄视频| 男人添女人高潮全过程视频| 一边亲一边摸免费视频| 国产精品香港三级国产av潘金莲 | 亚洲成国产人片在线观看| 久久狼人影院| 伦理电影免费视频| 十八禁网站网址无遮挡| 99久久综合免费| 国产福利在线免费观看视频| 精品少妇一区二区三区视频日本电影 | 国产成人精品无人区| 欧美日韩av久久| 免费高清在线观看视频在线观看| 婷婷色av中文字幕| 国产国语露脸激情在线看| 国产亚洲av片在线观看秒播厂|