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

    基于SSA優(yōu)化的GA-BP神經(jīng)網(wǎng)絡煤層底板突水預測模型與應用

    2022-01-04 12:10:56尹會永周鑫龍張歷峰王明麗
    煤田地質(zhì)與勘探 2021年6期
    關鍵詞:突水危險性底板

    尹會永,周鑫龍,郎 寧,張歷峰,王明麗,吳 燾,李 鑫

    基于SSA優(yōu)化的GA-BP神經(jīng)網(wǎng)絡煤層底板突水預測模型與應用

    尹會永1,周鑫龍1,郎 寧1,張歷峰2,王明麗3,吳 燾2,李 鑫4

    (1. 山東科技大學 地球科學與工程學院,山東 青島 266590;2. 棗莊礦業(yè)(集團)有限責任公司濱湖煤礦,山東 棗莊 277599;3. 山東省青島生態(tài)環(huán)境監(jiān)測中心,山東 青島 266003;4. 寧夏煤炭科學技術研究所有限公司,寧夏 銀川 750021)

    隨煤層開采深度的不斷增加,煤礦生產(chǎn)過程中面臨著復雜的突水機理和多變的突水主控因素,且各因素間相互聯(lián)系的不確定性,使底板突水預測的難度不斷增加。為準確預測底板突水危險性,針對底板突水的小樣本、非線性問題,首先利用遺傳算法(Genetic Algorithm,GA)將網(wǎng)絡隨機賦值的初始權值和閾值初次優(yōu)化,再選取搜索能力強、穩(wěn)定性較好的麻雀搜索算法(Sparrow Search Algorithm,SSA)對權值和閾值進行二次尋優(yōu),從而建立SSA-GA-BP神經(jīng)網(wǎng)絡底板突水預測模型。分析整理山東省濱湖煤礦地質(zhì)及水文地質(zhì)資料,選取含水層水壓、含水層厚度、隔水層厚度、斷層密度、斷層分維值、滲透系數(shù)、單位涌水量、底板破壞深度共8個因素,作為預測底板突水的主控因素,繪制各主控因素3D映射投影曲面圖;利用Surfer軟件中的克里金插值法提取50個數(shù)據(jù)點作為模型的輸入樣本(分為訓練集40個,測試集10個),對模型進行訓練學習,訓練誤差精度達到要求后,對濱湖煤礦3個未開采工作面的12個數(shù)據(jù)點進行突水危險性預測。為了驗證所建模型的準確性,利用BP、GA-BP、SSA-GA-BP這3種模型對測試集進行預測;為避免模型僅與BP網(wǎng)絡預測對比的片面性,同時選取以熵權法確定權重的模糊綜合評判法對測試集進行預測;將各網(wǎng)絡模型及方法的預測結果與實際值進行對比分析。結果表明:基于SSA優(yōu)化的GA-BP神經(jīng)網(wǎng)絡模型突水預測誤差較小,預測結果準確率更高,為礦井水害預測預報提供了科學的評價方法和理論依據(jù)。

    底板突水預測;麻雀搜索算法;遺傳算法;BP神經(jīng)網(wǎng)絡;熵權法;模糊綜合評判

    我國煤炭資源的開采強度及范圍不斷增大,煤田地質(zhì)及構造條件類型多樣。隨著開采深度的不斷增加,加之礦井復雜的突水機理與隱蔽的致災因素,導致煤層底板突水預測的難度不斷增加,礦井突水事故頻發(fā),給社會經(jīng)濟和人民安全造成嚴重威脅[1-3]。因此,使用準確有效的方法對底板突水危險性進行及時預測,有利于提出科學的風險預防決策,煤礦生產(chǎn)安全性進一步提升[4]。

    目前,傳統(tǒng)預測煤層底板突水危險性的方法主要有突水系數(shù)法、“下三帶”法、模糊證據(jù)理論法[5]、層次分析法[6]、多元信息復合法[7]等,為礦井的安全生產(chǎn)、水害防治提供了大量有效經(jīng)驗和參考。近年來,對于煤層底板突水預測出現(xiàn)了一些新方法,如集對分析法[8]、多決策樹法[9]、熵權耦合法[10]、Logistic回歸分析法[11]、突水概念場地模型法[12]、支持向量機法[13]、梯形模糊數(shù)法[14]、多模型融合評價法[15]等。上述方法各有特點,針對不同類型煤礦的突水預測也取得了較好的效果。但這些方法對現(xiàn)場地質(zhì)資料和參數(shù)的要求較高,且部分指標權重主觀性較強,因此,在應用時具有一定的局限性。面對上述問題,人工神經(jīng)網(wǎng)絡作為一種非線性動力系統(tǒng)具有實時適應性、自學習、并行性等功能,對復雜突水因素的處理能力較強,給底板突水預測提供了一種新思路,于小鴿等[16]提出了BP神經(jīng)網(wǎng)絡模型,該類方法所需樣本較少,對參數(shù)要求較低,但數(shù)據(jù)間聯(lián)系較弱,預測值與實際值有一定差距;陳建平等[17]建立了基于卷積神經(jīng)網(wǎng)絡預測模型,對數(shù)據(jù)的特征規(guī)律進行了有效提取,預測精度進一步提高,但由于隨機化的初始權值和閾值,使網(wǎng)絡容易陷入局部最優(yōu)解;施龍青等[18]利用灰狼算法改進網(wǎng)絡隨機賦值的初始權值、閾值,建立GWO-Elman預測模型,網(wǎng)絡搜索精度提升較大,但迭代收斂較慢。因此,可考慮用其他算法對神經(jīng)網(wǎng)絡進行優(yōu)化。薛建凱[19]于2020年提出一種新型群體智能優(yōu)化算法——麻雀搜索算法(Sparrow Search Algorithm,SSA),具有收斂速度快,搜索精度高,穩(wěn)定性好,魯棒性強等特點。針對煤層底板突水復雜非線性、多因素疊加問題,該算法有較好的適應性,對底板突水預測來說是一種新嘗試。

    筆者以山東棗莊市濱湖煤礦為研究背景,利用遺傳算法(Genetic Algorithm,GA)將網(wǎng)絡隨機化的初始權值、閾值初次優(yōu)化,再利用麻雀搜索算法對權值、閾值二次尋優(yōu),從而建立SSA-GA-BP網(wǎng)絡底板突水預測模型,進一步提高GA-BP網(wǎng)絡的非線性動態(tài)映射能力和突水預測的準確性,以期為礦井底板水害防治提供理論依據(jù)。

    1 研究區(qū)概況

    濱湖煤礦位于山東省棗莊市滕北礦區(qū)的中西部,區(qū)內(nèi)主要可采煤層有太原組12下、16煤層。其中,12下煤層全區(qū)分布,平均厚度1.27 m;研究區(qū)為162采區(qū),面積約6.37 km2,16煤為該采區(qū)主要開采煤層,平均厚度1.23 m。162采區(qū)位于藤縣背斜的北翼,屬單斜構造向西北傾斜,區(qū)內(nèi)揭露1 m以上斷層25條,構造復雜程度為中等。

    2 底板突水主控因素

    煤層底板突水是一種多因素復合疊加的非線性動態(tài)過程,各影響因素間有著密切的聯(lián)系[20]。本文以井田地質(zhì)、水文地質(zhì)、構造等資料為基礎,分析得到濱湖煤礦162采區(qū)煤層底板的主要充水水源為奧陶紀灰?guī)r(簡稱奧灰)含水層,充水通道主要有采動裂隙、斷層、裂隙帶,并根據(jù)以往煤層底板突水特征,選取含水層水壓(1)、含水層厚度(2)、隔水層厚度(3)、斷層密度(4)、斷層分維值(5)、滲透系數(shù)(6)、單位涌水量(7)、底板破壞深度(8)共8個因素作為研究煤層底板突水的具體依據(jù)。根據(jù)其具體資料,利用Origin軟件繪制各因素的3D映射投影曲面圖(圖1)。

    2.1 含水層水壓(A1)

    研究區(qū)16煤底板埋深為490~570 m。目前,奧灰水位高程為–168.53 m(2020年10月28日),煤層底板所承受水壓力較大,可以達到3.29~4.39 MPa,具有奧灰突水的危險性。含水層水壓3D映射投影曲面如圖1a所示。

    2.2 含水層厚度(A2)

    根據(jù)礦井實際突水情況及鉆孔揭露資料,引起研究區(qū)煤層底板突水的主要含水層為奧灰頂部上馬家溝組含水層,該段含水層淺部48~90 m富水性較強,與煤層底板具有直接水力聯(lián)系,因此,將該段作為奧灰含水層的等效厚度(圖1b)。

    2.3 隔水層厚度(A3)

    研究區(qū)16煤至奧灰間主要由泥巖、煤層、砂質(zhì)泥巖、細砂巖、粉砂巖及灰?guī)r組成,間距44.2~ 55.7 m,平均50.9 m,依據(jù)鉆孔資料利用Origin軟件繪制隔水層厚度3D映射投影曲面圖(圖1c)。

    2.4 斷層密度(A4)

    斷層密度指單位面積內(nèi)斷層跡線的條數(shù)。一般在構造復雜、密度較大的地方發(fā)育較多的儲水空間和導水通道,更容易發(fā)生突水[21]。斷層密度的計算公式為:

    式中:為網(wǎng)格中斷層(或褶皺軸部)跡線的條數(shù);為網(wǎng)格面積,km2。

    2.5 斷層分行維數(shù)值(A5)

    斷層分行維數(shù)值是描述斷裂構造復雜程度的合理參數(shù)。某區(qū)段的分行維數(shù)值越大,則表示該區(qū)斷層分支較多、密度及長度較大,對周邊巖體的剪切破壞也越嚴重,底板突水的危險性越大。定義集合()的相似維為:

    2.6 滲透系數(shù)(A6)

    滲透系數(shù)表明巖層的透水性強弱,是底板突水預測的重要因素。根據(jù)濱湖煤礦的放水試驗資料,利用Aquifer Test軟件中Theis配線法求取研究區(qū)不同范圍內(nèi)的滲透系數(shù)為0.10~0.29 m/d (圖1f)。

    2.7 單位涌水量(A7)

    單位涌水量是體現(xiàn)含水層排水能力的有效指標。根據(jù)濱湖煤礦奧灰水文鉆孔資料,研究區(qū)單位涌水量為0.05~0.27 L/(s?m),其3D映射投影曲面如圖1g所示。

    2.8 底板破壞深度(A8)

    底板破壞深度是受人工影響的充水通道,根據(jù)《建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規(guī)范》計算公式:

    3 模型的建立與應用

    3.1 麻雀搜索算法原理

    麻雀搜索算法(Sparrow Search Algorithm, SSA)是通過模擬自然界中麻雀生存機制進而提出的一種群體優(yōu)化算法[22]。麻雀群體可以分為發(fā)現(xiàn)者和跟隨者,發(fā)現(xiàn)者負責搜尋種群的食物,搜索能力強,屬于優(yōu)勢個體,跟隨者則跟隨發(fā)現(xiàn)者來搜尋食物。同時,為了提高自身的適應度,部分跟隨者會對發(fā)現(xiàn)者進行監(jiān)視,以便進行搜索和爭奪食物。此外,當天敵臨近時,種群會及時進行反捕食行為。發(fā)現(xiàn)者的位置更新公式如下:

    圖1 底板突水主控因素等值線及投影曲面

    跟隨者的位置更新公式如下:

    當受到天敵威脅時,偵查麻雀發(fā)出預警信號,種群迅速做出反捕食行為。偵查麻雀數(shù)量一般為種群的10%~20%,位置更新公式如下:

    3.2 神經(jīng)網(wǎng)絡模型優(yōu)化流程

    為避免出現(xiàn)網(wǎng)絡對初始權值和閾值隨機賦值而造成局部最優(yōu)解等缺陷,引入遺傳算法對其進行優(yōu)化,降低局部最優(yōu)概率,預測精度有一定的提高,但網(wǎng)絡全局尋優(yōu)能力尚且不足。針對上述問題,對已初次優(yōu)化的權值、閾值編碼后代入麻雀搜索算法中進行二次尋優(yōu),提高GA-BP網(wǎng)絡的尋優(yōu)能力和迭代收斂速度,優(yōu)化網(wǎng)絡的非線性動態(tài)映射能力,判別精度進一步提高。SSA優(yōu)化GA-BP神經(jīng)網(wǎng)絡算法流程如圖2所示。

    圖2 SSA優(yōu)化GA-BP神經(jīng)網(wǎng)絡算法流程

    3.3 實例分析

    3.3.1 數(shù)據(jù)來源及預處理

    整理分析棗莊市濱湖煤礦162采區(qū)的地質(zhì)勘探報告和鉆孔資料,根據(jù)突水系數(shù)法及高密度電法勘探結果對采區(qū)底板突水危險性的評價結果,將煤層底板的突水危險性劃分為3個不同等級,分別為安全區(qū)—0區(qū)、低風險區(qū)—1區(qū)和高風險區(qū)—2區(qū),利用Surfer軟件中的克里金插值法獲得采區(qū)50個數(shù)據(jù)點的突水危險性實際評價等級值(表1)。選取上述8種底板突水主控因素作為突水預測的指標依據(jù),利用Surfer軟件中的克里金插值法獲得各個數(shù)據(jù)點的主控因素數(shù)據(jù)值作為樣本數(shù)據(jù)集。由于樣本數(shù)據(jù)各指標值的大小差別較大,為避免低值數(shù)據(jù)被高值數(shù)據(jù)覆蓋,且防止節(jié)點輸出進入飽和區(qū),采用Matlab中的mapminmax函數(shù)將輸入數(shù)據(jù)進行歸一化處理,使其取值范圍為[0,1]。

    表1 突水危險性評價樣本數(shù)據(jù)集

    3.3.2 模型建立與驗證

    利用遺傳算法將網(wǎng)絡隨機化的權值、閾值初次優(yōu)化,具體的參數(shù)設定為:種群數(shù)目為50,最大迭代次數(shù)為80,變異、交叉率分別為0.1、0.2。繼而使用麻雀搜索算法對已優(yōu)化的網(wǎng)絡權值和閾值進行二次尋優(yōu),麻雀搜索算法參數(shù)設定為:種群數(shù)目為30,最大迭代次數(shù)為80,發(fā)現(xiàn)者和偵察者個數(shù)各取種群數(shù)目的20%,預警值2=0.8。

    為驗證模型是否準確有效,選用前40個樣本點數(shù)據(jù)作為SSA-GA-BP神經(jīng)網(wǎng)絡模型的訓練集,剩余10個樣本點數(shù)據(jù)作為模型測試數(shù)據(jù),網(wǎng)絡訓練次數(shù)達25次時,模型預測誤差可以達到7.938 1í10–9(圖3),低于目標誤差值1í10–8,模型對訓練集的預測結果準確率為100%,符合要求,因此,模型可用于濱湖煤礦162采區(qū)的底板突水預測。

    3.3.3 模型應用

    濱湖煤礦計劃在2021年開采162采區(qū)的3個工作面16205、16207和16208工作面,工作面位置如圖4所示。根據(jù)各工作面掘進期間所獲得的地質(zhì)數(shù)據(jù)資料,預測數(shù)據(jù)集由每個工作面各4個數(shù)據(jù)點組成,各點數(shù)據(jù)見表2,其中K1—K4屬于16205工作面,K5—K8屬于16207工作面,K9—K12屬于16208工作面。采用已訓練好且誤差精度達到要求的SSA-GA-BP預測模型對上述3個工作面的突水危險性進行預測,預測結果見表2。

    圖3 誤差收斂曲線

    圖4 濱湖煤礦162采區(qū)工作面分布

    3.4 其他方法對比

    3.4.1 BP神經(jīng)網(wǎng)絡方法對比

    為說明本文所建立的SSA-GA-BP網(wǎng)絡突水預測模型與以往神經(jīng)網(wǎng)絡突水預測模型的不同之處,用訓練集對BP網(wǎng)絡、GA-BP網(wǎng)絡進行訓練后,分別用于對測試集底板進行突水預測,網(wǎng)絡及算法參數(shù)設置與上文相同。測試集實際值與各方法預測值的對比如圖5所示,各方法的預測誤差如圖6所示。

    根據(jù)所得出的測試結果,預測精度最高的為SSA-GA-BP網(wǎng)絡,除對41號樣本預測誤差(–0.29)較大外,其余樣本誤差均較小,預測結果準確率為90%;GA-BP網(wǎng)絡對41號、47號、48號樣本的預測誤差較大,預測精度明顯較低,誤差精度較前者平均相差2個數(shù)量級,預測結果準確率為70%;而BP網(wǎng)絡的預測結果最差,對41號、42號、46號、47號和49號樣本的預測誤差均較大,預測結果準確率為50%。3種方法對測試集的預測結果及預測誤差詳見表3。從整體效果來看,上述方法的預測精度從高到低依次為:SSA-GA-BP、GA-BP、BP。

    圖7反映了隨迭代次數(shù)增加種群適應度的進化曲線,與GA算法相比,SSA算法在第5次迭代后種群適應度值的進化程度較為明顯,相同迭代次數(shù)下適應度遞減趨勢較快,經(jīng)過36次的迭代便達到了適應度的全局最優(yōu)值,且最終適應度值也得到了明顯的優(yōu)化,因此,本文所提出的SSA-GA優(yōu)化算法對BP神經(jīng)網(wǎng)絡的收斂速度及預測精度提升較大。

    表2 突水危險性評價預測數(shù)據(jù)集

    圖5 測試集預測值與實際值對比

    圖6 測試集預測誤差對比

    表3 各突水危險性方法預測評價等級值及預測誤差

    圖7 適應度進化曲線

    3.4.2 模糊綜合評判方法對比

    引入模糊綜合評判方法進行測試集的突水危險性評價,由于各因素對突水危險性影響的重要程度不同,因此,各因素權重的確定采用熵權法。

    熵權法依據(jù)不同因素間變異程度強弱從而對權重進行客觀賦值,因素的熵值越小,則該因素所包含的信息越重要,對應的權重就越大[24]。根據(jù)前文所述的8種底板突水主控因素,利用熵權法計算各因素的權重,公式如下。

    表4 模糊綜合評價等級分布

    結合樣本數(shù)據(jù),利用公式(11)計算得出各樣本所對應的突水危險等級模糊評價值表(表5),建立模糊綜合評判模型。經(jīng)過與實際值的對比,模糊綜合評價法對測試樣本中41號和47號的預測值分別為0和1,與實際值不符,模型預測準確率為80%;對預測樣本中K7、K10和K11的預測值分別為2、1和0,預測結果與SSA-GA-BP神經(jīng)網(wǎng)絡模型對應樣本預測值1、0和1不同。因后者模型對測試樣本的預測準確率更高,3個未開采工作面突水危險性最終預測結果以SSA-GA-BP網(wǎng)絡模型為準。

    4 結論

    a. 分析了影響煤層底板突水的8種主控因素,選取穩(wěn)定性好、全局搜索能力強的麻雀搜索算法(SSA)對GA-BP網(wǎng)絡權值與閾值進行二次尋優(yōu),建立SSA-GA-BP神經(jīng)網(wǎng)絡底板突水預測模型,經(jīng)測試集數(shù)據(jù)驗證,模型預測值與實際值較符合,預測準確率為90%。

    b. 利用SSA-GA-BP模型對山東棗莊濱湖煤礦162采區(qū)未開采的3個工作面底板突水危險性進行預測,結果顯示:16205工作面處于安全區(qū)無突水風險;16207工作面的K7樣本點處于突水低風險區(qū),有突水風險的可能;16208工作面的K9和K11樣本點分別處于危險區(qū)和低風險區(qū),其底板突水危險性較大,開采過程中需進行重點防治。

    表5 突水危險性評價等級綜合對比

    c. 建立GA-BP、BP預測模型,并對測試集進行預測,將預測結果與實際值進行對比,其預測準確率分別為70%、50%,預測精度較低。以熵權法確定主控因素權重的模糊綜合評判模型用于測試集數(shù)據(jù)驗證,結果顯示該模型預測準確率為80%。

    [1] 武強. 我國礦井水防控與資源化利用的研究進展、問題和展望[J]. 煤炭學報,2014,39(5):795–805.

    WU Qiang. Progress,problems and prospects of prevention and control technology of mine water and reutilization in China[J]. Journal of China Coal Society,2014,39(5):795–805.

    [2] 王皓,董書寧,喬偉,等. 礦井水害防控遠程服務云平臺構建與應用[J]. 煤田地質(zhì)與勘探,2021,49(1):208–216.

    WANG Hao,DONG Shuning,QIAO Wei,et al. Construction and application of remote service cloud platform for mine water hazard prevention and control[J]. Coal Geology & Exploration,2021,49(1):208–216.

    [3] 靳德武. 我國煤礦水害防治技術新進展及其方法論思考[J]. 煤炭科學技術,2017,45(5):141–147.

    JIN Dewu. New development of water disaster prevention and control technology in China coal mine and consideration on methodology[J]. Coal Science and Technology,2017,45(5):141–147.

    [4] 虎維岳,田干. 我國煤礦水害類型及其防治對策[J]. 煤炭科學技術,2010,38(1):92–96.

    HU Weiyue,TIAN Gan. Mine water disaster type and prevention and control counter measures in China[J]. Coal Science and Technology,2010,38(1):92–96.

    [5] 肖建于,童敏明,姜春露. 基于模糊證據(jù)理論的煤層底板突水量預測[J]. 煤炭學報,2012,37(增刊1):131–137.

    XIAO Jianyu,TONG Minming,JIANG Chunlu. Prediction of water inrush quantity from coal floor based on fuzzy evidence theory[J]. Journal of China Coal Society,2012,37(Sup.1):131–137.

    [6] 劉景,馮光俊,吳曉軍,等.基于層次分析法的富水區(qū)預測及其在頂板突水危險性評價中的作用[J]. 煤礦安全,2019,50(5):204–208.

    LIU Jing,F(xiàn)ENG Guangjun,WU Xiaojun,et al. Prediction of water rich area based on AHP and its role in risk assessment of roof water inrush[J]. Safety in Coal Mines,2019,50(5):204–208.

    [7] 尹會永,魏久傳,劉同彬,等. 基于多源信息復合的煤層底板突水評價[J]. 山東科技大學學報(自然科學版),2008,27(2):6–9.

    YIN Huiyong,WEI Jiuchuan,LIU Tongbin,et al. Evaluation of water inrush in seam floor based on multi-originated information complex[J]. Journal of Shandong University of Science and Technology(Natural Science),2008,27(2):6–9.

    [8] 韓承豪,魏久傳,謝道雷,等. 基于集對分析–可變模糊集耦合法的砂巖含水層富水性評價:以寧東礦區(qū)金家渠井田侏羅系直羅組含水層為例[J]. 煤炭學報,2020,45(7):2432–2443.

    HAN Chenghao,WEI Jiuchuan,XIE Daolei,et al. Water-richness evaluation of sandstone aquifer based on set pair analysis-variable fuzzy set coupling method:A case from Jurassic Zhiluo Formation of Jinjiaqu coal mine in Ningdong mining area[J]. Journal of China Coal Society,2020,45(7):2432–2443.

    [9] 李彥民,周晨陽,李鳳蓮. 基于代價敏感理論的多決策樹煤層底板突水預測模型[J]. 工礦自動化,2020,46(12):76–83.

    LI Yanmin,ZHOU Chenyang,LI Fenglian. Multi-decision tree prediction model for coal seam floor water inrush based on cost-sensitive theory[J]. Industry and Mine Automation,2020, 46(12):76–83.

    [10] 張曉亮. 熵權耦合層次分析賦權在煤層底板突水評價中的應用[J]. 煤田地質(zhì)與勘探,2017,45(3):91–95.

    ZHANG Xiaoliang. Application of entropy weight method and analytic hierarchy process in evaluation of water inrush from coal seam floor[J]. Coal Geology & Exploration,2017,45(3):91–95.

    [11] 劉偉韜,廖尚輝,劉士亮,等. 主成分logistic回歸分析在底板突水預測中的應用[J]. 遼寧工程技術大學學報(自然科學版),2015,34(8):905–909.

    LIU Weitao,LIAO Shanghui,LIU Shiliang,et al. Principal component logistic regression analysis in application of water outbursts from coal seam floor[J]. Journal of Liaoning Technical University(Natural Science),2015,34(8):905–909.

    [12] YIN Huiyong,ZHOU Wanfang,LAMOREAUX J W,et al. Water inrush conceptual site models for coal mines of China[J]. Environmental Earth Sciences,2018,77(22):746.

    [13] 趙琳琳,溫國鋒,邵良杉. 煤層底板突水危險性GSPCA-LSSVM評價模型[J]. 中國安全科學學報,2018,28(2):128–133.

    ZHAO Linlin,WEN Guofeng,SHAO Liangshan. GSPCA-LSSVM model for evaluating risk of coal floor groundwater bursting[J]. China Safety Science Journal,2018,28(2):128–133.

    [14] 劉晨雨,魏久傳,王杰,等. 基于AHP-TFN模型的底板突水危險性預測[J]. 中國礦業(yè),2019,28(8):124–129.

    LIU Chenyu,WEI Jiuchuan,WANG Jie,et al. Prediction of floor water inrush risk based on AHP-TFN model[J]. China Mining Magazine,2019,28(8):124–129.

    [15] 施龍青,曲興玥,韓進,等. 多模型融合評價煤層底板灰?guī)r巖溶突水危險性[J]. 煤炭學報,2019,44(8):2484–2493.

    SHI Longqing,QU Xingyue,HAN Jin,et al. Multi-model fusion for assessing the risk of inrush of limestone karst water through mine floor[J]. Journal of China Coal Society,2019,44(8):2484–2493.

    [16] 于小鴿,韓進,施龍青,等. 基于BP神經(jīng)網(wǎng)絡的底板破壞深度預測[J]. 煤炭學報,2009,34(6):731–736.

    YU Xiaoge,HAN Jin,SHI Longqing,et al. Forecast of destroyed floor depth based on BP neural networks[J]. Journal of China Coal Society,2009,34(6):731–736.

    [17] 陳建平,王春雷,王雪冬. 基于CNN神經(jīng)網(wǎng)絡的煤層底板突水預測[J]. 中國地質(zhì)災害與防治學報,2021,32(1):50–57.

    CHEN Jianping,WANG Chunlei,WANG Xuedong. Coal mine floor water inrush prediction based on CNN neural network[J]. The Chinese Journal of Geological Hazard and Control,2021,32(1):50–57.

    [18] 施龍青,張榮遨,徐東晶,等. 基于GWO-Elman神經(jīng)網(wǎng)絡的底板突水預測[J]. 煤炭學報,2020,45(7):2455–2463.

    SHI Longqing,ZHANG Rongao,XU Dongjing,et al. Prediction of water inrush from floor based on GWO-Elman neural network[J]. Journal of China Coal Society,2020,45(7):2455–2463.

    [19] 薛建凱. 一種新型的群智能優(yōu)化技術的研究與應用:麻雀搜索算法[D]. 上海:東華大學,2020.

    XUE Jiankai. Research and application of a novel swarm intelligence optimization technique:Sparrow search algorithm[D]. Shanghai:Donghua University,2020.

    [20] 祁春燕,邱國慶,張海榮. 底板突水預測模型的影響因素分析[J]. 武漢大學學報(信息科學版),2013,38(2):153–156.

    QI Chunyan,QIU Guoqing,ZHANG Hairong. Influencing factors analysis of floor water invasion prediction model[J]. Geomatics and Information Science of Wuhan University,2013,38(2):153–156.

    [21] YIN Huiyong,SHI Yongli,NIU Huigong,et al. A GIS-based model of potential groundwater yield zonation for a sandstone aquifer in the Juye coalfield,Shangdong,China[J]. Journal of Hydrology,2018,557:434–447.

    [22] 呂鑫,慕曉冬,張鈞,等. 混沌麻雀搜索優(yōu)化算法[J/OL]. 北京航空航天大學學報,2020:1–10. [2020-08-31]. https: //doi.org/10.137001j.bh.1001-5965.2020.0298

    LYU Xin,MU Xiaodong,ZHANG Jun,et al. Chaos sparrow search optimization algorithm[J/OL]. Journal of Beijing University of Aeronautics and Astronautics,2020:1–10. [2020-08-31]. https: //doi.org/10.137001j.bh.1001-5965.2020.0298

    [23] 王嶸冰,徐紅艷,李波,等. BP神經(jīng)網(wǎng)絡隱含層節(jié)點數(shù)確定方法研究[J]. 計算機技術與發(fā)展,2018,28(4):31–35.

    WANG Rongbing,XU Hongyan,LI Bo,et al. Research on method of determining hidden layer nodes in BP neural network[J]. Computer Technology and Development,2018,28(4):31–35.

    [24] 尹會永,趙涵,徐琳,等. 巖體質(zhì)量分級的改進模糊綜合評價法[J]. 金屬礦山,2020(7):53–58.

    YIN Huiyong,ZHAO Han,XU Lin,et al. Classification of rock mass in mine based on improved fuzzy comprehensive evaluation method[J]. Metal Mine,2020(7):53–58.

    [25] 李競贏,劉啟蒙,劉瑜,等. 基于GIS與熵值法的煤層頂板突水危險性評價[J]. 煤炭工程,2019,51(8):115–119.

    LI Jingying,LIU Qimeng,LIU Yu,et al. Risk assessment of water inrush from coal seam roof based on GIS and entropy method[J]. Coal Engineering,2019,51(8):115–119.

    Prediction model of water inrush from coal floor based on GA-BP neural network optimized by SSA and its application

    YIN Huiyong1, ZHOU Xinlong1, LANG Ning1, ZHANG Lifeng2, WANG Mingli3, WU Tao2, LI Xin4

    (1. College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao 266590, China; 2. Binhu Coal Mine, Zaozhuang Mining(Group) Co., LTD., Zaozhuang 277599, China; 3. Qingdao Eco-Environment Monitoring Center of Shandong Province, Qingdao 266003, China; 4. Ningxia Coal Science and Technology Research Institute Co., Ltd., Yinchuan 750021, China)

    With the increase of coal mining depth, coal production process is faced with complex water inrush mechanism and variable water inrush main control factors, and the uncertainties among the factors make the prediction of floor water inrush more difficult. In order to accurately predict the risk of floor water inrush, aiming at the small sample and non-linear problem of floor water inrush, firstly, genetic Algorithm is used to optimize the initial weights and thresholds of network random assignment, and then Sparrow Search Algorithm with strong search ability and good stability is selected to optimize the weights and thresholds for the second time, so as to establish the SSA-GA-BP neural network floor water inrush prediction model. Based on the analysis of geological and hydrological data of Binhu Coal Mine in Shandong Province, 8 factors including water pressure of aquifer, aquifer thickness, aquiclude thickness, fault density, fractal dimension value of fault, permeability coefficient, unit water inflow and floor failure depth are selected as the main control factors to predict floor water inrush, mapping the main controlling factors of 3D surface map projection. The Kriging interpolation method in surfer software is used to extract 50 data points as the input samples of the model(including 40 training sets and 10 test sets). The model is trained and studied. After the training error accuracy meets the requirements, the water inrush risk of 12 data points of 3 unmined working faces in Binhu Coal Mine is predicted. To verify the accuracy of the model, BP, GA-BP and SSA-GA-BP models are used to predict the test set; to avoid the one-sideness of comparing the model only with the prediction of BP network, the Fuzzy Comprehensive Evaluation Method, which determines the weight by Entropy Weight Method, is selected to predict the test set. The prediction results of each network model and method are compared with the actual values for analysis. The results show that the water inrush prediction error of GA-BP neural network model optimized by sparrow search algorithm is smaller, and the prediction accuracy is higher, which provides a scientific theoretical basis for mine water disaster prediction.

    prediction of water inrush from floor; Sparrow Search Algorithm; Genetic Algorithm; BP neural network; Entropy Weight Method; Fuzzy Comprehensive Evaluation

    語音講解

    TD745

    A

    1001-1986(2021)06-0175-11

    2021-04-20;

    2021-08-18

    國家重點研發(fā)計劃課題(2017YFC0804101);山東省自然科學基金項目(ZR2019MD013)

    尹會永,1979年生,男,河北故城人,博士,教授,碩士生導師,從事礦井水害防治方面的研究. E-mail:huiyongy@sdust.edu.cn

    尹會永,周鑫龍,郎寧,等. 基于SSA優(yōu)化的GA-BP神經(jīng)網(wǎng)絡煤層底板突水預測模型與應用[J]. 煤田地質(zhì)與勘探,2021,49(6):175–185. doi: 10.3969/j.issn.1001-1986.2021.06.021

    YIN Huiyong,ZHOU Xinlong,LANG Ning,et al. Prediction model of water inrush from coal floor based on GA-BP neural network optimized by SSA and its application[J]. Coal Geology & Exploration,2021,49(6):175–185. doi: 10.3969/j.issn.1001-1986.2021.06.021

    移動閱讀

    (責任編輯 周建軍)

    猜你喜歡
    突水危險性底板
    O-3-氯-2-丙烯基羥胺熱危險性及其淬滅研究
    礦井突水水源的判別方法
    危險性感
    礦井滯后突水機理探討
    輸氣站場危險性分析
    有色設備(2021年4期)2021-03-16 05:42:32
    基于AHP對電站鍋爐進行危險性分析
    李雅莊礦滯后突水主控因素分析
    軟土地基上超長底板節(jié)制閘沉降及應力分析
    大相嶺隧道高壓突水機理與預測分析
    底板巷一巷兩用的回采工作面防火實踐
    欧美黑人欧美精品刺激| 一卡2卡三卡四卡精品乱码亚洲| 露出奶头的视频| 久久久久久人人人人人| 丝袜在线中文字幕| 欧美日韩精品网址| 丝袜美足系列| 级片在线观看| 欧美国产日韩亚洲一区| 首页视频小说图片口味搜索| 操出白浆在线播放| 日韩av在线大香蕉| 嫩草影院精品99| 他把我摸到了高潮在线观看| 嫩草影视91久久| 国语自产精品视频在线第100页| 一级毛片女人18水好多| 成熟少妇高潮喷水视频| 国产av又大| 18禁美女被吸乳视频| 黄色视频不卡| 欧美黑人精品巨大| 女人高潮潮喷娇喘18禁视频| 亚洲人成网站在线播放欧美日韩| 精品卡一卡二卡四卡免费| 黄色丝袜av网址大全| av片东京热男人的天堂| 如日韩欧美国产精品一区二区三区| 亚洲av熟女| 一个人免费在线观看的高清视频| 天堂动漫精品| 日韩欧美一区视频在线观看| 韩国精品一区二区三区| 侵犯人妻中文字幕一二三四区| 欧美 亚洲 国产 日韩一| 看免费av毛片| 成人18禁在线播放| 十分钟在线观看高清视频www| 久9热在线精品视频| 亚洲色图 男人天堂 中文字幕| 中亚洲国语对白在线视频| 极品教师在线免费播放| 亚洲,欧美精品.| 免费少妇av软件| 在线视频色国产色| 精品电影一区二区在线| 激情视频va一区二区三区| 欧美国产精品va在线观看不卡| 黑人操中国人逼视频| 在线观看舔阴道视频| 91精品国产国语对白视频| 国产精品永久免费网站| 两人在一起打扑克的视频| 真人一进一出gif抽搐免费| 男女下面插进去视频免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 91大片在线观看| 国产精品亚洲av一区麻豆| 国产麻豆成人av免费视频| tocl精华| 国产精华一区二区三区| 久久久久久久精品吃奶| 亚洲人成77777在线视频| 国产三级黄色录像| 国产精品,欧美在线| 99re在线观看精品视频| 波多野结衣高清无吗| 亚洲精品中文字幕在线视频| 精品福利观看| 久久精品亚洲精品国产色婷小说| 欧美日本视频| 好看av亚洲va欧美ⅴa在| 亚洲中文日韩欧美视频| 久久 成人 亚洲| 黄色a级毛片大全视频| 亚洲欧洲精品一区二区精品久久久| 国产精品 欧美亚洲| 多毛熟女@视频| 午夜老司机福利片| 两性夫妻黄色片| 1024视频免费在线观看| 中文字幕人妻丝袜一区二区| 91成人精品电影| 国产精品影院久久| 欧美一区二区精品小视频在线| 亚洲欧美精品综合一区二区三区| 咕卡用的链子| 真人做人爱边吃奶动态| 午夜免费鲁丝| 最新美女视频免费是黄的| 国产亚洲欧美精品永久| 中出人妻视频一区二区| 精品无人区乱码1区二区| 亚洲色图 男人天堂 中文字幕| 免费看十八禁软件| 国产精品久久久久久人妻精品电影| 成在线人永久免费视频| 亚洲精品美女久久久久99蜜臀| 精品国产美女av久久久久小说| 午夜福利一区二区在线看| 日本 av在线| avwww免费| 午夜福利欧美成人| 成人三级黄色视频| 久久国产精品男人的天堂亚洲| 99re在线观看精品视频| 99re在线观看精品视频| avwww免费| 久久天堂一区二区三区四区| 午夜福利一区二区在线看| 日韩欧美一区视频在线观看| 婷婷六月久久综合丁香| 黄色女人牲交| 精品乱码久久久久久99久播| 亚洲男人天堂网一区| 国产一卡二卡三卡精品| 亚洲成人国产一区在线观看| 黄片小视频在线播放| 国产精品一区二区三区四区久久 | 免费无遮挡裸体视频| 变态另类成人亚洲欧美熟女 | 黄频高清免费视频| 久久九九热精品免费| 97超级碰碰碰精品色视频在线观看| 久久久久精品国产欧美久久久| 啦啦啦观看免费观看视频高清 | 中文字幕另类日韩欧美亚洲嫩草| 女警被强在线播放| 九色国产91popny在线| 9色porny在线观看| 91精品国产国语对白视频| 亚洲中文日韩欧美视频| 亚洲一区二区三区不卡视频| 搞女人的毛片| 国产精品一区二区精品视频观看| 在线十欧美十亚洲十日本专区| 在线观看免费视频日本深夜| 亚洲三区欧美一区| 国产精品久久久久久亚洲av鲁大| 最好的美女福利视频网| 午夜a级毛片| 嫩草影视91久久| 欧美日本亚洲视频在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 国产又爽黄色视频| 国产精品久久久人人做人人爽| 久久精品国产清高在天天线| 99在线人妻在线中文字幕| 一个人免费在线观看的高清视频| 日韩av在线大香蕉| 久久精品91无色码中文字幕| 给我免费播放毛片高清在线观看| 男人的好看免费观看在线视频 | 欧美性长视频在线观看| 免费在线观看视频国产中文字幕亚洲| 一级a爱视频在线免费观看| 亚洲成人免费电影在线观看| 国产高清视频在线播放一区| av中文乱码字幕在线| 18禁国产床啪视频网站| 免费高清视频大片| 黄色视频,在线免费观看| 国产亚洲精品久久久久5区| 亚洲国产欧美一区二区综合| 电影成人av| 亚洲一区高清亚洲精品| 中文字幕色久视频| 亚洲精品国产一区二区精华液| 日日夜夜操网爽| 亚洲一区中文字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 夜夜看夜夜爽夜夜摸| 日韩免费av在线播放| 亚洲精品在线观看二区| 一区福利在线观看| 国产成年人精品一区二区| 淫秽高清视频在线观看| 高清黄色对白视频在线免费看| 久久精品成人免费网站| www.www免费av| 不卡一级毛片| 午夜两性在线视频| 999久久久国产精品视频| 国产三级在线视频| 亚洲第一青青草原| 露出奶头的视频| 久久久久国内视频| 人人妻人人澡人人看| 多毛熟女@视频| 黑人操中国人逼视频| 久久久久国产精品人妻aⅴ院| 欧美亚洲日本最大视频资源| 亚洲av成人不卡在线观看播放网| 一区福利在线观看| АⅤ资源中文在线天堂| 亚洲av日韩精品久久久久久密| 免费看美女性在线毛片视频| 一级,二级,三级黄色视频| 欧美色欧美亚洲另类二区 | 欧美亚洲日本最大视频资源| 欧美一区二区精品小视频在线| 亚洲成人精品中文字幕电影| 亚洲av美国av| 日本a在线网址| 啪啪无遮挡十八禁网站| 多毛熟女@视频| 国产午夜福利久久久久久| 99精品欧美一区二区三区四区| 久久久精品欧美日韩精品| 精品久久蜜臀av无| 曰老女人黄片| 久热这里只有精品99| 香蕉丝袜av| 国产极品粉嫩免费观看在线| 午夜免费鲁丝| 亚洲欧美精品综合一区二区三区| 日韩 欧美 亚洲 中文字幕| 国产精品自产拍在线观看55亚洲| 两性夫妻黄色片| 久久久久精品国产欧美久久久| 黄色毛片三级朝国网站| 亚洲欧美精品综合一区二区三区| 午夜视频精品福利| 久久狼人影院| 久久中文字幕人妻熟女| 午夜福利视频1000在线观看 | 99久久久亚洲精品蜜臀av| 国产主播在线观看一区二区| 免费高清在线观看日韩| 久久亚洲真实| 国产成人欧美在线观看| 在线国产一区二区在线| 九色亚洲精品在线播放| 日韩欧美国产一区二区入口| 久久婷婷人人爽人人干人人爱 | 50天的宝宝边吃奶边哭怎么回事| 欧美日本视频| 美女国产高潮福利片在线看| 少妇粗大呻吟视频| 婷婷六月久久综合丁香| 国产黄a三级三级三级人| 欧美成人一区二区免费高清观看 | 久久久国产成人精品二区| 国产日韩一区二区三区精品不卡| 美国免费a级毛片| 18禁观看日本| 久久狼人影院| 久久精品人人爽人人爽视色| 麻豆av在线久日| 国产成人啪精品午夜网站| 91在线观看av| 亚洲欧美激情在线| 精品久久久久久,| 亚洲精品粉嫩美女一区| 18美女黄网站色大片免费观看| 色av中文字幕| 好看av亚洲va欧美ⅴa在| 亚洲中文字幕一区二区三区有码在线看 | 夜夜夜夜夜久久久久| 午夜免费鲁丝| 精品福利观看| 国产精品香港三级国产av潘金莲| 亚洲精品国产精品久久久不卡| 亚洲国产欧美一区二区综合| 色综合欧美亚洲国产小说| 一级毛片精品| 亚洲专区国产一区二区| 精品人妻在线不人妻| 禁无遮挡网站| 中文字幕av电影在线播放| 国产精品免费视频内射| 日本 av在线| 久久亚洲真实| 成熟少妇高潮喷水视频| 777久久人妻少妇嫩草av网站| 人人妻人人爽人人添夜夜欢视频| 一边摸一边抽搐一进一小说| 国产成人欧美在线观看| 我的亚洲天堂| 亚洲欧美日韩另类电影网站| 久久人人精品亚洲av| 久久中文字幕一级| 国产成人一区二区三区免费视频网站| 日韩大尺度精品在线看网址 | 18禁美女被吸乳视频| svipshipincom国产片| 一级黄色大片毛片| videosex国产| 婷婷丁香在线五月| 欧美亚洲日本最大视频资源| 国产精品一区二区三区四区久久 | 日韩成人在线观看一区二区三区| 91麻豆精品激情在线观看国产| 色精品久久人妻99蜜桃| 国产亚洲欧美精品永久| 国产成人精品久久二区二区免费| 日韩免费av在线播放| 欧美日本亚洲视频在线播放| 久久久精品欧美日韩精品| 日韩有码中文字幕| 制服人妻中文乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 97人妻精品一区二区三区麻豆 | 日韩欧美一区视频在线观看| 电影成人av| 级片在线观看| 亚洲中文字幕日韩| 国产一卡二卡三卡精品| 亚洲人成电影观看| 最近最新中文字幕大全电影3 | 国产亚洲精品综合一区在线观看 | aaaaa片日本免费| 中文字幕av电影在线播放| 国产一级毛片七仙女欲春2 | 国产精品二区激情视频| 中亚洲国语对白在线视频| 国产亚洲欧美精品永久| 一边摸一边抽搐一进一出视频| 久久国产乱子伦精品免费另类| 亚洲av日韩精品久久久久久密| 一区二区三区精品91| 亚洲欧洲精品一区二区精品久久久| 日韩精品中文字幕看吧| 国产又色又爽无遮挡免费看| 日韩有码中文字幕| 国产成+人综合+亚洲专区| 国产一区在线观看成人免费| 黄片大片在线免费观看| 亚洲一区二区三区不卡视频| 人人妻人人爽人人添夜夜欢视频| 无人区码免费观看不卡| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av日韩精品久久久久久密| 午夜成年电影在线免费观看| 久久久久亚洲av毛片大全| 亚洲五月天丁香| 一级毛片高清免费大全| 黄色女人牲交| 国产亚洲欧美98| 久热爱精品视频在线9| 日韩欧美免费精品| 制服人妻中文乱码| 久久国产精品男人的天堂亚洲| 国产欧美日韩综合在线一区二区| 99国产精品一区二区三区| 久久精品国产亚洲av香蕉五月| 老鸭窝网址在线观看| 精品久久久久久成人av| 777久久人妻少妇嫩草av网站| 久久精品国产清高在天天线| 久久精品成人免费网站| 正在播放国产对白刺激| 亚洲专区中文字幕在线| 一级片免费观看大全| 最新美女视频免费是黄的| 亚洲成人久久性| a级毛片在线看网站| 国产亚洲精品综合一区在线观看 | 久久久久久久久免费视频了| 免费看a级黄色片| 桃红色精品国产亚洲av| 国产成+人综合+亚洲专区| 啦啦啦 在线观看视频| 91大片在线观看| 一级a爱片免费观看的视频| 国产亚洲av嫩草精品影院| av电影中文网址| 中文字幕人妻丝袜一区二区| 婷婷丁香在线五月| 久久人人97超碰香蕉20202| 亚洲人成电影观看| 不卡av一区二区三区| 亚洲男人天堂网一区| 国产97色在线日韩免费| 成人av一区二区三区在线看| 搡老妇女老女人老熟妇| 亚洲男人的天堂狠狠| 曰老女人黄片| 久久天堂一区二区三区四区| 一进一出抽搐动态| 超碰成人久久| 国产一卡二卡三卡精品| 国产成人欧美在线观看| 久久久久久人人人人人| 亚洲精品久久国产高清桃花| 欧美绝顶高潮抽搐喷水| 国产欧美日韩一区二区三| 久久久精品国产亚洲av高清涩受| 亚洲欧美激情在线| 国产精品爽爽va在线观看网站 | 欧美激情 高清一区二区三区| e午夜精品久久久久久久| 国产成人欧美| 久久亚洲真实| 欧美一级a爱片免费观看看 | e午夜精品久久久久久久| 国产片内射在线| 亚洲欧美精品综合一区二区三区| 欧美av亚洲av综合av国产av| 国产黄a三级三级三级人| 日本vs欧美在线观看视频| 亚洲av五月六月丁香网| 一个人免费在线观看的高清视频| www.精华液| 亚洲专区中文字幕在线| 一级片免费观看大全| 久久 成人 亚洲| 99精品欧美一区二区三区四区| 欧美另类亚洲清纯唯美| 黄色 视频免费看| 免费看十八禁软件| 757午夜福利合集在线观看| 俄罗斯特黄特色一大片| 丁香六月欧美| 欧美黑人欧美精品刺激| 亚洲欧洲精品一区二区精品久久久| 美女免费视频网站| 日本vs欧美在线观看视频| 麻豆久久精品国产亚洲av| 久久天躁狠狠躁夜夜2o2o| 香蕉久久夜色| 亚洲熟妇熟女久久| 国内久久婷婷六月综合欲色啪| 人人妻人人澡欧美一区二区 | 在线观看66精品国产| 叶爱在线成人免费视频播放| 日日夜夜操网爽| 搡老岳熟女国产| e午夜精品久久久久久久| 波多野结衣高清无吗| 男人操女人黄网站| 又紧又爽又黄一区二区| 69av精品久久久久久| 国产免费男女视频| 黄色视频,在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 一进一出抽搐gif免费好疼| 一区二区三区精品91| av免费在线观看网站| 啦啦啦 在线观看视频| 天堂影院成人在线观看| 涩涩av久久男人的天堂| 一区在线观看完整版| 99国产精品一区二区三区| 国产私拍福利视频在线观看| 99香蕉大伊视频| 不卡av一区二区三区| 超碰成人久久| 亚洲欧美精品综合一区二区三区| 叶爱在线成人免费视频播放| 精品国产乱子伦一区二区三区| 午夜福利一区二区在线看| 亚洲国产中文字幕在线视频| 村上凉子中文字幕在线| 国产麻豆69| 亚洲中文字幕日韩| 少妇裸体淫交视频免费看高清 | 亚洲午夜理论影院| av欧美777| 亚洲人成伊人成综合网2020| 色婷婷久久久亚洲欧美| 国产成年人精品一区二区| 色av中文字幕| 国产在线观看jvid| 精品无人区乱码1区二区| 国产精品综合久久久久久久免费 | 国产亚洲欧美精品永久| 午夜视频精品福利| 免费人成视频x8x8入口观看| 欧美最黄视频在线播放免费| 亚洲自拍偷在线| 最近最新免费中文字幕在线| 久久热在线av| 色在线成人网| 男人的好看免费观看在线视频 | 日韩欧美一区视频在线观看| 国产av又大| 又黄又粗又硬又大视频| 十八禁网站免费在线| 国产成人欧美在线观看| 成人特级黄色片久久久久久久| 淫秽高清视频在线观看| svipshipincom国产片| www.自偷自拍.com| 日韩av在线大香蕉| 久久国产亚洲av麻豆专区| 国产不卡一卡二| 午夜精品久久久久久毛片777| 亚洲国产精品合色在线| 欧美日本中文国产一区发布| 色综合欧美亚洲国产小说| 每晚都被弄得嗷嗷叫到高潮| 午夜亚洲福利在线播放| 久久久久久久久免费视频了| 亚洲国产毛片av蜜桃av| 男人操女人黄网站| 又黄又爽又免费观看的视频| 19禁男女啪啪无遮挡网站| 亚洲欧美精品综合久久99| 国产亚洲精品久久久久久毛片| 午夜福利免费观看在线| 国产精品 欧美亚洲| 男人舔女人的私密视频| 一个人观看的视频www高清免费观看 | 88av欧美| 欧美午夜高清在线| 两人在一起打扑克的视频| 乱人伦中国视频| 国产麻豆69| 纯流量卡能插随身wifi吗| 亚洲人成电影观看| 在线十欧美十亚洲十日本专区| 亚洲精品粉嫩美女一区| 亚洲欧洲精品一区二区精品久久久| 精品久久久久久,| 久99久视频精品免费| 亚洲人成网站在线播放欧美日韩| 黄片大片在线免费观看| 国产一区二区三区视频了| 在线观看免费午夜福利视频| 99国产极品粉嫩在线观看| 国产亚洲欧美98| 亚洲成人精品中文字幕电影| 99久久综合精品五月天人人| 婷婷六月久久综合丁香| 黑人欧美特级aaaaaa片| 国产一区二区在线av高清观看| 9191精品国产免费久久| 51午夜福利影视在线观看| 成年版毛片免费区| 免费看十八禁软件| 不卡一级毛片| 免费在线观看黄色视频的| 婷婷精品国产亚洲av在线| 女警被强在线播放| 亚洲伊人色综图| 久久久国产成人精品二区| 日韩欧美国产一区二区入口| 久久这里只有精品19| 亚洲一区二区三区色噜噜| 69精品国产乱码久久久| 精品一区二区三区视频在线观看免费| 久久久久九九精品影院| 久久国产精品影院| 色播亚洲综合网| 久久精品亚洲熟妇少妇任你| 成在线人永久免费视频| 亚洲精品av麻豆狂野| 亚洲中文av在线| 在线国产一区二区在线| 嫩草影院精品99| 国产亚洲精品一区二区www| 日韩欧美免费精品| 手机成人av网站| 亚洲精品久久国产高清桃花| 麻豆成人av在线观看| 香蕉国产在线看| 12—13女人毛片做爰片一| 两性夫妻黄色片| 国产亚洲精品久久久久久毛片| 夜夜爽天天搞| 久久国产精品影院| 如日韩欧美国产精品一区二区三区| 国产一级毛片七仙女欲春2 | 人成视频在线观看免费观看| 国产精品1区2区在线观看.| av超薄肉色丝袜交足视频| 又紧又爽又黄一区二区| 日韩大码丰满熟妇| 在线观看www视频免费| 久久天躁狠狠躁夜夜2o2o| 女性被躁到高潮视频| 久久中文字幕人妻熟女| 亚洲精品av麻豆狂野| 亚洲中文av在线| 亚洲午夜理论影院| 国产精品av久久久久免费| 精品国产国语对白av| 久久中文字幕人妻熟女| 午夜福利影视在线免费观看| 久久久精品欧美日韩精品| 成人精品一区二区免费| 亚洲精品美女久久av网站| 黄片播放在线免费| 91麻豆精品激情在线观看国产| 国产av一区在线观看免费| 日本 av在线| 日韩欧美在线二视频| 男女下面进入的视频免费午夜 | 男女床上黄色一级片免费看| 国产在线精品亚洲第一网站| 他把我摸到了高潮在线观看| 国产片内射在线| 国产蜜桃级精品一区二区三区| 午夜成年电影在线免费观看| 免费少妇av软件| 精品午夜福利视频在线观看一区| 国产成人系列免费观看| 亚洲一区高清亚洲精品| 国产又爽黄色视频| 国产av一区在线观看免费| 中文字幕人妻丝袜一区二区| 夜夜躁狠狠躁天天躁| 长腿黑丝高跟| 一区二区三区国产精品乱码| 大陆偷拍与自拍| 国产免费男女视频| 久久亚洲真实| 国产精品亚洲一级av第二区| 国产成人影院久久av| 午夜视频精品福利| 久久午夜亚洲精品久久| 亚洲一区二区三区色噜噜| 真人一进一出gif抽搐免费| 丝袜在线中文字幕|