地震勘探成图方法速度研究ppt课件

1、之一 地震勘探成图方法速度研讨 研讨院勘探三室汇报人：覃新平 地震勘探成图方法速度研讨 一、概 况 二、构呵斥图方法速度分析的方法原理 三、速度计算软件的研制 四、成图效果分析及存在问题一、概略1. 成图方法速 度研讨现状2. 速度研讨缘由3. 速度研讨内 容及义务4. 速度研讨取 得的成果 针对台北凹陷大规模的速度场研讨是在1994年4月2019年7月由物探局一处为主体担任进展的，完成了台北凹陷除北部山前带和火焰山构造带之外的主要地震反射界面的速度场研讨任务，编制了二维地震解释工区的速度分析微机软件，然而软件的稳定性较差，但其研讨的思绪，采用的新技术、新方法对我们的速度分析任务依然有很大的指

2、点意义。 一、概略1. 成图方法速 度研讨现状2. 速度研讨缘由3. 速度研讨内 容及义务4. 速度研讨取 得的成果 目前，各家勘探单位都有自已的速度研讨人员和速度计算软件。一、概略1. 成图方法速 度研讨现状2. 速度研讨缘由3. 速度研讨内 容及义务4. 速度研讨取 得的成果 地震勘探需求进一步深化，那些复杂构造如台北凹陷北部山前带、火焰山七克台构造带的构造圈闭落实需求专题论证和研讨一、概略1. 成图方法速 度研讨现状2. 速度研讨缘由3. 速度研讨内 容及义务4. 速度研讨取 得的成果 另外，在横向速度变化猛烈的三维勘探区域，如何实现较准确的变速成图，寻觅低幅度圈闭，日益成为迫切需求处理

3、的问题一、概略1. 成图方法速 度研讨现状2. 速度研讨缘由3. 速度研讨内 容及义务4. 速度研讨取 得的成果 2、横向速度变化猛烈的三维地震工区的较准确变速构呵斥图方法。1、逆掩断层下盘构造形状的成图方法内 容一、概略1. 成图方法速 度研讨现状2. 速度研讨缘由3. 速度研讨内 容及义务4. 速度研讨取 得的成果任 务1、三塘湖条湖马朗凹陷，台北凹陷怡东地域的速度场研讨，提供构呵斥图方法2、神东、照南、葡萄沟三维工区的速度场计算任务3、TS94-158，ST98-501两条测线的模型正演任务4、软件研制方面，一是完成普通二维工区的速度场计算程序，二是完成三维变速成图计算软件一、概略1.

4、成图方法速 度研讨现状2. 速度研讨缘由3. 速度研讨内 容及义务4. 速度研讨取 得的成果 1. 利用20条测线，353个速度谱点，4口井测井资料，5层解释层位，完成台北凹陷怡东地域速度场研讨，提供分区块速度成图方案2. 利用43 条测线，735 个速度谱点，6口井 VSP资料，3口测井资料， 6层解释层位，划分条湖马朗凹陷为8条速度变化趋势带，以及相应的时深关系数据3. 完成神东、照南、葡萄沟三维工区的层位平均速度计算任务4. 对TS94-158，ST98-501两测线进展正演计算，初步得到两剖面的深度剖面 取 得 成 果一、概略1. 成图方法速 度研讨现状2. 速度研讨缘由3. 速度研讨

5、内 容及义务4. 速度研讨取 得的成果5. 完成了模型迭代法在二维工区的层速度计算程序，完成三维工区沿主测线进展了模型迭代法计算程序，完善了三维全自动变速成图程序6. 完成了报告7. 完成了报告 取 得 成 果 地震勘探成图方法速度研讨 一、概 况 二、构呵斥图方法速度分析的方法原理 三、速度计算软件的研制 四、成图效果分析及存在问题 二、构呵斥图方法速度分析的方法原理 1. 层速度计算方法 2. 二维、三维工区的构呵斥图方法 3. 数据处置方法 众所周知，层速度是一个重要的地层参数，在程度叠加、偏移、压力、储层预测、构呵斥图中得到广泛运用。 层速度小知识 层速度计算方法 地 下 地 质 模

