• <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ī)律研究
    中国国产av一级| 午夜福利视频精品| 久久久久久免费高清国产稀缺| 免费黄网站久久成人精品| 日韩 欧美 亚洲 中文字幕| 色网站视频免费| 中文字幕人妻熟女乱码| 少妇猛男粗大的猛烈进出视频| 妹子高潮喷水视频| 久久久久精品人妻al黑| 成人黄色视频免费在线看| www.精华液| 欧美 亚洲 国产 日韩一| 韩国精品一区二区三区| 国产成人精品福利久久| bbb黄色大片| av天堂久久9| 我要看黄色一级片免费的| 乱人伦中国视频| 国产欧美亚洲国产| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人免费无遮挡视频| 国产无遮挡羞羞视频在线观看| 午夜免费观看性视频| 天堂中文最新版在线下载| 性少妇av在线| 毛片一级片免费看久久久久| 99久久人妻综合| 国产在线免费精品| 亚洲欧洲国产日韩| 男男h啪啪无遮挡| 欧美成人精品欧美一级黄| 尾随美女入室| 别揉我奶头~嗯~啊~动态视频 | 香蕉丝袜av| 观看av在线不卡| 精品一区二区三区av网在线观看 | 亚洲在久久综合| 这个男人来自地球电影免费观看 | 国产毛片在线视频| 免费高清在线观看视频在线观看| 在线观看国产h片| 国产精品99久久99久久久不卡 | 成人影院久久| 亚洲五月色婷婷综合| 日韩欧美一区视频在线观看| 亚洲精品久久午夜乱码| 亚洲美女搞黄在线观看| 性高湖久久久久久久久免费观看| 日韩精品有码人妻一区| av电影中文网址| 伊人久久国产一区二区| 狠狠婷婷综合久久久久久88av| 久久这里只有精品19| 亚洲成人手机| av在线老鸭窝| 国产成人精品无人区| 亚洲av福利一区| 老司机影院成人| 在线观看一区二区三区激情| 制服丝袜香蕉在线| 日本猛色少妇xxxxx猛交久久| 丝袜脚勾引网站| 人妻 亚洲 视频| 精品福利永久在线观看| 中国国产av一级| 国产乱来视频区| 欧美变态另类bdsm刘玥| 日韩av不卡免费在线播放| av一本久久久久| 国产日韩欧美视频二区| 看免费成人av毛片| 成人国产av品久久久| 国产不卡av网站在线观看| 男女下面插进去视频免费观看| av网站在线播放免费| 美女国产高潮福利片在线看| 免费在线观看完整版高清| 亚洲精品,欧美精品| 9色porny在线观看| 男人舔女人的私密视频| 午夜福利在线免费观看网站| 亚洲综合色网址| 成年人午夜在线观看视频| 男人添女人高潮全过程视频| 日韩av在线免费看完整版不卡| 丰满乱子伦码专区| 久久久久久人人人人人| 在线天堂最新版资源| 亚洲成国产人片在线观看| 亚洲自偷自拍图片 自拍| 美女中出高潮动态图| 久久久久久久久久久久大奶| 天美传媒精品一区二区| 丝瓜视频免费看黄片| 99精国产麻豆久久婷婷| 国产成人欧美在线观看 | 老熟女久久久| 国产女主播在线喷水免费视频网站| 国产男女内射视频| 亚洲精品在线美女| 成人亚洲欧美一区二区av| 爱豆传媒免费全集在线观看| 麻豆乱淫一区二区| 精品酒店卫生间| 精品一区二区三区四区五区乱码 | 免费日韩欧美在线观看| 蜜桃在线观看..| 只有这里有精品99| 丝袜美足系列| 久久狼人影院| 男人操女人黄网站| 亚洲综合精品二区| 在线观看免费视频网站a站| 午夜日本视频在线| 一二三四中文在线观看免费高清| 国产色婷婷99| 午夜av观看不卡| 高清不卡的av网站| 巨乳人妻的诱惑在线观看| 一级片免费观看大全| 国产高清国产精品国产三级| 韩国精品一区二区三区| 精品福利永久在线观看| 久久97久久精品| 亚洲av在线观看美女高潮| 中文字幕人妻丝袜制服| 中文字幕最新亚洲高清| 国产深夜福利视频在线观看| 电影成人av| 黑人欧美特级aaaaaa片| 18禁观看日本| 纵有疾风起免费观看全集完整版| 两个人看的免费小视频| 精品人妻一区二区三区麻豆| 婷婷色麻豆天堂久久| 午夜免费鲁丝| 在线观看一区二区三区激情| 青草久久国产| 美国免费a级毛片| 亚洲精品久久成人aⅴ小说| 超色免费av| 日日爽夜夜爽网站| 精品一品国产午夜福利视频| av在线观看视频网站免费| 啦啦啦视频在线资源免费观看| 男女午夜视频在线观看| 国产女主播在线喷水免费视频网站| 成人影院久久| 亚洲伊人色综图| 最近最新中文字幕大全免费视频 | 亚洲熟女毛片儿| 国产在线视频一区二区| 爱豆传媒免费全集在线观看| 亚洲精品中文字幕在线视频| 久久久国产精品麻豆| 可以免费在线观看a视频的电影网站 | 蜜桃国产av成人99| 在线看a的网站| 制服人妻中文乱码| 国产亚洲精品第一综合不卡| 国产精品偷伦视频观看了| 国产免费视频播放在线视频| 秋霞伦理黄片| 精品久久久久久电影网| 亚洲欧美激情在线| videosex国产| 国产亚洲最大av| 中文字幕制服av| 肉色欧美久久久久久久蜜桃| 91老司机精品| 美女中出高潮动态图| 黄片无遮挡物在线观看| 亚洲欧美一区二区三区久久| 人妻一区二区av| 色婷婷av一区二区三区视频| 免费看av在线观看网站| 最近最新中文字幕大全免费视频 | 另类精品久久| 大话2 男鬼变身卡| 亚洲精品av麻豆狂野| 国产精品久久久av美女十八| 免费人妻精品一区二区三区视频| 超碰成人久久| 色94色欧美一区二区| 欧美日韩亚洲国产一区二区在线观看 | av卡一久久| 久久精品熟女亚洲av麻豆精品| 两性夫妻黄色片| 亚洲婷婷狠狠爱综合网| 人人妻人人爽人人添夜夜欢视频| av网站在线播放免费| 婷婷成人精品国产| 天天操日日干夜夜撸| 精品一区二区三区av网在线观看 | 69精品国产乱码久久久| 日本欧美国产在线视频| 看免费成人av毛片| 国产高清国产精品国产三级| 午夜影院在线不卡| 亚洲av成人精品一二三区| 国产片特级美女逼逼视频| 人人妻人人爽人人添夜夜欢视频| 亚洲精品乱久久久久久| 亚洲专区中文字幕在线 | 亚洲精品国产区一区二| 国产成人精品久久二区二区91 | 国产精品国产av在线观看| 少妇人妻 视频| 97人妻天天添夜夜摸| 91精品三级在线观看| 中文字幕色久视频| 男人舔女人的私密视频| 国产精品偷伦视频观看了| 亚洲综合精品二区| 伦理电影大哥的女人| 日韩av免费高清视频| 如何舔出高潮| 一级,二级,三级黄色视频| 亚洲成人手机| 久久精品国产亚洲av涩爱| 一本—道久久a久久精品蜜桃钙片| 丰满饥渴人妻一区二区三| 国产av精品麻豆| 中文字幕最新亚洲高清| 日韩 欧美 亚洲 中文字幕| 99香蕉大伊视频| 欧美日韩精品网址| 免费观看人在逋| 黑人欧美特级aaaaaa片| 精品一品国产午夜福利视频| 一本色道久久久久久精品综合| 激情视频va一区二区三区| 久久人妻熟女aⅴ| 校园人妻丝袜中文字幕| 亚洲精品一二三| 建设人人有责人人尽责人人享有的| 啦啦啦 在线观看视频| 国产精品三级大全| 亚洲欧美色中文字幕在线| 亚洲天堂av无毛| 国产老妇伦熟女老妇高清| 国产成人精品福利久久| 青青草视频在线视频观看| 亚洲第一av免费看| 亚洲伊人久久精品综合| 老汉色∧v一级毛片| 精品福利永久在线观看| 免费在线观看视频国产中文字幕亚洲 | 韩国av在线不卡| 国产熟女欧美一区二区| 欧美日韩视频高清一区二区三区二| tube8黄色片| 婷婷色综合大香蕉| 一级a爱视频在线免费观看| 大陆偷拍与自拍| 国产视频首页在线观看| 一区福利在线观看| 久久久久视频综合| 国产在线一区二区三区精| 一区二区三区精品91| 国产精品熟女久久久久浪| 精品酒店卫生间| 成人影院久久| 国产av一区二区精品久久| e午夜精品久久久久久久| 天堂俺去俺来也www色官网| 国产福利在线免费观看视频| av电影中文网址| 老汉色av国产亚洲站长工具| 欧美国产精品一级二级三级| 国产亚洲精品第一综合不卡| 菩萨蛮人人尽说江南好唐韦庄| 在线观看www视频免费| 久久99一区二区三区| 精品少妇内射三级| 在线天堂中文资源库| 激情五月婷婷亚洲| 久久婷婷青草| 又粗又硬又长又爽又黄的视频| 日韩制服丝袜自拍偷拍| av福利片在线| 国产精品久久久久久久久免| 中文欧美无线码| 97人妻天天添夜夜摸| 中文字幕另类日韩欧美亚洲嫩草| 咕卡用的链子| 波野结衣二区三区在线| 国产av国产精品国产| 18禁国产床啪视频网站| 男女高潮啪啪啪动态图| 成人午夜精彩视频在线观看| 可以免费在线观看a视频的电影网站 | 国产精品无大码| 日韩 亚洲 欧美在线| 免费女性裸体啪啪无遮挡网站| 91精品伊人久久大香线蕉| 亚洲成av片中文字幕在线观看| 男女边吃奶边做爰视频| 成年动漫av网址| 中国三级夫妇交换| 日本91视频免费播放| 人人妻人人澡人人看| 久久久久久免费高清国产稀缺| 久久久亚洲精品成人影院| 老司机深夜福利视频在线观看 | 欧美人与善性xxx| 免费av中文字幕在线| 国产日韩欧美在线精品| 久久久久国产一级毛片高清牌| 乱人伦中国视频| av福利片在线| 蜜桃国产av成人99| 成人黄色视频免费在线看| 国产日韩欧美视频二区| 十八禁人妻一区二区| 午夜激情av网站| 亚洲欧美一区二区三区黑人| 亚洲自偷自拍图片 自拍| 久久女婷五月综合色啪小说| 毛片一级片免费看久久久久| 90打野战视频偷拍视频| 嫩草影院入口| 国产欧美亚洲国产| 黄色毛片三级朝国网站| 国产一区二区激情短视频 | www.精华液| 午夜影院在线不卡| 色综合欧美亚洲国产小说| 伊人久久国产一区二区| 七月丁香在线播放| tube8黄色片| 日日撸夜夜添| 波多野结衣av一区二区av| 精品一区二区三区四区五区乱码 | 巨乳人妻的诱惑在线观看| 一级毛片 在线播放| 黑人巨大精品欧美一区二区蜜桃| 少妇人妻 视频| 久久久久久人人人人人| 视频区图区小说| 女的被弄到高潮叫床怎么办| 婷婷色综合大香蕉| 天天躁夜夜躁狠狠躁躁| 日本wwww免费看| 少妇被粗大的猛进出69影院| 国产精品免费大片| 80岁老熟妇乱子伦牲交| 久久午夜综合久久蜜桃| 男女床上黄色一级片免费看| 午夜免费鲁丝| 国产成人精品福利久久| 亚洲婷婷狠狠爱综合网| 人人妻人人澡人人爽人人夜夜| 91精品国产国语对白视频| 国产在线视频一区二区| 国产1区2区3区精品| 老汉色av国产亚洲站长工具| 中文字幕av电影在线播放| 亚洲国产日韩一区二区| 中文字幕av电影在线播放| 看免费成人av毛片| 美女福利国产在线| 大陆偷拍与自拍| 丰满少妇做爰视频| 亚洲欧美一区二区三区黑人| 9热在线视频观看99| 欧美精品一区二区大全| 精品少妇一区二区三区视频日本电影 | 777米奇影视久久| 欧美精品高潮呻吟av久久| 久久久久视频综合| 久久精品熟女亚洲av麻豆精品| 欧美日韩亚洲国产一区二区在线观看 | 国产一级毛片在线| 99热国产这里只有精品6| 中文字幕亚洲精品专区| 亚洲国产av新网站| 亚洲成人免费av在线播放| av片东京热男人的天堂| 搡老岳熟女国产| 中文字幕色久视频| 久久精品aⅴ一区二区三区四区| 国产激情久久老熟女| 一区福利在线观看| 国产精品一区二区在线不卡| 国产探花极品一区二区| 热re99久久国产66热| 国产高清不卡午夜福利| 久久精品亚洲熟妇少妇任你| 天天躁夜夜躁狠狠久久av| 久久久久国产精品人妻一区二区| 一区二区三区精品91| 久久精品久久久久久久性| 两性夫妻黄色片| 飞空精品影院首页| 亚洲美女视频黄频| 一个人免费看片子| 国产亚洲午夜精品一区二区久久| 1024香蕉在线观看| 人人妻人人添人人爽欧美一区卜| 国产在线免费精品| 亚洲精品国产av成人精品| 亚洲欧美精品自产自拍| 久热爱精品视频在线9| 国产黄色免费在线视频| 婷婷成人精品国产| 亚洲国产av新网站| 国产日韩欧美在线精品| 激情五月婷婷亚洲| 精品人妻熟女毛片av久久网站| xxxhd国产人妻xxx| 波野结衣二区三区在线| 老司机亚洲免费影院| 久久久久久久久久久久大奶| 在线观看www视频免费| 别揉我奶头~嗯~啊~动态视频 | 天堂俺去俺来也www色官网| 国产一区二区三区av在线| 中文字幕制服av| 日韩一卡2卡3卡4卡2021年| 一个人免费看片子| 黑人猛操日本美女一级片| 欧美激情 高清一区二区三区| 久久久久久久精品精品| 久久女婷五月综合色啪小说| 欧美国产精品va在线观看不卡| 成人毛片60女人毛片免费| 免费黄色在线免费观看| 制服诱惑二区| 亚洲国产精品国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 精品视频人人做人人爽| 中国国产av一级| 国产一区二区三区综合在线观看| 99九九在线精品视频| 亚洲第一区二区三区不卡| 国产 精品1| 久久99精品国语久久久| 亚洲专区中文字幕在线 | 1024香蕉在线观看| 黄色毛片三级朝国网站| 狠狠婷婷综合久久久久久88av| 亚洲国产精品成人久久小说| 大香蕉久久成人网| 日韩制服丝袜自拍偷拍| 啦啦啦在线观看免费高清www| 日韩一区二区视频免费看| 久久久久人妻精品一区果冻| 国产成人精品福利久久| 免费在线观看视频国产中文字幕亚洲 | 青春草视频在线免费观看| 亚洲欧洲精品一区二区精品久久久 | 又大又爽又粗| 久久久久久人人人人人| 少妇的丰满在线观看| 两个人免费观看高清视频| 欧美日韩视频高清一区二区三区二| 性高湖久久久久久久久免费观看| 免费日韩欧美在线观看| 一区二区三区精品91| 男人舔女人的私密视频| 中文精品一卡2卡3卡4更新| 中文字幕亚洲精品专区| 中文字幕色久视频| 亚洲成人国产一区在线观看 | 飞空精品影院首页| 最新在线观看一区二区三区 | 91精品三级在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲一码二码三码区别大吗| 国产精品一二三区在线看| bbb黄色大片| 天天影视国产精品| 国产精品一区二区在线不卡| 成人亚洲精品一区在线观看| 亚洲国产精品成人久久小说| 亚洲色图 男人天堂 中文字幕| 久久人人97超碰香蕉20202| 国产精品人妻久久久影院| 赤兔流量卡办理| 国产99久久九九免费精品| 色94色欧美一区二区| 亚洲精品日本国产第一区| 久久精品亚洲av国产电影网| 亚洲少妇的诱惑av| 欧美少妇被猛烈插入视频| 十八禁网站网址无遮挡| 又黄又粗又硬又大视频| 亚洲精品久久久久久婷婷小说| 亚洲国产成人一精品久久久| xxx大片免费视频| 18禁国产床啪视频网站| 好男人视频免费观看在线| 中文天堂在线官网| 婷婷成人精品国产| 亚洲免费av在线视频| 考比视频在线观看| 午夜福利视频精品| 伊人亚洲综合成人网| 黄色 视频免费看| 亚洲精品乱久久久久久| 国产毛片在线视频| 最近手机中文字幕大全| 免费高清在线观看视频在线观看| 亚洲伊人色综图| √禁漫天堂资源中文www| 国产一区二区三区综合在线观看| av不卡在线播放| 国产精品偷伦视频观看了| 日日啪夜夜爽| 中文字幕最新亚洲高清| 老司机影院毛片| 美女国产高潮福利片在线看| 欧美人与性动交α欧美精品济南到| av线在线观看网站| 国产97色在线日韩免费| 天天影视国产精品| 激情五月婷婷亚洲| 午夜福利在线免费观看网站| 中文乱码字字幕精品一区二区三区| 永久免费av网站大全| 亚洲欧美中文字幕日韩二区| av免费观看日本| 国产日韩欧美在线精品| 亚洲精品久久午夜乱码| 在线看a的网站| 国产一级毛片在线| 久久久久国产一级毛片高清牌| 男女床上黄色一级片免费看| 超碰成人久久| av免费观看日本| 久久国产精品男人的天堂亚洲| 我的亚洲天堂| 90打野战视频偷拍视频| 欧美日韩国产mv在线观看视频| 免费观看人在逋| 高清av免费在线| 午夜免费观看性视频| 成人漫画全彩无遮挡| 捣出白浆h1v1| 一区二区三区四区激情视频| 欧美日本中文国产一区发布| 欧美人与性动交α欧美软件| 亚洲美女黄色视频免费看| 丝袜喷水一区| 观看美女的网站| 中文字幕另类日韩欧美亚洲嫩草| 国产无遮挡羞羞视频在线观看| 可以免费在线观看a视频的电影网站 | 国产精品人妻久久久影院| 在线精品无人区一区二区三| 不卡av一区二区三区| 免费观看a级毛片全部| 男人添女人高潮全过程视频| 国产成人免费无遮挡视频| 母亲3免费完整高清在线观看| 国产欧美日韩一区二区三区在线| 日本猛色少妇xxxxx猛交久久| 久久久久久久久久久免费av| 亚洲精品成人av观看孕妇| 日日撸夜夜添| 亚洲久久久国产精品| 热99国产精品久久久久久7| 哪个播放器可以免费观看大片| 久久狼人影院| av卡一久久| 国产黄色免费在线视频| 精品亚洲成a人片在线观看| 9色porny在线观看| 大话2 男鬼变身卡| 亚洲专区中文字幕在线 | 女性生殖器流出的白浆| 免费黄色在线免费观看| 亚洲免费av在线视频| 亚洲成av片中文字幕在线观看| 在线 av 中文字幕| 日本一区二区免费在线视频| 国产精品久久久人人做人人爽| 好男人视频免费观看在线| 一区二区三区精品91| 国产日韩欧美在线精品| 啦啦啦视频在线资源免费观看| 日韩av不卡免费在线播放| 国产日韩欧美亚洲二区| 人人妻人人添人人爽欧美一区卜| 纵有疾风起免费观看全集完整版| 亚洲图色成人| 久久人人爽av亚洲精品天堂| 国产成人精品无人区| 丝瓜视频免费看黄片| 久久热在线av| 国产成人精品久久二区二区91 | 十八禁人妻一区二区| 国产精品久久久av美女十八| 最近中文字幕高清免费大全6| 一级片免费观看大全| 国产日韩欧美亚洲二区| 黄片播放在线免费| 纵有疾风起免费观看全集完整版| 老熟女久久久| 国产成人av激情在线播放| 高清黄色对白视频在线免费看| 免费黄网站久久成人精品| 成人漫画全彩无遮挡| 高清黄色对白视频在线免费看| 亚洲图色成人| 91老司机精品| 亚洲国产欧美一区二区综合| 日日啪夜夜爽| 校园人妻丝袜中文字幕| a级毛片在线看网站| 在线观看人妻少妇| xxxhd国产人妻xxx| 伊人亚洲综合成人网|