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

基于主成分分析法的城市人水和諧度評(píng)價(jià)

2018-03-16 01:56:19孟令爽唐德善史毅超
關(guān)鍵詞:評(píng)價(jià)模型

孟令爽, 唐德善, 史毅超

(河海大學(xué) 水利水電學(xué)院, 江蘇 南京 210098)

1 研究背景

隨著社會(huì)發(fā)展與人民生活水平不斷提高,人類各種活動(dòng)逐漸向城市轉(zhuǎn)移,促使城市快速發(fā)展與城市人口的急劇增長(zhǎng),水-經(jīng)濟(jì)-社會(huì)-生態(tài)-人之間的矛盾逐漸凸顯。近年來(lái)隨著生態(tài)文明理念的不斷深入,人水問(wèn)題逐漸被提升到城市發(fā)展的戰(zhàn)略地位,人水和諧成為新時(shí)期治水思路的本質(zhì)要求[1]。因此,綜合評(píng)價(jià)城市人水和諧度,剖析人水和諧發(fā)展中存在的問(wèn)題,對(duì)于加速城市生態(tài)文明建設(shè)具有極其重要的意義。

近年來(lái)國(guó)內(nèi)外學(xué)者相繼從多個(gè)角度開(kāi)展了一系列的人水關(guān)系研究。我國(guó)學(xué)者汪恕誠(chéng)[2]提出要破解中國(guó)水問(wèn)題,就要促進(jìn)人與自然的和諧相處;左其亭等[3]論述了人水和諧的內(nèi)涵并提出城市人水和諧評(píng)價(jià)的指標(biāo)體系;戴會(huì)超等[4]從人水主觀感受度、人水和諧發(fā)展度、人水和諧協(xié)調(diào)度3個(gè)方面對(duì)南京的人水和諧程度進(jìn)行了研究,表明南京人水和諧度逐年提高;張金鑫等[5]基于云模型理論,建立了新型的人水和諧評(píng)價(jià)指標(biāo)體系,提出了基于云模型的人水和諧評(píng)價(jià)方法;英國(guó)學(xué)者Tony Allan創(chuàng)造性提出“虛擬水”的概念,將水滲透到了人類發(fā)展的各個(gè)方面,從而更加全面地把握城市人水和諧。

本文在現(xiàn)有研究的基礎(chǔ)上,從水資源、社會(huì)、經(jīng)濟(jì)和生態(tài)4個(gè)角度構(gòu)建城市人水和諧評(píng)價(jià)指標(biāo)體系,建立基于主成分分析法[6]的城市人水和諧等級(jí)評(píng)價(jià)模型,旨在對(duì)城市人水和諧度做出客觀合理的評(píng)價(jià)。

2 城市人水和諧內(nèi)涵及評(píng)價(jià)指標(biāo)體系

2.1 城市人水和諧內(nèi)涵

人文社會(huì)系統(tǒng)與水系統(tǒng)相互關(guān)聯(lián),相互影響共同組成復(fù)雜的人-水系統(tǒng)。城市人水系統(tǒng)可以簡(jiǎn)單地理解為城市人文系統(tǒng)與城市水系統(tǒng)之間的關(guān)系問(wèn)題,它包含3方面的內(nèi)容:(1)人文系統(tǒng)可持續(xù)發(fā)展;(2)水系統(tǒng)不斷改善;(3)水系統(tǒng)是人文系統(tǒng)發(fā)展的基礎(chǔ),人文系統(tǒng)發(fā)展是水系統(tǒng)健康的保障。

2.2 城市人水和諧度評(píng)價(jià)指標(biāo)體系

評(píng)價(jià)指標(biāo)體系的建立遵循以下原則[7]:系統(tǒng)性、科學(xué)合理性、客觀性、無(wú)關(guān)性、可操作性等。在綜合考慮各種因素、原則的前提下,客觀合理地選擇獨(dú)立并且具有代表性的指標(biāo),進(jìn)而對(duì)評(píng)價(jià)對(duì)象進(jìn)行系統(tǒng)客觀科學(xué)地評(píng)價(jià)。

