Tesseral 2D正演模拟培训讲义 下载本文

FOLLOW ME(正演模拟培训讲义)

从一个bmp文件建立一个模型并进行全波场模拟: (一)粘贴一个bmp文件

1、给一个深度比例尺的速度—深度模型图在画图软件中把模型的比例尺和

范围做准,并输出一个bmp格式文件备用。 2、在工具条上点击 打开Tesseral主界面(图1),并填写正确的顶、底、

左、右坐标,并在Surface 组参数中选“Invisible”不可见地表 ,如果是地震(声波)波场模拟,地表一定不能出现。“Invisible”表示上行波将不再向下反射回到该模型按“确定”后,在删掉跳出的物理参数填写表后,出现一个空白的坐标网格图(图2)。

图1

图2

1

3、在工具条上选择 点左键光标变为多边形

,表示可以开始画图。

手工在屏幕上画多边形请看”Tutorial-Hinds&Kuzmiski May 2002”。现在讲如何从一个bmp文件描绘一个模型。从坐标外的一点开始,画一个矩形后双击,出现一个红色边界的矩形框后出现参数表,可以不填,点击“确定”,则完成了一个空白的矩形的建立(图3)。

4.点击

2

图3

打开文件,再点击Picture Files,选出准备好的bmp文件(图4)。

图4

5.文件显示的窗口中,确认所选的顶、底、左、右坐标无误后,点击OK(图5)。

图5

6.bmp文件就以正确的比例贴在Tesseral 的绘多边形的窗口上(图6)。

图6

3

(二)描绘多边形

1.开始描绘第一个多边形,从左侧的坐标线外起始,然后精细地用鼠标描绘地形线到坐标的右侧外点一下,然后下拉到底部,点一下,再向左下角拉线,点一下,然后双击左键矩形封闭。在出现的参数表上填写第一层的参数2800米,其余用缺省值(图7)。 图7

2.完成后点击OK,第一个多边形则做好了(图8)。 图8

4

3.为了看清楚原图的线,也可以选用 消色键(图9)。

图9

4.画以下的多边形,原则是整的多边形从坐标线外起始画,内部多边形在某一个多边形内部封闭,当两个多边形重叠以最上面呈现的部分有效,压在下面的失效(每个多边形的参数表,你都填写正确的P波参数,其余用缺省值,原则是从上到下,从整到零。并且在放大比例尺的条件下精细描绘(图10)。 图10

5

5. 当你将bmp文件上的所有层位都按一个个多边形描绘下来,就可以用 着色键将模型着上颜色看看描绘得是否与bmp图一致,如无误就可以撤掉bmp文件底图。方法如下:点击View,在出现的下拉菜单上选Pad Image,在它的侧出菜单上选Unload,点击后,这张底图就会消失。(图11)

6.编辑多边形,将光标指向要编辑的多边形点左键,击活了多边形的边框。移动光标,当坐标呈十字交叉的箭头状,即可用加点或拖拉已有点的办法,编辑多边形的各个边线(图12)。

6

图11

图12

7.需要加点则点右键出现下拉菜单选Insert Note,需要删除点将光标指向已有的点,点右键出现下拉菜单,选Delete Node即可删除。当看不清,想放大时,则点工具条上的 键,当 鼠标变成 + 按左键画框,可以放大细看,可在放大的窗内在鼠标呈+ 状下,再次画框放大,精细调整每一个多边形(图13)。

图13

8.输入岩性符号

在要输入岩性符号的多边形上击右键,出现的下拉菜单中选Edit Polygon点击左键,出现参数表,在左下方点Loal点击后出现许多岩性符号,点中要选的岩性符号则出现在Load上方的白板上,然后点OK则在该多边形中加入了这种岩性符号,用此方法可逐一加入每个多边形的岩性符号(图14)。

7

图14

(三)调整颜色

1.点击工具条上的色标调整钮,则色标出现在屏幕的右侧(图15)。 图15

2.点击工具条上的彩色按钮,可以在色标上直接编辑颜色(图16)。 图16

3.点击要编辑的速度区间小节点,即可出现一个对话框(图17)。

图17

8

当选中你要替换的颜色,点击后它就进入Select色块,点Additonal 模型上该区段上的颜色就会改变(图18)

4.调节满意后,可以存储该模型更换上一个模型。

图18

(四)定义观测系统

1.定义震源参数,在工具条上点击 Framework按钮,出现对话框(图19)。

图19

9

点击Source 按钮,出现震源定义会话框,首先定义放炮方式,该表中共有三方式可选Point 是最常用的模式,即对应于野外生产中的逐点放炮方式(图20)。

图20

1)如果我们的生产方式是等间距规则地移动炮点,选Free 即可将炮点安放到任意所需要的位置上。 2)Free 右侧的选项中,地面放炮选择Horizontal Line (井中放炮选择Vertical line)

