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

    利用連續(xù)重力觀測約束2014智利Iquique地震的震源機制解

    2016-07-29 02:04:43江穎李輝劉子維申重陽張曉彤
    地球物理學報 2016年3期
    關鍵詞:標量臺站震源

    江穎,李輝,劉子維,申重陽,張曉彤

    1 中國地震局地震研究所(中國地震局地震大地測量重點實驗室),武漢 430071 2 中國科學院測量與地球物理研究所大地測量與地球動力學國家重點實驗室,武漢 430077

    ?

    利用連續(xù)重力觀測約束2014智利Iquique地震的震源機制解

    江穎1,2,李輝1,劉子維1,申重陽1,張曉彤1

    1 中國地震局地震研究所(中國地震局地震大地測量重點實驗室),武漢4300712中國科學院測量與地球物理研究所大地測量與地球動力學國家重點實驗室,武漢430077

    摘要2014年4月1日,智利北部Iquique地區(qū)近海發(fā)生MW8.1地震,地震發(fā)生之后,國際上一些著名的地震科研機構和學者采用不同的數(shù)據(jù)和方法計算得到此次地震的震源機制解,但這些結果存在較大差異.地球長周期自由振蕩的振幅主要依賴于地震矩的大小及地震斷層的破裂方式,可以很好地約束地震震源機制.因此,本文根據(jù)2014年Iquique地震現(xiàn)有的6個不同震源機制解模擬計算了該地震激發(fā)的自由振蕩信號,并與全國連續(xù)重力臺網中16個彈簧重力儀的觀測結果進行比對,基于1.5~5.3 mHz的球型簡正模分析和約束了Iquique地震的震源機制解.研究發(fā)現(xiàn),基于美國地質調查局WPhase Moment Tensor Solution反演的震源機制解的自由振蕩模擬值與實際觀測符合最好,其相應的震級能較好反映Iquique地震釋放的總能量,而利用海嘯數(shù)據(jù)反演的標量地震矩偏小,聯(lián)合遠場和近場長周期觀測數(shù)據(jù)反演可顯著改善震源機制解.另外,還基于格爾木重力臺站的模擬值與觀測值定量分析了不同震源機制解參數(shù)對自由振蕩振幅的影響.結果表明地震的標量地震矩M0對自由振蕩振幅的影響最大,而斷層走向、傾角、滑動方向角和震源深度對自由振蕩的振幅影響相對較小.

    關鍵詞Iquique地震; 震源機制解; 地球自由振蕩; 重力觀測

    1引言

    智利西部沿海處于納斯卡板塊和南美板塊之間,由于納斯卡板塊以約8 cm·a-1的收斂速度向南美板塊俯沖(Angermann et al.,1999),該地區(qū)大地震頻發(fā).2014年4月1日23∶46 (UTC)在智利西北部Iquique近海地區(qū)發(fā)生了震級為MW8.1的地震,震中位為19.642°S,70.817°W(圖1,紅色五角星位置),震源深度25 km (USGS National Earthquake Information Center-NEIC: http:∥earthquake.usgs.gov/regional/neic/).地震發(fā)生后,Iquique地震的震源性質和發(fā)震構造吸引了許多學者的關注和研究(An et al.,2014;Lay et al.,2014;Ruiz et al.,2014;Schurr et al.,2014;Kato et al.,2014;Yagi et al.,2014;Bürgmann et al.,2014;Hayes et al.,2014;Liu et al.,2015).國際上一些著名的地震科研機構和學者采用不同的數(shù)據(jù)和方法計算得到Iquique地震的震源機制解,見表1.其中震源機制解1是GCMT(The Global Centroid-Moment-Tensor Project)通過計算單位矩張量元的格林函數(shù),將格林函數(shù)與震源時間函數(shù)做卷積,運用最小二乘求得最優(yōu)矩張量解.GCMT采用體波(截斷周期為50 s),長周期面波(截斷周期為50 s)和地幔波(截斷周期為200 s)計算了Iquique地震的震源機制解.震源機制解2、3是美國地質調查局(USGS)的結果,USGSCMT(USGS Centroid Moment Tensor Solution)采用體波和長周期面波反演震源機制,USGSWPHASE(USGS WPhase Moment Tensor Solution)采用W-phase(100~1000 s)的波形資料反演震源機制.對于中級地震,USGSBODY(USGS Body-Wave Moment Tensor Solution)采用遠場臺站記錄的體波數(shù)據(jù)反演震源機制,比如2013年的蘆山地震.單純用體波數(shù)據(jù)反演的標量地震矩偏小(Jiang et al.,2014a).通常情況下USGSBODY不提供大地震的震源機制解結果,如2011年日本Tohoku MW9.1地震等,包括此次Iquique MW8.1地震.Yagi 等(2014)利用48個臺站的遠震P波數(shù)據(jù),考慮了格林函數(shù)的不確定性反演了Iquique地震的震源機制解.Lay等(2014)利用遠程P波、SH波及海嘯數(shù)據(jù)進行有限斷層反演,得到的標量地震矩為1.7×1021N·m.An等(2014)利用3個深海評估和海嘯臺站的海嘯數(shù) 據(jù),反演得到了Iquique地震的滑動分布,并得到該地震的標量地震矩為1.22×1021N·m.從上述研究結果來看,不同的研究機構和學者給出的Iquique地震的震源機制解相差較大,特別是標量地震矩M0之間的差異顯著,最大值約為最小值的2倍,而這種差異不容忽視.

    表1 智利Iquique地震震源機制解一覽表

    注:文獻來源震源機制解1為Global-CMT Project的結果;震源機制解2和3分別為USGS提供2個不同的結果:USGSCMT、USGSWPHASE;4為Yagi et al.(2014)的結果;5為Lay et al.(2014)的結果;6為An et al.(2014)的結果.

    地震在瞬時釋放巨大能量、產生大面積地表破裂的同時,激發(fā)大量的地震波,并且導致地球整體的震動.地球的自由振蕩一旦被激發(fā),將以駐波形式穿過地球的內部介質.由于自由振蕩能量的大小與震源的破裂方式和破裂程度密切相關(Geller and Stein,1979;Stein et al.,2005),因此地震后的自由振蕩信號可用于約束地震震級和檢驗地震的震源機制解,它也是目前可對震源機制解進行總體檢驗和約束的非常有效的方法.利用自由振蕩信號約束MW9.0級以上大地震的震源機制解已經得到了不錯的效果,比如2004-12-26蘇門答臘大地震后,利用應變儀、長周期地震儀及超導重力儀記錄到的球型低頻模態(tài)0S0、0S2、0S3、0S4等,徑向模0S0、1S0等與理論模擬結果對比,發(fā)現(xiàn)相應的標量地震矩是GCMT結果的2.5~2.67倍,并且得出斷層的破裂長度、持續(xù)時間及破裂速度(Stein et al.,2005;Park et al.,2005;Lambotte et al.,2006).2011-03-11Tohoku大地震后,薛秀秀等(2012)利用0S0振型對Tohoku大地震的多個震源機制解進行分析和約束.對于中級地震來講低頻態(tài)的信號不易提取,Jiang等(2014)根據(jù)2013-04-20蘆山地震的4種不同震源機制解計算了蘆山地震頻段為2.3~5 mHz的自由振蕩信號,發(fā)現(xiàn)利用遠場體波反演的地震矩偏小.Hu等(2014)為了評估蘆山地震標量地震矩的正確性,對比了17個臺站和5個矩張量解的觀測值與模擬值,也得到了同樣的結論.目前,利用自由振蕩信號約束震源機制解所用到的觀測數(shù)據(jù)均來自超導重力儀(SG)和寬頻帶地震儀,當震級較大時,SG會發(fā)生限幅情況,因此在計算過程中必須去掉同震時的觀測段,對于MW8.1 Iquique地震來說,中國大陸的幾臺超導重力儀均出現(xiàn)了限幅情況,而寬頻帶地震儀主要用于記錄高頻信息.

    通過數(shù)字地震觀測網絡、中國大陸構造環(huán)境監(jiān)測網絡、中國地震背景場等國家重大科學工程項目,在中國大陸建成了一個集數(shù)據(jù)的自動匯集與數(shù)據(jù)庫式管理,儀器的遠程監(jiān)控和數(shù)據(jù)的圖形化分析處理為一體的連續(xù)重力觀測臺網,其平均空間分辨率約400 km.觀測儀器包括PET/gPhone、DZW、GS15、TRG-1和GWR超導等,除了SG以外,這些彈簧重力儀也可以觀測到微小的重力變化,分布廣泛且連續(xù)的彈簧重力儀觀測數(shù)據(jù)為檢測地球自由振蕩信號及利用自由振蕩信號約束震源機制解提供了很好的數(shù)據(jù)基礎.本文選取連續(xù)重力觀測臺網中分布廣泛的16個彈簧重力儀觀測數(shù)據(jù)(如圖1),利用自由振蕩信號觀測約束Iquique地震的震源機制解,并有效的評估哪一個震源機制解最符合地震實際釋放的能量,基于符合較好的重力站觀測資料,定量研究和分析影響自由振蕩振幅的因素.

    圖1 震中位置和16個連續(xù)重力觀測臺站的位置

    2數(shù)據(jù)和方法

    地球自由振蕩中的球型振蕩不僅有橫向分量,還有徑向分量,因此會引起地球的體積變化,從而導致重力發(fā)生變化.中國連續(xù)重力觀測臺網提供采樣率為1 Hz的重力觀測數(shù)據(jù)用于地球動力學研究.Iquique地震發(fā)生后,連續(xù)重力觀測臺網的重力儀清晰的記錄了重力變化.提取自由振蕩信號之前,首先對各臺站的觀測資料進行預處理,利用Tsoft(Vauterin,1998)重力潮汐數(shù)據(jù)預處理程序,去掉尖峰、突跳、間斷等.采用Eterna調和分析軟件,計算大氣導納因子,并在重力觀測中扣除大氣的影響(Wenzel,1996).

    r0處的點源地震激發(fā)臺站的簡正模位移可以利用雙力偶地震矩心矩張量理論和各項同性的地球模型模擬計算.在r處的位移u(r,t)為(Gilbert,1973):

    ×(M(r0):εn(r0)),

    (1)

    其中,Sn(r)為在某站點的標準化的n階模本征位移,ωn為本征頻率,Qn為質量因子,εn(r0)為應變張量,M(r0)為點源矩張量.二階6分量Mi(i=1,…,6),均勻地震矩張量M和εni組成εn(r0).則公式(1)可轉化為

    (2)

    (3)

    ai(r,t)為轉換函數(shù),Mi由標量地震矩M0、走向、傾角、滑動角組成(Aidetal.,1980).

    地球自由振蕩計算值的振幅主要取決于地震發(fā)生的位置和地震的震源機制解.模擬計算臺站自由振蕩需要的參數(shù)有:地震發(fā)生的時間、位置、走向、傾角、滑動角、標量地震矩、臺站的位置等.在參考的震源機制解中,地震的震源機制模型將震源等效為雙力偶點源,根據(jù)雙力偶地震矩心矩張量理論,由地球格林函數(shù)對矩張量加權求和,可計算矩張量源激發(fā)產生的地球自由振蕩位移(江穎等,2014b,2015).

    3結果分析

    3.1觀測值與模擬值比較

    將地震激發(fā)的地球自由振蕩信號的模擬值與相應的觀測值相比對,可評估震源機制解是否能正確反映地震的實際破裂情況.由于不同的地球模型對自由振蕩振幅的影響小于0.7%(薛秀秀等,2012),幾乎可以忽略.基于PREM地球模型(Dziewonski et al.,1981),結合震源機制解,將格林函數(shù)對矩張量加權求和,可以得到接收臺站的合成地震圖(Geller and Stein,1979).

    本文選取分布廣泛的16個彈簧重力儀觀測數(shù)據(jù)進行分析,臺站包括:北京(BJ)、福州(FZ)、高臺(GT)、格爾木(GEM)、海拉爾(HLE)、拉薩(LS)、蘭州(LZ)、牡丹江(MDJ)、憑祥(PX)、瓊中(QZ)、獅泉河(SQH)、炭山(TS)、烏加河(WJH)、烏什(WS)、西安(XA)、漳州(ZZ),見圖1.為了避免環(huán)境噪聲在簡正模振幅上的影響,選取16個臺站采樣間隔為1 s,時間長度為16h的重力觀測數(shù)據(jù),觀測數(shù)據(jù)和模擬數(shù)據(jù)在經過漢寧窗處理后進行離散傅里葉變換獲取線性振幅譜.連續(xù)觀測的彈簧重力儀會存在較大的零漂,選取的數(shù)據(jù)長度為16 h,零漂的影響較小,在數(shù)據(jù)處理過程中已經將該影響去掉.圖2顯示了16個重力臺站重力儀的自由振蕩觀測值和基于6個震源機制解計算的模擬值在1.5~5.3 mHz頻段的結果.在選取的頻段內,彈簧重力儀測得的自由振蕩信號有較高的信噪比.從圖2可以明顯看出基于震源機制解6的模擬值明顯偏小,而基于GCMT和USGS的震源機制解模擬計算的自由振蕩信號和觀測值符合的較好.

    3.2震源機制解整體評估

    為了更詳細的反映比較結果,基于16個重力儀臺站選用1.5~5.3 mHz頻段中所有的球型自由振蕩的模擬值和觀測值進行比較,利用殘余方差F(不符合度)和比例因子S(符合度)作為判斷標準:

    (4)

    (5)

    這里Oi和Si是觀測模的振幅和模擬值的振幅,n是所有振型的數(shù)量.當比例因子接近1的時候,其 震源機制解對應的M0是可靠的,可以描述震源的大小和強度.

    圖2 在1.5~5.3 mHz頻段內16個彈簧重力儀觀測數(shù)據(jù)與基于6個震源機制解算的模擬值

    另外,計算了基于各震源機制解結果的F和S的平均值及標準偏差,如表2.由表2可知,基于震源機制解3的自由振蕩模擬結果與實測結果最接近,其不符合度F的平均值為0.090,比例因子S為0.992,最接近1,F(xiàn)和S的標準偏差也最小,分別為0.055和0.065,最能反映地震釋放能量的大小.基于震源機制解1和2的模擬值與觀測值的殘余方差也較小,分別為0.098和0.091,顯示了較高的符合度.但是,基于震源機制解6的模擬值與觀測值相差較大,其不符合度F的平均值為0.214,比例因子S達到1.610.震源機制解中各個參數(shù)對地球自由振蕩信號振幅的影響需要進行進一步的計算與分析.

    表2 基于6個震源機制解的F與S的平均值和標準差

    3.3震源機制解影響定量分析

    分別計算震源機制解(標量地震矩、走向、傾角、滑動角和深度)對自由振蕩振幅的影響.基于USGSWPHASE的震源機制解計算的模擬值與觀測值最符合,由圖3的結果發(fā)現(xiàn),16個臺站中格爾木臺站觀測到的自由振蕩信號較接近模擬值,其對應的F為0.047,S為0.996.因此,基于格爾木臺站接收到的1.5~5.3 mHz頻段的自由振蕩信號平均振幅做定量分析.

    基于USGSWPHASE的震源機制解結果,分別給定標量地震矩、走向、傾角、滑動角和深度一個變化范圍,計算不同的震源機制解因素對自由振蕩振幅的影響.根據(jù)表1,標量地震矩的變化范圍為1.22×1021~2.35×1021N·m,保證走向、傾角、滑動角和深度不變,分別計算不同的標量地震矩對應的自由振蕩信號平均振幅,結果如圖4a,隨著標量地震矩的不斷增加,自由振蕩的平均振幅也在明顯增大,增量達到3.33×103nms-2/Hz;走向的變化范圍為240°~358°,保證M0、傾角、滑動角和深度不變,分別計算不同的走向對應的自由振蕩平均振幅,結果如圖4b,隨著走向增大,自由振蕩平均振幅也在慢慢增大,增加了2.25×103nms-2/Hz;傾角的變化范圍為12°~18°,保證M0、走向、滑動角和深度不變,分別計算不同的傾角對應的自由振蕩平均振幅,結果如圖4c,隨著傾角的不斷增加,振幅在慢慢減小,減小了3.08×103nms-2/Hz;滑動角變化范圍為74°~109°,保證M0、走向、傾角和深度不變,分別計算不同的滑動角對自由振蕩平均振幅的影響,結果如圖4d,隨著滑動角的不斷增加,振幅也在慢慢增大,增量達到2.13×103nms-2/Hz;深度變化范圍為20~33 km,保證M0、走向、傾角和滑動角不變,分別計算不同的深度對應的內核平動振蕩振幅的最大值,結果如圖4e,在20~33 km的范圍內,振幅有小范圍波動,差異僅為0.51×103nms-2/Hz.通過以上數(shù)值實驗說明M0對自由振蕩振幅的影響最大,而走向、傾角、滑動角和深度對內核平動振蕩振幅也有一定影響,但是影響相對較小.因此,將6個震源機制解中的M0均設置為相同的數(shù)值2.35×1021N·m,模擬計算格爾木臺站接收到的自由振蕩 信號,結果如圖5,當M0為2.35×1021N·m時,基于6個震源機制解得到的自由振蕩信號相差不大,且和觀測值均符合的較好.進一步證明了,M0是影響自由振蕩振幅的最大因素.An 等(2014)僅利用海嘯數(shù)據(jù)反演得到的震源機制解中M0偏小是造成模擬值和觀測值偏差較大的主要原因.

    圖4 震源機制解對自由振蕩振幅的影響

    圖5 格爾木臺重力觀測數(shù)據(jù)與基于改進的6個震源機制解解算的模擬值(其中,M0均為2.35×1021 N·m)

    4結論與討論

    通過比較分析Iquique地震激發(fā)的自由振蕩信號的實測值與模擬值,結果表明An 等(2014)等僅利用海嘯數(shù)據(jù)反演得到的震源機制解中M0偏小是造成模擬值和觀測值偏差較大的主要原因,說明對于本次的Iquique地震僅用海嘯數(shù)據(jù)做反演,有一定的局限性,這是由于海嘯數(shù)據(jù)對淺層地震的分辨率較高,而對于相對較深的逆沖型地震做數(shù)據(jù)分析的時候,并不一定能得到最佳解.利用海嘯數(shù)據(jù)反演震源機制主要是快速的得到破裂面如何破裂,在哪個方向上會有更大的破壞,而對于標量地震矩的估算會有一定偏差.另外,Lay 等(2014)從112個全球地震臺站選取200~1000 s的周期范圍的160個地面運動記錄的W-phase數(shù)據(jù)反演了Iquique地震,得到的標量地震矩為1.99×1021N·m,而另外的一組結合海嘯數(shù)據(jù)得到的震源機制解M0的為1.7×1021N·m,海嘯數(shù)據(jù)的加入造成了反演的標量地震矩結果偏小.Yagi等(2014)僅采用了遠震P波數(shù)據(jù),反演結果同樣偏小,可能是由于遠場體波記錄不到一些緩慢破裂的過程,且遠場體波容易受到其他震相的污染,特別是在主震后發(fā)生大量余震的情況下.對于大地震,USGS的地震檔案中一般不提供利用體波反演的震源機制解.USGSWPHASE選取長周期(周期為100~1000 s)方位角為5°~20°的近場數(shù)據(jù)和30°~90°的遠場數(shù)據(jù)聯(lián)合反演得到的震源機制解最佳.因此,利用長周期面波信號可以較好的約束Iquique地震的標量地震矩.

    本文基于地球自由振蕩分析了Iquique地震的震源機制解.結果表明利用彈簧重力儀記錄的1.5~5.3 mHz球型簡正模觀測可以有效地對Iquique地震的震源機制解進行總體評估.從研究結果來看,地震的斷層走向、傾角、滑動方向角和震源深度對自由振蕩的振幅有一定影響,但是影響較小,M0對自由振蕩振幅的影響最大.另外,基于USGSWPHASE反演得到的震源機制解的自由振蕩模擬值與觀測值符合最好,其反演的標量地震矩可以較好地反映Iquique地震釋放的能量,而僅利用海嘯數(shù)據(jù)反演的標量地震矩偏小,聯(lián)合近遠場觀測數(shù)據(jù)反演可顯著改善震源機制解.因此,利用自由振蕩的方法可以較好地判斷地震震源機制和有效約束地震的震級.

    另外,地球自由振蕩觀測只能對震源機制解進行總體綜合評估,而不能單獨對各個參數(shù)的影響單獨評估.如果能利用其他大地測量手段,如近場GPS觀測約束震源傾角、滑動方向角,結合地球自由振蕩的觀測和模擬可更細致地評估震源機制解.

    致謝感謝國家重力臺網中心提供的全國廣泛分布的連續(xù)重力觀測數(shù)據(jù).

    References

    AidK,Richards P G.1980.Quantitative Seismology : Theory and Methods.San Francisco: University Books.

    An C,Sepúlveda I,Liu P L F.2014.Tsunami source and its validation of the 2014 Iquique,Chile,earthquake.Geophysical Research Letters,41(11): 3988-3994.

    Angermann D,Klotz J,Reigber C.1999.Space-geodetic estimation of the Nazca-South America Euler vector.Earth and Planetary Science Letters,171(3): 329-334.

    Bürgmann R.2014.Earth science: Warning signs of the Iquique earthquake.Nature,512(7514): 258-259.

    Dziewonski A M,Anderson D L.1981.Preliminary reference Earth model.Phys.Earth Planet.Inter.,25(4): 297-356.

    Geller R J,Stein S.1979.Time-domain attenuation measurements for fundamental spheroidal modes (0S6to0S28) for the 1977 Indonesian earthquake.Bull.Seismol.Soc.Amer.,69(6):1671-1691.

    Gilbert F.1973.Derivation of source parameters from low-frequency spectra.Philosophical Transactions of the Royal Society A: Mathematical,Physical and Engineering Sciences,274(1239): 369-371.

    Hayes G P,Herman M W,Barnhart W D,et al.2014.Continuing megathrust earthquake potential in Chile after the 2014 Iquique earthquake.Nature,512(7514): 295-298.

    Hu X G,Jiang Y.2013.Evaluation of the seismic moment of the April 20,2013 Lushan earthquake.Earthq.Sci.,26(3-4): 169-177.

    Hu X G,Jiang Y,Sun H P.2014.Assessing the scalar moment of moderate earthquake and the effect of lateral heterogeneity on normal modes—An example from the 2013/04/20 Lushan earthquake,Sichuan,China.Physics of the Earth and Planetary Interiors,232: 61-71.

    Jiang Y,Hu X G,Liu C L,et al.2014.Constraining the focal mechanism of the Lushan earthquake with observations of the Earth′s free oscillations.Science China Earth Sciences,57(9): 2064-2070.

    Kato A,Nakagawa S.2014.Multiple slow-slip events during a foreshock sequence of the 2014 Iquique,Chile MW8.1 earthquake.Geophysical Research Letters,41(15): 5420-5427.

    Lambotte S,Rivera L,Hinderer J.2006.Rupture length and duration of the 2004 Aceh-Sumatra earthquake from the phases of the Earth′s gravest free oscillations.Geophysical Research Letters,33(3): L03307.

    Lay T,Yue H,Brodsky E E,et al.2014.The 1 April 2014 Iquique,Chile,MW8.1 earthquake rupture sequence.Geophysical Research Letters,41(11): 3818-3825.

    Liu C L,Zheng Y,Wang R J,et al.2015.Kinematic rupture process of the 2014 Chile MW8.1 earthquake constrained by strong-motion,GPS static offsets and teleseismic data.Geophysical Journal International,2: 1137-1145.

    Park J,Song T-R A,Tromp J,et al.2005.Earth′s free oscillations excited by the 26 December 2004 Sumatra-Andaman earthquake.Science,308(5725): 1139-1144.Ruiz S,Metois M,Fuenzalida A,et al.2014.Intense foreshocks and a slow slip event preceded the 2014 Iquique MW8.1 earthquake.Science,345(6201): 1165-1169.

    Schurr B,Asch G,Hainzl S,et al.2014.Gradual unlocking of plate boundary controlled initiation of the 2014 Iquique earthquake.Nature,512(7514): 299-302.

    Stein S,Okal E A.2005.Seismology: Speed and size of the Sumatra earthquake.Nature,308:1127-1139.

    Vauterin P.1998.Tsoft: Graphical and interactive software for the analysis of Earth tide data.∥ Proc.13th Int.Sympos.on Earth Tides,Brussels,Observatoire Royal de Belgique,Série Géophysique.481-486.

    Wenzel H G.1996.The nanogal software: Data processing package ETERNA 3.3.Bull.Inf.Marées Terrestres,124: 9425-9439.

    Yagi Y,Okuwaki R,Enescu B,et al.2014.Rupture process of the 2014 Iquique Chile Earthquake in relation with the foreshock activity.Geophysical Research Letters,41(12): 4201-4206.

    附中文參考文獻

    江穎,胡小剛,劉成利等.2014a.利用地球自由振蕩觀測約束蘆山地震的震源機制解.中國科學: 地球科學,44(12): 2689-2696.

    江穎,徐建橋,孫和平.2014b.深內部地球結構對內核平動振蕩本征周期的影響.地球物理學報,57(4): 1041-1048,doi: 10.6038/cjg20140403.

    江穎,徐建橋,孫和平.2015.地球內核平動振蕩的地震激發(fā).地球物理學報,58(1): 92-102,doi: 10.6038/cjg20150108.

    薛秀秀,胡小剛,郝曉光等.2012.利用重力觀測約束2011日本Tohoku大地震的震源機制.地球物理學報,55(9): 3006-3015,doi: 10.6038/j.issn.0001-5733.2012.09.019.

    (本文編輯汪海英)

    基金項目科研院所基本科研業(yè)務費(中國地震局地震研究所所長基金)(IS201526225),大地測量與地球動力學國家重點實驗室開放基金(SKLGED2015-1-2-EZ),國家自然科學基金(41404064,41374088)資助.

    作者簡介江穎,女,1986年生,現(xiàn)主要從事地球自由振蕩、重力數(shù)據(jù)分析及其應用研究.E-mail:jiangyingchen@126.com

    doi:10.6038/cjg20160316 中圖分類號P315

    收稿日期2015-04-08,2016-01-11收修定稿

    Constraining the focal mechanism of the Iquique earthquake with observations of the continuous gravity stations

    JIANG Ying1,2,LI Hui1,LIU Zi-Wei1,SHEN Chong-Yang1,ZHANG Xiao-Tong1

    1KeyLaboratoryofEarthquakeGeodesy,InstituteofSeismology,ChinaEarthquakeAdministration,Wuhan430071,China2StateKeyLaboratoryofGeodesyandEarth′sDynamics,InstituteofGeodesyandGeophysics,CAS,Wuhan430077,China

    AbstractOn 1 April 2014,a magnitude MW8.1 interplate thrust earthquake ruptured a densely instrumented region of Iquique seismic gap in northern Chile.After Earthquake,the focal mechanism solutions of Iquique earthquake were provided by well-respected international earthquake research institutions and scholars based on different data and methods,which were quite different.The amplitudes of the Earth′s free oscillations have a close relationship to earthquake focal mechanisms.Focal mechanisms of large earthquakes can be well analyzed and constrained with observations of long period free oscillations.We compared free oscillations observed by 16 spring gravimeters of continuous gravity stations with synthetic normal modes corresponding to six different focal mechanisms for the Iquique earthquake,and the focal mechanisms solutions of Iquique earthquake were analyzed and constrained by spherical normal modes in a 1.5 to 5.3 mHz frequency band.It′s found that the synthetic modes corresponding to the focal mechanism determined by the United States Geological Survey WPhase Moment Tensor Solution show agreement to the observed modes,suggesting that earthquake magnitudes predicted in this way can reflect the total energy released by the earthquake.The scalar seismic moment obtained by tsunami data inversion is significantly underestimated.Focal mechanism solutions can be improved by joint inversion of far- and near-field data.In addition,based on the synthetic and observed normal modes at Golmud station,quantitative analysis was carried out to research the influences of free oscillation′s amplitudes caused by different focal mechanism parameters.The results show that scalar seismic moment M0 has a major influence on the amplitudes of free oscillations,and the strike,dip,rake and depth of the hypocenter have minor influences.

    KeywordsIquique earthquake; Focal mechanism solutions; Earth′s free oscillations; Gravimeter observations

    江穎,李輝,劉子維等.2016.利用連續(xù)重力觀測約束2014智利Iquique地震的震源機制解.地球物理學報,59(3):941-951,doi:10.6038/cjg20160316.

    Jiang Y,Li H,Liu Z W,et al.2016.Constraining the focal mechanism of the Iquique earthquake with observations of the continuous gravity stations.Chinese J.Geophys.(in Chinese),59(3):941-951,doi:10.6038/cjg20160316.

    猜你喜歡
    標量臺站震源
    中國科學院野外臺站檔案工作回顧
    氣象基層臺站建設
    西藏科技(2021年12期)2022-01-17 08:46:38
    一種高效的橢圓曲線密碼標量乘算法及其實現(xiàn)
    一種靈活的橢圓曲線密碼并行化方法
    震源的高返利起步
    可控震源地震在張掖盆地南緣逆沖斷裂構造勘探中的應用
    華北地質(2015年3期)2015-12-04 06:13:25
    基層臺站綜合觀測業(yè)務管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    同步可控震源地震采集技術新進展
    單調Minkowski泛函與Henig真有效性的標量化
    MDOS平臺臺站級使用方法及技巧
    在线观看av片永久免费下载| 欧美人与善性xxx| 成人无遮挡网站| 国产精品三级大全| 成人毛片a级毛片在线播放| 国产精品人妻久久久久久| 91精品国产九色| 丰满少妇做爰视频| 亚洲国产精品专区欧美| 亚洲激情五月婷婷啪啪| 日韩,欧美,国产一区二区三区| 久久精品熟女亚洲av麻豆精品| 国产 精品1| 啦啦啦视频在线资源免费观看| 精品一区二区免费观看| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久久久电影| 中文字幕久久专区| 如何舔出高潮| 在线 av 中文字幕| 久久 成人 亚洲| 十分钟在线观看高清视频www | 久久综合国产亚洲精品| 国产精品爽爽va在线观看网站| 国产精品无大码| 又大又黄又爽视频免费| 交换朋友夫妻互换小说| 精品熟女少妇av免费看| 女性生殖器流出的白浆| av视频免费观看在线观看| 97在线人人人人妻| 亚洲国产精品国产精品| 精品99又大又爽又粗少妇毛片| 国产精品一二三区在线看| 亚洲国产精品成人久久小说| 啦啦啦视频在线资源免费观看| 在线观看美女被高潮喷水网站| 日韩强制内射视频| 国产精品国产三级国产av玫瑰| 精品人妻视频免费看| 亚洲av日韩在线播放| 国产有黄有色有爽视频| 黄色一级大片看看| 国产成人一区二区在线| 男人添女人高潮全过程视频| 亚洲精品乱码久久久v下载方式| 中文字幕亚洲精品专区| 中文字幕亚洲精品专区| 国产精品久久久久久精品古装| 日韩 亚洲 欧美在线| 日韩亚洲欧美综合| 亚洲精品国产av成人精品| 久久97久久精品| 国产成人免费无遮挡视频| 免费观看性生交大片5| 欧美另类一区| 日韩中文字幕视频在线看片 | 狂野欧美白嫩少妇大欣赏| 一二三四中文在线观看免费高清| 久久精品久久久久久噜噜老黄| 成人毛片a级毛片在线播放| 亚洲电影在线观看av| 国产精品一区二区三区四区免费观看| 亚洲第一av免费看| 熟妇人妻不卡中文字幕| 色视频www国产| 亚洲欧美清纯卡通| 亚洲欧洲日产国产| 少妇人妻 视频| 久久热精品热| 三级经典国产精品| 熟女av电影| 少妇精品久久久久久久| 久久精品国产亚洲av天美| 水蜜桃什么品种好| 精品久久久久久久末码| 色视频在线一区二区三区| 国产成人精品婷婷| 久久国产乱子免费精品| 中文天堂在线官网| 夫妻性生交免费视频一级片| 国产在线视频一区二区| h视频一区二区三区| 国产高清三级在线| 久久久久久久久久久丰满| 色婷婷久久久亚洲欧美| 亚洲精品乱久久久久久| 日本色播在线视频| 亚洲成人av在线免费| 97超碰精品成人国产| 国产深夜福利视频在线观看| 麻豆成人午夜福利视频| 欧美高清性xxxxhd video| 最近中文字幕高清免费大全6| 永久免费av网站大全| 高清午夜精品一区二区三区| 国产精品国产av在线观看| 女性生殖器流出的白浆| 永久免费av网站大全| 日日撸夜夜添| 国产精品久久久久久av不卡| 国产69精品久久久久777片| 免费黄频网站在线观看国产| 久久久久性生活片| 国产欧美日韩一区二区三区在线 | 亚洲国产精品专区欧美| 日韩中文字幕视频在线看片 | 亚洲三级黄色毛片| 乱码一卡2卡4卡精品| 51国产日韩欧美| 亚洲av.av天堂| 日本黄大片高清| 五月玫瑰六月丁香| 国产av一区二区精品久久 | 国产男人的电影天堂91| 秋霞伦理黄片| av线在线观看网站| 80岁老熟妇乱子伦牲交| 在线播放无遮挡| 亚洲欧美日韩东京热| 亚洲中文av在线| 亚洲av成人精品一区久久| 尾随美女入室| 老司机影院毛片| 久久久亚洲精品成人影院| 少妇人妻 视频| 多毛熟女@视频| 午夜福利视频精品| 免费黄色在线免费观看| 亚洲婷婷狠狠爱综合网| 亚洲国产精品一区三区| 黑人猛操日本美女一级片| 2018国产大陆天天弄谢| 国产精品一及| 91午夜精品亚洲一区二区三区| 婷婷色综合大香蕉| 18禁裸乳无遮挡动漫免费视频| 国产精品一区二区在线不卡| 亚洲一级一片aⅴ在线观看| 欧美精品一区二区免费开放| 两个人的视频大全免费| 国产免费福利视频在线观看| 成人毛片60女人毛片免费| 美女主播在线视频| 国产精品免费大片| 亚洲欧美日韩卡通动漫| 51国产日韩欧美| 老女人水多毛片| www.av在线官网国产| 国产免费一级a男人的天堂| 最新中文字幕久久久久| 精品久久国产蜜桃| 一本一本综合久久| 婷婷色av中文字幕| 精品视频人人做人人爽| 涩涩av久久男人的天堂| 免费人成在线观看视频色| 亚洲成人一二三区av| 欧美3d第一页| 亚洲成人手机| 人人妻人人添人人爽欧美一区卜 | 日本vs欧美在线观看视频 | 久久女婷五月综合色啪小说| 人妻制服诱惑在线中文字幕| 久久毛片免费看一区二区三区| 高清视频免费观看一区二区| 国产一区二区三区综合在线观看 | 国产成人精品福利久久| 男的添女的下面高潮视频| 高清日韩中文字幕在线| 日韩av免费高清视频| 亚洲高清免费不卡视频| 妹子高潮喷水视频| 亚洲av中文字字幕乱码综合| 看十八女毛片水多多多| 伦精品一区二区三区| 国产探花极品一区二区| 精品久久久噜噜| 午夜免费男女啪啪视频观看| 欧美区成人在线视频| av在线蜜桃| av免费观看日本| 黄片无遮挡物在线观看| 视频区图区小说| 日韩伦理黄色片| 日本vs欧美在线观看视频 | 久久国产乱子免费精品| 日韩不卡一区二区三区视频在线| 久久久久久久久大av| 自拍欧美九色日韩亚洲蝌蚪91 | 偷拍熟女少妇极品色| 91精品伊人久久大香线蕉| 嫩草影院入口| 亚洲国产色片| 国产一级毛片在线| 男人添女人高潮全过程视频| 免费大片黄手机在线观看| 午夜福利网站1000一区二区三区| 久久久久久久久久久免费av| 99精国产麻豆久久婷婷| 纵有疾风起免费观看全集完整版| 中国三级夫妇交换| 91精品国产国语对白视频| 妹子高潮喷水视频| 中文字幕久久专区| 久久久久久久久久人人人人人人| 欧美日韩视频精品一区| 中国国产av一级| xxx大片免费视频| 在线观看免费视频网站a站| 亚洲人成网站在线观看播放| 女人久久www免费人成看片| 国产黄频视频在线观看| 成人国产av品久久久| 久久精品夜色国产| 精品99又大又爽又粗少妇毛片| 欧美一级a爱片免费观看看| 少妇高潮的动态图| 2022亚洲国产成人精品| 777米奇影视久久| 好男人视频免费观看在线| 欧美高清性xxxxhd video| 免费看光身美女| 一区二区三区精品91| 成人国产麻豆网| 国产亚洲5aaaaa淫片| 国产永久视频网站| 中文字幕av成人在线电影| 身体一侧抽搐| 这个男人来自地球电影免费观看 | 欧美激情国产日韩精品一区| av福利片在线观看| 日日啪夜夜撸| 夫妻午夜视频| 2021少妇久久久久久久久久久| 欧美高清性xxxxhd video| 看免费成人av毛片| 少妇人妻一区二区三区视频| 有码 亚洲区| 在线免费十八禁| 日本欧美视频一区| 亚洲中文av在线| 精品国产乱码久久久久久小说| 草草在线视频免费看| 这个男人来自地球电影免费观看 | 熟女电影av网| 国产极品天堂在线| 在线观看一区二区三区| 亚洲高清免费不卡视频| 精品人妻熟女av久视频| 亚洲美女搞黄在线观看| 一级毛片aaaaaa免费看小| 欧美激情极品国产一区二区三区 | 大话2 男鬼变身卡| 美女主播在线视频| 国产精品一区二区在线不卡| 一本久久精品| a 毛片基地| 久久久久久久大尺度免费视频| 国产成人a区在线观看| 国产精品99久久久久久久久| 在线天堂最新版资源| 亚洲欧美精品专区久久| av又黄又爽大尺度在线免费看| 毛片一级片免费看久久久久| 内地一区二区视频在线| 国产精品一区二区三区四区免费观看| 日韩伦理黄色片| 天天躁日日操中文字幕| 国产深夜福利视频在线观看| 国产白丝娇喘喷水9色精品| 日本wwww免费看| 六月丁香七月| 99热这里只有精品一区| 成人一区二区视频在线观看| 国产精品一区二区三区四区免费观看| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩在线观看h| 亚洲av福利一区| 国产精品久久久久成人av| 一个人免费看片子| 成人亚洲精品一区在线观看 | 国产一区二区在线观看日韩| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品中文字幕在线视频 | 久久久久久久大尺度免费视频| 欧美老熟妇乱子伦牲交| 国产极品天堂在线| 欧美日韩视频高清一区二区三区二| 国产视频内射| 欧美+日韩+精品| 国产精品福利在线免费观看| 观看美女的网站| 婷婷色综合www| 亚洲精品国产av蜜桃| 99热全是精品| 在线免费观看不下载黄p国产| 久久人人爽人人片av| 建设人人有责人人尽责人人享有的 | 国产精品不卡视频一区二区| 国产伦精品一区二区三区四那| 国产成人精品福利久久| 美女脱内裤让男人舔精品视频| 精品人妻偷拍中文字幕| 日韩大片免费观看网站| av不卡在线播放| 男人添女人高潮全过程视频| 亚洲国产成人一精品久久久| 色5月婷婷丁香| 成人黄色视频免费在线看| 亚洲色图综合在线观看| 欧美 日韩 精品 国产| 蜜桃在线观看..| 亚洲欧美日韩无卡精品| 卡戴珊不雅视频在线播放| 欧美成人午夜免费资源| 伊人久久精品亚洲午夜| 五月伊人婷婷丁香| 色婷婷av一区二区三区视频| 成人漫画全彩无遮挡| 啦啦啦视频在线资源免费观看| 久久久久网色| 精品人妻一区二区三区麻豆| 国产精品人妻久久久久久| 两个人的视频大全免费| 最新中文字幕久久久久| 日本爱情动作片www.在线观看| 三级国产精品片| av专区在线播放| 亚洲精品aⅴ在线观看| 中文精品一卡2卡3卡4更新| 欧美成人一区二区免费高清观看| 国产爽快片一区二区三区| 国产一区二区三区综合在线观看 | 色网站视频免费| 男女下面进入的视频免费午夜| 青青草视频在线视频观看| 三级国产精品片| 尾随美女入室| 久久久亚洲精品成人影院| 国产男人的电影天堂91| 国产精品久久久久成人av| 国产成人a区在线观看| 国产乱人视频| 亚洲精品国产成人久久av| 亚洲精品中文字幕在线视频 | 九九爱精品视频在线观看| 亚洲中文av在线| 天天躁夜夜躁狠狠久久av| 精品久久久久久久久av| 亚洲欧美清纯卡通| 中文字幕免费在线视频6| 亚洲欧美日韩东京热| 18禁裸乳无遮挡动漫免费视频| 熟女人妻精品中文字幕| 久久av网站| 国产精品秋霞免费鲁丝片| 久久97久久精品| a级毛色黄片| 日本欧美视频一区| 高清午夜精品一区二区三区| 国产精品秋霞免费鲁丝片| 欧美少妇被猛烈插入视频| 一级片'在线观看视频| 国产无遮挡羞羞视频在线观看| 精品国产露脸久久av麻豆| 国国产精品蜜臀av免费| 一级黄片播放器| 大香蕉久久网| 国产一区亚洲一区在线观看| 国产 精品1| 人妻一区二区av| 久久99蜜桃精品久久| 日本av手机在线免费观看| 欧美高清性xxxxhd video| 国产在线视频一区二区| 七月丁香在线播放| 麻豆乱淫一区二区| 夫妻性生交免费视频一级片| 成人漫画全彩无遮挡| 精品亚洲成a人片在线观看 | 卡戴珊不雅视频在线播放| 亚洲三级黄色毛片| 日本av免费视频播放| 18禁在线无遮挡免费观看视频| 99久久精品一区二区三区| 日韩制服骚丝袜av| 亚洲国产精品成人久久小说| 日韩在线高清观看一区二区三区| 九九久久精品国产亚洲av麻豆| 亚洲四区av| 国产在线视频一区二区| 日韩中文字幕视频在线看片 | 亚洲精品亚洲一区二区| 国产免费一区二区三区四区乱码| 99久久精品一区二区三区| av国产免费在线观看| 午夜老司机福利剧场| 99热这里只有是精品50| 交换朋友夫妻互换小说| 综合色丁香网| 精品久久久噜噜| 欧美日韩精品成人综合77777| 最后的刺客免费高清国语| 亚洲丝袜综合中文字幕| 精品国产三级普通话版| 久久99热6这里只有精品| 亚洲国产高清在线一区二区三| 国产精品女同一区二区软件| 国产高清三级在线| 日本av免费视频播放| 一本一本综合久久| 国产在线免费精品| 国产亚洲欧美精品永久| 亚洲成人一二三区av| 插逼视频在线观看| 天美传媒精品一区二区| 男人添女人高潮全过程视频| 少妇被粗大猛烈的视频| 亚洲欧美中文字幕日韩二区| .国产精品久久| 成人影院久久| 熟妇人妻不卡中文字幕| av网站免费在线观看视频| 女人十人毛片免费观看3o分钟| 乱码一卡2卡4卡精品| 欧美区成人在线视频| 免费不卡的大黄色大毛片视频在线观看| 人妻少妇偷人精品九色| 日本黄色日本黄色录像| 国产精品三级大全| 久久久久久久国产电影| 久久国产亚洲av麻豆专区| 久久久久久久久久久丰满| 久久99蜜桃精品久久| 夜夜骑夜夜射夜夜干| 国产午夜精品一二区理论片| 搡女人真爽免费视频火全软件| 人妻少妇偷人精品九色| 夫妻性生交免费视频一级片| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久久av不卡| 欧美+日韩+精品| 精品久久久久久久久av| 国产一区有黄有色的免费视频| 一级毛片 在线播放| 性色av一级| 国产精品.久久久| 亚洲欧洲日产国产| 老女人水多毛片| 精品久久久久久久久亚洲| 狂野欧美激情性bbbbbb| 色综合色国产| 成人18禁高潮啪啪吃奶动态图 | 夜夜看夜夜爽夜夜摸| 日日啪夜夜撸| 亚洲va在线va天堂va国产| 不卡视频在线观看欧美| 精品一区在线观看国产| 成人午夜精彩视频在线观看| 日日撸夜夜添| 国产精品久久久久久精品古装| 国产女主播在线喷水免费视频网站| 99久久综合免费| 国产午夜精品一二区理论片| 欧美日韩精品成人综合77777| 乱系列少妇在线播放| 国产精品一区二区三区四区免费观看| 免费高清在线观看视频在线观看| 黄色怎么调成土黄色| 99热这里只有精品一区| 国产欧美亚洲国产| 精品午夜福利在线看| 晚上一个人看的免费电影| 交换朋友夫妻互换小说| 一二三四中文在线观看免费高清| 精品视频人人做人人爽| 少妇人妻 视频| 国产视频内射| 高清黄色对白视频在线免费看 | 国产亚洲午夜精品一区二区久久| 国产无遮挡羞羞视频在线观看| 熟女电影av网| 偷拍熟女少妇极品色| 老女人水多毛片| 欧美精品一区二区免费开放| 免费观看的影片在线观看| 欧美人与善性xxx| 国产男人的电影天堂91| 久久久a久久爽久久v久久| 少妇被粗大猛烈的视频| 欧美丝袜亚洲另类| 日本免费在线观看一区| 激情五月婷婷亚洲| 亚洲不卡免费看| 国内少妇人妻偷人精品xxx网站| 一区二区三区四区激情视频| 国产色婷婷99| 少妇被粗大猛烈的视频| 乱系列少妇在线播放| 国产有黄有色有爽视频| 青青草视频在线视频观看| 国产色婷婷99| 亚洲国产av新网站| 国产精品久久久久久久久免| 久久99精品国语久久久| 在线免费十八禁| 免费观看在线日韩| 欧美成人一区二区免费高清观看| 大香蕉久久网| 91久久精品国产一区二区三区| 高清在线视频一区二区三区| 免费人妻精品一区二区三区视频| 2021少妇久久久久久久久久久| 黄色日韩在线| 久久 成人 亚洲| 欧美激情极品国产一区二区三区 | 嫩草影院入口| 国产精品.久久久| 国产成人免费无遮挡视频| 国产欧美亚洲国产| 美女xxoo啪啪120秒动态图| 男女无遮挡免费网站观看| 99热网站在线观看| 国产乱人偷精品视频| 草草在线视频免费看| 国产乱人偷精品视频| 日日撸夜夜添| 老师上课跳d突然被开到最大视频| 少妇人妻久久综合中文| 99久久人妻综合| 亚洲国产精品专区欧美| 久久精品国产亚洲av天美| 亚洲熟女精品中文字幕| 国产乱来视频区| 男的添女的下面高潮视频| 午夜老司机福利剧场| 人妻制服诱惑在线中文字幕| 成人影院久久| 人妻 亚洲 视频| 精品久久久久久久久av| 精品一品国产午夜福利视频| 亚洲自偷自拍三级| av在线老鸭窝| 一本—道久久a久久精品蜜桃钙片| 精品人妻一区二区三区麻豆| 国内揄拍国产精品人妻在线| 如何舔出高潮| 国产真实伦视频高清在线观看| 亚洲欧美日韩另类电影网站 | 九九久久精品国产亚洲av麻豆| 特大巨黑吊av在线直播| 国产 一区精品| 久久精品久久久久久久性| 精品人妻偷拍中文字幕| 国产人妻一区二区三区在| 久久久亚洲精品成人影院| 久久人人爽人人片av| 亚洲精品亚洲一区二区| 秋霞在线观看毛片| 欧美亚洲 丝袜 人妻 在线| 一级黄片播放器| 在线观看免费日韩欧美大片 | 亚洲精品视频女| 少妇熟女欧美另类| 久久99蜜桃精品久久| 99热全是精品| 亚洲av国产av综合av卡| 最近的中文字幕免费完整| 91午夜精品亚洲一区二区三区| 人妻系列 视频| 日本av免费视频播放| 美女视频免费永久观看网站| 久久国产乱子免费精品| 特大巨黑吊av在线直播| 在线观看一区二区三区| 亚洲国产精品国产精品| 我的女老师完整版在线观看| 久久精品国产鲁丝片午夜精品| 精品久久久久久电影网| 特大巨黑吊av在线直播| 精华霜和精华液先用哪个| 国产精品久久久久久久久免| 少妇被粗大猛烈的视频| 欧美成人午夜免费资源| 99热网站在线观看| 最近最新中文字幕大全电影3| av又黄又爽大尺度在线免费看| 亚洲av国产av综合av卡| 丝袜脚勾引网站| 高清黄色对白视频在线免费看 | 毛片一级片免费看久久久久| h日本视频在线播放| 欧美成人a在线观看| 夫妻性生交免费视频一级片| 亚洲成人手机| 全区人妻精品视频| 国产精品嫩草影院av在线观看| 男的添女的下面高潮视频| 少妇熟女欧美另类| 91aial.com中文字幕在线观看| 国产无遮挡羞羞视频在线观看| 国产午夜精品久久久久久一区二区三区| 99热这里只有是精品在线观看| 亚洲精品,欧美精品| 啦啦啦啦在线视频资源| 日本与韩国留学比较| 韩国高清视频一区二区三区| 国产午夜精品久久久久久一区二区三区| 欧美3d第一页| 黄片wwwwww| 亚洲国产精品国产精品| 爱豆传媒免费全集在线观看| 国产日韩欧美在线精品| 久久久精品94久久精品|