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

    參考作物蒸散量計算方法在西寧的適用性分析

    2018-07-09 12:32:02孟昭睿張艷輝巴瑞壽祁澤學
    中國農(nóng)村水利水電 2018年6期
    關鍵詞:西寧市西寧青海省

    吳 萍,孟昭睿,吳 莉,張艷輝,巴瑞壽,祁澤學

    (1.青海省環(huán)境地質勘查局,西寧 810007;2.青海省環(huán)境地質重點實驗室,西寧 810007;3.邯鄲市環(huán)境保護局館陶縣分局,河北 邯鄲 057750;4.館陶縣安全生產(chǎn)監(jiān)督管理局,河北 邯鄲 057750)

    1948年Penman在英格蘭南部地區(qū)首先提出無水汽水平輸送情況下的參考作物蒸散量計算公式即1948 Penman 法(以下簡稱1948 Pen),1956年對公式進行了改進。Monteith在Penman等研究的基礎上提出了作物蒸散量計算的阻力模式,即Penman-Monteith(以下簡稱PM)公式,被認為是計算參考作物蒸散量(Reference Crop Evapotranspiration,ET0)[1]的標準化方法和ET0實測資料缺乏地區(qū)評價其他方法的標準。過去幾十年,學者已提出許多ET0簡便估算公式,如Hargreaves-Samani(Har)法[2]、McCloud(MC)法[3]和1948 Penman法 (1948 Pen)方法[4]等,由于輸入?yún)?shù)較少而得到廣泛應用[5-9],但對高海拔寒旱地區(qū)的研究較少。本文以西寧站1970-2016年的氣象資料為例,利用Penman-Monteith方法為衡量標準,分析對比Hargreaves-Samani法、McCloud法和1948 Penman 法在青海省西寧地區(qū)的適用性,旨在為青海西寧市計算ET0方法的選取提供理論參考。

    1 數(shù)據(jù)與方法

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

    氣象站西寧站(101°77′E,36°62′N)位于青藏高原東部,平均海拔2 261 m該處屬于大陸性高原半干旱氣候,其主要特點是:氣壓低,太陽輻射強,日照時間長,日夜溫差大,常年氣溫偏低,年降水量小,年蒸發(fā)量大,無霜期短,冰凍期長,冬無嚴寒,夏無酷暑,年平均降水量380 mm,蒸發(fā)量1 363.6 mm,年平均氣溫7.6 ℃,最高氣溫34.6 ℃,最低氣溫零下18.9 ℃,夏季平均氣溫17~19 ℃。主要氣象資料為:1970-2016年的逐日平均氣壓、最高氣溫、最低氣溫、平均相對濕度、日照時數(shù)、平均風速等,計算逐日參考作物蒸散量。

    1.2 參考作物蒸散量的計算

    (1)Penman-Monteith(PM)公式。具體計算公式為:

    (1)

    式中:ET0PM為Penman-Monteith公式計算得到的日參考作物蒸散量,mm/d;Δ為飽和水氣壓與溫度曲線的斜率,kPa/℃;Rn為參考作物冠層表面凈輻射,MJ/(m2·d);G為土壤熱通量,MJ/(m2·d);γ為干濕表常數(shù),kPa/℃;T為2 m高處的日平均氣溫,℃;U2為2 m高處的風速,m/s;es為飽和水汽壓,kPa;ea為實際水汽壓,kPa。

    (2)Hargreaves-Samani(H-S)公式。在美國西北部較干旱的氣候條件下建立的,它僅需要月最高、最低氣溫,具體計算公式為:

    (2)

    式中:ET0H為H-Sgreaves- Samani公式計算得到的日參考作物蒸散量,mm/d;Tmax為最高溫度,℃;Tmin為最低溫度,℃;T為平均溫度,℃;Ra為外空輻射,MJ/(m2·d);λ為水汽化潛熱,λ=2.45 MJ/kg;C0為轉換系數(shù),當Ra以mm/d為單位時,C0=0.002 3,以MJ/(m2·d)為單位時,C0=0.000 939。

    (3)McCloud(MC)公式?;谌掌骄鶜鉁兀旹T0為溫度的指數(shù)函數(shù),最初用于估算草坪草的潛在蒸散量,將原來公式的各變量單位換算成國際制得到的公式為:

    ET0M=KW1.8T

    (3)

    式中:ET0M為McCloud法計算得到的日參考作物蒸散量,mm/d;K=0.254;W=1.07。

    (4)1948 Penman(1948 Pen)公式。1948 Penman 法是依據(jù)能量平衡和紊流擴散原理導出的計算參考作物騰發(fā)量的方法。該方法具有堅實的理論基礎, 比以上的經(jīng)驗公式更科學合理, 目前仍為濕潤下墊面蒸散計算的主要方法。該方法需要氣溫、相對濕度、日照時數(shù)、風速等資料計算參考作物騰發(fā)量, 具體計算公式:

    (4)

    式中:ET048為1948 Penman 法計算得到的日參考作物蒸散量,mm/d。

    1.3 評價方法

    平均絕對偏差(Mean Absolute Error,MAE):

    (5)

    平均相對百分誤差(Mean Relative Percent Error,MRPE):

    (6)

    相對均方差(Relative Mean Square Error,RMSE):

    (7)

    式中:ET0PM(i)為PM公式計算的參考作物蒸散量系列值;ET0x(i)表示第二種方法計算的參考作物蒸散量系列值;MAE、MRPE、RMSE越接近0,效果越好。

    2 適用性分析

    2.1 逐日ET0序列評價分析

    對西寧氣象站點1970-2016年的氣象站資料采用3種方法進行逐日ET0計算,ET0以PM方法為比較標準,繪出多年平均逐日ET0變化曲線及其偏差(見圖1~圖3)。

    圖1 H-S方法與PM方法多年平均逐日ET0計算結果對比圖

    (1)H-S方法的逐日ET0變化趨勢近似為單調上升和下降的拋物線;與PM方法相比,H-S方法的ET0結果全年偏高于PM方法,上升段為1-7月,下降段為7-12月;尖峰7月12日的ET0極大值偏大1.16 mm/d,1-3月基本一致,絕對誤差在0.4 mm/d以上的天數(shù)在50%多集中在4-9月;整個ET0曲線的絕對偏差為-1.16~0.01 mm/d、相對偏差在10%以內的天數(shù)在15%多集中在1-3月。

    圖2 MC方法與PM方法多年平均逐日ET0計算結果對比圖

    圖3 1948 Pen方法與PM方法多年平均逐日ET0計算結果對比圖

    (2)MC方法的逐日ET0變化趨勢類似于H-S從1月單調上升至7月但峰值后單調下降;與PM方法相比,其ET0結果全年偏低于PM方法,尖峰7月20日的ET0極大值偏小1.28 mm/d,絕對誤差在1.5 mm/d以上的天數(shù)在50%多集中在4-7月;整個ET0曲線的絕對偏差為0.4~2.5 mm/d、相對偏差在全年都在25%以上,誤差較大。

    (3)1948 Pen方法的逐日ET0變化趨勢也近似為單調拋物線;與PM方法相比,1948 Pen方法的ET0結果全年偏高于PM方法結果從圖上看1948 Pen與PM全年變化趨勢有極好的一致性;尖峰7月15日的ET0極大值偏大0.68 mm/d,整個ET0曲線的絕對偏差為-0.36~-0.89 mm/d。相對偏差為10%~80%,其中40%以上的多分布于11-2月。

    2.2 逐月ET0序列評價分析

    由圖4可看出,3種方法各站的多年平均逐月ET0估算結果,年內都從1-12月呈拋物線形分布,最高值出現(xiàn)在7月,最低值出現(xiàn)在12月或1月;MC法估算結果與PM法各月ET0估算結果的絕對誤差最大,偏低12.5~73.0 mm/月,變化步調最不一致;H-S法月ET0估算結果無顯著差異的月份主要集中在12月-次年3月,4-10月H-S法月ET0估算結果高于PM法結果。H-S法月ET0估算結果在夏季明顯偏高,而且較PM法結果差異大,冬季計算結果普遍偏低,差異較小,這與很多研究得出的結論相似[10-13],絕對誤差為3.4~28.8 mm/月,MRPE和RMSE最大值為分別為35.8%、17.0;從圖4看1948 Pen法月ET0估算結果與PM法變化一致性良好,絕對誤差為11.1~20.0 mm/月,MRPE和RMSE最大值為分別為64.5%、14.1 mm/月;從總體看,1948 Pen法表現(xiàn)略優(yōu)于H-S法,MC法表現(xiàn)最差(見表1)。

    2.3 ET0相關性分析

    利用西寧站點的氣象資料,采用3種方法計算得到的多年平均逐日ET0值系列分別與PM方法的結果進行相關性分析,結果見圖5。由圖5可知,3種方法得到的趨勢線方程和相關系數(shù),其中H-S、1948 Pen法得相關系數(shù)達到0.98以上,遠遠超過了信度為0.01的臨界值,最差的是MC方法,只有0.819 2。

    圖4 西寧多年平均逐月ET0計算結果對比圖

    月份H-SMAEMRPERMSEMCMAEMRPERMSE1948 PenMAEMRPERMSE1 3.5 18.5 3.1 15.8 79.8 8.1 11.4 57.9 6.4 2 3.4 10.8 3.7 26.6 82.5 13.6 12.0 37.4 7.0 3 5.5 9.4 5.7 48.1 79.7 24.5 17.0 28.2 10.1 4 10.6 12.0 9.6 65.6 73.0 33.5 20.0 22.3 12.9 5 17.5 15.7 11.5 73.0 64.3 37.3 21.1 18.6 14.1 6 24.2 20.6 14.3 63.6 53.3 32.6 19.4 16.3 13.0 7 28.8 23.7 17.0 49.0 39.8 26.2 18.6 15.1 13.5 8 22.9 20.9 14.6 42.0 37.5 22.6 17.5 15.7 13.6 9 15.5 21.1 10.5 34.9 46.4 18.3 14.0 18.6 10.2 10 12.4 25.3 8.0 28.6 57.2 14.8 13.5 27.0 8.6 11 9.1 34.3 5.7 17.6 65.0 9.1 12.3 45.5 7.0 12 6.0 35.8 4.1 12.5 71.7 6.5 11.1 64.5 6.2

    圖5 不同計算方法與PM方法的多年平均逐日ET0相關關系圖

    3 方法修正及精度分析

    通過圖5相關性分析,在青海省西寧市使用1948 Penman法計算參考作物蒸散量的值乘以一個常數(shù)就能夠近似代替Penman-Monteith法的計算值,應根據(jù)不同時段對α重新修正,分別使用PM法和1948 Pen法對青海西寧市1970-2016年47 a實測到的日氣象數(shù)據(jù)進行計算后得到整年內的逐日平均ET048和ET0PM,對這些數(shù)據(jù)按月份分組再對每組中的ET048和ET0PM計算出相應的回歸方程從而得到合適的折算系數(shù),在青海省西寧市使用1948 Penman法估算短時段參考作物蒸散量時,應在不同的時段內使用不同的α值,具體的α范圍值見表2。

    表2 青海省西寧市1948 Penman法重新修正的α值

    西寧市站采用表2提供的修正α值得到1948 Penman法修正后的逐月ET048估算結果與PM法結果進行誤差分析,結果見表3。對比表3和表1可以看出,采用修正參數(shù)后的1948 Penman法在西寧市估算的ET0誤差明顯減小,如12月α值取0.61,MAE由原來的11.10 mm/月減小到0.02 mm/月,MRPE由64.5%減小到0.08%,RMSE由6.2 mm/月減小到0.03 mm/月;7月α值取0.87,MAE由原來的18.6 mm/月減小到0.21 mm/月,MRPE由15.1%減小到0.17%,RMSE由13.5 mm/月減小到2.7 mm/月,所以修正后精度較高。

    表3 修正后1948 Penman法與PM法計算結果誤差

    4 結 語

    (1)以PM方法作為標準,對其他方法進行評價,結果表明:在西寧站通過對3種方法計算的逐月ET0的比較,西寧地區(qū)3種方法變化趨勢基本一致,年內都從1-12月呈拋物線形分布,最高值出現(xiàn)在7月,最低值出現(xiàn)在12月

    或1月,與PM法結果相比,絕對誤差和相對均方差隨ET0值的增大而增大;1948 Pen和H-S法估算結果優(yōu)于MC法估算結果。

    (2)在西寧H-S法在4-10月誤差相對較大,當缺少輻射和風速資料時,按月修正公式的系數(shù),可得到精度相對較高的估算公式;1948 Pen法與PM法相關系數(shù)較高,一致性好,可以對1948 Pen修正式的計算值乘以一個常數(shù)就能夠近似代替PM公式的計算值。

    (3)1948 Penman法估算結果與標準值最接近,1948 Pen法計算參考作物蒸散量的值乘以一個常數(shù)就能夠近似代替Penman-Monteith法,本文按照不同時段對α進行修正,得到1948 Penman法適合青海省西寧市不同時段的α值,大大提高了青海省西寧市參考作物蒸散量的估算精度。

    參考文獻:

    [1] Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration guidelines for computing crop water requirements[C]∥ FAO Irrigation and Drainage Paper (No.56).Rome:FAO,1998:15-64.

    [2] Hargreaves G H,F ASCE,Allen RG. History and evaluation of Hargreaves evapotran-spiration equation[J].Journal of Irrigation and Drainage Engineering,2003,129(1):53-63.

    [3] Jennifer M Jacobs. Evaluation of reference evapotranspiration methodologies and AFSIRS crop water use simulation model(final report)[M].Florida: St.Johns River Water Management District Palatka, 2001:4-5.

    [4] Batchelo r C H.The accuracy of evapotranspiration estimated with the modified Penman equation[J]. Irrigation Science,1984,5(4):223-233.

    [5] 劉 鈺,蔡林根.參照騰發(fā)量的新定義及計算方法對比[J].水利學報,1997,(6):24-26.

    [6] 胡順軍,潘 渝,康紹忠,等.Penman-Monteith與Penman修正式計算塔里木盆地參考作物潛在騰發(fā)量比較[J].農(nóng)業(yè)工程學報,2005,(6):30-35.

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

    [8] 左德鵬,徐宗學,劉兆飛.基于氣溫的潛在蒸散發(fā)量估算方法在我國西北干旱地區(qū)的應用比較[J].干旱區(qū)資源與環(huán)境,2009,(10):123-131.

    [9] 丁加麗,彭世彰,徐俊增,等.基于溫度資料的參考作物蒸發(fā)蒸騰量計算方法[J].河海大學學報(自然科學版),2007,(6):633-637.

    [10] Karl Vanderlinden, J V Giraldez, M V Meirvenne. Assessing reference evapotranspiration by the Hargreaves Method in Southern Spain [J].Journal of Irrigation and Drainage Engineering, 2004,130(4):184-191.

    [11] 劉曉英,李玉中,王慶鎖.幾種基于溫度的參考作物蒸散量計算方法的評價[J].農(nóng)業(yè)工程學報,2006,(6):12-18.

    [12] 劉 鈺,L S Pereira.氣象數(shù)據(jù)缺測條件下參照騰發(fā)量的計算方法[J].水利學報,2001,(3):11-17.

    [13] 王新華,郭美華,徐中民.分別利用Hargreaves和PM公式計算西北干旱區(qū)ET0的比較[J].農(nóng)業(yè)工程學報,2006,(10):21-25.

    猜你喜歡
    西寧市西寧青海省
    西寧市人民政府大事記 2022年6月
    西寧市人民政府人事任免 2022年6月
    西寧市人民政府大事記 2022年7月
    西寧市人民政府人事任免 2022年7月
    Dynamical signatures of the one-dimensional deconfined quantum critical point
    輕輕松松聊漢語——“中國夏都”西寧
    金橋(2018年7期)2018-09-25 02:28:28
    青海西寧蘭州格爾木往來更暢通
    石油瀝青(2018年5期)2018-03-23 04:49:19
    落實“四個扎扎實實”展現(xiàn)開行使命擔當——國家開發(fā)銀行青海省分行
    青海政報(2017年20期)2017-03-16 02:38:19
    青海省交通運輸行業(yè)數(shù)據(jù)中心節(jié)能探索
    青海省人民政府辦公廳關于轉發(fā)省科協(xié)青海省自然科學優(yōu)犁匆圈獨砰審和勤加祛的通知
    青海政報(2016年21期)2016-02-28 12:26:23
    久久久精品区二区三区| 成在线人永久免费视频| 免费人妻精品一区二区三区视频| 成人影院久久| 久久久精品区二区三区| 亚洲久久久国产精品| 久久亚洲真实| 色综合欧美亚洲国产小说| 99久久99久久久精品蜜桃| 久久国产精品影院| 欧美午夜高清在线| 交换朋友夫妻互换小说| 国产成人啪精品午夜网站| 国产亚洲精品一区二区www | 午夜免费成人在线视频| 亚洲中文av在线| 夜夜爽天天搞| 午夜成年电影在线免费观看| 激情在线观看视频在线高清 | 午夜福利在线免费观看网站| 亚洲伊人久久精品综合| 99国产精品99久久久久| 日韩制服丝袜自拍偷拍| 欧美日本中文国产一区发布| 日本一区二区免费在线视频| 中文字幕高清在线视频| 国产日韩欧美亚洲二区| 久久人妻熟女aⅴ| 亚洲人成伊人成综合网2020| 天天躁日日躁夜夜躁夜夜| 另类精品久久| 国产欧美日韩一区二区三| 99国产精品一区二区蜜桃av | 在线十欧美十亚洲十日本专区| av天堂在线播放| 国产亚洲欧美在线一区二区| 精品久久久久久电影网| 王馨瑶露胸无遮挡在线观看| 丁香六月天网| 亚洲国产精品一区二区三区在线| 免费久久久久久久精品成人欧美视频| 国产精品98久久久久久宅男小说| 亚洲视频免费观看视频| 青草久久国产| 搡老乐熟女国产| 亚洲欧美日韩另类电影网站| 欧美老熟妇乱子伦牲交| 欧美日韩视频精品一区| 中文字幕高清在线视频| 大型黄色视频在线免费观看| 黄色视频,在线免费观看| 国产一区二区在线观看av| 91精品国产国语对白视频| 久久久久久人人人人人| 亚洲av日韩精品久久久久久密| 亚洲一卡2卡3卡4卡5卡精品中文| 一区二区三区国产精品乱码| 新久久久久国产一级毛片| 黑丝袜美女国产一区| 国产精品秋霞免费鲁丝片| 少妇猛男粗大的猛烈进出视频| 久久久久久久大尺度免费视频| 十八禁网站网址无遮挡| 免费不卡黄色视频| 亚洲va日本ⅴa欧美va伊人久久| videos熟女内射| 国产精品av久久久久免费| 99久久精品国产亚洲精品| 亚洲av片天天在线观看| 亚洲国产毛片av蜜桃av| 亚洲精华国产精华精| 国产精品免费大片| 国产精品1区2区在线观看. | 天天添夜夜摸| 少妇的丰满在线观看| 中文字幕精品免费在线观看视频| 欧美变态另类bdsm刘玥| 美女主播在线视频| 搡老熟女国产l中国老女人| 日本五十路高清| 久久国产精品大桥未久av| aaaaa片日本免费| 久久精品亚洲精品国产色婷小说| 色在线成人网| 午夜福利在线观看吧| 一级毛片电影观看| 变态另类成人亚洲欧美熟女 | 大片免费播放器 马上看| 日韩制服丝袜自拍偷拍| 99国产精品一区二区蜜桃av | 亚洲熟女毛片儿| 欧美激情 高清一区二区三区| 免费人妻精品一区二区三区视频| 日本av手机在线免费观看| 国产精品久久久久久精品古装| 精品人妻熟女毛片av久久网站| 午夜精品久久久久久毛片777| 91精品国产国语对白视频| 精品人妻1区二区| 亚洲九九香蕉| 中文字幕最新亚洲高清| 亚洲专区中文字幕在线| 黄频高清免费视频| 国产高清国产精品国产三级| 亚洲 国产 在线| 久久中文字幕一级| 最近最新中文字幕大全电影3 | 亚洲国产毛片av蜜桃av| 精品福利永久在线观看| 最新美女视频免费是黄的| 国产精品久久久久久精品古装| 美女高潮到喷水免费观看| 无人区码免费观看不卡 | 一边摸一边抽搐一进一小说 | 国产一区二区三区视频了| 在线天堂中文资源库| 99热国产这里只有精品6| 在线观看舔阴道视频| 久久久国产欧美日韩av| 一区二区三区精品91| 免费在线观看日本一区| 日韩成人在线观看一区二区三区| 视频在线观看一区二区三区| 一级毛片女人18水好多| 人人妻,人人澡人人爽秒播| 国内毛片毛片毛片毛片毛片| 欧美精品av麻豆av| 99国产极品粉嫩在线观看| 久久中文字幕人妻熟女| 久久久久精品国产欧美久久久| 一区二区三区激情视频| 国产男女超爽视频在线观看| 午夜久久久在线观看| 又紧又爽又黄一区二区| 一二三四社区在线视频社区8| 精品一区二区三区四区五区乱码| 97人妻天天添夜夜摸| 91麻豆精品激情在线观看国产 | 亚洲av国产av综合av卡| netflix在线观看网站| 日韩成人在线观看一区二区三区| 久久九九热精品免费| 国产成人精品久久二区二区免费| 国产精品久久久av美女十八| 国产福利在线免费观看视频| 国产片内射在线| 蜜桃国产av成人99| 国产成人精品无人区| 最新美女视频免费是黄的| 纯流量卡能插随身wifi吗| 国产精品自产拍在线观看55亚洲 | 国产精品久久久久久精品电影小说| 国产精品国产av在线观看| tocl精华| 免费在线观看视频国产中文字幕亚洲| 人妻一区二区av| 国产极品粉嫩免费观看在线| 国产精品久久久久久精品电影小说| 免费在线观看完整版高清| 亚洲人成电影免费在线| 亚洲午夜精品一区,二区,三区| 午夜福利,免费看| 一本综合久久免费| 人妻 亚洲 视频| 香蕉国产在线看| 黄色毛片三级朝国网站| 国产精品久久久av美女十八| 中文字幕人妻熟女乱码| 国产精品久久久久成人av| 日韩欧美一区二区三区在线观看 | 成年女人毛片免费观看观看9 | 国产成人欧美| 91av网站免费观看| 久久午夜综合久久蜜桃| 国产日韩欧美在线精品| 美女主播在线视频| 少妇 在线观看| 国产三级黄色录像| 欧美另类亚洲清纯唯美| 成人av一区二区三区在线看| 满18在线观看网站| 日日爽夜夜爽网站| 欧美午夜高清在线| 精品久久蜜臀av无| 狠狠精品人妻久久久久久综合| 亚洲七黄色美女视频| 国产精品一区二区免费欧美| 日韩视频在线欧美| 美女主播在线视频| 黑人巨大精品欧美一区二区mp4| 国产欧美亚洲国产| 正在播放国产对白刺激| 天天躁狠狠躁夜夜躁狠狠躁| 久久午夜亚洲精品久久| 一边摸一边抽搐一进一出视频| 亚洲男人天堂网一区| 精品一区二区三区四区五区乱码| 精品国产超薄肉色丝袜足j| 免费高清在线观看日韩| 怎么达到女性高潮| 天天躁狠狠躁夜夜躁狠狠躁| 国产欧美日韩综合在线一区二区| av免费在线观看网站| 欧美激情极品国产一区二区三区| 正在播放国产对白刺激| 精品亚洲成a人片在线观看| 中国美女看黄片| 国产1区2区3区精品| 国产又色又爽无遮挡免费看| 操美女的视频在线观看| 欧美日韩黄片免| 国内毛片毛片毛片毛片毛片| 亚洲精品av麻豆狂野| 国产真人三级小视频在线观看| 国产色视频综合| 中文字幕精品免费在线观看视频| 国产亚洲欧美在线一区二区| 久久人人爽av亚洲精品天堂| 少妇精品久久久久久久| 国产成人av教育| 国产精品一区二区精品视频观看| 中文字幕人妻丝袜制服| 日日摸夜夜添夜夜添小说| 久久精品熟女亚洲av麻豆精品| 麻豆国产av国片精品| 在线 av 中文字幕| 精品久久久久久久毛片微露脸| 国产片内射在线| 免费看a级黄色片| 桃红色精品国产亚洲av| 国产精品久久久久久精品电影小说| 国产三级黄色录像| 久久午夜综合久久蜜桃| 精品亚洲乱码少妇综合久久| 久久午夜亚洲精品久久| 女同久久另类99精品国产91| 黄色视频不卡| 国产精品一区二区免费欧美| 另类亚洲欧美激情| 夜夜夜夜夜久久久久| 欧美成人午夜精品| 欧美精品一区二区大全| 国产成人精品久久二区二区免费| 国产在线视频一区二区| av网站在线播放免费| 中文字幕制服av| 久久天躁狠狠躁夜夜2o2o| 精品少妇内射三级| 十八禁网站免费在线| 午夜激情久久久久久久| 国产男靠女视频免费网站| 香蕉国产在线看| 国产精品久久久久久精品电影小说| 亚洲av日韩精品久久久久久密| 伦理电影免费视频| 精品一区二区三卡| 国产精品美女特级片免费视频播放器 | 久久久精品免费免费高清| 久久久精品国产亚洲av高清涩受| 久久久久久久大尺度免费视频| 久久免费观看电影| 手机成人av网站| 视频区欧美日本亚洲| 高清在线国产一区| 国产高清激情床上av| 日本撒尿小便嘘嘘汇集6| 成人特级黄色片久久久久久久 | 中文字幕高清在线视频| 99在线人妻在线中文字幕 | 午夜两性在线视频| 男女边摸边吃奶| 91九色精品人成在线观看| 王馨瑶露胸无遮挡在线观看| 久久av网站| 99久久精品国产亚洲精品| 欧美日韩黄片免| 日本a在线网址| 777米奇影视久久| 国产无遮挡羞羞视频在线观看| 91成年电影在线观看| aaaaa片日本免费| av不卡在线播放| 1024香蕉在线观看| 高清黄色对白视频在线免费看| 人妻 亚洲 视频| 男女之事视频高清在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产伦理片在线播放av一区| 国产成人av教育| 丰满少妇做爰视频| 91成人精品电影| 亚洲七黄色美女视频| 国产高清videossex| 天堂俺去俺来也www色官网| 国产免费av片在线观看野外av| 久久性视频一级片| 淫妇啪啪啪对白视频| 美女视频免费永久观看网站| 女人久久www免费人成看片| 三级毛片av免费| 亚洲专区字幕在线| 日韩精品免费视频一区二区三区| 18在线观看网站| 丁香六月欧美| 天天影视国产精品| 久久久久久久久久久久大奶| 两个人看的免费小视频| 日韩大码丰满熟妇| 亚洲 欧美一区二区三区| 1024香蕉在线观看| 久久午夜综合久久蜜桃| 宅男免费午夜| 黄色丝袜av网址大全| 男女无遮挡免费网站观看| 欧美黄色淫秽网站| 国产男靠女视频免费网站| 操出白浆在线播放| 亚洲专区中文字幕在线| 女警被强在线播放| 91字幕亚洲| 一区二区三区精品91| 国产精品98久久久久久宅男小说| 日韩人妻精品一区2区三区| 夜夜骑夜夜射夜夜干| 91成年电影在线观看| 国产精品电影一区二区三区 | 欧美亚洲 丝袜 人妻 在线| 757午夜福利合集在线观看| 黑丝袜美女国产一区| 真人做人爱边吃奶动态| 久久精品人人爽人人爽视色| 亚洲国产欧美网| 国产一区二区三区综合在线观看| 50天的宝宝边吃奶边哭怎么回事| 久久av网站| 久久ye,这里只有精品| 两个人免费观看高清视频| 91麻豆av在线| 亚洲国产看品久久| av网站在线播放免费| 日本欧美视频一区| 精品欧美一区二区三区在线| 亚洲国产看品久久| 国产精品欧美亚洲77777| 搡老岳熟女国产| 免费不卡黄色视频| tube8黄色片| 成人国产av品久久久| 在线观看免费午夜福利视频| 天天躁夜夜躁狠狠躁躁| 老司机午夜十八禁免费视频| 极品少妇高潮喷水抽搐| 18禁裸乳无遮挡动漫免费视频| 极品少妇高潮喷水抽搐| 国产成人欧美在线观看 | 国产有黄有色有爽视频| 不卡av一区二区三区| 精品国产乱码久久久久久男人| 热re99久久国产66热| 视频区欧美日本亚洲| 久久国产亚洲av麻豆专区| 国产午夜精品久久久久久| 色视频在线一区二区三区| 亚洲精品一二三| 高清在线国产一区| 午夜福利,免费看| 国产精品一区二区免费欧美| 日本黄色日本黄色录像| 亚洲精品av麻豆狂野| 岛国毛片在线播放| 在线观看免费高清a一片| 在线天堂中文资源库| 国产av国产精品国产| 亚洲精品美女久久av网站| 亚洲 欧美一区二区三区| 亚洲,欧美精品.| 午夜久久久在线观看| 日韩欧美一区视频在线观看| 亚洲欧洲精品一区二区精品久久久| 国产欧美日韩一区二区三| 成人影院久久| 深夜精品福利| 久久 成人 亚洲| 亚洲精品国产精品久久久不卡| 精品少妇久久久久久888优播| 18禁国产床啪视频网站| 国产在线精品亚洲第一网站| 美女高潮喷水抽搐中文字幕| 亚洲精品久久午夜乱码| 国产成人精品在线电影| 亚洲欧美一区二区三区久久| 日韩视频一区二区在线观看| 丰满少妇做爰视频| 亚洲欧美一区二区三区黑人| 亚洲精品一卡2卡三卡4卡5卡| 免费看十八禁软件| 精品国产亚洲在线| 99香蕉大伊视频| 伊人久久大香线蕉亚洲五| 极品人妻少妇av视频| 十八禁高潮呻吟视频| 在线av久久热| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧洲日产国产| kizo精华| 啦啦啦在线免费观看视频4| 人人妻人人爽人人添夜夜欢视频| 最新美女视频免费是黄的| 黄片小视频在线播放| 日本av手机在线免费观看| 人人妻人人爽人人添夜夜欢视频| 国产1区2区3区精品| 精品久久久久久久毛片微露脸| 女人爽到高潮嗷嗷叫在线视频| 亚洲九九香蕉| 亚洲色图av天堂| 久久天躁狠狠躁夜夜2o2o| 精品少妇一区二区三区视频日本电影| 欧美日韩视频精品一区| 成人国产一区最新在线观看| 桃花免费在线播放| 亚洲国产看品久久| 国产一区二区 视频在线| 久久久水蜜桃国产精品网| 亚洲色图综合在线观看| 午夜福利在线免费观看网站| 免费在线观看视频国产中文字幕亚洲| 最近最新免费中文字幕在线| 精品一区二区三区视频在线观看免费 | 大型黄色视频在线免费观看| av片东京热男人的天堂| 国产成人系列免费观看| 法律面前人人平等表现在哪些方面| 亚洲自偷自拍图片 自拍| 国产三级黄色录像| 王馨瑶露胸无遮挡在线观看| 国产成人精品久久二区二区91| 久久精品亚洲熟妇少妇任你| 国产黄色免费在线视频| 国产成人欧美| 国产麻豆69| 免费少妇av软件| 法律面前人人平等表现在哪些方面| 国产91精品成人一区二区三区 | 一夜夜www| 久久精品亚洲熟妇少妇任你| 高清毛片免费观看视频网站 | 国产免费视频播放在线视频| 天天添夜夜摸| 黄色视频在线播放观看不卡| 午夜福利在线观看吧| 国产一区二区三区视频了| 久久香蕉激情| 国产野战对白在线观看| 午夜久久久在线观看| 国产免费av片在线观看野外av| 色视频在线一区二区三区| 国产午夜精品久久久久久| 99国产精品99久久久久| 国产成人精品久久二区二区91| 精品福利观看| 久久久国产成人免费| 伊人久久大香线蕉亚洲五| 精品久久蜜臀av无| 大型黄色视频在线免费观看| 大香蕉久久成人网| 香蕉国产在线看| netflix在线观看网站| 十八禁网站免费在线| 亚洲欧洲精品一区二区精品久久久| 巨乳人妻的诱惑在线观看| 三级毛片av免费| 999精品在线视频| 99香蕉大伊视频| 亚洲专区国产一区二区| 国产有黄有色有爽视频| 中文字幕av电影在线播放| 国产在线视频一区二区| 久久ye,这里只有精品| 另类亚洲欧美激情| 国产精品欧美亚洲77777| 一区福利在线观看| 一进一出抽搐动态| 国产日韩欧美视频二区| 免费观看人在逋| 王馨瑶露胸无遮挡在线观看| 91成人精品电影| 国产极品粉嫩免费观看在线| 久久国产精品影院| 交换朋友夫妻互换小说| 天天躁日日躁夜夜躁夜夜| 亚洲熟妇熟女久久| 高清在线国产一区| tube8黄色片| aaaaa片日本免费| 悠悠久久av| 最新的欧美精品一区二区| 少妇猛男粗大的猛烈进出视频| 亚洲精品美女久久av网站| 久久人妻av系列| 人妻 亚洲 视频| av一本久久久久| 久久久水蜜桃国产精品网| av不卡在线播放| 嫁个100分男人电影在线观看| 亚洲 国产 在线| 免费久久久久久久精品成人欧美视频| 亚洲av美国av| 天天躁日日躁夜夜躁夜夜| 视频在线观看一区二区三区| 91字幕亚洲| www.精华液| 国产三级黄色录像| 亚洲一区中文字幕在线| 99香蕉大伊视频| 丁香欧美五月| 国产在线精品亚洲第一网站| 99精国产麻豆久久婷婷| 亚洲专区国产一区二区| 精品一区二区三区四区五区乱码| 亚洲精品久久成人aⅴ小说| 国精品久久久久久国模美| 一夜夜www| 国产一区二区三区视频了| 国产区一区二久久| 国产伦理片在线播放av一区| 搡老熟女国产l中国老女人| 丝袜人妻中文字幕| 99国产极品粉嫩在线观看| av一本久久久久| 精品久久久久久久毛片微露脸| 精品福利观看| 久久久精品国产亚洲av高清涩受| 日韩欧美一区视频在线观看| 中文字幕高清在线视频| 1024视频免费在线观看| 国产深夜福利视频在线观看| 亚洲欧洲精品一区二区精品久久久| 国产精品.久久久| 曰老女人黄片| 无人区码免费观看不卡 | 久久久久国产一级毛片高清牌| 欧美老熟妇乱子伦牲交| 午夜福利影视在线免费观看| 久久久久国内视频| 99国产极品粉嫩在线观看| 中文字幕最新亚洲高清| 国产男女超爽视频在线观看| 亚洲熟女毛片儿| e午夜精品久久久久久久| 国产精品国产av在线观看| 99九九在线精品视频| 一边摸一边抽搐一进一小说 | 黄频高清免费视频| 精品一区二区三区视频在线观看免费 | 国产日韩欧美在线精品| 精品国产乱子伦一区二区三区| 久久中文字幕人妻熟女| 国产无遮挡羞羞视频在线观看| 首页视频小说图片口味搜索| 久久久久精品人妻al黑| 一边摸一边抽搐一进一出视频| 在线观看免费午夜福利视频| 亚洲精品久久成人aⅴ小说| 色综合婷婷激情| 久久精品成人免费网站| 99国产极品粉嫩在线观看| 中文字幕最新亚洲高清| 免费女性裸体啪啪无遮挡网站| 亚洲三区欧美一区| 亚洲精品乱久久久久久| 91麻豆精品激情在线观看国产 | 欧美日韩精品网址| 国产在线免费精品| 法律面前人人平等表现在哪些方面| 80岁老熟妇乱子伦牲交| 法律面前人人平等表现在哪些方面| 涩涩av久久男人的天堂| 一个人免费在线观看的高清视频| 桃花免费在线播放| 午夜福利欧美成人| 国产精品一区二区在线不卡| 亚洲精品美女久久久久99蜜臀| 黄色毛片三级朝国网站| 性高湖久久久久久久久免费观看| 啦啦啦中文免费视频观看日本| 久久久精品免费免费高清| 黑人欧美特级aaaaaa片| 国产在线免费精品| av不卡在线播放| 欧美精品一区二区大全| 一区二区三区激情视频| 精品熟女少妇八av免费久了| 午夜精品国产一区二区电影| 99久久精品国产亚洲精品| 啦啦啦免费观看视频1| 欧美中文综合在线视频| 亚洲精品一二三| 一级毛片精品| 精品久久久久久久毛片微露脸| 日韩免费高清中文字幕av| 中国美女看黄片| 欧美日韩黄片免| 性少妇av在线| 男女无遮挡免费网站观看| www.熟女人妻精品国产| 亚洲精品美女久久av网站| 91九色精品人成在线观看| 51午夜福利影视在线观看| 亚洲精品成人av观看孕妇| 男人操女人黄网站| 蜜桃国产av成人99| 亚洲国产av新网站| 日本欧美视频一区| 欧美国产精品一级二级三级|