3)需选中Cable interval则表示用户可以在测线上定义等间隔移动的炮点。 4)定义炮数则在Number 右侧的空格内输入要计算的总炮数。 5)定义炮间距在Interval 右侧的空格内输入参数。

6)定义要计算的炮数,缺省值从第一炮计算到第观测系统定义的最后一炮,也可以人工输入要计算的第一炮和最后一炮。

7)定义振源方式可选Compression(纵波)或Rotational (横波)。

8)定义子波类型,可选Single (单峰)Symmetric(对称单峰)Double(双峰)。

9)“Frequency”振源峰值频率其代表其将是生成子波最大振幅的频率。

2.定义检波点排列参数

1)在工具条上点击Framework按钮,出现对话框。再点击Observation按钮,

出现检波点会话窗口(图21)。

图21

10

2)同样在Free 的右侧选择中选Horizontal {表明检波点沿地面安放}Vertical 表明检波点在井下垂直安放。 3)“Free” 选中表明检波点可从任意位置开始。 4)“Cable Interval” 选中表明检波点按等投影间距安放在起伏的地面或斜井中。

5)Move with Sourc表明检波点随炮点一起整体移动。

6)定义排列位置From第一个检波点的位置,to到排列最后一个检波点的位置。

7)定义炮间距,在“Interval”输入。 8)“Margin”建议采用缺省值{用于定义实际的网格尺寸}。 9)“Start”定义起始时间,缺省值自动切除近炮点最大振幅值的切除时间。 10)“Stop”定义用户要计算的最大波传播时间。 11)“Sample”定义要计算合成记录的采样率,缺省值为一般确定的采样值。

3.定义计算波场快照的数量

1). EVERY:定义每隔几个炮点生成一个波场快照。 2). 选”End Truncation”则在计算结尾截断以适应我们原先定义的接收区和截止时间。

3).“Start”定义波场快照的起使时间。

4) ” Sample”定义隔多少时间拍一张波场快照.

4.将检波点和炮点移置地面

1). 将鼠标放在第一个检波点上,当鼠标变为检波点标志符时,将第一个检

波点拖至它的正确坐标位置。

2) 将鼠标放在第一炮点上,当鼠标变为炮点标志符时将炮点放在其正确位

置,也可用精调按钮,输入坐标的方式定位。 3). 我们现在希望用FINE TUNE ICON精调震源和检波点的

位置。在工具条上点击红圈标出的键,然后把指针放在你想要移动的目标体上, 目标体顶部就回出现一个模糊的

指针影象, 按左键点亮该目标体,并把它拖放到你想精调的位置上,在TUNE POSITION 窗口对该目标体键入新的坐标然后点击ENTER。 一的虚象目标体便飞到它的精调位置。关掉TUNE POSITION 窗口 则精调结束。

在这个例子中, 我们希望动整个震源的排列。将指针移动到第一个震源位置,第一个震源在X= -400 和 Z= -150 M, 如(图22)所示 随着震源被光标点亮并被拖放到新位置以后, 在精调菜单中键入X 为200 和Z 为 –150…

11 22 图

然后,在TUNE POSITION菜单中按的“INTER”,震源即移到正确的位置上(图23)。

图23

现在最左边震源的坐标 X=200 M 和 Z=-150 M。用点击最上角的X号 来关闭精调窗口。

(五)计算

