• <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
    美女 人体艺术 gogo| 欧美精品啪啪一区二区三区| 欧美色欧美亚洲另类二区| 国产成人av教育| xxxwww97欧美| 首页视频小说图片口味搜索| 三级毛片av免费| 亚洲午夜理论影院| 国产真实乱freesex| 十八禁人妻一区二区| 亚洲av第一区精品v没综合| 亚洲九九香蕉| 极品教师在线免费播放| 舔av片在线| 精品免费久久久久久久清纯| 亚洲性夜色夜夜综合| 香蕉丝袜av| 国产成人精品无人区| 日韩欧美免费精品| 禁无遮挡网站| 国产成人aa在线观看| 精品久久久久久成人av| 日本一本二区三区精品| 亚洲av中文字字幕乱码综合| 中文字幕av在线有码专区| 色综合婷婷激情| 黑人操中国人逼视频| 国产精品一区二区三区四区免费观看 | 国产午夜精品久久久久久| 国产亚洲欧美在线一区二区| 国产精品 国内视频| 日本一区二区免费在线视频| svipshipincom国产片| 看免费av毛片| 亚洲午夜理论影院| 757午夜福利合集在线观看| 妹子高潮喷水视频| 亚洲一区二区三区不卡视频| 国产久久久一区二区三区| а√天堂www在线а√下载| 欧美日韩亚洲综合一区二区三区_| 香蕉丝袜av| 老鸭窝网址在线观看| 十八禁人妻一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩一级在线毛片| 国产日本99.免费观看| 成年女人毛片免费观看观看9| 亚洲国产中文字幕在线视频| 亚洲欧洲精品一区二区精品久久久| 国产亚洲精品久久久久久毛片| 白带黄色成豆腐渣| 动漫黄色视频在线观看| ponron亚洲| 午夜免费成人在线视频| 国产又色又爽无遮挡免费看| 久久精品国产亚洲av香蕉五月| 91麻豆精品激情在线观看国产| 亚洲熟妇中文字幕五十中出| 亚洲国产欧洲综合997久久,| 99在线人妻在线中文字幕| 国产精品自产拍在线观看55亚洲| 极品教师在线免费播放| 色精品久久人妻99蜜桃| 两性夫妻黄色片| 亚洲人成电影免费在线| 精品无人区乱码1区二区| 久久草成人影院| 国产av一区在线观看免费| 欧美一级a爱片免费观看看 | 伊人久久大香线蕉亚洲五| 亚洲片人在线观看| 国产精品一区二区三区四区久久| 亚洲精品国产精品久久久不卡| 精品久久久久久久毛片微露脸| 国产精品自产拍在线观看55亚洲| 欧美成人性av电影在线观看| 高清毛片免费观看视频网站| 在线观看日韩欧美| 97超级碰碰碰精品色视频在线观看| 成人三级做爰电影| 怎么达到女性高潮| 两性午夜刺激爽爽歪歪视频在线观看 | 日本五十路高清| 久久精品国产清高在天天线| 亚洲无线在线观看| 欧美日韩乱码在线| 国模一区二区三区四区视频 | 青草久久国产| 在线看三级毛片| 亚洲男人的天堂狠狠| 久久午夜综合久久蜜桃| 成人欧美大片| 熟妇人妻久久中文字幕3abv| 国产亚洲av高清不卡| 一进一出抽搐gif免费好疼| 日韩欧美精品v在线| 日韩精品免费视频一区二区三区| 国产精品一区二区三区四区久久| 少妇裸体淫交视频免费看高清 | 亚洲国产精品久久男人天堂| 亚洲天堂国产精品一区在线| 国产精品野战在线观看| 9191精品国产免费久久| 亚洲最大成人中文| 夜夜看夜夜爽夜夜摸| 久久久久久九九精品二区国产 | 热99re8久久精品国产| 特大巨黑吊av在线直播| 99热这里只有是精品50| 国产亚洲精品综合一区在线观看 | 最近最新中文字幕大全免费视频| 欧美在线一区亚洲| 欧美色视频一区免费| 欧美中文日本在线观看视频| 最近最新中文字幕大全电影3| 亚洲国产日韩欧美精品在线观看 | a级毛片在线看网站| 精品欧美国产一区二区三| 狂野欧美激情性xxxx| 脱女人内裤的视频| 亚洲国产精品合色在线| 久久婷婷成人综合色麻豆| 精华霜和精华液先用哪个| 亚洲18禁久久av| 一个人免费在线观看电影 | 18禁裸乳无遮挡免费网站照片| 美女 人体艺术 gogo| 午夜福利成人在线免费观看| 91老司机精品| 国产伦在线观看视频一区| 国产精品亚洲一级av第二区| 91字幕亚洲| 国内揄拍国产精品人妻在线| 一二三四社区在线视频社区8| 亚洲欧美一区二区三区黑人| 两性夫妻黄色片| 国产v大片淫在线免费观看| 国产探花在线观看一区二区| 日韩国内少妇激情av| 熟女电影av网| 最好的美女福利视频网| 久久精品国产清高在天天线| 亚洲av成人av| 亚洲国产精品999在线| 丁香六月欧美| 在线观看美女被高潮喷水网站 | 在线观看一区二区三区| 看免费av毛片| 女人爽到高潮嗷嗷叫在线视频| 在线观看66精品国产| 国产成人系列免费观看| 十八禁网站免费在线| av超薄肉色丝袜交足视频| 无人区码免费观看不卡| 欧美日韩亚洲国产一区二区在线观看| xxx96com| 日韩精品中文字幕看吧| 神马国产精品三级电影在线观看 | 国产精品永久免费网站| 看片在线看免费视频| 亚洲熟女毛片儿| 欧美在线黄色| 亚洲国产欧美人成| a在线观看视频网站| 女警被强在线播放| 99久久综合精品五月天人人| 中文字幕久久专区| 亚洲精品国产一区二区精华液| 老司机午夜福利在线观看视频| 怎么达到女性高潮| 国产精品一区二区精品视频观看| 久久久久久久午夜电影| 听说在线观看完整版免费高清| 制服诱惑二区| 亚洲国产中文字幕在线视频| 国产精品亚洲美女久久久| 国产精品久久久久久久电影 | 丝袜人妻中文字幕| 国产视频内射| 国产成人av激情在线播放| 国内精品一区二区在线观看| 久久人人精品亚洲av| 视频区欧美日本亚洲| 日韩大码丰满熟妇| √禁漫天堂资源中文www| 国产伦一二天堂av在线观看| 人妻夜夜爽99麻豆av| 两个人免费观看高清视频| 亚洲欧美一区二区三区黑人| 欧美丝袜亚洲另类 | 国产精品精品国产色婷婷| 午夜福利视频1000在线观看| 亚洲精品国产一区二区精华液| 在线观看舔阴道视频| 久久久久久九九精品二区国产 | 久久久久精品国产欧美久久久| 国产伦一二天堂av在线观看| 伊人久久大香线蕉亚洲五| 国产精品一区二区三区四区久久| 嫁个100分男人电影在线观看| 午夜福利成人在线免费观看| 亚洲中文日韩欧美视频| 欧美精品亚洲一区二区| av国产免费在线观看| 一个人免费在线观看的高清视频| 给我免费播放毛片高清在线观看| 在线播放国产精品三级| 亚洲欧美日韩东京热| 97碰自拍视频| 哪里可以看免费的av片| 一边摸一边抽搐一进一小说| 精品午夜福利视频在线观看一区| 欧美乱色亚洲激情| 国产三级中文精品| 亚洲精品美女久久av网站| 99在线人妻在线中文字幕| 久久天躁狠狠躁夜夜2o2o| 成人国语在线视频| 成人永久免费在线观看视频| 日韩欧美 国产精品| 变态另类丝袜制服| 后天国语完整版免费观看| 精品久久久久久久久久免费视频| 国产精品久久久久久久电影 | 91大片在线观看| 色老头精品视频在线观看| 男女午夜视频在线观看| 深夜精品福利| 深夜精品福利| 香蕉国产在线看| 久久国产精品人妻蜜桃| 午夜激情av网站| av欧美777| 99久久无色码亚洲精品果冻| 丰满人妻熟妇乱又伦精品不卡| 午夜激情av网站| 国语自产精品视频在线第100页| 亚洲午夜理论影院| 91字幕亚洲| 亚洲av日韩精品久久久久久密| 免费在线观看黄色视频的| 一级毛片女人18水好多| 黑人巨大精品欧美一区二区mp4| 一a级毛片在线观看| cao死你这个sao货| 免费在线观看成人毛片| 久久精品aⅴ一区二区三区四区| 男女床上黄色一级片免费看| 丝袜美腿诱惑在线| 欧美日韩福利视频一区二区| 99国产精品99久久久久| 18禁观看日本| 99在线人妻在线中文字幕| 午夜精品久久久久久毛片777| 国内精品久久久久精免费| 精品一区二区三区视频在线观看免费| 国产主播在线观看一区二区| cao死你这个sao货| 中文字幕最新亚洲高清| 动漫黄色视频在线观看| 波多野结衣高清作品| 久久午夜亚洲精品久久| 午夜精品一区二区三区免费看| 国产高清视频在线播放一区| 国产精品爽爽va在线观看网站| 热99re8久久精品国产| 亚洲av成人一区二区三| 色在线成人网| 制服丝袜大香蕉在线| 少妇的丰满在线观看| 在线观看舔阴道视频| 亚洲国产精品sss在线观看| а√天堂www在线а√下载| 高清在线国产一区| 老司机靠b影院| 琪琪午夜伦伦电影理论片6080| 在线观看日韩欧美| 精品高清国产在线一区| 亚洲乱码一区二区免费版| 搡老妇女老女人老熟妇| 首页视频小说图片口味搜索| 欧美日韩瑟瑟在线播放| 久久这里只有精品中国| 毛片女人毛片| 一区二区三区国产精品乱码| 国产爱豆传媒在线观看 | 国产在线精品亚洲第一网站| 亚洲国产高清在线一区二区三| 日韩三级视频一区二区三区| 国产熟女午夜一区二区三区| 看免费av毛片| 日韩精品青青久久久久久| 极品教师在线免费播放| 亚洲免费av在线视频| 又黄又粗又硬又大视频| 美女大奶头视频| 中国美女看黄片| 亚洲欧美日韩高清在线视频| 午夜两性在线视频| 欧美一区二区国产精品久久精品 | 中文字幕高清在线视频| 99久久久亚洲精品蜜臀av| 亚洲国产欧洲综合997久久,| 99精品久久久久人妻精品| 99国产精品99久久久久| 村上凉子中文字幕在线| 欧美zozozo另类| 国产精品 欧美亚洲| 妹子高潮喷水视频| 久久久久久国产a免费观看| 中出人妻视频一区二区| 午夜两性在线视频| 久久香蕉国产精品| 桃色一区二区三区在线观看| 曰老女人黄片| 美女高潮喷水抽搐中文字幕| 精品一区二区三区四区五区乱码| 亚洲av中文字字幕乱码综合| 国产三级中文精品| 午夜亚洲福利在线播放| 最新美女视频免费是黄的| 少妇粗大呻吟视频| 狠狠狠狠99中文字幕| 两个人看的免费小视频| 欧美日本亚洲视频在线播放| 久久久久国内视频| 国产探花在线观看一区二区| 人人妻人人澡欧美一区二区| 夜夜爽天天搞| 久久国产精品人妻蜜桃| 亚洲av成人不卡在线观看播放网| 制服丝袜大香蕉在线| 级片在线观看| 男人舔女人的私密视频| 国产精品精品国产色婷婷| 老司机福利观看| 小说图片视频综合网站| 久久精品影院6| 欧美中文日本在线观看视频| 亚洲成人免费电影在线观看| 视频区欧美日本亚洲| 久久久久国产一级毛片高清牌| 又紧又爽又黄一区二区| 国产三级黄色录像| 欧美中文日本在线观看视频| 国产精品亚洲一级av第二区| 久久久水蜜桃国产精品网| 成人国语在线视频| 国产精品98久久久久久宅男小说| 亚洲国产欧洲综合997久久,| 99热6这里只有精品| 欧美绝顶高潮抽搐喷水| 黄色女人牲交| 麻豆国产av国片精品| 欧美高清成人免费视频www| 一a级毛片在线观看| 在线观看日韩欧美| 国产精品乱码一区二三区的特点| 他把我摸到了高潮在线观看| 黄色片一级片一级黄色片| 色老头精品视频在线观看| 波多野结衣高清作品| 女人被狂操c到高潮| 国产精品亚洲一级av第二区| 亚洲国产欧美人成| 特大巨黑吊av在线直播| 999久久久精品免费观看国产| 免费看美女性在线毛片视频| 日本在线视频免费播放| 久久久久久九九精品二区国产 | 看片在线看免费视频| www.999成人在线观看| 欧美黄色淫秽网站| 亚洲中文字幕日韩| 91av网站免费观看| 久久久久久久久免费视频了| 国产av在哪里看| 男女视频在线观看网站免费 | 亚洲精品美女久久av网站| 免费看日本二区| 小说图片视频综合网站| 午夜亚洲福利在线播放| 欧美另类亚洲清纯唯美| 午夜两性在线视频| 国产精品精品国产色婷婷| 久久久久国产一级毛片高清牌| 精品不卡国产一区二区三区| 久久久水蜜桃国产精品网| 女人被狂操c到高潮| 身体一侧抽搐| 一夜夜www| 夜夜躁狠狠躁天天躁| 国产一区在线观看成人免费| 制服诱惑二区| 亚洲一区中文字幕在线| 欧美一区二区国产精品久久精品 | 亚洲国产精品合色在线| 亚洲中文av在线| 伦理电影免费视频| 婷婷精品国产亚洲av在线| 午夜福利在线观看吧| 黑人操中国人逼视频| 丰满的人妻完整版| 国产亚洲欧美98| 国产亚洲精品久久久久久毛片| 国产精品精品国产色婷婷| 最近最新中文字幕大全电影3| 日本五十路高清| 久久天躁狠狠躁夜夜2o2o| 男男h啪啪无遮挡| 久久这里只有精品中国| 99热6这里只有精品| 在线观看www视频免费| 我要搜黄色片| 看免费av毛片| 国产主播在线观看一区二区| 欧美3d第一页| 成人高潮视频无遮挡免费网站| 哪里可以看免费的av片| 亚洲男人的天堂狠狠| 成人国产一区最新在线观看| 欧美黑人巨大hd| 在线观看舔阴道视频| 深夜精品福利| 欧美成人免费av一区二区三区| 久久久久久人人人人人| 一边摸一边抽搐一进一小说| 在线观看舔阴道视频| 亚洲男人天堂网一区| 婷婷丁香在线五月| 午夜福利高清视频| 中亚洲国语对白在线视频| 亚洲专区字幕在线| 国产成年人精品一区二区| 三级男女做爰猛烈吃奶摸视频| 久久精品亚洲精品国产色婷小说| 少妇被粗大的猛进出69影院| www国产在线视频色| 一进一出抽搐gif免费好疼| 久久精品91无色码中文字幕| 国产午夜精品论理片| 国产欧美日韩精品亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av成人不卡在线观看播放网| 日本一二三区视频观看| 不卡av一区二区三区| 亚洲av成人av| 亚洲男人的天堂狠狠| 亚洲美女黄片视频| 99在线视频只有这里精品首页| 色老头精品视频在线观看| 欧美黄色片欧美黄色片| 狂野欧美激情性xxxx| 伊人久久大香线蕉亚洲五| 十八禁人妻一区二区| 成人国产综合亚洲| 亚洲av电影在线进入| 精品一区二区三区av网在线观看| 国产又色又爽无遮挡免费看| 国产亚洲欧美98| 国产爱豆传媒在线观看 | 成年版毛片免费区| 91成年电影在线观看| 欧美一级a爱片免费观看看 | 欧美3d第一页| 午夜福利高清视频| 中亚洲国语对白在线视频| 欧美高清成人免费视频www| 免费在线观看完整版高清| 午夜福利成人在线免费观看| 人人妻人人澡欧美一区二区| 精品熟女少妇八av免费久了| 小说图片视频综合网站| av视频在线观看入口| 真人做人爱边吃奶动态| 亚洲av美国av| 久久久久精品国产欧美久久久| 巨乳人妻的诱惑在线观看| 免费在线观看黄色视频的| 亚洲成人免费电影在线观看| 色播亚洲综合网| 亚洲专区字幕在线| 一个人免费在线观看的高清视频| 一级毛片高清免费大全| 亚洲av成人av| 观看免费一级毛片| avwww免费| 中出人妻视频一区二区| 久久国产乱子伦精品免费另类| 国产成人av教育| 欧美+亚洲+日韩+国产| 亚洲成人久久爱视频| 欧美一区二区精品小视频在线| 免费在线观看黄色视频的| 日本免费一区二区三区高清不卡| 成人手机av| 日韩中文字幕欧美一区二区| 俺也久久电影网| 久久香蕉精品热| 久久久久九九精品影院| av福利片在线观看| av超薄肉色丝袜交足视频| 激情在线观看视频在线高清| 俄罗斯特黄特色一大片| 色综合婷婷激情| 国产亚洲欧美在线一区二区| 一级毛片女人18水好多| 亚洲国产看品久久| 日本撒尿小便嘘嘘汇集6| 国产一级毛片七仙女欲春2| 高清毛片免费观看视频网站| 黄色成人免费大全| 大型av网站在线播放| 人人妻人人看人人澡| 国产一区二区激情短视频| 99精品在免费线老司机午夜| 深夜精品福利| 又爽又黄无遮挡网站| 国产99白浆流出| 欧美黑人巨大hd| 亚洲精品一卡2卡三卡4卡5卡| av视频在线观看入口| 九色国产91popny在线| 又粗又爽又猛毛片免费看| 啪啪无遮挡十八禁网站| 两性午夜刺激爽爽歪歪视频在线观看 | 中出人妻视频一区二区| 欧美日本亚洲视频在线播放| 非洲黑人性xxxx精品又粗又长| 成人国产一区最新在线观看| 国产成人一区二区三区免费视频网站| 国产高清视频在线播放一区| 国产真实乱freesex| 欧美+亚洲+日韩+国产| 黄色片一级片一级黄色片| 亚洲av成人不卡在线观看播放网| 亚洲人与动物交配视频| ponron亚洲| 三级毛片av免费| 久久久久亚洲av毛片大全| 久久这里只有精品19| 女人高潮潮喷娇喘18禁视频| 久久久久性生活片| 2021天堂中文幕一二区在线观| 亚洲欧美激情综合另类| 女人高潮潮喷娇喘18禁视频| 国产单亲对白刺激| 亚洲精品在线观看二区| 国产区一区二久久| 久久人人精品亚洲av| 成年女人毛片免费观看观看9| 国产精品久久久av美女十八| 中出人妻视频一区二区| 一级a爱片免费观看的视频| 又爽又黄无遮挡网站| 美女午夜性视频免费| √禁漫天堂资源中文www| 黄色女人牲交| 两个人看的免费小视频| АⅤ资源中文在线天堂| 国产探花在线观看一区二区| 国产精品九九99| 欧美日本视频| 一夜夜www| 精品福利观看| 一边摸一边抽搐一进一小说| 熟女电影av网| 精品欧美国产一区二区三| ponron亚洲| 免费在线观看成人毛片| 亚洲成人精品中文字幕电影| 久久草成人影院| 国产成人系列免费观看| 午夜福利成人在线免费观看| 色噜噜av男人的天堂激情| 国产久久久一区二区三区| 欧美午夜高清在线| 久久人妻av系列| 床上黄色一级片| 欧美高清成人免费视频www| 制服丝袜大香蕉在线| 久久久久国产精品人妻aⅴ院| 久久久久久久精品吃奶| 欧美激情久久久久久爽电影| 成人午夜高清在线视频| 亚洲色图av天堂| 一边摸一边抽搐一进一小说| 黄色 视频免费看| 美女 人体艺术 gogo| 美女大奶头视频| 国产精品永久免费网站| 亚洲精品色激情综合| 91老司机精品| 久久久久久久午夜电影| 搡老熟女国产l中国老女人| 国产高清有码在线观看视频 | 在线观看一区二区三区| 国产人伦9x9x在线观看| 99在线人妻在线中文字幕| 一级毛片高清免费大全| 制服丝袜大香蕉在线| 又黄又粗又硬又大视频| 免费人成视频x8x8入口观看| 夜夜看夜夜爽夜夜摸| 免费在线观看黄色视频的| 男女午夜视频在线观看| 欧美另类亚洲清纯唯美| 国产主播在线观看一区二区| 精品久久久久久久久久免费视频| 五月玫瑰六月丁香| 黑人操中国人逼视频| 黄频高清免费视频| 国产精品1区2区在线观看.| 午夜精品在线福利|