CN110568490B - 一种高速层顶薄储层的识别方法 - Google Patents

一种高速层顶薄储层的识别方法 Download PDF

Info

Publication number
CN110568490B
CN110568490B CN201910774904.2A CN201910774904A CN110568490B CN 110568490 B CN110568490 B CN 110568490B CN 201910774904 A CN201910774904 A CN 201910774904A CN 110568490 B CN110568490 B CN 110568490B
Authority
CN
China
Prior art keywords
frequency
seismic
seismic data
time
thin
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
CN201910774904.2A
Other languages
English (en)
Other versions
CN110568490A (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
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
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, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201910774904.2A priority Critical patent/CN110568490B/zh
Publication of CN110568490A publication Critical patent/CN110568490A/zh
Application granted granted Critical
Publication of CN110568490B publication Critical patent/CN110568490B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/282—Application of seismic models, synthetic seismograms
    • 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles

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

一种高速层顶薄储层的识别方法
技术领域
本发明涉及一种高速层顶薄储层的识别方法，属于石油勘探技术领域。
背景技术
目前，识别高速层之上的薄储层(即高速层顶薄储层)主要有以下两种方法，方法一：该方法直接利用地震振幅反映高速层顶薄储层，但是由于薄储层形成的地震反射常常会淹没在高速层形成的强地震反射波形中，使得利用地震振幅不能够准确反映薄储层的分布；方法二：该方法基于地震资料中的高频成分与薄储层的响应密切相关，先采用分频信号处理方法(例如傅里叶变换、小波变换、S变换等)对地震资料进行分频处理，得到不同主频的地震数据体，然后从中选取较高频段的地震数据体进行地震解释以识别薄储层。然而，该方法通常选取的数据体频段较窄，这样会出现以下问题：1)使用较窄的高频段地震数据体进行地震解释，常常会由于频率太高而出现波形畸变，导致无法准确识别薄储层；2)由于高速层顶薄储层的岩性结构通常在横向上变化较快，较窄频率段的地震数据体并不能够准确反映这些横向变化。
综上所述，目前的薄储层识别方法，无法准确识别高速层之上的薄储层。
发明内容
本发明的目的是提供一种高速层顶薄储层的识别方法，用以解决目前的薄储层预测方法，无法准确识别高速层之上的薄储层的问题。
为实现上述目的，本发明提供了一种高速层顶薄储层的识别方法，该方法包括以下步骤：
对于勘探区内的任一钻井，根据测井资料和设定的一系列不同主频的地震子波，生成该钻井的一系列合成地震记录，分析该钻井的一系列合成地震记录，得到该钻井中高速层上的薄储层能够得到识别时的有效最低频率；
对于勘探区内的多口钻井，得到每口钻井的有效最低频率，结合每口钻井的有效最低频率，利用插值方法得到勘探区区域上的有效最低频率分布；
获取勘探区的地震资料，将所述地震资料进行时频分解，得到时频谱道集；
利用有效最低频率分布剔除时频谱道集中的低频部分，得到地震资料高频分量的时频谱；
对地震资料高频分量的时频谱进行时频分解的反变换，得到地震资料高频分量；
利用地震资料高频分量中的振幅值反映勘探区高速层顶薄储层的分布特征。
本发明的有益效果是：本发明用于反映勘探区高速层顶薄储层的分布特征的地震资料高频分量是通过勘探区区域上的有效最低频率分布得到的。第一：由于有效最低频率分布中的频率是高速层上的薄储层能够得到识别时的有效最低频率，那么利用有效最低频率分布剔除时，仅剔除了时频谱道集中的低于有效最低频率的低频部分，剩余的高频部分全部被保留了，就能保证得到的地震资料高频分量的频率段是较宽的，那么在利用频率段较宽的地震资料高频分量识别薄储层时，就不易出现波形畸变，从而能够准确识别高速层之上的薄储层；第二：由于有效最低频率分布中，不同位置对应的有效最低频率值各不相同，那么利用有效最低频率分布剔除时频谱道集中的低频部分时，就使得在不同位置所保留的高频段各不相同，进而保证了所得到的地震资料高频分量中包含多个高频段成分，这样在识别整个区域的薄储层时，就能够反映薄储层的横向变化特征。综上所述，本发明既能够准确识别高速层之上的薄储层，又能够反映薄储层的横向变化特征，由此得到的薄储层识别结果，能够为油气勘探提供依据。
进一步地，所述一系列不同主频的地震子波的起始主频值低于地震资料有效频率段的最低频率值。
进一步地，对于勘探区内的任一钻井，得到该钻井中高速层上的薄储层能够得到识别时的有效最低频率的步骤为：在该钻井的一系列地震子波的主频值依次增加的合成地震记录中，当高速层形成的强地震反射波形的上部刚开始出现明显的复波形式时，提取相应合成地震记录对应的地震子波的主频作为该钻井的有效最低频率。
利用高速层形成的强地震反射波形的上部刚出现明显的复波形式作为高速层上发育薄储层的标识，一方面利用波形变化作为识别标识，能够克服利用地震振幅不能够准确反映薄储层分布的缺陷；另一方面利用刚出现明显的复波形式时的合成地震记录对应的地震子波的主频作为有效最低频率，能够保证有效最低频率在地震资料的有效带宽内，进而实现利用地震资料识别高速层之上的薄储层。
进一步地，所述地震资料高频分量的时频谱利用分频公式得到，所述分频公式为：
Figure BDA0002174731860000031
其中，f、fe分别表示地震资料中某个线道对应的时频谱道集中的频率、有效最低频率分布中的有效最低频率，Ghigh(f)表示地震资料高频分量的时频谱。
利用有效最低频率分布中的有效最低频率，对地震资料中所有线道对应的时频谱道集中的频率进行一一筛选，剔除其中的低频部分，仅保留高频部分。一方面，由于有效最低频率分布中的频率是高速层上的薄储层能够得到识别时的有效最低频率，因此剔除时，仅剔除了时频谱道集中的低于有效最低频率的低频部分，剩余的高频部分全部被保留了，从而保证了所保留的高频部分的频率段较宽，即保证了地震资料高频分量的频率段是较宽的；另一方面，由于有效最低频率分布中，不同位置对应的有效最低频率值各不相同，那么在不同位置所保留的高频段也各不相同，从而保证了地震资料高频分量中包含多个高频段成分。
进一步地，所述地震子波的主频值从起始主频值起按照增加步长依次增加，所述增加步长为3Hz～5Hz。
进一步地，所述地震子波为雷克子波。
进一步地，所述插值方法为克里金插值方法。
附图说明
图1是本发明实施例中高速层顶薄储层的识别方法流程图；
图2-1是本发明实施例中钻井A的地震子波主频为20Hz～55Hz对应的合成地震记录图；
图2-2是本发明实施例中钻井A的地震子波主频为60Hz～95Hz对应的合成地震记录图；
图2-3是本发明实施例中钻井A的地震子波主频为100Hz～135Hz对应的合成地震记录图；
图3-1是本发明实施例中参考井的地质模型图；
图3-2是本发明实施例中参考井的正演地震剖面图；
图3-3是本发明实施例中参考井的正演地震剖面的高频分量图；
图4-1是本发明实施例中过钻井A、B和C的地震剖面图；
图4-2是本发明实施例中过钻井A、B和C的地震剖面的高频分量剖面图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及具体实施例对本发明进行进一步详细说明。
本实施例提供了一种高速层顶薄储层的识别方法(以下简称薄储层识别方法)，能够对勘探区的薄储层进行有效识别，并能揭示薄储层的横向变化特征。
本实施例中，勘探区的基底为石炭系凝灰岩，地震速度较高(为4500～5500m/s)。石炭系顶部发育一套薄储层，厚底较小，多小于10m，与石炭系高速层之间的泥岩夹层的厚度也较小。薄储层发育特征复杂，表现为横向连续性不强，变化较快。这些薄储层形成的地震反射，淹没在石炭系顶部强反射的伴生反射波谷内，共同形成强波谷地震反射，使得利用地震振幅不能够准确反映薄储层的分布。
如图1所示，利用本实施例的薄储层识别方法对勘探区高速层顶的薄储层进行识别，包括步骤1-步骤6：
步骤1、对于勘探区内的任一钻井，根据测井资料和设定的一系列不同主频的地震子波，生成该钻井的一系列合成地震记录，分析该钻井的一系列合成地震记录，得到该钻井中高速层上的薄储层能够得到识别时的有效最低频率。
本实施例中，以勘探区内三口钻井(钻井A、B和C)为例，对本实施例的薄储层识别方法进行介绍，钻井数目改变时该方法的实施过程与三口钻井类似，此处不再赘述。
其中，合成地震记录的制作过程为现有技术，本实施例中钻井A的一系列合成地震记录的制作过程包括：
1)利用钻井A的井旁地震道集资料，计算钻井A的井旁地震道主频fmain；
2)利用钻井A的测井资料(包括声波测井曲线和密度曲线)和与井旁地震道主频fmain一致的标准雷克子波，通过井震标定确定准确的时深关系；
3)从起始主频值开始，从低到高逐渐增加地震子波的主频值，得到一系列不同主频的地震子波；利用这一系列不同主频的地震子波，通过常规方法计算得到每个主频值对应的合成地震记录，从而得到钻井A的一系列合成地震记录。
本实施例中，地震资料主频较高，达到60Hz左右，且地震资料有效信号的频率段(即地震资料有效频率段)为30Hz～85Hz，由于地震资料有效频率段的最低频率值为30Hz，因而本实施例中，将地震子波的起始主频值设定为20Hz，使其稍低于最低频率值，并令地震子波主频值的增加步长为5Hz，所得到的钻井A的一系列合成地震记录如图2-1至图2-3所示，图中地震子波的最高主频值为135Hz，超出了地震资料有效频率段的最高频率值，是为了保证在高频部分有足够的范围，以确保能够准确确定薄储层的响应频率。作为其他实施方式，地震子波主频值的增加步长还可以从3Hz～5Hz中进行选择。
本实施例中的地震子波是雷克子波，作为其他实施方式，还可以根据勘探区的实际情况选择其他类型的地震子波。
本实施例中，得到钻井A的有效最低频率的过程如下：
如图2-1至图2-3所示，图中每个合成地震记录对应的地震子波的主频标注在图下边。从图中可以看出，最下部的强反射是高速层形成的。由于高速层上的薄储层一般不能够在地震有效频带范围内准确识别出顶底反射，因而通过波形变化来作为识别标志，认为在不同地震子波主频的合成地震记录上，对于高速层形成的强地震反射波形，当其上部出现明显的复波形式时，即可认为是能够识别该高速层之上的薄储层。
结合图2-1至图2-3，可以看出：随着地震子波主频不断增加，合成地震记录中的地震波形对薄储层的表征也越明显，只有当地震子波主频达到100Hz以上时，形成的地震波形才能够准确识别高速层顶部的薄储层。但这一频率已超过了地震资料的有效带宽(30Hz～85Hz)。由于地震资料频带宽度有限，为了利用地震资料识别高速层顶的薄储层，本实施例选择在合成地震记录中高速层形成的强地震反射波形的上部刚出现明显的复波形式时，作为识别薄储层的标志，如图2-2中的矩形框内所示，并将刚出现明显复波形式时的合成地震记录对应的地震子波的主频作为钻井A的有效最低频率fea。图2-2中刚出现明显复波形式时对应的地震子波的主频为80Hz，处于地震资料的有效带宽内。
其中，对不同的操作人员来讲，若他们对同一个合成地震记录中复波形式是否明显的定位不同，那他们所得到的有效最低频率的大小可能存在差异，但均需保证所得到的有效最低频率处于地震资料的有效带宽内。
步骤2、对于勘探区内的多口钻井，得到每口钻井的有效最低频率，结合每口钻井的有效最低频率，利用插值方法得到勘探区区域上的有效最低频率分布Farea。
利用步骤1可以得到勘探区内每口钻井的有效最低频率，本实施例中钻井A、B、C对应的有效最低频率用fea，feb，fec表示，结合钻井A、B和C的有效最低频率，利用克里金等常用插值方法，对勘探区各个位置的有效最低频率横向插值，得到勘探区区域上的有效最低频率分布Farea。
步骤3、获取勘探区的地震资料M0，将地震资料M0进行时频分解，得到时频谱道集G0。
步骤4、利用有效最低频率分布Farea剔除时频谱道集G0中的低频部分，得到地震资料高频分量的时频谱Ghigh。
其中，利用分频公式得到地震资料高频分量的时频谱Ghigh，分频公式表示为：
Figure BDA0002174731860000071
公式中，f表示地震资料中某个线道对应的时频谱道集中的频率，fe表示该线道对应的有效最低频率分布Farea中的有效最低频率，Ghigh(f)即地震资料高频分量的时频谱Ghigh。
步骤5、对地震资料高频分量的时频谱Ghigh进行时频分解的反变换，得到地震资料高频分量Mhigh。
步骤6、利用地震资料高频分量Mhigh中的振幅值反映待测钻井高速层顶薄储层的分布特征。
其中，本实施例的薄储层识别方法在实施过程中，各步骤之间的顺序可根据实际需要调整，例如得到有效最低频率分布Farea和得到时频谱道集G0之间的顺序可调整。
下面对本实施例的薄储层识别方法的有效性进行验证。
首先，通过构建地质模型，对地质模型进行地震正演得到地震剖面，并对地震剖面提取高频分量，分析其对高速层之上的薄储层的识别效果。其中，利用实际钻井作为参考井，建立接近实际地质特征的地质模型如图3-1所示，通过地震正演模拟得到的地震剖面图如图3-2所示，地震剖面的高频分量图如图3-3所示。
图3-1中参考井的曲线为自然电位曲线，能够反映砂岩储层、泥岩和凝灰岩等特征。建立的地质模型有三个薄储层：上部为薄层的砂岩储层发育在厚层泥岩之间，下部为两套超覆尖灭的薄砂岩储层；最下边的储层发育在凝灰岩顶部。
图3-2所示的地震剖面中，可见两套地震强反射和一套稍弱的地震反射：上部为薄储层形成的强反射；中部薄储层的反射由于调谐影响变弱；下部主要为凝灰岩基底形成的强反射，砂岩储层形成的反射淹没在其中。
结合图3-3，可以明显看到，利用地震剖面的高频分量能够准确确定超覆薄层的反射形态和尖灭点。
其次，在通过构建地质模型测试了利用地震资料高频分量识别高速层顶薄储层的识别效果后，将该方法应用在研究区，以得到能够满足生产需要的结果。
过钻井A、B和C的地震剖面及其高频分量剖面如图4-1和图4-2所示，图中下部的弱杂乱反射为石炭系凝灰岩地层的响应。在图4-1所示的地震剖面中，A、B、C三个钻井在石炭系高速层顶部均对应一套强波谷反射，图中已用虚线标出。这一套强波谷反射不能有效表征三口钻井中发育于石炭系高速层顶部的薄储层。而在图4-2所示的高频分量剖面中可见，石炭系顶部的强反射被削弱，保留的地震反射与三口钻井中石炭系顶部的薄储层能够对应。钻井A、B、C中的薄储层均不连通，且C井中发育两套薄储层，这些特征均在高频分量剖面中得到凸显，并且这些薄储层的横向尖灭变化也能够得到凸显。

