李洪麗, 劉財(cái),2, 田有,2,3* , 吳兆營(yíng), 朱洪翔

1 吉林大學(xué)地球信息探測(cè)儀器教育部重點(diǎn)實(shí)驗(yàn)室， 長(zhǎng)春 130026 2 吉林大學(xué)地球探測(cè)科學(xué)與技術(shù)學(xué)院， 長(zhǎng)春 130026 3 長(zhǎng)白山火山綜合地球物理教育部野外科學(xué)觀測(cè)研究站， 長(zhǎng)春 130021 4 吉林省地震局， 長(zhǎng)春 130021

0 引言

中國(guó)東北地區(qū)由多個(gè)微陸塊組成，主體位于西伯利亞地臺(tái)、華北克拉通與西太平洋板塊之間，構(gòu)造活動(dòng)復(fù)雜.該區(qū)的構(gòu)造演化在古生代期間主要受控于華北克拉通和西伯利亞地臺(tái)之間的古亞洲洋構(gòu)造體系演化的影響，以眾多微陸塊之間的碰撞-拼合和古亞洲洋的最終閉合為主要特征(唐杰等，2019).中生代以來(lái)，該地區(qū)經(jīng)歷了蒙古—鄂霍茨克閉合和中-新生代西太平洋板塊俯沖的雙重疊加和改造(田有等，2019).而東北亞大陸邊緣構(gòu)造演化(如長(zhǎng)白山火山、松遼盆地)自中生代以來(lái)主要受到西太平洋板塊向歐亞大陸巖石圈下大規(guī)模、多角度和多期次深俯沖作用的影響(唐克東等，2004；田有等，2019).

2013年10月，位于松原市的前郭爾羅斯地區(qū)發(fā)生了5.5級(jí)地震，震中位于查干花鎮(zhèn)東白音花村(東經(jīng)124.2°，北緯44.6°)，震源深度約11 km.長(zhǎng)春、哈爾濱等地有明顯震感，吉林市農(nóng)安縣有兩次明顯震感.截止到11月25日12時(shí)00分，共記錄到地震事件400多次，其中，10月31日又發(fā)生5.0級(jí)地震，震源深度6 km，11月23日發(fā)生5.8級(jí)地震，震源深度9 km，周邊市縣震感明顯，11月23日06時(shí)32分發(fā)生5.0級(jí)地震，震源深度8 km，4.0～4.9級(jí) 4次，3.0～3.9級(jí)10次.吉林省松原市官方2013年11月23日發(fā)布的統(tǒng)計(jì)，自10月31日至11月23日，松原前郭縣幾次地震已經(jīng)累計(jì)嚴(yán)重?fù)p毀房屋16210戶，直接經(jīng)濟(jì)損失201162萬(wàn)元人民幣.吉林油田就坐落在松原，而松原市是政治、經(jīng)濟(jì)、文化的中心，松原市一旦發(fā)生較大地震，將會(huì)對(duì)城市造成嚴(yán)重后果.因此，對(duì)于松原震區(qū)的發(fā)震機(jī)制研究具有重要的科學(xué)意義.

國(guó)內(nèi)自20世紀(jì)80年代開(kāi)始開(kāi)展地震層析成像研究.例如，眾多學(xué)者利用地震波走時(shí)數(shù)據(jù)對(duì)首都圈地區(qū)、西北造山帶、天山地區(qū)、昆侖造山帶、帕米爾等關(guān)鍵構(gòu)造區(qū)帶進(jìn)行了層析成像研究(金安蜀等，1980；胥頤等，2000；雷建設(shè)等，2002；趙文津等，2004；郭飚等，2004，2006；趙大鵬等2004；黃金莉和趙大鵬，2005；齊誠(chéng)等，2006；賀日政等，2006；錢(qián)輝等，2007，2011)，獲得了巖石圈尺度三維地震波速度結(jié)構(gòu)，探討了典型構(gòu)造區(qū)構(gòu)造演化特征，取得了重要研究成果.田有等(2011, 2019)對(duì)東北地區(qū)的地殼、上地幔三維結(jié)構(gòu)進(jìn)行研究，表明該區(qū)的構(gòu)造演化與西太平洋板塊的俯沖息息相關(guān).上述研究表明地震層析成像是研究深部殼幔結(jié)構(gòu)和發(fā)震構(gòu)造的最有效方法之一.

