CN112255673B - 一种基于地震反演的砂体顶界面自动追踪方法 - Google Patents
一种基于地震反演的砂体顶界面自动追踪方法 Download PDFInfo
- Publication number
- CN112255673B CN112255673B CN202011030754.3A CN202011030754A CN112255673B CN 112255673 B CN112255673 B CN 112255673B CN 202011030754 A CN202011030754 A CN 202011030754A CN 112255673 B CN112255673 B CN 112255673B
- Authority
- CN
- China
- Prior art keywords
- sand body
- sand
- top interface
- bodies
- target
- 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
- 239000004576 sand Substances 0.000 title claims abstract description 264
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000011161 development Methods 0.000 claims abstract description 12
- 238000003892 spreading Methods 0.000 claims abstract description 11
- 238000005260 corrosion Methods 0.000 claims description 39
- 230000007797 corrosion Effects 0.000 claims description 39
- 238000003708 edge detection Methods 0.000 claims description 24
- 238000001914 filtration Methods 0.000 claims description 6
- 241001270131 Agaricus moelleri Species 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 229910003460 diamond Inorganic materials 0.000 claims description 3
- 239000010432 diamond Substances 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 1
- 238000003384 imaging method Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 2
- 238000010276 construction Methods 0.000 abstract 1
- 238000004891 communication Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
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/282—Application of seismic models, synthetic seismograms
-
- 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
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/643—Horizon tracking
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及到一种基于地震反演的砂体顶界面自动追踪方法,解决了现有常规微幅构造识别没有落实的储层砂体顶界面且缺少高效的自动追踪手段的问题。包括以下步骤:通过三维地震反演,获得砂体在三维空间的展布模型;基于砂体的展布情况,进行砂体顶界面的自动追踪;基于砂体顶界面自动追踪方法进行砂体顶界面自动追踪;通过自动追踪的顶界面,得到顶界面三维展布特征;基于顶界面三维空间的展布进行微幅构造识别;进而识别出相对构造高点,从而为剩余油的发现和开发提供依据。该方法针对目的层储层砂体发育特征,建立了基于数字图像处理的砂体顶界面自动追踪技术,用它来反应砂体的微幅构造,用于提高微幅构造识别的精度和效率。
Description
技术领域
本发明涉及油田储层描述领域,特别涉及一种基于地震反演的砂体顶界面自动追踪方法。
背景技术
微幅构造(指的是受构造运动和沉积作用,地层局部产生的略高于临近区域的微小幅度的起伏,一般3~5米)发育是控制开发后期油田剩余油分布的重要因素之一,砂体内的石油受浮力的作用,容易优先运聚到微幅构造发育区域内形成油藏,因此精准预测和描述微幅构造的分布,有利于寻找剩余油,能够带来较大的经济效益。
目前相关的、公开发表的文章和发明专利,主要针对地层界面的起伏进行研究,而没有真正落实到储层砂体顶界面微幅度构造的研究。主要原因是缺少高效的自动追踪手段,手动解释误差和工作量较大。为此,在砂体反演预测模型的基础上,开展砂体顶界面自动追踪技术研发,以此提高砂体顶界面微幅构造的识别精度与效率,实现砂体微幅构造高效、精准表征。
发明内容
本发明在于克服背景技术中存在的现有常规微幅构造识别没有落实的储层砂体顶界面且缺少高效的自动追踪手段的问题,而提供一种基于地震反演的砂体顶界面自动追踪方法。该基于地震反演的砂体顶界面自动追踪方法,针对目的层砂体分布特征,建立了基于数字图像处理的砂体顶界面自动追踪技术,得到砂体顶界面微幅构造特征,适用于任何砂体储层条件下的微幅构造识别;大大提高了砂体顶界面微幅构造的识别精度与效率,实现砂体微幅构造高效、精准表征。
本发明解决其问题可通过如下技术方案来达到:一种基于地震反演的砂体顶界面自动追踪方法,包括以下步骤:
1)通过三维地震反演,获得砂体在三维空间的展布模型;
2)基于砂体的展布模型,进行砂体顶界面的自动追踪:
所述自动追踪过程主要包括以下步骤:
a、对整个砂体剖面进行腐蚀运算,分离少量连通的砂体;
b、上述将少量连通的砂体分离开后,基于现有砂体的连通关系,识别不同特征砂体;
c、根据步骤b识别出的不同特征砂体,去除离散的不连通的小砂体,保留连通的厚砂体,即目标砂体;
d、对保留的连通目标砂体进行边缘检测,利用边缘检测算法识别步骤c中提取的厚砂体的边界;
e、基于砂体边缘检测结果及识别出的砂体边界,拾取连通砂体顶界面;
f、为检验拾取连通砂体顶界面的准确性,消除腐蚀运算对界面位置的影响后,砂体的顶界面在原始剖面中进行叠加显示;
3)通过步骤(2f)得到的砂体顶界面,进行三维成图,得到顶界面三维展布特征;
4)通过所成的顶界面三维面图,进行微幅构造的识别,进而识别出相对构造高点,从而为剩余油的发现和开发提供依据。
进一步的,所述步骤(2a)对整个砂体剖面进行腐蚀运算是基于图像形态学进行腐蚀运算,其具体步骤为:腐蚀运算的定义如下:
(AθB)(x,y)=min{A(x+x',y+y')-B(x',y')|(x',y')∈DB}
其中,A(x,y)代表砂体,B(x’,y’)代表结构砂体,是面积较小的砂体,B砂体的形状可以是线形、矩形、正方形、菱形、椭圆形、圆形、球形等形状,x和x’表示横坐标,y和y’表示纵坐标,DB表示x’和y’的值域;腐蚀处理的结果是使原来的砂体减小一圈,且保证了少量连通的砂体会分离开。
进一步的,所述步骤(2b)提取特征砂体的方法为,将不连通的砂体确定为不同特征砂体,并用不同灰度表示,即不同灰度代表不同特征砂体。
进一步的,所述步骤(2c)识别厚砂体的方法为,考虑砂体的连通关系及面积大小,通过比较不同特征砂体面积大小,保留面积最大的联通砂体,去除其他的面积较小的零散砂体,保留连通的目标砂体,为下一步砂体顶界面的自动追踪做准备。
进一步的,所述步骤(2d)砂体边缘检测可以分为以下2个步骤:
(1)使用高斯滤波器,以平滑砂体剖面,滤除噪声:
为了平滑砂体,使用高斯滤波器与砂体进行卷积,以减少边缘检测时明显的噪声影响;
高斯滤波器计算公式如下:
其中,x代表横坐标,y代表纵坐标,π≈3.1415926,e表示以底数为欧拉数e的指数函数,σ表示高斯标准差,标准差公式为:公式中x1,x2,x3,……xN(N为砂体样品个数,i表示第几个数,取1~N)皆为实数,其平均数为μ;
(2)计算砂体的每个位置的梯度强度和方向:
对于滤除噪声、平滑后的砂体进行计算,计算砂体的每个位置的梯度强度和方向;
边缘检测的算子返回水平Gx和垂直Gy方向的一阶导数值,由此便可以确定砂体每个位置的梯度G和方向θ;
θ=arctan(Gy/Cx)
其中G为梯度强度,θ表示梯度方向,arctan为反正切函数;
应用双阈值检测来确定真实的和潜在的边缘,通过抑制孤立的弱边缘最终完成边缘检测。
进一步的,所述步骤(2f)消除腐蚀运算对界面位置的影响的方法为:
a、计算腐蚀运算前目标砂体与腐蚀运算后目标砂体的厚度差,即得出腐蚀运算使目标砂体减小的厚度;
b、将此厚度差校正到基于腐蚀运算后拾取的目标砂体顶界面上,得到实际位置的目标砂体顶界面;
c、将校正后的目标砂体顶界面在腐蚀运算前的砂体剖面中显示,验证校正结果的准确性。
进一步的,所述步骤(2a)图像形态学进行腐蚀运算的方法为:
⑴用结构砂体B,扫描砂体A的每一个位置;
⑵用B砂体与其覆盖的砂体做“与”操作;
⑶如果都为1,结果砂体A的该位置为1,否则为0。
本发明与上述背景技术相比较可具有如下有益效果:
(1)较高的精度
砂体微幅构造自动追踪方法避免了手动解释的误差,砂体顶界面的自动追踪的结果与后验井解释结果符合较好,通过示范区43口后验井的统计,自动追踪的砂体顶界面与井点砂体顶界面误差控制在0.1%以内。
(2)较高追踪效率
砂体微幅构造自动追踪方法避免了人工解释耗时耗力的问题,可以在短时间内完成砂体微幅构造的自动追踪,例如,解释870m*900m的三维工区的一个砂体顶界面,人工拾取大约需要2.5小时,而同样的工区,本方法识别出相同砂体顶界面用时约为30分钟,效率提高了5倍。在油田开发后期剩余油描述中具有较大的应用价值。
附图说明:
附图1是砂体微幅构造自动追踪方法流程图;
附图2是本发明实施例北东一区块反演的砂体剖面图;
附图3是本发明结构砂体B对砂体A腐蚀运算的示意图;
附图4是本发明实施例腐蚀运算后砂体剖面图;
附图5是本发明实施例基于砂体的连通关系确定不同砂体特征剖面图;
附图6是本发明实施例拾取的目标连通砂体剖面图;
附图7是本发明实施例砂体边缘检测识别出的砂体边界剖面;
附图8是本发明实施例拾取连通砂体顶界面剖面图;
附图9是本发明实施例拾取的连通砂体顶界面在原始砂体剖面叠加显示图;
附图10是本发明实施例拾取的连通砂体顶界面的三维显示。
具体实施方式
本发明提出了一种基于地震反演的砂体顶界面自动追踪方法,以下就本发明步骤及所能产生的效果,配合附图、实施例详细说明如下。
实施例
以大庆油田典型北东一区块为例,依据地震反演储层砂体三维模型(地震反演是石油地质行业常用的技术,可以用来建立砂体分布模型),发明砂体顶界面的自动追踪技术,实现精细描述砂体顶界面分布,从而进行微幅构造的表征,进而确定剩余油的挖潜有利部位。如图1所示,基于地震反演的砂体顶界面自动追踪方法,包括以下步骤:
1、通过三维地震反演,获得砂体在三维空间的展布模型;
2、基于获取的砂体展布模型,进行砂体顶界面的自动追踪:
S1、通过获取的砂体在三维空间的展布模型,读取反演的展布模型砂岩数据,本发明实施例中,三维砂体模型的大小为870m*900m*3000ms,在剖面上显示砂体,如图2所示,较厚且连续的砂体为目标砂体(如图2中S砂体),展布模型上存在厚度较薄且不连续的砂层为非目标砂体,对于上下虽有少部分连通的两个砂体,但是在开发上,考虑到驱油效果,这种连通没有实际意义,但在程序识别中两个砂体是连通的,因此就需要去除上部砂体对自动追踪界面的影响,如图2中箭头所在位置处。
S2、分离少量连通的砂体:为了分离少量连通的砂体,对整个砂体剖面进行腐蚀运算:
基于图像形态学进行腐蚀运算,腐蚀运算的定义如下:
(AθB)(x,y)=min{A(x+x',y+y')-B(x',y')|(x',y')∈DB}
其中,A(x,y)代表砂体,B(x’,y’)代表结构砂体,是面积较小的砂体,B砂体的形状可以是线形、矩形、正方形、菱形、椭圆形、圆形、球形等形状,本发明实施例中,定义的结构砂体B是长度为3,角度为60°的线形结构砂体。x和x’表示横坐标,y和y’表示纵坐标,DB表示x’和y’的值域。
B对A进行腐蚀的整个过程如图3所示(黑色区域代表砂体):⑴用结构砂体B,扫描砂体A的每一个位置;⑵用B砂体与其覆盖的砂体做“与”操作(“与”操作:如果第一个条件为假,则后面的不管真假,结果都为假。即两者同时为“1”,结果才为“1”,否则为“0”);⑶如果都为1,结果砂体A的该位置为1,否则为0。
腐蚀处理的结果是使原来的砂体减小一圈,且保证了少量连通的砂体会分离开。如图4所示,腐蚀运算后砂体面积缩减,保证了少量连通的砂体会分离开来,如图中箭头所示,上下砂体分离。
S3、上述将少量连通的砂体分离开后,基于现有砂体的连通关系,提取特征砂体,将不连通的砂体确定为不同特征砂体,并用不同灰度表示,即不同灰度代表不同特征砂体,如图5所示。
S4、识别厚砂体:根据步骤3提取的不同特征砂体,考虑砂体的连通关系及面积大小,本发明实施例中,最小连通砂体面积为3m2,最大面积为2774m2,其余连通砂体面积大小介于3~2774m2之间,通过比较不同特征砂体面积大小,保留面积最大的联通砂体,去除其他的面积较小的零散砂体,得到连通的目标砂体,为下一步砂体顶界面的自动追踪做准备。如图6所示。
S5、对保留的连通目标砂体进行边缘检测,利用边缘检测算法识别步骤4中提取的厚砂体的边界。砂体边缘检测可以分为以下2个步骤:
(1)使用高斯滤波器,以平滑砂体剖面,滤除噪声:
为了尽可能减少噪声对边缘检测结果的影响,必须滤除噪声以防止由噪声引起的错误检测。为了平滑砂体,使用高斯滤波器与砂体进行卷积,以减少边缘检测时明显的噪声影响。
高斯滤波器计算公式如下:
其中,x代表横坐标,y代表纵坐标,π≈3.1415926,e表示以底数为欧拉数e的指数函数,σ表示高斯标准差,标准差公式为:公式中x1,x2,x3,……xN(N为砂体样品个数,i表示第几个数,取1~N)皆为实数,其平均数为μ。
(2)计算砂体的每个位置的梯度强度和方向:
对于滤除噪声、平滑后的砂体进行计算,计算砂体的每个位置的梯度强度和方向。
边缘检测的算子返回水平Gx和垂直Gy方向的一阶导数值,由此便可以确定砂体每个位置的梯度G和方向θ。
θ=arctan(Gy/Gx)
其中G为梯度强度,θ表示梯度方向,arctan为反正切函数。
应用双阈值检测来确定真实的和潜在的边缘,通过抑制孤立的弱边缘最终完成边缘检测。砂体边缘检测算法识别的砂体边界的结果见图7。
S6、基于砂体边缘检测结果,拾取连通砂体顶界面,即拾取图7中砂体顶部的边界线。拾取的连通砂体顶界面见图8。
S7、由于拾取的砂体界面是在腐蚀运算后的砂体上获得的,现在的界面位置必然与实际界面位置有距离差,因此必须消除腐蚀运算对界面位置的影响。步骤为:
a、计算腐蚀运算前目标砂体与腐蚀运算后目标砂体的厚度差,即得出腐蚀运算使目标砂体减小的厚度,本发明实施例中,厚度差为3;
b、将此厚度差校正到基于腐蚀运算后拾取的目标砂体顶界面上,得到实际位置的目标砂体顶界面;
c、将校正后的目标砂体顶界面在腐蚀运算前的砂体剖面中显示,验证校正结果的准确性。
消除腐蚀运算对界面位置的影响后,砂体的顶界面在原始剖面中进行叠加显示,显示结果见图9。
3、将消除腐蚀运算对界面位置的影响后的连通砂体顶界面进行三维显示,制作顶界面三维面图。
4、通过所成的顶界面三维面图,进行微幅构造的识别,进而识别出相对构造高点,通过顶界面构造的起伏状况,寻找剩余油所在的相对构造高部位,从而为剩余油的发现和开发提供依据。(如图10所示,图中可见A、B、C处为相对高的微幅构造,为寻找剩余油所在的构造高部位。)
通过本实施例区块示范区43口后验井的统计,自动追踪的砂体顶界面与井点砂体顶界面误差控制在0.1%以内。解释870m*900m的三维工区的一个砂体顶界面,人工拾取大约需要2.5小时,而同样的工区,本方法识别出相同砂体顶界面用时约为30分钟,效率提高了5倍。在油田开发后期剩余油描述中具有较大的应用价值。
Claims (2)
1.一种基于地震反演的砂体顶界面自动追踪方法,包括以下步骤:
(1)通过三维地震反演,获得砂体在三维空间的展布模型;
(2)基于砂体的展布模型,进行砂体顶界面的自动追踪:
自动追踪过程包括以下步骤:
a、对整个砂体剖面进行腐蚀运算,分离少量连通的砂体;
所述对整个砂体剖面进行腐蚀运算是基于图像形态学进行腐蚀运算,其具体步骤为:腐蚀运算的定义如下:
(AθB)(x,y)=min{A(x+x’,y+y’)-B(x’,y’)|(x’,y’)∈DB}
其中,A(x,y)代表砂体,B(x’,y’)代表结构砂体,是面积较小的砂体,结构砂体B的形状为线形、矩形、正方形、菱形、椭圆形、圆形或球形形状,x和x’表示横坐标,y和y’表示纵坐标,DB表示x’和y’的值域;
腐蚀处理的结果是使原来的砂体减小一圈,且保证了少量连通的砂体会分离开;
b、将少量连通的砂体分离开后,基于现有砂体的连通关系,识别不同特征砂体;将不连通的砂体确定为不同特征砂体,并用不同灰度表示;
c、根据步骤b识别出的不同特征砂体,去除离散的不连通的小砂体,保留连通的厚砂体,即目标砂体;通过比较不同特征砂体面积大小,保留面积最大的连通砂体,去除其他的面积较小的零散砂体,得到连通的目标砂体;
d、对保留的连通的目标砂体进行边缘检测,利用边缘检测算法识别步骤c中提取的厚砂体的边界;砂体的边缘检测包括:
使用高斯滤波器,以平滑砂体剖面,滤除噪声:为了平滑砂体,使用高斯滤波器与砂体进行卷积,以减少边缘检测时明显的噪声影响;高斯滤波器计算公式如下:
,
其中,x代表横坐标,y代表纵坐标,π≈3.1415926,e表示以底数为欧拉数e的指数函数,σ表示高斯标准差,标准差公式为:
公式中x1,x2,x3,……xN,N为砂体样品个数,i表示第几个数,取1~N,皆为实数,其平均数为μ;
计算砂体的每个位置的梯度强度和方向:对于滤除噪声、平滑后的砂体进行计算,计算砂体的每个位置的梯度强度和方向;边缘检测的算子返回水平Gx和垂直Gy方向的一阶导数值,由此确定砂体每个位置的梯度G和方向θ;
θ=arctan(Gy/Gx)
其中G为梯度强度,θ表示梯度方向,arctan为反正切函数;
应用双阈值检测来确定真实的和潜在的边缘,通过抑制孤立的弱边缘最终完成边缘检测;
e、基于砂体边缘检测结果及识别出的砂体边界,拾取连通砂体顶界面;
f、为检验拾取连通砂体顶界面的准确性,消除腐蚀运算对界面位置的影响后,砂体的顶界面在原始剖面中进行叠加显示;
所述消除腐蚀运算对界面位置的影响的步骤为:计算腐蚀运算前目标砂体与腐蚀运算后目标砂体的厚度差,即得出腐蚀运算使目标砂体减小的厚度;
将此厚度差校正到基于腐蚀运算后拾取的目标砂体顶界面上,得到实际位置的目标砂体顶界面;
将校正后的目标砂体顶界面在腐蚀运算前的砂体剖面中显示,验证校正结果的准确性;
(3)通过步骤f得到的校正后的目标砂体顶界面,进行三维成图,制作顶界面三维面图;
(4)通过制作的顶界面三维面图,进行微幅构造的识别,进而识别出相对构造高点,从而为剩余油的发现和开发提供依据。
2.如权利要求1所述的一种基于地震反演的砂体顶界面自动追踪方法,其特征在于:步骤a图像形态学进行腐蚀运算的过程为:用结构砂体B,扫描砂体A的每一个位置;用结构砂体B与其覆盖的砂体做“与”操作;如果都为1,结果砂体A的该位置为1,否则为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011030754.3A CN112255673B (zh) | 2020-09-27 | 2020-09-27 | 一种基于地震反演的砂体顶界面自动追踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011030754.3A CN112255673B (zh) | 2020-09-27 | 2020-09-27 | 一种基于地震反演的砂体顶界面自动追踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112255673A CN112255673A (zh) | 2021-01-22 |
CN112255673B true CN112255673B (zh) | 2024-03-26 |
Family
ID=74233369
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011030754.3A Active CN112255673B (zh) | 2020-09-27 | 2020-09-27 | 一种基于地震反演的砂体顶界面自动追踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112255673B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101408624A (zh) * | 2007-10-08 | 2009-04-15 | 陶庆学 | 三维地震最佳时窗河道砂体储层预测评价技术 |
CN103675941A (zh) * | 2012-09-26 | 2014-03-26 | 中国石油天然气集团公司 | 一种获得复杂储层的砂体三维形态分布的方法 |
CN104636980A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 |
CN104715472A (zh) * | 2014-04-17 | 2015-06-17 | 中国石油化工股份有限公司 | 基于数学形态学的岩性体自动追踪方法 |
CN109541685A (zh) * | 2017-09-22 | 2019-03-29 | 中国石油天然气股份有限公司 | 一种河道砂体识别方法 |
CN109948291A (zh) * | 2019-04-01 | 2019-06-28 | 西南石油大学 | 一种砂体不连续界线方向自适应识别方法 |
-
2020
- 2020-09-27 CN CN202011030754.3A patent/CN112255673B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101408624A (zh) * | 2007-10-08 | 2009-04-15 | 陶庆学 | 三维地震最佳时窗河道砂体储层预测评价技术 |
CN103675941A (zh) * | 2012-09-26 | 2014-03-26 | 中国石油天然气集团公司 | 一种获得复杂储层的砂体三维形态分布的方法 |
CN104636980A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 |
CN104715472A (zh) * | 2014-04-17 | 2015-06-17 | 中国石油化工股份有限公司 | 基于数学形态学的岩性体自动追踪方法 |
CN109541685A (zh) * | 2017-09-22 | 2019-03-29 | 中国石油天然气股份有限公司 | 一种河道砂体识别方法 |
CN109948291A (zh) * | 2019-04-01 | 2019-06-28 | 西南石油大学 | 一种砂体不连续界线方向自适应识别方法 |
Non-Patent Citations (3)
Title |
---|
中国人工智能学会.《中国人工智能进展(2007)》.北京邮电大学出版社,2007,第656-657页. * |
港西油田六区河道砂体储层预测综合研究;邹拓;《地质学刊》;第39卷(第3期);第475-480页 * |
边缘检测技术砂体边界识别方法研究;许辉群 等;《石油天然气学报(江汉石油学院学报)》;第31卷(第5期);第75-77页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112255673A (zh) | 2021-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109389128B (zh) | 电成像测井图像特征自动提取方法及装置 | |
CN105865344A (zh) | 一种基于机器视觉的工件尺寸测量方法和装置 | |
CN110232682B (zh) | 一种基于图像的轨道异物探测方法 | |
CN106871901A (zh) | 一种基于地形特征匹配的水下地形匹配导航方法 | |
CN104778707A (zh) | 一种改进广义Hough变换的电解电容检测方法 | |
CN111274939A (zh) | 基于单目摄像头的道路路面坑洼破损自动提取方法 | |
CN111381275A (zh) | 一种地震数据的初至拾取方法和装置 | |
CN115953407A (zh) | 基于计算机视觉的半导体设备维修*** | |
CN114609675A (zh) | 基于高频旋回对碳酸盐岩地层沉积微地貌的定量恢复方法 | |
CN113848209B (zh) | 一种基于无人机与激光测距的堤坝裂缝检测方法 | |
CN112255673B (zh) | 一种基于地震反演的砂体顶界面自动追踪方法 | |
CN107610174A (zh) | 一种鲁棒的基于深度信息的平面检测方法及*** | |
CN111983678A (zh) | 一种快速评价深水砂体发育潜力的方法 | |
CN114265116B (zh) | 花岗岩类潜山风化壳沟脊幅度定量统计方法及装置 | |
CN112818501A (zh) | 基于动态监测数据校正碳酸盐岩油藏静态渗透率的方法 | |
CN115640546A (zh) | 一种图像与特征信息融合的岩性识别方法 | |
CN114418965A (zh) | 一种基于激光雷达技术的露头裂缝识别方法 | |
CN108564610B (zh) | 沉积体古水流方向角的计算方法和装置 | |
CN112684502A (zh) | 基于方位旅行时差异的裂缝预测方法及*** | |
CN112505752B (zh) | 基于剪切破裂共轭准则的叠后裂缝储层表征方法 | |
CN113325474B (zh) | 生物礁判别方法 | |
CN116721346B (zh) | 一种基于深度学习算法的岸线智能识别方法 | |
CN111060966B (zh) | 一种用于静校正技术的初至拾取方法 | |
CN113255155B (zh) | 一种储层不连续界限的阻渗能力等级的评价方法 | |
CN114624775B (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 |