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

    黑土區(qū)土壤水分特征曲線模擬及模型優(yōu)選

    2018-10-19 05:10:40王子龍常廣義姜秋香陳偉杰林百健
    東北農業(yè)大學學報 2018年9期
    關鍵詞:壤土吸力實測值

    王子龍,常廣義,姜秋香,付 強,陳偉杰,林百健

    (東北農業(yè)大學水利與土木工程學院,哈爾濱 150030)

    土壤水分特征曲線是模擬土壤水分和溶質運移重要輸入參數,廣泛應用于地下水補給、農業(yè)與土壤化學、水土保持等領域,可反映土壤孔隙分布狀況[1-3]。

    測定土壤水分特征曲線費時長、成本高、耗費人力,常用土壤水分特征曲線模型有Brooks-Corey(BC)[4],Gardner[5]、van Genuchten[6](VG)及Kosugi[7]、Fredlund 和 Xing[8]、Dual-porosity模型[9]等。土壤傳遞函數模型分為點預測模型和參數預測模型[10]。點預測模型利用田間含水量33 kPa、凋萎含水量1 500 kPa和土壤水力性質構建函數關系;參數預測模型基于模型參數與物理性質(如土壤質地、有機質、顆粒組成等)構建函數關系。土壤傳遞函數模型研究主要有線性和非線性擬合、人工神經網絡[11]等方法。土壤傳遞函數模型成果有Vereecken[12]模型,利用土壤中顆粒組成、有機碳等與VG模型參數建立關系式,Cosby[13]和Saxton[14]模型利用顆粒組成與實測土壤水分特征曲線參數作回歸分析,姚姣轉等在VG模型基礎上建立科爾沁沙地土壤傳遞函數等[15]。

    比較分析土壤水分特征曲線傳遞函數模型和經驗模型預測能力,如Woston等比較21種土壤傳遞函數,結果表明點預測模型模擬精度優(yōu)于參數預測模型[16];Kern分析對比6種土壤傳遞函數模型,發(fā)現Rawls模型預測精度最高[17];Tietje等評價13種土壤傳遞函數預測效果,結果表明Vereecken模型對土壤水分特征曲線預測值與實測值偏差最小,模擬精度最高[18];Matlan等對比4種經驗模型,結果表明MG模型對不同質地土壤土壤水分特征曲線模擬效果最好[19];丁新原等利用VG和Gardner模型模擬塔里木沙漠公路防護林SWCC,發(fā)現VG模型模擬精度高于Gardner模型[20]。目前,土壤水分特征曲線模型對比研究較多,但模型普適性和精確度依賴土壤類型,對比分析特定土壤適應性尤為必要。

    本文以松嫩平原黑土區(qū)黑土為例,利用VG、BC、LND、MG 4種土壤水分特征曲線模型模擬不同質地黑土土壤水分特征曲線,對比分析4種經驗模型模擬精度,選取最優(yōu)模型。研究結果為松嫩平原黑土區(qū)土壤水分特征曲線預測提供可靠模型和技術支持。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)處于松嫩平原黑土區(qū)南部黑龍江省哈爾濱市(見圖1),位于東經125°42'~130°10',北緯44°04'~46°40',全年平均降水量569.1 mm,冬長夏短,屬于溫帶大陸性季風氣候。哈爾濱市土壤類型較多,以黑土為主,包括壤土、黏土、砂壤土和砂黏壤土,少量砂土分布江河兩岸。黑土土壤養(yǎng)分含量較豐富,是主要耕作土壤,土壤有機質含量高,適宜農作物生長。

    圖1 研究區(qū)地理位置及采樣點分布Fig.1 Geographical location and sampling point distribution in the study area

    1.2 供試土樣采集與測定

    選取壤土、砂黏壤土、砂壤土、黏土、砂土5種土壤質地作為供試土樣,采樣點分布如圖1所示。

    每個土樣采集0~30 cm土層,分別用100 cm3環(huán)刀取原狀土,自封袋取散土。風干散土土樣、過篩2 mm,采用MS2000激光粒度儀分析土壤顆粒,按國際制分類,分為砂粒(粒徑0.02~2 mm)、粉粒(粒徑0.002~0.02 mm)和黏粒(粒徑<0.002 mm)。采用高速離心機法測定土壤水分特征曲線,分別測定土樣1、3、5、10、30、50、100、300、400、500、700、900、1 000、1 200 kPa吸力下體積含水率。干容重采用環(huán)刀法測定,土壤樣品質地分類如表1所示。

    表1 土壤樣品質地分類Table1 Textureclassification of soil samples

    1.3 土壤水分特征曲線經驗模型簡介

    土壤水分特征曲線經驗模型根據參數數量分為三參數模型、四參數模型和五參數模型。研究表明,模型參數越多,實測數據擬合精度更高。

    選取4種土壤水分特征曲線經驗模型van Genuchten(VG)、Brooks and Corey(BC)、Modified Gardner(MG)、Log-Normal Distribution(LND)模型模擬黑土區(qū)土壤水分特征曲線。其中VG模型是van Genuchten在1981年提出應用最廣泛模型;BC模型是由Brooks和Corey較早提出的土壤水分特征曲線模型,該模型簡單,應用較廣;MG模型是由Matlan提出的指數型經驗模型,模型參數有明確物理意義,應用范圍較廣;LND模型是Kosugi提出的對數正態(tài)分布模型。4種模型均為四參數模型,具有簡單性和普遍使用性等特點。

    ①VG模型表達式如下:

    式中,θ—體積含水率;θs—土壤飽和體積含水率;θr—滯留土壤體積含水率;h—土壤水吸力(cm);α、n、m為表征模型形狀參數,其中α數值上等于進氣值倒數。

    ②BC模型表達式如下:

    式中,Se—飽和度。

    λ—土壤孔隙尺寸分布參數

    ③MG模型表達式如下:

    式中,b、t為模型參數。

    ④LND模型表達式如下:

    式中,q為模型參數。

    4個模型除θs和θr參數外,還包含兩個模型參數,分別用P1和P2表示,其中P1單位為kPa、P2為分形參數,如表2所示。

    表2 土壤水分特征曲線模型參數Table 2 Parameters of soil water characteristic curve models

    1.4 模型參數優(yōu)化方法

    1.41 粒子群優(yōu)化算法

    采用粒子群優(yōu)化算法(Particle Swarm Optimization,PSO)優(yōu)化模型參數,模仿鳥群捕食行為提出粒子群優(yōu)化算法。首先初始化1個種群N,把每個種群中個體假設d維搜索空間中1個粒子,第i個粒子位置和速度分別為Xi=(xi1,xi2,…,xid)和Vi=(vi1,vi2,…,vid),i=1,2,3…N,每次迭代中,粒子通過自我更新找出最優(yōu)解,即個體最優(yōu)解pbest,Pi=(pi1,pi2…,pid);gbest是整個種群N種中最優(yōu)解,兩個最優(yōu)解確定后,粒子速度更新和位置更新方程如下:

    其中,w為慣性權因子,c1和c2為學習因子,rand1和rand2為0~1均勻分布隨機數。

    1.4.2 參數模型構建

    為精確獲得4種模型模型參數,其中飽和含水率θs采用實測數據,以模型擬合體積含水率值與土壤實測體積含水率殘差平方和最小為目標函數,求解該非線性最優(yōu)化問題,即表達式如下:

    式中,θobs為土壤體積含水率實測值,θmod為模型擬合值,N為樣本數。

    1.5 模型精度評價指標

    采用均方根誤差(RMSE)和Pearson相關系數R作為定量評價指標。

    式中,X為體積含水率實測值,Xi為經驗模型擬合值。RSME評估模型整體誤差;R值反映實測值與模擬值趨勢。

    2 結果分析

    2.1 不同土壤質地土壤水分特征曲線變化趨勢

    土壤水分特征曲線反映土壤持水能力。通過離心機測定各吸力段下對應土壤體積含水率,繪制不同質地吸力段下體積含水率與土壤水基質吸力變化趨勢如圖2所示。

    圖2 不同質地土壤水分特征曲線Fig.2 Soil moisture characteristic curve of different soil textures

    由圖2可知,吸力值小于100 kPa時的低吸力段,土壤水分特征曲線急劇下降,變化率大,因為該階段土壤大孔隙排水占主導;當土壤吸力大于100 kPa時,隨著吸力增大,曲線變化率逐漸趨于0,吸力值大于100 kPa時,中高吸力段土壤大孔隙水分在重力作用下排空,僅小孔隙存留水分,毛管力對水分產生作用力,土壤具有較好持水能力,土壤水分特征曲線緩慢降低并趨于平穩(wěn)。在相同吸力條件下,體積含水率黏土>壤土>砂黏壤土>砂壤土>砂土,說明黏土持水能力最強,壤土次之,砂土最弱。結合表1不同質地顆粒組成百分含量可知,黏粒含量越高,持水能力越強,由于黏粒粒徑較小,黏粒含量越高,小孔隙越多,比表面積越大,對水作用力越強;而砂土砂粒含量較高,粒徑較大,比表面積較小,水吸附能力弱。因此降水后砂土大孔隙較多易飽和,形成徑流,不利于水土保持。

    2.2 土壤水分特征曲線擬合

    本文采用VG、BC、MG、LND 4種經驗模型擬合5種土壤質地實測土壤含水率數據曲線,基于粒子群優(yōu)化算法對4種模型參數優(yōu)化值如表3所示,繪制5種土壤類型土壤水分特征曲線實測值與采用4種模型擬合值關系如圖3所示。由圖3可知,VG模型對5種土壤擬合值與實測值接近,其中實測值在VG模型曲線上。滯留土壤含水率為0.034~0.092,符合實際土壤滯留含水率。5種土壤質地中黏土飽和含水率和滯留含水率相對較高;砂土飽和含水率和滯留含水率最低,由土壤性質決定。參數n值反映土壤水分特征曲線坡度,土壤水分特征曲線坡度越緩,n值越大,反之越小[23]。砂土坡度最陡,n值不是最小,故認為n為表征曲線形狀系數[24]。VG模型參數P1值為0.0375~0.323,則土壤進氣值為3.096~26.667,與飽和含水率對應進氣值相符。BC模型對砂土和砂黏壤土擬合趨勢較好,壤土擬合趨勢較差,可較好擬合5種土壤質地曲線。BC模型優(yōu)化參數中滯留含水率為1.20~5.68×10-4,可知BC模型低估5種土壤質地滯留含水率;模型參數P1值為0.2051~0.9001,土壤進氣值為1.1109~4.8756??芍狟C模型低估飽和含水率對應進氣值。MG模型對5種類型土壤擬合效果較好,實測值在MG模型曲線上。其中除砂土外,對其他4種土壤擬合值更接近實測值,砂土實測值擬合偏差較大。MG模型擬合參數滯留含水率為0.0124~0.1469,與土壤中滯留含水率相近。LND模型優(yōu)化參數中滯留含水率為0.00012~0.0736。可知LND模型低估壤土和砂黏壤土滯留含水率。LND模型對5種類型土壤擬合值與實測值偏差較小,表明LND模型可較好擬合5種土壤質地土壤水分特征曲線。

    表3 不同土壤質地模型參數Table 3 Model parameters of different soil textures

    2.3 土壤水分特征曲線模型比較及優(yōu)選

    比較Brooks-Corey(BC模型)、Van Genuchten(VG模型)、Modified Garder(MG模型)和Log-Normal Distribution(LND模型)適用性,對比分析模擬結果。

    由表4可知,VG模型對5種類型土壤擬合值與實測值相關系數R為0.996~0.998;BC模型R為0.9547~0.9869;MG模型R為0.994~0.998;LND模型R為0.9820~0.9960,可見4種模型擬合值與實測值變化趨勢相似。其中,除壤土外,MG模型對其他4種土壤擬合效果優(yōu)于VG、BC、LND模型;另外,VG模型對壤土模擬效果最佳。由圖3可知,VG、MG、LND 3種模型在低吸力段(0~100 kPa)壤土、砂黏壤土、砂壤土、黏土、砂土擬合值與實測值重合,而BC模型擬合值明顯偏離實測值;在中高吸力段(>100 kPa)VG、MG、BC、LND 4種模型擬合值均在實測值上方。其中BC模型擬合值與實測值偏差較大,而VG、MG、LND 3種模型擬合值略高于實測值,可知VG、MG、LND3種模型對5種類型土壤擬合效果優(yōu)于BC模型。

    由殘差平方和值可知,4種模型均可較好擬合不同土壤質地含水量與水吸力間關系,VG、MG模型殘差平方和值顯著明顯低于BC和LND模型,其中BC模型殘差平方和值最高,可知MG和VG模型對5種類型土壤擬合效果最佳,BC模型擬合效果最差。從RMSE值來看,MG模型對5種土壤質地RMSE值為0.0033~0.0107;VG模型RMSE變化為0.0039~0.0092;LND 模型 RMSE 變化為 0.0096~0.0114;BC模型RMSE變化為0.0140~0.0250。結果表明,MG和VG模型RSME值明顯小于LND和BC模型,其中BC模型RMSE值比其他3種模型RMSE值高一個數量級,BC模型對5種土壤類型模擬效果最差。綜合比較擬合誤差分析可知,VG和MG模型對黑土區(qū)土壤模擬精度最高。

    表4 4種模型模擬土壤水分特征曲線R、RSME和殘差平方和Table4 R,RSMEand sum of squared residualsof soil water characteristic curvecompared to simulationsfor four models

    3 討 論

    本研究利用4種土壤水分特征曲線經驗模型擬合5種不同質地土壤曲線,4種模型擬合結果較佳,其中BC模型在擬合過程中低估土壤滯留含水率和土壤進氣值,由于實測飽和含水率存在誤差,導致出現低估土壤滯留含水率和土壤進氣值現象。模型模擬結果中高吸力段擬合大于實測值,由于高速離心機測定過程中土壤結構被破壞,容重改變,導致誤差[25]。

    對比分析4種經驗模型模擬結果可知,VG和MG模型模擬5種不同質地土壤效果最好。丁新原等比較VG和其他模型的土壤水分特征曲線模擬結果,VG模型模擬精度最高,與本研究結果一致[20,23,26],因VG和MG模型中各參數均有明確物理意義[19,27]。Matlan等利用4種模型模擬分析黏土、砂土、粉土3種土壤質地,表明BC模型對實測值擬合效果最差,與本文研究結果一致[19]。其他4種土壤質地,BC模型對砂土模擬效果最優(yōu),表明BC模型更適合模擬粗質地土壤水分特征曲線[28]。

    4 結 論

    黑土區(qū)不同質地土壤土壤體積含水率因土壤吸力變化持水能力存在差異,在相同吸力下,土壤黏粒含量越高,土壤體積含水率越高。在低吸力段,大孔隙排水,砂粒含量多,曲線變化率越快,釋水能力越強;在中高吸力段,土壤小孔隙比例較高,土壤對水吸附作用較強,因此土壤水分特征曲線逐漸趨于平穩(wěn)。

    本文利用粒子群優(yōu)化算法優(yōu)化VG、BC、MG、LND 4種土壤水分特征曲線模型參數。比較4種經驗模型擬合值與實測值,VG、MG模型對黑土區(qū)土壤水分特征曲線模擬精度較高。預測黑土區(qū)土壤水分特征曲線推薦使用VG和MG模型。

    猜你喜歡
    壤土吸力實測值
    深水大型吸力錨測試技術
    ROV在海上吸力樁安裝場景的應用及安裝精度和風險控制
    化工管理(2022年11期)2022-06-03 07:08:24
    ±800kV直流輸電工程合成電場夏季實測值與預測值比對分析
    深水吸力樁施工技術研究
    土壤質地及砧木影響蘋果根際微生物功能多樣性及其碳源利用
    常用高溫軸承鋼的高溫硬度實測值與計算值的對比分析
    哈爾濱軸承(2020年1期)2020-11-03 09:16:22
    左家林場核桃良種基地選址調查報告
    CONTENTS
    時代經貿(2019年3期)2019-11-29 11:15:16
    市售純牛奶和巴氏殺菌乳營養(yǎng)成分分析
    中國奶牛(2019年10期)2019-10-28 06:23:36
    一種基于實測值理論計算的導航臺電磁干擾分析方法
    電子制作(2018年23期)2018-12-26 01:01:22
    亚洲久久久国产精品| 夜夜骑夜夜射夜夜干| 免费看十八禁软件| 美女高潮到喷水免费观看| 视频区欧美日本亚洲| 亚洲国产毛片av蜜桃av| 黄色视频在线播放观看不卡| 色视频在线一区二区三区| 国产精品亚洲av一区麻豆| 丝袜人妻中文字幕| 久久精品人人爽人人爽视色| 久久狼人影院| 日本五十路高清| 不卡av一区二区三区| 亚洲情色 制服丝袜| 国产麻豆69| 视频区欧美日本亚洲| 久久久欧美国产精品| 欧美成人午夜精品| 国产成人精品久久二区二区免费| av福利片在线| 成年美女黄网站色视频大全免费| 欧美老熟妇乱子伦牲交| 美女午夜性视频免费| 亚洲国产欧美日韩在线播放| 在现免费观看毛片| 黄色a级毛片大全视频| 超色免费av| 亚洲国产日韩一区二区| 国产淫语在线视频| 亚洲 国产 在线| 一本综合久久免费| 一级毛片电影观看| 精品熟女少妇八av免费久了| 另类亚洲欧美激情| 免费在线观看视频国产中文字幕亚洲 | 日韩伦理黄色片| 99国产综合亚洲精品| 人妻 亚洲 视频| 亚洲专区中文字幕在线| 亚洲欧美激情在线| 女人爽到高潮嗷嗷叫在线视频| 久久ye,这里只有精品| 777米奇影视久久| 操出白浆在线播放| 国产黄色免费在线视频| 国产成人精品无人区| 欧美av亚洲av综合av国产av| 精品亚洲乱码少妇综合久久| 999精品在线视频| 一边摸一边做爽爽视频免费| 少妇猛男粗大的猛烈进出视频| 亚洲av电影在线观看一区二区三区| 日韩大码丰满熟妇| 你懂的网址亚洲精品在线观看| 国产精品九九99| 欧美av亚洲av综合av国产av| 亚洲专区中文字幕在线| 亚洲欧美日韩另类电影网站| 在线观看国产h片| 爱豆传媒免费全集在线观看| 中文字幕人妻熟女乱码| 丝袜人妻中文字幕| 欧美日韩成人在线一区二区| 亚洲精品自拍成人| 国产精品 国内视频| 五月开心婷婷网| 一二三四社区在线视频社区8| 日韩 欧美 亚洲 中文字幕| 少妇粗大呻吟视频| 国产精品国产三级国产专区5o| 久久 成人 亚洲| 色视频在线一区二区三区| 午夜激情久久久久久久| 欧美人与性动交α欧美精品济南到| 国产精品av久久久久免费| 欧美日本中文国产一区发布| 日韩 欧美 亚洲 中文字幕| 999精品在线视频| 狂野欧美激情性bbbbbb| 亚洲av男天堂| 91麻豆精品激情在线观看国产 | 91精品三级在线观看| 国产国语露脸激情在线看| 久久久欧美国产精品| 飞空精品影院首页| 下体分泌物呈黄色| 久久综合国产亚洲精品| 成人手机av| 中文字幕另类日韩欧美亚洲嫩草| www.av在线官网国产| 国产真人三级小视频在线观看| 久久久久网色| 真人做人爱边吃奶动态| 国产精品亚洲av一区麻豆| 丝袜人妻中文字幕| 国产黄频视频在线观看| 欧美精品一区二区大全| 满18在线观看网站| 色网站视频免费| 性色av乱码一区二区三区2| 国产亚洲欧美在线一区二区| 欧美在线黄色| 18禁观看日本| 免费少妇av软件| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲av日韩在线播放| 久久精品人人爽人人爽视色| 看十八女毛片水多多多| 国产免费现黄频在线看| 国产精品.久久久| 国产欧美日韩精品亚洲av| 欧美老熟妇乱子伦牲交| 久久精品熟女亚洲av麻豆精品| 黑人巨大精品欧美一区二区蜜桃| 热re99久久精品国产66热6| 亚洲中文字幕日韩| 午夜两性在线视频| 一级a爱视频在线免费观看| 曰老女人黄片| 男人爽女人下面视频在线观看| 国产片特级美女逼逼视频| 欧美日韩国产mv在线观看视频| 亚洲中文日韩欧美视频| 国产成人a∨麻豆精品| 国产人伦9x9x在线观看| 欧美少妇被猛烈插入视频| 国产精品 欧美亚洲| 黑人猛操日本美女一级片| 高清av免费在线| 国产国语露脸激情在线看| 精品久久久久久电影网| 色网站视频免费| 91成人精品电影| 国产99久久九九免费精品| 亚洲五月婷婷丁香| 老汉色av国产亚洲站长工具| 黑人巨大精品欧美一区二区蜜桃| 少妇猛男粗大的猛烈进出视频| 欧美久久黑人一区二区| 欧美少妇被猛烈插入视频| 欧美亚洲 丝袜 人妻 在线| 国产欧美日韩一区二区三 | 国产无遮挡羞羞视频在线观看| 91九色精品人成在线观看| 丰满少妇做爰视频| 精品卡一卡二卡四卡免费| 国产av国产精品国产| 一级毛片电影观看| 脱女人内裤的视频| 真人做人爱边吃奶动态| 一本大道久久a久久精品| 97人妻天天添夜夜摸| 国产亚洲av高清不卡| 国产精品一二三区在线看| 欧美成人精品欧美一级黄| 久久久久久人人人人人| 老司机靠b影院| 97在线人人人人妻| av在线app专区| 99国产精品免费福利视频| 久久久亚洲精品成人影院| 欧美日韩一级在线毛片| 18禁黄网站禁片午夜丰满| 亚洲熟女精品中文字幕| 最黄视频免费看| 欧美日韩一级在线毛片| 18禁裸乳无遮挡动漫免费视频| 国产精品亚洲av一区麻豆| xxx大片免费视频| 亚洲精品国产一区二区精华液| 丰满人妻熟妇乱又伦精品不卡| 在线观看免费午夜福利视频| 精品亚洲乱码少妇综合久久| 国产女主播在线喷水免费视频网站| 巨乳人妻的诱惑在线观看| 亚洲精品中文字幕在线视频| 久久国产精品人妻蜜桃| 欧美精品av麻豆av| 在线观看国产h片| 亚洲欧美日韩高清在线视频 | 国产真人三级小视频在线观看| 久久亚洲精品不卡| 亚洲成av片中文字幕在线观看| 老司机影院毛片| 国产伦人伦偷精品视频| 国产伦人伦偷精品视频| 在线观看免费日韩欧美大片| 亚洲欧美成人综合另类久久久| 一级毛片 在线播放| 激情五月婷婷亚洲| 王馨瑶露胸无遮挡在线观看| 亚洲欧美清纯卡通| 亚洲精品久久成人aⅴ小说| 天天影视国产精品| 极品少妇高潮喷水抽搐| 日韩电影二区| 一本一本久久a久久精品综合妖精| 91精品三级在线观看| 看十八女毛片水多多多| 你懂的网址亚洲精品在线观看| 亚洲av欧美aⅴ国产| 国产精品成人在线| 国产97色在线日韩免费| 国产精品欧美亚洲77777| 搡老乐熟女国产| 成人影院久久| 免费一级毛片在线播放高清视频 | 婷婷色综合www| 国产精品免费视频内射| 欧美 日韩 精品 国产| 精品亚洲乱码少妇综合久久| 成人午夜精彩视频在线观看| 人人妻人人爽人人添夜夜欢视频| 女性生殖器流出的白浆| 亚洲男人天堂网一区| 日韩熟女老妇一区二区性免费视频| 老司机影院毛片| 男女免费视频国产| 国产精品 国内视频| 亚洲欧洲国产日韩| 国产人伦9x9x在线观看| 中文字幕高清在线视频| 国产在视频线精品| 在现免费观看毛片| 在线观看国产h片| 国产成人精品久久二区二区免费| 中文字幕精品免费在线观看视频| 9色porny在线观看| 欧美久久黑人一区二区| 日韩中文字幕视频在线看片| 亚洲精品国产色婷婷电影| 亚洲中文日韩欧美视频| av电影中文网址| 亚洲欧美精品综合一区二区三区| 嫩草影视91久久| videos熟女内射| 午夜福利一区二区在线看| 深夜精品福利| 久久人妻福利社区极品人妻图片 | 超碰成人久久| 亚洲成人免费电影在线观看 | 亚洲人成电影免费在线| 精品一区二区三区四区五区乱码 | 50天的宝宝边吃奶边哭怎么回事| 一级毛片女人18水好多 | 女人高潮潮喷娇喘18禁视频| 日本欧美视频一区| 国产午夜精品一二区理论片| 欧美日本中文国产一区发布| 90打野战视频偷拍视频| 午夜福利免费观看在线| 最新的欧美精品一区二区| 18禁黄网站禁片午夜丰满| 一级毛片电影观看| 91九色精品人成在线观看| 午夜激情av网站| 亚洲精品一区蜜桃| 日韩欧美一区视频在线观看| 午夜福利,免费看| 精品国产一区二区三区久久久樱花| 只有这里有精品99| 中文字幕色久视频| xxxhd国产人妻xxx| 国产免费一区二区三区四区乱码| 亚洲欧美日韩另类电影网站| 97在线人人人人妻| 欧美激情 高清一区二区三区| 精品免费久久久久久久清纯 | 日韩av不卡免费在线播放| 成人手机av| 成人国产av品久久久| 亚洲精品国产av成人精品| 国产精品免费视频内射| 天天躁夜夜躁狠狠躁躁| 国产精品国产av在线观看| 国产伦人伦偷精品视频| 在线观看免费午夜福利视频| 男女国产视频网站| 欧美乱码精品一区二区三区| 黄网站色视频无遮挡免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 九草在线视频观看| 亚洲精品日韩在线中文字幕| 一本色道久久久久久精品综合| 免费在线观看日本一区| 老司机影院成人| 新久久久久国产一级毛片| 高清视频免费观看一区二区| 99国产综合亚洲精品| 日本av免费视频播放| 国产亚洲欧美精品永久| 亚洲五月色婷婷综合| 亚洲欧美日韩另类电影网站| 国产在线观看jvid| 人体艺术视频欧美日本| 青草久久国产| 亚洲欧洲国产日韩| 两个人看的免费小视频| 亚洲午夜精品一区,二区,三区| 十八禁高潮呻吟视频| 国产一区二区三区综合在线观看| 国产成人影院久久av| 日本猛色少妇xxxxx猛交久久| 热re99久久国产66热| 亚洲精品国产区一区二| 欧美老熟妇乱子伦牲交| 亚洲av国产av综合av卡| 丰满人妻熟妇乱又伦精品不卡| 老司机影院成人| 波多野结衣av一区二区av| 最近中文字幕2019免费版| 国产熟女午夜一区二区三区| 久久午夜综合久久蜜桃| 成年人午夜在线观看视频| 亚洲专区中文字幕在线| 九草在线视频观看| 欧美日韩一级在线毛片| 青春草亚洲视频在线观看| 精品少妇一区二区三区视频日本电影| 日韩视频在线欧美| 国产免费福利视频在线观看| 性少妇av在线| 精品高清国产在线一区| 国产一级毛片在线| 亚洲人成网站在线观看播放| 高清黄色对白视频在线免费看| 一级毛片 在线播放| 亚洲成人手机| 在线亚洲精品国产二区图片欧美| 精品少妇黑人巨大在线播放| 久久久亚洲精品成人影院| 国产激情久久老熟女| 我要看黄色一级片免费的| 一级黄片播放器| 午夜福利乱码中文字幕| 成年美女黄网站色视频大全免费| 亚洲av男天堂| 亚洲精品日本国产第一区| 精品国产一区二区三区四区第35| 你懂的网址亚洲精品在线观看| 国产麻豆69| 国产在视频线精品| 人人妻,人人澡人人爽秒播 | 大话2 男鬼变身卡| 日本vs欧美在线观看视频| 99久久99久久久精品蜜桃| 黑人巨大精品欧美一区二区蜜桃| 免费看不卡的av| 久久99精品国语久久久| 人人妻人人添人人爽欧美一区卜| 黄色 视频免费看| 啦啦啦 在线观看视频| 国产国语露脸激情在线看| 国产免费一区二区三区四区乱码| 精品亚洲成a人片在线观看| 国产精品av久久久久免费| 精品久久久久久电影网| av电影中文网址| 国产精品香港三级国产av潘金莲 | 大片免费播放器 马上看| 悠悠久久av| 亚洲av日韩精品久久久久久密 | 国产99久久九九免费精品| 另类精品久久| 国产精品.久久久| 亚洲欧美成人综合另类久久久| 美女福利国产在线| 丝瓜视频免费看黄片| 亚洲九九香蕉| 国产高清videossex| 国产精品免费视频内射| 国产主播在线观看一区二区 | 99九九在线精品视频| 午夜视频精品福利| 日韩熟女老妇一区二区性免费视频| 人人妻人人澡人人看| 欧美乱码精品一区二区三区| 亚洲欧美色中文字幕在线| 国产亚洲精品久久久久5区| 免费av中文字幕在线| 亚洲三区欧美一区| 精品人妻熟女毛片av久久网站| 国产成人影院久久av| 青青草视频在线视频观看| 午夜视频精品福利| 夫妻性生交免费视频一级片| 国产成人欧美在线观看 | 91精品伊人久久大香线蕉| 涩涩av久久男人的天堂| 国产伦人伦偷精品视频| 一级片'在线观看视频| 亚洲中文字幕日韩| 欧美日本中文国产一区发布| 制服诱惑二区| 男男h啪啪无遮挡| 无遮挡黄片免费观看| 亚洲国产精品一区二区三区在线| 啦啦啦在线观看免费高清www| 国产老妇伦熟女老妇高清| 久久av网站| 最黄视频免费看| 人妻人人澡人人爽人人| 久久亚洲国产成人精品v| 中文字幕人妻熟女乱码| 黄色毛片三级朝国网站| 国产野战对白在线观看| 亚洲欧美精品自产自拍| 亚洲欧洲日产国产| 丰满饥渴人妻一区二区三| 中文字幕人妻丝袜一区二区| 国产精品秋霞免费鲁丝片| 日本欧美国产在线视频| 人妻一区二区av| 少妇人妻久久综合中文| 国产精品久久久久久精品电影小说| 亚洲,欧美精品.| 久久久精品国产亚洲av高清涩受| av福利片在线| 中文字幕av电影在线播放| 丁香六月天网| 人成视频在线观看免费观看| 欧美精品av麻豆av| 91字幕亚洲| 国产精品 国内视频| 免费一级毛片在线播放高清视频 | 日韩电影二区| 青草久久国产| 无限看片的www在线观看| 赤兔流量卡办理| 国产极品粉嫩免费观看在线| 成年人午夜在线观看视频| www.999成人在线观看| 色婷婷av一区二区三区视频| 又紧又爽又黄一区二区| 另类亚洲欧美激情| 免费在线观看黄色视频的| 婷婷色综合www| 热99久久久久精品小说推荐| av片东京热男人的天堂| 国产日韩欧美视频二区| 亚洲av日韩精品久久久久久密 | 亚洲精品国产av成人精品| 精品人妻在线不人妻| 精品国产乱码久久久久久小说| 看免费av毛片| 中文字幕人妻丝袜制服| 人人妻人人澡人人看| 中文乱码字字幕精品一区二区三区| 在线观看免费视频网站a站| 操出白浆在线播放| 中文乱码字字幕精品一区二区三区| 精品少妇黑人巨大在线播放| 18禁裸乳无遮挡动漫免费视频| 韩国高清视频一区二区三区| 日韩熟女老妇一区二区性免费视频| 欧美另类一区| 久久人人97超碰香蕉20202| 99国产精品一区二区蜜桃av | 啦啦啦视频在线资源免费观看| 日本91视频免费播放| 只有这里有精品99| 亚洲中文日韩欧美视频| 午夜免费观看性视频| 欧美日韩福利视频一区二区| 国产精品99久久99久久久不卡| videosex国产| 成人18禁高潮啪啪吃奶动态图| 国产精品麻豆人妻色哟哟久久| 两性夫妻黄色片| 一本久久精品| 老汉色av国产亚洲站长工具| 国产成人欧美在线观看 | 波野结衣二区三区在线| 欧美激情 高清一区二区三区| 国产91精品成人一区二区三区 | 精品少妇久久久久久888优播| 亚洲av片天天在线观看| 伊人久久大香线蕉亚洲五| 最近最新中文字幕大全免费视频 | 久久亚洲精品不卡| 国产一区二区在线观看av| 一区二区三区激情视频| 在线看a的网站| 好男人视频免费观看在线| 18禁国产床啪视频网站| 色综合欧美亚洲国产小说| 丝瓜视频免费看黄片| 国产精品国产三级专区第一集| 2021少妇久久久久久久久久久| 国产欧美亚洲国产| www.精华液| 亚洲av在线观看美女高潮| 日韩大码丰满熟妇| 女性被躁到高潮视频| 人成视频在线观看免费观看| www日本在线高清视频| 香蕉丝袜av| 国精品久久久久久国模美| 男女床上黄色一级片免费看| 国产成人啪精品午夜网站| 国产野战对白在线观看| 99精品久久久久人妻精品| 国产亚洲av高清不卡| 一级毛片黄色毛片免费观看视频| av在线播放精品| 人妻一区二区av| 欧美国产精品一级二级三级| 涩涩av久久男人的天堂| 亚洲色图 男人天堂 中文字幕| 亚洲欧美精品综合一区二区三区| 黄频高清免费视频| 大片免费播放器 马上看| 亚洲精品美女久久久久99蜜臀 | 国产淫语在线视频| 一本—道久久a久久精品蜜桃钙片| 国产欧美日韩一区二区三区在线| 老司机影院毛片| 久久 成人 亚洲| 高清视频免费观看一区二区| 欧美精品av麻豆av| avwww免费| 亚洲国产欧美在线一区| 精品人妻一区二区三区麻豆| 久久精品久久久久久噜噜老黄| 亚洲精品国产区一区二| 男女边摸边吃奶| 高清欧美精品videossex| videosex国产| 午夜免费鲁丝| 熟女av电影| 亚洲情色 制服丝袜| 午夜福利一区二区在线看| 国产男女内射视频| 亚洲成人免费电影在线观看 | svipshipincom国产片| 免费少妇av软件| 青草久久国产| 人妻人人澡人人爽人人| 亚洲成人国产一区在线观看 | tube8黄色片| 国产极品粉嫩免费观看在线| 国产亚洲av高清不卡| 亚洲国产欧美一区二区综合| 夜夜骑夜夜射夜夜干| 国产男女超爽视频在线观看| 亚洲色图综合在线观看| 欧美黄色片欧美黄色片| 久久久久精品国产欧美久久久 | 国产无遮挡羞羞视频在线观看| 亚洲人成电影免费在线| 久久人人97超碰香蕉20202| 日本一区二区免费在线视频| 老汉色∧v一级毛片| 丝袜人妻中文字幕| 老司机在亚洲福利影院| 成年美女黄网站色视频大全免费| 好男人视频免费观看在线| 色94色欧美一区二区| 在线观看免费高清a一片| 国产三级黄色录像| 人人妻人人添人人爽欧美一区卜| 久热这里只有精品99| 亚洲综合色网址| 国产精品久久久av美女十八| 亚洲av男天堂| 丁香六月天网| 一本色道久久久久久精品综合| 新久久久久国产一级毛片| 国产成人91sexporn| 国产爽快片一区二区三区| 日本午夜av视频| 啦啦啦啦在线视频资源| 国产精品偷伦视频观看了| 亚洲欧美日韩另类电影网站| 免费高清在线观看视频在线观看| 久久人人爽人人片av| 嫁个100分男人电影在线观看 | 色94色欧美一区二区| 一区二区三区四区激情视频| 亚洲av日韩精品久久久久久密 | 久久久久久免费高清国产稀缺| 汤姆久久久久久久影院中文字幕| 亚洲五月色婷婷综合| 国产成人影院久久av| 日韩一区二区三区影片| 欧美亚洲 丝袜 人妻 在线| 啦啦啦啦在线视频资源| 欧美日韩视频精品一区| 飞空精品影院首页| 精品国产乱码久久久久久小说| 午夜精品国产一区二区电影| 免费在线观看黄色视频的| 国产在线一区二区三区精| 亚洲国产av新网站| 新久久久久国产一级毛片| 国产精品亚洲av一区麻豆| 国产主播在线观看一区二区 | 免费看不卡的av| 国产成人免费无遮挡视频| 中文字幕人妻丝袜制服| 又紧又爽又黄一区二区| 亚洲熟女精品中文字幕| 精品亚洲成a人片在线观看| av一本久久久久| 国产亚洲欧美在线一区二区| 欧美人与性动交α欧美软件| 男女边摸边吃奶| 国产成人欧美| 国产成人免费观看mmmm| 男女无遮挡免费网站观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲人成77777在线视频| 夜夜骑夜夜射夜夜干| 精品福利永久在线观看|