promax处理流程许辉群课件

1、实习要求1、按照计划时间完成实习任务；2、按照每项实习计划考核，每项9分，共计99分，考勤计分为：1分。其中：每个实习项目9分包括：a、软件操作+报告一、实习内容(1)定义观测系统使CMP道集分选正确(5)静校正检查静校量空间分布的合理性。(2)叠前去噪剔除不正常的炮、道和野值；压制强面波和其他不规则干 扰，使信噪比较高，波组特征清楚。(3)振幅补偿做相对保持振幅处理,使浅、中、深、道、炮间能量均衡(4)反褶积用好多道反褶积，实现零相位输出，达到压缩子波、提高地 震记录分辨率。(6)速度分析用速度谱或速度扫描获取速度,每1km不少于2个速度谱。(7)动校正和叠加叠加速度合理切除参数正确最终叠加

2、剖面质量较优 (8)倾角时差校正断面波、绕射波、水平及倾斜反射成像清楚。(9)叠后时间偏移偏移速度有效波归位合理,绕射波收敛,断点清楚。二、UNIX系统常用命令列出文件目录中内容改变文件目录拷贝文件创建目录删除文件及目录显示进程状态删除进程退出查阅命令三、启动ProMax系统1、启动仿真软件Exceed2、选择主机lade1 双击鼠标左键、输入用户名和口令用户名: pmxusr0116口令: 123456、设置环境变量setenv DISPLAY 192.168.1.21:0、启动ProMAX promax 或 promax3d 6、建立工区、测线和流程目录 ，即： Area 工区菜单Sele

3、ct 选择一个存在的工区； Add 加一新工区,即创建工区目录；Delete 删工区； Rename 工区更名；Copy 复制工区；Permission 改变工区的许可权 Line 测线菜单Select 选择一个存在的测线； Add 加一新测线,即创建测线目录；Delete 删测线； Rename 测线更名；Copy 复制测线；Permission 改变测线的许可权用来存放与地震测线有关的信息，如地震数据、几何库定义文件、参数表（速度、初至时间）、和有序参数文件( Ordered Parameter Files) Flow 流程菜单Add 加一流程，即创建流程目录，该目录 中有如下四个文件：

4、DescName:存放流程描述名 TypeName: 存放产品名(p、v、3) packet.job: 存处理列表和参数 job.output: 存放流程执行产生列表用以下作业核实/export/blade1/disk1/userdata/dty/dg_sqzh/1703文件中的FFIDDatasets数据集，击MB2可显示信息：数据 集的大小(字节)、数据类型、总道数、样点数、采样率、最小和最大 CDP数、炮点索引号等。击MB3可 显数据集的处理历史。Database数据库7、选择模块组织作业MB1MB2MB3四、处理流程1、定义观测系统 运行如下模块: 3D Land Geometry S

5、preadsheet* 弹出如下菜单:按从左到右顺序执行菜单项指定计算中点的方法:按地震道中现有的索引号来映射按震源表和班报表中的模式号匹配按震源表和班报表中的线号和桩号来匹配站号间隔:(标名)接收点站号间隔震源站号间隔横测线距离或间隔测区或工区方位角震源和接收点站号共线吗?震源类型：炮井 地面震源单位： 米 英尺坐标原点所有的x坐标减去X0(该值)所有的y坐标减去Y0(该值)例1 3D Land Geometry Spreadsheets Setup （Geometry Setup） Matching line and station numbers in the SIN and PAT s

6、preadsheets按炮点和模式表中的线号和站号来匹配。 Nominal Receiver Station Interval: 接收点间隔 50 Nominal Source Station Interval: 震源点间隔 100 Nominal Crossline Separation: 横测线间隔 25 Nominal Survey Azimuth: 工区方位角 335.67 Base Source Station co-ordinates upon a match between source and receiver station numbers? 炮线和接收线共线吗？ Yes N

