地震模型正演

地震模型正演
地震模型正演

地震模型正演与反演简介

一、地震模型正演(seismic forward modeling)的概念

如果我们已知地下的地质模型,它的地震响应如何?地震模型正演就是通过室内模拟得到地质模型对于地震波的响应。地震模型正演包括物理模拟和数值模拟,数值模拟就是应用相应的地球物理方程和数值计算求解已知的地质模型在假定激发源的作用下的地震相应。通常,我们针对特定的勘探区块,应用期望或实际的采集参数通过地震正演模拟野外地震采集,得到单炮记录,再通过速度分析、动校正、叠加、偏移等处理得到成像数据。图1为Marmousi速度模型,图2为正演得到的炮集记录,图3为正演得到的叠加剖面。

图1 Marmousi模型

图2正演炮集

图3 正演叠加剖面

二、数值模型正演方法

通常,我们提到的模型正演为数值模拟的模型正演,目前常用的数值模拟地震模型正演方法包括基于射线原理的射线追踪法,以及基于波动方程的有限差分法、有限元法、积分方程法、快速傅里叶变换法和拟谱法等。射线追踪法主要反映地震波的运动学特征,有限差分、有限元法则适合复杂地质构造的正演模拟,积分方程法涉及复杂的数学推导,快速傅里叶变换法在频率域计算得到正演数据。

三、数值模型正演的步骤

数值模拟求解地震模型正演问题的步骤主要包括以下三个方面:

1) 地质建模,根据研究对象和问题建立地球物理或地质模型;

2) 数学建模,根据应用的物理手段和地球物理模型建立相应的数学

模型;

3) 模拟计算,选择正演计算方法,编写计算程序进行数值模拟计算。

四、什么是地震反演

地震反演技术就是充分利用测井、钻井、地质资料提供的丰富的构造、层位、岩性等信息,从常规的地震剖面推导出地下地层的波阻抗、密度、速度、孔隙度、渗透率、沙泥岩百分比、压力等地球物理信息。反演就是由地震数据得到地质模型,进行储层、油藏研究。

地震资料反演可分为两部分:

1)通过有井(绝对)、无井(相对)波阻抗反演得到波阻抗、速度数据体。

2)利用测井、测试资料结合波阻抗、速度数据进行岩性反演,得到孔隙度、渗透率、砂泥百分比、压力等物理数据。

图4为包括火山锥、火山溢流相的特殊岩性体偏移地震剖面,从常规剖面上地质工作者很难进行岩性体的解释和追踪。图5为对应的反演剖面,显然,反演剖面更有利于地质体的解释、识别、和横向追踪。

图4偏移地震剖面

图5地震反演剖面

五、常用地震反演方法

当今的地震反演方法,一是建立在波动理论基础之上的。由其固有的不适定性,反演的效果尚不如以褶积模型为基础的方法,加之这种方法尚不够完善而未能得到普遍推广。二是建立在褶积模型基础上的,即叠后资料反演。叠后地震资料反演可分为三种基本的类型(Brain.h.Russell):带限反演,或称常规递推法、稀疏脉冲法以及以模型为基础的反演方法。如地震岩性模拟、广义线性反演、宽带约束反演等。

稀疏脉冲法包括最大似然反褶积、L1范数反褶积、最小熵反褶积、最大熵反褶积、同态反褶积等,稀疏脉冲反演是基于脉冲反褶积基础上的递推反演方法,其基本假设是地层的强反射系数是稀疏分布的。从地震道中根据稀疏的原则提取反射系数,与子波褶积后生成合成地震记录;利用合成地震记录与原始地震道残差的大小修改参与褶积的反射系数个数,再作合成地震记录;如此迭代,最终得到一个能最佳逼近原始地震道的反射系数序列。该方法适用于井数较少的地区,其主要优点是能够获得宽频带的反射系数,较好地解决地震反演的多解性问题,从而使反演结果更趋于真实。

常规递推法与稀疏脉冲反演法主要是利用反褶积方法来恢复反射系数序列,由经过标定的反射系数序列递推出相对波阻抗,然后加上从声波测井和地质模型中得到的低频分量,最终得到反演波阻抗。这两类方法的主要缺陷是选择可靠低频信息较为困难,由反射系数递推波阻抗过程中误差积累快,当反射系数存在较大误差时,递推出来的波阻抗剖面会面貌全非。

此外,经过反褶积处理的结果,并不代表真正的反射系数序列,稀疏脉冲法在地质结构复杂的条件下使用效果很差。其精度也难以满足储层预测、油藏描述的需要。

以模型为基础的反演方法以测井资料为约束条件,采用正、反演结合进行迭代,求取地下波阻抗将反演方法推向非线性问题。这种新方法利用了测井资料的高频和低频信息,大幅度拓宽了地震信号的频带,可以更好地获得薄层、薄互层的波阻抗信息。因而表现出强劲的发展势头,是目前国内外各软件公司重点发展的技术,也是油田开发阶段进行储层预测、油藏描述的主要应用技术。

叠前反演主要指AVO(AVA,振幅随炮检距或入射角变化)反演,是继亮点技术后出现的又一种油气直接识别技术,Zoeppritz方程描述了振幅与入射角的关系。精确的Zoeppritz方程是AVO烃类检测技术的基础。随着地球物理技术的飞速发展及利用地震信息进行岩性勘探与含油气储层的研究,

直接利用地震信息找油找气成为人们普遍关注的热点,亮点技术的出现,为直接找油、找气揭示了广阔的前景;同时真假亮点想象使人们陷入勘探的误区,如何识别真假亮点,即如何从亮点的振幅异常中识别出只与含气有关的亮点,多年来人们渴望找到一种方法来解决这一问题,八十年代AVO 技术的出现,对解决这一问题有了可能,近年来AVO技术日趋成熟,在国外AVO烃类检测技术的广泛应用,为油气田的开发带来了巨大的经济效益。

地震反射信号的振幅随炮检距的变化而变化。反射系数对反射角度的依赖关系反映了这种变化,在某些沉积环境中,这种振幅变化可能提供了烃类存在的重要信息,如当碎屑岩空隙中含气时,纵波速度降低,而横波速度不变,这造成地震反射波振幅与炮检距(或地震波入射角)之间的关系发生变化产生振幅异常。通常共中心点叠加压制了振幅随入射角变化的信息,因为叠加中每个同相轴的振幅表示所有不同炮检距上地震反射波的平均振幅。AVO分析处理技术试图研究振幅随入射角的变化关系,研究这种关系的变化,对储层的含油气性进行预测,特别适合于预测含气储层。求得的AVO属性以适当的方式显示,并对其进行解释、研究,找出区域上油气的分布规律。

