第1章地震偏移成像基础
- 格式:doc
- 大小:1.30 MB
- 文档页数:29
方法原理地震勘探中的反射共偏移方法,利用该方法进行有效的相位对比与追踪可获得反射界面的位置及厚度关系。
现对其方法原理作以简单介绍。
反射共偏移法又称陆地声纳法,依据反射波勘探原理,在单边排列的基础上选定最佳偏移距,即最佳反射窗口,采用单道或多道叠加小步长顺移前进观测系统(图1)。
“偏移成像”功能指在给定速度等参数后将地震时间剖面转换成空间剖面一种数据处理技术。
它最大限度地将反射同相轴归位到空间反射点上,同时消除了由于倾斜界面、尖灭点等引起的反射相位“偏移”现象,因而有利于地震剖面的解释和应用。
它所应用的原理就是多次覆盖原理。
所谓多次覆盖就是共反射点水平叠加,简称水平叠加,它是20世纪60年代初期发展起来的地震勘探方法。
它的实质是对反射界面上同一反射点进行重复观测,而激发点和接收点是在不断改变。
从而达到增强有效波、压制干扰波。
单边排列观测系统(如图1所示)指仅在接收点排列一侧激发的观测形式。
设某一单边排列接收道数为R 、道间距为I 、偏移距为O 、移动步距为P ,由几何地震学知当界面水平时其反射段长度为:2)1(-*=R I L 当整个排列移动步距P 小于反射段L 时出现反射段重复即多次覆盖。
图1 单边排列多次覆盖示意图资料的处理与解释反射共偏移探测数据在自行研制开发的KDZ2.8软件平台上进行,其中时间域里主要处理过程包括:信号录入、格式转换、预处理、数字滤波、修饰处理和偏移剖面形成与显示等内容,其中预处理包括道集重排、振幅平衡、静校正、二次采样等,修饰处理包括空间混波等。
处理的结果是由获得偏移时间剖面,根据剖面中反射相位同相轴的连续追踪与对比,结合已知地质资料及地质体的各种特征进行解释,最终形成地质剖面。
不同地质体在时间剖面中具有不同的反射波同相轴特征,即反射时间不同。
其界面的具体位置要根据每一组反射波旅行时间进行深度计算。
1、偏移技术分类【叠前/后偏移】可根据不同的标准对目前的地震偏移成像技术进行简单分类:按照所依据的理论基础,可以分为射、线类偏移成像和波动方程类偏移成像;根据输入数据类型,可以分为叠前偏移和叠后偏移;根据实现的时空域,可以分为时间偏移和深度偏移;按照维数,可以分为二维偏移以及三维偏移等;1.1叠前偏移使CSP道集记录或COF道集记录中的反射波归位,绕射波收敛。
●叠前偏移有椭圆切线法【手工方法,不适用】、Rockwell偏移叠加法【波前模糊法的拓展,计算量也很大】和Paturet-Tariel偏移叠加法【为了进行偏移,我们应当把的曲线上的地震能量(即采样点振幅)送到零炮检距绕射双曲线的顶点M上去叠加。
这样, 把各个相同炮检距的剖面偏移后叠加在一起即得偏移叠加剖面】等1.2叠后偏移基于水平叠加剖面,采用爆炸反射面的概念实现倾斜反射层归位和绕射波收敛。
●叠后偏移有波前模糊法、绕射曲线叠加法【两种方法原理简单,都是基于惠更斯原理提出的,前者将一个道上的波场值送到各个道上去叠加—输出道法,后者把各个道上的相应值取来在一道上叠加—输入道法,但是计算量很大】2、偏移成像特点●具有地震勘探本身的特征●计算机使其研究由地震波运动学特征过度到地震波动力学特征●提高地震空间分辨率和保真度●偏移成像是使反射界面最佳成像的一种技术●处理反射波,使之成为反映地下界面位置和反射系数值的反射界面的像3、偏移成像原理图偏移过程定量分析【Chun and Jacewitz ,1981】2(tan )/4t dx v t θ=221/2{1[1(tan )/4]}t dt t v θ=--221/2tan tan /[1(tan )/4]t t t v θθθ=-3.1 偏移前后的图例4、偏移方法分类5、实际中应用的一些偏移算法5.1 Kirchhoff 积分法【波场外推】适用条件:只满足均匀介质的情况。
[]111'1111(,,,)'4S R u u u x y z t u dS vR n t n R R n π⎧⎫-∂⎡∂⎤∂⎡∂⎤⎡⎤⎛⎫⎡⎤'''⎡⎤=-+⎨⎬ ⎪⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥∂∂∂∂⎣⎦⎝⎭⎣⎦⎣⎦⎣⎦⎩⎭⎰⎰式中的[[u]]不再是推迟场,而是超前场。
1、偏移技术分类【叠前/后偏移】可根据不同的标准对目前的地震偏移成像技术进行简单分类:按照所依据的理论基础,可以分为射、线类偏移成像和波动方程类偏移成像;根据输入数据类型,可以分为叠前偏移和叠后偏移;根据实现的时空域,可以分为时间偏移和深度偏移;按照维数,可以分为二维偏移以及三维偏移等;1.1叠前偏移使CSP道集记录或COF道集记录中的反射波归位,绕射波收敛。
●叠前偏移有椭圆切线法【手工方法,不适用】、Rockwell偏移叠加法【波前模糊法的拓展,计算量也很大】和Paturet-Tariel偏移叠加法【为了进行偏移,我们应当把的曲线上的地震能量(即采样点振幅)送到零炮检距绕射双曲线的顶点M上去叠加。
这样, 把各个相同炮检距的剖面偏移后叠加在一起即得偏移叠加剖面】等1.2叠后偏移基于水平叠加剖面,采用爆炸反射面的概念实现倾斜反射层归位和绕射波收敛。
●叠后偏移有波前模糊法、绕射曲线叠加法【两种方法原理简单,都是基于惠更斯原理提出的,前者将一个道上的波场值送到各个道上去叠加—输出道法,后者把各个道上的相应值取来在一道上叠加—输入道法,但是计算量很大】2、偏移成像特点●具有地震勘探本身的特征●计算机使其研究由地震波运动学特征过度到地震波动力学特征●提高地震空间分辨率和保真度●偏移成像是使反射界面最佳成像的一种技术●处理反射波,使之成为反映地下界面位置和反射系数值的反射界面的像3、偏移成像原理图偏移过程定量分析【Chun and Jacewitz ,1981】2(tan )/4t dx v t θ=221/2{1[1(tan )/4]}t dt t v θ=--221/2tan tan /[1(tan )/4]t t t v θθθ=-3.1 偏移前后的图例4、偏移方法分类5、实际中应用的一些偏移算法5.1 Kirchhoff 积分法【波场外推】适用条件:只满足均匀介质的情况。
[]111'1111(,,,)'4S R u u u x y z t u dS vR n t n R R n π⎧⎫-∂⎡∂⎤∂⎡∂⎤⎡⎤⎛⎫⎡⎤'''⎡⎤=-+⎨⎬ ⎪⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥∂∂∂∂⎣⎦⎝⎭⎣⎦⎣⎦⎣⎦⎩⎭⎰⎰式中的[[u]]不再是推迟场,而是超前场。
论偏移成像论文摘要地震偏移成像技术是现代地震勘探数据处理的三大基本技术之一,主要包括射线偏移和波动方程偏移两大类,主要目的是实现反射界面的空间归为和恢复反射界面空间的波场特征、振幅变化和反射系数,提高地震空间分辨率和地震保真度。
偏移技术具有地震勘探本身的特点,但是地震偏移方法本身由于使用计算机而引起了许多革命性的变化。
这就使得它从研究简单的探测目标的几何图形进而发展成研究反射界面空间的波场特征、振幅变化和反射率等,在本论文中主要介绍地震偏移成像技术的基本原理,地震剖面的偏移和叠加偏移,叠前部分偏移。
时间偏移和深度偏移等方面来介绍。
正文一、偏移成像的基本原理在水平叠加时间剖面上显示出来的反射点位置是沿地层下倾方向偏离了反射点的真实位置的,这种现象就称为偏移。
反射地震方法是根据在地面上以一定方式进行弹性波激发,并在地面的一定范围(孔径)内记录来自地下弹性分界面的反射波来研究地下地质岩层结构及其物性特征的一种方法。
因此,也可以把它看作一种反散射问题。
就反射地震观测方式的特点,它的成像问题要分作两步,第一步是按照一定的方式记录到达地面的反射波,第二步用计算机按一定的计算方法对观测数据进行处理,使之成为反映地下地质分层面位置及反射系数值的反射界面的像。
而地震偏移技术就是在第二步过程使反射界面最佳地成像的一种技术。
地震偏移可在叠前做也可在叠后做。
叠前偏移是把共炮点道集记录或共偏移距道集记录中的反射波归位到产生它们的反射界面上并使绕射波收敛到产生它的绕射点上。
在把反射波回投到反射界面上和绕射波收敛到绕射点上时要去掉传播过程的效应,如扩散与衰减等。
最后得到能够反映界面反射系数特点的并正确归位了的地震波形剖面,即偏移剖面。
叠后偏移是在水平叠加剖面的基础上进行的,针对水平叠加剖面上存在的倾斜反射层不能正确地归位和绕射波不能完全收敛的问题,采用了爆炸反射面的概念来实现倾斜反射层的正确归位和绕射波的完全收敛。
地震偏移的部分类型见表1-1。
《地震勘探原理与解释》复习要点第一章绪论(不作为考试内容)第二章地震波运动学理论§2.1 几何地震学基本概念1、掌握基本概念,如地震子波、波面、射线、振动图、波剖面、视速度、视波长、全反射、雷克子波。
2、掌握基本原理,如反射定律、透射定律、Snell定律、惠更斯原理、费马原理等。
3、地震波的分类。
§2.2 常速单界面的反射波特征及时距关系1、基本概念:时距曲线、时距曲面、时间场、自激自收、共激发点、偏移距、初至时间、纵测线、同相轴、正常时差、倾角时差、动校正等。
2、基本原理:虚震源原理、讨论时距曲线的实际意义、直达波时距曲线及方程、反射波时距曲线及方程、反射波时距曲线的主要特点。
§2.3 变速多界面的反射波特征及时距关系1、基本概念:均匀介质、层状介质、连续介质、参数方程、平均速度、射线方程、等时线方程、回折波、最大穿透深度等。
2、基本原理:水平层状介质和连续介质情况下讨论反射波时距曲线的基本思路;水平层状介质和连续介质情况下反射波时距曲线的主要特点。
§2.4 地震折射波运动学1、基本概念:折射波盲区、初至波、续至波、交叉时、信噪比等。
2、基本原理:产生折射波的条件;利用折射波法研究地下地层起伏的基本依据;折射波与反射波的主要差异。
3、分析理解:单界面(水平和倾斜)直达波、反射波与折射波时距曲线之间的关系;三层介质情况下折射波的时距曲线及其特点;折射波法在地震勘探中的应用。
§2.5 地震波动力学理论及应用本节不作为考试内容。
第三章地震资料采集方法与技术§3.1 野外工作概述1、掌握基本概念:低(降)速带、频散、群速度、相速度、多次波、虚反射、鸣震、交混回响。
2、掌握基本内容:试验工作内容、生产工作过程、激发条件、接收条件、调查干扰波的方法、干扰波的类型、各种干扰波的主要特点、面波特点、压制面波的方法、海上地震勘探的特点与特殊性、海上特殊干扰波、海上震源等。
地震成像引言在地面记录的未经偏移的反射地震能量(例如共偏移距记录)通常可提供地下构造的强相干图像。
尽管是强相干的,但该图象是不正确的,原因在于它受到几种畸变作用的影响,最显著的影响是来自地质界面截断的绕射波和倾斜界面上反射点与地面位置间能量的横向移动。
地震成像就是校正这些畸变影响的地震处理手段。
在常规处理中,地震成像是众多处理手段的最终阶段。
全部地震处理过程包括能量增益化、反褶积、静校正和速度分析、NMO和DMO、叠加,最后到主要用于构造成像的偏移。
根据偏移数据体绘制地震构造图,然后在构造图上布置钻探目的层的勘探井位。
这种做法一直延续至今,但是偏移后的数据已逐渐被用于后续进一步处理的输入数据。
例如,当地质学家和工程师们联合攻关已做过3D地震工作的油田的油气开发问题时,他们把偏移后的地震数据体用作地震属性分析处理的输入数据,诸如瞬时相位或振幅随偏移距变化的处理中。
通过已知井与建议井位处地震属性特征的对比,我们就能预测出建议井位处的岩性特征。
例如,倘若在几口井中某一属性可能都与孔隙度观测值的相关性很好,那么这个属性就可用于估算其它井位处的孔隙度。
多年来,地球物理学家都在用地震数据在比地震波长小得多的尺度上估算岩石物性特征,但直到最近他们才对偏移数据实施这种处理,它代表了偏移应用的一个全新领域。
这也促使我们在偏移中要非常重视象振幅这样的参数的应用。
时间和深度、叠后和叠前、2D和3D地球物理学家习惯使用时间域记录的地震数据。
以往,他们乐于解释时间记录,甚至是偏移后的时间记录。
这就与在深度域工作的地质学家和油藏工程师产生了根本性的矛盾,同时也给许多通过设计偏移算法来摸索地球物理学规律的数学家和物理学家带来困惑。
尽管深度域的解释结果明显会更正确,但迫使地球物理学家这样做会妨碍他们解释地震数据的能力,例如在资料质量差的地区,在偏移前后剖面上难以确定一一对应的反射波同相轴。
所以,就有了时间域偏移(时间偏移)和深度域偏移(深度偏移),它们之间的差异远不是表面上的。
地震成像方法
地震成像是一种重要的地震技术,它利用通过地壳中反射回来的弹性波形,来分析地壳中不同层次和结构的空间分布,以了解地质构造。
它可以将深层的隐藏地质体转换为与地表可见的形态相近的数据形态,这样它们可以更准确的被采集到和操作。
地震成像方法一般分为两类:逆时偏移(Reverse Time Migration,RTM)和差分时偏移(Difference Time Migration,DTM)。
RTM是一种基于弹性波传播的对称的正演方法,它可以根据多种类型的地形条件探测到积聚成图像的波形变化记录,以及地表表征信息。
它最大的优势是允许目标层次和尺度重构,使更多的反射特征能够被捕捉到地面或者测量位置。
DTM是一种采用差分法提取能量信息的非对称正演方法,它可以提取更小尺度下目标成像断层物理特征,这样就可以准确定位断层位置,及其关联的定量地震参数包括衰减,场设计及激发效果。
DTM和RTM方法的结合利用双水平(主和次)反射抛射,可以用于全方位的成层结构反演,以及分段连续的深部断层特征信息的反演和定位。
因此,这种技术可以保证地震成像在高普适性和精确表达的能力方面的效率与可靠性。
总之,地震成像方法是一种准确、可靠、通用的地震技术,它不仅可以检测出深层隐藏的地质形态,而且还可以进行全面的成层结构反演和断层特征的定位。
第一章地震偏移成像基础地震偏移技术是现代地震勘探数据处理的三大基本技术之一。
它是在过去的古典技术上发展起来的,其它两大技术都是从其它相关学科引进到地震中来的。
所以,偏移技术具有地震勘探本身的特征。
但是,地震偏移方法本身由于使用计算机而引起了许多革命性的变化。
这就是把它从研究简单的探测目标的几何图形进而发展成研究反射界面空间的波场特征、振幅变化和反射率等。
本章主要介绍地震偏移成像技术的基础知识。
首先给出偏移成像的概念;第二节介绍有限差分法的基础知识;第三节叙述基于波动方程的波场外推与地震成像原理;第四节讨论波场外推的Kirchhoff积分法;第五节简单分析Born近似和Rytov近似;最后阐述基于De Wolf近似、薄板近似、屏近似和相屏传播算子计算反向散射波场的方法。
§1.1 偏移成像的概念反射地震方法是根据在地面上以一定方式进行弹性波激发,并在地面的一定范围(孔径)内记录来自地下弹性分界面的反射波来研究地下地质岩层结构及其物性特征的一种方法。
因此,也可以把它看做是一种反散射问题。
就反射地震观测方式的特点,它的成像问题要分做两步,第一步是按照一定的方式记录到达地面的反射波,第二步用计算机按一定的计算方法对观测数据进行处理,使之成为反映地下地质分层面位置及反射系数值的反射界面的像。
而地震偏移技术就是在第二步过程使反射界面最佳地成像的一种技术。
地震偏移可在叠前做也可在叠后做。
叠前偏移是把共炮点道集记录或共偏移距道集记录中的反射波归位到产生它们的反射界面上并使绕射波收敛到产生它的绕射点上。
在把反射波回投到反射界面上和绕射波收敛到绕射点上时要去掉传播过程的效应,如扩散与衰减等。
最后得到能够反映界面反射系数特点的并正确归位了的地震波形剖面,即偏移剖面。
叠后偏移是在水平叠加剖面的基础上进行的,针对水平叠加剖面上存在的倾斜反射层不能正确地归位和绕射波不能完全收敛的问题,采用了爆炸反射面的概念来实现倾斜反射层的正确归位和绕射波的完全收敛。
第一章地震偏移成像基础地震偏移技术是现代地震勘探数据处理的三大基本技术之一。
它是在过去的古典技术上发展起来的,其它两大技术都是从其它相关学科引进到地震中来的。
所以,偏移技术具有地震勘探本身的特征。
但是,地震偏移方法本身由于使用计算机而引起了许多革命性的变化。
这就是把它从研究简单的探测目标的几何图形进而发展成研究反射界面空间的波场特征、振幅变化和反射率等。
本章主要介绍地震偏移成像技术的基础知识。
首先给出偏移成像的概念;第二节介绍有限差分法的基础知识;第三节叙述基于波动方程的波场外推与地震成像原理;第四节讨论波场外推的Kirchhoff积分法;第五节简单分析Born近似和Rytov近似;最后阐述基于De Wolf近似、薄板近似、屏近似和相屏传播算子计算反向散射波场的方法。
§1.1 偏移成像的概念反射地震方法是根据在地面上以一定方式进行弹性波激发,并在地面的一定范围(孔径)内记录来自地下弹性分界面的反射波来研究地下地质岩层结构及其物性特征的一种方法。
因此,也可以把它看做是一种反散射问题。
就反射地震观测方式的特点,它的成像问题要分做两步,第一步是按照一定的方式记录到达地面的反射波,第二步用计算机按一定的计算方法对观测数据进行处理,使之成为反映地下地质分层面位置及反射系数值的反射界面的像。
而地震偏移技术就是在第二步过程使反射界面最佳地成像的一种技术。
地震偏移可在叠前做也可在叠后做。
叠前偏移是把共炮点道集记录或共偏移距道集记录中的反射波归位到产生它们的反射界面上并使绕射波收敛到产生它的绕射点上。
在把反射波回投到反射界面上和绕射波收敛到绕射点上时要去掉传播过程的效应,如扩散与衰减等。
最后得到能够反映界面反射系数特点的并正确归位了的地震波形剖面,即偏移剖面。
叠后偏移是在水平叠加剖面的基础上进行的,针对水平叠加剖面上存在的倾斜反射层不能正确地归位和绕射波不能完全收敛的问题,采用了爆炸反射面的概念来实现倾斜反射层的正确归位和绕射波的完全收敛。
第一章地震偏移成像基础地震偏移技术是现代地震勘探数据处理的三大基本技术之一。
它是在过去的古典技术上发展起来的,其它两大技术都是从其它相关学科引进到地震中来的。
所以,偏移技术具有地震勘探本身的特征。
但是,地震偏移方法本身由于使用计算机而引起了许多革命性的变化。
这就是把它从研究简单的探测目标的几何图形进而发展成研究反射界面空间的波场特征、振幅变化和反射率等。
本章主要介绍地震偏移成像技术的基础知识。
首先给出偏移成像的概念;第二节介绍有限差分法的基础知识;第三节叙述基于波动方程的波场外推与地震成像原理;第四节讨论波场外推的Kirchhoff积分法;第五节简单分析Born近似和Rytov近似;最后阐述基于De Wolf近似、薄板近似、屏近似和相屏传播算子计算反向散射波场的方法。
§1.1 偏移成像的概念反射地震方法是根据在地面上以一定方式进行弹性波激发,并在地面的一定范围(孔径)内记录来自地下弹性分界面的反射波来研究地下地质岩层结构及其物性特征的一种方法。
因此,也可以把它看做是一种反散射问题。
就反射地震观测方式的特点,它的成像问题要分做两步,第一步是按照一定的方式记录到达地面的反射波,第二步用计算机按一定的计算方法对观测数据进行处理,使之成为反映地下地质分层面位置及反射系数值的反射界面的像。
而地震偏移技术就是在第二步过程使反射界面最佳地成像的一种技术。
地震偏移可在叠前做也可在叠后做。
叠前偏移是把共炮点道集记录或共偏移距道集记录中的反射波归位到产生它们的反射界面上并使绕射波收敛到产生它的绕射点上。
在把反射波回投到反射界面上和绕射波收敛到绕射点上时要去掉传播过程的效应,如扩散与衰减等。
最后得到能够反映界面反射系数特点的并正确归位了的地震波形剖面,即偏移剖面。
叠后偏移是在水平叠加剖面的基础上进行的,针对水平叠加剖面上存在的倾斜反射层不能正确地归位和绕射波不能完全收敛的问题,采用了爆炸反射面的概念来实现倾斜反射层的正确归位和绕射波的完全收敛。
地震偏移的效果见图1-1和图1-2。
地震偏移的类型见表1-1。
地震偏移技术在二十世纪六十年代以前是用手工操作的一种制图技术,只是用来求得反射点的空间位置,而不考虑反射波的特点。
它是一种古典的偏移方法。
早期的计算机偏移方法是在古典的偏移方法的基础上提出来的。
其中有的成功了,有的失败了。
成功的是那些符合波的传播特征的方法。
尽管这些方法使用了波前、绕射等地震波传播的惠更斯原理,但只是定性的、概念性的。
偏移剖面的质量虽然能够满足最基本的要求,但归位的精度和成像时的波形特征都不是很准确的。
因此,研究更有效的地震偏移方法是很迫切的。
二十世纪七十年代初J.Claerbout教授首先提出了用有限差分法解单程波动方程的近似式,用地面观测的地震数据重建地震波在地下传播过程中的波场,从这些传播过程的波场中提取使地震界面成像的那些数据,组成地震偏移剖面。
由于这种偏移方法在计算过程中要解波动方程或其近似式,所以被称为波动方程法偏移技术。
以后,French和Schneider等在绕射偏移法的基础上使用了波动方程解的Kirchhoff积分公式,发展为地震偏移的波动方程积分法。
使绕射偏移建立在可靠的波的基本原理上。
因而改善了偏移剖面,取得了良好的效果。
图1-1 (a)共中心点叠加剖面,(b)偏移剖面,(c)明显的绕射波D、偏移前的倾斜同相轴B和偏移后的倾斜同相轴A。
偏移使倾斜同相轴B归位到它的真实地下界面A,并使绕射波D收敛到其顶点P。
点画线指出了盐丘的边界。
图1-2 偏移前(a)及偏移后(b)弯曲反射界面(向斜和背斜)的形状。
详细情况见正文(模型据Union Oil Company)在二十世纪七十年代后期,Stolt和Gazdag等又先后提出了在频率-波数域解波动方程,外推地震波场的方法。
这种方法被称为F-K域偏移方法。
由于该方法计算简单,效率高,因而很快得到了推广。
上述三种波动方程偏移方法是同时并存的,因为它们各有自己的特点,因而不能用一个方法来取代其它方法。
使用时视具体条件和要求决定采用何种方法。
波动方程偏移方法在最近20年间迅速发展并不断完善,许多人对此做出了有益的贡献。
其中,Loewenthal等人的爆炸反射面的概念对于理解叠加剖面的偏移成像具有很大价值,Hubral,Larner等人提出的深度偏移的概念具有很大意义,Berkhout提出的偏移过程是一个空间褶积的概念对于偏移的横向分辨力的理解很有益处,马在田院士提出的高阶方程的分裂算法对提高有限差分法偏移的精度有很大贡献,Yilmaz等提出的双平方根法为解决叠前偏移奠定了基础。
现在仍有许多学者还在探索波动方程偏移技术,以期更加完善该方法。
表1-1 偏移方法分类§1.2 有限差分法的基础知识在计算机上进行数值运算,使用的是离散的和有限的数值,而不是连续的和无限的函数,为此要为离散数值的计算建立基本方法。
最基本和广泛使用的方法就是有限差分法,借助这一方法可以研究连续物理系统的性质,近似地、但相当精确地解出各种数理方程问题。
本节将概要地介绍有限差分法的基础知识,供以后各章解偏移方程使用。
一.差分方程的建立与求解1.有限差分的概念我们感兴趣的是用有限差分法解各类微分方程,因此,要把导数用有限差分来近似,所以我们这里只研究用有限差分近似导数的方法。
(1)一元函数差分法当一个函数u 和它的各阶导数是变量x 的单值的、有限的和连续的函数时,可以用泰勒原理展开为:⋅⋅⋅+∆+∆+∆+=∆+3332226121)()(dx u d x dx u d x dx du x x u x x u (1-1) 和⋅⋅⋅+∆-∆+∆-=∆-3332226121)()(dxu d x dx u d x dx du x x u x x u (1-2) 以上二式相加,得:)()(2)()(4222x dx u d x x u x x u x x u ∆O +∆+=∆-+∆+(1-3) 式中O(Δx 4)表示包含有Δx 的四阶和高于四阶以上的项。
从(1-1)~(1-3)式我们可以导出(图1-3)下列表示式:一阶向前差分:)]()([1x u x x u xu x -∆+∆=+δ )(x dxdu ∆O +=(1-4) 一阶向后差分: )]()([1x x u x u xu x ∆--∆=-δ )(x dx du ∆O +=(1-5) 一阶中心差分:)]()([21x x u x x u xu x ∆--∆+∆=δ)(2x dx du ∆O += (1-6)二阶差分:)]()(2)([122x x u x u x x u x u u x xx ∆-+-∆+∆==δδ)(222x dx ud ∆O +=(1-7) (2)多元函数差分法图1-3 一元函数差分法 图1-4 二元函数差分网格设有一个二元函数u(x, t),我们用网格把它们离散(图1-4),令x=i Δx ,t=j Δti 和j 为整数。
用u i, j =u(i Δx, j Δt)表示各网格点上的函数值。
现在我们用泰勒级数展开下面各点之值为:),(],)1[(,1j i j i t x x u t j x i u u ∆+=∆∆+=+⋅⋅⋅+∂∂∆+∂∂∆+∂∂∆+=j i j i j i j i x ux x u x x ux u ,333,222,,)(61)(21)((1-8)),(],)1[(,1j i j i t x x u t j x i u u ∆-=∆∆-=-⋅⋅⋅+∂∂∆-∂∂∆+∂∂∆-=j i j i j i j i x ux x u x x u x u ,333,222,,)(61)(21)((1-9)由此可以求出对x 的一阶和二阶差分为:)(1,,1j i j i x u u xu -∆=++δ )(x x u ∆O +∂∂=(1-10) )(1,1,j i j i x u u xu ---∆=δ )(x x u ∆O +∂∂=(1-11) )(21,1,1j i j i x u u xu -+-∆=δ )(2x xu ∆O +∂∂=(1-12) )2(21.1,,12j i j i j i x xx u u u xu u -++-∆==δδ )(222x xu ∆O +∂∂= (1-13) 同理可求出对t 的一阶和二阶差分:)(1,1,j i j i t u u tu -∆=++δ )(t tu ∆O +∂∂=(1-14) )(11,,---∆=j i j i t u u tu δ )(t tu ∆O +∂∂=(1-15) )(211,1,-+-∆=j i j i t u u tu δ )(2t tu ∆O +∂∂=(1-16) )2(11,,1,22-++-∆==j i j i j i t tt u u u tu u δδ )(222t tu ∆O +∂∂=(1-17) 2.建立差分方程建立差分方程的方法很多,有积分法,物理量守恒法,变分法和最小平方法等。
对于常系数的微分方程来说,积分法是最简便的。
用积分法构造差分方程的过程见参考文献[1]。
3.差分方程的格式差分方程的格式基本可分为两大类:即显式格式和隐式格式。
在实际工作中又可以衍生出许多格式,甚至可以用显式与隐式联合形式的差分格式。
在这里我们仅以抛物型偏微分方程为例说明常用的显式格式和Crank-Nicolson 格式。
(1)显式格式设有抛物型偏微分方程:022=∂∂-∂∂x u t u σ (1-18) 经推导可具体写出差分方程为:)2(,1,,12,1,j i j i j i j i j i u u u x tu u -+++-∆∆+=σ (1-19)所用的差分网格如图1-5所示。
令γσ=∆∆2x t,则(1-19)式可写为:)2(,1,,1,1,j i j i j i j i j i u u u u u -+++-+=γ (1-20)其中j i u ,,j i u ,1-和j i u ,1+为已知值。
从第j 时间层上的已知值,可用(1-20)式直接计算出第1+j 时间层的值。
以此类推,可解出全时间上的物理量值。
图1-5 显式差分格式 图1-6 隐式差分格式(2)Crank-Nicolson 隐式格式虽然显式法计算上很简单,但它有一个严重的缺点,即时间步长t ∆一定要很小,必须满足2102≤∆∆<x tσ,才能保持计算上稳定并达到必要的精度。
Crank-Nicolson 在1947年提出了一个使所有的有限r 值都满足计算要求(收敛性和稳定性)的方法。
他们把22/xu ∂∂项用第j 和第j+1时间行上的平均差分来逼近。
求出下列的差分方程(图1-6): ⎭⎬⎫⎩⎨⎧∆+-+∆+-=∆--++-++++2,1,,121,11,1,1,1,222x u u u x u u u t u u j i j i j i j i j i j i ji j i σ (1-21) 由此得出下列等式:j i j i j i j i j i j i u u u u u u ,1,,11,11,1,1)1(2)1(2+-++++-+-+=-++-γγγγγγ (1-22) 其中,2x t∆∆=σγ。