CN109283261A - 一种应用于超声相控阵井壁成像检测的回波信号处理方法 - Google Patents

一种应用于超声相控阵井壁成像检测的回波信号处理方法 Download PDF

Info

Publication number
CN109283261A
CN109283261A CN201710589633.4A CN201710589633A CN109283261A CN 109283261 A CN109283261 A CN 109283261A CN 201710589633 A CN201710589633 A CN 201710589633A CN 109283261 A CN109283261 A CN 109283261A
Authority
CN
China
Prior art keywords
signal
borehole wall
echo
ultrasonic phase
phase array
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
Application number
CN201710589633.4A
Other languages
English (en)
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.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
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 Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN201710589633.4A priority Critical patent/CN109283261A/zh
Publication of CN109283261A publication Critical patent/CN109283261A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4472Mathematical theories or simulation
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/002Survey of boreholes or wells by visual inspection
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/449Statistical methods not provided for in G01N29/4409, e.g. averaging, smoothing and interpolation
    • 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
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/0289Internal structure, e.g. defects, grain size, texture
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/10Aspects of acoustic signal generation or detection
    • G01V2210/14Signal detection
    • G01V2210/142Receiver location
    • G01V2210/1429Subsurface, e.g. in borehole or below weathering layer or mud line
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/40Transforming data representation
    • G01V2210/48Other transforms

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • General Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Analytical Chemistry (AREA)
  • Geophysics (AREA)
  • Pathology (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Immunology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Acoustics & Sound (AREA)
  • Remote Sensing (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明涉及一种应用于超声相控阵井壁成像检测的回波信号处理方法,该方法包括:步骤1)利用超声相控阵接收到的井壁回波信号生成相控波束;步骤2)对相控波束采用复小波变换算法进行处理,提取获得包络信号;步骤3)对提取的包络信号进行峰值判断,取包络信号中的最大值作为井壁回波的幅度数据。本发明的方法在进行相控波束形成后,采用复小波变换法对井壁回波信号进行包络提取,该方法能够通过选用与井壁回波信号中心频率相同的复小波基函数,更为精确地提取所需频率的信号,以获得更为平滑、稳定的井壁回波包络信号,使得井壁幅度数据的提取更为准确,从而有效地提高检测系统的信噪比,增强图像对比度,获得更好的井壁成像效果。

Description

