用户名: 密码: 验证码:
基于空间离散烈度点椭圆模型的历史地震参数估计方法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
由于地震仪器记录的历史较短,不足以在短时间、小尺度范围内弄清楚地震活动性规律,因而对历史地震的研究非常重要。地震史料只有转换成用现代地震参数(如震级、震中位置等)的形式来表达,才有可能用于地震科学的各种研究。已出版的各版本历史地震目录里地震参数估算方法主观经验性较强,不确定性较大;现有新方法多以圆烈度衰减关系为基础,忽略了烈度分布的椭圆特征。本文通过建立椭圆烈度分布模型,联立椭圆数学方程,直接带入原始烈度点信息来估算历史地震基本参数,并以蒙特卡洛方法对参数估定的不确定性进行定量分析,以提高历史地震参数估计的科学性,减少主观判断带入的不确定性。
     本文首先分析了历史地震资料的特点,认为我国历史地震记录虽然跨度时间长,资料丰富,但受历史条件所限,在时、空、强上分布均不均匀,地震记录与烈度点都存在着不完整性,而且历史地震参数不确定性较大。
     针对历史地震资料的特点,我们基于宏观影响场的图像特征来建立地震烈度椭圆分布模型。引入空间统计分析的概念,从空间数据的角度给出了离散烈度点的定义,阐述了空间离散烈度点建模理论。以原始烈度点空间分布为基础用最小二乘法拟合各地震各烈度区椭圆烈度估计线,进而对烈度估计线统计回归得到椭圆烈度分布模型。依据所能获取的地震烈度信息,本文建立了适用于中国南北地震带、东部地区和新疆地区的椭圆烈度分布模型。该模型的建立充分利用了大量原始离散烈度点信息,通过数学手段得到烈度估计线,避免了等震线圈定过程中的主观经验因素。
     根据所建立的椭圆烈度分布模型,联立考虑中心点和方向性的椭圆数学方程,建立了地震参数估计方程,给出了参数求解的基本条件,并用既有仪器测定记录又有宏观考察数据的现代地震进行了内符验算和外推验算,证明了算法的可行性。采用蒙特卡洛方法,通过大量的抽样模拟,定量分析了所得参数的不确定性,并进一步讨论了本文所提算法的稳定性和适用性。
     本文建立了一套估算地震震中与震级的方法,旨在用此方法来估算只有烈度记载且烈度数据点相对较少的历史地震参数。以云南地区和山西地区为例,将此方法应用于有真实历史记载的地震进行了试算。历史地震算例表明,此方法对烈度数据点较少、等震线不易勾画的历史地震颇为有效。处理过程直接明了,一定程度上减少了主观不确定性,可重复和再现。同时证明了该方法对烈度信息空间分布散乱的稳定性。
     由于单烈度点地震占了历史强震总数的近1/3,本文单独对这类地震进行了讨论分析。研究表明,历史强震目录里单烈度点地震的震级可能存在普遍偏小的问题,尤其是在新疆、西藏和青海等地广人稀而又地震活动频繁的地方,可能会造成大量历史大震的漏记或者错记。
     最后对研究工作进行了总结,讨论了存在的问题,并给出了使用历史地震资料的若干建议与研究展望。
     本文建立了椭圆烈度分布模型和地震参数估计方程,直接利用全部烈度点的空间分布信息来估算地震震级和震中,为历史地震参数的确定提供了新的思路与方法,进一步提升了地震史料的可用性,提高了历史地震参数的可靠性,对于历史地震的重建和历史地震参数的校核具有重要意义。
