基于光声泵浦成像的氧分压测量定量分析.pdf

基于光声泵浦成像的氧分压测量定量分析*何霄1)肖小舟1)何滨2)薛平2)肖嘉莹1)1)(中南大学基础医学院生物医学工程系,长沙410083)2)(清华大学物理系,低维量子物理国家重点实验室,北京100084)(2023年 6月 26 日收到;2023年 7月 21 日收到修改稿)光声泵浦成像是一种新型的高特异性光声分子影像技术.它可以避免常规光声成像中来自血液等强背景信号的干扰,实现组织中微弱目标分子的探测,并通过对泵浦-探测激光间的延时扫描,获得组织中的氧分压分布.本文采用亚甲基蓝作为分子探针,通过对血红蛋白溶液中氧分压变化的监测,开展了对光声泵浦成像的定量分析研究.本文采用高斯噪声模型,获得了三重态差分信号稳定性随着平均次数变化的规律,并在此基础上对氧分压测量的误差进行了分析.结果表明,在平均次数为 200 次条件下,氧分压在 300550mmHg(1mmHg=133Pa)区间内,所搭建系统的检测精度优于 33mmHg.本研究将对光声泵浦成像技术的进一步发展和应用起到重要的指导作用.关键词：光声泵浦成像,氧分压,瞬态三重差分,定量分析PACS：81.70.Cv,78.20.Pa,06.20.DkDOI:10.7498/aps.72.202310411引言光声成像(photoacousticimaging,PA)是一种新型的生物医学无损成像方式1.它通过收集脉冲激光照射组织所产生的超声波进行成像,因此兼具了光学成像丰富的分子影像功能和超声成像的高穿透深度的优势,在基础医学和临床医学中显示了出巨大的应用潜力.光声成像作为一种大深度分子影像手段,其主要依赖目标分子的特异性光吸收进行探测2.然而,由于光声成像中入射激光主要被血液中的血红蛋白所吸收,所以常规的光声成像通常用来对组织中的血管网络或血氧分布进行成像.当目标分子浓度较低时,其产生的微弱的光声信号就会被淹没在血液等强背景光声信号中3,4.因此,常规的光声成像仅能对高浓度的造影剂进行成像5.光声泵浦成像是一种基于泵浦-探测的新型光声分子影像技术6,7.由于生物组织分子的激发态寿命通常仅有皮秒或纳秒的量级,激发后会快速衰减到基态.而对于亚甲基蓝等卟啉类分子,其三重态寿命可达微秒量级.因此光声泵浦成像可以通过探测光信号的差分,排除血液等强背景光声信号的干扰,对处于三重态的该类分子实现高选择性的探测.因此,光声泵浦成像可以结合现代分子探针技术,实现大深度、高特异性的目标分子成像8.此外,光声泵浦成像还可以通过对三重态寿命的探测,在不依赖于探针浓度的情况下获得生物体内氧分压的分布9,10.在近期的工作中,本课题组提出了一种新颖的脉冲激光延迟调节技术11,实现了快速自动化的光声泵浦成像,并通过交错采集的方法,有效抑制了激光器能量的波动,显著地提升了光声泵浦成像的速度和精度,使得该技术在光动力治疗监测等临床中的应用成为了可能12.*低维量子物理国家重点实验室开放科研基金项目(批准号:KF202209)、湖南省自然科学基金(批准号:2022JJ30756)、中南大学创新驱动项目(批准号:2020CX004)和深圳市科技创新委员会自由探索基础研究计划(批准号:2021Szvup168)资助的课题.通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.21(2023)218101218101-1然而,目前光声泵浦成像技术的研究仍处于初始阶段.首先,光声泵浦成像需要进行大量的信号平均,以获得高信噪比的三重态差分图像,但是目前缺乏相关的图像稳定性的研究.此外,基于该技术的氧分压监测目前也仍处于定性研究阶段,同样缺乏相关的误差量化分析研究.为此,本文搭建了一套基于亚甲基蓝的光声泵浦氧分压检测系统,通过对牛血红蛋白溶液中氧分压的监测,展开对光声泵浦成像技术的定量研究.本文将对信号平均次数与三重态差分信号稳定性之间的关系进行量化分析,并进一步对所测氧分压的不确定度进行研究,从而对基于光声泵浦的氧分压测量的稳定性和准确性进行评估,为该技术的进一步发展和临床应用提供依据.2研究方法2.1 基于光声泵浦的氧分压测量原理光声泵浦技术原理如图 1 所示.采用两束脉冲激光照射样品,其中泵浦光用来将目标分子激发到激发态,由探测光对处于激发态的目标分子进行探测.当卟啉类分子(如亚甲基蓝)被泵浦光激发时,它会从基态 S0激发到激发单重态 S1,然后通过系间窜越(intersystemcrossing,ISC)转变为三重态(T1),由于自旋禁制的作用,这类 T1态分子相对稳定,寿命为微秒量级.因此在探测光照射时,仍然有大量分子处于 T1态.然而,生物组织内大多数分子激发态寿命仅为皮秒纳秒量级,泵浦光照射后会迅速衰减,处于激发态的分子很少.因此,可通过是否存在泵浦光时所得的光声信号的差分,选择性地获得 T1态分子的光声信号.这样,就可以排除组织内的血液等其他光吸收物质的背景信号,实现对卟啉等具有 T1态的分子的高特异性成像.上述光声差分检测方法称为瞬态三重差分(transienttripletdifferential,TTD),计算公式如下:STTD,t=Spump+probe,t=Spunp,t=Sprobe,t=,(1)STTD,t=Spump+probe,t=Spump,t=Sprobe,t=其中,是通过差分去除背景后得到的 TTD信号,是泵浦和探测光同时激发时的光声信号,是仅泵浦光激发产生的泵浦光信号,是仅探测光激发产生的探测光信号,表示探测激光相对于泵浦激光的延迟时间.pO2值得注意的是,T1态分子可与生物组织中的氧分子相互作用并发生动态淬灭(分子氧得到能量变成单线态氧,T1态分子失去能量返回基态 S0),使处于 T1态分子的寿命进行缩短,因此可通过 T1态的寿命计算得到组织中的氧分压(oxygenpar-tialpressure,).TTD 信号幅值和 T1态寿命的关系为A(t)=A0e(t/T),(2)A(t)A0pO2其中,代表延迟时间=t 时的 TTD 信号,代表 TTD 信号衰减前的幅度,而 T 则是估计得到的 T1态寿命.通过以下公式可将寿命 T 转换为生物组织中的 :T0/T=1+kQT0pO2,(3)T0pO2kQT0kQ其中,代表当 =0mmHg 时的 T1态寿命,而 则是淬灭速率常数.对于亚甲基蓝分子而言,=79.6s,=0.0036s1mmHg1.2.2 光声泵浦氧分压检测系统该系统结构组成见图 2(a),包括两台 OPO 可调谐激光器(SpitLightOPO600)、一台数字延迟/脉冲发生器(DG645,StanfordResearchSystems,CA)、一个 5MHz 水浸超声平探头(中心频率 5MHz、晶元尺寸 6mm、带宽 80%,常州道博超声电子有限公司定制)、脉冲接收/发射器(DPR500,Imagin-ant,NY)以及一张数据采集卡(NIPCI-5124).其中泵浦 OPO 激光器由 355nm 光泵浦,采用 665nm对亚甲基蓝激发,另一台 OPO 激光器由 532nm光泵浦,用于产生 830nm 的探测光.两台激光器重复频率均为 20Hz.泵浦和探测激光器通过一个磷光探针光声信号探测光泵浦光卟啉类分子游离氧分子(O2)时间磷光强度01图1光声泵浦成像原理示意图Fig.1.Illustrationoftheprincipleofphotoacousticpump-probeimaging.物理学报ActaPhys.Sin.Vol.72,No.21(2023)218101218101-2二向色镜耦合到12光纤束中.光纤束的输出端是两个 1.5mm3.5cm 的矩形状输出口,固定在线性超声阵列的两侧,如图 2(b)所示.泵浦激光器在样本上的激光能量密度为 15mJ/cm2,探测激光为 23mJ/cm2,均在人体安全阈值之内13.Spump+probe,t=Spump,t=Sprobe,t=采集数据时的序列流程如图 2(c)所示.在光声泵浦成像中,通过 DG645 在 50ms 的周期内产生 4 种相互之间有设定延时的脉冲,分别作为两台激光器的 Flash 和 Q-switch 触发信号.同时,探测脉冲的触发信号也作为 5124 采集卡的采集触发信号.当激光器出光后,照在样品上产生光声效应,通过超声探头采集信号并传输到 DPR500 将信号放大,接着传输到 5124 采集卡采集并保存到主机.通过上述步骤采集到 ,和 三种序列的数据,再经过差分后获得 TTD信号.由于 OPO 激光器需要稳定工作在 20Hz,所以泵浦光和探测光总会同时出光.为获取仅存在泵浦光或者探测光时候的光声信号,本文采用一种激光脉冲延迟调节技术,使用 DG645 将 OPO 泵浦光或探测光相对于采集卡触发信号向后调制一段时间(如 1ms 左右),这样就可实现三种采集序列的快速自动采集.在此基础上,三种序列采用交错采集的方式,从而避免了激光器能量长周期波动的影响,有利于降低平均次数,提高采集效率.pO2pO2本文中所测样品为掺有 100mg/mL 牛血红蛋白的亚甲基蓝溶液(浓度为 400mol/L).为改变溶液的 ,本文搭建了一套相对密闭的溶液循环装置,包括一个二口烧瓶、一个三口烧瓶、气体流量计、两个蠕动泵及其软管和一个与软管连通的透明塑料装样管.通过气体流量计将氧气和氮气分别通入三口瓶的溶液中,并借助气体流量计调节两种气体的通入速率来改变溶液的 ,如图 2(b)所示.2.3 材料与氧分压实验设计亚甲基蓝是一种芳香杂环化合物,其分子式为 C16H18ClN3S,常被用作化学指示剂、染料等14.作为一种临床上使用的注射药物,已经通过美国食品药品监督管理局的批准15,16,并在近年来被广泛应用于各种光动力治疗的相关研究中17,18.本文确定了亚甲基蓝的 S0态吸收峰为 665nm,其 T1态的最大吸收峰为 830nm 左右.为确定亚甲基蓝T1态的最大吸收峰,本文用 400mol/L 的亚甲基蓝溶液进行了光声光谱扫描实验(平均次数 n=500),并通过(1)式计算得到不同探测波长下的OPO脉冲激光器泵浦激光器探测激光器样品检测数字采集卡数字脉冲延时发射器脉冲接收放大器(a)气体流量计二口烧瓶三口烧瓶蠕动泵蠕动泵装样管探头及光纤固定模块(b)S1S2S3S4探测光泵浦光控制激光器触发信号采集1 ms50 ms50 ms50 ms1 ms(c)图2基于泵浦-探测的组织氧含量体外光声定量检测系统(a)光声泵浦成像系统图;(b)溶液循环装置图;(c)脉冲序列波形示意图Fig.2.Basedonpump-probetechnique,in vitroquantitativedetectionsystemfortissueoxygencontent:(a)Schematicdiagramofphotoacousticpump-probeimagingsystem;(b)diagramofsolutioncirculationdevice;(c)illustrationofpulsesequencewaveform.物理学报ActaPhys.Sin.Vol.72,No.21(2023)218101218101-3TTD 信号,从而获得其 T1态的瞬态光吸收谱.TTD 信号的波动性是光声泵浦成像系统稳定性的重要表征,直接影响到信号平均次数的选择与检测结果的准确性.为考察血红蛋白浓度以及亚甲基蓝浓度对 TTD 信号的影响,本文进行两组对比实验.在第 1 组实验中,样品溶液中的亚甲基蓝浓度为 400mol/L,共设置了 100,50 和 0mg/mL三种不同浓度的牛血红蛋白测试样品.其中,100mg/mL 与人体血红蛋白浓度相似,可用于模拟正常血液条件.第 2 组实验中,牛血红蛋白的含量为 0mg/mL,亚甲基蓝的浓度分别为 50,100,200 和 400mol/L.信号的平均次数为 200 次.pO2pO2pO2为获得系统在不同 下的测量误差,本文采集了 6 组不同 下的光声泵浦数据,其对应的氮气和氧气通入速率的比值分别为 10.25,10.5,10.75,11,11.5 和 12.在每个氮氧比下循环通气 30min 后,采集不同延迟下管中样品的 TTD信号,并通过数据拟合得到 T1态寿命和 .为了平衡采集时间和数据准确性,平均次数 n 设置为200 次.泵浦和探测光之间的延迟时间共有 10 个,分别为 0,0.5,0.75,1,1.5,2,2.75,3.5,4 和 6s.2.4 数据处理与分析方法pO2pO2在光声泵浦成像中,是通过 T1态寿命 T转化得到的,而 T 是通过对不同延时 下采集得到的 TTD 信号拟合得出,因此 测量的准确性和TTD 信号密切相关.TTD 信号既会受两台 OPO激光器脉冲能量波动的影响,也会受到组织内血液等背景信号的干扰,因此需要大量样本数据取平均值,才能获得可靠的数据.基于以上原因,研究平均次数对 TTD 的测量误差之间的关系显得尤为重要.本文采用交错采集的方法获取 TTD 数据,可认为单次测量得到的 TTD 信号符合随机独立分布19.因此,本文通过多次采集得到 TTD 的样本标准差s,来估算平均后 TTD 信号的标准误差(standarderrorofthemean,SEM),计算公式为20SE=s/n,(4)式中,n 为平均次数;TTD 信号的样本标准差 s 计算公式为21s=1n 1ni=1(xi x)2,(5)xi x其中,为样本中第 i 个检测值,为样本均值.pO2由(4)式可知,TTD 平均后的测量误差是随着平均次数 n 的增大而减小,但在实际的 测量中,由于采集的数据量较大,平均次数过高会导致采集时间过长,所以需要平衡采集时间和数据准确度之间的关系,合理地设置平均次数 n.pO2SEpO2pO2pO2As(,p)为对不同 下光声泵浦检测的误差进行评估,本文以 TTD 标准误差 为基础,采用高斯噪声模型,模拟生成不同 不同延时下的测量数据,然后拟合处理得到 数值,通过统计分析获得系统 测量的不确定度.模拟得到的 TTD 数据 如下式所示:As(,p)=Ae(,p)+SE(,p)RAD,(6)Ae(,p)SE(,p)pRADpO2其中,和 分别是氧分压为 、延迟为 条件下,实验测得的 TTD 信号幅值和标准误差,为采用 box-muller 方法生成的“服从标准高斯分布”的随机数22.对每个不同的 ,本文模拟生成了 10000 组测量数据,其中每组测量的信号平均次数为 200 次.3实验结果3.1 亚甲基蓝的瞬态吸收谱亚甲基蓝的基态(S0)光吸收谱可由分光光度计获得,其吸收峰在 665nm,如图 3(a)所示.因此,本文采用 665nm 作为激发波长,测量了亚甲基蓝 T1态的吸收光谱(如图 3(b).实验结果表明,在 700900nm 的波长范围内,MB 溶液的 TTD信号首先随波长的增大而增加,在 755nm 处上升趋势变缓慢,随后又逐渐上升,在 830nm 左右达到峰值,之后逐渐减小.因此,本文最终选择 830nm作为探测光波长.Spump+probe,=0,Spump,=0,Sprobe,=0STTD,=0光声泵浦成像能够通过信号差分的方法去除血液等强背景信号,获得 T1态分子的光声信号.图 4(a)为 400mol/L 亚甲基蓝溶液采集 100 次平均后得到的 和 信号.可以看出,亚甲基蓝的泵浦光和TTD 信号都明显大于探测光,这证明了亚甲基蓝有较高的基态和 T1态激发效率.对于单次测量的TTD 信号,提取其在图 4(a)中 TTD 平均信号峰值处(黑线)的值,可得 100 次测量 TTD 信号幅值的分布,如图 4(b)所示.可以看出,TTD 幅值围绕其均值(红线)在上下波动,其波动幅度代表了测量系统的稳定性.物理学报ActaPhys.Sin.Vol.72,No.21(2023)218101218101-43.2 TTD 信号波动性的影响因素研究图 5(a)显示了不同浓度牛血红蛋白溶液归一化后 TTD 标准差(橘红色虚线),以及归一化后的泵浦光信号的幅值(蓝色虚线).可以看出,随着样品中血红蛋白浓度的加大,其泵浦光信号逐渐增大,同时其 TTD 幅值的波动性也在加大.其中,泵浦光信号的增大是由于溶液光吸收系数的增大.由于各组测量相对独立,因此泵浦光信号的增加必然造成泵浦光信号标准差的增大.而由(1)式和(5)式可知,其也会造成 TTD 信号标准差的增大.另外,当溶液中血红蛋白含量一定,仅亚甲基蓝浓度发生变化时,可见 TTD 的标准差始终保持在 8%以下(图 5(b).这是因为在牛血红蛋白浓度一定的情况下,溶液中亚甲基蓝的浓度变化并不会引起其光吸收系数的太大改变.这两组实验结果表明,在亚甲基蓝浓度较低的条件下,光声泵浦成像的稳定性主要受到样品背景信号波动的影响,背景信号越强,TTD 信号的幅值波动越大.3.3 光声泵浦对氧分压的定量分析研究pO2R2图 6(a)显示了 3 组不同 下,所测量得到的归一化后的 TTD 衰减曲线,TTD 曲线上的不确定度为其标准差.可以看出,随着氮氧比的增大,亚甲基蓝的 T1态寿命明显延长,因此相应 TTD曲线的衰减速率会减慢.还可以看出,3 组 TTD 曲线都显示出非常好的指数衰减趋势,其拟合后曲线的拟合优度 均达到 0.95 以上.pO2pO2由于测量单组 TTD 衰减曲线所需时间(平均200 次)可达 10min,所以本文根据实验测量得到的 TTD 数值及其标准差,通过数值模拟的方法获得了 1 万组 TTD 衰减曲线,来对 测量的准确性进行评估.因此,对于单个 可拟合得到 1 万波长/nm70083090000.51.0峰值拟合曲线原始数据(b)归一化后光声差分信号波长/nm50060070080090000.21.21.00.80.60.4(a)光吸收度2.5 mmol/L5 mmol/L10 mmol/L20 mmol/L25 mmol/L图3亚甲基蓝的基态和激发态吸收光谱(a)基态吸收光谱图;(b)经能量校正后的 T1态吸收光谱图Fig.3.Ground-stateandexcited-stateabsorptionspectraofmethyleneblue:(a)Theabsorptionspectrumofthegroundstate;(b)theenergy-correctedabsorptionspectrumoftheT1state.采集序号02040608010000.250.500.751.00(b)归一化后TTD信号时间/ms01.4153.000-1.0-0.500.51.0(a)归一化后光声信号pump+probepumpprobeTTD图4通过光声泵浦成像得到的光声信号(a)归一化后的 MB 的光声信号示意图;(b)归一化后的 TTD 幅值分布图(基于平均值归一化)Fig.4.Photoacousticsignalsobtainedthroughphotoacousticpump-probeimaging:(a)SchematicrepresentationofthenormalizedphotoacousticsignalofMB;(b)normalizeddistributionofTTDamplitude(normalizedtothemean).物理学报ActaPhys.Sin.Vol.72,No.21(2023)218101218101-5R2个 T1态寿命.图 6(b)为氮氧比为 10.25 的条件下,对模拟所得的 1 万个 T1态寿命排序后得到的累计分布图.图中横坐标是排序后的样本序号,纵坐标为 T1态寿命.可以看出,所得 T1态寿命 T 分布在 0.850.97s 之间.此外,数据拟合结果表明,这 1 万条 TTD 衰减曲线的拟合优度 均在该组原实验数据的2%范围内,表明了在均次数200 次时获得了较好的测量准确性.pO2图 6(c)为根据图 6(b)所得的累积概率积分曲线.对该曲线进行微分即可得概率密度分布曲线,如图 6(d)中红色点状曲线所示.图 6(d)中蓝色曲线为概率密度分布曲线多项式拟合后所得平滑后的曲线,其能克服模拟次数的限制,更好地反映系统 测量值的概率密度分布.可以看出,该组氮氧比下最大似然的 T1态寿命为 0.9094s,相应的 95%置信区间为 0.88330.9375s.pO2pO2pO2pO2pO2通过(3)式将每组模拟所得的 T1态寿命转化为 值,即可得 检测值的统计分布.图 6(e)显示了 6 组不同氮氧比下最大似然的 数值,与其对应的 T1态寿命.图中曲线上标注的误差棒为T1态寿命的 95%置信区间.将图 6(e)中绿色圈所示数据(氮氧比为 11)放大可得图 6(f).可以看出,其 T1态寿命测量值在 0.6410.672s 范围内波动,对应的 在 409.7429.6mmHg 范围内波动,测量误差约为最优值的 4.75%.同理,也得到了其他组的测量误差,分别为 4.53%,5.01%,5.84%,4.75%,5.73%,5.69%.可见,本文中光声泵浦氧分压检测系统的测量误差总体保持在 6%以下,具有较高的测量准确性.4讨论传统的光声成像由于受到强血液背景信号的干扰,其分子影像能力较弱.相对地,光声泵浦成像可以通过对亚甲基蓝等具有长寿命 T1态的探针分子的泵浦-探测,实现深层组织中的微弱目标分子的高特异性成像.此外,光声泵浦成像还可以实现不依赖于探针浓度的组织氧分压的定量检测,在肿瘤、神经和心血管等的临床诊断中有广泛应用前景.然而,光声泵浦成像长期以来一直在成像速度、稳定性和准确性等方面面临挑战.虽然本课题组提出的激光脉冲延迟调节技术,以及其他一系列新型脉冲序列设计方法,有力地提升了其成像速度和质量,但是总体上光声泵浦成像的研究仍处于初级阶段,对其成像稳定性和准确性的定量分析研究还非常少.光声泵浦成像以 TTD 差分成像为基础.在 TTD信号波动性分析中,本文采用高斯噪声模型来评估该技术方法的检测性能.显然,随着平均次数的增加,TTD 信号的噪声会显著减少,进而提升图像的质量和可靠性.然而,需要注意的是,增加平均次数会导致时间分辨率的降低,并且过长时间的连续采集会导致激光出光能量逐渐衰减,进一步影响成像的准确性.因此,在选择平均次数时需要权衡稳定性和时间分辨率之间的关系,根据具体需求进行平衡.在本研究中,采集一组氧分压数据所需的时间为 5min.如果需要提高系统的时间分辨率,需要降低背景血红蛋白浓度,或者降低对系统检测精度的要求.对于氧分压变化本文研究还表明,与不同亚甲基蓝浓度的影响相比,背景信号的强度是影响 TTD 信号稳定性的关键因素.牛血红蛋白浓度/(mgSmL-1)05010000.20.40.60.81.0(a)归一化后的标准差pump,=0TTD,=0亚甲基蓝浓度/(mgSmL-1)5010020040001.00.9归一化后的标准差(b)pump,=0TTD,=0图5TTD 信号波动性分(a)含不同浓度牛血红蛋白的 TTD 标准差;(b)含不同浓度亚甲基蓝溶液的 TTD 标准差Fig.5.AnalysisofTTDsignalvariability:(a)TTDstandarddeviationwithdifferentconcentrationsofbovinehemoglobin;(b)TTDstandarddeviationwithdifferentconcentrationsofmethylenebluesolution.物理学报ActaPhys.Sin.Vol.72,No.21(2023)218101218101-6pO2pO2pO2pO2由于 检测依赖于 TTD 衰减的拟合,因此所测 的准确性与 TTD 信号波动性之间的关系较为复杂.为此,本文以实验数据为基础,通过数值仿真的方法来对 检测的误差进行分析.从结果可以看出,本文中的光声泵浦系统在不同的氧浓度下,其 检测的误差均在 6%之内.值得注意的是,该结果仅考虑了 TTD 信号波动性的影响,即本文主要考虑的是激光能量稳定性以及光声信号本身的涨落的影响.在实际检测中,样品的制备以及人员测量误差也会影响测量结果,这些因素的影响需要在以后的实验中进一步研究.pO2本文的研究结果将有助于对基于光声泵浦的氧分压测量方法的理解,并为该技术的进一步优化和实际应用提供指导.例如在光动力治疗中,光声泵浦成像技术有望对深层组织的光敏剂和 变化进行实时无损成像监测,从而通过术中治疗光的实时间/ms012345600.20.40.60.81.0(a)归一化后TTD信号N2:O2=1:0.25N2:O2=1:0.5N2:O2=1:1O2/mmHg2503203904605306000.50.60.70.80.9T1态寿命/ms(e)T1态寿命/ms0.850.88330.93750.90940.97(c)00.20.40.60.81.0累积概率/ms-195%的置信区间最优值个数02000400060008000100000.850.880.910.940.97(b)O2/mmHg419.40.570.660.730.6410.672409.7429.6T1态寿命/msT1态寿命/ms(f)T1态寿命/ms0102030概率密度0.850.88330.93750.90940.97(d)95%的置信区间最优值图6不同氧分压下的测量误差分析图(a)TTD 衰减曲线示意图;(b)T1态寿命累计分布图;(c)累积概率积分曲线图;(d)概率密度分布曲线图;(e)T1态寿命和氧分压的关系示意图;(f)氮氧比为 11 时 T1态寿命和氧分压波动性示意图Fig.6.Measurementerroranalysisunderdifferentoxygenpartialpressures:(a)SchematicofTTDdecaycurve;(b)cumulativedis-tributionofT1statelifetimes;(c)cumulativeprobabilityintegrationcurve;(d)probabilitydensitydistributioncurve;(e)schematicoftherelationshipbetweenT1statelifetimeandoxygenpartialpressure;(f)illustrationofT1statelifetimefluctuationwithoxygenpartialpressureatnitrogen-to-oxygenratioof11.物理学报ActaPhys.Sin.Vol.72,No.21(2023)218101218101-7pO2时调整,提高光动力治疗效果并减少副作用.因此,对光声泵浦成像方法成像精度和准确性的研究,将有助于对成像速度、泵浦-探测延迟时间,及其他相关监测参数的优化,以及对所得光敏剂和 图像的理解,从而有助于进行合理的治疗干预决策,助力个体化的精准诊疗.虽然本文仅通过仿体实验,对基于光声泵浦的氧分压监测方法进行初步的定量研究,但建立的相关理论框架和数据分析方法,将可用于后续的在体研究中.5结论光声泵浦成像通过差分探测的方法获得探针分子激发态信号,进而通过衰减分析获得样本中的氧分压.这个过程需要经过多个复杂的信号处理步骤.为了定量研究该方法,本文搭建了一套基于光声泵浦的氧分压监测系统,并通过高斯噪声模型,首次建立了基于光声泵浦的氧分压测量误差分析方法,并分析了系统对掺有亚甲基蓝的血红蛋白溶液在不同氮氧比下的氧分压测量误差.结果表明,该系统在 200 次样本平均条件下,在 300550mmHg的氧分压区间内,其检测误差范围保持在 6%以下.本文相关研究结果将推动光声泵浦技术的进一步发展和应用.参考文献 AttiaABE,BalasundaramG,MoothancheryM,DinishUS,Bi R,Ntziachristos V,Olivo M 2019 Photoacoustics 161100144Wang S,Lin J,Wang T F,Chen X Y,Huang P 2016Theranostics623942YaoJJ,WangLV2018Curr.Opin.Chem.Biol.451043LiL,WangLV2021BME Front.202198232684WeberJ,BeardPC,BohndiekSE2016Nat.Methods136395TanJWY,LeeCH,KopelmanR,WangXD2018Sci.Rep.892906SudD,ZhongW,BeerDG,MycekMA2006Opt.Express1444127ShaoQ,AshkenaziS2015J.Biomed.Opt.200360048AshkenaziS,HuangSW,HorvathT,KooYE,KopelmanR2018J.Biomed.Opt.130340239JoJ,LeeCH,FolzJ,TanJWY,WangXD,KopelmanR2019ACS Nano131402410WangB,XieY,HeX,JiangJS,XiaoJY,ChenZY2022Opt.Express303912911CorreiaJH,RodriguesJA,PimentaS,DongT,YangZC2021Pharmaceutics13133212LiL,ZhuLR,MaC,LinL,YaoJJ,WangLD,MaslovK,Zhang R Y,Chen W Y,Shi J H,Wang L V 2017 Nat.Biomed.Eng.1007113GaoL,ZhangC,LiCY,WangLV2013Appl.Phys.Lett.10219370514ZhaoWA,AliMM,BrookMA,LiYF2008Angew.Chem.Int.Ed.47633015GaoXY,LiJH2022Sci.China Chem.521609(inChinese)高晓怡,李景虹2022中国科学:化学52160916Orth K,Beck G,Genze F,Rck A 2000 J.Photochem.Photobiol.B Biol.5718617GrandeMPD,MiyakeAM,NagamineMK,LeiteJVP,da Fonseca I I M,Massoco C O,Dagli M L Z 2022Photodiagn.Photodyn.3710263518Al-TalibM,AlKadiriM,Al-MasriAQ2020Commun.Stat.Theory Methods49562719LeeDK,InJ,LeeS2015Korean.J.Anesthesiology6822020Schillaci M A,Schillaci M E 2022 Evol.Hum.Behav.171103230.21ThistletonWJ,MarshJA,NelsonK,TsallisC2007IEEETrans.Inf.Theory53480522物理学报ActaPhys.Sin.Vol.72,No.21(2023)218101218101-8Quantitative analysis of oxygen partial pressure measure-ments based on photoacoustic pump-probe imaging*HeXiao1)XiaoXiao-Zhou1)HeBin2)XuePing2)XiaoJia-Ying1)1)(Department of Biomedical Engineering,School of Basic Medical Science,Central South University,Changsha 410083,China)2)(State Key Laboratory of Low-Dimensional Quantum Physics,Department of Physics,Tsinghua University,Beijing 100084,China)(Received26June2023;revisedmanuscriptreceived21July2023)AbstractPump-probe-basedphotoacousticimagingisaninnovativetechniqueforhigh-specificitymolecularimagingin deep tissues.Compared with conventional photoacoustic imaging,this method effectively eliminates theinterferencefrombloodsignalandotherbackgroundsignal,enablingthedetectionofsubtletargetmolecules.Additionally,themanipulatingofthetimedelaybetweenthepumplaserandprobelasercanfacilitatenon-invasivemappingofoxygenpartialpressuredistributionwithintissues.Toquantifythephotoacousticpump-probeimaging,weusemethyleneblueasthemolecularprobetomonitorchangesinoxygenpartialpressurewithinahemoglobinsolution.UtilizingaGaussiannoisemodel,weinvestigatetherelationshipbetweenthestabilityofthetriplet-statedifferencesignalandtheaveragenumber,andalsoevaluatetheerrorassociatedwithmeasuringoxygenpartialpressure.Theresultsdemonstratethatthedetectionaccuracyofthesystemisbetterthan33mmHg(1mmHg=133Pa)intheoxygenpartialpressurerangeofabout300to550mmHgafter200 times of averaging.This research will play a significant role in guiding the further advancement andapplicationofpump-probe-basedphotoacousticimagingtechnology.Keywords:pump-probe-based photoacoustic imaging,oxygen partial pressure,transient triplet-differentialmethod,quantitativeanalysisPACS:81.70.Cv,78.20.Pa,06.20.DkDOI:10.7498/aps.72.20231041*ProjectsupportedbytheOpenResearchFundProgramoftheStateKeyLaboratoryofLow-DimensionalQuantumPhysics,China(GrantNo.KF202209),theNaturalScienceFoundationofHunanProvince,China(GrantNo.202