• <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
    国产女主播在线喷水免费视频网站 | 天天一区二区日本电影三级| 丰满乱子伦码专区| 男人舔女人下体高潮全视频| 自拍偷自拍亚洲精品老妇| 欧美+亚洲+日韩+国产| 全区人妻精品视频| 成人高潮视频无遮挡免费网站| 欧美日韩国产亚洲二区| 国产精品不卡视频一区二区| 日韩大尺度精品在线看网址| 日韩视频在线欧美| 可以在线观看毛片的网站| 直男gayav资源| 国产大屁股一区二区在线视频| 内地一区二区视频在线| 精品久久久噜噜| 一级毛片aaaaaa免费看小| 真实男女啪啪啪动态图| 99久久成人亚洲精品观看| 尤物成人国产欧美一区二区三区| 亚洲美女搞黄在线观看| 尤物成人国产欧美一区二区三区| 十八禁国产超污无遮挡网站| 国模一区二区三区四区视频| 菩萨蛮人人尽说江南好唐韦庄 | 成人高潮视频无遮挡免费网站| 免费观看的影片在线观看| 欧美日韩在线观看h| 22中文网久久字幕| 日本免费一区二区三区高清不卡| 国产午夜福利久久久久久| av福利片在线观看| h日本视频在线播放| 老司机福利观看| 免费观看的影片在线观看| 日韩 亚洲 欧美在线| 亚洲高清免费不卡视频| 一区二区三区免费毛片| 性色avwww在线观看| 一级毛片我不卡| 国产老妇伦熟女老妇高清| 免费在线观看成人毛片| 亚洲精品456在线播放app| 女人被狂操c到高潮| 久久99热6这里只有精品| 欧美丝袜亚洲另类| 久久久国产成人精品二区| 亚洲欧美精品综合久久99| 亚洲不卡免费看| 国产探花在线观看一区二区| 精品一区二区三区视频在线| 亚洲国产精品成人久久小说 | 少妇高潮的动态图| 亚洲经典国产精华液单| 美女脱内裤让男人舔精品视频 | 男女做爰动态图高潮gif福利片| 人人妻人人澡欧美一区二区| 国内精品一区二区在线观看| 熟妇人妻久久中文字幕3abv| 国产伦在线观看视频一区| 国产精品一区二区三区四区免费观看| 看免费成人av毛片| 久久久久久久久久成人| 成人高潮视频无遮挡免费网站| 亚洲中文字幕日韩| 久久久久久九九精品二区国产| 在线观看av片永久免费下载| 亚洲欧美日韩无卡精品| 菩萨蛮人人尽说江南好唐韦庄 | 高清午夜精品一区二区三区 | 最近的中文字幕免费完整| 麻豆精品久久久久久蜜桃| 国产亚洲av片在线观看秒播厂 | 国产午夜精品久久久久久一区二区三区| 亚洲精品乱码久久久久久按摩| 蜜臀久久99精品久久宅男| 91aial.com中文字幕在线观看| 国产精品福利在线免费观看| 久久99热这里只有精品18| 日韩欧美三级三区| 亚洲18禁久久av| 欧美性感艳星| 美女内射精品一级片tv| 亚洲va在线va天堂va国产| 国产精品美女特级片免费视频播放器| 久久99热6这里只有精品| 91精品国产九色| 精品人妻一区二区三区麻豆| 国产又黄又爽又无遮挡在线| 国产 一区 欧美 日韩| 国产精品电影一区二区三区| 久久亚洲国产成人精品v| 国产高清不卡午夜福利| 91久久精品电影网| 午夜a级毛片| 国产男人的电影天堂91| 久久精品综合一区二区三区| 日日摸夜夜添夜夜爱| 春色校园在线视频观看| 久久精品国产亚洲网站| 亚洲图色成人| 国产探花在线观看一区二区| 国产探花极品一区二区| 日本黄色视频三级网站网址| 此物有八面人人有两片| 久久久成人免费电影| 乱人视频在线观看| 校园人妻丝袜中文字幕| 免费人成视频x8x8入口观看| 亚洲最大成人手机在线| 极品教师在线视频| 大又大粗又爽又黄少妇毛片口| 尾随美女入室| 日本黄大片高清| 少妇人妻一区二区三区视频| АⅤ资源中文在线天堂| 亚洲欧美日韩卡通动漫| 99热6这里只有精品| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩东京热| 亚洲一区二区三区色噜噜| 3wmmmm亚洲av在线观看| 亚洲欧美成人综合另类久久久 | 亚洲精品色激情综合| 亚洲国产精品成人综合色| 亚洲成人久久爱视频| 亚洲第一电影网av| 少妇的逼水好多| 精品熟女少妇av免费看| 欧美日本亚洲视频在线播放| 我的女老师完整版在线观看| 欧美性感艳星| 熟女人妻精品中文字幕| 久久久国产成人精品二区| 一本精品99久久精品77| 国产精品不卡视频一区二区| 国产女主播在线喷水免费视频网站 | 国产黄a三级三级三级人| av在线播放精品| 久久久久网色| 特大巨黑吊av在线直播| 中文字幕av在线有码专区| 在线天堂最新版资源| av国产免费在线观看| 国产精品精品国产色婷婷| 午夜老司机福利剧场| 国产久久久一区二区三区| 国产精品99久久久久久久久| 欧美精品一区二区大全| 午夜a级毛片| 少妇熟女aⅴ在线视频| 日日摸夜夜添夜夜添av毛片| 永久网站在线| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲真实伦在线观看| 国产精品,欧美在线| 搡老妇女老女人老熟妇| 国产精品美女特级片免费视频播放器| 久久精品国产鲁丝片午夜精品| 尤物成人国产欧美一区二区三区| 精品日产1卡2卡| 看黄色毛片网站| 国产一级毛片七仙女欲春2| 成人三级黄色视频| 狂野欧美激情性xxxx在线观看| 插阴视频在线观看视频| 久久精品国产99精品国产亚洲性色| 亚洲av中文字字幕乱码综合| 成人亚洲欧美一区二区av| av在线观看视频网站免费| 亚洲av成人av| 亚洲高清免费不卡视频| 欧美3d第一页| 啦啦啦观看免费观看视频高清| 精品不卡国产一区二区三区| 91av网一区二区| 深爱激情五月婷婷| 国产人妻一区二区三区在| 国产成人a区在线观看| 国产黄片美女视频| av在线观看视频网站免费| 国产成人精品婷婷| 欧美成人一区二区免费高清观看| 亚洲欧美精品综合久久99| 日韩,欧美,国产一区二区三区 | 最好的美女福利视频网| 九色成人免费人妻av| 禁无遮挡网站| 久久精品国产亚洲av天美| 白带黄色成豆腐渣| 男人舔奶头视频| 一边亲一边摸免费视频| 日韩高清综合在线| 校园春色视频在线观看| 国产探花极品一区二区| 2022亚洲国产成人精品| 精品久久国产蜜桃| 99久国产av精品国产电影| 国产精品一二三区在线看| 成人午夜精彩视频在线观看| 国产免费一级a男人的天堂| 97超视频在线观看视频| 91久久精品国产一区二区三区| 免费搜索国产男女视频| 国产一级毛片在线| 国产午夜精品一二区理论片| 国内少妇人妻偷人精品xxx网站| 亚州av有码| 精品欧美国产一区二区三| av在线老鸭窝| 国产白丝娇喘喷水9色精品| 日韩精品有码人妻一区| 美女xxoo啪啪120秒动态图| 精品人妻熟女av久视频| 国产探花在线观看一区二区| 人体艺术视频欧美日本| 国产伦精品一区二区三区视频9| 色播亚洲综合网| 午夜精品国产一区二区电影 | 女的被弄到高潮叫床怎么办| 十八禁国产超污无遮挡网站| 婷婷精品国产亚洲av| 日韩中字成人| or卡值多少钱| 少妇熟女欧美另类| 人人妻人人看人人澡| 插逼视频在线观看| 国产成人91sexporn| 免费人成视频x8x8入口观看| 国产精品久久久久久久久免| 欧美日韩国产亚洲二区| 日本五十路高清| 亚洲最大成人av| 婷婷色av中文字幕| 国产午夜精品一二区理论片| 亚洲人成网站在线播放欧美日韩| 亚洲综合色惰| 99riav亚洲国产免费| 午夜福利成人在线免费观看| 久久久欧美国产精品| 12—13女人毛片做爰片一| 天堂√8在线中文| 能在线免费观看的黄片| 欧美变态另类bdsm刘玥| 国产三级在线视频| 狠狠狠狠99中文字幕| 亚洲美女搞黄在线观看| 久久久久国产网址| av国产免费在线观看| 六月丁香七月| 久久久精品大字幕| 两个人的视频大全免费| 国产成人精品婷婷| 国产美女午夜福利| 欧美日本视频| 日韩人妻高清精品专区| 欧美性猛交╳xxx乱大交人| 一级黄片播放器| 大又大粗又爽又黄少妇毛片口| 又粗又爽又猛毛片免费看| 国产成年人精品一区二区| 国产亚洲av嫩草精品影院| 国产探花在线观看一区二区| 日韩人妻高清精品专区| 18禁在线播放成人免费| 一级av片app| 亚洲国产精品合色在线| www.色视频.com| 给我免费播放毛片高清在线观看| 别揉我奶头 嗯啊视频| 在线天堂最新版资源| 亚洲第一区二区三区不卡| 国内精品一区二区在线观看| 我的老师免费观看完整版| 亚洲内射少妇av| 免费在线观看成人毛片| 天天一区二区日本电影三级| 国产不卡一卡二| 国产成人影院久久av| 国产色婷婷99| 久久热精品热| 色综合站精品国产| 啦啦啦啦在线视频资源| 少妇的逼水好多| 午夜免费男女啪啪视频观看| 自拍偷自拍亚洲精品老妇| 欧美日韩乱码在线| 欧美另类亚洲清纯唯美| 国产精品久久久久久久电影| 亚洲熟妇中文字幕五十中出| 国产精品久久久久久精品电影小说 | 国产 一区精品| 亚洲七黄色美女视频| 99久国产av精品| 网址你懂的国产日韩在线| 国产成人91sexporn| 亚洲最大成人手机在线| 久久久久免费精品人妻一区二区| 国产不卡一卡二| 能在线免费看毛片的网站| 啦啦啦观看免费观看视频高清| 此物有八面人人有两片| 日韩成人av中文字幕在线观看| 在线免费十八禁| 欧美日韩国产亚洲二区| 久久婷婷人人爽人人干人人爱| 午夜精品在线福利| 欧美日本视频| 久久婷婷人人爽人人干人人爱| 亚洲av免费高清在线观看| 国产亚洲精品av在线| 26uuu在线亚洲综合色| 男人的好看免费观看在线视频| 欧美最黄视频在线播放免费| 日产精品乱码卡一卡2卡三| 日韩成人av中文字幕在线观看| 国产久久久一区二区三区| 日本黄大片高清| 欧美性猛交黑人性爽| 午夜精品国产一区二区电影 | 丝袜喷水一区| 久久精品91蜜桃| 亚洲精品国产成人久久av| 久久精品久久久久久久性| 亚洲av二区三区四区| 久久久久久久久久黄片| 少妇熟女欧美另类| av福利片在线观看| 男女边吃奶边做爰视频| 国产91av在线免费观看| 麻豆久久精品国产亚洲av| 有码 亚洲区| 女人被狂操c到高潮| 三级男女做爰猛烈吃奶摸视频| 亚洲国产欧洲综合997久久,| 国产亚洲欧美98| 一个人看视频在线观看www免费| 亚洲成人精品中文字幕电影| 99国产精品一区二区蜜桃av| 亚洲图色成人| 男的添女的下面高潮视频| 一夜夜www| 亚洲天堂国产精品一区在线| 欧美三级亚洲精品| 我要看日韩黄色一级片| 中文字幕免费在线视频6| 一区二区三区高清视频在线| 校园人妻丝袜中文字幕| 观看免费一级毛片| 国产亚洲5aaaaa淫片| а√天堂www在线а√下载| 中文字幕免费在线视频6| 成年女人看的毛片在线观看| 一级毛片电影观看 | 欧美一区二区亚洲| 日本在线视频免费播放| 久久久久久九九精品二区国产| 青青草视频在线视频观看| 午夜久久久久精精品| 国产爱豆传媒在线观看| 免费人成视频x8x8入口观看| 青春草视频在线免费观看| 日产精品乱码卡一卡2卡三| 一级毛片久久久久久久久女| 直男gayav资源| 日本免费一区二区三区高清不卡| 国产黄色小视频在线观看| 老师上课跳d突然被开到最大视频| 久久久久久九九精品二区国产| 99国产精品一区二区蜜桃av| 成人美女网站在线观看视频| 人人妻人人澡欧美一区二区| 一级毛片我不卡| 国产免费男女视频| 亚洲欧洲日产国产| 九九爱精品视频在线观看| 联通29元200g的流量卡| 少妇丰满av| 男人和女人高潮做爰伦理| 亚洲,欧美,日韩| 91麻豆精品激情在线观看国产| 麻豆精品久久久久久蜜桃| 寂寞人妻少妇视频99o| 久久6这里有精品| 久久久久久久久中文| 亚洲精品色激情综合| 国产日本99.免费观看| 亚洲中文字幕一区二区三区有码在线看| 欧美性感艳星| av在线老鸭窝| 国产三级在线视频| 给我免费播放毛片高清在线观看| 内射极品少妇av片p| 国产精品一区二区三区四区久久| 国内精品宾馆在线| 日韩国内少妇激情av| 久久亚洲精品不卡| 韩国av在线不卡| 色哟哟哟哟哟哟| 欧美成人a在线观看| 亚洲第一电影网av| 国产三级在线视频| 国产亚洲av嫩草精品影院| 内射极品少妇av片p| 国产成人精品婷婷| 一卡2卡三卡四卡精品乱码亚洲| 亚洲美女搞黄在线观看| 国产精品.久久久| 国产成人精品婷婷| 日本在线视频免费播放| 2022亚洲国产成人精品| 能在线免费看毛片的网站| 亚洲性久久影院| 亚洲综合色惰| 又爽又黄无遮挡网站| av在线蜜桃| 亚洲高清免费不卡视频| 国产午夜精品论理片| 欧美精品国产亚洲| 免费黄网站久久成人精品| 亚洲真实伦在线观看| 好男人视频免费观看在线| 国产av一区在线观看免费| 亚洲欧美日韩无卡精品| 亚洲人与动物交配视频| 熟女电影av网| 国产黄片视频在线免费观看| 99久久精品一区二区三区| 男女下面进入的视频免费午夜| 人妻少妇偷人精品九色| 人体艺术视频欧美日本| 91久久精品国产一区二区三区| 男女那种视频在线观看| 热99在线观看视频| 久久精品久久久久久久性| 成人毛片a级毛片在线播放| 午夜福利在线观看吧| 国产伦在线观看视频一区| 欧美一区二区国产精品久久精品| 尤物成人国产欧美一区二区三区| 久久久久网色| 午夜精品国产一区二区电影 | 日本-黄色视频高清免费观看| 国产一区二区激情短视频| 国产乱人偷精品视频| 亚洲精品成人久久久久久| 韩国av在线不卡| 一进一出抽搐gif免费好疼| 国产黄片视频在线免费观看| 久久久欧美国产精品| 又粗又硬又长又爽又黄的视频 | 成人亚洲欧美一区二区av| av在线播放精品| 亚洲欧美成人精品一区二区| 中文精品一卡2卡3卡4更新| 国产精品久久久久久久久免| 日本五十路高清| 亚洲国产精品成人久久小说 | 成人无遮挡网站| 99久久九九国产精品国产免费| 日本av手机在线免费观看| 97超碰精品成人国产| 床上黄色一级片| 最近最新中文字幕大全电影3| 久久精品国产亚洲av香蕉五月| 欧美成人精品欧美一级黄| 国产三级在线视频| 最近的中文字幕免费完整| 悠悠久久av| 亚洲av免费高清在线观看| 高清午夜精品一区二区三区 | 亚洲色图av天堂| 日本爱情动作片www.在线观看| 久久亚洲国产成人精品v| 色综合站精品国产| 啦啦啦观看免费观看视频高清| 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区| 亚洲国产精品sss在线观看| 人妻系列 视频| 男女啪啪激烈高潮av片| 国产精品嫩草影院av在线观看| 九九在线视频观看精品| 亚洲五月天丁香| 成年版毛片免费区| 国产黄片视频在线免费观看| 在线播放无遮挡| 国产精品爽爽va在线观看网站| 国产伦精品一区二区三区四那| 亚洲18禁久久av| 日本一本二区三区精品| 国产精品,欧美在线| 天堂影院成人在线观看| 久久久久久久久大av| 蜜桃久久精品国产亚洲av| 精品久久久久久久久亚洲| 免费看日本二区| 久久欧美精品欧美久久欧美| 国产成人aa在线观看| 成人毛片60女人毛片免费| 亚洲精品国产av成人精品| 99久久九九国产精品国产免费| 久久精品影院6| 亚洲欧美日韩高清在线视频| av又黄又爽大尺度在线免费看 | 成人av在线播放网站| 永久网站在线| 男插女下体视频免费在线播放| 欧美日韩综合久久久久久| 一边摸一边抽搐一进一小说| 欧美最新免费一区二区三区| av卡一久久| av福利片在线观看| 国产精品久久视频播放| 三级经典国产精品| 狂野欧美激情性xxxx在线观看| 国产成人午夜福利电影在线观看| 国国产精品蜜臀av免费| 精品久久久久久久人妻蜜臀av| 哪里可以看免费的av片| 成人一区二区视频在线观看| 天堂√8在线中文| 成人欧美大片| 国产精品三级大全| 波多野结衣巨乳人妻| 色综合站精品国产| 人人妻人人看人人澡| 中文字幕熟女人妻在线| kizo精华| 听说在线观看完整版免费高清| 国产伦理片在线播放av一区 | 青春草亚洲视频在线观看| 久久精品国产99精品国产亚洲性色| 亚洲欧美精品专区久久| 中出人妻视频一区二区| 亚洲精品成人久久久久久| 欧美成人免费av一区二区三区| 97人妻精品一区二区三区麻豆| 中出人妻视频一区二区| 日本色播在线视频| 晚上一个人看的免费电影| 国产不卡一卡二| 高清毛片免费看| 免费av毛片视频| 国产高清三级在线| 久久人人爽人人爽人人片va| 国产精品电影一区二区三区| 男女边吃奶边做爰视频| 日韩欧美在线乱码| 少妇裸体淫交视频免费看高清| 日韩,欧美,国产一区二区三区 | 男女视频在线观看网站免费| 黄色日韩在线| 国产成人a∨麻豆精品| h日本视频在线播放| 免费看美女性在线毛片视频| 人妻制服诱惑在线中文字幕| 免费人成视频x8x8入口观看| 人体艺术视频欧美日本| 在现免费观看毛片| 午夜福利高清视频| 成年免费大片在线观看| 哪个播放器可以免费观看大片| 国产91av在线免费观看| 精品久久久噜噜| 热99re8久久精品国产| 国产乱人视频| 亚洲成a人片在线一区二区| 成人三级黄色视频| 内地一区二区视频在线| 久久久久久久久大av| 青春草视频在线免费观看| 国产黄a三级三级三级人| 国产精品不卡视频一区二区| 亚洲av第一区精品v没综合| 三级经典国产精品| 久久午夜福利片| 99久久久亚洲精品蜜臀av| 午夜激情福利司机影院| av在线观看视频网站免费| 一级av片app| 久久精品国产亚洲av香蕉五月| 又黄又爽又刺激的免费视频.| 国产高清视频在线观看网站| 国产伦精品一区二区三区四那| 18禁在线播放成人免费| 少妇高潮的动态图| 午夜老司机福利剧场| 色哟哟哟哟哟哟| a级毛色黄片| 18禁黄网站禁片免费观看直播| 精品久久久久久成人av| 亚洲精品色激情综合| 中文字幕av在线有码专区| 69人妻影院| 在线天堂最新版资源| 国产成人aa在线观看| 久久国产乱子免费精品| 高清毛片免费看| 长腿黑丝高跟| 亚洲欧美成人精品一区二区| 日日摸夜夜添夜夜爱| 欧美日韩一区二区视频在线观看视频在线 | 国产日韩欧美在线精品| 国产黄a三级三级三级人| 麻豆成人av视频| 亚洲五月天丁香| 99久久久亚洲精品蜜臀av| 日本黄色片子视频| 国产精品av视频在线免费观看| 国产综合懂色| 99在线人妻在线中文字幕| 欧美+日韩+精品| 欧美在线一区亚洲|