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

    基于聚類分析的華鎣山褶皺山系地下水水化學(xué)分類研究

    2013-12-14 02:52:36張明亮彭曉鳳
    地下水 2013年1期
    關(guān)鍵詞:華鎣山水化學(xué)水文地質(zhì)

    張明亮,許 模,彭曉鳳

    (1.成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國家重點實驗室,四川成都610059;2.四川省地質(zhì)礦產(chǎn)勘查開發(fā)局九一五水文地質(zhì)工程地質(zhì)隊,四川眉山620010)

    賦存于巖石圈中的地下水在運(yùn)移過程中不斷與圍巖發(fā)生各種化學(xué)反應(yīng),從而導(dǎo)致化學(xué)元素的遷移、聚集和分散[1]。由于巖性、巖石結(jié)構(gòu)、地下水徑流條件等的不同,地下水水化學(xué)組分差異較大。因此研究地下水的水化學(xué)組分,對掌握水文地質(zhì)條件有及其重要的作用。地下水系統(tǒng)的水化學(xué)分類及判別研究是水文地質(zhì)學(xué)的重要內(nèi)容之一[2]。為了更準(zhǔn)確和細(xì)致地劃分自然界各種水體的化學(xué)類型,前人做了許多相關(guān)的研究。美國科學(xué)家Hill[3]1940年首次使用三線圖對水化學(xué)類型進(jìn)行了描述。Piper[4]1944年對前人的三線圖進(jìn)行了修 改。隨 后,Durov[5]又 于 1948 年 推 出 了 Durov圖。1999Chadha[6]在充分總結(jié)前人經(jīng)驗的基礎(chǔ)上又推出了矩形圖。Piper圖和Durov圖能直觀展示區(qū)域水化學(xué)類型,但分類不精細(xì),特別是當(dāng)有大量的水化學(xué)數(shù)據(jù)時,對于水化學(xué)的定名和水化學(xué)分類不是很方便。在矩形水化學(xué)圖中SO42-與Cl-不分,Ca2+與Mg2+不分,而且當(dāng)各種離子相對質(zhì)量濃度為25%~75%時,次一級矩形中無法明確劃分水化學(xué)類型[7]。這幾種圖解在現(xiàn)今水化學(xué)的研究中都被廣泛應(yīng)用,但是,隨著學(xué)科的發(fā)展,傳統(tǒng)分類方法在科研和生產(chǎn)中的不適應(yīng)的現(xiàn)象越來越突出地表現(xiàn)了出來,尤其是在處理水污染問題的時候常常會有重金屬或者微量元素,甚至是有機(jī)物的混入[8]。而且由于分類標(biāo)準(zhǔn)大都是人為給定,例如舒氏分類中規(guī)定毫克當(dāng)量大于25%者才參與分類,從而會造成有時將兩個分析結(jié)果差別很小的水樣被列入不同類型,特別是對水化學(xué)成份差別不大的地區(qū),上述方法尚欠精細(xì)[9]。于是,數(shù)學(xué)工具就被逐漸引入了水化學(xué)分類研究,隨著后來多元分析的加入,便有了今天的水化學(xué)聚類分析分類方法。該方法排除了人為干擾,聚類精確,圖形顯示直觀,特別是水樣量大、組分極多的情況效果更佳;而且能根據(jù)研究的內(nèi)容調(diào)整并類距離,取得最好的解釋效果。

    1 聚類分析簡介

    聚類分析是一種考慮了多因素的數(shù)學(xué)分類方法,其基本原理是根據(jù)樣本自身的屬性,用數(shù)學(xué)的方法按照某些相似性指標(biāo),定量地確定樣本之間的親疏關(guān)系,并按這種親疏關(guān)系程度對樣本進(jìn)行聚類,與分類不同的是——聚類要求劃分的類是未知的[10]。

    系統(tǒng)聚類法是目前應(yīng)用最廣泛的一種聚類方法。其基本思想是:1、每個樣品為一類,計算各樣品之間的距離系數(shù);2、把距離最小的兩類合并為一類;3、計算新的類間的距離;4、重復(fù)2、3步一直到合并為一類為止[11]。這樣一種連續(xù)并類的過程可用一種類似于樹狀結(jié)構(gòu)的圖形即聚類譜系圖(俗稱樹狀圖)來表示,由聚類譜系圖可以清楚地看出全部樣本的聚類過程,從而作出對全部樣本的分類。系統(tǒng)聚類分為Q型聚類和R型聚類兩種:Q型聚類是對樣本進(jìn)行聚類,R型聚類是對變量進(jìn)行聚類。本次研究對于所取水樣的聚類采用Q型聚類分析。

    2 華鎣山研究區(qū)基本地質(zhì)概況

    川東隔檔式構(gòu)造位于四川盆地的東部,由一系列NENNE向的隔檔式褶皺組成,分布面積約5.3萬 km2。其西以華鎣山斷裂為界,與褶皺很不發(fā)育的川中平緩褶皺帶相鄰,其東與滇黔川鄂臺褶帶相隔于齊躍山斷裂,單個背斜呈線形延伸,左行斜列式展開,總體呈圍繞盆地東南邊界的弧形帚狀組合[12]。

    華鎣山背斜為一箱狀褶皺,處于隔檔式背斜的最西部,該復(fù)式背斜向南分岔在觀音峽、瀝鼻峽、溫塘峽三條隆起帶和北碚、壁山兩條沉降帶,形成川東褶皺的帚狀構(gòu)造。華鎣山背斜軸部開闊平緩,兩翼不對稱,北西翼陡,傾角60°~80°,南東翼緩,傾角 20°~40°,軸向 N10°~25°E,軸部向北傾伏,傾伏角2°~7°,見圖2。西槽谷發(fā)育,貫穿整個華鎣山背斜,東邊槽谷僅隧址區(qū)南段較發(fā)育,北段零星發(fā)育,見圖1研究區(qū)地貌圖。南段出露最老地層為三疊系嘉陵江組,北段出露的最老地層為三疊系雷口坡組,兩翼由須家河組、自流井組、新田溝組、沙溪廟組組成,研究區(qū)水文地質(zhì)圖見圖3。研究區(qū)范圍主要為州河向南擴(kuò)展約40 km的華鎣山褶皺山系,為華鎣山北段。研究區(qū)北側(cè)核部出露最老地層為雷口坡組(T2l)地層,核部可溶巖地層出露較窄。向南在張家溝-三叉河處開始出露嘉陵江組(T1j)地層,核部可溶巖地層變寬。受巖性控制,318國道南側(cè)巖溶地貌相對北側(cè)發(fā)育,順軸線發(fā)育連串大型洼地,溶溝、漏斗、落水洞、溶蝕槽谷、溶丘等地貌,北側(cè)巖溶地貌發(fā)育較少,一般為小型的洼地、溶溝,偶見落水洞和漏斗。同時受構(gòu)造控制,背斜核部西側(cè)可溶巖地層中沿軸向連串的巖溶洼地和溶溝形成一個連通的巖溶槽谷地貌,向北延伸至州河,高程246~605 m,州河為研究區(qū)排泄基準(zhǔn)面。核部東側(cè)在318國道南側(cè),九盤寺-梨竹寺-館子壩-小壩一帶沿軸向連串的巖溶洼地、漏斗、落水洞形成巖溶槽谷地形,但向北側(cè)延伸至核桃坪,洼地、漏斗逐漸減少,隧址區(qū)的中部和北側(cè)主要為巖溶山地地貌。東槽谷只在隧址區(qū)的南側(cè)發(fā)育,高程780~880 m,明顯西槽谷的切割地形遠(yuǎn)低于東槽谷的切割。因此,華鎣山隧道研究區(qū)南側(cè)為“一山三嶺二槽”地貌,北側(cè)為“一山二嶺一槽”地貌。隔檔式背斜獨(dú)特的地下水循環(huán)模式使得在背斜核部常形成高水位的槽谷,隧道穿越一旦揭露巖溶管道或大型溶洞,往往發(fā)生較大的涌突水和突泥災(zāi)害。所以對水化學(xué)組分的研究尤其重要,本文用Q型聚類方法對水樣進(jìn)行分析,以求結(jié)合分類更清楚的掌握研究區(qū)的水文地質(zhì)條件。

    圖1 華鎣山北段地貌圖

    圖2 華鎣山北段剖面示意圖

    圖3 華鎣山北段水文地質(zhì)簡圖

    3 水化學(xué)組分聚類分析

    取2010年8月在華鎣山背斜地區(qū)采集22個水樣水化學(xué)分析結(jié)果為研究對象,22組水樣中有8組水樣出露于T1j中,13組水樣出露于雷口坡T2l中,1組出露于須家河(T3xj)中。因地下水中8種常量離子總量可占水中TDS的95%~99%[13]。這些離子屬地下水中的標(biāo)型組分,其差異是地下水化學(xué)分類的主要依據(jù)。所以選取 Na+、K+、Mg2+、Ca2+、Cl-、、七大主離子作為特征變量,因為這七大離子是水的主要成分,水是什么樣的水化學(xué)性質(zhì)主要取決于這些離子的含量及組合情況,所以去這些離子作為特征變量是有理可依的,這些特征變量足以決定水類別。

    在進(jìn)行模糊聚類分析水化學(xué)類型時,按共原理設(shè)n個樣本(X1、X2、X3…Xn),因各個樣本有 m個分析指標(biāo),故任何一個樣本都是一個m維向量,可記為T,由此可構(gòu)成一個n×m階矩陣。

    1)數(shù)據(jù)標(biāo)準(zhǔn)化。由于各個變量的量綱不一致,即使統(tǒng)一量綱有事原始數(shù)據(jù)的大小也有懸殊,為了避免有些特征變量受到壓抑,在分類前先對數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化。

    2)計算類與類之間的距離,這里選取標(biāo)準(zhǔn)歐氏距離計算各組水樣的距離。

    3)找出最小值,然后選出次小值進(jìn)行連接。

    4)評價聚類結(jié)果。通過對不同距離計算方法及聚類方法的搭配,計算不同組合的相關(guān)系數(shù),結(jié)合實際分析,得到采用標(biāo)準(zhǔn)歐氏距離和離差平方和的聚類結(jié)果較符合實際的需要。聚類分析的譜系圖見圖4。

    在分類取多大水平值為宜,主要視具體地區(qū)的地貌,水文地質(zhì)條件,地下水的形成和運(yùn)動條件而定,以能反應(yīng)水化學(xué)成因規(guī)律為原則[14]。由聚類分析樹狀圖可以看出采用聚類分析對水樣進(jìn)行分類效果是非常明顯的。

    當(dāng)并類距離取d取25時,22組水樣分為兩類,即①②③④⑤為一類,⑥為一類,顯然這樣的分類不能表現(xiàn)出各組水樣的內(nèi)在聯(lián)系,不能反映該地區(qū)的水文地質(zhì)特征,分類過于粗糙,僅能反映22號水樣與其它水樣差異太大,可能受到污染。

    當(dāng)并類距離d=13時,22組水樣分為3類,即①②③④為一類,⑤為一類,⑥為一類。同樣,①②③④為一類,這一類中有19組水樣,這樣的分類不能表征它們的相互關(guān)系。⑤類中的礦化度較高,有一定的解釋意義。⑥類為22號水樣,能反映與其它水樣的差別。

    圖4 水化學(xué)組分聚類樹狀圖

    若將并類距離d設(shè)為10,則可將22組水樣分為4類,即①②③為一類,④、⑤、⑥分別自成一類。①②③類水樣都出露于雷口坡(T2l)和嘉陵江組(T1j)灰?guī)r地層當(dāng)中,具有巖溶水的性質(zhì),大都是 HCO3-Ca、HCO3-Ca·Mg型水,屬于巖溶水系統(tǒng)。④類由1、2、3、7、9號水樣組成,該組水樣的特點是徑流途徑短,礦化度較低。1、2號泉點出露于盧家河壩T2l地層中,該處T2l為T3xj地層中揭露的構(gòu)造天窗,水化學(xué)類型為HCO3·SO4-Ca·Mg型水,由于須家河地層中含硫,1、2號泉點在徑流過程中應(yīng)與須家河地層中的水發(fā)生混合,導(dǎo)致含量增高;7、9號水樣與1、2號水化學(xué)組分相似,僅含量低,因為其徑流途徑上穿越的須家河地層少,可以判斷7、9號水樣與1、2號補(bǔ)給源相同,反映出華鎣山背斜中的巖溶水受溝谷切割向兩側(cè)溝谷排泄的特征。⑤類由15和17號水樣組成,這兩個水樣都是暗河水,其特點是徑流途徑長,礦化度較高。如15號水樣的就達(dá) 330 ml/L,在22組水樣中最高。17號水樣取水點位于硐新店子(船石村)五組的一個地下暗河的出口,出露于嘉陵江(T1j)地層中,該地層主要的巖性為白云巖,按常理推斷水化學(xué)類型應(yīng)該為HCO3-Ca型或HCO3-Ca·Mg型水,而此處采樣的水卻是HCO3·SO4-Ca型水,而且的含量高達(dá)140 mg/L,該地層又無硫離子的來源,可以斷定該暗河的水遭到了附近工廠的污染。⑥類由22號單獨(dú)組成,該點的水為華鎣山核部匯入州河的支溝中的水,該處水樣的礦化度比其他水點的都高達(dá)744 mg/L,僅 SO42-的含量就達(dá)410 mg/L,可能是由于沿途接受有煤礦排水的緣故,可以斷定該組水樣已受污染。所以并類距離取為10較為合適。

    如將并類距離設(shè)為9,則可將水樣分為5類,即①②為一類,③、④、⑤、⑥分別自成一類。這樣在劇烈距離為10的基礎(chǔ)上再將①②③中的14組樣品細(xì)分兩類,即將①②和③分開??梢钥闯觫垲愔械?0、20、21號水樣都是出露于嘉陵江組(T1j)之中,該類水樣均屬于HCO3-Ca型水,其各組分指標(biāo)都較為接近,再根據(jù)其出露地點判斷,它們具有相同的補(bǔ)給徑流途徑,屬于同一個系統(tǒng)。①②類則全部屬于巖溶水。所以并類距離為9較10更好。

    如將并類距離再減小至8,則水樣分為6類,則在并類距離為9的基礎(chǔ)上把⑤中的15和17號分開,這樣的分類意義不大。并類距離減小到7時分為8類,分類的意義驟減,已不能反應(yīng)水樣的整體規(guī)律。若距離再減小,分類失去了意義。但是若從局部看,小距離的聚類還是有一定意義的,如當(dāng)并類距離d=1時,可以看出12、13為一類,結(jié)合實際分析他們?yōu)橥谎a(bǔ)給源。d=2時1、2可以為一類;它們的補(bǔ)給徑流都相同,只是排泄有所差異;d=5時①類(12、13、4、14、16)均為HCO3-Ca型水。所以應(yīng)根據(jù)研究的內(nèi)容不同選擇合適的并類距離,進(jìn)行分類的解釋。

    4 結(jié)語

    通過聚類分析,將在華鎣山地區(qū)取的22組水樣進(jìn)行了分類。取并類距離d=9將22組水樣分為5類,取得了較好的解釋結(jié)果。反映了該地區(qū)的巖溶水的水化學(xué)特征,結(jié)合分類能更清楚地掌握該地區(qū)同類泉點的補(bǔ)徑排特征,第三節(jié)中已詳述,對研究該地區(qū)的水文地質(zhì)條件有很重要的指導(dǎo)意義。

    應(yīng)根據(jù)研究需要選擇合適的并類距離,距離太大則聚類太粗糙,距離太小則失去聚類的意義。但是若從局部看,小距離的聚類還是有一定意義的,不容忽視。

    可以將未知水樣列入分類矩陣再進(jìn)行聚類,從而根據(jù)已知水樣的補(bǔ)徑排特征判斷未知水樣的補(bǔ)徑排特征。

    采用模糊數(shù)學(xué)聚類分析法進(jìn)行局部地區(qū)的地下水水化學(xué)分類,可以在樣本相近,具有模糊性的情況下進(jìn)行。能解決在單一地貌條件、單一水文地質(zhì)單元及小范圍內(nèi)的水化學(xué)分類。

    克服了Piper圖和Durov圖能直觀展示區(qū)域水化學(xué)類型,但分類不精細(xì),特別是當(dāng)有大量的水化學(xué)數(shù)據(jù)時,對于水化學(xué)的定名和水化學(xué)分類不是很方便的缺點。排除了舒氏分類人為的劃定以25%毫克當(dāng)量作為劃分水型的依據(jù)的人為性??梢匝杆俸Y選處異常數(shù)據(jù),特別是對地下水出露地層與其水化學(xué)性質(zhì)不符的水樣,能進(jìn)一步提高解決問題的速度和精度。

    聚類分析可以將分析結(jié)果直觀的用聚類樹狀圖展現(xiàn)出來,而且效率高,性能準(zhǔn)確可靠??梢愿鶕?jù)不同的研究目的來確定并類距離,合理解釋不同分類水樣的水文地質(zhì)意義。

    [1]王大純,張人權(quán),史毅虹,等.水文地質(zhì)學(xué)基礎(chǔ).北京:地質(zhì)出版社.1995.

    [2]江曉益,成春奇.礦區(qū)地下水系統(tǒng)水質(zhì)分類判別的多元統(tǒng)計分析[J].水文地質(zhì)工程地質(zhì).2009,(4):16-20.

    [3]Hill R A.1940.Geochemical Patterns in the Coachella Valley,California[J].Trans Am Geophys Union,21:46~ 49.

    [4]Piper A M.1944.A Graphic Procedure in GeochemicalⅠnterpretation of Water Analyses[J].Trans Am Geophys Union,25:914~ 923.

    [5]Durov S A.1948.Natural Waters and Graphic Representation of Their Compositions[J].Dokl Akad Nauk Ssssr,59:87—90.

    [6]Chadha D K.1999.A Proposed New Diagram for Geochemical Classification of Natural Waters and Ⅰnterpretation of Chemical Data[J].Hydrogeology Journal,7:431~ 439.

    [7]孫亞喬,錢會等.基于矩形圖的天然水化學(xué)分類和水化學(xué)規(guī)律研究[J].地球科學(xué)與環(huán)境學(xué)報.2007,29(1):75-79.

    [8]李慶生,伍坤宇.聚類分析在水化學(xué)分類研究中的應(yīng)用[J].化學(xué)工程與裝備.2010(11):107-109.

    [9]王文科,吳在寶.FuzzyⅠsodata聚類法在地下水水化學(xué)分類中的應(yīng)用.西安地質(zhì)學(xué)院學(xué)報[J].1991,13(3):59-66.

    [10]陽正熙,吳塹虹,彭直興,等.地學(xué)數(shù)據(jù)分析教程(試用本)[M].北京:科學(xué)出版社.2008:81-95.

    [11]胥澤銀,郭科編著.多元統(tǒng)計方法及其程序設(shè)計[M].四川:四川科學(xué)技出版社.

    [12]漆繼紅,許模等.川東銅鑼山背斜-南溫泉背斜溫泉水力聯(lián)系分析[J].人民長江.2011,(42)11:75-79.

    [13]沈照理,等.水文地球化學(xué)基礎(chǔ)[M].北京:科學(xué)出版社.1993.

    [14]蘇金貴.模糊聚類分析在水化學(xué)分類中的應(yīng)用[J].遼寧地質(zhì).1985,(3):251-258.

    猜你喜歡
    華鎣山水化學(xué)水文地質(zhì)
    幸福湖泉域地下水水化學(xué)特征及其演化機(jī)制
    人民黃河(2025年1期)2025-01-16 00:00:00
    基于抽水試驗計算水文地質(zhì)參數(shù)的分析與確定
    基于GPRS實現(xiàn)自動化水文地質(zhì)抽水試驗
    秦皇島北部地表水和地下水水化學(xué)特征及水質(zhì)評價
    基于轉(zhuǎn)移概率的三維水文地質(zhì)結(jié)構(gòu)劃分
    華鎣山游擊隊遺址
    黨史文苑(2020年1期)2020-03-16 03:17:30
    水文地質(zhì)在工程地質(zhì)勘察中的有效應(yīng)用
    我用一生在華鎣山愛著
    散文詩世界(2016年5期)2016-11-22 02:08:59
    改善華鎣山森林立地條件措施淺析
    衡水市新近系館陶組地?zé)崴瘜W(xué)特征與分布規(guī)律研究
    精品少妇黑人巨大在线播放| videossex国产| 各种免费的搞黄视频| 交换朋友夫妻互换小说| 母亲3免费完整高清在线观看 | 人妻一区二区av| 亚洲人成网站在线观看播放| 久久这里有精品视频免费| 久久精品国产亚洲av天美| 亚洲人与动物交配视频| 国产高清有码在线观看视频| 欧美性感艳星| 欧美另类一区| 国产一区亚洲一区在线观看| 十分钟在线观看高清视频www| 国产伦精品一区二区三区视频9| 看非洲黑人一级黄片| 成人黄色视频免费在线看| 久久热精品热| 国产成人精品无人区| 日本免费在线观看一区| 性高湖久久久久久久久免费观看| 免费看光身美女| a级片在线免费高清观看视频| 内地一区二区视频在线| 在线观看www视频免费| 狂野欧美白嫩少妇大欣赏| 下体分泌物呈黄色| 久久 成人 亚洲| 黄色毛片三级朝国网站| 久久97久久精品| 国产成人av激情在线播放 | 一区二区三区乱码不卡18| av电影中文网址| 日韩欧美精品免费久久| 精品卡一卡二卡四卡免费| 亚洲高清免费不卡视频| 午夜福利视频精品| 久久午夜福利片| 日韩免费高清中文字幕av| 最后的刺客免费高清国语| 一本一本综合久久| 精品一区二区三区视频在线| 免费日韩欧美在线观看| 曰老女人黄片| 韩国av在线不卡| 国产精品久久久久久久久免| 国产欧美日韩一区二区三区在线 | 久久99精品国语久久久| 在线观看免费日韩欧美大片 | 黄色怎么调成土黄色| 久久久久网色| 国产 精品1| 日本黄色片子视频| 蜜桃久久精品国产亚洲av| 大片免费播放器 马上看| 黑人欧美特级aaaaaa片| 男女高潮啪啪啪动态图| 99九九线精品视频在线观看视频| 国产成人av激情在线播放 | 亚洲欧洲国产日韩| 国产一区二区三区av在线| 精品久久久久久久久亚洲| 日韩制服骚丝袜av| 啦啦啦在线观看免费高清www| 少妇的逼好多水| 久久久久人妻精品一区果冻| av国产精品久久久久影院| 国产av国产精品国产| 一二三四中文在线观看免费高清| 美女内射精品一级片tv| 免费日韩欧美在线观看| 少妇 在线观看| 特大巨黑吊av在线直播| 亚洲av在线观看美女高潮| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久亚洲中文字幕| 在线播放无遮挡| 欧美激情极品国产一区二区三区 | 99国产精品免费福利视频| 99re6热这里在线精品视频| 人妻制服诱惑在线中文字幕| 亚洲人成77777在线视频| 国产成人精品久久久久久| 亚洲一区二区三区欧美精品| 亚洲色图 男人天堂 中文字幕 | 亚洲精品一二三| 免费大片黄手机在线观看| 狠狠婷婷综合久久久久久88av| 满18在线观看网站| 久久久久久久国产电影| 在线观看美女被高潮喷水网站| 国产免费一区二区三区四区乱码| 热99国产精品久久久久久7| 国产精品一区www在线观看| 女的被弄到高潮叫床怎么办| 内地一区二区视频在线| av视频免费观看在线观看| 十八禁网站网址无遮挡| 夫妻性生交免费视频一级片| 亚洲国产精品一区二区三区在线| 母亲3免费完整高清在线观看 | 亚洲精品成人av观看孕妇| 国产 精品1| 欧美日韩亚洲高清精品| 午夜影院在线不卡| 考比视频在线观看| 一个人看视频在线观看www免费| 久久久久久久久久人人人人人人| 狠狠婷婷综合久久久久久88av| 欧美亚洲 丝袜 人妻 在线| 女人精品久久久久毛片| 亚洲国产精品一区三区| 精品卡一卡二卡四卡免费| 成年人午夜在线观看视频| 国产一级毛片在线| 国产色爽女视频免费观看| 欧美激情 高清一区二区三区| 久久久久久久久久人人人人人人| 国产精品 国内视频| 男女啪啪激烈高潮av片| 国产精品成人在线| 日本wwww免费看| 午夜久久久在线观看| 久久韩国三级中文字幕| 丝袜脚勾引网站| 国产爽快片一区二区三区| 免费少妇av软件| 熟女av电影| 欧美日韩综合久久久久久| 亚洲av国产av综合av卡| 色婷婷av一区二区三区视频| 久久人人爽av亚洲精品天堂| 婷婷色综合www| 久久人人爽人人片av| 99久久人妻综合| 久热这里只有精品99| www.av在线官网国产| 天堂中文最新版在线下载| 只有这里有精品99| 日韩制服骚丝袜av| 欧美另类一区| 午夜91福利影院| 最近2019中文字幕mv第一页| 自线自在国产av| 大又大粗又爽又黄少妇毛片口| 在线精品无人区一区二区三| 日本av手机在线免费观看| 国产男人的电影天堂91| 中国国产av一级| 亚洲欧美清纯卡通| 人妻少妇偷人精品九色| 特大巨黑吊av在线直播| av专区在线播放| 久久久国产精品麻豆| 国产精品国产三级专区第一集| 日本vs欧美在线观看视频| 狂野欧美激情性xxxx在线观看| 国产成人a∨麻豆精品| 国产精品久久久久久精品电影小说| 香蕉精品网在线| 久久精品熟女亚洲av麻豆精品| 欧美激情极品国产一区二区三区 | 人成视频在线观看免费观看| 少妇被粗大猛烈的视频| 亚洲欧洲日产国产| videossex国产| 啦啦啦啦在线视频资源| 在线观看一区二区三区激情| 成年人免费黄色播放视频| 亚洲第一区二区三区不卡| 这个男人来自地球电影免费观看 | 一区二区三区免费毛片| 一边摸一边做爽爽视频免费| 亚洲不卡免费看| 午夜视频国产福利| 日本-黄色视频高清免费观看| 欧美bdsm另类| 成人18禁高潮啪啪吃奶动态图 | 精品久久久噜噜| 人人妻人人澡人人爽人人夜夜| 制服诱惑二区| 自拍欧美九色日韩亚洲蝌蚪91| 卡戴珊不雅视频在线播放| 99热全是精品| 亚洲精品久久午夜乱码| 久久久久久久久久久丰满| 最近中文字幕2019免费版| 高清av免费在线| 久久人人爽人人片av| 欧美激情极品国产一区二区三区 | 久久久久精品久久久久真实原创| 亚洲欧美中文字幕日韩二区| 亚洲中文av在线| 女人精品久久久久毛片| 啦啦啦啦在线视频资源| 三级国产精品欧美在线观看| 黄色毛片三级朝国网站| 午夜视频国产福利| 久久久久久久久大av| 水蜜桃什么品种好| 欧美日韩综合久久久久久| 多毛熟女@视频| 18禁动态无遮挡网站| 最近中文字幕2019免费版| 欧美另类一区| 十八禁高潮呻吟视频| a级毛片在线看网站| 国精品久久久久久国模美| 国产亚洲精品第一综合不卡 | 国产国拍精品亚洲av在线观看| 丰满饥渴人妻一区二区三| 大陆偷拍与自拍| 国国产精品蜜臀av免费| 国产亚洲精品第一综合不卡 | 久久午夜综合久久蜜桃| 又黄又爽又刺激的免费视频.| 亚洲精品久久久久久婷婷小说| 毛片一级片免费看久久久久| 尾随美女入室| 少妇的逼好多水| 久热久热在线精品观看| 天堂中文最新版在线下载| 美女xxoo啪啪120秒动态图| 高清毛片免费看| 久久久久久久久久久久大奶| 久久午夜福利片| 国产精品国产av在线观看| 少妇丰满av| 亚洲欧美中文字幕日韩二区| 国产av精品麻豆| 久久久久人妻精品一区果冻| 成人综合一区亚洲| 久久久久精品久久久久真实原创| 人妻系列 视频| 精品国产一区二区三区久久久樱花| 国产亚洲最大av| 国产爽快片一区二区三区| 成人免费观看视频高清| 久久久久久久久大av| 国产乱来视频区| 午夜视频国产福利| 69精品国产乱码久久久| 99热6这里只有精品| 国产男女超爽视频在线观看| 校园人妻丝袜中文字幕| 午夜免费男女啪啪视频观看| xxx大片免费视频| 美女视频免费永久观看网站| 美女大奶头黄色视频| 久久韩国三级中文字幕| 亚洲,欧美,日韩| 各种免费的搞黄视频| 999精品在线视频| 欧美xxⅹ黑人| 99九九在线精品视频| 亚洲人成网站在线观看播放| 亚洲av福利一区| 午夜精品国产一区二区电影| 国产黄频视频在线观看| 午夜激情av网站| 久久女婷五月综合色啪小说| 欧美日韩视频高清一区二区三区二| 日本欧美视频一区| 美女国产视频在线观看| 香蕉精品网在线| 国产精品 国内视频| 亚洲高清免费不卡视频| 国产一区二区三区综合在线观看 | 国产黄色视频一区二区在线观看| 国产熟女欧美一区二区| 久久国产亚洲av麻豆专区| 99久久精品一区二区三区| 亚洲伊人久久精品综合| 丝袜在线中文字幕| 久久久午夜欧美精品| 大香蕉97超碰在线| 波野结衣二区三区在线| 久久99热这里只频精品6学生| 免费人成在线观看视频色| 人人妻人人澡人人看| 久久精品国产亚洲网站| 国产永久视频网站| av不卡在线播放| 九草在线视频观看| 人人澡人人妻人| 男女无遮挡免费网站观看| 亚洲美女视频黄频| 中文字幕亚洲精品专区| 国产精品一国产av| 久久精品国产亚洲网站| 22中文网久久字幕| 在现免费观看毛片| av.在线天堂| 卡戴珊不雅视频在线播放| 在线播放无遮挡| 自拍欧美九色日韩亚洲蝌蚪91| 自线自在国产av| 狂野欧美白嫩少妇大欣赏| 国产av码专区亚洲av| xxxhd国产人妻xxx| 麻豆成人av视频| 国产精品三级大全| 青春草视频在线免费观看| www.色视频.com| 蜜桃国产av成人99| 亚洲美女黄色视频免费看| 丝袜脚勾引网站| 国产精品久久久久久久电影| av国产久精品久网站免费入址| 国产在线视频一区二区| 人人妻人人添人人爽欧美一区卜| 亚洲精华国产精华液的使用体验| 国产精品欧美亚洲77777| 国产一区有黄有色的免费视频| 少妇丰满av| 国产69精品久久久久777片| 久久热精品热| 欧美日韩av久久| 肉色欧美久久久久久久蜜桃| 能在线免费看毛片的网站| av有码第一页| 中文字幕av电影在线播放| 熟女av电影| 夫妻午夜视频| 午夜视频国产福利| 免费av中文字幕在线| 黄色欧美视频在线观看| 亚洲国产欧美在线一区| tube8黄色片| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久久久免| 日本与韩国留学比较| 日日摸夜夜添夜夜爱| 亚洲第一av免费看| 久久久久视频综合| 亚洲av成人精品一区久久| www.色视频.com| 久久久久国产精品人妻一区二区| 国产成人精品福利久久| 日本av手机在线免费观看| 日韩视频在线欧美| 日本午夜av视频| 亚洲欧美中文字幕日韩二区| 欧美精品一区二区大全| 国产片特级美女逼逼视频| 一级毛片黄色毛片免费观看视频| 精品久久久精品久久久| 精品国产一区二区三区久久久樱花| 母亲3免费完整高清在线观看 | 又粗又硬又长又爽又黄的视频| 精品人妻一区二区三区麻豆| 成人国产av品久久久| 欧美日韩视频精品一区| 嘟嘟电影网在线观看| 只有这里有精品99| 国产无遮挡羞羞视频在线观看| 亚洲精品久久久久久婷婷小说| 日韩 亚洲 欧美在线| 欧美精品人与动牲交sv欧美| 少妇 在线观看| 日韩av不卡免费在线播放| 曰老女人黄片| 亚洲精品视频女| 久久久久国产精品人妻一区二区| 国产极品天堂在线| 久久精品久久久久久久性| 蜜臀久久99精品久久宅男| 亚洲五月色婷婷综合| 午夜日本视频在线| 国产亚洲一区二区精品| av免费观看日本| 国产精品蜜桃在线观看| 久久精品国产自在天天线| 18禁动态无遮挡网站| 欧美日韩精品成人综合77777| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 内地一区二区视频在线| 亚洲av.av天堂| 在线观看人妻少妇| 亚洲精品色激情综合| 国产日韩欧美亚洲二区| 18禁裸乳无遮挡动漫免费视频| 久久久国产欧美日韩av| 大码成人一级视频| 精品一品国产午夜福利视频| av天堂久久9| 亚洲欧美成人精品一区二区| 校园人妻丝袜中文字幕| 看十八女毛片水多多多| 最近手机中文字幕大全| 纵有疾风起免费观看全集完整版| 亚洲五月色婷婷综合| 国产精品国产三级专区第一集| 久久精品国产亚洲网站| 免费观看a级毛片全部| 三上悠亚av全集在线观看| 国产国语露脸激情在线看| 乱人伦中国视频| www.色视频.com| 久久久国产欧美日韩av| 午夜久久久在线观看| 国产av一区二区精品久久| 伦理电影免费视频| 久久精品国产亚洲av涩爱| 亚洲精品av麻豆狂野| 国产精品一区二区三区四区免费观看| 18禁在线无遮挡免费观看视频| 黑丝袜美女国产一区| 亚洲国产最新在线播放| 寂寞人妻少妇视频99o| 久久久久精品性色| 最新的欧美精品一区二区| 热re99久久精品国产66热6| 亚洲精品日韩av片在线观看| 久久精品国产a三级三级三级| 十分钟在线观看高清视频www| 天天影视国产精品| 精品一区在线观看国产| 欧美激情国产日韩精品一区| 黄色毛片三级朝国网站| 日本免费在线观看一区| 日韩 亚洲 欧美在线| 只有这里有精品99| 18禁观看日本| 黑人高潮一二区| 在线观看免费视频网站a站| 黑人猛操日本美女一级片| 国产伦理片在线播放av一区| 91久久精品国产一区二区成人| 亚洲成人av在线免费| 夜夜骑夜夜射夜夜干| 大话2 男鬼变身卡| 久久人人爽人人片av| 91aial.com中文字幕在线观看| 亚洲精品国产av蜜桃| 亚洲av.av天堂| 精品一区二区三卡| 18禁在线播放成人免费| 色94色欧美一区二区| 高清欧美精品videossex| 中文天堂在线官网| 超碰97精品在线观看| 国产国语露脸激情在线看| 精品亚洲成国产av| 有码 亚洲区| 久久av网站| 亚洲欧美日韩卡通动漫| 亚洲av在线观看美女高潮| 视频区图区小说| 97在线人人人人妻| 免费久久久久久久精品成人欧美视频 | 一本久久精品| 日韩一本色道免费dvd| 热re99久久精品国产66热6| 久久精品久久久久久噜噜老黄| 丁香六月天网| 男女高潮啪啪啪动态图| 人妻制服诱惑在线中文字幕| 我要看黄色一级片免费的| 天天躁夜夜躁狠狠久久av| 亚洲成人一二三区av| 亚洲国产av影院在线观看| 国产色婷婷99| 国产精品女同一区二区软件| 亚洲人成77777在线视频| 女性被躁到高潮视频| 蜜臀久久99精品久久宅男| 久久久久久久国产电影| 亚洲av福利一区| 在线播放无遮挡| 91久久精品国产一区二区成人| 免费观看av网站的网址| 亚洲精品成人av观看孕妇| 少妇被粗大猛烈的视频| 亚洲第一av免费看| 美女xxoo啪啪120秒动态图| 成人亚洲精品一区在线观看| a级毛片黄视频| 久久99蜜桃精品久久| av线在线观看网站| 一边亲一边摸免费视频| 麻豆成人av视频| 精品国产国语对白av| 男的添女的下面高潮视频| 亚洲av综合色区一区| 久久久久久伊人网av| 久久这里有精品视频免费| 久久久久久久精品精品| 精品熟女少妇av免费看| 国产精品国产三级国产专区5o| 又大又黄又爽视频免费| 亚洲三级黄色毛片| 在线观看美女被高潮喷水网站| 久久久久久久久久人人人人人人| 亚洲人与动物交配视频| 人人妻人人澡人人看| 大又大粗又爽又黄少妇毛片口| 在线观看美女被高潮喷水网站| 99精国产麻豆久久婷婷| 99热全是精品| 色婷婷av一区二区三区视频| 热99久久久久精品小说推荐| 美女内射精品一级片tv| 91午夜精品亚洲一区二区三区| 午夜精品国产一区二区电影| 高清av免费在线| 国产成人freesex在线| 我要看黄色一级片免费的| 免费高清在线观看视频在线观看| 国产深夜福利视频在线观看| 免费av中文字幕在线| 女人久久www免费人成看片| 老司机亚洲免费影院| 大香蕉久久成人网| av在线app专区| 中文字幕最新亚洲高清| 蜜桃在线观看..| av天堂久久9| 成人综合一区亚洲| 精品少妇黑人巨大在线播放| 精品久久蜜臀av无| 国产免费又黄又爽又色| 一本—道久久a久久精品蜜桃钙片| 国语对白做爰xxxⅹ性视频网站| 亚洲av二区三区四区| xxxhd国产人妻xxx| 在线精品无人区一区二区三| 欧美性感艳星| 亚洲在久久综合| 国产综合精华液| 久久精品熟女亚洲av麻豆精品| 18禁在线播放成人免费| av视频免费观看在线观看| 啦啦啦在线观看免费高清www| 久久久久久久国产电影| 99热国产这里只有精品6| 免费观看无遮挡的男女| 亚洲欧美日韩另类电影网站| 伊人亚洲综合成人网| 久久国内精品自在自线图片| 国产精品久久久久久久电影| 国产精品熟女久久久久浪| 日韩人妻高清精品专区| 狂野欧美白嫩少妇大欣赏| 22中文网久久字幕| 简卡轻食公司| 两个人的视频大全免费| 免费日韩欧美在线观看| 精品国产一区二区三区久久久樱花| 精品国产国语对白av| 国产av精品麻豆| 一区二区三区乱码不卡18| 99视频精品全部免费 在线| 午夜福利影视在线免费观看| 九九在线视频观看精品| 免费观看的影片在线观看| 亚洲欧美精品自产自拍| 亚洲国产色片| 婷婷色麻豆天堂久久| 晚上一个人看的免费电影| 亚洲国产精品专区欧美| 国产日韩欧美亚洲二区| 爱豆传媒免费全集在线观看| 美女中出高潮动态图| 日韩中文字幕视频在线看片| 久久久久久久国产电影| 日本vs欧美在线观看视频| 国产高清三级在线| 色哟哟·www| 久久亚洲国产成人精品v| 在线观看国产h片| a级片在线免费高清观看视频| 国产精品欧美亚洲77777| 两个人免费观看高清视频| 国产亚洲午夜精品一区二区久久| 欧美97在线视频| 大片电影免费在线观看免费| 精品人妻一区二区三区麻豆| 一级毛片 在线播放| 亚洲av成人精品一区久久| 亚洲av福利一区| 99久久精品一区二区三区| 黑人猛操日本美女一级片| 在线播放无遮挡| 国产黄片视频在线免费观看| 精品酒店卫生间| 免费观看无遮挡的男女| 久久精品人人爽人人爽视色| 国产免费一区二区三区四区乱码| 夜夜看夜夜爽夜夜摸| 国产片特级美女逼逼视频| 日韩一本色道免费dvd| 婷婷色av中文字幕| 免费播放大片免费观看视频在线观看| 在线观看人妻少妇| 精品久久久久久久久av| 在线播放无遮挡| 黄色毛片三级朝国网站| 肉色欧美久久久久久久蜜桃| 国产精品99久久99久久久不卡 | 中文欧美无线码| 日韩电影二区| 99九九在线精品视频| 亚洲熟女精品中文字幕| 两个人免费观看高清视频| 亚洲精品国产av成人精品| 少妇被粗大猛烈的视频| 亚洲欧美成人精品一区二区| 婷婷色综合www| 最新中文字幕久久久久| 91精品一卡2卡3卡4卡| 成年人午夜在线观看视频| www.色视频.com|