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

    基于溫度的ET0模型改進及在海河流域的適用性

    2021-03-28 02:58:18張戴煒王春穎韓宇平
    人民黃河 2021年3期
    關(guān)鍵詞:模型

    張戴煒 王春穎 韓宇平

    摘 要:為提高HS模型估算海河流域參考作物蒸散量(ET0)精度,基于貝葉斯原理,利用海河流域8個地理分區(qū)1980—1999年29個氣象站的逐日氣象數(shù)據(jù)對其溫度指數(shù)、溫度系數(shù)和溫度常數(shù)進行改進,并以FAO Penman-Monteith(PM)模型為標準,利用1980—1999年和2000—2015年數(shù)據(jù)對HS改進模型適用性進行評價。海河流域HS改進模型的溫度系數(shù)、溫度常數(shù)和溫度指數(shù)的范圍分別為5×10-4~8×10-4、28.5~38.0和0.4~0.7。與PM 模型計算結(jié)果對比,HS改進模型均方根誤差(RMSE)下降,相對誤差(RE)明顯降低,1980—1999年RMSE由0.92下降到0.84,RE從-13%~3%降至-4%~1%,2000—2015年RMSE由1.00下降到0.96,RE從-18%~2%降至-16%~-1%。HS改進模型計算結(jié)果與PM計算結(jié)果的時空變化一致,具有較高的計算精度和較好的適應(yīng)性。HS改進模型有助于提高海河流域缺少風(fēng)速、輻射等氣象資料條件下ET0估算精度,可為海河流域ET0簡化計算、作物需水量計算及灌溉制度的確立提供參考。

    關(guān)鍵詞:參考作物蒸散量;Hargreaves-Samani 模型;貝葉斯理論;適用性分析;海河流域

    Abstract:In order to improve the accuracy of reference crop evapotranspiration (ET0) calculation with Hargreaves-Samani (HS) model, based on the Bayesian principle and the daily weather data of 29 meteorological stations in the eight geographical divisions of Haihe River basin from 1980 to 1999, this paper improved the temperature exponent (C), temperature coefficient (a) and temperature constant (m) of HS model. Furthermore, the FAO Penman-Monteith (PM) Model was regarded as the evaluation standard to estimate accuracy and adaptability of the improved HS model, data from 2000 to 2015 were used for model adaptability verification. The results show that the range of C, a and m coefficients are 5×10-4-8×10-4, 28.5-38 and 0.4-0.7 respectively. The absolute deviation, relative deviation and relative error of the improved HS model ET0 were significantly reduced, the relative error (RE) and the root mean square error (RMSE) of between the improved HS model ET0 and the PM model ET0 were reduced from -13%-3% to -4%-1%, 0.92 to 0.84 during 1980-1999, from -18%-2% to -16%-1%, 1 to 0.96 during 2000-2015. The spatial and temporal variation analysis shows that the improved HS model ET0 is basically consistent with the PM model ET0 and shows less difference. The study result will help to improve the estimation accuracy of ET0 in different climatic regions under limited meteorological data (without wind speed and radiation). The improved HS model can be used as a simplified recommended model for calculation of reference crop evapotranspiration in Haihe River basin, which provides a simple method for the determination of crop water demand and establishment of irrigation schedule in the area.

    Key words: evapotranspiration; Hargreaves-Samani (HS)Model; Bayesian Principles; adaptability analysis; Haihe River basin

    蒸散發(fā)是水文過程的重要環(huán)節(jié),參考作物蒸散量(reference crop evapotranspiration, ET0)是衡量蒸散發(fā)的綜合表現(xiàn)的通量。國際上普遍采用聯(lián)合國糧農(nóng)組織(Food and Agriculture Organization, FAO)于1998年提出的Penman-Monteith(PM)模型[1]計算ET0。PM模型需要大量氣象數(shù)據(jù),包括平均氣溫、最高氣溫、最低氣溫、風(fēng)速、輻射等。然而,缺乏氣象資料的地區(qū)很難獲取完整的氣象數(shù)據(jù),因此PM模型的應(yīng)用受到氣象數(shù)據(jù)完整性的限制。如何在缺乏大量氣象數(shù)據(jù)資料的情況下計算ET0,成為水文學(xué)、生態(tài)學(xué)及農(nóng)業(yè)水利領(lǐng)域的熱點課題[1]。眾多估算ET0的簡化公式不斷涌現(xiàn),如基于溫度的Blaney and Criddle、McCloud、Hamon、Thornthwaite、Hargreaves-Samani等公式,基于輻射的Jensen and Haise、Makkink、Priestley-Taylor、Irmak等公式[2]。基于輻射的方法所需要的參數(shù)較多,且輻射類的有關(guān)數(shù)據(jù)難以獲取,相比之下溫度類方法的應(yīng)用更為廣泛。其中,Hargreaves-Samani(HS)模型對于原始氣象數(shù)據(jù)資料的要求較低,在國內(nèi)外得到了廣泛應(yīng)用[3]。但是,以往研究結(jié)果表明HS模型計算結(jié)果存在誤差過大,甚至在部分地區(qū)出現(xiàn)蒸散量為負數(shù)的情況,因此HS模型需要進行修正以提高模型精度[4]。

    國內(nèi)外學(xué)者圍繞HS模型的參數(shù)校正、誤差影響因素分析、適用性評價等問題等開展了大量研究。楊永紅等[5]為了進一步提高HS模型在拉薩地區(qū)的適用性,引入了平均相對濕度因子,建立改進HS模型并進行了評價 。李晨等[6]基于貝葉斯方法改進的HS模型在川中丘陵區(qū)不同區(qū)域變異性較小,適應(yīng)性較強,具有較高的計算精度 。胡慶芳等[7]對HS模型在中國7個區(qū)域進行全局校正及適用性評價。范文波等[8]采用貝葉斯方法對HS模型在內(nèi)陸河流域進行修正及驗證,結(jié)果顯示修正后的HS模型滿足精度要求,且計算簡便,可為內(nèi)陸河流域ET0的計算提供參考。Xu等[9]在瑞士Vaud地區(qū)比較了HS模型和其他5種經(jīng)驗公式的計算精度,發(fā)現(xiàn)HS模型校準后能準確計算年ET0,但計算精度在所有方法中位居中等。Temesgen等[10]指出在濕潤低風(fēng)的環(huán)境下HS模型的計算結(jié)果優(yōu)于干燥強風(fēng)的環(huán)境,風(fēng)速是影響HS模型適用性的重要因素。此外,有研究指出在保留HS模型簡潔性的基礎(chǔ)上,需充分考慮區(qū)域溫度、降雨、輻射等氣象因子對HS模型參數(shù)和精度的影響[11-12]。

    HS模型簡捷、易算、對數(shù)據(jù)資料要求較低,但未考慮風(fēng)速、相對濕度等其他氣象因素影響,并且未充分考慮參數(shù)在不同區(qū)域的差異性。為提高HS模型估算海河流域參考作物蒸散量的計算精度和體現(xiàn)海河流域8個地理分區(qū)HS模型參數(shù)的空間差異性,本文基于1980—1999年海河流域29個氣象站的逐日氣象數(shù)據(jù)資料,引入貝葉斯理論對HS模型的溫度指數(shù)、溫度系數(shù)和溫度常數(shù)進行改進,以PM模型計算的蒸散量作為評價標準進行對比分析,并利用2000—2015年數(shù)據(jù)資料對改進的HS模型進行適用性評價。

    1 材料與方法

    1.1 研究區(qū)概況及數(shù)據(jù)來源

    海河流域位于北緯35°—43°、東經(jīng)112°—120°,包括北京市、天津市、河北省大部分地區(qū),山西省東部,河南省北部及山東、遼寧省、內(nèi)蒙古自治區(qū)的一部分地區(qū),總面積約31.8萬km2,約占全國面積的3%。海河流域是我國重要的糧食生產(chǎn)基地,其農(nóng)作物以冬小麥和夏玉米兩熟輪作制為主,其中夏玉米生育期為6—9月,冬小麥生育期為10月—次年5月。流域多年平均降水量540 mm,水面蒸發(fā)量1 100 mm,年平均氣溫1.5~14.0 ℃,多年平均相對濕度50%~70%,屬半濕潤半干旱地帶。

    本文根據(jù)水系分布和地貌特征將海河流域劃分為8個分區(qū)(見圖1)。海河流域有海河、灤河和徒駭馬頰河3大水系,其中:海河水系是主要水系,由北部的北運河、永定河和南部的大清河、子牙河、黑龍港運東、漳衛(wèi)河等組成;灤河水系包括灤河及冀東沿海諸河;徒駭馬頰河水系位于流域最南部。逐日氣象數(shù)據(jù)來源于國家氣象信息中心29個氣象站。

    1.2 蒸散發(fā)模型及貝葉斯理論

    1.2.1 PM模型

    聯(lián)合國糧農(nóng)組織推薦PM公式計算參考作物蒸散量,PM公式形式如下:

    1.2.2 HS模型

    1.2.3 基于貝葉斯理論的HS改進模型

    為提高HS模型估算海河流域ET0的精度,本文引入貝葉斯理論對HS模型3個參數(shù)進行修正與改進。貝葉斯原理是根據(jù)現(xiàn)有資料來估計該事物后期發(fā)生的可能性,并得到可能結(jié)果的一種統(tǒng)計推斷,這種理論方法包括先驗分布和后驗分布。當獲得新的樣本資料時,可以將樣本資料獲得的后驗分布作為新一輪評估的先驗信息,獲得新的后驗分布,從而獲得更加準確的估計參數(shù),提高模型的估計準確度[15]。本研究通過求解C、a、m的數(shù)學(xué)期望值μC、μa、μm改進HS模型。利用海河流域8個地理分區(qū)1980—1999年29個氣象站PM公式計算的ET0以及平均溫度、溫差、輻射獲得μm、μC、μa的后驗分布,取平均值作為HS改進模型的參數(shù)值。

    本文依據(jù)1980—1999年和2000—2015年海河流域29個氣象站的氣象數(shù)據(jù)資料進行PM和HS模型ET0計算,把PM模型作為參照標準,利用1980—1999年數(shù)據(jù)改進HS模型并用2000—2015年的數(shù)據(jù)驗證HS改進模型的適用性,分析HS模型改進前后計算ET0的精度。本研究選取均方根誤差RMSE和相對誤差RE作為模型精度的評價指標。

    2 結(jié)果與分析

    2.1 改進的HS模型參數(shù)分析

    海河流域不同分區(qū)HS改進模型的溫度系數(shù)、溫度常數(shù)和溫度指數(shù)最優(yōu)值見表1。海河流域HS改進模型的參數(shù)空間插值分布見圖2。由表1和圖2可見,3個參數(shù)在空間分布上呈較為明顯的地帶性特征。在海河流域邊緣地帶,如灤河、永定河、黑龍港及運東、漳衛(wèi)河等參數(shù)空間分布呈現(xiàn)兩極化。溫度系數(shù)C西南低東北高,最高值位于天津,最低值在五臺山、南宮,沿漳衛(wèi)河至灤河區(qū)域隨緯度增大而呈現(xiàn)不斷增大的特征。永定河和灤河的溫度系數(shù)C值分別為7.5×10-4、7.4×10-4,均高于溫度系數(shù)C在海河流域的平均值6.17×10-4。這二者均位于流域的北側(cè),可能與其地理位置、緯度、太陽輻射和溫度極值相關(guān)。溫度常數(shù)a在海河流域的擬合平均值為33.25,漳衛(wèi)河和北三河的數(shù)值呈現(xiàn)兩極化,擬合值分別為28.61和37.49,空間分布呈中部高南北低,高值位于五臺山和天津,低值位于承德、青龍、邢臺、蔚縣和新鄉(xiāng)。五臺山位于內(nèi)陸中部,天津位于入海口,而邢臺、新鄉(xiāng)等位于平原區(qū)域,這說明地形及氣候?qū)囟瘸?shù)a的影響相對較大。溫度指數(shù)m與溫度系數(shù)C的空間分布相反,呈現(xiàn)西南高東北低的態(tài)勢,溫度指數(shù)隨高程提高而減小。溫度指數(shù)m在海河流域的平均值為0.56,整體空間差異小。

    2.2 時空變化趨勢分析

    改進前后HS模型計算出的ET0與PM模型計算出的ET0年內(nèi)變化趨勢相同,且與氣溫的年度變化趨勢一致,均呈開口向下的拋物線狀,改進后的HS模型ET0計算值與PM模型的ET0更接近(見圖3)。但是,1985年北三河、灤河分區(qū)HS模型改進前ET0相比改進后更為接近PM標準值,1990年北三河、大清河分區(qū)HS模型改進前ET0相比改進后更接近PM標準值;全流域在1985年與1990年,改進前HS模型相比模型改進后效果更好。有些地區(qū)HS模型ET0計算值與PM模型ET0計算值差別較大,例如黑龍港及運東(見圖3)。其原因可能是冬季有些氣象站平均溫度低,導(dǎo)致HS模型ET0計算值為負值;也可能是夏秋季節(jié)雨水充足,空氣濕度較大,海拔較高地區(qū)和沿海地區(qū)風(fēng)速較大,而HS模型未考慮濕度項和風(fēng)速項對ET0的影響。

    PM模型和改進前后的HS模型計算出來的ET0進行空間對比,3個模型計算ET0的空間插值分布見圖4。由圖4可以看出,PM模型ET0與HS改進前ET0、改進后HS模型ET0的空間分布規(guī)律基本一致。PM模型ET0空間插值分布結(jié)果表明,北三河分區(qū)的天津氣象站附近蒸散量最高(1 360 mm),子牙河分區(qū)五臺山氣象站附近最低(880 mm)(見圖4(a))。HS模型改進后的ET0空間插值分布結(jié)果表明北三河分區(qū)的天津氣象站附近蒸散量最高(1 360 mm),子牙河分區(qū)五臺山氣象站附近最低(800 mm)(見圖4(c))。然而,HS模型改進前的ET0空間插值分布結(jié)果表明,黑龍港運東分區(qū)的南宮氣象站附近蒸散量最高(1 120 mm),子牙河分區(qū)的五臺山氣象站附近最低(560 mm)(見圖4(b))。因此,改進后的HS模型計算得出的ET0更加接近PM模型標準值,改進后的HS模型比改進前的HS模型從時間和空間分布上都有明顯提升,表明模型改進效果較為顯著。

    2.3 適用性分析

    1980—1999年海河流域 29個氣象站HS模型改進前計算ET0值與PM模型ET0值的差值(ET0(HS)-ET0(PM))見圖5(a)。結(jié)果表明,兩者平均偏差值為55.3 mm。該時間段29個氣象站HS模型改進后計算ET0值與PM模型差值(ET0(改進后HS)-ET0(PM))見圖5(b),兩者平均偏差值為11.9 mm。改進后HS模型相比改進前,同PM模型計算的ET0相比,空間差值顯著減小。

    選取海河流域8個分區(qū)1980—1999年和2000—2015年兩個時段,將PM模型計算的ET0與改進前和改進后HS模型計算的ET0對比,分別見表2和表3。由表2和表3可知,1980—1999年HS模型改進后均方根誤差RMSE由0.92降低至0.84,2000—2015年HS模型改進后RMSE值由1.00降低至0.96,RMSE整體呈現(xiàn)降低趨勢。其中大清河分區(qū)在不同時段RMSE改進后相比改進前均降低0.3左右,精度提高明顯。然而,1980—1999年漳衛(wèi)河分區(qū)和2000—2015年徒駭馬頰河分區(qū),HS模型改進后RMSE分別上升0.15、0.38,精度下降。1980—1999年改進后蒸散量的相對誤差RE(-4%~1%)明顯低于HS模型改進前相對誤差(-13%~3%);2000—2015年改進后蒸散量的相對誤差(-16%~-1%)明顯低于HS模型改進前相對誤差(-18%~2%),表明模型改進后精度有顯著提升。整體上,改進后的HS模型相比改進前計算結(jié)果更為精確與穩(wěn)定,改進后的HS模型在海河流域具有較好的適用性。

    3 討 論

    HS模型是基于溫度估算氣象數(shù)據(jù)有限地區(qū)的參考作物蒸散量的可靠方法,降低了ET0計算的氣象數(shù)據(jù)要求,簡化了計算過程。但HS模型經(jīng)驗性較強, Hargreaves 等[16] 推薦當Ra以MJ/(m2·d)為單位時,溫度系數(shù)C取值范圍為8.16×10-4~12.24×10-4。海河流域溫度系數(shù)C取值范圍為5×10-4~8×10-4,低于Hargreaves推薦值。李晨等[6]研究結(jié)果與本研究結(jié)果對比表明,川中丘陵區(qū)溫度系數(shù)、溫度常數(shù)和溫度指數(shù)的空間差異性較小,但海河流域HS改進模型3個參數(shù)的范圍分別為5×10-4~8×10-4、28.5~38.0和0.4~0.7,空間差異性更大??紤]海河流域HS模型參數(shù)的空間差異性,有利于提高HS模型計算ET0的精度。

    此外,為減小HS模型因忽略風(fēng)速、濕度、降雨等氣象因素影響對ET0計算所帶來的誤差,國內(nèi)外學(xué)者對HS 模型進行了改進,改進后HS模型形式更為復(fù)雜。Allen[17]通過研究認為在高程超過1 500 m的地區(qū)Hargreaves推薦值并不合理;Annandale等[18] 認為HS模型參數(shù)取值應(yīng)考慮大氣壓因素;Samani等[14] 認為在研究區(qū)域為北緯7°—50°時,修正參數(shù)應(yīng)依據(jù)研究區(qū)域的溫度梯度確定。本研究改進HS模型是在只考慮溫度的大背景下進行的,模型誤差存在空間差異性,對于HS改進模型誤差較大的地點,應(yīng)進一步考慮高程、風(fēng)速、濕度等因素改進HS模型。

    4 結(jié) 語

    為提高Hargreaves-Samani(HS)模型在海河流域的適用性,通過貝葉斯原理對HS模型的溫度系數(shù)、溫度常數(shù)和溫度指數(shù)進行改進,海河流域3個參數(shù)范圍分別為5×10-4~8×10-4、28.5~38.0和0.4~0.7。與PM 模型計算結(jié)果對比,HS改進模型的偏差減小,均方根誤差RMSE下降,相對誤差RE明顯降低。改進后的HS模型計算得到的ET0精度顯著提高,可作為海河流域ET0簡化計算的方法。本研究涉及的區(qū)域范圍相對較小,在海拔較高和沿海地區(qū)HS模型的適用性相對較差,這說明高程、風(fēng)速及濕度等因素也對ET0存在影響,有待于進一步研究。

    參考文獻:

    [1] ALLEN R G, PEREIRA L S, RAES D, et al. Crop Evapotranspiration - Guidelines for Computing Crop Water Requirements[M]. Rome: Food and Agriculture Organization of the United Nation, 1998:11-49.

    [2] 樊軍, 邵明安, 王全九.黃土區(qū)參考作物蒸散量多種計算方法的比較研究[J]. 農(nóng)業(yè)工程學(xué)報, 2008, 24(3):98-102.

    [3] HARGREAVES G H, ALLEN R G. History and Evaluation of Hargreaves Evapotranspiration Equation[J]. Journal of Irrigation and Drainage Engineering ASCE,2003, 192: 53-63.

    [4] 張本興,潘云,李小娟.中國不同氣候區(qū)域Hargreaves模型的修正[J].地理與地理信息科學(xué),2012,28(1):51-54.

    [5] 楊永紅,張展羽.改進Hargreaves方法計算拉薩參考作物蒸發(fā)蒸騰量[J].水科學(xué)進展,2009,20(5):614-618.

    [6] 李晨,崔寧博,魏新平,等.改進Hargreaves模型估算川中丘陵區(qū)參考作物蒸散量[J].農(nóng)業(yè)工程學(xué)報,2015,31(11):129-135.

    [7] 胡慶芳,楊大文,王銀堂,等.Hargreaves公式的全局校正及適用性評價[J].水科學(xué)進展,2011,22(2):160-167.

    [8] 范文波,吳普特,韓志全,等.瑪納斯河流域ET0影響因子分析及對Hargreaves法的修正[J].農(nóng)業(yè)工程學(xué)報,2012,28(8):19-24.

    [9] XU C Y, SINGH V P. Evaluation and Generalization of Temperature Based Methods for Calculating Evaporation[J]. Hydrological Processes, 2001, 15(2): 305-319.

    [10] TEMESGEN B, ECHINGS, DAVIDOFF B, et al. Comparison of Some Reference Evapotranspiration Equations for California[J]. Journal of Irrigation and Drainage Engineering, 2005, 131(1): 73-84.

    [11] GIUSEPPE M, SENATORE A. Regionalization of the Hargreaves Coefficient for the Assessment of Distributed Reference Evaportransportation in Southern Italy[J]. Journal of Irrigation and Drainage Engineering, 2013, 139(5): 349-362.

    [12] JABLOUN M, SAHLI A. Evaluation of FAO-56 Methodology for Estimating Reference Evaportransportation Using Limited Climatic Data Application to Tunisia[J]. Agric Water Manage, 2008, 95(6): 707-715.

    [13] HARGREAVES G H, ALLEN R G. History and Evaluation of Hargreaves Evapotranspiration Basin[J]. Irrig Drain Eng, 1985, 111(2):133-124

    [14] SAMANI Z A. Estimating Solar Radiation and Evapotransporation Using Minimum Climatological Data[J]. Irrig Drain Eng,2000,126(4):265-267

    [15] 任曉明,李章呂.貝葉斯決策理論的發(fā)展概況和研究動態(tài)[J].科學(xué)技術(shù)哲學(xué)研究,2013,30(2):1-7.

    [16] HARGREAVES G L, HARGREAVES G H, RILEY J P. Agricultural Benefits for Senegal River Basin[J]. Irrig Drain Eng,1985, 111: 113-124.

    [17] ALLEN R G. Self-Calibrating Method for Estimating Solarradiation from Air Temperature[J]. Hydrol Eng, 1997, 2(2): 56-67.

    [18] ANNANDALE J G, JOVANIC N Z, BENADE N, et al. Softwarefor Missing Data Error Analysis of Penman-Monteithreference Evaportransportation[J]. Irrig Sci, 2002, 21(2): 57-67.

    【責(zé)任編輯 許立新】

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    999精品在线视频| 看免费成人av毛片| 在线亚洲精品国产二区图片欧美| 黄色视频在线播放观看不卡| 久久国产亚洲av麻豆专区| 国产人伦9x9x在线观看 | 最近手机中文字幕大全| 丝袜美足系列| 国产精品 欧美亚洲| 亚洲国产av影院在线观看| 久久精品久久久久久噜噜老黄| 国产精品成人在线| 2022亚洲国产成人精品| 狠狠婷婷综合久久久久久88av| 精品人妻熟女毛片av久久网站| 欧美精品高潮呻吟av久久| 久久久a久久爽久久v久久| 自线自在国产av| 欧美亚洲日本最大视频资源| 人妻少妇偷人精品九色| 欧美国产精品va在线观看不卡| 久久精品久久精品一区二区三区| 欧美精品av麻豆av| 国产视频首页在线观看| 亚洲精品国产色婷婷电影| 男女国产视频网站| 日本猛色少妇xxxxx猛交久久| 美女国产视频在线观看| 成人国语在线视频| 亚洲美女视频黄频| 一区在线观看完整版| 视频区图区小说| 精品一区二区三区四区五区乱码 | 王馨瑶露胸无遮挡在线观看| 最近中文字幕高清免费大全6| xxxhd国产人妻xxx| 观看美女的网站| 人人澡人人妻人| 成人影院久久| 午夜免费鲁丝| 日韩成人av中文字幕在线观看| 国产一区二区 视频在线| 我的亚洲天堂| 国产视频首页在线观看| 成人影院久久| 综合色丁香网| 精品第一国产精品| 大片电影免费在线观看免费| 欧美国产精品一级二级三级| 一级毛片电影观看| 成人漫画全彩无遮挡| 久久国产亚洲av麻豆专区| 日日摸夜夜添夜夜爱| 丝袜美足系列| 午夜91福利影院| 美女国产高潮福利片在线看| www日本在线高清视频| 天堂俺去俺来也www色官网| 晚上一个人看的免费电影| 如何舔出高潮| 精品国产超薄肉色丝袜足j| 久久久精品区二区三区| 在线观看www视频免费| 女的被弄到高潮叫床怎么办| 亚洲精品一区蜜桃| 国产精品秋霞免费鲁丝片| 最新的欧美精品一区二区| 91午夜精品亚洲一区二区三区| 欧美亚洲 丝袜 人妻 在线| 欧美人与性动交α欧美精品济南到 | 在线 av 中文字幕| 久久久久网色| 一级黄片播放器| 久久久久久久大尺度免费视频| 欧美 日韩 精品 国产| 老熟女久久久| 免费黄频网站在线观看国产| 国产日韩欧美在线精品| 国产免费又黄又爽又色| 人妻 亚洲 视频| 亚洲精品一二三| 看免费成人av毛片| 亚洲一码二码三码区别大吗| 欧美bdsm另类| 自拍欧美九色日韩亚洲蝌蚪91| 啦啦啦中文免费视频观看日本| 欧美人与性动交α欧美软件| 18禁国产床啪视频网站| 狂野欧美激情性bbbbbb| 美女视频免费永久观看网站| 色吧在线观看| 日本vs欧美在线观看视频| 亚洲精品国产一区二区精华液| 中国国产av一级| 国产精品一区二区在线不卡| 久久人人爽av亚洲精品天堂| 丝袜脚勾引网站| 日韩制服丝袜自拍偷拍| 男女午夜视频在线观看| 日韩欧美一区视频在线观看| a级片在线免费高清观看视频| 国产一区二区激情短视频 | 久久国产精品大桥未久av| 天天操日日干夜夜撸| 日本av免费视频播放| 大码成人一级视频| 中国国产av一级| 热re99久久国产66热| 丝袜脚勾引网站| 久久久久久人人人人人| 亚洲国产毛片av蜜桃av| tube8黄色片| 另类精品久久| 国精品久久久久久国模美| 欧美中文综合在线视频| 亚洲经典国产精华液单| 十分钟在线观看高清视频www| 日韩一区二区三区影片| 最黄视频免费看| 亚洲国产精品999| 一区二区三区乱码不卡18| 国产一区二区三区av在线| 国产成人一区二区在线| 人人妻人人爽人人添夜夜欢视频| 啦啦啦中文免费视频观看日本| 久久免费观看电影| 老熟女久久久| 日日爽夜夜爽网站| 亚洲伊人久久精品综合| 在线精品无人区一区二区三| 午夜激情av网站| 七月丁香在线播放| 亚洲成色77777| 欧美日韩亚洲高清精品| 大片电影免费在线观看免费| 国产成人精品无人区| 国产亚洲av片在线观看秒播厂| 婷婷色综合www| 男人操女人黄网站| 免费在线观看完整版高清| 国产 精品1| 国产精品 欧美亚洲| 寂寞人妻少妇视频99o| 丝袜在线中文字幕| 国产精品欧美亚洲77777| 肉色欧美久久久久久久蜜桃| 国产一级毛片在线| 最近中文字幕高清免费大全6| 婷婷色麻豆天堂久久| 一个人免费看片子| 欧美人与性动交α欧美软件| 欧美日韩视频精品一区| 久久毛片免费看一区二区三区| 欧美日韩精品网址| 又粗又硬又长又爽又黄的视频| 久久久亚洲精品成人影院| 欧美97在线视频| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区 视频在线| 午夜日本视频在线| 美女午夜性视频免费| 欧美日韩av久久| 日韩三级伦理在线观看| 国产av精品麻豆| videos熟女内射| 久久国内精品自在自线图片| av又黄又爽大尺度在线免费看| 欧美激情高清一区二区三区 | 免费播放大片免费观看视频在线观看| 国产精品一二三区在线看| 婷婷色麻豆天堂久久| 一区福利在线观看| 如何舔出高潮| 热99国产精品久久久久久7| 日韩av免费高清视频| 国产精品成人在线| 大话2 男鬼变身卡| 国产探花极品一区二区| 国产免费又黄又爽又色| 人人妻人人爽人人添夜夜欢视频| 国产成人午夜福利电影在线观看| 国产在线视频一区二区| 看十八女毛片水多多多| 日韩一本色道免费dvd| 99香蕉大伊视频| 美女主播在线视频| 国产黄频视频在线观看| a级毛片在线看网站| 国产女主播在线喷水免费视频网站| 黑丝袜美女国产一区| 精品酒店卫生间| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久久国产精品人妻一区二区| 久久久久精品性色| 高清欧美精品videossex| 国产在视频线精品| 午夜福利影视在线免费观看| 精品少妇内射三级| www.精华液| 国产男女内射视频| 精品久久蜜臀av无| 18禁裸乳无遮挡动漫免费视频| 寂寞人妻少妇视频99o| 亚洲欧美一区二区三区黑人 | 9色porny在线观看| 国产男女内射视频| 在线观看国产h片| av一本久久久久| 国产男女内射视频| 国产激情久久老熟女| 国产男人的电影天堂91| av免费在线看不卡| 2018国产大陆天天弄谢| 中文精品一卡2卡3卡4更新| xxx大片免费视频| 日韩av免费高清视频| 午夜久久久在线观看| 亚洲精品乱久久久久久| 亚洲av男天堂| av福利片在线| 亚洲av成人精品一二三区| 人人妻人人爽人人添夜夜欢视频| 日韩精品免费视频一区二区三区| 亚洲三级黄色毛片| av国产久精品久网站免费入址| 欧美中文综合在线视频| 亚洲天堂av无毛| 如日韩欧美国产精品一区二区三区| 大香蕉久久成人网| 国产成人精品久久久久久| 永久网站在线| 男的添女的下面高潮视频| 免费人妻精品一区二区三区视频| 成年女人毛片免费观看观看9 | 欧美精品国产亚洲| 人人澡人人妻人| 欧美 亚洲 国产 日韩一| 18禁动态无遮挡网站| 男女无遮挡免费网站观看| 熟女电影av网| 两性夫妻黄色片| 欧美国产精品一级二级三级| 国产精品久久久久成人av| 精品亚洲乱码少妇综合久久| 免费在线观看完整版高清| 成人毛片a级毛片在线播放| 26uuu在线亚洲综合色| 国产一区二区 视频在线| 成年人午夜在线观看视频| 伊人久久国产一区二区| 不卡视频在线观看欧美| 亚洲av中文av极速乱| 国产成人精品婷婷| videossex国产| 三上悠亚av全集在线观看| 久久av网站| 婷婷色麻豆天堂久久| 高清av免费在线| 99久国产av精品国产电影| 99久久人妻综合| 国产av精品麻豆| 国产成人91sexporn| 少妇猛男粗大的猛烈进出视频| a 毛片基地| 在线观看免费视频网站a站| 亚洲第一区二区三区不卡| 久久这里只有精品19| 伦理电影免费视频| 水蜜桃什么品种好| 日韩在线高清观看一区二区三区| 人成视频在线观看免费观看| 大片电影免费在线观看免费| 99久国产av精品国产电影| av在线app专区| 国产有黄有色有爽视频| 亚洲欧洲国产日韩| 啦啦啦在线免费观看视频4| 中文字幕亚洲精品专区| 国产精品国产av在线观看| 大陆偷拍与自拍| 少妇人妻 视频| 国精品久久久久久国模美| 性色av一级| kizo精华| 亚洲av欧美aⅴ国产| 亚洲av福利一区| 免费在线观看视频国产中文字幕亚洲 | 国产男女超爽视频在线观看| 两性夫妻黄色片| 校园人妻丝袜中文字幕| 成人二区视频| 亚洲国产精品999| 亚洲激情五月婷婷啪啪| 又粗又硬又长又爽又黄的视频| 欧美人与性动交α欧美软件| 亚洲少妇的诱惑av| 欧美日韩综合久久久久久| 日本vs欧美在线观看视频| 9191精品国产免费久久| 精品99又大又爽又粗少妇毛片| 捣出白浆h1v1| 国产一区二区在线观看av| 丝袜美足系列| 亚洲少妇的诱惑av| 香蕉丝袜av| 国产精品人妻久久久影院| 亚洲成人手机| 亚洲五月色婷婷综合| 美女视频免费永久观看网站| 极品少妇高潮喷水抽搐| 黑丝袜美女国产一区| 黑丝袜美女国产一区| 中文天堂在线官网| 国产极品天堂在线| 免费高清在线观看视频在线观看| 日韩不卡一区二区三区视频在线| 中文字幕色久视频| 啦啦啦中文免费视频观看日本| 中文乱码字字幕精品一区二区三区| 黄网站色视频无遮挡免费观看| 欧美日韩亚洲高清精品| 国产免费又黄又爽又色| 女人久久www免费人成看片| 日韩一卡2卡3卡4卡2021年| 又黄又粗又硬又大视频| 天天影视国产精品| 麻豆av在线久日| 一区二区三区四区激情视频| 精品亚洲乱码少妇综合久久| 免费观看在线日韩| 欧美另类一区| 激情视频va一区二区三区| www.av在线官网国产| 久久精品久久久久久噜噜老黄| 丰满少妇做爰视频| 国产又爽黄色视频| 在线看a的网站| 91久久精品国产一区二区三区| 999精品在线视频| 久久精品人人爽人人爽视色| 精品久久久久久电影网| 亚洲精品日本国产第一区| 蜜桃在线观看..| 欧美日韩视频精品一区| 少妇猛男粗大的猛烈进出视频| 免费黄网站久久成人精品| 伦理电影大哥的女人| 十分钟在线观看高清视频www| 女性被躁到高潮视频| 免费在线观看完整版高清| 日本wwww免费看| 夜夜骑夜夜射夜夜干| 美女午夜性视频免费| 黄片无遮挡物在线观看| 9191精品国产免费久久| 在线观看免费视频网站a站| 美女主播在线视频| 欧美人与性动交α欧美软件| 另类亚洲欧美激情| 波多野结衣一区麻豆| 亚洲国产欧美网| 国产精品香港三级国产av潘金莲 | 日韩中文字幕欧美一区二区 | 人妻人人澡人人爽人人| 国产淫语在线视频| 国产熟女午夜一区二区三区| 又黄又粗又硬又大视频| 国产精品国产av在线观看| 亚洲av成人精品一二三区| 欧美日韩精品网址| 久久久久久久大尺度免费视频| 大香蕉久久网| 欧美精品国产亚洲| av女优亚洲男人天堂| 久久精品国产亚洲av涩爱| 桃花免费在线播放| 秋霞在线观看毛片| 欧美日本中文国产一区发布| 一级毛片 在线播放| 一边摸一边做爽爽视频免费| 欧美日本中文国产一区发布| 少妇的丰满在线观看| 久久久久久免费高清国产稀缺| 亚洲av电影在线进入| 男女高潮啪啪啪动态图| 亚洲成人一二三区av| 亚洲av日韩在线播放| 黄色怎么调成土黄色| 久久久a久久爽久久v久久| 亚洲精品自拍成人| 久久韩国三级中文字幕| 精品人妻在线不人妻| 亚洲av成人精品一二三区| 水蜜桃什么品种好| 久久久久国产网址| 亚洲内射少妇av| 国产精品蜜桃在线观看| 天天操日日干夜夜撸| 亚洲色图综合在线观看| 国产精品香港三级国产av潘金莲 | 久久久久久久大尺度免费视频| 黄频高清免费视频| 黑人欧美特级aaaaaa片| 自线自在国产av| 一级,二级,三级黄色视频| 色视频在线一区二区三区| 涩涩av久久男人的天堂| 丝袜人妻中文字幕| 王馨瑶露胸无遮挡在线观看| 免费看不卡的av| 精品久久久久久电影网| 啦啦啦啦在线视频资源| 色94色欧美一区二区| 国产激情久久老熟女| 欧美日韩亚洲高清精品| 欧美精品一区二区大全| 色播在线永久视频| 午夜免费男女啪啪视频观看| 国产成人免费观看mmmm| 国产精品久久久久久精品古装| 国产女主播在线喷水免费视频网站| av在线播放精品| 久久婷婷青草| 色婷婷av一区二区三区视频| 三级国产精品片| 男女免费视频国产| 亚洲精品久久午夜乱码| 精品亚洲成国产av| 香蕉国产在线看| 高清欧美精品videossex| 欧美+日韩+精品| 一边亲一边摸免费视频| 男人舔女人的私密视频| 亚洲av中文av极速乱| 精品卡一卡二卡四卡免费| 欧美黄色片欧美黄色片| 一本色道久久久久久精品综合| 另类精品久久| 免费观看av网站的网址| 亚洲av电影在线观看一区二区三区| 老熟女久久久| 91aial.com中文字幕在线观看| 日本欧美国产在线视频| 两个人免费观看高清视频| 亚洲综合色网址| 极品少妇高潮喷水抽搐| 国产 一区精品| 下体分泌物呈黄色| 亚洲成人av在线免费| 我的亚洲天堂| 在线亚洲精品国产二区图片欧美| 啦啦啦视频在线资源免费观看| 十八禁高潮呻吟视频| 最黄视频免费看| 三级国产精品片| 亚洲av欧美aⅴ国产| 精品久久久精品久久久| 女人久久www免费人成看片| 亚洲情色 制服丝袜| 美女中出高潮动态图| 色婷婷av一区二区三区视频| 丝袜在线中文字幕| 秋霞在线观看毛片| 日本欧美视频一区| 亚洲国产av新网站| 性色av一级| a级毛片在线看网站| 丝袜美足系列| 十八禁网站网址无遮挡| 精品一区在线观看国产| 亚洲第一av免费看| av免费在线看不卡| 久久人人97超碰香蕉20202| 只有这里有精品99| 亚洲欧美中文字幕日韩二区| 久久久久久久久久久免费av| 亚洲av电影在线观看一区二区三区| 国产xxxxx性猛交| 久久久久视频综合| 久久精品国产鲁丝片午夜精品| 如日韩欧美国产精品一区二区三区| 亚洲精品成人av观看孕妇| 天天操日日干夜夜撸| 热99国产精品久久久久久7| 最近中文字幕2019免费版| av在线播放精品| 亚洲,欧美,日韩| 精品一区在线观看国产| 精品视频人人做人人爽| 久久久久国产一级毛片高清牌| 亚洲av在线观看美女高潮| 中文字幕人妻熟女乱码| 免费看av在线观看网站| 亚洲精品日韩在线中文字幕| av女优亚洲男人天堂| 欧美xxⅹ黑人| 精品国产国语对白av| 亚洲人成77777在线视频| 最近中文字幕高清免费大全6| 宅男免费午夜| 男女国产视频网站| 久久精品熟女亚洲av麻豆精品| 可以免费在线观看a视频的电影网站 | 国产人伦9x9x在线观看 | 国产精品免费大片| 亚洲欧美色中文字幕在线| 一级毛片 在线播放| av线在线观看网站| 欧美精品一区二区免费开放| 性色av一级| 欧美精品一区二区免费开放| 欧美精品亚洲一区二区| 电影成人av| 不卡av一区二区三区| 中文字幕人妻丝袜制服| 天天躁日日躁夜夜躁夜夜| 久久久国产精品麻豆| 少妇猛男粗大的猛烈进出视频| 老汉色av国产亚洲站长工具| 我的亚洲天堂| 亚洲成av片中文字幕在线观看 | 久久久精品区二区三区| 亚洲精品,欧美精品| 日韩一区二区三区影片| 亚洲精品视频女| 久久久久国产精品人妻一区二区| 亚洲av欧美aⅴ国产| 色婷婷久久久亚洲欧美| 在线看a的网站| 国产老妇伦熟女老妇高清| 满18在线观看网站| 国产精品熟女久久久久浪| 国产福利在线免费观看视频| 春色校园在线视频观看| 香蕉丝袜av| 超碰成人久久| 最近最新中文字幕大全免费视频 | 汤姆久久久久久久影院中文字幕| 人妻系列 视频| 久久精品亚洲av国产电影网| 亚洲三区欧美一区| 日本av免费视频播放| 黄色 视频免费看| 大香蕉久久成人网| 国产精品久久久久久精品古装| 观看av在线不卡| 国产精品久久久久久精品古装| 亚洲美女黄色视频免费看| 咕卡用的链子| 国产精品二区激情视频| 精品第一国产精品| 在线观看免费视频网站a站| 久久这里只有精品19| 老熟女久久久| 菩萨蛮人人尽说江南好唐韦庄| 涩涩av久久男人的天堂| 亚洲成国产人片在线观看| av在线老鸭窝| 国产一区亚洲一区在线观看| 国产日韩欧美在线精品| 2018国产大陆天天弄谢| 午夜福利乱码中文字幕| 亚洲欧洲精品一区二区精品久久久 | 日韩av免费高清视频| 久久久精品免费免费高清| 搡老乐熟女国产| 久久亚洲国产成人精品v| 久久ye,这里只有精品| 国产淫语在线视频| 中文字幕色久视频| 国产乱人偷精品视频| 国产精品一国产av| 欧美日韩精品网址| 国产熟女欧美一区二区| 激情五月婷婷亚洲| 新久久久久国产一级毛片| 在线观看免费日韩欧美大片| 咕卡用的链子| 一区二区日韩欧美中文字幕| 男人操女人黄网站| 亚洲av综合色区一区| 亚洲熟女精品中文字幕| 中文字幕制服av| 七月丁香在线播放| 亚洲国产看品久久| 国产成人精品婷婷| 亚洲欧美成人精品一区二区| 日本爱情动作片www.在线观看| 亚洲中文av在线| 亚洲国产最新在线播放| 午夜福利网站1000一区二区三区| 九色亚洲精品在线播放| 青草久久国产| av线在线观看网站| 国产男人的电影天堂91| 亚洲精华国产精华液的使用体验| 黄色怎么调成土黄色| 色94色欧美一区二区| 色婷婷久久久亚洲欧美| 18在线观看网站| 欧美日韩亚洲高清精品| 国产一区二区在线观看av| 欧美av亚洲av综合av国产av | 制服人妻中文乱码| 欧美激情极品国产一区二区三区| 一区二区三区乱码不卡18| 老鸭窝网址在线观看| √禁漫天堂资源中文www| 欧美黄色片欧美黄色片| 男男h啪啪无遮挡| 国产成人精品久久久久久| 久久久久国产精品人妻一区二区| 国产熟女午夜一区二区三区| 国产日韩欧美亚洲二区| 国产亚洲欧美精品永久|