一种应用于超声相控阵井壁成像检测的回波信号处理方法
技术领域
本发明涉及超声成像检测技术领域,具体涉及一种应用于超声相控阵井壁成像检测的回波信号处理方法。
背景技术
超声井壁成像检测技术主要利用井壁或套管壁反射回波提取的幅度信息,识别多种类型的缺陷(裂缝和孔洞等),最终以回波幅度图的形式清晰直观地显示井壁的状况。超声相控阵井壁成像检测技术,相较于传统超声井壁成像技术具有避免探头旋转、聚焦点可调等优点,近年来得到了广泛的关注。然而,由于超声信号的非平稳特性,超声相控阵的聚焦特性和井壁回波信号容易受到井下复杂环境的影响,使得井壁回波信号信噪比低,从而造成成像结果的准确性较差。
超声回波信号处理方法是影响最终井壁检测系统成像质量好坏的重要因素之一,超声相控阵井壁回波信号处理主要包括两方面:1)将多路接收到的回波信号进行延时补偿并叠加合成,即实现相控波束形成过程;2)对合成后的回波进行包络提取,在包络信号的基础上提取井壁幅度数据,用于最终的成像显示。目前最常用的井壁回波包络提取方法是采用希尔伯特变换法。该方法简单实用,在信号规律性、周期性及频率特性较强的情况下,采用希尔伯特变换法获取的信号包络清晰、光滑、没有畸变,不需要作进一步的拟合处理。但是,该方法的一个弊端是它会将回波信号中的高次谐波也一同提取出来,产生不光滑的毛刺,降低井壁幅度数据提取的稳定性和准确性,使得最终的井壁成像结果中出现误检的情况,降低超声相控阵井壁成像检测的可靠性。
由上述内容可知,对于超声相控阵井壁成像检测,需要进一步改进井壁超声回波的信号处理方法,以保证超声相控阵井壁成像检测的可靠性,减小井壁成像结果中由于幅度数据提取不准确而造成的误检情况。
发明内容
本发明的目的在于,为解决由于井下复杂环境造成的超声相控阵井壁回波信号信噪比低、成像检测效果差的技术问题,本发明针对超声相控阵井壁成像检测,提出一种井壁超声回波信号处理方法。利用该方法对回波信号进行信号处理,能够获得更为平滑、稳定的井壁回波包络信号,使得井壁幅度数据的提取更为准确,从而有效地提高检测系统的信噪比,增强图像对比度,获得更好的井壁成像显示效果。
为解决上述问题,本发明提供的一种应用于超声相控阵井壁成像检测的回波信号处理方法,该方法包括:
步骤1)利用超声相控阵接收到的井壁回波信号生成相控波束;
步骤2)对相控波束采用复小波变换算法进行处理,提取获得包络信号;
步骤3)对提取的包络信号进行峰值判断,取包络信号中的最大值作为井壁回波的幅度数据。
作为上述技术方案的进一步改进,所述的步骤1)具体包括:
步骤101)利用插值因子对超声相控阵接收到的各路井壁回波信号进行一阶线性插值处理,并按照设定的延时规则对插值处理后的各路信号进行延时补偿;
步骤102)对延时补偿后的多个回波信号进行叠加计算,形成相控波束。
作为上述技术方案的进一步改进,所述步骤101)中经一阶线性插值处理的信号表示为:
其中,x(n)和x(n+1)分别表示第n个和第n+1个动态阵元插值前的回波信号,x′(n)表示第n个动态阵元插值后的回波信号,k和l分别是下一周期采样数据和当前采样数据的权重,且k+l等于系统采样精度与延时精度之比;
所述步骤101)中延时补偿的计算公式表示为:
τi=di/c
其中,τi为第i个阵元的延时,di为聚焦点到第i个阵元的声程,c为井下环境中的声速。
作为上述技术方案的进一步改进,所述步骤102)中叠加计算后形成的相控波束表示为:
其中,(x,z)为扫描平面内成像点位置坐标,si(t)为第i个阵元接收的回波信号。
作为上述技术方案的进一步改进,所述的步骤2)具体包括:
步骤201)对相控波束采用复小波变换算法进行处理得到复信号,复小波变换计算公式表示为:
其中,WTx(a,b)表示对信号x(t)基于参数a与b的小波变换,a和b分别为函数的尺寸及平移因子,φ(t)为小波基函数,t表示时间;根据小波分析的最大匹配原则,选取Morlet复小波作为井壁超声回波信号复小波变换的小波基函数,Morlet复小波函数的数学表达式为:
式中,β为形状控制参数,ω为角频率;
步骤202)对复信号进行求模运算获得包络信号,包络信号的计算公式表示为:
其中,E表示取包络,Re和Im分别表示取实部和取虚部。
本发明的一种应用于超声相控阵井壁成像检测的回波信号处理方法优点在于:本发明的方法在进行相控波束形成后,采用复小波变换法对井壁回波信号进行包络提取,该方法能够通过选用与井壁回波信号中心频率相同的复小波基函数,更为精确地提取所需频率的信号,以获得更为平滑、稳定的井壁回波包络信号,使得井壁幅度数据的提取更为准确,从而有效地提高检测系统的信噪比,增强图像对比度,获得更好的井壁成像显示效果。
附图说明
图1为本发明提供的应用于超声相控阵井壁成像检测的回波信号处理方法流程图。
具体实施方式
下面结合附图和实施例对本发明所述的一种应用于超声相控阵井壁成像检测的回波信号处理方法进行详细说明。
本发明针对超声相控阵井壁成像检测,提出一种井壁超声回波信号处理方法,用于解决由井下复杂环境造成的超声相控阵井壁回波信号信噪比低、成像检测效果差的技术问题。
本发明的方法在进行相控波束形成后,采用复小波变换法对井壁回波信号进行包络提取,以获得更为平滑、稳定的井壁回波包络信号,使得井壁幅度数据的提取更为准确,从而有效地提高检测系统的信噪比,增强图像对比度,获得更好的井壁成像显示效果。
本发明提供的一种应用于超声相控阵井壁成像检测的回波信号处理方法,该方法具体包括:
步骤1)利用超声相控阵接收到的井壁回波信号生成相控波束;
步骤2)对相控波束采用复小波变换算法进行处理,提取获得包络信号;
步骤3)对提取的包络信号进行峰值判断,取包络信号中的最大值作为井壁回波的幅度数据。
基于上述回波信号处理方法,如图1所示,在本实施例中,超声相控阵井壁成像检测的回波信号处理基本步骤包括:相控波束形成(即多路回波信号的延时、叠加)、对波束形成后的回波信号进行包络提取,以及在获得井壁回波信号包络的基础上提取井壁幅度数据,用于最终的成像显示。
下面对本发明方法中的相控波束形成、包络提取和幅度数据提取操作分别进行详细说明。
(1)相控波束形成
相控波束形成的具体步骤包括:
步骤101)根据实际的相控阵检测系统延时精度,利用插值因子对超声相控阵接收到的n路原始井壁回波信号进行一阶线性插值处理,只保留与延时有关的插值点上的波形数值,并按照设置的延时规则对n路数据进行延时补偿。
原始的采样数据x(n)经过插值处理后表示为:
上式中,x(n)和x(n+1)分别表示第n个和第n+1个动态阵元插值前的回波信号,x′(n)表示第n个动态阵元插值后的回波信号,k和l分别是下一周期采样数据和当前采样数据的权重,且k+l等于系统采样精度与延时精度之比。
然后按照设置的延时规则对n路数据进行延时补偿,延时规则由换能器孔径形状和聚焦点位置决定,在本实施例中,延时补偿的计算公式表示为:
τi=di/c
其中,τi为第i号阵元的延时,di为聚焦点到第i号阵元的声程,c为井下泥浆环境中的声速。
步骤102)对延时补偿后的多个回波信号进行叠加计算,叠加后的回波信号表达式为:
其中,I(x,z)表示叠加计算后形成的相控波束,(x,z)为扫描平面内成像点位置坐标,si(t)为第i号阵元接收的回波信号。
(2)包络提取
包络提取的具体步骤包括:
步骤201)对相控波束形成后的回波进行包络提取,这里采用复小波变换提取回波信号包络,得到复信号。对于回波信号x(t),其复小波变换计算公式表示为:
其中,WTx(a,b)表示对信号x(t)基于参数a与b的小波变换,a和b分别为函数的尺寸及平移因子,φ(t)为小波基函数,t表示时间。根据小波分析的最大匹配原则,在本实施例中,选取Morlet复小波作为井壁超声回波信号复小波变换的小波基函数,Morlet复小波函数的数学表达式为:
上式中,β为形状控制参数,ω为角频率。
步骤202)在复小波变换获得信号x(t)对应的复信号表示的基础上,对该复信号进行求模运算,即获得原信号x(t)的信号包络,计算公式为:
其中,E表示取包络,Re和Im分别表示取实部和取虚部。
(3)幅度数据提取
为获取井壁回波的幅度数值,具体处理方法是对提取的包络信号进行峰值判断,即取x(t)中的最大值作为井壁回波的幅度数据,完成信号处理操作。
在实际系统当中,假设采样后的数字序列为x(n),若存在点x(i),使得在序列x(n)当中,对于任意n(n≠i),都有x(i)>x(n),则取x(i)为该扫描点处的井壁回波幅度数据。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (5)

