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

    中國(guó)地震安全性評(píng)價(jià)中天然強(qiáng)震記錄選取

    2016-12-22 07:05:39溫瑞智任葉飛

    冀 昆,溫瑞智,任葉飛

    (中國(guó)地震局工程力學(xué)研究所 地震工程與工程振動(dòng)重點(diǎn)實(shí)驗(yàn)室, 哈爾濱 150080)

    ?

    中國(guó)地震安全性評(píng)價(jià)中天然強(qiáng)震記錄選取

    冀 昆,溫瑞智,任葉飛

    (中國(guó)地震局工程力學(xué)研究所 地震工程與工程振動(dòng)重點(diǎn)實(shí)驗(yàn)室, 哈爾濱 150080)

    為解決目前中國(guó)地震安全性評(píng)價(jià)中天然強(qiáng)震記錄選取的問(wèn)題,基于條件均值目標(biāo)譜的概念,從中國(guó)地震安全性評(píng)價(jià)的已有結(jié)果出發(fā),同時(shí)考慮震級(jí)、距離和地震動(dòng)離散性,對(duì)危險(xiǎn)性分析結(jié)果進(jìn)行設(shè)定地震解耦,構(gòu)建了反映結(jié)構(gòu)特征和地震動(dòng)特征的條件均值目標(biāo)譜,為天然強(qiáng)震記錄的選取提供匹配依據(jù).安全性評(píng)價(jià)實(shí)例的強(qiáng)震記錄結(jié)果表明,本文建議方法在保證選取物理意義與設(shè)定地震相符的前提下較好兼顧了全周期段的譜形匹配,均值偏差和離散性水平均較低.最后給出適用于中國(guó)地震安全性評(píng)價(jià)工作的天然強(qiáng)震記錄選取建議流程,所得結(jié)果可為結(jié)構(gòu)時(shí)程分析與抗震驗(yàn)算提供參考.關(guān)鍵詞: 安全性評(píng)價(jià);天然強(qiáng)震記錄;設(shè)定地震;條件均值譜;強(qiáng)震記錄選取

    地震安全性評(píng)價(jià)是根據(jù)建設(shè)工程場(chǎng)址周?chē)卣瓠h(huán)境調(diào)查、地質(zhì)條件勘測(cè)的結(jié)果,綜合評(píng)價(jià)和分析計(jì)算,最終按照工程類(lèi)型、性質(zhì)、重要性,科學(xué)合理地給出符合工程抗震設(shè)防要求的相應(yīng)地震動(dòng)參數(shù).中國(guó)相關(guān)標(biāo)準(zhǔn)規(guī)定對(duì)于重大建設(shè)工程,必須進(jìn)行地震安全性評(píng)價(jià),并根據(jù)該結(jié)果指導(dǎo)抗震設(shè)防.目前中國(guó)地震安全性評(píng)價(jià)工作一般以基于一致概率譜(UHS) 的人造地震動(dòng)作為基巖面輸入,采用一維土層等效線性化方法得到考慮場(chǎng)地反應(yīng)的地面人造地震動(dòng).但人造地震動(dòng)合成中包絡(luò)參數(shù)的確定和擬合尚沒(méi)有具體可靠的規(guī)定,這種不確定性導(dǎo)致地震安全性評(píng)價(jià)中給出不合理人造地震動(dòng)時(shí)程的問(wèn)題時(shí)有出現(xiàn)[1].而隨著結(jié)構(gòu)動(dòng)力時(shí)程分析的應(yīng)用日漸廣泛,選取天然強(qiáng)震記錄的需求變得越來(lái)越迫切.雖然現(xiàn)行規(guī)范GB 17741—2005《工程場(chǎng)地地震安全性評(píng)價(jià)》規(guī)定:“本地有強(qiáng)震動(dòng)記錄時(shí),宜充分利用其合成適合工程場(chǎng)地的基巖地震動(dòng)時(shí)程”,但由于地震動(dòng)本身的不可重復(fù)性與隨機(jī)性,在待建工程場(chǎng)地附近很難找到滿足設(shè)防要求的天然強(qiáng)震記錄,目前中國(guó)工程場(chǎng)地安全性評(píng)價(jià)工作中鮮有提供天然強(qiáng)震記錄.

    隨著國(guó)內(nèi)外強(qiáng)震記錄不斷積累和完善,為研究中國(guó)工程場(chǎng)地地震安全性評(píng)價(jià)中強(qiáng)震記錄選取提供了重要的基礎(chǔ)數(shù)據(jù).現(xiàn)有的天然強(qiáng)震記錄選取及修改方法眾多,但基本思路一致:首先確定滿足設(shè)定地震參數(shù)的目標(biāo)譜,進(jìn)而在強(qiáng)震記錄數(shù)據(jù)庫(kù)中尋找與目標(biāo)譜匹配程度滿足要求的結(jié)果.目標(biāo)譜是直接決定記錄選取結(jié)果是否合理的關(guān)鍵因素[2].不少學(xué)者已經(jīng)認(rèn)識(shí)到安全性評(píng)價(jià)中的一致概率譜是若干設(shè)定地震事件綜合后得到的包絡(luò)結(jié)果,并無(wú)法代表真實(shí)的地震動(dòng)譜形,并不適合指導(dǎo)強(qiáng)震記錄選取[3-4].針對(duì)上述目標(biāo)譜存在的問(wèn)題,2011年Baker[5]同時(shí)考慮設(shè)定地震信息和地震動(dòng)固有的離散性,提出了條件均值譜(CMS)概念:從目標(biāo)結(jié)構(gòu)特性出發(fā)將一致概率譜分解為若干條件均值目標(biāo)譜的組合.條件均值目標(biāo)譜具有更窄的譜帶寬,更明確的物理意義和更合理的形狀,近幾年在國(guó)外得到了廣泛實(shí)踐[6].國(guó)內(nèi)學(xué)者也驗(yàn)證了基于條件均值目標(biāo)譜的選波結(jié)果在評(píng)估長(zhǎng)周期結(jié)構(gòu)響應(yīng)離散性方面上更加合理[7].

    雖然條件均值目標(biāo)譜得到了不少學(xué)者的關(guān)注和青睞,但目前中國(guó)地震安全性評(píng)價(jià)工作中并沒(méi)有將地震動(dòng)離散性作為單獨(dú)參數(shù)予以考慮,因而該方法一直無(wú)法在中國(guó)實(shí)際應(yīng)用中推廣.為解決上述問(wèn)題,本文基于中國(guó)工程場(chǎng)地地震安全性評(píng)價(jià)結(jié)果給出設(shè)定地震信息解耦的思路方法;以某一工程場(chǎng)地作為案例,闡述解耦地震信息,構(gòu)建條件目標(biāo)譜和選取強(qiáng)震記錄的全過(guò)程;對(duì)比不同目標(biāo)譜選波結(jié)果差異后,給出適用于中國(guó)現(xiàn)行安評(píng)工作的天然強(qiáng)震記錄選取流程.

    1 設(shè)定地震解耦

    設(shè)定地震是條件均值目標(biāo)譜構(gòu)建以及強(qiáng)震記錄選取的基礎(chǔ).中國(guó)安全性評(píng)價(jià)工作給出的設(shè)定地震信息僅包含震級(jí)和震中距兩個(gè)參數(shù),主要用于標(biāo)定人造地震動(dòng)包絡(luò)函數(shù).實(shí)際觀測(cè)得到的強(qiáng)震記錄表明,即使同一震級(jí),同一震中距下的地震動(dòng)記錄,其反應(yīng)譜值都可能相差甚巨,這種地震動(dòng)自身的離散性會(huì)直接影響后面強(qiáng)震記錄選取和結(jié)構(gòu)響應(yīng)分析的結(jié)果.

    McGuire[8]于1995年最先提出了依據(jù)概率地震危險(xiǎn)性分析(PSHA)結(jié)果解耦設(shè)定地震的思路,在這一過(guò)程中引入了地震動(dòng)離散性參數(shù).所謂解耦,就是從危險(xiǎn)性分析的結(jié)果出發(fā),反推出同時(shí)滿足地震活動(dòng)性,潛源分布空間相容性的地震參數(shù)組合的貢獻(xiàn)比例.這種解耦的思路和結(jié)果已廣泛應(yīng)用于國(guó)外工程實(shí)踐,是其工程強(qiáng)震記錄選取的重要初選依據(jù).美國(guó)抗震設(shè)計(jì)規(guī)范(ASCE/SEI-7)[9]中已經(jīng)明確要求時(shí)程分析所選用的地震記錄應(yīng)以該解耦結(jié)果得到的最可信地震(MCE)為基礎(chǔ),歐洲抗震規(guī)范(Europe-8)[10]、新西蘭抗震規(guī)范(NZS)[11]等同樣規(guī)定在強(qiáng)震記錄選取時(shí)應(yīng)以該設(shè)定地震條件為初步篩選條件.

    1.1 標(biāo)準(zhǔn)差系數(shù)ε

    標(biāo)準(zhǔn)差系數(shù)ε定義為某設(shè)定地震下反應(yīng)譜目標(biāo)值與地震動(dòng)衰減關(guān)系均值之間的差值,用標(biāo)準(zhǔn)差的倍數(shù)度量,通常用來(lái)表征地震動(dòng)離散性.

    (1)

    式中:(M,R)為設(shè)定地震的震級(jí)-震中距組合,lnSa(T)為周期點(diǎn)T處的目標(biāo)反應(yīng)譜值,μlnSa(T)和σlnSa(T)為地震動(dòng)衰減關(guān)系在設(shè)定地震下的周期點(diǎn)T的均值和標(biāo)準(zhǔn)差.Baker等[12-13]認(rèn)為,標(biāo)準(zhǔn)差系數(shù)ε反映了目標(biāo)譜值與衰減關(guān)系中值在周期點(diǎn)T的差別程度,對(duì)地震記錄的反應(yīng)譜形狀具有重要影響.若ε為正,表明目標(biāo)周期處的目標(biāo)值要大于衰減關(guān)系中值,而其他周期點(diǎn)處的反應(yīng)譜值則要趨近于衰減關(guān)系中值,于是在周期點(diǎn)T形成一個(gè)“波峰”;反之,對(duì)于ε為負(fù)的反應(yīng)譜,會(huì)在T處形成“波谷”.

    1.2 基于PSHA結(jié)果解耦思路

    依據(jù)待建結(jié)構(gòu)特征和危險(xiǎn)性分析得到的一致概率譜,確定對(duì)應(yīng)的地震動(dòng)強(qiáng)度參數(shù)(IM, intensity measurement)目標(biāo)值IMtarget,一般選取目標(biāo)結(jié)構(gòu)自振周期T1處的一致概率譜值.該目標(biāo)值的超越概率為

    (2)

    式中:n為主要貢獻(xiàn)潛源的個(gè)數(shù);vi為第i個(gè)潛在震源區(qū)的地震年平均發(fā)生率,取決于該潛源所在的地震帶;fM(m)、fR(r)和fε(ε)分別為震級(jí)、震中距和標(biāo)準(zhǔn)差系數(shù)的概率密度函數(shù).

    將式(2)中變量離散化后的概率表達(dá)式為

    (3)

    式中:N1、N2、N3分別為震級(jí)、震中距和標(biāo)準(zhǔn)差系數(shù)的分檔數(shù);PM(mj)為某潛源內(nèi)震級(jí)為mj檔的概率[14],震級(jí)分檔時(shí)應(yīng)當(dāng)符合該潛源的震級(jí)上下限約束;將該潛源均勻離散為若干個(gè)微元后,各個(gè)距離檔的微元個(gè)數(shù)與微元總數(shù)之比即為某潛源內(nèi)震中距為rk檔的概率PR(rk);地震動(dòng)衰減關(guān)系的殘差可認(rèn)為服從對(duì)數(shù)正態(tài)分布,標(biāo)準(zhǔn)差系數(shù)服從標(biāo)準(zhǔn)正態(tài)分布,由正態(tài)分布概率密度函數(shù)積分即可得到標(biāo)準(zhǔn)差系數(shù)在εh檔的概率Pε(εh).

    P(IM≥IMtarget|mj,rk,εh)計(jì)算式為

    P(IM≥IMtarget|mj,rk,εh)=

    (4)

    根據(jù)貝葉斯全概率公式,目標(biāo)值條件下各個(gè)(mj,rk,εh)組合的條件概率P(mj,rk,εh|IM≥IMtarget)如式(5)所示,該條件概率即各設(shè)定地震對(duì)目標(biāo)值的貢獻(xiàn)比例,是解耦最終得到的結(jié)果.

    (5)

    上述分析過(guò)程表明,解耦所依據(jù)的目標(biāo)值可是任意一個(gè)特定周期點(diǎn)的加速度反應(yīng)譜值,具體的周期點(diǎn)由目標(biāo)結(jié)構(gòu)決定.

    1.3 工程實(shí)例分析

    通過(guò)一個(gè)工程場(chǎng)地安全性評(píng)價(jià)實(shí)例說(shuō)明設(shè)定地震解耦的過(guò)程(數(shù)據(jù)來(lái)源《港清三線輸氣管道工程地震安全性評(píng)價(jià)報(bào)告》).工程場(chǎng)址和所在區(qū)域主要潛源空間位置分布見(jiàn)圖1,涉及的潛源所屬地震帶均為華北平原地震帶,地震活動(dòng)參數(shù)a值為3.79,b值為0.83;地震動(dòng)衰減關(guān)系采用在中國(guó)地震安全性評(píng)價(jià)工程中廣泛使用的霍俊榮(1989)衰減關(guān)系[15].主要4個(gè)潛源對(duì)應(yīng)的震級(jí)上限以及空間分布函數(shù)列于表1,工程場(chǎng)址和所在區(qū)域主要潛源空間位置分布見(jiàn)圖1,安全性評(píng)價(jià)中通過(guò)概率危險(xiǎn)性分析得到的50 a超越概率2%、10%的一致概率譜見(jiàn)圖2.

    依據(jù)主要潛源的震級(jí)上下限和空間范圍,確定震級(jí)分檔范圍為[4.5~8],以0.5為分檔間隔,每檔震級(jí)的中值作為該檔震級(jí)的代表值;震中距r范圍為[1~100 km],以5 km作為分檔間隔;標(biāo)準(zhǔn)差系數(shù)ε區(qū)間為[-3~3].假設(shè)待建目標(biāo)結(jié)構(gòu)的自振周期為1.0 s,以加速度譜值Sa(T=1.0 s)為目標(biāo)值進(jìn)行設(shè)定地震解耦,見(jiàn)圖2.依據(jù)式(2)~(5)計(jì)算各個(gè)潛源下各個(gè)可能設(shè)定地震事件的貢獻(xiàn)比例.

    圖1 工程場(chǎng)址及所在區(qū)域主要影響潛源劃分

    Fig.1 Potential source distribution and location of target site

    圖2 不同超越概率下的一致概率譜

    Fig.2 Uniform spectrum for different exceeding probability

    表1 主要潛源震級(jí)上下限以及空間分布函數(shù)

    將各個(gè)潛源的計(jì)算結(jié)果進(jìn)行綜合疊加后,得到設(shè)定地震事件(M-R-ε)的貢獻(xiàn)比例分布見(jiàn)圖3.用不同位置帶顏色的三維柱狀圖表示設(shè)定地震組合,x軸代表震級(jí),y軸代表震中距,柱高度代表該事件對(duì)目標(biāo)值的貢獻(xiàn)比例,標(biāo)準(zhǔn)差系數(shù)用顏色表示,柱高為0表明該事件對(duì)目標(biāo)值沒(méi)有貢獻(xiàn).各工況下的主要貢獻(xiàn)地震集中在M=5.75~6.25,R=15~35 km的區(qū)域內(nèi),對(duì)應(yīng)的是36號(hào)和29號(hào)潛源,見(jiàn)圖1;30號(hào)潛源雖然距目標(biāo)場(chǎng)地較近,但由于震級(jí)上限僅為6.0,對(duì)解耦目標(biāo)并沒(méi)有顯著貢獻(xiàn),不同解耦目標(biāo)下的結(jié)果差異主要體現(xiàn)在標(biāo)準(zhǔn)差系數(shù)上.

    在求解出各可能地震事件的貢獻(xiàn)比例后,設(shè)定地震的設(shè)定思路有兩種:最大貢獻(xiàn)地震和平均影響地震.最大貢獻(xiàn)地震定義為貢獻(xiàn)比例最大的(M-R-ε)組合;平均影響地震則是將每個(gè)可能的設(shè)定地震組合按照貢獻(xiàn)比例進(jìn)行加權(quán)平均得到的結(jié)果.最大貢獻(xiàn)地震和平均影響地震均為危險(xiǎn)性分析解耦的重要產(chǎn)出.通常當(dāng)近、遠(yuǎn)場(chǎng)同時(shí)出現(xiàn)高貢獻(xiàn)比例的設(shè)定地震組合時(shí),應(yīng)當(dāng)給出不同區(qū)域?qū)?yīng)的最大貢獻(xiàn)地震作為設(shè)定地震[8].而本文的主要貢獻(xiàn)地震事件所在區(qū)域相對(duì)較為集中,為全面考慮整個(gè)潛源區(qū)內(nèi)所有微元對(duì)場(chǎng)地的影響特征,宜采用平均影響地震作為設(shè)定地震.

    不同超越概率的最大貢獻(xiàn)地震與平均影響地震結(jié)果列于表2,震級(jí)、震中距結(jié)果很接近,差異主要集中在標(biāo)準(zhǔn)差系數(shù)上,這從另一個(gè)角度說(shuō)明了僅靠震級(jí)和震中距并不足以完全描述設(shè)定地震,衰減關(guān)系離散性對(duì)結(jié)果的影響不可忽視.正如陶夏新等[16]指出的,PSHA中的目標(biāo)值已是進(jìn)行過(guò)不確定性校正的結(jié)果,如果直接借助衰減關(guān)系均值反推設(shè)定地震(M,R),得到的結(jié)果很容易與地震環(huán)境自相矛盾.

    圖3 50 a不同超越概率下設(shè)定地震(M-R-ε)的解耦結(jié)果及最大貢獻(xiàn)地震

    Tab.2 Maximum contribution and median contribution earthquake scenario forSa(T=1.0 s)

    超越概率最大貢獻(xiàn)地震平均影響地震MR/kmε貢獻(xiàn)比例/%MR/kmε50a2%6.25312.0~3.019.56.3931.822.0450a10%6.25311.0~1.515.76.3435.421.22

    2 強(qiáng)震記錄選取

    2.1 條件均值譜

    條件均值譜是以某超越概率下的目標(biāo)周期反應(yīng)譜值Sa(T*)為條件,利用地震動(dòng)衰減關(guān)系和相關(guān)系數(shù)公式計(jì)算非目標(biāo)周期處Ti的反應(yīng)譜值,從而推導(dǎo)出在整個(gè)周期段符合真實(shí)譜形規(guī)律的反應(yīng)譜,構(gòu)造公式為

    (6)

    式中:(M*-R*-ε(T*))為解耦得到的設(shè)定地震;ρ(Ti,T*)為ε(Ti)和ε(T*)的相關(guān)系數(shù),采用的是基于NGA-West1數(shù)據(jù)庫(kù)擬合得到的經(jīng)驗(yàn)公式[17],適用范圍0.01~10.0 s;μlnSa(Ti)為設(shè)定地震(M*,R*)條件下,地震動(dòng)衰減關(guān)系在周期點(diǎn)Ti處的反應(yīng)譜中值;σ(Ti)為對(duì)應(yīng)的標(biāo)準(zhǔn)差.由式(6)可知,條件均值譜介于一致概率譜和衰減關(guān)系均值之間;周期Ti與T*越遠(yuǎn),ε(Ti)和ε(T*)的相關(guān)性越低,條件反應(yīng)譜值越接近衰減關(guān)系的譜均值,否則條件反應(yīng)譜值越接近一致概率譜值.

    基于前文解耦結(jié)果構(gòu)建的條件均值目標(biāo)譜見(jiàn)圖4,在目標(biāo)周期處,條件反應(yīng)譜均值與目標(biāo)值一致,隨著周期點(diǎn)遠(yuǎn)離目標(biāo)周期點(diǎn),條件均值目標(biāo)譜逐漸偏離一致概率譜.從整個(gè)周期來(lái)看,一致概率譜整體偏于保守,這種差異在長(zhǎng)周期表現(xiàn)得尤為顯著.

    圖4 一致概率譜,條件均值目標(biāo)譜與衰減關(guān)系均值譜對(duì)比

    2.2 強(qiáng)震記錄選取

    分別選取50 a超越概率2%下Sa(T=1.0 s)的條件均值譜(CMS)、一致概率譜(UHS)作為目標(biāo)譜進(jìn)行強(qiáng)震記錄選取.強(qiáng)震記錄數(shù)據(jù)庫(kù)采用的是目前廣泛用于工程科研實(shí)踐且地震信息全面可靠的美國(guó)NGA-West1強(qiáng)震記錄數(shù)據(jù)庫(kù).首先結(jié)合設(shè)定地震信息對(duì)震級(jí)、震中距以及場(chǎng)地條件等重要地震信息參數(shù)加以初選.震級(jí)范圍為設(shè)定地震震級(jí)±0.5,即M5.7~M6.7;考慮到實(shí)際潛源分布和場(chǎng)地的相對(duì)位置,以震中距10~100 km為記錄篩選范圍;由于在PSHA計(jì)算過(guò)程中包括后面的解耦過(guò)程中所用的衰減關(guān)系均沒(méi)有考慮場(chǎng)地的土層放大,只適用于基巖場(chǎng)地,故規(guī)定待選記錄的臺(tái)站場(chǎng)地VS30(30 m覆蓋土層平均剪切波速)不得低于550 m/s.按上述要求進(jìn)行震級(jí)、震中距、場(chǎng)地條件篩選后,對(duì)滿足條件的加速度記錄以Sa(T=1.0 s)作為目標(biāo)值進(jìn)行線性放縮,調(diào)幅系數(shù)范圍為0.5~2.0倍.最后計(jì)算所有滿足條件記錄與目標(biāo)譜在全周期段的最小均方差,選取全部數(shù)據(jù)庫(kù)中均方差值最小的前7條記錄作為最終選取結(jié)果.這兩種不同目標(biāo)譜下的選波結(jié)果見(jiàn)圖5,用選波結(jié)果平均反應(yīng)譜與目標(biāo)譜的相對(duì)誤差評(píng)價(jià)匹配結(jié)果在各個(gè)周期點(diǎn)的準(zhǔn)確性,用選波結(jié)果在各個(gè)周期點(diǎn)對(duì)數(shù)值的標(biāo)準(zhǔn)差評(píng)價(jià)離散性.

    圖5 不同目標(biāo)譜下的選波結(jié)果

    由圖5可見(jiàn),基于CMS目標(biāo)譜的選取結(jié)果在絕大多數(shù)周期段相對(duì)誤差保持20%以?xún)?nèi),僅在5.0 s之后的相對(duì)誤差超過(guò)了20%,最大相對(duì)誤差低于30%,對(duì)數(shù)標(biāo)準(zhǔn)差在0.4以?xún)?nèi),均值和離散性的匹配結(jié)果均令人滿意.而基于UHS目標(biāo)譜的選取結(jié)果,在周期高于3.0 s后選取結(jié)果平均反應(yīng)譜已明顯低于目標(biāo)譜,最大誤差達(dá)到50%,而且長(zhǎng)周期處的對(duì)數(shù)標(biāo)準(zhǔn)差達(dá)到了0.6.從選波結(jié)果的地震信息來(lái)看,基于UHS目標(biāo)譜所選記錄的震中距整體偏近,7條記錄中一半以上的震中距低于15 km,與設(shè)定地震信息有明顯差距,也不滿足實(shí)際的潛源分布位置.而基于CMS目標(biāo)譜的選取結(jié)果地震信息整體上與目標(biāo)值吻合較好,震級(jí)與震中距均值與設(shè)定地震目標(biāo)很接近.

    3 地震動(dòng)選取建議流程

    將前文基于中國(guó)地震安全性評(píng)價(jià)PSHA分析結(jié)果進(jìn)行強(qiáng)震記錄選取的整體流程概括整理為圖6:首先,圍繞待建結(jié)構(gòu)的相關(guān)參數(shù)對(duì)PSHA結(jié)果進(jìn)行解耦,整理計(jì)算得到對(duì)應(yīng)的設(shè)定地震信息(M-R-ε);然后,利用設(shè)定地震信息對(duì)備選強(qiáng)震記錄進(jìn)行地震動(dòng)特性方面的篩選后,構(gòu)建條件均值目標(biāo)譜進(jìn)行強(qiáng)震記錄匹配.這一選取流程在解耦、篩選記錄、目標(biāo)譜構(gòu)建的整個(gè)過(guò)程均兼顧了結(jié)構(gòu)特性和地震動(dòng)離散性,選取結(jié)果作為基巖地震動(dòng)輸入進(jìn)行土層反應(yīng)計(jì)算后可用于待建結(jié)構(gòu)的時(shí)程分析或設(shè)計(jì)驗(yàn)算中,避免了傳統(tǒng)流程中構(gòu)建規(guī)范譜和重復(fù)匹配引入的誤差.

    圖6 中國(guó)地震安全性評(píng)價(jià)的天然地震動(dòng)選取建議流程

    Fig.6 Suggested procedure of recordings selection for Chinese seismic safety evaluation work

    4 結(jié) 論

    1)基于中國(guó)地震安全性評(píng)價(jià)PSHA的結(jié)果實(shí)現(xiàn)了同時(shí)考慮震級(jí)、震中距、標(biāo)準(zhǔn)差系數(shù)的地震信息解耦,將各主要潛源離散化得到各個(gè)可能的設(shè)定地震(M-R-ε)組合對(duì)目標(biāo)值的貢獻(xiàn),得到最大貢獻(xiàn)地震和平均影響地震,為構(gòu)建目標(biāo)譜和選取強(qiáng)震記錄提供考慮結(jié)構(gòu)特性與離散性的設(shè)定地震信息.

    2)在中國(guó)安全性評(píng)價(jià)解耦結(jié)果的基礎(chǔ)上,構(gòu)建了目前在國(guó)外已得到廣泛認(rèn)可和應(yīng)用的條件均值譜.依據(jù)該目標(biāo)譜選取得到的強(qiáng)震記錄結(jié)果在地震信息、均值和離散性上均較理想.

    3)在中國(guó)地震安全性評(píng)價(jià)現(xiàn)有工作基礎(chǔ)上,作為對(duì)人造地震動(dòng)時(shí)程的補(bǔ)充,提出了一套具有操作性的天然強(qiáng)震記錄選取流程,可與中國(guó)現(xiàn)有安全性評(píng)價(jià)工作較好銜接.

    [1] 李小軍. 工程場(chǎng)地地震安全性評(píng)價(jià)工作及相關(guān)技術(shù)問(wèn)題[J]. 震災(zāi)防御技術(shù), 2006, 1(1): 15-24. LI Xiaojun. Evaluation of seismic safety for engineering sites and relevant technical problems[J]. Technology for Earthquake Disaster Prevention, 2006, 1(1):15-24.

    [2] 李琳. 強(qiáng)震動(dòng)記錄選取中目標(biāo)譜的研究及應(yīng)用[D]. 哈爾濱: 中國(guó)地震局工程力學(xué)研究所, 2015. LI Ling. Research and application of target spectrum in the selection of strong motion records[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration, 2015.

    [3] BOMMER J J, SCOTT S G, SARMA S K. Hazard-consistent earthquakescenarios [J]. Soil Dynamics and Earthquake Engineering, 2000, 31:611-625.

    [4] NAEIM F, LEW M. On the use of design spectrum compatible timehistories [J]. Earthquake Spectra, 1995, 11(1):111-127.

    [5] BAKER J W. Conditional mean spectrum: tool for ground motion selection [J]. Journal of Structural Engineering, 2011, 137(3):322-331.

    [6] NEHRP(National Earthquake Hazards Reduction Program). Selecting and scaling earthquake ground motions for performing response-history analysis [R].

    [S.l.]:National Institute of Standards and Technology, 2011.

    [7] 陳波. 結(jié)構(gòu)非線性動(dòng)力分析中地震動(dòng)記錄的選擇和調(diào)整方法研究[D]. 北京:中國(guó)地震局地球物理研究所, 2013. CHEN Bo. Ground motion selection and modification methods for performing nonlinear dynamic analysis of buildings[D]. Beijing: Institute of Geophysics, China Earthquake Administration, 2013.

    [8] MCGUIRE R K. Probabilistic seismic hazard analysis and design earthquakes: closing the loop [J]. Bulletin of the Seismological Society ofAmerica, 1995, 85(5):1275-1284.

    [9] ASCE/SEI7 05. Minimum design loads for buildings and other structures[S].Reston: American Society of Civil, 2005.

    [10]Eurocode 8. European Standard [S]. Brussels: European Committee for Standardization, 2004.

    [11]NZS. 1170.5. Structural design actions[S]. Wellington: Standards New Zealand, 2004.

    [12]BAKER J W, CORNELL C A. A vector-valued ground motion intensity measure consisting of spectral acceleration and epsilon [J]. Earthquake Engineering and Structural Dynamics, 2005(34):1193-1217.

    [13]BAKER J W, CORNELL C A. Spectral shape, epsilon and record selection [J]. Earthquake Engineering and Structural Dynamics, 2006(35): 1077-1095.

    [14]高孟潭. 潛在震源區(qū)期望震級(jí)和期望距離及其計(jì)算方法[J].地震學(xué)報(bào), 1994, 16(3): 346-351. GAO Mengtan. Expected magnitude and distance calculation method for potential seismic source region[J]. Earthquake Science, 1994, 16(3):346-351.

    [15]霍俊榮. 近場(chǎng)強(qiáng)地面運(yùn)動(dòng)衰減規(guī)律的研究[D].哈爾濱:中國(guó)地震局工程力學(xué)研究所, 1989. HUO Junrong. Study on the near-field strong ground motion attenuation relationship[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration, 1989.

    [16]陶夏新,陶正如,師黎靜.設(shè)定地震——概率地震危險(xiǎn)性評(píng)估和確定性危險(xiǎn)性評(píng)估[J].地震工程與工程振動(dòng),2014, 34(4): 101-109. TAO Xiaxin, TAO Zhenru, SHI Lijing. Scenario earthquake-link of probabilistic seismic hazard assessment and deterministic seismic hazard assessment[J]. Journal of Earthquake Engineering and Engineering Vibration, 2014, 34(4): 101-109.

    [17]BAKER J W, JAYARAM N. Correlation of spectral acceleration values from NGA ground motion models [J]. Earthquake Spectra, 2008, 24(1):299-317.

    (編輯 趙麗瑩)

    Real strong ground-motion selection for Chinese seismic safety evaluation work

    JI Kun, WEN Ruizhi, REN Yefei

    (Key Laboratory of Earthquake Engineering and Engineering Vibration, Institute of Engineering Mechanics,China Earthquake Administration, Harbin 150080,China)

    To select real strong motions for Chinese seismic safety evaluation work, we disaggregated the result of Chinese seismic safety evaluation work into specified earthquake scenarios considering magnitude, distance and earthquake variation, and then the condition median spectrum was constructed as strong motion selection target. The results based on conditional mean spectrum was compared with uniform spectrum matching result for a real safety evaluation project. The comparison result indicates that the result from the proposed method is consistent with target earthquake scenario, and the proposed method preforms better in matching spectrum shape within the whole period. Finally, a practical real ground motion selection method was proposed for Chinese seismic safety evaluation work.

    seismic safety evaluation; real strong ground motions; earthquake scenario; conditional mean spectrum; strong ground motion selection

    10.11918/j.issn.0367-6234.2016.12.026

    2015-08-25

    地震行業(yè)科研專(zhuān)項(xiàng)(201508005); 黑龍江省科學(xué)基金(LC2015022);國(guó)家自然科學(xué)基金(51308515)

    冀 昆(1990—),男,博士研究生; 溫瑞智(1968—),男,研究員,博士生導(dǎo)師

    溫瑞智,ruizhi@iem.ac.cn

    P315.9

    A

    0367-6234(2016)12-0183-06

    国产免费福利视频在线观看| 能在线免费看毛片的网站| 丝袜脚勾引网站| 观看免费一级毛片| 久久女婷五月综合色啪小说 | 99久久精品国产国产毛片| 一级av片app| 久久精品熟女亚洲av麻豆精品| 国产熟女欧美一区二区| 免费av毛片视频| 视频中文字幕在线观看| 啦啦啦啦在线视频资源| 亚洲最大成人中文| 亚洲成人精品中文字幕电影| 黄色怎么调成土黄色| 乱系列少妇在线播放| 在线观看一区二区三区激情| 亚洲av日韩在线播放| 免费不卡的大黄色大毛片视频在线观看| 国产成人a∨麻豆精品| 男人狂女人下面高潮的视频| 黄色配什么色好看| 亚洲国产av新网站| 国产在视频线精品| 天堂网av新在线| 日本免费在线观看一区| 久久人人爽人人爽人人片va| 国产成人aa在线观看| 男人添女人高潮全过程视频| 亚洲精品视频女| 日韩av免费高清视频| 美女cb高潮喷水在线观看| 久久精品国产a三级三级三级| 精品酒店卫生间| av免费观看日本| 美女xxoo啪啪120秒动态图| 熟女人妻精品中文字幕| 亚洲图色成人| 午夜视频国产福利| 成人毛片60女人毛片免费| 男女边摸边吃奶| 亚洲国产欧美人成| 久久精品夜色国产| 亚洲性久久影院| 国产精品偷伦视频观看了| 狂野欧美激情性xxxx在线观看| 亚洲成人久久爱视频| 免费人成在线观看视频色| 国产成人精品久久久久久| 免费不卡的大黄色大毛片视频在线观看| 一级毛片我不卡| 国产乱人偷精品视频| av免费在线看不卡| 男人和女人高潮做爰伦理| 久久久a久久爽久久v久久| 欧美日韩综合久久久久久| 久热这里只有精品99| 成人亚洲欧美一区二区av| 国产精品秋霞免费鲁丝片| 亚洲aⅴ乱码一区二区在线播放| 晚上一个人看的免费电影| 久久鲁丝午夜福利片| 国产 一区精品| 深爱激情五月婷婷| 免费黄网站久久成人精品| 白带黄色成豆腐渣| 成人国产麻豆网| 欧美成人a在线观看| 亚洲图色成人| 国产免费一级a男人的天堂| 午夜激情久久久久久久| 国产淫片久久久久久久久| 日韩视频在线欧美| 永久免费av网站大全| 麻豆国产97在线/欧美| 黄色视频在线播放观看不卡| 欧美激情国产日韩精品一区| 美女国产视频在线观看| 国产 精品1| 另类亚洲欧美激情| av国产免费在线观看| 国产高清有码在线观看视频| 观看免费一级毛片| freevideosex欧美| 卡戴珊不雅视频在线播放| 熟妇人妻不卡中文字幕| 黄色日韩在线| 少妇的逼好多水| 狂野欧美激情性bbbbbb| 黄色配什么色好看| 精品午夜福利在线看| 成年人午夜在线观看视频| 成人欧美大片| 亚洲精品日韩av片在线观看| 丰满少妇做爰视频| 美女被艹到高潮喷水动态| 天天躁夜夜躁狠狠久久av| 日韩强制内射视频| 熟妇人妻不卡中文字幕| 国内精品美女久久久久久| 亚洲内射少妇av| 狂野欧美激情性bbbbbb| 少妇高潮的动态图| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产伦在线观看视频一区| 精品人妻视频免费看| 国产高清有码在线观看视频| 国产一区有黄有色的免费视频| 永久免费av网站大全| 国产亚洲最大av| 欧美丝袜亚洲另类| .国产精品久久| 亚洲av免费在线观看| 亚洲成人久久爱视频| 尾随美女入室| 中国三级夫妇交换| 久久久久性生活片| 久久女婷五月综合色啪小说 | 国产91av在线免费观看| 两个人的视频大全免费| 久久ye,这里只有精品| 婷婷色av中文字幕| 国产老妇女一区| 亚洲自偷自拍三级| 国产黄色视频一区二区在线观看| 国产毛片a区久久久久| 99热网站在线观看| 狂野欧美白嫩少妇大欣赏| 国产成人一区二区在线| 日日撸夜夜添| 久久人人爽人人片av| 99re6热这里在线精品视频| 日韩不卡一区二区三区视频在线| 亚洲欧美中文字幕日韩二区| 亚洲国产日韩一区二区| av又黄又爽大尺度在线免费看| 狂野欧美激情性bbbbbb| 人人妻人人澡人人爽人人夜夜| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 在线天堂最新版资源| 国产成人精品婷婷| 国产一级毛片在线| 三级国产精品欧美在线观看| 91在线精品国自产拍蜜月| 97在线人人人人妻| 亚洲一区二区三区欧美精品 | 国产av码专区亚洲av| 又粗又硬又长又爽又黄的视频| 免费黄色在线免费观看| 亚洲成色77777| 国产男人的电影天堂91| 久热久热在线精品观看| 人妻夜夜爽99麻豆av| 我要看日韩黄色一级片| 97在线人人人人妻| 日韩一本色道免费dvd| 色视频在线一区二区三区| 免费黄网站久久成人精品| 欧美xxⅹ黑人| 啦啦啦中文免费视频观看日本| videos熟女内射| 精品人妻偷拍中文字幕| 天天躁日日操中文字幕| 欧美变态另类bdsm刘玥| 3wmmmm亚洲av在线观看| 国产老妇伦熟女老妇高清| 中文资源天堂在线| 亚洲精品色激情综合| 在线免费观看不下载黄p国产| 一个人看视频在线观看www免费| 日韩制服骚丝袜av| 麻豆乱淫一区二区| 国产有黄有色有爽视频| 天美传媒精品一区二区| 男人添女人高潮全过程视频| 免费播放大片免费观看视频在线观看| 狂野欧美白嫩少妇大欣赏| 黄色一级大片看看| 中文精品一卡2卡3卡4更新| 日韩一区二区三区影片| 国产免费福利视频在线观看| 亚洲国产精品成人综合色| 亚洲人成网站高清观看| 日韩中字成人| 99热6这里只有精品| 精品酒店卫生间| 一边亲一边摸免费视频| 国产乱人偷精品视频| 日韩av不卡免费在线播放| 国产av国产精品国产| 99热这里只有是精品50| 看非洲黑人一级黄片| 久久久久网色| 国产高清不卡午夜福利| 色视频在线一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆成人av视频| 亚洲综合色惰| 嘟嘟电影网在线观看| 午夜福利在线观看免费完整高清在| 舔av片在线| 精品人妻一区二区三区麻豆| 欧美日韩视频精品一区| 丝袜脚勾引网站| 99久久九九国产精品国产免费| 国产成年人精品一区二区| 热re99久久精品国产66热6| 欧美区成人在线视频| 激情五月婷婷亚洲| 老女人水多毛片| 婷婷色麻豆天堂久久| 一区二区三区免费毛片| 小蜜桃在线观看免费完整版高清| 草草在线视频免费看| 人妻少妇偷人精品九色| 一级毛片黄色毛片免费观看视频| 一区二区三区四区激情视频| 亚洲丝袜综合中文字幕| 最近最新中文字幕免费大全7| 免费在线观看成人毛片| 亚洲欧美一区二区三区黑人 | 搡老乐熟女国产| 亚洲国产精品国产精品| 亚洲熟女精品中文字幕| 五月天丁香电影| 亚洲欧美日韩东京热| 在线观看三级黄色| 中文字幕亚洲精品专区| 国产极品天堂在线| av女优亚洲男人天堂| 免费大片黄手机在线观看| 国产成人精品婷婷| 国产大屁股一区二区在线视频| 亚洲欧美精品专区久久| 黄色视频在线播放观看不卡| 免费观看在线日韩| 国产av国产精品国产| 日韩欧美一区视频在线观看 | 亚洲精品日韩av片在线观看| 久久精品国产亚洲网站| 两个人的视频大全免费| 亚洲欧美清纯卡通| 亚洲av二区三区四区| 国产视频首页在线观看| 久久精品熟女亚洲av麻豆精品| 国内揄拍国产精品人妻在线| 又爽又黄无遮挡网站| 永久免费av网站大全| 亚洲精品影视一区二区三区av| 人人妻人人看人人澡| 日韩大片免费观看网站| 久久精品国产亚洲av天美| 亚洲精品视频女| 日本一二三区视频观看| 国产精品成人在线| 99久久精品国产国产毛片| 人妻制服诱惑在线中文字幕| 深爱激情五月婷婷| 国产精品人妻久久久久久| 丰满人妻一区二区三区视频av| 国产v大片淫在线免费观看| 国内少妇人妻偷人精品xxx网站| 男男h啪啪无遮挡| 尾随美女入室| 日韩强制内射视频| 精品视频人人做人人爽| 亚洲自偷自拍三级| 久久精品国产亚洲av天美| av播播在线观看一区| 在现免费观看毛片| 国产有黄有色有爽视频| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品999| 又粗又硬又长又爽又黄的视频| 亚洲精品乱码久久久v下载方式| 国产黄频视频在线观看| 国产淫片久久久久久久久| 欧美最新免费一区二区三区| 91久久精品国产一区二区成人| 欧美日本视频| 丰满乱子伦码专区| tube8黄色片| 婷婷色综合大香蕉| 久久久久国产精品人妻一区二区| 交换朋友夫妻互换小说| 一本久久精品| 插逼视频在线观看| 18禁在线无遮挡免费观看视频| 搞女人的毛片| 欧美一区二区亚洲| 欧美高清成人免费视频www| 五月玫瑰六月丁香| 男女那种视频在线观看| 最近最新中文字幕免费大全7| 欧美极品一区二区三区四区| 亚洲精品成人久久久久久| 最近2019中文字幕mv第一页| 国产免费福利视频在线观看| 毛片女人毛片| 欧美xxⅹ黑人| 国产精品熟女久久久久浪| 白带黄色成豆腐渣| 在线观看免费高清a一片| 中文资源天堂在线| 成人漫画全彩无遮挡| 秋霞伦理黄片| 中文精品一卡2卡3卡4更新| 亚洲天堂国产精品一区在线| 亚洲av.av天堂| 狠狠精品人妻久久久久久综合| 街头女战士在线观看网站| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| 国产美女午夜福利| 亚洲av欧美aⅴ国产| 99久国产av精品国产电影| 亚洲av二区三区四区| av国产精品久久久久影院| 亚洲性久久影院| 国产视频首页在线观看| 18+在线观看网站| 国产成人福利小说| 久热这里只有精品99| 啦啦啦在线观看免费高清www| 国产乱人偷精品视频| 26uuu在线亚洲综合色| 啦啦啦啦在线视频资源| 欧美日韩视频精品一区| 成人二区视频| 亚洲精品成人久久久久久| 热re99久久精品国产66热6| 日日啪夜夜撸| 成年免费大片在线观看| 街头女战士在线观看网站| 久久精品国产亚洲网站| 只有这里有精品99| 女人被狂操c到高潮| 亚洲一区二区三区欧美精品 | 精品久久久久久久久亚洲| 久久久久性生活片| 久久99蜜桃精品久久| 麻豆乱淫一区二区| 黄片无遮挡物在线观看| 干丝袜人妻中文字幕| 蜜臀久久99精品久久宅男| 亚洲精品成人av观看孕妇| 91久久精品电影网| 国产精品一区www在线观看| 麻豆久久精品国产亚洲av| 99久久精品国产国产毛片| 国产有黄有色有爽视频| 一本色道久久久久久精品综合| 国产v大片淫在线免费观看| 全区人妻精品视频| 亚洲国产精品成人综合色| 亚洲国产日韩一区二区| 夫妻午夜视频| 国产在线一区二区三区精| 最新中文字幕久久久久| 国产久久久一区二区三区| 美女主播在线视频| 久久久久久久国产电影| 成人国产麻豆网| 免费黄网站久久成人精品| 男女下面进入的视频免费午夜| 亚洲国产成人一精品久久久| 国产精品秋霞免费鲁丝片| 亚洲精品乱码久久久v下载方式| 久久久久精品久久久久真实原创| 九九爱精品视频在线观看| 精品国产乱码久久久久久小说| 欧美成人精品欧美一级黄| 少妇人妻一区二区三区视频| 日本黄色片子视频| 亚洲熟女精品中文字幕| 熟妇人妻不卡中文字幕| 亚洲av.av天堂| 国产亚洲av片在线观看秒播厂| 高清日韩中文字幕在线| 国产日韩欧美在线精品| 精品久久久久久久久av| 亚洲欧美成人精品一区二区| 欧美一区二区亚洲| 大香蕉久久网| 久久午夜福利片| 成人亚洲欧美一区二区av| 日韩欧美精品免费久久| 永久网站在线| 日韩 亚洲 欧美在线| 亚洲一区二区三区欧美精品 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品乱码久久久v下载方式| 久久99热这里只有精品18| 国产久久久一区二区三区| 亚洲国产精品专区欧美| 日本一本二区三区精品| 有码 亚洲区| 少妇的逼水好多| 亚洲精品国产色婷婷电影| 熟女av电影| 欧美日韩视频高清一区二区三区二| 久久精品国产亚洲av涩爱| 国产av国产精品国产| 一个人看的www免费观看视频| 97热精品久久久久久| 人体艺术视频欧美日本| 老司机影院成人| 欧美精品一区二区大全| 一级黄片播放器| 中文欧美无线码| 欧美日韩亚洲高清精品| 七月丁香在线播放| 中文字幕久久专区| 日韩强制内射视频| 看免费成人av毛片| 国语对白做爰xxxⅹ性视频网站| 肉色欧美久久久久久久蜜桃 | 少妇人妻一区二区三区视频| 高清毛片免费看| 亚洲国产精品专区欧美| 日韩电影二区| 一级二级三级毛片免费看| 欧美+日韩+精品| 丝袜美腿在线中文| 夜夜看夜夜爽夜夜摸| 日韩亚洲欧美综合| 免费高清在线观看视频在线观看| 国产v大片淫在线免费观看| 日韩不卡一区二区三区视频在线| 国产精品麻豆人妻色哟哟久久| 亚洲av福利一区| 欧美变态另类bdsm刘玥| 免费观看性生交大片5| 国产白丝娇喘喷水9色精品| 最后的刺客免费高清国语| 国产精品一区二区性色av| 少妇人妻精品综合一区二区| 婷婷色综合大香蕉| 久久久a久久爽久久v久久| 日本猛色少妇xxxxx猛交久久| 日韩在线高清观看一区二区三区| 免费观看无遮挡的男女| 一级毛片黄色毛片免费观看视频| 亚洲欧美一区二区三区国产| 哪个播放器可以免费观看大片| 免费黄网站久久成人精品| 丰满人妻一区二区三区视频av| 成人一区二区视频在线观看| 美女xxoo啪啪120秒动态图| 国产成人精品福利久久| 精品久久久精品久久久| 国产永久视频网站| 婷婷色av中文字幕| 国产成人免费观看mmmm| a级一级毛片免费在线观看| 一级毛片久久久久久久久女| 国产高清不卡午夜福利| 国产精品一区www在线观看| 亚洲自拍偷在线| 国产淫片久久久久久久久| 国产亚洲5aaaaa淫片| 欧美极品一区二区三区四区| 亚洲精品色激情综合| 赤兔流量卡办理| 亚洲精品久久午夜乱码| 插逼视频在线观看| 免费黄色在线免费观看| 欧美bdsm另类| 亚洲在线观看片| 国产视频首页在线观看| 亚洲色图av天堂| 禁无遮挡网站| 免费高清在线观看视频在线观看| 九九爱精品视频在线观看| 777米奇影视久久| 欧美老熟妇乱子伦牲交| 好男人视频免费观看在线| 国产女主播在线喷水免费视频网站| 久久久久性生活片| 亚洲国产av新网站| 精品视频人人做人人爽| 校园人妻丝袜中文字幕| 99视频精品全部免费 在线| 国产成人aa在线观看| 精华霜和精华液先用哪个| 观看免费一级毛片| 听说在线观看完整版免费高清| 国产欧美亚洲国产| 菩萨蛮人人尽说江南好唐韦庄| 大又大粗又爽又黄少妇毛片口| av在线天堂中文字幕| 女人久久www免费人成看片| 亚洲av免费高清在线观看| 亚洲一区二区三区欧美精品 | 啦啦啦在线观看免费高清www| 欧美日本视频| 丰满少妇做爰视频| freevideosex欧美| 久久久久久久久久久丰满| 亚洲欧洲日产国产| 国产成人精品福利久久| 亚洲精品日本国产第一区| 国产一区二区三区av在线| 中文在线观看免费www的网站| av一本久久久久| 国产精品一区二区三区四区免费观看| 热99国产精品久久久久久7| 日日啪夜夜撸| 97在线人人人人妻| 久久久精品欧美日韩精品| 国产白丝娇喘喷水9色精品| 在线免费十八禁| av国产免费在线观看| 国产日韩欧美在线精品| 午夜免费鲁丝| av.在线天堂| a级毛色黄片| 日日撸夜夜添| 尾随美女入室| 美女cb高潮喷水在线观看| 只有这里有精品99| 最近最新中文字幕大全电影3| 成人特级av手机在线观看| 五月开心婷婷网| 成年人午夜在线观看视频| 亚洲熟女精品中文字幕| 国产成人a区在线观看| 国产精品女同一区二区软件| 色网站视频免费| 国产黄色视频一区二区在线观看| 搡女人真爽免费视频火全软件| 精品一区二区三卡| 久久6这里有精品| 边亲边吃奶的免费视频| 中文字幕免费在线视频6| 有码 亚洲区| 成人免费观看视频高清| 亚洲欧美一区二区三区国产| 日韩制服骚丝袜av| 91精品一卡2卡3卡4卡| 好男人视频免费观看在线| 亚洲精品日韩在线中文字幕| 97在线人人人人妻| av在线天堂中文字幕| 一级毛片黄色毛片免费观看视频| 热re99久久精品国产66热6| 爱豆传媒免费全集在线观看| 欧美人与善性xxx| 99久国产av精品国产电影| 少妇人妻一区二区三区视频| 精品一区二区三卡| 色播亚洲综合网| 欧美精品人与动牲交sv欧美| 国产91av在线免费观看| a级毛片免费高清观看在线播放| 精品酒店卫生间| 亚洲人与动物交配视频| 亚洲综合色惰| 欧美+日韩+精品| 日韩亚洲欧美综合| 久久精品国产自在天天线| 午夜免费观看性视频| 男男h啪啪无遮挡| 国产黄色免费在线视频| 免费黄频网站在线观看国产| 日本av手机在线免费观看| 纵有疾风起免费观看全集完整版| 九九爱精品视频在线观看| 亚洲欧美日韩无卡精品| 99热这里只有是精品50| 国精品久久久久久国模美| 午夜福利高清视频| 高清在线视频一区二区三区| 青青草视频在线视频观看| 国产免费福利视频在线观看| 中文字幕免费在线视频6| av一本久久久久| 高清午夜精品一区二区三区| 免费看日本二区| 日本猛色少妇xxxxx猛交久久| .国产精品久久| 国产亚洲一区二区精品| 国产男人的电影天堂91| 看十八女毛片水多多多| 三级国产精品片| 精品午夜福利在线看| 色综合色国产| 人妻系列 视频| 大码成人一级视频| 久久久久久久久久久免费av| 国产乱来视频区| 久久久久久伊人网av| 亚洲在久久综合| 水蜜桃什么品种好| 七月丁香在线播放| 日韩,欧美,国产一区二区三区| 日韩 亚洲 欧美在线| 国产成人免费观看mmmm| 91精品一卡2卡3卡4卡| 亚洲欧美清纯卡通| eeuss影院久久| 亚洲av中文字字幕乱码综合| 插阴视频在线观看视频| 亚洲成人中文字幕在线播放| 小蜜桃在线观看免费完整版高清| 美女cb高潮喷水在线观看| 日韩伦理黄色片| 久久亚洲国产成人精品v| 久久久久网色| 少妇人妻久久综合中文| 不卡视频在线观看欧美| 亚洲成人久久爱视频| 免费不卡的大黄色大毛片视频在线观看| 国产成人精品福利久久| 国产视频内射|