CN105467447B - 相控趋势能量匹配的地震保幅评价方法 - Google Patents
相控趋势能量匹配的地震保幅评价方法 Download PDFInfo
- Publication number
- CN105467447B CN105467447B CN201410449069.2A CN201410449069A CN105467447B CN 105467447 B CN105467447 B CN 105467447B CN 201410449069 A CN201410449069 A CN 201410449069A CN 105467447 B CN105467447 B CN 105467447B
- Authority
- CN
- China
- Prior art keywords
- earthquake
- coefficient
- seismic
- imaging
- carried out
- 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
Links
Abstract
本发明提供一种利用相控趋势能量匹配来进行地震保幅性评价的方法。通过研究发现,地震(响应)属性与储层参数间通常存在一种软约束联系,即两者之间并非确定的数学关系,而是一定范围内的相关关系。本发明根据此原理,利用面向储层的正演模拟结果为(波形)保幅成像建模提供空间上连续的软约束信息,降低模型的不确定性,借以提高地震成像精度。本发明首次从沉积学、地震解释角度对地震资料保幅性进行评价。
Description
技术领域
本发明涉及勘探地球物理地震数据处理方法领域,具体是从沉积学、地震解释角度对地震资料保幅性进行评价的方法。
背景技术
随着勘探程度逐步加深,地震勘探的目标已由构造解释转变到岩性识别、油藏描述,面向岩性油藏刻画的地震数据保幅成像处理就显得异常重要。
如何对成像成果进行保幅性量化分析,一直以来争议较大。在生产实践中,研究人员也提出了一些保幅性评价方法,主要包括以下几个方面:①利用成像效果与时间切片分析成像的保幅性,即定性分析成像合理性及成像精度;②与已钻井的对比判定成像的保幅性,即利用合成记录分析储层的振幅变化关系与成像成果的振幅变化关系;③利用AVO等典型属性判别成像的保幅性,即分析处理后道集的AVO关系。这些方法都具有十分重要的现实意义,但这些方法都只能对地震资料保幅性进行定性评价,精度仍有待提高。
发明内容
本发明的目的是针对现有技术存在的问题,提供一种从沉积学、地震解释角度对地震资料保幅性评价,即相控趋势能量匹配的地震保幅评价方法,以提高储层预测的精度。
通过研究发现,地震(响应)属性与储层参数间通常存在一种软约束关系,即两者之间并非确定的数学关系,而是一定范围内的相关关系。本发明根据此原理,利用面向储层的正演模拟结果为(波形)保幅成像建模提供空间上连续的软约束信息,降低模型的不确定性,提高成像精度。纵向上地震成像结果与测井响应存在相关关系,符合压实趋势;横向上地震成像变化符合沉积规律,沉积相不同,地震反射特征各异。
本发明的目的可通过如下技术措施来实现:
步骤1,首先提出了相控趋势能量匹配的三个约束系数:
系数1,相关(相似)系数(R),指波形的相似系数。
公式为
系数2,幅度系数(S%),指两类数据等时采样点上数值的比值,其中A为采样点上的数值。公式为S%=A正演/A地震。
系数3,曲率系数(P%),指两类数据单位步长内拐点个数的比值,其中C为步长内拐点个数。公式为P%=C正演/C地震。
步骤2,根据三类约束系数对地震成像结果纵向进行定量约束。
利用单井测井地震响应与地震成像结果进行趋势能量匹配。将测井地震响应与地震成像结果进行相关、曲率及幅度系数运算。根据系数大小对地震成像结果纵向进行定量约束。
步骤3,根据三类约束系数对地震成像结果横向进行定量约束。
首先利用面向储层的二维正演模拟结果对地震成像结果横向进行定量约束。将二维地震响应与地震成像结果进行相关、曲率及幅度系数运算。根据系数大小对地震成像结果横向进行定量约束。
然后对面向储层的二维正演模拟结果及地震成像结果分别进行属性提取,利用正演提取的属性对成像结果属性横向进行定量约束。将正演提取的属性与地震成像结果属性进行相关、曲率及幅度系数运算。根据系数大小对地震成像结果横向进行定量约束。
附图说明
图1为相控趋势能量匹配的地震保幅评价方法实施流程图
图2为单井测井响应模拟结果
图3为地震成像与测井响应匹配图
图4为井旁道地震成像与测井地震响应相关分析图
图5为面向地质目标的相控地质模型
图6为地震正演模拟及属性提取结果
图7为地震成像及属性提取结果
图8为地震正演响应属性与地震成像属性叠合
图9为地震成像结果属性与地震正演模拟属性相关分析图
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文配合所附图式,作详细说明如下。
如图1所示,图1为本发明的相控趋势能量匹配的地震保幅评价方法的具体实施的流程图。
首先进行步骤101,对单井进行测井地震响应模拟(一维正演测井合成记录)(图2)。
步骤102,将测井地震响应与地震成像结果进行相关和幅度系数运算。根据系数大小对地震成像结果纵向进行定量约束(图3、图4)。幅度系数
70%<S%<130%并且相关系数R>60%时,说明地震成像数据纵向上达到了保幅成像的要求。
至此,相控趋势能量匹配的地震成像结果纵向评价完成。
步骤103,进行面向储层的相控地质模型建立(图5)。
步骤104,进行面向储层的二维正演模拟(图6)。
步骤105,对面向储层的二维正演模拟结果及地震成像结果分别进行属性提取(图6、图7)。
步骤106,将正演提取的属性与地震成像结果提取的属性进行相关、曲率及幅度系数运算。根据系数大小对地震成像结果横向进行定量约束。相关系数R>60%、幅度系数80%<S%<120%并且曲率系数70%<P%<130%时,说明地震成像数据在横向上达到了保幅成像的要求(图8、图9)。
至此,相控趋势能量匹配的地震保幅性评价完成。
Claims (1)
1.相控趋势能量匹配的地震保幅评价方法,其特征在于,包括步骤1:提出相控趋势能量匹配的三个约束系数;
该相控趋势能量匹配主要有三个约束系数:
系数1,相关系数或相似系数R,指波形的相似系数;
公式为;
系数2,幅度系数S%,指两类数据等时采样点上数值的比值,其中A为采样点上的数值;公式为S%=A正演/A地震;
系数3,曲率系数P%,指两类数据单位步长内拐点个数的比值,其中C为步长内拐点个数;公式为P%=C正演/C地震;
步骤2:根据三类约束系数对地震成像结果纵向进行定量约束;
利用单井测井地震响应与地震成像结果进行趋势能量匹配;将测井地震响应与地震成像结果进行相关、曲率及幅度系数运算;根据系数大小对地震成像结果纵向进行定量约束;
步骤3:根据三类约束系数对地震成像结果横向进行定量约束;
利用面向储层的二维正演模拟结果对地震成像结果纵向进行定量约束;将二维地震响应与地震成像结果进行相关、曲率及幅度系数运算;根据系数大小对地震成像结果横向进行定量约束;
对面向储层的二维正演模拟结果及地震成像结果分别进行属性提取,利用正演提取的属性对成像结果属性横向进行定量约束;将正演提取的地震属性与地震成像结果属性进行相关、曲率及幅度系数运算;根据系数大小对地震成像结果横向进行定量约束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410449069.2A CN105467447B (zh) | 2014-09-04 | 2014-09-04 | 相控趋势能量匹配的地震保幅评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410449069.2A CN105467447B (zh) | 2014-09-04 | 2014-09-04 | 相控趋势能量匹配的地震保幅评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105467447A CN105467447A (zh) | 2016-04-06 |
CN105467447B true CN105467447B (zh) | 2018-11-20 |
Family
ID=55605344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410449069.2A Active CN105467447B (zh) | 2014-09-04 | 2014-09-04 | 相控趋势能量匹配的地震保幅评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105467447B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108508489B (zh) * | 2018-03-19 | 2019-12-31 | 中国海洋石油集团有限公司 | 一种基于波形微变化匹配的地震反演方法 |
CN111090117B (zh) * | 2018-10-24 | 2021-11-19 | 中国石油化工股份有限公司 | 一种相控正演约束下的有效储层预测方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5719822A (en) * | 1996-10-04 | 1998-02-17 | Vector Seismic Data Processing, Inc. | Seismic data radon dip moveout method |
CN102176055B (zh) * | 2011-02-18 | 2012-10-17 | 中国石油化工股份有限公司 | 保幅处理分析评价方法 |
-
2014
- 2014-09-04 CN CN201410449069.2A patent/CN105467447B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5719822A (en) * | 1996-10-04 | 1998-02-17 | Vector Seismic Data Processing, Inc. | Seismic data radon dip moveout method |
CN102176055B (zh) * | 2011-02-18 | 2012-10-17 | 中国石油化工股份有限公司 | 保幅处理分析评价方法 |
Non-Patent Citations (3)
Title |
---|
地震资料保幅处理的讨论;郭树祥;《油气地球物理》;20090131;第7卷(第1期);全文 * |
地震资料处理中相对保幅性讨论;芮拥军;《物探与化探》;20110630;第35卷(第3期);全文 * |
地震资料处理保幅性评价方法综述与探讨;尚新民;《石油物探》;20140331;第53卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105467447A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dannemiller et al. | A new method for determination of magnetization direction | |
SA520412321B1 (ar) | تحديد تفاوتات لسرعة كثبان رملية | |
US9659252B2 (en) | Method to characterize heterogeneous anisotropic media | |
CN103954995B (zh) | 一种砂岩型铀矿勘查中的砂体识别方法 | |
US11099289B2 (en) | Multivariate analysis of seismic data, microseismic data, and petrophysical properties in fracture modeling | |
US10317554B2 (en) | Noise attenuation via thresholding in a transform domain | |
US10359529B2 (en) | Singularity spectrum analysis of microseismic data | |
CA3032780C (en) | Geostatistical analysis of microseismic data in fracture modeling | |
WO2015130441A1 (en) | Noise attenuation | |
CN106291682A (zh) | 一种基于基追踪方法的叠后声波阻抗反演方法 | |
CN107942404A (zh) | 一种确定裂缝的方法及装置 | |
Bourges et al. | A practical review of geostatistical processing applied to geophysical data: methods and applications | |
AU2016426667A1 (en) | Classifying well data using a support vector machine | |
NO20121473A1 (no) | System for modellering av geologiske strukturer | |
CN105467447B (zh) | 相控趋势能量匹配的地震保幅评价方法 | |
US9921324B2 (en) | Systems and methods employing upward beam propagation for target-oriented seismic imaging | |
CN106443791B (zh) | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 | |
Maity et al. | Reservoir characterization of an unconventional reservoir by integrating microseismic, seismic, and well log data | |
Miroslav et al. | Tutorial: The gravity-stripping process as applied to gravity interpretation in the eastern Mediterranean | |
CN105242317B (zh) | 一种纵波速度的确定方法及装置 | |
US10557957B2 (en) | System and method for improving resolution of digital seismic images | |
CN104166162B (zh) | 基于迭代三参数小波变换的缝洞发育带检测方法 | |
US10401515B2 (en) | Estimation of water properties from seismic data | |
Dubey et al. | Integrated deterministic and predictive discrete fracture network modeling for an Eocene carbonate reservoir, Bengal Basin, India | |
CN112649867B (zh) | 虚拟井构建方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |