郭 偉，劉甜甜 ， 范玲玲， 李貴明

(1.內(nèi)蒙古自治區(qū)地震局烏加河地震臺，內(nèi)蒙古 巴彥淖爾 015331；2.中國地震局地震預(yù)測研究所，北京 100036)

0 引言

地震定位是用各個臺站觀測到的地震波到時來反演地震的空間位置,地震精定位是開展地震學(xué)相關(guān)研究的基石。1912年，地震學(xué)家Geiger率先提出經(jīng)典地震定位方法[1]，但考慮到各個地區(qū)的地質(zhì)構(gòu)造不同、速度模型尚未完善及地震絕對到時拾取誤差較大等因素的影響，經(jīng)典地震定位法未能適用于所有地區(qū)。近年來，國內(nèi)外地震學(xué)者一直致力于提高地震定位精度的研究，傳統(tǒng)的定位方法有S-P到時差、單臺定位法、多臺定位法等。相對定位法主要是應(yīng)用時間之間的相互關(guān)系，常用的有主事件法和雙差定位法。2000年，Waldhauser和Ellsworth提出了雙差地震定位法[2]，應(yīng)用震源之間的相對位置消除速度模型對定位結(jié)果的影響。張海江在此基礎(chǔ)上提出了TomoDD方法[3]，實(shí)現(xiàn)了同時反演地震位置和速度結(jié)構(gòu)；楊智嫻等研究了雙差地震定位法在我國中西部地區(qū)地震精確定位中的應(yīng)用[4]；趙翠萍對1997—2003年新疆伽師震源區(qū)150多個地震進(jìn)行精定位，發(fā)現(xiàn)地震空間分布圖像顯著改善[5]，在空間上表現(xiàn)出較好的叢集；張廣偉等采用雙差定位法，對云南魯?shù)镸S6.5地震序列重新定位，結(jié)果顯示，主震的震源深度與破裂過程顯示的初始破裂深度較為接近，余震序列呈現(xiàn)出近東西—北西向的不對稱共軛狀分布[6]；張帆對岱海斷陷帶及周邊地區(qū)地震事件進(jìn)行雙差重定位后發(fā)現(xiàn)地震在空間分布上更加集中[7]；李翠平等對2017年11月23日重慶武隆MS5.0地震序列進(jìn)行重定位，發(fā)現(xiàn)余震集中分布在主震SW方向[8]；閻春恒等采用雙差定位法對2019年廣西北流MS5.2地震序列進(jìn)行重新定位，發(fā)現(xiàn)地震分布呈北西向優(yōu)勢展布[9]；王恩惠等對鹽城地區(qū)近10年的小震進(jìn)行重新定位，發(fā)現(xiàn)精度大幅提高[10]；李彬[11]、尹占軍[12]分別對內(nèi)蒙古阿拉善左旗5.8級地震序列和內(nèi)蒙古中部地區(qū)符合條件的地震事件進(jìn)行雙差定位，均得到較準(zhǔn)確的震源深度。

雙差定位法在實(shí)際應(yīng)用中要注意以下幾點(diǎn)：不需要主事件，適用于較大的叢集地震空間跨度；適用于震中分布集中的、相鄰地震間距離幾百米甚至更小的地震群；雙差定位算法的主要局限性在于臺站的幾何分布和可用性；近臺對雙差定位結(jié)果有較大的影響。

研究選取2009—2019年發(fā)生在鄂爾多斯塊體附近(北緯34°～42°,東經(jīng)104°～114°)的地震數(shù)據(jù)，重新量取地震波形，采用雙差定位法進(jìn)行重新定位，為進(jìn)一步研究該區(qū)域的地震活動性等奠定基礎(chǔ)。

1 地震地質(zhì)構(gòu)造背景

鄂爾多斯塊體位于華北克拉通塊體內(nèi)部，塊體內(nèi)部構(gòu)造活動微弱，周緣地震活動強(qiáng)烈，歷史地震較多[13]。鄂爾多斯塊體北緣主要位于內(nèi)蒙古中西部地區(qū)，在陰山隆起和鄂爾多斯隆起之間，東西走向[14]，中強(qiáng)地震較為活躍，主要包括河套斷陷系(有呼包凹陷、白彥花凹陷和臨河凹陷，這3個凹陷的北緣分別是大青山山前斷裂、烏拉山山前斷裂和狼山—色爾騰山山前斷裂)和山西斷陷系北部(主要有大同盆地，邊界斷裂有口泉斷裂、六棱山北麓斷裂和衡山北麓斷裂)[15]；鄂爾多斯塊體西緣位于南北地震帶北端，是青藏地塊、阿拉善地塊和鄂爾多斯地塊的交匯部位，新生代以來區(qū)域內(nèi)構(gòu)造活動強(qiáng)烈、南部弧形構(gòu)造帶以擠壓、走滑為特征，北部盆嶺構(gòu)造帶構(gòu)造活動以拉張為主，發(fā)育有規(guī)模較大的張性-扭性正斷層[16]。

