reflex

REFLEXW软件

使用指南

2

一说明

本手册主要用来指导用户如何使用REFLEXW软件，它并不包括所

有的细节。如果需要了解全部细节，请参阅英文版手册。

二安装及打开软件

安装软件时，只需点击tup即可。

另外，必须安装软件狗的驱动程序,安装时也是点击即可。

打开软件时，必须将软件狗插到USB口上，然后从“开始”→“程序”→“REFLEX”→

“REFLEXW”进入。当显示出“Projectdirectory”菜单时，即可使用该软件。（最好在桌面上

建立快捷键）。

进入菜单后，建议选择菜单中的“newproject”,出现“enternameofthenewproject”。输

入任务名（如a1等），点击OK,进入“Reflex_win”主菜单。

点击“Modules”,选择各相关模块。

三输入和第一次显示GPR数据

下面介绍如何输入和显示GPR数据，我们用MALA公司的RAMAC/GPR数据做例子，用

其它数据时仅需改变一些选项即可。

一、输入GPR数据

1.进入2D数据分析(2D-dataanalysis)

2.用选项File/Open/import(文件/打开/输入)进入import(输入)菜单，此时出现REFLEXW

DataImport（数据输入）菜单（见右图）

3

3.输入下列参数：

inputformat（输入格式）：RAMAC

outputformat（输出格式）：16bitinteger（16位整数）

filenamespecification（指定文件名）：originalname（原始名，举例）

选择X或Y作为剖面方向（ProfileDirection）及Y或X作为剖面常数（ProfileConstant）

选择道增量（traceincrement）和/或坐标（coordinates）是否从原始数据读出（从

ControlOption）

4.点击选项ConverttoReflex，出现一个文件打开菜单，并在任务目录下出现子目录ASCII。

你可以从该输入路径或从任何其它路径选择RAMAC文件（RD3或RAD文件），任何

情况下必须有RD3和RAD文件。选择好需要的RAMAC文件后，文件转换成REFLEXW

内置格式并保存在路径ROHDATA下。用选项PrimaryFile（第一文件）激活的输入数据

自动显示在第一窗口中

5.选择Exit退出

二、显示数据

1.当做完输入后,它将用标准绘图选项来显示数据,可以用选项Plot/Options(绘图/选项)来改

变绘图选项,点击该选项出现PlotOptions菜单(见右图)

2.用于GPR数据的主要PlotOptions包括:

-Plotmode(绘图模式)

-PointmodeScale(点模式比例)

-EnergyDecay(能量衰减)

4

-AmplitudeScale(振幅比例)

Plotmode(绘图模式):缺省时使用Pointmode(点模式)

PointmodeScale(点模式比例):

XYScaledPlot(XY比例绘图):数据完全显示在实际窗口中,两个比例选项Xscale和Yscale设

为1。该选项可用于小数据（几道）或将整个数据显示在主菜单中（可用Zoom选项），当

不做Zooming（缩放）或数据很大时，显示分辨率非常差。

PixelPerSample（每样点像素）：每个数据点的大小显示在给定的屏幕像素点中，只能在y方

向用缩放或重新改变比例。对大型数据该选项非常有用。

PixelPerTrace（每道像素）：两道间的距离给定在屏幕的像素点中，根据实际窗口尺寸绘制

每道的整个时间序列。该选项不能用缩放，如该选项可以用在大型数据中（很多道）。

EnergyDecay（能量衰减）：RAMAC/GPR数据是原始数据，它在时间方向没有任何增益，因

此为了补偿能量随时间的衰减，在没有使用增益滤波时，可以点击该选项。

AmplitudeScale（振幅比例）：该参数控制数据显示的对比度。

上图是上述参数显示的图象

注意：通常数据必须进行滤波，主要滤波包括：

（Staticcorrection）静校正补偿直达波的时间延迟

（y-gain）y（时间）增益

（dewow）抽取平均道

四用REFLEXW滤波

下面介绍用REFLEXW进行滤波，它包括单滤波（第二章）和批处理（第三章）。

一、概述

1.滤波选项仅用于REFLEXW格式数据，它意味着其它格式数据（像RAMAC，GSSI，

SEGY等）必须首先用打开File/import输入数据，输入数据保存在子目录rohdata

中。

2.滤波后的数据用相同的文件名和由处理标签给定的扩展名保存在各文件下，如：原

始文件：;滤波后的文件:profile1.00t(处理标签设为0)。滤波后的数据总

是保存在子目录procdata中。

5

3.滤波可以用于单数据剖面或许多剖面（所谓的批处理）。

二、单数据文件滤波

下面详细介绍单数据剖面滤波，滤波顺序用以下选项：

抽取平均道（subtract-mean（dewow））

静校正（staticcorrection）

手动增益（y）（manualgain（y））

背景去除（backgroundremoval）

它基本是GPR数据的标准滤波顺序，关于各滤波的步骤的帮助信息，按F1在滤波菜单中看

到。

1.进入模块二维数据分析（2D-dataanalysis）

2.用选项File/open/rawdata将数据文件载入

(load)到第一窗口,然后选择要处理的数据

文件。

3.选择的数据文件显示在第一窗口（用水平

分割），根据分割设置选项（见

选项Plot/options-选择垂直分割

或水平分割）,屏幕将分成垂直或水平的两个窗口(第一窗口和第二窗口),如果两

个选项都不选,屏幕将只有一个窗口。

4.用选项处理/1D滤波（Processing/1D-Filter）进行第一步处理

5.出现1D滤波窗口（1D-filterwindow），输入下列参数或选项：

-点击选项抽取平均道（subtract-mean（dewow））

-进入滤波参数时间窗口（timewindow）

-进入需要的处理标签（ProcessingLabel）

-开始（start）滤波

选则选项应用到举例道（applyonexampletrace），可以

控制输入滤波参数后的效果（滤波后的道filteredtrace和滤波后的频谱filteredspectrum），

