• <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打印中的模型分割與打包
    久久精品久久久久久噜噜老黄 | av天堂在线播放| 999久久久精品免费观看国产| 在线播放无遮挡| 欧美3d第一页| 午夜精品久久久久久毛片777| 亚洲五月婷婷丁香| 欧美xxxx黑人xx丫x性爽| 午夜老司机福利剧场| av女优亚洲男人天堂| 狂野欧美白嫩少妇大欣赏| 成人鲁丝片一二三区免费| 又爽又黄无遮挡网站| 搡老妇女老女人老熟妇| 熟女电影av网| 色吧在线观看| 国产真实伦视频高清在线观看 | 我的老师免费观看完整版| 人人妻,人人澡人人爽秒播| 亚洲,欧美精品.| 69人妻影院| 性色avwww在线观看| 欧美成狂野欧美在线观看| 十八禁人妻一区二区| 亚洲人成电影免费在线| 欧美丝袜亚洲另类 | 国产精品永久免费网站| 伦理电影大哥的女人| 国产午夜精品久久久久久一区二区三区 | 久久天躁狠狠躁夜夜2o2o| 免费看日本二区| 在线国产一区二区在线| 日韩精品青青久久久久久| АⅤ资源中文在线天堂| 美女cb高潮喷水在线观看| 欧美中文日本在线观看视频| 成人高潮视频无遮挡免费网站| av天堂在线播放| 亚洲无线在线观看| xxxwww97欧美| 亚洲一区二区三区不卡视频| 日本与韩国留学比较| 欧美日韩中文字幕国产精品一区二区三区| 少妇人妻精品综合一区二区 | 久久这里只有精品中国| 色播亚洲综合网| 男人狂女人下面高潮的视频| 特级一级黄色大片| av国产免费在线观看| 午夜亚洲福利在线播放| 中出人妻视频一区二区| 亚洲欧美日韩东京热| 精品久久久久久,| 中文字幕熟女人妻在线| 欧美黄色淫秽网站| 国产精品不卡视频一区二区 | 色av中文字幕| 激情在线观看视频在线高清| 热99在线观看视频| 成人无遮挡网站| 色尼玛亚洲综合影院| 老熟妇乱子伦视频在线观看| 久久这里只有精品中国| 亚洲五月婷婷丁香| 国产av麻豆久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 欧美黄色片欧美黄色片| 国产人妻一区二区三区在| 亚洲欧美日韩高清在线视频| 亚洲一区二区三区色噜噜| 久久久久久久久久黄片| 真人一进一出gif抽搐免费| 精品久久久久久久久久免费视频| 每晚都被弄得嗷嗷叫到高潮| 九九热线精品视视频播放| 精品人妻熟女av久视频| 午夜免费男女啪啪视频观看 | 不卡一级毛片| 国产精品日韩av在线免费观看| 国产真实乱freesex| 丁香六月欧美| 最好的美女福利视频网| 天美传媒精品一区二区| 亚洲不卡免费看| 伦理电影大哥的女人| 亚洲美女视频黄频| 久久午夜亚洲精品久久| 成人欧美大片| 人人妻,人人澡人人爽秒播| 俄罗斯特黄特色一大片| 免费看日本二区| 欧美激情久久久久久爽电影| 色综合亚洲欧美另类图片| 久久午夜福利片| 级片在线观看| 亚洲国产高清在线一区二区三| 亚洲,欧美精品.| 久久久色成人| 天堂av国产一区二区熟女人妻| 久久热精品热| 综合色av麻豆| 欧美中文日本在线观看视频| 国产精品98久久久久久宅男小说| 亚洲av美国av| 免费观看的影片在线观看| 亚洲成a人片在线一区二区| 男人舔奶头视频| 国产亚洲精品av在线| 99久久99久久久精品蜜桃| 伦理电影大哥的女人| 国产爱豆传媒在线观看| 精品久久久久久久久亚洲 | 日韩欧美一区二区三区在线观看| 精品久久国产蜜桃| 午夜老司机福利剧场| 成人鲁丝片一二三区免费| 99热精品在线国产| 白带黄色成豆腐渣| 午夜福利高清视频| 欧美精品国产亚洲| 非洲黑人性xxxx精品又粗又长| 99久久99久久久精品蜜桃| 国产真实伦视频高清在线观看 | 亚洲精品乱码久久久v下载方式| 久久久久久久久久成人| 亚洲七黄色美女视频| 老熟妇乱子伦视频在线观看| 国产乱人伦免费视频| 亚洲av五月六月丁香网| 亚洲 国产 在线| 亚洲最大成人av| 美女高潮的动态| 激情在线观看视频在线高清| 99热只有精品国产| 欧美日韩福利视频一区二区| 伊人久久精品亚洲午夜| 国产高清激情床上av| 嫁个100分男人电影在线观看| 怎么达到女性高潮| 亚洲国产精品999在线| 成人午夜高清在线视频| 午夜精品一区二区三区免费看| 99久久久亚洲精品蜜臀av| a级毛片免费高清观看在线播放| 亚洲成a人片在线一区二区| 国产精品乱码一区二三区的特点| 他把我摸到了高潮在线观看| 很黄的视频免费| 欧美精品啪啪一区二区三区| 在线免费观看的www视频| 一级作爱视频免费观看| 内射极品少妇av片p| 久久久成人免费电影| 麻豆成人午夜福利视频| 欧美日韩瑟瑟在线播放| 国产国拍精品亚洲av在线观看| 97超级碰碰碰精品色视频在线观看| 国产毛片a区久久久久| 国产精品,欧美在线| 别揉我奶头~嗯~啊~动态视频| 午夜免费男女啪啪视频观看 | 国产毛片a区久久久久| 免费电影在线观看免费观看| 亚洲人成电影免费在线| 伦理电影大哥的女人| 欧美黑人欧美精品刺激| 国产精品伦人一区二区| 日韩中文字幕欧美一区二区| 露出奶头的视频| 又爽又黄无遮挡网站| 欧美日韩综合久久久久久 | 少妇裸体淫交视频免费看高清| 亚洲精品乱码久久久v下载方式| 日本黄大片高清| 国产一区二区三区在线臀色熟女| 在线播放国产精品三级| 亚洲真实伦在线观看| 男女视频在线观看网站免费| 9191精品国产免费久久| 久久精品国产清高在天天线| 成年免费大片在线观看| 亚洲精品在线美女| 亚洲狠狠婷婷综合久久图片| 亚洲天堂国产精品一区在线| 神马国产精品三级电影在线观看| 午夜福利在线观看吧| 美女被艹到高潮喷水动态| 亚洲国产精品999在线| 午夜福利视频1000在线观看| 757午夜福利合集在线观看| 国产高清视频在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 97碰自拍视频| 亚洲五月婷婷丁香| 国产熟女xx| 欧美午夜高清在线| 亚洲一区高清亚洲精品| 色哟哟哟哟哟哟| 在线观看美女被高潮喷水网站 | 国产不卡一卡二| 久久午夜亚洲精品久久| 亚洲经典国产精华液单 | 亚洲欧美日韩东京热| 欧美黄色片欧美黄色片| 国产成人啪精品午夜网站| 99久久久亚洲精品蜜臀av| 一区二区三区四区激情视频 | 欧美日韩乱码在线| 夜夜躁狠狠躁天天躁| 成人三级黄色视频| 国产人妻一区二区三区在| 婷婷亚洲欧美| 国产黄a三级三级三级人| 色综合欧美亚洲国产小说| 女同久久另类99精品国产91| 观看美女的网站| 久久久久久久亚洲中文字幕 | 国产在视频线在精品| 此物有八面人人有两片| 亚洲五月婷婷丁香| 网址你懂的国产日韩在线| 亚洲国产欧美人成| 又黄又爽又刺激的免费视频.| 精品人妻1区二区| 男插女下体视频免费在线播放| 亚洲国产日韩欧美精品在线观看| 亚洲精品色激情综合| 成人鲁丝片一二三区免费| 亚洲第一区二区三区不卡| 亚洲av第一区精品v没综合| 啦啦啦观看免费观看视频高清| 2021天堂中文幕一二区在线观| www.www免费av| 国产熟女xx| 国产高清视频在线观看网站| 亚洲精品456在线播放app | 精品人妻偷拍中文字幕| 在线观看免费视频日本深夜| 少妇被粗大猛烈的视频| 国产伦人伦偷精品视频| 久久久久九九精品影院| 亚洲不卡免费看| 小说图片视频综合网站| 久久国产乱子伦精品免费另类| 国产欧美日韩精品亚洲av| 一边摸一边抽搐一进一小说| 国产日本99.免费观看| netflix在线观看网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲美女视频黄频| 我要搜黄色片| 久久亚洲真实| 国产av麻豆久久久久久久| 免费无遮挡裸体视频| 一区福利在线观看| 美女xxoo啪啪120秒动态图 | 少妇的逼水好多| 嫩草影视91久久| 亚洲国产欧洲综合997久久,| 性色av乱码一区二区三区2| 精品日产1卡2卡| 99国产精品一区二区三区| 国产av在哪里看| 97超视频在线观看视频| 免费大片18禁| 色综合欧美亚洲国产小说| 亚洲精品456在线播放app | 亚洲人成网站在线播| 色哟哟哟哟哟哟| 午夜福利免费观看在线| 啦啦啦观看免费观看视频高清| 日本五十路高清| 亚洲av免费高清在线观看| 窝窝影院91人妻| 久久国产乱子免费精品| 欧美高清成人免费视频www| 九色国产91popny在线| 又黄又爽又免费观看的视频| 草草在线视频免费看| 毛片女人毛片| 亚洲美女搞黄在线观看 | 深夜精品福利| 嫩草影院新地址| www.熟女人妻精品国产| 国产毛片a区久久久久| 亚洲激情在线av| 国产精品久久久久久人妻精品电影| 精品熟女少妇八av免费久了| 尤物成人国产欧美一区二区三区| 欧美激情久久久久久爽电影| 国产欧美日韩一区二区三| 69av精品久久久久久| 成年女人毛片免费观看观看9| 一本综合久久免费| 午夜免费成人在线视频| 欧美zozozo另类| 亚洲精品456在线播放app | 天堂av国产一区二区熟女人妻| 国产精品人妻久久久久久| 国产免费一级a男人的天堂| 亚州av有码| 丝袜美腿在线中文| 嫁个100分男人电影在线观看| 免费观看精品视频网站| 精品一区二区三区视频在线观看免费| www.999成人在线观看| 日本一二三区视频观看| 国产视频内射| 日韩国内少妇激情av| 国产av在哪里看| www.色视频.com| 天天躁日日操中文字幕| 夜夜夜夜夜久久久久| 中文字幕久久专区| 亚洲成人精品中文字幕电影| 成年女人永久免费观看视频| 国产亚洲欧美98| 欧美高清成人免费视频www| 麻豆久久精品国产亚洲av| 综合色av麻豆| 最近在线观看免费完整版| 天堂影院成人在线观看| 免费在线观看影片大全网站| 精品乱码久久久久久99久播| 可以在线观看的亚洲视频| 国产精品电影一区二区三区| 熟女电影av网| 午夜福利视频1000在线观看| 亚洲最大成人中文| 丝袜美腿在线中文| 亚洲av五月六月丁香网| 深爱激情五月婷婷| 真人做人爱边吃奶动态| 日本精品一区二区三区蜜桃| 中文字幕人妻熟人妻熟丝袜美| 成人国产一区最新在线观看| 美女 人体艺术 gogo| 国产黄片美女视频| 久久天躁狠狠躁夜夜2o2o| 赤兔流量卡办理| 亚洲av成人av| 99热这里只有精品一区| 午夜福利在线在线| 非洲黑人性xxxx精品又粗又长| 小说图片视频综合网站| 午夜免费激情av| 久久久久久久久久成人| 国产高清有码在线观看视频| 99精品在免费线老司机午夜| 九色成人免费人妻av| 欧美中文日本在线观看视频| av国产免费在线观看| 国产精品伦人一区二区| 欧美+亚洲+日韩+国产| 国产精品美女特级片免费视频播放器| 日韩中文字幕欧美一区二区| 男人和女人高潮做爰伦理| 免费看光身美女| 一夜夜www| 国产在线精品亚洲第一网站| 亚洲自偷自拍三级| 国产v大片淫在线免费观看| 99久久久亚洲精品蜜臀av| 嫩草影院新地址| 成人一区二区视频在线观看| 久久婷婷人人爽人人干人人爱| 成年女人毛片免费观看观看9| 成人av一区二区三区在线看| 一进一出抽搐动态| 欧美区成人在线视频| 日韩高清综合在线| 国产精品电影一区二区三区| 在线免费观看的www视频| 成人av在线播放网站| 少妇高潮的动态图| 国产精品一区二区三区四区久久| 天堂网av新在线| 一二三四社区在线视频社区8| 亚洲男人的天堂狠狠| 最新中文字幕久久久久| 亚洲人与动物交配视频| av中文乱码字幕在线| 夜夜夜夜夜久久久久| 国产免费一级a男人的天堂| 麻豆成人午夜福利视频| 亚洲 欧美 日韩 在线 免费| 成年人黄色毛片网站| 两性午夜刺激爽爽歪歪视频在线观看| 色播亚洲综合网| 国产白丝娇喘喷水9色精品| or卡值多少钱| 亚洲欧美日韩高清在线视频| 亚洲av.av天堂| 中文在线观看免费www的网站| 久久欧美精品欧美久久欧美| 在线看三级毛片| 国产野战对白在线观看| 热99re8久久精品国产| 男人舔奶头视频| 日韩欧美精品免费久久 | 久久精品综合一区二区三区| 十八禁国产超污无遮挡网站| 99热只有精品国产| 级片在线观看| 亚洲激情在线av| 非洲黑人性xxxx精品又粗又长| 亚洲成av人片在线播放无| 午夜福利成人在线免费观看| 日韩欧美在线乱码| 丰满的人妻完整版| 99热6这里只有精品| 亚洲电影在线观看av| 美女黄网站色视频| 欧美黄色淫秽网站| 不卡一级毛片| 淫秽高清视频在线观看| 观看免费一级毛片| 美女大奶头视频| 国产精品精品国产色婷婷| 欧美一级a爱片免费观看看| av国产免费在线观看| 免费观看的影片在线观看| 757午夜福利合集在线观看| 国产私拍福利视频在线观看| 午夜福利视频1000在线观看| 久久久久国产精品人妻aⅴ院| 美女高潮喷水抽搐中文字幕| 欧美极品一区二区三区四区| 永久网站在线| 色播亚洲综合网| 成人精品一区二区免费| 亚洲av.av天堂| 51国产日韩欧美| 天天一区二区日本电影三级| 赤兔流量卡办理| 日韩精品青青久久久久久| 一本一本综合久久| 久久久久免费精品人妻一区二区| 神马国产精品三级电影在线观看| 欧美乱妇无乱码| 亚洲,欧美,日韩| 亚洲国产欧美人成| 亚洲avbb在线观看| 午夜福利在线在线| 中国美女看黄片| 啦啦啦观看免费观看视频高清| 婷婷精品国产亚洲av在线| 国产欧美日韩一区二区精品| 色5月婷婷丁香| 婷婷六月久久综合丁香| 国产乱人视频| 最后的刺客免费高清国语| 国内精品一区二区在线观看| 极品教师在线免费播放| 少妇的逼水好多| 亚洲av第一区精品v没综合| 国产精品人妻久久久久久| 俄罗斯特黄特色一大片| 国产黄片美女视频| 丝袜美腿在线中文| 国产精品一区二区性色av| 亚洲欧美清纯卡通| 欧美一区二区精品小视频在线| 最新在线观看一区二区三区| 精品人妻1区二区| 欧美另类亚洲清纯唯美| 亚洲乱码一区二区免费版| 国产美女午夜福利| 丁香欧美五月| 欧美成人免费av一区二区三区| 国产91精品成人一区二区三区| 亚洲欧美精品综合久久99| 国产精品,欧美在线| 精品久久久久久久人妻蜜臀av| 长腿黑丝高跟| 亚洲一区二区三区不卡视频| 少妇裸体淫交视频免费看高清| 波野结衣二区三区在线| 国产人妻一区二区三区在| 亚洲 国产 在线| 久久草成人影院| 村上凉子中文字幕在线| 精品久久久久久久久久久久久| 国产成人福利小说| 真人做人爱边吃奶动态| 丝袜美腿在线中文| 搞女人的毛片| 最近中文字幕高清免费大全6 | 精品熟女少妇八av免费久了| 在线观看av片永久免费下载| 精品一区二区三区视频在线| 精品久久久久久,| 久久久久精品国产欧美久久久| 精品国产三级普通话版| 可以在线观看毛片的网站| av在线天堂中文字幕| 欧美+日韩+精品| 99riav亚洲国产免费| av天堂中文字幕网| 狂野欧美白嫩少妇大欣赏| 欧美日韩黄片免| 美女 人体艺术 gogo| 久久精品夜夜夜夜夜久久蜜豆| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一区二区三区四区久久| 久久久精品大字幕| 国产黄色小视频在线观看| 国产精品久久久久久久电影| 给我免费播放毛片高清在线观看| 神马国产精品三级电影在线观看| 欧美日韩福利视频一区二区| 久久精品久久久久久噜噜老黄 | 午夜激情欧美在线| 大型黄色视频在线免费观看| 亚洲人成网站在线播| 精品久久久久久久久亚洲 | 国内久久婷婷六月综合欲色啪| 一二三四社区在线视频社区8| 最近最新中文字幕大全电影3| 国产欧美日韩一区二区三| 中文字幕高清在线视频| 亚洲狠狠婷婷综合久久图片| 欧美日韩综合久久久久久 | 日本 欧美在线| 亚洲成人久久性| 麻豆国产av国片精品| 一级a爱片免费观看的视频| 十八禁国产超污无遮挡网站| 精品午夜福利在线看| 国产成+人综合+亚洲专区| 综合色av麻豆| 国产老妇女一区| 无人区码免费观看不卡| 午夜免费成人在线视频| 乱码一卡2卡4卡精品| 国产精品精品国产色婷婷| 国产午夜精品论理片| 亚洲五月婷婷丁香| 国产午夜福利久久久久久| 久久久国产成人免费| 色噜噜av男人的天堂激情| 内地一区二区视频在线| 国产精品日韩av在线免费观看| 婷婷亚洲欧美| 一边摸一边抽搐一进一小说| 自拍偷自拍亚洲精品老妇| xxxwww97欧美| 亚洲,欧美精品.| 国产免费男女视频| 男女视频在线观看网站免费| 亚洲av成人不卡在线观看播放网| 老司机福利观看| 国产蜜桃级精品一区二区三区| 99久久99久久久精品蜜桃| 黄色视频,在线免费观看| 一区二区三区激情视频| 一级av片app| 国产一区二区在线观看日韩| 免费一级毛片在线播放高清视频| 久久久久久久久久黄片| 好看av亚洲va欧美ⅴa在| 国产亚洲精品久久久com| 熟女电影av网| 亚洲一区高清亚洲精品| 啦啦啦观看免费观看视频高清| 成年女人永久免费观看视频| 国产精品1区2区在线观看.| 嫁个100分男人电影在线观看| 精品久久久久久,| 亚洲av电影不卡..在线观看| 成人特级av手机在线观看| 精品不卡国产一区二区三区| 午夜精品在线福利| 男女下面进入的视频免费午夜| 国产激情偷乱视频一区二区| 18禁在线播放成人免费| 国内精品美女久久久久久| 午夜精品久久久久久毛片777| 亚洲片人在线观看| 在线播放无遮挡| 岛国在线免费视频观看| 老熟妇仑乱视频hdxx| 亚洲男人的天堂狠狠| 高清日韩中文字幕在线| 99国产精品一区二区三区| 黄色视频,在线免费观看| 婷婷亚洲欧美| 久久久色成人| av福利片在线观看| 日韩亚洲欧美综合| 国产三级中文精品| 亚洲精品456在线播放app | 久久亚洲真实| 最好的美女福利视频网| 欧美成人性av电影在线观看| 国产精品久久电影中文字幕| 亚洲成av人片在线播放无| 自拍偷自拍亚洲精品老妇| 又爽又黄无遮挡网站| 亚洲精品一区av在线观看| 欧美丝袜亚洲另类 | 日本黄色片子视频| 国内揄拍国产精品人妻在线| 亚洲欧美精品综合久久99| 国产蜜桃级精品一区二区三区| 亚洲男人的天堂狠狠| 亚洲,欧美精品.| 久久人妻av系列| 亚洲精品粉嫩美女一区| 精品日产1卡2卡| 黄片小视频在线播放| 精品久久久久久久久久久久久| 欧美精品国产亚洲| 欧美日本视频| 久久久精品欧美日韩精品|