鄂爾多斯地塊東南緣地帶Moho深度變化特征研究

任 梟，徐志國，楊 輝，陳宏峰，鄒立曄

1 中國地震局地球物理研究所，北京 100081

2 中國地震臺網中心，北京 100045

1 引 言

自華北克拉通形成以來，經過了多期不同性質的構造活動，內部形成了構造特征不同的次級構造單元［1］.我們的研究區域位于鄂爾多斯東南緣，這個區域即是華北克拉通北東—北北東向構造與近東—西向構造的交接部位，又是新生代中國大陸西部碰撞擠壓到東部伸展拉張動力學轉換的關鍵地區［2］.鄂爾多斯地塊東部邊緣的核心是一個隆起區，以呂粱、太行山脈為特征；南緣是渭河盆地，是以下降接收沉積物為主的斷陷盆地.研究區區域的地質構造如圖1所示.

鄂爾多斯地塊東南緣邊界區域構造活動強烈［3］，地震活動頻繁，具有較強的地震活動性，多年來一直是地學家研究的熱點地區之一.歷史上在周邊斷裂系內共記錄到6級以上的淺源強烈地震22次，其中8級及以上地震2次，即東部邊緣的1303年9月25日洪洞8級地震和南部邊緣的1556年2月2日華縣8級地震；震級7≤M＜8級地震5次；其余則為6≤M＜7級的地震［4］.

為了解鄂爾多斯地塊東南緣區域的地殼深部結構特點，本文收集該區域內分布的寬頻帶數字地震固定臺站記錄的大量遠震體波波形資料，應用頻率域反褶積方法提取遠震P波接收函數，由H－κ方法測定了各臺站下方的Moho深度和Vp／Vs值，以研究鄂爾多斯東南緣地區的Moho深度及其變化特征.以期為研究鄂爾多斯地塊東南緣與周邊塊體之間的深部結構、地殼巖石動力學過程以及強震發生機制等提供可靠的地球物理學依據.

2 研究方法

接收函數是研究殼幔結構的主要手段之一，國內外學者均對其做了廣泛的研究與應用，并取得了一定的研究成果.近些年來，隨著大規模數字地震臺站的建設，臺站觀測密度顯著增強，積累了高質量的地震數據，使得接收函數研究向著高分辨率的方向發展，為獲取可靠的殼幔結構信息提供基礎.接收函數表示的是臺站下方入射P波的結構響應，當遠震P波通過地球內部，在速度不連續面P波轉化為S波，我們可以通過轉換波和多次波的走時提取，對Moho深度和Vp／Vs提供有效約束［5-8］.

本研究采用改進的頻率域反褶積技術，遠震P波波形的頻譜R（ω）除以源頻譜S（ω）：

S＊（ω）是S（ω）頻譜的共軛，e－ω2／4α2是高斯低通濾波器，用以去除高頻噪音，cσ20也叫水準量，用以壓制S（ω）頻譜的歧值點，可以發揮穩定反褶積的作用.由于在殼幔邊界存在著大的速度間斷面，P波在Moho面的轉化震相Ps經常是直達P波之后最主要的震相.理想情況下，初至P波的轉化波Ps和它的地殼多次反射震相PpPs和PpPs＋PsPs均清晰可見，并且具有相近的振幅.在給定地殼平均速度的條件下，P波接收函數中Ps震相的到時實際上指的是Ps震相和直達P波的到時差，它取決于間斷面的深度，射線參數以及臺站下方P波、S波速度結構.因此，利用P波接收函數獲得臺站下方地殼內部的S波速度結構是一個非線性反演的問題.公式（2）所給出的地殼厚度H和地殼平均速度之間是多解的關系.