本文在現(xiàn)有研究的基礎(chǔ)上,基于人水和諧的基本涵義,將指標(biāo)體系分為3個(gè)層面[8]:目標(biāo)層即城市人水和諧;準(zhǔn)則層為反映當(dāng)代城市人水關(guān)系的4大核心要素,即水資源情況、生態(tài)發(fā)展情況、經(jīng)濟(jì)發(fā)展情況、社會(huì)發(fā)展情況;指標(biāo)層為反映這4個(gè)核心方面的具體指標(biāo),通過(guò)DELPH法[9]、專家咨詢法[10]等對(duì)指標(biāo)進(jìn)行篩選整合,構(gòu)建了包含19項(xiàng)指標(biāo)的城市人水和諧度評(píng)價(jià)指標(biāo)體系[3](見(jiàn)表 1),其中C2、C10、C11、C12為定性指標(biāo),采用問(wèn)卷調(diào)查的方法確定;其余為定量指標(biāo),通過(guò)國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展公報(bào)等大數(shù)據(jù)資料確定。

表1 城市人水和諧度評(píng)價(jià)指標(biāo)體系

3 城市人水和諧度評(píng)價(jià)模型

城市人水和諧度評(píng)價(jià)是一個(gè)多指標(biāo)、高維度的問(wèn)題,其評(píng)價(jià)是通過(guò)數(shù)學(xué)的方法,利用已有的指標(biāo)數(shù)據(jù),建立相應(yīng)的評(píng)價(jià)模型,利用該模型與樣本數(shù)據(jù),對(duì)樣本的人水和諧度進(jìn)行評(píng)價(jià)。

主成分分析法是一種分析處理高維度數(shù)據(jù)的方法,它在分析和處理高維度的數(shù)據(jù)時(shí)可以極大程度降低原始數(shù)據(jù)維數(shù)。主成分分析法以其處理高維數(shù)據(jù)以及避免人為主觀成分、實(shí)現(xiàn)客觀評(píng)價(jià)的特點(diǎn),在多種評(píng)價(jià)中得到廣泛使用。因此本文采用基于Spass軟件的主成分分析法構(gòu)建城市人水和諧度評(píng)價(jià)模型,對(duì)城市人水和諧等級(jí)做出客觀評(píng)價(jià)。

3.1 主成分分析法基本思想

主成分分析法的原理即根據(jù)各指標(biāo)特征值以及指標(biāo)方差累計(jì)貢獻(xiàn)率確定指標(biāo)主成分,利用主成分代替原來(lái)的評(píng)價(jià)指標(biāo),從而降低了數(shù)據(jù)維度,使評(píng)價(jià)更加清晰、簡(jiǎn)便。通過(guò)對(duì)低維度空間中數(shù)據(jù)的研究,充分反映出高維數(shù)據(jù)的特點(diǎn),進(jìn)而對(duì)各個(gè)樣本進(jìn)行評(píng)價(jià)對(duì)比。

3.2 主成分分析法模型構(gòu)建過(guò)程

使用主成分分析法構(gòu)建城市人水和諧度評(píng)價(jià)模型,應(yīng)按以下步驟[11]進(jìn)行:

(1)步驟1:對(duì)評(píng)價(jià)指標(biāo)的原始數(shù)據(jù)進(jìn)行歸一化處理。

設(shè)共有m個(gè)評(píng)價(jià)指標(biāo),k個(gè)待評(píng)價(jià)樣本,列出待評(píng)樣本的指標(biāo)原始值矩陣,為了避免各個(gè)指標(biāo)因其單位、數(shù)量級(jí)等不同而造成數(shù)據(jù)操作以及數(shù)據(jù)對(duì)比的不便,可以采用公式(1)對(duì)評(píng)價(jià)指標(biāo)的原始數(shù)據(jù)進(jìn)行歸一化處理。

(1)

(2)步驟2:確定主成分個(gè)數(shù)。

對(duì)歸一化數(shù)據(jù)進(jìn)行主成分分析,當(dāng)成分特征值大于1且方差累計(jì)貢獻(xiàn)率大于85%時(shí),確定主成分個(gè)數(shù)p,方差累計(jì)貢獻(xiàn)率采用公式(2)確定。

(2)

式中:α為方差累計(jì)貢獻(xiàn)率;p為主成分個(gè)數(shù);m為指標(biāo)總個(gè)數(shù);αi為第i個(gè)指標(biāo)的方差貢獻(xiàn)率。

(3)步驟3:確定主成分表達(dá)式。

對(duì)于第k個(gè)評(píng)價(jià)樣本有:

(3)

