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

    灰壩工程領(lǐng)域巖體裂隙滲流浸潤(rùn)線及出漏點(diǎn)溶質(zhì)濃度數(shù)值模擬

    2017-06-19 15:08:10王婷婷徐文彬
    地下水 2017年3期
    關(guān)鍵詞:溶質(zhì)運(yùn)移滲流

    王婷婷,徐文彬

    (1.云南省電力投資有限公司,云南 昆明 650021;2.云南能投基礎(chǔ)設(shè)施投資開(kāi)發(fā)建設(shè)有限公司,云南 昆明 650021)

    ?

    灰壩工程領(lǐng)域巖體裂隙滲流浸潤(rùn)線及出漏點(diǎn)溶質(zhì)濃度數(shù)值模擬

    王婷婷1,徐文彬2

    (1.云南省電力投資有限公司,云南 昆明 650021;2.云南能投基礎(chǔ)設(shè)施投資開(kāi)發(fā)建設(shè)有限公司,云南 昆明 650021)

    灰壩工程中灰水滲漏易對(duì)地下水產(chǎn)生污染,為預(yù)防治理地下水污染,應(yīng)研究灰場(chǎng)污染物的運(yùn)移機(jī)理,在運(yùn)移工程中由于巖體裂隙的粗糙性和各向異性,滲流及溶質(zhì)運(yùn)移模型實(shí)際非常復(fù)雜。本文著眼于建立將粗糙裂隙巖體非連續(xù)滲流、溶質(zhì)運(yùn)移數(shù)學(xué)模型,應(yīng)用于陡立裂隙網(wǎng)絡(luò)巖體滲流浸潤(rùn)線及出漏點(diǎn)溶質(zhì)濃度數(shù)值模擬,為定量分析灰場(chǎng)對(duì)地下水的影響提供理論依據(jù)。

    裂隙巖體;滲流浸潤(rùn)線;出漏點(diǎn);數(shù)值模擬

    灰壩工程對(duì)環(huán)境的影響主要表現(xiàn)在灰水的滲漏上,分為對(duì)地下水和地表水的影響?;覉?chǎng)對(duì)地下水的影響主要表現(xiàn)在氟污染、微量元素滲出、PH值增大上。雖然土壤對(duì)氟離子有很強(qiáng)的吸附能力,但隨著氟離子的增加,將對(duì)地下水產(chǎn)生影響。為預(yù)防、治理地下水污染問(wèn)題,弄清灰場(chǎng)污染物運(yùn)移機(jī)理顯得十分重要。

    本文著眼于將裂隙巖體非連續(xù)滲流、溶質(zhì)運(yùn)移數(shù)學(xué)模型,應(yīng)用于陡立裂隙網(wǎng)絡(luò)巖體滲流浸潤(rùn)線及出漏點(diǎn)溶質(zhì)濃度數(shù)值模擬,為定量分析灰場(chǎng)對(duì)地下水的影響提供理論方法。

    1 裂隙巖體穩(wěn)定非連續(xù)滲流

    1.1 滲流數(shù)學(xué)模型

    這種模型忽略巖塊孔隙系統(tǒng)的滲透性,把巖體介質(zhì)系統(tǒng)看作單純的裂隙介質(zhì)系統(tǒng),每條裂隙的隙寬都是粗糙的、開(kāi)度變化的,滲流不滿足立方定律,具有6個(gè)特點(diǎn):

    (1)裂隙巖體中巖塊的滲透性一般很微弱,裂隙才是巖體中的水流通道,巖體裂隙網(wǎng)絡(luò)滲流具有定向性。

    (2)裂隙網(wǎng)絡(luò)滲流僅是一種網(wǎng)絡(luò)流。

    (3)裂隙網(wǎng)絡(luò)滲流是在巖體裂隙中非均勻的面狀流動(dòng),具有高度的非均勻性和各向異性。

    (4)裂隙網(wǎng)絡(luò)中的單一裂隙滲流不服從立方定律。

    (5)巖體的滲透性大小取決于巖體中結(jié)構(gòu)面的性質(zhì)及巖塊的巖性。

    (6)復(fù)雜裂隙系統(tǒng)中滲流,在裂隙交叉處,具有“偏流效應(yīng)”,即裂隙水流經(jīng)大小不等裂隙交叉處時(shí),水流偏向?qū)挻罅严兑粋?cè)流動(dòng)。

    假定某二維滲流區(qū)域有n個(gè)裂隙交叉點(diǎn)(節(jié)點(diǎn)),m個(gè)裂隙線單元,按流量平衡原理[1],節(jié)點(diǎn)i處的穩(wěn)定水流方程為

    (1)

    式中:[A]為裂隙網(wǎng)絡(luò)銜接矩陣:當(dāng)線單元j不銜接于節(jié)點(diǎn)i時(shí),aij=0;當(dāng)線單元j銜接于節(jié)點(diǎn)i方向時(shí),aij=1;當(dāng)線單元j銜接于節(jié)點(diǎn)i指向背離節(jié)點(diǎn)方向時(shí),aij=-1;C、b、l、H、Q分別為常數(shù)、線單元隙寬(裂隙粗糙)、線單元長(zhǎng)度、節(jié)點(diǎn)水頭、節(jié)點(diǎn)流量。

    1.2 溶質(zhì)運(yùn)移數(shù)學(xué)模型

    在二維裂隙網(wǎng)絡(luò)中,溶質(zhì)運(yùn)移過(guò)程主要受三個(gè)方面的控制:

    (1)裂隙網(wǎng)絡(luò)形態(tài),單裂隙粗糙程度;

    (2)水流狀態(tài);

    (3)溶質(zhì)本身的運(yùn)移擴(kuò)散、衰變和裂隙吸附作用等。

    在不考慮水流與巖塊骨架的耦合(裂隙壁水流通量為零),且溶質(zhì)巖壁吸附為線性平衡等溫吸附的前提下,溶質(zhì)運(yùn)移數(shù)值模型[2]為

    (2)

    其中:

    (3)

    式中:t為時(shí)間;e表示單元;i、j為單元節(jié)點(diǎn);l、b、u分別為單元長(zhǎng)、寬及水流流速;C為溶液中的溶質(zhì)濃度分布;qns是第二類(lèi)邊界上節(jié)點(diǎn)的外法線方向的溶質(zhì)通量;Nis是邊界單元上i節(jié)點(diǎn)在第二類(lèi)邊界上節(jié)點(diǎn)處的形函數(shù);αL為沿裂隙方向的彌散度;u水流流速;ξ是線單元局部橫坐標(biāo);kf為吸附分布系數(shù),實(shí)際上為各向異性;λ為衰變常數(shù)[3],等于ln(2/t1/2);t1/2為半衰期。

    2 陡立裂隙巖體浸潤(rùn)線及出漏點(diǎn)

    假定巖塊不透水,滲透水流只在巖體的裂隙網(wǎng)絡(luò)中流動(dòng),所以此時(shí)的滲流浸潤(rùn)線實(shí)際上是個(gè)裂隙中水位的連線。因而滲流浸潤(rùn)線的確定就轉(zhuǎn)化為各裂隙中水位的確定。

    同連續(xù)介質(zhì)滲流理論一樣,各裂隙中水位的確定同樣是一個(gè)非線性問(wèn)題,因而各裂隙中水位也是巖體裂隙網(wǎng)絡(luò)滲流的邊界。各裂隙中水位既要滿足其上的水頭等于其位置高度,又要滿足各裂隙的流量等于外界補(bǔ)給量。當(dāng)水流為穩(wěn)定流且浸潤(rùn)線上無(wú)補(bǔ)給排泄時(shí),浸潤(rùn)線上各裂隙流量等于零,即浸潤(rùn)線各裂隙兩端點(diǎn)水頭差等于零。這樣,各裂隙中浸潤(rùn)線邊界應(yīng)滿足以下方程

    (4)

    式中:Γ3為浸潤(rùn)線邊界,H和z分別表示水頭和位置高度。應(yīng)用非連續(xù)滲流數(shù)學(xué)模型確定浸潤(rùn)線上的裂隙的水位,可采用迭代法。具體實(shí)施步驟如下:

    (1)首先假定一個(gè)浸潤(rùn)線的位置,并假定此浸潤(rùn)線上各裂隙端點(diǎn)水頭等于其位置高度,此假定浸潤(rùn)線便成為第一類(lèi)邊界;

    (2)按巖體裂隙網(wǎng)絡(luò)滲流原理求解各節(jié)點(diǎn)的水頭和流量;

    (3)按式4中第二式調(diào)整浸潤(rùn)線位置,并再次假定調(diào)整后的浸潤(rùn)線上各裂隙端點(diǎn)水頭等于其位置高度;

    (4)重復(fù)2、3步,進(jìn)行迭代運(yùn)算,直到相鄰兩次迭代的浸潤(rùn)線位置相差在精度要求范圍內(nèi),此時(shí)的浸潤(rùn)線位置即為要確定的滲流浸潤(rùn)線位置。

    在滲流場(chǎng)確定的情況下,采用迭代法計(jì)算單元水流流速分布u;后將所有參數(shù)輸入到式2中,計(jì)算區(qū)域內(nèi)各時(shí)段節(jié)點(diǎn)溶質(zhì)濃度分布C及出漏電濃度值。

    3 實(shí)例分析

    某陡立巖體取50 cm×50 cm區(qū)域觀察,發(fā)現(xiàn)裂隙分布如圖1所示,有8條裂隙縱橫交錯(cuò),從下到上,從左到右,依次標(biāo)號(hào)為1~8。一次節(jié)點(diǎn)有32個(gè),線單元有40個(gè),若將一次線單元進(jìn)行二次剖分,將產(chǎn)生192個(gè)節(jié)點(diǎn),200個(gè)單元。裂隙隙寬服從均勻分布,每條裂隙的平均隙寬見(jiàn)表1所示。

    表1 裂隙基本參數(shù)

    模擬結(jié)果見(jiàn)圖1~圖3。

    圖1 裂隙巖體浸潤(rùn)線求解模型

    圖2 裂隙網(wǎng)絡(luò)浸潤(rùn)線及水頭分布(cm)

    從圖2中可以看出,裂隙巖體浸潤(rùn)線出現(xiàn)了階梯狀,這是由粗糙裂隙(線單元隙寬不同)下滲流的非均勻性造成的。需要說(shuō)明的是,裂隙巖體網(wǎng)絡(luò)滲流的浸潤(rùn)線表述的只是各縱向單裂隙中的水位,沒(méi)有實(shí)際的物理意義,水位之間的連線(浸潤(rùn)線)只是用于粗略估計(jì)線下巖體裂隙中水頭值。

    從圖3可以看出,污染物基本沿著裂隙進(jìn)行運(yùn)移,而且靠近浸潤(rùn)線的裂隙中,污染物運(yùn)移速度較快,這是因?yàn)榱严吨械乃髁魉佥^快。當(dāng)t=3 000 s后,整個(gè)裂隙網(wǎng)絡(luò)中污染物濃度分布基本保持穩(wěn)定,因?yàn)榇藭r(shí)由時(shí)間推移引起的濃度增加和因裂隙吸附、衰變作用引起的污染物濃度減少相平衡,進(jìn)入穩(wěn)定態(tài)。此時(shí),下游裂隙出流污染物通量為ΔF/C0=0.020 572 03 cm/s。為了弄清下游裂隙出流污染物通量隨時(shí)間的變化趨勢(shì),可以求出時(shí)間從t=0~3 000 s過(guò)程中,每間隔500 s時(shí)的出流通量值,然后繪出F/C0—T關(guān)系曲線,如圖4所示。

    圖3 裂隙巖體溶質(zhì)濃度分布C/C0

    圖4 下游裂隙出流污染物通量與時(shí)間關(guān)系

    從圖4中可以看出,隨著時(shí)間的推移,起始階段由于污染物尚未運(yùn)移到下游,因此,時(shí)間t<1 000 s時(shí)出流通量都為零,等到1 000 s

    特別注意的是,裂隙網(wǎng)絡(luò)滲流自由面出露點(diǎn)處的溶質(zhì)濃度是隨時(shí)間的變化量,需要采用迭代法求得:先假定滲流出露點(diǎn)處溶質(zhì)濃度和下游濃度一致為零,把它視為第一類(lèi)邊界條件上的節(jié)點(diǎn),參與整個(gè)區(qū)域的濃度計(jì)算;經(jīng)過(guò)一個(gè)Δt后,用與之最近節(jié)點(diǎn)上的濃度取代原來(lái)假定的濃度(零),再參與區(qū)域計(jì)算,如此反復(fù),經(jīng)過(guò)足夠長(zhǎng)的時(shí)間使之滿足迭代精度,即可求出整個(gè)區(qū)域的濃度分布和出露點(diǎn)處濃度。如圖3所示,t=3 000 s后,通過(guò)計(jì)算求得出露點(diǎn)處濃度為C/C0=0.005 447。

    4 結(jié)語(yǔ)

    (1)裂隙巖體滲流自由面與多孔連續(xù)介質(zhì)自由面不同,它僅僅表述的是裂隙中水位的連線,其滲流自由面出露點(diǎn)無(wú)論在何種流態(tài)下都只會(huì)出現(xiàn)在裂隙處。裂隙粗糙度對(duì)裂隙巖體網(wǎng)絡(luò)滲流自由面影響較大。

    (2)通過(guò)滲流出露點(diǎn)處濃度迭代法,可以推求裂隙巖體網(wǎng)絡(luò)滲流自由面下溶質(zhì)濃度分布和出露點(diǎn)處濃度值。

    (3)巖壁吸附能力(kf)隨時(shí)間推移會(huì)最終趨于飽和,吸附量和解析量持衡,換句話說(shuō),隨時(shí)間推移,巖體對(duì)水質(zhì)的凈化能力會(huì)慢慢消失。

    [1]柴軍瑞.大壩工程滲流力學(xué)[M].西藏:西藏人民出版社.2001.

    [2]徐文彬.巖體粗糙裂隙網(wǎng)絡(luò)非連續(xù)滲流及溶質(zhì)運(yùn)移研究[J].人民黃河.2014,93-95.

    [3]馮錫璋,等.放射性絕對(duì)測(cè)量和半衰期(Ⅱ)介紹一個(gè)非經(jīng)典的獲取衰變常數(shù)或半衰期的方法[J].高能物理與核物理.1982,6(1).

    The numerical simulation of the seepage line and the concentration of the solute in the seepage of rock mass in the area of the ash dam project

    WANG Ting-ting1,XU Wen-bin2

    (1.Yunnan Electric Power Investment Company Limited, Kunming, Yunnan, 650021 China;2.Take the Yunnan Infrastructure Investment Development Company Limited, Kunming, Yunnan, 650021 China)

    The ash water leakage is easy to pollute the groundwater. For the prevention and control of groundwater pollution, the mechanism of contaminant transport in ash field should be studied. Because of the roughness and anisotropy of rock fracture, the model of seepage and solute transport in nature is very complex in migration engineering. This paper focuses on applying mathematical model of discontinuous seepage flow in fractured rock mass and solute transport to numerical simulation of vertical fracture network of rock mass seepage lines as well as leakage concentration of solute. This can provide a theoretical basis for the quantitative analysis of the influence of ash field on groundwater.

    fractured rock mass;seepage line;leakage point;numerical simulation

    2016-12-28

    王婷婷(1982-),女,湖北宜昌人,工程師,主要從事水環(huán)境方面研究工作。

    徐文彬(1983-),男,湖北武漢人,高級(jí)工程師,主要從事滲流與環(huán)境方面工作。

    P641.135

    A

    1004-1184(2017)03-0030-03

    猜你喜歡
    溶質(zhì)運(yùn)移滲流
    有關(guān)溶質(zhì)質(zhì)量分?jǐn)?shù)的計(jì)算
    滴水成“冰”
    溶質(zhì)質(zhì)量分?jǐn)?shù)考點(diǎn)突破
    曲流河復(fù)合點(diǎn)壩砂體構(gòu)型表征及流體運(yùn)移機(jī)理
    東營(yíng)凹陷北帶中淺層油氣運(yùn)移通道組合類(lèi)型及成藏作用
    “溶質(zhì)的質(zhì)量分?jǐn)?shù)”計(jì)算歸類(lèi)解析
    開(kāi)采過(guò)程中上覆急傾斜巖層運(yùn)移規(guī)律模擬與研究
    川西坳陷孝泉-新場(chǎng)地區(qū)陸相天然氣地球化學(xué)及運(yùn)移特征
    簡(jiǎn)述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計(jì)算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    国产男靠女视频免费网站| 18禁裸乳无遮挡免费网站照片| 日韩av在线大香蕉| 国产单亲对白刺激| 亚洲国产欧美网| 精品国产超薄肉色丝袜足j| 在线视频色国产色| 俄罗斯特黄特色一大片| 国产精品美女特级片免费视频播放器 | 香蕉久久夜色| 国内少妇人妻偷人精品xxx网站 | 亚洲欧美日韩高清在线视频| 亚洲avbb在线观看| 亚洲无线观看免费| 国产乱人视频| 在线免费观看的www视频| 亚洲中文av在线| 99国产精品一区二区蜜桃av| 操出白浆在线播放| 久久久久久久久免费视频了| 精品福利观看| 美女免费视频网站| 九九热线精品视视频播放| 男女视频在线观看网站免费| av片东京热男人的天堂| 久久久精品大字幕| 国产精品99久久99久久久不卡| 国产三级中文精品| 久久久久久久久久黄片| 亚洲欧美日韩卡通动漫| 黄色片一级片一级黄色片| 国内久久婷婷六月综合欲色啪| 可以在线观看毛片的网站| 免费无遮挡裸体视频| 国产午夜福利久久久久久| 丁香欧美五月| 黄色视频,在线免费观看| 欧美乱妇无乱码| 久久精品国产99精品国产亚洲性色| 国产在线精品亚洲第一网站| 亚洲专区国产一区二区| 美女免费视频网站| 亚洲国产欧美人成| 免费观看的影片在线观看| 日本黄色片子视频| 精品久久久久久久人妻蜜臀av| 亚洲专区国产一区二区| 丁香欧美五月| 精品国产乱码久久久久久男人| 青草久久国产| 亚洲av电影不卡..在线观看| 国产伦一二天堂av在线观看| 国产亚洲av高清不卡| 2021天堂中文幕一二区在线观| 欧美黑人欧美精品刺激| 亚洲五月天丁香| e午夜精品久久久久久久| 亚洲人成网站高清观看| 久久香蕉精品热| 亚洲国产精品999在线| 哪里可以看免费的av片| 亚洲五月婷婷丁香| www.精华液| 亚洲国产欧美网| 欧美日韩黄片免| 国产又色又爽无遮挡免费看| 免费看美女性在线毛片视频| www.www免费av| 性欧美人与动物交配| 国产激情偷乱视频一区二区| av在线蜜桃| 久久精品aⅴ一区二区三区四区| 黄色视频,在线免费观看| 真人一进一出gif抽搐免费| 国产成人啪精品午夜网站| 在线观看舔阴道视频| 不卡一级毛片| 国产91精品成人一区二区三区| 在线免费观看的www视频| 69av精品久久久久久| 好男人在线观看高清免费视频| 99热这里只有精品一区 | 人人妻人人澡欧美一区二区| 在线永久观看黄色视频| 天天躁日日操中文字幕| 人妻久久中文字幕网| 久久精品国产综合久久久| 亚洲一区二区三区色噜噜| 久久亚洲精品不卡| 久久亚洲真实| 欧洲精品卡2卡3卡4卡5卡区| 一个人看视频在线观看www免费 | 他把我摸到了高潮在线观看| 欧美激情在线99| 97超级碰碰碰精品色视频在线观看| 中文字幕久久专区| 悠悠久久av| 亚洲av美国av| 亚洲成av人片免费观看| 亚洲精品乱码久久久v下载方式 | 亚洲成av人片在线播放无| 久久精品人妻少妇| a级毛片a级免费在线| 国产亚洲精品久久久com| 午夜成年电影在线免费观看| 一二三四在线观看免费中文在| 欧美在线黄色| 观看美女的网站| 欧美乱妇无乱码| 成年免费大片在线观看| 韩国av一区二区三区四区| 悠悠久久av| 在线免费观看的www视频| 午夜福利在线观看吧| 欧美另类亚洲清纯唯美| 一区二区三区国产精品乱码| 嫩草影视91久久| 国产精品一及| 日本一二三区视频观看| 久久久久久久久免费视频了| 不卡一级毛片| 美女黄网站色视频| 麻豆av在线久日| 蜜桃久久精品国产亚洲av| 在线a可以看的网站| 色综合亚洲欧美另类图片| 成人无遮挡网站| 一卡2卡三卡四卡精品乱码亚洲| 熟妇人妻久久中文字幕3abv| 精华霜和精华液先用哪个| 欧美性猛交╳xxx乱大交人| 国产麻豆成人av免费视频| 国产探花在线观看一区二区| 国产午夜福利久久久久久| 男女做爰动态图高潮gif福利片| 国产精品女同一区二区软件 | 蜜桃久久精品国产亚洲av| 国产精品国产高清国产av| 日韩人妻高清精品专区| 黄片小视频在线播放| 最近最新免费中文字幕在线| 老司机午夜福利在线观看视频| 国产免费男女视频| 午夜影院日韩av| 色综合站精品国产| 久久久久九九精品影院| 亚洲av成人精品一区久久| 亚洲成人免费电影在线观看| 国产私拍福利视频在线观看| 久久久色成人| 日韩大尺度精品在线看网址| ponron亚洲| 免费高清视频大片| aaaaa片日本免费| 亚洲 欧美一区二区三区| 九九在线视频观看精品| 中文字幕人妻丝袜一区二区| 亚洲人与动物交配视频| 99热精品在线国产| 999精品在线视频| 国产高清激情床上av| 亚洲熟妇熟女久久| 亚洲性夜色夜夜综合| 最近最新中文字幕大全免费视频| 麻豆成人午夜福利视频| av视频在线观看入口| 1024手机看黄色片| 亚洲第一欧美日韩一区二区三区| 亚洲人成伊人成综合网2020| 性欧美人与动物交配| 男女那种视频在线观看| 人人妻人人看人人澡| 亚洲欧美精品综合一区二区三区| 国产爱豆传媒在线观看| 男人舔奶头视频| www日本黄色视频网| 国产私拍福利视频在线观看| 欧美性猛交╳xxx乱大交人| 免费人成视频x8x8入口观看| 欧美成人一区二区免费高清观看 | 法律面前人人平等表现在哪些方面| 可以在线观看毛片的网站| 成人永久免费在线观看视频| 国产99白浆流出| 女人高潮潮喷娇喘18禁视频| 欧美午夜高清在线| 亚洲熟妇熟女久久| 69av精品久久久久久| 别揉我奶头~嗯~啊~动态视频| 亚洲精品456在线播放app | 久久午夜亚洲精品久久| 中文字幕高清在线视频| 国产av在哪里看| 亚洲av第一区精品v没综合| av片东京热男人的天堂| 麻豆成人av在线观看| 99热只有精品国产| 亚洲人成伊人成综合网2020| 亚洲午夜精品一区,二区,三区| 男女床上黄色一级片免费看| 国产麻豆成人av免费视频| 91九色精品人成在线观看| 久久人人精品亚洲av| 亚洲欧美一区二区三区黑人| 国产1区2区3区精品| 久久欧美精品欧美久久欧美| 亚洲黑人精品在线| 啦啦啦观看免费观看视频高清| 亚洲电影在线观看av| 亚洲精品国产精品久久久不卡| 黄色片一级片一级黄色片| 国产精品亚洲美女久久久| 午夜福利欧美成人| 欧美xxxx黑人xx丫x性爽| 一级毛片精品| 久久久久亚洲av毛片大全| 我要搜黄色片| 国产久久久一区二区三区| 久久精品91蜜桃| 九色成人免费人妻av| 国产成人一区二区三区免费视频网站| 老熟妇乱子伦视频在线观看| 久久精品影院6| 操出白浆在线播放| 夜夜看夜夜爽夜夜摸| 老汉色av国产亚洲站长工具| 免费在线观看成人毛片| 国产乱人伦免费视频| 色综合亚洲欧美另类图片| 亚洲欧美日韩高清专用| 色综合婷婷激情| 99久久久亚洲精品蜜臀av| 欧美一区二区国产精品久久精品| 国产单亲对白刺激| 午夜福利在线观看免费完整高清在 | 欧美成人一区二区免费高清观看 | 五月伊人婷婷丁香| 国产高清videossex| 亚洲专区国产一区二区| 免费看美女性在线毛片视频| 最新美女视频免费是黄的| 久久久久久人人人人人| 51午夜福利影视在线观看| 精品一区二区三区四区五区乱码| 免费av不卡在线播放| 国产精品一区二区三区四区久久| 久久精品国产99精品国产亚洲性色| 老熟妇乱子伦视频在线观看| 中文字幕人妻丝袜一区二区| x7x7x7水蜜桃| www日本在线高清视频| 色吧在线观看| 免费在线观看成人毛片| 国产精品亚洲一级av第二区| 日韩三级视频一区二区三区| 欧美中文日本在线观看视频| 日本成人三级电影网站| 中文字幕熟女人妻在线| 中出人妻视频一区二区| 精品久久久久久久人妻蜜臀av| 成人三级黄色视频| 青草久久国产| 波多野结衣巨乳人妻| 久久久成人免费电影| 操出白浆在线播放| 一个人看的www免费观看视频| 久久久久久久久久黄片| bbb黄色大片| 制服人妻中文乱码| 成人av一区二区三区在线看| 黄色视频,在线免费观看| 法律面前人人平等表现在哪些方面| 国产午夜福利久久久久久| 亚洲狠狠婷婷综合久久图片| 亚洲男人的天堂狠狠| 亚洲一区二区三区色噜噜| 久久精品国产99精品国产亚洲性色| 亚洲最大成人中文| 亚洲中文av在线| 久久午夜综合久久蜜桃| 精品国产亚洲在线| 露出奶头的视频| 国产成人福利小说| 久久久久国产一级毛片高清牌| 久久精品91蜜桃| 99国产精品一区二区蜜桃av| 99国产精品99久久久久| 国产欧美日韩精品亚洲av| 性欧美人与动物交配| 夜夜看夜夜爽夜夜摸| 黄色视频,在线免费观看| 欧美大码av| 色哟哟哟哟哟哟| 偷拍熟女少妇极品色| 高清在线国产一区| 国产三级在线视频| 中文字幕精品亚洲无线码一区| 日韩精品中文字幕看吧| 欧美一区二区国产精品久久精品| 欧美另类亚洲清纯唯美| 精华霜和精华液先用哪个| 男人的好看免费观看在线视频| 欧美日本视频| 18禁国产床啪视频网站| 美女 人体艺术 gogo| 国内精品一区二区在线观看| 日本 av在线| 亚洲成人久久爱视频| 色老头精品视频在线观看| 天堂动漫精品| 亚洲精品456在线播放app | 午夜福利高清视频| 观看美女的网站| www国产在线视频色| 亚洲成人免费电影在线观看| 午夜福利免费观看在线| 日本在线视频免费播放| 欧美日本亚洲视频在线播放| 国产淫片久久久久久久久 | 国产蜜桃级精品一区二区三区| 欧美黑人巨大hd| 国产1区2区3区精品| 国产精品香港三级国产av潘金莲| 一进一出抽搐动态| 免费无遮挡裸体视频| 久久精品综合一区二区三区| 国产高清videossex| 久久精品人妻少妇| 国模一区二区三区四区视频 | 国产主播在线观看一区二区| 成人性生交大片免费视频hd| 国产毛片a区久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 午夜精品在线福利| 曰老女人黄片| 国产一区二区三区在线臀色熟女| 中文在线观看免费www的网站| 人人妻,人人澡人人爽秒播| 夜夜爽天天搞| 91在线观看av| 亚洲精品国产精品久久久不卡| 99久久国产精品久久久| 熟女电影av网| а√天堂www在线а√下载| 亚洲av成人精品一区久久| 久久精品人妻少妇| bbb黄色大片| 日本成人三级电影网站| 丝袜人妻中文字幕| 亚洲欧美日韩东京热| 长腿黑丝高跟| 欧洲精品卡2卡3卡4卡5卡区| 丁香欧美五月| 欧美色欧美亚洲另类二区| 欧美日韩综合久久久久久 | 深夜精品福利| 精品99又大又爽又粗少妇毛片 | a级毛片a级免费在线| 久久伊人香网站| 日韩免费av在线播放| 中亚洲国语对白在线视频| 99久久99久久久精品蜜桃| 亚洲,欧美精品.| 可以在线观看的亚洲视频| 亚洲国产精品合色在线| 成在线人永久免费视频| 国产精品野战在线观看| 波多野结衣高清作品| 熟女少妇亚洲综合色aaa.| 国产午夜福利久久久久久| 日本a在线网址| 国产在线精品亚洲第一网站| 12—13女人毛片做爰片一| 色播亚洲综合网| www日本在线高清视频| 日日干狠狠操夜夜爽| 中文字幕高清在线视频| 精品乱码久久久久久99久播| 欧美另类亚洲清纯唯美| 18禁黄网站禁片免费观看直播| 少妇的丰满在线观看| 欧美+亚洲+日韩+国产| 亚洲国产色片| 老汉色∧v一级毛片| 后天国语完整版免费观看| 夜夜看夜夜爽夜夜摸| 亚洲精品久久国产高清桃花| 国产一级毛片七仙女欲春2| 国产精华一区二区三区| 国内精品美女久久久久久| 国产淫片久久久久久久久 | 一本一本综合久久| 亚洲熟女毛片儿| 黄频高清免费视频| 国产亚洲精品av在线| www.999成人在线观看| 国产精品久久电影中文字幕| 亚洲国产精品999在线| 午夜激情福利司机影院| 日韩免费av在线播放| 欧美乱妇无乱码| 黄色成人免费大全| 老鸭窝网址在线观看| 天堂av国产一区二区熟女人妻| 99在线人妻在线中文字幕| 国产人伦9x9x在线观看| 少妇丰满av| 美女 人体艺术 gogo| 悠悠久久av| 亚洲欧美日韩卡通动漫| 黑人巨大精品欧美一区二区mp4| 18禁美女被吸乳视频| 可以在线观看毛片的网站| 88av欧美| 中文在线观看免费www的网站| 美女黄网站色视频| 午夜成年电影在线免费观看| 国产精品综合久久久久久久免费| 这个男人来自地球电影免费观看| 黄色成人免费大全| 国产精品久久视频播放| 国产亚洲精品久久久久久毛片| 99国产精品一区二区蜜桃av| 国产毛片a区久久久久| 九九热线精品视视频播放| 99热精品在线国产| 国产精品一区二区三区四区久久| 中出人妻视频一区二区| 亚洲欧美日韩无卡精品| 久久久久性生活片| 97碰自拍视频| 黄色片一级片一级黄色片| 一区二区三区高清视频在线| 国产亚洲精品久久久com| 国产精品亚洲一级av第二区| 欧美在线一区亚洲| 国产精品亚洲美女久久久| 久久久精品欧美日韩精品| 中文字幕高清在线视频| 亚洲av美国av| 国产成人啪精品午夜网站| 国内少妇人妻偷人精品xxx网站 | 久久欧美精品欧美久久欧美| 夜夜爽天天搞| 一级毛片高清免费大全| 欧美一区二区国产精品久久精品| 一级毛片高清免费大全| 在线观看美女被高潮喷水网站 | 美女大奶头视频| 婷婷丁香在线五月| 日韩 欧美 亚洲 中文字幕| 美女高潮的动态| 在线观看66精品国产| 精品久久久久久久毛片微露脸| 日日夜夜操网爽| 亚洲av美国av| 夜夜躁狠狠躁天天躁| 18禁国产床啪视频网站| 国产av一区在线观看免费| 国产精品久久久久久亚洲av鲁大| 久久久国产成人免费| 黄色 视频免费看| 国产亚洲精品久久久com| 一级毛片精品| 亚洲av成人av| 中文在线观看免费www的网站| 九色成人免费人妻av| 精品一区二区三区av网在线观看| av欧美777| 又爽又黄无遮挡网站| 男人舔奶头视频| 欧美色欧美亚洲另类二区| 伊人久久大香线蕉亚洲五| 在线观看日韩欧美| 精品久久久久久,| 国产97色在线日韩免费| 少妇人妻一区二区三区视频| 99热这里只有精品一区 | 国产成人影院久久av| 一个人看的www免费观看视频| 一个人观看的视频www高清免费观看 | 成熟少妇高潮喷水视频| 男人舔女人下体高潮全视频| 亚洲国产精品久久男人天堂| 在线永久观看黄色视频| 国产av麻豆久久久久久久| 欧美一区二区国产精品久久精品| 9191精品国产免费久久| 亚洲精品在线观看二区| 成人av在线播放网站| 成人av一区二区三区在线看| 一级作爱视频免费观看| 午夜福利高清视频| 精品一区二区三区视频在线 | 亚洲国产精品999在线| 又爽又黄无遮挡网站| 在线视频色国产色| 成年女人毛片免费观看观看9| 国产精品久久久av美女十八| 99国产精品99久久久久| 舔av片在线| 国产精品99久久久久久久久| 亚洲av第一区精品v没综合| 国产精品一区二区三区四区免费观看 | 91麻豆av在线| 一区福利在线观看| 亚洲第一电影网av| 久久中文看片网| 中文在线观看免费www的网站| 久久精品夜夜夜夜夜久久蜜豆| 在线播放国产精品三级| 88av欧美| 亚洲av中文字字幕乱码综合| 一个人免费在线观看电影 | 亚洲国产精品成人综合色| 白带黄色成豆腐渣| 俺也久久电影网| 亚洲熟女毛片儿| 国内毛片毛片毛片毛片毛片| 久久久久久久久中文| 国产一区二区三区在线臀色熟女| 搡老熟女国产l中国老女人| 成人特级av手机在线观看| 欧美绝顶高潮抽搐喷水| 国产成人精品久久二区二区91| 亚洲欧美精品综合一区二区三区| 俺也久久电影网| 国产极品精品免费视频能看的| 欧美成狂野欧美在线观看| 国产成+人综合+亚洲专区| 天天添夜夜摸| 日本成人三级电影网站| 久久久国产成人免费| 久久精品人妻少妇| 最新美女视频免费是黄的| or卡值多少钱| 亚洲人成网站在线播放欧美日韩| 99热6这里只有精品| 国产免费av片在线观看野外av| 国产激情偷乱视频一区二区| 久久中文字幕人妻熟女| 国产免费男女视频| 免费观看的影片在线观看| 人妻夜夜爽99麻豆av| 99热6这里只有精品| 此物有八面人人有两片| 色综合欧美亚洲国产小说| 在线永久观看黄色视频| 中文字幕久久专区| 极品教师在线免费播放| 精品不卡国产一区二区三区| 国产欧美日韩一区二区精品| 免费看a级黄色片| 91麻豆av在线| 热99在线观看视频| 淫妇啪啪啪对白视频| 天堂√8在线中文| 国产伦精品一区二区三区四那| 欧美性猛交黑人性爽| 在线国产一区二区在线| 精品国产乱码久久久久久男人| 成年版毛片免费区| 给我免费播放毛片高清在线观看| 在线观看免费视频日本深夜| 非洲黑人性xxxx精品又粗又长| 搡老岳熟女国产| 精品久久久久久久久久久久久| 美女扒开内裤让男人捅视频| 岛国视频午夜一区免费看| 精品久久久久久久人妻蜜臀av| 亚洲九九香蕉| 精品一区二区三区视频在线观看免费| 午夜精品久久久久久毛片777| 国产精品影院久久| 国产精品久久电影中文字幕| 日韩三级视频一区二区三区| 在线观看舔阴道视频| 久久亚洲真实| 一区二区三区高清视频在线| 午夜激情欧美在线| 五月伊人婷婷丁香| 精品99又大又爽又粗少妇毛片 | 国产精品av视频在线免费观看| 两个人看的免费小视频| 国产精品九九99| 国产熟女xx| 午夜福利在线观看吧| 国产高潮美女av| 亚洲欧美激情综合另类| 国产精品免费一区二区三区在线| 看片在线看免费视频| 欧美成人一区二区免费高清观看 | 午夜免费观看网址| 小蜜桃在线观看免费完整版高清| 男女视频在线观看网站免费| 亚洲真实伦在线观看| 一本精品99久久精品77| 久久久精品欧美日韩精品| 精品一区二区三区视频在线观看免费| 99久国产av精品| 国产一区二区激情短视频| 欧美日韩综合久久久久久 | 亚洲熟女毛片儿| www国产在线视频色| 久久香蕉国产精品| 久久国产精品影院| 老司机午夜福利在线观看视频| 国产黄a三级三级三级人| av福利片在线观看| 99久久精品国产亚洲精品| 日本精品一区二区三区蜜桃| 成人国产综合亚洲| 午夜激情福利司机影院| 日本精品一区二区三区蜜桃| 亚洲午夜精品一区,二区,三区|