彈性與非彈性模型參數擾動對於走時與振幅異常敏感度算核之影響

以作者查詢圖書館館藏

、以作者查詢臺灣博碩士

、以作者查詢全國書目

、勘誤回報

、線上人數：10

、訪客IP：18.218.159.7

姓名

謝銘哲(Ming-Che Hsieh) 查詢紙本館藏

畢業系所

地球物理研究所

論文名稱

彈性與非彈性模型參數擾動對於走時與振幅異常敏感度算核之影響
(Study on the sensitivities of seismic traveltimes and amplitudes to elastic and anelastic model perturbations)

相關論文

★ 台灣地區中大型地震震源參數分析	★ 台灣北部地區之隱沒樣貌
★ 九二一集集地震之餘震(Mw≧6.0)震源破裂滑移分佈	★ 利用雙差分地震定位演算法重新定位過去十年台灣中、大型地震之餘震
★ 九二一集集地震三維震源過程與震波傳遞分析	★ 台灣弧陸碰撞構造之地殼及頂部地函的三維S波衰減模型
★ 集集地震之震前、同震及震後變形模式研究	★ 台灣地震震源尺度分析:2003年規模>6.0地震分析
★ 使用震源機制逆推台灣地區應力分區狀況	★ 地震水井水力學之理論模式改良與發展及同震水位資料分析
★ 台灣東北部外海地震之三維強地動模擬	★ 利用臨時寬頻地震網觀測嘉義地區淺層地殼之非均向性
★ 中大規模地震斷層參數之同步求解	★ 集集地震同震及震後應力演化與地震活動之相關性
★ 2005 年宜蘭雙主震之震源破裂滑移分析	★ 1999 集集地震後之黏彈性鬆弛效應

檔案

[Endnote RIS 格式]

[Bibtex 格式]

[相關文章]

[文章引用]

[完整記錄]

[館藏目錄]

[檢視]

[下載]

本電子論文使用權限為同意立即開放。
已達開放權限電子全文僅授權使用者為學術研究之目的，進行個人非營利性質之檢索、閱讀、列印。
請遵守中華民國著作權法之相關規定，切勿任意重製、散佈、改作、轉貼、播送，以免觸法。

摘要(中)

射線理論(Ray theory)為對震波作高頻近似假設，經常應用於三維速度構造逆推問題，該方法雖簡化問題且易於實行，但往往受限於兩個因素：(1)目標模型的解析力受限於菲涅爾帶(Fresnel zone)寬度；(2)因為高頻近似假設，致使觀測站所收集波形資訊完全無法應用，使大量寶貴資料被捨棄。然而採有限頻寬方法(finite-frequency approach)透過數值模擬與全波場(full-waveform)計算，我們可以計算參考模型下的理論地震圖，藉由比較理論與觀測地震圖二者差異，進而定量化到時異常(traveltime anomaly)與振幅異常(amplitude anomaly)等波形異常資訊，爾後經由波形異常資訊計算構造異質性對於有限頻寬(finite-frequency)震波的影響，即為敏感度算核。
　　本論文應用應變格林張量(strain Green’s tensor, or SGT)與三維交錯網格有限差分法(3-D staggered-grid finite-difference method)波形模擬，藉此探討模型參數擾動對於有限頻寬震波的影響及其敏感度算核。模型擾動可為彈性波速度擾動或非彈性介質之品質因子擾動，代表介質擾動對於波形的影響，並透過震源－觀測站的互換性(reciprocity)使用以測站位置為震源的應變格林張量，計算參考模型下震源－觀測站間的理論地震圖與到時和振幅異常的敏感度算核。首先在均質模型下計算速度擾動的敏感度算核，其震源－觀測站間的射線路徑為一直線，但實際震波為有限頻寬，因此速度擾動造成到時差異的算核為一中空橢球，或稱香蕉甜甜圈(banana-doughnut)，代表觀測站所收到的波形到時隨模型擾動之影響，不單單只是射線路徑(raypath)所穿越之速度構造所貢獻，而是含蓋更廣泛的區域。而品質因子的擾動可藉由該時間視窗的波形的希爾伯特轉換(Hilbert transform)表示，算核形貌如同前述為橢球型，但正負符號相異，即品質因子對於波形異常的貢獻與彈性波範疇下速度擾動大相逕庭。此外，震源機制與濾波器亦為影響算核形貌與分佈的因素，結果顯示能量輻射與散射深受上述物理機制影響。
　　若在三維模型下計算敏感度算核，則明顯可見在速度愈高的區域，其敏感度算核愈寬、影響波形的範圍愈廣，且相較於速度模型的側向變化，深度方向速度梯度較能影響算核形貌。若將速度擾動的算核推廣至非彈性衰減擾動，則同樣可見類似形貌的算核。

