三维（3D）反射地震学是成像地下火山的一种强大方法。通过使用来自塔里木盆地的工业三维地震数据，我们展示了如何从地震数据体中提取地下火山的岩床和通道结构。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
本文内容
摘要
摘要
管道系统的形态与结构可为玄武岩熔岩流的喷发速率和喷发方式提供关键信息。研究地下地质体最有效的方法是采用工业级三维 反射地震成像技术。然而，对地下火山体进行成像的策略与油气储层的成像方法存在显著差异。本研究中，我们处理了来自中国塔里木盆地北部的地震数据体，以说明如何通过不透明度渲染技术可视化岩床，以及如何利用时间切片技术对岩浆通道进行成像。在第一个案例中，我们利用标记岩床与围岩地层接触界面的地震层位隔离出目标区域，并应用不透明度渲染技术从地震数据体中提取岩床。所获得的精细岩床形态显示，岩浆流动方向从穹隆中心指向边缘。在第二个地震数据体中，我们采用时间切片技术对通道结构进行成像，这些通道对应于围岩中的明显不连续界面。一组在不同深度获取的时间切片表明，塔里木溢流玄武岩由中心式火山喷发形成，并由各自独立的管状通道供岩浆。
引言
大多数沉积盆地工业地震成像项目的目标是勘探油气储层。近年来，由于许多火山成因盆地含有大量火成岩且具备可观的油气储层，油气勘探已扩展至含有大量火成岩的盆地。然而，由于火山成因盆地中火成岩界面的存在，地震数据处理面临一系列由多种侵入体引起的挑战，例如能量透射减弱、固有衰减、干涉效应、折射和散射1。因此，油田公司正致力于减少此类"负面影响"对地震成像的影响2,3,4。
由于火成岩体与围岩之间存在显著的声阻抗差异，因此利用二维或三维地震反射成像技术可轻松识别沉积盆地中的火成岩体1,5,6。该方法能够清晰呈现火山管道系统中垂直与水平构造的精细图像7,8,9,10,11,12,13。然而，对地下火山体的成像策略与油气勘探中的方法存在显著差异8,14,15，这限制了工业地震数据在地下火山研究中的广泛应用，仅有少数成功案例报道10,15,16。本文报道了专为地下火山解释定制的地震数据处理详细流程。我们处理了两套地震数据体TZ47和YM2（图1），以展示如何在中国塔里木溢流玄武岩区可视化埋藏的火成岩体17。
访问受限。请登录或开始试用以查看此内容。
方案
注意：数据处理流程包括：合成地震记录计算、合成-实际地震道相关分析以及地质体提取。以下是每个步骤的详细说明。
1. 合成地震图的计算
- 计算井下测井曲线每个区间的声阻抗。
注意：声阻抗是“地震波速度”与“密度”的乘积（ρ*ν）。为了减少计算时间和混叠效应，数据通常会被平均到大于1英尺的采样间隔。 - 计算反射系数（R0) 通过使用声阻抗计算来确定每个界面处的值：
其中 ν1 和 ν2 分别为界面下方和上方各层的平均速度；ρ1 和 ρ2 是相应的平均密度。- 若井位未穿过火成岩体，则使用邻近已钻遇目标岩层的井获取相关参数（速度、密度， 等等。).
- 选择一个与邻近地震数据的振幅和相位谱相似的小波。
- 将合成子波与整个测井剖面的反射序列进行卷积，生成合成地震记录。最终的模拟地震记录 T(t) 可通过以下卷积模型进行描述：
其中 R0(tσ 是反射系数， w(t）是小波函数 n(t）是噪声。 - 若地震数据的频率在整个井段存在较大变化，则应使用具有不同相位和在不同深度区间主频不同的子波，重新计算合成地震记录。
- 如果合成记录与地震数据之间的匹配效果不理想，则重复该过程。
- 使用提供的软件进行计算（例如，Petrel E&P 软件平台
- 启动软件。选择 文件 | 打开项目 然后选择演示研究项目 tlm （用户可选择自己所需的项目）。项目应包含研究区域内的井数据、测井曲线、井口位置、地震数据体及解释层面。
- 点击 首页 | Windows | 2D 窗口 | 3D窗口 根据用户的偏好打开两个显示窗口以展示数据集。
- 在“输入窗格的 Wells 树”中，右键单击目标孔。打开 设置 小孔的窗口并选择 时间 按 Tab 键以创建新的时间日志。选择 速度函数，然后选择 DT 新时间日志中的数据。单击 好的 关闭设置窗口的按钮。一个新的单向时间日志将自动生成，并显示在“输入窗格的井树”中。
注意：单程时间录井是该井的时间-深度关系。有线录井数据域可转换为时间域，并在时间域窗口中显示。 - 激活现有 3D窗口 通过单击显示的窗口。如果没有 3D窗口 显示，创建一个新的 3D 窗口 通过点击 首页 | Windows | 3D窗口选择 TWT 在3D窗口的工具栏中选择以在时域中显示3D窗口。
- 选择具有代表性的测井曲线（如 'GR'、'DT' 或 'RT'） 威尔斯树 在3D窗口中显示它们；同时，在“输入”面板的“地震”树中选择地震剖面，以在同一3D窗口中显示。
- 使用 操作平面 使用 3D 窗口工具栏中的工具调整剖面位置，使其与井相交；用户将看到测井曲线已转换至时间域，并在同一 3D 窗口中与地震剖面一同显示。
- 点击 地震解释 | 地震测井关联 | 地震测井关联流程. 选择 综合地震井震联合 在“研究类型”行中选择相应类型，并在“井”行中添加目标井。在输入选项卡的 TDR 行中，选择校准的一维时间对数作为时间-深度关系；在地震行中选择地震数据体。在反射系数计算方法中任选一条测井数据。
- 点击 启动小波工具箱 创建一个Ricker子波用于此过程。单击 好的 将创建一个新的井段窗口和合成地震记录显示。
2. 将合成记录与实际地震反射层进行对比
- 使用平台中的自动化关联应用程序（如 Seismic Well Tie）将生成的合成地震道与地震剖面的垂直比例尺进行匹配。
- 调整合成地震图，以增强合成地震道与实际地震道中高振幅反射体的重叠程度。
- 反复调整合成地震图与实际地震道。当重叠程度达到最大时，解释人员即已实现合成地震图与实际地震道之间的“最佳匹配”。
- 重复该过程，直至关联结果达到所需水平。
- 使用提供的软件执行关联操作。
- 激活步骤 1.6.3 中创建的窗口，即从声波测井数据自动生成的单程时间测井。
注意：此自动生成的“单程时间测井”与实际地震反射界面并非完全对应。用户应校准单程时间测井与实际地震反射界面之间的关联关系。 - 为校准其关联关系，选择一个被井穿过且连续且具代表性的反射界面，然后手动调整测井深度。例如，要调整 DT 测井的深度，可在井树中右键单击 One-Way Time Log | 选择 Calculator 工具 | 然后在计算器工具的输入对话框中输入“DT=DT+10”，添加一个较小的时间增量（例如 10 ms）。
- 如果“10 ms”增量过大或过小，可在“Calculator”工具中将增量更改为其他时间值（可为负值）。反复检查测井与所选地震层位之间的关联性，并反复调整时间增量，直至关联性完全校准。
- 激活步骤 1.6.3 中创建的窗口，即从声波测井数据自动生成的单程时间测井。
3. 玄武岩岩床的提取
- 选择两个包围目标岩床的高振幅反射层。
注意：大多数侵入体在地震数据中表现为调谐反射包，即上、下侵入接触面的反射无法区分。当垂直侵入厚度介于 λ/4 与 λ/8 之间时（λ 为地震波长），会发生调谐现象19。因此，岩床在地震剖面中显示为一组强反射，其表观厚度是虚假的。 - 在对应于这两个高振幅反射层的层位之间提取探针。
注意：根据渲染技术的不同，存在多种工具可帮助解释人员更好地可视化目标体，例如“体块探针”、“层面探针”和“井探针”。然而，对于识别岩床与其围岩之间的接触关系，最佳工具是“层面探针”（“Surface probe, etc.” 是 Petrel 软件中的术语，软件用户应熟悉这些术语）。 - 通过调整体素连通性不透明度阈值，去除地质目标体周围区域。将默认阈值设为 20%。此处采用“不透明度渲染”的可视化方法，以显示玄武岩岩床提取结果（图 2C）。
注意：由于火成岩与沉积岩之间声阻抗差异显著，在二者界面处常出现高振幅反射。将低振幅部分设为透明，以突出火成体的形态。 - 由于隔离阈值可能高于 20–30%，需以较小增量逐步调整该值，以确保所有重要的火成体未被遗漏；阈值越大，越有可能丢失真实火成体的体积。
- 使用提供的软件执行操作：
- 点击 地震解释 面板，点击 插入层位探针。将在输入面板的地质体解释探针树中添加一个探针。双击所添加的层位探针，将弹出一个窗口。
- 在弹出窗口中点击 层位 选项卡，并选择两个圈定岩床区域的地震层面。点击 确定 以应用操作。
- 在输入面板中查看新添加的探针，三维窗口中将显示一个地震数据体。
- 双击该探针并选择 不透明度 选项卡。该选项卡将显示地震振幅的直方图。使用鼠标左键在直方图中绘制一条线，以控制地震数据体的不透明度。数据体中的低振幅部分应变为不可见，仅保留高振幅部分。
- 反复调整直方图，直至获得目标地质体的理想形态。
4. 取食通道的提取
- 在表层熔岩流下方不同深度处，选择连续且高能量的反射界面。
- 沿选定的地层界面进行时间切片，以识别与垂直通道相对应的不连续性。
- 反复调整双向时间（TWT），以获得导管不连续性的最佳成像效果。
注意：由于地震数据难以清晰成像垂直构造，因此需通过比较不同传播时间下的清晰度，从振幅体和方差体中选取更优的图像。 - 尝试不同的切片技术，然后选择能够更清晰显示不连续结构的方案。
注意：此处可使用不同的工具，例如方差体切片。其理论基础在于地震数据中每条地震剖面与其相邻地震道之间的相似性。另一种工具为方差体，它是通过对常规地震数据处理生成的一种新数据体，有助于识别构造与岩性变化、断层的平面组合等特征。 等等。20 - 将不同传播时间或深度的切片绘制到三维空间中。
- 使用提供的软件执行该操作。
- 双击 体积属性 在流程窗格的地球物理学分支中。检查 结构方法 在类别列中 方差 在属性列中，选择地震数据体输入到输入框，并在参数选项卡中调整其他参数。为获得更好的读取性能，勾选实现实列中的复选框。输入面板的地震树中将创建一个方差数据体。
- 右键单击方差立方体，然后单击 插入时间切片交点 在3D窗口中显示更多水平交点。使用 操控平面 使用3D窗口工具栏中的工具调整切片位置，以优化导管的显示效果。
- 右键单击地震振幅数据体并单击 插入时间切片交集 在3D窗口中显示更多水平交叉面。重复执行与步骤4.6.2相同的操作，以调整切片位置，优化导管的显示效果。
访问受限。请登录或开始试用以查看此内容。
结果
我们通过将上述技术应用于两种类型的火成岩体——水平岩床和垂直火山通道，来展示这些技术的实用性。采用不透明渲染技术提取岩床，而采用切片技术对火山通道进行解释。
提取火成岩床
在塔里木盆地北部英买2区，工业钻井已穿过多个岩床，但这些岩床的三维展布仍不清晰17。为了对岩床进行解释，我们对该区域的地震数据体进行了三维地震数据处理。首先，通过将合成地震记录与地震剖面进行对比，识别出地震数据体中与岩床相关的反射层位（图2A）；然后，在这些层位中插入表面测点（图2B），以约束岩床的横向展布范围；最后，采用透明度渲染技术从地震数据体中提取岩床的地质体（
访问受限。请登录或开始试用以查看此内容。
讨论
本文展示了两种用于描绘埋藏玄武岩火山管道系统形态与结构的方法：一种是不透明度渲染，另一种是时间切片。
不透明度渲染方法适用于与围岩地层具有连续且近水平界面的地质体。利用该方法，可提取岩浆凸起体的三维形态。通常情况下，流动方向应与岩浆凸起体的长轴方向一致。此外，界面处的反射系数（R0）需较高。若界面处的R0过低，则解释人员将无法将地表探针插入目标层位。例如，辉绿岩岩床的声波速度约为5500 m/s，而碳酸盐岩的速度与之相近，约为6,000 m/s12。因此，岩床与碳酸盐岩接触面处的反射系数过低，难以通过地表探针识别。使用该技术时，必须精确掌握目标岩石的速度信息。若速度数据缺失或估算不准确，该方法在地震数据体中的应用将受到极大限制。
时间切片法可应用于不具有连续且水平表面的地质体。当火成岩侵入体的声波速度与围岩相差较大（大多数情况下高于围岩）时，解释人员可利用时间切片技术来成像侵入体与周围岩石之间的边界。...
访问受限。请登录或开始试用以查看此内容。
披露
作者无任何利益冲突需要披露。
致谢
作者感谢国家自然科学基金对 WT（批准号：41272368）和 QKX（批准号：41630205）提供的经费支持。
访问受限。请登录或开始试用以查看此内容。
材料
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 佩特雷尔 E&P 软件平台 | Schlumberger | 软件版本：2014 |
参考文献
- Smallwood, J. R., Maresh, J. The properties, morphology and distribution of igneous sills: modelling, borehole data and 3D seismic from the Faroe-Shetland area. Geol. Soc. London Spec. Publ. 197 (1), 271-306 (2002).
- Millett, J. M., Hole, M. J., Jolley, D. W., Schofield, N., Campbell, E. Frontier exploration and the North Atlantic Igneous Province: new insights from a 2.6 km offshore volcanic sequence in the NE Faroe-Shetland Basin. J. Geol. Soc. 173 (2), 320-336 (2016).
- Lee, G. H., Kwon, Y. I., Yoon, C. S., Kim, H. J., Yoo, H. S. Igneous complexes in the eastern Northern South Yellow Sea Basin and their implications for hydrocarbon systems. Mar. Pet. Geol. 23 (6), 631-645 (2006).
- Rateau, R., Schofield, N., Smith, M. The potential role of igneous intrusions on hydrocarbon migration, West of Shetland. Pet. Geosci. 19 (3), 259-272 (2013).
- Magee, C., et al. Lateral magma flow in mafic sill complexes. Geosphere. 12 (3), 809-841 (2016).
- Magee, C., Jackson, C. A. L., Schofield, N. Diachronous sub-volcanic intrusion along deep-water margins: insights from the Irish Rockall Basin. Basin Res. 26 (1), 85-105 (2014).
- Symonds, P., Planke, S., Frey, O., Skogseid, J. Volcanic evolution of the Western Australian continental margin and its implications for basin development. The sedimentary basins of Western Australia. 2, 33-54 (1998).
- Thomson, K., Hutton, D. Geometry and growth of sill complexes: insights using 3D seismic from the North Rockall Trough. BVol. 66 (4), 364-375 (2004).
- Planke, S., Rasmussen, T., Rey, S., Myklebust, R. Petroleum Geology: North-West Europe and Global Perspectives-Proceedings of the 6th Petroleum Geology Conference. Doré, A. G., Vining, B. A. 6, Geological Society. London. 833-844 (2005).
- Magee, C., Hunt Stewart,, E,, Jackson, C. A. L. Volcano growth mechanisms and the role of sub-volcanic intrusions: Insights from 2D seismic reflection data. Earth Planet. Sci. Lett. 373, 41-53 (2013).
- Schofield, N. J., Brown, D. J., Magee, C., Stevenson, C. T. Sill morphology and comparison of brittle and non-brittle emplacement mechanisms. J. Geol. Soc. 169 (2), 127-141 (2012).
- Wang, L., Tian, W., Shi, Y. M., Guan, P. Volcanic structure of the Tarim flood basalt revealed through 3-D seismological imaging. Sci. Bull. 60 (16), 1448-1456 (2015).
- Sun, Q., et al. Neogene igneous intrusions in the northern South China Sea: Evidence from high-resolution three dimensional seismic data. Mar. Pet. Geol. 54, 83-95 (2014).
- Schofield, N., et al. Seismic imaging of 'broken bridges': linking seismic to outcrop-scale investigations of intrusive magma lobes. J. Geol. Soc. 169 (4), 421-426 (2012).
- Thomson, K. Volcanic features of the North Rockall Trough: application of visualisation techniques on 3D seismic reflection data. BVol. 67 (2), 116-128 (2005).
- Jackson, C. A. L. Seismic reflection imaging and controls on the preservation of ancient sill-fed magmatic vents. J. Geol. Soc. 169 (5), 503-506 (2012).
- Tian, W., et al. The Tarim picrite-basalt-rhyolite suite, a Permian flood basalt from northwest China with contrasting rhyolites produced by fractional crystallization and anatexis. CoMP. 160 (3), 407-425 (2010).
- Chen, M. -M., et al. Peridotite and pyroxenite xenoliths from Tarim, NW China: Evidences for melt depletion and mantle refertilization in the mantle source region of the Tarim flood basalt. Lithos. 204, 97-111 (2014).
- Magee, C., Maharaj, S. M., Wrona, T., Jackson, C. A. L. Controls on the expression of igneous intrusions in seismic reflection data. Geosphere. 11 (4), 1024-1041 (2015).
- Bahorich, M., Farmer, S. 3-D seismic discontinuity for faults and stratigraphic features: The coherence cube. The Leading Edge. 14 (10), 1053-1058 (1995).
访问受限。请登录或开始试用以查看此内容。
重印与许可
