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

    錢營(yíng)孜煤礦深部地下水水化學(xué)特征及來(lái)源解析

    2022-08-27 08:42:06劉啟蒙彭位華汪子濤
    煤田地質(zhì)與勘探 2022年8期
    關(guān)鍵詞:煤礦模型研究

    陳 凱,劉啟蒙,劉 瑜,彭位華,汪子濤,趙 祥

    (1.安徽理工大學(xué) 地球與環(huán)境學(xué)院,安徽 淮南 232001;2.安徽理工大學(xué) 深部煤礦采動(dòng)響應(yīng)與災(zāi)害防控國(guó)家重點(diǎn)實(shí)驗(yàn)室,安徽 淮南 232001;3.宿州學(xué)院 礦井水資源化利用安徽普通高校重點(diǎn)實(shí)驗(yàn)室,安徽 宿州 232000;4.中國(guó)科學(xué)院青海鹽湖研究所,青海 西寧 810008)

    在淺部煤炭資源逐漸枯竭的形勢(shì)下,如何為深部煤炭的安全開(kāi)采保駕護(hù)航已成為相關(guān)領(lǐng)域的研究熱點(diǎn)[1]。突水作為礦井深部開(kāi)采面臨的主要災(zāi)害之一,帶來(lái)的經(jīng)驗(yàn)教訓(xùn)是極其慘痛的。2021 年初,國(guó)家礦山安全監(jiān)察局公布的“2020 年全國(guó)煤礦事故十大典型案例”中便有2 起重大透水事故,共造成18 人死亡,直接經(jīng)濟(jì)損失高達(dá)6 054.03 萬(wàn)元。多次實(shí)踐證明,一旦發(fā)生井下突水事故,快速和精準(zhǔn)地判斷突水來(lái)源進(jìn)而采取相應(yīng)措施是將損失最小化的關(guān)鍵所在。因此,充分了解礦區(qū)地下水中化學(xué)組分的來(lái)源和演化機(jī)制顯得尤為重要,因?yàn)檫@是突水來(lái)源識(shí)別的必要前提。

    目前,學(xué)者們通過(guò)多種方法(如離子比值[2-3]、相關(guān)性分析[4]、聚類分析[5]和主成分分析[2-4]等)探討了不同礦區(qū)的含水層中地下水離子的來(lái)源問(wèn)題,這些研究方法已經(jīng)較為成熟,在一定程度上揭示了地下水化學(xué)特征形成機(jī)制,但無(wú)法進(jìn)一步給出水源貢獻(xiàn)的定量信息,這制約了對(duì)礦區(qū)深部地下水離子來(lái)源乃至水文地球化學(xué)演化的深入理解。近年來(lái),以Unmix 為典型代表的受體模型在大氣顆粒物中PM2.5[6]和PM10[7]、水溶性離子[8]以及土壤/沉積物中重金屬元素[9]和多環(huán)芳烴[10]的定量來(lái)源解析方面獲得成功應(yīng)用。Huang Fang 等[11]較早地將Unmix 模型引入地表水研究,并成功量化了工業(yè)廢水、農(nóng)業(yè)活動(dòng)和城市徑流對(duì)錢塘江中氨氮、六價(jià)鉻、石油類和鉛等污染物的具體貢獻(xiàn)。在地下水研究中,筆者[12]嘗試?yán)肬nmix 模型分離人為因素和地質(zhì)因素對(duì)研究區(qū)淺部含水層地下水中離子的貢獻(xiàn),A.Nagaraju 等[13]和N.S.Rao 等[14]基于Unmix 模型查明了不同巖石礦物組分的溶解對(duì)印度南部Chittoor 和Wanaparthy 地區(qū)高氟地下水中常規(guī)離子的貢獻(xiàn),這些研究證實(shí)了Unmix 模型的可遷移性。

    淮北煤田是我國(guó)重要的產(chǎn)煤基地,長(zhǎng)期的采煤活動(dòng)破壞了原本穩(wěn)定的水文地質(zhì)結(jié)構(gòu),隨著開(kāi)采深度的增加,煤層底部的高承壓灰?guī)r水對(duì)煤礦開(kāi)采的威脅也越來(lái)越大。在這種情況下,深入理解灰?guī)r水的水化學(xué)形成作用尤為重要。因此,筆者以安徽省淮北煤田的錢營(yíng)孜煤礦為研究示范,在系統(tǒng)測(cè)試灰?guī)r水中主要離子濃度的基礎(chǔ)上,通過(guò)多種分析方法(相關(guān)性分析、聚類分析和離子比值法)探討常規(guī)水化學(xué)組分的形成機(jī)制,并進(jìn)一步利用Unmix 模型量化不同因素對(duì)各個(gè)離子的具體貢獻(xiàn)信息,以期為研究區(qū)的礦井水害防治和水資源利用提供理論支持。

    1 材料與方法

    1.1 研究區(qū)概況

    錢營(yíng)孜煤礦位于安徽省宿州市西南部,距離宿州市中心15 km,與鄒莊礦、桃園礦、祁南礦相鄰。井田呈不規(guī)則多邊形展布,南北長(zhǎng)約8.3 km,東西寬約6 km。礦井揭露的地層有奧陶系、石炭系、二疊系、古近系、新近系和第四系,厚度近2 400 m,含煤地層為二疊系上、下石盒子組和山西組,厚度約1 266 m,自上而下含10 個(gè)煤組,其中32和82煤為錢營(yíng)孜煤礦的主要可采煤層,5、62、72和10 煤為次要開(kāi)采煤層。

    礦井含水層分布情況如圖1 所示。據(jù)礦井水文地質(zhì)劃分報(bào)告,四含(松散層第四含水層)不太發(fā)育且與煤層露頭直接接觸部位較少,因而對(duì)主采煤層影響不大;砂巖富水性弱,且經(jīng)多次疏放,雖存在涌水的可能但對(duì)開(kāi)采威脅程度較低。重點(diǎn)研究石炭系太原組灰?guī)r(簡(jiǎn)稱太灰)含水層,厚度為108~195 m,巖性由灰?guī)r、泥巖、粉砂巖、細(xì)砂巖及薄煤層互層組成,正常情況下對(duì)開(kāi)采10 煤無(wú)直接充水影響,但由于該含水層水壓大、水量豐富,且與開(kāi)采煤層距離較近(圖1),當(dāng)遇到斷層或陷落柱時(shí),可對(duì)礦坑直接充水,需引起高度重視。在構(gòu)造方面,錢營(yíng)孜煤礦處于宿北斷裂(北部)、固鎮(zhèn)?長(zhǎng)豐斷裂(東部)和板橋斷裂(南部)組成的斷塊內(nèi),總體構(gòu)造形態(tài)為一較寬緩向南仰起的向斜,并被一系列SN向、NNE 向斷層切割,受區(qū)域構(gòu)造的嚴(yán)格控制。

    圖1 錢營(yíng)孜煤礦地層巖性和含水層分布Fig.1 Hydrogeological column of Qianyingzi Coal Mine

    1.2 樣品收集與測(cè)試

    本次共收集灰?guī)r放水孔水樣22 個(gè),系統(tǒng)測(cè)試水中常規(guī)離子(Na+、Ca2+、Mg2+、Cl?、)濃度和pH 值,測(cè)試方法如下:pH 值利用便攜式水質(zhì)測(cè)試儀現(xiàn)場(chǎng)直接測(cè)試,濃度通過(guò)0.05 mol/L 的稀鹽酸滴定所得,其余陰陽(yáng)離子濃度通過(guò)離子色譜儀測(cè)試(ICS-600-900,Thermo Fischer Scientific,USA)。此外,溶解性總固體(TDS) 由所有離子濃度之和減去二分之一倍碳酸氫根濃度計(jì)算所得。水樣的陰陽(yáng)離子濃度平衡誤差均在±5%以內(nèi),表明測(cè)試結(jié)果可靠[15]。

    1.3 統(tǒng)計(jì)分析

    1.3.1 相關(guān)性分析

    相關(guān)性分析是傳統(tǒng)的數(shù)理統(tǒng)計(jì)方法之一,其通過(guò)相關(guān)系數(shù)來(lái)衡量變量之間的關(guān)聯(lián)程度,在各個(gè)學(xué)科中被廣泛應(yīng)用[4]。由于數(shù)據(jù)為非正態(tài)分布,因此,本研究選擇Spearman 秩相關(guān)系數(shù)進(jìn)行分析,可以降低數(shù)據(jù)的不確定性,從而獲得更穩(wěn)定的分析結(jié)果。2 個(gè)變量之間的Spearman 秩相關(guān)系數(shù)(ρ)計(jì)算公式如下:

    式中:n為等級(jí)個(gè)數(shù);di為2 個(gè)變量之間第i個(gè)樣品的等級(jí)差。

    1.3.2 R 型聚類分析

    系統(tǒng)聚類(包括Q 型和R 型聚類)是聚類分析中最常用的算法之一。在分組未知的情況下,劃分為同一組的樣品/變量被認(rèn)為具有相似的組成/來(lái)源特征[5]。在礦區(qū)水化學(xué)的研究中,前者多用于不同含水層的水力聯(lián)系分析或突水水源識(shí)別,而后者通常用于探討水化學(xué)組分的來(lái)源。在z-scores 標(biāo)準(zhǔn)化基礎(chǔ)下,本次使用歐式距離對(duì)地下水中主要離子進(jìn)行R 型聚類分析,以期獲得更可靠的聚類結(jié)果。二維空間的歐式距離計(jì)算公式如下:

    式中:L為點(diǎn)(x1,y1)與點(diǎn)(x2,y2)之間的歐式距離;|X|為點(diǎn)(x,y)到原點(diǎn)的歐氏距離。

    1.3.3 Unmix 模型

    Unmix 模型于2003 年由R.C.Henry[16]提出,該模型將物質(zhì)濃度歸于不同源貢獻(xiàn)的線性組合,通過(guò)自?;€分辨率的幾何技術(shù)對(duì)源貢獻(xiàn)進(jìn)行非負(fù)約束,原理如下:

    式中:Cij為樣品j中參數(shù)i的實(shí)際質(zhì)量濃度;aik為源k中參數(shù)i的質(zhì)量分?jǐn)?shù);Skj為源k中樣品j的總量;E為源組成的標(biāo)準(zhǔn)偏差。在深部地下水化學(xué)的研究中,離子濃度主要與水巖相互作用有關(guān),因此,源的類型可以考慮為不同種類的水巖相互作用,而源的貢獻(xiàn)可以理解為不同程度的水巖相互作用。

    Piper 三線圖[17]通過(guò)AqQA 軟件(version 1.5)繪制,相關(guān)性網(wǎng)絡(luò)圖解利用Gephi 軟件[18](version 0.9.2,https://gephi.org)完成,源解析通過(guò)Unmix 軟件開(kāi)展,其余統(tǒng)計(jì)分析和可視化均基于開(kāi)源的R 平臺(tái)(version 4.0.2,https://www.r-project.org)完成。

    2 結(jié)果與討論

    2.1 水文地球化學(xué)特征

    錢營(yíng)孜煤礦灰?guī)r含水層中主要離子質(zhì)量濃度、TDS 和pH 值的統(tǒng)計(jì)結(jié)果見(jiàn)表1。由表1 可知,水樣的TDS 值較高,為1 945~5 292 mg/L,平均3 754 mg/L,這可能與封閉條件下長(zhǎng)時(shí)間的水巖相互作用有關(guān)。pH 值為7.4~8.8,平均8.1,表明水質(zhì)呈中性至弱堿性。常規(guī)離子質(zhì)量濃度平均值排序如下:(2 305 mg/L)、Na+(718 mg/L)、Ca2+(270 mg/L)、(239 mg/L)、Cl?(186 mg/L)、Mg2+(153 mg/L)。變異系數(shù)(CV)表示數(shù)據(jù)的離散程度,在環(huán)境污染物溯源研究中,具有高變異系數(shù)的污染因子通常代表受人為作用影響,而對(duì)于深部地下水的化學(xué)組分而言,這種情況更多反映的是水巖相互作用類型和程度的差異[19]。在本研究中,除Cl?和具有偏高的CV 值(分別為0.57和0.50)外,其他參數(shù)的CV 值均低于0.4 (表1),可能與含水層中局部巖石礦物組成差異有關(guān)。此外,基于陰陽(yáng)離子的毫克當(dāng)量濃度劃分的水化學(xué)類型結(jié)果如圖2 所示,陰離子中占絕對(duì)優(yōu)勢(shì),陽(yáng)離子中大多數(shù)水樣以Na+為主導(dǎo),Mg2+和Ca2+次之,此離子濃度規(guī)律與陳松等[20]在本研究區(qū)附近礦區(qū)開(kāi)展的灰?guī)r水化學(xué)研究結(jié)果較為相似。據(jù)此,22 組水樣中有21 組為SO4-Na 型(占比95%),1 組為SO4-Mg 型(占比5%)。

    圖2 灰?guī)r含水層水Piper 圖Fig.2 Piper diagram of the groundwater samples in limestone aquifer

    表1 錢營(yíng)孜煤礦灰?guī)r含水層中水質(zhì)參數(shù)統(tǒng)計(jì)結(jié)果Table 1 Statistics of water quality parameters in the limestone aquifer of Qiangyingzi Coal Mine

    2.2 水化學(xué)控制機(jī)制

    2.2.1 離子比值分析

    多數(shù)研究[2-5]表明,離子比值分析是探討地下水水化學(xué)組成控制因素的有效手段。由于地下水中Cl?濃度相對(duì)穩(wěn)定,c(Na+)/c(Cl?)常用于解釋地下水中Na+的來(lái)源,若地下水樣品的c(Na+)/c(Cl?)等于1,說(shuō)明鹽巖的溶解是Na+的唯一來(lái)源。由圖3a 可看出,所有灰?guī)r水樣品均位于1∶1 線的上方,這種情況表明除鹽巖外,Na+還具有其他來(lái)源,可能包括:(1) 硅酸鹽礦物風(fēng)化溶解使水中Na+富集;(2) 陽(yáng)離子交替吸附作用使巖土顆粒中的Na+置換水中的Ca2+和Mg2+,從而導(dǎo)致水中Na+濃度升高。從圖3c 和圖3d 可進(jìn)一步證實(shí)含水層中硅酸鹽礦物風(fēng)化溶解作用的存在。用于揭示地下水中Ca2+和Mg2+的來(lái)源,若其值等于1,說(shuō)明受碳酸鹽和硫酸鹽礦物溶解的控制。在本研究區(qū)中(圖3b),所有水樣均位于1∶1 線右側(cè),說(shuō)明Ca2+和Mg2+相對(duì)虧缺,表明上文所述的陽(yáng)離子交換作用和硅酸鹽礦物風(fēng)化溶解作用的存在是合理的。

    圖3 地下水中主要離子比例關(guān)系Fig.3 Relationships of major ions in groundwater

    2.2.2 相關(guān)性分析

    通過(guò)式(1)獲得離子間相關(guān)系數(shù)后,在置信度(雙側(cè))0.01 顯著性水平下提取高于臨界值(22 個(gè)樣品的相關(guān)系數(shù)臨界值為0.54)的相關(guān)系數(shù),進(jìn)一步通過(guò)Gephi軟件生成網(wǎng)絡(luò)圖,可以更直觀地顯示多變量之間的內(nèi)在聯(lián)系(圖4a)。在水文地球化學(xué)的研究中,通常認(rèn)為在顯著性水平下相關(guān)的變量之間具有相似的來(lái)源[21]。如圖4a 所示,與Na+、Ca2+和Mg2+均具有顯著的正相關(guān)關(guān)系,相關(guān)系數(shù)依次為0.682、0.680 和0.661。據(jù)陳陸望等[3]的研究顯示,在還原條件下,礦區(qū)地下水中濃度與黃鐵礦的溶解有關(guān)。在本研究區(qū)中,地下水呈中性至弱堿性(表1),因此,高濃度的可能來(lái)源于硫酸鹽礦物溶解[14],代表性化學(xué)方程式如下:

    除硫酸鹽的溶解外,灰?guī)r水中Ca2+和Mg2+還可能來(lái)源于碳酸鹽礦物的溶解和陽(yáng)離子交替吸附作用[22]。從圖4a 中可看出,Ca2+和Mg2+之間也具有顯著性正相關(guān)關(guān)系(0.595),這可能與地層中普遍存在的碳酸鹽礦物的溶解有關(guān),例如白云石的溶解會(huì)向水中同時(shí)釋放Ca2+和Mg2+,相關(guān)方程式如下:

    圖4 灰?guī)r水主要離子相關(guān)性分析和R 型聚類分析結(jié)果Fig.4 Results of correlation analysis and R-type cluster analysis of major ions in groundwater from limestone aquifer

    氯堿指標(biāo)(Chloro-Alkaline indices:CAI-1,CAI-2)通常用于表征陽(yáng)離子交換作用的方向和強(qiáng)度[19,23],其計(jì)算公式為:

    式中:γ為離子的當(dāng)量濃度。

    如圖5 所示,與淮南煤田灰?guī)r水中離子作用關(guān)系相似[23],本研究中CAI-1 和CAI-2 均小于0,CAI-1 為?12.5~?1.1,CAI-2 為?0.7~?0.2,說(shuō)明該含水層地下水中Ca2+和Mg2+與Na+發(fā)生了一定程度的交替吸附作用,導(dǎo)致Ca2+和Mg2+濃度降低,而Na+含量升高,此結(jié)論與陳陸望等[3]基于主成分分析獲得的結(jié)論相一致。此外,和Cl?的相關(guān)系數(shù)為0.698,從表1 中可以觀察到,這2 種離子具有較高的CV 值(≥0.50),這種情況代表了局部地區(qū)存在氯鹽礦物溶解和硅酸鹽礦物風(fēng)化的貢獻(xiàn)。

    圖5 地下水樣氯堿指數(shù)圖Fig.5 Chloro–Alkaline indices diagram of groundwater samples

    2.2.3 聚類分析

    基于R 型聚類的主要離子濃度分析結(jié)果如圖4b所示,在距離為8 之內(nèi),6 種常規(guī)離子被劃分為2組(組1 和組2),組1 包括Mg2+、Ca2+、和Na+,組2包括和Cl?。前人的研究已查明,皖北礦區(qū)地下水中常規(guī)水化學(xué)組成主要受硅酸鹽礦物風(fēng)化、碳酸鹽、硫酸鹽和氯化物溶解的影響,且在不同含水層中水巖相互作用貢獻(xiàn)存在一定差異[24]。因此,組1 可以解釋為硫酸鹽和碳酸鹽礦物(含水層中的主要組成礦物)的溶解,因?yàn)槊⑾?、瀉鹽、石膏和白云石等礦物的溶解會(huì)向水中釋放大量的Na+、Ca2+和Mg2+。而組2代表氯鹽礦物的溶解和硅酸鹽類礦物的風(fēng)化,如鹽巖和長(zhǎng)石石英砂巖等的溶解。此分析結(jié)果與離子比值分析和相關(guān)性分析的結(jié)果可相互佐證。

    2.3 離子來(lái)源定量解析

    Unmix 模型具有較強(qiáng)的數(shù)據(jù)敏感性,若數(shù)據(jù)中缺失值或異常值過(guò)多可能會(huì)導(dǎo)致解析效果不佳。在本研究中,Unmix 模型運(yùn)行結(jié)果顯示Min Rsq 為0.95,Min Sig/Noise 為2.58,顯著高于模型的最低使用要求(Min Rsq >0.8 且Min Sig/Noise >2),表明解析結(jié)果合理可信[16]。

    從圖6 可看出,Unmix 模型成功分離了2 個(gè)源:源1 和源2。進(jìn)一步可以看出,2 個(gè)源對(duì)各個(gè)離子濃度的貢獻(xiàn):源1 對(duì)地下水中Ca2+、Mg2+、Na+和濃度的貢獻(xiàn)率較高,依次為73%、68%、63%和73%,相反,源2 對(duì)地下水中濃度的貢獻(xiàn)更高,分別為75%和66%。結(jié)合前文多元統(tǒng)計(jì)分析結(jié)果,源1 代表硫酸鹽和碳酸鹽礦物的溶解,源2 則反映了鹽巖的溶解和硅酸鹽礦物的風(fēng)化,且源1 對(duì)地下水離子濃度的貢獻(xiàn)(占比56%)要顯著高于源2(占比44%),這與地層中巖石礦物組成特征是基本一致的。此外,這種貢獻(xiàn)特征與筆者先前在安徽任樓煤礦開(kāi)展的研究結(jié)果相似(硫酸鹽礦物的溶解對(duì)灰?guī)r水中離子濃度貢獻(xiàn)高于50%)[22],證實(shí)了該結(jié)果的可靠性。

    圖6 基于Unmix 模型的主要離子濃度來(lái)源解析結(jié)果Fig.6 Source apportionment results for major ion concentrations based on Unmix model

    綜上所述,水巖相互作用是錢營(yíng)孜煤礦灰?guī)r水中常規(guī)離子的主要來(lái)源,其中硫酸鹽和碳酸鹽礦物的溶解導(dǎo)致了水中、Ca2+和Mg2+富集,隨著地下水徑流過(guò)程中陽(yáng)離子交替吸附作用的存在,水中Ca2+和Mg2+逐步被圍巖顆粒中Na+置換,最終形成了以SO4-Na 型為特征的高礦化度水體。

    3 結(jié) 論

    a.地下水中TDS 平均值達(dá)3 754 mg/L,水質(zhì)呈中至弱堿性,常規(guī)離子平均濃度從高到低排序?yàn)椋篘a+、Ca2+、、Cl?、Mg2+,SO4-Na 是主要的水化學(xué)類型。

    b.通過(guò)離子比值分析和多元統(tǒng)計(jì)分析查明了錢營(yíng)孜煤礦太灰含水層主要的水文地球化學(xué)作用為硫酸鹽和碳酸鹽礦物的溶解、鹽巖的溶解、硅酸鹽礦物的風(fēng)化和陽(yáng)離子交替吸附作用。

    c.應(yīng)用Unmix 模型對(duì)灰?guī)r水中離子濃度的反向溯源結(jié)果顯示,硫酸鹽和碳酸鹽礦物的溶解對(duì)Ca2+、Mg2+、Na+和濃度的貢獻(xiàn)率更高,而鹽巖的溶解和硅酸鹽礦物的風(fēng)化對(duì)Cl?和濃度的貢獻(xiàn)率更高,解析結(jié)果量化了不同水文地球化學(xué)作用的具體貢獻(xiàn)率。

    d.文中在傳統(tǒng)水化學(xué)成因解釋的基礎(chǔ)上進(jìn)一步結(jié)合受體模型探討了地下水中離子濃度的定量來(lái)源,研究成果為煤礦區(qū)深部水害防治和水資源利用提供理論支持。

    猜你喜歡
    煤礦模型研究
    一半模型
    FMS與YBT相關(guān)性的實(shí)證研究
    遼代千人邑研究述論
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    視錯(cuò)覺(jué)在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    3D打印中的模型分割與打包
    大型煤礦自動(dòng)化控制系統(tǒng)的設(shè)計(jì)與應(yīng)用
    上半年確定關(guān)閉煤礦名單513處
    97人妻天天添夜夜摸| tube8黄色片| 毛片一级片免费看久久久久| 亚洲 欧美一区二区三区| 制服人妻中文乱码| 日本猛色少妇xxxxx猛交久久| 亚洲av男天堂| 菩萨蛮人人尽说江南好唐韦庄| 草草在线视频免费看| 少妇人妻久久综合中文| 国产成人91sexporn| 少妇人妻 视频| 日韩精品有码人妻一区| 成人黄色视频免费在线看| 哪个播放器可以免费观看大片| 国产一区二区三区av在线| 超色免费av| 男女边摸边吃奶| 只有这里有精品99| a级毛片在线看网站| 午夜日本视频在线| 熟女电影av网| 国产高清不卡午夜福利| 欧美+日韩+精品| 欧美激情国产日韩精品一区| av视频免费观看在线观看| 午夜免费观看性视频| 亚洲丝袜综合中文字幕| 国产精品麻豆人妻色哟哟久久| 欧美成人午夜免费资源| 亚洲国产精品国产精品| 日韩视频在线欧美| 久久国产精品男人的天堂亚洲 | 九草在线视频观看| 国产白丝娇喘喷水9色精品| 成人亚洲精品一区在线观看| a级片在线免费高清观看视频| 宅男免费午夜| 成年女人在线观看亚洲视频| 精品一品国产午夜福利视频| 国产色婷婷99| 久久婷婷青草| 午夜福利,免费看| 欧美激情 高清一区二区三区| 考比视频在线观看| 国产午夜精品一二区理论片| 蜜桃在线观看..| 岛国毛片在线播放| 亚洲经典国产精华液单| 在线天堂最新版资源| 久久人妻熟女aⅴ| 永久免费av网站大全| 中文天堂在线官网| 男的添女的下面高潮视频| 这个男人来自地球电影免费观看 | 欧美日韩成人在线一区二区| 一级,二级,三级黄色视频| 日产精品乱码卡一卡2卡三| 一区二区三区精品91| av免费观看日本| 少妇人妻久久综合中文| av免费在线看不卡| 我的女老师完整版在线观看| 国产国拍精品亚洲av在线观看| 亚洲精品国产色婷婷电影| 青春草亚洲视频在线观看| 免费观看无遮挡的男女| videossex国产| 色婷婷久久久亚洲欧美| 一级黄片播放器| 欧美3d第一页| 高清毛片免费看| 久久精品国产亚洲av天美| 精品亚洲成a人片在线观看| 国产一区有黄有色的免费视频| 高清在线视频一区二区三区| 涩涩av久久男人的天堂| 考比视频在线观看| 捣出白浆h1v1| 一级毛片黄色毛片免费观看视频| 国产白丝娇喘喷水9色精品| 欧美日韩视频精品一区| 大片免费播放器 马上看| a级片在线免费高清观看视频| 99久久精品国产国产毛片| 亚洲精品国产av成人精品| 成年人午夜在线观看视频| 欧美成人午夜精品| 国产欧美日韩一区二区三区在线| 伊人久久国产一区二区| 在线观看一区二区三区激情| 少妇精品久久久久久久| 久热这里只有精品99| 日韩在线高清观看一区二区三区| 99国产精品免费福利视频| 国产免费现黄频在线看| 天堂中文最新版在线下载| 视频在线观看一区二区三区| 欧美精品一区二区免费开放| 最新的欧美精品一区二区| 国产 一区精品| 日韩伦理黄色片| 日韩大片免费观看网站| 精品一区二区三区四区五区乱码 | 亚洲欧美一区二区三区黑人 | 在线看a的网站| 婷婷色麻豆天堂久久| 国产欧美日韩综合在线一区二区| 草草在线视频免费看| 国产一区亚洲一区在线观看| 一级a做视频免费观看| 综合色丁香网| 91在线精品国自产拍蜜月| 日韩免费高清中文字幕av| 熟女人妻精品中文字幕| 久久精品久久久久久久性| 晚上一个人看的免费电影| 久久久国产欧美日韩av| 少妇的逼水好多| av有码第一页| 中文字幕另类日韩欧美亚洲嫩草| 精品福利永久在线观看| 日本午夜av视频| 亚洲成国产人片在线观看| 激情五月婷婷亚洲| 卡戴珊不雅视频在线播放| 一级毛片我不卡| 美女国产高潮福利片在线看| 超碰97精品在线观看| 亚洲av在线观看美女高潮| 免费黄频网站在线观看国产| 日韩在线高清观看一区二区三区| 看免费成人av毛片| 亚洲中文av在线| 精品国产露脸久久av麻豆| 一级片免费观看大全| 人人妻人人澡人人看| 五月玫瑰六月丁香| 亚洲情色 制服丝袜| 亚洲精品久久午夜乱码| 在现免费观看毛片| 超碰97精品在线观看| 巨乳人妻的诱惑在线观看| 亚洲精品美女久久av网站| 最近最新中文字幕大全免费视频 | 男女免费视频国产| freevideosex欧美| 精品久久久久久电影网| 最近中文字幕高清免费大全6| 热99国产精品久久久久久7| 黄网站色视频无遮挡免费观看| 一个人免费看片子| 免费看光身美女| 亚洲欧美精品自产自拍| 欧美日韩视频精品一区| 一本色道久久久久久精品综合| 欧美少妇被猛烈插入视频| 欧美日韩国产mv在线观看视频| 最近2019中文字幕mv第一页| 精品99又大又爽又粗少妇毛片| 午夜免费观看性视频| videosex国产| 国产永久视频网站| 国产一区亚洲一区在线观看| 天堂8中文在线网| 麻豆乱淫一区二区| 两个人免费观看高清视频| 十八禁网站网址无遮挡| 欧美日韩一区二区视频在线观看视频在线| av在线观看视频网站免费| 亚洲av男天堂| 亚洲情色 制服丝袜| 中文字幕精品免费在线观看视频 | 交换朋友夫妻互换小说| 国产 一区精品| 久久久精品94久久精品| 美女内射精品一级片tv| 秋霞在线观看毛片| 十八禁高潮呻吟视频| 国产又色又爽无遮挡免| 免费av不卡在线播放| 男女边摸边吃奶| 日本vs欧美在线观看视频| 一边亲一边摸免费视频| 久久人人爽av亚洲精品天堂| 免费av中文字幕在线| 亚洲欧美清纯卡通| 超碰97精品在线观看| av视频免费观看在线观看| 一级毛片黄色毛片免费观看视频| 亚洲色图 男人天堂 中文字幕 | 99久久综合免费| 麻豆乱淫一区二区| 精品亚洲成国产av| 国产一区二区三区综合在线观看 | 性色avwww在线观看| 久久精品久久久久久久性| 国产亚洲av片在线观看秒播厂| 国产成人一区二区在线| 国产免费福利视频在线观看| 亚洲欧美成人综合另类久久久| √禁漫天堂资源中文www| 午夜91福利影院| 青春草国产在线视频| 国产淫语在线视频| 亚洲成色77777| 高清视频免费观看一区二区| 亚洲欧美成人精品一区二区| 国产免费福利视频在线观看| 极品人妻少妇av视频| 丰满迷人的少妇在线观看| 黄色视频在线播放观看不卡| 丁香六月天网| 国产一区二区激情短视频 | 欧美xxxx性猛交bbbb| 国产精品久久久av美女十八| 亚洲精品美女久久av网站| 国产精品国产三级专区第一集| 午夜激情久久久久久久| av国产久精品久网站免费入址| 欧美成人精品欧美一级黄| 久久久欧美国产精品| av.在线天堂| 婷婷色综合大香蕉| 国产精品欧美亚洲77777| 亚洲av免费高清在线观看| 边亲边吃奶的免费视频| 久久婷婷青草| 久久人人爽人人片av| 老女人水多毛片| 捣出白浆h1v1| 丝袜人妻中文字幕| 色哟哟·www| 丝袜脚勾引网站| 新久久久久国产一级毛片| 岛国毛片在线播放| 黄片无遮挡物在线观看| 美女xxoo啪啪120秒动态图| 水蜜桃什么品种好| 精品亚洲成a人片在线观看| 99国产综合亚洲精品| www.色视频.com| 精品卡一卡二卡四卡免费| 亚洲av在线观看美女高潮| 亚洲一码二码三码区别大吗| 亚洲婷婷狠狠爱综合网| 成年动漫av网址| 日韩精品有码人妻一区| 热re99久久精品国产66热6| 亚洲成人av在线免费| av天堂久久9| 国产亚洲av片在线观看秒播厂| 久久99精品国语久久久| 久久久久精品性色| 90打野战视频偷拍视频| 在线观看国产h片| 91成人精品电影| 国产精品国产三级专区第一集| av电影中文网址| 久久精品国产亚洲av涩爱| 高清av免费在线| 九九爱精品视频在线观看| 最新的欧美精品一区二区| 亚洲精品国产av成人精品| 久久久久久久久久成人| 一级毛片电影观看| 精品一品国产午夜福利视频| 亚洲成国产人片在线观看| 亚洲国产欧美日韩在线播放| 国产极品天堂在线| 国产在视频线精品| 亚洲经典国产精华液单| 国精品久久久久久国模美| 成人亚洲欧美一区二区av| 精品人妻熟女毛片av久久网站| 91久久精品国产一区二区三区| 国产精品一二三区在线看| 午夜福利乱码中文字幕| 欧美亚洲 丝袜 人妻 在线| 青青草视频在线视频观看| 精品99又大又爽又粗少妇毛片| 久久久国产精品麻豆| 日韩三级伦理在线观看| 国产精品久久久久久精品电影小说| 国产在线一区二区三区精| av在线app专区| 纵有疾风起免费观看全集完整版| av在线播放精品| 丁香六月天网| 欧美精品人与动牲交sv欧美| 国产高清国产精品国产三级| 免费黄网站久久成人精品| 黑丝袜美女国产一区| 一边亲一边摸免费视频| 精品一区二区三区视频在线| 一区二区日韩欧美中文字幕 | 99热6这里只有精品| 成年美女黄网站色视频大全免费| 99九九在线精品视频| 99国产精品免费福利视频| 色哟哟·www| 国产一区有黄有色的免费视频| 久久精品国产自在天天线| 国产精品国产三级国产专区5o| 国产免费一级a男人的天堂| 18禁在线无遮挡免费观看视频| 熟女av电影| 有码 亚洲区| 一个人免费看片子| 色吧在线观看| 男女午夜视频在线观看 | 国产一区有黄有色的免费视频| 日韩不卡一区二区三区视频在线| 久久久久久人人人人人| 如日韩欧美国产精品一区二区三区| 国产精品99久久99久久久不卡 | 性高湖久久久久久久久免费观看| 欧美日韩视频高清一区二区三区二| 国产高清三级在线| 亚洲av欧美aⅴ国产| 国产男人的电影天堂91| 18禁动态无遮挡网站| 最新的欧美精品一区二区| 97人妻天天添夜夜摸| av.在线天堂| 国产1区2区3区精品| 国产毛片在线视频| 国产日韩一区二区三区精品不卡| 亚洲精华国产精华液的使用体验| 2018国产大陆天天弄谢| 精品国产一区二区三区四区第35| av又黄又爽大尺度在线免费看| 九九爱精品视频在线观看| av福利片在线| 男人舔女人的私密视频| 成人国产麻豆网| 国产一级毛片在线| 波多野结衣一区麻豆| 国产精品一二三区在线看| 这个男人来自地球电影免费观看 | 久久韩国三级中文字幕| 亚洲国产日韩一区二区| 久久久久久久精品精品| 亚洲精华国产精华液的使用体验| 精品第一国产精品| 久久精品国产亚洲av涩爱| 国产精品女同一区二区软件| 亚洲精品一二三| 夫妻午夜视频| 久久这里有精品视频免费| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产毛片av蜜桃av| 日韩三级伦理在线观看| 亚洲国产看品久久| 色哟哟·www| h视频一区二区三区| 成人国语在线视频| 久久精品国产亚洲av天美| 国产精品一国产av| 国产精品久久久久久久久免| 五月伊人婷婷丁香| 久久久久久人妻| 欧美bdsm另类| 中文字幕免费在线视频6| 成人黄色视频免费在线看| 久久精品久久久久久久性| 丝袜人妻中文字幕| 免费看av在线观看网站| 日本午夜av视频| 黄色 视频免费看| 亚洲精品乱码久久久久久按摩| 熟女人妻精品中文字幕| 蜜臀久久99精品久久宅男| 老司机影院成人| 你懂的网址亚洲精品在线观看| 人妻人人澡人人爽人人| 女人久久www免费人成看片| 最近中文字幕高清免费大全6| 精品人妻熟女毛片av久久网站| 久久99热这里只频精品6学生| 在线精品无人区一区二区三| 最近中文字幕高清免费大全6| 久久久久国产网址| 在线观看免费高清a一片| 国产在线视频一区二区| 国精品久久久久久国模美| 少妇高潮的动态图| 亚洲国产最新在线播放| 你懂的网址亚洲精品在线观看| 欧美 亚洲 国产 日韩一| 香蕉丝袜av| av电影中文网址| 欧美xxⅹ黑人| 午夜福利视频精品| 午夜视频国产福利| 久久久国产一区二区| 男人添女人高潮全过程视频| 黑人猛操日本美女一级片| 国产不卡av网站在线观看| 9色porny在线观看| 亚洲国产日韩一区二区| 亚洲国产最新在线播放| 亚洲av日韩在线播放| 精品一区二区三区视频在线| 这个男人来自地球电影免费观看 | 最近2019中文字幕mv第一页| 日产精品乱码卡一卡2卡三| 在线观看www视频免费| 国产 一区精品| 国产精品蜜桃在线观看| 免费高清在线观看视频在线观看| 国产片特级美女逼逼视频| 看十八女毛片水多多多| 亚洲色图 男人天堂 中文字幕 | 最后的刺客免费高清国语| 国产一区二区三区综合在线观看 | 久久亚洲国产成人精品v| 大片电影免费在线观看免费| 精品人妻一区二区三区麻豆| 成人国产麻豆网| 啦啦啦啦在线视频资源| 久久99一区二区三区| 亚洲,一卡二卡三卡| 蜜桃在线观看..| 久久午夜综合久久蜜桃| 国产av一区二区精品久久| 国产激情久久老熟女| 巨乳人妻的诱惑在线观看| 欧美人与性动交α欧美精品济南到 | 热99久久久久精品小说推荐| 精品久久蜜臀av无| 新久久久久国产一级毛片| 菩萨蛮人人尽说江南好唐韦庄| 女的被弄到高潮叫床怎么办| 在线观看免费日韩欧美大片| 建设人人有责人人尽责人人享有的| 精品福利永久在线观看| 午夜免费鲁丝| 精品国产一区二区三区四区第35| 日日撸夜夜添| 国产在线一区二区三区精| 26uuu在线亚洲综合色| 99热这里只有是精品在线观看| 男女边摸边吃奶| 最后的刺客免费高清国语| 免费av不卡在线播放| 一二三四中文在线观看免费高清| 亚洲国产精品999| 自拍欧美九色日韩亚洲蝌蚪91| 女性被躁到高潮视频| 欧美日韩亚洲高清精品| 中文天堂在线官网| 精品久久久精品久久久| 免费人成在线观看视频色| av在线观看视频网站免费| 亚洲欧美清纯卡通| 免费高清在线观看日韩| 黄色毛片三级朝国网站| 亚洲国产精品999| 激情五月婷婷亚洲| 天天影视国产精品| 七月丁香在线播放| 天天躁夜夜躁狠狠躁躁| 黄片无遮挡物在线观看| 99热6这里只有精品| 22中文网久久字幕| 又黄又粗又硬又大视频| 99视频精品全部免费 在线| 一二三四中文在线观看免费高清| 亚洲,欧美,日韩| 成人18禁高潮啪啪吃奶动态图| 精品一品国产午夜福利视频| 久久久久视频综合| 国产亚洲精品第一综合不卡 | 国产无遮挡羞羞视频在线观看| 国产黄色免费在线视频| 18禁动态无遮挡网站| 18禁裸乳无遮挡动漫免费视频| 久热这里只有精品99| 青青草视频在线视频观看| 夫妻性生交免费视频一级片| 国产亚洲精品久久久com| 亚洲精品日韩在线中文字幕| 精品99又大又爽又粗少妇毛片| av又黄又爽大尺度在线免费看| 国产精品人妻久久久久久| 免费观看无遮挡的男女| 久久99一区二区三区| 久久精品久久久久久久性| 精品一区二区三区四区五区乱码 | 亚洲精品乱久久久久久| av视频免费观看在线观看| 久久精品久久精品一区二区三区| 日本-黄色视频高清免费观看| 最近的中文字幕免费完整| av福利片在线| 久久狼人影院| av国产精品久久久久影院| 国产av码专区亚洲av| 少妇 在线观看| 免费观看a级毛片全部| 丝袜喷水一区| 中国国产av一级| 夫妻性生交免费视频一级片| 久久精品夜色国产| 一区二区av电影网| 久久精品熟女亚洲av麻豆精品| 欧美日韩视频高清一区二区三区二| 色网站视频免费| 少妇人妻精品综合一区二区| 草草在线视频免费看| 亚洲三级黄色毛片| 亚洲情色 制服丝袜| 欧美日韩综合久久久久久| 97精品久久久久久久久久精品| 欧美人与性动交α欧美精品济南到 | 在线观看www视频免费| 亚洲av男天堂| 国产精品不卡视频一区二区| videossex国产| 1024视频免费在线观看| 久久精品国产a三级三级三级| 国国产精品蜜臀av免费| 色94色欧美一区二区| 曰老女人黄片| 看非洲黑人一级黄片| 亚洲三级黄色毛片| 久久97久久精品| 纵有疾风起免费观看全集完整版| 中文字幕亚洲精品专区| 国产亚洲欧美精品永久| 免费看光身美女| av不卡在线播放| 国产 一区精品| 巨乳人妻的诱惑在线观看| 亚洲精品久久午夜乱码| 日本黄大片高清| 亚洲国产欧美日韩在线播放| 久久久国产一区二区| 国产在线免费精品| 水蜜桃什么品种好| 精品一区二区三区四区五区乱码 | 自拍欧美九色日韩亚洲蝌蚪91| 十八禁网站网址无遮挡| av不卡在线播放| 亚洲天堂av无毛| 巨乳人妻的诱惑在线观看| av在线播放精品| 亚洲av欧美aⅴ国产| 久久人人97超碰香蕉20202| 精品久久国产蜜桃| 亚洲中文av在线| 欧美xxxx性猛交bbbb| 精品一区二区免费观看| 热99久久久久精品小说推荐| 一区二区三区乱码不卡18| 18在线观看网站| 国产日韩欧美亚洲二区| 亚洲av电影在线观看一区二区三区| 国产熟女欧美一区二区| 丰满饥渴人妻一区二区三| 99久久中文字幕三级久久日本| 亚洲,欧美,日韩| 国产老妇伦熟女老妇高清| 中文字幕人妻熟女乱码| 美女中出高潮动态图| 久久人人爽av亚洲精品天堂| 最近最新中文字幕大全免费视频 | 亚洲综合色网址| 捣出白浆h1v1| 熟女人妻精品中文字幕| 极品人妻少妇av视频| 日韩在线高清观看一区二区三区| 亚洲国产最新在线播放| 国产一区二区三区av在线| 久久狼人影院| 丝袜美足系列| 2022亚洲国产成人精品| 国产亚洲精品第一综合不卡 | 免费女性裸体啪啪无遮挡网站| 丝袜在线中文字幕| 亚洲欧洲日产国产| 成人亚洲欧美一区二区av| 免费看av在线观看网站| 九色成人免费人妻av| 久久久久久久久久久免费av| 秋霞在线观看毛片| 亚洲一级一片aⅴ在线观看| 久久久久精品人妻al黑| 亚洲欧美中文字幕日韩二区| 精品久久久久久电影网| 18+在线观看网站| 美国免费a级毛片| 国产熟女午夜一区二区三区| 午夜视频国产福利| 最近手机中文字幕大全| 男人操女人黄网站| 成人午夜精彩视频在线观看| 最近手机中文字幕大全| 最黄视频免费看| 日韩一区二区视频免费看| 性色avwww在线观看| 欧美人与性动交α欧美精品济南到 | av在线app专区| 久久99热这里只频精品6学生| 精品人妻偷拍中文字幕| av黄色大香蕉| 在线观看美女被高潮喷水网站| 纵有疾风起免费观看全集完整版| 精品亚洲乱码少妇综合久久| 美女xxoo啪啪120秒动态图| 天天操日日干夜夜撸| 男女午夜视频在线观看 | 男女免费视频国产| 亚洲精品国产av蜜桃|