• <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)測量精度的評估
    我要飛翔
    少妇高潮的动态图| 少妇人妻精品综合一区二区 | 亚洲最大成人中文| 欧美日本视频| 国产毛片a区久久久久| 波野结衣二区三区在线| 国产一区二区三区视频了| 日韩中字成人| 亚洲七黄色美女视频| 在线看三级毛片| 色视频www国产| 又爽又黄无遮挡网站| 中出人妻视频一区二区| 日韩精品中文字幕看吧| avwww免费| 国产精品爽爽va在线观看网站| 欧美色视频一区免费| 九九爱精品视频在线观看| 久久精品国产亚洲av涩爱 | 亚洲人成网站在线播放欧美日韩| 高清日韩中文字幕在线| 嫁个100分男人电影在线观看| 蜜桃久久精品国产亚洲av| 伊人久久精品亚洲午夜| 国产午夜精品久久久久久一区二区三区 | 国产探花极品一区二区| 日本与韩国留学比较| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 最近最新中文字幕大全电影3| 欧美不卡视频在线免费观看| 99久久精品热视频| 久久久久国内视频| 观看免费一级毛片| 精品一区二区三区人妻视频| 最新中文字幕久久久久| 日韩欧美国产一区二区入口| 精品久久久久久成人av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲天堂国产精品一区在线| 久久久精品大字幕| 免费在线观看日本一区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲电影在线观看av| 在线观看av片永久免费下载| 国产男靠女视频免费网站| 国产蜜桃级精品一区二区三区| 长腿黑丝高跟| 久久亚洲真实| 一夜夜www| 日本-黄色视频高清免费观看| 在线国产一区二区在线| 日韩欧美一区二区三区在线观看| 亚洲精品亚洲一区二区| 国产精品免费一区二区三区在线| 免费电影在线观看免费观看| 国产三级在线视频| av国产免费在线观看| 欧美+亚洲+日韩+国产| 村上凉子中文字幕在线| 久久99热这里只有精品18| 久久人人爽人人爽人人片va| 一个人免费在线观看电影| 欧美zozozo另类| 成人鲁丝片一二三区免费| 精品久久久久久久人妻蜜臀av| 人妻夜夜爽99麻豆av| 长腿黑丝高跟| 国产亚洲欧美98| 久久久久免费精品人妻一区二区| 欧美日韩精品成人综合77777| 日本免费a在线| 一级av片app| 伊人久久精品亚洲午夜| 国产私拍福利视频在线观看| 男人狂女人下面高潮的视频| 国产伦一二天堂av在线观看| 草草在线视频免费看| 一个人观看的视频www高清免费观看| 精品日产1卡2卡| 成人高潮视频无遮挡免费网站| 欧美日韩精品成人综合77777| 亚洲欧美日韩卡通动漫| 久久久久久久久久久丰满 | 亚洲成人精品中文字幕电影| 亚洲黑人精品在线| 观看免费一级毛片| 97碰自拍视频| 99国产精品一区二区蜜桃av| 国产精品永久免费网站| 免费观看在线日韩| 91久久精品国产一区二区三区| 51国产日韩欧美| 日韩欧美在线乱码| 亚洲成人精品中文字幕电影| 成人无遮挡网站| 乱码一卡2卡4卡精品| 九九在线视频观看精品| 国产成人福利小说| 亚洲一区高清亚洲精品| 久久久久国产精品人妻aⅴ院| 国产不卡一卡二| 国产精品久久久久久av不卡| 亚洲在线观看片| 国产高潮美女av| 久久久久久国产a免费观看| 久久久久性生活片| 不卡视频在线观看欧美| 成人国产麻豆网| 久久午夜福利片| 久久人妻av系列| 日本黄大片高清| 十八禁网站免费在线| 国产黄色小视频在线观看| 亚洲av日韩精品久久久久久密| 中文字幕高清在线视频| 国产乱人伦免费视频| 免费观看在线日韩| 久久国产乱子免费精品| 国产亚洲91精品色在线| 日本免费一区二区三区高清不卡| 老熟妇仑乱视频hdxx| 国产精品三级大全| 午夜福利视频1000在线观看| 老师上课跳d突然被开到最大视频| 少妇高潮的动态图| 亚洲国产精品sss在线观看| 美女黄网站色视频| 一区二区三区免费毛片| 国产单亲对白刺激| 国产亚洲精品av在线| 日本免费一区二区三区高清不卡| 国产精品免费一区二区三区在线| 久久久精品欧美日韩精品| 少妇猛男粗大的猛烈进出视频 | 在现免费观看毛片| 国产成人影院久久av| 嫁个100分男人电影在线观看| 嫩草影院入口| a级毛片免费高清观看在线播放| 亚洲av美国av| 一进一出抽搐gif免费好疼| 动漫黄色视频在线观看| 久久精品国产自在天天线| 深爱激情五月婷婷| 能在线免费观看的黄片| 亚洲图色成人| 中文字幕高清在线视频| 久久精品91蜜桃| 国产不卡一卡二| 欧美性感艳星| 国产人妻一区二区三区在| 午夜久久久久精精品| 成人特级黄色片久久久久久久| 哪里可以看免费的av片| 97超级碰碰碰精品色视频在线观看| 一个人看视频在线观看www免费| 国内揄拍国产精品人妻在线| 国产麻豆成人av免费视频| 内射极品少妇av片p| 黄色女人牲交| 99热6这里只有精品| 久久精品国产亚洲av涩爱 | 亚洲自拍偷在线| 性插视频无遮挡在线免费观看| 十八禁国产超污无遮挡网站| 久久精品人妻少妇| 一区二区三区四区激情视频 | 成人亚洲精品av一区二区| 欧美不卡视频在线免费观看| 国产男人的电影天堂91| 婷婷丁香在线五月| 久久久久久九九精品二区国产| 制服丝袜大香蕉在线| 久久久久精品国产欧美久久久| 国产精品综合久久久久久久免费| 欧美激情在线99| 精品久久久久久久久久免费视频| 97人妻精品一区二区三区麻豆| 午夜福利在线观看吧| 国产伦精品一区二区三区四那| 18禁黄网站禁片免费观看直播| 99riav亚洲国产免费| 麻豆一二三区av精品| 18禁黄网站禁片免费观看直播| 春色校园在线视频观看| 日本熟妇午夜| 久久精品综合一区二区三区| 国产精品永久免费网站| 2021天堂中文幕一二区在线观| 嫩草影院入口| 国产亚洲精品综合一区在线观看| 深夜精品福利| 亚洲人成网站在线播| 亚洲中文字幕一区二区三区有码在线看| 舔av片在线| 久久久久久久久久久丰满 | 99热精品在线国产| 日韩中字成人| 亚洲人与动物交配视频| 亚洲av二区三区四区| 国产一区二区激情短视频| 国产精品久久久久久久久免| 22中文网久久字幕| 久久久久国内视频| 亚洲精华国产精华液的使用体验 | 精品人妻1区二区| 99国产精品一区二区蜜桃av| 色综合婷婷激情| 一区福利在线观看| 久99久视频精品免费| 国产毛片a区久久久久| 床上黄色一级片| 亚洲欧美精品综合久久99| 国产一区二区三区视频了| 舔av片在线| 如何舔出高潮| 成人国产综合亚洲| 高清日韩中文字幕在线| 久9热在线精品视频| 亚洲av中文字字幕乱码综合| 国产女主播在线喷水免费视频网站 | 久久久国产成人免费| 国产色婷婷99| 国产女主播在线喷水免费视频网站 | 看片在线看免费视频| 成人av一区二区三区在线看| 国产欧美日韩精品亚洲av| 亚洲欧美日韩高清专用| 人妻丰满熟妇av一区二区三区| 两个人的视频大全免费| av在线观看视频网站免费| 熟妇人妻久久中文字幕3abv| 免费观看在线日韩| 男女下面进入的视频免费午夜| 女人十人毛片免费观看3o分钟| 老司机福利观看| 精品福利观看| 色综合亚洲欧美另类图片| 亚洲欧美日韩卡通动漫| 国产精品久久电影中文字幕| .国产精品久久| 日日撸夜夜添| 午夜老司机福利剧场| 国产aⅴ精品一区二区三区波| 在线观看免费视频日本深夜| 97超视频在线观看视频| 国产欧美日韩一区二区精品| 久久久久久久久中文| 久久欧美精品欧美久久欧美| 免费人成视频x8x8入口观看| 欧美激情国产日韩精品一区| 亚洲熟妇熟女久久| 免费黄网站久久成人精品| 男女啪啪激烈高潮av片| 久久午夜亚洲精品久久| 亚洲va日本ⅴa欧美va伊人久久| 变态另类丝袜制服| 亚洲在线自拍视频| 日本 av在线| 久久精品国产亚洲网站| 悠悠久久av| 亚洲精品影视一区二区三区av| 午夜免费男女啪啪视频观看 | 国产成年人精品一区二区| 国产v大片淫在线免费观看| 22中文网久久字幕| 国产欧美日韩一区二区精品| 夜夜看夜夜爽夜夜摸| 搡老妇女老女人老熟妇| a级一级毛片免费在线观看| 亚洲综合色惰| 成人亚洲精品av一区二区| 久久欧美精品欧美久久欧美| 久99久视频精品免费| 在线观看免费视频日本深夜| 亚洲五月天丁香| 尾随美女入室| 淫秽高清视频在线观看| 99久久精品一区二区三区| 亚洲综合色惰| 欧美另类亚洲清纯唯美| 亚洲自拍偷在线| 乱系列少妇在线播放| 久久久久久久久久久丰满 | 美女cb高潮喷水在线观看| 内射极品少妇av片p| 久久精品综合一区二区三区| 可以在线观看毛片的网站| 高清日韩中文字幕在线| 99久国产av精品| 国产精品久久久久久亚洲av鲁大| 一进一出抽搐动态| 国产高潮美女av| 18禁裸乳无遮挡免费网站照片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av日韩精品久久久久久密| 十八禁国产超污无遮挡网站| 国产91精品成人一区二区三区| 欧美成人免费av一区二区三区| av天堂在线播放| 别揉我奶头~嗯~啊~动态视频| 国产视频一区二区在线看| 精品免费久久久久久久清纯| 亚洲性久久影院| 久久久久久国产a免费观看| 日本一本二区三区精品| 香蕉av资源在线| 成熟少妇高潮喷水视频| 欧美日韩精品成人综合77777| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区性色av| 午夜福利在线观看吧| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区三区在线臀色熟女| 在线国产一区二区在线| 琪琪午夜伦伦电影理论片6080| 十八禁国产超污无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 级片在线观看| 久久久成人免费电影| 免费看日本二区| 日日啪夜夜撸| 99热这里只有是精品在线观看| 日韩人妻高清精品专区| 久久久色成人| 999久久久精品免费观看国产| 中国美女看黄片| 婷婷亚洲欧美| 欧美黑人巨大hd| 不卡视频在线观看欧美| а√天堂www在线а√下载| 日韩欧美在线乱码| 亚洲性夜色夜夜综合| 赤兔流量卡办理| 波多野结衣高清作品| 最近最新中文字幕大全电影3| 一区福利在线观看| 别揉我奶头 嗯啊视频| 欧美3d第一页| 神马国产精品三级电影在线观看| 别揉我奶头~嗯~啊~动态视频| 国产午夜福利久久久久久| 午夜老司机福利剧场| 婷婷精品国产亚洲av| 久久6这里有精品| 搞女人的毛片| 成人午夜高清在线视频| 成人永久免费在线观看视频| 一个人免费在线观看电影| 看黄色毛片网站| 精品久久久久久久久av| 窝窝影院91人妻| 久久精品影院6| 午夜精品一区二区三区免费看| 91麻豆精品激情在线观看国产| 直男gayav资源| 最新中文字幕久久久久| 久久精品国产99精品国产亚洲性色| 99国产极品粉嫩在线观看| 狂野欧美白嫩少妇大欣赏| 一进一出抽搐gif免费好疼| 久9热在线精品视频| 身体一侧抽搐| 成人性生交大片免费视频hd| 午夜福利高清视频| 18禁黄网站禁片午夜丰满| 两个人的视频大全免费| 亚洲人与动物交配视频| av天堂在线播放| 亚洲欧美激情综合另类| 国产一区二区三区av在线 | 熟女人妻精品中文字幕| ponron亚洲| 久9热在线精品视频| 人妻夜夜爽99麻豆av| 最近在线观看免费完整版| 日本黄大片高清| 男人和女人高潮做爰伦理| 欧美三级亚洲精品| 简卡轻食公司| 嫩草影视91久久| 99热6这里只有精品| 欧美又色又爽又黄视频| 女人被狂操c到高潮| 国产av在哪里看| xxxwww97欧美| 国产 一区精品| 欧美中文日本在线观看视频| 精品久久久久久,| 在现免费观看毛片| 日本熟妇午夜| 少妇的逼好多水| 熟女人妻精品中文字幕| 欧美精品啪啪一区二区三区| 国产精品久久久久久精品电影| 国产真实乱freesex| 波野结衣二区三区在线| 狠狠狠狠99中文字幕| 窝窝影院91人妻| 黄色女人牲交| 又爽又黄无遮挡网站| av在线蜜桃| 免费搜索国产男女视频| 中文字幕高清在线视频| 国产精品免费一区二区三区在线| 亚洲精品456在线播放app | 真人一进一出gif抽搐免费| 成人高潮视频无遮挡免费网站| 国产精品三级大全| 又爽又黄a免费视频| 国产午夜精品久久久久久一区二区三区 | 91久久精品电影网| 在线看三级毛片| 国产真实伦视频高清在线观看 | 亚洲综合色惰| 中文字幕久久专区| 亚洲av日韩精品久久久久久密| 在线播放无遮挡| 午夜亚洲福利在线播放| 亚洲精品一卡2卡三卡4卡5卡| 久久久久国内视频| 两个人的视频大全免费| 麻豆精品久久久久久蜜桃| 日本与韩国留学比较| 极品教师在线视频| 精品一区二区免费观看| 亚洲三级黄色毛片| 美女黄网站色视频| 高清在线国产一区| 久久人人精品亚洲av| 亚洲成人久久爱视频| 国内少妇人妻偷人精品xxx网站| 日韩 亚洲 欧美在线| av中文乱码字幕在线| 如何舔出高潮| 香蕉av资源在线| 国产色婷婷99| 听说在线观看完整版免费高清| 亚洲成av人片在线播放无| 亚洲av一区综合| 亚洲av五月六月丁香网| www.www免费av| 免费在线观看影片大全网站| 国产视频内射| 搡女人真爽免费视频火全软件 | 成人精品一区二区免费| 国内毛片毛片毛片毛片毛片| 老熟妇乱子伦视频在线观看| 国产精品久久视频播放| 日韩欧美 国产精品| 国产日本99.免费观看| 精品久久久久久成人av| 国内少妇人妻偷人精品xxx网站| 51国产日韩欧美| 深夜a级毛片| 别揉我奶头 嗯啊视频| 亚洲自偷自拍三级| 亚洲av中文字字幕乱码综合| 男女之事视频高清在线观看| 99视频精品全部免费 在线| 男人和女人高潮做爰伦理| 亚洲精品亚洲一区二区| 色av中文字幕| 国产精品,欧美在线| 日日摸夜夜添夜夜添小说| 日韩欧美免费精品| 国产老妇女一区| 欧美色视频一区免费| 97碰自拍视频| 精品久久久久久成人av| 欧美精品国产亚洲| 国产精品亚洲美女久久久| 欧美黑人巨大hd| 三级国产精品欧美在线观看| 99热这里只有是精品在线观看| avwww免费| 亚洲在线自拍视频| 美女xxoo啪啪120秒动态图| 老熟妇乱子伦视频在线观看| 1024手机看黄色片| 午夜影院日韩av| 色尼玛亚洲综合影院| 五月玫瑰六月丁香| 亚洲性久久影院| 久久午夜福利片| 日韩欧美 国产精品| 又爽又黄无遮挡网站| www.色视频.com| 欧美潮喷喷水| 亚洲七黄色美女视频| 一区二区三区激情视频| 国内精品久久久久久久电影| avwww免费| 亚洲在线观看片| 深爱激情五月婷婷| 国产视频一区二区在线看| 18禁黄网站禁片午夜丰满| 亚洲精品成人久久久久久| 国产大屁股一区二区在线视频| 亚洲精品粉嫩美女一区| 人妻制服诱惑在线中文字幕| 禁无遮挡网站| 嫩草影视91久久| 亚洲人成网站在线播| 久久久久久久久大av| 少妇的逼水好多| 少妇猛男粗大的猛烈进出视频 | 三级男女做爰猛烈吃奶摸视频| 国产人妻一区二区三区在| 美女免费视频网站| 黄色欧美视频在线观看| 观看美女的网站| 夜夜看夜夜爽夜夜摸| 久久久久性生活片| 日日撸夜夜添| 国产精品不卡视频一区二区| 国产精品久久视频播放| 日韩欧美国产一区二区入口| 久久久久久久久久久丰满 | 国产高清三级在线| 偷拍熟女少妇极品色| 亚洲专区国产一区二区| 国产成人av教育| 国产高清激情床上av| 亚洲av一区综合| 欧美最新免费一区二区三区| 国产美女午夜福利| 婷婷六月久久综合丁香| 日日摸夜夜添夜夜添小说| 91久久精品国产一区二区成人| 国产女主播在线喷水免费视频网站 | 国产精品久久电影中文字幕| 麻豆国产97在线/欧美| 又黄又爽又免费观看的视频| 国产欧美日韩一区二区精品| 黄片wwwwww| 欧美潮喷喷水| 毛片一级片免费看久久久久 | 美女cb高潮喷水在线观看| 自拍偷自拍亚洲精品老妇| 日韩欧美国产在线观看| 听说在线观看完整版免费高清| 小说图片视频综合网站| 亚洲成a人片在线一区二区| 少妇丰满av| 国产欧美日韩精品一区二区| 免费搜索国产男女视频| 嫩草影视91久久| 国产av一区在线观看免费| 久久久精品欧美日韩精品| 九九在线视频观看精品| 干丝袜人妻中文字幕| 精品乱码久久久久久99久播| 欧美绝顶高潮抽搐喷水| 麻豆国产av国片精品| 少妇被粗大猛烈的视频| 成人一区二区视频在线观看| 一个人观看的视频www高清免费观看| 亚洲专区中文字幕在线| 久久久久久久久久久丰满 | 99精品在免费线老司机午夜| 久久6这里有精品| 久久久久久九九精品二区国产| 免费黄网站久久成人精品| 丰满乱子伦码专区| 99视频精品全部免费 在线| 国产69精品久久久久777片| 日本精品一区二区三区蜜桃| 欧美精品啪啪一区二区三区| 熟女人妻精品中文字幕| 久久久成人免费电影| 搡老熟女国产l中国老女人| 国产亚洲精品久久久com| 亚洲精品成人久久久久久| 国产精品野战在线观看| 国模一区二区三区四区视频| 村上凉子中文字幕在线| 免费在线观看日本一区| 日本 欧美在线| 日韩精品青青久久久久久| 又爽又黄a免费视频| 免费观看人在逋| 12—13女人毛片做爰片一| 亚洲人成网站在线播| 亚洲最大成人av| 亚洲第一区二区三区不卡| 欧美极品一区二区三区四区| 男女边吃奶边做爰视频| 国产老妇女一区| 久久草成人影院| 久久久久久久午夜电影| 国语自产精品视频在线第100页| 一本一本综合久久| 成人国产一区最新在线观看| 欧美丝袜亚洲另类 | 欧美xxxx性猛交bbbb| 日日啪夜夜撸| 老司机福利观看| 他把我摸到了高潮在线观看| 亚洲在线自拍视频| 美女高潮喷水抽搐中文字幕| 国产av麻豆久久久久久久| 国产精品女同一区二区软件 | 69av精品久久久久久| 夜夜看夜夜爽夜夜摸| 国产成人av教育| 精品久久久久久久久亚洲 | 免费不卡的大黄色大毛片视频在线观看 | 亚洲无线在线观看| 日本熟妇午夜| 午夜免费激情av| 久久久午夜欧美精品| 蜜桃久久精品国产亚洲av| 久久精品影院6|