劉雙慶　陳向軍　王熠熙　王俊菲

摘要：利用新疆和田臺(tái)陣首波Pn段前4s的高頻波動(dòng)形態(tài)，分析了2014年于田M_S7.3主震破裂初期幾個(gè)疑似破裂點(diǎn)的相對(duì)位置，給出兩種可能的初始破裂點(diǎn)運(yùn)動(dòng)學(xué)擴(kuò)展模型。并利用震中距約58km的于田臺(tái)站數(shù)據(jù)進(jìn)一步對(duì)該模型進(jìn)行闡述。分析結(jié)果推測(cè)：此次于田大地震，其起始破裂點(diǎn)開(kāi)始于一個(gè)水平各向異性較強(qiáng)的點(diǎn)位，在之后的3s內(nèi)破裂過(guò)程向各向異性相對(duì)較弱的區(qū)域發(fā)展，且破裂矢量略向南擴(kuò)展。

關(guān)鍵詞：和田臺(tái)陣；Pn段；破裂擴(kuò)展；障礙體

中圖分類號(hào)：P315.72 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1000-0666（2016）02-0295-08

O 引言

2014年2月12日17時(shí)19分48秒，在新疆于田（36.123°N，82.499°E）發(fā)生了M_S7.3大地震。這次地震震中位于距2008年3月21日于田M_S7.3地震的北東東向約90km處，在GoogleEarth地圖上顯示位于硝爾庫(kù)勒高原湖泊西南角，是阿爾金大斷裂的西南緣“爪”狀分枝斷裂交匯部位。該處北東東向的二級(jí)斷裂比較發(fā)育，如阿龍斷裂、庫(kù)牙克斷裂等（國(guó)家地震局《阿爾金活動(dòng)斷裂帶》課題組，1992；楊萬(wàn)志等，2005）。距離震中北西向約328km處的和田臺(tái)陣清晰地記錄了本次地震的波形，距離震中約58km的于田臺(tái)記錄到清晰的P波段波形（S波段記錄限幅）。于田M_S7.3地震發(fā)生后，中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司地震震源深度推進(jìn)研究小組參加了這次地震的室內(nèi)資料應(yīng)急處理工作，并得到一些相關(guān)的研究結(jié)果（陳向軍等，2014；宋秀青等，2014；王俊等，2014），但這些研究都未對(duì)臺(tái)陣記錄的初至的一些高頻成分展開(kāi)分析。實(shí)際上利用臺(tái)陣數(shù)據(jù)分析某些地震破裂過(guò)程時(shí)，由于其地震射線具有相近的傳播路徑，更利于反演震源處信息。因此，本文對(duì)和田臺(tái)陣記錄波形中的前幾秒的高頻成分作進(jìn)一步分析，以期獲得于田M_S7.3地震的一些震源運(yùn)動(dòng)過(guò)程特征。

1 和田臺(tái)陣入射線方位角分析

2008、2014年于田兩次M_S7.3地震，相關(guān)斷裂及和田臺(tái)陣位置見(jiàn)圖1。圖中Moho面深度數(shù)據(jù)來(lái)自滕吉文和曾融生（2002）的結(jié)果，震源機(jī)制來(lái)自萬(wàn)永革等（2010）及中國(guó)地震局地球物理研究所（2014）地震應(yīng)急處置的結(jié)果。其中和田臺(tái)陣位于和田地區(qū)皮山縣，其外形輪廓為圓環(huán)形，由一個(gè)中心臺(tái)和內(nèi)、外環(huán)構(gòu)成。臺(tái)陣各子臺(tái)都布設(shè)在同心圓環(huán)上，孔徑為3km，屬小孔徑臺(tái)陣。內(nèi)環(huán)有3個(gè)子臺(tái)，環(huán)半徑約為600m：外環(huán)有5個(gè)子臺(tái)，環(huán)半徑約為1500m（圖2a）。除中心臺(tái)同時(shí)安裝有一臺(tái)CTS-1E-120s寬頻帶地震計(jì)和CMG-40T-1型地震計(jì)外，其余8個(gè)子臺(tái)均安裝頻帶為0.5～50Hz的CMG-40T-1型地震計(jì)。根據(jù)新疆地震臺(tái)網(wǎng)編目的震中位置以及臺(tái)陣中心臺(tái)的坐標(biāo)，計(jì)算出震中相對(duì)中心臺(tái)的球面方位角為109.64°。為了對(duì)比研究，筆者利用臺(tái)陣波振面方法重新計(jì)算該方位角。計(jì)算過(guò)程分兩步：