式中:Fki為第k個(gè)樣本的第i個(gè)主成分表達(dá)式;bij為第i個(gè)主成分中,第j個(gè)指標(biāo)的權(quán)重,可以用第i個(gè)主成分與第j個(gè)指標(biāo)的相關(guān)系數(shù)除以主成分i的特征值開(kāi)平方根求得;xjk為第k個(gè)樣本中第j個(gè)指標(biāo)歸一化值。

(4)步驟4:建立主成分評(píng)價(jià)模型。

確定第k個(gè)樣本的評(píng)價(jià)函數(shù)Fk。

Fk=a1Fk1+a2Fk2+…+apFkp

(4)

通過(guò)與樣本結(jié)合,得到評(píng)價(jià)值,評(píng)價(jià)值越大說(shuō)明樣本越優(yōu):

F*=max(Fk)

(5)

式(4)、(5)中:Fk為第k個(gè)樣本的評(píng)價(jià)函數(shù);ap為第p個(gè)主成分方差貢獻(xiàn)率,可以通過(guò)第p個(gè)主成分方差與主成分方差累積之商表示;F*為最優(yōu)評(píng)價(jià)值,為所有樣本評(píng)價(jià)值的最大值。

4 實(shí)例研究

北京市是我國(guó)歷史文化古都,文化、政治、商業(yè)、科技發(fā)展中心,位于華北平原的西北邊緣,靠近天津并且被河北環(huán)繞。全市土地面積約為16 800km2,轄區(qū)內(nèi)有16個(gè)區(qū)、147個(gè)街道、38個(gè)鄉(xiāng)、144個(gè)鎮(zhèn)。北京市內(nèi)無(wú)天然湖泊,大、小水庫(kù)共85座,有五大自然水系匯入渤海。

北京市近年來(lái)發(fā)展迅速,人口激增,致使水荒現(xiàn)象不斷加重,城市發(fā)展對(duì)水的需求遠(yuǎn)遠(yuǎn)超過(guò)了水資源供給能力。多年來(lái),通常以犧牲地下水環(huán)境來(lái)保證城市的發(fā)展,地下水連年超采。北京市作為嚴(yán)重缺水的城市之一,人均水資源量遠(yuǎn)遠(yuǎn)達(dá)不到300m3,僅為我國(guó)人均水資源量的1/8,是全世界人均水資源占有量的1/30,因此對(duì)于北京市的人水和諧度應(yīng)該引起極大的重視。通過(guò)查閱北京市2011-2016年國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào),將北京市連續(xù)6年評(píng)價(jià)指標(biāo)原始數(shù)據(jù)匯總?cè)绫?所示。

表2 北京市2011-2016年評(píng)價(jià)指標(biāo)原始數(shù)據(jù)

注:表格中數(shù)據(jù)來(lái)自2011-2016年北京市國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)。

4.1 城市人水和諧度評(píng)價(jià)模型應(yīng)用

對(duì)原始數(shù)據(jù)進(jìn)行整理,根據(jù)前述模型構(gòu)建步驟,利用Spass軟件[12]計(jì)算得各成分初始特征值以及解釋的方差貢獻(xiàn)率,整理見(jiàn)表3。

表3 初始特征值及解釋的方差貢獻(xiàn)率

由表3可知,共有3個(gè)滿足特征值大于1且方差貢獻(xiàn)累計(jì)為93.05%的成分,由前文方差累計(jì)貢獻(xiàn)率大于85%的原則可知,這3個(gè)成分可以很好地反映原來(lái)19個(gè)評(píng)價(jià)指標(biāo)。所以根據(jù)主成分計(jì)算過(guò)程,選擇3個(gè)新的變量作為主成分,主成分矩陣整理見(jiàn)表4。

由表4可以看出,第1主成分在“生態(tài)發(fā)展情況”層面所包含指標(biāo)上的載荷極大,分別為0.978,0.978,0.987,0.989,0.990,0.995,所以F1主成分可以很好替代“生態(tài)發(fā)展情況”層面的指標(biāo)。第2主成分在“水資源情況”層面所包含指標(biāo)上的荷載極大,分別為0.902,0.719,0.925,0.902,F(xiàn)2可以較好替代“水資源情況”層面的指標(biāo)。第3主成分在“經(jīng)濟(jì)、社會(huì)發(fā)展情況”層面所包含指標(biāo)上的載荷較大,分別為0.277,0.027,0.152,0.768,0.398,所以F3主成分可以較好替代“經(jīng)濟(jì)、社會(huì)發(fā)展情況”層面的指標(biāo)。

