[工学]promax处理流程.ppt

《[工学]promax处理流程.ppt》由会员分享，可在线阅读，更多相关《[工学]promax处理流程.ppt（108页珍藏版）》请在三一文库上搜索。

1、地震资料数字处理,实验,学时：8 每组32人,指导教师：李谋杰,一、资料概况 1、工区概况 该工区位于天津市大港区千米桥以西沈青庄地区。工区北起北大港水库泄洪区，向南穿越红旗路，西起大中旺乡，东抵幸福小区采油作业一区一带。工区中北部几乎全部位于北大港水库、泄洪区内，东南部进入矿区。整个工区地表分布有河流、苇塘、作业区、生活区、交通干线、农田等，地表情况十分复杂。,2、施工信息,三维地震资料采集，设计满覆盖面积为249.17平方公里；共计17600炮。,测量工作进行了123天，共测三维线束31束，其中排列线128条；炮线260条；检波点定位点数38277个；炮点定位点数18148个。,每束线包括

2、8条检波线和8条炮线；每相邻两束线间重复4条检波线；线号和桩号均以50米为单位；1800以上为线号；1000以下为桩号；分别由西向东、由南向北递增，测线方位角为 。,坐标系：北京54坐标系，原点纬度 ,原点经度 。高程基准面：1956年黄海高程系。,3、实验资料,选择第17束线进行实验，测量资料和班报经整理后，存放在如下三个文本文件中： 17s.ref检波线文件；(线号：22562284；桩号498825) 17s.sin炮线文件； (线号：22632277；桩号498825) 17s.pat排列文件。 (136道/线；炮序号1352),解编后的地震数据存放在17031708(共6个)文件中，

3、其记录格式为SEGY格式；2ms采样；1088道/炮；3000样点/道；FFID：141560 (实际使用范围：142541),排列线距200米；道间距50米。 炮点线距100米；垂直炮线炮点距200米 纵测线(Inline)号的计算方法： (1) 计算(对应排列线号+对应炮线号)/2,其结果为Li； (2) 得到相同Li的次数为覆盖次数； (3) 对整个工区内的Li重新按整数编号或不重新编号。 横测线(crossline)号的计算方法： (1) 计算(每个接收点桩号+每条炮点桩号)/2,其结果为Xi； (2) 按道距之半整数编号或不重新编号，即CMP的编号。 面元边长 道距是接收线方向面元边

4、长的2倍，CMP点距=b.,2256 2260 2264 2268 2272 2276 2280 2284,498,825,横测线号: 824.25,横测线号: 523.75,炮线号: 2263 2265 2267 2269 2271 2273 2275 2277,2 2275 823.5,10 2263 823.5,共22条共中心线,(2256+2263)/2=2259.5,(2256+2265)/2=2260.5,(2284+2277)/2=2280.5,(2256+2275)/2=2265.5,(2256+2267)/2=2261.5,(2260+2263)/2=2261.5,移动200

5、米即4道,68道3350米,68道,450米9道,总覆盖次数 纵线方向覆盖次数：,横线方向覆盖次数：,4、实验内容,(1)定义观测系统使CMP道集分选正确 (2)叠前去噪剔除不正常的炮、道和野值；压制强面波和其他不规则干 扰，使信噪比较高，波组特征清楚。 (3)振幅补偿做相对保持振幅处理,使浅、中、深、道、炮间能量均衡 (4)反褶积用好多道反褶积，实现零相位输出，达到压缩子波、提高地 震记录分辨率。 (5)静校正检查静校量空间分布的合理性。 (6)速度分析用速度谱或速度扫描获取速度,每1km不少于2个速度谱。 (7)动校正和叠加叠加速度合理切除参数正确最终叠加剖面质量较优 (8)倾角时差校正断

