• <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久视频精品免费| 国产蜜桃级精品一区二区三区| 国产一区二区在线av高清观看| 男插女下体视频免费在线播放| 国产精品野战在线观看| 噜噜噜噜噜久久久久久91| 国产一级毛片七仙女欲春2| 国产激情欧美一区二区| 成人鲁丝片一二三区免费| 一本综合久久免费| 国产一区二区三区在线臀色熟女| 在线十欧美十亚洲十日本专区| 欧美色视频一区免费| 精品久久久久久,| 18+在线观看网站| 亚洲一区高清亚洲精品| 午夜老司机福利剧场| 久久香蕉国产精品| eeuss影院久久| 91九色精品人成在线观看| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久人妻精品电影| 国产精品亚洲一级av第二区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产成人aa在线观看| 叶爱在线成人免费视频播放| 黄片小视频在线播放| 成年版毛片免费区| 九色国产91popny在线| 在线播放国产精品三级| 日韩欧美在线二视频| 国产一区二区在线观看日韩 | h日本视频在线播放| 91在线观看av| 亚洲精品乱码久久久v下载方式 | 免费看光身美女| 亚洲精品乱码久久久v下载方式 | 欧美高清成人免费视频www| 亚洲男人的天堂狠狠| 日韩高清综合在线| 观看美女的网站| 变态另类丝袜制服| 老熟妇仑乱视频hdxx| 18+在线观看网站| 亚洲avbb在线观看| 日本一二三区视频观看| 哪里可以看免费的av片| 亚洲av成人不卡在线观看播放网| 国产伦一二天堂av在线观看| 老司机福利观看| 亚洲五月天丁香| 一二三四社区在线视频社区8| 日本免费一区二区三区高清不卡| 久久久久久久久中文| 婷婷丁香在线五月| 一本久久中文字幕| 综合色av麻豆| 色综合亚洲欧美另类图片| 熟女电影av网| 搞女人的毛片| 99热这里只有是精品50| 少妇的丰满在线观看| 嫁个100分男人电影在线观看| 亚洲av美国av| 亚洲美女黄片视频| 亚洲va日本ⅴa欧美va伊人久久| 成人18禁在线播放| 免费在线观看亚洲国产| 嫁个100分男人电影在线观看| 国产黄色小视频在线观看| 宅男免费午夜| 国产不卡一卡二| av天堂在线播放| 亚洲国产欧美网| 国产在视频线在精品| 成年女人毛片免费观看观看9| 51午夜福利影视在线观看| 欧美黄色片欧美黄色片| 亚洲国产精品sss在线观看| 久久久久国产精品人妻aⅴ院| 亚洲人成网站高清观看| 老司机午夜十八禁免费视频| 亚洲精品久久国产高清桃花| 国产高清视频在线观看网站| av中文乱码字幕在线| 网址你懂的国产日韩在线| 日韩大尺度精品在线看网址| 夜夜夜夜夜久久久久| 青草久久国产| 欧美性猛交╳xxx乱大交人| 99久久精品热视频| 日本五十路高清| 国产午夜福利久久久久久| 国产精品1区2区在线观看.| 国产亚洲精品av在线| 美女被艹到高潮喷水动态| 亚洲国产中文字幕在线视频| 桃色一区二区三区在线观看| 欧美色视频一区免费| 老鸭窝网址在线观看| 18禁在线播放成人免费| 亚洲国产欧洲综合997久久,| 国产三级中文精品| 九九在线视频观看精品| 韩国av一区二区三区四区| 成人18禁在线播放| 最新在线观看一区二区三区| 日本一二三区视频观看| 亚洲av一区综合| www.色视频.com| 一进一出抽搐gif免费好疼| 亚洲av五月六月丁香网| 国产色婷婷99| 久久草成人影院| 好男人在线观看高清免费视频| 日韩亚洲欧美综合| 日韩高清综合在线| 一a级毛片在线观看| 成人国产一区最新在线观看| 老熟妇乱子伦视频在线观看| 免费观看精品视频网站| 欧美成人免费av一区二区三区| 亚洲最大成人手机在线| 欧美在线一区亚洲| 国产一区二区激情短视频| 亚洲av美国av| av欧美777| 欧美日本亚洲视频在线播放| 在线观看免费视频日本深夜| 中文资源天堂在线| 免费在线观看日本一区| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区精品| 久久久久国内视频| 亚洲av熟女| 九色成人免费人妻av| 午夜福利欧美成人| 亚洲成人精品中文字幕电影| 国产真人三级小视频在线观看| 五月伊人婷婷丁香| 午夜福利18| 精品一区二区三区视频在线观看免费| 国产精品亚洲一级av第二区| 欧美日韩瑟瑟在线播放| 午夜日韩欧美国产| 中文字幕久久专区| 欧美极品一区二区三区四区| 少妇人妻精品综合一区二区 | 国产亚洲精品久久久com| 免费看美女性在线毛片视频| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久久久久久久| 久久久国产成人精品二区| 99久久精品国产亚洲精品| 国产淫片久久久久久久久 | 最好的美女福利视频网| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人av激情在线播放| 波多野结衣巨乳人妻| 一本精品99久久精品77| 亚洲av不卡在线观看| 桃色一区二区三区在线观看| 99热6这里只有精品| 性色avwww在线观看| 亚洲狠狠婷婷综合久久图片| 有码 亚洲区| 亚洲av成人av| 18禁国产床啪视频网站| 亚洲精品色激情综合| 国产欧美日韩一区二区三| 亚洲精品成人久久久久久| 精品久久久久久,| 丁香欧美五月| 国产不卡一卡二| 国内揄拍国产精品人妻在线| www日本黄色视频网| 在线视频色国产色| 大型黄色视频在线免费观看| 99视频精品全部免费 在线| 国产免费一级a男人的天堂| 国产爱豆传媒在线观看| 国产视频一区二区在线看| 网址你懂的国产日韩在线| 一本综合久久免费| 欧美日本视频| 国产亚洲精品综合一区在线观看| 免费观看精品视频网站| 国产在线精品亚洲第一网站| 天堂影院成人在线观看| 国产v大片淫在线免费观看| 国产日本99.免费观看| 国产欧美日韩一区二区精品| 午夜福利视频1000在线观看| 狠狠狠狠99中文字幕| 国产伦人伦偷精品视频| 听说在线观看完整版免费高清| 成人特级av手机在线观看| 亚洲精品粉嫩美女一区| 亚洲av成人av| 天堂动漫精品| 亚洲av免费在线观看| 免费看美女性在线毛片视频| 日韩高清综合在线| 一边摸一边抽搐一进一小说| 天天躁日日操中文字幕| 一夜夜www| 亚洲七黄色美女视频| 日本免费a在线| 天美传媒精品一区二区| 叶爱在线成人免费视频播放| 99国产综合亚洲精品| 听说在线观看完整版免费高清| 舔av片在线| 亚洲性夜色夜夜综合| 天天添夜夜摸| 国内毛片毛片毛片毛片毛片| 我要搜黄色片| 69人妻影院| 欧美在线一区亚洲| 啪啪无遮挡十八禁网站| 成人午夜高清在线视频| 一个人免费在线观看电影| 亚洲七黄色美女视频| 天天躁日日操中文字幕| 亚洲精品久久国产高清桃花| 日本一二三区视频观看| 午夜福利视频1000在线观看| 给我免费播放毛片高清在线观看| 给我免费播放毛片高清在线观看| 久久久色成人| 欧美日韩一级在线毛片| 亚洲精品乱码久久久v下载方式 | 亚洲av日韩精品久久久久久密| 亚洲内射少妇av| 女生性感内裤真人,穿戴方法视频| 国产国拍精品亚洲av在线观看 | 国产精品女同一区二区软件 | 欧美一级毛片孕妇| 欧美一级毛片孕妇| 女人被狂操c到高潮| 少妇的丰满在线观看| 亚洲人成网站在线播| 精品99又大又爽又粗少妇毛片 | 国产极品精品免费视频能看的| 身体一侧抽搐| av在线蜜桃| 国内精品一区二区在线观看| 91在线观看av| 无限看片的www在线观看| 一个人看视频在线观看www免费 | 深夜精品福利| 少妇的逼好多水| 九九热线精品视视频播放| 精品久久久久久久毛片微露脸| 亚洲av日韩精品久久久久久密| 亚洲狠狠婷婷综合久久图片| 亚洲欧美日韩卡通动漫| 日本一二三区视频观看| 亚洲第一欧美日韩一区二区三区| 亚洲男人的天堂狠狠| 精品国内亚洲2022精品成人| 欧美乱码精品一区二区三区| 手机成人av网站| h日本视频在线播放| 国内少妇人妻偷人精品xxx网站| 午夜影院日韩av| 国内揄拍国产精品人妻在线| 深夜精品福利| 九九热线精品视视频播放| 亚洲av电影在线进入| 日韩欧美精品免费久久 | 露出奶头的视频| 亚洲欧美一区二区三区黑人| 男插女下体视频免费在线播放| xxxwww97欧美| 国产精品永久免费网站| 午夜日韩欧美国产| 黄色日韩在线| 大型黄色视频在线免费观看| 久久婷婷人人爽人人干人人爱| 99在线人妻在线中文字幕| 欧美一区二区亚洲| 天天躁日日操中文字幕| 91久久精品国产一区二区成人 | www日本黄色视频网| 每晚都被弄得嗷嗷叫到高潮| 亚洲av二区三区四区| 免费在线观看日本一区| 精品日产1卡2卡| 日韩有码中文字幕| 午夜老司机福利剧场| 亚洲黑人精品在线| 俄罗斯特黄特色一大片| 少妇裸体淫交视频免费看高清| 亚洲av成人av| 五月玫瑰六月丁香| 婷婷六月久久综合丁香| 每晚都被弄得嗷嗷叫到高潮| 少妇熟女aⅴ在线视频| 成年版毛片免费区| 悠悠久久av| 男人舔女人下体高潮全视频| 十八禁网站免费在线| 欧美性猛交╳xxx乱大交人| 久久久久久久精品吃奶| 好男人在线观看高清免费视频| 成年女人毛片免费观看观看9| 天堂影院成人在线观看| 看黄色毛片网站| 白带黄色成豆腐渣| 一级a爱片免费观看的视频| 久久婷婷人人爽人人干人人爱| 国产亚洲欧美在线一区二区| 一个人观看的视频www高清免费观看| 99精品久久久久人妻精品| 免费看日本二区| 久久人妻av系列| 一级毛片高清免费大全| 岛国在线观看网站| 不卡一级毛片| 18禁黄网站禁片免费观看直播| 国产av不卡久久| 老鸭窝网址在线观看| 久久草成人影院| 夜夜躁狠狠躁天天躁| 亚洲一区二区三区不卡视频| 成年女人看的毛片在线观看| 国内少妇人妻偷人精品xxx网站| 女人高潮潮喷娇喘18禁视频| 在线免费观看的www视频| 叶爱在线成人免费视频播放| 最后的刺客免费高清国语| 女人被狂操c到高潮| 亚洲在线自拍视频| 黄色成人免费大全| 啦啦啦免费观看视频1| 色综合站精品国产| 在线观看午夜福利视频| 在线播放国产精品三级| 午夜两性在线视频| 99热只有精品国产| www.www免费av| 欧美日韩亚洲国产一区二区在线观看| 国内精品久久久久精免费| 亚洲 欧美 日韩 在线 免费| 丁香欧美五月| 淫秽高清视频在线观看| 麻豆国产av国片精品| 日本三级黄在线观看| 国产伦一二天堂av在线观看| 悠悠久久av| 男女做爰动态图高潮gif福利片| 国产精品久久视频播放| 狂野欧美白嫩少妇大欣赏| 国产麻豆成人av免费视频| 99久久精品国产亚洲精品| 内地一区二区视频在线| 国内精品久久久久精免费| 淫妇啪啪啪对白视频| 中文字幕av在线有码专区| 亚洲国产精品sss在线观看| 美女黄网站色视频| 精品欧美国产一区二区三| 在线播放国产精品三级| 日本熟妇午夜| 日韩免费av在线播放| 国产欧美日韩精品亚洲av| 久久精品91无色码中文字幕| 国内精品久久久久久久电影| 淫妇啪啪啪对白视频| 国产免费男女视频| 日本成人三级电影网站| 国产真人三级小视频在线观看| 丰满乱子伦码专区| 国产成人欧美在线观看| 看黄色毛片网站| 国产高清视频在线播放一区| 亚洲18禁久久av| 最近在线观看免费完整版| 啦啦啦韩国在线观看视频| 两人在一起打扑克的视频| 国产精品电影一区二区三区| 波多野结衣高清无吗| h日本视频在线播放| 给我免费播放毛片高清在线观看| 欧美成人性av电影在线观看| 精品乱码久久久久久99久播| 亚洲真实伦在线观看| 一本综合久久免费| 免费av观看视频| 看免费av毛片| 可以在线观看毛片的网站| 国产美女午夜福利| 国产亚洲精品一区二区www| av专区在线播放| 在线国产一区二区在线| 亚洲成人精品中文字幕电影| 日韩有码中文字幕| 两个人的视频大全免费| 夜夜爽天天搞| 亚洲av五月六月丁香网| 亚洲欧美一区二区三区黑人| 日本熟妇午夜| 日韩欧美免费精品| 欧美乱色亚洲激情| 国产精品精品国产色婷婷| 国产精品野战在线观看| 亚洲精品色激情综合| 成人18禁在线播放| 国产一区在线观看成人免费| 国产成人欧美在线观看| 村上凉子中文字幕在线| 老司机午夜十八禁免费视频| 国产精品综合久久久久久久免费| 久久午夜亚洲精品久久| 全区人妻精品视频| 亚洲五月天丁香| 婷婷亚洲欧美| 叶爱在线成人免费视频播放| 国产亚洲精品一区二区www| 深爱激情五月婷婷| 国产精品99久久久久久久久| 我要搜黄色片| 亚洲va日本ⅴa欧美va伊人久久| 极品教师在线免费播放| 一本综合久久免费| 日日干狠狠操夜夜爽| 无限看片的www在线观看| 国产精品精品国产色婷婷| 香蕉久久夜色| 亚洲av第一区精品v没综合| 无限看片的www在线观看| 午夜亚洲福利在线播放| 欧美日韩瑟瑟在线播放| 色播亚洲综合网| 少妇高潮的动态图| 亚洲精华国产精华精| 99久久九九国产精品国产免费| 18禁美女被吸乳视频| 一个人看视频在线观看www免费 | 成人av在线播放网站| 国产中年淑女户外野战色| 免费观看的影片在线观看| 国产成人av教育| 中出人妻视频一区二区| 午夜激情福利司机影院| 精品熟女少妇八av免费久了| 亚洲成a人片在线一区二区| 级片在线观看| 色综合婷婷激情| 好男人电影高清在线观看| 99精品久久久久人妻精品| 久久久久免费精品人妻一区二区| 最新美女视频免费是黄的| 午夜免费男女啪啪视频观看 | 琪琪午夜伦伦电影理论片6080| www.www免费av| 最新美女视频免费是黄的| 国产高清视频在线播放一区| 琪琪午夜伦伦电影理论片6080| 亚洲成a人片在线一区二区| 成人一区二区视频在线观看| 精品国产亚洲在线| 偷拍熟女少妇极品色| 亚洲久久久久久中文字幕| 91在线观看av| 搡老岳熟女国产| 一区二区三区免费毛片| 国产亚洲精品综合一区在线观看| 日韩欧美精品v在线| 国产精品一及| 国产麻豆成人av免费视频| 国产精品一区二区免费欧美| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 欧美日韩乱码在线| 51国产日韩欧美| 人妻丰满熟妇av一区二区三区| 又黄又爽又免费观看的视频| 亚洲中文日韩欧美视频| 国产成人av激情在线播放| 国产精品99久久99久久久不卡| 国产精品免费一区二区三区在线| 日韩有码中文字幕| 免费av观看视频| 色噜噜av男人的天堂激情| 久久久国产成人精品二区| 在线观看免费视频日本深夜| 亚洲乱码一区二区免费版| 变态另类丝袜制服| 麻豆国产av国片精品| 啦啦啦观看免费观看视频高清| 两个人的视频大全免费| 丰满人妻一区二区三区视频av | 国产成人福利小说| 夜夜夜夜夜久久久久| 一个人看的www免费观看视频| 久久香蕉国产精品| 久久草成人影院| 一个人看的www免费观看视频| 亚洲熟妇中文字幕五十中出| 欧美极品一区二区三区四区| 变态另类成人亚洲欧美熟女| 19禁男女啪啪无遮挡网站| 最新美女视频免费是黄的| 人人妻,人人澡人人爽秒播| av视频在线观看入口| 精品久久久久久久人妻蜜臀av| 亚洲自拍偷在线| 十八禁人妻一区二区| 亚洲va日本ⅴa欧美va伊人久久| www.色视频.com| 亚洲国产色片| 黄色成人免费大全| 丰满人妻一区二区三区视频av | 九九热线精品视视频播放| 国产美女午夜福利| 欧美日韩精品网址| 精品久久久久久久毛片微露脸| 午夜老司机福利剧场| 色视频www国产| 国产精品1区2区在线观看.| 亚洲精品美女久久久久99蜜臀| 美女高潮的动态| 我的老师免费观看完整版| 久久久国产精品麻豆| 婷婷丁香在线五月| 老司机在亚洲福利影院| 久久亚洲精品不卡| 国产亚洲精品久久久com| 国产又黄又爽又无遮挡在线| 国产99白浆流出| 午夜福利免费观看在线| 一本久久中文字幕| 国产黄片美女视频| 亚洲18禁久久av| 法律面前人人平等表现在哪些方面| 久久欧美精品欧美久久欧美| 两个人视频免费观看高清| 两个人看的免费小视频| 久久久成人免费电影| 亚洲人与动物交配视频| 日本黄大片高清| 一个人免费在线观看的高清视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲片人在线观看| 欧美另类亚洲清纯唯美| 男人舔女人下体高潮全视频| 少妇的丰满在线观看| 日韩成人在线观看一区二区三区| 亚洲精品456在线播放app | 制服丝袜大香蕉在线| 性欧美人与动物交配| 国产精品98久久久久久宅男小说| 琪琪午夜伦伦电影理论片6080| 特大巨黑吊av在线直播| 亚洲乱码一区二区免费版| 国产欧美日韩一区二区精品| 亚洲一区二区三区不卡视频| 国产高清视频在线观看网站| aaaaa片日本免费| 三级男女做爰猛烈吃奶摸视频| 成人特级av手机在线观看| 成人精品一区二区免费| 黄片小视频在线播放| 亚洲成a人片在线一区二区| 午夜福利18| 亚洲成人久久爱视频| 亚洲av成人不卡在线观看播放网| 国产精华一区二区三区| 色老头精品视频在线观看| 国产伦人伦偷精品视频| 国产高清三级在线| 亚洲国产精品999在线| 中亚洲国语对白在线视频| 欧美高清成人免费视频www| 99热只有精品国产| 日本一本二区三区精品| 国内揄拍国产精品人妻在线| e午夜精品久久久久久久| 神马国产精品三级电影在线观看| 国产成人影院久久av| 欧美激情久久久久久爽电影| 午夜福利免费观看在线| 免费人成在线观看视频色| 欧美午夜高清在线| 欧美最黄视频在线播放免费| 男女床上黄色一级片免费看| 成人一区二区视频在线观看| 亚洲av美国av| 国产成人欧美在线观看| 男女做爰动态图高潮gif福利片| 成熟少妇高潮喷水视频| 国产蜜桃级精品一区二区三区| 亚洲久久久久久中文字幕| 丰满人妻一区二区三区视频av | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 母亲3免费完整高清在线观看| 亚洲欧美日韩东京热| 一二三四社区在线视频社区8| 国产精品亚洲av一区麻豆| 免费无遮挡裸体视频| 国产黄色小视频在线观看| 高潮久久久久久久久久久不卡| 日韩精品青青久久久久久| 国产精品精品国产色婷婷| 日韩欧美在线乱码| 久久久久久久久久黄片| 亚洲一区高清亚洲精品| 亚洲成人精品中文字幕电影| 熟女人妻精品中文字幕| 国产视频一区二区在线看| 午夜福利视频1000在线观看| 国产毛片a区久久久久| 噜噜噜噜噜久久久久久91|