張致偉，程萬正，梁明劍，王曉山，龍 鋒，許 艷，陳文康，王世元

1 四川省地震局，成都 610041

2 河北省地震局，石家莊 050021

3 自貢市防震減災局，自貢 643000

四川自貢—隆昌地區(qū)注水誘發(fā)地震研究

張致偉1，程萬正1，梁明劍1，王曉山2，龍 鋒1，許 艷3，陳文康3，王世元1

1 四川省地震局，成都 610041

2 河北省地震局，石家莊 050021

3 自貢市防震減災局，自貢 643000

基于2007年1月—2010年8月四川自貢地方數字測震臺網和流動臺站記錄的地震觀測資料，并結合注水井（家33井）的加壓注水數據，研究了自貢—隆昌地區(qū)三個叢集區(qū)域（A、B、C）地震活動與家33井注水壓力及注水量的相關性，同時從注水區(qū)域及鄰區(qū)局部應力狀態(tài)呈現的差異研究其發(fā)震機理.結果表明：家33井加壓注水對其所在的B區(qū)地震活動具有明顯的調控作用，負壓注水階段，注水量對地震活動的影響不大，地震活動處于較弱狀態(tài)；加壓注水階段，在壓力持續(xù)升高的背景下，地震活動明顯增強，其頻次、強度與注水量呈現較好的相關性.研究區(qū)三次4級地震震源較淺、且均為逆沖型，震源機制節(jié)面與其震中附近斷層走向一致.基于大量小震震源機制解，統(tǒng)計獲得B區(qū)震源機制節(jié)面走向和P軸方位較A區(qū)相對離散；與川東南地區(qū)已有應力場結果相比，反演獲得的B區(qū)主壓應力方位較A區(qū)擾動偏大，可見在統(tǒng)一應力場的作用下，A、B區(qū)域小震震源機制解節(jié)面走向、P軸方位及反演得到的局部區(qū)域應力場呈現出的差異，與家33井加壓注水有一定的關系.

誘發(fā)地震，加壓注水，震源機制解，應力場，自貢—隆昌地區(qū)

1 引 言

隨著油氣工業(yè)的迅猛發(fā)展，國內外出現了多起因油氣田開采及高壓注水誘發(fā)地震活動的例子.20世紀20年代以來，美國德克薩斯州南部的鵝溪（Goose Creek）油田［1］、加利福尼亞州的威爾明頓（Wilmington）油田［2］、科羅拉多州丹佛東北地區(qū)的洛基山深井注入化學廢液［3］、中亞地區(qū)的烏茲別克葛茲力（Gazili）氣田及俄羅斯羅馬?；Z（Romashkino）油田［4］等開采及注水都不同程度的引起了鄰近地區(qū)的地表破裂和地震活動；而國內在70年代以來，也相繼觀測到華北任丘油田［5］、山東勝利油田［6－8］、重慶榮昌地區(qū)采氣注水［9］及四川長寧鹽礦井注水［10］等工業(yè)活動誘發(fā)地震的現象.高壓注水誘發(fā)地震源于孔隙壓力增加致使巖石的剪切強度降低，外加水的潤滑作用，從而在周圍高應力的作用下引起的斷層滑動［3，11－12］.

四川作為全國主要的油氣田和井礦鹽產區(qū)之一，存在多處采鹽、采氣及廢水回注井.其中位于自貢市大安區(qū)牛佛鎮(zhèn)與隆昌縣黃家場鎮(zhèn)交界的家33井目前為天然氣采空廢井，為避免地面環(huán)境污染，氣礦自貢作業(yè)區(qū)將其附近其他天然氣生產井所產生的工業(yè)廢水回注到家33井中.2009年1月上旬，該井出現容腔飽和，隨后開始人為加壓注水，家33井周邊地震活動也相應出現了明顯增強的異?，F象，并于2009年2月16日、5月22日發(fā)生ML4.4、ML4.2級地震，這一現象引起了當地政府的高度重視，隨后地震部門增設流動數字地震臺進行監(jiān)測，并開展了相應的研究工作.

以四川自貢—隆昌地區(qū)（29.0°N—29.8°N，104.4°E—105.4°E）為研究區(qū)域，本文基于2007年1月—2010年8月自貢地方數字測震臺網和流動臺站記錄的地震觀測資料，并結合家33井加壓注水數據，試圖分析注水區(qū)域及鄰區(qū)地震活動與家33井注水壓力、注水量的相關性，以及在統(tǒng)一區(qū)域應力場作用下，加壓注水對周邊局部區(qū)域應力場的擾動影響.

2 地質構造背景、臺站分布概況

2.1 區(qū)域地質構造及家33井位置

