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

    基于高分一號衛(wèi)星遙感影像的礦區(qū)生態(tài)安全評價研究
    ——以井陘礦區(qū)為例

    2015-06-07 11:09:03鋒,張慧,馬
    地理與地理信息科學 2015年5期
    關鍵詞:滑坡礦區(qū)指標

    劉 劍 鋒,張 可 慧,馬 文 才

    (1.河北省科學院地理科學研究所,河北 石家莊 050021;2.河北省地理信息開發(fā)應用工程技術研究中心,河北 石家莊 050021;3.河北省科學院,河北 石家莊 050081)

    ?

    基于高分一號衛(wèi)星遙感影像的礦區(qū)生態(tài)安全評價研究
    ——以井陘礦區(qū)為例

    劉 劍 鋒1,2,張 可 慧1,馬 文 才3

    (1.河北省科學院地理科學研究所,河北 石家莊 050021;2.河北省地理信息開發(fā)應用工程技術研究中心,河北 石家莊 050021;3.河北省科學院,河北 石家莊 050081)

    根據(jù)礦區(qū)生態(tài)環(huán)境的特點,以世界經(jīng)合組織(OCED)提出的“壓力-狀態(tài)-響應”(P-S-R) 框架理論為指導,構建了既適合遙感監(jiān)測技術定量評估又能較好表征礦區(qū)生態(tài)安全信息的評價指標體系,針對礦區(qū)開采地區(qū)空間地物的復雜性,提出了指標數(shù)據(jù)源遙感化的概念,基于“高分一號”衛(wèi)星遙感影像為主的多源、多時相、多分辨率遙感數(shù)據(jù)源,實現(xiàn)了多尺度目標監(jiān)測與評價技術。以井陘礦區(qū)為例,進行礦區(qū)生態(tài)安全評價研究,評價結果符合實際情況,表明該研究可為維護礦區(qū)及其周邊生態(tài)環(huán)境的原生平衡、引導礦區(qū)資源的合理開采與利用提供技術支持。

    高分一號;遙感;生態(tài)安全;礦區(qū);指標體系;評價

    0 引言

    隨著人類社會對地球資源的不斷開發(fā),生態(tài)環(huán)境的原生平衡不斷受到威脅,生態(tài)安全問題也日益凸顯。生態(tài)安全的研究愈來愈受到學術界、政府及公眾的廣泛關注,對生態(tài)安全概念的重新定義和擴展有了更多的研究和討論[1],從而開展了許多生態(tài)安全的相關研究[2-9]。礦產(chǎn)資源是人類社會賴以生存和發(fā)展必需的物質(zhì)基礎,而日益增速的礦山開發(fā)造成的生態(tài)環(huán)境問題也愈發(fā)嚴重,礦區(qū)生態(tài)安全已成為生態(tài)安全的一項重要研究內(nèi)容。

    礦山開發(fā)區(qū)域一般地形較為復雜、地物特征多樣,地物信息獲取不易,如果采用傳統(tǒng)的實地監(jiān)測方法,費時、費力又滿足不了實時監(jiān)測的要求。近年來,3S及空間數(shù)據(jù)庫技術的快速發(fā)展為礦區(qū)生態(tài)環(huán)境的相關研究提供了新的數(shù)據(jù)來源和技術手段,可及時、大范圍地監(jiān)測礦區(qū)生態(tài)環(huán)境狀況,具有傳統(tǒng)方法難以比擬的優(yōu)勢。相關研究中,馬玲[5]采用SPOT-5圖像和Quick Bird圖像對攀枝花釩鈦磁鐵礦-寶鼎煤礦區(qū)礦山開采引起的環(huán)境問題進行監(jiān)測,并利用模糊數(shù)學方法結合遙感和GIS技術對該區(qū)的生態(tài)環(huán)境質(zhì)量進行了評價;Pradhan等[6]基于遙感數(shù)據(jù)和GIS技術,采用神經(jīng)網(wǎng)絡模型對礦區(qū)滑坡區(qū)域生態(tài)安全進行評價;葉寶瑩等[7]通過MSS、TM(ETM)、中巴資源衛(wèi)星(CBERS)等多源多尺度遙感數(shù)據(jù),提取安太堡露天礦區(qū)的煤礦開采、土地破壞及土地復墾數(shù)據(jù),為大型露天礦區(qū)的生態(tài)安全評價奠定了基礎;朱運海等[8]在生態(tài)安全評價中采用動態(tài)評價因子,以發(fā)現(xiàn)生態(tài)安全在單個因子及整體上的發(fā)展趨勢和潛在問題。但由于礦區(qū)生態(tài)安全受礦業(yè)活動等的擾動性強,已有的評價方法需進一步完善和補充。

    針對礦山開采地區(qū)空間地物的復雜性,本文基于“高分一號”衛(wèi)星遙感影像為主的多源、多時相、多分辨率遙感數(shù)據(jù)源,探討不同因子的遙感反演技術方法,實現(xiàn)多尺度目標監(jiān)測與評價技術,采用“P-S-R”框架理論,構建快速有效的礦區(qū)生態(tài)安全評價指標體系,并對井陘礦區(qū)進行了實例研究,取得了較好的效果。本研究為維護礦區(qū)及其周邊生態(tài)環(huán)境的原生平衡、引導礦山資源的合理開采與利用提供技術支持。

    1 基于高分一號遙感影像的礦區(qū)生態(tài)安全評價方法

    1.1 礦區(qū)生態(tài)安全評價體系構建框架理論

    “P-S-R”(壓力-狀態(tài)-響應) 框架理論由世界經(jīng)合組織(OCED)提出,綜合考慮了人類活動、生態(tài)環(huán)境狀況及處理環(huán)境問題的政策措施三方面因素,可滿足礦區(qū)生態(tài)環(huán)境多目標、復雜的自然和人類活動評價要求。其同時考慮到社會經(jīng)濟、自然環(huán)境、人類活動等不同范疇的評價指標,是一種考慮比較全面、評價機制設計較完整的指導性評價框架理論,指標的選擇與評價對象選擇的理論依據(jù)、評價者的知識背景具有密切關系。

    本研究在“P-S-R”框架理論下展開,而利用遙感數(shù)據(jù)可以較好地提取相應環(huán)境因子數(shù)據(jù),因此研究采用我國首發(fā)的高時空分辨率遙感圖像——“高分一號”遙感影像為主的多源遙感數(shù)據(jù)源,同時輔以基礎數(shù)據(jù)庫中基礎、專題地理信息及社會經(jīng)濟統(tǒng)計數(shù)據(jù),結合儀器實地快速采樣分析法開展對研究區(qū)的監(jiān)測和評價研究,建立了基于3S技術的礦區(qū)生態(tài)安全P-S-R框架模型?;谶b感技術在空間、光譜、時間分辨率及信息提取技術方面的不斷進步,本研究提出了指標數(shù)據(jù)源遙感化的概念,即提高遙感數(shù)據(jù)在整個指標體系數(shù)據(jù)源中的占比,以期實現(xiàn)對研究區(qū)生態(tài)安全的逐年乃至逐月的監(jiān)測,為做到對整體生態(tài)安全或單因子的早期預警打下堅實基礎。

    1.2 礦區(qū)生態(tài)安全評價計算方法

    生態(tài)安全評價方法的提出經(jīng)歷了從簡單描述到精確定量計算的過程。目前常用的方法有灰色關聯(lián)度法、主成分投影法、層次分析法(AHP法)、生態(tài)足跡法、綜合指數(shù)法等。本文為更客觀地求取生態(tài)安全評價指數(shù),基于GIS技術,采用地學空間模型、AHP法和綜合指數(shù)法相結合的方法進行礦區(qū)生態(tài)安全評價。

    (1)評價指標的空間化與標準化。為更好地體現(xiàn)各評價指標在生態(tài)安全評價中的量化作用。需將其數(shù)據(jù)無量綱化,標準化處理方法如下:

    正相關指標的標準化(其指標值越大,則生態(tài)安全程度將越好,如植被蓋度):

    負相關指標的標準化(其指標值越大,則生態(tài)安全程度會越差,如礦業(yè)活動區(qū)面積):

    式中:xi′是指標i轉(zhuǎn)化后的無量綱值;xi為標準化前的原始指標值;xmax、xmin分別為區(qū)域內(nèi)最大、最小值[10]。

    其次,對非空間數(shù)據(jù)進行空間化處理。為了實現(xiàn)非空間數(shù)據(jù)與空間數(shù)據(jù)的疊加計算,需對其進行空間化處理。根據(jù)數(shù)據(jù)對應的行政區(qū)劃等級,分別將非空間數(shù)據(jù)值賦均值到鄉(xiāng)鎮(zhèn)或者行政村,然后進行柵格化處理,從而實現(xiàn)非空間數(shù)據(jù)的空間化處理。

    (2)礦區(qū)生態(tài)安全的綜合評價指數(shù)表達式為:

    Ai=∑KiBi

    式中:Ai代表綜合評價指數(shù);Ki為單項評價指標的權重值;Bi為單項評價指標的標準值。

    2 礦區(qū)生態(tài)安全評價指標體系構建

    礦區(qū)生態(tài)安全問題突出,其評價須利用已有評價標準并結合實際。指標構建原則如下:1)科學性原則。選擇能客觀、真實地反映礦區(qū)生態(tài)安全的現(xiàn)狀及演化特征的代表性指標。2)完整性與實用性原則。礦區(qū)屬自然和人工交互生態(tài)系統(tǒng),應全面考慮自然與社會的各項指標以完整反映礦區(qū)生態(tài)安全現(xiàn)狀及內(nèi)生機制等特征[11],同時考慮指標數(shù)據(jù)的易獲取性及可定量化性,并適應礦區(qū)發(fā)展水平。3)針對性和引導性原則。針對主要的生態(tài)安全問題設立指標,舍棄從屬目標,同時引導礦區(qū)生態(tài)安全的發(fā)展方向。4)動態(tài)性與穩(wěn)定性相結合原則。指標體系在一定時期后應做適當調(diào)整以適應礦區(qū)生態(tài)安全動態(tài)變化過程及社會經(jīng)濟的不同發(fā)展水平,以切實反映階段性特征,同時考慮評價的可對比性,指標體系在一段時間內(nèi)應保持相對穩(wěn)定性,原有指標和替代指標也可并存一段時間,以便于進行評價比較和趨勢分析。5)指標遙感化及快速獲取性原則。數(shù)據(jù)遙感化、數(shù)據(jù)庫支撐并采用其他先進數(shù)據(jù)采集技術可快速獲取客觀數(shù)據(jù),排除數(shù)據(jù)獲取周期長、實測數(shù)據(jù)缺乏等因素。6)數(shù)據(jù)空間化原則。各指標均具有空間屬性,反映其分布特征。

    根據(jù)上述原則,參考已有研究成果[9-13],構建包含壓力-狀態(tài)-響應3個層次的礦區(qū)生態(tài)安全評價指標體系,共涉及23個評價指標,并采用層次分析法結合專家打分賦予指標權重(表1)。

    3 礦區(qū)生態(tài)安全評價指標信息提取

    3.1 基于多源遙感數(shù)據(jù)的指標數(shù)據(jù)提取

    (1)遙感數(shù)據(jù)源選擇。考慮獲取遙感數(shù)據(jù)的經(jīng)費、時空分辨率及成像幅面寬度等因素,研究選用高分一號遙感影像數(shù)據(jù),并輔以中巴資源衛(wèi)星遙感數(shù)據(jù)、TM遙感數(shù)據(jù)以及GoogleEarth、天地圖等網(wǎng)絡遙感資源進行礦區(qū)相關數(shù)據(jù)的提取。高分一號衛(wèi)星是中國高分辨率對地觀測系統(tǒng)的第一顆衛(wèi)星,其突破了多光譜、高時空分辨率相結合的光學遙感技術,高精度、高穩(wěn)定度姿態(tài)控制技術,高分辨率數(shù)據(jù)處理與應用等關鍵技術,適合礦區(qū)生態(tài)環(huán)境關鍵因子的監(jiān)測。

    表1 評價指標體系及權重

    Table 1 Evaluation index system and weight

    目標層準則層[權重]要素層[權重]指標層[權重]數(shù)據(jù)來源評價指標體系壓力[0.44](Pressure)狀態(tài)[0.46](State)響應[0.10](Response)環(huán)境壓力[0.34]人文壓力[0.10]地質(zhì)構造[0.05]巖性[0.05]地形[0.10]礦區(qū)環(huán)境可治理度[0.05]土壤侵蝕度[0.06]植被指數(shù)[0.10]生物多樣性指標[0.05]社會響應指標[0.10]地質(zhì)災害數(shù)量[0.06]遙感數(shù)據(jù)地質(zhì)災害潛在危險源[0.05]遙感數(shù)據(jù)水體污染[0.09]遙感數(shù)據(jù)大氣粉塵污染[0.02]遙感數(shù)據(jù)大氣污染[0.04]監(jiān)測數(shù)據(jù)土壤化肥折純量[0.02]數(shù)據(jù)庫農(nóng)藥使用量[0.03]數(shù)據(jù)庫農(nóng)膜使用量[0.01]數(shù)據(jù)庫土壤Pb污染[0.02]實測數(shù)據(jù)礦業(yè)活動區(qū)面積[0.03]遙感數(shù)據(jù)>25°坡耕地面積[0.02]遙感數(shù)據(jù)人均耕地面積[0.01]數(shù)據(jù)庫人口密度[0.02]數(shù)據(jù)庫恩格爾系數(shù)[0.01]數(shù)據(jù)庫人口受教育程度[0.01]數(shù)據(jù)庫地質(zhì)構造[0.05]數(shù)據(jù)庫巖性[0.05]數(shù)據(jù)庫坡度[0.10]數(shù)據(jù)庫礦區(qū)環(huán)境可治理度[0.05]遙感數(shù)據(jù)土壤侵蝕度[0.06]數(shù)據(jù)庫NDVI[0.10]遙感數(shù)據(jù)生物豐度指數(shù)[0.05]遙感數(shù)據(jù)生態(tài)環(huán)境治理度[0.10]遙感數(shù)據(jù)

    (2)遙感數(shù)據(jù)預處理。一般情況下,礦區(qū)內(nèi)地貌類型復雜多樣,相對高差大,利用DEM數(shù)據(jù)對經(jīng)過粗校正的遙感圖像進行正射校正以消除高差等因素帶來的畸變;其后對高分一號2 m全色和8 m多光譜影像經(jīng)IHS變換后的I分量經(jīng)小波變換進行融合,得到兼具高空間分辨率和光譜分辨率的融合圖像。

    (3)指標信息提取。通過建立LUCC數(shù)據(jù)解譯標志并細化礦產(chǎn)活動占地解譯標志形成了融合圖像的LUCC遙感數(shù)據(jù)信息提取標志。其他信息解譯標志如下:1)地質(zhì)災害:主要包括滑坡、崩塌等?;乱话愠属せ巍⑸嘈?、梨形等,陡峭的滑坡壁和它們所形成的圍谷,遙感影像上表現(xiàn)為彎曲弧形。2)地質(zhì)災害危險性:采礦等活動造成掌子面、尾礦庫等固廢堆放,存在潰壩、崩塌等隱患,根據(jù)礦業(yè)采場及堆場位置結合DEM數(shù)據(jù),劃定隱患大致范圍。3)環(huán)境污染:選礦、洗煤廢水不經(jīng)處理或處理不達標排放,對水體或土壤等造成污染。采礦過程中產(chǎn)生的煙塵、粉塵主要對礦業(yè)活動區(qū)臨近區(qū)域及運輸?shù)缆穬蓚?cè)造成污染,表現(xiàn)在土壤或植被的顏色與周圍正常顏色的明顯差異。

    3.2 基于實地及其他來源的指標數(shù)據(jù)快速提取

    便攜式X射線熒光光譜儀對于Cu、Pb等重金屬元素的檢出限均低于區(qū)域元素背景值和土壤環(huán)境質(zhì)量標準一級標準值,說明便攜式X射線熒光光譜儀可用于土壤重金屬污染快速檢測和土地質(zhì)量評價[14]。實地采集的數(shù)據(jù)為土壤重金屬污染數(shù)據(jù)。采用X射線熒光分析儀,在野外按照1.5 km×1.5 km的網(wǎng)格布點,DEM坡度及地質(zhì)、斷裂帶、統(tǒng)計數(shù)據(jù)等數(shù)據(jù)來自河北省科學院地理科學研究所不斷更新的生態(tài)環(huán)境基礎數(shù)據(jù)庫。

    4 實例研究及結果分析

    4.1 研究區(qū)概況

    井陘礦區(qū)是河北省石家莊市的直轄區(qū),位于太行山區(qū)中段,北緯38°01′18″~38°08′03″,東經(jīng)113°58′50″~114°06′08″,周邊與井陘縣接壤,地勢西高東低,面積70.29 km2,轄兩鎮(zhèn)一鄉(xiāng),2個街道辦事處,30個行政村。礦區(qū)礦產(chǎn)資源以煤炭為主,具有一百多年的礦產(chǎn)資源開采史,該區(qū)經(jīng)濟曾因礦山開發(fā)而繁榮,在20世紀60年代,其焦煤產(chǎn)量位居全國第二位。

    井陘礦區(qū)由于面積較小,生態(tài)系統(tǒng)屬于人類社會經(jīng)濟系統(tǒng)和自然生態(tài)系統(tǒng)交互影響的復合生態(tài)系統(tǒng)。在利益等因素驅(qū)動下,資源遭到長時間、大規(guī)模、快速的開采利用,致使礦區(qū)礦產(chǎn)等資源枯竭、環(huán)境污染,極大地超過該區(qū)域生態(tài)承載能力,使當?shù)厣鷳B(tài)安全面臨巨大壓力,對區(qū)域生態(tài)安全構成嚴重威脅。而區(qū)域生態(tài)安全惡化后,又反過來制約社會經(jīng)濟的發(fā)展,形成惡性循環(huán)。井陘礦區(qū)煤炭資源已經(jīng)逐步枯竭,是一個典型的“因煤而建、因煤而興、因煤而衰”的資源型城區(qū),其礦山開發(fā)在給該區(qū)帶來經(jīng)濟上繁榮的同時也引發(fā)較多的生態(tài)環(huán)境問題,具有典型的代表性,故選擇該區(qū)作為研究區(qū)進行礦區(qū)生態(tài)安全評價,監(jiān)控和掌握該區(qū)域的生態(tài)安全現(xiàn)狀,并通過分析生態(tài)安全問題產(chǎn)生的原因為制定相應的對策打下基礎,最終達到使研究區(qū)生態(tài)安全處于良好狀態(tài)、實現(xiàn)可持續(xù)發(fā)展的目的。

    4.2 主要生態(tài)安全指標分析

    4.2.1 地質(zhì)災害數(shù)量分析 通過對高分一號衛(wèi)星數(shù)據(jù)的解譯判讀,輔助以Google earth等網(wǎng)絡高分辨率(1 m級)遙感影像、坡度數(shù)據(jù)及實地考察驗證。研究區(qū)地質(zhì)災害分布情況見表2。

    表2 研究區(qū)滑坡、崩塌分布

    Table 2 Distribution of landslide and avalanche

    分布區(qū)域地質(zhì)災害類型地質(zhì)災害數(shù)量賈莊鎮(zhèn)滑坡1崩塌2鳳山鎮(zhèn)滑坡7橫澗鄉(xiāng)--

    從地形上看,地質(zhì)災害主要分布在西部清涼山和云鳳山上;從行政區(qū)域看,主要分布在賈莊鎮(zhèn)和鳳山鎮(zhèn)。地質(zhì)災害的主要類型為滑坡,共有8個點;其次為崩塌,共有2個點。橫澗鄉(xiāng)由于僅有小部分較低坡度山地,因此沒有滑坡、崩塌發(fā)生。據(jù)四川攀西調(diào)查資料[15],按平均坡度分級進行統(tǒng)計如表3 所示,滑坡大多發(fā)生在10°~30°之間,崩塌多數(shù)出現(xiàn)在坡度大于30°的邊坡。

    表3 崩塌、滑坡與坡度的關系[16]

    Table 3 Relationship between avalanche,landslide and slope

    坡度(°)數(shù)量(個)比例(%)類型10沒有0滑坡10~20607.4滑坡21~3053265.2滑坡>3022427.4崩塌

    文獻[17]對過去30年間的2 238處崩塌歷史資料進行了統(tǒng)計,約80%的崩塌事例的坡度為30°~50°,其中邊坡坡度在40°時最多,并且一般土質(zhì)邊坡的坡度要緩于巖質(zhì)邊坡,且土質(zhì)邊坡崩塌時的坡度主要介于30°~40°之間,而巖質(zhì)邊坡崩塌時坡度在30°~50°之間[18]。本研究將滑坡、崩塌發(fā)生地點與坡度數(shù)據(jù)進行了疊加分析顯示,發(fā)生滑坡的地方坡度在12.5°~23.0°之間,崩塌一處坡度在31.0°~33.0°。影像判讀解譯結果與上述統(tǒng)計資料中滑坡、崩塌發(fā)生的坡度分布規(guī)律吻合。

    4.2.2 地質(zhì)災害規(guī)模分析 采空區(qū)在井陘礦區(qū)核心區(qū)域大面積分布(圖1)。各鄉(xiāng)鎮(zhèn)中采空區(qū)面積占總面積的比例最高的是鳳山鎮(zhèn),約25.7%;其次為賈莊鎮(zhèn),占21%;最少的是橫澗鄉(xiāng),占12.6%。采空區(qū)對其區(qū)域的各種土地利用類型均產(chǎn)生不同程度的影響。由表4可知,采空區(qū)影響面積最大的是耕地,約682.2hm2,占整個受影響區(qū)域面積的47.28%;其次是林地和中轉(zhuǎn)場地,分別占16.56%和12.62%;城鄉(xiāng)居民用地為5.58%。通過以上數(shù)據(jù)可知,采空區(qū)導致的地面塌陷、地裂縫等對農(nóng)業(yè)、林業(yè)、工業(yè)及居民生活均產(chǎn)生了較大影響。通過相關資料和實地調(diào)查得知,由于采空區(qū)地表應力發(fā)生改變,產(chǎn)生的塌陷、地裂縫等十分嚴重,導致居民及工礦企業(yè)建筑受損,棄耕撂荒現(xiàn)象嚴重。

    圖1 井陘礦區(qū)采空區(qū)分布

    Fig.1 Goaf distribution of Jingxing mining area

    表4 采空區(qū)土地利用類型及其面積

    Table 4 Area of land use in goaf area

    土地利用類型面積(m2)占采空區(qū)比例(%)耕地6827195.5847.275林地2391024.7416.557草地686982.454.757水域283826.331.965城鄉(xiāng)居民用地805681.575.579中轉(zhuǎn)場地1822223.0012.618固體廢棄物 196.780.001礦業(yè)建筑1249353.088.651交通工業(yè)用地374882.852.597總計14441366.39100

    4.2.3 地質(zhì)災害危險性分析 開采煤炭和地下水使地表產(chǎn)生不均勻沉降,從而引發(fā)塌陷、地裂縫等地質(zhì)災害,在煤炭資源幾乎枯竭的情況下,地方國有煤礦和鄉(xiāng)鎮(zhèn)村集體煤礦多次復采殘煤,進一步加大了地面塌陷、裂縫程度,造成地質(zhì)災害隱患;采石場等破壞了原有的山體平衡,形成掌子面等開采面,受雨水沖刷容易產(chǎn)生崩塌和滑坡等地質(zhì)災害;采石形成的尾礦、廢礦石如果貯存地點不適宜會導致堆積量逐年增大,在汛期容易產(chǎn)生泥石流;雨量年度分布極為不均,汛期強降雨易誘發(fā)地質(zhì)災害的發(fā)生。因此,采空區(qū)、礦山開采面、不適宜的尾礦庫和廢礦石堆放點,都具有地質(zhì)災害危險性。井陘礦區(qū)采礦歷史悠久,存在大面積的采空區(qū),經(jīng)遙感判讀在賈莊鎮(zhèn)和橫澗鄉(xiāng)存在大量礦山開采面,這些都是崩塌、滑坡、地面塌陷、地裂縫的潛在危險源。因此井陘礦區(qū)具有較高的地質(zhì)災害危險性。

    4.2.4 粉塵污染分析 井陘礦區(qū)存在大量煤炭開采、洗選、儲存、中轉(zhuǎn)等活動,加之大量采石場的存在,導致采場、中轉(zhuǎn)場地等礦業(yè)活動場地周圍及交通主干道兩側(cè)粉塵污染極為嚴重,對周邊區(qū)域土壤產(chǎn)生了嚴重影響,受污染農(nóng)作物導致減產(chǎn),受污染植被導致生長不良,并且視覺感官極差。相關區(qū)域在遙感影像上呈現(xiàn)以礦業(yè)活動場所為中心向外顏色由深變淺的規(guī)律,即離煤炭較近的地方受煤炭粉塵污染較重,隨著距離的增加影響逐漸減少;而采石場周邊則因為受淺色粉塵的影響呈現(xiàn)較周圍同類地物淺的色調(diào)(圖略)。

    4.3 礦區(qū)生態(tài)安全綜合評價結果分析

    依據(jù)公開出版的國家、地方和行業(yè)標準,結合已有研究成果,將研究區(qū)生態(tài)安全狀態(tài)分為安全狀態(tài)、較安全狀態(tài)、一般狀態(tài)、較不安全狀態(tài)和不安全狀態(tài)5個生態(tài)安全等級。根據(jù)綜合評價結果可知,研究區(qū)內(nèi)生態(tài)安全評價等級為一般狀態(tài)的面積最大,約占30.17%;其次為較安全級別,約占22.41%,安全狀態(tài)占17.24%;不安全狀態(tài)和較不安全狀態(tài)分別占15.81%和14.37%。由此可知該區(qū)域生態(tài)安全已受到威脅的區(qū)域面積占比為60.35%,生態(tài)安全問題應引起相關部門關注。

    由研究區(qū)生態(tài)安全狀態(tài)空間分布(圖2)可知,井陘的生態(tài)安全狀況總體格局呈現(xiàn)以下特點:1)研究區(qū)5種生態(tài)安全狀態(tài)呈現(xiàn)按鄉(xiāng)鎮(zhèn)集聚的特點。鳳山鎮(zhèn)生態(tài)安全狀態(tài)好于橫澗鄉(xiāng),而賈莊鎮(zhèn)生態(tài)安全狀態(tài)最差。時空分布上研究區(qū)北部處于較差的生態(tài)安全水平,南部整體處于較好的生態(tài)安全水平。2)研究區(qū)生態(tài)安全水平空間分布與人類活動分布和地形特點呈正相關。賈莊鎮(zhèn)人口最多,山地占比例最高,礦山開采活動密度大,工業(yè)較為發(fā)達,因此該區(qū)域生態(tài)安全承受了較大壓力,大部分區(qū)域處于預警狀態(tài)。

    綜合以上分析,大規(guī)模礦業(yè)活動是使得礦區(qū)內(nèi)部分區(qū)域生態(tài)安全處于較不安全或不安全狀態(tài)的主因。從生態(tài)安全受到威脅的區(qū)域看,很大程度上歸咎于該區(qū)域規(guī)模大、強度高、近百年的礦業(yè)開采歷史。礦區(qū)的資源開采型支柱產(chǎn)業(yè)單一,圍繞礦業(yè)開采的產(chǎn)業(yè)活動造成的地質(zhì)災害及隱患、景觀的破壞、土地資源占用、環(huán)境污染等負面影響導致該區(qū)域生態(tài)安全出現(xiàn)大面積預警狀態(tài)。

    圖2 生態(tài)安全空間分布格局

    Fig.2 Spatial distribution of ecological security

    5 結語

    本文首次將以高分一號衛(wèi)星遙感影像為主的多類型、多時相和不同空間分辨率的遙感圖像引入數(shù)據(jù)源,基于“壓力-狀態(tài)-響應”框架理論建立了一套適合礦區(qū)生態(tài)安全評價的遙感信息指標體系,涉及的23個指標涵蓋了自然、人類活動、社會經(jīng)濟等多方面的因素,是一個較全面的礦區(qū)生態(tài)安全評價體系,指標數(shù)據(jù)均來自遙感數(shù)據(jù)、數(shù)據(jù)庫、監(jiān)測數(shù)據(jù)及實地快測。遙感信息提取采用了計算機自動提取與目視解譯相結合的方法,既提高了信息提取的效率又保證了解譯精度,可滿足對礦區(qū)生態(tài)安全的快速監(jiān)測及評價的要求。經(jīng)實例驗證,本文方法能較好地評價實際生態(tài)安全狀況,可為相關研究提供借鑒。由于在指標權重的設定中采用了層次分析法與專家打分法相結合的方法,存在一定的主觀成分,如何取得更為客觀的權重體系,提高評價的普適性,將是下一步研究的重點。

    [1] ECHERSLEY R.Ecological security dilemmas[EB/OL].http://www.arts.monash.edu.au/ncas/teach/unit/pol/chpt08.html.2013-04-10.

    [2] 俞孔堅,高中潮.景觀生態(tài)和環(huán)境保護規(guī)劃的生態(tài)安全格式途徑[J].陜西環(huán)境,1999,6(4):42-45.

    [3] DOLBY S.Security and ecology in the age of globalization[R].ESCP Report 8, 2002.95-108.

    [4] 王讓虎,李曉燕,張樹文,等.東北農(nóng)牧交錯帶景觀生態(tài)安全格局構建及預警研究——以吉林省通榆縣為例[J].地理與地理信息科學,2014,30(2):111-115.

    [5] 馬玲.遙感與GIS技術在礦山環(huán)境監(jiān)測與質(zhì)量評價中的應用——以攀枝花釩鈦磁鐵礦寶鼎煤礦區(qū)為例[D].成都:成都理工大學,2008.

    [6] PRADHAN B, LEE S.Utilization of optical remote sensing data and GIS tools for regional landslide hazard analysis using an Artificial Neural Network Model[J].Earth Science Frontiers,2007,14(6):143-152.

    [7] 葉寶瑩,白中科,孔登魁,等.安太堡露天煤礦土地破壞與土地復墾動態(tài)變化的遙感調(diào)查[J].北京科技大學學報,2008,30(9):972-976.

    [8] 朱運海, 張百平, 曹銀璇,等.基于動態(tài)因子的區(qū)域生態(tài)安全評價——以河北省萬全縣為例[J].地理科學進展,2006,25(4):34-40.

    [9] 劉艷艷,吳大放,王朝暉.濕地生態(tài)安全評價研究進展[J].地理與地理信息科學, 2011,27(1):69-75.

    [10] 邵懷勇.攀西礦業(yè)開發(fā)集中區(qū)生態(tài)環(huán)境遙感信息提取及生態(tài)安全評價研究[D].成都:成都理工大學,2009.

    [11] 楊靜,王立芹.礦區(qū)生態(tài)安全評價指標體系的研究[J].山東科技大學學報(自然科學版),2005,24(l):36-39.

    [12] 左偉.基于RS、GIS和Models的區(qū)域生態(tài)環(huán)境系統(tǒng)安全綜合評價研究——以長江三峽庫區(qū)重慶市忠縣為例[D].南京:南京師范大學,2002.

    [13] DZ/T223-2007,礦山環(huán)境保護與綜合治理方案編制規(guī)范[S].中華人民共和國國土資源部,2007.

    [14] 冉景,王德建,王燦,等.便攜式X射線熒光光譜法與原子吸收/原子熒光法測定土壤重金屬的對比研究[J].光譜學與光譜分析,2010,1(83):127-129.

    [15] 譚萬沛,王成華.暴雨泥石流滑坡的區(qū)域預測與預報——以攀西地區(qū)為例[M].成都:四川科技出版社,1994.

    [16] 徐永年,匡尚富,李文武,等.邊坡形狀對崩塌的影響[J].泥沙研究,1999(5):67-73.

    [17] 申潤植.坡面框架工程設計與施工方法[M].日本三海堂,1995.

    [18] 唐邦興.山洪泥石流滑坡災害及防治[M].北京:科學出版社,1994.

    Research on Evaluation of Mining Area Ecological Security Based on Gaofen-1 Satellite Imagery:Take Jingxing Mining Area for Example

    LIU Jian-feng1,2,ZHANG Ke-hui1,MA Wen-cai3

    (1.InstituteofGeographicalSciences,HebeiAcademyofSciences,Shijiazhuang050021; 2.HebeiEngineeringResearchCenterforGeographicInformationApplication,Shijiazhuang050021; 3.HebeiAcademyofSciences,Shijiazhuang050081,China)

    According to the features of mining area environment,an evaluation index system of ecological security,representing mining area ecological information as well as suitable for quantified evaluation by remote sensing,is constructed based on the OCED(Organization for Economic Cooperation and Development) proposed "Pressure-State-Response" framework model.Aiming at the complexity of spatial ground objects in mining area,a concept of "remote sensing informatization for index data source" is put forward,which is to increase the proportion of remote sensing based index data,to realize a yearly or even monthly monitoring of mining area.Based on Gaofen-1 satellite imagery and other remote sensing data,remote sensing inversion methods of different index are discussed and multi-scale objects monitoring and evaluating technology are realized.Taking Jingxing mining area for example,the reasonable result of mining area ecological security evaluation shows that the study can provide support to maintain the original balance of mining area ecological environment and to guide the reasonable exploitation of mine resources.

    Gaofen-1 satellite;remote sensing;ecological security;mining area;index system;evaluation

    2014-10-11;

    2015-01-12

    河北省科學院科技計劃項目(14123、14130、15130)

    劉劍鋒(1974-),男,副研究員,研究方向為3S技術在生態(tài)環(huán)境中的應用。E-mail:ljf70s@126.com

    10.3969/j.issn.1672-0504.2015.05.024

    X37;TP75

    A

    1672-0504(2015)05-0121-06

    猜你喜歡
    滑坡礦區(qū)指標
    滑坡推力隱式解與顯式解對比分析——以河北某膨脹土滑坡為例
    加納Amanforom礦區(qū)Ⅲ號隱伏金礦帶的發(fā)現(xiàn)與評價
    加納Amanforom礦區(qū)Ⅲ號隱伏金礦帶的發(fā)現(xiàn)與評價
    湖北省??悼h堰邊上礦區(qū)發(fā)現(xiàn)超大型磷礦
    廣東省蕉嶺縣作壁坑礦區(qū)探明超大型銣礦
    最新引用指標
    莫讓指標改變初衷
    商周刊(2018年26期)2018-12-29 12:56:00
    淺談公路滑坡治理
    北方交通(2016年12期)2017-01-15 13:52:59
    基于Fluent的滑坡入水過程數(shù)值模擬
    “監(jiān)管滑坡”比“渣土山”滑坡更可怕
    山東青年(2016年3期)2016-02-28 14:25:50
    国产欧美日韩一区二区精品| 国产三级在线视频| 国产一区二区激情短视频| 一边摸一边抽搐一进一小说| 99热这里只有精品一区| 免费无遮挡裸体视频| 日韩成人在线观看一区二区三区| 欧美性猛交╳xxx乱大交人| 亚洲av免费高清在线观看| 亚洲第一电影网av| 三级男女做爰猛烈吃奶摸视频| 五月伊人婷婷丁香| 色播亚洲综合网| e午夜精品久久久久久久| 美女被艹到高潮喷水动态| 亚洲精品一卡2卡三卡4卡5卡| 五月伊人婷婷丁香| 欧美黑人巨大hd| 久久久国产成人精品二区| 亚洲成av人片在线播放无| 老熟妇仑乱视频hdxx| 中文字幕人妻丝袜一区二区| 欧美成狂野欧美在线观看| 午夜久久久久精精品| 黄色日韩在线| 国产精品av视频在线免费观看| 超碰av人人做人人爽久久 | 美女被艹到高潮喷水动态| 国内精品美女久久久久久| 免费电影在线观看免费观看| a在线观看视频网站| 欧美国产日韩亚洲一区| 亚洲人成网站在线播放欧美日韩| 国产成年人精品一区二区| 久久久久九九精品影院| 一本精品99久久精品77| 我的老师免费观看完整版| 99精品欧美一区二区三区四区| a在线观看视频网站| 精品人妻一区二区三区麻豆 | 90打野战视频偷拍视频| 成年女人永久免费观看视频| 噜噜噜噜噜久久久久久91| 国产成人aa在线观看| 99热精品在线国产| 久久欧美精品欧美久久欧美| 国产一区二区三区在线臀色熟女| 精品午夜福利视频在线观看一区| 国产免费男女视频| 草草在线视频免费看| 免费搜索国产男女视频| 在线播放国产精品三级| 国产极品精品免费视频能看的| 国产私拍福利视频在线观看| 日韩有码中文字幕| 人人妻人人看人人澡| 午夜福利高清视频| 免费看美女性在线毛片视频| 成人亚洲精品av一区二区| www日本在线高清视频| 美女大奶头视频| 午夜亚洲福利在线播放| 男人和女人高潮做爰伦理| 男女那种视频在线观看| 精品无人区乱码1区二区| 99久久综合精品五月天人人| 99热这里只有是精品50| 一级黄色大片毛片| 久久久久久九九精品二区国产| 神马国产精品三级电影在线观看| 999久久久精品免费观看国产| 亚洲va日本ⅴa欧美va伊人久久| 一本综合久久免费| 国产精品99久久99久久久不卡| 又紧又爽又黄一区二区| 亚洲成人免费电影在线观看| 不卡一级毛片| 悠悠久久av| 91av网一区二区| 久久99热这里只有精品18| 国产精品三级大全| 桃色一区二区三区在线观看| 在线观看舔阴道视频| 国产精品影院久久| 两个人视频免费观看高清| 亚洲av电影不卡..在线观看| 久久香蕉精品热| 亚洲,欧美精品.| 真人一进一出gif抽搐免费| 免费在线观看日本一区| 亚洲av成人不卡在线观看播放网| 狠狠狠狠99中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 白带黄色成豆腐渣| 国产精品免费一区二区三区在线| 国产真实乱freesex| 99久久精品一区二区三区| 日韩精品青青久久久久久| 亚洲性夜色夜夜综合| 国产69精品久久久久777片| 特级一级黄色大片| 久久精品国产自在天天线| 欧美成人a在线观看| 中文字幕人妻熟人妻熟丝袜美 | 亚洲国产色片| 又黄又粗又硬又大视频| 一区福利在线观看| 亚洲最大成人手机在线| 亚洲精品乱码久久久v下载方式 | 少妇人妻精品综合一区二区 | 久久精品91蜜桃| 亚洲国产欧洲综合997久久,| 国产中年淑女户外野战色| 成人精品一区二区免费| 少妇丰满av| 日韩欧美精品免费久久 | av在线蜜桃| 久久精品综合一区二区三区| 狂野欧美激情性xxxx| 老熟妇乱子伦视频在线观看| 色吧在线观看| 免费在线观看亚洲国产| 黄色日韩在线| 亚洲熟妇中文字幕五十中出| 欧美丝袜亚洲另类 | 亚洲中文日韩欧美视频| 国内揄拍国产精品人妻在线| 久9热在线精品视频| 色综合婷婷激情| 成人三级黄色视频| 欧美乱妇无乱码| 真人做人爱边吃奶动态| 一a级毛片在线观看| 日韩高清综合在线| 亚洲成av人片在线播放无| 亚洲精品成人久久久久久| 国产精品久久视频播放| 美女被艹到高潮喷水动态| 国产野战对白在线观看| 国产精品 国内视频| 两个人看的免费小视频| 婷婷精品国产亚洲av在线| а√天堂www在线а√下载| 法律面前人人平等表现在哪些方面| 午夜两性在线视频| 亚洲无线在线观看| 99久久成人亚洲精品观看| av在线蜜桃| av在线天堂中文字幕| 日韩欧美在线二视频| 国内精品美女久久久久久| 久久久精品欧美日韩精品| 色视频www国产| 亚洲中文日韩欧美视频| 亚洲无线在线观看| 国产伦人伦偷精品视频| 91久久精品国产一区二区成人 | 国产成人av教育| 婷婷丁香在线五月| 日韩欧美国产在线观看| 香蕉久久夜色| 亚洲国产高清在线一区二区三| 91av网一区二区| 嫩草影院精品99| 久久久色成人| 成年女人毛片免费观看观看9| 美女被艹到高潮喷水动态| 日本成人三级电影网站| 国产一区二区在线观看日韩 | 欧美中文综合在线视频| 97碰自拍视频| 最后的刺客免费高清国语| 亚洲av中文字字幕乱码综合| 美女免费视频网站| 亚洲人与动物交配视频| 一级黄色大片毛片| 少妇裸体淫交视频免费看高清| 国产高清videossex| 精品人妻偷拍中文字幕| 久久精品91无色码中文字幕| 91在线观看av| 黄片小视频在线播放| 国产美女午夜福利| 国产老妇女一区| 久久久久久久午夜电影| 国产精品永久免费网站| 免费看a级黄色片| 母亲3免费完整高清在线观看| 精品一区二区三区av网在线观看| 亚洲精品成人久久久久久| 禁无遮挡网站| 成人国产综合亚洲| 露出奶头的视频| 国产精品久久电影中文字幕| 亚洲电影在线观看av| 日本精品一区二区三区蜜桃| 欧美高清成人免费视频www| 亚洲无线观看免费| 亚洲 国产 在线| 琪琪午夜伦伦电影理论片6080| 亚洲人与动物交配视频| 丁香六月欧美| 在线视频色国产色| 男女视频在线观看网站免费| 国产视频内射| 国产av一区在线观看免费| 色老头精品视频在线观看| 国产视频一区二区在线看| 国产精品香港三级国产av潘金莲| 99在线人妻在线中文字幕| 一区二区三区免费毛片| 欧美丝袜亚洲另类 | 久久久久久久久大av| 亚洲欧美日韩卡通动漫| 色视频www国产| 亚洲第一电影网av| 脱女人内裤的视频| 亚洲电影在线观看av| 久久婷婷人人爽人人干人人爱| 男女之事视频高清在线观看| av片东京热男人的天堂| 国产在线精品亚洲第一网站| 亚洲国产欧美人成| 亚洲人成电影免费在线| 久久人妻av系列| 午夜福利在线观看免费完整高清在 | 亚洲在线观看片| 三级男女做爰猛烈吃奶摸视频| 人人妻人人澡欧美一区二区| 天堂动漫精品| 欧美成人一区二区免费高清观看| 99热6这里只有精品| 99热这里只有是精品50| 悠悠久久av| 日韩欧美三级三区| 亚洲黑人精品在线| 最新美女视频免费是黄的| 国产精品综合久久久久久久免费| av国产免费在线观看| 91久久精品国产一区二区成人 | 成年版毛片免费区| 一a级毛片在线观看| 老鸭窝网址在线观看| 99热这里只有是精品50| 在线天堂最新版资源| 女生性感内裤真人,穿戴方法视频| av片东京热男人的天堂| 亚洲国产精品合色在线| www日本黄色视频网| 在线观看午夜福利视频| 欧美最黄视频在线播放免费| 亚洲自拍偷在线| 欧美成人a在线观看| 精品国内亚洲2022精品成人| 国产成人系列免费观看| 午夜老司机福利剧场| 国产av麻豆久久久久久久| 欧美又色又爽又黄视频| 欧美黑人欧美精品刺激| 黄片大片在线免费观看| 国产又黄又爽又无遮挡在线| 精品不卡国产一区二区三区| 美女大奶头视频| 国产高清有码在线观看视频| 9191精品国产免费久久| 女生性感内裤真人,穿戴方法视频| 高清在线国产一区| 午夜a级毛片| 色视频www国产| 欧美一区二区精品小视频在线| 精品久久久久久久末码| 国产精品99久久久久久久久| 九九久久精品国产亚洲av麻豆| 身体一侧抽搐| 女同久久另类99精品国产91| 日韩精品青青久久久久久| 亚洲aⅴ乱码一区二区在线播放| 欧洲精品卡2卡3卡4卡5卡区| 黄色日韩在线| 搡老岳熟女国产| 欧美一区二区亚洲| 国产视频内射| 久久久久久久久中文| 伊人久久大香线蕉亚洲五| 久久精品综合一区二区三区| 国产视频内射| 日日摸夜夜添夜夜添小说| 欧美+亚洲+日韩+国产| 90打野战视频偷拍视频| aaaaa片日本免费| 亚洲国产欧美网| 在线视频色国产色| 亚洲五月天丁香| 成人一区二区视频在线观看| 亚洲av成人不卡在线观看播放网| 亚洲成人久久性| 熟女电影av网| 午夜福利欧美成人| 少妇裸体淫交视频免费看高清| 97超级碰碰碰精品色视频在线观看| 亚洲最大成人中文| 免费看美女性在线毛片视频| 欧美日韩国产亚洲二区| 国内毛片毛片毛片毛片毛片| 香蕉久久夜色| 日韩欧美免费精品| 国产亚洲精品一区二区www| 国产一区二区激情短视频| 精品一区二区三区av网在线观看| 美女高潮喷水抽搐中文字幕| 淫妇啪啪啪对白视频| www.色视频.com| 亚洲av电影不卡..在线观看| 少妇丰满av| 又爽又黄无遮挡网站| 99国产极品粉嫩在线观看| 操出白浆在线播放| 久久精品人妻少妇| 国产成人av教育| 12—13女人毛片做爰片一| 国产亚洲av嫩草精品影院| av天堂中文字幕网| 最后的刺客免费高清国语| 激情在线观看视频在线高清| 免费人成视频x8x8入口观看| 婷婷亚洲欧美| 国产精品 国内视频| 欧美激情久久久久久爽电影| 99国产精品一区二区蜜桃av| 欧美日韩精品网址| 最近最新免费中文字幕在线| 日本成人三级电影网站| 高清毛片免费观看视频网站| 亚洲国产高清在线一区二区三| 中亚洲国语对白在线视频| 国产毛片a区久久久久| 观看美女的网站| 91九色精品人成在线观看| 露出奶头的视频| 亚洲av中文字字幕乱码综合| 一个人免费在线观看的高清视频| 欧美日本视频| 深爱激情五月婷婷| 免费看光身美女| 在线观看美女被高潮喷水网站 | 国产av麻豆久久久久久久| 一二三四社区在线视频社区8| 日韩欧美国产在线观看| 亚洲av美国av| 精品国内亚洲2022精品成人| 国产精品日韩av在线免费观看| 亚洲欧美日韩高清在线视频| 亚洲精品粉嫩美女一区| 一区二区三区国产精品乱码| 噜噜噜噜噜久久久久久91| 午夜日韩欧美国产| 极品教师在线免费播放| 伊人久久精品亚洲午夜| 丝袜美腿在线中文| av福利片在线观看| 男人的好看免费观看在线视频| 在线a可以看的网站| 欧美成人a在线观看| 婷婷精品国产亚洲av| 欧美黑人巨大hd| 亚洲av二区三区四区| 最近最新免费中文字幕在线| 全区人妻精品视频| 国产精品嫩草影院av在线观看 | 免费人成在线观看视频色| 少妇的逼好多水| 日韩大尺度精品在线看网址| 国产黄a三级三级三级人| 欧美最黄视频在线播放免费| 90打野战视频偷拍视频| 国产免费男女视频| 一个人观看的视频www高清免费观看| 婷婷亚洲欧美| 日韩欧美精品免费久久 | 亚洲国产欧洲综合997久久,| 国产激情偷乱视频一区二区| a级毛片a级免费在线| 亚洲午夜理论影院| 午夜亚洲福利在线播放| 白带黄色成豆腐渣| 男人的好看免费观看在线视频| 成年女人看的毛片在线观看| 亚洲国产色片| av中文乱码字幕在线| aaaaa片日本免费| 午夜福利18| 色视频www国产| 美女免费视频网站| 黄色丝袜av网址大全| 少妇人妻精品综合一区二区 | 亚洲专区国产一区二区| 亚洲最大成人手机在线| 午夜福利在线观看免费完整高清在 | 狂野欧美白嫩少妇大欣赏| 伊人久久大香线蕉亚洲五| 一级黄片播放器| 国产野战对白在线观看| 欧美另类亚洲清纯唯美| 精品国产美女av久久久久小说| 国产 一区 欧美 日韩| 黄片小视频在线播放| 夜夜躁狠狠躁天天躁| 床上黄色一级片| 亚洲中文日韩欧美视频| 亚洲av第一区精品v没综合| 手机成人av网站| 亚洲五月天丁香| 在线天堂最新版资源| 精品欧美国产一区二区三| 12—13女人毛片做爰片一| 日韩亚洲欧美综合| 久久久久亚洲av毛片大全| 欧美黄色淫秽网站| 在线观看一区二区三区| 国产色爽女视频免费观看| 亚洲精品美女久久久久99蜜臀| 免费一级毛片在线播放高清视频| 叶爱在线成人免费视频播放| 成年版毛片免费区| 欧美成狂野欧美在线观看| 精品电影一区二区在线| www国产在线视频色| 少妇丰满av| 中文字幕人妻熟人妻熟丝袜美 | 成人永久免费在线观看视频| 天天躁日日操中文字幕| 熟女电影av网| 国产麻豆成人av免费视频| 日本五十路高清| 欧美一级毛片孕妇| 免费电影在线观看免费观看| 精品久久久久久久人妻蜜臀av| 偷拍熟女少妇极品色| www日本在线高清视频| 中国美女看黄片| 国内精品美女久久久久久| 午夜激情福利司机影院| 国产不卡一卡二| 国产在视频线在精品| 国产国拍精品亚洲av在线观看 | 一二三四社区在线视频社区8| 麻豆一二三区av精品| 欧美日韩乱码在线| 女同久久另类99精品国产91| 亚洲成人久久爱视频| 国产午夜精品久久久久久一区二区三区 | 色老头精品视频在线观看| 免费看十八禁软件| 日本免费一区二区三区高清不卡| 成年免费大片在线观看| 久久99热这里只有精品18| 少妇熟女aⅴ在线视频| 国产精品98久久久久久宅男小说| 极品教师在线免费播放| 99久久精品一区二区三区| 国产成年人精品一区二区| 亚洲男人的天堂狠狠| 亚洲美女黄片视频| 桃色一区二区三区在线观看| 黄色女人牲交| 成人欧美大片| 51午夜福利影视在线观看| 国产精品98久久久久久宅男小说| av女优亚洲男人天堂| 午夜影院日韩av| 亚洲久久久久久中文字幕| 国产成年人精品一区二区| 久久久久久久久中文| 成人午夜高清在线视频| 精品国产亚洲在线| 国产视频内射| 他把我摸到了高潮在线观看| 麻豆一二三区av精品| 亚洲欧美一区二区三区黑人| 日韩中文字幕欧美一区二区| 99热6这里只有精品| 欧美bdsm另类| 高潮久久久久久久久久久不卡| 午夜福利免费观看在线| 国产黄片美女视频| av女优亚洲男人天堂| 在线视频色国产色| 两个人看的免费小视频| 俄罗斯特黄特色一大片| 女同久久另类99精品国产91| 在线免费观看不下载黄p国产 | 99视频精品全部免费 在线| 别揉我奶头~嗯~啊~动态视频| 国产av在哪里看| 国产精品 欧美亚洲| 国产成人福利小说| 亚洲国产精品合色在线| 国产欧美日韩一区二区精品| 欧美av亚洲av综合av国产av| 精品人妻一区二区三区麻豆 | 精品国内亚洲2022精品成人| 性色av乱码一区二区三区2| 日本成人三级电影网站| 18禁裸乳无遮挡免费网站照片| 又紧又爽又黄一区二区| 日本a在线网址| 可以在线观看的亚洲视频| 男女那种视频在线观看| 久久香蕉国产精品| 窝窝影院91人妻| 国产亚洲欧美在线一区二区| 可以在线观看毛片的网站| 亚洲精品粉嫩美女一区| 99久久精品一区二区三区| 亚洲精品色激情综合| 精品人妻一区二区三区麻豆 | 国产精品一及| 99久久成人亚洲精品观看| 九九久久精品国产亚洲av麻豆| 亚洲人成网站在线播| 99久久精品热视频| 精品久久久久久久久久久久久| 免费电影在线观看免费观看| 淫妇啪啪啪对白视频| 亚洲精品国产精品久久久不卡| 波多野结衣巨乳人妻| 最近视频中文字幕2019在线8| a级一级毛片免费在线观看| 国产精品三级大全| netflix在线观看网站| 亚洲熟妇熟女久久| 成人永久免费在线观看视频| 中亚洲国语对白在线视频| 每晚都被弄得嗷嗷叫到高潮| 很黄的视频免费| 在线视频色国产色| 日本在线视频免费播放| 88av欧美| 国产精品野战在线观看| 免费观看人在逋| 变态另类成人亚洲欧美熟女| 日韩高清综合在线| 精品国产超薄肉色丝袜足j| 国产亚洲精品久久久com| 不卡一级毛片| 亚洲成av人片免费观看| 午夜两性在线视频| 亚洲人与动物交配视频| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆 | 丰满人妻一区二区三区视频av | 精华霜和精华液先用哪个| 精品久久久久久成人av| 国产精品,欧美在线| 免费观看精品视频网站| 久久精品影院6| 香蕉久久夜色| 亚洲精品影视一区二区三区av| 国产视频一区二区在线看| 国产淫片久久久久久久久 | 麻豆久久精品国产亚洲av| 99精品久久久久人妻精品| av中文乱码字幕在线| 中出人妻视频一区二区| 色吧在线观看| 亚洲中文日韩欧美视频| 99国产精品一区二区蜜桃av| 欧美日本视频| 一二三四社区在线视频社区8| 少妇的丰满在线观看| 日韩免费av在线播放| 成人性生交大片免费视频hd| 国产精品久久久久久亚洲av鲁大| 老司机福利观看| 国产精品av视频在线免费观看| 国产成人av激情在线播放| АⅤ资源中文在线天堂| 真实男女啪啪啪动态图| 欧美一区二区国产精品久久精品| 婷婷亚洲欧美| 日本一二三区视频观看| 欧美+亚洲+日韩+国产| 欧美日韩亚洲国产一区二区在线观看| 国产欧美日韩一区二区三| 免费一级毛片在线播放高清视频| 亚洲性夜色夜夜综合| 亚洲av二区三区四区| 亚洲成人久久爱视频| 欧美bdsm另类| 最新美女视频免费是黄的| 国产精品日韩av在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一区二区三区四区久久| 又黄又粗又硬又大视频| 男女视频在线观看网站免费| 麻豆久久精品国产亚洲av| 久久久久性生活片| 男人舔女人下体高潮全视频| 观看免费一级毛片| 丰满乱子伦码专区| 久久久久国产精品人妻aⅴ院| 午夜老司机福利剧场| 国产精品综合久久久久久久免费| 久久久国产成人精品二区| 免费电影在线观看免费观看| 国产欧美日韩精品一区二区| 深爱激情五月婷婷| 国产伦一二天堂av在线观看| 亚洲精华国产精华精| 男人舔奶头视频| 国产成人啪精品午夜网站| www.999成人在线观看| 无限看片的www在线观看| 亚洲五月天丁香| 在线观看av片永久免费下载|