当你建好一个模型以后可能用多种方法和方式试验、计算。为了将这些计算结果安全地保留,最好将模型对应不同的计算分别命名。(为了避免与程序加版本数字序号相区别,建议加字母序号加以区别,进行有效的文件管理。)

1. 第一个作业我们试验声波方程的单炮模型

(为此在模型的名称后面加一个大写的字母A,并复制一个模型备份。)

1) 点击该模型,模型即被调到主窗口。也可以用另一途径将模型调入主窗口,即打开Tesseral 2D窗口,将所出现的新建模型参数的子窗口一一取消关闭后,在空白主窗口的条件下,用打开文件---在模型文件中点击要打开的文件---然后打开。要调用的模型就出现在屏幕上了。 2) 在工具条上点击Framework窗口按钮 ,打开Framework窗口,点开source参数窗口,将计算组参数中的last和first的数字填写为同样的数,则表示只计算这一炮的炮记录和波场快照,此例中表示只计算81炮。 (图24)

图24

12

然后再打开observation 参数菜单检查修改炮记录的时间长度和采样率,同时检查并修改波场快照参数,注意为了提高计算速度,算波场快照的炮点数不要太多,波场快照的采样间隔可以选大,如现在我们选定的参数是每50ms照一张,每张8秒长的记录从20毫秒起照,每炮则照160张像。(图25)

图25

3)RUN,出现下拉菜单,选择声波方程计算。(图26)

图26

当计算引擎运行时计算对话框显示在屏幕上,它使用户能够看见当前处理状态。 (图27)

13

图27

4).当按下按钮 “Pause” 则停止该处理, 当处理暂停时,该按钮的标题变为“Resume”. 再次按下它即恢复处理。 5).用户按下窗口右上角的 “Hide” 能使对话窗最小化,如果计算对话框被最小化,按下在主窗口右下角的处理图标就能使它出现。按下“Terminate” 能终止处理。在终止处理时出现警告对话框:(图28)

图28

6).在计算时浏览器则启动监视该处理使你能够通过显示的图形,看到当前计算的进度。监视时,你可以看到随着计算而增长的炮集记录和波场快照(图29)。可以通过改变窗口的大小和点击 调出均衡、增益、平滑处理窗口来控制显示效果。(图30) 图29

14

图30

7).画图监视需要大量的时间,影响计算速度。如果你不需要图形监视,可采

用最小化浏览器窗口,这样处理器不进行画图处理,计算速度将加快。

8).当计算完成时,浏览器窗口显示最后计算的一个炮集和波场快照结果。

9).当你开始计算时,又想调整参数,可以用Pause键人为停止程序的执行,也可以终止程序,而且也可以从程序终止处或在程序崩溃后恢复计算。 10).如果计算没有结束时,你又调用建模器,你将见到下面的警告信息:

(图31)

不同的回答产生的结果不同:

a) Yes: 继续进行未完成的计算。

b) No: 不继续进行未完成的计算,也不保留已得到的计算的结果。

c) Cancel: 用户从建模器退出,计算结果没有改变。.

在程序计算运行期间和用户终止计算时,程序自动写入恢复数据。

15

11)报告窗口.(图32)

报告窗中显示关于计算和数据操作的辅助信息:

图32

在计算开始时报告窗中自动打开,也可以用“View” 的下拉菜单中的选项“Show Report”手动打开。 12) . 发送到报告窗口的数据,最后被存储在当前目录下,名为 “Tesseral.log”的文件中,供以后分析研究使用。 13)计算结束,计算完的第81炮声波模拟的炮集记录和第81炮的波场快照自动显示在浏览器窗口。(图33) 图33

14)该作业输出的文件都自动存在你原先存储模型文件的文件夹中。你可以从工具条上用打开 键查阅。输出的文件有:模型名+GatherAP-81.tgr, 为模拟的炮集记录, 炮集用Gather代表, AP表示用声波模拟算法,81是炮号, tgr 是表示文件类型的后缀; 模型名+SnapAP-81.tgr, 为模拟的炮集波场传播快照, 快照用Snap代表, AP表示用声波模拟算法,81是炮号, tgr 是表示文件类型的后缀;模型名+TimeAP-81.tgr, 为模拟的炮集波前传播的时间能量场, 时间能量场用Time代表;AP表示用声波模拟算法,81是炮号, tgr 是表示文件类型的后缀。

