新疆數字地震臺網地方性震級量規函數的初步研究

張志斌 朱皓清 李艷永

新疆維吾爾自治區地震局，烏魯木齊市新市區科學二街338號 830011

0 引言

地震震級是表征地震強弱的量度，是地震三要素的基本參數之一。第1個震級標度是李克特根據古登堡與達清夫的建議在1935年提出的。Richter（1958）在研究美國南加州的地震時引入的地方性震級標度ML，曾被認為“如果沒有進一步的證據，不能應用于加利福尼亞以外的地區”。

20世紀50年代中期，李善邦先生結合我國常用的62型（短周期地震儀）和基式（SK，中長周期地震儀）儀器特性，將Richter提出的震級標度變換為適用于我國的量規函數R1和R2，并在全國范圍內使用。嚴尊國等（1995）在收集了東北、華北、華東、華南、西南和西北等地區地動位移資料后，利用震級殘差統計分析和地方性地震最大震相平均衰減形態來研究我國東部和西部的地方性震級量規函數，并確定了使用范圍，從而進一步完善了ML震級測定方法。隨著全國臺網數字化改造的完成，各臺網之間的數據實現了共享，在大震速報和地震編目等日常工作中也逐漸發現不同區域臺網之間對于同一地方性地震事件所測震級有明顯差別，同一地震事件，不同震中距臺站所測震級差可達2.0級。與此同時，河北、甘肅和山西等區域臺網利用近幾年積累的豐富的數字地震資料研究了本區域的量規函數（趙明淳等，2005；陳繼鋒等，2013；董春麗等，2007），提高了其測定震級的精度。

新疆地處印度板塊向歐亞板塊俯沖的前緣地帶，是中國大陸主要的內陸地震活動區，地震災害頻發，每年記錄到大小地震2萬余次。新疆測震臺網有關本區域量規函數的研究還停留在由高國英等（1987）結合模擬測震數據利用統計方法得到的結果。本文旨在利用新疆測震臺網近幾年積累的數字地震資料，遵循我國地殼地區性特征的ML震級量規函數的繼承性、區域性、實用性原則（嚴尊國等，1992），依據區域臺網地理范圍的行政分區邊界，運用震級殘差統計和最大地動位移衰減特性2種方法，研究新疆臺網地方性震級量規函數，以期為提高新疆測震臺網測定地方性震級精度提供參考。

1 研究方法

本文采用震級殘差統計和最大地動位移衰減特性2種方法，分別計算新疆臺網地方性震級量規函數。震級殘差統計實際上采用的是均勻震級系統的思想，是一種自我完善的方法，是對現階段使用的量規函數的一種修正和完善。利用最大地動位移衰減特性研究地方性震級量規函數是建立在Richter關于零級地震的定義基礎上，兩者都非另立系統。

1.1 震級殘差統計方法

我國的地方性震級ML標度系統始建于1959年，建立的ML測定公式為

式中，AN、AE分別為NS、EW方向記錄最大地動位移（無須嚴格追蹤到同一時刻，單位以 μm計）；R（Δ）為適合各地區實際情況的量規函數。

1971年，針對我國常用的短周期和中長周期地震儀建立了地方性震級量規函數R1（Δ）和R2（Δ），進一步完善了 ML震級測定方法。1978年將 R1（Δ）列入《地震臺站觀測規范》（試行）。多年來的實際觀測和研究表明，R1（Δ）和R2（Δ）在反映我國地殼的地區性特征方面還存在一些問題。震級殘差統計方法遵循了地方性震級測定的繼承性，計算出的地方性震級量規函數是對現有的量規函數的一種修正。

臺網震級Mi，即對第j個臺站，按照地方性震級公式求得Mij，然后對第i個地震所有記錄到的子臺按照下式求震級平均值和標準偏差

對于第i個地震，第j個子臺的震級偏差采用下式計算

式中，k為間隔序數；N為間隔范圍內震級偏差樣本數（郭履燦等，1984；陳運泰等，2004）。

1.2 最大地動位移衰減特性方法

1935年，Richter使用美國南加州的地震資料驗證了“用同樣的儀器在相同的震中距離上記錄兩次地震所得的最大振幅之比是一個與震中距無關的常數”，同時定義了零級地震為“在震中距為100km時，標準短周期扭轉型地震儀（T＝0.8s，放大率 V＝2800，阻尼常數 D＝0.8）所測最大地震波振幅為1μm”（林仙坎，2006）。利用最大地動位移衰減特征研究地方性震級量規函數就是建立在上述理論基礎上。

