CN106597549A - 一种反射波提取方法及装置、反射声波成像测井仪器 - Google Patents

一种反射波提取方法及装置、反射声波成像测井仪器 Download PDF

Info

Publication number
CN106597549A
CN106597549A CN201611104623.9A CN201611104623A CN106597549A CN 106597549 A CN106597549 A CN 106597549A CN 201611104623 A CN201611104623 A CN 201611104623A CN 106597549 A CN106597549 A CN 106597549A
Authority
CN
China
Prior art keywords
receiving unit
acoustic signals
echo
sampling time
time point
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
CN201611104623.9A
Other languages
English (en)
Other versions
CN106597549B (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.)
Shandong Tianyuan Information Technology Co ltd
Shandong Tianyuan Information Technology Group Co ltd
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201611104623.9A priority Critical patent/CN106597549B/zh
Publication of CN106597549A publication Critical patent/CN106597549A/zh
Application granted granted Critical
Publication of CN106597549B publication Critical patent/CN106597549B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using 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)

Abstract

本发明的实施方式提供一种反射波提取方法及装置、反射声波成像测井仪器;该反射波提取方法包括:获取反射声波成像测井仪器中接收器包括的M个接收单元在N个采样时间点接收的声波信号;针对各接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的声波信号,并对m道声波信号进行处理;该处理是针对各采样时间点,将m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值,并将该接收单元接收的声波信号在该采样时间点的采样值与待剔除信号值的差值确定为该接收单元在该采样时间点对应的反射波信号值;输出各接收单元在各采样时间点对应的反射波信号值。本发明对于反射界面与井轴平行的情况仍然可获得良好的应用效果。

Description

