張 華，黃樹(shù)生，陸麗娟，黃惠寧，張 帆

(廣西壯族自治區(qū)地震局，南寧 530022)

廣西南丹大廠礦區(qū)震源參數(shù)研究

張 華，黃樹(shù)生，陸麗娟，黃惠寧，張 帆

(廣西壯族自治區(qū)地震局，南寧 530022)

利用廣西南丹縣大廠礦區(qū)地震監(jiān)測(cè)臺(tái)網(wǎng)記錄的地震波資料，采用Brune(1970)模型，將速度記錄譜歸算為震源位移譜，使用遺傳算法計(jì)算拐角頻率及零頻極限，然后計(jì)算地震矩、應(yīng)力降等小震震源參數(shù)，并系統(tǒng)分析各種參數(shù)之間的關(guān)系。大廠礦區(qū)38次較大地震的震源參數(shù)計(jì)算結(jié)果為：地震矩范圍在2.18×1011～7.89×1012N·m之間；震源破裂半徑78～439 m; 地震應(yīng)力降在0.01～1.79 MPa之間。表明：地震矩與震級(jí)、震源半徑、拐角頻率之間呈線性關(guān)系，而應(yīng)力降不依賴(lài)于地震矩，分布較為離散。大廠礦區(qū)應(yīng)力降值明顯偏低，這可能與該區(qū)域背景應(yīng)力低有關(guān)。

震源參數(shù)；定標(biāo)關(guān)系；大廠礦區(qū)

0 引言

通常認(rèn)為在地震孕育過(guò)程中，地殼介質(zhì)處于高應(yīng)力狀態(tài)，其輻射的地震波頻譜特征、衰減特性等和低應(yīng)力狀態(tài)是不同的。因此，對(duì)某一地區(qū)的地震進(jìn)行波譜分析，就可以得到表征震源破裂特征的震源參數(shù)(包括震源尺度、破裂方向和應(yīng)力降等)，并由此推斷震源區(qū)的構(gòu)造應(yīng)力環(huán)境，能夠使人們?cè)谝欢ǔ潭壬狭私庹鹪磪^(qū)的應(yīng)力狀態(tài)[1]。大量頻繁發(fā)生的中小地震震源參數(shù)和介質(zhì)參數(shù)攜帶了震源深度處應(yīng)力場(chǎng)和介質(zhì)性質(zhì)變化的信息，通過(guò)監(jiān)測(cè)這些參數(shù)的時(shí)空演化過(guò)程，可以對(duì)孕震區(qū)應(yīng)力場(chǎng)或介質(zhì)性質(zhì)變化過(guò)程進(jìn)行研究[2]。目前已有很多作者對(duì)構(gòu)造地震、水庫(kù)地震及火山地震的震源參數(shù)進(jìn)行了研究，但對(duì)礦區(qū)地震震源參數(shù)的研究并不多見(jiàn)。

震源參數(shù)的計(jì)算一般在頻率域進(jìn)行，在對(duì)臺(tái)站記錄的地震波資料進(jìn)行場(chǎng)地響應(yīng)、儀器響應(yīng)、幾何擴(kuò)散、路徑衰減、震源輻射方向性因子校正之后[3-5]，通過(guò)遺傳算法[6-8]可以獲得低頻波譜振幅極限、拐角頻率等參數(shù)，進(jìn)而根據(jù)一定的理論公式可計(jì)算應(yīng)力降、震源半徑、地震矩等[9]。