道数（tracenumber）用来决定要显示的道。

点击选项批处理（SequenceProc.），出现批处理菜单，我们将实际滤波步骤加到处理流程中。

6.滤波后的数据显示在第二窗口中,关闭(clo)滤波窗口,必须进入选项

File/ChangeSecondToPrimary(文件/将第二窗口改到第一窗口),以便把滤波后的数据用做

6

下一步处理。

7.用选项处理/静校正/切除（Processing/staticCorrection/muting）,进入第二步处理。

8.在静校正/切除窗口（StaticCorrection/muting）内选择静校正（static/correction）和子选项

移动到负时间（movetonegativetimes）。

9.出现右图的列表，你可以在该表中手动输入校准值，也可以用鼠标左键从图像中输入该

值。选项保存（save）保存校准文件，并用于以后各步骤中。

10.输入想要的处理标签（Processinglabel），按start开始静校正。注意：处理标签必须与

第一步不同，因为原始和处理后的数据文件文件名不能相同。

11.静校正后的数据显示在第二窗口，关闭（clo）静校正窗口，使用选项

File/ChangeSecondToPrimary，以便做下一步处理。

12.用选项处理/增益（Processing/Gain）做第三步处理。

7

13.在Gain窗口内选择手动增益（y）（manualgain（y）），出现一个列表用来输入增益值（类

似于以前的静校准），同时有一个新窗口用来交互输入增益值。

14.输入想要的处理标签（Processinglabel），按start开始手动增益。处理标签必须与以前

文件不同。

15.手动增益处理的数据显示在第二窗口，为了做下一步处理，必须关闭（clo）增益窗

口进入选项File/ChangeSecondToPrimary。

16.用选项处理/2D滤波（Processing/2D-filter）进入最后滤波步骤。

17.在2D-filter窗口内选择背景去除（backgroundremoval），此时用整个数据范围决定滤波

道（子选项整道wholeline），输入想要的处理标签（Processinglabel），按start开始手动

增益。处理标签必须与以前文件不同。

8

18.滤波已完成，你可以用选项File/ChangeSecondToPrimary，然后在plotoption菜单内关

闭选项,这样只显示最终滤波后的数据。整个处理流程用选项文件/处理流程

(File/ProcessingFlow)显示。实际文件的处理流程也可以用单独的文件保存，并用于以后

的批处理中。

三、用相同的批处理对其它数据做滤波

9

批处理可以对选定的许多剖面自动做滤波处理。使用时必须先定义批处理步骤，然后选择数

据线，最后开始进行批处理。

1.根据单个文件的滤波定义批处理

为了该目的,当进入任何滤波菜单时，你必须点击选项批处理(SequenceProc.)。打开批处

理菜单（SequenceProcessingmenu），用选项加入实际处理（AddActualProc.）,这样可以

把实际滤波步骤加到处理流程中。AddActualproc必须在指定滤波参数后进行，事先要

关闭实际滤波菜单。

做完所有的要滤波的步骤后，必须用选项保存批处理（SaveSequence）保存实际的批处

理。

该处理顺序文件可以在任何时间用载入批处理（LoadSequence）来载入。

通常实际处理步骤在处理顺序的最后加入，但也可以在任何位置加入。此时必须点击相

应的处理步骤，然后用选项删除处理（DeleteProc）来移去该处理步骤。

2.开始批处理

可以在处理单个文件时进入批处理菜单，也可以用选项处理/批处理

（Processing/SequenceProcessing）。

定义处理流程后，必须指定处理标签（ProcessingLabel），然后开始批处理。

点击选项StartActualLine将处理流程应用到实际原始文件中。

10

点击选项StartSeveralLines将处理流程应用到选定的许多文件中。

选项文件滤波（filefilter）来进入多文件选项的滤波，文件路径（filepath）指定要载入

文件的路径。

点击选项显示数据（plotdata）可以瞬时显示原始和滤波后的数据，它可以直接控制给

定的结果。

可以选择所谓的批模式或单处理模式（选项单处理SingleProcess），将各处理步骤加到

原始文件中。

五、滤波详述

1．软件产生的子目录：

ASCII：保存原始数据（如SEGY，RD3等）。

ROHDATA：保存转换后的原始数据。

PROCDATA：保存处理后的数据和时间切片。

1.11处理主菜单

1.11.1一维滤波（1D-Filter）

1.11.1.1平均值滤波（Meanfilter）

该滤波分别作用于每一道，它对选定的时间样点数做平均。

滤波参数平均范围（meanrange）是以样点形式决定窗口宽度，它是对各道的振幅值做平均。

如窗口宽度是4，做平均时计算前面的两个样点和后面的两个样点，对它们取平均后把平均

值放在该点。

该平均值的效果类似于时域中的低通滤波。它比低通滤波快得多，但没有低通滤波有效。

1.11.1.2中值滤波（Medianfilter）

该滤波分别作用于每一道，它对选定的时间样点数取中值。

滤波参数中值范围（medianrange）是以样点形式决定窗口宽度，它是对各道的振幅值取中

值。如窗口宽度是4，做中值时用前面的两个样点和后面的两个样点，取它们的中值放在该

点。

中值滤波可以从数据中去除毛刺，选择的窗口越宽，中值滤波的低通效果越明显。

1.11.1.3带通滤波（Bandpassfrequency）

该滤波分别作用于每一道，它是在频域中对各道做带通滤波。

频带由设置的四个频率值指定，第一点定义低切频率（low-cutfrequency）,第二点是平台的

开始（lowerplateau）。低切频率和平台开始之间，滤波是余弦窗口。第三点是平台结束（upper

plateau），第四点是高切频率（highcutfrequency），它们之间的滤波也是余弦窗口。低于低

切和高于高切的频谱都设为0。通过选择各点值，也可以近似实现低通或高通。

频率值可以手动设置，也可以从原始样点道的频谱中选取。

