CN113255155A - 一种储层不连续界限的阻渗能力等级的评价方法 - Google Patents

一种储层不连续界限的阻渗能力等级的评价方法 Download PDF

Info

Publication number
CN113255155A
CN113255155A CN202110654007.5A CN202110654007A CN113255155A CN 113255155 A CN113255155 A CN 113255155A CN 202110654007 A CN202110654007 A CN 202110654007A CN 113255155 A CN113255155 A CN 113255155A
Authority
CN
China
Prior art keywords
attribute
reservoir
discontinuous
boundary
discontinuous boundary
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
Application number
CN202110654007.5A
Other languages
English (en)
Other versions
CN113255155B (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.)
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
Original Assignee
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
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 Beijing Research Center of CNOOC China Ltd, CNOOC China Ltd filed Critical Beijing Research Center of CNOOC China Ltd
Priority to CN202110654007.5A priority Critical patent/CN113255155B/zh
Publication of CN113255155A publication Critical patent/CN113255155A/zh
Application granted granted Critical
Publication of CN113255155B publication Critical patent/CN113255155B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种储层不连续界限的阻渗能力等级的评价方法,包括以下步骤:建立多属性数据标准化模型;建立储层不连续界限的属性与多参数属性之间的数据栈,升维提取关联特征,并通过特征映射降维得到特征图库;根据特征矩阵w与属性图进行特征评价,得到评价等级。本发明公开的储层不连续界限的阻渗能力等级的评价方法,解决了常规不连续性预测方法地震特征弱差异性等现象导致的储层不连续界限阻渗能力等级认识不清晰的问题,同时减少了人工标注的人为干扰,提高了沉积岩性级别的不连续界限的阻渗能力等级评价的精度。

Description