Due to the short history of seismic instruments record, it is insufficient to know well of the law of seismic activity in a short time and within a scope of small-scale, thus it's very important to study on the historical earthquakes. Historical earthquake materials only converted into the form of modern seismic parameters (such as magnitude, epicenter, etc.), are possible for a variety of earthquake science research. The methods of earthquake parameters estimation in published historical earthquake catalogs partly are subjectively empirical with great uncertainty; Existing new methods are mostly based on circle intensity attenuation relationship, ignored the elliptical distribution characteristic of intensity. This paper established elliptical intensity distribution model, combined with ellipse mathematic equation, used the information of origial intensity points derectly to calculate historical earthquakes'basic parameters, and used the Monte Carlo method for quantitative analysis of the uncertainty of parameter estimation. These works are in order to improve the scientific and reduce the subjective uncertainty of the paramters'estimation of historical earthquakes.
     Firstly, this paper analyzed the characteristics of historical earthquake materials. Although Chinese historical earthquake records span a long time and are information rich, but due to the historical conditions, they are unevenly distributed in time, space and magnitude. The earthquake records and intensity data points are both not integrated, and the parameters of historical earthquakes have great uncertainty.
     According to the characteristics of historical earthquake materials, this paper established elliptical intensity distribution model based on the image feature of macro earthquake influence field. We introduced the concept of spatial statistical analysis, gave the definition of the discrete intensity points from the perspective of spatial data, and elaborated spatial modeling theory about discrete intensity point. Based on the spatial distribution of original intensity points, we used the least-square method to fit elliptical intensity estimation lines of different earthquakes of different intensity zones, then regressed the elliptical intensity distribution model. According to the seismic intensity information that can be obtained, this article established elliptical intensity distribution models that apply to the Chinese north-south seismic belt, the eastern region and the Xinjiang region. These models take advantage of the spatial distribution of large numbers of original discrete intensity points, use mathematical means to get intensity estimated lines to avoid subjective factors in the process of isoseismal plotting.
     Substituting the established elliptical intensity distribution model into the ellipse mathematic equation which considers the center point and the direction, we established the earthquake parameters estimation equation, and gave the basic conditions for solving, did internally checking and extrapolation checking on modern events with both instrument records and macro investigation records, proved the feasibility of the algorithm. Using the Monte Carlo method, we quantitatively analyzed the uncertainty of calculated parameters by simulating a large number of sampling and further discussed the stability and applicability of the algorithm.
     This paper has established a method to estimate earthquakes' epicenter and magnitude aimed on using this method to estimate the historical earthquake parameters that have only intensity information and the intensity data points relatively few. In Yunnan and Shanxi Province as an example, this method is applied to real historical earthquakes. Historical earthquake examples show that this method is quite effective for historical earthquakes with less intensity data points and difficulty to delineate isoseismal. The process of our method is straightforward. It can reduce the subjective uncertainty to some extent, and can be repeatable and reproduced. This work has also proved that the method is stability for the scattered spatial distribution of intensity information.
     As the historical earthquakes with a single intensity point account for nearly1/3of the total historical record, we discussed such earthquakes alone. Studies have shown that the magnitude of such earthquakes in historical earthquake catalog may generally smaller than the real, especially in Xinjiang, Tibet and Qinghai such frequent seismic activity but sparsely populated areas, it may cause a lot of large historical earthquakes omission or wrong noted.
     Finally, the research works are summarized, problems are discussed and a few recommendations are given on how to use historical earthquake materials.
     This paper has established elliptical intensity distribution model and seismic parameters estimation equation. This method directly uses the spatial distribution information of all intensity points to estimate seismic parameters, provides new idea for historical earthquake parameters determination, further enhances the availability of historical earthquakes and improves the reliability of the historical seismic parameters. For the reconstruction and parameters'checking of historical earthquakes, it has important significance.
