CN112649844B - 基于地震振幅差的反演方法及系统 - Google Patents

基于地震振幅差的反演方法及系统 Download PDF

Info

Publication number
CN112649844B
CN112649844B CN201910960266.3A CN201910960266A CN112649844B CN 112649844 B CN112649844 B CN 112649844B CN 201910960266 A CN201910960266 A CN 201910960266A CN 112649844 B CN112649844 B CN 112649844B
Authority
CN
China
Prior art keywords
seismic data
old
seis
new
inversion
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
CN201910960266.3A
Other languages
English (en)
Other versions
CN112649844A (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 Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical 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 Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201910960266.3A priority Critical patent/CN112649844B/zh
Publication of CN112649844A publication Critical patent/CN112649844A/zh
Application granted granted Critical
Publication of CN112649844B publication Critical patent/CN112649844B/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)
  • Geophysics And Detection Of Objects (AREA)

Abstract

公开了一种基于地震振幅差的反演方法及系统。该方法可以包括：针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据与归一化旧地震数据；计算归一化新地震数据与归一化旧地震数据的振幅差值；根据振幅差值建立关于反射系数的目标函数，计算目标函数趋于0时的反射系数；根据反射系数，针对每一道的振幅差值进行波阻抗反演，获得每一道的反演结果。本发明通过对地震振幅做差，然后对差异结果进行反演的方法，直接获得储层参数的变化信息，简化反演流程，通过对振幅差值的反演直接获得振幅差值结果的纵波阻抗反演结果。

Description