一种反射波提取方法及装置、反射声波成像测井仪器
技术领域
本发明的实施方式涉及石油地球物理勘探和声波信号处理领域,更具体地,本发明的实施方式涉及一种反射波提取方法及装置、反射声波成像测井仪器。
背景技术
本部分旨在为权利要求书中陈述的本发明的实施方式提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
随着世界范围内油气勘探开发程度的不断深入,油气勘探逐渐趋向于寻找中小型、隐蔽型和复杂构造油气藏。地震勘探和常规声波测井分别由于分辨率较低和探测深度较浅而无法满足这些油气藏的勘探需求。近年来新兴的反射声波成像测井技术以井中声源辐射到井外地层中的声波信号作为入射波,探测从井旁裂缝、地层界面或小构造等声阻抗不连续界面反射回来的声波信号,通过分析处理接收器接收到的全波列阵列信号,可以对井眼周围的地层构造进行声波成像,以了解井旁地质构造信息,这种新的测井技术在分辨率和探测深度方面恰好介于井间地震和常规声波测井之间,可以用于对井周围数米到数十米范围内的地层构造及地质体进行探测和精细描述,还可用于钻井地质导向,该技术具有良好的应用前景,因而受到国内外学者的广泛关注。
反射声波成像测井技术测量得到的波中包括反射波、井孔模式波及其他干扰模式波,其中,井孔模式波包括幅度较大的滑行纵波、滑行横波和斯通利波等,反射波的幅度相对较小。反射波提取是从反射声波成像测井技术测量得到的波中提取出反射波信号,是反射声波成像测井技术中极为关键的处理步骤之一,是后续处理的基础,反射波提取的效果决定了后续处理结果的准确性和可靠性。但由于反射波的幅度一般远小于井孔模式波的幅度,且反射波和其他干扰模式波在时域上常混叠在一起,给反射波信号的准确提取带来了很大的困难。因此,反射波提取是反射声波成像测井技术的重点和难点。
井孔模式波的幅度和波形特征与井孔模式传播的地层性质有关,也就是与所对应的测量地层有关,因而当源距一定时,反射声波成像测井技术在不同深度位置测得的波中,井孔模式波的幅度和波形特征是不一样的。传统反射波声波测井技术中的反射波提取是针对等源距情况下不同深度位置测得的波进行的,各道测量波中井孔模式波和反射波的特征都与不同深度的测量地层有关,这对于剔除各道测量波中的井孔模式波十分不利。
在过去20多年里,国内外众多学者在反射波提取方面做了大量的研究工作。
Hornby于1989年首次利用F-K滤波法从全波列数据中提取出了反射波,并对井旁倾斜地层界面进行了成像。1999年,Yamamoto等采用了一种新的处理方法,从BARS阵列声波数据中提取了反射波,该方法使用一个反褶积因子对不同道集数据进行处理,以聚焦反射波和直达波,再进行速度滤波和模式波消除。Coates等人于2000年根据反射波和噪声的时差差异进行了滤波处理。Chabot等人使用了径向倾角滤波法对偶极声波成像测井仪DSI测得的全波列数据进行处理以消除S、ST及反射ST波。2002年,Li等人使用VSP的处理方法对阿特拉斯公司的正交偶极阵列声波测井仪XMAC(Elite)在斜井中采集到的全波列数据进行了声反射成像处理,提取了其中的反射波并进行了偏移成像处理,但处理步骤较多,过程比较复杂。Tang等人利用带通滤波消除低频干扰,并指出可以根据反射波和直达波的传播路径不同,使用F-K滤波或者F-K滤波和中值滤波结合分离上、下行反射波,但中值滤波针对的是轴向深度方向的测量波形数据。
发明内容
但是,现有的反射波提取方法,如F-K滤波和中值滤波方法,所针对的数据都是不同深度点的相同源距的波形数据,井孔模式波的到达时间基本相当,若井旁的反射界面相对于井轴来说是倾斜的,则反射波的同相轴也是倾斜的,这种情况下,这些滤波方法是有效的。但是,当井旁的反射界面与井轴是平行时,反射波的同相轴也与井轴是平行,此时F-K滤波和中值滤波等现有的反射波提取方法就无效了。
为此,非常需要一种改进的反射波提取方法。
在本上下文中,本发明的实施方式期望提供一种反射波提取方法及装置、反射声波成像测井仪器。
在本发明实施方式的第一方面中,提供了一种反射波提取方法式,用于对反射声波成像测井仪器测得的声波信号进行处理;
所述反射声波成像测井仪器包括发射器和接收器,所述接收器包括M个按圆形排列且相互之间声电隔离的接收单元;所述接收单元用于接收所述发射器辐射出的在井眼和地层中传播的声波信号;所述发射器位于穿过所述M个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上;
所述反射波提取方法包括:
步骤1,获取所述M个接收单元在N个采样时间点所接收的声波信号;
步骤2,针对每一所述接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号,并对所述m道声波信号进行步骤21-步骤22的处理:
步骤21,针对每一所述采样时间点,将所述m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值;
步骤22,计算该接收单元所接收的声波信号在该采样时间点的采样值与所述待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值;
步骤3,输出每一所述接收单元在每一所述采样时间点对应的反射波信号值;
其中,M、m、N均为自然数,8≦M,3≦m≦M。
在本发明实施方式的第二方面中,提供了一种反射波提取装置,用于对反射声波成像测井仪器测得的声波信号进行处理;
所述反射声波成像测井仪器包括发射器和接收器,所述接收器包括M个按圆形排列且相互之间声电隔离的接收单元;所述接收单元用于接收所述发射器辐射出的在井眼和地层中传播的声波信号;所述发射器位于穿过所述M个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上;
所述反射波提取装置包括:
信号获取模块,用于获取所述M个接收单元在N个采样时间点所接收的声波信号;
信号道读取模块,用于针对每一所述接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号;
第一信号处理模块,用于针对每一所述采样时间点,将所述m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值;
第二信号处理模块,用于计算该接收单元所接收的声波信号在该采样时间点的采样值与所述待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值;
结果输出模块,用于输出每一所述接收单元在每一所述采样时间点对应的反射波信号值;
其中,M、m、N均为自然数,8≦M,3≦m≦M。
在本发明实施方式的第三方面中,提供了一种反射声波成像测井仪器,包括:发射器、接收器以及反射波提取装置;
所述接收器包括M个按圆形排列且相互之间声电隔离的接收单元;所述接收单元用于接收所述发射器辐射出的在井眼和地层中传播的声波信号;所述发射器位于穿过所述M个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上;
所述反射波提取装置包括:
信号获取模块,用于获取所述M个接收单元在N个采样时间点所接收的声波信号;
信号道读取模块,用于针对每一所述接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号;
第一信号处理模块,用于针对每一所述采样时间点,将所述m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值;
第二信号处理模块,用于计算该接收单元所接收的声波信号在该采样时间点的采样值与所述待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值;
结果输出模块,用于输出每一所述接收单元在每一所述采样时间点对应的反射波信号值;
其中,M、m、N均为自然数,8≦M,3≦m≦M。
借助于上述技术方案,本发明提供的反射波提取方法可以有效地提取出反射波信号,特别是对于反射界面与井轴平行的情况,仍然能获得良好的应用效果。本发明提供的反射波提取方法对于确定反射界面的所在方位具有良好的应用前景。
附图说明
通过参考附图阅读下文的详细描述,本发明示例性实施方式的上述以及其他目的、特征和优点将变得易于理解。在附图中,以示例性而非限制性的方式示出了本发明的若干实施方式,其中:
图1为反射声波成像测井的原理示意图;
图2示意性地示出了反射声波成像测井仪器中发射器和接收器的结构示意图;
图3为本发明提供的反射波提取方法;
图4为m道声波信号的采样值罗列成的N×m维数组;
图5示意性地示出了接收器中各个接收单元的编号;
图6为本发明提供的反射波提取装置;
图7为本发明提供的反射声波成像测井仪器的结构框架示意图;
在附图中,相同或对应的标号表示相同或对应的部分。
具体实施方式
下面将参考若干示例性实施方式来描述本发明的原理和精神。应当理解,给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明,而并非以任何方式限制本发明的范围。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
本领域技术技术人员知道,本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此,本公开可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
根据本发明的实施方式,提出了一种反射波提取方法及装置、反射声波成像测井仪器。
在本文中,需要理解的是,附图中的任何元素数量均用于示例而非限制,以及任何命名都仅用于区分,而不具有任何限制含义。
下面参考本发明的若干代表性实施方式,详细阐释本发明的原理和精神。
反射波提取方法
反射声波成像测井技术以井中声源辐射到井外地层中的声波信号作为入射波,探测从井旁裂缝、地层界面或小构造等声阻抗不连续界面反射回来的声波信号,通过分析处理接收器接收到的全波列阵列信号,对井眼周围的小地质构造进行声波成像,以了解井旁地质构造信息。
反射声波成像测井技术可通过反射声波成像测井仪器来实现。
反射声波成像测井仪器包括发射器和接收器。其中,接收器包括按圆形排列且相互之间声电隔离的多个接收单元;发射器位于穿过各个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上;接收单元用于接收发射器发射出的声波信号。
例如图1所示为反射声波成像测井的原理示意图。
在地层101的井眼102中,反射声波成像测井仪器的发射器103辐射声波信号,这些声波信号向井周360°范围内地层辐射,其中部分能量沿着井孔以及井壁外侧附近的地层传播至接收器105的各个接收单元106,承载这部分能量的波就是井孔模式波(包括滑行纵波、滑行横波、伪瑞利波和斯通利波,如图1中的虚线所示);还有部分能量透射到地层中并继续向远离井眼102的地层101中传播,当遇到声阻抗不连续的反射界面104时会产生反射波(如图1中的实线所示),反射波再传播回井中即被接收器105的各个接收单元106接收到。
由于发射器103位于穿过各个接收单元106所在圆形的圆心且与该圆形所在平面垂直的轴上,发射器103与各个接收单元106的距离相等。根据对称性可知,若地层性质相同,则各个接收单元106接收到的井孔模式波的到时、相位和幅度是基本一致的。而井外的反射界面104并不具有对称性,地层101中只有存在反射界面104的方位才会有反射波返回井眼102并被接收单元106接收到;而地层101中无反射界面104的方位,透射到地层中的声波能量会向远离井眼102的方向不断扩散,直至衰减消失。
可见,反射波只会从地层101中存在反射界面104的方位入射至接收器105,由于接收器105的各个接收单元106是按圆形排列的,其中,有些接收单元106的方位与反射界面104的所在方位较接近,有些接收单元106的方位则偏离反射界面104的所在方位,这种情况下,各个接收单元106接收到的反射波的幅度是不一致的。方位角越接近反射界面104的所在方位的接收单元106所接收到的反射波的幅度最大,方位角越偏离反射界面104的所在方位的接收单元106所接收到的反射波的幅度则越小。
在同一时刻,接收器105的各个接收单元106所接收的井孔模式波的幅度基本一致,而所接收的反射波的幅度则有大有小,其中幅度最大处对应的接收单元106即与反射界面104的所在方位最接近,因此,可以先通过确定该幅度最大处对应的接收单元106,然后确定该接收单元106在整个接收器105对应的圆形中所处的方位角,进而来确定反射界面104的所在方位。
反射声波成像测井仪器工作时,接收器105的各个接收单元106所接收的声波信号既包括到达时间、幅度基本一致的井孔模式波,又包括幅度大小不一的反射波,通过剔除其中的井孔模式波,提取反射波,就可以最终确定反射界面104的所在方位。
本发明提供一种反射波提取方法,用于对反射声波成像测井仪器测得的声波信号进行处理,剔除其中的井孔模式波,提取反射波。
如图2所示,反射声波成像测井仪器的接收器包括M个按圆形排列且相互之间声电隔离的接收单元,其发射器位于穿过这M个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上。
如图3所示,该反射波提取方法包括:
步骤S1,获取M个接收单元在N个采样时间点所接收的声波信号。其中,M、N均为自然数,8≦M。
本步骤中,每道声波信号包括N个采样值。
步骤S2,针对每一接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号,并对这m道声波信号进行步骤S21~步骤S22的处理。其中,m为自然数,3≦m≦M。
本步骤中,共读取m道声波信号,由于每道声波信号包括N个采样值,这m道声波信号的采样值可以罗列为如图4所示的N×m维数组。
如图4所示,行表示采样时间点,列表示声波信号的道数,元素aij表示第j道声波信号在采样时间点i的采样值。
步骤S21,针对每一采样时间点,将读取的m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值。
对于采样时间点i来说,读取的m道声波信号在该采样时间点对应的m个采样值分别为ai1…aij…aim,将这m个采样值的中位数标记为Med(i,m),计算该中位数的过程为:将这m个采样值ai1…aij…aim按照从大小顺序排序,若m为奇数,则Med(i,m)为排序结果中间的数值,若m为偶数,则Med(i,m)为排序结果中间的两个数值的平均数。
实际上,对数值序列ai1…aij…aim确定中位数这一过程是把该数值序列中较大或较小的数值过滤掉,最终确定的中位数是该数值序列中大小处于中间位置的数值。
相应的,从对应测量的物理参数上来说,数值序列ai1…aij…aim表示读取的m道声波信号在采样时间点i对应的m个采样值,即表示相应的m个接收单元所接收的声波信号的幅度,该数值序列的中位数表示这m个采样值中大小处于中间位置的数值,这一中位数即表示相应的m个接收单元所接收的声波信号中井孔模式波的幅度。
步骤S22,计算该接收单元所接收的声波信号在该采样时间点的采样值与待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值。
将第j道声波信号在采样时间点i对应的反射波信号值标记为Fij,则
Fij=aij-Med(i,m)
上述方程是从第j道声波信号在采样时间点i对应的采样值中减去待剔除信号值,其物理意义为:从第j道声波信号对应的接收单元所接收的声波信号中剔除井孔模式波,提取反射波。
步骤S3,输出每一接收单元在每一采样时间点对应的反射波信号值。
对所有的接收单元都执行步骤S2,最终可得到每一接收单元在每一采样时间点对应的反射波信号值,利用这些反射波信号值就可以进一步确定反射界面的所在方位。
现有的反射波提取方法,如F-K滤波和中值滤波方法,所针对的数据都是不同深度点的相同源距的波形数据,井孔模式波的到达时间基本相当,若井旁的反射界面相对于井轴来说是倾斜的,则反射波的同相轴也是倾斜的,这种情况下,这些滤波方法是有效的。但是,当井旁的反射界面与井轴是平行时,反射波的同相轴也与井轴是平行或近似平行,此时F-K滤波和中值滤波等现有的反射波提取方法就无效了。
而本发明提供的反射波提取方法可以有效地提取出反射波信号,特别是对于反射界面与井轴平行或近似平行的情况,仍然能获得良好的应用效果。本发明提供的反射波提取方法对于确定反射界面的所在方位具有良好的应用前景。
具体实施时,本发明提供的反射波提取方法可以有多种实施方式,例如:
可选的,可以采用单线程处理方式依次对接收器中的每个接收单元执行步骤S2,也可以多线程处理方式同时对接收器中的多个(例如全部)接收单元执行步骤S2。
可选的,在执行上述反射波提取方法时,可以配合改变反射声波成像测井仪器中发射器和接收器之间的源距,然后再执行步骤S1~步骤S3,这种实施方式可以在不同源距下提取反射波信号,有利于确定反射界面的所在方位。
可选的,在执行上述反射波提取方法时,可以配合改变所述反射声波成像测井仪器测量的深度,然后再执行步骤S1~步骤S3,这种实施方式可以在不同测量深度点提取反射波信号,继而确定不同测量深度点的反射界面的所在方位。
可选的,在执行上述反射波提取方法时,可以配合改变所述反射声波成像测井仪器测量的深度,以及改变所述发射器和所述接收器之间的源距,然后再执行步骤S1~步骤S3,这种实施方式既可以在不同测量深度点提取反射波信号,又可以在不同源距下提取反射波信号,可以有利于确定不同测量深度点的反射界面的所在方位和反射界面的延展。
上述反射波提取方法中的步骤S2中,m可以是奇数,也可以是偶数。
当m为奇数时,步骤S2中针对每一接收单元,读取该接收单元及其两侧各个相邻的接收单元所接收的声波信号。
如图5所示,假设接收器的各个接收单元的编号分别为R1,R2,R3…Rt-2,Rt-1,Rt,Rt+1,Rt+2…RM
例如m=5,读取接收单元Rt,以及接收单元Rt-2,Rt-1,Rt+1,Rt+2所接收的声波信号。
当m为偶数时,步骤S2中针对每一接收单元,读取该接收单元及其一侧个相邻的接收单元、另一侧个相邻的接收单元所接收的声波信号。
例如m=4,读取接收单元Rt,以及接收单元Rt-2,Rt-1,Rt+1所接收的声波信号,或者,读取接收单元Rt,以及接收单元Rt-1,Rt+1,Rt+2所接收的声波信号。
可见,当m为偶数时,步骤S2中针对每一接收单元,要在该接收单元两侧读取不同数目的接收单元的声波信号,考虑到需要针对接收器中的所有接收单元执行步骤S2,为了确保数据处理方式的一致性,在具体执行本发明提供的反射波提取方法时,若m为偶数,则对于每一个接收单元,分别将其两侧定义为第一预定侧和所述第二预定侧,并且,步骤S2中针对每一接收单元,读取该接收单元及其第一预定侧个相邻的接收单元、第二预定侧个相邻的接收单元所接收的声波信号。
例如图5中,假设定义第一预定侧为顺时针侧,第二预定侧为逆时针侧,m=4,则读取接收单元Rt,以及接收单元Rt-1,Rt+1,Rt+2所接收的声波信号。
应当注意,尽管在附图中以特定顺序描述了反射波提取方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。
反射波提取装置
本发明还提供一种反射波提取装置,用于从反射声波成像测井仪器测得的声波信号中提取反射波信号。
本发明提供的反射波提取装置可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
如图6所示,该反射波提取装置包括:
信号获取模块601,用于获取所述M个接收单元在N个采样时间点所接收的声波信号。其中,M、N均为自然数,8≦M。
信号道读取模块602,用于针对每一所述接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号。其中,m为自然数,3≦m≦M。
第一信号处理模块603,用于针对每一所述采样时间点,将所述m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值。
第二信号处理模块604,用于计算该接收单元所接收的声波信号在该采样时间点的采样值与所述待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值。
结果输出模块605,用于输出每一所述接收单元在每一所述采样时间点对应的反射波信号值。
可选的,该反射波提取装置还可以包括:源距控制模块,用于改变所述发射器和所述接收器之间的源距。
可选的,该反射波提取装置还可以包括:测量深度控制模块,用于改变所述反射声波成像测井仪器测量的深度。
可选的,该反射波提取装置还可以既包括以上的源距控制模块,又包括以上的测量深度控制模块。
可选的,m为奇数时,信号道读取模块602用于读取该接收单元及其两侧各个相邻的接收单元所接收的m道声波信号。
可选的,m为偶数时,信号道读取模块602用于读取该接收单元及其一侧个相邻的接收单元、另一侧个相邻的接收单元所接收的m道声波信号。
可选的,该反射波提取装置还可以包括:定义模块,用于第一预定侧和第二预定侧;并且,m为偶数时,信号道读取模块602用于读取该接收单元及其第一预定侧个相邻的接收单元、第二预定侧个相邻的接收单元所接收的m道声波信号,其中,所述第一预定侧与所述第二预定侧为该接收单元的不同侧。
上述反射波提取装置与本发明提供的反射波提取方法基于相同的发明思想实现,其工作原理可参照前述对反射波提取方法的介绍,此处不再赘述。
应当注意,尽管在上文详细描述中提及了反射波提取装置的若干模块,但是这种划分仅仅并非强制性的。实际上,根据本发明的实施方式,上文描述的两个或更多模块的特征和功能可以在一个模块中具体化。反之,上文描述的一个模块的特征和功能可以进一步划分为由多个模块来具体化。
反射声波成像测井仪器
本发明还提供一种反射声波成像测井仪器,如图7所示为该反射声波成像测井仪器的结构框架示意图,该反射声波成像测井仪器包括:发射器701、接收器702以及反射波提取装置703。
其中,发射器用于发射出声波信号。接收器包括M个按圆形排列且相互之间声电隔离的接收单元;每个接收单元用于接收发射器辐射出的在井眼和地层中传播的声波信号。发射器位于穿过这M个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上。
具体实施时,发射器和接收器的实施方式可以有如下形式:
(1)发射器和接收器可以直接设置于钻铤上,以便在随钻环境下进行反射声波成像测量;
(2)发射器和接收器的数量均可以为一个或多个,当发射器和/或接收器的数量为多个时,通过触发不同的发射器和/或接收器,即可完成源距的改变工作;
(3)与接收器包括多个按圆形排列且相互之间声电隔离的接收单元类似,发射器也可以包括多个按圆形排列且相互之间声电隔离的发射单元,通过触发不同的发射单元和接收单元执行相应的功能,可以实现相同源距下对地层的不同方位进行反射声波成像测量;
(4)发射单元/接收单元可以是按圆形均匀排列,以便确定每个发射单元/接收单元的方位;
(5)发射单元/接收单元可以设置于钻铤外表面的刻槽中,以便尽量少占用空间;
(6)发射单元/接收单元可以呈圆弧片状,且其弧度与钻铤外表面的弧度一致,以便适合在钻铤上安装和使用;
(7)在发射单元/接收单元与钻铤之间可以设置隔声材料(例如聚四氟乙烯等)进行声隔离,以尽量避免声波沿钻铤传播。
需要说明的是,在实施本发明时,可以根据实际情况去设置发射器和接收器的装设位置、数量、形状等实施方式,本发明对发射器和接收器的装设位置、数量、形状等实施方式不作具体限定,即以上说明仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,其它任何可行的实施方式均应包含在本发明的保护范围之内。
反射波提取装置包括:
信号获取模块,用于获取所述M个接收单元在N个采样时间点所接收的声波信号;
信号道读取模块,用于针对每一所述接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号;
第一信号处理模块,用于针对每一所述采样时间点,将所述m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值;
第二信号处理模块,用于计算该接收单元所接收的声波信号在该采样时间点的采样值与所述待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值;
结果输出模块,用于输出每一所述接收单元在每一所述采样时间点对应的反射波信号值;
其中,M、m、N均为自然数,8≦M,3≦m≦M。
上述反射波提取装置与本发明提供的反射波提取方法基于相同的发明思想实现,其工作原理可参照前述对反射波提取方法的介绍,此处不再赘述。
具体实施时,反射波提取装置可以具体实现为完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
应当注意,尽管在上文详细描述中提及了反射波提取装置的若干模块,但是这种划分仅仅并非强制性的。实际上,根据本发明的实施方式,上文描述的两个或更多模块的特征和功能可以在一个模块中具体化。反之,上文描述的一个模块的特征和功能可以进一步划分为由多个模块来具体化。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrative logical block),单元,和步骤可以通过电子硬件、电脑软件,或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability),上述的各种说明性部件(illustrative components),单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块,或单元,或装置都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。

