超声层析成像检测系统的研究与实现

    该系统由超声发射接收电路、信号采集、控制电路、信号处理与显示部分组成，如图2所示。整个超声CT成像检测系统由数据(声时)采集、数据处理，以及图像处理显示3部分构成。所需设备和元件包括：若干超声发射探头和接收探头；由单片机组成的脉冲控制电路；2个读写存储器RAM；2个由DSP组成的数字信号处理电路(分别内置射线追踪程序和反演迭代程序)，及DSP驱动的显示阵列LCD。

    超声发射装置是由脉冲控制电路以一定周期发射超声波，接收探头接收信号后，通过高精度计时器得到最小走时矩阵，并存储于存储器1，从而完成数据采集；由存储器1中的最小走时矩阵初步建立(假定超声直线传播)介质内部的慢度矩阵，南慢度矩阵按照射线追踪得到走时路径，存储于存储器2，按照存储器1，2中的信息，经DSP处理，完成方程求解，得到最终精确的慢度矩阵，完成该系统的核心功能一数据处理；重建的慢度矩阵转化为灰度值，由DSP控制LCD显示阵列。最后由工程技术人员对照标准的无缺陷的介质，分析试件，确定有无缺陷，缺陷的位置，大小以及严重情况。

4 图像重建算法

    实现超声层析成像有正演和反演两个基本过程。正演是反演的基础，其模型选取与求解精度直接影响反演变精度；而反演过程实际是对问题的最优化过程。反演方法一般分为变换重建法和级数展开法两大类。由于级数展开法通过离散重建区域，降低原有问题的非线性，适合投影数据不够精确、不能完全且以射线为曲线的层析成像。根据反演的理论基础，层析成像分为以射线理论为基础的射线层析和以波动理论为基础的波动层析。目前CT技术从直射层析成像向弯曲射线层析成像发展。反演方法由最小二乘法发展到各种约束条件下的加权阻尼最小二乘法以及统计法，如最大熵法。观测参数由单一走时数据向多参数数据发展。反演方法则从线性(代数重建法、联合跌代重建法等方法)向非线性(梯队法、模糊神经网络算法和遗传算法等方法)方向发展。射线追踪是确定波从激发点传至接收点的射线路径及其走时的主要技术。因为层析面上各点的波速不相等，波传播的路径是一个与材料不均匀性相关的曲线，其真实路径未知，所以可采用射线追踪理论求取其真实路径。而射线追踪算法分为向前处理与向后处理两部分。

    (1)向前处理先计算离散网格模型上所有边界点的最小走时，然后计算发射点单元所在列各单元边界上所有离散点的最小走时。并记录对应的次级源，接着对发射点单元所在行各单元边界上所有点的计算与第二步相同。在求出每网格内局部走时的基础上，对发射源发出的射线在整个模型上的走时相加便得到全局最小走时。

    (2)向后处理根据向前处理中得到的各单元边界上所有计算点的最小走时及相应次级源，追踪所有发射点到接收点间具有最小走时的射线路径。首先求出接收点所在单元边界上走时最小的离散点，以走时最小的离散点为新的接收点，重复上一步骤向前推，直到发射点所在单元为止。将发射点与最后的射线交点相连。即完成全部向后处理。

5 层析结果

    图4为假设模型的网格分布图，其背景区域的速度为v0=4．000 m／s，异常区(即黑色区域)速度为v1=3 000 m／s。探头布置采用环绕方式，利用线性插值射线追踪算法和联合迭代重建算法得到速度层析图像，结果如图5所示，从图5中可清晰发现低速区。因此，解决了图像检测中图像清晰度不高的困扰。此设计方案可运用到各种需要利用超声波图像检测的领域。

6 结论

    本文所采用的阵列检测方法，在射线追踪算法和SIRT重建方法的基础上，所提取的走时数和超声波数目大大增加，这样在相同的迭代次数条件下，得到结果更精确，重建图像结果更清晰准确