一种储层不连续界限的阻渗能力等级的评价方法
技术领域
本发明涉及地质勘探技术领域,具体涉及一种储层不连续界限的阻渗能力等级的评价方法。
背景技术
储层非均质性研究贯穿了油田的整个勘探开发历程,储层的“不连续界限”是指由于沉积作用、成岩改造及构造运动等因素,在储层中形成的岩性尖灭、砂体叠置、物性变化及小断层等对流体流动产生影响的各类界限。储层不连续界限主要包括构造级别和沉积岩性级别这两大类型。
通过常规的地震属性分析技术可以准确识别构造级别的不连续界限,但对沉积岩性级别的不连续界限检测效果不佳。为了更好的识别沉积岩性级别的不连续界限,近年来借鉴于图像处理领域而来的平面算法能够准确检测到沉积岩性级别的不连续界限。然而,在实际应用中发现,这类算法所拾取的是一个各等级界限的综合,包括了岩性界限、小断层、以及砂体内部的界限,并不能体现出不同等级的界限往往具有不同程度的阻渗能力。总之,由于储层不连续界限阻渗能力等级认识不清,现有技术中缺乏能对储层不连续界限的阻渗能力等级进行有效评价的方法。
发明内容
本发明的目的在于提供一种储层不连续界限的阻渗能力等级的评价方法,用以解决现有技术中不对储层不连续界限的阻渗能力等级进行有效评价的问题。
本发明提供一种储层不连续界限的阻渗能力等级的评价方法,包括以下步骤:
步骤S1:选取储层不连续界限的属性及基于地震数据提取的阻渗参数,并对所选属性进行值域标准化处理,建立多属性数据标准化模型;
步骤S2:建立储层不连续界限的属性与多参数属性之间的数据栈,升维提取关联特征,并通过特征映射降维得到特征图库;
步骤S3:对得到的离散特征向量集合w'进行线性压缩得到特征矩阵w;采用定位数值编码技术对得到的离散特征向量集合w'进行线性压缩,得到特征矩阵w;根据特征矩阵w与属性图进行特征评价,得到评价等级。
优选地,所述步骤S1具体包括以下步骤:
步骤S11:选取不连续界限属性,包括砂体厚度A1、砂体孔隙度A2,不连续界限宽度A3,不连续界限两侧砂体高程A4以及阻渗参数In
步骤S12:分别对选取的不连续界限属性进行值域标准化处理,建立多属性数据标准化模型:
分别选取储层不连续界限属性砂体厚度A1、砂体孔隙度A2,不连续界限宽度A3,不连续界限两侧砂体高程A4,阻渗参数In进行值域标准化处理,标准化处理公式为:
Figure BDA0003111896260000021
式中,yi表示标准化处理后的属性值,xi表示原始属性值域的某个属性值;n表示原始属性值域的属性值个数;
Figure BDA0003111896260000022
表示原始属性值域xi的期望值,
Figure BDA0003111896260000023
优选地,所述步骤S2具体包括以下步骤:
步骤S21:根据储层不连续界限的属性逐点建立多参数属性数据栈,得到储层不连续界限的数据栈;
步骤S22:遍历储层不连续界限数据栈,对数据栈中数据进行升维操作,求取关联特征,按照关联特征计算公式对储层不连续界限的属性的进行特征关联,该关联特征计算公式为:
Figure BDA0003111896260000024
式中,T(x,y)表示属性x和属性y的关联特征;m表示参数属性个数;k表示储层不连续界限的属性的组合个数,0<k≤m;
步骤S23:通过特征映射的降维方法求取关联特征,通过关联特征得到特征图库。
优选地,所述步骤S23具体包括以下步骤:
将所求取的关联特征标签映射到第一空白域w1',清除原位置中对应内容,根据单变量特征选择来突出弱差异特征,弱差异特征的表达公式为:
Figure BDA0003111896260000025
式中,X2(x,y)表示属性x和属性y的差异特征;
再次映射到第二空白域w”1
重复迭代前操作,直至数据栈中内容全部出栈。
优选地,所述步骤S3具体包括以下步骤:
步骤S31:采用定位数值编码方法对离散特征向量做线性压缩处理,线性压缩表达式如下:
Figure BDA0003111896260000026
其中,N表示属性个数;px、py分别表示空白域中数据坐标;c表示当前坐标点位置的数值;
步骤S32:将属性数据与特征矩阵w进行全局视觉相似计算,特征矩阵w包含关联特征,其计算公式如下:
Figure BDA0003111896260000031
其中,C为常数;GV(i)表示全局视觉相似度,i表示像素点;
步骤S33:获取评价等级的最终结果。
与现有技术相比,本发明的有益效果是:
本发明公开了一种储层不连续界限的阻渗能力等级的评价方法,通过研究油气田储层非均质性和精细描述油藏剩余油对对沉积岩性级别的不连续界限进行等级评价,具体表现在:属性图形式特征强化处理后,使得弱差异内容突出,特征更为标准,增大了属性可用率,同时形成一体化地震数据挖掘处理流程,减少了对人工标注的依赖,加工后的特征可用于标签建立,做到特征与模型一对一;采用全局视觉相似计算作为评价等级方法,对每一个特征计算得分,根据等级分值易于判断特征是否适合所用学习方法,使得学习模型更为精准。本发明公开的储层不连续界限的阻渗能力等级的评价方法,解决了常规不连续性预测方法地震特征弱差异性等现象导致的储层不连续界限阻渗能力等级认识不清晰的问题,同时减少了人工标注的人为干扰,提高了沉积岩性级别的不连续界限的阻渗能力等级评价的精度。
具体实施方式
以下实施例用于说明本发明,但不用来限制本发明的范围。
实施例1
实施例1提供一种储层不连续界限的阻渗能力等级的评价方法,包括以下步骤:
步骤S1:选取储层不连续界限的属性及基于地震数据提取的阻渗参数,并对所选属性进行值域标准化处理,建立多属性数据标准化模型;具体包括以下步骤:
步骤S11:选取不连续界限属性,包括砂体厚度A1、砂体孔隙度A2,不连续界限宽度A3,不连续界限两侧砂体高程A4以及阻渗参数In
步骤S12:分别对选取的不连续界限属性进行值域标准化处理,建立多属性数据标准化模型:
分别选取储层不连续界限属性砂体厚度A1、砂体孔隙度A2,不连续界限宽度A3,不连续界限两侧砂体高程A4,阻渗参数In进行值域标准化处理,标准化处理公式为:
Figure BDA0003111896260000032
式中,yi表示标准化处理后的属性值,xi表示原始属性值域的某个属性值;n表示原始属性值域的属性值个数;
Figure BDA0003111896260000041
表示原始属性值域xi的期望值,
Figure BDA0003111896260000042
步骤S2:建立储层不连续界限的属性与多参数属性之间的数据栈,升维提取关联特征,并通过特征映射降维得到特征图库。
其中,“升维提取”指的是将二维的原始数据的维度提升至多维,详细步骤见下述;“特征”指的是能突显阻渗能力;“特征图库”指的是数据格式为图片,将所提取的特征集合即为特征图库。
具体包括以下步骤:
步骤S21:根据储层不连续界限的属性逐点建立多参数属性数据栈,得到储层不连续界限的数据栈;
其中,“逐点建立”指的是初始所选取的属性数据格式为图片即二维数据,图片由像素点构成,逐点建立即扫描每个像素点。
步骤S22:遍历储层不连续界限数据栈,对数据栈中数据进行升维操作,求取关联特征,按照关联特征计算公式对储层不连续界限的属性的进行特征关联,该关联特征计算公式为:
Figure BDA0003111896260000043
式中,T(x,y)表示属性x和属性y的关联特征;m表示参数属性个数;k表示储层不连续界限的属性的组合个数,0<k≤m;
这里的“升维操作”指的是对单个属性提取出的特征进行堆叠;
步骤S23:通过特征映射的降维方法求取关联特征,通过关联特征得到特征图库。具体包括以下步骤:
将所求取的关联特征标签映射到第一空白域w1',清除原位置中对应内容,根据单变量特征选择来突出弱差异特征,弱差异特征的表达公式为:
Figure BDA0003111896260000044
式中,X2(x,y)表示属性x和属性y的差异特征;
再次映射到第二空白域w”1
重复迭代前操作,直至数据栈中内容全部出栈。
步骤S3:对得到的离散特征向量集合w'进行线性压缩得到特征矩阵w;采用定位数值编码技术对得到的离散特征向量集合w'进行线性压缩,得到特征矩阵w;根据特征矩阵w与属性图进行特征评价,得到评价等级。
具体包括以下步骤:
步骤S31:采用定位数值编码方法对离散特征向量做线性压缩处理,线性压缩表达式如下:
Figure BDA0003111896260000052
其中,N表示属性个数;px、py分别表示空白域中数据坐标;c表示当前坐标点位置的数值;
步骤S32:将属性数据与特征矩阵w进行全局视觉相似计算,特征矩阵w包含关联特征,其计算公式如下:
Figure BDA0003111896260000051
其中,C为常数;GV(i)表示全局视觉相似度,i表示像素点;
步骤S33:获取评价等级的最终结果。
其中,评价等级分为3类:连通、半连通、不连通。同时,还可根据评价等级高低确定是否选用该属性。
虽然,上文中已经用一般性说明及具体实施例对本发明作了详尽的描述,但在本发明基础上,可以对之作一些修改或改进,这对本领域技术人员而言是显而易见的。因此,在不偏离本发明精神的基础上所做的这些修改或改进,均属于本发明要求保护的范围。

