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

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx

    2014年魯?shù)榈卣?MS=6.5)靜態(tài)庫侖應力變化及其影響

    2016-09-02 06:12:07朱守彪
    地震地質(zhì) 2016年1期
    關(guān)鍵詞:魯?shù)?/a>主震庫侖

    繆 淼 朱守彪

    1)中國地震局地殼應力研究所、地殼動力學重點實驗室、北京 100085 2)中國地震局地球物理研究所、北京 100081

    ?

    2014年魯?shù)榈卣?MS=6.5)靜態(tài)庫侖應力變化及其影響

    繆淼1,2)朱守彪1)*

    1)中國地震局地殼應力研究所、地殼動力學重點實驗室、北京1000852)中國地震局地球物理研究所、北京100081

    2014年8月3日在云南省昭通市魯?shù)榭h發(fā)生了MS6.5左旋走滑型地震。為了解本次地震所造成的影響、文中利用前人反演的震源模型計算了主震產(chǎn)生的同震靜態(tài)庫侖應力變化、考察其對1個月內(nèi)余震的觸發(fā)效果以及對周圍斷層的影響。經(jīng)計算后發(fā)現(xiàn)、82.43%的后續(xù)余震發(fā)生在靜態(tài)庫侖應力增加>0.01MPa的區(qū)域;說明大部分余震是由主震觸發(fā)產(chǎn)生的。然后、使用周圍斷層作為接收斷層、計算了魯?shù)榈卣鹪谥苓厰鄬由袭a(chǎn)生的庫侖應力變化:發(fā)現(xiàn)昭通-魯?shù)閿嗔?西支)的東北端及該斷層東支的東北段庫侖應力增大、發(fā)生地震的危險性會增大;而該區(qū)的大涼山斷裂、蓮峰斷裂、則木河斷裂、小江斷裂和馬邊-鹽津斷裂的庫侖應力都減小、即未來發(fā)生地震的危險性下降。另外、計算還發(fā)現(xiàn)、當選擇最優(yōu)破裂面作為庫侖應力變化的投影面時、無論使用不同的震源模型還是改變有效摩擦系數(shù)、魯?shù)榈卣饘ζ浜罄m(xù)余震始終具有較好的觸發(fā)作用。

    2014年魯?shù)榈卣痨o態(tài)庫侖應力變化余震觸發(fā)

    0 引言

    北京時間2014年8月3日16時30分、云南省昭通市魯?shù)榭h發(fā)生MS6.5地震(以下簡稱 “魯?shù)榈卣稹?;震中位置為27.08°N、103.37°E;震源深度13.3km(張廣偉等、2014;圖1)。本次地震發(fā)生在滇東北、南北地震帶中南段、那里屬于青藏高原東南部川滇菱形塊體與華南地塊之間的Ⅰ級塊體邊界構(gòu)造帶(Xuetal.、2003; Zhangetal.、2003)。雖然魯?shù)榈卣饹]有發(fā)現(xiàn)明顯的地表破裂帶、但據(jù)現(xiàn)有的資料可知、發(fā)震斷層并非該區(qū)的主要構(gòu)造昭通-魯?shù)閿嗔选⒍荖W向的包谷垴-小河斷裂、屬于鮮水河-小江斷裂系東側(cè)大涼山斷裂南端的組成部分、與馬邊-鹽津斷裂一起組成了由青藏高原最東緣與華南地塊相互作用形成的前緣最新構(gòu)造變形帶、分解了鮮水河-小江斷裂系安寧河斷裂和則木河斷裂的左旋走滑分量(徐錫偉等、2014)。

    圖1 2014云南魯?shù)镸S6.5地震震源位置、周圍強震分布(a)與震源機制解(b)(來源于:Global CMT)Fig. 1 (a)Location of hypercenter of the 2014 MS6.5 Ludian earthquake and strong shocks around the area; (b)Focal mechanism of the 2014 MS6.5 Ludian earthquake(from: Global CMT)

    一次地震之后、主震對余震是否具有觸發(fā)作用以及周邊斷層的地震危險性是人們關(guān)注的問題、分析一次地震前后周圍地區(qū)的應力變化對判斷未來的震情具有一定的作用。近年來大量學者使用庫侖應力變化(ΔCFS)來考察震后的應力狀態(tài)、從而考察對余震的觸發(fā)作用以及周邊斷層的加、卸載情況(Kingetal.、1994; Harris、1998; Stein、1999; Todaetal.、2000; 萬永革等、2000; 沈正康等、2003; Linetal.、2004; McCloskeyetal.、2005; 陳連旺等、2008; 張竹琪等、2008; 單斌等、2012; 繆淼等、2012、2013)。King等(1994)認為使用 “最優(yōu)破裂面”作為靜態(tài)庫侖應力變化的投影面、得到的結(jié)果可以用來解釋余震的分布、預測未來的地震危險性。對此、他們計算了1992年美國LandersMW7.3地震造成的破裂面附近最優(yōu)方向上的庫侖應力變化、發(fā)現(xiàn)余震廣泛分布于庫侖應力增加0.01MPa的區(qū)域、而庫侖應力降低的區(qū)域余震活動較少。Mallman等(2007)進一步使用各種不同的模型對Landers地震和1995年日本Kobe地震進行驗證、最終發(fā)現(xiàn)以最優(yōu)破裂面投影得到的庫侖應力上升區(qū)域、分布了更多的余震。Ma等(2005)計算了1999年臺灣集集ML7.3地震后庫侖應力變化在最優(yōu)破裂面上的投影、發(fā)現(xiàn)余震活動與其有極好的對應關(guān)系。2008年汶川MS8.0地震后的研究成果(Nalbantetal.、2010; Wanetal.、2010; Xuetal.、2010; 解朝娣等、2010)也表明、主震產(chǎn)生的靜態(tài)庫侖應力變化對后續(xù)余震有較好的觸發(fā)效果。

    那么、本次魯?shù)榈卣饘ζ溆嗾鹗欠窬哂杏|發(fā)作用呢?未來哪些地區(qū)的地震危險性上升了呢?周圍斷層受到的應力擾動情況又是如何呢?為此、本研究利用前人反演的震源斷層模型、計算魯?shù)榈卣甬a(chǎn)生的靜態(tài)庫侖應力變化、考察其對余震的觸發(fā)作用、并計算周圍斷層面的庫侖應力變化、判斷該區(qū)未來的地震危險性。

    1 靜態(tài)庫侖應力變化

    根據(jù)庫侖破裂假設、巖石趨近于破裂的庫侖破裂應力σf為(Kingetal.、1994)

    (1)

    式(1)中、τ為地震破裂面上剪應力的大小、σn為正應力、P為孔隙流體壓力、μ為斷層面介質(zhì)的摩擦系數(shù)。定義壓應力為正。

    然而、精確確定地下應力張量是極為困難的、通常使用庫侖破裂應力變化(Harris、1998)。當μ不隨時間變化時、由式(1)庫侖應力變化為

    (2)

    孔隙流體壓力變化ΔP控制著斷層面上的有效正應力、當巖石應力的改變遠遠快于巖石中的流體壓力擴散時、流體壓力變化ΔP可以通過Skemptons系數(shù)B在式(2)中得到反映。取視摩擦系數(shù)μ′=μ(1-B)、它給出了孔隙流體和斷層面上的介質(zhì)特性、范圍為0~1。那么式(2)變?yōu)?/p>

    (3)

    庫侖破裂應力變化需要定義于具體的斷層面。通常、以最優(yōu)破裂面投影得到的庫侖應力變化、可以解釋余震的分布情況、預測未來的后續(xù)地震活動(Kingetal.、1994; Harrisetal.、1996; Maetal.、2005; Todaetal.、2005; Xuetal.、2010)。而所謂最優(yōu)破裂面是指計算庫侖應力變化時、某一產(chǎn)狀接收斷層面上計算得到的庫侖應力變化幅值大于同一地點其他任意產(chǎn)狀的接收斷層面得到的應力變化幅值。

    (4)

    圖2 最優(yōu)取向斷層面上應力的坐標系示意圖(King et al.、1994)Fig. 2 The axis system used for calculations of Coulomb stresses on optimum failure planes(after King et al.、1994).

    (5)

    (6)

    (7)

    求得σ33和τ13的相對變化Δσ33和Δτ13、進而最優(yōu)取向投影方向上的庫侖應力為

    (8)

    2 資料與計算結(jié)果

    2.1余震時空分布特征

    圖3 余震的空間展布、深度分布(a)與隨時間衰減情況及擬合曲線(b)Fig. 3 (a)Spatial and depth distribution of aftershocks; (b)Temporal decay of aftershock rates

    魯?shù)榈卣鸢l(fā)生后、斷層及其周圍地區(qū)有大量的余震發(fā)生。截至2014年9月2日、共記錄到M≥0的余震活動2,694次、其中4.0~4.9級地震12次、3.0~3.9級地震62次、2.0~2.9級地震505次(來源于:中國地震臺網(wǎng)中心)。圖3a展示了余震序列分布的平面圖、可以看出、余震呈現(xiàn)共軛型分布、優(yōu)勢分布方向為NW、空間中延展約22km、穿過了該區(qū)的昭通-魯?shù)閿嗔?;近EW向延展較短、約為17km;地震主要分布于26km深度范圍內(nèi)、主要集中在4~17km。另外、還可以看到余震分布由主震處沿2個相互共軛的斷層分別向SE和近EW向擴展、且逐漸變淺。

    余震序列的衰減可以使用修正的大森公式來描述(Utsuetal.、1995):

    (9)

    式(9)中、n為后續(xù)余震的次數(shù);K、c、p是常數(shù);t表示主震后的時間。p是描述序列特征的重要參數(shù)、其值通常接近于1。利用ZMAP軟件(Wiemer、2001)擬合本次余震序列的衰減曲線(MC=2.2)、得到的結(jié)果如圖3b所示、其中p=0.70、c=0.01、K=23.8。以往的研究表明、中國大陸的p值在0.63~1.54之間變化(蔣海昆等、2006)、因此、本次魯?shù)榈卣鸬挠嗾痣m然衰減較慢、但仍然是1個正常的序列。

    2.2靜態(tài)庫侖應力變化

    本文通過計算主震產(chǎn)生的靜態(tài)庫侖應力變化、并對比余震活動在空間中的分布、考察主震對余震的觸發(fā)效果。

    圖4 魯?shù)榈卣鹩邢迶鄬幽P?張勇等、2014)Fig. 4 Dislocation model of the Ludian earthquake(after Zhang Yong et al.、2014).

    表1 區(qū)域構(gòu)造應力場的主要參數(shù)

    Table1 Main parameters of the regional tectonic stress field

    方位角/(°)傾角/(°)應力值/MPa最大主應力25410中間主應力161852最小主應力29540

    注使用魯?shù)榈卣鹫鹪礄C制解顯示的P、T軸方位(來源于:Global CMT)作為區(qū)域應力場的方位、應力值參考前人的做法。

    魯?shù)榈卣鸬恼鹪礄C制解顯示的2個節(jié)面分別為:節(jié)面Ⅰ:走向70°、傾角85°、滑動角180°;節(jié)面Ⅱ:走向160°、傾角90°、滑動角5°(來源于:Global CMT)。由前人的研究成果可知、節(jié)面Ⅱ為實際破裂面。張勇等(2014)在魯?shù)榈卣鸷?、利用地表臺站接收到的地震波形數(shù)據(jù)、反演了魯?shù)榈卣鸬钠屏堰^程模型(圖4)、使用該模型計算同震庫侖應力變化。該模型由19×12個滑移單元組成、斷層的最大深度為24km、最大位移為0.48m。物性參數(shù)方面、參照King等(1994)和Stein等(1997)的方法、巖石的泊松比和有效摩擦系數(shù)分別取0.25和0.4、剪切模量選為3.2×104MPa。計算時采用最優(yōu)破裂面進行投影、由于最優(yōu)破裂面的取向與區(qū)域應力狀態(tài)相關(guān)、應使用該地區(qū)的絕對應力場方位及傾角進行計算(Wan、2010)。但在許多情況下、某些地區(qū)的真實應力場很難獲得、因此、前人多以主震的震源機制所示的P、T軸方位近似代表區(qū)域應力場方向、選擇最優(yōu)破裂面(Kingetal.、1994; Robinsonetal.、2000; Maetal.、2005; Xuetal.、2010)。本文參考前人的做法、選擇魯?shù)榈卣鹫鹪礄C制顯示的P、T軸方位角和傾角作為區(qū)域應力場方位(表1)進行相關(guān)計算、構(gòu)造應力場的大小參照前人的選法(Kingetal.、1994; Steinetal.、1997; Todaetal.、2005)。將庫侖應力變化投影至最優(yōu)破裂面、并與1個月內(nèi)的余震空間分布進行對比。利用Coulomb程序(http: ∥www.coulombstress.org/)進行庫侖應力變化的計算。

    圖5 魯?shù)榈卣鸬撵o態(tài)庫侖應力變化在最優(yōu)破裂面的投影及余震的分布Fig. 5 Mapviews of static Coulomb stress change of Ludian earthquake and aftershock distribution.a 顯示深度為5km;b 顯示深度為10km

    圖5 顯示了魯?shù)榈卣甬a(chǎn)生的靜態(tài)庫侖應力變化在最優(yōu)破裂面上的投影、顯示的深度分別為5km和10km;黑色空心圓圈為余震活動。由圖5 可以看到、大部分的余震落入主震產(chǎn)生的庫侖應力增大區(qū)、將余震活動的分布與庫侖應力變化進行全空間范圍內(nèi)的對比、發(fā)現(xiàn)共有89.55%的余震位于庫侖應力上升的區(qū)域、而被主震所觸發(fā)(ΔCFS>0.01MPa)的余震占總數(shù)的82.43%??梢娛褂敏?shù)榈卣鸬淖顑?yōu)破裂面進行投影、得到的靜態(tài)庫侖應力變化可以解釋余震的分布。因此、本次魯?shù)榈卣饘罄m(xù)余震具有觸發(fā)作用。

    3 周圍斷層的影響

    將庫侖應力變化投影至具體斷層面、可以得到該斷層受到的加、卸載情況、判斷其未來危險性(Todaetal.、2008;Nalbantetal.、2010)。因此、本文對該區(qū)幾條主要斷層的庫侖應力變化情況進行了計算、包括昭通-魯?shù)閿嗔选⑸彿鍞嗔?、大涼山斷裂、則木河斷裂、小江斷裂和馬邊-鹽津斷裂、具體斷層參數(shù)如表2 所示(鄧起東等、2002)、物性參數(shù)與前文相同。計算中分別選擇各個斷層作為庫侖應力變化的投影面、結(jié)果如圖6 所示??梢钥吹?、距離震源較近的昭通-魯?shù)閿嗔咽艿捷^大的影響、西支東北端的庫侖應力增大較多、整體庫侖應力變化范圍為-2.330~3.905MPa;東支的東北段整體呈庫侖應力上升、庫侖應力變化為-0.917~1.147MPa、未來的地震危險性增加;對于該區(qū)的其他斷裂、由于距離相對較遠、受到的影響較小、且整體呈庫侖應力減小的趨勢、地震危險性下降。

    圖6 魯?shù)榈卣饘χ車饕顒訑鄬拥挠绊慒ig. 6 The diagram of Coulomb stress change on the surrounding major active faults.紅色區(qū)域表示該條斷裂的庫侖應力增加、危險性增加;藍色區(qū)域表示該條斷裂的庫侖應力減小、危險性降低

    表2 魯?shù)橹車饕顒訑嗔训膮?shù)與受到的庫侖應力變化

    Table2 Parameters of main active faults around epicenter and the Coulomb stress changes resolved on them

    斷裂名稱傾角/(°)滑動角/(°)ΔCFS范圍/MPa昭通-魯?shù)閿嗔?西支)57115-2.330~3.905昭通-魯?shù)閿嗔?東支)57115-0.072~0.012蓮峰斷裂79115-0.005~0.010大涼山斷裂900-0.003~0.006則木河斷裂900-0.001~0.004小江斷裂900-0.009~0.004馬邊-鹽津斷裂900-0.0006~-0.0001

    4 討論

    本文利用已發(fā)布的有限斷層模型、計算了2014年魯?shù)榈卣鸬撵o態(tài)庫侖應力變化在最優(yōu)破裂面上的投影、并對比了1個月內(nèi)的余震空間分布、可以看出魯?shù)榈卣饘罄m(xù)余震影響較大、超過80%的余震為主震所觸發(fā)。但在具體計算靜態(tài)庫侖應力變化時、有幾個問題值得我們思考:

    (1)計算1次地震所產(chǎn)生的靜態(tài)庫侖應力變化時、地震的震源模型是極為重要的。現(xiàn)有的研究中、通常使用地震波反演得到的有限斷層模型進行計算。然而、由于反演的不惟一性、得到的有限斷層模型也有多種可能、進而導致觸發(fā)效果的差異。本文使用了張勇等公布的模型進行計算、那么、對于本次魯?shù)榈卣?、使用其他的反演模型是否會得到不同的觸發(fā)效果呢?

    為了避免震源模型的獨特性所帶來的特殊觸發(fā)情況、本文使用其他學者反演的有限斷層模型進行相同的計算。劉成利等(2014)和郝金來等(2014)在本次魯?shù)榈卣鸷笸瑯永貌ㄐ螖?shù)據(jù)反演了發(fā)震斷層模型(圖7)、利用上述2個模型計算魯?shù)榈卣甬a(chǎn)生的靜態(tài)庫侖應力變化。物性參數(shù)與前面保持一致、同樣投影至最優(yōu)破裂面、得到的結(jié)果如圖8 所示。與前面的計算結(jié)果相似、利用劉成利等反演的斷層模型計算的靜態(tài)庫侖應力變化、66.86%的余震分布在庫侖應力上升超過0.01MPa的區(qū)域。而對于郝金來模型、使用同樣的參數(shù)與方法、余震觸發(fā)率為73.69%。另外、使用上面2個模型考察魯?shù)榈卣饘χ車鷶鄬拥挠绊?、如圖9 所示。對比圖6 可以看出、不同的震源模型產(chǎn)生的加、卸載作用雖然有一些細節(jié)上的變化、但是整體趨勢并沒有太大變化。值得注意的是、2個模型的計算結(jié)果中、昭通-魯?shù)閿嗔?西支)的西南段由原來的應力下降變?yōu)閼ι仙⒁虼?、該斷裂的地震危險性需要進一步關(guān)注。

    圖7 魯?shù)榈卣鹩邢迶鄬幽P虵ig. 7 Dislocation model of the Ludian earthquake.a 劉成利等(2014)提供; b 郝金來等(2014)提供

    圖8 不同震源模型計算的庫侖應力變化與余震的分布Fig. 8 Coulomb stress change calculated with different source models and aftershocks distribution.a、b 震源模型:劉成利等(2014)、顯示深度分別為5km和10km;c、d 震源模型:郝金來等(2014)、顯示深度分別為5km和10km

    圖9 不同震源模型計算的靜態(tài)庫侖應力變化對周圍斷層的影響Fig. 9 The influence of Coulomb stress change on the surrounding faults calculated by different source models.a 震源模型:劉成利等(2014); b 震源模型:郝金來等(2014)

    表3 使用不同有效摩擦系數(shù)計算的庫侖應力變化對余震的觸發(fā)效果

    Table3 Triggering effect of the Coulomb stress changes on the aftershocks calculated with different effective friction coefficient

    震源模型余震處于正庫侖應力變化區(qū)/%μ=0.0μ=0.4μ=0.8張勇等,201479.2682.4383.45劉成利等,201466.2166.8667.55郝金來等,201476.1773.6970.23

    表4 使用主震破裂面進行投影計算的庫侖應力變化對余震的觸發(fā)效果

    Table4 Triggering effect of Coulomb stress changes on the aftershocks resolved from the mainshock failure plane

    震源模型余震處于正庫侖應力變化區(qū)/%μ=0.0μ=0.4μ=0.8張勇等,20145.7724.9337.45劉成利等,201426.3930.8742.01郝金來等,201423.9130.7851.61

    另外、由式(4)、(8)可以看出、不同的有效摩擦系數(shù)對于最優(yōu)破裂面的選擇、靜態(tài)庫侖應力變化的計算都有直接的影響。因此、本文進一步研究了斷層面上的不同摩擦系數(shù)對地震觸發(fā)效果的影響。保持剪切模量和泊松比不變、分別考察有效摩擦系數(shù)變化為0.0~0.8、步長為0.1時、各個模型的庫侖應力變化與余震的分布關(guān)系。表3 給出了有效摩擦系數(shù)取0.0、0.4和0.8時庫侖應力變化對余震的觸發(fā)效果??梢钥吹?、不論斷層模型和有效摩擦系數(shù)如何變化、主震對后續(xù)余震總是具有觸發(fā)效果的。因此、對于魯?shù)榈卣?、使用主震產(chǎn)生的同震靜態(tài)庫侖應力變化、可以很好地反映本次地震產(chǎn)生的影響、解釋余震的分布。

    (2)靜態(tài)庫侖應力變化需要定義于具體斷層面、在以往的應力觸發(fā)研究中、使用各種類型的斷層作為接收斷層已有較為廣泛的研究、對于周邊斷層的危險性判斷具有較好的效果(Caskeyetal.、1997; Wanetal.、2010; 李玉江等、2013)。然而、如何選擇合適的接收斷層考察余震的觸發(fā)效果是一個較為復雜的問題。通常會使用已有的余震震源機制的節(jié)面進行投影(Steinetal.、1997)、但由于已有余震震源機制的多樣性、而且無法確定哪個節(jié)面代表真實的斷層、利用該方法考察數(shù)量龐大的余震、顯然無法得到準確的結(jié)果;另外、許多研究認為主震震源附近的余震、具有與主震相同或平行的破裂面、因此、使用與主震破裂面平行的斷層面進行投影、可以得到庫侖應力變化的最大值、用來預測未來余震的發(fā)生(Andersonetal.、1999; Chietal.、2006; Zhanetal.、2011)。

    但是、實際的計算中、使用主震破裂面進行投影得到的靜態(tài)庫侖應力變化并不能很好地觸發(fā)余震。運用本文上面提到的模型、保持物性參數(shù)不變、選擇主震破裂面作為靜態(tài)庫侖應力變化的投影面進行計算、考察對余震的觸發(fā)效果、得到的結(jié)果如表4 所示。從表中可以看到、無論震源模型和有效摩擦系數(shù)如何變化、由主震破裂面投影得到的靜態(tài)庫侖應力變化、并不能很好地解釋余震的分布。另外、對比表3 和表4 中的數(shù)據(jù)可以看到、在各種情況下、使用最優(yōu)破裂面投影的結(jié)果總是比主震破裂面投影的觸發(fā)效率高;因此、對于本次魯?shù)榈卣?、使用最?yōu)破裂面投影、計算得到的靜態(tài)庫侖應力變化可以用來解釋余震的分布、結(jié)果是比較可靠的。

    圖10 蘆山地震產(chǎn)生的靜態(tài)庫侖應力變化在最優(yōu)破裂面的投影以及余震的分布Fig. 10 Mapview of static Coulomb stress change of Lushan earthquake and aftershocks distribution.

    (3)通常情況下、以最優(yōu)破裂面投影得到的同震靜態(tài)庫侖應力變化、都可以較好地解釋余震的分布情況。那么、這樣的計算是否對于所有地震都適用呢?如2013年四川蘆山MS7.0地震、震后3個月內(nèi)共有6,798個余震發(fā)生、這些余震是否為主震所觸發(fā)呢?

    王衛(wèi)民等(2013)在蘆山地震后相繼發(fā)布了震源斷層模型、利用該模型計算靜態(tài)庫侖應力變化、選擇最優(yōu)破裂面作為投影面、其他物性參數(shù)等與上文所述相同。將計算結(jié)果與余震的空間分布相比較、考察蘆山地震對后續(xù)余震的觸發(fā)作用。然而、得到的結(jié)果并不理想、如圖10 所示??梢钥吹?、即使采用最優(yōu)破裂面進行投影、蘆山地震主震對后續(xù)余震的觸發(fā)效果也并不好、大部分余震處在靜態(tài)庫侖應力下降的區(qū)域。統(tǒng)計結(jié)果表明、僅有15.77%的余震為主震所觸發(fā)。這與魯?shù)榈卣鸬那闆r截然不同、主震對余震的觸發(fā)率明顯下降。顯然、僅使用蘆山地震的同震靜態(tài)庫侖應力變化、并不能很好地解釋余震的空間分布。因此、對于地震的靜態(tài)觸發(fā)問題、應進行更加深入的研究、包括使用不同的介質(zhì)參數(shù)計算、考慮動態(tài)庫侖應力變化和孔隙流體的作用等。

    5 結(jié)論

    通過以上數(shù)值計算與分析、得到以下初步結(jié)論:

    2014年魯?shù)镸S6.5地震對其后續(xù)余震具有觸發(fā)作用、82.43%的余震發(fā)生在靜態(tài)庫侖應力上升的區(qū)域。由于本次地震震級不大、對周圍斷層的影響較?。荷彿鍞嗔选⒋鬀錾綌嗔?、馬邊斷裂、則木河斷裂和小江斷裂整體呈卸載作用、庫侖應力變化較小;昭通-魯?shù)閿嗔?西支)的東北端受到本次地震的加載作用、昭通-魯?shù)閿嗔?東支)也呈現(xiàn)少量的庫侖應力上升、未來的地震危險性上升。另外、在考察主震對余震的觸發(fā)作用時、應使用 “最優(yōu)破裂面”進行投影;而無論震源模型與有效摩擦系數(shù)如何變化、本次魯?shù)榈卣鹗冀K對余震具有觸發(fā)作用;但同震靜態(tài)觸發(fā)對2013年蘆山地震并不適用、有待進一步的研究。

    致謝張勇、劉成利、王衛(wèi)民為本研究提供了震源模型數(shù)據(jù)、在此表示感謝。

    陳連旺、張培震、陸遠忠、等. 2008. 川滇地區(qū)強震序列庫侖破裂應力加卸載效應的數(shù)值模擬 [J]. 地球物理學報、51(5): 1411—1421.

    CHEN Lian-wang、ZHANG Pei-zhen、LU Yuan-zhong、etal. 2008. Numerical simulation of loading/unloading effect on Coulomb failure stress among strong earthquakes in Sichuan-Yunnan area [J]. Chinese Journal of Geophysics、51(5): 1411—1421(in Chinese).

    鄧起東、張培震. 2002. 中國活動構(gòu)造基本特征 [J]. 中國科學(D輯)、32(12): 1020—1030.

    DENG Qi-dong、ZHANG Pei-zhen. 2002. Basic characteristics of active tectonics of China [J]. Science in China(Ser D)、32(12): 1020—1030(in Chinese).

    郝金來、王衛(wèi)民、姚振興. 2014. 2014年8月3日云南魯?shù)镸S6.5地震震源破裂過程反演初步結(jié)果[EB/OL] http: ∥www.igg.cas.cnxwzxkyjz/201408/t20140808_4172104.html.

    HAO Jin-lai、WANG Wei-min、YAO Zhen-xing. 2014. Preliminary result for rupture process of Aug. 3、2014、MS6.5 earthquake、Ludian、Yunnan[EB/OL]. http: ∥www.igg.cas.cnxwzxkyjz/201408/t20140808_4172104.html.

    蔣海昆、曲延軍、李永莉、等. 2006. 中國大陸中強地震余震序列的部分統(tǒng)計特征 [J]. 地球物理學報、49(4): 1110—1117.

    JIANG Hai-kun、QU Yan-jun、LI Yong-li、etal. 2006. Some statistic features of aftershock sequences in Chinese mainland [J]. Chinese Journal of Geophysics、49(4): 1110—1117(in Chinese).

    李玉江、陳連旺、陸遠忠、等. 2013. 汶川地震的發(fā)生對周圍斷層穩(wěn)定性影響的數(shù)值模擬 [J]. 地球科學、38(2): 398— 410.

    LI Yu-jiang、CHEN Lian-wang、LU Yuan-zhong、etal. 2013. Numerical simulation on influences of Wenchuan earthquake on the stability of faults in the neighborhood [J]. Scientific Journal of Earth Science、38(2): 398— 410(in Chinese).

    劉成利、鄭勇、熊熊、等. 2014. 利用區(qū)域?qū)掝l帶數(shù)據(jù)反演魯?shù)镸S6.5地震震源破裂過程 [J]. 地球物理學報、57(9): 3028—3037.

    LIU Cheng-li、ZHENG Yong、XIONG Xiong、etal. 2014. Rupture process ofMS6.5 Ludian earthquake constrained by regional broadband seismograms [J]. Chinese Journal of Geophysics、57(9): 3028—3037(in Chinese).

    繆淼、朱守彪. 2012. 俯沖帶上特大地震靜態(tài)庫侖應力變化對后續(xù)余震觸發(fā)效果的研究 [J]. 地球物理學報、55,(9): 2982—2993.

    MIAO Miao、ZHU Shou-biao. 2012. A study of the impact of static Coulomb stress changes of megathrust earthquakes along subduction zone on the following aftershocks [J]. Chinese Journal of Geophysics、55(9): 2982—2993(in Chinese).

    繆淼、朱守彪. 2013. 2013年蘆山MS7.0地震產(chǎn)生的靜態(tài)庫侖應力變化及其對余震空間分布的影響 [J]. 地震學報、35(5): 619— 631.

    MIAO Miao、ZHU Shou-biao. 2013. The static Coulomb stress change of the 2013 LushanMS7.0 earthquake and its impact on the spatial distribution of aftershocks [J]. Acta Seismologica Sinica、35(5): 619— 631(in Chinese).

    單斌、李佳航、韓立波、等. 2012. 2010年MS7.1玉樹地震同震庫侖應力變化以及對2010年MS5.2囊謙地震的影響 [J]. 地球物理學報、55(9): 3028—3042.

    SHAN Bin、LI Jia-hang、HAN Li-bo、etal. 2012. Coseismic Coulomb stress change caused by 2010MS=7.1 Yushu earthquake and its influence to 2011MS=5.2 Nangqên earthquake [J]. Chinese Journal of Geophysics、55(9): 3028—3042(in Chinese).

    沈正康、萬永革、甘衛(wèi)軍、等. 2003. 東昆侖活動斷裂帶大地震之間的黏彈性應力觸發(fā)研究 [J]. 地球物理學報、46(6): 786—795.

    SHEN Zheng-kang、WAN Yong-ge、GAN Wei-jun、etal. 2003. Viscoelastic triggering among large earthquakes along the east Kunlun Fault system [J]. Chinese Journal of Geophysics、46(6): 786—795(in Chinese).

    萬永革、吳忠良、周公威、等. 2000. 幾次復雜地震中不同破裂事件之間的 “應力觸發(fā)”問題 [J]. 地震學報、22,(6): 568—576.

    WAN Yong-ge、WU Zhong-liang、ZHOU Gong-wei、etal. 2000. “Stress triggering” between different rupture events in several earthquakes [J]. Acta Seismologica Sinica、22(6): 568—576(in Chinese).

    王衛(wèi)民、郝金來、姚振興. 2013. 2013年4月20日四川蘆山地震震源破裂過程反演初步結(jié)果 [J]. 地球物理學報、56(4): 1412—1417.

    WANG Wei-min、HAO Jin-lai、YAO Zhen-xing. 2013. Preliminary result for rupture process of Apr. 20、2013、Lushan earthquake、Sichuan、China [J]. Chinese Journal of Geophysics、56(4): 1412—1417(in Chinese).

    解朝娣、朱元清、Xinglin Lei、等. 2010.MS8.0汶川地震產(chǎn)生的應力變化空間分布及其對地震活動性的影響 [J]. 中國科學(D輯)、40(6): 688— 698.

    XIE Chao-di、ZHU Yuan-qing、LEI Xing-lin、etal. 2010. Pattern of stress change and its effect on seismicity rate caused byMS8.0 Wenchuan earthquake [J]. Science in China(Ser D)、40(6): 688— 698(in Chinese).

    徐錫偉、江國焰、于貴華、等. 2014. 魯?shù)?.5級地震發(fā)震斷層判定及其構(gòu)造屬性討論 [J]. 地球物理學報、57(9): 3060—3068.

    XU Xi-wei、JIANG Guo-yan、YU Gui-hua、etal. 2014. Discussion on seismogenic fault of the LudianMS6.5 earthquake and its tectonic attribution [J]. Chinese Journal of Geophysics、57(9): 3060—3068(in Chinese).

    張廣偉、雷建設、梁姍姍、等. 2014. 2014年8月3日云南魯?shù)镸S6.5地震序列重定位與震源機制研究 [J]. 地球物理學報、57(9): 3018—3027.

    ZHANG Guang-wei、LEI Jian-she、LIANG Shan-shan、etal. 2014. Relocations and focal mechanism solutions of the 3 August 2014 Ludian、YunnanMS6.5 earthquake sequence [J]. Chinese Journal of Geophysics、57(9): 3018—3027(in Chinese).

    張勇、許力生、陳運泰、等. 2014. 2014年8月3日云南魯?shù)镸W6.1(MS6.5)地震破裂過程 [J]. 地球物理學報、57(9): 3052—3059.

    ZHANG Yong、XU Li-sheng、CHEN Yun-tai、etal. 2014. Rupture process of the 3 August 2014 Ludian、Yunnan、MW6.1(MS6.5)earthquake [J]. Chinese Journal of Geophysics、57(9): 3052—3059(in Chinese).

    張竹琪、陳永順、林間. 2008. 1997年伽師震群中相鄰正斷層和走滑斷層之間相互應力作用 [J]. 中國科學 (D輯)、38(3): 334—342.

    ZHANG Zhu-qi、CHEN Yong-shun、LIN Jian. 2008. Stress interaction between normal fault and strike-slip fault in Jiashi earthquake group 1997 [J]. Science in China(Ser D)、38(3): 334—342(in Chinese).

    Anderson G、Johnson H. 1999. A new statistical test for static stress triggering: Application to the 1987 Superstition Hills earthquake sequence [J]. Journal of Geophysical Research、104(B9): 20153—20160、20168.Caskey S J、Wesnousky S G. 1997. Static stress changes and earthquake triggering during the 1954 Fairview Peak and Dixie Valley earthquakes、central Nevada [J]. Bulletin of the Seismological Society of America、87(3): 521—527.

    Chi W C、Hauksson E. 2006. Fault-perpendicular aftershock clusters following the 2003MW=5.0 Big Bear、California、earthquake [J]. Geophysical Research Letters、33(7): L07301.

    Harris R A. 1998. Introduction to special section: Stress triggers、stress shadows、and implications for seismic hazard [J]. Journal of Geophysical Research、103(B10): 24347—24358.

    Harris R A、Simpson R W. 1996. In the shadow of 1857-the effect of the Great Ft. Tejon earthquake on subsequent earthquakes in southern California [J]. Geophysical Research Letters、23(3): 229—232.

    King G C P、Stein R S、Lin J. 1994. Static stress changes and the triggering of earthquakes [J]. Bulletin of the Seismological Society of America、84(3): 935—953.

    Lin J、Stein R S. 2004. Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults [J]. Journal of Geophysical Research、109(B2): B02303.

    Ma K、Chan C、Stein R S. 2005. Response of seismicity to Coulomb stress triggers and shadows of the 1999MW=7.6 Chi-Chi、Taiwan、earthquake [J]. Journal of Geophysical Research、110(B5): B05S19.

    Mallman E P、Zoback M D. 2007. Assessing elastic Coulomb stress transfer models using seismicity rates in southern California and southwestern Japan [J]. Journal of Geophysical Research: Solid Earth、112(B3): B03304.

    McCloskey J、Nalbant S S、Steacy S. 2005. Indonesian earthquake: Earthquake risk from co-seismic stress [J]. Nature、434: 291—291.

    Nalbant S S、McCloskey J. 2010. Stress evolution before and after the 2008 Wenchuan、China earthquake [J]. Earth and Planetary Science Letters、307(1-2): 222—232.

    Okada Y. 1992. Internal deformation due to shear and tensile faults in a half-space [J]. Bulletin of the Seismological Society of America、82(2): 1018—1040.

    Robinson R、McGinty P J. 2000. The enigma of the Arthur’s Pass、New Zealand、earthquake 2. The aftershock distribution and its relation to regional and induced stress fields [J]. Journal of Geophysical Research、105(B7): 16139—16150.

    Stein R S. 1999. The role of stress transfer in earthquake occurrence [J]. Nature、402: 605— 609.

    Stein R S、Barka A A、Dieterich J H. 1997. Progressive failure on the North Anatolian Fault since 1939 by earthquake stress triggering [J]. Geophysical Journal International、128(3): 594— 604.

    Toda S, Lin J, Meghraoui M,etal. 2008. 12 May 2008M=7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems [J]. Geophys Res Lett,35(17):L17305.

    Toda S、Stein R S. 2000. Did stress triggering cause the large off-fault aftershocks of the 25 March 1998MW=8.1 Antarctic Plate earthquake? [J]. Geophysical Research Letters、27(15): 2301—2304.

    Toda S、Stein R S、Richards-Dinger K、etal. 2005. Forecasting the evolution of seismicity in southern California: Animations built on earthquake stress transfer [J]. Journal of Geophysical Research、110(B5): B05S16.

    Utsu T、Ogata Y、Matsu’ura R S. 1995. The centenary of the Omori formula for a decay law of aftershock activity [J]. Journal of Physics of the Earth、43(3): 1—33.

    Wan Y G. 2010. Contemporary tectonic stress field in China [J]. Earthquake Science、23: 377—386.

    Wan Y Z、Shen Z K. 2010. Static Coulomb stress changes on faults caused by the 2008MW7.9 Wenchuan、China earthquake [J]. Tectonophysics、491(1-4): 105—118.Wiemer S. 2001. A software package to analyze seismicity: ZMAP [J]. Seismological Research Letters、72(3): 373—382.

    Xu C、Wang J、Li Z、etal. 2010. Applying the Coulomb failure function with an optimally oriented plane to the 2008MW7.9 Wenchuan earthquake triggering [J]. Tectonophysics、491(1-4): 119—126.

    Xu X W、Wen X Z、Zheng R Z、etal. 2003. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region、China [J]. Science in China(Ser D)、46(2): 210—226.

    Zhan Z、Jin B、Wei S、etal. 2011. Coulomb stress change sensitivity due to variability in mainshock source models and receiving fault parameters: A case study of the 2010-2011 Christchurch、New Zealand、earthquakes [J]. Seismological Research Letters、82(6): 800—814.

    Zhang P、Deng Q、Zhang G、etal. 2003. Active tectonic blocks and strong earthquakes in the continent of China [J]. Science in China (Ser D)、46(2): 13—24.

    Abstract

    On Aug. 3rd、2014、aMS6.5 earthquake struck Ludian County、Yunnan Province. It is a typical left-lateral strike-slip event. With the purpose of understanding the influence of the Ludian earthquake、this paper firstly calculates the co-seismic Coulomb failure stress changes of the mainshock with the employment of the finite dislocation source model inversed by other researchers and studies the triggering effect to the aftershocks within a month. We find that 82.43% of the aftershocks are located in the Coulomb stress increasing area(ΔCFS>0.01MPa)、therefore、most of the aftershocks are triggered by the mainshock. Then、regarding the surrounding active faults as the receive faults、the Coulomb stress changes of the mainshock are calculated to investigate the impact on the faults nearby. The result shows that only the northeast end of the west branch and northeast part of the east branch of Zhaotong-Ludian faults have been brought to failure. However、the other faults such as Daliangshan Fault、Lianfeng Fault、Zemuhe Fault、Xiaojiang Fault and Mabian-Yanjin Fault are unloaded after the Luidian event、so the possibility of future earthquake is decreased around these faults. Besides、when the optimal failure plane is chosen as the receive fault of the Coulomb stress changes、the Ludian earthquake always has good triggering effect to the aftershocks no matter which source models and effective friction coefficients are chosen.

    THE STATIC COULOMB STRESS CHANGE OF THE 2014 LUDIAN EARTHQUAKE AND ITS INFLUENCE ON THE AFTERSHOCKS AND SURROUNDING FAULTS

    MIAO Miao1、2)ZHU Shou-biao1)

    1)KeyLaboratoryofCrustalDynamics、InstituteofCrustalDynamics、ChinaEarthquakeAdministration,Beijing100085、China2)InstituteofGeophysics、ChinaEarthquakeAdministration、Beijing100081、China

    2014 Ludian earthquake、static Coulomb stress change、aftershock triggering

    10.3969/j.issn.0253-4967.2016.01.013

    2014-11-01收稿、2015-10-01改回。

    地震動力學國家重點實驗室開放基金項目(LED2012B01)與中央級科研院所科研業(yè)務專項(ZDJ2014-01、ZDJ2013-21)共同資助。

    朱守彪、男、研究員、E-mail: zhushoubiao@gmail.com。

    P315.72+7

    A

    0253-4967(2016)01-0169-13

    繆淼、男、1987年生、2012年畢業(yè)于中國地震局地殼應力研究所固體地球物理學專業(yè),獲碩士學位、助理研究員、現(xiàn)主要進行地球動力學的研究、電話010-62842659、E-mail:miaom1987@gmail.com。

    猜你喜歡
    魯?shù)?/a>主震庫侖
    1976年唐山強震群震后庫侖應力演化及其與2020年古冶5.1級地震的關(guān)系
    地震研究(2021年1期)2021-04-13 01:04:46
    基于粘彈庫侖應力變化的后續(xù)最大地震震級估計及2008、2014年于田2次7.3級地震之間關(guān)系的討論
    中國地震(2015年1期)2015-11-08 11:11:18
    多塔斜拉橋在主震-余震序列波下地震位移研究
    讓我輕輕地把魯?shù)閾肀Вü?jié)選)
    目擊
    龍卷流旋轉(zhuǎn)與地震成因
    一種周期庫侖作用勢優(yōu)化法的改進
    計算物理(2014年1期)2014-03-11 17:01:03
    長程庫侖勢對高溫超導渦旋電荷的影響
    從魯?shù)榈卣鹂赐话l(fā)事件應急管理機制的進步
    大壩與安全(2014年4期)2014-02-28 02:37:04
    2014年云南魯?shù)?.5級地震烈度分布與房屋震害特征
    地震研究(2014年4期)2014-02-27 09:31:09
    404 Not Found

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx
    404 Not Found

    404 Not Found


    nginx
    欧美bdsm另类| 久久精品熟女亚洲av麻豆精品 | 免费观看的影片在线观看| 国产极品天堂在线| 国产探花在线观看一区二区| 毛片一级片免费看久久久久| 国产极品天堂在线| 日韩成人伦理影院| 成年人午夜在线观看视频 | 久久99蜜桃精品久久| 国产精品久久久久久精品电影| 91久久精品电影网| 联通29元200g的流量卡| 丝袜美腿在线中文| 男人舔奶头视频| 97精品久久久久久久久久精品| 又大又黄又爽视频免费| 久久久精品免费免费高清| 亚洲欧美日韩卡通动漫| 内射极品少妇av片p| 国模一区二区三区四区视频| 男女边摸边吃奶| 国产精品一区二区三区四区免费观看| or卡值多少钱| 亚洲电影在线观看av| 久久精品国产亚洲av天美| 国产黄色小视频在线观看| 国产av码专区亚洲av| 秋霞伦理黄片| 能在线免费看毛片的网站| 黄片无遮挡物在线观看| kizo精华| 99热全是精品| 夫妻午夜视频| 精品酒店卫生间| 婷婷色综合www| 久久99热这里只有精品18| 亚洲av二区三区四区| 亚洲av不卡在线观看| 搡老乐熟女国产| 亚洲精品aⅴ在线观看| 久久精品熟女亚洲av麻豆精品 | 亚洲怡红院男人天堂| av专区在线播放| 人妻少妇偷人精品九色| 久久人人爽人人爽人人片va| 一个人看的www免费观看视频| 一级av片app| 一级二级三级毛片免费看| 国产黄频视频在线观看| 啦啦啦韩国在线观看视频| 一区二区三区乱码不卡18| 久久久成人免费电影| 色综合色国产| 色尼玛亚洲综合影院| 卡戴珊不雅视频在线播放| 美女主播在线视频| 久久久久久久久久久免费av| 尤物成人国产欧美一区二区三区| 久久久久久久久大av| 国产精品一区二区三区四区免费观看| 99久国产av精品| 日韩av在线大香蕉| 久久久久精品久久久久真实原创| 欧美+日韩+精品| 国产片特级美女逼逼视频| 欧美另类一区| 2021天堂中文幕一二区在线观| 精品国内亚洲2022精品成人| 亚洲国产精品成人久久小说| 精品久久久久久久久av| 午夜福利网站1000一区二区三区| 精品亚洲乱码少妇综合久久| 亚洲人与动物交配视频| av在线天堂中文字幕| 日韩伦理黄色片| 人妻夜夜爽99麻豆av| av在线亚洲专区| 嘟嘟电影网在线观看| 亚洲综合精品二区| 午夜激情久久久久久久| 亚洲精品久久久久久婷婷小说| 一级毛片aaaaaa免费看小| 亚洲怡红院男人天堂| 日日啪夜夜爽| 最近的中文字幕免费完整| 卡戴珊不雅视频在线播放| 久久99蜜桃精品久久| 成人高潮视频无遮挡免费网站| 女人被狂操c到高潮| 青青草视频在线视频观看| 美女国产视频在线观看| 国产精品综合久久久久久久免费| 91精品一卡2卡3卡4卡| 免费黄频网站在线观看国产| 国产中年淑女户外野战色| 成人午夜高清在线视频| 久久久久久久久久黄片| 国内少妇人妻偷人精品xxx网站| 国产伦精品一区二区三区四那| 免费看日本二区| 日韩欧美三级三区| 精品一区二区三区人妻视频| 亚洲人成网站在线观看播放| 久久精品熟女亚洲av麻豆精品 | 国产精品国产三级国产av玫瑰| 欧美人与善性xxx| 午夜久久久久精精品| 久久久色成人| 日本免费a在线| 一级毛片 在线播放| 91午夜精品亚洲一区二区三区| 看十八女毛片水多多多| 一级毛片aaaaaa免费看小| 女的被弄到高潮叫床怎么办| 国语对白做爰xxxⅹ性视频网站| 中国国产av一级| 搡女人真爽免费视频火全软件| 日本-黄色视频高清免费观看| 91狼人影院| av黄色大香蕉| 日本欧美国产在线视频| 淫秽高清视频在线观看| 日韩av不卡免费在线播放| 黄色日韩在线| 人妻一区二区av| 狂野欧美白嫩少妇大欣赏| 国产高清有码在线观看视频| 亚洲欧洲日产国产| 亚洲欧美精品专区久久| 成人亚洲欧美一区二区av| 不卡视频在线观看欧美| 久久精品国产亚洲av天美| 三级国产精品片| 国内精品宾馆在线| 国产成人a∨麻豆精品| 久久久久久久午夜电影| 少妇丰满av| 日韩精品青青久久久久久| 看十八女毛片水多多多| 欧美一区二区亚洲| 欧美xxⅹ黑人| 波野结衣二区三区在线| 成人欧美大片| 国产一级毛片在线| 亚洲婷婷狠狠爱综合网| 男女边吃奶边做爰视频| 免费播放大片免费观看视频在线观看| 日韩人妻高清精品专区| .国产精品久久| 麻豆成人av视频| 中文字幕人妻熟人妻熟丝袜美| 一本一本综合久久| 校园人妻丝袜中文字幕| 七月丁香在线播放| 亚洲精品亚洲一区二区| 夫妻性生交免费视频一级片| 欧美人与善性xxx| 久久这里有精品视频免费| 极品教师在线视频| 老司机影院成人| 国产综合精华液| 国产精品日韩av在线免费观看| 日韩伦理黄色片| av播播在线观看一区| 直男gayav资源| 能在线免费观看的黄片| 亚洲va在线va天堂va国产| 久久97久久精品| a级毛片免费高清观看在线播放| 一级片'在线观看视频| 国产精品av视频在线免费观看| 国产乱人视频| av.在线天堂| 欧美最新免费一区二区三区| 亚洲欧美日韩东京热| 日韩成人伦理影院| 免费av观看视频| 赤兔流量卡办理| 80岁老熟妇乱子伦牲交| 观看免费一级毛片| 又爽又黄无遮挡网站| 精品人妻熟女av久视频| 亚洲国产精品成人久久小说| 婷婷色综合www| 亚洲精品国产av成人精品| 91精品伊人久久大香线蕉| 一夜夜www| 超碰av人人做人人爽久久| 久久人人爽人人爽人人片va| 国产成人91sexporn| 神马国产精品三级电影在线观看| 汤姆久久久久久久影院中文字幕 | 日本三级黄在线观看| 国产精品嫩草影院av在线观看| 最后的刺客免费高清国语| 好男人视频免费观看在线| 天天躁日日操中文字幕| 国产精品久久久久久av不卡| 国产免费视频播放在线视频 | 亚洲成人一二三区av| 日韩一区二区三区影片| 国产精品美女特级片免费视频播放器| 一个人免费在线观看电影| 一级毛片电影观看| 少妇的逼好多水| 丰满人妻一区二区三区视频av| 日韩不卡一区二区三区视频在线| 男人舔奶头视频| 日韩av在线免费看完整版不卡| 欧美三级亚洲精品| 成人亚洲欧美一区二区av| 久久鲁丝午夜福利片| 国产伦理片在线播放av一区| 91精品一卡2卡3卡4卡| 国产欧美日韩精品一区二区| videos熟女内射| 秋霞在线观看毛片| 青春草国产在线视频| 国产片特级美女逼逼视频| 日本欧美国产在线视频| 在线 av 中文字幕| 男女视频在线观看网站免费| 91在线精品国自产拍蜜月| 国产精品三级大全| 亚洲自偷自拍三级| 久久久久久久久久黄片| 丰满人妻一区二区三区视频av| av免费在线看不卡| 亚洲自拍偷在线| 亚洲精品久久久久久婷婷小说| 久久精品人妻少妇| 亚洲激情五月婷婷啪啪| 人人妻人人澡人人爽人人夜夜 | 亚洲精品一区蜜桃| 能在线免费观看的黄片| 久久久久九九精品影院| 99热网站在线观看| 欧美日韩在线观看h| 日产精品乱码卡一卡2卡三| 久久久午夜欧美精品| av在线观看视频网站免费| 日本与韩国留学比较| 熟女电影av网| 日本熟妇午夜| 99久久九九国产精品国产免费| 天天一区二区日本电影三级| 少妇高潮的动态图| 少妇熟女aⅴ在线视频| 老司机影院毛片| 欧美精品一区二区大全| 欧美精品一区二区大全| 最近最新中文字幕免费大全7| 床上黄色一级片| 好男人视频免费观看在线| 99久久精品国产国产毛片| 乱系列少妇在线播放| 99热这里只有是精品在线观看| 欧美成人精品欧美一级黄| 日本熟妇午夜| 国产探花极品一区二区| 三级经典国产精品| 你懂的网址亚洲精品在线观看| av.在线天堂| 国产男女超爽视频在线观看| 22中文网久久字幕| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧洲日产国产| 美女内射精品一级片tv| 午夜精品在线福利| 日韩一区二区三区影片| 亚洲成人一二三区av| 观看美女的网站| 听说在线观看完整版免费高清| av卡一久久| 男插女下体视频免费在线播放| av女优亚洲男人天堂| 嘟嘟电影网在线观看| 看黄色毛片网站| 九草在线视频观看| 老司机影院毛片| 中文天堂在线官网| 婷婷色麻豆天堂久久| 国产午夜精品久久久久久一区二区三区| 婷婷色综合大香蕉| 国产成人福利小说| 国产黄色小视频在线观看| 天天躁日日操中文字幕| 午夜福利视频1000在线观看| 亚洲成人一二三区av| 午夜福利成人在线免费观看| 嘟嘟电影网在线观看| 熟妇人妻不卡中文字幕| 十八禁国产超污无遮挡网站| 啦啦啦韩国在线观看视频| kizo精华| 国产毛片a区久久久久| 伦理电影大哥的女人| 99re6热这里在线精品视频| 最近中文字幕2019免费版| 一区二区三区高清视频在线| 亚洲精品乱码久久久久久按摩| 免费观看a级毛片全部| av免费在线看不卡| 麻豆国产97在线/欧美| 美女xxoo啪啪120秒动态图| 日本免费a在线| 亚洲婷婷狠狠爱综合网| 日韩 亚洲 欧美在线| 美女高潮的动态| 亚洲精品一区蜜桃| 国产久久久一区二区三区| a级一级毛片免费在线观看| 国产69精品久久久久777片| av黄色大香蕉| 国产精品不卡视频一区二区| 久久久久久久久久久免费av| 日韩人妻高清精品专区| 男人舔女人下体高潮全视频| 亚洲av成人av| 亚洲欧美日韩卡通动漫| 精品人妻熟女av久视频| 九九爱精品视频在线观看| 国产免费又黄又爽又色| av在线亚洲专区| 肉色欧美久久久久久久蜜桃 | 一二三四中文在线观看免费高清| 深爱激情五月婷婷| 亚洲国产最新在线播放| 2021少妇久久久久久久久久久| 麻豆精品久久久久久蜜桃| 亚洲av电影不卡..在线观看| 国产在线一区二区三区精| 欧美变态另类bdsm刘玥| 亚洲欧美日韩东京热| 永久免费av网站大全| 日韩精品有码人妻一区| a级毛片免费高清观看在线播放| 美女内射精品一级片tv| videossex国产| 亚洲精品成人av观看孕妇| 午夜福利成人在线免费观看| 成人午夜精彩视频在线观看| 老司机影院毛片| 国产成人精品久久久久久| 日韩欧美精品v在线| 亚洲熟妇中文字幕五十中出| 嫩草影院精品99| 欧美xxxx黑人xx丫x性爽| 亚洲欧美中文字幕日韩二区| 又爽又黄无遮挡网站| 久久精品国产亚洲av天美| 18+在线观看网站| 夫妻性生交免费视频一级片| 国内少妇人妻偷人精品xxx网站| 深夜a级毛片| 亚洲内射少妇av| 国语对白做爰xxxⅹ性视频网站| 国产精品蜜桃在线观看| 中国美白少妇内射xxxbb| 欧美三级亚洲精品| 国产 一区精品| 国产精品一及| 啦啦啦中文免费视频观看日本| 久久久久精品性色| 精品人妻一区二区三区麻豆| 春色校园在线视频观看| 男女视频在线观看网站免费| 女人被狂操c到高潮| 国产成人午夜福利电影在线观看| 午夜福利在线观看免费完整高清在| 少妇被粗大猛烈的视频| 亚洲av一区综合| 国产精品一区www在线观看| 精品久久久久久久久久久久久| 伊人久久精品亚洲午夜| 国产免费一级a男人的天堂| 九九久久精品国产亚洲av麻豆| 免费大片18禁| 亚洲欧美一区二区三区国产| 国产日韩欧美在线精品| 看十八女毛片水多多多| 美女大奶头视频| 狂野欧美激情性xxxx在线观看| 日韩欧美精品免费久久| 99热这里只有精品一区| 精品久久久久久久久久久久久| 国产探花极品一区二区| 国产精品国产三级国产av玫瑰| 午夜激情福利司机影院| 少妇高潮的动态图| av免费在线看不卡| 国产精品一区www在线观看| 国产成人freesex在线| 在线免费观看的www视频| 国产中年淑女户外野战色| 岛国毛片在线播放| 精品久久久久久电影网| 人妻少妇偷人精品九色| 国产精品一区二区三区四区免费观看| 26uuu在线亚洲综合色| 午夜精品在线福利| 亚洲,欧美,日韩| 午夜日本视频在线| 精品久久久久久久久av| 国产精品嫩草影院av在线观看| 亚洲国产成人一精品久久久| 日韩成人伦理影院| 日产精品乱码卡一卡2卡三| 久久久国产一区二区| 国产精品熟女久久久久浪| 精品久久久久久久久久久久久| 黄色配什么色好看| 国产一区二区在线观看日韩| 免费看美女性在线毛片视频| 久久精品国产亚洲网站| 亚洲三级黄色毛片| 人人妻人人澡人人爽人人夜夜 | 国产精品.久久久| 婷婷色麻豆天堂久久| 国产一区二区亚洲精品在线观看| 美女cb高潮喷水在线观看| 中文字幕免费在线视频6| 久久久久免费精品人妻一区二区| 男女边摸边吃奶| 又粗又硬又长又爽又黄的视频| 99久久精品国产国产毛片| 又爽又黄a免费视频| 99视频精品全部免费 在线| 亚洲欧洲国产日韩| av在线播放精品| 欧美成人精品欧美一级黄| 欧美不卡视频在线免费观看| 深夜a级毛片| 乱人视频在线观看| 久久久久性生活片| 精品国内亚洲2022精品成人| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产精品成人久久小说| 亚洲人成网站在线播| 岛国毛片在线播放| 日韩伦理黄色片| 天美传媒精品一区二区| 免费看日本二区| 一个人看视频在线观看www免费| 国产精品国产三级专区第一集| 久久久国产一区二区| 婷婷色综合www| 国产色婷婷99| 日韩,欧美,国产一区二区三区| 国产伦一二天堂av在线观看| 久久久久精品性色| 日韩成人av中文字幕在线观看| 神马国产精品三级电影在线观看| 国产av码专区亚洲av| 午夜精品在线福利| 久久久久久久国产电影| 精品少妇黑人巨大在线播放| 少妇人妻一区二区三区视频| 视频中文字幕在线观看| 免费看日本二区| 伦精品一区二区三区| av在线观看视频网站免费| 乱人视频在线观看| 色视频www国产| 国产伦在线观看视频一区| 亚洲精品视频女| 亚洲最大成人av| 日韩一区二区三区影片| 久久久久免费精品人妻一区二区| 国产亚洲5aaaaa淫片| 最近视频中文字幕2019在线8| 亚洲国产色片| 在线播放无遮挡| 深夜a级毛片| 午夜福利成人在线免费观看| 国产免费又黄又爽又色| 午夜福利在线观看免费完整高清在| 亚州av有码| 好男人在线观看高清免费视频| 国产免费又黄又爽又色| 欧美精品一区二区大全| 看非洲黑人一级黄片| 欧美日韩国产mv在线观看视频 | 亚洲一级一片aⅴ在线观看| 日本一二三区视频观看| 真实男女啪啪啪动态图| 成人毛片60女人毛片免费| 日本猛色少妇xxxxx猛交久久| 51国产日韩欧美| 午夜视频国产福利| 免费看a级黄色片| 最近中文字幕高清免费大全6| 九九爱精品视频在线观看| av天堂中文字幕网| 日韩大片免费观看网站| 99久国产av精品国产电影| 欧美97在线视频| 日韩人妻高清精品专区| 国产成人福利小说| 亚洲精品第二区| 内射极品少妇av片p| 五月伊人婷婷丁香| 成人无遮挡网站| 国产淫语在线视频| 婷婷色麻豆天堂久久| 美女xxoo啪啪120秒动态图| 麻豆乱淫一区二区| 日韩成人伦理影院| av国产免费在线观看| 午夜福利视频精品| 国产日韩欧美在线精品| 国产av在哪里看| 日韩av在线大香蕉| 亚洲美女搞黄在线观看| 精品国产一区二区三区久久久樱花 | 国产精品一区二区三区四区久久| 99热全是精品| 欧美不卡视频在线免费观看| 青春草亚洲视频在线观看| 久99久视频精品免费| 赤兔流量卡办理| 久久久久久久久久黄片| 乱码一卡2卡4卡精品| 久久久久精品久久久久真实原创| 又爽又黄无遮挡网站| 亚洲自拍偷在线| 99热这里只有是精品在线观看| 在线观看av片永久免费下载| 老司机影院成人| 国内精品美女久久久久久| 国产亚洲av嫩草精品影院| 在线播放无遮挡| 亚洲欧洲国产日韩| 成人亚洲精品av一区二区| 91午夜精品亚洲一区二区三区| 欧美性感艳星| 国产午夜精品一二区理论片| 有码 亚洲区| 久久精品熟女亚洲av麻豆精品 | 激情 狠狠 欧美| 在线播放无遮挡| 成人漫画全彩无遮挡| 中文字幕av在线有码专区| 国产亚洲最大av| 中文字幕av在线有码专区| 国产亚洲最大av| 国产亚洲精品av在线| 亚洲av一区综合| 搞女人的毛片| 欧美日韩视频高清一区二区三区二| 色5月婷婷丁香| 国产精品爽爽va在线观看网站| 午夜视频国产福利| 国产黄a三级三级三级人| 啦啦啦韩国在线观看视频| 国产一区亚洲一区在线观看| 又黄又爽又刺激的免费视频.| 国产 亚洲一区二区三区 | 极品教师在线视频| av免费观看日本| 如何舔出高潮| 亚洲欧洲国产日韩| 直男gayav资源| 久久99热这里只有精品18| 日本三级黄在线观看| 一级毛片 在线播放| 又黄又爽又刺激的免费视频.| 亚洲精品乱久久久久久| 真实男女啪啪啪动态图| 国产精品人妻久久久影院| 免费黄频网站在线观看国产| 亚洲自偷自拍三级| 久久精品熟女亚洲av麻豆精品 | 亚洲av免费高清在线观看| 身体一侧抽搐| 亚洲国产色片| 亚洲精品亚洲一区二区| 国产亚洲一区二区精品| 少妇人妻一区二区三区视频| 特大巨黑吊av在线直播| 又大又黄又爽视频免费| 日日啪夜夜爽| 国产男人的电影天堂91| 午夜视频国产福利| 舔av片在线| 亚洲性久久影院| 少妇猛男粗大的猛烈进出视频 | 国产老妇女一区| 国产亚洲精品久久久com| 国产女主播在线喷水免费视频网站 | 蜜桃久久精品国产亚洲av| 国产成人91sexporn| 日本一本二区三区精品| 国产男人的电影天堂91| 美女被艹到高潮喷水动态| 久久久精品94久久精品| 国产成人一区二区在线| 色播亚洲综合网| 熟女电影av网| 精品久久久精品久久久| 精品一区二区三区人妻视频| 国产美女午夜福利| 日韩欧美精品v在线| 欧美高清性xxxxhd video| 久久久成人免费电影| 国产精品麻豆人妻色哟哟久久 | av在线老鸭窝| 久久久久久久午夜电影| 久久久久久久久久成人| 亚洲婷婷狠狠爱综合网| 大话2 男鬼变身卡| 99热6这里只有精品| 婷婷色综合大香蕉| 天堂影院成人在线观看|