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

    大氣降水與濟南泉水動態(tài)變化的定量關(guān)系研究

    2017-03-08 05:10:56遲光耀邢立亭主恒祥侯新宇邢學(xué)睿
    地下水 2017年1期
    關(guān)鍵詞:泉域雨量站興隆

    遲光耀,邢立亭,主恒祥,侯新宇,相 華,邢學(xué)睿

    (1.濟南大學(xué) 資源與環(huán)境學(xué)院,山東 濟南 250022;2.山東省濟南市水文局,山東 濟南 250013;3.河北地質(zhì)大學(xué),河北 石家莊 050031)

    大氣降水與濟南泉水動態(tài)變化的定量關(guān)系研究

    遲光耀1,邢立亭1,主恒祥1,侯新宇1,相 華2,邢學(xué)睿3

    (1.濟南大學(xué) 資源與環(huán)境學(xué)院,山東 濟南 250022;2.山東省濟南市水文局,山東 濟南 250013;3.河北地質(zhì)大學(xué),河北 石家莊 050031)

    濟南保泉歷經(jīng)40多年歷史,歷年枯水期泉水一直受到斷流威脅。針對濟南泉水動態(tài)變化現(xiàn)狀,采用相關(guān)分析法,研究不同時間尺度下大氣降水與泉水動態(tài)的關(guān)系。研究表明:濟南泉域巖溶地下水的主要的補給來源是大氣降水,泉水位高低主要受降雨大小的影響;大氣降水對泉水位的影響具有滯后性特點;不同時間尺度下累計降水量回歸分析得出興隆雨量站、燕子山雨量站對泉水位影響滯后時間較短;在現(xiàn)狀開采條件下,當(dāng)旬降雨量累計值大于17.58 mm時對泉水位回升有積極的影響。綜合濟南泉域巖溶水流向、各雨量站距泉的距離、降雨對泉水位影響的滯后時間、各雨量站降水大小與泉水位變幅的關(guān)系程度等方面,得出興隆雨量站所在區(qū)域降雨對泉水位的影響較大,是適宜的保泉補源地,能達(dá)到科學(xué)補源保泉的目的。

    保泉;大氣降水;水位變幅;面降雨量;濟南

    地下水的補給來源有大氣降水、地表水、凝結(jié)水、來自其他含水層或含水系統(tǒng)的水等[1],從水文循環(huán)角度來看大氣降水是地下水最主要的補給源。

    我國巖溶分布廣泛[2],南北方的巖溶發(fā)育程度、巖溶水動態(tài)特征都存在一系列差異,其巖溶水對降雨的響應(yīng)也比較敏感[3]。在我國北方典型的巖溶泉域——濟南巖溶泉域,巖溶水是重要的供水水源[4]。根據(jù)巖溶水的動態(tài)特征、氫氧同位素等證實巖溶水的補給來源是大氣降水[5]。受人類活動的影響,地下水動態(tài)特征與大氣降水的相關(guān)性發(fā)生了明顯變化。例如大量開采地下水會引起水位大幅度下降,形成地下水超采漏斗[6],人工回灌又會引起地下水位上升[7],人為工程活動改變?nèi)奶烊粍討B(tài)特征,并且地下水位動態(tài)存在年內(nèi)季節(jié)變化、年際周期性變化及隨機波動等一定的規(guī)律性[8],近些年來,國內(nèi)外很多專家學(xué)者通過各種方法研究大氣降水與地下水動態(tài)的變化規(guī)律。例如,通過建立模型計算得出在不同水文地質(zhì)條件下,降雨形成的徑流流速快慢也影響著年累計地下水補給總量[9];研究在氣候變化的因素下,大氣降水對沖積層地下水位帶來的影響[10];采用小波分析方法可定量評價泉域地下水位動態(tài)及其與降水的關(guān)系[11,12];采用多元線性回歸模型模擬降水對地下水動態(tài)影響進行了分析,為泉水保護提供了依據(jù)[13];還有采用時間序列法、頻譜分析、人工神經(jīng)網(wǎng)絡(luò)、黑箱模型法、灰色分析、系統(tǒng)動態(tài)神經(jīng)網(wǎng)絡(luò)模式、核函數(shù)法等[14-20],這些研究提高了濟南泉域的研究精度,但是濟南泉水依然受到泉水?dāng)嗔鞯耐{,從工程實踐角度,現(xiàn)有人工補源工程效果不明顯、具有很大盲目性。本文在前人工作的基礎(chǔ)上,研究不同時間尺度下大氣降雨與泉水動態(tài)的定量關(guān)系,找出合適的回灌補源區(qū)域,達(dá)到科學(xué)保泉的目的。

    1 研究區(qū)概況

    濟南市地處魯中山地的北緣,地形南高北低,南部為前震旦系片麻巖組成的結(jié)晶基底,寒武系~奧陶系地層由老到新向北傾斜出露,燕山期巖漿巖體分布在北部,在這一特定的地形、地質(zhì)構(gòu)造條件下,巖溶水在南部山區(qū)得到大氣降水的補給后,其運動方向與地形坡向和地層傾向方向大體一致,自南向北運移,巖溶水向北運動遇到巖漿巖體而受阻,地下水富集,在適宜的地形和有利構(gòu)造部位承壓水以泉的形式出露(圖1)。寒武紀(jì)張夏組石灰?guī)r含水層與奧陶紀(jì)石灰?guī)r含水層之間存在崮山組和長山組泥質(zhì)灰?guī)r與頁巖,削弱了寒武紀(jì)與奧陶紀(jì)含水層的水力聯(lián)系,單斜自流構(gòu)造內(nèi)部,有多條斷裂構(gòu)造切割,加快地下水的徑流,因此,不同地段降水對泉水位的影響程度不一。

    注:1.碳酸鹽巖類裂隙巖溶水,>5 000(單井涌水量,m3/d,下同);2.碳酸鹽巖類裂隙巖溶水,1 000~5 000;3.碳酸鹽巖類裂隙巖溶水水,500~1 000;4.碳酸鹽巖類裂隙巖溶水,<500;5.變質(zhì)巖、巖漿巖;6.石炭、二疊砂頁巖;7.間接補給區(qū);8.直接補給區(qū)與間接補給區(qū)分區(qū)線;9.斷層

    圖1 研究區(qū)水文地質(zhì)圖及雨量站分布位置圖

    2 數(shù)據(jù)采集

    2.1 雨量站

    選擇泉域直接和間接補給區(qū)內(nèi)15個雨量站,主要分布在紹而、東紅廟、興隆、燕子山、黃臺橋、邱家莊、臥虎山、窩鋪、南高而、崮山、東梧、劉家莊、王家莊、棗林、西營等地。降水資料時長為2009—2015年,觀測頻率為日觀測數(shù)據(jù)。

    2.2 泉水位

    選取2009—2015年的趵突泉水位觀測資料,觀測頻率為日觀測數(shù)據(jù)。數(shù)據(jù)處理采用excel軟件。

    3 結(jié)果分析與討論

    3.1 各站降水特征

    濟南市降雨表現(xiàn)在空間上的不均一性,位于研究區(qū)東南部的邱家莊、西營、東梧雨量站2015年平均降雨量約466.67 mm,研究區(qū)中部的東紅廟、紹而、興隆雨量站2015年平均降雨量570 mm左右,研究區(qū)西北部的劉家莊雨量站2015年平均降雨量481 mm(圖2),補給區(qū)的降水量多于排泄區(qū)。由于濟南泉域不同地點的地層巖性入滲能力不同,因此不同站點的雨量對趵突泉水位影響具有差異性。

    圖2 2015年研究區(qū)降雨量等值線圖

    3.2 泉水位變化特征

    根據(jù)歷年觀測資料,泉水位的升高與降雨量的關(guān)系比較明顯(圖3和圖4),在豐水期7~9月份水位上升,在枯水期3~6月份水位下降,泉水位受大氣降水的影響比較顯著,泉水位變化具有典型的季節(jié)性。

    圖3 2009—2015年旬降雨量值與泉旬末水位的關(guān)系

    圖4 2009—2015年五日降雨量值與泉五日末水位的關(guān)系

    為進一步提高研究精度選用旬作為研究降雨量與泉水位差關(guān)系的時間尺度,通過旬時間尺度降雨量和泉水位差的散點圖得出線性方程與x軸的交點為17.58 mm,即在現(xiàn)狀開采條件下旬降雨量累計值大于17.58 mm時,可引起泉水位回升;線性方程斜率為0.003 3(圖5),這一結(jié)果與前人在年尺度和月尺度條件下得出的降雨量和泉水位的線性方程的斜率0.003 3和0.003 8基本接近[7],而五日時間尺度降雨量與泉水位差的相關(guān)系數(shù)較差(圖6)。

    圖5 2009~2015年水位旬變化量與旬降雨量的關(guān)系

    圖6 2009~2015年水位五日變化量與五日降雨量的關(guān)系

    3.3 水位與降水補給的滯后時間分析

    多年水位監(jiān)測表明泉水動態(tài)相對于降水有一定滯后(圖3、圖4),但通過不同時間尺度大小的降水產(chǎn)生滯后效應(yīng)是有差別的,以年尺度、月尺度降水來研究水位變化會掩蓋水位的微觀變化,因此,從旬尺度、五日和次降雨量來探究降雨量對泉水位的影響。

    以旬尺度為例設(shè)計滯后時間分別為1 d、2 d、3 d、4 d、5 d、6 d,泉水位變幅采用旬初日至第11、12、13、14、15、16日的水位差值,分別計算旬時間尺度降水量與6個泉水位變幅的相關(guān)系數(shù)。五日降雨和次降雨采用與旬降雨相同的計算方法。得到不同時間尺度下水位滯后時間較短、對泉水位影響較大的雨量站(表1)。

    表1 不同時間尺度下對泉水位影響較大的雨量站

    (注:雨量站排列按照相關(guān)系數(shù)由大到小、滯后時間由短到長的順序排列)

    從表1可以看出次降雨量尺度下相關(guān)系數(shù)整體最好,所以選取次降雨量尺度下對水位與降水補給的滯后時間的影響因素進行分析。在次降雨量尺度下,將各雨量站的最佳滯后時間與距泉的距離做相關(guān)性分析(圖7),可以看出兩者總體上是正相關(guān)的關(guān)系,即距泉的距離越遠(yuǎn)滯后時間越長,但也存在異常情況如黃臺橋和劉家莊雨量站,二者滯后時間相同,但是距離不同,由此可以推斷出距離只是滯后時間大小的主導(dǎo)因素之一,還與各雨量站的巖溶、裂隙發(fā)育程度,巖溶水的流向及與泉群的連通性有關(guān)。雨量站所處的地層不同也是影響其滯后原因之一 (圖1),雨量站分別分布在奧陶紀(jì)區(qū)、寒武紀(jì)區(qū),通過圖1及表2,可以看出在次時間尺度下,滯后時間短且相關(guān)系數(shù)大的雨量站如興隆、燕子山雨量站都在奧陶紀(jì)區(qū)內(nèi),距離泉群較近,說明奧陶紀(jì)區(qū)內(nèi)雨量站與泉群的關(guān)系比較密切。

    表2 次降雨尺度下各雨量站滯后時間及地層分布統(tǒng)計表

    (注:括號內(nèi)為相關(guān)系數(shù))

    圖7 次降雨量下滯后時間與距泉距離的關(guān)系

    圖8 研究區(qū)雨量站的泰森多邊形

    3.4 泉水位與各站面雨量之間的關(guān)系

    利用泰森多邊形法[21]對研究區(qū)進行分區(qū)后進行面雨量與水位關(guān)系研究(圖8)。以旬尺度為例,已知旬尺度下各雨量站的滯后時間,為了消除降雨滯后的影響,即選取各個雨量站的最佳滯后時間,計算各雨量站旬面雨量累計值與其滯后時間水位差的相關(guān)關(guān)系,分析各雨量站降水與泉水位的關(guān)系,五日面雨量與次面雨量采用與旬面雨量相同的計算方法。計算得到不同時間尺度下對泉水位影響最相關(guān)的面雨量站為興隆站和東紅廟站。

    表3 不同時間尺度下的最佳面雨量站統(tǒng)計表

    回歸分析表明(表3),在相同降雨量的條件下,這兩個雨量站的降水對泉水位的增幅較大。通過比較不同時間尺度下各雨量站面雨量與泉水位差的關(guān)系,發(fā)現(xiàn)興隆雨量站對泉水位的影響比較大,適宜作為回灌補源地。

    3.5 降水與泉水位關(guān)系揭示的水文地質(zhì)問題

    根據(jù)前人對濟南泉域巖溶水的流向研究分析,泉域巖溶水的徑流方向和徑流強度受地形、地貌、巖性和地質(zhì)構(gòu)造等因素控制,其運動方向與地形及巖層的傾斜方向大體一致,在接受大氣降水滲入和地表水滲漏補給后總體方向由南向北運動。研究區(qū)的各雨量站距泉的位置不同、接受降雨補給的面積不同、各雨量站與泉群的連通性不同,所以降水對泉水位的影響不一。

    通過上述資料數(shù)據(jù)分析得出:(1)位于寒武紀(jì)張夏地層內(nèi)的臥虎山、邱家莊、南高而等雨量站距離泉群較遠(yuǎn),滯后時間較長;奧陶紀(jì)區(qū)內(nèi)的興隆、燕子山雨量站距離泉群較近,滯后時間短,因此濟南回灌補源地應(yīng)當(dāng)選在奧陶紀(jì)區(qū)內(nèi);(2)在相同降雨強度和距離泉群相同距離的前提下,東紅廟和興隆雨量站的滯后時間不同,并且興隆雨量站的滯后時間短于東紅廟,由此確定興隆雨量站與泉群的連通性好于東紅廟,因此興隆地區(qū)更利于回灌;(3)從面降雨量相關(guān)系數(shù)來看,在相同降雨強度和距離泉群相同距離的前提下,興隆雨量站的相關(guān)系數(shù)大于東紅廟雨量站,也能準(zhǔn)確地揭示興隆雨量站的巖溶發(fā)育程度、與泉群的連通性好于東紅廟。

    興隆雨量站位于興隆強滲漏帶范圍內(nèi),是濟南南部山區(qū)重要的生態(tài)和水源涵養(yǎng)區(qū)、泉水直接補給區(qū),并且興隆強滲漏帶是距離四大泉群直線距離最近、面積最大的強滲漏帶,所以綜合面降雨量與泉水位相關(guān)系數(shù)、滯后時間和距泉群的距離等數(shù)據(jù),可以得出興隆地區(qū)宜作為泉水回灌補源地。

    4 結(jié)語

    (1) 濟南泉域巖溶水的最主要補給來源是大氣降水,泉水位的高低受降雨量的影響較大。

    (2) 綜合旬、五日和次降雨量以及濟南泉域地下水的流動方向、各雨量站距泉的距離、降雨滯后時間、面降雨量和水位差的相關(guān)關(guān)系等方面,得出興隆雨量站、燕子山雨量站較其他雨量站降水滯后時間較短,有利于回灌。

    (3) 通過研究不同區(qū)域雨量站與泉水位的關(guān)系,得出奧陶紀(jì)區(qū)內(nèi)雨量站與泉群的關(guān)系密切,濟南回灌補源應(yīng)當(dāng)選在奧陶紀(jì)區(qū)內(nèi)而非寒武紀(jì)區(qū),因此興隆雨量站區(qū)域宜作為泉水回灌補源地。

    (4) 在枯水期維持濟南泉水動態(tài)穩(wěn)定的旬降雨量累計值大于17.58 mm時,與前人研究相比,提高了濟南降水與泉水位響應(yīng)的研究精度。

    [1]王大純.水文地質(zhì)學(xué)基礎(chǔ)[M]. 北京:地質(zhì)出版社.1995:100-101.

    [2]袁道先.中國巖溶學(xué)[M]. 北京:地質(zhì)出版社.1993:1-8.

    [3]Ellen K H, Laura T, William B W. Quantifying the place of karst aquifers in the groundwater to surf-ace water continuum:A time series analysis study of storm behavior in Pennsylvania water resources[J]. Journal of Hydrology, 2009,376:307-317.

    [4]祁曉凡,李文鵬,李海濤,等.濟南巖溶泉域地下水位、降水、氣溫與大尺度氣象遙相關(guān)[J]. 水文地質(zhì)工程地質(zhì).2015,42(6):18-28.

    [5]萬利勤.濟南泉域巖溶地下水的示蹤研究[D].北京:中國地質(zhì)大學(xué).2008.

    [6]Yong Yang,Guo-Min Li,Yan-Hui Dong,et al.Influece of South to North Water Transfer on groundwater dynamic change in Beijing plain[J].Environ Earth S-ci,2012,(65):1323-1331.

    [7]周娟,邢立亭,滕朝霞.制約濟南巖溶大泉持續(xù)噴涌的主因素閾值研究[J].華東師范大學(xué)學(xué)報(自然科學(xué)版).2015,(3):145-156.

    [8]FAN Y H,HUO X L,HAO Y H,et al.An assembled extreme value statistical model of karst spring dis-charge[J].Journal of Hydrology,2013,504(11):57-68.

    [9]Orhan Gunduz,Celalettin Simsek. Influence of Clim-ate Change on Shallow Groundwater Resources:The Link Between Precipitation and Groundwater Leve-ls in Alluvial Systems[J].Climate Change and its E-ffects on Water Resources,2011,3:225-233.

    [10]Noam Zach Dvory,Yakov Livshitz, Eilon Adar, et al.The effect of hydrogeological conditions on variabi-lity and dynamic of groundwater recharge in a car-bonate aquifer at local scale[J]. Journal of Hydrology,2016,535:480-494.

    [11]祁曉凡,楊麗芝,韓曄.濟南泉域地下水位動態(tài)及其對降雨響應(yīng)的交叉小波分析[J].地球科學(xué)進展.2012,27(9):970-978.

    [12]DAMIR J,VESNA D J.A frequency domain approa-ch to groundwater recharge estimation in karst[J]. Journal of Hydrology,2004,289:95-110.

    [13]張建芝,邢立亭.回歸分析法在地下水動態(tài)分析中的應(yīng)用[J].地下水.2010,32(4):88-90.

    [14]Wen Xi Lu,Ying Zhao,Hai Bo Chu,et al. The anal-ysis of groundwater levels influenced by dual facto-rs in western Jilin Province by using time series a-nalysis method[J]. Applied Water Science, 2014,3(4):251-260.

    [15]陳葆仁,洪再吉,汪福炘,等.地下水動態(tài)及其預(yù)測[M].北京:科學(xué)出版社.1988.

    [16]張建芝.濟南泉域地下水動態(tài)特征及監(jiān)測網(wǎng)優(yōu)化[D].濟南:濟南大學(xué).2011.

    [17]Jan C D,Chen T H,Huang H M.Analysis of rainfal-linduced quick groundwater-level response by using a Kernel function[J]. Paddy Water Environment, 2013(11): 135-144.

    [18]Fi-John Chang,Pin-An Chen,Chen-Wuing Liu,et al. Regional estimation of groundwater arsenic concent-rations through systematical dynamic-neural modeling[J].Journal of Hydrology,2013,499:265-274.

    [19]胡克禎,張建芝,邢立亭.基于時間序列分析的地下水動態(tài)研究[J].水科學(xué)與工程技術(shù).2011(5):32-34.

    [20]Mao X M,Shang S H,Liu X.Groundwater level pre-dictions using artificial neural networks[J]. Tsinghua Science and Technology,2002,7(6):574-579.

    [21]龍胤慧,廖梓龍.基于泰森多邊形法的慶陽市面雨量計算[J].河北工程大學(xué)學(xué)報(自然科學(xué)版).2012,29(3):64-67.

    The study of quantitative relationship between the spring water and the dynamic change of the atmospheric precipitation in Ji'nan

    CHI Guang-yao1,XING Li-ting1,ZHU Heng-xiang1,HOU Xin-yu1,XIANG Hua2,XING Xue-rui3

    (1.Resource and Environment College of Jinan University, Jinan 250022,China;2.Hydrographic Office of Jinan City, Jinan, 250013, Shandong ,China;3.Hebei Geo university, Shijiazhuang 050031,China)

    The protection of Ji'nan spring has history more than 40 years, and the spring during dry season over the years has been threatened by cutoff. According to current situation of dynamic change, this paper studies by the correlation analysis method that the relationship between precipitation and spring water in different time scales. Research shows that: The main recharge sources of karst groundwater in Ji'nan spring is precipitation, and the water level of the spring is mainly affected by rainfall. The influence of precipitation on the water level of the spring has the characteristic of hysteresis. The regression analysis of the cumulative precipitation in different time scales shows that the lag time of the Xinglong rainfall station and the Yanzishan rainfall station affect on the water level of the spring is short. Under the condition of the current exploitation, it has a positive effect on the recovery of spring water level when the cumulative value of natural precipitation is more than 17.58mm. The study that the flow direction of the karst water in the spring of Ji'nan, the distance between the spring and each rainfall station, the lag time of the rainfall affect on the spring water level, the relationship between the precipitation and the change of water level of the springs shows that, it is in the region of the Xinglong rainfall station that the influence of precipitation on the water level of the spring is great. Therefore, the Xinglong rainfall station is a suitable recharge place for spring protection and can achieve the purpose of scientific compensation for spring protection.

    hoizumi;atmospheric precipitation;range of water level;surface rainfall;Jinan

    2016-10-17

    國家自然科學(xué)基金:內(nèi)陸平原區(qū)粘土對淺層地下水與中深層咸水循環(huán)交替的阻滯作用(41172222);淺層大

    遲光耀(1993-),男,山東濰坊人,在讀碩士研究生。主攻方向:地下水動態(tài)演化。

    邢立亭(1966-),男,山東青島人,教授,主要從事地下水環(huán)境演化研究。

    P641.2

    A

    1004-1184(2017)01-0008-04

    陸演化咸水水文地質(zhì)參數(shù)變異機制研究(41472216)

    猜你喜歡
    泉域雨量站興隆
    遼西山洪災(zāi)害重點防治區(qū)遙測雨量站點優(yōu)化布設(shè)研究
    山西省泉域水資源保護條例
    山西水利(2022年9期)2022-03-12 09:03:35
    山西省人大常委會舉辦《山西省泉域水資源保護條例》新聞發(fā)布會
    山西水利(2022年9期)2022-03-12 09:03:35
    信息熵方法在遼寧省不同分區(qū)雨量站網(wǎng)布設(shè)的應(yīng)用研究
    黑龍洞泉域子單元劃分及特征分析
    特別的生日禮物
    中小流域雨量站網(wǎng)密度規(guī)劃與研究
    敦煌月牙泉水補給及泉域地下水動態(tài)研究
    雨量站網(wǎng)測量精度的評估
    我要飛翔
    亚洲av男天堂| 国产成人精品福利久久| 卡戴珊不雅视频在线播放| 中国美白少妇内射xxxbb| 高清av免费在线| 免费av不卡在线播放| 日本91视频免费播放| 亚洲图色成人| 久久6这里有精品| 日本vs欧美在线观看视频 | 国产亚洲91精品色在线| freevideosex欧美| 看免费成人av毛片| 下体分泌物呈黄色| 欧美日韩视频高清一区二区三区二| 精品国产乱码久久久久久小说| 秋霞伦理黄片| 中文字幕精品免费在线观看视频 | 久久久a久久爽久久v久久| 久久午夜福利片| 少妇的逼好多水| 欧美精品一区二区免费开放| 波野结衣二区三区在线| 在线观看免费视频网站a站| 91精品国产九色| 91久久精品国产一区二区成人| 一级二级三级毛片免费看| 亚洲国产精品一区三区| 亚洲国产精品一区三区| 男人舔奶头视频| 国产亚洲精品久久久com| a级毛片免费高清观看在线播放| 我要看日韩黄色一级片| 一区在线观看完整版| 看免费成人av毛片| 少妇精品久久久久久久| 99热6这里只有精品| 男男h啪啪无遮挡| 国产亚洲91精品色在线| 五月伊人婷婷丁香| 永久网站在线| 亚洲国产精品999| 精品亚洲成a人片在线观看| 欧美3d第一页| videossex国产| 99久久综合免费| 大香蕉97超碰在线| 在线免费观看不下载黄p国产| 亚洲欧美中文字幕日韩二区| 久久99热这里只频精品6学生| 色吧在线观看| 亚洲av成人精品一二三区| 久久久久久伊人网av| 久久久久视频综合| 这个男人来自地球电影免费观看 | 亚洲成人手机| 国模一区二区三区四区视频| 精品国产国语对白av| 久久国产亚洲av麻豆专区| 一区二区三区乱码不卡18| 免费看日本二区| 美女视频免费永久观看网站| 青春草国产在线视频| 国产69精品久久久久777片| 欧美日韩亚洲高清精品| 新久久久久国产一级毛片| 精品人妻一区二区三区麻豆| 日日啪夜夜爽| 免费黄频网站在线观看国产| 亚洲不卡免费看| 国产乱来视频区| 热re99久久精品国产66热6| 亚洲美女搞黄在线观看| 久久久久久久亚洲中文字幕| 精品人妻一区二区三区麻豆| 国产精品一二三区在线看| 日韩亚洲欧美综合| 啦啦啦中文免费视频观看日本| 国产亚洲欧美精品永久| xxx大片免费视频| 日韩一本色道免费dvd| 日韩av在线免费看完整版不卡| 少妇人妻一区二区三区视频| 少妇被粗大猛烈的视频| 我的老师免费观看完整版| 国产精品不卡视频一区二区| 在线精品无人区一区二区三| 国产精品99久久久久久久久| .国产精品久久| 国产av码专区亚洲av| 大话2 男鬼变身卡| 日产精品乱码卡一卡2卡三| 少妇被粗大猛烈的视频| 各种免费的搞黄视频| 少妇的逼好多水| 男女国产视频网站| 天美传媒精品一区二区| 久久97久久精品| 精品国产国语对白av| 国产淫片久久久久久久久| 97精品久久久久久久久久精品| 国产精品伦人一区二区| 99久久人妻综合| 老熟女久久久| 国产亚洲一区二区精品| 国产美女午夜福利| 久久久久国产精品人妻一区二区| 建设人人有责人人尽责人人享有的| 熟妇人妻不卡中文字幕| av在线老鸭窝| 欧美人与善性xxx| 日日摸夜夜添夜夜添av毛片| 青春草视频在线免费观看| 丁香六月天网| 免费黄网站久久成人精品| 免费看av在线观看网站| 少妇人妻精品综合一区二区| 国产精品久久久久久精品电影小说| 午夜老司机福利剧场| 美女福利国产在线| 免费观看无遮挡的男女| 亚洲精品久久午夜乱码| 午夜福利视频精品| 国产亚洲精品久久久com| 性色av一级| 狂野欧美白嫩少妇大欣赏| 欧美 亚洲 国产 日韩一| 人人妻人人添人人爽欧美一区卜| 观看美女的网站| 偷拍熟女少妇极品色| 五月玫瑰六月丁香| 日韩不卡一区二区三区视频在线| 日韩av在线免费看完整版不卡| 成人亚洲精品一区在线观看| 国产伦理片在线播放av一区| 免费看不卡的av| kizo精华| 午夜免费男女啪啪视频观看| 久久毛片免费看一区二区三区| 蜜桃在线观看..| 日韩,欧美,国产一区二区三区| 国产av精品麻豆| 亚洲性久久影院| 赤兔流量卡办理| 国产淫语在线视频| 日韩,欧美,国产一区二区三区| 99re6热这里在线精品视频| 国产乱人偷精品视频| av线在线观看网站| 亚洲国产成人一精品久久久| 极品教师在线视频| 特大巨黑吊av在线直播| 丝袜在线中文字幕| 最近的中文字幕免费完整| 女性被躁到高潮视频| 国产视频内射| 日韩亚洲欧美综合| 亚洲精品第二区| 亚洲欧美一区二区三区国产| 免费观看的影片在线观看| 99久久精品国产国产毛片| 高清欧美精品videossex| 国产精品熟女久久久久浪| 国产片特级美女逼逼视频| 成年人午夜在线观看视频| 熟女人妻精品中文字幕| 欧美xxxx性猛交bbbb| 国产精品国产三级国产专区5o| 国产欧美日韩一区二区三区在线 | 日本av手机在线免费观看| 久久久久网色| 亚洲图色成人| 国产色爽女视频免费观看| 日本91视频免费播放| 免费少妇av软件| 人妻 亚洲 视频| 一级毛片黄色毛片免费观看视频| 婷婷色麻豆天堂久久| 观看av在线不卡| 我的女老师完整版在线观看| 国产亚洲最大av| 女性被躁到高潮视频| 91精品国产九色| av线在线观看网站| 18+在线观看网站| 在线看a的网站| 欧美日韩一区二区视频在线观看视频在线| 国产黄色免费在线视频| 一边亲一边摸免费视频| 国产亚洲最大av| 黄色毛片三级朝国网站 | 亚洲成色77777| 亚洲精品日韩av片在线观看| 噜噜噜噜噜久久久久久91| 久久综合国产亚洲精品| 国精品久久久久久国模美| 国产免费一区二区三区四区乱码| 一本大道久久a久久精品| 亚洲欧美中文字幕日韩二区| 秋霞在线观看毛片| 一级av片app| 中文字幕制服av| 成人综合一区亚洲| 永久网站在线| 九九久久精品国产亚洲av麻豆| 一区二区av电影网| 欧美日韩视频高清一区二区三区二| 欧美 亚洲 国产 日韩一| 国产91av在线免费观看| 高清视频免费观看一区二区| 日韩成人伦理影院| 亚洲经典国产精华液单| 免费观看在线日韩| 国产亚洲5aaaaa淫片| 一区二区av电影网| 亚洲欧洲日产国产| 午夜免费鲁丝| 亚洲四区av| 久久精品国产亚洲网站| 内地一区二区视频在线| 日韩成人伦理影院| 久久久久久久久久久免费av| 一级毛片黄色毛片免费观看视频| 伊人亚洲综合成人网| 国产真实伦视频高清在线观看| 日本wwww免费看| 日韩av免费高清视频| 91久久精品国产一区二区三区| 欧美性感艳星| 99九九在线精品视频 | 久久久久国产网址| 全区人妻精品视频| 久久免费观看电影| 日本色播在线视频| av免费在线看不卡| 久久久久人妻精品一区果冻| 91精品国产国语对白视频| 大香蕉97超碰在线| 久久99热6这里只有精品| 老熟女久久久| 午夜视频国产福利| 国产精品一区二区三区四区免费观看| 国产av精品麻豆| 亚洲天堂av无毛| 高清av免费在线| av不卡在线播放| 久久久国产欧美日韩av| 美女脱内裤让男人舔精品视频| 寂寞人妻少妇视频99o| 欧美精品人与动牲交sv欧美| 五月伊人婷婷丁香| 黑人高潮一二区| 中文字幕精品免费在线观看视频 | 日日撸夜夜添| 久久国产亚洲av麻豆专区| 99九九在线精品视频 | 日本黄大片高清| 少妇被粗大猛烈的视频| 亚洲精品久久午夜乱码| 国产老妇伦熟女老妇高清| 99热这里只有精品一区| 五月开心婷婷网| 亚洲欧美成人综合另类久久久| 激情五月婷婷亚洲| 乱系列少妇在线播放| 欧美97在线视频| 在线看a的网站| 日韩精品有码人妻一区| 免费av中文字幕在线| av不卡在线播放| 又黄又爽又刺激的免费视频.| 精品人妻熟女毛片av久久网站| 国产永久视频网站| 桃花免费在线播放| 亚洲婷婷狠狠爱综合网| 高清毛片免费看| 少妇高潮的动态图| 如日韩欧美国产精品一区二区三区 | 日韩不卡一区二区三区视频在线| 一级毛片电影观看| 一区二区三区免费毛片| 99热6这里只有精品| 中文在线观看免费www的网站| 尾随美女入室| 午夜福利在线观看免费完整高清在| 久久国产乱子免费精品| 2018国产大陆天天弄谢| 国产精品久久久久久av不卡| 免费在线观看成人毛片| 老女人水多毛片| 欧美日本中文国产一区发布| 美女大奶头黄色视频| www.色视频.com| 欧美日韩视频高清一区二区三区二| 男人和女人高潮做爰伦理| 久久精品国产亚洲av涩爱| 亚洲欧美日韩东京热| 欧美精品人与动牲交sv欧美| 亚洲国产av新网站| 国产日韩一区二区三区精品不卡 | 国产黄频视频在线观看| 欧美精品国产亚洲| 国产成人aa在线观看| 亚洲欧美日韩东京热| 少妇人妻 视频| 精品国产国语对白av| 99热6这里只有精品| 午夜影院在线不卡| 伦理电影大哥的女人| 黄片无遮挡物在线观看| 91在线精品国自产拍蜜月| 欧美成人午夜免费资源| 亚洲欧美日韩东京热| av天堂中文字幕网| 3wmmmm亚洲av在线观看| 久久午夜综合久久蜜桃| 国产精品一区www在线观看| 涩涩av久久男人的天堂| 免费大片黄手机在线观看| 亚洲精品自拍成人| 一级片'在线观看视频| 精品久久久精品久久久| 国产91av在线免费观看| 高清视频免费观看一区二区| 男女啪啪激烈高潮av片| 黑人高潮一二区| 亚洲va在线va天堂va国产| 国产一级毛片在线| 99国产精品免费福利视频| 婷婷色综合大香蕉| 亚洲婷婷狠狠爱综合网| 丝袜喷水一区| 天堂中文最新版在线下载| 精品国产乱码久久久久久小说| 亚洲国产精品一区三区| 麻豆乱淫一区二区| 成人午夜精彩视频在线观看| 日本与韩国留学比较| av国产精品久久久久影院| 51国产日韩欧美| 亚洲国产毛片av蜜桃av| 插逼视频在线观看| 九九久久精品国产亚洲av麻豆| 一区二区三区四区激情视频| 麻豆精品久久久久久蜜桃| 久久狼人影院| 午夜福利在线观看免费完整高清在| 日本vs欧美在线观看视频 | 人妻 亚洲 视频| 欧美精品人与动牲交sv欧美| 国产男人的电影天堂91| 一本大道久久a久久精品| 又爽又黄a免费视频| 波野结衣二区三区在线| 亚洲欧洲国产日韩| 少妇猛男粗大的猛烈进出视频| 婷婷色av中文字幕| 久久久久久久久久成人| 欧美亚洲 丝袜 人妻 在线| 在线观看美女被高潮喷水网站| 国产爽快片一区二区三区| 在线看a的网站| a级毛片免费高清观看在线播放| 在线亚洲精品国产二区图片欧美 | 蜜桃久久精品国产亚洲av| 国产欧美亚洲国产| 日韩不卡一区二区三区视频在线| av不卡在线播放| 午夜福利,免费看| 免费观看a级毛片全部| 国产精品女同一区二区软件| 韩国高清视频一区二区三区| 国产精品国产av在线观看| 国产黄色免费在线视频| 国产高清国产精品国产三级| 韩国av在线不卡| 欧美老熟妇乱子伦牲交| 2022亚洲国产成人精品| 国产视频首页在线观看| 人妻系列 视频| 高清视频免费观看一区二区| 午夜av观看不卡| av国产精品久久久久影院| 中文字幕制服av| 纵有疾风起免费观看全集完整版| 国产精品免费大片| 久久精品国产亚洲av天美| 日韩熟女老妇一区二区性免费视频| 精品国产一区二区久久| 精品国产国语对白av| 91成人精品电影| 国产免费福利视频在线观看| 日韩av免费高清视频| 狂野欧美白嫩少妇大欣赏| 日本黄色日本黄色录像| 最新的欧美精品一区二区| 久久人人爽av亚洲精品天堂| 特大巨黑吊av在线直播| 亚洲国产精品一区二区三区在线| 国产黄片视频在线免费观看| 久久久a久久爽久久v久久| 国产成人精品无人区| 99久久中文字幕三级久久日本| 国语对白做爰xxxⅹ性视频网站| 国产乱人偷精品视频| 交换朋友夫妻互换小说| 人人妻人人添人人爽欧美一区卜| 九九爱精品视频在线观看| 黄色日韩在线| 日韩大片免费观看网站| 99久国产av精品国产电影| av专区在线播放| 一区二区三区精品91| 亚洲欧美日韩卡通动漫| 高清在线视频一区二区三区| 亚洲欧美一区二区三区黑人 | 久久久久久伊人网av| 色婷婷久久久亚洲欧美| 午夜免费男女啪啪视频观看| 日本av免费视频播放| 午夜福利视频精品| 国产精品麻豆人妻色哟哟久久| av网站免费在线观看视频| 欧美三级亚洲精品| 亚洲一级一片aⅴ在线观看| 亚洲欧美中文字幕日韩二区| 久久亚洲国产成人精品v| 国产一区二区在线观看日韩| 美女内射精品一级片tv| 精品一区二区三卡| av在线观看视频网站免费| 国精品久久久久久国模美| 老女人水多毛片| 日韩av不卡免费在线播放| 久久久久久久精品精品| 十八禁高潮呻吟视频 | 国产成人免费无遮挡视频| 高清毛片免费看| freevideosex欧美| 高清欧美精品videossex| 亚洲欧美成人精品一区二区| 亚洲综合色惰| 男人狂女人下面高潮的视频| 波野结衣二区三区在线| 一本色道久久久久久精品综合| 亚洲人与动物交配视频| 五月开心婷婷网| 自拍欧美九色日韩亚洲蝌蚪91 | 22中文网久久字幕| 美女内射精品一级片tv| 色5月婷婷丁香| 午夜久久久在线观看| 高清毛片免费看| 中文天堂在线官网| 综合色丁香网| 国产伦精品一区二区三区四那| www.色视频.com| 日本黄色片子视频| 美女脱内裤让男人舔精品视频| av一本久久久久| 久久久国产欧美日韩av| 3wmmmm亚洲av在线观看| 亚洲第一区二区三区不卡| 日韩欧美 国产精品| 丰满迷人的少妇在线观看| 如何舔出高潮| 精品国产一区二区三区久久久樱花| 多毛熟女@视频| 国产爽快片一区二区三区| 国产有黄有色有爽视频| 有码 亚洲区| 99热这里只有精品一区| 午夜福利在线观看免费完整高清在| 午夜免费鲁丝| 精品亚洲成a人片在线观看| 欧美日韩视频精品一区| 亚洲欧美中文字幕日韩二区| 黑丝袜美女国产一区| 亚洲情色 制服丝袜| 亚洲国产毛片av蜜桃av| 国产极品天堂在线| 99九九在线精品视频 | 一级毛片 在线播放| 狠狠精品人妻久久久久久综合| 91aial.com中文字幕在线观看| 黄色视频在线播放观看不卡| 国产精品.久久久| 精品一区在线观看国产| 精品亚洲乱码少妇综合久久| freevideosex欧美| 美女内射精品一级片tv| 日韩在线高清观看一区二区三区| www.色视频.com| 国产精品一二三区在线看| 一级毛片黄色毛片免费观看视频| 久久综合国产亚洲精品| 国产精品麻豆人妻色哟哟久久| 如日韩欧美国产精品一区二区三区 | 国产伦在线观看视频一区| 亚洲精品色激情综合| 如何舔出高潮| 精品人妻偷拍中文字幕| 亚洲婷婷狠狠爱综合网| 色5月婷婷丁香| 啦啦啦在线观看免费高清www| 亚洲av不卡在线观看| 日韩一本色道免费dvd| 能在线免费看毛片的网站| 欧美区成人在线视频| 成人无遮挡网站| 色94色欧美一区二区| 男男h啪啪无遮挡| 18+在线观看网站| 久久久久视频综合| 精品99又大又爽又粗少妇毛片| 亚洲激情五月婷婷啪啪| 中文字幕精品免费在线观看视频 | 人人澡人人妻人| 亚洲丝袜综合中文字幕| 久久久久久伊人网av| 丝袜在线中文字幕| 日韩一本色道免费dvd| 日韩免费高清中文字幕av| 欧美日本中文国产一区发布| 一级毛片电影观看| 日韩视频在线欧美| 国内揄拍国产精品人妻在线| 热re99久久精品国产66热6| 日韩一本色道免费dvd| 丝袜脚勾引网站| 寂寞人妻少妇视频99o| 99久久精品国产国产毛片| 夫妻性生交免费视频一级片| 国产免费视频播放在线视频| 久久鲁丝午夜福利片| 一级毛片久久久久久久久女| 精品一区二区三卡| 亚洲一级一片aⅴ在线观看| 国产深夜福利视频在线观看| 69精品国产乱码久久久| 三级国产精品片| 亚洲精品自拍成人| 曰老女人黄片| 制服丝袜香蕉在线| 亚洲精品国产av成人精品| 国产精品一区二区在线观看99| 久久97久久精品| 十八禁网站网址无遮挡 | 国产高清三级在线| 人人妻人人澡人人看| 91精品国产九色| 91aial.com中文字幕在线观看| 国产无遮挡羞羞视频在线观看| 国产 一区精品| a级毛色黄片| 亚洲精品日韩在线中文字幕| 曰老女人黄片| 大话2 男鬼变身卡| 一本久久精品| 国产成人午夜福利电影在线观看| 亚洲av福利一区| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 免费看日本二区| 国产黄频视频在线观看| 男人爽女人下面视频在线观看| 插逼视频在线观看| videossex国产| 成人免费观看视频高清| 午夜免费鲁丝| 国产一区二区三区综合在线观看 | 国产精品99久久久久久久久| 亚洲精品国产成人久久av| 国产色婷婷99| 各种免费的搞黄视频| av在线老鸭窝| 黑人高潮一二区| 一级爰片在线观看| 国产亚洲欧美精品永久| 日韩欧美一区视频在线观看 | 日韩电影二区| 男女国产视频网站| 尾随美女入室| 最近中文字幕高清免费大全6| 99国产精品免费福利视频| 久久青草综合色| 亚洲中文av在线| 美女内射精品一级片tv| 插逼视频在线观看| 国产精品熟女久久久久浪| 国产免费福利视频在线观看| 在现免费观看毛片| 亚洲av欧美aⅴ国产| 在线免费观看不下载黄p国产| 免费人妻精品一区二区三区视频| 国产69精品久久久久777片| 少妇丰满av| 又爽又黄a免费视频| 国产精品人妻久久久影院| h视频一区二区三区| 国产高清国产精品国产三级| 中国美白少妇内射xxxbb| 视频中文字幕在线观看| 欧美日本中文国产一区发布| av国产精品久久久久影院| av有码第一页| 亚洲高清免费不卡视频| 欧美日韩综合久久久久久| 插阴视频在线观看视频| 欧美三级亚洲精品| 精品一区二区三区视频在线| 亚洲av欧美aⅴ国产| 久久精品夜色国产| 一区二区三区精品91| 国产精品久久久久久精品古装|