廣西南丹縣大廠錫多金屬礦田是世界聞名的超大型礦床，是我國(guó)錫多金屬礦業(yè)的重要基地，礦產(chǎn)資源非常豐富。作為有上千年開(kāi)采歷史的有色金屬礦區(qū)，大廠礦區(qū)地下開(kāi)采坑道縱橫交錯(cuò)，已形成約83 km2范圍內(nèi)的地下作業(yè)采空區(qū)，從而導(dǎo)致地震活動(dòng)日趨頻繁。特別是2007年以來(lái)，礦區(qū)連續(xù)發(fā)生多次有明顯震感的地震活動(dòng)，對(duì)當(dāng)?shù)厝罕娚a(chǎn)生活秩序和礦山安全生產(chǎn)造成了較大的影響[10]。深入了解大廠礦區(qū)的震源破裂特征，對(duì)礦山安全生產(chǎn)和社會(huì)秩序穩(wěn)定具有重要的現(xiàn)實(shí)意義。本文利用廣西大廠地震臺(tái)網(wǎng)記錄的數(shù)字地震資料對(duì)該礦區(qū)礦震的震源參數(shù)進(jìn)行了計(jì)算研究。

1 計(jì)算方法

地面運(yùn)動(dòng)的位移譜可以表示為：

Aij(f)=Ai0(f)·G(Rij)·Sj(g)·e-k(f)Rij

(1)

式中：Aij(f)為第i個(gè)地震在第j個(gè)臺(tái)站的觀測(cè)譜；Ai0(f)為第i個(gè)地震的震源譜；Rij為第i個(gè)地震到第j個(gè)臺(tái)站的震源距；G(Rij)為幾何擴(kuò)散函數(shù)；Sj(f)為第j個(gè)臺(tái)站的場(chǎng)地響應(yīng)；k(f)為非彈性衰減系數(shù)。對(duì)于某次地震的某個(gè)臺(tái)站記錄，從振幅譜中扣除儀器響應(yīng)、噪聲影響、傳播路徑影響、場(chǎng)地響應(yīng)，就可以得到該臺(tái)站記錄的震源位移譜A0j(f)。這樣，對(duì)于某次地震的每個(gè)臺(tái)站記錄都求出震源位移譜A0j(f)，再通過(guò)式(2)求平均震源譜來(lái)作為該地震的觀測(cè)震源譜

(2)

根據(jù)Brune模型，理論震源位移譜可以表示為：

|Oth(f)|=Ω0·[1+(f/f0)2]-1

(3)

式中：Oth(f)表示理論震源譜，Ω0表示震源譜低頻極限值，f0表示拐角頻率。Ω0和f0為地震位移譜的2個(gè)獨(dú)立參數(shù)。通過(guò)觀測(cè)數(shù)據(jù)計(jì)算模型參數(shù)有多種方法，本研究采用的是遺傳算法，以零頻極限和拐角頻率為獨(dú)立變量，確定下式為目標(biāo)函數(shù)[8]。根據(jù)儀器性能及采樣率，限定在1～18 Hz之間搜索頻率，使目標(biāo)函數(shù)最小，進(jìn)而獲得Ω0和f0。

(4)

式中：p為頻率點(diǎn)數(shù)，ε表示目標(biāo)函數(shù)值，Oob為同一地震的波譜轉(zhuǎn)化為震源位移譜后的多臺(tái)平均值，Oth為根據(jù)ω平方模型計(jì)算的震源譜，k為頻率點(diǎn)序號(hào)。

并由下式獲得地震矩M0，

M0=b·Ω0

(5)

b=4πβ3ρ·[Rs(φ,θ)·Re]-1

(6)

式中：Rs(φ,θ)為輻射方向性因子，Re為自由表面反射系數(shù)，ρ表示地殼介質(zhì)密度，β表示S波速度，應(yīng)力降Δσ、震源半徑r0分別使用(7)、(8)式計(jì)算：

(7)

(8)

2 大廠礦區(qū)地質(zhì)構(gòu)造與地震活動(dòng)

