跳到主要內容

簡易檢索 / 詳目顯示

研究生: 王顥鈞
Haw-Chun Wang
論文名稱: Pre-Stack Diffraction Stack Depth Migration of Active Source Short-offset Marine and Long-offset Seismic Data
指導教授: 陳浩維
How-Wei Chen
口試委員:
學位類別: 碩士
Master
系所名稱: 地球科學學院 - 地球科學學系
Department of Earth Sciences
論文出版年: 2019
畢業學年度: 107
語文別: 英文
論文頁數: 185
中文關鍵詞: TAIGER震測資料偏移疊前深度偏移
外文關鍵詞: Depth Migration, Diffraction Sum, Fresnel Zone
相關次數: 點閱:16下載:0
分享至:
查詢本校圖書館目錄 查詢臺灣博碩士論文知識加值系統 勘誤回報
  • 震波在通過不同阻抗之材質時,波傳速度的改變會體現於震測資料的走時資訊。 震測資料偏移(seismic migration)即是以此對波形資料進行重定位與回歸,所得成像為材質阻抗的在空間中的幾何分布,其有助人們判斷地下之材質岩性、構造,或是作為進一步逆推分析之初始模型以及用於偏移速度分析(Migration Velocity Analysis, MVA)等。
    本研究使用的疊前繞射疊加深度偏移(PreSDSDM)具有時間偏移法之部分優點,此方法對速度資訊之誤差有較高的容忍度,且理論上亦不需針對高程進行靜態修正,也不需顧慮資料之聚排方式,適用三維測線,且輸出成果為深度域,不須再將成像結果由時間域轉深度域,是為一種具有潛力的偏移方式。
    以此為出發點,本研究將對此方法進行原理分析、測試以及實際應用。所選用的震測資料分別為TAIGER計畫中收集的長支距陸上炸測與海上多頻道反射震測(MCS)資料。海上資料測線為TAIGER 2009 MGL0908-TST,資料量足夠且炸測點位相當密集;而陸上震測資料則是TAIGER 2008於南部Transect T4b測線施作之廣角反射震測測線,是一條橫跨133公里蜿蜒測線,有四個炸點與609個TEXAN單分量接收器,成像條件相對嚴苛。陸上資料所使用之速度模型參考自Kuo-Chen et al.(2012)的Local及Regional速度模型,以及Huang et al.(2014)所提出之速度模型,另有Roma (2018)以tau-p波場逆推所得之一維速度。
    研究成果方面,海上資料的應用非常成功,可以觀察到Lin et al. (2012)等人所提出的部分構造,並於海面下1.7公里處附近發現明顯的海底仿擬反射(BSR),或有甲烷水合物蘊藏於該地之可能。
    陸上資料方面,受限於資料的數量及品質,炸點間距過大,淺部反射訊號難以收斂。另外西部平原區淺層沖積物造成之Multiple Diving Wave (Huang et al., 2013),也對研究構成一定的難度。目前以挑選中遠支距初達波進行偏移之結果,可以觀察到西部平原區基盤之明顯象徵,以及炸點SP2東側疑似觸口斷層之構造,惟炸點SP3與SP4之間接收器密度不足,難以判斷中央山脈以下之構造。此外,使用Kuo-Chen et al.(2012)及Huang et al.(2014)速度模型進行偏移之結果,發現淺層成像收斂均出現”X”型圖樣,意味著淺層速度有高估的可能。
    本研究所用之方法仍有許多成長空間,將來可朝使用天然震源進行偏移作為研究方向,亦或是引入三維速度模型,以三維的測線幾何作為成像考量等等。


    The Pre-Stack Diffraction-Stack Depth Migration (PreSDSDM) algorithm was implemented with short-offset Marine Multichannel Seismic (MCS) data, the MGL0908-TST survey, and Long-offset active source land seismic from Transect T4b Line of TAIGER Project. Tomographic velocity models from Kuo-Chen et al. (2012), Huang et al. (2014), and Roma (2018) are also utilized in this study to derive 1D velocity for implementing PreSDSDM.
    For deep structure imaging with seismic data. The imaging principle follows generalized imaging conditions for 3D depth migration. Resolving power of proposed approach depends on propagating velocity and frequency bandwidth of waveform, geometry of source-receiver pair and spatial aperture coverage of available data. Such kinds of details were discussed in this study.
    Current results shows outstanding output for MGL0908-TST survey. Several detailed structural were distinguishable and supports several interpretation from other study. On the other hand, the results from active source land seismic survey from TAIGER project surfers from insufficient source and receiver coverage, crooked line, multiple diving wave (Huang et al., 2013), and low signal to noise ration. However, we’re still able to distinguish features that indicates certain structures, such as basement under Coastal Plain and Chukou Fault which is on the eastern side of shot point 2. In addition, apparent migration smiles could be found on the image migrated with different tomographic velocity model. Such fact implies that, current tomographic velocity models appear to overestimate velocity at the shallow part of subsurface.
    Further improvement could be made in the future for better imaging results, such as adopting full 3D depth migration, improving steps for data conditioning, and to uses passive source for imaging.

    摘要 i Abstract iii Acknowledgement 致謝 iv Table of Contents v List of Figures vii List of Tables xvi Chapter 1 Introduction 1 1.1 Earth’s Subsurface Structure Imaging 1 1.2 Overview of TAIGER Project and Used Dataset 2 1.3 Tectonic Setting and Geological Background of Research Area 3 1.4 Objectives of Study 7 1.5 Arrangement of Thesis 7 References of Chapter 1 18 Chapter 2 Principle of Pre-stack Diffraction Stack Depth Migration 20 2.1 Principle of Depth Imaging 20 2.1.1 Principles of Diffraction Stack and Pre-stack Depth Migration 21 2.2 Migration Kernel and Imaging Condition 27 2.2.1 Finite-difference Travel-time Calculations 28 2.2.2 Stretching Effect of Fresnel Zone along Offset 29 2.3 Migration of Synthetic Acoustic and Elastic Data 32 2.3.1 2D Synthetic Seismic Profile along CPC-84-PSN-VA Seismic Line 34 2.3.2 Synthetic Application for Locating and Imaging Earthquake Source 36 References of Chapter 2 56 Chapter 3 Long-offset Seismic Line MGL0908-TSTGL0908-TST (TAIGER 2009 Marine - Leg 4, R/V Langseth) 61 3.1 Depth Image Reconstruction with Practical PreSDSDM Approach 61 3.2 Data Description 67 3.3 Methods and Results 68 3.3.1 Data Pre-Processing 68 3.3.2 VRMS, VInterval and VStacking Distribution Determined From Semblance Analysis 70 3.3.3 Offset Dependent Imaging Condition and Migration Strategy 70 3.4 Results, Discussions and Interpretation 76 References of Chapter 3 96 Chapter 4 TAIGER 2008 Land Experiment: South Main Array (Transect T4b Line) 99 4.1 The TAIGER (Taiwan Integrated Geodynamic Research) Project 99 4.2 Data Description 100 4.2.1 Survey Location and Geometry 100 4.2.2 Spectral Analysis of Four Common Shot Gathers 101 4.2.3 Multiple diving waves 104 4.3 Pre-Processing Strategy 105 4.3.1 Universal Filtering 106 4.3.2 Deconvolution (DECON) 107 4.3.3 Adaptive Filtering 107 4.3.4 Time-Frequency Domain Noise Suppression 107 4.4 Available Reference Velocity Models and Evaluations 108 4.4.1 Migration Velocity Model Setup 110 4.4.2 Offset-Dependent 1D Velocity Models from Roma (2018) 110 4.5 Migration Results and Interpretations 111 References of Chapter 4 148 Chapter 5 Discussions and Conclusions 149 5.1 Over-all Scopes on Data Pre-conditioning, Model Preparation and Migration Features 149 5.2 Discussions and Conclusions 153 5.3 Limitation of this Study and Suggestions for Future Work 158 References 160

    Abbott, R. E. and J.N. Louie, 2001, Geophysical confirmation of low-angle normal slip on the historically active Dixie Valley fault, Nevada
    Abma, R., J. Sun, N. Bernitsas, 1999, Antialiasing methods in Kirchhoff migration, Geophysics, Vol. 64, No. 6, 1783-1792.
    Angelier, J., Lee, J.C., Hu, J.C. and Chu, H.T. 2003. Three-dimensional deformation along the rupture trace of the September 21st, 1999, Taiwan earthquake: a case study in the Kuangfu school. J. Struct. Geol.
    Bancroft, J.C., 1997, A Practical Understanding of Pre/Poststack Vol.1.SEG.
    Bancroft, J.C., 2007, A Practical Understanding of Pre/Poststack Vol.2 (Pre stack). SEG.
    Basler-Reeder, K., J. Louie, S. Pullammanappallil, G. Kent, 2016, Joint optimization of vertical component gravity and P-wave first arrivals by simulated annealing, Geophysics, Vol. 81, No. 4, 10.1190/GEO2015-0643.1, ID59-ID71.
    Bauer K., T. Ryberg, G. S. Fuis, and S. Lüth, 2013, Seismic Imaging of the Waltham Canyon Fault, California Comparison of Ray-Theoretical and Fresnel Volume Prestack Depth Migration, Bull. Seismol. Soc. Am. 78, 2062–2076.
    Baysal, E., Kosloff, D.D. and Sherwood, J.W.C., 1983, Reverse time migration. Geophysics, 48, 1514-1524.
    Bednar, J.B., Yoon, K., Shin, C. and Lines, L.R., 2003, One Way vs Two Way Wave Equation Imaging – Is Two-Way Worth It? 65th EAGE Conference & Exhibition, Extended Abstracts, B11.
    Berkhout, A. J., 1984, Seismic resolution—A quantitative analysis of resolving power of acoustical echo techniques: Geophysical Press.
    Bernitsas, N., Sun, J. and Sicking, C., 1997, Prism waves – an explanation for curved seismic horizons below the edge of salt bodies. 59th EAGE Conference & Exhibition, Extended Abstracts.
    Bevc, D. and D. E. Lumley, 1994, When is anti-aliasing needed in Kirchhoff migration?, Stanford exploration project, Report 80, 467-476.
    Bevc, D. and Jon, Claerbout, 1997, Fast anti-aliased Kirchhoff migration and modeling, Short Note, Stanford exploration project, Report 75, Nov. 17, 97-103.
    Beylkin, G., 1985, Imaging of discontinuities in the inverse scattering problem by inversion of a causal generalized Radon transform: J. Math. Phys. 26, 99–108.
    Beylkin, G., Oristaglio, M., and Miller, D., 1985, Spatial resolution of migration algorithms, in Berkhout, A. J., Ridder, J., and van der Waals, L. F., Eds., Proc. 14th Internat. Symp. on Acoust. Imag.,155–167.
    Bleistein, N., 1987, On the imaging of reflectors in the earth: Geophysics, 52, 931–942.
    Cavalca, M., and Lailly, P., 2005, Prismatic reflections for the delineation of salt bodies. 72nd Annual International Meeting, SEG, Expanded Abstracts.
    Chávez-Pérez, S. and J. Louie, 1998, Crustal imaging in southern California using earthquake sequences, Tectonophysics, 286, 223-236.
    Chiang, C.S., Yu H.S. and Chou, Y.W. 2004. Characterisics of wedge-top depozone of the southern Taiwan foreland basin. Basin Research, 16: 65-78.
    Chi, W.-C., Reed, D.L., Liu, C.-S., Lundberg, N., 1998. Distribution of the bottom stimulating reflector in the offshore Taiwan collision zone. Terr. Atmos. Ocean. Sci. 9(4), 779-794.
    Claerbout, J.F., 1971, Toward a unifned theory of reflector mapping. Geophysics, 36, no. 3, 467-481.
    Claerbout, J. E., 1997, Anti aliasing, Stanford exploration project, Report 73, Nov. 18, 395-415.
    Claerbout, J. F., 2010, Basic Earth Imaging, online free book, Stanford Exploration Project.
    Cohen, J. K., Hagin, F. G., and Bleistein, N., 1986, Three-dimensional Born inversion with an arbitrary reference: Geophysics, 51, 1552–1558.
    Eaton, D. W., 2010, Microseismic Moment Tensors the Good, the Bad and the Ugly, CSEG Recorder, Expanded Abstracts, Nov. 2010, 44-47.
    Effrianto, P.A., 2018, Master Thesis: AVO/AVA Modeling, Analysis and Inversion of CMP Gathers for Gas Hydrates Investigation in Lower Fangliao Basin, SW-Taiwan, National Central University.
    Eyre, T. S., 2015, Overview of moment-tensor inversion of microseismic events, The Leadding Edge, Aug. 2015,
    Francis T. Wua, H. Kuo-Chen, K.D. McIntosh, 2014, Subsurface imaging, TAIGER experiments and tectonic models of Taiwan. Journal of Asian Earth Sciences 90 (2014) 173–208
    Gray, S.H., 1992, Frequency-selective design of the Kirchhoff migration operator. Geophysical Prospecting, 40, 565-572.
    Gray, S., 2004, Gaussian beam migration of common-shot records. 71st Annual International Meeting, SEG, Expanded Abstracts, 953-956.
    Hagedoorn, J. G., 1954, A Process of Seismic Reflection Interpretation, Geophysical Prospecting, II, 85-127.
    Hale, D., 1991, Stable explicit depth extrapolation of seismic wavefields. Geophysics, 56, 1770-1777.
    Hale, D., Hill, N.R. and Stefani, J., 1992, Imaging salt with turning seismic waves. Geophysics, 57, no.11, 1453-1462.
    Hemon, C., 1978, Equations d’onde et modeles. Geophysical Prospecting, 26, 790–821.
    Hill, N.R., 1990, Gaussian beam migration. Geophysics, 55 (11), 1416–1428.
    Ho, C. S., 1975, An Introduction to the Geology of Taiwan: Explanatory Text on the Geologic Map of Taiwan, 1975, The Ministry of Economic Affairs, Republic of China, Taipei, Taiwan, 277 pages with 16 plates.
    Ho, C. S., 1979, Geologic and tectonic framework of Taiwan. Geol. Soc. China Mem. 3, 57-72.
    Huang, B.S., Wang, C.Y., Okaya, D., Lee, S.J., Lai, Y.C., Wu, F., Liang, W.T., Huang, W.G., 2013. Multiple diving waves and steep velocity gradients in the western Taiwan coastal plain: an investigation based on the TAIGER experiment, Bulletin of the Seismological Society of America, Vol. 103, No. 2A, pp. 925–935.
    Huang, T. C., 1980, Oligocene to Pleistocene calcareous .nannofossil biostratigraphy of the Hsuehshan Range and western foothills in Taiwan: Geol. Paleont, SE Asia, vol. 21, pp. 191-210.
    Huang, T. Y., 1978, Significant new look on the Tertiary stratigraphy of Taiwan: Petroleum Geology of Taiwan, no.15, p. 167-179.
    Huang, Y. and G. T. Schuster, 2014, Resolution limits for wave equation imaging, J. Applied Geophysics, 107, 137-148.
    Ian F. Jones, Robert I. Bloor, Biondo L. Biondi, and John T. Etgen, 2008, Prestack Depth Migration and Velocity Model Building. Society of Exploration Geophysicists
    Jones, I.F. and Lambaré, G., 2003, Wave equation versus ray based imaging. First Break, 21 (2), 11–13.
    Koren, Z., Ravve, I., Bartana, A. and Kosloff, D., 2007, Local angle domain in seismic imaging. 69th EAGE Conference & Exhibition, Extended Abstracts, P287.
    Kuo-Chen, H., Wu, F.T., Roecker, S.W., 2012. Three-dimensional P velocity structures of the lithosphere beneath Taiwan from the analysis of TAIGER and related seismic data sets. Journal Of Geophysical Research, Vol. 117, B06306, doi:10.1029/2011JB009108.
    Kvasnička, M. and V. Červený, 1996, Analytical expression for fresnel volumes and interface fresnel zones of seismic body waves. Part 1 Direct and unconverted reflected waves, Studia Geoph. Et Geod., 40, 136-155.
    Kvasnička, M. and V. Červený, 1996, Analytical expressions for Fresnel volumes and interface fresnel zones of seismic body waves. Part 2 Transmitted and converted waves. Head waves, Studia Geoph. Et Geod., 40, 381-397.
    Kwiatek, G., B.Z. Yehuda, 2013, Assessment of P and S wave energy radiated from very small shear-tensile seismic events in a deep south African mine, J. Geophys. Res., 118, 3630-3641, doi:10.1002/jgrb.50274.
    Lin, A. T., Watts, A. B., and Hesselbo, S. P., 2003. Cenozoic stratigraphy and subsidence history of the South China Sea margin in the Taiwan region. Basin Research (2003) 15, 453–478, doi:10.1046/j.1365-2117.2003.00215.x
    Lin, C.-C., 2005, Master Thesis: Structure and Sedimentary Characteristics of Gas Hydrate-Bearing Sediments off Southwest Taiwan, National Central University.
    Lin, C.-C., Lin, A.T., Liu, C.-S., Schnurle, P., Chen, G.-Y., Liao, W.-Z., 2008, Geological controls on BSR occurrences in the incipient arc-continent collision zone off southwest Taiwan, Marine and Petroleum Geology (2008), doi:10.1016/j.marpetgeo.2008.11.002.
    Lin, C.C., 2010, Geological controls of gas hydrate occurances and gas hydrate resource assessment, offshore southwest Taiwan. National Cental University, PhD thesis.
    Liu, C.-S., Schnurle, P., Wang, Y., Chuang, S.-H., Chen, S.-C., Hsiuan, T.-H. 2006. Distribution and characters of gas hydrate offshore of southwestern Taiwan. Terr. Atmos. Ocean. Sci. 17, 615-644.
    Lin, Y. P., Li Zhao, and Hung, S. H., 2011. Assessment of Tomography Models of Taiwan Using First-Arrival Times from the TAIGER Active-Source Experiment
    Leveille, J.P., Jones, I.F., Zhou, Z-Z., Wang, B., Liu, F., 2011, Subsalt Imaging for Exploration, Production and Development. A Review. Geophysics, 76, WB3–WB20.
    Liu, F., Zhang, G., Morton, S.A. and Leveille, J.P., 2011, An effective imaging condition for reverse-time migration using wavefield decomposition. Geophysics, 76, S29–S39.
    Louie, J. N., R. W. Clayton, R. J. LeBras, 1988, Three-dimensional imaging of steeply dipping structure near the San Andreas fault, Parkfield, California, Geophysics, Vo 53, No. 2, 176-185.
    Louie, J. N. and J. Qin, 1991, Subsurface Imaging of the Garlock Fault, Cantil Valley, California, 1979-1983, J. Geophys. Res., 96, No B9, 14461-14479, Aug. 10.
    Louie, J. N. et al., 2002, Multimode migration of scattered and converted waves for the structure of the Hikurangi slab interface, New Zealand, Tectonophysics, 355, 227-246.
    Lumley, D. E., J. F. Claerbout, D. Bevc, 1994, Anti-aliasing Kirchhoff 3D migration, Stanford exploration project, Report 80, May. 15, 447-490.
    Lumley, D. E., J. F. Claerbout, D. Bevc, 1994, Anti-aliasing Kirchhoff 3D migration, Society of Exploration Geophysics, Expanded abstract, SM4.6.
    Lüth, S. Buske,R. Giese and A. Goertz, 2005, Fresnel volume migration of multicomponent data, Geophysics, Vol 70, No. 6, Nov-Dec, S121-S129.
    Malavieille, J., and Trullenque, G., 2007. Consequences of continental subduction on forearc basin and accretionary wedge deformation in SE Taiwan: Insights from analogue modeling. Tectonophysics 466 (2009) 377–394.
    McMechan, G.A., 1983, Migration by extrapolation of time-dependent boundary values. Geophysical Prospecting, 31, 413–420.
    Nowack, R. L., 1992, Wavefront and solutions of the eikonal equation, Geophys. J. Int., 110, 55-62,
    Ou, G. B., 2008, Seismological Studies for Tensile Faults, Terr. Atmos. Ocean. Sci., Vol. 19, No. 5, 463-471.
    Popov, M.M., 1982, A new method of computation of wave fields using Gaussian beams. WaveMotion, 4, 85–97.
    Pisen, S., 2014, Master Thesis: Reservoir Characterization via Amplitude versus Offset Analysis and Impedance Inversion in the Thrace Basin of Northwest Turkey, University of Houston
    Reed, D.L., Lundberg, N., Liu, C.-S., Kuo, B.-Y. 1992. Structural relations along the margins of the offshore Taiwan accretionary wedge: implication for accretion and crustal kinematics. Acta Geol. Taiwan 30, 105-122
    Roecker, S., C. Thurber, K. Roberts, L. Powell, 2006, Refining the image of the San Andreas Fault near Parkfield, California using a finite difference travel time computation technique, Tectonophysics, 426, 189-205.
    Roma Widiyansari., 2018, Master Thesis: Body- and Surface-waves Wavefield Inversion in Southern Taiwan: Application to TAIGER T4b Wide Angle Refraction/Reflection Data, National Central University.

    Sava, P. and Fomel, S., 2006, Time-shift imaging condition in seismic migration. Geophysics, 71, S209–S218.
    Schuster, G. T., Decomposition of migration kernel into upgoing and downdoging components, Society of Exploration Geophysics, Expanded abstract, SM4.7.
    Schuster, G. T., 2010, Basics of Seismic Imaging, Cambridge University Press. Copyright G. Schuster 2010.
    Sergio Chávez-Pérez and John N.Louie, 1998, Crustal imaging in southern California using earthquake sequences. Tectonophysics 286 (1998) 223-236
    Shui-Beih Yu, Horng-Yue Chen, Long-Chen Kuo, 1997, Velocity field of GPS stations in the Taiwan area. Tectonophysics 274 (1997) 41-59
    Suppe, J., 1981. Mechanics of mountain-building and metamorphism in Taiwan, Mem.Geol. Soc. China 4, 67–89
    Teng, L.S., 1990. Geotectonic evolution of late Cenozoic arc-continent collision in Taiwan. In: J. Angelier (Editor), Geodynamic Evolution of the Eastern Eurasian Margin. Tectonophysics, 183: 57-76.
    Tian, N., T.E. Fan, Z.J. Wang, W. T. Cai, Lateral wave field stacking of seismic Fresnel zones for the generalized-offset case, Applied Geophysics, Vol. 12, No. 2, June, 2015, DOI:10.1007/s11770-015-0488, 236-243.
    Vaillant, L. and S. Fomel, 1999, On Stolt stretch time migration, Stanford exploration project, Report 102, Oct. 25, 1-11.
    Vidale, J. , 1988, Finite-difference calculation of travel times, Bull. Seismol. Soc. Am., 78, 6, 2062–2076.
    Vidale, J., 1990, Finite-difference calculation of traveltimes in three dimensions, Geophysics, Vol 55, No. 5, 521-526.
    Welayaturromadhona, 2016, Master Thesis: AVO Analysis of Detecting Submarine Gas Hydrate in Lower Fangliao Basin, SW-Taiwan, National Central University.
    Williamson, P., Wang, B., Bevc, D., Jones, I.F., 2010, Full wave-equation methods for complex imaging challenges. The Leading Edge, 29 (3), 264–268
    Whitmore, N.D., 1983, Iterative depth migration by backward time propagation. 53rd Annual International Meeting, SEG, Expanded Abstracts.
    Woodward, M. J. , 1992, Wave-equation tomography, Geophysics, Vol. 57, No 1, Jan., 1992, 15-26.
    Wu, Q.J., Y. H. Li, R.Q. Zhang, N. L. Zhang, 2007, 2D Kirchhoff migration for receiver function, Vol. 50, No. 2, 476-482.
    Wu, R. and K. Aki, 1985, Scattering characteristics of elastic waves by an elastic heterogeneity, Geophysics, Vol. 50, No. 4, 582-595.
    Yilmaz. O., 2001. Seismic Data Analysis, Society of Exploration Geophysicists
    Yoon, K. and Marfurt, K.J., 2006, Reverse-time migration using the Poynting vector. Exploration Geophysics, 37, 102–107.
    Yu, C. Y., Yang, C. H., Yang, K., Cheng, R.C., & Huang, S.T., 2009, A Study of Complex Substructures in a Foothill Belt of Northwestern Taiwan Using Two-Dimensional Gravity Multiple-Source Model Inversion, Terr. Atmos. Ocean. Sci., 20, 791-798, doi: 10.3319/TAO.2009.02.02.01(T)
    Zahradník, J., 2015, Analysis of the source scanning algorithm with a new P-wave picker, J. Seismol., 19, 423-441.
    Zhu, J., R Lines and H Gray, 1997, Smiles and frowns in migration velocity analysis, CREWES Research Report, Vol. 9, 33, 1-21.
    Zhu, T., 2012, Complex-beam migration and land depth imaging, CSEG Recorder, Expanded Abstracts, Dec. 2012, 22-26.

    QR CODE
    :::