折射地震资料处理软件研制及应用.pdf

折射地震资料处理软件研制及应用 梁志强 刘国华 (长安大学地质工程与测绘工程学院 西安710054) 摘要一直以来，折射地震勘探数据的解释都是通过人工对比．人工读取数据解释戒图。这种 方法效率低．解释周期长，数据凄取精度不高，为此．我们研发了一套基于t。法原理的折射地震 数据自动解释和成像系统．本系统基于VC和fortran两种语言设计而成，具有界面良好，功能齐 全。使用方便，快捷．运算速度快等优点．通过有限元数值模拟的模型资料和实际野外采集的地 震资料验证与应用．都得到了良好的结果。 Fortran 关键诃to法折射地震Vc 引言 折射地震勘探技术主要应用于地质填图、探测基底深度、确定第四系及第三系厚度，以 及反射地震勘探中的低速带调查等，有时也用于地学大断面中的深部Moho面及上伏构造的 勘探研究。以往的折射地震处理解释软件大都采用人工对比折射渡，读取数据解释成图，这 种方法不仅计算速度慢，而用解释精度低。我们用VC和Fortran两种语言进行混合编程研制 了一套基于相遇时距曲线法原理的折射地震数据自动解释和成像系统。由于以下五太模块组 成：①折射地震剖面成图；②数据的预处理及编辑；③图形显示；④反演计算；⑤偏移及成 果输出。通过本系统对有限元数值模拟的地震记录和实际勘探地震数据资料的处理及解释，说 明本系统具有界面友好，功能齐全，使用方便、快捷，运算速度较快等优点。 该系统也可用于工程勘探中铁路公路选线、灾害地质体调查、隧道勘察、岩石硬度估计 和供水水源勘察等方面。 1基本原理 浅层地震勘探折射波法中常用的解释方法有t。法、准旅行时法、时间场法、哈莱斯法、共 轭点法、时间项法、波路计算法等，本系统采用t。法反演地质剖面，其基本原理如图1所示： 图1中的了1。．和r，。为直达波时距曲线，丁“和丁：n是界面R的折射波时距曲线。 设折射界面的曲率半径比其埋深大得多，△SⅣQ可近似看作等腰三角形，高h(法线深 度)可用下式表示： h=Vlto／(2cosi)=Kto (1) 式中，fo一72^+T2日一T=T。^一(了1一丁2B)=丁¨一△f (2) ●______________一 (3) K=V1／2cosi=VlVz／(2√明一V}) 利用丁¨和丁：口，在测点T量出丁一丁2s(工)=zlt(z)，然后从T^：b)减去z3tb)，连接 一52— 各点的t。Cr)即成to(z)关系曲线，如图1所示。 y，可由直达波时距曲线丁。。和丁。。的斜率求取，为 了求取K值，还必须由y。求出n，为此可类似求t。 b)的方法，作出差数时距曲线，令 口(工)一T2^+丁z口一丁一T2^+△￡(4) 利用了’2B量出丁一T2B(z)=出(z)，将At(T)和71：^ (z)相加得到目(z)，连接这些点即为差数时距曲线日 o)。当界面倾角不大时．V2_22tr／At．由0(z)曲线 可求出V。。 求出界面速度yz后，可根据公式(3)求出K值，然 后由t。(z)曲线和公式(1)求出各点的法线深度ho)。 以各接收点为圆心，各点的法线深度h(z)为半径画弧， 各弧的包络线即为所要求的折射界面。 图l二层构造殛其时距曲线 2实现方法 “法解释示意围 在地球物理数据处理中，多使用Fortran语言编写 程序．Fortran语言是一种面向过程的语言，是目前流行较广的适用于科学计算的高级语言， 但它在实现人机对话、交互处理方面很不理想，界面也不够美观，图形处理功能也比较欠缺。 VC是一种面向对象的高级编程语言，它不但界面美观，图形处理功能强大。而且在实现人机 对话、交互处理方面都很方便，与硬件底层打交道更显示出它的强大优