紫坪鋪水庫(kù)水位變化對(duì)剪切波分裂參數(shù)的影響

劉莎， 吳朋

1 中國(guó)地震局地球物理研究所, 北京 100081 2 四川省地震局, 成都 610041

?

紫坪鋪水庫(kù)水位變化對(duì)剪切波分裂參數(shù)的影響

劉莎1， 吳朋2

1 中國(guó)地震局地球物理研究所, 北京 100081 2 四川省地震局, 成都 610041

本文通過(guò)對(duì)2006—2009年四川紫坪鋪水庫(kù)庫(kù)區(qū)8個(gè)地震臺(tái)站記錄的地震事件，采用剪切波分裂方法獲得了水庫(kù)庫(kù)區(qū)剪切波分裂參數(shù)，并結(jié)合地震活動(dòng)性與水庫(kù)水位之間的變化關(guān)系，分析了紫坪鋪水庫(kù)庫(kù)區(qū)地殼應(yīng)力的變化特征.剪切波分裂結(jié)果顯示該研究區(qū)域快波偏振方向有兩個(gè)，分別為北東向和北西向，充分體現(xiàn)了紫坪鋪水庫(kù)地區(qū)地殼應(yīng)力是由北西向的區(qū)域主壓應(yīng)力與南東走向的龍門(mén)山斷裂帶綜合作用的結(jié)果.慢波延遲時(shí)間平均值為5.8 ms·km-1，慢波延遲時(shí)間較大的地區(qū)位于庫(kù)壩和庫(kù)尾，分別是水庫(kù)蓄水排水引起地殼應(yīng)力變化最大的區(qū)域.對(duì)比慢波延遲時(shí)間的變化和水庫(kù)水位的變化顯示了慢波延遲時(shí)間與水庫(kù)水位之間的一致變化關(guān)系，揭示了水庫(kù)的蓄水排水對(duì)地殼應(yīng)力的影響.

快波偏振方向； 慢波延遲時(shí)間； 水位； 地殼應(yīng)力

1 引言

近些年來(lái)，隨著國(guó)內(nèi)外大型高壩水庫(kù)的不斷興建，水庫(kù)誘發(fā)地震事件也隨之增多.水庫(kù)蓄水或排水通常會(huì)引起地殼巖體中應(yīng)力集中和能量釋放而產(chǎn)生地震(秦嘉政等，2009).這些地震位置一般鄰近重要的水庫(kù)工程設(shè)施，而且震源較淺，震中烈度相對(duì)較高，往往具有很大的破壞性，造成大壩及附近建筑物的破壞和人員傷亡，甚至引起滑坡、坍塌等嚴(yán)重的次生災(zāi)害.因此，水庫(kù)地區(qū)地震的研究不僅是水利水電工程研究的內(nèi)容，也是地震學(xué)研究的重要內(nèi)容之一.

剪切波分裂作為一個(gè)有效的觀(guān)測(cè)手段，可用于對(duì)地殼應(yīng)力、流體壓力等方面的探索(張永久等，2010; Gao et al., 2011; Tang et al., 2005).剪切波分裂參數(shù)快波偏振方向和慢波延遲時(shí)間是用于分析地殼應(yīng)力場(chǎng)的狀態(tài)和各向異性參數(shù)的靜態(tài)和動(dòng)態(tài)探討的重要依據(jù)(Crampin, 1984; 高原等，1999).剪切波分裂對(duì)地殼中各向異性介質(zhì)較為敏感，與地殼應(yīng)力場(chǎng)和斷裂線(xiàn)的分布密切相關(guān).快波偏振方向與地殼微裂隙排列方向一致，與原地主壓應(yīng)力的方向一致(Crampin, 1978; Zhang et al., 2000; Musumeci et al.,2005).慢波延遲時(shí)間通常是評(píng)價(jià)各向異性程度的物理量，對(duì)地殼中微裂隙的幾何形態(tài)和密度分布非常敏感，可以反映地震發(fā)生前后應(yīng)力的積累與釋放過(guò)程，而且也可以很好地反映地殼中水壓變化所引起的微裂隙變化.

紫坪鋪水庫(kù)位于四川省都江堰市北西9 km處的岷江上游，是一座壩高156 m庫(kù)容11.12億立方米的大型水利工程.紫坪鋪水庫(kù)已經(jīng)建立了專(zhuān)門(mén)用來(lái)監(jiān)測(cè)水庫(kù)地區(qū)中小地震的數(shù)字遙測(cè)地震臺(tái)網(wǎng)，該地震臺(tái)網(wǎng)由7個(gè)臺(tái)站組成.在紫坪鋪水庫(kù)庫(kù)區(qū)還有油榨坪臺(tái)站，屬于四川地震臺(tái)網(wǎng).自2005年9月紫坪鋪水庫(kù)下閘蓄水之后，庫(kù)區(qū)弱震活動(dòng)頻繁，為紫坪鋪水庫(kù)庫(kù)區(qū)地震相關(guān)研究提供了豐富的數(shù)據(jù)資料(胡先明，2007；雷興林等，2008).