当噪声信号的频率范围与有效信号不同时，带通滤波可以抑制噪声。滤波参数设置道时，建

议它既包括噪声的主要部分，也包括信号的主要部分。

1.11.1.4巴特沃斯带通滤波（Bandpassbutterworth）

该滤波分别作用于每一道，它是在时域中用递归滤波进行带通滤波。

11

滤波带通过设置两个频率值指定，第一点决定低切频率（low-cutfrequency），第二点决定高

切频率(highcutfrequency)。低于低切和高于高切的频谱都设为0。通过选择各点值，也可以

近似实现低通或高通。

频率值可以手动设置，也可以从原始样点道的频谱中选取。

当噪声信号的频率范围与有效信号不同时，带通滤波可以抑制噪声。滤波参数设置道时，建

议它既包括噪声的主要部分，也包括信号的主要部分。

1.11.1.5滤波/时间相关（Filter/timedependent）

该滤波分别作用于每一道，它是在时域中用递归滤波进行与时间相关的带通滤波。

滤波带通过设置两个频率值对指定。第一个频率值对时间开始有效，第二个频率值对时间结

束有效。每对的第一点决定低切频率(low-cutfrequency),每对的第二点决定高切频率(high

cutfrequency)。低于低切和高于高切的频谱都设为0。根据输入的各对值（线形插值），对每

一时间值计算新的滤波边界。使用这些计算的滤波边界进行巴特沃斯滤波。

例：由于深部高频噪声较强，如果你想随时间增加而增大低通滤波：

设置第一滤波对导致大带宽（如f1=100MHz,f2=500MHz）.

设置第二滤波对导致小带宽（如f1=100MHz,f2=300MHz）.

使用这些值，浅部的高频部分保存下来，而深部的高频信号被滤掉。

与巴特沃斯带通滤波一致，频率值可以手动设置，也可以从原始样点道的频谱中选取。

1.11.1.6陷波/频率（Norchfilter/frequ.）

该滤波分别作用于每一道，它是在频域中对每一道进行陷波。

滤波带通过设置两个频率值指定。第一点决定低通平台(lowerplateau)，第二点是低切频率

(low-cutfrequency)。低频平台的终点和低切频率之间，滤波是余弦窗口。第三点是高切频率

（highcutfrequency），第四点是高频平台开始（upperplateau），它们之间的滤波也是余弦

窗口。切除频率之间的频谱都设为0。

频率值可以手动设置，也可以从原始样点道的频谱中选取。

陷波特别适用于抑制近似单频噪声。处理时应选择包括很强噪声的道，这将易于设置滤波频

带。

1.11.1.7反褶积（Deconvolution）

该滤波分别作用于每一道，该滤波允许对剖面应用反褶积，这就是所谓的毛刺反褶积，即它

可以创立较高的带宽和较平滑的谱（根据毛刺的谱）。该程序采用Levinson递归算子（Wiener

滤波法），对每一道予先设定的样点数做相关计算，用做Levison算子的输入信号，来计算

递归的反褶积滤波道。然后滤波道褶积为原始道。

反褶积需要四个参数（自相关范围的开始和结束，滤波道长度，白噪声百分比）。用自相关

范围开始的指标（）和结束的指标()，可以控制反褶积效果。当需要

对深部信号进行反褶积时，自相关范围必须取大的时间范围，而小数值范围仅在需要将原始

信号转换成一个毛刺时有效。

滤波道长度(filterlength)通常取小的数值，它必须小于或等于自相关范围。

为了稳定反褶积处理，必须把白噪声（noi（perc.））加到数据中。输入参数噪声（百分比）

决定白噪声。数值太大会导致反褶积执行不充分，数值太小会产生稳定性问题。

1.11.1.8反褶积/成型

该统计反褶积是应用Toeplitz结构的线性方程中的Levinson递归算子。它分别作用于每一道。

对每一道，对预定义的样点数做自相关，以及对输入和期望输出做自相关计算。

12

这些相关是输入Levinson算子，然后递归计算滤波道反褶积，然后该滤波道褶积成原始道。

根据希望选择输出的反褶积的目的，可以使用以下反褶积：

毛刺-反褶积

毛刺反褶积的目的是集中时域中的信号子波，或在频域的反褶积道中使频谱更宽和更平滑。

我们期望输出的是一个毛刺，即一个能量最集中的子波。如果选择0延迟毛刺（0lagspike），

能量集中在子波开始处；如果选择毛刺（spiking），则能量集中在右参数延迟（lag）定义的

样点处。实际上输入信号在混合相位处。由于毛刺反褶积使振幅谱变平坦，反褶积会造成高

频噪声，因此需要在反褶积后做带通滤波。我们通过选择参数滤波谱（filterspectrum）来进

行处理。

子波成型

子波成型的目的是形成输出信号的子波，即改善输出信号子波的最小相位特性。输入信号子

波的波形可以转换成自由选择的子波，如Richer子波。反褶积类型零相位（zeropha）和

最小相位（）可以将输入信号子波转换成零相位或最小相位子波。实际上，两种

方法相差很小，由于转换成最小相位子波需要较长的运算时间，我们一般选择零相位反褶积。

预测反褶积

预测反褶积的目的是抑制多次震荡或反射。为了抑制多次震荡，根据多次震荡双程走时选择

延迟；为了抑制多次反射，选择较小的延迟。

如果各道之间的信号波形变化很大，反褶积会降低反射的空间一致性，它是由于各道间自相

关的强烈变化引起的。为解决这一问题，采用参数平均自相关（meanauto）来计算反褶积

滤波。所有各道用相同反褶积滤波做反褶积，空间一致性也会保持。

反褶积需要三个参数（自相关范围的开始和结束，滤波道的长度），用自相关范围开始

（）和()可以控制反褶积的效果。为了抑制多次震荡需要选择大的时

间范围，当需要把原始信号转换成毛刺时需要小的范围。滤波道长度（filterlength）通常很