Claims (5)

1.一种储层不连续界限的阻渗能力等级的评价方法,其特征在于,包括以下步骤:
步骤S1:选取储层不连续界限的属性及基于地震数据提取的阻渗参数,并对所选属性进行值域标准化处理,建立多属性数据标准化模型;
步骤S2:建立储层不连续界限的属性与多参数属性之间的数据栈,升维提取关联特征,并通过特征映射降维得到特征图库;
步骤S3:对得到的离散特征向量集合w'进行线性压缩得到特征矩阵w;采用定位数值编码技术对得到的离散特征向量集合w'进行线性压缩,得到特征矩阵w;根据特征矩阵w与属性图进行特征评价,得到评价等级。
2.如权利要求1所述的储层不连续界限的阻渗能力等级的评价方法,其特征在于,所述步骤S1具体包括以下步骤:
步骤S11:选取不连续界限属性,包括砂体厚度A1、砂体孔隙度A2,不连续界限宽度A3,不连续界限两侧砂体高程A4以及阻渗参数In
步骤S12:分别对选取的不连续界限属性进行值域标准化处理,建立多属性数据标准化模型:
分别选取储层不连续界限属性砂体厚度A1、砂体孔隙度A2,不连续界限宽度A3,不连续界限两侧砂体高程A4,阻渗参数In进行值域标准化处理,标准化处理公式为:
Figure FDA0003111896250000011
式中,yi表示标准化处理后的属性值,xi表示原始属性值域的某个属性值;n表示原始属性值域的属性值个数;
Figure FDA0003111896250000012
表示原始属性值域xi的期望值,
Figure FDA0003111896250000013
3.如权利要求1所述的储层不连续界限的阻渗能力等级的评价方法,其特征在于,所述步骤S2具体包括以下步骤:
步骤S21:根据储层不连续界限的属性逐点建立多参数属性数据栈,得到储层不连续界限的数据栈;
步骤S22:遍历储层不连续界限数据栈,对数据栈中数据进行升维操作,求取关联特征,按照关联特征计算公式对储层不连续界限的属性的进行特征关联,该关联特征计算公式为:
Figure FDA0003111896250000014
式中,T(x,y)表示属性x和属性y的关联特征;m表示参数属性个数;k表示储层不连续界限的属性的组合个数,0<k≤m;
步骤S23:通过特征映射的降维方法求取关联特征,通过关联特征得到特征图库。
4.如权利要求3所述的储层不连续界限的阻渗能力等级的评价方法,其特征在于,所述步骤S23具体包括以下步骤:
将所求取的关联特征标签映射到第一空白域w′1,清除原位置中对应内容,根据单变量特征选择来突出弱差异特征,弱差异特征的表达公式为:
Figure FDA0003111896250000021
式中,X2(x,y)表示属性x和属性y的差异特征;
再次映射到第二空白域w″1
重复迭代前操作,直至数据栈中内容全部出栈。
5.如权利要求1所述的储层不连续界限的阻渗能力等级的评价方法,其特征在于,所述步骤S3具体包括以下步骤:
步骤S31:采用定位数值编码方法对离散特征向量做线性压缩处理,线性压缩表达式如下:
Figure FDA0003111896250000023
其中,N表示属性个数;px、py分别表示空白域中数据坐标;c表示当前坐标点位置的数值;
步骤S32:将属性数据与特征矩阵w进行全局视觉相似计算,特征矩阵w包含关联特征,其计算公式如下:
Figure FDA0003111896250000022
其中,C为常数;GV(i)表示全局视觉相似度,i表示像素点;
步骤S33:获取评价等级的最终结果。
CN202110654007.5A 2021-06-11 2021-06-11 一种储层不连续界限的阻渗能力等级的评价方法 Active CN113255155B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110654007.5A CN113255155B (zh) 2021-06-11 2021-06-11 一种储层不连续界限的阻渗能力等级的评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110654007.5A CN113255155B (zh) 2021-06-11 2021-06-11 一种储层不连续界限的阻渗能力等级的评价方法

