殷常陽 石永祥 伍晗 寧杰遠(yuǎn),2,?

基于合成數(shù)據(jù)對高鐵震源虛擬道集生成方法的驗證

殷常陽1石永祥1伍晗1寧杰遠(yuǎn)1,2,?

1.北京大學(xué)地球與空間科學(xué)學(xué)院, 北京 100871; 2.河北紅山地球物理國家野外科學(xué)觀測研究站, 邢臺 054000; ?通信作者, E-mail: njy@pku.edu.cn

根據(jù)穩(wěn)相點分析的疊加方法原理, 提出一種波形疊加方法, 消除臺站記錄中高鐵波場的串?dāng)_噪聲, 生成虛擬道集, 并用合成的高鐵數(shù)據(jù)進行驗證。通過建立高鐵橋梁振動方程, 獲得精細(xì)的高鐵震源時間函數(shù), 然后使用有限元方法,在分層介質(zhì)模型中進行高速鐵路波場正演模擬計算, 并對不同速度的高鐵列車正演模擬記錄進行時移疊加, 獲得虛擬道集, 從而驗證疊加方法的正確性。

高鐵; 干涉場; 震源時間函數(shù); 虛擬道集

高鐵列車經(jīng)過高架橋時會引起橋體振動, 并通過橋墩傳導(dǎo)到地下。關(guān)于高鐵列車經(jīng)過橋體的振動模型和高鐵通過時引起周圍地表振動的響應(yīng)及頻率, 已有大量的研究[1?9]。然而, 高鐵振動信號復(fù)雜, 非常難以處理。溫景充等[10]基于新的密集觀測資料, 指出高鐵震源激發(fā)的波場是穩(wěn)定的、重復(fù)性強的地震干涉波場。高鐵相鄰橋墩的震源之間存在串?dāng)_噪聲[11], 如果采用傳統(tǒng)地震勘探中常用的互相關(guān)等混響分離方法, 無法得到正確的虛擬道集。最近, 石永祥等[12]通過分析高鐵波場的特征, 得到其相干振幅和相干相位的解析表達(dá)式, 發(fā)現(xiàn)高鐵波場有干涉波場的特性, 存在主瓣, 可以利用這個特征進行地下介質(zhì)成像。通過穩(wěn)相點分析, 不同速度的列車記錄疊加可以消除源與源之間的串?dāng)_噪聲, 也可以生成高鐵的虛擬道集[13]。

為了驗證生成虛擬道集的原理, 需要正演出高速鐵路通過高架橋的波場數(shù)據(jù)。本文參照 Wu 等[5]模擬列車通過橋梁時引發(fā)周圍地表位移的振動模型, 對高架橋采用歐拉?伯努利橋梁模型, 對樁體和土體采用開爾文體模型, 以便得到更貼近真實情況的高鐵震源時間函數(shù)。利用得到的震源時間函數(shù), 通過譜元法正演彈性條件下的高速鐵路波場, 得到波場合成數(shù)據(jù)?；趯铣蓴?shù)據(jù)的處理, 合成虛擬道集, 并用合成的高鐵數(shù)據(jù)驗證該方法的正確性。

1 地震勘探中混合震源數(shù)據(jù)分離

其中,

將式(6)代入式(7):

式(8)的兩邊同時出現(xiàn) [?]-1, 是一個調(diào)節(jié)所有頻率下振幅的因子[16]。將式(8)的右邊?一項展開:

如果震源編碼矩陣中的震源激發(fā)時間1,2, …,t是隨機的, 即滿足下式[17]:

其中, ???表示集平均(ensemble average),表示單位 矩陣。

若震源激發(fā)時間1,2, …,t是隨機的, 那么式(9)是對角線為1, 非對角線為0的單位矩陣, 即滿足式(10), 那么式(10)就能保證串?dāng)_噪聲的壓制。因此, 震源激發(fā)時間的隨機性保證了地震勘探中多震源數(shù)據(jù)的分離, 能夠壓制串?dāng)_噪聲, 這是地震勘探中混合震源數(shù)據(jù)分離的關(guān)鍵[14?18]。

實際上, 高速鐵路波場的震源激發(fā)時間延遲不是隨機的, 而是依次等間距線性延遲激發(fā), 不滿足式(10)。所以, 串?dāng)_噪聲不能通過傳統(tǒng)地震勘探方法壓制, 而是要基于高速鐵路波場的特點提取虛擬格林函數(shù), 并合成虛擬道集。

2 基于穩(wěn)相點疊加的虛擬道集生成方法