大廠礦區(qū)位于江南古陸西南緣的NW 向(南)丹-(河)池坳褶帶中段。從晚古生代至今，礦區(qū)附近地區(qū)經(jīng)歷過(guò)不同的地質(zhì)發(fā)展時(shí)期和多期構(gòu)造運(yùn)動(dòng)，地質(zhì)構(gòu)造較為復(fù)雜。礦區(qū)內(nèi)出露地層主要為一套泥盆系-二疊系的碎屑巖-碳酸鹽巖-硅質(zhì)巖組合，泥盆系為主要賦礦層位[11]。自中生代以來(lái)，區(qū)內(nèi)構(gòu)造經(jīng)歷了印支期強(qiáng)烈擠壓及其后的伸展剪切作用，前者表現(xiàn)為緊閉的線形褶皺及逆沖斷層，后者表現(xiàn)為層內(nèi)伸展剪切、拉斷的石香腸構(gòu)造及張扭性斷層，成巖、成礦作用明顯受后期構(gòu)造所控制。區(qū)域性斷裂——丹池?cái)嗔炎员倍洗┻^(guò)近場(chǎng)區(qū)，將區(qū)內(nèi)斷裂系統(tǒng)分成兩大斷裂構(gòu)造域，東部斷裂構(gòu)造主要以北東向?yàn)橹?，西部斷裂?gòu)造主要以NW、NNW向?yàn)橹鳎鞑康拇髲S斷裂正是礦田內(nèi)的主干構(gòu)造，長(zhǎng)度約10 km，走向?yàn)镹W320°，傾向NE，傾角28°～72°，產(chǎn)狀上陡下緩，具有“犁式”逆沖斷裂特征[12]。

為更好地監(jiān)測(cè)礦區(qū)地震活動(dòng)，自2009年8月以來(lái)相繼在大廠礦區(qū)建設(shè)了8個(gè)遙測(cè)臺(tái)站，臺(tái)站間距3 km左右，可監(jiān)測(cè)震級(jí)下限低至ML-0.5級(jí)。該臺(tái)網(wǎng)的建立為礦震研究提供了豐富的數(shù)據(jù)，自2009年8月至2013年底，該臺(tái)網(wǎng)共記錄到ML0.0級(jí)以上地震415次，其中ML0.0～0.9級(jí)316次，ML1.0～1.9級(jí)95次，ML2.0～2.9級(jí)4次，最大地震為2010年4月20日6時(shí)34分大廠鎮(zhèn)山草溝(24.83°N，107.60°E)ML2.7級(jí)地震。

在計(jì)算震源參數(shù)前，首先采用遺傳算法對(duì)礦區(qū)地震進(jìn)行重新定位。為了驗(yàn)證遺傳算法在大廠礦區(qū)的定位效果，我們收集到5次典型的人工爆破資料，并將遺傳算法定位結(jié)果與實(shí)際爆破位置進(jìn)行了對(duì)比(表1)。從表1給出的大廠臺(tái)網(wǎng)對(duì)5次網(wǎng)內(nèi)爆破定位結(jié)果對(duì)比可知，用遺傳算法對(duì)該區(qū)域爆破的定位精度最小水平誤差為78 m，最大不超過(guò)280 m，平均定位水平誤差在146 m范圍內(nèi)，優(yōu)于臺(tái)網(wǎng)地震定位精度(500 m)的標(biāo)準(zhǔn)。定位深度方面，臺(tái)網(wǎng)報(bào)告(hypo81定位)深度大部分在5 km以上，最大深度達(dá)到11.2 km；通過(guò)遺傳算法得到的最終定位結(jié)果是震源深度為1.5～6.5 km，平均深度為2.8 km。圖1為基于遺傳算法重新定位后的結(jié)果[13]。

表1 爆破實(shí)際位置與遺傳算法定位結(jié)果對(duì)比表

圖1 重新定位后的震中分布圖

3 資料處理

3.1 資料選取

