折射静校正.docx

折射静校正.docx

《折射静校正.docx》由会员分享，可在线阅读，更多相关《折射静校正.docx（12页珍藏版）》请在冰豆网上搜索。

折射静校正

目录

一、Fathmom：

计算折射静校正的基本流程及批量拾取

1.计算基本流程

2.初至批量拾取

二、GII：

观测系统设置及质量控制

1、GII的窗口环境及鼠标响应

2、编辑功能（Edit）

3、对排列表的操作

4、各种质量控制图件的显示

5、实用程序

6、接口

三、BIO：

地震数据文件格式转换

四、Picker：

初至拾取

五、Branch：

折射界面的定义

六、Fathanal：

折射界面分析

七、Fathmodl：

近地表模型及折射静校正计算

一、Fathom：

计算折射静校正的基本流程及批量拾取

1、用Fathom计算折射静校正的基本流程

1）GeoScribeII

置好观测系统，特别是变观部分一定要置正确。

高程和井深（如果大的话）一定要输入，井深较大时最好用低降速带速度计算出井口时间供校正用。

炮、检点的X，Y坐标最好用实测坐标而不要用内插值或任意给Y值，否则容易造成就结果误差较大。

用plot功能检查各项参数。

2）BIO　数据转换

准备将SEGY数据文件格式转换成cpt内部格式文件。

用BIO生成绿山的.cpt内部格式数据文件时需要有能正确反映数据格式的描述文件（.tdf）。

由于所提供SEGY文件各各有异，应先用头块功能卸出搞清差异再对.tdf文件作编辑。

注意：

道头字里面的野外记录号（FFID）必须与GII数据库里置的记录号（RecordNo.）一致，绝不允许有重复号。

如果用SEGY文件生成观测系统，则必须在道头中置入检波点桩号并编辑相应的.tdf文件。

3）Picker模块：

拾取初至

在“ShotMap”窗口中活化要拾取的炮（把“Shotincreament”改成１即可活化所有炮），此窗口中也可用来删掉不好的预测值。

按一下“PickShot”按钮就可选择要拾取的单炮并显示在“Pick”窗口中。

在“TracePrep”窗口中可调整显示时窗，增益等参数并可以对数据作自动增益控制、滤波、相位旋转等简单处理以提高初至的信噪比，便于拾取（可视需要选择炮检距拾取）。

按“Vel”按钮可让炮集按绝对炮检距重排便于鼠标确定线形动校正速度。

按“PickOpts”按钮选择拾取参数。

算法——大部分情况用“Vel.tracking”选件。

拾取修正——拾取波峰、波谷还是零交线。

缺省为“none”，此时按“atright”按钮，小红“＋”会自动排成一水平横线。

确定后不得在拾取中途变更。

移动选件——用于加线性动校，加折射静校正值，计算基准面静校正及加基准面静校正值。

各拾取按钮的功能：

“zoom”方式，左上角的按钮一个也没被按下。

“AtRight”：

按下此按钮表示鼠标所指处右边的各道被自动拾取。

“AtLeft”：

按下此按钮表示鼠标所指处左边的各道被自动拾取。

“Pick”　：

单道人工拾取。

两个箭头按钮：

分别表示炮号前进一个间隔或后退一个间隔。

具体取决于“Map”窗口中的增量。

利用“Displaysetup”窗口下的“PickGuide”功能可把前一炮的初至拾取时间以一个小短红杠显示在本炮中以指导初至的拾取。

拾取完成后在“File”下退出模块即可（不须“save”），程序自动保存拾取值。

４）Branch定义分支点

启动Branch模块出现两个窗口“CMP”窗口和“BranchPt”窗口；三维项目此时要求给出面元网格的大小（面元的整数倍）。

CMP窗口：

显示整条线的平面位置，“＋”代表一个CMP点，灰色表示无覆盖次数，黑色表示有覆盖。