主成分中各指標(biāo)得分結(jié)果整理見(jiàn)表5。

表4 主成分矩陣

表5 主成分得分系數(shù)

由前述計(jì)算過(guò)程與表3、表5可得北京市人水和諧度評(píng)價(jià)模型:

F1=0.016X1+0.080X2+0.008X3-0.038X4+

0.006X5+0.016X6+0.082X7+0.082X8+

0.083X9+0.083X10+0.083X11+0.083X12+

0.083X13-0.081X14-0.075X15-0.083X16+

0.079X17-0.035X18+0.027X19

F2=0.214X1-0.059X2+0.170X3-0.126X4+

0.219X5+0.214X6-0.039X7+0.028X8+

0.010X9-0.007X10-0.029X11-0.002X12-

0.007X13+0.039X14-0.073X15-0.017X16-

0.062X17-0.004X18+0.186X19

F3=0.189X1+0.077X2-0.359X3-0.167X4+

0.182X5+0.189X6+0.019X7-0.006X8-

0.099X9-0.016X10-0.013X11-0.017X12+

0.038X13-0.107X14+0.181X15+0.018X16+

0.100X17+0.052X18+0.260X19

F=0.6747F1+0.2388F2+0.0865F3

將北京市2011-2016年的指標(biāo)歸一化數(shù)值代入該模型,得2011-2016年的綜合得分依次為:-0.9466,-0.1141,-0.3251,-0.2154,0.4491,1.1521。 2011-2016年北京市人水和諧度結(jié)果匯總見(jiàn)表6。

表6 2011-2016年各項(xiàng)得分及排名

4.2 評(píng)價(jià)結(jié)果分析

4.2.1 計(jì)算結(jié)果分析 根據(jù)表6中各項(xiàng)分?jǐn)?shù)以及排名繪制“各項(xiàng)得分-年份”折線圖(圖1),從圖1中可以直觀看出:

(1)2011-2016年北京市人水和諧度總體呈現(xiàn)上升趨勢(shì),其中2012年水資源總量為近幾年最大,達(dá)到39.5×108m3,所以2012年排在第3名,其余年份排名呈現(xiàn)持續(xù)上升趨勢(shì),表明由于最嚴(yán)格水資源管理措施的實(shí)施以及2014年南水北調(diào)南水正式入京,北京市人水和諧程度逐漸得到改善。

(2)由F1主成分排名可以看出北京市生態(tài)建設(shè)情況逐年上升,表明北京市政府采取的一系列生態(tài)治理措施如減少污染排放、促進(jìn)低碳發(fā)展、改善生態(tài)環(huán)境質(zhì)量、提升城市生態(tài)環(huán)境質(zhì)量和承載能力等收到了極大的成效,但在城市綠化率、水文化建設(shè)等方面仍有較大的提升空間。

(3) 由F2主成分排名可以看出, 2012年降水較多,水資源量豐富,人均水資源量最高為190.89m3,因此其水資源情況最優(yōu),排名第1;2013年北京市水資源總量較2012年下降較大,但其人口增長(zhǎng)率較高,因此使其水資源情況出現(xiàn)大幅度下降;2014年人均水資源量最低,為100.4m3,由于北京市加強(qiáng)了水資源管理、實(shí)施水務(wù)統(tǒng)管、開(kāi)展污水回用、雨洪利用、鼓勵(lì)市民節(jié)約用水以及南水北調(diào)南水入京,北京市水資源情況從2014年開(kāi)始有了大幅度改善。

(4)由F3主成分排名看出,北京市社會(huì)經(jīng)濟(jì)發(fā)展情況發(fā)展不夠理想,人口增長(zhǎng)較快,人口密度大,2016年人口密度為1 324人/km2,對(duì)人水和諧造成極大威脅。北京市應(yīng)嚴(yán)格控制城市人口增長(zhǎng),加速產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型,形成以第三產(chǎn)業(yè)為主導(dǎo)的現(xiàn)代化節(jié)水型城市。綜合4項(xiàng)排名來(lái)看,2016年的4項(xiàng)得分中,F(xiàn)3主成分排名較低,所以今后北京市應(yīng)重點(diǎn)從社會(huì)經(jīng)濟(jì)層面改善其人水和諧狀況,促進(jìn)“人-水-經(jīng)濟(jì)-社會(huì)-生態(tài)”五位一體化發(fā)展。