7、o Surface Seismic source Units : Meter Feet Receivers (填写接收点参数表） 由此表生成接收点有序参数文件SRF Ordered Parameter FileStation接收点站号(桩号)本例为498825/线Pt index站点识别器(19)全 1x接收点x坐标为测量数据y接收点y坐标为测量数据Elev接收点高程为测量数据Line接收点线号本例线号为19Static接收点静校量全 01。击MB12。击MB1import3。选择文件路径4。击MB15。(1)用MB1选始行 (2)用MB2选末行 (3)Ctrel+d 去掉行6。击MB1，弹出

8、7。定义文件名，再击ok8。击station9。按MB1 拖动10。用MB1击,则将数据输入表中。 Sources (填写炮点参数表）Source用户定义的震源编号线号*10000+站号Line震源线号1060Station震源站号1153St Index同一炮点识别器(19) 全1 (每个站点只放一炮)x震源点的X坐标 y震源点的Y坐标z震源点的高程FFID野外文件号按班报填写Time放炮时间可以不填Date放炮日期可以不填Offset炮点垂直炮线的偏移距,+右、-左可以不填Skid炮点平行炮线的偏移距,+大、-小号可以不填Shot fold*接收道数变化的Pat Shift排列滚动的站点值

9、全0Static用户定义的静校量全0注意：有意将表中坐标与文件中坐标填写不一致 Patterns (填写排列或班报参数表）Src Line炮线号1060Src Point炮站号1-153Src Index炮索引号与St Index一致全1Min Chan最小道号Max/Gap Chan最大道号Chan Inc道号增量全1Rcvr Line接收线号19Rcvr MinChan最小道号所在的桩号Rcvr MaxChan最大道号所在的桩号Rcvr Index接收点索引号与Pt Index一致全1 Bin (数据面元化) 必须执行下面1、3、5项1、指定中点：按SIN和PAT表中的线号和站号来匹配2、

10、定义面元网格3、面元中点4、面元数据质量控制5、数据库定案计算炮点到检波点间的距离；计算炮点到检波点的中点坐标；计算炮点到检波点的方位角CDP=炮点+检波点本例网格原点位置（min inline & crosslineXorg=532051.50Yorg=4282645.50网格的方位角面元边长dx面元边长dy网格的坐标原点Inline 间的距离crossline 间的距离存放本表的文件名Cdp、inline、crossline号形成CMP 道的最小/最大炮检距炮检距面元之间距离震源-接收点间方位角方位角公差当定义面元化参数后，计算网格维数沿网格Y轴的方位角面元在x方向的边长，即dx面元在Y方

11、向的边长，即dy3D网格的坐标原点网格在x方向的距离(/dx= x方向上的线数)网格在Y方向的距离(/dy= y方向上的线数)存放网格参数的文件名起始CDP号起始纵线号起始横线号距面元中心的最小炮检距距面元中心的最大炮检距面元化时炮检距增量震源接收点的方位角震源接收点的方位角公差纵线平行于网格的x轴纵线平行于网格的y轴Dx=50mDy=25mBin=面元击Database,弹出TRC记录道SRF地面位置SIN炮点索引CDP共深度点CHN道OFB偏移距面元ILN纵测线XLN横测线PAT排列模式浏览数据库属性：get1、检查覆盖次数的变化情况：CDPFOLD2、检查每炮的道数变化情况：SINNCH

12、ANS3、检查每炮的有效检波点：View-Predefined-SIN-SRF-offset4、检查CDP上的炮检距分布：View-Predefined-offset-CDP-SIN5 、检查SRF上的高程分布：SRF-Geometry ELEV6 、检查SIN上的高程分布：SIN-Geometry ELEV覆盖次数按Inline-xline显示的覆盖次数和CDP号：2、将几何信息装入道头、地震数据写入库内 SEG-Y Input IF （选择FFID） Inline Geometry Header Load Disk Data Output ( 注意 New 和 Append的切换) END

13、IF40张卡片映像(共3200字节)卷头信息(共400字节)第一道的道头信息(共240字节)第一道的地震数据(常为IBM浮点格式)第二道的道头信息(共240字节)第二道的地震数据(常为IBM浮点格式)1任务号（4字节)2测线号（4字节)3卷号（4字节)4道数/炮或总道数5采样间隔采样间隔us6样点数/道7数据格式码：1浮点（4字节) 2定点（4字节） 3定点（2字节） 8记录类型码：0叠后道 1炮集 2CDP道集 913备用14计量系统代码：1米 2英尺 备用100备用SEG_Y 格式：卷头(4字节/字,共100字）:1测线内道序号 (字节1 4)2卷内道序号 (字节5 8)3FFID & I

14、LN (字节9 12)4道号 & XLN (字节1316)5震源点号 (字节1720)6CMP号 & XLN (字节2124)7CMP集内道号 (字节2528)8道识别码： (字节2930) 1=地震数据；2=死道；3=空道 4 =爆炸信号；5 井口道；垂直叠加道数： (字节3132)9水平叠加道数 ： (字节3334) 数据类型码： (字节3536) 1=生产；2=试验64 (字节237240) 道头(4字节/字,共60字）:SEG-Y 类型使用的存储类别输入磁盘文件名当采样间隔不正确时,可忽略输入数据的采样间隔,给出其正确值.每道的样点数(在忽略2#道头时给出)处理历史存入卷头吗?输入辅助

15、道吗?从道头中得到道号吗?输入道格式显示道集信息吗输入的最大时间是叠加数据吗?每集的最大道数主分选道头字输入第一选择吗?输入第二选择吗输入全局的XY参考坐标吗?若为yes,则道头中xy减此值.使用坐标比例重置道头吗?选道: 1-10(2)“(2)”表示隔一道读一道.道选择方式:Include包含Exclude拒绝列出指定的道连机几何信息道头装入与数据库对应的主道头字用有效道号来匹配吗？去掉道头中为空CDP的道吗？去掉道头中为空接收点的道吗？冗余诊断吗？道头运算选择方式定义道头运算表达式固定等式方式顺序重编号方式旋转与平移方式改变数据属性磁盘数据输出指定输出数据集的文件名是新文件还是已存在的文件

16、？输出的记录长度指定样值存储格式跳过主磁盘存储区吗？8、16位 为整型数，每一百个样值有一 比例因子；32位 为IEEE格式的浮点数。3、叠前参数测试(1)获取初至切除、去废道废炮参数 Disk Data Input Automatic Gain Control Trace Display磁盘数据输入从其他工区、测线中读取数据吗？选择数据集繁衍输入文件的处理历史吗？读道选件: get all/sort交互数据存取吗？选择进入的主、次和第三道头字对数据库分选顺序表在内存或磁盘分检多次读数据吗？仅处理道头吗？忽略输入数据的采样间隔吗？分选顺序表填写举例：( 所用符号：“()”,” ”,”,”/ ”

17、, ” * ”, ”:”, “,”,)120(4)表示取 1，5，9，13，17; “(4)”为增量; “”表示范围。3(4) 表示取3,3,3,3; “(4)”表示前一值的重复次数。 1(2), 4-10(2) 表示取1，1，4，6，8，1043 表示取顺序值,即4,5,6; 又如130(10)3,即1,2,3,11,12,13,21,22,23 * 通配符，其范围-999999.0999999.0; / 分组符如10-14(2)/15-21(3)表示取10,12,14,15,18,21; : 分级符如：1100：116 (前为主头字，后次道头字)应用方式: 见增益因子的类型, 见增益因子长

18、度 (ms)因子的应用基准, 见拒绝硬化零吗?鲁棒比例吗?(强行求因子)只应用增益因子,不保存于道头中应用增益因子,并将因子保存于道头中去掉增益控制,并删除道头中增益因子输出每道的增益因子自动增益控制应用基准绝对平均值均方根值中值(按滑动时窗求上面三值)中心头部尾部道显示选择显示设备显示的起止时间(ms)每屏的最大道数每屏间的重叠道数每屏的道集(线段)数使用变道间隔吗？输出方式道显示方式变面积基线道头显示参数自动保存屏幕吗？保存屏幕图像的最大数按彩色方式保存屏幕吗？保存屏幕图像到何处收集的屏幕数道显示方向/极性标注的主/次道头字次道注释方式/及增量道比例方式(习惯的)剪切的道偏离值样值所乘的标

19、量值道比例选件(单/多道)道显示方向(垂直/水平)SOURCE:有效炮号(用户定义的)FFID:野外文件标识号SIN: 炮号(内部号)REC_SLOC: 检波点号SOU_SLOC: 外部震源位置号SRF_SLOC: 外部检波点位置号CDP: CDP面元号OFFSET: 带符号的炮检距AOFFSET: 绝对炮检距CHAN: 记录道号REPEAT: 重复数据的复制号CS_STAT: 炮点相关_自动统计静校正量CR_STAT: 检波点相关_自动统计静校正量PS_STAT: 炮点最大能量_自动统计静校正量PR_STAT: 检波点最大能量_自动统计静校正量TRIMSTAT: 修饰静校正量FB_PICK:

20、 初至拾取时间TR_FOLD: 实际道的覆盖次数DMOOFF: DMO偏移距面元ILINE_NO: 纵测线号XLINE_NO: 横测线号CDP_SLOC: 外部CDP位置号删除道极性反转拾取顶部切除拾取底部切除拾取手术切除拾取各种各样的时窗编辑道头值编辑数据库值拾取其他层位道切除重新应用以前的切除吗?切除时间参考切除类型斜坡时间(ms)外推切除时间从数据库得到切除文件吗?选择切除参数文件(1) 参数表拾取 disk Data Input Read datafrom other lines/surveys?.No Select datasetShots-with-geometry trace r

21、ead option.sort interaction Data Access?.yes select primary trace header entry.FFID select secondary trace header entry.OFFSET select tertiary trace header entryNONE sort order for dataset.*:*/Automatic Gain Control AGC operator length1000Trace Display Number of ENSEMBLES(line segments)/screen2 Do y

22、ou want to use variable trace spacing.yes运行此流程拾取如下参数表：Kill tracePick Top mutePick Miscellaneous Time GatesPick Reverse trace(2)真振幅恢复参数测试 DISK Data Input Interactive Data Access?.No Select primary trace header entry.SIN Sort order list for dataset.16/ Trace Kill/Reverse Trace editing MODEkill Get edi

23、ts form the DATABASE?.yes Trace Kill/Reverse Trace editing MODEReverse Get edits form the DATABASE?.yes Trace Muting Parameter Test Enter parameter VALUES:.12|9|6 Trace grouping to reproduce.Ensembles True Amplitude Recovery Apply spherical divergence corrections?.yes Basis for spherical spreading1/

24、dist Apply inelastic attenuation corrections?.No Get TAR velocity function from database?.No Should the velocity be treated as space variableNo Specify TAR velocity function: 0-1500,500-2000, Apply dB/sec correction?.Yes dB/sec correction constant:.99999 Apply time raised to a power correction?.No A

25、pply function to data or Remove effect of amplitude correction?.Apply Maximum appliction TIME:5500Trace Display Number of Ensembles(line segments)/screen:.2 Trace scaling option:Entire Screen Number of display panels:.2说明：球面扩散(spherical divergence)地震波的振幅随传播距离的增大而减小的现象，称为波前或球面扩散。非弹性吸收(inelastic atten

26、uation)由于岩石是非完全弹性体，地震波在传播过程中一部分能量将转化为热能而消耗掉，从而使波的能量变弱的现象，称为地层的吸收衰减作用。真振幅恢复(True Amplitude Recovery)就是从地面检波器记录的原始振幅中，消除波前扩散和地层吸收衰减的影响，使其恢复到仅与地下反射界面的反射系数有关的真振幅值，使恢复后的地震波振幅正比于反射系数的大小。这种保持反射波振幅与反射系数呈正比的剖面称为相对振幅保持剖面。(4)反褶积参数测试Disk Data InputTrace kill/reverseTrace mutingTrue Amplitude Recovery dB/sec cor

27、rection constant .9Reproduce Traces Trace grouping to reproduce.Ensembles Total Number of datasets3IF Trace selection MODE .Include Select primary trace header wordREPEAT Select secondary trace header word.NONE Specify trace list1 Trace Display Label Trace labelOriginal InputELSEIF Trace select MODE

28、.Include Select Primary trrace header word.REPEAT Selectr secondary trace header wordNONE Specify trace list.2 Spiking/Predictive Decon Trace Display LabelELSE Spiking/Predictive Decon Trace Display LabelENDIFTrace Display 处参数说明： Type of deconvolutionMinimum phase spiking Decon operator length(ms).1

29、60 Operator white noise level.0.1 Window rejection factor2 Time gate reference.time Get decon gates from the DATABASE?.yes Select decon gate parameter filedecon gate Output traces or filersNormal decon output Apply a bandpass filter after decon?.No Re-apply trace mute after decon?.yesTrace Display L

30、abel Trace label:.Decon 处参数说明：Spiking/Predictive Decon Apply a bandpass filter after decon?.yes Bandpass filter frequency values:6-12-70-80 其余参数见Trace Display Label Trace Label Decon with filterENDIFTrace Display Number of Ensembles/screen.3 Number of display panels.1 trace scaling option.Individual

31、(3)F-K分析和滤波F-K Analysis Panel width in traces(面板宽度-道数) .136 Starting time for analysis(分析的起始时间) .0 Ending time for analysis (分析的终止时间) .0 Distance between input trace (输入的道间距) .50 Starting display configuration (初始显示设置) .TX-TK-FX-FK Position of zero wavenumber in display(零波数位置) .Center Position of ze

32、ro frequency in display (零波数位置) .Top Plot FK,TK,TX panels in DB or Linear?(以分贝/线性绘) .DBSCALE Initial TX gain setting(percentile) (设置TX初始增益) .98 Initial FK maximum gain setting(db) (设置FK最大增益) 0 Initial FK minimum gain setting(db) (设置FK最小增益) .0 percent flat for trace ramping (道斜坡平坦百分比) .100 (不斜坡) perc

33、ent flat for time ramping (时间斜坡平坦百分比) 100 (不斜坡) Select mute polygon table (选择切除多边形表) .fk mute Mode of F-K filter operation (F-K滤波器操作方式) .REJECT Percent flat for F-K filter windowing (FK滤波器窗口平坦百分比) .90 Time length of F-K filter(ms) (FK滤波器时间长度)500 Spatial extent of F-K filter(traces) (FK滤波器空间范围）.50F-K

34、滤波Distance between input trace (输入道的间距).50Panel width in trace (以道计的面板宽度)136Test the filter impulse response?(测试滤波器的脉冲响应吗)NoPercent flat for time ramping (时间斜坡平坦百分比) .100Percent flat for offset ramping (偏移斜坡平坦百分比).100Get polygon mute file from the database(从数据库中得到多边形切除文件)yes select mute parameter fi

35、le (选择切除参数文件).fk_muteMode of F-K filter operation.RejectPercent flat for F-K filter windowing.90Time length of F-K filter(ms).500Spatial extent of F-K filter(traces)50Re-apply T-X trace mute after filter?.(滤波后重新道切除吗).yesPercentage of K-space to keep around k=0(在K=0附近保留的k空间百分比)0 Trace Display Label T

36、race labelFK FilterELSE Trace Display LabelENDIFTrace Display Number of Ensembles(line segments)/Screen.2(5)频谱分析确定有效波主频；检查反褶积效果。Disk Data Input Select dataset.Shots with geometry Trace Read option.Sort Interactive Data Access?.No Select Primary trace header entry.SIN Select secondary trace header en

37、try.NONE Sort order for dataset.16/Automatic Gain ControlInteractive Spectral Analysis Data select methodSimple Display data by traces or ensembles.Ensembles检查反褶积效果Disk Data InputAutomatic Gain ControlReproduce Traces Trace grouping to reproduceEnsembles Total number of datasets.2IF Select Primary t

38、race header wordRepeat Specify trace list.1ELSEIF Select Primary trace header word.Repeat Specify trace list.2 Trace Muting Select mute parameter file.FB mute Spiking/Predictive Decon Decon operator length(ms)160 Select decon gate parameter file.decon gateENDIFInteractive Spectral Analysis Data sele

39、ct methodMuliple Subset(多子集方式) Freeze the selected subset?(冻结所选子集).yes Display data by traces or ensemblesEnsembles 4、高程静校正四种基准面的定义统一基准面：是一个由用户提供的与海平面平行的水平面，最终叠加和偏移结果必须校正到统一基准面。CMP参考面：是一个以CMP为单元统计出的等时面，其静校正量为该CMP道集中的炮点静校正量和接收点静校正量的算术平均值。速度分析、动校正和叠加在此面上进行。辅助基准面：是由用户临时定义的接近地表的一个水平面，一般在偏移处理时使用。浮动基准面：对C

40、MP高程经过平滑后所得到一个平滑面(非水平面)。(2)约定校正量定义：校正量应用：当校正量为正值时，记录往上提；当校正量为负值时，记录往下拉。CMP参考面校正量(3)计算高程静校正量 Datum Statics Calculation* Elevation or Refraction(高程或折射).Elevation Final datum elevation(最终基准面高程)-10 Replacement velocity(替换速度)1600 Database math method.Shot Holes Using Uphole Info NMO static method .Elevat

41、ion Length of smoother.51 Processing DATUM.NMO DATUM Run ID.01 用DBTools 选择SRF，观察下列属性：R_ATAT01 、 F_DATUM DATUMVEL 、ELEV，注意静校正量与高程呈反向关系。用DBTools 选择CDP ，观察ELEV和N-DATUM替换速度=基岩速度=折射波时间线斜率倒数；风化层速度=直达波时间线斜率倒数。(4)应用高程静校正Disk Data Input Select dataset.Shot-with geometry Trace read optionGet allTrace kill/reverseTrace MutingTrue Amplitude