地震检测模型

楚雄师范学院 2014年“雁峰杯”数学建模竞赛论文 题目地震检测 姓名杨子月 学院数学与统计学院 专业数学与应用数学 2014年5月28日

地震检测模型 摘要 继2008年5月12日在四川汶川大地震之后,2013年4月22日四川雅安又发生了一次7.0级地震,这些重大自然灾害,给我们每一位中国人带来了巨大的伤痛,痛定思痛,我们应该为减少震后灾害做些事情。当地震发生时,震中位置的快速确定对第一时间展开抗震救灾起到非常重要的作用,而震中位置可以通过多个地震观测站点接收到地震波的时间推算得到。 现已采集到某地观测的30个指标的数据,和该地区该时期内已发生地震的经纬度、地震波到达的时间的数据。科学地截取这些数据的有用片段,对数据进行合理地预测处理,用数学方法计算出地震的中心位置。 关键词:地震检测经纬度地震波到达时间震源中心

一、问题重述 假设你是一位地震学家,在某地部署了30座地震台。这些地震台装备了测量和记录地质运动的设备。现已采集了这30座地震台的坐标和某次地震时这些的地震台测得的地震运动到达时间t,现在我们需要建立一个数学模型求出这次地震中心的坐标M(x,y)。 二、模型假设 1、假设震源在地下,发生地震之后地震波沿着各个方向匀速传播,且在传播过程中速度保持不变。 2、假设地震波在各种介质中的传播速度相等。 3、假设地震发生的区域范围内时差为零。 4、、假设由于其他因素而引起10多个指标数据的变化以及非正常波动可以忽略不计。 5、假设地震的前兆指标的数据特征符合一定的概率统计分布。 6、地形各观测点没有剧烈变化。 通过以上条件虽然不能精确求出地震发生的地点,但是可以建立一种在空间和时间上准确模拟地震发生以及预测的模型机制,对于地震预报及防治有很大的现实意义。地震源可能在地下,地震发生之后,地震波从震源点开始以球面方式沿各个方向传播,在空间和时间上是一个三维的立体模型结构。 三、符号说明及名词解释 3.1符号说明 震中位置 M(x,y) 经度 x(度) 纬度 y(度) 震源深度 h(千米) 地震波在各种介质中的传播速度v(千米/秒) 地震波到达时间 t(秒) 3.2 名词解释 地震波:地震被按传播方式分为三种类型:纵波、横波和面波。纵波是推进波,地壳中传播速度为5.5~7千米/秒,最先到达震中,又称P波,它使地面发生上下振动,破坏性较弱。横波是剪切波:在地壳中的传播速度为3.2~4.0千米/秒,第二个到达震中,又称S波,它使地面发生前后、左右抖动,破坏性较强。面波又称L波,是由纵波与横波在地表相遇后激发产生的混合波。其波长大、振幅强,只能沿地表面传播,是造成建筑物强烈破坏的主要因素。[1]

全球地震采集技术发展趋势分析

全球地震采集技术发展趋势分析 2014-12-04 15:16:00 1 地震采集声波多普勒传感器波场旋转分量多源地震采集双圆形采集吴伟分布式震源组合dsa 文|吴伟等 中石化石油物探技术研究院

地震技术结合多学科高新技术发展最新成果,在油气勘探开发中发挥了日益重要的核心支撑作用。在遵循采集、处理、解释一体化的发展思路下,借助于先进仪器装备和各种采集新技术的不断推出,地震采集技术近年来的发展明显加快,正向着适应更恶劣地表条件、更复杂地下构造和更隐蔽含油气圈闭勘探需求的精细采集方向发展。 采集理念由过去的追求共中心点叠加次数向“以记录波场为中心”转变,采用24位超万道地震仪、数字检波器加网络技术支撑,精细表层调查和模型驱动的采集设计,进行单点接收、大动态范围、无线化传输、超多道记录、小面元网格、高覆盖次数、高品质震源、多分量接收、全方位信息、环保型作业的高密度三维地震全波场采集,不断提高地震资料的纵、横向分辨率和有效信息的精确度。多源地震等高效地震采集技术的出现,提高了采集效率并给数据处理技术带来变革。本文从采集理念、新型传感器、多源高效地震、海上地震新技术等方面总结了近年来地震采集技术的发展动态,展望采集技术的发展方向。限于篇幅,没有详细介绍目前在采集中已普遍使用的一些方法技术,如采集前的模拟技术(局部照明分析、振幅分析等),激发接收条件的分析与改进,地震采集脚印的模拟与采集效应的消除,山前带、沙漠、滩海等特殊环境的采集等。 一些新的技术在采集中的应用,如稀疏采样(包括随机采样、压缩采样等在空间、时间域中突破采样定理要求的采集)、光纤传输等,文中也没有涉及。 1 “以记录波场为中心”的采集新理念

仿真地震模拟

Part 01公司简介

北京盛世民安科技发展有限公司 ↘北京盛世民安科技发展有限公司始于2009年,总部位北京经济技术发区(国 家级),是一家集的设计、研发、制造 的集成产业公司。公司主要技术人员主 要来自各大公司、各大院校,在研制、 生产、销售上有多年的综合管理经验, 多年的产品销售和优质售后服务,让很 多的行业对我公司产品有了更好的认识 和认可,在我公司自己的不断努力下和 客户们的拥戴下,在上海、南京、西安、成都等地区分别设立了办事处,为全国 24小时提供优质服务机制。公司涉及软 件、动漫、自动化、机械、液压、电子 等多个行业。↘公司主要产品有:模拟灭火、模拟报警、烟雾逃生、地震体验平台、地震体验小 屋、地震科普展品、仿真动感平台、三 自由度平 ↘台、六自由度平台等。客户已经遍布各大科技馆、展览馆、房地产、政府、企 业、文化影视等诸多领域和各行各业。

part02建设背景

