一类反问题—生物体发光断层成像—的理论分析和数值模拟.pdf

摘 要 近年来，随着科学技术的发展，生物医学成像得到了长足的进步．分子成像 是在细胞和分子水平上研究活体生理和病理过程的一个快速发展着的生物医 学成像方法．它的同标是非入侵地准确描述活体细胞和分子过程．分子成像主 要基于三种技术：核素成像技术，核磁共振成像技术和光学成像技术．与传统的 成像方法，如计算机层析成像(CT)、正电子发射断层成像(PET)、核磁共振成像 (MRI)以及它们的合成等，相比，作为光学成像模型，生物放光成像(BLI)有它独 特的优点．对于活体人体研究，由于组织的吸收与发散特性，光学成像主要局限 于表层位置，此时MRj或PET是比较合适的选择．然而，对于小动物来说，由于 较短的路径长度，一大部分光子可以到达动物的外表，从而能以比MRI和PET 低得多的成本探测到生物光和荧光信号．对于活体应用来说，用生物光源进行成 像特别具有吸引力，因为它并不需要外部的刺激光源，从而背景噪音低而灵敏性 高． 生物发光断层成像(BLT)是一个非常有前途的BLI，因为他能实时地揭示 生物体的分子和细胞活动．几年来，人们就BLT问题做了大量的理论分析和数值 仿真．BLT的主要任务就是利用边界测量的光信号来决定小动物体内或一些大 器官的表层的光子密度分布．利用它，我们能在三维空间中对小动物内部的光源 进行定位并定量计算出生物光源的强度密度分布．进行BLT的第一步是确定组 织的光学特性，这正是光学断层摄影(DOT)问题的的主要任务．因为光在任何实 体中的传送都受到吸收和散射的影响，所以我们需要知道光子在生物组织中传 送的精确表达式．一般来说，生物光光子在组织中的传播过程能被漫射转移方程 (RTE)或蒙特卡罗模型(MCM)准确描述．然而，就当前而言，两者在计算上都是 不可行的．因此，当光的波长在600nnl附近或更大时，人们通常用扩散近似方程 来代替RTE．从数学上来讲，BLT是一个反问题，通常需要对它进行正则化处理． 在这篇论文中，我们讨论了BIT和多谱BLT的理论分析和数值求解等问题．全 文共分五章． 在第一章里，我们首先简单介绍一下BLT问题的生物医学背景．然后，在 1．2节中，我们给出了BLT问题的数学形式．在1．3节，我们给出了一些主要符号， 它们将会在接下来的章节里被反复地用到． 第二章主要是关于经典BLT的理论和数值方法的研究．2．1节陈述了经典 BLT的形式．然后，在2．2节里，在一些假定下，利用共轭方法，我们将一个iEN 化的最小化问题转化为一个偏微分方程组．这为快速重构光源函数提供了可能． 在第2．3节中，利用正则化方法，我们研究了BLT问题的一个更为广义框架．在 一类反问题一生物体发光断层成像一 这个框架下，我们建立了问题的适定性结果，考察了它与以前别人文章之间的关 系．在第2．4节中，我们发展了BLT问题的一个新的变分形式．利用该变分形式， 我们可以从理论上证明光源函数并不连续依赖于测量数据．在这一节里，我们还 给出了正则化有限元解最优收敛阶的证明．另外，利用伴随方程，我们构造了一 个简单而有效的迭代格式． 第三章研究的是多谱BLT问题．首先，在3．2节，基于一个惩罚策略，我们 提出了关于多谱BLT的全新的方法．这个新的数学框架的特点是其数值测量的 使用涉及到两个相关但并不相同的边值问题(BVPs)．该数学框架包含了通常的 多谱BLT模型．然后，在第3．2节，我们改进【l】中的误差阶．接着，在第3．3．2小 节里，我们用线性而不是通常的常数有限元来离散光源函数空间．利用有效集策 略，我们获得更高阶的有限元误差估计． 第四章我们继续深入讨论BLT问题．首先，在4．2节，我们研究了具有随空 间改变折射率的BLT问题．为此，利用Tikhonov正则化方法，我们介绍一个广义 的框架，证明了它的适定性，并利用有限元方法(FEMs)，建立了数值解的误差估 计式．在第43节，我们提出了一个数学模型，它将BLT和DOT结合在一起，也 就是说，同时进行两种类型的重构而不是一前一后．凶为在实际中，光学参数通 常被假定是常数或分片光滑函数，因此我们设定光学参数的求解范围是有界变 差(BV)函数空间．我们证明了解的存在性和数值解的收敛性并给出了一个数值 求解格式． 最后，在第五章，我们作了总结并展望未来进一步要做的工作． 关键词：生物发光断层成像，反问题，适定性，有限元方法，误差估计