• <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ū)代表性觀測點地下水位變化特征分析
    基于升降溫全曲線的鋼筋混凝土梁溫度場分析
    亚洲无线观看免费| 久久久a久久爽久久v久久| 夜夜夜夜夜久久久久| 免费看日本二区| 精品久久久久久成人av| 国产91av在线免费观看| 欧美成人免费av一区二区三区| 午夜福利成人在线免费观看| 国产综合懂色| 免费av不卡在线播放| 久久国产乱子免费精品| 乱人视频在线观看| 日韩 亚洲 欧美在线| 黄色欧美视频在线观看| 麻豆久久精品国产亚洲av| 亚洲av.av天堂| 日韩av在线大香蕉| 亚洲图色成人| 小说图片视频综合网站| 在线免费观看的www视频| 天天躁日日操中文字幕| 欧美极品一区二区三区四区| 久久久久久久久久成人| 黄片wwwwww| 少妇裸体淫交视频免费看高清| 在线免费观看不下载黄p国产| 黄色欧美视频在线观看| 91午夜精品亚洲一区二区三区| 一边亲一边摸免费视频| 一进一出抽搐动态| 亚洲性久久影院| 国产精品久久久久久亚洲av鲁大| 乱人视频在线观看| 亚洲av电影不卡..在线观看| 亚洲经典国产精华液单| 久久中文看片网| 午夜福利在线观看吧| 欧美极品一区二区三区四区| 人人妻人人澡人人爽人人夜夜 | 亚洲美女搞黄在线观看| 全区人妻精品视频| 又粗又硬又长又爽又黄的视频 | 一个人看视频在线观看www免费| av天堂在线播放| 九九热线精品视视频播放| 97超碰精品成人国产| 三级男女做爰猛烈吃奶摸视频| 欧美日韩精品成人综合77777| 超碰av人人做人人爽久久| 久久精品国产亚洲网站| 久久精品国产亚洲av涩爱 | 精品日产1卡2卡| 内射极品少妇av片p| 村上凉子中文字幕在线| 亚洲国产精品成人综合色| 如何舔出高潮| 亚洲经典国产精华液单| 精品午夜福利在线看| 给我免费播放毛片高清在线观看| 一本精品99久久精品77| 国产免费一级a男人的天堂| 亚洲熟妇中文字幕五十中出| 欧美xxxx性猛交bbbb| 国产在线男女| 亚洲欧美精品专区久久| 国产黄片视频在线免费观看| 蜜桃久久精品国产亚洲av| 国产毛片a区久久久久| 少妇丰满av| 亚洲aⅴ乱码一区二区在线播放| 成人亚洲精品av一区二区| 久久久久国产网址| 久久久久久伊人网av| 国产老妇伦熟女老妇高清| 高清在线视频一区二区三区 | 精品人妻偷拍中文字幕| 亚洲内射少妇av| 老司机影院成人| 成人综合一区亚洲| 日产精品乱码卡一卡2卡三| 中国国产av一级| 成年免费大片在线观看| 黄色一级大片看看| 免费电影在线观看免费观看| 精品人妻熟女av久视频| 成人永久免费在线观看视频| 国产伦精品一区二区三区视频9| 美女内射精品一级片tv| 天堂网av新在线| 我的老师免费观看完整版| 少妇的逼好多水| 午夜激情欧美在线| 老熟妇乱子伦视频在线观看| 日韩三级伦理在线观看| 人妻系列 视频| 国产伦理片在线播放av一区 | 最近2019中文字幕mv第一页| 欧美激情国产日韩精品一区| 成人国产麻豆网| 午夜老司机福利剧场| 久久久a久久爽久久v久久| 寂寞人妻少妇视频99o| АⅤ资源中文在线天堂| 少妇人妻精品综合一区二区 | 亚洲内射少妇av| 日韩 亚洲 欧美在线| 亚洲成av人片在线播放无| 只有这里有精品99| 嫩草影院入口| 联通29元200g的流量卡| 国产真实伦视频高清在线观看| 校园人妻丝袜中文字幕| 中文在线观看免费www的网站| 日本-黄色视频高清免费观看| 久久人人爽人人爽人人片va| 麻豆av噜噜一区二区三区| 亚洲在久久综合| 国产不卡一卡二| 最后的刺客免费高清国语| 精品久久久久久久久久免费视频| 2022亚洲国产成人精品| 国产黄片视频在线免费观看| 亚洲人成网站高清观看| 久久99蜜桃精品久久| 欧美zozozo另类| 精品久久久久久成人av| 色5月婷婷丁香| 99久久九九国产精品国产免费| 一个人看视频在线观看www免费| 久久久久国产网址| 1000部很黄的大片| 97超视频在线观看视频| 69av精品久久久久久| 午夜福利视频1000在线观看| 久久精品国产99精品国产亚洲性色| 大又大粗又爽又黄少妇毛片口| 伦理电影大哥的女人| 亚洲av中文av极速乱| 黑人高潮一二区| 国产精品久久久久久av不卡| 亚洲成a人片在线一区二区| 国产成人精品一,二区 | 我要搜黄色片| 精品人妻熟女av久视频| 国产久久久一区二区三区| 国产亚洲精品av在线| 国产一区二区在线观看日韩| 在线免费十八禁| 国产伦一二天堂av在线观看| 午夜a级毛片| 国内精品宾馆在线| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 极品教师在线视频| 国产熟女欧美一区二区| a级毛片a级免费在线| 国产精品久久久久久久久免| 三级国产精品欧美在线观看| 久久6这里有精品| 国产色爽女视频免费观看| 久久久久久九九精品二区国产| 综合色丁香网| 国产日韩欧美在线精品| 日本色播在线视频| 一本久久中文字幕| 国产成人freesex在线| 久久99热这里只有精品18| 老司机影院成人| 久久久久久久久久成人| 天堂√8在线中文| 亚洲成a人片在线一区二区| 国产白丝娇喘喷水9色精品| 亚州av有码| 久久久久性生活片| 高清午夜精品一区二区三区 | 欧美不卡视频在线免费观看| 男女边吃奶边做爰视频| 直男gayav资源| 婷婷精品国产亚洲av| 国产国拍精品亚洲av在线观看| a级毛色黄片| av天堂中文字幕网| 欧美性感艳星| 国产日韩欧美在线精品| 成人国产麻豆网| 国产精品久久久久久精品电影小说 | 观看免费一级毛片| 亚洲欧美日韩无卡精品| 男人和女人高潮做爰伦理| 亚洲性久久影院| 我的女老师完整版在线观看| 极品教师在线视频| 男人舔奶头视频| 亚洲国产精品成人综合色| 亚洲三级黄色毛片| 草草在线视频免费看| 可以在线观看的亚洲视频| 欧美激情久久久久久爽电影| 精品久久久久久久久av| 亚洲第一电影网av| 97人妻精品一区二区三区麻豆| 最后的刺客免费高清国语| 免费av观看视频| 久久久久久伊人网av| 免费在线观看成人毛片| 天堂中文最新版在线下载 | 国产激情偷乱视频一区二区| 国产成人91sexporn| 最近的中文字幕免费完整| 国产精品电影一区二区三区| 国产三级在线视频| 久久国内精品自在自线图片| 国产精品电影一区二区三区| 欧美日韩在线观看h| 亚洲欧美清纯卡通| 伦理电影大哥的女人| 男女边吃奶边做爰视频| 久久久久久久亚洲中文字幕| 男女做爰动态图高潮gif福利片| 国产人妻一区二区三区在| 黄色一级大片看看| 久久综合国产亚洲精品| 99久久精品一区二区三区| 热99在线观看视频| 婷婷色av中文字幕| 中文字幕av成人在线电影| 免费人成视频x8x8入口观看| 日日摸夜夜添夜夜添av毛片| 国产国拍精品亚洲av在线观看| 成人亚洲精品av一区二区| 国内少妇人妻偷人精品xxx网站| 99国产精品一区二区蜜桃av| 国产午夜精品论理片| 国产成年人精品一区二区| 中出人妻视频一区二区| av天堂在线播放| 最近手机中文字幕大全| 亚洲色图av天堂| 亚洲无线在线观看| 亚洲欧美日韩高清在线视频| 日韩制服骚丝袜av| 嘟嘟电影网在线观看| 日本黄大片高清| а√天堂www在线а√下载| 亚洲三级黄色毛片| 99热只有精品国产| 久久精品国产亚洲av天美| 2022亚洲国产成人精品| 亚洲成人久久爱视频| 夫妻性生交免费视频一级片| 99国产精品一区二区蜜桃av| 欧美+日韩+精品| 亚洲国产欧洲综合997久久,| 日韩欧美国产在线观看| 一区福利在线观看| 如何舔出高潮| 亚洲欧美日韩无卡精品| 青春草国产在线视频 | 日韩av在线大香蕉| 不卡视频在线观看欧美| 国产精品女同一区二区软件| 国产成人a区在线观看| 一区福利在线观看| 欧美成人免费av一区二区三区| 91午夜精品亚洲一区二区三区| 插阴视频在线观看视频| 欧美最新免费一区二区三区| 欧美zozozo另类| 免费av观看视频| 亚洲国产欧美在线一区| 最近最新中文字幕大全电影3| 中文欧美无线码| 国产精品一区www在线观看| 亚洲人成网站在线播放欧美日韩| 免费搜索国产男女视频| 黄色配什么色好看| 欧美日韩乱码在线| 欧美日本视频| 欧美一区二区国产精品久久精品| 国产探花在线观看一区二区| 爱豆传媒免费全集在线观看| 日本黄色片子视频| 久久久久久久久久成人| 国产女主播在线喷水免费视频网站 | 日韩欧美精品免费久久| 午夜久久久久精精品| 免费av不卡在线播放| 国产精品,欧美在线| 国产老妇伦熟女老妇高清| 欧美+日韩+精品| 精品不卡国产一区二区三区| 深夜精品福利| 久久久久性生活片| 黄色配什么色好看| 精品免费久久久久久久清纯| 91精品一卡2卡3卡4卡| 国产国拍精品亚洲av在线观看| 精品久久久久久久久久久久久| av免费观看日本| 在线免费观看不下载黄p国产| 身体一侧抽搐| 好男人在线观看高清免费视频| 久久精品国产亚洲av涩爱 | 少妇熟女欧美另类| 中文字幕人妻熟人妻熟丝袜美| 在线播放国产精品三级| 色综合色国产| 乱人视频在线观看| 日韩视频在线欧美| 亚洲欧美成人综合另类久久久 | 青春草国产在线视频 | 69av精品久久久久久| 午夜免费男女啪啪视频观看| 亚洲精品粉嫩美女一区| 啦啦啦观看免费观看视频高清| 久久精品国产亚洲av香蕉五月| 麻豆精品久久久久久蜜桃| 三级国产精品欧美在线观看| 亚洲精品国产av成人精品| 岛国毛片在线播放| 综合色丁香网| 99久久成人亚洲精品观看| 自拍偷自拍亚洲精品老妇| 在线免费十八禁| 亚洲人成网站在线播放欧美日韩| 美女内射精品一级片tv| 国产白丝娇喘喷水9色精品| 青青草视频在线视频观看| 免费人成在线观看视频色| 日韩欧美一区二区三区在线观看| 嘟嘟电影网在线观看| 日日啪夜夜撸| 国产黄色视频一区二区在线观看 | 精华霜和精华液先用哪个| 国产麻豆成人av免费视频| 女同久久另类99精品国产91| 久久韩国三级中文字幕| 国内精品一区二区在线观看| 青青草视频在线视频观看| 亚洲天堂国产精品一区在线| 精品国产三级普通话版| 91精品国产九色| 久久99精品国语久久久| 欧美色欧美亚洲另类二区| 插逼视频在线观看| 国产精品久久久久久av不卡| h日本视频在线播放| 我要搜黄色片| 亚洲欧美成人精品一区二区| 久久久久久国产a免费观看| av在线老鸭窝| 99热网站在线观看| 国产av一区在线观看免费| 26uuu在线亚洲综合色| 99热全是精品| 男女啪啪激烈高潮av片| 1024手机看黄色片| 欧美日韩一区二区视频在线观看视频在线 | 熟女人妻精品中文字幕| 悠悠久久av| 亚洲欧美清纯卡通| 免费观看在线日韩| 国产精华一区二区三区| 在线天堂最新版资源| 少妇猛男粗大的猛烈进出视频 | 天天躁日日操中文字幕| 欧美不卡视频在线免费观看| 99久久无色码亚洲精品果冻| 伦理电影大哥的女人| 欧美日韩乱码在线| 亚洲欧美日韩东京热| 18禁在线播放成人免费| 亚洲无线观看免费| 日本免费a在线| 美女被艹到高潮喷水动态| 草草在线视频免费看| 特大巨黑吊av在线直播| 3wmmmm亚洲av在线观看| 国产精品久久久久久久久免| 亚洲乱码一区二区免费版| 中文字幕精品亚洲无线码一区| 国产一级毛片七仙女欲春2| 18禁在线无遮挡免费观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 色综合站精品国产| 可以在线观看毛片的网站| 变态另类成人亚洲欧美熟女| 国产老妇女一区| 国产高潮美女av| 久久人人精品亚洲av| 少妇人妻一区二区三区视频| 人体艺术视频欧美日本| 日本色播在线视频| 久久久久久国产a免费观看| 卡戴珊不雅视频在线播放| 一边摸一边抽搐一进一小说| 在线免费观看不下载黄p国产| 麻豆av噜噜一区二区三区| 国产精品美女特级片免费视频播放器| 两个人视频免费观看高清| 国产不卡一卡二| 嫩草影院新地址| 能在线免费看毛片的网站| 美女cb高潮喷水在线观看| 又粗又硬又长又爽又黄的视频 | 国产成人精品久久久久久| 国内精品美女久久久久久| 国产熟女欧美一区二区| ponron亚洲| 日本色播在线视频| 精品午夜福利在线看| 日本爱情动作片www.在线观看| 日本黄色片子视频| 少妇的逼水好多| 一级黄色大片毛片| av在线观看视频网站免费| 免费电影在线观看免费观看| 成人特级av手机在线观看| 日本一二三区视频观看| 欧美极品一区二区三区四区| 在线播放无遮挡| 97在线视频观看| 国产成人a区在线观看| av又黄又爽大尺度在线免费看 | 国产成人a区在线观看| 国产私拍福利视频在线观看| 精品久久久久久久久久久久久| 欧美bdsm另类| 亚洲美女搞黄在线观看| 18+在线观看网站| 国产白丝娇喘喷水9色精品| 国产精品一区二区在线观看99 | 18禁在线播放成人免费| 亚洲一级一片aⅴ在线观看| 亚洲av不卡在线观看| 亚洲精品成人久久久久久| www.av在线官网国产| 亚洲成a人片在线一区二区| 久久人人爽人人片av| 日韩精品有码人妻一区| 亚洲一区二区三区色噜噜| 精品人妻视频免费看| 亚洲av.av天堂| 欧美bdsm另类| 国产精品.久久久| 看十八女毛片水多多多| 国产精品99久久久久久久久| 边亲边吃奶的免费视频| 久久久精品欧美日韩精品| 观看免费一级毛片| 亚洲欧美日韩无卡精品| 国产精品久久久久久精品电影小说 | 少妇丰满av| 最近的中文字幕免费完整| 在线a可以看的网站| 青青草视频在线视频观看| 精品久久久久久久久久免费视频| 哪个播放器可以免费观看大片| 麻豆成人午夜福利视频| 成人高潮视频无遮挡免费网站| 日本免费a在线| 国产精品.久久久| 插阴视频在线观看视频| 1000部很黄的大片| 51国产日韩欧美| 国产真实伦视频高清在线观看| 久久草成人影院| 欧美又色又爽又黄视频| 热99re8久久精品国产| 亚洲第一区二区三区不卡| av在线播放精品| 亚洲在线自拍视频| 日本与韩国留学比较| 亚洲aⅴ乱码一区二区在线播放| av黄色大香蕉| 中文字幕av成人在线电影| 一个人免费在线观看电影| 人妻系列 视频| 亚洲色图av天堂| 亚洲国产欧美人成| 亚洲国产精品国产精品| 亚洲欧美日韩无卡精品| 国产精品美女特级片免费视频播放器| 人体艺术视频欧美日本| 天天躁日日操中文字幕| 少妇的逼水好多| 国产精品一区二区三区四区久久| 欧美另类亚洲清纯唯美| 日本五十路高清| 麻豆国产av国片精品| 三级男女做爰猛烈吃奶摸视频| 国产黄片视频在线免费观看| 亚洲人成网站在线观看播放| 亚洲自偷自拍三级| 亚洲图色成人| 亚洲美女视频黄频| 女人被狂操c到高潮| 久久欧美精品欧美久久欧美| 久久久国产成人免费| 美女脱内裤让男人舔精品视频 | 亚洲四区av| 亚洲欧美日韩东京热| 高清日韩中文字幕在线| 国产成人a∨麻豆精品| 久久久国产成人精品二区| 尾随美女入室| 国产成人精品婷婷| 晚上一个人看的免费电影| 1000部很黄的大片| 给我免费播放毛片高清在线观看| 国产精品一区二区三区四区免费观看| 麻豆国产97在线/欧美| 亚洲国产精品久久男人天堂| 国产三级中文精品| 插阴视频在线观看视频| 91精品一卡2卡3卡4卡| 99在线人妻在线中文字幕| 99热这里只有是精品在线观看| 黄色日韩在线| 久久综合国产亚洲精品| 联通29元200g的流量卡| av免费在线看不卡| 99久久人妻综合| 久久久成人免费电影| 女的被弄到高潮叫床怎么办| 欧美+亚洲+日韩+国产| 久久九九热精品免费| 夜夜看夜夜爽夜夜摸| 久久韩国三级中文字幕| 99精品在免费线老司机午夜| av国产免费在线观看| 老司机影院成人| 国语自产精品视频在线第100页| 成人漫画全彩无遮挡| 亚洲久久久久久中文字幕| 午夜福利在线在线| 国产亚洲av嫩草精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 男女边吃奶边做爰视频| 联通29元200g的流量卡| 久久亚洲精品不卡| 亚洲av不卡在线观看| 久久精品夜夜夜夜夜久久蜜豆| 男人狂女人下面高潮的视频| 久久人人精品亚洲av| 欧美色欧美亚洲另类二区| 国产女主播在线喷水免费视频网站 | a级毛片免费高清观看在线播放| 精品久久久久久成人av| 少妇人妻一区二区三区视频| 亚洲成人久久爱视频| 亚洲成人av在线免费| 一本久久精品| 国产精品1区2区在线观看.| 久久韩国三级中文字幕| 99久久九九国产精品国产免费| 一区二区三区四区激情视频 | 成人特级av手机在线观看| 男女那种视频在线观看| 少妇被粗大猛烈的视频| 午夜福利成人在线免费观看| 在线观看一区二区三区| 国产单亲对白刺激| 久久人人精品亚洲av| 日韩一区二区三区影片| 又粗又爽又猛毛片免费看| 麻豆国产av国片精品| 国产乱人视频| 一级av片app| 亚洲内射少妇av| 欧美不卡视频在线免费观看| 能在线免费观看的黄片| 亚洲美女视频黄频| 人人妻人人看人人澡| 欧美在线一区亚洲| 国产精品综合久久久久久久免费| 校园春色视频在线观看| 哪里可以看免费的av片| 99久久九九国产精品国产免费| 熟女人妻精品中文字幕| 午夜a级毛片| 欧美高清成人免费视频www| 69av精品久久久久久| 22中文网久久字幕| 中文字幕人妻熟人妻熟丝袜美| 亚洲最大成人av| 高清午夜精品一区二区三区 | 久久精品国产清高在天天线| 97超碰精品成人国产| 欧美成人一区二区免费高清观看| 男的添女的下面高潮视频| 99久久人妻综合| 亚洲无线在线观看| 国产亚洲精品久久久久久毛片| 久久精品国产99精品国产亚洲性色| 国产精品国产三级国产av玫瑰| 久久鲁丝午夜福利片| 可以在线观看的亚洲视频| 日韩精品有码人妻一区| 国产亚洲精品av在线| 亚洲va在线va天堂va国产| 成人三级黄色视频| 成人一区二区视频在线观看| 日日干狠狠操夜夜爽| 亚洲av中文av极速乱| 级片在线观看| av在线播放精品| 男女啪啪激烈高潮av片| 乱码一卡2卡4卡精品| 69人妻影院| 边亲边吃奶的免费视频| 日韩欧美国产在线观看| 欧美一区二区亚洲| 国产成人精品久久久久久| 九九爱精品视频在线观看|