（1）利用中心臺(tái)CMG-40T-1型地震計(jì)記錄的Z分量中首波Pn初至及其后10s的數(shù)據(jù)，對(duì)其余8個(gè)子臺(tái)記錄的波形進(jìn)行互相關(guān)，計(jì)算這8個(gè)子臺(tái)相對(duì)中心臺(tái)的超越或延遲時(shí)間。

（2）以中心臺(tái)為參考坐標(biāo)，按0.01°步長(zhǎng)進(jìn)行全方位掃描，計(jì)算各個(gè)子步方位角所對(duì)應(yīng)的入射波射線到8個(gè)子臺(tái)的相對(duì)走時(shí)，并轉(zhuǎn)換成與中心臺(tái)的時(shí)差，與所計(jì)算出的互相關(guān)時(shí)差進(jìn)行對(duì)比，找出二者總時(shí)差最小的入射線方位角。

利用互相關(guān)方法計(jì)算8個(gè)子臺(tái)與中心臺(tái)的時(shí)差時(shí)，分別對(duì)不同波形掃描時(shí)間長(zhǎng)度進(jìn)行了嘗試，顯示掃描時(shí)間長(zhǎng)度為9～11s計(jì)算出的互相關(guān)時(shí)差比較穩(wěn)定，其中掃描時(shí)間為10s的結(jié)果如圖2a所示。由圖可見(jiàn)，在震中方位上的子臺(tái)，其相對(duì)到時(shí)為負(fù)，反方位上為正。0～360°掃描得出方位角為110.46°時(shí)，波射線到達(dá)各臺(tái)的相對(duì)時(shí)差與互相關(guān)波形的時(shí)差相差最小，約0.0493s（圖2b）。圖2a中方括號(hào)內(nèi)的數(shù)字為子臺(tái)站號(hào)，方括號(hào)前數(shù)字為時(shí)間差。圖2b中方括號(hào)內(nèi)的數(shù)字為方位角，方括號(hào)上方數(shù)字為時(shí)間差，不同方位角上的半徑長(zhǎng)度對(duì)應(yīng)于不同波射線時(shí)差與互相關(guān)時(shí)間差的絕對(duì)差值。本文得出的方位角110.46°與109.64°相差0.82°。進(jìn)一步利用Matlab軟件提供的reckon函數(shù)，求出距離中心子臺(tái)327.952km（該值為中心子臺(tái)與新疆地震臺(tái)網(wǎng)所給的震中之間的球面大圓距離），方位角110.46°所對(duì)應(yīng)的震中位置為（36.084°N，82.479°E），比新疆地震臺(tái)網(wǎng)編目的結(jié)果偏西南4.714km，但它滿足地震編目所要求的震中誤差小于5km的一類精度水平。因此利用波形互相關(guān)方法可以較好地獲取地震震中方位角，當(dāng)震中距給定時(shí)，可以得出震中位置。從圖2a還可以看出，震中與中心子臺(tái)的連線幾乎穿過(guò)4、6號(hào)子臺(tái)。7、8、2號(hào)子臺(tái)偏南，3、9、5號(hào)子臺(tái)偏北。

2 臺(tái)陣Pn段的高頻時(shí)間不同步性分析

2.1 破裂過(guò)程穿透的兩種可能介質(zhì)模型