自貢—隆昌地區(qū)位于四川盆地東南腹地，盆地基底為前震旦系結晶基底，一套強磁性的深變質巖、混合巖及中、基性火山雜巖建造，蓋層為比較穩(wěn)定的海陸相沉積建造［13］（圖1）.四川盆地的基本構造格局成型于晚三疊世，印支運動導致了松潘—甘孜海槽的閉合、龍門山構造帶的推覆，以及山前前陸盆地的形成［14－15］.新生代以來，伴隨著印—亞板塊碰撞和青藏高原的強烈隆升，四川盆地西界的龍門山構造帶主要表現為逆沖－右旋走滑運動特性，盆地的構造變形樣式則受到了構造運動的疊加影響，形成盆地內斷褶構造格局，產生了許多油氣構造及鹽拱構造.其地表構造樣式則呈一系列平緩、低陡的背斜與向斜相間排列，地表斷裂較少，常與背斜構造相伴生，為斷裂彎曲背斜或斷層擴展背斜構造成因［16－17］，地下蓋層中常存在地腹構造.圖1還標注了家33注水井及三次4級地震的位置.

家33井位于自貢市大安區(qū)牛佛鎮(zhèn)與隆昌縣黃家場鎮(zhèn)交界，其所在區(qū)域的主要構造為牛佛渡—普臺寺背斜（圖2），長約20km，軸向N80°E.從黃家場氣田的二疊系下統(tǒng)（P1）頂界埋深等值線圖來看，蓋層中褶皺構造的軸向呈NE－NEE，與地表出露的構造軸向基本一致.牛佛渡斷層正好呈NW走向展布于P1頂界高點向低點傾斜過渡的部位，傾向NE，主要表現為擠壓逆斷層性質.家33井的注水深度在2500m左右，層位為二疊系茅口組，該產氣層位的斷層和節(jié)理發(fā)育，利于回注廢水.

2.2 臺站分布概況

圖3 研究區(qū)域臺站分布圖Fig.3 The distribution map of stations in study area

本研究使用了兩個臺網（圖3）的地震觀測資料.其中自貢地方數字測震臺網（灰色三角）運行于1994年，由椅子灣（YZW）、石佛寺（SFS）、碾子山（NZS）、三元井（SYJ）和趙家坡（ZJP）5個子臺組成.2009年1月，該臺網的東側富順—隆昌交界地區(qū)出現了小震活動明顯增強的異?，F象，隨后四川省地震局在小震增強區(qū)域周圍布設了5套流動測震臺（黑色三角），包括何市臺（HES）、牛佛臺（NIF）、王達山（WDS）、黃家鎮(zhèn)（HUJ）和沱江臺（TUJ），該流動臺網于2009年7月正式開始運行.這10個子臺均勻地展布于小震增強區(qū)域周圍，平均臺距為10km，監(jiān)測能力震級下限可達ML0.1級.

3 地震活動特征

圖4給出了自貢—隆昌地區(qū)2007年1月—2010年8月2706次地震的精定位結果及家33井的位置.地震主要集中分布在自貢（A區(qū)）、富順—隆昌交界（B區(qū)）及隆昌（C區(qū)）三個叢集區(qū)域，其中A區(qū)小震呈長軸NNE向展布，B、C兩區(qū)小震則沿長軸NW向展布.三次ML4級地震時空分布顯示，2008年10月10日自貢ML4.6級地震發(fā)生在A區(qū)，震中位于家33井西側約20km，該地震發(fā)生在汶川8.0級地震之后，家33井容腔飽和加壓注水之前；2009年2月16日、5月22日富順ML4.4、ML4.2級地震發(fā)生在家33井所在的B區(qū)，且發(fā)生于家33井容腔飽和后的加壓注水時段.

3.1 地震活動與注水壓力的關系

圖4 研究區(qū)域中小地震及注水井位置分布圖（2007－01—2010－08）Fig.4 The distribution map of water injection well and moderate－small earthquakes from January 2007to August 2010in study area

為了進一步了解上述三個叢集區(qū)域地震活動與家33井注水壓力的關系，本文分別給出了A、B、C三個區(qū)域地震活動頻次與家33井加壓數據的關系（圖5）.其中A區(qū)地震活動增強時段為汶川8.0級地震之后、家33井加壓之前，加壓之后地震活動頻次顯示正常水平（圖5a）.分析認為2008年5月12日汶川8.0級地震的發(fā)生在一定程度上改變了四川盆地地下應力場環(huán)境，致使盆地東南緣出現地震活動增強現象.家33井所在的B區(qū)地震活動明顯增強發(fā)生在2009年1月注水井加壓之后（圖5b），地震活動頻次與家33井注水壓力有很好的相關性，可見家33井注水壓力對B區(qū)地震活動具有明顯的調控作用；C區(qū)地震活動頻次在家33井加壓之后也有所增強（圖5c），但因距注水井相對B區(qū)較遠，地震活動沒有B區(qū)明顯.

研究發(fā)現：B、C兩區(qū)小震活動增強均為本區(qū)初次現象，而且均以微小地震為主，家33井注水壓力對這兩個區(qū)域的地震活動具有不同程度的調控作用.其中B區(qū)小震震中靠近注水井，持續(xù)加壓導致地下巖層節(jié)理、裂隙擴張，外加水的潤滑作用產生新的微破裂，從而使得該地區(qū)中小地震活動增強.微震活動與水流通道的關系也是比較復雜的，并非想象中的沿一個平面分布，而與注水所致裂隙有相當密切的關系［18］.鑒于C區(qū)小震活動的時空分布特征，認為該區(qū)小震活動可能是家33井加壓致使B、C區(qū)域之間地下裂隙貫通，流體向C區(qū)深部延伸產生微裂隙的結果.反之，也可根據地震活動深部空間分布來推斷流體在地下的流向及裂隙的延伸方向.