小，它必须小于或等于自相关范围。延迟（lag）必须小于滤波长度。参数,

,filterlength和lag用实际时间尺寸给出。

1.11.1.9抽取平均道（subtract-mean）

该滤波分别作用于每一道，用它对每道的每个值做滑动平均计算，该滑动平均值从中心点中

减去。

计算时必须输入滤波参数时间窗口（timewindow）值。

该滤波用于消除低频部分，因此窗口范围应设为大约一个主周期。

1.11.1.10去直流漂移（subtract-DC-shift）

该滤波分别作用于每一道，它是所谓的零平均，即去除时间常数漂移。

必须输入两个时间值（和），在该时间范围内对每道的各样点做平均，然后减去

该平均值。时间一在时间窗口的约三分之二处，时间二在离最后样点约5个样点处。

1.11.1.11互相关（crosscorrelation）

该选项用于计算指定时间道与剖面所有其它道的互相关。

通过选择由滤波参数给定的道/时间来做互相关，互相关道在参考道窗口中显示。第一个滤

波参数是开始参数（），第二个是结束参数（）,第三个参数给出开始

13

距离范围（ce），第四个是结束距离范围（ce）。它对选定道的选定

时间范围进行自动平均。如选定的道/时间范围可以是直达波范围。

互相关用于改善信噪比，它的另一作用是缩短原始波列（当环境干扰较大时）。

1.11.1.12样点重组（resampling）

该滤波用来在时间方向重组样点，它是在频率范围内进行。

滤波参数timeincr.在时间方向上决定新的采样间隔。

该滤波可以把用不同时间间隔采集的许多剖面变成等时间间隔，注意：降低时间间隔并不能

增加垂直分辨率。

1.11.1.13Walsh带通（Walshbandpass）

该滤波分别作用于每一道，你可以用Walsh变换对每一道进行带通滤波。该变换是基于阶跃

函数的标准系统而不是用于付立叶变换的正弦函数。因此它更适用于渐变数据（如时间切

片）。

滤波时要指定Walsh系数的低切（lowercutoff）和高切（uppercutoff），最大Walsh系数等

于样点数的一半。锥度可用于低切和高切系数（参数lowertaper和uppertaper）。

1.11.2增益（Gain）

1.11.2.1AGC增益（AGC-Gain）

该增益是用预先指定的窗口在y方向对信号振幅进行放大。

滤波参数窗口长度（windowlength）决定在给定方向的窗口尺寸，该程序首先计算每一道的

整个时间范围的平均振幅，然后按比例对每个振幅值进行放大，在选定窗口内使平均振幅具

有相同的值。窗口尺寸决定放大效果，窗口尺寸只包含一个样点时，意味着每个样点的振幅

都相同；窗口尺寸包含整道时表示没有放大。小窗口尺寸使各点振幅较均衡，大窗口振幅变

化大。

该选项是对信号的振幅放大。

1.11.2.2能量衰减（Energydecay）

该选项是根据平均振幅衰减曲线对y方向进行增益。

首先从所有现存的道中决定平均衰减曲线，然后对每一点的值用衰减曲线的值去除。

所有数据点都乘以一个比例因子（scalingfactor），你可以用以下方式选择比例因子：如果文

件格式是16位整数，那么所有振幅都在16数据位内；否则那些振幅自动限制在32000或

-32000内。

与AGC增益类似，它是对信号的振幅进行放大。

1.11.2.3增益函数（Gainfunction）

该滤波是对个点数据分别乘以给定的函数g（y）或g（t）。

函数g（t）包括线性和指数部分：g（t）=（1+a*t）*exp（b*t）其中a=a‘/脉宽

b=b’*v/8.69v=0.1m/ns或1.0m/ms,脉宽由主频得到.

必须输入两个滤波参数a‘（线性增益）和b’（指数增益），a‘没有单位，b’单位是dB/m。

另外必须输入开始时间（starttime）和最大增益（）。

数据乘以该函数是为了补偿衰减或散射造成的损失。

1.11.2.4Removeheadergain

14

该增益没用。

1.11.2.5手动增益（y）（Manualgain(y)）

该选项交互式定义y方向的数字式增益曲线，并将其应用到数据中。

激活该项后，出现一个表格和一个窗口，他们都可以输入增益曲线（time在给定的时间坐标，

factor用db）。选定的值自动输入到表中，反过来也可以到增益窗口中。可以输入20-80db，

在两个值之间线性插值。输入增益值的效果自动显示在滤波后的道窗口中。

选项save用ASCII文件在任务子目录下保存实际编辑的值，文件名用y增益，扩展名任选。

ASCII文件行包括两栏：第一栏是时间；第二栏是因子（factor），单位db。

选项load载入要编辑数值的ASCII文件，列表和窗口可自动更新。

选项ret重新设置当前编辑的值。

第一个和最后一个增益值保存在处理流程（processingflow）中。

1.11.2.6手动增益（x）（Manualgain(x)）

该选项交互式定义x方向的数字式增益曲线，并将其应用到数据中。

激活该项后，出现一个表格，可以输入新的增益曲线。必须输入距离（distance）和增益因

子（factor）（线性），在两个值之间线性插值。可以在剖面希望处理的位置上点机鼠标左键

来决定距离值，乘积因子必须在表中手工输入。

选项save用ASCII文件在任务子目录下保存实际编辑的值，文件名用x增益，扩展名任选。

ASCII文件行包括两栏：第一栏是时间；第二栏是因子（factor），单位db。

选项load载入要编辑数值的ASCII文件，列表和窗口可自动更新。

选项ret重新设置当前编辑的值。

该选项用来对剖面的某一部分进行增益。

1.11.2.7比例窗口增益（x）（Scaledwindowgain(x)）

该增益对指定的y窗口在x方向对振幅进行增益。

选定该项时，程序需要输入时间窗口（开始窗口startwindow和结束窗口endwindow），在