Claims (6)

1.一种高速层顶薄储层的识别方法，其特征在于，该方法包括以下步骤：
对于勘探区内的任一钻井，根据测井资料和设定的一系列不同主频的地震子波，生成该钻井的一系列合成地震记录，分析该钻井的一系列合成地震记录，得到该钻井中高速层上的薄储层能够得到识别时的有效最低频率；对于勘探区内的任一钻井，得到该钻井中高速层上的薄储层能够得到识别时的有效最低频率的步骤为：在该钻井的一系列地震子波的主频值依次增加的合成地震记录中，当高速层形成的强地震反射波形的上部刚开始出现明显的复波形式时，提取相应合成地震记录对应的地震子波的主频作为该钻井的有效最低频率；
对于勘探区内的多口钻井，得到每口钻井的有效最低频率，结合每口钻井的有效最低频率，利用插值方法对勘探区各个位置的有效最低频率横向插值，得到勘探区区域上的有效最低频率分布；勘探区不同位置对应的有效最低频率值各不相同；
获取勘探区的地震资料，将所述地震资料进行时频分解，得到时频谱道集；
利用有效最低频率分布剔除时频谱道集中的低频部分，得到地震资料高频分量的时频谱；
对地震资料高频分量的时频谱进行时频分解的反变换，得到地震资料高频分量；
利用地震资料高频分量中的振幅值反映勘探区高速层顶薄储层的分布特征。
2.根据权利要求1所述的高速层顶薄储层的识别方法，其特征在于，所述一系列不同主频的地震子波的起始主频值低于地震资料有效频率段的最低频率值。
3.根据权利要求1所述的高速层顶薄储层的识别方法，其特征在于，所述地震资料高频分量的时频谱利用分频公式得到，所述分频公式为：
Figure FDA0003127232470000021
其中，f、fe分别表示地震资料中某个线道对应的时频谱道集中的频率、有效最低频率分布中的有效最低频率，Ghigh(f)表示地震资料高频分量的时频谱，G0为时频谱道集。
4.根据权利要求1所述的高速层顶薄储层的识别方法，其特征在于，所述地震子波的主频值从起始主频值起按照增加步长依次增加，所述增加步长为3Hz～5Hz。
5.根据权利要求1-4任一项所述的高速层顶薄储层的识别方法，其特征在于，所述地震子波为雷克子波。
6.根据权利要求1-4任一项所述的高速层顶薄储层的识别方法，其特征在于，所述插值方法为克里金插值方法。
CN201910774904.2A 2019-08-21 2019-08-21 一种高速层顶薄储层的识别方法 Active CN110568490B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910774904.2A CN110568490B (zh) 2019-08-21 2019-08-21 一种高速层顶薄储层的识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910774904.2A CN110568490B (zh) 2019-08-21 2019-08-21 一种高速层顶薄储层的识别方法

