《地球物理解释基础.ppt》由会员分享,可在线阅读,更多相关《地球物理解释基础.ppt(89页珍藏版)》请在金锄头文库上搜索。
1、东方地球物理公司2007年5月地地 球球 物物 理理 解解 释释 基基 础础 (2)东方地球物理公司2007年5月第第1414章章解释与盐构造有关的圈闭解释与盐构造有关的圈闭 许多重要油田许多重要油田 和盐圈闭联系在一起和盐圈闭联系在一起著名的墨西哥湾、美国几个洲、加拿大、北海、北非、德国 、里海地区 都存在盐圈闭。塔里木盆塔里木盆地的克拉地的克拉2 2气田也和盐圈闭有关气田也和盐圈闭有关盐与众不同盐与众不同 有较低的密度和较高的地震速度较低的密度和较高的地震速度 “ “漂浮漂浮” ”状侵入到沉积物之下状侵入到沉积物之下侵入体产生各种盐体形状 ,盐体侧翼成倾斜状 ;盐侵入体之上形成断层圈闭;岩
2、盖上呈垂直盐株状 ;古老的盐丘有厚层堆积物 (石膏、碳酸盐岩)盐丘的地震勘探成像问题是关键盐丘的地震勘探成像问题是关键盐丘的3D形状,通常需要3D偏移。盐通常具有比围岩要高很多的P-波速度,围绕盐和周围沉积之间横向速度差大成为成像主要问题速度横向变化、三维形状的盐丘和陡倾角足以值得应用三维速度横向变化、三维形状的盐丘和陡倾角足以值得应用三维叠前深度偏移叠前深度偏移 东方地球物理公司2007年5月Subslat imaging via target-oriented 3D prestack depth migration By D. Ratcliff, C. A. Jacewitz, and S
3、. H. Gray 通过针对目标的通过针对目标的3D叠前深度偏移盐丘成像叠前深度偏移盐丘成像(墨西哥湾Vermilion构造的盐丘)东方地球物理公司2007年5月 (引自引自Ratcliff 等人,等人,1994 ) 盐丘数据的叠前深度偏移剖面盐丘数据的叠前深度偏移剖面,112次覆盖次覆盖东方地球物理公司2007年5月3D叠前深度偏移的必要性叠前深度偏移的必要性(a )2D叠后时间偏移反映出一个不完整和畸变的TOS成像(箭头)(b)3D叠后时间偏移,3D偏移消除了畸变偏移消除了畸变,盐顶清晰的成像盐顶清晰的成像 (引自Ratcliff 等人,1994 ) 2D、3D叠后时间偏移的比较叠后时间偏
4、移的比较东方地球物理公司2007年5月(a)2D叠前时间偏移显示了不正确的盐底( BOS )位置,并缺少盐的反射(b)2D叠前深度偏移剖面,盐的成像有相当改进 (引自Ratcliff 等人,1994 ) 2D叠前时间偏移叠前时间偏移与与2D叠前深度偏移的比较叠前深度偏移的比较东方地球物理公司2007年5月用用3D叠后偏移叠后偏移建立建立3D速度场速度场 (a)3D沉积层速度场沉积层速度场横剖面横剖面,横向速度变化很小(b)3D叠后深度偏移叠后深度偏移第一次迭代第一次迭代盐的 成像(TOS以上正确像)(c)用TOS修正的修正的3D速度场速度场(d)3D叠后深度偏移叠后深度偏移第二次迭代第二次迭代
5、(BOS以上正确成像)(e)用BOS修正的修正的3D速度场速度场(f)3D叠后深度偏移叠后深度偏移第三次迭代第三次迭代(盐下叠加后还保留的同相轴的正确成像)东方地球物理公司2007年5月建立建立3D速度模型速度模型沉沉积岩速度岩速度场 井的控制井的控制 3D DMO速度场 3D MBS(叠前偏移)速度分析 3D 叠后深度偏移3D 盐和沉和沉积层速度速度场 3D 叠后深度偏移叠后深度偏移 GOCAD 3D 速度包应用井的信息、3D DMO(倾角动校正)速度信息、2D叠前偏移速度分析信息和3D叠后深度偏移,来建立3D沉积层速度场沉积层速度场用3D叠后深度偏移,应用3D设计软件来建立盐和沉积层的盐和
6、沉积层的3D速度场速度场东方地球物理公司2007年5月3D叠前深度偏移流程图叠前深度偏移流程图(在建立了在建立了3D速度场后应用)速度场后应用) 野外数据 重采样和编辑 球面扩散校正 切除 反褶积 滤波 数据与导航数据合并 3D共炮检距选排 3D共炮检距偏移 输出纵、横测线子集东方地球物理公司2007年5月覆盖次数对比覆盖次数对比3D叠前深度偏移叠前深度偏移(a)单次覆盖单次覆盖(1325m)(b)9次覆盖次覆盖(1275-1400m)(c)47次覆盖次覆盖(900-2100m)(d)112次覆盖次覆盖(引自(引自 Ratcliff 等人等人, 1994 ) 东方地球物理公司2007年5月炮检
7、距对比炮检距对比 3D叠前偏移,用以下列炮检距叠前偏移,用以下列炮检距:(a)炮检距值范围)炮检距值范围1300-2000m,(b)炮检距值范围炮检距值范围375-2000m。包含了近炮检距,改进了。包含了近炮检距,改进了TOS的成像的成像(引自(引自 Ratcliff 等人等人, 1994 ) 东方地球物理公司2007年5月TOS面面AVO合成记录合成记录 TOS交界面上的交界面上的AVO合成记录。合成记录。TOS的响应是来自所有炮检距;的响应是来自所有炮检距;向右随炮检距增加振幅减小(引自向右随炮检距增加振幅减小(引自 Ratcliff 等人等人, 1994 )(盐顶的反射能量大部分来自近
8、炮检距盐顶的反射能量大部分来自近炮检距) 东方地球物理公司2007年5月炮检距比较炮检距比较 3D3D叠前深度偏移叠前深度偏移,50,50次覆盖次覆盖, , ,比较炮检距范围对盐成像的影响比较炮检距范围对盐成像的影响 炮检距:炮检距:( a ) 375-1600 m 和和 (b ) 375- 3000 m (来自来自Ratcliff 等人等人, 1994 ) 东方地球物理公司2007年5月2D 、3D叠前深度偏移比较叠前深度偏移比较2D叠前深度偏移,显示了剖面叠前深度偏移,显示了剖面平面外的平面外的TOS,BOS不好不好 3D叠前深度偏移叠前深度偏移 TOS和和BOS都能正都能正确成像钻井穿过
9、清晰成像的盐背斜确成像钻井穿过清晰成像的盐背斜 东方地球物理公司2007年5月第第15章章地地 震震 模模 拟拟地震模型和地震模拟的概念地震模型和地震模拟的概念(Seismic Model 、Seismic Modeling)Seismic Model是名词, Seismic Modeling是动名词模拟模拟是去建立模型或模型响应的方式方法和过程是去建立模型或模型响应的方式方法和过程模型模型“:它可以推演出能与观测结果比较效果的一种概念,用于更好地:它可以推演出能与观测结果比较效果的一种概念,用于更好地理解观测结果。分为概念模型、物理模型或数学模型。理解观测结果。分为概念模型、物理模型或数学模
10、型。” ” (Sheriff, 1991) 地震模型能以一维(地震模型能以一维(1D)、二维()、二维(2D)或三维()或三维(3D)形式变化。这些模型与)形式变化。这些模型与实际情形的精确度取决与地质环境的吻合程度实际情形的精确度取决与地质环境的吻合程度地震模拟试图模拟地下的岩石性质和波在地下传播时地震波的地震模拟试图模拟地下的岩石性质和波在地下传播时地震波的传播响应传播响应地震地震- -模拟方法除维数外还包括不同的方法模拟方法除维数外还包括不同的方法地震模拟也可分为正演模拟和反演模拟地震模拟也可分为正演模拟和反演模拟模型的选择是在成本和模型的有效性之间取得平衡模型的选择是在成本和模型的有效
11、性之间取得平衡东方地球物理公司2007年5月 模模 拟拟 方方 法法模拟类型模拟类型数学模型数学模型一般性一般性费用费用法线入射反射系法线入射反射系数数1D反射系数值由下式给出 2211R= 22 + 11 对水平层和垂直旅行的波是有效的。保留了多次波对计算反射率很便宜,如果有多次波稍微贵一些 振幅随炮检距振幅随炮检距(或入射角)变(或入射角)变化化 “1.5D”地下模型是1D加非零(2D),用Zoeppritz方程严格地讲,对水平层是有效的,一般不包括多次波AVO模型比法线入射反射率花费多的多,但比波动方程求解便宜 射线追踪射线追踪按照Snell定理2D、3D求解。包括通过渐近线射线追踪的振
12、幅当非均质体的尺度与Fresnel带相比很大时,一般是可以应用的。通常忽略了绕射 多数情况下中等花费。通常计算射线路径很便宜,但计算振幅要增加费用波动方程有限差波动方程有限差分(分(FD)或有限)或有限元(元(FE)求解)求解1D、2D、3D数字表达为 1 2u2u = 2 t2FD 一般是矩形网格FE 更是一般用的,用网格算法费用一般性的价格物理模拟物理模拟1D, 2D, 3D需要将物理模型材料校正为成比例的模型 这种模型用物理材料提供建立模型本身是昂贵的,但模拟的运行通常比数值模拟要便宜 东方地球物理公司2007年5月地震模拟的用途地震模拟的用途 设计激发设计激发- -接收的观测系统接收的
13、观测系统偏离盐丘两种炮检距的偏离盐丘两种炮检距的VSP反射模型反射模型(引自(引自Whitmore 和和 Lines, 1986 ) 东方地球物理公司2007年5月地震模拟的用途地震模拟的用途解释工作结果预测(解释工作结果预测(用正演模拟和反演模拟用正演模拟和反演模拟)加强解释基础工作(加强解释基础工作(合成记录合成记录)数据处理大量应用反演数据处理大量应用反演(反褶积、静校正和速度估算是1D 模型,层析成像速度分析方法,是通过2D或3D模型 ,地震偏移,可以认为它是一种构造反演 )测试地震处理算法的正确性测试地震处理算法的正确性噪音影响测试噪音影响测试东方地球物理公司2007年5月零偏移距波
14、场映射(引自零偏移距波场映射(引自Whitmore 和和 Lines, 1986 ) 波动方程全解可以逼真地得到所有的波至,包括直达波、波动方程全解可以逼真地得到所有的波至,包括直达波、折射波、反射波和绕射波,所有的一次波和多次波折射波、反射波和绕射波,所有的一次波和多次波 东方地球物理公司2007年5月第第16章章地地 震震 反反 演演 正演模拟和反演的关系正演模拟和反演的关系 正演模拟正演模拟用一个数学关系式数学关系式,对给 出的一组模型参数模型参数合成地下响应。 反演或反演或“ “反演模拟反演模拟” ” 与正演 模拟 “ 相反”的过程。对一个给出的数据集,寻 求定义一个与观测数据相符的地
15、质模型地质模型 从数学上讲,反问题由于比方程式更多 未知数的存在,能够引起不确定性,产生 多解,所以反演的多解性是固有的反演的多解性是固有的 东方地球物理公司2007年5月反演的多解性(非唯一性、不确定性)反演的多解性(非唯一性、不确定性)(a)褶积模型的基本的数学多解性道褶积模型的基本的数学多解性道 (b)模型)模型1:震源子波(左)和脉冲:震源子波(左)和脉冲 响应响应 (右右)(c)模型)模型2:虚反射震源子波(左)和:虚反射震源子波(左)和 脉冲响应脉冲响应 (右)右) 噪音会引起大的变化或估算模型参数的噪音会引起大的变化或估算模型参数的不稳定,破坏解答的正确性不稳定,破坏解答的正确性
16、东方地球物理公司2007年5月1D 1D 模型的地震反演模型的地震反演 地震数据处理大部分的是基于近似水平层状地层地震数据处理大部分的是基于近似水平层状地层1D模型的假设模型的假设 ,包括动校正、水平叠加等包括动校正、水平叠加等 包括有密度、速度和厚度特征的一系列水平层的1D地质模型东方地球物理公司2007年5月反演技术反演技术地震反演技主要分四类:地震反演技主要分四类:(1)、基于地震数据的声波阻抗反演基于地震数据的声波阻抗反演(2)、基于模型的测井属性反演基于模型的测井属性反演(3)、基于地质统计的随机模拟与随机反演基于地质统计的随机模拟与随机反演(4)、叠前地震反演叠前地震反演 东方地球
17、物理公司2007年5月常用的反演常用的反演地震波阻抗估算地震波阻抗估算算法:算法:递归反演递归反演(早期的地震反演算法)(早期的地震反演算法) 可以可以从反射系数和上面层的阻抗推断下面地层的阻抗从反射系数和上面层的阻抗推断下面地层的阻抗。这个反演常常。这个反演常常叫作叫作Seislog反演反演 也可用也可用密度和速度之间的密度和速度之间的Gardner关系式关系式 将密度替换为速度,反演结将密度替换为速度,反演结果就变成速度函数果就变成速度函数合成声波测井曲线合成声波测井曲线 声波测井曲线与合成声波测井曲线之间的主要差别,是地震数据中缺声波测井曲线与合成声波测井曲线之间的主要差别,是地震数据中
18、缺少低频带宽(典型的是少低频带宽(典型的是0-5 Hz ) 另一个主要问题是另一个主要问题是缺少高频成分缺少高频成分,这是因为地震数据也缺少高频(有,这是因为地震数据也缺少高频(有代表性的是代表性的是100 Hz至至Nyquist频率)频率) 通常是通常是用现有的声波测井信息或用估算的层速度来重新获得低频用现有的声波测井信息或用估算的层速度来重新获得低频 东方地球物理公司2007年5月声波测井曲线可表示成速度声波测井曲线可表示成速度函数函数( 0-5 Hz )和精细的速度和精细的速度函数函数( 6-250 Hz)之和之和 (引自引自 Lidseth,1979)东方地球物理公司2007年5月从声
19、波曲线上去除高频成从声波曲线上去除高频成分导致降低分辨率的例子分导致降低分辨率的例子 ( (引自引自 Lidseth,1979)Lidseth,1979)东方地球物理公司2007年5月第第1717章章地震旅行时层析成像地震旅行时层析成像层析成像层析成像(tomography)“tomo”是希腊字,切片的意思,层析成像的意思是一个物体的切片图像医疗诊断的医疗诊断的CT技术技术原理是通过沿各个方向穿过人体的X射线,测量X射线的强度,确定人体不同部位的吸收性质地震旅行时层析成像地震旅行时层析成像是一种利用大量炮点和检波点综合是一种利用大量炮点和检波点综合观测结果求取速度与反射系数分布的方法观测结果求
20、取速度与反射系数分布的方法层析成像技术有两个假定前提条件层析成像技术有两个假定前提条件 假定物性是位置的连续函数假定物性是位置的连续函数 假定介质可离散化成有限数量均匀的面元假定介质可离散化成有限数量均匀的面元在地震旅行时层析成像中在地震旅行时层析成像中, ,地下介质被分解为面元层析的地下介质被分解为面元层析的目标是求解每个面元的速度目标是求解每个面元的速度东方地球物理公司2007年5月从炮点到接收点的射线路径是由位于不同面元中的射线段组成,根据各个面元射线段的长度和各个小面元的速度来计算旅行时由初始模型计算波至时间与观测值进行对比(正演),根据两者的时间差对模型进行修改,模型正演、测量时间差
21、、修改模型正演、测量时间差、修改模型这一迭代过程一直到时差模型这一迭代过程一直到时差小于给定值(最小平方差)小于给定值(最小平方差)东方地球物理公司2007年5月层析成像技术层析成像技术层析的目标是求解每个面元的速度层析的目标是求解每个面元的速度层析成像方法的第一步是从未叠加的地震资料上或直接从野层析成像方法的第一步是从未叠加的地震资料上或直接从野外观测值拾取旅行时外观测值拾取旅行时 建立初始模型作射线追踪,由初始模型计算的波至时间与观建立初始模型作射线追踪,由初始模型计算的波至时间与观测值进行比较测值进行比较根据模型值与观测值之差对模型进行修改根据模型值与观测值之差对模型进行修改拾取旅行时拾
22、取旅行时- -建立模型建立模型- -模型正演模型正演- -测量时间差测量时间差- -修改模型修改模型这一迭代过程一直进行到时差最小(最小平方差)这一迭代过程一直进行到时差最小(最小平方差)此时的速度就是要求取的速度,这个过程也是一个反演过程此时的速度就是要求取的速度,这个过程也是一个反演过程层析成像技术中既有正演也有反演层析成像技术中既有正演也有反演东方地球物理公司2007年5月地震旅行时层析成像的应用地震旅行时层析成像的应用 利用直达波和折射波的信息确定近地表速度结构,深度利用直达波和折射波的信息确定近地表速度结构,深度建立建立 速度模型速度模型 叠前深度偏移建立速度场叠前深度偏移建立速度场
23、 井间勘测,用地震旅行时井间层析成像估算井与井之间井间勘测,用地震旅行时井间层析成像估算井与井之间的地层的速度结构的地层的速度结构 垂直地震剖面法(垂直地震剖面法(VSP),地面到井中的层析成像),地面到井中的层析成像 天然地震学中也用广泛的应用天然地震学中也用广泛的应用 东方地球物理公司2007年5月第第18章章 3D反射地震反射地震解释过程建立解释过程建立3D空间概念的空间概念的重要性重要性东方地球物理公司2007年5月两个绕射点的绕射双曲线 (引自Wu 等人, 1996 ) ,(a) 和 (b)分别代表y =400m平面和 y = 800m 平面东方地球物理公司2007年5月对两个点绕射
24、模型2D偏移(上面)和3D偏移(下面)的对比 ( 引自 Wu等人,1996 ) 东方地球物理公司2007年5月未偏移、2D偏移、和3D偏移剖面的比较 说明需要3D成像(引自Brown, 1991 )东方地球物理公司2007年5月第第19章章介绍介绍AVO方法方法什么是什么是AVO?研究研究CMP道集内相对振幅,称作振幅随炮检距道集内相对振幅,称作振幅随炮检距变化的分析(变化的分析(AVO)。研究相对振幅随反射角的变化,这种方法称作)。研究相对振幅随反射角的变化,这种方法称作振幅随入射角变化的分析(振幅随入射角变化的分析(AVA) AVO分析能解决什么地质分析能解决什么地质问题问题 碎屑岩气藏直
25、接烃类指示碎屑岩气藏直接烃类指示 在碳酸盐油藏中可能识别孔隙发育带在碳酸盐油藏中可能识别孔隙发育带 随炮检距变化的反射系数随炮检距变化的反射系数 AVO的计算的计算东方地球物理公司2007年5月AVO、岩石物性和孔隙流体、岩石物性和孔隙流体弹性模量弹性模量与岩石性质的相关性优于速度和岩石性质的相关性与岩石性质的相关性优于速度和岩石性质的相关性弹性模量与地震弹性模量与地震- -波速度有关,由以下方程式表达波速度有关,由以下方程式表达式中的弹性模量:式中的弹性模量:k k是体积模量是体积模量 、 是拉每常数(表示不可压缩性)是拉每常数(表示不可压缩性), 是是剪切模量(表示刚性)剪切模量(表示刚性
26、) ,和,和 是密度是密度泊松比是泊松比是随随P-P-波和波和S-S-波速度比函数变化而变化的一个波速度比函数变化而变化的一个重要的重要的岩性参数岩性参数 东方地球物理公司2007年5月泊松比泊松比-东方地球物理公司2007年5月碎屑岩的碎屑岩的AVOAVO原理原理 当当P-P-波以非法线入射到达岩石界面时,入射波以非法线入射到达岩石界面时,入射P-P-波能量的一小部份转换波能量的一小部份转换为为S-S-波的能量。波的能量。P-P-波和波和S-S-波速度的差别将会因含气砂岩波速度的差别将会因含气砂岩/ /页岩接触和页岩接触和含气砂岩含气砂岩/ /含水砂岩接触引起不同的反射响应含水砂岩接触引起不
27、同的反射响应 泊松比是随泊松比是随P-P-波和波和S-S-波速度比函数变化而变化的一个重要的岩性参数波速度比函数变化而变化的一个重要的岩性参数AVOAVO依赖于炮检距的作用,它是一个附加的岩石物性的函数依赖于炮检距的作用,它是一个附加的岩石物性的函数 在在碳酸盐碳酸盐油藏中,没有证实油藏中,没有证实AVOAVO本身可识别孔隙流体,但可能识别孔本身可识别孔隙流体,但可能识别孔隙发育带。隙发育带。充满气的多孔石灰岩相与含水的石灰岩相相比可能具有较充满气的多孔石灰岩相与含水的石灰岩相相比可能具有较低的低的VpVp/Vs/Vs值值 对碳酸盐油藏对碳酸盐油藏AVOAVO的应用不能凭直觉进行的,其准确性取
28、决于优质的的应用不能凭直觉进行的,其准确性取决于优质的井控制和该地区计算的井控制和该地区计算的S-S-波速度。波速度。 东方地球物理公司2007年5月P-P-波和波和S-S-波对孔隙流体有不同的敏感度波对孔隙流体有不同的敏感度 少量的气体注入碎屑沉积岩孔隙空间中,岩石的少量的气体注入碎屑沉积岩孔隙空间中,岩石的P-P-波速度将剧烈减小波速度将剧烈减小S-S-波速度可以随气体进入孔隙空间而略有增加波速度可以随气体进入孔隙空间而略有增加 (a)(a)多孔固体注入少量气多孔固体注入少量气P-P-波速度迅速减小。波速度迅速减小。S-S-波速度随含气波速度随含气饱和度增加呈线性增加关系。(饱和度增加呈线
29、性增加关系。(b b)泊松比随含气饱和度增加呈)泊松比随含气饱和度增加呈现减小的效应现减小的效应(Allen (Allen 和和PeddyPeddy, 1993 ), 1993 )东方地球物理公司2007年5月 非法线入射时界面上波的分离非法线入射时界面上波的分离P波倾斜入射交界面上波型转换波倾斜入射交界面上波型转换入射的入射的P-P-波遇到两种介质的交界面时,波遇到两种介质的交界面时,分离成分离成4 4个分量:反射个分量:反射P-P-波、反射波、反射S-S-波、透射波、透射P-P-波和透射波和透射S-S-波波 分界面上入射角、反射角和透射角分界面上入射角、反射角和透射角 符合符合Snell定
30、理的关系定理的关系 反射和透射系数随入射角变化(随炮检距反射和透射系数随入射角变化(随炮检距)是是AVO分析的基础分析的基础 东方地球物理公司2007年5月AVO和和Zoeplpritz方程方程AVO分析一般考虑若干形式的分析一般考虑若干形式的Zoeppritz方程方程 对一个入射到交界面的平面波,对一个入射到交界面的平面波, Zoeppritz方程描述反射和透射的方程描述反射和透射的P- -波和波和S- -波波 许多研究推导出一些近似的Zoeppritz方程 Aki 和 Richard (1980 );Shuey (1985 );Hilterman (1990 );Smith 和 Gidlo
31、w (1987 ) ;Fatti 等人 (1994)简化了反射系数和入射角之间的关系式 Shuey (1985 )近似关系式是常用的近似关系式是常用的 式中式中R0是法线入射是法线入射P- P-波反射系数,或波反射系数,或“ “截距截距” ”,G是是“ “梯度梯度” ”项。梯度的定义项。梯度的定义是入射角的是入射角的函数,函数,CDP道集上在每个时间样点振幅的变化率。梯度包含了全道集上在每个时间样点振幅的变化率。梯度包含了全部的部的AVO效应效应 东方地球物理公司2007年5月部分其它部分其它Zoeplpritz方程方程HiltermanHilterman ( 1989 ) ( 1989 )提
32、出了提出了ShueyShuey (1985 ) (1985 ) 关系式的近似方程式关系式的近似方程式 式中式中NI NI = = 法线入射反射系数,法线入射反射系数, = = 上下介质间的泊松比之差,上下介质间的泊松比之差, = = 入射角入射角P-P-波和波和S-S-波的通用关系式波的通用关系式用用P-P-波和波和 S-S-波反射系数表示的流体因子波反射系数表示的流体因子 东方地球物理公司2007年5月AVO分析的陷阱产生的因素分析的陷阱产生的因素客观条件客观条件人为因素人为因素地质因素地质因素1、地表不一致性2、环境噪声3、震源噪声4、球面扩散5、非弹性衰减及子波弥散6、多次波7、各向异性
33、1、增益控制2、检波器组合3、最大炮检距太小4、记录道的不一致性5、设备产生的高低频干扰6、处理程序或参数不当7、模型参数不准1、地层倾斜或弯曲2、浅层气3、火成岩屏蔽4、盐丘遮挡5、岩性组合尽管如此,尽管如此,AVO在检测气藏,碳酸盐岩勘探中还是值得应用的技术在检测气藏,碳酸盐岩勘探中还是值得应用的技术东方地球物理公司2007年5月第第20章章油油 藏藏 表表 征征在在2121世纪,全世界多于世纪,全世界多于95%的石油将来自现有的油田的石油将来自现有的油田 通过有效的油藏表征增加产量是可能的通过有效的油藏表征增加产量是可能的油藏表征描述油藏和含有烃类的岩石特性油藏表征描述油藏和含有烃类的岩
34、石特性油藏表征是多学科领域的油藏表征是多学科领域的,油藏表征依赖于油藏工程、地质和地球物,油藏表征依赖于油藏工程、地质和地球物理技术。综合研究来自这些领域有关油藏的各个方面的信息理技术。综合研究来自这些领域有关油藏的各个方面的信息油藏表征项目组的构成需要工程师、地质家和地球物理家的相互协作油藏表征项目组的构成需要工程师、地质家和地球物理家的相互协作 油藏表征的基本目标油藏表征的基本目标 烃类的存在烃类的存在 储层的孔隙度储层的孔隙度 储层的渗透率储层的渗透率 东方地球物理公司2007年5月岩石物理岩石物理什么是岩石物理学什么是岩石物理学 “用岩石、测井资料研究岩石物理性质它们的相互关系的科学”
35、。(Robert E. Sherff,应用地球物理百科辞典孔隙度、渗透率、密度、地层的有效压力及地震波速度等孔隙度、渗透率、密度、地层的有效压力及地震波速度等岩石特性和它们之间的关系岩石特性和它们之间的关系是我们在储层地球物理领域关注的孔隙度是一个关键因素孔隙度是一个关键因素它将决定油藏在岩石中的补给。地震速度和孔隙度有关。通过应用穿过岩石骨架和充满流体孔隙的地震旅行时间平均值,Wyllie时间-平均方程使速度与孔隙度联系起来。Wyllie方程式对砂岩相当有效 最重要的物性参数渗透率最重要的物性参数渗透率,它与孔隙的连通性有关,它与孔隙的连通性有关,渗透渗透率与孔隙度通常有粗略的比例关系,有时
36、可根据地震振幅率与孔隙度通常有粗略的比例关系,有时可根据地震振幅和速度推断和速度推断另外,裂缝和裂缝的方向性影响渗透率 。东方地球物理公司2007年5月Wyllie时间时间- -平均方程平均方程(地震振幅解释)速度速度- -孔隙度转换方程孔隙度转换方程V V = =整个岩石的速度,整个岩石的速度, VmaVma = =基质的速度,基质的速度,VflVfl= = 孔隙流体的速度孔隙流体的速度 = = 孔隙度孔隙度 t 表示相应的旅行时间表示相应的旅行时间全部传播时间是分别通过全部传播时间是分别通过多孔隙物质和基质时间之和多孔隙物质和基质时间之和 东方地球物理公司2007年5月综合的油藏表征综合的
37、油藏表征 刻度问题刻度问题油藏采样数据类型的比较油藏采样数据类型的比较地面观测的波长是很粗的地面观测的波长是很粗的 、井间数据的分辨率按数量级井间数据的分辨率按数量级远超过地面地震数据远超过地面地震数据 、测井曲线和岩芯是高、测井曲线和岩芯是高- -频的采样,频的采样,但井与井之间的采样是低频的但井与井之间的采样是低频的 所有的油藏信息都是有用的。我们需要综合所有的资料所有的油藏信息都是有用的。我们需要综合所有的资料来描述油藏来描述油藏数据类型数据类型 油藏体油藏体 调查的波长调查的波长 地面3D地震数据 整个油藏数据体 10-100m ( 取决于速度和频率 )井中地震数据 井和震源之间的范围
38、(体积) VSP 5-50m 井间地震井间地震1-10m井的数据 井筒附近的岩体 10 cm 1 m东方地球物理公司2007年5月地面地震数据和井间成像(速度和反射系数)、井间速度、声波测井、地面地震数据和井间成像(速度和反射系数)、井间速度、声波测井、测井和岩芯的波长刻度比较测井和岩芯的波长刻度比较(引自引自Harris等人等人, 1995 ) 东方地球物理公司2007年5月井间反射与地面地震成像的比较井间反射与地面地震成像的比较 (a) 与地面地震反射测线对应的井间反射成像 b)西德克萨斯北Cowden油田的地面地震剖面,阴影区标出对应的井间勘测范围(引自Lines 等人, 1995 )
39、东方地球物理公司2007年5月第第21章章时间时间- -推移地震推移地震时间时间- -推移地震(有时称作推移地震(有时称作4D地震地震)最主要最普遍的应用就是向重)最主要最普遍的应用就是向重- -油油藏注蒸气的监测油油藏注蒸气的监测 当含油砂岩受热,温度从当含油砂岩受热,温度从250 C 上升至上升至1400C时,时,P-波速度明显降低波速度明显降低(大约(大约30%)。)。地震响应发生引人注目的变化地震响应发生引人注目的变化 地震速度图。通过画白圆圈的地震速度图。通过画白圆圈的几口井注入蒸汽导致地震速度几口井注入蒸汽导致地震速度减小减小(引自引自Lines等人等人, 1990 )。速度单位是
40、速度单位是km/s,注蒸汽前速,注蒸汽前速度为度为2.40 km/s 东方地球物理公司2007年5月蒸汽带的地震监测蒸汽带的地震监测一个实例一个实例加拿大Saskatchewan省的Pikes Peak重-油油田 。用蒸汽-驱提高采收率生产了4200万桶原油 。由于油藏内注入高温高压蒸汽,油的黏度降低,流(动)度增加。即可从临近的井筒也可以在用于循环注汽的同一口井中生产原油早白垩世Waseca组是产油层,产油深度约450m。是一个填满河口底部切割河谷的均匀砂岩组、砂泥岩组互层和一个页岩组盖层1991年勘测使用可控震源14-110 Hz扫描,2000年勘测使用14-150 Hz带宽扫描(三分量、
41、四分量)。要求第二次勘测处理用14-110 Hz的高截频滤波针对地震探测蒸汽前缘试验了5种技术: 监测勘测与基础勘测反射系数的差异反射系数的差异 监测勘测与基础勘测声阻抗估算的差异声阻抗估算的差异 监测勘测与基础勘测P-波旅行时时间间隔波旅行时时间间隔的比较 从多分量数据中估算Vp/Vs的变化 Q值(衰减系数的倒数)随温度的变化 东方地球物理公司2007年5月P- -波反射系数剖面对比波反射系数剖面对比 1991测线测线 H2000测线测线 东方地球物理公司2007年5月地震反射系数地震反射系数- -差异剖面差异剖面在圆圈内见到一个大的差异。这个范围处在注汽在圆圈内见到一个大的差异。这个范围处
42、在注汽/生产井位置,生产井位置,我们感兴趣的地段。差异是与注蒸汽有关联时间下凹的结果我们感兴趣的地段。差异是与注蒸汽有关联时间下凹的结果 (引自引自Watson等人等人, 2002 ) 反射系数差异是子波处理后获得的。对反射系数差异是子波处理后获得的。对2000测线应用子波整形让它与测线应用子波整形让它与1991测线匹配测线匹配 东方地球物理公司2007年5月声阻抗差异剖面声阻抗差异剖面圆形标出了阻抗减小最明显的范围。在圆形标出了阻抗减小最明显的范围。在Waseca- Sparky层层之间差异最大。在注汽井之间差异最大。在注汽井/生产井区域没有出现较低的阻抗生产井区域没有出现较低的阻抗 (引自
43、引自Watson等人等人, 2002 ) 东方地球物理公司2007年5月P-波旅行时比波旅行时比 Waseca 层层H2000/H1991线线P-波旅行时比波旅行时比(引自引自Watson等人等人, 2002 ) 东方地球物理公司2007年5月岩芯样品,温度对纵波和横波速度岩芯样品,温度对纵波和横波速度的影响(来源:岩芯实验室)的影响(来源:岩芯实验室)注气后纵波比横波速度下降更快,注气后纵波比横波速度下降更快,所以所以Vp/Vs比下降更快比下降更快Mannville- Lower Mannville层的层的Vp/Vs比绘图比绘图(引自引自Watson等人等人, 2002 ) 东方地球物理公司
44、2007年5月结结 论论时移地震是油藏表征的主要工具之一时移地震是油藏表征的主要工具之一时移地震主要的目标是显示由于油藏物理变化引起的地震时移地震主要的目标是显示由于油藏物理变化引起的地震响应差异响应差异 通过这些地震响应差异圈定油藏的动态范围通过这些地震响应差异圈定油藏的动态范围主要用于是蒸汽带的绘图主要用于是蒸汽带的绘图也许在也许在COCO2 2 和注水监测中还可应用和注水监测中还可应用 东方地球物理公司2007年5月第第22章章多分量地震多分量地震弹性波理论对地下特性更完全的描述需要考虑弹性波理论对地下特性更完全的描述需要考虑P-波和波和S-波波 P-波是纵向可压缩波,质点位移与波的传播
45、是同一方向波是纵向可压缩波,质点位移与波的传播是同一方向 S-波是横波,因为质点的位移是垂直波的运动方向波是横波,因为质点的位移是垂直波的运动方向 多分量地震采用垂直的、径向的和横向的三分量接收装置多分量地震采用垂直的、径向的和横向的三分量接收装置 为什么应用为什么应用S-波?波? 碎屑岩中要有效地区分砂岩和页岩,需要碎屑岩中要有效地区分砂岩和页岩,需要Vp/Vs比比 可以用可以用S-波调查裂缝的方位和方向波调查裂缝的方位和方向 描述横波在遇到垂直裂缝时分裂成快、描述横波在遇到垂直裂缝时分裂成快、慢波方式慢波方式(引自引自Sheriff, 1991 ) 快横波具有与裂缝平行的位移,快横波具有与
46、裂缝平行的位移,慢波具有与裂缝垂直的位移。慢波具有与裂缝垂直的位移。裂裂缝的密度和缝的密度和方向方向将决定将决定S-波的方波的方位以及快慢波之间的旅行时差。位以及快慢波之间的旅行时差。所以,这个方法能帮助所以,这个方法能帮助刻画储层刻画储层中的裂缝和描述裂缝的渗透性中的裂缝和描述裂缝的渗透性。 东方地球物理公司2007年5月转换转换- -波勘探波勘探 横波方式的主要障碍是使用横横波方式的主要障碍是使用横- -波震源的价格和对地面的破坏。因此,波震源的价格和对地面的破坏。因此,大部分的应用是使用转换波震源多分量记录大部分的应用是使用转换波震源多分量记录 P-S转换波的描述(引自Stewart 等
47、人,2002) 东方地球物理公司2007年5月转换转换- -波波 的用途的用途:用于岩性识别或流体检测;消除气:用于岩性识别或流体检测;消除气- -云的勘探云的勘探 Tommeliten油田海-底记录P-波叠加(上)和P-S波叠加(下)的比较(引自Granli等人, 1999 ) 东方地球物理公司2007年5月第第2323章章垂直地震剖面法垂直地震剖面法什么是垂直地震剖面法垂直地震剖面法(VSP)?垂直地震剖面法垂直地震剖面法(VSP)是在是在地面激发,测量不同深度推靠在钻井井壁的检波器记录的地震信号地面激发,测量不同深度推靠在钻井井壁的检波器记录的地震信号” ” VSP的用途的用途提供地震与
48、井中之间的连接;形成详细的速度剖面;提供地震与井中之间的连接;形成详细的速度剖面;预测钻头下面地层深度;精确地估算各向异性参数预测钻头下面地层深度;精确地估算各向异性参数VSP的类型:常规的类型:常规VSP放炮分成两类放炮分成两类 : 零井源距零井源距VSP:包括:包括: 零井源距零井源距VSP-相距很近的接收点并记录全部波场相距很近的接收点并记录全部波场 激发点位于井口的校验炮激发点位于井口的校验炮- -接受点稀梳地布置在井中并且仅使用初至接受点稀梳地布置在井中并且仅使用初至 非零井源距非零井源距VSP:激发点距离井在有效间隔内激发点距离井在有效间隔内 ,包括:多井源距的包括:多井源距的 V
49、SP。多方位。多方位VSP。3DVSP既是多井源距的也是多方位的既是多井源距的也是多方位的 非常规的非常规的逆逆VSP是通过地面是通过地面接收井中激发获得的接收井中激发获得的 东方地球物理公司2007年5月VSP 示意图(示意图(a)零井源距)零井源距 、近井源距或校验炮,(、近井源距或校验炮,(b)非零井源距)非零井源距 东方地球物理公司2007年5月逆逆VSP采样射线路径示意图采样射线路径示意图 (a)2D勘测,(勘测,(b)3D勘测勘测 东方地球物理公司2007年5月井间勘测井间勘测使用两口井,在一口井中接收,而在另一口井中激发。能量直接通过地下传播,用获得的地下速度剖面转换为初至波旅行
50、时 井间勘测示意图井间勘测示意图 东方地球物理公司2007年5月VSP 采集和处理采集和处理三分量检波器三分量检波器 。许多检波器放置在一个单独的许多检波器放置在一个单独的“ “用具用具” ”内内 (例如有(例如有5 5个检波器的个检波器的5-级用具)。在用具内检波器间距可以变化检波器间距和级数将级用具)。在用具内检波器间距可以变化检波器间距和级数将根据勘测如何实施决定。已开发了根据勘测如何实施决定。已开发了80级,长度级,长度3048 m和和400级,长度级,长度7630 m 的用具的用具 当当VSP勘测要求勘测要求VSP用具内检波器间距用具内检波器间距相等时,在井中相等时,在井中VSP用具
51、每次上提一个用具每次上提一个用具长度。用具长度。( a ) 和和 (b) 图解说明了连续图解说明了连续的两炮,每炮用具提升一个用具长度。的两炮,每炮用具提升一个用具长度。当当VSP用具内检波器间距大于用具内检波器间距大于VSP勘测勘测要求,用具在井中交叉接收点位置提升,要求,用具在井中交叉接收点位置提升,(c )和和(d) 图解说明了对需要的连续两炮图解说明了对需要的连续两炮用具在井中交叉提升用具在井中交叉提升 东方地球物理公司2007年5月VSP 采集和处理采集和处理VSP震源震源:VSP地面震源和地面地震勘测震源是一样的地面震源和地面地震勘测震源是一样的 重垂、重垂重垂、重垂震源、弹性震源
52、、弹性- -波发生器(波发生器(EWG)、炸药和可控震源。震源的类型由勘测施)、炸药和可控震源。震源的类型由勘测施工决定工决定 ; 井下震源包括压电震源、可控震源和气枪,井下震源包括压电震源、可控震源和气枪, 要没有破坏性。压电震源和气枪震源都有要没有破坏性。压电震源和气枪震源都有 套管波噪音问题套管波噪音问题VSP处理处理:了解给出的显示很重要了解给出的显示很重要 原始原始VSP炮记录炮记录,标出了,标出了上行波和下行波能量上行波和下行波能量 VSP资料沿资料沿x-轴代表检波点深度轴代表检波点深度,沿沿z-轴代表旅行时轴代表旅行时东方地球物理公司2007年5月 零零- -井源距处理的第一步是
53、将波场分离成上行波和下行波井源距处理的第一步是将波场分离成上行波和下行波(a)全波场记录,分成)全波场记录,分成 (b)下行波能量和)下行波能量和 (c) 上行波能量。上行波能量。(c)上面的箭头指出了某些主要的阻抗界面上面的箭头指出了某些主要的阻抗界面(引自引自Hinds 和和Kuzmiski, 2000) 初至曲线(或初至曲线(或一次下行波能一次下行波能量)给出了井量)给出了井中详细的垂直中详细的垂直-速度剖面速度剖面 上行一次波能量上行一次波能量给出了主要阻抗给出了主要阻抗界面的深度。位界面的深度。位于阻抗边界的检于阻抗边界的检波点同时记录了波点同时记录了下行波能量的初下行波能量的初至和
54、上行波能量至和上行波能量 波场分离后很容波场分离后很容易从一次波能量易从一次波能量中识别多次波。中识别多次波。初至能量仅仅是初至能量仅仅是下行波一次波能下行波一次波能量。所有的其它量。所有的其它下行波场是多次下行波场是多次波波 东方地球物理公司2007年5月整个记录减去直达波时间会把下行波拉直,可用于估算反可用于估算反褶积算子来提高褶积算子来提高VSP成像质量成像质量 (a)下行波场下行波场 (b) TT显示显示下行波拉直下行波拉直 (引自引自Hinds 和和Kuzmiki, 2000) 东方地球物理公司2007年5月 整个记录加直达波时间会把上行波拉直,突出上行波整个记录加直达波时间会把上行
55、波拉直,突出上行波(a)上行波场上行波场 (b) +TT显示显示校直上行波。箭头指出转换校直上行波。箭头指出转换 前、后的两个层前、后的两个层(引自引自Hinds 和和Kuzmiki, 2000)。多次波不与多次波不与初至能量交初至能量交叉,很容易叉,很容易将多次波与将多次波与一次波区分一次波区分 东方地球物理公司2007年5月走廊切除、走廊叠加走廊切除、走廊叠加 (a)应用反褶积的应用反褶积的+TT显显示示(b)走廊切除走廊切除 (c)内部走廊叠加内部走廊叠加 (d)外部走廊叠加外部走廊叠加 (引自引自Hinds 和和 Kuzmiki, 2000) 东方地球物理公司2007年5月第第24章章
56、地球物理数据的综合反演地球物理数据的综合反演 综合反演综合反演两种或更多类型独立观测值的同时反演单一类型数据的反演往往是不确定性单一类型数据的反演往往是不确定性(多解多解),通过综合通过综合(联合)反演不同类型地球物理观测资料,可以减少不确(联合)反演不同类型地球物理观测资料,可以减少不确定性定性 “ “综合反演综合反演” ”的目的是要获得与所有可利用的地面和井中的目的是要获得与所有可利用的地面和井中地球物理数据一致的模型地球物理数据一致的模型综合反演包括联合反演和连续反演综合反演包括联合反演和连续反演 联合反演联合反演两个独立的反演过程两个独立的反演过程 ,仅限于两种类型独立观测值的同时反演
57、 连续反演连续反演联合反演和连续反演步骤的差别在于它们对联合反演和连续反演步骤的差别在于它们对观测资料处置不同。联合反演加权观测的数据集,并把它观测资料处置不同。联合反演加权观测的数据集,并把它们放置到一个数据矢量中,然而们放置到一个数据矢量中,然而连续反演单独对待观测的连续反演单独对待观测的数据集数据集东方地球物理公司2007年5月地震和重力数据联合反演地震和重力数据联合反演 联合反演流程图联合反演流程图用地震和重力数据迭代模拟用地震和重力数据迭代模拟联合反演的一般方法联合反演的一般方法S =地震数据,地震数据,MS =地震模地震模型数据,型数据,G = 重力数据,重力数据,MG = 重力模
58、型数据,重力模型数据,i = 地层地层i的密度,的密度,Vi = 地层地层i的速度,的速度,di (x ) =地层地层i的厚度的厚度 (引自引自Lines等人等人., 1988 ) 东方地球物理公司2007年5月顺序反演的实际数据顺序反演的实际数据 用实际数据的例子来说明连续反演,包括以下数据集连续反演,包括以下数据集:地震地震CDP叠加剖面叠加剖面井中声波测井数据井中声波测井数据垂直地震剖面垂直地震剖面(VSPs)井中重力仪井中重力仪(BHGM )数据数据地面布伽重力图地面布伽重力图所有的数据反演之中,通常连续反演是首选所有的数据反演之中,通常连续反演是首选 东方地球物理公司2007年5月地
59、震数据集应用连续地震数据集应用连续- -反演反演 使用可获得的地震和声波数据估算的速度-深度模型 连续-反演问题步骤1地震旅行时数据的反演 东方地球物理公司2007年5月地震旅行时反演流程图地震旅行时反演流程图(引自引自Lines等人等人, 1988 ) 东方地球物理公司2007年5月用地震旅行时反演获得的初始地震模型。速度单位是用地震旅行时反演获得的初始地震模型。速度单位是m/s。 井的位置在井的位置在2700-m处处(引自引自Lines等人等人, 1988 ) 东方地球物理公司2007年5月连续反演类的重力反演流程图连续反演类的重力反演流程图(引自引自Lines等人等人, 1988 ) 东
60、方地球物理公司2007年5月结结 论论完全自动的综合反演将是很困难完全自动的综合反演将是很困难 ,在所有的数据反演之中,通常连续,在所有的数据反演之中,通常连续反演是首选反演是首选 反演过程中数据集之间的联结是由层的几何形态产生的反演过程中数据集之间的联结是由层的几何形态产生的 重力数据能够使我们估算那些层的边界,它不如地震数据确定的好。重力数据能够使我们估算那些层的边界,它不如地震数据确定的好。然而,地震数据确定模型的层速度,地面和井中重力数据集确定层的然而,地震数据确定模型的层速度,地面和井中重力数据集确定层的密度密度重力数据模型拟合中,约束的最小重力数据模型拟合中,约束的最小- -平方算
61、法是有用的平方算法是有用的 经验的速度经验的速度- -密度关系式在一些地质区域证实是有成效的密度关系式在一些地质区域证实是有成效的 合作反演的最终目标是产生一个最终的地下模型,这个模型满足所有合作反演的最终目标是产生一个最终的地下模型,这个模型满足所有可得到的地球物理资料可得到的地球物理资料几种地球物理数据集的反演近似地描绘了地下的同一部分,几种地球物理数据集的反演近似地描绘了地下的同一部分,比任何单一数据集的反演减小了不确定性比任何单一数据集的反演减小了不确定性 东方地球物理公司2007年5月第第25章章地质统计学地质统计学地质统计学是模型预测的工具并能指出预测的概率地质统计学是模型预测的工
62、具并能指出预测的概率 单单- -变量型的空间预测变量型的空间预测 克里金(克里金(Kriging)地质统计法)地质统计法 这种方法是在一个单点上对数据点周围数据用加权和来估算变量值 假如估算一个变量的值为u, 在图上点的特定的位置x0,可以用下面表达式表示变量: x0 和xi 表示2D图上定义的向量。变量u可以表示某些物理量(如孔隙度或地层深度) 在一个特定位置的模型估算值是用一些权wi值对相邻点进行加权给出的。计算这些权值,使得在这些点上模型值和观测值之间的预测误差最小化。在这个极小化中,用最小-平方准则。要使估算值不偏离,可以增加约束,这个约束是所有权值之和等于1 东方地球物理公司2007
63、年5月 多变量预测多变量预测协克里金法协克里金法 数值数值 由由w 和和W分别加权的和分别加权的和 来描写。数学表达如下来描写。数学表达如下 新的变量新的变量 (xj )可以表示地震振幅可以表示地震振幅 某种意义上,从有足够变量的协克里金法估算的模型可以使用全部来某种意义上,从有足够变量的协克里金法估算的模型可以使用全部来 自油藏数据测量值。这个公式可以达到前面油藏表征所规定的目标,自油藏数据测量值。这个公式可以达到前面油藏表征所规定的目标, 并用所有可利用的数据定义综合反演模型并用所有可利用的数据定义综合反演模型东方地球物理公司2007年5月 用克里金法得到的砂岩用克里金法得到的砂岩-河道深
64、度图河道深度图坐标轴代表距离(坐标轴代表距离(m),误差误差10m (引自引自 Gosse, 1999 ) 用协克里金法对海绿石构造深度估算值用协克里金法对海绿石构造深度估算值 误差误差2.5m(引自引自 Gosse, 1999 ) 东方地球物理公司2007年5月第第26章章绘等值线的艺术和科学绘等值线的艺术和科学石油勘探工作由一系列描绘地下地质情况的等值线图随伴 地下构造图 上面的地层被剥后的地形 基准面每种图件都有一个起始深度的基准面最初拾取是通过测井曲线与层的顶部连接。然后通过地震数据体对比反射波至逐道拾取时间或深度鉴别地震拾取质量是好、无效或是不好。质量通常是信噪比的函数并影响某个点的
65、可信度 等厚图表示地层厚度的变化。等时线图表示层间地震-反射时间的变化绘等值线应该把一些的地质概念放在心上 绘等值线 解释的过程东方地球物理公司2007年5月绘等值线的简单的准则绘等值线的简单的准则 在数据点采集时,要孤立最高和最低的值。一条等值线应该在点与点之间通过,那些点的数值比等值线上的数值要高或低。一条等值线决不能穿过它本身或其它的等值线,但逆断层和倒转褶皱除外。一条等值线不能并入不同值的等值线。对山谷和河流等值线总是尖头指向上游。当控制点不足时用虚线。尽管是可靠点,要保持图件简明、可行。东方地球物理公司2007年5月关于机械绘等值线的警告关于机械绘等值线的警告计算机程序对简单的地质构
66、造能够有效地绘制等值线,但是计算机算法不包含地质概念,并且没有手工绘制的断层的不连续性好。计算机算法对噪音拾取过于敏感,使高频等值线摆动,在地质上是不合的。(a)计算机绘制的等值线和几个内在的问题计算机绘制的等值线和几个内在的问题 (b)手工绘制的等值线图(引自手工绘制的等值线图(引自Tucker, 1988)东方地球物理公司2007年5月机械绘等值线绘图中的缺点机械绘等值线绘图中的缺点以为等值线值是以为等值线值是“ “真理真理” ”(也就是绝对正确);地震拾取(也就是绝对正确);地震拾取值有误差并且应该不要过于严格,使得每一条等值线围绕值有误差并且应该不要过于严格,使得每一条等值线围绕这个值
67、弯曲一下。这个值弯曲一下。使用的等值线间隔太小,因此造成图上突出了噪声,胜过使用的等值线间隔太小,因此造成图上突出了噪声,胜过地质趋势地质趋势绘制等值线与地质类型不相称;一个地区地质背景的知识绘制等值线与地质类型不相称;一个地区地质背景的知识将帮助避免这点将帮助避免这点仅仅绘一层图仅仅绘一层图箭尾型等值线,通常由于地震测线位置错误造成箭尾型等值线,通常由于地震测线位置错误造成产生方形等值线产生方形等值线用太多的断层,导致象一块用太多的断层,导致象一块“ “碎裂的玻璃碎裂的玻璃” ” 形状形状 东方地球物理公司2007年5月绘等值线的艺术和科学绘等值线的艺术和科学地震地震- -等值线图是勘探工作
68、非常重要的部分。它对发展石等值线图是勘探工作非常重要的部分。它对发展石油勘探是必不可少的油勘探是必不可少的 一个地球物理解释员绘等值线的艺术和科学是关键技巧一个地球物理解释员绘等值线的艺术和科学是关键技巧 图件应该清楚地描绘地质构造类型和井位的依据。图件应图件应该清楚地描绘地质构造类型和井位的依据。图件应该综合定位、地质、地球物理和油藏工程的观念该综合定位、地质、地球物理和油藏工程的观念 专业知识和实际经验对等值线绘图是十分重要的专业知识和实际经验对等值线绘图是十分重要的东方地球物理公司2007年5月两张图对比两张图对比 (a)有许多断层的有许多断层的“破碎的玻璃破碎的玻璃”外观外观 (b)描
69、绘勘探区域的图描绘勘探区域的图 (引自(引自Tucker, 1988)东方地球物理公司2007年5月第第27章章结结 束束 语语覆盖了若干基础的地球物理原理覆盖了若干基础的地球物理原理 解释是一个综合的过程解释是一个综合的过程 。有效地解释必须利用全部已知的资料,来提高解有效地解释必须利用全部已知的资料,来提高解释精度。释精度。 在过去的十年,地球物理解释发生了许多变化在过去的十年,地球物理解释发生了许多变化 3D地震勘探现在已是油藏表征的标准工具和常规手段地震勘探现在已是油藏表征的标准工具和常规手段 地震解释工作站使巨大数量的地震解释工作站使巨大数量的3D数据能够容易处理和有效地进行分析数据
70、能够容易处理和有效地进行分析 可视化工具已经得到重大的进展可视化工具已经得到重大的进展 通过相干技术的发展,断层探测已经取得了很大的进步。通过相干技术的发展,断层探测已经取得了很大的进步。地震解释的缺陷将始终伴随着我们地震解释的缺陷将始终伴随着我们 多年来,一个最普遍的缺陷来自对地震剖面时间异常当作深度异常的误解多年来,一个最普遍的缺陷来自对地震剖面时间异常当作深度异常的误解 各向异性有关的陷阱,例如,在各向异性介质中以各向同性方法成像,将各向异性有关的陷阱,例如,在各向异性介质中以各向同性方法成像,将 引起成像位置的错误引起成像位置的错误 东方地球物理公司2007年5月 对于开发石油和矿产勘探解释是最重要对于开发石油和矿产勘探解释是最重要 解释是艺术和科学二者的结合解释是艺术和科学二者的结合 解释需要熟练人员的干预和洞察力解释需要熟练人员的干预和洞察力 使用先进的计算机硬件和软件使用先进的计算机硬件和软件