《工程与环境物探》课程设计.doc

工程与环境物探课程设计学生姓名： 学 号： 学 院： 专业班级： 指导教师： 二零一零年一月十九日目 录第一章 折射波法基本原理和方法31.1折射波法基本理论31.2 折射波法观测系统41.3 折射波分层解释的法4第二章 实验仪器设备及实验方法72.1 仪器设备72.2 采用的观测系统9第三章 现场资料采集实验步骤10第四章 软件说明和数据处理过程104.1 软件功能104.2 数据处理过程11第四章 数据处理结果及图形绘制23参考书目24第一章 折射波法基本原理和方法1.1 折射波法基本理论以水平界面的两层介质进行简要的说明，假设地下深度为h，有一个水平的速度分界面r，上、下两层的速度分别为v1和v2，且v2v1。如图11所示。从激发点o至地面某一接收点d的距离为x，折射波旅行的路程为ok、ke、ed之和，则它的旅行时t为： （11）图1-1 水平两层介质折射波时距曲线为了简便起见，先作如下证明：从o，d两点分别作界面r的垂线，则oadgh，再自a、g分别作ok，ed的垂线，几何上不难证明bakegfi，因已知，所以： （12）即 和 （13）上式说明，波以速度v1旅行bk（或ef）路程与以速度v2旅行ak（或ec）路程所需的时间是相等的。将式（13）的关系和式（11）作等效置换，并经变换后可得： （14）这就是水平两层介质的折射波时距曲线方程。它表示时距曲线是一条直线，若令x0，则可得时距曲线的截距时间t0（时距曲线延长与t轴相交处的时间值） （15）式（15）表示出界面深度h和截距时间t0之间的关系，当已知v1和v2时，可以求出界面的深度h。1.2 折射波法观测系统根据折射波场形成条件和特征，折射波观测系统必须避开盲区，且要把接收部分尽量放到待测地层折射波区范围。当水平层状介质满足折射条件的前提下，固定一个激发点，将排列沿测线由近及远进行时距观测，将得到由浅入深各层介质的地震波信息。在地震记录上可观察到各层介质折射波的动力学特点， 从而判别层间干涉、波形置换特征。在时距曲线上将反映出各层介质折射波运动学的空间分布规律及介质的物理力学性质。为了消除表层不均匀及界面起伏的影响，往往采用相遇时距曲线观测系统。相遇时距曲线观测系统如图2-2所示，同一观测地段分别在两端o1 和o2 点激发， 图2-2此观测系统采集的地震记录，可得两支方向相反的时距相遇曲线s1和s2。ae段折射， o1 o2接收， ea段折射，o1o2接收。其优点可弥补单一时距曲线的不足，可以从不同方向反映界面变化。1.3 折射波分层解释的法折射波解释法是常用的地震折射波解释方法，它是针相遇时距曲线观测系统采集发展起来的解释方法。t0法解释的主要原理与方法如下：t0法又称为t0差数时距曲线法，是解释折射波相遇时距曲线最常用的方法之一。当折射界面的曲率半径比其埋深大得很多的情况下，t0法通常能取得很好的效果，且具有简便快速的优点。如图13所示，设有折射波相遇的时距曲线s1和s2，两者的激发点分别是o1和o2， (a) (b)图1-3 t0法折射界面示意图若在剖面上任意取一点d，则在两条时距曲线中可以分别得到其对应的走时t1和t2，从图中可以得到： （16）且在o1和o2点，时距曲线s1和s2的走时是相等的，称之为互换时，用t表示，则有： （17）当界面的曲率半径远大于其埋深时，图中的bdc可以近似地看作为等腰三角形，若自d点作bc的垂直平分线dm（dm即为该点的界面深度h），于是有： 和 （18）将公式（16）中的t1和t2相加，并且减去（17）式，再将（18）式代入后可以得到： （19）式（19）便是任意点d的t0值公式，由此可得d点的折射界面法线深度h为： （110）令和则式（110）可以写为： （111）因此只要从相遇时距曲线中分别求出各个观测点的t0值和k值，就可以得出各个点的界面深度h。从上述的公式可以看出，只要从时距曲线上读取t1，t2和互换时t，就可以算出各个点的t0值，并可以在图上绘制相应的t0(x)曲线（13(b)中所示）。关于k值的求取：根据斯奈尔定律可将k值表达式写成： （112）由公式（112）可以看出，只要求得波速v1和v2则很容易得出k值。其中v1通常可以根据表层的直达波速度来确定，因此关键就是v2值的求取，为此引出参数时距曲线方程：令 （113）对式（113）两边对x求导，可得： （114）式中和分别为上倾方向时距曲线s1和下倾方向时距曲线s2的斜率（即视速度v*的倒数）。根据公式： 和 因为是同样o1o2内观测段，设上倾方向x为正，下倾方向x为负，则： 和 它们分别对x求导有如下形式： 和 （115）将（115）代入（114）式中，经一些变换后可得： （116）于是可以求得波速v2为： （117）当折射界面倾角小于15时，可以近似的写成 （118）因此只要根据（113）式在相遇时距曲线图上构置(x)曲线，并求取其斜率的倒数，则可以根据（113）式得出波速v2，进而代入（112）式中求得k值。知道了k值和各个观测点的t0值后，可以根据（111）式计算出各个点的界面深度h。然后，以各观测点为圆心，以其对应的h为半径画弧，可以得出一系列的圆弧，作这些圆弧的包络线即为折射界面的位置。或由、和计算h和i，再用公式计算得到真深度h0。第二章 实验仪器设备及实验方法2.1 仪器设备仪器设备主要包括: 振源、检波器、工程检测仪、触发开关、电源、大线等,具体布置见图2-1。sws-1a型多功能检测仪就是陆地进行折射波勘察的主要仪器设备之一，见图2-2。 图2-1 仪器设备布置折射波法振源类型：炸药(或雷管)，重锤和垫板。 检波器：竖向检波器，其自振频率约为38hz，所测地震波主频率的20-25。一般折射波法的自振频率较高。触发开关：分为外触发和内触发两种。折射波勘察使用外触发，用安置在震源附近触发传感器或检波器进行触发。图2-3 sws多功能检测仪电源：可用12伏的直流电或12伏60安时的蓄电池。工程检测仪：不同型号的检测仪其功能不尽相同，但目前来看新型的检测仪都具有如下的基本性能和功能： 性能好坏主要取决于a/d变换器和储存器。ad 转换器性能要求是转换速率快，分辨率高（数字化输出的位数高）和可靠性高。存储器的主要要求是存储速度快、存储容量大。体积小和耗能小。 功能主要实现：（1）多通道；（2）ad转换；（3）记录和显示；（4）信号分析功能（软件系统）； （5）打印和数据输出系统。 2.2采用的观测系统根据折射波场形成条件和特征，折射波观测系统必须避开盲区，且要把接收部分尽量放到待测地层折射波区范围。当水平层状介质满足折射条件的前提下，固定一个激发点，将排列沿测线由近及远进行时距观测，将得到由浅入深各层介质的地震波信息。在地震记录上可观察到各层介质折射波的动力学特点， 从而判别层间干涉、波形置换特征。在时距曲线上将反映出各层介质折射波运动学的空间分布规律及介质的物理力学性质。为了消除表层不均匀及界面起伏的影响，实验中我们采用相遇时距曲线观测系统。相遇时距曲线观测系统如图2-2所示，同一观测地段分别在两端o1 和o2 点激发，图2-2此观测系统采集的地震记录，可得两支方向相反的时距相遇曲线s1和s2。ae段折射， o1 o2接收， ea段折射，o1o2接收。其优点可弥补单一时距曲线的不足，可以从不同方向反映界面变化。本实验观测系统设计如下：检波道数：24道采样间隔：1ms采样总数：1024偏移量：5m道间距：1m采样长度：1024ms第三章 现场资料采集实验步骤：布置地震测线。测量确定测线起点和终点位置，布设测线标志，确定激发点和接收点标志。：按测线起点方向，将检测仪器放在适当的位置，并寻找激发点和相应的接收点。在激发点选择合适的震源类型，并在激发点附近布置触发的传感器或检波器，将其与仪器连接，用于仪器记录时触发；在接收点排列方向上布置大线和检波器，并将大线与仪器连接；实验我们主要采用锤击震源，将振源和一组检波器布置在一条直线上（纵测线），排列相对较长。：仪器系统通电后启动采集软件，按观测系统布设的参数设置各采集参数。检波道数主要由有效传播范围和仪器通道数决定，实验中仪器采取24道接收，每道采样数1024，采样时间间隔由采样定理确定，实验中我们采用时间间隔为0.500ms，道间距为2m，偏移距为6m。采样排列长度为（24-1）*2=46m。：检查激发、接收和仪器是否正常，还要注意周边是否有人为震动源，各部分正常，且无强烈的震动源后，提醒周边人员注意，保持相对静止，如果是炸药震源让他们尽快撤离危险地带或采取防护措施。保证安全无误后，仪器操作员发出指令震源激振，仪器接收地震记录。并检查地震记录的质量。接收记录满足要求，可以保存；如不正常，及时查明原因，重新采集。 第四章 软件说明和数据处理过程4.1 软件功能数据主要采用的软件有：(1)记录编辑（rcdedtv.exe）(2)记录滤波和增益处理（rcdprcv.exe）(3)噪音及干扰波切除（mutnoiv.exe）(4)地震波拾取（pktm.exe）软件。记录编辑（rcdedtv.exe）可以对记录进行修饰和拼结，截短等，记录滤波和增益处理（rcdprcv.exe）可以对记录进行频谱分析，用不同的滤波类型进行滤波处理，噪音及干扰波切除（mutnoiv.exe）是对记录上存在的触发噪音，声波，面波等干扰波进行切除，突出有效波，以减小干扰波的影响。4.2 数据处理过程（一） 实验采集到的原始数据1、第一次激发波形图d20.dat2、第二次激发波形图d13.dat3、第三次激发波形图d21.dat4、第四次激发波形图d33.dat（二） 记录编辑在系统主界面下选择记录编辑（rcdedtv.exe）回车，即进入“记录编辑”主菜单load(调入记录)按地震记录文件的驱动器号、路径和文件名输入，回车确认。选择modify record(记录修饰处理)。1、d02切除后图形2、d13切除后图形3、d21切除后图形4、d33切除后图形（三） 记录滤波和增益处理进入记录滤波和增益处理（rcdprcv.exe）功能后，主菜单显示见图选择spectrum做频谱分析,进入该功能后，在显示当前记录的屏幕右上角出现子菜单，见图。以 选择子菜单中的各功能项，回车确认。进入spectrum(谱分析)后在视频右上角出现选用时间窗口的菜单，见图。同时在视频左上角第一道顶端有一时段标。以 移动光标，、和pgup、pgdn改变标内数值，与此同时左上角时段标也随之移动，帮助观察所选的道及时窗。确认后按回车键，屏幕呈现频谱图形。在彩显上，绿色为全道谱成分；浅蓝色为时窗内的谱成分。在该时窗下，进行谱分析。各数据各道频谱分析见下图：1、d02第一道频谱分析2、d02第六道频谱分析3、d02第二十四道频谱分析4、d13第一道频谱分析5、d13第六道频谱分析6、d13第二十四道频谱分析7、d21第一道频谱分析8、d21第六道频谱分析9、d21第二十四道频谱分析10、d33第一道频谱分析11、d33第六道频谱分析12、d33第二十四道频谱分析进入filtering(滤波)后在视频右上角出现提示菜单，见图，同时光标在屏幕左上角第一道始端。按s键，向右移动光标至具有代表性的道上，回车，屏幕展示该道的频谱图，并在屏幕右上角显示滤波开关及可以设置滤波类型及响应参数的窗口。按、箭头键可以改变滤波的关断和滤波类型。滤波器选择有4档，我们选用的是band pass filter带通滤波进行滤波。利用、箭头移动光标选择条目，利用、和pgup、pgdn改变数值，以此确定滤波参数，在此过程中按f2键屏幕上的频谱内根据所确定参数适时显示该道滤波后的频谱和被压制的谱成分，以此方便调整滤波参数的准确性。确认滤波参数后，回车即实现滤波操作，屏幕上显示出滤波后的记录面貌，同时右上角出现一个窗口，该窗口提示enter键是显示滤波前与滤