研究 · 本科毕设Research · Undergraduate thesis

微血管光声三维重建与剪切波粘弹性估计3D photoacoustic reconstruction and shear-wave viscoelastic estimation of microvasculature

光声感知与智能成像实验室（OSSII Lab），指导教师马海钢教授，2025.12 – 2026.06Optoacoustic Sensing and Intelligent Imaging Laboratory (OSSII Lab), advised by Prof. Haigang Ma, Dec 2025 – Jun 2026

光声成像Photoacoustic imaging 剪切波Shear wave COMSOL k-Wave

概览Overview

微血管的几何构型与力学特性，对肿瘤、糖尿病等疾病的早期诊断有参考价值。但目前缺少能在真值已知的条件下，系统验证成像与粘弹性反演算法的仿真平台。我搭建了光声（PA）与剪切波（SW）双模态联合仿真框架，在受控条件下评估这些算法的表现。

The geometry and mechanics of microvessels carry information for early diagnosis of diseases such as cancer and diabetes, yet there is no simulation platform that can systematically test imaging and viscoelastic inversion algorithms against a known ground truth. I built a joint photoacoustic (PA) and shear-wave (SW) simulation framework to evaluate these algorithms under controlled conditions.

方法Method

  1. 三维微血管建模：用骨架-体素场方法生成单管、随机网络和真实网络三类模型。3D vessel modelling: skeleton-voxel field models of a single vessel, a random network and a real network.

  2. 多物理场仿真：在 COMSOL 中分别设置压力声学瞬态模块（PA）和固体力学模块（SW）。Multiphysics simulation: transient pressure acoustics (PA) and solid mechanics (SW) in COMSOL.

  3. 传感器阵列：PA 为四侧 256 路线阵，SW 为底面 64 路近场线阵。Sensor arrays: a four-sided 256-channel linear array for PA, a 64-channel near-field array at the bottom for SW.

  4. PA

    光声重建：用 k-Wave 时间反演算法重建三维初始声压分布。Photoacoustic reconstruction: 3D initial pressure recovered with k-Wave time reversal.

  5. SW

    粘弹性估计：滑窗互相关渡越时间法，配合 Maxwell 模型估计剪切波速、模量和粘度。Viscoelastic estimation: sliding-window cross-correlation transit time with a Maxwell model to estimate shear-wave speed, modulus and viscosity.

结果Results

  • PSNR9.70dBPA 重建，随机网络，ROI 区域内PA reconstruction, random network, within the ROI
  • SSIM0.72PA 重建，随机网络，ROI 区域内PA reconstruction, random network, within the ROI
  • 剪切波速Shear-wave speed3.077m/s相对误差 0.76%relative error 0.76%
  • 剪切模量 GShear modulus G9.48kPa相对误差 5.87%relative error 5.87%
  • 粘度 ηViscosity η0.0534Pa·s相对误差 6.3%relative error 6.3%
噪声鲁棒性Noise robustness
在无噪声、40、30、20 dB 四级信噪比下，PA 重建的 SSIM 波动小于 0.01；波速误差在 30 dB 及以上低于 1.5%，20 dB 时约 2.16%。Across noise-free, 40, 30 and 20 dB conditions, PA reconstruction SSIM varied by less than 0.01; wave-speed error stayed below 1.5% at 30 dB and above, and was about 2.16% at 20 dB.

局限Limits

  • 以上均为仿真结果，尚未经过实验验证。All results are from simulation and have not yet been validated experimentally.
  • PSNR 与 SSIM 仅在 ROI 区域内计算。PSNR and SSIM are computed within the ROI only.

材料Materials

  • 论文Thesis
    基于 COMSOL 与 k-Wave 联合仿真的微血管光声三维重建及剪切波粘弹性估计研究3D Photoacoustic Reconstruction and Shear Wave Viscoelastic Estimation of Microvasculature Based on COMSOL and k-Wave Co-Simulation
    本科毕业论文，天津医科大学，2026.05，全文暂不公开Undergraduate thesis, Tianjin Medical University, May 2026. Full text not public.
  • 海报Poster
    基于仿真的微血管光声/剪切波粘弹性三维重建算法研究Simulation-Based PA/SW Dual-Modality 3D Reconstruction for Microvascular Viscoelasticity
    查看海报 PDFView poster (PDF)

← 返回研究← Back to research