梁建文 劉中憲 巴振寧 任葉飛

1） 中國天津 300350 天津大學建筑工程學院

2） 中國天津 300384 天津城建大學土木工程學院

3） 中國哈爾濱 150080 中國地震局工程力學研究所

強地震動模擬分析涉及地震學、巖石力學、土動力學等多個學科，是工程地震學的核心內容之一，國內外學者在該領域已取得較多重要進展.近年來，隨著國家韌性城鄉(xiāng)建設的需要，區(qū)域復雜場地地震區(qū)劃及重大工程建設對地震動場的準確模擬提出了更高要求.綜合考慮震源、路徑和場地效應的區(qū)域尺度地震動模擬以及場地-結構（群）地震相互作用等已成為當前研究的熱點和難點.隨著計算機軟硬件技術的飛速提升及三維地質信息、地震觀測數(shù)據(jù)的積累，目前已可實現(xiàn)從斷層動力破裂到城市建筑群地震破壞的全過程模擬（陸新征，2020）.

基于物理模型的數(shù)值模擬方面，有限差分法、有限元法、譜元法、邊界元法等數(shù)值方法取得了長足的進步，對實際巖土層復雜的介質特性和幾何特征的適應性越發(fā)良好（賀春暉等，2017）；針對近斷層地震動模擬，考慮實際復雜幾何特征及局部參數(shù)的斷層運動學或動力學模型逐步完善.由于核電站等抗震設計對地震動頻帶的要求較高，近年來地震動模擬有效頻率范圍也在不斷提升（Poursartipet al，2020）.各類數(shù)值模擬方法均有其優(yōu)勢和局限性，仍需不斷發(fā)展完善，如有限元、有限差分法在處理高頻波動時存在的數(shù)值頻散問題，譜元法在處理跨尺度、復雜幾何形狀等問題中存在的不均勻網(wǎng)格協(xié)調性問題，以及邊界元方法存在的介質非線性的高效積分難題等.另外，由于地震動形成的全過程存在不確定性，隨機有限斷層法、經(jīng)驗格林函數(shù)法等隨機方法也在不斷發(fā)展.綜合確定性數(shù)值法和隨機法的優(yōu)勢，采用隨機有限斷層法模擬高頻地震動、數(shù)值法模擬低頻地震動的確定性與隨機性混合的模擬方法在區(qū)域尺度場地寬頻地震動模擬中被廣泛應用.

國內外有大量城市坐落于高地震烈度的沉積盆地或河谷之中，盆地或河谷的地震動放大效應研究已取得豐碩成果（Triffunac，2016）.但對城市直下型地震，對于坐落在沉積盆地或河谷之中的城市而言，需考慮近斷層效應和盆地或河谷效應的耦合作用，極大地增加了問題的復雜性，如何進一步提高現(xiàn)有數(shù)值方法的有效模擬頻率、計算效率、穩(wěn)定性，以及科學考慮盆地表層松軟土層非線性效應，還需要進行地震學、工程學、計算科學等的跨學科探索與攻關.

場地-結構（群）相互作用方面，目前關于區(qū)域場地如盆地、山體地形等與建筑群的地震相互作用已有較多進展，揭示了密集建筑群對地震動的明顯影響.另外，城市地下空間大規(guī)模開發(fā)利用，地下結構（群）對地震動的影響及臨近結構物抗震問題亦引起廣泛關注，特別是超大型地下結構（群）對地震波的散射效應需充分考慮結構形式、上覆土層厚度以及鄰近結構相互作用等因素.

此外，隨著近幾年人工智能技術的推進，基于機器學習的地震動模擬與預測成為工程地震學領域的一個重要探索方向.典型應用為采用實測數(shù)據(jù)或數(shù)值模擬結果構建數(shù)據(jù)集，借助人工神經(jīng)網(wǎng)絡等構建地震動預測方程.為進一步揭示地震動形成的全過程物理機制，機器學習算法在近斷層地震動和寬頻地震動模擬中也有所發(fā)展（Khosravikia，Clayton，2021）.

為了集中展示強地震動模擬及工程應用的最新進展，特組織出版 《 強地震動模擬及工程應用專輯 》 .本專輯論文共16 篇，研究內容主要涵蓋斷層效應模擬、場地地震響應分析、地震動模擬方法以及場地-結構（群）相互作用等方面，具體如下：

1） 斷層效應模擬

盆地直下型斷層發(fā)震常導致坐落于盆地之上的城市發(fā)生嚴重震害.于彥彥等采用譜元法模擬不同位置直下型有限斷層破裂時盆地內的強地震動分布場，重點研究了斷層位置對三維沉積盆地地震動影響的機理，可為盆地區(qū)域地震災害風險評估提供參考.

為拓展單臺工作站模擬大尺度復雜場地地震動的能力，巴振寧等基于CUDA 編程平臺提出了工作站級的CPU-GPU 并行異構方法.研究表明，CPU-GPU 異構并行方法可以有效地提高譜元法模擬的計算效率，具有良好的應用前景.

