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

    高潛水位采煤沉陷區(qū)積水范圍動態(tài)演化規(guī)律

    2020-06-05 02:11:34陳曉謝張文濤朱曉峻孫愛國楊曉玉
    煤田地質(zhì)與勘探 2020年2期
    關(guān)鍵詞:增長量采動積水

    陳曉謝,張文濤,朱曉峻,劉 輝,4,張 敏,孫愛國,楊曉玉

    高潛水位采煤沉陷區(qū)積水范圍動態(tài)演化規(guī)律

    陳曉謝1,2,張文濤3,朱曉峻1,2,劉 輝1,2,4,張 敏1,2,孫愛國5,楊曉玉6

    (1. 安徽大學(xué) 資源與環(huán)境工程學(xué)院,安徽 合肥 230601;2. 安徽省礦山生態(tài)修復(fù)工程實驗室,安徽 合肥 230601;3. 安徽理工大學(xué) 地球與環(huán)境學(xué)院,安徽 淮南 232001;4. 河北工程大學(xué) 礦業(yè)與測繪學(xué)院,河北 邯鄲 056038;5. 安徽恒源煤電股份有限公司五溝煤礦,安徽 淮北 235131;6. 合肥財經(jīng)職業(yè)學(xué)院,安徽 合肥 230601)

    為了研究高潛水位采煤沉陷區(qū)積水范圍動態(tài)演化規(guī)律,以安徽淮北五溝煤礦1031工作面為研究區(qū)域,基于Landsat-8遙感數(shù)據(jù)解譯結(jié)果,總結(jié)了2013—2017年1031工作面上方地表積水范圍演化規(guī)律,分析凈降水量、地下水埋深、工作面推進(jìn)距離等因素對積水范圍演化的影響,并提出積水邊界角的概念,建立了積水邊界角隨采動時間變化的函數(shù)關(guān)系式。研究發(fā)現(xiàn):高潛水位采煤沉陷區(qū)地表積水演化分為4個階段:未形成期、同步增長期、殘余增長期、相對穩(wěn)定期;在同步增長期和殘余增長期,工作面推進(jìn)距離、凈降水量與積水面積日增長量呈正相關(guān),地下水埋深與積水面積日增長量呈負(fù)相關(guān);工作面推進(jìn)距離是積水范圍演化的關(guān)鍵影響因素,在工作面推進(jìn)距離為476 m左右時,地表產(chǎn)生積水,工作面停采后5個月左右,積水面積趨于穩(wěn)定;積水邊界角先呈減小趨勢,再趨于穩(wěn)定,在同步增長期和殘余增長期,積水邊界角總體呈減小趨勢,但受氣象、地下水埋深等因素影響而波動;相對穩(wěn)定期,積水邊界角整體趨于90°。本研究為高潛水位采煤沉陷積水區(qū)土地利用規(guī)劃、土地復(fù)墾、水陸復(fù)合生態(tài)系統(tǒng)建立等提供理論依據(jù)。

    高潛水位;采煤沉陷區(qū);遙感監(jiān)測;積水演化階段;積水邊界角;五溝煤礦

    煤炭資源在我國能源體系中占非常重要的地位,占國家能源消費總量的64%[1],煤炭的大規(guī)模開采造成了嚴(yán)重的地表塌陷、地裂縫等地質(zhì)災(zāi)害[2-3]。兩淮地區(qū)是我國重要的煤炭基地,該地區(qū)地下潛水位高、松散層厚、下沉系數(shù)大,煤炭資源的開采引起地表下沉形成盆地,加上大氣降水與地下水補給等因素的共同影響,沉陷盆地內(nèi)形成大面積的積水,導(dǎo)致沉陷區(qū)水土流失嚴(yán)重,村莊、耕地被積水淹沒,原有的陸地生態(tài)系統(tǒng)轉(zhuǎn)化為水陸復(fù)合生態(tài)系統(tǒng)[4-5],對當(dāng)?shù)厣鷳B(tài)環(huán)境、社會經(jīng)濟造成了嚴(yán)重的影響。

    關(guān)于沉陷區(qū)積水演化研究,目前主要集中在積水面積演化和水循環(huán)兩個方面。在積水演化方面,彭蘇萍等[6]以淮南市為例,利用多時相TM遙感影像監(jiān)測地表積水塌陷區(qū)動態(tài)變化信息;魏礦靈等[7]以淮南礦區(qū)為例,通過監(jiān)測不同時期SPOT-5遙感影像,確定開采沉陷的位置及積水變化范圍;馬建威等[8]利用Landsat系列衛(wèi)星數(shù)據(jù),研究了武漢地區(qū)1973—2015年的城區(qū)水面范圍變化,分析多種因素對水體面積的影響等。在水循環(huán)方面,陸垂裕等[9-10]提出了面向?qū)ο竽K化的分布式水文模型MODCYCLE,并借助相關(guān)水文推理和數(shù)值模擬分析,對淮南典型高潛水位采煤沉陷區(qū)的地下水機制進(jìn)行辨析,得到高潛水位采煤沉陷區(qū)積水的補給來源和排泄渠道;李慧等[11-12]利用分布式水文模型,深入研究東部高潛水位采煤沉陷積水區(qū)庫容變化的主要影響因素及其相互關(guān)系。

    針對高潛水位采煤沉陷區(qū)積水演化規(guī)律研究,多針對地表或淺層地表的積水范圍,結(jié)合工作面采動影響的積水動態(tài)演化規(guī)律方面研究較少。本文以五溝煤礦1031工作面為研究區(qū),基于Landsat-8衛(wèi)星遙感數(shù)據(jù),對1031工作面地表積水范圍進(jìn)行長時序的動態(tài)監(jiān)測,提取積水邊界矢量并統(tǒng)計積水面積,分析積水范圍動態(tài)演化規(guī)律及其影響因素,以期為高潛水位采煤沉陷區(qū)土地利用規(guī)劃、土地復(fù)墾、水陸復(fù)合生態(tài)系統(tǒng)建立提供理論依據(jù)。

    1 研究區(qū)概況

    1031工作面位于安徽省濉溪縣五溝煤礦,是典型的高潛水位采煤區(qū),該地區(qū)降水量豐富,年均降水量827 mm,地下水位埋深1.5~4.0 m,多年平均埋深2.48 m[13]。1031工作面采動時間為2013年10月至2015年3月,工作面所在位置及現(xiàn)狀如圖1所示,底圖為2019年11月11日研究區(qū)遙感圖像,其中1031工作面全長1 411 m,寬180 m,煤層傾角為3°~16°,平均傾角8°,煤層平均厚度3.8 m,平均埋深(0)約364 m。工作面上覆巖層以砂巖和泥巖為主,厚度約94 m,松散層厚度約270 m,鉆孔柱狀如圖2所示。工作面所在區(qū)域地表地勢平坦,標(biāo)高+26.50~+27.56 m,無大流量河流經(jīng)過。由于該地區(qū)地下水埋深淺,降水量大[14],且地表下沉系數(shù)大,工作面采動后地表很快將形成積水[15]。

    圖1 五溝煤礦研究區(qū)所在位置及現(xiàn)狀

    2 積水范圍遙感監(jiān)測

    礦區(qū)地表積水周期長、范圍大,且所處位置交通不便,常規(guī)的觀測手段難以獲取大范圍、長時序、完整統(tǒng)一的積水范圍演化數(shù)據(jù);遙感觀測手段具有觀測周期短、數(shù)據(jù)量大、精度高等特點,因此,選擇遙感資料作為數(shù)據(jù)源,提取采煤沉陷區(qū)地表積水范圍信息[16]。本次選用Landsat-8系列衛(wèi)星的遙感數(shù)據(jù),影像空間分辨率為30 m,衛(wèi)星回歸周期為16 d。

    2.1 遙感圖像獲取與預(yù)處理

    五溝煤礦1031工作面采動時間為2013年10月至2015年3月,地表積水盆地不受其他工作面采動影響的截止時間為2016年10月,故選取2013—2016年Landsat-8獲取的遙感影像。剔除云量較大、能見度低的遙感影像,最終選取16期滿足水體提取基本要求的遙感影像。對選取的每幅影像進(jìn)行預(yù)處理,步驟包括:輻射定標(biāo)、大氣校正和圖像裁剪,得到最終影像[17]。

    圖2 五溝煤礦1031工作面鉆孔柱狀

    2.2 積水范圍提取

    目前,較高效的水體提取方法為水體指數(shù)法,常用的水體指數(shù)包括歸一化差異水體指數(shù)NDWI(Normalized Difference Water Index)[18]、改進(jìn)后的歸一化差異水體指數(shù)MNDWI[19](Modified NDWI)及自動水體提取指數(shù)AWEI(Automated Water Extrac-tion Index)[20]。NDWI利用近紅外波段和綠波段的反射率差值提取水體,該方法提取時可能有部分不透水層像素會被誤判為水體;故提出了MNDWI指數(shù),將NDWI中的近紅外波段換成短波紅外波段,增加水體的提取精度[19];AWEI針對Landsat-5影像設(shè)計,通過最大化水體和非水體區(qū)分度而得到最優(yōu)的系數(shù)組合[21]。結(jié)合研究區(qū)條件,對比各種水體指數(shù)特點,選擇MNDWI指數(shù)用于該研究區(qū)的水體提取。

    選擇Landsat-8影像中的綠波段(Band 3)和短波紅外波段(Band 6)反射值進(jìn)行計算,公式如下:

    式中:Green和MIR分別為綠光波段和短波紅外波段的反射值。

    通過MNDWI指數(shù)對遙感圖像進(jìn)行處理,根據(jù)影像內(nèi)的水體特征,統(tǒng)計積水面積并將積水邊界轉(zhuǎn)成矢量文件輸出。如圖3所示,底圖為Landsat-8標(biāo)準(zhǔn)彩色合成影像,白框部分為同期研究區(qū)積水矢量邊界,積水矢量邊界與水體邊界基本吻合,可知,通過MNDWI提取的積水邊界與水體邊界基本吻合,具有較高的準(zhǔn)確性,能夠滿足研究需求。

    圖3 2016年9月2日積水矢量邊界提取效果

    3 積水動態(tài)演化規(guī)律

    通過MNDWI方法,得到16幅積水范圍矢量并統(tǒng)計積水面積,結(jié)合工作面推進(jìn)位置,監(jiān)測積水范圍動態(tài)變化,分析地表積水范圍演化規(guī)律。

    圖4為1031工作面推進(jìn)位置與地表積水范圍變化疊加圖像,由圖4可以看出,地表積水形成明顯滯后于地下工作面采動時間,當(dāng)工作面推進(jìn)476 m后,即推進(jìn)距離接近采深(0)時,地表積水才開始形成;工作面推進(jìn)過程中,積水范圍演化方向與工作面推進(jìn)方向一致,隨著工作面推進(jìn),積水范圍沿工作面推進(jìn)方向擴大,積水面積日增長量約440 m2/d;2015年3月3日(工作面停止采動)至2015年7月30日,積水范圍仍然沿工作面推進(jìn)方向擴大,積水面積日增長量約120 m2/d;工作面停采后5個月左右,積水范圍不再擴大,趨于穩(wěn)定。

    為更方便研究積水范圍演化過程中的影響因素,通過工作面采動狀態(tài)及積水范圍的演化規(guī)律分析,將高潛水位采煤沉陷區(qū)積水范圍演化劃分為4個時期,如圖5所示。

    a. 未形成期 該階段為工作面開始采動至地表積水開始形成期間,地表積水尚未形成;

    圖4 2013—2016年工作面推進(jìn)距離與地表積水范圍變化

    圖5 五溝煤礦1031工作面積水面積變化

    b.同步增長期 該階段為地表首次觀測到積水至工作面停止采動期間,工作面處于采動狀態(tài),地表積水范圍隨工作面采動同步增長;

    c. 殘余增長期 該階段為工作面停止采動至停采后5個月左右,地表受采動殘余沉降影響仍呈動態(tài)變化,積水范圍沿工作面推進(jìn)方向呈增長趨勢,積水面積呈緩慢增長趨勢;

    d. 相對穩(wěn)定期 該階段為工作面停采5個月后,地表基本趨于穩(wěn)定,積水邊界不再沿采動方向擴大,積水面積整體趨于穩(wěn)定。

    4 地表積水演化的關(guān)鍵影響因素

    沉陷區(qū)積水范圍演化中,首先受工作面采動影響產(chǎn)生地表沉陷;另外,大氣降水和地下水補給是沉陷區(qū)積水的主要來源,占積水區(qū)水量來源90%以上,水分蒸發(fā)是主要的排泄方式[10],其他影響因素對積水范圍影響較小,非主要影響因素。

    為研究主要影響因素對積水范圍的作用,首先,根據(jù)2013—2016年國家氣象觀測數(shù)據(jù)計算研究區(qū)的降水量、蒸發(fā)量[22]及兩者差值,得到凈降水量,分析凈降水量對積水范圍演化各階段的影響;然后,根據(jù)濉溪縣地下水位埋深觀測數(shù)據(jù),分析地下水位埋深對積水演化各時期的影響;再次,結(jié)合工作面推進(jìn)距離,分析工作面推進(jìn)距離對積水范圍演化各時期的影響;最后,根據(jù)上述分析得到積水范圍演化的關(guān)鍵影響因素。

    4.1 凈降水量

    大氣降水、蒸發(fā)作用同步進(jìn)行,均對積水范圍產(chǎn)生影響,因此,應(yīng)綜合分析兩者對沉陷區(qū)積水范圍演化的作用。通過計算研究時期內(nèi)日降水量、日蒸發(fā)量差值之和,得到研究區(qū)凈降水量累計值。圖6為自首次觀測到積水以后的凈降水量與積水面積變化圖。

    圖6 凈降水量與積水面積相關(guān)關(guān)系

    積水未形成期,該沉陷區(qū)未產(chǎn)生積水。同步增長期和殘余增長期,積水面積整體呈增長趨勢,凈降水量與積水面積日增長量呈正相關(guān),隨著凈降水量增大,積水面積日增長量越大;凈降水量減小,積水面積日增長量相對減小。據(jù)圖6可知,2014年5月8日至2014年12月2日,凈降水量由–57 mm左右先減小至–170 mm左右,再增長到200 mm左右,同期積水面積由1.8 hm2增長到14.0 hm2,平均日增長量約584 m2;2014年12月2日至2015年3月3日(工作面停止采動)期間,凈降水量由200 mm減小至105 mm左右,同期積水面積由14.0 hm2增長到15.3 hm2,平均日增長量約116 m2,積水面積日增長量較前者明顯減小。相對穩(wěn)定期,凈降水量與積水面積日增長量負(fù)相關(guān),如2015年7月30日至2016年3月10日,凈降水量由180 mm減小至75 mm,由于凈降水量整體為正,對積水仍呈補充狀態(tài),同期積水面積呈小幅度增長,平均日增長量約為72 m2;2016年6月4日至2016年8月1日,凈降水量整體由15 mm增長至39 mm,但過程中凈降水量出現(xiàn)負(fù)值,同期積水面積呈小幅度減小,平均日減小量約170 m2。

    4.2 地下水位埋深

    淺層地下水在沉陷區(qū)形成積水后,與地表水體產(chǎn)生水力聯(lián)系,對地表水有一定的補充作用,是沉陷區(qū)積水主要來源之一[23],圖7為研究時期地下水位埋深與積水面積變化圖。

    圖7 地下水埋深與積水面積相關(guān)關(guān)系

    積水未形成期,地表未觀測到積水,該階段地表下沉量較小,未與地下水產(chǎn)生水力聯(lián)系;同步增長期與殘余增長期,地表下沉量變大,沉陷區(qū)形成積水區(qū)域,該階段地下水位埋深與積水面積日增加量呈負(fù)相關(guān),地下水位埋深快速減小,積水面積日增長量增大;地下水位埋深增加,地下水對沉陷區(qū)補給減少,積水面積日增長量減小。據(jù)圖7可知,2014年3月21日至2014年12月2日,地下水位埋深由3.7 m 左右減小到2.4 m左右,同期積水面積平均日增長量約504 m2;2015年3月24日至2015年7月30日,地下水位埋深由2.9 m左右減小至1.8 m,積水面積平均日增長量為139 m2;2014年12月2日至2015年3月24日,地下水位埋深由2.4 m左右增大至2.9 m,同期凈降水量由約210 mm減小到115 mm左右,積水面積平均日增長量約116 m2,較前兩區(qū)間增長速度明顯減??;相對穩(wěn)定期,沉陷區(qū)基本穩(wěn)定,地下水與沉陷區(qū)積水已產(chǎn)生穩(wěn)定水力聯(lián)系,淺層地下水為積水區(qū)域提供涵養(yǎng)條件[12],維持地表水的相對穩(wěn)定,該階段地下水位埋深與積水范圍演化相關(guān)性不明顯。

    4.3 工作面推進(jìn)距離

    工作面采動后,采區(qū)周圍巖體原始應(yīng)力平衡被破壞,在應(yīng)力達(dá)到重新平衡的過程中,上覆巖層產(chǎn)生連續(xù)的移動、變形和非連續(xù)破壞,地表產(chǎn)生移動盆地。當(dāng)移動盆地達(dá)到一定深度且與地下水產(chǎn)生水力聯(lián)系時,地表將形成積水區(qū)域,隨著地下工作面采動,地表積水區(qū)域隨之同步擴大。由積水范圍演化規(guī)律可知,地表積水范圍沿著工作面采動方向演化[24],可見工作面推進(jìn)距離與積水范圍演化具有明顯的相關(guān)性。

    圖8為工作面推進(jìn)距離與積水面積變化關(guān)系圖。據(jù)圖8可知,地表積水形成明顯滯后于工作面采動,在積水未形成期,受工作面采動影響巖層發(fā)生破壞,并逐漸傳遞至地表,地表形成下沉盆地,但下沉量較小,地表尚未形成穩(wěn)定積水區(qū)域;同步增長期,隨著下沉量增大,地表移動盆地已形成積水,截至工作面停止采動,積水面積隨工作面采動增長至15.3 hm2左右;殘余增長期,工作面停采后,采動導(dǎo)致的巖體移動尚未穩(wěn)定,該階段地表積水面積仍有小范圍增加,由15.3 hm2增長至17.1 hm2;相對穩(wěn)定期,該階段巖體應(yīng)力平衡基本穩(wěn)定,沉陷區(qū)地表趨于穩(wěn)定,但受降水量、蒸發(fā)量、地下水等因素影響,積水面積仍然存在小幅度波動。

    圖8 工作面推進(jìn)距離與地表積水面積相關(guān)關(guān)系

    4.4 關(guān)鍵影響因素確定

    綜上可知,積水未形成期,地表主要受工作面采動影響形成下沉盆地,由于下沉較小,地表尚未形成積水;同步增長期,地表下沉量增大,地表形成積水,該階段凈降水量、采動距離與積水面積日增長量呈正相關(guān),地下水埋深與積水面積日增長量呈負(fù)相關(guān);殘余增長期,該階段工作面停止采動,受殘余沉降影響,地表仍然有少量下沉,該階段積水面積仍然呈上升趨勢,凈降水量與積水面積日增長量呈正相關(guān),地下水埋深與積水面積日增長量呈負(fù)相關(guān);相對穩(wěn)定期,該階段地表下沉盆地基本穩(wěn)定,受降水量、蒸發(fā)量和地下水埋深的波動,積水區(qū)面積整體穩(wěn)定呈小范圍波動。

    凈降水量、地下水埋深、工作面推進(jìn)距離均與地表積水范圍演化有明顯相關(guān)性,但凈降水量、地下水位埋深受氣候條件影響均呈周期性變化[25],根據(jù)積水演化規(guī)律分析可知,工作面推進(jìn)距離增加為積水范圍演化的關(guān)鍵影響因素。

    5 工作面推進(jìn)距離與積水邊界角的關(guān)系

    為量化工作面采動距離與積水范圍演化的關(guān)系,提出積水邊界角的概念。工作面推進(jìn)距離是積水范圍演化的關(guān)鍵影響因素,與積水面積變化呈正相關(guān)。根據(jù)積水邊界角隨工作面采動距離的變化,建立工作面推進(jìn)方向上推進(jìn)位置與地表積水范圍之間的動態(tài)關(guān)系。

    如圖9所示,工作面推進(jìn)位置,工作面推進(jìn)方向剖面正上方積水邊界,連接,與工作面推進(jìn)方向的夾角定義為積水邊界角。圖中,在水平面投影距離為,2為工作面到積水水面的高差,0為采深,1為積水水面到原地表的距離,由于2遠(yuǎn)大于1,故2約等于0。

    圖9 積水邊界角示意圖

    積水邊界角計算公式為:

    式中:為積水邊界角,(°);0為采深,m;為在水平面的投影距離,m。

    根據(jù)式(2),計算每個觀測時間段積水邊界角。2013年10月至2014年3月21日地表未形成積水,積水邊界角記為0,從2014年3月21日開始統(tǒng)計,并利用統(tǒng)計數(shù)據(jù)繪制積水邊界角變化曲線,如圖10所示。

    圖10 積水邊界角變化曲線

    為便于研究積水演化周期內(nèi)積水邊界角的變化特征,根據(jù)工作面推進(jìn)速度與推進(jìn)距離,建立推進(jìn)距離與時間的關(guān)系,最終建立整個演化過程中積水邊界角與時間的關(guān)系。

    a. 未形成期 地表未形成積水,積水邊界角記為0。

    c. 殘余增長期 建立積水邊界角與工作面停采時間的關(guān)系。

    d. 相對穩(wěn)定期 該階段積水范圍整體動態(tài)平衡,積水邊界角小范圍波動,趨于90°。

    積水演化各階段中,積水邊界角與工作面采動時間的關(guān)系式為:

    由圖10、式(3)可得積水演化周期內(nèi)積水邊界角與工作面采動時間的關(guān)系。積水未形成期,不觀測積水邊界角;同步增長期和殘余增長期,積水邊界角整體呈減小趨勢,其中,同步增長期積水邊界角與時間呈三次函數(shù)關(guān)系,積水邊界角整體隨時間變化呈減小趨勢,即積水邊界角隨著工作面采動距離的增大而減??;殘余增長期,積水邊界角與采動時間呈一次函數(shù)關(guān)系,積水邊界角隨工作面停采時長的增加而減小;相對穩(wěn)定期,積水邊界角趨于穩(wěn)定,接近90°。

    6 結(jié)論

    a.高潛水位采煤沉陷區(qū)地表積水的形成分為4個階段:未形成期,開始采動至采動距離476 m;同步增長期,采動距離476 m至工作面停采;殘余增長期,工作面停采至停采后5個月左右;相對穩(wěn)定期,停采5個月后。

    b.積水范圍演化受到大氣降水、蒸發(fā)、地下水位埋深、工作面采動等因素的影響,在同步增長和殘余增長期,凈降水量、工作面采動距離與積水面積日增長量呈正相關(guān),地下水位埋深與積水面積日增長量呈負(fù)相關(guān);工作面推進(jìn)距離為積水范圍演化的關(guān)鍵影響因素。

    c.積水邊界角在同步增長期和殘余增長期整體呈減小趨勢;同步增長期的積水邊界角隨工作面推進(jìn)距離的增加而減??;殘余增長期的積水邊界角隨工作面停采時長的增加而減??;積水穩(wěn)定期,積水邊界角趨于穩(wěn)定,接近90°。

    請聽作者語音介紹創(chuàng)新技術(shù)成果等信息,歡迎與作者進(jìn)行交流

    [1] 英國石油公司.《BP世界能源統(tǒng)計年鑒》[EB/OL]. (2019- 06-14)[2019-10-10] https://www.bp.com/zh_cn/china/home/news/ reports.html. British Petroleum Company.《BP Statistical Review of World Energy》[EB/OL]. (2019-06-14)[2019-10-10] https://www.bp. com/zh_cn/china/home/news/reports.html.

    [2] 劉輝,鄧喀中,雷少剛,等. 采動地裂縫動態(tài)發(fā)育規(guī)律及治理標(biāo)準(zhǔn)探討[J]. 采礦與安全工程學(xué)報,2017,34(5):884–890. LIU Hui,DENG Kazhong,LEI Shaogang,et al. Dynamic developing law and governance standard of ground fissures caused by underground mining[J]. Journal of Mining & Safety Engineering,2017,34(5):884–890.

    [3] LIU Hui,DENG Kazhong,ZHU Xiaojun,et al. Effects of mining speed on the developmental features of mining-induced ground fissures[J]. Bulletin of Engineering Geology and the Environment,2019,78(8):6297–6309.

    [4] ZHU Xiaojun,GUO Guangli,ZHA Jianfeng,et al. Surface dynamic subsidence prediction model of solid backfill mining[J]. Environmental Earth Sciences,2016,75(6):1–9.

    [5] 孔令健,姜春露,鄭劉根,等. 淮北臨渙礦采煤沉陷區(qū)不同水體水化學(xué)特征及其影響因素[J]. 湖泊科學(xué),2017,29(5):1158–1167. KONG Lingjian,JIANG Chunlu,ZHENG Liugen,et al. Characters of hydrochemistry and their influenced factors of different waters in the Linhuan coal mining subsidence area of Huaibei City[J]. Journal of Lake Sciences,2017,29(5):1158–1167.

    [6] 彭蘇萍,王磊,孟召平,等. 遙感技術(shù)在煤礦區(qū)積水塌陷動態(tài)監(jiān)測中的應(yīng)用:以淮南礦區(qū)為例[J]. 煤炭學(xué)報,2002,27(4):374–378. PENG Suping,WANG Lei,MENG Zhaoping,et al. Monitoring the seeper subside in coal district by the remote sensing:Examples from Huainan coal district[J]. Journal of China Coal Society,2002,27(4):374–378.

    [7] 魏礦靈,王啟春,郭廣禮,等. 利用遙感數(shù)據(jù)監(jiān)測礦區(qū)開采沉陷積水變化[J]. 煤礦安全,2014,45(1):13–16. WEI Kuangling,WANG Qichun,GUO Guangli,et al. Using remote sensing data to monitor the change of mining subsidence water accumulation[J]. Safety in Coal Mines,2014,45(1):13–16.

    [8] 馬建威,黃詩峰,許宗男. 基于遙感的1973—2015年武漢市湖泊水域面積動態(tài)監(jiān)測與分析研究[J]. 水利學(xué)報,2017,48(8):903–913. MA Jianwei,HUANG Shifeng,XU Zongnan. Satellite remote sensing of lake area in Wuhan from 1973 to 2015[J]. Journal of Hydraulic Engineering,2017,48(8):903–913.

    [9] 陸垂裕,秦大庸,張俊娥,等. 面向?qū)ο竽K化的分布式水文模型MODCYCLEⅠ:模型原理與開發(fā)篇[J]. 水利學(xué)報,2012,43(10):1135–1145. LU Chuiyu,QIN Dayong,ZHANG Jun’e,et al. MODCYCLE-an object oriented modularized hydrological model Ⅰ:Theory and development[J]. Journal of Hydraulic Engineering,2012,43(10):1135–1145.

    [10] 陸垂裕,陸春輝,李慧,等. 淮南采煤沉陷區(qū)積水過程地下水作用機制[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(10):122–131. LU Chuiyu,LU Chunhui,LI Hui,et al. Groundwater influencing mechanism in water-ponding process of Huainan coal mining subsidence area[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2015,31(10):122–131.

    [11] 李慧,陸垂裕,孫青言,等. 分布式“河道–沉陷區(qū)–地下水”水循環(huán)耦合模型—Ⅰ模型原理與開發(fā)[J]. 水科學(xué)進(jìn)展,2016,27(2):214–223.LI Hui,LU Chuiyu,SUN Qingyan,et al. Research on the distributed “river-subsidence area-groundwater” coupling simulation of coal mining subsidence areaⅠ:Theory and development[J]. Advances in Water Science,2016,27(2):214–223.

    [12] 李慧. 采煤沉陷區(qū)分布式水循環(huán)模型研究[D]. 北京:中國水利水電科學(xué)研究院,2016. LI Hui. Research on the distributed water cycle model of coal mining subsidence area[D]. Beijing:China Institute of Water Resources and Hydropower Research,2016.

    [13] 王發(fā)信,柏菊. 淮北平原淺層地下水埋深區(qū)域分布特點[J]. 地下水,2014,36(5):51–53. WANG Faxin,BO Ju. Distribution characteristics of shallow groundwater depth in Huaibei plain[J]. Ground Water,2014,36(5):51–53.

    [14] 喬叢林,史明禮,蘇婭,等. 淮北平原地區(qū)水文特征[J]. 水文,2000,20(3):55–58. QIAO Conglin,SHI Mingli,SU Ya,et al. Hydrological features of Huaibei plain[J]. Journal of China Hydrology,2000,20(3):55–58.

    [15] 徐良驥,嚴(yán)家平,高永梅. 安徽省兩淮礦區(qū)地質(zhì)環(huán)境治理技術(shù)[J]. 煤田地質(zhì)與勘探,2007,35(6):37–40. XU Liangji,YAN Jiaping,GAO Yongmei. Comprehensive treatment technique of geological environment in Huainan & Huaibei mining area[J]. Coal Geology & Exploration,2007,35(6):37–40.

    [16] 王海波,馬明國. 基于遙感的湖泊水域動態(tài)變化監(jiān)測研究進(jìn)展[J]. 遙感技術(shù)與應(yīng)用,2009,24(5):674–684. WANG Haibo,MA Mingguo. A review of monitoring change in lake water areas based on remote sensing[J]. Remote Sensing Technology and Application,2009,24(5):674–684.

    [17] 祝佳. Landsat8衛(wèi)星遙感數(shù)據(jù)預(yù)處理方法[J]. 國土資源遙感,2016,28(2):21–27. ZHU Jia. Analysis of Landsat8 satellite remote sensing data preprocessing[J]. Remote Sensing for Land & Resources,2016,28(2):21–27.

    [18] MCFEETERS S K. The use of the Normalized Difference Water Index(NDWI) in the delineation of open water features[J]. International Journal of Remote Sensing,1996,17(7):1425–1432.

    [19] XU Hanqiu. Modification of Normalised Difference Water Index(NDWI) to enhance open water features in remotely sensed imagery[J]. International Journal of Remote Sensing,2006,27(14):3025–3033.

    [20] FEYISA G L,MEILBY H,F(xiàn)ENSHOLT R,et al. Automated water extraction index:A new technique for surface water mapping using Landsat imagery[J]. Remote Sensing of Environment,2014,140:23–35.

    [21] 計璐艷,尹丹艷,宮鵬. Landsat長時間序列的陽澄湖湖面圍網(wǎng)時空變化[J]. 遙感學(xué)報,2019,23(4):717–729.JI Luyan,YIN Danyan,GONG Peng. Temporal-spatial study on enclosure culture area in Yangcheng lake with long-term Landsat time series[J]. Journal of Remote Sensing,2019,23(4):717–729.

    [22] 張晶,李妍清. 基于地統(tǒng)計學(xué)的柬埔寨年降水量數(shù)據(jù)空間插值[J]. 人民長江,2018,49(22):100–103. ZHANG Jing,LI Yanqing. Spatial interpolation of annual precipitation data in Cambodia based on geostatistics[J]. Yangtze River,2018,49(22):100–103.

    [23] 張磊,秦小光,劉嘉麒,等. 淮南采煤沉陷區(qū)積水來源的氫氧穩(wěn)定同位素證據(jù)[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版),2015,45(5):1502–1514. ZHANG Lei,QIN Xiaoguang,LIU Jiaqi,et al. Characters of hydrogen and oxygen stable isotope of different water bodies in Huainan coal mining area[J]. Journal of Jilin University(Earth Science Edition),2015,45(5):1502–1514.

    [24] 何國清,楊倫,凌賡娣,等. 礦山開采沉陷學(xué)[M]. 徐州:中國礦業(yè)大學(xué)出版社,1994:84–92. HE Guoqing,YANG Lun,LING Gengti,et al. Mining subsidence[M]. Xuzhou:China University of Mining and Technology Press,1994:84–92.

    [25] 陸小明,楊秒,王振龍,等. 淮北平原氣象因素對裸地潛水蒸發(fā)的影響[J]. 灌溉排水學(xué)報,2019,38(4):84–91. LU Xiaoming,YANG Miao,WANG Zhenlong,et al. Metrological factors affecting evaporation of shallow groundwater in the absence of plants in Huaibei plain[J]. Journal of Irrigation and Drainage,2019,38(4):84–91.

    Dynamic evolution law of water accumulation range in coal mining subsidence area with high-level groundwater

    CHEN Xiaoxie1,2, ZHANG Wentao3, ZHU Xiaojun1,2, LIU Hui1,2,4, ZHANG Min1,2, SUN Aiguo5, YANG Xiaoyu6

    (1. College of Resources and Environmental Engineering, Anhui University, Hefei 230601, China; 2. Anhui Province Engineering Laboratory for Mine Ecological Remediation, Hefei 230601, China; 3. School of Earth and Environment, Anhui University of Science & Technology, Huainan 232001, China; 4. School of Mining and Geomatics, Hebei University of Engineering, Handan 056038, China; 5. Wugou Coal Mine of Anhui Hengyuan Coal Power Co. Ltd.,Huaibei 235131, China; 6. Hefei College of Finance & Economics, Hefei 230601, China)

    In order to study the dynamic evolution law of the water accumulation range with high-level groundwater, the working face 1031 of Wugou coal mine was taken as the research area.Based on Landsat-8 remote sensing data,the evolution law of water accumulation range above working face 1031 during 2013-2017 was summarized and the influence of factors such as net rainfall, groundwater depth and working surface propulsion distance on the evolution of water accumulation range was analyzed.The concept of the boundary angle of the water accumulation range was proposed, and the functional relationship between the angle and the mining time was established.The study found that the evolution of surface water in the coal mining subsidence area was divided into four periods: unformed, synchronous growth, residual growth, stable. The groundwater depth is negatively correlated with the daily growth of the water accumulation range. The working face advancement distance is the driving factor for the change of the water accumulation range. When the working distance is about 476 m, water is accumulated on the surface, about 5 months after stopping mining, the area of water accumulation range tends to be stable, in the synchronous growth and residual growth period, the angle generally decreases, in the stable period the angle tends to be 90°as a whole, but it fluctuates due to other factors.This study provides a theoretical basis for land use planning, land reclamation, and establishment of a land-water complex ecosystem in the coal mining subsidence area with high-level groundwater.

    high-level groundwater; coal mining subsidence area; remote sensing monitoring; evolution of subsidence water accumulation range; the boundary angle of water accumulation range; Wugou coal mine

    TD803

    A

    10.3969/j.issn.1001-1986.2020.02.020

    1001-1986(2020)02-0126-08

    2019-10-18;

    2020-02-22

    國家自然科學(xué)基金項目(51804001);安徽省自然科學(xué)基金項目(1808085QE147);安徽省高校自然科學(xué)研究重點項目(KJ2018A0003)

    National Natural Science Foundation of China(51804001);Natural Science Foundation of Anhui Province, China(1808085QE147);Key Project of Natural Science Research of Universities in Anhui Province(KJ2018A0003)

    陳曉謝,1993年生,男,安徽合肥人,碩士研究生,從事礦山開采沉陷方向研究. E-mail:chenxiaoxie99@163.com

    劉輝,1982年生,男,山東肥城人,博士,副教授,從事礦山開采沉陷及生態(tài)環(huán)境治理等方面的教學(xué)與科研工作. E-mail:lhui99@aliyun.com

    陳曉謝,張文濤,朱曉峻,等. 高潛水位采煤沉陷區(qū)積水范圍動態(tài)演化規(guī)律[J]. 煤田地質(zhì)與勘探,2020,48(2):126–133.

    CHEN Xiaoxie,ZHANG Wentao,ZHU Xiaojun,et al. Dynamic evolution law of water accumulation range in coal mining subsidence area with high-level groundwater[J]. Coal Geology & Exploration,2020,48(2):126–133.

    (責(zé)任編輯 周建軍)

    猜你喜歡
    增長量采動積水
    原來是輸卵管積水惹的禍
    不同配比的基質(zhì)對不同多肉植物生長的影響
    小熊當(dāng)當(dāng)玩積水
    不同坡度和不同坡向?qū)鸹ú枭L量的影響
    采動影響下淺埋輸氣管道與土體耦合作用機理
    原來是輸卵管積水惹的禍
    金都1號火龍果果實發(fā)育過程中品質(zhì)的變化規(guī)律分析
    施用新型肥料對閩南桉樹生長的影響
    深部采動巷道頂板穩(wěn)定性分析與控制
    采動巖體滲流力學(xué)研究進(jìn)展
    免费不卡黄色视频| 日本91视频免费播放| 亚洲第一av免费看| 国产有黄有色有爽视频| 久久久亚洲精品成人影院| 精品国产露脸久久av麻豆| 日韩一区二区视频免费看| 国产精品99久久99久久久不卡 | 亚洲av电影在线观看一区二区三区| 欧美黑人欧美精品刺激| 91成人精品电影| 国产精品一区二区精品视频观看| 国产激情久久老熟女| 亚洲精品自拍成人| 老司机影院毛片| 成年动漫av网址| 日本欧美国产在线视频| av在线播放精品| 中国国产av一级| 亚洲一级一片aⅴ在线观看| 亚洲精品视频女| 一区二区日韩欧美中文字幕| 国产成人欧美| 999精品在线视频| 男人爽女人下面视频在线观看| 男女边吃奶边做爰视频| 久久 成人 亚洲| av国产久精品久网站免费入址| 高清不卡的av网站| 免费不卡黄色视频| 久久99精品国语久久久| 久久天堂一区二区三区四区| 国产成人欧美| av卡一久久| 性少妇av在线| 国产精品成人在线| 国产精品.久久久| 巨乳人妻的诱惑在线观看| 嫩草影视91久久| 丰满饥渴人妻一区二区三| 成年人午夜在线观看视频| 人妻人人澡人人爽人人| 中文字幕高清在线视频| 天堂俺去俺来也www色官网| 宅男免费午夜| 欧美在线一区亚洲| 少妇的丰满在线观看| 一本久久精品| 三上悠亚av全集在线观看| 久久久精品免费免费高清| 亚洲av日韩精品久久久久久密 | 日本91视频免费播放| 黄片小视频在线播放| 色网站视频免费| 欧美亚洲 丝袜 人妻 在线| 亚洲av男天堂| 国产一区二区 视频在线| 我要看黄色一级片免费的| 婷婷色av中文字幕| 国产午夜精品一二区理论片| 国产精品一国产av| 51午夜福利影视在线观看| 岛国毛片在线播放| 日韩,欧美,国产一区二区三区| 亚洲伊人色综图| 久久久久网色| 丰满乱子伦码专区| 国产成人91sexporn| 母亲3免费完整高清在线观看| 国产av精品麻豆| 久久鲁丝午夜福利片| 亚洲人成网站在线观看播放| 男人添女人高潮全过程视频| 纯流量卡能插随身wifi吗| 美女国产高潮福利片在线看| 久久韩国三级中文字幕| 91老司机精品| 久久99热这里只频精品6学生| 日本爱情动作片www.在线观看| 热99国产精品久久久久久7| 欧美 日韩 精品 国产| 午夜免费观看性视频| 精品福利永久在线观看| 最近最新中文字幕免费大全7| 欧美日韩视频高清一区二区三区二| av不卡在线播放| 97在线人人人人妻| 最近中文字幕2019免费版| 欧美人与性动交α欧美精品济南到| 日韩伦理黄色片| 黄色毛片三级朝国网站| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻丝袜一区二区 | 欧美精品av麻豆av| 午夜激情av网站| 在线观看人妻少妇| 国产精品一区二区在线观看99| 亚洲精品,欧美精品| 一边亲一边摸免费视频| av不卡在线播放| 韩国精品一区二区三区| 色94色欧美一区二区| av在线观看视频网站免费| 国产男人的电影天堂91| 热99国产精品久久久久久7| 国产精品久久久久久精品电影小说| av在线app专区| 亚洲国产精品999| 天堂中文最新版在线下载| 天天躁夜夜躁狠狠久久av| 午夜日本视频在线| 国产av码专区亚洲av| 欧美最新免费一区二区三区| 国产欧美亚洲国产| 久久毛片免费看一区二区三区| 波多野结衣av一区二区av| 欧美成人精品欧美一级黄| 免费看av在线观看网站| 国产一区二区 视频在线| 大片电影免费在线观看免费| 日本一区二区免费在线视频| 捣出白浆h1v1| 免费黄网站久久成人精品| 丁香六月天网| 精品一区二区三卡| 国产精品久久久人人做人人爽| 99re6热这里在线精品视频| 亚洲欧美成人精品一区二区| 午夜福利一区二区在线看| 中文字幕亚洲精品专区| 久久久国产欧美日韩av| 亚洲第一av免费看| 男女午夜视频在线观看| av电影中文网址| 两个人看的免费小视频| 99热网站在线观看| 爱豆传媒免费全集在线观看| 久久久久国产一级毛片高清牌| 久久ye,这里只有精品| av国产精品久久久久影院| 黄色视频在线播放观看不卡| 久久鲁丝午夜福利片| 日韩一本色道免费dvd| 又大又爽又粗| a 毛片基地| 国产精品久久久久久精品电影小说| 好男人视频免费观看在线| 国产一区二区 视频在线| 国产片内射在线| 国产精品欧美亚洲77777| 亚洲精品国产一区二区精华液| 高清av免费在线| 亚洲免费av在线视频| 婷婷色综合大香蕉| 国产精品无大码| 日韩制服骚丝袜av| 成年女人毛片免费观看观看9 | 免费看不卡的av| 男的添女的下面高潮视频| 黄色视频不卡| e午夜精品久久久久久久| 一区二区三区四区激情视频| www.自偷自拍.com| 99精国产麻豆久久婷婷| 精品一区二区免费观看| 热99国产精品久久久久久7| 国产精品久久久久久精品电影小说| 亚洲欧美成人精品一区二区| 国产精品久久久av美女十八| 精品一区二区三卡| 国产色婷婷99| 精品亚洲乱码少妇综合久久| 97人妻天天添夜夜摸| 精品国产国语对白av| 久久精品国产综合久久久| 欧美日韩精品网址| 黄色 视频免费看| 一级爰片在线观看| 19禁男女啪啪无遮挡网站| 日本色播在线视频| 又粗又硬又长又爽又黄的视频| 9191精品国产免费久久| 日日啪夜夜爽| 99国产综合亚洲精品| 99久久精品国产亚洲精品| 美女福利国产在线| 亚洲精品国产av蜜桃| 在线精品无人区一区二区三| 久久精品久久久久久噜噜老黄| 午夜日韩欧美国产| 欧美精品高潮呻吟av久久| 在线天堂最新版资源| av片东京热男人的天堂| 成年人午夜在线观看视频| 在线观看人妻少妇| 两性夫妻黄色片| 久久精品久久精品一区二区三区| 自线自在国产av| 日韩欧美一区视频在线观看| 伊人亚洲综合成人网| 成人三级做爰电影| 少妇 在线观看| 国产成人午夜福利电影在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦视频在线资源免费观看| 亚洲欧美精品综合一区二区三区| 欧美日韩亚洲综合一区二区三区_| 七月丁香在线播放| 我要看黄色一级片免费的| 午夜影院在线不卡| 在线观看免费日韩欧美大片| 国产成人免费无遮挡视频| 另类精品久久| 九草在线视频观看| 国产亚洲午夜精品一区二区久久| 国产亚洲av高清不卡| 一本一本久久a久久精品综合妖精| 青草久久国产| 免费观看性生交大片5| 女人被躁到高潮嗷嗷叫费观| 90打野战视频偷拍视频| 久久热在线av| 亚洲三区欧美一区| 另类精品久久| 国产精品一区二区在线不卡| 午夜福利视频精品| 久久久久精品国产欧美久久久 | 久久99热这里只频精品6学生| 亚洲五月色婷婷综合| 欧美成人精品欧美一级黄| 久久狼人影院| av有码第一页| 老司机影院毛片| 亚洲熟女精品中文字幕| 一级毛片电影观看| 一区二区三区乱码不卡18| 国产精品一二三区在线看| 午夜激情av网站| 十八禁高潮呻吟视频| 深夜精品福利| 一级,二级,三级黄色视频| 免费观看a级毛片全部| 午夜久久久在线观看| 丁香六月欧美| 九色亚洲精品在线播放| 在线观看人妻少妇| 99国产精品免费福利视频| 天天操日日干夜夜撸| 制服诱惑二区| 搡老乐熟女国产| 国产在视频线精品| 亚洲欧美精品综合一区二区三区| 亚洲国产毛片av蜜桃av| 午夜福利,免费看| 视频在线观看一区二区三区| 日本猛色少妇xxxxx猛交久久| 国产黄色视频一区二区在线观看| 丝袜喷水一区| 夫妻性生交免费视频一级片| 别揉我奶头~嗯~啊~动态视频 | 国产精品久久久久成人av| 中文字幕亚洲精品专区| 老司机影院毛片| 亚洲精品,欧美精品| 夫妻午夜视频| av不卡在线播放| 精品人妻一区二区三区麻豆| 叶爱在线成人免费视频播放| 日韩制服骚丝袜av| av国产精品久久久久影院| 丝袜脚勾引网站| 99久久人妻综合| 赤兔流量卡办理| 水蜜桃什么品种好| 久久 成人 亚洲| 中文字幕制服av| 午夜免费观看性视频| 街头女战士在线观看网站| 美女福利国产在线| 天天影视国产精品| 日韩一本色道免费dvd| 欧美国产精品一级二级三级| 我要看黄色一级片免费的| 成年动漫av网址| 青草久久国产| 丝袜美足系列| 亚洲欧美色中文字幕在线| 18在线观看网站| 久久ye,这里只有精品| 三上悠亚av全集在线观看| 成人亚洲精品一区在线观看| 亚洲一级一片aⅴ在线观看| 丝袜美足系列| 亚洲三区欧美一区| 国产成人系列免费观看| 免费女性裸体啪啪无遮挡网站| 满18在线观看网站| 精品一区二区三区四区五区乱码 | 亚洲av成人精品一二三区| 亚洲av成人精品一二三区| 中国三级夫妇交换| 国产精品免费大片| 一本久久精品| 久久久久久久国产电影| 久久av网站| 中文天堂在线官网| 成人国语在线视频| 91精品伊人久久大香线蕉| 国产精品一二三区在线看| 欧美日本中文国产一区发布| 国产男女内射视频| 又黄又粗又硬又大视频| 久久人人97超碰香蕉20202| 日韩欧美一区视频在线观看| www.av在线官网国产| 色网站视频免费| 男女之事视频高清在线观看 | 欧美黄色片欧美黄色片| 考比视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 免费人妻精品一区二区三区视频| √禁漫天堂资源中文www| 一级毛片我不卡| 亚洲情色 制服丝袜| 老司机在亚洲福利影院| 国产精品久久久久久久久免| 中文字幕亚洲精品专区| 国产淫语在线视频| 99国产综合亚洲精品| 男的添女的下面高潮视频| 亚洲国产av新网站| 色综合欧美亚洲国产小说| √禁漫天堂资源中文www| 精品久久久精品久久久| 免费在线观看视频国产中文字幕亚洲 | 在线观看免费高清a一片| 老司机亚洲免费影院| 自线自在国产av| 亚洲中文av在线| 久久久久精品久久久久真实原创| 午夜日韩欧美国产| 男人添女人高潮全过程视频| 99精品久久久久人妻精品| 亚洲av中文av极速乱| 国产片内射在线| 亚洲美女搞黄在线观看| 51午夜福利影视在线观看| 看非洲黑人一级黄片| 欧美人与性动交α欧美软件| 欧美黄色片欧美黄色片| 日韩中文字幕视频在线看片| 日韩中文字幕欧美一区二区 | 亚洲第一av免费看| 不卡视频在线观看欧美| 母亲3免费完整高清在线观看| 国产精品人妻久久久影院| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美一区二区三区黑人| 欧美日韩视频精品一区| 成人国语在线视频| 毛片一级片免费看久久久久| 黄色视频在线播放观看不卡| 嫩草影视91久久| 欧美日本中文国产一区发布| 国产精品.久久久| 老鸭窝网址在线观看| 丰满乱子伦码专区| 最近手机中文字幕大全| 精品人妻一区二区三区麻豆| 亚洲成色77777| 99精品久久久久人妻精品| 久久青草综合色| 免费看av在线观看网站| 亚洲在久久综合| 国产男女超爽视频在线观看| 人人妻人人澡人人爽人人夜夜| 99热国产这里只有精品6| 咕卡用的链子| 午夜影院在线不卡| 成人毛片60女人毛片免费| 亚洲综合色网址| 色播在线永久视频| 国产亚洲精品第一综合不卡| 人体艺术视频欧美日本| 一区在线观看完整版| 精品国产国语对白av| 久久99热这里只频精品6学生| 悠悠久久av| 亚洲一区二区三区欧美精品| 国产老妇伦熟女老妇高清| 人人妻人人澡人人爽人人夜夜| 国产精品一区二区在线不卡| 无遮挡黄片免费观看| 美女主播在线视频| 九色亚洲精品在线播放| 成人亚洲欧美一区二区av| 青青草视频在线视频观看| 美女中出高潮动态图| 少妇人妻久久综合中文| 亚洲成人一二三区av| 国产黄色免费在线视频| 美女视频免费永久观看网站| 国产精品秋霞免费鲁丝片| 熟妇人妻不卡中文字幕| 操美女的视频在线观看| 青春草视频在线免费观看| 国产亚洲午夜精品一区二区久久| 免费在线观看完整版高清| 黄片播放在线免费| 久久久久久久久久久免费av| 国产精品久久久久成人av| 永久免费av网站大全| 亚洲欧美精品综合一区二区三区| 国产精品无大码| 久久久久久久国产电影| 亚洲一级一片aⅴ在线观看| 国产免费又黄又爽又色| av天堂久久9| 精品少妇黑人巨大在线播放| 亚洲国产欧美网| 免费女性裸体啪啪无遮挡网站| 男人爽女人下面视频在线观看| 久久免费观看电影| 亚洲精品中文字幕在线视频| 极品人妻少妇av视频| 国产成人系列免费观看| 伦理电影免费视频| 黄色视频不卡| 高清欧美精品videossex| 人成视频在线观看免费观看| 国产亚洲av高清不卡| 成人黄色视频免费在线看| 熟妇人妻不卡中文字幕| 久久久久国产一级毛片高清牌| 一区二区日韩欧美中文字幕| 日韩成人av中文字幕在线观看| 久久久久久久久久久免费av| 在现免费观看毛片| 超色免费av| 国产欧美日韩一区二区三区在线| 日日撸夜夜添| 精品一品国产午夜福利视频| 色视频在线一区二区三区| 国产免费现黄频在线看| 老司机深夜福利视频在线观看 | 一个人免费看片子| 国产激情久久老熟女| 亚洲,欧美精品.| 中文字幕av电影在线播放| 国产av码专区亚洲av| 亚洲欧美中文字幕日韩二区| 大香蕉久久成人网| 亚洲av在线观看美女高潮| 少妇的丰满在线观看| av视频免费观看在线观看| 日本av手机在线免费观看| 精品酒店卫生间| 又黄又粗又硬又大视频| 国产黄色免费在线视频| 久久人人97超碰香蕉20202| 中文字幕色久视频| 黄片小视频在线播放| 午夜福利视频精品| 永久免费av网站大全| 性色av一级| 精品一区二区三区四区五区乱码 | 精品人妻熟女毛片av久久网站| 99久久99久久久精品蜜桃| 亚洲欧美一区二区三区黑人| 韩国高清视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 日韩精品免费视频一区二区三区| 亚洲欧洲日产国产| 久久鲁丝午夜福利片| 女人爽到高潮嗷嗷叫在线视频| 大话2 男鬼变身卡| 国产无遮挡羞羞视频在线观看| 又大又黄又爽视频免费| 丝袜喷水一区| 成人漫画全彩无遮挡| 午夜福利免费观看在线| 欧美日韩精品网址| 亚洲熟女精品中文字幕| 久久久国产欧美日韩av| 久久影院123| 女性生殖器流出的白浆| 久久久久精品性色| 观看av在线不卡| 亚洲国产毛片av蜜桃av| 9色porny在线观看| 精品国产露脸久久av麻豆| 少妇 在线观看| 99热国产这里只有精品6| 日本欧美视频一区| 我的亚洲天堂| 欧美激情高清一区二区三区 | 亚洲国产av影院在线观看| 日韩精品免费视频一区二区三区| 天美传媒精品一区二区| 男女之事视频高清在线观看 | 男女午夜视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 麻豆av在线久日| 午夜精品国产一区二区电影| 2018国产大陆天天弄谢| 啦啦啦 在线观看视频| 日本欧美国产在线视频| 欧美日韩视频精品一区| 免费少妇av软件| 欧美少妇被猛烈插入视频| 精品酒店卫生间| 亚洲国产精品成人久久小说| 在线观看免费午夜福利视频| 成年动漫av网址| 国精品久久久久久国模美| 高清黄色对白视频在线免费看| 国产精品国产三级国产专区5o| 国产国语露脸激情在线看| 成年动漫av网址| 丁香六月欧美| 亚洲七黄色美女视频| 亚洲精品中文字幕在线视频| xxx大片免费视频| 十八禁网站网址无遮挡| 女性被躁到高潮视频| 免费观看a级毛片全部| 巨乳人妻的诱惑在线观看| 精品少妇久久久久久888优播| 国产伦人伦偷精品视频| 999精品在线视频| 精品午夜福利在线看| 亚洲精品美女久久久久99蜜臀 | 美女午夜性视频免费| 久久人人爽av亚洲精品天堂| 免费女性裸体啪啪无遮挡网站| 性色av一级| 狂野欧美激情性xxxx| 亚洲在久久综合| 国产精品嫩草影院av在线观看| 哪个播放器可以免费观看大片| 国产xxxxx性猛交| 国产99久久九九免费精品| 欧美xxⅹ黑人| 国产精品久久久久久精品古装| 精品午夜福利在线看| 亚洲精品美女久久av网站| 美女国产高潮福利片在线看| 久久久久精品久久久久真实原创| 80岁老熟妇乱子伦牲交| 午夜91福利影院| 一区二区三区精品91| 亚洲少妇的诱惑av| 久久综合国产亚洲精品| 成年美女黄网站色视频大全免费| 飞空精品影院首页| 男女边吃奶边做爰视频| 国产男女超爽视频在线观看| 国产精品免费视频内射| 亚洲国产av影院在线观看| 亚洲成人国产一区在线观看 | 精品国产一区二区三区四区第35| 精品少妇久久久久久888优播| 国产成人精品久久二区二区91 | 欧美精品一区二区免费开放| 久久国产精品男人的天堂亚洲| 校园人妻丝袜中文字幕| 欧美精品一区二区大全| www.自偷自拍.com| 久久影院123| 欧美xxⅹ黑人| 男人舔女人的私密视频| 国产精品一区二区在线不卡| 男女之事视频高清在线观看 | 精品一区二区三卡| 久久韩国三级中文字幕| 亚洲人成电影观看| 亚洲第一av免费看| 久久韩国三级中文字幕| 日韩大码丰满熟妇| 无遮挡黄片免费观看| 亚洲 欧美一区二区三区| 久久免费观看电影| 在线观看人妻少妇| 一区二区三区四区激情视频| 亚洲精品视频女| av在线老鸭窝| 丁香六月欧美| 18在线观看网站| 天天影视国产精品| 欧美激情高清一区二区三区 | 狂野欧美激情性xxxx| 九草在线视频观看| 99热全是精品| 亚洲综合色网址| 亚洲av电影在线进入| 热re99久久精品国产66热6| 高清av免费在线| 天天添夜夜摸| 国产成人欧美| 丰满乱子伦码专区| 狠狠婷婷综合久久久久久88av| 韩国高清视频一区二区三区| 大码成人一级视频| 人人妻人人澡人人爽人人夜夜| 久久久精品区二区三区| 国产精品免费大片| 人体艺术视频欧美日本| 晚上一个人看的免费电影| 亚洲欧洲精品一区二区精品久久久 | 日本vs欧美在线观看视频| 国产av一区二区精品久久| 高清欧美精品videossex| 久久精品久久久久久噜噜老黄| 熟女av电影| 久久久久精品性色|