引文
陈培善,刘家森.1975.用位错模型研究震级与烈度的关系.地球物理学报,18(3):183-195
    陈达生,刘汉兴.地震烈度椭圆衰减关系.华北地震科学,1989,7(8):31-42
    陈棋福,郑大林,车时.2002,中国震例(1992-1994).北京:地震出版社,1-489
    陈棋福,郑大林,刘桂萍等.2002,中国震例(1995-1996).北京:地震出版社,1-429
    陈棋福,郑大林,高荣胜.2003,中国震例(1997-1999).北京:地震出版社,1-468
    陈棋福,郑大林,车时等.2008,中国震例(2000-2002).北京:地震出版社,1-570
    陈运泰,刘瑞丰,2004,地震的震级,地震地磁观测与研究,25(6),1-12
    陈念英,王树华,1994,震中区地面建筑物震害的力学特性与地震震中判别,成都大学学报[J],13(2),36-43
    曹学锋,阎秀杰,时振梁,2004,关于疑难历史地震研究,中国地震,20(3):294-298
    崔鑫,2010,地震烈度衰减模型研究与地震震动图理论初探,中国科学技术大学硕士学位论文,4
    常银辉,2012,基于断层距的地震烈度衰减关系研究,中国地震局工程力学研究所
    董瑞树,冉洪流,任国强,1998,历史地震资料在核动力工程地震安全性评价中的应用,中国地震,14(3),32-37
    董瑞树,周庆,冉洪流,2007.历史地震考察中历史资料辨伪与查证.震灾防御技术,2(1):50-57
    董瑞树,周庆,陈晓利,等.2009.1631年湖南省常德地震的再考证.地震地质,31(1):162-173
    顾功叙(主编),1983,中国地震目录,北京:地震出版社,1-663
    国家地震局震害防御司,1995,中国历史强震目录(公元前23世纪-公元1911年),北京:地震出版社,1-471
    国家地震局地球物理研究所,1989,复旦大学中国历史地理研究所,中国历史地震图集(远古至元时期),北京:中国地图出版社,1-174
    国家地震局地球物理研究所,1989,复旦大学中国历史地理研究所.中国历史地震图集(明时期).北京:地图出版社,1-194
    国家地震局地球物理研究所,1990,复旦大学中国历史地理研究所.中国历史地震图集(清时期).北京:中国地图出版社,1-240
    国家地震局,1981,中国地震烈度区划工作报告,北京:地震出版社,11
    国家核安全局,1998,核安全导则汇编(上),北京:中国法制出版社,343-374
    龚双辉,2007,空间统计分析及在区域经济中的应用,华中科技大学硕士学位论文,1-4
    胡聿贤.1999,地震安全性评价技术教程.北京:地震出版社,97,174-176
    胡聿贤,2003,地震工程学,北京:地震出版社,53
    胡聿贤.2006,地震工程学.北京:地震出版社,50-51,69-70
    霍俊荣.1989,近场强地面运动衰减规律的研究[博士论文].黑龙江:国家地震局工程力学研究所
    何彬彬,陈翠华,方涛,郭达志,2006,空间数据关联规则挖掘的不确定性处理及度量,地理与地理信息科学,22(6),5-8
    皇甫岗,罗荣联,2008,运用历史地震资料重在考证,地震研究,31(4),304-307
    黄永林,高俊锁,2009,历史地震资料整理与利用的现状与发展——以江苏省及其邻近地区为例[J],灾害学,24(3),103-108
    黄芳,2005,空间统计学及其在空间模式分析中的应用,华中师范大学硕士学位论文
    郝敏,谢礼立,徐龙军,2005,关于地震烈度物理标准研究的若干思考,地震学报,,27(2),230-234
    蒋溥,1993,工程地震学概论,北京,地震出版社,21
    金畅,2006,蒙特卡洛方法中随机数发生器和随机抽样方法的研究,大连理工大学硕士学位论文
    贾素娟,1996,历史地震震中的甄别与修订及其构造意义,中国地震学会第六次学术大会论文摘要集,199
    李善邦,1960,中国地震门录.北京:科学出版社,1-790
    李善邦,1981,中国地震,北京:地震出版社,207
    刘昌森,1989,对某些历史地震参数的修订意见,中国历史地震研究文集1,北京:地震出版社,62-67
    刘昌森,1993,历史有感地震处理的原则与方法,东北地震研究,9(3),16-24
    刘百篪,郑文俊,郭华等,2002,活断层工作方法在中早期历史地震研究中的作用,中国地震,18(3),283-288
    刘湘南,黄方,王平,2008,GIS空间分析原理与方法(第二版),北京:科学出版社,185-186
    卢荣俭,宋雅桐,陈达生,1981,等震线与震级的统计关系,地震工程研究报告集(第四集),北京:地震出版社,14-34
    雷建成,高孟潭,俞言祥,2006,西南地区近代地震的震中烈度与有感半径的统计研究,震灾防御技术,1(2):137-145
    林国良,王健,2012,基于烈度点的中国历史地震资料数据系统试编制,地震学报,34(1),118-124
    毛定一,1984,地震史料在地震预报中的实用价值,华南地震,4(4):87-89
    齐书勤,卢振恒,吴戈,孟宪东,1991,历史地震学简论,中国历史地震研究文集2,北京:地震出版社,1-5
    齐书勤,1991,中国历史地震研究十年,中国历史地震研究文集2,北京:地震出版社,10-14
    柔洁,伊力亚尔·阿不力孜,陈建波,2006,关于地震烈度的讨论,中国地震,22(1),94-102
    任雪梅,2011,中国地震局地球物理研究所博士学位论文,20-36
    时振梁,李裕彻,1989,中国历史地震研究的回顾和展望,中国历史地震研究文集1,北京:地震出版社,1
    时振梁,何淑韵,曹学峰,1991,应用历史地震资料研究地震活动性,中国历史地震研究文集2,北京:地震出版社,86-91
    时振梁,曹学锋,1998,历史地震目录修订问题,中国地震,14(3),1-5
    沈建文,华宜平,1989,关于地震烈度衰减模型的系统偏差.地震学报,11(1),38-45
    桑永胜,2010,空间数据分析的神经计算方法,电子科技大学博士学位论文,1
    谭其骧,1982,中国历史地图集,第七册:元、明时期,上海:地图出版社,42-43
    谭毅培,陈棋福,2012,利用宏观烈度资料估计历史地震震源参数的研究概述,地球物理学进展,27(6):2384-2394
    孙寿成,朱书俊,历史地震学及其研究,1991,《中国历史地震研究文集2》,北京:地震出版社,6-9
    孙继浩,2011,川滇及邻区中强地震烈度衰减关系的适用性研究,中国地震局地震预测研究所硕士学位论文
    汪素云,1993,地震重新定位在核电厂厂址地震安全性评价中的重要性,地震地质,15(4),389-394
    汪素云,武宦英,1993,华北地区地震烈度的衰减关系,中国地震区划文集,北京:地震出版社,185-191
    汪素云,俞言祥,高阿甲等,2000,中国分区地震动衰减关系的确定.中国地震,16(2):99-106
    汪素云,张晓东,时振梁,2002,历史地震资料的应用与开发,中国地震,18(4),371-376
    王泽皋,1985,大震区的“长期活动”,地震学报,7(3),254-266
    王健.1999.地震活动性震级结构特征及其应用研究.北京:中国地震局地球物理研究所:61-62
    王健.2007.渤海海域历史地震和地震海啸.地震学报,29(5):549-557.
    王健,吴宣,高阿甲,陈鲲,邵磊,2010,历史地震参数校订方法探讨—以北京延庆公元294年地震为例,地震学报,32(1):108-117
    王光远,结构软设计理论初探,1987,哈尔滨:哈尔滨建工学院出版社,1-18
    王虎栓,1991,地震烈度及其物理标准的研究,世界地震工程,4,14-21
    王越,1998,北京历史地震资料汇编,北京:专利文献出版社,1-186
    吴清,高孟潭,徐伟进,2012,历史强震震中精度统计特征及其对地震危险性研究的影响,地震学报,34(4),537-548
    吴戈,刘昌森,瞿文杰,江在雄,吴镝等,2001,黄海及其沿岸历史地震编目与研究,北京:地震出版社,1-152
    武汉大学测绘学院测量平差学科组,2003,误差理论与测量平差基础,武汉:武汉大学出版社,6-22
    谢毓寿,蔡美彪(主编),1983,中国地震历史资料汇编(第1卷).北京:科学出版社,1-227
    谢毓寿,蔡美彪(主编),1985a,中国地震历史资料汇编(第2卷).北京:科学出版社,1-949
    谢毓寿,蔡美彪(主编),1985b,中国地震历史资料汇编(第4卷)(上).北京:科学出版社,1-729
    谢毓寿,蔡美彪(主编),1986,中国地震历史资料汇编(第4卷)(下).北京:科学出版社,1-258
    谢毓寿,蔡美彪(主编),1987a,中国地震历史资料汇编(第3卷)(上).北京:科学出版社,1-540
    谢毓寿,蔡美彪(主编),1987b,中国地震历史资料汇编(第3卷)(下).北京:科学出版社,541-1427
    谢毓寿,蔡美彪(主编),1988,中国地震历史资料汇编(第5卷).北京:科学出版社,1-782
    谢毓寿,1982,历史地震震级的评定,地震地磁观测与研究,3,10-13
    谢毓寿,1986,中国的历史地震研究,中国地震,2(2),8-15
    谢毓寿,地震烈度,1988,北京:地震出版社,81-82
    谢毓寿,1991,历史地震参数的估定,东北地震研究,7(1),1-6
    徐吉廷,1994,阿坝州历史地震的烈度评估定与震级估算,阿坝科技,(1),48-53
    徐爱信,2007,历史地震档案资料在防震减灾中的应用,国际地震动态,(2),43-46
    肖亮,2011,水平向基岩强地面运动参数衰减关系研究,中国地震局地球物理研究所
    许卫晓,2011,烈度分布快速评估方法研究,中国地震局工程力学研究所硕士学位论文,9-12
    鄢家全,张志中,王健,温增平,俞言祥,刘爱文,潘华,吕红山,郝玉芹,2011,中国历史地震烈度表研究,地震学报,33(4):515-531
    游扬声,2005,一般分布模式下GIS位置数据的不确定性研究,武汉大学博士学位论文
    张四昌,1989,宏观震中及其定位精度的探讨,中国历史地震研究文集1,北京:地震出版 社,127-130
    张扬,马干,史保平,张健,杨勇.2009.华北地区烈度衰减模型建立及其用于震中区域和震级的定量估算.地震学报,31(3):290-306.
    中国科学院地震工作委员会历史地震工作组,1956.中国地震资料年表,北京:科学出版社,1-1653
    中央地震工作小组,1971,中国地震目录,北京:科学出版社,1-663
    中国科学院地球物理研究所,1976,地震学基础,北京:科学出版社,42-45
    中国地震局震灾应急救援司,2010,2001-2005年中国大陆地震灾害损失评估汇编,北京:地震出版社,1-560
    Ambraseys, N. N.,1985. Intensity-attenuation and magnitude-intensity relationships for northwest European earthquakes. Earthq. Eng. Struct. D.,13:733-778.
    Bakun W H, Wentworth C M. Estimating earthquake location and magnitude from seismic intensity data. Bull Seism Soc Amer,1997,87(6):1502-1521.
    Bakun W H. Seismic Activity of the San Francisco Bay Region. Bull Seism Soc Amer,1999,89(3): 764-784.
    Bakun W H, A C Johnston, M G Hopper. Estimating Locaitons and Magnitudes of Earthquakes in Eastern North America from Modified Mercalli Intensities. Bull Seism Soc Amer,2003,93(1): 190-202.
    Bakun W H. Magnitude and location of historical earthquakes in Japan and implications for the 1855 Ansei Edo earthquake. J.Geophys.Res,2005,110:B02304
    Bakun W H. Estimating Location and Magmotudes of Earthquakes in the Southern California from Modified Mercalli Intensities. Bull Seism Soc Amer,2006,96(4):1278-1295.
    Bakun W H, Oona Scotti. Regional intensity attenuation models for France and the estimation of magnitude and location of historical earthquakes. Geophysics J Int,2006,164:596-610.
    Diane I Doser. Estimating Magnitude and Location of Alaskan Earthquakes Using Intensity Data. Bull Seism Soc Amer,2009,99(6):3430-3453.
    Frankel A. Implications of felt area-magnitude relations for earthquake scaling and the average frequency of perceptible ground motion. Bulletin of the Seismological Society of America, 1994,84:462-465.
    Gutenberg, B. and Richter, C. F.,1956, Earthquake Magnitude, Intensity, Energy, and Acceleration (Second Paper), BSSA,46(2),105-146
    Gasperini P, Bernardini F, Valensise G, eta.1999, Defining seismogenic sources from historical earthquake felt reports. Bull Seism Soc Amer,89(1),94-110
    Gasperini, P. and Ferrari, G.,2000. Deriving numerical estimates from descriptive information: The computation of earthquake parameters. Ann. Geofis.,43,729-746.
    Hinzen K G, Oemisch M.2001. Location and magnitude from seismic intensity data of recent and historic earthquake in the northern Rhine area, central Europe. Bull Seism Soc Amer,91(1): 40-56
    Johnston, A. C.,1996a. Seismic moment assessment of earthquakes in stable continental regions-Ⅰ. Instrumental seismicity. Geophys. J. Int.,124:381-414
    Johnston, A. C.,1996b. Seismic moment assessment of earthquakes in stable continental regions-Ⅱ. Historical seismicity. Geophys. J. Int.,125:639-678
    Johnston, A. C.,1996c. Seismic moment assessment of earthquakes in stable continental regions-Ⅲ. New Madrid 1811-1812, Charleston 1886, and Lisbon 1755. Geophys. J. Int.,126: 314-344
    Kovesligethy, R D.,1906, A makroszeizmikus rengesek feldolgozasa. Mathematikai es Termeszettudomanyi Ertesito,24,349-368
    Musson, R. M. W. and Jemenes, M. J.,2008. Macroseismic estimation of earthquake parameters. Sixth Framework Programme report.
    Mantyniemi, P.,2008. Earthquake of 4 November 1898 in northern Europe:New insights. J. Geophys. Res.,113:B11303, doi:10.1029/2007JB005461.
    Peruzza L.,1992, Procedure of macroseismic epicentre evaluation for seismic hazard purposes. Proceedings of the XXIII General Assembly of the ESC, Prague.,2:434-437.
    Richter, C. F.,1958. Intensities and isoseismals, in elementary seismology. San Francisco:W. H. Freeman,135-149.
    Sondergeld C H, Estey L H,1981, Acoustic emission study of microfracturing during the cyclic loading of westerly granite. J Geophys Res,86(B4):2915-2924.
    Toppozada, T. R.,1975. Earthquake magnitude as a function of intensity data in California and Western Nevada. Bull. Seism. Soc. Am.,65,1223-1238
    Teves-Costa, P. and Batllo, J.,2011. The 23 April 1909 Benavente earthquake (Portugal): macroseismic field revision. J. Seismol.,15:59-70. DOI 10.1007/s10950-010-9207-6.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700