发明名称 一种基于高斯回波模型测量薄层材料厚度的方法
摘要 本发明公开了一种基于高斯回波模型测量薄层材料厚度的方法。其步骤包括:1)将薄层材料放置于基体材料表面,使超声波探头分别位于基体材料和薄层材料正上方,测得基体材料和薄层材料的超声波回波信号s<sub>1</sub>(t)、s<sub>2</sub>(t);2)利用高斯回波模型对基体材料回波信号s<sub>1</sub>(t)进行迭代拟合,得到去除噪声后的拟合信号h(t);3)利用脉冲回波模型对薄层材料回波信号s<sub>2</sub>(t)进行迭代拟合,得到拟合结果;4)选取薄层材料回波信号s<sub>2</sub>(t)拟合结果中的前两项回波,两项回波分别对应的脉冲回波模型参数空间中抵达时间分量为τ<sub>1</sub>,τ<sub>2</sub>,已知薄层材料中超声波传播速度为c,则薄层材料的厚度D=(τ<sub>2</sub>-τ<sub>1</sub>)·c/2。本发明可以实现对薄层材料厚度的测量,且可以获得小于采样间距的测量精度。
申请公布号 CN103234494B 申请公布日期 2015.07.15
申请号 CN201310121229.6 申请日期 2013.04.09
申请人 浙江大学 发明人 居冰峰;孙泽青;白小龙;孙安玉;范有成;张建刚
分类号 G01B17/02(2006.01)I;G01N29/07(2006.01)I 主分类号 G01B17/02(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 张法高
主权项 一种基于高斯回波模型测量薄层材料厚度的方法,采用扫描超声波显微镜,扫描超声波显微镜包括超声波探头(1)、三维直线电机(2)、导轨(3)、基体材料(4)、薄层材料(5)、水槽(6)、电机控制器(7)、超声波发射接收器(8)、计算机(9),水槽(6)底部放有基体材料(4),基体材料(4)上设有薄层材料(5),基体材料(4)上方设有超声波探头(1),超声波探头(1)上端与三维直线电机(2)固定,导轨(3)上设有三维直线电机(2),超声波探头(1)与超声波发射接收器(8)相连,三维直线电机(2)与电机控制器(7)相连,计算机(9)分别与电机控制器(7)、超声波发射接收器(8)相连;其特征在于方法的步骤如下:1)将薄层材料(5)放置于基体材料(4)表面,并置于盛有水的水槽(6)中,开启扫描超声波显微镜;2)调节扫描超声波显微镜的三维直线电机(2)的Y轴电机使超声波探头(1)位于基体材料(4)正上方,测量基体材料(4)表面的超声波回波信号s<sub>1</sub>(t);3)建立高斯回波模型,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>;</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>&beta;e</mi><mrow><mo>-</mo><mi>&alpha;</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&tau;</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><mi>&phi;</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000712919270000011.GIF" wi="1058" he="108" /></maths>θ=[ατf<sub>c</sub>φβ]式中f<sub>1</sub>(θ;t)代表高斯回波的时域信号,其中t为时间参量,θ代表高斯回波模型的参数空间,参数空间具体包括衰减参数α,抵达时间参数τ,中心频率参数f<sub>c</sub>,相位参数φ,幅度参数β;4)将M个步骤3)中的高斯回波模型进行线性叠加后记作h(t),以ν(t)表示高斯白噪声,利用h(t)拟合基体材料(4)表面的超声波回波信号s<sub>1</sub>(t):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>m</mi></msub><mo>;</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>h</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>m</mi></msub><mo>;</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000712919270000012.GIF" wi="1186" he="190" /></maths>式中<img file="FDA0000712919270000013.GIF" wi="110" he="85" />代表超声波回波信号s<sub>1</sub>(t)的拟合估计值,θ<sub>m</sub>=[α<sub>m</sub> τ<sub>m</sub> f<sub>cm</sub> φ<sub>m</sub> β<sub>m</sub>]为M个被叠加高斯回波中第m项回波所对应的参数空间取值;5)寻找到薄层材料(5)伸出基体材料(4)表面的部分,即薄层材料(5)上下表面完全浸在水中的部分,并调节扫描超声波显微镜的三维直线电机(2)的Y轴电机使超声波探头(1)位于薄层材料(5)此部分的正上方,测量薄层材料(5)的超声波回波信号s<sub>2</sub>(t);6)基于超声波回波信号s<sub>1</sub>(t)拟合结果h(t)建立超声探头脉冲回波模型,f<sub>2</sub>(ψ;t)=β<sup>*</sup>h(t‑τ<sup>*</sup>)ψ=[β<sup>*</sup> τ<sup>*</sup>]式中f<sub>2</sub>(ψ;t)代表脉冲回波的时域信号,其中t为时间参量,ψ为脉冲回波模型的参数空间,参数空间具体包括幅度参数β<sup>*</sup>,抵达时间参数τ<sup>*</sup>;7)将N个步骤6)中描述的超声探头脉冲回波模型进行线性叠加后记作g(t),以ν<sup>*</sup>(t)表示高斯白噪声,利用g(t)拟合薄层材料(5)的超声波回波信号s<sub>2</sub>(t):<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mover><mi>s</mi><mo>^</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>&psi;</mi><mi>n</mi></msub><mo>;</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>v</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><mi>h</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>*</mo><mo>{</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>&beta;</mi><mi>n</mi><mo>*</mo></msubsup><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msubsup><mi>&tau;</mi><mi>n</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mo>}</mo><mo>+</mo><msup><mi>v</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000712919270000021.GIF" wi="1045" he="411" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>g</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>&psi;</mi><mi>n</mi></msub><mo>;</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000712919270000022.GIF" wi="524" he="196" /></maths>又记:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>g</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>&beta;</mi><mi>n</mi><mo>*</mo></msubsup><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msubsup><mi>&tau;</mi><mi>n</mi><mo>*</mo></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000712919270000023.GIF" wi="610" he="190" /></maths>式中<img file="FDA0000712919270000024.GIF" wi="125" he="82" />代表超声波回波信号s<sub>2</sub>(t)的拟合估计值,<img file="FDA0000712919270000025.GIF" wi="281" he="92" />为N个被叠加脉冲回波中第n项回波所对应的参数空间取值;8)已知薄层材料中超声波传播速度为c,薄层材料的厚度D=(τ<sub>2</sub>‑τ<sub>1</sub>)·c/2。
地址 310027 浙江省杭州市西湖区浙大路38号