前人對(duì)于東北地區(qū)地震層析成像研究工作重點(diǎn)關(guān)注巖石圈尺度構(gòu)造演化、活火山深部起源等，而對(duì)東北地區(qū)的淺部精細(xì)速度結(jié)構(gòu)、破壞性地震發(fā)震機(jī)制及盆地內(nèi)部中上地殼高精度速度結(jié)構(gòu)等方面并未進(jìn)行深入研究.針對(duì)松原地震區(qū)震群的發(fā)震機(jī)制與該區(qū)速度結(jié)構(gòu)、波速比結(jié)構(gòu)等諸多科學(xué)問(wèn)題，本文收集松原前郭地震區(qū)及鄰區(qū)地震臺(tái)站采集的地震震相到時(shí)資料，采用雙差層析成像方法構(gòu)建松原震區(qū)中、上地殼三維速度結(jié)構(gòu)和泊松比模型，并基于反演得到的成像結(jié)果進(jìn)一步計(jì)算了該區(qū)地下介質(zhì)的裂縫密度和飽和度分布特征，揭示震區(qū)未來(lái)的地震危險(xiǎn)性，并且為地震監(jiān)測(cè)規(guī)劃、民居震害預(yù)防工作和震后恢復(fù)重建提供抗震指導(dǎo)建議，最終達(dá)到最大限度減輕地震災(zāi)害的目的.

1 數(shù)據(jù)資料與成像方法

1.1 數(shù)據(jù)資料

本次地震數(shù)據(jù)源自吉林省地震局布設(shè)在松原地區(qū)的17個(gè)地震臺(tái)站在2013年10月—2018年12月期間記錄的地震事件(圖1).震中位置比較集中(圖2)，其中絕大多數(shù)地震事件分布在東經(jīng)123.7°～124.5°、北緯44.5°～45.4°范圍內(nèi).在深度方向上，大部分事件分布在0～19 km范圍內(nèi)，有兩個(gè)事件震源深度大于20 km，最大震源深度為24 km.地震事件在空間上的分布形態(tài)近似為柱狀分布.

圖1 (a)松原地震區(qū)構(gòu)造、地震臺(tái)站和震群位置分布圖. 黑色三角代表地震臺(tái)站；圓點(diǎn)代表地震；穿過(guò)震群的實(shí)線代表速度剖面位置；虛線表示斷層；(b)研究區(qū)位置圖Fig.1 (a) Map showing tectonics map, seismic stations and earthquakes in the Songyuan area. Black triangles are seismic stations. Small circles are earthquakes. The solid lines passing through seismic events are locations of velocity profiles. Dashed solid lines are faults； (b) The location of the study area

圖2 松原地區(qū)地震重定位前(a—c)后(d—f)震源位置分布Fig.2 Hypocenters before (a—c) and after (d—f) relocation of earthquakes in Songyuan area

本研究共采用2926條P波走時(shí)數(shù)據(jù)和2665條S波走時(shí)數(shù)據(jù).P波和S波射線分布如圖3所示.從射線分布圖上看，在研究區(qū)中心位置，即震群位置射線非常密集，后面對(duì)于該成像分辨率檢測(cè)結(jié)果也證明在震群集中區(qū)域分辨率較高.雙差成像方法主要的優(yōu)點(diǎn)在于可以采用時(shí)間對(duì)進(jìn)行反演運(yùn)算，所以依據(jù)地震之間的距離對(duì)515個(gè)地震事件進(jìn)行了配對(duì)，其中共有415個(gè)事件滿足配對(duì)條件，獲得47870個(gè)P波事件對(duì)和46350個(gè)S波事件對(duì).

