秦滿忠 李順成 劉旭宙 張淑珍 張元生，

1）中國蘭州730000中國地震局蘭州地震研究所

2）中國蘭州730000中國地震局地震預(yù)測研究所蘭州科技創(chuàng)新基地

3）中國北京100029中國地震局地質(zhì)研究所地震動力學(xué)國家重點實驗室

引言

核幔界面衍射波Pdiff是一個極遠(yuǎn)震弱震相，單脈沖型，振幅較小，僅為PP的1／10—1／5，一般在當(dāng)MS≥7.0且干擾背景較小時才能夠在長周期地震記錄中得以分辨（李克，2007）.丹麥地震學(xué)家萊曼1936年首次在地震圖上識別出核幔界面衍射波Pdiff震相（Lehmann，1936）（圖1），并討論了震中距為100°—140°的核幔界面上的衍射波Pdiff問題.Rial和Cormier定量地研究了地球?qū)︴劈c（178°≤Δ≤180°）附近的Pdiff，PKIKP及PP等震相的基本特征，發(fā)現(xiàn)對跖點附近的Pdiff和PP震相具有聚焦效應(yīng)（振幅放大現(xiàn)象），而 PKIKP 沒有聚焦效應(yīng)（Rial，1978；Rial，Cormier，1980）.

Pdiff震相在全球長周期地震數(shù)據(jù)疊加的觀測走時曲線中的最大記錄震中距約150°（Astiz et al，1996）.唐燕娟（1997）使用1980—1995年全國地震臺網(wǎng)極遠(yuǎn)震資料得到Pdiff震相記錄范圍在100°—158°附近.據(jù)趙榮國（1987）觀測研究表明，使用長周期763型地震儀可記錄到弱幅度的Pdiff震相，其最大震中距為175.6°（成都臺）.許健生等（2012）使用寬頻帶數(shù)字地震記錄對蘭州地震臺記錄的極遠(yuǎn)震進行了分析，得到Pdiff震相的記錄范圍為114.2°—175.7°.

2010年2月27日在智利中南部發(fā)生MS8.8特大地震，甘肅數(shù)字測震臺網(wǎng)和甘東南野外觀測流動臺陣記錄到該地震豐富的（210個臺站）高質(zhì)量觀測波形，其記錄震中距范圍在170°—180°之間.本文通過讀取智利地震初至震相的觀測走時，計算該震相的觀測慢度，并與Pdiff震相理論慢度相比較來判定初至震相是否為核幔界面衍射波震相，進一步探討Pdiff震相的記錄特征、最大記錄震中距以及不同方位傳播的Pdiff震相觀測走時與核幔界面復(fù)雜形態(tài)之間的關(guān)系.

圖1 核幔界面衍射波Pdiff震相和PKIKP震相射線路徑傳播示意圖Fig.1 Schematic diagram of ray paths of seismic phases Pdiffand PKIKP

1 臺網(wǎng)資料

甘肅數(shù)字測震臺網(wǎng)和甘東南野外觀測流動臺陣均布設(shè)于青藏高原東北緣.甘肅“十五”數(shù)字測震臺網(wǎng)自2008年6月正式運行，由蘭州、高臺、安西、嘉峪關(guān)和天水等5個有人值守的國家數(shù)字地震臺和39個區(qū)域遙測數(shù)字地震臺組成（馮建剛等，2012），目前擁有包括鄰省（寧夏、陜西、四川、青海、內(nèi)蒙）在內(nèi)的76個數(shù)字測震臺站；甘東南野外觀測流動臺陣由中國地震局地質(zhì)研究所和中國地震局蘭州地震研究所共同合作在甘東南地區(qū)勘址架設(shè)，共設(shè)有7條測線，150個流動臺站，平均臺間距約10km（圖2）.使用REF TEK-130型數(shù)據(jù)采集器和Guralp公司生產(chǎn)的CMG-3ESPC寬頻帶地震計，其頻帶范圍為60—0.02s.甘東南野外觀測流動臺陣運行期間（2009年11月—2011年12月）積累了連續(xù)、可靠、高質(zhì)量的地震觀測波形數(shù)據(jù)，為深入研究地球內(nèi)部物理結(jié)構(gòu)提供了寶貴的資料.

