1、浅层地震课程教案 第一部分 教学大纲一、本课程的适应层次、专业与参考学时 浅层地震方法及应用是适用于地球探测与信息、矿产普查与勘探、油气田开发工程及其它相关专业本科和研究生的一门选修课程。该课程参考学时74学时,其中实验实习10学时。二、课程的性质、目的和基本要求浅层地震方法及应用是适用于地球探测与信息、矿产普查与勘探、油气田开发工程及其它相关专业的一门选修课程。它是根据地球物理学原理来解决矿产、工程、石油、环境及国防等领域中存在的地质问题的一门学科,是地勘、工程、水文及环境勘查的重要方法之一。通过本课程的学习,要求学生了解本课程的基本原理和野外工作方法技术,特别在解决工程与环境地质问题时的抗
2、干扰、提高分辨率的措施,提高学生应用地震新方法解决相应地质问题的实际能力。这就必须加强实验实习课的份量,跟上目前工程建设和环境地质灾害探测对人才的要求。三、本课程的基本及重点内容根据教学计划要求和课程性质,本课程的教学内容主要包括:物探分类、浅层地震特点及发展史;弹性波的基本理论;地震波的时距曲线;地震勘探野外数据采集技术及方法;地震数据处理技术;地震资料的解释与应用;工程地震新技术。本课程的重点是地震勘探野外数据采集技术、处理技术、解释技术及工程地震新方法。地震勘探方法始于上世纪40年代处,经过半个多世纪的深入研究和运用,地震勘探技术得到了长足发展和进步。目前,它已成功应用于矿产资源勘探、石
3、油勘察、工程地质调查、水文地质调查、环境地质调查等诸多领域。涉及到矿产、石油、交通、建筑、考古、国防等与国计民生密切相关的各个行业。通过本课程的系统学习,学生应掌握地震勘探的基本原理、数据处理技术及解释方法,在此基础上进一步拓宽,适当了解近年来发展起来的一些新的相关技术和方法,争取做到学生走出校门后,可以独立从事施工设计并带队从事野外勘探、资料处理和成果解释的水平。为我国国民经济建设输送短缺人才。本课程的难点在于该学科涉及的范围很广,学生需要有扎实的数理基础和较好的计算机知识来理解和学习该课程。四、课程主要内容及学时分配本课程总学时74学时,其中理论教学64学时,实验实习10学时,实验、实习预
4、计安排校内进行。学时安排如下:浅层地震学时分配表章次内容合计讲授学时分配实验教学实习机动绪论22一弹性波的基本理论1414二地震波的时距曲线44三地震勘探野外数据采集技术及方法16142四地震数据处理1010五地震资料的解释与应用14104六地震新技术14104合 计746410注 课外作业不计入教学时数,任课教师可酌情增减。 实验课内容初步拟定为:实验课1:地震仪认识与操作实习,掌握地震仪器参数选择,和各功能键的使用实验课2:对折射波进行识别,绘制时距曲线,掌握从地震记录图上读取初至时间的方法,并进行界面速度和厚度的计算。实验课3:反射波测桩实验或瑞雷面波及折射波资料自动解五、考核办法 本课
5、程采用期末考试和平时作业成绩相结合的办法评定学生的学习成绩,其中期末考试占总成绩的70,平时作业成绩占总成绩的30,平时作业将在每一节课程之后布置给学生。六、教材及参考用书1教材:工程与环境地震勘探技术王俊茹编著,地质出版社,2002年。浅层地震勘探应用技术王振东编著,地质出版社,1988年。2参考书目应用地球物理教程何樵登主编,地质出版社,1991年。勘探地震学美R.E.谢里夫 加L.P.吉尔达特编,石油工业出版社,1999年总的来讲,所选教材是针对“工程探测”的,这与我们的培养方案是紧密相扣的。目的是要让大家对相应的、必须的、常用的方法技术有较深入的了解。在方法原理、理论、数理方面,难度不
6、是太大,比较浅显,但希望大家认真对待。考试将安排在课表上安排的本门课程的最后两个学时,采用开卷笔试形式进行。第二部分 课程具体讲授内容绪 论一、 学习意义l 热门课程:这是由国家经济建设的需要决定的。由于国家经济建设和社会持续发展的需要,蓬勃发展起来的工程与环境地震勘探技术与国民经济的方方面面密切相关,使之成为热门研究课题致意。如:高层建筑、重型厂房、桥梁、隧道、机场、水坝、电站、港口、码头、高速公路、铁路及其它巨型工程建设项目的地质基础资料的调查。勘探目标包括:岩溶塌陷、表层软泥层、地面沉降、地裂缝、岩体滑坡、泥石流等,甚至于煤层储量、采矿及隧道掘进过程中的超前探测等。另外,由于人们环保意识
7、的加强及重视,如对地下隐蔽污染源、地质灾害、考古等方面的探测不能靠传统的挖掘或钻探取样方法来解决,而应发展新型的无损探测技术,地震勘探就是其中的重要技术之一。总的来说,一门课程必须有需求,才有设计的必要,才有生源,才能有所为,才能有发展。二战时,很多从事海军声纳技术的专家,由于战争停止而转向石油勘探,从而发展了地震勘探技术。现在,它已发展为与国民经济的众多领域密切相关的一门学科。l 是重要的物探方法之一:对于地质或物探专业的学生来讲,地震方法是常用的、也是行之有效的方法之一。近年来,高分辨率地震勘探、VSP(Vertical Seismic Profile,垂直地震剖面)、瞬态瑞雷波法、反射波
8、测桩技术、横波勘探技术及常时微动观测技术等迅猛发展,自会有它们发展的必然性。另外,从国家对物探工作的总体投入来看,地震占到了95%,也可以说明问题。因此,希望在坐的各位有必要掌握好这们方法技术。二、 地震勘探方法概述1. 地球物理勘探(简称物探)方法是以物理学原理为基础,利用电子学、计算机数字处理技术、信息论等科学领域中的新技术,所建立起来的一整套勘探地下矿产或结构的方法。借助于各种物探仪器,根据地质任务,采用相应的观测方法(装置或系统)来观测地下岩石的各种物理响应,从而推断和解释地下岩石的构造特点、岩石性质等等,以达到勘察地下矿产或结构的目的。物探方法主要有: 重力勘探利用岩石的密度差异;
9、磁法勘探利用岩石的磁性差异; 电法勘探利用岩石的电性差异; 放射性勘探利用岩石的放射性差异; 地震勘探利用岩石的弹性差异。【注意:占5%的投入;占95%的投入】2. 地震勘探是物探方法中的重要方法之一。它是以地质介质的弹性为基础,利用地震仪器观测地震波信号,结合测区的地质、钻探以及其它已知的物探资料,对观测记录进行处理、解释,并作出地质推断的一门学科。 应用范围:在油气、煤田勘探、工程地质勘察等领域发挥着重要作用,已成为最为有效的勘探方法。一般来讲,地震勘探较其它物探方法的勘探精度高。缺点是:野外作业及室内资料处理工作复杂,工作量大,因而成本高。 方法分类: 观测系统(以反射波法为例): 激发
10、方式:炸药、锤击、震源枪、电火花等。 波型: 纵波(P波、压缩波、胀缩波) 横波(S波、剪切波,分为SH波和SV波) 面波(Rayleigh波、Love波)3. 浅层地震 主要方法:反射波法:不受速度场限制;主要应用于中、深层;分辨能力高。折射波法:条件是下部地层波速大于上部地层波速;主要用于表层;分辨能力相对较低。 面波法:条件是自由表面附近介质非均匀;主要用于表层分层及计算表层介质的横波速度;分辨能力较强。 应用:在工程地质勘察中,工程地震勘探的工作量已占工程物探总量的65%以上,主要用于:地矿、能源、水利、铁道、交通、城建、考古、环境、国防等。 主要解决的常见地质问题:A. 目的层埋深、
11、厚度(包括风化层);B. 断层、破碎带、洞穴、地下埋设物(包括考古、城市地下管网等);C. 矿产资源(煤、砂体、石油、);D. 地质灾害(由波速地压或地应力的变化);E. 地基、桩基、坝体质量评估(面波或频率探测)。三、 地震勘探方法的发展历程1. 地震勘探历史A. 采集记录方面:模拟信号数字信号(提高动态范围,室内可处理)B. 仪器方面: 单道多道; 信号有线传输无线传输(采集站); 仪器功放检波器集成功放+仪器功放。C. 数字处理方面: 单道 水平叠加 叠前处理力度加大(多道滤波) (AVO)(偏移技术) (叠前偏移) 弹性介质模型双向介质粘弹性介质各向异性介质研究 构造功能构造+岩性功能
12、(油藏描述)2. 勘探方法A. 针对应用领域: 石油地震 煤田地震(煤层储量、断层、瓦斯超前预报) 工程地震(表层结构、承载力)B采集方面: 单道二维(2D)三维(3D)四维(4D)C针对波型:D针对观测系统: 地面观测:反射、折射、面波 井中观测:透射波法(CT) VSP:零偏非零偏Walk Away第一章 弹性波理论基础地震勘探是研究人工激发的弹性波在岩石中的传播规律,以探测地下地质构造或地质赋存状态的地球物理勘探方法。因此,可将弹性力学的基本理论引入到地震勘探中来。第一节 弹性理论概述一、基本概念l 弹性:物体在外力作用下会发生形变。当外力消失时,由于物体具有阻止形变的内力,它会使物体恢
13、复原状,物体的这种属性被称为弹性。l 塑性:如果外力超过物体的弹性极限,或外力作用时间太长,当外力消失时,物体不能恢复原状,物体的这种性质被称为塑性。l 观点1:弹性与塑性是两个极限概念,是为了解决或讨论问题的方便常常要作的假设。多数情形下,介质的性质介于二者之间,没有绝对的弹性和塑性。l 观点2:一种材料或物体的弹、塑性与内因和外因两方面的因素有关。内因:材料本身的性质;外因:作用力的大小,作用时间的长短。l 弹性介质:产生弹性形变的介质。l 地震弹性波:在弹性介质中传播的地震波。l 弹性介质的几种类型: 根据弹性性质与空间方向的变化关系分:A. 各向同性介质:弹性体的弹性性质与空间方向无关
14、。B. 各向异性介质:弹性体的弹性性质与空间方向相关。 根据地震波速与空间位置的变化关系分:A均匀介质:地震波速不随空间坐标变化而变化。B非均匀介质:地震波速随空间坐标变化而变化。C连续介质:地震波速在介质中连续变化。D层状介质:在非均匀介质中,介质的地震波速变化表现出成层性。(显然,当层状介质的层厚趋于零时就是连续介质)实际地下介质多表现出各向异性,但当横向变化非常缓慢时,可以忽略,同时,也是为了讨论问题的方便,可将其假设为各向同性介质。同样,实际介质多表现为非均匀介质,对于沉积岩来讲,同时,也为了讨论问题的方便,我们可将介质看成层状介质,且往往假设层内为均匀介质,当沉积旋迴非常发育时,层数
15、增多,层厚变薄,此时的层状介质就过渡到连续介质。二、弹性系数弹性性质用弹性系数(或弹性模量)来度量。对于各向同性介质,弹性性质可由等5个弹性系数来确定,而其中独立的参数只有2个,即已知其中的2个,就可计算出另外3个,一般通常用拉梅()常数和剪切模量这两个系数来表述。对于各向异性介质,描述其弹性性质的弹性常数有21个。 (一)杨氏模量() 定义:表示物体抗拉伸或挤压的力学参数。越大,抗拉伸或挤压的阻力越大。(二)剪切模量()定义:是表示物体阻止剪切应变的力学参数,单位与应力相同。越大,切应变越小。液体中=0。 (三)泊松比()定义: 反映的是物体的横向拉伸(或压缩)对纵向的压缩(或拉伸)的影响,
16、越大,影响越小。 一般未胶结的砂土较高,而坚硬岩石的较小。 自然界中,值的取值范围,当=0.25时的介质被称为泊松固体。 流体的=0.5。 (四)体变模量或体积模量()定义: 表示物体的抗压性质,所以又称为抗压缩系数。 (五)拉梅常数()定义: 我们知道,在弹性力学中,对于弹性介质,应力与应变满足虎克(Hook)定律:(条件:均匀各向同性介质)式中: 此时,上述5个弹性模量的数学意义可表示为:5个参数之间的关系:三、波动方程波动是指弹性体内相邻质点间存在应力梯度时产生质点的相对位移。波动方程便是研究应力不平衡时的状态方程。由固体弹性理论可知,在均匀各向同性介质中,三维波动方程为:该方程又被称做
17、Navier方程。式中:为外力;是介质质点受外力作用后产生的位移向量;统称为拉梅常数;为介质的密度;标量是体变系数,。l 基础知识:梯度:散度:旋度:对于Navier方程两边分别取散度(div)和旋度(rot),则可得:式中:;为标量,表示的是一种胀缩力(因为,当时);为矢量,表示的是一种旋转力(因为)。上式表明,介质在两种不同得力得作用下,产生两种不同性质的扰动。在胀缩力的作用下,将产生仅由体变系数决定的介质体积相对胀缩的扰动,称为纵波;而在旋转力作用下,则产生由决定的角度转动的扰动,称为横波。纵波又称胀缩波,横波又称剪切波。它们分别以速度传播。即:一般地,位移矢量和力矢量既存在散度,也存在
18、旋度。因此,可表示为:式中:分别为位移场的标量位和矢量位;分别为力场的标量位和矢量位。此时,用位函数表示的波动方程可表示位:在远离震源区,震源力已不复存在,此时,上述方程便变成以下齐次方程:【这表明:当我们今后用上述方程去描述波的传播问题时,所求解出来的并不是位移,而是位移的位函数】第二节 弹性波的形成在外力的作用下,弹性介质中存在两种扰动。胀缩力的扰动对应,即介质中产生了体积形变,体积形变的传播形成纵波;旋转力的扰动对应,介质中质点产生了旋转形变(切应变),切应变的传播形成横波。它们都属于体波。此外,还有沿自由表面传播的面波。下面将讨论各种波的形成及传播特点。为考虑问题方便,假设介质为均匀各
19、向同性介质,则其弹性参数具有球形对称性。因此,可用球面坐标系来讨论问题。球面坐标系与直角坐标系的关系为:在图中:r:球面上任意一点P到直角坐标原点O的距离;:r所在圆面与X轴的夹角;:OP与Z轴的夹角。则有:一、纵波(P波、胀缩波、疏密波、压缩波)在球面坐标系()中,纵波的波动方程可表示为:令:,则得其达朗贝尔解:式中,为两个任意常数,它们与震源得性质有关。若令:则表示了波动在某一个固定时刻得状态,即波阵面。上式中,第一项表示波动随时间增加远离震源方向传播,第二项则相反。所以前者被称为发散波,称后者为会聚波。会聚波不符合实际情况,是物理不可实现得,因此,不予考虑。此时,其达朗贝尔解可进一步修改
20、为:相应的位移为:式中:为震源强度函数;为单位向量。根据纵波位移表达式,可以总结出纵波具有以下特点:(1) 质点位移方向与的方向一致,即质点振动方向与波的传播方向一致。因此,纵波又称为线性极化波。由于有正有负,因此在纵波扰动带内将会间隔出现膨胀带和压缩带,所以又称为疏密波或压缩波;(2) 纵波的传播速度:(3) 当一定时,质点位移决定于,即质点位移的大小决定于离震源的距离、震源强度函数及其变化率;(4) 振动强度随传播距离增大而减小,这一现象被称为球面扩散。二、横波(剪切波、S波,包括SH波和SV波)横波由旋转力产生。位移位函数所表达的波动方程为:在球面坐标系()中,其解为:因此,横波具有以下
21、特点:(1) 传播速度:(2) 横波在传播方向上的质点位移,而在垂直于传播方向的和方向上具有位移,说明横波质点振动方向与传播方向正交。横波也是线性极化波。振动方向为水平时称SH波,振动方向为垂直时称SV波。(3) 振动强度决定于旋转激发力的强度函数及其变化率。(4) 横波同样具有球面扩散特征。(5) 在液体和气体中,由于,所以不存在横波。(6) 根据可以求泊松比:由此可见,的取值范围为。一般岩石可看作泊松固体,即,此时,可见横波速度纵波速度,即横波传播得比纵波慢一些,这就是为什么当地震发生后,往往会感觉到两次强烈震动得原因。(7) 只要测出,便可计算其它弹性摸量值:用这种方法测得的弹性摸量称为
22、动态弹性摸量,二实验室内测得的则称为静态弹性摸量。三、面波(Rayleigh波、Love波)除体波外,还有另一类波,从能量上讲,它们只分布在弹性分界面附近,因此,统称为面波,如Rayleigh波和Love波。A. Rayleigh波(又称地滚波)这是一种在自由表面(空气与地球表面形成得弹性分界面)附近大约一个波长范围内传播的波。这种波是英国学者Rayleigh于1887年在理论上首先证明出来而得名。它具有以下特点:(1) 质点振动轨迹沿与波传播方向成反方向的椭圆轨道运动,即沿逆时针方向的椭圆轨迹运动,因此它是椭圆极化波;(2) 这种椭圆轨迹是由相位相差的两个相互垂直的振动分量合成而得;(3)
23、振动能量沿垂直方向衰减快,而沿水平方向(近似于传播方向)衰减慢。因此只在地表附近一个波长范围内传播。由于体波球面扩散时,而Rayleigh波,所以,在远离震源时,面波能量往往强于体波能量,这就是称之为地滚波得由来;(4) 传播速度:由Rayleigh方程求解该方程得:例如:0.250.330.400.500.920.9330.9430.956一般土的泊松比,因此,可以近似看作,即可通过计算面波速度而近似获得横波速度。(5) 存在频散现象,即波速是频率的函数。利用这一特征,可以从事工程勘察,即面波勘探。B. Love波是一种类似于SH型的面波,质点在水平方向振动并垂直于波的传播方向,存在于表层和
24、地下层界面附近。其形成条件是:上部地层下部地层这种波的能量往往比较弱(当用垂直检波器检测时),因此很少讨论。于之相反,而Rayleigh波的振动椭圆轨迹的长轴在垂直方向,所以,垂直检波器很容易检测到它的存在。四、有效波和干扰波有效波工作时用来解决测区地质任务的波;干扰波有碍于提取并识别有效波的其它波。l 有效波与干扰波的概念是相对的。l 但有些干扰波(如随机干扰、工业电干扰等)绝对属于干扰波,是必须去除或压制的 ,从而发展了很多相应的数字处理技术去压制它们。l 信噪比(S/N)定义为:有效信号强度/干扰信号强度。第三节 弹性波的描述一、地震波场的基本特征地震波场的基本特征包括:l 运动学特针地
25、震波传播的时间(旅行时或走时)与空间位置之间所展现出的特征。主要与波速及地下构造因素有关。l 动力学特征地震波传播时,其波形、频率、振幅随空间和时间变化而变化所展现出的特征。主要与激发源、介质吸收及构造形态等有关。假设地下介质是弹性介质,地震波场即可视为弹性波场。因此,下面从更广义的角度来表述弹性波场描述所需的几个基本概念。二、振动图和波剖面图波在弹性介质中传播时,既是时间t,又是空间r的函数,即位移。1. 振动图对于固定的r,即只观察某一固定的质点,其振动随时间t的变化关系,即,该函数的图象就称为振动图。对于简谐振动来说,如我们知道:表示振幅;称为周期;表示频率等。对于一般意义的弹性波,由于
26、介质的粘滞性等原因,其频率常表现为:由于它不是精确的周期振动,所以:视振幅,它描述振动的强弱视周期,它描述振动的快慢;视频率初至时间,它表示某质点刚刚开始振动的时间延续度,它表示质点从开始振动到停止振动的时间间隔。它直接影响到地震勘探的分辨能力(因此,要压缩它,即发展了所谓的反褶积技术或去子波技术)。2. 波剖面图对于固定的时间t,各质点的位移的图象就称为波剖面图。它是表示在某一时刻各质点之间振动强弱关系的图象。视波长,即相邻两个波峰之间的距离;视波数,即单位距离内视波长的数目。为表示方便,以后取消“视”或“”号。关系式:,表示视速度,表示视周期。地震传播既表现出运动学特征,又表现出动力学特征
27、。我们知道:A. 运动学特征即描述波在传播过程中,波前、射线、等时面、视速度等参数随空间变化而变化的关系。即用波剖面图来表述;B. 动力学特征即描述波在传播过程中,波形、波幅或波谱(频率)、能量等参数随时间变化而变化的关系。即用振动图来表述。3. 波前面与等相位面以均匀各向同性介质中的点震源为例,其球面波传播示意图为:对于某个时刻,介质振动将形成三个区域:振动结束区域、扰动区域、振动未到区域。波前面(波前)某一时刻介质中刚刚开始振动的质点所连成的面。波尾面(波尾)某一时刻介质中刚刚停止振动的质点所连成的面。振动带波前与波尾之间的介质区域。此时,其中所有质点正处于振动状态。等相位面在某一时刻,具
28、有相同相位状态的质点所连成的面。显然,波前面和波尾面都是等相位面。随着球面波前的进一步扩大,当扩大到非常大时,可以把球面上的局部近似看成平面。这就是到以后可将远离点震源的波作为平面波来研究的原因。4. 时间场与等时面波前面上的每一点都有与之对应的旅行时间t,因此,旅行时间t也是空间的函数。这种反映波场面时空关系的物理场称为时间场,即。在地震波传播区域内,若将旅行时间相同的各点连成曲面,就构成时间等值面(简称等时面)。等时面与波前面形状是相似的,都表示在某一时刻刚刚开始振动的点所连成的面,但二者的含义不同。等时面为时间量纲,而波前面为空间量纲。射线是用来表示波的传播路径与方向的几何线。射线原理(
29、费马原理)波沿射线路径传播的时间要比沿其它路径传播的时间小,即射线肯定是垂直波前面的。因此,射线原理也称为最小时间原理。由于波前面的形状与等时面的形状相一致,因此,射线也垂直于等时面。l 思考:在均匀和非均匀介质中,射线和等时面的形状分别是什么样?5. 地震波的频谱振动图是地震信号随时间t而变化的图象,可用“时域”函数f(t)表示。有时,为了讨论问题的方便,在地震资料处理中,需要在“频率域”来研究地震信号的特点。如反射波、面波、工业电干扰等,虽然在时域内它们是交织在一起的,但在频率域内,它们便可以较好地区分。因此,我们有必要在此讨论一下地震波的频谱。设所对应的频谱为,则为一复数函数。由付氏变换
30、的单一对应关系我们知道,二者是可以互换的,即:由付氏逆变换可见:地震信号是由一系列频率的简谐振动复合而成的。换句话说,不同频率的简谐振动可以复合成非常复杂的地震信号。可用复变函数来表示。式中:振幅谱相位谱所以,该式可以更清楚地表明:信号可以看作由无限个不同频率、不同初始相位的简谐振动的叠加。各个频率的简谐振动的振幅就构成了振幅谱;各个频率的简谐振动的相位就构成了相位谱。如:振幅谱和相位谱合起来,统称为“频谱”。“频谱分析”将时域地震信号通过付氏变换得到其频谱曲线,进而进行分析的过程。通过频谱分析,我们可以看出各种频率成分对地震波能量的贡献的大小。其中振幅最大的频率称为“主频”。主频频带(大部分
31、能量集中于此)频带宽度在地震勘探过程种,我们所记录到的各种类型的波,其主频及频带是有所区分的。如:因此,通过设计合适的滤波器,就可达到压制不期望的波(干扰波),保留期望保留的波(有效波),从而提高地震信号的信噪比的目的。这也是频谱分析的出发点之一。以下是一些典型(或理想)波形函数的频谱。 c从上图可以看出,只要时域信号具有轴对称性,即,则其相位谱恒等于,即。证明如下:式中:(振幅谱)(相位谱)又因:由于为偶函数,而sin为奇函数,可见:,从而。第四节 弹性波的传播质点振动是产生波动的根源,是基础,而波动是振动的传播过程。其形成机理:A点振动,势必对B点形成作用力,使B点振动;同时,B点有阻止A
32、点振动的反作用力,使A点再回复,向相反方向振动,到一定时候,C点对A点形成作用力,使A点在其平衡位置来回振动。质点间振动相互传递,由近及远,从而形成波动。一、波动传播原理(一) 惠更斯原理在弹性介质中,已知时刻的波前,求时刻的波前时可将原时刻的波前面上的每个质点看作新的震源,新震源以传播形成一系列圆圈,所有这些圆圈的包络就是时刻新的波前的位置。 惠更斯原理同时适用于均匀和非均匀介质。尽管如此,但在作图时应注意:对于均匀介质,所有子圆的大小都相同;而对于非均匀介质,由于时刻的波前面上的每个质点处的速度不一样,所以,各子圆的大小也将发生变化。(二) 费马原理(最小时间原理)波动总是沿传播时间最小的
33、路径传播,这些路径就是射线。在均匀介质中,射线为直线;在均匀层状介质中,射线满足斯奈尔定律,即波动沿满足斯奈尔定律的路径传播时所用时间最小。同时,折射波的存在也有力地说明了费马原理。(三) 视速度定理真速度V 沿射线方向估计出的波的传播速度。出射角 射线与地面法线之间的夹角。视速度 由于观测方向偏离射线方向,此时,由观测数据估计出的波的传播速度将不等于真速度,而称为视速度。V。视速度定理:二、地震波反射、透射和折射波的形成在弹性分界面处,产生反射、透射波是波动的共性。本节讨论地震波在什么条件下才会发生反射、透射及折射,其能量的变化和传播规律等。(一) 斯奈尔定律(Snells Law)Snel
34、l定律的描述: 入射线、反射线、透射线在同一平面内(即射线平面); 入射角=反射角 透射角取决于入射角和界面上、下介质的波速比值,即:即:,P为射线参数,用于确定反射线与透射线轨迹的参量。另外请注意:斯奈尔定律也称为反射透射定律;此时的透射就是光学中的折射。关于斯奈尔定律的证明参见教材中P.14。(二) 垂直入射(或法向入射)时波的反射和透射假设:A.各向同性弹性介质; B.波通过分界面时满足两个连续性条件: 弹性位移的法向分量和切向分量连续; 应力的法向分量与切向分量连续。依据:A.斯奈尔定律; B.地震波可分解为无限个简谐振动。由此可推导出用于揭示入射、反射、透射波能量分配关系的佐普瑞兹(
35、Zoepritz)方程。在该方程中,令入射角,并求解得:式中:反射系数;透射系数;反射波振幅;透射波振幅;入射波振幅;上、下介质的密度;上、下介质的速度。从上式可以看出:即入射波在经过弹性分界面时,能量只分解为反射波和透射波。另外,还可看出: ; 可0,也可0,与上、下介质的弹性参数有关:若,则,此时,反射波的质点振动与入射波的质点振动有相位延迟,延迟量为,这种现象称为半波损失。我们把: 叫做波阻抗; 具有波阻抗差异的分界面称为反射界面或反射面。因此,反射波形成的条件是:界面上、下介质必须存在波阻抗差异。讨论: 若,则即使存在地质界面,也无法通过反射地震波法观测到反射信息。这就是往往地质界面与
36、地震相(或地震界面)有时不一致的客观原因。 若,则将产生强反射。但此时,表明透射能量非常弱,这种现象不利于利用反射地震波法观测其下部介质的信息。此时,我们称该层为高阻屏蔽层。(三) 非法向入射时波的反射和透射 将产生四种波型:反射纵波、反射横波、透射纵波、透射横波。其射线仍满足斯奈尔定律,即: 设为的反射系数;为的反射系数;为的透射系数;为的透射系数。则其能量分配满足以下方程(在法向入射时提出的前提条件下得到的):写称矩阵方程:这就是著名的Zoepritz方程。若令,则根据斯奈尔定律有,我们便可得到前述法向入射时纵波的反射系数和透射系数的表达式,而且有。因此,如果纵波法向入射,则只产生反射纵波
37、和透射纵波,它们是同类波。若非法向入射,则除产生反射、透射纵波以外,还将产生反射、透射横波。因为它们与P波类型不同,所以又叫转换波(又称P-SV波)。若已知介质弹性参数及入射角,则可根据上述Zoepritz方程计算出P、S波的反射系数、及透射系数、。同时,我们可以看出,当地下介质弹性参数一定时,随着入射角发生变化,反射波的能量也将发生变化,即反射强度将随入射角的变化而变化。下图分别是对应和两种情形的反射、透射波之间的能量分配曲线。根据能量守恒原则有:由于入射角的变化实际上是由于地面上炮检距的变化造成的,因此,反射强度的变化又称AVO(Amplitude Versus Offset, 振幅随炮检
38、距变化而变化)。由此发展了所谓的AVO技术:即如果我们能从观测记录中找出反射波振幅随炮检距的变化关系曲线,那么,我们便可根据Zoepritz方程反演出地下介质的弹性参数等,从而得到所谓的“P波剖面”、“密度剖面”、“速度比剖面”、“泊松比剖面”、“横波剖面”、“振幅比剖面”等。通过大量的模型计算可以发现:上、下介质的速度差异对反射强度变化的影响明显,而密度的变化对反射波强度变化的影响相对较小。(四) 折射波的形成与传播1. 折射波的形成条件:V2V1根据斯奈尔定律:若,则一定有一个入射角,使得。令,则: 此时的入射角称为临界角。当入射角达到临界角时,透射波将沿界面滑行,产生滑行波。 当入射角(
39、临界角)时,透射波转化为滑行波。滑行波沿界面滑行时,界面上的每个质点都相当于新的扰动源(根据惠更斯原理),新的扰动源产生的波再到达地面,被布置在地面上的检波器检测到,此时记录下的波即为折射波。由于折射波沿界面滑行时,是以下伏介质中的波速传播的,当炮检距比较大时,它可能先于直达波到达检波器,所以有时又叫“初至”。2. 折射波的接收条件:xxm只有当时,才能接收到折射波。对于水平层状介质,式中:。在区域内,无法接收到折射波,因此该区域被称为“折射波盲区”。影响折射波盲区大小的因素:A. 界面上、下介质的速度比,越大,盲区越大;B. 折射面的埋深h,h越大,盲区越大;C. 折射面的倾角(如下图),显
40、然,。即在折射面的上倾方向接收折射波比在折射面的下倾方向接收折射波,其盲区范围要小。当折射面倾角增大到时,即在下倾方向就再也无法接收到折射波了。三、地震波的绕射和散射地震波的绕射与光学、声学中的衍射现象雷同。地震光学、声学绕射衍射透射折射当地震波遇到几何尺寸与其波长基本相当的不均匀地质体时,如断层错位造成的棱、角点(简称断点)、地层尖灭点、洞穴等往往都能满足绕射波的形成条件。因此,绕射波是识别断层、地层尖灭点、洞穴等突变地质体的重要标志。绕射波在水平叠加记录上的形态是一条规则的双曲线,因此,在剖面上很容易识别。对于任意点C,它离开震源的水平距离为x绕射现象同样可用惠更斯原理来解释。事实上,界面
41、上的每个点都可看作是新震源,在地面上观测到的一个连续界面上的反射波是由所有这些新震源产生的绕射波的叠加。此时,我们称之为“广义绕射”,由于地质体不连续造成的间断点而引起的绕射波被称为“狭义绕射”。如不特殊申明,一般所说的绕射是指“狭义绕射”。有时,地层由于侵蚀,造成层面不光滑,粗糙不平,当地震波传播到这种凹凸不平的界面时会产生散射现象。散射也称为漫射,将造成地震剖面中反射面不连续,断断续续,时隐时现的情况。第五节 弹性波的衰减概念:能通量密度(波的强度)单位时间通过垂直于波动传播方向(即波前面)单位面积的能量。用I表示。一般,即I与振幅A的平方成正比。地震波在传播过程中,其能量的变化由多种因素
42、引起。主要有:波前扩散、介质吸收造成能量衰减、分界面造成反射、透射及层间多次反射。一、波前扩散(球面扩散、几何扩散)以球面波为例,随着波前面的不断扩大,通过单位面积的能量将逐渐变小。设激发总能量为E,则:考虑到,所以,这表明,球面波的强度与传播距离的平方成反比,振幅与传播距离成反比。同样,可以推导出柱面波、平面波所对应的通式:式中,分别对应平面波、柱面波、球面波的扩散。二、吸收衰减介质对地震波的吸收实际地质介质并非完全弹性介质,弹性波在非完全弹性介质中传播时,介质中质点振动的能量因质点之间相互摩擦,有一部分能量要转化为热能而损失掉,这种现象称为介质对地震波的吸收。粘滞介质具有吸收性能的非理想弹
43、性介质。或叫“粘弹性介质”。吸收衰减可用以下方程表示:式中,为处的振幅;为吸收系数,表示单位距离内振幅的衰减率。吸收系数是频率f的函数,即。 当频率很低时,满足,则:式中,为粘滞系数。 当频率教高时,满足,则:由此可见,f越高,越大。所以,大地具有低通滤波作用。同时,波的传播速度在频率比较高时,亦会随f变化,这种现象称为“波的频散特征”。一般在地震勘探所能记录到的有效地震信号的频带内,体波的频散现象很弱,可以忽略。吸收的强弱随岩石的性质变化而变化。一般地,疏松岩石(如风化层)比硬岩石的吸收要大。因此,地震波在通过表层风化层时,会大大地被吸收掉,尤其是高频成分,很快被吸收掉,导致地震勘探的分辨能
44、力下降。为尽量避免这一现象发生,人们往往将激发点置于风化层以下,并将检波器挖坑埋置。最后,综合球面扩散和地层吸收:三、透射损失为考虑问题简单,只考虑一次反射,未考虑多次反射。式中,表示由界面反射回来的波经过由到介质中的界面时的透射系数。若假设入射角,则依次类推:这表明,弹性波每透过一个弹性分界面时,必使入射波的振幅乘以一个因子,这个因子称为透射损失因子。由于,所以,入射波的能量由于透射损失会耗损一部分。但由于,所以这部分损失往往被忽略。第六节 地震反射波记录道的形成设表示入射波的振幅,表示吸收系数,r表示波的传播距离,表示地震子波(它由震源函数和介质性质决定),则由地下某一反射界面反射到地面的波的位移为:若令则1. 地震记录的形成:就是将地下所有反射界面所获得的各位移量相加形成的振动图。绘图表示如下:对应各界面的波形受多种因素影响,如:、等;振幅受r、的影响;极性受反射系数的影响。若各界面反射项的子波长度(为相邻两层之间的双程旅行时),则叠加后,互不干涉,可以分开,即可以分辨各层的存在。反之,则无法分开,从而无法分辨各层的存在。在地震勘探中,将分辨率定义为T/4,即由相邻层之间的厚度所引起的双程旅行时至少应大于或等于T/4,最好是T/2,便可分辨。即:其中f为地震波的主频。 优势波由强反射系数界面所引起的反射记录中