為了對(duì)比和田臺(tái)陣各子臺(tái)記錄的Pn段的一些高頻特征，前文已給出了不同子臺(tái)的時(shí)差以及中心臺(tái)入射線的方位角。將子臺(tái)按到時(shí)的先后順序顯示Z分量（圖3），圖3中x軸時(shí)間跨度近4s。由于臺(tái)陣距離震中約330km，地殼內(nèi)顯著的直達(dá)Pg波尚未到達(dá)，該時(shí)間段波形成分基本為首波。圖3中除了t₁初至的起始點(diǎn)外，在1s左右大周期成分上存在約0.4s的小周期成分。由于數(shù)據(jù)采樣率為100點(diǎn)/s，利用MSDP軟件，將波形放大到可以識(shí)別0.01s的極限顯示精度。根據(jù)t₁、t₂、t₃、t₄（t_i分別為幾個(gè)壓縮波形的起始時(shí)刻）附近的曲線拐點(diǎn)情況、曲線變化形態(tài)及9個(gè)臺(tái)站間波形對(duì)比，逐點(diǎn)分析局部時(shí)間的幅值變化，從而精細(xì)提取各個(gè)子臺(tái)每個(gè)小周期成分的起始時(shí)刻t_i，j（其中i為1～4不同時(shí)刻，j為1～9子臺(tái)號(hào)），并用“+”標(biāo)記在相應(yīng)的波形上。圖4中將子臺(tái)號(hào)標(biāo)注在橫坐標(biāo)軸上，且圖4的臺(tái)號(hào)順序是以垂直入射線（圖2）從西南往東北按距入射線遠(yuǎn)近依次排列，以利于從入射線切線方向分析到時(shí)差異。t₄之后約1s，Pg震相族抵達(dá)臺(tái)陣，波形變得復(fù)雜化。由于Pg初至對(duì)應(yīng)的破裂點(diǎn)早于t₄處Pn對(duì)應(yīng)的破裂點(diǎn)，Pg、Pn地殼速度模型的可靠性對(duì)其分析結(jié)果有較大影響，因此未考慮后續(xù)的細(xì)節(jié)震相。

由圖3可以發(fā)現(xiàn)，不同時(shí)刻點(diǎn)上的“+”有一定的時(shí)間交錯(cuò)性。為此，每個(gè)時(shí)刻點(diǎn)上9個(gè)子臺(tái)都減去子臺(tái)6的到時(shí)，將結(jié)果繪于圖4a。由圖4a可看出，不同震中距具有一定的走時(shí)差異，且后3個(gè)時(shí)刻的高頻成分走時(shí)整體略低。圖4b是在圖4a的基礎(chǔ)上再將t₁時(shí)刻的各個(gè)子臺(tái)時(shí)差對(duì)應(yīng)相減。經(jīng)過(guò)兩次減法，基本扣除了每個(gè)子臺(tái)GPS授時(shí)誤差及臺(tái)陣各子臺(tái)之間距離的影響。圖4b顯示3條曲線幾乎都小于0，而且沿垂直于波射線排列的子臺(tái)所獲取的差值呈“n”型分布，表明后3組波動(dòng)到6號(hào)子臺(tái)的相對(duì)時(shí)間比第1組到達(dá)的相對(duì)時(shí)間減少了。引起這種現(xiàn)象可能有3種因素：①人工提取各個(gè)壓縮波形初至到時(shí)引入了人為誤差；②從震源傳播到臺(tái)陣的射線路徑上存在著橫切路徑的強(qiáng)烈的橫向速度擾動(dòng)，引起一些非均勻相移；③震源破裂過(guò)程中，不同時(shí)刻破裂點(diǎn)所處的介質(zhì)波速參數(shù)存在空間分布差異而引起相移。

對(duì)于第①種因素，筆者認(rèn)為肯定存在拾取誤差，但該誤差應(yīng)呈現(xiàn)非相關(guān)特征，不會(huì)出現(xiàn)“∩”型分布。對(duì)比圖2a和圖3，發(fā)現(xiàn)人工提取的初始時(shí)間和互相關(guān)方法給出的序列并不一致，特別是4、5號(hào)子臺(tái)的順序。這種差異可能是做互相關(guān)時(shí)使用了10s的時(shí)間，在這個(gè)時(shí)段中震源破裂點(diǎn)已出現(xiàn)位置變化，5號(hào)子臺(tái)更早接收到后續(xù)破裂點(diǎn)的震動(dòng)。2、3號(hào)子臺(tái)也有類似的現(xiàn)象。而發(fā)生這個(gè)現(xiàn)象的臺(tái)站主要分布于相對(duì)中心臺(tái)與震中連線近似對(duì)稱的位置。對(duì)于第②種因素，在2.2節(jié)進(jìn)一步分析。