该窗口内对振幅进行归一化。

它既可以固定归一化值（选fixmeanvalueactivated），也可以自动从数据中决定该值（选fix

meanvaluedeactivated）。

用选项fixmeanvalueactivated，程序先在该窗口的所有道计算平均振幅（它意味着振幅也三

个滤波参数中显示）。

用选项fixmeanvaluedeactivated，必须在三个滤波参数内手动输入平均振幅。如果由于地面

耦和原因，其振幅特性不同，你需要比较不同的数据文件时它非常有用（如用时间切片）。

一道的每个时间样点都乘以一个比例因子，它是由每道y窗口内的平均振幅决定及以前计算

和确定的平均总振幅决定。在所有道的y窗口内，其结果是对信号振幅进行增益。

当地面耦合不同时，用该选项进行补偿。

1.11.2.8x-距离衰减（x-distancedecay(db)）

雷达不用。

1.11.3静校正/切除（Staticcorrection/muting）

1.11.3.1静校正(staticcorrection)

15

略

1.11.3.2动校正（Dynamiccorrection）

略

1.11.3.3移动开始时间（movestarttime）

该选项可以移动开始时间。

使用时选择手动输入（manualinput）。在参考道中确定直达波的开始时间t，然后在输入框

内输入-t，这样直达波的开始位置移到0点。

1.11.3.4切除（Muting）

该选项在指定点处将上部数据或下部数据置为0。

用选项上切除（mutedataabovemutingcurve）将指定点上部的数据置为0。用选项下切除

（mutedatabelowmutingcurve）将指定点下部的数据置为0。

激活选项切除（muting）后，将出现带距离和时间参数的列表，它可以交互式输入新的切除

曲线。距离和时间必须按相关的单位输入。在两点之间线性插入各值，在x的边缘处，属于

最小或最大x坐标的指定值将外插。在原始剖面中可以用鼠标左键在选定的位置确定数值，

选定的值自动放到表中，反过来也可加到剖面图像中。

选项应用到整个剖面（applyonwholeline）可以将预先选定的实际切除点重复加到整个剖面

中。

选项保存（save）在任务路径下用ASCII文件保存实际编辑的值，文件扩展名可以任选。

ASCII的每一行由两栏组成：1栏：距离；2栏：时间。

选项载入（load）从ASCII文件中载入要编辑的值，列表和剖面图像自动升级。

选项重置（ret）重新设置当前编辑值。

1.11.3.5时间切除（Timecut）

该选项可以把各道限定在预先设定的最大时间内。

在给定时间维度内的滤波参数最大时间（）决定要限定的各道时间，它不用压缩或

分开图像，只是把时间范围减小到预先设定的范围。

在数据采集时，如果时间范围设置得太大时可以选择该项，这样你可以节省硬盘空间，在随

后的数据处理中加快计算速度。

1.11.3.6校准最大相位（）

该选项可以在给定的时间窗口内自动发现最大振幅的相位，并使该相位变平坦。

滤波参数开始窗口（startwindow）和结束窗口（endwindow）决定时间窗口，用来发现要

校正的相位，极性（polarity）用来决定相位的极性，它在给定窗口内决定要选择振幅的相

位。

参数平均（average（traces））用于做滚动平均，输入的值应大于2。

用该选项可以对直达波做快速静校正或自动改变开始时间。

1.11.3.7校准拾取相位（Correctpickedpha）

它可以对拾取的相位做静校正，使拾取的相位变平坦。

必须在滤波参数组中设置目标时间（targettime），基于平坦相位的时间设置成目标时间。

根据进入的最大时间和最大拾取时间的差自动改变开始时间。

16

用该选项可以使直达波等相位变平坦。

1.11.3.8校准两层（Correctfor2layers）

如果数据是在崎岖不平的地面上方的平坦介质（如雪）时，可以使用该选项。

介绍略。

1.11.3.9校准三维层析成像（Correct3D-tomography）

用该选项做三维层析校准。

1.11.4去剪切/算术函数（Declipping/arithmeticfunctions）

1.11.4.1去剪切（Declipping）

根据采集的数据，你可以选择不同的去剪切方式，以便自动指定哪个数据剪切，哪个数据不

剪切：

当数据文件中剪切掉的值有最大值时，应选择去剪切/最大值（declipping/）。

如果剪切掉的数据有一个平台，即这些信号的所有数据都相同，应选择去剪切/平台

（declipping/plateau）。

如果剪切掉的数据有扰动，此时只能应用去剪切/域值（declipping/threshold）。如剪切的相位

从正相位跳到负相位，可以使用该选项。

如果选择去剪切/最大值和去剪切/平台，当至少有两个样点振幅值相同时，就使用剪切功能。

选择子菜单最大值，加入的最大值必须等于整个剖面中的最大值；选择子菜单平台，给定的

振幅值没有限制；选择子菜单域值，必须输入域值，所有大于或等于该值的数，都剪切掉。

如果数值达不到平台值，或予处理的步骤（如增益函数或AGC）使振幅没有很大改善，就

可以选择该值。

去剪切后，所有振幅值自动乘以比例因子，比例因子的选择原则是：如果REFLEX格式是

16位整数，则振幅值在16位内；否则它自动限定到32000~-32000内。

1.11.4.2算术函数

使用该处理步骤可以对当前道的样点做简单的算术计算。

加（addition）需要输入一个值，并将该值加到每个点上。

乘（mulyiplication）需要输入一个值，并将每个点都乘以该值。

平方（square），平方根（root），自然对数（naturallog）

绝对值（absolute），指数（expx），样点范围设置为0（o0）

1.11.5合成道-分析/谱分析

1.11.5.1合成道-分析(complextrace-analysis)

该滤波是用希尔波特变换来计算瞬时特性如瞬时振幅,瞬时频率和瞬时相位.