之前有學(xué)者(張永久等，2010)已經(jīng)對(duì)紫坪鋪水庫(kù)庫(kù)區(qū)剪切波分裂進(jìn)行探討，獲得了水庫(kù)庫(kù)區(qū)各向異性特征.該結(jié)果采用的是紫坪鋪水庫(kù)庫(kù)區(qū)2004年8月—2008年5月汶川地震之前的地震數(shù)據(jù)，本文選擇2006—2009年紫坪鋪水庫(kù)庫(kù)區(qū)的地震波形資料，獲得水庫(kù)庫(kù)區(qū)的各向異性特征，進(jìn)而分析水庫(kù)水位的變化與地殼應(yīng)力的變化關(guān)系，并對(duì)比張永久等(2010)的結(jié)果，探討兩次剪切波分裂獲得的地震各向異性特征有何異同.

2 地震和水位數(shù)據(jù)

紫坪鋪水庫(kù)位于龍門(mén)山斷裂帶中段，龍門(mén)山斷裂帶處于中國(guó)大陸南北地震帶的中部，東南側(cè)為揚(yáng)子準(zhǔn)地臺(tái)，西側(cè)為松潘—甘孜地塊.地震震源機(jī)制解結(jié)果表明龍門(mén)山斷裂帶所處地區(qū)的現(xiàn)代應(yīng)力場(chǎng)主壓應(yīng)力軸方位為北西西向(鄧起東等，1994).紫坪鋪水庫(kù)自2005年9月下閘蓄水之后弱震活動(dòng)較為活躍，為剪切波分裂研究提供了豐富的數(shù)據(jù)資料.本文選取了2006—2009年紫坪鋪水庫(kù)庫(kù)區(qū)8個(gè)地震臺(tái)站記錄的地震數(shù)據(jù)，包括紫坪鋪水庫(kù)數(shù)字地震臺(tái)網(wǎng)的7個(gè)臺(tái)站(八角BAJ，白巖BAY，桂花樹(shù)GHS，靈巖寺LYS，廟子坪MZP，桃子坪TZP，鉆洞子ZDZ)和四川數(shù)字地震臺(tái)網(wǎng)的油榨坪YZP臺(tái)站，8個(gè)臺(tái)站平均臺(tái)間距10 km左右，采樣率為100 sps (samples per second)(圖1).

在2006—2009年整個(gè)觀(guān)測(cè)時(shí)間段內(nèi)，紫坪鋪水庫(kù)水位起伏較大，水庫(kù)庫(kù)區(qū)弱震活動(dòng)較強(qiáng)(圖2).M-T圖顯示水庫(kù)庫(kù)區(qū)沒(méi)有大于4.0級(jí)的地震，地震震級(jí)主要集中在3.0級(jí)以下.而且水庫(kù)庫(kù)區(qū)地震活動(dòng)與水庫(kù)水位的升降變化有一定的相關(guān)性.2006年上半年水庫(kù)水位處于830 m左右，8月13日開(kāi)始，水庫(kù)水位開(kāi)始出現(xiàn)大幅度上升，到10月7日累計(jì)上升54 m，此后水位開(kāi)始下降，地震活動(dòng)從2006年8月開(kāi)始逐漸增強(qiáng).2007年5月20日水庫(kù)水位再次上升，12月14日水位達(dá)到高值后開(kāi)始下降，直到2008年6月降至820 m.2008年2月庫(kù)區(qū)地震活動(dòng)性較強(qiáng).2008年11月，水庫(kù)水位再次達(dá)到峰值隨后下降.2009年6月水庫(kù)水位重新開(kāi)始上升，地震活動(dòng)性在2009年下半年也表現(xiàn)出逐漸增強(qiáng)的趨勢(shì).

3 剪切波分裂方法

本文采用偏振分析法獲得剪切波分裂參數(shù)，偏振分析法雖然計(jì)算速度較慢，但是準(zhǔn)確率較高，直觀(guān)性強(qiáng)(Crampin et al., 2006).

