CN103376469B - 一种基于超声成像测井的裂缝定量评价方法 - Google Patents
一种基于超声成像测井的裂缝定量评价方法 Download PDFInfo
- Publication number
- CN103376469B CN103376469B CN201210125518.9A CN201210125518A CN103376469B CN 103376469 B CN103376469 B CN 103376469B CN 201210125518 A CN201210125518 A CN 201210125518A CN 103376469 B CN103376469 B CN 103376469B
- Authority
- CN
- China
- Prior art keywords
- crack
- wave amplitude
- sonic wave
- sliding window
- width
- 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
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
根据本发明的一方面,提供了一种用于超声成像测井中的裂缝定量评价方法,包括步骤:a)利用超声成像测井仪来获得井周的声波振幅图像,在所述声波振幅图像中包括至少一条裂缝;b)选择所述至少一条裂缝中的一条裂缝作为评价对象;c)确定选择的裂缝在井周上的轨迹长度;d)针对所述选择的裂缝,确定搜索范围;e)设定声波振幅阈值以便在所述搜索范围内统计声波振幅低于所述声波振幅阈值的声波振幅异常点的数目;f)估计与所述数目的声波振幅异常点对应的异常面积;以及g)基于所述轨迹长度和所述异常面积,评价所述选择的裂缝的平均视宽度。
Description
技术领域
本发明大体上涉及石油地质勘探与开发的领域,尤其是涉及一种基于超声成像测井的裂缝定量评价方法。
背景技术
测井,也被称为地球物理测井或石油测井,是利用岩层的电化学特性、导电特性、声学特性、放射性等地球物理特性来测量地球物理参数的方法。在油田勘探与开发过程中,测井是确定和评价油层、气层的重要手段之一,也是解决一系列地质问题的重要手段。它能直接为石油地质和工程技术人员提供各项资料和数据。
随着油气勘探程度的提高,新发现油气藏在规模上趋于小型化,在储层物性及构造形态上趋于复杂化,因此勘探难度越来越大。常规测井在储层描述、油气层解释等方面已不能满足要求。因此在这一背景下,成像测井应运而生。
成像测井在复杂油气储层的解释评价中发挥着越来越重要的作用。成像测井具有可视性和直观性的特点,可以进行储层结构特征、构造特征、沉积特征、地应力方向等方面的评价分析,可以解决常规测井难以解决的测井地质解释问题。尤其在裂缝性储层评价中,它能直观地提供裂缝的形态、裂缝发育程度等地质信息。
成像测井可以包括超声成像测井和电成像测井。裂缝定量评价是成像测井解释的一个重要方面。裂缝定量评价包括对各种裂缝参数(包括但不限于裂缝平均视宽度)的评价。目前国内应用的超声成像测井仪例如包括斯伦贝谢公司的UBI、阿特拉斯公司的CBIL等等。斯伦贝谢公司利用有限元法对裂缝进行定量评价,适用于电成像FMS、FMI测井资料的处理而不适用于超声成像测井资料的处理。
因此,需要一种基于超声成像测井的裂缝定量评价方法。
发明内容
本发明的目的是提供一种基于超声成像测井的裂缝定量评价方法,该方法能够相对准确地确定各种裂缝参数,包括裂缝平均视宽度、裂缝平均水动力宽度及裂缝视面积孔隙度等等。
根据本发明的一方面,提供了一种基于超声成像测井的裂缝定量评价方法,包括步骤:a)利用超声成像测井仪来获得井周的声波振幅图像,在所述声波振幅图像中包括至少一条裂缝;b)选择所述至少一条裂缝中的一条裂缝作为评价对象;c)确定选择的裂缝在井周上的轨迹长度;d)针对所述选择的裂缝,确定搜索范围;e)设定声波振幅阈值以便在所述搜索范围内统计声波振幅低于所述声波振幅阈值的声波振幅异常点的数目;f)估计与所述数目的声波振幅异常点对应的异常面积;以及g)基于所述轨迹长度和所述异常面积,评价所述选择的裂缝的平均视宽度。
根据本发明的另一方面,所述确定所述选择的裂缝在井周上的轨迹长度的步骤包括按照下式确定所述轨迹长度:
L=π*CAL*(1+cosθ),
其中L为所述轨迹长度;CAL为井径;以及θ为所述选择的裂缝的视倾角。
根据本发明的另一方面,所述估计与所述数目的声波振幅异常点对应的异常面积的步骤包括估计每个声波振幅异常点在声波振幅图像上所占的面积与所述数目的乘积。
根据本发明的另一方面,每个声波振幅异常点在声波振幅图像上所占的面积按照下式得出:
其中S为每个声波振幅异常点在声波振幅图像上所占的面积;BS为钻头直径;m为超声成像测井仪在井周方向上旋转一周扫描的点数;以及d为超声成像测井仪在纵向方向上的采样间隔。
根据本发明的另一方面,所述评价所述选择的裂缝的平均视宽度的步骤包括将所述平均视宽度确定为所述异常面积除以所述轨迹长度。
根据本发明的另一方面,所述搜索范围是以所述选择的裂缝为中心在10mm-40mm内。
根据本发明的另一方面,所述设定声波振幅阈值的步骤包括设定声波振幅阈值以使得声波振幅异常点的数目占统计总数70%-90%。
根据本发明的另一方面,提供了一种基于超声成像测井的裂缝定量评价方法,包括步骤:a)利用超声成像测井仪来获得井周的声波振幅图像,在所述声波振幅图像中包括至少一条裂缝;b)在所述声波振幅图像上选择一个滑动窗并且针对所述滑动窗内的每条裂缝执行以下步骤:确定所述每条裂缝在井周上的轨迹长度;针对每条裂缝,确定相应的搜索范围;设定声波振幅阈值以便在所述搜索范围内统计声波振幅低于所述声波振幅阈值的声波振幅异常点的数目;估计与所述数目的声波振幅异常点对应的异常面积;以及基于所述轨迹长度和所述异常面积,评价所述每条裂缝的平均视宽度。
根据本发明的另一方面,所述裂缝定量评价方法进一步包括:基于所述每条裂缝的平均视宽度,按照下式确定在滑动窗内的裂缝平均视宽度:
其中VDA为在滑动窗内的裂缝平均视宽度,
Wi为第i条裂缝的平均视宽度,以及n为在滑动窗内裂缝的条数。
根据本发明的另一方面,裂缝定量评价方法进一步包括:基于所述每条裂缝的平均视宽度,按照下式确定在滑动窗内的裂缝平均水动力宽度(VAH):
其中VAH为在滑动窗内的裂缝平均水动力宽度,Wi为第i条裂缝的平均视宽度,以及n为在滑动窗内裂缝的条数。
根据本发明的另一方面,裂缝定量评价方法进一步包括:基于所述每条裂缝的平均视宽度,按照以下操作来确定在滑动窗内的裂缝视面积孔隙度:
A)确定在滑动窗内裂缝的面积
其中SS为在滑动窗内裂缝的面积,Wi为第i条裂缝的平均视宽度,n为滑动窗内裂缝的条数,以及L’为滑动窗内裂缝的长度,当滑动窗包含整条裂缝时为裂缝在井周上的轨迹长度;当滑动窗包含部分裂缝时按下式确定:
其中Dw为滑动窗的深度;wl为滑动窗长;Df为裂缝深度;CAL为井径;θ为裂缝的视倾角;L为裂缝在井周上的轨迹长度;
B)利用在滑动窗内裂缝的面积来确定裂缝视面积孔隙度:
其中VPA为在滑动窗内的裂缝视面积孔隙度,SS为在滑动窗内裂缝的面积,wl为滑动窗长,以及CAL为井径。
以上相当宽泛地概括了本发明的特征和技术优势,以便可以更好地理解本发明的以下详细描述。在下文中将描述本发明的附加特征和优势。本领域技术人员应当认识到,所公开的概念和具体实施例可以容易地用作用于修改或设计用于实现本发明的相同目的的其他结构的基础。本领域技术人员也应当认识到,这种等效的构造并不脱离所附权利要求中限定的本发明的精神和范围。
通过结合附图来阅读后面的具体实施方式,可以更好地理解本发明的特征和优点。
附图说明
现在将参照附图来解释本发明的实施例。应当注意,这些实施例用于图示基本原理,使得仅图示为了理解基本原理而必需的那些特征。附图未按比例。另外,相似标号在附图中通篇表示相似特征。
图1是根据本发明实施例的超声成像测井过程的示意图。
图2a是根据本发明实施例的裂缝平面与井眼相交的示意图。
图2b是根据本发明实施例的图2a按照井眼圆周方向展开的示意图。
图3是根据本发明实施例的针对裂缝的搜索范围的示意图。
图4a是根据本发明实施例的叠加有图3所示的搜索范围的声波振幅图像。
图4b是根据本发明实施例的在图4a所示的声波振幅图像中的裂缝的搜索范围内声波振幅数据的直方图。
图5是根据本发明实施例的滑动窗的示意图。
具体实施方式
以下具体描述涉及附图,这些附图通过图示方式示出了可以实施本发明的具体细节和实施例。充分具体描述这些实施例以使本领域的技术人员能够实施本发明。可以利用其它实施例并且可以进行结构、逻辑和电改变而不脱离本发明的范围。各种实施例未必互斥,因为一些实施例可以与一个或者多个其它实施例组合以形成新实施例。
在下面的详细描述中参照形成本详细描述的一部分的附图,在所述附图中通过说明的方式示出了其中可以实践本发明的具体实施例。在这方面,参照所描述的附图的取向使用了诸如“顶”、“底”、“前”、“后”、“首”、“尾”等等的方向术语。由于实施例的组件可以被定位在若干种不同的取向中,因此所述方向术语被用于说明的目的而决不是进行限制。要理解的是,在不背离本发明的范围的情况下可以利用其他实施例并且可以做出结构或逻辑的改变。因此,不要将下面的详细描述视为限制性意义,并且本发明的范围由所附权利要求书限定。
要理解的是,除非具体另行声明,否则这里所描述的各个示例性实施例的特征可以彼此组合。
如在本说明书中所采用的术语“裂缝”可以指的是储层裂缝,一般为张开缝。
如在本说明书中所采用的术语“裂缝定量评价”可以指的是对裂缝各种参数的评价,包括但不限于裂缝平均视宽度、裂缝平均水动力宽度及裂缝视面积孔隙度。
如在本说明书中所采用的术语“井周和/或井壁”可以指的是井眼的四周和/或墙壁。
如在本说明书中所采用的术语“井径”可以指的是井眼的直径。
如在本说明书中所采用的术语“方位”和“方位角”可互换地使用,除非上下另外明确指出。
如在本说明书中所采用的术语“滑动窗”可以指的是在声波振幅图像的深度方向上占据一定宽度的窗口,所述窗口的形状例如可以包括矩形等等常规的几何形状。
下面描述本发明的思想:根据超声成像测井数据和/或井周声波振幅图像,确定裂缝的几何特性;然后根据裂缝对声波振幅的影响,确定因裂缝引起的声波振幅异常点;接着通过统计裂缝引起的声波振幅异常点的数目,从而确定单条裂缝的平均视宽度;最后根据各条裂缝的平均视宽度,确定在滑动窗内的裂缝平均视宽度、裂缝平均水动力宽度及裂缝视面积孔隙度(VPA),从而达到通过超声成像测井对裂缝进行定量评价的目的。
接下来,具体参考附图来解释本发明的实施例。
参照图1,图1是根据本发明实施例的超声成像测井过程的示意图。超声成像测井过程是采用超声成像测井仪120进行的,所述超声成像测井仪120比如采用斯伦贝谢公司的UBI、阿特拉斯公司的CBIL、哈里伯顿公司的CAST和中油测井有限公司的BHTV等等。
如图1所示,所述超声成像测井过程包括:超声成像测井仪120向井周和/或井壁110发射100-1000kHz频率范围优选地是200-600kHz频率范围的超声波束130;超声波束130经过井壁110反射回来以形成超声回波140,从而再被超声成像测井仪120(具体地是内部的换能器,图1中未示出)接收;超声成像测井仪120以一定的角速度ω1围绕轴线125沿着扫描迹线150连续旋转,该角速度ω1的范围为0-100转/分钟,优选0-50转/分钟,更优选5-10转/分钟,不同超声成像测井仪可以采用不同的转数;同时超声成像测井仪120以一定的提升速度v1上提,该提升速度v1的范围为100-1000英尺/小时,优选600英尺/小时;因此超声成像测井仪120沿井周方向(由箭头指示)和纵向方向(由箭头y指示)连续记录经反射的超声回波140,从而产生相应的声波振幅数据及其图像,如图4所示。
根据本发明的实施例,基于用于钻开井眼的钻头(图1中未示出)直径、超声成像测井仪120在纵向方向y上的采样间隔(如图1中的符号d表示)以及超声成像测井仪120在井周方向上旋转一周扫描的声波振幅数据点数目(在下文中简称为点数),可以按照如下公式确定单个声波振幅数据点在声波振幅图像上所占的面积S:
其中BS为钻头直径;
m为超声成像测井仪120在井周方向上旋转一周扫描的点数;以及
d为超声成像测井仪120在纵向方向y上的采样间隔。
接下来参照图2a和图2b,图2a是根据本发明实施例的裂缝平面与井眼相交的示意图;图2b是图2a按照井眼圆周方向的展开图,示出了裂缝在井周上的形状。
如图2a所示,井眼220例如具有圆柱形的形状,其中井眼220的直径(在下文中也称为井径)由双向箭头符号CAL表示。当裂缝平面210与井眼220相交时,相交的是一个椭圆250。裂缝平面210与垂直于井眼220的轴线(在下文中也称为井轴)230的平面240的夹角称为裂缝的视倾角,如图2a中的角度θ所指示。裂缝平面210的法向量在垂直于井轴230的平面240上的投影与该平面240正北的夹角称为裂缝的视方位,如图2a中的角度α所指示。
如图2b所示,将圆柱形的井眼220按照井眼的圆周方向(即,北(N)-东(E)-南-西-北,图2b中的水平轴x所示)展开,此时所述椭圆250示为一条正弦波曲线260,即裂缝在井周280上的轨迹为正弦波曲线260。在图2b中的纵轴y表示井眼的深度,其中裂缝的深度由虚线270指示。
另外在其他实施例中,当裂缝平面210与圆柱形的井眼220垂直相交时,上述的椭圆250就退化成圆并且相应地所述正弦波曲线就退化成一条直线。在下文中主要以前一种情况(即,椭圆-正弦波曲线的裂缝)来解释本发明,但是本发明同样也适用于后一种情况(即,圆-直线的裂缝)。
根据如图2a和图2b所示的裂缝的产状数据(包括裂缝的深度、视倾角和视方位)和井径,可以按照如下公式确定裂缝260在井周上的轨迹长度L,如图2b所示:
L=π*CAL*(1+cosθ) (2)
其中CAL为井径;θ为裂缝的视倾角。
接下来参照图3,图3是根据本发明实施例的针对裂缝的搜索范围的示意图。如图3所示,井眼也按圆周方向(即,北-东-南-西-北,如图3中的水平轴x所示)展开。在图3中纵轴y表示井眼的深度而水平轴x表示0到360°之间的方位。与上文针对图2b所述的一样,在图3中裂缝在井周380上的形状也将是正弦波曲线的。因此在图3中实正弦线300也将代表裂缝,Df表示该裂缝300的深度,而上和下虚正弦线310和320所包围的区域代表该裂缝的搜索范围SL,即将在该搜索范围SL内对声波振幅数据点进行统计以得到由裂缝引起的声波振幅低于预定值的声波振幅异常点的数目,这将在下面进一步详细解释。
接下来参照图4a,图4a是根据本发明实施例的叠加有图3所示的搜索范围的声波振幅图像,其中所述声波振幅图像是通过如图1所述的超声成像测井过程而获得的。如图4a所示,纵轴y表示井眼的深度而水平轴x表示0到360°之间的方位,深度比例为1∶20。在图4a中,在裂缝400上面和下面的正弦波曲线410和420表示针对该裂缝的搜索范围的上界和下界,即要对声波振幅数据进行统计的区域的上界和下界。在图4a中可以根据各种技术从声波振幅图像识别出所述裂缝,例如神经网络、模糊识别、小波变换法等等,在此本发明将不再赘述。
参照图4b,图4b是根据本发明实施例的在图4a所示的声波振幅图像中的裂缝的搜索范围内声波振幅数据的直方图。如图4b所示,绘制了搜索范围内声波振幅数据的直方图。在图4b中纵轴是出现频率(即百分比,%)而横轴是声波振幅,曲线450是累计出现频率。作为实例,曲线450给出了70%所对应的声波振幅A,这意味着声波振幅小于A的声波振幅数据点的数目占统计总数(搜索范围内所有声波振幅数据点的数目)的70%,即由声波振幅小于A的数据点所包围的阴影面积占总阴影面积的70%。
现在结合图3和图4a及4b,说明根据本发明的实施例的统计过程。
根据本发明的实施例,由于上述的超声成像测井过程所得到的声波振幅数据是一个二维数组,因此根据该数组的行数可以确定出每一行所在的深度而根据该数组的列数及测井时所测的方位曲线可确定出每一列所在的方位。第i行第j列的声波振幅数据点AMP(i,j)的坐标位置(yi,xj)可按如下方式确定:
yi=D+(i-1)*d i=1,2,3,...n
(3)
j=1,2,3,...,m
其中D为声波振幅数据的起始深度;
d为超声成像测井仪在纵向方向上的采样间隔;
P1AZ为方位曲线;
n为声波振幅数据的总行数;以及
m为超声成像测井仪在井周方向上旋转一周扫描的点数。
然后以所选裂缝的轮廓线为中心,上下各加半个搜索范围以形成搜索区域,如图3所示。作为实例,每一条裂缝的搜索区域可按下式来确定:
i=1,2,3,...,72 (4)
若xi>360,则xi=xi-360
其中θ为裂缝的视倾角;α为裂缝的视方位;
CAL为井径;Df为裂缝的深度;以及
SL为搜索范围,通常在10mm-40mm内。
接着利用公式(3)和(4),对与落在公式(4)内的坐标对应的声波振幅进行统计以绘制直方图,如图4b所示。然后将搜索区域内的声波振幅数据与预定振幅值(即振幅阈值)进行比较,确定因裂缝引起的比预定振幅值小的声波振幅异常点的数目,从而利用上面的公式(1)来确定裂缝的异常面积。
根据本发明的实施例,一般而言裂缝引起的声波振幅异常点的数目占统计总数70%-90%,因此裂缝的异常面积可按下式确定:
A=S*TN*P (5)
其中A为裂缝的异常面积;
S为声成像单点在声波振幅图像上所占的面积,由公式(1)给出;
TN为搜索区域内声波振幅数据的数目;以及
p为声波振幅异常点的数目占统计总数的百分比,通常所述百分比在70%-90%之间。
根据本发明的实施例,例如利用上面的公式(2)和(5),可以确定单条裂缝平均视宽度W如下:
其中A为裂缝的异常面积;以及
L为裂缝在井周上的轨迹长度。
图5是根据本发明实施例的滑动窗550的示意图。如图5所示,纵轴y表示井眼的深度而水平轴x表示0到360°之间的方位,深度比例为1∶40。在图5中可以根据各种技术从声波振幅图像识别出裂缝,例如神经网络、模糊识别、小波变换法等等,在此本发明将不再赘述。在声波振幅图像580上正弦线500是根据图像所识别的裂缝而虚线方框550代表滑动窗并且wl表示滑动窗长。滑动窗的深度指的是滑动窗中心点所在的深度。
注意,在确定裂缝视面孔隙度(VPA)时,滑动窗内裂缝的长度指的是落在滑动窗内正弦线的那些部分的长度。
根据如图5所示的滑动窗550内各条裂缝的平均视宽度,确定所述滑动窗550内裂缝平均视宽度(VDA)、裂缝平均水动力宽度(VAH)及裂缝视面积孔隙度(VPA)。
作为实例,裂缝平均视宽度(VDA)可按下式确定:
其中Wi为第i条裂缝的平均视宽度;n为滑动窗内裂缝的条数。
作为实例,裂缝的裂缝平均水动力宽度(VAH)可按下式确定:
其中Wi为单条裂缝的平均视宽度;n为滑动窗内裂缝的条数。
作为实例,裂缝视面积孔隙度VPA可按如下步骤确定:
(1)确定滑动窗内裂缝的面积(SS)
其中Wi为单条裂缝的平均视宽度;n为滑动窗内裂缝的条数;
L’为滑动窗内裂缝的长度,即落在滑动窗内正弦线的那些部分的长度,如图5所示;当滑动窗包含整条裂缝时,即为裂缝在井周上的轨迹长度;当滑动窗包含部分裂缝时,按下式确定:
其中Dw为滑动窗的深度;wl为滑动窗长;Df为裂缝深度;CAL为井径;θ为裂缝的视倾角;L为裂缝在井周上的轨迹长度;
(2)利用滑动窗内裂缝的面积SS来确定裂缝视面积孔隙度(VPA):
其中wl为滑动窗长;CAL为井径。
综上所述,提供一种基于超声成像测井的裂缝定量评价方法,该方法能够相对准确地确定各种裂缝参数,包括裂缝平均视宽度、裂缝平均水动力宽度及裂缝视面积孔隙度等等。
例如,根据本发明的裂缝定量评价方法具有如下优点以及效果:
为储层结构特征、构造特征、沉积特征、地应力方向等方面的评价分析提供了基础。尤其在裂缝性储层评价中提供了裂缝的形态、裂缝发育程度等地质信息,从而为裂缝性油气藏的发现和勘探提供了技术支持。
虽然以上参考附图示出并描述了本发明的各种示例性实施例,但是本领域技术人员将清楚可以做出将实现本发明一些优点的各种改变和修改而未脱离本发明的精神和范围。因此,只要本发明的这些修改、变型和替换属于本发明的权利要求书及其等同技术方案的范围之内,则本发明也意图包含这些修改、变型和替换。此外,就在详细描述或权利要求书中使用术语“包括”、“具有”、“带有”或它们的其他变体来说,这样的术语意图以与术语“包含”类似的方式是包含性的。
应当理解的是,本发明的实施例可以被实施在分立电路、部分集成电路或完全集成电路或者编程装置中。此外,术语“示例性”仅仅意味着作为一个实例,而不是最佳的或最优的。还要认识到,这里所描绘的特征和/或元素是出于简单及易于理解的目的而以相对于彼此的特定尺寸示出的,实际的尺寸可能与这里所示出的尺寸显著不同。
本领域合理技术人员将明白,可以用执行相同功能的其它部件适当地进行替换。应当提到,即使在尚未明确提到这一点的那些情况下,参照具体图所解释的特征可以与其它图的特征组合。另外,可以在使用适当处理器指令的全软件实施方其中或者在利用硬件逻辑与软件逻辑的组合以实现相同结果的混合实施方其中实现本发明的方法。这样对发明概念的修改旨在于由所附权利要求覆盖。
此外,本申请的范围不旨在限于说明书中描述的过程、机器、制造、物质组成、装置、方法和步骤的具体实施例。如本领域普通技术人员从本发明的公开中容易认识到的,根据本发明可以利用与这里描述的对应实施例执行实质上相同的功能或实现实质上相同的结果的、当前存在或以后待开发的过程、机器、制造、物质组成、装置、方法或步骤。据此,所附权利要求旨在在其范围内包括这些过程、机器、制造、物质组成、装置、方法或步骤。
Claims (10)
1.一种基于超声成像测井的裂缝定量评价方法,包括步骤:
a)利用超声成像测井仪来获得井周的声波振幅图像,在所述声波振幅图像中包括至少一条裂缝;
b)选择所述至少一条裂缝中的一条裂缝作为评价对象;
c)确定选择的裂缝在井周上的轨迹长度;
d)针对所述选择的裂缝,确定搜索范围;
e)设定声波振幅阈值以便在所述搜索范围内统计声波振幅低于所述声波振幅阈值的声波振幅异常点的数目;
f)估计与所述数目的声波振幅异常点对应的异常面积;以及
g)基于所述轨迹长度和所述异常面积,评价所述选择的裂缝的平均视宽度;其中所述确定所述选择的裂缝在井周上的轨迹长度的步骤包括按照下式确定所述轨迹长度:
其中L为所述轨迹长度;
CAL为井径;以及
θ为所述选择的裂缝的视倾角。
2.根据权利要求1的裂缝定量评价方法,其中所述估计与所述数目的声波振幅异常点对应的异常面积的步骤包括估计每个声波振幅异常点在声波振幅图像上所占的面积与所述数目的乘积。
3.根据权利要求2的裂缝定量评价方法,其中每个声波振幅异常点在声波振幅图像上所占的面积按照下式得出:
其中S为每个声波振幅异常点在声波振幅图像上所占的面积;
BS为钻头直径;
m为超声成像测井仪在井周方向上旋转一周扫描的点数;以及
d为超声成像测井仪在纵向方向上的采样间隔,其中所述纵向方向与所述井周方向垂直。
4.根据权利要求1的裂缝定量评价方法,其中所述评价所述选择的裂缝的平均视宽度的步骤包括将所述平均视宽度确定为所述异常面积除以所述轨迹长度。
5.根据权利要求1的裂缝定量评价方法,其中确定搜索范围的步骤包括以所述选择的裂缝的轮廓线为中心在纵向方向上上下平移半个搜索范围,其中所述搜索范围在10mm-40mm内。
6.根据权利要求1的裂缝定量评价方法,其中所述设定声波振幅阈值的步骤包括设定声波振幅阈值以使得声波振幅异常点的数目占统计总数70%-90%。
7. 一种基于超声成像测井的裂缝定量评价方法,包括步骤:
利用超声成像测井仪来获得井周的声波振幅图像,在所述声波振幅图像中包括至少一条裂缝;和
在所述声波振幅图像上选择一个滑动窗并且针对所述滑动窗内的每条裂缝执行以下步骤:
确定所述每条裂缝在井周上的轨迹长度;
针对每条裂缝,确定相应的搜索范围;
设定声波振幅阈值以便在所述搜索范围内统计声波振幅低于所述声波振幅阈值的声波振幅异常点的数目;
估计与所述数目的声波振幅异常点对应的异常面积;以及
基于所述轨迹长度和所述异常面积,评价所述每条裂缝的平均视宽度;
其中所述确定所述每条裂缝在井周上的轨迹长度的步骤包括按照下式确定所述轨迹长度:
其中L为所述轨迹长度;
CAL为井径;以及
θ为所述每条裂缝的视倾角。
8.根据权利要求7的裂缝定量评价方法,进一步包括:基于所述每条裂缝的平均视宽度,按照下式确定在滑动窗内的裂缝平均视宽度:
其中VDA为在滑动窗内的裂缝平均视宽度,
Wi为第i条裂缝的平均视宽度,以及
n为在滑动窗内裂缝的条数。
9.根据权利要求7的裂缝定量评价方法,进一步包括:基于所述每条裂缝的平均视宽度,按照下式确定在滑动窗内的裂缝平均水动力宽度VAH:
其中VAH为在滑动窗内的裂缝平均水动力宽度,
Wi为第i条裂缝的平均视宽度,以及
n为在滑动窗内裂缝的条数。
10. 根据权利要求7的裂缝定量评价方法,进一步包括:基于所述每条裂缝的平均视宽度,按照以下操作来确定在滑动窗内的裂缝视面积孔隙度:
A)确定在滑动窗内裂缝的面积
其中SS为在滑动窗内裂缝的面积,
Wi为第i条裂缝的平均视宽度,
n为滑动窗内裂缝的条数,以及
L’为滑动窗内裂缝的长度,当滑动窗包含整条裂缝时为裂缝在井周上的轨迹长度;当滑动窗包含部分裂缝时按下式确定:
其中Dw为滑动窗的深度;wl为滑动窗长;Df为裂缝深度;CAL为井径;为裂缝的视倾角;L为裂缝在井周上的轨迹长度;
B)利用在滑动窗内裂缝的面积来确定裂缝视面积孔隙度:
其中VPA为在滑动窗内的裂缝视面积孔隙度,
SS为在滑动窗内裂缝的面积,
wl为滑动窗长,以及
CAL为井径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210125518.9A CN103376469B (zh) | 2012-04-26 | 2012-04-26 | 一种基于超声成像测井的裂缝定量评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210125518.9A CN103376469B (zh) | 2012-04-26 | 2012-04-26 | 一种基于超声成像测井的裂缝定量评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103376469A CN103376469A (zh) | 2013-10-30 |
CN103376469B true CN103376469B (zh) | 2017-09-26 |
Family
ID=49461864
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210125518.9A Active CN103376469B (zh) | 2012-04-26 | 2012-04-26 | 一种基于超声成像测井的裂缝定量评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103376469B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103791866A (zh) * | 2014-01-20 | 2014-05-14 | 河北建设勘察研究院有限公司 | 一种大直径竖井孔形检测方法 |
CN105422089B (zh) * | 2015-11-30 | 2018-07-27 | 长江大学 | 一种利用井周超声图像生成孔隙度谱的方法 |
CN105842751B (zh) * | 2016-03-17 | 2019-02-22 | 成都创源油气技术开发有限公司 | 页岩储层裂缝评价方法 |
CN106443782B (zh) * | 2016-09-20 | 2018-10-09 | 中国地质大学(北京) | 一种断层和裂缝发育密度、均匀性以及组合样式评价方法 |
CN108362780B (zh) * | 2017-12-07 | 2020-09-22 | 中航复合材料有限责任公司 | 一种评价复合材料结构r区缺陷超声检出概率的方法 |
CN108362252A (zh) * | 2018-01-19 | 2018-08-03 | 上海理工大学 | 窗户缝隙高度的测试方法 |
CN113075748B (zh) * | 2021-03-10 | 2024-04-16 | 中石化石油工程技术服务有限公司 | 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 |
CN113820317B (zh) * | 2021-08-25 | 2024-05-03 | 绍兴文理学院 | 一种基于水动力学的裂缝检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101929973A (zh) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | 裂缝储层含油气饱和度定量计算方法 |
US8147514B2 (en) * | 1998-02-06 | 2012-04-03 | P Tech, Llc. | Apparatus and method for securing a portion of a body |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2324810C2 (ru) * | 2006-05-31 | 2008-05-20 | Шлюмберже Текнолоджи Б.В. | Способ определения размеров трещины гидроразрыва пласта |
-
2012
- 2012-04-26 CN CN201210125518.9A patent/CN103376469B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8147514B2 (en) * | 1998-02-06 | 2012-04-03 | P Tech, Llc. | Apparatus and method for securing a portion of a body |
CN101929973A (zh) * | 2009-06-22 | 2010-12-29 | 中国石油天然气股份有限公司 | 裂缝储层含油气饱和度定量计算方法 |
Non-Patent Citations (5)
Title |
---|
Fracture apertures from electrical borehole scans;S.M.Luthi 等;《Geophysics》;19900731;第55卷(第7期);第821-833页 * |
利用模拟井研究用声波成像资料计算裂缝宽度问题;乔德新 等;《石油勘探与开发》;20050228;第32卷(第1期);第76-79页 * |
基岩油气藏裂缝性储层的成像测井评价及裂缝预测;牛虎林 等;《石油学报》;20100331;第31卷(第2期);第264-269页 * |
声、电成像测井资料裂缝识别技术及其应用;王鹏 等;《测井技术》;20001231;第24卷;第487-490页 * |
成像测井资料定量计算方法研究及软件开发;乔德新;《中国优秀博硕士学位论文全文数据库(博士)·基础科学辑》;20051015(第6期);第A011-75页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103376469A (zh) | 2013-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103376469B (zh) | 一种基于超声成像测井的裂缝定量评价方法 | |
CN104635269B (zh) | 一种基于叠前方位角道集预测火成岩裂缝型储层的方法 | |
CN104977617B (zh) | 储层裂缝识别方法和成像测井储层裂缝识别方法 | |
US11474271B2 (en) | Methods and systems for automated sonic imaging | |
US11604301B2 (en) | Methods and systems for automated sonic imaging | |
CN103149589A (zh) | 火成岩油气勘探方法及装置 | |
CN104520733A (zh) | 地震正交分解属性 | |
US9529109B2 (en) | Methods and apparatus for determining slowness of wavefronts | |
CN104459768A (zh) | 一种基于可视化的三维空间目标地质体追踪方法 | |
CN109425899A (zh) | 一种碳酸盐岩断层破碎带分布的预测方法及装置 | |
CN107728214A (zh) | 一种裂缝预测方法 | |
CN109655908A (zh) | 基于HTI介质的方位AVAz的正演模拟方法及系统 | |
CN105425299A (zh) | 确定地层裂缝分布的方法和装置 | |
US11719836B1 (en) | Methods of oil and gas exploration using digital imaging | |
CN103399345B (zh) | 一种潜山裂缝分布的勘测方法与装置 | |
CN108957527A (zh) | 岩层网状裂缝的地震预测方法 | |
CN107976717B (zh) | 碎屑岩油藏等级的确定方法及装置 | |
US20150109887A1 (en) | Sonic adaptor for converting sonic or ultrasonic waveform data for use with a seismic-based computer program | |
CN110531413A (zh) | 一种小断层超前三维可视化建模方法 | |
Rush et al. | Geostatistical facies modeling trends for oolitic tidal sand shoals | |
CN114017014B (zh) | 声波测井波动声学方法 | |
Verma et al. | Calibration of attribute anomalies through prestack seismic modeling | |
CN110579798B (zh) | 一种等反射角度间隔的地震采集观测方法及系统 | |
Seller et al. | Field performance of a new borehole televiewer tool and associated image processing techniques | |
Rocco et al. | Steam Development Areal Surveillance Programme in Petroleum Development Oman |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201106 Address after: No. 6, liupukang, Xicheng District, Beijing 100120 Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: CHINA NATIONAL LOGGING Corp. Patentee after: China National Petroleum Corporation Great Wall Drilling Engineering Co.,Ltd. Address before: 100101 Beijing Chaoyang District Anli Road No. 101 Patentee before: China National Petroleum Corporation Great Wall Drilling Engineering Co.,Ltd. |
|
TR01 | Transfer of patent right |