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

    基于二級(jí)修正的區(qū)域水資源需求量預(yù)測(cè)模型

    2010-05-12 07:31:24吳鳳平陳艷萍
    水資源保護(hù) 2010年3期
    關(guān)鍵詞:需水量需求量城區(qū)

    吳 丹,吳鳳平,陳艷萍

    (河海大學(xué)商學(xué)院,江蘇南京 210098)

    基于二級(jí)修正的區(qū)域水資源需求量預(yù)測(cè)模型

    吳 丹,吳鳳平,陳艷萍

    (河海大學(xué)商學(xué)院,江蘇南京 210098)

    根據(jù)區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展的用水現(xiàn)狀以及水資源需求的總體變化趨勢(shì),將Logistic生長(zhǎng)曲線模型和GM(1,1)等維新息模型相結(jié)合,建立組合預(yù)測(cè)模型,組合預(yù)測(cè)區(qū)域水資源需求量,并利用ARIMA模型對(duì)區(qū)域水資源需求量的組合預(yù)測(cè)結(jié)果進(jìn)行一級(jí)修正;在此基礎(chǔ)上,根據(jù)區(qū)域節(jié)水規(guī)劃的總體目標(biāo),剖析區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展綜合節(jié)水潛力,對(duì)區(qū)域水資源需求量的一級(jí)修正結(jié)果進(jìn)行二級(jí)修正,最終得到區(qū)域水資源需求量的綜合預(yù)測(cè)結(jié)果。算例分析結(jié)果表明,模型具有較高的預(yù)測(cè)精度。

    水資源需求量;Logistic生長(zhǎng)曲線模型;GM(1,1)等維新息模型;ARIMA模型;二級(jí)修正

    隨著區(qū)域社會(huì)經(jīng)濟(jì)的快速發(fā)展,水資源需求量不斷增長(zhǎng),水資源的稀缺性已嚴(yán)重制約區(qū)域社會(huì)、經(jīng)濟(jì)以及生態(tài)環(huán)境之間的協(xié)調(diào)持續(xù)發(fā)展。區(qū)域水資源需求量預(yù)測(cè)是初始水權(quán)配置以及水資源優(yōu)化配置的前提和基礎(chǔ),通過(guò)水資源需求量預(yù)測(cè)研究,提高區(qū)域水資源需求量預(yù)測(cè)精度,有利于協(xié)調(diào)區(qū)域水資源供需平衡。

    目前,許多學(xué)者對(duì)區(qū)域水資源需求量預(yù)測(cè)模型進(jìn)行了系統(tǒng)研究,王好芳等[1]建立了GM(1,1)模型預(yù)測(cè)區(qū)域社會(huì)經(jīng)濟(jì)系統(tǒng)的水資源需求量;方浩等[2]將神經(jīng)網(wǎng)絡(luò)和模糊理論相結(jié)合,建立了模糊神經(jīng)網(wǎng)絡(luò)模型,從模糊神經(jīng)網(wǎng)絡(luò)角度并運(yùn)用灰色系統(tǒng)理論預(yù)測(cè)區(qū)域水資源需求量;白雪華等[3]利用改進(jìn)的RBF神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法預(yù)測(cè)區(qū)域水資源需求量。郝轉(zhuǎn)等[4]利用GM(1,1)等維新息模型預(yù)測(cè)區(qū)域水資源需求量;羅利民等[5]將神經(jīng)網(wǎng)絡(luò)和小波分析理論相結(jié)合,建立了小波神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)區(qū)域水資源需求量;陳坤[6]在系統(tǒng)闡述各類水資源需求量預(yù)測(cè)方法優(yōu)缺點(diǎn)的基礎(chǔ)上,通過(guò)分析區(qū)域各行業(yè)需水量的影響因素,建立了回歸模型預(yù)測(cè)區(qū)域水資源需求量;左其亭等[7]提出了社會(huì)經(jīng)濟(jì)安全量化方法,將社會(huì)經(jīng)濟(jì)安全條件與需水量預(yù)測(cè)方法相結(jié)合,建立了社會(huì)經(jīng)濟(jì)安全條件下的區(qū)域水資源需求量預(yù)測(cè)模型。續(xù)繽等[8]將城區(qū)用水量的非平穩(wěn)時(shí)序分解為趨勢(shì)項(xiàng)和殘差項(xiàng)兩部分,采用Logistic成長(zhǎng)曲線函數(shù)擬合趨勢(shì)項(xiàng),采用ARIMA模型擬合殘差項(xiàng),建立了城市水資源需求量的Logistic-ARIMA耦合預(yù)測(cè)模型。

    鑒于國(guó)內(nèi)外相關(guān)研究成果,針對(duì)區(qū)域水資源需求量的預(yù)測(cè),用水定額法、人均綜合用水量法、指標(biāo)分析法、回歸分析法、時(shí)間序列法、灰色預(yù)測(cè)法以及人工神經(jīng)網(wǎng)絡(luò)法等理論方法都得到了廣泛應(yīng)用。但是,用水定額法、人均綜合用水量法以及指標(biāo)分析法的用水定額難以精確計(jì)量;回歸分析法僅考慮了預(yù)測(cè)對(duì)象與影響因素之間的因果關(guān)系;時(shí)間序列法僅考慮了預(yù)測(cè)對(duì)象隨時(shí)間變化的發(fā)展趨勢(shì);灰色預(yù)測(cè)法未考慮系統(tǒng)內(nèi)各變量之間的相互關(guān)系和相互作用;人工神經(jīng)網(wǎng)絡(luò)法等其他方法往往分析過(guò)程比較復(fù)雜。鑒于此,基于各類水資源需求量預(yù)測(cè)方法存在的不足之處,為提高預(yù)測(cè)的準(zhǔn)確性和精度,可結(jié)合2種或多種預(yù)測(cè)方法,建立組合預(yù)測(cè)模型,對(duì)區(qū)域水資源需求量進(jìn)行組合預(yù)測(cè)。本文根據(jù)區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展的用水現(xiàn)狀以及水資源需求的總體變化趨勢(shì),將Logistic生長(zhǎng)曲線模型和GM(1,1)等維新息模型相結(jié)合,建立組合預(yù)測(cè)模型,組合預(yù)測(cè)區(qū)域水資源需求量。在此基礎(chǔ)上,利用ARIMA模型對(duì)區(qū)域水資源需求量進(jìn)行一級(jí)修正;并根據(jù)區(qū)域節(jié)水規(guī)劃的總體目標(biāo),剖析區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展綜合節(jié)水潛力,對(duì)區(qū)域水資源需求量一級(jí)修正結(jié)果進(jìn)行二級(jí)修正,以確定區(qū)域水資源需求量的綜合預(yù)測(cè)結(jié)果。

    1 基于組合預(yù)測(cè)模型的區(qū)域水資源需求量預(yù)測(cè)

    組合預(yù)測(cè)是用2個(gè)或2個(gè)以上不同預(yù)測(cè)模型對(duì)同一預(yù)測(cè)對(duì)象進(jìn)行預(yù)測(cè),然后將各個(gè)單獨(dú)的預(yù)測(cè)結(jié)果進(jìn)行加權(quán),作為組合預(yù)測(cè)結(jié)果。組合預(yù)測(cè)可適度減少預(yù)測(cè)風(fēng)險(xiǎn),提高預(yù)測(cè)精度,其關(guān)鍵是恰當(dāng)確定單個(gè)預(yù)測(cè)模型的權(quán)重系數(shù)。基于國(guó)內(nèi)外相關(guān)研究成果,本文將Logistic生長(zhǎng)曲線模型和GM(1,1)等維新息模型相結(jié)合,建立組合模型,對(duì)區(qū)域水資源需求量進(jìn)行組合預(yù)測(cè)。

    區(qū)域水資源需求量預(yù)測(cè)的組合預(yù)測(cè)模型可表述為

    1.1 組合預(yù)測(cè)模型y1t預(yù)測(cè)值的確定

    Logistic生長(zhǎng)曲線是比利時(shí)數(shù)學(xué)家維哈爾斯特(P.F.Verhulst)發(fā)現(xiàn)的一種特殊曲線。皮爾(R.Pearl)和里德(L.J.Reed)曾經(jīng)應(yīng)用這一曲線研究人口增長(zhǎng)規(guī)律,所以Logistic生長(zhǎng)曲線又稱為Pearl-Reed曲線。由于區(qū)域水資源需求量的總體變化趨勢(shì)與Logistic生長(zhǎng)曲線模型的變化趨勢(shì)相吻合,因此,可采用Logistic生長(zhǎng)曲線模型預(yù)測(cè)區(qū)域水資源需求總量的發(fā)展趨勢(shì),提高區(qū)域水資源需求總量預(yù)測(cè)的準(zhǔn)確性。

    Logistic生長(zhǎng)曲線模型確定區(qū)域水資源需求量y1t預(yù)測(cè)值的表達(dá)式為

    式中:k、a、b為L(zhǎng)ogistic生長(zhǎng)曲線模型的參數(shù)。

    1.2 組合預(yù)測(cè)模型y2t預(yù)測(cè)值的確定

    灰色系統(tǒng)預(yù)測(cè)理論對(duì)于信息不完整或不完全的實(shí)際情況具有良好的適用性,其中GM(1,1)模型具有充分利用“少數(shù)據(jù)”進(jìn)行預(yù)測(cè)的優(yōu)點(diǎn),在區(qū)域水資源需求量預(yù)測(cè)中得到了較為廣泛的運(yùn)用。GM(1,1)等維新息模型在不斷補(bǔ)充新信息的同時(shí),及時(shí)去掉老化信息,更好地揭示了區(qū)域水資源需求量的發(fā)展趨勢(shì),由于GM(1,1)等維新息模型能夠獲得較高的區(qū)域水資源需求量預(yù)測(cè)精度。因此,采用GM(1,1)等維新息模型對(duì)區(qū)域水資源需求量進(jìn)行預(yù)測(cè)。

    GM(1,1)等維新息模型確定區(qū)域水資源需求量y2t預(yù)測(cè)值的具體步驟為:

    步驟1,記區(qū)域現(xiàn)狀用水量按照時(shí)間的先后順序排列而成的原始數(shù)列為x(0):

    根據(jù)灰色系統(tǒng)理論對(duì)區(qū)域現(xiàn)狀用水量的原始數(shù)列進(jìn)行一階累加(1-AGO)生成后,得生成列x(1),即

    1.3 組合預(yù)測(cè)模型權(quán)重w1與w2的確定

    式(1)中權(quán)重w1與w2可采用誤差平方和或誤差絕對(duì)值和最小準(zhǔn)則確定。本文采用誤差平方和最小準(zhǔn)則,即

    式中yt為第t年區(qū)域?qū)嶋H用水量,目標(biāo)函數(shù)Zt是一個(gè)二次函數(shù),并且是一個(gè)凸函數(shù),由Kuhn-Tucker條件可知,它的任意一個(gè)Kuhn-Tucker點(diǎn)就是整體最小點(diǎn)。建立廣義拉格朗日函數(shù):

    該廣義拉格朗日函數(shù)的Kuhn-Tucker點(diǎn)應(yīng)滿足ΔZt(w,λ)=0,因此,廣義拉格朗日函數(shù)的解為

    將權(quán)重w1、w2代入式(1),可確定第t年區(qū)域水資源需求量的組合預(yù)測(cè)值y^t。

    2 基于二級(jí)修正的區(qū)域水資源需求量預(yù)測(cè)

    2.1 區(qū)域水資源需求量的一級(jí)修正結(jié)果

    基于組合預(yù)測(cè)模型,組合預(yù)測(cè)得到區(qū)域水資源需求量。為進(jìn)一步減少預(yù)測(cè)誤差,提高預(yù)測(cè)精度,利用ARIMA模型[9]對(duì)組合預(yù)測(cè)值t進(jìn)行一級(jí)修正,修正與實(shí)際水資源需求量的差值Yt=yt-t,從而確定組合預(yù)測(cè)模型和ARIMA模型的疊加預(yù)測(cè)值。

    ARIMA模型由美國(guó)學(xué)者George Box和英國(guó)統(tǒng)計(jì)學(xué)家Gwilym Jenkinsy于20世紀(jì)70年代首次提出。它是利用自回歸模型(auto regressive模型,簡(jiǎn)稱AR模型)與移動(dòng)平均模型(moving average模型,簡(jiǎn)稱MA模型)搭配形成的隨機(jī)型時(shí)序ARIMA模型;ARIMA模型中有多個(gè)參數(shù),無(wú)季節(jié)成分時(shí)可記為ARIMA(p,d,q)模型,其中,p為自回歸過(guò)程階數(shù),d為差分的階數(shù),q為移動(dòng)平均過(guò)程階數(shù)。ARIMA(p,d,q)模型經(jīng)d階差分之后成為ARMA(p,q)模型,采用極大似然法、矩估計(jì)法和最小二乘法等參數(shù)估計(jì)方法,對(duì)ARMA(p,q)模型進(jìn)行參數(shù)估計(jì),將估計(jì)參數(shù)和最佳模型階數(shù)代入隨機(jī)差分方程,可確定ARMA(p,q)模型。

    設(shè)(Y1,Y2,…,Yt,…,Yn)(t=1,2,…,n)表示組合預(yù)測(cè)值與實(shí)際水資源需求量結(jié)果的差值Yt經(jīng)d階差分后組成的平穩(wěn)時(shí)間序列,(Y1,Y2,…,Yt,…,Yn)可擬合成具有p階自回歸部分、q階移動(dòng)平均部分的ARMA(p,q)模型隨機(jī)差分方程:

    式中:φi(i=1,2,…,p)為自回歸系數(shù);θi(i=0,1,…,q)為移動(dòng)平均系數(shù);εt為零均值、方差為 δ2pq的白噪聲序列。ARMA(p,q)模型參數(shù)估計(jì)的具體計(jì)算步驟可參考文獻(xiàn)[8-9]。

    2.2 區(qū)域水資源需求量的二級(jí)修正結(jié)果

    基于區(qū)域水資源需求量的一級(jí)修正結(jié)果,根據(jù)區(qū)域節(jié)水規(guī)劃的總體目標(biāo),采用區(qū)域各用水行業(yè)節(jié)水潛力的計(jì)量方法[10],剖析區(qū)域生活、生態(tài)環(huán)境以及生產(chǎn)等“三生”用水行業(yè)發(fā)展的節(jié)水潛力以及區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展綜合節(jié)水潛力,通過(guò)專家咨詢,從定性和定量相結(jié)合的角度分析規(guī)劃年區(qū)域綜合節(jié)水潛力的取值范圍,對(duì)區(qū)域水資源需求量一級(jí)修正結(jié)果進(jìn)行二級(jí)修正,得到規(guī)劃年區(qū)域水資源需求量的取值范圍。

    3 算例分析

    以2000—2005年H市城區(qū)用水量數(shù)據(jù)(表1)為依據(jù),基于二級(jí)修正的水資源需求量預(yù)測(cè)模型,預(yù)測(cè)規(guī)劃年H市城區(qū)水資源需求量。

    表1 2000~2005年H市城區(qū)用水量 萬(wàn)m3

    根據(jù)表1中的數(shù)據(jù),利用Logistic生長(zhǎng)曲線模型與GM(1,1)等維新息模型,組合預(yù)測(cè)2010年、2015年、2020年H市城區(qū)水資源需求量,見表2。

    表2 規(guī)劃年H市城區(qū)需水量 萬(wàn)m3

    針對(duì)組合預(yù)測(cè)模型得到的H市城區(qū)需水量,采用ARIMA(1,0,1)模型對(duì)其差值進(jìn)行擬合分析,差值擬合后的變化趨勢(shì)見圖1。

    圖1 ARIMA(1,0,1)模型需水量差值擬合

    圖1表明,利用ARIMA(1,0,1)模型對(duì)H市城區(qū)需水量的差值進(jìn)行擬合時(shí),2000—2010年差值的變化表現(xiàn)出大幅度的波動(dòng)現(xiàn)象,而2010—2020年差值的變化趨于平穩(wěn)。這說(shuō)明H市城區(qū)的需水量預(yù)測(cè)值基本符合Logistic生長(zhǎng)曲線的變化趨勢(shì),H市城區(qū)需水量的變化仍然處于Logistic生長(zhǎng)曲線的上升階段,即在2010—2020年H市城區(qū)需水量將繼續(xù)保持持續(xù)增長(zhǎng)的趨勢(shì)。將組合預(yù)測(cè)模型和ARIMA預(yù)測(cè)模型進(jìn)行耦合,可進(jìn)一步提高H市城區(qū)需水量預(yù)測(cè)的精度,得到H市城區(qū)各規(guī)劃年需水量的一級(jí)修正結(jié)果;根據(jù)規(guī)劃年H市城區(qū)水資源需求量一級(jí)修正結(jié)果,結(jié)合H市城區(qū)節(jié)水規(guī)劃的總體目標(biāo),分析H市城區(qū)規(guī)劃年綜合節(jié)水潛力,得到H市城區(qū)各規(guī)劃年需水量的二級(jí)修正結(jié)果,見表3。

    表3 規(guī)劃年H市城區(qū)需水量的修正結(jié)果 萬(wàn)m3

    表3結(jié)果表明,H市城區(qū)水資源需求量預(yù)測(cè)值與其水資源需求量的發(fā)展趨勢(shì)相吻合。通過(guò)與文獻(xiàn)[8]建立的Logistic-ARIMA耦合模型的預(yù)測(cè)結(jié)果進(jìn)行對(duì)比可知,基于二級(jí)修正的區(qū)域水資源需求量預(yù)測(cè)模型的預(yù)測(cè)結(jié)果相對(duì)偏低,根本原因在于,模型是在采用Logistic生長(zhǎng)曲線模型和GM(1,1)等維新息模型對(duì)H市城區(qū)水資源需求量進(jìn)行組合預(yù)測(cè)的基礎(chǔ)上,采用ARIMA模型,并結(jié)合H市城區(qū)節(jié)水規(guī)劃總體目標(biāo),分析綜合節(jié)水潛力,對(duì)H市城區(qū)水資源需求量預(yù)測(cè)值進(jìn)行二級(jí)修正,因此,其需水量預(yù)測(cè)結(jié)果更加符合發(fā)展趨勢(shì)。

    4 結(jié) 語(yǔ)

    結(jié)合Logistic生長(zhǎng)曲線模型和GM(1,1)等維新息模型,建立了區(qū)域水資源需求量組合預(yù)測(cè)模型,在此基礎(chǔ)上,采用ARIMA模型,根據(jù)區(qū)域節(jié)水規(guī)劃的總體目標(biāo),對(duì)區(qū)域水資源需求量組合預(yù)測(cè)值進(jìn)行修正,確定基于二級(jí)修正的區(qū)域水資源需求量綜合預(yù)測(cè)值。結(jié)合H市城區(qū)社會(huì)經(jīng)濟(jì)用水量數(shù)據(jù)進(jìn)行算例分析,研究結(jié)果表明,預(yù)測(cè)模型具有較高的預(yù)測(cè)精度。

    [1]王好芳,董增川.區(qū)域經(jīng)濟(jì)-人口系統(tǒng)需水量的灰色預(yù)測(cè)方法研究[J].西北水資源與水工程,2001,12(2):53-54.

    [2]方浩,李蓓,石娜,等.基于模糊神經(jīng)網(wǎng)絡(luò)的區(qū)域需水預(yù)測(cè)計(jì)算模型[J].西北水資源與水工程,2003,14(4):1-3.

    [3]白雪華,郭旭穎.改進(jìn)的RBF網(wǎng)絡(luò)在區(qū)域需水預(yù)測(cè)中的應(yīng)用[J].青島建筑工程學(xué)院學(xué)報(bào),2005,26(3):87-89.

    [4]郝轉(zhuǎn),張文鴿.GM(1,1)等維新息模型在區(qū)域需水量預(yù)測(cè)中的應(yīng)用[J].東北水利水電,2006,24(4):6-8.

    [5]羅利民,方浩,仲躍,等.小波神經(jīng)網(wǎng)絡(luò)算法在區(qū)域需水預(yù)測(cè)中的應(yīng)用[J].計(jì)算機(jī)工程與應(yīng)用,2006(3):200-201.

    [6]陳坤.上海水資源可持續(xù)利用的經(jīng)濟(jì)學(xué)研究[M].上海:上海人民出版社,2007:49-55.

    [7]左其亭,李琳.社會(huì)經(jīng)濟(jì)安全條件下區(qū)域需水量預(yù)測(cè)方法[J].水資源保護(hù),2008,24(1):6-11.

    [8]續(xù)繽,郭立山.城市水資源需求量預(yù)測(cè)Logistic-ARIMA耦合模型[J].水利經(jīng)濟(jì),2008,26(5):14-16.

    [9]寧宣熙,劉思峰.管理預(yù)測(cè)與決策方法[M].北京:科學(xué)出版社,1997:98-100.

    [10]宋長(zhǎng)虹,高淑紅,胡志亮.節(jié)水潛力的計(jì)算方法與應(yīng)用[J].黑龍江水利科技,2008,36(1):18.

    Regional water demand forecast based on second-order modifying model

    WU Dan,WU Feng-ping,CHEN Yan-ping
    (Business School,Hohai University,Nanjing 210098,China)

    The Logistic model and GM(1,1)constant dimension new information model were combined to establish a comprehensive forecasting model to forecast regional water demand based on current water use and the total development trend of water demand in the course of regional socio-economic sustainable development.The forecast of regional water demand was modified using the ARIMA model for first-order modification.According to the overall objective of regional water resources planning and analysis of the comprehensive regional socio-economic water-saving potential,the first-order modification results of regional water demand were refined by second-order modification,and the final results of regional water demand were obtained.The results showed that the model had high predictive accuracy and actual effect.

    water demand;Logistic model;GM(1,1)constant dimension new information model;ARIMA model;second-order modification

    TV213

    A

    1004-6933(2010)03-0001-04

    國(guó)家“十一五”科技支撐計(jì)劃重點(diǎn)項(xiàng)目(2007BAB28B03);江蘇省國(guó)家重點(diǎn)學(xué)科培育建設(shè)點(diǎn)技術(shù)經(jīng)濟(jì)及管理項(xiàng)目。

    吳丹(1986—),男,江西撫州人,博士研究生,主要從事資源配置與規(guī)劃管理研究。E-mail:wu_daniel@hhu.edu.cn

    (收稿日期:2008-12-04 編輯:熊水斌)

    猜你喜歡
    需水量需求量城區(qū)
    長(zhǎng)沙市望城區(qū)金地三千府幼兒園
    幼兒畫刊(2022年8期)2022-10-18 01:44:10
    從數(shù)學(xué)角度看“彈性”
    南京城市中心區(qū)生態(tài)需水量探析
    金霞早油蟠在保定滿城區(qū)的表現(xiàn)及栽培技術(shù)
    河北果樹(2020年2期)2020-05-25 06:58:26
    組合預(yù)測(cè)模型在區(qū)域需水量預(yù)測(cè)中的應(yīng)用
    福田要建健康城區(qū)
    帕滿灌區(qū)供需水量平衡分析
    大洋河生態(tài)需水量分析與評(píng)價(jià)
    對(duì)淮安市城區(qū)河道生態(tài)建設(shè)的思考
    2017年我國(guó)汽車軟管需求量將達(dá)6.4億m
    橡膠科技(2015年3期)2015-02-26 14:45:02
    一本—道久久a久久精品蜜桃钙片| 亚洲av在线观看美女高潮| 1024香蕉在线观看| 亚洲av日韩在线播放| 久久精品久久久久久噜噜老黄| 日本午夜av视频| 成人国产av品久久久| 不卡av一区二区三区| 婷婷成人精品国产| 91精品三级在线观看| 纵有疾风起免费观看全集完整版| 亚洲国产中文字幕在线视频| 老司机在亚洲福利影院| 午夜免费观看性视频| 国产成人精品久久久久久| 免费不卡黄色视频| 亚洲国产最新在线播放| 国产野战对白在线观看| 国产精品av久久久久免费| 欧美日韩一级在线毛片| 人人妻人人添人人爽欧美一区卜| 婷婷色av中文字幕| 亚洲,欧美精品.| 精品福利永久在线观看| 男女边吃奶边做爰视频| av女优亚洲男人天堂| 午夜影院在线不卡| 最近最新中文字幕大全免费视频 | 午夜91福利影院| 国产欧美日韩一区二区三区在线| 久久久国产精品麻豆| 国产熟女午夜一区二区三区| 久久97久久精品| 国产无遮挡羞羞视频在线观看| 久久天堂一区二区三区四区| 青草久久国产| 日日爽夜夜爽网站| 免费黄网站久久成人精品| 大香蕉久久网| 亚洲国产欧美一区二区综合| 亚洲欧美一区二区三区黑人| 女人久久www免费人成看片| 精品一品国产午夜福利视频| 国产精品久久久久久人妻精品电影 | 高清不卡的av网站| 国产深夜福利视频在线观看| 丝袜美腿诱惑在线| 啦啦啦在线观看免费高清www| 亚洲精品国产av成人精品| 精品人妻一区二区三区麻豆| 亚洲国产精品一区二区三区在线| 国产亚洲午夜精品一区二区久久| 日韩中文字幕欧美一区二区 | 久久99精品国语久久久| 免费观看性生交大片5| 久久久久久免费高清国产稀缺| 18禁国产床啪视频网站| 国产熟女午夜一区二区三区| 天天躁日日躁夜夜躁夜夜| 亚洲国产精品国产精品| 人体艺术视频欧美日本| 午夜激情久久久久久久| 在线免费观看不下载黄p国产| 国产精品久久久人人做人人爽| 日韩人妻精品一区2区三区| 午夜久久久在线观看| av女优亚洲男人天堂| 午夜福利乱码中文字幕| 久久免费观看电影| 午夜免费鲁丝| 午夜精品国产一区二区电影| 叶爱在线成人免费视频播放| 国产精品亚洲av一区麻豆 | 十八禁人妻一区二区| av一本久久久久| 久久久久国产精品人妻一区二区| 中文字幕人妻丝袜一区二区 | av在线老鸭窝| 99国产综合亚洲精品| e午夜精品久久久久久久| 精品国产国语对白av| 女人高潮潮喷娇喘18禁视频| 欧美精品一区二区大全| 青青草视频在线视频观看| 免费看av在线观看网站| 久久鲁丝午夜福利片| 国产欧美日韩综合在线一区二区| 亚洲婷婷狠狠爱综合网| 秋霞伦理黄片| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利影视在线免费观看| 国产xxxxx性猛交| 国产成人精品久久二区二区91 | 午夜福利视频在线观看免费| 永久免费av网站大全| 免费久久久久久久精品成人欧美视频| 亚洲欧美精品综合一区二区三区| 久久天堂一区二区三区四区| 日本色播在线视频| 亚洲欧美一区二区三区久久| 日韩 欧美 亚洲 中文字幕| 美女大奶头黄色视频| 亚洲精品国产色婷婷电影| 国产乱来视频区| 在线精品无人区一区二区三| 欧美日韩福利视频一区二区| 免费在线观看视频国产中文字幕亚洲 | 我要看黄色一级片免费的| 天天影视国产精品| 久久99精品国语久久久| av国产久精品久网站免费入址| 亚洲欧美色中文字幕在线| 纵有疾风起免费观看全集完整版| 看免费av毛片| 天天躁狠狠躁夜夜躁狠狠躁| 国产女主播在线喷水免费视频网站| 在线观看三级黄色| 99久久99久久久精品蜜桃| 亚洲精品aⅴ在线观看| 国产精品久久久久久人妻精品电影 | 99精国产麻豆久久婷婷| 黄片无遮挡物在线观看| 熟女少妇亚洲综合色aaa.| 欧美日韩国产mv在线观看视频| 性色av一级| 十分钟在线观看高清视频www| 在线观看国产h片| 国产精品久久久久久精品电影小说| 18禁观看日本| 久久av网站| 麻豆av在线久日| 日韩精品有码人妻一区| 欧美中文综合在线视频| 亚洲国产精品一区三区| 男男h啪啪无遮挡| 黄网站色视频无遮挡免费观看| 国产一区有黄有色的免费视频| 亚洲久久久国产精品| 大香蕉久久成人网| av天堂久久9| 黄色毛片三级朝国网站| 69精品国产乱码久久久| 这个男人来自地球电影免费观看 | av国产久精品久网站免费入址| 啦啦啦啦在线视频资源| 亚洲四区av| 搡老乐熟女国产| 午夜免费鲁丝| 老司机影院成人| 久久久久久久国产电影| 国产色婷婷99| 十八禁网站网址无遮挡| 亚洲国产欧美一区二区综合| 狂野欧美激情性xxxx| 欧美日韩视频高清一区二区三区二| 日本欧美国产在线视频| 午夜精品国产一区二区电影| 精品一品国产午夜福利视频| 国产av精品麻豆| 国产极品粉嫩免费观看在线| 久久久精品94久久精品| 亚洲av成人不卡在线观看播放网 | 99精国产麻豆久久婷婷| 十八禁网站网址无遮挡| 欧美在线一区亚洲| 99热国产这里只有精品6| √禁漫天堂资源中文www| 欧美 亚洲 国产 日韩一| 中文字幕制服av| h视频一区二区三区| 一级毛片黄色毛片免费观看视频| 少妇被粗大猛烈的视频| 久久热在线av| 国产精品久久久久久精品电影小说| 十八禁高潮呻吟视频| 久久精品国产亚洲av涩爱| 麻豆av在线久日| 啦啦啦在线观看免费高清www| 亚洲自偷自拍图片 自拍| 国产深夜福利视频在线观看| 久久av网站| 国产成人91sexporn| 欧美日韩精品网址| 国产精品av久久久久免费| 欧美日韩亚洲国产一区二区在线观看 | 五月天丁香电影| 乱人伦中国视频| 十八禁网站网址无遮挡| 你懂的网址亚洲精品在线观看| av片东京热男人的天堂| 久久久久久免费高清国产稀缺| 大陆偷拍与自拍| 久久久精品94久久精品| 熟妇人妻不卡中文字幕| av在线观看视频网站免费| netflix在线观看网站| 午夜免费男女啪啪视频观看| 最新在线观看一区二区三区 | 日韩制服丝袜自拍偷拍| 电影成人av| 午夜福利视频精品| avwww免费| 日本色播在线视频| 侵犯人妻中文字幕一二三四区| 久久女婷五月综合色啪小说| 搡老乐熟女国产| 亚洲精品成人av观看孕妇| 黑人猛操日本美女一级片| 又粗又硬又长又爽又黄的视频| 99久久精品国产亚洲精品| 各种免费的搞黄视频| 2021少妇久久久久久久久久久| 亚洲伊人色综图| 女性生殖器流出的白浆| 日本一区二区免费在线视频| av在线播放精品| 午夜精品国产一区二区电影| videosex国产| 视频区图区小说| 大片免费播放器 马上看| 国产野战对白在线观看| 国产免费现黄频在线看| 18在线观看网站| 日本wwww免费看| 国产激情久久老熟女| 日韩免费高清中文字幕av| 嫩草影视91久久| 成人亚洲精品一区在线观看| 国产精品一国产av| 男女高潮啪啪啪动态图| 亚洲国产欧美网| 女的被弄到高潮叫床怎么办| 国产一卡二卡三卡精品 | 日本vs欧美在线观看视频| 美女福利国产在线| 国产精品久久久久久久久免| 久久久亚洲精品成人影院| 精品国产国语对白av| 国产伦理片在线播放av一区| 欧美另类一区| 亚洲综合色网址| 亚洲婷婷狠狠爱综合网| 一本一本久久a久久精品综合妖精| 中国三级夫妇交换| 男人舔女人的私密视频| 黄片无遮挡物在线观看| 在线 av 中文字幕| 亚洲精华国产精华液的使用体验| 国产成人91sexporn| 水蜜桃什么品种好| 国产在线视频一区二区| 日韩av免费高清视频| 熟妇人妻不卡中文字幕| 欧美人与性动交α欧美精品济南到| 97人妻天天添夜夜摸| 亚洲成色77777| 免费av中文字幕在线| 久久久久久久国产电影| 大陆偷拍与自拍| 亚洲av福利一区| 我要看黄色一级片免费的| 男人添女人高潮全过程视频| 国产黄色免费在线视频| 日韩av免费高清视频| 另类精品久久| 国产免费现黄频在线看| 一本色道久久久久久精品综合| 午夜91福利影院| 国产亚洲精品第一综合不卡| 最近中文字幕2019免费版| 亚洲国产av新网站| 色播在线永久视频| 精品亚洲乱码少妇综合久久| 成人免费观看视频高清| 国产麻豆69| 一区二区三区乱码不卡18| av在线老鸭窝| 免费女性裸体啪啪无遮挡网站| 自拍欧美九色日韩亚洲蝌蚪91| 国产片特级美女逼逼视频| 久久毛片免费看一区二区三区| 男女边摸边吃奶| 国产一区二区 视频在线| 亚洲美女搞黄在线观看| 国产高清不卡午夜福利| 亚洲成国产人片在线观看| 精品少妇久久久久久888优播| 韩国高清视频一区二区三区| 国产午夜精品一二区理论片| 最新的欧美精品一区二区| 国产片内射在线| 丰满饥渴人妻一区二区三| 十八禁网站网址无遮挡| 国产欧美日韩综合在线一区二区| 国产成人a∨麻豆精品| 嫩草影视91久久| 欧美日韩一区二区视频在线观看视频在线| 大码成人一级视频| 国产精品免费视频内射| 亚洲精华国产精华液的使用体验| 国产精品一国产av| 精品少妇内射三级| 韩国精品一区二区三区| 日韩精品有码人妻一区| 欧美日韩亚洲高清精品| 哪个播放器可以免费观看大片| 亚洲欧美精品综合一区二区三区| 久久99一区二区三区| 国产免费一区二区三区四区乱码| 黄色一级大片看看| 久久国产亚洲av麻豆专区| 人妻一区二区av| 亚洲欧洲日产国产| 成人午夜精彩视频在线观看| 又粗又硬又长又爽又黄的视频| 亚洲精品一区蜜桃| 99久久精品国产亚洲精品| 999久久久国产精品视频| 国产精品无大码| 韩国av在线不卡| 纵有疾风起免费观看全集完整版| 曰老女人黄片| 男的添女的下面高潮视频| 美国免费a级毛片| 国产在线免费精品| 亚洲,一卡二卡三卡| 亚洲国产精品999| 伦理电影免费视频| 熟女av电影| 亚洲精品国产色婷婷电影| 欧美在线黄色| 日本猛色少妇xxxxx猛交久久| 日本欧美视频一区| 欧美日韩亚洲国产一区二区在线观看 | 国产片特级美女逼逼视频| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片我不卡| 老司机影院毛片| 男男h啪啪无遮挡| 岛国毛片在线播放| 亚洲精品av麻豆狂野| 捣出白浆h1v1| 日本av手机在线免费观看| 一本一本久久a久久精品综合妖精| 欧美精品av麻豆av| 午夜免费男女啪啪视频观看| 搡老乐熟女国产| 亚洲视频免费观看视频| av天堂久久9| 秋霞伦理黄片| 99久久精品国产亚洲精品| 99九九在线精品视频| av一本久久久久| 久久久精品区二区三区| 午夜福利免费观看在线| 黄片播放在线免费| 亚洲国产毛片av蜜桃av| 亚洲国产精品成人久久小说| 日本av手机在线免费观看| 一边摸一边做爽爽视频免费| 亚洲成人国产一区在线观看 | 高清欧美精品videossex| 亚洲欧美成人精品一区二区| 激情五月婷婷亚洲| 超碰97精品在线观看| 国产精品久久久久久人妻精品电影 | 不卡av一区二区三区| 亚洲,欧美,日韩| 一级毛片 在线播放| 一区二区三区激情视频| 亚洲一区二区三区欧美精品| 国产精品99久久99久久久不卡 | 少妇 在线观看| 欧美精品一区二区大全| 中文字幕精品免费在线观看视频| 免费女性裸体啪啪无遮挡网站| 免费观看性生交大片5| 深夜精品福利| 亚洲情色 制服丝袜| 精品亚洲成a人片在线观看| 亚洲美女黄色视频免费看| 日韩一区二区三区影片| 免费人妻精品一区二区三区视频| 丝袜脚勾引网站| 国产av码专区亚洲av| 又黄又粗又硬又大视频| 国产男女超爽视频在线观看| 国产人伦9x9x在线观看| 黄片小视频在线播放| 亚洲精品国产色婷婷电影| av视频免费观看在线观看| 国产男人的电影天堂91| 精品一区二区三卡| 欧美中文综合在线视频| 午夜免费观看性视频| 免费av中文字幕在线| 精品视频人人做人人爽| 男女高潮啪啪啪动态图| 国产又爽黄色视频| 久久99精品国语久久久| 国产成人av激情在线播放| 伊人久久大香线蕉亚洲五| 亚洲人成77777在线视频| 2018国产大陆天天弄谢| 国产淫语在线视频| 国产精品 欧美亚洲| 欧美黑人精品巨大| 久久久久精品国产欧美久久久 | 麻豆精品久久久久久蜜桃| 伊人亚洲综合成人网| 国产 精品1| 国产精品 国内视频| 最近的中文字幕免费完整| 亚洲精品自拍成人| 另类精品久久| 免费在线观看黄色视频的| 男人舔女人的私密视频| av片东京热男人的天堂| 亚洲欧美色中文字幕在线| 777米奇影视久久| 精品国产一区二区三区四区第35| 中文精品一卡2卡3卡4更新| 亚洲国产看品久久| 成人三级做爰电影| 另类亚洲欧美激情| 成人影院久久| 久久久久精品人妻al黑| av国产久精品久网站免费入址| 精品国产一区二区三区久久久樱花| 精品久久久精品久久久| 纵有疾风起免费观看全集完整版| 久久久久国产精品人妻一区二区| 在线观看免费日韩欧美大片| 久久久久久久大尺度免费视频| bbb黄色大片| 午夜精品国产一区二区电影| 色婷婷av一区二区三区视频| 精品亚洲成a人片在线观看| 国产男人的电影天堂91| 宅男免费午夜| 国产av国产精品国产| videosex国产| 国产精品久久久久成人av| 精品免费久久久久久久清纯 | 久久久久久久精品精品| 久久鲁丝午夜福利片| 在线观看免费高清a一片| 日本午夜av视频| 性色av一级| 69精品国产乱码久久久| 亚洲av成人不卡在线观看播放网 | 日本一区二区免费在线视频| 亚洲少妇的诱惑av| 成人国产av品久久久| 日韩不卡一区二区三区视频在线| 涩涩av久久男人的天堂| 欧美国产精品一级二级三级| 亚洲,欧美,日韩| 夫妻午夜视频| 一本一本久久a久久精品综合妖精| xxxhd国产人妻xxx| 99精国产麻豆久久婷婷| 在线 av 中文字幕| 国产极品天堂在线| 丰满饥渴人妻一区二区三| 丝袜美腿诱惑在线| 久久午夜综合久久蜜桃| 国产成人系列免费观看| 免费人妻精品一区二区三区视频| 老司机亚洲免费影院| 亚洲成色77777| 国产亚洲最大av| 18禁国产床啪视频网站| 热re99久久国产66热| 国产欧美亚洲国产| 久久久久久久精品精品| 久久精品人人爽人人爽视色| 一边摸一边抽搐一进一出视频| 免费观看av网站的网址| 精品久久久精品久久久| 满18在线观看网站| 亚洲av综合色区一区| 久久久久精品久久久久真实原创| 久久久国产欧美日韩av| 爱豆传媒免费全集在线观看| 亚洲精品久久成人aⅴ小说| 搡老乐熟女国产| 日韩,欧美,国产一区二区三区| 亚洲国产欧美网| 在线看a的网站| 精品少妇一区二区三区视频日本电影 | 激情五月婷婷亚洲| 国产精品.久久久| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品久久二区二区91 | 精品福利永久在线观看| 另类精品久久| 国产精品亚洲av一区麻豆 | 久久97久久精品| 久久人妻熟女aⅴ| 亚洲欧美日韩另类电影网站| 在线观看免费午夜福利视频| 精品国产一区二区三区四区第35| 久久久久久久精品精品| 精品免费久久久久久久清纯 | 久久久久精品人妻al黑| 亚洲自偷自拍图片 自拍| 1024视频免费在线观看| 午夜精品国产一区二区电影| 2018国产大陆天天弄谢| 国产欧美亚洲国产| a级毛片黄视频| 大香蕉久久成人网| 亚洲伊人色综图| 一级片'在线观看视频| 人妻一区二区av| 国产免费一区二区三区四区乱码| 美女大奶头黄色视频| 亚洲精品中文字幕在线视频| 色吧在线观看| 深夜精品福利| av不卡在线播放| 免费高清在线观看视频在线观看| 国产深夜福利视频在线观看| www.精华液| 五月开心婷婷网| 亚洲精品日本国产第一区| 久久性视频一级片| 涩涩av久久男人的天堂| 久久热在线av| 欧美日韩综合久久久久久| 国产精品免费视频内射| 视频区图区小说| 日韩电影二区| av国产精品久久久久影院| 免费女性裸体啪啪无遮挡网站| 性色av一级| 国产成人欧美在线观看 | 亚洲色图 男人天堂 中文字幕| 波多野结衣av一区二区av| 国产精品一区二区在线不卡| 国产精品一区二区在线观看99| 黄色视频在线播放观看不卡| 别揉我奶头~嗯~啊~动态视频 | 久久久国产欧美日韩av| 久久综合国产亚洲精品| 国产熟女欧美一区二区| 久久久久久免费高清国产稀缺| 国产成人欧美在线观看 | 99热国产这里只有精品6| 欧美97在线视频| 人人妻人人添人人爽欧美一区卜| av.在线天堂| 精品国产超薄肉色丝袜足j| 日日撸夜夜添| 亚洲,欧美,日韩| 老司机亚洲免费影院| 中文字幕人妻熟女乱码| 精品视频人人做人人爽| 91精品伊人久久大香线蕉| 亚洲激情五月婷婷啪啪| 老熟女久久久| svipshipincom国产片| 国产精品一区二区在线观看99| 国产亚洲av高清不卡| 精品国产一区二区久久| 另类亚洲欧美激情| 日韩免费高清中文字幕av| 香蕉国产在线看| 国产一区有黄有色的免费视频| 看免费av毛片| 纵有疾风起免费观看全集完整版| 九九爱精品视频在线观看| 最近最新中文字幕免费大全7| 新久久久久国产一级毛片| 久久久精品国产亚洲av高清涩受| 涩涩av久久男人的天堂| 国产成人av激情在线播放| 国产成人精品久久二区二区91 | 两性夫妻黄色片| 日韩 欧美 亚洲 中文字幕| av国产精品久久久久影院| 国产精品 国内视频| 女人高潮潮喷娇喘18禁视频| 亚洲av电影在线观看一区二区三区| 亚洲第一区二区三区不卡| 亚洲精品中文字幕在线视频| 日本猛色少妇xxxxx猛交久久| 国产精品.久久久| 欧美成人午夜精品| 一本色道久久久久久精品综合| 欧美人与善性xxx| 五月天丁香电影| 在线观看免费午夜福利视频| 久久天堂一区二区三区四区| 欧美 亚洲 国产 日韩一| 亚洲天堂av无毛| avwww免费| 欧美日韩国产mv在线观看视频| 亚洲美女黄色视频免费看| 成年动漫av网址| 日韩欧美精品免费久久| 黄片无遮挡物在线观看| 国产xxxxx性猛交| 国产免费视频播放在线视频| 久久99精品国语久久久| 美女中出高潮动态图| 中文字幕高清在线视频| 汤姆久久久久久久影院中文字幕| 国产又爽黄色视频| 色播在线永久视频| 国产黄频视频在线观看| 国产免费又黄又爽又色| 国产成人免费无遮挡视频| 婷婷色麻豆天堂久久|