瞬时振幅(envelope)计算反射强度,它正比于信号能量的平方根,瞬时振幅给出各道的能量

特性,另外它可以决定直达波信号。

瞬时相位（instantaneouspha）是计算剖面的连续性，瞬时相位的跳跃表示到达信号的叠加。

该相位可以用选项相位/常数（pha/cont）或相位/2倍象素（pha/2pi）连续绘出。显示是

最好用波形方式。

瞬时相位的时间率的改变是瞬时频率（instantaneousfrequency）。瞬时频率的计算是用相位

的中心差，瞬时频率的绘制可以分辨直达波的叠加。此时波形显示比灰度显示更有用，瞬时

频率的变化等级更高。

17

1.11.5.2道频谱（tracespectrum）

道频谱可以计算各道的振幅谱。

可以将振幅谱计算限定到指定的时间范围，此时需要输入滤波参数组开始时间（starttime）

和结束时间（endtime），缺省时开始时间设为各道的开始时间，结束时间设为各道的最大时

间。

Y轴标签自动改为频率。

1.11.5.3滑动窗口谱（Movingwindowspectrum）

滑动窗口谱对所有道的所有平均时间样点做滑动窗口频率分析。

对每个时间值，该滤波对指定的窗口做振幅谱（滤波参数窗口宽度windowwidth），窗口用

每n个样点的参数沿y（时间）轴连续滑动（滤波参数每n个样点）。对参数

每n个样点用大于1的指定的值，来进行时间损耗分析计算，它不是用于每个时间值，而是

仅用于每n个时间样点选一个。振幅谱暂时加到选定的道数中（滤波参数第一道firsttrace

和最后道lasttrace）。因此得到时间-频率相关性，其距离相关性被消除。

滑动谱分析可以得到随时间移动的频率范围。它可以应用于衰减分析或识别与频率相关的反

射或传播。

1.11.5.4频谱白噪声（Spectralwhitening）

该滤波用来补偿频率随时间的衰减。

该滤波的目的是使给定频带的谱变平坦（参数开始频率startfrequency和结束频率end

frequency）。为这一目的，对每一道做窄重叠的带通滤波，随后计算振幅衰减曲线，得到每

个频带的衰减率。将这些衰减曲线应用到每个带通滤波中，并将其相加。为了保存原始道的

振幅特性，将数据乘以该振幅衰减曲线。比例因子（scaling）将滤波道放大到给定的值。

频谱白噪声用来补偿频率随时间的衰减。

1.11.5.5一次导数（tive）

该滤波用来计算每道时间方向的一次导数。其方法是在频域乘以频率。

一次导数可以增强各道的高频部分，因此在某些情况下可增强垂直分辨率。

1.11.5.6积分（Intergration）

该滤波是对各道在时间方向积分。其方法是在频域除以频率。

积分可以增强各信号的低频部分。

1.11.6插道/拆分(Traceinterpolation/resorting)

1.11.6.1插入标记（Markerinterpol）

根据手工设置的标记或自动从数据中抽取的标记（子选项读标记readmarker），用该选项可

以在X方向插入数据。

激活该选项后，出现一列表，可以用道（道数tracenr.）和相邻标记的距离(标记增量markerinc.)

来输入标记位置,根据希望的道增量(traceincrement)和实际设置的两个标记间的距离,可以用

以下公式重新计算道数:计算的道数=标记增量/道增量。标记增量（markerincrement）必须

在滤波参数组中予先指定，并可以在输入的表中改变每一个标记位置。

道增量必须是每个标记增量的倍数，该程序自动计算整个标记范围，然后根据实际道增量确

定需要的道数，如果该数与实际道数不一致，将会出现警告信息。

18

标记位置可以用鼠标左键在原始剖面中指定，也可在输入列表中输入道数。各个标记增量必

须在列表中输入。另外可以从原始文件的道头上读标记位置（选项读标记ReadMarker），也

可以用选项注释标记（commentmarker）激活注释标记。如果该选项用于批处理中，标记位

置必须从道头文件读出。激活选项注释标记，就能用距离和注释标记；否则只能用距离标记。

1.11.6.2道增量－样点重组（Traceincr-resampling）该选项是用指定的新道增量在X方向对

数据进行重组。

如果由新道增量产生的新道数小于原始文件的道数，一些道将会略去；如果由新道增量产

生的新道数大于原始文件的道数，将会增加一些道。

如果将不同道增量的数据和起来，该选项会很有用。

1.11.6.3创建等距离道（）该选项可以插入非等距离数据，使其变成等距

离数据。

1.11.6.4Y方向翻转剖面（YFlipProfile）

该选项在时间（y）方向翻转剖面。

1.11.6.5X方向翻转剖面（XFlipProfile）

该选项可以在X方向翻转剖面。

1.11.6.6求助道头（Resort-traceheader）

该选项根据道头文件中保存的距离，来在X方向求助数据。

1.11.6.7求助组（Resort-group）

该软件对给定的道组来求助数据。

1.11.6.8插道-3D文件（Traceinterpol-3DFile）

该选项在X方向重新将2D剖面插入3D文件中。

1.11.6.9X翻转－3D文件（XFlip-3DFile）

该选项在X方向翻转各个3D文件。

1.11.7编辑道/道范围

1.11.7.1移去（Remove）

该选项可以删除各道或道范围。

1.11.7.2移去范围（Removerange）

该选项在剖面的开始和结束处删除若干道。

1.11.7.3插道（Interpolate）

该软件可以在相邻两道之间插入新道。当选择“marginaltraces”时，它将复制相邻道。

1.11.7.4抽道（Extract）

该选项从原始文件中分别抽取各道或道范围。

19

1.11.7.4替换（Replace）

该选项用相等的其它文件的相关道来替换某一道或道范围。

1.11.7.6反极性（Reverpolarity）

该选项可以对某一道或道范围反极性。

1.11.7.7设置为零（Settozero）

该选项将某一道或道范围设置为零。

