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

    采煤驅(qū)動(dòng)下復(fù)雜井田含水層化學(xué)特征與水力聯(lián)系辨識(shí)

    2022-06-29 08:53:40侯澤明張圣微2李鋼柱劉志強(qiáng)
    中國(guó)環(huán)境科學(xué) 2022年6期
    關(guān)鍵詞:直羅第四系圖版

    黃 磊,侯澤明,韓 萱,許 磊,張圣微,2李鋼柱,劉志強(qiáng)

    采煤驅(qū)動(dòng)下復(fù)雜井田含水層化學(xué)特征與水力聯(lián)系辨識(shí)

    黃 磊1*,侯澤明1,韓 萱1,許 磊1,張圣微1,2李鋼柱1,劉志強(qiáng)1

    (1.內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院,內(nèi)蒙古自治區(qū)水資源保護(hù)與利用重點(diǎn)實(shí)驗(yàn)室,內(nèi)蒙古 呼和浩特 010018;2.內(nèi)蒙古自治區(qū)農(nóng)牧業(yè)大數(shù)據(jù)研究與應(yīng)用重點(diǎn)實(shí)驗(yàn)室,內(nèi)蒙古 呼和浩特 010018)

    通過(guò)分析復(fù)雜井田不同含水層的常規(guī)元素、微量元素、氘氧同位素與氚同位素的水化學(xué)特征來(lái)判斷各含水層的水力聯(lián)系,并分別建立了Piper識(shí)別圖版、Durov識(shí)別圖版、氘氧同位素識(shí)別圖版,借此可以快速甄別礦井突水的來(lái)源.結(jié)果表明:研究區(qū)地下水主要為大氣降水補(bǔ)給,第四系含水層與直羅組含水層存在顯著聯(lián)系,直羅組含水層與延安組存在有限的聯(lián)系.利用各含水層不同巖性導(dǎo)致的水化學(xué)離子特征差異和氘氧同位素、氚同位素的示蹤特性建立識(shí)別圖版,可有助于快速識(shí)別補(bǔ)連塔礦區(qū)突水水源,并對(duì)不同含水層的突水事故提出針對(duì)性的解決措施.

    復(fù)雜井田;多層含水層;水力聯(lián)系;礦井突水;水化學(xué)特征;同位素示蹤

    當(dāng)前,國(guó)內(nèi)外學(xué)者對(duì)采煤驅(qū)動(dòng)影響地下水系統(tǒng)變化機(jī)制[1-2]進(jìn)行了一定研究:如錢鳴高等[3]通過(guò)研究采煤驅(qū)動(dòng)作用對(duì)巖層原有平衡狀態(tài)的破壞,認(rèn)為采煤活動(dòng)導(dǎo)致裂隙增加、地下水水位下降與地表塌陷;冀瑞君等[4]通過(guò)研究神東礦區(qū)采煤活動(dòng)對(duì)窟野河流域地下水循環(huán)的影響機(jī)制,認(rèn)為采動(dòng)破壞了隔水層的隔水性,使含水盆地匯水面積減小,泉流量減小或干涸,同時(shí)一部分地下水轉(zhuǎn)化為礦井水;范立民等[5]通過(guò)對(duì)比1994年與2015年毛烏素沙漠與黃土高原接壤區(qū)泉數(shù)量、流量變化趨勢(shì)等,認(rèn)為混合入滲補(bǔ)給點(diǎn)受隔水層結(jié)構(gòu)破壞與側(cè)向補(bǔ)給截?cái)嗍菍?dǎo)致近20a來(lái)泉點(diǎn)大量消失的主要原因.

    鄂爾多斯市伊金霍洛旗境內(nèi)有70多座煤礦,其中大規(guī)模大產(chǎn)能煤礦主要分布在黃河的一級(jí)支流烏蘭木倫河兩岸,是典型的沿河礦區(qū).連年持續(xù)高強(qiáng)度開采導(dǎo)致地下水問(wèn)題日益嚴(yán)重,已經(jīng)成為制約當(dāng)?shù)厣鐣?huì)發(fā)展的重要因素.本文選取伊金霍洛旗補(bǔ)連塔煤礦為研究區(qū),該煤礦坐落于烏蘭木倫河西岸,是世界第一大單井井工礦.采煤驅(qū)動(dòng)作用直接影響著烏蘭木倫河流域引發(fā)地下水水位下降、淺水含水層疏干等一系列環(huán)境地質(zhì)問(wèn)題[6].

    本研究區(qū)范圍內(nèi)前人的研究多集中于烏蘭木倫河流域地下水化學(xué)特征與成因機(jī)制分析[7]或某一離子的空間分布變化規(guī)律[8].多手段、多角度對(duì)復(fù)雜地質(zhì)條件下井田各含水層化學(xué)特征與各含水層之間是否存在水力聯(lián)系的相關(guān)問(wèn)題研究較少. 因此,選取補(bǔ)連塔礦區(qū)地下水為研究對(duì)象,通過(guò)繪制主要含水層地下水化學(xué)特征識(shí)別圖版[9],揭示采煤驅(qū)動(dòng)下地下水化學(xué)特征與含水層聯(lián)通關(guān)系,為礦區(qū)快速定位礦井突水來(lái)源提供科學(xué)參考.

    1 材料與方法

    1.1 研究區(qū)概況

    補(bǔ)連塔礦區(qū)位于(109°33¢~110°10¢E,39°01¢~ 39°30¢N)內(nèi)蒙古西部干旱半干旱地區(qū)(圖1),年降雨量396.8mm,年蒸發(fā)量2772mm,面積約為130.9km2,地勢(shì)較為平坦,整體呈東高西低態(tài)勢(shì),黃河一級(jí)支流烏蘭木倫河位于礦區(qū)東側(cè)邊界,受當(dāng)?shù)夭擅菏韪捎绊?烏蘭木倫河及研究區(qū)內(nèi)三條支溝為季節(jié)性流水.

    圖1 研究區(qū)概況

    1.2 水文地質(zhì)條件

    補(bǔ)連塔礦區(qū)地下水類型包括松散巖類孔隙水、松散巖類裂隙水與基巖類裂隙水,含水層順序由上至下分別為第四系含水層、直羅組含水層與延安組含水層.

    松散巖類孔隙水的含水層為第四系含水層,厚度較薄,厚度約在10~30m,主要賦存于烏蘭木倫河兩側(cè)及其一級(jí)支流、河漫灘、一級(jí)階地及局部沙丘區(qū)域,降水入滲是孔隙水的主要補(bǔ)給來(lái)源,此外還接受灌溉回流與山區(qū)基巖側(cè)向補(bǔ)給.松散巖類裂隙水的含水層為直羅組含水層與下覆延安組含水層,直羅組含水層由西向東逐漸增厚,厚度約為10~120m之間,含水層在研究區(qū)西北側(cè)存在出露,出露處直接受風(fēng)化作用侵蝕;延安組為主要煤系地層,受采動(dòng)作用強(qiáng)烈.

    主要隔水層為兩層,一:是介于第四系含水層與直羅組含水層的砂質(zhì)泥巖隔水層;二:是介于直羅組含水層與延安組含水層的砂質(zhì)泥巖隔水層.[10]

    圖2 研究區(qū)水文地質(zhì)示意

    1.3 樣品處理與分析

    根據(jù)研究區(qū)氣候與水文地質(zhì)特征,在2021年8月進(jìn)行了水樣采集,采樣過(guò)程嚴(yán)格遵守《地下水環(huán)境監(jiān)測(cè)技術(shù)規(guī)范》[11](HJ/T164-2004),以烏蘭木倫河及其上游呼和烏素支流、下游活雞兔溝支流為界限,進(jìn)行網(wǎng)格狀采樣,并適當(dāng)向研究區(qū)外延伸,按照水樣采集深度分為4類:地表水采樣點(diǎn)6個(gè)、第四系含水層采樣點(diǎn)11個(gè)、直羅組含水層采樣點(diǎn)6個(gè)、延安組含水層采樣點(diǎn)6個(gè)(表1).

    水樣裝于500mL水樣瓶?jī)?nèi),每個(gè)點(diǎn)采集3瓶,其中一瓶加入硝酸保存劑測(cè)試重金屬離子,并使用Parafilm封口膜密封保存.后將水樣送往內(nèi)蒙古自治區(qū)生態(tài)環(huán)境科學(xué)研究院進(jìn)行測(cè)試,共測(cè)試Na+、K+、Ca2+、Mg2+、F-、Cl-、SO42-、NO3-、CO32-、HCO3-、TDS、pH值、D、18O、3H共15項(xiàng)指標(biāo).

    表1 各層含水層采樣點(diǎn)編號(hào)

    1.4 研究方法與數(shù)據(jù)處理

    對(duì)水樣測(cè)試數(shù)據(jù)進(jìn)行數(shù)理統(tǒng)計(jì)分析[12-13];繪制不同含水層paper三線圖來(lái)判斷地下水化學(xué)主導(dǎo)類型;通過(guò)Schoeller圖[14-15]判斷各含水層離子置換與微量元素賦存情況;借助Durov圖[16-17]對(duì)地下水整體進(jìn)行水化學(xué)分類;根據(jù)以上分析結(jié)果,建立該研究區(qū)常規(guī)水化學(xué)組分識(shí)別圖版.

    2 結(jié)果與分析

    2.1 常規(guī)離子分析

    2.1.1 地表水與第四系含水層對(duì)比 由圖3可看出:地表水中Na+、K+、SO42-、Cl-富集,其中SO42-、Na+離子占比達(dá)70%以上;第四系含水層中Na+、K+有所下降,Ca2+、Mg2+離子濃度開始上升至22%左右,略少于Na+、K+的26%,依舊為Na+、SO42-為主導(dǎo),研究區(qū)的地下水主要為大氣降水補(bǔ)給,降水入滲后,逐漸與周圍巖石發(fā)生溶濾作用和陽(yáng)離子交換吸附作用[18],導(dǎo)致Na+、K+離子濃度降低而Ca2+、Mg2+離子濃度升高.

    地表水采樣點(diǎn)地表水03-地表水01的化學(xué)類型由SO42-×Cl-—Ca2+×Mg2+型轉(zhuǎn)為SO42-×Cl--Mg2+型,鄰近采樣路徑第四系01-第四系02下,第四系含水層采樣點(diǎn)的化學(xué)類型由SO42-×Cl--Ca2+型轉(zhuǎn)為SO42-×Cl--Mg2+,地表水與第四系含水層水樣空間變化趨勢(shì)一致,說(shuō)明地表水與第四系含水層之間存在水力聯(lián)系.

    2.1.2 第四系含水層與直羅組含水層對(duì)比 由圖4可知:直羅組含水層較第四系含水層的采樣點(diǎn)中的SO42-濃度較高且數(shù)值較為穩(wěn)定,HCO3-、Mg2+離子濃度增高,而Cl-離子濃度降低,且降幅較大,地下水化學(xué)類型為Na+-SO42-主導(dǎo)型,兩個(gè)含水層的水化學(xué)特征具有相同的趨勢(shì).考慮到第四系含水層與直羅組含水層之間的隔水層為粉砂巖與粉砂質(zhì)泥巖,具有較好的隔水性,所以這兩層含水層本應(yīng)屬于兩個(gè)相對(duì)獨(dú)立的含水系統(tǒng),但直羅組地層受風(fēng)化侵蝕作用[19],導(dǎo)致地層厚度不一,東部剝蝕嚴(yán)重,逐漸向西部增厚,最厚處可達(dá)160m左右.所以含水的砂巖巖層厚度變化較大,在東部邊緣存在尖滅現(xiàn)象,且上部裂隙較為發(fā)育,導(dǎo)致第四系含水層與直羅組含水層存在水力聯(lián)系.

    圖3 地表水與第四系含水層piper三線圖

    圖4 第四系含水層與直羅組含水層piper三線圖

    2.1.3 直羅組含水層與延安組含水層對(duì)比 延安組為研究區(qū)主要煤系地層,由圖5可知:延安組含水層部分水樣落點(diǎn)與直羅組含水層水樣落點(diǎn)出現(xiàn)重合,而另一部分落點(diǎn)與直羅組含水層水樣落點(diǎn)相差較大,在數(shù)值上整體具有較高SO42-與HCO3-濃度,Cl-離子濃度大幅下降,地下水化學(xué)類型為HCO3--Ca2+主導(dǎo)型.piper三線圖中延安組2號(hào)、3號(hào)水樣依舊與直羅組含水層落點(diǎn)具有較強(qiáng)的關(guān)聯(lián)性,但延安組1, 4, 5, 6號(hào)水樣落點(diǎn)獨(dú)立于直羅組含水層水樣落點(diǎn),推測(cè)出現(xiàn)這種情況的原因?yàn)?采煤驅(qū)動(dòng)作用下煤層上覆巖層裂隙發(fā)育,直羅組含水層中的地下水滲入延安組含水層[25-28],導(dǎo)致延安組含水層的地下水化學(xué)組分變得更為復(fù)雜,位于工作面上方的延安組水樣點(diǎn)2、3受擾動(dòng)影響較大;而位于工作面下方的地下水受采煤帶來(lái)的擾動(dòng)作用較小[20-22],所以延安組水化學(xué)特征與直羅組含水層的差異性更為顯著,直羅組含水層與延安組含水層是否存在水力聯(lián)系,需通過(guò)氫氧同位素進(jìn)一步判斷.

    圖5 直羅組含水層與延安組含含水層piper三線圖

    2.2 氘氧同位素分析

    以標(biāo)準(zhǔn)海洋水的千分差(SMOW)表示所采水樣的氘氧值:

    式中::代表樣品;:標(biāo)準(zhǔn)樣品;:同位素比率.

    以鄂爾多斯市的氘氧同位素?cái)?shù)據(jù)建立當(dāng)?shù)氐慕邓€方程(LMWL)為[23-24]:

    dD=6.5d18O-4.2(2)

    以烏蘭木倫河及其季節(jié)性支流所采水樣的氘氧值計(jì)算當(dāng)?shù)睾恿髡舭l(fā)線方程(EL)和第四系含水層蒸發(fā)線方程分別為:

    dD=3.85d18O-32.21(3)

    dD=4.76d18O-23.41(4)

    通過(guò)克雷格溫度效應(yīng)公式進(jìn)行誤差驗(yàn)證:

    d18O=0.695t(℃)-13.6(‰)(5)

    dD=5.6t(℃)-100(‰)(6)

    式中:為當(dāng)?shù)仄骄鶜鉁?補(bǔ)連塔平均氣溫6.2℃,將帶入(5)、(6)中,得到18O=-9.291、D=65.28,將18O=-9.21帶入(3)式,得D=67.42,與通過(guò)(6)式得到的D=65.28誤差僅為2.14‰,說(shuō)明(2)式可以代表研究區(qū)降水線.

    圖6 樣品dD、d18O關(guān)系

    由圖6可以看出,在研究區(qū)內(nèi),滲入水與地下水充分混合后,有28例樣品值仍位于當(dāng)?shù)亟邓€右下方,D、18O沿當(dāng)?shù)亟邓€呈線性分布,說(shuō)明研究區(qū)地下水的主要補(bǔ)給來(lái)源為大氣降水入滲補(bǔ)給,第四系含水層部分水樣出現(xiàn)18O漂移現(xiàn)象[31],可能受到蒸發(fā)作用影響.

    地表水水樣與第四系含水層水樣的D、18O分布范圍相互重疊,反映了地表水入滲補(bǔ)給過(guò)程;第四系含水層水樣的D則普遍貧于地表水,但豐于直羅組含水層與延安組含水層水樣,這也符合第四系含水層為較年輕地層的實(shí)際情況[25];第四系含水層水樣與直羅組含水層水樣的落點(diǎn)相互重疊,說(shuō)明兩者間具有一定的水力聯(lián)系,部分水樣的落點(diǎn)接近重合,說(shuō)明兩者間可能存在垂向補(bǔ)給關(guān)系,進(jìn)一步證明了第四系含水層地下水與直羅組含水層地下水存在水力聯(lián)系;直羅組含水層與部分延安組含水層落點(diǎn)相接近,說(shuō)明直羅組含水層與延安組含水層之間存在一定的水力聯(lián)系,但另一部分水樣點(diǎn)D、18O出現(xiàn)貧化,說(shuō)明這部分采樣點(diǎn)可能位于延安組含水層較深位置或采樣點(diǎn)位于冒落、導(dǎo)水裂隙帶影響區(qū)域邊緣或影響區(qū)域外,受采動(dòng)作用影響較小,含水層貯水條件更為封閉,導(dǎo)致這部分水樣與上覆含水層的聯(lián)系較小.

    2.3 3H元素分析

    氚元素的更新周期要慢于氘元素[26],借此可以測(cè)定氚元素來(lái)判斷水體的形成年代,第四系含水層水樣(埋深小于30m)3H值介于12.4~22.5TU,均值為17.1TU,說(shuō)明第四系含水層所含地下水均來(lái)自現(xiàn)代大氣降水;直羅組含水層水樣(埋深位于20~ 120m)3H介于6.3~22.4TU,均值為14.1TU,在數(shù)值上略小于第四系含水層,說(shuō)明直羅組含水層所含地下水為次現(xiàn)代-近代混合補(bǔ)給水;延安組含水層水樣3H介于4.3~18.6TU,均值為11.4TU,整體小數(shù)值上于直羅組含水層,說(shuō)明延安組含水層地下水的貯藏條件相對(duì)封閉,但與直羅組含水層水樣在數(shù)值上仍然存在交叉,判斷延安組含水層與直羅組含水層存在水力聯(lián)系,所含地下水為次現(xiàn)代-近代混合地下水.由圖7可以看出:數(shù)值偏小的水樣采樣點(diǎn)應(yīng)該位于受采煤驅(qū)動(dòng)擾動(dòng)較小的區(qū)域,處于相對(duì)封閉的狀態(tài),與上覆含水層的聯(lián)系較小,所含地下水主要為次現(xiàn)代補(bǔ)給水.

    圖7 3H含量

    綜上所述,研究區(qū)地下水的主要補(bǔ)給來(lái)源為大氣降水,在烏蘭木倫河沿岸地區(qū),由于煤層埋深較淺,采空區(qū)冒落、裂隙帶發(fā)育高度甚至可以直達(dá)地表,在長(zhǎng)期采動(dòng)作用下,煤系地層上覆含水層原有裂隙進(jìn)一步發(fā)育至貫穿第四系含水層底板,使第四系含水層、直羅組含水層與延安組含水層之間產(chǎn)生水力聯(lián)系.

    2.4 微量離子分析

    2.4.1 NO3-的變化趨勢(shì)分析 選取NO3-和F-離子作為參考元素進(jìn)行分析,由圖8可知,在地表采樣點(diǎn)所采水樣的NO3-離子濃度較為穩(wěn)定,入滲補(bǔ)給進(jìn)入第四系含水層后濃度逐漸升高,直羅組、延安組含水層中NO3-濃度整體上低于第四系含水層,推測(cè)第四系含水層NO3-濃度增高的原因?yàn)?研究區(qū)地表存在大面積塌陷回填后形成松樹種植區(qū),所需施氮肥含量遠(yuǎn)高于周邊稀疏分布的耕地,所以研究區(qū)第四系含水層硝酸的濃度主要受林地影響[27-29].結(jié)合研究區(qū)高程分布情況來(lái)看,第四系07采樣點(diǎn)海拔較高,不利于NO3-在地下水中富集,因此該點(diǎn)的濃度低于第四系含水層其他采樣點(diǎn).

    2.4.2 F-的變化趨勢(shì)分析 如圖8所示,F-離子濃度在地表水05點(diǎn)數(shù)值較大,地表水05點(diǎn)所在位置西北600m處為補(bǔ)連塔煤礦的露天礦堆渣場(chǎng),導(dǎo)致F-離子濃度在地表水05點(diǎn)數(shù)值較高,受煤矸石堆放場(chǎng)影響導(dǎo)致F-含量在地表水中數(shù)值整體略高于第四系含水層水樣.

    圖8 微量元素濃度折線

    直羅組地下水在進(jìn)入延安組含水層后,F-離子濃度明顯上升,可能有3個(gè)原因[30]:(1)延安組含水層為補(bǔ)連塔礦區(qū)的主要采煤層,巖層內(nèi)富含有機(jī)礦物,與地下水發(fā)生溶濾作用后,導(dǎo)致F-濃度升高[35]; (2)少量F-會(huì)與Ca2+反應(yīng)生成極難溶于水的CaF2后達(dá)到水巖平衡(圖9),所以F_易隨水體遷移;(3)采動(dòng)產(chǎn)生降落漏斗導(dǎo)致水力梯度增大,直接影響地下水徑流,并形成局部地下水水流系統(tǒng),上覆含水層地下水通過(guò)導(dǎo)水裂隙溶解了更多的F-,并最終匯入到延安組含水層中,同時(shí)延安為主要煤系地層,煤層開采使得地下水與圍巖充分接觸,導(dǎo)致延安組F-離子濃度增高[31].

    由圖9可見(jiàn):在水體呈弱堿性的情況下,F-、Cl-濃度偏高,水體TDS整體偏高.一般情況下,高TDS、低Ca2+環(huán)境為高F-的生成條件,對(duì)于補(bǔ)連塔礦區(qū)延安組含水層,高TDS的同時(shí),Ca2+濃度并不低,但F-濃度依舊較高,因此F-存在多個(gè)來(lái)源,延安組含水層與上覆含水層存在較密切的水力聯(lián)系.

    圖9 各化學(xué)組分Scholler圖

    長(zhǎng)時(shí)間的采動(dòng)作用會(huì)導(dǎo)致:(1)上覆巖層的裂隙發(fā)育,裂隙連通性增強(qiáng);(2)局部第四系含水層與直羅組含水層會(huì)通過(guò)裂隙向下層延安組含水層進(jìn)行補(bǔ)給,形成新的匯流中心.綜上所述,第四系含水層地下水與直羅組相近的元素變化趨勢(shì)(圖9)說(shuō)明:第四系含水層與直羅組含水層存在密切的水力聯(lián)系.延安組含水層在變化趨勢(shì)上與直羅組含水層整體相近,且出現(xiàn)F-離子的逐層累計(jì)現(xiàn)象(圖8),但部分點(diǎn)在常規(guī)元素上與直羅組含水層相差較大說(shuō)明:直羅組含水層與延安組含水層存在有限的水力聯(lián)系.

    3 水化學(xué)特征識(shí)別圖版

    3.1 識(shí)別流程

    研究區(qū)各含水層相互存在水力聯(lián)系,導(dǎo)致不同含水層水樣整體在數(shù)值上出現(xiàn)一定程度共同性.但受不同地質(zhì)環(huán)境、貯藏條件與煤礦開采條件的影響,各含水層的地下水化學(xué)特征同時(shí)存在明顯的差異性,差異主要反映在HCO3-、Mg2+、Cl-、SO42-、Ca2+、Na+、K+濃度的變化上.通過(guò)各層地下水的差異性將識(shí)別圖版劃分為4個(gè)區(qū)域,分別對(duì)應(yīng)地表水、第四系、直羅組、延安組含水層,由此建立礦區(qū)突水水源識(shí)別圖版.在礦區(qū)發(fā)生突水事故時(shí),通過(guò)分析突水水樣的化學(xué)組分特征值以及在識(shí)別圖版上的落點(diǎn)分區(qū),來(lái)確定突水來(lái)源,同時(shí)結(jié)合不同類型的識(shí)別圖版,提高結(jié)果可靠性,以達(dá)到快速確定礦區(qū)突水來(lái)源的目的.

    3.2 常規(guī)元素識(shí)別圖版

    由圖10可知,地表水主要集中在Cl--SO42-- HCO3-三角形的SO42-端,落點(diǎn)集中于①區(qū);第四系含水層主要集中在Cl--SO42--HCO3-三角形的Ca2++ Mg2+端,落點(diǎn)集中于②區(qū);直羅組含水層主要集中在Mg2+-Ca2+-K+三角形的SO42-端,落點(diǎn)集中于③區(qū);延安組含水層主要集中在Na++K+端,落點(diǎn)集中于④區(qū).通過(guò)常規(guī)元素實(shí)現(xiàn)的判斷過(guò)程,主要由Cl--SO42-- HCO3-三角形和Mg2+-Ca2+-K+三角形完成.由圖11可知:在各含水層原有的地下水化學(xué)類型特征下,地表水與第四系含水層落點(diǎn)集中于TDS較低的①、②兩區(qū),直羅組含水層與延安組含水層落點(diǎn)集中于TDS較高的③、④兩區(qū),同時(shí)查看Piper識(shí)別圖版與Durov圖版可更準(zhǔn)確地判斷地下水來(lái)源.

    圖10 Piper三線圖識(shí)別圖版

    圖11 Durov識(shí)別圖版

    3.2 氘氧同位素識(shí)別圖版

    由圖12可知:地表水的D、18O值最大,處于偏重端,落點(diǎn)集中于①區(qū);第四系含水層D、18O值豐于地表水且落點(diǎn)較分散,多集中于②區(qū);直羅組含水層D、18O值豐于第四系含水層,落點(diǎn)多集中于③區(qū);延安組含水層D、18O同位素貧化,明顯處于偏輕端,落點(diǎn)多集中于④區(qū),該圖版可清晰快速地辨別水源是否來(lái)自延安組含水層.

    圖12 dD、d18O同位素識(shí)別圖版

    4 采煤驅(qū)動(dòng)下地下水循環(huán)模式構(gòu)建

    通過(guò)常規(guī)元素與氫氧同位素判斷含水層之間是否存在水力聯(lián)系,通過(guò)微量元素判斷含水層之間水力聯(lián)系的強(qiáng)弱,構(gòu)建了采動(dòng)作用下地表水、地下水循環(huán)模式,分別為孔隙地下水流系統(tǒng)(第四系含水層)、裂隙地下水流系統(tǒng)(直羅組含水層與延安組含水層)與局部地下水流系統(tǒng)三個(gè)層級(jí).采煤驅(qū)動(dòng)下形成的局部地下水流系統(tǒng)接受大氣降水、部分地表水、孔隙水與裂隙水的集中補(bǔ)給,改變了地表水、孔隙水、裂隙水原有的循環(huán)形式,并在局部形成新的匯流中心(如圖13所示).

    圖13 補(bǔ)連塔地下水循環(huán)概念模型

    5 結(jié)論

    5.1 地表水化學(xué)類型以Mg2+-SO42-×Cl-為主導(dǎo)型,第四系含水層水化學(xué)類型同樣為Mg2+-SO42-×Cl-主導(dǎo)型,二者化學(xué)類型轉(zhuǎn)化在空間上具有一致性;直羅組含水層水化學(xué)類型為Na+-SO42-主導(dǎo)型,與第四系含水層相比,HC03-、Mg2+濃度上升,Cl-離子濃度下降,且減幅較大;延安組含水層水化學(xué)類型為HCO3--Ca2+主導(dǎo)型,與直羅組含水層相比,SO42-與HCO3-濃度較高,Cl-離子濃度大幅下降.受采動(dòng)作用影響,3個(gè)含水層中F-呈逐層遞增的趨勢(shì),延安組含水層F-濃度偏高;受植被施肥影響,第四系含水層中NO3-濃度偏較高,研究區(qū)植被為第四系含水層NO3-含量大幅升高的主要原因.

    5.2 地表水與第四系含水層存在相互轉(zhuǎn)換的關(guān)系,存在密切水力聯(lián)系;第四系含水層與直羅組含水層存在密切水力聯(lián)系;直羅組含水層與延安組含水層存在有限的水力聯(lián)系.

    5.3 通過(guò)常規(guī)離子、微量元素、氘氧同位素與氚同位素等多種分析手段對(duì)復(fù)雜地質(zhì)條件下井田不同層位水體進(jìn)行識(shí)別與分析,總結(jié)了不同含水層地下水化學(xué)類型與水力聯(lián)系,并建立Piper識(shí)別圖版、Durov圖版、氘氧同位素圖版,通過(guò)聯(lián)合比對(duì)可快速開展礦區(qū)突水溯源,并借此提出針對(duì)性的補(bǔ)救措施.

    [1] 楊永剛,孟志龍,秦作棟,等.采礦對(duì)山西水資源破壞過(guò)程的同位素示蹤研究[J]. 中國(guó)環(huán)境科學(xué), 2013,33(8):1447-1453.

    Yang Y G, Meng Z L, Qin Z D, et al. Isotope tracer study on the destruction processs of water resources by miningi n Shanxi Province [J]. China Environmental Science, 2013,33(8):1447-1453.

    [2] 陳 雪,許丹丹,錢雅慧,等.淮北劉橋二礦煤矸石多環(huán)芳烴污染特征及毒性評(píng)價(jià)[J]. 中國(guó)環(huán)境科學(xué), 2022,42(2):753-760.

    Chen X, Xu D D, Qian Y H, et al. Pollution characteristics and toxicity evaluation of coal gangue polycyclic aromatic hydro carbonsin Liuqiao No.2 Coal Mine, Huaibei [J]. China Environmental Science, 2022,42(2):753-760.

    [3] 錢鳴高,許家林.煤炭開采與巖層運(yùn)動(dòng) [J]. 煤炭學(xué)報(bào), 2019,44(4): 973-984.

    Qian M G, XuJ L. Behaviorsofstratamovementincoalmining [J]. Journalof China Coal Society, 2019,44(4):973-984.

    [4] 冀瑞君,彭蘇萍,范立民,等.神府礦區(qū)采煤對(duì)地下水循環(huán)的影響——以窟野河中下游流域?yàn)槔?[J]. 煤炭學(xué)報(bào), 2015,40(4):938-943.

    Ji R J, Peng S P, Fan L M, et al. Effectofcoalexploita- tiononground watercirculationinthe Shenfuminearea: Anexam-plefrommiddle and lowerreachesofthe Kuye Riverbasin [J]. Journalof China Coal Society, 2015,40(4):938-943.

    [5] 范立民,向茂西,彭 捷,等.毛烏素沙漠與黃土高原接壤區(qū)泉的演化分析 [J]. 煤炭學(xué)報(bào), 2018,43(1):207-218.

    Fan L M, Xiang M X, Peng J, et al. Evolutionanaly sisonspring sincontiguo usareaof Maowusu Desertand Loess Plateau [J]. Journalof China Coal Society, 2018,43(1):207-218.

    [6] 梁蓉蓉.孔隙含水層底板參數(shù)變化條件下煤礦開采對(duì)松散含水層影響規(guī)律的數(shù)值模擬研究 [D]. 太原:太原理工大學(xué).

    Liang R R. Numerical simulation studyontheinfluence of coal miningonlooseaquiferundertheconditionofchangeofporeaquiferfloorparameters [D]. Taiyuan: Taiyuanuniversityoftechnology.

    [7] 侯澤明,黃 磊,韓 萱,等.采煤驅(qū)動(dòng)下神東礦區(qū)地下水化學(xué)特征及控制因素[J]. 中國(guó)環(huán)境科學(xué), 2022,42(5):2250-2259.

    Hou Z M, Huang L, Han X, et al. Hydrochemical character istics and controlling factors of groundwater driven by coal mining in Shendong mining area [J]. 2022,42(5):2250-2259.

    [8] 顧大釗,張建民,王振榮,等.神東礦區(qū)地下水變化觀測(cè)與分析研究 [J]. 煤田地質(zhì)與勘探, 2013,41(4):35-39.

    Gu D Z, Zhang J M. Wang Z R, et al. Observation and analysis of groundwater change in Shendong mining area [J]. Coal Geology & Exploration, 2013,41(4):35-39.

    [9] 喬元棟,孟召平,張 村,等.復(fù)雜構(gòu)造井田含水層特征及其水力聯(lián)系辨識(shí) [J/OL]. 煤炭學(xué)報(bào): 1-11[2021-12-30].http://kns.cnki.net/ kcms/detail/11.2190.td.20211105.1522.005.html.

    Qiao Y D, Meng Z P, Zhang C, et. Characteristics of mine aqui fer with complex structure and identification of its hydraulic connection [J/OL]. Journal of China Coal Society:1-11[2021-12-30]. http://kns. cnki.net/kcms/detail/11.2190.td.20211105.1522.005.html.

    [10] 劉利軍.補(bǔ)連塔礦區(qū)水文地質(zhì)報(bào)告[R]. 內(nèi)蒙古自治區(qū):神東煤炭集團(tuán)地測(cè)公司, 2007..

    Liu L J. Hydrogeological Report of Bulianta Mining Area [R]. Inner Mongolia: Shendong Coal Group Geological Survey Company, 2007.

    [11] HJ/T164-2004 地下水環(huán)境監(jiān)測(cè)技術(shù)規(guī)范 [S].

    HJ/T164-2004 Groundwaterenvironmentalmonitoringtechnicalspeci fication [S].

    [12] 房麗晶,高瑞忠,賈德彬,等.草原流域地下水化學(xué)時(shí)空特征及環(huán)境驅(qū)動(dòng)因素——以內(nèi)蒙古巴拉格爾河流域?yàn)槔?[J]. 中國(guó)環(huán)境科學(xué), 2021,41(5):2161-2169.

    Fang L J, Gao R Z, Jia D B, et al. Spatial- temporal characteristics of groundwater quality and its environment driving faxtors of Steppe Basin- taken Balager river basin of Inner Mongolia for instance [J]. China Environmental Science, 2021,41(5):2161-2169.

    [13] 王 建,韓海東,許君利,等.塔里木河流域出山徑流水化學(xué)特征研究 [J]. 中國(guó)環(huán)境科學(xué), 2021,41(4):1576-1587.

    Wang J, Han H D, Xu J L, et. Hydro chemical characteristics of the mountain runoff in Tarim River Basin, China [J]. China Environmental Science, 2021,41(4):15,76-1587.

    [14] 劉 鑫,向 偉,司炳成.汾河流域淺層地下水水化學(xué)和氫氧穩(wěn)定同位素特征及其指示意義 [J]. 環(huán)境科學(xué), 2021,42(4):1739-1749.

    Liu X, Xiang W, Si B C, Hydrochemical and isotopic characteristics in the shallow groundwater of the Fenhe River Basin and indicative significance [J]. Environmental Science, 2021,42(4):1739-1749.

    [15] 仁增拉姆,羅 珍,陳虎林,等.西藏年楚河流域水化學(xué)特征分析 [J]. 地球與環(huán)境, 2021,49(4):358-366.

    Renzeng Lamu, Luo Z, Chen H L, et al. Hydrochemical characteristics of Nyangchu River Basinin Tibet [J]. Earthand Environment, 2021, 49(4):358-366.

    [16] 王文祥,李文鵬,蔡月梅,等.黑河流域中游盆地水文地球化學(xué)演化規(guī)律研究 [J]. 地學(xué)前緣, 2021,28(4):184-193.

    Wang W X, Li W P, Cai Y M, et al. Thehydrogeochemical evolution ofgroundwaterinthemiddlereachesofthe Heihe River Basin [J]. Earth Science Frontiers, 2021,28(4):184-193.

    [17] 雷 米,周金龍,張 杰,等.新疆博爾塔拉河流域平原區(qū)地表水與地下水水化學(xué)特征及轉(zhuǎn)化關(guān)系[J]. 環(huán)境科學(xué), 2022,43(4):1873-1884

    Lei M, ZHou J L, Zhang J, et al. Hydrochemical characteristics and transformation relationship of surface water and groundwater in the plain area of Bortala River Basin, Xinjiang [J]. Environmental Science, 2022,43(4):1873-1884.

    [18] 房滿義,李雪妍,張 根,等.大柳塔煤礦地下水庫(kù)水巖作用機(jī)理分析 [J/OL]. 煤炭科學(xué)技術(shù), 1-8[2021-12-31].http://kns.cnki.net/kcms/ detail/11.2402.TD.20200502.0813.002.html.

    Fang M Y, Li X Y, Zhang G, et al. Discussion on water-rock interaction mechanism in underground reservoir of Daliuta coal mine [J/OL]. Coal Scienceand Technology, 1-8[2021-12-31]. http://kns.cnki.net/kcms/ detail/11.2402.TD.20200502.0813.002.html.

    [19] 李 偉..神東礦區(qū)2012年水文地質(zhì)報(bào)告[R]. 內(nèi)蒙古自治區(qū):神東煤炭集團(tuán)地測(cè)公司, 2012.

    Li W. Shendong Mining Area hydrogeological report [R]. Inner Mongolia: Shendong Coal Group Geological Survey Company, 2012.

    [20] 楊雅琪.辛置礦井主要充水含水層地下水化學(xué)特征及成因機(jī)制研究 [D]. 江蘇:中國(guó)礦業(yè)大學(xué), 2020.

    Yang Y Q. Research on groundwater chemical characteristics and gene sis mechanism of main water-filled aquifers in Xinzhi Coal Mine [D]. China Universityof Mining & Technology.

    [21] 樊 燕.煤礦開采對(duì)上覆含水層影響的數(shù)值模擬研究 [D]. 太原:太原理工大學(xué), 2011.

    Fan Y. Numerical simulation study on the influence of coal mining on overlying aquifer [D]. Taiyuan: Taiyuan Universityof Technology.

    [22] 王甜甜,張 雁,趙 偉,等.伊敏礦區(qū)地下水水化學(xué)特征及其形成作用分析 [J]. 環(huán)境化學(xué), 2021,40(5):1480-1489.

    Wang T T, Zhang Y, Zhao W, et al. Hydrogeochemical characteristics and formation process of groundwater in Yimin mining area [J].. Environmental Chemistry, 2021,40(5):1480-1489.

    [23] 賈艷琨,王經(jīng)蘭,王東升.環(huán)境同位素在水文地質(zhì)和環(huán)境地質(zhì)研究中的應(yīng)用 [J]. 地球?qū)W報(bào), 2005,26(Z1):307-308.

    Jia Y K, Wang J L, Wang D S. The application of environmental isotopesto hydrogeology and environmental geolo-gy [J]. Acta Geoscientica Sinica, 2005,26(Z1): 307-308.

    [24] Ayenew T, Kebede S, Alemyahu T. Environmental isotopes and hydrochemical study applied to surface water and groundwater in teraction in the Awash Riverbasin [J]. Hydrological Processes, 2010,22(10): 1548-1563.

    [25] 徐秋娥,劉澄靜,角媛梅,等.穩(wěn)定氫氧同位素示蹤水汽來(lái)源對(duì)哈尼梯田降水補(bǔ)給的影響 [J]. 生態(tài)學(xué)報(bào), 2020,40(5):1709-1717.

    Xu Q E, Liu C J, Jiao Y M, et al. Impacts of stable isotopic composition and moisture sources of precipitation on precipitation recharge of Hani Rice Terraces during the dry season [J]. Acta Ecologica Sinica, 2020,40(5):1709-1717.

    [26] 郭洋楠,楊俊哲,張 政,等.神東礦區(qū)礦井水的氫氧同位素特征及高氟礦井水形成的水-巖作用機(jī)制 [J/OL]. 煤炭學(xué)報(bào):1-15[2021-10- 09].https://doi.org/10.13225/j.cnki.jccs.2021.0388.

    Guo Y N, Yang J Z, Zhang Z, et al. Hydrogenand oxygenlsotope characteristicsof mine water in the Shendong Mine area and water- rock reactions mechanism of the formation of high-fluoride mine water [J]. Journal of China Coal Society, 1-15[2021-10-09].

    [27] 付昌昌,李向全,馬劍飛,等.窟野河流域中游煤礦區(qū)地下水質(zhì)量及補(bǔ)給來(lái)源研究 [J]. 水文, 2018,38(6):42-47.

    Fu C C, Li X Q, Ma J F, et al. Study on groundwater quality and recharge sources in middler eaches of Kuye River Basin [J]. hydrology, 2018,38(6):42-47.

    [28] 梁向陽(yáng),方 剛,黃 浩.榆神礦區(qū)曹家灘井田水文地球化學(xué)特征研究 [J]. 干旱區(qū)資源與環(huán)境, 2020,34(5):102-108.

    Liang X Y, Fang G, Huang H. Hydrogeochemical characteristics of Caojiatan mine field in Yushen Mining Area [J]. Journalof Arid Land Resourcesand Environment, 2020,34(5):102-108.

    [29] 黃平華,祝金峰,鄧 勇,等.地下水中氚同位素分布模型及其應(yīng)用 [J]. 煤炭學(xué)報(bào), 2013,38(S2):448-452.

    Huang P H, Zhu J F, Deng Y, et al. Distribution model and application of groundwatert ritium isotope [J]. Journal of the China Coal Society, 2013,38(S2):448-452.

    [30] 張 鑫,張 妍,畢直磊,等.中國(guó)地表水硝酸鹽分布及其來(lái)源分析 [J]. 環(huán)境科學(xué), 2020,41(4):1594-1606.

    Zhang X, Zhang Y, Bi Z L, et al. Distribution and source anal ysis of nitratein surfacewaters of China [J]. Environmental Earth Sciences, 2020,41(4):1594-1606.

    [31] 蘇 賀,康衛(wèi)東,楊永康.基于水化學(xué)和穩(wěn)定同位素的黃土區(qū)地下水硝酸鹽來(lái)源示蹤 [J]. 太原理工大學(xué)學(xué)報(bào), 2021,52(5):775-788.

    Su H, Kang W D, Yang Y K. Source tracing of ground water nitrate in loess area based on hydrochemistry and stable isotope [J]. Journal of Taiyuan University of Technology, 2021,52(5):775-788.

    [32] 王 磊.地下水中硝酸鹽氮污染源解析 [D]. 北京:中國(guó)地質(zhì)大學(xué)(北京), 2016.

    Wang L. Analysis of nitrate nitrogen pollution sources in groundwater [D]. Beijing: China University of Geosciences (Beijing).

    [33] 郝春明,張 偉,何瑞敏,等.神東礦區(qū)高氟礦井水分布特征及形成機(jī)制 [J]. 煤炭學(xué)報(bào), 2021,46(6):1966-1977.

    Hao C M, Zhang W, He R M, et al. Formation mechanisms for elevated fluoride in the minewater in Shendong coal-mining district. [J] Journalof China Coal Society, 2021,46(6):1966-1977.

    [34] 唐春雷,鄭秀清,梁永平.龍子祠泉域巖溶地下水水化學(xué)特征及成因 [J]. 環(huán)境科學(xué), 2020,41(5):2087-2095.

    Tang C L, Zheng X Q, Liang Y P. Hydrochemical characteristics and formation causes of ground karstwater systems in the Longzici Spring Catchment [J]. Environmental Science, 2020,41(5):2087-2095.

    Identification of chemical characteristics and hydraulic connection of each aquifer in complex mine field driven by coal mining.

    HUANG Lei1*, HOU Ze-ming1, HAN Xuan1, XU Lei1, ZHANG Sheng-wei1,2, LI Gang-zhu1, LIU Zhi-qiang1

    (1. Inner Mongolia Water Resource Protection and Utilization Key Laboratory, Water Conservancy and Civil Engineering College, Inner Mongolia Agricultural University, Hohhot 010018, China;2.Inner Mongolia Autonomous Region Key Laboratory of Big Data Research and Application of Agriculture and Animal Husbandry, Hohhot 010018, China)., 2022,42(6):2697~2706

    By analyzing the numerical characteristics of conventional elements, trace elements, deuterium, oxygenand tritium isotope of different aquifers, the groundwater chemical characteristics and hydraulic connection of each aquifer were judged. The three identification charts (piper and durov identification charts of conventional hydrochemical components and deuterium tritium isotope identification charts) were established to realize the rapid identification of aquifer water source. The results showed that the main recharge source of groundwater in the study area was atmospheric precipitation. There was a significant indigenous connection between Quaternary aquifer and Zhiluo Formation aquifer, limited indigenous connection between the Zhiluo Formation aquifer and the Yan'an Formation aquifer. The three identification charts were established to determine the water inrush source and put forward targeted solutions for water inrush from different aquifers in Bulianta mining area based on the differences in hydrochemical ion characteristics, deuterium-oxygen and tritium isotopes.

    complex mine fields;multi-layer aquifer;hydraulic connection;mine water inrush;hydrochemical characteristics;isotopetracing

    X523

    A

    1000-6923(2022)06-2697-10

    黃 磊(1983-)男,內(nèi)蒙古呼和浩特人,博士,副教授,主要從事水文地質(zhì)與環(huán)境地質(zhì)領(lǐng)域的教學(xué)與科研工作.發(fā)表論文20余篇.

    2022-01-05

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51969023);內(nèi)蒙古自治區(qū)科技計(jì)劃項(xiàng)目(2020GG0076)

    * 責(zé)任作者, 副教授, lei_huang@yeah.net

    猜你喜歡
    直羅第四系圖版
    巨厚充水含水層疏放水層位優(yōu)化研究
    鄂爾多斯盆地北部直羅組含水層研究進(jìn)展與水害防治建議
    第四系膠結(jié)砂層水文地質(zhì)鉆進(jìn)技術(shù)研究初探
    河北平原區(qū)第四系(200m以淺)巖土體綜合熱導(dǎo)率特征綜述
    北疆第四系覆蓋層上明渠結(jié)構(gòu)優(yōu)化研究
    陸相頁(yè)巖氣資源評(píng)價(jià)初探:以延長(zhǎng)直羅——下寺灣區(qū)中生界長(zhǎng)7段為例
    圖版 Ⅰ Plate Ⅰ
    圖版II Plate II
    直羅油田三疊系延長(zhǎng)組長(zhǎng)8段油氣成藏條件分析
    圖版
    国产一区二区激情短视频| 一级a爱片免费观看的视频| 97人妻精品一区二区三区麻豆| 欧美黑人欧美精品刺激| 亚洲精品久久国产高清桃花| 精品久久久久久成人av| 亚洲色图av天堂| 最好的美女福利视频网| or卡值多少钱| 97超级碰碰碰精品色视频在线观看| 精品国内亚洲2022精品成人| 三级国产精品欧美在线观看| 丰满人妻一区二区三区视频av | 人妻丰满熟妇av一区二区三区| 成年版毛片免费区| 国内揄拍国产精品人妻在线| 亚洲七黄色美女视频| 99久久精品国产亚洲精品| 精品久久久久久久久久久久久| 天堂影院成人在线观看| 高清在线国产一区| 亚洲人成网站高清观看| 男女做爰动态图高潮gif福利片| 久久人人精品亚洲av| 国产精品av视频在线免费观看| 国产一区二区在线av高清观看| 欧美性感艳星| 男人的好看免费观看在线视频| 女警被强在线播放| 国产三级在线视频| 欧美一区二区精品小视频在线| 91av网一区二区| 欧美bdsm另类| 男人舔女人下体高潮全视频| avwww免费| 男人舔女人下体高潮全视频| 国产精品一及| 人妻夜夜爽99麻豆av| 成人国产一区最新在线观看| 精品电影一区二区在线| 欧美一级a爱片免费观看看| 午夜免费激情av| 亚洲片人在线观看| 蜜桃久久精品国产亚洲av| 长腿黑丝高跟| 又黄又爽又免费观看的视频| 国产极品精品免费视频能看的| 国产爱豆传媒在线观看| 亚洲美女视频黄频| 国产av不卡久久| 亚洲av一区综合| 久99久视频精品免费| 性欧美人与动物交配| 精品一区二区三区av网在线观看| 一级毛片高清免费大全| 婷婷精品国产亚洲av| 精品国产亚洲在线| 国产色爽女视频免费观看| 男人舔女人下体高潮全视频| 在线国产一区二区在线| 高清在线国产一区| 人妻久久中文字幕网| 非洲黑人性xxxx精品又粗又长| 亚洲精品久久国产高清桃花| 天天躁日日操中文字幕| 色播亚洲综合网| or卡值多少钱| 国产一级毛片七仙女欲春2| 偷拍熟女少妇极品色| 久久伊人香网站| 免费高清视频大片| 欧美绝顶高潮抽搐喷水| 综合色av麻豆| 俄罗斯特黄特色一大片| 搡老妇女老女人老熟妇| 女人十人毛片免费观看3o分钟| 丝袜美腿在线中文| 欧美xxxx黑人xx丫x性爽| 国产精品久久久久久人妻精品电影| 少妇的逼好多水| 久久久久久国产a免费观看| 国产欧美日韩一区二区精品| 亚洲av中文字字幕乱码综合| 一区二区三区高清视频在线| 九九热线精品视视频播放| 中文字幕精品亚洲无线码一区| 看片在线看免费视频| 又粗又爽又猛毛片免费看| 国产黄a三级三级三级人| 精品熟女少妇八av免费久了| 老司机深夜福利视频在线观看| 午夜福利18| 欧美一区二区国产精品久久精品| xxxwww97欧美| 国产乱人伦免费视频| 国产精品久久久久久精品电影| 久久精品91无色码中文字幕| 观看免费一级毛片| 大型黄色视频在线免费观看| 制服丝袜大香蕉在线| 欧美3d第一页| 免费av不卡在线播放| 久久精品国产清高在天天线| 国产精品嫩草影院av在线观看 | 国产精品 国内视频| 欧美国产日韩亚洲一区| av天堂中文字幕网| 国产成年人精品一区二区| 国产欧美日韩一区二区精品| 欧美绝顶高潮抽搐喷水| 宅男免费午夜| 亚洲欧美日韩高清专用| 亚洲国产日韩欧美精品在线观看 | 亚洲国产中文字幕在线视频| 在线天堂最新版资源| 五月玫瑰六月丁香| 夜夜夜夜夜久久久久| 无限看片的www在线观看| 国产不卡一卡二| 精品久久久久久久末码| 国产成人a区在线观看| 成人av一区二区三区在线看| 九九热线精品视视频播放| e午夜精品久久久久久久| 特级一级黄色大片| 亚洲国产精品sss在线观看| 两个人视频免费观看高清| 亚洲精品乱码久久久v下载方式 | 99久国产av精品| 国产毛片a区久久久久| 一a级毛片在线观看| 18美女黄网站色大片免费观看| 国产色爽女视频免费观看| 每晚都被弄得嗷嗷叫到高潮| 国产精品自产拍在线观看55亚洲| 国产精品 欧美亚洲| 亚洲中文字幕日韩| 香蕉丝袜av| 欧美性猛交╳xxx乱大交人| 美女高潮的动态| 岛国在线观看网站| 国产精品免费一区二区三区在线| 一进一出抽搐动态| 国产在线精品亚洲第一网站| 亚洲人成网站高清观看| 啦啦啦观看免费观看视频高清| 亚洲成人久久爱视频| 欧美成人a在线观看| 久久婷婷人人爽人人干人人爱| 免费无遮挡裸体视频| 法律面前人人平等表现在哪些方面| 久久久久久国产a免费观看| 嫩草影视91久久| 久久精品国产综合久久久| 日日夜夜操网爽| 麻豆一二三区av精品| 又黄又粗又硬又大视频| 成人国产综合亚洲| 特级一级黄色大片| 欧美精品啪啪一区二区三区| 精品人妻偷拍中文字幕| 欧美+亚洲+日韩+国产| 久久中文看片网| 丰满乱子伦码专区| 国产精品嫩草影院av在线观看 | 精品人妻1区二区| 亚洲av美国av| 嫩草影院入口| 少妇的丰满在线观看| 亚洲色图av天堂| 国产精品一区二区三区四区久久| 国产成人欧美在线观看| 午夜福利在线观看免费完整高清在 | 97碰自拍视频| 18禁裸乳无遮挡免费网站照片| 国产精品电影一区二区三区| 久久精品国产亚洲av香蕉五月| 99在线视频只有这里精品首页| 久久久久免费精品人妻一区二区| 国产成年人精品一区二区| 很黄的视频免费| 亚洲av五月六月丁香网| 国产一区二区激情短视频| 99久国产av精品| 日本黄色视频三级网站网址| 乱人视频在线观看| 国产成人啪精品午夜网站| 免费在线观看影片大全网站| 久久久久久久精品吃奶| 国产激情欧美一区二区| 亚洲国产精品999在线| 在线观看午夜福利视频| 国产伦在线观看视频一区| 他把我摸到了高潮在线观看| 中文字幕精品亚洲无线码一区| 激情在线观看视频在线高清| a级毛片a级免费在线| 熟女少妇亚洲综合色aaa.| 免费在线观看影片大全网站| 一夜夜www| 亚洲人成网站高清观看| 男人舔女人下体高潮全视频| 国产精品美女特级片免费视频播放器| 国产乱人伦免费视频| 国产一区二区在线av高清观看| 99久久综合精品五月天人人| 国产乱人视频| 午夜免费激情av| 欧美乱码精品一区二区三区| www.色视频.com| bbb黄色大片| 激情在线观看视频在线高清| 天天添夜夜摸| 欧美午夜高清在线| 亚洲五月婷婷丁香| 99久久九九国产精品国产免费| 国产免费av片在线观看野外av| 69人妻影院| 欧美3d第一页| 国产午夜精品论理片| 国产乱人视频| 偷拍熟女少妇极品色| 99热6这里只有精品| 亚洲天堂国产精品一区在线| 淫妇啪啪啪对白视频| 欧美国产日韩亚洲一区| 一本精品99久久精品77| 99国产综合亚洲精品| 色噜噜av男人的天堂激情| 成人欧美大片| 亚洲精华国产精华精| 亚洲精品色激情综合| 男人舔女人下体高潮全视频| 国产男靠女视频免费网站| 国内精品美女久久久久久| 真实男女啪啪啪动态图| 看黄色毛片网站| 成人特级黄色片久久久久久久| avwww免费| 午夜福利在线观看吧| 欧美大码av| 动漫黄色视频在线观看| 中文字幕人妻丝袜一区二区| 在线观看免费午夜福利视频| 脱女人内裤的视频| 内射极品少妇av片p| 午夜影院日韩av| 亚洲精品粉嫩美女一区| 国产精品美女特级片免费视频播放器| 国产精品一区二区免费欧美| 国内少妇人妻偷人精品xxx网站| 亚洲无线在线观看| 久久午夜亚洲精品久久| 国产色爽女视频免费观看| 精品国内亚洲2022精品成人| 少妇熟女aⅴ在线视频| 日韩成人在线观看一区二区三区| 日韩欧美三级三区| 91久久精品国产一区二区成人 | 国产欧美日韩精品一区二区| 亚洲va日本ⅴa欧美va伊人久久| 一进一出抽搐gif免费好疼| 久久99热这里只有精品18| 精品久久久久久久末码| 国产av一区在线观看免费| 在线观看免费视频日本深夜| 日韩欧美一区二区三区在线观看| 亚洲av二区三区四区| 中文资源天堂在线| 国内精品一区二区在线观看| 国产精品久久久久久久电影 | 中文字幕人成人乱码亚洲影| 国产真实伦视频高清在线观看 | 制服丝袜大香蕉在线| 成年女人永久免费观看视频| 51午夜福利影视在线观看| 少妇人妻精品综合一区二区 | 久久99热这里只有精品18| 色精品久久人妻99蜜桃| 色视频www国产| 少妇熟女aⅴ在线视频| 国产精品免费一区二区三区在线| 非洲黑人性xxxx精品又粗又长| 国产探花在线观看一区二区| 高清日韩中文字幕在线| 搡老妇女老女人老熟妇| 亚洲av成人精品一区久久| 波多野结衣高清无吗| 久久久久久国产a免费观看| 亚洲国产中文字幕在线视频| 伊人久久大香线蕉亚洲五| 中亚洲国语对白在线视频| 久久久久久久亚洲中文字幕 | 一级黄色大片毛片| 亚洲内射少妇av| 久久久国产成人精品二区| 三级国产精品欧美在线观看| 久久久久久久精品吃奶| 看片在线看免费视频| 999久久久精品免费观看国产| 叶爱在线成人免费视频播放| 天堂影院成人在线观看| 村上凉子中文字幕在线| 久久午夜亚洲精品久久| 无限看片的www在线观看| 88av欧美| 中出人妻视频一区二区| 一级黄色大片毛片| 亚洲av免费在线观看| 欧美成狂野欧美在线观看| 欧美大码av| 国产爱豆传媒在线观看| 嫩草影院精品99| 狂野欧美白嫩少妇大欣赏| 九色国产91popny在线| 波野结衣二区三区在线 | 91在线精品国自产拍蜜月 | 久久亚洲精品不卡| 嫩草影院精品99| 亚洲第一电影网av| www日本在线高清视频| 国产一级毛片七仙女欲春2| 十八禁网站免费在线| 亚洲欧美日韩卡通动漫| 亚洲久久久久久中文字幕| 欧美日韩黄片免| 听说在线观看完整版免费高清| 高清在线国产一区| 日韩精品青青久久久久久| 欧美日韩乱码在线| 法律面前人人平等表现在哪些方面| av天堂中文字幕网| 少妇的逼好多水| 国内精品久久久久久久电影| 免费看光身美女| www.999成人在线观看| 无人区码免费观看不卡| 欧美日韩国产亚洲二区| 久久这里只有精品中国| 国产成人影院久久av| 人人妻人人看人人澡| 热99在线观看视频| 国产精品美女特级片免费视频播放器| 99精品在免费线老司机午夜| 丁香欧美五月| 99riav亚洲国产免费| 国产成人a区在线观看| 老熟妇仑乱视频hdxx| 亚洲国产精品成人综合色| 国产精品永久免费网站| 亚洲av成人不卡在线观看播放网| 免费看十八禁软件| 精品久久久久久久人妻蜜臀av| 国产一区二区在线观看日韩 | 嫩草影视91久久| 久久香蕉国产精品| 婷婷丁香在线五月| 成人av一区二区三区在线看| 老鸭窝网址在线观看| 黑人欧美特级aaaaaa片| 搡老岳熟女国产| 男插女下体视频免费在线播放| 波多野结衣巨乳人妻| 精品国产美女av久久久久小说| 黄色片一级片一级黄色片| 亚洲av二区三区四区| 神马国产精品三级电影在线观看| 国产精品日韩av在线免费观看| 老司机午夜福利在线观看视频| 手机成人av网站| 免费av毛片视频| 亚洲成人久久爱视频| av福利片在线观看| 国产三级黄色录像| 国产av一区在线观看免费| 日韩精品青青久久久久久| 老鸭窝网址在线观看| 99久久99久久久精品蜜桃| 国产精品,欧美在线| 最新在线观看一区二区三区| 成年女人看的毛片在线观看| 国产一区二区在线av高清观看| 欧美激情在线99| 欧美一级a爱片免费观看看| 午夜影院日韩av| 男人和女人高潮做爰伦理| 91九色精品人成在线观看| 日韩欧美 国产精品| 亚洲,欧美精品.| 午夜福利成人在线免费观看| 亚洲av成人不卡在线观看播放网| 美女cb高潮喷水在线观看| 日韩精品中文字幕看吧| 亚洲狠狠婷婷综合久久图片| 欧美大码av| 日韩 欧美 亚洲 中文字幕| 国产成人aa在线观看| av天堂中文字幕网| 舔av片在线| 99热这里只有是精品50| 在线观看一区二区三区| 日韩有码中文字幕| 一个人看的www免费观看视频| 亚洲精品在线美女| 女人高潮潮喷娇喘18禁视频| 欧美bdsm另类| 97超视频在线观看视频| 嫁个100分男人电影在线观看| e午夜精品久久久久久久| 啦啦啦观看免费观看视频高清| 日本熟妇午夜| 好男人电影高清在线观看| 十八禁网站免费在线| 别揉我奶头~嗯~啊~动态视频| 国产午夜精品论理片| 婷婷六月久久综合丁香| 操出白浆在线播放| 老司机在亚洲福利影院| 国产真人三级小视频在线观看| 18美女黄网站色大片免费观看| 亚洲无线观看免费| 在线观看一区二区三区| 日本 欧美在线| 小蜜桃在线观看免费完整版高清| 亚洲成人精品中文字幕电影| 男女视频在线观看网站免费| 欧美日韩国产亚洲二区| www.999成人在线观看| 色尼玛亚洲综合影院| 日韩欧美精品v在线| 美女cb高潮喷水在线观看| 在线播放无遮挡| 久久人人精品亚洲av| 嫁个100分男人电影在线观看| 国产高清三级在线| 日韩 欧美 亚洲 中文字幕| 热99re8久久精品国产| 国产精品影院久久| 中出人妻视频一区二区| av福利片在线观看| 成人特级av手机在线观看| 精品久久久久久久久久免费视频| 久久亚洲精品不卡| 女人被狂操c到高潮| 嫩草影视91久久| 亚洲国产色片| 欧美日韩黄片免| 18禁裸乳无遮挡免费网站照片| 激情在线观看视频在线高清| 欧美一区二区亚洲| 久久九九热精品免费| 久久精品亚洲精品国产色婷小说| 国产成年人精品一区二区| 精品免费久久久久久久清纯| 熟妇人妻久久中文字幕3abv| 少妇裸体淫交视频免费看高清| 亚洲国产欧美网| 国产一区二区在线av高清观看| 日日夜夜操网爽| 国产精品日韩av在线免费观看| 国产伦精品一区二区三区四那| 欧美黑人欧美精品刺激| 在线播放国产精品三级| 亚洲片人在线观看| 亚洲av不卡在线观看| 亚洲欧美日韩无卡精品| 97超级碰碰碰精品色视频在线观看| 波多野结衣高清无吗| 亚洲人与动物交配视频| 国产高潮美女av| 天天躁日日操中文字幕| 日韩国内少妇激情av| 噜噜噜噜噜久久久久久91| 久久人妻av系列| 色噜噜av男人的天堂激情| 成人鲁丝片一二三区免费| 国产精品综合久久久久久久免费| 男插女下体视频免费在线播放| 亚洲va日本ⅴa欧美va伊人久久| 哪里可以看免费的av片| 人妻久久中文字幕网| 国产av在哪里看| 成人性生交大片免费视频hd| 国产亚洲精品久久久久久毛片| 51午夜福利影视在线观看| 搡老熟女国产l中国老女人| 丰满乱子伦码专区| 99精品在免费线老司机午夜| 亚洲国产欧洲综合997久久,| 国产色爽女视频免费观看| 久久精品国产自在天天线| 中出人妻视频一区二区| 丁香欧美五月| 人妻久久中文字幕网| 国产不卡一卡二| 内地一区二区视频在线| 草草在线视频免费看| 国产视频内射| 亚洲 国产 在线| 99精品久久久久人妻精品| 九九在线视频观看精品| 国内毛片毛片毛片毛片毛片| 天堂av国产一区二区熟女人妻| 男女做爰动态图高潮gif福利片| 亚洲欧美精品综合久久99| 2021天堂中文幕一二区在线观| 日韩有码中文字幕| 身体一侧抽搐| 午夜免费激情av| 成人国产综合亚洲| 国产高清有码在线观看视频| 欧美日本视频| 亚洲av一区综合| 一边摸一边抽搐一进一小说| 国内久久婷婷六月综合欲色啪| 国产免费av片在线观看野外av| 亚洲国产中文字幕在线视频| 3wmmmm亚洲av在线观看| 精品无人区乱码1区二区| 中文字幕av成人在线电影| 美女 人体艺术 gogo| 中文字幕av成人在线电影| 黄色片一级片一级黄色片| 日本一本二区三区精品| 18禁国产床啪视频网站| 国产成人av激情在线播放| 亚洲精品粉嫩美女一区| 中文资源天堂在线| 亚洲七黄色美女视频| 国产精品乱码一区二三区的特点| 日韩欧美免费精品| 亚洲国产日韩欧美精品在线观看 | 国产精品久久久久久久久免 | 国产精品99久久久久久久久| 听说在线观看完整版免费高清| 啦啦啦观看免费观看视频高清| 欧美午夜高清在线| 99热这里只有是精品50| 久久亚洲精品不卡| 中出人妻视频一区二区| 97人妻精品一区二区三区麻豆| 女警被强在线播放| 欧美又色又爽又黄视频| 久久久久性生活片| 亚洲精品日韩av片在线观看 | 国产乱人伦免费视频| 欧美xxxx黑人xx丫x性爽| 午夜两性在线视频| 国产亚洲av嫩草精品影院| 国产精品影院久久| 国产又黄又爽又无遮挡在线| 观看美女的网站| 动漫黄色视频在线观看| 亚洲第一欧美日韩一区二区三区| 99在线人妻在线中文字幕| 国产精品国产高清国产av| 在线观看免费午夜福利视频| 九九久久精品国产亚洲av麻豆| 在线国产一区二区在线| 久久久久九九精品影院| 97人妻精品一区二区三区麻豆| 日韩大尺度精品在线看网址| 69人妻影院| 国产精品综合久久久久久久免费| 亚洲不卡免费看| 亚洲狠狠婷婷综合久久图片| 亚洲电影在线观看av| 首页视频小说图片口味搜索| 国产欧美日韩精品亚洲av| 欧美不卡视频在线免费观看| 亚洲人成伊人成综合网2020| 97超视频在线观看视频| 亚洲最大成人中文| bbb黄色大片| 精品国产美女av久久久久小说| 亚洲一区二区三区不卡视频| 真人一进一出gif抽搐免费| avwww免费| 日韩欧美在线乱码| 久久性视频一级片| 国产不卡一卡二| 18+在线观看网站| 久久性视频一级片| 中文字幕人成人乱码亚洲影| 中文字幕高清在线视频| 美女cb高潮喷水在线观看| 五月伊人婷婷丁香| 亚洲中文日韩欧美视频| 亚洲五月婷婷丁香| 久久精品国产自在天天线| 午夜免费男女啪啪视频观看 | 国产三级在线视频| 日韩大尺度精品在线看网址| 母亲3免费完整高清在线观看| 亚洲真实伦在线观看| 国内精品久久久久精免费| 欧美成人a在线观看| 2021天堂中文幕一二区在线观| 国产色婷婷99| 精品久久久久久久末码| 亚洲欧美日韩无卡精品| 国产激情偷乱视频一区二区| 国产成人影院久久av| www日本黄色视频网| 成人精品一区二区免费| 99久久无色码亚洲精品果冻| 三级国产精品欧美在线观看| 桃色一区二区三区在线观看| 免费av毛片视频| 国产欧美日韩精品一区二区| 听说在线观看完整版免费高清| 白带黄色成豆腐渣| 亚洲欧美日韩高清专用| 亚洲美女黄片视频| 亚洲天堂国产精品一区在线|