Publications (2)

Publication Number Publication Date
CN113255155A true CN113255155A (zh) 2021-08-13
CN113255155B CN113255155B (zh) 2024-05-28

Family

ID=77187820

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110654007.5A Active CN113255155B (zh) 2021-06-11 2021-06-11 一种储层不连续界限的阻渗能力等级的评价方法

Country Status (1)

Country Link
CN (1) CN113255155B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6041017A (en) * 1996-08-05 2000-03-21 Fred L. Goldsberry Method for producing images of reservoir boundaries
US6865521B1 (en) * 2003-01-29 2005-03-08 Michael John Padgett Method for horizon binning for an area of interest
US20110015909A1 (en) * 2009-07-16 2011-01-20 Gang Zhao Reservoir modeling method
CN103792576A (zh) * 2014-01-28 2014-05-14 中国石油天然气股份有限公司 一种基于梯度结构张量的储层非均质检测方法及设备
CN105093303A (zh) * 2015-08-11 2015-11-25 中国海洋石油总公司 一种河流相砂体不连续界线精细表征的地震解释方法
WO2016041189A1 (zh) * 2014-09-19 2016-03-24 杨顺伟 一种评价页岩气储层及寻找甜点区的方法
CN106569272A (zh) * 2016-11-11 2017-04-19 东北石油大学 一种基于数据性质空间升维的地震属性融合方法
CN107894615A (zh) * 2017-11-13 2018-04-10 中国石油化工股份有限公司华北油气分公司勘探开发研究院 一种定量化评价三维地震属性预测储层参数有效性的方法
CN110807544A (zh) * 2019-10-08 2020-02-18 北京科技大学 一种基于机器学习的油田剩余油饱和度分布的预测方法
CN110929413A (zh) * 2019-12-04 2020-03-27 西南石油大学 一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法
CN112782760A (zh) * 2020-12-11 2021-05-11 中国海洋石油集团有限公司 一种利用地震储层不连续边界解剖辫状河储层结构的方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6041017A (en) * 1996-08-05 2000-03-21 Fred L. Goldsberry Method for producing images of reservoir boundaries
US6865521B1 (en) * 2003-01-29 2005-03-08 Michael John Padgett Method for horizon binning for an area of interest
US20110015909A1 (en) * 2009-07-16 2011-01-20 Gang Zhao Reservoir modeling method
CN103792576A (zh) * 2014-01-28 2014-05-14 中国石油天然气股份有限公司 一种基于梯度结构张量的储层非均质检测方法及设备
WO2016041189A1 (zh) * 2014-09-19 2016-03-24 杨顺伟 一种评价页岩气储层及寻找甜点区的方法
CN105093303A (zh) * 2015-08-11 2015-11-25 中国海洋石油总公司 一种河流相砂体不连续界线精细表征的地震解释方法
CN106569272A (zh) * 2016-11-11 2017-04-19 东北石油大学 一种基于数据性质空间升维的地震属性融合方法
CN107894615A (zh) * 2017-11-13 2018-04-10 中国石油化工股份有限公司华北油气分公司勘探开发研究院 一种定量化评价三维地震属性预测储层参数有效性的方法
CN110807544A (zh) * 2019-10-08 2020-02-18 北京科技大学 一种基于机器学习的油田剩余油饱和度分布的预测方法
CN110929413A (zh) * 2019-12-04 2020-03-27 西南石油大学 一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法
CN112782760A (zh) * 2020-12-11 2021-05-11 中国海洋石油集团有限公司 一种利用地震储层不连续边界解剖辫状河储层结构的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
吴佳忆;王志强;王庆宇;邵帅;孙姝婷;: "基于砂体储渗能力与集散程度的储层评价方法", 大庆石油地质与开发, no. 05, 1 October 2018 (2018-10-01) *

