CN107728198A - 地震道集成像方法及系统 - Google Patents

地震道集成像方法及系统 Download PDF

Info

Publication number
CN107728198A
CN107728198A CN201610665776.4A CN201610665776A CN107728198A CN 107728198 A CN107728198 A CN 107728198A CN 201610665776 A CN201610665776 A CN 201610665776A CN 107728198 A CN107728198 A CN 107728198A
Authority
CN
China
Prior art keywords
processing
common image
image gathers
stacked section
domain common
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
Application number
CN201610665776.4A
Other languages
English (en)
Other versions
CN107728198B (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.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
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 China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201610665776.4A priority Critical patent/CN107728198B/zh
Publication of CN107728198A publication Critical patent/CN107728198A/zh
Application granted granted Critical
Publication of CN107728198B publication Critical patent/CN107728198B/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

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)
  • Image Processing (AREA)

Abstract

公开了一种地震道集成像方法及系统。该方法可以包括：基于共中心点道集与地震速度场，通过叠前偏移，获得角度域共成像点道集与偏移距域共成像点道集；基于角度域共成像点道集与偏移距域共成像点道集，通过拉伸切除处理，获得处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面；基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面，通过系数加权与差值计算，获得差值剖面；以及基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面加权的组合和差值剖面的计算，获得最终的共成像点道集。

Description

