CN106968668B - 一种得分式测井岩性优选解释方法 - Google Patents

一种得分式测井岩性优选解释方法 Download PDF

Info

Publication number
CN106968668B
CN106968668B CN201710220259.0A CN201710220259A CN106968668B CN 106968668 B CN106968668 B CN 106968668B CN 201710220259 A CN201710220259 A CN 201710220259A CN 106968668 B CN106968668 B CN 106968668B
Authority
CN
China
Prior art keywords
lithology
logging
different
interpretation
database
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.)
Expired - Fee Related
Application number
CN201710220259.0A
Other languages
English (en)
Other versions
CN106968668A (zh
Inventor
刘敬寿
丁文龙
杨海盟
周学慧
肖子亢
黄昌杰
谷阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201710220259.0A priority Critical patent/CN106968668B/zh
Publication of CN106968668A publication Critical patent/CN106968668A/zh
Application granted granted Critical
Publication of CN106968668B publication Critical patent/CN106968668B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及油气田勘探开发领域，尤其是一种得分式测井岩性优选解释方法。本发明专利通过岩心归位、成像测井识别以及阵列声波测井等建立不同的常规测井曲线与岩性的识别标准数据库；依据每种岩性的有效判别数据库的大小，建立岩性判别顺序，在不同阈值条件下，分别计算不同岩性的测井识别数据库对应的得分率；通过循环迭代，输出不同井、不同层段测井解释结果。本发明专利解决了传统的测井解释软件往往解释结果单一，在不同的地区、不同岩性段适用性差的问题；它实现了测井岩性的多元合理解释，为研究人员筛选解释结果提供了合理的依据。

Description

一种得分式测井岩性优选解释方法
技术领域
本发明涉及油气田勘探开发领域，尤其是一种得分式测井岩性优选解释方法。
背景技术
测井解释岩性具有多解性，不同的解释人员的解释结论往往不一，这取决于研究者对工区的了解程度、专业背景等等；传统的测井解释软件往往解释结果单一，给研究人员选择的余地较小，很难兼顾地质资料的特殊性以及科研人员的专业知识背景，测井软件在不同的地区、不同岩性段适用性差。如何充分利用岩心、成像测井识别以及阵列声波测井等较为准确的岩性解释结论，同样也是测井解释软件面临的一大难题。本发明专利通过岩心归位、成像测井识别以及阵列声波测井等建立不同的常规测井曲线与岩性的识别标准数据库；依据每种岩性的有效判别数据库的大小，建立岩性判别顺序，在不同阈值条件下，分别计算不同岩性的测井识别数据库对应的得分率；通过循环迭代，输出不同井、不同层段测井解释结果。
发明内容
本发明旨在解决上述问题，提供了一种得分式测井岩性优选解释方法，它实现了测井岩性的优选解释。
本发明的技术方案为：一种得分式测井岩性优选解释方法，具体步骤如下：
第一步 建立岩性识别数据库。
通过岩心归位、成像测井识别以及阵列声波测井建立不同的常规测井曲线与岩性的识别标准数据库；所述的常规测井曲线与岩性的识别标准数据库是指建立不同的岩性对应的测井曲线响应数据库，利用每次岩心归位、成像测井识别以及阵列声波测井解释结果，建立岩性与对应的一套测井相应数据库一次响应，进而建立不同岩性对应的测井识别数据库：W1、W2、...Wi...Wm-1、Wm，其中，1≤i≤m；m为岩性的数目；对于第i类岩性，共有ni组有效判别数据库，标记为：Yi，1、Yi，2、...Wi，j...Yi，ni-1、Yi，ni；对于第i类岩性的第j组有效判别数据库，共有qi,j条测井曲线判别条件，依次标记为：Xi，j，1、Xi，j，2、...Xi，j，k...Xi，j，qi,j-1、Xi，j，qi,j(图2)。
第二步 确定岩性判别顺序。
依据每种岩性的有效判别数据库ni的大小，建立岩性判别顺序，优选的是，ni数值大者的优先识别为该种岩性。
第三步 测井岩性输出阈值设定(得分标准的判断)。
所述的岩性输出阈值是指对于第i类岩性的第j组有效判别数据库的qi,j条测井曲线，确定测井岩性输出阈值(Φu)0≤Φu≤1；例如，u＝5时，设定为Φ1＝0.2，Φ2＝0.4，Φ3＝0.6，Φ4＝0.8，Φ5＝1，在测井解释过程中，对于某一岩性段共有s条测井曲线测井响应值，某种岩性的测井解释得分标准记为：
Figure GDA0002417420040000021
ei,j为测井解释过程中对于第i类岩性的第j组有效判别数据库，满足qi,j条测井曲线判别条件的曲线条数。
对于不同的阈值(Φu)：
ηi≥Φu (2)
当ηi满足公式(2)时，即在阈值(Φu)条件下，该某一岩性段的第i类岩性解释可能性加一分。
第四步 不同阈值条件下，岩性得分率计算。
在不同阈值条件下，分别计算测井识别数据库：W1、W2、...Wi...Wm-1、Wm对应的得分率：
Figure GDA0002417420040000022
εu,i为第u阈值对应的i岩性得分率；依次计算不同阈值对应的不同岩性得分率，在同一阈值下，得分最高的判定为该阈值下测井岩性识别结果；优选的是，相同的得分率时优先识别为ni数值大的所对应的岩性。
第五步循环迭代，输出不同井、不同层段测井解释结果。
通过阈值(Φu)的循环迭代，输出不同井、不同层段测井解释结果。
本发明的有益效果是：本发明专利通过岩心归位、成像测井识别以及阵列声波测井等建立不同的常规测井曲线与岩性的识别标准数据库；依据每种岩性的有效判别数据库的大小，建立岩性判别顺序，在不同阈值条件下，分别计算不同岩性的测井识别数据库对应的得分率；通过循环迭代，输出不同井、不同层段测井解释结果。它实现了测井岩性的多元合理解释，为研究人员筛选解释结果提供了合理的依据。
附图说明
图1为一种得分式测井岩性优选解释方法的流程图。
图2为不同岩性对应的测井识别数据库。
图3为研究区构造位置图。
图4研究区某丼不同阈值条件的岩性解释结果。
具体实施方式
下面结合附图说明本发明的具体实施方式：
以鄂尔多斯地区东仁沟区块为例说明本发明的具体实施方式，研究区区域构造处于鄂尔多斯盆地伊陕斜坡带的西部边缘，构造简单，主要为西倾单斜背景上由差异压实作用形成的一系列由东向西倾没的低幅鼻状隆起，鼻状隆起轴线近于东西向(图3)。东西向构造与砂体有效配合，形成良好的鼻隆背斜型圈闭，有利于油气聚集。从长7砂体顶面构造等值线图上看，该区为西倾单斜背景上的差异压实成因的鼻隆背斜型圈闭构造，全区构造相对平缓，沿倾斜方向坡降10m/km左右。
东仁沟区块定4115井区长7地层厚度一般85～110m左右。岩性主要为油页岩、暗色泥岩、砂质泥岩夹黑灰色粉细砂岩。长7总体表现为顶部和底部各发育一套油页岩或黑色泥岩，中部发育两套砂体，砂体间为一层砂泥岩互层。砂层总体表现为由细到粗的反韵律性为主，局部有正韵律。砂岩中主要发育平行层理、板状斜层理。局部发育沙纹层理、透镜状层理、块状构造和斜层理。水平井开发目的层长71为本区新发现的油层，为一单砂体，一般砂体厚度10～25m。研究区长7砂岩岩石类型为细、细—中粒灰绿色岩屑长石砂岩、长石岩屑砂岩为主，粒径粒度中值平均0.10mm，粒度范围一般0.1～0.25mm，最大粒径0.25～0.45mm。颜色大都呈灰色、灰绿色、褐灰色及灰白色，薄层状、中层状、中厚层状至块状均可见。砂岩中的石英含量平均为24.7％；长石含量平均46.4％；岩屑含量平均16.8％。根据研究区储层取心井岩心样品分析物性资料及测井孔渗解释结果统计，长7砂岩储层孔隙度最大值14.6％，最小值6.2％，平均10.3％；渗透率最大4.1mD，最小0.2mD，平均0.54mD。研究区长7水平井实施区砂体呈北东南西向展布，区内砂体发育，连片性较好，砂体厚度变化不大，一般15～25m。河道中心砂体发育最好，局部砂体厚度大于30m。研究区油层发育情况主要受砂体分布影响，其发育特征也与砂体发育特征类似。油层发育与砂体匹配关系较好，一般在河道中心砂体厚度大、物性好的地方，油层也较发育，有效厚度一般10～15m，局部达20m以上。
研究区测井岩性解释流程为：
第一步 建立岩性识别数据库。通过岩心归位、成像测井识别以及阵列声波测井等建立不同的常规测井曲线与岩性的识别标准数据库；所述的常规测井曲线与岩性的识别标准数据库是指建立不同的岩性对应的测井曲线响应数据库，利用每次岩心归位、成像测井识别以及阵列声波测井解释结果，建立岩性与对应的一套测井相应数据库一次响应，进而建立不同岩性对应的测井识别数据库：W1、W2、...Wi...Wm-1、Wm，其中，1≤i≤7；对于第i类岩性，共有ni组有效判别数据库，标记为：Yi，1、Yi，2、...Wi，j...Yi，ni-1、Yi，ni；对于第i类岩性的第j组有效判别数据库，共有qi,j条测井曲线判别条件，依次标记为：Xi，j，1、Xi，j，2、...Xi，j，k...Xi，j，qi,j-1、Xi，j，qi,j(图2).
第二步 确定岩性判别顺序。
依据每种岩性的有效判别数据库ni的大小，建立岩性判别顺序，优选的是，ni数值大者的优先识别为该种岩性；依次判别：油页岩、凝灰岩、泥岩、粉砂岩、细砂岩、钙质细砂岩、中砂岩。
第三步 测井岩性输出阈值设定。
所述的岩性输出阈值是指对于第i类岩性的第j组有效判别数据库的qi,j条测井曲线，确定测井岩性输出阈值u＝4，设定为Φ1＝0.2，Φ2＝0.4，Φ3＝0.6，Φ4＝0.8在测井解释过程中，对于某一岩性段共有s条测井曲线测井响应值，利用公式(1)-(2)确定该某一岩性段的得分。
第四步 不同阈值条件下，岩性得分率计算。
在不同阈值条件下，利用公式(3)，得分最高的判定为该阈值下测井岩性识别结果；优选的是，相同的得分率时依据岩性判别顺序确定对应的岩性。
第五步 循环迭代，输出不同井、不同层段测井解释结果。
通过阈值(Φu)的循环迭代，输出不同井、不同层段测井解释结果(图4)。
上面以举例方式对本发明进行了说明，但本发明不限于上述具体实施例，凡基于本发明所做的任何改动或变型均属于本发明要求保护的范围。

Claims (1)

1.一种得分式测井岩性优选解释方法，步骤如下：
第一步 建立岩性识别数据库；
通过岩心归位、成像测井识别以及阵列声波测井建立不同的常规测井曲线与岩性的识别标准数据库；所述的常规测井曲线与岩性的识别标准数据库是指建立不同的岩性对应的测井曲线响应数据库，利用每次岩心归位、成像测井识别以及阵列声波测井解释结果，建立岩性与对应的一套测井相应数据库一次响应，进而建立不同岩性对应的测井识别数据库：W1、W2、...Wi...Wm-1、Wm，其中，1≤i≤m；m为岩性的数目；对于第i类岩性，共有ni组有效判别数据库，标记为：Yi，1、Yi，2、...Wi，j...Yi，ni-1、Yi，ni；对于第i类岩性的第j组有效判别数据库，共有qi,j条测井曲线判别条件，依次标记为：Xi，j，1、Xi，j，2、...Xi，j，k...Xi，j，qi,j-1、Xi，j，qi,j；
第二步 确定岩性判别顺序；
依据每种岩性的有效判别数据库ni的大小，建立岩性判别顺序，ni数值大者的优先识别为该种岩性；
第三步 测井岩性输出阈值设定；
所述的岩性输出阈值是指对于第i类岩性的第j组有效判别数据库的qi,j条测井曲线，确定测井岩性输出阈值Φu；所述的Φu满足0≤Φu≤1；在测井解释过程中，对于某一岩性段共有s条测井曲线测井响应值，某种岩性的测井解释得分标准记为：
Figure FDA0002417420030000011
ei,j为测井解释过程中对于第i类岩性的第j组有效判别数据库，满足qi,j条测井曲线判别条件的曲线条数；
对于不同的阈值Φu：
ηi≥Φu (2)
当ηi满足公式(2)时，即在阈值Φu条件下，该某一岩性段的第i类岩性解释可能性加一分；
第四步 不同阈值条件下，岩性得分率计算；
在不同阈值条件下，分别计算测井识别数据库：W1、W2、...Wi...Wm-1、Wm对应的得分率：
Figure FDA0002417420030000012
εu,i为第u阈值对应的i岩性得分率；依次计算不同阈值对应的不同岩性得分率，在同一阈值下，得分最高的判定为该阈值下测井岩性识别结果；相同的得分率时，依据岩性判别顺序判别岩性；
第五步 循环迭代，输出不同井、不同层段测井解释结果。
CN201710220259.0A 2017-04-06 2017-04-06 一种得分式测井岩性优选解释方法 Expired - Fee Related CN106968668B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710220259.0A CN106968668B (zh) 2017-04-06 2017-04-06 一种得分式测井岩性优选解释方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710220259.0A CN106968668B (zh) 2017-04-06 2017-04-06 一种得分式测井岩性优选解释方法

Publications (2)

Publication Number Publication Date
CN106968668A CN106968668A (zh) 2017-07-21
CN106968668B true CN106968668B (zh) 2020-08-14

Family

ID=59337366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710220259.0A Expired - Fee Related CN106968668B (zh) 2017-04-06 2017-04-06 一种得分式测井岩性优选解释方法

Country Status (1)

Country Link
CN (1) CN106968668B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108227036A (zh) * 2018-01-22 2018-06-29 中国石油大港油田勘探开发研究院 一种细粒沉积岩岩心归位的方法
CN113945992B (zh) * 2020-07-15 2024-06-04 中国石油化工股份有限公司 泥岩、油页岩识别方法、装置、电子设备及介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040158406A1 (en) * 2003-01-31 2004-08-12 Harrison Christopher J. Method for computing complexity, confidence and technical maturity indices for reservoir evaluations
US20110231164A1 (en) * 2010-03-18 2011-09-22 Schlumberger Technology Corporation Generating facies probablity cubes
CN103529474A (zh) * 2012-07-02 2014-01-22 中国石油化工股份有限公司 采用岩性细分实现岩相精细描述的方法
CN104632206A (zh) * 2015-02-05 2015-05-20 中国石油天然气股份有限公司 一种石油测井标准层确定方法及装置
CN104914482A (zh) * 2014-03-13 2015-09-16 中国石油化工股份有限公司 一种复杂砂砾岩岩相组合类型定量识别方法
CN104989392A (zh) * 2015-07-10 2015-10-21 中国石油天然气股份有限公司 一种岩性识别方法
CN105201492A (zh) * 2015-09-07 2015-12-30 山东科技大学 综合识别巨厚煤层内部沉积间断面的方法
CN105201490A (zh) * 2014-06-05 2015-12-30 中国石油化工股份有限公司 一种页岩层段岩石相分析方法
CN105697002A (zh) * 2014-11-24 2016-06-22 中国石油化工股份有限公司 一种用于识别煤系地层岩性的方法
CN105986811A (zh) * 2015-02-11 2016-10-05 中国石油化工股份有限公司 基于常规测井资料五维度岩性判识方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040158406A1 (en) * 2003-01-31 2004-08-12 Harrison Christopher J. Method for computing complexity, confidence and technical maturity indices for reservoir evaluations
US20110231164A1 (en) * 2010-03-18 2011-09-22 Schlumberger Technology Corporation Generating facies probablity cubes
CN103529474A (zh) * 2012-07-02 2014-01-22 中国石油化工股份有限公司 采用岩性细分实现岩相精细描述的方法
CN104914482A (zh) * 2014-03-13 2015-09-16 中国石油化工股份有限公司 一种复杂砂砾岩岩相组合类型定量识别方法
CN105201490A (zh) * 2014-06-05 2015-12-30 中国石油化工股份有限公司 一种页岩层段岩石相分析方法
CN105697002A (zh) * 2014-11-24 2016-06-22 中国石油化工股份有限公司 一种用于识别煤系地层岩性的方法
CN104632206A (zh) * 2015-02-05 2015-05-20 中国石油天然气股份有限公司 一种石油测井标准层确定方法及装置
CN105986811A (zh) * 2015-02-11 2016-10-05 中国石油化工股份有限公司 基于常规测井资料五维度岩性判识方法
CN104989392A (zh) * 2015-07-10 2015-10-21 中国石油天然气股份有限公司 一种岩性识别方法
CN105201492A (zh) * 2015-09-07 2015-12-30 山东科技大学 综合识别巨厚煤层内部沉积间断面的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
应用测井资料自动确定沉积相的计算机程序;肖义越等;《地质科学》;19930131;第28卷(第1期);第68-78页 *

Also Published As

Publication number Publication date
CN106968668A (zh) 2017-07-21

Similar Documents

Publication Publication Date Title
CN105201490B (zh) 一种页岩层段岩石相分析方法
CN104991274A (zh) 多种地质因素定量约束下的单一圈闭级的有利区优选方法
CN103437760B (zh) 一种使用阵列感应数据快速评价油水层的方法
CN113283108B (zh) 一种定量评价页岩油储层可压裂性的方法及系统
CN110837117B (zh) 一种含油气盆地中洼陷的综合评价方法
CN106968668B (zh) 一种得分式测井岩性优选解释方法
CN111624677B (zh) 一种致密砂岩气藏单期河道砂体边界刻画方法
CN113608278A (zh) 沉积盆地红杂色层中的砂岩型铀矿定位方法
CN116122801A (zh) 一种页岩油水平井体积压裂可压性综合评价方法
CN109555518B (zh) 一种基于聚类和判别算法的冲积扇单井构型识别方法
CN112034526A (zh) 基于岩相组合的灰质泥岩发育区薄层浊积砂体的地震识别方法
CN110656922B (zh) 基于笔石带特征的页岩等时地层测井划分方法及系统
CN105089659B (zh) 一种砾岩油藏渗流单元识别方法
CN108241785A (zh) 一种非均质性储层饱和度场精细表征方法
CN106285659B (zh) 一种低渗岩性非均质气藏开发单元划分方法
CN115201912B (zh) 基于聚类分析的岩石矿物体积含量曲线转换岩性剖面方法
CN111198406A (zh) 一种红层储层的因子分析测井岩性识别方法
Raasch Upper Paleozoic Section at Highwood Pass, Alberta: Carboniferous
CN113673771B (zh) 一种页岩气水平井压裂分段方法
CN110888164B (zh) 一种确定三角洲前缘席状砂储层砂体尖灭点的方法
CN114488298A (zh) 一种预测微生物岩甜点储层的方法
Garlough et al. Hugoton Gas Field, Grant, Haskell, Morton, Stevens, and Seward Counties, Kansas, and Texas County, Oklahoma
CN114076992A (zh) 基于岩相的曲流河储层微观孔喉分类评价方法
CN105569650A (zh) 基于马尔科夫链数学模型的井间沉积序列研究方法
Srivastava et al. Palynology of Kamthi Formation in Godavari Graben

Legal Events

Date Code Title Description
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20190805

Address after: 266580 Qingdao Changjiang Road, Huangdao District, Shandong, No. 66

Applicant after: China Petroleum University (East China)

Address before: 100083 Xueyuan Road, Beijing, China University of Geosciences, (Beijing), 29

Applicant before: China University of Geosciences (Beijing)

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200814

Termination date: 20210406

CF01 Termination of patent right due to non-payment of annual fee