The Fault Geometry and Co-seismic Surface Deformation Around the Northern Taichung Basin

台灣中部地區於1999年受到集集地震之重創，造成兩千餘位民眾喪生，上千億的財物損失。造成此地震之車籠埔斷層上盤北段於同震時有著近十公尺的位移量，此等變形規模不僅為台灣百年唯一，在世界上亦不多見。由於車籠埔斷層北段有著異常巨大的變形量，因此本研究即針對該地區進行斷層面構造幾何之探討，並建立了一區域尺度之大型三維地質構造模型，以及一公里尺度的小型三維地震斷層面幾何模型，試圖找出此等規模之變形與其斷層構造之關聯，以及在此等級的上盤位移下，伴隨主斷層活動的次生構造幾何形貌與運動方式，最後使用數值方法模擬近地表斷層面上的滑移分布。 西部麓山帶於苗栗的褶皺區塊至台中的覆瓦排列斷層區塊間，於地表上有明顯的分界，大致平行於后里台地以及大安、大甲兩溪，而車籠埔斷層上盤之最大抬升地區的即座落於本轉換帶之南緣，為找出本構造區塊轉換帶與車籠埔斷層上盤抬升之關聯，本研究建立一個橫跨苗栗以及台中地區的三維地質構造模型。此模型使用八條構造平衡剖面建構，而其中一條為本研究自行製作之構造平衡剖面。該剖面之建構是依據中油公司之震測剖面及地表地層位態，配合重定位之震源控制剖面深處而製作完成。 此三維地質構造模型顯示，苗栗與台中構造區塊交界帶似乎受控於一大型向南傾斜之斷坡，橫跨苗栗與台中兩區的構造均受到此斷坡之影響，其中包括出磺坑－關刀山背斜之隱沒、三義斷層上盤大安背斜之生成、以及車籠埔斷層之走向轉折等現象均受到該斷坡之影響。由於車籠埔斷層下盤之大安背斜阻擋車籠埔斷層向北的運動，因此於該背斜處會有較大的同震抬升產生，而此區即為本次集集地震垂直抬升之最大地區。 於集集地震同震最大抬升區之內，本研究使用豐原市的都市計畫樁作為工具，分析該市東側之同震變形。結果顯示該區域內之上盤變形並非單一塊體運動，而是數個塊體以1 m內的位移差共同造成上盤的位移。各塊體依其岩性不同而呈現類似剛體或有10 ~ 20 cm的內部變形之運動模式。本研究結果顯示豐原市東側集集地震地表破裂上盤區域仍有許多次生活動構造，於集集地震時產生大小不一之位移。 本研究使用計算自都計樁之空間位移向量，建構一個公里級尺度的地震斷層構造幾何模型。該模型之形貌顯示公里級尺度之下的斷層傾角可以有相當大的變化，在短短數百公尺之內造成區域性的斷坡與斷坪。而利用鑽井資料測試本研究建構之斷層面幾何模型，可得到僅15 m的深度差異，本研究認為該模型確實可代表集集地震的同震地震斷層形貌。 最後，本研究使用此公里級之地震斷層幾何模型進行數值模擬，以了解斷層上盤於此尺度之下的地表形變模式、材料特性以及斷層面上同震滑移量分布。模擬結果顯示，斷層上盤的地表形變完全受控於其下之斷層形貌，而公里級尺度之下的斷層面上滑移量分布則顯示出量值大小亦受控於斷層面之幾何形貌。以本研斷層面滑移分布為例，滑移最大之位置即位在兩斷坡之交界處，且位在利用地表位移觀察到的一剛性運動塊體之邊界上，此現象顯示斷層面上的滑移量分布與地表觀測到的位移模式有絕對之關聯。 本論文使用數種不同之工具建構台中盆地北側於不同尺度之下的構造形貌，希望經由三維構造模型之呈現，釐清台中盆地北側之構造形貌，以及車籠埔斷層上盤活動模式與本區域之其他構造的關聯。而本論文利用數值方法討論公里級尺度之集集地震同震變形之研究，解析出了以往區域尺度相關研究無法解析之近地表淺層斷層上盤運動模式以及斷層面上滑移量的分布。希望未來有更多地球物理、鑽井以及測量資料檢測本研究建構之模型，並期望在不久的將來有更多研究結果補強本模型之不足，可以更清楚的解析出台中盆地北側的地下地質，特別是可能的發震斷層構造形貌，以期減少本區可能之災害性地震發生時所帶來之傷害。The northern Taichung basin is a transition zone between two structural domains in central Taiwan: folds-dominated domain in the north and imbricate thrusts- dominated domain in the south. Both of the north and south domains are highly active area. Series of large earthquakes occurred in 1935 in the north domain, and the devastating Mw 7.6 Chi-Chi earthquake occurred in 1999 in the south domain. The main purposes of this study are: to construct the structure geometry in the highly active structural transition area; to reveal detailed deformation pattern in the hanging wall of Chi-Chi earthquake and to construct a kilometer scale geometry model of the earthquake fault; and to simulate the slip distribution on the proposed fault plane. A 3D geological model was established based on seven published balanced cross sections plus one cross section constructed in the first part of the thesis. The constructed cross section is based on the surface geological data and relocated earthquake hypocenters. This 3D geological model indicated the high co-seismic displacements in northern section of Chi-Chi earthquake fault is result from a fault- bend fold in its footwall. The fold is caused by a lateral ramp underneath the Daan and Dacha Rivers. The anticlinal fold is also responsible for the changing direction of the rupture trace in the northern Chi-Chi earthquake fault. The second part is to characterize the deformational behaviors of the hanging wall block due to earthquake fault slip in northern Chi-Chi earthquake fault. In the case of the 1999 Chi-Chi earthquake, hanging wall block of the earthquake fault showed complex deformation pattern at the kilometer scale. Because previous studies mainly characterize the fault at the regional scale, it is of interest and also a challenge to characterize the fault at a smaller scale with a higher resolution. In this study, the geometry of a kilometer-scale patch of the fault plane is reconstructed using the displacement data collected from the densely distributed city planning benchmarks. The study area is approximately 4 km by 8 km in size, and contains as many as 924 benchmarks. Among the benchmarks, 62 have both horizontal and vertical displacement data, and the rest of the benchmarks have the horizontal displacement data. Based on the assumption of constant volume, the earthquake fault geometry is built by using the 62 slip vectors. The derived fault geometry model is rather consistent with the borehole data from the nearby 450 m well. The last part of the thesis applies numerical model to evaluate the fault geometry and to simulate the slip distribution on the fault plane. The fault geometry can be evaluated using the concept of fault parallel flow. The method of fault parallel flow can restore the topography before the earthquake and the deformation controlled by the fault geometry. Comparing pre-earthquake and post-earthquake cross sections of topography shows very high correlation coefficient, above 0.99 in average. Half-space elastic dislocation model is then used to simulate the slip distribution on the fault plane, which is inversed from the displacement data of city planning benchmarks. The result shows significant slip variation at the kilometer scale and is interpreted to be related to the deformation zone. The thesis uses different methods to reveal the spatial distribution and 3D geometry of seismogenic structures in northern Taichung basin at different scale. The results indicate better resolution of the seismogenic can apply to the earthquake hazards mitigation in this area in the future.目 錄 口試委員會審定書...........................................I 致 謝....................................................II 摘 要...................................................III Abstract...................................................V 目 錄...................................................VII 圖目錄.....................................................X 表目錄...................................................XII 第一章 序論................................................1 1.1 前言...................................................1 1.2台中盆地北側三維地質構造模型研究概述....................1 1.3上盤同震位移與公里級尺度地震斷層形貌建構概述............2 1.4上盤同震變位與斷層面位移分布數值模擬概述................3 第二章 台中盆地北側三維地質構造模型........................4 2.1 前言...................................................4 2.2 苗栗與台中地區之構造概述...............................6 2.2.1 苗栗地區構造概述.....................................6 2.2.2 台中地區構造概述.....................................7 2.3構造平衡剖面之建構與回復................................8 2.3.1構造平衡剖面之建構工具及方法簡述......................8 2.3.2構造平衡剖面之回復....................................9 2.3.3建構大甲－后里－白冷之地質構造平衡剖面...............11 2.3.4 不同基底滑移面形貌對構造平衡剖面建構之比較..........14 2.4 苗栗與台中構造區塊交界帶之三維地下構造模型............15 2.4.1研究區域內已發表之構造平衡剖面.......................15 2.4.2研究區域內之立體構造模型.............................19 2.5 小結..................................................31 第三章 地震斷層上盤同震位移及地震斷層面建構－以豐原市東側集 集地震斷層為例.....................................32 3.1 前言..................................................32 3.2前人研究...............................................33 3.2.1研究區域地形與地質概述...............................33 地形概述.................................................33 地質概述.................................................38 3.2.2集集地震大尺度同震地表變形相關研究...................40 3.2.3集集地震斷層北段之野外測量地表變形相關研究...........41 3.3研究使用之都市計劃樁介紹...............................43 3.4豐原市東側集集地震同震變形模式.........................46 3.4.1橫跨上下盤之水平變形模式.............................46 3.4.2水平位移場以及上盤內部之水平方向相對變形.............51 3.4.3垂直位移場...........................................56 3.5豐原市東側集集地震斷層幾何模型之建構與測試.............58 3.6討論：集集地震斷層之上盤運動模式.......................64 3.7小結...................................................66 第四章 同震地形變位與斷層面位移分布數值模擬...............67 4.1前言...................................................67 4.1.1 FPF回復法簡介.......................................67 4.1.2彈性半空間斷層錯位模型簡介...........................71 4.2地震斷層幾何模型之數值模擬結果.........................75 4.2.1以FPF法模擬上盤之同震地表形變........................75 4.2.2以斷層錯位法評估研究區域之地質材料特性...............80 4.2.3以斷層錯位法模擬分布於斷層面之同震滑移量.............86 4.3討論：影響上盤運動模擬之因素...........................87 4.3.1影響FPF法模擬地表形變結果之因素......................88 4.3.2地質材料特性與彈性半空間斷層錯位模擬結果之關聯.......88 4.4小結...................................................90 第五章 結論...............................................91 參考文獻..................................................93 中文部分..................................................93 地質圖....................................................93 期刊、專書及報告..........................................93 英文部分..................................................95 附錄A 都市計畫樁測量誤差.................................102 附錄B 研究區域北段於三種不同包松比環境下模擬都計樁位移值一覽 表.................................................104 圖 目 錄 圖2-1研究區域地質圖........................................5 圖2-2苗栗台中地區地下斷層模型立體圖........................6 圖2-3曲滑褶皺恢復平衡法地層回復示意圖......................9 圖2-4楊耿明等人所製作之Line 20剖面........................10 圖2-5楊耿明等人所製作之Line 20剖面回復後之形貌............10 圖2-6基底滑移面於10 km以下之構造平衡剖面及剖面回復圖......12 圖2-7基底滑移面於10 km以上之構造平衡剖面及剖面回復圖......13 圖2-8 Namson所製作之Line 1剖面............................15 圖2-9楊耿明等人所製作之22號剖面...........................16 圖2-10楊耿明等人所製作之21號剖面..........................17 圖2-11楊耿明等人所製作之20號剖面..........................18 圖2-12洪日豪與Suppe所製作之剖面...........................18 圖2-13楊耿明等人所製作之18號剖面..........................19 圖2-14研究區域立體構造模型俯視圖..........................21 圖2-15研究區域立體構造與地層模型俯視圖。..................22 圖2-16模型內之震源分布....................................23 圖2-17南北走向之模型剖面與解釋............................24 圖2-18三義斷層下盤木山層層頂模型圖........................25 圖2-19三義斷層下盤東坑層層頂模型圖........................25 圖2-20三義斷層下盤桂竹林層層頂模型圖......................26 圖2-21三義斷層下盤錦水頁岩層頂模型圖......................26 圖2-22三義斷層下盤卓蘭層層頂模型圖........................27 圖2-23三義斷層以及其上盤木山層層頂模型圖..................28 圖2-24三義斷層上盤東坑層層頂模型圖........................28 圖2-25三義斷層上盤桂竹林層層頂模型圖......................29 圖2-26車籠埔斷層模型圖....................................29 圖2-27三義斷層上盤錦水頁岩層頂模型圖......................30 圖2-28三義斷層上盤卓蘭層層頂模型圖........................30 圖3-1台灣中部地區地形與地質簡圖...........................34 圖3-2研究區域地形分區簡圖.................................35 圖3-3豐原市水源路斷層崖照片...............................36 圖3-4豐原市中正公園內的單斜崖.............................36 圖3-5南陽路階地旁於集集地震同震時產生的斷層崖.............37 圖3-6研究區域地質簡圖.....................................39 圖3-7集集地震六年後拍攝野外露頭之相片.....................39 圖3-8豐原市千分之一都計圖.................................45 圖3-9a集集地震前後都計圖疊合比較圖........................46 圖3-9b豐原市東側都計圖重測區域圖..........................46 圖3-10水源路測線西向與北向水平位移分量圖..................48 圖3-11南陽路測線西向與北向水平位移分量圖..................49 圖3-12南陽路階地構造分布圖................................51 圖3-13研究區域上盤水平位移向量圖與角度、淨向量長度直方圖..53 圖3-14研究區域都計樁位移場圖..............................54 圖3-15上盤相對位移向量圖..................................57 圖3-16地震斷層幾何面製作流程..............................59 圖3-17立體地震斷層幾何模型................................60 圖3-18集集地震地震斷層模型等深圖..........................62 圖3-19地震斷層鑽井測試示意圖..............................63 圖4-1 FPF法回復地層示意圖.................................69 圖4-2角狀錯動原理示意圖...................................73 圖4-3順推位移觀測點與斷層面幾何關係示意圖.................75 圖4-4研究區域北段剖面位置圖...............................76 圖4-5未回復之地形以及斷層剖面圖...........................77 圖4-6回復震後地形剖面與震前地形剖面比較圖.................77 圖4-7研究區域北段斷層面逆推用向量分布圖...................82 圖4-8包松比與六個滑移方向的最大位移量之關聯圖.............82 圖4-9模形於不同包松比環境下模擬出之向量與實測值誤差數量分布 直方圖...............................................85 圖4-10順推模擬結果與觀測值比較圖..........................86 圖4-11位移向量逆推位移量至斷層面上之滑移分布圖............87 表 目 錄 表4-1 研究區域內上盤滑移點淨滑移長數量表..................78 表4-2 研究區域北段上盤地形剖面回復長度、單剪角度與相關係數關 係表................................................79 表4-3 Poly3D測試結果......................................81 表4-4模形於不同包松比環境下模擬出之向量與實測值誤差數量分布 圖..................................................8