6、面波、绕射波、水平及倾斜反射成像清楚。 (9)叠后时间偏移偏移速度有效波归位合理,绕射波收敛,断点清楚。,二、UNIX系统常用命令,三、启动ProMax系统 1、启动仿真软件Exceed,2、选择主机lade1 双击鼠标左键,、输入用户名和口令,用户名: pmxusr0116,口令: 123456,、设置环境变量 setenv DISPLAY 192.168.1.21:0 、启动ProMAX promax 或 promax3d 6、建立工区、测线和流程目录 ，即：, Area 工区菜单 Select 选择一个存在的工区； Add 加一新工区,即创建工区目录； Delete 删工区； Renam

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

8、) packet.job: 存处理列表和参数 job.output: 存放流程执行产生列表,用以下作业核实/export/blade1/disk1/userdata/dty/dg_sqzh/1703文件中的FFID,Datasets数据集，击MB2可显示信息：数据 集的大小(字节)、数据类型、总道 数、样点数、采样率、最小和最大 CDP数、炮点索引号等。击MB3可 显数据集的处理历史。 Database数据库 7、选择模块组织作业,MB1,MB2,MB3,四、处理流程 1、定义观测系统 运行如下模块: 3D Land Geometry Spreadsheet* 弹出如下菜单:,按从左到右顺序执

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

10、s按炮点和模式表中的线号和站号来匹配。 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 No Surface