按下鼠标左键拉开一个范围可选择一定数量的CMP作分析，该段CMP的中心点会出现一个红圈，表示该范围的CMP显示在Pick窗口中。

CMP范围的大小以及分段个数视低、降速带的横向变化情况而定，复杂的话，可多选几个点。

Pick窗口：

在Pick窗口中用鼠标沿点子的趋势确定折射层的启始和最终炮检距。

注意以下几点：

●如果点子出现很散或出现趋势类似的平行线，说明拾取有误或定义观测系统有问题，返回GII或Picker进一步检查。

●发现存在一个以上折射界面时，定义第一层的炮检距范围与第二层的不要重叠，要分开一个点。

●Branch允许最多定义十个折射层，但实际上计算静校正值一般一个折射界面即可，最多二个界面（三层结构）。

●二维数据一般选一组CMP分析已可以，三维数据至少要分析三个以上的CMP组才能实现内插。

模块新增的“Rejection”窗口中可计算并显示拾取时间的残余量（Residuals）供了解拾取值的质量并可把超出门槛的拾取值排除在建立的层位表之外。

５）分析折射初至（Fathanal）

启动本模块。

一旦选择好.mas文件后，模块即自动加载层位表以及初至拾取值等数据。

等“AnalysisOption”窗口打开后，注意观察模块建议你要做的工作（Recommendation）。

一般是先算速度然后分析延迟时。

具体步骤：

●在“Analysis”下拉菜单中选择速度分析可打开折射速度分析框，然后在所提供的三种速度分析算法中选取一种并按“OK”。

●在随之打开的“VelocityAnalysisOptions”框中给出预计的速度范围（Branch已有了解）以及炮数增量（取决于已拾取的炮数以及初至质量）。

按“OK”后，速度计算立即进行并很快显示出来，速度结果须经过编辑和平滑，使之更符合地质概念后才能用于延迟时计算。

下一步为延迟时分析。

●在“Analysis”下用鼠标选“Delaytime”。

●选择“iterativeanalysis”，按“process”。

●对接着出现的“iterativerefratorananlysis”窗口选择当前分析的折射界面

●如果不是初次分析延迟时，或是修改速度后想重算延迟时，要按“Deletecurrentdelay-timesolution”按钮把自己有的延迟时删掉。

●迭代次数以２—４次为宜，高于４次耗时且无意义。

计算结束后延迟时窗口中即显示炮、检点的延迟时曲线，必要时做编辑。

注意：

如存在没有拾取的炮，则没有拾取的炮就没有延迟时值，需要对延迟时进行内插以便填补空缺值。

如果在“Branch”中已定义了多个折射界面，可逐个界面进行，没有次序的限制，分析结束退出分析模块时可自己决定保存或不保存算出的延迟时和速度。

一般按缺省文件保存，模型模块自动加载；自己给文件名保存，供多方案效果比较。

6）反演近地表模型并计算静校正值（Fathmodl）

●在打开项目.mas文件后，即出现“ModelBuild”窗口，由此可以选择速度和延迟时的来源以及确定生成几个折射界面的模型。

注意：

要分清layer和refractor之间的含义区别。

分析模块中的“reftactor”指的是底界的形态，而模块中的“layer”则指的是有速度和时间厚度的层段，一定要分清两者的不同。

●加载信息后，出现“layerSpecification”窗口，要求你提供第一层（layer1）的速度（低速带速度），这就是为什么要求即使不做低速带调查但也要打几口微测井控制低速带速度的原因。

●低速带速度可按常数或变速给出，按“OK”后，模块很快完成并给出模块已完成的信息。

●反演出的界面通常比较粗糙，而实际上远比它光滑的多，故此时务必对反演出的折射界面作平滑。

平滑应自浅层向深层逐界面进行（多界面），模块会相应调整与速度、延迟时的关系。

注意：

静校正值在Fathmodl中是以模型为基础计算出的。

