各向异性深度偏移
VTI各向异性叠前深度偏移技术应用研究——以JZ工区三维地震资料处理为例
狄 帮让 ， 培成． 震偏 移成 像分 辨率 的定 量分 析 顾 地 ＿］ 石 油大学学报 ，０ ５ ２ （ ） ２—６ Ｊ． ２ ０ ，９ ５ ：３２
邱天林． 沙漠地 区地震 资料处 理方法 研究 ［］ 石 油 Ｊ．
物探 ，０ ４ ４ （ ）４ ８４ １ ２ ０ ，３ ５ ：５—６
１ — ９７— ２０１
［３ 王丽， ］ 屠世杰 ， 赵传雪 ， 叠前深度偏 移技术在 江苏 等．
火成岩地区的应用 口］ 石油物探 ，０６ ４ （） ７—２ ． ２ ０ ，５ １ ：９８
［ 孙耀华 ， ４］ 宋建平． 叠前深度偏移技术在焉耆盆地本东 地 区构 造研 究中的应用Ｅ］ 石油物 探 ，０ ３ ４ （ ） Ｊ． ２０ ， ２ ３ ：
（ 接 第 １ ６页） 上 ８
［４ Ｓｈｏｄｒ 著． ］ ｃｒｅｅ Ｗ Ｆ 覆盖次数和面元尺寸怎样影 响数 据 可解释 ｌ［］高林译． 生 Ｊ． 石油物探译丛 ，９９１ ：９３ １９（）１—２ ［５ 张伟 ， ］ 尹成 ， 田继东 ， 三维观测 系统参数 的退化性 等． 处理试验 ［］ 石油物探 ，０ ７４ （ ）６ —３ Ｊ． ２ ０ ，６ １ ：９７ ［６ 张军华 ， 帆 ， ］ 张 郑旭 刚 ， 地 震采 集 脚 印综 合评 述 等． Ｉ］ 石油仪器 ，０ ７ ２ （） Ｉ４ －． Ｊ ２ ０ ，１ ５ ：－ ［７ 侯 成福 ， 连斌 ， ］ 蒋 高书琴． 维观测 系统与 采集脚 印 三 ＿］ 石油地球物理勘探 ，０ ７４ （ ） ６ １６ ５ Ｊ． ２ ０ ，２ ６ ：１ —１ ［ ８］ 张军华 ， 藏胜涛 ， ， 地震资料信 噪 比定量计算 徐辉 等． 及方法 比较 ［］ 石 油地 球 物理 勘 探 ， ０ ９ ４ （ ） Ｊ． ２０ ，４ ４ ：
VTI介质克希霍夫叠前深度偏移及其应用
VTI介质克希霍夫叠前深度偏移及其应用夏常亮;王永明;夏密丽;刘红久;胡浩;王祥春【摘要】针对叠前深度偏移速度反演多解性及层位标定和偏移结果不匹配,以伊拉克某构造复杂区块地震资料为例,详细介绍了垂直对称轴横向各向同性(VTI)介质的克希霍夫(Kirchhoff)叠前深度偏移及其应用和注意事项.提出利用剥层层速度修正方法反演层速度和测井曲线趋势约束联合解决速度反演多解性问题;利用叠前时间偏移均方根速度场通过约束速度反演(CVI)获得初始沿层层速度,从而保证初始层速度场的准确性和有效减少剥层层速度修正方法反演层速度的迭代次数;通过VTI介质的偏移解决偏移结果与层位标定不匹配问题.实际应用表明,前述Kirchhoff叠前深度偏移流程,能够有效提高叠前深度偏移工作效率,获得可靠性更强的深度域层速度模型,有效提高速度反演精度,获得与井上层位一致的地震层位,满足勘探开发的需求.【期刊名称】《新疆石油地质》【年(卷),期】2018(039)006【总页数】5页(P732-736)【关键词】Kirchhoff叠前深度偏移;横向各向同性介质;Thomsen参数;层剥离;层速度【作者】夏常亮;王永明;夏密丽;刘红久;胡浩;王祥春【作者单位】中国地质大学地球物理与信息技术学院,北京100083;中国石油东方地球物理勘探有限责任公司研究院长庆分院,西安710021;中国石油东方地球物理勘探有限责任公司研究院长庆分院,西安710021;中国石油东方地球物理勘探有限责任公司研究院长庆分院,西安710021;中国石油东方地球物理勘探有限责任公司研究院长庆分院,西安710021;中国石油东方地球物理勘探有限责任公司研究院长庆分院,西安710021;中国地质大学地球物理与信息技术学院,北京100083【正文语种】中文【中图分类】P631.4叠前深度偏移是解决复杂地质构造条件下地震波场成像的有效工具，偏移算法及其理论依据是决定偏移效果的关键[1]。
VSP_随钻地震地质导向建模技术应用研究
442023年6月上 第11期 总第407期信息技术与应用China Science & Technology Overview0 引言随着油气勘探开发的不断深入，石油地震勘探目标向尺度小、细、深及复杂特征的趋势发展，高质量勘探开发难度日益加大，若要资源突破，物探先行已成为目前高质量勘探、少井高产和效益开发的主要方法，因此，地震资料处理速度建模技术得到了迅速发展，已经从常规的时间域速度建模向深度域速度建模发展，并研发了地下非均质地层介质的各向异性速度建模，应用国内外很多建模学者正在研究的FWI 全波形反演速度建模。
当前，高精度速度建模还处于瓶颈期，超深层碳酸盐岩地质复杂，速度精度不足导致断裂、缝洞体等成像不清晰、构造归位不准确等问题，FWI 全波形反演速度建模成为其追求的目标，但是由于FWI 建模对资料要求很高，陆上地震超深层勘探目标效果还不显著，所以，要获得更多高品质地震资料，实现“多做物探少打井，打高产井”的目标，迫切需要效率和精度都高的速度建模技术来支撑，保证地震反演速度精度与钻井速度高度一致，再通过高精度成像结果为钻井提供更好的分析资料，实时指导钻井轨迹的调整。
井中地震勘探作为能够快速获得垂向地震剖面和最初了解地下信息的技术，随着采集装备、处理技术的发展，从1917年至今，井中地震勘探技术已经在垂直地震剖面（VSP）基础上形成了零井源距VSP、非零井源距VSP、变井源距VSP、井间地震、三维VSP、随钻地震等系列，成为不可或缺的勘探方法，在油气勘探开发中被广泛应用[1]。
其中，零井源距VSP、非零井源距VSP 被广泛应用于层位与深度标定、速度求取、地震波吸收衰减因子求取和提高分辨率井控处理；随钻导向技术逐步用于钻头前地层深度预测、地层压力预测、钻头导向和提高储层钻遇率方面，并与零井源距VSP 一起应用于高精度速度建模，通过获得的高精度速度对地震资料进行高精度成像，从而帮助随钻井调整轨迹方向。
各向异性叠前深度偏移技术识别低幅度构造
轴在观 测坐 标 系中具有 倾 角时就 会形成 倾斜 横 向各 向同性 先假设 ５＝ ，先进行初始基于 Ｖ Ｔ Ｉ 介质的 目 标线叠前深 介 质。垂直 对称 轴的横 向各 向同性介 质模型 是非 常重要 的 度 偏移 ，然后 进行 迭代 ，求 取准确 的 参数 。 同样 ，如 各 向异性 模型 ，因为 实际中约 ７ ０ ％ 的沉积岩 石展现 为垂直 果我 们有准 确的 ， ７ 值 和 值 ，也可利用 公式计算 ： 对称 轴 的横 向各向 同性 介 质的各 向异性 ，常用来 描述 由周 期性 的薄互 层 、岩石 内部结 构和 平行排 列 的微裂 隙 引起 的
气生 、储 、 盖 方 面 的 作 用 。
［ （ Ｖ ｏ ／ Ｖ ｏ ） 一 １ ］ ／ ２
对于 不同 的 值 ，Ｐ波 的相 速度 曲线在 近水 平时 相差
较大 ，说 明 值决定 了 Ｐ波 的水 平相速 度的大 小 ，反 映 了
横 向各 向同性 介质 的对称 轴为 铅垂 方 向 ，称其 为垂 直 Ｐ波水 平传 播 的相速 度与垂 直传 播的 相速 的差异 。地 震资 对称 轴 的横 向各 向 同性 （ ＶＴＩ介 质 ） ；横 向各 向 同性 介 料分 析 中的各 向异性 ，也 是基 于弱各 向异性 的假 设 ，纵波 质的对 称轴 为水 平时 ，称其 为水 平对称 轴的 横 向各 向同性 速 度取决 于参 数 和 ， 当 ＝ 时 ，简 化为 椭 圆各 向异 （ ＨＴ Ｉ介质 ） ，当垂直 对称 轴的横 向各 向同性 介质的对 称 性 ，所以 一般 我们 在获取 了各 向异性 参数 后 ，我们 可以
同时 ，我 们知 道对 于 各 向同性 速 度 和 各 向异性 速 度 有 Ｖ ｏ ／ Ｖ ａ 就等于 △ ／ △ Ｚ Ａ ，由此改公式也可以表示为 ：
深度偏移培训总结-高现俊
各向异性叠前深度偏移
各向同性叠前深度偏移速度 远道拉平
N
300近道是否拉 平 优化各向同性速度V0
Y
从井资料和各向同性资料计算Delta 把各向同性层速度V0比例转换为 各向异性速度Va Epsilon＝Delta
Y
N
优化 Epsilon
运行VTI各向异 性叠前深度偏移
结 束
运行VTI各向异性叠前深度偏移
9建立剩余深度延迟模型
10层析成像更新层速度模型
11 重复7、8步骤直到道集 拉平和剩余深度延迟归零 12创建最终深度域层速度体 13最终深度偏移
处理中心深度偏移学习班第7期
优化深度域层速度体
初始深度域层速度体
最终深度域层速度体
处理中心深度偏移学习班第7期
叠前深度偏移
初始深度偏移剖面
最终深度偏移剖面
处理中心深度偏移学习班第7期
学习体会 1、初学者正确理解深度域的概念对学会深度偏移 很重要，例如：面、图、模型、实体模型、体等。
2、丰富的构造解释知识对做好深度偏移很关键。
3、上手操作是牢记纷繁操作的根本。
处理中心深度偏移学习班第7期
感谢左老师的悉心指导！ 感谢处理中心为我们提供的学习机会！
感谢华北分院为我们所作的周到安排！
构造模型
弯曲射线时间偏移剖面
处理中心深度偏移学习班第7期
复杂构造叠前深度偏移流程：
1 建立工区 2 数据加载 3 时间偏移域层位解释 4 时间偏移域层位建模 5 抽取沿层层速度 ６编辑沿层层速度并建模 ７ 创建深度域层位模型 8 创建初始深度层速度剖面 9 叠前深度偏移
10 深度域构造模型调整并建模
处理中心深度偏移学习班第7期
各向异性叠前深度偏移
VTI介质角度域叠前深度偏移
VTI介质角度域叠前深度偏移李江;李庆春【摘要】研究了VTI介质角度域偏移方法,以各向同性双平方根方程的角度域偏移方法为基础,从VTI介质qP波频散关系出发,推导出VTI介质角度域偏移的波场延拓算子;在频率—波数域处理以横向均匀速度传播的波场,在空间域处理具有速度扰动特征的波场以提高波场延拓精度.模型试算和实际资料处理结果表明:各向同性偏移方法由于未考虑各向异性参数的影响,绕射波不能完全收敛,波场聚焦效果差,降低了成像剖面的分辨率和信噪比,不能对地质构造精确成像;VTI介质角度域偏移可对断层、盐丘、小尺度地质体精确成像.对于角度域偏移产生的角度域共成像点道集(ADCIG)而言,各向同性偏移的ADCIG同相轴无法校平,残留断点绕射波,波场无法正确聚焦,不能正确反映局部地质特征;VTI介质偏移的ADCIG同相轴较平直,角度范围更宽,波场归位准确,精度较高.因此,利用VTI介质角度域偏移方法可对复杂构造精确成像.%The angle domain migration in VTI media is discussed in this paper.Based on the angle domain migration of isotropic double-square-root equation,a wavefield extension operator of angle domain migration in VTI media is derived from the qP wave dispersion equation. The complex velocity and anisotropic parameter field are divided into two parts,one is lateral uniform background field,and the other is the disturbances of velocity and anisotropic parameters.Then the wave propagated at a uniform velocity and anisotropy field is processed in the frequency wavenumber domain,and the wave propagated at a velocity and anisotropy disturbance field is corrected by time-shift in the spatial domain.So the accuracy of wave field extension is greatly improved.Basedon our model and real data tests,the following observation are obtained:A.Because the isotropic migration method does not take into account the influence of anisotropy parameters, the diffraction wave cannot be completely converged and the wavefield is misfocused,which causes low resolution and low signal-to-noise ratio.So conventional migration methods cannot accurately image geological structures,while the VTI media migration can accurately image faults,salt mounds,and small scale geological bodies;B.For the angle domain common imaging gathers(ADCIGs)generated by the prestack migration, the event of the isotropic method cannot be equalized, the residual fault diffraction is wound and the wavefield cannot be properly focused,which cannot correctly reflect local geological characteristics,while the ADCIGs from VTI media migration are relatively straight,the angle range is wider,the wavefield is accurately positioned.So the image accuracy is higher. Therefore,the angle domain migration in VTI media can be suitable for complex-structure accurate imaging.【期刊名称】《石油地球物理勘探》【年(卷),期】2019(054)002【总页数】12页(P330-340,前插3)【关键词】VTI介质;角度域偏移;双平方根方程;频散关系;波场延拓【作者】李江;李庆春【作者单位】中国煤炭科工集团西安研究院有限公司,陕西西安 710077;长安大学地质工程与测绘学院,陕西西安 710054【正文语种】中文【中图分类】P6310 引言地球介质一般具有各向异性特征，其中VTI介质(具有垂直对称轴的横向各向同性介质)是一种常见的各向异性介质。
各向异性介质中的叠前高斯波束深度偏移
摘要 ：高斯波束深度偏 移 （ ｌ １９ ，０ １ Ｈｉ，９ ０ ２ ０ ）是 一种可替代 Ｋ ｒｈ ｏ 偏移和波动方程偏移的方法。它不但克服 了 ｌ ｉ ｈｆ ｃ
Ｋｉ ｈ ｏ 偏移不 能对 多次波成像的缺 陷， ｒ ｈｆ ｃ 同时又保 留 Ｋ ｒｈ ｏ 偏 移利用 回转波对陡倾角构造成像 的优点。它通过 ｉ ｈｆ ｃ 直接 向下外推局部平面波代替 单道散射波 ， 避免 了 Ｋｉ ｈ ｏ 偏移 中固有 的对 不稳定噪 音的偏移， ｒ ｈｆ ｃ 可获 得与波动方
Ｈｌ ２０ ） 用 公 式 阐 明 了 共 偏 移 距 域 叠 前 ｉ ｌ（０ １ Ｇ Ｍ， Ｂ 本文 的共 炮 点 Ｇ Ｍ 对他 的公 式 进 行 了修 改 。 Ｂ
一
成像问题 ， 在构造复杂地区成像效果差。另一方面 ， 波动方程偏移可 以对多次波正确成像 ， 但在陡倾角 构造 区成像受 限制 ， 同时 ， 因花费太 高也很难将波
动方 程偏 移延 伸 到各 向异性 介 质 。
高斯波束偏移 （ Ｂ 是近年来提 出的一种替 Ｇ Ｍ） 代 叠 前 Ｋｒ ｈｆ 和 波 动 方 程 偏 移 方 法 （ ｉ， ｉｈ ｏ ｃ Ｈｌ ｌ
１９ ，０ １。该方 法保 留了 Ｋｒｈ ｏ 偏 移 方法 的许 ９ ０２０ ） ｉ ｈｆ ｃ 多 优 点 ， 能 够有 效 地利 用 回转波 对 陡 倾 角构 造 成 如
程偏 移相 当的地下成像效果 。
关键词 ： 向异性介质 ； 各 叠前偏移成像 ； ｒｈ ｏ 偏移； Ｋｉ ｈ ｆ ｃ 高斯波束偏移
在 过去 １ 里 ，ｉ ｈ ｏ 偏移 一直 是 叠前 地震 Ｏ年 Ｋｒ ｈｆ ｃ 成 像 的主要 方法 。它可 以有 效地 利用 回转 波对 陡倾
１ 共炮点叠前高斯波束偏移
211212153_各向同性与各向异性逆时偏移成像对比分析
价值工程0引言随着勘探目标的日趋复杂化，如复杂盐丘及盐下、逆掩推覆构造、复杂高陡断裂及碳酸盐岩缝洞储集体等，逆时偏技术受到越来越多关注并在上述复杂构造成像方面取得了相对于单程波和射线偏移更好的效果。
但逆时偏移本身计算代价要远高于单程波和射线偏移，进一步考虑各向异性会大大增加方程复杂度，计算和存储量也随之大幅增加。
得益于计算机硬件和针对性优化算法的快速发展，上述制约各向异性逆时偏移广泛应用的难题不断得到有效解决。
由于在地质演化过程中地层挤压形变等作用，倾斜横向各向同性（Tilted TransverseIsotropy ，TTI ）甚至倾斜正交各向异性（Tilted Orthorhombic Isotropy ，TORT ）普遍发育。
以TTI 介质模型为例，当前TTI 介质逆时偏移多以声学近似展开，基于声学近似的TTI 介质耦合拟声波逆时偏移在实际应用过程中存在倾角剧烈变化时数值求解不稳定、菱形伪横波干扰以及求解效率低等难题。
为此，Fletcher 等[1]分析了波场在倾角剧烈变化的复杂构造区传播时出现的不稳定性问题，并给出了最大化降低横波反射能量的控制参数，将有限横波速度代入方程实现了拟声波波场的稳定传播。
Zhang 等[2]同样对波场传播的稳定性问题进行了原理上的深入分析，指出其稳定性与空间差分算子具有密切联系，通过对微分算子应用自共轭性质以保证其传播的稳定性，得到了TTI 介质拟声波方程。
Duveneck 等[3]直接从胡克定律出发，通过对应力和应变旋转亦推导出了TTI 介质拟声波方程，其波场物理意义明确且保守振幅能量。
由于基于声学假设得到的耦合声学近似方程难以避免特有的稳定性、菱形横波干扰等问题，近年来各向异性纯P 波方程也得到越来越多关注。
Crawley 等[4]将Etgen 等[5]提出的伪解析法VTI 介质纯P 波正演模拟算法拓展到TTI 介质，实现了TTI 介质的纯P 波逆时偏移。
各向异性弹性波叠后逆时深度偏移
ｐ ｐ ｒ ｓ ｕ ｉｄ ｉ ｓ ｅ ｆｒ ｖ ｒ ｅ ｔ ｅ ｃ ｎ ｉ ｕ ｔｏ ｌ ｏ ｉ ｍ ｎ ｂ ｏ ｂ ｎ ｏ ｎ ａ ｙ ｃ ｎ ｉｉ ｎ ｆ ａ ｉｏ ｒ ｐ ｃ ｅ ａ — ａ ｅ ｔ ｄ ｅ ｓ ｕ ｓｏ ｅ ｅ ｓ － ｉ ｏ ｔ ａ ｉ ｎ ａ ｇ ｒｔ ｍ ｎ ｈ ａ ｄ ａ ｓ ｒ ｉ ｇ ｂ ｕ ｄ ｒ ｏ ｄ ｔ ｓ ｏ ｎ ｓ ｔ ｏ ｉ ｌ ｓ ｏ
维普资讯
第 ４ １卷 第 ３ 期 ２０ ０ ２年 ９月
石
油
物
探
Ｖｏ ． １。 ． Ｉ ４ Ｎｏ ３
Ｓ ｐ． ０ ２ ｅ ｔ２ ０
ＧＥ０ＰＨ ＹＳＩ ＣＡＩ ＲＯＳ Ｐ ＰＥＣＴＩ ＮＧ ＿ Ｒ Ｆ ） ＰＥＴＲ０ Ｉ （ ＥＵＭ
ｇ ａ ｅ ｅ ｕ ｔ ｆｔ ｄ ｌ ｈ ｗ ｈ ｔ ｔ ｅ ｓ ｓ ｅ ｃ ｎ ｐｒ ｄｕ ｅ ｃ ｅ ｒ ａ ｄ ａ ｃ ｒ ｔ ｍａ ｓ ｏ ｏ ｐ ｅ ｅ ｌ ｇ ｃ ｌ ｒ ｔ ｄ ｒ ｓ ｌｓ ｏ ｗｏ ｍｏ ｅ ｓ ｓ ｏ ｔ ａ ｈ ｙ ｔ ｍ ａ ｏ ｃ ｌ ａ ｎ ｃ ｕ ａ ｅ ｉ ｇｅ ｆ ｃ ｍ ｌ ｘ ｇ ｏ ｏ ｉａ ｓ ｒ ｃ ｕ ｅ ，ｗ ｈ ｒ ａ ｈ ｓ ｔｏ ｉａ ｌ ｉ ｒ ｔ ｄ ｒ ｓ ｌｓｆ ｒ ｔ ｅ ｓ ｍｅ ｍｏ ｅ ｓ ｓ ｇ ｅ ｔｔ ａ ｒ ｄ ｔｏ ａ ｓ ｔ ｏ ｉ ｍｉ ｒ — ｔｕ ｔ ｒ ｓ ｅ ｅ ｓｔ ｅ ｉ ｏ ｒ ｐ ｃ ｌ ｍ ｇ ａ ｅ ｅ ｕ ｔ ｏ ｈ ａ ｄ ｌ ｕ ｇ ｓ ｈ ｔｔ ａ ｉｉ ｎ ｌ ｏ ｒ ｐ ｃ ｙ ｉ ｇ ａ ｔ ｎ ｃ ｎ ｃ ｕ ｅ ｃ ｎ ｉ ｅ ａ ｌ ｒ ｏ ｓ ｉ ａ ａ ｓ ｏ ｓｄ ｒ ｂ ｅ ｅ ｒ ｒ ． ｏ
TTI各向异性叠前深度偏移技术在库车复杂山地的应用
TTI各向异性叠前深度偏移技术在库车复杂山地的应用吴超;许安明;尚江伟;陈维力;朱婧;章国威【摘要】库车坳陷复杂山地地表地下双重复杂、地表高差大、山体发育、沟壑纵横,古近系膏盐岩挤压变形严重,厚度变化大,盐上浅层高陡,盐下目的层逆冲叠瓦断块发育,导致地震资料信噪比低,成像效果较差.通过多年复杂山地地震资料处理研究,形成了起伏地表TTI各向异性叠前深度偏移技术.叠前深度偏移速度建模中,采用了小平滑基准面,通过微测井约束层析反演计算静校正厚度、时间和速度,建立较高精度的浅表层速度模型;综合应用地表露头、重磁电、地质和钻测井资料约束建立了合理的中深层速度模型;采用井控TTI各向异性参数提取及网格层析成像技术提高了速度模型和成像的精度.通过TTI各向异性叠前深度偏移处理,库车坳陷复杂山地地震资料信噪比和成像质量明显提高,为区带研究及圈闭落实奠定了基础.【期刊名称】《新疆石油地质》【年(卷),期】2019(040)001【总页数】5页(P103-107)【关键词】库车坳陷;复杂山地;TTI各向异性;叠前深度偏移;层析反演;速度模型【作者】吴超;许安明;尚江伟;陈维力;朱婧;章国威【作者单位】中国石油塔里木油田分公司勘探开发研究院,新疆库尔勒841000;中国石油塔里木油田分公司勘探开发研究院,新疆库尔勒841000;中国石油塔里木油田分公司勘探开发研究院,新疆库尔勒841000;中国石油塔里木油田分公司勘探开发研究院,新疆库尔勒841000;中国石油塔里木油田分公司勘探开发研究院,新疆库尔勒841000;中国石油塔里木油田分公司勘探开发研究院,新疆库尔勒841000【正文语种】中文【中图分类】P631.4库车坳陷位于塔里木盆地北缘的山前复杂构造区，是在古生代前陆盆地之上叠加的新生代前陆盆地。
发育三叠系和侏罗系两套煤系烃源岩，生烃强度大、资源量大；同时发育以白垩系巴什基奇克组为储集层，古近系库姆格列木群膏盐岩为盖层的优质储盖组合。
地震深度成像
地震深度成像目前，有许多勘探目标无法使用常规地震成像方法进行识别，而利用叠前深度成像技术，作业公司可以对包括最复杂构造在内的地质特征进行清晰成像。
这种准确的结果可以降低风险并帮助确定储量。
在上个世纪，地震解释人员普遍采纳在时间域处理和显示的地震图像。
在目前许多热点勘探地区，尤其是由于断裂或盐体侵入导致构造复杂和地震速度突变的地区，时间域处理方法可能产生容易使人误解的结果，只有深度成像技术可以确定地下构造特征的形态和真实的位置。
在有些情况下，深度域和时间域图像之间的差异，可能会形成或否定一个远景目标，导致构造扩大或缩小圈闭范围、目标偏离数百英尺或数百米以及储量增多或减少等，这种差异可能换来一口代价昂贵的干井而不是一个发现。
本文描述了深度成像技术演变为一种复杂地层特征成像的地震数据处理特殊技术的过程。
文中的研究实例介绍了在美国陆上、墨西哥湾和北海地区作业的石油公司是如何利用深度成像技术来改善其钻探成功率的。
地震勘探技术发展大事记在已经过去的20世纪里，有一些显著的、里程碑式的事件反映出地震勘探技术的进步。
虽然许多新技术从引入到实际采纳的成熟期要花费10年时间，但每项技术最终又都创造出新的勘探机会。
从上世纪20年代开始，人们引入单次覆盖模拟记录方法检测倾斜地层（下一页）。
[1] 在30年代，这项新技术是在盐丘周围获得发现的关键技术，成为标准应用方法。
到50年代，出现了共深度点（CDP）叠加技术实现的多次覆盖地震数据，使信噪比得到明显改善。
在60年代，人们引入数字数据采集和处理技术，取代了早期的模拟和光点法，使地震数据质量发生了重大改进，并在世界许多地区获得新的发现。
在整个70年代，数字数据和二维勘探成为常规技术，这些技术一起打开了北海和其它富有挑战性地区的勘探局面。
虽然时间域处理技术是标准方法，但人们引入和测试了二维叠后深度偏移方法。
首批小规模三维勘探数据在一些已开发油田采集，以便改善油藏划分的能力。
叠前深度偏移技术
叠前深度偏移技术一、技术原理及主要技术内容叠前深度偏移技术已由克希霍夫积分法发展到波动方程法，同时还发展了其它的偏移方法，如：高斯束（Beam）偏移、相移屏偏移技术、转换波叠前深度偏移、各向异性叠前深度偏移等，现把上述各种方法分述如下：（1）克希霍夫积分法叠前深度偏移：该偏移方法一般由两部分组成：一部分是旅行时计算，另外一部分是克希霍夫积分处理。
偏移的精度主要取决于旅行时的精度。
旅行时计算建立在费马原理的基础上，即地下两点间的一切可能路径中实际路径对应于最小旅行时间。
它遵循倒转射线追踪机制，大多数情况下使用对应于体波而不是首波的射线，这样减少了偏移成像的畸变，且输出轨迹是灵活的。
新方法主要改进了原方法中单波至、不保幅的缺点，现在是计算多波至旅行时，并且具有振幅与相位保持特性，最具代表性的方法是由以色列PARADIGM公司发展的共反射角克希霍夫积分法，其原理与方法是：由成像点到地面采用照明式射线追踪；在每个射线均计算旅行时、观测位置、相位旋转因子、慢度；在特定倾角每对射线均是潜在反射；求和某成像点同一层的所有反射形成共反射成像道集；所有到达时的振幅与相位都是保持的。
高斯射线束（Gaussian Beam）偏移方法有别于常规的克希霍夫积分法深度偏移方法，目前只有Chevron公司使用它，它分多组射线束进行研究，采用Gaussian法振幅衰减与相位抛物线近似等。
具体讲它是将震源和接受点波场局部分解成“束”，并利用精确的射线追踪将这些束返回地下。
一个地面位置能发出几个束，不同的束对应不同的初始传播方向，每个束独立于其他束传播，且受单个射线管引导。
射线管可以重叠，所以能量能在成像位置、震源位置及接受点位置间以多个路径传播，因此高斯射线束偏移可处理多路径。
该种方法部分解决了常规克希霍夫积分法精度不高的问题。
（2）波动方程法叠前深度偏移：该种方法研究多波至，易振幅与相位保持，精度高，但费机时，主要方法有有限差分法（FD）与相移校正法（PSPC），它们均基于单程波动方程、平方根算子向下延拓，并使用多个参考速度。
叠前深度偏移培训总结汇报
观覆 测盖 系次 统数
偏 移 距 分 布
垂 直 函 数
浮 动 面
注意：CMP加载中有6个必须需道头，浮动面 处理中心叠前深度偏移培训班-第3期
CMP道集的加载与检查
道集类型
道头信息检查
二、
培 训 内 容
处理中心叠前深度偏移培训班-第3期
二、
培 训 内 容
CMP道集的加载与检查
地震属性图
处理中心叠前深度偏移培训班-第3期
注意：这是基础性工作，很重要!!!
处理中心叠前深度偏移培训班-第3期
设置层位表
模型建立
层位解释
解释也可调用解释层位
解释原则：1、尽量解释强反射界面； 2、一般200-500米解释一层； 3、在最后一层有效层后加一层兜底的层，作为控制层， 以免软件出错。
处理中心叠前深度偏移培训班-第3期
PICKS
各向异性PSDM
处理中心叠前深度偏移培训班-第3期
二、
培 训 内 容
剩余延迟
剩余延迟收敛前相干谱（3D)
处理中心叠前深度偏移培训班-第3期
二、
培 训 内 容
实体模型的剩余延迟收敛前后相干谱 初始剩余延迟
迭代后剩余延迟
处理中心叠前深度偏移培训班-第3期
二、
培 训 内 容
初始深度偏移剖面（2D)
迭代后的深度偏移剖面(2D) 处理中心叠前深度偏移培训班-第3期
模型建立
T_SURFACE
层位模型的建立
SOLID_MODEL
通过3D CANVAS用拾取的PICKS创建时间层位修复时间层位直至创建实体模型
处理中心叠前深度偏移培训班-第3期
模型建立
2、建立初始沿层层速度(MAP)
轴对称各向异性介质波动方程深度偏移的对称非平稳相移方法
第５ Ｏ卷 第 １ 期
２ ０ ０ ７年 １月
地
球
物
理学报 源自Ｖｏ １ ． ５ ０．Ｎｏ． １
ＣＨ Ｉ ＮＥ Ｓ Ｅ Ｊ ＯＵ ＲＮ Ａ Ｌ Ｏ Ｆ Ｇ Ｅ ＯＰ Ｈ ＹＳ Ｉ Ｃ Ｓ
Ｊ ａ ｎ ． ，２ ０ ０ ７
刘礼农 ， 陈 小宏 ， 张剑锋等 ． 轴 对 称 各 向异 性 介 质 波 动 方 程 深度 偏 移 的对 称 非 平 稳 相 移 方 法 ． 地球物理学报 ， ２ ０ ０ ７ ， ｓ ｏ （ １ ） ： ２ ７ ７～
性介质中的相移算子 ， 进 而 将 以 相 移算 子 为基 础 的 对 称 非 平 稳 相 移 方 法 推 广 到 各 向 异性 介 质 ， 发 展 了 一 个 三 维 各
向异 性介 质 的 深 度 偏 移 方 法 ．文 中 使 用 的 各 向 异 性 介 质 的速 度 模 型 与 现 行 的 各 向 异 性 构 造 的 速 度 估 计 方 法 一 致 ，
适用于三维 Ｔ Ｉ 介质 ．
关键 词 波 动 方程 ， 各 向 异性 介 质 ， 相移算子 ， 对 称 非 平 稳 相 移延 拓 算 子
文章编 号 ０ ０ ０ １ — ５ ７ ３ ３ （ ２ ０ 町） Ｏ ｌ — Ｏ ２ ７ ７— ０ ８ 中图分 类号 Ｐ ６ ３ １ 收 稿 日期 ２ Ｏ Ｏ ５ — ０ ８ — ０ ４ ， ２ ０ Ｏ ６ — ０ ８ — ２ １ 收 修 定稿
２ ８ ４
ｌ Ｊ ｉ ｕ Ｌ Ｎ， Ｃ ｈ ｅ ｎ Ｘ Ｈ， Ｚ ｈ ａ “ ｇ Ｊ Ｆ ， ｅ ｔ ａ １ ．Ｗａ ｖ ｅ ｅ ｑ ｕ ａ ｔ ｉ 。 ｎ ｄ ｅ ｐ ｔ ｈ ｍ ｉ ｇ ｒ ａ ｔ ｉ 。 ｎ ｗ ｉ ｔ ｈ ｙ ｍｍｅ ｔ ｒ ｉ ｃ ｎ ｏ ｎ － ｓ ｔ ａ ｔ ｉ 。 ｎ ａ ｒ ｙ ｐ ｈ ａ ｓ ｅ ｓ ｈ ｉ ｉ ｆ ｍｅ ｔ ｈ 。 ｄ ｉ ｎ ａ ｘ ｉ — ｓ ｙ ｍｍｅ ｃ ａ ｌ
晶体各向异性及工艺参数对蓝宝石切片加工面形偏差的影响
晶体各向异性及工艺参数对蓝宝石切片加工面形偏差的影响*王 超1， 葛培琪1,2， 贺基凯3， 王新辉3（1. 山东大学 机械工程学院， 济南 250061）（2. 山东大学，高效洁净机械制造教育部重点实验室， 济南 250061）（3. 青岛高测科技股份有限公司， 青岛 266114）摘要 在多线锯切片加工中，由于蓝宝石晶体材料的各向异性，不同加工位置处的力学性能不同，导致金刚石线锯在垂直进给方向上产生偏移，从而造成切片加工晶片面形偏差。
为深入研究晶体各向异性对切片加工晶片面形偏差的影响机制，通过分析蓝宝石晶体材料的特性，并计算常用晶面弹性模量分布情况，结合线锯受力情况模拟，计算晶片的面形偏差，分析工艺参数对面形偏差的影响。
结果表明：C 面、A 面和M 面晶片面形偏差不受蓝宝石各向异性的影响；R 面晶片切片加工时，可选择切片进给角度为90°或270°以获得较小晶片面形偏差；减小比进给速度或采用变速进给方法可降低晶片面形偏差。
关键词 蓝宝石；金刚石线锯；晶体各向异性；面形偏差；工艺参数中图分类号 TG74 文献标志码 A 文章编号 1006-852X(2023)05-0612-09DOI 码 10.13394/ki.jgszz.2022.0207收稿日期 2022-11-24 修回日期 2023-02-08蓝宝石因其优异的物理性能、化学性能和光学性能，已成为现代工业中极为重要的材料，在衬底材料等领域中广泛应用[1-2]。
蓝宝石一般加工流程为：切片、研磨和抛光，切片后晶片表面质量将直接决定后续加工效率和成本[3-4]。
蓝宝石硬度高且脆性大，加工难度较高，晶体点阵为六方密排堆积的方式排列，沿不同晶向材料性能不同，为典型的各向异性材料，因此晶体取向会影响蓝宝石材料性能[5]。
蓝宝石制备方法不同，其弹性模量和断裂韧性不同，不同晶面蓝宝石临界载荷不同，其硬度也不相同，蓝宝石晶体取向同样会影响其纳米塑性等[6-9]。
叠前深度偏移
7、沿层剩余延迟迭代
东方公司研究院处理中心
叠前时间偏移
8、处理效果 叠前深度偏移
东方公司研究院处理中心
培训内容
叠前深度偏移基本原理 三维各向异性叠前深度偏移流程 二维叠前深度偏移基本流程 几点体会
东方公司研究院处理中心
几点体会与认识
1、GeoDepth软件的速度模型建立以及层速度迭代 优化方法具有行业领先优势； 2、GeoDepth软件产品多、应用模块多、模块中选 件多，需要在今后的项目中逐步消化应用； 3、此次培训收获很多，一是更系统地了解深度偏 移的工作流程与操作，二是掌握了Geodepth的更多 功能，为下一步做好深度偏移、用好深度偏移提供 了保障。
东方公司研究院处理中心
3、建立初始沿层层速度
基本步骤： • 层位网格化,建模 (MAP) • 生成时间域层速度 (CVI) • 沿层抽取沿层层速度并建模 （MAP） 1.编辑层速度平面图并建模 （MAP）很重要
东方公司研究院处理中心
4、更新沿层层速度
初始叠前深度偏移
CRP道集是否拉平
Y
N
计算沿层剩余延迟
理论模型
叠加剖面 时间偏移
叠加剖面 时间偏移
深度偏移
东方公司研究院深处度理偏移中心
各向异性叠前深度偏移工作思路
1
2
工区建立数据加载
道集的加载与检查 速度的加载与检查
初始速度模型的建立
层位模型的建立 创建初始深度域层速度体
4
各向异性深度偏移
井约束求取各向异性参数 各向异性叠前深度偏移
东方公司研究院处理中心
二维叠前深度偏移基本流程
工区建立
1.数据加载
TTI各向异性逆时偏移技术及应用
TTI各向异性逆时偏移技术及应用王咸彬【摘要】地下介质广泛存在各向异性,传统各向同性地震偏移成像技术往往会导致成像精度不高甚至深度偏差问题,宽方位采集技术和高精度逆时偏移(RTM)成像技术的应用更是突显了各向异性的影响.从弱各向异性弹性波波动方程出发,首先采用拟声波近似得到VTI各向异性伪声波控制方程,然后引入交叉导数项进行坐标旋转得到TTI各向异性伪声波控制方程,再由高阶有限差分方法得到TTI-RTM偏移算子,最后采用波场校正消除横波分量影响,提高各向异性偏移算子的精度.模型试算和实际资料处理结果表明,该技术在处理各向异性介质地震资料时具有更高的精度,是高精度地震成像理想的技术手段.%Anisotropy is widespread in the subsurface formation medium.The conventional isotropic seismic migration imaging technique often leads to low imaging accuracy and even depth error.The application of wide azimuth seismic acquisition technique and high precision RTM imaging technology have highlighted the influence of anisotropy.In this paper we begin with the weakly anisotropic elastic wave equation,First,the VTI anisotropic pseudo acoustic wave control equation is obtained by quasi acoustic approximation.Then,the cross derivative term is introduced to coordinate rotation to get the TTI anisotropic pseudo acoustic wave control equation.Next,the TTI-RTM migration operator is obtained by high order finite difference.Finally,the influence of the S-wave component is eliminated by wave field correction to improve the accuracy of the anisotropic migration operator.The model test and field data processing results show that,the technique is more beneficial to theanisotropic medium with higher precision and it's an ideal technique for high precision seismic imaging.【期刊名称】《石油物探》【年(卷),期】2017(056)004【总页数】9页(P534-542)【关键词】各向异性;拟声波近似;控制方程;TTI逆时偏移【作者】王咸彬【作者单位】中国石油化工股份有限公司石油物探技术研究院,江苏南京211103【正文语种】中文【中图分类】P631地下介质具有广泛的各向异性特性,如我国东部陆相砂泥岩薄互层具有长波长各向异性特征,西部海相岩溶—裂缝型碳酸盐岩储层具有裂隙诱导各向异性特征。
Tesseral 培训讲义1-概述
X,m
в
A=A(x) X-comp. Arefl
г
当水平层状介质之上有近地表的复杂介质时，把激发点、接收点都放在地表时，AVO振幅的规律变化已经无法分析 和辨认了，模型正演的结果表明，对带有近地表复杂介质的陆上地震资料进行AVO分析是难以取得理想结果的。
TESSERAL 的特色 各向异性介质模拟,也说明了各向异性介质对偏移精度的影响
通讯地址：北京海淀区北三环中路39号418室，邮编：100088，电话：（010）62374349，
传真：（010）62374349，E-mail:ws9663@.
Tesseral 2-D产品开发的简短里程碑： • 从1997-1999集中力量开发Tesseral 2D的第一版本， 进行市场调研和宣传。
左图为VSP-CDP 转换后的时间剖面，右图为偏移剖面
TESSERAL 的特色各向异性介质模拟,下图为在各向同性和各向异性介质中模拟 的VSP 波场对比
a
左图展示了TTI 各向异性介质 中有关透射转换横波的例子： (a)模型中下面一层为各向异性 介质，它的倾斜对称轴为45º 。 模拟的垂直VSP 炮集(波场 Z分量)示于图(a) ，其中的绿色箭 头指示在各向异性介质中传播 的转换波。 (b)在模型下层为各向同性介质 中模拟，得到的垂直VSP 炮集( 波场 Z-分量)，可见在各向同性 介质中，这样的转换波是不存 在的。
Exploding surface和Exploding reflectors层爆炸震源方允许用户很快地得到近似算 法的地震剖面。
浏览器：
可以采用适当的方式,浏览由计算引擎计算的合成地震记录和波场传播的快照。并输出计 算结果，以便在Tesseral 软件包中，或以外的其他系统，进行后续处理。它也可以非常 方便地对比不同的波场分量,改变信号的振幅值等。 解释人员可以：
各向异性弹性波高阶有限差分法叠前逆时深度偏移
１
：＝ ｃ ／ Ｐ 薹 ＋ （ ２ Ｉ 一 ） ｃ 州 ） ；ｋ） ＋ ２ ￣－ ¨一 ，１ ／ ｊ） －／ （２ ｚ ＋ ｃ＇ｋ１；ｌ１ 薹ｌ ；ｊ ２ｋ ；２ ｉ，＋ ｃ＇－－ｉ－ Ｃ，－ －ｊ Ⅳ ｉ＋ ；ｊ） ￣ｊ ｋ２ ｉ － ＋－ ＋ ＋ ＋ ） ’＋ｋ ＋ ＋ ｚ１ －１ ２ ， ｌ ，
＿ Ｌ ４ 鲁 ４ －４ － 乙 ｃ ４
ａ ８ Ｉ ｒ
ｌ Ｄ 为介质密度。Ｔ ｏ ｓｎ ｈｍ ｏ 参数符号定义为
ｃ：
Ｉ 百 Ｄ 百
：
十百 采用 高阶交错 网格有 限差分法［２ １４ ６１ （） －将式 １ 离 散为任意２ Ｎ为任意正整数） Ｎ（ 阶精度 ， 有
』 鲁＋ Ｄ 十 警
相对较少 。为此 ， 本文以各 向异性全弹性波动方程
为基础 ， 在前人的研究基础之上 ， 构造了适用于各 向 异性介质中弹性波多分量双程波联合偏移的叠前逆 时偏移算子 ， 以各 向异性介质 中最具代表性的横 向 各 向同性介质为例 ， 用数值实例来验证本文介绍 的 逆时偏移算子的可实用性和有效性 。
随 着弹性 波 动理 论研 究 和野外 地 震勘探 的不 断 深入 ， 量科 学研 究结 果 表 明 ， 大 地球 介 质具有 普 遍 的
前 有 限 元 法叠 前 历 史 深 度 偏移 ， 秉 铭 等ｎ研 究 过 张 。 １ 各 向 同性 介 质 有 限差 分 法 弹性 波 叠 后 逆 时偏 移 问
像 更加清晰准确 ， 偏移成像 精度较高 ， 并且 因此开展各 向异性叠前逆 时成像 可为 当前 高精度岩性地 震勘探提供 方
法指导。
关键词 ： 阶交错 网格 有 限差分 法 ； 向异性弹性 波动 方程 ； 高 各 叠前逆 时深度偏 移 ； 最大 绝对 振幅能 量； 多分量逆 时
201209各向异性克希霍夫叠前深度偏移(石油物探)
文章编号 %...%''%#$.%$$.4.'@+"%.收稿日期 $.%$.#.*'改回日期 $.%$.@%.(作者简介 周巍#%*+#)$"女"高级工程师"现主要从事各向异性克希霍夫叠前深度成像工作("各向异性克希霍夫叠前深度偏移周"巍 王鹏燕 杨勤勇 方伍宝 潘宏勋 刘旭跃 郭书娟中国石油化工股份有限公司石油物探技术研究院 江苏南京$%%%.#摘要 随着地震勘探精度的提高"各向异性对于常规克希霍夫叠前深度偏移效果的影响已不容忽视(研究了利用1波非双曲旅行时求取各向异性参数的方法"当道集的排列长度远大于勘探目的层深度时"1波走时资料包含的信息可以用来确定C (V 介质的各向异性参数(利用C (V 介质的!\1道集"研究了二步法提取,参数的方法"即先利用短排列地震数据求取均方根动校正速度"再利用长排列数据和水平速度扫描法得到水平速度"通过转换得到各向异性,参数'采用基于声学近似的C (V 介质程函方程"通过各向异性射线追踪计算旅行时(模型及实际资料处理结果表明"利用上述方法提取各向异性参数"进行C (V 介质各向异性克希霍夫叠前深度偏移"使成像精度和质量得到了明显提高(关键词 各向异性参数'叠前深度偏移'速度分析'射线追踪'旅行时Z %#%%.3#*+*!a 36H H >3%...;%''%3$.%$3.43..*中图分类号 1+#%3'文献标识码 -""克希霍夫叠前深度偏移成像技术具有简便*高效*适用性强等优点"在复杂地区地震勘探中发挥了重要作用(但是近年来随着勘探技术的不断发展"对成像精度的要求也越来越高"基于各向同性理论的克希霍夫叠前深度偏移已不能满足实际生产的需要"有必要开展各向异性克希霍夫叠前深度偏移研究(目前"地震资料处理中所指的各向异性介质主要是横向各向同性介质"如上覆地层中广泛存在的页岩和薄互层"会产生具有垂直对称轴的极性各向异性#C (V "横向各向同性$(影响横向各向同性介质的克希霍夫叠前深度偏移的一项关键技术是各向异性射线追踪"因为克希霍夫叠前深度偏移的主要计算量是射线追踪计算旅行时(各向异性克希霍夫叠前深度偏移与各向同性克希霍夫叠前深度偏移的不同之处就体现在旅行时的计算上"幸运的是这方面的研究工作近年来取得了很大进展+%;4,(影响各向异性介质精确成像的另一项关键技术是各向异性参数的可靠估计(实际生产中深度偏移预测的目标位置不准确并不表明该方法本质上存在错误"而是速度估算方法存在不足(因此"在克希霍夫各向异性叠前深度偏移中"建立正确的地下介质速度场和各向异性参数场#,参数$是非常重要的(对于横向各向同性介质"%*&+年(0D F H 8>提出了具有明确物理意义的各向异性参数-*!和."并分别给出了基于这些参数的相速度和时差速度在任意强度以及弱各向异性条件下的表达式++,"为各向异性研究奠定了基础(之后有不少学者对横向各向同性介质的速度特性及各向异性参数估计方法作了大量研究"但主要针对纵波非双曲线时差速度分析和旅行时反演+@;%%,(我们基于-7O 0=76K =0和B 98J 0O =等提出的非双曲线时差各向异性参数估算方法+@;*,"对,参数提取和深度域模型建立进行了讨论'提出应用!\1道集数据按偏移距大小分段处理的方法提取动校正速度J >F D 和,8K K 等效参数"并通过P 6T 公式获得层间,#&$值和层速度模型'通过深度偏移不断修正速度模型"实现了各向异性克希霍夫叠前深度偏移处理(%"各向异性克希霍夫叠前深度偏移方法""各向异性克希霍夫叠前深度偏移的旅行时计算采用各向异性介质的射线追踪方法(我们采用-7O 0=76K =0从N1波标量波动方程出发推导出的基于声学近似的C (V 介质程函方程+%,#%M $,$J $>F D%).%#$0$M%).%#$2+,$M J $1.%).%#$C $N $,J $>F D J $1.%).%#$0$M%).%#$2+,$%).%#$C $K %#%$通过特征值法"推导出射线追踪方程组+%,+@'第4%卷第4期$.%$年*月石"油"物"探B )X 15S A V !-Y1Q X A 1)!(V "BZ X Q1)(Q X Y )[\CD 734%""D 34A 8E3"$.%$I 0I TK J $>F D Q 0+%M $,#%N Q $C J $1.$,I C I TK #%N $J $>F D $,Q $R Q C J $1.I Q 0I T KN J >F D #%M $,$Q $R %J %0N J $>F D Q $R ,0M J >F D Q $R Q $C J $1.$,%J %0M J >F D ,#$0N #%N $J $>F D ,Q $R $Q $C J 1.%J 1.%0I Q C I T KN J >F D #%M $,$Q $R %J %C N J $>FD Q $R ,C M J >F D Q $R Q $C J $1.$,%J %C M J >F D ,#$C N #%N $J $>FD ,Q $R $Q $C J 1.%J 1.%0I G I TK J $>F D #%M $,$Q $0M #%N 'J $>F D ,Q $0$Q $C J $1.I Q 2I T KN J >F D #%M $,$Q $R %J >F D %2N J $>F D Q $R ,2M J >F D Q $R Q $C J $1.$,%J >F D %2M J >F D ,#$2N #%N $J $>F D ,Q $R $Q $C J 1.%J 1.%2#$$式中%J >F D 为"\X 速度'J 1.为垂直速度',为各向异性参数",i -e !%d $!'Q R "D C "D 0分别为R "C "0方向的慢度(从公式#%$中可以看出"用于各向异性克希霍夫叠前深度偏移的参数有#个%时差速度*垂直速度和,参数(除了垂直速度由井资料获得外"其它两个参数都可以直接在地震数据上拾取得到($"C(V 介质各向异性速度分析对于C (V 介质"1波反射旅行时的非双曲性非常明显(-7O 0=76K =0和(H M =>O 6>给出了大偏移距单一界面1波非双曲线反射旅行时公式+&,%G $#0$K G $.M 0$J $>F D N $,0'J $>F D +G $.J $>F D M #%M $,$0$,##$##$式表明"在各向异性介质中"反射1波的旅行时不再符合双曲线轨迹(图%是根据##$式计算的单一界面反射波时距曲线"界面深度是%...F (从图%看出%#当偏移距大于界面深度时#0-%...F $"各向异性对反射波旅行时的影响非常明显"偏移距越大影响也越大"所以大偏移距资料适合各向异性参数提取'$,只对大偏移距时距曲线有影响",越大表明各向异性越强"对旅行时曲线影响越大'%当偏移距小于界面深度时",影响不明显(当,i .时"公式##$变成各向同性介质的双曲线形式#图%中红色曲线$(图%"单一界面各向异性时距曲线#界面深度为%...F ""\X 速度为$...F !H $""根据,参数的定义"1波水平速度可以表示为J 0D 9K J 1.J >F D#'$把#'$式代入##$式得G $#0$K G $.M 0$J $>F D N #J $0D 9N J $>F D $0'J $>F D #G $.J '>F D M J $0D 90$$#4$为了更好地适合过大偏移距的情况#%L 4,0!C ,$L 4$"B 98J 0O =和(H M =>O 6>在公式#4$中J $0D 90$项前面引入常系数)i %L $"改善了对过大大偏移距的适用性+*,"即G $#0$K G $.M 0$J $>F D N #J $0D 9N J $>F D $0'J $>F D #G $.J '>F D M )J $0D 90$$"#+$公式#+$的应用效果通常好于公式##$"因为水平速度J 0D 9和时差速度J >F D 是同一个量纲"水平速度谱比,谱的能量更聚焦"有利于拾取速度函数(求取水平速度有两种方法%#一步法"即采用J >F D 和J 0D 9联合扫描的方法得到J >F D 和J 0D 9'$两步法"先采用近偏移距数据确定J >F D "然后用远偏移@@'第4期周"巍等3各向异性克希霍夫叠前深度偏移距数据获得J 0D 9(对于多层介质"1波反射旅行时公式为+@,"G $#0"-$K G $.#-$M 0$J $>F D #-$N J $0D 9#-$N J $>F D #-+,$0'J $>F D -G $.-J '>F D -M )J $0D 9-0$#@$其中"J 0D 9#-$是层状介质的有效水平速度(根据C (V 介质动校正公式+@,"有,#-$K %&%J '>F D #-$G .#-3$+2&-&K%#J #&$>F D$'#%M &,#&$$G #&$,.N 4%#&$定义,#-$K %G .-&-&K%#J #&$>F D $'#%M &,#&$$G #&$.K %G .-&-&K %#J #&$>F D $$+'#J #&$0D 9$$N "##J #&$>F D $$,G #&$.#*$应用P 6T 公式"得到层间J #&$0D 9和,#&$"即J #&$0D 9K J #&$>F D#%.$,#&$K %&#J #&$>F D $',#&$G .#&$N ,#&N %$G .#&N %$G .#&$N G .#&N %+$N ""#J #&$>F D $,'#%%$""对每层的速度和,参数进行时深转换"得到深度域的层速度和层间,参数"再由声波测井或C A 1测井得到垂直速度"进行各向异性叠前深度偏移处理(速度分析方法与各向同性叠前深度偏移类似"也需要进行多次迭代(图$是各向异性克希霍夫叠前深度偏移和速度分析的流程(图$"各向异性克希霍夫叠前深度偏移和速度分析流程#"模型试算Q H N "`(#介质理论模型各向异性参数提取设计一C (V 介质各向异性模型"根据公式#+$和公式##$"采用二步法计算了该模型的速度谱和,谱#图#$"模型参数如表%所示(图#"C(V 模型#参数见表%$速度谱估计=小偏移距正常时差速度谱',大偏移距水平速度谱'J 各向异性,参数谱&@'石"油"物"探第4%卷表%"C(V模型各向异性参数J>F D F H e% ,' F%....&..%4...3.#.%&..%*...3%$4$@..$'...3%'44@+. ""从图#可见 水平速度谱的能量聚焦要好于,参数谱 更便于拾取速度 图'是各向同性和各向异性动校正结果 由于介质是各向异性的 进行各向同性动校正后 远偏移距处的同相轴没有拉平 图', 而采用水平速度谱中拾取的水平速度做非双曲时差动校正 则能够使远偏移距处的同相轴拉平 图'J 表$给出了各向异性参数,的提取结果与理""""论值 对比可见 相对误差都在%.c以内 精度较高 说明提取各向异性参数的方法是有效的 表$",参数理论值与计算值的对比分析,理论值,计算值误差第%层...第$层.3.#..3.#$+3@c第#层.3%$4.3%$4.第'层.3%'4.3%'.e#3'cQ H L"薄互层模型各向异性叠前深度偏移对于厚度远小于地震波长的薄互层介质 可以将其近似为具有一定方向的横向各向同性 (V 介质 因此设计了一个水平各向同性薄互层模型 模拟C(V介质 如图4所示图'"各向同性和各向异性动校正=!\1道集记录 ,各向同性动校正 J各向异性动校正*@'第4期周"巍等3各向异性克希霍夫叠前深度偏移""模型共分'层 第% # '层是各向同性介质 厚度均为4..F 速度依次为为$... '... +...F H 第$层是薄互层 小层厚度是$F 总厚度是4..F 低速层速度是$&..F H 高速层速度是#$..F H 由于波长远大于小层厚度 所以薄互层应该具有较强的各向异性 模拟采集了%**炮记录 炮间距是%..F 每炮+.%道道间距是%.F 采样间隔$F H 记录长度$H 子波主频是#.5/采用本文介绍的方法对图4所示薄互层模型!P 1%#..道集进行了速度分析图+ 图@ 薄互层被当作一个厚层处理 相当于针对#个厚层进行速度分析 4.. @.. %%4.F H 附近 图+是小偏移距正常时差速度分析结果 在速度谱中可以看到#个强能量团 蓝色 分别对应#个厚层 图+J """"是正常时差动校正结果 可以看到 小偏移距道集被拉平 但因为各向异性的影响 远偏移距道集未拉平 存在较大的剩余动校正量 图@是远偏移距道集的非双曲时差速度分析结果 在水平速度谱上有#个明显的能量团蓝色 可以方便地拾取水平速度 通过非双曲时差动校正 远偏移距道集被拉平 图@中未被拉平的同相轴是由于受到低速薄层的影响所致根据拾取得到的水平速度和公式 &至公式 %% 得到层间各向异性, & 参数 图&=是建立的深度域层间,参数场 可以看出 薄互层的各向异性参数,高达.3.#. 属于强各向异性 图&,是建立的深度域层速度剖面图*是克希霍夫各向同性叠前深度偏移和各向异性叠前深度偏移结果对比 可以看出 各向同""""图4"薄互层模型图+"小偏移距正常时差速度分析=速度谱 ,!P 1%#..道集 J 正常时差动校正结果.&'石"油"物"探第4%卷图@"大偏移距非双曲时差速度分析=速度谱 ,!P 1%#..道集 J非双曲时差动校正结果图&"薄互层模型的深度域层速度和各向异性参数剖面=深度域,参数剖面 ,深度域层速度剖面%&'第4期周"巍等3各向异性克希霍夫叠前深度偏移图*"薄互层模型各向同性= 与各向异性 , 克希霍夫叠前深度偏移结果性偏移成像由于受薄互层的强各向异性影响 在薄互层下方%4..F 附近的各向同性反射界面 第#个反射界面 成像很不清楚 不仅如此 薄互层底界面 %...F 附近 也没有各向异性偏移成像清楚 分辨率低 而通过各向异性克希霍夫叠前深度偏移成像后 薄互层底界面及下方的各向同性反射界面的成像都很清楚 分辨率明显提高 说明各向异性速度分析方法和各向异性克希霍夫叠前深度偏移方法是正确的'"实际资料应用效果分析图%.是新疆某地区实际单炮记录 每炮""""有#$.道 道间距是4.F 属于大偏移距资料 有利于做各向异性分析 我们对该资料进行了各向异性速度分析和参数提取 图%%是建立的深度域层速度场 图%$是深度域水平速度场 图%#是深度域层间,参数场 利用这些参数进行了各向同性克希霍夫叠前深度偏移 图%' 和各向异性克希霍夫叠前深度偏移 图%4 对此图%'和图%4可以看出 在浅层 特别是在深度小于$...F 的浅层 矩形框所示 各向异性叠前深度偏移明显好于各向同性叠前深度偏移 同相轴的横向连续性有明显改善 在图%'和图%4中箭头所示位置各向同性叠前深度偏移结果出现杂乱反射 而各向异性叠前深度偏移结果出现强能量且连续性较好的同相轴 各向异""""图%."新疆某地区实际单炮记录$&'石"油"物"探第4%卷图%%"深度域层速度剖面图%$"深度域水平速度剖面图%#"深度域层间/参数剖面#&'第4期周"巍等3各向异性克希霍夫叠前深度偏移性叠前深度偏移之所以对浅层成像效果改善明显 是因为浅层的入射角范围更大 即偏移距与深度的比值 0 C更大 地层的各向异性影响更大的缘故 图%+是各向同性叠前深度偏移和各向异性叠前深度偏移的深度对比图 可以看到两者存在明显的深度差图%'"各向同性克希霍夫叠前深度偏移剖面图%4"各向异性克希霍夫叠前深度偏移剖面图%+"各向同性和各向异性叠前深度偏移结果对比'&'石"油"物"探第4%卷4"结束语各向异性克希霍夫叠前深度偏移效果明显优于各向同性克希霍夫叠前深度偏移(本文针对横向各向同性介质"研究了基于非双曲时差旅行时的各向异性参数提取*深度域速度模型和各向异性参数模型建立的方法"给出了各向异性克希霍夫叠前深度偏移的处理流程(通过模型和实际资料处理试验"证明上述方法及其处理流程是正确有效的"可以在实际生产中推广应用(本研究仅限于C(V 介质"没有涉及倾斜地层"这方面有待开展进一步的研究工作(致谢 本文研究过程中得到了中国石油化工股份有限公司石油物探技术研究院蔡杰雄 王于静的很大帮助 在此一并表示感谢参"考"文"献+%,"-7O0=76K=0(3->=J D<H:6JL=M88N<=:6D>K D9=>6H D:9D E6JF8I6=+`,3B8D E0G H6J H"$..."+4#'$%%$#*;%$4.+$,"]0<("B9=G A5"W=>?P3R6>8F=:6J=>II G>=F6J 9=G:9=J6>?6>=>6H D:9D E6JF8I6=%(08D9G=>I=E E76J=;:6D>+`,3)T E=>I8I-,H:9=J:H D K@4:0->><=7V>:89>=:A)B\:?"$..4"*+;**+#,"Y=L:D>P!"V H==J`5"C8H:9<FQ W"8:=73A76E;H76I6>=L=G;H D F8E9=J:6J=76F E76J=:6D>H D K H86H F6J M87D J;6:G=>6H D:9D E G D>I8E:06F=?6>?+`,3(08Y8=I6>?)I?8"$..%"$.#%$%@.;@#+',"段鹏飞"程玖兵3(V介质射线追踪与局部角度域叠前深度偏移+`,3A1B!A)B深圳$.%%国际地球物理会议论文集"$.%%"*'.;*'4P<=>1Z"!08>?`U3Q=G:9=J6>?6>(VF8I6==>I7D J=7=>?78I D F=6>R69J00D K K E98H:=J OI8E:0F6?9=;:6D>+`,3A1B!A)BA08>/08>$.%%V>:89>=:6D>=7B8;D E0G H6J=7!D>K898>J8A G FE D H6<F"$.%%"*'.;*'4+4,"(H M=>O6>V"(0D F H8>Y3V>M89H6D>D K98K78J:6D>:9=M;87:6F8HK D9:9=>H M89H8=>6H D:9D E G+`,3B8D E0G H6J H%**4"+.#'$%%.*4;%%.@++,"(0D F H8>Y-3W8=O87=H:6J=>6H D:9D E G+`,3B8D;E0G H6J H"%*&+"4%#%.$%%*4';%*+++@,"-7O0=76K=0(3C87D J6:G=>=7G H6H<H6>?>D>0G E89,D76JF D M8D<:6>:9=>H M89H87G6H D:9D E6J F8I6=+`,3B8D;E0G H6J H"%**@"+$#+$%%&#*;%&4'+&,"-7O0=76K=0("(H M=>O6>V3C87D J6:G=>=7G H6H6>:9=>H;M89H87G6H D:9D E6J F8I6=+`,3B8D E0G H6J H"%**4"+.#4$%%44.;%4+++*,"B98J0O=C"(H M=>O6>V3Z8=H6,676:G D K>D>0G E89,D76JF D M8D<:6>M89H6D>6>:9=>H M89H87G6H D:9D E6J F8I6=+`,3B8D E0G H6J H"%**&"+###$%*4@;*+*+%.,"B98J0O=C"18J0-"(H M=>O6>V"8:=73C87D J6:G=>=7;G H6H K D9:67:8I:9=>H M89H87G6H D:9D E6JF8I6="=E0G H6J=7F D I8776>?8T=F E78+`,3B8D E0G H6J H"$..%"++##$%*.';*%.+%%,"Y8A:<>K KS"B98J0O=C"(H M=>O6>V3P8E:0I D F=6> M87D J6:G=>=7G H6H6>C(VF8I6=<H6>?H<9K=J81;L=M8I=:="6H6:K8=H6,786+`,3B8D E0G H6J H"$..%"++##$%&*@;*.#编辑 戴春秋/////////////////////////////////////////////上接第'+&页"""!=6bY3->8K K8J:6M8F8:0D I:DH<E E98H H=J D<H:6J L=M8=>I06?08>89?G>D6H8K98N<8>J G;I6M6H6D>=77G=>I=I=E:6M87G+`,3X67B8D E0G H6J=719D H E8J:6>?"%***"#'#'$%#@#;#&.+*,"王君"周兴元"恽春华3地表一致性约束下异常振幅衰减技术的实现+`,3石油物探"$.%%"4.#4$%'*#;'*&W=>?`"]0D<b S"5<6!53Q8=76/=:6D>D K H<9K=J8J D>H6H:8>:J D>H:9=6>8I=>D F=7G=F E76:<I8=::8><=;:6D>:8J0>D7D?G+`,3B8D E0G H6J=719D H E8J:6>?K D918;:9D78<F"$.%%"4.#4$%'*#;'*&+%.,"刘伟"曹思远"王征"等3随机噪声的多尺度多方向域压制方法研究+`,3石油物探"$.%%"4##$%#.%;#.+Y6<W"!=D AS"W=>?]"8:=73\<7:6;H J=78_F<7:6;I698J:6D>H<E E98H H6>?F8:0D I K D99=>I D F>D6H8+`,3B8D E0G H6J=719D H E8J:6>?K D918:9D78<F"$.%%"4##$%#.%;#.++%%,"刘明乾"于相海"周夏丽"等3非纵测线地震资料相干噪声压制+`,3石油物探"$.%%"4.##$%#.@;#.*"#%&Y6<\^"S<b5"]0D<bY"8:=73!D0898>:>D6H8=::8><=:6D>6>,9D=I H6I8H86H F6JI=:=+`,3B8D E0G H6;J=719D H E8J:6>?K D918:9D78<F"$.%%"4.##$%#.@;#.*"#%&编辑 陈"杰4&'第4期周"巍等3各向异性克希霍夫叠前深度偏移。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Don Lawton University of Calgary Calgary, Alberta, Canada
Fold-fault Research Project
1580 m
Sideslip
480 m
Migrated section (isotropic)
Depth (m)
Migrated section (anisotropic)
Depth (m)
Lateral position error in target
hinterland foreland
V
...you know the closer your destination, the more you’re slipslidin’ away ...
(Paul Simon, 1977)
Vertical Transverse Isotropy (VTI)
V90 (fast)
ε = (V90 - V0)/V0 4[V δ = 4[ 45/V0 - 1] - [V90/V0 - 1]
V
90
(fa st
V45
V0 lo (s w )
)
Conclusions
• • • • •
Clastic rocks are anisotropic Anisotropy impacts depth imaging Depth and position errors will occur if isotropic migration velocities are used Vεδθ models are required for migration Values of anisotropy parameters are needed
V90 V0
Tilted transverse isotropy (TTI)
wavefront
ow Sl
Fa st
Group
Phase
Implementation – prestack Kirchhoff
source
midpoint
receiver
Aplanatic surface
• Raytrace through velocity model at all sources and receivers • Create gridded traveltime tables • For each trace, sum source and receiver tables • Distribute weighted trace amplitudes over aplanatic surfaces
V0 (slow)
Photo courtesy D. Spratt
Isotropic wavefront ε = 0.0 δ = 0.0
V V
Anisotropic wavefront (VTI) ε = 0.2 δ = 0.05
V90 V0
Anisotropic wavefront (VTI) ε = 0.2 δ = −0.2
Time migration“imaging velocity”
Anisotropic depth migration
[V0, ε, δ, θ]
Azimuthal refraction line layout
Dip
St rik e
45
o
Multioffset VSP
A VSP well B
) st fa (
V w) lo (s
DEPTH
Actual position
Imaged position (isotropic)
Lawton et al., 2001
Thrust fault model (layered fibreglass, homogeneous, anisotropic)
Anadarko Canada Corp. BG Canada Exploration Burlington Resources ConocoPhillips Devon Canada Encana Corporation Ecopetrol Unocal/Northrock Midland Valley Husky Energy Nexen Energy PetroCanada Resources Shell Canada Resources Suncor Energy Talisman Inc Occidental Oil and Gas Landmark Graphics
Vestrum et al., 1999
Anisotropic depth migration
Data courtesy Husky & Talisman.
Vestrum et al., 1999
Tarata Thrust Zone, Taranaki Basin
from Blackburn (1998)
Acknowledgements
Rob Vestrum, Jennifer Leslie, Graziella Kirtland Grech, Rachel Newrick, Helen Isaac,
Sponsors of the Fold-fault Research Project (FRP)
0.0 0.0 Distance (km) 5.0
Depth (km)
1.0
2.0
Migrated section (isotropic)
Migrated section (anisotropic)
Isotropic depth migration
Data courtesy Husky & Talisman.
Natural Sciences and Engineering Research Council
Base step
Isotropic migration operator
midpoint
Anisotropic migration operator
midpoint
FRP TTI step model
Shotpoint 0m midpoint geophone 5100 m
1500 m
Anisotropic
FRP phenolic TTI step model
FRP TTI step model
0m 5100 m
1500 m
Anisotropic Vo = 2925 m/s ε = 0.22 δ = 0.10
1580 m
480 m
Multichannel record, step model
Base of anisotropic layer Top step