11、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击,则将数据输入表中。, Sources (填写炮点参数表）,注意： 有意将 表中 坐标与文件中 坐标填写不一致, Patterns (填写排列或班报参

12、数表）, Bin (数据面元化) 必须执行下面1、3、5项,1、指定中点：按SIN和PAT表中的线号和站号来匹配 2、定义面元网格 3、面元中点 4、面元数据质量控制 5、数据库定案,计算炮点到检波点间的距离； 计算炮点到检波点的中点坐标； 计算炮点到检波点的方位角,CDP=炮点+检波点,本例网格原点位置（min inline & crossline,Xorg=532051.50 Yorg=4282645.50,网格的方位角 面元边长dx 面元边长dy,网格的坐标原点,Inline 间的距离,crossline 间的距离,存放本表的文件名,Cdp、inline、crossline号,形成CMP

13、 道的最小/最大炮检距,炮检距面元之间距离,震源-接收点间方位角,方位角公差,当定义面元化参数后，计算网格维数,沿网格Y轴的方位角 面元在x方向的边长，即dx 面元在Y方向的边长，即dy 3D网格的坐标原点 网格在x方向的距离(/dx= x方向上的线数) 网格在Y方向的距离(/dy= y方向上的线数) 存放网格参数的文件名 起始CDP号 起始纵线号 起始横线号 距面元中心的最小炮检距 距面元中心的最大炮检距 面元化时炮检距增量 震源接收点的方位角 震源接收点的方位角公差 纵线平行于网格的x轴 纵线平行于网格的y轴,Dx=50m,Dy=25m,Bin=面元,击Database,弹出,TRC记录道

14、 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 6 、检查SIN上的高程分布：SIN-Geometry ELEV,覆盖次数,按Inline-xline显示的覆盖次

15、数和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 格式：,卷头(4字节/字,共100字）:,道头(4字节/字,共60字）:,SEG-Y 类型 使用的存储类别 输入磁盘文件名

16、 当采样间隔不正确时,可忽略输入数据的采样间隔,给出其正确值. 每道的样点数(在忽略2#道头时给出) 处理历史存入卷头吗? 输入辅助道吗? 从道头中得到道号吗? 输入道格式 显示道集信息吗 输入的最大时间 是叠加数据吗? 每集的最大道数 主分选道头字,输入第一选择吗? 输入第二选择吗 输入全局的XY参考坐标吗?若为yes,则道头中xy减此值. 使用坐标比例 重置道头吗?,选道: 1-10(2) “(2)”表示隔一道读一道.,道选择方式: Include包含 Exclude拒绝 列出指定的道,连机几何信息道头装入 与数据库对应的主道头字 用有效道号来匹配吗？ 去掉道头中为空CDP的道吗？ 去掉道

17、头中为空接收点的道吗？ 冗余诊断吗？,道头运算 选择方式 定义道头运算表达式,固定等式方式 顺序重编号方式 旋转与平移方式 改变数据属性,磁盘数据输出 指定输出数据集的文件名 是新文件还是已存在的文件？ 输出的记录长度 指定样值存储格式 跳过主磁盘存储区吗？,8、16位 为整型数，每一百个样值有一 比例因子； 32位 为IEEE格式的浮点数。,3、叠前参数测试,(1)获取初至切除、去废道废炮参数 Disk Data Input Automatic Gain Control Trace Display,磁盘数据输入 从其他工区、测线中读取数据吗？ 选择数据集 繁衍输入文件的处理历史吗？ 读道选件

18、: 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，10 43 表示取顺序值,即4,5,6; 又如130(10)3,即1,2,3,11,12,13,21,22,23 *

19、 通配符，其范围-999999.0999999.0; / 分组符如10-14(2)/15-21(3)表示取10,12,14,15,18,21; : 分级符如：1100：116 (前为主头字，后次道头字),应用方式: 见 增益因子的类型, 见 增益因子长度 (ms) 因子的应用基准, 见 拒绝硬化零吗? 鲁棒比例吗?(强行求因子),只应用增益因子,不保存于道头中 应用增益因子,并将因子保存于道头中 去掉增益控制,并删除道头中增益因子 输出每道的增益因子,自动增益控制,应用基准,绝对平均值 均方根值 中值 (按滑动时窗求上面三值),中心 头部 尾部,道显示 选择显示设备 显示的起止时间(ms) 每

20、屏的最大道数 每屏间的重叠道数 每屏的道集(线段)数 使用变道间隔吗？ 输出方式 道显示方式 变面积基线 道头显示参数 自动保存屏幕吗？ 保存屏幕图像的最大数 按彩色方式保存屏幕吗？ 保存屏幕图像到何处 收集的屏幕数 道显示方向/极性 标注的主/次道头字 次道注释方式/及增量 道比例方式(习惯的),剪切的道偏离值 样值所乘的标量值 道比例选件(单/多道) 道显示方向(垂直/水平),SOURCE:有效炮号(用户定义的) FFID:野外文件标识号 SIN: 炮号(内部号) REC_SLOC: 检波点号 SOU_SLOC: 外部震源位置号 SRF_SLOC: 外部检波点位置号 CDP: CDP面元号

21、 OFFSET: 带符号的炮检距 AOFFSET: 绝对炮检距 CHAN: 记录道号 REPEAT: 重复数据的复制号 CS_STAT: 炮点相关_自动统计静校正量 CR_STAT: 检波点相关_自动统计静校正量 PS_STAT: 炮点最大能量_自动统计静校正量 PR_STAT: 检波点最大能量_自动统计静校正量 TRIMSTAT: 修饰静校正量 FB_PICK: 初至拾取时间 TR_FOLD: 实际道的覆盖次数 DMOOFF: DMO偏移距面元 ILINE_NO: 纵测线号 XLINE_NO: 横测线号 CDP_SLOC: 外部CDP位置号,删除道 极性反转 拾取顶部切除 拾取底部切除 拾取

22、手术切除 拾取各种各样的时窗 编辑道头值 编辑数据库值 拾取其他层位,道切除 重新应用以前的切除吗? 切除时间参考 切除类型 斜坡时间(ms) 外推切除时间 从数据库得到切除文件吗? 选择切除参数文件,(1) 参数表拾取 disk Data Input Read datafrom other lines/surveys?.No Select datasetShots-with-geometry trace read optionsort interaction Data Access?yes select primary trace header entryFFID select second

23、ary 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.yes,运行此流程拾取如下参数表： Kill trace Pick Top mute Pick Miscellan

24、eous 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/Reverse Trace editing MODEReverse Get edits form the DAT

25、ABASE?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 corrections?.No Get TAR velocity function from database?No Sho

26、uld 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 of amplitude correction?.Apply Maximum appliction TIME:55

27、00 Trace Display Number of Ensembles(line segments)/screen:2 Trace scaling option:Entire Screen Number of display panels:2,说明： 球面扩散(spherical divergence)地震波的振幅随传播距离的增大而减小的现象，称为波前或球面扩散。,非弹性吸收(inelastic attenuation)由于岩石是非完全弹性体，地震波在传播过程中一部分能量将转化为热能而消耗掉，从而使波的能量变弱的现象，称为地层的吸收衰减作用。,真振幅恢复(True Amplitude Rec

28、overy)就是从地面检波器记录的原始振幅中，消除波前扩散和地层吸收衰减的影响，使其恢复到仅与地下反射界面的反射系数有关的真振幅值，使恢复后的地震波振幅正比于反射系数的大小。这种保持反射波振幅与反射系数呈正比的剖面称为相对振幅保持剖面。,(4)反褶积参数测试 Disk Data Input Trace kill/reverse Trace muting True Amplitude Recovery dB/sec correction constant .9 Reproduce Traces Trace grouping to reproduce.Ensembles Total Number

29、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 trrace header wordREPEAT Selectr secondary trace header wordN

30、ONE 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 level.0.1 Window rejection factor2 Time gate referencetime Ge

31、t 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,处参数说明： Spiking/Predictive Decon Apply a bandpass filter after d

32、econ?.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 Panel width in traces(面板宽度-道数) .136 Starting time for ana

33、lysis(分析的起始时间) 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 (零波数位置) Top Plot FK,TK,TX panels in DB or Linear?(以分贝/线性绘) DB

34、SCALE 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 (时间斜坡平坦百分比) 100 (不斜坡) Select mute polygon table (选择切除多边形表) fk

35、 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 (输入道的间距)50 Panel width in trace (以道计的面板宽度)136 Test the f

36、ilter 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 of F-K filter operationReject Percent flat for F-K filter

37、 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 Trace Display Label ENDIF Trace Display Number of Ensemb

38、les(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 dataset.16/ Automatic Gain Control Interactive Spectral A

39、nalysis 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 wordRepeat Specify trace list.1 ELSEIF Select Primary trace head

40、er 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 Subset(多子集方式) Freeze the selected subset?(冻结所选子集)yes Displ

41、ay data by traces or ensemblesEnsembles,4、高程静校正,四种基准面的定义 统一基准面：是一个由用户提供的与海平面平行的水平面，最终叠加和偏移结果必须校正到统一基准面。 CMP参考面：是一个以CMP为单元统计出的等时面，其静校正量为该CMP道集中的炮点静校正量和接收点静校正量的算术平均值。速度分析、动校正和叠加在此面上进行。 辅助基准面：是由用户临时定义的接近地表的一个水平面，一般在偏移处理时使用。 浮动基准面：对CMP高程经过平滑后所得到一个平滑面(非水平面)。 (2)约定 校正量定义：,校正量应用：当校正量为正值时，记录往上提；当校正量为负值时，记录往

42、下拉。 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 of smoother51 Processing DATUMNMO DATUM Run ID01 用DBTools

43、选择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 Muting True Amplitude Recovery F-K filter Spiking/Predict

44、ive 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 Dataset FilenameShot-decon/elev statics,5、速度分析 (1)数据准备 Disk Da

45、ta Input (按CDP分选) Bandpass Filter (窄带滤波) Automatic Gain Control(l300-500) Disk Data Output (2) 常速扫描 CVS Strip Analysis*,(2) 速度分析 Disk Data Input,Velocity Analysis,选择显示设备,指定存储速度拾取的表,为计算变量置初值吗？,求和成道集的CDP数,是否对面元应用动校正？,分析数据预先计算吗？,是否应用静校正？,相似速度采样率,相似速度计算窗口(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 p