↘地震模拟的宗旨就是真实的再现地震过程并使观众如身临其境般感受地震,其关键就是“真实”。基于地震模拟的宗旨,该产品通过模拟普通小屋模拟地震发生前、发生过程中和震后的整个过程,并提示观众在地震中如何躲避、如何逃生。 ↘地震体验模拟器包括整体建筑及屋内陈设、由虚拟显示技术和装饰模拟出的房屋窗体、房屋地板下面的六自由度液压运动平台及控制台。整体建筑基于戈壁滩小屋的构造。 展品通过小屋实体、液压运动平台、控制台并辅以视频及音响系统组成一个立体的地震环境,使观众与整个环境融为一体,从而全方位展示地震的全过程,使观众真实感受地震。其中液压运动平台按控制台的操作能提供相当于3-8级地震效果,包括横波与纵波,可以体验地震发生前、地震中及余震。

地震成因及风险模型

摘要首先,介绍了地震形成的自然因素和非自然因素,并对其发生原因进行了分析和研究;其次,通过对1999-2009年间的地震现场的灾害调查资料和损失评估的资料进行研究,对其进行分析与处理,采取以实际烈度区作为分配单元,建立适用于县级区域小尺度的地震风险分析模型;最后基于烈度的地震分析模型与基于建筑物易损性地震分析模型预测结果进行比较,由此说明在县级区域小尺度上,地震风险分析模型的适用性。 关键词分析模型烈度地震灾害损失评估 A Study on the Analysis Mode for the Causes and Risks of Earthquakes//Wang Xingyu[1],Chen Peng[2]* Abstract We first introduced the earthquake causes which con-sist of the natural factors and unnatural factors.Second,we use data from the investigation of earthquake disasters and the inf-ormation of disaster losses1999-2009,gathering respective sec-tions strength data of the population,the per capita GDP and the area of land,from evaluation information of earthquake disa-ster and statistics of the earthquake province when the earthqu-ake happens.By analyzing and handling the above data,use virtual broken-level areas as allocation units and establish risk analysis model for medium and small scale earthquake of https://www.360docs.net/doc/9414397488.html,parison the two forecast results from the earthquake analysis model based on strength and building damage to prove that the earthquake analysis model based on strength is more suitable for medium and small scale county.Finally,we introduced some knowledge about self-protection when the earthquake happens. Key words analysis model;strength;earthquake disaster;evalu-ation of losses First-author's address Changkou Middle School of Fuyang City,311400,Hangzhou,Zhejiang,China 1引言 地震作为中国灾害中破坏力最强,损失最严重的灾种,而被研究者所重视。因此,对地震成因的分析以及如何减轻地震风险是这篇文章的主要内容。随着城市化进程日益加快,承灾数量不断增加,但是,灾害评价分析模型还不能满足现代应急的需求。以往的地震评价研究模式,主要是针对地震风险分析模型在地震减灾中长期规划研究,但不适用于突发性地震事件。如何快速、准确地对突发性地震事件发生前或者发生时做出应急管理是决策者面临的重大难点。本文提出了一个地震风险分析模型,该模型能够很好地满足地震应急需求,同时在一定程度上满足抗震风险分析需求。 2地震发生的自然原因 地震是地壳运动的一种特殊表现形式,也是极为常见的地质现象。地震有多种成因,根据其成因分为构造地震、火山地震和陷落地震三种主要类型。 2.1构造运动 构造地震是由地壳运动所引起的地震。一般而言,地壳运动是长期的、缓慢的,一旦地壳所积累的地应力超过了组成地壳岩石极限强度时,岩石就要发生断裂而引起地震,也就是说地应力从逐渐积累到突然释放时才发生地震。构造地震是一种活动频繁、影响范围大、破坏力强的地震,世界上最多(90%以上)和最大的地震都属于构造地震[1]。 2.2火山运动 由于火山活动时岩浆喷发的冲击力或热力作用而引起的地震,称为火山地震。火山地震一般较小,数量约占地震总数的7%左右,地震和火山往往存在关联,火山爆发可能会激发地震,而发生在火山附近的地震也可能引起火山爆发,通常发生在板块的生长边界。其特点是震源常限于火山活动地带,一般深度不超过10公里的浅源地震,震级较大,多属于没有主震的地震群型,影响范围小。 2.3陷落运动 陷落地震是由于岩层大规模崩塌或陷落而引起的地震。这种地震为数很少,只占地震总数的3%左右,一般震级较小,影响范围不大,地震能量主要来自重力作用。陷落运动主要发生在石灰岩或其他易溶的岩石地区,由于地下溶洞不断扩大,洞顶崩塌,引起震动,导致矿洞塌陷或大规模山崩、滑坡等也可能导致这类地震发生。 3地震发生的非自然原因 在特定地区由于某种地壳外界的非自然因素而引起的地震,称为诱发地震。这些外界因素可能是地下核爆炸、陨石坠落、油井灌水等,其中最常见的是水库地震。水库蓄水,石油和天然气、盐卤、地下热(汽)储的开发,废液处理和油田开采中的深井注水,钻进过程中的井漏,矿山抽排水,固 ([1]杭州市富阳市场口中学浙江·杭州311400; [2]吉林师范大学旅游与地理科学学院吉林·四平136000) 中图分类号:K909文献标识码:A文章编号:1672-7894(2012)15-0081-03 81

地震预测模型doc

精心整理2011年赣南师院数学建模竞赛选拔赛 题目地震预测模型 摘要: 本文前三个任务主要考虑是各指标的变化对地震发生问题的影响,通过对各指标数据量的分析建立相应的模型,并对任务四和任务五给出了合理的解答。 针对任务一:我们从原始数据中计算出各项指标的日均值,绘制出各指标分年度的时间序列图, 磁波幅度 。 关键词: 一·问题的重述 1.1背景分析 地震是地壳快速释放能量过程中造成的振动。虽然预测地震是世界性难题,但迄今科学界普遍认为,有可能反映地震前兆特征的指标可能不少于10个。已经有专业仪器在多个定点实时按秒记录这些指标的数据,期望通过对记录数据的分析研究找到地震的前兆特征。 现已采集到某地2005年1月1日至2010年6月30日按小时观测的10多个指标的数据,和该地区该时期内已发生地震的时刻、经纬度、震级及震源深度的数据。这些数据中隐藏着地震发生的前兆特征。科学地截取这些数据的有用片段,对数据进行合理地预处理,用数学方法揭示地震前兆