地震道集成像方法及系统
技术领域
本发明涉及油汽地球物理技术领域，更具体地，涉及一种利用差值剖面的地震道集成像方法及系统。
背景技术
地震成像是地震勘探的关键技术，不论是在研究领域还是在工业应用中，都是反射地震方法有史以来的热点问题。一般来讲，叠前偏移除了可以输出成像结果，还可以输出未完全叠加的部分成像数据。把所有部分成像数据中成像点横向位置相同的道组合起来，就形成了共成像道集(CIG，Common Image Gather)。其有三方面的用途：其一，基于共成像道集的偏移速度分析；其二，在振幅保真程度较高的共成像道集上进行振幅随偏移距变化(AVO，Amplitude Versus Offset)或者振幅随角度变化(AVA，Amplitude Versus Angle)分析；其三，对共成像道集进行适当的切除、去噪处理、剩余曲率校正，然后再叠加成像，会进一步提高构造图像的质量。
同其他叠前偏移方法一样，共炮集叠前偏移也可输出共成像道集，得到偏移距域共成像点道集(ODCIG，Offset Domain Common Image Gather)，但在强横向变速情况下，偏移距域共成像点道集存在运动学和动力学上的假象。为了避开这些假象，在现有技术中，从角度域进行成像，得到角度域共成像点道集(ADCIG，Angle Domain Common ImageGather)。
发明人发现，无论是偏移距域共成像点道集还是角度域共成像点道集，两者在成像中都各自具备相对的优劣：偏移距域共成像点道集的浅层数据与深层比相对覆盖次数较少，浅层相对深层的畸变较大，在成像时会进行较大的切除；角度域共成像点道集则相反，深层数据相对浅层数据覆盖次数较多。因此，有必要开发一种高精度的地震道集成像方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种地震道集成像方法及系统，其通过多种地震道集联合成像，实现高精度的地震道集成像。
根据本发明的一方面，提出了一种地震道集成像方法，所述方法可以包括：基于共中心点道集与地震速度场，通过叠前偏移，获得角度域共成像点道集与偏移距域共成像点道集；基于所述角度域共成像点道集与所述偏移距域共成像点道集，通过拉伸切除处理，获得处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面；基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面，通过系数加权与差值计算，获得差值剖面；以及基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面加权的组合和所述差值剖面的计算，获得最终的共成像点道集。
根据本发明的另一方面，提出了一种地震道集成像系统，所述系统可以包括：用来基于共中心点道集与地震速度场，通过叠前偏移，获得角度域共成像点道集与偏移距域共成像点道集的单元；用来基于所述角度域共成像点道集与所述偏移距域共成像点道集，通过拉伸切除处理，获得处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面的单元；用来基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面，通过系数加权与差值计算，获得差值剖面的单元；以及用来基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面加权的组合和所述差值剖面的计算，获得最终的共成像点道集的单元。
根据本发明的地震道集成像方法，通过地震道集联合成像，引入差值剖面，突出角度域共成像点道集与偏移距域共成像点道集中各自的有效信号，从而使最终的共成像点道集的精度显著提高。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的地震道集成像方法的步骤的流程图。
图2示出了角度域共成像点道集与偏移距域共成像点道集成像的区别的示意图。
图3a与图3b分别示出了根据现有技术对偏移距域共成像点道集进行拉伸切除处理前和处理后的示意图。
图4a与图4b分别示出了根据现有技术的处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面的示意图。
图5示出了根据本发明的一个实施方式的最终的共成像点道集的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
实施方式1
图1示出了根据本发明的地震道集成像方法的步骤的流程图。
在该实施方式中，根据本发明的地震道集成像方法包括：步骤101，基于共中心点道集与地震速度场，通过叠前偏移，获得角度域共成像点道集与偏移距域共成像点道集；步骤102，基于角度域共成像点道集与偏移距域共成像点道集，通过拉伸切除处理，获得处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面；步骤103，基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面，通过系数加权与差值计算，获得差值剖面；以及步骤104，基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面加权的组合和差值剖面的计算，获得最终的共成像点道集。
该实施方式通过多种地震道集结合成像，实现高精度的地震道集成像。
下面详细说明根据本发明的地震道集成像方法的具体步骤。
在一个示例中，基于共中心点道集与地震速度场，通过叠前偏移，可以获得角度域共成像点道集与偏移距域共成像点道集。然后，基于角度域共成像点道集与偏移距域共成像点道集，通过拉伸切除处理，获得处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面。
在一个示例中，拉伸切除处理可以包括：消除远道畸变引入的噪音。
图2示出了角度域共成像点道集与偏移距域共成像点道集成像的区别的示意图。图3a和图3b分别示出了根据现有技术对偏移距域共成像点道集进行拉伸切除处理前和处理后的示意图，其中，图3a示出了根据现有技术对偏移距域共成像点道集进行拉伸切除处理前的示意图，图3b示出了根据现有技术拉伸切除处理后得到的处理后的偏移距域共成像点道集叠加剖面的示意图，其中，横轴表示共深度点序号，纵轴表示时间，单位为秒。
如图2所示，A为浅层大角度小偏移距处的能量，B为深层小角度大偏移距处的能量，在实际地震数据处理中，一般浅层因覆盖次数多，大角度的浅层能量对最终的成像效果贡献权重较小，且因为偏移距偏大的地方速度不是很精确而容易使得这部分能量没有被正确归位而形成噪音，即可以视为无效能量，如A能量所示；而深层覆盖次数少，能量发散较浅层严重，小角度的地震数据则在成像的效果中权重较高，是一种有效能量，如B能量所示。
如图3a与图3b所示，在共成像点道集的切除中，偏移距域共成像点道集对拉伸部分的切除将无效能量A保留，而将有效能量B切除，这对成像效果不利，而角度域共成像点道集的切除能克服此类缺陷，将A切除且保留B。可见，偏移距域共成像点道集和角度域共成像点道集两者在成像中都各自具备相对的优劣，而将偏移距域共成像点道集和角度域共成像点道集两者结合成像，能够显著提高最终的共成像点道集的精度。
在一个示例中，基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面，通过系数加权与差值计算，可以获得差值剖面。
在一个示例中，系数加权与差值计算可以为：
C1A(tr1)-C2O(tr2)＝Cha (1)
其中，A(tr1)为处理后的角度域共成像点道集叠加剖面，C1为处理后的角度域共成像点道集叠加剖面的加权系数，O(tr2)为处理后的偏移距域共成像点道集叠加剖面，C2为处理后的偏移距域共成像点道集叠加剖面的加权系数，Cha为差值剖面。
具体地，角度域共成像点道集截取tr1道后获得处理后的角度域共成像点道集叠加剖面A(tr1)，偏移距域共成像点道集截取tr2道后获得处理后的偏移距域共成像点道集叠加剖面O(tr2)，将A(tr1)与O(tr2)分别乘以加权系数C1和加权系数C2，达到能量归一化效果，使A(tr1)与O(tr2)的各点能量值属于同一能量级别，之后再相减得到差值剖面Cha。
在一个示例中，基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面加权的组合和差值剖面的计算，可以获得最终的共成像点道集。
在一个示例中，处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面加权的组合可以为：C4A(tr3)、C5O(tr4)、C6A(tr5)+C7O(tr6)。其中，A(tr3)与A(tr5)为处理后的角度域共成像点道集叠加剖面，C4与C6为处理后的角度域共成像点道集叠加剖面的加权系数，O(tr4)与O(tr6)为处理后的偏移距域共成像点道集叠加剖面，C5与C7为处理后的偏移距域共成像点道集叠加剖面的加权系数。
在一个示例中，基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面加权的组合和差值剖面的计算可以为：
C3Cha+C4A(tr3)＝F1 (2)
C3Cha+C5O(tr4)＝F2 (3)
C3Cha+C6A(tr5)+C7O(tr6)＝F3 (4)
其中，C3为差值剖面的加权系数，F1、F2、F3为最终的共成像点道集。
其中，A(tr1)、A(tr3)、A(tr5)可以满足：A(tr3)＝A(tr5)＝C1A(tr1)，O(tr2)、O(tr4)、O(tr6)可以满足：O(tr4)＝O(tr6)＝C2O(tr2)。
最终的共成像点道集F可以用三种方法得到，分别为：
C3Cha+C4A(tr3)＝F1 (2)
C3Cha+C5O(tr4)＝F2 (3)
C3Cha+C6A(tr5)+C7O(tr6)＝F3 (4)
其中，C3为差值剖面的加权系数，F1、F2、F3为最终的共成像点道集，也即，最终的共成像点道集F可以表示为F1、F2或F3。本领域的技术人员可以根据需要研究的目标能量来做参数调节和公式选择(公式(2)、(3)或(4))，优化计算，进而获得最终的共成像点道集F。
具体地，可以根据(2)将Cha与A(tr3)分别乘以加权系数C3与加权系数C4，达到能量归一化效果，使Cha与A(tr3)的各点能量值属于同一能量级别，之后再相加得到最终的共成像点道集F1。
可以根据(3)将Cha与O(tr4)分别乘以加权系数C3与加权系数C5，达到能量归一化效果，使Cha与O(tr4)的各点能量值属于同一能量级别，之后再相加得到最终的共成像点道集F2。
可以根据(4)将Cha、A(tr5)与O(tr6)分别乘以加权系数C3、加权系数C6与加权系数C7，达到能量归一化效果，使C3Cha的值域与C6A(tr5)+C7O(tr6)的值域在同一能量范围，之后再相加得到最终的共成像点道集F3。
由此可见，C1-C7都是为了使所在参数公式中的各求和项或者求差项的各项能量归一化，因此C1-C7也不一定是唯一值，需要本领域的技术人员在归一化的基础上根据需要目标能量来做调节，有做一定的灵活选择的空间。
根据本发明的地震道集成像方法，通过地震道集联合成像，引入差值剖面，突出角度域共成像点道集与偏移距域共成像点道集中各自的有效信号，从而使得最终的共成像点道集的精度显著提高。
应用示例
为便于理解本发明实施方式的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
图4a与图4b分别示出了根据现有技术的处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面的示意图，其中，图4a示出了根据现有技术的处理后的角度域共成像点道集叠加剖面的示意图，图4b示出了根据现有技术的处理后的偏移距域共成像点道集叠加剖面的示意图，其中，横轴表示共深度点序号，纵轴表示时间，单位为毫秒。对比发现：处理后的角度域共成像点道集叠加剖面比处理后的偏移距域共成像点道集叠加剖面深层成像清晰，处理后的偏移距域共成像点道集叠加剖面比处理后的角度域共成像点道集叠加剖面浅层成像清晰。
图5示出了根据本发明的一个实施方式的最终的共成像点道集的示意图，其中，横轴表示共深度点序号，纵轴表示时间，单位为毫秒。将处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面进行系数加权与差值计算，获得差值剖面，并将差值剖面加回角度域共成像点道集叠加剖面与偏移距域共成像点道集叠加剖面加权的组合，获得最终的共成像点道集。
由此可见，最终的共成像点道集结合角度域共成像点道集与偏移距域共成像点道集各自的优势，突出中各自的有效信号，不论浅、中、深层，成像效果均优于角度域共成像点道集与偏移距域共成像点道集，从而使精度显著提高。
本领域技术人员应理解，上面对本发明的实施方式的描述的目的仅为了示例性地说明本发明的实施方式的有益效果，并不意在将本发明的实施方式限制于所给出的任何示例。
实施方式2
根据本发明的实施方式，提供了一种地震道集成像系统，所述系统可以包括：用来基于共中心点道集与地震速度场，通过叠前偏移，获得角度域共成像点道集与偏移距域共成像点道集的单元；用来基于角度域共成像点道集与偏移距域共成像点道集，通过拉伸切除处理，获得处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面的单元；用来基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面，通过系数加权与差值计算，获得差值剖面的单元；以及用来基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面加权的组合和差值剖面的计算，获得最终的共成像点道集的单元。
该实施方式通过地震道集联合成像，引入差值剖面，突出角度域共成像点道集与偏移距域共成像点道集中各自的有效信号，从而使得最终的共成像点道集的精度显著提高。
在一个示例中，系数加权与差值计算可以为：
C1A(tr1)-C2O(tr2)＝Cha (1)
其中，A(tr1)为处理后的角度域共成像点道集叠加剖面，C1为处理后的角度域共成像点道集叠加剖面的加权系数，O(tr2)为处理后的偏移距域共成像点道集叠加剖面，C2为处理后的偏移距域共成像点道集叠加剖面的加权系数，Cha为差值剖面。
在一个示例中，处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面加权的组合可以为：C4A(tr3)、C5O(tr4)、C6A(tr5)+C7O(tr6)，其中，A(tr3)与A(tr5)为处理后的角度域共成像点道集叠加剖面，C4与C6为处理后的角度域共成像点道集叠加剖面的加权系数，O(tr4)与O(tr6)为处理后的偏移距域共成像点道集叠加剖面，C5与C7为处理后的偏移距域共成像点道集叠加剖面的加权系数。
在一个示例中，基于处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面加权的组合和差值剖面的计算可以为：
C3Cha+C4A(tr3)＝F1 (2)
C3Cha+C5O(tr4)＝F2 (3)
C3Cha+C6A(tr5)+C7O(tr6)＝F3 (4)
其中，C3为差值剖面的加权系数，F1、F2、F3为最终的共成像点道集。
在一个示例中，拉伸切除处理可以包括：消除远道畸变引入的噪音。
本领域技术人员应理解，上面对本发明的实施方式的描述的目的仅为了示例性地说明本发明的实施方式的有益效果，并不意在将本发明的实施方式限制于所给出的任何示例。
以上已经描述了本发明的各实施方式，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施方式。在不偏离所说明的各实施方式的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施方式的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施方式。