圖5 各區(qū)地震頻次N（柱狀圖）與家33井注水壓力曲線（灰色曲線）圖（a）區(qū)域A；（b）區(qū)域B；（c）區(qū)域C.Fig.5 The earthquake frequency of each region and the curve graph of injection pressure of Jia 33well（a）Region A；（b）Region B；（c）Region C.

3.2 B區(qū)地震活動與注水量的關系

家33井注水壓力與上述三個區(qū)域地震活動的關系表明：家33井注水壓力對B區(qū)地震活動具有明顯調控作用.下面將詳細討論B區(qū)地震活動在不同加壓階段與注水量的關系.圖6包含了注水壓力（MPa）、注水量（m3）、M－T和N－T圖.我們將其分為負壓（圖中1段）和加壓（圖中2段）兩個階段分別進行討論.

第一階段為負壓注水階段，即2007年1月—2008年12月，此階段家33井沒有加壓，但注水量不斷變化，2008年3月以前，注水量一直控制在315m3內，2008年4、5月份注水量有所增加，最高值達到了630m3，在經歷了2008年5月份短期高值后，注水量于6月份迅速下降為約200m3，持續(xù)兩個月之后于8月份又急劇上升到500m3左右，且一直持續(xù)到了年底.從B區(qū)的震級、頻次時間圖來看，該時段地震活動一直處于較弱狀態(tài).顯然在負壓注水階段，注水量對地震活動的影響不大.

第二階段為加壓注水階段，2009年1月上旬，家33井出現容腔飽和，隨后開始加壓注水.容腔飽和致使家33井注水量總體呈現下降趨勢，但仍然存在有起伏的波動.2009年1月13日開始加壓，B區(qū)地震活動相應出現了顯著增強現象，2月10日壓力迅速升高到2.7MPa后于2月16日在家33井周邊發(fā)生ML4.4級地震，5月上旬至中旬，壓力持續(xù)為2.82～2.96MPa，又于5月22日發(fā)生ML4.2級地震.6月注水壓力已經達到4MPa以上，小震月頻度也達到了160多次.這一現象引起了地震部門的高度重視，隨后與企業(yè)進行溝通，該作業(yè)井于7月份將注水量下降到了初始注水量（315m3）以內，隨后的兩個月內月頻次及3級以上地震強度均有所下降.2009年1月—2010年7月期間，在注水壓力持續(xù)升高的背景下，地震活動強度、頻次與注水量呈現出較好的對應關系，即注水量下降—升高—下降與地震活動減弱—增強—減弱的顯著關系.任丘油田845號井加壓注水與地震活動［19］同樣呈現與本研究相似的形態(tài).可見，加壓注水階段，地震活動與注水量呈現了較好的相關性.

4 震源機制解及應力場特征

4.1 初始速度模型及求解方法

基于趙珠等［20］獲得的四川東部盆地地殼上地幔平均速度模型，并結合近些年該區(qū)域的地殼速度研究結果［21－22］，最終采用表1作為本研究的初始速度模型.

圖6 區(qū)域B地震活動與家33井加壓注水的曲線圖Fig.6 The seismicity of region B and the curve graph of injection pressure of Jia 33well

表1 研究區(qū)域地殼P波平均速度模型Table 1 The average model of P wave velocity in study area

本文針對不同震級的地震采用不同的震源機制求取方法.對于研究區(qū)內的中等地震，采用最近十多年國際上不斷發(fā)展和完善的CAP波形反演方法［23－28］，其綜合利用了近震中體波和面波信息，將寬頻帶數字地震波形記錄分解為體波Pnl和面波兩部分，計算并搜索理論地震波形與真實地震波形之間擬合誤差函數最小的機制解.該方法具有計算臺站數量少、反演結果對地殼速度結構模型及橫向變化的依賴性相對較小等優(yōu)點［25－28］，在獲得震源機制解的同時還能給出最佳擬合震源深度.

對于小震震源機制解，梁尚鴻等［29］提出了一種利用區(qū)域地震臺網地震波的ˉP、ˉS垂直分量振幅比資料的求解方法，以層狀介質中一點源位錯震源模型，采用廣義透射系數的快速算法和理論地震圖擬合直達波最大振幅比來求取小震震源參數.胡新亮等［30］就該方法測定小震震源機制解的可靠性進行了分析，認為小震的發(fā)生具有很強的隨機性，不便逐一進行具體分析，統(tǒng)計研究應該更加客觀和恰當.基于上述思路，本文針對研究區(qū)獲得的大量小震震源機制解，采取了統(tǒng)計的方法開展研究工作.

4.2 研究區(qū)域3次4級地震震源機制解分析