通過假定lg A－Δ的關系可以表示為

式中，i＝1～N；A為最大地動位移；Δ為震中距。通過“有進有出”的多元回歸分析方法找出lg A－Δ的關系式（中國科學院計算中心概率統計組，1979），從而擬合出最大地動位移衰減關系。利用零級地震的定義，進而求出地方性震級量規函數。

2 資料選取與處理結果

隨著“十五”數字化地震觀測網絡項目和中國地震背景場探測項目的建設，新疆測震臺網基本建成了一個覆蓋全疆、擁有84個數字測震臺站的地震觀測網絡。本文采用了新疆數字地震臺網2009年1月～2014年12月所記錄的地震事件。依托省級臺網之間數據共享，即新疆測震臺網與甘肅、西藏和青海相鄰臺網共享原始地震記錄，使新疆測震臺網部分記錄事件由網外變為網內。同時對所選的地震事件逐條都進行仔細地綜合分析，地震定位采用由廣東省地震局開發的MSDP交互處理軟件，測定震源基本參數時采用單純型法定位程序進行定位，并利用3400走時表求出發震時刻。

2.1 震級殘差統計方法

為了使被研究事件盡可能地覆蓋研究區域且綜合考慮近臺和遠臺的權重比，選取了2009年1月～2014年12月間2.0≤ML≤5.4的地震事件作為研究對象，同時要求每個地震事件至少有5個觀測臺站記錄，且震源深度小于 50km，最終共得到16269個地震樣本（圖1），震中距范圍取0～800km。在測定震級時，將速度記錄仿真成位移記錄后（消除周期對結果的影響），量取兩水平向最大振幅，然后根據式（1）求出各子臺震級值，繼而求出多臺震級平均值。

對所選用的新疆數字地震臺網84個子臺所記錄的179561個地震數據進行整理統計，按照式（2）、（4）計算得出單臺震級與平均震級的偏差，并做統計分析，結果見圖2。由圖2可見，偏差值基本呈正態分布，集中度較高，在被分析的179561個單臺震級中，偏差值（絕對值）在0.2以內的有106615次，占 59.3%；0.3～0.5有 60460次，占 33.7%；0.6～0.8有 11645次，占6.5%；大于0.9的841次，占0.5%。

采用式（3）～（5）得出各子臺震級相對于臺網震級偏差、平均偏差和標準偏差等，計算結果見表1。每個臺站的平均偏差范圍從－0.40～0.79，除了WCW（五彩灣）臺偏差較大外，剩余臺站都控制在±0.5的范圍內。標準偏差值反映了震級偏差值的離散程度，其中WCW（五彩灣）和STZ（石梯子）2個臺站的標準偏差較大，其余都在±0.5范圍內，說明震級偏差離散程度較小。其中WCW臺的臺基類型為砂巖，砂巖對地震具有放大效應，計算出來的震級也偏大；而STZ臺的臺基為礫巖，由其場地條件計算出的場地響應值較大，對地震震級也起到放大作用（劉建明等，2014）。

為了研究單臺震級隨震中距的變化，利用式（6）將前述179561個單臺震級偏差數據按照震中距由近到遠排序，在0～800km范圍內，以10km為間隔，分為80個間隔，求出每個間隔內其平均偏差值ΔMk。計算結果見圖（3）。

從圖3可以看出，當震中距在0～150km以內時，單臺震級與平均震級的偏差基本都超過了－0.18，在個別震中距段甚至接近－0.3，表明此震中距范圍量規函數偏小。當震中距超過150km時，單臺震級與平均震級的偏差基本都在±0.1的范圍內，震級偏差較小，取震級偏差的負值為校正值，對現有的量規函數進行校正，得出較正后的量規函數。

2.2 最大地動位移衰減特性方法

利用多元回歸分析研究地方性震級量規函數時，選取的“母本”資料對量規函數的實用性和科學性有很大影響。黃瑋瓊等（1982）在研究最大地動位移的衰減規律時，提出震級選取應相差不大，如若選取的震級偏大，體波震級已不能反映地震的真實大小，地方性震級量規函數已失去意義；若選取的震級偏小，新疆區域面積大，會造成遠臺無法記錄到地震，容易造成近臺權重過大，而失去地方性震級量規函數的真實性。

圖3 震級偏差值隨震中距的分布

