CN109387873A - 一种缝洞储集体反演方法及系统 - Google Patents
一种缝洞储集体反演方法及系统 Download PDFInfo
- Publication number
- CN109387873A CN109387873A CN201710662577.2A CN201710662577A CN109387873A CN 109387873 A CN109387873 A CN 109387873A CN 201710662577 A CN201710662577 A CN 201710662577A CN 109387873 A CN109387873 A CN 109387873A
- Authority
- CN
- China
- Prior art keywords
- fracture
- data
- post
- value
- continuity
- 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.)
- Pending
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
公开了一种缝洞储集体反演方法和系统。该方法包括:1)获取叠后地震数据;2)基于井资料以及从井资料获取的地震子波信息,对所述叠后地震数据进行反演,获得波阻抗数据体;3)基于缝洞储集体的非均质特征得到缝洞发育优势区域体的数值;将步骤2)得到的波阻抗体数据体和步骤3)得到的缝洞发育优势区域体的数值相乘,得到修正后的波阻抗数据体。本发明利用缝洞储层的非均质特征,屏蔽反演结果中连续性较强的低阻区域,优选出符合非均质空间规律的低阻抗体,用于地下溶洞的刻画,可以更加准确地对地下溶洞进行刻画,计算储集体规模,为储量评估以及后续油藏模拟提供基础。
Description
技术领域
本发明涉及地震勘探技术领域,更具体地,涉及一种缝洞储集体反演方法及系统。
背景技术
缝洞型储集体在塔河油田广泛分布,在地震剖面上呈现出串珠状强反射的特征。在该类储层预测中,平面上主要使用振幅变化率属性进行刻画,而在纵向上利用能量体属性进行刻画。然而因子波旁瓣的影响,地震反射波在纵向上具有拉伸效应。因此,直接在利用地震数据或者振幅类属性进行储集体空间计算会带来偏差,影响储集体的精细刻画及后续储量评估的可靠性。常用的解决方案是,利用井震标定估算出子波,再进行叠后地震反演去除子波旁瓣影响,并获得波阻抗信息。这样一方面可以消除子波旁瓣影响,使得纵向上更为聚焦;另一方面波阻抗具有明确的物理意义:例如缝洞储集体因为当中充填了油气,具有降低的阻抗;而围岩一般是灰岩,阻抗值相对较高;利用这种规律,我们可以使用阻抗结果进行实现缝洞储集体的精细刻画。但是在实际的生产过程中经常碰到一个问题:在地下溶洞的周边如果同时沉积一套砂泥质岩层,在阻抗上也会体现为一套连续的低阻体,其值域与洞穴对应的阻抗一致。此时,单一利用阻抗值大小对缝洞体刻画则会带来问题,尤其在利用阻抗值进行储集体规模计算的时候,计算机很可能会把这套连续的低阻层也归于缝洞体,这样就导致储量评估结果偏大。目前业内的反演方法有多种多样,但是基本上都是利用子波信息消除旁瓣影响并获得阻抗结果,尚无对应的技术去解决反演结果中连续低阻抗的问题。因此,有必要提出一种能够进行缝洞预测的方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
针对常规波阻抗反演在预测缝洞储集体中碰到的连续低阻抗问题,本发明提出了一种缝洞储集体反演方法和系统,利用溶洞储集体的空间特性对反演结果进行自动筛选,在储集体刻画中剔除连续性强的低阻地质体影响,使得利用阻抗进行缝洞预测更加准确,并可以更好地应用于储量评估工作。
根据本发明的一方面,提出一种缝洞储集体反演方法。该方法主要包括:
1)获取叠后地震数据;
2)基于井资料以及从井资料获取的地震子波信息,对所述叠后地震数据进行反演,获得波阻抗数据体;
3)基于缝洞储集体的非均质特征得到缝洞发育优势区域体的数值;
4)将步骤2)得到的波阻抗体数据体和步骤3)得到的缝洞发育优势区域体的数值相乘,得到修正后的波阻抗数据体。
优选地,所述步骤3)包括以下子步骤:
3.1)获取所述叠后地震数据的连续性属性;
3.2)对所述连续性属性进行平滑处理;
3.3)对进行平滑处理后的叠后地震数据设置阈值a,将连续性属性的值大于a的区域赋值为1,将连续性属性的值小于a的区域赋值为0。
优选地,在子步骤3.1)中,通过以下方式获取所述连续性属性:
根据所述叠后地震数据的空间梯度建立结构张量,提取相应的特征值以反映连续性属性。
优选地,在子步骤3.2)中,采用中值滤波对所述连续性属性进行平滑处理。
根据本发明的另一方面,提出一种缝洞储集体反演系统。该系统主要包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
1)获取叠后地震数据;
2)基于井资料以及从井资料获取的地震子波信息,对所述叠后地震数据进行反演,获得波阻抗数据体;
3)基于缝洞储集体的非均质特征得到缝洞发育优势区域体的数值;
4)将步骤2)得到的波阻抗体数据体和步骤3)得到的缝洞发育优势区域体的数值相乘,得到修正后的波阻抗数据体。
优选地,所述步骤3)包括以下子步骤:
3.1)获取所述叠后地震数据的连续性属性;
3.2)对所述连续性属性进行平滑处理;
3.3)对进行平滑处理后的叠后地震数据设置阈值a,将连续性属性的值大于a的区域赋值为1,将连续性属性的值小于a的区域赋值为0。
优选地,在子步骤3.1)中,通过以下方式获取所述连续性属性:
根据所述叠后地震数据的空间梯度建立结构张量,提取相应的特征值以反映连续性属性。
优选地,在子步骤3.2)中,采用中值滤波对所述连续性属性进行平滑处理。
本发明利用缝洞储层的非均质特征,屏蔽反演结果中连续性较强的低阻区域,优选出符合非均质空间规律的低阻抗体,用于地下溶洞的刻画,可以更加准确地对地下溶洞进行刻画,计算储集体规模,为储量评估以及后续油藏模拟提供基础。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的,或者将在并入本文中的附图和随后的具体实施例中进行详细陈述,这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1为根据本发明的示例性实施方案的缝洞储集体反演方法的流程图;
图2显示了示例性叠后地震数据;
图3显示了反演所得到的波阻抗数据体;
图4显示了从叠后地震数据计算得到的连续性属性值;
图5显示了对连续性属性进行平滑处理之后的结果;
图6显示了最终得到的波阻抗数据体。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
以下参照图1详细描述根据本发明示例性实施方案的缝洞储集体反演方法。
该方法主要包括:
步骤1:获取叠后地震数据。
原始的叠后地震数据是后续反演的数据基础。图2显示了示例性叠后地震数据,由图2可见,短轴的串珠状地震信号对应着地下发育的溶洞。
步骤2:基于井资料以及从井资料获取的地震子波信息,对所述叠后地震数据进行反演,获得波阻抗数据体。
反演方法有多种,可以采用常用的基于模型反演方法,利用井资料建立起地下初始模型,然后修正初始模型,使得在模型正演的数据与实际的叠后地震数据匹配。
从图3所示的反演结果中可以看出,串珠对应的区域具有低阻特征,同时纵向上会更加收敛;此外,图2中标识连续的地震信号,反演结果中对应着一套连续的低阻抗数据体;这种数据体可能是沉积过程中岩性变化导致,其阻抗值与溶洞对应的阻抗值接近,但非溶蚀导致。为了利用阻抗对溶洞进行精细刻画,需要将连续的低阻体去除。
步骤3:基于缝洞储集体的非均质特征得到缝洞发育优势区域体的数值。
在一个示例中,所述步骤3)包括以下子步骤:
3.1)获取所述叠后地震数据的连续性属性;
3.2)对所述连续性属性进行平滑处理;
3.3)对进行平滑处理后的叠后地震数据设置阈值a,将连续性属性的值大于a的区域赋值为1,将连续性属性的值小于a的区域赋值为0。
对叠后地震数据提取某种能够衡量空间连续性的属性,该类属性有多种计算方法,可以使用计算其结构张量,该属性根据地震数据的空间梯度建立结构张量,并提取相应的特征值,用于定量刻画地震数据的空间变化程度,如图4所示:属性值越大,代表空间连续性小,而图2中连续的地震反射轴,对应的张量属性值很小。
为了更好地刻画溶洞的优势区域,对非连续性属性进行平滑处理,例如可以采用中值滤波对所述连续性属性进行平滑处理。
如图5所示,可以看到在经过平滑处理后非连续性强的区域基本上成片状分布。
设定阈值a,对于图5中大于a的区域(地震空间变化较大)赋值为1,小于a的区域(地震空间连续性较好)赋值为0。其中阈值a需要根据实际情况进行判定,例如可以将阈值取值为50。
步骤4:将步骤2得到的波阻抗体数据体和步骤3得到的缝洞发育优势区域体的数值相乘,得到修正后的波阻抗数据体。
将图3所示的反演结果与所对应的属性相乘,滤掉空间连续性较好的区域,最终得到修正后的波阻抗结果,如图6所示。
本发明还提出了一种缝洞储集体的反演系统。该系统主要包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
1)获取叠后地震数据;
2)基于井资料以及从井资料获取的地震子波信息,对所述叠后地震数据进行反演,获得波阻抗数据体;
3)基于缝洞储集体的非均质特征得到缝洞发育优势区域体的数值;
4)将步骤2)得到的波阻抗体数据体和步骤3)得到的缝洞发育优势区域体的数值相乘,得到修正后的波阻抗数据体。
在一个示例中,所述步骤3)包括以下子步骤:
3.1)获取所述叠后地震数据的连续性属性;
3.2)对所述连续性属性进行平滑处理;
3.3)对进行平滑处理后的叠后地震数据设置阈值a,将连续性属性的值大于a的区域赋值为1,将连续性属性的值小于a的区域赋值为0。
在一个示例中,在子步骤3.1)中,通过以下方式获取所述连续性属性:
根据所述叠后地震数据的空间梯度建立结构张量,提取相应的特征值以反映连续性属性。
在一个示例中,在子步骤3.2)中,采用中值滤波对所述连续性属性进行平滑处理。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (8)
1.一种缝洞储集体反演方法,其特征在于,包括:
1)获取叠后地震数据;
2)基于井资料以及从井资料获取的地震子波信息,对所述叠后地震数据进行反演,获得波阻抗数据体;
3)基于缝洞储集体的非均质特征得到缝洞发育优势区域体的数值;
4)将步骤2)得到的波阻抗体数据体和步骤3)得到的缝洞发育优势区域体的数值相乘,得到修正后的波阻抗数据体。
2.根据权利要求1所述的缝洞储集体反演方法,其中,所述步骤3)包括以下子步骤:
3.1)获取所述叠后地震数据的连续性属性;
3.2)对所述连续性属性进行平滑处理;
3.3)对进行平滑处理后的叠后地震数据设置阈值a,将连续性属性的值大于a的区域赋值为1,将连续性属性的值小于a的区域赋值为0。
3.根据权利要求2所述的缝洞储集体反演方法,其中,在子步骤3.1)中,通过以下方式获取所述连续性属性:
根据所述叠后地震数据的空间梯度建立结构张量,提取相应的特征值以反映连续性属性。
4.根据权利要求2所述的缝洞储集体反演方法,其中,在子步骤3.2)中,采用中值滤波对所述连续性属性进行平滑处理。
5.一种缝洞储集体反演系统,其特征在于,所述系统包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
1)获取叠后地震数据;
2)基于井资料以及从井资料获取的地震子波信息,对所述叠后地震数据进行反演,获得波阻抗数据体;
3)基于缝洞储集体的非均质特征得到缝洞发育优势区域体的数值;
4)将步骤2)得到的波阻抗体数据体和步骤3)得到的缝洞发育优势区域体的数值相乘,得到修正后的波阻抗数据体。
6.根据权利要求5所述的缝洞储集体反演系统,其中,所述步骤3)包括以下子步骤:
3.1)获取所述叠后地震数据的连续性属性;
3.2)对所述连续性属性进行平滑处理;
3.3)对进行平滑处理后的叠后地震数据设置阈值a,将连续性属性的值大于a的区域赋值为1,将连续性属性的值小于a的区域赋值为0。
7.根据权利要求6所述的缝洞储集体反演系统,其中,在子步骤3.1)中,通过以下方式获取所述连续性属性:
根据所述叠后地震数据的空间梯度建立结构张量,提取相应的特征值以反映连续性属性。
8.根据权利要求6所述的缝洞储集体反演系统,其中,在子步骤3.2)中,采用中值滤波对所述连续性属性进行平滑处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710662577.2A CN109387873A (zh) | 2017-08-04 | 2017-08-04 | 一种缝洞储集体反演方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710662577.2A CN109387873A (zh) | 2017-08-04 | 2017-08-04 | 一种缝洞储集体反演方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109387873A true CN109387873A (zh) | 2019-02-26 |
Family
ID=65413802
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710662577.2A Pending CN109387873A (zh) | 2017-08-04 | 2017-08-04 | 一种缝洞储集体反演方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109387873A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112379435A (zh) * | 2020-10-30 | 2021-02-19 | 中国石油天然气集团有限公司 | 相控岩溶型缝洞集合体刻画方法及装置 |
CN113970787A (zh) * | 2020-07-22 | 2022-01-25 | 中国石油化工股份有限公司 | 物性参数反演方法、装置、计算机设备和存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103116186A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油天然气集团公司 | 一种小尺度非均质储集体的容积确定方法 |
CN103792576A (zh) * | 2014-01-28 | 2014-05-14 | 中国石油天然气股份有限公司 | 一种基于梯度结构张量的储层非均质检测方法及设备 |
CN103926617A (zh) * | 2014-04-21 | 2014-07-16 | 中国石油化工股份有限公司 | 一种缝洞储集体的检测方法和检测装置 |
CN104181598A (zh) * | 2014-09-04 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 计算地层不连续性属性值的方法及装置 |
CN105301647A (zh) * | 2014-06-10 | 2016-02-03 | 中国石油化工股份有限公司 | 区分灰质泥岩和砂岩的方法 |
CN105353407A (zh) * | 2015-10-28 | 2016-02-24 | 中国石油化工股份有限公司 | 一种叠后地震波阻抗反演方法 |
CN105589098A (zh) * | 2014-10-29 | 2016-05-18 | 中国石油天然气股份有限公司 | 碳酸盐岩去除沉积泥质横向影响的储层反演方法和系统 |
-
2017
- 2017-08-04 CN CN201710662577.2A patent/CN109387873A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103116186A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油天然气集团公司 | 一种小尺度非均质储集体的容积确定方法 |
CN103792576A (zh) * | 2014-01-28 | 2014-05-14 | 中国石油天然气股份有限公司 | 一种基于梯度结构张量的储层非均质检测方法及设备 |
CN103926617A (zh) * | 2014-04-21 | 2014-07-16 | 中国石油化工股份有限公司 | 一种缝洞储集体的检测方法和检测装置 |
CN105301647A (zh) * | 2014-06-10 | 2016-02-03 | 中国石油化工股份有限公司 | 区分灰质泥岩和砂岩的方法 |
CN104181598A (zh) * | 2014-09-04 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 计算地层不连续性属性值的方法及装置 |
CN105589098A (zh) * | 2014-10-29 | 2016-05-18 | 中国石油天然气股份有限公司 | 碳酸盐岩去除沉积泥质横向影响的储层反演方法和系统 |
CN105353407A (zh) * | 2015-10-28 | 2016-02-24 | 中国石油化工股份有限公司 | 一种叠后地震波阻抗反演方法 |
Non-Patent Citations (1)
Title |
---|
胡修全: "涪陵地区长兴组地震相分析及储层预测研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113970787A (zh) * | 2020-07-22 | 2022-01-25 | 中国石油化工股份有限公司 | 物性参数反演方法、装置、计算机设备和存储介质 |
CN112379435A (zh) * | 2020-10-30 | 2021-02-19 | 中国石油天然气集团有限公司 | 相控岩溶型缝洞集合体刻画方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2823710C (en) | Methods and systems regarding models of underground formations | |
CN107367757A (zh) | 碳酸盐岩断溶体的刻画方法 | |
US10921477B2 (en) | Method, apparatus, and system for eliminating frequency dispersion effect | |
CN103792576A (zh) | 一种基于梯度结构张量的储层非均质检测方法及设备 | |
CN105005080A (zh) | 一种利用振幅比率属性识别地层圈闭尖灭线的方法 | |
CN106291677A (zh) | 一种基于匹配追踪方法的叠后声波阻抗反演方法 | |
CN103941286A (zh) | 一种基于平面导航的快速三维层位解释方法 | |
CN108873065B (zh) | 砂岩优质储层预测方法及装置 | |
CN107861162B (zh) | 基于微电极测井资料的天然裂缝识别方法及系统 | |
CN103869358A (zh) | 一种基于直方图均衡的断层识别方法及设备 | |
CN109387873A (zh) | 一种缝洞储集体反演方法及系统 | |
CN108629459B (zh) | 储层含烃孔隙的检测方法及装置 | |
CN105301651A (zh) | 地震地质一体化低级序小断层确定方法及装置 | |
CN114462113A (zh) | 基于振幅比的隐蔽河道刻画方法、装置、电子设备及介质 | |
CN113627607A (zh) | 碳酸盐岩储层沉积相识别方法、装置、电子设备及介质 | |
CN110954947B (zh) | 时域线性低频融合方法及系统 | |
CN112649893B (zh) | 一种面向薄储层的多资料多参数融合建模方法及系统 | |
CN112147687B (zh) | 一种储层含气性预测方法及预测系统 | |
CN110308488B (zh) | 确定洞穴充填程度的方法及系统 | |
CN112394398A (zh) | 一种基于岩相与地震相的地层尖灭线识别方法 | |
CN107831537B (zh) | 高砂地比浊积岩储层描述方法 | |
CN105467447A (zh) | 相控趋势能量匹配的地震保幅评价方法 | |
CN116165709B (zh) | 一种基于缝洞体识别的油气勘探方法、装置及设备 | |
NO20190217A1 (en) | Correcting biases in microseismic-event data | |
CN113495293B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190226 |