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

    基于GIS與BP神經(jīng)網(wǎng)絡(luò)的礦區(qū)塌陷易發(fā)性預(yù)測

    2021-10-17 14:00:18王紳皓謝婉麗奚家米
    煤礦安全 2021年9期
    關(guān)鍵詞:易發(fā)礦區(qū)神經(jīng)網(wǎng)絡(luò)

    王紳皓,謝婉麗,奚家米,井 旭

    (1.西北大學(xué) 地質(zhì)學(xué)系,陜西 西安 710069;2.大陸動力學(xué)國家重點實驗室,陜西 西安 710069)

    地面塌陷是礦山常見的地質(zhì)災(zāi)害之一,礦區(qū)的地面塌陷為采礦工作及當(dāng)?shù)鼐用駧砹藝?yán)重的安全隱患及經(jīng)濟損失[1]。因此,礦區(qū)塌陷的預(yù)測是保障生產(chǎn)工作正常運行,保護居民生命財產(chǎn)安全的重要工作[2]。針對塌陷易發(fā)性的預(yù)測,已經(jīng)有眾多學(xué)者進行了大量的研究,并提出了各種判別方法。劉寶琛等[3]提出將概率積分法運用到地下開采等實際工程引發(fā)地面塌陷的預(yù)測工作中;謝和平等[4]使用相似材料作了巷道和采場的巖層塌陷實驗,并基于非線性的連續(xù)介質(zhì)模型對地下工程的塌陷進行了非線性力學(xué)有限元分析??讘椓⒌萚5]論述了采空區(qū)地表為斜坡條件下,上覆巖層發(fā)生多種變形、破壞的力學(xué)特征;陳學(xué)軍[6]、包惠明[7]等分別將GIS技術(shù)、模糊數(shù)學(xué)理論應(yīng)用在巖溶塌陷評價模型中,對其易發(fā)性進行定量評價。以上方法在塌陷易發(fā)性預(yù)測中均有較好的效果,但由于巖土體的空間特性特殊,建立力學(xué)模型的方法復(fù)雜且難以簡化;數(shù)值模擬難以全面考慮礦山塌陷中的復(fù)雜因素,選取合理參數(shù),模糊評判存在超模糊現(xiàn)象。為補充以往研究不足,為此選用BP神經(jīng)網(wǎng)絡(luò)進行礦區(qū)塌陷的預(yù)測研究。BP神經(jīng)網(wǎng)絡(luò)具有主動學(xué)習(xí)、多層次,可以通過學(xué)習(xí)大量的樣本來獲得適合的連接權(quán)值,反映出各種復(fù)雜的非線性關(guān)系[8]。綜合GIS系統(tǒng)空間分析功能在礦區(qū)塌陷預(yù)測突出的優(yōu)越性,采用柵格為評價單元,提高預(yù)測精度,最大限度地體現(xiàn)出礦區(qū)內(nèi)部塌陷易發(fā)性的空間差異[9],為模型預(yù)測提供足夠的樣本數(shù)。

    1 礦區(qū)塌陷預(yù)測評價模型

    BP神經(jīng)網(wǎng)絡(luò)算法是一種信息處理體系,利用神經(jīng)網(wǎng)絡(luò)能夠進行復(fù)雜的邏輯運算,實現(xiàn)多因子相互作用的非線性關(guān)系。BP神經(jīng)網(wǎng)絡(luò)由輸入層、隱含層和輸出層3部分組成[10]。網(wǎng)絡(luò)訓(xùn)練過程主要信號的正向傳播和誤差反向傳播[11],通過BP算法的迭代依次對隱含層到輸出層、輸入層到隱含層的權(quán)重和偏置進行調(diào)節(jié),使權(quán)重系數(shù)和閾值參數(shù)逐漸收斂于合適的值,當(dāng)相鄰2次誤差值達到目標(biāo)精度時,模型趨于合理,訓(xùn)練完成。

    1.1 信號的正向傳播

    設(shè)xi為人工神經(jīng)網(wǎng)絡(luò)隱含層的輸入值,則隱含層輸出值Hk為:

    式中:n為輸入節(jié)點數(shù);wik為隱含層間連接權(quán)值;ak為隱含層閾值;g為隱含層激勵函數(shù)。

    激活函數(shù)的選取一般認(rèn)為理想的激活函數(shù)是階躍函數(shù),但階躍函數(shù)具有不連續(xù)、不可導(dǎo)的缺點。塌陷預(yù)測采用Sigmoid函數(shù):

    Sigmoid函數(shù)的優(yōu)勢在于其不同于階躍函數(shù)的二值化輸出結(jié)果,其輸出值為(0,1)區(qū)間內(nèi)的可控數(shù)值,該數(shù)值可以用于表示塌陷易發(fā)性的大小。

    1.2 誤差反向傳播

    通過連接權(quán)值和偏置的更新使誤差函數(shù)取值逐漸減小,當(dāng)誤差取最小值時,認(rèn)為這時的權(quán)重參數(shù)趨于最優(yōu)解[12]。BP神經(jīng)網(wǎng)絡(luò)算法基于梯度下降法來求解修正權(quán)值。進行參數(shù)的傳遞分析可以發(fā)現(xiàn)存在這樣的影響傳遞鏈條關(guān)系,即wkj影響輸出層輸出值,最終影響誤差的大小。反向傳播中權(quán)值更新公式:

    式中:wkj為連接權(quán)值;ej為網(wǎng)絡(luò)預(yù)測誤差;p為輸入節(jié)點數(shù);η為學(xué)習(xí)速率。

    當(dāng)相鄰2次之間的誤差值小于目標(biāo)值時認(rèn)為算法收斂,算法迭代結(jié)束[14]。

    2 礦區(qū)概況及塌陷影響因子分析

    2.1 礦區(qū)地質(zhì)概況

    研究礦區(qū)地處陜西省廟哈孤礦區(qū)東南部,該礦區(qū)為典型的黃土梁峁地貌,溝壑縱橫,梁面寬緩平坦,場地內(nèi)地勢較為平坦。區(qū)內(nèi)地層主要為三迭系細(xì)粒石英砂巖,厚度約80~200 m;侏羅系粉砂巖,泥巖夾煤層,厚度約89.78 m,新近系紅色黏土,平均厚度3.1 m;第四系黃土覆蓋大部分山梁、緩坡,平均厚度24.15 m。

    礦區(qū)降水較少,年內(nèi)降水量變化大,水系較發(fā)育,流經(jīng)礦區(qū)的水系主要為2條:孤山川河自井田西側(cè)邊界穿過,府谷縣內(nèi)流長57 km,流域面積1 018 km2,年平均流量3.48 m3/s,溝谷寬緩;流小南川溝流域面積約6 km2,河流受區(qū)域獨特的地質(zhì)地貌條件及氣象因素的控制,雨季洪峰多。

    井田內(nèi)設(shè)計開采煤層為22、31、42、51、52上、52煤,煤層平均可采總厚度12.22 m。以52煤為例,煤層厚度1.47~3.32 m,平均2.23 m,煤厚變化較小,埋深淺,煤層埋藏深度在0~225.66 m,煤層頂板巖性為1套厚層粉砂巖、局部夾薄層粉砂質(zhì)泥巖,厚16 m。

    2.2 塌陷發(fā)育情況

    地質(zhì)災(zāi)害類型主要為滑坡、崩塌、塌陷和泥石流,其中因煤礦開采,該地區(qū)塌陷顯著增多。截至2019年10月,研究區(qū)塌陷已增至10余處,目前塌陷有形成貫通的趨勢。其中塌陷面積最大達到3 km2,嚴(yán)重影響到礦區(qū)煤礦的繼續(xù)開采以及當(dāng)?shù)鼐用袢粘I?,礦區(qū)塌陷分布如圖1。

    圖1 礦區(qū)塌陷分布圖Fig.1 Distribution map of mining area collapse

    以二采區(qū)形成的地面塌陷為例,其形成時間為2013—2018年,位于井田中部,走向長度1.91 km,傾斜長度1.06 km,該采區(qū)52號煤層平均厚度2.38 m,煤層結(jié)構(gòu)較簡單。該地面塌陷區(qū)面積較大,塌陷區(qū)東側(cè)較寬,向西逐漸變窄,地面裂縫發(fā)育,裂縫長約7~50 m。裂縫成組平行發(fā)育,間距5~20 m不等,部分裂縫可見深度達1.5 m,兩側(cè)錯落10~50 cm,在斜坡地帶局部發(fā)生了滑塌,地面塌陷裂縫方向與采取工作面方向相一致。

    2.3 塌陷影響因子選取

    2.3.1 數(shù)據(jù)源

    研究區(qū)塌陷預(yù)測評價建模的數(shù)據(jù)源主要有:①塌陷編錄信息和野外調(diào)查獲取的相關(guān)資料,包括研究區(qū)1∶2 000井田勘查報告、礦井建設(shè)工程地質(zhì)災(zāi)害易發(fā)性評估報告、礦井設(shè)計書等;②DEM高程數(shù)據(jù),用于提取地形地貌等基礎(chǔ)環(huán)境因子;③研究區(qū)地質(zhì)調(diào)查資料。

    2.3.2 影響因子

    礦區(qū)地面塌陷大多數(shù)是由于頂板的破壞而引起的,即礦體被開挖、采掘后,失去下部巖層支持的頂板在自重和次生應(yīng)力作用下發(fā)生剪切或拉裂破壞[13],該礦區(qū)采用走向長壁后退式綜合機械化采煤方法,全部垮落法管理頂板。研究區(qū)垂直剖面由上至下可簡化為松散覆蓋物、堅硬巖層、含水層等,區(qū)內(nèi)分布構(gòu)造裂隙,采空區(qū)垂直剖面如圖2。

    圖2 采空區(qū)垂直剖面圖Fig.2 Vertical section view of goaf

    在煤層回采后,上覆巖體完全垮落,垮落帶可達松散覆蓋物底部,在松散覆蓋層形成裂隙,使松散層本身發(fā)生垮落,導(dǎo)致地表塌陷[14]。

    影響礦區(qū)塌陷的因素多種多樣,各個采區(qū)情況具有獨特性[15],眾多研究表明塌陷的發(fā)生是基礎(chǔ)環(huán)境和工程影響共同作用的結(jié)果。采煤方法和頂板管理方法在該礦區(qū)不同工作面中相同,在預(yù)測模型中不予考慮。塌陷預(yù)測模型易發(fā)性因素包括:①覆巖強度是反映巖體工程地質(zhì)特征的最根本因素,直接影響著巖體的物理力學(xué)性質(zhì)及其受力變形的全過程;②水文地質(zhì)特征,地表水、地下水對巖土體的作用是一種地質(zhì)應(yīng)力作用,水與巖土體相互作用改變著巖土體的物理、化學(xué)及力學(xué)性質(zhì);③地質(zhì)構(gòu)造主要包括不連續(xù)結(jié)構(gòu)面、構(gòu)造應(yīng)力,地質(zhì)構(gòu)造活動使得巖體中聚集大量的彈性變形能,產(chǎn)生構(gòu)造應(yīng)力,受到開挖擾動的影響,構(gòu)造應(yīng)力釋放使巖體中出現(xiàn)局部應(yīng)力集中現(xiàn)象,誘發(fā)地面塌陷[16];④巖層的變形和破壞是從直接頂開始,自下而上擴散直至地面的漸變過程,終采時間越長,塌陷的易發(fā)性越高;⑤采礦工程將巖體中原始應(yīng)力平衡狀態(tài)破壞,并使其周圍一定范圍的巖層應(yīng)力發(fā)生重新分布,而開采深度和跨度是影響著擾動后原巖應(yīng)力的大小、方向與分布狀態(tài)的重要因素;⑥采煤高度影響著塌陷的易發(fā)程度,當(dāng)采煤高度較高或多煤層迭置開采時,覆巖經(jīng)多次塌陷分級傳至地面造成地表沉陷,產(chǎn)生大量裂縫。

    結(jié)合研究區(qū)的實際情況,本研究選取研究區(qū)3個基礎(chǔ)環(huán)境因子和5個工程影響因子,其中基礎(chǔ)環(huán)境因子包括水文特征、地層構(gòu)造和覆巖強度[17],工程影響因子包括終采時間、頂板跨度、開采深度、采煤高度以及空間迭置層數(shù)。礦區(qū)塌陷影響因子量化方法及易發(fā)性分級見表1。

    表1 礦區(qū)塌陷影響因子量化方法及易發(fā)性分級Table 1 Quantification method and stability classification of mining area collapse influence factor

    2.4 塌陷易發(fā)性單元劃分

    利用預(yù)測模型開展塌陷區(qū)易發(fā)性預(yù)測的1個關(guān)鍵問題是如何進行預(yù)測單元的劃分[18]。常用的預(yù)測單元多為柵格[19]、同期開采區(qū)域、行政區(qū)劃單元等。其中,柵格具有表達目標(biāo)明確,模型計算效率高,數(shù)據(jù)后處理易操作等優(yōu)點[20]。GIS系統(tǒng)具有范圍廣闊且功能強大的空間分析能力,使用其核心工具集可以將數(shù)據(jù)快速處理為柵格化,為數(shù)據(jù)空間分析提供了豐富、便利的環(huán)境和算子[21]。

    研究利用GIS的數(shù)據(jù)管理將離散型、連續(xù)型數(shù)據(jù)整合為面數(shù)據(jù),利用GIS的空間分析功能對各影響因子進行柵格化,根據(jù)礦區(qū)現(xiàn)有資料及GIS系統(tǒng)有效精度選擇適宜的網(wǎng)格大小為20 m×20 m,將48.9 km2研究區(qū)按大小劃分為119 762個單元,通過最鄰近分配法分別對各影響因子單元進行指標(biāo)采樣,采樣后的柵格數(shù)據(jù)中每1個像元點為包含8個塌陷影響因子的樣本組。礦區(qū)的8個影響因子分區(qū)如圖3。

    圖3 影響因子分區(qū)圖Fig.3 Impact factor classification

    3 礦區(qū)塌陷預(yù)測

    采用GIS系統(tǒng)對礦區(qū)影響因子提取466×257個單元,對某些數(shù)據(jù)偏離平均值過大的畸點取相鄰柵格單元平均值,同時為防止由于數(shù)據(jù)數(shù)量級的不統(tǒng)一,對各類評價指標(biāo)數(shù)據(jù)進行歸一化處理[20]。通過Matlab軟件利用礦區(qū)中部及南部466×173個單元作為樣本進行訓(xùn)練構(gòu)建BP神經(jīng)網(wǎng)絡(luò)模型,并保留礦區(qū)北部466×84個單元檢驗預(yù)測結(jié)果。

    3.1 神經(jīng)網(wǎng)絡(luò)模型塌陷預(yù)測

    神經(jīng)網(wǎng)絡(luò)隱含層的確定是建模的1個關(guān)鍵問題,隱含層節(jié)點過少,網(wǎng)絡(luò)難以識別樣本,從而無法建立復(fù)雜的映射關(guān)系,導(dǎo)致模型訓(xùn)練能力降低,預(yù)測效果差;隱含層節(jié)點過多則會導(dǎo)致網(wǎng)絡(luò)的過擬合以及訓(xùn)練效率的降低。目前常用的方法是通過多次實驗以確定隱含層的最優(yōu)節(jié)點數(shù)[22]。參考式(5)~式(7)進行隱含層節(jié)點選?。?/p>

    式中:l為隱含層節(jié)點數(shù);n為輸入層節(jié)點數(shù);m為輸出層節(jié)點數(shù);a為0~10之間的常數(shù)。

    最終確定隱含層節(jié)點數(shù)為5。

    利用上文生成的輸入、輸出訓(xùn)練樣本在Matlab環(huán)境下進行神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練和檢驗工作。訓(xùn)練使用newff神經(jīng)網(wǎng)絡(luò)函數(shù)。為保證預(yù)測誤差較小,隱含層節(jié)點傳遞函數(shù)選用S型函數(shù)tansig函數(shù),輸出層節(jié)點傳遞選用purelin函數(shù),控制學(xué)習(xí)速率為0.1,選擇Levenberg_Marquardt的BP算法訓(xùn)練函數(shù)trainlm,性能分析選擇均方差性能分析函數(shù)mse,誤差允許限度設(shè)置為1×10-2。訓(xùn)練過程發(fā)現(xiàn)迭代次數(shù)為2 253時誤差達到預(yù)期目標(biāo),可以認(rèn)為此時模型訓(xùn)練完成。

    將研究區(qū)柵格化數(shù)據(jù)組代入訓(xùn)練完成的神經(jīng)網(wǎng)絡(luò)模型,將模型輸出結(jié)果使用image函數(shù)進行圖像化,礦區(qū)塌陷預(yù)測分區(qū)圖如圖4。

    圖4 礦區(qū)塌陷預(yù)測分區(qū)圖Fig.4 Mine collapse prediction classification

    3.2 塌陷預(yù)測結(jié)果檢驗

    從BP神經(jīng)網(wǎng)絡(luò)模型的預(yù)測分區(qū)圖來看,預(yù)測區(qū)塌陷情況整體符合礦區(qū)采空塌陷的實際情況:礦區(qū)塌陷主要分布于下部有采空的地區(qū),同時塌陷受到地下巷道支護控制,具備煤柱支護的地區(qū)塌陷易發(fā)性最低,且支護具備一定區(qū)域范圍內(nèi)的塌陷控制能力,將塌陷分割為多個區(qū)域且塌陷易發(fā)性明顯降低。地層巖性變化較小時,上覆巖層硬質(zhì)巖層占比越大,采空區(qū)塌陷的易發(fā)性越低。采空區(qū)頂板越厚、跨度越小,越不容易發(fā)生離層破壞或頂板折斷,對應(yīng)區(qū)域塌陷易發(fā)性越低。塌陷的易發(fā)性在一定程度上受到水文特征的影響,其中,中南部礦區(qū)內(nèi)有水系發(fā)育的地區(qū)塌陷易發(fā)性較高。

    為檢驗BP神經(jīng)網(wǎng)絡(luò)的預(yù)測結(jié)果,根據(jù)礦區(qū)實際發(fā)生的地面塌陷作為檢驗標(biāo)準(zhǔn),采用ROC曲線完成神經(jīng)網(wǎng)絡(luò)預(yù)測模型的結(jié)果驗證[23]。大量研究表明,ROC曲線可以作為驗證預(yù)測模型的良好指標(biāo)[24]。選取礦區(qū)北部466×84個單元采用繪制ROC曲線,計算得到曲線下面積AUC=0.902,預(yù)測精度可達90.2%,BP模型ROC曲線測試圖如圖5。

    圖5 BP模型ROC曲線測試圖Fig.5 ROC curve test diagram of BP model

    驗證結(jié)果表明:采用BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型對礦區(qū)地下開采引起的地面塌陷可以很好地達到預(yù)期的預(yù)測結(jié)果。

    4 結(jié)語

    1)采礦區(qū)是塌陷地質(zhì)災(zāi)害高發(fā)區(qū)之一,通過選取3個基礎(chǔ)環(huán)境因子和5個工程影響因子構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型,BP神經(jīng)網(wǎng)絡(luò)具有很強的主動學(xué)習(xí)、適應(yīng)能力,反映出礦區(qū)地面塌陷發(fā)育過程中各種復(fù)雜的非線性關(guān)系。

    2)結(jié)合GIS技術(shù)將礦區(qū)數(shù)據(jù)柵格化和人工神經(jīng)網(wǎng)絡(luò)模型應(yīng)用于塌陷易發(fā)性分區(qū),可以很好地滿足預(yù)測精度,同時在訓(xùn)練計算過程中的迭代次數(shù)和耗時較少,可以達到快速、優(yōu)良的收斂效果。

    3)將預(yù)測模型與實際調(diào)查資料進行對比,BP神經(jīng)網(wǎng)絡(luò)模型用于礦區(qū)易發(fā)性預(yù)測的效果較好,適用于礦區(qū)開采結(jié)束后地面塌陷的預(yù)測,為環(huán)境保護、土地復(fù)墾方案以及礦區(qū)工程建設(shè)提供可行性、科學(xué)性及準(zhǔn)確性的防治建議。

    4)目前建立塌陷預(yù)測模型所采集的塌陷實際發(fā)生樣本均發(fā)生在終采6年之內(nèi),時間尺度較短,采空區(qū)頂板的殘余變形、破裂并未完全停止,后期深入研究可根據(jù)時間尺度更長的塌陷樣本完成礦區(qū)地面塌陷預(yù)測模型。

    猜你喜歡
    易發(fā)礦區(qū)神經(jīng)網(wǎng)絡(luò)
    機用鎳鈦銼在乳磨牙根管治療中的應(yīng)用
    貴州省地質(zhì)災(zāi)害易發(fā)分區(qū)圖
    夏季羊易發(fā)疾病及防治方法
    冬季雞腸炎易發(fā) 科學(xué)防治有方法
    加納Amanforom礦區(qū)Ⅲ號隱伏金礦帶的發(fā)現(xiàn)與評價
    加納Amanforom礦區(qū)Ⅲ號隱伏金礦帶的發(fā)現(xiàn)與評價
    湖北省??悼h堰邊上礦區(qū)發(fā)現(xiàn)超大型磷礦
    廣東省蕉嶺縣作壁坑礦區(qū)探明超大型銣礦
    神經(jīng)網(wǎng)絡(luò)抑制無線通信干擾探究
    電子制作(2019年19期)2019-11-23 08:42:00
    基于神經(jīng)網(wǎng)絡(luò)的拉矯機控制模型建立
    重型機械(2016年1期)2016-03-01 03:42:04
    亚洲,欧美精品.| 午夜精品一区二区三区免费看| 欧美日韩中文字幕国产精品一区二区三区| 久久久久久九九精品二区国产| 国产成+人综合+亚洲专区| 欧美极品一区二区三区四区| 欧美zozozo另类| 久久99热这里只有精品18| 大型黄色视频在线免费观看| 观看免费一级毛片| 亚洲av电影不卡..在线观看| 性色avwww在线观看| 一个人观看的视频www高清免费观看| 日韩欧美免费精品| 国产精品精品国产色婷婷| 国产国拍精品亚洲av在线观看 | 亚洲 国产 在线| 99国产极品粉嫩在线观看| 91在线精品国自产拍蜜月 | 超碰av人人做人人爽久久 | 99riav亚洲国产免费| 变态另类丝袜制服| 亚洲精品成人久久久久久| 97碰自拍视频| 婷婷精品国产亚洲av在线| 成人国产综合亚洲| 国产午夜福利久久久久久| 亚洲人成网站高清观看| 国产亚洲欧美98| 久久九九热精品免费| 男女之事视频高清在线观看| 亚洲av电影在线进入| 国产精品一区二区三区四区久久| www日本黄色视频网| 国产高清视频在线播放一区| 老司机在亚洲福利影院| 国产亚洲精品一区二区www| 亚洲在线观看片| 亚洲av不卡在线观看| av女优亚洲男人天堂| 在线免费观看不下载黄p国产 | 国产高清视频在线播放一区| 国语自产精品视频在线第100页| 999久久久精品免费观看国产| 国产毛片a区久久久久| 亚洲欧美日韩东京热| 老司机深夜福利视频在线观看| 亚洲午夜理论影院| 国产av麻豆久久久久久久| 日日摸夜夜添夜夜添小说| 日韩欧美国产在线观看| 国产成人啪精品午夜网站| 国产激情偷乱视频一区二区| 日韩大尺度精品在线看网址| 亚洲精品久久国产高清桃花| 制服丝袜大香蕉在线| 男女之事视频高清在线观看| 美女黄网站色视频| 五月伊人婷婷丁香| 91麻豆精品激情在线观看国产| 久久精品国产亚洲av香蕉五月| 热99在线观看视频| 久久99热这里只有精品18| 人妻夜夜爽99麻豆av| 成年女人看的毛片在线观看| 999久久久精品免费观看国产| 小说图片视频综合网站| 色老头精品视频在线观看| 一区二区三区国产精品乱码| xxx96com| 九色国产91popny在线| 国产欧美日韩精品一区二区| 观看免费一级毛片| 日韩欧美国产在线观看| 一卡2卡三卡四卡精品乱码亚洲| 很黄的视频免费| 久久性视频一级片| 91久久精品国产一区二区成人 | or卡值多少钱| 亚洲av美国av| 日韩欧美一区二区三区在线观看| 一进一出好大好爽视频| 亚洲成人中文字幕在线播放| 美女免费视频网站| 欧美av亚洲av综合av国产av| 国产午夜精品久久久久久一区二区三区 | 夜夜躁狠狠躁天天躁| 日本熟妇午夜| 国产高潮美女av| 亚洲欧美日韩无卡精品| 国产欧美日韩一区二区三| 男女下面进入的视频免费午夜| 亚洲精品日韩av片在线观看 | 亚洲av第一区精品v没综合| 在线观看午夜福利视频| 欧美日韩福利视频一区二区| 国产精品亚洲av一区麻豆| 成人亚洲精品av一区二区| 在线看三级毛片| 久久久精品大字幕| 成年版毛片免费区| 人妻夜夜爽99麻豆av| 色在线成人网| 黄色女人牲交| 久久久久久大精品| 精品免费久久久久久久清纯| 欧美激情在线99| 最近最新免费中文字幕在线| 国产淫片久久久久久久久 | 极品教师在线免费播放| 亚洲精品影视一区二区三区av| 国产中年淑女户外野战色| 国产精品永久免费网站| 国产亚洲精品av在线| 一进一出抽搐gif免费好疼| avwww免费| 有码 亚洲区| 成人无遮挡网站| 青草久久国产| 老司机午夜十八禁免费视频| 亚洲天堂国产精品一区在线| 日本熟妇午夜| 一区二区三区高清视频在线| www日本黄色视频网| 欧美性猛交黑人性爽| а√天堂www在线а√下载| 日韩高清综合在线| aaaaa片日本免费| 国产在视频线在精品| 好看av亚洲va欧美ⅴa在| 成人特级黄色片久久久久久久| av中文乱码字幕在线| 免费电影在线观看免费观看| 亚洲狠狠婷婷综合久久图片| 午夜老司机福利剧场| 亚洲成a人片在线一区二区| 久久伊人香网站| 日日摸夜夜添夜夜添小说| 亚洲色图av天堂| 伊人久久大香线蕉亚洲五| eeuss影院久久| 国产精品免费一区二区三区在线| 超碰av人人做人人爽久久 | 女人被狂操c到高潮| 怎么达到女性高潮| 精品不卡国产一区二区三区| 久久99热这里只有精品18| 亚洲七黄色美女视频| 别揉我奶头~嗯~啊~动态视频| 国产成人av激情在线播放| 12—13女人毛片做爰片一| 97超视频在线观看视频| 久久6这里有精品| 亚洲一区高清亚洲精品| 亚洲成av人片在线播放无| 国产精品综合久久久久久久免费| 久久精品国产自在天天线| 欧美成人性av电影在线观看| 99精品久久久久人妻精品| 色尼玛亚洲综合影院| 国产欧美日韩精品一区二区| 亚洲精品一卡2卡三卡4卡5卡| 2021天堂中文幕一二区在线观| 哪里可以看免费的av片| 免费观看的影片在线观看| 色尼玛亚洲综合影院| 欧美日韩精品网址| 亚洲片人在线观看| 国产高清激情床上av| 波野结衣二区三区在线 | 欧美日韩国产亚洲二区| 悠悠久久av| 久久精品国产清高在天天线| а√天堂www在线а√下载| 成年人黄色毛片网站| 人妻久久中文字幕网| 中文资源天堂在线| 日韩精品青青久久久久久| 老汉色av国产亚洲站长工具| 免费搜索国产男女视频| 一进一出抽搐动态| 国产综合懂色| 成年版毛片免费区| 国产男靠女视频免费网站| 有码 亚洲区| 毛片女人毛片| 欧美一级毛片孕妇| 天天一区二区日本电影三级| 国产三级在线视频| 亚洲精品美女久久久久99蜜臀| 亚洲五月天丁香| 亚洲一区高清亚洲精品| 一个人看视频在线观看www免费 | 国产久久久一区二区三区| 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| 日韩亚洲欧美综合| 精品久久久久久久人妻蜜臀av| xxx96com| 国内精品美女久久久久久| 亚洲av电影在线进入| 特大巨黑吊av在线直播| 全区人妻精品视频| 久久精品影院6| 18禁黄网站禁片午夜丰满| 高潮久久久久久久久久久不卡| 久久香蕉精品热| 国产私拍福利视频在线观看| 毛片女人毛片| 精品久久久久久久久久免费视频| 全区人妻精品视频| 中文字幕人妻丝袜一区二区| 国产高清有码在线观看视频| 少妇人妻精品综合一区二区 | 国产精品香港三级国产av潘金莲| 母亲3免费完整高清在线观看| 成熟少妇高潮喷水视频| 91九色精品人成在线观看| 日韩成人在线观看一区二区三区| 欧美性猛交╳xxx乱大交人| 色播亚洲综合网| 国产成人aa在线观看| 国内久久婷婷六月综合欲色啪| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利18| 久久久久久国产a免费观看| 久99久视频精品免费| 精品久久久久久久人妻蜜臀av| 日本黄大片高清| 亚洲激情在线av| 最后的刺客免费高清国语| 精品久久久久久久久久免费视频| 精品人妻偷拍中文字幕| 88av欧美| 久久午夜亚洲精品久久| 国产精品野战在线观看| 久久欧美精品欧美久久欧美| 哪里可以看免费的av片| 亚洲aⅴ乱码一区二区在线播放| 一级作爱视频免费观看| 国产成人a区在线观看| 伊人久久大香线蕉亚洲五| АⅤ资源中文在线天堂| 久久久国产成人免费| 男女那种视频在线观看| 999久久久精品免费观看国产| 久久精品国产亚洲av香蕉五月| 女人十人毛片免费观看3o分钟| 久久草成人影院| 午夜视频国产福利| 国内少妇人妻偷人精品xxx网站| 一进一出好大好爽视频| 91在线精品国自产拍蜜月 | 久久久久久九九精品二区国产| 99精品欧美一区二区三区四区| 久久性视频一级片| 国产成人a区在线观看| 欧美最新免费一区二区三区 | 亚洲成人免费电影在线观看| 少妇丰满av| 叶爱在线成人免费视频播放| 在线免费观看不下载黄p国产 | 免费看十八禁软件| 亚洲国产中文字幕在线视频| 又紧又爽又黄一区二区| 亚洲国产欧美人成| 成人鲁丝片一二三区免费| 国产亚洲精品av在线| 午夜免费成人在线视频| 偷拍熟女少妇极品色| 日韩欧美在线二视频| 亚洲人成网站在线播放欧美日韩| 最新中文字幕久久久久| 欧美一区二区亚洲| 日韩亚洲欧美综合| 亚洲午夜理论影院| 国产精品久久久久久亚洲av鲁大| 久久久久久大精品| 黄色片一级片一级黄色片| 国产伦精品一区二区三区四那| 久久精品人妻少妇| 日韩欧美免费精品| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 中文字幕人妻熟人妻熟丝袜美 | 欧美xxxx黑人xx丫x性爽| 精品人妻一区二区三区麻豆 | 99视频精品全部免费 在线| 国产高清视频在线观看网站| 亚洲成a人片在线一区二区| 欧美日韩国产亚洲二区| av天堂在线播放| 欧美一区二区亚洲| 国产精品久久电影中文字幕| 亚洲 欧美 日韩 在线 免费| 国产伦人伦偷精品视频| 法律面前人人平等表现在哪些方面| 国产老妇女一区| 天堂动漫精品| 老司机福利观看| 999久久久精品免费观看国产| 最近视频中文字幕2019在线8| 国产单亲对白刺激| 欧美日韩精品网址| 88av欧美| 成人av在线播放网站| 成人永久免费在线观看视频| 搡老岳熟女国产| 成人无遮挡网站| 国产成人av激情在线播放| 三级国产精品欧美在线观看| 黄片小视频在线播放| 久久九九热精品免费| 在线国产一区二区在线| 亚洲自拍偷在线| 欧美黑人巨大hd| 最近最新中文字幕大全电影3| 国产一级毛片七仙女欲春2| 国产成年人精品一区二区| 在线观看日韩欧美| 国产探花在线观看一区二区| 91字幕亚洲| 嫩草影院精品99| 亚洲第一欧美日韩一区二区三区| 一个人看的www免费观看视频| 国产精品一及| tocl精华| 午夜免费男女啪啪视频观看 | 一个人免费在线观看的高清视频| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 极品教师在线免费播放| 最近最新中文字幕大全电影3| 亚洲av成人不卡在线观看播放网| 国产亚洲欧美在线一区二区| 少妇的逼好多水| 成人特级av手机在线观看| 国产精品 欧美亚洲| 午夜日韩欧美国产| 国产三级中文精品| 极品教师在线免费播放| 亚洲男人的天堂狠狠| 夜夜看夜夜爽夜夜摸| 深爱激情五月婷婷| 免费在线观看亚洲国产| 欧美乱码精品一区二区三区| 丁香欧美五月| 极品教师在线免费播放| 两个人看的免费小视频| 欧美+亚洲+日韩+国产| 嫩草影院精品99| 亚洲国产精品久久男人天堂| 亚洲专区国产一区二区| 黄片小视频在线播放| 成人性生交大片免费视频hd| 国产亚洲精品综合一区在线观看| 少妇丰满av| 99在线视频只有这里精品首页| 午夜免费男女啪啪视频观看 | 丰满人妻一区二区三区视频av | netflix在线观看网站| 一级a爱片免费观看的视频| 99riav亚洲国产免费| 波多野结衣巨乳人妻| 久久精品国产亚洲av香蕉五月| 一级黄片播放器| 伊人久久精品亚洲午夜| 深夜精品福利| 国产v大片淫在线免费观看| 精品电影一区二区在线| 性色av乱码一区二区三区2| 看免费av毛片| 午夜福利欧美成人| 国产精华一区二区三区| 怎么达到女性高潮| 国产真实伦视频高清在线观看 | 国产精品影院久久| 俺也久久电影网| 国产熟女xx| 色噜噜av男人的天堂激情| 黄片大片在线免费观看| 成人欧美大片| av国产免费在线观看| 精品欧美国产一区二区三| 亚洲成av人片在线播放无| 亚洲精品日韩av片在线观看 | 日本一本二区三区精品| 最好的美女福利视频网| 久久6这里有精品| 99热精品在线国产| 欧美av亚洲av综合av国产av| 国产精品嫩草影院av在线观看 | 男女之事视频高清在线观看| 又黄又爽又免费观看的视频| 中亚洲国语对白在线视频| 黄色成人免费大全| 男插女下体视频免费在线播放| 动漫黄色视频在线观看| 国产亚洲精品久久久com| 国产99白浆流出| 精品欧美国产一区二区三| 99久久久亚洲精品蜜臀av| 日韩有码中文字幕| 欧美日本亚洲视频在线播放| 人妻久久中文字幕网| 午夜a级毛片| 亚洲第一电影网av| 美女高潮喷水抽搐中文字幕| 日日摸夜夜添夜夜添小说| 综合色av麻豆| 内射极品少妇av片p| 69人妻影院| 老司机在亚洲福利影院| 日本一二三区视频观看| 国产精品综合久久久久久久免费| 久久久久性生活片| 国产一区二区在线观看日韩 | 无遮挡黄片免费观看| 老司机午夜福利在线观看视频| 少妇丰满av| 精品久久久久久久久久久久久| 日韩欧美三级三区| 丁香六月欧美| 国内毛片毛片毛片毛片毛片| 国产精品国产高清国产av| 欧美成人一区二区免费高清观看| 欧美日韩黄片免| 别揉我奶头~嗯~啊~动态视频| 97碰自拍视频| 老司机午夜福利在线观看视频| 男女下面进入的视频免费午夜| 亚洲第一电影网av| 国产国拍精品亚洲av在线观看 | 欧美黑人欧美精品刺激| 久久久久久久久中文| 欧美性感艳星| 免费无遮挡裸体视频| 日韩高清综合在线| 99久久久亚洲精品蜜臀av| 一级黄片播放器| 久久久久久九九精品二区国产| 国产精华一区二区三区| 国产亚洲精品久久久久久毛片| 久久亚洲真实| 18+在线观看网站| 女生性感内裤真人,穿戴方法视频| 黑人欧美特级aaaaaa片| 久久久久国产精品人妻aⅴ院| av欧美777| 禁无遮挡网站| 亚洲中文日韩欧美视频| 99久久精品热视频| 一卡2卡三卡四卡精品乱码亚洲| 热99re8久久精品国产| 国产午夜福利久久久久久| 久久天躁狠狠躁夜夜2o2o| 天天一区二区日本电影三级| 极品教师在线免费播放| 中文字幕久久专区| www.999成人在线观看| 黄色丝袜av网址大全| 免费高清视频大片| 韩国av一区二区三区四区| а√天堂www在线а√下载| 欧美日韩综合久久久久久 | 日韩欧美三级三区| 看片在线看免费视频| 一夜夜www| 三级男女做爰猛烈吃奶摸视频| 最近在线观看免费完整版| 亚洲精品一区av在线观看| 99热这里只有精品一区| 长腿黑丝高跟| 国产探花极品一区二区| 国产精华一区二区三区| 亚洲精品一区av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲内射少妇av| 久久精品亚洲精品国产色婷小说| www.www免费av| 日韩免费av在线播放| 亚洲精品456在线播放app | 久久精品人妻少妇| 国产一级毛片七仙女欲春2| 欧美国产日韩亚洲一区| 91麻豆av在线| 怎么达到女性高潮| 每晚都被弄得嗷嗷叫到高潮| 久久久久久久久久黄片| 在线免费观看不下载黄p国产 | 欧美一级毛片孕妇| 999久久久精品免费观看国产| 制服丝袜大香蕉在线| 国产成人av激情在线播放| 3wmmmm亚洲av在线观看| 99在线视频只有这里精品首页| 国内精品久久久久精免费| 精品人妻1区二区| 日韩亚洲欧美综合| 搡老熟女国产l中国老女人| 亚洲av免费在线观看| 国产97色在线日韩免费| 亚洲精品在线观看二区| 麻豆成人午夜福利视频| 欧美日韩综合久久久久久 | 国产亚洲欧美98| 757午夜福利合集在线观看| 午夜影院日韩av| 久久6这里有精品| 亚洲第一电影网av| 男女午夜视频在线观看| 午夜亚洲福利在线播放| 亚洲第一电影网av| 在线观看av片永久免费下载| 成人特级av手机在线观看| 亚洲第一欧美日韩一区二区三区| 搡老熟女国产l中国老女人| 国产精品影院久久| 99精品欧美一区二区三区四区| 免费搜索国产男女视频| 亚洲精品在线观看二区| 国产精品自产拍在线观看55亚洲| 一本一本综合久久| 欧美大码av| 亚洲av五月六月丁香网| 丁香欧美五月| 99久国产av精品| 久久精品国产亚洲av涩爱 | 尤物成人国产欧美一区二区三区| 亚洲成av人片在线播放无| 好男人电影高清在线观看| 亚洲激情在线av| 一夜夜www| 99热这里只有是精品50| 黄色女人牲交| 最近在线观看免费完整版| 两个人的视频大全免费| 最后的刺客免费高清国语| 国产精品国产高清国产av| 国产私拍福利视频在线观看| 精品无人区乱码1区二区| 国产伦一二天堂av在线观看| 国产高清激情床上av| 午夜福利视频1000在线观看| xxx96com| 亚洲成人久久性| 男女视频在线观看网站免费| 国产毛片a区久久久久| 成人特级av手机在线观看| netflix在线观看网站| 精品乱码久久久久久99久播| 精品久久久久久久久久久久久| 热99在线观看视频| 亚洲av电影不卡..在线观看| 99国产极品粉嫩在线观看| aaaaa片日本免费| 小说图片视频综合网站| 免费搜索国产男女视频| 成年版毛片免费区| а√天堂www在线а√下载| 亚洲国产欧美网| www.999成人在线观看| 精品久久久久久成人av| 亚洲精品亚洲一区二区| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久久电影 | 在线观看舔阴道视频| 欧美性感艳星| 成人欧美大片| 长腿黑丝高跟| 国产精品 欧美亚洲| 免费av不卡在线播放| 岛国视频午夜一区免费看| 久久草成人影院| 又粗又爽又猛毛片免费看| 免费搜索国产男女视频| 一本综合久久免费| 琪琪午夜伦伦电影理论片6080| 淫妇啪啪啪对白视频| 少妇人妻一区二区三区视频| 亚洲国产色片| 日韩中文字幕欧美一区二区| 久久亚洲真实| 757午夜福利合集在线观看| а√天堂www在线а√下载| 久久精品国产亚洲av涩爱 | 一卡2卡三卡四卡精品乱码亚洲| 嫩草影视91久久| 性色avwww在线观看| 波野结衣二区三区在线 | 熟女人妻精品中文字幕| 国产一级毛片七仙女欲春2| 精品久久久久久久末码| 日韩欧美在线乱码| 国产91精品成人一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲专区中文字幕在线| 美女cb高潮喷水在线观看| 国产乱人伦免费视频| 欧美性猛交╳xxx乱大交人| 精品久久久久久久久久久久久| 日本与韩国留学比较| 黄色片一级片一级黄色片| 国产av在哪里看| 成人亚洲精品av一区二区| 伊人久久大香线蕉亚洲五| 美女cb高潮喷水在线观看| 亚洲av五月六月丁香网| 十八禁人妻一区二区| 1000部很黄的大片| 麻豆成人午夜福利视频| 每晚都被弄得嗷嗷叫到高潮| а√天堂www在线а√下载| 成人欧美大片| 男女做爰动态图高潮gif福利片| 中出人妻视频一区二区|