對(duì)于第③種因素，給出以下模型。因?yàn)楹吞锱_(tái)陣屬小孔徑臺(tái)陣，直徑約3km，震中距328km，臺(tái)陣對(duì)震源的張角很小（約0.524。）。射線在震源體之外的傳播過(guò)程中幾乎穿透相同的介質(zhì)，會(huì)引起后3組時(shí)差變小，最有可能源于震源破裂點(diǎn)附近的介質(zhì)參數(shù)在空間分布上的變化。假定3種情況：（1）震源破裂點(diǎn)下部的各向異性差異減少；（2）震源破裂點(diǎn)上部的各向異性差異減少；（3）震源破裂點(diǎn)上下部的各向異性差異不變。

圖5中，保持模型c參數(shù)不變，以作對(duì)比。破裂過(guò)程都由黑點(diǎn)向白點(diǎn)擴(kuò)展；每方格的高度為1，方格上的數(shù)字對(duì)應(yīng)速度比值；黑色箭頭為Pn傳播方向。模型b中，破裂點(diǎn)從下往上破裂，而首波都是朝下傳播，因此這里的箭頭使用虛線，以示區(qū)別。3種模型的Pn穿透不同參數(shù)的兩個(gè)方格厚度后的時(shí)差分別為：模型a：1.3409；模型b：1.8333；模型c：2。

因此只要在破裂過(guò)程中，破裂點(diǎn)擴(kuò)展的后期介質(zhì)參數(shù)比前期的參數(shù)（這里指P波波速）更趨于相近時(shí)，則從不同時(shí)刻傳出的首波震動(dòng)到達(dá)臺(tái)陣各子臺(tái)的時(shí)間將趨近。如果破裂后期的介質(zhì)參數(shù)差異大于前期，那將大于模型c的值，不滿足圖4b中所顯示的現(xiàn)象。由于在各向同性線彈性模型中，介質(zhì)波速與介質(zhì)參數(shù)關(guān)系式為V_p=？其中，λ、μ為拉梅常數(shù)，分別反映壓縮膨脹效應(yīng)和剪切形變效應(yīng)；ρ為介質(zhì)密度，對(duì)于地殼介質(zhì)，λ、μ的垂向變化一般較ρ明顯。因此推測(cè)于田M_S7.3地震很可能符合（1）、（2）其中的一種情況，破裂點(diǎn)都有由破裂面上各向異性較強(qiáng)的點(diǎn)位向各向異性趨同的點(diǎn)位擴(kuò)展的過(guò)程，在此過(guò)程中出現(xiàn)了疑似障礙體或凹凸體。障礙體指高臨界強(qiáng)度區(qū)，也可以是斷層面上的一些二級(jí)凹凸面（Aki，1979），凹凸體指應(yīng)力集中區(qū)（宮武隆，1993；Dragoni，Piombo，1998）的特征。在筆者研究的這4s時(shí)間內(nèi)，破裂點(diǎn)由各向異性較強(qiáng)（如水平差異導(dǎo)致的垂向組合分布）的點(diǎn)位向較弱點(diǎn)位過(guò)渡，每個(gè)破裂點(diǎn)的子震源在這個(gè)時(shí)間段內(nèi)依次產(chǎn)生并傳播首波，但后續(xù)破裂點(diǎn)的首波所經(jīng)過(guò)的射線上的介質(zhì)差異減弱，從而出現(xiàn)后續(xù)子破裂的首波走時(shí)更接近。實(shí)際上，在精細(xì)的破裂過(guò)程反演計(jì)算中，障礙體破裂過(guò)程的滑動(dòng)幅度與其斷層相應(yīng)位置處的剛度有一定的模型關(guān)系（何玉梅，鄭天愉，1998）。這種反演方法上對(duì)剛度的差異性設(shè)置，與圖5的速度垂向異性簡(jiǎn)化模型都屬于刻畫(huà)震源體內(nèi)部的各向異性特征。至于于田主震破裂方式是從上往下，還是從下往上，筆者沒(méi)有搜集到更多的同精度的資料供分析。由于時(shí)差呈“∩”型分布，不是單斜曲線，不支持該短暫破裂過(guò)程只朝單側(cè)進(jìn)行的情況，因此筆者著重討論由水平分異引起的垂向介質(zhì)參數(shù)差異。而張廣偉等（2014）反演的結(jié)果也顯示此次大地震破裂面比較集中，水平兩側(cè)擴(kuò)展的尺度不明顯。