使用多次波震相可以有效地壓制其多解性，通過公式（3）和公式（4）可以獲得地殼厚度H 和Vp／Vs的精確估計［9-10］.接收函數的疊加可以有效地壓制隨機噪聲，增強信號的信噪比，這就需要在時間域內對一簇事件進行疊加［11］.由于從實際的接收函數中拾取Pms震相和PpPms震相是非常困難的，Zhu和 Kanamori［12］提出用 H－Vp／Vs域網格搜索和疊加的方法來確定Moho界面深度和Vp／Vs，并且這一方法得到了廣泛的應用.使用Ps，PpPs以及PpSs＋PsPs三個震相的振幅構造目標函數：

圖1 鄂爾多斯東南緣地質構造圖（引自“鄧起東，基于GIS的地震分析預報系統”）Fig.1 Geophysics of Ordos Southeastern areas

其中r（ti）表示接收函數的振幅，t1，t2，t3是分別相應于依據地殼厚度H和Vp／Vs計算的Ps，PpPs和PpSs＋PsPs的到時，ω是每個震相的權值.Ps震相的振幅最大，所以賦予的權值也最大（ω1＞ω2＋ω3）.如果多次波清楚的話，那么最大值所對應的地殼厚度H和Vp／Vs可以通過在給定的取值范圍內搜索得到.當三個震相（Ps，PpPs和PpSs＋PsPs）振幅值的疊加對應正確的H 和κ時，三個震相相交，S（H，κ）趨于最大值.

通常情況下，為了增強信噪比和消除不均勻性造成的散射波、背景噪音以及計算誤差的影響，每個臺站的不同震中距的接收函數進行疊加得到平均的接收函數.

3 資料分析與處理

本文利用山西省地震臺網、陜西省地震臺網的固定測震臺站2009—2010年間記錄到的遠震體波資料.臺站分布如圖1所示，這些臺站所配置的地震儀類型多達數十種，均為三分向寬頻帶儀器，頻帶寬度分別為20s～40Hz、60s～40Hz或120s～40Hz，采樣率均為100sps，臺站的分布覆蓋了鄂爾多斯地塊東南緣邊界區域.我們選取了858個震級Mb≥5.4，震中距范圍30°～90°的遠震觀測數據，事件分布情況見圖2，計算得到對應單臺下方的Moho面深度和Vp／Vs值.

圖2 選用的遠震事件分布圖Fig.2 Locations of events used in the study

在數據預處理工作中，要考慮具有較完整的方位角分布和清晰記錄，以利于獲得較好的平均結果；參照目錄篩選每個臺站的原始記錄，力求震相清楚、初動尖銳、波形簡單；接著對選取的數據作坐標旋轉、去除均值和傾斜分量、加時間窗等預處理；之后采用頻率域反褶積計算方法提取單臺單個事件的接收函數，在此基礎上，選取波形數據質量較好、信噪比較高、多次反射震相清晰的接收函數用于計算Moho深度和波速比.

我們采用Zhu［12］提出的 H－Vp／Vs域網格搜索和疊加的方法來確定Moho界面深度和Vp／Vs值，其優點是：無需人工標定Pms、PpPms和PpSms＋PsPms震相的到時.接收函數掃描法計算Moho深度和波速比過程中，給定的地殼P波平均速度為6.3km／s，對于與 Moho界面相關的 Pms、PpPms和PpSms＋PsPms震相采用0.7，0.2，0.1的加權值，H 和Vp／Vs的變化范圍為20～60km和1.6～1.9搜索確定最優解.

4 結果與討論

通過研究區內各臺站接收函數得到各臺站下方的Moho深度和波速比.文中以臺站ANZ為例說明獲取其下方Moho深度和波速比.圖3表示臺站ANZ記錄遠震資料提取的徑向接收函數.橫軸代表P波初至后的時間，縱軸代表反方位角.由圖可見，臺站ANZ的接收函數在直達P波之后5.3s處可以非常清楚地識別出Moho界面的Ps震相.臺站ANZ接收函數徑向分量極性和振幅隨著方位角的變化而變化，可能是由于地殼各向異性或殼內傾斜界面所致.此外，在圖中P波初至震相至Moho界面Ps轉換震相范圍間，存在若干個地殼內部次一級速度界面產生的Ps轉換波，這可能由于臺站ANZ位于臨汾盆地邊緣，臺站下方地質結構復雜，地殼演化過程后期由于強烈的地質作用影響，盆地大量的沉積物沉積成層，地殼內的多層界面產生Ps轉換波.