剪切波分裂方法研究地殼各向異性要求選取的地震波記錄必須在剪切波窗口內(nèi).利用剪切波分裂探討地殼介質(zhì)的各向異性，還要求地震數(shù)據(jù)有較高的信噪比.根據(jù)淺層勘探給出的速度模型(表1)(趙珠等，1997)，該地區(qū)深度12 km，震中距20 km的地震入射角大約58°時(shí)，實(shí)際地震射線(xiàn)的入射角僅35°.因此本文選取58°的入射角作為剪切波窗口.對(duì)數(shù)據(jù)進(jìn)行進(jìn)一步的篩選，選取地震事件震源深度在5～12 km范圍內(nèi)，而且剪切波信噪比大于6.0.最后符合剪切波分裂計(jì)算條件的地震數(shù)據(jù)共512條(表2).

圖1 紫坪鋪水庫(kù)庫(kù)區(qū)地震(2006—2009年)和臺(tái)站分布Fig.1 Distribution of earthquakes (2006—2009) and stations in Zipingpu reservoir

圖2 紫坪鋪水庫(kù)庫(kù)區(qū)水庫(kù)水位變化(a)，地震M-T(b)及N-T(c)圖Fig.2 The water level (a), M-T (b) and N-T (c) of Zipingpu reservoir

圖3 BAJ(八角臺(tái))剪切波分裂分析示意圖

圖4 剪切波偏振方向等面積極射投影與等面積玫瑰圖

表1 剪切波速度模型Table 1 The velocity model of shear wave

在確定剪切波分裂參數(shù)時(shí)，文中選取剪切波到時(shí)前20個(gè)采樣點(diǎn)和到時(shí)后40個(gè)采樣點(diǎn)的波形數(shù)據(jù)窗口(時(shí)長(zhǎng)0.6 s)，將這60個(gè)采樣點(diǎn)數(shù)據(jù)分為三段分別繪制質(zhì)點(diǎn)偏振圖進(jìn)行分析(圖3).畫(huà)出所有滿(mǎn)足數(shù)據(jù)要求的質(zhì)點(diǎn)偏振圖，然后從這些圖中挑選出快剪切波質(zhì)點(diǎn)振動(dòng)為線(xiàn)性或近于線(xiàn)性的質(zhì)點(diǎn)偏振圖.由于各向異性介質(zhì)的存在，剪切波偏振方向第一個(gè)突然轉(zhuǎn)變的位置被認(rèn)為是快剪切波的到達(dá)，慢剪切波的到達(dá)同樣也是極化偏振方向的突然轉(zhuǎn)變，快、慢剪切波到時(shí)差即為慢波延遲時(shí)間(圖3).

4 快波偏振方向結(jié)果

快波偏振方向通常反映了原地主壓應(yīng)力的方向，與地殼中應(yīng)力誘發(fā)的定向排列的微裂隙面的走向一致(賴(lài)院根等，2006；高原等，1995；石玉濤等，2009).從表2可以看出，8個(gè)臺(tái)站均得到了30條以上的可靠結(jié)果，特別是BAJ臺(tái)有172個(gè)結(jié)果.剪切波分裂結(jié)果統(tǒng)計(jì)8個(gè)臺(tái)站的快波偏振方向均有兩個(gè)，本文分別計(jì)算剪切波分裂參數(shù)平均值以及對(duì)應(yīng)的均方根誤差.快波偏振方向的誤差結(jié)果都小于30°，表明各個(gè)臺(tái)站的快波偏振方向優(yōu)勢(shì)取向都較為明顯.結(jié)果顯示研究區(qū)域內(nèi)快波偏振方向主要有兩個(gè)，分別為北西向和北東向.

表2 四川紫坪鋪水庫(kù)庫(kù)區(qū)地震臺(tái)站剪切波分裂參數(shù)Table 2 The parameters of shear wave splitting in Zipingpu reservoir

快波偏振方向等面積投影圖(圖4)顯示，研究區(qū)域內(nèi)8個(gè)臺(tái)站快波偏振方向有兩個(gè)，但是其優(yōu)勢(shì)取向都較為明顯.BAJ臺(tái)，BAY臺(tái)，GHS臺(tái)，LYS臺(tái)，MZP臺(tái)，TZP臺(tái)和ZDZ臺(tái)7個(gè)臺(tái)站快波偏振方向顯示為北西向和北東向，而YZP臺(tái)的結(jié)果為北東向和近似的南北向.北東向的快波偏振方向與龍門(mén)山斷裂帶的走向一致，北西向的快波偏振方向與研究區(qū)域主壓應(yīng)力的方向一致.

5 慢波延遲時(shí)間與水庫(kù)水位