1.11.7.8复制（Duplicate）

该选项复制某一道或道范围。

1.11.7.9移动（Move）

该项可以在时间方向移动各道（为了修正可能的触发错误）。

1.11.7.10插入剖面（Inrtprofile）

该选项可以在特殊的位置插入剖面。

1.11.7.11混合剖面（Mixprofile）

该选项可以在原来剖面上加、减、乘或除第二个剖面。

1.11.7.12联合文件（）

该选项将原始文件和其它几个文件合并成一个新文件。

1.11.7.13合并文件（mergefiles）

该选项将原始文件和其它几个文件合并成一个新文件。

1.11.7.14插入0道（inrtzerotraces）

该选项可以插入零振幅值的道。

1.11.82D滤波（2D-Filter）

1.11.8.1滑动平均(Runningaverage)

该滤波作用于选定的道数，它对每个时间段中选定的道数进行滑动平均。

滑动平均作用于选定的道（参数平均道averagetraces），最大道数为256。如平均道为4

时，水平方向左边两个点，右边两个点，即对每一个时间值都考虑5个点。

可以限制滤波作用的时间范围，此时必须输入参数开始时间（starttime）和结束时间（end

time）。数据的其它部分不进行滤波。作为缺省值每道开始时间设为数据的起点，结束时间

设为每道最大时间处，

该滤波的作用是抑制噪声，它的效果表现在水平方向的能量一致性。

1.11.8.2抽取平均（Subtractingaverage）

该滤波作用于选定的道数，它对每个时间段中选定的道数进行抽取平均。滤波执行所谓的

滑动背景去除。

20

抽取平均作用于选定的道（参数平均道averagetraces），最大道数为256。如平均道为4

时，水平方向左边两个点，右边两个点，即对每一个时间值都考虑5个点。

可以限制滤波作用的时间范围，此时必须输入参数开始时间（starttime）和结束时间（end

time）。数据的其它部分不进行滤波。作为缺省值每道开始时间设为数据的起点，结束时间

设为每道最大时间处，

该滤波的作用是抑制水平一致能量，它的效果表现在有差异的信号部分。

1.11.8.3背景去除（Backgroundremoval）

该滤波作用于选定的道数，它在实际剖面中选定的时间/距离范围内抽取一个平均道（道

范围）。它执行所谓的背景去除。

必须定义平均道的时间/距离范围，该平均道显示在参考道窗口中。时间范围由滤波参数

开始时间（starttime）和结束时间（endtime）来定义；距离范围由滤波参数开始距离（start

distance）和结束距离(enddistance)来定义。也可以在原始剖面中交互式定义数据范围。

用按钮可以指定参考道是由整个剖面（wholeline）决定，或指定的距离范围（distancerange）

决定，或仅用距离范围（nly）,或由剖面部分（lineparts）自动决定。当使用整个

剖面或剖面部分时，滤波参数开始距离和结束距离不再重要。

采用该选项可以从整个剖面中消除一致性的噪声，使有用信号更加清楚。该滤波也可以抑

制水平一致的能量。其效果也表现在有差异的信号部分。

注意：该滤波也可能产生假信号。

1.11.8.4平均xy滤波（averagexy-filter）

该滤波对每个时间段的选定的xy区域计算平均。

该平均可以在许多道之间进行（参数道数es(x)），也可以在许多样点之间进行（参

数样点数les(y)）。当选择较大区域时，该滤波效果较好。

也可以限制滤波的时间范围，此时必须输入开始时间（starttime）和结束时间(endtime)。

该滤波可以抑制道和时间的噪声，其效果相当于在x方向和y方向的低通滤波。

1.11.8.5中值xy滤波（medianxy-filter）

该滤波对每个时间段的选定的xy区域计算中值。

该中值可以在许多道之间进行（参数道数es(x)），也可以在许多样点之间进行（参

数样点数les(y)）。当选择较大区域时，该滤波效果较好。

也可以限制滤波的时间范围，此时必须输入开始时间（starttime）和结束时间(endtime)。

该滤波可以抑制道和时间的噪声和毛刺，其效果相当于在x方向和y方向的低通滤波。

1.11.8.6压缩（Compress）

该选项可以在时间（Y）方向和距离（X）方向进行压缩。

滤波参数factor(x-dir)和factor(y-dir)决定在x方向和y方向压缩窗口的尺寸。

当应用计算时间加强处理步骤（如fk滤波）前，该压缩非常有效。

1.11.8.7扩展（Expand）

该选项可以在时间（Y）方向和距离（X）方向进行扩展。

滤波参数factor(x-dir)和factor(y-dir)决定在x方向和y方向扩展窗口的尺寸。

请注意应用该选项分辨能力不会增加。

21

1.11.8.8叠加道（Stacktraces）

该选项可以对选定的道进行瞬时叠加。

滤波参数时输入要叠加的道数（nroftraces）。

该滤波可以增加信噪比。

1.11.8.9抽道（Subtracttraces）

该选项可以抽取每一组道的第一道。

滤波参数es是确定道数组，该滤波将每一组的第一道抽去。

1.11.8.9压缩3D文件（Compress3D-file）

该选项可以在两个方向压缩3D文件，一个3D文件包括许多平行的2D剖面。

当应用计算时间加强处理步骤（如fk滤波）前，该压缩非常有效。

1.11.9偏移（Migration）

1.11.9.1衍射叠加（Diffractionstack）

它只对零偏移距（即收发天线在同一点）数据有效。

1.11.9.2克希霍夫偏移（Kirchhoffmigration）

（待续）

六拾取不同反射体，层显示和

创建报告

下面介绍拾取不同反射体，定义深度转换速度，将不同的反射体合并到所谓的层显示图像。

一、一般流程

1.拾取反射体-反射体不一定是连续的,它也可以由不同的孤立体组成。

2.用特殊的名字（如layer1）保存拾取的一个反射体-输入层数（layernumber），平均

