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

    兩淮采煤沉陷積水區(qū)水體水化學特征及影響因素

    2014-06-07 05:55:14孫鵬飛易齊濤許光泉
    煤炭學報 2014年7期
    關(guān)鍵詞:水化學淮北淮南

    孫鵬飛,易齊濤,許光泉

    (安徽理工大學地球與環(huán)境學院,安徽淮南 232001)

    兩淮采煤沉陷積水區(qū)水體水化學特征及影響因素

    孫鵬飛,易齊濤,許光泉

    (安徽理工大學地球與環(huán)境學院,安徽淮南 232001)

    為分析沉陷積水區(qū)主要補給水源及相應(yīng)的溶質(zhì)來源,在淮南“潘謝”礦區(qū)和淮北“朱-楊莊”礦區(qū)各選取3個研究站點,基于水化學基本理論和原理,分析了各研究水域主要離子質(zhì)量濃度、組成及類型,采用Gibbs圖及因子分析方法,對水化學特征的影響因素進行了討論。結(jié)果表明:水體總?cè)芙庑怨腆w(TDS)質(zhì)量濃度在豐、平、枯水期總體上逐漸升高,淮北沉陷積水區(qū)水體的礦化度、堿度和總硬度均高于淮南沉陷積水區(qū)。水化學組成在區(qū)域范圍具有一定的空間變異性,HCO和Na+在陰陽離子中占有最大比例;淮南3個站點水化學類型主要為Na+-HCO-Cl-和Na+-Ca2+-HCO型,淮北3個站點則為Na+-Mg2+-HCO型。Gibbs圖和因子分析共同揭示了淮南沉陷積水區(qū)水化學特征主要受淺層地下水和地表徑流的雙重影響,而淮北沉陷積水區(qū)離子組成主要體現(xiàn)了區(qū)域淺層地下水化學的特征。

    采煤沉陷區(qū);水化學;因子分析;淮南礦區(qū);淮北礦區(qū)

    Key words:coalmines subsidence area;water chemistry;factor analysis;Huainan mining area;Huaibeimining area

    兩淮地區(qū)煤礦開采導致了廣泛的地表破壞和積水,現(xiàn)已導致300 km2左右的土地沉陷,未來10 a內(nèi)將擴展至500 km2以上,水資源儲存潛力可達109m3[1],對礦區(qū)可持續(xù)發(fā)展及生態(tài)環(huán)境等均產(chǎn)生重大影響。大面積的地表沉陷與積水改變了礦區(qū)原有的地形地貌,并使區(qū)域水交換條件和循環(huán)過程等發(fā)生變化。一方面,沉陷區(qū)形成巨大的儲水空間,攔蓄地表徑流;另一方面,水位的季節(jié)變化在有些地段導致地表水和淺層地下水發(fā)生水力聯(lián)系,進而影響水化學物質(zhì)濃度及平衡、污染物質(zhì)輸運和累積。對一般的河流、湖庫等水體而言,水化學組成及特征能夠揭示流域蒸發(fā)、巖石風化、大氣沉降及人為輸入的相互聯(lián)系[2-4],是研究水體環(huán)境特征和流域水文地球化學過程的重要方法。國內(nèi)已對長江、黃河[5-6]等主要河流,以及丹江口水庫、太湖[2,7]等不同湖泊水庫類水體主要離子質(zhì)量濃度、組成、來源及影響因素等方面進行了較為細致的研究。

    此前兩淮采煤沉陷積水區(qū)的研究主要集中于水體污染評價[8]、富營養(yǎng)化特征[1,9-10]以及水域沉積物[11-12]等方面,而關(guān)于水化學組分及特征缺乏相關(guān)研究報道。筆者將借助于河流、湖庫水化學的研究理論和方法,對兩淮采煤沉陷積水區(qū)的水化學特征進行科學表征,結(jié)合礦區(qū)特有的區(qū)域水文地質(zhì)條件,解析沉陷區(qū)水體主要的補給水源及相應(yīng)的溶質(zhì)來源,并討論相關(guān)影響因素,以期為礦區(qū)水資源規(guī)劃、水域管理與利用及生態(tài)重建與恢復等提供基礎(chǔ)資料和科學依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域概況

    兩淮礦區(qū)是全國14個大型煤炭基地之一,其中淮南礦區(qū)面積約3 600 km2,位于我國南北地理分界線附近,平均氣溫15.6℃,年均降雨量926 mm,礦區(qū)地表水系豐富,屬于淮河水系,淺層地下水位在1.5~3.0 m,與地表水之間水力聯(lián)系十分緊密。淮北礦區(qū)面積約9 600 km2,平均氣溫14.8℃,年均降雨量869.6 mm,區(qū)內(nèi)沱河、澮河、岱河、渦河等水系自西北向東南流過,最終匯入洪澤湖,淺層地下水埋深為1.0~3.0 m,水化學類型以HCO型為主。

    首先在兩淮選取兩個代表性礦區(qū),即淮南“潘謝”礦區(qū)和淮北“朱-楊莊”礦區(qū)(圖1)進行研究。

    圖1 兩淮礦區(qū)地理位置與水系分布Fig.1 Location and water system distribution of the Huainan&Huaibeimining areas

    “潘謝”礦區(qū)是現(xiàn)階段淮南井田分布最為密集、積水區(qū)面積最大的礦區(qū),在水生態(tài)環(huán)境重建與恢復中具有重要地位;“朱-楊莊”礦區(qū)是淮北生態(tài)城市建設(shè)的重要組成部分,是地質(zhì)生態(tài)環(huán)境治理投入最大、最具系統(tǒng)性的工程,該研究區(qū)域水體分別執(zhí)行著水源保護地、生態(tài)濕地、城市景觀等重要功能。筆者選取的兩個研究區(qū)域具有較強的對比性,能夠體現(xiàn)區(qū)域水化學的特征與差異。

    (1)淮南“潘謝”采煤沉陷積水區(qū)。

    “潘謝”礦區(qū)位于淮河中游北岸,淮北平原南緣,東西長70 km,南北寬25 km,面積約1 571 km2。區(qū)域地表水系發(fā)達,有西淝河、黑河、泥河和架河等河流,最終均匯入淮河,現(xiàn)沉陷積水區(qū)總面積約70 km2,平均水深3.0~5.0 m。在礦區(qū)的東、中、西部共設(shè)置3個研究站點,分別為潘集(PXS-1)、顧橋(PXS-3)和謝橋(PXS-5),具體如圖2(a)所示。

    (2)淮北“朱-楊莊”采煤沉陷積水區(qū)。

    淮北市現(xiàn)有采煤沉陷區(qū)30多處,總水域面積為32 km2左右,平均每年以3%~5%的速度增加[9]?!爸?楊莊”沉陷積水區(qū)位于淮北市東南方向,東西寬6 km,南北長10 km,根據(jù)“淮北市資源枯竭城市礦山地質(zhì)環(huán)境治理重點項目”,擬將東湖、中湖、南湖聯(lián)通,形成串珠狀湖泊系列,綜合整治水域面積超過10 km2,調(diào)蓄庫容2 476萬m3。區(qū)域內(nèi)的岱河(Dai River)自北向南與龍河(Long River)交匯至雷河(Lei River)。在該區(qū)設(shè)置3個研究站點,分別為東湖(HBDH)、中湖(HBZH)和南湖(HBNH),具體如圖2(b)所示。

    1.2 水質(zhì)采樣及分析方法

    結(jié)合水文條件及季節(jié)變化,分別于2011年8月(豐水期)、11月(平水期),2012年3月(枯水期),在PXS-1,PXS-3和PXS-5站點各布置5~9個采樣點采集表層水樣,具體位置同文獻[12]中的潘集(PJ)、顧橋(GQ)和謝橋(XQ)采樣點;于2012年8月(豐水期)、11月(平水期),2013年3月(枯水期)在HBDH,HBZH和HBNH站點設(shè)5~8個采樣點, HBNH采樣點同文獻[1],同時在Dai River,Long River,Lei River設(shè)3個斷面進行水質(zhì)采樣,方法參照相關(guān)規(guī)范進行[13]。為對比分析,于淮南市環(huán)境保護監(jiān)測站調(diào)取淮河淮南段國控斷面水質(zhì)監(jiān)測數(shù)據(jù),及淮南市2012年雨水監(jiān)測數(shù)據(jù)。

    現(xiàn)場使用野外便攜式Y(jié)SI儀測定水溫、pH等參數(shù)。采集的水樣分成兩份:一份采用HCl滴定法分析水體中HCO-3含量;另一份用孔徑0.45μm的醋酸纖維素濾膜進行過濾,濾樣裝入50 mL聚乙烯瓶中,用離子色譜儀(ICS1000,美國戴安公司)測定陽離子(K+,Na+,Ca2+,Mg2+)和陰離子(F-,Cl-,SO, NO)的質(zhì)量濃度。TDS質(zhì)量濃度(水體總?cè)芙夤腆w)采用各離子質(zhì)量濃度之和減去1/2 HCO質(zhì)量濃度的方法計算得到??傆捕群坎捎肅a2+,Mg2+含量換算成CaO含量的方法計算。實驗所用藥品均為優(yōu)級純,實驗用水為超純水,所有器皿均用0.3%HCl溶液浸泡24 h后,再用去離子水清洗。陰陽離子測試過程中,空白、標準同時測定,標準溶液為國家標準物質(zhì)研究中心標準溶液。為保證樣品測定準確度,對主要陰、陽離子平衡進行計算,其當量濃度比值在1.00±0.25[14]時,認為離子平衡,否則重新測定。

    1.3 數(shù)據(jù)處理及統(tǒng)計分析

    對各個研究站點,首先統(tǒng)計各離子指標(質(zhì)量濃度和毫克當量濃度)的均值和標準差,用SigmaPolt 10.0繪制水化學三角圖和Gibbs圖,水化學類型用AquaChem 3.7進行判斷。利用SPSS 19.0中的“因子分析”模塊,經(jīng)主成分分析及Varimax旋轉(zhuǎn),對各站點的離子來源進行解析。

    2 研究結(jié)果與分析

    2.1 水化學特征分析

    2.1.1 主要離子組成及含量

    淮南“潘謝”沉陷積水區(qū)主要離子組成見表1,3個研究站點的離子質(zhì)量濃度隨季節(jié)變化明顯,豐水期(FS)、平水期(PS)、枯水期(KS)的離子質(zhì)量濃度依次增大。TDS在上述3個時期的質(zhì)量濃度分別為192.5~453.0,303.9~527.2,308.4~684.6 mg/L,總體上在300~1 000 mg/L,屬于中等或較高礦化度水;3個時期總硬度在55~150 mg/L,屬于軟水;pH呈現(xiàn)弱堿性,均值為8.3。

    3個站點陽離子質(zhì)量濃度由大到小依次為:Na+, Ca2+,Mg2+,K+,豐水期均值分別為58.1,28.3,14.2和3.8 mg/L,平水期為72.8,38.9,17.5和3.9 mg/ L,枯水期則為110.5,42.1,21.7和4.8 mg/L。陰離子質(zhì)量濃度由大到小依次為:HCOCl-,SO,豐水期均值分別為168.7,63.4和53.5 mg/L,平水期為217.3,92.6和72.2 mg/L,枯水期則為235.3,127.4和116.7 mg/L。潘謝礦區(qū)3個研究站點陽離子以Na+,Ca2+為主,占陽離子總量的57.3%和26.6%,陰離子以HCO,Cl-為主,占陰離子總量的54.7%和24.2%。

    由于PXS-1站點和泥河連通,接受礦井水的排污,PXS-5站點和周圍的農(nóng)業(yè)主干渠道濟河連通,接納周圍農(nóng)業(yè)排水,兩個站點汛期受雨水沖淡作用明顯,平、枯水期濃度高于豐水期;較為封閉的PXS-3站點3個時期TDS含量和總硬度變化不大。同時沉陷積水區(qū)TDS含量和總硬度顯著高于臨近的淮河(僅給出枯水期離子濃度作為參考)。

    淮北“朱-楊莊”沉陷積水區(qū)主要離子組成見表2,3個湖泊陽離子質(zhì)量濃度總體順序大到小依次為: Na+,Mg2+,Ca2+,K+,均值分別為139.7,32.7,27.4和8.7 mg/L;陰離子由大到小為:HCO,SO,Cl-,總體均值分別為306.2,168.4和135.2 mg/L。其中,陽離子以Na+,Ca2+為主,占陽離子總量的59.9%和20.1%,陰離子則以HCO,SO為主,占陰離子總量的53.8%和27.5%。

    表1 淮南潘謝礦區(qū)沉陷積水區(qū)水體主要離子組成Table 1 The com position of major ions in the subsided waters of the Panxiem ining area mg/L

    表2 淮北朱-楊莊礦區(qū)沉陷積水區(qū)水體主要離子組成Table 2 The composition ofmajor ions in the subsided waters of the Zhu-Yangzhuangm ining area mg/L

    3個時期TDS質(zhì)量濃度均值空間變化為: HBDH(791.6 mg/L)>HBZH(643.9 mg/L)>HBNH(564.5 mg/L),Lei River(853.8 mg/L)>Long River(841.1 mg/L)>Dai River(615.8 mg/L),而總硬度的空間變化亦是HBNH最小。

    就兩淮礦區(qū)而言,淮北沉陷區(qū)水體TDS質(zhì)量濃度和總硬度明顯高于淮南沉陷積水區(qū)。而同一礦區(qū)小范圍內(nèi)各積水區(qū)之間的離子組成也存在一定的空間變異性。

    2.1.2 水化學類型特征

    水化學三角圖可以表明水體的主離子組成變化和特征,并辨別其控制端元[15]。在以Ca2+-Mg2+-(K++ Na+)組成的陽離子三角圖中(圖3(a)),受碳酸鹽巖影響的水體,其組分點會落在Ca2+-Mg2+線上,靠近Ca2+端元;受蒸發(fā)鹽巖影響,其組分點偏向(K++Na+)端元。在陰離子三角圖中也有類似反映(圖3(b)),流經(jīng)碳酸鹽巖地區(qū)的水體以HCO為主導,組分點落在HCO一端;流經(jīng)蒸發(fā)鹽巖地區(qū)的水體其組分點會落在SO-Cl-線上,遠離HCO一端。

    主要離子組成中,PXS-1,PXS-3,PXS-5站點之間存在一定差異(圖3),其中PXS-1受其與泥河相連通的影響,豐、平、枯水期離子變化范圍較大,而PXS-3和PXS-5相對封閉,其水化學類型在各個時期均較為穩(wěn)定。3個站點陽離子主要分布在(K++ Na+)-Ca2+線上,并靠近(K++Na+)端元,Mg2+含量較穩(wěn)定,占陽離子總量的18%~25%。陰離子分布差別顯著,PXS-1和PXS-3站點主要分布在HCO-Cl-線上,偏向HCO一端,SO所占比例較小,均值為14.6%;PXS-5站點主要分布在三角圖中央?yún)^(qū)域, 3種離子質(zhì)量濃度差別較小。對比淮河淮南段中下游枯水期數(shù)據(jù),其陽離子組成與3個研究站點較為接近,陰離子組成則和PXS-5站點類似,而淮南雨水中陰、陽離子以SO和Ca2+為主?;茨铣料莘e水區(qū)3個研究站點的水化學類型主要歸為Na+-HCO-Cl-型、Na+-Ca2+-HCO型。

    淮北3個湖泊站點的陽離子主要分布在(K++ Na+)-Mg2+線上(圖4),靠近(K++Na+)一端,Ca2+所占比例較小,均值為15.7%。其中HBDH與HBNH站點分布較為相似,(K++Na+)含量占陽離子總量的77.7%,而HBZH站點的(K++Na+)含量相對較小??傮w上,3個湖泊站點的水化學類型為Na+-Mg2+-HCO型,3條河流的水化學類型主要為Na+-Ca2+-HCO-Cl-型,Na+-HCO-SO型和Na+-Mg2+-HCO-SO型。

    圖3 淮南潘謝礦區(qū)沉陷積水區(qū)水體的主要離子組成三角圖Fig.3 Fraction ofmajor ions in the subsided waters from the Panxiemining area

    圖4 淮北朱-楊莊礦區(qū)沉陷積水區(qū)水體的主要離子組成三角圖Fig.4 Fraction ofmajor ions in the subsided waters from the Zhu-Yangzhuangmining area

    2.2 水化學特征影響因素

    2.2.1 Gibbs圖

    利用Gibbs的半對數(shù)坐標圖可以對影響水體離子特征的因素進行分析,其縱坐標以對數(shù)表示TDS質(zhì)量濃度,橫坐標以算術(shù)值表示ρ(Na+)/ρ(Na++ Ca2+)或ρ(Cl-)/ρ(Cl-+HCO)的比值[16]??梢詫⒌乇硭x子起源的自然影響因素劃分為3種機制:大氣降水、風化和蒸發(fā)-結(jié)晶作用[7]。一般認為,TDS質(zhì)量濃度很低(10 mg/L)且具有較大的ρ(Na+)/ρ(Na++ Ca2+)或ρ(Cl-)/ρ(Cl-+HCO)比值(接近于1)的水體,主要分布在圖的右下角,反映了海洋起源的大氣降水影響;TDS質(zhì)量濃度較高(70~300 mg/L)的水體且ρ(Na+)/ρ(Na++Ca2+)或ρ(Cl-)/ρ(Cl-+HCO)比值在0.5左右或小于0.5,主要分布在圖的左側(cè)中央,反映了巖石風化的影響;TDS質(zhì)量濃度最高的水體且具有較高的ρ(Na+)/ρ(Na++Ca2+)或ρ(Cl-)/ρ(Cl-+ HCO)比值(接近于1),與海水相似,主要分布在圖的右上角,體現(xiàn)了干旱地區(qū)蒸發(fā)作用的影響[17]。

    圖5表明潘謝沉陷積水區(qū)3個研究站點的TDS質(zhì)量濃度中等,ρ(Na+)/ρ(Na++Ca2+)的比值在0.55~0.81,ρ(Cl-)/ρ(Cl-+HCO)的比值在0.20~0.44,同時對比淮河河水和淮南雨水的分布特征,說明潘謝沉陷區(qū)水體離子特征主要受“巖石風化”及“蒸發(fā)-沉淀”作用的影響,而ρ(Cl-)/ρ(Cl-+HCO)值較ρ(Na+)/ρ(Na++Ca2+)小,在Gibbs圖上略偏向于左上方,可能受人為活動輸入的影響。

    圖5 淮南潘謝礦區(qū)沉陷積水區(qū)水化學的吉布斯分布模式Fig.5 Plots of themajor ionswithin the Gibbs boomerang envelope forwaters in the Panxiemining area

    圖6 淮北朱-楊莊礦區(qū)沉陷積水區(qū)水化學的吉布斯分布模式Fig.6 Plots of themajor ions within the Gibbs boomerang envelope for waters in the Zhu-Yangzhuangmining area

    淮北研究水體(圖6)中TDS質(zhì)量濃度較高, ρ(Na+)/ρ(Na++Ca2+)的比值亦較大(0.52~0.89), ρ(Cl-)/ρ(Cl-+HCO)的比值較小(0.16~0.42),且湖泊的比值普遍高于河流,體現(xiàn)了朱-楊莊沉陷區(qū)水體的封閉性,離子特征主要受“蒸發(fā)-沉淀”的作用,另一方面與區(qū)域水文地質(zhì)條件有關(guān)。另外可以看出,各個站點的離子組成大部分都落在了Gibbs圖中,只有朱-楊莊礦區(qū)的部分河流超出邊界,這是與Gibbs圖主要基于大型湖泊和河流繪制,小型水體以及干旱/半干旱地區(qū)小型湖泊人為干擾活動較大的因素有關(guān)[15,17]。

    2.2.2 離子來源判別

    對離子毫克當量濃度數(shù)據(jù)按研究站點進行KMO和Bartlett球形度檢驗,6個站點的KMO統(tǒng)計量均值為0.627,大于0.5,Bartlett球形度檢驗概率均為0,說明在99%的置信區(qū)間內(nèi)顯著相關(guān)。因此,本研究數(shù)據(jù)適合作因子分析。因子分析中抽取方法選擇主成分法,采用Varimax(最大正交旋轉(zhuǎn)法)對因子載荷進行旋轉(zhuǎn),目的在于使因子載荷矩陣結(jié)構(gòu)簡化,即每個變量在同一因子載荷的方差向最大和最小兩方面轉(zhuǎn)化。得到的旋轉(zhuǎn)載荷矩陣、特征值和貢獻率見表3。

    表3 兩淮采煤沉陷積水區(qū)水化學因子分析的旋轉(zhuǎn)載荷矩陣和貢獻率Table 3 Rotated com ponentmatrix and factor contribution to water solutes in term s of factor analysis in studied m ining areas

    PXS-1站點只提取出1個因子,貢獻率為91.09%,由于該積水區(qū)和泥河連通,常年受納周圍地區(qū)包括礦井水排放在內(nèi)的地表徑流,同時在不同季節(jié)和淺層地下水存在一定的水力聯(lián)系,各離子含量相關(guān)性強,因子分析很難區(qū)分各來源的具體貢獻;PXS-3站點提取出2個因子,由于其水體較為封閉,主要補給水源為淺層地下水,可將其視為第1因子,貢獻率為61.68%,第2因子與Cl-,K+,HCO相關(guān)性較大,貢獻率為24.62%,代表沉陷區(qū)周圍農(nóng)業(yè)活動影響; PXS-5站點同樣提取出2個因子,但各因子的貢獻與PXS-3站點存在一定差異,由于該站點周圍存在階段性調(diào)水的濟河,受納周圍及上游地區(qū)來水,將2個因子初步判斷為淺層地下水和人為活動輸入,其貢獻率相當。

    從淮北因子分析可以看出,HBDH只分離出1個因子,基于該站點已建設(shè)成為濕地公園,有良好的生態(tài)緩沖帶,周圍地表徑流輸入較小,其主要水源為地下水,各離子濃度在不同季節(jié)相關(guān)度較高(R均值為0.944)。HBZH和HBNH兩個站點較為相似,淺層地下水為主要的水源補給,但從因子分析中可以提取2個因子,分別與Ca2+,Cl-,F-和K+,Na+相關(guān)性較高,淮北地區(qū)地下水中F-含量相對較高,主要源于含氟礦石(CaF2)的溶解[18-19]。第1因子可能體現(xiàn)了這部分離子的來源,第2因子則體現(xiàn)了淮北地下水K+, Na+含量較高的特點。

    3 討 論

    水化學研究理論主要用于揭示世界各大河流溶質(zhì)含量與組成、來源,以及與流域地理氣候、水文地質(zhì)、生態(tài)環(huán)境和人為活動等方面的聯(lián)系。然而,水化學理論主要起源于大江、大河及其所對應(yīng)的大型流域范圍甚至全球尺度,和筆者研究的采煤沉陷區(qū)內(nèi)數(shù)十至數(shù)百平方公里的空間范圍相去甚遠[17],因此基于水化學一般原理辨別出小區(qū)域的特殊規(guī)律至關(guān)重要。由于本研究區(qū)域范圍較小(類似于小流域),沉陷區(qū)水體離子來源及組成主要反映了小區(qū)域水文地質(zhì)條件和生態(tài)環(huán)境的相關(guān)特征。從TDS含量來看,沉陷區(qū)水體遠高于我國地理氣候位置相似的主要大型河流和湖泊,如長江中下游[20]、淮河淮南段、鄱陽湖[21]等,屬于中等或較高礦化度水;從水化學三角圖中體現(xiàn)的離子組成來看,同一礦區(qū)小范圍內(nèi)亦存在一定的差異。

    淮南和淮北水文地質(zhì)條件、地表水系及土壤植被存在較大的差異?;茨系V區(qū)土壤黏土層較厚,具有隔水作用,積水區(qū)與淺層地下水可能主要通過采礦活動產(chǎn)生的裂隙進行溝通;淮北沉陷區(qū)干旱,土壤中砂質(zhì)成分多,導致其和淺層地下水的水力聯(lián)系更加直接?;茨系V區(qū)部分水體和周圍河流水系貫通(PXS-1)或發(fā)生水量交換(PXS-5),使得離子組成較為復雜,體現(xiàn)了其離子來源主要為淺層地下水和人為活動影響的地表水體;淮北礦區(qū)3個研究站點較為封閉,與周圍河流沒有直接的水力聯(lián)系,其補給水源主要為地下水,受外界影響相對較小。

    另一方面,沉陷積水區(qū)面積普遍較小且多為封閉,地下水及礦區(qū)地表徑流輸入的溶質(zhì)在“蒸發(fā)-沉淀”作用下,水體的礦化度將會逐漸升高,水化學組成發(fā)生改變,最終會導致水體咸化或堿化。特別是淮北沉陷積水區(qū),更加干旱缺水的氣候水文條件,導致其河流水系及沉陷區(qū)水體堿度顯著高于淮南地區(qū),并已經(jīng)開始影響水體功能。對此可以考慮在沉陷區(qū)設(shè)置節(jié)制閘,將其和周圍河流水系連通,在豐水期利用雨洪資源的沖淡作用。

    由于研究過程中缺乏周圍淺層地下水化學長期詳實的數(shù)據(jù),未來可以對沉陷區(qū)地表水體與周圍淺層地下水進行長期觀測,研究二者之間具體的水力聯(lián)系過程;由于影響水化學性質(zhì)的因素很多,各因素導致離子間的相互關(guān)系較為復雜,加之研究區(qū)域較小,下一步應(yīng)采取更為科學的定量化方法,如質(zhì)量平衡法和溶解物質(zhì)同位素示蹤法進行水體水源及離子來源判別。相關(guān)研究不僅可以為沉陷區(qū)水資源量預(yù)測及水化學變化提供基礎(chǔ)資料,而且可以為區(qū)域生態(tài)環(huán)境修復與重建等提供更加科學的依據(jù)。

    4 結(jié) 論

    (1)兩淮采煤沉陷區(qū)水體具有較高的礦化度,淮南潘謝礦區(qū)3個研究站點豐水期、平水期、枯水期TDS平均質(zhì)量濃度分別為306.6,408.4和544.9 mg/ L,淮北朱-楊莊礦區(qū)3個站點均值則為665.0,629.3和861.1 mg/L,枯水期質(zhì)量濃度高于豐水期,淮北沉陷積水區(qū)水體礦化度、堿度和總硬度均高于淮南沉陷水體。

    (2)沉陷區(qū)水體化學組成在礦區(qū)小范圍內(nèi)具有一定的空間變異性,HCO和Na+分別在陰陽離子中占有最大比例。ρ(Na+)/ρ(Na++Ca2+)比值遠高于ρ(Cl-)/ρ(Cl-+HCO),礦區(qū)水體中較高的HCO含量是導致以上兩組比值出現(xiàn)不平衡的主要原因?;茨铣料莘e水區(qū)的水化學類型主要為Na+-HCO-Cl-, Na+-Ca2+-HCO型,淮北3個湖泊站點的水化學類型則為Na+-Mg2+-HCO型。

    (3)Gibbs圖和因子分析結(jié)果均揭示了兩淮礦區(qū)水文地質(zhì)條件上的差異,淮南沉陷積水區(qū)由于水體補給來源的多樣,離子組成較為復雜,淺層地下水和地表徑流為水化學離子組成的主要影響因素;淮北沉陷積水區(qū)的研究站點較封閉,補給水源為地下水,離子組成主要體現(xiàn)了區(qū)域淺層地下水化學的相關(guān)特征。

    [1] 王婷婷,易齊濤,胡友彪,等.兩淮采煤沉陷區(qū)水域水體富營養(yǎng)化及氮、磷限制模擬實驗[J].湖泊科學,2013,25(6):916-926.

    Wang Tingting,Yi Qitao,Hu Youbiao,et al.Eutrophication and nutrientenrichmentbioassays in thewatersof the Huainan and Huaibei coal mining subsidence areas,Anhui Province[J].Journal of Lake Sciences,2013,25(6):916-926.

    [2] 李思悅,譚 香,徐志方,等.湖北丹江口水庫主要離子化學季節(jié)變化及離子來源分析[J].環(huán)境科學,2008,29(12):3353-3359.

    Li Siyue,Tan Xiang,Xu Zhifang,et al.Seasonal variation in themajor ion chemistry and their sources in the HubeiDanjiangkou Reservoir,China[J].Environmental Science,2008,29(12):3353-3359.

    [3] Timofeeff M N,Lowenstein T K,Silva M A M,et al.Secular variation in the major-ion chemistry of seawater:Evidence from fluid inclusions in Cretaceous Halites[J].Geochimica et Cosmochimica Acta, 2006,70:1977-1994.

    [4] Brennan SK,Lowenstein TK.Themajor ion composition of Silurian seawater[J].Geochimica et Cosmochimica Acta,2002,66:2683-2700.

    [5] 夏學齊,楊忠芳,王亞平,等.長江水系河水主要離子化學特征[J].地學前緣,2008,15(5):194-202.

    Xia Xueqi,Yang Zhongfang,Wang Yaping,etal.Major ion chemistry in the Yangtze River[J].Earth Science Frontiers,2008,15(5): 194-202.

    [6] 陳靜生,王飛越,何大偉.黃河水質(zhì)地球化學[J].地學前緣, 2006,13(1):58-73.

    Chen Jingsheng,Wang Feiyue,He Dawei.Geochemistry of water quality of the Yellow River basin[J].Earth Science Frontiers,2006, 13(1):58-73.

    [7] 葉宏萌,袁旭音,葛敏霞,等.太湖北部流域水化學特征及其控制因素[J].生態(tài)環(huán)境學報,2010,19(1):23-27.

    Ye Hongmeng,Yuan Xuyin,Ge Minxia,et al.Water chemistry characteristics and controlling factors in the northern rivers in the Taihu basin[J].Ecology and Environmental Sciences,2010,19(1):23-27.

    [8] 何春桂,劉 輝,桂和榮.淮南市典型采煤塌陷區(qū)水域環(huán)境現(xiàn)狀評價[J].煤炭學報,2005,30(6):754-758.

    He Chungui,Liu Hui,GuiHerong.Environmentevaluation on typical water area resulting from coal-mine subsidence in Huainan[J].Journal of China Coal Society,2005,30(6):754-758.

    [9] 曲喜杰,易齊濤,胡友彪,等.兩淮采煤沉陷積水區(qū)水體營養(yǎng)鹽時空分布及富營養(yǎng)化進程[J].應(yīng)用生態(tài)學報,2013,24(11): 3249-3258.

    Qu Xijie,Yi Qitao,Hu Youbiao,et al.Nutrient spatiotemporal distribution and eutrophication process in subsided waters of Huainan and Huaibeimining areas,China[J].Chinese Journal of Applied Ecology,2013,24(11):3249-3258.

    [10] 桂和榮,王和平,方文慧,等.煤礦塌陷區(qū)水域環(huán)境指示微生物——藍藻的研究[J].煤炭學報,2007,32(8):848-853.

    Gui Herong,Wang Heping,Fang Wenhui,et al.Study of the cyanophyta which is one of the environment indicator microorganism in mining subsidencewater area[J].Journalof China Coal Society, 2007,32(8):848-853.

    [11] 謝 凱,易齊濤,孫鵬飛.采煤沉陷積水區(qū)土壤覆水初期氧化還原反應(yīng)過程與磷遷移轉(zhuǎn)化的耦合研究[J].環(huán)境科學學報, 2013,33(11):3101-3110.

    Xie Kai,Yi Qitao,Sun Pengfei.The early processes of redox and phosphorusmigration through simulated experiments of inundated soil column in the subsidence waters of coal mine areas[J].Acta Scientiae Circumstantiae,2013,33(11):3101-3110.

    [12] 易齊濤,孫鵬飛,謝 凱,等.區(qū)域水化學條件對淮南采煤沉陷區(qū)水域沉積物磷吸附特征的影響研究[J].環(huán)境科學,2013, 34(10):3894-3903.

    Yi Qitao,Sun Pengfei,Xie Kai,et al.Impact of regional water chemistry on the phosphorus isothermal adsorption of the sediments in three subsidencewaters of the Huainan Mine Areas[J].Environmental Science,2013,34(10):3894-3903.

    [13] 金相燦,屠清瑛.湖泊富營養(yǎng)化調(diào)查規(guī)范(第二版)[M].北京:中國環(huán)境科學出版社,1990.

    [14] 林雨霏,劉素美,紀 雷,等.黃海春夏季降水中常量離子化學特征分析[J].海洋環(huán)境科學,2007,26(2):116-120.

    Lin Yufei,Liu Sumei,Ji Lei,et al.Chemical characteristic analysis on major ions of rainwater at Qianliyan island in spring and summer [J].Marine Environmental Science,2007,26(2):116-120.

    [15] 侯昭華,徐 海,安芷生.青海湖流域水化學主離子特征及控制因素初探[J].地球與環(huán)境,2009,37(1):11-19.

    Hou Zhaohua,Xu Hai,An Zhisheng.Major ion chemistry of waters in lake Qinghai catchment and the possible controls[J].Earth and Environment,2009,37(1):11-19.

    [16] 陳靜生,何大偉.珠江水系河水主要離子化學特征及成因[J].北京大學學報(自然科學版),1999,35(6):786-793.

    Chen Jingsheng,He Dawei.Chemical characteristics and genesis ofmajor ions in the Pearl River basin[J].Acta Scientiarum Naturalium Universitatis Pekinensis,1999,35(6):786-793.

    [17] Kilham P.Mechanisms controlling the chemical composition of lakes and rivers:Data from Africa[J].Limnology and Oceanography,1990,35(1):80-831.

    [18] 丁 丹,許光泉,何曉文,等.淮北平原淺層地下水氟的水化學特征及影響因素分析[J].水資源保護,2009,25(2):64-68.

    Ding Dan,Xu Guangquan,He Xiaowen,etal.Chemical characteristics and influential factors of fluorine ions in shallow groundwater of Huaibei plain[J].Water Resources Protection,2009,25(2):64-68.

    [19] 吳泊人,錢家忠.安徽省淮北平原地下水氟、溶解性總固體空間變異特征[J].地理科學,2011,31(4):453-458.

    Wu Boren,Qian Jiazhong.Spatial variation characteristics of fluorine and total dissolved solid in groundwater in Huaibei plain,Anhui Province[J].Scientia Geographica Sinica,2011,31(4):453-458.

    [20] 陳靜生.河流水質(zhì)原理及中國河流水質(zhì)[M].北京:科學出版社,2007.

    [21] 胡春華,周文斌,夏思奇.鄱陽湖流域水化學主離子特征及其來源分析[J].環(huán)境化學,2011,30(9):1620-1626.

    Hu Chunhua,Zhou Wenbin,Xia Siqi.Characteristics ofmajor ions and the influence factors in Poyang Lake catchment[J].Environmental Chemistry,2011,30(9):1620-1626.

    Characteristics of water chem istry and their influencing factors in subsidence waters in the Huainan and Huaibeim ining areas,Anhui Province

    SUN Peng-fei,YIQi-tao,XU Guang-quan

    (School of Earth and Environment,Anhui University of Science and Technology,Huainan 232001,China)

    :In order to know water sources and total dissolved solids(TDS)ofmining subsidencewaters,six siteswere selected in the mining areas of Huainan“Panxie”and Huaibei“Zhu-Yang Zhuang”.The contents,compositions and types ofmajor ions were analyzed based on the general theories and principles of water chemistry.Using the Gibbs graph and themethod of factor analysis,the influencing factors were further discussed.The results indicate that the TDS are elevated in all sites during the period ofwet season,median water period,and dry season.TDS,alkalinity and hardness of the subsidencewaters in Huainbeiare obviously higher than those in Huainan.Meanwhile,Na+and HCOaccounts for a great part of the total cations and anions respectively,although their water chemistry varied depending on the sampling site.The Huainan mining subsidence watersmainly belongs to water chemistry type of“Na+-HCO-Cl-”and“Na+-Ca2+-HCO”,and it is the type of“Na+-Mg2+-HCO”in the subsidence waters of Huaibei.Over all,the Gibbs graph and factor analysis show thatwater chemistry in the Huainan mining subsidencewaters aremainly related to local shallow underground water and surface runoff,whereas it reflects the relationship to underground water for the compositions of ions in the Huaibeimining subsidence waters.

    孫鵬飛,易齊濤,許光泉.兩淮采煤沉陷積水區(qū)水體水化學特征及影響因素[J].煤炭學報,2014,39(7):1345-1353.

    10.13225/j.cnki.jccs.2013.0959

    Sun Pengfei,Yi Qitao,Xu Guangquan.Characteristics ofwater chemistry and their influencing factors in subsidence waters in the Huainan and Huaibeim ining areas,Anhui Province[J].Journal of China Coal Society,2014,39(7):1345-1353.doi:10.13225/j.cnki.jccs.2013.0959

    2013-07-08 責任編輯:韓晉平

    國家自然科學基金資助項目(41202242);教育部留學回國人員科研啟動基金資助項目(第45批)

    孫鵬飛(1991—),男,安徽安慶人,碩士研究生。E-mail:superfell@126.com。通訊作者:易齊濤(1979—),男,湖北天門人,副教授,博士。E-mail:yiqitao@163.com

    P592;X142

    A

    0253-9993(2014)07-1345-09

    猜你喜歡
    水化學淮北淮南
    幸福湖泉域地下水水化學特征及其演化機制
    人民黃河(2025年1期)2025-01-16 00:00:00
    《淮南師范學院學報》投稿須知
    《淮北師范大學學報》(自然科學版)征稿簡則
    秦皇島北部地表水和地下水水化學特征及水質(zhì)評價
    《淮北師范大學學報》(自然科學版)征稿簡則
    《淮北枳》
    學生天地(2016年10期)2016-04-16 05:14:49
    衡水市新近系館陶組地熱水水化學特征與分布規(guī)律研究
    淮北 去產(chǎn)能的黑色面孔
    能源(2016年10期)2016-02-28 11:33:25
    CRADLE OF TOFU BY DAVID dawson
    筠連煤層氣井產(chǎn)出水化學特征及意義
    中國煤層氣(2015年5期)2015-08-22 03:26:32
    a级一级毛片免费在线观看| 亚洲色图av天堂| 日韩欧美国产在线观看| 乱码一卡2卡4卡精品| 成年av动漫网址| 久久精品国产亚洲av天美| 一卡2卡三卡四卡精品乱码亚洲| 成人亚洲精品av一区二区| 亚洲av不卡在线观看| 亚洲人成网站在线播| 久久精品国产99精品国产亚洲性色| 国产69精品久久久久777片| 欧美3d第一页| 国产又黄又爽又无遮挡在线| 赤兔流量卡办理| 97超碰精品成人国产| 久久亚洲国产成人精品v| av在线亚洲专区| 伦精品一区二区三区| 国产精品99久久久久久久久| 欧美区成人在线视频| 你懂的网址亚洲精品在线观看 | 一个人看视频在线观看www免费| 国产精品不卡视频一区二区| 熟妇人妻久久中文字幕3abv| 日本成人三级电影网站| 97超视频在线观看视频| 国产淫片久久久久久久久| 国产又黄又爽又无遮挡在线| 国产精品一区二区性色av| 国产成人91sexporn| 夜夜爽天天搞| 国产亚洲av片在线观看秒播厂 | 午夜福利成人在线免费观看| 亚洲av二区三区四区| 看非洲黑人一级黄片| 日本av手机在线免费观看| 青春草国产在线视频 | 久久久久久久久久黄片| 男人和女人高潮做爰伦理| 亚洲国产精品国产精品| 天堂中文最新版在线下载 | 婷婷色av中文字幕| 黄色一级大片看看| 久久久久久伊人网av| 精品久久久久久久末码| 亚洲在线自拍视频| 日韩精品青青久久久久久| 国产一区二区激情短视频| 亚洲国产精品成人久久小说 | 99九九线精品视频在线观看视频| 国产精品久久久久久久久免| 国产精品一区二区在线观看99 | 人人妻人人澡人人爽人人夜夜 | 在线观看午夜福利视频| 真实男女啪啪啪动态图| 卡戴珊不雅视频在线播放| 久久6这里有精品| 亚洲最大成人手机在线| 久久草成人影院| 国产精品美女特级片免费视频播放器| 国产麻豆成人av免费视频| 亚洲不卡免费看| 只有这里有精品99| 成人漫画全彩无遮挡| 麻豆精品久久久久久蜜桃| 在线播放无遮挡| 亚洲精品久久久久久婷婷小说 | 99久久精品热视频| 精品欧美国产一区二区三| 久久久精品大字幕| 三级国产精品欧美在线观看| 五月玫瑰六月丁香| 久久精品国产亚洲网站| 色5月婷婷丁香| 免费看a级黄色片| 18禁在线播放成人免费| 国产精品蜜桃在线观看 | 亚洲精品国产成人久久av| 久久精品夜色国产| 日韩,欧美,国产一区二区三区 | 日韩av在线大香蕉| 别揉我奶头 嗯啊视频| 中国美白少妇内射xxxbb| 久久人人精品亚洲av| 观看美女的网站| 秋霞在线观看毛片| 亚洲最大成人中文| 国产不卡一卡二| 男女视频在线观看网站免费| 99久久人妻综合| 成人亚洲欧美一区二区av| 亚洲av.av天堂| 最近手机中文字幕大全| av专区在线播放| 男女啪啪激烈高潮av片| 嫩草影院新地址| 欧美成人免费av一区二区三区| www.av在线官网国产| 成人国产麻豆网| 好男人视频免费观看在线| 午夜精品在线福利| 哪里可以看免费的av片| av在线蜜桃| 精品熟女少妇av免费看| 国内少妇人妻偷人精品xxx网站| 老熟妇乱子伦视频在线观看| 久久鲁丝午夜福利片| 12—13女人毛片做爰片一| 91麻豆精品激情在线观看国产| 美女cb高潮喷水在线观看| 黄色欧美视频在线观看| 此物有八面人人有两片| 日本免费a在线| 校园春色视频在线观看| 人妻夜夜爽99麻豆av| 长腿黑丝高跟| 午夜激情福利司机影院| 人妻制服诱惑在线中文字幕| 国产私拍福利视频在线观看| 人妻系列 视频| 天天一区二区日本电影三级| 国产伦精品一区二区三区视频9| 好男人视频免费观看在线| 男女那种视频在线观看| 国内精品美女久久久久久| 亚洲国产欧美人成| 亚洲欧美精品综合久久99| 99久久九九国产精品国产免费| 久久婷婷人人爽人人干人人爱| 春色校园在线视频观看| 亚洲精华国产精华液的使用体验 | 亚洲国产色片| 国产免费男女视频| 黄片无遮挡物在线观看| 亚洲欧美清纯卡通| 一级毛片aaaaaa免费看小| 亚洲精品亚洲一区二区| а√天堂www在线а√下载| 欧洲精品卡2卡3卡4卡5卡区| 99久久久亚洲精品蜜臀av| 久久久精品欧美日韩精品| 又爽又黄无遮挡网站| 国产精品国产高清国产av| 热99在线观看视频| 三级毛片av免费| 成人永久免费在线观看视频| 99久久中文字幕三级久久日本| 亚洲国产色片| 亚洲无线观看免费| 高清日韩中文字幕在线| 三级毛片av免费| 免费看光身美女| 婷婷六月久久综合丁香| 亚洲av免费高清在线观看| 国产成人午夜福利电影在线观看| 亚洲欧美成人综合另类久久久 | 国产在线精品亚洲第一网站| 国产乱人视频| 日韩av在线大香蕉| 免费观看精品视频网站| avwww免费| 亚洲精品国产av成人精品| 可以在线观看毛片的网站| 亚洲真实伦在线观看| 成人午夜高清在线视频| 精品人妻一区二区三区麻豆| 在线观看av片永久免费下载| 日本av手机在线免费观看| 亚洲久久久久久中文字幕| 亚洲电影在线观看av| 欧美一区二区国产精品久久精品| 高清毛片免费看| 1000部很黄的大片| 欧美激情在线99| 国产亚洲av片在线观看秒播厂 | 国产精品永久免费网站| 此物有八面人人有两片| 六月丁香七月| 国产伦在线观看视频一区| 久久精品夜夜夜夜夜久久蜜豆| 春色校园在线视频观看| 精品少妇黑人巨大在线播放 | 男插女下体视频免费在线播放| 可以在线观看毛片的网站| 国产午夜精品久久久久久一区二区三区| 99久久无色码亚洲精品果冻| 日韩在线高清观看一区二区三区| 久久精品久久久久久噜噜老黄 | av天堂在线播放| 小蜜桃在线观看免费完整版高清| av黄色大香蕉| 国产日韩欧美在线精品| 国产探花极品一区二区| 边亲边吃奶的免费视频| 日本免费一区二区三区高清不卡| 搡女人真爽免费视频火全软件| 偷拍熟女少妇极品色| 一卡2卡三卡四卡精品乱码亚洲| 日本黄色视频三级网站网址| 国产黄片视频在线免费观看| 美女被艹到高潮喷水动态| av国产免费在线观看| 国产精华一区二区三区| 久久鲁丝午夜福利片| 爱豆传媒免费全集在线观看| 欧美另类亚洲清纯唯美| 神马国产精品三级电影在线观看| 高清毛片免费观看视频网站| 中文资源天堂在线| 精品国产三级普通话版| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区在线观看99 | 国产探花极品一区二区| 五月玫瑰六月丁香| 日本在线视频免费播放| 2022亚洲国产成人精品| 99久久精品国产国产毛片| 日日摸夜夜添夜夜添av毛片| 国产精品久久久久久精品电影| 熟女人妻精品中文字幕| 久久精品91蜜桃| 国产精品精品国产色婷婷| av黄色大香蕉| 亚洲真实伦在线观看| 97热精品久久久久久| 波多野结衣巨乳人妻| 日韩精品有码人妻一区| 最近2019中文字幕mv第一页| 欧美日韩一区二区视频在线观看视频在线 | 两个人视频免费观看高清| 久久久国产成人精品二区| 国产视频内射| 色噜噜av男人的天堂激情| 国产精品女同一区二区软件| 国产日韩欧美在线精品| 搡女人真爽免费视频火全软件| 51国产日韩欧美| h日本视频在线播放| 有码 亚洲区| 亚洲人成网站在线播| ponron亚洲| 丝袜美腿在线中文| 中国美白少妇内射xxxbb| 国产精品一及| 欧美另类亚洲清纯唯美| 亚洲激情五月婷婷啪啪| 午夜精品在线福利| 成人高潮视频无遮挡免费网站| 国产精品人妻久久久影院| 成人二区视频| 人妻夜夜爽99麻豆av| 91久久精品国产一区二区成人| 欧美潮喷喷水| 亚洲18禁久久av| 蜜桃亚洲精品一区二区三区| 国产精品久久久久久亚洲av鲁大| 国产真实乱freesex| 99久久九九国产精品国产免费| 成年女人永久免费观看视频| 啦啦啦观看免费观看视频高清| 国产一区二区在线观看日韩| 欧美色欧美亚洲另类二区| 午夜福利在线观看免费完整高清在 | 少妇人妻精品综合一区二区 | 91久久精品国产一区二区三区| 久久精品国产鲁丝片午夜精品| 亚洲精品自拍成人| 国产精品一区www在线观看| h日本视频在线播放| 99九九线精品视频在线观看视频| 免费人成在线观看视频色| 日韩精品有码人妻一区| 国产成人aa在线观看| 成人亚洲欧美一区二区av| 观看美女的网站| 亚洲三级黄色毛片| av.在线天堂| 免费av毛片视频| 国产黄a三级三级三级人| 国产精品嫩草影院av在线观看| 国产成人a区在线观看| 午夜福利在线观看免费完整高清在 | 99riav亚洲国产免费| 又黄又爽又刺激的免费视频.| 寂寞人妻少妇视频99o| 精品久久久久久久久av| 一进一出抽搐动态| 在线播放无遮挡| 国产真实伦视频高清在线观看| 亚洲av熟女| 直男gayav资源| 你懂的网址亚洲精品在线观看 | 久久精品综合一区二区三区| 国产亚洲av嫩草精品影院| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久久电影| 亚洲18禁久久av| 人妻少妇偷人精品九色| 久久久久免费精品人妻一区二区| 亚洲自拍偷在线| 我要搜黄色片| 欧美成人精品欧美一级黄| 九草在线视频观看| 草草在线视频免费看| 中国美女看黄片| av在线老鸭窝| 日韩欧美 国产精品| 欧美性猛交╳xxx乱大交人| 人妻系列 视频| 中文资源天堂在线| 男女边吃奶边做爰视频| 亚洲在线观看片| 欧美bdsm另类| 成人二区视频| 国产极品精品免费视频能看的| 国产一区二区亚洲精品在线观看| 99视频精品全部免费 在线| 在线观看66精品国产| 黄片无遮挡物在线观看| 国产精品女同一区二区软件| 99热这里只有是精品在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 欧美bdsm另类| 国模一区二区三区四区视频| 中文字幕人妻熟人妻熟丝袜美| 不卡一级毛片| 久久久久久大精品| 1000部很黄的大片| 偷拍熟女少妇极品色| 国产91av在线免费观看| 在线a可以看的网站| 精品午夜福利在线看| 亚洲成人久久爱视频| 麻豆成人av视频| 亚洲精品成人久久久久久| 国产亚洲av嫩草精品影院| 色综合色国产| 久久综合国产亚洲精品| 亚洲无线观看免费| 欧美性猛交╳xxx乱大交人| 午夜免费男女啪啪视频观看| 亚洲18禁久久av| av在线老鸭窝| 久久精品国产清高在天天线| 观看免费一级毛片| 美女黄网站色视频| 亚洲在线观看片| 色吧在线观看| 丝袜喷水一区| 中文精品一卡2卡3卡4更新| 中文在线观看免费www的网站| avwww免费| 午夜福利成人在线免费观看| 看非洲黑人一级黄片| 国产一级毛片在线| 18禁黄网站禁片免费观看直播| 麻豆精品久久久久久蜜桃| 日本免费一区二区三区高清不卡| 中文字幕精品亚洲无线码一区| 边亲边吃奶的免费视频| 国产真实乱freesex| 1000部很黄的大片| 成人av在线播放网站| 永久网站在线| 精品欧美国产一区二区三| 国产极品精品免费视频能看的| 一级毛片我不卡| 变态另类成人亚洲欧美熟女| 男插女下体视频免费在线播放| av天堂中文字幕网| 国产成人影院久久av| 中文亚洲av片在线观看爽| 成人无遮挡网站| 大香蕉久久网| 亚洲av熟女| 国产毛片a区久久久久| 精品免费久久久久久久清纯| 亚洲国产精品合色在线| 欧美在线一区亚洲| 国产精品一区二区在线观看99 | 国产成人精品久久久久久| 日本一二三区视频观看| 亚洲一区二区三区色噜噜| 精品欧美国产一区二区三| 精品熟女少妇av免费看| 精品欧美国产一区二区三| 亚洲欧美精品自产自拍| 国产淫片久久久久久久久| 麻豆久久精品国产亚洲av| 久久九九热精品免费| 亚洲人成网站在线播放欧美日韩| 亚洲精品久久国产高清桃花| 熟妇人妻久久中文字幕3abv| 有码 亚洲区| 亚洲人成网站高清观看| 精品午夜福利在线看| 国产色爽女视频免费观看| 97超视频在线观看视频| 欧美最黄视频在线播放免费| 久久久精品欧美日韩精品| 18+在线观看网站| 一级毛片电影观看 | 欧美人与善性xxx| 在线观看美女被高潮喷水网站| 国产精品乱码一区二三区的特点| 欧美在线一区亚洲| 黄色日韩在线| 春色校园在线视频观看| 12—13女人毛片做爰片一| 三级毛片av免费| 男女啪啪激烈高潮av片| 三级毛片av免费| 欧美+亚洲+日韩+国产| 一级黄片播放器| 亚洲成人精品中文字幕电影| 一级av片app| 久久国产乱子免费精品| 九九爱精品视频在线观看| 国产不卡一卡二| 国产久久久一区二区三区| 国产私拍福利视频在线观看| 成人综合一区亚洲| 国产视频首页在线观看| 国产精品99久久久久久久久| 99久久无色码亚洲精品果冻| 国产午夜福利久久久久久| 精品一区二区三区视频在线| 午夜福利在线在线| 成人亚洲精品av一区二区| 国产高清激情床上av| 国产在视频线在精品| 亚洲av成人av| 免费一级毛片在线播放高清视频| av在线天堂中文字幕| 亚洲国产精品合色在线| 亚洲高清免费不卡视频| 永久网站在线| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品sss在线观看| 亚洲av电影不卡..在线观看| 日韩欧美三级三区| 国产激情偷乱视频一区二区| 乱人视频在线观看| h日本视频在线播放| 91久久精品国产一区二区三区| 欧美最黄视频在线播放免费| 日产精品乱码卡一卡2卡三| 亚洲人成网站在线播放欧美日韩| 成人性生交大片免费视频hd| 久久亚洲国产成人精品v| 国产一区亚洲一区在线观看| 成年av动漫网址| 久久热精品热| 免费观看的影片在线观看| 黄色配什么色好看| 丰满乱子伦码专区| 国产av一区在线观看免费| 亚洲乱码一区二区免费版| 老司机影院成人| 免费av不卡在线播放| 在线观看美女被高潮喷水网站| 大又大粗又爽又黄少妇毛片口| 久久国产乱子免费精品| 亚洲最大成人手机在线| 99久久九九国产精品国产免费| 九色成人免费人妻av| 高清在线视频一区二区三区 | 91午夜精品亚洲一区二区三区| 韩国av在线不卡| 哪个播放器可以免费观看大片| 亚洲人成网站在线播| 久久久久久伊人网av| 国产日韩欧美在线精品| 97在线视频观看| 联通29元200g的流量卡| 久久亚洲精品不卡| 美女国产视频在线观看| 三级毛片av免费| 网址你懂的国产日韩在线| 国产精品精品国产色婷婷| 成人午夜精彩视频在线观看| 麻豆国产97在线/欧美| 校园春色视频在线观看| 亚洲欧美精品专区久久| 精品人妻熟女av久视频| 成人特级黄色片久久久久久久| 好男人视频免费观看在线| 国内精品一区二区在线观看| 成人漫画全彩无遮挡| 赤兔流量卡办理| 国产亚洲91精品色在线| 高清在线视频一区二区三区 | 亚洲精华国产精华液的使用体验 | 国产精品电影一区二区三区| 丰满人妻一区二区三区视频av| 国内久久婷婷六月综合欲色啪| eeuss影院久久| 综合色av麻豆| 亚洲av中文字字幕乱码综合| 校园人妻丝袜中文字幕| 高清午夜精品一区二区三区 | 亚洲成人精品中文字幕电影| 高清在线视频一区二区三区 | 日本爱情动作片www.在线观看| 一级黄色大片毛片| 最近2019中文字幕mv第一页| 成人国产麻豆网| 国产色婷婷99| 欧美性猛交╳xxx乱大交人| 看免费成人av毛片| 日本欧美国产在线视频| 深夜a级毛片| 国产精品永久免费网站| 亚洲熟妇中文字幕五十中出| 看黄色毛片网站| 丝袜美腿在线中文| 1000部很黄的大片| 国产精品伦人一区二区| 久久中文看片网| 18+在线观看网站| 亚洲高清免费不卡视频| 精品人妻偷拍中文字幕| 一卡2卡三卡四卡精品乱码亚洲| av免费在线看不卡| 97超碰精品成人国产| 欧美成人一区二区免费高清观看| 91午夜精品亚洲一区二区三区| 成熟少妇高潮喷水视频| 天堂√8在线中文| 人妻制服诱惑在线中文字幕| 最近手机中文字幕大全| 亚洲精品自拍成人| 日本三级黄在线观看| 久久久久久久午夜电影| 久久久久久久久久久免费av| 午夜爱爱视频在线播放| 黄片无遮挡物在线观看| av黄色大香蕉| 久久久精品欧美日韩精品| 在线a可以看的网站| 国产人妻一区二区三区在| 深夜a级毛片| 久久久久久久亚洲中文字幕| 嫩草影院新地址| 高清午夜精品一区二区三区 | 26uuu在线亚洲综合色| 国产综合懂色| 偷拍熟女少妇极品色| 天堂中文最新版在线下载 | 日韩人妻高清精品专区| 女的被弄到高潮叫床怎么办| 深夜精品福利| 一个人看视频在线观看www免费| 男人和女人高潮做爰伦理| 亚洲乱码一区二区免费版| 国产乱人偷精品视频| 日韩强制内射视频| 全区人妻精品视频| 啦啦啦啦在线视频资源| 色吧在线观看| 少妇的逼水好多| 亚洲高清免费不卡视频| 国产亚洲5aaaaa淫片| 亚洲天堂国产精品一区在线| 色综合站精品国产| 国产精品一区二区在线观看99 | 亚洲精品粉嫩美女一区| 不卡视频在线观看欧美| 久久久国产成人免费| 亚洲欧洲国产日韩| 国产精品野战在线观看| 91久久精品国产一区二区成人| 国产爱豆传媒在线观看| 高清在线视频一区二区三区 | 夜夜爽天天搞| ponron亚洲| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品粉嫩美女一区| 久久久久久久久久成人| 久久韩国三级中文字幕| 熟女电影av网| 欧美色欧美亚洲另类二区| 国产成人91sexporn| 在线国产一区二区在线| 好男人视频免费观看在线| 欧美又色又爽又黄视频| 国产精品永久免费网站| 国产高清视频在线观看网站| 亚洲人成网站在线观看播放| 夜夜爽天天搞| 久久久午夜欧美精品| 成人亚洲精品av一区二区| 国产蜜桃级精品一区二区三区| 国产 一区 欧美 日韩| 日韩强制内射视频| 国产精品久久久久久久电影| 女的被弄到高潮叫床怎么办| 国产国拍精品亚洲av在线观看| 看免费成人av毛片| 国产黄色视频一区二区在线观看 | 欧美成人a在线观看| ponron亚洲| 久久精品国产鲁丝片午夜精品| 在线观看美女被高潮喷水网站| 麻豆久久精品国产亚洲av| 欧美成人免费av一区二区三区| 哪里可以看免费的av片| 婷婷精品国产亚洲av| 亚洲一区高清亚洲精品| 国产精品人妻久久久影院| www.av在线官网国产| 午夜激情欧美在线| 日韩三级伦理在线观看| 岛国毛片在线播放|