慢波延遲時(shí)間對(duì)地殼中微裂隙的幾何形態(tài)較為敏感，因此它在一定程度上反映了地殼中應(yīng)力的分布特點(diǎn)(Gao and Crampin, 2008; Liu et al., 2014).研究區(qū)域內(nèi)512條地震事件的慢波延遲時(shí)間平均值為5.8 ms·km-1.慢波延遲時(shí)間的最大值位于八角臺(tái)(BAJ)，兩個(gè)方向上的慢波延遲時(shí)間分別是(7.7±3.8) ms·km-1，(8.1±3.1) ms·km-1.其次是靈巖寺臺(tái)(LYS)，兩個(gè)方向上的慢波延遲時(shí)間分別是(6.5±2.9) ms·km-1，(5.4±2.0) ms·km-1.最小值位于油榨坪臺(tái)(YZP)，兩個(gè)方向上的慢波延遲時(shí)間分別是(3.8±1.6) ms·km-1，(4.2±1.7) ms·km-1.其余5個(gè)臺(tái)站慢波延遲時(shí)間多集中在5 ms·km-1左右(表2).慢波延遲時(shí)間較大的靈巖寺臺(tái)和八角臺(tái)分別位于紫坪鋪水庫(kù)的庫(kù)壩區(qū)和庫(kù)尾區(qū)(圖1)，正好是處于水庫(kù)蓄水引起地殼應(yīng)力變化較大的兩個(gè)區(qū)域.慢波延遲時(shí)間最小的油榨坪臺(tái)距離水庫(kù)最遠(yuǎn)，地殼應(yīng)力受水庫(kù)水位變化影響最小.可見(jiàn)，水庫(kù)的庫(kù)壩和庫(kù)尾地區(qū)地殼應(yīng)力受水庫(kù)水位的影響較大.

由于慢波延遲時(shí)間的離散性，本文采用 7點(diǎn)滑動(dòng)平均來(lái)表示慢波延遲時(shí)間在時(shí)間上的變化.慢波延遲時(shí)間的7點(diǎn)滑動(dòng)平均曲線(xiàn)顯示慢波延遲時(shí)間變化較為明顯，特別是BAJ，BAY, GHS和LYS臺(tái)站(圖5).BAJ臺(tái)站慢波延遲時(shí)間從2006年初開(kāi)始升高直到10月底升至12 ms·km-1，隨后慢波延遲時(shí)間開(kāi)始下降直到2007年4月.之后，慢波延遲時(shí)間再次升高直到2007年10月.2007年年底，慢波延遲時(shí)間一直維持在11 ms·km-1，隨后降低至7 ms·km-1.2008年2月，慢波延遲時(shí)間再次升高直到2008年10月，12月慢波延遲時(shí)間又一次降低，2009年7月降至4 ms·km-1，隨后開(kāi)始持續(xù)升高.BAY臺(tái)站在2006年和2007年，慢波延遲時(shí)間變化不明顯，2008年年初開(kāi)始稍有降低.2008年7月慢波延遲時(shí)間上升直到2009年初，隨后開(kāi)始降低，7月降至最低3.5 ms·km-1.GHS臺(tái)站慢波延遲時(shí)間也可以觀(guān)測(cè)到變化的趨勢(shì)，在2006年底和2008年11月分別可以觀(guān)測(cè)到慢波延遲時(shí)間的兩個(gè)峰值，達(dá)到7 ms·km-1.在LYS臺(tái)站，慢波延遲時(shí)間的變化較小，在2007年底慢波延遲時(shí)間出現(xiàn)低谷，低至4 ms·km-1.

圖5 慢波延遲時(shí)間與水庫(kù)水位對(duì)比圖

慢波延遲時(shí)間對(duì)地殼應(yīng)力場(chǎng)所引起的微裂隙變化非常敏感.紫坪鋪水庫(kù)庫(kù)區(qū)慢波延遲時(shí)間在觀(guān)測(cè)時(shí)間段內(nèi)表現(xiàn)出了很好的一致性變化.為了分析水庫(kù)水位對(duì)地殼應(yīng)力變化的影響，文中將慢波延遲時(shí)間變化與水庫(kù)水位變化進(jìn)行對(duì)比分析(圖5).結(jié)果顯示慢波延遲時(shí)間的變化和水庫(kù)水位之間存在很好的一致性.每當(dāng)水庫(kù)水位出現(xiàn)大幅度上升，慢波延遲時(shí)間也同步出現(xiàn)相應(yīng)的增大；水位下降之后，慢波延遲時(shí)間隨之減小.這一結(jié)果與之前的結(jié)果(張永久等，2010)較為吻合.水庫(kù)水位與慢波延遲時(shí)間的一致變化關(guān)系在Gao和Crampin(2003)關(guān)于慢剪切波時(shí)間延遲變化與壓力變化的試驗(yàn)中也得到了證實(shí).水庫(kù)的蓄水與排水通常會(huì)引起地殼流體壓力的變化，剪切波分裂慢波延遲時(shí)間可能會(huì)由于地殼應(yīng)力的微小變化而產(chǎn)生顯著的變化，該結(jié)論也與Crampin和Zatsepin(1997)所提出的APE(Anisotropic Poro-elasticity, 各向異性孔隙彈性)理論相吻合.