摘要(英)

Ray theory is based on infinitely-high frequency approach and has been widely used in seismic tomography. It is convenient and is therefore commonly used in seismology. However, ray theory tomography is constrained by two important shortcomings: (1) the resolving power of the tomography model is limited by the width of the first Fresnel zone; (2) the full waveform information (e.g. the time history of the seismic response) is discarded because of the high frequency assumption. The recently developed finite-frequency tomography approach provides a way of fully utilizing the spatial and temporal responses of seismic wavefield. The finite-frequency approach adopts numerical or semi-analytical methods to compute the Green’s functions in reference models. Sensitivity kernels can be expressed in terms of the waveform perturbations due to point sources generated by structural heterogeneities, and these waveform perturbations result in traveltime and amplitude anomalies.
　　In this thesis, we apply the strain Green’s tensor (SGT) approach and 4th-order staggered-grid finite-difference method to solve the wave equations and calculate the sensitivity kernels for both elastic and anelastic perturbations of the medium properties. The ellipsoid and hollow appearance (in the shape of banana-doughnut) of the traveltime anomaly kernels for velocity (νp and νs) perturbations will be shown first in the homogeneous elastic medium for each phase. Sensitivities are close to zero on the raypath and strongly influenced by the free surface. This phenomenon is due to finite-frequency effect on the full-waveform that cannot be seen in ray theory. In addition, the anelastic sensitivity kernels of the quality factor (Qp and Qs) perturbations can be expressed in terms of the original waveform’s Hilbert transform. The shapes of kernels are similar to those for the elastic parameters, but the patterns of kernels in terms of sign changes indicate that the contributions to the waveform by the anelastic perturbation are quite different than velocity perturbations of the reference models. We also discuss other effects, such as filtering and focal mechanism, on the sensitivity kernels. The results show that these effects have great influence on the sensitivity kernel.Finally, we demonstrate the same calculation procedure of a scenario earthquake in southern Taiwan in a 3-D velocity model. The result shows that the appearance of sensitivity kernels strongly depends on vertical velocity gradients. Unlike in homogeneous model, the traveltime anomaly kernel is generally non-zero on the raypath, which can be explained by the multi-pathing effect of seismic wave propagation in a heterogeneous structure.

關鍵字(中)

★ 應變格林張量
★ 敏感度算核
★ 有限差分法

關鍵字(英)

★ sensitivity kernel
★ finite-difference method
★ strain Green's tensor

論文目次