Publications (2)

Publication Number Publication Date
CN110568490A CN110568490A (zh) 2019-12-13
CN110568490B true CN110568490B (zh) 2021-09-17

Family

ID=68775789

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910774904.2A Active CN110568490B (zh) 2019-08-21 2019-08-21 一种高速层顶薄储层的识别方法

Country Status (1)

Country Link
CN (1) CN110568490B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111444154B (zh) * 2020-03-26 2020-11-27 涵涡智航科技(玉溪)有限公司 一种在单片机***中高效存储日志文件的方法
CN112180464B (zh) * 2020-11-03 2024-05-24 中国石油化工股份有限公司 一种储层物性的识别方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2683106A1 (en) * 2001-04-30 2002-11-07 The Regents Of The University Of California Frequency-dependent processing and interpretation (fdpi) of seismic data for identifying, imaging and monitoring fluid-saturated underground reservoirs
CN104991270A (zh) * 2015-06-08 2015-10-21 中国石油化工股份有限公司 一种利用地震正演模拟的叠置砂体识别方法
CN106707339B (zh) * 2016-11-22 2018-09-25 中国石油天然气集团有限公司 一种深度域相控反演方法
CN106405647B (zh) * 2016-12-22 2018-04-13 成都晶石石油科技有限公司 一种关于沉积地层厚度的调谐反演方法