1.一种应用于超声相控阵井壁成像检测的回波信号处理方法,其特征在于,该方法包括:
步骤1)利用超声相控阵接收到的井壁回波信号生成相控波束;
步骤2)对相控波束采用复小波变换算法进行处理,提取获得包络信号;
步骤3)对提取的包络信号进行峰值判断,取包络信号中的最大值作为井壁回波的幅度数据。
2.根据权利要求1所述的应用于超声相控阵井壁成像检测的回波信号处理方法,其特征在于,所述的步骤1)具体包括:
步骤101)利用插值因子对超声相控阵接收到的各路井壁回波信号进行一阶线性插值处理,并按照设定的延时规则对插值处理后的各路信号进行延时补偿;
步骤102)对延时补偿后的多个回波信号进行叠加计算,形成相控波束。
3.根据权利要求2所述的应用于超声相控阵井壁成像检测的回波信号处理方法,其特征在于,所述步骤101)中经一阶线性插值处理的信号表示为:
其中,x(n)和x(n+1)分别表示第n个和第n+1个动态阵元插值前的回波信号,x′(n)表示第n个动态阵元插值后的回波信号,k和l分别是下一周期采样数据和当前采样数据的权重,且k+l等于系统采样精度与延时精度之比;
所述步骤101)中延时补偿的计算公式表示为:
τi=di/c
其中,τi为第i个阵元的延时,di为聚焦点到第i个阵元的声程,c为井下环境中的声速。
4.根据权利要求3所述的应用于超声相控阵井壁成像检测的回波信号处理方法,其特征在于,所述步骤102)中叠加计算后形成的相控波束表示为:
其中,(x,z)为扫描平面内成像点位置坐标,si(t)为第i个阵元接收的回波信号。
5.根据权利要求1所述的应用于超声相控阵井壁成像检测的回波信号处理方法,其特征在于,所述的步骤2)具体包括:
步骤201)对相控波束采用复小波变换算法进行处理得到复信号,复小波变换计算公式表示为:
其中,WTx(a,b)表示对信号x(t)基于参数a与b的小波变换,a和b分别为函数的尺寸及平移因子,φ(t)为小波基函数,t表示时间;根据小波分析的最大匹配原则,选取Morlet复小波作为井壁超声回波信号复小波变换的小波基函数,Morlet复小波函数的数学表达式为:
式中,β为形状控制参数,ω为角频率;
步骤202)对复信号进行求模运算获得包络信号,包络信号的计算公式表示为:
其中,E表示取包络,Re和Im分别表示取实部和取虚部。
CN201710589633.4A 2017-07-19 2017-07-19 一种应用于超声相控阵井壁成像检测的回波信号处理方法 Pending CN109283261A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710589633.4A CN109283261A (zh) 2017-07-19 2017-07-19 一种应用于超声相控阵井壁成像检测的回波信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710589633.4A CN109283261A (zh) 2017-07-19 2017-07-19 一种应用于超声相控阵井壁成像检测的回波信号处理方法