圖3 P波和S波射線分布Fig.3 Ray paths of P-wave and S-wave

1.2 震源定位與成像方法

雙差層析成像方法是基于雙差地震定位方法(Waldhauser and Ellsworth, 2000)發(fā)展而來(lái)，該方法利用地震對(duì)相對(duì)于同一臺(tái)站之間相對(duì)走時(shí)差的差來(lái)同時(shí)反演速度結(jié)構(gòu)和高精度震源位置(Zhang and Thurber, 2003).雙差層析成像方法在強(qiáng)震區(qū)精細(xì)結(jié)構(gòu)構(gòu)建、俯沖帶的高分辨率成像、火山的速度結(jié)構(gòu)反演等方面具有較好的效果.由于雙差層析成像方法是利用相鄰事件的相對(duì)走時(shí)進(jìn)行反演，所以在地震多發(fā)區(qū)域，如板塊俯沖帶，大斷層區(qū)域，發(fā)生的地震比較集中，容易產(chǎn)生配對(duì)而成像效果好.

圖4 松原震群重定位前(a)后(b)走時(shí)殘差統(tǒng)計(jì)直方圖Fig.4 Statistical histograms of travel time residuals before (a) and after (b) relocation of earthquakes in Songyuan area

作為相對(duì)定位方法的雙差法，由于地震對(duì)到臺(tái)站之間的路徑幾乎是一致的，這可以有效消除地殼速度結(jié)構(gòu)帶來(lái)的影響，進(jìn)而提高震源定位精度.本研究對(duì)地震重新定位前后P波走時(shí)殘差的分布特征進(jìn)行了統(tǒng)計(jì)(圖4)，從圖上看出定位之前雖然走時(shí)殘差值呈正態(tài)分布特征，但更多的數(shù)據(jù)是集中在-0.5～+0.5，殘差較大(圖4a)；定位后，更多的走時(shí)殘差集中在-0.2～+0.2，說(shuō)明定位后震源位置更接近于實(shí)際地震位置，震源定位精度較高(圖4b).定位后震源位置分布更加集中(圖2d—f)，而從定位后的地震深度直方圖上(圖5)可以看出，地震分布的深度也主要集中在5～14 km，而5～10 km深度的上部地殼區(qū)域地震震源更為集中.圖2a—c表示未定位的震源分布，不難看出震源分布較為凌亂，深度方向上也沒(méi)有規(guī)律可循，但定位后地震成簇分布，表明發(fā)震斷層的主體走向?yàn)楸蔽?南東向(圖2d—f).

圖5 松原震群重定位前(a)后(b)震源深度統(tǒng)計(jì)直方圖Fig.5 Statistical histograms of focal depths before (a) and after (b) relocation of the Songyuan earthquake cluster

1.3 O′Connell-Budiansky理論與裂縫密度、飽和度的計(jì)算

O′Connell-Budiansky理論是一種基于各向同性介質(zhì)條件下，利用地震波速度、泊松比參數(shù)估計(jì)地下巖石裂縫密度和飽和度的方法(O′Connell and Budiansky, 1974; Zhao and Mizuno, 1999).該理論假設(shè)巖石介質(zhì)中變平的圓形裂縫隨機(jī)排列，并充滿了流體.依據(jù)該理論，裂縫密度參數(shù)ε可以定義為

ε=N〈α3〉，

(1)

其中，N為單位體積含有的裂縫數(shù)量，〈α〉表示圓形裂縫的平均半徑.如果假設(shè)Ndry表示干巖石的裂縫數(shù)量，則充滿流體的裂縫數(shù)量為Nsad=N-Ndry，所以飽和度ξ可以定義為

(2)

