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

    數(shù)值法在卡莫亞礦區(qū)水源地地下水補(bǔ)給資源量計(jì)算中的應(yīng)用

    2017-08-28 22:49:13李淼清李貴仁
    地下水 2017年4期
    關(guān)鍵詞:資源量水源地水文地質(zhì)

    李淼清,鄔 立,李貴仁

    (華北有色工程勘察院有限公司,河北 石家莊 050021)

    數(shù)值法在卡莫亞礦區(qū)水源地地下水補(bǔ)給資源量計(jì)算中的應(yīng)用

    李淼清,鄔 立,李貴仁

    (華北有色工程勘察院有限公司,河北 石家莊 050021)

    卡莫亞礦區(qū)位于剛果(金)東南部的加丹加高原,區(qū)域水資源比較匱乏,為解決礦山供水問題,尋找滿足礦山正常生產(chǎn)用水的供水水源地,在對礦區(qū)水源地水文地質(zhì)條件綜合分析的基礎(chǔ)上,分別采用斷面法及數(shù)值法對供水水源地地下水資源進(jìn)行評價(jià)。重點(diǎn)對數(shù)值法進(jìn)行詳細(xì)論述,通過對邊界條件、含水層及地下水流動(dòng)狀態(tài)進(jìn)行概化,建立水源地的水文地質(zhì)概念模型,運(yùn)用Visual MODFLOW數(shù)值模擬軟件建立地下水?dāng)?shù)值模擬模型,通過均衡分析計(jì)算,最終計(jì)算獲得區(qū)內(nèi)地下水補(bǔ)給資源量為13 551.89 m3/d。與采用斷面法計(jì)算的結(jié)果13 328.81 m3/d相比較,兩種方法計(jì)算結(jié)果基本吻合,而數(shù)值法在計(jì)算過程紅充分運(yùn)用了群孔抽水資料進(jìn)行模型識(shí)別,結(jié)果相對合理,因此計(jì)算結(jié)果以數(shù)值法為準(zhǔn)。

    Visual MODFOW;地下水補(bǔ)給資源量;卡莫亞礦區(qū);水均衡分析

    卡莫亞礦區(qū)位于剛果(金)東南部的加丹加高原,由于水資源比較匱乏,礦區(qū)的正常生產(chǎn)活動(dòng)已受到一定的影響,為徹底解決礦山的供水問題,擬在礦區(qū)范圍內(nèi)尋找滿足礦山正常生產(chǎn)用水的供水水源地[1]。本文在對水源地的地質(zhì)、水文地質(zhì)條件進(jìn)行綜合分析的基礎(chǔ)上,采用斷面法及數(shù)值法對供水水源地的地下水資源進(jìn)行了評價(jià),計(jì)算出水源地的補(bǔ)給資源量,為礦區(qū)未來的供水方向提供了依據(jù)。數(shù)值法采用Visual MODFLOW軟件進(jìn)行三維數(shù)值模型的建立,該軟件是在整合MODFLOW軟件模塊的基礎(chǔ)上應(yīng)用可視化技術(shù)開發(fā)研制,無縫集成了MODFLOW-2000、WinPEST、MT3D、MODPATH等軟件[2-4]。

    1 研究區(qū)概況

    卡莫亞銅鈷礦區(qū)位于剛果(金)加丹加省坎博韋鎮(zhèn),與著名的“坎博韋銅鈷礦選廠”相鄰,交通較為便利。礦區(qū)至利卡西(Likasi)市區(qū)距離約25 km,由碎石公路相連;礦區(qū)南東方向距離加丹加省省會(huì)盧本巴希市(Lubumbashi)約130 km,利卡西市至盧本巴希市有100 km較好的公路。

    工作區(qū)處于剛果(金)東南部的加丹加高原,海拔標(biāo)高在1 375~1 560 m之間,地理環(huán)境隸屬加丹加高原上的低矮丘陵區(qū),山巒起伏不大,總體呈南西高、北東低之勢。礦區(qū)所在地屬熱帶草原氣候,年均氣溫為27℃,氣候溫和。按年降雨量分旱、雨兩季:5月至9月為旱季,降雨量極少;10月至翌年4月為雨季,大雨乃至暴雨不斷。

    區(qū)內(nèi)出露的地層主要是第四系(Q)及中元古界加丹加群的孔德龍古組(K)、羅安組(R),由新到老依次為第四系、孔德龍古組(K)、羅安組(R)。第四系主要為沖洪積及殘坡積物,巖性多為粘性土、角礫以及砂礫石等,多分布在地勢低洼的溝谷地帶、坡麓及平緩的山坡等處,厚度一般0~12 m左右??椎慢埞沤M(K)地層為區(qū)域范圍出露的主要地層,在礦區(qū)的北側(cè)大面積出露,另外坎博韋推覆體下亦為本組地層,出露的巖性主要為砂礫巖及淺紫紅色條帶狀粉砂巖。羅安組(R)是本區(qū)及區(qū)域最主要的含礦層位,為一套碎屑巖-泥質(zhì)巖-碳酸鹽巖組合的沉積建造。

    礦區(qū)內(nèi)分布著多個(gè)含水層,巖性結(jié)構(gòu)差異導(dǎo)致含水層透水性強(qiáng)弱不均一。在垂向上,含水層強(qiáng)弱相間,但由于區(qū)內(nèi)構(gòu)造發(fā)育,巖層破碎,弱含水層也具有一定的透水性,各含水層之間存在水力聯(lián)系,含水層間水位未發(fā)生明顯變化,含水層整體為潛水含水層,各含水層構(gòu)成一個(gè)統(tǒng)一的地下含水系統(tǒng)。從礦區(qū)整體空間上來看,表現(xiàn)為內(nèi)強(qiáng)外弱的總體特征,大氣降水入滲補(bǔ)給為地下水重要充水來源,內(nèi)部R2礦山段為地下水儲(chǔ)存和運(yùn)移通道,其中R2.3(CMN)巖溶裂隙含水層、R2.1(RSC、RSF)巖溶裂隙含水層為主要含水層。

    2 水文地質(zhì)概念模型

    2.1 計(jì)算區(qū)域

    計(jì)算區(qū)邊界基本與卡莫亞和坎博韋兩大礦區(qū)邊界吻合,具體為西部邊界為西南礦體、西礦體外緣,接近西側(cè)分水嶺;南部邊界為兩礦區(qū)南側(cè)分水嶺;東部邊界為坎博韋礦東緣,處于東側(cè)分水嶺上;北部邊界為兩礦羅安組與孔德龍古組地層交界處,地勢向東北逐漸降低;整個(gè)計(jì)算區(qū)為區(qū)域內(nèi)一個(gè)半封閉的次級水文地質(zhì)單元。

    2.2 含水層概化

    區(qū)內(nèi)含水層巖組具有很強(qiáng)的不均一性,垂向和水平向上強(qiáng)弱含水層呈相間分布。區(qū)內(nèi)主要含水層組為羅安組第二亞組下段R2.1、上段R2.3、羅安組第四亞組下段R4.1以及部分構(gòu)造破碎帶Br,巖性主要為白云巖、硅化白云巖;弱含水層主要為羅安組第二亞組中段R2.2、羅安組第一亞組R1、羅安組第四亞組上段R4.2、上孔德龍古組下段Ks1.3及部分構(gòu)造破碎帶,巖性主要為頁巖、砂巖、粉砂巖、礫巖、完整的白云巖;各含水層之間存在水力聯(lián)系,總體為潛水含水層,局部承壓。綜上所述,將含水層概化為非均質(zhì)潛水含水層。

    2.3 邊界條件概化

    卡莫亞礦區(qū)西部、北部邊界為上孔德龍古組弱含水層,概化為弱透水Ⅱ類邊界;卡莫亞礦南部和坎博韋礦東部邊界恰位于分水嶺上,概化為零流量邊界;上邊界為大氣降水補(bǔ)給邊界,下邊界為上孔德龍古組弱含水層,概化為弱透水Ⅱ類邊界。

    3 斷面法計(jì)算地下水補(bǔ)給資源量

    3.1 均衡區(qū)范圍

    東界為坎博韋鎮(zhèn)東側(cè)地表分水嶺附近,南界為礦區(qū)南部地表分水嶺附近山脊一線,西界為西南礦體西南側(cè)山脊-西礦體山脊-生活營區(qū)一線,北界為北礦體北側(cè)邊界孔德龍古組地層一線;面積約12.18 km2。

    3.2 地下水補(bǔ)給量計(jì)算

    在現(xiàn)狀條件下勘探區(qū)內(nèi)地下水的補(bǔ)給來源有:南側(cè)、東側(cè)為地表分水嶺,地下水無側(cè)向補(bǔ)給,北側(cè)、西側(cè)為透水性較弱的孔德龍古組弱含水層,側(cè)向補(bǔ)給量較小,勘探區(qū)內(nèi)地下水的主要補(bǔ)給來源為大氣降水的垂向入滲補(bǔ)給。

    計(jì)算公式:

    Q補(bǔ)=Q降+Q東+Q西+Q南+Q北

    (1)

    式中:Q補(bǔ)為地下水補(bǔ)給資源量(m3/d);Q降為區(qū)內(nèi)降水垂向補(bǔ)給資源量(m3/d);Q東為東側(cè)側(cè)向徑流補(bǔ)給資源量(m3/d);Q南為南側(cè)側(cè)向徑流補(bǔ)給資源量(m3/d);Q西為西側(cè)側(cè)向徑流補(bǔ)給資源量(m3/d);Q北為北側(cè)側(cè)向徑流補(bǔ)給資源量(m3/d)。由于東側(cè)和南側(cè)均為地表分水嶺,無地下水側(cè)向補(bǔ)給資源量,因此

    Q東=0,Q南=0。

    區(qū)內(nèi)西側(cè)、北側(cè)為孔德龍古組地下水補(bǔ)給邊界,根據(jù)Darcy公式計(jì)算其側(cè)向排泄量,選擇礦區(qū)東部地下水等水位線圖上的1390m地下水位等勢線作為計(jì)算斷面,地下水徑流量按下式計(jì)算:

    (2)

    式中:Q徑為通過計(jì)算斷面的天然地下水徑流量(m3/d);Ki為第i 個(gè)斷面含水層滲透系數(shù)(m/d);Ii為第i個(gè)斷面地下水平均水力坡度(%);Mi為第i個(gè)斷面含水層厚度(m);Li為第i個(gè)斷面長度(m)。

    由此可計(jì)算出,Q北=2 606.67 m3/d,Q西=673.64 m3/d。

    勘探區(qū)內(nèi)降水量較大,且區(qū)內(nèi)沖溝較發(fā)育,多見羅安組含水層出露地表,溶蝕裂隙發(fā)育,有利于降水入滲補(bǔ)給下部含水層,因此其降水補(bǔ)給資源量計(jì)算如下:

    Q降=A*α*W

    (3)

    式中:Q降為地下水補(bǔ)給資源量(m3/d);A為均衡區(qū)面積,取值12.18×106m2;α為降水入滲系數(shù),參考剛果(金)其他類似礦區(qū),取值0.25;W為區(qū)內(nèi)日平均降水量,區(qū)內(nèi)年平均降水量為1 197.9 mm,則日平均降水量為0.003 3 m;由此可計(jì)算出Q降=10 048.50 m3/d

    綜上可計(jì)算出勘探區(qū)地下水補(bǔ)給量為:Q補(bǔ)=13 328.81 m3/d。

    4 數(shù)值法計(jì)算地下水補(bǔ)給資源量

    4.1 數(shù)學(xué)模型

    根據(jù)前述水文地質(zhì)條件分析,數(shù)學(xué)模型為:

    (4)

    4.2 數(shù)值模型

    4.2.1 空間離散

    采用有限差分法對物理模型進(jìn)行矩形剖分,模型范圍(261250,290750)~(266250,293500),對地層分界線、抽水孔、觀測孔附近進(jìn)行加密處理,共剖分成78行、158列、25層。共計(jì)78×158×25=308 100個(gè)單元格[5]。

    4.2.2 水文地質(zhì)參數(shù)的選取

    通過對研究區(qū)地層巖性、地質(zhì)構(gòu)造、鉆孔巖芯采取率、抽水試驗(yàn)以及地下水水位恢復(fù)試驗(yàn)等資料的分析,將模擬區(qū)劃分為24個(gè)水文地質(zhì)參數(shù)分區(qū)(表1、圖1、圖2)。

    4.2.3 初始流場

    根據(jù)群孔抽水試驗(yàn)前的鉆孔地下水位觀測資料,繪制初始流場。

    表1 數(shù)值模擬水文地質(zhì)參數(shù)分區(qū)表

    圖1 滲透性和給水度參數(shù)分區(qū)圖(1 336 m高程)

    4.2.4 模型識(shí)別

    利用模型反演求參的方法,對模型的初始流場、水文地質(zhì)參數(shù)初始值和邊界條件進(jìn)行識(shí)別,識(shí)別的主要依據(jù)是將抽水試驗(yàn)階段各觀測孔水位觀測數(shù)據(jù)作為目標(biāo)函數(shù),進(jìn)行參數(shù)的反演識(shí)別。經(jīng)過反復(fù)調(diào)整初始參數(shù),通過對各觀測孔水位模擬值與地下水流場實(shí)際觀測值的擬合,獲得各區(qū)較為合理的水文地質(zhì)參數(shù),作為地下水天然資源量計(jì)算和未來水源地運(yùn)行預(yù)測的基礎(chǔ)。經(jīng)過反復(fù)調(diào)試,地下水位動(dòng)態(tài)曲線擬合結(jié)果較好,令人較為滿意。擬合出的較合理的水文地質(zhì)參數(shù)如表2所示,地下水位動(dòng)態(tài)曲線擬合結(jié)果見圖3、圖4。

    圖2 滲透性和給水度參數(shù)分區(qū)圖(Y=291 905剖面)

    圖3 西區(qū)各孔水位擬合曲線圖

    圖4 東區(qū)各孔水位擬合曲線圖

    ZOneKx[m/d]Ky[m/d]Kz[m/d]Sy1234567891011121314151617181920212223240.0150.0150.00150.10110.10.10550.50.15220.10.15550.50.20110.10.102.52.50.250.10101010.200.010.010.010.053.53.50.350.100.10.10.010.101.51.50.150.100.0010.0010.0070.05220.20.15220.20.15880.80.15110.10.150.50.50.050.050.50.50.050.050.10.10.010.10440.40.10440.40.05110.10.10110.10.05

    4.3 地下水補(bǔ)給量計(jì)算

    在模型中建立均衡區(qū),均衡區(qū)范圍為模型計(jì)算區(qū)域,將模型中DCK01、DCK02井抽水量設(shè)置為零,計(jì)算現(xiàn)有狀態(tài)下均衡區(qū)內(nèi)地下水補(bǔ)給資源量。經(jīng)計(jì)算,地下水補(bǔ)給資源量為13 551.89 m3/d,天然狀態(tài)下地下水資源均衡量差為-5 648.11 m3/d,已處于超采狀態(tài)。均衡區(qū)內(nèi)地下水資源量如表3所示。

    表3 水源地天然狀態(tài)下地下水資源均衡量表 m3/d

    5 結(jié)語

    (1)在水文地質(zhì)條件分析的基礎(chǔ)上,建立礦區(qū)水文地質(zhì)概念模型,模型概化合理,邊界清晰,構(gòu)成了一個(gè)相對獨(dú)立的水文地質(zhì)單元。單元內(nèi)主要含水層為R2.3、R2.1和R4巖溶裂隙含水層組,含水層間具有較好的水力聯(lián)系,地下水運(yùn)動(dòng)以水平運(yùn)動(dòng)為主。地下水補(bǔ)給以大氣降水入滲補(bǔ)給為主,側(cè)向徑流補(bǔ)給為輔,地下水排泄主要為人工開采。

    (2)地下水資源計(jì)算以斷面法為基礎(chǔ),以數(shù)值法為重點(diǎn),不同方法的計(jì)算結(jié)果相互對比、驗(yàn)證,更能客觀地反映地下水資源量。按斷面法計(jì)算,地下水補(bǔ)給資源量為13 328.81 m3/d;按數(shù)值法計(jì)算,地下水補(bǔ)給資源量為13 551.89 m3/d。兩種方法計(jì)算結(jié)果基本吻合,而數(shù)值法在計(jì)算過程紅充分運(yùn)用群孔抽水資料進(jìn)行模型識(shí)別,結(jié)果相對合理,因此計(jì)算結(jié)果以數(shù)值法為準(zhǔn),水源地地下水補(bǔ)給資源量為13 551.89 m3/d。

    [1]華北有色工程勘察院有限公司.剛果(金)加丹加省卡莫亞銅鈷礦區(qū)供水水文地質(zhì)勘探報(bào)告[R].石家莊:華北有色工程勘察院有限公司.2014.11.

    [2]李貴仁,趙珍,陳植華.?dāng)?shù)值模擬在反演礦區(qū)水文地質(zhì)條件中的應(yīng)用[J].水文地質(zhì)工程地質(zhì).2013.40(2):19-23.

    [3]王慶永,賈忠華.Visual MODFLOW及其在地下水模擬中的應(yīng)用[J].水資源與水工程學(xué)報(bào).2007.18(5):90-92.

    [4]李貴仁,趙珍,陳植華.復(fù)雜巖溶礦區(qū)疏干條件下的地下水?dāng)?shù)值模擬—以福建省馬坑鐵礦為例[J].中國巖溶.2012.31(4):382-387.

    [5]馮更辰,郝俊杰,譚?。甐isual MODFLOW模型在白澗鐵礦區(qū)礦井涌水量預(yù)測中的應(yīng)用[J].中國巖溶.2011.30(3):271-276.

    Application of Numerical Method in Calculating the Amount of Groundwater Recharge Resources at Water Source Area of Mining Area

    LI Miao-qing,WU Li,LI Gui-ren

    (Huabei Nonferrous Engineering Investigation Institute Co Ltd,Hebei Shijiazhuang 050021)

    The Karmo mining area is located in the Katanga Plateau in the southeastern part of the Congo (DRC), and water resources are scarce in this region. In order to solve the problem of water supply, and look for the source area for water supply which can satisfy the normal usage of water, based on the comprehensive analysis of water and geological conditions in this region, the groundwater resources at the water supply area are evaluated by cross-section method and numerical method respectively. This paper focuses on the numerical method, the concept model of hydrogeological geology is established through the generalization of boundary conditions, aquifer and flowing status of groundwater, and then uses Visual MODFLOW numerical simulation software to establish the groundwater numerical simulation model. After the equilibrium analysis and calculation, the amount of groundwater recharge resources in the area is 13 551.89 m3/d. Compared with the results calculated by the cross-section method 13 328.81 m3/ d, the results of the two methods are basically consistent, but numerical method fully uses the hole pumping data for model identification in the calculation process, the result is relatively reasonable, so the results should be from the numerical method.

    Visual MODFOW;the amount of groundwater recharge resources; the Karmo mining area; water equilibrium analysis

    2017-04-07

    李淼清(1978-),男,湖南平江人,高級工程師,主要從事水文地質(zhì)方面工作。

    P641.4+3

    B

    1004-1184(2017)04-0036-04

    猜你喜歡
    資源量水源地水文地質(zhì)
    江埡庫區(qū)魚類群落組成和資源量評估
    基于抽水試驗(yàn)計(jì)算水文地質(zhì)參數(shù)的分析與確定
    基于GPRS實(shí)現(xiàn)自動(dòng)化水文地質(zhì)抽水試驗(yàn)
    鈾礦數(shù)字勘查資源量估算方法應(yīng)用與驗(yàn)證
    基于轉(zhuǎn)移概率的三維水文地質(zhì)結(jié)構(gòu)劃分
    塞拉利昂通戈金剛石礦資源量上升
    水文地質(zhì)在工程地質(zhì)勘察中的有效應(yīng)用
    基于因果分析圖法的飲用水源地保護(hù)探討
    南水北調(diào)中線水源地土壤侵蝕經(jīng)濟(jì)損失估算
    吉林省供水水源地問題探討
    男人舔女人的私密视频| 午夜福利,免费看| 岛国毛片在线播放| 亚洲av中文av极速乱| 精品一区二区三区四区五区乱码 | 国产精品一区二区在线观看99| av有码第一页| 午夜福利,免费看| 国产精品久久久久久av不卡| 啦啦啦中文免费视频观看日本| 制服人妻中文乱码| 中文字幕人妻丝袜制服| 午夜福利在线免费观看网站| av有码第一页| 久久久久国产网址| 18+在线观看网站| 亚洲男人天堂网一区| 亚洲,欧美精品.| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91精品三级在线观看| 国产精品不卡视频一区二区| 国产成人精品福利久久| 肉色欧美久久久久久久蜜桃| 国产黄色免费在线视频| 亚洲av成人精品一二三区| 高清av免费在线| 免费人妻精品一区二区三区视频| 日本欧美视频一区| 精品人妻在线不人妻| 女人高潮潮喷娇喘18禁视频| 国产精品久久久av美女十八| 精品一品国产午夜福利视频| 亚洲欧美成人精品一区二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产高清不卡午夜福利| 国产男女超爽视频在线观看| 欧美xxⅹ黑人| 国产黄色免费在线视频| 精品国产乱码久久久久久小说| 色吧在线观看| 青春草亚洲视频在线观看| 亚洲成人手机| 国产老妇伦熟女老妇高清| 卡戴珊不雅视频在线播放| 成人毛片60女人毛片免费| av.在线天堂| 美女中出高潮动态图| 国产毛片在线视频| 99re6热这里在线精品视频| 国产伦理片在线播放av一区| 在线精品无人区一区二区三| 国产精品一区二区在线观看99| 久久精品国产亚洲av高清一级| 亚洲精品自拍成人| 99香蕉大伊视频| 色吧在线观看| 国产片内射在线| 天美传媒精品一区二区| 9191精品国产免费久久| 日韩 亚洲 欧美在线| 久热久热在线精品观看| 午夜影院在线不卡| 亚洲第一av免费看| 成人毛片a级毛片在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 欧美国产精品va在线观看不卡| 国产又爽黄色视频| 午夜日本视频在线| 亚洲,一卡二卡三卡| 中文字幕制服av| 捣出白浆h1v1| 亚洲综合色惰| 成年人午夜在线观看视频| 欧美97在线视频| 人成视频在线观看免费观看| 亚洲国产精品一区三区| 国产精品蜜桃在线观看| 免费高清在线观看日韩| 大码成人一级视频| 欧美日韩av久久| 80岁老熟妇乱子伦牲交| 一区在线观看完整版| 91在线精品国自产拍蜜月| 久热久热在线精品观看| 亚洲人成77777在线视频| 极品人妻少妇av视频| 日韩欧美精品免费久久| 久久狼人影院| 免费高清在线观看视频在线观看| 高清不卡的av网站| 波多野结衣av一区二区av| 亚洲人成77777在线视频| 99久久中文字幕三级久久日本| 亚洲精品国产av成人精品| 国产精品久久久av美女十八| 2022亚洲国产成人精品| 啦啦啦在线免费观看视频4| 免费黄网站久久成人精品| 久久久久久久久久久免费av| 久久久久视频综合| 丝袜脚勾引网站| 国产精品国产av在线观看| 大香蕉久久成人网| 男男h啪啪无遮挡| 精品人妻一区二区三区麻豆| 亚洲激情五月婷婷啪啪| 青青草视频在线视频观看| 老汉色∧v一级毛片| 国产欧美亚洲国产| 久久久精品国产亚洲av高清涩受| 一本大道久久a久久精品| 久久精品熟女亚洲av麻豆精品| 大话2 男鬼变身卡| 久久综合国产亚洲精品| 中文字幕制服av| 99国产综合亚洲精品| 亚洲一区二区三区欧美精品| 亚洲一级一片aⅴ在线观看| 男女国产视频网站| 国产成人免费无遮挡视频| 久久国内精品自在自线图片| 天堂俺去俺来也www色官网| 丝瓜视频免费看黄片| 国产av一区二区精品久久| 女人久久www免费人成看片| 亚洲av欧美aⅴ国产| 日本wwww免费看| 在线观看国产h片| 热99国产精品久久久久久7| 女人被躁到高潮嗷嗷叫费观| 熟女av电影| 电影成人av| 日韩中文字幕视频在线看片| 亚洲欧美中文字幕日韩二区| 岛国毛片在线播放| 国产亚洲午夜精品一区二区久久| 狠狠精品人妻久久久久久综合| 久久精品人人爽人人爽视色| 一本—道久久a久久精品蜜桃钙片| 亚洲一区中文字幕在线| 99国产精品免费福利视频| 又大又黄又爽视频免费| 9色porny在线观看| 国产精品女同一区二区软件| 免费观看在线日韩| 久久久久久久久久久久大奶| 中文欧美无线码| 亚洲精品国产一区二区精华液| av卡一久久| 如何舔出高潮| 亚洲国产精品999| 最近最新中文字幕大全免费视频 | 久久久久久免费高清国产稀缺| 久久这里有精品视频免费| 免费观看av网站的网址| 国产一区二区在线观看av| 日韩成人av中文字幕在线观看| 麻豆乱淫一区二区| 午夜福利一区二区在线看| 国产成人精品久久久久久| 一级毛片电影观看| 97在线人人人人妻| 欧美精品一区二区大全| 丝袜美足系列| 人妻一区二区av| 少妇被粗大猛烈的视频| 曰老女人黄片| 夫妻性生交免费视频一级片| 制服人妻中文乱码| 97在线人人人人妻| 制服诱惑二区| 成年人午夜在线观看视频| 伦精品一区二区三区| 日日撸夜夜添| 人妻系列 视频| 久久女婷五月综合色啪小说| 中国国产av一级| 亚洲av在线观看美女高潮| 母亲3免费完整高清在线观看 | 日韩一区二区视频免费看| 超色免费av| av国产久精品久网站免费入址| 亚洲精品久久成人aⅴ小说| 国产 一区精品| 国产在视频线精品| 熟女av电影| 亚洲综合色惰| 香蕉精品网在线| 欧美激情高清一区二区三区 | 国产成人aa在线观看| 91精品三级在线观看| 国产97色在线日韩免费| 女的被弄到高潮叫床怎么办| freevideosex欧美| 日本欧美国产在线视频| videosex国产| 成人18禁高潮啪啪吃奶动态图| 亚洲伊人色综图| 永久网站在线| 亚洲国产欧美网| 人妻 亚洲 视频| 免费在线观看黄色视频的| 久久久久网色| 国产精品一区二区在线观看99| 欧美国产精品一级二级三级| tube8黄色片| 久久99蜜桃精品久久| 精品国产国语对白av| av视频免费观看在线观看| 免费播放大片免费观看视频在线观看| 99热全是精品| 制服丝袜香蕉在线| 男男h啪啪无遮挡| 少妇的逼水好多| 日本wwww免费看| 99热全是精品| 精品福利永久在线观看| 国精品久久久久久国模美| 午夜福利乱码中文字幕| 18禁裸乳无遮挡动漫免费视频| 精品国产超薄肉色丝袜足j| 欧美激情高清一区二区三区 | 久久久久久久大尺度免费视频| 满18在线观看网站| 国产成人免费无遮挡视频| 久久精品aⅴ一区二区三区四区 | 中文字幕制服av| 纵有疾风起免费观看全集完整版| av卡一久久| 少妇的逼水好多| 久久久久久伊人网av| 韩国高清视频一区二区三区| 久久免费观看电影| 欧美人与性动交α欧美软件| 女性生殖器流出的白浆| 午夜免费鲁丝| 亚洲精品一二三| 各种免费的搞黄视频| 国产精品久久久久久精品古装| 亚洲一码二码三码区别大吗| 日韩欧美精品免费久久| 成人免费观看视频高清| 久久久久久久精品精品| 波多野结衣av一区二区av| 永久网站在线| 一区二区三区精品91| 热99久久久久精品小说推荐| 亚洲av日韩在线播放| 国产精品秋霞免费鲁丝片| 国产熟女欧美一区二区| 精品午夜福利在线看| 亚洲欧美精品综合一区二区三区 | 你懂的网址亚洲精品在线观看| 色吧在线观看| av不卡在线播放| 国产av一区二区精品久久| 街头女战士在线观看网站| 久久精品久久精品一区二区三区| 成人二区视频| 亚洲成色77777| 自拍欧美九色日韩亚洲蝌蚪91| 少妇人妻精品综合一区二区| 91精品三级在线观看| 国产欧美日韩综合在线一区二区| 精品少妇久久久久久888优播| 丝袜在线中文字幕| 免费观看在线日韩| 人人妻人人添人人爽欧美一区卜| 日韩在线高清观看一区二区三区| 日韩av免费高清视频| 国产在线一区二区三区精| 超色免费av| 国产片特级美女逼逼视频| 黄色配什么色好看| 人人澡人人妻人| 国产精品香港三级国产av潘金莲 | 国产精品偷伦视频观看了| 久久国内精品自在自线图片| 熟女电影av网| 一级毛片我不卡| 看免费成人av毛片| 国产精品女同一区二区软件| 1024香蕉在线观看| 色哟哟·www| 亚洲精品国产av蜜桃| 精品国产一区二区三区久久久樱花| 国产精品一区二区在线观看99| 日韩精品免费视频一区二区三区| 欧美黄色片欧美黄色片| 如日韩欧美国产精品一区二区三区| av国产精品久久久久影院| 九九爱精品视频在线观看| 午夜福利在线观看免费完整高清在| 久久久精品94久久精品| 热re99久久精品国产66热6| 成人亚洲精品一区在线观看| 午夜激情久久久久久久| 秋霞在线观看毛片| 欧美日本中文国产一区发布| 桃花免费在线播放| 曰老女人黄片| 国产精品二区激情视频| av片东京热男人的天堂| 久久久久精品性色| 久久精品国产亚洲av涩爱| 国产国语露脸激情在线看| 男人爽女人下面视频在线观看| 亚洲av电影在线进入| 日本爱情动作片www.在线观看| 久久久欧美国产精品| 精品久久蜜臀av无| 免费少妇av软件| 国产国语露脸激情在线看| 婷婷成人精品国产| 亚洲天堂av无毛| 少妇熟女欧美另类| 免费高清在线观看视频在线观看| 日韩一区二区三区影片| 汤姆久久久久久久影院中文字幕| 国产乱来视频区| 日韩一本色道免费dvd| 在线 av 中文字幕| 最近最新中文字幕大全免费视频 | 一本—道久久a久久精品蜜桃钙片| 熟妇人妻不卡中文字幕| 精品国产乱码久久久久久小说| 国产精品一二三区在线看| 日韩av不卡免费在线播放| a级毛片黄视频| 夫妻性生交免费视频一级片| 亚洲国产最新在线播放| 国产av精品麻豆| 国产免费视频播放在线视频| 波野结衣二区三区在线| 日韩欧美精品免费久久| 免费观看av网站的网址| www日本在线高清视频| 在线天堂最新版资源| 精品卡一卡二卡四卡免费| 看免费成人av毛片| 亚洲 欧美一区二区三区| 中文字幕人妻丝袜一区二区 | 少妇 在线观看| 永久免费av网站大全| 最近手机中文字幕大全| 亚洲欧美一区二区三区黑人 | 最新的欧美精品一区二区| 国产片特级美女逼逼视频| 国产免费福利视频在线观看| www日本在线高清视频| 成年美女黄网站色视频大全免费| 亚洲综合精品二区| 亚洲久久久国产精品| 精品久久蜜臀av无| 一区在线观看完整版| 精品一区二区免费观看| 91精品国产国语对白视频| 久久久久久久国产电影| 叶爱在线成人免费视频播放| 欧美中文综合在线视频| 亚洲国产精品成人久久小说| 免费在线观看完整版高清| 国产黄色免费在线视频| 国产高清国产精品国产三级| 免费看av在线观看网站| 国产精品.久久久| 精品午夜福利在线看| 亚洲精品成人av观看孕妇| 岛国毛片在线播放| 桃花免费在线播放| 国产日韩欧美视频二区| 99精国产麻豆久久婷婷| √禁漫天堂资源中文www| 精品久久久精品久久久| 精品亚洲乱码少妇综合久久| 日韩电影二区| 日韩av在线免费看完整版不卡| 色哟哟·www| 国产精品一二三区在线看| 中文字幕人妻丝袜制服| 国产亚洲午夜精品一区二区久久| 久久久久国产精品人妻一区二区| 亚洲,欧美,日韩| 精品久久久久久电影网| 成年人午夜在线观看视频| 久久鲁丝午夜福利片| 欧美中文综合在线视频| 亚洲国产成人一精品久久久| 中文天堂在线官网| 卡戴珊不雅视频在线播放| 一区二区日韩欧美中文字幕| 又大又黄又爽视频免费| 国产亚洲一区二区精品| 大话2 男鬼变身卡| 啦啦啦视频在线资源免费观看| 精品视频人人做人人爽| 精品国产乱码久久久久久小说| 亚洲婷婷狠狠爱综合网| 激情视频va一区二区三区| 亚洲av电影在线观看一区二区三区| 久久人人爽人人片av| 人妻一区二区av| 亚洲综合色惰| 人妻系列 视频| av有码第一页| 大片电影免费在线观看免费| 性色av一级| 天天操日日干夜夜撸| 亚洲国产欧美网| 国产精品不卡视频一区二区| 色婷婷av一区二区三区视频| 黄网站色视频无遮挡免费观看| 欧美激情极品国产一区二区三区| 国产精品偷伦视频观看了| 黑人欧美特级aaaaaa片| 高清在线视频一区二区三区| 欧美精品人与动牲交sv欧美| 日本av手机在线免费观看| 国产日韩欧美在线精品| 国产黄频视频在线观看| 亚洲四区av| 在线天堂中文资源库| 日本猛色少妇xxxxx猛交久久| 亚洲精品中文字幕在线视频| 在线天堂中文资源库| 国产极品粉嫩免费观看在线| 欧美国产精品va在线观看不卡| 尾随美女入室| 国产男女内射视频| 777久久人妻少妇嫩草av网站| 黄色一级大片看看| 欧美另类一区| 一个人免费看片子| 日本色播在线视频| 一二三四在线观看免费中文在| 国产白丝娇喘喷水9色精品| 丰满迷人的少妇在线观看| 日韩电影二区| av免费在线看不卡| 国产爽快片一区二区三区| 亚洲av国产av综合av卡| 美国免费a级毛片| 日本黄色日本黄色录像| 国产精品无大码| 国产免费视频播放在线视频| 日本欧美视频一区| 亚洲美女黄色视频免费看| 欧美人与性动交α欧美软件| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成人手机av| 超碰成人久久| 啦啦啦在线观看免费高清www| 国产爽快片一区二区三区| 在线观看免费日韩欧美大片| 美女福利国产在线| 宅男免费午夜| 国产成人91sexporn| 777久久人妻少妇嫩草av网站| 免费高清在线观看日韩| 久久久久人妻精品一区果冻| 国产精品香港三级国产av潘金莲 | 麻豆精品久久久久久蜜桃| 亚洲欧美一区二区三区黑人 | 午夜福利在线免费观看网站| 国产成人av激情在线播放| 久久毛片免费看一区二区三区| 国产精品免费大片| 午夜福利视频在线观看免费| 99热网站在线观看| 午夜精品国产一区二区电影| 国产熟女午夜一区二区三区| 国产精品久久久久久精品电影小说| 久久女婷五月综合色啪小说| 青春草亚洲视频在线观看| 久久人人97超碰香蕉20202| av.在线天堂| 久久ye,这里只有精品| 一区福利在线观看| 日韩制服骚丝袜av| 侵犯人妻中文字幕一二三四区| 美女脱内裤让男人舔精品视频| 久久精品久久精品一区二区三区| 亚洲国产精品999| 久久久久人妻精品一区果冻| av视频免费观看在线观看| 大片电影免费在线观看免费| 在线观看一区二区三区激情| 色婷婷av一区二区三区视频| 伊人久久大香线蕉亚洲五| 久久99精品国语久久久| 亚洲精品国产av成人精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av网站免费在线观看视频| 国产成人av激情在线播放| 久久久国产欧美日韩av| 91国产中文字幕| 国产日韩欧美视频二区| 母亲3免费完整高清在线观看 | 免费少妇av软件| 18禁裸乳无遮挡动漫免费视频| av在线app专区| av片东京热男人的天堂| 国产乱人偷精品视频| 日本vs欧美在线观看视频| 十八禁高潮呻吟视频| 性少妇av在线| 亚洲av在线观看美女高潮| 欧美bdsm另类| 欧美成人午夜精品| 国产有黄有色有爽视频| 黄片播放在线免费| 十八禁高潮呻吟视频| 大香蕉久久成人网| 久久精品夜色国产| 久久免费观看电影| 亚洲欧美清纯卡通| 国产免费现黄频在线看| 久久97久久精品| 亚洲国产日韩一区二区| 黄色怎么调成土黄色| 高清视频免费观看一区二区| 9色porny在线观看| 亚洲精品美女久久久久99蜜臀 | 美女脱内裤让男人舔精品视频| 9191精品国产免费久久| 成人免费观看视频高清| 欧美日韩视频高清一区二区三区二| 亚洲国产精品国产精品| 超色免费av| 最近2019中文字幕mv第一页| 国产在线视频一区二区| 观看美女的网站| 国产伦理片在线播放av一区| 欧美另类一区| 久久精品aⅴ一区二区三区四区 | av女优亚洲男人天堂| 麻豆av在线久日| 少妇 在线观看| 欧美亚洲 丝袜 人妻 在线| 国产免费视频播放在线视频| 亚洲欧美日韩另类电影网站| 亚洲精品美女久久久久99蜜臀 | 国产成人精品一,二区| 99香蕉大伊视频| 女的被弄到高潮叫床怎么办| a级毛片在线看网站| a 毛片基地| 久久久久久久久久久久大奶| 少妇熟女欧美另类| 中文欧美无线码| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧洲国产日韩| 午夜福利视频在线观看免费| 一本—道久久a久久精品蜜桃钙片| 亚洲国产精品国产精品| 精品久久久久久电影网| 国产精品久久久久久久久免| 最新的欧美精品一区二区| 久久国产亚洲av麻豆专区| 日本黄色日本黄色录像| 黄色毛片三级朝国网站| 国产精品一区二区在线观看99| 国语对白做爰xxxⅹ性视频网站| 国产av精品麻豆| 99久久综合免费| 精品国产一区二区三区久久久樱花| 精品国产国语对白av| 精品国产乱码久久久久久男人| 五月天丁香电影| 午夜日本视频在线| 啦啦啦啦在线视频资源| 波多野结衣一区麻豆| 婷婷成人精品国产| 波多野结衣一区麻豆| 一区二区三区激情视频| 1024视频免费在线观看| 亚洲欧美精品综合一区二区三区 | 日本色播在线视频| 国产精品久久久久久久久免| 欧美日韩精品成人综合77777| 国产无遮挡羞羞视频在线观看| 免费在线观看黄色视频的| 亚洲美女搞黄在线观看| 亚洲综合色惰| 两个人免费观看高清视频| 少妇猛男粗大的猛烈进出视频| 精品国产乱码久久久久久男人| 高清欧美精品videossex| 纵有疾风起免费观看全集完整版| 一级黄片播放器| 国产一区二区三区av在线| 日韩一区二区三区影片| 国产黄色免费在线视频| 国产亚洲一区二区精品| 亚洲美女黄色视频免费看| 亚洲国产欧美日韩在线播放| 成人国产麻豆网| 欧美精品高潮呻吟av久久| 亚洲av福利一区| 成年女人在线观看亚洲视频| av一本久久久久| 捣出白浆h1v1| 亚洲精品国产av成人精品| 麻豆av在线久日| 亚洲精品国产av蜜桃| 在线观看国产h片| 亚洲一码二码三码区别大吗| 久久精品国产a三级三级三级| 人人澡人人妻人| 老鸭窝网址在线观看| 精品亚洲成a人片在线观看| 丰满乱子伦码专区| 看免费成人av毛片| 欧美人与性动交α欧美软件| 母亲3免费完整高清在线观看 |