Also Published As

Publication number Publication date
CN110568490A (zh) 2019-12-13

Similar Documents

Publication Publication Date Title
CN105334535B (zh) 一种识别薄储层隐蔽岩性油气藏的方法
CN106842307B (zh) 一种基于正演约束下波形分类再检索的储层精细预测方法
CN111596364B (zh) 基于高精度层序地层格架的地震沉积微相组合分析方法
CN101349764B (zh) 一种地震旋回分析方法
Yanhu et al. A method of seismic meme inversion and its application
CN108020863A (zh) 一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法
CN110568490B (zh) 一种高速层顶薄储层的识别方法
CN114994758B (zh) 碳酸盐岩断控储层的波阻抗提取与结构表征方法和***
CN110927817B (zh) 一种薄层叠置河道水平井轨迹优化设计及控制的方法
CN108254785B (zh) 储层确定方法和装置
CN113109875A (zh) 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法
CN107703548A (zh) 基于沉积物品质因子和回波损失级曲线峰谷的浅地层层界划分方法
CN111090117B (zh) 一种相控正演约束下的有效储层预测方法及***
CN112180464B (zh) 一种储层物性的识别方法
CN113589365B (zh) 基于时频域信息的储层尖灭线描述方法
CN115629417A (zh) 基于地震沉积学的多尺度融合及相控颗粒滩刻画方法
CN110673211B (zh) 一种基于测井与地震数据的品质因子建模方法
CN110632660B (zh) 基于地震数据体的薄砂体表征方法及装置
CN110568491B (zh) 一种品质因子q的估算方法
Jafari et al. Estimation of acoustic-impedance model by using model-based seismic inversion on the Ghar Member of Asmari Formation in an oil field in southwestern Iran
CN110389381B (zh) 一种基于地震属性的砂体储层预测方法及装置
CN110941028B (zh) 一种对碳酸盐岩溶蚀孔洞型地热能储层位置进行定位的方法和***
CN113806674A (zh) 古河道纵向尺度的量化方法、装置、电子设备及存储介质
Butorin Application of wavelet spectral decomposition for geological interpretation of seismic data
Sadeghi et al. Global stochastic seismic inversion using turning bands simulation and co-simulation

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