中文摘要...............................................i
英文摘要.............................................iii
致謝..................................................iv
目錄...................................................v
圖目錄...............................................vii
第一章緒論............................................1
1.1 研究動機與目的.....................................1
1.2 文獻回顧...........................................2
1.3 本文範疇與內容.....................................4
第二章研究方法與原理..................................6
2.1 量化波形異常.......................................6
2.2 敏感度算核.........................................8
2.3 應變格林張量......................................13
2.4 以應變格林張量表示敏感度算核......................15
2.5 品質因子擾動之敏感度算核..........................16
2.6 三維有限差分法波傳模擬............................19
2.6.1交錯網格有限差分法...............................20
2.6.2 震源參數........................................22
2.6.3 吸收邊界條件....................................23
2.6.4 以有限差分法計算應變格林張量....................24
第三章參考模型中敏感度算核評估.......................30
3.1 均質模型建置與波傳模擬............................30
3.2 時間偏移量修正....................................31
3.3 均質模型中的敏感度算核............................33
3.3.1 P波相...........................................33
3.3.2 pP波相..........................................34
3.3.3 pS波相..........................................35
3.4 品質因子擾動之敏感度算核實例......................36
第四章濾波器、震源機制、速度模型對於敏感度算核影響...58
4.1 濾波器效應........................................58
4.2 震源機制..........................................59
4.3 三維速度模型中的敏感度算核........................61
第五章結論...........................................86
參考文獻..............................................88

參考文獻

