基于射线理论的各向异性煤层地震响应.pdfVIP

计 算 物 理 Vo1．31．NO．4 第 31卷 第 4期 2014年 7月 CHINESE JOURNALOFCOMPUTATIONALPHYSICS Ju1．，2014 文章编号：1001—246X(2014)O4_0444_05 基于射线理论的各 向异性煤层地震响应 李 勤 ， 李庆春 ， 张 林 (1．西安科技大学 地质与环境学院，陕西 西安 710054；2．长安大学 地质工程与测绘学院，陕西 西安 710054 3．陕西省地质调查中心，陕西 西安 710016) 摘 要：以VTI型煤层为研究对象，从介质中的地震波的射线方向、走时计算等多个方面展开讨论，分析VTI型煤 层的试射法射线追踪技术，充分利用VC++的面向对象类框架泛型编程的特点和Fo~ran数值计算方面的优势 ，通 过混合编程实现VTI型煤层中的地震波试射法射线追踪 ；通过各向同性煤层和VTI型煤层的射线追踪试算对 比， 论证对VTI型煤层进行各向异性射线追踪的必要性，同时验证本文VTI型煤层射线追踪方法的正确性．最后，通过 计算时间成本比较，表明该方法计算速度较快，能满足地震正演模拟的要求．该方法适用于各向异性煤层的多波射 线追踪． 关键词 ：VTI型煤层 ；射线追踪 ；各向异性 中图分类号 ：P631 文献标志码 ：A 0 引言 煤层裂隙是 国内外煤矿地震勘探研究的重点和难点，裂隙发育复杂 ¨ ，其中一种较为典型的裂隙系统 可近似为平行于地表 的裂隙系统 ，这种裂隙系统称为 VTI型，而这种煤层就称之为 VTI型煤层 ．VTI型煤 层作为各 向异性煤层的一种常见模型，裂隙发育导致的各 向异性会对地震波的记录产生较大影响 ．但 是，由于煤层相对地震子波来说较薄，各 向异性煤层地震波正演模拟也相对复杂．目前，对于煤层进行各向异 性的研究主要集中在实验室分析。。卜 ，波动方程法地震波模拟 及 AVO特征分析 ¨ ”．射线追踪法和波 动方程法各有优缺点．波动方程法相对来说，所生成的地震记录波形丰富，信息量大，比较接近实际地震记 录；而射线追踪法适合于追踪单一波形，所得模拟地震记录可直接用于地震波各向异性响应特征研究．射线 追踪法地震波模拟由于其具有显示直观、高效计算的特点，在石油地震勘探方面取得 了显著成效 ¨ ，然而 却一直未见用于各向异性煤层中．试射法是一种基于射线的算法 ，也是一种经典的射线追踪算法 ，有很大的 发展空间 ．目前的地震记录正演模拟软件很多，如Tesseral等 ，它们大多是对层厚较大的模型进行模拟 ，还 没有对 VTI型煤层的正演模拟模块． 本文以各向异性煤层为研究对象，采用试射法对煤层中地震波进行射线追踪，分析VTI型煤层 中走时及 射线路径的传播特征 ；并将VC++与 Fortran语言有机结合 ，实现混合编程 ，研发基于VTI型煤层的地震记录 正演模拟模块 ，将其作为一个开放式模块 ，可加载到其它软件系统 中，提高正演模拟的精度和效率． 1 各 向异性煤层 中试射法射线追踪 相对于一般 VTI型介质来说 ，VTI型煤层 的结构更为复杂，射线方 向的计算精度要求更高．本文将 VTI 型煤层的试射法射线追踪分为以下几个步骤进行 ：给定群角进行射线试射 ，进行群相关系转换 ，寻找射线 出 射角度 ，计算走时． 1．1 群相关系 各 向同性介质中，相速度与群速度传播方 向相同，大小相等，而在各向异性介质 中，波的相速度指波前传 收稿 日期 ：2013—07—27；修回 日期 ：2013—11—25 基金项 目：国家 自然科学基金41374145)；高等学校博士学科点专项科研基金 (20120205130002)；陕西省教育厅专项科研计 划 (12JK0470)资助项 目 作者简介 ：李勤 (1979一)，女，博士，讲师，从事应用地球物理信息技术工作，E-mail：eriliqin@126．tom 第 4期 李 勤等 ：基于射线理论的各向异