的数据特征,是一项很有意义的研究工作。 题给数据中的这10多个指标,究竟哪些与地震的发生有关,有何种关系,是单一关系还是复合关系;除这10多个指标外还有哪些因素及含题给指标在内的哪些指标的哪种数学模型更能反映地震的前兆特征等等,人们迄今仍不很清楚,需要进行深入地研究。地震数据的观测是持续进行的,随着时间的推移数据的规模会不断扩大。从中挖掘地震的前兆特征,必须有合理的数学模型,也必须有科学高效的算法分析平台。因此,需要我们结合附件中给出的实际记录数据,尝试完成以下任务。 1.2任务的提出 任务一:分析数据特征,建立数学模型以度量各指标对地震发生的敏感程度。 越大 任务三:中要结合题给数据,建立数学模型来研究地震发生前的数量特征。主要运用贝叶斯判别分析法进行建模,对已给数据进行先验信息、后验信息分析。 任务四:要将计算程序集结成地震数据分析平台,能够完成其它地震数据的分析,并能自动输出前任务的重要分析结果。 任务五:是针对进一步的研究设想写一篇切实可行的报告。 三·问题的基本假设 (1)地震监测点的监测设施能正常运转; (2)地震监测设施周围不存在影响其工作效能的干扰源,如飞机场、发电厂等;

叠加地震记录的相移波动方程正演模拟数值模拟实验共22页

《地震数值模拟》实验报告 一、实验题目 叠加地震记录的相移波动方程正演模拟

二、实验目的 1.掌握各向同性介质任意构造、水平层状速度结构地质模型的相移波动方程正演模拟基本理论 2.实现方法与程序编制 3.由正演记录初步分析地震信号的分辨率。 三、实验原理 1、地震波传播的波动方程 设(x,z)为空间坐标,t为时间,地震波传播速度为v(x,z),则二位介质中任意位置、任意时刻的地震波场为p(z,x,t):压缩波——纵波。则二维各向同性均匀介质中地震波传播的遵循声波方程为 2、傅里叶变换的微分性质 p(t)与其傅里叶变换的P(w)的关系: 3、地震波传播的相移外推公式 令速度v不随x变化,只随z变化,则利用傅里叶变换微分性质把波动方程(变换到频率-波数域,得: 4、初始条件和边界条件 按照爆炸界面理论,反射界面震源在t=0时刻同时起爆,此时刻的波场就是震源。根据不同情况,可直接使用反射系数脉冲或子波作震源。如果直接使用反射系数作震源脉冲,则初始条件可表示为: 5、边界处理

(1)边界反射问题 把实际无穷空间区域中求解波场的问题化为有穷区域求解时,左右两边使用零边界条件。物理上假设探区距Xmin与Xmax两个端点很远,在两个端点上收到的反射波很弱。但是,上述条件在实际中不能成立,造成零边界条件反而成为绝对阻止波通过的强反射面。在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。 (2)边界强反射的处理 镶边法、削波法、吸收边界都能有效消除边界强反射。 削波法就是在波场延拓过程中,没延拓一次,在其两侧均匀衰减到零,从而消除边界强反射的影响。假设横向总长度为NX,以两边Lx道吸波为例,有以下吸波公式: 四、实验内容

高阶统计量地震子波估计建模

2006年10月 第41卷 第5期  3山东省东营市中国石油大学(华东)信息与控制工程学院,257061本文于2005年12月21日收到,修改稿于2006年5月12日收到。 本项研究受高等学校博士学科点专项科研基金(No.20020008004)部分资助。 ?处理方法? 高阶统计量地震子波估计建模 戴永寿3①② 郑德玲① 魏 磊② 霍志勇② (①北京科技大学信息工程学院;②中国石油大学(华东)信息与控制工程学院) 摘 要 戴永寿,郑德玲,魏磊,霍志勇.高阶统计量地震子波估计建模.石油地球物理勘探,2006,41(5):514~518,540 本文在反射系数序列为非高斯、平稳和统计独立的随机过程,地震子波为非因果、混合相位的假设条件下,分别应用滑动平均(MA )和自回归滑动平均(ARMA )模型对地震记录进行建模,并采用运算代价较小的基于高阶累积量的线性化求解方法———累积量矩阵方程法进行了子波提取和模型适应性的研究。数值模拟结果和实际地震数据处理结果表明:自回归滑动平均(ARMA )模型比滑动平均(MA )模型具有参数节省、模型更为高效的特点;累积量矩阵方程法可以有效地压制加性高斯噪声,但对累积量样本估计的准确性要求较高;如果累积量样本估计的误差和方差适度,结合自回归滑动平均(ARMA )模型描述的累积量矩阵方程法可以高效、准确地估计出地震子波。 关键词 高阶累积量 子波 自回归滑动平均(ARMA ) 滑动平均(MA ) 建模 1 引言 作为地震资料反褶积处理、波阻抗反演以及正演模拟的基础工作,准确的地震子波估计对于高分辨率、高信噪比、高保真度的地震勘探数据处理具有极为重要的意义。统计性子波提取方法的基本原理是首先对反射系数序列的分布做某种假设,然后利用地震记录的统计信息进行子波估计。在没有任何先验知识的情况下,通常假设反射系数序列为一个非高斯、平稳和统计独立的随机过程,假设子波为一个非因果、非最小相位系统,加性噪声为高斯色噪声。因此在利用地震记录的统计信息进行子波估计时,其高阶累积量不仅能保留系统的相位信息,而且能较好地压制高斯色噪声,显示出此法的优越性。 近年来,基于高阶累积量的参数化子波估计方法得到了快速发展。Lazear [1]首先引入滑动平均(MA )模型描述地震记录,然后将子波四阶矩和地震资料的四阶累积量在最小均方误差意义下进行拟合,并用梯度下降法求解目标函数。随后,Velis 等人[2]及尹成等人[3]试图应用特性更好的全局最优化 方法解拟合函数,但求解效率普遍较低。石殿祥等 人[4]基于高阶累积量研究了非最小相位子波提取问题,虽取得了一定的成果,但依然沿用了滑动平均(MA )模型来描述地震记录。 本文分别采用滑动平均(MA )模型和自回归滑动平均(ARMA )模型来描述地震记录,并借助基于高阶累积量的线性化参数估计方法———矩阵方程法求解模型参数,最终精确估计了地震子波。 2 地震记录的滑动平均(MA)模型描 述及矩阵方程法子波提取 地震记录y (n )可视为一个零均值的平稳随机过程,且符合如下褶积模型 y (n )= ∑q i =0 w (i )r (n - i )+v (n ) =w (n )3r (n )+v (n ) (1) 式中:w (n )为地震子波;r (n )为反射系数序列;v (n )为环境噪声。显然,式(1)符合典型的滑动平均(MA )模型表达式,因此可以把地震记录看作是有限脉冲响应(FIR )系统的含噪输出。对于上述模型有如下假设:

实用文库汇编之数学建模地震预测模型

*实用文库汇编之 * 题目:地震预测数学建模 姓名:张志鹏 学号:12291233 学院:电气工程学院 姓名: 赵鑫 学号:10291033 学院:电气工程学院 数学建 模竞赛 论文

姓名:张书铭学号:12291232 学院:电气工程学院 目录 摘要 (3) 一、问题重述 (4) 二、问题的分析 (4) 三、建模过程 (5) 问题1:地震时间预测 (5) 1、问题假设 (5) 2、参数定义 (6) 3、求解 (6) 问题2:地震地点预测 (7) 1、问题假设: (7) 2、参数定义 (7) 3、求解过程: (7) 四、模型的评价与改进 (10) 参考文献 (11)

摘要 大地振动是地震最直观、最普遍的表现。在海底或滨海地区发生的强烈地震,能引起巨大的波浪,称为海啸。在大陆地区发生的强烈地震,会引发滑坡、崩塌、地裂缝等次生灾害。对人们的生产生活成巨大影响,严重威胁人们的生命和财产安全,所以,对地震的预测是十分必要的。 本文根据从1900年以来中国发生的八级以上地震的时间和地点分析,利用合理的数学建模方法,对下一次中国可能发生的八级以上地震的和时间和地点进行合理的预测。建模方法分为对于时间的预测和地点的预测两个方面。 问题1:对于时间的预测 采用的方法为指数平滑法,它是通过计算指数平滑值,配合一定的时间序列预测模型对现象的未来进行预测。其原理是任一期的指数平滑值都是本期实际观察值与前一期指数平滑值的加权平均。 问题2:对于地点的预测 根据长久的数据表明,八级以上地震主要发生在东经70°——110°,北纬20°——50°这个范围内,据此将整个地震带划分为100个区域,按顺序进行编号。建立时间与地震区域编号的数学模型,利用线性回归的方法对下次地震地点预测。

地震模型正演

地震模型正演与反演简介 一、地震模型正演(seismic forward modeling)的概念 如果我们已知地下的地质模型,它的地震响应如何?地震模型正演就是通过室内模拟得到地质模型对于地震波的响应。地震模型正演包括物理模拟和数值模拟,数值模拟就是应用相应的地球物理方程和数值计算求解已知的地质模型在假定激发源的作用下的地震相应。通常,我们针对特定的勘探区块,应用期望或实际的采集参数通过地震正演模拟野外地震采集,得到单炮记录,再通过速度分析、动校正、叠加、偏移等处理得到成像数据。图1为Marmousi速度模型,图2为正演得到的炮集记录,图3为正演得到的叠加剖面。 图1 Marmousi模型 图2正演炮集

图3 正演叠加剖面 二、数值模型正演方法 通常,我们提到的模型正演为数值模拟的模型正演,目前常用的数值模拟地震模型正演方法包括基于射线原理的射线追踪法,以及基于波动方程的有限差分法、有限元法、积分方程法、快速傅里叶变换法和拟谱法等。射线追踪法主要反映地震波的运动学特征,有限差分、有限元法则适合复杂地质构造的正演模拟,积分方程法涉及复杂的数学推导,快速傅里叶变换法在频率域计算得到正演数据。 三、数值模型正演的步骤 数值模拟求解地震模型正演问题的步骤主要包括以下三个方面: 1) 地质建模,根据研究对象和问题建立地球物理或地质模型; 2) 数学建模,根据应用的物理手段和地球物理模型建立相应的数学 模型; 3) 模拟计算,选择正演计算方法,编写计算程序进行数值模拟计算。 四、什么是地震反演 地震反演技术就是充分利用测井、钻井、地质资料提供的丰富的构造、层位、岩性等信息,从常规的地震剖面推导出地下地层的波阻抗、密度、速度、孔隙度、渗透率、沙泥岩百分比、压力等地球物理信息。反演就是由地震数据得到地质模型,进行储层、油藏研究。

