CN111142159A - 中深层异常体地震速度建模方法及系统 - Google Patents
中深层异常体地震速度建模方法及系统 Download PDFInfo
- Publication number
- CN111142159A CN111142159A CN201811301562.4A CN201811301562A CN111142159A CN 111142159 A CN111142159 A CN 111142159A CN 201811301562 A CN201811301562 A CN 201811301562A CN 111142159 A CN111142159 A CN 111142159A
- Authority
- CN
- China
- Prior art keywords
- speed
- model
- depth domain
- abnormal body
- velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000002159 abnormal effect Effects 0.000 title claims abstract description 132
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000003384 imaging method Methods 0.000 claims abstract description 112
- 238000013508 migration Methods 0.000 claims abstract description 22
- 230000005012 migration Effects 0.000 claims abstract description 22
- 230000008859 change Effects 0.000 claims description 26
- 239000000203 mixture Substances 0.000 claims description 9
- 238000005056 compaction Methods 0.000 claims description 8
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 8
- 239000011707 mineral Substances 0.000 claims description 8
- 239000011435 rock Substances 0.000 claims description 8
- 238000003325 tomography Methods 0.000 claims description 8
- 230000000877 morphologic effect Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 abstract description 30
- 238000004458 analytical method Methods 0.000 description 20
- 230000000694 effects Effects 0.000 description 8
- 235000010755 mineral Nutrition 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000012937 correction Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 235000002639 sodium chloride Nutrition 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000002547 anomalous effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000007598 dipping method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002674 ointment Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
公开了一种中深层异常体地震速度建模方法及系统。该方法可以包括：步骤1：根据原始模型获得初始深度域速度模型；步骤2：根据初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集；步骤3：根据深度域CRP道集，拾取剩余速度，进而获得剩余延迟；步骤4：判断成像剖面是否存在异常体，若否，获得最终的深度域速度模型；若是，以更新的深度域速度模型为初始深度域速度模型，重复步骤2‑4。本发明根据高信噪比区拾取深层与低信噪比区的剩余速度，判断异常体并填充，获得成像剖面清晰、速度谱剩余能量小、道集平整、沿层速度谱剩余小、数据可信度高的深度域速度模型。
Description
技术领域
本发明涉及油汽地球物理技术领域，更具体地，涉及一种中深层异常体地震速度建模方法及系统。
背景技术
精确地震速度的求取是地球物理的核心问题之一，无论是对地震资料偏移、时深转换、静校正，还是对岩性、岩相和储层预测，均起着重要作用。传统的深度域层速度方法是根据叠加速度求均方根速度，再由DIX公式计算层速度，其精度较低，特别是对倾斜界面，即使做倾角校正，其精度也很难满足实际需要。常规速度分析已经可以满足常规处理的需要，而无法解决复杂构造的速度分析问题。近些年发展的各种地震成像方法和速度分析方法主要解决的是地下复杂的地质构造，针对更加符合实际的地下地质情况，解决实际问题。产生了如发展的常规叠加速度分析用于解决倾斜地层的倾角时差校正，高阶动校正优化来解决大偏移距问题，以及层析技术等。到目前为止，1990年，Deregowski(1990)提出的叠前时间域环速度分析方法；1987年到1989年，Al-yahya提出了剩余曲率叠前深度偏移速度分析；1991年到1992年，Mackay和Abma深度聚焦分析叠前深度偏移速度分析；1997年，Berkhout基于双聚焦思想提出的CFP道集偏移速度分析等方法都是叠前偏移速度分析方法中经典的方法。现在很多叠前速度分析方法都或多或少的有这些思想的印记。
不同的速度分析方法速度准确性的判别准则是不同的，针对于主流的叠加速度分析方法来说，主要是以共中心点道集的同相轴是否拉平作为速度正确与否的判别准则。对于偏移速度分析方法来说，一般通过偏移成像的结果的最优化可以判断用于偏移的速度场是否正确。它是将偏移成像和速度分析结合起来，根据偏移成像的结果的好坏来对初始的速度场模型进行速度的更新，针对不同的偏移速度分析方法，有不同的判别准则，目前主流的地震速度建立方法主要为以下三类，具体介绍如下：
(1)1986年，Faye和Jeannot提出分别通过零时间和零偏移距成像条件得到的成像深度是否相等来判断用于偏移的速度场是否正确。他们提出，用于成像的偏移速度场正确时，在地下成像点位置处零时间成像深度点和零偏移距成像聚焦点位置相同，否则不一致。用于偏移成像的速度较高时，成像深度点大于聚焦深度点，否则，则相反。
(2)1989年，Al-Yahya提出了CIG拉平准则。该速度分析方法的判别准则是如果用于成像的速度场模型正确，则从偏移成像结果中提取出的CIG是拉平的，否则偏移成像中提取的CIG道集表现为向下或者向上弯曲，可以在该同相轴上提取剩余曲率，然后用于更新速度。当用于偏移的速度大于真实速度时，CIG向下弯曲，当用于偏移的速度小于真实速度时，CIG向上弯曲。
(3)基于CFP道集的等时原理。1997年，Berkhout提出的基于双聚焦理论的叠前深度偏移速度分析方法，速度正确与否的判别准则是使用等(旅行)时原理。该准则认为在用于偏移的速度场正确时，聚焦点响应与其相应的逆时聚焦算子两者的时差为零，差异时移面板上道集的同相轴是拉平的。否则，则不是拉平的。当速度场偏大时，DTS面板上的道集表现为向下弯曲，当速度场偏小时，DTS面板表现为向上弯曲。根据DTS面板上的非零时差进行速度分析，这点和基于剩余曲率的偏移速度分析方法相似。
目前科研领域主要是基于方法(3)进行探讨，但此种方法虽然精度略高，李璐更先进，但十分耗时，且对应的偏移方法对速度模型要求很高；而工业界尤其大部分商业软件，主要基于方法(2)进行速度分析和建模，此种方法简单易行，实用效率高，且灵活方便，但是没有针对异常地质体的建模方法。因此，有必要开发一种中深层异常体地震速度建模方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种中深层异常体地震速度建模方法及系统，其能够根据高信噪比区拾取深层与低信噪比区的剩余速度，判断异常体并填充，获得成像剖面清晰、速度谱剩余能量小、道集平整、沿层速度谱剩余小、数据可信度高的深度域速度模型。
根据本发明的一方面，提出了一种中深层异常体地震速度建模方法。所述方法可以包括：步骤1：根据原始模型获得初始深度域速度模型；步骤2：根据所述初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集；步骤3：根据所述深度域CRP道集，拾取剩余速度，进而获得剩余延迟，其中，针对深层与低信噪比区域，沿同层的高信噪比区域的趋势进行拾取，保证同层速度平滑；步骤4：根据所述剩余延迟，进行层析成像，获得成像剖面，判断所述成像剖面是否存在异常体，若否，则根据所述成像剖面获得最终的深度域速度模型；若是，则对所述异常体进行填充，获得更新的深度域速度模型，以所述更新的深度域速度模型为初始深度域速度模型，重复步骤2-4，获得最终的深度域速度模型。
优选地，根据原始模型获得初始深度域速度模型包括：针对所述原始模型设定形态权重因子和速度变化权重因子；根据所述形态权重因子和所述速度变化权重因子，对所述原始模型进行转化，获得所述初始深度域速度模型。
优选地，所述原始模型为RMS均方根速度模型。
优选地，所述形态权重因子和所述速度变化权重因子的关系为：
A2>2A1 (1)
其中，A1为形态权重因子，A2为速度变化权重因子。
优选地，判断所述成像剖面是否存在异常体包括：判断所述成像剖面的清晰度和能量聚焦程度是否符合预设要求；判断所述成像剖面的道集拉平程度和剩余能量是否符合预设要求；判断所述成像剖面的剩余速度能量是否符合预设要求；判断所述成像剖面的速度与井上速度的对应性是否符合预设要求；若上述判断均为是，则所述成像剖面不存在异常体，否则，所述成像剖面不存在异常体。
优选地，对所述异常体进行填充包括：根据所述异常体的尺度将异常体分为小尺度异常体与大尺度异常体；针对所述小尺度异常体，通过与小尺度异常体同层相近的速度进行填充；针对所述大尺度异常体，若存在构造相同或相近的已知井，则通过所述已知井的速度进行填充。
优选地，针对所述大尺度异常体进行填充还包括：若不存在构造相同或相近的已知井，则根据地质认识，通过岩石矿物成分构成、压实程度与同层速度进行填充。
优选地，针对所述大尺度异常体进行填充还包括：若不存在地质认识，则以同层相近的速度进行填充。
根据本发明的另一方面，提出了一种中深层异常体地震速度建模系统，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：步骤1：根据原始模型获得初始深度域速度模型；步骤2：根据所述初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集；步骤3：根据所述深度域CRP道集，拾取剩余速度，进而获得剩余延迟，其中，针对深层与低信噪比区域，沿同层的高信噪比区域的趋势进行拾取，保证同层速度平滑；步骤4：根据所述剩余延迟，进行层析成像，获得成像剖面，判断所述成像剖面是否存在异常体，若否，则根据所述成像剖面获得最终的深度域速度模型；若是，则对所述异常体进行填充，获得更新的深度域速度模型，以所述更新的深度域速度模型为初始深度域速度模型，重复步骤2-4，获得最终的深度域速度模型。
优选地，根据原始模型获得初始深度域速度模型包括：针对所述原始模型设定形态权重因子和速度变化权重因子；根据所述形态权重因子和所述速度变化权重因子，对所述原始模型进行转化，获得所述初始深度域速度模型。
优选地，所述原始模型为RMS均方根速度模型。
优选地，所述形态权重因子和所述速度变化权重因子的关系为：
A2>2A1 (1)
其中，A1为形态权重因子，A2为速度变化权重因子。
优选地，判断所述成像剖面是否存在异常体包括：判断所述成像剖面的清晰度和能量聚焦程度是否符合预设要求；判断所述成像剖面的道集拉平程度和剩余能量是否符合预设要求；判断所述成像剖面的剩余速度能量是否符合预设要求；判断所述成像剖面的速度与井上速度的对应性是否符合预设要求；若上述判断均为是，则所述成像剖面不存在异常体，否则，所述成像剖面不存在异常体。
优选地，对所述异常体进行填充包括：根据所述异常体的尺度将异常体分为小尺度异常体与大尺度异常体；针对所述小尺度异常体，通过与小尺度异常体同层相近的速度进行填充；针对所述大尺度异常体，若存在构造相同或相近的已知井，则通过所述已知井的速度进行填充。
优选地，针对所述大尺度异常体进行填充还包括：若不存在构造相同或相近的已知井，则根据地质认识，通过岩石矿物成分构成、压实程度与同层速度进行填充。
优选地，针对所述大尺度异常体进行填充还包括：若不存在地质认识，则以同层相近的速度进行填充。
本发明具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的中深层异常体地震速度建模方法的步骤的流程图。
图2示出了根据本发明的一个实施例的成像剖面的示意图。
图3示出了根据本发明的一个实施例的速度谱剩余能量谱的示意图。
图4示出了根据本发明的一个实施例的深度域CRP道集的示意图。
图5示出了根据本发明的一个实施例的沿层剩余速度能量谱的示意图。
图6示出了根据本发明的一个实施例的初始深度域速度模型的成像剖面的示意图。
图7示出了根据本发明的一个实施例的初始深度域速度模型的一个CMP点的垂直剩余速度谱的示意图。
图8示出了根据本发明的一个实施例的根据初始深度域速度模型获得的CRP道集的示意图。
图9示出了根据本发明的一个实施例的初始深度域速度模型的沿层速度谱的示意图。
图10示出了根据本发明的一个实施例的最终深度域速度模型的沿层速度谱的示意图。
图11示出了根据本发明的一个实施例的最终深度域速度模型的成像剖面的示意图。
图12示出了根据本发明的一个实施例的最终深度域速度模型的速度谱的示意图。
图13示出了根据本发明的一个实施例的根据最终深度域速度模型获得的CRP道集的示意图。
图14a、14b分别示出了根据本发明的一个实施例的膏盐填充速度为4500m/s的速度模型与对应的地震成像的示意图。
图15a、15b分别示出了根据本发明的一个实施例的膏盐填充速度为6000m/s的速度模型与对应的地震成像的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的中深层异常体地震速度建模方法的步骤的流程图。
在该实施例中，根据本发明的中深层异常体地震速度建模方法可以包括：步骤1：根据原始模型获得初始深度域速度模型；步骤2：根据初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集；步骤3：根据深度域CRP道集，拾取剩余速度，进而获得剩余延迟，其中，针对深层与低信噪比区域，沿同层的高信噪比区域的趋势进行拾取，保证同层速度平滑；步骤4：根据剩余延迟，进行层析成像，获得成像剖面，判断成像剖面是否存在异常体，若否，则根据成像剖面获得最终的深度域速度模型；若是，则对异常体进行填充，获得更新的深度域速度模型，以更新的深度域速度模型为初始深度域速度模型，重复步骤2-4，获得最终的深度域速度模型。
在一个示例中，根据原始模型获得初始深度域速度模型包括：针对原始模型设定形态权重因子和速度变化权重因子；根据形态权重因子和速度变化权重因子，对原始模型进行转化，获得初始深度域速度模型。
在一个示例中，原始模型为RMS均方根速度模型。
在一个示例中，形态权重因子和速度变化权重因子的关系为：
A2>2A1 (1)
其中，A1为形态权重因子，A2为速度变化权重因子。
在一个示例中，判断成像剖面是否存在异常体包括：判断成像剖面的清晰度和能量聚焦程度是否符合预设要求；判断成像剖面的道集拉平程度和剩余能量是否符合预设要求；判断成像剖面的剩余速度能量是否符合预设要求；判断成像剖面的速度与井上速度的对应性是否符合预设要求；若上述判断均为是，则成像剖面不存在异常体，否则，成像剖面不存在异常体。
在一个示例中，对异常体进行填充包括：根据异常体的尺度将异常体分为小尺度异常体与大尺度异常体；针对小尺度异常体，通过与小尺度异常体同层相近的速度进行填充；针对大尺度异常体，若存在构造相同或相近的已知井，则通过已知井的速度进行填充。
在一个示例中，针对大尺度异常体进行填充还包括：若不存在构造相同或相近的已知井，则根据地质认识，通过岩石矿物成分构成、压实程度与同层速度进行填充。
在一个示例中，针对大尺度异常体进行填充还包括：若不存在地质认识，则以同层相近的速度进行填充。
具体地，根据本发明的中深层异常体地震速度建模方法可以包括：
步骤1：针对原始模型设定形态权重因子和速度变化权重因子；根据形态权重因子和速度变化权重因子，对原始模型进行转化，获得初始深度域速度模型，其中，原始模型为RMS均方根速度模型，形态权重因子和速度变化权重因子的关系为公式(1)。
步骤2：根据初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集。
图5示出了根据本发明的一个实施例的沿层剩余速度能量谱的示意图。
步骤3：根据深度域CRP道集，拾取剩余速度，进而获得剩余延迟，如图5中横坐标CMP(Common Middle Point)对应的2100至2400点段，沿层速度谱的纵坐标轴为剩余能量值，从图中可以发现2100至2400点段为低信噪比段，因为剩余能量表征为能量十分不集中的散点段，几乎无法拾取，传统技术的对应策略基本为将剩余能量视为0值对待或者拾取对应点的最大值，而实际上，此两种方法前者仅保留了初始速度，而与同层修正后速度差异很大，后者使得速度跳动很大且不保证速度正确。本方法针对深层与低信噪比区域，沿同层的高信噪比区域的趋势进行拾取，保证同层速度平滑，符合客观的地震沉积学基本原理。
图2示出了根据本发明的一个实施例的成像剖面的示意图。
图3示出了根据本发明的一个实施例的速度谱剩余能量谱的示意图。
图4示出了根据本发明的一个实施例的深度域CRP道集的示意图。
步骤4：根据剩余延迟，进行层析成像，获得成像剖面，判断成像剖面是否存在异常体，通过四要素来判断，判断成像剖面的清晰度和能量聚焦程度是否符合预设要求，如图2所示；判断成像剖面的道集拉平程度和剩余能量是否符合预设要求，其中，道集不发生较大畸变即可，因为在低信噪比区域对应的道集点基本没有可见和可信的水平同相轴，如图4所示；判断成像剖面的剩余速度能量是否符合预设要求，其中，每个单CMP点的速度谱剩余能量趋势变小，沿层速度谱的剩余能量趋势变小，低信噪比区的剩余速度能量趋势(沿剩余能量散点的趋势线)向0值方向变小，而不是和高信噪比区域一样，要求值为0，因为低信噪比区的剩余能量不聚焦，如图3所示；判断成像剖面的速度与井上速度的对应性是否符合预设要求，即与井上对应速度基本对应，与同层高信噪比的速度大小和变化在趋势上相对应，如图5所示；若上述判断均为是，则成像剖面不存在异常体，否则，成像剖面不存在异常体，其中，本领域技术人员可以根据具体情况设定各个预设要求。
若成像剖面不存在异常体，则根据成像剖面获得最终的深度域速度模型；若成像剖面存在异常体，则对异常体进行填充，根据异常体的尺度将异常体分为小尺度异常体与大尺度异常体，小尺度异常体的尺度在半个地震波长以及该尺度以下，大尺度异常体的尺度大于地震波长尺度，尤其大于3倍波长且大于纵向采样20倍以上；针对小尺度异常体，因为速度本身相对其他地震数据而言是低频量，小尺度异常体的异常速度对成像影响较小，通过与小尺度异常体同层相近的速度进行填充；针对大尺度异常体需要进行单独匹配速度的填充，若存在构造相同或相近的已知井，则通过已知井的速度进行填充，若不存在构造相同或相近的已知井，则根据地质认识，通过岩石矿物成分构成、压实程度与同层速度进行填充，若不存在地质认识，则以同层相近的速度进行填充，获得更新的深度域速度模型。
常规速度建模时，不考虑地质认识，仅按照常规流程，在异常体对应的层位处的层速度直接作为异常体速度后进行优化，但常规流程此处层位速度一般与横向上邻近点的层速度大小相近，速度变化平缓，体现不出异常体的存在，以这样的层速度为基础进行优化一方面迭代轮次较多，一方面最终也很难优化得到较为精确的速度结果。此时借鉴相应地质认识，对异常体单独的进行速度填充。一般工区在地震速度建模前就对异常体有一定的地质认识，包括异常体的岩性构成、大致深度范围、大致厚度情况、临区相应异常体的测井速度等等，进而对异常体的大致速度有个最大最小值的估算。结合这些已有认识，采用速度可能的极值进行填充，用不同极值填充后的速度进行成像，对成像后的剖面和剩余能量进行对比，进而选择效果较好的对应速度值，以此速度值去填充异常体，再以此异常体速度为基础去优化，便能更快更好的得到最终迭代的速度结果。
现有技术中，大尺度异常体的速度填充一般采用人工主观设定上下限值进行参数实验，例如上下限速度一般为上下层速度极值，再通过偏移进行成像效果对比和道集拉平效果对比以及速度剩余能量对比。这种参数实验一方面测试时间成本高，如上下限速度为4000m/s到7000m/s，间隔250m/s，需要进行12轮耗时的偏移流程和成果对比，且如果异常体速度不在上下限速度值范围内，则直接导致速度填充的错误结果。通过本方法进行异常体填充，一方面不需要多轮参数实验，从而节省时间成本，另一方面打破人为的主观的上下限速度测试限制，符合实际地质沉积基本客观规律，保证速度结果的准确可靠。
以更新的深度域速度模型为初始深度域速度模型，重复步骤2-4，获得最终的深度域速度模型。
本方法根据高信噪比区拾取深层与低信噪比区的剩余速度，判断异常体并填充，获得成像剖面清晰、速度谱剩余能量小、道集平整、沿层速度谱剩余小、数据可信度高的深度域速度模型。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
根据本发明的中深层异常体地震速度建模方法可以包括：
步骤1：针对原始模型设定形态权重因子和速度变化权重因子；根据形态权重因子和速度变化权重因子，对原始模型进行转化，获得初始深度域速度模型，其中，原始模型为RMS均方根速度模型，形态权重因子和速度变化权重因子的关系为公式(1)。
步骤2：根据初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集。
步骤3：根据深度域CRP道集，拾取剩余速度，进而获得剩余延迟，其中，针对深层与低信噪比区域，沿同层的高信噪比区域的趋势进行拾取，保证同层速度平滑。
步骤4：根据剩余延迟，进行层析成像，获得成像剖面，判断成像剖面是否存在异常体，判断成像剖面的清晰度和能量聚焦程度是否符合预设要求；判断成像剖面的道集拉平程度和剩余能量是否符合预设要求，其中，道集不发生较大畸变即可，因为在低信噪比区域对应的道集点基本没有可见和可信的水平同相轴；判断成像剖面的剩余速度能量是否符合预设要求，其中，每个单CMP点的速度谱剩余能量趋势变小，沿层速度谱剩余能量趋势变小，低信噪比区的剩余速度能量趋势(沿剩余能量散点的趋势线)向0值方向变小，而不是和高信噪比区域一样，要求值为0，因为低信噪比区的剩余能量不聚焦；判断成像剖面的速度与井上速度的对应性是否符合预设要求，即与井上对应速度基本对应，与同层高信噪比的速度大小和变化在趋势上相对应；若上述判断均为是，则成像剖面不存在异常体，否则，成像剖面不存在异常体，其中，本领域技术人员可以根据具体情况设定各个预设要求。
若成像剖面不存在异常体，则根据成像剖面获得最终的深度域速度模型；若成像剖面存在异常体，则对异常体进行填充，根据异常体的尺度将异常体分为小尺度异常体与大尺度异常体，小尺度异常体的尺度在半个地震波长以及该尺度以下，大尺度异常体的尺度大于地震波长尺度，尤其大于3倍波长且大于纵向采样20倍以上；针对小尺度异常体，因为速度本身相对其他地震数据而言是低频量，小尺度异常体的异常速度对成像影响较小，通过与小尺度异常体同层相近的速度进行填充；针对大尺度异常体需要进行单独匹配速度的填充，若存在构造相同或相近的已知井，则通过已知井的速度进行填充，若不存在构造相同或相近的已知井，则根据地质认识，通过岩石矿物成分构成、压实程度与同层速度进行填充，若不存在地质认识，则以同层相近的速度进行填充，获得更新的深度域速度模型。以更新的深度域速度模型为初始深度域速度模型，重复步骤2-4，获得最终的深度域速度模型。
图6示出了根据本发明的一个实施例的初始深度域速度模型的成像剖面的示意图。
图7示出了根据本发明的一个实施例的初始深度域速度模型的一个CMP点的垂直剩余速度谱的示意图。
图8示出了根据本发明的一个实施例的根据初始深度域速度模型获得的CRP道集的示意图。
图9示出了根据本发明的一个实施例的初始深度域速度模型的沿层速度谱的示意图。
图10示出了根据本发明的一个实施例的最终深度域速度模型的沿层速度谱的示意图，与图9对比，可以发现剩余能量更加平直的集中在了剩余能量为0值的附近，代表剩余能量更小，有效能量更加准确，对应的速度更加准确。
图11示出了根据本发明的一个实施例的最终深度域速度模型的成像剖面的示意图，与图6相比，地震成像层位更加清晰，归位深度更加准确。
图12示出了根据本发明的一个实施例的最终深度域速度模型的速度谱的示意图，与图7相比，同一CMP点对应的剩余速度能量更加平直的集中在了剩余能量为0值的附近，代表剩余能量更小，有效能量更加准确，对应的速度更加准确。
图13示出了根据本发明的一个实施例的根据最终深度域速度模型获得的CRP道集的示意图，与图8相比，同一CMP点的CRP成像道集同相轴更加平整，表明速度更加准确。
如图6和图11所示，从剖面中可以发现，本例中，5000米上下存在一个明显地质层位，这里称为No1层，6000米上下也存在一地质层位，这里称为No2层。No2层位目标层位，信噪比较差。且其存在圈闭幅度高低以及该圈闭的真假难以判断。计划在此圈闭打井，成本较高。难点还在于：
(1)发育推覆构造，断裂系统很发育。目标圈闭幅度较小约200米，且可能受上层盐下速度不准的影响，造成成像失真，如果该圈闭不存在，造成的损失将十分巨大，因此，希望进行高精度沿层深度速度建模，再进行偏移，进一步确定真实构造情况；
(2)发育膏岩盐，难以准确成像，速度难于确定；
(3)目的层信噪比低，难于基于数据的速度建模。
由No1层和No2层的沿层速度谱可以发现，No1层的高信噪比段较多，完成高信噪此段建模，通过图9和图10的对比可以发现此段建模效果良好，优化后的剩余速度能量基本为0，且其他要求均符合四要素控制的原则和方法。而该段的低信噪比段为本次任务的关键段，即，地段信噪比很低，没有集中地剩余速度能量可以拾取，但对下层圈闭构造的真伪起着至关重要的决定性作用。采用沿同层的高信噪比区域的趋势进行拾取后，其成像效果符合本发明的四个要素要求的原则标准和方法，且能量更加集中，且趋近于0值了。而目标层No2段，在优化前基本全段没有集中地剩余能量，基本无法进行常规速度建模。采用由浅入深、逐层迭代、低信噪比四要素控制的方法原则后，出现了可以参考的剩余能量，数据质量发生的本质上的变化。即成像效果得到了很大改进，数据可信度得到了极大的提升，更加有效的指导亿元级别井位的部署。
图14a、14b分别示出了根据本发明的一个实施例的膏盐填充速度为4500m/s的速度模型与对应的地震成像的示意图。
图15a、15b分别示出了根据本发明的一个实施例的膏盐填充速度为6000m/s的速度模型与对应的地震成像的示意图。
还需要注意的是，针对上层段大规模推覆构造，尤其目的层上方的膏盐岩发育情况，需要进行速度填充实验优化。通过异常体填充，因为2000米范围外存在一井位钻穿该目的层段，测得井点处的速度值为5800m/s，且已有地质认识为，该膏盐填充矿物组分构成压实程度对应速度范围在4500m/s至6000m/s之间。因此结合井位和地质认识，该区该深度进行了速度填充试验，分别填充较为极端的4500m/s和6000m/s进行试验验证。如图14a、图14b、图15a、图15b所示。
由填充试验可以发现，图14b的白框中圈出的膏盐构造成像效果不如图15b，后者膏盐边界更加清晰，地层更加平滑，符合深层地质认识，也更加贴近同层附近井数据的对应速度。符合本发明的标准和方法原则，因此可以在6000m/s速度填充膏盐的基础上再进一步优化。
综上所述，本发明根据高信噪比区拾取深层与低信噪比区的剩余速度，判断异常体并填充，获得成像剖面清晰、速度谱剩余能量小、道集平整、沿层速度谱剩余小、数据可信度高的深度域速度模型。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的中深层异常体地震速度建模系统，其上存储有计算机程序，其中，程序被处理器执行时实现以下步骤：步骤1：根据原始模型获得初始深度域速度模型；步骤2：根据初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集；步骤3：根据深度域CRP道集，拾取剩余速度，进而获得剩余延迟，其中，针对深层与低信噪比区域，沿同层的高信噪比区域的趋势进行拾取，保证同层速度平滑；步骤4：根据剩余延迟，进行层析成像，获得成像剖面，判断成像剖面是否存在异常体，若否，则根据成像剖面获得最终的深度域速度模型；若是，则对异常体进行填充，获得更新的深度域速度模型，以更新的深度域速度模型为初始深度域速度模型，重复步骤2-4，获得最终的深度域速度模型。
在一个示例中，根据原始模型获得初始深度域速度模型包括：针对原始模型设定形态权重因子和速度变化权重因子；根据形态权重因子和速度变化权重因子，对原始模型进行转化，获得初始深度域速度模型。
在一个示例中，原始模型为RMS均方根速度模型。
在一个示例中，形态权重因子和速度变化权重因子的关系为：
A2>2A1 (1)
其中，A1为形态权重因子，A2为速度变化权重因子。
在一个示例中，判断成像剖面是否存在异常体包括：判断成像剖面的清晰度和能量聚焦程度是否符合预设要求；判断成像剖面的道集拉平程度和剩余能量是否符合预设要求；判断成像剖面的剩余速度能量是否符合预设要求；判断成像剖面的速度与井上速度的对应性是否符合预设要求；若上述判断均为是，则成像剖面不存在异常体，否则，成像剖面不存在异常体。
在一个示例中，对异常体进行填充包括：根据异常体的尺度将异常体分为小尺度异常体与大尺度异常体；针对小尺度异常体，通过与小尺度异常体同层相近的速度进行填充；针对大尺度异常体，若存在构造相同或相近的已知井，则通过已知井的速度进行填充。
在一个示例中，针对大尺度异常体进行填充还包括：若不存在构造相同或相近的已知井，则根据地质认识，通过岩石矿物成分构成、压实程度与同层速度进行填充。
在一个示例中，针对大尺度异常体进行填充还包括：若不存在地质认识，则以同层相近的速度进行填充。
本系统根据高信噪比区拾取深层与低信噪比区的剩余速度，判断异常体并填充，获得成像剖面清晰、速度谱剩余能量小、道集平整、沿层速度谱剩余小、数据可信度高的深度域速度模型。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种中深层异常体地震速度建模方法，包括：
步骤1：根据原始模型获得初始深度域速度模型；
步骤2：根据所述初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集；
步骤3：根据所述深度域CRP道集，拾取剩余速度，进而获得剩余延迟，其中，针对深层与低信噪比区域，沿同层的高信噪比区域的趋势进行拾取，保证同层速度平滑；
步骤4：根据所述剩余延迟，进行层析成像，获得成像剖面，判断所述成像剖面是否存在异常体，若否，则根据所述成像剖面获得最终的深度域速度模型；若是，则对所述异常体进行填充，获得更新的深度域速度模型，以所述更新的深度域速度模型为初始深度域速度模型，重复步骤2-4，获得最终的深度域速度模型。
2.根据权利要求1所述的中深层异常体地震速度建模方法，其中，根据原始模型获得初始深度域速度模型包括：
针对所述原始模型设定形态权重因子和速度变化权重因子；
根据所述形态权重因子和所述速度变化权重因子，对所述原始模型进行转化，获得所述初始深度域速度模型。
3.根据权利要求2所述的中深层异常体地震速度建模方法，其中，所述原始模型为RMS均方根速度模型。
4.根据权利要求2所述的中深层异常体地震速度建模方法，其中，所述形态权重因子和所述速度变化权重因子的关系为：
A2>2A1 (1)
其中，A1为形态权重因子，A2为速度变化权重因子。
5.根据权利要求1所述的中深层异常体地震速度建模方法，其中，判断所述成像剖面是否存在异常体包括：
判断所述成像剖面的清晰度和能量聚焦程度是否符合预设要求；
判断所述成像剖面的道集拉平程度和剩余能量是否符合预设要求；
判断所述成像剖面的剩余速度能量是否符合预设要求；
判断所述成像剖面的速度与井上速度的对应性是否符合预设要求；
若上述判断均为是，则所述成像剖面不存在异常体，否则，所述成像剖面不存在异常体。
6.根据权利要求1所述的中深层异常体地震速度建模方法，其中，对所述异常体进行填充包括：
根据所述异常体的尺度将异常体分为小尺度异常体与大尺度异常体；
针对所述小尺度异常体，通过与小尺度异常体同层相近的速度进行填充；
针对所述大尺度异常体，若存在构造相同或相近的已知井，则通过所述已知井的速度进行填充。
7.根据权利要求6所述的中深层异常体地震速度建模方法，其中，针对所述大尺度异常体进行填充还包括：
若不存在构造相同或相近的已知井，则根据地质认识，通过岩石矿物成分构成、压实程度与同层速度进行填充。
8.根据权利要求7所述的中深层异常体地震速度建模方法，其中，针对所述大尺度异常体进行填充还包括：
若不存在地质认识，则以同层相近的速度进行填充。
9.一种中深层异常体地震速度建模系统，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现以下步骤：
步骤1：根据原始模型获得初始深度域速度模型；
步骤2：根据所述初始深度域速度模型与共中心点道集，进行目标线叠前深度偏移，获得深度域CRP道集；
步骤3：根据所述深度域CRP道集，拾取剩余速度，进而获得剩余延迟，其中，针对深层与低信噪比区域，沿同层的高信噪比区域的趋势进行拾取，保证同层速度平滑；
步骤4：根据所述剩余延迟，进行层析成像，获得成像剖面，判断所述成像剖面是否存在异常体，若否，则根据所述成像剖面获得最终的深度域速度模型；若否，则对所述异常体进行填充，获得更新的深度域速度模型，以所述更新的深度域速度模型为初始深度域速度模型，重复步骤2-4。
10.根据权利要求9所述的中深层异常体地震速度建模系统，其中，根据原始模型获得初始深度域速度模型包括：
针对所述原始模型设定形态权重因子和速度变化权重因子；
根据所述形态权重因子和所述速度变化权重因子，对所述原始模型进行转化，获得所述初始深度域速度模型；
其中，所述形态权重因子和所述速度变化权重因子的关系为：
A2>2A1 (1)
其中，A1为形态权重因子，A2为速度变化权重因子。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811301562.4A CN111142159B (zh) | 2018-11-02 | 2018-11-02 | 中深层异常体地震速度建模方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811301562.4A CN111142159B (zh) | 2018-11-02 | 2018-11-02 | 中深层异常体地震速度建模方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111142159A true CN111142159A (zh) | 2020-05-12 |
| CN111142159B CN111142159B (zh) | 2022-04-08 |
Family
ID=70516267
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811301562.4A Active CN111142159B (zh) | 2018-11-02 | 2018-11-02 | 中深层异常体地震速度建模方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111142159B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113960658A (zh) * | 2020-07-20 | 2022-01-21 | 中国石油天然气股份有限公司 | 一种基于地质地震模型的测井约束速度建模方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102879814A (zh) * | 2011-07-15 | 2013-01-16 | 中国石油天然气集团公司 | 一种准确的深度域层速度更新方法 |
| CN103544185A (zh) * | 2012-07-16 | 2014-01-29 | 中国石油化工股份有限公司 | 一种测井数据文件存储方法 |
| CN104536043A (zh) * | 2014-12-26 | 2015-04-22 | 中国石油天然气股份有限公司 | 一种深度域整体速度模型融合方法及装置 |
| US20150177401A1 (en) * | 2013-12-19 | 2015-06-25 | Chevron U.S.A. Inc. | Method of Correcting Velocity for Complex Surface Topography |
| CN104950330A (zh) * | 2014-03-25 | 2015-09-30 | 中国石油化工股份有限公司 | 气顶油藏深度域成像的速度建模方法及系统 |
| CN106597533A (zh) * | 2016-11-17 | 2017-04-26 | 中国石油化工股份有限公司 | 一种用于山前带地震资料处理的深度域速度建模方法 |
-
2018
- 2018-11-02 CN CN201811301562.4A patent/CN111142159B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102879814A (zh) * | 2011-07-15 | 2013-01-16 | 中国石油天然气集团公司 | 一种准确的深度域层速度更新方法 |
| CN103544185A (zh) * | 2012-07-16 | 2014-01-29 | 中国石油化工股份有限公司 | 一种测井数据文件存储方法 |
| US20150177401A1 (en) * | 2013-12-19 | 2015-06-25 | Chevron U.S.A. Inc. | Method of Correcting Velocity for Complex Surface Topography |
| CN104950330A (zh) * | 2014-03-25 | 2015-09-30 | 中国石油化工股份有限公司 | 气顶油藏深度域成像的速度建模方法及系统 |
| CN104536043A (zh) * | 2014-12-26 | 2015-04-22 | 中国石油天然气股份有限公司 | 一种深度域整体速度模型融合方法及装置 |
| CN106597533A (zh) * | 2016-11-17 | 2017-04-26 | 中国石油化工股份有限公司 | 一种用于山前带地震资料处理的深度域速度建模方法 |
Non-Patent Citations (1)
| Title |
|---|
| 杜宏钧: "逆时偏移技术在八岭山－花园地区的应用", 《江汉石油职工大学学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113960658A (zh) * | 2020-07-20 | 2022-01-21 | 中国石油天然气股份有限公司 | 一种基于地质地震模型的测井约束速度建模方法及装置 |
| CN113960658B (zh) * | 2020-07-20 | 2024-01-30 | 中国石油天然气股份有限公司 | 一种基于地质地震模型的测井约束速度建模方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111142159B (zh) | 2022-04-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| CN106094032B (zh) | 一种构建地层速度模型的方法 | |
| CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
| CN110361781B (zh) | 一种基于地震反射振幅分析的沉积期古地貌恢复方法 | |
| CN102879819B (zh) | 保持地震波场运动学特征的地震数据处理方法及装置 | |
| CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
| CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
| CN108663713B (zh) | 一种建立深度域构造模型的方法 | |
| CN102901985A (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
| WO2005026776A1 (en) | Wide-offset-range pre-stack depth migration method for seismic exploration | |
| CN107390266A (zh) | 基于角道集的速度更新方法及叠前深度偏移速度建模方法 | |
| CN111123359A (zh) | 随钻测井与地层格架约束的井周地震成像探测方法及装置 | |
| CN105629300B (zh) | 提高复杂构造偏移数据信噪比的方法 | |
| CN103149588A (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
| CN103852789B (zh) | 用于地震数据的非线性层析方法及其装置 | |
| CN106443791B (zh) | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 | |
| CN103076628A (zh) | 一种孔径优化的叠前时间偏移的处理方法 | |
| CN111142159B (zh) | 中深层异常体地震速度建模方法及系统 | |
| CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
| CN104614762B (zh) | 疏松砂岩气藏边界确定方法及装置 | |
| CN109932749B (zh) | 一种井震标定方法、装置 | |
| CN108693556B (zh) | 一种构造约束的速度平滑方法和系统 | |
| CN109975876A (zh) | 一种基于构造层位的井震融合速度模型的建模方法 | |
| CN113589375B (zh) | 基于倾斜层约束旅行时计算的vsp层速度反演方法 | |
| CN109839659B (zh) | 一种对叠前深度偏移剖面进行迭代优化的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |