三维土_结构动力相互作用体系分析的两步时域显式波动有限元过程.docVIP

doi: 10． 3969 / j． issn． 0253 － 4967． 2012． 03 doi: 10． 3969 / j． issn． 0253 － 4967． 2012． 03． 009 三维弯曲界面炮集 波反射时距特征分析 P 郝重涛* 宋利虎 姚 陈 蔡明刚 ( 中国地震局地质研究所，地震动力学国家重点实验室，北京 100029) 对于三维反射地震，文中将三维弯曲界面视为倾向和倾角变化的三维倾斜界面的包络，并 摘 要 依据三维倾斜界面反射提出计算三维弯曲界面 P 波反射时距的一种新方法。该方法基于倾角 CDP 理论，从界面反射点出发，给出地表 CDP 网格点，由检波点位置计算炮点位置，得到三维地震的炮点 网格和检波点网格的分布关系，最终计算 P 波的三维反射时距。炮集数值模拟的结果表明，该方法有 较高的计算效率，便于分析反射时距随炮点位置的变化，特别能揭示界面反射点位置和时距的关系。 三维弯曲界面和三维倾斜界面的反射时距特征对比表明，仅从同一测线方位的时距难以将这两种类 型的界面反射区分开来，但比较不同方位的时距差异可以解决这一问题。 关键词 三维地震勘探 三维弯曲界面 P 波反射时距 界面反射点 中图分类号: 文献标识码: A 文章编号: 0253 － 4967( 2012) 03 － 0477 － 11 P35． 63 0 引言 随着三维地震勘探技术的深入发展，地震勘探面临着越来越复杂的勘探目标和越来越高的 勘探精度的挑战( 张中杰，2002) 。目标体的复杂性突出表现在地下发育复杂构造，如倾斜界面 和背斜、向斜等，地层倾角变化大。为了对复杂构造的关键部位进行精确成像，这就要求我们提 高对地震波在复杂构造中传播规律的认识，进而采用一些特殊的、有针对性的勘探方法，比如突 破传统的建立在水平层状介质理论基础上的地震勘探技术的限制 ( 佘德平，2004 ) 。为了达到 这一目标，需应用地震数值模拟来了解地震波在三维复杂构造中的传播规律。 数值模拟技术是研究各种地震地质条件下构造形态、岩性等各种地质因素与地震波响应特 征之间关系的一门技术，是在假定地下介质构造模型和相应的物理参数已知的情况下，模拟研 究地震波在地下介质中的传播规律，并计算在各观测点所应观测到的数值地震记录的一种地震 模拟方法( 裴正林等，2004; 佘德平，2004 ) 。反射时距模拟是地震波场数值模拟的重要组成部 分，是研究复杂地区地震资料采集、处理和解释的有效辅助手段。在地震资料采集方面，数值模 拟可用于指导观测系统设计，并对野外观测系统进行评估和优化; 在地震数据处理方面，数值 模拟可以检验数据处理方法的正确性，为反演提供理论基础 ( 裴 正 林 等，2004; 佘 德 平，2004; 刘洋等，2008) 。 〔收稿日期〕 2012 － 04 － 09 收稿，2012 － 05 － 11 改回。 〔基金项目〕 地震动力学国家重点实验室开放基金( LED2011B02) 和国家自然科学基金( 41104076、 共同资助。 * 通讯作者: 郝重涛，副研究员，E － mail: haoct@ ies． ac． cn。 三维地震技术兴起于 20 世纪 70 年代，基于地下构造本身是一个三维实体，用三维的方法进行研究，才能得到对地质结构的全面正确认识( French 三维地震技术兴起于 20 世纪 70 年代，基于地下构造本身是一个三维实体，用三维的方法 进行研究，才能得到对地质结构的全面正确认识( French，1974) 。三维地震勘探中的 “三维”包 含两方面内容: 一是三维观测系统，同时布设多条测线，即在地表采用网格面接收，炮点和检波 点按一定规律分布在平面上; 二是地质体的空间位置和射线传播路径的描述是三维的。近年 来，三维地震数值模拟方法发展迅速。运动学方面，传统的有打靶法和弯曲法，在此基础上逐渐 发展起来的有基于程函方程、费马原理和惠更斯原理的波前扩展类方法( Vidale，1988; 徐昇等， 1996; 张霖 斌 等，1996 ) ，基于图论和费马原理的最短路径法 ( Nakanishi et al．，1986; Moser， 1991; 王辉等，2000; 赵爱华等，2003，2004; 张美根等，2006) 。上述方法用网格剖分来模拟地 下构造，当网格稀疏时，计算结果会出现较大偏差; 而为了提高计算精度做加密剖分网格等修 正，会显著增加计算量，降低效率。因此受计算机内存和运算速度的影响，这些方法很难应用到 实际工作中; 且三维地震中涉及到多个炮点，炮点位置改变，相应的追踪到的反射点分布也发 生变化，不利于后续的叠加处理。 针对界面反