• <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)展
    在线观看舔阴道视频| 视频区图区小说| 国产精品免费一区二区三区在线 | tocl精华| 国产精品.久久久| h视频一区二区三区| 无限看片的www在线观看| av网站在线播放免费| 三级毛片av免费| 亚洲七黄色美女视频| 亚洲久久久国产精品| 在线观看人妻少妇| 国产精品免费视频内射| 国产aⅴ精品一区二区三区波| 国产精品美女特级片免费视频播放器 | 巨乳人妻的诱惑在线观看| 日本黄色视频三级网站网址 | 99香蕉大伊视频| 视频在线观看一区二区三区| 视频区欧美日本亚洲| 99精品欧美一区二区三区四区| 人人妻,人人澡人人爽秒播| 侵犯人妻中文字幕一二三四区| 美女福利国产在线| 成人18禁高潮啪啪吃奶动态图| 免费久久久久久久精品成人欧美视频| 国产色视频综合| 亚洲免费av在线视频| av又黄又爽大尺度在线免费看| 成年人黄色毛片网站| 亚洲七黄色美女视频| 黄色视频不卡| 国产精品免费一区二区三区在线 | 女警被强在线播放| 色精品久久人妻99蜜桃| 一个人免费在线观看的高清视频| 99在线人妻在线中文字幕 | 妹子高潮喷水视频| 精品福利观看| 丝袜美足系列| 天堂8中文在线网| 国产激情久久老熟女| 中文字幕色久视频| 多毛熟女@视频| 少妇的丰满在线观看| 国产国语露脸激情在线看| a级毛片在线看网站| 久久人妻熟女aⅴ| 中文字幕另类日韩欧美亚洲嫩草| 露出奶头的视频| 亚洲成人国产一区在线观看| 国产在线一区二区三区精| 欧美成人午夜精品| 成年人午夜在线观看视频| 亚洲 欧美一区二区三区| 久久热在线av| 中文字幕最新亚洲高清| 久久影院123| 99久久99久久久精品蜜桃| 757午夜福利合集在线观看| 岛国在线观看网站| 亚洲国产成人一精品久久久| 久久狼人影院| 国产精品二区激情视频| 欧美在线一区亚洲| 欧美乱妇无乱码| 十八禁人妻一区二区| 69精品国产乱码久久久| 2018国产大陆天天弄谢| 黑人猛操日本美女一级片| 伦理电影免费视频| 精品少妇内射三级| 99国产精品一区二区三区| 亚洲色图av天堂| 黄色 视频免费看| 老熟女久久久| 久久精品aⅴ一区二区三区四区| 精品人妻熟女毛片av久久网站| 五月开心婷婷网| 国产免费福利视频在线观看| 亚洲成a人片在线一区二区| 久久青草综合色| 一级黄色大片毛片| 亚洲少妇的诱惑av| 亚洲欧美色中文字幕在线| 十八禁网站免费在线| 一夜夜www| 久久久精品国产亚洲av高清涩受| 高清在线国产一区| 精品福利永久在线观看| 国产精品美女特级片免费视频播放器 | 国产在线一区二区三区精| 美女高潮喷水抽搐中文字幕| 夜夜爽天天搞| e午夜精品久久久久久久| 人人澡人人妻人| netflix在线观看网站| 久久婷婷成人综合色麻豆| 一区二区日韩欧美中文字幕| 国产精品av久久久久免费| 成人黄色视频免费在线看| 男女高潮啪啪啪动态图| 国产成+人综合+亚洲专区| 国产深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 别揉我奶头~嗯~啊~动态视频| 久久人人爽av亚洲精品天堂| 日韩大片免费观看网站| 三级毛片av免费| 啦啦啦视频在线资源免费观看| 亚洲精品中文字幕在线视频| 少妇被粗大的猛进出69影院| 免费观看av网站的网址| 国产麻豆69| 丁香六月欧美| 老汉色∧v一级毛片| kizo精华| 亚洲精品自拍成人| 免费av中文字幕在线| 欧美人与性动交α欧美软件| av天堂在线播放| av有码第一页| 女警被强在线播放| 成人精品一区二区免费| 在线亚洲精品国产二区图片欧美| 成年人免费黄色播放视频| 成人手机av| 欧美中文综合在线视频| 一边摸一边抽搐一进一出视频| 新久久久久国产一级毛片| 国产男女超爽视频在线观看| 一区在线观看完整版| 国产aⅴ精品一区二区三区波| 在线观看66精品国产| 一区二区日韩欧美中文字幕| av网站免费在线观看视频| 成人国产av品久久久| 亚洲午夜精品一区,二区,三区| 91字幕亚洲| 国产在线免费精品| 制服诱惑二区| av片东京热男人的天堂| 高潮久久久久久久久久久不卡| 精品久久久精品久久久| 99久久国产精品久久久| 黄色丝袜av网址大全| 国产老妇伦熟女老妇高清| 在线观看人妻少妇| 90打野战视频偷拍视频| 久久性视频一级片| 久久久精品94久久精品| 在线 av 中文字幕| videos熟女内射| av有码第一页| bbb黄色大片| 91精品国产国语对白视频| 精品国产超薄肉色丝袜足j| 国产精品麻豆人妻色哟哟久久| 在线永久观看黄色视频| 国产免费福利视频在线观看| 制服诱惑二区| 国产单亲对白刺激| 国产成人av激情在线播放| 热99国产精品久久久久久7| 黄色成人免费大全| 少妇粗大呻吟视频| 丝袜喷水一区| 91麻豆精品激情在线观看国产 | 亚洲精品国产色婷婷电影| 国产成人免费无遮挡视频| 一区二区三区国产精品乱码| 亚洲午夜理论影院| 免费日韩欧美在线观看| 午夜激情av网站| 热re99久久精品国产66热6| 女人被躁到高潮嗷嗷叫费观| 国产亚洲一区二区精品| 色综合婷婷激情| 我的亚洲天堂| 91精品国产国语对白视频| 久热爱精品视频在线9| 亚洲伊人色综图| 精品国产一区二区三区四区第35| 国产av一区二区精品久久| 深夜精品福利| 19禁男女啪啪无遮挡网站| 日韩视频一区二区在线观看| 亚洲七黄色美女视频| 久久久久国内视频| 波多野结衣一区麻豆| 菩萨蛮人人尽说江南好唐韦庄| 操美女的视频在线观看| 国产亚洲欧美精品永久| 伦理电影免费视频| 窝窝影院91人妻| 又黄又粗又硬又大视频| 999久久久国产精品视频| 蜜桃在线观看..| 大香蕉久久网| 久久中文字幕一级| 电影成人av| 日本av免费视频播放| 国产成人av教育| 亚洲精品乱久久久久久| 高清在线国产一区| 午夜福利影视在线免费观看| 妹子高潮喷水视频| 丁香欧美五月| 十八禁网站免费在线| 最近最新中文字幕大全免费视频| 久久香蕉激情| 午夜久久久在线观看| 多毛熟女@视频| 亚洲自偷自拍图片 自拍| 蜜桃国产av成人99| h视频一区二区三区| 在线天堂中文资源库| 高清黄色对白视频在线免费看| 欧美精品一区二区免费开放| 熟女少妇亚洲综合色aaa.| 中文字幕精品免费在线观看视频| 久久影院123| 免费一级毛片在线播放高清视频 | √禁漫天堂资源中文www| 免费在线观看日本一区| 国产精品一区二区精品视频观看| 日韩欧美国产一区二区入口| 久久亚洲精品不卡| 一个人免费在线观看的高清视频| 国产成人影院久久av| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久久免费视频了| 搡老熟女国产l中国老女人| 亚洲欧美精品综合一区二区三区| 在线十欧美十亚洲十日本专区| 久久国产精品影院| 一本综合久久免费| 午夜91福利影院| 一区二区三区激情视频| 热99久久久久精品小说推荐| 亚洲久久久国产精品| 亚洲精品中文字幕一二三四区 | 一区二区三区国产精品乱码| 国产精品久久久久成人av| 建设人人有责人人尽责人人享有的| 人妻 亚洲 视频| 叶爱在线成人免费视频播放| 一区在线观看完整版| 91麻豆av在线| 亚洲自偷自拍图片 自拍| 亚洲av成人不卡在线观看播放网| 99国产精品免费福利视频| 嫁个100分男人电影在线观看| 大码成人一级视频| 又黄又粗又硬又大视频| 男女无遮挡免费网站观看| 亚洲精品自拍成人| www.精华液| av不卡在线播放| 久久精品国产a三级三级三级| 中文字幕人妻丝袜一区二区| 国产成人精品久久二区二区91| 国产成人av教育| 久久99热这里只频精品6学生| 亚洲色图av天堂| 美女高潮喷水抽搐中文字幕| 纵有疾风起免费观看全集完整版| svipshipincom国产片| 丰满饥渴人妻一区二区三| 一本一本久久a久久精品综合妖精| 美女主播在线视频| 国产成人精品久久二区二区免费| 在线观看人妻少妇| 国产免费福利视频在线观看| 亚洲欧洲精品一区二区精品久久久| 国产精品1区2区在线观看. | 精品少妇一区二区三区视频日本电影| 一本一本久久a久久精品综合妖精| 欧美人与性动交α欧美精品济南到| 少妇猛男粗大的猛烈进出视频| 人人妻,人人澡人人爽秒播| 丝袜美腿诱惑在线| 亚洲精品乱久久久久久| 国产精品99久久99久久久不卡| 亚洲精品一二三| 视频在线观看一区二区三区| 国产精品欧美亚洲77777| 国产精品国产高清国产av | 涩涩av久久男人的天堂| 深夜精品福利| 一边摸一边抽搐一进一小说 | 国产成人免费观看mmmm| 国产成+人综合+亚洲专区| 亚洲国产欧美日韩在线播放| 国产高清激情床上av| 国内毛片毛片毛片毛片毛片| 国产精品美女特级片免费视频播放器 | 免费在线观看影片大全网站| 亚洲精品美女久久av网站| 亚洲熟女精品中文字幕| 久久精品国产综合久久久| 国产日韩欧美亚洲二区| 欧美日韩亚洲综合一区二区三区_| 极品人妻少妇av视频| 99九九在线精品视频| 国产精品麻豆人妻色哟哟久久| 久久99一区二区三区| 久久久久视频综合| 岛国在线观看网站| 国产精品国产av在线观看| 大香蕉久久成人网| 人人澡人人妻人| 少妇粗大呻吟视频| 久久av网站| 久久香蕉激情| 成人永久免费在线观看视频 | 一进一出抽搐动态| 制服人妻中文乱码| 国产日韩欧美视频二区| 国产av一区二区精品久久| 午夜福利在线观看吧| 怎么达到女性高潮| 乱人伦中国视频| 成人免费观看视频高清| 99热国产这里只有精品6| 天堂8中文在线网| 久久久久久久大尺度免费视频| 国产xxxxx性猛交| 叶爱在线成人免费视频播放| 国产精品免费大片| 亚洲精品国产区一区二| 国产欧美日韩一区二区三| 久久久久久久精品吃奶| 国产区一区二久久| 精品少妇内射三级| 国产极品粉嫩免费观看在线| 国产一区二区 视频在线| 在线观看66精品国产| 精品一区二区三卡| 中文字幕人妻丝袜制服| aaaaa片日本免费| 国产免费视频播放在线视频| 欧美日韩黄片免| 久久av网站| 日韩大码丰满熟妇| 黑人猛操日本美女一级片| 亚洲欧美精品综合一区二区三区| 国产又色又爽无遮挡免费看| 亚洲欧美一区二区三区黑人| 国产精品亚洲av一区麻豆| 欧美日韩av久久| 亚洲视频免费观看视频| 超色免费av| 久久ye,这里只有精品| 999精品在线视频| 五月开心婷婷网| bbb黄色大片| 免费在线观看视频国产中文字幕亚洲| 日韩精品免费视频一区二区三区| 欧美在线一区亚洲| 色视频在线一区二区三区| 老熟妇仑乱视频hdxx| 90打野战视频偷拍视频| 亚洲伊人久久精品综合| 热re99久久国产66热| 色在线成人网| 纵有疾风起免费观看全集完整版| 多毛熟女@视频| 动漫黄色视频在线观看| 欧美黑人欧美精品刺激| av福利片在线| 免费一级毛片在线播放高清视频 | 久久久精品区二区三区| 日韩欧美三级三区| 波多野结衣一区麻豆| 我的亚洲天堂| 久久久国产欧美日韩av| 在线看a的网站| 首页视频小说图片口味搜索| 99九九在线精品视频| 精品人妻1区二区| 亚洲av日韩精品久久久久久密| 精品久久久久久电影网| 亚洲va日本ⅴa欧美va伊人久久| 日本五十路高清| 少妇裸体淫交视频免费看高清 | 免费人妻精品一区二区三区视频| 国产片内射在线| 亚洲自偷自拍图片 自拍| 久久国产精品男人的天堂亚洲| 国产免费福利视频在线观看| 国产精品亚洲av一区麻豆| 亚洲欧美一区二区三区久久| 99国产极品粉嫩在线观看| 天天添夜夜摸| 夫妻午夜视频| 国产亚洲一区二区精品| 日韩三级视频一区二区三区| 99香蕉大伊视频| 99久久99久久久精品蜜桃| 伊人久久大香线蕉亚洲五| 黄色视频,在线免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 大片电影免费在线观看免费| 欧美精品啪啪一区二区三区| 国产又爽黄色视频| 亚洲精华国产精华精| 啦啦啦中文免费视频观看日本| 免费看a级黄色片| 亚洲一码二码三码区别大吗| 国产高清国产精品国产三级| 欧美+亚洲+日韩+国产| 狠狠狠狠99中文字幕| 脱女人内裤的视频| 国产日韩欧美亚洲二区| 国产熟女午夜一区二区三区| 一本大道久久a久久精品| 国产精品香港三级国产av潘金莲| 桃红色精品国产亚洲av| 丁香六月欧美| 天天影视国产精品| 亚洲黑人精品在线| 美女扒开内裤让男人捅视频| 亚洲成人免费电影在线观看| 国产午夜精品久久久久久| 十分钟在线观看高清视频www| 9热在线视频观看99| 高清视频免费观看一区二区| 精品国产乱码久久久久久小说| 国产成+人综合+亚洲专区| 久久亚洲精品不卡| 一区在线观看完整版| 美国免费a级毛片| 久久精品熟女亚洲av麻豆精品| 黄片大片在线免费观看| 久久国产亚洲av麻豆专区| 一个人免费看片子| 久久精品成人免费网站| 日日爽夜夜爽网站| 无遮挡黄片免费观看| 又大又爽又粗| 美女福利国产在线| 日韩 欧美 亚洲 中文字幕| 久久性视频一级片| 午夜日韩欧美国产| 精品视频人人做人人爽| 久久午夜亚洲精品久久| 美国免费a级毛片| 欧美午夜高清在线| 国产成人av教育| 黄色视频在线播放观看不卡| 老汉色av国产亚洲站长工具| 久久久久视频综合| 蜜桃国产av成人99| 51午夜福利影视在线观看| 精品一区二区三区av网在线观看 | 国产一卡二卡三卡精品| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产欧美日韩在线播放| 免费观看人在逋| 男女高潮啪啪啪动态图| 人人妻,人人澡人人爽秒播| a在线观看视频网站| 精品久久久精品久久久| 一二三四社区在线视频社区8| 美女福利国产在线| 亚洲欧美激情在线| 亚洲九九香蕉| 青青草视频在线视频观看| 亚洲国产看品久久| 国产麻豆69| 午夜两性在线视频| 免费观看av网站的网址| 天堂8中文在线网| 日本欧美视频一区| 亚洲精品中文字幕一二三四区 | 日韩中文字幕欧美一区二区| 1024香蕉在线观看| 国产精品秋霞免费鲁丝片| 国产精品1区2区在线观看. | 大型黄色视频在线免费观看| 成人手机av| 757午夜福利合集在线观看| 亚洲av电影在线进入| 国产日韩一区二区三区精品不卡| 精品亚洲成a人片在线观看| 日本av免费视频播放| a级毛片黄视频| 国产免费现黄频在线看| 99精品在免费线老司机午夜| 国产高清激情床上av| 精品一区二区三区四区五区乱码| 天堂俺去俺来也www色官网| 国产有黄有色有爽视频| 看免费av毛片| 国产91精品成人一区二区三区 | 老司机靠b影院| 1024视频免费在线观看| 国产精品久久久久成人av| 国产三级黄色录像| 国产精品一区二区在线不卡| tube8黄色片| 国产熟女午夜一区二区三区| 一边摸一边抽搐一进一出视频| 亚洲中文av在线| 久久热在线av| avwww免费| 极品人妻少妇av视频| 999久久久精品免费观看国产| 满18在线观看网站| 黄色毛片三级朝国网站| 亚洲视频免费观看视频| av又黄又爽大尺度在线免费看| 久久久久久人人人人人| 色精品久久人妻99蜜桃| 91九色精品人成在线观看| 91精品三级在线观看| 一级毛片电影观看| 波多野结衣av一区二区av| 欧美日韩福利视频一区二区| 久久国产精品影院| 亚洲七黄色美女视频| 大片免费播放器 马上看| 99国产精品免费福利视频| 9色porny在线观看| 欧美黄色片欧美黄色片| 国产片内射在线| 久热这里只有精品99| 婷婷丁香在线五月| 天天影视国产精品| 日韩大片免费观看网站| 日韩制服丝袜自拍偷拍| 中文字幕人妻熟女乱码| 亚洲精品中文字幕在线视频| 18在线观看网站| 精品少妇一区二区三区视频日本电影| 国产免费福利视频在线观看| 咕卡用的链子| 欧美日韩中文字幕国产精品一区二区三区 | 免费少妇av软件| 亚洲免费av在线视频| 999久久久国产精品视频| 久久精品亚洲av国产电影网| 日本黄色视频三级网站网址 | 老鸭窝网址在线观看| av不卡在线播放| 亚洲人成77777在线视频| 99精品在免费线老司机午夜| 十八禁网站免费在线| 国产在线精品亚洲第一网站| 在线永久观看黄色视频| √禁漫天堂资源中文www| 两个人免费观看高清视频| 熟女少妇亚洲综合色aaa.| 成年动漫av网址| 亚洲欧洲精品一区二区精品久久久| 国产91精品成人一区二区三区 | 欧美在线黄色| 午夜激情久久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 夜夜夜夜夜久久久久| 伊人久久大香线蕉亚洲五| 久久久国产成人免费| 另类精品久久| 69精品国产乱码久久久| 999久久久国产精品视频| av片东京热男人的天堂| 亚洲国产毛片av蜜桃av| 大码成人一级视频| 日韩欧美一区视频在线观看| 久久久久精品国产欧美久久久| 亚洲成人国产一区在线观看| 两个人看的免费小视频| 美女高潮喷水抽搐中文字幕| 99精品在免费线老司机午夜| 人人妻人人澡人人爽人人夜夜| 中文字幕精品免费在线观看视频| 最近最新中文字幕大全免费视频| 777久久人妻少妇嫩草av网站| 青草久久国产| 飞空精品影院首页| 精品一区二区三区四区五区乱码| 国产精品影院久久| 亚洲免费av在线视频| 一区福利在线观看| 精品第一国产精品| 人妻一区二区av| 91老司机精品| 日本精品一区二区三区蜜桃| 久久人人爽av亚洲精品天堂| 午夜福利,免费看| 黄片播放在线免费| 免费在线观看完整版高清| 国产一区二区 视频在线| 久久中文看片网| 天堂8中文在线网| 99国产精品免费福利视频| 午夜免费鲁丝| 日本黄色日本黄色录像| 丁香六月欧美| 久久国产精品大桥未久av| 一本久久精品| 男女免费视频国产| 激情在线观看视频在线高清 | 女人久久www免费人成看片| 精品高清国产在线一区| 国产精品自产拍在线观看55亚洲 | 一级毛片精品| 露出奶头的视频| 午夜福利影视在线免费观看| 日本a在线网址| 欧美日韩精品网址| 中国美女看黄片| 动漫黄色视频在线观看| 久久久国产成人免费|