CN113589380A - 地震资料断层储层信息空间融合方法 - Google Patents
地震资料断层储层信息空间融合方法 Download PDFInfo
- Publication number
- CN113589380A CN113589380A CN202010371016.9A CN202010371016A CN113589380A CN 113589380 A CN113589380 A CN 113589380A CN 202010371016 A CN202010371016 A CN 202010371016A CN 113589380 A CN113589380 A CN 113589380A
- Authority
- CN
- China
- Prior art keywords
- fault
- seismic data
- seismic
- data volume
- data
- 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
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 26
- 238000005070 sampling Methods 0.000 claims abstract description 23
- 230000002159 abnormal effect Effects 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000010276 construction Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000007499 fusion processing Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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
本发明提供一种地震资料断层储层信息空间融合方法,该地震资料断层储层信息空间融合方法包括:步骤1:选择地震数据体作为待融合数据体;步骤2:计算步骤1中地震数据体对应的断层属性数据体;步骤3:对指示断层位置数据采样点进行赋异常值,对非指示断层位置的采样点进行赋小值处理,从而得到赋值断层属性数据体;步骤4:将步骤1中的地震数据体和步骤3中赋值断层属性数据体进行代数相加,得到融合后的地震数据体。该地震资料断层储层信息空间融合方法将断层信息做增强处理后与原始地震数据体进行融合,增强了地震数据中的断层信息,在做地震属性分析时可以更清晰的显示断层附加信息,提高了地质分析的效率和准确度。
Description
技术领域
本发明涉及勘探地球物理地震数据处理方法领域,特别是涉及到一种地震资料断层储层信息空间融合方法。
背景技术
随着我国东部油田进入开发后期,剩余油多富集于“薄、小、深、散”的地质目标中,其中由低序级断层控制下的复杂断块油藏占到了相当大的比例。目前对低序级断层的刻画主要依赖于地震数据。地震数据中包含丰富的构造、储层及流体信息,其中构造信息包括地层的层位信息和断层信息。如何在保持原有地震数据储层、流体信息的情况下有效突出断层信息,对于准确刻画低序级断层是地震构造解释的关键。现有地震技术通常是对地震数据进行断层增强处理,同时参考不同的地震数据体进行构造解释,或是利用地震数据的断层属性从平面或是剖面上进行低序级断层的解释和组合。这种低序级断层的解释方式存在以下几个问题,1、需要选择不同的地震数据体进行对比分析,数据的选择标准不统一,构造解释的随意性比较大;2、原有方法中的无论是地震数据中的地层层位信息、储层信息和流体信息与地震断层属性中的断层指示信息是相互割裂,二者没有有机的融合起来,造成剖面上构造解释、平面上储层预测结果不高效。3、常规的复杂断块岩性构造油藏储层预测需要首先通过精细的构造解释准确描述断层形态,进一步通过沿层振幅属性结合地震振幅特征开展储层的预测,该方法储层预测周期长,不能适应复杂断块岩性油藏快速工区筛选、有利目标区快速定位等前期地层研究的需求。
在申请号:201510172789.3的中国专利申请中,涉及到一种改善时移地震资料一致性的融合处理方法,所述方法在基准地震数据和监测地震数据的资料处理过程中,在其中一个或多个一致性处理步骤中是将基准地震数据和监测地震数据融合在一起,构成一块三维数据,然后进行一致性处理;其中,所述一致性处理包括层析静校正、地表一致性振幅补偿、地表一致性反褶积、去噪、速度分析、地表一致性剩余静校正;所述方法包括:S1,数据输入:将时移地震资料的基准地震数据和监测地震数据分别进行输入;S2,预处理:根据时移地震资料的条件,对时移地震资料进行预处理;S3,融合数据:按统一的处理网格将经过预处理后的基准地震数据和监测地震数据融合成一套数据;S4,对融合后的数据进行一致性处理,得到处理好的记录;S5,数据输出:将处理好的记录输出。该专利是改善时移地震资料一致性的融合处理方法,无法加强或改善地震资料中断层以及储层的特征信息,因此无法帮助地质技术人员更有效地识别描述断层储层信息及分布特征。
为此我们发明了一种新的地震资料断层储层信息空间融合方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种将原始的地震数据与地震断层数据有机融合起来,提高了前期地质选区的效率和后期断层构造解释、储层分布预测的精度的地震资料断层储层信息空间融合方法。
本发明的目的可通过如下技术措施来实现:地震资料断层储层信息空间融合方法,该地震资料断层储层信息空间融合方法包括:步骤1:选择地震数据体作为待融合数据体;步骤2:计算步骤1中地震数据体对应的断层属性数据体;步骤3:对指示断层位置数据采样点进行赋异常值,对非指示断层位置的采样点进行赋小值处理,从而得到赋值断层属性数据体;步骤4:将步骤1中的地震数据体和步骤3中赋值断层属性数据体进行代数相加,得到融合后的地震数据体。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,选择现有地震数据体或者根据储层描述需求对现有地震数据体开展提高分辨率处理,作为待融合数据体。
在步骤2中,以步骤1中的待融合数据体或步骤1中的待融合数据体通过断层增强处理得到的数据体为基础,计算得到能够指示断层的属性体、方差体以及进一步优化得到蚂蚁体这些属性体。
在步骤2中,首先对步骤1中待融合数据体的进行断层增强处理,进一步根据断层增强处理结果计算方差体,最后利用方差体进行蚂蚁追踪得到蚂蚁体。
在步骤3中,针对步骤2中选择的断层属性数据体确定断层指示的门槛值,对指示断层位置数据采样点进行赋异常值,异常值须远大于步骤1中地震采样点的最大值或远小于地震采样点的最小值,对非指示断层位置的采样点进行赋小值处理,从而得到赋值断层属性数据体。
在步骤3中,异常值及小值的计算方法包含以下步骤:
步骤31,统计步骤1中地震数据的最大值Vmax及最小值Vmin;
步骤32,异常值选择远大于Vmax的数值或远小于Vmin;
小值是选择绝对值远小于Vmax和Vmin绝对值的数值,取零。
本发明中的地震资料断层储层信息空间融合方法,将断层信息做增强处理后与原始地震数据体进行融合,增强了地震数据中的断层信息,在做地震属性分析时可以更清晰的显示断层附加信息,提高了地质分析的效率和准确度。
附图说明
图1为本发明的一具体实施例中处理过程的流程图;
图2为本发明的一具体实施例中原始地震数据地震剖面的示意图;
图3为本发明的一具体实施例中断层强化处理地震剖面的示意图;
图4为本发明的一具体实施例中方差体地震剖面的示意图;
图5为本发明的一具体实施例中蚂蚁体地震剖面的示意图;
图6为本发明的一具体实施例中蚂蚁体地震赋值处理剖面的示意图;
图7为本发明的一具体实施例中融合数据体地震剖面的示意图;
图8为本发明的一具体实施例中原始地震数据沿层振幅属性的示意图;
图9为本发明的一具体实施例中融合数据体沿层振幅属性的示意图;
图10为本发明的地震资料断层储层信息空间融合方法的一具体实施例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图10所示,图10为本发明的地震资料断层储层信息空间融合方法的流程图。
步骤101:选择地震数据体作为待融合数据体,或根据储层描述需求对该地震数据体开展提高分辨率处理并作为新的待融合数据体;
步骤102:计算步骤101中地震数据体对应的断层属性数据体;步骤102中断层属性数据体可以通过步骤101中的地震数据或步骤101中的地震数据通过断层增强处理得到的数据体为基础,计算得到的能够指示断层的属性体、方差体以及进一步优化得到的蚂蚁体等属性体。
步骤103:针对步骤102中选择的断层属性数据体确定断层指示的门槛值,对指示断层位置数据采样点进行赋异常值(异常值须远大于步骤一中地震采样点的最大值或远小于地震采样点的最小值),对非指示断层位置的采样点进行赋小值处理,从而得到赋值断层属性数据体;
异常值及小值的计算方法包含以下步骤:
步骤1,统计步骤101中地震数据的最大值Vmax及最小值Vmi n;
步骤2,异常值通常选择为远大于Vmax的数值或远小于Vmin,可以对Vmax或Vmin分别乘以10倍作为异常值;
小值是指通常选择绝对值远小于Vmax和Vmin绝对值的数值,通常情况下可以取零。
步骤104:将步骤101中的地震数据体和步骤103中赋值断层属性数据体进行代数相加,得到融合后的地震数据体。
在应用本发明的一具体实施例中,如图1所示,包括了以下处理流程:
步骤一:选择胜利油田某A区块,其地震剖面如图2所示;
步骤二:首先对地震数据进行断层增强处理,其处理结果如图3所示;进一步对断层增强处理结果计算方差体,其结果如图4所示;最后利用方差体进行蚂蚁追踪得到蚂蚁体,其结果如图5所示。
步骤三:通过统计发现地震数据采样点的最大vmax=18000,最小值Vmin=-17900,本实施例中采样将最大值Vmax乘以10作为异常值;通过统计分析蚂蚁体在地震剖面断层位置处的数值,设定断层指示门槛为0.8,;利用门槛值和异常值对蚂蚁体进行赋值处理,赋值规则为当蚂蚁体采样点数值大于门槛值0.8时,采样点赋异常值,当采样点数值小于门槛值时赋零值处理,其处理结果如图6所示。
步骤四:将步骤一中的地震数据体和步骤三中赋值地震属性数据体进行代数相加,得到融合后的地震数据体,其结果如图7所示。
通过对比原始地震数据和地震断层体地震剖面可以看出,地震断层体剖面上可清晰准确的刻画低序级断层,有利于构造解释中低序级断层的解释和组合,提高了构造解释人员的解释效率。进一步对原始地震数据体和地震断层体取沿层地震振幅属性进行分析(如图8、图9所示),通过图9可以看出基于地震断层体的地震沿层振幅属性不但显示了地震振幅的相对强弱指示了砂体的发育位置,地震断层及低序级断层在平面上的展布也清晰的刻画出来,为构造解释人员开展储层、流体快速预测提供了可能。
Claims (6)
1.地震资料断层储层信息空间融合方法,其特征在于,该地震资料断层储层信息空间融合方法包括:
步骤1:选择地震数据体作为待融合数据体;
步骤2:计算步骤1中地震数据体对应的断层属性数据体;
步骤3:对指示断层位置数据采样点进行赋异常值,对非指示断层位置的采样点进行赋小值处理,从而得到赋值断层属性数据体;
步骤4:将步骤1中的地震数据体和步骤3中赋值断层属性数据体进行代数相加,得到融合后的地震数据体。
2.根据权利要求1所述的地震资料断层储层信息空间融合方法,其特征在于,在步骤1中,选择现有地震数据体或者根据储层描述需求对现有地震数据体开展提高分辨率处理,作为待融合数据体。
3.根据权利要求1所述的地震资料断层储层信息空间融合方法,其特征在于,在步骤2中,以步骤1中的待融合数据体或步骤1中的待融合数据体通过断层增强处理得到的数据体为基础,计算得到能够指示断层的属性体、方差体以及进一步优化得到蚂蚁体这些属性体。
4.根据权利要求1所述的地震资料断层储层信息空间融合方法,其特征在于,在步骤2中,首先对步骤1中待融合数据体的进行断层增强处理,进一步根据断层增强处理结果计算方差体,最后利用方差体进行蚂蚁追踪得到蚂蚁体。
5.根据权利要求1所述的地震资料断层储层信息空间融合方法,其特征在于,在步骤3中,针对步骤2中选择的断层属性数据体确定断层指示的门槛值,对指示断层位置数据采样点进行赋异常值,异常值须远大于步骤一中地震采样点的最大值或远小于地震采样点的最小值,对非指示断层位置的采样点进行赋小值处理,从而得到赋值断层属性数据体。
6.根据权利要求1所述的地震资料断层储层信息空间融合方法,其特征在于,在步骤3中,异常值及小值的计算方法包含以下步骤:
步骤31,统计步骤1中地震数据的最大值Vmax及最小值Vmin;
步骤32,异常值选择远大于Vmax的数值或远小于Vmin;
小值是选择绝对值远小于Vmax和Vmin绝对值的数值,取零。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010371016.9A CN113589380A (zh) | 2020-04-30 | 2020-04-30 | 地震资料断层储层信息空间融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010371016.9A CN113589380A (zh) | 2020-04-30 | 2020-04-30 | 地震资料断层储层信息空间融合方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113589380A true CN113589380A (zh) | 2021-11-02 |
Family
ID=78237960
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010371016.9A Pending CN113589380A (zh) | 2020-04-30 | 2020-04-30 | 地震资料断层储层信息空间融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113589380A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459801A (zh) * | 2014-12-10 | 2015-03-25 | 中国石油天然气集团公司 | 用于识别断层的相干增强处理方法 |
KR101656862B1 (ko) * | 2016-03-15 | 2016-09-13 | 한국지질자원연구원 | 추계학적 지진 단층 파열 모델링 장치 및 방법 |
WO2017020267A1 (zh) * | 2015-08-05 | 2017-02-09 | 深圳朝伟达科技有限公司 | 一种水力压裂微地震事件的解释方法 |
CN108931814A (zh) * | 2017-05-26 | 2018-12-04 | 中国石油化工股份有限公司 | 一种基于多属性融合的基岩裂缝预测的方法 |
US20190227185A1 (en) * | 2017-05-26 | 2019-07-25 | Chevron U.S.A. Inc. | System and method for predicting fault seal from seismic data |
CN110261905A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 基于倾角控制的复值相干微断层识别方法 |
CN110568493A (zh) * | 2019-08-21 | 2019-12-13 | 中国石油化工股份有限公司 | 一种复杂断块盆地隐蔽性断层的识别方法 |
-
2020
- 2020-04-30 CN CN202010371016.9A patent/CN113589380A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459801A (zh) * | 2014-12-10 | 2015-03-25 | 中国石油天然气集团公司 | 用于识别断层的相干增强处理方法 |
WO2017020267A1 (zh) * | 2015-08-05 | 2017-02-09 | 深圳朝伟达科技有限公司 | 一种水力压裂微地震事件的解释方法 |
KR101656862B1 (ko) * | 2016-03-15 | 2016-09-13 | 한국지질자원연구원 | 추계학적 지진 단층 파열 모델링 장치 및 방법 |
CN108931814A (zh) * | 2017-05-26 | 2018-12-04 | 中国石油化工股份有限公司 | 一种基于多属性融合的基岩裂缝预测的方法 |
US20190227185A1 (en) * | 2017-05-26 | 2019-07-25 | Chevron U.S.A. Inc. | System and method for predicting fault seal from seismic data |
CN110261905A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 基于倾角控制的复值相干微断层识别方法 |
CN110568493A (zh) * | 2019-08-21 | 2019-12-13 | 中国石油化工股份有限公司 | 一种复杂断块盆地隐蔽性断层的识别方法 |
Non-Patent Citations (5)
Title |
---|
卢美月;汤子余;张永健;: "蚂蚁追踪技术在Plutonio油田X区块中的应用", 科学技术与工程, vol. 20, no. 08, pages 2992 - 2995 * |
周赏;王永莉;韩天宝;李楠;史德莲;成鹏;方海飞;张颖;马意红;: "小断层综合解释技术及其应用", 石油地球物理勘探, vol. 47, no. 1, pages 50 - 54 * |
赵虎;尹成;朱仕军;: "多属性融合技术研究", 勘探地球物理进展, vol. 32, no. 02, pages 119 - 121 * |
郭彪;李勇;何剑;丁昭恒;: "基于小波变换的地震属性多源信息融合预测技术研究", 矿物岩石, vol. 35, no. 03, pages 115 - 120 * |
马春华;: "老地震资料再处理后井震结合断层识别方法", 油气田地面工程, vol. 32, no. 08, pages 3 - 4 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2600944C1 (ru) | Формирование моделей распознаваемых геологических структур на основании набора узловых точек | |
EP3217192B1 (en) | Method for looking for oil and gas reservoir on the basis of trap-3d software | |
US20160123121A1 (en) | Method of searching for oil-gas reservoir based on trap-3d software | |
US10732310B2 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
CN103777243A (zh) | 砂泥岩薄互层储层厚度预测方法 | |
US10890688B2 (en) | Method for generating secondary data in geostatistics using observed data | |
AU2013402201A1 (en) | In-situ wellbore, core and cuttings information system | |
CN105626057A (zh) | 一种确定陆相页岩气有机碳含量的测井优化方法 | |
Li et al. | Quantitative evaluation of the carbonate reservoir heterogeneity based on production dynamic data: A case study from Cretaceous Mishrif formation in Halfaya oilfield, Iraq | |
US10969507B2 (en) | Seismic channel geobody extraction for basin fairway simulation | |
CN106777707B (zh) | 一种利用改进的蜘蛛网图进行测井岩性定量识别的方法 | |
CN111239809B (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
CN106468171A (zh) | 多信息耦合的递进式砂砾岩岩相判别方法 | |
CN113589380A (zh) | 地震资料断层储层信息空间融合方法 | |
Huang | Integrating time-lapse seismic with production data: A tool for reservoir engineering | |
CN109444960A (zh) | 用于断层构造的判定方法 | |
CN113419274B (zh) | 基于高精度层序格架模型的三维地震切片属性体提取方法 | |
CN112649893B (zh) | 一种面向薄储层的多资料多参数融合建模方法及系统 | |
Masoud et al. | Reservoir Characterization and Geostatistical Model of the Cretaceous and Cambrian-Ordovician Reservoir Intervals, Meghil Field, Sirte Basin, Libya | |
CN105717541A (zh) | 储层流动单元划分方法 | |
CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
CN112347707A (zh) | 油藏流动非均质性的表征方法及装置 | |
CN108318936A (zh) | 一种地层划分处理方法和装置 | |
Ma et al. | Multi-level ultra-deep fault-controlled karst reservoirs characterization methods for the Shunbei field | |
CN115903026B (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 |