Claims (10)

1.一种地震道集成像方法，包括：
基于共中心点道集与地震速度场，通过叠前偏移，获得角度域共成像点道集与偏移距域共成像点道集；
基于所述角度域共成像点道集与所述偏移距域共成像点道集，通过拉伸切除处理，获得处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面；
基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面，通过系数加权与差值计算，获得差值剖面；以及
基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面加权的组合和所述差值剖面的计算，获得最终的共成像点道集。
2.根据权利要求1所述的地震道集成像方法，其中，所述系数加权与差值计算为：
C1A(tr1)-C2O(tr2)＝Cha (1)
其中，A(tr1)为所述处理后的角度域共成像点道集叠加剖面，C1为所述处理后的角度域共成像点道集叠加剖面的加权系数，O(tr2)为所述处理后的偏移距域共成像点道集叠加剖面，C2为所述处理后的偏移距域共成像点道集叠加剖面的加权系数，Cha为所述差值剖面。
3.根据权利要求2所述的地震道集成像方法，其中，所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面加权的组合为：C4A(tr3)、C5O(tr4)、C6A(tr5)+C7O(tr6)，
其中，A(tr3)与A(tr5)为所述处理后的角度域共成像点道集叠加剖面，C4与C6为所述处理后的角度域共成像点道集叠加剖面的加权系数，O(tr4)与O(tr6)为所述处理后的偏移距域共成像点道集叠加剖面，C5与C7为所述处理后的偏移距域共成像点道集叠加剖面的加权系数。
4.根据权利要求3所述的地震道集成像方法，其中，所述基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面加权的组合和所述差值剖面的计算为：
C3Cha+C4A(tr3)＝F1 (2)
C3Cha+C5O(tr4)＝F2 (3)
C3Cha+C6A(tr5)+C7O(tr6)＝F3 (4)
其中，C3为所述差值剖面的加权系数，F1、F2、F3为所述最终的共成像点道集。
5.根据权利要求1所述的地震道集成像方法，其中，所述拉伸切除处理包括：消除远道畸变引入的噪音。
6.一种地震道集成像系统，包括：
用来基于共中心点道集与地震速度场，通过叠前偏移，获得角度域共成像点道集与偏移距域共成像点道集的单元；
用来基于所述角度域共成像点道集与所述偏移距域共成像点道集，通过拉伸切除处理，获得处理后的角度域共成像点道集叠加剖面与处理后的偏移距域共成像点道集叠加剖面的单元；
用来基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面，通过系数加权与差值计算，获得差值剖面的单元；以及
用来基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面加权的组合和所述差值剖面的计算，获得最终的共成像点道集的单元。
7.根据权利要求6所述的地震道集成像系统，其中，所述系数加权与差值计算为：
C1A(tr1)-C2O(tr2)＝Cha (1)
其中，A(tr1)为所述处理后的角度域共成像点道集叠加剖面，C1为所述处理后的角度域共成像点道集叠加剖面的加权系数，O(tr2)为所述处理后的偏移距域共成像点道集叠加剖面，C2为所述处理后的偏移距域共成像点道集叠加剖面的加权系数，Cha为所述差值剖面。
8.根据权利要求7所述的地震道集成像系统，其中，所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面加权的组合为：C4A(tr3)、C5O(tr4)、C6A(tr5)+C7O(tr6)，
其中，A(tr3)与A(tr5)为所述处理后的角度域共成像点道集叠加剖面，C4与C6为所述处理后的角度域共成像点道集叠加剖面的加权系数，O(tr4)与O(tr6)为所述处理后的偏移距域共成像点道集叠加剖面，C5与C7为所述处理后的偏移距域共成像点道集叠加剖面的加权系数。
9.根据权利要求8所述的地震道集成像系统，其中，所述基于所述处理后的角度域共成像点道集叠加剖面与所述处理后的偏移距域共成像点道集叠加剖面加权的组合和所述差值剖面的计算为：
C3Cha+C4A(tr3)＝F1 (2)
C3Cha+C5O(tr4)＝F2 (3)
C3Cha+C6A(tr5)+C7O(tr6)＝F3 (4)
其中，C3为所述差值剖面的加权系数，F1、F2、F3为所述最终的共成像点道集。
10.根据权利要求6所述的地震道集成像系统，其中，所述拉伸切除处理包括：消除远道畸变引入的噪音。
CN201610665776.4A 2016-08-12 2016-08-12 地震道集成像方法及系统 Active CN107728198B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610665776.4A CN107728198B (zh) 2016-08-12 2016-08-12 地震道集成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610665776.4A CN107728198B (zh) 2016-08-12 2016-08-12 地震道集成像方法及系统