Claims (16)

1.一种反射波提取方法,其特征在于,用于对反射声波成像测井仪器测得的声波信号进行处理;
所述反射声波成像测井仪器包括发射器和接收器,所述接收器包括M个按圆形排列且相互之间声电隔离的接收单元;所述接收单元用于接收所述发射器辐射出的在井眼和地层中传播的声波信号;所述发射器位于穿过所述M个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上;
所述反射波提取方法包括:
步骤1,获取所述M个接收单元在N个采样时间点所接收的声波信号;
步骤2,针对每一所述接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号,并对所述m道声波信号进行步骤21-步骤22的处理:
步骤21,针对每一所述采样时间点,将所述m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值;
步骤22,计算该接收单元所接收的声波信号在该采样时间点的采样值与所述待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值;
步骤3,输出每一所述接收单元在每一所述采样时间点对应的反射波信号值;
其中,M、m、N均为自然数,8≦M,3≦m≦M。
2.根据权利要求1所述的反射波提取方法,其特征在于,所述反射波提取方法还包括:改变所述发射器和所述接收器之间的源距,并继续执行所述步骤1~步骤3。
3.根据权利要求1所述的反射波提取方法,其特征在于,所述反射波提取方法还包括:改变所述反射声波成像测井仪器测量的深度,并继续执行所述步骤1~步骤3。
4.根据权利要求1所述的反射波提取方法,其特征在于,所述反射波提取方法还包括:改变所述反射声波成像测井仪器测量的深度,以及改变所述发射器和所述接收器之间的源距,并继续执行所述步骤1~步骤3。
5.根据权利要求1所述的反射波提取方法,其特征在于,m为奇数;
所述的读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号,包括:读取该接收单元及其两侧各个相邻的接收单元所接收的m道声波信号。
6.根据权利要求1所述的反射波提取方法,其特征在于,m为偶数;
所述的读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号,包括:读取该接收单元及其一侧个相邻的接收单元、另一侧个相邻的接收单元所接收的m道声波信号。
7.根据权利要求6所述的反射波提取方法,其特征在于,所述的读取该接收单元及其一侧个相邻的接收单元、另一侧个相邻的接收单元所接收的m道声波信号,包括:
读取该接收单元及其第一预定侧个相邻的接收单元、第二预定侧个相邻的接收单元所接收的m道声波信号,其中,所述第一预定侧与所述第二预定侧为该接收单元的不同侧。
8.一种反射波提取装置,其特征在于,用于对反射声波成像测井仪器测得的声波信号进行处理;
所述反射声波成像测井仪器包括发射器和接收器,所述接收器包括M个按圆形排列且相互之间声电隔离的接收单元;所述接收单元用于接收所述发射器辐射出的在井眼和地层中传播的声波信号;所述发射器位于穿过所述M个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上;
所述反射波提取装置包括:
信号获取模块,用于获取所述M个接收单元在N个采样时间点所接收的声波信号;
信号道读取模块,用于针对每一所述接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号;
第一信号处理模块,用于针对每一所述采样时间点,将所述m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值;
第二信号处理模块,用于计算该接收单元所接收的声波信号在该采样时间点的采样值与所述待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值;
结果输出模块,用于输出每一所述接收单元在每一所述采样时间点对应的反射波信号值;
其中,M、m、N均为自然数,8≦M,3≦m≦M。
9.根据权利要求8所述的反射波提取装置,其特征在于,所述反射波提取装置还包括:源距控制模块,用于改变所述发射器和所述接收器之间的源距。
10.根据权利要求8所述的反射波提取装置,其特征在于,所述反射波提取装置还包括:测量深度控制模块,用于改变所述反射声波成像测井仪器测量的深度。
11.根据权利要求8所述的反射波提取装置,其特征在于,所述反射波提取装置还包括:
源距控制模块,用于改变所述发射器和所述接收器之间的源距;
测量深度控制模块,用于改变所述反射声波成像测井仪器测量的深度。
12.根据权利要求8所述的反射波提取装置,其特征在于,m为奇数;
所述的信号道读取模块用于读取该接收单元及其两侧各个相邻的接收单元所接收的m道声波信号。
13.根据权利要求8所述的反射波提取装置,其特征在于,m为偶数;
所述的信号道读取模块用于读取该接收单元及其一侧个相邻的接收单元、另一侧个相邻的接收单元所接收的m道声波信号。
14.根据权利要求13所述的反射波提取装置,其特征在于,所述的反射波提取装置还包括:定义模块,用于第一预定侧和第二预定侧;
所述的信号道读取模块用于读取该接收单元及其第一预定侧个相邻的接收单元、第二预定侧个相邻的接收单元所接收的m道声波信号,其中,所述第一预定侧与所述第二预定侧为该接收单元的不同侧。
15.一种反射声波成像测井仪器,其特征在于,包括:发射器、接收器以及反射波提取装置;
所述接收器包括M个按圆形排列且相互之间声电隔离的接收单元;所述接收单元用于接收所述发射器辐射出的在井眼和地层中传播的声波信号;所述发射器位于穿过所述M个接收单元所在圆形的圆心且与该圆形所在平面垂直的轴上;
所述反射波提取装置包括:
信号获取模块,用于获取所述M个接收单元在N个采样时间点所接收的声波信号;
信号道读取模块,用于针对每一所述接收单元,读取包括该接收单元在内的m个相邻的接收单元所接收的m道声波信号;
第一信号处理模块,用于针对每一所述采样时间点,将所述m道声波信号在该采样时间点对应的m个采样值的中位数确定为待剔除信号值;
第二信号处理模块,用于计算该接收单元所接收的声波信号在该采样时间点的采样值与所述待剔除信号值的差值,并将该差值确定为该接收单元在该采样时间点对应的反射波信号值;
结果输出模块,用于输出每一所述接收单元在每一所述采样时间点对应的反射波信号值;
其中,M、m、N均为自然数,8≦M,3≦m≦M。
16.根据权利要求15所述的反射声波成像测井仪器,其特征在于,所述接收器中的接收单元按圆形均匀排列。
CN201611104623.9A 2016-12-05 2016-12-05 一种反射波提取方法及装置、反射声波成像测井仪器 Active CN106597549B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611104623.9A CN106597549B (zh) 2016-12-05 2016-12-05 一种反射波提取方法及装置、反射声波成像测井仪器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611104623.9A CN106597549B (zh) 2016-12-05 2016-12-05 一种反射波提取方法及装置、反射声波成像测井仪器