Publications (1)

Publication Number Publication Date
CN109283261A true CN109283261A (zh) 2019-01-29

Family

ID=65184880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710589633.4A Pending CN109283261A (zh) 2017-07-19 2017-07-19 一种应用于超声相控阵井壁成像检测的回波信号处理方法

Country Status (1)

Country Link
CN (1) CN109283261A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112834623A (zh) * 2021-01-06 2021-05-25 浙江大学 一种基于压缩感知技术的声学覆盖层回声降低测量方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102507739A (zh) * 2011-11-15 2012-06-20 北京理工大学 碳纤维复合材料冲击损伤的超声导波缺陷成像方法
CN102539542A (zh) * 2011-11-15 2012-07-04 北京理工大学 复小波变换法提取陶瓷基芳纶纤维增强复合板中单模式超声导波信息

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102507739A (zh) * 2011-11-15 2012-06-20 北京理工大学 碳纤维复合材料冲击损伤的超声导波缺陷成像方法
CN102539542A (zh) * 2011-11-15 2012-07-04 北京理工大学 复小波变换法提取陶瓷基芳纶纤维增强复合板中单模式超声导波信息

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李杏华 等: "高精度波束形成的相控阵超声系统研究", 《传感技术学报》 *
董晗 等: "超声相控阵井壁成像检测的回波信号处理", 《无损检测》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112834623A (zh) * 2021-01-06 2021-05-25 浙江大学 一种基于压缩感知技术的声学覆盖层回声降低测量方法

Similar Documents

Publication Publication Date Title
US11255708B2 (en) Extended range ADC flow meter
EP2137547B1 (en) System and method for positioning
WO2017167191A1 (zh) 地震数据处理方法和装置
CN106569196B (zh) 一种基于压缩感知的地面雷达多目标检测方法
CN107179535A (zh) 一种基于畸变拖曳阵的保真增强波束形成的方法
US7512518B2 (en) Method for measuring thin layers in solid state devices
CN103760522A (zh) 用于时差估计与多站时钟误差校准的方法及系统
CN105675126B (zh) 一种用于检测多频多源复杂稳定声场声压的新方法
CN110231620A (zh) 一种噪声相关系统跟踪滤波方法
CN108186050B (zh) 一种基于超声通道数据的多普勒血流速度成像方法和系统
CN104748704A (zh) 薄壁结构超声共振测厚频谱分析内插校正方法
CN106569180B (zh) 一种基于Prony方法的方位估计算法
CN112230212B (zh) 雷达测距信号处理方法和装置
CN110048795A (zh) 一种地震仪采集数据时钟的方法及装置
Handegard et al. Tracking individual fish from a moving platform using a split-beam transducer
CN113009572B (zh) 一种基于横波偏振分析预测裂缝方位角的方法
CN102636773B (zh) 基于信道多途特性的单基元抗距离模糊方法
CN102540153A (zh) 一种基于站间直达波干扰的阵列幅相误差校准方法
CN109283261A (zh) 一种应用于超声相控阵井壁成像检测的回波信号处理方法
CN110196407A (zh) 一种基于频率预估的单矢量水听器信号来波方向估计方法
CN103196451B (zh) 一种基于Fisher信息矩阵的导航脉冲星选择方法
CN110850386A (zh) 一种基于分数阶域特征的旋翼类无人机深度学习识别方法
CN105527624B (zh) 一种雷达回波动态估算噪声的方法和气象雷达系统
CN106093941A (zh) 一种基于fmcw雷达系统的冰层截面成像方法
Russell et al. OBSrange: A new tool for the precise remote location of ocean‐bottom seismometers

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190129

RJ01 Rejection of invention patent application after publication