Publications (2)

Publication Number Publication Date
CN107728198A true CN107728198A (zh) 2018-02-23
CN107728198B CN107728198B (zh) 2019-02-19

Family

ID=61201534

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610665776.4A Active CN107728198B (zh) 2016-08-12 2016-08-12 地震道集成像方法及系统

Country Status (1)

Country Link
CN (1) CN107728198B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110297274A (zh) * 2018-03-21 2019-10-01 王高成 一种利用偏移公式校正共中心点道集实现叠加成像的方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102944894A (zh) * 2012-11-26 2013-02-27 中国科学院地质与地球物理研究所 一种地震叠前偏移成像方法
CN103777242A (zh) * 2012-10-24 2014-05-07 中国石油化工股份有限公司 一种深度聚焦和道集同相轴拉平联合的速度判别方法
US20140293744A1 (en) * 2013-04-02 2014-10-02 Bp Corporation North America Inc. Specular filter (sf) and dip oriented partial imaging (dopi) seismic migration
EP2673662B1 (en) * 2011-02-08 2014-12-17 Total SA Method of analyzing seismic data
CN104570077A (zh) * 2013-10-18 2015-04-29 中国石油化工股份有限公司 一种逆时偏移偏移距域共成像点道集提取方法
CN104635266A (zh) * 2015-01-14 2015-05-20 中国石油天然气集团公司 一种基于偏移距-方位角域的多维道集显示方法以及系统
US20150293248A1 (en) * 2014-04-09 2015-10-15 John Theodore Etgen Seismic Adaptive Focusing

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2673662B1 (en) * 2011-02-08 2014-12-17 Total SA Method of analyzing seismic data
CN103777242A (zh) * 2012-10-24 2014-05-07 中国石油化工股份有限公司 一种深度聚焦和道集同相轴拉平联合的速度判别方法
CN102944894A (zh) * 2012-11-26 2013-02-27 中国科学院地质与地球物理研究所 一种地震叠前偏移成像方法
US20140293744A1 (en) * 2013-04-02 2014-10-02 Bp Corporation North America Inc. Specular filter (sf) and dip oriented partial imaging (dopi) seismic migration
CN104570077A (zh) * 2013-10-18 2015-04-29 中国石油化工股份有限公司 一种逆时偏移偏移距域共成像点道集提取方法
US20150293248A1 (en) * 2014-04-09 2015-10-15 John Theodore Etgen Seismic Adaptive Focusing
CN104635266A (zh) * 2015-01-14 2015-05-20 中国石油天然气集团公司 一种基于偏移距-方位角域的多维道集显示方法以及系统

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110297274A (zh) * 2018-03-21 2019-10-01 王高成 一种利用偏移公式校正共中心点道集实现叠加成像的方法