6 地震各向異性的討論

紫坪鋪水庫(kù)位于龍門(mén)山斷裂帶上，龍門(mén)山斷裂帶作為2008年汶川地震的發(fā)震斷裂，在構(gòu)造上極其復(fù)雜.快波偏振方向通常與區(qū)域主壓應(yīng)力方向或者活動(dòng)斷裂的走向表現(xiàn)一致.研究區(qū)域中，8個(gè)臺(tái)站均顯示有兩個(gè)快波偏振方向，BAJ臺(tái)，BAY臺(tái)，GHS臺(tái)，LYS臺(tái)，MZP臺(tái)，TZP臺(tái)和 ZDZ臺(tái)7個(gè)臺(tái)站剪切波分裂快波偏振方向顯示北東和北西兩個(gè)方向，而YZP臺(tái)快波偏振方向?yàn)楸睎|向和近似南北向.YZP臺(tái)站北東向的快波偏振方向與斷裂帶的走向一致，南北向的快波偏振方向與其他臺(tái)站結(jié)果存在較大差異，與斷裂帶走向和最大主壓應(yīng)力方向也不一致.南北向的快波偏振方向結(jié)果與張永久等(2010)的結(jié)果相一致，證實(shí)了本文結(jié)果的可靠性.YZP臺(tái)站位于紫坪鋪水庫(kù)東南方向，慢波延遲時(shí)間的結(jié)果顯示該臺(tái)站受到水庫(kù)影響較小，該臺(tái)站南北向快波偏振方向反映了YZP臺(tái)站附近復(fù)雜的應(yīng)力場(chǎng)，該區(qū)域應(yīng)力方向并不僅僅是受到龍門(mén)山活動(dòng)斷裂帶單一的影響，也暗示了附近隱伏斷層的情況.以往的研究結(jié)果(高原等，1995；賴(lài)院根等，2006；石玉濤等2009)也表明快波偏振方向也揭示了新的構(gòu)造.

對(duì)比張永久等(2010)的結(jié)果(圖6)，BAJ臺(tái)，LYS臺(tái)，ZDZ臺(tái)，YZP臺(tái)和BAY臺(tái)快波偏振方向均有兩個(gè)方向.本文的結(jié)果中MZP臺(tái)，TZP臺(tái)和GHS臺(tái)仍然有兩個(gè)快波偏振方向，而之前的結(jié)果(張永久等，2010)顯示只有一個(gè)快波偏振方向.本文的地震數(shù)據(jù)采用的是2006—2009年，跨越汶川地震.通常大地震的發(fā)生會(huì)使得震區(qū)地殼應(yīng)力進(jìn)行調(diào)整，進(jìn)而快波偏振方向表現(xiàn)的較為零散，所以本文觀(guān)測(cè)到這三個(gè)臺(tái)站快波偏振方向有兩個(gè).較之前的結(jié)果，MZP臺(tái)和TZP臺(tái)快波偏振方向多了北東向，與龍門(mén)山活動(dòng)斷裂的走向一致.而GHS臺(tái)多了北西向的快波偏振方向.

紫坪鋪水庫(kù)8個(gè)地震臺(tái)站慢波延遲時(shí)間平均值為5.8 ms·km-1，與張永久等(2010)的結(jié)果5.38 ms·km-1近似.位于庫(kù)尾的白巖臺(tái)(BAY)和庫(kù)壩的靈巖寺臺(tái)(LYS)慢波延遲時(shí)間明顯大于其他幾個(gè)臺(tái)站.水庫(kù)的庫(kù)壩和庫(kù)尾相比于水庫(kù)其他地區(qū)地殼水壓較大，使得微裂隙孔隙流體壓力增強(qiáng)，慢波延遲時(shí)間將表現(xiàn)出較大的值，這一結(jié)果表明紫坪鋪水庫(kù)庫(kù)壩和庫(kù)尾附近是地殼應(yīng)力相對(duì)較大的區(qū)域.對(duì)比慢波延遲時(shí)間與水庫(kù)水位之間的關(guān)系表明兩者之間的變化有一定的相關(guān)性，水庫(kù)水位的增加往往伴隨著慢波延遲時(shí)間的增大，這是由于水庫(kù)蓄水導(dǎo)致了地殼中孔隙壓力的擴(kuò)散(Talwani,1997；陶瑋等，2014)，使得地殼裂隙孔隙壓強(qiáng)增大，進(jìn)而觀(guān)測(cè)到慢波延遲時(shí)間的增大.

圖6 紫坪鋪水庫(kù)庫(kù)區(qū)快剪切波偏振方向