16

炮集记录、波场快照示于图34;时间能量场示于图35。 图34

图35

2.第二个作业我们试验弹性波方程的单炮模型(仍然计算同一个模型的第81炮) 1) 打开模型的方法与第一个作业相同,但是当你仍选这同一个模型打开文件时,屏幕上只有你上次浏览的炮集和快照,模型图被隐藏在后面的窗口中,刚才计算完毕后,模型在上一排的第一个窗口上,炮集和波场快照在下一排的第一和第二个窗口中显示。已经在窗口中的图,不可能在已有图形的窗口中重调。所以你要再次使

17

用它,只需将表面的窗口关掉,而把模型窗口呈现出来。

2)模型参数未变,只需重复上述步奏点击RUN,在出现的下拉菜单中,选择弹性波方程计算(Elastic Modeling)

3)其他的做法都与上面相同。只是最后输出的结果是弹性波方程模拟的结果呈现在浏览器窗口。炮集和波场快照(图36),时间能量场(图37)。

图36

图37

18

4)该作业输出的文件都自动存在你原先存储模型文件的文件夹中。你可以从工具条上用打开 键查阅。输出的文件有:模型名+GatherEP-81.tgr, 为模拟的炮集记录, 炮集用Gather代表, EP表示用弹性波模拟算法,81是炮号, tgr 是表示文件类型的后缀; 模型名+SnapEP-81.tgr, 为模拟的炮集波场传播快照, 快照用Snap代表, EP表示用声波模拟算法,81是炮号, tgr 是表示文件类型的后缀;模型名+TimeEP-81.tgr, 为模拟的炮集波前传播的时间能量场, 时间能量场用Time代表;EP表示用声波模拟算法,81是炮号, tgr 是表示文件类型的后缀。

3. 第三个作业我们试验弹性波方程模拟产生50炮中点放炮记录。1)检查你要做的模型文件,将其更名放入一个新建的文件夹内。同时检查新建文件夹所在的磁盘空间是否充足。(必须准备5-10G)的磁盘空间才能准备这一个大作业。 2)打开Tesseral 2D窗口,将所出现的新建模型参数的子窗口一一取消关闭后,在空白主窗口的条件下,用打开文件---在模型文件中点击要打开的文件---然后打开。要调用的模型就出现在屏幕上了。

3)在工具条上点击Framework窗口 按钮 ,打开Framework窗口,点开source参数窗口,将计算组参数中的last和first的数字填写为50和1的数,则表示在这个已定义的观测系统参数下只计算第一炮到第50炮的记录。 (图38)

图38

4)点开observation参数窗口,该参数表中表示每个排列240道检波器,道间距30米,排列总长度7170米,每炮记录从0.02秒开始记录,记录长度为7秒,采样间隔为4毫秒;每隔10炮照一个波场快照,每个波场快照,每隔50毫秒照一张像。(为了提高计算速度可以不要照波场快照,即将Every填写为0即可)(图39)。

19

5)参数全部检查核实完后在主菜单上选Run在其下拉菜单上选Elastic Modeling点击,程序开始处理,这一处理时间需要一整天。计算完毕,在浏览器窗口呈现出模型、计算的最后一个炮集和最后一个波场快照。

6)该作业的输出文件都自动存在你原先存储模型文件的文件夹中。你可以从工具条上用打开 键再按Viewer fieles查阅。输出的文件有:模型名+GatherEP-炮号.tgr, 为模拟的炮集记录; 模型名+SnapEP-炮号.tgr, 为模拟的炮集波场传播快照;模型名+TimeEP-炮号.tgr, 为模拟的炮集波前传播的时间能量场。其他的文件还有模型名+WaveEP-1.tgr 其中记录的是模拟用的子波。Tesseral文本文件,其中记录了该作业的使用参数,和每一炮记录的运行过程。模型名-0.avo是这一个作业产生的文件但不列在可打开的文件之中。