O′Connell和Budiansky(1974)推導(dǎo)獲得地震波速度(VP,VS)、泊松比(σ)、裂縫密度(ε)和飽和度(ξ)的關(guān)系式為

(3)

(4)

(5)

2 反演計(jì)算與分辨率分析

依據(jù)松原地震區(qū)地震與臺(tái)陣的分布，本次研究模型的選取范圍為東經(jīng)119°～128.5°，北緯43.5°～46°，深度范圍為0～25 km.但模型格點(diǎn)設(shè)置范圍適當(dāng)減小.本次研究采用規(guī)則網(wǎng)格，在經(jīng)度和緯度方向上節(jié)點(diǎn)距離為0.25°，深度方向在0～20 km范圍內(nèi)，每5 km設(shè)置一個(gè)節(jié)點(diǎn).

本次研究采用的初始速度模型為層狀模型，各個(gè)深度層位上的初始速度值參考iasp91全球一維速度模型給出，輸入分別為各層(節(jié)點(diǎn))縱波速度和縱橫波速度比.在深度方向上的5個(gè)節(jié)點(diǎn)(0 km·s-1、5 km·s-1、10 km·s-1、15 km·s-1、20 km·s-1)縱波速度依次為5.8 km·s-1、5.8 km·s-1、5.95 km·s-1、6.15 km·s-1和6.5 km·s-1，波速比設(shè)置為1.73.

地震層析成像結(jié)果的分辨率分析是十分重要的環(huán)節(jié).偏導(dǎo)加權(quán)之和(DWS，Derivative Weight Sum)是衡量反演分辨率和反演精度的有效手段之一(Thurber, 1983).DWS越高，說(shuō)明地震波走時(shí)數(shù)據(jù)對(duì)于該地區(qū)的速度約束越強(qiáng)，反演得到的速度結(jié)構(gòu)精度就越高.所以本研究給出了深度平面上的DWS分布(圖6).從圖中可以看出，白色區(qū)域速度反演的分辨率和精度均較高，這與該區(qū)地震和臺(tái)站分布的密集程度相對(duì)應(yīng).震群附近及研究區(qū)東北部區(qū)域成像精度較高.波速比結(jié)構(gòu)是通過(guò)將反演得到的P波和S波三維速度值直接相除獲得，所以對(duì)于波速比結(jié)構(gòu)的成像分辨率參考了P波和S波速度成像中分辨率較低的S波成像分辨率，即S波成像分辨率在一定程度上代表了波速比結(jié)構(gòu)的分辨率.

圖6 各水平深度層位上P波(a)和S波(b)的DWS(Derivative Weight Sum，偏導(dǎo)加權(quán)之和)分布Fig.6 DWS (Derivative Weight Sum) distribution of P-wave (a) and S-wave (b) at every horizontal depth

鑒于以上DWS分布特征表明，本研究的地震波P波和S波成像的分辨率較好的區(qū)域均分布在地震集中的震群附近，因此本研究針對(duì)P波、S波與波速比結(jié)構(gòu)的深入探討也集中在震群區(qū)域.

3 成像結(jié)果與討論

3.1 介質(zhì)結(jié)構(gòu)異常解釋

圖7為5個(gè)典型深度的P波、S波速度與VP/VS結(jié)構(gòu)圖像及利用反演的速度結(jié)構(gòu)計(jì)算獲得的裂縫密度和飽和度分布圖, 并標(biāo)出主要斷層和震源分布.從圖中可以看出P和S波速度分布具有很好的相似性.總體來(lái)看在近地表低速速度結(jié)構(gòu)呈現(xiàn)近南北向，向深部轉(zhuǎn)變?yōu)楸北睎|向.震區(qū)附近的P波速度較低并延伸至20 km，低S波速異常也延伸至20 km深度附近.VP/VS分布表明該區(qū)8～10 km以淺表現(xiàn)較高異常；在深部區(qū)域，震群附近VP/VS值較低，而裂縫密度和飽和度在10 km以淺區(qū)域也表現(xiàn)為高值.