地震資料選取廣西大廠礦區(qū)地震監(jiān)測(cè)臺(tái)網(wǎng)的波形數(shù)據(jù)，大廠礦區(qū)位于107°30′～107°40′E， 24°49′～24°54′N(xiāo)范圍之間，地面面積約83 km2。大廠臺(tái)網(wǎng)包括8個(gè)遙測(cè)子臺(tái)，全部采用數(shù)字無(wú)線遙測(cè)方式進(jìn)行組網(wǎng)(圖2)。臺(tái)網(wǎng)采用FSS-3M地震計(jì)，采用大動(dòng)態(tài)24位數(shù)據(jù)采集器，動(dòng)態(tài)范圍不低于110 dB，采樣頻率為每秒100點(diǎn)，這些地震儀在1～40 Hz之間具有速度平坦的響應(yīng)。在重新定位的基礎(chǔ)上，選取大廠礦區(qū)2009年1月至2013年12月發(fā)生的1.0級(jí)以上地震，然后按照每個(gè)地震至少被3個(gè)臺(tái)站記錄到，每個(gè)臺(tái)站至少記錄到3個(gè)地震的原則，通過(guò)信噪比計(jì)算共挑選出8個(gè)臺(tái)站，38次地震用于震源波譜參數(shù)的計(jì)算。

圖2 大廠礦區(qū)地震臺(tái)站分布圖

3.2 地震波譜計(jì)算

取水平向S波記錄數(shù)據(jù)計(jì)算付氏譜，要求所取波段包含90%以上S波能量[14-15]，同時(shí)具有較高的信噪比。為了保證地震波譜數(shù)據(jù)的穩(wěn)定性，在計(jì)算付氏譜時(shí)采用了Chael(1987)提出的延時(shí)窗技術(shù)[16]。該技術(shù)先把S波段分割為256點(diǎn)的小數(shù)據(jù)段，計(jì)算付氏譜，之后以128點(diǎn)數(shù)據(jù)段前移，前次計(jì)算數(shù)據(jù)段的一半被重疊，進(jìn)行第2次計(jì)算，依次類(lèi)推，這樣就獲得了整個(gè)數(shù)據(jù)段分割、重疊的n段譜數(shù)據(jù)，用下式把n小段數(shù)據(jù)歸算為全部S波段波譜數(shù)據(jù)

(9)

式中：O(f)為觀測(cè)位移譜，n為數(shù)據(jù)段總數(shù)，T為整段S波窗長(zhǎng)，t為數(shù)據(jù)分段窗長(zhǎng)，f代表頻率。使用(9)式分別計(jì)算東西、南北分量震源譜，并根據(jù)(10)式合成最終的震源觀測(cè)譜[4]：

Oob(f)=[(Oew(f))2+(Ons(f))]1/2

(10)

式中：Oob(f)為觀測(cè)譜，Oew(f)、Ons(f)為東西向、南北向觀測(cè)譜。在計(jì)算S波位移譜的同時(shí)計(jì)算噪聲譜，噪聲太大的地震波記錄不參與計(jì)算。對(duì)同一地震使用多臺(tái)記錄的振幅譜均值作為有效振幅譜(圖3)。

a 拉甲臺(tái)記錄波形經(jīng)校正后的觀測(cè)譜與噪聲譜對(duì)比圖 b 譜數(shù)據(jù)擬合圖 注：a圖，藍(lán)色曲線表示南北向和東西向譜結(jié)果，紅色曲線表示平均結(jié)果，黑色曲線表示噪聲譜結(jié)果； b圖，藍(lán)色曲線表示各臺(tái)觀測(cè)譜數(shù)據(jù)結(jié)果，綠色曲線表示多臺(tái)平均結(jié)果，紅色曲線表示數(shù)據(jù)擬合結(jié)果。

4 計(jì)算結(jié)果

根據(jù)Brune模型理論及公式，計(jì)算得到了大廠礦區(qū)38次較大地震的震源波譜參數(shù)。通過(guò)對(duì)大廠礦區(qū)礦震的計(jì)算得到，其地震矩M0在2.18×1011～7.89×1012N·m之間；地震的拐角頻率在3～17 Hz 之間，主要集中在5～15 Hz之間；震源破裂半徑范圍78～439 m; 地震應(yīng)力降在0.01～1.79 MPa 之間。在地震矩與震級(jí)、震源半徑及拐角頻率的定標(biāo)關(guān)系方面，本文的研究結(jié)果同國(guó)內(nèi)外研究結(jié)果基本一致[17-20](表2)。

4.1 地震矩與近震震級(jí)、震源半徑的關(guān)系