2.2 關(guān)于震源區(qū)以外速度界面影響的分析

進(jìn)一步討論2.1節(jié)中第②種因素的影響。不只傳播路徑上的速度間斷面可以衍生震相，地震計(jì)的瞬態(tài)反應(yīng)如果沒(méi)有有效的阻尼作用，也將出現(xiàn)拖尾現(xiàn)象。對(duì)于地震儀器的因素，由于同臺(tái)的瞬態(tài)響應(yīng)不變，因此如果存在后續(xù)“拖尾”震相，這些震相應(yīng)該是等間隔的，不會(huì)出現(xiàn)“∩”型分布。而對(duì)于速度界面的影響，引入震中距58km的于田（YUT）臺(tái)站波形進(jìn)行分析。同時(shí)為了進(jìn)一步對(duì)比，還引入了一個(gè)震中距與之相近（56km）、位錯(cuò)為走滑兼正斷性質(zhì)的、且位于初動(dòng)向上的小震波形。該波形為北京南山村（NSC）臺(tái)的記錄的2012年6月18日天津?qū)氎鍹_S4.0地震的波形（震源機(jī)制解：走向142°，傾角87°，滑動(dòng)角-15°）。由于天津有近50個(gè)強(qiáng)震加速度計(jì)獲取了相應(yīng)的烈度分布，可以確定該地震為北西向破裂，所以與之共軛的北東向震源機(jī)制解舍棄。使用這個(gè)對(duì)比地震同時(shí)還考慮到4級(jí)地震破裂過(guò)程簡(jiǎn)單。

圖6中將原始記錄的EW、NS分量旋轉(zhuǎn)成T、R分量，以確認(rèn)于田臺(tái)站記錄的直達(dá)剪切波sg初至。對(duì)比圖6b、f，NSC臺(tái)記錄的Sg-Pg比YUT的臺(tái)記錄的略小，但都接近8s，反映出二者震中距的大小。由于寶坻Ms4.0地震破裂尺度小，在圖6e～g中可以明顯看出，在Pg之后，Sg之前，波形幅度明顯比Pg前2個(gè)周期的振幅小。而且在這個(gè)震中距，反射震相P^*、P₁₁尚不發(fā)育，因此Pg振幅占絕對(duì)優(yōu)勢(shì)。而于田M_S7.3地震破裂面尺度可達(dá)50×15km（陳運(yùn)泰，顧浩鼎，2004），若以0.8倍的橫波速度破裂，破裂10km尺度也需約3.125s（破裂速度假定為3.2km/s）。因此筆者在局部放大圖6d中發(fā)現(xiàn)，波幅強(qiáng)度呈增長(zhǎng)趨勢(shì)，特別是圖6c中的Px處出現(xiàn)了一個(gè)更顯著的介質(zhì)破裂。由此推測(cè)，即使在于田地震破裂的前3s內(nèi)也存在多個(gè)級(jí)聯(lián)子破裂，與永善大震破裂過(guò)程類似（曾融生，師潔珊，1978）。而位于震中距328km處的和田臺(tái)陣記錄的Pn波形中捕捉了這些級(jí)聯(lián)破裂的某些初動(dòng)曲線拐點(diǎn)。