速度（meanvelocity）和层编码（layercode）。

3.重复步骤1和2，直到拾取所有的反射体为止。

4.进入层显示（layer-show）

5.如果用层速度代替平均速度，则创建速度文件

6.创建层显示

7.创建报告

二、拾取反射体

每个反射体必须分别拾取，反射体不一定是连续的,它也可以由不同的孤立体组成。下面描

述不同步骤：

1.进入模块2D-dataannalysis

2.载入(load)想要的GPR或地震ZO(0偏移距)剖面(GPR剖面一般是0偏移距的)

3.点击选项拾取(pick)

22

4.用下述拾取功能之一拾取反射体(pickthereflector)

相位追踪(phafollower):该选项用于长剖面（建议用显示选项PixelsPerTrace）而且反

射体非常连续，用鼠标右键点击停止自动相位追踪，然后可以在不连续处强制进行相位

追踪，此时已存在的追踪线将用新追踪线来代替。

连续拾取（continuouspick）：用该选项进行连续手动拾取。

手动拾取（manualpick）：用该选项进行单道拾取或改变各个拾取点。

自动拾取（autopick）：用于连续反射体。

5.用选项保存(save)用REFLEXW格式来保存一个反射体的拾取,输入下列参数:

编码(globalcode):任何编码

层数(layernumber):用于合并不同拾取反射体,最小层数是0,最大层数是100,每个反射体

有自己的层数。

平均速度（meanvelocity）：输入平均速度用于深度转换。

文件名（filename）：用于拾取的文件名（扩展名自动定为pck）

6.重新设置(ret):当你要拾取新的反射体时,重新所有的拾取

7.重复步骤4-6知道所有反射体拾取完为止,注意：对每个不同反射体必须改变层数，

最大层数是100。

8.点击选项层显示（layer-show）或在analy/layer-show下（必须载入相应的剖面），

出现新的层显示参数菜单。

三、产生用于时深转换的速度文件（打开层显示菜单）

用选项createLayerShow做时深转换时（time-depthconversion），必须考虑速度。

一种是对每个拾取文件的整个剖面用常数平均速度（constantmeanvelocities），它是用拾取

处理来设置的。

另一种是用ASCII速度数据文件（ASCIIvelocitydatafile），此时每层速度也许会非常倾斜，

并将其作为层速度。该ASCII文件包括以下信息：

注释行commentline

注释行commentline

第一速度行tyline

……..

第n速度行tyline

每个速度行包括:

distancev-layer1v-layer2……v-layer100

23

层速度可以不定义在希望的距离位置,单程序在两个给定值间线性插入。

该ASCII速度数据文件可用任何编辑器手动产生，也可用核文件（corefiles）或在layerShow

中用选项产生速度（createvelocity）。最后选项详述如下：

1.在layershow面板内点击产生速度（createvelocity），打开createvelocity菜单。

2.有两种不同方式产生ASCII速度数据文件：层速度手动设置，它意味着一层内速度没有

倾斜变化；或速度用拾取层边界的振幅校准（amplitudecalibration）来计算。

3.对个层输入速度（velocity）和文件名（filename）（速度文件的扩展名自动取为vla）。

4.点击选项开始（start）产生ASCII文件，该文件显示在窗口中；点击选项关闭（clo）

关闭文件显示。

5.点击选项关闭（clo）来关闭产生速度菜单。

文件可以用任何编辑器手动改变，以考虑速度的变化。

四、将不同拾取的反射体合并到一层模型中（打开层显示菜单）

该选项将各拾取文件合并保存到一个文件中，用曲线文件显示，并将报告文件输出到文件或

打印机。每层最大合并拾取数等于每个剖面最大道数，最大层数是100。

反射体的拾取传输时间用平均速度或速度文件转换成深度。

以下是产生层显示的详细步骤：

1.在层显示菜单内击活选项产生(create),打开产生层显示菜单(见右图)

2.在层显示菜单内点击选项拾取文件(pickfiles),选择想要的拾取文件（用shift键做多个选

择），选择的文件列在右边的窗口内。

3.选择其它参数如开始/结束位置（start/endposition），增量（increment）和时深转换的速

24

度选项（velocitychoice）。

如果不用平均速度，而是用速度文件（velocityfile），首先必须创立速度文件，然后点击

选项velocityfile载入想要的速度文件，选项插入层（interpolatelayer）控制如何处理孤

立的反射体。

4.输入用于层显示的文件名（filename，扩展名是lay），点击start创立层显示。

5.点击clo关闭创立层显示菜单。

6.层显示显示在下面的窗口，传输时间拾取与剖面一起显示在上面的窗口。

7.用层显示菜单的load选项栽入任何显示。

注意：如果平均拾取速度（meanpickvelocities）用于时间深度转换，每层的时深转换与上

层转换无关。与之对比，如果用ASCII速度数据文件（ASCIIvelocitydatafile），每层的时深

转换与上层转换有关。因此，如果上面某层的拾取不连续，必须控制特殊层的时深转换，此

时应用选项插入层（interpolatelayer）。

如果点击该选项，则假定上面每层都连续。

如果不点击该选项，则用紧邻的上层速度。它将导致在上层不连续处的实际层有尖锐台阶。

五、产生层显示的输出报告（打开层显示菜单）

当前层显示的结果可以保存在报告中，可以用打印机打印，也可以保存成文件。根据各设置，

报告包括单层信息如厚度，速度，振幅和坐标。另外，也可以报告出实际载入的核文件和标

记注释。它也可以减去单点信息如代表两层结构的厚度。

1.在层显示菜单内点击选项产生报告(createreport),打开产生报告菜单。前提是必须载

入层显示或创立层显示。

2.输入必要的输出参数：

-开始/结束位置（start/endposition）和增量（increment）

-层（layer）：输入要报告的层数：如1-5或1-5，8，9

3．点击选项开始（start）。

（3D待续）