本文使用的地震數(shù)據(jù)資料跨越汶川地震，對(duì)比汶川地震之前的剪切波分裂結(jié)果(張永久等，2010)發(fā)現(xiàn)，紫坪鋪庫(kù)區(qū)8個(gè)地震臺(tái)站快波偏振方向結(jié)果基本一致，慢波延遲時(shí)間反映的各向異性程度也相當(dāng).水庫(kù)周?chē)穆ㄑ舆t時(shí)間主要受到水庫(kù)水位的影響，與地殼水壓的變化相關(guān).本文得到的剪切波分裂結(jié)果初步反映了水庫(kù)庫(kù)區(qū)的區(qū)域應(yīng)力情況以及水庫(kù)的蓄水排水對(duì)地殼應(yīng)力的影響.

7 結(jié)論

紫坪鋪水庫(kù)自2005年9月下閘蓄水之后，水庫(kù)庫(kù)區(qū)弱震活動(dòng)頻繁. 而且紫坪鋪水庫(kù)已經(jīng)建立了專(zhuān)門(mén)用于監(jiān)測(cè)水庫(kù)地區(qū)中小地震的數(shù)字遙測(cè)地震臺(tái)網(wǎng)，積累了大量的地震數(shù)據(jù)資料. 本文利用2006—2009年紫坪鋪水庫(kù)庫(kù)區(qū)8個(gè)地震臺(tái)站記錄的地震事件，采用剪切波分裂方法分析了水庫(kù)庫(kù)區(qū)地殼應(yīng)力的分布特征，同時(shí)也探討了水庫(kù)水位的變化對(duì)剪切波分裂參數(shù)的影響.

紫坪鋪水庫(kù)庫(kù)區(qū)剪切波分裂快波偏振方向主要有兩個(gè)，北東向和北西向，反映了研究區(qū)域內(nèi)地殼各向異性的分布特征. 本文的結(jié)果與之前的剪切波分裂結(jié)果表現(xiàn)一致，紫坪鋪水庫(kù)地區(qū)快波偏振方向和區(qū)域主壓應(yīng)力的方向、龍門(mén)山斷裂的走向一致. 紫坪鋪水庫(kù)8個(gè)臺(tái)站慢波延遲時(shí)間平均值為5.8 ms·km-1. 水庫(kù)的庫(kù)壩和庫(kù)尾相比于水庫(kù)其他地區(qū)慢波延遲時(shí)間要大，表明庫(kù)壩和庫(kù)尾是地殼應(yīng)力相對(duì)較大的地區(qū).水庫(kù)庫(kù)區(qū)慢波延遲時(shí)間與水庫(kù)水位之間表現(xiàn)出了很好的一致性變化. 水庫(kù)的蓄水和排水通常會(huì)引起地殼水壓的變化，進(jìn)而使得地殼裂隙孔隙壓強(qiáng)發(fā)生相應(yīng)的變化.

致謝 感謝四川省地震局為本文的研究提供數(shù)據(jù)的支持和幫助；感謝兩位匿名評(píng)審專(zhuān)家提出的寶貴意見(jiàn).

Crampin S. 1978. Seismic-wave propagation through a cracked solid: polarization as a possible dilatancy diagnostic.Geophys.J.Int., 53(3): 467-496. Crampin S. 1984. An introduction to wave propagation in anisotropic media.Geophys.J.Int., 76(1): 17-28.

Crampin S, Gao Y. 2006. A review of techniques for measuring shear-wave splitting above small earthquakes.Phys.Earth.Planet.Inter., 159: 1-14.

Crampin S, Zatsepin S V. 1997. Modelling the compliance of crustal rock—II. Response to temporal changes before earthquakes.Geophys.J.Int., 129(3): 495-506.

Deng Q D, Chen S F, Zhao X L. 1994. Tectonics, seismicity and dynamics of Longmenshan Mountains and its adjacent regions.SeismologicalandGeology(in Chinese), 16(4): 389-403.

Gao Y, Zheng S H, Sun Y. 1995. Crustal crack anisotropy in Tangshan.ActaSeismologicaSinica(in Chinese), 17(3): 283-293.

Gao Y, Zheng S H, Zhou H L. 1999. Polarization patterns of fast shear wave in Tangshan region and their variations.ChineseJ.Geophys. (in Chinese), 42(2): 228-232.

Gao Y, Crampin S. 2003. Temporal variations of shear-wave splitting in field and laboratory studies in China.J.Appl.Geophys., 54(3-4): 279-287.

Gao Y, Crampin S. 2008. Shear-wave splitting and earthquake forecasting.TerraNova, 20(6): 440-448.

Gao Y, Wu J, Fukao Y, et al. 2011. Shear wave splitting in the crust in North China: stress, faults and tectonic implications.Geophys.J.Int., 187(2): 642-654.

