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

    中國不同氣候區(qū)參考作物蒸散量計(jì)算模型適用性評價(jià)

    2022-04-06 05:06:34王君勤葉生進(jìn)
    節(jié)水灌溉 2022年3期
    關(guān)鍵詞:氣候區(qū)溫帶適用性

    王君勤,葉生進(jìn),樊 毅

    (1.四川省水利科學(xué)研究院,成都610065;2.四川省都江堰東風(fēng)渠管理處,成都610065)

    0 引 言

    參考作物蒸散量(potential evapotranspiration,ET0)是作物需水量的關(guān)鍵部分,其準(zhǔn)確計(jì)算對精準(zhǔn)灌溉預(yù)報(bào)、灌溉水資源優(yōu)化調(diào)度和區(qū)域水資源規(guī)劃有重要意義。蒸滲儀可實(shí)測ET0,具有較高的精度,但該類方法價(jià)格昂貴且技術(shù)復(fù)雜,實(shí)際應(yīng)用性不強(qiáng)。因此,通常采用經(jīng)驗(yàn)?zāi)P凸浪鉋T0,F(xiàn)AO56 Penman-Monteith(PM)公式因計(jì)算精確且適用范圍廣,一般被作為標(biāo)準(zhǔn)公式。然而PM 模型需要?dú)夂蛞蜃虞^多,嚴(yán)重限制了它的應(yīng)用,同時(shí)許多地區(qū)并沒有標(biāo)準(zhǔn)氣象站提供數(shù)據(jù)支持,因此ET0簡化計(jì)算模型開發(fā)與適用性評價(jià)對于氣象資料缺失地區(qū)的ET0估算尤為重要[1]。ET0簡化計(jì)算模型可分為溫度法、輻射法、質(zhì)量傳遞法和綜合法4種類型,其中溫度法和輻射法由于所需數(shù)據(jù)較少且結(jié)構(gòu)相對簡單在世界范圍內(nèi)廣泛應(yīng)用,然而這兩類模型區(qū)域性強(qiáng)在使用時(shí)必須進(jìn)行適用性評價(jià)。

    Almorox 等[2]在全球研究了11 種基于溫度資料的ET0模型,認(rèn)為Hargreaves 模型精度最高,而McCloud模型和Thornthwaite模型精度最低;Mohammad[3]對22 種輻射模型在伊朗地區(qū)的適用性進(jìn)行分析,得出Stephens 模型在10 個(gè)省份精度較高,回歸分析校正模型后所有省份精度提高,可見輻射模型精度與參數(shù)變化密切;不同類ET0經(jīng)驗(yàn)?zāi)P途C合評價(jià),Azhar等[4]在澳大利亞東南部研究了10個(gè)不同類型的ET0模型,表明任何一個(gè)模型都不能在所有氣候區(qū)有好的表現(xiàn);Tabari 等[5]在伊朗濕潤氣候下分析31 種計(jì)算模型,發(fā)現(xiàn)輻射模型表現(xiàn)較好,基于溫度的模型次之。褚榮浩等[6]評估了13 種典型ET0模型在安徽省的適用性,發(fā)現(xiàn)在氣象數(shù)據(jù)完整時(shí),VA3模型可作為PM 模型的最佳替代公式;在4-10月推薦使用基于溫度和輻射的模型以及VA1和VA2模型,而在其他月份使用質(zhì)量傳輸模型。

    ET0經(jīng)驗(yàn)?zāi)P突径际墙⒃谝欢俣ɑA(chǔ)上的半經(jīng)驗(yàn)半理論公式,具有很強(qiáng)的地區(qū)性,同時(shí)受復(fù)雜的地理、水文及氣候環(huán)境等的影響,目前仍無法解決ET0計(jì)算公式普適性這一科學(xué)難題。現(xiàn)有的ET0計(jì)算模型評價(jià)地區(qū)性均較強(qiáng),且針對中國不同氣候區(qū)基于溫度和輻射的ET0經(jīng)驗(yàn)?zāi)P偷倪m應(yīng)性評價(jià)還未有系統(tǒng)研究。因此,本研究綜合分析了28個(gè)ET0經(jīng)驗(yàn)?zāi)P驮谥袊? 個(gè)氣候區(qū)的適用性,以FAO56 Penman-Monteith 公式為標(biāo)準(zhǔn)模型比較中國不同氣候區(qū)基于溫度和輻射的ET0經(jīng)驗(yàn)?zāi)P偷挠?jì)算精度,評價(jià)不同類型ET0計(jì)算模型在不同氣候區(qū)的適用性;提出不同氣候區(qū)ET0計(jì)算最佳推薦模型,為提高不同氣候區(qū)的ET0估算精度提供支撐。

    1 材料與方法

    1.1 研究區(qū)域和氣象數(shù)據(jù)采集

    中國氣候類型主要包括熱帶季風(fēng)氣候、亞熱帶季風(fēng)氣候、溫帶季風(fēng)氣候、溫帶大陸性氣候和高原山地氣候5 類[7]。Feng等[8]考慮到熱帶季風(fēng)區(qū)面積小、氣象站點(diǎn)少,將該區(qū)域納入到亞熱帶季風(fēng)區(qū)進(jìn)行研究。Yang等[9]根據(jù)柯本氣候分類法將中國大陸地區(qū)劃分為高原山丘區(qū)(HG)、溫帶大陸區(qū)(Dfc)、溫帶季風(fēng)區(qū)(Dwa)和亞熱帶季風(fēng)區(qū)(Cwa)4 個(gè)氣候區(qū)域。因此本研究從溫帶季風(fēng)區(qū)(TMZ)、溫帶大陸區(qū)(TCZ)、高原山丘區(qū)(MPZ)、亞熱帶季風(fēng)區(qū)(SMZ)4 個(gè)區(qū)域分析研究ET0模型的適用性。溫帶季風(fēng)區(qū)夏季高溫多雨,冬季寒冷干燥,屬半濕潤區(qū);溫帶大陸區(qū)夏季高溫,冬季寒冷,全年降水較少,屬于干旱區(qū);高原山丘區(qū),輻射較強(qiáng),但全年低溫,屬半干旱區(qū);亞熱帶季風(fēng)區(qū)夏季高溫多雨,冬季溫和少雨,屬濕潤區(qū)。

    本研究所使用的氣象數(shù)據(jù)來自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://cdc.cma.gov.cn/home.do),由于設(shè)備及天氣等因素的影響,部分站點(diǎn)缺乏輻射及風(fēng)速數(shù)據(jù)或沒有足夠序列的數(shù)據(jù),考慮到數(shù)據(jù)資料的可用性,在全國范圍內(nèi)選取了具有代表性的95 個(gè)站點(diǎn),高原山丘區(qū)、溫帶大陸區(qū)、溫帶季風(fēng)區(qū)和亞熱帶季風(fēng)區(qū)站點(diǎn)數(shù)量分別為11個(gè)、24個(gè)、23個(gè)和37個(gè)。日氣象資料包括日照時(shí)數(shù)(n)、相對濕度(RH)、最高氣溫(Tmax)、最低氣溫(Tmin)、降水量(P)、風(fēng)速(u)和太陽輻射(Rs)。由于實(shí)測輻射數(shù)據(jù)的限制,研究時(shí)段為1994-2016年,其中1994-2010年的數(shù)據(jù)用于模擬校正模型的參數(shù),2011-2016年的數(shù)據(jù)用于驗(yàn)證模型。本文研究區(qū)劃分及氣象站點(diǎn)分布如圖1所示。

    圖1 中國不同氣候區(qū)劃分及氣象站點(diǎn)分布圖Fig.1 Division of different climatic regions and distribution of meteorological stations in China

    1.2 參考作物蒸散量計(jì)算模型

    1.2.1 FAO56 Penman-Monteith模型(PM)

    1998年FAO-56 推薦Penman-Monteith 模型[10]為計(jì)算ET0的標(biāo)準(zhǔn)方法,該方法以空氣動(dòng)力學(xué)和能量平衡原理為基礎(chǔ),理論嚴(yán)謹(jǐn)、精確度較高,在世界范圍內(nèi)廣泛使用。具體公式如下:

    式中:ET0的單位為mm/d;Rn為輸入冠層的凈輻射量,MJ/(m2·d);Ta為日平均氣溫,℃;G為土壤熱通量,MJ/(m2·d);es為飽和水氣壓,kPa;ea為實(shí)際水氣壓,kPa;Δ為飽和水汽壓與溫度曲線的斜率,kPa/℃;γ為干濕溫度計(jì)常數(shù),kPa/℃;u2為距地面2 m高處的風(fēng)速,m/s。

    1.2.2 基于溫度的模型

    溫度測量精度較高,溫度模型對數(shù)據(jù)要求低,應(yīng)用較多[11]。本文選用了13種溫度模型公式,具體見表1。

    表1 參考作物蒸散量溫度模型計(jì)算公式Tab.1 The calculation formulas of temperature-based ET0 method

    1.2.3 基于輻射的模型

    輻射法主要輸入輻射數(shù)據(jù),但是太陽輻射數(shù)據(jù)精度一般都較低,全國僅有123 個(gè)站點(diǎn)實(shí)測輻射數(shù)據(jù)資料可用[25]。本文選用了15種輻射計(jì)算模型,見表2。

    表2 參考作物蒸散量輻射模型計(jì)算公式Tab.2 The calculation formulas of radiation-based ET0 method

    1.2.4 模型評價(jià)指標(biāo)

    本文以PM 模型為標(biāo)準(zhǔn),采用決定系數(shù)(R2)、均方根誤差(RMSE)、相對均方根誤差(RRMSE)、模型有效性系數(shù)(ENS)、平均絕對誤差(MAE)和一致指數(shù)(IOA)6 個(gè)指標(biāo)來評價(jià)模型的表現(xiàn),利用全球績效指標(biāo)(Global Performance Indicator,GPI)綜合分析指標(biāo)效果[38],公式如下:

    式中:Yi為模型計(jì)算的第i日ET0值;Xi為標(biāo)準(zhǔn)模型計(jì)算的第i日ET0值;為Yi的平均值;為Xi的平均值;n為數(shù)據(jù)樣本數(shù)量。R2越接近1,表明模型吻合度越高;RMSE越小,表明模型偏差越?。籖RMSE越小,模型精度越高;ENS越接近1,表明模型效率越高;MAE越接近0,表明誤差越小;IOA越接近1,一致性越好。

    式中:Ti為R2、RMSE、RRMSE、ENS、MAE、IOA歸一化后的值;為對應(yīng)參數(shù)的中位數(shù);αi在Ti為R2、ENS、IOA時(shí)取-1,Ti為RMSE、RRMSE、MAE時(shí)取1。GPI越大即模 型 準(zhǔn)確度越高。

    2 結(jié)果與分析

    2.1 溫度模型適用性評價(jià)

    圖2 為基于溫度法的ET0估算模型在不同氣候區(qū)的排名,溫度模型平均精度表明13 個(gè)溫度模型區(qū)域表現(xiàn)為:TCZ 區(qū)表現(xiàn)最好,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.937、0.56 mm/d、0.446 mm/d、18.73%、0.905、0.978;SMZ區(qū)次之,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.874、 0.502 mm/d、 0.387 mm/d、 18.34%、 0.642、 0.947;MPZ 區(qū)排名第三,R2、RMSE、MAE、RRMSE、ENS和IOA分別 為0.887、 0.517 mm/d、 0.435 mm/d、 19.54%、 0.731、0.948;TMZ 區(qū)表現(xiàn)最差,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.87、0.659 mm/d、0.501 mm/d、25.75 %、0.779、0.951。因此,基于溫度的模型在溫帶大陸區(qū)(TCZ)適用性最好。該類模型所需輸入數(shù)據(jù)少,數(shù)據(jù)精確度高,當(dāng)區(qū)域數(shù)據(jù)資料不足時(shí)可供選擇。

    圖2 基于溫度的ET0模型在不同氣候區(qū)的區(qū)域平均精度Fig.2 Regional average accuracy of the ET0 methods based on temperature in different climatic regions

    圖3 為溫度模型不同氣候區(qū)日ET0比較。除HS3 和Rom 模型外,基于溫度的模型ET0日估算值與PM 模型誤差較小,這可能是由于HS3 模型輸入了日波動(dòng)較大的降雨量值,且本文使用的最小二乘法無法有效優(yōu)化非線性結(jié)構(gòu)Rom 模型的經(jīng)驗(yàn)系數(shù)。4 個(gè)氣候區(qū)PM 模型ET0均值分別為2.558、2.990、2.644、2.736 mm/d,TCZ 區(qū)ET0最大,可見ET0值越大,溫度模型精確度越高。根據(jù)表3溫度模型平均表現(xiàn),TMZ區(qū)B-R模型有最高的GPI 值,R2、RMSE、MAE、RRMSE、ENS和IOA分 別 為0.926、0.425 mm/d、0.337 mm/d、16.62%、0.923、0.980(ET0均值為2.482 mm/d,圖3)。HScor 模型表現(xiàn)與B-R模型相近,排名第二(ET0均值為2.517 mm/d);TCZ區(qū)M-B模型表現(xiàn)最佳,GPI 為0.181,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.968、0.354 mm/d、0.249 mm/d、11.82%、0.968、0.992(ET0均值為2.956 mm/d);MPZ 區(qū)B-R 模型可得到更精確的ET0,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.952、0.271 mm/d、0.208 mm/d、10.36%、0.95、0.987(ET0均值為2.581 mm/d);SMZ 區(qū)HScor 模型排名第1,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.945、0.255 mm/d、0.195 mm/d、9.32%、0.944、0.985(ET0均值為2.702 mm/d)??傊琓MZ 和MPZ 區(qū)B-R 模型適用性最 好,TCZ 和SMZ 區(qū)分別推薦MB、HScor模型,4個(gè)區(qū)域最適宜溫度模型輸入氣象因子為Tmax和Tmin,與其他溫度模型相比包含了大氣頂層輻射Ra,可見Ra對于ET0估算具有積極影響。

    表3 基于溫度法的ET0計(jì)算模型在不同氣候區(qū)計(jì)算精度Tab.3 The calculation accuracy of ET0 methods based on temperature in different climate zones

    圖3 溫度模型不同氣候區(qū)日ET0比較Fig.3 Comparison of daily ET0 values of temperature methods in different climatic zones

    2.2 輻射模型適用性評價(jià)

    圖4 為基于輻射法的ET0估算模型在不同氣候區(qū)的表現(xiàn)排名,輻射模型平均精度表明15 個(gè)輻射模型區(qū)域表現(xiàn)為:SMZ區(qū)模型精度最高,R2、RMSE、MAE、RRMSE、ENS和IOA分別 為0.965、 0.269 mm/d、 0.222 mm/d、 9.844%、 0.923、0.981;TCZ 區(qū)次之,R2、RMSE、MAE、RRMSE、ENS和IOA分 別 為0.964、0.515 mm/d、0.437 mm/d、17.24%、0.903、0.976;TMZ 區(qū)排名第三,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.944、0.464 mm/d、0.389 mm/d、18.14%、0.872、0.97;MPZ 區(qū)表現(xiàn)最差,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.946、0.49 mm/d、0.441 mm/d、18.53%、0.645、0.952。濕潤地區(qū)水分條件充足,空氣動(dòng)力學(xué)的影響也較小,可用能量輻射值的多少是影響ET0的關(guān)鍵,因而以輻射作為主要輸入的輻射模型在亞熱帶季風(fēng)區(qū)(SMZ)表現(xiàn)最好。

    圖5 為輻射模型不同氣候區(qū)日ET0比較。除Turc 模型,基于輻射的模型ET0日估算值與PM 模型誤差較小。Turc 模型表現(xiàn)較差的原因可能是該模型為非線性結(jié)構(gòu),最小二乘法無法高效地提高其精度。由表4 可知,TMZ 區(qū)和TCZ 區(qū)FAO24Rad模型GPI值最高,R2、RMSE、MAE、RRMSE、ENS和IOA分別 為0.982、0.99,0.245、0.237 mm/d,0.21、0.199 mm/d,9.566%、7.913%,0.975、0.986,0.993、0.996 (ET0均值為2.611、3.008 mm/d,圖4);MPZ 和SMZ 區(qū)PT 模型表現(xiàn)最佳,GPI 分別為0.118、0.224,R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.984、0.979,0.179、0.159 mm/d,0.133、0.125 mm/d,6.772%、5.802%,0.978、0.978,0.995、0.995(ET0均值為2.553、2.733 mm/d)??偠灾現(xiàn)AO24Rad 模型在北方地區(qū)適用性較好,PT模型在南方地區(qū)有突出表現(xiàn)。FAO24Rad模型是當(dāng)風(fēng)速和相對濕度無法獲得時(shí)提出的一種PM 模型的替代公式,劉曉英等[39]也表明FAO24Rad 模型在半干旱區(qū)具有較高的ET0精度。PT 模型在澳大利亞利用濕潤路面的資料提出,以平衡蒸發(fā)為基礎(chǔ),引進(jìn)經(jīng)驗(yàn)常數(shù)導(dǎo)出的估算無平流條件蒸發(fā)的模型,因而其在南方地區(qū)表現(xiàn)更好[40]。

    圖5 輻射模型不同氣候區(qū)ET0比較Fig.5 Comparison of daily ET0 values of radiation methods in different climate zones

    2.3 不同氣候區(qū)模型推薦

    圖6通過GPI值及排名對28個(gè)模型表現(xiàn)做綜合評價(jià),輻射類模型總體表現(xiàn)比溫度模型好,GPI 排名前十的輻射模型4 個(gè)氣候區(qū)分別有10、8、9、10個(gè)。TMZ區(qū)ET0估算首選輻射模型的FAO24Rad 模型,最優(yōu)的溫度模型為B-R 模型,排名為12;TCZ 區(qū)ET0推薦模型分別為輻射模型FAO24Rad(第1 名)、溫度模型MB(2005)(第7 名);MPZ 區(qū)最優(yōu)推薦為輻射模型PT,此外B-R 是表現(xiàn)最佳的溫度模型(排名9);SMZ 區(qū),輻射法中PT 模型排名第1,表現(xiàn)優(yōu)于溫度模型,溫度法推薦HScor(排名12)。綜上所述,4 個(gè)氣候區(qū)表現(xiàn)最優(yōu)的為輻射模型FAO24Rad(TMZ、TCZ 區(qū))和PT(MPZ、SMZ)。然而輻射模型雖表現(xiàn)較好,但輻射數(shù)據(jù)有時(shí)很難獲得,特別是在一些發(fā)展中國家。當(dāng)輻射資料無法獲取時(shí),表現(xiàn)較好的溫度模型是一個(gè)不錯(cuò)的選擇,TMZ 和MPZ 區(qū)可選用B-R 模型,TCZ 區(qū)選用MB(2005)模型,且SMZ 區(qū)選用HScor 模型。上述結(jié)果表明溫度模型相較于輻射模型表現(xiàn)較差,這可能是由于該類模型輸入的溫度參數(shù)比較容易受云層厚度、天氣和風(fēng)速變化的影響,在日尺度上波動(dòng)較大,產(chǎn)生了更多的誤差[41]。此外本文以實(shí)測輻射數(shù)據(jù)為輸入變量,這可能是輻射模型表現(xiàn)較好的另一原因。

    圖6 經(jīng)驗(yàn)?zāi)P突赑M模型在不同氣候區(qū)表現(xiàn)的GPI值及排名Fig.6 GPI values and ranking of the empirical methods Based on the PM method in different climate zones

    3 結(jié) 論

    本文以FAO56 Penman-Monteith 公式為標(biāo)準(zhǔn)模型,研究28個(gè)ET0經(jīng)驗(yàn)?zāi)P停?3 個(gè)溫度模型和15 個(gè)輻射模型)在溫帶季風(fēng)區(qū)(TMZ)、溫帶大陸區(qū)(TCZ)、高原山丘區(qū)(MPZ)和亞熱帶季風(fēng)區(qū)(SMZ)4 個(gè)氣候區(qū)的適用性。通過統(tǒng)計(jì)指標(biāo)R2、RMSE、MAE、RRMSE、ENS、IOA及GPI的分析,得出以下結(jié)論:

    (1)ET0經(jīng)驗(yàn)?zāi)P驮? 個(gè)氣候區(qū)的適用性差異較大。溫度模型在TCZ 區(qū)有最好表現(xiàn)(R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.937、0.56 mm/d、0.446 mm/d、18.73 %、0.905、0.978);輻射模型在SMZ 區(qū)模型精度最高(R2、RMSE、MAE、RRMSE、ENS和IOA分別為0.965、0.269 mm/d、0.222 mm/d、9.844%、0.923、0.981)。

    (2)全國范圍內(nèi)基于輻射的ET0模型表現(xiàn)優(yōu)于溫度模型。不同氣候區(qū)ET0估算首先推薦輻射模型,TMZ 和TCZ 區(qū)為FAO24Rad 模型,MPZ 和SMZ 區(qū)PT 模型與FAO56 PM 的ET0值最接近。當(dāng)輻射數(shù)據(jù)不可獲得時(shí),4 個(gè)氣候區(qū)(TMZ、TCZ、MPZ 和SMZ)分別推薦基于溫度的B-R 模型,MB(2005)模型、B-R模型、HScor模型。

    猜你喜歡
    氣候區(qū)溫帶適用性
    第31期 參考答案
    Facts of Yellowstone
    強(qiáng)調(diào)簡潔和適用性 MICHI by Rotel X5/X3合并功放
    環(huán)保技術(shù)在土木工程領(lǐng)域的適用性探討
    小議階級分析理論的適用性
    各氣候區(qū)被動(dòng)建筑節(jié)能技術(shù)實(shí)用性分析
    科技資訊(2017年23期)2017-09-09 13:08:43
    冷涼氣候區(qū)不同形態(tài)氮肥對蘋果根系氮代謝的影響
    濕熱氣候區(qū)建筑防熱研究進(jìn)展
    新一代車上轉(zhuǎn)換裝置在廠礦企業(yè)鐵路運(yùn)輸中的適用性研究
    中國不同氣候區(qū)土壤濕度特征及其氣候響應(yīng)
    超碰成人久久| 午夜两性在线视频| 亚洲精品在线观看二区| 长腿黑丝高跟| av福利片在线| 黄色a级毛片大全视频| 国产精品1区2区在线观看.| 亚洲自拍偷在线| videosex国产| 操美女的视频在线观看| 9191精品国产免费久久| 久久久久亚洲av毛片大全| 长腿黑丝高跟| bbb黄色大片| av国产精品久久久久影院| 水蜜桃什么品种好| 亚洲视频免费观看视频| 久久性视频一级片| 久久影院123| 精品福利观看| 一a级毛片在线观看| 日韩欧美免费精品| 国产色视频综合| 中出人妻视频一区二区| 欧美人与性动交α欧美软件| 99re在线观看精品视频| 一级片免费观看大全| 99久久综合精品五月天人人| 国产熟女xx| 99精品久久久久人妻精品| 日本 av在线| 男女午夜视频在线观看| 国产男靠女视频免费网站| 色综合站精品国产| 十分钟在线观看高清视频www| 窝窝影院91人妻| 身体一侧抽搐| 亚洲av日韩精品久久久久久密| 国产熟女午夜一区二区三区| 中文字幕高清在线视频| 99国产极品粉嫩在线观看| av有码第一页| 99久久精品国产亚洲精品| 免费一级毛片在线播放高清视频 | 50天的宝宝边吃奶边哭怎么回事| 999精品在线视频| 亚洲视频免费观看视频| 国产熟女午夜一区二区三区| 亚洲欧美一区二区三区久久| 亚洲国产毛片av蜜桃av| 久久中文字幕人妻熟女| 夜夜夜夜夜久久久久| 咕卡用的链子| 免费看a级黄色片| 久久久久久大精品| 精品无人区乱码1区二区| 看黄色毛片网站| 97碰自拍视频| 欧美日韩乱码在线| 国内久久婷婷六月综合欲色啪| 黄色成人免费大全| 亚洲成人国产一区在线观看| 午夜福利在线观看吧| 高清在线国产一区| 每晚都被弄得嗷嗷叫到高潮| 天天添夜夜摸| av网站在线播放免费| 色老头精品视频在线观看| 亚洲国产看品久久| 大型av网站在线播放| 美女扒开内裤让男人捅视频| 中文亚洲av片在线观看爽| 亚洲精品一二三| 夜夜爽天天搞| 极品人妻少妇av视频| 亚洲精华国产精华精| 91成人精品电影| 久久国产乱子伦精品免费另类| 黑人巨大精品欧美一区二区mp4| 久久婷婷成人综合色麻豆| 天天躁夜夜躁狠狠躁躁| 国产一区在线观看成人免费| 日韩欧美一区视频在线观看| 日韩三级视频一区二区三区| 18禁美女被吸乳视频| 欧美精品一区二区免费开放| 久久精品国产99精品国产亚洲性色 | 制服诱惑二区| 青草久久国产| 久久久久国产一级毛片高清牌| 久久 成人 亚洲| 一级黄色大片毛片| 国产高清视频在线播放一区| 大码成人一级视频| www.熟女人妻精品国产| 在线观看免费高清a一片| 一级a爱视频在线免费观看| 国产精品综合久久久久久久免费 | 一区二区三区精品91| 中文字幕人妻熟女乱码| 国产激情欧美一区二区| 亚洲精品一区av在线观看| 男人操女人黄网站| 香蕉久久夜色| 精品国产美女av久久久久小说| 国产成人影院久久av| 男女高潮啪啪啪动态图| 久久婷婷成人综合色麻豆| av有码第一页| 亚洲国产中文字幕在线视频| 久久精品91蜜桃| 91精品三级在线观看| 三级毛片av免费| av电影中文网址| 啦啦啦免费观看视频1| 满18在线观看网站| √禁漫天堂资源中文www| 97人妻天天添夜夜摸| 亚洲视频免费观看视频| 老司机午夜十八禁免费视频| 国产精品98久久久久久宅男小说| 国产黄a三级三级三级人| bbb黄色大片| 欧美乱色亚洲激情| 欧美黑人欧美精品刺激| 国产成人av教育| 国产成+人综合+亚洲专区| av欧美777| 嫁个100分男人电影在线观看| 日韩大尺度精品在线看网址 | 国产免费现黄频在线看| e午夜精品久久久久久久| 久久人人爽av亚洲精品天堂| 在线观看日韩欧美| 久久中文字幕人妻熟女| 最近最新中文字幕大全免费视频| av福利片在线| 日韩精品中文字幕看吧| 精品乱码久久久久久99久播| 久久久国产成人精品二区 | 一个人免费在线观看的高清视频| 色尼玛亚洲综合影院| 热re99久久精品国产66热6| 五月开心婷婷网| 夜夜夜夜夜久久久久| 欧美日韩精品网址| 亚洲av日韩精品久久久久久密| www国产在线视频色| 水蜜桃什么品种好| 亚洲狠狠婷婷综合久久图片| 热re99久久精品国产66热6| 五月开心婷婷网| 精品久久蜜臀av无| 又紧又爽又黄一区二区| 国产精品自产拍在线观看55亚洲| 久久久国产一区二区| 一区二区三区激情视频| 国产精品av久久久久免费| 欧美另类亚洲清纯唯美| 亚洲精品国产区一区二| 人成视频在线观看免费观看| 色综合婷婷激情| 国产一区在线观看成人免费| 日韩有码中文字幕| 欧美 亚洲 国产 日韩一| 80岁老熟妇乱子伦牲交| 成人黄色视频免费在线看| 丝袜美足系列| 欧美在线黄色| 欧美色视频一区免费| 午夜91福利影院| 午夜a级毛片| 欧美成狂野欧美在线观看| 亚洲av片天天在线观看| 男女下面进入的视频免费午夜 | 91麻豆精品激情在线观看国产 | 在线观看免费视频网站a站| 欧美日韩视频精品一区| 热99re8久久精品国产| netflix在线观看网站| 亚洲国产看品久久| 国产国语露脸激情在线看| 热99国产精品久久久久久7| 18禁国产床啪视频网站| 18禁裸乳无遮挡免费网站照片 | 免费在线观看影片大全网站| 亚洲成人免费av在线播放| 亚洲男人天堂网一区| 极品教师在线免费播放| 久久香蕉激情| 日韩欧美在线二视频| 一本大道久久a久久精品| 最好的美女福利视频网| 亚洲av熟女| 最新美女视频免费是黄的| 悠悠久久av| 亚洲av成人av| tocl精华| 窝窝影院91人妻| 夜夜躁狠狠躁天天躁| 在线播放国产精品三级| 久久精品亚洲精品国产色婷小说| 天天躁狠狠躁夜夜躁狠狠躁| 在线观看一区二区三区| 久久久久久久久久久久大奶| 久久久精品国产亚洲av高清涩受| 少妇的丰满在线观看| 黄色毛片三级朝国网站| 一级a爱视频在线免费观看| 亚洲精品在线美女| 日本黄色日本黄色录像| 久9热在线精品视频| 免费看a级黄色片| 99精品欧美一区二区三区四区| 欧美在线一区亚洲| 这个男人来自地球电影免费观看| 后天国语完整版免费观看| 十分钟在线观看高清视频www| 中文字幕av电影在线播放| 久久国产精品影院| 欧美老熟妇乱子伦牲交| 少妇裸体淫交视频免费看高清 | 999久久久精品免费观看国产| 国产成年人精品一区二区 | 丁香六月欧美| 日韩精品青青久久久久久| 中文字幕人妻丝袜一区二区| 午夜福利在线观看吧| 91大片在线观看| 国产精品亚洲一级av第二区| 国产精品免费一区二区三区在线| 亚洲国产精品一区二区三区在线| 亚洲人成伊人成综合网2020| 欧美最黄视频在线播放免费 | 天天躁狠狠躁夜夜躁狠狠躁| 高清黄色对白视频在线免费看| 国产熟女午夜一区二区三区| 99国产精品一区二区蜜桃av| 欧美丝袜亚洲另类 | 国产又色又爽无遮挡免费看| 日韩三级视频一区二区三区| 免费高清在线观看日韩| 国产单亲对白刺激| 9热在线视频观看99| 80岁老熟妇乱子伦牲交| 成年人黄色毛片网站| 成人精品一区二区免费| 最好的美女福利视频网| 亚洲狠狠婷婷综合久久图片| 性少妇av在线| 久久久国产成人精品二区 | 色在线成人网| 国产成年人精品一区二区 | 精品一区二区三区av网在线观看| 欧美成人午夜精品| 欧美激情 高清一区二区三区| 久久久国产成人免费| 亚洲欧美精品综合久久99| 啦啦啦在线免费观看视频4| 亚洲欧美激情在线| 婷婷六月久久综合丁香| 欧美在线一区亚洲| 免费在线观看亚洲国产| 欧美大码av| 丝袜在线中文字幕| 日韩视频一区二区在线观看| 高清毛片免费观看视频网站 | 日日爽夜夜爽网站| 老汉色∧v一级毛片| 两个人免费观看高清视频| 亚洲一码二码三码区别大吗| 中文欧美无线码| 天堂中文最新版在线下载| 亚洲 欧美 日韩 在线 免费| 婷婷六月久久综合丁香| 欧美成人性av电影在线观看| 一a级毛片在线观看| 成在线人永久免费视频| 十八禁人妻一区二区| 久久伊人香网站| 女同久久另类99精品国产91| av欧美777| 97热精品久久久久久| 九九热线精品视视频播放| 午夜激情福利司机影院| 成人高潮视频无遮挡免费网站| 热99在线观看视频| 国产淫片久久久久久久久 | 毛片女人毛片| 中文字幕高清在线视频| 精品一区二区三区视频在线观看免费| 亚洲av.av天堂| 亚洲,欧美精品.| 一区福利在线观看| 欧美黑人巨大hd| 成人鲁丝片一二三区免费| 伦理电影大哥的女人| 一区福利在线观看| 国模一区二区三区四区视频| 又爽又黄a免费视频| 国产精品综合久久久久久久免费| 欧美性猛交╳xxx乱大交人| 欧美国产日韩亚洲一区| 国产高潮美女av| 免费看美女性在线毛片视频| 9191精品国产免费久久| 舔av片在线| av专区在线播放| 亚洲自偷自拍三级| 亚洲在线自拍视频| 日韩欧美三级三区| 亚洲第一电影网av| 欧美午夜高清在线| 国产成人欧美在线观看| 亚洲成人免费电影在线观看| 性色av乱码一区二区三区2| 色av中文字幕| 久久热精品热| 国产麻豆成人av免费视频| 成年女人看的毛片在线观看| 日韩有码中文字幕| 免费在线观看影片大全网站| 久久久色成人| 国产精品久久久久久久电影| 性色avwww在线观看| 老熟妇乱子伦视频在线观看| 久久亚洲真实| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩中文字幕国产精品一区二区三区| 一级黄片播放器| 91麻豆精品激情在线观看国产| 亚洲熟妇熟女久久| 国产亚洲欧美98| 网址你懂的国产日韩在线| 久久精品国产亚洲av香蕉五月| 日韩欧美精品免费久久 | 国产午夜精品久久久久久一区二区三区 | 久99久视频精品免费| 国产成+人综合+亚洲专区| 婷婷六月久久综合丁香| 亚洲国产色片| 脱女人内裤的视频| 欧美一级a爱片免费观看看| 免费高清视频大片| 99久久成人亚洲精品观看| 国产男靠女视频免费网站| 亚洲国产色片| 好男人电影高清在线观看| 午夜福利视频1000在线观看| 国产视频一区二区在线看| 一区二区三区免费毛片| 在线天堂最新版资源| 亚洲精品影视一区二区三区av| 成人高潮视频无遮挡免费网站| 观看美女的网站| 欧美+日韩+精品| 国产精品女同一区二区软件 | 亚洲国产高清在线一区二区三| a在线观看视频网站| 日本成人三级电影网站| 亚洲精品久久国产高清桃花| 亚洲中文字幕日韩| 真人一进一出gif抽搐免费| 国产一区二区三区视频了| 午夜免费激情av| 一级黄色大片毛片| a在线观看视频网站| 亚洲成人免费电影在线观看| 亚洲欧美日韩高清专用| 中文字幕av成人在线电影| 精品欧美国产一区二区三| 亚洲一区二区三区色噜噜| 男人的好看免费观看在线视频| 一个人看视频在线观看www免费| 18禁在线播放成人免费| 人妻制服诱惑在线中文字幕| 日韩av在线大香蕉| 亚洲av二区三区四区| 一个人观看的视频www高清免费观看| 欧美激情久久久久久爽电影| 亚洲精华国产精华精| 在线观看66精品国产| 18美女黄网站色大片免费观看| 精品福利观看| 嫁个100分男人电影在线观看| 国产私拍福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 丰满的人妻完整版| 91字幕亚洲| 尤物成人国产欧美一区二区三区| 免费av毛片视频| 国产aⅴ精品一区二区三区波| 国产私拍福利视频在线观看| 中文字幕久久专区| 国产精品爽爽va在线观看网站| 99久久99久久久精品蜜桃| 99国产综合亚洲精品| 特级一级黄色大片| 免费人成视频x8x8入口观看| 色吧在线观看| 日日夜夜操网爽| 成人特级av手机在线观看| 午夜日韩欧美国产| 亚洲av一区综合| 久久精品国产亚洲av涩爱 | 国产不卡一卡二| 中出人妻视频一区二区| 亚洲色图av天堂| 一级黄色大片毛片| 午夜影院日韩av| 午夜福利在线观看吧| 在线十欧美十亚洲十日本专区| 久久婷婷人人爽人人干人人爱| 欧美丝袜亚洲另类 | 人人妻人人看人人澡| 久久久成人免费电影| 成人鲁丝片一二三区免费| 一区二区三区高清视频在线| 夜夜爽天天搞| 特级一级黄色大片| or卡值多少钱| 国产伦在线观看视频一区| 一个人免费在线观看电影| 国产精品乱码一区二三区的特点| 成人无遮挡网站| 精品一区二区三区视频在线观看免费| 18禁黄网站禁片午夜丰满| 欧美日本亚洲视频在线播放| 精品99又大又爽又粗少妇毛片 | 毛片一级片免费看久久久久 | 亚洲av成人不卡在线观看播放网| 美女高潮的动态| 男人舔女人下体高潮全视频| 又黄又爽又刺激的免费视频.| 免费高清视频大片| 久久久久久久久中文| 欧美日韩乱码在线| 一级黄色大片毛片| 成人av在线播放网站| av中文乱码字幕在线| 免费看日本二区| 一区福利在线观看| 天堂网av新在线| 一个人免费在线观看的高清视频| 久久久久国产精品人妻aⅴ院| 欧美成人a在线观看| 动漫黄色视频在线观看| 亚洲人成伊人成综合网2020| 午夜a级毛片| 久久国产精品人妻蜜桃| 淫秽高清视频在线观看| 在线十欧美十亚洲十日本专区| 欧美成人一区二区免费高清观看| 一进一出抽搐gif免费好疼| 一级黄片播放器| 3wmmmm亚洲av在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲午夜理论影院| 搡老妇女老女人老熟妇| 熟女电影av网| 少妇人妻精品综合一区二区 | 亚洲欧美激情综合另类| 精品久久久久久久久亚洲 | 成年女人看的毛片在线观看| 亚洲人成伊人成综合网2020| av欧美777| 精品一区二区三区av网在线观看| 听说在线观看完整版免费高清| 国产精品国产高清国产av| 男插女下体视频免费在线播放| 中文字幕av成人在线电影| 欧美3d第一页| 香蕉av资源在线| 99久久精品一区二区三区| 一个人免费在线观看电影| 麻豆成人av在线观看| 成人av在线播放网站| 精品久久久久久,| 久久久久精品国产欧美久久久| 亚洲国产欧洲综合997久久,| 国产精品亚洲美女久久久| 18美女黄网站色大片免费观看| 日韩欧美 国产精品| 国产单亲对白刺激| 日韩欧美三级三区| 欧美日本视频| 性插视频无遮挡在线免费观看| 99久久久亚洲精品蜜臀av| 欧美中文日本在线观看视频| 久久精品国产清高在天天线| 国产精品人妻久久久久久| 日韩av在线大香蕉| 12—13女人毛片做爰片一| 最新在线观看一区二区三区| 国产私拍福利视频在线观看| netflix在线观看网站| 天天躁日日操中文字幕| 久久热精品热| 国产精品久久电影中文字幕| 亚洲av成人av| 免费人成在线观看视频色| 亚洲内射少妇av| 美女被艹到高潮喷水动态| 天堂√8在线中文| 99热精品在线国产| 波多野结衣高清作品| 最新在线观看一区二区三区| 色视频www国产| 日本五十路高清| 听说在线观看完整版免费高清| 国产成+人综合+亚洲专区| 99riav亚洲国产免费| 国产乱人视频| 精品久久久久久成人av| 女生性感内裤真人,穿戴方法视频| 亚洲最大成人中文| 国产精品嫩草影院av在线观看 | 久久精品夜夜夜夜夜久久蜜豆| 精品国内亚洲2022精品成人| av视频在线观看入口| 色哟哟哟哟哟哟| 男人狂女人下面高潮的视频| 精品久久国产蜜桃| 婷婷亚洲欧美| 亚洲成人久久爱视频| 最新中文字幕久久久久| 啦啦啦韩国在线观看视频| 国产高清有码在线观看视频| 日韩欧美 国产精品| 人妻久久中文字幕网| 国产av不卡久久| 亚洲国产精品合色在线| 好男人在线观看高清免费视频| 美女大奶头视频| 男女那种视频在线观看| 免费大片18禁| 久久国产精品影院| 日本撒尿小便嘘嘘汇集6| 高清日韩中文字幕在线| 国产伦在线观看视频一区| 午夜福利视频1000在线观看| 深夜a级毛片| 天堂网av新在线| 91在线精品国自产拍蜜月| av在线蜜桃| 亚洲av电影不卡..在线观看| 亚洲国产欧美人成| 欧美最新免费一区二区三区 | 麻豆久久精品国产亚洲av| 极品教师在线视频| 国产国拍精品亚洲av在线观看| 伊人久久精品亚洲午夜| 一级作爱视频免费观看| 日本 av在线| 午夜日韩欧美国产| 天天躁日日操中文字幕| 久久精品91蜜桃| 欧美精品国产亚洲| 午夜精品久久久久久毛片777| 人人妻人人澡欧美一区二区| 中出人妻视频一区二区| 国产一区二区在线观看日韩| 最近最新中文字幕大全电影3| 此物有八面人人有两片| 首页视频小说图片口味搜索| 久久亚洲精品不卡| 国产在线精品亚洲第一网站| 国产一区二区激情短视频| 真人做人爱边吃奶动态| 长腿黑丝高跟| 亚洲av中文字字幕乱码综合| 色哟哟·www| 欧美日本视频| 欧美zozozo另类| 欧美日韩综合久久久久久 | 免费电影在线观看免费观看| 窝窝影院91人妻| 舔av片在线| 日韩国内少妇激情av| 女人十人毛片免费观看3o分钟| 女同久久另类99精品国产91| 青草久久国产| 中文字幕av成人在线电影| 极品教师在线视频| 91麻豆精品激情在线观看国产| 日韩欧美在线乱码| 久久午夜福利片| 嫩草影视91久久| 国产久久久一区二区三区| АⅤ资源中文在线天堂| 国产探花极品一区二区| 男人的好看免费观看在线视频| 男人舔奶头视频| 别揉我奶头 嗯啊视频| 麻豆国产97在线/欧美| 欧美激情久久久久久爽电影| 久9热在线精品视频| 深夜a级毛片| 国产色爽女视频免费观看| 色哟哟·www| 日韩精品中文字幕看吧| 麻豆国产av国片精品| 亚洲精品在线观看二区| 久久精品国产亚洲av天美| АⅤ资源中文在线天堂| 国产午夜精品论理片| 嫩草影院入口| 亚洲专区中文字幕在线| 久久久国产成人精品二区| 精品国产三级普通话版| 国产色婷婷99| 亚洲av第一区精品v没综合| 成人特级黄色片久久久久久久| 美女 人体艺术 gogo| 日本 av在线|