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

    地面沉降-回彈及地下水位波動(dòng)的InSAR長(zhǎng)時(shí)序監(jiān)測(cè)——以德州市為例

    2014-09-26 02:23:26葛大慶殷躍平郭小方
    自然資源遙感 2014年1期
    關(guān)鍵詞:德州深層基線

    葛大慶,殷躍平,王 艷,張 玲,郭小方,王 毅

    (1.中國(guó)地質(zhì)大學(xué)(北京)水資源與環(huán)境學(xué)院,北京 100083;2.中國(guó)國(guó)土資源航空物探遙感中心,北京 100083;3.中國(guó)地質(zhì)環(huán)境監(jiān)測(cè)院,北京 100081)

    0 引言

    長(zhǎng)期過(guò)量開(kāi)采地下水,特別是深層地下水,是目前我國(guó)一些平原和三角洲地區(qū)產(chǎn)生地面沉降的主要原因[1]。地下水位的持續(xù)下降改變了原有的補(bǔ)徑排特征,形成了區(qū)域性地下水位降落漏斗,地面沉降隨之產(chǎn)生并擴(kuò)大。我國(guó)已形成了以華北平原、長(zhǎng)江三角洲及汾渭斷陷盆地為主的三大地面沉降嚴(yán)重區(qū),華北平原已出現(xiàn)近20個(gè)大型復(fù)合水位降落漏斗,地面沉降呈現(xiàn)出區(qū)域性且連片分布的特征[2-3]。

    雷達(dá)干涉測(cè)量(InSAR)技術(shù)為區(qū)域性地面沉降的快速準(zhǔn)確和連續(xù)監(jiān)測(cè)提供了重要手段。以永久性散射體干涉測(cè)量(permanent scatterer InSAR,PSIn-SAR)[4]為代表的相干目標(biāo)(coherent target,CT)時(shí)間序列分析技術(shù)克服了差分InSAR技術(shù)所面臨的失相干和大氣相位延遲等問(wèn)題,通過(guò)對(duì)相干目標(biāo)的差分干涉相位序列進(jìn)行時(shí)序分析,逐個(gè)估計(jì)和分離地形相位、形變相位、大氣相位和噪聲等信息,進(jìn)而解算相干目標(biāo)的形變速率和累積形變量。與之類似的還有短基線集(small baseline subset,SBAS)[5]以及干涉測(cè)量點(diǎn)目標(biāo)分析(interferometric point target analysis,IPTA)[6]等技術(shù),這些技術(shù)已在地面沉降、滑坡、開(kāi)采沉陷、地震和冰川滑移等多尺度緩慢地表形變研究中得以發(fā)展和應(yīng)用[7-9]。

    德州市位于山東省北部,是華北平原地面沉降和深層地下水位降落漏斗極為典型的地區(qū)之一。其地下水位降落漏斗在1978年前已經(jīng)形成,是目前華北平原水位降落深度最大的漏斗,與衡水、滄州等漏斗連成一片。該區(qū)的地面沉降在1989年以前就已經(jīng)出現(xiàn),至2010年累積沉降量超過(guò)1 m。與地面沉降對(duì)應(yīng)的是,德州市地下水位呈現(xiàn)出多年連續(xù)下降和年內(nèi)季節(jié)性波動(dòng)的特征[10]。由于地面沉降與地下水位下降相關(guān),因此在德州市深、淺層地下水位均呈現(xiàn)出季節(jié)性波動(dòng)的情況下,地面沉降場(chǎng)的時(shí)空演變特征值得研究,這包括:沉降場(chǎng)的空間變化特征,即沉降中心的移動(dòng)及其變化幅度;沉降場(chǎng)的時(shí)間變化特征,即季節(jié)性波動(dòng),包括是否會(huì)出現(xiàn)地面回彈,回彈幅度多大,地面沉降與地下水位波動(dòng)是否同步等方面。針對(duì)上述問(wèn)題,本文以相干目標(biāo)InSAR時(shí)序分析方法為主,利用2004年1月—2010年10月近7 a的ENVISAT衛(wèi)星數(shù)據(jù)研究分析德州地區(qū)地面沉降-回彈的動(dòng)態(tài)變化特征及其與地下水位波動(dòng)的對(duì)應(yīng)關(guān)系,以揭示該地區(qū)地面沉降季節(jié)性變化的產(chǎn)生原因和誘發(fā)機(jī)制。

    1 相干目標(biāo)InSAR時(shí)序分析方法

    InSAR時(shí)序分析主要有PSInSAR和SBAS兩類方法。根據(jù)干涉像對(duì)的組合模式,前者以單一主影像構(gòu)成差分干涉圖序列,后者則利用短基線原則以多個(gè)主影像構(gòu)成若干個(gè)干涉紋圖子集。單一主影像的干涉紋圖序列處理主要顧及大氣相位和非線性形變的統(tǒng)計(jì)模型,有利于對(duì)大氣相位的準(zhǔn)確估計(jì),但對(duì)數(shù)據(jù)量要求過(guò)高,且受相干性影響而對(duì)大變形的監(jiān)測(cè)能力不足;多個(gè)主影像則降低了對(duì)數(shù)據(jù)量的要求,將時(shí)空基線小于一定閾值的干涉像對(duì)組合,生成干涉紋圖序列,增加對(duì)同一信號(hào)的采樣密度,有利于準(zhǔn)確求解形變速率,在數(shù)據(jù)量較小的情況下也可以實(shí)現(xiàn)。在選擇相干目標(biāo)的方式上,前者以點(diǎn)目標(biāo)為主,后者兼顧點(diǎn)目標(biāo)和分布式目標(biāo)[5]。本文集成了2種方法的優(yōu)點(diǎn),以短基線為準(zhǔn)則構(gòu)成差分干涉相位圖,利用點(diǎn)目標(biāo)識(shí)別算法提取相干目標(biāo)。以Delanay三角網(wǎng)連接相鄰的相干像元,用二維周期圖估計(jì)點(diǎn)間形變速率和高程誤差改正。以此為基礎(chǔ),根據(jù)大氣相位特征,利用奇異值分解算法解算空間濾波后的殘余相位,并對(duì)解算結(jié)果進(jìn)行時(shí)域?yàn)V波,以求解非線性形變,得到不同時(shí)刻的累積變形量,再利用線性回歸法求解每個(gè)相干目標(biāo)的變形速率。

    1.1 干涉相位時(shí)序分析的基本模型

    對(duì)于給定的M景雷達(dá)數(shù)據(jù),在短基線條件下生成N個(gè)差分干涉圖(一般情況下,N>M),對(duì)于其中的任一差分干涉圖k(k=1,2,…,N),其相位模型為

    顯然,待求解量d與InSAR參數(shù)不直接相關(guān),但受其他分量的影響。對(duì)于大氣分量,其本質(zhì)上是2次大氣延遲影響的疊加,在空間上具有一定的相關(guān)性[4]。大氣波動(dòng)在時(shí)間上的隨機(jī)特性決定了2次互差也為隨機(jī)量。需要重點(diǎn)區(qū)分的是大氣相位與軌道誤差,雖然二者都具有空間低頻的特性,但各自的影響范圍有所不同。基線誤差多是全局性的,而大氣影響半徑有限,對(duì)于大范圍處理,二者應(yīng)區(qū)別對(duì)待。高程誤差與垂直基線相關(guān),垂直基線是變量。變形量與干涉圖的時(shí)間基線相關(guān),時(shí)間間隔是變量,相對(duì)變形的速率與時(shí)間相關(guān)。

    1.2 地表形變參數(shù)的提取算法

    1.2.1 線性分量的估計(jì)

    考慮到大氣的空間相關(guān)性,如果對(duì)相鄰2點(diǎn)(小于大氣相關(guān)距離)求差,則可削弱大氣影響。相鄰相干目標(biāo)i和j差分干涉相位的互差△為

    由于相鄰目標(biāo)間的相對(duì)變形速率與時(shí)間間隔相關(guān),干涉圖的時(shí)間基線為自變量,因而可將形變量d分解為線性形變速率和非線性形變量2部分。如此,時(shí)序分析過(guò)程則轉(zhuǎn)為參數(shù)估計(jì),即

    式中:CB為與垂直基線相關(guān)的系數(shù);T為時(shí)間基線;△ε為相對(duì)高程誤差;△υ為相對(duì)形變速率;μN(yùn)L為非線性形變量;α為大氣相位;n為噪聲。

    對(duì)于纏繞相位,其周期數(shù)由相對(duì)高程誤差和線性形變速率共同決定,采用二維參數(shù)估計(jì)方法可實(shí)現(xiàn)對(duì)相位周期的估算。由此,構(gòu)建目標(biāo)函數(shù)

    并將其從相位互差中減去,得到殘余相位,即

    顯然,在準(zhǔn)確估計(jì)目標(biāo)函數(shù)的參數(shù)△ε和△υ時(shí),殘余相位將最小化。

    由于相鄰2點(diǎn)存在多個(gè)干涉圖,即存在N維B和T,因而,對(duì)△ε和△v的估計(jì)使模型相關(guān)系數(shù)最大化,即

    在纏繞相位條件下,△ε和△v為周期函數(shù)的二維頻率,可利用二維周期圖估計(jì)使模型相關(guān)系數(shù)最大化;在解纏相位條件下,則轉(zhuǎn)換為二維線性函數(shù),利用參數(shù)估計(jì)方法可實(shí)現(xiàn)對(duì)△ε和△v的求解。

    在相鄰2點(diǎn)解算時(shí)需對(duì)所有相干目標(biāo)建立連接關(guān)系,整體求解形變場(chǎng)的速率和高程誤差改正量。這一目標(biāo)的實(shí)現(xiàn)可通過(guò)建立Delanay三角網(wǎng)的方法完成,也可以利用冗余網(wǎng)構(gòu)建更為復(fù)雜的連接關(guān)系,強(qiáng)化對(duì)待解算方程組的約束。利用鄰近法則將所有距離滿足大氣相關(guān)距離的相干目標(biāo)連接起來(lái),在求解完成相鄰點(diǎn)間的互差后,通過(guò)最小二乘或加權(quán)平均的方法求解每個(gè)目標(biāo)相對(duì)于參考點(diǎn)的變形速率,即總體形變速率場(chǎng)。

    1.2.2 非線性分量的估計(jì)

    對(duì)于具有顯著非線性形變過(guò)程,仍需對(duì)殘余相位進(jìn)行更為復(fù)雜的處理,以提取非線性形變量。從差分相位中去除高程誤差估計(jì)值和線性速率后,殘余相位為

    由于已解算出整周相位,因而殘余相位為相位主值,其大小

    這里,首先要實(shí)現(xiàn)對(duì)單個(gè)差分干涉圖中殘余相位的濾波處理。由于大氣表現(xiàn)出空間低頻特性,而非線性變形的空間變化范圍較小,相對(duì)大氣而言表現(xiàn)為高通特性。因而,對(duì)殘余相位圖進(jìn)行空間低通濾波可以進(jìn)一步弱化大氣影響。在短基線集條件下求解的非線性形變量,不同于PSInSAR經(jīng)典處理策略,不直接對(duì)大氣進(jìn)行估計(jì),而是通過(guò)濾波減弱大氣的影響。殘余相位經(jīng)空間高通濾波后,大氣分量已經(jīng)減弱。此時(shí),進(jìn)一步的處理是求解與雷達(dá)數(shù)據(jù)對(duì)應(yīng)的不同時(shí)刻的非線性變形量,包含噪聲等影響。根據(jù)主輔影像的關(guān)系,可將解纏后的殘余相位分解為

    式中:p和q分別表示生成第k(k=1,…,N)景差分干涉圖的主輔影像的獲取時(shí)間。由于采用短基線原則,在求解每個(gè)時(shí)刻對(duì)應(yīng)的殘余相位時(shí)出現(xiàn)秩虧方程組,解決這一問(wèn)題的辦法就是奇異值分解法(singular value decomposition,SVD)[3]。在此基礎(chǔ)上,對(duì)M景雷達(dá)數(shù)據(jù)的殘余相位進(jìn)行時(shí)域低通濾波處理,以提取最終的非線性形變量,可表示為

    在求解出非線性形變量后,每個(gè)相干目標(biāo)對(duì)應(yīng)的相變序列為

    式中vest為線性形變速率。

    1.2.3 平均形變速率的重新估計(jì)

    對(duì)于顯著的非線性形變過(guò)程,利用前述二維線性模型估計(jì)出來(lái)的線性速率往往偏小。在恢復(fù)整個(gè)形變序列的基礎(chǔ)上,可根據(jù)對(duì)應(yīng)的每個(gè)目標(biāo)的形變序列,利用最小二乘法重新估計(jì)形變速率,即

    式中d和T分別為相對(duì)起始時(shí)刻的形變量和時(shí)間間隔。

    2 雷達(dá)數(shù)據(jù)及InSAR時(shí)序分析處理

    研究中共獲取2004年1月—2010年10月近7 a的ENVISAT衛(wèi)星ASAR數(shù)據(jù)。該數(shù)據(jù)覆蓋了德州地區(qū)主要沉降區(qū),總數(shù)據(jù)量為46景(Track-447,F(xiàn)rame-2853),其基線分布如圖1所示。

    圖1 差分干涉圖序列的短基線關(guān)系Fig.1 Small baseline of the differential interferograms series

    需要說(shuō)明的是,利用幅度離散指數(shù)[1]、相干系數(shù)[5]以及子視相關(guān)[6]均可以很好地識(shí)別出相干目標(biāo)候選點(diǎn)。對(duì)于候選點(diǎn),只有在滿足模型相關(guān)函數(shù)最大化約束條件下才能稱為時(shí)序分析中的相干目標(biāo)。實(shí)際上,一些目標(biāo)后向散射特性穩(wěn)定,但非線性變化強(qiáng)烈,不滿足該條件。因而,可通過(guò)降低相關(guān)系數(shù)的閾值,將其納入到形變解算網(wǎng)絡(luò)中,以增加監(jiān)測(cè)點(diǎn)的密度。

    時(shí)空基線閾值的設(shè)置取決于SAR系統(tǒng)參數(shù)。空間基線主要考慮雷達(dá)系統(tǒng)臨界基線,而時(shí)間基線則應(yīng)根據(jù)變形幅度大小,短時(shí)間基線相當(dāng)于降低了形變場(chǎng)的梯度,便于相位解纏和變形參數(shù)估計(jì)。因而,本研究中設(shè)定的短基線條件為:時(shí)間基線<450 d,空間基線<300 m。

    圍繞德州地區(qū)地面沉降和地下水位漏斗的基本分布特征,分別針對(duì)區(qū)域性大范圍沉降和重點(diǎn)沉降中心進(jìn)行了2類處理:利用2007—2010年4 a間的數(shù)據(jù)提取了德州地區(qū)100 km×100 km平均地面沉降速率;分析全區(qū)地面沉降和地下水位變化關(guān)系,并確定時(shí)序分析過(guò)程中穩(wěn)定的地面參考點(diǎn)。在此基礎(chǔ)上,對(duì)德州城區(qū)內(nèi)20 km×20 km范圍的主要沉降中心的7 a數(shù)據(jù)進(jìn)行了時(shí)序分析處理,提取每個(gè)相干目標(biāo)的沉降率和累積沉降量序列,用以詳細(xì)分析多年和年度內(nèi)沉降中心的時(shí)空變化特征。

    3 區(qū)域性地面沉降與地下水位變化

    3.1 區(qū)域地面沉降分布

    圖2示出德州地區(qū)約50 km×60 km范圍內(nèi)2007年初—2010年底4 a間平均地面沉降速率。

    圖2 德州地區(qū)2007—2010年CTInSAR監(jiān)測(cè)平均地面沉降率Fig.2 Average subsidence velocity of Dezhou area from 2007 to 2010 derived by CTInSAR data

    從區(qū)域分布上來(lái)看,德州市區(qū)沉降較為顯著,沉降區(qū)向西北方向延伸,與河北省景縣和故城相連,構(gòu)成連片沉降區(qū)。區(qū)內(nèi)景縣沉降向西北方向擴(kuò)展,是目前衡水—德州沉降區(qū)的重要組成部分,沉降速率均達(dá)到40~50 mm/a。德州沉降區(qū)向東南方向,分布著3個(gè)明顯的沉降中心,分別為陵縣、臨邑縣和平原縣,其中以陵縣最為突出,中心位于縣城西北菜園村,最大沉降速率達(dá)75 mm/a。臨邑縣沉降中心位于臨盤(pán)采油廠,中心處沉降速率達(dá)50~60 mm/a。

    3.2 與深層地下水位漏斗的比較分析

    InSAR監(jiān)測(cè)結(jié)果準(zhǔn)確揭示了德州地區(qū)區(qū)域性沉降和單個(gè)漏斗的分布。連續(xù)多年的地下水位變化監(jiān)測(cè)結(jié)果表明,德州市深層地下水位降落漏斗是目前魯北地區(qū)最大的漏斗,與周邊的小漏斗相連接,引起區(qū)域性地面沉降[8-9]。深層地下水開(kāi)采(第Ⅳ,Ⅴ含水層組)是該地區(qū)地面沉降的主要原因。以臨邑為例,該沉降中心深層地下水降落漏斗的-30 m等水壓線與德城區(qū)深層漏斗相連,已成為德州漏斗的一部分;-40 m等水壓線包圍臨邑縣城和臨盤(pán)鎮(zhèn),圈閉面積達(dá)192 km2。漏斗中心在臨盤(pán)采油廠一帶,水位埋深超過(guò)59.70 m,水位標(biāo)高-41.82 m左右。圖2中的沉降中心準(zhǔn)確地反映了目前該漏斗的分布位置和影響范圍。由于本區(qū)地下水水位動(dòng)態(tài)屬于較穩(wěn)定的連續(xù)開(kāi)采消耗型,年水位變化趨勢(shì)與開(kāi)采量大小密切相關(guān)。臨邑深層地下水開(kāi)采量達(dá)3.56×104m3/a。其中,第Ⅳ和第Ⅴ含水層組地下水開(kāi)采量約占90%,地下水位多年平均降速2.0 m/a,從而使臨邑深層地下水降落漏斗迅速擴(kuò)展,成為較大的區(qū)域性降落漏斗。

    4 德城區(qū)沉降中心的沉降與回彈

    4.1 德城區(qū)主要沉降中心

    2010年之前連續(xù)7 a的InSAR時(shí)序分析結(jié)果(圖3)表明,德城區(qū)有3個(gè)典型沉降中心,分別位于德城區(qū)陳莊(德棉一廠)、長(zhǎng)莊(華源紡織廠)和肖何莊。其中陳莊沉降中心呈南北向分布,長(zhǎng)約2.2 km,寬1.5 km,中心最大沉降速率達(dá)到55 mm/a;長(zhǎng)莊沉降中心呈北東向,長(zhǎng)為2 km,寬為1.5 km,中心最大速率也超過(guò)50 mm/a。連片分布沉降出現(xiàn)在德城區(qū)東部的付莊和宋官屯鎮(zhèn),這一地區(qū)為近年來(lái)快速發(fā)展的新區(qū)。根據(jù)1991—2010年水準(zhǔn)測(cè)量結(jié)果[11],德城區(qū)沉降中心累計(jì)沉降量為-1 186.9~-636.9 mm,位于德棉一廠(即陳莊漏斗)位置,中心多年平均速率為59.35 mm/a。

    圖3 德城區(qū)2004—2010年CTInSAR監(jiān)測(cè)平均地面沉降率Fig.3 Average subsidence velocity of Decheng district from 2004 to 2010 derived by CTInSAR data

    其中,2005—2006年沉降量為-56 mm;2006—2007年沉降量為-89.0 mm;2007—2010年總沉降量為93.2 mm,年平均沉降量為31.1 mm,與 InSAR監(jiān)測(cè)反映的德城區(qū)陳莊沉降中心的變化完全一致。

    4.2 地面沉降與回彈的時(shí)間效應(yīng)

    由2004—2010年InSAR監(jiān)測(cè)結(jié)果發(fā)現(xiàn):德城區(qū)多個(gè)沉降中心年內(nèi)均表現(xiàn)出季節(jié)性沉降與回彈特征(圖4)。陳莊、長(zhǎng)莊、肖何莊以及德棉二廠4個(gè)典型沉降中心均出現(xiàn)回彈,回彈最顯著的位于陳莊、長(zhǎng)莊和德棉二廠沉降中心。其中陳莊的沉降和回彈幅度均為最大,對(duì)應(yīng)的年份為2007年,2—8月間沉降幅度超過(guò)80 mm,至翌年2月回彈15~20 mm。除德棉二廠外,其他3個(gè)中心2008年沉降幅度均有所減緩,范圍縮小,2009年繼續(xù)減緩,而市區(qū)東部在8月份之后沉降加快,2010年繼續(xù)增加。

    圖4 德城區(qū)典型沉降中心的季節(jié)性沉降與回彈Fig.4 Seasonal subsidence and rebound of the major center in Decheng district

    圖4 所示5個(gè)年度的變化特征表明,以陳莊為代表的沉降中心具有明顯的季節(jié)性變化特征,表現(xiàn)為每年的2—3月開(kāi)始快速沉降,至8—9月達(dá)到最大,此后逐步回彈,至翌年2—3月重新開(kāi)始沉降,沉降與回彈之比約為4∶1。

    5 地面沉降及回彈與地下水位變化的關(guān)系

    5.1 持續(xù)沉降與深層地下水位的關(guān)系

    德州地區(qū)地下水系統(tǒng)分為4個(gè)孔隙含水巖組:第Ⅰ含水巖組潛水層,是主要開(kāi)采層;第Ⅱ含水巖組屬于咸水層,基本不開(kāi)采;第Ⅲ含水巖組,底板埋深450~500 m;第Ⅳ含水巖組,底板埋深800~950 m。第Ⅲ,Ⅳ含水巖組屬于承壓含水層,水頭埋深70~100 m,水質(zhì)普遍較好,是城鎮(zhèn)居民生活和工業(yè)用水的主要開(kāi)采層,德州地下水位降落漏斗屬于該層位,其中心水位(位于陳莊國(guó)棉一廠)埋深超過(guò)120 m[10]。根據(jù)德州深層地下開(kāi)采引發(fā)地面沉降變化閾值的研究[11]:持續(xù)開(kāi)采深層地下水是地面沉降的主要原因,地面沉降與深層水位降落呈現(xiàn)出顯著的線性關(guān)系,深層地下水位每下降1 m,產(chǎn)生的地面沉降約為17.5 mm。該經(jīng)驗(yàn)值與近4 a來(lái)地面測(cè)量結(jié)果一致,速率為30~40 mm/a。

    圖5所示為德城區(qū)202長(zhǎng)觀孔(位于陳莊漏斗區(qū))1996—2010年的深層水位變化情況,表現(xiàn)出持續(xù)下降和年內(nèi)波動(dòng)的特點(diǎn)。德城區(qū)淺層地下水位變化曲線如圖6所示。

    圖5 德城區(qū)202長(zhǎng)觀孔深層地下水位變化曲線Fig.5 Dynamic curve of deep groundwater level of No.202 hole in Decheng district

    圖6 德城區(qū)淺層地下水位變化曲線Fig.6 Dynamic curve of shallow groundwaterlevel in Decheng district

    5.2 季節(jié)性沉降-回彈與地下水位波動(dòng)的關(guān)系

    以陳莊漏斗為例,其動(dòng)態(tài)變化(圖3和圖4)與深、淺層地下水位的季節(jié)性波動(dòng)密切相關(guān)。每年1—5月降水少,農(nóng)業(yè)用水量大,淺層地下水位急劇下降;7—8月雨量增多、河道水量補(bǔ)給和開(kāi)采量相對(duì)減少,地下水埋深逐漸回升,至年底達(dá)到次年最高點(diǎn)。深層承壓水位變化主要取決于開(kāi)采量,年水位變化與開(kāi)采量關(guān)系密切。每年3—7月,深層地下水連續(xù)大量開(kāi)采,漏斗內(nèi)壓力水頭急劇下降;8月—翌年2月由于開(kāi)采量減小,加之大量引入黃河水,水位有所回升,回升幅度一般小于2.0 m。

    開(kāi)采量決定了地面沉降的范圍和強(qiáng)度。目前,在深層水位降落漏斗中心區(qū)(陳莊),開(kāi)采井密度為334眼/km2。2007年以前,由于深層地下水開(kāi)采量大,漏斗內(nèi)壓力水逐年下降,年下降幅度較大,一般2~4 m;2008年以后,由于采取封井措施,深層地下水開(kāi)采量減少,水位呈現(xiàn)回升趨勢(shì)[11]。

    陳莊、長(zhǎng)莊等中心的沉降與回彈的波動(dòng)與上述變化十分吻合。圖7分別給出了陳莊(國(guó)棉一廠)、長(zhǎng)莊(華源紡織廠)、國(guó)棉二廠、肖何莊和宋官屯5個(gè)沉降中心連續(xù)7 a的沉降過(guò)程(監(jiān)測(cè)日期為每年的1月30日),均表現(xiàn)出季節(jié)性的非線性波動(dòng)(紅色線表示初次估計(jì)沉降速率)。以陳莊為例,其年度內(nèi)沉降與回彈與近年深層水位波動(dòng)相對(duì)應(yīng),沉降速率在2008年后出現(xiàn)回落,小于2007年之前。

    圖7 德城區(qū)典型沉降中心2004—2010年地面沉降及深層水位變化序列Fig.7 Subsidence history of coherent target in the typical subsidence center of Decheng district

    從地下水天然補(bǔ)給而言,春旱夏澇的氣候特征決定了德州地區(qū)春季地下水補(bǔ)給量的不足[11],而同時(shí)用水量又進(jìn)一步增大,特別是陳莊、長(zhǎng)莊和國(guó)棉二廠3個(gè)沉降中心,均為棉紡、化工等耗水企業(yè),需要大量采水,而汛期澇季時(shí),地表徑流補(bǔ)給及時(shí),地下水開(kāi)采總量減小,水位出現(xiàn)回升,地面沉降也隨之減緩,并在部分時(shí)段內(nèi)出現(xiàn)回彈的現(xiàn)象。

    6 結(jié)論

    1)2004—2010年間的連續(xù)InSAR監(jiān)測(cè)有效揭示了德州地區(qū)區(qū)域性地面沉降和典型漏斗季節(jié)性地面沉降-回彈的發(fā)展與變化特征。德州市及其周邊主要地下水降落漏斗與地面沉降中心相互一致,證明了深層地下水開(kāi)采是本地區(qū)地面沉降的主要原因。

    2)德城區(qū)陳莊、長(zhǎng)莊等典型漏斗的季節(jié)性沉降-回彈與該地區(qū)深淺層地下水位波動(dòng)相關(guān),地面沉降與地下水位變化表現(xiàn)出良好的相關(guān)關(guān)系,變化時(shí)間基本同步,滯后期約1個(gè)月。地面沉降漏斗呈現(xiàn)出3—8月快速沉降,沉降幅度達(dá)50~80 mm,8月—翌年2月地面回彈,回彈幅度達(dá)10~20 mm的變化特點(diǎn)。

    志謝:感謝山東省魯北工程地質(zhì)勘察院鄒祖光、陳松、楊亞賓以及河北水文地質(zhì)工程地質(zhì)四隊(duì)杜興明、田小偉等專家的指導(dǎo)和深入討論。

    [1]殷躍平,張作辰,張開(kāi)軍.我國(guó)地面沉降現(xiàn)狀及防治對(duì)策研究[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2005,16(2):1-8.Yin Y P,Zhang Z C,Zhang K J.Land subsidence and countermeasures for its prevention in China[J].The Chinese Jounal of Geological Hazard and Control,2005,16(2):1- 8.

    [2]葛大慶,郭小方,張 玲,等.華北平原地面沉降區(qū)InSAR監(jiān)測(cè)[R].北京:中國(guó)國(guó)土資源航空物探遙感中心,2010.Ge D Q,Guo X F,Zhang L,et al.Subsidence monitoring with SAR interferometry in the North China Plain[R].Beijing:China Aero Geophysical Surveying and Remote Sensing Center for Land and Resources,2010.

    [3]王 艷,葛大慶,郭小方,等.長(zhǎng)江三角洲地區(qū)地面沉降InSAR監(jiān)測(cè)[R].北京:中國(guó)國(guó)土資源航空物探遙感中心,2010.Wang Y,Ge D Q,Guo X F,et al.Subsidence monitoring with SAR interferometry in the Yangtze River Delta area[R].Beijing:China Aero Geophysical Surveying and Remote Sensing Center for Land and Resources,2010.

    [4]Ferretti A,Prati C,Rocca F.Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry[J].IEEE Trans Geosci Remote Sensing,2001,38(5):2202-2212.

    [5]Mora O,Mallorqui J J,Broquetas A.Linear and nonlinear terrain deformation maps from a reduced set of interferometric SAR images[J].IEEE Trans Geosci Remote Sensing,2003,41(10):2243-2253.

    [6]Charles W,Wegmuller U,Andreas W,et al.Interferometric point target analysis with JERS-1 L-band SAR data[C]//Proceedings of the IEEE Transactfions on Geoscience and Remote Sensing.Toulouse,F(xiàn)rance,2003:4359-4361.

    [7]Hoffmann J,Zebker H A,Galloway D L,et al.Seasonal subsidence and rebound in Las Vegas Valley,Nevada,observed by synethetic aperture Radar interferometry[J].Water Resources Research,2001,37(6):1551-1566.

    [8]王 艷,廖明生,李德仁,等.利用長(zhǎng)時(shí)間序列相干目標(biāo)獲取地面沉降場(chǎng)[J].地球物理學(xué)報(bào),2007,50(2):598-604.Wang Y,Liao M S,Li D R,et al.Subsidence velocity retrieveal from long term coherent targets in Radar interferometric stacks[J].Chinese Jounal Geophys,2007,50(2):598-604.

    [9]http://www.terrafirma.eu.com/documents.htm

    [10]王小剛,陳 松,王秀芹,等.德州市地面沉降成因及防治對(duì)策淺析[J].地質(zhì)災(zāi)害與環(huán)境保護(hù),2006,17(3):62-66.Wang X G,Chen S,Wang X Q,et al.A Discussion on the causes of the ground subsidence and its countermeasures in Dezhou City[J].Journal of Geological Hazards and Environment Preservation,2006,17(3):62-66.

    [11]楊麗芝.德州深層地下水位降落漏斗演變機(jī)制與可調(diào)控性研究[D].北京:中國(guó)地質(zhì)科學(xué)院,2009.Yang L Z.Form principle and controlling- adjusting research about deep ground water depression cone in Dezhou[D].Beijing:Chinese Academy of Geologieal Seienee,2009.

    [12]王勝嶺,宋 波,王德生,等.德州市臨盤(pán)采油區(qū)地面沉降監(jiān)測(cè)[J].山東國(guó)土資源,2009,25(1):25-28,32.Wang S L,Song B,Wang D S,et al.Land subsidence monitoring of Linpan oil mining area in Dezhou City[J].Land and Resources in Shandong Province,2009,25(1):25-28,32.

    [13]何國(guó)榮,馮在敏,馮克印,等.水準(zhǔn)測(cè)量網(wǎng)監(jiān)測(cè)在德州市地面沉降中的應(yīng)用[J].山東國(guó)土資源,2012,28(6):28-30.He G R,F(xiàn)eng Z M,F(xiàn)eng K Y,et al.Application of leveling surveying for subsidence monitoring in Dezhou City[J].Land and Resources in Shandong Province,2012,28(6):28-30.

    猜你喜歡
    德州深層基線
    德州大陸架石油工程技術(shù)有限公司
    鉆采工藝(2022年4期)2022-10-22 10:24:20
    適用于MAUV的變基線定位系統(tǒng)
    航天技術(shù)與甚長(zhǎng)基線陣的結(jié)合探索
    科學(xué)(2020年5期)2020-11-26 08:19:14
    德州魯源貨場(chǎng)信號(hào)聯(lián)鎖設(shè)備關(guān)鍵技術(shù)的應(yīng)用
    考慮各向異性滲流的重力壩深層抗滑穩(wěn)定分析
    SAM系統(tǒng)對(duì)TDCS數(shù)據(jù)的優(yōu)化處理與深層應(yīng)用
    一種改進(jìn)的干涉儀測(cè)向基線設(shè)計(jì)方法
    在德州,電力經(jīng)紀(jì)人幫你選電!
    能源(2016年3期)2016-12-01 05:10:48
    德州地區(qū)懸鈴木方翅網(wǎng)蝽的綜合防治措施
    對(duì)“醫(yī)患失去信任”的深層憂慮
    99国产精品一区二区蜜桃av| 全区人妻精品视频| 午夜福利18| 国产欧美日韩精品亚洲av| 国产黄色小视频在线观看| 99在线视频只有这里精品首页| 直男gayav资源| 欧美人与善性xxx| 麻豆久久精品国产亚洲av| 色综合亚洲欧美另类图片| 我要搜黄色片| 美女 人体艺术 gogo| 亚洲在线观看片| 在线观看免费视频日本深夜| 一进一出好大好爽视频| 黄色视频,在线免费观看| 国产亚洲精品久久久com| 色尼玛亚洲综合影院| 日韩,欧美,国产一区二区三区 | 欧美极品一区二区三区四区| 成人永久免费在线观看视频| 美女大奶头视频| 激情 狠狠 欧美| 亚洲最大成人手机在线| av在线老鸭窝| 亚洲人成网站在线播放欧美日韩| 高清日韩中文字幕在线| 亚洲精品国产av成人精品 | 97超视频在线观看视频| 最后的刺客免费高清国语| 午夜福利成人在线免费观看| 日韩,欧美,国产一区二区三区 | 三级国产精品欧美在线观看| 99热网站在线观看| 99riav亚洲国产免费| 99热这里只有精品一区| 一个人看的www免费观看视频| 美女大奶头视频| 俄罗斯特黄特色一大片| 91精品国产九色| 在线播放无遮挡| 久久久久久九九精品二区国产| 亚洲一级一片aⅴ在线观看| 一区二区三区四区激情视频 | 中文资源天堂在线| 大型黄色视频在线免费观看| 国产黄a三级三级三级人| 日韩高清综合在线| 天堂影院成人在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品sss在线观看| 国国产精品蜜臀av免费| 国产成人aa在线观看| 99热这里只有是精品50| 麻豆久久精品国产亚洲av| 精品乱码久久久久久99久播| 少妇猛男粗大的猛烈进出视频 | 欧美精品国产亚洲| 久久人人精品亚洲av| 3wmmmm亚洲av在线观看| 搡老岳熟女国产| 波多野结衣高清作品| 综合色丁香网| 99久国产av精品国产电影| 成人性生交大片免费视频hd| 又爽又黄无遮挡网站| 国产色婷婷99| 国产精品99久久久久久久久| 黄色一级大片看看| 日韩一本色道免费dvd| 免费观看在线日韩| 可以在线观看毛片的网站| 男女那种视频在线观看| 国产色婷婷99| 国产毛片a区久久久久| 我的老师免费观看完整版| 国产精品久久久久久亚洲av鲁大| 午夜福利在线在线| 久久这里只有精品中国| 床上黄色一级片| 国产一区二区亚洲精品在线观看| 一级毛片电影观看 | а√天堂www在线а√下载| 日本色播在线视频| 成人特级黄色片久久久久久久| 美女被艹到高潮喷水动态| 久久草成人影院| 亚洲国产色片| 日本欧美国产在线视频| 少妇猛男粗大的猛烈进出视频 | 国产综合懂色| 亚洲欧美清纯卡通| 一进一出好大好爽视频| 国产精品福利在线免费观看| 国产极品精品免费视频能看的| 精品久久久久久久人妻蜜臀av| 蜜桃久久精品国产亚洲av| 国产美女午夜福利| 午夜老司机福利剧场| 毛片女人毛片| 最近的中文字幕免费完整| 秋霞在线观看毛片| 如何舔出高潮| 午夜久久久久精精品| 亚洲中文日韩欧美视频| 免费无遮挡裸体视频| 亚洲电影在线观看av| 午夜福利高清视频| av在线观看视频网站免费| 91av网一区二区| 免费高清视频大片| 99riav亚洲国产免费| 一卡2卡三卡四卡精品乱码亚洲| 91在线观看av| 免费在线观看影片大全网站| 亚洲精品日韩av片在线观看| 亚洲国产欧美人成| 日韩亚洲欧美综合| 99久久成人亚洲精品观看| 久久精品国产亚洲av涩爱 | 亚洲国产精品sss在线观看| 乱人视频在线观看| 中文字幕精品亚洲无线码一区| 啦啦啦韩国在线观看视频| 日本黄色视频三级网站网址| 久久婷婷人人爽人人干人人爱| 日日干狠狠操夜夜爽| 国产单亲对白刺激| 亚洲性久久影院| 22中文网久久字幕| 亚洲av美国av| 欧美在线一区亚洲| 非洲黑人性xxxx精品又粗又长| 欧美国产日韩亚洲一区| 久久亚洲国产成人精品v| 国产老妇女一区| 欧美日韩在线观看h| 精品久久久久久久久av| 六月丁香七月| 国产精品伦人一区二区| av专区在线播放| 在线a可以看的网站| 深夜a级毛片| 国产一区二区三区在线臀色熟女| 美女cb高潮喷水在线观看| 99久国产av精品| 精品久久久噜噜| 亚洲无线观看免费| 身体一侧抽搐| 超碰av人人做人人爽久久| 国产v大片淫在线免费观看| 黄色配什么色好看| 国产精品免费一区二区三区在线| 可以在线观看的亚洲视频| 人人妻人人看人人澡| 亚洲真实伦在线观看| 欧美一区二区亚洲| 欧美成人精品欧美一级黄| 国产毛片a区久久久久| 久久久成人免费电影| 99riav亚洲国产免费| 欧美一级a爱片免费观看看| 亚洲成人久久性| 床上黄色一级片| 免费在线观看成人毛片| 两性午夜刺激爽爽歪歪视频在线观看| 看非洲黑人一级黄片| 天天躁日日操中文字幕| 大又大粗又爽又黄少妇毛片口| 中文字幕av成人在线电影| 我要搜黄色片| 亚洲人成网站高清观看| av在线观看视频网站免费| 嫩草影院新地址| 亚洲av电影不卡..在线观看| 免费av不卡在线播放| 色综合亚洲欧美另类图片| 国产淫片久久久久久久久| 亚洲人成网站在线观看播放| avwww免费| 午夜福利在线在线| 日日干狠狠操夜夜爽| 久久久久精品国产欧美久久久| av.在线天堂| 97在线视频观看| 午夜福利视频1000在线观看| 伊人久久精品亚洲午夜| 国产av不卡久久| 美女高潮的动态| 日韩精品中文字幕看吧| 人人妻,人人澡人人爽秒播| 久久久精品94久久精品| 噜噜噜噜噜久久久久久91| 欧美三级亚洲精品| 十八禁网站免费在线| 麻豆成人午夜福利视频| 久久久久久久亚洲中文字幕| 女同久久另类99精品国产91| 久久婷婷人人爽人人干人人爱| 欧美精品国产亚洲| 精品免费久久久久久久清纯| 国产人妻一区二区三区在| 又爽又黄无遮挡网站| 欧美日韩国产亚洲二区| 亚洲成人精品中文字幕电影| 中文字幕熟女人妻在线| 久久九九热精品免费| 99精品在免费线老司机午夜| 午夜a级毛片| 国产亚洲av嫩草精品影院| 18+在线观看网站| 一级a爱片免费观看的视频| 97热精品久久久久久| 亚洲丝袜综合中文字幕| 在线观看一区二区三区| 51国产日韩欧美| 欧美绝顶高潮抽搐喷水| 国产精品三级大全| 亚洲不卡免费看| 国产午夜福利久久久久久| 亚洲无线在线观看| 色尼玛亚洲综合影院| 一区二区三区高清视频在线| 亚洲精品影视一区二区三区av| 国产高清三级在线| 人人妻,人人澡人人爽秒播| 毛片一级片免费看久久久久| 国产精品不卡视频一区二区| 香蕉av资源在线| 乱系列少妇在线播放| 婷婷亚洲欧美| 激情 狠狠 欧美| 久久精品国产亚洲av天美| 日韩欧美 国产精品| 一级毛片久久久久久久久女| 亚洲三级黄色毛片| 午夜影院日韩av| 性色avwww在线观看| 日韩成人伦理影院| 国产欧美日韩精品亚洲av| 免费看日本二区| 国产精品久久久久久久电影| 一区二区三区四区激情视频 | 小蜜桃在线观看免费完整版高清| 91麻豆精品激情在线观看国产| 三级毛片av免费| 日韩精品中文字幕看吧| 日韩强制内射视频| 丝袜喷水一区| 看非洲黑人一级黄片| 亚洲国产精品成人久久小说 | 91av网一区二区| 日本欧美国产在线视频| 亚洲av五月六月丁香网| 我要搜黄色片| 99久久精品一区二区三区| 五月伊人婷婷丁香| 国产成人aa在线观看| 国产白丝娇喘喷水9色精品| 国产成人福利小说| 久久精品国产鲁丝片午夜精品| 亚洲真实伦在线观看| 日韩高清综合在线| 久久久久久久久中文| 热99在线观看视频| 啦啦啦啦在线视频资源| 夜夜看夜夜爽夜夜摸| 偷拍熟女少妇极品色| 亚洲人成网站高清观看| 亚洲美女视频黄频| 久久99热6这里只有精品| 国产视频一区二区在线看| 听说在线观看完整版免费高清| 国产av一区在线观看免费| 亚洲精品日韩在线中文字幕 | 欧洲精品卡2卡3卡4卡5卡区| 人人妻人人澡欧美一区二区| 国产人妻一区二区三区在| 亚洲av五月六月丁香网| 天美传媒精品一区二区| 人人妻人人看人人澡| 97超视频在线观看视频| 国产av一区在线观看免费| 亚洲自拍偷在线| 国产精品一区www在线观看| 三级男女做爰猛烈吃奶摸视频| 欧美性感艳星| 午夜日韩欧美国产| 欧美性猛交黑人性爽| 一本精品99久久精品77| avwww免费| 麻豆精品久久久久久蜜桃| 成人特级黄色片久久久久久久| 偷拍熟女少妇极品色| 国内精品久久久久精免费| 亚洲精品成人久久久久久| 午夜爱爱视频在线播放| 国产亚洲精品久久久久久毛片| 欧美日韩乱码在线| 亚洲无线观看免费| 久久久精品欧美日韩精品| 国产高清三级在线| 久久精品国产自在天天线| 久99久视频精品免费| 在线a可以看的网站| 亚洲av一区综合| 人妻丰满熟妇av一区二区三区| 国产不卡一卡二| 午夜激情欧美在线| 成人av一区二区三区在线看| 久久午夜亚洲精品久久| 亚洲电影在线观看av| 毛片一级片免费看久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲无线在线观看| 大又大粗又爽又黄少妇毛片口| 久久亚洲精品不卡| 免费看av在线观看网站| 又粗又爽又猛毛片免费看| 18禁裸乳无遮挡免费网站照片| 成人永久免费在线观看视频| 美女xxoo啪啪120秒动态图| av专区在线播放| 亚洲av美国av| 天堂影院成人在线观看| 草草在线视频免费看| 在线看三级毛片| 亚洲第一区二区三区不卡| 色播亚洲综合网| 久久久久久久久中文| 国产成年人精品一区二区| 国产精品永久免费网站| 国产不卡一卡二| 97超级碰碰碰精品色视频在线观看| 亚洲电影在线观看av| 久久韩国三级中文字幕| 久久久精品大字幕| 亚州av有码| 免费观看人在逋| 日韩欧美一区二区三区在线观看| 亚洲熟妇中文字幕五十中出| 亚洲av免费高清在线观看| 人妻丰满熟妇av一区二区三区| 97碰自拍视频| 欧美性感艳星| 一卡2卡三卡四卡精品乱码亚洲| 天堂av国产一区二区熟女人妻| 91av网一区二区| 亚洲四区av| 欧美激情国产日韩精品一区| av免费在线看不卡| 能在线免费观看的黄片| 又黄又爽又免费观看的视频| 99久久精品热视频| 好男人在线观看高清免费视频| 伦精品一区二区三区| 最新中文字幕久久久久| 99久国产av精品国产电影| 国产aⅴ精品一区二区三区波| 成人特级黄色片久久久久久久| 国产精品三级大全| 国产淫片久久久久久久久| 在线国产一区二区在线| 亚洲激情五月婷婷啪啪| 舔av片在线| 国产精品一区www在线观看| 精品久久久久久久久久免费视频| 小说图片视频综合网站| 最新中文字幕久久久久| 亚洲av第一区精品v没综合| 国产麻豆成人av免费视频| 精品一区二区三区av网在线观看| 国产成人freesex在线 | 久久久精品94久久精品| 毛片女人毛片| 最近最新中文字幕大全电影3| 少妇高潮的动态图| 一级a爱片免费观看的视频| 美女xxoo啪啪120秒动态图| 一级黄片播放器| 少妇的逼水好多| 国产乱人视频| 淫秽高清视频在线观看| 狂野欧美白嫩少妇大欣赏| 联通29元200g的流量卡| 3wmmmm亚洲av在线观看| 欧美色视频一区免费| 欧美日韩国产亚洲二区| 丰满人妻一区二区三区视频av| 免费看av在线观看网站| 国产女主播在线喷水免费视频网站 | 亚州av有码| 精品日产1卡2卡| 国产乱人视频| 成人特级av手机在线观看| 久久午夜亚洲精品久久| 国产毛片a区久久久久| 男女边吃奶边做爰视频| 久久久久久久久久久丰满| 男人舔奶头视频| 午夜福利视频1000在线观看| 日韩国内少妇激情av| 欧美一区二区国产精品久久精品| 久久久a久久爽久久v久久| 久久精品国产自在天天线| 美女 人体艺术 gogo| 国产精品亚洲一级av第二区| 观看美女的网站| 日韩大尺度精品在线看网址| 精品人妻视频免费看| 99riav亚洲国产免费| 搡老熟女国产l中国老女人| 日韩欧美一区二区三区在线观看| 三级毛片av免费| 中文亚洲av片在线观看爽| 乱人视频在线观看| 欧美另类亚洲清纯唯美| 精品一区二区三区人妻视频| 久久6这里有精品| 中文资源天堂在线| 亚洲av一区综合| 欧美+日韩+精品| 国产精品爽爽va在线观看网站| 国产精品野战在线观看| 国产一区二区在线av高清观看| 99精品在免费线老司机午夜| 舔av片在线| 久久久久免费精品人妻一区二区| 嫩草影院新地址| 欧美绝顶高潮抽搐喷水| 性插视频无遮挡在线免费观看| 在线观看66精品国产| 亚洲最大成人av| 丝袜喷水一区| 国产精品国产三级国产av玫瑰| 淫妇啪啪啪对白视频| 亚洲内射少妇av| av专区在线播放| 舔av片在线| 91在线观看av| 69人妻影院| 欧美激情久久久久久爽电影| 中国美女看黄片| 熟女电影av网| 好男人在线观看高清免费视频| 国产免费一级a男人的天堂| 别揉我奶头~嗯~啊~动态视频| 国产免费一级a男人的天堂| 国产真实乱freesex| 日韩av不卡免费在线播放| 51国产日韩欧美| 日产精品乱码卡一卡2卡三| 夜夜夜夜夜久久久久| 最好的美女福利视频网| 成人二区视频| 一区福利在线观看| 久久人妻av系列| 亚洲18禁久久av| 欧美最黄视频在线播放免费| 日韩强制内射视频| 成人永久免费在线观看视频| 国产精品女同一区二区软件| 亚洲性久久影院| 亚洲av第一区精品v没综合| 免费看光身美女| 黄色视频,在线免费观看| 欧美最黄视频在线播放免费| 热99在线观看视频| 最近2019中文字幕mv第一页| 久久鲁丝午夜福利片| 国产一区二区在线av高清观看| 欧美另类亚洲清纯唯美| 91久久精品电影网| 天天躁日日操中文字幕| 国产成人一区二区在线| 黑人高潮一二区| 中文字幕免费在线视频6| 久久久午夜欧美精品| 亚洲国产色片| 国产高清不卡午夜福利| 一个人观看的视频www高清免费观看| 国产av麻豆久久久久久久| 日本精品一区二区三区蜜桃| 亚洲国产精品国产精品| 亚洲电影在线观看av| 欧美性感艳星| 亚洲成av人片在线播放无| 久久久久久国产a免费观看| 大又大粗又爽又黄少妇毛片口| 丰满人妻一区二区三区视频av| av黄色大香蕉| 日韩欧美一区二区三区在线观看| 精品午夜福利视频在线观看一区| 此物有八面人人有两片| 国产精品嫩草影院av在线观看| 国产毛片a区久久久久| 久久综合国产亚洲精品| 欧美潮喷喷水| 亚洲一级一片aⅴ在线观看| 人妻久久中文字幕网| av在线观看视频网站免费| 午夜激情福利司机影院| 久久精品人妻少妇| 简卡轻食公司| 1000部很黄的大片| 亚洲欧美清纯卡通| 亚洲人成网站在线播放欧美日韩| 欧美xxxx黑人xx丫x性爽| 亚洲精品一卡2卡三卡4卡5卡| 国产成人a∨麻豆精品| 三级经典国产精品| 中国美女看黄片| 99视频精品全部免费 在线| 欧美xxxx黑人xx丫x性爽| 小说图片视频综合网站| 亚洲va在线va天堂va国产| 夜夜看夜夜爽夜夜摸| 亚洲欧美清纯卡通| 欧美三级亚洲精品| 一本一本综合久久| 亚洲成人久久性| 人妻少妇偷人精品九色| 亚洲av中文av极速乱| 秋霞在线观看毛片| 人人妻人人澡欧美一区二区| 国产在视频线在精品| 亚洲国产欧美人成| 网址你懂的国产日韩在线| 能在线免费观看的黄片| 欧美日韩精品成人综合77777| 午夜精品一区二区三区免费看| 欧美日韩乱码在线| 欧美高清成人免费视频www| 听说在线观看完整版免费高清| 欧美日韩国产亚洲二区| 国产老妇女一区| 欧美性感艳星| 国产午夜精品论理片| 草草在线视频免费看| 少妇的逼水好多| 一区二区三区免费毛片| 欧美最黄视频在线播放免费| 成人特级av手机在线观看| 免费av不卡在线播放| 久久人人爽人人片av| av在线观看视频网站免费| 国产又黄又爽又无遮挡在线| 在线观看美女被高潮喷水网站| 精品久久久久久久久久免费视频| 国产黄色小视频在线观看| 精品一区二区三区人妻视频| 人妻少妇偷人精品九色| 草草在线视频免费看| 午夜日韩欧美国产| 久久综合国产亚洲精品| 亚洲一区二区三区色噜噜| 男人狂女人下面高潮的视频| www日本黄色视频网| 国产欧美日韩精品一区二区| 久久人人爽人人爽人人片va| 久久久久国产精品人妻aⅴ院| 精品一区二区三区视频在线| 九九在线视频观看精品| 亚洲欧美清纯卡通| 国产免费男女视频| 91在线精品国自产拍蜜月| 中文字幕久久专区| 波多野结衣高清作品| 亚洲国产精品国产精品| 国产精品免费一区二区三区在线| 晚上一个人看的免费电影| 少妇被粗大猛烈的视频| 97超级碰碰碰精品色视频在线观看| 亚洲成人精品中文字幕电影| 国产精品美女特级片免费视频播放器| 日韩人妻高清精品专区| 一本久久中文字幕| 成人特级黄色片久久久久久久| 超碰av人人做人人爽久久| 日日摸夜夜添夜夜添av毛片| 99久久精品热视频| 日本一本二区三区精品| 欧美日韩在线观看h| 国产 一区 欧美 日韩| 搡老妇女老女人老熟妇| 亚洲性久久影院| 麻豆国产av国片精品| 亚洲欧美成人综合另类久久久 | 搡女人真爽免费视频火全软件 | 久久精品国产亚洲av涩爱 | 禁无遮挡网站| 免费av不卡在线播放| 午夜福利视频1000在线观看| 看非洲黑人一级黄片| 一级av片app| 成年女人永久免费观看视频| 黄色日韩在线| 搡老妇女老女人老熟妇| 中文字幕av成人在线电影| 乱码一卡2卡4卡精品| 国产成人91sexporn| 久久久精品大字幕| 国产熟女欧美一区二区| 久久中文看片网| 国产淫片久久久久久久久| 人妻制服诱惑在线中文字幕| 97热精品久久久久久| 人妻夜夜爽99麻豆av| 欧美xxxx黑人xx丫x性爽| 国内精品一区二区在线观看| av免费在线看不卡| 国产乱人偷精品视频| 日韩三级伦理在线观看| 秋霞在线观看毛片| 国产午夜精品久久久久久一区二区三区 |