4. 第四个作业我们试验使用弹性波方程模拟产生50炮中点放炮记录,做切除. 1) 合并单炮记录为一个整的炮集记录文件供后续处理使用。点击主菜单上的Run,然后在它的下拉菜单上选Grid Merge,程序自动将单个的炮集记录合并为一整体炮集录文件,存储在同一个文件夹内,文件名为模型名+GatherEP.tgr。同时程序出现提示窗口,询问是否删除原有的炮集记录?如果点Yes则删除原有炮集,如果点No则将原单个炮集文件都保留。(如果计算机的硬盘不够大,建议采用Yes)。

1)打开一个空白窗格,使之活化,然后打开“模型名+GatherEP.tgr”合并的炮集文件。

2)做初至切除和底部切除。点击主菜单上的Edit,在它的下拉菜单上选Muting,在Muting的侧向菜单上点Upper Line, Upper Line被激活状态为其左侧的小方格内打一个小对勾。此时可以利用鼠标左键在炮集录上点切除线,当点完一个炮集以后,每一个炮集上都画有相同的切除线。如果这条切除线适合于所有的炮你可以不用逐炮一一手工定义。但是如果每炮不同于这个个参考线。

20

39 图

则需要逐炮定义。当你想再次编辑切除线,则必须再次点击主菜单上的Edit,在它的下拉菜单上选Muting,在Muting的侧向菜单上点Edit Upper Line,当出现编辑菜单时,点OK此时编辑功能被激活,切除线变为黄色,即可进行编辑。

当所有炮都编辑完以后,再次点击主菜单上的Edit,在它的下拉菜单上选Muting,在Muting的侧向菜单上点Upper Line使之消除活化状态, 方格及小对勾消失,即完成了所有炮集的初至切除(图40)。

图40

与上面情况类似可以进行底部切除,即把上面步奏中的Upper Line,变为Lower Line和Edit Lower Line。另外如果没一炮的底切除线如果都一样也不用逐炮绘,直接再次点击主菜单上的Edit,在它的下拉菜单上选Muting,在Muting的侧向菜单上点LowerLine使之消除活化状态, 方格及小对勾消失,即完成了所有炮集的底切除(图41)。 图41

21

5. 第五个作业我们试验在切除后的炮记录上直接做速度谱分析 1)点击主菜单上的Edit,在它的下拉菜单上选Velocity,在Velocity的侧向菜单上点Create Velocity Spectra,然后出现速度谱参数菜单(图42)。

图42

在该菜单中输入如下适应速度分析的参数: ? Start Velocity –适宜的最小速度 (m/s); ? End Velocity -适宜的最大速度(m/s); ? Step Velocity –合适的步长 (m/s);

? Gradient – 时延斜率– 通过输入一个数值 或用移动滑棒确定。该斜率定义信号在最大排列处的动校正量(NMO );在值较大的情况下信号保留较多;在值较小的情况下,信号超过给定的界限,就会在速度谱分析中漏掉。

? number 在本次分析中所包含的炮数。

o “Every” 选项 –适用于对该剖面中的所有炮进行分析,这将需

要一会儿时间; o “First&Last” 选项 – 适用于仅对第一和最后的一些道进行分

析; o “Number” 选项 – 适用于对给定一定数量的道进行分析,从当

前所在的道开始算起。

沿着剖面的炮点坐标也被作为信息显示。

点击 OK 键 开始计算速度谱(图43)。

22

图43

计算完毕速度谱被显示在相邻的窗口(图44)。为改变速度谱的显示,从主菜单条选择 “View” → “Visualization”,或点击 按钮然后选择 “Basic”或 “Advanced” (详见 Tesseral 2-D 用户手册6.6节). 建议用10-20ms的时窗对速度谱进行平滑处理。 (“Smoothing选项”)。

图44

23

6. 第六个作业我们试验由炮记录抽取CMP道集

1)输入记录是合并后的炮集录文件“模型名+GatherEP.tgr”。