地震预测研究发展战略几点思考

地震预测研究发展战略几点思考 岳明生 (2005年1月6日) 同志们:  2004年12月13—14日,中国地震局在北京召开了《地震预报发展战略研讨会》,近80位院士、专家参加了会议。会议分二个阶段进行,第一阶段是专题报告,第二阶段是自由发言。专题报告观点明确,内容丰富,论述详实。自由发言踊跃,提出了很多很好的建议。陈建民局长、宋瑞祥同志分别到会发表了重要讲话。与会院士、专家一致认为,这是一次学术空气、讨论氛围都很好的会议。院士、专家们的很多建议对制定我国地震预报发展战略都具有十分重要的参考价值。 会议形成了三点共识,一是地震预报是十分困难的,要比我们料想的困难程度难得多。但近四十年来的观测事实说明,地震是有规律、有前兆的,最终实现地震预测是可能的。为此,我们必须正视困难,增强信心,坚定开展地震预测研究和预报工作的决心。二是我国近四十年来的地震预测研究和地震预报工作是有成绩的,其在全世界的先进地位是有目共睹的,但监测、预报、科研越来越脱节、基础研究落后、创新不足的问题已经严重地影响到地震预报工作的发展,必须引起高度重视。三是突出重点,重建地震预测研究实验场,带动我国地震预测研究和地震预报工作向纵深发展。正确处理继承与发展的关系,已有的监测预报成果是我们的宝贵财富,是发展的基础,既要继承好,又要不断向纵深发展,同时必须不断地引进新理论、新技术、新思路;全面加强地震预测研究,建立地震预测研究可持续发展和监测、预报、科研有机结合机制,既要重视发震机理的研究,又要重视地震前兆现象的研究;实行开放式的地震预测研究,欢迎国内外、行业内外和广大人民群众参与我国的地震预测研究,鼓励多路探索,百花齐放,百家争鸣。  下面我对会议形成的三点共识谈点认识。 一、正视困难,增强信心,坚定开展地震预测研究的决心 地震预测科学是一门以科学假说为理论依据,以观测为基础的预测科学。根据马克思主义认识论,我们可以将假说(模型)、观测、地震建立如下关系:  如图1所示,人们根据对地震成因和地震前兆机理的已有认识,建立模型;再根据模型选择观测方法,制定观测方案,确定预期的地震前兆指标;地震发生后把实际观测到的地震前兆现象与预期的地震前兆指标进行对比,找出差距;根据差距再去修改模型,然后再重复模型-观测-地震的循环过程,以至无穷,一步步缩小预期的地震前兆指标与实际观测到的地震前兆现象的差距,最后实现地震的准确预测。   图1 地震震源形成及其演化的模型,自Reid1911年提出地震弹性回跳理论后,又相继提出了地震包体理论(B.T.Brady)、红肿理论(傅承义)等多种理论模型,但都不能全面解释震源形成和演化过程。现在开展的地震前兆观测有地震活动性、地壳形变、地下水、地电、地磁、应力-应变、气象、诱发因子、宏观异常等十几类近百种方法,也发现了大量的异常现

