亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于濕潤(rùn)區(qū)分層假定對(duì)Green-Ampt模型的改進(jìn)

        2018-10-20 06:31:44李泳霖王仰仁通信作者孫小平鄭志偉
        關(guān)鍵詞:非飽和濕潤(rùn)幅度

        李泳霖,王仰仁,通信作者,孫小平,鄭志偉

        ?

        基于濕潤(rùn)區(qū)分層假定對(duì)Green-Ampt模型的改進(jìn)

        李泳霖1,2,王仰仁1,2,通信作者,孫小平3,鄭志偉1,2

        (1. 天津農(nóng)學(xué)院 水利工程學(xué)院,天津 300384;2. 天津市節(jié)水灌溉技術(shù)與裝備校企協(xié)同創(chuàng)新實(shí)驗(yàn)室,天津 300384; 3. 山西省水利水電科學(xué)研究院,太原 030012)

        假定入滲過(guò)程分為前期、后期兩個(gè)階段,入滲前期土壤飽和區(qū)與非飽和區(qū)厚度之比為一常數(shù)2,非飽和區(qū)厚度隨著入滲的進(jìn)行而逐漸增大;入滲后期非飽和區(qū)的厚度不再隨時(shí)間變化,為一定值,只有飽和區(qū)厚度隨時(shí)間變化,基于此對(duì)Green-Ampt模型進(jìn)行了改進(jìn)。用傳統(tǒng)Green-Ampt模型和改進(jìn)模型對(duì)6組入滲測(cè)試資料進(jìn)行擬合,傳統(tǒng)模型擬合的相關(guān)系數(shù)為0.982(6組數(shù)據(jù)的平均值),相比傳統(tǒng)模型,改進(jìn)模型的擬合精度明顯提高,相關(guān)系數(shù)提高了0.016,達(dá)到0.998(對(duì)應(yīng)=15);其中兩組數(shù)據(jù)的效果尤為明顯,相關(guān)系數(shù)分別提高了0.07和0.04,均達(dá)到了0.99以上。改進(jìn)模型與傳統(tǒng)模型同樣具有較強(qiáng)的物理意義,但與傳統(tǒng)模型相比,改進(jìn)模型擬合得到的穩(wěn)定入滲率更接近于測(cè)試值,具有更廣泛的適用性。

        Green-Ampt入滲模型;累計(jì)入滲量;非飽和區(qū);入滲鋒;入滲時(shí)間

        Green-Ampt入滲模型是Green和Ampt在1911年提出的基于毛管理論的積水入滲模型[1],該模型的基本假定是,入滲時(shí)存在明確的水平濕潤(rùn)鋒面,將濕潤(rùn)和未濕潤(rùn)區(qū)域明顯分開(kāi),濕潤(rùn)區(qū)為飽和含水率,濕潤(rùn)鋒前即為初始含水率,即土壤水分剖面為階梯狀分布,故該模型又稱(chēng)活塞(或打氣筒)模型[2]。由于該模型具有較強(qiáng)的物理意義、參數(shù)少、求解簡(jiǎn)單,從而被廣泛應(yīng)用。為適用不同積水條件[3-5]、不同土質(zhì)土層條件、不同地形條件,研究者對(duì)模型進(jìn)行了不同程度的改進(jìn)。其中,劉姍姍等[6]對(duì)模型關(guān)鍵參數(shù)簡(jiǎn)化求解及與土壤物理參數(shù)的關(guān)系進(jìn)行了研究;國(guó)外學(xué)者提出均質(zhì)土壤在入滲期間從上至下依次可分為飽和層、半飽和層、濕潤(rùn)層和干土層[7]。Ma等[8]提出,濕潤(rùn)鋒后土壤導(dǎo)水率近似用飽和導(dǎo)水率和飽和系數(shù)之積代替;Bouwer[9-10]建議濕潤(rùn)鋒后導(dǎo)水率為飽和導(dǎo)水率的0.5倍;毛麗麗等[11]通過(guò)分段擬合直線(xiàn)方法描述含水量的分布;王文焰等[12]通過(guò)設(shè)計(jì)入滲試驗(yàn)測(cè)得濕潤(rùn)層厚度約為濕潤(rùn)峰后土層厚度的一半,并在此基礎(chǔ)上假定濕潤(rùn)層內(nèi)含水率分布為橢圓曲線(xiàn),對(duì)模型進(jìn)行了相應(yīng)改進(jìn);彭振陽(yáng)等[13]采用 Richards 模型,分析了不同入滲條件下各層厚度以及各層內(nèi)含水率和導(dǎo)水率的變化規(guī)律,并據(jù)此規(guī)律改進(jìn) Green-Ampt 模型。本文針對(duì)傳統(tǒng)模型擬合參數(shù)s過(guò)小甚至為負(fù)值的情況[14],對(duì)傳統(tǒng)模型進(jìn)行改進(jìn),將入滲階段分為前期和后期,將土壤由上至下依次分為飽和區(qū)、非飽和區(qū)和未濕潤(rùn)區(qū),并結(jié)合試驗(yàn)數(shù)據(jù)進(jìn)行驗(yàn)證。結(jié)果表明,改進(jìn)模型的計(jì)算精度顯著提高,與傳統(tǒng)模型相比,穩(wěn)定入滲率更接近實(shí)際情況。

        1 基本假定與模型建立

        1.1 傳統(tǒng)Green-Ampt模型

        Green-Ampt模型研究的是初始干燥的土壤在薄層積水時(shí)的入滲問(wèn)題,基本假定是入滲時(shí)存在明顯的水平濕潤(rùn)鋒面,將濕潤(rùn)和未濕潤(rùn)的區(qū)域截然分開(kāi),濕潤(rùn)區(qū)水分為飽和含水率θ,濕潤(rùn)鋒前為初始含水率θ。如圖1。

        地表積水深度記為,不隨時(shí)間改變,濕潤(rùn)鋒的位置為z,隨時(shí)間前移,濕潤(rùn)鋒處的土壤水吸力為s,該模型的主要任務(wù)是得出入滲量,入滲率及濕潤(rùn)鋒面位置與入滲時(shí)間的關(guān)系。

        由達(dá)西定律可求出入滲水量由地表進(jìn)入土壤的通量,即入滲率

        式中:K為飽和導(dǎo)水率,也稱(chēng)為穩(wěn)定入滲率。

        根據(jù)模型假定,由水量平衡原理,可得出累計(jì)入滲量和濕潤(rùn)鋒z的關(guān)系:

        由入滲率與入滲量的關(guān)系可得:

        聯(lián)立式(1)、式(3)可得:

        對(duì)式(4)整理、積分,并利用時(shí)z=0,有:

        式(1)、式(2)、式(5)便是Green-Ampt模型的主要入滲關(guān)系式。由(5)式可得z的關(guān)系,代入式(2)可得的關(guān)系:

        1.2 Green-Ampt模型的改進(jìn)

        1.2.1 基本假定

        假定入滲過(guò)程中土壤由上到下依次為飽和區(qū)、非飽和區(qū)和未濕潤(rùn)區(qū),并依據(jù)非飽和區(qū)厚度的變化將整個(gè)入滲過(guò)程分為前期、后期兩個(gè)階段。

        入滲前期(圖2-a),非飽和區(qū)厚度隨著入滲時(shí)間的增加而增加,并假定非飽和區(qū)厚度與飽和區(qū)厚度之比為一定值2,即:

        式中:為一常數(shù),Δ為非飽和區(qū)厚度,z為飽和區(qū)厚度。

        入滲后期(圖2-b),非飽和區(qū)厚度不再隨入滲時(shí)間的增加而增加,非飽和區(qū)厚度Δ為一定值Δz

        入滲前期結(jié)束時(shí),非飽和區(qū)厚度達(dá)到最大值Δz,此時(shí)飽和區(qū)與非飽和區(qū)依然有式(7)的關(guān)系,記此時(shí)的飽和區(qū)厚度為Δz。

        圖2 改進(jìn)Green-Ampt模型假設(shè)示意圖

        1.2.2 模型建立

        入滲前期,由達(dá)西定律和水量平衡原理,分別可得累計(jì)入滲率與累計(jì)入滲量:

        由累計(jì)入滲量與入滲率的關(guān)系,并結(jié)合式(7)有:

        聯(lián)立式(8)可得:

        對(duì)式(11)積分,并利用0時(shí)z=0,則有:

        將式(7)、式(9)代入式(12),可得~的關(guān)系:

        或者:

        入滲前期結(jié)束時(shí),對(duì)應(yīng)總?cè)霛B量為Δ,相應(yīng)的時(shí)間記為t,則有:

        進(jìn)入入滲后期,>t,累計(jì)入滲量與入滲率的關(guān)系:

        對(duì)式(15)積分,并利用z=Δz時(shí)=t,可得:

        式中:Δzt時(shí)刻對(duì)應(yīng)的飽和區(qū)厚度,式(9)中,當(dāng)=Δ時(shí),s=Δs,然后聯(lián)合式(7)、式(9)和式(16),可得:

        則改進(jìn)后的入滲時(shí)間與總?cè)霛B量的關(guān)系式為:

        式(18)中參數(shù)、、ΔI、K需要通過(guò)擬合來(lái)確定,確定了4個(gè)參數(shù)的值后可由式(14)得的值。

        2 模型檢驗(yàn)

        為分析改進(jìn)效果,使用雙環(huán)入滲儀,在天津農(nóng)學(xué)院西校區(qū)農(nóng)田水循環(huán)試驗(yàn)基地進(jìn)行入滲試驗(yàn),實(shí)驗(yàn)基地土壤類(lèi)型為中壤土,在試驗(yàn)過(guò)程中人工加水,積水厚度保持在5 cm,共獲得6組試驗(yàn)數(shù)據(jù),用于測(cè)試模型。其中,4組數(shù)據(jù)(第1組,第2組,第4組,第5組)由小麥(2016年3月20日)地測(cè)得,兩組數(shù)據(jù)(第3組,第6組)由玉米(2017年5月18日)地測(cè)得。

        2.1 模型參數(shù)的擬合

        分別采用傳統(tǒng)模型和改進(jìn)模型對(duì)6組入滲試驗(yàn)數(shù)據(jù)進(jìn)行擬合。傳統(tǒng)模型使用式(6)進(jìn)行擬合,擬合目標(biāo)是使累計(jì)入滲時(shí)間的觀測(cè)值和模擬值的誤差平方和最小。首先對(duì)參數(shù)、K設(shè)置相應(yīng)的初始值[15],然后對(duì)參數(shù)進(jìn)行規(guī)劃求解,使誤差平方和最小,計(jì)算出傳統(tǒng)Green-Ampt模型累計(jì)入滲時(shí)間的模擬值與觀測(cè)值的相關(guān)系數(shù),以及相應(yīng)的參數(shù)、K值。根據(jù)實(shí)測(cè)資料,由各時(shí)段的累計(jì)入滲量除以相應(yīng)的時(shí)間,得到各時(shí)段的入滲率,選取各組入滲數(shù)據(jù)最后一個(gè)時(shí)段的入滲率,近似作為測(cè)試的穩(wěn)定入滲率K,然后,計(jì)算測(cè)點(diǎn)間參數(shù)、K的變差系數(shù),見(jiàn)表1。

        改進(jìn)模型的擬合過(guò)程與傳統(tǒng)模型類(lèi)似,不同的是,對(duì)常數(shù)從小到大設(shè)置10個(gè)不同的數(shù)(1,2,3,4,5,7,10,15,20,25),對(duì)于不同的分別進(jìn)行擬合,并計(jì)算相關(guān)系數(shù)以及各參數(shù)的變差系數(shù),結(jié)果見(jiàn)表2。

        表1 傳統(tǒng)模型的參數(shù)值

        表2 不同β時(shí)各組參數(shù)值

        2.2 擬合結(jié)果對(duì)比分析

        2.2.1值的確定

        改進(jìn)模型的6組相關(guān)系數(shù)均隨的增大而增加,6組相關(guān)系數(shù)平均值從0.993 0增大到0.998 8,但各組相關(guān)系數(shù)增大的幅度不同(表2)。其中,第2組相關(guān)系數(shù)值隨增大而增加的幅度最明顯,從0.97提高到0.99;第4組數(shù)據(jù)次之;第1、3、6組相關(guān)系數(shù)隨增大而增加的幅度不明顯;第5組數(shù)據(jù),當(dāng)較小時(shí),相關(guān)系數(shù)有增大的趨勢(shì),但當(dāng)>15時(shí)相關(guān)系數(shù)出現(xiàn)減小的趨勢(shì)。綜合來(lái)看,的取值越大,相關(guān)系數(shù)越高,但不宜過(guò)大。

        Δ表示當(dāng)非飽和區(qū)達(dá)到最大厚度時(shí),對(duì)應(yīng)滲入的總水量,它與土壤初始含水率和土壤類(lèi)型有很大關(guān)系。整體來(lái)看(圖3),隨著的增大,參數(shù)Δ整體呈增大的趨勢(shì),但各組參數(shù)Δ的擬合值范圍不同。當(dāng)較小時(shí),參數(shù)Δ隨增大的變化幅度較明顯;當(dāng)≥10時(shí),參數(shù)Δ趨于平穩(wěn),增加幅度不明顯。

        圖3 參數(shù)ΔI隨β的變化趨勢(shì)

        式(18)中的參數(shù)K(表2),隨著的增大而增大,整體呈增加趨勢(shì),增加幅度有所不同,如第2組數(shù)據(jù),參數(shù)K從0.018增大到0.225,增加幅度較大;第6組數(shù)據(jù)從2.081增大到2.098,增加幅度較小。當(dāng)>15時(shí),各組數(shù)據(jù)增加幅度均較小,參數(shù)K趨于平穩(wěn)。

        參數(shù)為模型中的(+s)(θ-θ),由表(2)可以看出,隨著的增大,各組入滲數(shù)據(jù)對(duì)應(yīng)的參數(shù)逐漸減小,同樣減小的程度有所不同,其中,第4組參數(shù)減小的程度最大,從147 0減小到34;其次是第2組數(shù)據(jù),從717減小到8.3;其他4組數(shù)據(jù)中,值也逐漸減小。當(dāng)>15時(shí),參數(shù)趨于平穩(wěn)。

        綜合來(lái)看,取不同值時(shí),不同入滲數(shù)據(jù)所對(duì)應(yīng)的相關(guān)系數(shù)均隨值的增大而增加,但當(dāng)>15時(shí),增加幅度明顯減小;各參數(shù)(K、A、Δ)隨的變化趨勢(shì)和變化幅度均有所不同,有的隨增大而增大,有的隨增大而減小,但當(dāng)>15時(shí),各參數(shù)均趨于穩(wěn)定,因此,建議?。?5。

        2.2.2 擬合效果

        整體來(lái)看,改進(jìn)模型較傳統(tǒng)模型擬合效果明顯提高,傳統(tǒng)模型6組相關(guān)系數(shù)平均值為0.979 0,改進(jìn)模型計(jì)算的相關(guān)系數(shù)增大到0.998 5。其中,第2、第4組數(shù)據(jù)效果更為顯著,如圖5。傳統(tǒng)模型和改進(jìn)模型對(duì)入滲后期的擬合效果理想,但入滲前期改進(jìn)模型的擬合效果要好于傳統(tǒng)模型,模擬的累計(jì)入滲時(shí)間與實(shí)測(cè)的累計(jì)入滲時(shí)間的相關(guān)系數(shù)分別從0.92和0.96提高到0.993 0和0.999 6;其次為第1組數(shù)據(jù),其相關(guān)系數(shù)從0.989提高到0.997;第3、第5、第6組數(shù)據(jù),改進(jìn)模型的相關(guān)系數(shù)較傳統(tǒng)模型也均有所提高。

        2.2.3 參數(shù)穩(wěn)定性分析

        改進(jìn)模型的參數(shù)(=15)比傳統(tǒng)模型的參數(shù)更為穩(wěn)定,并且更接近于實(shí)際測(cè)試的結(jié)果(表1、表2)。其中,改進(jìn)模型參數(shù)K的變化范圍(0.170~2.096)明顯小于傳統(tǒng)模型(0.001~2.068),對(duì)應(yīng)的變差系數(shù)分別為0.92和1.59;改進(jìn)模型參數(shù)的變化范圍(0.31~41.97)明顯小于傳統(tǒng)模型(3.4~15 705.3),對(duì)應(yīng)的變差系數(shù)分別為1.43和1.61;改進(jìn)模型參數(shù)Δ的變化范圍為26.42~119.08,變差系數(shù)為0.539,顯著小于參數(shù)K和的變差系數(shù),不同測(cè)試點(diǎn)的參數(shù)Δ較參數(shù)K和更為穩(wěn)定。

        3 結(jié)論

        1)基于濕潤(rùn)區(qū)分層的假設(shè),對(duì)傳統(tǒng)Green- Ampt入滲模型進(jìn)行了改進(jìn),得到了一個(gè)以分段函數(shù)表示的入滲模型。

        2)基于實(shí)測(cè)入滲資料,以累計(jì)入滲時(shí)間測(cè)試值與模擬值誤差平方和最小為目標(biāo),擬合確定了傳統(tǒng)Green-Ampt模型和改進(jìn)模型的參數(shù)值。結(jié)果表明,改進(jìn)模型顯著提高了擬合精度,其參數(shù)值較傳統(tǒng)模型更為合理和穩(wěn)定,穩(wěn)定入滲率更接近于測(cè)試值。

        3)值是改進(jìn)模型的一個(gè)重要參數(shù)值,分析表明,改進(jìn)模型的擬合精度(2)隨值的增大而增加,增加幅度快速減小,相應(yīng)參數(shù)隨值的增大,其變化幅度均有明顯減小的趨勢(shì),鑒于此,建議?。?5。

        [1] Green W H,Ampt G A.Studies on soil physics: Flow of air and water through soils[J].Journal of Agriculture Science,1911,4(1):1-24.

        [2] 雷志棟,楊詩(shī)秀,謝森傳. 土壤水動(dòng)力學(xué)[M]. 北京:清華大學(xué)出版社,1998.

        [3] 趙偉霞,張振華,蔡煥杰,等.恒定水頭井入滲Green-Ampt模型的改進(jìn)與驗(yàn)證[J].水利學(xué)報(bào),2010,41(4):464-470.

        [4] 郭向紅,孫西歡,馬娟娟,等.不同入滲水頭條件下的Green-Ampt模型[J].農(nóng)業(yè)工程學(xué)報(bào),2010,26(3):64-68.

        [5] 馬娟娟,孫西歡,郭向紅. 基于Green-Ampt模型的變水頭積水入滲模型建立及其參數(shù)求解[J]. 水利學(xué)報(bào),2010,41(01):61-67.

        [6] 劉姍姍,白美建,許迪,等. Green-Ampt模型參數(shù)簡(jiǎn)化及土壤物理參數(shù)的關(guān)系[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(1):106-110.

        [7] Coleman E A,Bodman G B.Moisture and energy conditions during downward entry of water into moist and layered soils[J].Soil Sci Soc Amer Proc,1945,9(C):3-11.

        [8] Ma Y,F(xiàn)eng S Y,Su D Y,et al.Modeling water infiltration in a large layered soil column with a modifies green-ampt model and HYDRUS-1D[J]. Computers and Electronics in Agriculture,2010,71(s1):S40-S47.

        [9] Bouwer H. Rapid field measurement of air-entry value and hydraulic conductivity of soil as significant parameters in flow system analysis [J].Water Resources Research, 1966,2(4):729-738.

        [10] Bouwer H. Infiltration of water into non-uniform soil[J].Journal of Irrigation and Drainage Engineering,1969,95 (IR4):451-462.

        [11] 毛麗麗,雷廷武,劉汗,等. 用水平土柱和修正的Green-Ampt模型確定土壤的入滲性能[J].農(nóng)業(yè)工程學(xué)報(bào),2009,25(11):35-40.

        [12] 王文焰,王全九,張建豐,等. 甘肅秦王川地區(qū)土壤水分運(yùn)動(dòng)參數(shù)及相關(guān)性[J].水土保持學(xué)報(bào),2002,16(3):110-113.

        [13] 彭振陽(yáng),黃介生,伍靖?jìng)?,? 基于分層假設(shè)的Green-Ampt模型改進(jìn)[J]. 水科學(xué)進(jìn)展,2012,23(1):59-66.

        [14] 周鑫洋,王仰仁,劉群昌. 四個(gè)入滲公式對(duì)農(nóng)田潮土的適應(yīng)性研究[J]. 灌溉排水學(xué)報(bào),2015,34(8):96-100.

        [15] 沈洪政,王仰仁,韓娜娜. 非線(xiàn)性入滲模型求參中參數(shù)初始值確定方法研究[J]. 中國(guó)農(nóng)村水利水電,2017 (10):183-187.

        責(zé)任編輯:宗淑萍

        Improvement of Green-Ampt model based on layered assumptions in wetting regions

        LI Yong-lin1,2, WANG Yang-ren1,2,Corresponding Author, SUN Xiao-ping3, ZHENG Zhi-wei1,2

        (1. College of Water Conservancy Engineering, Tianjin Agricultural University, Tianjin 300384, China; 2. University-Enterprise Collaborative Innovation Laboratory of Water-Saving Irrigation Technology and Equipment, Tianjin 300384, China; 3. Shanxi Institute of Water Resources and Hydropower Research, Taiyuan 030012, China)

        It is assumed that the infiltration process is divided into two stages of pre-infiltration and post-infiltration. The ratio of soil saturated zone to unsaturated zone is a constant 2in the pre-infiltration stage, and the thickness of unsaturated zone increases with infiltration. The thickness of unsaturated zone in the post-infiltration stage no longer changes with time, which is a fixed value,and only the thickness of saturated zone changes with time. Based on this, Green-Ampt model is improved. The traditional Green-Ampt model and improved model were used to fit six groups of infiltration test data. The correlation coefficient of traditional model fitting was 0.982(the average of six groups of data). Compared with the traditional model, the fitting accuracy of the improved model was obvious. The correlation coefficient was increased by 0.016, reaching 0.998 (corresponding to=15). The data of two groups were especially effective, which were increased by 0.07 and 0.04 respectively, both reaching over 0.99. The improved model has the same strong physical meaning as the traditional model. However, compared with the traditional model, the stable infiltration rate obtained by the improved model is closer to the test value and the improved model has wider applicability.

        Green-Ampt infiltration model; accumulated infiltration capacity; unsaturated zone; infiltration front; infiltration time

        1008-5394(2018)03-0078-06

        10.19640/j.cnki.jtau.2018.03.017

        S275

        A

        2018-05-01

        國(guó)家自然科學(xué)基金項(xiàng)目(51779174);天津市農(nóng)業(yè)科技成果轉(zhuǎn)化與推廣項(xiàng)目(201701150);天津市科技計(jì)劃項(xiàng)目(17YFZC SF00930);山西省水利科學(xué)技術(shù)研究與推廣項(xiàng)目(201618)

        李泳霖(1994-),男,碩士在讀,主要從事灌溉排水理論與新技術(shù)研究。E-mail:335985474@qq.com。

        王仰仁(1962-),男,教授,博士,主要從事灌溉排水技術(shù)研究。E-mail:wyrf@163.com。

        猜你喜歡
        非飽和濕潤(rùn)幅度
        The Desert Problem
        非飽和原狀黃土結(jié)構(gòu)強(qiáng)度的試驗(yàn)研究
        微波超寬帶高速數(shù)控幅度調(diào)節(jié)器研制
        海邊的沙漠
        非飽和多孔介質(zhì)應(yīng)力滲流耦合分析研究
        基于ANSYS的四連桿臂架系統(tǒng)全幅度應(yīng)力分析
        非飽和土基坑剛性擋墻抗傾覆設(shè)計(jì)與參數(shù)分析
        非飽和地基土蠕變特性試驗(yàn)研究
        2014年中期預(yù)增(降)幅度最大的50家上市公司
        對(duì)“進(jìn)步幅度分公式”的實(shí)踐思考
        亚洲无码视频一区:| 男女深夜视频网站入口| 国产人妖网站在线视频| 亚洲最大中文字幕在线| 国产农村熟妇videos| 国内露脸中年夫妇交换| 丝袜人妻无码中文字幕综合网| 国产丝袜美腿一区二区三区| 国产老熟妇精品观看| 成人黄色网址| 一级二级中文字幕在线视频| 99re6久精品国产首页| 久久午夜一区二区三区| 免费人成小说在线观看网站 | 亚洲国产精品区在线观看| 狠狠色噜噜狠狠狠777米奇小说| 中文字幕一区二区三区久久网站| 激情综合五月天开心久久| 青青草成人在线播放视频| 免费网站看v片在线18禁无码| 日本免费人成视频播放| 极品少妇被后入内射视| 亚洲成a人一区二区三区久久| 精品国产乱码久久久久久婷婷 | 国内精品久久人妻互换| 午夜精品免费视频一区二区三区| 国产二级一片内射视频播放| 亚洲欧洲无码一区二区三区| 99riav精品国产| 中文字幕亚洲精品专区| 狠狠色噜噜狠狠狠777米奇| 精品久久人人妻人人做精品| 中文国产成人精品久久一区| 日韩一区二区av极品| 久久不见久久见中文字幕免费| 亚洲国产成人无码影院| 日本一道本加勒比东京热| 国内精品久久久久久99| 国产精品半夜| 国产男女乱婬真视频免费| 蜜桃av人妻精品一区二区三区|