基于地震振幅差的反演方法及系统
技术领域
本发明涉及地球物理勘探领域，更具体地，涉及一种基于地震振幅差的反演方法及系统。
背景技术
地震反演技术是地震响应与储层参数之间最有效的桥梁，将地震反演技术应用到时移地震领域可以直接解决地震响应与油藏性质变化之间的关系。目前时移地震反演多采用对时移前后的地震数据分别进行反演，通过两次反演结果的差异得出储层参数的变化信息。而时移差异性反演可以直接对时移数据差异进行反演并获得储层参数的变化信息，与常规方法相比有着明显的计算效率优势，并可以简化反演流程。在该反演方法中，常用的手段则是对每次时移的地震数据进行反演，然后用反演的结果做差，对得到的差异结果进行分析。目前时移地震反演多采用对时移前后的地震数据分别进行反演，通过两次反演结果的差异得出储层参数的变化信息。现有的方法需要对每期的数据都做一次反演，计算量大。因此，有必要开发一种基于地震振幅差的反演方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种基于地震振幅差的反演方法及系统，其能够通过对地震振幅做差，然后对差异结果进行反演的方法，直接获得储层参数的变化信息，简化反演流程，通过对振幅差值的反演直接获得振幅差值结果的纵波阻抗反演结果。
根据本发明的一方面，提出了一种基于地震振幅差的反演方法。所述方法可以包括：针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据与归一化旧地震数据；计算所述归一化新地震数据与所述归一化旧地震数据的振幅差值；根据所述振幅差值建立关于反射系数的目标函数，计算所述目标函数趋于0时的反射系数；根据所述反射系数，针对每一道的振幅差值进行波阻抗反演，获得每一道的反演结果。
优选地，所述归一化新地震数据为：
其中，为第t个归一化新地震数据，seis(t)new为第t个新地震数据，seis(min)new为新地震数据的最小值，seis(max)new为新地震数据的最大值。
优选地，所述归一化旧地震数据为：
其中，为第t个归一化旧地震数据，seis(t)old为第t个旧地震数据，seis(min)old为旧地震数据的最小值，seis(max)old为旧地震数据的最大值。
优选地，所述目标函数为：
J＝(W(t)R(t)-seis(t)*)T(W(t)R(t)-seis(t)*)+λ(R(t)TR(t)-C) (3)
其中，J为目标函数，R(t)为反射系数，W(t)为给定子波，seis(t)*为振幅差值，λ为阻尼系数，C为常数。
优选地，通过公式(4)进行波阻抗反演：
其中，Zn+1为波阻抗反演，Z0为该道第一个振幅差值，N为每个地震道数值的总个数。
根据本发明的另一方面，提出了一种基于地震振幅差的反演系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据与归一化旧地震数据；计算所述归一化新地震数据与所述归一化旧地震数据的振幅差值；根据所述振幅差值建立关于反射系数的目标函数，计算所述目标函数趋于0时的反射系数；根据所述反射系数，针对每一道的振幅差值进行波阻抗反演，获得每一道的反演结果。
优选地，所述归一化新地震数据为：
其中，为第t个归一化新地震数据，seis(t)new为第t个新地震数据，seis(min)new为新地震数据的最小值，seis(max)new为新地震数据的最大值。
优选地，所述归一化旧地震数据为：
其中，为第t个归一化旧地震数据，seis(t)old为第t个旧地震数据，seis(min)old为旧地震数据的最小值，seis(max)old为旧地震数据的最大值。
优选地，所述目标函数为：
J＝(W(t)R(t)-seis(t)*)T(W(t)R(t)-seis(t)*)+λ(R(t)TR(t)-C) (3)
其中，J为目标函数，R(t)为反射系数，W(t)为给定子波，seis(t)*为振幅差值，λ为阻尼系数，C为常数。
优选地，通过公式(4)进行波阻抗反演：
其中，Zn+1为波阻抗反演，Z0为该道第一个振幅差值，N为每个地震道数值的总个数。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的基于地震振幅差的反演方法的步骤的流程图。
图2示出了根据本发明的一个实施例的归一化旧地震数据的示意图。
图3示出了根据本发明的一个实施例的归一化新地震数据的示意图。
图4示出了根据本发明的一个实施例的振幅差值的示意图。
图5示出了根据本发明的一个实施例的反演结果的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的基于地震振幅差的反演方法的步骤的流程图。
在该实施例中，根据本发明的基于地震振幅差的反演方法可以包括：步骤101，针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据与归一化旧地震数据；步骤102，计算归一化新地震数据与归一化旧地震数据的振幅差值；步骤103，根据振幅差值建立关于反射系数的目标函数，计算目标函数趋于0时的反射系数；以及步骤104，根据反射系数，针对每一道的振幅差值进行波阻抗反演，获得每一道的反演结果。
在一个示例中，归一化新地震数据为：
其中，为第t个归一化新地震数据，seis(t)new为第t个新地震数据，seis(min)new为新地震数据的最小值，seis(max)new为新地震数据的最大值。
在一个示例中，归一化旧地震数据为：
其中，为第t个归一化旧地震数据，seis(t)old为第t个旧地震数据，seis(min)old为旧地震数据的最小值，seis(max)old为旧地震数据的最大值。
在一个示例中，目标函数为：
J＝(W(t)R(t)-seis(t)*)T(W(t)R(t)-seis(t)*)+λ(R(t)TR(t)-C) (3)
其中，J为目标函数，R(t)为反射系数，W(t)为给定子波，seis(t)*为振幅差值，λ为阻尼系数，C为常数。
在一个示例中，通过公式(4)进行波阻抗反演：
其中，Zn+1为波阻抗反演，Z0为该道第一个振幅差值，N为每个地震道数值的总个数。
具体地，根据本发明的基于地震振幅差的反演方法可以包括：
针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据为公式(1)，归一化旧地震数据为公式(2)。
计算归一化新地震数据与归一化旧地震数据的振幅差值，根据振幅差值建立关于反射系数的目标函数为公式(3)，计算目标函数趋于0时的反射系数，根据反射系数，针对每一道的振幅差值，通过公式(4)进行波阻抗反演，获得每一道的反演结果。
本方法通过对地震振幅做差，然后对差异结果进行反演的方法，直接获得储层参数的变化信息，简化反演流程，通过对振幅差值的反演直接获得振幅差值结果的纵波阻抗反演结果。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
图2示出了根据本发明的一个实施例的归一化旧地震数据的示意图。
图3示出了根据本发明的一个实施例的归一化新地震数据的示意图。
针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据为公式(1)，如图2所示，归一化旧地震数据为公式(2)，如图3所示。
计算归一化新地震数据与归一化旧地震数据的振幅差值，如图4所示，根据振幅差值建立关于反射系数的目标函数为公式(3)，计算目标函数趋于0时的反射系数，根据反射系数，针对每一道的振幅差值，通过公式(4)进行波阻抗反演，获得每一道的反演结果，如图5所示。
综上所述，本发明通过对地震振幅做差，然后对差异结果进行反演的方法，直接获得储层参数的变化信息，简化反演流程，通过对振幅差值的反演直接获得振幅差值结果的纵波阻抗反演结果。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例，提供了一种基于地震振幅差的反演系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据与归一化旧地震数据；计算归一化新地震数据与归一化旧地震数据的振幅差值；根据振幅差值建立关于反射系数的目标函数，计算目标函数趋于0时的反射系数；根据反射系数，针对每一道的振幅差值进行波阻抗反演，获得每一道的反演结果。
在一个示例中，归一化新地震数据为：
其中，为第t个归一化新地震数据，seis(t)new为第t个新地震数据，seis(min)new为新地震数据的最小值，seis(max)new为新地震数据的最大值。
在一个示例中，归一化旧地震数据为：
其中，为第t个归一化旧地震数据，seis(t)old为第t个旧地震数据，seis(min)old为旧地震数据的最小值，seis(max)old为旧地震数据的最大值。
在一个示例中，目标函数为：
J＝(W(t)R(t)-seis(t)*)T(W(t)R(t)-seis(t)*)+λ(R(t)TR(t)-C) (3)
其中，J为目标函数，R(t)为反射系数，W(t)为给定子波，seis(t)*为振幅差值，λ为阻尼系数，C为常数。
在一个示例中，通过公式(4)进行波阻抗反演：
其中，Zn+1为波阻抗反演，Z0为该道第一个振幅差值，N为每个地震道数值的总个数。
本系统通过对地震振幅做差，然后对差异结果进行反演的方法，直接获得储层参数的变化信息，简化反演流程，通过对振幅差值的反演直接获得振幅差值结果的纵波阻抗反演结果。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (6)

