按：转帖不代表本人赞同或反对文章观点。 他们预测了08年512大地震“国内要是不要,我们就给国外。”两位爱国者长长的叹息 我相信,我的民族能接受一切真相。——笔者翟明磊 中国地震局发言人称: 地震不可预测,这次汶川地震情况尤为特殊,震前没有发现任何前兆信息,没有前震,没有收到任何预测意见。我历时四个月的调查证明这是彻头彻尾的谎言。 但本文并不满足于证伪，我想带给读者的不是一次寻找小丑与英雄的旅行，而是面对科学与我们人性的惨痛见证，汶川地震预报失败在我们民族性格与体制中的腐败点：不合作，不独立，不负责，胆怯。我对自己的要求是不哗众取宠，不偏不倚，恰如其分,我对读者唯一的要求：对科学问题更多耐心与理解。 让我们开始寻找汶川预测真相的旅行。 ——笔者翟明磊 半个台测出汶川地震 2008年5月11晚上，北京民族学院南路一间屋子的灯通宵未熄，一对白发苍苍的老夫妻紧紧盯着计算机屏幕。老先生的老花镜只有一条腿，歪歪地架在鼻梁上有点滑稽。——他们是全世界唯一知道地震马上就要发生的人。 凌晨两点时，他们知道的是：距四川红格550公里到800公里的环带中,12日至13日将发生七级到八级的大地震。 他们不知道的是：震中的方位与具体地址。 因为他们已经弹尽粮绝，他们只剩下半个台站。
钱复业不为公众所知，这位七十四岁的老人在地震界却是一位实力人物。人称中国洋地电的“祖师奶奶”。当年她还是一位俊俏姑娘时就奉周总理之命，将苏联的地电方法引入中国，这位当年的风火的巾帼标兵在邢台试验场一干就八年。提起当年，最难忘的是邢台百姓拔了他们架下的电线,三十岁钱复业说“我们是为人民服务的。”邢台人对没报出地震的钱说：“你们为人民屁服务。”八年没有休息日，没有回家，孩子管父亲叫“叔叔。”活得象野人一样，常年住帐篷吃干粮，15分钟观测一次仪器.在特大洪水中是老百姓救了她们。地电的方法是将一公里长的电线埋入地下，当地层受压时，地下三百立方米的体积电阻变大，电阻率变小，从而测出地震前兆。得到李四光首肯后，这一方法在全国推广，成为预测地震的主要手段。目前全国仍有110个地电台，负责人大多是钱复业的徒子徒孙。 2003年，“地电祖师奶奶”闹革命了。因为她发现全国的地电台站大部分失灵了。原因是自从全国花了二十七个亿更新数字化台站后，这位创始人发现自己报不出地震了。这一年，她上报了四次地震，三次成功，一次失败。前三次用的是手动的传统地电仪，虚报的一次用的是地震局数字台站的数字仪器。 图片说明下图为汶川地震HRT波前兆红色部分即为异常区域。上图为唐山松潘等三个地震异常比较
地震资料反演技术概论(波阻抗、岩性反演处理技术) 一九九八年九月 辽河油藏工程培训班材料 编写人：钟俊
地震资料反演技术概论 前言 一.反演的概念、目的 二.反演的发展历史及趋势 三.反演的基本方法 四.反演的限制条件 五.反演的基本流程 六.反演实例
前言 地震、测井、钻井是石油工作者认识地下地质构造、地层、岩性、物性、含油气性的最重要的信息来源。虽然测井、钻井仅能提供井孔附近的有关信息，尤其是有关岩性、物性、含油气性的信息，但是这些信息往往具有很高的分辨率，可信度、准确性，能确切地指出含油气层的位置，定量化分析与储层、油藏有关的参数。然而一个油气田勘探、开发方案的设计、实施、调整仅靠测井、钻井资料是远远不够的，必须与地震资料相结合进行综合分析才能取得良好效果。 地震资料的分辨率虽然远远不及测井、钻井，但是随着地震勘探技术的发展，从光电记录、模拟记录到数字记录，从二维到三维，地震资料的信噪比、分辨率、成像的准确性都获得了极大的提高，由于地震资料包含大量地下地质信息，覆盖面积广，具有三维特性，所以这项技术的使用越来越受到石油工作者的重视，如何利用地震资料研究地下地质构造、地层？如何进行储层预测、油藏描述？如何进行油藏、含油气层的预测？ 这些问题促使地球物理学家、地质学家开发应用了一系列地震资料特殊处理技术，如地震资料反演技术、地震属性分析技术、AVO分析技术，这些技术充分利用测井、钻井、地震的长处，使人们对地下储层、油藏的研究从点到面、从二维到三维、从三维可视化研究到油藏动态监测、从定性研究到定量化研究，大大提高了钻探成功率，有效地指导了油田开发，
我国未来地震预报和防灾工作开展的几点建议 摘要：四川汶川大地震对我国人民生命财产造成了巨大损失，通过分析该次地震的突然性和成因以及地震发生后中暴露出的建筑物大量倒塌，较高的震亡率等问题等对我国地震预防和防灾机制进行了分析，并提出了个人的一些见解 关键词：地震预防抗震减灾城市生命线 引言 我国由于东邻环太平洋地震带，南接欧亚地震带，地震情况相当复杂，具有频度高、强度大、分布广、震源浅的特征，同时由于人口稠密、建筑物抗震能力较低，因此我国是世界上地震活动水平最高、地震灾害最重的国家，震后建筑物和城市生命线遭到严重破坏，火，水，气，泥石流等次生灾害也时有发生，人民的安全受到严重威胁，因此对于地震的预防和采取相应的防震减灾措施已经成为当务之急。 2008年5月12日在我国四川省汶川发生的8.0级大地震属于大陆内部地震，是大型板块的断层活动，为浅源地震，由北川断层的逆冲-右旋错动导致的。同时由于聚集能量巨大，在突然释放时能量沿着板块裂缝传递，对各板块产生挤压，地震破裂尺度较大，导致其他省份也产生明显震感，最终造成数万人死亡，建筑物大量倒塌，同时引发包括堰塞湖，泥石流等次生灾害。汶川大地震也是中国一九四九年以来破坏性最强、波及范围最大的一次地震，在断层错动时间，地震张量指数，以及地震的强度、烈度上均超过了1976年唐山大地震。 地震发生前夕，我国地震监测台网没有做出地震预报，但由于我国防灾应急机制启动迅速，以及各省市对受灾地区的迅速有效支援最大限度的减小了人员伤亡。 1 地震预报的现状和困难 严重的地震灾害及其加速发展的震灾形势给社会公众带来了对防震减灾的强烈社会需求。而通过地震预报和在预报基础上的震灾防御是实现地震减灾的最基本途径。 1.1 地震预报的研究现状 从６０年代中期开始，世界一些地震频繁的国家相继开展有计划的地震预报研究，美国和日本都是多地震国家，也是目前世界上地震预报技术最先进的两个国家，其中美国于1964年组织了一批有声望的地震科学家拟定了地震预报的研究规划，开展了与地震孕育、发生相关的地震活断层调查、地震前兆观测和地震孕育理论等地震预报研究，并于８０年代在加利福尼亚州一个名叫帕克菲尔德的地震区建立了地震预报实验场。日本政府从1964年开始推行地震预报研究的第一个５年计划。1994年已进入第７个地震预报５年计划，其重点是地震预报实用化和确定地震预报方法、提高地震预报精度的观测研究，并加强地震预报的基础研究和新技术开发。前苏联则从６０年代初开始，在中亚远东地区建立一系列地震预报实验场，开展地震预报的现场研究和基础性的实验论研究。但从总体上看，３０多年的科学进展与实现地震预报的科学目标之间还存在很大的距离。正如美国地震学会会长、地震预报评估委员会主席、加州理工学院教授克拉伦斯．艾伦在评定地震预报进展情况时所说：地震预报的进展要比初期预料的缓慢得多，地震预报的科学难度要比原先预料的困难得多。 1.2 困难所在 地震预报由于涉及到大陆地震成因和孕震理论，地震前兆机理，地震前兆探测中的基础性研究问题等因素。考虑到地球的不可入性，以及在不同的地理构造环境、不同的时间阶段，不同震级的地震显示出的地震孕律的复杂性和地震发生的小概率性使得当前地震预报工作整体进展不快。当前的地震预报总体水平是很低的。准确的短临预报意见也是非常少的，因此目前世界地震学界最主流的学术观点仍然是地震无法预测。
竭诚为您提供优质的服务，优质的文档，谢谢阅读/双击去除 地震能否预测 地震能否预测 其实，“地震能否预测”作为一个基本的科学问题，一直拷问着地球科学界。 美国密苏里州立大学教授刘勉认为，要回答这一问题，需要将长期预测和短期预测分开来谈。长期预测目前已经有了很多进展。但对于短期预测来说，人类仍有很长的路要走。当然，一些新技术的出现，比如全球定位系统，会为地震预测研究带来很大的帮助。 伊利诺伊州立大学的副教授宋晓东（音译）的回答是：“依据人们目前对地震发生过程的了解程度，还不足以得出能否进行短期预测的结论。中国是世界上唯一成功进行过地震预报的国家，但这属于经验型预测范畴，其中包括年成功预测海城大地震。” 日本ｅｈｉｍｅ大学的赵大鹏教授的答案更为保守；“根据我们目前的知识和能力，还不可能够预测地震。只有在人们不断加深对地震的了解之后，才可能得出能否预测地震的明确结论。” 多年来，对这一问题国际上一直存在着激烈争论。其中，一个很
有影响的事例是，ｋｏｇａｎ、ｇｅｌｌｅｒ、ｊａｃｋｓｏｎ三位教授于上世纪九十年代末联名在《科学》杂志上发表了《地震无法预测》的论文。他们经历了近年苦心研究之后认为，地震是无法预测的，应打消可能会在几小时、几天或几个月之前预测地震的念头。他们的主要论点是：总结的地震前兆都不可靠；现行的经验预报方法是走不通的；按照地震过程的物理实质，准确预报是不可能的。他们认为“这是一项毫无希望的工作！” 在中国地震局的官方网页上，则有这样的介绍：“我们做出的较大时间尺度的中长期预报已有一定的可信度，但短临预报的成功率还相对较低……从世界范围说，地震预报仍处于探索阶段，尚未完全掌握地震孕育发展的规律，我们的预报主要是根据多年积累的观测资料和震例，进行经验性预报。” 进行预测有无必要？ 在《简明大不列颠百科全书》上，有对地震这样的描述：“任何天灾都比不上地震，能在如此短促的时间，如此广大的范围，造成如此巨大的损失。”正因如此，许多人也一直试图预测地震的发生。 其中一类就是一些痴迷的“民间爱好者”。这些人没有接受过基本的科学教育，往往通过一些神奇的手段，就自认为成功预测了地震的发生。不久前，一位河北的农民曾找到记者说，他自己早就预测到了伊朗巴姆大地震将要发生，其唯一的证据就是他手写的一堆谁也看不懂的数学公式。 在，一个经营私人天文台的日本人串田嘉男一夜成名。据报道，
No.13,2010 现代商贸工业 Modern Bus iness Trade Industry2010年第13期 地震波层析成像反演方法及其研究综述 冯 微 (长江大学物理科学与技术学院,湖北荆州434025) 摘 要:通过研究利用初至波走时的层析反演方法建立近地表速度模型,提供近地表地下介质的速度信息,进一步为静校正或浅层工程勘探服务。 关键词:速度建模;层析成像;初至波 中图分类号:TB 文献标识码:A 文章编号:1672 3198(2010)13 0368 01 地震勘探是利用人工在地表激发和接收地震波,再对地震波作分析处理以及解释而得到地下构造信息和岩性信息的一种方法。在整个地震勘探过程中,精确的求取地震波在地下介质中的传播速度,一直是地震勘探的核心问题之一。尤其在地表条件较复杂的区域,地表速度的横向剧烈变化会严重影响中深层目的层的成像效果。近地表速度不准确,将会直接影响到速度分析、偏移成像的质量以及静校正的精度等地震勘探的各个环节和最终的勘探成果。 1 地震面波及波形反演 利用面波进行结构反演一直是了解地球介质结构的重要途径。近几年来,在面波理论和面波反演方面做了大量工作。陈蔚天和陈晓非(2001)提出了一种求解水平层状海洋-地球模型中面波振型问题的新算法,它简洁、高效,彻底消除了高频情况下数值计算的精度失真问题。张碧星等(2000,2002)对瑞利波勘探中 之字形频散曲线形成的物理机理和多模性问题进行了理论分析,研究了诸波模的传播特性及相互关系,以及地表下低速层介质的位置、厚度及其它参数对 之字形频散曲线的相互影响.在面波反演理论方面,朱良保等(2001)通过保角变换,把面波群速度的反演变成了球谐系数的线性化反演,使其计算速度快,等值线光滑,构造界限清晰。众多研究者根据从面波资料求出的频散曲线,对不同地区的地下速度结构作了反演,揭示了横向结构差异的广泛存在。 根据走时反演地下结构是获取结构信息的经典做法。刘伊克等(2001)根据三维地震观测的初至走时数据,利用最小平方与QR分解相结合的算法,在三维空间重建近地表低降速带速度模型。同时,采用分形算法克服了初至波波形差异以及折射波相位反转导致的拾取误差,实现了三维初至拾取的大规模全自动化运算。李录明等(2000)针对地震勘探中的复杂地表问题,提出了一套地震初至波表层模型层析反演方法.它利用地震直达波、回折波、折射波以及三者组合的初至波和层析反演方法具有的纵、横向变速优势,实现适应速度任意变化的复杂表层模型反演。 在利用远震体波接收函数反演地下结构方面。钱辉等(2001)对接收函数反演地壳结构速度的算法作了分析,使之适应正演参数的变化,并利用天然地震接收函数揭示了青藏高原东部地壳结构。 近年来,非线性反演越来越受到重视,许多研究者把新的最优化理论引入地震学反演中。孟洪鹰和刘贵忠(1999)提出了多尺度地震波形反演的小波变换方法。对于一维非线性地震波形反演问题,此方法和已有的简单迭代法及多重网格法比较表明,此方法更为有效。杨峰和聂在平(2000)提出了用于二维轴对称非均匀介质结构的反演和成像的一种新的反演迭代方法变分玻恩迭代方法.与传统的玻恩迭代方法相比,其收敛速度和成像质量均有较大改善。 2 地震勘探、测井问题中的地震波研究及其它 在地震勘探和测井方面,许多研究者针对实际问题,提出了新的方法。沈建国和张海澜(2000)计算了井内靠近井壁的偏心声源激发的声场,得到了在井壁不同位置的接收波形,分析了直达波、井壁反射波、纵波、横波和面波在这些波形中的反映。为了处理横向强变速介质中的深度成像问题,程玖兵等(2001)提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,在基于反射系数估算的成像条件下,可实现叠前深度偏移成像。陈生昌等(2001)实现了一种基于拟线性Born近似的叠张海明等:地震波研究前深度偏移方法,扩大了拟线性Born近似的应用范围,使其能够适应更强的横向速度变化。张美根和王妙月(2001)利用有限元法和最小走时射线追踪的界面点法,实现了各向异性弹性波的叠前逆时偏移.陈志德等(2002)利用叠前深度域地震成像对速度模型变化的敏感性,采用偏移迭代逐次逼近最佳成像速度,研究开发了一套快捷有效的三维叠前深度偏移深度域速度模型建立技术。顾汉明等(2002)在频率-波数域中采用解析法,解出多层条件下海底实测的多分量地震数据分解成上行和下行P波和S波的算法,导出海底各层地震反射系数随入射角变化(简称RVA)的递推计算公式。金胜汶等(2002)给出了一种高效率、高精度的炮检距域叠前深度偏移方法,并得到各个不同照射角下的成像结果。 3 讨论和结论 地震波理论是固体地球物理学研究的重要基础.地震波研究领域的任何实质性进展都会促进固体地球物理学的发展.在过去的4年里,中国地球物理学家在该领域做了很多有意义的研究工作,其中不乏创新性的理论工作.当前地震波研究领域的重要课题包括: (1)复杂地球介质中地震波激发与传播理论; (2)高效计算三维介质中地震波传播的数值方法; (3)利用先进的地震波数值模拟方法,开展设定地震与强地面运动的数值模拟研究,为精细的地震危险分析与预测奠定基础。 参考文献 [1]周庆凡.我国天然气发展前景广阔[J].中国石化,2009. [2]刘英祥.我国天然气价格与天然气发展问题研究[J].企业经济, 2009. [3]牛建娣.我国天然气市场供需状况及发展对策分析[D].对外经济 贸易大学,2007. ! 368 !
精心整理2011年赣南师院数学建模竞赛选拔赛 题目地震预测模型 摘要： 本文前三个任务主要考虑是各指标的变化对地震发生问题的影响，通过对各指标数据量的分析建立相应的模型，并对任务四和任务五给出了合理的解答。 针对任务一：我们从原始数据中计算出各项指标的日均值，绘制出各指标分年度的时间序列图， 磁波幅度 。 关键词： 一·问题的重述 1.1背景分析 地震是地壳快速释放能量过程中造成的振动。虽然预测地震是世界性难题，但迄今科学界普遍认为，有可能反映地震前兆特征的指标可能不少于10个。已经有专业仪器在多个定点实时按秒记录这些指标的数据，期望通过对记录数据的分析研究找到地震的前兆特征。 现已采集到某地2005年1月1日至2010年6月30日按小时观测的10多个指标的数据，和该地区该时期内已发生地震的时刻、经纬度、震级及震源深度的数据。这些数据中隐藏着地震发生的前兆特征。科学地截取这些数据的有用片段，对数据进行合理地预处理，用数学方法揭示地震前兆
的数据特征，是一项很有意义的研究工作。 题给数据中的这10多个指标，究竟哪些与地震的发生有关，有何种关系，是单一关系还是复合关系；除这10多个指标外还有哪些因素及含题给指标在内的哪些指标的哪种数学模型更能反映地震的前兆特征等等，人们迄今仍不很清楚，需要进行深入地研究。地震数据的观测是持续进行的，随着时间的推移数据的规模会不断扩大。从中挖掘地震的前兆特征，必须有合理的数学模型，也必须有科学高效的算法分析平台。因此，需要我们结合附件中给出的实际记录数据，尝试完成以下任务。 1.2任务的提出 任务一：分析数据特征，建立数学模型以度量各指标对地震发生的敏感程度。 越大 任务三：中要结合题给数据，建立数学模型来研究地震发生前的数量特征。主要运用贝叶斯判别分析法进行建模，对已给数据进行先验信息、后验信息分析。 任务四：要将计算程序集结成地震数据分析平台，能够完成其它地震数据的分析，并能自动输出前任务的重要分析结果。 任务五：是针对进一步的研究设想写一篇切实可行的报告。 三·问题的基本假设 （1）地震监测点的监测设施能正常运转； （2）地震监测设施周围不存在影响其工作效能的干扰源，如飞机场、发电厂等；
地震波阻抗反演方法综述 一、地震反演技术研究现状 地震反演方法是一门综合运用数学、物理、计算机科学等学科发展起来的新技术新方法，每当数学方法、物理理论有了新的认识和发展时，就会有新的地震反演技术、方法的提出。随着计算机技术的不断发展、硬件设施的不断升级，这些方法技术得到了实践验证和提升，反过来地震反演技术运用中出现的新问题、新思路又不断促使数学方法、地球物理学理论的再次发展。时至今日，地震反演技术仍然是一个不断发展、不断成熟、不断丰富着的领域。 反演是正演的逆过程，在地震勘探中正演是已知地下的地质构造情况、岩性物性分布情况，根据地震波传播规律和适当的数学计算方法模拟地震波在地下传播以及接收地震波传输到地表信息的过程。地球物理反演就是使用已知的地震波传播规律和计算方法，将地表接收到的地震数据通过逆向运算，预测地下构造情况、岩性物性分布情况的过程。地震波阻抗正演是对反演的理论基础和实现手段。 1959年美国人Edwin Laurentine Drake在宾夕法尼亚州开凿的第一口钻井揭开了世界石油工业的序幕。从刚开始的查看地质露头、寻找构造高点寻找石油，到通过地震剖面的亮点技术寻找石油，再到现在运用多种科学技术手段进行油气资源的预测，石油勘探经历了一个飞速的发展历程。 声波阻抗（AI）是介质密度和波在介质中传播速度的乘积，它能够反映地下地质的岩性信息。声波阻抗反演技术是20世纪70年代加拿大Roy Lindseth博士提出的，通过反演能够将反映地层界面信息的地震数据变为反映岩性变化的波阻抗（或速度）信息。由于波阻抗与地下岩石的密度、速度等信息紧密联系，又可以直接与已知地质、钻井测井信息对比，因此广泛应用于储层的预测和油藏描述中，深受石油工作者的喜爱。70年代后期，从地震道提取声波资料的合成声波技术得到了快速发展，以此为基础发展的基于模型的一维有井波阻抗反演技术，提高了反演结果的可靠性。进入80年代，Cooke等人将数学中的广义线性方法运用于地震资料反演，提出了广义线性地震反演。此后Seymour等人又提出了测井声波资料和地震数据正反演相结合求取地下声波阻抗的测井约束反演，大大拓宽了反演结果的纵向分辨能力。 90年代，在基于前人对地质统计学研究的基础上Bortoli和Haas提出了地质统计学反演，Dubrule等人对该方法进行了改进和推广。在国内随着油田对地震反演技术的广泛应用，以周竹生为主提出的地震、地质和测井资料联合反演方法，将地质信息引入地震反演中，提高的反演结果与地质认识的联系，克服了线性反演存在的缺陷。1996年，李宏兵等人将宽频带约束方法应用于递推反演并对其进行改进，减弱了噪音对反演结果的影响。 1999年，任职于英国石油公司的Connolly在《弹性波阻抗》一文中介绍了弹性波阻抗（EI）的概念和计算方法，阐述了不同入射角度（偏移距）地震道集部分叠加反演波阻抗随入射角之间的关系，但是该方法求取的弹性阻抗随入射角变化很大，无法与常规叠后反演波阻抗直接比较，因此推广应用较为困难。2002年，Whitcombe通过修正Patrick Connolly的计算公式，得到了弹性波阻抗的归一化求取方法，消除了弹性阻抗随入射角变化大的难题。2003年，西北大学马劲风教授从Zoeppritz方程简化出发提出了广义弹性波阻抗的概念，克服了以往波阻抗反演要求地震波垂直入射到地表的假设条件，推导出了任意入射角下纵波反射系数的递推公式，提高了中等入射角度下弹性波阻抗反演的精度。
【全文】地震属性分析技术综述 [摘要] 地震属性是从地震资料中提取的隐藏有用信息，因而地震属性分析技术近几年在油气勘探开发中得到了广泛的应用与研究。本文对地震属性分析技术的发展状况进行了归纳、总结，简单阐述了地震属性分析技术的在不同时期所用到的基本原理和方法。特别对新地震属性进行了具体介绍。最后对该技术进一步的研究工作进行了总结和展望。 摘要:在勘探和开发周期的各个阶段，地震资料在复杂油藏系统的解释过程中，扮演着至关重要的角色。然而，缺少一种有效地将地质知识应用于地震解释中的上具。随着一系列属性新技术的出现，对地震属性进行充分研究，就给地质家提供了快速地从三维地震数据中获得地质信息的能力。尤其在用常规解释手段难以识别日的储层的情况下，属性分析技术更是给地质上作人员指出了新的方向。 [关键词] 地震属性储层预测叠前数据叠后数据 关键词:储层;波形分析;地震属性 1.引言 地震属性是指叠前或叠后的地震数据经过数学变换而导出的有关地震波的几何形态、运动学特征、动力学特征和统计学特征的特殊度量值。地震属性的发展大致从20世纪60年代的直接烃类检测和亮点、暗点、平点技术开始，经历了70年代的瞬时属性(主要是振幅属性)和复数道分析，90年代的多维属性(特别是相干体属性)分析，21世纪的地震相分析等阶段[1一SJ。随着地震属性分析技术的发展与研究，该技术已广泛应用于储层预测、油气藏动态监测、油气藏特征描述等领域，并取得了很好的效果。总之，地震属性分析技术可以从地震资料中提取隐藏其中的多种有用信息，这为油气勘探与开发提供了丰富宝贵的资料，也为解决复杂地质体评价提供了实用的分析手段。因此，对该技术进行深人调查研究具有很强的现实意义。 地震属性是指从地震数据中导出的关于儿何学、运动学、动力学及统计特性的特殊度量值。它可包括时问属性、振幅属性、频率属性和吸收衰减属性，不同的属性可指示不同的地质现象。地震属性分析则是从地震资料中提取其中的有用信息，并结合钻井资料，从不同角度分析各种地震信息在纵向和横向上的变化，以揭示出原始地震剖面中不易被发现的地质异常现象及含油气情况。 地震属性分析技术的研究已由线、面信息扩展到三维体信息，从分类提取扰化发展为一项系统的应用技术。随着地震技术的日趋成熟，地震属性技术近儿年也发展迅速，其中有多属性联合解释技术、波形分析技术、吸收滤波技术等。应用地震属性分析技术去完善勘探生产中的油藏描述工作，已经成为油藏地球物理的核心内容。利用地震属性分析技术预测岩性和有利储集体，描述油藏特征及孔隙度变化，寻找难以发现的隐蔽油区，以至于监测流体运动和进行其它综合研究，一直是石油工作人员追求的目标。 1波形分析技术的研究与应用 通常的层段属性只是表示了某儿个地震信号的物理参数(振幅、相位、频率等)，但它们没有一个能够单独描述地震信号的异常，而地震信号的任何物理参数的变化总是对应着反映地震道形状的变化，所以，研究和分析地震资料中代表各种属性总体特征的地震道形状(波形)，应该能有非常不错的效果[，]。 1. 1波形分析技术的原理及处理过程
地震预测 根据所认识的地震发生规律，用科学方法对未来地震发生的时间、地点和强度作预先的估计。地震预报则是在具备一定可靠程度的前提下将地震预测的意见向公众宣布。有实用价值的地震预报必须同时报出时间、地点和强度。科学的地震预测是将来实现地震预报的基础。 地震预测是第二次世界大战结束以后开展的探索性研究项目，特别是中、短期或临震前的预测尚处于探索阶段，远远没有到可以实用的程度。一些学者对实现地震预报抱有怀疑，对于用行政手段组织地震预测持保留态度。 地震预测的科学前提是认识地震孕育和发生的物理过程，包括地球介质物理、力学性质的异常变化。但是人类对地震成因和地震发生的规律还知之甚少，主要是因为地震是宏观自然界中大规模的深层的变动过程，不同于实验室中单纯的可控条件下进行的样品试验，其影响因素过于复杂，还可能有人类未知的因素存在。人们所能做的是在地面上观测某些物理量，这种观测通常是不完全的和不完善的。在当代科学技术条件下，人们还不能深入地球内部直接或间接观测深层介质的物理状态，因为测量过程本身就将打破原有的状态。而所能观测的物理量异常变化是否与地震的发生真正相关还不能确知。这就是地震预测研究所以进展缓慢的真实原因。 地震预测研究有3种不同的思路: ①地震地质。地震发生在地壳中上层，故认定地震应属于地质过程。研究已发生的大地震的地质构造特点，应有助于今后判定何处具备发生大地震的地质背景。但有些地震发生前，地质构造往往不甚明朗，震后才发现有某个断层，认为与地震有关。 ②地震统计。对过去已发生的地震，运用数理统计方法，从中发现地震发生的规律，特别是时间序列的规律，根据过去以推测未来。此法把地震问题归结为数学问题。因需要对大量地震资料作统计，研究的区域往往过大，所以判定地震的地点有困难，而且外推常常不准确。 ③地震前兆。地震是地球介质的破裂，故认定地震应属于物理过程。观测地球物理场各种参量以及地下水等异常变化，可能找到有用的地震前兆。前兆研究中的最大困难是，观测中常遇到各种天然的和人为的干扰，而所谓的前兆与地震的对应往往也是经验性的。尚未找到一种普遍适用的可靠的前兆。 以上3种思路都有片面性，都不能独立地解决地震预测问题。实际采取的是综合的办法，把3种不同思路所得放在一起对比参照，力求对未来的地震活动作出估计。 地震预测是世界难题，第一，地球的不可入性。大家知道上天容易入地难，我们对地下发生的变化，只能通过地表的观测来推测；第二，地震孕律的复杂性。通过专家多年的研究，现在逐渐认识到地震孕育、发生、发展的过程十分复杂，在不同的地理构造环境、不同的时间阶段，不同震级的地震都显示出相当复杂的孕律过程；第三，地震发生的小概率性。大家可能都感觉到，全球每年都有地震发生，有些还是比较大的地震。但是对于一个地区来说，地震发生的重复性时间是很长的，几十年、几百年、
第一章反演理论 第一节基本概念 一．反演和正演 1．反演 反演是一个很广的概念，根据地震波场、地球自由振荡、交变电磁场、重力场以及热学等地球物理观测数据去推测地球内部的结构形态及物质成分，来定量计算各种有关的物理参数，这些都可以归结为反演问题。在地震勘探中，反演的一个重要应用就是由地震记录得到波阻抗。 有反演，还有正演。要正确理解反演问题，还要知道正演的概念。 2．正演 正演和反演相反，它是对一个假设的地质模型，给定某些参数（如速度、层数、厚度）用理论关系式（数学模型）推导出某种可测量的量（如地震波）。在地震勘探中，正演的一个重要应用就是制作合成地震记录。 3．例子 考虑地球内部的温度分布，假定地球内部的温度随深度线性增加，其关系式可表示成：T(z)=a+bz 正演：给定a和b，求不同深度z的对应温度T(z) 反演：已经在不同点z测得T(z)，求a和b。 二．反演问题描述和公式表达的几个重要问题 1．应用哪种参数化方式——离散的还是连续的？ 2．地球物理数据的性质是什么？观测中的误差是什么？ 3．问题能不能作为数学问题提出，如果能够，它是不是适定的？ 4．对问题有无物理约束？ 5．能获得什么类型的解，达到什么精度？要求得到近似解、解的范围、还是精确解？ 6．问题是线性的还是非线性的？ 7．问题是欠定的、超定的、还是适定的？ 8．什么是问题的最好解法？ 9．解的置信界限是什么？能否用其它方法来评价？ 第二节反演的数学基础
一．解超定线性反问题 1．简单线性回归 可利用最小平方法确定参数a 、b 使误差的平方和最小。 ??? ? ?? ?∑-∑∑∑-∑= -=∑∑-=2 2)()(x x n y x xy n b x b y n x b y a （1-2-1） 拟合公式为： bx a y +=? （1-2-2） 该方法的公式原来只适用于解超定问题，但同样适用于欠定问题，当我们有多个参数时，称为多元回归，在地球物理领域广泛采用这种方法。此过程用矩阵形式表示，则称为广义最小平方法矩阵方演。 2．非约束最小平方法反演——广义矩阵方法 由前面讨论可知，参数估计的最小平方方法用矩阵公式表示，所得到的算法等价于一个或多个模型参数的一个或多个数据集反演，步骤为： 问题定义→矩阵公式→最小平方解 线性问题采用广义矩阵形式 d=Gm （1-2-3） 对于精确的数据模型，参数m 为 m=G -1d （1-2-4） 但是由于试验误差，实际数据将不能精确拟合获得，故采用最小平方法求解。解的矩阵表示式为 d G G G m T T 1][?-= （1-2-5） 上式具体计算时可用奇异值分解方法 G=U ∧V T 最后，得 m ?=（G T G ）-1G T d=V ∧-1U T d （1-2-6）
地震属性体处理 1、分频处理属性 分频处理属性可将地震振幅和属性数据转换成更为清晰的地下地质图像，识别薄层或能量衰减区。将各地震道分解成不同的频带成分，有助于突出复杂的断裂体系以及储层的分布特征。分频处理的技术主要是通过 “Gabor-Morlet” 子波对复数地震道进行谱 分解，类似于小波变换。用来帮助地质家和解 释人员进行如下的勘探研究工作： （1）薄层检测以及薄层厚度估计； （2）衰减分析——直接进行油气检测 （3）提高地震分辨率 该方法通过连续的时频分析来描述时间-- 频率的瞬时信号能量密度。与以往常规的谱分 解使用离散傅立叶变换不同，该方法使用 Gabor-Morley 子波来提高时间-频率的分辨率。 提供了两种计算瞬时能量的方法：等空间中心频率和倍频程频率。输出结果可以分解成多种属性体：时间-频率体、时间切片，然后进行分析。 2、地震属性分析 地震属性分析使我们获得更多极有价值的多方位信息，从而使油藏的描述更准确、更细致。帕拉代姆地震属性库包括丰富的地震属性，如振福包络、瞬时频率、吸收系数以及相对波阻抗等20多 种复地震道（Hilbert ）属性、多道几何属性，谱分解属 性和用户自定义属性见图。这些地震属性可分别表征地 震影像的不同特征，从而使解释人员以少量的工作即可 获得大量的地质信息，其中多地震道几何属性包括倾角 体、方位角体、非连续性和照明体。这些属性旨在强化 地震影像的非连续性特征，因此对识别地质体的构造特 征（如断层）、地层边界、河道和地质体的几何样式十 分有效。在这些属性体提取的基础上，利用PCA 主组分 分析技术进行属性优化分析，同时也可借助多属性体交 会VXPLOT 识别异常体。通过多属性体交汇、神经网络 测井参数反演、多属性体的波形分类以及变时窗/等时 窗的地震相划分等综合技术，并借助多属性体立体可视 化浏览技术实现对地下构造、地层和储层岩性的综合解 释。 光照体属性 常用提取的地震属性有信号包络、瞬时频率、瞬时相位、相对波阻抗、分频处理等。
大地震来临之前预测方法 副标题： 作者：本站整理文章来源：网络点击数：715 更新时间：2003-12-10 切事物都有过去，一切事物都有未来。未来是一个令人神往的字眼，它给人以希望.未来又是个充满未知的字眼，它给人以迷茫。 要创造光辉的未来，就要正确地预测未来。凡事预则立，不预则废。未来可知吗？科学表明，在一定意义上，未来可知，未来可测。未来是怎么预测的呢？ 日月有常，星辰有行，日月星辰的运行是有一定的周期，它们周而复始，循环往复。月亮围着地球转，地球围着太阳转，年年月月，月月年年，运行不息。月亮的盈亏引起了潮汐的变化。大海永不止息地潮涨潮落，大潮小潮，早潮晚潮，都有一定的规律。 根据月球运动周期，我们可以预测月相变化，预测海潮涨潮退潮的时间。 根据太阳地球和月亮的运动周期和相对位置，天文学家能够预测日食和月食。科学家们早就预测到，97年3月9日，将发生本世纪最后一次日全食。就在这一天9时09分左右，一个巨大的黑影准时地遮没了太阳。月晕而风，础润而雨。 任何事物的发生，常常有一个蕴酿的过程。蕴酿过程中，总会有某些前兆发生，认识了前兆现象我们就可以预测未来。 地震是地球内部运动的一种自然的现象。它造成极大的破坏力。地震来临时，地光闪烁，天崩地裂。高楼夷为平地，城市变为废墟。地震是地质运动的激烈变化，它的形成，伴随着能量的积累，这个过程会引起大地的形变，会引起声、光、电磁等物理、气象和生物的变化。认识了这些变化，就能预测地震。 １９７３年，中国地震工作者发现了辽宁营口地区的地质形变、地磁等等的变化，根据这些地震前兆，加强了观测的频度和密度。１９７４年１２月间，人们发现了井水变浑、冒泡等水文反常现象，人们又发现蝴蝶、蚂蚁等昆虫活动反常，家禽和家畜也活动异常。２月初，地震工作者又发现了水中氡的含量剧变、地电反应异常，有的井水自喷，有的井水则猛然下降。对这些现象进行了综合分析后，地震工作者认为地震在即，２月４日发出了紧急预报，不久地震发生了。由于预报及时，大大减少了人员的死伤。现在，地震专家应用了灵敏的仪器，不断地收集着地质变化的各种信息，帮助预测和预报地震。 人们常说，天有不测风云，其实，地震可测，风云也可测。风云变幻，但是，它也会有各种征兆可寻。云，是大气流动和水汽多少的反映。因此，观云可以测天。天高云淡，说明水汽稀少，秋高气爽，干燥无雨。如勾的卷云，表明一场风雨就要来临，正如民谚所说的，天上勾勾云，地下雨淋淋。天上布满鱼鳞斑的高积云，表明气流稳定，预示着一个大好的晴天。农民会据此预测，“天上布满鱼鳞斑，明日晒谷不用翻。”事物之间常常有着内在的联系。 １９８２年２月底，墨西哥的爱尔.基琼火山爆发.科学家们预测着火山喷发会引起什么
地震反演的类型 1．1 反演的分类 1）从所利用的地震资料来分可分两类：叠前反演和叠后反演; 2）从测井资料在其中所起作用大小可分为四类：地震直接反演，测井控制下的地震反演，测井—地震联合反演和地震控制下的测井内插外推； 3）从实现方法上可分三类：直接反演、基于模型反演和地震属性反演。 4）从反演模型参数来分主要有：储层特性(如:孔隙度、渗透率、饱和度等)反演、岩石物性反演、地质结构反演、各向异性参数反演、阻抗反演以及速度反演等; 5）从使用的数学方法可分为：最优化拟合反演、遗传算法反演、蒙特卡罗反演、Born近似反演、统计随机反演以及基于神经网络的反演等。 1．2几种主要反演方法的概述 叠前反演尚处于研究试验阶段，而叠后地震反演近年来快速发展，形成了多种技术。下面简要介绍几种主要反演方法：直接反演（递推反演和道积分反演）、基于模型反演、地震属性反演、测井约束反演和叠前AVO反演。 1．2．1直接反演 两种基本做法：递推反演和道积分反演。 1）递推反演：递推反演是一种基于反射系数递推计算地层波阻抗的直接地震反演方法。它完全依赖于地震资料本身的品质，地震资料噪音对反演结果敏感，影响大，地震带宽窄会导致分辨率相对较低，难以满足储层描述的要求。典型的有Seislog,Glog，稀疏脉冲反演(实现方法又有MED,AR,MLD,BED方法等)等;Seislog,CLOG等使用测井信息后，只获得剖面上关键点的低频分量，整个剖面上的低频信息要靠内插来求得。 优点：计算简单，递推列累计误差小。其结果直接反映岩层的速度变化，可以以岩层为单元进行地质解释。缺点：由于受地震固有频率的限制，分辨率低，无法适应薄层解释的需要；其次，无法求得地层的绝对波阻抗和绝对速度，不能用于定量计算储层参数。这种方法在处理过程中不能用地质或测井资料对其进行约束控制，因而其结果比较粗略。 2）道积分反演：是以反褶积为基础的地震直接反演法。道积分是利用叠后地震资料计算相对波阻抗的直接反演方法，它无需测井资料控制，计算简单，其结果直接反映了岩层的速度变化，但受地震资料固有频宽的限制，分辨率低，无法适应薄层解释的需要，无法求得地层的绝对波阻抗和绝对速度，不能用于定量计算储层参数。 优点：能比较完整地保留地震反射的基本特征（断层、产状），不存在基于模型方法的多解性问题，能够明显地反映岩相、岩性的空间变化，在岩性相对稳定的条件下，能较好地反映储层的物性变化。 缺点：由于受地震频带宽度的限制，递推反演资料的分辨率相对较低，不能满足薄储层的研究需要。 1．2．2基于模型的反演 1）基于模型的反演：就是从地质模型出发，采用模型优选迭代扰动算法(广义线性或非线性最优化算法)，通过不断修改更新模型，使模型正演合成地震资料与实际地震数据最佳吻合，最终的模型数据便是反演结果。 实现方法有广义线性反演(GLI)(Cooke,1983);宽带约束反演(BCI)(Martinez,1988);地震岩性模拟(SLIM)(Ge lfand,1984);具有全局优化特点的遗传算法、模拟退火法(Smith等1992:Sen和Stoffa,1995);蒙特卡罗搜索法(Cary和Chapman,19 98)以及人工神经网络法(Ca lderron-Macias 等，1998)等。 目前,以模型为基础的反演方法一般都是依据测井及地质资料建立初始模型,通过广义线性反
地下流体异常地震预测法 ——李宇杰土木1202 20121167 引言 多年的地震监测预报实践表明(在强震前可观测到十分明显的水位、水氡等测项的中期和中短期异常变化，但实践表明只采用看图识字方法来识别测项异常是不够的，往往会造成异常判定的失误。在前人研究的基础上按照同一化技术思路、差分从属函数、变差和变化率等方法来判定地下流体的中期及中短期前兆异常。由不同方法所给出的异常变化彼此之间具有较大的可比较性，从而有利于地震的综合分析预报。 一理论提出 最早提出系统发展地下流体观测的国家是苏联，但直到1966年4月塔什干地震发生后，苏联才开始大规模开展以水文地球化学为基础地震预测方法的研究。 美国的地下流体研究方向主要是地下水水位、水温、孔隙压力大小和水氡等。正在实施的“地球透镜计划”可用于研究大地构造和地壳中流体间的关系。 日本对地下流体的观测包括水位、水温及地下水的化学性质，重点发展深度(>1000 m)探测及对深度流体的直接取样以避开地表因素的干扰。 内陆最早由刘耀炜（现任中国地震局地壳应力研究所监测预报研究室主任、研究员）等实践并发展。对1966年3月26日宁晋百尺口6.2
级地震提出了较好的预测，这是中国第一次对强余震成功预测，而主要依据就是震中区水位急剧涨落、井水翻花、发浑等临震异常. 二理论 作为地下流体前兆，观测项目主要有水位、水温(地热)、水氡与水汞等，近年来开始并得到推广的还有气氡、气汞和H2、He、CO2等.根据布雷斯等人得到的有效应力定律关系可知，孔隙压力的增加会使有效正应力减小，而对剪应力分量没有影响.孔隙压力增大，断层面上的抗剪强度下降，当断层面上积累的剪切力变大时，断层发生滑动；另外，流体压力所产生的应力转移会使剪应力增加，进一步加强岩石不稳定性（王博等，2008）。Matthai等认为流体压力近似于静岩压力时，断层膨胀或蠕动重复发生，会使各封闭的孔隙连通，形成网络.进而产生流体，导致较弱的断层岩石局部发生水压破裂。另有学者指出，水的扩散会在岩体结构面内产生相应的孔隙压力，从而造成弱化断层的粘滑不稳定，导致岩体破裂并诱发地震（施行觉等，1986；颜玉定等，2005)。陆明勇认为孔隙压力及动水压力的变化异常可以作为中短期的地震预测(陆明勇等，2005). 三异常判定方法简介
地震反演方法概述 地震反演：由地震信息得到地质信息的过程。 地震反射波法勘探的基础在于：地下不同地层存在波阻抗差异，当地震波传播有波阻抗差异的地层分界面时，会发生反射从而形成地震反射波。地震反射波等于反射系数与地震子波的褶积，而某界面的法向入射发射系数就等于该界面上下介质的波阻抗差与波阻抗和之比。也就是说，如果已知地下地层的波阻抗分布，我们可以得到地震反射波的分布，即地震反射剖面。即由地层波阻抗剖面得到地震反射波剖面的过程称为地震波阻抗正演，反之，由地震反射剖面得到地层波阻抗剖面的过程称为地震波阻抗反演。 叠前反演主要是指AVO反演，通过AVO反演，可以获得全部的岩石参数，如：岩石密度、纵横波速度、纵横波阻抗、泊松比等。叠前反演与叠后反演的根本区别在于叠前反演使用了未经叠加的地震资料。多道叠加虽然能够改善资料的品质，提高信噪比，但是另一方面，叠加技术是以东校正后的地震反射振幅、波形等特征不随炮检距变化的假设为基础的。实际上，来自同一反射点的地震反射振幅在不同炮检距上是不同的，并且反射波形也随炮检距的变化而发生变化。这种地震反射振幅、波形特征随炮检距的变化关系很复杂，主要原因就在于不同炮检距的地震波经过的地层结构、弹性性质、岩性组合等许多方面都是不同的。叠加破坏了真实的振幅关系，同时损失了横波信息。叠前反演通过叠前地震信息随炮检距的变化特征，来揭示岩性和油气的关系。叠前反演的理论基础是地震波的反射和透射理论。理论上讲，利用反射振幅随入射角的变化规律可以实现全部岩性参数的反演，提取纵波速度、横波速度、纵横波速度比、岩石密度、泊松比、体积模量、剪切模量等参数。 叠后地震剖面相当于零炮检距的自激自收记录。与叠前反演不同，叠后反演只能得到纵波阻抗。虽然叠后反演与叠前反演想必有很多不足之处，但由于其技术方法成熟完备，到目前为止，叠后反演仍然是主流的反演类型，是储层预测的核心技术。 介绍几种叠后反演方法： 1）道积分：利用叠后地震资料计算地层相对波阻抗（速度）的直接反演方法。因为它是在地层波阻抗随深度连续可微的条件下推导出来的，因而又称为连续反演。 原理简述： 上述公式表示，反射系数的积分正比于波阻抗Z的自然对数，这是一种简单的相对波阻抗概念。 适用条件及优缺点 与绝对波阻抗反演相比，道积分的优点：1.递推时累积误差较小；2.计算简单，不需要反射系数标定；3.无需钻井控制，在勘探储气即可推广使用。 缺点：1.由于这种方法受到地震固有频宽的限制，分辨率低，无法适用于薄层解释的需要；2.需要地震记录经过子波零相位化处理；3.无法求得地层的绝对波阻抗和绝对速度，不能用于定量计算储层参数；4.这种方法在处理过程中不能用地质或测井资料对其进行约束控制，因而结果比较粗略。 2）递推反演方法：根据反射系数进行递推计算地层波阻抗或层速度，其关键在于由原始地震记录估算反射系数和波阻抗，测井资料不直接参入反演，只起到标定和质量控制的作用。因此又称为直接反演。 原理简述： 利用以上公式，可以从声波时差曲线及密度曲线上（没有密度曲线时可以利用Gardnar 公式进行换算）选择标准层波阻抗作为基准波阻抗，将反褶积得到的反射系数转为波阻抗。