亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于溫度的ET0模型改進(jìn)及在海河流域的適用性

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

        張戴煒 王春穎 韓宇平

        摘 要:為提高HS模型估算海河流域參考作物蒸散量(ET0)精度,基于貝葉斯原理,利用海河流域8個地理分區(qū)1980—1999年29個氣象站的逐日氣象數(shù)據(jù)對其溫度指數(shù)、溫度系數(shù)和溫度常數(shù)進(jìn)行改進(jìn),并以FAO Penman-Monteith(PM)模型為標(biāo)準(zhǔn),利用1980—1999年和2000—2015年數(shù)據(jù)對HS改進(jìn)模型適用性進(jìn)行評價。海河流域HS改進(jìn)模型的溫度系數(shù)、溫度常數(shù)和溫度指數(shù)的范圍分別為5×10-4~8×10-4、28.5~38.0和0.4~0.7。與PM 模型計算結(jié)果對比,HS改進(jìn)模型均方根誤差(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改進(jìn)模型計算結(jié)果與PM計算結(jié)果的時空變化一致,具有較高的計算精度和較好的適應(yīng)性。HS改進(jìn)模型有助于提高海河流域缺少風(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)模型對于原始?xì)庀髷?shù)據(jù)資料的要求較低,在國內(nèi)外得到了廣泛應(yīng)用[3]。但是,以往研究結(jié)果表明HS模型計算結(jié)果存在誤差過大,甚至在部分地區(qū)出現(xiàn)蒸散量為負(fù)數(shù)的情況,因此HS模型需要進(jìn)行修正以提高模型精度[4]。

        國內(nèi)外學(xué)者圍繞HS模型的參數(shù)校正、誤差影響因素分析、適用性評價等問題等開展了大量研究。楊永紅等[5]為了進(jìn)一步提高HS模型在拉薩地區(qū)的適用性,引入了平均相對濕度因子,建立改進(jìn)HS模型并進(jìn)行了評價 。李晨等[6]基于貝葉斯方法改進(jìn)的HS模型在川中丘陵區(qū)不同區(qū)域變異性較小,適應(yīng)性較強,具有較高的計算精度 。胡慶芳等[7]對HS模型在中國7個區(qū)域進(jìn)行全局校正及適用性評價。范文波等[8]采用貝葉斯方法對HS模型在內(nèi)陸河流域進(jìn)行修正及驗證,結(jié)果顯示修正后的HS模型滿足精度要求,且計算簡便,可為內(nèi)陸河流域ET0的計算提供參考。Xu等[9]在瑞士Vaud地區(qū)比較了HS模型和其他5種經(jīng)驗公式的計算精度,發(fā)現(xiàn)HS模型校準(zhǔn)后能準(zhǔn)確計算年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ù)進(jìn)行改進(jìn),以PM模型計算的蒸散量作為評價標(biāo)準(zhǔn)進(jìn)行對比分析,并利用2000—2015年數(shù)據(jù)資料對改進(jìn)的HS模型進(jìn)行適用性評價。

        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改進(jìn)模型

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

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

        2 結(jié)果與分析

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

        海河流域不同分區(qū)HS改進(jìn)模型的溫度系數(shù)、溫度常數(shù)和溫度指數(shù)最優(yōu)值見表1。海河流域HS改進(jìn)模型的參數(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 時空變化趨勢分析

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

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

        2.3 適用性分析

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

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

        3 討 論

        HS模型是基于溫度估算氣象數(shù)據(jù)有限地區(qū)的參考作物蒸散量的可靠方法,降低了ET0計算的氣象數(shù)據(jù)要求,簡化了計算過程。但HS模型經(jīng)驗性較強, Hargreaves 等[16] 推薦當(dāng)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改進(jìn)模型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 模型進(jìn)行了改進(jìn),改進(jìn)后HS模型形式更為復(fù)雜。Allen[17]通過研究認(rèn)為在高程超過1 500 m的地區(qū)Hargreaves推薦值并不合理;Annandale等[18] 認(rèn)為HS模型參數(shù)取值應(yīng)考慮大氣壓因素;Samani等[14] 認(rèn)為在研究區(qū)域為北緯7°—50°時,修正參數(shù)應(yīng)依據(jù)研究區(qū)域的溫度梯度確定。本研究改進(jìn)HS模型是在只考慮溫度的大背景下進(jìn)行的,模型誤差存在空間差異性,對于HS改進(jìn)模型誤差較大的地點,應(yīng)進(jìn)一步考慮高程、風(fēng)速、濕度等因素改進(jìn)HS模型。

        4 結(jié) 語

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

        參考文獻(xiàn):

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

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

        [7] 胡慶芳,楊大文,王銀堂,等.Hargreaves公式的全局校正及適用性評價[J].水科學(xué)進(jìn)展,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在帕金森病模型中的表達(dá)及分布
        函數(shù)模型及應(yīng)用
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        3D打印中的模型分割與打包
        日韩一级黄色片一区二区三区| 久久www色情成人免费观看| 精品国产18久久久久久| 无码丰满少妇2在线观看| 国产人在线成免费视频麻豆| 蜜桃一区二区免费视频观看| 日韩国产一区二区三区在线观看| 亚洲中文字幕精品乱码2021| 亚洲精品无人区| 国产亚洲欧美日韩综合一区在线观看| 美女叉开双腿让男人插| 日本少妇又色又紧又爽又刺激| 成人影院在线观看视频免费| 老熟女的中文字幕欲望 | 玖玖资源站无码专区| 精品人妻av一区二区三区不卡| 日本人妻三级在线观看| 国产女优一区在线观看| 久久精品国产免费观看三人同眠| 人人爽人人爽人人片av| 久久精品国产亚洲av大全| 日产精品久久久久久久| 视频女同久久久一区二区三区 | 中文国产乱码在线人妻一区二区| 日本一本之道高清不卡免费| 丁香综合网| 中文字幕视频二区三区| 国产黑丝美女办公室激情啪啪 | 国产精品无码无卡无需播放器| 国产成人精品无码一区二区老年人| 两个人看的www高清视频中文| 一区二区三无码| 亚洲av男人的天堂在线| 亚洲男人天堂黄色av| 海角国精产品一区一区三区糖心| 亚洲色图视频在线观看网站| 国产精品不卡免费版在线观看| 亚洲丁香婷婷久久一区二区| 国产精品51麻豆cm传媒| 亚洲饱满人妻视频| 国产成人久久精品77777综合|