選取四川區(qū)域地震臺網震中距在250km以內臺站的寬頻帶波形記錄，采用CAP方法反演獲得了自貢—隆昌地區(qū)發(fā)生的三次ML4.0級以上地震震源機制解，詳細參數列于表2.圖7給出了上述三次4級地震震源機制反演中的深度擬合誤差分布.可見，這三次地震的震源深度分別在3、5km和4km處具有最小的反演誤差，表明這三次地震均屬淺源地震.事件b、c的震源機制比較相似，且有別于事件a，事件a的震源機制反演結果不穩(wěn)定，隨著不同的深度變化較大，而事件b、c則除了大于14km的震源機制解稍有不同外，其他深度解的圖案幾乎一致.

圖8給出了自貢—隆昌地區(qū)三次ML4.0級以上地震震源機制解分布.三次地震均表現為逆沖型錯動類型，震源機制參數顯示，事件a的節(jié)面1呈NNE走向，與震中附近的長巖斷層走向一致，震源機制P軸方位為近NWW向.事件b、c的震源機制參數比較一致，節(jié)面1呈NW向，與其震中附近的牛佛渡斷層走向一致，震源機制P軸方位為近NEE向.已有研究結果表明自貢—隆昌所在的川東南地區(qū)區(qū)域應力場方位為近NW向［31－32］，距家33井西約20km的事件a震源機制P軸方位與川東南地

區(qū)應力場方向接近，然而發(fā)生在家33井附近的事件b、c震源機制解P軸方位卻與區(qū)域應力場存在一定的偏差.趙根模等［19］研究發(fā)現任丘油田局部區(qū)域幾次4級地震的等震線形態(tài)同樣存在差異，反映了誘發(fā)地震不僅受區(qū)域應力場影響，還要受到注水井局部原生構造的影響.

表2 研究區(qū)域ML4.0級以上地震的震源機制解參數Table 2 The focal mechanism solution parameters of ML≥4.0earthquakes in study area

這里需要說明的是，本文將中等地震震源機制P軸方位與已有區(qū)域應力場比較分析，旨在觀測現今地震活動所呈現力學機制的一致或變化，而不認為單個或幾個中等地震震源機制解的P軸就代表地塊主壓應力方向.Julian等［33］認為，理想化的雙力偶力系能表征均勻各向同性介質中的平面剪切斷裂，但不穩(wěn)定液體流，環(huán)狀構造上的剪切破裂和張性斷裂等過程可能產生與此機制不同的地震機制，在地熱和火山環(huán)境中是特別容易出現的.本研究的局部區(qū)域存在人為造成的不穩(wěn)定液體流，對相距20km的兩區(qū)域中等地震震源機制P軸方向與區(qū)域應力場進行這種比較，目的在于嘗試分析家33井加壓注水對周邊不同局部區(qū)域應力場的擾動影響.

4.3 A、B區(qū)域小震震源機制及應力場分析

利用2007年1月—2010年8月自貢地方數字測震臺網和四川流動臺站記錄的地震波形資料，采用振幅比方法測定了該研究區(qū)域小震震源機制解.其中區(qū)域A的地震活動較弱，震級下限選為ML1.5，區(qū)域B的地震活動較強，震級下限選為ML2.5.最終獲得了72次小震的震源機制，圖9給出了集中分布在區(qū)域A、B的小震震源機制解及其局部應力場反演結果，在求得的震源機制解中，A區(qū)絕大多數小震發(fā)生在2007和2008兩年，而B區(qū)小震則發(fā)生在2009年家33井容腔飽和后的加壓注水階段.

基于大量小震震源機制解，我們按10°間隔進行統(tǒng)計并計算歸一頻數，分別給出了區(qū)域A、B的小震震源機制解的節(jié)面和力軸參數玫瑰圖（圖10）.因為小震無法區(qū)分斷層面和輔助面，故在統(tǒng)計中同等看待，合在一起進行分析.

A區(qū)參數統(tǒng)計結果顯示（圖10a）：小震節(jié)面（strike）主要集中在NW和NEE兩個優(yōu)勢方向，其次為近NS走向，節(jié)面傾角（dip）近直立，主要分布在60°—90°范圍內，其中80°—90°所占比例最大，根據滑動角（slip）分析震源力學作用方式，主要呈現為左旋和右旋走向滑動，其次為正傾型錯動類型，逆沖型地震相對較少，但依然獨立存在，震源機制解P軸方位（paz）一致性較好，以近NWW－SEE為優(yōu)勢方向，傾角（pdip）分布在0°—60°范圍內，其中30°—60°所占比例較大，T軸的優(yōu)勢方向（taz）主要集中在NNE－SSW向，傾角（tdip）在30°—50°范圍內較多.用同樣的方法統(tǒng)計得到B區(qū)小震震源機制解的參數玫瑰圖（見圖10b），節(jié)面走向分布比較分散，各個角度都有，小有優(yōu)勢的節(jié)面走向分別為NNE和EW向，節(jié)面傾角分布在50°—80°，主要呈現為走滑和逆沖型地震，正斷型地震相對較少，小震震源機制P軸方位相對比較凌亂，但優(yōu)勢方位為NE和NW向，P軸傾角以近水平的低傾角（＜10°）占優(yōu)勢，T軸優(yōu)勢方向為NW－SE向，也表現出較低的傾角.從統(tǒng)計結果來看，B區(qū)小震震源機制節(jié)面及P軸方位較A區(qū)結果相對比較離散，分析認為可能是由于家33井加壓注水導致地下巖層節(jié)理、裂隙擴張所致.