標(biāo)量地震矩M0計(jì)算結(jié)果介于2.18×1011～7.89×1012N·m 之間。在單對(duì)數(shù)坐標(biāo)下，地震矩M0與震級(jí)ML之間呈現(xiàn)較好的相關(guān)性，隨著震級(jí)的增大，地震矩也相應(yīng)地增大。從圖4中看出，在所研究的震級(jí)范圍內(nèi)，地震矩與震級(jí)有很好的線性關(guān)系。地震矩與震級(jí)之間的線性關(guān)系式能反映不同地區(qū)的地下介質(zhì)狀態(tài)和地震活動(dòng)特點(diǎn)。用最小二乘法擬合出地震矩M0與震級(jí)ML之間關(guān)系式為：

LgM0=1.48ML+9.08

(11)

計(jì)算得出，38 次ML1.0級(jí)以上地震的破裂半徑在78～439 m之間，地震矩與震源半徑之間呈線性關(guān)系，并給出了用最小二乘法擬合得到的地震矩與震源半徑之間的關(guān)系式為：

LgM0=0.003R+11.63

(12)

4.2 地震矩與應(yīng)力降、拐角頻率的關(guān)系

38 次ML1.0級(jí)以上地震的應(yīng)力降計(jì)算結(jié)果介于 0.01～1.79 MPa之間，研究區(qū)應(yīng)力降值整體偏低，這種現(xiàn)象可能與地下介質(zhì)相對(duì)破碎、非均勻性特征突出、震源機(jī)制復(fù)雜及該區(qū)域背景應(yīng)力低有關(guān)。圖5給出了應(yīng)力降Δσ與地震矩M0之間的關(guān)系，可以看出應(yīng)力降不依賴(lài)于地震矩，分布較為離散。

從地震矩和拐角頻率的關(guān)系來(lái)看，二者明顯有依賴(lài)關(guān)系，地震矩越大，拐角頻率越小，這與以往的研究結(jié)果一致。根據(jù) Brune( 1970) 模型，用最小二乘法擬合出 38次ML1.0級(jí)以上地震拐角頻率和地震矩M0之間的關(guān)系式為

Lgf0=-0.338 LgM0+5.07

(13)

表2 相關(guān)參數(shù)擬合數(shù)據(jù)結(jié)果比較

圖4 大廠礦區(qū)地震矩與近震震級(jí)、震源半徑的關(guān)系

圖5 地震矩與應(yīng)力降、拐角頻率的關(guān)系

5 結(jié)論與討論

本文用廣西大廠礦區(qū)地震臺(tái)網(wǎng)觀測(cè)的數(shù)字波形資料，估算了該區(qū)域的震源參數(shù)，得到了一些有意義的結(jié)果：地震矩與震級(jí)、震源半徑之間呈線性關(guān)系；地震矩和拐角頻率有明顯的依賴(lài)關(guān)系，地震矩越大，拐角頻率越小，這與以往的研究結(jié)果一致；應(yīng)力降不依賴(lài)于地震矩，分布較為離散。

大廠礦區(qū)地震的拐角頻率在3～17 Hz之間，主要集中在5～15 Hz。而張永久等[9]給出紫坪鋪水庫(kù)的拐角頻率在2～13 Hz之間，劉學(xué)軍等[21]給出騰沖火山拐角頻率在4.37～11.87 Hz之間，姚安構(gòu)造地震的拐角頻率在3.10～11.40 Hz之間。相對(duì)一些地區(qū)，大廠礦區(qū)的拐角頻率相對(duì)較高，這可能與大廠礦區(qū)地震多為淺源地震，觀測(cè)臺(tái)網(wǎng)密集，記錄到的地震高頻成份較多有關(guān)。