2 資料選取

選取2009—2019年發(fā)生在鄂爾多斯塊體附近(北緯34°～42°,東經(jīng)104°～114°)的地震事件為研究對象，采用雙差定位法重新定位。研究區(qū)域內(nèi)臺站和地震觀測報(bào)告中的地震分布如圖1所示。可以看出，地震基本上都分布在鄂爾多斯塊體西北緣，塊體內(nèi)部較少。

圖1 鄂爾多斯塊體西北緣臺站及地震分布Fig.1 Stations and earthquake distribution in the northwest margin of Ordos block

3 原理方法

雙差定位算法是調(diào)整事件對的相對位置，使得觀測值和理論值的到時差最小。對于臺站k，記錄到的第i個地震的體波走時t滿足:

式中:τi為事件的初始時刻;u為慢度;ds為積分路徑長度。

假設(shè)地震事件i、j相距足夠近，這兩個地震事件到同一臺站k的路徑幾乎相同,則:

4 雙差定位結(jié)果及分析

4.1 速度模型

地震學(xué)中經(jīng)典的一維速度模型有J-B模型、PREM、IASP91和AK135模型?？紤]到不同區(qū)域的地質(zhì)構(gòu)造不同，多年來學(xué)者們對不同區(qū)域分別開展了相應(yīng)的速度模型研究。

劉芳[17]建立了內(nèi)蒙古地區(qū)地殼速度模型(簡稱“2015模型”)，趙艷紅[18]在“2015模型”的基礎(chǔ)上分別建立了內(nèi)蒙古西部、中部和東部3個區(qū)域的地殼速度模型。賈昕曄[19]在劉芳、趙艷紅研究的基礎(chǔ)上開展了內(nèi)蒙古中西部地區(qū)一維速度模型檢驗(yàn)，結(jié)果顯示，分區(qū)速度模型更符合內(nèi)蒙古區(qū)域構(gòu)造特征，所以，該文采用的速度模型是賈昕曄的研究結(jié)果(見表1)。其中：H表示地殼厚度，VP表示P波的速度，VS表示S波的速度。

表1 內(nèi)蒙古中部地區(qū)水平層狀速度模型Table 1 Horizontal layered velocity model in central Inner Mongolia

4.2 結(jié)果分析

應(yīng)用雙差定位法，對鄂爾多斯塊體附近的地震進(jìn)行重新定位，定位前后的地震分布如圖2所示?？梢钥闯?，地震基本上都分布在鄂爾多斯塊體西北緣，塊體內(nèi)部較少。由于內(nèi)蒙古地區(qū)地域遼闊且地質(zhì)構(gòu)造背景復(fù)雜，臺站分布不均勻，西部地區(qū)臺站分布較稀疏。雙差定位對于包圍較好的地震叢集有優(yōu)勢，在圖2的基礎(chǔ)上考慮構(gòu)造背景分區(qū)，分別對銀川吉蘭泰斷陷帶(鄂爾多斯塊體西緣)和河套斷陷帶(鄂爾多斯塊體北緣)的地震進(jìn)行重新定位，結(jié)果如第13頁圖3、4所示。雙差定位后的地震分布更集中，西緣和北緣均呈條帶狀分布，朝斷裂帶走向聚集，基本上與斷層走向一致，塊體西緣的地震沿銀川吉蘭泰斷陷帶分布，塊體北緣的地震沿河套斷陷帶分布，顯示出研究區(qū)內(nèi)地震活動與活動構(gòu)造的密切關(guān)系。

圖2 鄂爾多斯塊體雙差定位前后地震分布Fig.2 Seismic distribution before and after double difference location in Ordos block

進(jìn)行雙差定位前后鄂爾多斯塊體及其西緣、北緣的均方根殘差如第13頁表2所示，可以看出，雙差定位后均方根殘差較之前降低了一個數(shù)量級。由于臺站密集，地震分布集中的區(qū)域叢集事件的殘差較小，反演過程組成叢集事件較多，說明雙差定位的精度較目前測震臺網(wǎng)的定位精度有提高。

