• <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电影在线播放| 中文字幕人妻熟女乱码| 欧美精品一区二区大全| 国产免费福利视频在线观看| 亚洲人成77777在线视频| 夫妻性生交免费视频一级片| 岛国毛片在线播放| 尾随美女入室| 老汉色av国产亚洲站长工具| 国产伦理片在线播放av一区| 免费看不卡的av| 久久人人97超碰香蕉20202| 久久精品久久久久久久性| 国产精品国产三级国产专区5o| 在现免费观看毛片| 高清不卡的av网站| 在线观看美女被高潮喷水网站| 制服人妻中文乱码| 啦啦啦中文免费视频观看日本| 国产精品一区二区在线不卡| 大香蕉久久成人网| 国产极品天堂在线| 最新的欧美精品一区二区| 亚洲一区二区三区欧美精品| 成人毛片a级毛片在线播放| 永久免费av网站大全| 亚洲国产精品一区二区三区在线| 成人毛片60女人毛片免费| 男女下面插进去视频免费观看| 大话2 男鬼变身卡| 国产成人精品久久久久久| 成年美女黄网站色视频大全免费| 国产精品.久久久| videosex国产| 制服丝袜香蕉在线| 国产女主播在线喷水免费视频网站| 欧美xxⅹ黑人| 熟女电影av网| 亚洲av电影在线进入| 捣出白浆h1v1| 在线精品无人区一区二区三| 黄色配什么色好看| 大话2 男鬼变身卡| 国产一区二区激情短视频 | 欧美97在线视频| 日韩电影二区| 秋霞伦理黄片| 亚洲视频免费观看视频| 免费在线观看完整版高清| 国产午夜精品一二区理论片| 亚洲三级黄色毛片| 亚洲精品久久成人aⅴ小说| 少妇被粗大的猛进出69影院| 一区二区日韩欧美中文字幕| 秋霞伦理黄片| 欧美激情 高清一区二区三区| 欧美+日韩+精品| 女人久久www免费人成看片| 久久女婷五月综合色啪小说| 你懂的网址亚洲精品在线观看| 亚洲欧美一区二区三区黑人 | 婷婷色综合www| 欧美日韩av久久| 老司机亚洲免费影院| 不卡av一区二区三区| 亚洲精品日本国产第一区| 亚洲欧美一区二区三区国产| 狂野欧美激情性bbbbbb| 国产精品久久久久久久久免| 少妇精品久久久久久久| 欧美变态另类bdsm刘玥| kizo精华| 欧美激情 高清一区二区三区| 宅男免费午夜| 岛国毛片在线播放| 精品国产一区二区三区久久久樱花| 91精品三级在线观看| 美女午夜性视频免费| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 狠狠婷婷综合久久久久久88av| 久久99一区二区三区| 亚洲精品视频女| 亚洲av国产av综合av卡| 国产极品粉嫩免费观看在线| 少妇的丰满在线观看| 欧美中文综合在线视频| 欧美另类一区| 91久久精品国产一区二区三区| 多毛熟女@视频| 青青草视频在线视频观看| 好男人视频免费观看在线| 不卡av一区二区三区| 久久人人97超碰香蕉20202| 精品国产一区二区三区四区第35| 不卡视频在线观看欧美| 久久久久久久久久久久大奶| 制服诱惑二区| 久久精品国产自在天天线| 国产97色在线日韩免费| 成人国产av品久久久| 在线观看www视频免费| 中文字幕亚洲精品专区| 丝袜喷水一区| 国产一区二区三区综合在线观看| 久久久久久久久久久免费av| 欧美日韩一区二区视频在线观看视频在线| 亚洲激情五月婷婷啪啪| 99香蕉大伊视频| 亚洲精品乱久久久久久| 日韩av免费高清视频| 亚洲,一卡二卡三卡| 婷婷色综合www| 热re99久久精品国产66热6| 国产野战对白在线观看| 久久精品aⅴ一区二区三区四区 | 五月伊人婷婷丁香| 亚洲综合色网址| 在线观看美女被高潮喷水网站| 精品久久久精品久久久| 少妇 在线观看| 午夜福利视频在线观看免费| 日韩在线高清观看一区二区三区| 国产av码专区亚洲av| 18在线观看网站| 一本色道久久久久久精品综合| 久久综合国产亚洲精品| 丰满少妇做爰视频| 欧美精品亚洲一区二区| 久久人人97超碰香蕉20202| 秋霞伦理黄片| www.av在线官网国产| 国产成人精品久久久久久| 最新的欧美精品一区二区| 少妇被粗大猛烈的视频| 国产精品人妻久久久影院| 亚洲欧洲日产国产| 亚洲欧美一区二区三区久久| 最近的中文字幕免费完整| av在线app专区| 91精品三级在线观看| 亚洲欧美精品自产自拍| 中文字幕人妻熟女乱码| 中国国产av一级| 多毛熟女@视频| 久久午夜综合久久蜜桃| 色播在线永久视频| 亚洲精品中文字幕在线视频| 国产免费又黄又爽又色| 国产免费现黄频在线看| 青草久久国产| 大片电影免费在线观看免费| 99久久综合免费| 少妇 在线观看| 亚洲av在线观看美女高潮| 亚洲视频免费观看视频| 五月开心婷婷网| 视频在线观看一区二区三区| 精品一区二区三区四区五区乱码 | 十八禁高潮呻吟视频| 久久鲁丝午夜福利片| 午夜福利网站1000一区二区三区| 久久久精品国产亚洲av高清涩受| 校园人妻丝袜中文字幕| 男人爽女人下面视频在线观看| 国产爽快片一区二区三区| 1024香蕉在线观看| 亚洲精品视频女| 日韩三级伦理在线观看| 国产1区2区3区精品| 侵犯人妻中文字幕一二三四区| av在线老鸭窝| 免费大片黄手机在线观看| 日韩人妻精品一区2区三区| 十分钟在线观看高清视频www| 久久青草综合色| 黄色视频在线播放观看不卡| 午夜影院在线不卡| 国产精品.久久久| 青草久久国产| www.精华液| 王馨瑶露胸无遮挡在线观看| 性色av一级| 99国产综合亚洲精品| 国产成人精品久久二区二区91 | 日韩在线高清观看一区二区三区| 不卡视频在线观看欧美| 久久精品亚洲av国产电影网| 男女边吃奶边做爰视频| 夜夜骑夜夜射夜夜干| 国产成人欧美| 热99久久久久精品小说推荐| 午夜福利在线观看免费完整高清在| 如日韩欧美国产精品一区二区三区| 夫妻午夜视频| 欧美国产精品va在线观看不卡| 韩国精品一区二区三区| 精品人妻熟女毛片av久久网站| 人成视频在线观看免费观看| 色婷婷av一区二区三区视频| 中国国产av一级| 女性被躁到高潮视频| 高清在线视频一区二区三区| 亚洲在久久综合| 亚洲精品日韩在线中文字幕| 日本猛色少妇xxxxx猛交久久| 亚洲欧美清纯卡通| 午夜精品国产一区二区电影| 另类亚洲欧美激情| 久久99一区二区三区| 两性夫妻黄色片| 国产精品 欧美亚洲| 高清欧美精品videossex| 亚洲图色成人| 人人妻人人澡人人爽人人夜夜| 69精品国产乱码久久久| 男女高潮啪啪啪动态图| 一级毛片我不卡| 日本91视频免费播放| videosex国产| 欧美bdsm另类| 卡戴珊不雅视频在线播放| 国产欧美日韩一区二区三区在线| 伦精品一区二区三区| 一区在线观看完整版| 精品一区在线观看国产| 最近2019中文字幕mv第一页| 久久这里有精品视频免费| 天堂中文最新版在线下载| 国产精品亚洲av一区麻豆 | 午夜91福利影院| 国产深夜福利视频在线观看| 女人久久www免费人成看片| 国产极品粉嫩免费观看在线| 久久 成人 亚洲| 亚洲av在线观看美女高潮| 美女高潮到喷水免费观看| 精品国产乱码久久久久久小说| 9191精品国产免费久久| 国产乱人偷精品视频| 精品酒店卫生间| 99香蕉大伊视频| 色视频在线一区二区三区| 黄色视频在线播放观看不卡| 男女无遮挡免费网站观看| 成人毛片60女人毛片免费| 中国三级夫妇交换| 亚洲成人手机| 国产深夜福利视频在线观看| 亚洲中文av在线| 免费在线观看黄色视频的| 久久久久久久精品精品| 纯流量卡能插随身wifi吗| 国产成人精品福利久久| 美女主播在线视频| 啦啦啦视频在线资源免费观看| 中文字幕亚洲精品专区| 国产成人精品福利久久| 亚洲av成人精品一二三区| 亚洲国产欧美网| 久久久国产欧美日韩av| 免费观看av网站的网址| 国产精品人妻久久久影院| 777久久人妻少妇嫩草av网站| 三上悠亚av全集在线观看| 18禁国产床啪视频网站| 91国产中文字幕| 国产人伦9x9x在线观看 | 国产男女超爽视频在线观看| 久久久久久伊人网av| 丰满少妇做爰视频| 亚洲美女搞黄在线观看| 丝袜美腿诱惑在线| 大码成人一级视频| 欧美人与性动交α欧美软件| 男女无遮挡免费网站观看| a级毛片在线看网站| 久久久久人妻精品一区果冻| 黄网站色视频无遮挡免费观看| 赤兔流量卡办理| 国产精品成人在线| 亚洲国产欧美网| 亚洲精品一区蜜桃| 精品国产露脸久久av麻豆| 亚洲精品乱久久久久久| 中文乱码字字幕精品一区二区三区| 天天影视国产精品| 国产精品 国内视频| av一本久久久久| 午夜激情av网站| 日韩在线高清观看一区二区三区| 如日韩欧美国产精品一区二区三区| 日日爽夜夜爽网站| 国产在线一区二区三区精| 免费女性裸体啪啪无遮挡网站| xxx大片免费视频| 久久久久国产精品人妻一区二区| 蜜桃在线观看..| 欧美日韩视频精品一区| 啦啦啦视频在线资源免费观看| 丝袜在线中文字幕| 欧美激情极品国产一区二区三区| 亚洲欧洲国产日韩| 久久午夜福利片| 97在线人人人人妻| 满18在线观看网站| 国产精品三级大全| 国产精品国产三级国产专区5o| 亚洲精品美女久久av网站| 日韩不卡一区二区三区视频在线| 午夜91福利影院| 亚洲av欧美aⅴ国产| 性少妇av在线| 女性被躁到高潮视频| 另类亚洲欧美激情| 99香蕉大伊视频| 成年人午夜在线观看视频| 秋霞伦理黄片| av免费观看日本| 午夜福利影视在线免费观看| 午夜免费观看性视频| 两性夫妻黄色片| 久久97久久精品| 老女人水多毛片| 欧美精品人与动牲交sv欧美| 中国国产av一级| 国产精品久久久久久精品古装| 久久av网站| 国产女主播在线喷水免费视频网站| 中国三级夫妇交换| av不卡在线播放| 天天影视国产精品| 久久女婷五月综合色啪小说| 亚洲av男天堂| 久久国产亚洲av麻豆专区| 另类精品久久| www.av在线官网国产| av免费在线看不卡| 亚洲精品第二区| 免费大片黄手机在线观看| 精品久久蜜臀av无| 啦啦啦啦在线视频资源| 亚洲经典国产精华液单| 午夜激情久久久久久久| 在线免费观看不下载黄p国产| 黄色配什么色好看| 大码成人一级视频| 亚洲三区欧美一区| 女人久久www免费人成看片| 有码 亚洲区| 中文字幕制服av| 大话2 男鬼变身卡| 观看美女的网站| 亚洲少妇的诱惑av| 男女啪啪激烈高潮av片| av一本久久久久| 欧美日韩亚洲高清精品| 免费观看av网站的网址| 免费久久久久久久精品成人欧美视频| 97精品久久久久久久久久精品| 国产白丝娇喘喷水9色精品| 亚洲欧洲日产国产| 婷婷色av中文字幕| 欧美最新免费一区二区三区| 午夜激情久久久久久久| 丰满饥渴人妻一区二区三| 国产男女超爽视频在线观看| 午夜激情久久久久久久| 日韩 亚洲 欧美在线| 国产白丝娇喘喷水9色精品| 亚洲综合色网址| 老汉色∧v一级毛片| 色网站视频免费| 国产精品 欧美亚洲| 女的被弄到高潮叫床怎么办| 亚洲成人av在线免费| 国产精品欧美亚洲77777| 黄色 视频免费看| 国产成人精品无人区| 日韩成人av中文字幕在线观看| 日日啪夜夜爽| 天天躁夜夜躁狠狠久久av| 最近中文字幕高清免费大全6| 日本欧美视频一区| 成人黄色视频免费在线看| 欧美国产精品va在线观看不卡| 亚洲在久久综合| 精品久久蜜臀av无| 热re99久久精品国产66热6| 夫妻性生交免费视频一级片| 国产精品三级大全| 国产欧美日韩一区二区三区在线| 精品国产乱码久久久久久男人| 亚洲欧洲日产国产| 久久久久久免费高清国产稀缺| 免费观看在线日韩| 国产成人免费观看mmmm| 视频区图区小说| 中文字幕人妻熟女乱码| 成人亚洲精品一区在线观看| 看十八女毛片水多多多| 老汉色av国产亚洲站长工具| 一区福利在线观看| 日韩电影二区| 纵有疾风起免费观看全集完整版| 免费av中文字幕在线| 又黄又粗又硬又大视频| 另类精品久久| 纵有疾风起免费观看全集完整版| 天美传媒精品一区二区| 巨乳人妻的诱惑在线观看| 少妇被粗大猛烈的视频| 99久久精品国产国产毛片| 亚洲在久久综合| 亚洲精品国产一区二区精华液| 韩国av在线不卡| 国产精品.久久久| 亚洲av电影在线观看一区二区三区| 亚洲经典国产精华液单| av免费观看日本| 亚洲一区中文字幕在线| 成人免费观看视频高清| 国产男女内射视频| 老熟女久久久| 亚洲美女黄色视频免费看| 9191精品国产免费久久| 另类亚洲欧美激情| 国产一区二区在线观看av| 亚洲欧美日韩另类电影网站| 国产有黄有色有爽视频| 精品少妇久久久久久888优播| 欧美97在线视频| 日韩制服骚丝袜av| 啦啦啦啦在线视频资源| 母亲3免费完整高清在线观看 | 激情视频va一区二区三区| 国产女主播在线喷水免费视频网站| 日日爽夜夜爽网站| 18禁国产床啪视频网站| 欧美精品国产亚洲| www.自偷自拍.com| 一区二区三区精品91| 欧美 亚洲 国产 日韩一| 黄色配什么色好看| 久久热在线av| 国产一区有黄有色的免费视频| av福利片在线| 国产成人aa在线观看| 国产av码专区亚洲av| 一区二区三区四区激情视频| 在线观看免费视频网站a站| 一级毛片 在线播放| 国产成人欧美| 亚洲激情五月婷婷啪啪| 中文字幕精品免费在线观看视频| 26uuu在线亚洲综合色| av又黄又爽大尺度在线免费看| 亚洲精品国产av成人精品| 纯流量卡能插随身wifi吗| 男女边摸边吃奶| 久久久久网色| 中文字幕人妻丝袜制服| 久久97久久精品| 蜜桃在线观看..| 欧美中文综合在线视频| 丝袜喷水一区| 亚洲欧美成人综合另类久久久| 老汉色∧v一级毛片| 18在线观看网站| av不卡在线播放| av又黄又爽大尺度在线免费看| 成人漫画全彩无遮挡| 久久婷婷青草| 美女国产视频在线观看| 飞空精品影院首页| 日产精品乱码卡一卡2卡三| 下体分泌物呈黄色| 中文字幕另类日韩欧美亚洲嫩草| 99久久综合免费| 大香蕉久久网| 校园人妻丝袜中文字幕| 18+在线观看网站| 在线观看三级黄色| 久久久a久久爽久久v久久| 日韩制服骚丝袜av| 熟妇人妻不卡中文字幕| 69精品国产乱码久久久| 男人爽女人下面视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 少妇猛男粗大的猛烈进出视频| 中文字幕人妻丝袜一区二区 | 两个人看的免费小视频| 另类亚洲欧美激情| 男女边吃奶边做爰视频| 日韩精品免费视频一区二区三区| xxx大片免费视频| 国产精品久久久av美女十八| 一区二区av电影网| 国产欧美日韩一区二区三区在线| 制服诱惑二区| 亚洲精品在线美女| 免费人妻精品一区二区三区视频| 精品少妇一区二区三区视频日本电影 | 免费久久久久久久精品成人欧美视频| tube8黄色片| 免费不卡的大黄色大毛片视频在线观看| 久久精品国产鲁丝片午夜精品| 免费人妻精品一区二区三区视频| 一级a爱视频在线免费观看| 纵有疾风起免费观看全集完整版| 在现免费观看毛片| 日本免费在线观看一区| videossex国产| 午夜福利在线免费观看网站| av卡一久久| 国产一区二区激情短视频 | 亚洲三区欧美一区| 亚洲一区二区三区欧美精品| a 毛片基地| 交换朋友夫妻互换小说| 9热在线视频观看99| av网站免费在线观看视频| 欧美日韩av久久| 美女大奶头黄色视频| 亚洲国产精品成人久久小说| 国产成人精品久久二区二区91 | 我的亚洲天堂| 久久亚洲国产成人精品v| 在线观看国产h片| 久久毛片免费看一区二区三区| 国产亚洲最大av| a级片在线免费高清观看视频| www.精华液| 日韩电影二区| 中文精品一卡2卡3卡4更新| 日本欧美视频一区| 少妇的丰满在线观看| 丝袜喷水一区| 国产97色在线日韩免费| 丁香六月天网| 久久久久国产一级毛片高清牌| 国产日韩欧美在线精品| 国产精品人妻久久久影院| 国产精品麻豆人妻色哟哟久久| 久久这里只有精品19| 中文字幕亚洲精品专区| 亚洲美女视频黄频| 极品少妇高潮喷水抽搐| 色吧在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 中文精品一卡2卡3卡4更新| 天天影视国产精品| 国产激情久久老熟女| 午夜福利一区二区在线看| 久久影院123| 亚洲精品美女久久av网站| 亚洲精品av麻豆狂野| 90打野战视频偷拍视频| 丰满少妇做爰视频| 欧美xxⅹ黑人| 久久久久网色| 1024视频免费在线观看| 国产日韩欧美视频二区| 国产精品一国产av| 午夜精品国产一区二区电影| 亚洲国产av影院在线观看| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 久久精品aⅴ一区二区三区四区 | 欧美国产精品va在线观看不卡| 狠狠婷婷综合久久久久久88av| 男人爽女人下面视频在线观看| 寂寞人妻少妇视频99o| 国产片特级美女逼逼视频| 亚洲精品日韩在线中文字幕| 国产无遮挡羞羞视频在线观看| 日韩伦理黄色片| 国产亚洲精品第一综合不卡| 熟女少妇亚洲综合色aaa.| 久久亚洲国产成人精品v| 日日摸夜夜添夜夜爱| 熟女电影av网| 国产精品久久久久久精品电影小说| 亚洲激情五月婷婷啪啪| 婷婷色av中文字幕| 国产精品秋霞免费鲁丝片| 久久久a久久爽久久v久久| 日韩av免费高清视频| 高清视频免费观看一区二区| 校园人妻丝袜中文字幕| 亚洲成色77777| 热99国产精品久久久久久7| 日韩av不卡免费在线播放| 性色av一级| 搡女人真爽免费视频火全软件| 午夜福利视频在线观看免费| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕亚洲精品专区| 99热网站在线观看| 美女中出高潮动态图| 青春草亚洲视频在线观看| 中文字幕人妻丝袜一区二区 | 精品久久久精品久久久| 高清黄色对白视频在线免费看| 青春草国产在线视频| 国产国语露脸激情在线看| 欧美精品高潮呻吟av久久| 国产精品免费大片| 人人妻人人爽人人添夜夜欢视频| 97在线人人人人妻| 精品午夜福利在线看| 亚洲人成网站在线观看播放| 美女xxoo啪啪120秒动态图| 免费大片黄手机在线观看| 最近手机中文字幕大全|