圖4 最大地動位移衰減關系選用的地震震中分布

為了使不同震中距范圍內的觀測資料盡可能地等權，在選取資料時，要求每個地震記錄的臺數不少于10個，且震源深度不能超過50km，同時每個臺站記錄的地震條目也不能低于10個。依據上述條件選取了2009～2014年746個ML3.6～4.5地震作為研究對象（圖4），用以研究新疆地區最大地動位移衰減特性，進而得出地方性震級量規函數。

本研究所量取的最大地動位移都是由廣東省地震局開發的MSDP經W·A仿真處理后得到的。根據每個地震定位的震中位置，由各臺站的經緯度算出地震與臺站的距離-震中距Δ，再對各臺站記錄的EW向和NS向最大地動位移值平均數取對數，這樣就得到了lg A－Δ的數據對，共計算得到17011個這樣的數據對。將17011個數據對按照震中距從0～800km的范圍排序，由于數據較多，各臺最大地動位移信息包含的物理意義也相差較大，若不經任何形式的預處理直接進行回歸分析，會使回歸效果不理想。預處理時按照每10km的間隔，從小到大順序分割出一系列震中距段（i＝1～80），把每一段中的 lg Aj－Δj（j＝1～n，n為該震中距段 lg A－Δ的數量）進行平均，得到 lg Ai－Δi（i＝1～80）。采用式（7）對這 80對數據進行“有進有出”的多元回歸分析，回歸分析過程中顯著性水平為0.05，多元回歸分析模型標準誤差為0.07，回歸模型成立。

擬合出來的最大地動位移衰減關系式為

其中，常數項為震源信息項，與衰減特征無關。

根據李氏震級關于零級地震的定義，由lg2800＝3.45，依據最大地動位移衰減關系式，得出地方性震級量規函數的表達式為

3 討論和結論

經前文的討論，可得出如下認識和結論：

（1）由震級殘差統計和最大地動位移衰減關系得出的地方性震級量規函數如圖5所示。

圖5 兩種方法計算出的量規函數與現階段使用的量規函數的對比

經分析圖5可知，現階段使用的量規函數在震中距0～200km的范圍內，地方性震級量規函數偏小，導致近臺震級偏小，由震級的偏差值分布特征也可反映出在0～200km的范圍內震級偏差皆小于平均震級；而由最大地動位移衰減關系得出的地方性震級量規函數在400～800km范圍內比現階段使用的量規函數和由震級殘差得出量規函數都偏?。辉?00～800km范圍內測算出來的震級值偏差值保持在±0.1的范圍內，計算地方性震級時，應盡量選用此范圍內的臺站進行測算。

（2）新疆測震臺網擁有84個數字測震臺站，再加上共享鄰省的17個臺站數據，能較好地覆蓋全區域的地震活動，對于邊境地區和塔里木盆地中部，由于臺站分布偏向一側，故可忽略此區域的震源輻射方位。此外，本文所選用的均為ML≥2.0地震，保證至少5個以上臺站的觀測資料對震中的圍繞，基本消除了震源輻射方位的影響。

（3）由震級殘差統計得出的量規函數中儀器、方位等其他影響因素：通過對選用的新疆測震臺網數字地震記錄資料，對75個臺站的數據進行分析，除了 TCH、WCW、XKR、STZ和JIG這5個臺的震級偏差較大外，有54個臺站的震級偏差都在±0.2的范圍內，占所有臺數的72%。其中JIG臺是臨時臺，臺基條件為砂巖，且風化嚴重；WCW臺在礦區當中，記錄的地震事件常常受到干擾；對TCH、XKR和STZ臺，由于這3個臺計算出的場地響應值較大（劉建明等，2014），臺基對地震波有放大作用，計算出的震級會普遍偏大。

（4）分別由震級殘差統計和地動位移衰減特性得出的量規函數之間有一定的差距，主要表現在所使用的方法的不同。震級殘差統計方法得出的量規函數是在原有量規函數的基礎上一種取平均的思想，而由地動位移衰減特性得出的量規函數是反映真實的地震波的衰減，但由地動位移計算量規函數需依賴于各臺站精確的儀器響應。由于現階段使用的臺站地震儀器響應都存在一定的偏差，故在今后的研究中應兩者結合對比，以便得出更加精確的新疆地區地方性震級量規函數。

致謝：感謝劉瑞豐研究員對本文修改提出的寶貴意見。