震源參數(shù)研究中，大多數(shù)研究結(jié)果中提到了拐角頻率與地震大小的負(fù)相關(guān)關(guān)系[22-24]，即地震震級(jí)越大，地震矩越大，拐角頻率越小。大廠礦區(qū)地震矩和拐角頻率存在明顯負(fù)相關(guān)關(guān)系，這與以往的研究結(jié)果一致。但亦有研究表明，這種負(fù)相關(guān)關(guān)系在紫坪鋪水庫(kù)地區(qū)不成立[9]。姚立珣等[25]認(rèn)為：水庫(kù)誘發(fā)地震震源區(qū)的介質(zhì)條件變化與一般構(gòu)造地震不同，與水的作用有密切關(guān)系。巖體的擴(kuò)容、飽和以及浸水弱化對(duì)地震波傳播速度、破裂速度及破裂尺度有較大影響，而這些也是決定拐角頻率等震源參數(shù)的主要因素之一。受觀測(cè)資料樣本限制，本文未對(duì)大廠礦區(qū)地震拐角頻率的變化特征及原因進(jìn)行深入分析。

應(yīng)力降表征地震瞬間斷層錯(cuò)動(dòng)時(shí)位錯(cuò)面上的應(yīng)力變化，通過(guò)分析研究地震的應(yīng)力降可以認(rèn)識(shí)地震發(fā)生地區(qū)的構(gòu)造背景應(yīng)力。國(guó)內(nèi)外一些學(xué)者對(duì)不同類(lèi)型的地震應(yīng)力降進(jìn)行了研究，如華衛(wèi)等[26]從水庫(kù)誘發(fā)地震與構(gòu)造地震應(yīng)力降比較結(jié)果發(fā)現(xiàn)，在兩者可比較震級(jí)范圍內(nèi)(ML2.0～4.0)，雖然應(yīng)力降有一定的波動(dòng)范圍，但三個(gè)水庫(kù)地區(qū)(龍灘、三峽、新豐江)中相同震級(jí)的水庫(kù)地震，其應(yīng)力降值整體上比構(gòu)造地震明顯偏低，大約小10倍。而且，震級(jí)越小這種差異似乎越大。一些學(xué)者對(duì)火山地區(qū)的應(yīng)力降進(jìn)行了研究，發(fā)現(xiàn)幾乎所有活火山的火山地震的地震應(yīng)力降都比較低(如日本大島火山、意大利埃特納火山、美國(guó)圣安德烈斯火山和夏威夷火山及“九五”騰沖火山監(jiān)測(cè)研究等)，火山地震低應(yīng)力降現(xiàn)象的物理解釋是因?yàn)槠湔鹪礄C(jī)制的力學(xué)特性與構(gòu)造的不同，由于火山地震主要是火山巖漿活動(dòng)引起的，而熔融或半熔融狀態(tài)的高溫物質(zhì)是難以積累較高的破裂應(yīng)力的[27]。本文研究中大廠礦區(qū)礦震應(yīng)力降值整體明顯偏低，其原因?yàn)楹?？一些研究認(rèn)為，構(gòu)造地震平均來(lái)說(shuō)似乎比水壓破裂和礦震等誘發(fā)地震有較高的應(yīng)力降(大約為10倍)。這是因?yàn)榉亲匀坏沫h(huán)境局部地影響原來(lái)位置處的應(yīng)力水平，從而導(dǎo)致流體(或采礦)誘發(fā)的地震可能發(fā)生在比自然地震更低的構(gòu)造應(yīng)力地區(qū)[28-29]。結(jié)合地質(zhì)環(huán)境及地震活動(dòng)水平分析認(rèn)為，大廠礦區(qū)應(yīng)力降值整體偏低的現(xiàn)象可能與地下介質(zhì)相對(duì)破碎、非均勻性特征突出、震源機(jī)制復(fù)雜及該區(qū)域背景應(yīng)力低有關(guān)。

大廠礦區(qū)具有密集的臺(tái)網(wǎng)，對(duì)震級(jí)較小地震的記錄有很高的信噪比，對(duì)于更接近點(diǎn)源地震模型假設(shè)的小地震，能得出更為真實(shí)可信的地震震源參數(shù)，是礦震研究的天然試驗(yàn)場(chǎng)，可進(jìn)一步研究震源參數(shù)的時(shí)空演化過(guò)程，結(jié)合震源機(jī)制，反演其區(qū)域應(yīng)力場(chǎng)，探討應(yīng)力場(chǎng)變化與較大礦震的關(guān)系，以期對(duì)礦震發(fā)震成因研究進(jìn)行深入探索。