2） 場地地震響應分析

高玉峰等基于1992 年在我國臺灣翡翠河谷上觀測到的六條地形影響臺陣記錄，利用線源SH 波入射下非對稱V 形河谷地震波傳播解析理論，模擬了河谷臺陣各點的地震動，從理論上揭示了翡翠河谷地形效應對地震動影響的規(guī)律.

劉啟方和陳長龍利用模擬和觀測記錄分析了場地非線性時變特征，研究表明，相較于移動窗譜比法，基于移動窗解卷積法可以較好地揭示場地非線性隨震動水平的變化過程，識別非線性發(fā)生的閾值、非線性的變化程度及強震動后的恢復程度.

楊彩紅等給出了圓弧形凹陷在SH 波入射下瞬態(tài)反應的解析解，研究了凹陷在脈沖輸入與實際地震動輸入下瞬態(tài)反應的基本特征，并重點分析了凹陷對地震動的放大效應，該研究對于深入認識圓弧形凹陷對地震動放大效應具有參考價值.

韓璐等以1979 年帝王谷地震地震為例，分別使用小波變換法、峰點法和零點法等三種方法識別近斷層強震觀測記錄的速度脈沖，并借助FK 法合成的地震動，探討了近斷層脈沖型地震動的產(chǎn)生機制.

孟思博等以半圓形、V 形沉積河谷為例，以入射波條件、沉積內外介質屬性和河谷形狀為特征參數(shù)，以沉積河谷地震動放大系數(shù)為預測目標參數(shù)，建立了基于差分進化—人工神經(jīng)網(wǎng)絡的沉積河谷地震動放大效應預測模型，具有工程參考價值.

3） 地震動模擬方法

波動輸入是場地地震反應分析的重要環(huán)節(jié)，由于近海工程場地存在向海洋傾斜的單體邊坡，且邊坡一側直接與海水耦合，其特殊的場地條件導致難以直接求解得到自由場分布，對現(xiàn)有的基于自由波場的場地地震波動輸入方法提出了挑戰(zhàn).為此，寶鑫等提出了一種基于混合波場地震動輸入技術的近海場地地震反應分析方法，解決了這一工程難題.

邢浩潔等將其最近發(fā)展的多人工波速優(yōu)化透射邊界應用于高精度譜元法的地震波動模擬當中，通過與經(jīng)典的人工邊界進行比較分析，證明該方法是一種便捷、高效的人工邊界條件實現(xiàn)手段，該研究為實現(xiàn)簡單高效的地震波動模擬提供了一種新途徑.

李鴻晶和王競雄研究了構建時域譜單元質量特性模型的數(shù)學機制，針對時域切比雪夫譜單元和勒讓德譜單元建立了一種直接導出譜單元一致質量矩陣和集中質量矩陣的統(tǒng)一數(shù)學方法，從物理角度探討了譜單元質量特性模型的合理性.

地震動合成方法的應用前景與其所能表達的有效頻帶范圍緊密相關.曹澤林等研究基于頻率波數(shù)域格林函數(shù)的地震動合成FK 方法，分析了頻帶范圍的關鍵影響要素及處理措施.該研究表明，地震動合成FK 方法具有合成寬頻帶地震動的能力.

王競雄等提出了一種求解水平成層場地地震反應的時域集中質量切比雪夫譜元法，通過節(jié)點積分法嚴格地導出集中質量矩陣，克服了傳統(tǒng)切比雪夫譜元法由于質量矩陣為非對角矩陣形式所帶來的計算效率不高的問題，具有良好應用前景.

4） 場地-結構（群）相互作用

針對地下結構抗震設計中埋深這一重要因素，趙密等從一維場地地震響應進行分析，研究了面向地下結構的最優(yōu)地震動峰值指標隨埋深變化規(guī)律，確定不同埋深下的最優(yōu)地震動強度指標，可為結構抗震性能評價提供合理的地震動指標參考.

董瑞等給出了一種模擬地下結構和地面結構耦聯(lián)災變過程的模擬方法，并通過對地鐵車站-橋梁相互作用模型地震災變過程的仿真模擬，給出了地鐵車站和鄰近橋梁的耦聯(lián)災變機理，該研究對于地下結構選址及抗震設計具有參考價值.

輸入地震動的選取是地下結構抗震設計與分析的首要環(huán)節(jié)，禹海濤等為探討最不利地震動在地下結構抗震設計中的適用性，開展了軟土隧道基于地震響應的輸入地震動排序研究，該研究對軟土隧道的輸入地震動選擇具有參考價值.

禹海濤等基于Biot 飽和兩相介質理論，采用間接邊界積分方程法分別建立水平和豎向雙線隧道動力作用分析模型，研究了SV 波垂直入射下飽和地層淺埋平行隧道動力響應機制，該研究對飽和地層淺埋平行隧道抗震設計具有參考價值.

值此專輯出版之際，感謝審稿專家對本專輯文章提出的寶貴意見，感謝各位作者給予的大力支持.