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

    SV波斜入射下成層盆地地震動時-頻域放大特征研究

    2022-05-11 08:33:02魏成前于彥彥丁海平
    地震工程與工程振動 2022年2期
    關(guān)鍵詞:入射角觀測點震動

    魏成前,于彥彥,丁海平

    (蘇州科技大學(xué)江蘇省結(jié)構(gòu)工程重點實驗室,江蘇蘇州 215011)

    引言

    國內(nèi)外大量的地震觀測記錄、數(shù)值模擬及震害調(diào)查研究結(jié)果表明,沉積盆地作為一種特殊的場地,會顯著放大其內(nèi)部長周期地震動,并延長地震動的持續(xù)時間,這對于盆地內(nèi)的高層建筑、大跨度橋梁和其它自振周期較長的工程結(jié)構(gòu)非常不利。1985年墨西哥8.1級地震中,遠(yuǎn)離震中400 km的墨西哥城中高層建筑破壞嚴(yán)重,震后研究認(rèn)為這與該城市坐落于深厚沉積盆地有關(guān)[1-2]。其后在1995年日本神戶地震中的大阪盆地[3]、1999年中國臺灣集集地震中臺北盆地[4]、2008年四川汶川地震中的渭河盆地[5]及2014年魯?shù)榈卣瘕堫^山鎮(zhèn)盆地[6]等地出現(xiàn)的震害,均表明盆地對震害會有加重作用。

    現(xiàn)有研究表明,入射地震波的角度、幅值和頻譜特性等參數(shù)對盆地地震效應(yīng)具有顯著的影響。張建經(jīng)等[7]研究了SH波斜入射時,地震波的入射角、卓越頻率和盆地傾角對盆地地表地震響應(yīng)的影響。結(jié)果表明,隨著地震波入射角的增大,加速度反應(yīng)譜的峰值有減小的趨勢;地震波的卓越頻率越大,影響盆地地表加速度峰值的地震波入射角的變化范圍就越小。Khanbabazadehand等[8]通過輸入不同幅值的強震地震波研究了盆地介質(zhì)分層和土層參數(shù)的變化對盆地動力響應(yīng)的影響,發(fā)現(xiàn)雙層介質(zhì)盆地的共振周期普遍大于單層盆地,雙層軟黏土盆地的譜放大倍數(shù)最大。Zhu等[9]采用顯式有限差分法詳細(xì)研究了SH波垂直下入射波主頻和盆地結(jié)構(gòu)幾何形狀之間的相互作用及盆地內(nèi)外介質(zhì)阻抗比對盆地地震響應(yīng)的影響,結(jié)果表明盆地地面運動的振幅和分布特點是由入射波與盆地動力特性以及盆地幾何結(jié)構(gòu)相互作用決定的。陳學(xué)良等[10]研究了Rayleigh面波作用下盆地介質(zhì)參數(shù)、盆地尺寸對盆地場地地震動特性的影響規(guī)律,發(fā)現(xiàn)當(dāng)盆地內(nèi)介質(zhì)由軟變硬時,盆地對地表地震動的放大作用由大變小。劉中憲等[11]基于高精度譜元法,研究了不同斷層傾角及不同震源頻率下盆地的地震響應(yīng)規(guī)律,發(fā)現(xiàn)隨著斷層傾角的增大,地表加速度的峰值有減小的趨勢。震源的主頻較高時,盆地的地震動放大效應(yīng)顯著增強。梁佳利等[12]基于粘彈性邊界和等效線性化方法研究了二維均勻盆地的非線性地震反應(yīng)特征,結(jié)果顯示入射波的頻譜和幅值對盆地地震響應(yīng)影響非常明顯。但以往對于成層盆地模型在SV波斜入射下地震動的時、頻域綜合放大特性的研究相對較少。此外,當(dāng)前對于盆地場點的工程地震效應(yīng)評價中,采用較多的方法是根據(jù)場點處的鉆孔資料建立一維土層模型,進而分析地震波垂直入射下地表點的地震動響應(yīng),而忽略了入射角度的影響以及二維盆地模型橫向上的波速不均勻分布所引起的盆地放大作用。

    文中以SV波為入射地震波,結(jié)合顯式有限元法和透射邊界技術(shù),詳細(xì)研究了梯形形狀的二維成層盆地模型,相比盆地各場點處厚度不等的一維水平成層模型在相同輸入波、不同入射角度下地表地震動的放大特征,并以地震動峰值大小、地震動峰值放大系數(shù)和傅里葉譜譜比作為定量分析的指標(biāo)。

    1 計算方法及計算模型

    1.1 計算方法

    模擬地震波斜入射條件下盆地場地地震響應(yīng)時,首先需要計算自由波場,以此作為斜入射條件下盆地模型的波動輸入。均勻或成層半空間自由場可以根據(jù)波動傳播規(guī)律在頻域內(nèi)計算得到[13]。當(dāng)使用有限元方法分析地震波散射問題時,有必要從半無限空間介質(zhì)中切出有限的計算區(qū)域,同時,必須在邊界處引入適當(dāng)?shù)娜斯み吔鐥l件模擬區(qū)域內(nèi)外行波的傳播過程,以達到消除截斷邊界處人工反射波的目的。文中采用透射邊界[14],在模型的左、右兩側(cè)邊界和底邊界設(shè)置二階透射邊界。計算區(qū)域內(nèi)的節(jié)點可以分為內(nèi)節(jié)點和邊界節(jié)點,位于人工邊界上的節(jié)點u0稱為人工邊界節(jié)點,其余的節(jié)點均為內(nèi)節(jié)點u,全部內(nèi)節(jié)點的運動方程為

    式中:M為質(zhì)量矩陣;C為阻尼矩陣;K為剛度矩;P為外力矢量,文中P為0。

    人工邊界節(jié)點運動采用多次透射公式計算,即

    式中:N為透射階數(shù)(文中取2階);0表示人工邊界點;j為與節(jié)點0相鄰的計算點;u0p+1為邊界節(jié)點p+1時刻的位移;u p+1-jj為p+1時刻j計算點的位移。

    關(guān)于散射波場的輸入問題,其輸入是來自半無限域的彈性波。在人工邊界附近區(qū)域,既有穿過人工邊界向外傳播的外行波,又有向計算區(qū)內(nèi)部傳播的內(nèi)行波。由于透射邊界只能模擬外行波,所以必須從總波場中分離出外行波。假設(shè)分離出的外行波位移為us,u為總波場位移,ur為參考波場位移(一般取自由場位移),三者存在以下關(guān)系,

    將式(4)代入式(2)得到邊界節(jié)點p+1時刻的總場位移為

    利用式(5)即可實現(xiàn)散射問題的人工邊界節(jié)點運動方程的求解,此方程也同時實現(xiàn)了地震波的輸入。

    1.2 計算模型

    文中所采用的盆地計算模型如圖1(a)所示,模型尺寸為5 000 m×1 000 m,盆地頂部寬4 000 m,底部寬3 600 m,深度H為200 m,盆地傾角α=45°。盆地內(nèi)土層分3層,厚度從上到下依次為30,70,100 m。模型介質(zhì)參數(shù)見表1。在地表x=100~4 900 m范圍內(nèi)選取55個觀測點,其中x=500~700 m,x=4 300~4 500 m范圍內(nèi)觀測點間隔為40 m,其余間隔均為100 m。盆地內(nèi)最小網(wǎng)格尺寸為2 m,基巖最大網(wǎng)格尺寸為10 m。模擬時間步長0.001 s。為了與盆地模型對比來研究盆地的放大效應(yīng),建立與圖1(a)各觀測點處的土層分布和厚度完全相同的一維水平成層模型作為參考模型,共計55個一維土層模型。以測點7和12為例,其對應(yīng)的一維土層模型分別如圖1(b)、(c)所示。各土層參數(shù)與盆地模型相同,進而模擬相同地震波輸入下一維土層的反應(yīng)。采用寬度0.4 s的SV脈沖波作為基巖輸入地震波,從模型左側(cè)傾斜入射,設(shè)地震波入射角(入射方向與垂直方向的夾角)為θ,其位移時程曲線及其傅里葉譜見圖2。

    圖2 輸入波時程及其傅里葉譜Fig.2 Time history of input wave and its spectrum

    表1 計算模型介質(zhì)參數(shù)Table 1 Physical parameters of the calculation model

    對于該盆地模型,采用剪切波速加權(quán)平均法計算盆地最深處的等效一維土層模型(土層分布如圖1(c))的基階自振頻率f0,計算公式為:

    圖1 計算模型示意圖Fig.1 Schematic of the calculation model

    式中:vsi和hi分別為各土層的剪切波速和厚度;H為土層總厚度(200 m)。根據(jù)式(6)計算得f0=0.7 Hz。2.2節(jié)中將基于該自振頻率將模擬結(jié)果對應(yīng)的頻率無量綱化。

    2 計算結(jié)果

    為研究地震波入射角對成層盆地地表地震動放大特征的影響,利用圖1所示盆地模型和盆地不同場點處相同土層厚度的一維水平成層模型,在保持其余參數(shù)不變的基礎(chǔ)上,依次改變地震波的入射角θ為0°、10°、20°和30°,然后分別模擬得到相應(yīng)入射角下模型的地震動響應(yīng)。需要指出的是,由于文中研究的是平面SV波入射下的盆地場地效應(yīng),當(dāng)SV波的入射角大于一定角度時,反射波就變成了非均勻平面波[15],不在文中討論范圍之內(nèi)。因此,文中入射SV波將具有臨界入射角,即SV波的入射角需小于按以下公式確定的角度:

    式中,β為P波的反射角,當(dāng)β為90°時,入射角θ達到臨界值。由表1可知基巖半空間介質(zhì)Vs=1 650 m/s,Vp=2 858 m/s,代入式(7)計算可得SV波臨界入射角θ約為35o,因此文中最大入射角取為30°。

    2.1 入射角對盆地地震動時域放大特征的影響

    不同入射角下盆地模型模擬得到的地表觀測點水平和垂直分量的位移時程如圖3所示??梢钥闯?,盆地內(nèi)觀測點的位移時程相比基巖更加復(fù)雜、振蕩,地震動持時顯著延長。隨著入射角θ的增大,地表各觀測點直達體波的到時會有滯后的現(xiàn)象。水平分量直達體波的強度逐漸減弱,垂直分量上則逐漸增強。不同入射角度下盆地兩側(cè)邊緣均有明顯的次生面波產(chǎn)生,但隨入射角度的增大,盆地左側(cè)次生面波的強度逐漸大于右側(cè)。不同入射角度下盆地左側(cè)的面波強度有增強的趨勢,右側(cè)面波強度有減弱趨勢,垂直分量的次生面波強度受入射角的影響相對較小。

    圖3 不同入射角下盆地地表觀測點水平(左)和垂直分量(右)的位移時程Fig.3 Displacement time histories of horizontal(left)and vertical(right)components of surface observation points corresponding to different incidence angles

    圖4為不同入射角下盆地地表觀測點的位移峰值(PGD)分布。對于水平分量,當(dāng)入射角θ=0°即垂直入射時,地表位移峰值基本呈對稱分布。隨著入射角的增大,盆地左側(cè)邊緣觀測點的位移峰值普遍大于右側(cè)。盆地邊緣處觀測點(6~10號,45~50號觀測點)的PGD受入射角的影響十分顯著:6~10號觀測點范圍內(nèi),PGD隨入射角的增加而增大;45~50號觀測點范圍內(nèi),PGD隨入射角的增大而減小。盆地中央?yún)^(qū)域PGD隨入射角的增大而減小,但固定入射角下其值始終比較穩(wěn)定,可能和盆地較大的寬度有關(guān)。同時,入射角對PGD在盆地內(nèi)的分布具有重要影響,其中主要影響盆地邊緣效應(yīng)出現(xiàn)的位置。以盆地左側(cè)最強烈地震動區(qū)域為例,當(dāng)入射角θ=0°或10°時,PGD有先增大后至平穩(wěn)的趨勢,PGD最大值發(fā)生在x=900 m處;而入射角θ=20°或30°時,PGD有先增大后減小至平穩(wěn)的趨勢,PGD最大值發(fā)生在x=620 m處??傮w而言,PGD最大值出現(xiàn)的位置有隨入射角的增大而向盆地左側(cè)移動的趨勢。

    圖4 不同入射角度下盆地地表觀測點的位移峰值分布Fig.4 PGD distributions of observation points on the surface of the basin under different incidence angles

    對于垂直分量,與盆地角點相同距離的左側(cè)觀測點的位移峰值明顯大于右側(cè),但最大PGD值兩側(cè)基本相同。在盆地左側(cè)和右側(cè)邊緣處(6~15號,40~50號觀測點范圍)均表現(xiàn)出明顯的盆地邊緣效應(yīng),PGD在此處變化劇烈。PGD最大值與入射角之間具有明顯的規(guī)律性,即隨入射角的增大而增大。盆地中央?yún)^(qū)域的PGD也隨入射角度的增加而增大(10°~30°時增大最顯著),這與水平分量的趨勢相反。不同入射角下盆地左側(cè)和右側(cè)觀測點垂直分量的PGD最大值均出現(xiàn)在6號和50號觀測點處。

    55個觀測點對應(yīng)的一維水平成層模型在不同角度SV波入射下的位移峰值分布如圖5所示。可以看出,同一模型(固定觀測點)相同厚度覆蓋土層下,隨著入射角的增大水平分量PGD呈減小趨勢,垂直分量PGD則呈增大趨勢。同一入射角度下,水平分量上覆蓋土層厚度越大,其PGD值也越大;垂直分量上在測點7處(覆蓋層厚度100 m)地震動幅值最大。此外,垂直入射下水平成層模型垂直分量PGD為0,即垂直入射下,模型沒有產(chǎn)生垂直方向的振動。由此可見,垂直入射下二維盆地模型產(chǎn)生的垂直方向上的位移是由盆地特殊的幾何形狀引起的。

    圖5 不同入射角度下55個一維模型地表觀測點的位移峰值分布Fig.5 PGD distributions of surface observation points of 55 1D models under different incidence angles

    定義盆地地震動放大系數(shù)(AF)為盆地模型模擬得到的地表位移峰值與對應(yīng)場點處的一維水平成層模型模擬得到的地表位移峰值的比值,其分布如圖6所示,圖7為圖6(a)在盆地邊緣附近區(qū)域的局部放大??梢钥闯?,對于水平分量,當(dāng)?shù)卣鸩ù怪比肷鋾r,盆地邊緣區(qū)域觀測點的地震動放大系數(shù)大多小于1。隨著入射角度的增大,放大系數(shù)也逐漸增大,當(dāng)入射角θ=30°時,8號觀測點放大系數(shù)達到了最大值1.43。盆地中央?yún)^(qū)域水平分量地震動放大系數(shù)穩(wěn)定在1左右。由于地震波從左側(cè)入射,盆地右側(cè)區(qū)域?qū)θ肷浣堑淖兓鄬Σ幻舾校卣饎臃糯笙禂?shù)均小于1,且入射角越大,放大系數(shù)總體越小。此外,水平分量上盆地外側(cè)4~5號測點區(qū)域的放大系數(shù)也大于1,應(yīng)主要受盆地構(gòu)造的影響。對于垂直分量,放大系數(shù)隨入射角的增加急劇減小,其放大效應(yīng)在盆地邊緣區(qū)域(6號及50號測點處)最顯著,最大放大系數(shù)可達4.5。相比而言,盆地垂直分量的放大系數(shù)顯著大于水平分量。因此平面波入射下除考慮盆地構(gòu)造的影響外,還需考慮入射角度的影響。

    圖6 不同入射角度下盆地地表觀測點的放大系數(shù)分布Fig.6 Distributions of amplification factor of surface observation points of the basin under different incidence angles

    2.2 入射角對盆地地震動頻域放大特征的影響

    選取盆地內(nèi)6、8、10、12、14和20號觀測點(距盆地角點540~1 700 m范圍,位置如圖1所示),分別定義為點A、B、C、D、E和F。圖8為A~F觀測點在不同入射角下水平分量位移時程的傅里葉譜(本小節(jié)只研究水平分量地震動),其中頻率f基于盆地最深處一維等效土層模型的自振頻率f0進行了無量綱化??梢钥闯觯瑢ν挥^測點,不同入射角度下其頻譜特征基本相似,但卓越頻率(傅里葉譜最大幅值對應(yīng)的頻率)隨入射角度的增大而略有增加。從盆地邊緣到盆地內(nèi)部,觀測點的卓越頻率有從高頻向低頻移動的趨勢,如入射角θ=30°時,A~F號觀測點的卓越頻率分別等于2.87f0、1.91f0、1.93f0、1.31f0、1.13f0、1.39f0。此外,A~E點傅氏譜幅值隨入射角增大而增大,而位于盆地中央的F號觀測點的傅氏譜幅值對入射角的變化相對不敏感,且其幅值小于盆地邊緣處觀測點。

    圖8 不同入射角度下A~F號觀測點的傅里葉譜Fig.8 Fourier spectra of observation points A~F under different incidence angles

    圖8(續(xù))Fig.8(Continued)

    圖9給出了A~F號觀測點在不同入射角下位移時程傅里葉譜比曲線。本小節(jié)定義頻域放大倍數(shù)(AFF)為盆地模型觀測點的位移時程傅里葉譜與一維土層模型位移時程傅里葉譜之比,用來定量描述頻域內(nèi)盆地場地相比一維土層模型的附加放大效應(yīng)。可以看出,地震波從左側(cè)斜入射時,入射角越大,盆地地表觀測點的譜比值越大。如觀測點A在4種入射角下的最大放大倍數(shù)分別為2.12、2.38、2.65和2.97,這與時域結(jié)果對應(yīng)。盆地邊緣處的放大倍數(shù)要大于盆地其它區(qū)域,在C號觀測點(x=700 m)處,頻域放大倍數(shù)達到了最大值3.5,而盆地內(nèi)部F號觀測點的頻域放大系數(shù)總體最小。值得注意的是,A點在f0~2.5f0頻段內(nèi)的AFF值基本都小于1,主要放大頻段集中分布在3f0附近;B~E觀測點的譜比曲線整體呈雙峰分布,一個大致位于f0~2.5f0的范圍(如C點和E點在此頻段內(nèi)的AFF最大值分別達到了2.26和2.12),另一個(也是最主要放大頻段)集中在4f0左右;而F點在較寬頻帶范圍內(nèi)的AFF值均較為均衡。表明相比等效的一維土層模型,盆地邊緣土層較薄區(qū)域?qū)ο鄬Α案哳l”的地震動放大較顯著,且主要放大頻段集中,而靠近盆地內(nèi)部土層較厚,在較寬的頻帶范圍內(nèi)均表現(xiàn)出明顯放大效應(yīng)。此外,固定觀測點處譜比曲線的分布特征相近,基本不受入射角度的影響。

    圖9 不同入射角度下A~F觀測點的譜比曲線Fig.9 Fourier spectra ratio curves of observation points A~F under different incidence angles

    不同入射角下盆地內(nèi)6~50號觀測點譜比最大值分布如圖10所示,可以看出譜比最大值A(chǔ)FFmax與入射角之間有很好的相關(guān)性。總體來說,盆地邊緣區(qū)域的AFFmax大于盆地其它區(qū)域,從6號點到35號觀測點,AFFmax有隨入射角增大而增大的趨勢,但35號觀測點之后(逐漸靠近盆地右側(cè)邊界),AFFmax卻隨入射角度的增大而減小。

    圖10 不同入射角度下盆地內(nèi)觀測點的譜比最大值分布Fig.10 Distributions of the maximum spectral ratio of observation points in the basin under different incidence angles

    3 成層盆地與等效均勻盆地地震動及其放大特征比較

    考慮盆地內(nèi)介質(zhì)分層對不同入射角度的影響,建立和圖1(a)相同尺寸的均勻盆地模型,即盆地內(nèi)土層為均一介質(zhì),其厚度仍為200 m。剪切波速取成層盆地模型的等效剪切波速[16],即

    進而模擬相同入射波(圖2)下均勻盆地的地震響應(yīng)與放大效應(yīng)。

    圖11給出了均勻盆地模型地表觀測點的水平分量位移峰值及相比一維等效土層模型的放大系數(shù)分布??梢钥闯?,相比成層盆地,均勻盆地模型盆地邊緣效應(yīng)遠(yuǎn)小于成層盆地,其PGD最大值為0.374 m,而成層盆地模型PGD最大值為0.619 m,且地震動最大值的位置基本不隨入射角的變化而改變,恒定在12號觀測點處,這與成層盆地模型入射角≤10°時的PGD最大值位置相同。此外,盆地內(nèi)PGD最大值隨入射角的增大而一致減小。從峰值放大系數(shù)來看,在給定入射波下,入射角度較小時(≤10°)成層盆地內(nèi)部的放大系數(shù)(圖6(a))小于等效均勻盆地,入射角大于10°后則大于等效均勻盆地,如均勻盆地模型在30°入射角時盆地內(nèi)放大系數(shù)最大值為1.14(顯著小于成層盆地模型的1.45),出現(xiàn)在11號測點附近。同時,均勻盆地模型的放大系數(shù)分布曲線相比成層模型結(jié)果簡單,盆地內(nèi)最大放大系數(shù)相近,基本不隨入射角的變化而改變。

    圖11 不同入射角度下均勻盆地模型地表觀測點的位移峰值及放大系數(shù)分布Fig.11 Distributions of PGD and amplification factor of observation points on the surface of the single layer basin

    4 結(jié)論

    文中采用有限元法結(jié)合透射邊界模擬SV波斜入射下成層盆地的地震動響應(yīng)。通過改變地震波的入射角度,研究了入射角對成層盆地相比等效一維土層模型的時、頻域放大特征的影響,得到以下結(jié)論:

    (1)入射角θ對盆地地表地震動峰值影響顯著。水平分量上,盆地邊緣區(qū)域的位移峰值隨θ增加而增大,中間及右側(cè)區(qū)域的位移峰值隨之減?。淮怪狈至可?,盆地所有區(qū)域的位移峰值均隨θ增加而增大。

    (2)入射角度對盆地邊緣效應(yīng)區(qū)域的位置有影響。對水平分量,θ大于20°時出現(xiàn)在盆地斜邊區(qū)域范圍內(nèi),小于20°時則位于盆地斜邊范圍右側(cè)的區(qū)域;對垂直分量,邊緣效應(yīng)始終出現(xiàn)在緊鄰盆地角點的位置,不受入射角的影響。

    (3)不同入射角度下,盆地左側(cè)邊緣區(qū)域水平分量地震動放大系數(shù)顯著大于右側(cè),其值在此區(qū)域達到最大,且隨θ增大而增大(最大放大系數(shù)1.43),右側(cè)則反之;垂直分量上左右兩側(cè)放大系數(shù)基本相同,其值隨入射角的增大而減小,最大放大系數(shù)4.5左右。兩分量上放大系數(shù)的分布特征基本不受入射角度的影響,最大放大系數(shù)的位置與最大位移峰值位置基本相同。

    (4)固定觀測點在不同入射角度下的頻譜特征類似,但峰值頻率有隨θ的增大略變大的趨勢。入射角變化時不同位置觀測點的譜比曲線分布特征相似,但譜比值均隨入射角度的增大而增大。不同θ角下成層盆地邊緣區(qū)域以相對“高頻”的放大為主,放大頻段較為集中,向盆地內(nèi)部則在較寬頻帶范圍內(nèi)均有放大。

    (5)盆地內(nèi)介質(zhì)波速的縱向變化對其地震動分布影響顯著。盆地內(nèi)分層介質(zhì)與盆地斜邊的相互作用可能是導(dǎo)致大角度入射時成層盆地地震動放大效應(yīng)相比等效均勻盆地顯著的主要原因,因此在斜入射下盆地地震動的估計中,必須同時考慮盆地內(nèi)介質(zhì)橫向和縱向波速的變化。

    猜你喜歡
    入射角觀測點震動
    一般三棱鏡偏向角與入射角的關(guān)系
    高速公路網(wǎng)連續(xù)式交通量調(diào)查觀測點布設(shè)方法研究
    智能城市(2021年3期)2021-04-12 04:40:50
    震動減脂儀可以減肥?
    洛陽市老城區(qū)西大街空間形態(tài)與熱環(huán)境耦合關(guān)系實測研究
    綠色科技(2019年12期)2019-07-15 11:13:02
    預(yù)制圓柱形鎢破片斜穿甲鋼靶的破孔能力分析*
    水電工程場地地震動確定方法
    振動攪拌 震動創(chuàng)新
    中國公路(2017年18期)2018-01-23 03:00:38
    用經(jīng)典定理證明各向異性巖石界面異常入射角的存在
    張掖市甘州區(qū)代表性觀測點地下水位變化特征分析
    基于升降溫全曲線的鋼筋混凝土梁溫度場分析
    中文字幕制服av| 爱豆传媒免费全集在线观看| 日本精品一区二区三区蜜桃| 女性被躁到高潮视频| 欧美日韩一级在线毛片| 一本一本久久a久久精品综合妖精| 成年人午夜在线观看视频| 欧美一级毛片孕妇| 91av网站免费观看| 男人舔女人的私密视频| 精品第一国产精品| 免费看十八禁软件| 久久精品熟女亚洲av麻豆精品| netflix在线观看网站| 各种免费的搞黄视频| 91成年电影在线观看| 三上悠亚av全集在线观看| 欧美少妇被猛烈插入视频| 亚洲av成人不卡在线观看播放网 | 纵有疾风起免费观看全集完整版| 亚洲伊人色综图| 亚洲中文日韩欧美视频| 桃红色精品国产亚洲av| 日本91视频免费播放| 两个人看的免费小视频| 大香蕉久久成人网| 国产又色又爽无遮挡免| 无限看片的www在线观看| 欧美人与性动交α欧美精品济南到| 亚洲国产精品999| 亚洲精品国产av成人精品| 国产亚洲欧美精品永久| 交换朋友夫妻互换小说| kizo精华| 国产精品国产三级国产专区5o| 亚洲国产中文字幕在线视频| 69av精品久久久久久 | 久久久精品94久久精品| 美女脱内裤让男人舔精品视频| 精品人妻一区二区三区麻豆| 热re99久久国产66热| 一区在线观看完整版| 国产精品二区激情视频| 国产精品99久久99久久久不卡| av一本久久久久| 日韩人妻精品一区2区三区| 国产淫语在线视频| 丝瓜视频免费看黄片| 99国产精品一区二区三区| 狠狠婷婷综合久久久久久88av| 丝袜在线中文字幕| 黄色毛片三级朝国网站| 精品国产乱码久久久久久小说| 国产男女超爽视频在线观看| 一级a爱视频在线免费观看| √禁漫天堂资源中文www| 久久99一区二区三区| 一区二区日韩欧美中文字幕| 女性被躁到高潮视频| 777米奇影视久久| 久久人妻福利社区极品人妻图片| 好男人电影高清在线观看| 国产日韩一区二区三区精品不卡| 亚洲激情五月婷婷啪啪| 一个人免费看片子| a级毛片黄视频| 国产精品 欧美亚洲| 男女下面插进去视频免费观看| 极品少妇高潮喷水抽搐| 大片免费播放器 马上看| 日韩制服骚丝袜av| 国产有黄有色有爽视频| 亚洲专区中文字幕在线| 窝窝影院91人妻| 又大又爽又粗| 亚洲av美国av| 亚洲成人免费电影在线观看| 成人18禁高潮啪啪吃奶动态图| 俄罗斯特黄特色一大片| 国内毛片毛片毛片毛片毛片| 少妇精品久久久久久久| 久久久精品区二区三区| 99国产精品一区二区蜜桃av | 黄色视频不卡| 日本五十路高清| 大片免费播放器 马上看| 国产一区二区在线观看av| 国产亚洲欧美在线一区二区| 纯流量卡能插随身wifi吗| 日日摸夜夜添夜夜添小说| 日本av免费视频播放| 一区二区av电影网| 90打野战视频偷拍视频| 丁香六月天网| 亚洲精品成人av观看孕妇| 亚洲情色 制服丝袜| 视频区图区小说| 一级,二级,三级黄色视频| 一二三四社区在线视频社区8| 999久久久国产精品视频| 黑人操中国人逼视频| 亚洲国产欧美在线一区| 69精品国产乱码久久久| 俄罗斯特黄特色一大片| 久久精品国产亚洲av香蕉五月 | 大片电影免费在线观看免费| 中文字幕人妻丝袜一区二区| 国产亚洲精品一区二区www | 亚洲精品日韩在线中文字幕| 黑人巨大精品欧美一区二区mp4| 国产精品偷伦视频观看了| 一个人免费看片子| 精品高清国产在线一区| 亚洲av日韩在线播放| 久久久久久人人人人人| 手机成人av网站| 大片免费播放器 马上看| 国产亚洲午夜精品一区二区久久| 精品熟女少妇八av免费久了| 中国国产av一级| 亚洲欧美精品自产自拍| 欧美日韩中文字幕国产精品一区二区三区 | 99香蕉大伊视频| 纵有疾风起免费观看全集完整版| 99精品欧美一区二区三区四区| 亚洲精品成人av观看孕妇| 久久久久久久精品精品| 国产成人欧美在线观看 | 国产精品 欧美亚洲| 在线观看免费午夜福利视频| 性色av乱码一区二区三区2| 老汉色av国产亚洲站长工具| 国产一卡二卡三卡精品| av在线app专区| 九色亚洲精品在线播放| 免费在线观看日本一区| 日本撒尿小便嘘嘘汇集6| 国产黄色免费在线视频| 亚洲精品美女久久av网站| 久久精品国产a三级三级三级| 亚洲专区国产一区二区| 中文字幕av电影在线播放| 一级片'在线观看视频| 国产欧美日韩综合在线一区二区| 亚洲激情五月婷婷啪啪| 大型av网站在线播放| 精品久久久久久久毛片微露脸 | 在线天堂中文资源库| 91成人精品电影| 精品福利观看| 一区二区三区四区激情视频| av又黄又爽大尺度在线免费看| 欧美成狂野欧美在线观看| 母亲3免费完整高清在线观看| 免费女性裸体啪啪无遮挡网站| 国产伦理片在线播放av一区| 亚洲第一av免费看| 日韩一区二区三区影片| 又黄又粗又硬又大视频| 午夜久久久在线观看| 99久久精品国产亚洲精品| 人人妻人人澡人人看| 大香蕉久久成人网| 久久免费观看电影| 精品国产乱码久久久久久男人| www.av在线官网国产| av在线播放精品| 黑人欧美特级aaaaaa片| 亚洲avbb在线观看| 国产一区二区在线观看av| 国产色视频综合| 国产精品久久久久久人妻精品电影 | 黄片播放在线免费| 国产精品 欧美亚洲| 狂野欧美激情性xxxx| 久久久久久亚洲精品国产蜜桃av| 法律面前人人平等表现在哪些方面 | 日韩人妻精品一区2区三区| 久久精品亚洲熟妇少妇任你| www.熟女人妻精品国产| 国产麻豆69| 操出白浆在线播放| 亚洲激情五月婷婷啪啪| 丝袜人妻中文字幕| 99re6热这里在线精品视频| 久久国产精品男人的天堂亚洲| 巨乳人妻的诱惑在线观看| 午夜精品国产一区二区电影| 国产精品香港三级国产av潘金莲| 午夜影院在线不卡| 欧美av亚洲av综合av国产av| 涩涩av久久男人的天堂| 国产亚洲av高清不卡| 国产精品99久久99久久久不卡| 亚洲va日本ⅴa欧美va伊人久久 | 97人妻天天添夜夜摸| 国产日韩欧美视频二区| 超碰成人久久| 99九九在线精品视频| 国产1区2区3区精品| 日日爽夜夜爽网站| 精品少妇内射三级| 国产免费一区二区三区四区乱码| 午夜福利视频精品| 久久久久久久久免费视频了| 天堂8中文在线网| 一级毛片电影观看| xxxhd国产人妻xxx| 久久 成人 亚洲| 国产免费一区二区三区四区乱码| 亚洲成人手机| 午夜福利视频精品| 丝袜美足系列| 久久精品国产综合久久久| 亚洲色图综合在线观看| 日韩中文字幕视频在线看片| 亚洲国产av影院在线观看| 久久久久久亚洲精品国产蜜桃av| 狂野欧美激情性xxxx| 亚洲国产欧美在线一区| 色综合欧美亚洲国产小说| 日本vs欧美在线观看视频| 国产区一区二久久| 国产成人啪精品午夜网站| 自拍欧美九色日韩亚洲蝌蚪91| 可以免费在线观看a视频的电影网站| 国产成人免费观看mmmm| 成年女人毛片免费观看观看9 | 精品欧美一区二区三区在线| 夜夜夜夜夜久久久久| √禁漫天堂资源中文www| 国产高清国产精品国产三级| 久久国产精品人妻蜜桃| 两个人看的免费小视频| 啦啦啦在线免费观看视频4| 99久久99久久久精品蜜桃| 搡老熟女国产l中国老女人| 精品久久久精品久久久| 热99国产精品久久久久久7| 精品福利观看| 亚洲精品中文字幕在线视频| 18禁观看日本| 国产淫语在线视频| 日韩欧美一区视频在线观看| 日本精品一区二区三区蜜桃| 日韩熟女老妇一区二区性免费视频| 国产免费现黄频在线看| 亚洲人成77777在线视频| 久久99一区二区三区| 欧美亚洲 丝袜 人妻 在线| 在线观看免费午夜福利视频| 性少妇av在线| 欧美激情久久久久久爽电影 | av视频免费观看在线观看| 日韩电影二区| 欧美在线黄色| 91精品伊人久久大香线蕉| 国产亚洲精品久久久久5区| 日韩一卡2卡3卡4卡2021年| 免费在线观看日本一区| 人成视频在线观看免费观看| 男女免费视频国产| 在线观看免费日韩欧美大片| 最黄视频免费看| 精品熟女少妇八av免费久了| av片东京热男人的天堂| www.精华液| 欧美一级毛片孕妇| 一级片'在线观看视频| 免费av中文字幕在线| 人妻久久中文字幕网| 亚洲avbb在线观看| 黑人猛操日本美女一级片| 丰满人妻熟妇乱又伦精品不卡| 日本一区二区免费在线视频| 各种免费的搞黄视频| 免费久久久久久久精品成人欧美视频| 欧美国产精品一级二级三级| 久久久久久久精品精品| 亚洲国产欧美网| 天天躁日日躁夜夜躁夜夜| 亚洲专区字幕在线| 亚洲国产精品成人久久小说| 精品少妇一区二区三区视频日本电影| 一区二区三区激情视频| 超碰97精品在线观看| 天天躁夜夜躁狠狠躁躁| 欧美老熟妇乱子伦牲交| 高清黄色对白视频在线免费看| 亚洲激情五月婷婷啪啪| 亚洲色图 男人天堂 中文字幕| 老司机午夜十八禁免费视频| 欧美大码av| 国产免费视频播放在线视频| 亚洲自偷自拍图片 自拍| 欧美日韩黄片免| 精品一区二区三卡| 午夜福利在线免费观看网站| 后天国语完整版免费观看| 亚洲第一欧美日韩一区二区三区 | 国产一区二区三区综合在线观看| 香蕉丝袜av| 亚洲熟女精品中文字幕| 一级片'在线观看视频| 青春草亚洲视频在线观看| 91精品伊人久久大香线蕉| 麻豆乱淫一区二区| av线在线观看网站| 国产精品久久久人人做人人爽| 欧美日韩亚洲综合一区二区三区_| 亚洲精品久久久久久婷婷小说| e午夜精品久久久久久久| 两个人看的免费小视频| 色老头精品视频在线观看| 亚洲av成人不卡在线观看播放网 | 久久久久久亚洲精品国产蜜桃av| 国产精品 欧美亚洲| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品国产av蜜桃| 一区二区日韩欧美中文字幕| 亚洲九九香蕉| 精品国产一区二区久久| 看免费av毛片| 狂野欧美激情性bbbbbb| 极品少妇高潮喷水抽搐| www.999成人在线观看| 久久久国产精品麻豆| 视频在线观看一区二区三区| 日本vs欧美在线观看视频| 久久久久久久久免费视频了| 亚洲午夜精品一区,二区,三区| 中文字幕制服av| 国产不卡av网站在线观看| 纵有疾风起免费观看全集完整版| 成人国产av品久久久| 国产成人av教育| 亚洲精品av麻豆狂野| 蜜桃国产av成人99| 国产精品久久久人人做人人爽| 日本a在线网址| 桃花免费在线播放| 99久久99久久久精品蜜桃| 在线亚洲精品国产二区图片欧美| 91精品三级在线观看| 嫁个100分男人电影在线观看| 777米奇影视久久| 国产一卡二卡三卡精品| 9热在线视频观看99| 999久久久精品免费观看国产| 久久久欧美国产精品| 电影成人av| 老司机午夜十八禁免费视频| 中文字幕最新亚洲高清| 国产精品麻豆人妻色哟哟久久| 久久久久精品人妻al黑| 肉色欧美久久久久久久蜜桃| 国产国语露脸激情在线看| 99热网站在线观看| 国产精品秋霞免费鲁丝片| 青青草视频在线视频观看| 丁香六月欧美| 日韩一区二区三区影片| 日韩精品免费视频一区二区三区| 91成年电影在线观看| 汤姆久久久久久久影院中文字幕| 最新在线观看一区二区三区| 久久久久久久久免费视频了| 无限看片的www在线观看| 成年av动漫网址| www.999成人在线观看| 国产精品成人在线| 一本色道久久久久久精品综合| 亚洲欧美日韩高清在线视频 | 午夜精品久久久久久毛片777| 男女国产视频网站| av免费在线观看网站| 天天影视国产精品| 欧美精品人与动牲交sv欧美| 777米奇影视久久| 亚洲一区中文字幕在线| 青春草视频在线免费观看| 精品福利永久在线观看| 亚洲avbb在线观看| 久热爱精品视频在线9| 一区二区日韩欧美中文字幕| 视频在线观看一区二区三区| 亚洲精品乱久久久久久| 老司机亚洲免费影院| 青春草视频在线免费观看| 亚洲欧美一区二区三区黑人| 午夜久久久在线观看| 97人妻天天添夜夜摸| 亚洲国产精品999| 亚洲人成电影观看| 久久国产精品人妻蜜桃| 久久久久久免费高清国产稀缺| 久久狼人影院| 国产成人精品无人区| 十分钟在线观看高清视频www| 欧美人与性动交α欧美软件| 久久久欧美国产精品| 中文欧美无线码| 91精品伊人久久大香线蕉| 日本猛色少妇xxxxx猛交久久| 日韩有码中文字幕| 一级毛片女人18水好多| 国产亚洲av高清不卡| 亚洲第一青青草原| 日本av免费视频播放| 国产av又大| 美女大奶头黄色视频| 国产精品一区二区在线不卡| 老司机靠b影院| 成在线人永久免费视频| 99热国产这里只有精品6| 99国产精品99久久久久| 亚洲 国产 在线| kizo精华| 大香蕉久久网| 在线av久久热| 不卡av一区二区三区| av在线app专区| 亚洲国产精品一区二区三区在线| 免费不卡黄色视频| 免费高清在线观看视频在线观看| 久久午夜综合久久蜜桃| 少妇 在线观看| 亚洲精品久久成人aⅴ小说| 老鸭窝网址在线观看| 欧美精品亚洲一区二区| 午夜影院在线不卡| 中文字幕人妻丝袜一区二区| 国产精品二区激情视频| 亚洲国产欧美一区二区综合| 性色av一级| 男人添女人高潮全过程视频| 啪啪无遮挡十八禁网站| av天堂久久9| 久久国产精品影院| 美女视频免费永久观看网站| 久久人妻福利社区极品人妻图片| 久久久久国产一级毛片高清牌| 又紧又爽又黄一区二区| 人人妻人人爽人人添夜夜欢视频| av片东京热男人的天堂| 免费在线观看视频国产中文字幕亚洲 | 中文字幕高清在线视频| 搡老岳熟女国产| 亚洲黑人精品在线| 人人澡人人妻人| 精品人妻1区二区| av网站免费在线观看视频| 国产人伦9x9x在线观看| 国产一级毛片在线| 黄色视频在线播放观看不卡| 亚洲av成人不卡在线观看播放网 | 亚洲精品自拍成人| 青春草亚洲视频在线观看| 亚洲欧美激情在线| 一二三四社区在线视频社区8| 99热网站在线观看| 不卡av一区二区三区| 高清视频免费观看一区二区| 免费人妻精品一区二区三区视频| 50天的宝宝边吃奶边哭怎么回事| 悠悠久久av| 国产一卡二卡三卡精品| 国产熟女午夜一区二区三区| 精品少妇内射三级| 日韩中文字幕欧美一区二区| 国产精品二区激情视频| 性少妇av在线| 欧美 亚洲 国产 日韩一| 宅男免费午夜| 最黄视频免费看| netflix在线观看网站| 日韩电影二区| 久久久久国内视频| 午夜91福利影院| 国产日韩欧美亚洲二区| a级毛片在线看网站| 久久综合国产亚洲精品| 丝袜喷水一区| 1024视频免费在线观看| 久久久久久久久免费视频了| 午夜免费成人在线视频| 亚洲七黄色美女视频| 亚洲情色 制服丝袜| 黑人巨大精品欧美一区二区蜜桃| 91av网站免费观看| 高潮久久久久久久久久久不卡| 涩涩av久久男人的天堂| 青草久久国产| 精品国产国语对白av| 国产又色又爽无遮挡免| 最新的欧美精品一区二区| 国产激情久久老熟女| 91精品国产国语对白视频| 欧美人与性动交α欧美软件| 97人妻天天添夜夜摸| 午夜激情av网站| 黄色 视频免费看| 成人av一区二区三区在线看 | 久久精品亚洲av国产电影网| 亚洲av成人一区二区三| 午夜福利,免费看| 日韩有码中文字幕| 老熟妇乱子伦视频在线观看 | 国产一区二区三区在线臀色熟女 | 精品少妇黑人巨大在线播放| 首页视频小说图片口味搜索| 日韩一区二区三区影片| 黑丝袜美女国产一区| 在线av久久热| 亚洲欧美精品综合一区二区三区| 国产欧美日韩精品亚洲av| 女人精品久久久久毛片| 性色av一级| 中文字幕另类日韩欧美亚洲嫩草| 午夜免费鲁丝| 亚洲国产看品久久| 香蕉国产在线看| 日日爽夜夜爽网站| 国产免费av片在线观看野外av| 久久精品久久久久久噜噜老黄| 国产亚洲精品一区二区www | 啪啪无遮挡十八禁网站| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av香蕉五月 | 午夜日韩欧美国产| 欧美黄色片欧美黄色片| 日韩大片免费观看网站| 欧美+亚洲+日韩+国产| 亚洲国产精品一区二区三区在线| av在线播放精品| 日韩欧美一区视频在线观看| 最黄视频免费看| 一区二区三区激情视频| 亚洲欧美精品自产自拍| 窝窝影院91人妻| 欧美变态另类bdsm刘玥| 日韩制服骚丝袜av| 午夜福利影视在线免费观看| 人妻一区二区av| 国产精品免费大片| 国产欧美日韩一区二区三区在线| 青草久久国产| 日韩一区二区三区影片| 成在线人永久免费视频| 久久99热这里只频精品6学生| 午夜福利乱码中文字幕| 日韩大片免费观看网站| 黄网站色视频无遮挡免费观看| 国产福利在线免费观看视频| 午夜福利视频精品| 大香蕉久久网| 国产区一区二久久| 久久人人97超碰香蕉20202| 久久精品成人免费网站| 精品福利观看| 51午夜福利影视在线观看| 国产精品熟女久久久久浪| 国产精品.久久久| 成人影院久久| 精品一品国产午夜福利视频| 在线观看免费日韩欧美大片| 法律面前人人平等表现在哪些方面 | 精品第一国产精品| 欧美午夜高清在线| 高清av免费在线| av欧美777| www.熟女人妻精品国产| 亚洲精品一二三| 国产成人欧美在线观看 | 亚洲午夜精品一区,二区,三区| 日韩欧美一区二区三区在线观看 | 成年女人毛片免费观看观看9 | 国产高清视频在线播放一区 | 青春草亚洲视频在线观看| 日韩电影二区| 国产成人啪精品午夜网站| 男人添女人高潮全过程视频| 99久久精品国产亚洲精品| 男人舔女人的私密视频| 久久天躁狠狠躁夜夜2o2o| av电影中文网址| 十八禁网站网址无遮挡| 黑人猛操日本美女一级片| 国产亚洲av片在线观看秒播厂| 精品一区二区三区四区五区乱码| 亚洲成人免费电影在线观看| 日韩大码丰满熟妇| 一边摸一边抽搐一进一出视频| 日韩欧美免费精品| 丁香六月欧美| av有码第一页| 久久影院123| 久久久欧美国产精品| 婷婷丁香在线五月| 一个人免费看片子| 亚洲欧美成人综合另类久久久| 久久久久久人人人人人| 丝袜人妻中文字幕| 久久人人爽人人片av| 最近最新免费中文字幕在线| 少妇 在线观看| 高清黄色对白视频在线免费看| 在线十欧美十亚洲十日本专区| 俄罗斯特黄特色一大片| 91精品伊人久久大香线蕉| 亚洲成av片中文字幕在线观看| 伦理电影免费视频| 一区二区日韩欧美中文字幕| 人妻久久中文字幕网| 国产成人系列免费观看|