数学建模地震预测模型

数学建模竞赛论文 题目:地震预测数学建模 姓名:张志鹏学号:12291233 学院:电气工程学院 姓名:赵鑫学号:10291033 学院:电气工程学院

姓名:张书铭学号:12291232 学院:电气工程学院

目录 摘要 (3) 一、问题重述 (4) 二、问题的分析 (4) 三、建模过程 (5) 问题1:地震时间预测 (5) 1、问题假设 (5) 2、参数定义 (6) 3、求解 (6) 问题2:地震地点预测 (7) 1、问题假设: (7) 2、参数定义 (8) 3、求解过程: (8) 四、模型的评价与改进 (12) 参考文献 (13)

摘要 大地振动是地震最直观、最普遍的表现。在海底或滨海地区发生的强烈地震,能引起巨大的波浪,称为海啸。在大陆地区发生的强烈地震,会引发滑坡、崩塌、地裂缝等次生灾害。对人们的生产生活成巨大影响,严重威胁人们的生命和财产安全,所以,对地震的预测是十分必要的。 本文根据从1900年以来中国发生的八级以上地震的时间和地点分析,利用合理的数学建模方法,对下一次中国可能发生的八级以上地震的和时间和地点进行合理的预测。建模方法分为对于时间的预测和地点的预测两个方面。 问题1:对于时间的预测 采用的方法为指数平滑法,它是通过计算指数平滑值,配合一定的时间序列预测模型对现象的未来进行预测。其原理是任一期的指数平滑值都是本期实际观察值与前一期指数平滑值的加权平均。 问题2:对于地点的预测 根据长久的数据表明,八级以上地震主要发生在东经70°——110°,北纬20°——50°这个范围内,据此将整个地震带划分为100个区域,按顺序进行编号。建立时间与地震区域编号的数学模型,利用线性回归的方法对下次地震地点预测。

基于MATCOM 平台的VC++与MATLAB混合编程地震正演模拟(完整版)

基于MATCOM 平台的VC++与MATLAB混合编程地震正演模拟 摘要:正演模拟是地震勘探的重要方法之一。为解决基于Matcom4.5的VC++与Matlab接口程序的地震正演模拟问题,通过用Matcom4.5与VC++结合编程的新方法实现EDA介质中的地震波正演模拟的实例,结果显示通过Matcom实现VC++与Matlab混合编程的可行性并提高了编程效率。 关键词: Matcom;VC++; Matlab;混合编程;正演模拟 MIXED PROGRAMMING OF SEISMIC FORW ARD MODELING BETWEEN VC++ Abstract: Forward modeling is one of the important means in seismic exploration. This paper in order to solve the interface program of seismic forward modeling between Matlab and VC++ based on Matcom4.5,by using an example which achieved the results of simulate the wave equation in EDA media by a new method which is mixed programming between Matcom4.5 and VC++,the results show the advantage of efficiency by mixed programming between Matlab and VC++ based on Matcom4.5. Keywords: Matcom;VC++;Matlab;Mixed programming;Forward modeling 0引言 近年来,随着勘探程度深入化和勘探目标的复杂化,地震勘探无论在采集、处理还是解释技术方面都面临越来越多的挑战。复杂构造油气田、岩性油气藏和裂缝性油气藏等隐蔽性油气藏的勘探与开发越来越受到地球物理工作者的重视。为此,需要发展复杂介质的地震波传播理论,使其进一步接近地下介质的实际情况。正演模拟研究对复杂介质中波的传播具有指导意义[1]。 目前在地震数据处理中应用比较广泛的软件是MATLAB,但由于MATLAB本身的程序编写的一些局限,使得它在处理许多应用问题时灵活性不够。本文引入MATCOM软件平台,它的功能相当强大,可提供近千个MATLAB的基本功能函数,通过必要的设置,可以直接实现与C++的混合编程,而不必再依赖MATLAB;笔者基于MATCOM4.5的VC++与MATLAB的接口技术实现EDA介质中的地震波正演模拟,此方法提高了代码的使用效率,说明MATCOM与VC++混合编程的可行性,取得较好的效果,为技术人员提供了另一种实现正演模拟的途径,同时也为开展相应其他研究工作奠定了基础[2]。 1VC++与MATLAB接口 Visual C++是Windows平台下主要的应用程序开发环境之一,它能方便地实现软件开发,开发的系统具有界面友好、执行速度快、易维护和升级等优点。C++

基于MATLAB的地震正演模型实现[1]