由于一個(gè)大地震破裂的過(guò)程是不均勻的，每次子破裂點(diǎn)的位置在橫向和縱向上都不一樣。破裂過(guò)程中震源區(qū)介質(zhì)的性質(zhì)（彈性、均勻性、各向異性）都將發(fā)生變化，關(guān)于臺(tái)陣子臺(tái)時(shí)差現(xiàn)象的分析，除了上述3種情況外，是否還有另外的解釋？這里引入兩個(gè)概念：一是震源區(qū)靜態(tài)的各向異性分布；二是地震發(fā)生過(guò)程伴隨時(shí)間變化的局部性質(zhì)瞬態(tài)改變。第一種情況為圖5所示情形。對(duì)于隨著破裂過(guò)程而產(chǎn)生的介質(zhì)性質(zhì)瞬態(tài)變化，如果破裂點(diǎn)外的介質(zhì)均勻或水平層狀均勻，則這個(gè)子破裂產(chǎn)生的波源傳到臺(tái)陣各子臺(tái)的時(shí)間變化較少，難以觀測(cè)，或呈單斜曲線。因?yàn)閺钠湔鹪礄C(jī)制解可知，斷層面傾角約82°，而本文研究的子破裂間隔約1s，如果是位錯(cuò)在入射面上的水平破裂分量引發(fā)的波速擾動(dòng)所致（水平層狀均勻，垂直分量不產(chǎn)生此效應(yīng)），當(dāng)破裂速度取3.2km/s，震源體附近地殼P波速度6.2km/s，速度擾動(dòng)為10%時(shí)，兩次破裂的各子臺(tái)對(duì)應(yīng)時(shí)差為0.0018s。這個(gè)量值在采樣率為100點(diǎn)/s觀察不到。

而對(duì)于水平成層介質(zhì)，垂直入射面的位錯(cuò)破裂分量引發(fā)的介質(zhì)P波速度變化，對(duì)各子臺(tái)的相對(duì)時(shí)差無(wú)影響，即沿入射面法向的介質(zhì)瞬態(tài)變化對(duì)各子臺(tái)相對(duì)走時(shí)無(wú)影響（處于同一平面波陣面同等相移不影響各子臺(tái)的相對(duì)到時(shí)）。對(duì)于水平單向破裂，則需進(jìn)一步分析。根據(jù)圖1斷裂走向及震源機(jī)制解，取破裂面走向239°，且假定本文研究的短時(shí)間破裂過(guò)程中，破裂面走向不變，根據(jù)震中相對(duì)各子臺(tái)的方位角，求出圖7a中各子臺(tái)的θ，進(jìn)而可以利用破裂長(zhǎng)度d₁、d'₁，計(jì)算出P₂、P₃到子臺(tái)i的走時(shí)與原路徑L₀走時(shí)差。先假定破裂由P₁向西南破裂10km到達(dá)P₂，破裂過(guò)程介質(zhì)力學(xué)性質(zhì)發(fā)生瞬時(shí)改變，可以產(chǎn)生10%的P波速度擾動(dòng)，則以0.9倍的P波速度穿過(guò)L₀-L₁的時(shí)間差，如圖7b的Linel曲線所示（減去了最小變化的7號(hào)子臺(tái)走時(shí)差，括號(hào)里的數(shù)值是7號(hào)的相對(duì)走時(shí)差），以1.1倍P波速度穿過(guò)L₀-L₁的時(shí)間差，如圖7b的Line2曲線所示。而由P₁向東北破裂10km到達(dá)P₃的兩個(gè)速度計(jì)算的走時(shí)差曲線分別為L(zhǎng)ine3、Line4。計(jì)算表明，在10%速度擾動(dòng)的情況下，要使該臺(tái)陣子臺(tái)的相對(duì)走時(shí)差達(dá)到0.01s的識(shí)別水平，破裂長(zhǎng)度應(yīng)不小于10km。因此對(duì)于本文研究的3～4s的波形長(zhǎng)度，即使存在單側(cè)破裂，從P波信息里也僅僅能略有反映，但不明顯。不過(guò)對(duì)于單側(cè)破裂，各子臺(tái)的相對(duì)走時(shí)是單斜型曲線。從Line1、Line2等后括號(hào)內(nèi)數(shù)字對(duì)比還可顯示，向西南破裂引起的相對(duì)走時(shí)差略小于向東北破裂，這是因?yàn)榕_(tái)陣位于震中西北，而破裂走向?yàn)楸睎|，向東南破裂時(shí)，新破裂點(diǎn)到各子臺(tái)的距離都減小了（圖7a）。因此，對(duì)比圖4b與圖7b的走時(shí)差量級(jí)，筆者認(rèn)為震源區(qū)應(yīng)存在“靜態(tài)”的各向異性（圖5a、b），才能出現(xiàn)較符合子臺(tái)相對(duì)走時(shí)差減少的特征。如果彈性模量值較大區(qū)與較小區(qū)的P波速度差異為10%，則要產(chǎn)生本文可觀測(cè)時(shí)差（如0.04s）的彈性異常體尺度約1.23km（與入射面平行的方向）。如果彈性模量值較大區(qū)與較小區(qū)的P波速度差異為2%，則產(chǎn)生本文可觀測(cè)時(shí)差的彈性異常體尺度約6.19km（條件同上）。對(duì)于于田地震，當(dāng)震源體附近存在各向異性時(shí)，產(chǎn)生同等速度擾動(dòng)的各向異性體尺度要比單側(cè)破裂的長(zhǎng)度小若干倍。本文觀測(cè)到的時(shí)差量級(jí)更傾向于震源體附近存在各向異性體的結(jié)論。