圖9 區(qū)域A（a）、B（b）中小地震震源機制和應力場應力主軸的投影Fig.9 The projection of focal mechanism and stress axis of moderate－small earthquakes in region A（a）and region B（b）

圖10 區(qū)域A（a）、B（b）小震震源機制解的節(jié)面和力軸參數玫瑰圖Fig.10 Rose diagrams of nodal plane and stress axis parameters of small earthquakes focal mechanism solutions in region A（a）and region B（b）

本文采用FMSI方法［34］，基于A、B區(qū)域的小震震源機制解分別反演獲得了各區(qū)應力場，結果如圖9所示，并將參數列于表3.A區(qū)應力場的結果顯示：最大主壓應力軸σ1的方位為254°、傾角34°，σ2的方位為157°、傾角10°，最小主壓應力軸σ3的方位為52°、傾角54°；然而B區(qū)應力場的最大主壓應力軸σ1的方位為27°、傾角55°，σ2的方位為285°、傾角9°，最小主壓應力軸σ3的方位為189°、傾角34°.與川東南地區(qū)已有應力場研究結果［32－33］相比，B區(qū)應力場主壓應力方位較A區(qū)的擾動更大，表明在統(tǒng)一應力場的作用下，自貢—隆昌地區(qū)注水區(qū)域及鄰區(qū)局部應力場的非一致性.

表3 A、B區(qū)域應力場反演結果Table 3 Inversion results of stress field of region A and B

大區(qū)域板塊或地塊的應力場是基本穩(wěn)定的，而且具有統(tǒng)一的應力場作用，但這并不排除地塊內局部小區(qū)域應力場出現的擾動現象.Seeber［35］認為地震可作為應力變化的標志，雖然小地震對總的變形貢獻很小，但它們在空間、時間和運動上的分布對應力變化非常敏感.影響地殼力學狀態(tài)的各種現象，包括諸如斷層蠕動或巖漿注入的自然現象，水庫蓄水或油田注水的大型工程活動.本研究區(qū)域介質因注水已發(fā)生了變化，震源力學機制同樣也反映出了這種變化.

5 結論與討論

結合注水井的加壓注水數據，分析了四川自貢—隆昌地區(qū)三個叢集區(qū)域地震活動與家33井注水壓力及注水量的相關性，并且基于中小地震震源機制解，進一步探討了加壓注水對周邊局部區(qū)域應力場的擾動影響.獲得的主要認識如下：

（1）距家33井西側約20km的A區(qū)地震活動增強主要受汶川8.0級地震的影響，與家33井加壓注水無關；而家33井注水壓力卻對B、C區(qū)域的地震活動具有不同程度的調控作用，其中B區(qū)小震震中靠近注水井，調控作用較為顯著，注水井持續(xù)加壓導致地下巖層節(jié)理、裂隙擴張，外加水的潤滑作用，從而產生更多微破裂，致使小震增強；位于家33井東南方向約15km的C區(qū)小震活動可能源于加壓致使B、C區(qū)域之間地下裂隙貫通，流體向C區(qū)深部延伸產生微裂隙的結果.B區(qū)地震活動在不同加壓時段與注水量的關系表明：負壓注水階段，注水量對地震活動的影響不大，地震活動處于較弱狀態(tài)；加壓注水階段，在壓力持續(xù)升高的背景下，地震活動頻次、強度與注水量呈現較好的相關性.

（2）研究區(qū)三次4級地震震源較淺、且均為逆沖型，震源機制節(jié)面與各自震中附近斷層走向一致.發(fā)生在A區(qū)的地震事件a震源機制P軸方位與川東南地區(qū)已有應力場結果一致，而家33井所在區(qū)域的事件b、c震源機制解P軸方位卻存在一定的偏差；基于大量小震震源機制解，統(tǒng)計獲得B區(qū)小震震源機制節(jié)面及P軸方位較A區(qū)結果離散，可能源于家33井加壓注水導致地下巖層節(jié)理、裂隙擴張所致，與川東南地區(qū)已有應力場結果相比，反演得到B區(qū)應力場主壓應力方位較A區(qū)的擾動更大，表明在統(tǒng)一應力場的作用下，自貢—隆昌地區(qū)注水區(qū)域及鄰區(qū)局部應力場的非一致性.

注水與地震的關系是相當復雜的，分析認為自貢—隆昌注水區(qū)小震活動增強及局部應力場擾動均與家33井的加壓注水有關，加壓注水是誘發(fā)地震的主要因素，同時除了受區(qū)域應力場作用以外，還與局部淺層構造、注水方式及注水量等因素有關.本文僅對自貢—隆昌地區(qū)地震活動與家33井注水壓力及注水量的關系，震源機制及局部小區(qū)域應力場擾動特征做了初步探討，在今后的研究中有必要從物理意義上開展更進一步的深入研究.

