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

    丘陵山區(qū)地下水流動特征及核素遷移數(shù)值模擬

    2019-05-31 01:54:38朱君陳超李婷張艾明
    山西大學學報(自然科學版) 2019年2期
    關鍵詞:沖溝丘陵第四系

    朱君,陳超,李婷,張艾明

    (中國輻射防護研究院 核環(huán)境模擬與評價技術重點實驗室,山西 太原 030006)

    0 引言

    據(jù)全球核電機組的數(shù)量和分布統(tǒng)計,一半都建設在內(nèi)陸,美國、俄羅斯、加拿大等國家則更多。目前,我國的核電機組主要分布在東南沿海地區(qū),未來發(fā)展到一定規(guī)模,沿海沒有足夠的地方,內(nèi)陸核電的發(fā)展難以避免。2008年,湖南桃花江、湖北咸寧和江西彭澤三大內(nèi)陸核電項目正式獲得批準并開展前期準備工作,但是三年后的日本福島核事故,出于安全因素考慮內(nèi)陸核電項目全部進入擱淺期。

    根據(jù)取排水特點,內(nèi)陸核電廠址大部分將處于大江、河流等丘陵山區(qū)。氚(H-3)是核反應堆釋放的特征污染物,1994-2005年,大亞灣核電站共排放液態(tài)氚5.6×105GBq。作為一種弱β放射性核素,不會對人體造成外照射危害,但是由于半衰期長,直接或者間接進入人體后,易與DNA等有機質(zhì)結(jié)合,對組織和器官造成內(nèi)照射,誘發(fā)突變。更重要的是H-3不受地下水環(huán)境巖土介質(zhì)的吸附和阻滯作用,進入含水層后其遷移速度與地下水流速基本上一致,危害極大。因此,合理評估H-3在丘陵山區(qū)地下水中的遷移、轉(zhuǎn)化行為,對內(nèi)陸核電的選址、設計、建造以及事故應急措施的制定是十分迫切和必要的。

    丘陵山區(qū)在沖溝、河谷階地、山梁等不同的地貌單元組合下,地下水系統(tǒng)的補給與徑流變得十分復雜,并且第四系松散層受侵蝕作用,僅堆積在沖溝底部和河谷,受地形阻隔,地層的不連續(xù)決定了沖溝間第四系松散孔隙地下水的相對獨立性,如何描述第四系松散孔隙地下水,在數(shù)值模擬中合理的建立丘陵山區(qū)的水文地質(zhì)概念模型是難點。目前相關模型通常分為兩類,第一類是將第四系松散層與基巖的交界處作為邊界,將沖溝或者河谷單獨分割出來[1-2],該方法局限于第四系松散孔隙單一含水層且范圍受限,對研究松散孔隙地下水與碎屑巖裂隙地下水之間的轉(zhuǎn)化,以及多地層三維數(shù)值模擬顯然不適用。如,吳樂[3]采用GMS數(shù)值模擬軟件,建立了北京西山地區(qū)地下水流三維模型,分析南水北調(diào)條件下,不同開采方案的區(qū)域含水層系統(tǒng)響應特征,雖然研究是針對西山地區(qū),但模型以山前平原分界線為界,實質(zhì)上只考慮了山前沖洪積扇地區(qū),而將山區(qū)地下水作為側(cè)向徑流補給賦予模型;龔繼文[4]基于GMS平臺以重慶某山區(qū)為例,構(gòu)建三維水文地質(zhì)模型,得到山區(qū)的地下水流場空間分布,但是模擬結(jié)果只是砂巖裂隙含水層、碳酸鹽巖溶含水層的地下水流場,沒有耦合分布在沖溝與河漫灘的第四系松散孔隙含水層。另一類是將松散孔隙含水層與碎屑巖裂隙含水層整體考慮,得到統(tǒng)一的地下水流場[5-6],簡單概化可能造成地下水水位線在沖溝上“懸空”,或者是增加沖溝第四系松散孔隙地下水的水力梯度,造成污染物遷移過快,無法真實地反映丘陵山區(qū)地下水流動特征下的污染物遷移規(guī)律。

    地下水數(shù)值模擬已成為解決地下水資源開發(fā)、定量評價和溶質(zhì)遷移等各領域問題的重要手段,從地形地貌特征劃分,國內(nèi)研究主要集中在平原區(qū)和盆地區(qū)[7-17]。例如,胡立堂[18]應用TOUGH2軟件建立了北京市平原區(qū)飽和~非飽和地下水三維流模型,研究地下水供水方案的合理性,其結(jié)果為南水北調(diào)入京后北京市地下水補給和壓采方案提供了有效的技術支持。董少剛[19]分析了太原盆地地下水三維流動特征,建立松散沉積層的地質(zhì)結(jié)構(gòu)模型和地下水的三維流動數(shù)學模型,利用數(shù)值模擬方法再現(xiàn)了太原盆地地下水流動的時空變化規(guī)律。但是真正能夠刻畫丘陵山區(qū)地下水系統(tǒng)補給、徑流和排泄特點的模型,鮮有報道。

    國外已經(jīng)將丘陵山區(qū)的地下水數(shù)值模型成功用于評估乏燃料后處理廠、核廢物處置場對地下水環(huán)境的影響。EUREX后處理廠始建于20世紀60年代,位于意大利北部山前沖積傾斜平原,2006年監(jiān)測發(fā)現(xiàn)地下水中存在放射性核素90Sr, Sanctis等[20]應用FEFLOW建立了三維地下水核素遷移模型,通過監(jiān)測數(shù)據(jù)對模型進行了驗證和率定,成功的反映了90Sr在地下水系統(tǒng)中的遷移、轉(zhuǎn)化規(guī)律。Gmünder等[21]應用FEFLOW建立了瑞士北部山區(qū)深地質(zhì)高放廢物處置場的區(qū)域三維水文地質(zhì)模型,評估了地下水滲流場對處置場的安全影響;并詳細論述了如何在模型中精細概化地層剝蝕出露、斷層等復雜地質(zhì)構(gòu)造[22],以及對地下水滲流場的影響。Luo等[23]采用近、遠場三維嵌套模型的方法,同時建立了瑞士北部山區(qū)6個中、低放廢物處置場的三維地下水數(shù)值模型,刻畫了地下水系統(tǒng)的補給、徑流和排泄特征,為選址及安全評估提供科學依據(jù);并實現(xiàn)了在模型中概化構(gòu)造及斷層[24-26],闡述優(yōu)先通道對核素遷移、擴散的影響。

    內(nèi)陸核電的地形地貌特征有別于平原區(qū)和盆地區(qū),如何準確的描述丘陵山區(qū)地下水流動的復雜特性,是現(xiàn)階段亟須解決的問題。本文將選擇山西汾河河谷的典型丘陵山區(qū)作為研究對象,擬從三維地形地質(zhì)模型的構(gòu)建著手,應用GOCAD地質(zhì)建模軟件,刻畫地層的分布、剝蝕以及傾向特點;在結(jié)合地下水數(shù)值模擬軟件FEFLOW,構(gòu)建三維水文地質(zhì)模型,精細概化地下水系統(tǒng)的補給、徑流和排泄特征,以期解決丘陵山區(qū)復雜的地下水流動問題。

    1 建模方法

    1.1 地質(zhì)建模

    選擇具有代表性的山西寧武縣西馬坊鄉(xiāng)~新堡鄉(xiāng)的汾河河谷作為研究對象,位于寧靜向斜西翼,傾角5~78°,下伏基巖地層剝蝕出露,呈狹長的條狀分布,地質(zhì)構(gòu)造復雜,屬于下靜游巖溶泉域??傮w地勢西北高、東南低,最高海拔1 681 m,位于西部山梁;最低海拔1 320 m,位于汾河階地,由黃土峁梁、沖溝、河谷階地組成,屬于典型的丘陵山區(qū)地貌。

    地層由新到老依次分為第四系、侏羅系、三疊系。第四系主要分布于山間沖溝和汾河河谷,上部為灰黃色亞砂土、細砂土,下部為礫石層。侏羅系由中統(tǒng)云崗組(J2y)紫紅色砂質(zhì)泥巖、灰綠色含礫長石砂巖、黃綠色粉砂巖;大同組(J2d)不等厚互層的灰白色石英砂巖、灰色砂質(zhì)泥巖組成。三疊系由中統(tǒng)銅川組(T2t)灰紅色長石砂巖、黑色頁巖和二馬營組(T2er)肉紅色長石砂巖、紫色砂質(zhì)巖組成。

    為了能夠刻畫丘陵山區(qū)地下水系統(tǒng)的補給、徑流和排泄條件。必須在模型中盡可能真實的反映研究區(qū)的地形地質(zhì)條件,因此需要建立三維數(shù)學模型,并在模型中刻畫第四系松散層的分布特點,以及侏羅系云崗組和大同組、三疊系銅川組和二馬營組地層的出露、剝蝕以及走向、傾向特點。

    應用GOCAD軟件建立三維地形地質(zhì)模型,地形數(shù)據(jù)從1:1 000圖上矢量化獲取,模型自上而下分為第四系松散層、侏羅系云崗組和大同組、三疊系銅川組和二馬營組地層。采用三角網(wǎng)格剖分,除第四系松散層,其余地層均剖分為2層,因此空間上共剖分為9層,在沖溝和汾河河谷區(qū)加密,共產(chǎn)生節(jié)點912 780個,單元格1 634 139個。

    1.2 水文地質(zhì)建模

    研究區(qū)地下水類型主要分為第四系松散孔隙地下水和碎屑巖裂隙地下水。松散孔隙含水層主要分布在較大的沖溝和汾河河谷,由亞砂土、細砂土、沙礫石層組成,水位埋深淺,接受大氣降水、地表水的入滲補給和基巖風化裂隙水的側(cè)向補給,地層的不連續(xù)決定了沖溝間第四系松散孔隙地下水的相對獨立性,不可能穿越基巖山體發(fā)生水力聯(lián)系,所有沖溝間的第四系松散孔隙地下水沿溝谷向下游排泄,最后匯入汾河河谷產(chǎn)生聯(lián)系。碎屑巖裂隙含水層與上覆第四系含水層有較好的水力聯(lián)系,在出露區(qū)接受大氣降水補給,地下水順巖層傾斜由向斜兩翼向軸部運動,具承壓性。

    將GOCAD建立的地形地質(zhì)模型導入地下水數(shù)值模擬軟件FEFLOW,同時模擬區(qū)地下水主要賦存于松散孔隙含水層和碎屑巖類裂隙含水層兩類。自上而下劃分第四系松散孔隙含水層、侏羅系云崗組砂巖裂隙含水層、侏羅系大同組砂巖裂隙含水層、三疊系銅川組砂巖裂隙含水層和二馬營組砂巖裂隙含水層。

    模型的西北部(AB段)以三疊系二馬營組出露處為界,處理為流量邊界;東南部(CD段)以汾河為界,處理為定水頭邊界;東北部(BC段)以地表分水嶺為界,西南部(AD段)以沖溝為界,均處理為通量邊界。整個模擬區(qū)東西長約9.0 km,南北寬約11.0 km,面積約82.0 km2。各含水層系統(tǒng)在地層出露處接受大氣降水補給。另外,地下水系統(tǒng)符合質(zhì)量守恒定律、能量守恒定律和達西定律;考慮污染物運移的特點,地下水運動概化成空間三維穩(wěn)定流;參數(shù)隨空間變化,體現(xiàn)了系統(tǒng)的非均質(zhì)性。因此,將整個模型概化為非均質(zhì)各向異性、空間三維結(jié)構(gòu)、穩(wěn)定地下水流系統(tǒng)(見圖2)。

    水文地質(zhì)參數(shù)主要包括含水層x、y、z方向的主滲透系數(shù)Kxx、Kyy、Kzz和大氣降水補給入滲系數(shù)。滲透系數(shù)Kxx、Kyy取值主要參考煤礦水文地質(zhì)勘查成果中的注水、壓水和抽水實驗獲取的實測值,并在模型識別過程中進行修正,垂直方向的滲透系數(shù)Kzz取水平方向滲透系數(shù)的1/10。該地區(qū)多年平均降雨量為468.1 mm,降雨入滲系數(shù)α參考經(jīng)驗值,松散巖類取0.06,碎屑巖類取0.24(見表1)。

    表1 含水層概化及參數(shù)取值

    1.3 核素遷移溶質(zhì)建模

    從污染物遷移途徑考慮,第四系松散孔隙含水層作為首要收納水體,最先受到污染,且與地表水聯(lián)系最為緊密,研究核素在丘陵山區(qū)的遷移、轉(zhuǎn)化規(guī)律,必定不可忽略第四系松散孔隙含水層。如何精細的描述丘陵山區(qū)復雜地貌條件的地下水流動特征,是研究核素遷移運動規(guī)律的關鍵問題。

    對三維地質(zhì)、水文地質(zhì)模型進行識別和驗證后,加入對流~彌散模塊,建立核素遷移溶質(zhì)模型。H-3作為核反應堆釋放的特征污染物,重要的性質(zhì)是不受地下水巖土介質(zhì)的吸附和阻滯作用,遷移速度與地下水流速基本上一致。因此,在溶質(zhì)遷移模型中不考慮吸附、化學反應,仍以對流與彌散為主,加入衰變。

    (1)

    式中:t為核素運移時間,d;C為溶解相濃度,Bq/cm3;λ為衰變常量,d-1;Dxx、Dyy、Dzz為x、y、z主方向的彌散系數(shù),cm2/d;qxx、qyy、qyy為x、y、z主方向的流速。

    另外,彌散度是影響H-3遷移的一個關鍵參數(shù),參考文獻和經(jīng)驗值第四系松散孔隙含水層縱向彌散度αL取5.0 m,碎屑巖類裂隙含水層縱向彌散度αL取0.3 m,橫向彌散度αT取縱向彌散度的1/10。假設H-3每天泄漏1 000 Bq進入地下水,按連續(xù)泄漏考慮,計算30 a和60 a H-3污染暈隨時間和空間的變化情況。

    2 結(jié)果和討論

    2.1 地下水水流數(shù)值模擬結(jié)果

    將邊界條件、水文地質(zhì)參數(shù)等代入模型,并且以穩(wěn)定流運行得到模型的地下水水位等值線,見圖2。由計算水位線分布可知,第四系松散孔隙地下水在接受兩側(cè)基巖裂隙水的補給后,向沖溝匯集,最終由西向東排泄于汾河。碎屑巖裂隙地下水受寧靜向斜地質(zhì)構(gòu)造的影響,地下水順巖層傾斜由向斜兩翼向軸部運動,總體流向由西北向東南徑流,具有承壓性。模擬結(jié)果與水文地質(zhì)調(diào)查所描述的地下水流動特征基本一致。

    Fig.2 Groundwater flow of model圖2 模型地下水流場圖

    模型的識別與驗證是模擬過程中重要的一步工作,通過擬合實測地下水水位,識別水文地質(zhì)參數(shù)、邊界條件,使建立的模型更加符合研究區(qū)的水文地質(zhì)條件,以便更精確的定量模擬區(qū)地下水的補、徑、排。

    調(diào)查研究區(qū)內(nèi)的18個村莊的飲用水水井資料,包括地下水水位標高和含水層類型。一方面用于分析研究區(qū)的地下水流動特點;另一方面可以作為水位實測資料,驗證模型的合理性。共有飲用水水井26個,23個取自第四系松散孔隙含水層,分布在沖溝和汾河河谷;其余3個取自碎屑巖裂隙含水層。同時,收集研究區(qū)煤礦的4個水文地質(zhì)鉆孔資料(見表2)。

    表2 水位實測數(shù)據(jù)

    由擬合結(jié)果可知,除了16#和21#點位外,地下水實測水位與預測水位的殘差平方和SSE為26.02,見圖3。所建立的水文地質(zhì)模型、邊界條件以及水文地質(zhì)參數(shù)的選取能夠反映丘陵山區(qū)復雜地形地貌條件下的地下水流運動特征,可以利用模型進一步預測核素遷移、轉(zhuǎn)化規(guī)律。

    Fig.3 Fitting of groundwater level圖3 地下水水位擬合

    2.2 核素遷移數(shù)值模擬結(jié)果

    計算30年后H-3的濃度暈中心值為23.0 Bq/L,0.1 Bq/L濃度暈沿溝谷向下游遷移了約700 m,第四系松散孔隙地下水沿程水頭變化范圍1 443 m~1 431 m,水力梯度約0.017,流速0.057 m/d(見圖4)。

    60年后H-3的濃度暈中心值為25.0 Bq/L,基本上處于穩(wěn)定,0.1 Bq/L濃度暈沿溝谷向下游遷移了約1 650 m,沿程水頭變化范圍1 443 m~1 407 m,水力梯度約0.022,流速0.073 m/d(見圖5)。

    污染物遷移速度取決于地下水的水力梯度,而丘陵山區(qū)沖溝與河谷中的第四系松散孔隙地下水水力梯度并不偏大。針對丘陵山區(qū)復雜的地下水流動特征,也有將松散孔隙含水層與碎屑巖裂隙含水層統(tǒng)一考慮,得到統(tǒng)一的地下水流場[5-6],這種簡單的概化造成結(jié)果嚴重的失真,一方面可能造成地下水水位線在沖溝上“懸空”;另一方面增加了沖溝第四系松散孔隙地下水的水力梯度,造成污染物遷移過快,參考上述文獻,地下水流速會到達6.2 m/d,無法真實的反映丘陵山區(qū)地下水流動特征下的污染物遷移規(guī)律。

    Fig.4 Migration distribution of H-3 after 30 years圖4 30年H-3遷移分布圖

    Fig.5 Migration distribution of H-3 after 60 years圖5 60年H-3遷移分布圖

    3 結(jié)論

    本文選擇山西汾河河谷典型丘陵山區(qū)作為研究對象,首先應用GOCAD軟件建立三維地形地質(zhì)模型,刻畫了第四系松散層、侏羅系云崗組和大同組、三疊系銅川組和二馬營組地層的分布、剝蝕以及傾向特點;然后結(jié)合地下水數(shù)值模擬軟件FEFLOW,構(gòu)建三維水文地質(zhì)模型。計算結(jié)果可以精細的概化丘陵山區(qū)第四系松散孔隙水、碎屑巖裂隙水的補給、徑流和排泄特征,第四系松散孔隙地下水在接受兩側(cè)基巖裂隙水的補給后,向沖溝匯集,最終由西向東排泄于汾河;碎屑巖裂隙地下水受寧靜向斜地質(zhì)構(gòu)造的影響,地下水順巖層傾斜由向斜兩翼向軸部運動,總體流向由西北向東南徑流。

    以不受巖土介質(zhì)吸附和阻滯的H-3作為內(nèi)陸核電的特征污染物,按每天連續(xù)泄漏1 000 Bq考慮,在第四系松散孔隙地下水中,60年后H-3的0.1 Bq/L濃度線沿溝谷向下游遷移了約1 650 m,遷移速度取決于地下水的水力梯度,并沒有想象中的偏大,而將松散孔隙含水層與碎屑巖裂隙含水層地下水流場統(tǒng)一考慮,人為增加了污染物遷移速度,造成結(jié)果失真。

    猜你喜歡
    沖溝丘陵第四系
    第四系膠結(jié)砂層水文地質(zhì)鉆進技術研究初探
    加快推進丘陵山區(qū)農(nóng)業(yè)生產(chǎn)托管服務
    風機基礎設計在沖溝發(fā)育地區(qū)的應用
    胡丘陵《栽樹》
    散文詩(2017年22期)2017-06-09 07:55:21
    關于推進丘陵山區(qū)農(nóng)業(yè)機械化發(fā)展的思考
    河北平原區(qū)第四系(200m以淺)巖土體綜合熱導率特征綜述
    關于丘陵山區(qū)農(nóng)機化發(fā)展的幾點思考
    北疆第四系覆蓋層上明渠結(jié)構(gòu)優(yōu)化研究
    沖溝不同部位土壤機械組成及抗沖性差異①
    土壤(2016年6期)2016-02-08 07:53:18
    基于高分辨率影像的沖溝發(fā)育及其時空變化研究
    日韩免费高清中文字幕av| 最新中文字幕久久久久| 少妇人妻久久综合中文| 亚洲高清免费不卡视频| 特大巨黑吊av在线直播| 国产又色又爽无遮挡免| 91久久精品电影网| 777米奇影视久久| 一二三四中文在线观看免费高清| 亚洲精品456在线播放app| 日本一本二区三区精品| 国产伦理片在线播放av一区| 国产高潮美女av| 国产高潮美女av| 国产精品人妻久久久影院| 丰满少妇做爰视频| videos熟女内射| 国产黄片美女视频| 国产亚洲5aaaaa淫片| 久久热精品热| 免费看日本二区| 在线免费十八禁| 男人爽女人下面视频在线观看| av女优亚洲男人天堂| 麻豆精品久久久久久蜜桃| 久久久久久久国产电影| 亚洲av成人精品一区久久| 另类亚洲欧美激情| 一个人看视频在线观看www免费| 国产熟女欧美一区二区| 中文字幕av成人在线电影| 少妇被粗大猛烈的视频| 日韩三级伦理在线观看| 午夜亚洲福利在线播放| 国产精品久久久久久精品古装| 国产精品不卡视频一区二区| 国产 一区精品| 国产综合精华液| 青春草国产在线视频| 赤兔流量卡办理| 80岁老熟妇乱子伦牲交| 男人舔奶头视频| 国产一区二区亚洲精品在线观看| 美女高潮的动态| 国产精品不卡视频一区二区| 日韩中字成人| 插阴视频在线观看视频| 青春草视频在线免费观看| 夜夜看夜夜爽夜夜摸| 亚洲av.av天堂| 男男h啪啪无遮挡| 熟妇人妻不卡中文字幕| av福利片在线观看| 国产亚洲91精品色在线| 免费黄色在线免费观看| 特级一级黄色大片| 一区二区av电影网| 国产精品嫩草影院av在线观看| 午夜精品一区二区三区免费看| 国产精品一区二区性色av| 久久久国产一区二区| 中文乱码字字幕精品一区二区三区| 国产成人a区在线观看| 寂寞人妻少妇视频99o| 日韩大片免费观看网站| 2021少妇久久久久久久久久久| av在线老鸭窝| 中文字幕av成人在线电影| 男女国产视频网站| 蜜桃亚洲精品一区二区三区| 亚洲国产日韩一区二区| 色吧在线观看| 97精品久久久久久久久久精品| 人人妻人人看人人澡| 男人舔奶头视频| 一级二级三级毛片免费看| 建设人人有责人人尽责人人享有的 | 久久国产乱子免费精品| 六月丁香七月| 精品久久久噜噜| 午夜福利网站1000一区二区三区| 成人漫画全彩无遮挡| 性色av一级| 十八禁网站网址无遮挡 | 久久久亚洲精品成人影院| 国产亚洲最大av| 久久精品国产亚洲网站| 看十八女毛片水多多多| 精品午夜福利在线看| 亚洲欧美中文字幕日韩二区| 成年女人在线观看亚洲视频 | 香蕉精品网在线| 国产欧美亚洲国产| 免费播放大片免费观看视频在线观看| 国产亚洲精品久久久com| 一级爰片在线观看| 2022亚洲国产成人精品| 国精品久久久久久国模美| 成人二区视频| 美女脱内裤让男人舔精品视频| 我的女老师完整版在线观看| 在线免费观看不下载黄p国产| 午夜免费观看性视频| 日韩人妻高清精品专区| 午夜亚洲福利在线播放| 在线观看人妻少妇| 男的添女的下面高潮视频| 丰满乱子伦码专区| 亚洲国产成人一精品久久久| 亚洲第一区二区三区不卡| 久久久久久伊人网av| 欧美日韩视频高清一区二区三区二| 91在线精品国自产拍蜜月| 只有这里有精品99| 男女下面进入的视频免费午夜| 国产免费福利视频在线观看| 国内揄拍国产精品人妻在线| 亚洲国产欧美人成| 欧美xxxx黑人xx丫x性爽| 婷婷色av中文字幕| 99热这里只有是精品50| 啦啦啦啦在线视频资源| 另类亚洲欧美激情| 亚洲av中文av极速乱| 日韩大片免费观看网站| 深爱激情五月婷婷| 日韩免费高清中文字幕av| 久久精品久久久久久久性| 国产综合懂色| 久久鲁丝午夜福利片| 久久久午夜欧美精品| 爱豆传媒免费全集在线观看| 性插视频无遮挡在线免费观看| 青春草国产在线视频| 夜夜爽夜夜爽视频| 久久女婷五月综合色啪小说 | 国产男女内射视频| 久久午夜福利片| 嘟嘟电影网在线观看| 欧美高清成人免费视频www| 日韩人妻高清精品专区| 午夜激情久久久久久久| 国产高清有码在线观看视频| 亚洲婷婷狠狠爱综合网| 91久久精品电影网| 婷婷色综合大香蕉| 国产黄a三级三级三级人| 亚洲综合精品二区| 久久久久久久国产电影| 亚洲,欧美,日韩| 22中文网久久字幕| 大码成人一级视频| 好男人视频免费观看在线| 久热久热在线精品观看| 久久亚洲国产成人精品v| 亚洲一区二区三区欧美精品 | 亚洲av不卡在线观看| 亚洲,欧美,日韩| 婷婷色麻豆天堂久久| 色视频在线一区二区三区| 国产乱人偷精品视频| 国产成人精品久久久久久| 激情五月婷婷亚洲| 国产日韩欧美亚洲二区| 欧美区成人在线视频| xxx大片免费视频| 性色av一级| 99精国产麻豆久久婷婷| 男女国产视频网站| 欧美成人精品欧美一级黄| 插逼视频在线观看| 久久久久久久午夜电影| 精品人妻一区二区三区麻豆| 亚洲成人精品中文字幕电影| 亚洲精品日本国产第一区| 成人欧美大片| 精品一区二区免费观看| 性色av一级| 美女脱内裤让男人舔精品视频| 91久久精品国产一区二区三区| 女的被弄到高潮叫床怎么办| 国产高潮美女av| 日本爱情动作片www.在线观看| 日本一本二区三区精品| 男女下面进入的视频免费午夜| 中文乱码字字幕精品一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 三级国产精品欧美在线观看| 卡戴珊不雅视频在线播放| 国产一区有黄有色的免费视频| 免费看a级黄色片| 18禁在线无遮挡免费观看视频| 联通29元200g的流量卡| 少妇人妻 视频| 亚洲欧美中文字幕日韩二区| 亚洲成人av在线免费| 欧美国产精品一级二级三级 | 欧美3d第一页| av网站免费在线观看视频| 日本爱情动作片www.在线观看| 一级爰片在线观看| 草草在线视频免费看| 在线免费十八禁| 九草在线视频观看| 国产精品爽爽va在线观看网站| 亚洲av在线观看美女高潮| 亚洲精品乱码久久久久久按摩| 欧美zozozo另类| 精品视频人人做人人爽| 久久久成人免费电影| 国产黄色视频一区二区在线观看| 国产熟女欧美一区二区| 欧美激情在线99| 联通29元200g的流量卡| 日日撸夜夜添| 国产爽快片一区二区三区| 日韩欧美 国产精品| 亚洲国产精品成人综合色| 男女无遮挡免费网站观看| 少妇丰满av| 麻豆成人午夜福利视频| 美女脱内裤让男人舔精品视频| 欧美性感艳星| 波多野结衣巨乳人妻| 亚洲欧美成人综合另类久久久| 国产黄色视频一区二区在线观看| 肉色欧美久久久久久久蜜桃 | 精品国产露脸久久av麻豆| 深爱激情五月婷婷| 赤兔流量卡办理| 亚洲av.av天堂| 免费看av在线观看网站| 久久热精品热| 五月天丁香电影| 成年版毛片免费区| 超碰av人人做人人爽久久| 国产精品麻豆人妻色哟哟久久| 久久国产乱子免费精品| 麻豆乱淫一区二区| 国产精品一区二区性色av| 亚洲国产精品成人综合色| 午夜免费男女啪啪视频观看| 亚洲国产欧美在线一区| 新久久久久国产一级毛片| av播播在线观看一区| 国产精品嫩草影院av在线观看| 精品酒店卫生间| 97在线视频观看| 新久久久久国产一级毛片| 欧美日韩视频精品一区| 三级国产精品片| 亚洲电影在线观看av| 嫩草影院新地址| 成人欧美大片| 午夜激情久久久久久久| 国产色爽女视频免费观看| 日韩中字成人| 亚洲精品456在线播放app| 国产爱豆传媒在线观看| 国内精品宾馆在线| 在线天堂最新版资源| 亚洲欧美日韩卡通动漫| 亚洲国产精品999| 国产一区亚洲一区在线观看| 精品人妻一区二区三区麻豆| 亚洲国产av新网站| 晚上一个人看的免费电影| 久久久亚洲精品成人影院| 亚洲va在线va天堂va国产| 久久精品国产亚洲av涩爱| 日韩av免费高清视频| 亚洲国产精品国产精品| 欧美bdsm另类| 成人二区视频| 黄片无遮挡物在线观看| 国产中年淑女户外野战色| 日本-黄色视频高清免费观看| 亚洲欧美日韩卡通动漫| 不卡视频在线观看欧美| 亚洲av电影在线观看一区二区三区 | 欧美亚洲 丝袜 人妻 在线| 免费观看的影片在线观看| 色综合色国产| 丝袜美腿在线中文| 国产一区二区三区综合在线观看 | 天堂俺去俺来也www色官网| 国产成人精品久久久久久| 嫩草影院入口| 少妇高潮的动态图| 免费观看无遮挡的男女| 黄色配什么色好看| 欧美性感艳星| 精品一区二区三卡| 亚洲成人精品中文字幕电影| 丰满少妇做爰视频| 精品亚洲乱码少妇综合久久| 国产精品.久久久| 亚洲国产欧美在线一区| tube8黄色片| 日韩不卡一区二区三区视频在线| 免费看av在线观看网站| 中文欧美无线码| 最近的中文字幕免费完整| 制服丝袜香蕉在线| 日本一二三区视频观看| 全区人妻精品视频| 成人无遮挡网站| 日韩,欧美,国产一区二区三区| 精品一区在线观看国产| 男人添女人高潮全过程视频| 国产欧美另类精品又又久久亚洲欧美| 国产 一区 欧美 日韩| 国产精品秋霞免费鲁丝片| 干丝袜人妻中文字幕| 男女啪啪激烈高潮av片| 黄色视频在线播放观看不卡| 边亲边吃奶的免费视频| 丝袜脚勾引网站| 天堂中文最新版在线下载 | 国产成人freesex在线| 日韩强制内射视频| 亚洲欧美成人精品一区二区| 欧美性猛交╳xxx乱大交人| 免费人成在线观看视频色| 亚洲自偷自拍三级| 国产免费又黄又爽又色| 高清av免费在线| 久久99热这里只有精品18| 国产视频内射| 免费看a级黄色片| 少妇人妻精品综合一区二区| 亚洲精品日韩在线中文字幕| 少妇熟女欧美另类| 2021天堂中文幕一二区在线观| 狂野欧美白嫩少妇大欣赏| 国产免费一级a男人的天堂| 在线观看免费高清a一片| 在线精品无人区一区二区三 | 国产一级毛片在线| 亚洲天堂av无毛| 干丝袜人妻中文字幕| 一区二区三区精品91| 亚洲精品第二区| 日韩一本色道免费dvd| 老司机影院成人| 全区人妻精品视频| videossex国产| 亚洲经典国产精华液单| 能在线免费看毛片的网站| 久久久久久久久久人人人人人人| 王馨瑶露胸无遮挡在线观看| 大陆偷拍与自拍| av国产久精品久网站免费入址| .国产精品久久| 在线观看免费高清a一片| 黄色一级大片看看| 国产成人免费无遮挡视频| 国产精品一区二区三区四区免费观看| 色播亚洲综合网| 亚洲色图av天堂| 69av精品久久久久久| 人妻制服诱惑在线中文字幕| 免费观看无遮挡的男女| 国产男女超爽视频在线观看| 亚洲精品国产av成人精品| 欧美区成人在线视频| 一级毛片我不卡| 插逼视频在线观看| 老司机影院成人| 永久免费av网站大全| 色婷婷久久久亚洲欧美| 男女啪啪激烈高潮av片| 国产亚洲一区二区精品| 亚洲自拍偷在线| 国产69精品久久久久777片| 欧美亚洲 丝袜 人妻 在线| 麻豆精品久久久久久蜜桃| 亚洲精品,欧美精品| 性色av一级| 欧美日韩在线观看h| 色视频www国产| 欧美日韩在线观看h| 男人狂女人下面高潮的视频| 看非洲黑人一级黄片| 久久久久久久大尺度免费视频| 亚洲精品日韩av片在线观看| 高清毛片免费看| 亚洲欧美成人综合另类久久久| 别揉我奶头 嗯啊视频| 免费av毛片视频| 在线观看人妻少妇| 日韩av免费高清视频| 人妻 亚洲 视频| 久久ye,这里只有精品| 亚洲国产精品999| 少妇裸体淫交视频免费看高清| 亚洲,欧美,日韩| 国产高清有码在线观看视频| 国产男女超爽视频在线观看| 国产精品嫩草影院av在线观看| 交换朋友夫妻互换小说| 三级国产精品片| 亚洲欧美日韩无卡精品| 蜜桃久久精品国产亚洲av| 久久久成人免费电影| 免费黄频网站在线观看国产| 麻豆国产97在线/欧美| 看黄色毛片网站| a级毛片免费高清观看在线播放| 波多野结衣巨乳人妻| av在线天堂中文字幕| 伦精品一区二区三区| 女人十人毛片免费观看3o分钟| 看十八女毛片水多多多| 久久国内精品自在自线图片| 中文天堂在线官网| 久久精品久久久久久久性| 美女主播在线视频| 亚洲天堂国产精品一区在线| 九色成人免费人妻av| 久久久精品免费免费高清| 国产精品精品国产色婷婷| 伦精品一区二区三区| 69av精品久久久久久| 欧美日韩综合久久久久久| 亚洲av免费高清在线观看| 边亲边吃奶的免费视频| 美女高潮的动态| 欧美bdsm另类| 国产精品麻豆人妻色哟哟久久| 九色成人免费人妻av| kizo精华| 国产免费视频播放在线视频| 天堂中文最新版在线下载 | 成年女人看的毛片在线观看| 少妇 在线观看| 黄片无遮挡物在线观看| 直男gayav资源| 黑人高潮一二区| 成人午夜精彩视频在线观看| 超碰av人人做人人爽久久| 五月开心婷婷网| 精品久久久久久久末码| 边亲边吃奶的免费视频| 成人鲁丝片一二三区免费| 男女下面进入的视频免费午夜| 中文在线观看免费www的网站| 亚洲av福利一区| 看黄色毛片网站| 国产黄片视频在线免费观看| 成人亚洲欧美一区二区av| 大香蕉97超碰在线| 免费高清在线观看视频在线观看| 成年av动漫网址| 国产欧美另类精品又又久久亚洲欧美| 国产精品人妻久久久久久| 91在线精品国自产拍蜜月| 99精国产麻豆久久婷婷| 国产美女午夜福利| 欧美日韩精品成人综合77777| av在线app专区| 精品久久久久久久久av| 伊人久久精品亚洲午夜| 日本一二三区视频观看| 汤姆久久久久久久影院中文字幕| 蜜桃久久精品国产亚洲av| 亚洲自偷自拍三级| av在线亚洲专区| 丝袜美腿在线中文| 身体一侧抽搐| 成人二区视频| 97超视频在线观看视频| 另类亚洲欧美激情| 卡戴珊不雅视频在线播放| 不卡视频在线观看欧美| 性插视频无遮挡在线免费观看| 永久网站在线| 热re99久久精品国产66热6| 草草在线视频免费看| 国产探花在线观看一区二区| 成人漫画全彩无遮挡| 亚洲成色77777| 国产淫片久久久久久久久| 亚洲真实伦在线观看| 一级av片app| av在线观看视频网站免费| 国产成人aa在线观看| 黄色日韩在线| 久久97久久精品| www.av在线官网国产| 一级毛片我不卡| 日韩av在线免费看完整版不卡| 亚洲精品456在线播放app| 一级毛片 在线播放| 少妇的逼水好多| av国产精品久久久久影院| 精品熟女少妇av免费看| 成年免费大片在线观看| 亚洲伊人久久精品综合| 另类亚洲欧美激情| 大片免费播放器 马上看| 成人毛片a级毛片在线播放| 美女视频免费永久观看网站| 18禁裸乳无遮挡动漫免费视频 | 69av精品久久久久久| 麻豆久久精品国产亚洲av| 成人鲁丝片一二三区免费| 国产成人免费观看mmmm| 欧美高清成人免费视频www| 天天躁日日操中文字幕| 久久精品人妻少妇| 高清毛片免费看| 国产亚洲午夜精品一区二区久久 | 国产老妇伦熟女老妇高清| 国产精品偷伦视频观看了| av在线app专区| 久久ye,这里只有精品| 欧美最新免费一区二区三区| 久久精品夜色国产| 精品久久久久久电影网| 80岁老熟妇乱子伦牲交| 亚洲婷婷狠狠爱综合网| 国产黄色免费在线视频| 日日啪夜夜爽| 一级二级三级毛片免费看| 精品国产露脸久久av麻豆| 国产精品国产三级国产专区5o| 亚洲美女视频黄频| 国产午夜精品一二区理论片| 王馨瑶露胸无遮挡在线观看| 午夜福利在线在线| 在线观看av片永久免费下载| 成年av动漫网址| 久久久久久国产a免费观看| 精品酒店卫生间| 亚洲成人久久爱视频| 亚洲精品日韩av片在线观看| 午夜精品一区二区三区免费看| 校园人妻丝袜中文字幕| 身体一侧抽搐| 中文字幕亚洲精品专区| 男女边吃奶边做爰视频| 免费av不卡在线播放| 久久久久网色| 婷婷色综合大香蕉| 国产欧美日韩精品一区二区| 少妇猛男粗大的猛烈进出视频 | 国产91av在线免费观看| 久久久久久久午夜电影| 在线a可以看的网站| 久久精品夜色国产| 美女国产视频在线观看| 亚洲性久久影院| 亚洲最大成人手机在线| 欧美日韩精品成人综合77777| 联通29元200g的流量卡| 一个人观看的视频www高清免费观看| 91久久精品电影网| 亚洲图色成人| 久久影院123| 精品熟女少妇av免费看| 欧美成人精品欧美一级黄| 有码 亚洲区| 久久鲁丝午夜福利片| 在线 av 中文字幕| 日韩精品有码人妻一区| 免费av观看视频| 美女cb高潮喷水在线观看| 亚洲天堂av无毛| av专区在线播放| av国产久精品久网站免费入址| 自拍偷自拍亚洲精品老妇| 国产一区有黄有色的免费视频| 日本爱情动作片www.在线观看| 日韩欧美 国产精品| 99久久人妻综合| 美女xxoo啪啪120秒动态图| 天堂中文最新版在线下载 | 一区二区三区精品91| 国产日韩欧美在线精品| 亚洲国产精品成人综合色| 久久久午夜欧美精品| 国产精品麻豆人妻色哟哟久久| 亚洲欧美日韩卡通动漫| freevideosex欧美| 国产成人a区在线观看| 丝袜美腿在线中文| 日日啪夜夜撸| 精品久久久噜噜| 毛片一级片免费看久久久久| 久久久午夜欧美精品| 国产一区有黄有色的免费视频| 免费观看无遮挡的男女| 日日啪夜夜撸| 午夜精品一区二区三区免费看| 一个人观看的视频www高清免费观看| 亚洲精品久久久久久婷婷小说| 一级爰片在线观看| 91精品伊人久久大香线蕉| 神马国产精品三级电影在线观看| 国产色爽女视频免费观看| 一二三四中文在线观看免费高清| 汤姆久久久久久久影院中文字幕| 国产日韩欧美在线精品| 天天躁日日操中文字幕| 91aial.com中文字幕在线观看| 国产探花在线观看一区二区| 欧美xxxx黑人xx丫x性爽| 亚洲熟女精品中文字幕| 国产欧美日韩精品一区二区| 一级毛片我不卡| 97人妻精品一区二区三区麻豆| 久热这里只有精品99| 午夜老司机福利剧场| av免费观看日本|