CN110703354B - 一种基于井震结合的岩性边界识别方法 - Google Patents
一种基于井震结合的岩性边界识别方法 Download PDFInfo
- Publication number
- CN110703354B CN110703354B CN201911013467.9A CN201911013467A CN110703354B CN 110703354 B CN110703354 B CN 110703354B CN 201911013467 A CN201911013467 A CN 201911013467A CN 110703354 B CN110703354 B CN 110703354B
- Authority
- CN
- China
- Prior art keywords
- seismic
- well
- reservoir
- research area
- main
- 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
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
Abstract
本发明涉及石油勘探开发领域,特别是一种基于井震结合的岩性边界识别方法。该方法获取研究区的地震资料和测井资料,根据声波曲线和密度曲线制作合成地震记录,通过合成地震记录标定,开展主力油层地震响应特征分析,并且综合利用地层体切片技术、地震属性分析技术与砂体展布特征研究来识别岩性边界,先通过分析来说明研究区可以开展地震属性刻画的研究,并在主力油层的砂岩平面展布范围内进行地震属性提取,弥补了因砂泥岩速度叠置而无法开展叠后波阻抗反演进行储层预测的不足,提高了岩性边界的预测精度,降低了勘探开发成本。
Description
技术领域
本发明涉及石油勘探开发领域,特别是一种基于井震结合的岩性边界识别方法。
背景技术
随着石油勘探开发的不断深入,构造油藏越来越难发现,构造油藏依赖于对断层-岩性油藏的发现。然而,针对断层-岩性圈闭的识别,不仅需要精细的构造解释,同样,也需要开展相应的储层预测工作,但是对于部分老油区,因各种原因,有效的电子版测井资料参差不齐,尤其是准确标定储层需要用到的密度曲线缺失,进而导致相应的物探分析工作受到限制。另外,通过对有限的电子版数据进行统计分析,发现部分研究区目的层段砂泥岩速度叠置,对其进行叠后波阻抗反演难以准确预测储层边界。此外,如果仅通过开展砂体展布特征研究来预测储层边界,可能会因为部分研究区储层相变快,单砂体准确描述难度大,而导致对断层-岩性圈闭的识别造成偏差。
发明内容
本发明的目的是提供一种基于井震结合的岩性边界识别方法,用以解决现有对砂泥岩速度叠置区域进行储层边界刻画不准确的问题。
为了实现上述目的,本发明提供一种基于井震结合的岩性边界识别方法,包括以下步骤:
1)获取研究区的地震资料和测井资料,根据声波曲线和密度曲线制作合成地震记录,根据所述合成地震记录标定出所述研究区的主力油层;
2)对所述研究区开展地层对比,将对比结果与标定出的主力油层结合确定准确的主力油层;
3)对已钻井准确的主力油层开展加砂、去砂试验,判断所述研究区时间剖面上的振幅与砂体厚度是否存在正相关关系;
4)若所述研究区时间剖面上的振幅与砂体厚度存在正相关关系,则在所述研究区的任一主力油层的砂岩平面展布范围内提取地震属性,地震属性包括均方根振幅属性和瞬间Q因子属性,开展地震属性分析得到储层的岩性边界,其中,所述研究区内主力油层的砂岩平面展布范围为利用地层体切片技术预测得到。
有益效果是,通过合成地震记录标定,开展主力油层地震响应特征分析,并且综合利用地层体切片技术、地震属性分析技术与砂体展布特征研究来识别岩性边界,先通过分析来说明研究区可以开展地震属性刻画的研究,并在主力油层的砂岩平面展布范围内进行地震属性提取,弥补了因砂泥岩速度叠置而无法开展叠后波阻抗反演进行储层预测的不足,提高了岩性边界的预测精度,降低了勘探开发成本。
进一步地,对于有效的电子版测井资料参差不齐,尤其是准确标定储层需要用到的密度曲线缺失,为了避免相应的物探分析工作受到限制,提高合成地震记录质量,步骤1)中若研究区中存在缺失密度曲线的井,则根据缺失密度曲线的井的声波曲线通过似加德纳公式转换得到该井的密度曲线。
进一步地,为了准确得到实际的密度曲线,所述似加德纳公式为根据研究区的地震资料和测井资料运用最小二乘法拟合得到。
进一步地,为了降低对准确的主力油层的误判率,步骤2)中若对比结果与标定出的主力油层存在不一致,则对存在不一致的主力油层进行去砂试验,若发现该主力油层下部的波峰存在地震响应,则认为该套波峰为准确的主力油层。
附图说明
图1是本发明的研究区的地层体切片第6层图;
图2是本发明的研究区的地层体切片第7层图;
图3是本发明的研究区的地层体切片第8层图;
图4是本发明的研究区的地层体切片第13层图;
图5是本发明的研究区主力油层的砂岩等值线与构造叠合图;
图6是本发明的研究区主力油层的均方根振幅剖面属性图;
图7是本发明的研究区主力油层的瞬时Q因子剖面属性图;
图8是本发明的关键井主力油层沿层均方根振幅属性与砂岩等厚图叠合图。
具体实施方式
下面结合附图对本发明做进一步详细的说明。
本发明提供一种基于井震结合的岩性边界识别方法,包括以下步骤:
1)获取研究区的地震资料和测井资料,根据声波曲线和密度曲线制作合成地震记录,根据合成地震记录标定出研究区的主力油层。
对研究区内的井搜集测井数据,4口井有声波曲线和密度曲线,29口井只有声波曲线,有伽马曲线的只有5口。
若研究区中存在缺失密度曲线的井,则根据缺失密度曲线的井的声波曲线通过似加德纳公式转换得到该井的密度曲线。
高质量的合成地震记录需要同时用到声波曲线和密度曲线,在缺少密度曲线的情况下,密度曲线通常是用landmark软件里面的加德纳公式转换得来。但是,加德纳公式并不完全适用于我们的研究区,因此,在研究区选取了一个断层少,测井资料较为齐全,且地层相对平缓的区域,拟合了一个似加德纳公式。
研究区的似加德纳公式的原型为:
ρ=aVb
其中,ρ为研究区地层密度,V为研究区地层速度。
对工区内的井求取线性关系ρ=aVb,求取过程如下:
第一步:由ρ=aVb可推导出公式Ln(ρ)=Ln(a)+b*Ln(V);
第二步:建立上步公式的直线方程,既求解Ln(a),b的过程;
第三步:应用最小二乘法原理,将Ln(ρ)的实测值与计算值的离差的平方和的最小值做为“优化判据”,即可确定Ln(a),b;
第四步:运用EXCEL中的LINEST函数,即可求出以(Ln(V),Ln(ρ))为变量的直线方程的斜率(b)与截距Ln(a);
第五步:依据Ln(V)和LN(ρ)数列,应用最小二乘法,使用LINEST函数计算得到b=0.0792,a=1.1626。
最终,似加德纳公式如下:
ρ=1.1626V0.0792。
2)对研究区开展地层对比,将对比结果与标定出的主力油层结合确定准确的主力油层。
建立覆盖整个研究区的6条骨架剖面,并以此对全区开展地层对比及统层工作。
若对比结果与标定出的主力油层存在不一致,则对存在不一致的主力油层进行去砂试验,若发现该主力油层下部波峰出现地震响应,则认为该存在不一致的主力油层为准确的主力油层。本实施例中,通过对主力油层的去砂试验,发现其下部一套波峰存在地震响应,因此,可对这一套波峰进行追踪刻画。
3)对已钻井标定的准确的主力油层开展加砂、去砂试验,判断研究区时间剖面上的振幅与砂体厚度存在正相关关系。
本实施例中,通过对关键井(即主力油层标定的准确的已钻井)的加砂、去砂试验,发现加砂后,主力油层所在同相轴振幅有所增强,且相位关系不变;去砂后,主力油层所在同相轴振幅有所减弱,且相位关系不变。因此,分析认为,研究区时间剖面上的振幅与砂体厚度存在正相关的关系,可以通过提取振幅属性来刻画本工区的河道砂体。
4)利用地层体切片技术预测研究区内主力油层的砂岩平面展布范围,在研究区的任一主力油层的砂岩平面展布范围内提取地震属性,开展地震属性分析得到储层的岩性边界。
根据研究区的准确的主力油层,建立目的层段的层段模型,在层序界面的控制下,对夹在其间的地层采用顶底均衡的切片方式提取15个均方根振幅属性切片。如图1至图4所示,可以看出,层段内第6、7、8和13层的地层体切片,刻画了层段内不同砂体平面展布特征和演化过程。其中,第7层地层体切片较好的刻画了物源的展布特征,展布方向整体表现为由北向南展布,且形成东、西两个方向的分支砂体,其中研究区砂体主要受北东向物源控制(图5)。
获得研究区内一个主力油层的砂岩平面展布范围,在此基础上,开展地震属性分析来预测储层的岩性边界。
地震属性分析包括剖面属性分析和平面属性分析。
剖面属性分析:如图6和图7所示,通过均方根振幅属性和瞬间Q因子属性对地震剖面的分析,发现关键井的主力油层向高部位逐渐出现一个属性变化,推测可能是存在岩性变化,这与连井剖面图的分析结果基本一致。
平面属性分析:通过单层位属性提取方式,对合成地震记录和地层对比得到的准确的主力油层,沿层上下各开2ms时窗提取均方根振幅属性。如图8所示,可以看出主力油层沿层均方根振幅属性对该区砂体刻画的较清晰,边界明显,与自北向南展布的砂体基本吻合。
综合分析认为该研究区储层与北东向断层配置易形成断层-岩性圈闭。
以上给出了本发明涉及的具体实施方式,但本发明不局限于所描述的实施方式。在本发明给出的思路下,采用对本领域技术人员而言容易想到的方式对上述实施例中的技术手段进行变换、替换、修改,并且起到的作用与本发明中的相应技术手段基本相同、实现的发明目的也基本相同,这样形成的技术方案是对上述实施例进行微调形成的,这种技术方案仍落入本发明的保护范围内。
Claims (3)
1.一种基于井震结合的岩性边界识别方法,其特征在于,包括以下步骤:
1)获取研究区的地震资料和测井资料,根据声波曲线和密度曲线制作合成地震记录,根据所述合成地震记录标定出所述研究区的主力油层;
2)对所述研究区开展地层对比,将对比结果与标定出的主力油层结合确定准确的主力油层;若对比结果与标定出的主力油层存在不一致,则对存在不一致的主力油层进行去砂试验,若发现该主力油层下部的波峰存在地震响应,则认为该套波峰为准确的主力油层;
3)对已钻井准确的主力油层开展加砂、去砂试验,判断所述研究区时间剖面上的振幅与砂体厚度是否存在正相关关系;
4)若所述研究区时间剖面上的振幅与砂体厚度存在正相关关系,则在所述研究区的任一主力油层的砂岩平面展布范围内提取地震属性,地震属性包括均方根振幅属性和瞬间Q因子属性,开展地震属性分析得到储层的岩性边界,其中,所述研究区内主力油层的砂岩平面展布范围为利用地层体切片技术预测得到。
2.根据权利要求1所述的基于井震结合的岩性边界识别方法,其特征在于,步骤1)中若研究区中存在缺失密度曲线的井,则根据缺失密度曲线的井的声波曲线通过似加德纳公式转换得到该井的密度曲线。
3.根据权利要求2所述的基于井震结合的岩性边界识别方法,其特征在于,所述似加德纳公式为根据研究区的地震资料和测井资料运用最小二乘法拟合得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911013467.9A CN110703354B (zh) | 2019-10-23 | 2019-10-23 | 一种基于井震结合的岩性边界识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911013467.9A CN110703354B (zh) | 2019-10-23 | 2019-10-23 | 一种基于井震结合的岩性边界识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110703354A CN110703354A (zh) | 2020-01-17 |
CN110703354B true CN110703354B (zh) | 2021-11-16 |
Family
ID=69201184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911013467.9A Active CN110703354B (zh) | 2019-10-23 | 2019-10-23 | 一种基于井震结合的岩性边界识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110703354B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111340379B (zh) * | 2020-03-02 | 2022-06-17 | 西南石油大学 | 一种稀井网区深海水道砂体建筑结构的解剖方法 |
CN111505720B (zh) * | 2020-05-21 | 2023-03-28 | 中国海洋石油集团有限公司 | 一种岩性圈闭的刻画方法 |
CN111736216B (zh) * | 2020-05-25 | 2023-06-16 | 恒泰艾普集团股份有限公司 | 一种基于地震的油水界面地震层位确定方法 |
CN112180465B (zh) * | 2020-11-03 | 2022-06-21 | 中国石油化工股份有限公司 | 一种井震联合的地层展布确定方法 |
CN113514883B (zh) * | 2021-06-18 | 2023-03-17 | 中国石油化工股份有限公司 | 一种断层-岩性油藏刻画方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330828A (zh) * | 2014-10-27 | 2015-02-04 | 中国石油天然气股份有限公司 | 甜点储层的预测方法及预测装置 |
CN105334535A (zh) * | 2015-11-13 | 2016-02-17 | 中国石油化工股份有限公司 | 一种识别薄储层隐蔽岩性油气藏的方法 |
CN106371139A (zh) * | 2016-08-16 | 2017-02-01 | 中国石油化工股份有限公司 | 一种应用层段振幅属性识别岩性油藏的方法及装置 |
CN107121699A (zh) * | 2017-02-28 | 2017-09-01 | 中国石油化工股份有限公司 | 一种地震相控制下的沉积微相识别方法 |
CN107728230A (zh) * | 2017-09-11 | 2018-02-23 | 中国石油大港油田勘探开发研究院 | 一种基于地球物理三联技术的岩性油藏预测方法 |
CN109212627A (zh) * | 2017-07-06 | 2019-01-15 | 中国石油化工股份有限公司 | 一种基于细分层系的薄储层描述及增储规模确定方法 |
CN109441438A (zh) * | 2018-10-18 | 2019-03-08 | 中国石油化工股份有限公司 | 一种构造复杂区岩性边界表征方法 |
-
2019
- 2019-10-23 CN CN201911013467.9A patent/CN110703354B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330828A (zh) * | 2014-10-27 | 2015-02-04 | 中国石油天然气股份有限公司 | 甜点储层的预测方法及预测装置 |
CN105334535A (zh) * | 2015-11-13 | 2016-02-17 | 中国石油化工股份有限公司 | 一种识别薄储层隐蔽岩性油气藏的方法 |
CN106371139A (zh) * | 2016-08-16 | 2017-02-01 | 中国石油化工股份有限公司 | 一种应用层段振幅属性识别岩性油藏的方法及装置 |
CN107121699A (zh) * | 2017-02-28 | 2017-09-01 | 中国石油化工股份有限公司 | 一种地震相控制下的沉积微相识别方法 |
CN109212627A (zh) * | 2017-07-06 | 2019-01-15 | 中国石油化工股份有限公司 | 一种基于细分层系的薄储层描述及增储规模确定方法 |
CN107728230A (zh) * | 2017-09-11 | 2018-02-23 | 中国石油大港油田勘探开发研究院 | 一种基于地球物理三联技术的岩性油藏预测方法 |
CN109441438A (zh) * | 2018-10-18 | 2019-03-08 | 中国石油化工股份有限公司 | 一种构造复杂区岩性边界表征方法 |
Non-Patent Citations (4)
Title |
---|
井震结合识别砂体上倾尖灭岩性圈闭;柯钦 等;《中国地球科学联合学术年会2018》;20181231;第2012-2013页 * |
利用地震属性预测河道砂体;关达 等;《石油物探》;20040430;第43卷;第56-58页 * |
岔河集构造带综合研究及预探目标评价;朱正平等;《内蒙古石油化工》;20051130(第11期);第83页 * |
泌阳凹陷杜坡地区断层-岩性油藏识别技术研究;刘娟 等;《石油地质与工程》;20170331;第31卷(第2期);第41-43、47页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110703354A (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110703354B (zh) | 一种基于井震结合的岩性边界识别方法 | |
CN100557464C (zh) | 基于叠前波场模拟的地震勘探层位标定方法 | |
CN111596364B (zh) | 基于高精度层序地层格架的地震沉积微相组合分析方法 | |
CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
Noorlandt et al. | Characterisation of ground motion recording stations in the Groningen gas field | |
CN111399056B (zh) | 一种基于划分方位角滤波预测裂缝强度的方法 | |
CN103592680B (zh) | 一种基于正反演的测井数据和深度域地震剖面合成方法 | |
CN104570067A (zh) | 一种地球物理勘探中相控地震反演方法 | |
Kianoush et al. | Compressional and shear interval velocity modeling to determine formation pressures in an oilfield of SW Iran | |
CN110954944A (zh) | 一种断层圈闭含油高度地震预测方法 | |
CN103777245A (zh) | 基于地震资料的油气成藏条件定量评价方法 | |
RU2145100C1 (ru) | Способ поиска, разведки и исследования нефтегазовой залежи | |
Chen et al. | Geophysical interpretation of a subsurface landslide in the Southern Qinshui Basin | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
CN103628866A (zh) | 获取地层的地层压力系数的方法 | |
Dickmann | 3D Tunnel Seismic Prediction: A next generation tool to characterize rock mass condition ahead of the tunnel face | |
Bang et al. | Mean refracted ray path method for reliable downhole seismic data interpretations | |
CN112394392B (zh) | 对烃源岩分布情况进行评价的方法和装置 | |
CN113279746A (zh) | 一种套管变形风险区域的确定方法及应用 | |
CN113589365A (zh) | 基于时频域信息的储层尖灭线描述方法 | |
CN112180464A (zh) | 一种储层物性的识别方法 | |
Song et al. | The impact of local velocity anomaly on structure confirming and the varying velocity mapping strategy | |
CN112180443A (zh) | 页岩气二维地震甜点区优选方法及装置 | |
Gargouri | Multicomponent 3D seismic interpretation of the Marcellus shale Bradford county, Pennsylvania | |
CN110795513A (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 |