圖7 不同深度層位上的P波、S波、波速比、裂縫密度和飽和度成像結(jié)果Fig.7 Tomographic results of P waves, S waves, VP/VS, crack density and saturation degree at different depths

為了深入探討介質(zhì)結(jié)構(gòu)變化與地震活動(dòng)性的關(guān)系，本研究沿著震群中心區(qū)域截取了3個(gè)剖面(圖8)，剖面位置見(jiàn)圖1.松原震區(qū)附近共發(fā)育3條隱伏斷層(鄧起東，2007)分別是NE向扶余—肇東斷裂、克山—大安斷裂及NW向查干泡—道字井?dāng)嗔?從截取的剖面上看低P波延伸至20 km深度附近，低S波速度結(jié)構(gòu)向下延伸至10 km深度附近；VP/VS結(jié)構(gòu)在5 km以淺區(qū)域表現(xiàn)為高值，深部VP/VS較低；而高裂縫密度和高飽和度異常亦延伸至10 km附近.從震源分布來(lái)看，地震主要發(fā)生在P波和S波低度、低VP/VS區(qū)域.水平剖面上(圖7)地震分布總體表現(xiàn)為沿北西-南東向，即松原震區(qū)的地震主要發(fā)生在NW-SE向的查干泡—道字井?dāng)嗔焉?；縱向上，地震分布近于垂直，但從垂直扶余—肇東斷裂的剖面上(圖8中BB′剖面)可以看出該斷裂是具有向北東方向傾斜的特征.

圖8 震群附近地震波速度、波速比結(jié)構(gòu)、裂縫密度和飽和度剖面圖(剖面位置見(jiàn)圖1)Fig.8 Vertical tomographic profiles of VP, VS, VP/VS, crack density and saturation degree (locations of profiles are shown in Fig.1)

3.2 流體注入與誘發(fā)地震的可能性

地殼巖石的地震波速度結(jié)構(gòu)、波速比結(jié)構(gòu)與壓力、溫度、巖石組成、裂縫密度和流體含量等密切相關(guān)，實(shí)驗(yàn)研究表明流體對(duì)巖石地震波速度的影響取決于流體類型(水、部分熔融)及巖石孔隙的形狀等(Takei，2002).對(duì)較小的巖石孔隙縱橫比，流體可以降低巖石地震波速度，提高波速比；對(duì)較大的巖石孔隙縱橫比，流體導(dǎo)致地震波速度和波速比均降低(Nugraha et al., 2013; Takei, 2002).

圖8為松原震區(qū)速度結(jié)構(gòu)與震群分布剖面圖，圖中的地震為剖面兩側(cè)2 km范圍內(nèi)的震源位置投影，經(jīng)過(guò)雙差法重定位后的地震分布表明地震位置與速度結(jié)構(gòu)之間有很強(qiáng)的依賴關(guān)系.淺層(<5 km)速度結(jié)構(gòu)表現(xiàn)為低P波和S波速度異常、高波速比異常、且具有高裂縫密度和高飽和度異常，顯然該深度速度異常與淺層流體作用密切相關(guān)，而這一深度發(fā)生的地震很少；在5～15 km深度附近，主體表現(xiàn)為低P波速度、低S波速度和低波速比異常，10 km以下區(qū)域表現(xiàn)為低裂縫密度和低飽和度異常，可能表征了較小的巖石孔隙比，所以流體致使地震波速降低，同時(shí)也導(dǎo)致波速比降低.