地质上越是合理的模型，得出的静校正值也越正确，微测井分层及速度也为你评价模型提供了依据。

●修改层参数可调整模型，关键是低速带的速度。

比如，发现低速带底界的形态与地表高程线恰好相反，很可能所给的低速带速度太高，需要适度的降低并观察模型的变化。

●模型满意后，在“model”菜单下选择“staticsoptions”可根据需要按浮动基准面或按最终基准面计算静校正值。

通常是：

给出最终基准面和替换速度，OK；静校正计算瞬时完成并显示在窗口中；在“File”菜单下可存入GeoScribe数据库或存入ASCII文件，通常都存一下，GeoScribere中的静校正值供调“Picker”模块检查静校正值的质量，满意的话即可将已存入的ASCII格式的静校正文件拷出。

以上为用Fathmom计算静校正值的基本流程，具体操作请参看各模块说明。

需要着重强调的几点是：

a）修改单炮初至拾取值后，一定要再建（rebuild）层位表

b）修改速度场后一定要重算延迟时，但重算延迟时前要把原来延迟时删掉

c）拾取值问题较大时（“Branch”中可发现点子分散、乱、飘）可在“Branch”中计算“Rejection”以便用门槛值将拾取有误的点子剔除在外，保证层位表的质量。

d）“Picker”中拾取初至时一定要注意不要把直达波、反射波全部拾取近来，这样会严重破坏计算效果。

注意：

只拾取折射初至。

2、批量拾取：

1）批量拾取流程图

2）

批量拾取的注意要点：

a）在Picker里按“Batchpcik”按钮时，这时将修改预测拾取值而不涉及人工拾取的那些空炮（控制点），故人工拾取初至时要认真、仔细。

批量拾取结束时，模块会告诉你批量拾取供拾取了多少炮（人工拾取炮除外）。

b）检查批量拾取的质量可以通过Picker里的“MapOption”进行，给一个增量，比如是2，即每2炮检查一炮，如果质量满意，可退出Picker。

c）如果数据信噪比很低，会发现批量拾取的炮中有相当数量不行，需要把某些炮的预测拾取值删掉（在“ShotMap”窗口中的“MapOption”对话框中），再在删掉的范围中人工多拾取几炮加强控制。

接着在Branch中修改层位表，返回Fathanal重算速度和延迟时并重新生成预测拾取值。

这时需进入Picker再运行批量拾取并检查那些原来有问题炮的拾取质量。

d）拾取没有问题后，有两条路径可以选择：

●调整分支点，并重新生成层位表，因为原来的层位表是依据每隔几炮拾取一炮而生成的，故炮全部拾取后必须加以调整。

如果第一次分析得到的延迟时已经存在，则可在Branch中加上延迟时使拾取时间点更集中，这对近地表条件复杂或高程变化大的地方特别有用。

●不调整Branch而直接进入Fathanal，但模块初始化时一定要选择“Rebuildthelayertables”（重建层位表），另外，分析时不要再用原来算出的速度和延迟时，必须要逐层分别计算速度（平滑）、延迟时。

分析结束后，在“File”菜单下选择“Save”，把各层的速度和延迟时以二进制缺省文件保存起来。

二、GII：

观测系统设置及质量控制

目的：

为观测系统的输入及质量控制提供灵活、直观、简便易用的工具，具有以下特点。

●支持多窗口显示，显示内容与排列表内的值密切相关，可快速发现错误，加以纠正；

●QC图件可直接形成绘图文件输出到打印机；

●排列表具有多种输入及编辑功能，使观测系统的设置简单易行；

●允许多个ASCII文件同时输入；

●与多种常用的处理系统有输出接口，便于数据直接存取；

●可按标准工业格式如：

SEG—P1,UKOOA以及ASCII列表输出观测系统的有关信息。

1、GII的窗口环境及鼠标响应