Aki, K., and P. G. Richards. Quantitative Seismology, Second Ed., University Science Books, Sausalito, California, 700 pp., 2002.
Cerjan, C., D. Kosloff, R. Kosloff, and M. Reshef, A non-reflecting boundary condition for direct acoustic and elastic wave equations, Geophys., 50, 705-708, 1985.
Chen, P., A unified methodology for seismic waveform analysis and inversion, Ph.D. Thesis, Univ. of Southern California, Los Angeles, California, 201 pp., 2005.
Chen, W.-F., and A. F. Saleeb, Constitutive Equations for Engineering Materials, Revised Ed., Elsevier, 1096pp., 1994.
Clayton, R., and B. Engquist, Absorbing boundary conditions for acoustic and elastic wave equations, Bull. Seismol. Soc. Am., 67, 1529-1540, 1977.
Dahlen, F. A., Finite-frequency sensitivity kernels for boundary topography perturbations, Geophys. J. Int., 162, 525-540, 2005.
Dahlen, F. A., and J. Tromp, Theoretical Global Seismology, Princeton University Press, Princeton, 944pp., 1998.
Dahlen, F. A., S.-H. Hung, and G. Nolet, Fréchet kernels for finite-frequency traveltimes, I. Theory, Geophys. J. Int., 141, 157-174, 2000.
Dahlen, F. A., and Y. Zhou, Surface-wave group-delay and attenuation kernels, Geophys. J. Int., 165, 545-554, 2006.
Furumura, T., B. L. N. Kennett, and H. Takenaka, Parallel 3-D pseudospectral simulation of wave propagation by using a Connection machine and workstation-clusters, Geophys., 63, 279-288, 1998.
Graves, R. W., Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences, Bull. Seism. Soc. Am., 86, 1091–1106, 1996.
Gee, L. S., and T. H. Jordan, Generalized seismological data functionals, Geophys. J. Int., 111, 363-390., 1992.
Hung, S.-H., F. A. Dahlen, and G. Nolet., Fréchet kernels for finite-frequency traveltimes, II. Examples, Geophys. J. Int., 141, 175-203, 2000.
Kim, K. H., J. M. Chiu, J. Pujol, K. C. Chen, B. S. Huang, Y. H. Yeh, and P. Shen, Three-dimensional Vp and Vs structural model associated with the active subduction and collision tectonics in the Taiwan region, Geophys. J. Int., 162, 204-220., 2005.
Komatitsch, D., and J. P. Vilotte, The spectral-element method: an efficient tool to simulate the seismic reponse of 2-D and 3D geological structures, Bull. seism. Soc. Am., 88(2), 368–392., 1998.
Lay, T. and T. C. Wallace, Modern Global Seismology, Academic Press, San Diego, California, 535pp., 1995.
Levander, A. R., Fourth-order finite-difference P-SV seismograms, Geophys., 53, 1425-1436, 1988.
Liu, Q., and J. Tromp, Finite-frequency sensitivity kernels for global seismic wave propagation based upon adjoint methods, Geophys. J. Int., 174, 265-286, 2008.
Liu, H.-P., D. L. Anderson, and H. Kanamori, Velocity dispersion due to anelasticity: implications for seismology and mantle composition, Geophys. J. R. Astr. Soc. 47, 41-58., 1976.
Ma, K. F., J. H. Wang, and D. Zhao, Three-dimensional seismic velocity structure of the crust and uppermost mantle beneath Taiwan, J. Phys. Earth, 44, 85-105, 1996.
Magistrale, H., S. Day, R. W. Clayton, and R. Graves, The SCEC Southern California reference three-dimensional seismic velocity model Version 2, Bull. Seism. Soc. Am., 90, S65-S76, 2000.
Marcinkovich, C., and K. Olsen, On the implementation of perfectly matched layers in a three-dimensional fourth-order velocity-stress finite difference scheme, J. Geophys. Res., 108(B5), 2276, 2003.
Marquering, H., F. A. Dahlen, and G. Nolet, Three-dimensional sensitivity kernels for finite-frequency traveltimes: the banana-doughnut paradox, Geophys. J. Int., 137, 805-815, 1999.
Olsen, K. B., Simulation of three-dimensional wave propagation in the Salt Lake Basin, Ph.D. Thesis, Univ. of Utah, Salt Lake City, 157 pp., 1994.
Rau, R.-J., and F. T. Wu, Tomographic imaging of lithospheric structures under Taiwan, Earth Planet. Sci. Lett., 133, 517- 532, 1995.
Roecker, S. W., Y. H. Yeh, and Y. B. Tsai, Three-dimensional P and S wave velocity structures beneath Taiwan; deep structure beneath an arccontinent collision, J. Geophys. Res., 92, 10547-10570, 1987.
Tromp, J., C. Tape, and Q. Liu, Seismic tomography, adjoint methods, time reversal, and banana-doughnut kernels, Geophys. J. Int., 160, 195-216, 2005.
Virieux, J., P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method, Geophys., 51, 889– 901, 1986.
Wu, Y.-M., C.-H. Chang, L. Zhao, J. B. H. Shyu, Y.-G. Chen, K. Sieh, and J.-P. Avouac, Seismic tomography of Taiwan: Improved constraints from a dense network of strong motion stations, J. Geophys. Res., 112, B08312, 2007.
Zhang, Z., Y. Shen, and L. Zhao, Finite-frequency sensitivity kernels for head waves, Geophys. J. Int., 171, 847-856, 2007.
Zhao, L., T. H. Jordan, and C. H. Chapman, Three-dimensional Fréchet differential kernels for seismic delay times, Geophys. J. Int., 141, 558-576, 2000.
Zhao, L., T. H. Jordan, K. B. Olsen, and P. Chen, Fréchet kernels for imaging regional earth structure based on three-dimensional reference models, Bull. Seism. Soc. Am., 95, 2066-2080, 2005.
Zhao, L., and T. H. Jordan, Structural sensitivities of finite-frequency seismic waves: a full-wave approach, Geophys. J. Int., 165, 981-990, 2006.
Zhao, L., P. Chen, and T. H. Jordan, Strain Green’s tensors, reciprocity, and their applications to seismic source and structure studies, Bull. Seism. Soc. Am., 96, 1753-1763, 2006.
趙里，知識天地－全波場波形分析：地震學研究的新方向，中央研究院週報，第1115期，2007。
楊欣穎，模擬波與波線走時：有限頻寬與波線層析成像的極限，國立台灣大學地質科學研究所碩士論文，台灣台北，57頁，2005。
王乾盈，台灣大地動力學國際合作研究計畫（TAIGER）陸上炸測野外工作手冊，國立中央大學地球物理研究所淺層反射震測實驗室，台灣中壢，35頁，2008。

指導教授

馬國鳳、趙里
(Kuo-Fong Ma、Li Zhao)

審核日期

2009-8-23

推文