致 謝 四川省地震局監(jiān)測研究所、自貢市防震減災局為本研究提供了地震觀測資料；CAP波形反演程序來源于美國圣路易斯大學的朱露培博士，振幅比方法采用河北省地震局刁桂苓研究員提供的程序；研究中還得到了四川省地震局杜方研究員、張永久副研究員等的大力支持和幫助；兩位評審專家提出了非常中肯、有價值的修改意見.在此一并致謝！

（References）

［1］ Pratt W E，Johnson D W.Local subsidence of the Goose Creek oil field.Journal of Geology，1926，34（1）：.577－590.

［2］ Yerkes R F，Castle R O.Seismicity and faulting attributable to fluid extraction.Engineering Geology，1976，10（24）：151－167.

［3］ Healy J H，Rubey W W，Griggs D T，et al.The Denver earthquakes.Science，1968，161（3848）：1301－1310.

［4］ 刁守中，周煥鵬，郭愛香.山東角07井注（漏）水誘發(fā)地震—震源參數與地震序列特征.地震地質，1987，9（3）：84－89.Diao S Z，Zhou H P，Guo A X.Induced earthquakes by injecting and leaking water in Jiao 7oil well，Shandong Province—characteristics of the focal parameters and earthquake sequences.Seismology and Geology（in Chinese），1987，9（3）：84－89.

［5］ 刁桂苓，張德元，趙根模.任邱油田注水地震的初步研究.西北.地震學報，1982，4（3）：89－95.Diao G L，Zhang D Y，Zhao G M.Preliminary study on water injection earthquake of Renqiu oil－fleld.Northwestern Seismological Journal（in Chinese），1982，4（3）：89－95.

［6］ 刁守中，周煥鵬，郭愛香.山東角07井注（漏）水誘發(fā)地震—震源參數與地震序列特征.地震地質，1987，9（3）：84－89.Diao S Z，Zhou H P，Guo A X.Induced earthquakes by injecting and leaking water in Jiao 7oil well，Shandong Province—characteristics of the focal parameters and earthquake sequences.Seismology and Geology（in Chinese），1987，9（3）：84－89.

［7］ 刁守中.山東角07井注水誘發(fā)地震序列的時間分布特征.地震，1989，（1）：27－31.Diao S Z.Time distribution feature of seismic sequence induced by water flooding in Jiao 07well in Shandong Province.Earthquake（in Chinese），1989，（1）：27－31.

［8］ 刁守中，蔣海昆，徐學炎.山東勝利油田角07井注水地震序列的演化及其機制.地震學報，1990，12（4）：399－406.Diao S Z，Jiang H K，Xu X Y.The evolution and mechanism of earthquake sequences induced by water injection in Jiao 7 oil well of Shengli oil field，Shandong Province.Acta Seismologica Sinica（in Chinese），1990，12（4）：399－406.

［9］ 朱麗霞，黃世源，魏紅梅.榮昌地區(qū)注水地震研究.大地測量與地球動力學，2007，27（6）：86－90.Zhu L X，Huang S Y，Wei H M.On fluid－injection induced earthquake in Rongchang area.Journal of Geodesy and Geodynamics（in Chinese），2007，27（6）：86－90.

［10］ 阮祥，程萬正，張永久等.四川長寧鹽礦井注水誘發(fā)地震研究，中國地震，2008，24（3）：226－234.Ruan X，Cheng W Z，Zhang Y J，et al.Research of the earthquake induced by water injections in salt mines in Changning，Sichuan.Earthquake Research in China（in Chinese），2008，24（3）：226－234.

［11］ Segall P.Earthquakes triggered by fluid extraction.Geology，1989，17（10）：942－946.

［12］ Raleigh C B，Healy J H，Bredehoeft J D.Faulting and crustal stress at Rangely，Colorado.／／Heard H C，Borg I Y，Carter N L，et al eds.Flow and Fracture of Rocks（Griggs Volume）：American Geophysical Union Geophysical Monograph，1972，16：275－284.

［13］ 許志琴，候立偉，王宗秀等.中國松潘—甘孜造山帶的造山過程.北京：地質出版社，1992：1－190.Xu Z Q，Hou L W，Wang Z X，et al.The Orogenic Process of the Songpan－Garze Orogenic Belt（in Chinese）.Beijing：Beijing Publishing House，1992：1－190.

［14］ 周榮軍，唐榮昌，雷建成.四川盆地潛在震源區(qū)的細致劃分.四川地震，2005，（3）：1－6.Zhou R J，Tang R C，Lei J C.Attentive seismic hazard mapping for Sichuan Basin.Earthquake Research in Sichuan（in Chinese），2005，（3）：1－6.