采用網格搜索和疊加方法得到臺站ANZ下方的 Moho深度為43.70km，Vp／Vs為1.84，如圖4所示，圖中漸變的顏色表示的是s（H，κ）值的變化，當顏色最深時，該數值達到最大，最大值對應的Moho深度值及波速比即為疊加結果.

深反射人工測深（Deep Seismic Sounding）是依靠多次記錄疊加獲取Pm＊震相去獲取Moho深度.因為DSS是人工源，所以源多，接收點密，適于獲取精細結構，但是由于造價高，即便能較好地測得每條測線上的結構，測線數量也是有限的.本研究使用固定臺站記錄的天然地震事件，通過反演接收函數去找Moho面.這種方法優點是造價低，但只能大尺度地對Moho形態進行研究，不太容易獲取其精細結構.近年來，隨著固定臺站的加密建設，其精細程度也會隨之增加.下面將利用接收函數方法獲取的Moho深度與DSS方法獲取的結果進行對比.按照臺站ANZ的經緯度，在Li［13］于2006年利用DSS方法得到的“中國大陸Moho深度等值線分布圖”中找到其對應點的Moho深度為44.00km.我們隨機抽取了若干臺站，按照同樣方法對應比較了Moho深度值.可以看出：在誤差范圍內，本研究結果與Li［13］研究結果基本一致.

表1 Moho深度值對比表Table 1 Comparison of Moho depth

利用鄂爾多斯地塊東南緣區域范圍內的58個臺站的遠震波形資料，采用接收函數地殼厚度和Vp／Vs掃描方法，得到該區域的Moho深度和波速比，詳見表2.

鄂爾多斯地塊東南緣地區Moho界面深度界于31.0～53.1km之間，平均的 Moho深度為40.7km，Moho最厚為陜西安康臺的53.10km±1.50km，最薄為陜西藍田臺的31.00km±1.76km，二者相差22.1km，研究區域整體存在明顯的起伏變化.其東緣的Moho深度介于33.4～45km之間，以太原斷陷盆地附近的Moho深度較淺，最淺處為33.4km；太原斷陷盆地東部北段的延懷盆地、蔚縣盆地、陽原盆地和南段的臨汾盆地附近Moho深度變化不大，平均深度為40km.而在東部邊緣的東側，因存在著北北東走向的山西斷陷帶，導致塊體邊緣的Moho深度要大于塊體內部的Moho深度.鄂爾多斯地塊南緣的Moho深度介于31～53.1km之間，自東段向西段Moho深度逐漸變大，從渭河盆地附近的31km增厚至秦嶺造山帶地段的53.1km.

圖3 臺站ANZ的徑向接收函數Fig.3 Radial receiver function of ANZ

上述數值表明，鄂爾多斯地塊東南緣地殼結構的橫向變化劇烈，殼內結構復雜.將結果匯總使用差值法繪制了鄂爾多斯地塊東南緣及其周邊連續的Moho深度等值線分布圖，詳見圖5.由圖可見，鄂爾多斯地塊東南緣的Moho深度等值線分布趨勢為NNE方向，Moho深度與構造環境有關：在凹陷和盆地區，中心處的Moho深度小于邊緣處，例如：東緣地區的太原盆地、南緣地區的渭河盆地；相反在地形抬升地區，例如：南緣的秦嶺造山帶，中心區域的Moho深度大于邊緣地帶.