Also Published As

Publication number Publication date
CN113255155B (zh) 2024-05-28

Similar Documents

Publication Publication Date Title
Lorentzen et al. History matching channelized reservoirs using the ensemble Kalman filter
GB2444167A (en) Method for stratigraphic interpretation of seismic images
CN110598748B (zh) 基于卷积神经网络融合的异质图像变化检测方法及装置
CN101901343A (zh) 基于立体约束的遥感影像道路提取方法
CN113516650B (zh) 基于深度学习的电路板塞孔缺陷检测方法及装置
Canaz Sevgen et al. An improved RANSAC algorithm for extracting roof planes from airborne lidar data
CN110400371B (zh) 一种水平构造地貌实体的三维模型构建方法
CN113703045B (zh) 一种基于轻量级网络的地震相识别方法
CN112731522B (zh) 地震地层智能识别方法、装置、设备及存储介质
Borgos et al. Automated geometry extraction from 3D seismic data
CN115471467A (zh) 一种高分辨率光学遥感影像建筑物变化检测方法
CN113223170B (zh) 一种基于致密砂岩ct图像三维重构的孔隙识别方法
CN111381275A (zh) 一种地震数据的初至拾取方法和装置
CN111639067A (zh) 多特征融合卷积自编码的多元地球化学异常识别方法
CN116630802A (zh) 一种基于SwinT与尺寸自适应卷积的电力设备锈蚀缺陷图像检测方法
CN107563296B (zh) 基岩海岸岸线的提取方法及系统
CN109146888A (zh) 一种基于差商分析的土裂隙分形维数计算方法
CN114937023A (zh) 一种承压设备的缺陷识别方法及装置
CN106897683B (zh) 一种遥感图像的地物检测方法及系统
CN113255155B (zh) 一种储层不连续界限的阻渗能力等级的评价方法
US6829537B2 (en) Contour-based attribute seismic mapping
CN117593601A (zh) 一种基于深度学习的水尺验潮方法
CN107657615A (zh) 基于增量cae的高分辨sar图像变化检测方法
CN106570524A (zh) 一种储层流体类型的识别方法及装置
CN114581556B (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