《promax处理流程-许辉群.ppt》由会员分享,可在线阅读,更多相关《promax处理流程-许辉群.ppt(102页珍藏版)》请在三一文库上搜索。
1、地震资料数字处理,实习,学时:3 周,指导教师:许辉群,实习要求,1、按照计划时间完成实习任务; 2、按照每项实习计划考核,每项9分,共计99分, 考勤计分为:1分。 其中:每个实习项目9分包括:a、软件操作+报告,一、实习内容,(1)定义观测系统使CMP道集分选正确 (5)静校正检查静校量空间分布的合理性。 (2)叠前去噪剔除不正常的炮、道和野值;压制强面波和其他不规则干 扰,使信噪比较高,波组特征清楚。 (3)振幅补偿做相对保持振幅处理,使浅、中、深、道、炮间能量均衡 (4)反褶积用好多道反褶积,实现零相位输出,达到压缩子波、提高地 震记录分辨率。 (6)速度分析用速度谱或速度扫描获取速度
2、,每1km不少于2个速度谱。 (7)动校正和叠加叠加速度合理切除参数正确最终叠加剖面质量较优 (8)倾角时差校正断面波、绕射波、水平及倾斜反射成像清楚。 (9)叠后时间偏移偏移速度有效波归位合理,绕射波收敛,断点清楚。,二、UNIX系统常用命令,三、启动ProMax系统 1、启动仿真软件Exceed,2、选择主机lade1 双击鼠标左键,、输入用户名和口令,用户名: pmxusr0116,口令: 123456,、设置环境变量 setenv DISPLAY 192.168.1.21:0 、启动ProMAX promax 或 promax3d 6、建立工区、测线和流程目录 ,即:, Area 工区
3、菜单 Select 选择一个存在的工区; 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文件中的FFID,Datasets数据集,击MB2可显示信息:数据 集的大小(字节)、数据类型、总道 数、样点数、采样率、最小和最大 CDP数、炮点索引号等。击MB3可 显数据集的处理历史。 Database数据库 7、选择模块组织作业,MB1,MB2,MB3,四、处理流程 1、定义观测
5、系统 运行如下模块: 3D Land Geometry Spreadsheet* 弹出如下菜单:,按从左到右顺序执行菜单项,指定计算中点的方法: 按地震道中现有的索引号来映射 按震源表和班报表中的模式号匹配 按震源表和班报表中的线号和桩号来匹配 站号间隔: (标名)接收点站号间隔 震源站号间隔 横测线距离或间隔 测区或工区方位角 震源和接收点站号共线吗? 震源类型: 炮井 地面震源 单位: 米 英尺 坐标原点 所有的x坐标减去X0(该值) 所有的y坐标减去Y0(该值),例1 3D Land Geometry Spreadsheets Setup (Geometry Setup) Matchin
6、g line and station numbers in the SIN and PAT spreadsheets按炮点和模式表中的线号和站号来匹配。 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
7、 and receiver station numbers? 炮线和接收线共线吗? Yes No Surface Seismic source Units : Meter Feet, Receivers (填写接收点参数表),由此表生成接收点有序参数文件SRF Ordered Parameter File,1。击MB1,2。击MB1 import,3。选择文件路径,4。击MB1,5。(1)用MB1选始行 (2)用MB2选末行 (3)Ctrel+d 去掉行,6。击MB1,弹出,7。定义文件名,再击ok,8。击station,9。按MB1 拖动,10。用MB1击,则将数据输入表中。, Source
8、s (填写炮点参数表),注意: 有意将 表中 坐标与文件中 坐标填写不一致, Patterns (填写排列或班报参数表), Bin (数据面元化) 必须执行下面1、3、5项,1、指定中点:按SIN和PAT表中的线号和站号来匹配 2、定义面元网格 3、面元中点 4、面元数据质量控制 5、数据库定案,计算炮点到检波点间的距离; 计算炮点到检波点的中点坐标; 计算炮点到检波点的方位角,CDP=炮点+检波点,本例网格原点位置(min inline & crossline,Xorg=532051.50 Yorg=4282645.50,网格的方位角 面元边长dx 面元边长dy,网格的坐标原点,Inline
9、 间的距离,crossline 间的距离,存放本表的文件名,Cdp、inline、crossline号,形成CMP 道的最小/最大炮检距,炮检距面元之间距离,震源-接收点间方位角,方位角公差,当定义面元化参数后,计算网格维数,沿网格Y轴的方位角 面元在x方向的边长,即dx 面元在Y方向的边长,即dy 3D网格的坐标原点 网格在x方向的距离(/dx= x方向上的线数) 网格在Y方向的距离(/dy= y方向上的线数) 存放网格参数的文件名 起始CDP号 起始纵线号 起始横线号 距面元中心的最小炮检距 距面元中心的最大炮检距 面元化时炮检距增量 震源接收点的方位角 震源接收点的方位角公差 纵线平行于
10、网格的x轴 纵线平行于网格的y轴,Dx=50m,Dy=25m,Bin=面元,击Database,弹出,TRC记录道 SRF地面位置 SIN炮点索引 CDP共深度点 CHN道 OFB偏移距面元 ILN纵测线 XLN横测线 PAT排列模式,浏览数据库属性:,get,1、检查覆盖次数的变化情况:CDPFOLD 2、检查每炮的道数变化情况:SINNCHANS 3、检查每炮的有效检波点:View-Predefined-SIN-SRF-offset 4、检查CDP上的炮检距分布:View-Predefined-offset-CDP-SIN 5 、检查SRF上的高程分布:SRF-Geometry ELEV
11、6 、检查SIN上的高程分布:SIN-Geometry ELEV,覆盖次数,按Inline-xline显示的覆盖次数和CDP号:,2、将几何信息装入道头、地震数据写入库内 SEG-Y Input IF (选择FFID) Inline Geometry Header Load Disk Data Output ( 注意 New 和 Append的切换) ENDIF,40张卡片映像(共3200字节),卷头信息 (共400字节),第一道的道头信息(共240字节),第一道的地震数据 (常为IBM浮点格式),第二道的道头信息(共240字节),第二道的地震数据 (常为IBM浮点格式),SEG_Y 格式:,
12、卷头(4字节/字,共100字):,道头(4字节/字,共60字):,SEG-Y 类型 使用的存储类别 输入磁盘文件名 当采样间隔不正确时,可忽略输入数据的采样间隔,给出其正确值. 每道的样点数(在忽略2#道头时给出) 处理历史存入卷头吗? 输入辅助道吗? 从道头中得到道号吗? 输入道格式 显示道集信息吗 输入的最大时间 是叠加数据吗? 每集的最大道数 主分选道头字,输入第一选择吗? 输入第二选择吗 输入全局的XY参考坐标吗?若为yes,则道头中xy减此值. 使用坐标比例 重置道头吗?,选道: 1-10(2) “(2)”表示隔一道读一道.,道选择方式: Include包含 Exclude拒绝 列出
13、指定的道,连机几何信息道头装入 与数据库对应的主道头字 用有效道号来匹配吗? 去掉道头中为空CDP的道吗? 去掉道头中为空接收点的道吗? 冗余诊断吗?,道头运算 选择方式 定义道头运算表达式,固定等式方式 顺序重编号方式 旋转与平移方式 改变数据属性,磁盘数据输出 指定输出数据集的文件名 是新文件还是已存在的文件? 输出的记录长度 指定样值存储格式 跳过主磁盘存储区吗?,8、16位 为整型数,每一百个样值有一 比例因子; 32位 为IEEE格式的浮点数。,3、叠前参数测试,(1)获取初至切除、去废道废炮参数 Disk Data Input Automatic Gain Control Trac
14、e Display,磁盘数据输入 从其他工区、测线中读取数据吗? 选择数据集 繁衍输入文件的处理历史吗? 读道选件: get all/sort 交互数据存取吗? 选择进入的主、次和第三道头字 对数据库分选顺序表 在内存或磁盘分检 多次读数据吗? 仅处理道头吗? 忽略输入数据的采样间隔吗?,分选顺序表填写举例:( 所用符号:“()”,” ”,”,”/ ”, ” * ”, ”:”, “,”,) 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,1
15、0 43 表示取顺序值,即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 (前为主头字,后次道头字),应用方式: 见 增益因子的类型, 见 增益因子长度 (ms) 因子的应用基准, 见 拒绝硬化零吗? 鲁棒比例吗?(强行求因子),只应用增益因子,不保存于道头中 应用增益因子,并将因子保存于道头中 去掉增益控制,并删除道头中增益因子 输出每道的增益因子,自动增益控制,应用基准,绝对平
16、均值 均方根值 中值 (按滑动时窗求上面三值),中心 头部 尾部,道显示 选择显示设备 显示的起止时间(ms) 每屏的最大道数 每屏间的重叠道数 每屏的道集(线段)数 使用变道间隔吗? 输出方式 道显示方式 变面积基线 道头显示参数 自动保存屏幕吗? 保存屏幕图像的最大数 按彩色方式保存屏幕吗? 保存屏幕图像到何处 收集的屏幕数 道显示方向/极性 标注的主/次道头字 次道注释方式/及增量 道比例方式(习惯的),剪切的道偏离值 样值所乘的标量值 道比例选件(单/多道) 道显示方向(垂直/水平),SOURCE:有效炮号(用户定义的) FFID:野外文件标识号 SIN: 炮号(内部号) REC_SL
17、OC: 检波点号 SOU_SLOC: 外部震源位置号 SRF_SLOC: 外部检波点位置号 CDP: CDP面元号 OFFSET: 带符号的炮检距 AOFFSET: 绝对炮检距 CHAN: 记录道号 REPEAT: 重复数据的复制号 CS_STAT: 炮点相关_自动统计静校正量 CR_STAT: 检波点相关_自动统计静校正量 PS_STAT: 炮点最大能量_自动统计静校正量 PR_STAT: 检波点最大能量_自动统计静校正量 TRIMSTAT: 修饰静校正量 FB_PICK: 初至拾取时间 TR_FOLD: 实际道的覆盖次数 DMOOFF: DMO偏移距面元 ILINE_NO: 纵测线号 XL
18、INE_NO: 横测线号 CDP_SLOC: 外部CDP位置号,删除道 极性反转 拾取顶部切除 拾取底部切除 拾取手术切除 拾取各种各样的时窗 编辑道头值 编辑数据库值 拾取其他层位,道切除 重新应用以前的切除吗? 切除时间参考 切除类型 斜坡时间(ms) 外推切除时间 从数据库得到切除文件吗? 选择切除参数文件,(1) 参数表拾取 disk Data Input Read datafrom other lines/surveys?.No Select datasetShots-with-geometry trace read optionsort interaction Data Acces
19、s?yes select primary trace header entryFFID select secondary trace header entry.OFFSET select tertiary trace header entryNONE sort order for dataset.*:*/,Automatic Gain Control AGC operator length1000 Trace Display Number of ENSEMBLES(line segments)/screen2 Do you want to use variable trace spacing.
20、yes,运行此流程拾取如下参数表: Kill trace Pick Top mute Pick Miscellaneous Time Gates Pick Reverse trace,(2)真振幅恢复参数测试 DISK Data Input Interactive Data Access?No Select primary trace header entrySIN Sort order list for dataset16/ Trace Kill/Reverse Trace editing MODEkill Get edits form the DATABASE?yes Trace Kill
21、/Reverse Trace editing MODEReverse Get edits form the DATABASE?yes Trace Muting Parameter Test Enter parameter VALUES:.12|9|6 Trace grouping to reproduceEnsembles,True Amplitude Recovery Apply spherical divergence corrections?.yes Basis for spherical spreading1/dist Apply inelastic attenuation corre
22、ctions?.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 Apply function to data or Remove effect
23、of amplitude correction?.Apply Maximum appliction TIME:5500 Trace Display Number of Ensembles(line segments)/screen:2 Trace scaling option:Entire Screen Number of display panels:2,说明: 球面扩散(spherical divergence)地震波的振幅随传播距离的增大而减小的现象,称为波前或球面扩散。,非弹性吸收(inelastic attenuation)由于岩石是非完全弹性体,地震波在传播过程中一部分能量将转化为
24、热能而消耗掉,从而使波的能量变弱的现象,称为地层的吸收衰减作用。,真振幅恢复(True Amplitude Recovery)就是从地面检波器记录的原始振幅中,消除波前扩散和地层吸收衰减的影响,使其恢复到仅与地下反射界面的反射系数有关的真振幅值,使恢复后的地震波振幅正比于反射系数的大小。这种保持反射波振幅与反射系数呈正比的剖面称为相对振幅保持剖面。,(4)反褶积参数测试 Disk Data Input Trace kill/reverse Trace muting True Amplitude Recovery dB/sec correction constant .9 Reproduce T
25、races Trace grouping to reproduce.Ensembles Total Number of datasets3 IF Trace selection MODE .Include Select primary trace header wordREPEAT Select secondary trace header word.NONE Specify trace list1 Trace Display Label Trace labelOriginal Input,ELSEIF Trace select MODE.Include Select Primary trra
26、ce header wordREPEAT Selectr secondary trace header wordNONE Specify trace list2 Spiking/Predictive Decon Trace Display Label ELSE Spiking/Predictive Decon Trace Display Label ENDIF Trace Display,处参数说明: Type of deconvolutionMinimum phase spiking Decon operator length(ms)160 Operator white noise leve
27、l.0.1 Window rejection factor2 Time gate referencetime 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?yes Trace Display Label Trace label:.Decon,处参数说明:
28、 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 filter ENDIF Trace Display Number of Ensembles/screen.3 Number of display panels.1 trace scaling option.Individual,(3)F-K分析和滤波 F-K Analysis
29、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 zero frequency in display (零波数
30、位置) 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 (不斜坡) percent flat for time ramping (时间斜
31、坡平坦百分比) 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滤波器空间范围).50,F-K滤波,Distance between input trace
32、 (输入道的间距)50 Panel width in trace (以道计的面板宽度)136 Test the filter impulse response?(测试滤波器的脉冲响应吗)No Percent flat for time ramping (时间斜坡平坦百分比) .100 Percent flat for offset ramping (偏移斜坡平坦百分比).100 Get polygon mute file from the database(从数据库中得到多边形切除文件)yes select mute parameter file (选择切除参数文件)fk_mute Mode
33、of F-K filter operationReject Percent flat for F-K filter windowing.90 Time length of F-K filter(ms)500 Spatial extent of F-K filter(traces)50 Re-apply T-X trace mute after filter?.(滤波后重新道切除吗)yes Percentage of K-space to keep around k=0(在K=0附近保留的k空间百分比)0,Trace Display Label Trace labelFK Filter ELSE
34、 Trace Display Label ENDIF Trace Display Number of Ensembles(line segments)/Screen2,(5)频谱分析确定有效波主频;检查反褶积效果。,Disk Data Input Select datasetShots with geometry Trace Read optionSort Interactive Data Access?No Select Primary trace header entrySIN Select secondary trace header entry.NONE Sort order for
35、dataset.16/ Automatic Gain Control Interactive Spectral Analysis Data select methodSimple Display data by traces or ensemblesEnsembles,检查反褶积效果,Disk Data Input Automatic Gain Control Reproduce Traces Trace grouping to reproduceEnsembles Total number of datasets2 IF Select Primary trace header wordRep
36、eat Specify trace list.1 ELSEIF Select Primary trace header wordRepeat Specify trace list.2,Trace Muting Select mute parameter file.FB mute Spiking/Predictive Decon Decon operator length(ms)160 Select decon gate parameter filedecon gate ENDIF Interactive Spectral Analysis Data select methodMuliple S
37、ubset(多子集方式) Freeze the selected subset?(冻结所选子集)yes Display data by traces or ensemblesEnsembles,4、高程静校正,四种基准面的定义 统一基准面:是一个由用户提供的与海平面平行的水平面,最终叠加和偏移结果必须校正到统一基准面。 CMP参考面:是一个以CMP为单元统计出的等时面,其静校正量为该CMP道集中的炮点静校正量和接收点静校正量的算术平均值。速度分析、动校正和叠加在此面上进行。 辅助基准面:是由用户临时定义的接近地表的一个水平面,一般在偏移处理时使用。 浮动基准面:对CMP高程经过平滑后所得到一个
38、平滑面(非水平面)。 (2)约定 校正量定义:,校正量应用:当校正量为正值时,记录往上提;当校正量为负值时,记录往下拉。 CMP参考面校正量,(3)计算高程静校正量 Datum Statics Calculation* Elevation or Refraction(高程或折射)Elevation Final datum elevation(最终基准面高程)-10 Replacement velocity(替换速度)1600 Database math methodShot Holes Using Uphole Info NMO static method Elevation Length o
39、f smoother51 Processing DATUMNMO DATUM Run ID01 用DBTools 选择SRF,观察下列属性:R_ATAT01 、 F_DATUM DATUMVEL 、ELEV,注意静校正量与高程呈反向关系。 用DBTools 选择CDP ,观察ELEV和N-DATUM,替换速度=基岩速度=折射波时间线斜率倒数;风化层速度=直达波时间线斜率倒数。 (4)应用高程静校正 Disk Data Input Select datasetShot-with geometry Trace read optionGet all Trace kill/reverse Trace
40、Muting True Amplitude Recovery F-K filter Spiking/Predictive Decon Datum Static Apply Source datum statics database parameter.SIN GEOMETRY S_STAT01 Receiver datum statics database parameter.SRF GEOMETRY R_STAT01 CDP datum statics database parameter.CDP GEOMETRY C_STAT01 Disk Data Output Output Datas
41、et FilenameShot-decon/elev statics,5、速度分析 (1)数据准备 Disk Data Input (按CDP分选) Bandpass Filter (窄带滤波) Automatic Gain Control(l300-500) Disk Data Output (2) 常速扫描 CVS Strip Analysis*,(2) 速度分析 Disk Data Input,Velocity Analysis,选择显示设备,指定存储速度拾取的表,为计算变量置初值吗?,求和成道集的CDP数,是否对面元应用动校正?,分析数据预先计算吗?,是否应用静校正?,相似速度采样率,
42、相似速度计算窗口(ms),相似速度计算个数,叠加速度函数个数,每个叠加条带的CDP数,速度计算方式,速度在时间0处的变化,速度在最大时间处的变化,拷贝拾取到下一位置?,动校正最大拉伸百分比,设置相似比例和自动拾取参数?,从存在的参数表中得到指导函数?,层速度,是否设置可见项?,(3) 3D Velocity Viewer/Editor*,Disk data Input select datasetcdp sort100-1500(50)3,5-1500(50)3,5 Ensemble Redefine Band pass filter Automatic Gain Control Veloci
43、ty Analysis is the data precomputed.yes,速度分析:(1)预计算;(2)应用预计算结果。,(1)流程,(2)应用预计算结果,速度扫描,6、一次叠加、拾取层位、计算剩余静校正量 Disk Data Input Normal Moveout Correction Automatic Gain Control CDP/Ensemble Stack Disk Data Output Disk Data Input Bandpass Filter Auotomatic Gain Control Stack Display 2D/3D Max Power Autostatic*,7、二次速度分析 Disk Data Input Apply Residual Statics Bandpass Filter Automatic Gain Control Disk Data Output CVS Strip Analysis*,7、二次速度分析 Disk Data Input Apply Residual Statics Bandpass Filter Automatic Gain Control Disk Data Output CVS Strip Analysis*,
链接地址:https://www.31doc.com/p-2772892.html