Also Published As

Publication number Publication date
CN107728198B (zh) 2019-02-19

Similar Documents

Publication Publication Date Title
CN107505651B (zh) 地震初至波和反射波联合斜率层析成像方法
CN106405651B (zh) 一种基于测井匹配的全波形反演初始速度模型构建方法
CN105301641B (zh) 一种方位各向异性速度反演的方法及装置
Trad et al. Fast and robust deblending using apex shifted Radon transform
CN102778693A (zh) 一种基于反射波层拉平提取并消除的绕射波分离处理方法
CN106970417B (zh) 椭圆展开转换波速度分析方法与系统
WO2009105964A1 (zh) 一种三维地震记录的叠前似二维变换的方法
CN107817526A (zh) 叠前地震道集分段式振幅能量补偿方法及系统
CN104483705A (zh) 一种三维剩余静校正方法
CN110780351B (zh) 纵波和转换波叠前联合反演方法及系统
CN103744116B (zh) 一种叠前道集的全时间域相位一致性校正方法
CN103217710A (zh) 获取共偏移距共方位角域道集的方法及系统
CN105629300B (zh) 提高复杂构造偏移数据信噪比的方法
CN105510975B (zh) 提高地震数据信噪比的方法及装置
CN106249292A (zh) 一种共反射点道集优化处理方法及系统
CN108957553B (zh) 动校正量递推修正的无拉伸畸变动校正方法及装置
CN107728198A (zh) 地震道集成像方法及系统
CN106990434B (zh) 椭圆展开转换波成像方法及系统
CN109917461A (zh) 一种提高侏罗系地层及基岩成像质量的方法及装置
CN106814394B (zh) 台站式检波器与常规检波器联合低频恢复方法
CN105572741B (zh) 一种计算3d高频静校正量的方法
CN109143327B (zh) 一种基于块状约束的q值初始建模方法及系统
CN109425892B (zh) 地震子波的估计方法及系统
CN106932821A (zh) 地震层析反演中的一种目标射线追踪技术
CN112946742B (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