CN107797143B - 串珠状地震响应特征聚焦寻优方法及系统 - Google Patents

串珠状地震响应特征聚焦寻优方法及系统 Download PDF

Info

Publication number
CN107797143B
CN107797143B CN201610807449.8A CN201610807449A CN107797143B CN 107797143 B CN107797143 B CN 107797143B CN 201610807449 A CN201610807449 A CN 201610807449A CN 107797143 B CN107797143 B CN 107797143B
Authority
CN
China
Prior art keywords
replacement
velocity
optimal
beading
migration
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.)
Active
Application number
CN201610807449.8A
Other languages
English (en)
Other versions
CN107797143A (zh
Inventor
许凯
唐金良
王世星
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sinopec Geophysical Research Institute
China Petrochemical Corp
Original Assignee
Sinopec Geophysical Research Institute
China Petrochemical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sinopec Geophysical Research Institute, China Petrochemical Corp filed Critical Sinopec Geophysical Research Institute
Priority to CN201610807449.8A priority Critical patent/CN107797143B/zh
Priority to RU2017131249A priority patent/RU2760102C2/ru
Priority to US15/697,745 priority patent/US10288752B2/en
Publication of CN107797143A publication Critical patent/CN107797143A/zh
Application granted granted Critical
Publication of CN107797143B publication Critical patent/CN107797143B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/301—Analysis for determining seismic cross-sections or geostructures
    • 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

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

公开了一种串珠状地震响应特征聚焦寻优方法及系统。该方法可以包括：基于叠前时间偏移速度模型，进行速度替换，获得多个替换偏移速度体；基于所述多个替换偏移速度体，进行叠前时间偏移，获得多个叠后数据体；基于地震资料中串珠状地震响应特征，选取三维立方体窗口；基于所述多个叠后数据体，获得多个叠后数据体在所述三维立方体窗口中的统计瞬时振幅；以及基于所述统计瞬时振幅，获得最优偏移速度和最优地震道。

Description

串珠状地震响应特征聚焦寻优方法及系统
技术领域
本发明涉及地震勘探领域，更具体地，涉及一种串珠状地震响应特征聚焦寻优方法及系统。
背景技术
在地震勘探领域，针对碳酸盐岩非常规储层，常发育大型溶洞，在地震剖面上表现为串珠状地震响应特征，“串珠”成像与偏移速度模型密切相关，当偏移速度偏大时，“串珠”往上翘且能量较弱；当偏移速度偏小时，“串珠”则往下弯且能量较弱；当偏移速度正好时，则“串珠”最收敛且能量最强。而由于常规速度谱的拾取，网格较大，导致偏移速度模型不准确，因此“串珠”也常常出现上翘或下弯的现象，“串珠”成像不准确，不利于后续的钻井及进一步的储层预测。针对碳酸盐岩非常规储层“串珠”不收敛的问题，现有技术一般采用最小二乘偏移进行成像，以期达到“串珠”收敛。
发明人发现采用最小二乘偏移进行成像技术计算耗费较大，十分依赖初始模型且收敛极易陷入局部极值，导致串珠状地震响应特征始终不收敛；全波形反演亦面临同样的问题。因此，有必要开发一种能够提高储层预测精度的串珠状地震响应特征聚焦寻优方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种串珠状地震响应特征聚焦寻优方法及系统，其能够通过对目的层的偏移速度替换进行连续叠前时间偏移，得到不同成像效果的偏移剖面，进一步根据振幅能量差异，得到最优偏移速度和最优地震道，从而提高储层预测的精度和钻井的成功率。
根据本发明的一方面，提出了一种串珠状地震响应特征聚焦寻优方法。所述方法可以包括：基于叠前时间偏移速度模型，进行速度替换，获得多个替换偏移速度体；基于所述多个替换偏移速度体，进行叠前时间偏移，获得多个叠后数据体；基于地震资料中串珠状地震响应特征，选取三维立方体窗口；基于所述多个叠后数据体，获得多个叠后数据体在所述三维立方体窗口中的统计瞬时振幅；以及基于所述统计瞬时振幅，获得最优偏移速度和最优地震道。
根据本发明的另一方面，提出了一种串珠状地震响应特征聚焦寻优系统，所述系统可以包括：用于基于叠前时间偏移速度模型，进行速度替换，获得多个替换偏移速度体的单元；用于基于所述多个替换偏移速度体，进行叠前时间偏移，获得多个叠后数据体的单元；用于基于地震资料中串珠状地震响应特征，选取三维立方体窗口的单元；用于基于所述多个叠后数据体，获得多个叠后数据体在所述三维立方体窗口中的统计瞬时振幅的单元；以及用于基于所述统计瞬时振幅，获得最优偏移速度和最优地震道的单元。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的串珠状地震响应特征聚焦寻优方法的步骤的流程图。
图2示出了根据本发明的一个实施例的沿层速度替换示意图。
图3a-图3i分别示出了根据本发明的一个实施例的沿层速度替换叠前时间偏移剖面示意图。
图4示出了根据本发明的一个实施例的三维立方体窗口扫描示意图。
图5示出了根据本发明的一个实施例的某成像点振幅能量随偏移速度变化示意图。
图6示出了根据本发明的一个实施例的某测线寻优剖面示意图。
图7a和图7b分别示出了根据本发明的一个实施例的某测线寻优剖面放大示意图和某测线原始剖面放大示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
实施例1
图1示出了根据本发明的串珠状地震响应特征聚焦寻优方法的步骤的流程图。
在该实施例中，根据本发明的串珠状地震响应特征聚焦寻优方法可以包括：S1、基于叠前时间偏移速度模型，进行速度替换，获得多个替换偏移速度体；S2、基于所述多个替换偏移速度体，进行叠前时间偏移，获得多个叠后数据体；S3、基于地震资料中串珠状地震响应特征，选取三维立方体窗口；S4、基于所述多个叠后数据体，获得多个叠后数据体在所述三维立方体窗口中的统计瞬时振幅；以及S5、基于所述统计瞬时振幅，获得最优偏移速度和最优地震道。
该实施例通过对目的层的偏移速度替换进行连续叠前时间偏移，得到不同成像效果的偏移剖面，进一步根据振幅能量差异，得到最优偏移速度和最优地震道，从而实现提高储层预测的精度和钻井的成功率。
下面详细说明根据本发明的串珠状地震响应特征聚焦寻优方法的具体步骤。
在一个示例中，可以基于叠前时间偏移速度模型，进行速度替换，获得多个替换偏移速度体。
在一个示例中，获得所述多个替换偏移速度体可以包括：基于所述叠前时间偏移速度模型，设定最小替换速度、最大替换速度以及替换速度变化间隔，对已有偏移速度进行沿标志层速度替换，进而获得所述多个替换偏移速度体。具体地，可以首先对已有的叠前时间偏移速度模型进行沿标志层速度替换，设已有偏移速度为v，设定最小替换速度、最大替换速度以及合适的速度变化间隔，以得到不同替换偏移速度体。
在一个示例中，可以基于所述多个替换偏移速度体，进行叠前时间偏移，获得多个叠后数据体。可以分别采用工区已有的叠前时间偏移处理流程进行速度替换连续叠前时间偏移，从而获得不同的叠后数据体。
在一个示例中，可以基于地震资料中串珠状地震响应特征，选取三维立方体窗口；然后，可以基于所述多个叠后数据体，获得多个叠后数据体在所述三维立方体窗口中的统计瞬时振幅。
在一个示例中，获得所述统计瞬时振幅可以包括：根据L1范数准则，统计每个叠后数据体在所述三维立方体窗口中的瞬时振幅，进而获得多个叠后数据体的统计瞬时振幅。
具体地，根据实际资料中的“串珠”大小，可以选取合适的三维立方体窗口(每个样点属性由周围多道和多个采样点属性表示)，对地震资料进行滑动扫描，并根据L1范数准则，统计每个数据三维立方体窗口中的瞬时振幅属性大小，得到9个不同数据体的统计瞬时振幅大小。
其中，叠后地震道s(t)和相对应的瞬时振幅Amp(t)之间的关系满足：
s(t)＝Re(t)+i·Im(t) (1)
其中，s(t)表示叠后地震道，Amp(t)表示叠后地震道的瞬时振幅，Re(t)表示叠后地震道的实部，Im(t)表示叠后地震道的虚部。
在一个示例中，可以基于所计算的瞬时振幅，获得最优偏移速度和最优地震道。
在一个示例中，获得最优偏移速度和最优地震道可以包括：基于所述统计瞬时振幅，求取振幅极值对应的下标，进而获得所述最优偏移速度和所述最优地震道。优选地，可以通过最小二乘极值拟合算法对所述统计瞬时振幅进行拟合，获得所述振幅极值对应的下标(也即，最优下标)，然后，可以提取出相应的最优地震道和最优偏移速度。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
图2示出了根据本发明的一个实施例的沿层速度替换示意图，其中，根据新疆某工区实际地层的速度，设定最小替换速度为0.99v，最大替换替换速度为1.01v，速度替换变化间隔为0.0025v。
图3a-3i分别示出了沿层速度替换叠前时间偏移剖面示意图。其中，根据不同替换速度进行连续叠前时间偏移，得到不同的叠后数据体。图3e为初始偏移速度v，图3a-图3d中的替换速度分别为0.99v、0.9925v、0.995v、0.9975v；图3f-图3i中的替换速度分别为1.0025v、1.005v、1.0075v、1.01v。
图4示出了根据本发明的一个实施例的三维立方体窗口扫描示意图。其中，x轴为联络测线(Xline)，y轴为主测线(Inline)。根据地震资料中串珠状地震响应特征的大小，选取合适的三维立方体窗口，工区选取7×7×15的窗口，即中心点分别由纵测线和联络测线7道，纵向时间采样点15点表示；对同一样点不同替换速度立方体窗口内的振幅属性进行L1范数统计。
图5示出了根据本发明的一个实施例的某成像点振幅能量随偏移速度变化示意图。其中，以某一成像点为例，不同数据间均存在最优振幅值及其对应的下标，因此利用最小二乘拟合算法寻找振幅极值及其对应的下标；最终提取相应的最优地震道及最优偏移速度。
图6示出了根据本发明的一个实施例的某测线寻优剖面示意图。图7a和图7b分别示出了根据本发明的一个实施例的某测线寻优剖面放大示意图和某测线原始剖面放大示意图。其中，图7a的寻优剖面相较于图7b的原始剖面，双串珠更收敛，且分开度较好，单串珠“尾巴”较小，且更收敛。可以得出通过实施本发明的串珠状地震响应特征聚焦寻优方法得到的剖面，串珠成像更收敛，更精确，从而有效提高了储层预测的精度和钻井的成功率。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
实施例2
根据本发明的实施例，提供了一种串珠状地震响应特征聚焦寻优系统，所述系统可以包括：用于基于叠前时间偏移速度模型，进行速度替换，获得多个替换偏移速度体的单元；用于基于所述多个替换偏移速度体，进行叠前时间偏移，获得多个叠后数据体的单元；用于基于地震资料中串珠状地震响应特征，选取三维立方体窗口的单元；用于基于所述多个叠后数据体，获得多个叠后数据体在所述三维立方体窗口中的统计瞬时振幅的单元；以及用于基于所述统计瞬时振幅，获得最优偏移速度和最优地震道的单元。
该实施例通过对目的层的偏移速度替换进行连续叠前时间偏移，得到不同成像效果的偏移剖面，进一步根据振幅能量差异，得到最优偏移速度和最优偏移成像剖面，从而实现提高储层预测的精度和钻井的成功率。
在一个示例中，获得所述多个替换偏移速度体可以包括：基于所述叠前时间偏移速度模型，设定最小替换速度、最大替换速度以及替换速度变化间隔，对已有偏移速度进行沿标志层速度替换，进而获得所述多个替换偏移速度体。
在一个示例中，获得所述统计瞬时振幅可以包括：根据L1范数准则，统计每个叠后数据体在所述三维立方体窗口中的瞬时振幅，进而获得多个叠后数据体的统计瞬时振幅。
在一个示例中，获得最优偏移速度和最优地震道可以包括：基于所述统计瞬时振幅，求取振幅极值对应的下标，进而获得所述最优偏移速度和所述最优地震道。
在一个示例中，通过最小二乘极值拟合算法对所述统计瞬时振幅进行拟合可以获得所述振幅极值对应的下标。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (10)

1.一种串珠状地震响应特征聚焦寻优方法，包括：
基于叠前时间偏移速度模型，进行速度替换，获得多个替换偏移速度体；
基于所述多个替换偏移速度体，进行叠前时间偏移，获得多个叠后数据体；
基于地震资料中串珠状地震响应特征，选取三维立方体窗口；
基于所述多个叠后数据体，获得多个叠后数据体在所述三维立方体窗口中的统计瞬时振幅；以及
基于所述统计瞬时振幅，获得最优偏移速度和最优地震道。
2.根据权利要求1所述的串珠状地震响应特征聚焦寻优方法，其中，获得所述多个替换偏移速度体包括：
基于所述叠前时间偏移速度模型，设定最小替换速度、最大替换速度以及替换速度变化间隔，对已有偏移速度进行沿标志层速度替换，进而获得所述多个替换偏移速度体。
3.根据权利要求1所述的串珠状地震响应特征聚焦寻优方法，其中，获得所述统计瞬时振幅包括：
根据L1范数准则，统计每个叠后数据体在所述三维立方体窗口中的瞬时振幅，进而获得多个叠后数据体的统计瞬时振幅。
4.根据权利要求1所述的串珠状地震响应特征聚焦寻优方法，其中，获得所述最优偏移速度和所述最优地震道包括：
基于所述统计瞬时振幅，求取振幅极值对应的下标，进而获得所述最优偏移速度和所述最优地震道。
5.根据权利要求4所述的串珠状地震响应特征聚焦寻优方法，其中，通过最小二乘极值拟合算法对所述统计瞬时振幅进行拟合而获得所述振幅极值对应的下标。
6.一种串珠状地震响应特征聚焦寻优系统，包括：
用于基于叠前时间偏移速度模型，进行速度替换，获得多个替换偏移速度体的单元；
用于基于所述多个替换偏移速度体，进行叠前时间偏移，获得多个叠后数据体的单元；
用于基于地震资料中串珠状地震响应特征，选取三维立方体窗口的单元；
用于基于所述多个叠后数据体，获得多个叠后数据体在所述三维立方体窗口中的统计瞬时振幅的单元；以及
用于基于所述统计瞬时振幅，获得最优偏移速度和最优地震道的单元。
7.根据权利要求6所述的串珠状地震响应特征聚焦寻优系统，其中，获得所述多个替换偏移速度体包括：
基于所述叠前时间偏移速度模型，设定最小替换速度、最大替换速度以及替换速度变化间隔，对已有偏移速度进行沿标志层速度替换，进而获得所述多个替换偏移速度体。
8.根据权利要求6所述的串珠状地震响应特征聚焦寻优系统，其中，获得所述统计瞬时振幅包括：
根据L1范数准则，统计每个叠后数据体在所述三维立方体窗口中的瞬时振幅，进而获得多个叠后数据体的统计瞬时振幅。
9.根据权利要求6所述的串珠状地震响应特征聚焦寻优系统，其中，获得所述最优偏移速度和所述最优地震道包括：
基于所述统计瞬时振幅，求取振幅极值对应的下标，进而获得所述最优偏移速度和所述最优地震道。
10.根据权利要求9所述的串珠状地震响应特征聚焦寻优系统，其中，通过最小二乘极值拟合算法对所述统计瞬时振幅进行拟合而获得所述振幅极值对应的下标。
CN201610807449.8A 2016-09-07 2016-09-07 串珠状地震响应特征聚焦寻优方法及系统 Active CN107797143B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201610807449.8A CN107797143B (zh) 2016-09-07 2016-09-07 串珠状地震响应特征聚焦寻优方法及系统
RU2017131249A RU2760102C2 (ru) 2016-09-07 2017-09-05 Способ и система автоматического распознавания центра залежи в карстовой пещере
US15/697,745 US10288752B2 (en) 2016-09-07 2017-09-07 Automatic focus identification method and system for Karst cave reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610807449.8A CN107797143B (zh) 2016-09-07 2016-09-07 串珠状地震响应特征聚焦寻优方法及系统

Publications (2)

Publication Number Publication Date
CN107797143A CN107797143A (zh) 2018-03-13
CN107797143B true CN107797143B (zh) 2019-11-12

Family

ID=61530724

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610807449.8A Active CN107797143B (zh) 2016-09-07 2016-09-07 串珠状地震响应特征聚焦寻优方法及系统

Country Status (1)

Country Link
CN (1) CN107797143B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110967754B (zh) * 2018-09-29 2022-03-08 中国石油化工股份有限公司 一种基于偏移速度寻优的缝洞储层充填与流体识别方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080137480A1 (en) * 2006-12-07 2008-06-12 Woodside Energy Limited Method of Building a Subsurface Velocity Model
CN101839999B (zh) * 2009-03-20 2011-12-07 中国石油集团东方地球物理勘探有限责任公司 一种确定叠前时间偏移最佳速度剖面的方法
CN102590864A (zh) * 2011-12-31 2012-07-18 中国石油集团西北地质研究所 两步法层析反演近地表建模方法
CN104459782B (zh) * 2014-12-10 2017-02-22 中国石油天然气股份有限公司 使用薄层的沿层层析反演地震速度建模方法以及建模单元

Also Published As

Publication number Publication date
CN107797143A (zh) 2018-03-13

Similar Documents

Publication Publication Date Title
CN101839999B (zh) 一种确定叠前时间偏移最佳速度剖面的方法
CN102928873B (zh) 基于四维能量聚焦的地面微地震定位方法
CN102565857B (zh) 自动剩余动校正方法
CN102841373A (zh) 基于方位角约束的微地震事件定位方法
CN110133715A (zh) 一种基于初至时差和波形叠加的微地震震源定位方法
CN104483705A (zh) 一种三维剩余静校正方法
CN106970417B (zh) 椭圆展开转换波速度分析方法与系统
CN101551467A (zh) 一种基于边界检测的初至波自动拾取方法
CN107797140B (zh) 串珠状地震响应特征聚焦寻优方法及系统
CN107797143B (zh) 串珠状地震响应特征聚焦寻优方法及系统
CN104267432B (zh) 一种基于规则化的转换波共转换点道集高精度抽取方法
CN106199704A (zh) 一种三维三分量海底电缆地震资料速度建模方法
CN106324682A (zh) 一种应用于冻土层地区的表层结构调查方法
CN102338885A (zh) 三分量vsp资料初至时间自动拾取方法
CN106837313B (zh) Lwf存储式声波测井慢度提取方法
CN106761715B (zh) Lwf存储式声波测井首波检测方法
CN107870359B (zh) 微地震事件识别方法及装置
CN103777242A (zh) 一种深度聚焦和道集同相轴拉平联合的速度判别方法
CN106894762A (zh) 一种灌注桩冲击成孔的控制方法
CN104765066B (zh) 地震三维波速扫描聚焦成像方法
CN102053275A (zh) 一种用于单点地震室内组合的相对静校正量计算方法
CN104536048B (zh) 一种自动切除地震波波至的方法
CN106226815A (zh) 转换波等效速度比谱的获取方法及转换波速度分析方法
CN106443829A (zh) 一种近地表模型构建方法及装置
CN105223630A (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