CN114485746B - 基于时分复用多载波探测光干涉型传感器的光声成像系统 - Google Patents
基于时分复用多载波探测光干涉型传感器的光声成像系统 Download PDFInfo
- Publication number
- CN114485746B CN114485746B CN202111603250.0A CN202111603250A CN114485746B CN 114485746 B CN114485746 B CN 114485746B CN 202111603250 A CN202111603250 A CN 202111603250A CN 114485746 B CN114485746 B CN 114485746B
- Authority
- CN
- China
- Prior art keywords
- signal
- type sensor
- interference type
- photoacoustic
- laser
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 46
- 238000001514 detection method Methods 0.000 title claims abstract description 34
- 230000003287 optical effect Effects 0.000 claims abstract description 77
- 238000005070 sampling Methods 0.000 claims description 57
- 150000004770 chalcogenides Chemical class 0.000 claims description 41
- 230000010287 polarization Effects 0.000 claims description 39
- 230000004044 response Effects 0.000 claims description 35
- 239000000523 sample Substances 0.000 claims description 25
- 238000001228 spectrum Methods 0.000 claims description 24
- 238000000034 method Methods 0.000 claims description 13
- 239000000463 material Substances 0.000 claims description 11
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 4
- 230000001678 irradiating effect Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 2
- 238000005305 interferometry Methods 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 230000035945 sensitivity Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 206010034972 Photosensitivity reaction Diseases 0.000 description 3
- 238000012634 optical imaging Methods 0.000 description 3
- 230000036211 photosensitivity Effects 0.000 description 3
- 239000006096 absorbing agent Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000012879 PET imaging Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000000411 transmission spectrum Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/35306—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
- A61B5/0095—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying light and detecting acoustic waves, i.e. photoacoustic measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
- G01N2021/458—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods using interferential sensor, e.g. sensor fibre, possibly on optical waveguide
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Acoustics & Sound (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明提出基于时分复用多载波探测光干涉型传感器的光声成像系统，涉及光声成像的技术领域，以光干涉型传感器作为基本结构，首先利用信号发生器输出若干正弦波射频信号，利用可调激光器产生探测激光，IQ调制器基于正弦波射频信号对可调激光器产生的探测激光进行单边带调制，生成多路载波，基于光干涉型传感器实现光干涉型传感器的输出光谱在10GHz范围内，与IQ调制器和光电探测器的带宽限制配合，提高了探测系统的时域利用率，而且光干涉型传感器的输出光谱中同时存在多个谐振峰，同时对输出光谱中的多个谐振峰通道进行光强度监测，进而获得高信噪比的光声信号，实现了高灵敏度的光声成像。
Description
技术领域
本发明涉及光声成像的技术领域，更具体地，涉及一种基于时分复用多载波探测光干涉型传感器的光声成像系统。
背景技术
光声成像是一种非标记、无损伤的成像技术，具有很高的光学对比度和成像分辨率，因此在生物、医学等领域发挥着重要作用。当利用短脉冲激光照射生物组织时，生物组织中的分子吸收光子后，能量会以热量的形式释放能量。释放的热量导致吸收体局部温度升高，温度升高后导致热膨胀而产生压力波，进而产生光声信号。光声成像结合了光学成像和声学成像的优点。一方面，比纯光学成像穿透更深，通过聚焦光束的方式，可激发组织深层的光声信号；另一方面，比传统的MRI以及PET成像拥有更高的分辨率，其图像分辨率可达到亚微米、微米量级。
超声波传感器属于光声成像系统的核心器件，传感器的性能直接决定了成像系统的效果。光干涉型传感器为光学成像系统最有前景的核心器件，原因在于如下几个因素:光干涉型传感器制备与目前的CMOS技术是兼容的、其器件尺寸可以做到极小、同时具有很高的品质因数。基于光干涉型传感器的光声成像系统具有极高的灵敏度、大的信噪比和很高的成像精度，同时也避免了压电传感器因为射频串扰而面临的高频局限的问题。当声波作用于光干涉型传感器时，由于形变效应和弹光效应的影响，微腔模式的有效折射率改变引起谐振条件变化，导致输出光谱中谐振峰位置漂移。目前主流系统中一般使用单频激光器探测并解调光谱中的超声脉冲响应信息：通过强度反馈波长跟踪技术，将激光器波长锁定于谐振某一边带的一半光强的波长处，以获取超声信号的最大线性响应；当微腔受外界超声信号的作用，发生光谱谐振峰漂移时，具有极高斜率的边带，可以将微弱的谐振漂移量转化为更明显的光强信号。如现有技术中公开了一种基于数字光频梳与微腔阵列的高灵敏超声探测方法，利用探测信号发射单元发出合适的探测信号，由连续激光器发射频率为f0的单频连续光信号，经过频移器，频移器受任意信号发生器调制；然后，探测光进入传输链路进行传输，经过调制的、具有 BW带宽的光信号经过逛耦合器，耦合进入光路系统，光信号经过传输链路入射到微腔阵列传感单元；光信号以光速快速依次分别经历多个具有高品质因子的微腔，得到微腔阵列的透过光谱，该方法可用于光声成像中。然而，一方面，使用单频激光器探测探测得到的光声信号的信噪比取决于传感元件的灵敏度和系统噪声等，在传感元件以及系统中各种设备固定的情况下，通常很难直接对系统信噪比进行提升；另一方面，探测光强变化的光电探测器和示波器的带宽可达到几 GHz甚至几十GHz，远高于光声信号的频率范围，探测系统的时域利用率较低。
发明内容
为解决传统单载波探测并解调光谱中的超声脉冲响应信息，从而用于光声成像的方式时域利用率低、实现成本高且探测光声信号信噪比低的问题，本发明提出一种基于时分复用多载波光干涉型传感器的光声成像系统，通过时分复用的多载波探测技术，同时对谐振范围的多个波长通道进行光强度监测，有效控制成本，提高时域利用率，进而获得高信噪比的光声信号，实现了灵敏度的光声成像。
为了达到上述技术效果，本发明的技术方案如下：
一种基于时分复用多载波探测光干涉型传感器的光声成像系统，包括：信号发生器、可调激光器、IQ调制器、偏振控制器、光干涉型传感器、光电探测器及示波器；信号发生器在每个采样周期内，按n个等时间采样间隔依次输出n个频率的多载波正弦波射频信号：f1，…，fn，n个正弦波射频信号输入IQ调制器，可调激光器产生中心频率为f0的探测激光，IQ调制器基于n个正弦波射频信号对可调激光器产生的探测激光进行单边带调制，在每个采样周期内产生中心频率为f0+f1、…、f0+fn的信号激光，并输入偏振控制器进行偏振态调整，通过带隙光聚焦照射，调整光干涉型传感器输出光谱中同时产生若干个谐振峰，且谐振峰中心频率为f0，调整偏振态后的信号激光作为探针，监测光干涉型传感器在每个采样周期内的n个等时间采样间隔对生物组织内激发的光声信号的时域响应信号；之后，作为探针的信号激光进入光电探测器，光强度转化为电压强度信号后传输至示波器，携带光干涉型传感器对光声信号的时域响应信号被示波器接收。
在技术方案中，以光干涉型传感器作为基本结构，首先利用信号发生器输出若干正弦波射频信号，利用可调激光器产生探测激光，IQ调制器基于正弦波射频信号对可调激光器产生的探测激光进行单边带调制，生成了多路子载波，即信号激光，有效控制了成本，避免了多路光源复用与解复用的复杂系统的应用。基于光干涉型传感器的硫系材料光敏性，实现光干涉型传感器的输出光谱在10GHz 范围内，与IQ调制器和光电探测器的带宽限制配合，提高了探测系统的时域利用率，而且光干涉型传感器的输出光谱中同时存在多个谐振峰，通过时分复用的多路载波探测技术，同时对谐振峰范围的多个波长通道进行光强度监测，进而获得高信噪比的光声信号，实现了高灵敏度的光声成像。
优选地，设每个采样周期时长表示为T，按n个等时间间隔依次输出n个正弦波射频信号时，每一个正弦波射频信号与采样间隔的对应关系为：在(0～1/n)T 时间间隔内，输出频率为f1的正弦波射频信号，依次类推，在(n-1/n～1)T时间间隔内，输出频率为fn的正弦波射频信号。
优选地，所述可调激光器为单一的c波段可调激光器。在此，后续通过调制的方式生成多路载波信号激光，有效控制了成本，也避免了多路光源复用与解复用的复杂系统的应用。
优选地，光干涉型传感器为基于硫系材料的片上光学微腔阵列传感器件，包括n个硫系微环以及一根将n个硫系微环串联耦合的直波导，所述n个硫系微环的尺寸相同，利用带隙光分别聚焦照射每一个硫系微环，调节该硫系微环产生的谐振峰的位置，使得光干涉型传感器的输出光谱中所有谐振峰的中心频率集中于 10GHz范围以内，可以避免在光声成像过程中需要的额外芯片电极工艺以及持续的电压驱动。
优选地，偏振控制器调整的偏振态为TE偏振或TM偏振。对于一个光干涉型传感器件，输入不同偏振态的激光，会得到不同的光谱；如果输入混合偏振态的光源，会得到混乱的光谱，这对于传感探测是不好的；因此一般控制输入单一的偏振态。
优选地，生物组织内激发的光声信号由激光器聚焦激发光源照射生物组织产生。
优选地，监测光干涉型传感器在每个采样周期内的n个等时间采样间隔对生物组织内激发的光声信号的时域响应时，提取n个硫系微环中每一个硫系微环对光声信号的时域响应信号。
优选地，在提取n个硫系微环中每一个硫系微环对光声信号的时域响应信号后，得到n个时域响应信号，对n个时域响应信号作采样周期对齐，并做平均处理，以降低随机噪声。
一种根据所述的基于时分复用多载波探测光干涉型传感器的光声成像系统实现的光声成像方法，所述方法包括：
S1.确定采样周期，利用信号发生器在每个采样周期内，按n个等时间采样间隔依次输出n个频率的多载波正弦波射频信号：f1，…，fn，将n个正弦波射频信号输入IQ调制器；
S2.利用可调激光器产生中心频率为f0的探测激光，IQ调制器基于n个正弦波射频信号对可调激光器产生的探测激光进行单边带调制，在每个采样周期内产生中心频率为f0+f1、…、f0+fn的信号激光，并输入偏振控制器调整偏振态；
S3.通过带隙光聚焦照射，调整光干涉型传感器的输出光谱同时产生若干个谐振峰，且谐振峰中心频率为f0；
S4.调整偏振态后的信号激光作为探针，监测光干涉型传感器在每个采样周期内的n个等时间采样间隔对生物组织内激发的光声信号的时域响应信号
S5.作为探针的信号激光进入光电探测器，光强度转化为电压强度信号后传输至示波器，携带光干涉型传感器对生物组织内激发的光声信号的时域响应信号被示波器接收。
优选地，步骤S5之后还包括：对时域响应信号作采样周期对齐，并做平均处理，以降低随机噪声。
与现有技术相比，本发明技术方案的有益效果是：
本发明提出一种基于时分复用多载波探测光干涉型传感器的光声成像系统，以光干涉型传感器作为基本结构，首先利用信号发生器输出若干正弦波射频信号，利用可调激光器产生探测激光，IQ调制器基于正弦波射频信号对可调激光器产生的探测激光进行单边带调制，生成了多路载波，采用调制方式有效控制了成本，避免了多路光源复用与解复用的复杂系统的应用。基于光干涉型传感器的硫系材料光敏性，实现光干涉型传感器的输出光谱在10GHz范围内，与IQ调制器和光电探测器的带宽限制配合，提高了探测系统的时域利用率，而且光干涉型传感器的输出光谱中同时存在多个谐振峰，通过时分复用的多路载波探测技术，同时对谐振峰范围的多个波长通道进行光强度监测，进而获得高信噪比的光声信号，实现了高灵敏度的光声成像。
附图说明
图1表示本发明实施例1中提出的基于时分复用多载波探测光干涉型传感器的光声成像系统结构图；
图2表示本发明实施例2中提出的光干涉型传感器结构示意图；
图3表示本发明实施例2中提出的光干涉型传感器输出光谱的示意图；
图4表示本发明实施例2中提出的信号发生器输出的多载波正弦射频信号频率随采样时间间隔的变化关系图；
图5表示本发明实施例1中提出的示波器接收的时域响应信号原始波形图；
图6表示图5中A处的局部放大图；
图7表示平均处理后得到的光声信号的示意图；
图8表示本发明实施例3中提出的基于实施例1所述光声成像系统实现的光声成像方法的流程图。
其中，1.信号发生器；2.可调激光器；3.IQ调制器；4.偏振控制器；5.光干涉型传感器；6.光电探测器；7.示波器；51.硫系微环；52.直波导。
具体实施方式
附图仅用于示例性说明，不能理解为对本专利的限制；
为了更好地说明本实施例，附图某些部位会有省略、放大或缩小，并不代表实际尺寸；
对于本领域技术人员来说，附图中某些公知内容说明可能省略是可以理解的。
附图中描述位置关系的仅用于示例性说明，不能理解为对本专利的限制；
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
传统单载波探测并解调光谱中的超声脉冲响应信息，从而用于光声成像的方式，一方面，得到的光声信号的信噪比取决于传感元件的灵敏度和系统噪声等，在传感元件以及系统中各种设备固定的情况下，通常很难直接对系统信噪比进行提升；另一方面，探测光强变化的光电探测器和示波器的带宽可达到几GHz甚至几十GHz，远高于光声信号的频率范围，探测系统的时域利用率较低，因此，综合存在时域利用率低、实现成本高且探测光声信号信噪比低的问题，为解决以上问题，本实施例提出一种基于时分复用多载波探测光干涉型传感器的光声成像系统，该系统的结构图如图1所示，参见图1，光声成像系统包括：信号发生器1、可调激光器2、IQ调制器3、偏振控制器4、光干涉型传感器5、光电探测器 6及示波器7；信号发生器1在每个采样周期内，按n个等时间采样间隔依次输出n个频率的多载波正弦波射频信号：f1，…，fn，n个正弦波射频信号输入IQ 调制器3，可调激光器2产生中心频率为f0的探测激光，IQ调制器3基于n个正弦波射频信号对可调激光器2产生的探测激光进行单边带调制，在每个采样周期内产生中心频率为f0+f1、…、f0+fn的信号激光，并输入偏振控制器4进行偏振态调整，通过带隙光聚焦照射，调整光干涉型传感器5输出光谱中同时产生若干个谐振峰，且谐振峰中心频率为f0，调整偏振态后的信号激光作为探针，监测光干涉型传感器5在每个采样周期内的n个等时间采样间隔对生物组织内激发的光声信号(生物组织内激发的光声信号由激光器聚焦激发光源照射生物组织产生，参见图1，基于532nm激光器聚焦激发光源)的时域响应信号；之后，作为探针的信号激光进入光电探测器6，光强度转化为电压强度信号后传输至示波器 7，携带光干涉型传感器5对光声信号的时域响应信号被示波器7接收。
整体上，以光干涉型传感器5作为基本结构，首先利用信号发生器1输出若干正弦波射频信号，利用可调激光器2产生探测激光，IQ调制器3基于正弦波射频信号对可调激光器2产生的探测激光进行单边带调制(“单边带调制”属于IQ 调制器的一种工作方式，本领域熟知，具体的，“单边带调制”对于中心频率为f0的激光和fn的正弦波射频信号，会根据设置偏压产生f0-fn或f0+fn的信号激光，而“双边带调制”则会同时输出f0-fn和f0+fn的激光，这里选用单边)，生成了多路载波，即多路信号激光，有效控制成本，避免多路光源复用与解复用的复杂系统的应用。基于光干涉型传感器5的硫系材料光敏性，实现光干涉型传感器5的输出光谱在10GHz范围内，与IQ调制器3和光电探测器6的带宽限制配合，提高探测系统的时域利用率，而且光干涉型传感器5的输出光谱中同时存在多个谐振峰，通过时分复用的多路载波探测技术，同时对谐振峰范围的多个波长通道进行光强度监测，进而获得高信噪比的光声信号，实现了高灵敏度的光声成像。
在本实施例中，光干涉型传感器5为基于硫系材料的片上光学微腔阵列传感器件，包括n个硫系微环51以及一根将n个硫系微环51串联耦合的直波导52， n个硫系微环51的尺寸相同。参见图2，在本实施例中，光干涉型传感器5的结构采用3个硫系微环51以及一根将3个硫系微环51串联耦合的直波导52，组成片上光学微腔阵列传感器件，根据硫系材料的特性，经过合理设计，可使其在输出光谱中实现10GHz范围内(这样可以避免在光声成像过程中需要的额外芯片电极工艺以及持续的电压驱动)同时存在多个谐振峰，分别由每一个硫系微环 51对输入光进行相位干涉产生。
正常的硫系微环的数量越多，传感信噪比提升越大，本实施例仅以3个硫系微环为例，对系统工作过程作简要说明。在本实施例中，每个硫系微环51的半径为60um，圆心间距为65um，硫系微环51的宽度2um；3个硫系微环51通过一根宽度为1.5um的直波导52串联耦合。由于硫系材料具有光敏特性(即光致折变效应，带隙光照射在材料上引起材料分子结构变化，从而改变材料折射率)，在片上光学微腔阵列传感器件加工完成后，可使用带隙光分别聚焦照射单个硫系微环51，调节其对应产生的谐振峰位置，使得输出光谱中三个谐振峰的中心频率集中在10GHz范围以内，具体如图3中的曲线所示。图3中横坐标表示频率，纵坐标表示光强，3个硫系微环51发生谐振产生谐振峰的点可参见图3，以片上光学微腔阵列传感器件为基本结构，若要监测三个谐振的时域变化，需要产生中心频率分别为f0+f1、f0+f2、f0+f3(这里是指每个谐振的左边带一半功率处的频率) 的激光。
在本实施例中，偏振控制器4调整的偏振态为TE偏振或TM偏振，以配合系统中光传输。
下面基于图2所示的光干涉型传感器结构对系统的工作过程进行说明，结合图1，在系统实际工作时，利用电脑控制信号发生器1在每个采样周期内，按n 个等时间采样间隔依次输出n个频率的多载波正弦波射频信号：f1，…，fn，n 个正弦波射频信号输入IQ调制器3，设每个采样周期时长表示为T，按n个等时间间隔依次输出n个正弦波射频信号时，每一个正弦波射频信号与采样间隔的对应关系为：在(0～1/n)T时间间隔内，输出频率为f1的正弦波射频信号，依次类推，在(n-1/n～1)T时间间隔内，输出频率为fn的正弦波射频信号，针对于图2所示的光干涉型传感器结构时，设每个采样周期T＝5ns，对应采样频率为200MHz 内，0～1/3T时间间隔内输出频率为f1的正弦波射频信号；同样的，在1/3 T～2/3 T和 2/3 T～T的时间内分别输出频率为f2、f3正弦波射频信号，信号发生器1输出的正弦波射频信号频率随时间的变化关系整体如图4所示，输出的正弦波射频信号接入IQ调制器3，将对可调激光器2产生的中心频率为f0的激光进行单边带调制，采样周期内的每个时间间隔内，分别产生中心频率为f0+f1、f0+f2、f0+f3的信号激光多路载波，之后信号激光经过偏振控制器调整所需偏振态(TE或TM偏振)，进入图2所示的光干涉型传感器结构，作为探针分别监测微腔阵列对生物组织内激发的光声信号的时域响应信号，作为探针的信号激光进入光电探测器6，光强度转化为电压强度信号后传输至示波器7，携带光干涉型传感器5对光声信号的时域响应信号被示波器7接收，接收的时域响应信号原始波形图如图5所示，监测光干涉型传感器5在每个采样周期内的3个等时间采样间隔对生物组织内激发的光声信号的时域响应时，提取3个硫系微环51中每一个硫系微环对光声信号的时域响应信号，对图5中A进行局部放大，按照每个采样周期内分别对应的时间段，提取出三个信号激光探针分别的强度变化波形，也即三个硫系微环分别的时域响应信号如图6所示，由于所设计图2所示的微腔阵列的各硫系微环到声源距离的差值远小于信号波长(对于生物光声信号的典型中心频率30MHz，水中声速1500m/s，对应的超声波长为50um)，因此信号传播至各微腔的相位差可忽略不计。
对3个时域响应信号作采样周期对齐，并做平均处理，以降低随机噪声，得到图7中随机噪声更低的光声时域信号，在此，采样周期对齐时，假设原本采样时间为：0 1 2 3；0.33 1.33 2.33 3.33；0.66 1.66 2.66 3.66的三路信号，对齐为：0 1 2 3；0 1 2 3；0 12 3，对齐了才能对相同采样时间的三组电压做数值平均，平均信号的信噪比由原始波形的4.10提升至7.09，提升约1.72倍，噪比提升，说明对信号的采集更加清晰，可传感更微小的超声信号，成像质量更高，且符合N次信号平均对噪声的理论抑制比进一步说明可通过继续增加硫系微环的数量(信号数)，提高信噪比。
使用成熟的二维步进电机控制系统，用532nm聚焦光源逐点扫描生物组织，对于组织内部不同生物结构以及成分激发出不同强度的光声信号。基于本发明实施例提出的时分复用多载波探测结合光干涉型传感器的超声传感，即可将微弱光声信号转化为极高信噪比的时域电信号，从而生成出具有高灵敏度以及高对比度的组织二维图像。
实施例2
在本实施例中，可调激光器2为单一的c波段可调激光器。单一个c波段可调激光器2作为探测光源，通过调制的方式生成多路载波信号激光，有效控制了成本，也避免了多路光源复用与解复用的复杂系统的应用。
实施例3
本实施例根据实施例1所述的基于时分复用多载波探测光干涉型传感器的光声成像系统实现的光声成像方法，该方法的流程图如图8所示，所述方法包括：
S1.确定采样周期，利用信号发生器1在每个采样周期内，按n个等时间采样间隔依次输出n个频率的多载波正弦波射频信号：f1，…，fn，将n个正弦波射频信号输入IQ调制器3；
S2.利用可调激光器2产生中心频率为f0的探测激光，IQ调制器3基于n个正弦波射频信号对可调激光器2产生的探测激光进行单边带调制，在每个采样周期内产生中心频率为f0+f1、…、f0+fn的信号激光，并输入偏振控制器4调整偏振态；
S3.通过带隙光聚焦照射，调整光干涉型传感器5的输出光谱同时产生若干个谐振峰，且谐振峰中心频率为f0；
S4.调整偏振态后的信号激光作为探针，监测光干涉型传感器5在每个采样周期内的n个等时间采样间隔对生物组织内激发的光声信号的时域响应信号
S5.作为探针的信号激光进入光电探测器6，光强度转化为电压强度信号后传输至示波器7，携带光干涉型传感器5对生物组织内激发的光声信号的时域响应信号被示波器7接收。
步骤S5之后还包括：对时域响应信号作采样周期对齐，并做平均处理，以降低随机噪声。
显然，本发明的上述实施例仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。
Claims (9)
1.一种基于时分复用多载波探测光干涉型传感器的光声成像系统，其特征在于，包括：信号发生器（1）、可调激光器（2）、IQ调制器（3）、偏振控制器（4）、光干涉型传感器（5）、光电探测器（6）及示波器（7）；信号发生器（1）在每个采样周期内，按n个等时间采样间隔依次输出n个频率的多载波正弦波射频信号：f1，…，fn，n个正弦波射频信号输入IQ调制器（3），可调激光器（2）产生中心频率为f0的探测激光，IQ调制器（3）基于n个正弦波射频信号对可调激光器（2）产生的探测激光进行单边带调制，在每个采样周期内产生中心频率为f0+f1、…、f0+fn的信号激光，并输入偏振控制器（4）进行偏振态调整，通过带隙光聚焦照射，调整光干涉型传感器（5）输出光谱中同时产生若干个谐振峰，且谐振峰中心频率为f0，调整偏振态后的信号激光作为探针，监测光干涉型传感器（5）在每个采样周期内的n个等时间采样间隔对生物组织内激发的光声信号的时域响应信号；之后，作为探针的信号激光进入光电探测器（6），光强度转化为电压强度信号后传输至示波器（7），携带光干涉型传感器（5）对光声信号的时域响应信号被示波器（7）接收；
光干涉型传感器（5）为基于硫系材料的片上光学微腔阵列传感器件，包括n个硫系微环（51）以及一根将n个硫系微环（51）串联耦合的直波导（52），所述n个硫系微环（51）的尺寸相同，利用带隙光分别聚焦照射每一个硫系微环，调节该硫系微环产生的谐振峰的位置，使得光干涉型传感器（5）的输出光谱中所有谐振峰的中心频率集中于10GHz范围以内。
2.根据权利要求1所述的基于时分复用多载波探测光干涉型传感器的光声成像系统，其特征在于，设每个采样周期时长表示为T，按n个等时间间隔依次输出n个正弦波射频信号时，每一个正弦波射频信号与采样间隔的对应关系为：在(0~1/n)T时间间隔内，输出频率为f1的正弦波射频信号，依次类推，在(n-1/n~1)T时间间隔内，输出频率为fn的正弦波射频信号。
3.根据权利要求1所述的基于时分复用多载波探测光干涉型传感器的光声成像系统，其特征在于，所述可调激光器（2）为单一的c波段可调激光器。
4.根据权利要求1所述的基于时分复用多载波探测光干涉型传感器的光声成像系统，其特征在于，偏振控制器（4）调整的偏振态为TE偏振或TM偏振。
5.根据权利要求1所述的基于时分复用多载波探测光干涉型传感器的光声成像系统，其特征在于，生物组织内激发的光声信号由激光器聚焦激发光源照射生物组织产生。
6.根据权利要求1所述的基于时分复用多载波探测光干涉型传感器的光声成像系统，其特征在于，监测光干涉型传感器（5）在每个采样周期内的n个等时间采样间隔对生物组织内激发的光声信号的时域响应时，提取n个硫系微环（51）中每一个硫系微环对光声信号的时域响应信号。
7.根据权利要求6所述的基于时分复用多载波探测光干涉型传感器的光声成像系统，其特征在于，在提取n个硫系微环（51）中每一个硫系微环对光声信号的时域响应信号后，得到n个时域响应信号，对n个时域响应信号作采样周期对齐，并做数值平均处理，以降低随机噪声。
8.一种根据权利要求1所述的基于时分复用多载波探测光干涉型传感器的光声成像系统实现的光声成像方法，其特征在于，所述方法包括：
S1.确定采样周期，利用信号发生器（1）在每个采样周期内，按n个等时间采样间隔依次输出n个频率的多载波正弦波射频信号：f1，…，fn，将n个正弦波射频信号输入IQ调制器（3）；
S2.利用可调激光器（2）产生中心频率为f0的探测激光，IQ调制器（3）基于n个正弦波射频信号对可调激光器（2）产生的探测激光进行单边带调制，在每个采样周期内产生中心频率为f0+f1、…、f0+fn的信号激光，并输入偏振控制器（4）调整偏振态；
S3.通过带隙光聚焦照射，调整光干涉型传感器（5）的输出光谱同时产生若干个谐振峰，且谐振峰中心频率为f0；
S4.调整偏振态后的信号激光作为探针，监测光干涉型传感器（5）在每个采样周期内的n个等时间采样间隔对生物组织内激发的光声信号的时域响应信号
S5.作为探针的信号激光进入光电探测器（6），光强度转化为电压强度信号后传输至示波器（7），携带光干涉型传感器（5）对生物组织内激发的光声信号的时域响应信号被示波器（7）接收。
9.根据权利要求8所述的光声成像方法，其特征在于，步骤S5之后还包括：对时域响应信号作采样周期对齐，并做平均处理，以降低随机噪声。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111603250.0A CN114485746B (zh) | 2021-12-24 | 2021-12-24 | 基于时分复用多载波探测光干涉型传感器的光声成像系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111603250.0A CN114485746B (zh) | 2021-12-24 | 2021-12-24 | 基于时分复用多载波探测光干涉型传感器的光声成像系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114485746A CN114485746A (zh) | 2022-05-13 |
| CN114485746B true CN114485746B (zh) | 2023-10-31 |
Family
ID=81496282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111603250.0A Active CN114485746B (zh) | 2021-12-24 | 2021-12-24 | 基于时分复用多载波探测光干涉型传感器的光声成像系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114485746B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117338254B (zh) * | 2023-12-06 | 2024-04-16 | 中山大学附属第一医院 | 一种基于气管导管光声成像的肺泡水氧含量监测系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101799608A (zh) * | 2010-04-02 | 2010-08-11 | 上海交通大学 | 基于硅基微环谐振腔的电控宽带光子射频移相器 |
| CN106249354A (zh) * | 2016-09-27 | 2016-12-21 | 华中科技大学 | 一种基于微环与马赫曾德尔干涉结构的微波光子带阻滤波器 |
| CN109507116A (zh) * | 2018-10-15 | 2019-03-22 | 哈尔滨工业大学 | 一种基于聚合物薄膜微腔的光声光谱气体传感装置及方法 |
| CN110082591A (zh) * | 2019-03-29 | 2019-08-02 | 华东师范大学 | 一种大带宽高精度微波频率测量装置及方法 |
| CN113552071A (zh) * | 2021-08-26 | 2021-10-26 | 天津大学 | 一种光声成像系统 |
| CN113703246A (zh) * | 2021-08-30 | 2021-11-26 | 中山大学 | 基于负热光系数硫系微腔的孤子频梳自产生系统及方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7049843B2 (en) * | 2004-03-10 | 2006-05-23 | Tektronix, Inc. | Signal acquisition probing system using a micro-cavity laser capable of sensing DC voltages |
| US20070015992A1 (en) * | 2005-06-30 | 2007-01-18 | General Electric Company | System and method for optoacoustic imaging |
| CN102322880B (zh) * | 2011-08-18 | 2013-06-05 | 天津大学 | 偏振敏感的分布式光频域反射扰动传感装置和解调方法 |
| CN105099570B (zh) * | 2014-05-22 | 2018-08-17 | 复旦大学 | 正交多载波光源及pdm-qpsk信号发射装置 |
| CN110927121B (zh) * | 2019-12-05 | 2022-07-05 | 深圳大学 | 一种基于白光干涉光谱的相位型spr检测装置及方法 |
| CN112902861B (zh) * | 2021-01-26 | 2022-08-02 | 哈尔滨工程大学 | 一种基于超大测量范围pdh传感的应变测量装置 |
-
2021
- 2021-12-24 CN CN202111603250.0A patent/CN114485746B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101799608A (zh) * | 2010-04-02 | 2010-08-11 | 上海交通大学 | 基于硅基微环谐振腔的电控宽带光子射频移相器 |
| CN106249354A (zh) * | 2016-09-27 | 2016-12-21 | 华中科技大学 | 一种基于微环与马赫曾德尔干涉结构的微波光子带阻滤波器 |
| CN109507116A (zh) * | 2018-10-15 | 2019-03-22 | 哈尔滨工业大学 | 一种基于聚合物薄膜微腔的光声光谱气体传感装置及方法 |
| CN110082591A (zh) * | 2019-03-29 | 2019-08-02 | 华东师范大学 | 一种大带宽高精度微波频率测量装置及方法 |
| CN113552071A (zh) * | 2021-08-26 | 2021-10-26 | 天津大学 | 一种光声成像系统 |
| CN113703246A (zh) * | 2021-08-30 | 2021-11-26 | 中山大学 | 基于负热光系数硫系微腔的孤子频梳自产生系统及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114485746A (zh) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111474158B (zh) | 一种二维谱成像系统和二维成像方法 | |
| CN114485746B (zh) | 基于时分复用多载波探测光干涉型传感器的光声成像系统 | |
| Maxwell et al. | Polymer microring resonators for high-frequency ultrasound detection and imaging | |
| WO2016111862A1 (en) | Method and apparatus for photoacoustic tomography using orbital angular momentum(oam) | |
| CN111277338B (zh) | 一种产生宽带混沌激光的装置 | |
| US7657129B2 (en) | Optical focal plane data coupler | |
| WO2005062090A2 (en) | High energy photonic arbitrary waveform source | |
| JP7086225B2 (ja) | 時間差を有するポンプ光とプローブ光の位相コード変調を使用する光ファイバｂｏｃｄａセンサ | |
| JP5489730B2 (ja) | 波長可変光源装置 | |
| WO2014099962A2 (en) | Wavelength stabilization | |
| JP4845145B2 (ja) | 広帯域多波長光源 | |
| CN111637909A (zh) | 一种偏振正交型双泵浦脉冲botda传感装置及方法 | |
| CN111007154A (zh) | 柔性超声换能器、制作方法及全光超声发射与检测方法 | |
| KR20120082473A (ko) | 광 정합 필터 | |
| Fujisue et al. | Demodulation of acoustic signals in fiber Bragg grating ultrasonic sensors using arrayed waveguide gratings | |
| WO2018164895A1 (en) | Optical data modulators with feedback wavelength control | |
| CN114279569B (zh) | 多载波跟踪探测光干涉型传感器光谱的反馈控制方法 | |
| CN111458950B (zh) | 基于石墨烯与xpm作用的空间双相全光开关装置及调制方法 | |
| CN101320188A (zh) | 基于硅基环形谐振腔的光子可调宽带射频移相器 | |
| CN114190974B (zh) | 一种超声成像系统及其成像方法 | |
| CN115374828A (zh) | 一种张量卷积核加速芯片及方法 | |
| CN114777901A (zh) | 干涉式光纤水听器系统及其线性调频方法 | |
| CN115406521B (zh) | 一种基于温控反馈的阵列式超声探测装置 | |
| RU2623710C1 (ru) | Способ определения центральной частоты симметричной оптической структуры (варианты) и устройство для его реализации | |
| KR101466749B1 (ko) | 광원 장치 및 이를 이용한 광 영상 진단 장치 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |