三维复杂层状介质中多震相走时联合反演成像
|
摘要
采用新近提出的改进型不规则最短路径多次波射线追踪正演技术,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,讨论了三维复杂层状模型中利用多震相走时资料进行联合反演成像的技术方法.考虑到不同震相种类走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数;另外,考虑到同时反演速度模型和反射界面起伏中不同参数变化对走时影响程度的不同,Jacobi偏导矩阵元素中引入了不同种类参数的归一化因子.几种数值模拟实例表明:多震相走时的联合或同时反演成像是一种提高走时成像空间分辨率,进而降低重建速度模型失真度行之有效的方法技术.
By exploring a newly developed multiple(reflected,transmitted and converted)arrival tracking algorithm with modified irregular shortest-path method,and combining the conjugated gradient method to solve the constrained,damped least square problem,we discuss and analyze a simultaneous inversion algorithm combining different travel time datasets in 3D complex layered media.In the inversion process we introduce different weighting factors according to the different picking errors of different sets of arrival times.In order to balance the influence on the two different model parameters in simultaneous inversion,we normalize two different model parameters in Jacobi matrix in the inversion process.The numerical results show that it is a practical and efficient way to improve the spatial resolution and reduce the artifact in reconstructing velocity distribution and reflected subsurface interface.
By exploring a newly developed multiple(reflected,transmitted and converted)arrival tracking algorithm with modified irregular shortest-path method,and combining the conjugated gradient method to solve the constrained,damped least square problem,we discuss and analyze a simultaneous inversion algorithm combining different travel time datasets in 3D complex layered media.In the inversion process we introduce different weighting factors according to the different picking errors of different sets of arrival times.In order to balance the influence on the two different model parameters in simultaneous inversion,we normalize two different model parameters in Jacobi matrix in the inversion process.The numerical results show that it is a practical and efficient way to improve the spatial resolution and reduce the artifact in reconstructing velocity distribution and reflected subsurface interface.
引文
[1] Aki K,Lee W H K.Determination of three dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes.J Geophys Res,1976,81:4381～4399
[2] Zhao D.Global tomographic images of mantle plumes and subducting slabs,insight into deep Earth dynamics.Phys Earth Planet Inter,2004,146:3～34
[3] Lei J,Zhao D.A new insight into the Hawaiian plume.Earth Planet Sci Lett,2006,241 :438～453
[4] Lei J,Zhao D.Global P-wave tomography;on the effect of various mantle and core phases.Phys Earth Planet Inter,2006,154:44～69
[5] Zhao D,Hasegawa A,Horiuchi S.Tomographic imaging of P and S wave velocity structure beneath northeastern Japan.J Geophys Res ,1992,97:19909～19928
[6] Lei J,Xie F,Lan C,et al.Seismic images under the Beijing region inferred from P and PmP data.Phys Earth Planet Inter,2008,168:134～146
[7] Zhao D,Todo S,Lei J.Local earthquake reflection tomography of the Landers aftershock area.Earth Planet Sci Lett,2005,235:623～631
[8] 刘福田.震源位置和速度结构的联合反演(Ⅰ)--理论和方法.地球物理学报,1984,27(2) :167～175 Liu F T.Simultaneous inversion of earthquake hypocenters and velocity structure (Ⅰ)-Theory and method.Chinese J Geophys (in Chinese) ,1984,27(2) :167～175
[9] 华标龙,刘福田.界面和速度反射联合成像--理论与方法.地球物理学报,1 995,38(6) :750～756 Hua B L,Liu F T.Joint reflection imaging of velocity and interface-Theory and method.Chinese J Geophys (in Chinese),1995,38(6) :750-756
[10] 周龙泉,刘福田,陈晓非.三维介质中速度结构和界面的联合成像.地球物理学报,2006,49(4) :1062～1067 Zhou L Q,Liu F T,Chen X F.Simultaneous tomography of 3-D velocity structure and interface.Chinese J Geophys (in Chinese),2006,49(4) :1062～1067
[11] 李松林,吴宁远,宋占隆等.速度分布和界面位置的联合反演.地震学报,1997,19(4) :383～392 Li S L,Wu N Y,Song Z L,et al.Simultaneous inversion of velocity distribution and interface positions.Acta Seismologica Sinica (in Chinese),1997,19(4) :383～392
[12] Pavlis G L,Booker J R.The mixed discrete-continuous inversion problem:application to the simultaneous determination of earthquake hypocenter and velocity structure.J Geophys Res,1980,85:4801～4810
[13] Kennett B L N,Sambridge M S,Williamson P R.Subspace methods for large inverse problem with multiple parameter classes.Geophysical Journal,1988,94:237～247
[14] Williamson P R.Tomographic inversion in reflection seismology.Geophys J Int,1990,100:255～274
[15] Zelt C A,Hojka A M,Flueh E R,et al.3D simultaneous seismic refraction and reflection tomography of wide-angle data from the central Chilean margin.Geophys Res Lett,1999,26:2577～2580
[16] Zelt C A,Sain K,Naumenko J V,et al.Assessment of crustal velocity models using seismic refraction and reflection tomography.Geophys J Int,2003,153:609～626
[17] Zelt C A,Ellis R M,Zelt B C.Three-dimensional structure across the Tintina strike-slip fault,northern Canadian Cordillera,from seismic refraction and reflection tomography.Geophys J Int,2006,167:1292～1308
[18] McCaughey M,Singh S C.Simultaneous velocity and interface tomography of normal-incidence and wide-aperture seismic traveltime data.Geophys J Int,1997,131:87～99
[19] James W,Hobro D,Satish C,et al.Three-dimensional tomographic inversion of combined reflection and refraction seismic traveltime data.Geophys J Int,2002,152:79～93
[20] Zhang J,ten Brink U,Toksoz M.Nonlinear refraction and reflection travel time tomography.J Geophys Res ,1998,103 (B12) :29743～29757
[21] Korenaga J,Holbrook W,Kent G,et al.Crustal structure of the southeast Greenland margin from joint refraction and reflection seismic tomography.J Geophys Res,2000,105 (B9) :21591～21614
[22] 毛伟健,Stua G W.多震相的层析成像与在介质中的多波型射线跟踪.CT理论与应用研究,1996,5(4) :42～48 Mao W J,Stua G W.On multi-phase tomography and wavetype ray tracing in earth-media.CT Theory and Application (in Chinese),1996,5(4) :42～48
[23] Rawlinson N,Sambridge M.Multiple reflection and transmission phases in complex layered media using a multistage fast marching method.Geophysics,2004,69:1338～1350
[24] De Kool M,Rawlinson N,Sambridge M.A practical grid-based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous media.Geophys J Int,2006,167:253～270
[25] Bai C Y,Tang X P,Zhao R.2D/3D multiply transmitted,converted and reflected arrivals in complex layered media with the modified shortest path method.Geophys J Int,2009,179:201～214
[26] 唐小平,白超英.最短路径算法下三维层状介质中多次波追踪.地球物理学报,2009,52(10) :2635～2643 Tang X P,Bai C Y.Multiple ray tracing within 3-D layered media with the shorted path method.Chinese J Geophys (in Chinese),2009,52(10) :2635～2643
[27] 唐小平,白超英.利用改进后的最短路径法追踪二维复杂层 状介质中的多次波.地球物理学进展,2009,24(6) :2087～2096 Tang X P,Bai C Y.Multiple ray tracing within 2D layered media with the shorted path method.Progress in Geophys (in Chinese),2009,24(6) :2087～2096
[28] 赵瑞,白超英.复杂层状模型中多次波快速追踪--一种基于非规则网格的最短路径算法.地震学报,2010,32(4) :433～444 Zhao R,Bai C Y.Fast multiple ray tracing within complex layered media:the shortest path method based on irregular grid cells.Acta Seismologica Sinica (in Chinese),2010,32(4) :433～444
[29] Bai C Y,Huang G J,Zhao R.2D/3D irregular shortest-path raytracing for multiple arrivals and its applications.Geophys J Int,2010,183:1596～1612
[30] Menke W.Geophysical Data Analysis:Discrete Inverse Theory.San Diego,California:Academic Press Inc.,1984． 40～126
[31] Meyerholtz K A,Szpakowski S A.Convolutional quelling in seismic tomography.Geophysics,1989,54:570～580
[32] Shaw P R,Orcutt J A.Waveform inversion of seismic refraction data and applications to young Pacific crust.Geophys J R Astr Soc,1985,82:375～414
[33] Tarantola A,Valette B.Generalized non-linear inverse problems solved using the least squares criterion.Rev Geophys Space Phys,1982,20:219～232
[34] Zhou B,Sinadinovski C,Greenhalgh S A.Iterative algorithm for the damped minimum norm,least-squares and constrained problem in seismic tomography.Explor Geophys,1992,23:497～505
[35] Bai C Y,Greenhalgh S A.3-D non-linear travel time tomography:imaging high contrast velocity anomalies.Pure Appl Geophys,2005,162:2029～2049
[36] Bishop T N,Bube K P,Cutler R T,et al.Tomographic determination of velocity and depth in laterally varying media.Geophysics,1985,50:903～923
[37] Bai C Y,Greenhalgh S.3-D local earthquake hypocenter determination with an 'irregular' shortest-path method.Bull Seism Soc Am,2006,96:2257～2268
[38] Scales J.Tomographic inversion via the conjugate gradient method.Geophysics,1987,52:179～185
[39] 黄国娇,白超英.二维复杂层状介质中地震多波走时联合反演成像.地球物理学报,2010,53(12) :2972～2981 Huang G J,Bai C Y.Simultaneous inversion with multiple traveltimes within 2-D complex layered media.Chinese J Geophys (in Chinese),2010,53(12) :2972～2981
[2] Zhao D.Global tomographic images of mantle plumes and subducting slabs,insight into deep Earth dynamics.Phys Earth Planet Inter,2004,146:3～34
[3] Lei J,Zhao D.A new insight into the Hawaiian plume.Earth Planet Sci Lett,2006,241 :438～453
[4] Lei J,Zhao D.Global P-wave tomography;on the effect of various mantle and core phases.Phys Earth Planet Inter,2006,154:44～69
[5] Zhao D,Hasegawa A,Horiuchi S.Tomographic imaging of P and S wave velocity structure beneath northeastern Japan.J Geophys Res ,1992,97:19909～19928
[6] Lei J,Xie F,Lan C,et al.Seismic images under the Beijing region inferred from P and PmP data.Phys Earth Planet Inter,2008,168:134～146
[7] Zhao D,Todo S,Lei J.Local earthquake reflection tomography of the Landers aftershock area.Earth Planet Sci Lett,2005,235:623～631
[8] 刘福田.震源位置和速度结构的联合反演(Ⅰ)--理论和方法.地球物理学报,1984,27(2) :167～175 Liu F T.Simultaneous inversion of earthquake hypocenters and velocity structure (Ⅰ)-Theory and method.Chinese J Geophys (in Chinese) ,1984,27(2) :167～175
[9] 华标龙,刘福田.界面和速度反射联合成像--理论与方法.地球物理学报,1 995,38(6) :750～756 Hua B L,Liu F T.Joint reflection imaging of velocity and interface-Theory and method.Chinese J Geophys (in Chinese),1995,38(6) :750-756
[10] 周龙泉,刘福田,陈晓非.三维介质中速度结构和界面的联合成像.地球物理学报,2006,49(4) :1062～1067 Zhou L Q,Liu F T,Chen X F.Simultaneous tomography of 3-D velocity structure and interface.Chinese J Geophys (in Chinese),2006,49(4) :1062～1067
[11] 李松林,吴宁远,宋占隆等.速度分布和界面位置的联合反演.地震学报,1997,19(4) :383～392 Li S L,Wu N Y,Song Z L,et al.Simultaneous inversion of velocity distribution and interface positions.Acta Seismologica Sinica (in Chinese),1997,19(4) :383～392
[12] Pavlis G L,Booker J R.The mixed discrete-continuous inversion problem:application to the simultaneous determination of earthquake hypocenter and velocity structure.J Geophys Res,1980,85:4801～4810
[13] Kennett B L N,Sambridge M S,Williamson P R.Subspace methods for large inverse problem with multiple parameter classes.Geophysical Journal,1988,94:237～247
[14] Williamson P R.Tomographic inversion in reflection seismology.Geophys J Int,1990,100:255～274
[15] Zelt C A,Hojka A M,Flueh E R,et al.3D simultaneous seismic refraction and reflection tomography of wide-angle data from the central Chilean margin.Geophys Res Lett,1999,26:2577～2580
[16] Zelt C A,Sain K,Naumenko J V,et al.Assessment of crustal velocity models using seismic refraction and reflection tomography.Geophys J Int,2003,153:609～626
[17] Zelt C A,Ellis R M,Zelt B C.Three-dimensional structure across the Tintina strike-slip fault,northern Canadian Cordillera,from seismic refraction and reflection tomography.Geophys J Int,2006,167:1292～1308
[18] McCaughey M,Singh S C.Simultaneous velocity and interface tomography of normal-incidence and wide-aperture seismic traveltime data.Geophys J Int,1997,131:87～99
[19] James W,Hobro D,Satish C,et al.Three-dimensional tomographic inversion of combined reflection and refraction seismic traveltime data.Geophys J Int,2002,152:79～93
[20] Zhang J,ten Brink U,Toksoz M.Nonlinear refraction and reflection travel time tomography.J Geophys Res ,1998,103 (B12) :29743～29757
[21] Korenaga J,Holbrook W,Kent G,et al.Crustal structure of the southeast Greenland margin from joint refraction and reflection seismic tomography.J Geophys Res,2000,105 (B9) :21591～21614
[22] 毛伟健,Stua G W.多震相的层析成像与在介质中的多波型射线跟踪.CT理论与应用研究,1996,5(4) :42～48 Mao W J,Stua G W.On multi-phase tomography and wavetype ray tracing in earth-media.CT Theory and Application (in Chinese),1996,5(4) :42～48
[23] Rawlinson N,Sambridge M.Multiple reflection and transmission phases in complex layered media using a multistage fast marching method.Geophysics,2004,69:1338～1350
[24] De Kool M,Rawlinson N,Sambridge M.A practical grid-based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous media.Geophys J Int,2006,167:253～270
[25] Bai C Y,Tang X P,Zhao R.2D/3D multiply transmitted,converted and reflected arrivals in complex layered media with the modified shortest path method.Geophys J Int,2009,179:201～214
[26] 唐小平,白超英.最短路径算法下三维层状介质中多次波追踪.地球物理学报,2009,52(10) :2635～2643 Tang X P,Bai C Y.Multiple ray tracing within 3-D layered media with the shorted path method.Chinese J Geophys (in Chinese),2009,52(10) :2635～2643
[27] 唐小平,白超英.利用改进后的最短路径法追踪二维复杂层 状介质中的多次波.地球物理学进展,2009,24(6) :2087～2096 Tang X P,Bai C Y.Multiple ray tracing within 2D layered media with the shorted path method.Progress in Geophys (in Chinese),2009,24(6) :2087～2096
[28] 赵瑞,白超英.复杂层状模型中多次波快速追踪--一种基于非规则网格的最短路径算法.地震学报,2010,32(4) :433～444 Zhao R,Bai C Y.Fast multiple ray tracing within complex layered media:the shortest path method based on irregular grid cells.Acta Seismologica Sinica (in Chinese),2010,32(4) :433～444
[29] Bai C Y,Huang G J,Zhao R.2D/3D irregular shortest-path raytracing for multiple arrivals and its applications.Geophys J Int,2010,183:1596～1612
[30] Menke W.Geophysical Data Analysis:Discrete Inverse Theory.San Diego,California:Academic Press Inc.,1984． 40～126
[31] Meyerholtz K A,Szpakowski S A.Convolutional quelling in seismic tomography.Geophysics,1989,54:570～580
[32] Shaw P R,Orcutt J A.Waveform inversion of seismic refraction data and applications to young Pacific crust.Geophys J R Astr Soc,1985,82:375～414
[33] Tarantola A,Valette B.Generalized non-linear inverse problems solved using the least squares criterion.Rev Geophys Space Phys,1982,20:219～232
[34] Zhou B,Sinadinovski C,Greenhalgh S A.Iterative algorithm for the damped minimum norm,least-squares and constrained problem in seismic tomography.Explor Geophys,1992,23:497～505
[35] Bai C Y,Greenhalgh S A.3-D non-linear travel time tomography:imaging high contrast velocity anomalies.Pure Appl Geophys,2005,162:2029～2049
[36] Bishop T N,Bube K P,Cutler R T,et al.Tomographic determination of velocity and depth in laterally varying media.Geophysics,1985,50:903～923
[37] Bai C Y,Greenhalgh S.3-D local earthquake hypocenter determination with an 'irregular' shortest-path method.Bull Seism Soc Am,2006,96:2257～2268
[38] Scales J.Tomographic inversion via the conjugate gradient method.Geophysics,1987,52:179～185
[39] 黄国娇,白超英.二维复杂层状介质中地震多波走时联合反演成像.地球物理学报,2010,53(12) :2972～2981 Huang G J,Bai C Y.Simultaneous inversion with multiple traveltimes within 2-D complex layered media.Chinese J Geophys (in Chinese),2010,53(12) :2972～2981
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心