致謝：感謝浙江省地震局朱新運(yùn)博士提供計(jì)算程序。

[1]夏愛(ài)國(guó),趙翠萍.天山中東段地區(qū)震源參數(shù)的初步應(yīng)用研究[J].內(nèi)陸地震,2006,20(3):245-251.

[2]趙翠萍,陳章立,華衛(wèi),等.中國(guó)大陸主要地震活動(dòng)區(qū)中小地震震源參數(shù)研究[J].地球物理學(xué)報(bào),2011,54(6):1481-1489.

[3]Atkinson G M,R F Mereu.The shape of ground motion attenuation curves in Southeastern Canada[J]. Bull. Seism. Soc. Am. ,1992,82(5):2014-2031．

[4]Moya A,Jorge A.Inverson of source parameters and site effects from strong ground motion record using genetic algorithms[J].Bull.Seism.Soc.Am.,2000,90(4):977-992

[5]Bindi D,Spallarossa D，Augliera P，et al．．Source parameters estimated from the aftershock of the 1997 Umbria Marche(Italy)seismic sequence[J].Bull. Seis. Soc. Am. ,2001,91(3):448-455．

[6]劉杰,鄭斯華,黃玉龍.利用遺傳算法反演非彈性衰減系數(shù)、震源參數(shù)和場(chǎng)地響應(yīng)[J].地震學(xué)報(bào),2003,25(2)：211-218．

[7]郭曉,張?jiān)?莘海亮,等.祁連山中東段地區(qū)非彈性衰減系數(shù)、震源參數(shù)和場(chǎng)地響應(yīng)研究[J].西北地震學(xué)報(bào),2007,29(4):319-325.

[8]朱新運(yùn),于俊誼.基于MATLAB的小震震源參數(shù)計(jì)算軟件研制[J].西北地震學(xué)報(bào),2008,30(4):380-384．

[9]張永久,喬慧珍.紫坪鋪水庫(kù)庫(kù)區(qū)震源波譜參數(shù)研究[J].西北地震學(xué)報(bào),2011,33(2):117-122.

[10]姚宏,郭培蘭,龍政強(qiáng),等.廣西南丹大廠礦區(qū)地震活動(dòng)特征分析[J].地震,2011,31(3):112-120.

[11]蔡明海,梁婷,吳德成,等.廣西大廠礦田花崗巖地球化學(xué)特征及其構(gòu)造環(huán)境[J].地質(zhì)科技情報(bào),2004,23(2):57-62.

[12]陳毓川,黃民智,徐鈺,等.大廠錫礦地質(zhì)[M].北京:地質(zhì)出版社, 1993.

[13]張華,張忠利,姚宏,等.基于遺傳算法的大廠礦區(qū)地震定位研究[J].工程地球物理學(xué)報(bào),2014,11(2):260-265.

[14]蘭從欣,劉杰,鄭斯華,等.北京地區(qū)中小地震震源參數(shù)反演[J].地震學(xué)報(bào),2005,27(5):498-507.

[15]秦嘉政,葉建慶,錢(qián)曉東,等. 2000年姚安地震的震源參數(shù)[J].地球物理學(xué)報(bào),2003,46(5):633-641.

[16]Chael E. Spectral sealing of earthquakes in the Miramiehi region of New Brunswick[J].BSSA,1987,77,347-365

[17] Bakun W H.Seismic moment，local magnitudes and coda-duration magnitude for earthquakes in central California[J].BSSA,1984,74:439-458．

[18]于俊誼,朱新運(yùn).浙江珊溪水庫(kù)地震震源參數(shù)研究[J].中國(guó)地震,2008,24(4):379-387．