大量的監(jiān)測(cè)表明非常規(guī)油氣藏勘探開(kāi)發(fā)過(guò)程中地下深處注水、注氣會(huì)誘發(fā)天然地震(Lei et al., 2013, 2017；Walsh III and Zoback, 2015; 孫小龍等，2018).當(dāng)有大量的流體被壓入地下時(shí)，地下巖石的孔隙壓力會(huì)隨之增大，改變地下原有的應(yīng)力分布情況，進(jìn)而可能導(dǎo)致地質(zhì)力學(xué)變化、斷層復(fù)活、微地震活動(dòng)，甚至?xí)衅茐男缘卣鸢l(fā)生(Lei et al., 2017).2017年1月在四川筠連縣發(fā)生了最大震級(jí)為4.9級(jí)的震群.Lei等(2017) 分別從統(tǒng)計(jì)分析、矩張量反演、震源重新定位、數(shù)值計(jì)算的庫(kù)侖破裂應(yīng)力結(jié)果以及地震發(fā)生指數(shù)等不同方面進(jìn)行研究，認(rèn)為該震群中的大多數(shù)地震是因大量注水再激活預(yù)先存在的斷層而發(fā)生的.2003—2013年間，美國(guó)俄克拉荷馬州中部由于地下深層注入大量廢水影響，地震活動(dòng)明顯增強(qiáng)(Walsh III and Zoback, 2015)，俄克拉荷馬州波尼發(fā)生最強(qiáng)烈的一次地震(5.8級(jí)).Pei等 (2018)對(duì)俄克拉荷馬州地區(qū)的地震構(gòu)造進(jìn)行了各向異性的二維P波層析成像，結(jié)果顯示大多數(shù)中尺度地震都發(fā)生在高速體區(qū)域內(nèi)部或高速體與低速體的邊界附近，表明它們跟不同基底巖石性質(zhì)的地質(zhì)邊界和上地殼中的物質(zhì)性質(zhì)有關(guān).

上述研究表明工業(yè)開(kāi)采期間的注水、注氣會(huì)不同程度的誘發(fā)該區(qū)域天然地震的發(fā)生.而松原震區(qū)同樣處于吉林油田范圍內(nèi).油氣田的開(kāi)采后期或達(dá)到增產(chǎn)的目的，都可能會(huì)采用注水驅(qū)油技術(shù)，這勢(shì)必導(dǎo)致地層裂縫和流體含量的變化，從而使得地層地質(zhì)結(jié)構(gòu)發(fā)生改變.但對(duì)比上述研究結(jié)構(gòu)，我們可以看出松原地區(qū)上地殼的速度結(jié)構(gòu)和震群發(fā)生位置具有較為明顯的區(qū)別.四川筠連震區(qū)和美國(guó)俄克拉荷馬州波尼震群發(fā)生的較大地震均發(fā)生在高速體或高速體與低速體過(guò)渡區(qū)域，而松原震群均發(fā)生在低P波和S波速度異常區(qū)域，只是淺層(<5 km)和中深層(5～15 km)的波速比異常具有較大差異.吉林油田前郭爾多斯地區(qū)油層為泉頭組四段的扶余油層，深度在350～450 m深度左右(王艷麗, 2018)，但松原前郭爾多斯震區(qū)主震群發(fā)生在5～15 km深度，所以推測(cè)該區(qū)地震的發(fā)生可能與油田注水活動(dòng)并無(wú)直接的關(guān)系.

3.3 速度結(jié)構(gòu)不均一性與地震發(fā)生的關(guān)系

前人的研究表明流體與強(qiáng)震發(fā)生有著密切關(guān)系(Zhao et al.，2000，Tian and Liu, 2013；繆淼和朱守彪，2012).Hickman等(1995)研究表明流體廣泛地存在于地球內(nèi)部.日本神戶7.3級(jí)地震震源區(qū)域介質(zhì)表現(xiàn)為低P波和S波速度、高泊松比異常(Zhao et al.，1996)，表明該地震的發(fā)生與殼內(nèi)流體的作用息息相關(guān).流體可能來(lái)自深部俯沖的太平洋板塊脫水作用.流體作用導(dǎo)致了巖石有效應(yīng)力的下降和孔隙壓力的增加，使斷層弱化而觸發(fā)地震(田有等，2007).上述的研究多基于大洋和大陸俯沖帶地區(qū)，如太平洋俯沖帶在日本島弧下面的俯沖作用.由于俯沖的太平洋板塊富含流體，帶入到上部地幔后釋放到地幔楔內(nèi)部，使得島弧地殼富含流體；但大陸內(nèi)部強(qiáng)震的形成機(jī)制可能有不同認(rèn)識(shí).