圖4 臺站ANZ的徑向接收函數時間剖面（a）HK搜索，橫軸是波速比，H是深度；（b）根據計算的深度得到的的理論P波和轉換震相的位置，橫軸為時間，單位s；（c）所有疊加的接收函數和理論的震相位置，橫軸是時間，單位s，縱軸是歸一化的幅度.Fig.4 Radial receiver function profile of ANZ（a）The example of H－κstacking results obtained from Station ANZ.Horizontal axis means Vp／Vsratio；vertical axis means depth；（b） The position of theoretical converted wave.Horizontal axis means time；vertical axis means ray parameter；（c）The receiver functions for Station ANZ.Horizontal axis means time；vertical axis means normalized amplitude.

圖5 研究區域總的Moho深度等值線分布圖Fig.5 Moho depth isoline map in this area

表2 Moho深度與Vp／Vs值Table 2 Moho depth and Vp／Vsratio

關于鄂爾多斯地塊Moho深度的資料，由于資料來源和解釋方法上的不同，具體數據存在差異.如Li［13］于2006年研究所得“中國大陸 Moho深度等值線分布圖”中鄂爾多斯塊體對應的Moho深度為42～44km；陳九輝等［14］于2005年在研究鄂爾多斯地塊上地幔S波速度結構時得出鄂爾多斯地塊的地殼厚度為40km左右；而韓忠東、楊峰杰等［15-16］得出鄂爾多斯地塊南緣環形區域內部地殼厚度為40～42km，鄂爾多斯地塊北部地殼厚度為47～49km.經對比分析，本研究結果具有較高可信度.

由圖6可以看出鄂爾多斯東南緣區域整體的Vp／Vs分布具有明顯的分區特征.整個研究區域地殼波速比介于1.60～1.90之間，平均波速比為1.77.地殼的泊松比σ可以直接由Vp／Vs求取，研究區域的泊松比介于0.18～0.31之間.陜西富平臺站、藍田臺站和山西的臨汾臺站、鎮川臺站的波速比最高，均為1.90，陜西商南臺站、山西的右玉臺站和大同臺站波速比相對低，分別為1.60、1.63、1.68，對應泊松比為0.18、0.20、0.23.相對高的波速比多存在于沉積盆地斷陷中心地帶及附近凹陷區域；而較低的波速比則分布在鄂爾多斯地塊東部邊緣太行山脈與南部邊緣秦嶺造山帶等隆起區域附近.不同臺站的波速比反映了研究區域的地殼分層結構明顯，具有很強的橫向不均勻性.根據已有研究成果表明，地殼泊松比與介質的巖石組成有關，常見巖石的泊松比值為0.25左右，而石英含量的提高會顯著地降低巖石的泊松比值，例如砂巖、粉砂巖及石英巖的泊松比值都很低；而灰巖、白云巖等碳酸鹽巖以及片麻巖等變質巖的泊松比值則相對較高［17］，如輝長巖的為0.30.因此，研究區的泊松比值偏低，這些特征對應了鄂爾多斯塊體自新生代以來即以垂直升降為主的演化歷史［18］，也符合關于鄂爾多斯盆地塊體在上地殼較淺部位主要以砂巖或變質鎂鐵巖為主的實際情況［19］.

圖6 研究區域總的Vp／Vs等值線分布圖Fig.6 Vp／Vsratio isoline map in this area

從Vp／Vs值分布體現了：相對高的波速比多存在于沉積盆地斷陷中心地帶及附近凹陷區域，例如臨汾盆地與渭河盆地等裂谷斷陷盆地附近存在較高的Vp／Vs值，其等值線分布圖較好地勾勒出塊體東南緣的主要沉積盆地.較低的波速比則分布在鄂爾多斯地塊東部邊緣太行山脈與南部邊緣秦嶺造山帶等隆起區域附近.塊體東緣邊界帶的Vp／Vs值分布具有顯著的南、北分段特征，顯現出凹陷與突起相間分布的特征.總之，不同臺站的波速比反映了研究區域的地殼分層結構明顯，具有很強的橫向不均勻性，計算結果很好地印證了鄂爾多斯地塊地殼結構的復雜多變性.