基于MATLAB的地震正演模型实现 贾跃玮 (中国地质大学(北京) 北京100083) 摘 要 人工合成地震正演模型是进行三维模型计算的基础。针对地震勘探的原理,本文运用MATLAB强大数学计算和图像可视化功能,对一个三层介质模型制作了人工合成地震记录。文章首先说明了地震记录形成的物理机制,然后介绍了地质模型的构造及参数选择,最后针对该具体地质模型制作了合成地震记录。 关键词 地震;MATLAB;正演 0引 言 地震勘探就是利用地下介质弹性和密度的差异,通过观测和分析大地对人工激发地震波的响应,推断地下岩层的性质和形态的地球物理方法。地震勘探是钻探前勘测石油与天然气资源的重要手段,在煤田和工程地质勘查、区域地质研究和地壳研究等方面,也得到广泛应用。 人工合成二维地震模型记录是各种复杂地震模型正演计算的基础,是对地震勘探经典理论的忠实实现。在实际工作中,针对具体地质构造进行二维地震模拟能够有效帮助地球物理工作者在地震剖面上识别各种地质现象。MATLAB环境集编程、画图于一体,特别适合人工合成地震记录的快速实现。因此,我们在MATLAB环境下设计了一个三层地质模型,并对该模型模拟了地震记录,旨在可视化地观察地震波场记录特征并验证地震褶积模型。 1地震记录形成的物理机制 在地震记录上看到的波形是地震子波叠加的结果,从地下许多反射界面发生反射时形成的地震子波,振幅大小决定于反射界面反射系数的绝对值,极性的正负决定于反射系数的正负,到达时间的先后取决于界面深度和覆盖层的波速。若地震子波波形用S(t)表示,反射系数是双程垂直反射旅行时t的函数,用R(t)表示,地震记录f(t)形成的物理过程在数学上就可表示为:f(t)=S(t)3R(t)=∫0T S(τ)R(t-τ)dτ 地震子波和反射系数资料常常不易取得,因此计算时常做这样一些假设: (1)地质模型的建立是来自大量观察实际地质结构的经验性归纳总结。 (2)为了模型建立和计算过程中突出理论数值,去除了一些干扰因素,对一切衰减、噪声都不进行考虑。 (3)地层在横向上均匀,纵向上是由大量具有不同弹性性质的薄层构成。 (4)地震子波以平面波形式垂直入射到界面,各薄层的反射子波与地震子波形状相同,只是振幅及极性不同。 (5)所有波的转换、吸收及绕射等能量损失都不考虑。 基于以上这些假设条件进行地震记录合就必须已知地震子波以及地层的反射系数,而反射系数又主要由地层的波阻抗反映,所以必须首先获取地层的速度和密度资料。 速度资料可通过连续速度测井获得,密度资料可从密度测井获得,得不到密度资料时,可近似假定密度不变,以速度曲线代替波阻抗曲线来计算反射系数。加德纳根据实际资料提出了一个由速度推算密度的经验公式: ρ=0.23V0.25 (速度单位:英尺/秒) 或 ρ=0.31V0.25 (速度单位:米/秒)

地震模型1

地震模型 摘要 地震预报是当代自然科学领域里一个直接关系到人民生命安全和社会经济的发展,是科学难度很大的前沿课题。近几十年来,人民的努力虽然取得了丰硕的成果,预报的实践有些有限的成功。 在地震发生前有很多前兆性指标,如磁场强度、氡值、电压、雨量、水温等,这些指标都从不同的侧面反映了地震活动的各项特征。在正常情况下,也常常有些指标有明显的异常,而令一些指标并不出现异常。这些都给实际的前兆指标数据特征分析以及地震预报工作带来了困难。 本文针对地震影响因素多,数据多,联系弱的特点,分别建立了针对各个任务的数学模型,首先,处理初始数据时选取日平均变化来消除一天中数据的随机因素的影响,然后更具这些数据建立了地震发生前后各个指标的不同程度度量模型解决了地震对指标的敏感度的度量问题,进而找出了EW,气氡,水位,电压,雨量这些衡量地震的主要因素,并分析出了这些指标的重要程度,在解决找统计量时,利用上次任务所得的主要因素的组合来衡量一次地震的地震前兆指标数量特征,并得出了电压达到26.8754,气氡达到17.685左右时就有地震发生的可能性,考虑全面周到,效果较好,在模型改进中所提出的判别分析法,科学有效,对数据利用率较高。而除此之外的其他几个指标与震级有关,这样就这些指标的分析在任务三中分别建立了地磁前兆异常的动态从属函数模型以及非线性综合模型,进而找出了地震级数与这些因素之间的非线性关系,对于上述数学方法在地震预测预报中的应用没有达

到预期效果,须进一步研究改进。综合这些工作,我们应用易语言编写了分析地震数据的平台,通过分析地震数据生成了报告。 一、问题重述 地震是地壳快速释放能量过程中造成的振动。虽然预测地震是世界性难题,但迄今科学界普遍认为,有可能反映地震前兆特征的指标可能不少于10个。已经有专业仪器在多个定点实时按秒记录这些指标的数据,期望通过对记录数据的分析研究找到地震的前兆特征。

地震预测问题数学建模论文

2015年北京交通大学数学建模竞赛 承诺书 我们仔细阅读了数学建模竞赛的竞赛规则。 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 所属学校(请填写完整的全名): 参赛队员(打印并签名) :1. 2. 3. 日期:年月日赛区评阅编号(由赛区组委会评阅前进行编号):

目录 一、摘要 (2) 二、问题的重述 (3) 三、模型的假设与符号说明 (3) (1)、模型的假设 (3) (2)、符号说明 (4) 四、模型的建立 (4) (1). 由原始数据计算一次累加数据 (4) (2). 建立矩阵 (4) (3). 求u和a (2) (4). 时间响应方程 (2) 五、模型的求解 (3) (1). 矩阵求解 (3) (2).求u和a (3) (3). 时间响应方程 (4) 六、模型检验或误差分析 (7) (1)、残差检验: (7) (2)、后验差检验: (8) (3)、结论 (8) (4)、误差分析 (8) 七、模型评价 (8) 八、参考文献 (9)

地震预报问题 一、摘要 本文通过建立模型,解决了地震预报问题,简述如何利用正确的建模,从几组数据中预测下一次地震的时间与地点。 20世纪以来发生在中国8级以上大地震共四次,数据较少,很难找到规律,用神经网络等建模形式来“训练数据”较为困难,且到目前为止,地震的发生尚无规律可言,即使是相邻的两次地震,彼此之间也很难说有什么影响。对于地震来说,筛选出的四个数据,它们的时间地点没有必然的联系,所以可以将时间和地点的预测分开进行建模。 在时间上,由于数据较少,内部数据规律不能很好地表示出来,可以选择灰色预测,它是一种对含有不确定因素的系统进行预测的方法,可以弥补地震事件中的不确定性。通过关联分析几次地震时间,并对原始数据进行生成处理来寻找系统变动的规律,生成有较强规律性的数据序列,然后建立相应的微分方程模型,从而预测事物未来发展趋势的状况,预测未来某一时刻的特征量,或达到某一特征量的时间。最终运用matlab进行编程得出结果。 在地点上,根据板块漂移假说和地震带理论,发生在地震带上的可能性较高。对于地点,主要通过经纬度来确定,假设经度和纬度之间没有关系,通过分别建模预测下次地震的经度和纬度。由于数据较少,可以运用曲线拟合法,通过图像的类比得到近似的函数。最小化问题的精度,依赖于所选择的函数模型,用q个相关变量或p个附加的相关变量去拟和。通常人们将一个可能的、对不相关变量的构成都无困难的函数类型充作函数模型(本题中运用正弦函数)。求出参数后,用参数使所选择的函数模型同观测值Y 相匹配。最终得到可预测下次地震地点的函数。在求解方面,利用了mathematica软件求解。 经过建模求解,并进行误差分析之后,预测下次8级及以上的地震时间是2087年10月15日,地点是N36.55?,E106.76?(甘肃省内) 论文的末尾给出了模型优缺点的分析和评价,并提出了改进方向,如果数据再多一些,或者是要进行7级以上地震预测的话,会更准确一些。 关键词:地震预报时间地点灰色理论模型非线性回归模型