地殼結(jié)構(gòu)的橫向不均勻性與地震分布密切相關(guān)(田有等，2007；田有，2008；曲中黨等，2018)，尤其是大陸板塊內(nèi)部區(qū)域絕大多數(shù)大地震(>6.0級(jí))發(fā)生與高地震波速度異常密切相關(guān)(Huang and Zhao，2004；田有等，2007)；近來(lái)對(duì)華北1960年以來(lái)發(fā)生的大于6級(jí)的地震研究發(fā)現(xiàn)，這些地震多數(shù)也都分布在高P波、S波速度，低泊松比區(qū)域，或者發(fā)生在高速和低速異常體的分界區(qū)域，表明高速區(qū)多屬地殼脆性介質(zhì)，應(yīng)力容易集中，因而導(dǎo)致大地震；相反，低速度區(qū)域則可能代表了破裂程度高、富含流體或溫度較高的區(qū)域，因而傾向于發(fā)生無(wú)震變形或小震群(田有等，2007；田有，2008).這與松原震區(qū)非常相似，該區(qū)地震均發(fā)生在低P波、S波與高波速比區(qū)域，而且最近幾年尚未發(fā)生過(guò)7級(jí)左右的大地震，表明該區(qū)的低速物質(zhì)不容易積累產(chǎn)生大地震的應(yīng)力，所以推測(cè)松原震區(qū)發(fā)生較大破壞性地震的可能性較低.松原前郭震區(qū)震源機(jī)制研究顯示該區(qū)4.5級(jí)以上的地震具有逆沖斷層兼走滑性質(zhì)(吳微微等，2014)，而震群又發(fā)生在NW走向的斷層上，即與太平洋北西向俯沖方向一致.

4 結(jié)論

(1)地震重新精確定位結(jié)果表明松原地區(qū)地震主要發(fā)生在NW-SE的查干泡—道字井?dāng)嗔焉?，該斷裂傾向?yàn)楸睎|向，近垂直斷裂.

(2)松原震群區(qū)主要表現(xiàn)為低P波、低S波、低波速比結(jié)構(gòu)，中國(guó)大陸內(nèi)部，尤其是從華北地區(qū)大地震(>6.0級(jí))的研究來(lái)看，大地震主要發(fā)生在高P波和S波速度異常，低泊松比異常區(qū)，所以推測(cè)松原地震區(qū)發(fā)生更大地震的可能性較低.

(3)松原震群區(qū)淺地表區(qū)域(<5 km)整體表現(xiàn)為低P波、低S波、高波速比、高裂縫密度和飽和度異常，淺層(<5 km)高波速比結(jié)構(gòu)可能與流體有關(guān)，但該區(qū)產(chǎn)油的扶余油層較淺(350～450 m深度)，推測(cè)地震發(fā)生可能與油田驅(qū)油注水并無(wú)直接相關(guān).從區(qū)域構(gòu)造出發(fā)表明松原地震發(fā)震構(gòu)造與東北地區(qū)的區(qū)域構(gòu)造應(yīng)力場(chǎng)密切相關(guān),即可能與太平洋板塊俯沖遠(yuǎn)程效應(yīng)有關(guān).

致謝感謝張海江教授提供的雙差地震層析成像計(jì)算程序(TOMODD)，感謝中國(guó)國(guó)家地震局和吉林省地震局提供的地震波走時(shí)數(shù)據(jù)，感謝兩位評(píng)審專家提出的詳細(xì)的、寶貴的修改意見(jiàn)和建議.本文中圖件均用GMT軟件繪制(Wessel and Smith，1998).