圖2 甘肅數(shù)字測震臺網(wǎng)和甘東南野外觀測流動臺陣分布圖Fig.2 Distribution of the Gansu Digital Seismic Network（blue triangles）and southeastern Gansu temporary observation array（black triangles），where the red square is the position of the Lanzhou small aperture seismic array

2 數(shù)據(jù)處理

本文選取甘肅數(shù)字測震臺網(wǎng)和甘東南野外觀測流動臺陣共同記錄的2010年2月27日發(fā)生在智利中南部 MS8.8特大地震（36.12°S、72.90°W，h＝23km）數(shù)字觀測波形數(shù)據(jù).

首先截取智利特大地震數(shù)字觀測波形數(shù)據(jù)，對觀測波形數(shù)據(jù)進行去傾斜、去平均處理，然后作濾波處理，去掉記錄畸形、信噪比差的觀測波形，最終選取震中距為170°—180°的155條垂直分量的觀測波形數(shù)據(jù).為了易于讀取初至震相走時，更好地顯示震相（Pdiff和PKIKP）和觀測波形數(shù)據(jù)，我們使用了30s的低通濾波（圖3）.

圖3 甘肅數(shù)字測震臺網(wǎng)和甘東南野外觀測流動臺陣觀測到的智利中南部MS8.8地震波形Fig.3 Waveforms of a MS＝8.8event that occurred in south-central Chile on February 27，2010，recorded by the Gansu Digital Seismic Network and the southeastern Gansu temporary observation array

結(jié)合IASP91理論走時曲線，從圖3中1 215s附近很容易識別出極遠(yuǎn)震PKIKP震相.同時，在PKIKP震相之前存在一個較弱的初至震相，具有單脈沖型、振幅較小，周期約42s.Pdiff震相在IASP91理論模型的最大記錄震中距為159°，而本文中臺網(wǎng)記錄初至震相的震中距范圍為170°—180°.為了確定該初至震相是否為核幔界面衍射波震相Pdiff，我們分別讀取了初至震相的起始觀測走時和峰值觀測走時，并分別與對應(yīng)的震中距作了線性擬合（圖4）：

式中，t1為初至震相起始觀測走時，t2為初至震相峰值觀測走時，t3為Pdiff震相IASP91理論走時，單位為s；Δ為震中距，單位為度；相關(guān)系數(shù)分別為R1＝0.962，R2＝0.975，R3＝1.000.通過線性擬合初至震相起始觀測走時t1和峰值觀測走時t2，得到對應(yīng)的平均觀測慢度值分別為4.104s／°和3.980s／°，與Pdiff震相的理論慢度4.439s／°基本一致.由此可以確定，甘肅數(shù)字測震臺網(wǎng)和甘東南野外觀測流動臺陣記錄的初至震相為核幔界面衍射波Pdiff震相，其最大記錄震中距為179.7°.圖3中Pdiff震相紅色地震觀測波形的峰值點隨震中距的增大存在明顯的觀測走時超前和滯后現(xiàn)象，這與核幔界面凹凸不平的復(fù)雜形態(tài)有關(guān).

圖4 初至震相（Pdiff）起始觀測走時（R1，t1）和峰值觀測走時（R2，t2）的線性擬合Fig.4 Linear regression of origin and peak observational travel times for the first arrival phase（Pdiff）.R1，t1and R2，t2are correlation coefficients and linear regression equations of the origin observational travel time and the peak observational travel time，respectively

3 討論與結(jié)論

本文選取了2010年2月27日智利MS8.8特大地震，保證其在全球所有地震臺均具有清晰、良好的地震觀測波形記錄特征；甘肅數(shù)字測震臺網(wǎng)和甘東南野外觀測流動臺陣擁有大量、密集的地震臺站（210個），對此次特大地震有連續(xù)、可靠、高質(zhì)量的地震觀測波形，初至震相在臺網(wǎng)記錄中的震中距為170°—180°.