6、型T1: V=2000m/sT3: V=4500m/sT4: V=6000m/sT5: V=4800m/sT6: V=3700m/sT2: V=4000m/s 层速度计算方法：模型正演 利用野外单炮 单 炮 射 线 路 径 炮点 层速度计算方法：模型正演 利用野外单炮 层速度计算方法：模型正演 利用野外单炮 理 论 时 距 曲 线 实 际 野 外 接 收 记 录反射轴相关最好，层速度模型确立 层速度计算方法：模型正演 利用途置剖面红层为实践剖面层位黑层为正演结果层位正演深度剖面 实践时间剖面与正演时间剖面对比图 层速度计算方法 模型反演DIX公式 假设地层为程度层状介质，层速度Vi与均方根速度

7、Vrms有如下关系Vi,n2 = Vrms,n2To,n - Vrms,n-12 To,n-1To,n - To,n-1 层速度计算方法 模型反演双曲斜率法V1V2 SA=x1, SB=x2V2 2 = x2 - x1T2 - T1dxdT 层速度计算方法 波阻抗反演法 目前，另一个子课题“地球物了解释方法研讨之储层反演正在运用 层速度计算方法 模型反演模型迭代法 恣意倾斜层状介质射线途径表示图 又称为射线追踪法，是建立在地震波运动学根底上的。该方法的根本思绪是将地下地层笼统为倾斜层状介质，用射线追踪的方法求出初始层速度模型及其共中心点时距曲线，然后模拟做速度谱的方法，用实际曲线对它进展拟合求

8、出叠加速度，把它和实践叠加速度进展比较，误差小于给定的范围以内，即以为所给定的模型是正确的，此点的层速度值因此确定。然后进展下一个谱点的计算。Xk 的递推公式Zk 的递推公式k 的逆递推公式T 0 , k 的递推公式M. Shah (沙赫 ) 公式 层速度计算方法 模型反演模型迭代法 层速度计算方法 以葡1井为例，阐明DIX公式与模型迭代法计算的结果对比 葡1井叠加速度谱图解释叠加速度曲线 过葡1井地震测线：T90140B500600700800900100011001200130014001500160017001800190020002500260027002800290030003100

9、3200330034003500图 例DIX 公式平均速度模型迭代平均速度VSP 平均速度时间平均速度（ 米/ 秒） 葡1井三种平均速度对比效果地震层位To时间VSP深度DIX公式计算深度模型迭代法计算深度DIX公式相对误差模型迭代法相对误差N1t1080 1574159915761.60%0.13%Esh1215 1800183718152.06%0.83%K1280 1925196119351.87%0.52%J3q1370 2100214321202.05%0.95%J2q1425 2210225422301.99%0.91%J2s11540 2425248324502.39%0.99%

10、J2x4+J2x31600 2550260625752.20%0.98%J2x2+J2x11700 2760283628032.75%1.56% 葡1井计算深度误差表模型迭代法计算的平均速度更接近与VSP速度 平均速度对比图 层速度计算方法 模型迭代法计算精度高于DIX公式 计算方法比较 层速度计算方法单炮质量好、信噪比较高可操作性强，实际可靠 适用程度、平缓地层 另外一个子课题运用适用各种地层倾角，精度较高 层速度计算方法 选用 选用 二、构呵斥图方法速度分析的方法原理 1. 层速度计算方法 2. 二维、三维工区的构呵斥图方法 3. 数据处置方法 构呵斥图方法速度分析的方法原理二维、三维工区

11、的构呵斥图流程图 二维、三维工区的构呵斥图方法 自激自收时间对比接收点双程旅行时间(ms)A328.57B312.92C248.57A328.57B305.71C248.57抽象后模型实际模型 运用模型正演计算层速度的表示图简化了射线途径，呵斥时间计算的误差 二维、三维工区的构呵斥图方法 逆掩断层下盘 利用途置剖面，进展模型正演 平缓、倾斜地层 模型迭代法 二、构呵斥图方法速度分析的方法原理 1. 层速度计算方法 2. 二维、三维工区的构呵斥图方法 3. 数据处置方法 数据处置方法 正确解释速度谱 VSP测井及声波测井资料对速度场的校正 层速度的平滑 叠加速度谱解释 构呵斥图方法速度分析的方法

12、原理数据处置方法 A 能量团B 叠加效果 C 阅历 构呵斥图方法速度分析的方法原理数据处置方法 层速度平滑二维平面上五点三次公式平滑 ST94-L-26测线P2h层速度平滑图平滑后层速度平滑前层速度 构呵斥图方法速度分析的方法原理数据处置方法 层速度平滑 方格网平滑 马朗凹陷P2h平滑前层速度平面图 马朗凹陷P2h平滑后层速度平面图构呵斥图方法速度分析的方法原理 数据处置方法 层速度平滑 可见，层速度数据经过平滑后，其规律性得到加强 构呵斥图方法速度分析的方法原理数据处置方法 VSP速度对平均速度校正表示图校正前速度曲线校正后速度曲线 构呵斥图方法速度分析的方法原理数据处置方法 声波测井曲线校