圖1 各項(xiàng)得分-年份折線圖

4.2.2 與投影尋蹤法評(píng)價(jià)結(jié)果對(duì)比分析 投影尋蹤法[13-14]也是高維度問(wèn)題的一種有效處理方法,其基本原理就是通過(guò)最大化投影指標(biāo)函數(shù)來(lái)確定最優(yōu)的投影方向,從而將高維度數(shù)據(jù)投影到低維度子空間上,通過(guò)對(duì)低維度空間中數(shù)據(jù)的分析,使研究結(jié)果更加清晰,從而達(dá)到對(duì)高維原始數(shù)據(jù)進(jìn)行研究的目的。通過(guò)投影尋蹤法計(jì)算,最佳投影方向?yàn)?0.2644, 0.2579, 0.1361, 0.2601, 0.2785, 0.2466, 0.2077, 0.2941, 0.2468, 0.1524, 0.2682, 0.1377, 0.2534, 0.2135, 0.3117, 0.2321, 0.0928, 0.2467, 0.0634),2011年-2016年投影值依次為:0.4477,1.4770,1.4766,1.4777,2.7881,3.5444。北京市2011-2016年人水和諧度從高到低依次為:2016、2015、2014、2012、2013、2011年。主成分分析法評(píng)價(jià)結(jié)果與投影尋蹤法評(píng)價(jià)結(jié)果基本一致,結(jié)果見(jiàn)表7,評(píng)價(jià)結(jié)果與北京市基本情況相符,基于主成分分析法的人水和諧度評(píng)價(jià)切實(shí)可行。

表7 2011-2016年兩種方法評(píng)價(jià)結(jié)果排名對(duì)比

5 結(jié) 論

通過(guò)對(duì)城市的人水和諧度進(jìn)行評(píng)價(jià),可以反映出城市人水問(wèn)題,本文利用主成分分析法對(duì)城市人水和諧度問(wèn)題做出了客觀評(píng)價(jià),取得以下兩點(diǎn)結(jié)論:

(1)本文立足于北京市的人水關(guān)系現(xiàn)狀,從城市水資源情況、生態(tài)發(fā)展情況、經(jīng)濟(jì)發(fā)展情況、社會(huì)發(fā)展情況4個(gè)方面提出了包含19項(xiàng)評(píng)價(jià)指標(biāo)的城市人水和諧度評(píng)價(jià)指標(biāo)體系。

(2)利用主成分分析法提出了城市人水和諧度評(píng)價(jià)模型的構(gòu)建方法,可避免人為主觀因素對(duì)評(píng)價(jià)結(jié)果的影響,并利用該模型對(duì)北京市2011-2016年的人水和諧度做出評(píng)價(jià),結(jié)果與投影尋蹤法所得結(jié)果基本一致,說(shuō)明了將主成分分析法與城市人水和諧評(píng)價(jià)問(wèn)題結(jié)合的可行性,并根據(jù)評(píng)價(jià)結(jié)果,為北京市人水和諧建設(shè)提出了改進(jìn)措施。

[1] 左其亭,張 云.人水和諧量化研究方法及應(yīng)用[M].北京:中國(guó)水利水電出版社, 2009.

[2] 汪恕誠(chéng).人與自然和諧相處-破解中國(guó)水問(wèn)題的核心理念[J].今日國(guó)土,2004(Z2):6-9.

[3] 左其亭,張 云,林 平.人水和諧評(píng)價(jià)指標(biāo)及量化方法研究[J].水利學(xué)報(bào),2008,39(4):440-447.

[4] 戴會(huì)超,唐德善,張范平,等.城市人水和諧度研究[J].水利學(xué)報(bào),2013,44(8):973-978+986.

[5] 張金鑫,唐德善,丁億凡,等.基于云模型的流域人水和諧評(píng)價(jià)方法[J].水電能源科學(xué), 2015,33(9):155-158+127.

[6] 賈喜午,周 堅(jiān),呂慶云,等.基于主成分分析法的復(fù)合營(yíng)養(yǎng)米評(píng)價(jià)模型的建立[J].中國(guó)糧油學(xué)報(bào),2015, 30(7):123-127.

[7] 郭 瀟,方國(guó)華,章哲愷.跨流域調(diào)水生態(tài)環(huán)境影響評(píng)價(jià)指標(biāo)體系研究[J].水利學(xué)報(bào), 2008, 39(9):1125-1130+1135.