3 結(jié)論與討論

本文利用和田臺(tái)陣的地震資料對(duì)2014年于田M_S7.3地震前4s的破裂過(guò)程進(jìn)行了細(xì)致分析。利用中心子臺(tái)與震中的連線確定臺(tái)陣中各個(gè)子臺(tái)的相對(duì)位置。在此基礎(chǔ)上對(duì)子臺(tái)波形進(jìn)行特定排序，并提取了不同時(shí)刻點(diǎn)上的成組高頻壓縮波初至。這些初至走時(shí)之間的相互差異反映了一些震源破裂點(diǎn)處的介質(zhì)特性。為進(jìn)一步討論這些特性的可能分布，利用寶坻地震的NSC臺(tái)資料與于田地震的于田臺(tái)資料進(jìn)行對(duì)比，結(jié)果反映和田臺(tái)陣記錄的Pn段內(nèi)確實(shí)存在多級(jí)級(jí)聯(lián)破裂現(xiàn)象。如果早期破裂點(diǎn)處的各向異性差異大于后續(xù)破裂點(diǎn)（只分析了3s）的差異，則可以更好解釋本文提取到的“∩”型時(shí)間分布現(xiàn)象。這也說(shuō)明了和田臺(tái)陣記錄到了這個(gè)大地震最初破裂區(qū)的一些級(jí)聯(lián)破裂的凹凸體或障礙體特征。若利用圖4b中的3、9、5子臺(tái)時(shí)差比其他子臺(tái)更小的細(xì)微差異，則可進(jìn)一步推測(cè)后續(xù)破裂點(diǎn)將向南偏移，使得位于入射線北部的子臺(tái)經(jīng)過(guò)更長(zhǎng)路徑的各向異性較弱的區(qū)域，從而降低相對(duì)時(shí)差。這與中國(guó)地震局地球物理研究所應(yīng)急處置給出的破裂過(guò)程反演分布中往南往深部擴(kuò)展的特征相吻合，同時(shí)也與張廣偉等（2014）復(fù)雜反演的往南破裂結(jié)果吻合。

由于和田臺(tái)陣距離震中還是較遠(yuǎn)，一些近場(chǎng)高頻成分未能有效記錄。對(duì)高頻成分的人工提取位置也不能認(rèn)定其必然很可靠。因此本文的結(jié)論還需要收集更多的近場(chǎng)臺(tái)陣資料或相關(guān)類似資料進(jìn)行判定。但利用臺(tái)陣獲取大震破裂過(guò)程的細(xì)微特征，這項(xiàng)工作是非常值得去研究的。

非常感謝在2014年于田地震應(yīng)急工作中，新疆維吾爾自治區(qū)地震局蔚曉利研究員及唐明帥博士給予的幫助。同時(shí)感謝審稿老師對(duì)初稿的修改建議。