1.一种基于地震振幅差的反演方法，其特征在于，包括：
针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据与归一化旧地震数据；
计算所述归一化新地震数据与所述归一化旧地震数据的振幅差值；
根据所述振幅差值建立关于反射系数的目标函数，计算所述目标函数趋于0时的反射系数；
根据所述反射系数，针对每一道的振幅差值进行波阻抗反演，获得每一道的反演结果；
其中，所述目标函数为：
J＝(W(t)R(t)-seis(t)*)T(W(t)R(t)-seis(t)*)+λ(R(t)TR(t)-C) (3)
其中，J为目标函数，R(t)为反射系数，W(t)为给定子波，seis(t)*为振幅差值，λ为阻尼系数，C为常数；
其中，通过公式(4)进行波阻抗反演：
其中，Zn+1为波阻抗反演，Z0为该道第一个振幅差值，N为每个地震道数值的总个数。
2.根据权利要求1所述的基于地震振幅差的反演方法，其中，所述归一化新地震数据为：
其中，为第t个归一化新地震数据，seis(t)new为第t个新地震数据，seis(min)new为新地震数据的最小值，seis(max)new为新地震数据的最大值。
3.根据权利要求1所述的基于地震振幅差的反演方法，其中，所述归一化旧地震数据为：
其中，为第t个归一化旧地震数据，seis(t)old为第t个旧地震数据，seis(min)old为旧地震数据的最小值，seis(max)old为旧地震数据的最大值。
4.一种基于地震振幅差的反演系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
针对时移地震前后两次的地震数据分别进行归一化处理，获得归一化新地震数据与归一化旧地震数据；
计算所述归一化新地震数据与所述归一化旧地震数据的振幅差值；
根据所述振幅差值建立关于反射系数的目标函数，计算所述目标函数趋于0时的反射系数；
根据所述反射系数，针对每一道的振幅差值进行波阻抗反演，获得每一道的反演结果；
其中，所述目标函数为：
J＝(W(t)R(t)-seis(t)*)T(W(t)R(t)-seis(t)*)+λ(R(t)TR(t)-C) (3)
其中，J为目标函数，R(t)为反射系数，W(t)为给定子波，seis(t)*为振幅差值，λ为阻尼系数，C为常数；
其中，通过公式(4)进行波阻抗反演：
其中，Zn+1为波阻抗反演，Z0为该道第一个振幅差值，N为每个地震道数值的总个数。
5.根据权利要求4所述的基于地震振幅差的反演系统，其中，所述归一化新地震数据为：
其中，为第t个归一化新地震数据，seis(t)new为第t个新地震数据，seis(min)new为新地震数据的最小值，seis(max)new为新地震数据的最大值。
6.根据权利要求4所述的基于地震振幅差的反演系统，其中，所述归一化旧地震数据为：
其中，为第t个归一化旧地震数据，seis(t)old为第t个旧地震数据，seis(min)old为旧地震数据的最小值，seis(max)old为旧地震数据的最大值。
CN201910960266.3A 2019-10-10 2019-10-10 基于地震振幅差的反演方法及系统 Active CN112649844B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910960266.3A CN112649844B (zh) 2019-10-10 2019-10-10 基于地震振幅差的反演方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910960266.3A CN112649844B (zh) 2019-10-10 2019-10-10 基于地震振幅差的反演方法及系统

Publications (2)

Publication Number Publication Date
CN112649844A CN112649844A (zh) 2021-04-13
CN112649844B true CN112649844B (zh) 2023-08-15

Family

ID=75343271

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910960266.3A Active CN112649844B (zh) 2019-10-10 2019-10-10 基于地震振幅差的反演方法及系统

Country Status (1)