［15］ 宋鴻彪，羅志立.四川盆地基底及深部地質結構研究的進展.地學前緣，1995，2（3－4）：231－237.Song H B，Luo Z L.The study of the basement and deep geological structures of Sichuan basin，China.Earth Science Frontiers（in Chinese），1995，2（3－4）：231－237.

［16］ 錢洪.四川斷裂活動的區(qū)域性差異及其與區(qū)域地殼運動的關系.地震研究，1995，18（1）：49－55.Qian H.Regional variation of fault movement in Sichuan and its relation to regional crustal movement.Journal of Seismological Research（in Chinese），1995，18（1）：49－55.

［17］ 李勇，曾允孚.試論龍門山逆沖推覆作用的沉積響應——以成都盆地為例.礦物巖石，1994，14（1）：58－66.Li Y，Zeng Y F.On the sedimentary response to thrusting of Longmenshan thrust belt in Chengdu Basin.Journal of Mineralogy and Petrology（in Chinese），1994，14（1）：58－66.

［18］ 張寶紅，邱澤華.關于注水地震研究的幾個問題.現代地質，1994，8（3）：329－333.Zhang B H，Qiu Z H.On study of fluid injection earthquake.Geoscience（in Chinese），1994，8（3）：329－333.

［19］ 趙根模，楊港生.注水誘發(fā)地震的震源應力場和尾波衰減特征.地震地質，1990，12（4）：303－310.Zhao G M，Yang G S.Focal stress field by water injectioninduced earthquakes and characteristics of coda wave attenuation.Seismology and Geology（in Chinese），1990，12（4）：303－310.

［20］ 趙珠，張潤生.四川地區(qū)地殼上地幔速度結構的初步研究.地震學報，1987，9（2）：154－166.Zhao Z，Zhang R S.Primary study of crustal and upper mantle velocity structure of Sichuan Province.Acta SeismologicaSinica（in Chinese），1987，9（2）：154－166.

［21］ 王椿鏞，Mooney W D，王溪莉.川滇地區(qū)地殼上地幔三維速度結構研究.地震學報，2002，24（1）：1－16.Wang C Y，Mooney W D，Wang X L.Study on 3－D velocity structure of crust and upper mantle in Sichuan－Yunnan region，China.Acta Seismologica Sinica（in Chinese），2002，24（1）：1－16.

［22］ 馬宏生，張國民，聞學澤等.川滇地區(qū)三維P波速度結構反演與構造分析.地球科學（中國地質大學學報），2008，33（5）：591－602.Ma H S，Zhang G M，Wen X Z，et al.3－D P wave velocity structure tomographic inversion and its tectonic interpretation in Southwest China.Earth Science（Journal of China University of Geosciences）（in Chinese），2008，33（5）：591－602.

［23］ Zhao L S，Helmberger D V.Source estimation from broadband regional seismograms.Bull.Seism.Soc.Amer.，1994，84（1）：91－104.

［24］ Zhu L P，Helmberger D V.Advancement in source estimation techniques using broadband regional seismograms.Bull.Seism.Soc.Amer.，1996，86（5）：1634－1641.

［25］ 鄭勇，馬宏生，呂堅等.汶川地震強余震（MS≥5.6）的震源機制解及其與發(fā)震構造的關系.中國科學D輯，2009，39（4）：413－426.Zheng Y，Ma H S，LüJ，et al.Source machanism of strong aftershocks（MS≥5.6）of the 2008／05／12Wenchuan earthquake and the implication for seismotectonics.ScienceinChina（Series D），2009，52（6）：739－753.

［26］ 呂堅，鄭勇，倪四道等.2005年11月26日九江—瑞昌MS5.7、MS4.8地震的震源機制解與發(fā)震構造研究.地球物理學報，2008，51（1）：158－164.LüJ，Zheng Y，Ni S D，et al.Focal mechanisms and seismogenic structures of the MS5.7and MS4.8Jiujiang－Ruichang earthquakes of Nov.26，2005.Chinese J.Geophys.（in Chinese），2008，51（1）：158－164.

［27］ 龍鋒，張永久，聞學澤等.2008年8月30日攀枝花—會理6.1級地震序列ML≥4.0事件的震源機制解.地球物理學報，2010，53（12）：2852－2860.Long F，Zhang Y J，Wen X Z，et al.Focal mechanism solutions of ML≥4.0events in the MS6.1Panzhihua－Huili earthquake sequence of Aug 30，2008.Chinese J.Geophys.（in Chinese），2010，53（12）：2852－2860.

［28］ 易桂喜，龍鋒，張致偉.汶川MS8.0地震余震震源機制時空分布特征.地球物理學報，2012，55（4）：1213－1227.Yi G X，Long F，Zhang Z W.Spatial and temporal variation of focal mechanisms for aftershocks of the 2008 MS8.0 Wenchuan earthquake.Chines J.Geophys（in Chinese），2012，55（4）：1213－1227.

［29］ 梁尚鴻，李幼銘，束沛鎰等.利用區(qū)域地震臺網P、S振幅比資料測定小震震源參數.地球物理學報，1984，27（3）：247－257.Liang S H，Li Y M，Shu P Y，et al.On the determining of source parameters of small earthquakes by suing amplitude ratios of P and S from regional network observations.Chinese J.Geophys.（in Chinese），1984，27（3）：249－257.

