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

    MIKE SHE 模型在杭錦旗引黃灌區(qū)的應(yīng)用研究

    2022-10-19 03:45:54韓金旭張翔宇劉姝芳谷曉偉
    人民黃河 2022年10期
    關(guān)鍵詞:杭錦旗引黃灌區(qū)非飽和

    韓金旭,張翔宇,劉姝芳,李 舒,張 丹,谷曉偉

    (1.中國水利水電科學(xué)研究院,北京 100038; 2.黃河水利科學(xué)研究院,河南 鄭州 450003)

    杭錦旗引黃灌區(qū)為自流灌區(qū),灌區(qū)受南部鄂爾多斯臺地補(bǔ)給和北側(cè)黃河影響,地下水水位較高,加之當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)習(xí)慣重灌輕排,灌區(qū)排水不暢,灌區(qū)土壤飽受鹽漬化危害。 鹽漬化土壤要靠大水漫灌壓鹽,漫灌又進(jìn)一步抬高地下水水位,加劇土壤鹽漬化,從而形成惡性循環(huán)。 因此,科學(xué)排水是灌區(qū)除澇、防漬、防止土壤鹽漬化以及改善耕作條件的有效途徑[1-2]。 目前,已有一些關(guān)于引黃灌區(qū)灌溉排水退水規(guī)律的研究,如河套灌區(qū)農(nóng)田排水過程模擬[3]、甘肅景電引黃灌區(qū)退水規(guī)律與模擬[4]等,但有關(guān)內(nèi)蒙古黃河南岸灌區(qū)灌溉排水過程的研究較少,缺少水文參數(shù)的確定研究,同時限制了水文模型的模擬應(yīng)用[5]。 本文以杭錦旗引黃灌區(qū)為研究對象,采用MIKE SHE 模型建立符合灌區(qū)實(shí)際的水文模型,旨在為灌區(qū)排水規(guī)律研究提供模型基礎(chǔ)支撐,為灌區(qū)水文參數(shù)的確定提供參考依據(jù)。

    1 概 況

    杭錦旗位于鄂爾多斯高原西北部,東鄰達(dá)拉特旗、東勝區(qū),西南與烏審旗、鄂托克旗接壤,北靠黃河與巴彥淖爾市隔河相望,屬中溫帶半干旱高原大陸性氣候區(qū),長年干旱少雨,冬季漫長而寒冷,夏季溫和而短促[6]。 杭錦旗引黃灌區(qū)西起三盛公水利樞紐,東至毛不拉孔兌,北臨黃河右岸防洪大堤,南接庫布其沙漠邊緣,沿黃河呈狹長帶狀分布,東西長約170 km,南北寬5~10 km,總灌溉面積2.79 萬hm2。 該灌區(qū)從三盛公水利樞紐上游南岸2.7 km 總干渠進(jìn)水閘引水自流灌溉,總干渠全長204.67 km。

    灌區(qū)受南部鄂爾多斯臺地補(bǔ)給和北側(cè)黃河影響,地下水水位較高,埋深為0.5 ~4.0 m,地下水總體運(yùn)動方向由西南向東北。 為解決灌區(qū)排水及地下水位過高的問題,灌區(qū)建設(shè)了總排干溝,全長195.768 km,貫通全灌區(qū)。 由于灌區(qū)總體地勢較低,排水無法順利排出,因此灌區(qū)配備了6 座常規(guī)排水泵站和2 座臨時排水泵站,將排水及時抽排灌區(qū)。 在歷次續(xù)建配套與節(jié)水改造中由于建設(shè)資金投入不足等,因此排水系統(tǒng)普遍存在建設(shè)標(biāo)準(zhǔn)偏低、配套程度低、抗御自然水旱災(zāi)害能力不強(qiáng)等問題,排水不能及時快速排除,影響了農(nóng)田排水的正常抽排。 根據(jù)杭錦旗引黃灌區(qū)多年灌水?dāng)?shù)據(jù)統(tǒng)計(jì)結(jié)果,2010—2019 年年均灌溉引水量22 574 萬m3,年均排水量3 653 萬m3,排水量總體呈上升趨勢(見圖1),其中2018 年排水量最大,達(dá)到4 450 萬m3。

    圖1 灌區(qū)多年排水量變化情況

    2 模型構(gòu)建

    2.1 模型結(jié)構(gòu)

    MIKE SHE 是一款較為全面的分布式水文模型,它可以模擬陸相水循環(huán)中主要的水文過程,包括降雨、蒸散發(fā)、地表坡面流、地表河道流、地下徑流以及各個水文過程之間的交互作用[7]。 MIKE SHE 模型由數(shù)個獨(dú)立的、相互聯(lián)系的、基于過程的模塊構(gòu)成,主要有坡面流(OL)、河流和湖泊(OC)、非飽和帶/包氣帶(UZ)、蒸散發(fā)(ET)與飽和帶(SZ)等5 個子模塊。 每個子模塊用于一個主要水文過程的描述,根據(jù)不同的模擬要求,這些子模塊可以互相分離開來也可以綜合起來應(yīng)用,分離開來可分別描述流域水文循環(huán)的各個過程,綜合起來可描述整個流域的水文循環(huán)過程[8]。模型中坡面流模擬采用二維圣維南方程組,河流和湖泊模塊為MIKE 11 圣維南一維水動力學(xué)模型,飽和帶即地下水的模擬采用集總式的線性水庫法,非飽和帶模擬采用水量平衡法或理査德方程[9]。 MIKE SHE 相較于其他水文模型具有明顯的優(yōu)勢,可以靈活地選用必要的模塊以匹配灌區(qū)特點(diǎn),能夠清晰地模擬灌區(qū)水分的運(yùn)動和轉(zhuǎn)化過程,在此基礎(chǔ)上建立灌區(qū)水文循環(huán)模型,進(jìn)而研究灌區(qū)排水形成過程是切實(shí)可行的。

    2.2 數(shù)據(jù)庫搭建

    灌區(qū)排水過程影響因素繁多,受氣象因素、地下水運(yùn)動條件、土壤屬性、地形地貌等自然因素和引水量、灌溉制度、作物種植結(jié)構(gòu)、渠系構(gòu)造等人為因素共同影響,其形成過程及變化規(guī)律較為復(fù)雜[10-11]。 結(jié)合杭錦旗引黃灌區(qū)水循環(huán)特點(diǎn),搭建灌區(qū)水文模型數(shù)據(jù)庫,主要包括模型范圍及地表高程、氣象數(shù)據(jù)、土地利用分布、灌溉數(shù)據(jù)、排干溝設(shè)置、坡面流數(shù)據(jù)、土壤數(shù)據(jù)、地下水條件等方面的基礎(chǔ)數(shù)據(jù)及參數(shù)。 因文章篇幅有限,在此選擇模型范圍及地表高程、土地利用類型分布、土壤數(shù)據(jù)等基礎(chǔ)數(shù)據(jù)進(jìn)行重點(diǎn)說明,其他數(shù)據(jù)僅簡要介紹。

    本次研究收集到灌區(qū)2017—2019 年實(shí)測逐日排水資料,在此確定模型的模擬時段為2017—2019 年,其中2017—2018 年為模型數(shù)據(jù)率定校準(zhǔn)階段、2019年為模型數(shù)據(jù)驗(yàn)證階段。

    2.2.1 模型范圍及地表高程

    根據(jù)灌區(qū)分布狀況,在ArcGIS 中確定灌區(qū)范圍shp 面文件,并利用ArcGIS 中投影轉(zhuǎn)換工具將坐標(biāo)系轉(zhuǎn)換為CGCS2000_3_Degree_GK_CM_108E,模型范圍總面積為1 620 km2。 灌區(qū)高程采用30 m 分辨率數(shù)字高程數(shù)據(jù)集。 根據(jù)模型運(yùn)行精度需求及現(xiàn)實(shí)情況,模型計(jì)算空間網(wǎng)格大小選取500 m×500 m,模型范圍起始x、y坐標(biāo)點(diǎn)為(416 000,4 460 000),將模型范圍沿x和y方向劃分為350×140 個網(wǎng)格。

    2.2.2 土地利用類型分布

    根據(jù)國家土地資源分類標(biāo)準(zhǔn)并結(jié)合當(dāng)?shù)貙?shí)際情況,將模型范圍內(nèi)土地資源分為耕地、草地、林地、裸地、居民區(qū)、水體及其他用地。 模型范圍內(nèi)的各土地利用類型占地面積統(tǒng)計(jì)見表1。

    表1 土地利用類型統(tǒng)計(jì)

    土地利用模塊植被屬性子項(xiàng)中需要在植被屬性文件中輸入各類用地的葉面積指數(shù)(LAI)及根系深度(RD)時間序列。 根據(jù)FAO 文件及查閱相關(guān)文獻(xiàn)資料確定不同土地類型植被根系、葉面積指數(shù)隨時間變化規(guī)律,得到研究區(qū)域的作物參數(shù)值,作為模型的輸入項(xiàng)。 灌區(qū)播種作物主要為玉米、葵花,2019 年兩種作物種植面積占總耕地面積的90%,此外還種植有少量小麥、雜糧、牧草等。 模型中按玉米、葵花兩種作物考慮,玉米、葵花種植面積比為4 ∶6。 查閱相關(guān)文獻(xiàn)資料得到作物生育期及對應(yīng)的LAI、Kc(作物系數(shù))與RD時間序列[12],并根據(jù)兩種作物種植面積權(quán)重比得到耕地類型對應(yīng)的植被屬性參數(shù),時間序列見圖2。

    2.2.3 土壤數(shù)據(jù)

    非飽和帶即包氣帶是指以地表為上邊界、地下水位為下邊界的非飽和水流區(qū)域。 非飽和帶土壤含水量和土壤水勢的變化受非飽和水流運(yùn)動過程的影響,此外,地表水下滲至非飽和帶土壤中的水量以及非飽和帶與飽和帶之間的水量交換過程也受到非飽和帶水流運(yùn)動過程的控制。 非飽和帶上部受降水、蒸發(fā)和下滲的影響,土壤含水量波動較大,下邊界受地下水位的波動發(fā)生相應(yīng)的改變。 本次模擬采用兩層水量平衡法(2 Layer UZ),根據(jù)國家土壤質(zhì)地分類標(biāo)準(zhǔn),將模型范圍內(nèi)的土壤類型分為砂質(zhì)壤土、砂質(zhì)黏壤土、黏壤土3種類型(見表2),主要類型為砂質(zhì)壤土。

    表2 不同土壤質(zhì)地顆粒含量 %

    2.2.4 其他數(shù)據(jù)

    氣象數(shù)據(jù)采用臨河國家氣象站2017—2019 年數(shù)據(jù),主要包括風(fēng)速、降水量、氣溫、日照時數(shù)和相對濕度等數(shù)據(jù),用Penman-Monteith 公式,計(jì)算參考作物蒸散發(fā)量ET0。 該灌區(qū)采用地面灌溉方式,灌溉水源為黃河,因此在模型中灌溉水源設(shè)定為外部水源,灌溉過程采用2017—2019 年灌區(qū)實(shí)際引水?dāng)?shù)據(jù)。 灌區(qū)總排干溝全長195.768 km,溝深2 m,內(nèi)邊坡按1 ∶3 考慮,縱比降為1/4 500~1/7 500,排干采用一維河道模型MIKE 11模擬,將MIKE 11 模型與MIKE SHE 耦合形成完整的地表-地下水循環(huán)模擬模型。 坡面流是地表積水的運(yùn)動,一般采用有限差分法模擬坡面流,需要輸入3 個參數(shù),即曼寧系數(shù)作為率定參數(shù)、滯蓄水深設(shè)置為田埂高度500 mm、初始地表水深為0。

    2.3 模型輸入

    在輸入數(shù)據(jù)之前,需先對模型進(jìn)行初始設(shè)置。 首先,根據(jù)研究目的以及灌區(qū)實(shí)際情況選定模型的子模塊,選定MIKE 模型中的坡面流(OL)、河流和湖泊(OC)、非飽和帶/包氣帶(UZ)、蒸散發(fā)(ET)、飽和帶(SZ)5 個子模塊;然后,對模擬的時段進(jìn)行設(shè)置,根據(jù)基礎(chǔ)數(shù)據(jù)詳細(xì)程度,本次模擬時段為2017-01-01 至2019-12-31,以a 為單位,分別進(jìn)行模擬,MIKE SHE模型的模擬時間步長為24 h,MIKE 11 河道模擬時間步長為10 min,其他模擬參數(shù)的設(shè)置采用模型默認(rèn)參數(shù);最后,在相關(guān)設(shè)置完成后,導(dǎo)入已生成的地表高程、氣象數(shù)據(jù)、土地利用類型分布、灌溉數(shù)據(jù)、坡面流參數(shù)、土壤參數(shù)、地下水分布條件等相關(guān)文件并設(shè)置相關(guān)參數(shù)。

    3 結(jié)果分析

    3.1 參數(shù)率定及驗(yàn)證

    模型中許多參數(shù)由于缺乏實(shí)測數(shù)據(jù),因此基于理論參考范圍和經(jīng)驗(yàn)賦予了初始值,需要根據(jù)模擬值與實(shí)測值的擬合程度進(jìn)行參數(shù)調(diào)整率定。 在此選取對模型模擬結(jié)果有較大影響的參數(shù)作為率定參數(shù),主要包括河道水動力、坡面流、非飽和帶土壤屬性、含水層水文地質(zhì)等參數(shù),并選取相關(guān)系數(shù)(R)、均方差(RMSE)、納什效率系數(shù)(E)3 個逐日排水過程指標(biāo)及年排水量相對誤差評價模擬結(jié)果的準(zhǔn)確程度[13]。 模型率定選取的參數(shù)及率定結(jié)果見表3。

    表3 模型率定選取的參數(shù)及率定結(jié)果

    3.2 模擬結(jié)果分析

    根據(jù)以上參數(shù)率定結(jié)果,對2019 年排水過程進(jìn)行驗(yàn)證模擬。 最終率定期及驗(yàn)證期的模擬結(jié)果評價見表4,驗(yàn)證期2019 年排水流量模擬結(jié)果見圖3。

    圖3 2019 年排水流量模擬結(jié)果

    表4 模型模擬結(jié)果評價

    可以看出,模擬結(jié)果存在兩個問題:一是平均誤差較大,總體評價指標(biāo)不夠優(yōu)化;二是2019 年作為驗(yàn)證期,模擬結(jié)果明顯優(yōu)于率定期。 下面針對這兩個問題進(jìn)行分析討論。

    (1)模型率定驗(yàn)證所采用實(shí)測值為灌區(qū)實(shí)際抽排量,排水閘門及泵站需人工開啟關(guān)閉,實(shí)測排水量受到監(jiān)測手段及人為因素影響;模型模擬值為受降水、灌溉及地下水側(cè)向徑流影響的天然排水量。 本次模擬值與實(shí)測值本質(zhì)上有一定差別,模擬值均大于實(shí)測值,因此平均誤差較大,評價指標(biāo)看似不夠優(yōu)化。

    (2)與2017 年、2018 年相比,2019 年排水設(shè)備相對完善,人員配備到位,具有較好的監(jiān)測計(jì)量條件,且增加了臨時排水站,排水能及時地抽出溝道,受人為因素影響較小,因此驗(yàn)證期2019 年模擬效果相對較好,相關(guān)性較高、誤差較小。

    4 結(jié) 語

    (1)對于水文模型來說,所需基礎(chǔ)資料多且復(fù)雜,部分?jǐn)?shù)據(jù)又受到人為因素影響,給模型模擬帶來了不確定性,因此模型有一定誤差是允許的;同時2019 年的模擬結(jié)果表明,在實(shí)測數(shù)據(jù)受干擾影響較小的情況下,模擬值與實(shí)測值擬合程度較高,模型可以較好地模擬灌區(qū)天然排水過程,認(rèn)為模型模擬結(jié)果是較為合理的,近3 a 天然排水量應(yīng)接近4 961 萬m3。

    (2)本次模型以2017—2018 年為參數(shù)率定期、2019 年為驗(yàn)證期,主要依據(jù)杭錦旗引黃灌區(qū)實(shí)測排水?dāng)?shù)據(jù)進(jìn)行率定驗(yàn)證,率定參數(shù)具有參考性,但在時間、空間及率定要素選取上也有一定片面性,相關(guān)模型研究可以在本次模擬結(jié)果的基礎(chǔ)上進(jìn)行不斷率定修正,以期不斷完善相關(guān)水文模型及參數(shù)研究。

    (3)灌區(qū)現(xiàn)狀排水直接排入黃河干流,由于水質(zhì)不達(dá)標(biāo),該部分水量不被算作退水量,即仍然占用黃河耗水指標(biāo),若將灌區(qū)排水納入?yún)^(qū)域水資源統(tǒng)一配置,則既可以減少向黃河排污、又能夠?yàn)閰^(qū)域經(jīng)濟(jì)社會發(fā)展開源。 在灌區(qū)節(jié)水及渠系改造的背景下,下一步可以針對不同模擬情景預(yù)測未來排水量,為灌區(qū)排水資源化利用提供支撐。

    猜你喜歡
    杭錦旗引黃灌區(qū)非飽和
    趙口引黃灌區(qū)邁向高質(zhì)量發(fā)展之路
    人民黃河(2022年11期)2022-11-28 15:10:36
    非飽和原狀黃土結(jié)構(gòu)強(qiáng)度的試驗(yàn)研究
    杭錦旗地區(qū)辮狀河定量地質(zhì)知識庫建立及應(yīng)用
    從姓氏的角度看杭錦旗名稱演變
    民族地區(qū)本土音樂文化教育傳承現(xiàn)狀與思考——以杭錦旗“古如歌”教育傳承基地為例
    人民音樂(2017年7期)2017-07-19 13:03:04
    非飽和多孔介質(zhì)應(yīng)力滲流耦合分析研究
    梁山引黃灌區(qū)渠道防滲效果與灌溉用水效率分析
    非飽和土基坑剛性擋墻抗傾覆設(shè)計(jì)與參數(shù)分析
    引黃灌區(qū)泥沙處理與利用技術(shù)發(fā)展現(xiàn)狀及分析
    科技視界(2016年10期)2016-04-26 12:31:19
    非飽和地基土蠕變特性試驗(yàn)研究
    免费人妻精品一区二区三区视频| 精品久久久久久久久亚洲| 免费高清在线观看视频在线观看| 日本-黄色视频高清免费观看| 最近中文字幕2019免费版| 性高湖久久久久久久久免费观看| 大香蕉久久网| 久久鲁丝午夜福利片| 午夜精品国产一区二区电影| 蜜桃久久精品国产亚洲av| 亚洲精品一区蜜桃| 高清欧美精品videossex| 亚洲aⅴ乱码一区二区在线播放| 亚洲av.av天堂| 99久久精品热视频| 26uuu在线亚洲综合色| 美女主播在线视频| 26uuu在线亚洲综合色| 国产伦精品一区二区三区视频9| 日本wwww免费看| 久久久久久久精品精品| 如何舔出高潮| 亚洲中文av在线| 国产真实伦视频高清在线观看| 国产精品爽爽va在线观看网站| 久久久色成人| 国产欧美日韩精品一区二区| 成人毛片a级毛片在线播放| 日日撸夜夜添| 国产亚洲精品久久久com| 高清在线视频一区二区三区| 国产白丝娇喘喷水9色精品| 久久国产乱子免费精品| 超碰av人人做人人爽久久| 国产精品一区www在线观看| 老女人水多毛片| 国产在线免费精品| 亚洲成人中文字幕在线播放| 啦啦啦中文免费视频观看日本| 成年免费大片在线观看| 熟女av电影| 免费观看在线日韩| 成人毛片60女人毛片免费| 亚洲精品久久久久久婷婷小说| 国产真实伦视频高清在线观看| 国产一区二区三区av在线| 在线观看一区二区三区| 有码 亚洲区| 最近中文字幕高清免费大全6| 久久97久久精品| 国产成人精品福利久久| 国产av国产精品国产| 爱豆传媒免费全集在线观看| 在线观看美女被高潮喷水网站| videossex国产| 亚洲国产色片| 欧美高清性xxxxhd video| 成年女人在线观看亚洲视频| 99久久精品一区二区三区| 久久这里有精品视频免费| 久久毛片免费看一区二区三区| 日韩,欧美,国产一区二区三区| 欧美三级亚洲精品| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美中文字幕日韩二区| 久久毛片免费看一区二区三区| 这个男人来自地球电影免费观看 | 日韩,欧美,国产一区二区三区| 亚洲熟女精品中文字幕| 91精品国产国语对白视频| 男人狂女人下面高潮的视频| 国产一区二区三区综合在线观看 | 亚洲婷婷狠狠爱综合网| 极品少妇高潮喷水抽搐| 欧美xxⅹ黑人| 一本—道久久a久久精品蜜桃钙片| 美女xxoo啪啪120秒动态图| 国产精品.久久久| 欧美日本视频| 国产av精品麻豆| 婷婷色麻豆天堂久久| 国产成人免费观看mmmm| 一级a做视频免费观看| 女人十人毛片免费观看3o分钟| 大香蕉97超碰在线| 国产无遮挡羞羞视频在线观看| 国产乱人偷精品视频| 在线观看免费视频网站a站| 久久精品国产亚洲av涩爱| 国产免费福利视频在线观看| 成人黄色视频免费在线看| 在线观看三级黄色| av在线播放精品| 少妇人妻精品综合一区二区| 中文字幕制服av| 亚洲丝袜综合中文字幕| 欧美激情极品国产一区二区三区 | 国产精品福利在线免费观看| 九九久久精品国产亚洲av麻豆| 亚洲欧洲日产国产| 亚洲,一卡二卡三卡| 国产淫片久久久久久久久| 免费看日本二区| av在线app专区| 亚洲国产精品999| 天堂中文最新版在线下载| 五月伊人婷婷丁香| 多毛熟女@视频| 午夜福利影视在线免费观看| 成人二区视频| 久久毛片免费看一区二区三区| 久久热精品热| 丰满人妻一区二区三区视频av| 哪个播放器可以免费观看大片| 联通29元200g的流量卡| 少妇的逼水好多| 在线观看av片永久免费下载| 国产欧美另类精品又又久久亚洲欧美| 亚州av有码| 日韩制服骚丝袜av| 国产精品国产三级专区第一集| 国产大屁股一区二区在线视频| 亚洲国产毛片av蜜桃av| 国产精品人妻久久久久久| 黄色一级大片看看| 精品熟女少妇av免费看| 国产黄频视频在线观看| 色吧在线观看| 一二三四中文在线观看免费高清| 女的被弄到高潮叫床怎么办| 亚洲欧美一区二区三区黑人 | 久久久久久久精品精品| 成人美女网站在线观看视频| 欧美亚洲 丝袜 人妻 在线| 国产精品女同一区二区软件| 老师上课跳d突然被开到最大视频| 亚洲电影在线观看av| 一级a做视频免费观看| 久久久久性生活片| 一级毛片aaaaaa免费看小| 国产黄片美女视频| 久久久久久九九精品二区国产| 久久久午夜欧美精品| 久久久久精品久久久久真实原创| 国产精品免费大片| 婷婷色综合www| 国产色婷婷99| 黄片wwwwww| 一级毛片我不卡| 联通29元200g的流量卡| 欧美日韩视频精品一区| 久久久色成人| 黄色配什么色好看| 欧美日韩精品成人综合77777| 国产精品秋霞免费鲁丝片| 国产精品99久久久久久久久| 国产精品一区二区在线观看99| 国产亚洲5aaaaa淫片| 精品一区二区三区视频在线| 国产日韩欧美亚洲二区| 欧美最新免费一区二区三区| 亚洲av成人精品一二三区| 一级a做视频免费观看| 中文字幕精品免费在线观看视频 | 99热国产这里只有精品6| 亚洲国产最新在线播放| 中国美白少妇内射xxxbb| 男女下面进入的视频免费午夜| 午夜福利影视在线免费观看| 十八禁网站网址无遮挡 | 我要看黄色一级片免费的| 人人妻人人看人人澡| 日韩中字成人| 熟女人妻精品中文字幕| 人人妻人人澡人人爽人人夜夜| 亚洲av福利一区| 久久精品久久久久久久性| 老司机影院毛片| 国产男女内射视频| 成年免费大片在线观看| 亚洲人成网站在线观看播放| 国产伦精品一区二区三区四那| 一级黄片播放器| 99久久人妻综合| 亚洲国产日韩一区二区| 免费观看a级毛片全部| 观看美女的网站| 啦啦啦啦在线视频资源| 97超碰精品成人国产| 99热国产这里只有精品6| 日本欧美国产在线视频| www.av在线官网国产| 国产一区亚洲一区在线观看| 欧美日韩亚洲高清精品| 国产女主播在线喷水免费视频网站| 青青草视频在线视频观看| 亚洲av不卡在线观看| 尾随美女入室| 两个人的视频大全免费| 国产女主播在线喷水免费视频网站| 国内揄拍国产精品人妻在线| 天堂8中文在线网| 国产亚洲精品久久久com| 99热这里只有是精品在线观看| 在线观看免费视频网站a站| 欧美日韩一区二区视频在线观看视频在线| 国产精品伦人一区二区| 日韩视频在线欧美| 免费观看av网站的网址| 18禁裸乳无遮挡动漫免费视频| 久久韩国三级中文字幕| 在线 av 中文字幕| 观看免费一级毛片| 亚洲成色77777| 国产黄频视频在线观看| 欧美国产精品一级二级三级 | 九色成人免费人妻av| av不卡在线播放| 国产一区有黄有色的免费视频| 美女脱内裤让男人舔精品视频| 亚洲伊人久久精品综合| 六月丁香七月| 亚洲美女搞黄在线观看| 五月开心婷婷网| 成人国产麻豆网| 卡戴珊不雅视频在线播放| 精品99又大又爽又粗少妇毛片| 国产一区二区三区综合在线观看 | 免费观看无遮挡的男女| 国产91av在线免费观看| 网址你懂的国产日韩在线| 亚洲av成人精品一二三区| 99久久综合免费| 最近2019中文字幕mv第一页| 国产在线一区二区三区精| 日韩,欧美,国产一区二区三区| 精华霜和精华液先用哪个| 日本猛色少妇xxxxx猛交久久| 午夜免费鲁丝| 免费黄色在线免费观看| 亚洲av不卡在线观看| 精品酒店卫生间| 日韩不卡一区二区三区视频在线| 久久6这里有精品| 中文字幕亚洲精品专区| 欧美激情极品国产一区二区三区 | 国产成人a∨麻豆精品| 丝瓜视频免费看黄片| 女人久久www免费人成看片| 亚洲伊人久久精品综合| 最后的刺客免费高清国语| 最近最新中文字幕大全电影3| 寂寞人妻少妇视频99o| 国产男女内射视频| 最近2019中文字幕mv第一页| 日韩视频在线欧美| av网站免费在线观看视频| 十分钟在线观看高清视频www | 美女xxoo啪啪120秒动态图| 狂野欧美白嫩少妇大欣赏| 午夜免费观看性视频| 午夜福利网站1000一区二区三区| 一级毛片 在线播放| 免费观看的影片在线观看| 亚洲图色成人| 看十八女毛片水多多多| 欧美日韩视频精品一区| 青青草视频在线视频观看| 99精国产麻豆久久婷婷| 男女无遮挡免费网站观看| 熟妇人妻不卡中文字幕| 国产伦理片在线播放av一区| 日本免费在线观看一区| 婷婷色综合www| 亚洲不卡免费看| 色网站视频免费| 精品国产三级普通话版| 精品久久久久久久末码| 中文字幕制服av| 亚洲中文av在线| 一本色道久久久久久精品综合| 午夜福利视频精品| 高清毛片免费看| 99精国产麻豆久久婷婷| 国产午夜精品久久久久久一区二区三区| 久久ye,这里只有精品| 国产日韩欧美亚洲二区| 人体艺术视频欧美日本| av播播在线观看一区| 啦啦啦视频在线资源免费观看| 国产日韩欧美亚洲二区| av网站免费在线观看视频| 国产精品精品国产色婷婷| 中文字幕亚洲精品专区| 最近中文字幕2019免费版| 日本猛色少妇xxxxx猛交久久| 美女高潮的动态| 黄色日韩在线| 自拍欧美九色日韩亚洲蝌蚪91 | 美女cb高潮喷水在线观看| 欧美一级a爱片免费观看看| 午夜福利在线在线| 卡戴珊不雅视频在线播放| 国产成人精品一,二区| 日日啪夜夜撸| 亚洲成人手机| 男女无遮挡免费网站观看| 亚洲熟女精品中文字幕| 高清欧美精品videossex| 国产精品国产av在线观看| 亚洲真实伦在线观看| 久久 成人 亚洲| 黄色欧美视频在线观看| 国产精品福利在线免费观看| 99久久综合免费| 在线观看免费视频网站a站| 成人黄色视频免费在线看| 国语对白做爰xxxⅹ性视频网站| h视频一区二区三区| 亚洲成人中文字幕在线播放| 一本久久精品| 亚洲av中文字字幕乱码综合| 黄片wwwwww| 成年人午夜在线观看视频| 亚洲在久久综合| 最近最新中文字幕免费大全7| 久久人人爽人人片av| 女人久久www免费人成看片| 久久人人爽人人爽人人片va| 亚洲av综合色区一区| 全区人妻精品视频| 亚洲av电影在线观看一区二区三区| 亚洲不卡免费看| 精品人妻偷拍中文字幕| 亚洲aⅴ乱码一区二区在线播放| 久久国产精品男人的天堂亚洲 | 国产黄片视频在线免费观看| 亚洲国产毛片av蜜桃av| av播播在线观看一区| 亚洲精品国产av成人精品| 人人妻人人看人人澡| 精品久久久久久电影网| 水蜜桃什么品种好| 成人亚洲精品一区在线观看 | 99热这里只有是精品在线观看| 久久国产精品大桥未久av | 天天躁日日操中文字幕| 男人添女人高潮全过程视频| 亚洲欧美精品专区久久| 中文字幕精品免费在线观看视频 | 久热这里只有精品99| 黄片wwwwww| 熟女人妻精品中文字幕| www.av在线官网国产| xxx大片免费视频| 久久久久性生活片| 少妇的逼水好多| 久久久午夜欧美精品| 久久韩国三级中文字幕| 久久精品熟女亚洲av麻豆精品| 亚洲精品乱码久久久久久按摩| 美女福利国产在线 | 亚洲精品乱码久久久久久按摩| 天美传媒精品一区二区| 亚洲美女搞黄在线观看| 国内少妇人妻偷人精品xxx网站| 国产av国产精品国产| 中文字幕久久专区| 日本-黄色视频高清免费观看| 超碰av人人做人人爽久久| 亚洲国产毛片av蜜桃av| 观看av在线不卡| 国产精品熟女久久久久浪| 人人妻人人爽人人添夜夜欢视频 | 在线天堂最新版资源| 欧美性感艳星| 亚洲av成人精品一区久久| 日韩伦理黄色片| 国产 精品1| 国产精品一区二区在线观看99| 久久ye,这里只有精品| 97在线人人人人妻| 欧美精品一区二区免费开放| 精品一区二区三区视频在线| 少妇高潮的动态图| 国内揄拍国产精品人妻在线| 亚洲国产精品成人久久小说| 18+在线观看网站| 久久久久久伊人网av| 最近2019中文字幕mv第一页| 色综合色国产| 十分钟在线观看高清视频www | 久久6这里有精品| 精品视频人人做人人爽| 日韩不卡一区二区三区视频在线| 免费黄色在线免费观看| 日韩电影二区| 看非洲黑人一级黄片| 国产美女午夜福利| 一个人看的www免费观看视频| 亚洲精品国产成人久久av| 久久久久久久久久久丰满| 久久99精品国语久久久| 国产精品熟女久久久久浪| 看十八女毛片水多多多| 我要看日韩黄色一级片| 一个人免费看片子| 成人黄色视频免费在线看| 精品人妻熟女av久视频| 色网站视频免费| 黄色怎么调成土黄色| 久久久精品免费免费高清| 日韩大片免费观看网站| 亚洲精品456在线播放app| 两个人的视频大全免费| 欧美成人一区二区免费高清观看| 99re6热这里在线精品视频| 精品久久国产蜜桃| videos熟女内射| 人妻夜夜爽99麻豆av| 日韩欧美精品免费久久| 黄色配什么色好看| 丰满乱子伦码专区| 一个人免费看片子| 五月伊人婷婷丁香| 亚洲aⅴ乱码一区二区在线播放| 国产色婷婷99| 日本vs欧美在线观看视频 | 高清毛片免费看| 亚洲欧美日韩无卡精品| 亚洲欧美清纯卡通| 九草在线视频观看| 激情五月婷婷亚洲| 国产成人精品一,二区| 少妇被粗大猛烈的视频| 日韩 亚洲 欧美在线| 美女脱内裤让男人舔精品视频| 欧美xxⅹ黑人| 亚洲国产日韩一区二区| 夫妻午夜视频| 成人美女网站在线观看视频| 亚洲熟女精品中文字幕| 少妇 在线观看| h日本视频在线播放| 一边亲一边摸免费视频| 亚洲美女搞黄在线观看| 永久免费av网站大全| 久久午夜福利片| 亚洲最大成人中文| 国产精品女同一区二区软件| av专区在线播放| 亚洲怡红院男人天堂| 欧美日韩亚洲高清精品| 久久99精品国语久久久| 毛片女人毛片| 久久女婷五月综合色啪小说| 色视频在线一区二区三区| 国产午夜精品一二区理论片| 亚洲自偷自拍三级| 久久久久久久大尺度免费视频| 日产精品乱码卡一卡2卡三| 久久国内精品自在自线图片| 女性生殖器流出的白浆| 黄色怎么调成土黄色| 伦精品一区二区三区| 99九九线精品视频在线观看视频| 亚洲精品一二三| 成人18禁高潮啪啪吃奶动态图 | 在线观看一区二区三区| 国产精品久久久久久久久免| 亚洲精品自拍成人| 久久精品夜色国产| 精品人妻偷拍中文字幕| 亚洲欧美一区二区三区国产| 22中文网久久字幕| av不卡在线播放| 大话2 男鬼变身卡| 日韩免费高清中文字幕av| 联通29元200g的流量卡| 国产精品国产三级国产av玫瑰| 国产一级毛片在线| 久久99蜜桃精品久久| a级毛色黄片| 国产精品福利在线免费观看| 妹子高潮喷水视频| 国产精品无大码| 韩国高清视频一区二区三区| 久久久午夜欧美精品| 免费黄色在线免费观看| 亚洲精品亚洲一区二区| 国产美女午夜福利| 成人毛片60女人毛片免费| 亚洲欧美清纯卡通| 老司机影院毛片| 欧美zozozo另类| 国产毛片在线视频| 搡老乐熟女国产| 亚洲av中文字字幕乱码综合| 国产成人精品婷婷| 免费少妇av软件| 中国三级夫妇交换| 极品少妇高潮喷水抽搐| 亚洲伊人久久精品综合| 高清毛片免费看| 视频区图区小说| 女的被弄到高潮叫床怎么办| 日本午夜av视频| 干丝袜人妻中文字幕| 国产精品国产三级专区第一集| 狂野欧美激情性bbbbbb| 永久免费av网站大全| 一二三四中文在线观看免费高清| 日韩在线高清观看一区二区三区| 亚洲精品国产色婷婷电影| 在线观看免费视频网站a站| 又粗又硬又长又爽又黄的视频| 国产亚洲av片在线观看秒播厂| 久久久成人免费电影| 久久精品国产亚洲av天美| 久久精品久久久久久噜噜老黄| 亚洲精品456在线播放app| 日本av免费视频播放| 成年免费大片在线观看| 亚洲熟女精品中文字幕| 亚洲精品456在线播放app| 少妇人妻一区二区三区视频| 一级av片app| 日本wwww免费看| 欧美成人精品欧美一级黄| 最新中文字幕久久久久| 久久久久久久国产电影| 黄片wwwwww| 熟妇人妻不卡中文字幕| 国产精品av视频在线免费观看| 九九爱精品视频在线观看| 精品国产乱码久久久久久小说| 女的被弄到高潮叫床怎么办| av国产免费在线观看| 亚洲国产色片| 如何舔出高潮| 亚洲精品视频女| 熟女人妻精品中文字幕| 久久久久久九九精品二区国产| 久久精品夜色国产| 日本与韩国留学比较| 亚洲精品国产av蜜桃| 日韩制服骚丝袜av| 纯流量卡能插随身wifi吗| 少妇被粗大猛烈的视频| 精品99又大又爽又粗少妇毛片| 成人国产麻豆网| 国产色爽女视频免费观看| 亚洲一级一片aⅴ在线观看| 日韩电影二区| 男女边吃奶边做爰视频| 麻豆成人av视频| 超碰av人人做人人爽久久| 国产精品人妻久久久久久| 国产高清国产精品国产三级 | 最黄视频免费看| 午夜激情福利司机影院| 亚洲真实伦在线观看| 国产69精品久久久久777片| 精品久久久久久久久亚洲| 亚洲国产精品国产精品| 亚洲欧美成人综合另类久久久| 亚洲图色成人| 亚洲国产精品专区欧美| 伊人久久国产一区二区| 哪个播放器可以免费观看大片| 久久久国产一区二区| 精品一区二区三卡| 丰满乱子伦码专区| 亚洲,欧美,日韩| 国产成人免费无遮挡视频| 老司机影院毛片| 日本欧美国产在线视频| 天天躁夜夜躁狠狠久久av| 建设人人有责人人尽责人人享有的 | 欧美老熟妇乱子伦牲交| 亚洲精品久久午夜乱码| 久久久久久久久久久免费av| 亚洲va在线va天堂va国产| 午夜福利网站1000一区二区三区| 免费黄色在线免费观看| 哪个播放器可以免费观看大片| 99热这里只有是精品在线观看| 亚洲美女黄色视频免费看| 国产爽快片一区二区三区| 成人国产av品久久久| 欧美激情极品国产一区二区三区 | 国产一区有黄有色的免费视频| 九色成人免费人妻av| 中文资源天堂在线| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产高清在线一区二区三| 丝袜脚勾引网站| 国产黄频视频在线观看| 成年av动漫网址| av又黄又爽大尺度在线免费看| 亚洲人成网站在线观看播放| 女人十人毛片免费观看3o分钟| 欧美三级亚洲精品| 三级国产精品片| 国内揄拍国产精品人妻在线| 夜夜爽夜夜爽视频| 嫩草影院新地址| 你懂的网址亚洲精品在线观看| 简卡轻食公司| 日本午夜av视频| 国产成人精品一,二区| 一级二级三级毛片免费看| .国产精品久久| 最近2019中文字幕mv第一页| 3wmmmm亚洲av在线观看|