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

    MODFLOW在新疆灌區(qū)地下水資源開發(fā)評價中的應用

    2011-02-23 07:05:10董新光
    地下水 2011年4期
    關鍵詞:開采量水源地含水層

    劉 磊,董新光,吳 斌

    (1.新疆農業(yè)大學 水利與土木工程學院,新疆 烏魯木齊830052;2.新疆水利廳,新疆 烏魯木齊830000)

    新和縣位于新疆維吾爾自治區(qū)天山南麓,塔里木盆地北緣,東隔渭干河與庫車相望,北隔確勒塔格山與拜城縣相鄰,南連沙雅縣,東西長136 km,南北寬91 km,全縣總面積為8 823 km2,耕地面積49.35×m 104萬畝。地處歐亞大陸腹地,塔里木盆地北部,遠離海洋,北部和西部受天山屏障的阻隔,大西洋和印度洋暖濕氣流部分可翻越帕米爾高原或天山進入本縣,同時東南面又受塔克拉瑪干大沙漠的干熱氣流影響,從而使新和縣呈現(xiàn)暖溫帶大陸性干旱氣候的特征。

    新疆的農業(yè)屬于灌溉型農業(yè),尤其是在南疆地表水缺乏的情況下,有限的且不斷減少的地表水資源已無法保證農業(yè)經濟可持續(xù)發(fā)展。為改善新和縣灌區(qū)耕地農業(yè)用水季節(jié)性和高峰期水資源供需矛盾,提高灌溉保證程度,保證農業(yè)經濟可持續(xù)發(fā)展,科學合理的開發(fā)利用該區(qū)域地下水資源是一條有效途徑。

    1 研究區(qū)水文地質概況

    評價區(qū)處于渭干河沖洪積傾斜平原中上部,渭干河是該沖洪積平原區(qū)第四紀松散層孔隙地下水形成的唯一來源,渭干河在出山口后形成了巨大的沖洪積扇。渭干河沖洪積扇的中上部主要為地下水的徑流區(qū),地下水徑流條件由北向南,由扇頂至扇前緣逐漸變差,隨著地層巖性的逐漸變細,地形坡度變緩,潛水水位埋深也逐漸變淺。在沖洪積扇的下部,由于地形坡度變緩,地層顆粒變細,從而減緩了地下水向下游排泄的速度,垂向排泄逐漸成為主要的排泄方式。

    評價區(qū)地下水的補給方式主要有上游地下水側向流入補給和區(qū)內地表水的垂向入滲補給,大氣降水入滲補給鑒于區(qū)內降雨量小且次降雨量很少達到10 mm以上,在此可予以忽略。評價區(qū)地層結構簡單,地層顆粒較粗大,滲透性好,地下水在接受了渭干河及上游渠系水入滲補給后,以側向徑流流入的方式補給本區(qū)。

    評價區(qū)地處渭干河西岸,渠系發(fā)育,尤魯都斯干渠、塔什艾日克干渠、依其艾日克干渠和沙雅總干渠等多條干渠均在工作區(qū)分布或穿過,同時,區(qū)內支、斗渠配套齊全,形成縱橫交錯的水系網,除干渠做防滲外,支、斗渠多防滲效果差或未防滲,從而使渠道滲漏成為測區(qū)地下水的另一主要補給源。此外,做為渭干河流域的古老灌區(qū),評價區(qū)內農田分布廣泛,耕作層主要為粉土和粉砂,田間灌溉水大量入滲,從而形成了本區(qū)地下水的另一重要補給源。

    項目區(qū)地層巖性在314國道以北主要以砂礫石夾中粗砂、粉土、粉質粘土為主,以南逐漸變?yōu)橹写稚皧A粉土、粉質粘土層,第四系松散堆積物在工作區(qū)厚度大于150 m。巨厚的松散堆積層為地下水提供了良好的賦存空間,地下水類型主要為潛水,局部存在微承壓水。

    在項目區(qū)東北部,含水層結構為單一的砂礫石潛水賦存區(qū),水位埋深多為3~5 m,在縣城北部和東北部大于5 m,在西南部多在1~3 m,局部地段由于受地形影響在3~5 m之間。地下水流向,總體上與地形坡降一致,即向南、西南方向徑流,但局部因受引水渠及排水渠的影響而稍有改變。

    地下水的徑流條件主要受地形條件和地層巖性控制,測區(qū)地形平坦開闊,地勢北高南低,地形坡降1.0‰ ~2.5‰。含水層巖性主要為礫石、中粗砂、細砂和粉砂,自北向南,顆粒逐漸變細,徑流條件也相應逐漸變差。評價區(qū)地下水的排泄方式有潛水蒸發(fā)蒸騰排泄、地下水側向排泄、排渠排泄和人工開采等。

    2 模型的建立

    2.1 模型研究范圍

    地下水模型的研究范圍與地下水調查范圍一致,在平面上呈半扇形沿東北-西南方向展布,總面積為295 km2,東北端角點坐標為(640000,4620000),東南端角點坐標為(635500,4590000),西北角點的坐標為(621700,4600000)(見圖1)。

    圖1 地下水模型研究范圍

    2.2 含水層的概化

    根據模型范圍內含水層介質的成因、物質組成、空間分布規(guī)律,將模擬對象概化為非均質各項同性含水層。平面上根據含水層分布規(guī)律、結合勘探孔抽水試驗資料進行了參數分區(qū);在垂直方向大部分地段含水層變化不明顯,僅作為一層,模擬深度為150 m。

    2.3 數學模型

    根據含水層埋藏分布、巖性及水流特征,模型區(qū)數學模型選用一層二維非穩(wěn)定流數值模型,其具體數學模型表達式如下。

    式中,Kxx、Kyy分別為 X、Y方向的滲透系數,L/T;H 為水頭值,L;M為含水層厚度、ε為源匯項,L/T;S為給水度或彈性釋水系數,潛水區(qū)取重力給水度,承壓區(qū)取彈性釋水系數;Ω為模擬范圍;n為邊界面的外法線方向;Γ為側邊界;B為底邊界。

    對上述地下水流數學模型,采用正交網格三維有限差分數值方法求解。為保證計算結果的收斂性,采用交替方向隱式差分格式,利用強隱式(SIP)求解差分方程。采用美國地調局(USGS)三維地下水流模擬軟件包Processing MODFLOW進行數值計算。

    研究區(qū)以2000年8月實測地下水位作為模擬預測模型的初始水位。利用MODFLOW自帶的工具Digitizer對實測地下水等水位線圖數字化,再用插值模塊Field Interpolator進行插值,對每一個單元格賦初始水頭值。

    2.4 邊界條件的處理

    依據模擬區(qū)水文地質條件,平面上邊界條件概化如下(見圖2),將底邊界處理為隔水邊界,上邊界作為開放邊界,考慮渠系入滲和人工開采,分別用Processing Modflow中的Recharge(面狀補給)、Well(水井)程序包來處理;側向補給邊界作為定流量邊界,用Well程序包來處理;側向排泄邊界用General Head Boundary(一般水頭邊界)程序包來處理;排水渠道用 Drain(排水)程序包來處理;蒸發(fā)用Evapotranspiration(蒸發(fā))程序包來處理。其中AB邊界為第二類補給流量邊界;BC和FA邊界為零流量邊界;CD、DE、EF邊界為第二類排泄流量邊界。

    2.4.1 空間與時間離散

    空間上采用等間距有限差分的離散方法,在地下水模型中采用自動剖分,模擬范圍內將含水層離散為1層、81行、69列,差分網格的大小為400 m×400 m,網格平均單元面積0.16 km2。模型計算區(qū)單元數(有效單元數)為1 842個,面積為 294.72 km2。

    時間上以2000年1月作為初始時刻,根據地下水位的觀測時間,模擬時段以月為單位,將模型調試與參數率定過程的7年分成84個時段,將2000年的初始水位運行到2007年,再以模擬好的07年的流場作為規(guī)劃年的初始流場,規(guī)劃年模型調試與參數率定過程的10年分為120個時段,模擬計算以天為最小單位。

    圖2 模型邊界條件概化及規(guī)劃開采井平面分布圖

    2.4.2 模型中數據選用

    地面高程數據從1:50000比例尺的地形圖上提取,然后利用MODFLOW自帶的插值模塊Field Interpolator插值,再賦值到各個單元格,作為地面高程。底部高程值用地面高程減去150 m。結合含水層成因、空間分布規(guī)律及勘探孔抽水試驗資料進行了參數分區(qū)。依據抽水試驗資料確定含水層的滲透系數和給水度。

    面狀補給量用Recharge程序包進行模擬,紅旗農場的面狀補給量為3 373×104m3/a,試驗林場的面狀補給量為1 685×104m3/a。側向補給量及側向排泄量用well程序包進行模擬,灌區(qū)邊界BC、CD和DE的側向補給量分別為316×104m3/a、773 ×104m3/a和621 ×104m3/a,灌區(qū)邊界 FI、IJ和JK的側向排泄量分別為、3 405×104m3/a、2 790×104m3/a和945×104m3/a。

    地下水開采量用well程序包進行模擬。規(guī)劃機井開采期每月的地下水開采量依據灌區(qū)機井開采方案中的數據進行賦值。研究區(qū)排水渠道的排水量用Drain程序包進行模擬,排水量為 2419×104m3/a。研究區(qū)蒸發(fā)量用Evapotranspiration程序包進行模擬,蒸發(fā)量為689×104m3/a。

    3 模型參數率定

    3.1 率定依據和方法

    2007年8月實測等水位線圖作為模型檢驗的流場圖,2000年的等水位線經模型運行7年后,與2007年的等水位線擬合較好的等水位線作為規(guī)劃年的初始流場。2000年到2007年評價區(qū)地下水動態(tài)觀測資料作為率定含水層參數的依據。以地下水均衡計算的地下水補給量、排泄量作為模型檢驗的輸入,以控制性勘探孔地下水位作為模型率定的檢驗標準。

    采用自動與手動相結合的方法,通過計算水位和實際水位的擬合分析,不斷地修改反復調整參數,當兩者之間誤差“最小時”,即認為此時的參數值代表含水層的參數。

    3.2 平面流場的擬合

    經調試計算,2000年的等水位線經模型運行7年后,得到2007年8月流場與實測流場比較,趨勢基本一致,將2007年8月的流場作為規(guī)劃年的初始流場,模擬潛水初始流場與實測流場比較見圖3。計算值與實際值偏差最大為0.5 m,大部分小于0.1 m。地下水位擬合誤差為0.0~0.2 m的單元數占總單元數的82%,0.2~0.5 m單元數占總單元數的18%。故認為擬合效果較好。

    圖3 07年8月模擬流場與實際流場擬合圖

    3.3 數值模擬結果與地下水均衡結果對比分析

    通過數值法計算的地下水補給量和排泄量,與地下水均衡法計算的補給量和排泄量結果比較分析(見表1),可進一步提高地下水資源的評價精度。

    從全區(qū)計算結果分析:地下水均衡法計算的地下水總補給量為12 671×104m3/a,模型計算為12 511×104m3/a,相差160×104m3/a,相對誤差為1.26%。地下水均衡法計算的總排泄量為13 224×104m3/a,模型計算為13 113×104m3/a,相差111×104m3/a,相對誤差為0.84%。地下水均衡法計算的均衡差為-552×104m3/a,模型計算為-602×104m3/a,相差 49 ×104m3/a。

    地下水模型進行參數和各項補排量調試和調整的主要依據是將2000年流場運行7年后要與2007年實測的流場基本一致,以2007年的水均衡計算的各項地下水補排量作為模型調試的初始值。在率定參數同時,調整各項補排量,使地下水模擬的初始流場與2007年的實測流場擬合最佳為止,得出評價區(qū)模型模擬各項補給和排泄量與地下水均衡法算的量擬合最佳。本次運用兩種評價方法,從流場上講,研究區(qū)地下水模擬的初始流場與2007年的實測流場基本一致;從補排關系上,水均衡法與地下水模型計算的總補給量、總排泄量和均衡結果差異不大。故認為此計算成果是可信的。

    表1 地下水均衡法與數值法計算結果比較 104m3

    4 應用

    4.1 開采方案布置

    依據《新和縣布喀塔木水源地(擴建)工程報告》,水源地擬布機井50眼,計算區(qū)原水源地現(xiàn)狀開采量為6 786.14×104m3/a,原水源地如按設計開采量開采時計算區(qū)總開采量為(6 786.14~1 880) ×104m3/a=4 906.14 ×104m3/a,計算區(qū)可開采量為6 351.81×104m3/a,則尚可開(6 351.81~4 906.14) ×104m3/a=1 445.67×104m3/a。即計算區(qū)在今后的地下水開發(fā)利用過程中,尚可開采1 445.67×104m3/a的地下水量。

    設計開采量是在可開采量滿足的條件下綜合灌區(qū)缺水量、缺水時間及考慮本水源地置換地表水方案進行確定的。根據這50眼機井的不同的規(guī)劃開采量,擬定了3種方案對評價區(qū)在項目運行5年后和10年后的水量與水位進行預測。

    方案一:按照原水源地的開采井數和開采量繼續(xù)運行。

    方案二:根據《新和縣布喀塔木水源地(擴建)工程報告》,布喀塔木水源地預計新建機井50眼,總開采量為1 000×104m3/a。原有機井的開采量保持不變,仍然為6 786.14×104m3/a,總的開采量為 7 786.14 ×104m3/a。

    方案三:根據《新和縣布喀塔木水源地(擴建)工程報告》,布喀塔木水源地預計新建機井50眼,總開采量為1 000×104m3/a。原有機井的開采量減小到設計開采量4 906.14×104m3/a,總的開采量為 5 906.14 ×104m3/a。

    4.2 各方案模型預測

    通過不同開采方案預測,方案一與方案二的水井開采量大于水源地的可開采量,運行5年后(2012年)潛水區(qū)大部分地段地下水位在現(xiàn)狀年基礎上下降0.3~2 m,集中開采區(qū)下降1~2.5 m。從方案一到方案二,隨開采量增大,地下水位降幅增大,同等降幅面積在擴大。方案三開采量小于水源地的可開采量,原水源地水位在現(xiàn)狀年基礎上有所上升0.1~0.5 m,新建開采井在現(xiàn)狀年基礎上水位下降0.1~1.1 m。

    各方案運行10年后(2017年)方案一與方案二的水井潛水區(qū)大部分地段地下水位在現(xiàn)狀年基礎上下降0.4~3 m,集中開采區(qū)下降2~3 m。從方案2到方案1,隨開采量增大,地下水位降幅增大。方案三開采量小于水源地的可開采量,原水源地水位在現(xiàn)狀年基礎上有所上升0.2~1.0 m,新建開采井在現(xiàn)狀年基礎上水位下降0.2~1.2 m。

    5 結語

    對新和縣老灌區(qū)通過MODFLOW軟件建立的地下水數值模擬通過驗證符合研究區(qū)地下水含水層的邊界條件和實際真實情況,模型運行穩(wěn)定可靠。利用模型對新灌區(qū)的地下水的三個開采方案進行了預測分析評價,得出按照原方案或加大開采量都超過了改研究區(qū)的可開采量,運行5年和10年后該區(qū)地下水的水量和水位都有下降趨勢,故在增加開采井的情況下減少開采量是比較符合地下水資源合理利用的方案。模擬說明這一優(yōu)化方案可以較好的緩解該區(qū)域地下水水資源與農業(yè)發(fā)展的矛盾,為該區(qū)域的農業(yè)經濟發(fā)展和地下水開采提供了參考。

    [1]Mc Donald G,Michael and Arlen W,Harbaugh.A modular threedimensional finite-difference ground-water flow model[M].United State Government Printing Office,Washington,1998.

    [2]盧文喜.地下水運動數值模擬過程中邊界條件問題探討[J].水利學報,2003(3):33-36.

    [3]吳劍鋒,朱學愚.由MODFLIOW淺談地下水流數值模擬軟件的發(fā)展趨勢[J].工程勘察,2000,(2):12-15.

    [4]王智,高瑾,董新光.MODFLOW在三工河流域地下水資源評價中的應用[J].新疆農業(yè)大學報,2002,25(4):29 -34.

    [5]卞玉梅,盧文喜,馬洪云.Visual MODFLOW在水源地地下水數值模擬中的應用[J].東北水利水電,2006,(3).

    [6]Sun Minzheng,Chunmiao.Long-term groundwater management by a MODFLOW based dynamic optimization tool[J].Journal of the American Water Resources Association,1999,35(1):99-111.

    猜你喜歡
    開采量水源地含水層
    全球多個含水層里的水正快速流失
    再談河北省灤平縣馬營子鄉(xiāng)高鍶天然礦泉水特征與開采量估算
    礦產勘查(2020年3期)2020-12-28 00:27:22
    美國西部奧加拉拉含水層水位下降原因初探
    基于因果分析圖法的飲用水源地保護探討
    利用統(tǒng)計分析法預測地熱水可開采量應注意的問題
    河北地質(2016年1期)2016-03-20 13:51:59
    中國新疆石油開采量總額增長
    中亞信息(2016年1期)2016-03-07 01:00:42
    全球地下含水層下降驚人:要被抽干了
    地理教學(2015年14期)2015-03-31 20:04:53
    巖溶含水層水流模型研究進展
    南水北調中線水源地土壤侵蝕經濟損失估算
    吉林省供水水源地問題探討
    吉林地質(2014年4期)2014-03-11 16:47:59
    久久人人爽av亚洲精品天堂| 亚洲欧美清纯卡通| 我要看黄色一级片免费的| 日韩 亚洲 欧美在线| 欧美亚洲日本最大视频资源| 热99久久久久精品小说推荐| 久久99一区二区三区| 如日韩欧美国产精品一区二区三区| 国产欧美日韩综合在线一区二区| 国产一区二区 视频在线| 啦啦啦中文免费视频观看日本| 国产精品久久久久久精品电影小说| 一区二区三区激情视频| 在线观看人妻少妇| 啦啦啦 在线观看视频| 久久中文字幕一级| 各种免费的搞黄视频| 天堂8中文在线网| 欧美人与性动交α欧美软件| av线在线观看网站| 久久av网站| 99国产精品一区二区蜜桃av | 亚洲欧洲精品一区二区精品久久久| 一级毛片精品| 丰满饥渴人妻一区二区三| 无遮挡黄片免费观看| 99精品欧美一区二区三区四区| www.999成人在线观看| 国产亚洲精品久久久久5区| 制服诱惑二区| 国产精品一区二区在线不卡| 久久久久网色| 两个人免费观看高清视频| 午夜福利在线免费观看网站| 操出白浆在线播放| 黄色怎么调成土黄色| 老司机影院成人| 午夜视频精品福利| 久久亚洲国产成人精品v| 免费观看av网站的网址| 操美女的视频在线观看| 亚洲中文日韩欧美视频| 久久久久久久大尺度免费视频| 欧美黑人精品巨大| 亚洲七黄色美女视频| 欧美日韩福利视频一区二区| 亚洲综合色网址| 免费一级毛片在线播放高清视频 | 12—13女人毛片做爰片一| 在线 av 中文字幕| 成年人午夜在线观看视频| 一级黄色大片毛片| 精品福利观看| 老汉色∧v一级毛片| 久久久久久人人人人人| 男女午夜视频在线观看| 国产一卡二卡三卡精品| 捣出白浆h1v1| 亚洲欧美一区二区三区黑人| 久久久久国内视频| 久久九九热精品免费| 香蕉国产在线看| 日本一区二区免费在线视频| 中文精品一卡2卡3卡4更新| 免费人妻精品一区二区三区视频| 高清黄色对白视频在线免费看| 免费高清在线观看日韩| 亚洲色图 男人天堂 中文字幕| 日本猛色少妇xxxxx猛交久久| 国产免费一区二区三区四区乱码| 91字幕亚洲| 宅男免费午夜| 一级,二级,三级黄色视频| 国产精品秋霞免费鲁丝片| 久久人人爽av亚洲精品天堂| 欧美日韩精品网址| 在线观看人妻少妇| 久久人妻福利社区极品人妻图片| 色婷婷av一区二区三区视频| 91九色精品人成在线观看| 国产高清视频在线播放一区 | 狂野欧美激情性xxxx| 黄色a级毛片大全视频| 我的亚洲天堂| 亚洲午夜精品一区,二区,三区| 亚洲精品成人av观看孕妇| 国产成人系列免费观看| 叶爱在线成人免费视频播放| 国产一区二区在线观看av| 成人18禁高潮啪啪吃奶动态图| 欧美国产精品va在线观看不卡| 亚洲成人免费电影在线观看| 国产伦理片在线播放av一区| 十八禁人妻一区二区| 亚洲综合色网址| 亚洲专区字幕在线| 国产在线观看jvid| 满18在线观看网站| 美女扒开内裤让男人捅视频| 999久久久国产精品视频| 青春草视频在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| bbb黄色大片| 搡老乐熟女国产| 精品高清国产在线一区| 欧美日韩黄片免| 午夜影院在线不卡| 久久99一区二区三区| 国产精品亚洲av一区麻豆| 丝袜美足系列| 亚洲精品一区蜜桃| 亚洲人成77777在线视频| 欧美日韩视频精品一区| av国产精品久久久久影院| 欧美亚洲日本最大视频资源| 一二三四社区在线视频社区8| 中国美女看黄片| 国产亚洲欧美精品永久| 美女国产高潮福利片在线看| 精品视频人人做人人爽| 高清黄色对白视频在线免费看| 性少妇av在线| 国产亚洲精品一区二区www | 一级毛片精品| 韩国精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 91国产中文字幕| 成人影院久久| 欧美成人午夜精品| avwww免费| 性色av乱码一区二区三区2| 一进一出抽搐动态| 99九九在线精品视频| 在线亚洲精品国产二区图片欧美| 亚洲精品美女久久久久99蜜臀| 欧美亚洲 丝袜 人妻 在线| 国产精品 欧美亚洲| 亚洲三区欧美一区| 中文字幕人妻丝袜一区二区| 18禁黄网站禁片午夜丰满| 国产精品免费大片| 黄色视频在线播放观看不卡| 国产淫语在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产免费一区二区三区四区乱码| 他把我摸到了高潮在线观看 | 99香蕉大伊视频| 黑人欧美特级aaaaaa片| 国产精品一二三区在线看| 亚洲视频免费观看视频| 亚洲色图 男人天堂 中文字幕| 日本av手机在线免费观看| 欧美日韩黄片免| 亚洲av美国av| 亚洲欧美成人综合另类久久久| 久久青草综合色| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区三区在线臀色熟女 | 一级毛片精品| 亚洲成人手机| 亚洲av电影在线观看一区二区三区| 国产有黄有色有爽视频| 最近中文字幕2019免费版| 日本a在线网址| 亚洲一区二区三区欧美精品| 免费在线观看黄色视频的| 成人三级做爰电影| 国产99久久九九免费精品| 老司机影院成人| 激情视频va一区二区三区| 亚洲全国av大片| 午夜福利乱码中文字幕| 国产亚洲欧美在线一区二区| 欧美精品人与动牲交sv欧美| 欧美日韩视频精品一区| 美女福利国产在线| 无遮挡黄片免费观看| 亚洲国产精品一区二区三区在线| 成人18禁高潮啪啪吃奶动态图| 久久久久久免费高清国产稀缺| 99国产精品免费福利视频| 久久香蕉激情| 亚洲精品av麻豆狂野| 亚洲专区中文字幕在线| 久久人妻熟女aⅴ| 国产亚洲欧美精品永久| 蜜桃国产av成人99| 国产av精品麻豆| √禁漫天堂资源中文www| 九色亚洲精品在线播放| 十八禁网站网址无遮挡| 后天国语完整版免费观看| 亚洲欧美日韩高清在线视频 | 亚洲精品中文字幕在线视频| 国产一区二区三区在线臀色熟女 | 中文字幕最新亚洲高清| 少妇粗大呻吟视频| 国产一区二区激情短视频 | 欧美人与性动交α欧美软件| 久热这里只有精品99| 美女大奶头黄色视频| 久久久久久久精品精品| tube8黄色片| av超薄肉色丝袜交足视频| 成人影院久久| 国产深夜福利视频在线观看| 精品福利永久在线观看| 国精品久久久久久国模美| 久久精品国产a三级三级三级| 国产伦人伦偷精品视频| 国产高清国产精品国产三级| 每晚都被弄得嗷嗷叫到高潮| 最近中文字幕2019免费版| 日韩有码中文字幕| 一本大道久久a久久精品| 天天躁夜夜躁狠狠躁躁| 男女下面插进去视频免费观看| 国产亚洲午夜精品一区二区久久| 热re99久久精品国产66热6| 欧美激情极品国产一区二区三区| 国产成人欧美在线观看 | 亚洲中文字幕日韩| 18禁国产床啪视频网站| 国产熟女午夜一区二区三区| 少妇猛男粗大的猛烈进出视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧美一区二区三区黑人| 黑人欧美特级aaaaaa片| 999久久久精品免费观看国产| 黄色视频不卡| 精品久久蜜臀av无| 欧美成人午夜精品| 19禁男女啪啪无遮挡网站| 亚洲免费av在线视频| www日本在线高清视频| 国产淫语在线视频| 最黄视频免费看| 国产亚洲精品一区二区www | 亚洲欧洲精品一区二区精品久久久| 999久久久精品免费观看国产| 人人妻人人澡人人爽人人夜夜| 久久久久久亚洲精品国产蜜桃av| 97人妻天天添夜夜摸| 另类亚洲欧美激情| 日本一区二区免费在线视频| 高清欧美精品videossex| 两个人免费观看高清视频| 日本av免费视频播放| 午夜日韩欧美国产| 亚洲国产看品久久| 一区二区三区精品91| 91精品伊人久久大香线蕉| 无遮挡黄片免费观看| 久久精品成人免费网站| 久久99热这里只频精品6学生| www.自偷自拍.com| 国产成人精品无人区| 亚洲欧洲精品一区二区精品久久久| 少妇粗大呻吟视频| 高清黄色对白视频在线免费看| 久久热在线av| 黑人巨大精品欧美一区二区蜜桃| 色精品久久人妻99蜜桃| 国精品久久久久久国模美| 国产在线观看jvid| a在线观看视频网站| 丁香六月欧美| 国产真人三级小视频在线观看| 日本一区二区免费在线视频| 国产日韩欧美视频二区| 叶爱在线成人免费视频播放| 99国产精品99久久久久| 欧美老熟妇乱子伦牲交| 视频区图区小说| 国产老妇伦熟女老妇高清| 亚洲男人天堂网一区| 精品久久久久久电影网| 侵犯人妻中文字幕一二三四区| 19禁男女啪啪无遮挡网站| 国产男人的电影天堂91| 伊人久久大香线蕉亚洲五| 亚洲七黄色美女视频| 色视频在线一区二区三区| 无遮挡黄片免费观看| 桃红色精品国产亚洲av| 欧美 亚洲 国产 日韩一| 国产伦理片在线播放av一区| av又黄又爽大尺度在线免费看| 亚洲成人免费电影在线观看| 久久女婷五月综合色啪小说| 中国国产av一级| 丝袜喷水一区| 免费观看人在逋| 中文欧美无线码| 中文字幕人妻熟女乱码| 老司机深夜福利视频在线观看 | 狠狠狠狠99中文字幕| 青青草视频在线视频观看| 欧美日韩黄片免| 少妇裸体淫交视频免费看高清 | 国产欧美日韩一区二区精品| 一区二区三区乱码不卡18| 国产精品久久久人人做人人爽| 午夜老司机福利片| 秋霞在线观看毛片| 亚洲精品中文字幕一二三四区 | 亚洲人成电影观看| 窝窝影院91人妻| 精品国产一区二区三区四区第35| 欧美日韩一级在线毛片| 久热爱精品视频在线9| 国产男女超爽视频在线观看| 精品高清国产在线一区| 成人亚洲精品一区在线观看| 亚洲欧美日韩高清在线视频 | 黄频高清免费视频| 国产精品久久久久成人av| 国产免费现黄频在线看| 91精品国产国语对白视频| 99国产精品免费福利视频| 搡老乐熟女国产| 久热这里只有精品99| 婷婷丁香在线五月| 亚洲第一av免费看| 亚洲色图 男人天堂 中文字幕| 欧美 亚洲 国产 日韩一| 亚洲成人国产一区在线观看| 狂野欧美激情性xxxx| 精品国产超薄肉色丝袜足j| 91av网站免费观看| 久久久久久久精品精品| kizo精华| 久久久久久亚洲精品国产蜜桃av| 不卡一级毛片| 男女边摸边吃奶| 亚洲精品av麻豆狂野| 日韩 欧美 亚洲 中文字幕| 欧美日韩一级在线毛片| 国产视频一区二区在线看| 亚洲av美国av| 777米奇影视久久| 在线亚洲精品国产二区图片欧美| 这个男人来自地球电影免费观看| 久久精品亚洲熟妇少妇任你| 欧美激情 高清一区二区三区| 女人精品久久久久毛片| 不卡一级毛片| 亚洲一码二码三码区别大吗| 王馨瑶露胸无遮挡在线观看| 一二三四社区在线视频社区8| 美女扒开内裤让男人捅视频| 欧美+亚洲+日韩+国产| 最新的欧美精品一区二区| av电影中文网址| 久久亚洲精品不卡| 日本wwww免费看| 秋霞在线观看毛片| av网站免费在线观看视频| 国产亚洲精品一区二区www | 黑人操中国人逼视频| 黄片小视频在线播放| 老司机午夜福利在线观看视频 | 悠悠久久av| 欧美激情极品国产一区二区三区| 法律面前人人平等表现在哪些方面 | 热99国产精品久久久久久7| 十八禁网站网址无遮挡| 欧美少妇被猛烈插入视频| 新久久久久国产一级毛片| 国产精品成人在线| 黑人巨大精品欧美一区二区mp4| 91麻豆精品激情在线观看国产 | 久久 成人 亚洲| 精品视频人人做人人爽| 各种免费的搞黄视频| 在线亚洲精品国产二区图片欧美| 久久久久国内视频| 美女国产高潮福利片在线看| 人人妻人人澡人人爽人人夜夜| 9色porny在线观看| 日韩电影二区| 亚洲精品乱久久久久久| 国产深夜福利视频在线观看| 一区在线观看完整版| 欧美精品亚洲一区二区| 亚洲欧美一区二区三区黑人| 日韩免费高清中文字幕av| 欧美激情极品国产一区二区三区| 午夜免费鲁丝| 国产91精品成人一区二区三区 | 永久免费av网站大全| 免费在线观看日本一区| 久热这里只有精品99| av又黄又爽大尺度在线免费看| 欧美成人午夜精品| 99国产精品一区二区蜜桃av | 最近中文字幕2019免费版| 亚洲欧美日韩高清在线视频 | 19禁男女啪啪无遮挡网站| www.熟女人妻精品国产| 精品国产一区二区三区四区第35| 久久久久视频综合| 伊人久久大香线蕉亚洲五| 亚洲天堂av无毛| 国产精品久久久久久精品电影小说| 日韩欧美免费精品| 欧美日韩国产mv在线观看视频| 日韩,欧美,国产一区二区三区| 一级毛片电影观看| 巨乳人妻的诱惑在线观看| 精品久久久久久电影网| 乱人伦中国视频| 后天国语完整版免费观看| 免费在线观看视频国产中文字幕亚洲 | 满18在线观看网站| 欧美性长视频在线观看| 国产一区二区三区在线臀色熟女 | 一本久久精品| 午夜福利在线免费观看网站| 看免费av毛片| 亚洲,欧美精品.| 黑人猛操日本美女一级片| 97在线人人人人妻| 99久久精品国产亚洲精品| 超碰97精品在线观看| 中国美女看黄片| 国产精品99久久99久久久不卡| 日本a在线网址| av欧美777| 十分钟在线观看高清视频www| 久久精品国产a三级三级三级| 亚洲人成电影免费在线| 999精品在线视频| 久久性视频一级片| 丰满饥渴人妻一区二区三| 1024视频免费在线观看| 丝袜脚勾引网站| 91av网站免费观看| 亚洲欧美清纯卡通| 亚洲九九香蕉| 亚洲欧美一区二区三区久久| 国产精品自产拍在线观看55亚洲 | 国产免费现黄频在线看| 久久青草综合色| 99re6热这里在线精品视频| 黑人猛操日本美女一级片| av在线老鸭窝| 亚洲第一欧美日韩一区二区三区 | 精品乱码久久久久久99久播| 法律面前人人平等表现在哪些方面 | 90打野战视频偷拍视频| 日日爽夜夜爽网站| 国产精品久久久人人做人人爽| 国产亚洲欧美精品永久| 在线观看免费视频网站a站| 黄色片一级片一级黄色片| 超碰97精品在线观看| 精品人妻熟女毛片av久久网站| 久久久久精品国产欧美久久久 | 99热网站在线观看| 一区福利在线观看| 人人澡人人妻人| 五月开心婷婷网| 亚洲精品国产av成人精品| 欧美精品人与动牲交sv欧美| 欧美日韩黄片免| 精品一区二区三区四区五区乱码| 十八禁高潮呻吟视频| 久久久国产成人免费| 丝袜在线中文字幕| 久久久久久亚洲精品国产蜜桃av| 欧美av亚洲av综合av国产av| 一本—道久久a久久精品蜜桃钙片| 新久久久久国产一级毛片| 熟女少妇亚洲综合色aaa.| 亚洲av美国av| 天天躁日日躁夜夜躁夜夜| 欧美在线一区亚洲| 老司机靠b影院| 亚洲av男天堂| 另类亚洲欧美激情| 91精品国产国语对白视频| 黄频高清免费视频| 亚洲中文字幕日韩| 精品久久久久久电影网| 一区福利在线观看| av国产精品久久久久影院| 久久人妻福利社区极品人妻图片| 亚洲精品第二区| 国产极品粉嫩免费观看在线| 久久久久精品人妻al黑| 午夜福利影视在线免费观看| 亚洲av国产av综合av卡| 18禁观看日本| 精品久久久久久电影网| 久久久精品免费免费高清| av国产精品久久久久影院| 精品国产乱码久久久久久男人| 亚洲专区国产一区二区| bbb黄色大片| 精品一区二区三卡| 精品少妇久久久久久888优播| 黄色视频,在线免费观看| 亚洲午夜精品一区,二区,三区| 香蕉国产在线看| 精品亚洲成a人片在线观看| 亚洲欧洲日产国产| 国产在线视频一区二区| 欧美国产精品va在线观看不卡| 国产黄频视频在线观看| 国产一级毛片在线| 亚洲第一欧美日韩一区二区三区 | 日韩三级视频一区二区三区| 99国产精品一区二区三区| 免费看十八禁软件| 乱人伦中国视频| 欧美另类一区| 国产三级黄色录像| 一本大道久久a久久精品| 超色免费av| 精品人妻1区二区| av一本久久久久| 国产免费现黄频在线看| 国产精品1区2区在线观看. | 99re6热这里在线精品视频| 人妻人人澡人人爽人人| 免费人妻精品一区二区三区视频| videos熟女内射| 成人影院久久| 免费高清在线观看日韩| 亚洲精品成人av观看孕妇| 热99re8久久精品国产| 日韩 欧美 亚洲 中文字幕| 精品少妇久久久久久888优播| 黄色视频,在线免费观看| 狂野欧美激情性bbbbbb| 亚洲欧美成人综合另类久久久| 久久亚洲精品不卡| 国产1区2区3区精品| 51午夜福利影视在线观看| 动漫黄色视频在线观看| 日韩中文字幕视频在线看片| 国产欧美日韩精品亚洲av| 黑人操中国人逼视频| 国产精品 国内视频| 自线自在国产av| 中文精品一卡2卡3卡4更新| 亚洲av美国av| 午夜福利影视在线免费观看| 中文欧美无线码| 国产精品麻豆人妻色哟哟久久| 国产片内射在线| 精品久久久久久久毛片微露脸 | 少妇裸体淫交视频免费看高清 | 国产伦理片在线播放av一区| 精品国产国语对白av| 最近最新免费中文字幕在线| 亚洲国产成人一精品久久久| 久久久国产一区二区| 国产一区二区激情短视频 | 国产黄频视频在线观看| 国产野战对白在线观看| 黄片播放在线免费| 精品国产乱码久久久久久男人| 国产精品国产三级国产专区5o| 欧美黄色淫秽网站| 国产不卡av网站在线观看| 午夜福利视频精品| 亚洲黑人精品在线| 18禁观看日本| √禁漫天堂资源中文www| 亚洲中文av在线| 欧美精品亚洲一区二区| 日本撒尿小便嘘嘘汇集6| 欧美老熟妇乱子伦牲交| 天天躁日日躁夜夜躁夜夜| 国产色视频综合| 久久国产精品男人的天堂亚洲| 色视频在线一区二区三区| 亚洲精品中文字幕一二三四区 | 捣出白浆h1v1| 欧美性长视频在线观看| 国产熟女午夜一区二区三区| 日韩视频在线欧美| 精品国产乱子伦一区二区三区 | 国产精品1区2区在线观看. | 国产精品自产拍在线观看55亚洲 | 久久 成人 亚洲| av超薄肉色丝袜交足视频| 97在线人人人人妻| 亚洲欧美日韩高清在线视频 | av又黄又爽大尺度在线免费看| 亚洲国产精品999| 中文字幕精品免费在线观看视频| 久久香蕉激情| 午夜免费观看性视频| 欧美乱码精品一区二区三区| 超色免费av| 美女午夜性视频免费| 精品高清国产在线一区| 18在线观看网站| 天天影视国产精品| 黄色怎么调成土黄色| 交换朋友夫妻互换小说| 精品亚洲乱码少妇综合久久| 精品一区在线观看国产| 美女脱内裤让男人舔精品视频| 丝袜在线中文字幕| 久久久精品94久久精品| 国内毛片毛片毛片毛片毛片| 亚洲精品国产区一区二| 国产精品欧美亚洲77777| 精品国产乱码久久久久久男人| 久久久久久久大尺度免费视频| 国产高清videossex|