台北都会区地震波传递之地形效应.PDF

台北都會區地震波傳遞之地形效應 李憲忠 中央研究院地球科學研究所 摘 要 台北都會區位在北台灣之鬆軟沈積盆地上 。過去二十幾年來，地震一直是 盆地內造成嚴重建築物損害與人員傷亡的主要天然災害，特別是 1999 年的 921 大地震(M =7.3)與 2002 年 3 月 31 日的花蓮外海地震(M =6.8) 。過去的研究顯示， L L 盆地內部的低速物質與盆地幾何形狀會影響震波能量的傳播行為並延長地振動 時間。然而目前學界對於地表地形所造成的地震波異常波傳行為的研究並不多， 特別是在真實地形之響應上。台北盆地四周為群山所環繞，對比於相對較淺的盆 地構造，盆地附近的地形對地震波傳遞所造成的影響可以預期。由於台北附近山 區之地表觀測並不多，因此透過數值計算進行理論模擬將是最直接而快速有效的 方法。此研究透過譜元素法(Spectral-Element Method, SEM)針對台北盆地附近區 域進行三維全波場彈性波模擬，研究中將模擬的重點放在台北都會區之地表地 形，分析盆地周圍之主要地表構造對地震波之傳遞可能造成的影響。 地表數值模型與模擬方法 盆地附近之地表地形由台灣地區數位地形資料(40 公尺 DEM)所建立。圖一 為台北盆附近之地表地形。主要的構造有台北盆地、北邊之大屯火山群、東南側 之雪山山脈以及西側的林口台地。相較於目前已知的盆地深度(約 1000m) ，台北 盆地周圍的群山高度可達 1000m 以上，且在短距離內有很大的高程差，在這樣 陡峭的地表起伏下，地形效應對地震波傳遞的影響有待仔細加以檢驗。目前譜元 素法在三維波傳數值模擬上最具彈性與發展潛力，然而由於模型網格的建立具相 當的複雜與困難度，特別是在地表地形的整合上，因此譜元素法在普遍應用上仍 有所困難。於台北盆地譜元素法的研究中，我們嘗試在網格中加入一個緩衝層 (Buffer Layer)用於吸收陡峭地形所造成的網格扭曲來避免品質不好的網格造成 的計算不穩定。目前我們已有能力將真實地表地形完整地建置於網格模型中，透 過這個方法結合 SEM 的特性，本研究以 40 公尺的高解析度建構台北盆地附近之 地表地形網格，並藉而提高合成地震波之穩定計算頻率至 1.5 Hz 。研究中以2004 年 10 月 23 日台北四獸山地震為模擬對象，假設其震源持續時間為 0.9 秒來模擬 1 Hz 以上的地震波傳行為，並探討地震波能量在真實地表地形構造下的傳播現 象。 模擬結果與討論 分析模擬結果中地形對最大速度值(PGV)的影響，研究中 PGV 放大值定義為： [(具地形響應之 PGV – 無地形響應之 PGV ) / 無地形響應之 PGV] * 100% 如圖二所示，大致上地表地形會造成山頂區域的 PGV 放大；相反地山谷地區則 出現 PGV 減小的現象，也就是谷地地形之 PGV 值較完全水平的地形來的小。圖 二中並展示的部分地點的合成波形，可以清楚看到山區位置的紀錄其最大速度值 比未考量地形之模型增加許多，特別是在 S 波的震幅上。由此模擬結果看來，地 形所造成的 PGV 變化值約在 ± 100% ，但若考量山峰與山谷間的最大差值則可能 達到兩倍以上。另一個值得注意的現象是，在純粹考量地形而沒有盆地內部低速 物質的情況下，台北盆地東區仍可以看到 PGV 值有少許增加的趨勢，推斷這個 現象是導因於地震所發生的位置與附近之地表地形相互作用所致。 於本研究中地表地形對地震波能量傳遞的影響是顯而易見的，特別是在局 部區域，如山頂的 PGV 放大、山谷地區的 PGV 減小、以及地形產生的散射波 造成振動持續時間延長與總累積動能增加等現象。這些地形效應在台北盆地周圍 陡峭多變的地表構造中尤為顯著。主要的 PGV 放大出現在山區，然而在沒有考 量盆地內部低速物質的情況下，受到震源位置與地形間交互作用的影響，台北盆 地東區仍出現較大的 PGV 分佈。除了部分山谷地區，整個研究區域的振動總累 加動能亦受到地形的影響而增加。這些現象均指出，台北盆地周圍複雜的地表形 貌亦是造成盆地內部異常波傳的