為了確定初至震相Pdiff起始觀測走時的可靠性，對智利特大地震甘肅臺網(wǎng)SBC臺（Δ＝170.54°）垂直分量波形分別選用了2—10s、10—30s、30—100s帶通濾波和30s低通濾波，并與原始觀測波形進行了比較（圖5）.結(jié)果表明，通過30s低通和30—100s帶通濾波處理后的觀測波形與原始觀測波形在頻率、相位和記錄波形特征上具有很好的一致性，30s低通濾波結(jié)果的初至震相Pdiff起始更為尖銳，更有利于讀取初至震相Pdiff起始觀測到時，其走時為1 151.1s.

許建生等（2013）利用中國數(shù)字地震臺網(wǎng)11個臺的寬頻帶長周期數(shù)字記錄對極遠(yuǎn)震（震源深度≤100km）震相Pdiff走時與震中距關(guān)系作了線性擬合，我們利用其擬合關(guān)系式計算得到SBC臺Pdiff震相觀測走時為1 147.3s，與我們讀取的起始觀測走時僅差3.8s.由于本次智利地震震源深度為23km，這種差別應(yīng)是震源深度造成的，因為震源深度越淺，走時越大.由本文式（3）計算的SBC臺Pdiff震相IASP91理論走時為1 136.1s，與許健生等（2013）線性擬合關(guān)系計算的走時1 147.3s也相差11.2s.因此，本文中Pdiff震相的起始觀測到時應(yīng)該是可靠的.

通過分別進行線性擬合初至震相的起始觀測走時、峰值觀測走時與震中距的關(guān)系，其相關(guān)系數(shù)分別為0.961 9和0.975 2，獲得的平均觀測慢度值分別為4.104s／°和3.980s／°，與Pdiff震相IASP91理論模型的理論慢度值4.439s／°基本一致，故可確定甘肅數(shù)字測震臺網(wǎng)和甘東南野外觀測流動臺陣記錄的初至震相為核幔界面衍射波Pdiff震相.由于該震相傳播距離很大，高頻成分在傳播過程中逐漸被吸收，其周期達到42s，平均幅度不到PP震相的百分之一.當(dāng)Δ＞179°時，由于PKIKP震相與Pdiff震相的相互疊加作用，Pdiff震相表現(xiàn)為約1／2的孤峰（圖3）.本文中Pdiff震相最大記錄震中距為179.7°，為目前該震相的最大記錄震中距，記錄臺站位于甘肅省慶陽市西南部的鎮(zhèn)原臺（ZYT）.

以PKIKP震相最大幅度（1 200—1 230s）歸一化觀測波形，與Pdiff震相最大幅度進行比較.結(jié)果表明，Pdiff與PKIKP的相對最大振幅比在震中距為170°—177°時變化不大，在震中距為177°—180°時，最大振幅比隨震中距的增大而增大（圖6）.這可能與Pdiff震相在地球?qū)︴劈c的聚焦效應(yīng)有關(guān)（Rial，1978；Rial，Cormier，1980）.

圖5 原始觀測波形與不同頻率濾波后的觀測波形的比較Fig.5 Comparison of the original observation waveform （thick line）with four filtered observation waveforms（thin lines）

圖6 Pdiff震相與PKIKP震相的最大振幅比Fig.6 The maximum amplitude ratio of Pdiffto PKIKP

表1 不同方位角Pdiff震相起始觀測走時tPdiff與震中距Δ（單位：度）的線性擬合Table 1 Linear regression between origin observational travel time（tPdiff）and its epicentral distanceΔ （in unit of degree）for Pdiffphases with different azimuth range

表2 不同方位角Pdiff震相峰值觀測走時tpeak與震中距Δ（單位：度）的線性擬合Table 2 Linear regression between peak observational travel time（tpeak）and its epicentral distanceΔ （in unit of degree）for Pdiffphase with different azimuth range

