摘要
讨论了全局选择震源初始位置下的矩阵反演求取全局解的问题.与流行的全局优化搜寻直接定位方法相比,该反演算法是一种行之有效的定位方法,具有定位精度高、计算时间少以及对噪声数据不敏感等特点.其突出的优点是在不增加计算难度和计算时间的前提下,用矩阵反演的方法确保得到全局最小值解,可适应于地震早期预警,海啸早期预警,以及大震速报等实际工作.
A new technique of matrix inversion with initial global sampling was developed to seek a global solution.Comparing with the global grid search method it is a practical selection for earthquake location procedure in the sense of solution accuracy,computational efficiency and non-sensitivity to noise data.The dominant feature of the new method is that a global solution can be found with matrix inversion without extra computational time and mathematic complexity required.It is a good choice to the earthquake early warning,tsunami early warning and rapid hazard assessment and emergency response after strong earthquake occurrence.
引文
刘福田.1984.震源位置和速度结构的联合反演(Ⅰ):理论和方法[J].地球物理学报,27(2):167--175.
王椿镛,王溪莉,颜其中.1993.昆明地震台网多事件定位问题的初步研究[J].地震学报,15(2):136--145.
杨智娴,陈运泰.2004.用双差地震定位法再次精确测定1998年张北—尚义地震序列的震源参数[J].地震学报,26(2):115--120.
周仕勇,许忠淮,韩京.1999.主地震定位方法分析以及1997年新疆伽师震群高精度定位[J].地震学报,21(3):258--265.
Bai C Y,Greenhalgh S A,Zhou B.2007.3-D ray tracing by a modified shortest path method[J].Geophysics,97:T27--T36.
Bai C Y,Greenhalgh S A.2006.3-Dlocal earthquake hypocenter determination with an‘irregular’shortest-path method[J].Bull Seism Soc Amer,96:2257--2268.
Bai Chao Ying.2004.Three-dimensional seismic kinematic inversion with application to reconstruction of the velocitystructure of Rabaul volcano[D].Adelaide:Adelaide University:127--146.
Crosson R S.1976.Crustal structure modeling of earthquake data,Simultaneous least squares estimation of hypocenterand velocity parameters[J].J Geophys Res,81:3036--3046.
Douglas A.1976.Joint epicenter determination[J].Nature,215:45--48.
Julian B R,Gubbins D.1977.Three-dimensional seismic ray tracing[J].J Geophys,43:95--113.
Lomax A,Curtis A.2001.Fast,probabilistic earthquake location in3D models using oct-tree importance sampling[J/OL].Geophys Res Abstr,3:955(www.alomax.net/nlloc/octtree).
Moser TJ,Van Eck T,Nolet G.1992.Hypocenter determinationin strongly heterogeneous earth models usingthe shor-test path method[J].J Geophys Res,97:6563--6572.
Nelson G D,Vidale J E.1990.Earthquake location by3-Dfinite difference travel times[J].Bull Seism Soc Amer,80:395--410.
Sambridge MS,Kennett B L N.2001.Seismic event location:Nonlinear inversion using a neighbourhood algorithm[J].Pure Appl Geophys,158:241--257.
Spencer C,Gubbins D.1980.Travel-time inversion for simultaneous earthquake location and velocity structure determi-nationin laterally varying media[J].Geophys J R astr Soc,63:95--116.
UmJ,Thurber C.1987.Afast algorithmfor two-point seismic ray tracing[J].Bull Seism Soc Amer,77:972--986.
Waldhauser F,Ellsworth WL.2000.Adouble-difference earthquake location algorithm:Method and application to theNorthern Hayward Fault,California[J].Bull Seism Soc Amer,90:1353--1368.
Yong Z X,Waldhauser F,Chen Y T,Richards P G.2005.Double-difference relation of earthquakes in central-westernChina[J].J Seismology,9:241--264.