第13頁表3是雙差定位后震源位置三分向誤差，可以看出，雙差定位后三分向誤差基本均在1 km以內(nèi)。

從第13頁圖5、6、7(a圖均為重新定位前，b圖均為重新定位后)可以看出，重新定位前震源深度分布主要集中在5～20 km，定位后基本呈正態(tài)分布。尹占軍[12]對內(nèi)蒙古中部地區(qū)地震事件進(jìn)行雙差定位，結(jié)果表明，震源深度分布在6～10 km。張帆[20]通過Pn-Pg走時差計(jì)算內(nèi)蒙古地區(qū)地震震源深度發(fā)現(xiàn)，原始地震目錄深度在10 km以內(nèi)，重新測定的深度分布在5～40 km范圍內(nèi)。為此，研究了鄂爾多斯塊體重新定位前后震源深度隨經(jīng)緯度變化的情況，結(jié)果如第14頁圖8所示?？梢钥闯?，重新定位前(見圖8a、8b)震源深度在5～10 km范圍內(nèi)，水平分層現(xiàn)象明顯，這與我國地震目錄的處理方式有關(guān)。目前，測震臺網(wǎng)多采用人工分圖方式，有些震源由于臺站記錄到的震相不足等導(dǎo)致無法獲得準(zhǔn)確的震源深度，從而根據(jù)以往的分圖經(jīng)驗(yàn)人為設(shè)定一個震源深度(采用“地區(qū)優(yōu)勢”深度代替實(shí)際的震源深度[21])，與實(shí)際的深度存在誤差，進(jìn)一步說明對地震進(jìn)行重新定位是必要的。

鄂爾多斯塊體重新定位后(見圖8c、8d)震源深度隨經(jīng)緯度變化情況顯示，重定位后消除了水平分層現(xiàn)象，與張帆[20]通過Pn-Pg走時差重新測定的內(nèi)蒙古地區(qū)地震震源深度分布結(jié)果相吻合。在垂直方向震源呈條帶狀分布，地震優(yōu)勢分布在斷裂帶附近，說明雙差定位算法可以較好地應(yīng)用震源的相對位置，對深度有一定的約束力，能夠更真實(shí)地反映當(dāng)?shù)氐牡刭|(zhì)構(gòu)造情況。

圖3 鄂爾多斯塊體西緣雙差定位結(jié)果Fig.3 Double difference location results of the western margin of Ordos block

圖4 鄂爾多斯塊體北緣雙差定位結(jié)果Fig.4 Double difference location results of the northern margin of Ordos block

表2 雙差定位前后均方根殘差

表3 雙差定位后震源位置三分向誤差

圖5 鄂爾多斯塊體重新定位前后震源深度分布Fig.5 Focal depth distribution before and after relocation of Ordos block

圖6 北緣重新定位前后震源深度分布Fig.6 Distribution of focal depth before and after relocation in the northern margin

圖7 西緣重新定位前后震源深度分布Fig.7 Focal depth distribution berore and after relocation of the western margin

5 結(jié)論與討論

采用雙差定位法對2009—2019年發(fā)生在鄂爾多斯塊體西北緣的地震事件進(jìn)行重新定位，結(jié)果表明：

(1) 雙差定位的精度較高。雙差定位后均方根殘差平方和RMS顯著降低，從精定位前的1.79～2.04降到0.06～0.15，精定位后震源位置三分向誤差基本都在1 km以內(nèi)，垂直定位誤差略高于水平誤差。

圖8 震源深度隨經(jīng)緯度分布圖Fig.8 Distribution of focal depth with latitude and longitude

(2) 震源深度有所改善。消除了重新定位前震源深度水平分層的現(xiàn)象，重新定位后震源深度基本呈正態(tài)分布，在垂直方向上地震呈條帶狀分布。

(3) 震中分布更集中。定位精度明顯提高，鄂爾多斯塊體西緣的地震沿銀川吉蘭泰斷陷帶分布，北緣的地震沿河套斷陷帶分布，與塊體周緣復(fù)雜的地質(zhì)構(gòu)造背景和動力學(xué)過程密切相關(guān)。

經(jīng)過雙差定位的地震精度、震源位置及深度都有改善，為研究該區(qū)域的地震活動性等奠定基礎(chǔ)，震源深度的精確測定在今后的研究中仍需不斷完善。