甘東南野外觀測流動臺陣關(guān)于智利MS8.8特大地震的記錄為我們提供了頗為寶貴的研究資料，結(jié)合甘肅數(shù)字測震臺網(wǎng)觀測記錄波形本文研究了此次地震不同方位（190°—340°）、震中距在175°—180°之間的148條地震記錄波形對Pdiff震相起始觀測走時（表1）和峰值觀測走時（表2）的影響.將觀測走時按照方位角（190°—340°）劃分為7檔（表1，表2），分別作了不同方位角走時曲線的線性擬合，如圖7所示.可以看出，隨著方位角的逐漸增大，Pdiff震相起始觀測走時和峰值觀測走時均存在較明顯的滯后現(xiàn)象.核幔界面是物理變化極其復(fù)雜的地帶，是一個凹凸不平的界面（Wysession et al，1992；Arthur，John，1993；Sze，van der Hilst，2003），下地幔普遍存在超低速區(qū)，其P波速度和S波速度分別驟降約10%和30%（Mori，Helmberger，1995；Helmberger et al，1998；McNamara et al，2010），超低速區(qū)之上D″區(qū)速度結(jié)構(gòu)都表現(xiàn)為低速，且厚度明顯加厚（唐群署，李麗紅，2006）.相關(guān)研究結(jié)果表明，富含鐵的氧化物可能分離結(jié)晶形成原始海洋巖漿停留在核幔界面之上，并在該界面之上形成高度約100km的穹狀結(jié)構(gòu)，這些少量富含鐵的氧化物的存在會明顯降低地震波速度（Labrosse et al，2007；Sun et al，2013）.因此，方位角的增大所引起的Pdiff震相觀測走時差異可能與核幔界面的復(fù)雜形態(tài)有關(guān).

與IASP91理論模型Pdiff震相走時曲線的延伸部分相比，本文線性擬合的Pdiff震相起始觀測走時曲線存在明顯的滯后現(xiàn)象，滯后時間約15s（圖7）.而利用蘭州小孔徑臺陣（圖2）資料疊加觀測走時曲線結(jié)果（秦滿忠等，2014）表明，遠(yuǎn)震P震相觀測走時滯后于理論P震相走時，并且這種差異在臺陣的東部和西部明顯不同，西部的觀測走時滯后現(xiàn)象更為明顯，其中蘭州小孔徑臺陣記錄的2003年1月22日墨西哥極遠(yuǎn)震（18.77°N、104.10°W，h＝24km，Δ＝119.5°，MS7.6，Az＝30.5°）Pdiff震相起始觀測走時（916.4s）相對于理論觀測走時（909.5s）滯后6.9s.而智利MS8.8特大地震震中位于甘肅數(shù)字測震臺網(wǎng)西部，其記錄震中距（175°—180°）遠(yuǎn)大于墨西哥地震，且在核幔界面衍射傳播的距離更長（75°左右）.因此，Pdiff震相起始觀測走時的滯后不僅與地震波在低速的上、下地幔介質(zhì)中傳播有關(guān)（Soldati et al，2006；秦滿忠等，2014），同時與核幔界面存在超低速區(qū)及凹凸不平的復(fù)雜形態(tài)有關(guān).

圖7 不同方位角Pdiff震相起始觀測走時曲線和峰值觀測走時曲線與IASP91理論走時曲線延伸部分對比Fig.7 Comparison of Pdifforigin and peak observed travel time curves in different azimuth ranges with theoretical travel time curve of the IASP91model

馮建剛，蔣長勝，韓立波，陳繼鋒.2012.甘肅測震臺網(wǎng)監(jiān)測能力及地震目錄完整性分析［J］.地震學(xué)報，34（5）：646－－658.

Feng J G，Jiang C S，Han L B，Chen J F.2012.Analysis on the monitoring capability of seismic networks and completeness of earthquake catalogues in Gansu region［J］.Acta Seismologica Sinica，34（5）：646－－658（in Chinese）.

李克.2007.地震學(xué)與地震觀測［M］.北京：地震出版社：172.

Li K.2007.Seismology and Earthquake Observation［M］.Beijing：Seismological Press：172（in Chinese）.

秦滿忠，沈旭章，張元生，劉旭宙.2014.利用蘭州小孔徑地震臺陣資料疊加觀測走時曲線［J］.地震學(xué)報，36（1）：59－－69.