a）G2提供下述几种窗口

图表窗口——用于设置观测系统

输入窗口——交互定义要输入的文本文件

显示窗口——以图方式检查观测系统，包括高程、工区位置、迭加图（二维），面元属性（二维弯线及三维）等。

注意：

输入窗口（import）只在输入时出现；排列表窗口每一种只允许出现一个，而面元属性窗口和面元统计窗口可同时打开多个以供比较。

如果想使窗口中的内容不要刷新，将鼠标置于该窗口中然后按一下“ESC”键。

b）打印功能

在File菜单下选择Print可生成活化窗口的屏幕卸出供打印

●在装有Windows的微机上，可将屏幕卸出直接送到打印机

●在UNIX工作站上，生成PS格式绘图文件。

在微机上打印时，步骤如下：

PrintSetup选择打印机options：

选择打印机

驱动程序OK，关闭setup窗口OK，打印

在UNIX工作站上，步骤如下：

选择图件摆放位置（横打：

Landscape或竖打：

Portrait）

选择纸页的大小（见手册39页）

选择比例尺

选择背景图象为彩色、灰度或黑色

在File下面要给出输出的绘图文件名，缺省名为g2.ps

用lpr命令把生成的PS格式文件送到PS打印机去。

例如对PH打印机则用命令：

lpr—phpg2.ps

c）鼠标与响应

窗口

鼠标动作

响应

排列表

“shift”+双击左键

转入平面图中响应的炮、检点位置

平面图显示，高程等

“shift”+单击左键

识别炮、检点

平面图（只适用三维）

“ctrl”+单击左键

面元属性的编号

面元属性显示窗口

“shift”+单击左键

识别面元

“ctrl”+单击左键

显示蜘蛛图

面元属性表窗口

“shift”+单击左键

识别面元

面元统计窗口

“shift”+单击左键

给出统计报告（具体内容视当时显示选择而定

2、编辑（Edit）功能

对于排列表中不同类型的空格分别提供不同的编辑功能，大部分编辑功能是面向列的，因为一列正好对应于某个数据域。

下表总结了各种编辑功能的应用。

编辑功能

格的选择方法

cut,copy,paste

要贴上（paste）数据的格数（行数x列数）必须恰好等于cut或copy时所选的格数。

删除行

单行或多行

排列表扩展

检波点排列表；station 列炮点排列表；source列

填入

所选取的填入范围顶部一定要有已填入数字的格

公式操作

单列

转向菜单

对活化排列起作用，无需选择

插入行

任何行或一行中的任意一格

内插

一列中进行，但所选范围上下一定要包含已填数的格

算术操作

单列或多列

查询

单列或多列

平滑

单列或多列

重排

单列，即按某一列中的内容进行整个数据表重排

绝对/相对坐标转换

多活化排列表起作用，无须选择

GII有三维内插/外推功能，在做三维空间内插/外推时，列中的全部已有值（包括零）按棋盘格排列，而要填值的空格则位于三角形中，其值按三角上的值算出。

注意：

只有位于已经活化的范围内的空格才填入新值。

新增加的空间平滑算子只适用于三维，其平滑原理类似与空间内插，但要注意的是一旦进行了三维平滑后值再也不能恢复，故建议对要平滑的内容先备份一下或暂时存入空列中以备不测，否则非常容易造成整个观测系统的混乱。

改变排列表中的列的位置及列的宽度

改变列的位置：

将光标指到要变动的那列，按下鼠标左键及键盘上的“ctrl”键不放的同时把该列拉到需要的位置上后放开鼠标键和“ctrl”键。

改变列宽度，把光标指到一列右边，按下鼠标左键和“ctrl”键不放将该列拉到所需宽度后再放松鼠标和“ctrl”键。

3、对排列表的操作

在着手填排列表前，一定要通过“General”菜单填好响应项目参数，如果填好排列表后再改变这些参数有可能丢失数据。