13、正方法 通常，利用声波测井资料对速度场进展校正之前，要进展环境校正 环境校正 泥岩侵蚀校正 周波腾跃校正 构呵斥图方法速度分析的方法原理数据处置方法 声波测井曲线校正方法 岩层受泥浆侵蚀后声波测井表示图 当接纳探头首先接纳到的是来自界面上的滑行波时，没有时差误差 当接纳探头首先接纳到的是来自界面上的滑行波时，产生时差误差值到达最大 当接纳探头接纳到的是来自界面至界面之间的滑行波时，产生时差误差 接纳情况1 接纳情况2 接纳情况3 构呵斥图方法速度分析的方法原理数据处置方法 声波测井曲线校正方法 泥质岩石 声波的周波泥质岩石的声波时差校正公式：tc = to + Q ( C1to + C2 )

14、- to受泥浆侵蚀 腾跃景象 反复中值滤波 构呵斥图方法速度分析的方法原理数据处置方法 声波测井曲线校正方法 玉东1井声波校正前的合成记录 玉东1井声波校正后的合成记录 构呵斥图方法速度分析的方法原理数据处置方法 声波测井曲线校正方法 杜1井声波校正前的合成记录 杜1井声波校正后的合成记录 构呵斥图方法速度分析的方法原理数据处置方法 声波测井曲线校正方法 玉东1井各种速度对比图 校正后的声波时差曲线速度更接近于 VSP 速度 构呵斥图方法速度分析的方法原理数据处置方法 声波测井曲线校正反复中值滤波 中值滤波后的声波时差曲线 中值滤波所滤掉的声波时差曲线 构呵斥图方法速度分析的方法原理数据处置方

15、法 声波测井曲线校正方法 校正后的声波时差速度更接近于地层的实践速度，这样，大大提高了声波时差曲线的可信程度 地震勘探成图方法速度研讨 一、概 况 二、构呵斥图方法速度分析的方法原理 三、速度计算软件的研制 四、成图效果分析及存在问题 速度计算软件的研制 二维工区计算程序 二维工区速度计算及配套软件的功能 1. 叠加速度数据的检错、阅读 2. 计算剖面倾角视差，反射波射线追踪 3. 计算层速度，对层速度进展平滑、校正 4. 对声波资料进展校正，并绘制合成地震记录 速度计算软件的研制 三维工区计算程序 三维工区速度计算及软件的功能 1. 沿提取速度谱方向进展模型迭代计算，在井数据约束 下，得到层

16、速度、平均速度、校正后平均速度、速度校正量等数据文件 2. 沿地震反射层提取各种速度 3. 地震反射层位的时深转换 4. 沿时间切片提取各种速度，提取十字剖面速度，提取井旁道速度 速度计算软件的研制 程序运转环境 二维速度计算只需求满足：80386 或80386以上微机、单色或彩色显示器、绘图仪、数字化仪，硬盘容量为 50 MKB以上，内存容量为4MB以上 硬件配置： 三维速度计算在SUN4或SUN4以上任务站上运转 软件配置： DOS6.0以上版本的操作系统，QB4.5编译器，任务站上运用UNIX-FORTRAN言语编译器 地震勘探成图方法速度研讨 一、概 况 二、构呵斥图方法速度分析的方法

17、原理 三、速度计算软件的研制 四、成图效果分析及存在问题 成图效果分析 条湖马朗凹陷 钻 井 层 速 度 对 比 图高值区主要是石炭系、火成岩低值区浅层中低、中高值区各系组 成图效果分析 条湖马朗凹陷 马朗凹陷P2h层速度平面图 速度梯度为 50-100米/秒*千米，变化较为平缓 成图效果分析 条湖马朗凹陷 钻 井 平 均 速 度 对 比 图 1秒时速度相差达1400m/s 1.5秒时速度相差达800m/s 东北部速度梯度为10-150米/秒*千米 西南部速度梯度为200-350米/秒*千米 成图效果分析 条湖马朗凹陷 成图效果分析 条湖马朗凹陷 平均速度=5000米/秒区间 成图效果分析 条

