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

    空間插值法在熱環(huán)境流動觀測中的應用

    2012-11-30 02:33:14楊永川楊軻王志浩盧軍
    中南大學學報(自然科學版) 2012年9期
    關鍵詞:冷島插值法克里

    楊永川,楊軻,王志浩,盧軍

    (1.重慶大學 城市建設與環(huán)境工程學院,重慶,400045;2.重慶大學 三峽庫區(qū)生態(tài)環(huán)境教育部重點實驗室,重慶,400045)

    1 研究方法

    1.1 數(shù)據(jù)來源

    本研究所用的溫度是2009年9月9日在重慶都市區(qū)華巖公園進行流動觀測所測得的 3組溫度(觀測時間分別為 9∶00,11∶00 和 13∶00),流動范圍約 1.5 km×2.5 km,每組溫度有500個左右。

    1.2 插值方法

    1.2.1 反距離加權插值法

    反距離加權插值法是最常用的插值方法之一,其依據(jù)是空間數(shù)值與距離成反比的空間特性,公式為:

    式中:Z(x0)為x0點處的估算值;N為用于插值的鄰近點的個數(shù);Z(xi)為樣點xi處的實測值;λi為第i個樣點對估值點貢獻的權重,它由下式計算[12]:

    式中:di0為待估點與樣點間的距離;p為距離的冪,它控制著隨距離增大權系數(shù)的變化,同時,也決定著內(nèi)插結(jié)果的平滑效果。其選擇標準是最小平均誤差。

    1.2.2 普通克里格法(Ordinary Kriging)

    普通克里格法是地統(tǒng)計學中的一種重要的插值方法,它是以區(qū)域變量理論為基礎,以變差函數(shù)為主要工具,在保證估計值滿足無偏性和有效性的前提下求得估計值。設區(qū)域化變量f(x)滿足二階平穩(wěn)假設或本征假設,則待插點 P的估計值為(式中fi是n個已知點的函數(shù)值,ω是n個已知點的全系數(shù),并滿足:再根據(jù)估計的方差最小的條件:∑ωiγ(xj?xi)+μ=γ(xp?xi)(i=1,2,…,n)(其中,μ 為拉格朗日算子,γ(xj?xi)為已知點間的變差函數(shù)值,γ(xp?xi)為已知點與待插點間的變差函數(shù)值。)求出待插點P的估計值[13]。簡單地說,克里格插值法就是一種特定的滑動加權平均法。

    在臨床實踐中發(fā)現(xiàn),由于疾病和透析引起的焦慮抑郁情緒使得透析患者囿于家中,懶于活動,并且不愿意與社會聯(lián)系。本研究結(jié)果顯示,動機性訪談可以促使患者積極參與運動,適當?shù)倪\動不僅可以改善患者的身體狀態(tài),而且可以對其心理狀態(tài)產(chǎn)生積極影響,從而改善患者的心境和情緒,患者的社會活動和參與感增強,生活能力顯著提高,生活質(zhì)量明顯改善[16]。而生活質(zhì)量提高的PD患者,其運動依從性也明顯好于對照組,原因可能是生活質(zhì)量提高的腹膜透析患者,他們的身體和心理狀態(tài)都處于較高水平,更易于接受運動訓練。

    普通克里格法根據(jù)擬合半方差模型的不同可分為球形模型(Spherical)、環(huán)狀模型(Circular)、指數(shù)模型(Exponential)、高斯模型(Gaussian)和線性模型(Linear)。普通克里格法使用數(shù)學函數(shù)與方法指定參數(shù)以適應在不同數(shù)據(jù)中的不同直線或曲線,普通克里格法假設數(shù)據(jù)的任何部分都是自由變化的,利用這5種模型,以確保提供的變化可以滿足數(shù)據(jù)的變化,但這并不可能總是與實際的數(shù)據(jù)相適應。

    1.2.3 樣條函數(shù)法(Spline)

    Spline是在空間插值時準確地通過實測樣點擬合出連續(xù)光滑表面。其表達式如下[14]:

    式中:Z為要預測的值;n為參與插值的實測點數(shù);λi為一系列線性方程解所確定的系數(shù);γi為估測點到第i點的距離;R(γi)和T(x,y)表達式如下:

    式中:τ2為權重系數(shù);γ為已知點與采樣點之間的距離;k0為改正后的貝塞爾函數(shù);c為常數(shù);a為線性方程的系數(shù)。

    1.3 數(shù)據(jù)分析

    首先對流動觀測所得到的溫度進行誤差修正和時間修正[11],將所測溫度按流動觀測路線排列,從排列好的溫度中每隔2行提出1行(約33.3%)用于驗證插值結(jié)果,用剩下的(約 66.6%)溫度進行插值,插值在ARCGIS9.2中完成,再利用空間分析中的柵格計算工具在插值結(jié)果圖層和驗證點圖層之間進行乘法運算提取出驗證點的插值結(jié)果,并對驗證點的插值結(jié)果圖層和驗證點的實測圖層進行比較,計算均絕對誤差(MAE)、平均相對誤差(MRE)和平均標準差(RMSIE),考察使用不同的插值方法時3種誤差的變化情況。通過分析誤差確定在流動觀測溫度中的最佳插值方法。

    2 結(jié)果分析

    分別對3組66.6%的溫度進行插值,得到的結(jié)果有所不同,圖1所示為上午9∶00的溫度插值結(jié)果,圖2所示為3組溫度在3種插值方法不同參數(shù)下的誤差計算結(jié)果。下面對 3種插值方法結(jié)合溫度進行具體分析。

    2.1 反距離加權法的分析

    在3組數(shù)據(jù)中,IDW的平均相對誤差和平均絕對誤差隨著冪指數(shù)的增加而減小,但這并不能說明IDW法的插值結(jié)果就是隨著冪指數(shù)的升高而減小,因為平均標準差同時在增大,這說明估值的靈敏度和極值效應隨著冪指數(shù)的升高而增加。另外,從IDW插值的原理來看,隨著冪指數(shù)增大,較遠的點對插入的數(shù)據(jù)影響越小。由于所采用的數(shù)據(jù)驗證點非常接近,驗證點受鄰近點的影響很大,所以,IDW的插值結(jié)果中平均絕對誤差和平均相對誤差在一定范圍內(nèi)必然會隨著冪指數(shù)的升高而減小,但是,冪指數(shù)的升高并不能保證在插值范圍內(nèi)所有值的準確性,因此,在數(shù)據(jù)密度較小或局部數(shù)據(jù)密度遠大于總體數(shù)據(jù)密度時不建議采用IDW法,在實測數(shù)據(jù)分布均勻且密度較大的情況下采用 IDW 法較好,但要根據(jù)數(shù)據(jù)的具體情況驗證出 1個合適的冪指數(shù)以保證插值結(jié)果的準確性。用IDW法在插值時是在一定的搜索半徑內(nèi)的數(shù)據(jù)加權平均得到未知點的數(shù)據(jù),當每個點的數(shù)據(jù)變化較小時,插值結(jié)果變化也相對較小,只有在采樣點屬性變化劇烈和頻繁的地方才會出現(xiàn)較大的變化。

    綜上所述,用IDW對GPS+溫濕度自動記錄儀這種流動觀測得到的溫度進行插值,可以得到較好的插值結(jié)果,但是,如果觀測點不均勻理論上會導致局部插值結(jié)果嚴重失真。這種方法可以作為一種簡單快速的方法使用。

    2.2 普通克里格(Ordinary Kriging)插值法的分析

    從圖2中的3項誤差分析可以看出在不同溫度中普通克里格的各種模型表現(xiàn)的效果有所不同,相比較來說指數(shù)模型的各項誤差值始終為最小,其他幾種模型的驗證參數(shù)的誤差在3組溫度中表現(xiàn)不穩(wěn)定,高斯模型的3項誤差都較大,在2組溫度中環(huán)狀模型和球形模型誤差略小于高斯模型的誤差,但在第2組溫度中最大,因此,這2種模型的誤差表現(xiàn)不穩(wěn)定,不建議選用。從3組溫度的各種誤差參數(shù)的平均值來看指數(shù)模型的標準差、相對誤差、絕對誤差均為最小,如圖3所示。

    相對于其他2種插值方法,普通克里格插值法的各項參數(shù)雖然不是最好的,但普通克里格插值法中指數(shù)模型在3組溫度中各項誤差總處于穩(wěn)定的位置,因此,采用普通克里格插值法中指數(shù)模型來對 GPS+溫度濕度自動記錄儀的流動觀測溫度進行插值可以保證插值結(jié)果的精度。

    普通克里格插值法是利用半方差圖的模型進行測定點之間的內(nèi)插,該方法已經(jīng)廣泛用于各種地理數(shù)據(jù)處理。劉加平等[15]在西安冬季城市熱島調(diào)查研究時使用了普通克里格插值法;李軍龍等[16]在對全國 2 114個氣象站的氣象數(shù)據(jù)進行交叉分析驗證時也發(fā)現(xiàn)普通克里格插值法對溫度插值來說是首選的插值方法。

    2.3 樣條函數(shù)(Spline)插值法的分析

    3組數(shù)據(jù)在不同插值方法下的各種誤差平均值比較值見圖3。從圖3可以看出:樣條函數(shù)插值法3個驗證參數(shù)變化沒有規(guī)律,時大時小,這主要是由其本身的算法決定的。樣條函數(shù)插值法是用多項式擬合的方法產(chǎn)生平滑的插值曲線,由于多項式的階數(shù)較低,對溫度誤差的響應不敏感,具有較好的保凸性、逼真性和平滑性,但這種變形會對原始溫度空間插值產(chǎn)生過于理想化的偏差且這個偏差不能估計。

    李軍龍等[16]在研究中指出樣條函數(shù)插值法可以較好地反映高山區(qū)域的年均溫度空間分布特征,在地形復雜的山區(qū)使用樣條函數(shù)插值法較好。由于流動觀測只是用在微氣候的觀測上,觀測范圍很小,因此,在對溫度插值時用樣條函數(shù)插值法的意義不大。

    圖1 上午9:00溫度在不同插值法下的插值結(jié)果Fig.1 Results of 9:00 am data interpolated in different methods

    圖2 3組溫度在不同插值方法下的誤差比較Fig.2 Comparison of accuracy about three groups of data in different interpolation methods

    圖3 3組溫度在不同插值方法下的各種誤差平均值比較Fig.3 Comparison of the average of accuracy about the three groups of data in different interpolation method

    3 重慶都市區(qū)華巖公園的熱環(huán)境分析評價

    冷島效應是相對于城市熱島效應提出的一個概念,在水體和綠地覆蓋率高的地區(qū),溫度要比周圍溫度低,這種現(xiàn)象被稱為冷島效應(Cool-island effect)[17?19]。Chang等[20]在 2007年對臺灣的 61個都市區(qū)公園的冷島效應進行了觀測;李東海等[21]在2008年利用遙感技術分析了河流的冷島效應影響范圍。本研究以華巖公園的冷島效應觀測溫度為例對普通克里金中指數(shù)模型的插值結(jié)果進行了簡單分析。

    華巖公園上午9∶00的溫度場(見圖4)。從圖4可以看出:在整個觀測范圍內(nèi)最低溫度出現(xiàn)在公園的中心位置為34.5 ℃,最高溫度出現(xiàn)在觀測范圍的南部(公園的南偏東方向)為 36.3 ℃。公園南部的溫度場隨離公園距離的增加均勻升高,在距公園1.5 km的距離內(nèi)溫度升高了1.5 ℃左右,然而,在公園北部、東部和西部溫度變化梯度較大,在不到0.5 km內(nèi)升高了1.5℃左右,可以看出公園內(nèi)溫度要比其周圍的都市區(qū)環(huán)境溫度明顯要低。由該溫度場可以看出公園在對四周的冷島效應并不相等,在公園的南部溫度梯度要明顯比其他幾個方向小的多,主要是因為此時有北偏西小于1.5 m/s的微風。

    圖4 華巖公園上午9∶00的溫度場Fig.4 Temperature field obtained at 9:00 am in Huayan Park

    華巖公園上午11∶00的溫度場(見圖5)。從圖5可以看出:在整個觀測范圍內(nèi)最低溫度出現(xiàn)在公園的中心位置為37.1 ℃,最高溫度出現(xiàn)在觀測范圍的地南部(公園的南偏東方向)為38.2 ℃。此時,公園外的風速小于1.5 m/s,風向為北偏西。公園南部的溫度場隨離公園的距離的增大均勻升高,在距公園1.5 km的距離內(nèi)溫度升高不到1.5 ℃。然而,在公園北部、東部和西部,溫度變化梯度較大,公園對周圍的冷島效應范圍非常有限,在不到 300 m的范圍內(nèi)溫度就達到38.5 ℃,溫升超過1.5 ℃。此時,冷島效應受風向的影響非常明顯,并且比 9∶00的冷島效應的輻射范圍大。

    圖5 華巖公園上午11∶00的溫度場Fig.5 Temperature field obtained at 11:00 am in Huayan Park

    華巖公園13∶00的溫度場(見圖6)。從圖6可以看出:在整個觀測范圍內(nèi),最低溫度出現(xiàn)在公園的中心位置為37.5 ℃,最高溫度出現(xiàn)在觀測范圍的地西部(公園的正西方向)為 40.2 ℃,此時,公園外的風速小于0.5 m/s,風向不穩(wěn)定。公園南部的溫度場從西到東均勻升高,在距離公園0.3 km的距離內(nèi)溫度升高1.2 ℃,在公園東部距離公園0.5 km處溫度升高了1.3 ℃,在公園西部距離公園0.3 km處溫度升高了1.9 ℃,在公園北部距離公園0.3 km處溫度升高了1.9 ℃,因此,此時公園對周圍的冷島效應范圍非常有限。

    從上述3個溫度場來看,上午9∶00公園下風向處溫度呈均勻升高趨勢,但此時受太陽輻射影響較小,在人口密集區(qū)溫度受人為排熱的影響較大;在上午11∶00,公園下風向溫度仍呈均勻升高趨勢,在1.5 km時溫度升高不到1.5 ℃,但溫度升高的梯度比9∶00的溫度場要小;13∶00時的溫度場可以看出在1.5 km的范圍內(nèi)沿風向方向(風速小于 0.5 m/s)溫度升高不到1.0 ℃。因此,該公園的冷島效應在1.5 km的范圍內(nèi)13∶00時要比9∶00和11∶00的要大,但在人口密集區(qū)受人為排熱的影響,公園冷島效應對周圍環(huán)境的冷卻效應有限。

    4 結(jié)論

    (1)對于 GPS+溫濕度自動記錄儀的流動觀測溫度的插值來說,克里格插值法中指數(shù)模型是首選的插值方法;克里格插值法中指數(shù)模型對于這種類型的溫度進行插值的結(jié)果較為理想,且精度穩(wěn)定,而其他模型的插值結(jié)果的精度不太穩(wěn)定。

    (2)IDW法適用于測點密度大且分布均勻的數(shù)據(jù)處理,而這種類型的數(shù)據(jù)分布不均勻,使用IDW法不能保證其精度,而且要根據(jù)數(shù)據(jù)的密度通過驗證來確定一個合適的冪指數(shù)。

    (3)樣條插值法在 GPS+溫濕度自動記錄儀的流動觀測溫度中插值的精度不穩(wěn)定且不能估計,因此,使用該方法的意義不大。

    (4)由 3個用普通克里金插值法中指數(shù)模型得到的溫度場分析得出:當風速小于0.5 m/s時,在距公園0.3 km的人口密集區(qū)溫度比公園內(nèi)溫度高2.0 ℃;冷島效應在13∶00時較為明顯;公園的冷島效應范圍受風向和太陽輻射強度的影響較大。

    [1]CHEN Xiao-long,ZHAO Hong-mei,LI Ping-xiang,et al.Remote sensing image-based analysis of the relationship between urban heat island and land use/cover changes[J].Remote Sensing of Environment,2006,104(2):133?146.

    [2]Yuan F,Bauer M E.Comparison of impervious surface area and normalized difference vegetation index as indicators of surface urban heat island effects in Landsat imagery[J].Remote Sensing of Environment,2006,11(3):58?65.

    [3]MA Ya,KUANG Yao-qiu,HUANG Ning-sheng.Coupling urbanization analyses for studying urban thermal environment and its interplay with biophysical parameters based on TM/ETM+imagery[J].International Journal of Applied Earth Observation and Geoinformation,2010,12(2):110?118.

    [4]Imhoff M L,Zhang Ping,Wolfe R E,et al.Remote sensing of the urban heat island effect across biomes in the continental USA[J].Remote Sensing of Environment,2010,114(3):504?513.

    [5]Kolokotroni M,Giannitsaris I,Watkins R.The effect of the London urban heat island on building summer cooling demand and night ventilation strategies[J].Solar Energy,2006,80:383?392.

    [6]Lee S H,Baik J J.Statistical and dynamical characteristics of the urban heat island intensity in Seoul[J].Theoretical and Applied Climatology,2010,100(1/2):227?237.

    [7]Sofer M,Potchter O.The urban heat island of a city in an arid zone:the case of Eilat,Israel[J].Theoretical and Applied Climatology,2006,85(1/2):81?88.

    [8]郭勇,龍步菊,劉偉東,等.北京城市熱島效應的流動觀測和初步研究[J].氣象科技,2006,34(6):656?661.GOU Yong,LONG Bu-ju,LIU Wei-dong,et al.Primary study of urban heat island effect in Beijing using mobile measurements[J].Meteorological Science and Technology,2006,34(6):656?661.

    [9]Lindberg F.Modelling the urban climate using a local governmental geo-database[J].Meteorological Applications,2007,14(3):263?273.

    [10]HUANG Liang-mei,ZHAO De-hua,WANG Jia-zhen,et al.Scale impacts of land cover and vegetation corridors on urban thermal behavior in Nanjing,China[J].Theoretical and Applied Climatology,2008,94(3/4):241?257.

    [11]小澤行雄,吉野正敏.小氣候調(diào)查方法[M].郭殿福譯.南京:廣西人民出版社,1982:40?42.Ozawa Y,Yoshino M(Translated by GUO Dian-fu).Local Climatology Study Method[M].GUO Dian-fu transl.Nanning:Guangxi People Press,1982:40?42.

    [12]Watson D F,Philip G M.A refinement of inverse distance weighted interpolation[J].Geoprocessing,1985,2:315?327.

    [13]Oliver M A.Kriging:A method of interpolation for geographical information systems[J].International Journal of Geographic Information Systems,1990,4:313?332.

    [14]Franke R.Smooth interpolation of scattered data by local thin plate splines[J].Comp &Maths,with Appls,1982,8(4):237?281.

    [15]劉加平,林憲德,劉艷峰,等.西安冬季城市熱島調(diào)查研究[J].太陽能學報,2007,28(8):912?917.LIU Jia-ping,LIN Xian-de,LIU Yan-feng,et al.Survey on winter urban heat island in Xi’an[J].Acta Energiae Solaris Sinica,2007,28(8):912?917.

    [16]李軍龍,張劍,張叢,等.氣象要素空間插值方法的比較分析[J].草業(yè)科學,2006,23(8):6?11.LI Jun-long,ZHANG Jian,ZHANG Cong,et al.Analyze and compare the spatial interpolation methods for climate factor[J].Pratacultural Science,2006,23(8):6?11.

    [17]Hamada S,Ohta T.Seasonal variations in the cooling effect of urban green areas on surrounding urban areas[J].Urban Forestry&Urban Greening,2010,9(1):15?24.

    [18]Marie-Claude F,Konrad G.Holocene climate change and its effect on lake ecosystem production on Northern Victoria Island,Canadian Arctic[J].Journal of Paleolimonology,2010,43(2):219?234.

    [19]Jansson C,Jansson P E,Gustafsson D.Near surface climate in an urban vegetated park and its surroundings[J].Theoretical and Applied Climatology,2007,89:185?193.

    [20]Chang C R,Lib M H,Chang S D.A preliminary study on the local cool-island intensity of Taipei city parks[J].Landscape and Urban Planning,2007,80(4):386?395.

    [21]李東海,艾彬,黎夏.基于遙感和GIS的城市水體緩解熱島效應的研究[J].熱帶地理,2008,28(9):414?418.LI Dong-hai,AI Bin,LI Xia.Urban water Body alleviating heat island effect based on RS and GIS:A case study of Dongguan City[J].Tropical Geography,2008,28(9):414?418.

    猜你喜歡
    冷島插值法克里
    今晚不能去你家玩啦!
    知識窗(2023年12期)2024-01-03 01:38:55
    我可以咬一口嗎?
    知識窗(2023年2期)2023-03-05 11:28:27
    極熱條件下空冷島永磁電動機散熱能力研究
    防爆電機(2022年4期)2022-08-17 06:00:22
    高寒條件下空冷島風機集群永磁直驅(qū)系統(tǒng)批量改造分析
    防爆電機(2022年4期)2022-08-17 05:59:10
    高寒條件下空冷島變頻永磁直驅(qū)改造及應用
    防爆電機(2022年3期)2022-06-17 01:41:20
    你今天真好看
    《計算方法》關于插值法的教學方法研討
    智富時代(2019年7期)2019-08-16 06:56:54
    你今天真好看
    讀者(2018年24期)2018-12-04 03:01:34
    空冷島的熱態(tài)沖洗及防凍措施
    電站輔機(2017年3期)2018-01-31 01:46:41
    基于二次插值法的布谷鳥搜索算法研究
    人妻 亚洲 视频| 色94色欧美一区二区| 欧美激情极品国产一区二区三区| 国产福利在线免费观看视频| 久久久久久亚洲精品国产蜜桃av| 久久 成人 亚洲| 少妇 在线观看| 一个人免费看片子| 国产精品一区二区在线观看99| 欧美精品一区二区免费开放| 精品第一国产精品| 老司机靠b影院| 18禁国产床啪视频网站| 1024视频免费在线观看| 成人影院久久| av免费在线观看网站| 国产高清国产精品国产三级| 午夜精品久久久久久毛片777| 少妇猛男粗大的猛烈进出视频| 99久久精品国产亚洲精品| 久9热在线精品视频| 欧美日韩成人在线一区二区| 别揉我奶头~嗯~啊~动态视频| 一区二区三区精品91| 亚洲国产精品一区二区三区在线| 1024视频免费在线观看| 亚洲天堂av无毛| 免费不卡黄色视频| 成人三级做爰电影| 成人影院久久| 国产1区2区3区精品| 又黄又粗又硬又大视频| 精品卡一卡二卡四卡免费| 老熟妇仑乱视频hdxx| 日韩欧美三级三区| 欧美精品人与动牲交sv欧美| 2018国产大陆天天弄谢| 在线 av 中文字幕| avwww免费| 国产成人精品无人区| 久久久久国产一级毛片高清牌| 亚洲专区字幕在线| 亚洲精品久久成人aⅴ小说| 五月开心婷婷网| 久久香蕉激情| 王馨瑶露胸无遮挡在线观看| 亚洲精品在线观看二区| 狠狠狠狠99中文字幕| 亚洲精品乱久久久久久| 精品福利观看| 一区二区日韩欧美中文字幕| 久久性视频一级片| 黄色毛片三级朝国网站| 无遮挡黄片免费观看| 久久中文字幕一级| 高清av免费在线| 在线十欧美十亚洲十日本专区| 国产男靠女视频免费网站| 亚洲av日韩精品久久久久久密| 最黄视频免费看| 精品第一国产精品| 在线观看www视频免费| 国产精品免费一区二区三区在线 | 少妇精品久久久久久久| 成在线人永久免费视频| 精品久久久久久电影网| 久久精品亚洲精品国产色婷小说| 搡老岳熟女国产| 国产伦人伦偷精品视频| 精品欧美一区二区三区在线| 老司机在亚洲福利影院| 国产伦理片在线播放av一区| 操出白浆在线播放| 午夜成年电影在线免费观看| 日韩视频一区二区在线观看| 久久香蕉激情| 亚洲中文av在线| 亚洲人成77777在线视频| 久久亚洲真实| 精品国产乱码久久久久久小说| 亚洲第一欧美日韩一区二区三区 | 丰满迷人的少妇在线观看| 精品卡一卡二卡四卡免费| 欧美精品亚洲一区二区| 久热爱精品视频在线9| 婷婷成人精品国产| 久久久国产欧美日韩av| 亚洲第一青青草原| 老司机影院毛片| 精品国产乱码久久久久久小说| 久久青草综合色| 999精品在线视频| 国产高清激情床上av| 国产在线视频一区二区| 俄罗斯特黄特色一大片| 国产免费现黄频在线看| 极品教师在线免费播放| 夜夜夜夜夜久久久久| 女人高潮潮喷娇喘18禁视频| 国产精品.久久久| 91老司机精品| 欧美国产精品va在线观看不卡| 妹子高潮喷水视频| 老熟妇仑乱视频hdxx| 蜜桃国产av成人99| 亚洲第一av免费看| 另类精品久久| 久久久国产欧美日韩av| 91麻豆精品激情在线观看国产 | 亚洲精品久久午夜乱码| 久久久久视频综合| 久久久国产一区二区| 天天躁日日躁夜夜躁夜夜| 91九色精品人成在线观看| 在线亚洲精品国产二区图片欧美| 成年人午夜在线观看视频| 免费不卡黄色视频| 午夜精品久久久久久毛片777| 久久午夜亚洲精品久久| 亚洲伊人久久精品综合| 无人区码免费观看不卡 | www.999成人在线观看| 最新美女视频免费是黄的| xxxhd国产人妻xxx| 怎么达到女性高潮| 国产成+人综合+亚洲专区| 中文字幕精品免费在线观看视频| 午夜91福利影院| 国精品久久久久久国模美| 亚洲欧洲日产国产| 无人区码免费观看不卡 | 9热在线视频观看99| 美女午夜性视频免费| 日本欧美视频一区| 久久久精品区二区三区| 午夜精品国产一区二区电影| 99国产极品粉嫩在线观看| 亚洲,欧美精品.| 欧美亚洲 丝袜 人妻 在线| 中文字幕最新亚洲高清| 丝袜在线中文字幕| 日韩欧美一区二区三区在线观看 | 婷婷成人精品国产| 午夜精品国产一区二区电影| 日本五十路高清| 大片免费播放器 马上看| 在线观看舔阴道视频| 又紧又爽又黄一区二区| 国产野战对白在线观看| 精品人妻熟女毛片av久久网站| 国产精品成人在线| 国产高清国产精品国产三级| 亚洲午夜理论影院| 中文亚洲av片在线观看爽 | 亚洲国产看品久久| 黑丝袜美女国产一区| 在线观看免费日韩欧美大片| 老司机福利观看| 丁香六月欧美| 成人特级黄色片久久久久久久 | 久久久久国内视频| 久久久久国产一级毛片高清牌| 香蕉久久夜色| www.精华液| 免费女性裸体啪啪无遮挡网站| 精品国产乱码久久久久久小说| 一级毛片精品| 女性被躁到高潮视频| 欧美 日韩 精品 国产| 两个人免费观看高清视频| 少妇裸体淫交视频免费看高清 | 久久精品国产综合久久久| 97在线人人人人妻| 高清视频免费观看一区二区| 高清毛片免费观看视频网站 | 777久久人妻少妇嫩草av网站| 激情在线观看视频在线高清 | 欧美av亚洲av综合av国产av| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品一区二区www | 亚洲欧美色中文字幕在线| 大陆偷拍与自拍| av福利片在线| 久久青草综合色| 超碰97精品在线观看| 嫩草影视91久久| 欧美人与性动交α欧美精品济南到| 黄片大片在线免费观看| 美女扒开内裤让男人捅视频| 高清欧美精品videossex| 精品亚洲乱码少妇综合久久| 美女高潮到喷水免费观看| 桃红色精品国产亚洲av| 欧美久久黑人一区二区| 手机成人av网站| 在线观看一区二区三区激情| 国产有黄有色有爽视频| 十分钟在线观看高清视频www| 国产99久久九九免费精品| 亚洲三区欧美一区| 国产欧美日韩一区二区三| 国产91精品成人一区二区三区 | 午夜福利一区二区在线看| 亚洲精品美女久久av网站| 成人手机av| 亚洲欧洲精品一区二区精品久久久| av免费在线观看网站| h视频一区二区三区| 国产精品久久久久久人妻精品电影 | 午夜福利视频精品| 久久久久久久大尺度免费视频| 久久久久精品人妻al黑| 成人黄色视频免费在线看| 人人妻人人添人人爽欧美一区卜| 日韩中文字幕欧美一区二区| 亚洲男人天堂网一区| 免费人妻精品一区二区三区视频| 国产成人免费无遮挡视频| 精品一区二区三卡| 久久久久久亚洲精品国产蜜桃av| 久久久久久久国产电影| 高清毛片免费观看视频网站 | 一本—道久久a久久精品蜜桃钙片| 又紧又爽又黄一区二区| 国产成人精品久久二区二区91| 免费在线观看日本一区| kizo精华| 自拍欧美九色日韩亚洲蝌蚪91| 婷婷成人精品国产| 日韩大码丰满熟妇| 美女国产高潮福利片在线看| 人人妻人人澡人人爽人人夜夜| 色尼玛亚洲综合影院| 日韩免费av在线播放| 一级,二级,三级黄色视频| 精品国产国语对白av| 久久av网站| 少妇的丰满在线观看| 999精品在线视频| cao死你这个sao货| 男男h啪啪无遮挡| 两个人免费观看高清视频| 免费不卡黄色视频| 亚洲精品中文字幕在线视频| 男女免费视频国产| 欧美乱码精品一区二区三区| videosex国产| xxxhd国产人妻xxx| 国产伦理片在线播放av一区| 久久久久久亚洲精品国产蜜桃av| www.自偷自拍.com| 国产av又大| 亚洲一区中文字幕在线| 成在线人永久免费视频| 美女福利国产在线| 日本黄色日本黄色录像| 国产高清videossex| 日韩成人在线观看一区二区三区| 精品国产一区二区久久| 老汉色∧v一级毛片| 亚洲黑人精品在线| 夜夜骑夜夜射夜夜干| 99在线人妻在线中文字幕 | 一边摸一边做爽爽视频免费| 欧美精品人与动牲交sv欧美| 一本久久精品| 大香蕉久久网| av国产精品久久久久影院| 亚洲精品久久成人aⅴ小说| 夜夜夜夜夜久久久久| 日韩一区二区三区影片| 69av精品久久久久久 | av超薄肉色丝袜交足视频| 在线观看舔阴道视频| 久久久久精品人妻al黑| 国产成人av激情在线播放| 久久亚洲精品不卡| 欧美日韩亚洲国产一区二区在线观看 | av免费在线观看网站| 香蕉丝袜av| 黄色毛片三级朝国网站| 国产一区二区三区综合在线观看| 国产精品一区二区免费欧美| 久久久欧美国产精品| 国产97色在线日韩免费| 亚洲少妇的诱惑av| 老司机午夜十八禁免费视频| 黑人巨大精品欧美一区二区蜜桃| 99热国产这里只有精品6| 中文亚洲av片在线观看爽 | 国产精品久久电影中文字幕 | 欧美黑人精品巨大| 亚洲精品国产一区二区精华液| 伊人久久大香线蕉亚洲五| 亚洲av国产av综合av卡| 在线观看66精品国产| 性色av乱码一区二区三区2| 亚洲精品av麻豆狂野| 国产精品国产av在线观看| 菩萨蛮人人尽说江南好唐韦庄| 天堂俺去俺来也www色官网| 老鸭窝网址在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲 国产 在线| 成年人黄色毛片网站| 一级黄色大片毛片| 亚洲av欧美aⅴ国产| 久久精品国产综合久久久| 久久久久精品人妻al黑| 看免费av毛片| 另类亚洲欧美激情| av福利片在线| 中文字幕av电影在线播放| 日本一区二区免费在线视频| 亚洲中文日韩欧美视频| 999久久久精品免费观看国产| 亚洲精品美女久久av网站| 99热网站在线观看| 亚洲欧美精品综合一区二区三区| 在线天堂中文资源库| 十分钟在线观看高清视频www| 又紧又爽又黄一区二区| 国产亚洲欧美在线一区二区| 一级片'在线观看视频| 丰满饥渴人妻一区二区三| 大陆偷拍与自拍| 久久久精品94久久精品| 国产欧美日韩一区二区三区在线| 国产免费福利视频在线观看| av天堂久久9| 一个人免费看片子| av又黄又爽大尺度在线免费看| 国产亚洲精品第一综合不卡| 日韩大片免费观看网站| 在线观看一区二区三区激情| 午夜激情久久久久久久| 久久精品aⅴ一区二区三区四区| 777米奇影视久久| 波多野结衣av一区二区av| 极品人妻少妇av视频| 丝瓜视频免费看黄片| 久久 成人 亚洲| 日本黄色日本黄色录像| 国产免费视频播放在线视频| 大香蕉久久网| 色综合婷婷激情| 亚洲中文av在线| 日日摸夜夜添夜夜添小说| 亚洲成人国产一区在线观看| 日本撒尿小便嘘嘘汇集6| 中文字幕高清在线视频| 首页视频小说图片口味搜索| 热99国产精品久久久久久7| 在线永久观看黄色视频| 亚洲国产av影院在线观看| 亚洲色图av天堂| 国产精品九九99| 国产男女内射视频| 狠狠精品人妻久久久久久综合| 69av精品久久久久久 | 国产xxxxx性猛交| 两性夫妻黄色片| 国产精品熟女久久久久浪| 捣出白浆h1v1| 国产激情久久老熟女| 中文字幕制服av| 亚洲黑人精品在线| 亚洲成人手机| 两个人免费观看高清视频| 日本撒尿小便嘘嘘汇集6| 精品国产亚洲在线| 久久中文字幕一级| 国产成+人综合+亚洲专区| 久久人妻福利社区极品人妻图片| 一本大道久久a久久精品| 免费观看av网站的网址| 国产精品久久久久久精品电影小说| cao死你这个sao货| 在线观看免费视频日本深夜| 18在线观看网站| 天堂俺去俺来也www色官网| 国产色视频综合| 亚洲精品成人av观看孕妇| 国产黄色免费在线视频| 蜜桃国产av成人99| 无人区码免费观看不卡 | 色尼玛亚洲综合影院| 国产精品国产av在线观看| 欧美大码av| 美女高潮喷水抽搐中文字幕| 日日摸夜夜添夜夜添小说| 女人久久www免费人成看片| 大片免费播放器 马上看| 99精国产麻豆久久婷婷| 在线 av 中文字幕| 最新的欧美精品一区二区| 亚洲欧美激情在线| 国产av又大| 亚洲专区中文字幕在线| 久久久久久亚洲精品国产蜜桃av| 午夜精品国产一区二区电影| 中文欧美无线码| 激情在线观看视频在线高清 | 午夜成年电影在线免费观看| 久久毛片免费看一区二区三区| 亚洲国产av影院在线观看| 亚洲成人免费电影在线观看| 午夜精品久久久久久毛片777| 亚洲精品在线美女| 久久免费观看电影| 欧美日韩av久久| 欧美精品一区二区大全| 中文字幕高清在线视频| 亚洲成av片中文字幕在线观看| 制服人妻中文乱码| 久久亚洲精品不卡| 美女扒开内裤让男人捅视频| videos熟女内射| 欧美亚洲 丝袜 人妻 在线| 黄色成人免费大全| 精品国产一区二区三区久久久樱花| 99精品在免费线老司机午夜| 国产伦理片在线播放av一区| 午夜久久久在线观看| 啦啦啦免费观看视频1| 人人妻人人澡人人看| 激情视频va一区二区三区| 亚洲精品自拍成人| 777久久人妻少妇嫩草av网站| 51午夜福利影视在线观看| 亚洲国产欧美网| 亚洲av美国av| www.熟女人妻精品国产| 色在线成人网| 日本av免费视频播放| 下体分泌物呈黄色| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦中文免费视频观看日本| 国产91精品成人一区二区三区 | 日韩中文字幕欧美一区二区| 亚洲精品中文字幕一二三四区 | av视频免费观看在线观看| 汤姆久久久久久久影院中文字幕| 丰满少妇做爰视频| kizo精华| 国产精品亚洲av一区麻豆| 999久久久国产精品视频| 满18在线观看网站| 啪啪无遮挡十八禁网站| 在线观看免费高清a一片| 成年人免费黄色播放视频| 中文字幕最新亚洲高清| 久久精品国产亚洲av香蕉五月 | www.精华液| 99re在线观看精品视频| 日本vs欧美在线观看视频| 免费在线观看影片大全网站| 国产成人影院久久av| 亚洲国产欧美日韩在线播放| 日韩大片免费观看网站| 久久毛片免费看一区二区三区| 精品国产乱子伦一区二区三区| 成人亚洲精品一区在线观看| 757午夜福利合集在线观看| 国产成人影院久久av| 欧美日韩亚洲高清精品| 中文字幕色久视频| 精品一品国产午夜福利视频| 精品少妇久久久久久888优播| 大陆偷拍与自拍| 欧美日韩av久久| 香蕉国产在线看| 午夜免费成人在线视频| 国产成人啪精品午夜网站| 宅男免费午夜| 真人做人爱边吃奶动态| 亚洲一码二码三码区别大吗| 精品国产超薄肉色丝袜足j| 美女高潮喷水抽搐中文字幕| 亚洲成人手机| 母亲3免费完整高清在线观看| 国产一区二区三区视频了| 免费不卡黄色视频| 久久毛片免费看一区二区三区| 窝窝影院91人妻| 黄频高清免费视频| 日韩免费高清中文字幕av| 亚洲精品中文字幕一二三四区 | 男女免费视频国产| 国产成+人综合+亚洲专区| 女人被躁到高潮嗷嗷叫费观| 99精国产麻豆久久婷婷| 亚洲av成人一区二区三| 日韩欧美一区视频在线观看| 大片电影免费在线观看免费| 国产真人三级小视频在线观看| 一本色道久久久久久精品综合| 91成人精品电影| 99热网站在线观看| 国产在线一区二区三区精| 久久久精品国产亚洲av高清涩受| 久久人人97超碰香蕉20202| 欧美日韩福利视频一区二区| 老熟妇仑乱视频hdxx| 香蕉国产在线看| 成人特级黄色片久久久久久久 | 久久国产精品人妻蜜桃| 精品乱码久久久久久99久播| 久久国产精品人妻蜜桃| 99精品欧美一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 女性生殖器流出的白浆| 亚洲色图av天堂| 久久久国产成人免费| 国产野战对白在线观看| 最新美女视频免费是黄的| 欧美日韩中文字幕国产精品一区二区三区 | 老汉色∧v一级毛片| 亚洲精品乱久久久久久| 亚洲成人国产一区在线观看| 精品久久久精品久久久| 久久久久久久精品吃奶| 在线观看www视频免费| 性少妇av在线| 丁香欧美五月| 久久青草综合色| 天堂中文最新版在线下载| 精品国产超薄肉色丝袜足j| 久久久久久亚洲精品国产蜜桃av| 黄片小视频在线播放| av电影中文网址| 1024视频免费在线观看| 午夜视频精品福利| 两人在一起打扑克的视频| 欧美变态另类bdsm刘玥| 大片免费播放器 马上看| 午夜福利在线观看吧| 丝瓜视频免费看黄片| 高清在线国产一区| 成在线人永久免费视频| 少妇 在线观看| 视频区图区小说| 亚洲国产毛片av蜜桃av| 少妇的丰满在线观看| 美女午夜性视频免费| av福利片在线| 久久久久久久国产电影| 香蕉国产在线看| 久久久久久人人人人人| 精品熟女少妇八av免费久了| 黄色 视频免费看| 国产黄频视频在线观看| 成人特级黄色片久久久久久久 | 性少妇av在线| www.自偷自拍.com| 大香蕉久久成人网| 在线观看免费视频日本深夜| 亚洲人成电影观看| 999久久久精品免费观看国产| 少妇裸体淫交视频免费看高清 | 亚洲九九香蕉| 精品高清国产在线一区| 一个人免费在线观看的高清视频| 咕卡用的链子| 狠狠狠狠99中文字幕| 热99国产精品久久久久久7| 视频区欧美日本亚洲| 啦啦啦视频在线资源免费观看| 亚洲av电影在线进入| 精品少妇内射三级| 国产男女超爽视频在线观看| 国产高清国产精品国产三级| 免费在线观看完整版高清| 黄色怎么调成土黄色| 国产在线视频一区二区| 怎么达到女性高潮| 一二三四社区在线视频社区8| 免费在线观看日本一区| 国产免费现黄频在线看| 无限看片的www在线观看| 男女下面插进去视频免费观看| 曰老女人黄片| av网站免费在线观看视频| 亚洲成人免费电影在线观看| 久9热在线精品视频| 另类精品久久| 一级片免费观看大全| 国产精品亚洲一级av第二区| 国产色视频综合| 嫩草影视91久久| 国产区一区二久久| 亚洲国产毛片av蜜桃av| 亚洲国产欧美在线一区| 巨乳人妻的诱惑在线观看| 美女扒开内裤让男人捅视频| 国产深夜福利视频在线观看| 老司机午夜十八禁免费视频| 波多野结衣av一区二区av| 老司机福利观看| a级片在线免费高清观看视频| 久久影院123| 久久久水蜜桃国产精品网| 男男h啪啪无遮挡| 在线看a的网站| 亚洲专区国产一区二区| 青草久久国产| 亚洲中文字幕日韩| 91国产中文字幕| 日本黄色视频三级网站网址 | 黄频高清免费视频| 手机成人av网站| 最新的欧美精品一区二区| 欧美日韩黄片免| 亚洲国产av新网站| 巨乳人妻的诱惑在线观看|