地震波数值模拟方法研究综述

地震波数值模拟方法研究综述 在地学领域,对于许多地球物理问题,人们已经得到了它应遵循的基本方程(常微分方程或偏微分方程)和相应的定解条件,但能用解析方法求得精确解的只是少数方程性质比较简单,且几何形状相当规则的问题。对于大多数问题,由于方程的非线性性质,或由于求解区域的几何形状比较复杂,贝y 不能得到解析解。这类问题的解决通常有两种途径。一是引入简化假设,将方程和几何边界简化为能够处理的情况,从而得到问题在简化状态下的解答。但这种方法只是在有限的 情况下是可行的,过多的简化可能导致很大的误差甚至错误的解答。因此人们多年来寻找和发展了另一种求解方法一- 数值模拟方法。 地震数值模拟(SeismicNumericalModeling) 是地震勘探和地震学的基础,同时也是地震反演的基础。所谓地震数值模拟,就是在假定地下介质结构模型和相应的物理参数已知的情况下,模拟研究地震波在地下各种介质中的传播规律,并计算在地面或地下各观测点所观测到的数值地震记录的一种地震模拟方法。地震波场数值模拟是研究复杂地区地震资料采集、处理和解释的有效辅助手段,这种地震数值模拟方法已经在地震勘探和天然地震领域中得到广泛应用。 地震数值模拟的发展非常迅速,现在已经有各种各样的 地震数值模拟方法在地震勘探和地震学中得到广泛而有效

的应用。这些地震波场数值模拟方法可以归纳为三大类,即几何射线法、积分方程法和波动方程法。波动方程数值模拟方法实质上是求解地震波动方程,因此模拟的地震波场包含 了地震波传播的所有信息,但其计算速度相对于几何射线法要慢。几何射线法也就是射线追踪法,属于几何地震学方法,由于它将地震波波动理论简化为射线理论,主要考虑的是地震波传播的运动学特征,缺少地震波的动力学信息,因此该方法计算速度快。因为波动方程模拟包含了丰富的波动信息,为研究地震波的传播机理和复杂地层的解释提供了更多的佐证,所以波动方程数值模拟方法一直在地震模拟中占有重要地位。 1地震波数值模拟的理论基础 地震波数值模拟是在已知地下介质结构的情况下,研究地震波在地下各种介质中传播规律的一种地震模拟方法,其理论基础就是表征地震波在地下各种介质中传播的地震波传播理论。上述三类地震波数值模拟方法相应的地震波传播理论的数学物理表达方式不尽相同。射线追踪法是建立在以射线 理论为基础的波动方程高频近似理论基础上的,其数学表形式为程函方程和传输方程。积分方程法是建立在以惠更斯原理为基础的波叠加原理基础上的,其数学表达形式为波动方程的格林函数域积分方程表达式和边界积分方程表达式动方程数值解法是建立在以弹性或粘弹性理论和牛顿力学

地震预测模型

龙源期刊网 https://www.360docs.net/doc/9414397488.html, 地震预测模型 作者:宫唐小恒林建明苏静周杰 来源:《电脑知识与技术·学术交流》2008年第32期 摘要:文章旨在研究地震预测问题,讨论了全球16个强震活动区的强震活动次数与其后一年我国强震发生与否的关系。首先,由模糊聚类和变量规格化的思想,结合前15年的数据求出全球16个区发生强震次数与中国后一年发生强震的相关系数矩阵,通过对前15年16个区发生强震次数加权求和,得到各地的影响值,取其中的最小值作为影响值的下限值,建立预测预测中国后一年是否发生强震模型;其次,从数据中筛选出云南地区每月最高震级组成固定时间间隔的地震序列,利用气象学预测中的均生函数方法生成一列相应的时间序列,经周期性外延、残差值逐步拟合得到对当地地震最大震级的预测模型。 关键词:地震预测;相关系数;加权和;均生函数;时间序列 中图分类号:P31文献标识码:A文章编号:1009-3044(2008)32-1252-05 The Model of Earthquake Prediction GONG TANG Xiao-heng1, LIN Jian-ming2, SU Jing3, ZHOU Jie2 (1.School of Computer Science, Sichuan University,Chengdu 610065,China;2.College of Mathematics, Sichuan University,Chengdu 610065,China;3.School of Electrical Engineering and Information,Sichuan University,Chengdu 610065,China) Abstract: This paper was mainly about the earthquake prediction on the basis of the relationship between the frequency of strong earthquakes in 16 seismically active zones all over the world and the probability of strong earthquakes which may happen in China the next year. Firstly, according to fuzzy clustering and variable normalizing, we obtained a relevant coefficient matrix about the relationship between the frequency of strong earthquakes in 16 seismically active zones all over the world and the probability of strong earthquakes which may happen in China the next year with data of the previous 15 years. Through weighted summation of the frequency of strong earthquakes in 16 seismically active zones all over the world during the previous 15 y ears, we obtained the “influencing value” of each zone and set the minimum value as the lower limit value, and then established the prediction model of strong earthquakes in China the next year. Secondly, we screened out the highest earthquake magnitudes each month in Yunnan to compose an earthquake sequence. Next we applied to mean generating function from Meteorological prediction to generate a time sequence. Thus we established the prediction model of the highest earthquake magnitudes in Yunnan by using cyclic extension and adjusted residual gradual fitting. Key words: earthquake prediction; relevant coefficient; weighted summation; mean generating function; time sequence

相关文档
最新文档