Qin M Z，Shen X Z，Zhang Y S，Liu X Z.2014.Observed travel-time curves by stacking records from Lanzhou small aperture seismic array［J］.Acta Seismologica Sinica，36（1）：59－－69（in Chinese）.

唐群署，李麗紅.2006.核幔邊界D″區(qū)的地震學(xué)研究進展［J］.地學(xué)前緣，13（2）：213－－223.

Tang Q S，Li L H.2006.The Earth’s lowermost mantle and its seismological research progress［J］.Earth Science Frontiers，13（2）：213－－223（in Chinese）.

唐燕娟.1997.對P波衍射波Pdif震相記錄范圍的討論［J］.地震地磁觀測與研究，18（2）：32－－36.

Tang Y J.1997.A discussion on recorded distance for diffraction wave Pdif of wave P［J］.Seismological and Geomagnetic Observation and Research，18（2）：32－－36（in Chinese）.

許健生，隗永剛，張淑珍.2012.Pdif震相在CDSN蘭州地震臺的記錄特征［J］.西北地震學(xué)報，34（3）：239－－244.

Xu J S，Wei Y G，Zhang S Z.2012.The recording characters of Pdif phase in seismic wave at CDSN Lanzhou seismic station［J］.Northwestern Seismological Journal，34（3）：239－－244（in Chinese）.

許健生，隗永剛，張淑珍.2013.Pdif波在中國CDSN臺站的運動學(xué)特征［J］.地震工程學(xué)報，35（3）：610－－617.

Xu J S，Wei Y G，Zhang S Z.2013.The kinematic characteristics of Pdif-waves at CDSN stations［J］.China Earthquake Engineering Journal，35（3）：610－－617（in Chinese）.

趙榮國.1987.中國763長周期臺網(wǎng)觀測到Δ＝175.6°的幔－核界面衍射波［J］.國際地震動態(tài)，（3）：14－－16.

Zhao R G.1987.Waves diffracted at the mantle-core boundary observed by type-763long－period seismograph network of China［J］.Recent Developments in World Seismology，（3）：14－－16（in Chinese）.

Arthur R，John W.1993.Inference of core-mantle boundary topography from ISC PcP and PKP traveltimes［J］.Geophys J Int，115（3）：991－－1011.

Astiz L，Earle P，Shearer P.1996.Global stacking of broadband seismograms［J］.Seismol Res Lett，67（4）：8－－18.

Helmberger D V，Wen L，Ding X.1998.Seismic evidence that the source of the Iceland hotspot lies at the core-mantle boundary［J］.Nature，396（6708）：251－－255.

Labrosse S，Hernlund J W，Coltice N.2007.A crystallizing dense magma ocean at the base of the Earth’s mantle［J］.Nature，450（7171）：866－－869.

McNamara A K，Garnero E J，Rost S.2010.Tracking deep mantle reservoirs with ultra-low velocity zones［J］.Earth Planet Sci Lett，299（1／2）：1－－9.

Mori J，Helmberger D V.1995.Localized boundary layer below the mid-Pacific velocity anomaly identified from a PcP precursor［J］.J Geophys Res，100（B10）：20359－－20365.

Rial J A.1978.On the focusing of seismic body waves at the epicentre’s antipode［J］.Geophys J Int，55（3）：737－－743.Rial J A，Cormier V F.1980.Seismic waves at the epicenter’s antipode［J］.J Geophys Res，85（B5）：2661－－2668.

Soldati G，Boschi L，Piersanti A.2006.Global seismic tomography and modern parallel computers［J］.Ann Geophys，49（4／5）：977－－986.

Sun D Y，Helmberger D V，Jackson J M，Clayton R W，Bower D J.2013.Rolling hills on the core-mantle boundary［J］.Earth Planet Sci Lett，361：333－－342.

Sze E K M，van der Hilst R D.2003.Core mantle boundary topography from short period PcP，PKP，and PKKP data［J］.Physics Earth Planet Inter，135（1）：27－－46.

Wysession M E，Okal E A，Bina C R.1992.The structure of the core-mantle boundary from diffracted waves［J］.J Geophys Res，97（B6）：8749－－8764.