CN110794476A - 一种基于断溶体相控的反演方法 - Google Patents

一种基于断溶体相控的反演方法 Download PDF

Info

Publication number
CN110794476A
CN110794476A CN201810864177.4A CN201810864177A CN110794476A CN 110794476 A CN110794476 A CN 110794476A CN 201810864177 A CN201810864177 A CN 201810864177A CN 110794476 A CN110794476 A CN 110794476A
Authority
CN
China
Prior art keywords
wave impedance
inversion
solution
properties
frequency model
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
CN201810864177.4A
Other languages
English (en)
Other versions
CN110794476B (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 Northwest Oil Field Co
Original Assignee
China Petroleum and Chemical Corp
Sinopec Northwest Oil Field 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, Sinopec Northwest Oil Field Co filed Critical China Petroleum and Chemical Corp
Priority to CN201810864177.4A priority Critical patent/CN110794476B/zh
Publication of CN110794476A publication Critical patent/CN110794476A/zh
Application granted granted Critical
Publication of CN110794476B publication Critical patent/CN110794476B/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
    • G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本发明公开了一种基于断溶体相控的反演方法，所述方法包括：提取针对断溶体储集体的属性描述；基于所述属性描述，解析获取断溶体属性与波阻抗体间的转换关系；基于所述转换关系建立波阻抗反演低频模型；将所述波阻抗反演低频模型加入到常规反演流程中，得到断溶体地震相控制下的波阻抗体。根据本发明的方法，可以针对以溶蚀孔洞、裂缝及断裂带等组成的断溶体储集体进行波阻抗反演，从而实现对复杂碳酸盐岩储集体——“断溶体”储层的定量或半定量刻画，为波阻抗转换到孔隙度体做前期准备，进而完成有效储集体量化。

Description

一种基于断溶体相控的反演方法
技术领域
本发明涉及地质勘探领域，具体涉及一种基于断溶体相控的反演方法。
背景技术
在地质勘探领域，断控岩溶区是一种较为常见的地质构造区，其断控规模储集体为溶洞、断裂及裂缝综合体，一般称之为断溶体。
针对断控岩溶区，目前已经形成“振幅梯度属性+蚂蚁体+杂乱结构反射”多属性融合对断溶体立体结构进行刻画的技术，但这些属性刻画只能定性描述断溶体，要对断溶体进行定量或半定量刻画，目前常用方法是波阻抗反演技术。然而常规波阻抗反演只能反映溶洞型储层，且只能确定溶洞的顶面，对于溶洞内幕刻画及断裂和裂缝型储集体刻画稍显无力。
发明内容
本发明提供了一种基于断溶体相控的反演方法，所述方法包括：
提取针对断溶体储集体的属性描述；
基于所述属性描述，解析获取断溶体属性与波阻抗体间的转换关系；
基于所述转换关系建立波阻抗反演低频模型；
将所述波阻抗反演低频模型加入到常规反演流程中，得到断溶体地震相控制下的波阻抗体。
在一实施例中，解析获取断溶体属性与波阻抗体间的转换关系，包括：
将所述断溶体属性与已知井的阻抗进行交汇分析，获取属性体转换成波阻抗体的转换关系式。
在一实施例中，解析获取断溶体属性与波阻抗体间的转换关系，还包括：
根据钻井标定调整所述转换关系式。
在一实施例中，解析获取断溶体属性与波阻抗体间的转换关系，还包括：
如果没有钻井曲线，通过在常规阻抗反演结果上提取伪井资料或借鉴邻区相似井进行拟合。
在一实施例中，解析获取断溶体属性与波阻抗体间的转换关系，还包括：
拟合所述断溶体属性与孔隙度参数。
在一实施例中，提取针对断溶体储集体的属性描述，其中，提取针对断溶体储集体的结构张量属性。
在一实施例中：
在工区内实钻井的井点处沿纵向时间轴提取结构张量属性，作为该井的结构张量曲线；
利用波阻抗曲线与结构张量曲线之间的交汇图，拟合它们之间的转换关系式。
在一实施例中，基于所述转换关系建立波阻抗反演低频模型，其中，利用所述转换关系，结合碳酸盐岩储层的非均质性，编程获取作为低频模型的波阻抗体。
在一实施例中，将所述波阻抗反演低频模型加入到常规反演流程中，其中：
按照常规波阻抗反演方法进行井震标定及提取子波；
在选择低频模型趋势时，采用所述波阻抗反演低频模型。
本发明还提出了一种存储介质，所述存储介质上存储有可实现如本发明所述方法的程序代码。
根据本发明的方法，可以针对以溶蚀孔洞、裂缝及断裂带等组成的断溶体储集体进行波阻抗反演，从而实现对复杂碳酸盐岩储集体——“断溶体”储层的定量或半定量刻画，为波阻抗转换到孔隙度体做前期准备，进而完成有效储集体量化。
本发明的其它特征或优点将在随后的说明书中阐述。并且，本发明的部分特征或优点将通过说明书而变得显而易见，或者通过实施本发明而被了解。本发明的目的和部分优点可通过在说明书、权利要求书以及附图中所特别指出的步骤来实现或获得。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例共同用于解释本发明，并不构成对本发明的限制。在附图中：
图1是根据本发明一实施例的方法流程图；
图2是根据本发明一实施例测井波阻抗与断溶体属性之一的结构张量之间的拟合关系图；
图3是根据本发明一实施例由结构张量得到的低频模型剖面与对应的结构张量剖面的对比图；
图4是根据本发明一实施例断溶体相控波阻抗反演结果与常规反演结果对比图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式，借此本发明的实施人员可以充分理解本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程并依据上述实现过程具体实施本发明。需要说明的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。
在地质勘探领域，断控岩溶区是一种较为常见的地质构造区，其断控规模储集体为溶洞、断裂及裂缝综合体，一般称之为断溶体。
针对断控岩溶区，目前已经形成“振幅梯度属性+蚂蚁体+杂乱结构反射”多属性融合对断溶体立体结构进行刻画的技术，但这些属性刻画只能定性描述断溶体，要对断溶体进行定量或半定量刻画，目前常用方法是波阻抗反演技术。然而常规波阻抗反演只能反映溶洞型储层，且只能确定溶洞的顶面，对于溶洞内幕刻画及断裂和裂缝型储集体刻画稍显无力。
针对上述问题，本发明提出了一种基于断溶体相控的反演方法。在本发明的方法中，在常规波阻抗反演方法的基础上加入了可以作为低频模型的波阻抗体，从而可以针对以溶蚀孔洞、裂缝及断裂带等组成的断溶体储集体进行波阻抗反演。根据本发明的方法，可以实现对复杂碳酸盐岩储集体——“断溶体”储层的定量或半定量刻画，为波阻抗转换到孔隙度体做前期准备，进而完成有效储集体量化。
接下来基于附图详细描述根据本发明实施例的方法的详细流程，附图的流程图中示出的步骤可以在包含诸如一组计算机可执行指令的计算机系统中执行。虽然在流程图中示出了各步骤的逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
如图1所示，在一实施例中，方法包括：
提取针对断溶体储集体的属性描述(S110)；
基于提取到的属性描述，解析获取断溶体属性与波阻抗体间的转换关系(S120)；
基于获取的转换关系建立波阻抗反演低频模型(S130)；
将建立的波阻抗反演低频模型加入到常规反演流程中，得到断溶体地震相控制下的波阻抗体(S140)。
进一步的，在一实施例中，将波阻抗反演低频模型加入到常规反演流程中，其中：
按照常规波阻抗反演方法进行井震标定及提取子波；
在选择低频模型趋势时，采用波阻抗反演低频模型。
在上述步骤中，关键点之一在于获取断溶体属性与波阻抗体间的转换关系。具体的，在一实施例中，解析获取断溶体属性与波阻抗体间的转换关系，包括：
将断溶体属性与已知井的阻抗进行交汇分析，获取属性体转换成波阻抗体的转换关系式。
具体的，在一实施例中，提取能够描述断溶体储集体的多种属性，并将该属性与已知井的阻抗进行交汇分析，得到属性体转换成波阻抗体的关系式。
进一步的，在一实施例中，在建立属性体转换成波阻抗体的转换关系式的过程中，根据钻井标定调整转换关系式。
进一步的，在一实施例中，如果没有钻井曲线，通过在常规阻抗反演结果上提取伪井资料或借鉴邻区相似井进行拟合。
进一步的，在一实施例中，拟合过程还包括：拟合断溶体属性与孔隙度参数。
进一步的，在实际应用场景中，结构张量也叫梯度结构张量法，是目前雕刻断溶体轮廓最为有效的属性。因此，在一实施例中，在提取针对断溶体储集体的属性描述的过程中，提取针对断溶体储集体的结构张量属性。
具体的，在一实施例中：
在工区内实钻井的井点处沿纵向时间轴提取结构张量属性，作为该井的结构张量曲线；
利用波阻抗曲线与结构张量曲线之间的交汇图，拟合它们之间的转换关系式。
进一步的，在本发明的方法中，关键点还在于建立波阻抗反演低频模型。在一实施例中，利用断溶体属性与波阻抗体间的转换关系，结合碳酸盐岩储层的非均质性，编程获取作为低频模型的波阻抗体。
本发明的方法能够实现对复杂碳酸盐岩储集体——“断溶体”储层的定量或半定量刻画，为波阻抗转换到孔隙度体做前期准备，进而完成有效储集体量化。与常规反演方法对比，本发明的反演方法能够实现对断溶体的综合刻画与描述，并且能够实现对断溶体的定量雕刻，指导工区储量计算，并最终提高钻井成功率，提高勘探开发经济效益。
进一步的，在一实施例中，在反演得到的波阻抗数据体上，通过设定15000g/cm3*m/s这个门槛值，将小于该值的波阻抗提取出来，便能方便快捷的实现断溶体的空间展布雕刻。
基于本发明的方法，本发明还提出了一种存储介质，该存储介质上存储有可实现如本发明所述方法的程序代码。
接下来基于具体的应用实例详细描述根据本发明一实施例的方法的具体执行效果。
在本实施例中，应用的反演方法为约束稀疏脉冲反演。这里需要说明的是，在本发明其它实施例中，可以根据具体需要，采用其他的反演方法。
具体的，在本实施例中，反演方法目标函数为：
F＝Lp(r)+λLq(s-d)+a-1L1ΔZ (1)
目标函数等式右边分别是反射系数的绝对值和(Lp(r))，合成声波记录与原始地震数据的差值(λLq(s-d))，趋势约束项(a-1L1ΔZ)。约束稀疏脉冲反演基于地震数据，用波阻抗趋势进行约束，以期用最少的反射稀疏脉冲，对合成记录与地震道进行最佳匹配。波阻抗趋势即波阻抗反演低频模型。
在本实施例中，利用结构张量属性建立波阻抗反演低频模型。具体的，在工区内实钻井的井点处沿纵向时间轴提取结构张量属性，作为该井的结构张量曲线，然后利用波阻抗曲线与结构张量曲线之间的交汇图，拟合它们之间的关系式，如图2所示，图2为测井波阻抗与断溶体属性之一的结构张量之间的拟合关系图，灰色部分201是原始散点，黑线202是拟合的递归关系。纵轴y代表波阻抗，横轴x代表结构张量，拟合公式为：
y＝-507lnx+17550 (2)
在该关系式基础上，把结构张量x转换成波阻抗体y，即得到断溶体控制下的波阻抗低频模型。如图3所示，图3为由结构张量得到的低频模型剖面(图3上)与对应的结构张量剖面(图3下)的对比图。
需要注意的一点是，由于结构张量属性和波阻抗属性之间的关系是一个拟合过程，并不是一个确定关系式，在操作中需要反复测试，选取最优化的拟合关系式进行转换。
利用反演软件按常规反演步骤进行反演操作，提取子波，制作合成记录，在选择低频模型趋势时，选择由上一步得到的基于断溶体相控的低频模型，最终得到断溶体相控反演结果。如图4所示，图4为断溶体相控波阻抗反演结果(图4上)与常规反演结果(图4下)对比图。图中可以看出两种反演趋势都符合地震展布，但细节方面有所不同，结构张量约束的反演结果更能体现断溶体的模式。
虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。本发明所述的方法还可有其他多种实施例。在不背离本发明实质的情况下，熟悉本领域的技术人员当可根据本发明做出各种相应的改变或变形，但这些相应的改变或变形都应属于本发明的权利要求的保护范围。

Claims (10)

1.一种基于断溶体相控的反演方法，其特征在于，所述方法包括：
提取针对断溶体储集体的属性描述；
基于所述属性描述，解析获取断溶体属性与波阻抗体间的转换关系；
基于所述转换关系建立波阻抗反演低频模型；
将所述波阻抗反演低频模型加入到常规反演流程中，得到断溶体地震相控制下的波阻抗体。
2.根据权利要求1所述的方法，其特征在于，解析获取断溶体属性与波阻抗体间的转换关系，包括：
将所述断溶体属性与已知井的阻抗进行交汇分析，获取属性体转换成波阻抗体的转换关系式。
3.根据权利要求2所述的方法，其特征在于，解析获取断溶体属性与波阻抗体间的转换关系，还包括：
根据钻井标定调整所述转换关系式。
4.根据权利要求2所述的方法，其特征在于，解析获取断溶体属性与波阻抗体间的转换关系，还包括：
如果没有钻井曲线，通过在常规阻抗反演结果上提取伪井资料或借鉴邻区相似井进行拟合。
5.根据权利要求2所述的方法，其特征在于，解析获取断溶体属性与波阻抗体间的转换关系，还包括：
拟合所述断溶体属性与孔隙度参数。
6.根据权利要求2～5中任一项所述的方法，其特征在于，提取针对断溶体储集体的属性描述，其中，提取针对断溶体储集体的结构张量属性。
7.根据权利要求6所述的方法，其特征在于：
在工区内实钻井的井点处沿纵向时间轴提取结构张量属性，作为该井的结构张量曲线；
利用波阻抗曲线与结构张量曲线之间的交汇图，拟合它们之间的转换关系式。
8.根据权利要求2～7中任一项所述的方法，其特征在于，基于所述转换关系建立波阻抗反演低频模型，其中，利用所述转换关系，结合碳酸盐岩储层的非均质性，编程获取作为低频模型的波阻抗体。
9.根据权利要求1～8中任一项所述的方法，其特征在于，将所述波阻抗反演低频模型加入到常规反演流程中，其中：
按照常规波阻抗反演方法进行井震标定及提取子波；
在选择低频模型趋势时，采用所述波阻抗反演低频模型。
10.一种存储介质，其特征在于，所述存储介质上存储有可实现如权利要求1～9中任一项所述方法的程序代码。
CN201810864177.4A 2018-08-01 2018-08-01 一种基于断溶体相控的反演方法 Active CN110794476B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810864177.4A CN110794476B (zh) 2018-08-01 2018-08-01 一种基于断溶体相控的反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810864177.4A CN110794476B (zh) 2018-08-01 2018-08-01 一种基于断溶体相控的反演方法

Publications (2)

Publication Number Publication Date
CN110794476A true CN110794476A (zh) 2020-02-14
CN110794476B CN110794476B (zh) 2021-10-19

Family

ID=69425494

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810864177.4A Active CN110794476B (zh) 2018-08-01 2018-08-01 一种基于断溶体相控的反演方法

Country Status (1)

Country Link
CN (1) CN110794476B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112394416A (zh) * 2020-10-20 2021-02-23 中国石油天然气集团有限公司 非均质断裂控制储层预测方法及装置
CN112505760A (zh) * 2020-10-15 2021-03-16 中国石油天然气股份有限公司 一种基于虚拟井的断溶体储层刻画方法
CN115421214A (zh) * 2022-10-10 2022-12-02 成都理工大学 基于断裂识别和易溶岩相组合的断溶体刻画方法
CN116299701A (zh) * 2023-02-23 2023-06-23 北京阳光吉澳能源技术有限公司 断溶体储层识别方法、装置及电子设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101980053A (zh) * 2010-10-18 2011-02-23 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN103527184A (zh) * 2013-10-28 2014-01-22 北京大学 一种白云岩储层的预测方法和系统
US20150301208A1 (en) * 2014-04-22 2015-10-22 Westerngeco L.L.C. Seismic data processing
CN105223617A (zh) * 2015-09-21 2016-01-06 中国石油化工股份有限公司江汉油田分公司物探研究院 利用地震属性控制初始波阻抗的建模方法及系统
CN107390264A (zh) * 2017-07-20 2017-11-24 中国石油化工股份有限公司 碳酸盐岩断溶体内部结构的表征方法
US20180203146A1 (en) * 2015-07-28 2018-07-19 Schlumberger Technology Corporation Seismic Constrained Discrete Fracture Network

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101980053A (zh) * 2010-10-18 2011-02-23 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN103527184A (zh) * 2013-10-28 2014-01-22 北京大学 一种白云岩储层的预测方法和系统
US20150301208A1 (en) * 2014-04-22 2015-10-22 Westerngeco L.L.C. Seismic data processing
US20180203146A1 (en) * 2015-07-28 2018-07-19 Schlumberger Technology Corporation Seismic Constrained Discrete Fracture Network
CN105223617A (zh) * 2015-09-21 2016-01-06 中国石油化工股份有限公司江汉油田分公司物探研究院 利用地震属性控制初始波阻抗的建模方法及系统
CN107390264A (zh) * 2017-07-20 2017-11-24 中国石油化工股份有限公司 碳酸盐岩断溶体内部结构的表征方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HANQING LIU ET AL.: "Removing strong reflections of coal seams based on the fast exponential pursuit algorithm", 《SEG DENVER 2014 ANNUAL MEETING》 *
RADMILA PEDERSEN-TATALOVIC ET AL.: "Event-based low-frequency impedance modeling using well logs and seismic attributes", 《THE LEADING EDGE》 *
叶云飞等: "基于立体层析反演的低频模型构建在深水区储层反演中的应用:以南海深水 W 构造为例", 《吉林大学学报(地球科学版)》 *
鲍典等: "塔河油田碳酸盐岩断溶体油藏分隔性描述方法研究", 《新疆石油天然气》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112505760A (zh) * 2020-10-15 2021-03-16 中国石油天然气股份有限公司 一种基于虚拟井的断溶体储层刻画方法
CN112505760B (zh) * 2020-10-15 2024-05-28 中国石油天然气股份有限公司 一种基于虚拟井的断溶体储层刻画方法
CN112394416A (zh) * 2020-10-20 2021-02-23 中国石油天然气集团有限公司 非均质断裂控制储层预测方法及装置
CN115421214A (zh) * 2022-10-10 2022-12-02 成都理工大学 基于断裂识别和易溶岩相组合的断溶体刻画方法
CN115421214B (zh) * 2022-10-10 2023-10-20 成都理工大学 基于断裂识别和易溶岩相组合的断溶体刻画方法
CN116299701A (zh) * 2023-02-23 2023-06-23 北京阳光吉澳能源技术有限公司 断溶体储层识别方法、装置及电子设备
CN116299701B (zh) * 2023-02-23 2024-02-13 北京阳光吉澳能源技术有限公司 断溶体储层识别方法、装置及电子设备

Also Published As

Publication number Publication date
CN110794476B (zh) 2021-10-19

Similar Documents

Publication Publication Date Title
CN110794476B (zh) 一种基于断溶体相控的反演方法
RU2669948C2 (ru) Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности
Ray et al. Building more robust low-frequency models for seismic impedance inversion
CN108562938B (zh) 一种消除频散效应的方法、装置及系统
US10087733B2 (en) Fracture mapping using vertical seismic profiling wave data
CN103913774A (zh) 基于微地震事件的储层地质力学参数反演方法
CN110095811B (zh) 膏岩层速度模型构建处理方法及装置
CN108629459B (zh) 储层含烃孔隙的检测方法及装置
EP3014310A1 (en) System and method for remaining resource mapping
CN112379439A (zh) 地震数据中纵波与横波匹配的方法及装置
CN114428356A (zh) 相控反演刻画碳酸盐岩断溶体内部结构的方法及电子设备
CN117388924A (zh) 一种地质导向方法及装置
US9014984B2 (en) Resolution matched nonlinear resolution enhancement of well logs
CN112147700A (zh) 速度异常区的低频模型构建方法及系统
CN106291713B (zh) 一种确定地层岩性的处理方法及装置
CN112394416B (zh) 非均质断裂控制储层预测方法及装置
CN112415581B (zh) 缝洞储层反演方法及系统
CN112363244B (zh) 波阻抗反演方法与碳酸盐岩非均质储层预测方法及系统
US10180513B2 (en) Method for mapping the propagation of earth fractures
CN113009579B (zh) 地震数据反演方法及装置
CN118625386A (zh) 一种煤系地层的砂岩储层预测方法和装置及设备
CN116629140A (zh) 基于人工智能和实钻数据的地层四压力预测方法及装置
CN118091744A (zh) 一种碳酸盐岩非层状介质的低频模型构建方法及装置
CN115542399A (zh) 一种碳酸盐岩缝洞型储集体厚度分布预测方法
CN114135264A (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