5 結 論

本文使用遠震接收函數得到了鄂爾多斯地塊東南緣地帶的Moho深度和Vp／Vs值分布，分析表明：

（1）應用固定寬頻帶遠震接收函數是一種有效的研究臺站下方地殼上地幔結構的方法，其結果與DSS結果相對應，具有一定的可信度.隨著固定臺站波形記錄的增加，可以進一步提高反演結果的可靠性.

（2）鄂爾多斯地塊東南緣地帶的Moho深度分布具有明顯的分塊特征.塊體內部結構比較穩定，盆地中部地殼厚度變化較為平緩，塊體外側地殼結構與塊體東緣內側地殼深度具有明顯的差異性，自東向西地殼厚度呈現逐漸增厚的趨勢，很好地對應了其地質構造特點.

（3）研究區域Vp／Vs值分布也具有明顯的分塊特征，不同臺站的波速比反映了研究區域的地殼分層結構明顯，具有很強的橫向不均勻性.相對高的波速比多存在于沉積盆地斷陷中心地帶及附近凹陷區域，較低的波速比則分布在鄂爾多斯地塊東部邊緣太行山脈與南部邊緣秦嶺造山帶等隆起區域附近.

（4）原本計算了59個臺站的數據，但在本文中未使用山西昔陽臺的結果.原因是：從該臺的接收函數圖看出：來自北西-南東方向比較聚焦，震相穩定；來自北東-南西方向聚焦較差.這可能是臺站的基址問題造成的，該臺站處在山西地塹的東南緣邊界，可能是在具有一定寬度的一級或者二級邊界斷裂帶內，例如：太行山斷裂帶.所以對應的莫霍轉換震相不清晰，得不到確定的Moho深度值.這個結果正說明：地表和深部地質構造條件決定了數據質量.該地帶作為一級構造特征和重要的布格重力異常帶，這種現象值得深入研究，我們下一步工作將分別選用位于平原、過渡帶和斷陷盆地內的臺站資料進行計算對比，結合分析地表地質構造現象來認識此類問題.

（References）

［1］嘉世旭，張先康.華北不同構造塊體地殼結構及其對比研究.地球物理學報，2005，48（3）：611－620.Jia S X，Zhang X K.Crustal structure and comparison of different tectonic blocks in North China.Chinese J.Geophys.（in Chinese），2005，48（3）：611－620.

［2］陸一鋒，徐鳴潔，王良書等.相位加權疊加方法在探測鄂爾多斯東南緣地殼結構中的應用.高校地質學報，2011，17（4）：562－568.Lu Y F，Xu M J，Wang L S，et al.An application of phaseweighted stacks to study the crustal structures in the southeastern margin of the Ordos block.Geological Journal of China Universities （in Chinese），2011，17（4）：562－568.

［3］張學民，刁桂苓，束沛鎰.鄂爾多斯塊體及其東南緣剪切波速度結構與波速比研究.中國地震，2004，20（1）：53－63.Zhang X M，Diao G L，Shu P Y.Study on shear wave velocity structure and velocity ratio under the Ordos Plate and its southeastern marginal areas.Earthquake Research in China （in Chinese），2004，20（1）：53－63.

［4］任雪梅，高孟潭，馮靜.鄂爾多斯地塊周緣歷史強震的影響烈度研究.中國地震，2010，26（4）：450－456.Ren X M，Gao M T，Feng J.Study on historical strong earthquakes affected intensity of counties and cities in the Erdos area.Earthquake Research in China （in Chinese），2010，26（4）：450－456.

［5］Ducker K G，Sheehan A F.Mantle discontinuity structure from midpoint stacks of converted P to S waves across the Yellow stone hotspot track.Journal Geophysical Research，1997，102（B4）：8313－8327.