Hu X M. 2007. Seismicity background around Zipingpu Dam before impounding.EarthquakeResearchinSichuan(in Chinese), (2): 16-21.

Lai Y G, Liu Q Y, Chen J H, et al. 2006. Shear wave splitting and the features of the crustal stress field in the Capital circle.ChineseJ.Geophys. (in Chinese), 49(1): 189-196.

Lei X L, Ma S L, Wen X Z, et al. 2008. Integrated analysis of stress and regional seismicity by surface loading—A case study of Zipingpu reservoir.SeismologyandGeology(in Chinese), 30(4): 1046-1064.

Liu S, Crampin S, Luckett R, et al. 2014. Changes in shear wave splitting before the 2010 Eyjafjallaj?kull eruption in Iceland.Geophys.J.Int., 199(1): 102-112.

Musumeci C, Patanè D, Scarfi L, et al. 2005. Stress directions and shear-wave anisotropy: observations from local earthquakes in southeastern Sicily, Italy.Bull.Seismol.Soc.Am., 95(4): 1395-1374.

Qin J Z, Liu L F, Qian X D. 2009. Research on characteristics and prediction of reservoir induced seismicity.JournalofSeismologicalResearch(in Chinese), 32(2): 105-113.

Shi Y T, Gao Y, Zhao C P, et al. 2009. A study of seismic anisotropy of Wenchuan earthquake sequence.ChineseJ.Geophys. (in Chinese), 52(2): 398-407.

Talwani P. 1997. On the nature of reservoir-induced seismicity.PureAppl.Geophys., 150(3-4): 473-492, doi: 10.1007/s000240050089.

Tang C H, Rial J A, Lees J M. 2005. Shear-wave splitting: a diagnostic tool to monitor fluid pressure in geothermal fields.Geophys.Res.Lett., 32(21): L21317, doi: 10.1029/2005GL023551.

Tao W, Masterlark T, Shen Z K, et al. 2014. Triggering effect of the Zipingpu Reservoir on the 2008MW7.9 Wenchuan, China, Earthquake due to poroelastic coupling.ChineseJ.Geophys. (in Chinese), 57(10): 3318-3331, doi: 10.6038/cjg20141019.

Zhang Y J, Gao Y, Shi Y T, et al. 2010. The shear-wave splitting study of Sichuan Zipingpu reservoir region.ChineseJ.Geophys. (in Chinese), 53(9): 2091-2101, doi: 10.3969/j.issn.0001-5733.2010.09.009.

Zhang Z J, Li Y K, Lu D Y, et al. 2000. Velocity and anisotropy structure of the crust in the Dabieshan orogenic belt from wide-angle seismic data.Phys.EarthPlanet.Inter., 122(1-2): 115-131.

Zhao Z, Fan J, Zheng S H, et al. 1997. Precision determination of the crustal structure and hypocentral locations in the Longmenshan thrust nappe belt.ActaSeismologicaSinica(in Chinese), 19(6): 615-622.

附中文參考文獻(xiàn)

鄧起東, 陳社發(fā), 趙小麟. 1994. 龍門(mén)山及其鄰區(qū)的構(gòu)造和地震活

動(dòng)及動(dòng)力學(xué). 地震地質(zhì), 16(4): 389-403.

高原, 鄭斯華, 孫勇. 1995. 唐山地區(qū)地殼裂隙各向異性. 地震學(xué)報(bào), 17(3): 283-293.

高原, 鄭斯華, 周惠蘭. 1999. 唐山地區(qū)快剪切波偏振圖像及其變化. 地球物理學(xué)報(bào), 42(2): 228-232.

胡先明. 2007. 紫坪鋪水庫(kù)蓄水前天然地震活動(dòng). 四川地震, (2): 16-21.

賴(lài)院根, 劉啟元, 陳九輝等. 2006. 首都圈地區(qū)橫波分裂與地殼應(yīng)力場(chǎng)特征. 地球物理學(xué)報(bào), 49(1): 189-196.

雷興林, 馬勝利, 聞學(xué)澤等. 2008. 地表水體對(duì)斷層應(yīng)力與地震時(shí)空分布影響的綜合分析——以紫坪鋪水庫(kù)為例. 地震地質(zhì), 30(4): 1046-1064.

秦嘉政, 劉麗芳, 錢(qián)曉東. 2009. 水庫(kù)誘發(fā)地震活動(dòng)特征及其預(yù)測(cè)方法研究. 地震研究, 32(2): 105-113.

石玉濤, 高原, 趙翠萍等. 2009. 汶川地震余震序列的地震各向異性. 地球物理學(xué)報(bào), 52(2): 398-407.

