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

    從同震和震后形變分析1668年M8.5級郯城地震對周邊地震活動性的影響

    2016-12-05 02:28:57孫玉軍吳中海
    地質力學學報 2016年3期
    關鍵詞:郯城同震庫侖

    徐 昊,孫玉軍,吳中海

    (中國地質科學院地質力學研究所,北京 100081)

    ?

    從同震和震后形變分析1668年M8.5級郯城地震對周邊地震活動性的影響

    徐 昊,孫玉軍,吳中海

    (中國地質科學院地質力學研究所,北京 100081)

    以山東郯城1668年大地震為例,以前人地表地質調查結果為約束,利用彈性位錯理論初步獲取了該地震的同震破裂模型;在此基礎上,基于粘彈性分層模型分析了該地震的同震和震后形變,同時以主震斷層為接收斷層計算了庫侖應力分布,進一步討論了地幔不同粘滯性系數(shù)對地表形變和庫侖應力變化的影響。計算結果顯示,該地震是一個右旋走滑為主兼有一定逆沖性質的地震,其同震位移巨大,能量釋放較徹底;同震破裂造成震中郯城縣西北、東北和南部部分斷層庫侖應力增加,而震后形變使得這些斷層庫侖應力進一步增加,在單縣、宿遷和日照等地,地震后350 a庫侖應力變化量達到+1bar—+1MPa量級;地幔粘滯性系數(shù)不同,形變量和庫侖應力變化達到穩(wěn)定的時間不同,但最終趨于穩(wěn)定的數(shù)值基本一致。

    郯城大地震;同震和震后形變;庫侖應力變化;粘滯性

    0 引言

    郯廬斷裂帶是中國東部一條重要的活動斷裂帶,歷史上沿該斷裂帶發(fā)生過多次大地震。1668年7月25日晚,在山東南部發(fā)生了一次曠古未有的特大地震,根據歷史資料得到該地震的震級約為Ms8.5級,極震區(qū)位于山東省郯城、臨沭、臨沂交界處,震中位置(34.8°N,118.5°E),極震區(qū)烈度達Ⅻ級[1]。由于極震區(qū)大部分位于郯城縣境內,故稱為郯城地震。這是我國大陸東部板塊內部最強烈的地震,造成巨大財產損失和人員傷亡。

    關于1668郯城地震,我國學者已經對其進行了相關研究,分析了這次地震發(fā)震斷層的構造特征和破裂機制[2~3],采用不同手段得到了此次地震震源參數(shù)[4~6],估算了郯城8.5級地震的復發(fā)周期[7~8]。但關于此次地震造成的應力場變化研究很少。近年來,越來越多的研究表明,地震發(fā)生會隨之釋放斷層上積聚的應力,但應力并不會憑空消失,部分應力轉移到其他地區(qū),導致應力積聚,進而誘發(fā)后續(xù)地震,此即地震的應力觸發(fā)理論[9]。根據應力觸發(fā)理論,庫侖應力增強相當于斷層額外負荷的增加,容易誘發(fā)后續(xù)地震;反之,庫侖應力減弱的區(qū)域——應力影區(qū),斷層負荷部分卸載,地震發(fā)生的概率降低。1668年郯城8.5級特大地震的同震及震后應力特征如何?其造成的應力場變化對周圍地區(qū)的地震活動性有何影響?本文圍繞這些問題,利用彈性位錯理論和分層巖石圈模型,計算了郯城1668年地震同震和震后形變以及庫侖應力的變化,探討本次地震對郯廬斷裂帶及附近區(qū)域地震活動性的影響,以期為進一步研究該區(qū)域的地震活動性提供參考。

    1 區(qū)域構造背景

    地震的發(fā)生與活動斷層聯(lián)系密切,斷層活動是地震產生的原因。郯廬斷裂帶是中國大陸東部一條北北東向的深大斷裂,也是中國東部有巨大影響的強震活動帶,史料記載曾發(fā)生過公元前70年安丘7級地震、1597年渤海7級地震、1668年郯城8.5級地震、1888年渤海7.5級地震、1969年渤海7.4級地震、1975年海城7.3級地震等[10]。其中1668年郯城8.5級地震的地震斷層為一高角度右旋走滑逆斷層,呈北北東向展布(見圖1),斷層面向東傾[6]。

    圖1 研究區(qū)構造背景圖(橢圓內加粗斷裂為此次地震主震斷層)Fig.1 Tectonic background map of Tanlu fault and its surrounding regions

    2 計算方法

    2.1 分層模型

    用于地震研究的彈性位錯模型主要可分為3類,第一類主要基于均勻彈性半空間介質模型[11],第二類主要基于分層介質模型[12],第三類基于成層球對稱地球模型的球形位錯模型[13]。本文基于地殼的分層模型采用格林函數(shù)解求解這次地震引起的同震位移以及引起的應力場變化。震后很長時間,有許多影響震后應力變化的因素,如震后余滑、孔隙效應和粘性松弛等,其中震后粘性松弛效應對震后形變和應力起到非常重要的作用[14]。Wang[15]提出了利用正交歸一法計算地震應力場Green函數(shù)方法,并在此基礎上建立了粘彈松弛分層模型下的地震同震及震后形變模型,編寫了數(shù)值計算程序——PSGRN/PSCMP軟件,用該程序基于分層、重力作用下的彈性模型來研究由于斷層錯動引起的同震和震后形變。假設郯城及周邊地區(qū)震后彈性松弛的過程主要發(fā)生在30 km上地殼以下。震后粘彈松弛應力變化是指一定深度下,受高溫高壓環(huán)境的影響,巖石層的力學性質逐漸由脆性向粘塑性轉化。由于地震發(fā)生時破裂的速度很快,在短時間內粘彈性的中下地殼以及上地幔力學性質表現(xiàn)為彈性體,然后震后隨時間推移,中下地殼和上地幔所積累的應力和應變逐漸釋放,并向上傳輸?shù)綇椥缘纳系貧ぶ?,導致上地殼中應力狀態(tài)發(fā)生改變[16]。對于地殼和上地幔來說,粘滯性系數(shù)不同,對震后形變和庫侖應力變化影響不同。粘滯性系數(shù)越低,松弛時間越短,震后短期內形變和庫侖應力增加的都比較快,很快能達到穩(wěn)定狀態(tài);相反,粘滯性系數(shù)越高,震后短期形變和庫侖應力增加的都比較慢,但最終的穩(wěn)定值相對較高[17]。模型假定地殼厚度為30 km,視為完全彈性體,粘滯性系數(shù)取無窮大;研究區(qū)缺乏長期的形變觀測資料和相關的震后形變研究成果,模型復雜化需要更多的參數(shù)來約束,因此采用簡單的Maxwell體分析地幔粘彈性變形。由于未獲取研究區(qū)地幔準確粘滯性系數(shù),本文模型中計算所需的粘滯性系數(shù)取1.0×1018Pa·s,最后會分析不同地幔粘滯性系數(shù)條件下,計算得到的形變量和庫侖應力的差異,并結合李勇紅等[5]給出的郯城8.5級地震震源區(qū)附近地殼上地幔速度結構,給出最終的分層模型及參數(shù)(見表1)。

    表1 分層模型及參數(shù)

    2.2 庫侖應力變化的計算

    隨著斷層上應力不斷累積,最終剪應力能夠克服斷裂上的靜摩擦力,鎖閉的斷層會發(fā)生破裂,斷層的破裂符合庫倫破裂準則[18]。利用彈性位錯理論,可以計算地震同震破裂所引起的應力場變化[19]。所研究的接收斷層庫侖應力變化為:

    (1)

    式中:Δσf為計算斷層面上庫侖應力的變化,Pa;Δτs是斷層面上剪應力的變化,Pa;μ′為有效摩擦系數(shù);Δσn是計算斷層面上正應力的變化,Pa。關于μ′的取值,不同學者有所差別。King等[20]認為有效摩擦系數(shù)對庫侖應力變化的結果將產生中等大小的影響,本文在計算過程中取中間值0.4。

    通過主震斷層模型計算靜態(tài)庫侖應力不但能夠有效反映余震的發(fā)展態(tài)勢,同時也能夠獲取對周圍活動斷層的影響[17]。由于此次地震發(fā)生時間距今已有350 a,本文僅選取一條主震斷層計算庫侖應力。采用之前學者研究所獲得的斷層相關參數(shù),走向21.57°、傾角89.51°、滑移角148.9°;起始破裂坐標為34.5°N、118.35°E,斷層面起始深度及其實際寬度取0 km和50 km[21],斷層破裂長度為150 km[6]。

    3 計算結果與分析

    3.1 同震和震后形變

    1668年郯城8.5級主震斷層為一高角度右旋走滑兼有逆沖性質的斷層,呈北北東向展布,斷層面東傾,震后觀察到水平最大形變量和垂直最大形變量分別為10 m和3 m[6]。盡管計算過程是在半無限空間假設模型下得到的,未考慮相關地形的影響,但通過圖像反映出的形變,依然可以判斷主震斷層為一右旋走滑逆斷層。計算得到同震時刻地表水平變形量和350 a后地表水平變形量(見圖2a,2c),同震時刻地表垂直變形量和350 a后地表垂直變形量(見圖2b,2d),其中350 a以后水平方向和垂直方向變形量分別為8.33 m和2.73 m,與王華林等[6]觀測數(shù)據存在一定誤差,由于在模型構建中缺少斷層面上精準的滑移數(shù)據,因此認為斷層模型依然能夠在一定程度上反映當時的地震破裂變形。

    a—同震時刻地表水平方向形變量;b—同震時刻地表垂直方向形變量(向上為正);c—350 a后地表水平方向形變量;d—350 a后地表垂直方向形變量(向上為正)圖2 同震及震后0 km處地表形變量Fig.2 The coseismic and postseismic deformation of the surface

    對比地表處水平方向的形變可以發(fā)現(xiàn),同震時刻和350 a以后斷層右旋的性質沒有改變,但350 a后水平形變量稍大于同震時刻水平形變量,距離主震斷層越遠增量越大。對比地表處垂直方向的形變可以發(fā)現(xiàn),同震時刻和350 a以后斷層逆沖的性質也未改變,但局部影區(qū)形變方向發(fā)生了倒轉,以右下角區(qū)域為例,同震時刻位移方向垂直向上,350 a后則向下變形。推測以上兩種變化原因與地殼巖石粘滯性有關:同震時刻斷層兩盤發(fā)生水平方向錯動,從而使距離破裂帶較近的區(qū)域瞬間獲得較大的形變量,而震后距破裂帶較遠的區(qū)域水平形變在巖石粘滯性作用下緩慢增加;已知主震斷層滑移角為149.5°,同震時刻斷層上下盤沿著走向錯動,必然會引起所在盤沿滑移方向一側的抬升,因此同震時刻右下和左上形變均為向上方向;隨著時間的推移,由于巖石粘滯性影響,變形逐漸恢復甚至引起反向形變,從而使圖中原先向上方向的變形轉變?yōu)橄蛳伦冃?。關于這種變化的量化分析,會在后來的部分進一步論述。

    3.2 庫侖應力變化

    每一次大地震發(fā)生都會引起周圍區(qū)域的應力增加或降低。一般來說,如果一個斷層的庫侖應力增加,則該斷層發(fā)生地震的危險性增強,反之則降低。同震破裂造成震中郯城縣西北、東北和南部部分地區(qū)庫侖應力增加,而震后形變使得這些地區(qū)庫侖應力進一步增加,在單縣、宿遷和日照等地,地震350 a后庫侖應力變化量達到+1 bar—+1 MPa量級(見圖3)。觀察研究區(qū)內震后350 a間Ms1—Ms7地震事件發(fā)現(xiàn),地震主要發(fā)生在這些庫侖應力增加的區(qū)域內。雖然地震的發(fā)生與構造運動直接相關,但這樣的特大型地震發(fā)生后對周圍斷層的應力積累影響十分明顯,而且計算結果顯示震后350 a該影響依然存在,只是目前已經趨于穩(wěn)定。

    同震庫侖應力增加最明顯區(qū)域主要位于斷層破裂帶的兩端及同震破裂不完全的區(qū)域。斷層破裂帶應力釋放明顯,而破裂帶的兩端則為應力集中的區(qū)域。本次地震的震源深度為20 km左右,因此,20 km深度同震應力釋放較地表同震應力釋放更為明顯(見圖3a,3b),沿著破裂帶主要為應力影區(qū);而在地表有部分滑動位移小,破裂不完全,應力增加的區(qū)域(見圖3a,3c)。對比同震時刻和350 a的同震庫侖應力分布,可以看出震后350 a在地表和20 km處應力增加的區(qū)域在面積和數(shù)值上都明顯增大。將本次地震后該區(qū)域的所有地震投影在平面后發(fā)現(xiàn),該區(qū)域的大部分地震事件也主要位于應力增加的區(qū)域。

    為了更清楚地觀察各應力影區(qū)的形變和庫侖應力變化趨勢,筆者選取南京、宿遷、單縣和日照4座城市計算并繪制其在深度20 km處350 a間位移和庫侖應力變化曲線。南京和日照在地震發(fā)生后庫侖應力都減少,不同的是南京在地震發(fā)生約20 a后庫侖應力逐漸增加,至今增加到-25 kPa左右(見圖4a);而日照在地震反生后庫侖應力稍微減少,隨后一直增加,到后來變化值穩(wěn)定在+1000 kPa左右(見圖4b)。宿遷和單縣則在地震發(fā)生后庫侖應力一直增加,最終宿遷庫侖應力增加穩(wěn)定在+1300 kPa左右(見圖4c),而單縣僅增加了+150 kPa左右(見圖4d)。庫侖應力增量表明宿遷和日照等地與郯廬斷裂帶主震斷層相同性質的斷層地震危險性增強。

    考慮同震和震后粘彈松弛效應對應力場變化的影響時,地幔的粘滯性系數(shù)不同取值可能會改變研究區(qū)的應力場分布。由于U(x)反映南北方向形變量,斷層走向為21.57°且為右旋斷層,U(x)的值能清楚反映斷層的運動性質。本文選取不同地幔粘滯性系數(shù)計算了郯城地震發(fā)生后350 a內,南京市深度20 km處的U(x)和庫侖應力變化值(見圖5),結果表明,震后形變量和庫侖應力變化量初始值相同,受粘滯性系數(shù)影響較小。隨著時間的推移,不管粘滯性大小,最終形變量穩(wěn)定值都近似相同,均為-1 m,但粘滯性系數(shù)越小,U(x)越快達到穩(wěn)定值,當粘滯性系數(shù)為1×1017Pa·s時,10 a左右就可達到穩(wěn)定值,當粘滯性系數(shù)為7×1018Pa·s時,則需要300 a左右才能達到穩(wěn)定值(見圖5a)。對于庫侖應力變化來說,震后初值變化也都相同,均為-17.5 kPa,曲線變化模式類似,都是庫侖應力先減少后增加最后趨于穩(wěn)定;當粘滯性系數(shù)為1×1017Pa·s時,50 a左右即可達到穩(wěn)定值,約為-24 kPa;當粘滯性系數(shù)為1×1018Pa·s時,150 a左右達到穩(wěn)定,約為-22.5 kPa;當粘滯性系數(shù)為7×1018Pa·s時,350 a時仍未達到穩(wěn)定狀態(tài),但可以估計其達到穩(wěn)定狀態(tài)的值基本相同(見圖5b)。

    a—地表處同震庫侖應力變化分布;b—20 km深處同震庫侖應力變化分布;c—地表處震后350 a庫侖應力變化分布;d—20 km深處震后350 a庫侖應力變化分布圖3 同震庫侖應力變化分布圖(灰色圓圈為1668年郯城地震后至今研究區(qū)內發(fā)生的M1—M7地震事件)Fig.3 The distribution map of coseismic Coulomb stress change

    a—南京市;b—日照市;c—宿遷市;d—單縣市△CFS—庫侖應力變化量;U(x)—南北向位移(北方向為正);U(y)—東西向位移(東方向為正);U(z)—垂向位移(向上為正)圖4 研究區(qū)內不同城市震后350 a內形變量和庫侖應力變化Fig.4 Deformation and Coulomb stress change in different cities in study regions 350 years after the earthquake

    圖5 不同地幔粘滯性系數(shù)對應的震后形變和庫侖應力變化(南京市)Fig.5 The postseismic deformation and Coulomb stress change of different mantle viscosities(Nanjing)

    4 結論

    本文基于粘彈松弛分層模型下的地震同震及震后形變模型,采用數(shù)值計算程序——PSGRN/PSCMP軟件,利用前人地表觀測和反演的地震數(shù)據計算了主震斷層同震以及震后的地表形變和庫侖應力變化量。

    通過對比同震和350 a后震中郯城縣東南部分位移形變量,從垂直向上的位移量轉變?yōu)橄蛳挛灰屏?,可知震后最終位移方向在受地殼粘彈性效應的影響下會發(fā)生倒轉。

    以主震斷層為接收斷層計算了同震和震后庫侖應力。同震庫侖應力增加的區(qū)域主要位于斷層破裂帶的兩端及同震破裂不完全的區(qū)域,具體分布于震中郯城縣西北、東北和東南的部分區(qū)域,震后形變使得這些地區(qū)庫侖應力進一步增加,在單縣、宿遷和日照等地,地震350 a后庫侖應力變化量達到+1 bar—+1 MPa量級。其他地區(qū)庫侖應力降低,南京地震350 a后庫侖應力變化量為-25 kPa。震后該區(qū)域發(fā)生的地震主要分布于庫侖應力增加的區(qū)域。

    通過對比地幔不同粘滯性系數(shù)對震后形變和庫侖應力的影響,發(fā)現(xiàn)對于同一區(qū)域,粘滯性系數(shù)越小,震后形變量和庫侖應力變化達到穩(wěn)定的時間越短,反之則時間越長;不同的地幔粘滯性系數(shù)對最終穩(wěn)定的形變量和庫侖應力影響不大。

    在討論此次地震與區(qū)域內周邊地震活動性關系時,主要依據庫侖應力變化量,雖然地震的發(fā)生與構造運動直接相關,但這樣的特大型地震發(fā)生后對周圍斷層的應力積累影響十分明顯,而且計算結果顯示震后350 a后該影響依然存在,但目前已經趨于穩(wěn)定,實際地震的發(fā)生還要綜合考慮構造活動等其他因素。此外由于1668—1970年之間區(qū)域地震事件缺少深度數(shù)據,所以在考慮后續(xù)地震在0 km和20 km處分布與庫侖應力變化的關系時誤差較大。

    [1] 馬玉香,鐘普裕.1668年山東郯城8.5級地震綜述[J].國際地震動態(tài),2009,(2):9~18.

    MA Yu-xiang,ZHONG Pu-yu.Summarization of the Tancheng Earthquake with 8.5 in 1668[J]. Recent developments in World Seismology, 2009, (2): 9~18.

    [2] 高維明,鄭朗蓀,李家靈,等.1668年郯城8.5級地震的發(fā)震構造[J].中國地震,1988,4(3):9~15.

    GAO Wei-ming,ZHENG Lang-sun,LI Jia-ling.The triggering structure of theMs8.5 earthquake of Tan Cheng in 1668[J].Earthquake Research in China,1988, 4(3): 9~15.

    [3] 鞠林雪.1668年郯城8.5級地震發(fā)震斷層構造特征與應力場分析[D].合肥:合肥工業(yè)大學,2012.

    JU Lin-xue.The research on the structure of triggring faults and stress field in 1668: The Tancheng earthquake[D]. Hefei: Hefei University of Technology, 2012.

    [4] 李清河,張元生,鮑海英,等.據地殼速度結構推測1668年山東郯城8.5級大地震震中[J].中國地震,2014,30(1):30~42.

    LI Qing-he, ZHANG Yuan-sheng, BAO Hai-ying, et al. Crust veolcity structure of the epicenter ofM8.5 earthquake at Tancheng, Shandong, China, in 1668[J]. Earthquake Research in China, 2014, 30 (1): 30~42.

    [5] 李永紅,胡新亮,許萍,等.以雙差定位方法對郯城8.5級地震震中附近現(xiàn)代小震重新定位[J].西北地震學報,2011,33(1):71~75.

    LI Yong-hong, HU Xin-liang, XU Ping, et al. Relocation of modern small earthquakes near the epicenter of Tancheng 8.5 earthquake using the Double Difference Earthquaake Location Algorithm[J]. Northwestern Seismological Journal, 2011, 33(1): 71~75.

    [6] 王華林,耿杰.1668年郯城8.5級地震震源參數(shù)及其討論[J].地震學刊,1996,(4):27~33.

    WANG Hua-lin, GENG Jie. Discussion about focus parameters of Tancheng earthquake ofMs8.5 in 1668[J]. Journal of Seismology, 1996, (4): 27~33.

    [7] 王華林.1668年郯城8.5級地震斷裂的全新世滑動速率,古地震和強震復發(fā)周期[J].西北地震學報,1995,17(4):1~12,22.

    WANG Hua-lin. Holocene displacement rate,paleoearthquakes and recurrence intervals of strong earthquakes along the 1668 Tancheng Earthquake (Ms=8.5) Fault[J]. Northwestern Seismological Journal, 1995, 17(4): 1~12.

    [8] 王啟鳴.三維有限元方法估算郯城8.5級大震重復周期[J].中國地震,1988,4(3):54~58.

    WANG Qi-ming. The calculation of recurrence intervals of theMs8.5 earthquake in Tancheng by the three dimension finite element method[J]. Earthquake Research in China, 1988, 4(3): 54~58.

    [9] 萬永革,吳忠良,周公威,等.地震應力觸發(fā)研究[J].地震學報,2002,24(5):533~551.

    WAN Yong-ge, WU Zhong-liang, Zhou Gong-wei, et al. Research on seismic stress triggering[J]. Acta Seismologica Sinica, 2002, 24(5): 533~551.

    [10] 魏光興.郯廬帶地震活動性研究[M].北京:地震出版社,1993.

    WEI Guang-xing. Research on the seismic activity of Tanlu earthquake belt[M]. Beijing: Seismological Press, 1993.

    [11] Okada Y. Surface deformation to shear and tensile faults in a halfspace[J]. Bulletin of the Seismological Society of America, 1985, 75 (4): 1135~1154.

    [12] Pollitz F F. Postseismic relaxation theory on a spherical Earth[J]. Bulletin of the Seismological Society of America, 1992, 82(1): 422~453.

    [13] 付廣裕,孫文科.球體位錯理論計算程序的總體設計與具體實現(xiàn)[J].地震,2012,32(2):73~87.

    FU Guang-yu, SUN Wen-ke. Overall design and special strctures of the computing codes for coseismic deformations on a layered sperical earth[J]. Earthquake, 2012, 32(2): 73~87.

    [14] 石耀霖.關于應力觸發(fā)和應力影概念在地震預報中應用的一些思考[J].地震,2001, 21(3):1~7.

    SHI Yao-lin. Sress triggers and stress shadows How to apply these concepts to earthquake prediction[J]. Earthquake, 2001, 21(3): 1~7.

    [15] Wang R, Lorenzo-Martín F, Roth F. PSGRN/PSCMP: A new code for calculating co-and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory[J]. Computers & Geosciences, 2006, 32: 527~541.

    [16] Freed A M. Earthquake triggering by static, dynamic, and postseismic stress transfer[J]. Annual Review of Earth & Planetary Sciences, 2005, 33: 335~367.

    [17] 孫玉軍,董樹文,范桃園,等.從同震和震后形變分析日本東北Mw9.0級大地震對近場地震活動性的影響[J].地球物理學進展,2013,28(3):1131~1139.

    SUN Yu-jun, DONG Shu-wen, FAN Tao-yuan, et al. The effect of TohotuMw9.0 earthquake on the near field seismic activity from coseismic and postseismic deformation[J]. Progress in Geophys, 2013, 28(3): 1131~1139.

    [18] 江國焰.利用庫侖破裂準則評估活動斷層地震危險性[D].武漢:武漢大學,2013.

    JIANG Guo-yan. Assessing the seismic hazards of active faults using Coulomb Failure Criterion[D]. Wuhan: Wuhan University, 2013.

    [19] Okada Y. Internal deformation due to shear and tensile fault in a half space[J]. Bulletin of the Seismological Society of America, 1992, 92: 1018~1040.

    [20] King G C P, Stein R S, Lin J. Static stress changes and the triggering of earthquakes[J]. Bull.seismol.soc.amer, 1994, 84: 935~953.

    [21] 周翠英,刁桂苓,耿杰,等.1668年郯城大地震震源斷層三維特征反演[J].地球物理學進展,2013,28(6):2814~2824.

    ZHOU Cui-ying, DIAO Gui-ling, GENG jie, et al. 3-D characteristics inversion of hypocenter fault-plane of the 1668 Tancheng great earthquake[J]. Progress in Geophys, 2013, 28(6): 2814~2824.

    THE EFFECT OF 1668 TANCHENG M8.5 EARTHQUAKE ON THE SEISMIC ACTIVITY OF THE VICINITY FROM COSEISMIC AND POSTSEISMIC DEFORMATION

    XU Hao, SUN Yu-jun, WU Zhong-hai

    (Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100081, China)

    The Tanlu fault zone is an important active fault zone in east China. Many great earthquakes have occurred along it. The 1668M8.5 Tancheng earthquake is the largest recorded event. With the constraint of surface rupture, we get a coseismic rupture model based on the elastic dislocation theory. In addition, we calculate the coseismic, postseismic deformation and Coulomb stress change of this earthquake with a viscoelastic multilayered model. The results show that the rupture fault of the earthquake is adextral strike-slip fault, which has slight thrust. The coseismic displacement of this event is very large. The accumulated energy is released thoroughly. The coseismic rupture of the 1668 Tancheng earthquake increased the Coulomb stress on the northwest, northeast and southeast faults of Tanlu fault zone. The Coulomb stress increases further due to the postseismic deformation. After 350 years of the Tancheng earthquake, the Coulomb stress changes in Shanxian, Suqian and Rizhaoget to +1bar-+1MPa. When the mantle is imposed different viscosities, the duration time which the deformation and Coulomb changes become stable will be different. However, the stable values of them are generally same.

    1668 TanchengM8.5 earthquake; coseismic and postseismic deformation; Coulomb stress change; viscosity

    1006-6616(2016)03-0568-09

    2016-04-24

    中國地質調查局地質調查項目“長江經濟帶活動構造與區(qū)域地殼穩(wěn)定性調查”(DD20160268)

    徐昊(1993-),男,碩士研究生,研究方向地球動力學。E-mail:xhdky2011@163.com

    孫玉軍(1983-),男,博士,副研究員,主要從事地球動力學方面的研究。E-mail:sunyujunabc@163.com

    P315.5

    A

    猜你喜歡
    郯城同震庫侖
    1976年唐山強震群震后庫侖應力演化及其與2020年古冶5.1級地震的關系
    地震研究(2021年1期)2021-04-13 01:04:46
    My Mother Taught Me About……
    山東郯城:多肉種植成花卉苗木產業(yè)“新寵”
    云南思茅大寨井水位地震同震響應特征分析*
    地震研究(2018年4期)2018-11-23 02:29:36
    2種模式下2001—2016年郯城主要糧食作物氣象適宜度對比分析
    基于粘彈庫侖應力變化的后續(xù)最大地震震級估計及2008、2014年于田2次7.3級地震之間關系的討論
    中國地震(2015年1期)2015-11-08 11:11:18
    蘆山地震前后介質波速變化與GPS應變場相關性研究?
    蘆山Ms7.0地震引起的水位同震響應特征分析
    四川地震(2014年2期)2014-12-02 04:16:30
    一種周期庫侖作用勢優(yōu)化法的改進
    計算物理(2014年1期)2014-03-11 17:01:03
    長程庫侖勢對高溫超導渦旋電荷的影響
    男人的好看免费观看在线视频| 操出白浆在线播放| 午夜福利在线在线| 亚洲成人中文字幕在线播放| 日韩欧美 国产精品| 午夜福利在线在线| 男女做爰动态图高潮gif福利片| 欧美黑人欧美精品刺激| 国产精品一区二区三区四区久久| 十八禁网站免费在线| 九九热线精品视视频播放| 黄色视频,在线免费观看| 天天躁日日操中文字幕| 国产精品女同一区二区软件 | 亚洲国产中文字幕在线视频| 亚洲av五月六月丁香网| 午夜福利欧美成人| 69av精品久久久久久| xxxwww97欧美| 日韩欧美在线乱码| 亚洲精品亚洲一区二区| 美女高潮喷水抽搐中文字幕| a在线观看视频网站| 欧美中文综合在线视频| 国产乱人伦免费视频| 91av网一区二区| 欧美性感艳星| 国产真实乱freesex| 亚洲在线自拍视频| 12—13女人毛片做爰片一| 欧美黄色片欧美黄色片| 久久久久国内视频| 亚洲人成网站高清观看| 色哟哟哟哟哟哟| 97超视频在线观看视频| 给我免费播放毛片高清在线观看| 国产色婷婷99| 精品电影一区二区在线| 女生性感内裤真人,穿戴方法视频| 真人一进一出gif抽搐免费| 午夜久久久久精精品| 一个人看的www免费观看视频| av国产免费在线观看| 国产伦人伦偷精品视频| 色av中文字幕| h日本视频在线播放| 少妇丰满av| 男女之事视频高清在线观看| 欧美极品一区二区三区四区| 波野结衣二区三区在线 | 99久久精品国产亚洲精品| 小蜜桃在线观看免费完整版高清| 欧美日本亚洲视频在线播放| av在线天堂中文字幕| 亚洲av二区三区四区| 男人舔奶头视频| 久久精品综合一区二区三区| 桃色一区二区三区在线观看| 成人高潮视频无遮挡免费网站| 搡老岳熟女国产| 国内毛片毛片毛片毛片毛片| 亚洲人成网站在线播| 国产毛片a区久久久久| 3wmmmm亚洲av在线观看| 午夜老司机福利剧场| 国产又黄又爽又无遮挡在线| 久久久久久人人人人人| 国产精品精品国产色婷婷| 免费在线观看成人毛片| 欧美日韩瑟瑟在线播放| avwww免费| 欧美成人性av电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩综合久久久久久 | 国产一级毛片七仙女欲春2| 99精品欧美一区二区三区四区| 午夜免费激情av| 日本在线视频免费播放| 国内精品美女久久久久久| 十八禁网站免费在线| 熟女人妻精品中文字幕| 99国产精品一区二区三区| 亚洲精品成人久久久久久| 最近最新中文字幕大全免费视频| 国产精品久久电影中文字幕| 动漫黄色视频在线观看| 国产高清视频在线观看网站| 午夜激情福利司机影院| 亚洲最大成人中文| 久久婷婷人人爽人人干人人爱| 欧美xxxx黑人xx丫x性爽| 日韩欧美国产在线观看| 一进一出抽搐gif免费好疼| 久久精品综合一区二区三区| 欧美成人免费av一区二区三区| 手机成人av网站| 热99在线观看视频| 男人和女人高潮做爰伦理| 日日夜夜操网爽| 亚洲精品亚洲一区二区| 一个人免费在线观看电影| 精品日产1卡2卡| 亚洲av五月六月丁香网| 亚洲精品美女久久久久99蜜臀| 亚洲人成电影免费在线| 少妇的逼水好多| 久久人人精品亚洲av| 日日干狠狠操夜夜爽| 国产精品永久免费网站| 国产高清视频在线播放一区| 亚洲狠狠婷婷综合久久图片| 亚洲av不卡在线观看| 黄色日韩在线| 麻豆国产97在线/欧美| 亚洲人成伊人成综合网2020| 精品福利观看| 亚洲av免费在线观看| 国产美女午夜福利| 天天一区二区日本电影三级| 欧美黄色片欧美黄色片| 亚洲成a人片在线一区二区| 亚洲欧美日韩东京热| 久久婷婷人人爽人人干人人爱| 国产亚洲av嫩草精品影院| 美女被艹到高潮喷水动态| 神马国产精品三级电影在线观看| 色播亚洲综合网| 久久人妻av系列| 久久久久久人人人人人| 最近最新中文字幕大全免费视频| 国产色婷婷99| 亚洲精品成人久久久久久| 久久久久国内视频| 国产欧美日韩精品一区二区| 五月伊人婷婷丁香| 精品欧美国产一区二区三| 午夜视频国产福利| 亚洲片人在线观看| 成人av在线播放网站| 国产亚洲精品一区二区www| 午夜福利欧美成人| 日韩有码中文字幕| 婷婷亚洲欧美| 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| 在线观看免费视频日本深夜| 天天添夜夜摸| 一本精品99久久精品77| 日本与韩国留学比较| 久久久久久久精品吃奶| 免费人成在线观看视频色| 国产av在哪里看| 久久久久久久精品吃奶| 精品久久久久久久人妻蜜臀av| 日韩欧美精品v在线| 亚洲五月天丁香| 亚洲欧美日韩高清在线视频| 亚洲精品日韩av片在线观看 | 国产精品乱码一区二三区的特点| 国产av不卡久久| 国产精品亚洲av一区麻豆| 免费av不卡在线播放| 十八禁网站免费在线| 美女高潮的动态| 白带黄色成豆腐渣| 国产精品爽爽va在线观看网站| 99久久综合精品五月天人人| 亚洲成人久久性| 天美传媒精品一区二区| 淫妇啪啪啪对白视频| 日本五十路高清| av视频在线观看入口| 最近最新中文字幕大全电影3| 日韩人妻高清精品专区| 精品久久久久久久末码| 啦啦啦韩国在线观看视频| 村上凉子中文字幕在线| 久久久色成人| 免费无遮挡裸体视频| 免费高清视频大片| 最新在线观看一区二区三区| 午夜两性在线视频| 天天添夜夜摸| 午夜精品一区二区三区免费看| 精品人妻1区二区| 一个人观看的视频www高清免费观看| 久久久精品大字幕| 男人舔女人下体高潮全视频| 在线观看免费午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久| 色吧在线观看| 免费看日本二区| 熟女电影av网| 51午夜福利影视在线观看| 亚洲久久久久久中文字幕| 国产免费一级a男人的天堂| 亚洲av一区综合| 内地一区二区视频在线| 日本五十路高清| 国模一区二区三区四区视频| 免费av观看视频| 午夜日韩欧美国产| 99久久精品热视频| 色视频www国产| av天堂中文字幕网| 久久久国产精品麻豆| 中出人妻视频一区二区| 国产高潮美女av| 欧美zozozo另类| 真实男女啪啪啪动态图| 美女大奶头视频| 嫁个100分男人电影在线观看| 天堂av国产一区二区熟女人妻| 一级a爱片免费观看的视频| 俺也久久电影网| 91九色精品人成在线观看| 色综合亚洲欧美另类图片| 国产中年淑女户外野战色| 婷婷精品国产亚洲av| 久久久久久久午夜电影| 观看免费一级毛片| 亚洲美女视频黄频| 首页视频小说图片口味搜索| 色哟哟哟哟哟哟| 女人十人毛片免费观看3o分钟| 国产色婷婷99| 18禁国产床啪视频网站| 三级毛片av免费| 欧美区成人在线视频| 久久久国产精品麻豆| 在线播放国产精品三级| 美女大奶头视频| 精品熟女少妇八av免费久了| 国产老妇女一区| 偷拍熟女少妇极品色| 国产精品一区二区免费欧美| 久久婷婷人人爽人人干人人爱| 非洲黑人性xxxx精品又粗又长| 99视频精品全部免费 在线| 免费观看人在逋| 看免费av毛片| 宅男免费午夜| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品日韩av片在线观看 | 亚洲人成网站在线播放欧美日韩| 亚洲aⅴ乱码一区二区在线播放| 国产激情欧美一区二区| 欧美黄色片欧美黄色片| 精品久久久久久久久久久久久| 99国产极品粉嫩在线观看| 性色avwww在线观看| 国产伦在线观看视频一区| 天堂动漫精品| 国产色爽女视频免费观看| 91av网一区二区| 亚洲人与动物交配视频| 身体一侧抽搐| 午夜免费激情av| 久久久久国产精品人妻aⅴ院| 国产黄a三级三级三级人| 欧洲精品卡2卡3卡4卡5卡区| 3wmmmm亚洲av在线观看| 狂野欧美激情性xxxx| 国产高清视频在线播放一区| 国产精品嫩草影院av在线观看 | 日韩欧美国产在线观看| 九九久久精品国产亚洲av麻豆| 日本撒尿小便嘘嘘汇集6| 成人特级黄色片久久久久久久| 国产高清有码在线观看视频| 欧美绝顶高潮抽搐喷水| 免费大片18禁| 中文字幕人妻熟人妻熟丝袜美 | 国产精品精品国产色婷婷| 老司机深夜福利视频在线观看| 亚洲人成网站在线播放欧美日韩| 免费大片18禁| 三级毛片av免费| 免费看光身美女| 亚洲片人在线观看| 亚洲精品影视一区二区三区av| 国产 一区 欧美 日韩| 99热精品在线国产| 国产真人三级小视频在线观看| 久9热在线精品视频| 婷婷亚洲欧美| 99在线人妻在线中文字幕| 午夜免费激情av| 午夜久久久久精精品| 亚洲精品一区av在线观看| 18美女黄网站色大片免费观看| 国产高清三级在线| 国产三级中文精品| 两个人视频免费观看高清| 真实男女啪啪啪动态图| 男人舔奶头视频| 久久久久久久精品吃奶| 久久久精品大字幕| 亚洲乱码一区二区免费版| 午夜福利在线观看吧| 窝窝影院91人妻| 一级a爱片免费观看的视频| 1024手机看黄色片| 一级毛片女人18水好多| 色在线成人网| 69av精品久久久久久| 亚洲在线观看片| 亚洲精品一区av在线观看| 亚洲欧美一区二区三区黑人| 在线播放国产精品三级| 欧美一区二区精品小视频在线| 性色av乱码一区二区三区2| 国产在视频线在精品| 欧美一区二区精品小视频在线| 男女之事视频高清在线观看| 日韩欧美国产在线观看| 亚洲精品一区av在线观看| 亚洲精品影视一区二区三区av| 亚洲中文日韩欧美视频| av欧美777| 精品一区二区三区视频在线 | 色综合亚洲欧美另类图片| 99热6这里只有精品| 欧美日韩瑟瑟在线播放| 午夜福利免费观看在线| 中文字幕久久专区| 亚洲18禁久久av| 最近最新中文字幕大全电影3| 九九久久精品国产亚洲av麻豆| 成年女人毛片免费观看观看9| 91字幕亚洲| 午夜福利视频1000在线观看| 老司机在亚洲福利影院| 色播亚洲综合网| 又黄又爽又免费观看的视频| 亚洲国产精品成人综合色| 老司机深夜福利视频在线观看| 亚洲性夜色夜夜综合| 少妇的逼水好多| 久久久久久国产a免费观看| 国产精品电影一区二区三区| 国产精品日韩av在线免费观看| 天堂网av新在线| 一进一出好大好爽视频| 久久国产精品人妻蜜桃| 国产成年人精品一区二区| 久久国产精品人妻蜜桃| 亚洲国产欧洲综合997久久,| 高潮久久久久久久久久久不卡| 国产精品免费一区二区三区在线| 国产高清三级在线| 天堂√8在线中文| 天堂网av新在线| 波多野结衣高清作品| 国产在线精品亚洲第一网站| 国产av不卡久久| 啦啦啦韩国在线观看视频| 免费高清视频大片| 窝窝影院91人妻| 国产成人福利小说| 亚洲av中文字字幕乱码综合| 乱人视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 天美传媒精品一区二区| 日韩欧美国产一区二区入口| 99热6这里只有精品| 人人妻,人人澡人人爽秒播| 别揉我奶头~嗯~啊~动态视频| 每晚都被弄得嗷嗷叫到高潮| 午夜精品久久久久久毛片777| 日本黄大片高清| 午夜福利视频1000在线观看| 悠悠久久av| 1024手机看黄色片| 国产三级中文精品| 久久精品国产99精品国产亚洲性色| 最好的美女福利视频网| 特大巨黑吊av在线直播| 黄色女人牲交| 国产精品一区二区免费欧美| 国产精品三级大全| 香蕉久久夜色| 免费人成视频x8x8入口观看| 欧美区成人在线视频| 亚洲五月婷婷丁香| 两个人的视频大全免费| 亚洲成av人片免费观看| 日韩欧美 国产精品| 国产美女午夜福利| 午夜福利18| 精品午夜福利视频在线观看一区| 99久久精品热视频| 亚洲欧美日韩高清在线视频| 亚洲无线在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 色播亚洲综合网| 亚洲专区国产一区二区| 禁无遮挡网站| 好男人电影高清在线观看| 一级黄片播放器| 国产三级在线视频| 亚洲国产中文字幕在线视频| 精品乱码久久久久久99久播| 看免费av毛片| 亚洲性夜色夜夜综合| 老司机深夜福利视频在线观看| 亚洲欧美日韩卡通动漫| 一级黄色大片毛片| 久久久国产成人免费| 日韩人妻高清精品专区| 在线看三级毛片| 美女免费视频网站| 国产av一区在线观看免费| 日韩欧美三级三区| 国产欧美日韩精品一区二区| 免费看美女性在线毛片视频| 网址你懂的国产日韩在线| a在线观看视频网站| 老汉色∧v一级毛片| 51国产日韩欧美| 女人高潮潮喷娇喘18禁视频| 亚洲第一欧美日韩一区二区三区| 日本 av在线| 久久人妻av系列| 国产黄a三级三级三级人| 成人精品一区二区免费| 性色av乱码一区二区三区2| 老司机在亚洲福利影院| 男人和女人高潮做爰伦理| 伊人久久大香线蕉亚洲五| 亚洲一区高清亚洲精品| 久久这里只有精品中国| 蜜桃亚洲精品一区二区三区| 亚洲精品亚洲一区二区| 麻豆国产av国片精品| 亚洲无线观看免费| 国产一区二区在线av高清观看| 国产成人欧美在线观看| 91字幕亚洲| 欧美色视频一区免费| 变态另类成人亚洲欧美熟女| 国产精品爽爽va在线观看网站| 免费人成在线观看视频色| 国产高清视频在线观看网站| 欧美国产日韩亚洲一区| 国产真实伦视频高清在线观看 | 精品人妻1区二区| 国产美女午夜福利| 三级国产精品欧美在线观看| 最后的刺客免费高清国语| 在线十欧美十亚洲十日本专区| 天天添夜夜摸| 国产激情偷乱视频一区二区| 男女床上黄色一级片免费看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲av二区三区四区| 欧美成人免费av一区二区三区| 日韩欧美精品v在线| 国产成人啪精品午夜网站| 亚洲,欧美精品.| 国产av麻豆久久久久久久| 超碰av人人做人人爽久久 | 亚洲av一区综合| 婷婷精品国产亚洲av| tocl精华| 午夜影院日韩av| 欧美国产日韩亚洲一区| 国产成人福利小说| 亚洲欧美日韩东京热| a在线观看视频网站| 精品人妻偷拍中文字幕| 精品午夜福利视频在线观看一区| 日韩欧美三级三区| 91九色精品人成在线观看| 夜夜爽天天搞| 99久久久亚洲精品蜜臀av| 97超级碰碰碰精品色视频在线观看| 啦啦啦免费观看视频1| 两个人看的免费小视频| 中文字幕熟女人妻在线| 亚洲专区国产一区二区| 无人区码免费观看不卡| 午夜福利欧美成人| 一级毛片女人18水好多| 精品免费久久久久久久清纯| 一区二区三区激情视频| 老司机午夜十八禁免费视频| 久久99热这里只有精品18| 窝窝影院91人妻| 18禁美女被吸乳视频| 亚洲人成网站在线播放欧美日韩| 99热6这里只有精品| av欧美777| 免费人成在线观看视频色| 日本免费a在线| 精品不卡国产一区二区三区| 精品无人区乱码1区二区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲人成网站高清观看| 香蕉丝袜av| 亚洲人成网站在线播放欧美日韩| 老司机午夜十八禁免费视频| 国产精品久久久人人做人人爽| 99久久精品一区二区三区| 久久6这里有精品| 三级国产精品欧美在线观看| 亚洲第一欧美日韩一区二区三区| 五月玫瑰六月丁香| 伊人久久精品亚洲午夜| 亚洲成av人片免费观看| 亚洲国产精品999在线| 九色国产91popny在线| 乱人视频在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲欧美激情综合另类| 国产精品一及| 一本综合久久免费| 国产亚洲精品久久久com| 亚洲精品一区av在线观看| 精品人妻一区二区三区麻豆 | 成人一区二区视频在线观看| 女人十人毛片免费观看3o分钟| 老汉色∧v一级毛片| 成年女人毛片免费观看观看9| 国产成年人精品一区二区| 国产精品99久久久久久久久| 亚洲av免费高清在线观看| 亚洲国产精品999在线| 淫妇啪啪啪对白视频| 天美传媒精品一区二区| 精品一区二区三区人妻视频| 特大巨黑吊av在线直播| 麻豆国产97在线/欧美| 老熟妇乱子伦视频在线观看| bbb黄色大片| 午夜a级毛片| 精品久久久久久久人妻蜜臀av| 母亲3免费完整高清在线观看| 日韩人妻高清精品专区| www.熟女人妻精品国产| 欧美性猛交╳xxx乱大交人| 亚洲成人免费电影在线观看| 国产精品亚洲av一区麻豆| 欧美黑人欧美精品刺激| 欧美精品啪啪一区二区三区| 欧美日韩福利视频一区二区| 男女做爰动态图高潮gif福利片| 少妇的逼水好多| 国模一区二区三区四区视频| 国产三级在线视频| 久久久精品大字幕| 99国产精品一区二区蜜桃av| 国产高潮美女av| 内地一区二区视频在线| 精品久久久久久久毛片微露脸| 精品不卡国产一区二区三区| 欧美黄色淫秽网站| 中文字幕人妻丝袜一区二区| 又黄又爽又免费观看的视频| 狂野欧美白嫩少妇大欣赏| 草草在线视频免费看| 午夜免费男女啪啪视频观看 | 波野结衣二区三区在线 | ponron亚洲| svipshipincom国产片| 淫秽高清视频在线观看| 国产精品久久久久久久电影 | 宅男免费午夜| 日韩欧美在线乱码| 久久性视频一级片| 午夜福利欧美成人| 午夜日韩欧美国产| 给我免费播放毛片高清在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲av中文字字幕乱码综合| 亚洲精品久久国产高清桃花| 色在线成人网| 99热精品在线国产| 精品久久久久久久末码| 午夜福利18| 成人一区二区视频在线观看| 中亚洲国语对白在线视频| 神马国产精品三级电影在线观看| 真实男女啪啪啪动态图| 欧美激情在线99| 欧美又色又爽又黄视频| xxx96com| 色播亚洲综合网| 日本免费一区二区三区高清不卡| 亚洲第一电影网av| 欧美一区二区亚洲| 一进一出抽搐gif免费好疼| 99久久成人亚洲精品观看| 在线观看美女被高潮喷水网站 | 亚洲精品成人久久久久久| 一级毛片高清免费大全| 俺也久久电影网| 中文字幕人成人乱码亚洲影| 一区二区三区国产精品乱码| 俺也久久电影网| 国产一区二区在线观看日韩 | 九色成人免费人妻av| 婷婷精品国产亚洲av| 性欧美人与动物交配| 成年女人永久免费观看视频| 日本成人三级电影网站| 91麻豆av在线| 噜噜噜噜噜久久久久久91| 身体一侧抽搐| 亚洲第一电影网av| 69人妻影院| 又粗又爽又猛毛片免费看| 黑人欧美特级aaaaaa片| 国产精品久久久久久人妻精品电影| 国产男靠女视频免费网站|