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

    利用背景噪聲研究紫坪鋪水庫水位加卸載對庫區(qū)地下介質(zhì)波速的影響

    2015-09-03 11:04:08安艷茹張曉東
    中國地震 2015年4期
    關(guān)鍵詞:紫坪波速蓄水

    安艷茹 張曉東

    1)中國地震局地震預(yù)測研究所,北京市復(fù)興路63號 100036

    2)中國地震臺(tái)網(wǎng)中心,北京市西城區(qū)三里河南橫街5號 100045

    0 引言

    中國是世界上水庫數(shù)量最多的國家,截至2011年底,全國已建成水庫97246座,總庫容8104.1億m3,在建水庫756座,總庫容 1219.0億 m3(孫振剛等,2013)。中國也是世界上水庫地震多發(fā)地之一,尤其是近些年在地震多發(fā)的川滇地區(qū)建設(shè)了大量高壩、大庫容的梯級電站,水庫地震時(shí)有發(fā)生,大型水庫的水庫地震問題亦引起了地震學(xué)者的高度重視。汶川地震后,由于紫坪鋪水庫與其特殊的時(shí)空關(guān)系,更是引發(fā)廣泛關(guān)注(雷興林等,2008;Shemin et al,2009;馬文濤等,2011;Deng et al,2010;盧顯等,2010)。雷興林等(2008)認(rèn)為,紫坪鋪水庫在蓄水過程中對其地下的龍門山中央斷層和山前斷層有明顯的作用;Shemin等(2009)認(rèn)為蓄水后的庫侖應(yīng)力增量足以引發(fā)地震;馬文濤等(2011)認(rèn)為汶川地震與紫坪鋪水庫在時(shí)空上均有一定的聯(lián)系;Deng等(2010)認(rèn)為紫坪鋪水庫周邊的M-t圖像沒有增強(qiáng),庫水不能滲透到汶川地震深度14km的震源處;盧顯等(2010)認(rèn)為都江堰震群與紫坪鋪水庫無關(guān)。無論如何,水庫的蓄水對地下介質(zhì)產(chǎn)生影響都是個(gè)不爭的事實(shí),而利用數(shù)字地震資料研究庫區(qū)介質(zhì)和應(yīng)力場變化則是十分有效的方法。目前有關(guān)庫區(qū)水位變化與庫區(qū)地震的關(guān)系、庫區(qū)介質(zhì)變化與庫區(qū)地震的關(guān)系等方面的研究較多(周斌等,2010;王惠琳等,2012a、2012b;盧顯等,2013),但有關(guān)庫區(qū)水位變化與庫區(qū)介質(zhì)變化的研究還相對較少。本文基于紫坪鋪水庫數(shù)字地震臺(tái)網(wǎng)的連續(xù)波形數(shù)據(jù),利用噪聲互相關(guān)技術(shù)研究水位的加卸載和滲透過程對庫區(qū)地下介質(zhì)波速的影響,以求深化對庫區(qū)地下結(jié)構(gòu)、狀態(tài)和動(dòng)力學(xué)過程的認(rèn)識。

    1 紫坪鋪水庫及地震臺(tái)網(wǎng)概況

    紫坪鋪水庫位于四川省都江堰市區(qū)西北方向9km處的岷江上游麻溪鄉(xiāng),距成都市區(qū)約60km,距汶川8.0級大地震微觀震中最近處僅6km左右。紫坪鋪水庫于2001年3月29日正式動(dòng)工興建,2006年12月竣工,最大壩高156m,總庫容11.12×108m3。本文利用四川省地震局水庫地震研究所提供的水位數(shù)據(jù),繪制了紫坪鋪水庫水位變化圖(圖1)。由圖1可見,2005年9月30日下閘蓄水,數(shù)天內(nèi)水位迅速升高至死水位817m以上,蓄水后最高水位為875.18m。2005年9月30日~2007年底,紫坪鋪水庫共經(jīng)歷了3次大規(guī)模蓄水、2次泄水過程。2005年12月5日水位達(dá)835.91m,2006年10月14日水位上升到最高值875.18m,2007年12月12日水位再次上升到873.39m。

    圖1 紫坪鋪水庫水位變化

    紫坪鋪水庫位于青藏高原東緣的龍門山造山帶的中段。庫區(qū)及鄰區(qū)被主干斷裂分隔,可劃分為茂汶韌性剪切帶、中央推覆構(gòu)造帶、龍門山前緣拆離帶、前陸擴(kuò)展變形帶和川西前陸盆地等5個(gè)地質(zhì)構(gòu)造單元,它們在物質(zhì)組成、構(gòu)造層次和變形樣式等方面具有明顯的差異。紫坪鋪水庫庫體坐落于龍門山前緣拆離帶內(nèi),鄰近中央推覆構(gòu)造帶和前陸擴(kuò)展變形帶(周斌等,2010)。庫區(qū)地下水分為碎屑巖類裂隙孔隙水和碳酸鹽巖巖溶水等2大類(圖2),如圖2所示,庫區(qū)西南主要巖性為碳酸鹽巖,東北主要巖性為碎屑巖和部分碳酸鹽巖。

    紫坪鋪水庫數(shù)字遙測地震臺(tái)網(wǎng)于2004年8月16日正式采集地震信息,2005年6月27日通過驗(yàn)收(胡先明等,2006)。臺(tái)網(wǎng)包括7個(gè)數(shù)字遙測地震臺(tái),使用短周期地震儀,頻帶寬度1~40Hz。此外,該臺(tái)網(wǎng)共享四川區(qū)域數(shù)字臺(tái)網(wǎng)中YZP(油榨坪)臺(tái)的數(shù)據(jù)。8個(gè)地震臺(tái)大致均勻地分布在庫區(qū)周邊,平均臺(tái)間距約為14km,臺(tái)站分布如圖3所示,臺(tái)站及儀器參數(shù)見表1。

    圖2 庫壩區(qū)地下水類型(據(jù)王云基(2001))

    圖3 紫坪鋪水庫臺(tái)網(wǎng)臺(tái)站分布

    表1 紫坪鋪水庫臺(tái)網(wǎng)臺(tái)站及儀器參數(shù)表

    2 觀測資料

    本研究使用紫坪鋪臺(tái)網(wǎng)及YZP臺(tái)垂直分量的連續(xù)波形數(shù)據(jù),時(shí)間范圍2005年1月1日~2008年1月1日。圖4為各臺(tái)站數(shù)據(jù)連續(xù)記錄情況,可用數(shù)據(jù)的天數(shù)為1003天,最長連續(xù)空缺時(shí)間為27天。

    圖4 各臺(tái)數(shù)據(jù)連續(xù)情況

    3 理論基礎(chǔ)與數(shù)據(jù)處理

    利用背景噪聲測量地下介質(zhì)波速連續(xù)變化的主要技術(shù)路徑是:比較經(jīng)驗(yàn)格林函數(shù)(CCFs)與參考格林函數(shù)(REF)。其中,CCFs表征的是一段時(shí)間內(nèi)地下介質(zhì)的狀態(tài),而 REF代表的是地下介質(zhì)的背景狀態(tài)。計(jì)算過程分為4步:①計(jì)算臺(tái)站對間不同時(shí)間段的噪聲互相關(guān)函數(shù),即經(jīng)驗(yàn)格林函數(shù)CCFs;②獲取參考格林函數(shù)REF,并計(jì)算每個(gè)CCFs與REF的走時(shí)延遲Δt;③對不同臺(tái)站對間不同時(shí)間錯(cuò)動(dòng)的走時(shí)延遲進(jìn)行平均,并使用一個(gè)相對速度均勻一致變化(Δv/v=常數(shù))的簡單模型(Lecocq et al,2014)來表征庫區(qū)地下介質(zhì)波速的變化。

    本文使用MSNoise軟件(Lecocq et al,2014)處理數(shù)據(jù)。該軟件已應(yīng)用于奧克蘭火山地區(qū)(Kasper,2014)及地震前后(Berkeley Seismological Laboratory,2014)的分析研究中,并得到了穩(wěn)定可靠的結(jié)果。

    3.1 噪聲互相關(guān)提取經(jīng)驗(yàn)格林函數(shù)CCFs

    3.1.1 單臺(tái)數(shù)據(jù)預(yù)處理

    本文研究的是庫區(qū)地下介質(zhì)微小的波速變化,因此能否得到穩(wěn)定可靠的CCFs非常關(guān)鍵。首先,我們必須通過數(shù)據(jù)預(yù)處理,獲得較為純凈的背景噪聲,以保證CCFs的質(zhì)量。將每個(gè)臺(tái)站的連續(xù)波形數(shù)據(jù)合并成每天一個(gè)文件,然后進(jìn)行預(yù)處理,包括:去平均、尖滅處理,高通、低通濾波,降采樣,時(shí)間域歸一化,譜白化處理等。通常,波形數(shù)據(jù)總會(huì)存在一個(gè)非零的均值,這會(huì)影響到數(shù)據(jù)的分析;為了使數(shù)據(jù)標(biāo)準(zhǔn)化,適用于各種標(biāo)準(zhǔn)公式,必須在數(shù)據(jù)分析前去平均值(鄭治真,1979)。另外,在對數(shù)據(jù)進(jìn)行譜域操作,如FFT、濾波時(shí),若數(shù)據(jù)的兩端不為零,則會(huì)出現(xiàn)譜域假象,因而實(shí)際數(shù)據(jù)經(jīng)常需要做尖滅處理,以使數(shù)據(jù)兩端在短時(shí)間窗內(nèi)逐漸變成零值。本文研究的是庫區(qū)淺層介質(zhì)波速的變化,關(guān)注的是相對高頻的信息,所以我們選擇了0.1~2.0Hz的濾波參數(shù)。為了減少計(jì)算量和存儲(chǔ)量,將數(shù)據(jù)降采樣到10Hz。為了減少地震信號、儀器異常信號和臺(tái)站附近非穩(wěn)定噪聲源的影響,我們對數(shù)據(jù)進(jìn)行了時(shí)間域歸一化處理(Tukey,1962)。為了使信號在不同頻率上的能量更為均衡,我們對數(shù)據(jù)進(jìn)行了譜白化處理。

    3.1.2 互相關(guān)計(jì)算提取經(jīng)驗(yàn)格林函數(shù)

    利用背景噪聲提取地下結(jié)構(gòu)信息設(shè)想的探索最初可追溯到20世紀(jì)50、60年代,Aki(1957)認(rèn)為可以用背景噪聲提取地下結(jié)構(gòu)中的面波頻散性質(zhì);Claerbout(1968)提出可以利用背景噪聲恢復(fù)一維層狀介質(zhì)的反射響應(yīng)。類似想法的首次成功應(yīng)用是在太陽地震學(xué)研究中,即通過對太陽表面噪聲進(jìn)行互相關(guān)計(jì)算,成功地提取出了聲波時(shí)距曲線(Duvall et al,1993)。隨后,噪聲互相關(guān)技術(shù)在超聲學(xué)領(lǐng)域獲得重大進(jìn)展。Weaver等(2001)在聲學(xué)上給出了隨機(jī)波場的互相關(guān)特性,即通過從發(fā)散的或隨機(jī)的波場中提取格林函數(shù)來計(jì)算地球的彈性響應(yīng)。證明該特性的實(shí)例就是基于一個(gè)彈性體中有發(fā)散場的模型

    式中,x為點(diǎn)的位置;t為時(shí)間;i為虛數(shù)單位;an為模態(tài)激發(fā)函數(shù);un和ωn為地球的本征方程和本征頻率。發(fā)散場的一個(gè)重要性質(zhì)就是模態(tài)振幅是不相關(guān)的隨機(jī)變量

    式中,δnm為克羅內(nèi)克函數(shù);F為頻譜能量密度。因式(2)的交叉項(xiàng)在取平均時(shí)消失,所以場中x、y點(diǎn)的相互關(guān)系簡化為

    式中,τ為時(shí)間錯(cuò)動(dòng)。比較發(fā)現(xiàn),式(3)與x、y點(diǎn)間真實(shí)的格林函數(shù)相比,僅差一個(gè)振幅因子F。因此,根據(jù)足夠長時(shí)間內(nèi)得到的場和場的關(guān)系,可從發(fā)散場中提取2點(diǎn)的格林函數(shù),這正是噪聲成像的理論基礎(chǔ)。之后,該方法在海洋聲學(xué)(Roux et al,2004)和地震學(xué)(Shapiro et al,2004、2005)等領(lǐng)域得到了迅速的發(fā)展。

    數(shù)學(xué)上,頻率域與時(shí)間域的互相關(guān)計(jì)算是等價(jià)的。若2個(gè)序列x(t)、y(t)的傅里葉變換分別是X(f)、Y(f),則頻率域的互相關(guān)公式為

    其中,X*(f)為X(f)的復(fù)共軛;而C(f)的傅里葉反變換正是時(shí)間域的互相關(guān)c(t)。

    本文將1天的數(shù)據(jù)分成48段×30min(50%重疊),在頻率域?qū)γ總€(gè)臺(tái)站對的每段數(shù)據(jù)做互相關(guān)運(yùn)算,并疊加形成每天的經(jīng)驗(yàn)格林函數(shù)CCF。在測量波速時(shí),會(huì)因?yàn)榕cREF相關(guān)度較差的CCF而產(chǎn)生強(qiáng)烈的變化;因此我們分別以10、30、50天作為滑動(dòng)窗,疊加獲得 CCFs,以提高相關(guān)度,獲得更為可靠的波速變化。圖5中,以BAJ-GHS臺(tái)站對為例,展示了互相關(guān)提取的經(jīng)驗(yàn)格林函數(shù)以50天為滑動(dòng)窗疊加后的能量干涉圖,并繪出了不同長度的滑動(dòng)窗獲得的CCFs與REF的互相關(guān)系數(shù)曲線。由圖5可見,10天滑動(dòng)窗的互相關(guān)系數(shù)曲線抖動(dòng)較明顯,30、50天滑動(dòng)窗的互相關(guān)系數(shù)曲線較為平滑,能夠突出顯著的變化。

    圖5 BAJ-GHS臺(tái)站對50天疊加的CCFs的能量干涉圖(a)、紫坪鋪水庫水位變化(b)、CCFs與REF互相關(guān)系數(shù)的變化(c)

    3.2 計(jì)算走時(shí)延遲Δt

    為了研究庫區(qū)水位對地下不同深度介質(zhì)波速的影響,分析水位壓力及滲透的作用過程,我們分別對 1~2s(0.5~1.0Hz)、2~4s(0.25~0.50Hz)、4~8s(0.125~0.250Hz)等 3個(gè)周期進(jìn)行計(jì)算。根據(jù)瑞雷面波的特性(圖6),上述3個(gè)周期的波分別對深度0~2、1~4、1~8km的介質(zhì)變化敏感。

    圖6 瑞雷面波波速在不同周期的敏感性(據(jù)Liu等(2014))

    3.2.1 獲取參考格林函數(shù)REF

    獲取參考格林函數(shù)REF有2種方法,即絕對疊加法和相對疊加法(Lecocq et al,2014)。由于水庫蓄水導(dǎo)致的地下介質(zhì)波速變化并不劇烈,因此我們選擇絕對疊加法,即疊加所有的CCF成為REF。圖7為獲得的28對臺(tái)站的參考格林函數(shù),此處進(jìn)行了反序疊加,旨在展示的方便。

    3.2.2 計(jì)算走時(shí)延遲

    本文使用移動(dòng)窗口互譜方法計(jì)算相對走時(shí)變化,該方法由 Ratdomopurbo等(1995)提出,最早被Poupinet等(1984)用于研究地震對之間的相對速度變化,其優(yōu)點(diǎn)是在頻率域進(jìn)行計(jì)算,可以明確相關(guān)函數(shù)中相關(guān)信號的帶寬。Brenguier等(2008a、2008b)利用該方法先后研究了法國富爾奈斯火山地下介質(zhì)的相對波速變化及2004年P(guān)arkfield地區(qū)6.0級地震前后沿圣安德烈斯斷層地下介質(zhì)的相對波速變化。Clarke等(2011)詳細(xì)闡述了移動(dòng)窗口互譜方法的原理和步驟,并驗(yàn)證了該方法監(jiān)測地殼波速變化的分辨率和準(zhǔn)確率。

    首先,我們將序列CCFs與REF分成許多重疊的窗,去平均、尖滅處理,然后進(jìn)行傅里葉變換得到Fcur(f)和Fref(f)。互譜X(f)被定義為

    式中,f為頻率,星號表示復(fù)共軛性。在頻率域,兩序列的相似程度通過能量密度間的相關(guān)系數(shù)C(f)來評估。由式(5)推導(dǎo)得出

    圖7 28對臺(tái)站的參考格林函數(shù)(反序疊加)

    式中,φ(f)為解纏相位,它與頻率f成線性比例關(guān)系

    對于不同的頻率f和相位φ,可通過最小二乘加權(quán)線性回歸獲得m,對應(yīng)的誤差為em,對應(yīng)的權(quán)重為wj

    式中,Cj是相關(guān)系數(shù);Xj為互譜。這種權(quán)重公式的優(yōu)勢在于,能夠在相關(guān)系數(shù)相對為常數(shù)而互譜能量變化的情況下給出較為不同的權(quán)重。根據(jù)式(7)中Δt與m的關(guān)系,將m與em分別除以2π,即可獲得臺(tái)站對間對應(yīng)于不同時(shí)間錯(cuò)動(dòng)的走時(shí)延遲Δt及誤差eΔt。

    3.3 計(jì)算平均相對波速變化

    假設(shè)相對波速變化Δv/v在空間上是均勻的,那么它與相對走時(shí)變化Δt/t是相反數(shù),即Δv/v=-Δt/t(Ratdomopurboet al,1995)??紤]到臺(tái)站間距、信號的能量強(qiáng)弱、CCFs與 REF的相關(guān)程度等因素,本文選擇了互相關(guān)正負(fù)錯(cuò)動(dòng)的8~20s部分(圖5),此部分為面波及尾波。我們將28組臺(tái)站對間不同時(shí)間錯(cuò)動(dòng)的走時(shí)延遲Δt及對應(yīng)的誤差eΔt進(jìn)行平均,得到庫區(qū)地下介質(zhì)平均的走時(shí)延遲及平均誤差對于3個(gè)周期,我們分別通過最小二乘加權(quán)線性回歸獲得為了減小誤差,擬合時(shí)要求直線強(qiáng)制通過原點(diǎn)(Lecocqet al,2014),則

    其中,pi為權(quán)重

    b對應(yīng)的方差為

    最終,根據(jù)式(9),我們獲得了庫區(qū)地下介質(zhì)平均的相對波速變化(圖 8)。

    4 結(jié)論與討論

    本文的研究時(shí)間段包含了紫坪鋪水庫3次大規(guī)模蓄水和2次泄水的過程,通過噪聲互相關(guān)提取了庫區(qū)28組臺(tái)站對間的經(jīng)驗(yàn)格林函數(shù)CCFs。圖5中以BAJ-GHS臺(tái)站對為例展示了臺(tái)站對間的互相關(guān)結(jié)果,這對臺(tái)站跨越了西南至東北的整個(gè)庫區(qū)(圖3),十分具有代表性。從圖5中可以看出,水位加卸載時(shí)CCFs與REF的互相關(guān)系數(shù)明顯降低,并在時(shí)間上有一定延遲,這定性地表征了水位影響了庫區(qū)地下介質(zhì)波速的變化。本文利用互相關(guān)結(jié)果,通過移動(dòng)窗口互譜方法,分別對 1~2s(0.5~1.0Hz)、2~4s(0.25~0.50Hz)、4~8s(0.125~0.250Hz)等3個(gè)周期,計(jì)算了28組臺(tái)站對間不同時(shí)間錯(cuò)動(dòng)的走時(shí)延遲,3個(gè)周期分別對應(yīng)著0~2、1~4、1~8km的深度。最終,使用面波及尾波部分,通過加權(quán)線性回歸捕捉到了庫區(qū)地下介質(zhì)平均的相對波速變化(圖8)。可以看出,庫區(qū)地下介質(zhì)相對波速變化不超過0.1%。水庫水位的加卸載主要通過壓力和滲透2種作用影響地下介質(zhì)。Biot(1956)認(rèn)為,孔隙流體會(huì)增加壓縮波的速度并降低剪切波的速度。庫區(qū)以碳酸鹽巖為主,具備較好的滲透條件。圖8的3個(gè)周期中,1~2s對應(yīng)的2km內(nèi)的介質(zhì)相對波速變化最大,變化幅度超過了0.05%,這是由于在同樣的水位壓力下,淺層的滲透作用更強(qiáng)烈。

    圖8 庫區(qū)地下介質(zhì)平均相對波速變化

    4.1 第1次蓄水及泄水

    由圖8可見,在2005年9月30日第1次大規(guī)模蓄水前,水位與相對波速變化沒有明顯關(guān)系。蓄水后數(shù)10天內(nèi),水位由760.36m迅速升高至835.91m,3個(gè)周期對應(yīng)的3個(gè)深度的介質(zhì)相對波速隨之同時(shí)迅速降低,說明此時(shí)庫水壓力作用是導(dǎo)致相對波速變化的主要原因,而滲透作用只影響了2km以內(nèi)的淺層介質(zhì)。之后,水庫開始泄水,對應(yīng)的相對波速也逐漸升高。

    4.2 第2次蓄水及泄水

    第2次蓄水開始于2006年8月15日,2006年10月14日水位上升到最高值875.18m,前述3個(gè)深度的介質(zhì)相對波速隨之降低,存在明顯的時(shí)間延遲。其中,1~2、2~4s對應(yīng)的相對波速幾乎同時(shí)降到最小值,而4~8s的相對波速在近1個(gè)月后才降到最小值。這說明在承受同樣的庫水壓力下,滲透作用已成為影響相對波速變化的主要因素,且已影響至深度4km左右的斷層。在第2次泄水過程中,3個(gè)深度對應(yīng)的相對波速變化略有延遲,但幾乎同時(shí)升高至最大值,這說明此時(shí)滲透作用已影響至深度8km左右的斷層。

    4.3 第3次蓄水

    第3次蓄水于2007年12月12日上升到873.39m,3個(gè)周期的相對波速幾乎同時(shí)降到最小值,與水位變化非常一致,互相關(guān)系數(shù)分別達(dá)-0.81,-0.66,-0.76(圖9)。此時(shí)庫水壓力和滲透作用共同影響著介質(zhì)相對波速的變化,水已滲透至深度為8km以下的斷層。

    圖9 第3次蓄水時(shí)的水位與相對波速變化的最小二乘擬合曲線

    綜上,本文分別對紫坪鋪水庫3次蓄水及2次泄水過程進(jìn)行了分析,討論了不同階段影響相對波速變化的主要因素,追蹤了水在不同深度介質(zhì)中的滲透過程。結(jié)果表明:在紫坪鋪水庫大規(guī)模蓄水前,相對波速變化與水位變化沒有明顯的相關(guān)性;而在之后的水位加卸載過程中,表現(xiàn)出了明顯的負(fù)相關(guān)性,即水位升高,相對波速降低,水位降低,相對波速升高。分析認(rèn)為,第1次蓄水時(shí),壓力是導(dǎo)致介質(zhì)相對波速變化的主要因素;而后2次蓄水時(shí),滲透作用成為主要因素。在第2次蓄水至最高點(diǎn)時(shí),滲透作用影響至深度4km左右的斷層,在第2次泄水至最低點(diǎn)時(shí),滲透作用已影響至深度8km左右的斷層。

    盧顯等(2010)通過研究水庫水位與水庫區(qū)域地震頻度的關(guān)系發(fā)現(xiàn),在2005年9月、2006年9月和2007年9月這3個(gè)水位快速上升的月份,地震頻次也相應(yīng)增加,且地震的震源深度絕大部分集中在5km處。根據(jù)本文的研究,這3個(gè)月份對應(yīng)著相對波速由0值左右迅速降低的過程。因此對于地下淺層介質(zhì)相對波速連續(xù)變化的研究,有可能作為水庫地震預(yù)測的依據(jù),具有十分重要的意義。

    致謝:本文使用的大量連續(xù)波形及水庫資料均由四川省地震局提供。研究過程中,本文得到了馬延路博士、楊志高博士、周龍泉研究員、杜方研究員、戴仕貴高工、謝蓉華高工和韓進(jìn)高工的支持與幫助,在此一并表示衷心的感謝。同時(shí),非常感謝審稿專家提出的寶貴意見。

    猜你喜歡
    紫坪波速蓄水
    2022年巴西水電站水庫蓄水率達(dá)50%
    基于實(shí)測波速探討地震反射波法超前預(yù)報(bào)解譯標(biāo)志
    西藏在建大型水電站工程通過蓄水驗(yàn)收
    金馬河安全泄流管理體系的構(gòu)建及建議
    四川水利(2017年1期)2017-03-31 09:25:13
    強(qiáng)震下紫坪鋪壩前大型古滑坡體變形破壞效應(yīng)
    吉林地區(qū)波速比分布特征及構(gòu)造意義
    紫坪鋪電廠3號水輪發(fā)電機(jī)組下導(dǎo)振擺超標(biāo)分析與處理
    基于分位數(shù)回歸的剪切波速變化規(guī)律
    北京20座下凹式立交橋汛前可蓄水
    雨水集流微蓄水技術(shù)的實(shí)施
    欧美绝顶高潮抽搐喷水| 最新中文字幕久久久久 | 精品国产乱码久久久久久男人| 中国美女看黄片| 99热6这里只有精品| 成人欧美大片| 每晚都被弄得嗷嗷叫到高潮| 丰满人妻一区二区三区视频av | 色精品久久人妻99蜜桃| 精品久久蜜臀av无| 又黄又爽又免费观看的视频| 最新中文字幕久久久久 | 在线国产一区二区在线| 免费看a级黄色片| 成年女人永久免费观看视频| 国产成人av教育| 免费观看精品视频网站| 成人av一区二区三区在线看| 床上黄色一级片| 午夜福利视频1000在线观看| 成熟少妇高潮喷水视频| 可以在线观看毛片的网站| 国产成人精品久久二区二区免费| 免费人成视频x8x8入口观看| 国产午夜精品论理片| 亚洲国产日韩欧美精品在线观看 | 国产免费男女视频| 大型黄色视频在线免费观看| 男女做爰动态图高潮gif福利片| a级毛片a级免费在线| 日韩av在线大香蕉| 18美女黄网站色大片免费观看| 免费在线观看成人毛片| 中文字幕最新亚洲高清| 成人永久免费在线观看视频| 搡老妇女老女人老熟妇| or卡值多少钱| 欧美成人性av电影在线观看| 午夜福利视频1000在线观看| 欧美中文日本在线观看视频| 精品国产乱子伦一区二区三区| 又黄又爽又免费观看的视频| 久久精品人妻少妇| 国产精品一及| 51午夜福利影视在线观看| 国产淫片久久久久久久久 | 免费无遮挡裸体视频| 国产高清视频在线播放一区| 老汉色∧v一级毛片| 在线观看免费午夜福利视频| x7x7x7水蜜桃| 国产亚洲欧美98| 亚洲九九香蕉| a在线观看视频网站| av片东京热男人的天堂| 99久国产av精品| 国产精品一区二区免费欧美| 中文字幕av在线有码专区| 欧美不卡视频在线免费观看| 国内少妇人妻偷人精品xxx网站 | 成人av在线播放网站| svipshipincom国产片| 国产精品一区二区精品视频观看| 亚洲激情在线av| 成人无遮挡网站| 特大巨黑吊av在线直播| 国产精品1区2区在线观看.| 国产主播在线观看一区二区| 午夜福利视频1000在线观看| 成人精品一区二区免费| 国产三级中文精品| АⅤ资源中文在线天堂| 天天躁日日操中文字幕| 桃色一区二区三区在线观看| 日本一本二区三区精品| 一个人看视频在线观看www免费 | 日本黄色视频三级网站网址| 欧美av亚洲av综合av国产av| 色吧在线观看| 国产av不卡久久| 露出奶头的视频| 制服丝袜大香蕉在线| 啪啪无遮挡十八禁网站| 欧美日韩福利视频一区二区| 亚洲国产中文字幕在线视频| 亚洲人成网站高清观看| 免费在线观看视频国产中文字幕亚洲| 成人特级黄色片久久久久久久| 国产午夜福利久久久久久| 欧美日韩黄片免| 黄色 视频免费看| 无限看片的www在线观看| 熟妇人妻久久中文字幕3abv| 亚洲中文字幕一区二区三区有码在线看 | 天堂网av新在线| 欧美高清成人免费视频www| 国产黄片美女视频| 少妇丰满av| 欧美成人性av电影在线观看| 人妻久久中文字幕网| 一个人看视频在线观看www免费 | 欧美成人一区二区免费高清观看 | 长腿黑丝高跟| bbb黄色大片| 国产高清三级在线| 精品电影一区二区在线| 日韩欧美免费精品| 国产精品 欧美亚洲| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久人人做人人爽| 一边摸一边抽搐一进一小说| 国产av麻豆久久久久久久| 999久久久国产精品视频| 亚洲无线在线观看| 中出人妻视频一区二区| 麻豆成人av在线观看| 国产蜜桃级精品一区二区三区| 十八禁人妻一区二区| 黑人欧美特级aaaaaa片| av天堂在线播放| 我要搜黄色片| 级片在线观看| 日韩免费av在线播放| 九九热线精品视视频播放| 亚洲人成伊人成综合网2020| 国产aⅴ精品一区二区三区波| 99热这里只有是精品50| 久久香蕉精品热| 亚洲18禁久久av| 一级作爱视频免费观看| 国产伦在线观看视频一区| 国产野战对白在线观看| 日日干狠狠操夜夜爽| 亚洲国产中文字幕在线视频| 99久久久亚洲精品蜜臀av| 琪琪午夜伦伦电影理论片6080| 97超级碰碰碰精品色视频在线观看| 欧美乱色亚洲激情| 午夜免费观看网址| 中文字幕精品亚洲无线码一区| 搡老熟女国产l中国老女人| 无遮挡黄片免费观看| 91在线精品国自产拍蜜月 | 国产黄片美女视频| 变态另类丝袜制服| 波多野结衣高清无吗| 观看美女的网站| 法律面前人人平等表现在哪些方面| 亚洲最大成人中文| 人人妻,人人澡人人爽秒播| 精品国产美女av久久久久小说| 欧美性猛交黑人性爽| 国产精品 国内视频| 美女cb高潮喷水在线观看 | av天堂在线播放| 宅男免费午夜| a级毛片在线看网站| 在线免费观看不下载黄p国产 | 全区人妻精品视频| 国产私拍福利视频在线观看| 男女之事视频高清在线观看| 欧美中文日本在线观看视频| 超碰成人久久| 老司机在亚洲福利影院| 女同久久另类99精品国产91| 欧美日韩综合久久久久久 | 九九在线视频观看精品| 18禁黄网站禁片免费观看直播| 亚洲av电影在线进入| 后天国语完整版免费观看| 国产av一区在线观看免费| 成年版毛片免费区| 男女下面进入的视频免费午夜| 国产伦人伦偷精品视频| 在线免费观看不下载黄p国产 | 午夜福利在线在线| 九色成人免费人妻av| 后天国语完整版免费观看| 国产aⅴ精品一区二区三区波| 一二三四在线观看免费中文在| 啦啦啦免费观看视频1| 亚洲欧美一区二区三区黑人| 国内精品久久久久久久电影| 日韩国内少妇激情av| 色在线成人网| 97超视频在线观看视频| 日韩国内少妇激情av| 久久久久亚洲av毛片大全| 久久中文字幕一级| 法律面前人人平等表现在哪些方面| 国产亚洲av高清不卡| 免费av毛片视频| 香蕉国产在线看| 国产激情欧美一区二区| av天堂在线播放| 亚洲黑人精品在线| 免费在线观看日本一区| 国产99白浆流出| 成人特级av手机在线观看| 久久这里只有精品19| 国产伦精品一区二区三区四那| av在线蜜桃| 亚洲无线在线观看| 欧美丝袜亚洲另类 | 日本黄大片高清| 亚洲中文av在线| 国产69精品久久久久777片 | 欧美zozozo另类| 精品99又大又爽又粗少妇毛片 | 成人av一区二区三区在线看| 成年女人永久免费观看视频| 免费观看的影片在线观看| 中文字幕最新亚洲高清| 在线观看免费视频日本深夜| 99久久精品热视频| 久久久久国内视频| 淫妇啪啪啪对白视频| 久久久久亚洲av毛片大全| 51午夜福利影视在线观看| 日本在线视频免费播放| 欧美精品啪啪一区二区三区| 色综合亚洲欧美另类图片| 国产淫片久久久久久久久 | 村上凉子中文字幕在线| 亚洲熟妇中文字幕五十中出| 最近最新免费中文字幕在线| 色尼玛亚洲综合影院| 免费看光身美女| avwww免费| 亚洲黑人精品在线| 久久久久精品国产欧美久久久| 精品国产亚洲在线| 久久久久久国产a免费观看| 国产97色在线日韩免费| 无限看片的www在线观看| 亚洲av成人不卡在线观看播放网| 午夜亚洲福利在线播放| 俺也久久电影网| 美女扒开内裤让男人捅视频| 欧美日韩国产亚洲二区| 神马国产精品三级电影在线观看| 欧美日韩乱码在线| 日韩欧美在线二视频| 国产午夜福利久久久久久| 国产精品日韩av在线免费观看| 久久这里只有精品19| 脱女人内裤的视频| 人人妻,人人澡人人爽秒播| 91av网站免费观看| 久久国产精品人妻蜜桃| 午夜视频精品福利| 成人特级av手机在线观看| 欧美成人免费av一区二区三区| 国产精品国产高清国产av| 精品国产美女av久久久久小说| 日本一二三区视频观看| 怎么达到女性高潮| 免费看a级黄色片| 天天躁日日操中文字幕| 国产久久久一区二区三区| 一区二区三区国产精品乱码| 日韩欧美 国产精品| 日韩欧美免费精品| 色尼玛亚洲综合影院| 精品一区二区三区视频在线观看免费| 天堂√8在线中文| 久久午夜亚洲精品久久| 美女被艹到高潮喷水动态| 亚洲av成人不卡在线观看播放网| 偷拍熟女少妇极品色| 国产成人啪精品午夜网站| 国产精品亚洲一级av第二区| 99精品欧美一区二区三区四区| 熟女人妻精品中文字幕| 日本熟妇午夜| 国产欧美日韩精品亚洲av| 三级毛片av免费| 国产亚洲精品一区二区www| 久久精品aⅴ一区二区三区四区| 一个人观看的视频www高清免费观看 | 99久久精品国产亚洲精品| 国产成+人综合+亚洲专区| 国产乱人视频| 舔av片在线| 精品99又大又爽又粗少妇毛片 | 国产av不卡久久| 国产视频内射| 国产毛片a区久久久久| 亚洲av免费在线观看| 两性夫妻黄色片| netflix在线观看网站| 成人三级黄色视频| 天天躁日日操中文字幕| 午夜久久久久精精品| 亚洲中文字幕一区二区三区有码在线看 | 俺也久久电影网| 脱女人内裤的视频| 亚洲真实伦在线观看| 欧美性猛交╳xxx乱大交人| 成人鲁丝片一二三区免费| av天堂中文字幕网| 国产91精品成人一区二区三区| 欧美日韩一级在线毛片| 国产精品99久久久久久久久| 草草在线视频免费看| 久9热在线精品视频| 亚洲自偷自拍图片 自拍| 精品国产美女av久久久久小说| 三级毛片av免费| 久久国产精品影院| 欧美性猛交黑人性爽| 在线观看免费午夜福利视频| 美女大奶头视频| 亚洲欧美日韩高清专用| 国产亚洲精品av在线| av欧美777| 美女 人体艺术 gogo| 国产精品久久久久久久电影 | 亚洲熟妇熟女久久| 国产视频内射| 亚洲av五月六月丁香网| 九九久久精品国产亚洲av麻豆 | 真实男女啪啪啪动态图| 日韩免费av在线播放| 一夜夜www| 一本精品99久久精品77| 91字幕亚洲| 99国产精品一区二区蜜桃av| 亚洲欧洲精品一区二区精品久久久| 我的老师免费观看完整版| 国产黄色小视频在线观看| 在线播放国产精品三级| ponron亚洲| 最近最新中文字幕大全免费视频| 夜夜躁狠狠躁天天躁| 亚洲美女黄片视频| 国产亚洲av嫩草精品影院| 国产精品99久久久久久久久| 宅男免费午夜| 97超级碰碰碰精品色视频在线观看| 久久国产精品影院| 听说在线观看完整版免费高清| 91麻豆av在线| 成人三级黄色视频| 亚洲av电影不卡..在线观看| 国产主播在线观看一区二区| 高清在线国产一区| 国产精品久久久久久精品电影| 99热这里只有精品一区 | 99久久精品国产亚洲精品| 此物有八面人人有两片| e午夜精品久久久久久久| 欧美又色又爽又黄视频| 两人在一起打扑克的视频| av中文乱码字幕在线| 国产精品美女特级片免费视频播放器 | 免费观看人在逋| 国产熟女xx| 中文字幕熟女人妻在线| 国产精品99久久久久久久久| 成人午夜高清在线视频| 99国产精品一区二区三区| 黄频高清免费视频| 少妇熟女aⅴ在线视频| 狠狠狠狠99中文字幕| 不卡av一区二区三区| 日本免费a在线| 欧美乱色亚洲激情| 国产成人精品久久二区二区91| 少妇熟女aⅴ在线视频| 狠狠狠狠99中文字幕| 久久久久国产精品人妻aⅴ院| 狠狠狠狠99中文字幕| 成人高潮视频无遮挡免费网站| 精品免费久久久久久久清纯| 日韩精品中文字幕看吧| 亚洲精华国产精华精| 亚洲欧美精品综合一区二区三区| 性色avwww在线观看| 中文亚洲av片在线观看爽| 国模一区二区三区四区视频 | tocl精华| 老熟妇仑乱视频hdxx| 亚洲国产色片| 日韩精品青青久久久久久| 精品无人区乱码1区二区| 91字幕亚洲| 免费在线观看亚洲国产| 亚洲成人精品中文字幕电影| 性欧美人与动物交配| 亚洲国产精品999在线| 婷婷丁香在线五月| 丰满的人妻完整版| 国产成人av激情在线播放| 99在线人妻在线中文字幕| 99久久精品一区二区三区| 国产精品日韩av在线免费观看| 人人妻人人看人人澡| 嫩草影院精品99| 脱女人内裤的视频| 亚洲激情在线av| 日本一二三区视频观看| 亚洲欧美一区二区三区黑人| xxx96com| 日韩成人在线观看一区二区三区| 香蕉国产在线看| av天堂中文字幕网| 日本在线视频免费播放| 大型黄色视频在线免费观看| 香蕉国产在线看| 国语自产精品视频在线第100页| 俄罗斯特黄特色一大片| svipshipincom国产片| 好男人电影高清在线观看| 欧美性猛交╳xxx乱大交人| 观看免费一级毛片| av在线蜜桃| 99国产精品一区二区蜜桃av| 国产激情欧美一区二区| 一进一出抽搐gif免费好疼| 黄频高清免费视频| 亚洲 欧美一区二区三区| 日本与韩国留学比较| 亚洲美女黄片视频| 在线观看日韩欧美| 一区二区三区激情视频| 亚洲黑人精品在线| 国产精品99久久久久久久久| 久久久国产精品麻豆| 夜夜躁狠狠躁天天躁| 国产免费av片在线观看野外av| 免费看日本二区| 亚洲真实伦在线观看| 成人性生交大片免费视频hd| 在线观看66精品国产| 岛国在线免费视频观看| 香蕉久久夜色| 九九热线精品视视频播放| 啦啦啦观看免费观看视频高清| 午夜免费观看网址| 日本在线视频免费播放| 欧美在线黄色| 国产日本99.免费观看| 亚洲 国产 在线| 欧美色视频一区免费| 亚洲中文字幕日韩| 久久天躁狠狠躁夜夜2o2o| www日本在线高清视频| 少妇裸体淫交视频免费看高清| 人人妻人人看人人澡| 国产一区二区激情短视频| 日本撒尿小便嘘嘘汇集6| 午夜亚洲福利在线播放| 国产精品电影一区二区三区| 国产私拍福利视频在线观看| 法律面前人人平等表现在哪些方面| www.自偷自拍.com| 亚洲欧美日韩卡通动漫| 午夜福利在线观看吧| netflix在线观看网站| 变态另类丝袜制服| 久久热在线av| 婷婷六月久久综合丁香| 日本三级黄在线观看| 久久婷婷人人爽人人干人人爱| 欧美日韩亚洲国产一区二区在线观看| 国产不卡一卡二| 欧美另类亚洲清纯唯美| 午夜福利欧美成人| 亚洲欧美精品综合久久99| 亚洲欧美日韩无卡精品| 夜夜夜夜夜久久久久| 久久国产乱子伦精品免费另类| 18禁裸乳无遮挡免费网站照片| 日本成人三级电影网站| 最近最新免费中文字幕在线| 一二三四在线观看免费中文在| 亚洲五月天丁香| 一个人看的www免费观看视频| 欧美成人性av电影在线观看| 免费看光身美女| 国产成人精品久久二区二区91| 老熟妇乱子伦视频在线观看| 国产精品 国内视频| 国产精品1区2区在线观看.| 久久精品91无色码中文字幕| 国产精品99久久99久久久不卡| 亚洲五月天丁香| 在线观看免费视频日本深夜| 午夜免费激情av| 国产99白浆流出| 别揉我奶头~嗯~啊~动态视频| 12—13女人毛片做爰片一| 欧美日韩福利视频一区二区| 成人国产一区最新在线观看| 精品一区二区三区四区五区乱码| 日韩人妻高清精品专区| 99热6这里只有精品| 久久久久亚洲av毛片大全| 在线免费观看不下载黄p国产 | 神马国产精品三级电影在线观看| 后天国语完整版免费观看| 亚洲成人久久性| 可以在线观看的亚洲视频| 嫩草影视91久久| 日本黄色片子视频| 18禁黄网站禁片免费观看直播| 久久久久国产一级毛片高清牌| 丁香欧美五月| 校园春色视频在线观看| 日本 av在线| 亚洲av电影不卡..在线观看| 国产黄片美女视频| 男女视频在线观看网站免费| 国产成+人综合+亚洲专区| 亚洲aⅴ乱码一区二区在线播放| 一二三四在线观看免费中文在| 99riav亚洲国产免费| av在线蜜桃| 狂野欧美白嫩少妇大欣赏| 日本免费一区二区三区高清不卡| 青草久久国产| 免费在线观看影片大全网站| 午夜激情福利司机影院| 无遮挡黄片免费观看| 美女扒开内裤让男人捅视频| 国产精品久久久久久精品电影| 日韩欧美在线二视频| 午夜福利欧美成人| 午夜免费激情av| 一个人看的www免费观看视频| 丰满人妻一区二区三区视频av | 亚洲成av人片免费观看| 色老头精品视频在线观看| 亚洲七黄色美女视频| 天天一区二区日本电影三级| 亚洲国产欧美一区二区综合| 亚洲色图av天堂| 午夜福利在线观看吧| 动漫黄色视频在线观看| 亚洲成av人片免费观看| 1000部很黄的大片| а√天堂www在线а√下载| 亚洲午夜精品一区,二区,三区| 中亚洲国语对白在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 欧美在线一区亚洲| 亚洲国产精品合色在线| 免费人成视频x8x8入口观看| av天堂中文字幕网| 国产精华一区二区三区| 国产精品综合久久久久久久免费| 中文字幕精品亚洲无线码一区| 色在线成人网| 精品久久久久久成人av| 一进一出好大好爽视频| 淫妇啪啪啪对白视频| 亚洲av日韩精品久久久久久密| 日韩大尺度精品在线看网址| 1024香蕉在线观看| 99久久精品一区二区三区| 国产精品精品国产色婷婷| 久久精品国产综合久久久| 午夜亚洲福利在线播放| 精品福利观看| 嫁个100分男人电影在线观看| 亚洲一区二区三区色噜噜| 亚洲成人久久爱视频| 日韩成人在线观看一区二区三区| 精品久久久久久成人av| 精品乱码久久久久久99久播| 欧美最黄视频在线播放免费| 露出奶头的视频| 又黄又爽又免费观看的视频| 欧美xxxx黑人xx丫x性爽| 麻豆久久精品国产亚洲av| 午夜激情福利司机影院| 黄频高清免费视频| 一级毛片精品| 久久久精品欧美日韩精品| 欧美精品啪啪一区二区三区| 在线观看午夜福利视频| 久久国产精品影院| 精品一区二区三区视频在线 | 亚洲国产精品999在线| 亚洲国产精品久久男人天堂| 深夜精品福利| 亚洲av片天天在线观看| 岛国在线免费视频观看| 熟女电影av网| 国产成人一区二区三区免费视频网站| 精品国产美女av久久久久小说| 又黄又粗又硬又大视频| 亚洲人成电影免费在线| 久久精品国产清高在天天线| 亚洲在线自拍视频| 国产精品爽爽va在线观看网站| 亚洲电影在线观看av| 亚洲av电影不卡..在线观看| 精品久久久久久久毛片微露脸| 日韩中文字幕欧美一区二区| 成熟少妇高潮喷水视频| 亚洲精品粉嫩美女一区| www日本在线高清视频| 美女 人体艺术 gogo| 国产成年人精品一区二区| 精品福利观看| 首页视频小说图片口味搜索| 九色成人免费人妻av| 九九在线视频观看精品| 很黄的视频免费| 婷婷六月久久综合丁香| 国产精品野战在线观看| av视频在线观看入口| 色视频www国产|