CN103529473A - 一种地震反演过程中差异压实的补偿方法 - Google Patents
一种地震反演过程中差异压实的补偿方法 Download PDFInfo
- Publication number
- CN103529473A CN103529473A CN201310450291.XA CN201310450291A CN103529473A CN 103529473 A CN103529473 A CN 103529473A CN 201310450291 A CN201310450291 A CN 201310450291A CN 103529473 A CN103529473 A CN 103529473A
- Authority
- CN
- China
- Prior art keywords
- data
- wave impedance
- seismic
- sample spot
- numerical value
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种地震反演过程中差异压实的补偿方法，其特征在于利用工区内的地震测井录井数据完成地震反演，提取波阻抗平面图；在进行小层对比后，统计样品点坐标与样品点之间的深度差异数据；将沿地震解释层的波阻抗数值按数值格式导出，根据波阻抗数值的x,y大地坐标与小层对比样品点的大地坐标建立对应关系，得到一个波阻抗数值随深度增加而增加的关系式；运用关系式对导出的波阻抗数据体进行二次运算，重新平面成图，从而得到补偿压实作用之后的平面波阻抗矢量图。本发明最大限度的体现了地层由于岩性变化所引起波阻抗值的变化，避免了对地震体分块而导致的数据丢失，具有保证数据完整、保证反演工作质量和准确性、提高工作效率的特点。
Description
技术领域：
本发明涉及一种地震反演过程中差异压实的补偿方法，属地震勘探资料处理技术领域。
背景技术：
地震反演技术是充分利用钻井、测井、地质资料所提供的有关构造、岩性、层位等方面的丰富信息，由常规地震剖面来反演推导出地下地层的密度、波阻抗、速度、渗透率、孔隙度、压力、砂泥岩百分比等信息。波阻抗反演是指利用地震资料反演推导出地层波阻抗(或速度)的地震特殊处理解释技术。与地震模式识别预测油气、神经网络预测地层参数、振幅拟合预测储集层厚度等统计性方法相比,波阻抗反演处理解释技术具有十分明确的物理意义,是储集层岩性预测、油藏特征描述的确定性方法,在实际应用中取得了显著的地质效果,因此地震反演通常特指波阻抗反演。
地震反演通常分为叠前反演和叠后反演两大类。近20年来,叠后地震反演取得了巨大的进展,已形成了多种成熟技术。按测井资料在其中所起作用的大小又可分成4类:地震直接反演、测井控制下的地震反演、测井-地震联合反演和地震控制下的测井内插外推。
地震波阻抗反演主要是根据波阻抗数值即速度与密度的乘积，构造作用中的压实作用能通过影响岩石密度而影响最终的波阻抗数值，由于地势的差异性所导致的差异压实作用则使得波阻抗数值不能很好的反映出不同的岩性，所以对地震反演过程中的差异压实进行补偿则显得非常必要。
目前国内外用于地震反演中针对差异压实所采取的方法主要是将研究区进一步细分为若干个小工区，以小工区为单元来进行地震反演，但这种处理方法存在如下不足之处：1.将研究工区分为若干个小工区来进行地震反演无法在整体上准确把握区域构造与沉积背景；2.分为若干个小工区来进行地震反演的成果很难完全拼接起来而得到整个工区的整体地震反演成果，从而增加了地质学科和地球物理学科综合勘探的难度；3.分为若干个小工区来进行地震反演时所能够利用到的基础数据相比整体反演来说无疑是不足的，因此，在缺乏足够多的控制因素的前提下，其反演结果的随机性也增大，从而降低了反演成果的准确度和可信度；4.分为若干个小工区来进行地震反演，相比整体反演在反演处理的工作时间上相对要加长，工作强度上相对要增加，从而造成了资源和工作量的浪费。
发明内容：
为了克服现有技术的不足，本发明的目的在于提供一种地震反演过程中差异压实的补偿方法，最大限度的体现地层由于岩性变化所引起阻抗值的变化，从总体上把握整个研究区的构造与沉积背景，避免对地震体分块而导致的数据丢失，具有保证数据的完整性与精度、最大限度的保证反演工作的质量和准确性、提高工作效率的特点。
本发明是通过如下技术方案来实现上述目的的。
本发明提供的一种地震反演过程中差异压实的补偿方法，包括如下步骤：
（1）、提取工区内的地震体segy数据、地震解释层位数据、测井数据与录井数据；进行地震稀疏脉冲反演，得到研究区的波阻抗数据体，绘制得到地震解释层的波阻抗平面图；
（2）、收集工区内的岩心数据、录井解释数据，并进行小层对比，挑选出岩石类型相同、地层沉积厚度相同，但埋深不同的样品点；统计出样品点的大地坐标与样品点之间的深度差异数据；
（3）、依据工区内绘制的地震解释层的波阻抗平面图，按数值格式导出一组沿地震解释层的波阻抗数值，按其的x,y大地坐标与小层对比样品点的大地坐标建立对应关系；
所述的坐标建立对应关系，为波阻抗数值点的大地坐标与小层对比样品点的大地坐标两者尽可能的接近甚至重合，要求平面上距离最短，且要考虑方位角和偏移距等可能引起误差的各种因素；
（4）、根据工区研究对象的精度要求，设定补偿压实作用差值的最小计算单元，按所对应的样品点将深度差异数据与波阻抗数值差异数据进行交汇分析，得到一个波阻抗数值随深度增加而增加的关系式；
y=Ax+B （1）
其中：
Y；为样品点的波阻抗值；
x：为样品点的深度差值；
A、B：为数据统计常数；
（5）运用得到的线性关系式（1）对导出的波阻抗数据体进行二次运算，对全区不同深度点的波阻抗数据根据线性关系式（1）乘以不同系数之后得出新的波阻抗值，然后重新平面成图，从而得到补偿压实作用之后的沿地震解释层位的平面波阻抗矢量图。
本发明与现有的技术相比具有如下有益效果：
1、本发明的差异压实的补偿方法，最大限度的体现了地层由于岩性变化所引起波阻抗值的变化，排除由于构造压实作用影响岩层密度从而引起的反演结果不准确，而不能真实反映地层发育特征与分布范围。
2、本发明的差异压实的补偿方法，从总体上把握整个研究区的构造与沉积背景，避免对地震体分块而导致的数据丢失，保证了数据的完整性与精度，最大限度的保证了反演工作的质量和准确性，提高了工作效率。
具体实施方式：
下面结合具体实施例，对本发明作进一步说明。
本发明提供的一种地震反演过程中差异压实的补偿方法，包括如下步骤：
（1）、提取工区内的地震体segy数据、地震解释层位数据、测井数据与录井数据；根据测井、录井数据完成井上层位标定；根据井上层位标定结合地震反射特征完成地震解释,根据地震数据体进行地震稀疏脉冲反演；得到研究区的波阻抗数据体，根据地震解释层位对地震反演后得到的波阻抗数据体进行沿地震解释层切片，得到沿地震解释层位波阻抗平面图；绘制得到地震解释层的波阻抗平面图；
沿地震解释层提取波阻抗数值平面图，该步工作要求地震层位以及断层解释，然后合成记录制作以及时深转换公式。沿地震解释层提取波阻抗平面图中的反演工作要求较为精确的地震层位以及断层解释，然后合理的合成记录制作以及时深转换；
（2）、收集工区内的岩心数据、录井解释数据，进行小层对比；根据对取心井的观察与数据统计，挑选出岩石类型相同、地层沉积厚度相同，但埋深不同的样品点；根据样品点的选择，统计其大地坐标和样品点的埋深数据；统计出样品点的大地坐标与样品点之间的埋深数据之间的差值；
统计样品点坐标与样品点之间的深度差异数据，该步数据提取工作要求合理的小层对比方案；提取样品的时候需要在不同岩性区域，不同构造单元中提取样品点；统计样品点坐标与样品点之间的深度差异数据中的数据提取工作要求合理的小层对比方案，尽可能多的提取数据点，保证样品点的代表性与非偶然性；
（3）、依据工区内绘制的地震解释层的波阻抗平面图，按数值格式导出一组沿地震解释层的波阻抗数值，得到一组研究区沿地震解释层的包括大地坐标的波阻抗数值；根据波阻抗数值中的x,y大地坐标与小层对比样品点的大地坐标建立对应关系，计算其波阻抗数值差异；
所述的坐标建立对应关系，为波阻抗数值点的大地坐标与小层对比样品点的大地坐标两者尽可能的接近甚至重合，要求平面上距离最短，且要考虑方位角和偏移距等可能引起误差的各种因素；
（4）、根据工区研究对象的精度要求，设定补偿压实作用差值的最小计算单元，按所对应的样品点将深度差异数据与波阻抗数值差异数据进行交汇分析，得到一个关于两组差值数据之间的关系式，其波阻抗数值随深度增加而增加；关系式如下：
y=Ax+B （1）
其中：
Y；为样品点的波阻抗值；
x：为样品点的深度差值；
A、B：为数据统计常数；
为了得到地震反演波阻抗数值随深度增加而增加的关系式，需根据研究目标精度要求设定数据单元，且数据单元应该不小于研究目标要求精度；
为了得到补偿压实作用之后的平面波阻抗矢量图必须最大限度的排除由于压实作用对波阻抗数值的干扰，尽可能的反映由于岩性变化而引起的波阻抗数值变化的空间展布和变化规律；
（5）运用得到的线性关系式（1）对导出的波阻抗数据体进行二次运算，对全区不同深度点的波阻抗数据根据线性关系式（1）乘以不同系数之后得出新的波阻抗值，然后重新平面成图，从而得到补偿压实作用之后的沿地震解释层位的平面波阻抗矢量图。
Claims (2)
1.一种地震反演过程中差异压实的补偿方法，其特征在于包括如下步骤：
（1）、提取工区内的地震体segy数据、地震解释层位数据、测井数据与录井数据；进行地震稀疏脉冲反演，得到研究区的波阻抗数据体，绘制得到地震解释层的波阻抗平面图；
（2）、收集工区内的岩心数据、录井解释数据，并进行小层对比，挑选出岩石类型相同、地层沉积厚度相同，但埋深不同的样品点；统计出样品点的大地坐标与样品点之间的深度差异数据；
（3）、依据工区内绘制的地震解释层的波阻抗平面图，按数值格式导出一组沿地震解释层的波阻抗数值，按其的x,y大地坐标与小层对比样品点的大地坐标建立对应关系；
（4）、根据工区研究对象的精度要求，设定补偿压实作用差值的最小计算单元，按所对应的样品点将深度差异数据与波阻抗数值差异数据进行交汇分析，得到一个波阻抗数值随深度增加而增加的关系式；
y=Ax+B （1）
其中：
Y；为样品点的波阻抗值；
x：为样品点的深度差值；
A、B：为数据统计常数；
（5）运用得到的线性关系式（1）对导出的波阻抗数据体进行二次运算，对全区不同深度点的波阻抗数据根据线性关系式（1）乘以不同系数之后得出新的波阻抗值，然后重新平面成图，从而得到补偿压实作用之后的沿地震解释层位的平面波阻抗矢量图。
2.根据权利要求1所述的一种地震反演过程中差异压实的补偿方法，其特征在于建立波阻抗数值中的x,y大地坐标与小层对比样品点的大地坐标的对应关系，为波阻抗数值点的大地坐标与小层对比样品点的大地坐标两者尽可能的接近甚至重合，要求平面上距离最短，且要考虑方位角和偏移距等可能引起误差的各种因素。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310450291.XA CN103529473B (zh) | 2013-09-27 | 2013-09-27 | 一种地震反演过程中差异压实的补偿方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310450291.XA CN103529473B (zh) | 2013-09-27 | 2013-09-27 | 一种地震反演过程中差异压实的补偿方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103529473A true CN103529473A (zh) | 2014-01-22 |
| CN103529473B CN103529473B (zh) | 2016-07-06 |
Family
ID=49931626
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310450291.XA Expired - Fee Related CN103529473B (zh) | 2013-09-27 | 2013-09-27 | 一种地震反演过程中差异压实的补偿方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103529473B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105527653A (zh) * | 2015-12-01 | 2016-04-27 | 中国地质大学（北京） | 一种基于地质信息的虚拟井构建方法 |
| CN106405642A (zh) * | 2016-07-19 | 2017-02-15 | 西安石油大学 | 基于脱压实声波速度的地震反演储层预测方法 |
| CN111859631A (zh) * | 2020-07-01 | 2020-10-30 | 南京师范大学 | 岩浆岩体侵入过程的反演方法及装置 |
| CN113703042A (zh) * | 2020-05-21 | 2021-11-26 | 中国石油天然气股份有限公司 | 一种确定层间碎屑岩砂体累加厚度的方法及装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6058073A (en) * | 1999-03-30 | 2000-05-02 | Atlantic Richfield Company | Elastic impedance estimation for inversion of far offset seismic sections |
| EP1914568A2 (en) * | 2006-10-20 | 2008-04-23 | BHP Billiton Innovation Pty Ltd | Method for determining impedence coefficients of a seismic trace |
| CN102116871A (zh) * | 2011-01-24 | 2011-07-06 | 中国石油天然气股份有限公司 | 一种提高地震反演数据预测储层能力的方法 |
| CN102841378A (zh) * | 2011-06-23 | 2012-12-26 | 中国石油天然气股份有限公司 | 一种地震反演数据预测储层的方法 |
-
2013
- 2013-09-27 CN CN201310450291.XA patent/CN103529473B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6058073A (en) * | 1999-03-30 | 2000-05-02 | Atlantic Richfield Company | Elastic impedance estimation for inversion of far offset seismic sections |
| EP1914568A2 (en) * | 2006-10-20 | 2008-04-23 | BHP Billiton Innovation Pty Ltd | Method for determining impedence coefficients of a seismic trace |
| CN102116871A (zh) * | 2011-01-24 | 2011-07-06 | 中国石油天然气股份有限公司 | 一种提高地震反演数据预测储层能力的方法 |
| CN102841378A (zh) * | 2011-06-23 | 2012-12-26 | 中国石油天然气股份有限公司 | 一种地震反演数据预测储层的方法 |
Non-Patent Citations (3)
| Title |
|---|
| 刘成斋，等: "稀疏脉冲和基于模型反演在王家岗沙四段砂岩油藏精细勘探中的应用", 《石油地球物理勘探》, vol. 37, no. 4, 31 August 2002 (2002-08-31), pages 401 - 406 * |
| 张铭，等: "AVO约束稀疏脉冲反演技术在濮城油田储层预测中的应用", 《石油天然气学报》, vol. 34, no. 4, 30 April 2012 (2012-04-30), pages 78 - 82 * |
| 李进: "霸县凹陷文安斜坡沙二段－东营组地震反演研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》, no. 01, 15 January 2013 (2013-01-15), pages 20 - 55 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105527653A (zh) * | 2015-12-01 | 2016-04-27 | 中国地质大学（北京） | 一种基于地质信息的虚拟井构建方法 |
| CN105527653B (zh) * | 2015-12-01 | 2018-01-12 | 中国地质大学（北京） | 一种基于地质信息的虚拟井构建方法 |
| CN106405642A (zh) * | 2016-07-19 | 2017-02-15 | 西安石油大学 | 基于脱压实声波速度的地震反演储层预测方法 |
| CN113703042A (zh) * | 2020-05-21 | 2021-11-26 | 中国石油天然气股份有限公司 | 一种确定层间碎屑岩砂体累加厚度的方法及装置 |
| CN113703042B (zh) * | 2020-05-21 | 2023-09-29 | 中国石油天然气股份有限公司 | 一种确定层间碎屑岩砂体累加厚度的方法及装置 |
| CN111859631A (zh) * | 2020-07-01 | 2020-10-30 | 南京师范大学 | 岩浆岩体侵入过程的反演方法及装置 |
| CN111859631B (zh) * | 2020-07-01 | 2024-02-23 | 南京师范大学 | 岩浆岩体侵入过程的反演方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103529473B (zh) | 2016-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104977618B (zh) | 一种评价页岩气储层及寻找甜点区的方法 | |
| WO2016041189A1 (zh) | 一种评价页岩气储层及寻找甜点区的方法 | |
| CN102692645B (zh) | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 | |
| CN105759310B (zh) | 复杂非均质储层介质中地震波衰减及速度频散预测方法 | |
| CN102466815B (zh) | 三叠系碎屑岩油气藏识别方法 | |
| CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
| CN103675907A (zh) | 一种基于岩相约束的avo反演油气检测方法 | |
| CN103713319B (zh) | 一种基于地震约束建模的叠前反演方法 | |
| CN105044770A (zh) | 致密砂砾岩气藏储层定量预测方法 | |
| CN101634717A (zh) | 基于测井和叠前道集地震数据的精细横波阻抗求取技术 | |
| CN104345345A (zh) | 一种页岩储层总有机碳toc含量预测方法 | |
| CN103454685A (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
| CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
| CN105629325B (zh) | 前陆盆地冲积扇精细刻画与预测方法 | |
| CN104678434A (zh) | 一种预测储层裂缝发育参数的方法 | |
| CN104678438B (zh) | 一种co2地质封存中四维地震资料co2分布预测的方法 | |
| CN103257363A (zh) | 一种探测地下裂缝性储层中裂缝倾角的方法 | |
| CN105549082A (zh) | 超深碳酸盐储层三维地质力学场建立方法和系统 | |
| CN103777245A (zh) | 基于地震资料的油气成藏条件定量评价方法 | |
| CN105242307A (zh) | 复杂碳酸盐储层地震孔隙度获取方法及装置 | |
| CN103529473B (zh) | 一种地震反演过程中差异压实的补偿方法 | |
| CN102288996B (zh) | 叠前afi反演定量预测含气概率的方法 | |
| CN114609675A (zh) | 基于高频旋回对碳酸盐岩地层沉积微地貌的定量恢复方法 | |
| CN103487830B (zh) | 一种地震振幅类属性勘探中补偿差异压实的方法 | |
| CN103197348B (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 | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160706 Termination date: 20210927 |