• <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视频| 久久人人爽人人爽人人片va| 婷婷色综合www| 成人二区视频| 日日摸夜夜添夜夜添av毛片| 免费高清在线观看视频在线观看| 久久久久视频综合| 2021少妇久久久久久久久久久| 欧美日韩亚洲高清精品| 国产成人精品福利久久| 欧美日韩精品成人综合77777| 欧美少妇被猛烈插入视频| 国产片特级美女逼逼视频| 在线观看人妻少妇| 中文字幕人妻丝袜制服| 亚洲精品乱久久久久久| 一本大道久久a久久精品| 亚洲国产成人一精品久久久| 少妇高潮的动态图| 草草在线视频免费看| 久久青草综合色| 建设人人有责人人尽责人人享有的| 一级毛片aaaaaa免费看小| 免费观看av网站的网址| 在线观看三级黄色| 亚洲欧美中文字幕日韩二区| 日韩欧美一区视频在线观看 | 噜噜噜噜噜久久久久久91| 精品一区二区三区视频在线| 精品一品国产午夜福利视频| 大陆偷拍与自拍| 少妇丰满av| 亚洲国产精品一区三区| 男人狂女人下面高潮的视频| 黄色配什么色好看| 国产精品福利在线免费观看| 嫩草影院新地址| 亚洲精品成人av观看孕妇| 在线观看免费视频网站a站| 天堂中文最新版在线下载| 日日撸夜夜添| 国产老妇伦熟女老妇高清| tube8黄色片| 国产黄片美女视频| 日韩伦理黄色片| 搡老乐熟女国产| 最近的中文字幕免费完整| 秋霞伦理黄片| av在线观看视频网站免费| av一本久久久久| 精品人妻熟女av久视频| 久久狼人影院| 欧美亚洲 丝袜 人妻 在线| 丰满乱子伦码专区| 欧美丝袜亚洲另类| 少妇高潮的动态图| 久久久久视频综合| 亚洲综合精品二区| 日韩精品免费视频一区二区三区 | 久久精品国产亚洲av天美| 五月开心婷婷网| 3wmmmm亚洲av在线观看| av卡一久久| 免费黄网站久久成人精品| 黄色一级大片看看| 极品人妻少妇av视频| 最后的刺客免费高清国语| 国产在线一区二区三区精| 三级经典国产精品| av免费在线看不卡| 99久久精品国产国产毛片| 日本欧美视频一区| 国产女主播在线喷水免费视频网站| 亚洲国产欧美在线一区| 亚洲精品成人av观看孕妇| 久久久欧美国产精品| 亚洲精品日韩在线中文字幕| 欧美性感艳星| 成年av动漫网址| 三级经典国产精品| 亚洲综合色惰| 高清欧美精品videossex| 少妇精品久久久久久久| 精品人妻偷拍中文字幕| 在线观看美女被高潮喷水网站| 久久精品久久久久久久性| 亚洲,欧美,日韩| 少妇人妻一区二区三区视频| 中文字幕av电影在线播放| 亚洲图色成人| 国产av码专区亚洲av| 97超视频在线观看视频| 51国产日韩欧美| 国产午夜精品一二区理论片| 美女主播在线视频| 亚洲综合精品二区| 亚洲欧美成人精品一区二区| 又大又黄又爽视频免费| 国内精品宾馆在线| 欧美老熟妇乱子伦牲交| 亚洲精品,欧美精品| 久久精品久久精品一区二区三区| 国产日韩欧美亚洲二区| 日日啪夜夜撸| av在线播放精品| 国产爽快片一区二区三区| 国产黄片美女视频| 亚洲欧洲日产国产| 三级国产精品片| 亚洲,欧美,日韩| 久久久久精品久久久久真实原创| 国模一区二区三区四区视频| 免费黄频网站在线观看国产| 亚洲精品乱码久久久久久按摩| 精品国产乱码久久久久久小说| 一级毛片久久久久久久久女| 亚洲av综合色区一区| 国产免费又黄又爽又色| 国产精品无大码| 校园人妻丝袜中文字幕| 久久久久久久国产电影| 免费观看在线日韩| 中国三级夫妇交换| 国产精品嫩草影院av在线观看| 午夜免费观看性视频| 最新的欧美精品一区二区| 久久综合国产亚洲精品| 久久久久久人妻| 国产av精品麻豆| 2018国产大陆天天弄谢| 嘟嘟电影网在线观看| 国产午夜精品一二区理论片| 美女主播在线视频| 色婷婷久久久亚洲欧美| 欧美变态另类bdsm刘玥| 成人国产麻豆网| 亚洲va在线va天堂va国产| 成人综合一区亚洲| 免费观看a级毛片全部| 国产精品久久久久久av不卡| 少妇的逼水好多| 午夜av观看不卡| 91精品伊人久久大香线蕉| 中文字幕人妻熟人妻熟丝袜美| 成人毛片60女人毛片免费| 国产日韩一区二区三区精品不卡 | 99热这里只有是精品50| 一区在线观看完整版| 热re99久久国产66热| av在线老鸭窝| 日本av免费视频播放| 在线观看免费高清a一片| 丁香六月天网| 蜜臀久久99精品久久宅男| 国产高清有码在线观看视频| 国产日韩欧美亚洲二区| 国产永久视频网站| 国产免费一区二区三区四区乱码| 日韩欧美一区视频在线观看 | 老女人水多毛片| 国产女主播在线喷水免费视频网站| 美女xxoo啪啪120秒动态图| 国产爽快片一区二区三区| 国产免费福利视频在线观看| 久久 成人 亚洲| 少妇精品久久久久久久| 精品久久久久久久久亚洲| 国语对白做爰xxxⅹ性视频网站| 大又大粗又爽又黄少妇毛片口| 亚洲精品乱码久久久久久按摩| 亚洲欧美日韩卡通动漫| 国产伦精品一区二区三区四那| av免费观看日本| 美女cb高潮喷水在线观看| av视频免费观看在线观看| 国产 一区精品| 能在线免费看毛片的网站| 日韩精品免费视频一区二区三区 | 欧美激情极品国产一区二区三区 | 精品国产一区二区久久| 男女边摸边吃奶| 一区二区三区四区激情视频| 日韩一区二区视频免费看| 欧美日韩在线观看h| 插阴视频在线观看视频| 国产精品99久久久久久久久| 亚洲自偷自拍三级| 欧美成人午夜免费资源| 久久久久久久久久久久大奶| 亚洲高清免费不卡视频| 免费看不卡的av| 最新中文字幕久久久久| 亚洲精品第二区| 看免费成人av毛片| 欧美少妇被猛烈插入视频| 国产91av在线免费观看| 国产黄色免费在线视频| 日韩免费高清中文字幕av| 另类亚洲欧美激情| 视频区图区小说| 亚洲精品日本国产第一区| 精品人妻熟女毛片av久久网站| 亚洲av男天堂| 美女视频免费永久观看网站| 美女大奶头黄色视频| 国产高清有码在线观看视频| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 乱系列少妇在线播放| 性高湖久久久久久久久免费观看| 成人影院久久| 丝袜喷水一区| 看非洲黑人一级黄片| 国产片特级美女逼逼视频| 人妻制服诱惑在线中文字幕| 成人特级av手机在线观看| 亚洲不卡免费看| 在线观看av片永久免费下载| 午夜91福利影院| 国产精品国产av在线观看| 少妇人妻 视频| 久久久久精品性色| 夜夜看夜夜爽夜夜摸| 视频中文字幕在线观看| 久热这里只有精品99| 最新中文字幕久久久久| 免费看av在线观看网站| a级毛片在线看网站| 亚洲婷婷狠狠爱综合网| 一个人看视频在线观看www免费| 爱豆传媒免费全集在线观看| 精品99又大又爽又粗少妇毛片| 纵有疾风起免费观看全集完整版| videossex国产| 亚洲性久久影院| 自拍欧美九色日韩亚洲蝌蚪91 | 最近最新中文字幕免费大全7| 麻豆乱淫一区二区| 欧美精品人与动牲交sv欧美| 亚洲欧美清纯卡通| 国产毛片在线视频| 亚洲精品国产成人久久av| 久久精品国产鲁丝片午夜精品| 亚洲精品视频女| 成人无遮挡网站| 99re6热这里在线精品视频| 一边亲一边摸免费视频| 中文在线观看免费www的网站| 亚洲真实伦在线观看| 99久国产av精品国产电影| 免费看光身美女| 夫妻午夜视频| 欧美精品高潮呻吟av久久| 亚洲熟女精品中文字幕| 国产精品一二三区在线看| 黑丝袜美女国产一区| 中国国产av一级| 丝袜喷水一区| 亚洲av男天堂| 亚洲精品一区蜜桃| 国内揄拍国产精品人妻在线| 热re99久久精品国产66热6| 欧美三级亚洲精品| 好男人视频免费观看在线| 麻豆成人av视频| 国产亚洲91精品色在线| 精品亚洲乱码少妇综合久久| 高清黄色对白视频在线免费看 | 亚洲中文av在线| 九九久久精品国产亚洲av麻豆| 一级毛片我不卡| 欧美人与善性xxx| 高清av免费在线| 精品人妻一区二区三区麻豆| 免费观看a级毛片全部| 午夜激情福利司机影院| 国产永久视频网站| 下体分泌物呈黄色| 免费观看在线日韩| 亚洲精品国产av成人精品| 黄色配什么色好看| 亚洲精品久久午夜乱码| 亚洲,欧美,日韩| 免费av中文字幕在线| av在线app专区| 亚洲美女黄色视频免费看| 精品人妻熟女av久视频| 中文天堂在线官网| 国产精品熟女久久久久浪| 日韩熟女老妇一区二区性免费视频| 欧美性感艳星| 日日摸夜夜添夜夜添av毛片| 国产高清国产精品国产三级| 久久这里有精品视频免费| 免费不卡的大黄色大毛片视频在线观看| 国产一级毛片在线| av播播在线观看一区| 狂野欧美激情性bbbbbb| 亚洲精品aⅴ在线观看| 一本久久精品| 国产有黄有色有爽视频| 国产免费一区二区三区四区乱码| 亚洲成色77777| 国产综合精华液| 99久久精品一区二区三区| 色吧在线观看| 精品人妻一区二区三区麻豆| 亚洲av男天堂| videos熟女内射| 一区二区av电影网| 黑人巨大精品欧美一区二区蜜桃 | 久久热精品热| 水蜜桃什么品种好| 91午夜精品亚洲一区二区三区| 高清欧美精品videossex| 欧美最新免费一区二区三区| 久久韩国三级中文字幕| 岛国毛片在线播放| 一区二区三区四区激情视频| 黄色配什么色好看| 卡戴珊不雅视频在线播放| 国产成人一区二区在线| 99九九线精品视频在线观看视频| 尾随美女入室| 天堂8中文在线网| av视频免费观看在线观看| 2018国产大陆天天弄谢| 在线观看免费视频网站a站| 精品午夜福利在线看| 97在线人人人人妻| 亚洲综合色惰| 丝袜喷水一区| 777米奇影视久久| 乱码一卡2卡4卡精品| 精品国产乱码久久久久久小说| 另类精品久久| 国产精品无大码| 国产亚洲一区二区精品| 在线观看免费日韩欧美大片 | 尾随美女入室| 欧美日韩视频高清一区二区三区二| 国产在线男女| 91久久精品国产一区二区三区| 秋霞伦理黄片| 成人影院久久| 国产伦在线观看视频一区| 国产高清国产精品国产三级| 亚洲国产精品999| 激情五月婷婷亚洲| 日韩成人伦理影院| 亚洲情色 制服丝袜| 久久毛片免费看一区二区三区| 久久久国产一区二区| 亚洲美女视频黄频| 亚洲美女搞黄在线观看| 国产伦精品一区二区三区视频9| 嘟嘟电影网在线观看| 亚洲天堂av无毛| 国产精品伦人一区二区| 哪个播放器可以免费观看大片| 99热这里只有是精品在线观看| 一级,二级,三级黄色视频| 亚洲欧美一区二区三区国产| 蜜臀久久99精品久久宅男| 一级黄片播放器| 这个男人来自地球电影免费观看 | 日韩欧美一区视频在线观看 | 大香蕉久久网| 一区二区三区免费毛片| 久久女婷五月综合色啪小说| 国模一区二区三区四区视频| 欧美三级亚洲精品| 欧美日韩av久久| 国产成人精品无人区| 日本欧美国产在线视频| 国内少妇人妻偷人精品xxx网站| 在线 av 中文字幕| 纵有疾风起免费观看全集完整版| 51国产日韩欧美| 亚洲精品日韩在线中文字幕| 免费看av在线观看网站| 久久精品久久久久久噜噜老黄| 国产精品蜜桃在线观看| 亚洲欧洲国产日韩| 日本91视频免费播放| 桃花免费在线播放| 午夜影院在线不卡| 国国产精品蜜臀av免费| 亚洲国产精品999| 青青草视频在线视频观看| 色5月婷婷丁香| 午夜免费鲁丝| 色94色欧美一区二区| 在线观看免费高清a一片| 国产精品一区二区在线观看99| 在线观看免费高清a一片| 国产精品一区二区在线观看99| videossex国产| 高清午夜精品一区二区三区| av女优亚洲男人天堂| 日韩人妻高清精品专区| 麻豆乱淫一区二区| 日本午夜av视频| 欧美成人精品欧美一级黄| 国产极品粉嫩免费观看在线 | 日韩欧美一区视频在线观看 | 亚洲经典国产精华液单| 日韩欧美精品免费久久| 国产成人精品久久久久久| 国产精品久久久久成人av| 午夜福利在线观看免费完整高清在| 国产老妇伦熟女老妇高清| 男的添女的下面高潮视频| 亚洲欧美日韩另类电影网站| 中文欧美无线码| 久久久久久久大尺度免费视频| 3wmmmm亚洲av在线观看| 国产日韩欧美亚洲二区| 九色成人免费人妻av| 久久久欧美国产精品| 人人澡人人妻人| 不卡视频在线观看欧美| 啦啦啦视频在线资源免费观看| 国产男女超爽视频在线观看| 欧美 亚洲 国产 日韩一| 制服丝袜香蕉在线| av在线观看视频网站免费| 三级国产精品片| 成人黄色视频免费在线看| 亚洲国产成人一精品久久久| 色吧在线观看| 国产av一区二区精品久久| 欧美少妇被猛烈插入视频| kizo精华| 青青草视频在线视频观看| 精品酒店卫生间| av国产久精品久网站免费入址| 高清午夜精品一区二区三区| 国产成人aa在线观看| 国产黄频视频在线观看| 亚洲不卡免费看| 日韩中文字幕视频在线看片| 久久久精品免费免费高清| 久久影院123| 一个人看视频在线观看www免费| 午夜福利视频精品| 大话2 男鬼变身卡| 秋霞在线观看毛片| 又大又黄又爽视频免费| 成人美女网站在线观看视频| 中国国产av一级| 国产在线免费精品| 人妻人人澡人人爽人人| 日韩av不卡免费在线播放| 天堂俺去俺来也www色官网| av免费观看日本| 极品少妇高潮喷水抽搐| 国产亚洲欧美精品永久| a级片在线免费高清观看视频| 天美传媒精品一区二区| 国产男女超爽视频在线观看| 丝袜脚勾引网站| 亚洲欧美精品专区久久| 女的被弄到高潮叫床怎么办| 少妇丰满av| av在线老鸭窝| 成人18禁高潮啪啪吃奶动态图 | 国产女主播在线喷水免费视频网站| 免费在线观看成人毛片| 欧美日韩精品成人综合77777| videos熟女内射| 夫妻午夜视频| 国产成人精品福利久久| 色婷婷久久久亚洲欧美| 街头女战士在线观看网站| 久久精品国产a三级三级三级| 国产成人精品婷婷| 99国产精品免费福利视频| 毛片一级片免费看久久久久| 国产成人精品无人区| 欧美丝袜亚洲另类| av免费在线看不卡| 欧美高清成人免费视频www| 日日摸夜夜添夜夜爱| 最近中文字幕2019免费版| 亚洲欧美一区二区三区国产| 看非洲黑人一级黄片| 国产日韩欧美在线精品| 国产综合精华液| 大陆偷拍与自拍| 最近2019中文字幕mv第一页| 精华霜和精华液先用哪个| tube8黄色片| 久久午夜综合久久蜜桃| av福利片在线观看| 51国产日韩欧美| 日韩伦理黄色片| 日本欧美视频一区| 午夜免费鲁丝| 久久人妻熟女aⅴ| 亚洲精品成人av观看孕妇| 久久久久精品久久久久真实原创| 久久精品熟女亚洲av麻豆精品| 少妇裸体淫交视频免费看高清| 国产av精品麻豆| 亚洲欧美日韩另类电影网站| 一区二区av电影网| 两个人免费观看高清视频 | 三级国产精品欧美在线观看| 国产亚洲一区二区精品| 十分钟在线观看高清视频www | 自拍欧美九色日韩亚洲蝌蚪91 | 人妻人人澡人人爽人人| 日产精品乱码卡一卡2卡三| 国产爽快片一区二区三区| 2018国产大陆天天弄谢| 黄色一级大片看看| 丝袜脚勾引网站| 日韩伦理黄色片| 日韩精品免费视频一区二区三区 | 搡女人真爽免费视频火全软件| 国产免费又黄又爽又色| 久久精品久久精品一区二区三区| 三级国产精品欧美在线观看| 欧美成人精品欧美一级黄| 日韩中文字幕视频在线看片| 精品少妇黑人巨大在线播放| 老司机亚洲免费影院| 免费观看a级毛片全部| 久久国产精品男人的天堂亚洲 | 成人影院久久| 色网站视频免费| 国产乱来视频区| 国产一区二区在线观看av| 另类亚洲欧美激情| 丰满少妇做爰视频| 在线观看免费日韩欧美大片 | 男女免费视频国产| 婷婷色av中文字幕| 亚洲精品色激情综合| 亚洲性久久影院| 大香蕉97超碰在线| 91aial.com中文字幕在线观看| 亚洲国产毛片av蜜桃av| 精品少妇久久久久久888优播| 肉色欧美久久久久久久蜜桃| 日本黄色片子视频| 大码成人一级视频| 黄色配什么色好看| 亚洲人与动物交配视频| 色婷婷av一区二区三区视频| 免费看不卡的av| 亚洲va在线va天堂va国产| av播播在线观看一区| 一级爰片在线观看| 国产av国产精品国产| 高清黄色对白视频在线免费看 | 精品卡一卡二卡四卡免费| 亚洲一区二区三区欧美精品| 亚洲精品国产av蜜桃| 欧美另类一区| .国产精品久久| 国产乱来视频区| 久久久久视频综合| www.av在线官网国产| 麻豆成人午夜福利视频| 日韩伦理黄色片| 少妇被粗大的猛进出69影院 | 亚洲精品乱久久久久久| 国产在线视频一区二区| 欧美精品亚洲一区二区| 亚洲美女视频黄频| 久久亚洲国产成人精品v| 一区二区三区乱码不卡18| 国产一区二区三区av在线| 美女中出高潮动态图| 国产精品久久久久久av不卡| 国产亚洲5aaaaa淫片| 在线观看人妻少妇| 91在线精品国自产拍蜜月| 国产女主播在线喷水免费视频网站| 熟妇人妻不卡中文字幕| av一本久久久久| 一级毛片黄色毛片免费观看视频| 极品教师在线视频| 日本猛色少妇xxxxx猛交久久| 国产精品一区二区在线观看99| a级毛色黄片| 免费不卡的大黄色大毛片视频在线观看| 日本黄色日本黄色录像| 丝袜脚勾引网站| 亚洲av男天堂| 日韩,欧美,国产一区二区三区| 亚洲欧美一区二区三区国产| 99久久中文字幕三级久久日本| 久久久久久久久久成人| 人妻少妇偷人精品九色| 日韩制服骚丝袜av| 国产伦理片在线播放av一区| 亚洲四区av| 日韩伦理黄色片| 在线播放无遮挡| 亚洲精品日本国产第一区| 国产伦精品一区二区三区视频9| 亚洲性久久影院| 国产精品一区www在线观看| 亚洲国产av新网站| 亚洲国产精品999| 久久久久精品性色| 亚洲av.av天堂| 国产欧美亚洲国产| 久久久久国产精品人妻一区二区| 一级毛片 在线播放| 精品人妻熟女毛片av久久网站|