• <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)槊⑾酢a鹽、石膏和白云石等礦物的溶解會(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處
    夜夜爽夜夜爽视频| 狠狠狠狠99中文字幕| a级毛片免费高清观看在线播放| 爱豆传媒免费全集在线观看| 亚洲精品乱久久久久久| 亚洲国产欧美在线一区| 国产成人精品婷婷| 国产亚洲91精品色在线| 色哟哟·www| 一边亲一边摸免费视频| 国产成人精品婷婷| 色哟哟·www| 青春草亚洲视频在线观看| 欧美一区二区亚洲| 亚洲国产欧美人成| 舔av片在线| 最近中文字幕高清免费大全6| 岛国在线免费视频观看| 精品久久久久久久末码| 亚洲综合精品二区| 在线免费十八禁| 久久人人爽人人片av| 久久久久性生活片| 嘟嘟电影网在线观看| 亚洲av不卡在线观看| 亚洲国产精品久久男人天堂| videos熟女内射| 中文字幕久久专区| 国产精品野战在线观看| 在线播放无遮挡| 久久久久国产网址| 天堂中文最新版在线下载 | 欧美zozozo另类| 国产国拍精品亚洲av在线观看| 国产亚洲91精品色在线| 日韩制服骚丝袜av| 国产伦精品一区二区三区视频9| 国产精品一及| 精品久久久久久久久亚洲| 日本色播在线视频| 91aial.com中文字幕在线观看| 日韩人妻高清精品专区| 91午夜精品亚洲一区二区三区| 看十八女毛片水多多多| 色网站视频免费| 国产精品久久久久久久电影| 一级黄色大片毛片| 岛国毛片在线播放| 精品国产露脸久久av麻豆 | 中文字幕亚洲精品专区| 国产v大片淫在线免费观看| 51国产日韩欧美| 天天躁日日操中文字幕| 亚洲国产最新在线播放| 能在线免费看毛片的网站| 精品久久久久久久末码| 亚洲欧美清纯卡通| 99久久精品国产国产毛片| 久99久视频精品免费| 午夜爱爱视频在线播放| 午夜福利视频1000在线观看| 男插女下体视频免费在线播放| 男插女下体视频免费在线播放| 免费观看的影片在线观看| 中文乱码字字幕精品一区二区三区 | 国语自产精品视频在线第100页| 午夜日本视频在线| 激情 狠狠 欧美| 亚洲最大成人手机在线| 欧美潮喷喷水| 亚洲av中文av极速乱| 久久人人爽人人爽人人片va| 日韩 亚洲 欧美在线| 精品久久久噜噜| 国产白丝娇喘喷水9色精品| 97人妻精品一区二区三区麻豆| 国产精品日韩av在线免费观看| 99久国产av精品国产电影| 天堂影院成人在线观看| av在线蜜桃| 久久久久久久久久黄片| 日韩亚洲欧美综合| 大又大粗又爽又黄少妇毛片口| 国产精品久久久久久久电影| 国产高清视频在线观看网站| 能在线免费观看的黄片| 联通29元200g的流量卡| 精品一区二区三区人妻视频| 国产亚洲91精品色在线| 51国产日韩欧美| 国产一区二区在线av高清观看| 老司机影院成人| 国产色爽女视频免费观看| 在线观看66精品国产| 国产亚洲5aaaaa淫片| 国产亚洲午夜精品一区二区久久 | 一夜夜www| 成年免费大片在线观看| 天天一区二区日本电影三级| 国产av在哪里看| 成年女人永久免费观看视频| 久久99热6这里只有精品| 小蜜桃在线观看免费完整版高清| 男的添女的下面高潮视频| 中文精品一卡2卡3卡4更新| 久久午夜福利片| 亚洲久久久久久中文字幕| 亚洲国产精品合色在线| 亚洲国产精品专区欧美| 欧美成人一区二区免费高清观看| 国产日韩欧美在线精品| 日本-黄色视频高清免费观看| 日日撸夜夜添| 中文乱码字字幕精品一区二区三区 | 免费观看精品视频网站| 在线观看av片永久免费下载| 男女国产视频网站| 久热久热在线精品观看| 日韩,欧美,国产一区二区三区 | 日韩欧美国产在线观看| 成人毛片60女人毛片免费| 久久精品影院6| 久久久久网色| 18+在线观看网站| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 最近的中文字幕免费完整| 成人欧美大片| 国产爱豆传媒在线观看| 一级毛片久久久久久久久女| 亚洲不卡免费看| 少妇高潮的动态图| 亚洲av男天堂| 亚洲国产精品成人久久小说| 国产欧美日韩精品一区二区| 村上凉子中文字幕在线| 99国产精品一区二区蜜桃av| av在线老鸭窝| 精品人妻偷拍中文字幕| 天美传媒精品一区二区| 色尼玛亚洲综合影院| 黄色配什么色好看| 黑人高潮一二区| 亚洲av免费高清在线观看| 夜夜爽夜夜爽视频| 日韩av在线免费看完整版不卡| 欧美性猛交黑人性爽| 国产人妻一区二区三区在| 韩国高清视频一区二区三区| 国产精品久久久久久精品电影| 九九热线精品视视频播放| 国产精品久久久久久精品电影小说 | 日韩一区二区三区影片| 美女被艹到高潮喷水动态| 国产白丝娇喘喷水9色精品| 国产一区二区三区av在线| 最近视频中文字幕2019在线8| 一本一本综合久久| 级片在线观看| 国语自产精品视频在线第100页| 成人三级黄色视频| 五月玫瑰六月丁香| 99国产精品一区二区蜜桃av| 亚洲国产精品成人久久小说| 直男gayav资源| 亚洲精华国产精华液的使用体验| 三级经典国产精品| 欧美三级亚洲精品| 日韩人妻高清精品专区| 人体艺术视频欧美日本| 亚洲在久久综合| 一级爰片在线观看| 欧美变态另类bdsm刘玥| 国产精品.久久久| 国产精品国产三级专区第一集| 美女黄网站色视频| 日韩av在线大香蕉| 男人和女人高潮做爰伦理| 国产亚洲一区二区精品| 久久精品熟女亚洲av麻豆精品 | 欧美又色又爽又黄视频| 精品久久久久久久末码| 精品久久久久久久人妻蜜臀av| 自拍偷自拍亚洲精品老妇| 免费观看性生交大片5| 在线播放国产精品三级| 国产三级中文精品| 天天躁日日操中文字幕| 国产在线男女| 亚洲伊人久久精品综合 | 黄色日韩在线| 一级毛片电影观看 | 91午夜精品亚洲一区二区三区| 一级毛片久久久久久久久女| 亚洲国产精品国产精品| 男人舔奶头视频| 天美传媒精品一区二区| 欧美97在线视频| 亚洲av中文av极速乱| 69人妻影院| 精品久久久久久久久久久久久| 三级国产精品欧美在线观看| 日本爱情动作片www.在线观看| 中文字幕免费在线视频6| 久久久久久久久久久免费av| 国产精品三级大全| 啦啦啦啦在线视频资源| 人体艺术视频欧美日本| 我要搜黄色片| 国产色婷婷99| 成人国产麻豆网| 毛片女人毛片| 一个人看视频在线观看www免费| 久久精品夜色国产| 欧美+日韩+精品| 国产精品av视频在线免费观看| 禁无遮挡网站| 高清视频免费观看一区二区 | 一区二区三区免费毛片| 成人午夜精彩视频在线观看| 亚洲美女视频黄频| 国产中年淑女户外野战色| 一个人看视频在线观看www免费| www.色视频.com| 1024手机看黄色片| 欧美激情国产日韩精品一区| 又粗又硬又长又爽又黄的视频| 人妻系列 视频| 少妇高潮的动态图| 国产视频首页在线观看| 久久99热6这里只有精品| 国产午夜精品一二区理论片| 大话2 男鬼变身卡| 亚洲一级一片aⅴ在线观看| 色播亚洲综合网| 麻豆精品久久久久久蜜桃| 赤兔流量卡办理| 亚洲自拍偷在线| 99国产精品一区二区蜜桃av| 亚洲精品乱码久久久v下载方式| av天堂中文字幕网| 免费电影在线观看免费观看| 久久国内精品自在自线图片| 国产亚洲av片在线观看秒播厂 | 久久草成人影院| 日韩成人伦理影院| 国产极品精品免费视频能看的| 91av网一区二区| 欧美性猛交黑人性爽| 亚洲最大成人手机在线| 免费看av在线观看网站| 一级爰片在线观看| 老司机影院毛片| 国产 一区精品| 九九热线精品视视频播放| 久久久久九九精品影院| 亚洲人成网站在线播| 国产精品不卡视频一区二区| 亚洲av日韩在线播放| 国产伦精品一区二区三区视频9| 一个人观看的视频www高清免费观看| 中文字幕制服av| av在线蜜桃| 狠狠狠狠99中文字幕| 欧美一级a爱片免费观看看| 国产一区二区三区av在线| 毛片女人毛片| 网址你懂的国产日韩在线| 成年女人永久免费观看视频| 亚洲天堂国产精品一区在线| 国产真实伦视频高清在线观看| 性色avwww在线观看| 国产亚洲精品av在线| 国产成人91sexporn| 久久人人爽人人片av| 一级黄色大片毛片| 久久精品久久精品一区二区三区| 啦啦啦观看免费观看视频高清| 亚洲成人中文字幕在线播放| 中文字幕精品亚洲无线码一区| 久久久久久久久久久免费av| 一级毛片久久久久久久久女| a级一级毛片免费在线观看| 欧美bdsm另类| 99热这里只有是精品在线观看| 插逼视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美潮喷喷水| 国产精品国产三级国产av玫瑰| 国产69精品久久久久777片| 亚洲国产高清在线一区二区三| 18禁在线播放成人免费| 免费观看精品视频网站| 美女国产视频在线观看| 国产69精品久久久久777片| 亚洲第一区二区三区不卡| 69人妻影院| av专区在线播放| 久久婷婷人人爽人人干人人爱| 欧美色视频一区免费| 中文字幕制服av| 综合色丁香网| 亚洲国产最新在线播放| 中文字幕制服av| 天天一区二区日本电影三级| 国产三级中文精品| 国产国拍精品亚洲av在线观看| 亚洲国产精品成人久久小说| 国产成人精品一,二区| 午夜福利在线观看免费完整高清在| 精品久久久久久久久亚洲| 深爱激情五月婷婷| a级毛色黄片| 免费av观看视频| 只有这里有精品99| 偷拍熟女少妇极品色| 晚上一个人看的免费电影| 91久久精品电影网| 午夜免费激情av| 不卡视频在线观看欧美| 激情 狠狠 欧美| 精品久久久噜噜| 免费播放大片免费观看视频在线观看 | 青春草国产在线视频| 国产精品一区二区性色av| 天天躁日日操中文字幕| 人妻制服诱惑在线中文字幕| 少妇裸体淫交视频免费看高清| 国产精品电影一区二区三区| 免费电影在线观看免费观看| 国产精品永久免费网站| 亚洲乱码一区二区免费版| 亚洲av免费高清在线观看| 久久久久久久久久久免费av| 欧美性感艳星| 18+在线观看网站| av又黄又爽大尺度在线免费看 | 国产精品福利在线免费观看| 热99re8久久精品国产| 久久精品91蜜桃| 亚洲国产精品成人综合色| 波多野结衣巨乳人妻| 黑人高潮一二区| 亚洲精品乱码久久久v下载方式| 村上凉子中文字幕在线| 久久精品国产亚洲网站| 国产女主播在线喷水免费视频网站 | 一本一本综合久久| 日本黄色视频三级网站网址| 精品久久久久久久末码| 两性午夜刺激爽爽歪歪视频在线观看| 狂野欧美激情性xxxx在线观看| 精品久久久噜噜| 国产成人91sexporn| 18禁动态无遮挡网站| 亚洲真实伦在线观看| 国产淫片久久久久久久久| 老司机影院毛片| av在线播放精品| 国产熟女欧美一区二区| av卡一久久| 中文资源天堂在线| 欧美性猛交黑人性爽| 国产日韩欧美在线精品| 国产一区亚洲一区在线观看| 日韩欧美国产在线观看| www.色视频.com| 亚洲av电影不卡..在线观看| av播播在线观看一区| 最近的中文字幕免费完整| 国产精品一区二区性色av| 国语对白做爰xxxⅹ性视频网站| 亚洲无线观看免费| 国产精品精品国产色婷婷| 可以在线观看毛片的网站| 亚洲自偷自拍三级| 老女人水多毛片| 日韩制服骚丝袜av| 欧美成人午夜免费资源| 国产不卡一卡二| 成人二区视频| 欧美日韩在线观看h| 午夜激情福利司机影院| 亚洲精品亚洲一区二区| 亚洲欧美一区二区三区国产| 国产熟女欧美一区二区| 午夜免费激情av| 亚洲av二区三区四区| 亚洲欧美日韩卡通动漫| 黄色一级大片看看| av在线播放精品| 搞女人的毛片| 日本免费a在线| 尾随美女入室| 偷拍熟女少妇极品色| 国产精品人妻久久久久久| 亚洲欧洲日产国产| 午夜视频国产福利| 久久久精品94久久精品| a级一级毛片免费在线观看| 久久精品影院6| 亚洲av成人精品一区久久| 亚洲av男天堂| 最近最新中文字幕大全电影3| 国产成人91sexporn| 内地一区二区视频在线| 日韩制服骚丝袜av| 精品人妻一区二区三区麻豆| 亚洲综合色惰| 高清毛片免费看| 亚洲美女视频黄频| 精品久久久久久久久久久久久| 国产精品av视频在线免费观看| 亚洲av电影不卡..在线观看| 久久精品久久久久久噜噜老黄 | av又黄又爽大尺度在线免费看 | 婷婷六月久久综合丁香| 国产在线一区二区三区精 | 超碰97精品在线观看| 蜜臀久久99精品久久宅男| 狂野欧美白嫩少妇大欣赏| 国产精品综合久久久久久久免费| 亚洲欧美一区二区三区国产| 久久欧美精品欧美久久欧美| 久久99热6这里只有精品| 午夜福利在线在线| 99热这里只有是精品在线观看| 汤姆久久久久久久影院中文字幕 | 亚洲精品乱码久久久v下载方式| 在线观看一区二区三区| 99热这里只有精品一区| 亚洲怡红院男人天堂| av线在线观看网站| 国产私拍福利视频在线观看| 国产高清三级在线| 亚洲国产欧美在线一区| 亚洲三级黄色毛片| 亚洲国产欧美人成| 久久亚洲国产成人精品v| 欧美性感艳星| 久久久久久久久中文| 最近最新中文字幕免费大全7| 久久精品国产99精品国产亚洲性色| 在线播放国产精品三级| 一级毛片aaaaaa免费看小| 欧美日韩精品成人综合77777| 亚洲成人久久爱视频| 日韩,欧美,国产一区二区三区 | 亚洲成人久久爱视频| 九九在线视频观看精品| 99久久无色码亚洲精品果冻| 精品少妇黑人巨大在线播放 | 赤兔流量卡办理| 99在线人妻在线中文字幕| 亚洲国产精品国产精品| av国产久精品久网站免费入址| 蜜桃亚洲精品一区二区三区| 欧美三级亚洲精品| av视频在线观看入口| 亚洲av中文字字幕乱码综合| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品人妻久久久影院| 好男人视频免费观看在线| videos熟女内射| 2021天堂中文幕一二区在线观| 欧美成人一区二区免费高清观看| 97超碰精品成人国产| av免费在线看不卡| 蜜桃亚洲精品一区二区三区| 久久久久久九九精品二区国产| 狠狠狠狠99中文字幕| 免费一级毛片在线播放高清视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产免费男女视频| 久久精品影院6| 一卡2卡三卡四卡精品乱码亚洲| 亚洲激情五月婷婷啪啪| 久久久久九九精品影院| 日本熟妇午夜| 亚洲国产欧洲综合997久久,| 国产精品久久电影中文字幕| 欧美精品一区二区大全| 国产成人aa在线观看| av在线观看视频网站免费| 超碰97精品在线观看| 哪个播放器可以免费观看大片| 高清毛片免费看| 亚洲一级一片aⅴ在线观看| 大又大粗又爽又黄少妇毛片口| 久久精品综合一区二区三区| 日日摸夜夜添夜夜添av毛片| 最近中文字幕高清免费大全6| 国语自产精品视频在线第100页| 一本久久精品| 天天躁夜夜躁狠狠久久av| 成年免费大片在线观看| 日韩精品有码人妻一区| 国产 一区 欧美 日韩| 少妇的逼好多水| 免费av观看视频| 91在线精品国自产拍蜜月| 三级男女做爰猛烈吃奶摸视频| 99久国产av精品国产电影| 国内少妇人妻偷人精品xxx网站| 中文字幕精品亚洲无线码一区| 少妇人妻一区二区三区视频| 亚洲精品乱码久久久v下载方式| 亚洲精品乱久久久久久| 欧美变态另类bdsm刘玥| 日本爱情动作片www.在线观看| 国产精品1区2区在线观看.| 久久国内精品自在自线图片| 午夜日本视频在线| 晚上一个人看的免费电影| 99热网站在线观看| 欧美潮喷喷水| 午夜亚洲福利在线播放| 色网站视频免费| 成人av在线播放网站| 2021天堂中文幕一二区在线观| 99热这里只有精品一区| 亚洲国产欧洲综合997久久,| 国产亚洲精品av在线| 亚洲aⅴ乱码一区二区在线播放| 高清午夜精品一区二区三区| 亚洲欧美一区二区三区国产| 日本午夜av视频| 国产精品永久免费网站| 99在线人妻在线中文字幕| 国产伦理片在线播放av一区| 天天躁日日操中文字幕| 人妻夜夜爽99麻豆av| 一区二区三区四区激情视频| 在线观看66精品国产| 欧美色视频一区免费| 只有这里有精品99| 亚洲精品一区蜜桃| 床上黄色一级片| 麻豆成人av视频| 好男人视频免费观看在线| 边亲边吃奶的免费视频| 成人二区视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜免费激情av| 国产三级中文精品| 免费黄网站久久成人精品| 亚洲精品乱久久久久久| 在线播放无遮挡| 国产精品久久久久久av不卡| 午夜福利网站1000一区二区三区| 超碰97精品在线观看| 99热网站在线观看| 午夜视频国产福利| 亚洲国产精品成人综合色| 日本黄大片高清| 欧美3d第一页| 国产精品,欧美在线| 波多野结衣巨乳人妻| 69av精品久久久久久| 国产免费福利视频在线观看| 亚洲一区高清亚洲精品| 久久久久久久久大av| 波多野结衣巨乳人妻| 一二三四中文在线观看免费高清| 国产精品一区二区性色av| 欧美丝袜亚洲另类| 七月丁香在线播放| 91精品一卡2卡3卡4卡| 嫩草影院入口| 亚洲激情五月婷婷啪啪| 欧美高清性xxxxhd video| 精品久久久久久久久av| 岛国毛片在线播放| 亚洲国产色片| 黄片wwwwww| 夜夜爽夜夜爽视频| 日本av手机在线免费观看| 国产精品久久久久久久电影| 日韩一区二区三区影片| 国产成人91sexporn| 日韩欧美在线乱码| 亚洲五月天丁香| 国产乱人视频| 九九爱精品视频在线观看| 有码 亚洲区| 白带黄色成豆腐渣| 国产黄色视频一区二区在线观看 | 免费人成在线观看视频色| 久久精品久久精品一区二区三区| 国产免费福利视频在线观看| 国产单亲对白刺激| 亚洲自偷自拍三级| 日韩国内少妇激情av| 老司机影院成人| 亚洲成人中文字幕在线播放| 免费观看精品视频网站| 亚洲欧美日韩卡通动漫| 国产爱豆传媒在线观看| 免费黄网站久久成人精品| 男女边吃奶边做爰视频| h日本视频在线播放| 狂野欧美白嫩少妇大欣赏| 久久婷婷人人爽人人干人人爱| 国产午夜精品一二区理论片| 国产精品久久视频播放| 日日啪夜夜撸| 日韩欧美国产在线观看| 欧美性感艳星| 天堂中文最新版在线下载 | 卡戴珊不雅视频在线播放| 纵有疾风起免费观看全集完整版 | 99热精品在线国产| 丝袜美腿在线中文| 精品人妻偷拍中文字幕| 久久久久久久久久久丰满|