Country Link
CN (1) CN112649844B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008140655A1 (en) * 2007-05-09 2008-11-20 Exxonmobil Upstream Research Company Inversion of 4d seismic data
CN102508294A (zh) * 2011-10-20 2012-06-20 西北大学 一种利用时移地震勘探资料进行差异avo分析的方法
WO2013159810A1 (en) * 2012-04-24 2013-10-31 Statoil Petroleum As Processing data representing a physical system
CN104062680A (zh) * 2013-03-22 2014-09-24 中国石油化工股份有限公司 一种计算波阻抗反演目标函数梯度的方法
CN106324665A (zh) * 2015-07-01 2017-01-11 中国石油化工股份有限公司 反演裂缝密度的方法及系统
CN106646603A (zh) * 2017-01-04 2017-05-10 中海石油(中国)有限公司 一种实际时移地震资料处理差异的可靠性判断方法
CN110068860A (zh) * 2019-05-09 2019-07-30 中国石油大学(北京) 一种时移地震数据处理方法、装置及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7636275B2 (en) * 2007-02-06 2009-12-22 Conocophillips Company Direct time lapse inversion of seismic data
FR2933499B1 (fr) * 2008-07-03 2010-08-20 Inst Francais Du Petrole Methode d'inversion conjointe de donnees sismiques representees sur des echelles de temps differentes
US9442204B2 (en) * 2012-08-06 2016-09-13 Exxonmobil Upstream Research Company Seismic inversion for formation properties and attenuation effects

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008140655A1 (en) * 2007-05-09 2008-11-20 Exxonmobil Upstream Research Company Inversion of 4d seismic data
CN102508294A (zh) * 2011-10-20 2012-06-20 西北大学 一种利用时移地震勘探资料进行差异avo分析的方法
WO2013159810A1 (en) * 2012-04-24 2013-10-31 Statoil Petroleum As Processing data representing a physical system
CN104062680A (zh) * 2013-03-22 2014-09-24 中国石油化工股份有限公司 一种计算波阻抗反演目标函数梯度的方法
CN106324665A (zh) * 2015-07-01 2017-01-11 中国石油化工股份有限公司 反演裂缝密度的方法及系统
CN106646603A (zh) * 2017-01-04 2017-05-10 中海石油(中国)有限公司 一种实际时移地震资料处理差异的可靠性判断方法
CN110068860A (zh) * 2019-05-09 2019-07-30 中国石油大学(北京) 一种时移地震数据处理方法、装置及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于稀疏优化的薄层反射系数反演及应用研究;张倩;《中国优秀硕士学位论文全文数据库基础科学辑》;20170215;第7-33页 *

Also Published As

Publication number Publication date
CN112649844A (zh) 2021-04-13

Similar Documents

Publication Publication Date Title
Dong et al. Denoising the optical fiber seismic data by using convolutional adversarial network based on loss balance
US10127905B2 (en) Apparatus and method for generating acoustic model for speech, and apparatus and method for speech recognition using acoustic model
US20200102821A1 (en) Reducing resonant noise in seismic data acquired using a distributed acoustic sensing system
CN109633753B (zh) 地表噪声压制方法及装置
Ma et al. Automated arrival-time picking using a pixel-level network
US10288754B2 (en) Seismic noise removal
CN106324665A (zh) 反演裂缝密度的方法及系统
CN112649844B (zh) 基于地震振幅差的反演方法及系统
CN113222029B (zh) 基于数据融合的地震成像方法、装置、设备及存储介质
CN109584895A (zh) 语音降噪方法及装置
Yuan et al. Attenuation of linear noise based on denoising convolutional neural network with asymmetric convolution blocks
CN104781698A (zh) 用于隔离地震数据中的信号的系统和方法
US20210011185A1 (en) Fiber optic noise mitigation for distributed acoustic sensing
CN111435174B (zh) 强反射地区地震资料振幅补偿方法及装置
CN113219539A (zh) 光纤声波传感地震数据中井筒波干扰的去除方法及装置
CN104749627A (zh) 基于相似性的微地震信号凸显方法
CN107346035B (zh) 识别断裂的方法及装置
CN112649857A (zh) 基于叠前地震数据的流体因子直接反演方法及系统
CN112114364B (zh) 偶极横波反射波补偿方法及装置
CN116610937B (zh) 在隐式空间上进行低频信息延拓的方法、装置及电子设备
CN117849875B (zh) 一种地震信号分析方法、系统、装置及存储介质
CN111435173B (zh) 岩性油藏叠后地震资料去噪方法及装置
CN110618448B (zh) 多波动态匹配方法及系统
CN113031071B (zh) 地震波长波长静校正方法及装置
US20230298326A1 (en) Image augmentation method, electronic device and readable storage medium

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