18、湖马朗凹陷 三塘湖条湖马朗凹陷速度分析阐明： 东北部速度变化较小，采用恒速塘参1井平均速度成图 西南缘速度变化很大，采用分区块变速成图 成图效果分析 条湖马朗凹陷 平面上划分8条速度变化趋势线，各条速度线之间为线性变化速度变化趋势线图 成图效果分析 条湖马朗凹陷 趋势线平均速度曲线 8条速度趋势线的平均速度与时间关系 成图效果分析 条湖马朗凹陷 马6井为例阐明变速成图的效果 成图效果分析 条湖马朗凹陷 趋势线平均速度曲线马5井声波速度恒速成图速度 变速成图时，马6井位于趋势线L5与L6之间 成图效果分析 条湖马朗凹陷层位底界深度(m)地层厚度(m)解释时间(ms)层速度(m/s)预测深度(m)

19、绝对误差(m)相对误差C127912795304826693-545.81%J2t1377985803920782-743.21%J2x15491726703822949-1438.73%T16479881014001227-19425.50%火成岩顶215851197561941570327.25%火成岩底2777619131036962317816.56%P2t30642871465370326651413.02% 马6深度预测误差表恒速成图 成图效果分析 条湖马朗凹陷层位底界深度(m)地层厚度(m)解释时间(ms)层速度(m/s)预测深度(m)绝对误差(m)相对误差C1279127953

20、048261284-50.40%J2t13779858039201384-70.50%J2x154917267038221563-140.90%T16479881014001841-19411.80%火成岩顶21585119756194215530.10%火成岩底277761913103696276980.30%P2t3064287146537033050140.50% 马6 井 位于速度趋势线L 5 与L 6 的中间 T层层速度为1400米/秒，偏小 马6深度预测误差表变速后 成图效果分析 条湖马朗凹陷 变速成图后，成图精度大大提高 成图效果分析 神东三维 crossline511 速度 I

21、nline 221 速度 成图效果分析 神东三维测线西北向平均速度略有升高 成图效果分析 神东三维西北向平均速度抬升 成图效果分析 神东三维 西北部平均速度场较高 南部平均速度较低 成图效果分析 神东三维 神东2井速度分析 成图效果分析 神东三维 神1井、神东2井平均速度曲线神1井VSP平均速度1 0 02 0 03 0 04 0 05 0 06 0 07 0 08 0 09 0 01 0 0 01 1 0 01 2 0 01 3 0 01 4 0 01 5 0 01 6 0 01 7 0 01 8 0 01 9 0 02 0 0 02 1 0 02 2 0 02 3 0 02 4 0 02

22、5 0 01800190020002100220023002400250026002700280029003000310032003300340035003600To时间平 均 速 度 （ 米 / 秒 ）98年成图平均速度神东2井平均速度 成图效果分析 神东三维 神1井、神东2井层速度曲线神1井层速度曲线神东2井层速度曲线 成图效果分析 神东三维 神东地震反射神东地震反射TEshTEsh层构造图层构造图9999年度年度神东2井预测深度为例，阐明成图效果 神东地震反射神东地震反射TEshTEsh层构造图层构造图9898年度年度 成图效果分析 神东三维层位钻井深度(m)98年成图深度(m)99年成

23、图深度(m)98年成图误差99年成图误差TEsh195219861961 1.74% 0.46%TK1s1228123422310 2.67% 1.27% 神东2井98、99年深度预测误差表 成图精度比去年提高1个百分点左右 葡萄沟地域地震反射葡萄沟地域地震反射Tj2qTj2q层构造图层构造图(97(97年度年度) ) 成图效果分析 葡萄沟三维 葡萄沟地域葡萄沟地域9999年度七克台组构造图年度七克台组构造图上盘上盘下盘下盘 成图效果分析 葡萄沟三维 葡萄沟地域成图精度对比表葡萄沟地域成图精度对比表99年成图精度年成图精度 1.1%97年成图精度年成图精度1.1%7.6%精度提高0.5%|1.5% 成图效果分析 葡萄沟三维 成图效果分析 今年，我们加强了速度成图方法的研讨，使成图精度比去年有所提高，到达或超越了部颁规范 利用途置程度剖面，进展模型正演 成图效果分析 模型正演初步运用 逆掩断层下盘为平缓时的射线途径 模型正演初步运用 逆掩断层下盘为平缓时的正演剖面 下盘时间抬升，产生假断层、假构造 模型正演初步运用 逆掩断层下盘为上倾时的正演剖面 出现空白段 模型正演初步运用 逆掩断层下盘为下倾时的正演剖面 出现小幅度断层，时间稍有抬