［30］ 胡新亮，刁桂苓，馬瑾等.利用數字地震記錄的ˉP，ˉS振幅比資料測定小震震源機制解的可靠性分析.地震地質，2004，26（2）：347－354.Hu X L，Diao G L，Ma J，et al.Reliability analysis of focal mechanism solutions of micro－earthquakes determined from amplitude ration ofˉPandˉSrecorded by digital seismograph.Seismology and Geology（in Chinese），2004，26（2）：347－354.

［31］ 崔效鋒，謝富仁，張紅艷.川滇地區(qū)現代構造應力場分區(qū)及動力學意義.地震學報，2006，28（5）：451－461.Cui X F，Xie F R，Zhang H Y.Recent tectonic stress field zoning in Sichuan－Yunnan region and its dynamic interest.Acta Seismologica Sinica（in Chinese），2006，28（5）：451－461.

［32］ 程萬正，阮祥，張永久.川滇次級地塊震源機制解類型與一致性參數.地震學報，2006，28（6）：561－573.Cheng W Z，Ruan X，Zhang Y J.Types of focal mechanism solutions and parameter consistency of the sub－blocks in Sichuan and Yunnan Provinces.Acta Seismologica Sinica（in Chinese），2006，28（6）：561－573.

［33］ Julian B R，Miller A D，Foulger G R.Non－double－couple earthquakes 1.theory.Reviews of Geophysics，1998，36（4）：525－549.

［34］ Gephart J W.FMSI：A fortran program for inverting fault／slickenside and earthquake focal mechanism data to obtain the regional stress tensor.Computers ＆Geosciences，1990，16（7）：953－989.

［35］ Seeber L，Armbruster J G.Earthquakes as beacons of stress change.Nature，2000，407（6800）：69－72.

Study on earthquakes induced by water injection in Zigong－Longchang area，Sichuan

ZHANG Zhi－Wei1，CHENG Wan－Zheng1，LIANG Ming－Jian1，WANG Xiao－Shan2，LONG Feng1，XU Yan3，CHEN Wen－Kang3，WANG Shi－Yuan1
1 Earthquake Adiministration of Sichuan Province，Chengdu610041，China
2 Earthquake Adiministration of Hebei Province，Shijiazhuang050021，China
3 Earthquake Adiministration of Zigong Municipality，Zigong643000，China

Based on seismic data recorded by Zigong local digital seismic network and mobile stations from January 2007to August 2010，the correlation between pressure，water injection rate of Jia 33well and seismic activities of three cluster regions（A，B，and C）in Zigong－Longchang area is studied，and the seismogenic mechanism of induced earthquakes is studied by the space concentration distribution of shallow seismic activity and the difference of local stress state.The result shows that the injection pressure of Jia 33well has obvious regulatory effect on seismic activity of region B.At the stage of negative pressure injection，the effect of water injection rate on seismic activity is small，the seismic activity of region B is very weak，at thestage of injection pressure，when the pressure is continuously elevated，the seismicity is obviously increased，and the frequency and intensity present a good correlation with the water injection rate.According to the focal mechanism solutions and stress field characteristic of earthquakes in research area，the focal depth of three ML≥4earthquakes are all shallow，and the dislocation types all present reverse faulting，the strikes of the nodal plane and the fault near epicenter are consistent.Based on a lot of small earthquake focal mechanism solutions，the statistical results of nodal strike and P－axis direction in region B are relatively disperse than the results of region A.Compared with the existing results of regional stress field in the southeast of Sichuan Province，the stress field disturbance of region B is bigger than the result of region A.Under the action of a unified stress field，the nodal and P－axis direction of focal mechanism solutions and stress field in region A and B present difference，which have a certain relationship with injection pressure of Jia 33well.

Induced earthquake，Injection pressure，Focal mechanism solution，Stress field，Zigong－Longchang area

10.6038／j.issn.0001－5733.2012.05.021

P315

2011－12－06，2012－04－18收修定稿

地震科技星火計劃項目（XH1021Y）、國家自然科學基金項目（41090294）及中國地震局監(jiān)測預報司震情跟蹤定向工作任務（2012020107）聯(lián)合資助.

張致偉，男，1983年生，助理工程師，2009年碩士畢業(yè)于中國地震局蘭州地震研究所固體地球物理專業(yè)，主要從事數字地震學和地震活動性等方面的研究.E－mail：zzw1983107＠163.com

張致偉，程萬正，梁明劍等.四川自貢—隆昌地區(qū)注水誘發(fā)地震研究.地球物理學報，2012，55（5）：1635－1645，

10.6038／j.issn.0001－5733.2012.05.021.

Zhang Z W，Cheng W Z，Liang M J，et al.Study on earthquakes induced by water injection in Zigong－Longchang area，Sichuan.Chinese J.Geophys.（in Chinese），2012，55（5）：1635－1645，doi：10.6038／j.issn.0001－5733.2012.05.021.

（本文編輯 胡素芳）