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

    多源數(shù)據(jù)挖掘下潘謝礦區(qū)深部灰?guī)r水突水預(yù)警研究

    2022-03-18 06:29:14陳永春安士凱徐燕飛
    煤田地質(zhì)與勘探 2022年2期
    關(guān)鍵詞:松散層突水太原

    畢 波 ,陳永春 ,謝 毫 ,安士凱 ,徐燕飛

    (1.煤礦瓦斯治理國(guó)家工程研究中心,安徽 淮南 232001;2.平安煤炭開(kāi)采工程技術(shù)研究院有限責(zé)任公司,安徽 淮南 232001)

    我國(guó)煤礦地質(zhì)條件復(fù)雜、開(kāi)采深度不斷加大導(dǎo)致礦井突水的威脅不斷增加,水害事故造成的經(jīng)濟(jì)損失一直高居首位[1];另一方面,高新技術(shù)與傳統(tǒng)技術(shù)融合,國(guó)內(nèi)各產(chǎn)業(yè)轉(zhuǎn)型升級(jí)不斷深入[2-4]。武強(qiáng)[5]展望了礦井水防控工作創(chuàng)新研究的重大科學(xué)問(wèn)題和關(guān)鍵核心技術(shù)與方法;儲(chǔ)婷婷[6]研究灰?guī)r成分與其中水成分之間的相關(guān)性,揭示了水巖相互作用的機(jī)理;李垣志等[7]、徐星等[8]通過(guò)降低樣本指標(biāo)維度和K 折交叉驗(yàn)證算法,簡(jiǎn)化BP 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu),同時(shí)提高突水水源判斷的精度。楊建等[9]、閆鵬程[10]通過(guò)光譜預(yù)處理對(duì)水樣進(jìn)行熒光光譜分析,確定水源類(lèi)型。劉劍民等[11]基于水化學(xué)分析、王心義等[12]基于熵權(quán)?模糊可變集理論、黃平華等[13]基于多元統(tǒng)計(jì)分析,分別建立了礦井突水水源識(shí)別模型。

    隨著地殼運(yùn)動(dòng)、人為采動(dòng)等因素的影響,深層灰?guī)r各含水層之間一直進(jìn)行著極其復(fù)雜的物質(zhì)、能量交換。導(dǎo)水陷落柱、斷層構(gòu)造都會(huì)形成地下水系統(tǒng)循環(huán)通道[14]。只簡(jiǎn)單研究各含水層水化學(xué)指標(biāo)之間的“表象差異”,無(wú)法準(zhǔn)確判斷復(fù)雜深層灰?guī)r水突水預(yù)兆。

    淮南煤田位于華北板塊東南側(cè),煤田的潘謝礦區(qū)上覆有巨厚的新生界松散層,又受南北兩翼逆沖斷層和斜切斷層的切割,構(gòu)成一個(gè)較為封閉的水文地質(zhì)單元。礦區(qū)自上而下充水水源主要為新生界松散砂層孔隙水、二疊系煤系砂巖裂隙水、底板石炭?二疊系太原組灰?guī)r(簡(jiǎn)稱(chēng)太灰)和奧陶系灰?guī)r(簡(jiǎn)稱(chēng)奧灰)巖溶裂隙水。

    隨著煤礦逐漸向深部開(kāi)采,潘謝礦區(qū)下一階段的主采煤層為A 組煤。對(duì)A 組煤層開(kāi)采威脅最大的因素來(lái)自底板高承壓灰?guī)r水,底板各層灰?guī)r含水層的原始水壓在5 MPa 以上。目前,直接影響煤層安全開(kāi)采的為太原組C2I 組石灰?guī)r含水層,間接影響的為奧陶系和寒武系灰?guī)r(簡(jiǎn)稱(chēng)寒灰)含水層。2017 年5 月25 日,淮南礦業(yè)集團(tuán)潘二礦12123 工作面底板聯(lián)絡(luò)巷發(fā)生隱伏陷落柱突水淹礦,最初出水水質(zhì)特征為太灰水,當(dāng)數(shù)小時(shí)后轉(zhuǎn)變?yōu)閵W灰水水質(zhì)特征時(shí)突水已無(wú)法控制。而往往重大的突水事故都是深部高壓灰?guī)r水以淺層含水層為通道突入礦井發(fā)生的。通過(guò)探明淺層含水層與深層灰?guī)r水的水力聯(lián)系,查找礦井突水風(fēng)險(xiǎn)區(qū)域。

    筆者通過(guò)研究深層灰?guī)r水(奧灰水、寒灰水)對(duì)淺層含水層(上覆石炭–二疊系太原組灰?guī)r、新生界松散層)的影響程度,來(lái)分析發(fā)生奧灰水突水的風(fēng)險(xiǎn),可提高突水預(yù)警的準(zhǔn)確性,縮小預(yù)判時(shí)間,以期為后續(xù)采取防治措施爭(zhēng)取寶貴時(shí)間。

    1 水文觀測(cè)孔水位變化數(shù)據(jù)分析

    1.1 潘謝礦區(qū)的補(bǔ)徑排條件

    潘謝礦區(qū)地質(zhì)條件極其復(fù)雜,構(gòu)造斷層十分發(fā)育。探明斷距大于3 m 的斷層140 多條,見(jiàn)斷點(diǎn)226 個(gè),其中正斷層73 條,占總數(shù)52%;兩淮礦區(qū)發(fā)現(xiàn)近30 個(gè)陷落柱。導(dǎo)水構(gòu)造發(fā)育(例如導(dǎo)水陷落柱)區(qū)域各含水層補(bǔ)給充足容易發(fā)生突水事故,而各含水層通常在采掘工作面前方幾十米甚至幾百米的深處,很難探查,風(fēng)險(xiǎn)較大。如何辨別回采區(qū)域內(nèi)淺層含水層與深層高壓灰?guī)r含水層的水力聯(lián)系,是防治水工作的重點(diǎn)和難點(diǎn)。

    目前,潘謝礦區(qū)正在使用的地面水文觀測(cè)孔182 個(gè),其中新生界松散層含水層70 個(gè)、太灰含水層76 個(gè)、奧灰含水層25 個(gè)、寒灰含水層7 個(gè)、陷落柱探測(cè)孔3 個(gè)、火成巖1 個(gè)。潘謝礦區(qū)地面水文觀測(cè)孔布置及突水后水位變化如圖1 所示。圖1 對(duì)潘二礦突水后主要的奧灰、寒灰鉆孔水位下降高度進(jìn)行標(biāo)注,藍(lán)色標(biāo)注鉆孔受斷層影響水位下降高度相對(duì)較高;紅色標(biāo)注鉆孔為潘三礦的(十二O)鉆孔,因受董崗郢次級(jí)向斜影響,與潘二礦突水點(diǎn)無(wú)水力聯(lián)系。可以推斷,各區(qū)域含水層的補(bǔ)給與斷層密切相關(guān),巖層被拉伸的背斜和富水性良好的斷層都可以作為很好的連通通道。

    表1 為潘二礦突水后3 個(gè)月各地面水文觀測(cè)孔水位變化統(tǒng)計(jì),其中水位高程為礦內(nèi)同一含水層突水前的平均高程。由圖1、表1 可以看出,潘謝礦區(qū)各礦奧灰、寒灰含水層鉆孔突水后水位下降深度隨距突水點(diǎn)的距離(平面距離)增加而減少,說(shuō)明奧灰、寒灰水對(duì)潘二礦突水點(diǎn)進(jìn)行補(bǔ)給。潘謝礦區(qū)各礦的太灰?guī)r溶發(fā)育不均一,具有淺部發(fā)育、斷層帶附近發(fā)育和厚層灰?guī)r發(fā)育的特征。突水后水位也在下降,但水位變化和平均高程差異性很大,并無(wú)規(guī)律可循。說(shuō)明灰?guī)r水也對(duì)突水點(diǎn)進(jìn)行補(bǔ)給,而且是局部通過(guò)奧灰、寒灰間接進(jìn)行補(bǔ)給。新生界松散層含水層地面水文觀測(cè)孔水位幾乎沒(méi)有下降,說(shuō)明新生界松散層含水層對(duì)突水點(diǎn)沒(méi)有明顯補(bǔ)給,且與太原組、奧陶系等灰?guī)r含水層水力聯(lián)系弱。

    表1 突水前后地面水文觀測(cè)孔水位變化統(tǒng)計(jì)Table 1 Statistical table of water level changes in surface hydrological observation holes before and after water inrush

    圖1 潘謝礦區(qū)地面水文觀測(cè)孔布置及突水后水位變化Fig.1 Water level changes of hydrological observation holes in Ordovician and Cambrian after water inrush

    以前普遍認(rèn)為,由于斷層構(gòu)造的影響,潘謝礦區(qū)為一個(gè)較為封閉的水文地質(zhì)單元,地下水以消耗型為主[15-17]。從各含水層水位高程和突水后的水位變化可以看出,礦井主要充水水源為儲(chǔ)量豐富的奧灰和寒灰水。

    1.2 基于數(shù)據(jù)挖掘的突水危險(xiǎn)點(diǎn)識(shí)別

    為確定各灰?guī)r含水層之間的水力聯(lián)系,分析潘二礦突水后各水文觀測(cè)孔水位變化,如圖2 所示。從圖2a 可以看出,潘二礦突水一個(gè)月后在距突水點(diǎn)25 km以?xún)?nèi),奧陶系、寒武系水位下降與距突水點(diǎn)距離成冪函數(shù)(y=94 397x–0.795,R2=0.911 2)變化。說(shuō)明距突水點(diǎn)25 km 以?xún)?nèi)的奧陶系、寒武系、太原組C23 組灰?guī)r和部分太原組C22 灰?guī)r水力聯(lián)系密切。從圖2b 可以看出,距突水點(diǎn)25 km 以外的水文觀測(cè)孔,在突水1 a 后水位也普遍下降,除張集礦外(同期對(duì)灰?guī)r水進(jìn)行抽排)水位下降高度與距突水點(diǎn)距離成反比,而圖2c 顯示突水前1 a 水文觀測(cè)孔水位變化并無(wú)明顯規(guī)律。

    圖2 灰?guī)r水文觀測(cè)孔突水前后水位變化與距突水點(diǎn)距離散點(diǎn)圖Fig.2 Scatter plots of water level changes and distances from water inrush point before and after water inrush in the limestone hydrological observation holes

    為進(jìn)一步確認(rèn)距突水點(diǎn)較遠(yuǎn)礦井各含水層水力聯(lián)系,對(duì)2017 年潘二礦突水5 個(gè)月后顧橋礦、張集礦、謝橋礦的水文觀測(cè)孔水位變化進(jìn)行分析。由于各水文觀測(cè)孔水位變化不大,為消除隨機(jī)誤差,加入2018 年同期的水位變化和突水前水位高程2 類(lèi)數(shù)據(jù)進(jìn)行層次聚類(lèi),來(lái)評(píng)估與奧灰水水位變化相近的鉆孔。為了便于分析,增加變化為0 的55 號(hào)點(diǎn),標(biāo)記為水位變化很小的樣本。標(biāo)號(hào)1?18 為奧陶系、寒武系灰?guī)r鉆孔,19?54 為太原組灰?guī)r鉆孔(19?28 第三組、29?39 第二組、40?54 第一組)、56?89 為松散層鉆孔(56?64 為下含、65?78 為中含、79?89 為上含)、55 號(hào)為零參考點(diǎn)。

    聚類(lèi)分析結(jié)果如圖3 所示,主要分為4 類(lèi)。

    圖3 水文觀測(cè)孔聚類(lèi)分析Fig.3 Cluster analysis of hydraulic observational pores

    第1 組聚類(lèi)結(jié)果為水位在–200 m 左右的太原組灰?guī)r鉆孔,與其他灰?guī)r鉆孔水力聯(lián)系不大。

    第2 組聚類(lèi)結(jié)果為水位變化很平穩(wěn)的鉆孔,其中謝橋礦3 個(gè)奧陶系灰?guī)r孔、2 個(gè)太原組第三組灰?guī)r孔,說(shuō)明謝橋礦太原組第三組灰?guī)r與奧陶系灰?guī)r水力聯(lián)系緊密。包含顧北礦2 個(gè)太原組第一組灰?guī)r、2 個(gè)松散層下含鉆孔,顧北礦太原組灰?guī)r水與松散層下含水有水力聯(lián)系。包含張集礦1 個(gè)西風(fēng)井松散層鉆孔,說(shuō)明張集礦西風(fēng)井松散層與太原組灰?guī)r有一定的水力聯(lián)系。

    第3 組聚類(lèi)結(jié)果為水位變化相對(duì)較大的鉆孔,其中包含顧北礦6 個(gè)奧陶系、寒武系灰?guī)r鉆孔和6 個(gè)太原組灰?guī)r鉆孔,說(shuō)明顧北礦太原組灰?guī)r與奧陶系、寒武系灰?guī)r水力密切。包含顧橋礦全部的3 個(gè)奧陶系、寒武系灰?guī)r鉆孔和2 個(gè)太原組第二組、第三組灰?guī)r鉆孔,并沒(méi)有太原組第一組鉆孔,說(shuō)明顧橋礦奧陶系、寒武系灰?guī)r與太原組第二組、第三組灰?guī)r有水力聯(lián)系。包含謝橋礦3 個(gè)奧陶系、寒武系灰?guī)r鉆孔和2 個(gè)太原組第二組、第三組灰?guī)r鉆孔,并沒(méi)有太原組第一組鉆孔,說(shuō)明謝橋礦與顧橋礦相同,奧陶系、寒武系灰?guī)r與太原組第二組、第三組灰?guī)r有水力聯(lián)系。包含張集礦2 個(gè)奧陶系、寒武系灰?guī)r鉆孔、2 個(gè)太原組灰?guī)r鉆孔和1 個(gè)西風(fēng)井松散層鉆孔。說(shuō)明張集礦奧陶系、寒武系灰?guī)r與太原組灰?guī)r有水力聯(lián)系,且西風(fēng)井松散層與太灰組灰?guī)r有一定的水力聯(lián)系。

    第4 組聚類(lèi)結(jié)果包含為55 號(hào)零參考點(diǎn)為水位變化不大的鉆孔,其中大部分的松散層鉆孔和小部分灰?guī)r鉆孔。

    通過(guò)聚類(lèi)分析,識(shí)別出各礦突水風(fēng)險(xiǎn)區(qū)。謝橋礦太原組灰?guī)r突水風(fēng)險(xiǎn)區(qū)有4 個(gè),其中第三組灰?guī)r孔3 個(gè)、第二組灰?guī)r孔1 個(gè)。顧北礦灰?guī)r突水風(fēng)險(xiǎn)區(qū)有6 個(gè),其中第三組灰?guī)r孔3 個(gè)、第一組灰?guī)r孔3 個(gè),且有下含與太原組第一組灰?guī)r孔水力聯(lián)系密切的鉆孔3 個(gè)。顧橋礦灰?guī)r突水風(fēng)險(xiǎn)區(qū)有2 個(gè),同一位置其第二、第三組灰?guī)r孔各一個(gè)。張集礦灰?guī)r突水風(fēng)險(xiǎn)區(qū)太原組第二、第三組灰?guī)r孔各一個(gè)、西風(fēng)井松散層鉆孔2 個(gè)。

    2 水質(zhì)化驗(yàn)數(shù)據(jù)挖掘

    2.1 水源識(shí)別模型建立

    由于各含水層的賦存年代、巖石巖性、補(bǔ)給水源、埋藏深度等條件不同,會(huì)造成礦井水水化學(xué)成分有所不同,通常采用PiPer 三線圖、特征離子等方法來(lái)分析含水層的水化學(xué)特征,來(lái)判別井下巷道出水來(lái)源。但這種方法的廣度和深度都具有局限性,對(duì)于整個(gè)礦區(qū)多對(duì)礦井、多個(gè)含水層的水源識(shí)別難以實(shí)現(xiàn)。尤其是地質(zhì)條件復(fù)雜區(qū)域,各含水層的水力聯(lián)系復(fù)雜,如潘謝礦區(qū)灰?guī)r水都是同一補(bǔ)給源,只是在長(zhǎng)時(shí)間水巖耦合作用下各含水層的水化學(xué)性質(zhì)發(fā)生了微小變化,水化學(xué)成分差別不大。如此復(fù)雜的地下水賦存條件下,通過(guò)傳統(tǒng)手段難以對(duì)特定區(qū)域含水層進(jìn)行精確探查。

    為識(shí)別各含水層的水化學(xué)特征,收集整理潘謝礦區(qū)各礦井1980?2018 年的7 000 條水質(zhì)化驗(yàn)數(shù)據(jù)。由于研究區(qū)內(nèi)各礦井水質(zhì)化驗(yàn)信息的記錄存在不一致性,且在長(zhǎng)時(shí)間的保存過(guò)程中存在信息偏差導(dǎo)致樣本噪聲大,一般的算法很難準(zhǔn)確對(duì)水源進(jìn)行區(qū)分識(shí)別。隨機(jī)森林算法[18]通過(guò)生成大量獨(dú)立的決策樹(shù)對(duì)結(jié)果進(jìn)行分類(lèi),具有魯棒性[19]和較強(qiáng)的泛化能力,適用于各含水層水質(zhì)化驗(yàn)數(shù)據(jù)的分類(lèi)分析。但隨機(jī)森林算法也有不足,當(dāng)分類(lèi)樣本偏差很大時(shí),分類(lèi)結(jié)果會(huì)偏向于樣本數(shù)據(jù)多的類(lèi)別。通過(guò)對(duì)數(shù)量少的樣本類(lèi)別進(jìn)行合并或人為增加樣本數(shù)量,提高模型的準(zhǔn)確率。對(duì)比選擇不同參數(shù)后的準(zhǔn)確率和各參數(shù)對(duì)各含水層分類(lèi)的貢獻(xiàn)值見(jiàn)表2。其中MDA(Mean Decrease Accuracy)表示把一參數(shù)變隨機(jī)數(shù)后模型準(zhǔn)確率下降的程度,數(shù)值越大說(shuō)明參數(shù)越重要。MDG(Mean Decrease Gini)表示參數(shù)對(duì)決策樹(shù)分支觀測(cè)值異質(zhì)性的影響,數(shù)值越大說(shuō)明參數(shù)越重要。通過(guò)比較模型準(zhǔn)確率,優(yōu)選參數(shù)為Cl?、礦井類(lèi)別、Na++K+、Ca2+、采樣年限、為參考分類(lèi)依據(jù)。對(duì)于奧灰水分類(lèi)貢獻(xiàn)值比較高的常規(guī)水化學(xué)離子為

    表2 水樣各參數(shù)對(duì)于分類(lèi)的貢獻(xiàn)值Table 2 Contribution rate of each parameter of water samples to classification

    2.2 基于隨機(jī)森林的突水危險(xiǎn)點(diǎn)識(shí)別

    對(duì)于建立的水源識(shí)別模型采用10 折交叉驗(yàn)證方法[7]。把樣本按比例隨機(jī)分成10 份,8 份作為訓(xùn)練樣本,2 份作為測(cè)試樣本。計(jì)算10 次后,取計(jì)算結(jié)果組成混淆矩陣見(jiàn)表3,可以近似評(píng)估各含水層分類(lèi)的準(zhǔn)確率。結(jié)果顯示:奧灰水準(zhǔn)確率90%、砂巖水準(zhǔn)確率95%、太灰水準(zhǔn)確率90%、松散層水準(zhǔn)確率92%。

    表3 隨機(jī)森林分類(lèi)結(jié)果混淆矩陣Table 3 Confusion matrix of random forest classification results

    有一部分太灰水與奧灰水連通,導(dǎo)致太灰水分類(lèi)準(zhǔn)確性較低。通過(guò)多次隨機(jī)森林計(jì)算,輸出錯(cuò)誤的樣本,篩查樣本錯(cuò)誤率在70%以上,確定為突水風(fēng)險(xiǎn)區(qū)。對(duì)每一個(gè)礦井的錯(cuò)誤樣本進(jìn)行統(tǒng)計(jì)分析,結(jié)果如下:

    (1) 丁集礦篩查出5 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有10% 砂巖水被錯(cuò)分成松散層,有11% 的松散層水被錯(cuò)分成砂巖水,有4%松散層水被錯(cuò)分成灰?guī)r水,說(shuō)明丁集礦頂板砂巖水主要補(bǔ)給來(lái)自松散層,少部分來(lái)自太灰水。

    (2) 顧北礦篩查出7 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有25%奧灰水被錯(cuò)分成太灰水,說(shuō)明顧北礦奧灰水和太灰水有很強(qiáng)的水力聯(lián)系;有5%的太灰水被錯(cuò)分成松散層水,且有4%松散層水被錯(cuò)分成灰?guī)r水,說(shuō)明顧北礦太灰水松散層有一定水力聯(lián)系。

    (3) 顧橋礦篩查出1 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有25%奧灰水被錯(cuò)分成太灰水,說(shuō)明顧北礦奧灰水和太灰水有很強(qiáng)的水力聯(lián)系。

    (4) 潘北礦篩查出12 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有18%砂巖水被錯(cuò)分成太灰水,說(shuō)明有部分頂板砂巖與太灰有水力聯(lián)系;有15%奧灰水被錯(cuò)分成太灰水,說(shuō)明潘北礦奧灰水和太灰水有一定的水力聯(lián)系。

    (5) 潘二礦篩查出1 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有5%太灰水被錯(cuò)分成奧灰水,由于潘二礦發(fā)生過(guò)奧灰水突水事故,奧灰水與太灰水的水質(zhì)基本相同。

    (6) 潘三礦篩查出3 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有3%砂巖水被錯(cuò)分成松散層水,說(shuō)明潘三礦的頂板砂巖水與松散層水有水力聯(lián)系;有2 個(gè)第三組灰?guī)r水全部被錯(cuò)分成奧灰水,說(shuō)明潘三礦第三組灰?guī)r水與奧灰水水力聯(lián)系很強(qiáng)。

    (7) 潘一礦篩查出3 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有60%松散層水被錯(cuò)分成砂巖水,說(shuō)明潘一礦的松散層水和砂巖水有很強(qiáng)水力聯(lián)系。

    (8) 謝橋礦篩查出8 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有23% 奧灰水被錯(cuò)分成太灰水,有22% 太灰水被錯(cuò)分成奧灰水,說(shuō)明謝橋礦奧灰水和太灰水有很強(qiáng)的水力聯(lián)系。

    (9) 張集礦篩查出4 個(gè)突水風(fēng)險(xiǎn)區(qū)。統(tǒng)計(jì)分析得出有3%太灰水被錯(cuò)分成砂巖水,有4%砂巖水被錯(cuò)分成太灰水,說(shuō)明有一小部分太灰水通過(guò)裂隙補(bǔ)給煤系砂巖水。

    3 快速預(yù)警系統(tǒng)構(gòu)建

    3.1 突水風(fēng)險(xiǎn)區(qū)及敏感指標(biāo)確定

    運(yùn)用隨機(jī)森林算法篩選出44 個(gè)突水風(fēng)險(xiǎn)區(qū),采用水文觀測(cè)孔水位變化分析出19 個(gè)突水風(fēng)險(xiǎn)區(qū),可以相互驗(yàn)證、互補(bǔ)。

    由于各含水層的巖性不同,其地下水水化學(xué)成分也有所不同。找出含水層離子特征,突水危險(xiǎn)源智能檢測(cè)技術(shù)可以為實(shí)施防治水措施提供快速、可靠的依據(jù),極大地避免礦井發(fā)生突水事故和減少突水事故產(chǎn)生的損失。

    為探尋快速識(shí)別灰?guī)r水水源的敏感指標(biāo),分析含水層水化學(xué)垂直分布特征。分別對(duì)奧灰水、太原組3組灰?guī)r水(三灰水)、太原組1 組灰?guī)r水(一灰水)、新生界松散層四含水(下含水)、新生界松散層三含水(中含水)、新生界松散層一、二含水(上含水)進(jìn)行常規(guī)七大離子分布對(duì)比分析[20]。各含水層常規(guī)離子箱線圖如圖4 所示。

    可以看出,奧灰水和三灰水的常規(guī)離子水化學(xué)成分非常接近,說(shuō)明奧灰水和三灰水的水力聯(lián)系緊密。圖4中橙黃線為Cl?含量變化,可以看出與含水層的埋藏深度成正比。綠線為含量變化,與紅線分布成反比。在灰?guī)r含水層SO42?離子與埋藏深度成正比,離子與埋藏深度成反比。因此,可以作為區(qū)分奧灰水、太灰水的敏感指標(biāo)。

    圖4 含水層的常規(guī)離子分布Fig.4 Conventional ion distribution of the aquifer

    TDS、Cl?相關(guān)性很強(qiáng),只需要監(jiān)測(cè)TDS 就可以達(dá)到監(jiān)測(cè)其他離子的效果,且TDS 在線監(jiān)測(cè)技術(shù)已經(jīng)很成熟。結(jié)合建立水源識(shí)別模型時(shí)分析出的離子指標(biāo),綜合考慮確定常規(guī)水化學(xué)監(jiān)測(cè)指標(biāo)為T(mén)DS、同時(shí)需對(duì)水文地質(zhì)條件分析中經(jīng)常采用的溫度、流量、相近水文觀測(cè)孔的水位變化進(jìn)行在線監(jiān)測(cè)。

    3.2 快速預(yù)警系統(tǒng)建立

    對(duì)于傳統(tǒng)的突水危險(xiǎn)預(yù)警,主要依靠現(xiàn)場(chǎng)征兆來(lái)判斷,很難快速、準(zhǔn)確對(duì)突水危險(xiǎn)進(jìn)行評(píng)估、預(yù)測(cè)。如果錯(cuò)過(guò)采取最佳措施的窗口期,很容易造成重大的突水事故。通過(guò)確定各層水源敏感指標(biāo),運(yùn)用連續(xù)在線監(jiān)測(cè)系統(tǒng),對(duì)出水地點(diǎn)進(jìn)行突水預(yù)警[21](圖5)。

    圖5 灰?guī)r水突水危險(xiǎn)識(shí)別與預(yù)警系統(tǒng)Fig.5 Risk identification and early warning system for limestone water inrush

    當(dāng)前期巷道發(fā)現(xiàn)出水時(shí)水量一般很小,根據(jù)附近是否存在斷層等構(gòu)造、是否為灰?guī)r水突水危險(xiǎn)區(qū)來(lái)確定對(duì)應(yīng)的措施。當(dāng)附近存在斷層等構(gòu)造或在判定的灰?guī)r水突水危險(xiǎn)區(qū),就有可能發(fā)生突水事故。需要設(shè)置在線監(jiān)測(cè)設(shè)備,對(duì)出水點(diǎn)進(jìn)行實(shí)時(shí)監(jiān)測(cè),同時(shí)定期進(jìn)行取樣化驗(yàn)分析。當(dāng)出水地點(diǎn)附近不存在斷層等構(gòu)造且不在判別的灰?guī)r水突出危險(xiǎn)區(qū)域,發(fā)生灰?guī)r水突水事件的幾率很小。通過(guò)對(duì)出水溫度、水量、附近地面水文觀測(cè)水位變化等常規(guī)人工監(jiān)測(cè),并定期對(duì)出水點(diǎn)進(jìn)行取樣化驗(yàn)分析。在線監(jiān)測(cè)設(shè)備將巷道出水的TDS、含量、含量、溫度、流量和附近水文觀測(cè)孔的水位變化等數(shù)據(jù)傳輸?shù)降孛婵刂剖?,通過(guò)大數(shù)據(jù)分析技術(shù)對(duì)突水危險(xiǎn)進(jìn)行評(píng)估;取樣化驗(yàn)分析結(jié)果運(yùn)用水源識(shí)別模型進(jìn)行分析,確定出水水源;常規(guī)的人工監(jiān)測(cè)根據(jù)多年防治水經(jīng)驗(yàn)進(jìn)行判別是否有突水危險(xiǎn),3 種分析方法任意一種判別有突水危險(xiǎn)應(yīng)立即采取防治水措施。

    構(gòu)建快速準(zhǔn)確突水預(yù)警體系對(duì)礦井出水點(diǎn)進(jìn)行監(jiān)測(cè),為實(shí)施防治水措施提供快速、可靠的依據(jù),可以極大地避免礦井發(fā)生突水事故和減少突水事故產(chǎn)生的損失。

    4 結(jié)論

    a.通過(guò)分析各礦水質(zhì)化驗(yàn)資料、地面水文觀測(cè)水位變化數(shù)據(jù)及潘二礦突水后各水文觀測(cè)孔水位變化的時(shí)空規(guī)律,得出水文觀測(cè)孔的水位變化數(shù)據(jù)比水位高程數(shù)據(jù)更靈敏,潘謝礦區(qū)深層灰?guī)r水由下向上對(duì)淺層灰?guī)r水進(jìn)行補(bǔ)給。

    b.獲得潘謝礦區(qū)各含水層水化學(xué)空間分布特征,確定各層含水層分類(lèi)顯著因子;基于改進(jìn)隨機(jī)森林算法,建立潘謝礦區(qū)的水源識(shí)別模型。

    c.構(gòu)建快速準(zhǔn)確突水預(yù)警系統(tǒng),通過(guò)在線監(jiān)測(cè)系統(tǒng)對(duì)礦井出水點(diǎn)進(jìn)行監(jiān)測(cè),基于大數(shù)據(jù)分析、水源識(shí)別模型對(duì)灰?guī)r水突水危險(xiǎn)進(jìn)行預(yù)警,為實(shí)施防治水措施提供快速、可靠的依據(jù)。

    猜你喜歡
    松散層突水太原
    礦井突水水源的判別方法
    厚松散層薄基巖工作面導(dǎo)水裂隙帶高度預(yù)計(jì)方法
    太原清廉地圖
    礦井滯后突水機(jī)理探討
    津南區(qū)分層標(biāo)組地質(zhì)鑒別孔松散層取芯施工技術(shù)與效果分析
    除夜太原寒甚
    凍結(jié)松散層中未凍水的變化規(guī)律
    李雅莊礦滯后突水主控因素分析
    大相嶺隧道高壓突水機(jī)理與預(yù)測(cè)分析
    基載比對(duì)厚松散層薄基巖煤層安全開(kāi)采的數(shù)值模擬研究
    河南科技(2014年23期)2014-02-27 14:18:47
    欧美性感艳星| 一本综合久久免费| 夜夜躁狠狠躁天天躁| 精品99又大又爽又粗少妇毛片 | a在线观看视频网站| 亚洲欧美日韩高清专用| 在线免费观看不下载黄p国产 | 97超级碰碰碰精品色视频在线观看| 亚洲国产精品久久男人天堂| 国产一区二区亚洲精品在线观看| a级毛片免费高清观看在线播放| 久久久精品欧美日韩精品| 中国美女看黄片| 欧美xxxx性猛交bbbb| 久久99热6这里只有精品| 日韩欧美国产在线观看| 男女下面进入的视频免费午夜| 人妻丰满熟妇av一区二区三区| 中文字幕av在线有码专区| 久久久久国产精品人妻aⅴ院| 亚洲成人免费电影在线观看| 人人妻人人澡欧美一区二区| 成人毛片a级毛片在线播放| 亚洲精品一区av在线观看| 国产主播在线观看一区二区| 精品久久久久久久久av| 一个人免费在线观看电影| 69av精品久久久久久| 国产欧美日韩一区二区精品| 色av中文字幕| 丁香六月欧美| 精品福利观看| 日韩人妻高清精品专区| 免费电影在线观看免费观看| 淫秽高清视频在线观看| 男人的好看免费观看在线视频| 美女高潮的动态| 看免费av毛片| 人妻制服诱惑在线中文字幕| 久久精品国产亚洲av涩爱 | 久久久久亚洲av毛片大全| 久久99热这里只有精品18| av在线天堂中文字幕| 午夜精品在线福利| 亚洲不卡免费看| 好看av亚洲va欧美ⅴa在| 我的老师免费观看完整版| 3wmmmm亚洲av在线观看| 久久婷婷人人爽人人干人人爱| 有码 亚洲区| 欧美成人性av电影在线观看| 男人和女人高潮做爰伦理| 中文字幕精品亚洲无线码一区| 色视频www国产| 如何舔出高潮| 观看免费一级毛片| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲人成网站在线播放欧美日韩| 99国产极品粉嫩在线观看| 美女高潮喷水抽搐中文字幕| 午夜亚洲福利在线播放| 久9热在线精品视频| 欧美国产日韩亚洲一区| 欧美高清成人免费视频www| 一区二区三区激情视频| 自拍偷自拍亚洲精品老妇| 日本 av在线| 精品久久国产蜜桃| 欧美成人a在线观看| 久久性视频一级片| 久久久国产成人精品二区| 国产成人a区在线观看| 色5月婷婷丁香| 国语自产精品视频在线第100页| 国产精品一区二区性色av| 神马国产精品三级电影在线观看| 久久国产乱子伦精品免费另类| 丰满人妻一区二区三区视频av| 久久久久久国产a免费观看| 久久久久久久精品吃奶| 熟女人妻精品中文字幕| 99热这里只有精品一区| 99久久99久久久精品蜜桃| 日本黄色视频三级网站网址| 精品久久久久久久久亚洲 | 国产私拍福利视频在线观看| 国产一区二区三区在线臀色熟女| 最近中文字幕高清免费大全6 | 床上黄色一级片| 最近中文字幕高清免费大全6 | 天美传媒精品一区二区| 国产亚洲精品久久久com| 国产精品伦人一区二区| 中文字幕久久专区| 久久精品国产亚洲av涩爱 | 亚洲av成人av| 男女之事视频高清在线观看| 69人妻影院| 俺也久久电影网| 一进一出抽搐动态| 男人狂女人下面高潮的视频| 天堂av国产一区二区熟女人妻| 亚洲精品在线美女| 在线天堂最新版资源| 亚洲最大成人中文| 伦理电影大哥的女人| 最好的美女福利视频网| 亚洲18禁久久av| 天天躁日日操中文字幕| 91麻豆av在线| 欧美+日韩+精品| 色在线成人网| 乱人视频在线观看| 国产成人aa在线观看| 夜夜躁狠狠躁天天躁| 国产精品亚洲一级av第二区| 深夜a级毛片| 久久精品夜夜夜夜夜久久蜜豆| 久久精品影院6| 99久久无色码亚洲精品果冻| 一级黄片播放器| 夜夜躁狠狠躁天天躁| 欧美+日韩+精品| 免费av不卡在线播放| 精品人妻视频免费看| 国产av不卡久久| 欧美成人性av电影在线观看| 国产欧美日韩一区二区精品| 亚洲av不卡在线观看| 欧美潮喷喷水| 欧美不卡视频在线免费观看| 12—13女人毛片做爰片一| 精品一区二区免费观看| 欧美性感艳星| 天天躁日日操中文字幕| 桃色一区二区三区在线观看| 观看免费一级毛片| 激情在线观看视频在线高清| 99国产极品粉嫩在线观看| 欧美精品国产亚洲| 亚洲精品影视一区二区三区av| 久久久国产成人精品二区| 午夜视频国产福利| 永久网站在线| 久久亚洲精品不卡| 免费大片18禁| 国产麻豆成人av免费视频| 变态另类丝袜制服| 十八禁人妻一区二区| 亚洲在线观看片| 91字幕亚洲| 天美传媒精品一区二区| 免费在线观看日本一区| av欧美777| 给我免费播放毛片高清在线观看| 国产成+人综合+亚洲专区| 最新中文字幕久久久久| 51国产日韩欧美| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲熟妇中文字幕五十中出| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av| 免费高清视频大片| 女人被狂操c到高潮| 变态另类成人亚洲欧美熟女| 久久精品影院6| 欧美黄色淫秽网站| 国产一级毛片七仙女欲春2| 中文字幕熟女人妻在线| 中文在线观看免费www的网站| 99久久久亚洲精品蜜臀av| 99国产综合亚洲精品| 免费黄网站久久成人精品 | 国产精品日韩av在线免费观看| 国产精品亚洲av一区麻豆| 久久国产精品影院| 亚洲精品粉嫩美女一区| 国产久久久一区二区三区| 1000部很黄的大片| 69人妻影院| 亚洲在线自拍视频| 国产不卡一卡二| 久久久色成人| 亚洲av二区三区四区| 99在线人妻在线中文字幕| 欧美丝袜亚洲另类 | 色噜噜av男人的天堂激情| 国产亚洲av嫩草精品影院| 在线国产一区二区在线| 日韩欧美三级三区| 欧美另类亚洲清纯唯美| 精品人妻1区二区| 成人无遮挡网站| 午夜视频国产福利| 久久99热这里只有精品18| 18禁黄网站禁片午夜丰满| 婷婷亚洲欧美| 我要看日韩黄色一级片| 日本一本二区三区精品| 18禁黄网站禁片免费观看直播| a在线观看视频网站| 91在线观看av| 男女床上黄色一级片免费看| 国产精品久久电影中文字幕| 久久欧美精品欧美久久欧美| 国产精华一区二区三区| 人妻制服诱惑在线中文字幕| or卡值多少钱| 看免费av毛片| 香蕉av资源在线| xxxwww97欧美| 亚洲成人久久性| 99国产精品一区二区蜜桃av| 麻豆一二三区av精品| 久久精品夜夜夜夜夜久久蜜豆| 国产人妻一区二区三区在| 少妇熟女aⅴ在线视频| 99在线人妻在线中文字幕| 欧美日韩亚洲国产一区二区在线观看| 久久欧美精品欧美久久欧美| 日韩亚洲欧美综合| 亚洲成人免费电影在线观看| 俄罗斯特黄特色一大片| 成人欧美大片| 在线免费观看的www视频| 国产69精品久久久久777片| 亚洲 欧美 日韩 在线 免费| 欧美日韩中文字幕国产精品一区二区三区| 99久久精品国产亚洲精品| 日韩中文字幕欧美一区二区| 有码 亚洲区| 成人特级av手机在线观看| 精品午夜福利视频在线观看一区| 国产乱人视频| 久久午夜亚洲精品久久| av天堂中文字幕网| 欧美激情在线99| 波多野结衣巨乳人妻| 亚洲综合色惰| 又爽又黄无遮挡网站| 国产一级毛片七仙女欲春2| 亚洲无线在线观看| 午夜精品在线福利| 欧美一区二区亚洲| 国产欧美日韩一区二区精品| 真实男女啪啪啪动态图| 看免费av毛片| 美女高潮的动态| 国产探花极品一区二区| 欧美日本亚洲视频在线播放| 我要看日韩黄色一级片| 久久久国产成人精品二区| 又黄又爽又免费观看的视频| 草草在线视频免费看| 18禁裸乳无遮挡免费网站照片| 久久久精品大字幕| 91字幕亚洲| 制服丝袜大香蕉在线| 中文字幕免费在线视频6| 国语自产精品视频在线第100页| 亚洲综合色惰| 精品一区二区三区视频在线| 免费一级毛片在线播放高清视频| 99国产极品粉嫩在线观看| 精品免费久久久久久久清纯| 久久精品国产亚洲av天美| 亚洲 欧美 日韩 在线 免费| 亚洲无线在线观看| 亚洲av不卡在线观看| 国产极品精品免费视频能看的| 亚洲aⅴ乱码一区二区在线播放| 看片在线看免费视频| 国产高清激情床上av| 97超级碰碰碰精品色视频在线观看| 欧美三级亚洲精品| 一区二区三区免费毛片| 国产精品1区2区在线观看.| 黄色视频,在线免费观看| 免费观看精品视频网站| 欧美黑人欧美精品刺激| 男女那种视频在线观看| 久久中文看片网| 中文字幕免费在线视频6| 国产男靠女视频免费网站| 别揉我奶头~嗯~啊~动态视频| 日本成人三级电影网站| 啪啪无遮挡十八禁网站| 欧美日韩瑟瑟在线播放| 最近中文字幕高清免费大全6 | 成人美女网站在线观看视频| 他把我摸到了高潮在线观看| 成人鲁丝片一二三区免费| 亚洲 国产 在线| 在现免费观看毛片| 一区二区三区激情视频| 日本一二三区视频观看| 小说图片视频综合网站| 欧美成人a在线观看| 久久精品国产亚洲av涩爱 | 欧美最新免费一区二区三区 | 国产久久久一区二区三区| 精品国产亚洲在线| 性色av乱码一区二区三区2| 精品人妻一区二区三区麻豆 | 日韩av在线大香蕉| 国产色爽女视频免费观看| 欧美区成人在线视频| 国产精品三级大全| 免费看日本二区| 熟女电影av网| 淫秽高清视频在线观看| 老司机福利观看| 不卡一级毛片| 亚洲av一区综合| 直男gayav资源| 一区二区三区高清视频在线| 国产精品久久久久久人妻精品电影| 欧美一级a爱片免费观看看| 国产乱人视频| 一个人观看的视频www高清免费观看| 日韩大尺度精品在线看网址| 男人舔奶头视频| 亚洲五月天丁香| 欧美潮喷喷水| av欧美777| 午夜老司机福利剧场| 色综合婷婷激情| 五月伊人婷婷丁香| 午夜精品久久久久久毛片777| 亚洲va日本ⅴa欧美va伊人久久| 亚洲,欧美精品.| 亚洲成a人片在线一区二区| 日韩精品中文字幕看吧| 国产精品一区二区三区四区久久| 亚洲成a人片在线一区二区| 欧美三级亚洲精品| 一夜夜www| 人妻丰满熟妇av一区二区三区| 国产三级在线视频| 欧美xxxx黑人xx丫x性爽| 免费一级毛片在线播放高清视频| 久久精品夜夜夜夜夜久久蜜豆| 国产美女午夜福利| 国产黄片美女视频| 变态另类丝袜制服| 老司机福利观看| 成人av一区二区三区在线看| 国产精品精品国产色婷婷| 日本在线视频免费播放| av天堂在线播放| 久久久久九九精品影院| 亚洲无线观看免费| 国产精品一及| 最新在线观看一区二区三区| 身体一侧抽搐| 欧美zozozo另类| 日日夜夜操网爽| 男女下面进入的视频免费午夜| 国产真实乱freesex| 国产精品乱码一区二三区的特点| 精品午夜福利在线看| 欧美色欧美亚洲另类二区| 亚洲人成网站在线播放欧美日韩| 听说在线观看完整版免费高清| 亚洲不卡免费看| 97碰自拍视频| 日韩欧美精品免费久久 | 亚洲无线在线观看| 国内久久婷婷六月综合欲色啪| 久久人人精品亚洲av| 桃红色精品国产亚洲av| 国产精品亚洲美女久久久| 亚州av有码| 美女高潮的动态| 不卡一级毛片| 69人妻影院| 欧美潮喷喷水| 日韩欧美一区二区三区在线观看| 99热这里只有是精品在线观看 | 老司机午夜福利在线观看视频| 日本黄色视频三级网站网址| 亚洲欧美日韩卡通动漫| 国产在线男女| 欧美午夜高清在线| 首页视频小说图片口味搜索| 全区人妻精品视频| 亚洲av一区综合| 国产精品电影一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕精品亚洲无线码一区| 欧美日韩综合久久久久久 | 国产精品,欧美在线| 美女高潮的动态| 两个人的视频大全免费| 老鸭窝网址在线观看| 韩国av一区二区三区四区| 亚州av有码| 亚洲人成网站在线播| 此物有八面人人有两片| 高清毛片免费观看视频网站| 久久精品国产亚洲av涩爱 | 成人av在线播放网站| 国产黄片美女视频| 午夜福利免费观看在线| 国产在视频线在精品| 最好的美女福利视频网| 亚洲欧美日韩高清专用| 国产白丝娇喘喷水9色精品| or卡值多少钱| 亚洲在线观看片| 久久亚洲精品不卡| 久久精品国产99精品国产亚洲性色| 大型黄色视频在线免费观看| 亚洲五月婷婷丁香| 看十八女毛片水多多多| 亚洲性夜色夜夜综合| 精品久久久久久久久久免费视频| www日本黄色视频网| 国产人妻一区二区三区在| 看十八女毛片水多多多| 免费大片18禁| 99久久精品一区二区三区| 久久99热这里只有精品18| 国产精品亚洲av一区麻豆| 久久久久久九九精品二区国产| 最后的刺客免费高清国语| 国产精品乱码一区二三区的特点| a级毛片免费高清观看在线播放| 十八禁国产超污无遮挡网站| 亚洲国产精品sss在线观看| 99在线视频只有这里精品首页| 女同久久另类99精品国产91| 女人十人毛片免费观看3o分钟| 日韩大尺度精品在线看网址| 日本 av在线| 欧美丝袜亚洲另类 | 欧美日韩黄片免| 日韩高清综合在线| 黄色配什么色好看| 国产精品自产拍在线观看55亚洲| 99热这里只有是精品50| 国产不卡一卡二| 又爽又黄无遮挡网站| 国产三级中文精品| 午夜精品在线福利| 女同久久另类99精品国产91| 亚洲av熟女| 搡老熟女国产l中国老女人| 国产一区二区在线av高清观看| 亚洲三级黄色毛片| 舔av片在线| 又黄又爽又免费观看的视频| 国产精华一区二区三区| 日韩欧美国产一区二区入口| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 国产欧美日韩精品一区二区| 天美传媒精品一区二区| 国产亚洲精品久久久久久毛片| 真人做人爱边吃奶动态| 精品一区二区三区人妻视频| 女人十人毛片免费观看3o分钟| 内射极品少妇av片p| 中文字幕熟女人妻在线| 免费在线观看日本一区| 丰满人妻一区二区三区视频av| 国产在视频线在精品| 他把我摸到了高潮在线观看| 免费看日本二区| 国产野战对白在线观看| 久久久久性生活片| 成人国产一区最新在线观看| 日韩欧美精品免费久久 | 色综合亚洲欧美另类图片| 欧美黄色淫秽网站| 成人av一区二区三区在线看| 欧美成人一区二区免费高清观看| 亚洲,欧美,日韩| 色综合站精品国产| 国产精品亚洲美女久久久| 最近最新免费中文字幕在线| 国产高清视频在线播放一区| 国产人妻一区二区三区在| av在线蜜桃| 国产精品嫩草影院av在线观看 | 麻豆一二三区av精品| 欧美高清成人免费视频www| 精品国产三级普通话版| 国产视频内射| 国产乱人视频| 国产成人影院久久av| 亚洲国产精品合色在线| 嫩草影院精品99| 乱码一卡2卡4卡精品| 欧美在线黄色| 97超级碰碰碰精品色视频在线观看| 日韩人妻高清精品专区| 国产欧美日韩一区二区精品| 亚洲 欧美 日韩 在线 免费| 色精品久久人妻99蜜桃| 人妻制服诱惑在线中文字幕| 搡老妇女老女人老熟妇| 天堂av国产一区二区熟女人妻| 久久这里只有精品中国| 如何舔出高潮| 2021天堂中文幕一二区在线观| 51国产日韩欧美| 动漫黄色视频在线观看| 十八禁人妻一区二区| 国产老妇女一区| 亚洲七黄色美女视频| 伊人久久精品亚洲午夜| 性色av乱码一区二区三区2| 18+在线观看网站| 成人特级av手机在线观看| 免费在线观看日本一区| 欧美午夜高清在线| 亚洲av电影不卡..在线观看| 亚洲自拍偷在线| 中文字幕精品亚洲无线码一区| 成年女人看的毛片在线观看| 亚洲精品影视一区二区三区av| 少妇高潮的动态图| 少妇裸体淫交视频免费看高清| 国产熟女xx| 18禁在线播放成人免费| 精品久久久久久久人妻蜜臀av| 精品久久国产蜜桃| 18+在线观看网站| 偷拍熟女少妇极品色| 韩国av一区二区三区四区| 琪琪午夜伦伦电影理论片6080| 亚洲人成网站高清观看| 久久99热这里只有精品18| 国产一区二区三区在线臀色熟女| 精品久久久久久久久久免费视频| 美女高潮喷水抽搐中文字幕| 最近最新中文字幕大全电影3| 九九久久精品国产亚洲av麻豆| 在线播放国产精品三级| 高清日韩中文字幕在线| 国产精品综合久久久久久久免费| 毛片女人毛片| 无人区码免费观看不卡| 国产高潮美女av| 很黄的视频免费| 无遮挡黄片免费观看| 欧美日韩亚洲国产一区二区在线观看| 久久精品国产自在天天线| 欧美黄色淫秽网站| 3wmmmm亚洲av在线观看| 人妻久久中文字幕网| 国内揄拍国产精品人妻在线| 直男gayav资源| 国产免费av片在线观看野外av| 很黄的视频免费| 亚洲乱码一区二区免费版| 丝袜美腿在线中文| 成人鲁丝片一二三区免费| 久久热精品热| 亚洲成人久久爱视频| 大型黄色视频在线免费观看| 精品熟女少妇八av免费久了| 国产高清视频在线播放一区| 最近在线观看免费完整版| 在线a可以看的网站| 亚洲黑人精品在线| 波野结衣二区三区在线| 在线观看美女被高潮喷水网站 | a级一级毛片免费在线观看| 欧美精品国产亚洲| 国产大屁股一区二区在线视频| 久久久久久久久久黄片| 99热这里只有是精品50| 国产黄a三级三级三级人| 赤兔流量卡办理| 在线观看美女被高潮喷水网站 | 国产免费av片在线观看野外av| 又黄又爽又刺激的免费视频.| 国产一区二区亚洲精品在线观看| 国产午夜精品论理片| 99久久无色码亚洲精品果冻| 亚洲中文字幕日韩| 熟女人妻精品中文字幕| 美女大奶头视频| 欧洲精品卡2卡3卡4卡5卡区| 香蕉av资源在线| 午夜精品一区二区三区免费看| 欧美黑人欧美精品刺激| 亚洲三级黄色毛片| 亚洲国产欧洲综合997久久,| 免费人成在线观看视频色| 变态另类丝袜制服| 免费av观看视频| 观看美女的网站| 看十八女毛片水多多多| 日韩高清综合在线| 国产高清视频在线播放一区| 日本撒尿小便嘘嘘汇集6| 欧美在线黄色| 国产 一区 欧美 日韩| 久久精品综合一区二区三区| 在线十欧美十亚洲十日本专区| 午夜福利视频1000在线观看| 又粗又爽又猛毛片免费看| 欧美色欧美亚洲另类二区| 搡女人真爽免费视频火全软件 | 国产伦精品一区二区三区视频9| 精品一区二区三区av网在线观看| 九九热线精品视视频播放| 国内精品一区二区在线观看| 亚洲欧美清纯卡通| 成人国产一区最新在线观看| x7x7x7水蜜桃| 久久精品夜夜夜夜夜久久蜜豆|