2)在主菜单上选Run,在它的下拉菜单上选Gathering(GATHER)。见(图45) 点击后屏幕上即出现一个有关道集参数的菜单,它是根据你定义的观测系统而得出的。(图46) 图45

图46

24

上面标明的输出文件名及存储的地点。右侧的三个参数Start X对应于首CDP号所应该对应的X坐标;Stop X对应于末CDP号所应该对应的X坐标;X Step为 CDP间距。

3)检查首尾CDP号与CDP间距都无误后,点击小窗口上的Advanced选项,选取道集分选类型参数。这里我们选择共中心点道集。(图47) 图47

4)参数选完,点击OK后程序开始计算,屏幕上有作业动态窗口供我们监视处理,道集分选处理的速度是比较快的,但此后还要进行网格化处理,这需要较长的时间。(图48) 图48

25

5)网格化处理后的CDP道集记录见(图49)。

图49

6)(图50)中下面的信息框中示道集文件的名字为:模型名+GatherEP-GATHER.rec 左面的文件列表中可见为带有卷盘标记的rec文件。

图50

26

7. 第七个作业我们试验由炮记录做CMP叠加剖面 1)激活一个显示面板,输入tgr炮集文件。

2)激活另一个显示面板,输入速度谱文件。做进一步解释效验。 3)激活炮集记录面板,对炮集记录做切除。

4)在主菜单上选RUN,在其下拉菜单上选CMP Stack(SUM) 5)屏幕上出现CMP 叠加会话框即显示在。

在会话框中有三个选项表: “General”, “Advanced” 和 “Velocity”。 点击 “General”显示通用参数: 有三个显示参数格:

? Start Time, Stop time – 以毫秒表示的处理的开始时间和结束时间; ? Time Step –时间步长以 ms为单位。

这些参数是程序根据做模型时所输入的参数记算的,是不可改变的。 ? Start X, Stop X – 时间剖面起始和终止坐标 单位为m; ? X Step – CMP道间距单位为 m;

建议不要改变表中已由程序根据做模型输入的参数自动记算的Start X, Stop X 和X Step参数值 。

6)点击“Advanced”显示参数表:

在“Advanced” 会话参数表中的处理参数:

? Gradient –动校正时延梯度(NMO) 输入一个数值或用滑标确定。 ? Normalization –CMP 时间剖面道的规格化系数。 (Normalization=0, 1, 2)。值为0 输道的振幅用覆盖次数的和规格化,值为 1 –用覆盖次数的均方根规

27

格化,值为 2 –用覆盖次数的四次方根规格化。

7)点击“Velocity” 如果是用速度谱分析得到的速度,则该会话框将如下所示,显示全部的 CMP 叠加速度。

点击OK 键,即可产生 CMP时间剖面。

如果没做速度谱,且速度曲线必须用手工输入,则显示的是一个空白的对话框:

你可以从一个 text 文件把所需要的值输入到此空白会话框内。该程序所用的值时间以ms为单位,速度 – 以m/s为单位。

8)输入完所有需要的参数,点击 OK 键,然后等待处理结束。如果你需要对所产生的结果做进一步处理,将有提示语问你:你想保存前一步产生的结果吗?

图51

28

所得到的CMP 时间剖面(图51)。

8. 第八个作业我们试验由炮记录做DMO叠加剖面(DMO) 1) 激活一个显示面板,输入tgr炮集文件。 2) 激活炮集记录面板,对炮集记录做切除。

3) 在主菜单上选RUN,在其下拉菜单上选Dip Moveout Stack(DMO)。 4) 屏幕上出现DMO 叠加会话框。

5) “General”, “Advanced”参数的含义和填写方法与CMP叠加相同。

6) “Velocity”必须使用平均速度。在空白会话框的条件下用手工编辑的方法拷贝、粘贴的方法将速度输入空白表格,并存储备用。

7)当点击表中的Edit键时,即出现下面的速度编辑表。

8)当点击表中上方的Add即出现X 坐标输入表。

29

