CN111983680A - 基于嵌入三维卷积算子的小尺度岩性边界刻画方法 - Google Patents
基于嵌入三维卷积算子的小尺度岩性边界刻画方法 Download PDFInfo
- Publication number
- CN111983680A CN111983680A CN202010894942.4A CN202010894942A CN111983680A CN 111983680 A CN111983680 A CN 111983680A CN 202010894942 A CN202010894942 A CN 202010894942A CN 111983680 A CN111983680 A CN 111983680A
- Authority
- CN
- China
- Prior art keywords
- gradient
- scale
- small
- convolution operator
- edge
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000000694 effects Effects 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 230000002787 reinforcement Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 4
- 238000011161 development Methods 0.000 abstract description 3
- 239000003208 petroleum Substances 0.000 abstract description 2
- 238000007670 refining Methods 0.000 abstract description 2
- 239000004576 sand Substances 0.000 description 4
- 238000003708 edge detection Methods 0.000 description 3
- 230000002708 enhancing effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于嵌入三维卷积算子的小尺度岩性边界刻画方法,属于石油地震勘探开发技术领域。首先利用设计的X、Y、Z三个方向的三维卷积算子求取变化梯度数据体;在变化梯度数据体的基础上,再利用局部增强算法对微弱的岩性变化带实现局部增强,放大数据体内部微小的岩性变化信息;最后利用边缘锐化算法实现细化边缘效果,使边缘信息更精细。利用该方法有效提取数据体内部小尺度岩性变化细节,实现对地下地质储层小尺度的河道识别、三级构型研究及后续精细建模。
Description
技术领域
本发明属于石油地震勘探开发技术领域,具体涉及一种利用嵌入三维卷积算子的小尺度岩性边界刻画方法,增强对岩性边界、小尺度河道、废弃河道响应特征的方法。
背景技术
地震信息是一种大尺度的信息,对垂向厚度小于5米的地质体刻画较为困难。随着油田精细化开发生产,常规的刻画砂体边界信息的属性已经不能满足生产研究的需要,因此准确识别小尺度点坝边缘、砂体内部的末期河道水线、废弃河道的位置对砂体精细研究具有重要意义。常用的最小距离法、协方差法等方法对断裂体系的刻画、小尺度断层的研究具有明显效果,但对点坝边部的岩性变化点、末期河道水线、废弃河道的刻画显得无能为力。
本发明提出的三维卷积算子实现对地震数据体X、Y、Z三个方向求取地震数据体内部变化梯度,然后利用边缘增强、边缘细化等方法实现对微弱异常信息的增强处理,达到刻画微小地质体或砂体内部微小变化的目的。
发明内容
为了解决现有技术的问题,本发明提供了一种基于嵌入三维卷积算子的小尺度岩性边界刻画方法,解决现有技术中对岩性边界、小尺度河道、废弃河道的地质体刻画较为困难的问题。
本发明技术方案为:
一种基于嵌入三维卷积算子的小尺度岩性边界刻画方法,包括如下步骤:
(1)设计一种三维卷积算子,实现对地震数据体X、Y、Z三个方向求取变化梯度,并按照一定的比例因子求取均方根。
(2)设计一种三维数据体局部梯度增强算法,实现对步骤(1)中得到的梯度数据体进行局部梯度较弱的细节增强。
(3)设计一种锐化异常边缘效果的算法,实现对步骤(2)中得到数据体局部异常边缘精细化。
(4)选取合适截止值对数据体二值化,得到小尺度岩性边缘数据体。
所述步骤(1)中所用计算方法如下:
其中z1、z2、z3、……z26、z27为3x3x3的三维时窗里的数值,其各自空间位置如附图所示,为三维时窗内X方向的变化梯度,为三维时窗内Y方向的变化梯度,为三维时窗内Z方向的变化梯度,G(x,y,z)为卷积运算后的新数据体,m、n、r值域范围[0,1]。
所述步骤(2)中,所用计算方法如下:
E(x,y,z)为梯度边缘增强体;M[G(x,y,z)]为梯度体窗口内中值;sd为梯度体窗口内标准差;t为调节参数,值域范围[1,10]。
所述步骤(3)中,所用计算方法如下:
R(x,y,z)为边缘细化后的数据体,k为系数参数,值域范围[0.5,0.9],依据实际地震资料进行调整。
附图说明
图1为本发明三维卷积算子模型;(a)为原始地震数据;(b)为X方向算子;(c)为Y方向算子;(d)为Z方向算子;
图2中(a)为常规算法边缘检测切片案例一;(b)为案例一的本发明算法边缘检测切片效果;(c)为常规算法边缘检测切片案例二;(d)为案例二的本发明算法边缘检测切片效果。
具体实施方式
下面对本发明进行详细的说明。具体实施方式包括以下步骤:
(1)对原始地震资料利用三维卷积算子求取X、Y、Z三个方向的变化梯度,并利用公式1计算出一个梯度体,其中调节参数m、n、r根据实际地震资料情况和地下小尺度地质体展布特征进行调节,本次应用中均取0.5,应用效果最佳。
(2)利用公式2对步骤(1)中所得梯度体进行局部增强,其中M[G(x,y)]为梯度体窗口内中值;sd为梯度体窗口内标准差;t为调节参数,其值域范围为[1,10],根据不同地震资料选取不同值,本次应用中取值7,应用效果最佳;
(3)利用公式3对步骤(2)中所得增强梯度体进化边缘锐化处理,使得异常边界更加精细,其中k为系数参数,值域范围为[0.5,0.9],依据实际地震资料进行调整,本次应用中,根据实际地震资料进行参数调节,确定选取0.7时效果最佳。
(4)步骤4为灰度处理模式,可以根据实际需要进行运算。
所述步骤(1)中所用计算方法如下:
其中m、n、r值域范围[0,1],本次应用中三个参数均取值0.5,效果最佳。
所述步骤(2)中所用计算方法如下:
E(x,y,z)为梯度边缘增强体;M[G(x,y,z)]为梯度体窗口内中值;sd为梯度体窗口内标准差;t为调节参数,值域范围(1,10),本次应用中取值7,效果最佳。
所述步骤(3)中,所用计算方法如下:
其中k为系数参数,值域范围[0.5,0.9],依据实际地震资料进行调整,本次应用中取值0.7,效果最佳。
本发明基于嵌入三维卷积算子的小尺度岩性边界刻画方法,实现地震数据体内部小尺度信息刻画,达到识别点坝边界、小尺度河道的目的。首先利用设计的X、Y、Z三个方向的三维卷积算子求取地震数据体的变化梯度体;在变化梯度体基础上,再利用局部增强算法对微弱的岩性变化带实现局部增强,放大数据体内部微小的岩性变化信息;最后利用边缘锐化算法实现细化边缘效果,使边缘信息更精细。利用该方法有效提取数据体内部小尺度岩性变化细节,实现对地下地质储层小尺度的河道识别、三级构型研究及精细建模。
上述实例对本发明做了进一步描述,各方法的实施步骤都可以修改,凡是在本发明上述技术方案的基础上做出的各种修改,都属于本发明所保护的范围。
Claims (4)
1.一种基于嵌入三维卷积算子的小尺度岩性边界刻画方法,其特征在于所述方法包括以下几个步骤:
(1)设计一种三维卷积算子,实现对地震数据体X、Y、Z三个方向求取变化梯度,并按照一定的比例因子求取均方根;
(2)设计一种三维数据体局部梯度增强算法,实现对步骤(1)中得到的梯度数据体进行局部梯度较弱的细节增强;
(3)设计一种锐化异常边缘效果的算法,实现对步骤(2)中得到数据体局部异常边缘精细化;
(4)选取合适截止值对数据体二值化,得到小尺度岩性边缘数据体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010894942.4A CN111983680A (zh) | 2020-08-31 | 2020-08-31 | 基于嵌入三维卷积算子的小尺度岩性边界刻画方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010894942.4A CN111983680A (zh) | 2020-08-31 | 2020-08-31 | 基于嵌入三维卷积算子的小尺度岩性边界刻画方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111983680A true CN111983680A (zh) | 2020-11-24 |
Family
ID=73441509
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010894942.4A Pending CN111983680A (zh) | 2020-08-31 | 2020-08-31 | 基于嵌入三维卷积算子的小尺度岩性边界刻画方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111983680A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2008241460A1 (en) * | 2007-04-17 | 2008-10-30 | Saudi Arabian Oil Company | Enhanced isotropic 2D and 3D gradient method |
US20120257477A1 (en) * | 2011-04-06 | 2012-10-11 | Ahmed Adnan Aqrawi | Amplitude contrast seismic attribute |
CN104020492A (zh) * | 2013-07-01 | 2014-09-03 | 西安交通大学 | 一种三维地震资料的保边滤波方法 |
CN104459769A (zh) * | 2013-09-22 | 2015-03-25 | 中国石油化工股份有限公司 | 一种地震图像增强方法 |
CN105158795A (zh) * | 2015-08-27 | 2015-12-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用地层叠前纹理属性值来检测缝洞的方法 |
US20190235114A1 (en) * | 2015-05-08 | 2019-08-01 | King Fahd University Of Petroleum And Minerals | Method for processing marine seismic data |
-
2020
- 2020-08-31 CN CN202010894942.4A patent/CN111983680A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2008241460A1 (en) * | 2007-04-17 | 2008-10-30 | Saudi Arabian Oil Company | Enhanced isotropic 2D and 3D gradient method |
US20120257477A1 (en) * | 2011-04-06 | 2012-10-11 | Ahmed Adnan Aqrawi | Amplitude contrast seismic attribute |
CN104020492A (zh) * | 2013-07-01 | 2014-09-03 | 西安交通大学 | 一种三维地震资料的保边滤波方法 |
CN104459769A (zh) * | 2013-09-22 | 2015-03-25 | 中国石油化工股份有限公司 | 一种地震图像增强方法 |
US20190235114A1 (en) * | 2015-05-08 | 2019-08-01 | King Fahd University Of Petroleum And Minerals | Method for processing marine seismic data |
CN105158795A (zh) * | 2015-08-27 | 2015-12-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用地层叠前纹理属性值来检测缝洞的方法 |
Non-Patent Citations (6)
Title |
---|
ASJAD AMIN ET AL.: "A new approach for salt dome detection using a 3D multidirectional edge detector", 《APPLIED GEOPHYSICS》 * |
HUAJING CHEN ET AL.: "Application of seismic multi-attributes in the prediction of the fractures of granite buried hill", 《2016 SEG INTERNATIONAL EXPOSITION AND 86TH ANNUAL MEETING》 * |
SALEH AL-DOSSARY ET AL.: "Fault Detection and Characterization Using a 3D Multidirectional Sobel Filter", 《SPE SAUDI ARABIA SECTION ANNUAL TECHNICAL SYMPOSIUM AND EXHIBITION HELD IN KHOBAR 2013》 * |
孙夕平等: "小尺度边缘特征地震检测技术研究", 《石油地球物理勘探》 * |
沈德海等: "一种基于Sobel 算子梯度增强的边缘检测算法", 《电子设计工程》 * |
马佳国等: "精细识别小尺度废弃河道对海上低油柱油藏水平井实施的意义", 《复杂油气藏》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104047597B (zh) | 含油气泥页岩地层测井曲线标准化方法 | |
CN103838936B (zh) | 适用于浊积砂低渗透储层的高精度构造应力场模拟方法 | |
CN103513286B (zh) | 一种地质模型约束下的滩坝结构单元判别方法 | |
CN106569272B (zh) | 一种基于数据性质空间升维的地震属性融合方法 | |
CN108121009B (zh) | 开发后期复杂断块变速构造成图方法 | |
CN103775075A (zh) | 一种全井段岩性识别方法 | |
CN107966732A (zh) | 基于空间结构导向的地震属性变化率求取方法 | |
CN104863574A (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
CN105447762A (zh) | 一种流体替换的低渗透油藏水淹信息的计算方法 | |
CN110320574B (zh) | 基于缓坡三角洲薄层砂体刻画的方法 | |
CN108008117A (zh) | 裂缝发育模式与储层基质-裂缝渗透率敏感性预测方法 | |
CN114609675A (zh) | 基于高频旋回对碳酸盐岩地层沉积微地貌的定量恢复方法 | |
CN109324171A (zh) | 一种基于岩性统计的沉积相定量识别方法 | |
Fang et al. | Fine description of structure and sedimentary microfacies of Li32 block of Lijin oilfield, Dongying depression, China | |
Li et al. | Application of the Wavelet Transform and INPEFA in Sequence Stratigraphy | |
Tao et al. | A de-flat roundness method for particle shape quantitative characterization | |
CN113361111A (zh) | 一种低渗透油藏储层特征模型的划分方法 | |
CN111983680A (zh) | 基于嵌入三维卷积算子的小尺度岩性边界刻画方法 | |
CN111768503A (zh) | 一种基于三维地质模型的海砂资源量估算方法 | |
CN101625420B (zh) | 一种油藏描述的方法 | |
CN103984015A (zh) | 基于正演系数滤波的位场分离方法和位场分离系统 | |
CN113552623A (zh) | 曲流河沉积砂体边界的确定方法和确定装置 | |
CN113325474B (zh) | 生物礁判别方法 | |
Li et al. | A case study on statistical wireline log parameters in identifying shallow-water delta microfacies of Late Dongying Formation, northern Liaozhong Depression, Bohai Bay Basin | |
CN114594518B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20201124 |