在頻率域中, 所有橋墩產(chǎn)生的震源時間函數(shù)all()激發(fā)的波場被臺站記錄到的垂向數(shù)據(jù)信號(g,)為

石永祥等[12]把方程(12)右端進行如下變換:

殷常陽等[13]基于穩(wěn)相點性質(zhì), 提出對不同的速度進行疊加:

圖1 相干相位在不同頻率、介質(zhì)速度和列車速度下的分布

3 高速鐵路震源時間函數(shù)的正演

為了驗證上述方法, 需要模擬行駛在高架橋上高鐵的地震波場。首先, 需要得到一個貼近真實的震源時間函數(shù)。

為了方便推導(dǎo), 取橋梁的微元段 d, 使用如圖2所示的歐拉?伯努利橋梁模型。其中,(,)是單位長度上外力,(,)為垂直撓度,為橋梁轉(zhuǎn)動慣量,為橋梁彈性模量,為橋梁單位長度的質(zhì)量。歐拉?伯努利橋梁方程[5]為

圖2 歐拉?伯努利橋梁模型[5]

高鐵由多節(jié)車廂組成, 每節(jié)車廂都可以視為負(fù)載, 且每節(jié)車廂都有前后兩組車輪, 由于列車的速度恒定, 所以每節(jié)車廂到達(dá)橋墩的時間間隔可以視為恒定[5?6], 則外力(,)表示為

橋梁振動方程的初始條件和邊界條件[5?7]為

式(20)中,表示橋墩支撐座的等效彈性系數(shù),表示橋梁的轉(zhuǎn)動慣量,表示橋梁的楊氏模量。

求解式(17)～(20), 可以得到橋梁在外力作用下的垂直撓度(,)。

前后輪一起對橋墩彈性支座施加的力為

求得total()后, 還需要求解高鐵移動對橋墩的影響, 從而計算橋墩與土體之間的相互作用, 最后求出激發(fā)波場的震源時間函數(shù)s()。

根據(jù)文獻[5?8]的研究結(jié)果, 列車在高架橋引起的地面振動模型如圖3所示。彈性系數(shù)為的彈簧支座下方, 連接著質(zhì)量為的橋梁支柱(可以視為剛體)。為了保證支柱穩(wěn)定, 底部有樁帽托住。樁帽有5個樁體支撐, 每個樁體可以視為開爾文體模型。同時, 土壤地基也可視為開爾文體模型。所以在頻率域中, 圖3所示模型的定解方程[5?8]為

由式(24)可以得到樁體的位移:

圖3 高鐵?橋梁?支撐住?樁體?土體模型[5?8]

其中,

得到頻率域的震源時間函數(shù)s():

把式(28)從頻率域轉(zhuǎn)換到時間域:

s()即為時間域的震源時間函數(shù), 也是由高鐵移動引起的震源作用于橋墩時, 激發(fā)波場的震源時間函數(shù)。不同車速時橋梁長度=32m 的震源時間函數(shù)如圖4所示, 其中高鐵、橋梁、樁體和土體的參數(shù)來自文獻[5]。

all(,,,)的加載方向豎直向下。

利用譜元法, 由震源時間函數(shù)all(,,,), 可以進行高速鐵路波場正演: 頂層是自由表面, 其余是 Clayton-Engquist 吸收邊界[19]。本文采用圖5所示的介質(zhì)速度模型: 第一層密度1=2.2kg/m3, P波速度1=1600m/s, S 波速度1=400m/s, 厚度1=150m; 第二層密度2=2.4kg/m3, P 波速度2=2400m/s, S波速度2=750m/s, 厚度2=250m。

圖6是高速鐵路平面分布圖, 行駛在鐵路高架橋(黑色實線)上高鐵的速度為, 橋墩間距=32m, 共有63個橋墩。

4 合成數(shù)據(jù)結(jié)果

根據(jù)式(15)和(16), 合成虛擬道集的步驟如下。

1)消除震源影響。高速列車速度變化造成震源時間函數(shù)變化, 所以需要反卷積震源時間函數(shù), 要注意保留震源時間函數(shù)時間延遲的信息。

圖4 橋梁長度L=32 m時不同車速下的震源時間函數(shù)

圖5 介質(zhì)的速度模型示意圖

實線表示高速鐵路, 由高架橋支撐; 五角星表示其中兩個橋墩, 三角形表示布設(shè)在地表的臺站; 虛線表示 x 軸從 100 m 至1900 m, y軸30 m處, 臺間距為100 m的測線