9)在X右侧的空格内输入速度点的X坐标值点OK。速度点的X坐标值就加到 编辑窗口的X表的序列中。然后点选Activ edit,小白框内出现对勾,即可编辑时间、速度对。将你准备好的平均速度文件TXT 格式,用拷贝、粘贴将数据置入表中右下方的空格中。注意这里的时间单位是秒,速度单位是米/秒。

10)输入完所有需要的参数,点击 OK 键,然后等待处理结束。

30

52图

所得到的DMO时间剖面(图52)。

9. 第九个作业我们试验由炮记录做克希霍夫叠前偏移 (PSM)

该程序不用做CMP时间剖面和叠后偏移,而是直接从炮集做动校正偏移。为此偏移,必须输入平均速度。

为做叠前偏移,选择“Edit” → ”Pre-Stack Kirchhoff Migration”。叠前偏移会话框出现。

在会话框中有四个列表键:“General”, “Advanced “Velocity” 和, “Aperture”。 在 “General”表中指示:

? Start Time, Stop Time – 处理的开始和结束时间单位为ms; ? Time Step – 采样间隔 ,单位为ms。 这些参数是不改变的。

? Start X, Stop X – 偏移剖面的起始和终止坐标。单位为m; ? X Step –偏移剖面的道间距单位为m。

建议不要改变这些参数 Start X, Stop X 和 X Step .

在Advanced 表中 所指示的处理参数 (类似于1.3):

? Gradient –(NMO) 动校正时延梯度,输入一个数值或用滑标确定。 ? Normalization –偏移时间剖面道的规格化系数。 (Normalization=0, 1, 2)。值为0 输道的振幅用覆盖次数的和规格化,值为 1 –用覆盖次数的均方根规格化,值为 2 –用覆盖次数的四次方根规格化。

在 “Velocity” 表中输入偏移速度值,该值等于平均速度。

31

在下一个版本中计划由模型本身生成速度值(平均速度)。

Aperture”表是为了输入偏移孔径,描绘介质区域大小,在其中保留来自没一个炮道的波场。孔径的定义好象偏移距,单位为米。从炮点到右侧和左侧,即必须定义一组三个值: T –时间 (ms), Al – 左孔径(m), Ar – 右孔径 (m)。手工编辑方法,类似于速度表编辑的方法。

在每一个节点上左孔径和右孔径值都不能为负值。孔径的大小可根据经验公式定义:L = Z * tg ?,

这里 L – 孔径大小,Z – 目标层深度,? – 目标层倾角,处理时间与孔径大小成正比。

输入所有需要的参数后, 然后点击 OK 键并等待处理结束。

叠前克希霍夫偏移的实例

图53显示的叠前偏移结果不甚令人满意,这是因为,孔径参数的选择不好 这里为了通试程序只选了一个节点参数,T=0 为100,T=7 为3000。

32

53

由于时间的关系,其他处理自己通试。发现问题有待专家授课时解决。下面展示专家处理的结果。图54、55、56、57、58。

图54是由中国专家给出的模型

图55是由乌克兰专家根据上述模型做的炮集记录,然后做切除、叠前深度偏移,加时间能量克希霍夫算子得到的偏移剖面的时间域显示。绿色箭头标出的位置同

33

54

图55图

相轴发生畸变,这是由于以上层位的速度变化引起的,而这样的畸变在深度域剖面上就不存在 。

图56

图56由乌克兰专家根据上述模型做的炮集记录,然后做切除、叠前深度偏移,加时间能量克希霍夫算子得到的偏移剖面加上了波形的时间域显示。

图57

图57是由乌克兰专家根据上述模型做的炮集记录经过切除,叠前深度偏移加时

34

间能量克希霍夫算子的深度域显示剖面。

图58

图58是由乌克兰专家根据上述模型做的炮集记录经过切除,叠前深度偏移加时间能量克希霍夫算子的得到剖面。显示方式为带有波形的深度域显示剖面。

图59

图59是由乌克兰专家根据上述模型做的炮集记录经过切除,叠前深度偏移加时间能量克希霍夫算子的深度偏移剖面。与原模型边界线叠合,位置和振幅精度都相当吻合,绿色箭头所示之处两边速度无差异,所以就没有反射能量。

35