由于觀測成本低、深部信号信噪比高，接收函數目前是殼幔結構探測最常用的方法。但其分辨率遠遠低于深反射，通常僅能夠提供地殼變形的基本特征，無法在殼内變形、發震構造、控礦構造等調查研究中發揮關鍵作用。雖然近年來，短周期密集台陣通過加密觀測提高了接收函數成像的分辨率，但由于地震波頻率的限制，這種提高是有限的。分辨率不僅與台間距有關，也與地震波的頻率或波長有關。由于接收函數使用的遠震體波信号傳播了3000至1萬公裡，是以其高頻成分衰減殆盡，低頻限制了接收函數空間分辨率的提高。
為了提高天然地震方法探測地殼速度結構的分辨率，中國科學院地質與地球實體研究所岩石圈演化國家重點實驗室碩士研究所學生張良雨和導師田小波研究員提出了Pn波接收函數方法。地震學中常見的折射波是由殼内地震産生的P波以臨界角由地殼向下進入地幔頂部并沿地幔頂部滑行，滑行的同時部分能量傳回地殼、向上傳播并被台站記錄，也被成為Pn震相（如圖1所示）。Pn主要出現在2°-15°震中距範圍内，傳播距離遠小，衰減少，如圖2所示，相對于遠震體波（低于3Hz），Pn波富含高頻成分（20-30Hz）。
圖1 Pn波傳播路徑示意圖。不同顔色虛線表示不同震中距Pn波的傳播路徑
圖2 遠震P波和Pn波的波形頻率範圍比較。藍色和紅色垂直虛線分别是長度為15秒和20秒的波形信号截取時窗。寬頻帶台站遠震P波是西安台震中距為41.6°（約4500公裡）的5.6級地震的P波記錄。寬頻帶台站Pn波是該台震中距6.2°（約700公裡）的5.6級地震的Pn波記錄。由于寬頻帶台陣資料采樣率限制，可分析的頻率最高位10赫茲。短周期台站Pn波是震中距11.1°（約1200公裡）的4.8級地震的Pn波記錄
傳統的接收函數通過遠震波形反褶積獲得，需要滿足平面波、無其它震相幹擾和近垂直入射3個假設。而該項研究注意到Pn波有如下特點：（1）Pn折射波傳回地殼時均以臨界角進入地殼，在很長滑行傳播距離内，傳回地殼的折射波傳播角度保持不變，也就是說折射波由上地幔傳回地殼傳播時是滿足平面波要求的；（2）Pn波形在震中距5°-15°範圍内可以避免其它震相的幹擾；（3）通過波形記錄的質點運動極化分析可以将Pn波和它的轉換S波分離，不需要滿足近垂直入射的要求。因而可以使用震中距5°-15°範圍内的Pn波波形記錄，通過合适的時窗長度截取波形，通過極化分析和坐标旋轉分離Pn波分量和轉換S波分量，再利用這兩個分量的反褶積，可獲得Pn波接收函數。如圖3所示，通過雙層地殼模型的波場數值模拟合成Pn波記錄，通過反褶積運算獲得的Pn波接收函數包含了莫霍面和康拉德界面的轉換波。
通過理論合成模型分析（圖3）和實際寬頻帶台站的Pn波接收函數與遠震接收函數的對比（圖4），表明Pn波接收函數能夠可靠的反應地殼速度結構。同時還可以看到，由于Pn波向上傳播的角度遠大于遠震P波，導緻Pn波接收函數中轉換波的振幅明顯大于在遠震接收函數中的轉換波振幅。這一特點，有助于通過Pn波接收函數對殼内速度差小的界面進行成像或速度反演。将新方法應用于跨陰山的南北向短周期密集台陣，獲得Pn波接收函數時間剖面（如圖5所示）。通過不同頻率範圍的比較可以看出，由于Pn波富含高頻成分，可以極大地提升接收函數的波形頻率。
圖3 理論模型合成Pn波及其Pn波接收函數。(a)為簡單的雙層地殼模型，包含莫霍界面(Moho)和康拉德界面(Conrad)；(b)為利用波場數值模拟合成的震中距12°的Pn波3分量波形，其中300秒以後的振幅被縮小了20倍；(c)為根據偏振方向進行坐标系旋轉投影的Pn波分量和轉換S波分量，這兩個分量進行反褶積獲得Pn波接收函數；(d)其中Pcs和Pms分别為來自康拉德界面和莫霍界面的轉換波
圖4 Pn波接收函數與遠震接收函數對比。(a)和(b)分别是西安台和昆明台。Pn波接收函數(PnRF)與遠震接收函數(tel-RF)的時間軸均已動校正的統一的射線參數。圖中垂直的紅色虛線和藍色虛線是基于前人地殼速度模型的預測莫霍面轉換波到時
圖5 短周期密集台陣的Pn波接收函數時間剖面。(a)、(b)和(c)分别使用了高斯系數為5.0、10.0和15.0的低通濾波。紅色虛線訓示可連續追蹤的轉換波震相，其中6秒左右為莫霍面的轉換波震相
由于頻率高，該方法可以将天然地震地殼結構探測的分辨能力提升到接近人工源深地震反射的水準（如圖6所示）。目前的天然地震轉換波方法（遠震接收函數）頻率低（<3.0Hz），而震中距為5°-15°的Pn波頻率可以高達20Hz以上，該頻率的Pn波接收函數在40公裡以淺的橫向和垂向分辨率分别為1-3公裡和0.2公裡，非常接近深地震反射的探測水準。是以，該方法有望成開展造山帶精細結構、大型成礦帶和大型地震帶地殼精細結構探測的 “利器”。
圖6 不同探測方法的分辨率。黑色實線、虛線和點線為遠震接收函數(tele-RF)不同頻率波形的分辨率随深度的變化；紅色實線、虛線和點線為Pn波接收函數(PnRF)不同頻率波形的分辨率随深度的變化；藍色實線為深反射(DSR)的分辨率随深度的變化
研究成果發表在國際地學權威學術期刊JGR-Solid Earth（張良雨，田小波*. Pn-wave receiver function [J]. Journal of Geophysical Research: Solid Earth, 2024, 129: e2023JB028318. DOI: 10.1029/2023JB028318.）。研究受到國家自然科學基金（42030308,41974053）的資助。
美編：傅士旭(華東師大）
校對：萬鵬