陶瑋, Masterlark T, 沈正康等. 2014. 紫坪鋪水庫(kù)造成孔隙彈性耦合變化及其對(duì)2008年汶川地震觸發(fā)作用. 地球物理學(xué)報(bào), 57(10): 3318-3331, doi: 10.6038/cjg20141019.

張永久, 高原, 石玉濤等. 2010. 四川紫坪鋪水庫(kù)庫(kù)區(qū)地震剪切波分裂研究. 地球物理學(xué)報(bào), 53(9): 2091-2101, doi: 10.3969/j.issn.0001-5733.2010.09.009.

趙珠, 范軍, 鄭斯華等. 1997. 龍門(mén)山斷裂帶地殼速度結(jié)構(gòu)和震源位置的精確修定. 地震學(xué)報(bào), 19(6): 615-622.

(本文編輯 何燕)

The effect of water level changes in Zipingpu reservoir on the parameters of shear wave splitting

LIU Sha1, WU Peng2

1InstituteofGeophysics,ChinaEarthquakeAdministration,Beijing100081,China2EarthquakeAdministrationofSichuanProvince,Chengdu610041,China

The water storage and drainage of the reservoir usually cause the stress change in the crust, which can induce earthquakes. The Zipingpu reservoir is located in the middle section of Longmenshan fault with complex structure. The digital seismic network has been established in the reservoir to monitor moderate and small earthquakes, and accumulated a wealth of data, facilitating the study of earthquakes in the reservoir.

Shear wave splitting can be used to detect the stress change in the crust, and obtain the anisotropic static and dynamic gauge. The polarization direction of fast shear waves is consistent with the direction of cracks in the crust, indicating the regional stress direction in the crust. The time delay of the slow shear waves is related to the distribution of the micro-cracks in the crust which implies the change of water pressure in the crust. And the delay time can be used to study the effect of reservoir water level on the crustal anisotropy. This paper obtained the shear wave splitting parameters by the polarization analysis with high accuracy and strong visuality, using seismic data recorded by 8 seismic stations in the Zipingpu reservoir area from 2006 to 2009. According to the velocity model of the shallow exploration, the shear wave window chosen in this paper is 58 degree. The focal depths of the earthquakes are in the range of 5～12 km, and the signal to noise ratio of shear waves is greater than 6.0. 512 seismic events are justified by the data selection criteria. On the basis of the particle moving figures, we choose the position where the particle trajectory deflects vertically to ascertain the arrival time of fast and slow shear waves.

Every seismic station recorded more than 30 events in the study area, and the polarization direction at each station is obvious. Combining with the relationship between the seismic activity and water level of the Zipingpu reservoir, this paper obtains the variation characteristics of the stress in the crust. The results of shear wave splitting show two polarization directions of fast shear waves, which are northeast and northwest. The northeast polarization direction of fast shear waves is consistent with the strike of the Longmenshan fault. And the northwest direction agrees with the direction of the regional principal compressive stress. The average delay time of the slow shear waves in the reservoir is 5.8 ms·km-1. The large delay time of the slow shear waves is from the station BAJ and LYS at the dam and reservoir tail, respectively, where the stress change is large because of the water storage and drainage in the reservoir. The small delay time is from the station YZP far from the reservoir, where the stress in the crust is affected least by the water level of the reservoir.

Comparison with the previous results of shear wave splitting in the Zipingpu reservoir shows agreement with our results. The polarization direction of fast shear waves is consistent with the direction of the regional principal compressive stress and the strike of active faults. The relationship between the delay time of the slow shear waves and the water level of the reservoir suggests the increase of water level in the reservoir is often accompanied with the increase of the delay time. The reservoir storage causes the diffusion of the pore pressure in the crust, which makes the aspect ratio of the microcracks in the crust increase.

Polarization direction of fast shear waves; Delay time of slow shear waves; Water level; Stress in the crust

10.6038/cjg20151118.Liu S, Wu P. 2015. The effect of water level changes in Zipingpu reservoir on the parameters of shear wave splitting.ChineseJ.Geophys. (in Chinese),58(11):4106-4114,doi:10.6038/cjg20151118.

中國(guó)地震局基本科研業(yè)務(wù)費(fèi)(DQJB13B07)，“2013年臺(tái)站三結(jié)合——四川地區(qū)S波分裂研究”(2013-48)共同資助.

劉莎，女，1985年生，助理研究員，主要從事地球介質(zhì)地震各向異性研究.E-mail：liusha@cea-igp.ac.cn

10.6038/cjg20151118

P315

2015-04-28，2015-09-28收修定稿

劉莎，吳朋. 2015. 紫坪鋪水庫(kù)水位變化對(duì)剪切波分裂參數(shù)的影響.地球物理學(xué)報(bào)，58(11):4106-4114,