4)疊加不同速度的波形。上述每個臺站獲得16 個不同高速列車速度產(chǎn)生的波形, 每個速度對應(yīng)的波形起始時間為1,2, …,16。首先需要對齊波形, 即把不同速度波形的起始時間1,2, …,16平移到零時刻, 然后進行疊加, 即獲得第號橋墩到臺站的格林函數(shù)。

選取第20號橋墩(608m, 500m)和第30號(928m, 500m)橋墩, 求取這兩個橋墩到臺陣的格林函數(shù)。圖7顯示兩個橋墩分別到坐標(biāo)為(1000m, 30m)的臺站(圖6中三角形所示)分量虛擬道集中的單道, 反射S波的理論到時滿足

圖8是第20號橋墩和第30號橋墩分別到軸從100m至1900m,軸30m處, 臺間距為100m的測線臺站(圖6中虛線所示)的虛擬道集, 反射S波的理論到時仍然滿足式(31)。

從圖7和8可以看出, 在反射S波震相到時合成數(shù)據(jù)的分量中, 彈性波存在多個震相的相互干擾。與聲波方程的模擬計算相比, 彈性波方程的計算更加復(fù)雜[13]。同時, 反射波能量強, 首波和直達(dá)波能量弱, 疊加后仍然存在部分串?dāng)_噪音, 相對來說, 反射波的到時能更好地體現(xiàn)出來。

(a)第 20 號橋墩; (b)第 30 號橋墩; 黑色虛線示意反射S波的理論到時

(a)第 20 號橋墩; (b)第 30 號橋墩; 黑色虛線示意反射 S 波的理論到時

本文只考慮震源時間函數(shù)向下時的合成記錄情況。真實的震源時間函數(shù)不僅有垂直向下的, 也有平行于鐵路方向的。在進行虛擬道集合成時, 疊加速度范圍的選取至關(guān)重要, 必須保證疊加速度項覆蓋相干相位的極大值點。本文采取的是70～100m/s的速度范圍疊加,的極大值點在此速度范圍內(nèi)是可以取到的。如果疊加速度沒有覆蓋的極大值, 則會造成虛擬道集的生成失敗, 因此需要確定的極大值點的對應(yīng)速度后再進行疊加。

高鐵列車速度在大范圍的變化可以保證有效地消除串?dāng)_噪聲。實際上, ±5m/s 的變化范圍已經(jīng)能夠滿足穩(wěn)相點求和的需求(如圖1所示, 主瓣區(qū)以外相位變化劇烈)。觀測結(jié)果顯示, 車速的實際變化幅度為±5 m/s[20], 可以滿足計算需求。

虛擬道集中, 有時橋墩附近臺站的反射S波到時不清晰(圖8)。原因是如果接收臺站過于接近橋墩, 橋墩震源振動的近場項會對波場特征產(chǎn)生較大的影響, 導(dǎo)致到時特征不顯著。因此, 在疊加不同高鐵列車的速度時, 臺站與橋墩震源應(yīng)當(dāng)保持一定的距離, 避免橋墩震源振動的近場項對虛擬道集的生成造成影響。

5 結(jié)論

本文提出一種基于穩(wěn)相點分析的波形疊加方法, 可以消除高鐵列車通過高架橋時引起的多橋墩震源的串?dāng)_噪聲, 生成特定頻率高鐵波場的虛擬道集。高架橋的橋梁采用歐拉?伯努利模型, 土體和樁體采用開爾文體模型, 通過正演模擬得到高鐵列車通過高架橋時的震源時間函數(shù), 并使用有限元方法給出分層介質(zhì)模型中高鐵波場的正演模擬結(jié)果。在此基礎(chǔ)上, 對不同速度的高鐵列車正演模擬記錄進行疊加, 獲得虛擬道集。通過對比, 發(fā)現(xiàn)結(jié)果與反射S波的理論走時相符, 驗證了本文高鐵震源虛擬道集生成方法的正確性。

[1]高廣運, 謝偉, 陳娟, 等. 高鐵運行引起的高架橋群樁基礎(chǔ)地面振動衰減分析. 巖土力學(xué), 2019, 40 (8): 3197?3206

[2]Chen Feng. Prediction and mitigation analyses of ground vibrations induced by high speed train with 3-dimensional finite element method and substructure method. Journal of Vibration and Control, 2011, 17 (11): 1703?1720

[3]Krylov V V. Generation of ground vibrations by su-perfast trains. Applied Acoustics, 1995, 44(2): 149? 164