[8] 趙洪杰,唐德善.黑河中游灌區(qū)水資源綜合效益評(píng)價(jià)研究[J].節(jié)水灌溉, 2006(6):58-60+62.

[9] 馮彩云,許新橋,孫振元.北京近自然園林綠地植物群落綜合評(píng)價(jià)指標(biāo)體系研究[J].安徽農(nóng)業(yè)大學(xué)學(xué)報(bào), 2014, 41(6):950-955.

[10] 孔 越,陳 娟.社區(qū)衛(wèi)生服務(wù)滿意度指標(biāo)體系中專家咨詢法的可靠性分析[J].解放軍預(yù)防醫(yī)學(xué)雜志,2007, 25(4):259-261.

[11] 王林鵬. 基于主成分分析法的深基坑支護(hù)方案優(yōu)選模型[J].河南科學(xué),2011,29(6):699-702.

[12] 楊 競(jìng), 童禎恭, 劉玉哲.SPSS軟件對(duì)飲用水水質(zhì)進(jìn)行主成分分析評(píng)價(jià)的運(yùn)用[J]. 環(huán)境科學(xué)與技術(shù), 2011, 34(7):171-174.

[13] 黃顯峰,賈永樂(lè),方國(guó)華.基于投影尋蹤法的城市水生態(tài)文明建設(shè)評(píng)價(jià)[J].水資源保護(hù),2016,32(6):117-122.

[14] 方國(guó)華. 多目標(biāo)決策理論、方法及其應(yīng)用[M].北京:科學(xué)出版社, 2011.

猜你喜歡
評(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à)
主站蜘蛛池模板: 美女无遮挡被啪啪到高潮免费| 亚洲A∨无码精品午夜在线观看| 欧美综合成人| 99伊人精品| 日本www在线视频| 日韩欧美国产中文| 国产AV无码专区亚洲A∨毛片| 欧美亚洲日韩中文| 日本精品αv中文字幕| 国产福利影院在线观看| 不卡国产视频第一页| 日本不卡视频在线| 国产精品林美惠子在线播放| 99久久精品视香蕉蕉| 久久亚洲国产最新网站| 久久久久国色AV免费观看性色| 超清无码熟妇人妻AV在线绿巨人| m男亚洲一区中文字幕| 国产在线精彩视频二区| 波多野衣结在线精品二区| 亚洲国产天堂在线观看| av尤物免费在线观看| 国产精品污视频| 夜精品a一区二区三区| 四虎成人精品在永久免费| 无码中文AⅤ在线观看| 欧美午夜小视频| 青青极品在线| 一级毛片免费观看不卡视频| 97视频精品全国在线观看| 国产中文在线亚洲精品官网| 制服丝袜在线视频香蕉| 色欲综合久久中文字幕网| 亚洲欧美激情小说另类| 毛片视频网址| 免费久久一级欧美特大黄| 精品无码视频在线观看| 一区二区三区四区日韩| 日本久久久久久免费网络| 毛片基地视频| 国产又粗又爽视频| 亚洲视频无码| 欧美在线一级片| 久久国产精品影院| 国产色伊人| 一本大道香蕉久中文在线播放 | 久久亚洲天堂| 亚洲欧洲日韩综合| 欧美日韩一区二区三区四区在线观看| 中文一级毛片| 色综合天天娱乐综合网| 又大又硬又爽免费视频| 久久久久免费精品国产| 色偷偷一区二区三区| 就去吻亚洲精品国产欧美| 国产成人免费观看在线视频| 午夜毛片免费观看视频 | 国产麻豆福利av在线播放| 国产精品短篇二区| 国产乱人视频免费观看| 午夜不卡视频| 久草青青在线视频| 99视频在线免费观看| 久久国产乱子| 99精品免费在线| 国产精品页| 91福利免费视频| AV天堂资源福利在线观看| 国产一区二区精品福利| 国产成人在线无码免费视频| 亚洲综合激情另类专区| 亚洲国产亚综合在线区| 国产午夜精品鲁丝片| 中文国产成人精品久久| 国产SUV精品一区二区| 亚洲性一区| 国产免费自拍视频| 视频国产精品丝袜第一页| 成人综合网址| 日韩成人在线视频| 久久国产精品影院| 国产免费人成视频网|