Publications (2)

Publication Number Publication Date
CN106597549A true CN106597549A (zh) 2017-04-26
CN106597549B CN106597549B (zh) 2017-11-14

Family

ID=58595298

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611104623.9A Active CN106597549B (zh) 2016-12-05 2016-12-05 一种反射波提取方法及装置、反射声波成像测井仪器

Country Status (1)

Country Link
CN (1) CN106597549B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108035708A (zh) * 2017-11-20 2018-05-15 中国石油天然气股份有限公司 一种去除地层界面反射波的方法及装置
CN109116415A (zh) * 2018-10-11 2019-01-01 中国石油天然气股份有限公司 地震波数据分离方法、装置及存储介质
CN109116424A (zh) * 2018-10-11 2019-01-01 中国石油天然气股份有限公司 地震波数据的低波数噪音分离方法、装置及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102373923A (zh) * 2010-08-20 2012-03-14 中国石油天然气股份有限公司 一种储层识别方法
CN104950333A (zh) * 2007-05-21 2015-09-30 普拉德研究及开发股份有限公司 用于处理声学波形数据的方法和系统
CN105589103A (zh) * 2014-10-20 2016-05-18 中国石油大学(华东) 井周地质反射体环向扫描成像的探测方法
US20160291189A1 (en) * 2014-11-07 2016-10-06 Halliburton Energy Services, Inc. Apparatus and methods of extracting reflections from acoustic array data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104950333A (zh) * 2007-05-21 2015-09-30 普拉德研究及开发股份有限公司 用于处理声学波形数据的方法和系统
CN102373923A (zh) * 2010-08-20 2012-03-14 中国石油天然气股份有限公司 一种储层识别方法
CN105589103A (zh) * 2014-10-20 2016-05-18 中国石油大学(华东) 井周地质反射体环向扫描成像的探测方法
US20160291189A1 (en) * 2014-11-07 2016-10-06 Halliburton Energy Services, Inc. Apparatus and methods of extracting reflections from acoustic array data

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108035708A (zh) * 2017-11-20 2018-05-15 中国石油天然气股份有限公司 一种去除地层界面反射波的方法及装置
CN108035708B (zh) * 2017-11-20 2021-04-30 中国石油天然气股份有限公司 一种去除地层界面反射波的方法及装置
CN109116415A (zh) * 2018-10-11 2019-01-01 中国石油天然气股份有限公司 地震波数据分离方法、装置及存储介质
CN109116424A (zh) * 2018-10-11 2019-01-01 中国石油天然气股份有限公司 地震波数据的低波数噪音分离方法、装置及存储介质
CN109116424B (zh) * 2018-10-11 2020-06-09 中国石油天然气股份有限公司 地震波数据的低波数噪音分离方法、装置及存储介质