[4]Wu W H, Smith H A. Efficient modal analysis for structures with soil-structure interaction. Earthquake Engineering & Structural Dynamics, 1995, 24(2): 283–299

[5]Wu Y S, Yang Y B. A semi-analytical approach for analyzing ground vibrations caused by trains moving over elevated bridges. Soil Dynamics and Earthquake Engineering, 2004, 24(12): 949–962

[6]Yang Y B, Hung H, Chang D W. Train-induced wave propagation in layered soils using finite/infinite element simulation. Soil Dynamics and Earthquake Engineering, 2003, 23(4): 263–278

[7]Yang Y B, Lin C L, Yau J D, et al. Mechanism of resonance and cancellation for train-induced vib-rations on bridges with elastic bearings. Journal of Sound and Vibration, 2004, 269(1/2): 345–360

[8]Yang Y B, Yau J D, Hsu L C. Vibration of simple beams due to trains moving at high speeds. Enginee-ring Structures, 1997, 19(11): 936–944

[9]陳棋福, 李麗, 李綱, 等. 列車振動的地震記錄信號特征. 地震學(xué)報, 2004, 26(6): 651–659

[10]溫景充, 寧杰遠(yuǎn), 張獻兵. 高鐵地震波場特點的理論分析. 北京大學(xué)學(xué)報(自然科學(xué)版), 2019, 55(5): 813–822

[11]張喚蘭, 王保利, 寧杰遠(yuǎn), 等. 高鐵地震數(shù)據(jù)干涉成像技術(shù)初探. 地球物理學(xué)報, 2019, 62(6): 2321–2327

[12]石永祥, 溫景充, 寧杰遠(yuǎn). 高鐵震源地下介質(zhì)成像的理論分析. 中國科學(xué): 地球科學(xué), 2022, 52(5): 893?902

[13]殷常陽, 石永祥, 伍晗, 等. 利用高鐵震源生成虛擬道集的理論分析. 地球物理學(xué)報, 出版中

[14]Mahdad A, Doulgeris P, Blacquiere G. Separation of blended data by iterative estimation and subtraction of blending interference noise. Geophysics, 2011, 76(3): Q9–Q17

[15]Wapenaar K, Neut J, Van der Thorbecke J. Deblen-ding by direct inversion. Geophysics, 2012, 77(3): A9–A12

[16]Wu S, Blacquière G, van Groenestijn G-J A. Shot repetition: an alternative seismic blending code in marine acquisition. Geophysics, 2018, 83(6): P29–P37

[17]Schuster G T, Wang X, Huang Y, et al. Theory of multisource crosstalk reduction by phase-encoded sta-tics: theory of multisource crosstalk reduction. Geo-physical Journal International, 2011, 184(3): 1289–1303

[18]Moore I, Dragoset B, Ommundsen T, et al. Simulta-neous source separation using dithered sources // SEG Technical Program Expanded Abstracts. Las Vegas, 2008: 2806–2810

[19]Clayton R, Engquist B. Absorbing boundary condi-tions for acoustic and elastic wave equations. Bulletin of the Seismological Society of America, 1977, 67(6): 1529–1540

[20]蔣一然, 鮑鐵釗, 寧杰遠(yuǎn), 等. 高架橋下方高鐵地震信號頻譜特征研究. 北京大學(xué)學(xué)報(自然科學(xué)版), 2019, 55(5): 829?838

Validation of the Virtual Gather Generation Method from Synthetic Data for High-Speed Rail Seismic Source

YIN Changyang1, SHI Yongxiang1, WU Han1, NING Jieyuan1,2,?

1. School of Earth and Space Sciences, Peking University, Beijing 100871; 2. Hebei Hongshan National Geophysical Observation and Research Station, Xingtai 054000; ? Corresponding author, E-mail: njy@pku.edu.cn

Based on the principle of superposition method which can eliminate the crosstalk noise of the high-speed rail seismic wavefield recorded by the stations and generate virtual gathers, the operation ways of the superposition method are given and the method is verified on the synthetic high-speed rail data. The precise source time function of the high-speed rail is obtained by establishing the vibration equation of the high-speed rail viaduct, and the forward modeling of the high-speed rail wavefield in layered medium is calculated by using the finite element methodand time-shift superposition of forward modeling of the high-speed rail wavefield with different speeds is carried out. The virtual gathers are obtained, which verifies the correctness of the stacking method.

high-speed rail; interference field; source time function; virtual gather

10.13209/j.0479-8023.2022.047

2021–06–01;

2022–05–22

河北省地震科技星火計劃(DZ20200827054)和國家自然科學(xué)基金(41874071)資助