［6］Li A B，Fischer K M，Wysession M E，et al.Mantle discontinuities and temperature under the North American continental keel.Nature，1998，395（6698）：160－163.

［7］Jones C H，Phinney R A.Seismic structure of the lithosphere from teleseismic converted arrivals observed at small arrays in the southern Sierra Nevada and vicinity，California.Journal of Geophysical Research，1998，103（B5）：10065－10090.

［8］Kosarev G，Kind R，Sobolev S V，et al.Seismic evidence for a detached Indian lithospheric mantle beneath Tibet.Science，1999，283（5406）：1306－1309.

［9］Zandt G，Myers S C，Wallace T C.Crust and mantle structure across the Basin and Range－Colorado Plateau boundary at 37°N latitude and implications for Cenozoic extensional mechanism.Journal of Geophysical Research，1995，100（B6）：10529－10548.

［10］Zandt G， Ammon C J. Continental crust composition constrained by measurements of crusta1Poissson′s ratio.Nature，2002，374（6518）：152－154.

［11］Owans T J，Zandt G，Taylor S R.Seismic evidence for an ancient rift beneath the Cumberland plateau，Tennessee：A detailed analysis of broadband teleseismic P waveforms.Journal of Geophysical Research，1984，89（B9）：7783－7795.

［12］Zhu L P，Kanamori H.Moho depth variation in southern California from teleseismic receiver functions.Journal of Geophysical Research，2000，105（B2）：2969－2980.

［13］Li S L，Mooney W D，Fan J C.Crustal structure of mainland China from deep seismic sounding data.Tectonophysics，2006，420（1－2）：239－252.

［14］陳九輝，劉啟元，李順成等.青藏高原東北緣——鄂爾多斯地塊地殼上地幔S波速度結構.地球物理學報，2005，48（2）：333－342.Chen J H，Liu Q Y，Li S C，et al.Crust and upper mantle S－wave velocity structure across Northeastern Tibetan P1ateau and Ordos block.Chinese J.Geophys.（in Chinese），2005，48（2）：333－342.

［15］韓忠東，陳顒，楊鋒杰等.鄂爾多斯地臺南部邊緣環形構造的GPS證據及深部地球動力學解釋.中國西部科技，2008，7（3）：1－5.Han Z D，Chen Y，Yang F J，et al.The GPS evidences and depth earth dynamics explanations of the ring structure in south Ordos platform.Science and Technology of West China （in Chinese），2008，7（3）：1－5.

［16］楊鋒杰，王明鎮，李增學等.鄂爾多斯盆地南部環形影像特征及地質意義.地球科學與環境學報，2006，28（3）：37－41.Yang F J，Wang M Z，Li Z X，et al.Ring structure in south Ordos Basin and its geologic genesis.Journal of Earth Sciences and Environment，2006，28（3）：37－41.

［17］李永華，吳慶舉，安張輝等.青藏高原東北緣地殼S波速度結構與泊松比及其意義，地球物理學報，2006，49（5）：1359－1368.Li Y H，Wu Q J，An Z H，et al.The Poisson ratio and crustal structure across the NE Tibetan Plateau determined from receiver functions.Chinese J.Geophys.（in Chinese），2006，49（5）：1359－1368.

［18］廖昌珍，張岳橋，溫長順.鄂爾多斯盆地東緣邊界帶構造樣式及其區域構造意義，地質學報，2007：81（4）：466－474.Liao C Z，Zhang Y Q，Wen C S.Structural styles of the Eastern Boundary zone of the Ordos Basin and its regional tectonic significance.Acta Geologica Sinica （in Chinese），2007，81（4）：466－474.

［19］江為為，郝天珧，宋海斌.鄂爾多斯盆地地質地球物理場特征與地殼結構.地球物理學進展，2000，15（3）：45－53.Jiang W W，Hao T Y，Song H B.Crustal structure and geological and geophysical features of Ordos basin.Progress in Geophysics（in Chinese），2000，15（3）：45－53.