Also Published As

Publication number Publication date
CN106597549B (zh) 2017-11-14

Similar Documents

Publication Publication Date Title
Esmersoy et al. Acoustic imaging of reservoir structure from a horizontal well
CN107810431A (zh) 校准地震检波器阵列的深度的方法
CN105093281B (zh) 一种反演框架下的地震多波建模方法
CN107642114B (zh) 桩基浇注前桩底隐患探查方法及其探查装置
GB2426052A (en) Method for estimating the horizontal shear modulus of a formation while taking into account the presence of an acoustic tool in the borehole.
CN110348135A (zh) 一种随钻声波测井评价地层渗透率的方法
CN106597549B (zh) 一种反射波提取方法及装置、反射声波成像测井仪器
CN103901465A (zh) 全息三维地震勘探观测系统的设计方法
CN108646292B (zh) 裂缝密度预测方法、装置及计算机存储介质
CN105093296A (zh) 一种优化观测系统的方法及装置
CN106593420B (zh) 一种多源组合声反射测井仪及井旁界面探测方法
CN103362502A (zh) 在声波测井中消除直达波干扰的方法、系统及声波测井仪
Pan et al. Reflection signals and wellbore scattering waves in acoustic logging while drilling
CN108594301B (zh) 一种具有差异特征的地震数据融合的方法及处理终端
CN105044779B (zh) 基于相控接收指向性的反射界面方位定量判定方法及装置
CN105403917A (zh) 一种裂缝检测的方法及装置
CN104267434B (zh) 一种三维多分量地震观测系统炮检距分布获取方法及装置
AU2011332299A1 (en) Methods of hydrocarbon detection using spectra dominant frequency and measures of energy decay on the low side and high side of spectra dominant frequency
Morris et al. A new sonic array tool for full waveform logging
CN103941289B (zh) 一种绕射点成像的方法和设备
CN108345032A (zh) 一种弱照明区域高信噪比偏移成像方法
CN104297786B (zh) 检测地层裂缝走向方位角度的方法和装置
CN102914797A (zh) 一种获得地层各向异性系数的方法及装置
CN107229069B (zh) 一种对共检波点数据道集进行速度分析的方法
CN106019372B (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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Che Xiaohua

Inventor after: Zhao Huagang

Inventor after: Qiao Wenxiao

Inventor after: Ju Xiaodong

Inventor after: Zeng Lianbo

Inventor after: Zhao Honglin

Inventor after: Lu Junqiang

Inventor after: Menbaiyong

Inventor before: Che Xiaohua

Inventor before: Qiao Wenxiao

Inventor before: Ju Xiaodong

Inventor before: Zeng Lianbo

Inventor before: Zhao Honglin

Inventor before: Lu Junqiang

Inventor before: Menbaiyong

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20170911

Address after: 257091 South Road, Dongying District, Shandong, No. 228, No.

Applicant after: SHANDONG TIANYUAN INFORMATION TECHNOLOGY CO.,LTD.

Applicant after: China University of Petroleum (Beijing)

Address before: 102249 Beijing city Changping District Road No. 18

Applicant before: China University of Petroleum (Beijing)

GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 257091 No. 228 Nanyi Road, Dongying District, Dongying City, Shandong Province

Patentee after: Shandong Tianyuan Information Technology Group Co.,Ltd.

Patentee after: China University of Petroleum (Beijing)

Address before: 257091 No. 228 Nanyi Road, Dongying District, Dongying City, Shandong Province

Patentee before: Shandong Tianyuan Information Technology Co.,Ltd.

Patentee before: China University of Petroleum (Beijing)

Address after: 257091 No. 228 Nanyi Road, Dongying District, Dongying City, Shandong Province

Patentee after: Shandong Tianyuan Information Technology Co.,Ltd.

Patentee after: China University of Petroleum (Beijing)

Address before: 257091 No. 228 Nanyi Road, Dongying District, Dongying City, Shandong Province

Patentee before: SHANDONG TIANYUAN INFORMATION TECHNOLOGY Co.,Ltd.

Patentee before: China University of Petroleum (Beijing)