CN108375785B - 裂缝带位置校正方法及装置 - Google Patents
裂缝带位置校正方法及装置 Download PDFInfo
- Publication number
- CN108375785B CN108375785B CN201810029983.XA CN201810029983A CN108375785B CN 108375785 B CN108375785 B CN 108375785B CN 201810029983 A CN201810029983 A CN 201810029983A CN 108375785 B CN108375785 B CN 108375785B
- Authority
- CN
- China
- Prior art keywords
- data
- depth
- stoneley wave
- crack
- position information
- 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
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明的目的是提供一种裂缝带位置校正方法及装置。所述方法包括:获得声波测井数据中的斯通利波信号;根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。通过本发明的实施方式以校正所述裂缝带位置信息,包括所述裂缝带的深度数据和宽度数据。
Description
技术领域
本发明涉及地质勘探领域,尤其涉及一种裂缝带位置校正方法及装置。
背景技术
碳酸盐岩储层是我国陆上石油后备资源主要的接替领域和储量增长的主要目标,该储层的主要特点为岩石非均质性强,基质孔隙度偏低。对于这样复杂储层,裂缝不仅是主要的渗流通道,而且还控制着溶孔、溶洞的发育,以及地层中原始流体的分布状况和泥浆侵入特性等。裂缝在地下通常是以裂缝簇的形式存在的,即在某个深度段内存在多条裂缝,形成一个裂缝带,如何利用测井资料准确识别和定量评价裂缝带发育情况对碳酸盐岩油气勘探非常重要。
现有技术中,1989年Hornby发现当井旁地层中存在裂缝时,声波测井测量的斯通利波波形中会出现“V”字形的上行反射波和下行反射波,利用上行反射波和下行反射波的交叉点确定裂缝所在位置。1998年余春昊利用中值滤波和频率-波数域变换等技术将斯通利波分解为斯通利波直达波、斯通利波上行反射波和斯通利波下行反射波,通过反射系数极大值的深度确定裂缝所在位置。但实际情况下,可能会存在多条裂缝形成的裂缝带,上述方法因裂缝之间的干扰而导致在识别裂缝带位置时适应性较差。
发明内容
本说明书实施方式的目的是提供一种裂缝带位置校正方法、装置、电子设备及存储器介质,能够识别和定量评价裂缝带发育情况。
本说明书实施方式提供一种裂缝带位置校正方法,所述方法包括:获得声波测井数据中的斯通利波信号;根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。
本说明书实施方式还提供一种裂缝带位置校正装置,所述装置包括:第一处理单元,用于获得声波测井数据中的斯通利波信号;第二处理单元,用于根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;第三处理单元,用于根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。
本说明书实施方式还提供一种电子设备,包括:数据接收设备,处理器;所述数据接收设备,用于接收声波测井数据;所述处理器,用于获得所述声波测井数据中的斯通利波信号;根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;根据所述斯通利波能量表征值在深度维度的极大值和极小值,对裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。
本说明书实施方式还提供一种计算机存储介质,所述计算机存储介质存储有计算机程序指令,在所述计算机程序指令被执行时实现:获得声波测井数据中的斯通利波信号;根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;根据所述斯通利波能量表征值在深度维度的极大值和极小值,对裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。
由以上本说明书实施方式提供的技术方案可见,本说明书实施方式在大致确定水平裂缝带发育位置的基础之上,利用斯通利波幅度获得所述裂缝带的第二位置信息,所述第二位置信息包括所述裂缝带上边界深度数据和所述裂缝带下边界深度数据,用以评价水平裂缝带起始位置及其宽度。以上做法实现了提供了一种新方法来评价裂缝带发育情况,方便后续的孔隙度等研究。
附图说明
为了更清楚地说明本说明书实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书中记载的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本说明书提供的一种裂缝带位置校正方法的流程图;
图2是本说明书提供的一个测井结果示意图;
图3是本说明书提供的一个获得斯通利波信号的流程图;
图4是本说明书提供的一个获得第二位置信息的流程图;
图5是本说明书提供的一个模拟单条裂缝的斯通利波幅度曲线示意图;
图6是本说明书提供的一个获得校正数据的流程图;
图7是本说明书提供的一个电子设备示意图。
具体实施方式
下面将结合本说明书实施方式中的附图,对本说明书实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本说明书一部分实施方式,而不是全部的实施方式。基于本说明书中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式,都应当属于本申请保护的范围。
请参阅图1。本说明书提供的一种裂缝带位置校正方法。所述裂缝带位置校正方法可以包括以下步骤。
在本实施方式中,执行所述裂缝带位置校正方法的客体可以是具有逻辑运算功能的电子设备。所述电子设备可以是服务器和客户端。所述客户端可以为台式电脑、平板电脑、笔记本电脑、智能手机、数字助理等。当然,客户端并不限于上述具有一定实体的电子设备,其还可以为运行于上述电子设备中的软体。还可以是一种通过程序开发形成的程序软件,该程序软件可以运行于上述电子设备中。所述服务器可以为具有运算和网络交互功能的电子设备;也可以为运行于该电子设备中,为数据处理和网络交互提供业务逻辑的软体。
在本实施方式中,所述裂缝可以是指岩石受成岩作用或构造作用等产生破裂,破裂两侧的岩石沿破裂面没有发生明显的相对位移,或仅有微量位移的断裂构造的裂缝组成的裂缝带,也可以是指人为施力或诱导下造成岩石破裂形成。多个裂缝组成所述裂缝带。所述裂缝带的宽度是指所述裂缝带最上条裂缝的顶部与最下条裂缝底部之间的宽度。所述裂缝带与周边地质的物理特性不同,例如电阻率、以及对声波的传递、反射等作用。所述裂缝带是油气储层中一种重要的储渗空间,对裂缝带的评价,有利于油气的勘探和开采。
步骤S10:获得声波测井数据中的斯通利波信号.
在本实施方式中,声波测井是指在一定深度下的充液井孔内布置声源并产生井孔中的多种模式波,然后利用各种模式波的声学信息来评价井旁地层性质的一种技术。所述声波测井是用来检测井旁裂缝的重要测井依据。
在本实施方式中,所述声波测井数据是指勘探过程中声波测井获取的数据,当然还包括初始数据。所述声波测井获取的数据可以是波速、时间、振幅等,所述初始数据可以是井径大小等。
在本实施方式中,所述斯通利波是指在两种不同介质的半空间体的交界面上传播的波。具体的,例如,当单极声源在井孔中激发声波时,通常会产生四种沿着井壁传播的模式波:折射纵波、折射横波、伪瑞利波以及斯通利波。其中斯通利波作为一种导波,它的传播速度略低于井孔流体的纵波速度,并且它对穿过井孔的水平裂缝非常敏感,因为此时斯通利波的能量会沿着水平裂缝泄漏到地层中去。本说明书采用斯通利波对所述裂缝带的响应,识别所述裂缝带。
在本实施方式中,从接收器中分离出的斯通利波信号包括斯通利波对应的接收器深度维度的位置数据,以及接收器中分离出的斯通利波的时间维度信息。例如,在某个时间段中,接收器深度位置数据为1114米,接收器接收的声波时间段为时间△t,则所述斯通利波在该深度的信号对应于此时的深度1114米,所述斯通利波信息分布在时间△t上。
在本实施方式中,根据声波测井数据,获得斯通利波信号,可以是利用所述斯通利波在声波测井数据中,与其他波性质不同来分离获得所述斯通利波。所述性质可以是,声波的速度、声波的幅度、声波的周期、声波的能量等。例如,在声波测井中,接收器记录声波的全波列波形图,所述全波列波形图包括滑行纵波、滑行横波、伪瑞利波、斯通利波等,根据斯通利波的速度与其他波形速度不同,接收器接收到的声波的时间信号不同,分离所述斯通利波等方法。具体的,可以是根据声波测井数据,获取单极波形数据;所述单极波形数据进行低通滤波,所述低通滤波的滤波频率范围可以是200Hz-2000Hz,获取所述斯通利波信号,有利于利用斯通利波对所述裂缝带的敏感性分析裂缝带的位置和宽度。
步骤S12:根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值。
在本实施方式中,时间数据表征斯通利波信号在深度位置处对应的时间,所述斯通利波幅度对应于该时间范围内斯通利波的平均振幅,是一种声波的能量表示。
在本实施方式中,斯通利波能量表征值用于表示接收器在不同位置处接收的斯通利波的能量。获取所述斯通利波能量表征值可以通过求取对应时间范围内所述斯通利波的平均振幅来获取所述斯通利波能量表征值。具体的,例如,在所述斯通利波信号的第一个周期和/或第二个周期对应的时间区间内,计算获得斯通利波的平均振幅值作为所述斯通利波能量的表征值对应于该接收器深度位置。请图2(a)是本说明书提供的一个单级波形示意图。参阅图2 (b),对低通滤波后的波形进行以开时间窗的形式获取信息,所述开时间窗可以是针对斯通利波的第一个周期和/或第二个周期开窗,此时干扰波信号较少,往后周期的波形中还包括其他的干扰波,时间窗长度可以是一到两个波长对应的时间,一般为1ms。在该时间窗内取振幅平均值,作为所述斯通利波能量表征值。图2中3ft表示声源T与接收器之间的距离。请参阅图2(d),根据每个位置处获得的斯通利波能量表征值,获得在深度位置处,斯通利波幅度曲线。
步骤S14:根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。
请参阅图2,在本实施方式中,所述裂缝带的初始位置信息可以是根据利用斯通利波波形中的“V”字形特征识别方法或电阻率成像技术图中的正弦曲线来确定水平裂缝带的大致位置。所述斯通利波波形中的“V”字形特征是指低频斯通利波中包含直达波和反射波信息,在低频斯通利波的变密度图上,经常可以分辨出明显的反射斯通利波,即“V”字形的干涉条纹。根据“V”字形的干涉条纹确定裂缝带的位置数据。“V”字形特征对于单条裂缝效果最佳,多条裂缝的裂缝带也适用,但是对于裂缝带效果不一定好。例如,请参阅图2(c)在出现“V”字形的干涉条纹的斯通利波变密度图,不一定对应存在裂缝带位置,存在裂缝带的地方,在斯通利波变密度图也不一定有明显的“V”字形特征。所述电阻率成像技术是指根据不同介质电阻率不同原理进行对地层分析,来寻找裂缝带的技术。与斯通利波“V”字形特征识别方法相比,利用电阻率成像图获取的水平裂缝带所在位置更为准确一些。但是,对于裂缝带宽度较小的裂缝带等情况下,所述电阻率成像技术适用性较差。本说明书将所述斯通利波信号结合位置数据,对所述裂缝带的位置数据进行校正,得到第二位置信息。
在本实施方式中,所述第二位置信息包括所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。具体的,是指所述裂缝带最上一条裂缝的顶部对应的深度数据,和所述裂缝带最下一条裂缝的底部对应的深度数据。
在本实施方式中,根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正。请参阅图2(d),其中,所述斯通利波能量表征值在深度维度的极大值和极小值分别在XXX5.2和XXX4.75处,则XXX4.75作为校正后的上边界深度数据,XXX5.2作为校正后的下边界深度数据。
上述实施方式,通过获取每个深度位置的斯通利波能量表征值,利用斯通利波对裂缝带的敏感特性,以极小值对应的深度数据作为所述裂缝带的上边界深度,极大值对应的深度数据作为所述裂缝带的下边界深度,获得较为准备的裂缝带位置信息。
在一个具体的场景示例中,执行裂缝带位置校正方法的客体可以是客户端或服务器。所述客户端或服务器根据接收到的声波测井数据进行滤波、正演等操作,获得所述裂缝带的宽度数据和深度数据。所述客户端或服务器将所述宽度数据和深度数据以文字或图像的形式展示给用户,以方便用户进行分析。
请参阅图2。在本场景示例中,所述裂缝带包括裂缝带的位置数据。所述位置数据可以是根据斯通利波波形中的“V”字形特征或电阻率成像技术来确定水平裂缝带的大致位置。图2(c)展示的是斯通利波变密度图,所述变密度图是声波全波测井的一种记录方式。它记录的是较长一段声波波形的变化,是将井下仪器接收到的声波信号中各波峰幅度的变化变换为亮度的变化后记录在胶片上。因为地层较薄,反射波并没有展示出清晰的“V”字形,但在大多数情况下可以根据“V”字形特征获得对应的裂缝带深度数据。图2(e)是电阻率成像图。电阻率成像技术包括对电阻率扫描成像测井数据进行处理,得到反应井周的电阻率井成像图像。地层中,不同介质对应不同电阻率,对不同电阻率进行标识,获得对应的不同颜色深浅电阻率成像图。在电阻率成像图中观察到了几条清晰的黑色水平条带,这说明该井段存在一个水平裂缝带。
在本场景示例中,根据声波测井资料,获得斯通利波信号中的用于表征所述裂缝带的信号。在声波测井的背景基础上,根据获得的所述声波测井的测井资料,对所述声波测井的测井资料进行处理,获取单极波形并进行低通滤波,获取较纯净的斯通利波信号,然后通过开时间窗的方式计算出一条斯通利波幅度曲线。具体的,可以是对于实际测量的单极波形采用数值滤波器或模拟滤波器的滤波方法进行低通滤波。滤波频率范围可以选择为200Hz-2000Hz,即只通过200Hz-2000Hz的波,此时能有效滤除其他的井孔模式波,比如纵波、横波以及伪瑞利波等,滤波结果请参阅图2(a)。
在本场景示例中,根据斯通利波信号中表征所述裂缝带的信号,获得对应的表征值,以识别所述裂缝带。具体的,可以是计算所述时间窗内的斯通利波平均振幅值,获得所述斯通利波能量的表征值。再根据每个表征值对应深度数据,获得在深度—幅度坐标系中的斯通利波幅度曲线。
在本场景示例中,通过所述表征值校正所述位置数据,获得所述第二位置信息,所述第二位置信息包括所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。具体的,根据斯通利波波形中的“V”字形特征或电阻率成像技术来确定水平裂缝带的大致位置,再根据所述表征值确定所述裂缝带的具体深度数据和宽度数据。在实际情况下,斯通利波波形中的“V”字形特征因干扰波存在,可能会出现在不同位置深度,将所述斯通利波幅度曲线与所述斯通利波波形中的“V”字形特征进行对比,可以精确确定具体的裂缝带位置。所述斯通利波幅度曲线中,极小值的深度数据对应所述裂缝带的上边界深度数据,极大值的深度数据对应所述裂缝带的下边界深度数据。请参阅图2(d),即所述裂缝带第二数据中,上边界深度数据位XXX4.75,下边界深度数据位XXX5.2。当然,若位置数据是由电阻率成像获得的,在获取第二位置信息时同理,此处不再复述。当然,所述位置数据还可以是根据斯通利波波形中的“V”字形特征和所述电阻率成像技术共同获取的。
在本场景示例中,根据声波测井数据进行正演,得到第一校正数据和第二校正数据,用以校正所述第二位置信息。可以是根据所述井径大小、目标地层纵波速度、目标地层横波速度建立数值模型,来获得所述第一校正数据和所述第二校正数据。例如,采用与实际声波测井的目标井大致同等井况条件(井径216mm,地层纵波速度6500m/s,横波速度3500m/s)进行模拟,根据力的平衡建立平衡方程,所述平衡方程包括内源,即声源对应的力。根据实际声波测井,设定边界条件等。根据所述平衡方程和边界条件,通过二阶有限差分法等,求解所述平衡方程,获得所述斯通利波幅度与所述深度的数据集,并映射至深度—幅度坐标系中。在一个数值模型中,所述模型中的裂缝带为1mm,此时可以提高精度和减少运算量,获得对应的斯通利波幅度曲线。其中,根据所述斯通利波信号的数据集可以获得第一校正数据和第二校正数据。具体的,在图形上的表现形式是,曲线中有多个峰值和谷值,并且水平裂缝左侧幅度最大值与所述模型中裂缝带距离h1,右侧幅度最小值与所述模型中裂缝带距离h2。所述h1就是在此场景下的第二校正数据,所述h2就是在此场景中的第一校正数据。请参阅图5,第一校正数据h2为0.12m,第二校正数据h1为0.12m。
在本场景示例中,还可以是通过数值拟合的方式校正所述第二位置信息。例如,通过数值模型进行正演,修改所述数值模型中的裂缝带宽度,获得不同裂缝带宽度对应的上边界深度校正数据和下边界深度校正数据。将所述第二位置信息与不同裂缝带宽度下进行正演的上边界深度校正数据和下边界深度校正数据,进行匹配,获得校正后第二位置信息。例如,第二位置信息中上边界1111.3m深度,下边界1111.5m深度;所述数值模型中,调整预设的裂缝带宽度,当设定裂缝带宽度为0.18时,对应的极大值与极小值的深度数据差值为0.2m,其中,极大值的深度数据与所述模型中裂缝带的底部深度数据差值为0.009,极小值的深度数据与所述模型中裂缝带的顶部深度数据差值为0.011。根据上述数据进行校正,得到所述裂缝带的上边界数据位1111.309m,所述裂缝带的下边界数据位1111.489m,所述裂缝带的宽度数据为 0.18m。当然,在本场景示例中,可以基于上述思想通过拟合公式的方式来识别所述裂缝带。具体的,可以是通过正演,获得几组不同裂缝带宽度对应的斯通利波信号,将多组数据进行公式拟合,可以是线性拟合也可以是多项式拟合,得到关联式。最后,将所述第二位置信息或相关数据带入所述关联式,获得对应的校正值,校正所述第二位置信息。
在本场景示例中,根据所述第一校正数据和第二校正数据,校正所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。具体的,例如,在所述深度—幅度坐标系中,根据正演获得的斯通利波信号,第一校正数据为0.12m,第二校正数据为0.12m。所述第二位置信息中,裂缝带下侧最近的幅度峰值对应深度为XXX5.2m,减去第二校正数据0.12m,即为所述裂缝带底部深度XXX5.08m;距离裂缝带上侧最近的幅度谷值对应深度为XXX4.75m,加上第一校正数据0.12m,即为所述裂缝带顶部深度XXX4.87m。所述裂缝带的宽度为,所述裂缝带底部深度与所述裂缝带顶部深度的差值,即为0.21m。
请参阅图3,在一个实施方式中,根据声波测井数据,获得斯通利波信号包括以下步骤。
步骤S20:根据声波测井数据,获取单极波形数据。
在本实施方式中,所述单极波形为所述数据都位于基准值的一侧。具体的,在波形图中,可以所述数据都位于是基准线上方,即单极性。
在本实施方式中,获取单级波形数据,用于分析所述斯通利波的幅度,有利于获得所述斯通利波的表征值。所述波形中,包括时间信号△t和所述声波的幅度,以及所述波形对应的该点的深度数据。
步骤S22:对所述单极波形数据进行低通滤波,获取所述斯通利波信号。
在本实施方式中,通过低通滤波,提取所述斯通利波信号。所述斯通利波信号相较于其他信号,所述斯通利波信号处于低频区,尤其处于200Hz至2000Hz之间。
在本实施方式中,获取的所述斯通利波信号包括深度维度、时间维度、幅度等。具体的,在每个接收器处,接收的声波信号对应该点的深度位置。所述的时间维度用于表征所述斯通利波在时间长度上的能量分布。所述斯通利波幅度用于表征声波的能量,根据预设的时间范围,取该时间范围的平均振幅,作为该深度点的斯通利波能量表征值。根据所述斯通利波信号的深度维度和每个深度维度的表征值,形成所述斯通利波信号数据集。根据所述斯通利波信号数据集,映射至深度—幅度坐标系,获得所述斯通利波幅度曲线,用以分析所述裂缝带的位置。
通过上述方式,可以求取较为纯净的斯通利波信号,方便后续数据分析。
在一个实施方式中,所述低通滤波的滤波频率范围为大于等于200Hz且小于等于2000Hz。
在本实施方式中,所述斯通利波作为一种导波,它的传播速度略低于井孔流体的纵波速度,与其他波对比,所述斯通利波幅度较大,对穿过井孔的水平裂缝非常敏感。所述低通滤波的滤波频率范围为大于等于200Hz且小于等于2000Hz时,即只通过200Hz-2000Hz的波,能有效滤除其他井孔模式波,比如纵波、横波以及伪瑞利波等。
通过上述方式,可以求取较为纯净的斯通利波信号,方便后续数据分析。
在一个实施方式中,根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值时包括:计算指定时间范围内的所述斯通利波信号的振幅平均值;以所述振幅平均值作为该深度位置的斯通利波能量表征值;其中,所述指定时间为所述斯通利波信号的第一个周期和/或第二个周期对应的时间区间。
在本实施方式中,通过表征值来表征对应深度点的斯通利波能量信号,以用于分析所述斯通利波对所述裂缝带的能量变化响应。具体的,所述斯通利波信号包括深度维度、时间维度、和振幅。在所述时间维度上,声波信号包括多个周期,在第一个周期和/或第二个周期对应的时间区间中,可以获得较为纯净的斯通利波信号,往后的斯通利波信号中包括了其他干扰波,不宜采用。将所述斯通利波信号的振幅平均值作为所述斯通利波能量的表征值,用于表征接收器接收得到的斯通利波能量。当然所述方法还包括根据振幅平均值的平方来作为所述斯通利波能量的表征值。
在本实施方式中,所述第一个周期与第二个周期表征声波波包对应的时间范围,第一个周期可以是所述斯通利波第一个波包对应的时间范围,第二个周期可以是所述斯通利波第二个波包对应的时间范围。
在本实施方式中,通过求取斯通利波在指定时间范围内幅度的平均值,作为所述斯通利波能量表征值,以表征在对应深度位置处,接收器接收到的斯通利波能量。
请参阅图4,在一个实施方式中,所述方法还可以包括以下步骤。
步骤S30:根据声波测井数据进行正演,得到第一校正数据和第二校正数据;所述第一校正数据用于校正所述裂缝带上边界深度数据,所述第二校正数据用于校正所述裂缝带下边界深度数据。
在本实施方式中,所述正演是指由源的属性推导出场的分布属性;具体的,可以是已知地球介质的性质,如(地震波传播速度等),获得波的走时等(即地震波在地球中的传播时间等)。例如,在本实施方式中,可以是根据所述井径大小、目标地层纵波速度、目标地层横波速度建立数值模型,来获得所述第一校正数据和所述第二校正数据。所述数值模型是对特定的地质、地球物理问题作适当的简化,采用数值计算的方法获取声波响应。
在本实施方式中,可以通过数值模拟方法进行正演。所述数值模拟可以选择有限差分方法或有限元方法。例如,根据力的平衡建立平衡方程,所述平衡方程包括内源,即声源对应的力。根据实际声波测井,设定边界条件,如,所述边界条件可以是完全吸收边界条件等。根据所述平衡方程和边界条件,通过二阶有限差分法等,求解平衡方程。
在本实施方式中,可以通过正演获得不同裂缝带宽度对应的斯通利波信号。具体的,通过数值模型进行正演,获得数值模型中,不同深度位置对应的斯通利波能量,即斯通利波幅度,建立深度和幅度的坐标系,用于表征所述模型中的斯通利波信号。具体的,可以根据不同的初始参数获得不同的信号。例如,可以通过调节裂缝带宽度,获得模型中,对应不同裂缝带宽度的斯通利波信号。在一种优选的实施方式中,所述数模模型中的裂缝带宽度选取为 1mm。
请参阅图5。在本实施方式中,在正演获得的斯通利波信号中,根据斯通利波幅度曲线,取得对应所述模型中的裂缝带的第一校正数据和第二校正数据。图5中,在所述模型中,裂缝带宽度为1mm,深度位置大约在1.2m。根据正演,获得斯通利波幅度曲线,在距所述裂缝带左侧0.12m处达到极大值,所述极大值为所述第二校正数据。根据斯通利波幅度曲线,在距所述裂缝带又侧0.12m处达到极小值,所述极小值为所述第一校正数据。
步骤S32:根据所述第一校正数据和第二校正数据,校正所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。
在本实施方式中,通过所述第一校正数据校正所述裂缝带上边界深度数据,所述第二校正数据校正所述裂缝带下边界深度数据。具体的,所述运算方式根据所述正演模型中所述斯通利波信号与裂缝带的宽度位置的关系确定。例如,根据所述第二位置数据中所述裂缝带上边界深度数据和所述裂缝带下边界深度数据,将所述裂缝带上边界深度数据加上所述第一校正数据,将所述裂缝带下边界深度数据减去所述第二校正数据,得到校正后的所述裂缝带的上边界深度数据和所述裂缝带下边界深度数据。
在本实施方式中,通过正演,进一步对所述裂缝带第二位置信息进行校正,获得较为准确的所述裂缝带位置信息。
请参阅图6。在一个实施方式中,根据声波测井数据正演包括以下步骤。
步骤S40:根据声波测井数据,建立数值模型;其中,所述声波测井数据至少包括井径大小、目标地层纵波速度、目标地层横波速度;所述数值模型中的斯通利波信号至少包括预设裂缝带宽度数据和深度数据。
步骤S42:根据所述数值模型中的斯通利波信号,获得所述第一校正数据和所述第二校正数据;其中,所述第一校正数据为所述数值模型中的斯通利波信号极小值对应的所述深度数据与所述预设裂缝带顶部对应的所述深度数据的差值;所述第二校正数据为所述数值模型中的斯通利波信号极大值对应的所述深度数据与所述预设裂缝带底部对应的所述深度数据的差值。
在本实施方式中,通过正演来校正所述裂缝带的宽度数据和深度数据。具体的,可以采用与实际声波测井的目标井大致同等井况条件,如井径216mm,地层纵波速度6500m/s,横波速度3500m/s,来进行模拟,根据力的平衡建立平衡方程,所述平衡方程包括内源,即声源对应的力。根据实际声波测井,设定边界条件,如,所述边界条件可以是完全吸收边界条件等。根据所述平衡方程和边界条件,通过二阶有限差分法等,求解所述平衡方程,获得所述斯通利波幅度与所述深度的数据集,并映射至深度—幅度坐标系中。在一个数值模型中,所述模型中的裂缝带为1mm,此时可以提高精度和减少运算量,获得对应的斯通利波幅度曲线。其中,根据所述斯通利波信号的数据集可以获得第一校正数据和第二校正数据。具体的,在图形上的表现形式是,曲线中有多个峰值和谷值,并且水平裂缝左侧幅度最大值与所述模型中裂缝带距离h1,右侧幅度最小值与所述模型中裂缝带距离h2。所述h1就是在此场景下的第二校正数据,所述h2就是在此场景中的第一校正数据。
在本实施方式中,还可以是通过数值拟合的方式校正所述第二位置信息。例如,通过数值模型进行正演,修改所述数值模型中的裂缝带宽度,获得不同裂缝带宽度对应的上边界深度校正数据和下边界深度校正数据。将所述第二位置信息与不同裂缝带宽度下进行正演的上边界深度校正数据和下边界深度校正数据,进行匹配,获得校正后第二位置信息。例如,第二位置信息中上边界XXX1.3m深度,下边界XXX1.5m深度;所述数值模型中,当设定裂缝带宽度为0.18时,对应的极大值与极小值的深度数据差值为0.2m,其中,极大值的深度数据与所述模型中裂缝带的底部深度数据差值为0.009,极小值的深度数据与所述模型中裂缝带的顶部深度数据差值为0.011。根据上述数据进行校正,得到所述裂缝带的上边界数据位 XXX1.309m,所述裂缝带的上边界数据位XXX1.489m,所述裂缝带的宽度数据为0.18m。当然,在本场景示例中,可以基于上述思想通过拟合公式的方式来识别所述裂缝带。具体的,可以是通过正演,获得几组不同裂缝带宽度对应的斯通利波信号,将多组数据进行公式拟合,可以是线性拟合也可以是多项式拟合,得到关联式。最后,将所述第二位置信息或相关数据带入所述关联式,获得对应的校正值,校正所述第二位置信息。
在一个实施方式中,所述预设裂缝带宽度数据为1mm。
在本实施方式中,所述模型中的裂缝带可以预设为1mm,此时运算时,可以忽略所述裂缝带的宽度,可以提高精度和减少运算量,获得对应的斯通利波幅度曲线。
通过本说明书实施方式,可以根据声波测井资料,利用斯通利波幅度曲线获取水平裂缝带的深度数据和宽度数据。通过本说明书实施方式还可以进一步校正所述深度数据和所述宽度数据,为后续分析提供准确依据。例如,在一个场景中,通过本说明书实施方式,获得水平裂缝带顶部深度数据为XXX4.87m,底部深度数据为XXX5.09m,由此计算的水平裂缝带宽度数据为0.21m,当通过电阻率成像,也能获得水平裂缝带宽度,其值为0.2m,两者非常接近,所述方法进一步校正了电阻率成像获得的水平裂缝带的深度数据和宽度数据。需要指出的是,实际声波测井的深度采样间隔为0.05m,若缩小深度采样间隔,应该能进一步提高所述方法的准确性。
本说明书实施方式中还提供了一种裂缝带识别装置,如下面的实施方式所述。由于一种裂缝带识别装置解决问题的原理与一种裂缝带位置校正方法相似,因此一种裂缝带识别装置的实施可以参见一种裂缝带位置校正方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。该装置具体可以包括:第一处理单元、第二处理单元、第三处理单元。下面对该结构进行具体说明。
第一处理单元,用于获得声波测井数据中的斯通利波信号。
在本实施方式中,根据声波测井数据,获得斯通利波信号,可以是利用所述斯通利波在声波测井数据中,与其他波性质不同来分离获得所述斯通利波。所述性质可以是,声波的速度、声波的幅度、声波的周期、声波的能量等。例如,在声波测井中,接收器记录声波的全波列波形图,所述全波列波形图包括滑行纵波、滑行横波、伪瑞利波、斯通利波等,更具斯通利波的速度与其他波形速度不同,接收器接收到的声波的时间信号不同,分离所述斯通利波等方法。具体的,可以是根据声波测井数据,获取单极波形数据;所述单极波形数据进行低通滤波,所述低通滤波的滤波频率范围可以是200Hz-2000Hz,获取所述斯通利波信号,有利于利用斯通利波对所述裂缝带的敏感性分析裂缝带的位置和宽度。
第二处理单元,用于根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值。
在本实施方式中,斯通利波能量表征值用于表示不同位置处斯通利波的能量。获取所述斯通利波能量表征值可以通过求取对应时间范围内所述斯通利波的平均振幅来获取所述斯通利波能量表征值。具体的,例如,在所述斯通利波信号的第一个周期和/或第二个周期对应的时间区间内,计算获得斯通利波的平均振幅值作为所述斯通利波能量的表征值对应于该接收器深度位置。请参阅图2(b),对低通滤波后的波形进行以开时间窗的形式获取信息,所述开时间窗可以是针对斯通利波的第一个周期和第二个周期开窗,此时干扰波信号较少,往后周期的波形中还包括其他的干扰波,时间窗长度可以是一到两个波长对应的时间,一般为 1ms。在该时间窗内取振幅平均值,作为所述斯通利波能量表征值。图中3ft表示声源T与接收器之间的距离。请参阅图2(d),根据每个位置出获得的斯通利波能量表征值,获得在深度位置处,斯通利波幅度曲线。
第三处理单元,用于根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。
在本实施方式中,所述裂缝带的初始位置信息可以是根据利用斯通利波波形中的“V”字形特征识别方法或电阻率成像技术图中的正弦曲线来确定水平裂缝带的大致位置。所述斯通利波波形中的“V”字形特征是指低频斯通利波中包含直达波和反射波信息,在低频斯通利波的变密度图上,经常可以分辨出明显的反射斯通利波,即“V”字形的干涉条纹。根据“V”字形的干涉条纹确定裂缝带的位置数据。“V”字形特征对于单条裂缝效果最佳,多条裂缝的裂缝带也适用,但是对于裂缝带效果不一定好。例如,请参阅图2(c)在出现“V”字形的干涉条纹的斯通利波变密度图,不一定对应存在裂缝带位置,存在裂缝带的地方,在斯通利波变密度图也不一定有明显的“V”字形特征。所述电阻率成像技术是指根据不同介质电阻率不同原理进行对地层分析,来寻找裂缝带的技术。与斯通利波“V”字形特征识别方法相比,利用电阻率成像图获取的水平裂缝带所在位置更为准确一些。但是,对于裂缝带宽度较小的裂缝带等情况下,所述电阻率成像技术适用性较差。本说明书将所述斯通利波信号结合位置数据,对所述裂缝带的位置数据进行校正,得到第二位置信息。
在本实施方式中,所述第二位置信息包括所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。具体的,是指所述裂缝带最上一条裂缝的顶部对应的深度数据,和所述裂缝带最下一条裂缝的底部对应的深度数据。
在本实施方式中,根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正。请参阅图2(d),其中,所述斯通利波能量表征值在深度维度的极大值和极小值分别在XXX5.2和XXX4.75处,则XXX4.75作为校正后的上边界深度数据,XXX5.2作为校正后的下边界深度数据。
上述实施方式,通过获取每个深度位置的斯通利波能量表征值,利用斯通利波对裂缝带的敏感特性,以极小值对应的深度数据作为所述裂缝带的上边界深度,极大值对应的深度数据作为所述裂缝带的下边界深度,获得较为准备的裂缝带位置信息。
在一个实施方式中,所述装置还可以包括以下单元。
第四处理单元,用于根据声波测井数据进行正演,得到第一校正数据和第二校正数据;所述第一校正数据用于校正所述裂缝带上边界深度数据,所述第二校正数据用于校正所述裂缝带下边界深度数据;
第五处理单元,用于根据所述第一校正数据和第二校正数据,校正所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。
在本实施方式中,所述正演是指由源的属性推导出场的分布属性;具体的,可以是已知地球介质的性质,如(地震波传播速度等),获得波的走时等(即地震波在地球中的传播时间等)。例如,在本实施方式中,可以是根据所述井径大小、目标地层纵波速度、目标地层横波速度建立数值模型,来获得所述第一校正数据和所述第二校正数据。所述数值模型是对特定的地质、地球物理问题作适当的简化,采用数值计算的方法获取声波响应。
在本实施方式中,可以通过数值模拟方法进行正演。所述数值模拟可以选择有限差分方法或有限元方法。例如,根据力的平衡建立平衡方程,所述平衡方程包括内源,即声源对应的力。根据实际声波测井,设定边界条件,如,所述边界条件可以是完全吸收边界条件等。根据所述平衡方程和边界条件,通过二阶有限差分法等,求解平衡方程。
在本实施方式中,可以通过正演获得不同裂缝带宽度对应的斯通利波信号。具体的,通过数值模型进行正演,获得数值模型中,不同深度位置对应的斯通利波能量,即斯通利波幅度,建立深度和幅度的坐标系,用于表征所述模型中的斯通利波信号。具体的,可以根据不同的初始参数获得不同的信号。例如,可以通过调节裂缝带宽度,获得模型中,对应不同裂缝带宽度的斯通利波信号。在一种优选的实施方式中,所述数模模型中的裂缝带宽度选取为 1mm。
请参阅图5。在本实施方式中,在正演获得的斯通利波信号中,根据斯通利波幅度曲线,取得对应所述模型中的裂缝带的第一校正数据和第二校正数据。图5中,在所述模型中,裂缝带宽度为1mm,深度位置大约在1.2m。根据正演,获得斯通利波幅度曲线,在距所述裂缝带左侧0.12m处达到极大值,所述极大值为所述第二校正数据。根据斯通利波幅度曲线,在距所述裂缝带又侧0.12m处达到极小值,所述极小值为所述第一校正数据。
在本实施方式中,通过所述第一校正数据校正所述裂缝带上边界深度数据,所述第二校正数据校正所述裂缝带下边界深度数据。具体的,所述运算方式根据所述正演模型中所述斯通利波信号与裂缝带的宽度位置的关系确定。例如,在图5中,根据所述第二位置数据中所述裂缝带上边界深度数据和所述裂缝带下边界深度数据,将所述裂缝带上边界深度数据加上所述第一校正数据,将所述裂缝带下边界深度数据减去所述第二校正数据,得到校正后的所述裂缝带的上边界深度数据和所述裂缝带下边界深度数据。
在本实施方式中,通过正演,进一步对所述裂缝带第二位置信息进行校正,获得较为准确的所述裂缝带位置信息。
请参阅图7。本说明书实施方式中还提供了一种电子设备,包括:数据接收设备,处理器;所述数据接收设备,用于接收声波测井数据;所述处理器,用于获得所述声波测井数据中的斯通利波信号;根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;根据所述斯通利波能量表征值在深度维度的极大值和极小值,对裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。
所述数据接收设备是用户和计算机系统之间进行信息交换的主要装置之一。所述数据接收设备包括键盘、鼠标、摄像头、扫描仪、光笔、手写输入板、语音输入装置等;所述数据接收设备用于把原始数据和处理这些数的程序输入到计算机中。
所述处理器可以按任何适当的方式实现。例如,处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。
在本实施方式中,该电子设备具体实现的功能和效果,可以与其它实施方式对照解释,在此不再赘述。
本说明书实施方式中还提供了一种计算机存储介质,所述计算机存储介质存储有计算机程序指令,在所述计算机程序指令被执行时实现:获得声波测井数据中的斯通利波信号;根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;根据所述斯通利波能量表征值在深度维度的极大值和极小值,对裂缝带的位置信息进行校正,得到第二位置信息;其中,所述极小值对应的深度位置数据为所述第二位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述第二位置信息中的所述裂缝带下边界深度数据。
在本实施方式中,所述存储器包括但不限于随机存取存储器(Random AccessMemory, RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive, HDD)或者存储卡(Memory Card)。所述存储器可以用于存储计算机程序指令。网络通信单元可以是依照通信协议规定的标准设置的,用于进行网络连接通信的接口。
在本实施方式中,该计算机存储介质存储的程序指令具体实现的功能和效果,可以与其它实施方式对照解释,在此不再赘述。
本说明书实施方式提出识别裂缝带最本质的特点是,根据斯通利波对所述裂缝带敏感的特征,通过所述斯通利波的幅度曲线获取所述裂缝带的第二位置信息;再通过正演的方式,校正所述第二位置信息,获得较为精确的所述裂缝带位置信息。
尽管本申请内容中提到一种裂缝带位置校正方法、装置、电子设备及存储器介质,但是,本申请并不局限于必须是行业标准或实施例所描述的情况等,某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例,仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的装置或模块等,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述模块的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个模块或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,移动终端,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络 PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请。
Claims (9)
1.一种裂缝带位置校正方法,其特征在于,所述方法包括:
获得声波测井数据中的斯通利波信号;
根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;
根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正,得到校正后的位置信息;其中,所述极小值对应的深度位置数据为所述校正后的位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述校正后的位置信息中的所述裂缝带下边界深度数据;
根据所述声波测井数据进行正演,得到第一校正数据和第二校正数据;所述第一校正数据用于校正所述裂缝带上边界深度数据,所述第二校正数据用于校正所述裂缝带下边界深度数据;
根据所述第一校正数据和第二校正数据,校正所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。
2.如权利要求1所述的方法,其特征在于,根据声波测井数据,获得斯通利波信号包括:
根据声波测井数据,获取单极波形数据;
对所述单极波形数据进行低通滤波,获取所述斯通利波信号。
3.如权利要求2所述的方法,其特征在于,所述低通滤波的滤波频率范围为大于等于200Hz且小于等于2000Hz。
4.如权利要求1所述的方法,其特征在于,根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值时包括:
计算指定时间范围内的所述斯通利波信号的振幅平均值;以所述振幅平均值作为该深度位置的斯通利波能量表征值;其中,所述指定时间为所述斯通利波信号的第一个周期和/或第二个周期对应的时间区间。
5.如权利要求1所述的方法,其特征在于,根据声波测井数据正演包括:
根据声波测井数据,建立数值模型;其中,所述声波测井数据至少包括井径大小、目标地层纵波速度、目标地层横波速度;所述数值模型中还至少包括预设裂缝带宽度数据及对应的深度数据;
根据所述数值模型中的斯通利波信号,获得所述第一校正数据和所述第二校正数据;其中,所述第一校正数据为所述数值模型中在深度方向上斯通利波平均幅度极小值对应的所述深度数据与所述预设裂缝带顶部对应的所述深度数据的差值;所述第二校正数据为所述数值模型中在深度方向上斯通利波平均幅度极大值对应的所述深度数据与所述预设裂缝带底部对应的所述深度数据的差值。
6.如权利要求5所述的方法,其特征在于,所述预设裂缝带宽度数据为1mm。
7.一种裂缝带位置校正装置,其特征在于,所述装置包括:
第一处理单元,用于获得声波测井数据中的斯通利波信号;
第二处理单元,用于根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;
第三处理单元,用于根据所述斯通利波能量表征值在深度维度的极大值和极小值,对所述裂缝带的位置信息进行校正,得到校正后的位置信息;其中,所述极小值对应的深度位置数据为所述校正后的位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述校正后的位置信息中的所述裂缝带下边界深度数据;
第四处理单元,用于根据所述声波测井数据进行正演,得到第一校正数据和第二校正数据;所述第一校正数据用于校正所述裂缝带上边界深度数据,所述第二校正数据用于校正所述裂缝带下边界深度数据;
第五处理单元,用于根据所述第一校正数据和第二校正数据,校正所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。
8.一种电子设备,其特征在于,包括:数据接收设备,处理器;
所述数据接收设备,用于接收声波测井数据;
所述处理器,用于获得所述声波测井数据中的斯通利波信号;根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;根据所述斯通利波能量表征值在深度维度的极大值和极小值,对裂缝带的位置信息进行校正,得到校正后的位置信息;其中,所述极小值对应的深度位置数据为所述校正后的位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述校正后的位置信息中的所述裂缝带下边界深度数据;根据所述声波测井数据进行正演,得到第一校正数据和第二校正数据;所述第一校正数据用于校正所述裂缝带上边界深度数据,所述第二校正数据用于校正所述裂缝带下边界深度数据;根据所述第一校正数据和第二校正数据,校正所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。
9.一种计算机存储介质,其特征在于,所述计算机存储介质存储有计算机程序指令,在所述计算机程序指令被执行时实现:
获得声波测井数据中的斯通利波信号;
根据所述斯通利波信号在每个深度位置处的时间数据与斯通利波幅度,得到该深度位置的斯通利波能量表征值;
根据所述斯通利波能量表征值在深度维度的极大值和极小值,对裂缝带的位置信息进行校正,得到校正后的位置信息;其中,所述极小值对应的深度位置数据为所述校正后的位置信息中的所述裂缝带上边界深度数据,所述极大值对应的深度位置数据为所述校正后的位置信息中的所述裂缝带下边界深度数据;根据所述声波测井数据进行正演,得到第一校正数据和第二校正数据;所述第一校正数据用于校正所述裂缝带上边界深度数据,所述第二校正数据用于校正所述裂缝带下边界深度数据;根据所述第一校正数据和第二校正数据,校正所述裂缝带上边界深度数据和所述裂缝带下边界深度数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810029983.XA CN108375785B (zh) | 2018-01-12 | 2018-01-12 | 裂缝带位置校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810029983.XA CN108375785B (zh) | 2018-01-12 | 2018-01-12 | 裂缝带位置校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108375785A CN108375785A (zh) | 2018-08-07 |
CN108375785B true CN108375785B (zh) | 2020-03-10 |
Family
ID=63016654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810029983.XA Active CN108375785B (zh) | 2018-01-12 | 2018-01-12 | 裂缝带位置校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108375785B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110687607B (zh) * | 2019-09-18 | 2021-11-26 | 南方科技大学 | 一种斯通利波探测方法及系统 |
CN113884580B (zh) * | 2020-07-02 | 2023-08-22 | 中国石油天然气股份有限公司 | 激波管实验斯通利波幅度提取方法及装置 |
CN112505780B (zh) * | 2020-10-27 | 2024-05-28 | 中国石油天然气集团有限公司 | 地层深度数据的校正方法及装置 |
CN112666554A (zh) * | 2020-12-17 | 2021-04-16 | 江苏中路工程技术研究院有限公司 | 一种沥青路面雷达振幅特征裂缝宽度识别方法 |
CN117211764B (zh) * | 2023-10-18 | 2024-04-05 | 西南石油大学 | 一种致密气层裂缝宽度随钻录井解释方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104345346B (zh) * | 2014-10-30 | 2017-10-17 | 中国石油天然气股份有限公司 | 一种获取裂缝宽度的方法 |
CN106503353A (zh) * | 2016-10-25 | 2017-03-15 | 中国海洋石油总公司 | 阵列声波测井多孔介质模型定量计算裂缝孔隙度方法 |
CN107587871B (zh) * | 2017-08-07 | 2020-05-08 | 中国石油天然气股份有限公司 | 确定水平裂缝宽度的方法及装置 |
-
2018
- 2018-01-12 CN CN201810029983.XA patent/CN108375785B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108375785A (zh) | 2018-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108375785B (zh) | 裂缝带位置校正方法及装置 | |
Grana | Joint facies and reservoir properties inversion | |
US8427902B2 (en) | Method for monitoring a subsoil zone, particularly during simulated fracturing operations | |
US20150088424A1 (en) | Identifying geological formation depth structure using well log data | |
US6721661B2 (en) | Method of distinguishing types of geologic sedimentation | |
CA2931435A1 (en) | Method for developing a geomechanical model based on seismic data, well logs and sem analysis of horizontal and vertical drill cuttings | |
CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
EA005692B1 (ru) | Зависящая от частоты обработка и интерпретация (fdpi) сейсмических данных для идентификации, изображения и мониторинга насыщенных флюидом подземных месторождений | |
US10024991B2 (en) | Method and device for processing seismic signals | |
US10983235B2 (en) | Characterizing depositional features by geologic-based seismic classification | |
Xu et al. | Method for calculating the fracture porosity of tight-fracture reservoirs | |
US11215035B2 (en) | Method to predict reservoir formation permeability using combined acoustic and multi-frequency dielectric measurements | |
AU2018272313A1 (en) | System and method for predicting fault seal from seismic data | |
Zaree et al. | Estimating fracture intensity in hydrocarbon reservoir: an approach using DSI data analysis | |
CN109100803B (zh) | 微断裂的确定方法和装置 | |
KR101324285B1 (ko) | 물리 검층 데이터의 모델링 방법 | |
RU2716757C2 (ru) | Анализ дифференциальной энергии при дипольном акустическом измерении | |
CN108490490B (zh) | 裂缝带表征方法、装置、电子设备及计算机存储介质 | |
US11346968B2 (en) | Estimation of reservoir flow properties from seismic data | |
US11474272B2 (en) | Methods and systems for identifying and plugging subterranean conduits | |
CN113806674A (zh) | 古河道纵向尺度的量化方法、装置、电子设备及存储介质 | |
CN114755740A (zh) | 岩石分布确定方法、装置、设备及介质 | |
CN111897004A (zh) | 一种基于大数据分析技术的测井预测方法 | |
CN113970786B (zh) | 暗河厚度预测方法、系统、存储介质及电子设备 | |
CN111624662A (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 |