[19]劉麗芳,秦嘉政.會(huì)澤5.3級(jí)和文山5.3級(jí)地震序列震源參數(shù)研究[J].地震研究,2006,29(1):43-49.

[20]師海闊,朱新運(yùn),賀永忠,等.寧夏及鄰區(qū)地震震源參數(shù)研究[J].大地測(cè)量與地球動(dòng)力學(xué),2010,30(S1):38-43.

[21]劉學(xué)軍,葉建慶,蔡紹平,等. 騰沖火山區(qū)地震的震源參數(shù)[J].防災(zāi)減災(zāi)工程學(xué)報(bào),2004,24(1):22-27.

[22]劉學(xué)軍,葉建慶,蔡紹平,等.姚安地震序列的震源參數(shù)研究[J].地震研究,2002,25(3):240 -245．

[23]張永久,程萬(wàn)正.用S波研究雅江地震序列震源波譜[J].中國(guó)地震,2003,19(4):340-350．

[24]張永久,彭立國(guó),程萬(wàn)正.馬爾康地震序列震源參數(shù)研究[J].中國(guó)地震,2006,22(1):85-91．

[25]姚立珣,鐘羽云,張震峰,等.巖體裂隙充水后對(duì)拐角頻率的影響[J].西北地震學(xué)報(bào),2004,26(4):315-321

[26]華衛(wèi),陳章立,鄭斯華,等. 水庫(kù)誘發(fā)地震與構(gòu)造地震震源參數(shù)特征差異性研究——以龍灘水庫(kù)為例[J].地球物理學(xué)進(jìn)展,2012,27(3):924-935.

[27]李立平,李一正,葉建慶,等.瀾滄一耿馬地震震源動(dòng)力學(xué)參數(shù)研究[J].地震研究,1990,13(3):223-233.

[28] Abercrombie R E . Leary P.Source parameters of small earthquakes recorded at 2.5 km depth,Cajon Pass Southern California: implications for earthquake scaling[J]．Geophys．Res．Lett．,1993,20(14):1511-1514．

[29]董春麗,王向陽(yáng),鄭樹(shù)平,等. 山西省內(nèi)連續(xù)5次ML≥4.0地震序列震源區(qū)附近應(yīng)力降的初步研究[J]. 華北地震科學(xué),2013, 31(3):14-22.

Study on Source Parameters of Earthquakes in Dachang Mining Area

ZHANG Hua, HUANG Shu-sheng, LU Li-juan, HUANG Hui-ning, ZHANG fan

(Earthquake Administration of Guangxi Autonomous Region, Nanning 530022, China)

Using the Brune model (1970) and the digital seismic wave data of the earthquake occurred in Dachang mining area, which was recorded by Guangxi Dachang seismic network; we transferred the velocity spectrum into source displacement spectrum. The corner frequency and zero frequency were calculated by using genetic algorithm, further, the source parameters including seismic moment and stress drop were calculated and the relationships among them were analyzed systematically. This paper studies 38 bigger regional earthquakes occurred in this area, the results show that they seismic moments are 2.18×1011～7.89×1012Nom and rupture radius are 78～439 m, stress drop are 0.01～1.79 MPa. The relations between seismic moment, body-wave magnitude, rupture radius and corner frequency are near log-linear．The distribution of stress drop is dispersed and show little obvious correlation with seismic moment. Dachang mining stress drop values significantly lower, this may be related to the low regional background stress.

source parameter; scaling relation; Dachang mining area

10.3969/j.issn.1003-1375.2014.03.006

2014-03-31

地震科技星火計(jì)劃項(xiàng)目(項(xiàng)目編號(hào)：XH12036Y)；廣西科學(xué)研究與技術(shù)開(kāi)發(fā)計(jì)劃課題(項(xiàng)目編號(hào)：桂科攻1140004-3)

張華(1978-)，男，工程師，主要從事數(shù)字地震學(xué)及地震預(yù)測(cè)研究.E-mail:huazhang1222@163.com

P315.31

A

1003-1375(2014)03-0024-07