CN113406838A - 实现双光频梳系统重频倍频的方法及系统 - Google Patents
实现双光频梳系统重频倍频的方法及系统 Download PDFInfo
- Publication number
- CN113406838A CN113406838A CN202110635935.7A CN202110635935A CN113406838A CN 113406838 A CN113406838 A CN 113406838A CN 202110635935 A CN202110635935 A CN 202110635935A CN 113406838 A CN113406838 A CN 113406838A
- Authority
- CN
- China
- Prior art keywords
- optical
- frequency
- frequency comb
- double
- interference
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000001228 spectrum Methods 0.000 claims abstract description 51
- 230000003595 spectral effect Effects 0.000 claims abstract description 12
- 230000003068 static effect Effects 0.000 claims abstract description 12
- 239000006185 dispersion Substances 0.000 claims abstract description 11
- 238000012805 post-processing Methods 0.000 claims abstract description 10
- 230000003287 optical effect Effects 0.000 claims description 123
- 210000001520 comb Anatomy 0.000 claims description 28
- 230000009977 dual effect Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 8
- 239000000835 fiber Substances 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 238000000862 absorption spectrum Methods 0.000 claims description 4
- 230000035559 beat frequency Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 15
- 238000012935 Averaging Methods 0.000 abstract description 12
- 230000001427 coherent effect Effects 0.000 abstract description 10
- 230000035945 sensitivity Effects 0.000 abstract description 6
- 230000007774 longterm Effects 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 239000013307 optical fiber Substances 0.000 description 16
- 238000005070 sampling Methods 0.000 description 10
- 230000010287 polarization Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000012952 Resampling Methods 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 238000000411 transmission spectrum Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 235000013405 beer Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002789 length control Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
- G02F1/35—Non-linear optics
- G02F1/353—Frequency conversion, i.e. wherein a light beam is generated with frequency components different from those of the incident light beams
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/02—Details
- G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
- G01J3/0218—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using optical fibers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/12—Generating the spectrum; Monochromators
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/45—Interferometric spectrometry
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
- G02F1/01—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour
- G02F1/21—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour by interference
- G02F1/212—Mach-Zehnder type
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
- G02F1/01—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour
- G02F1/21—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour by interference
- G02F1/225—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour by interference in an optical waveguide structure
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
- G02F1/01—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour
- G02F1/21—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour by interference
- G02F1/225—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour by interference in an optical waveguide structure
- G02F1/2252—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour by interference in an optical waveguide structure in optical fibres
Landscapes
- Physics & Mathematics (AREA)
- Nonlinear Science (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
一种实现双光频梳系统重频倍频的方法及系统，双光频梳包括重复频率不同且拥有重叠的光谱范围的信号光频梳光源和本振光频梳光源，其中，在双光频梳的输出端接入n阶级联的马赫曾德干涉仪，实时采集双光频梳产生的干涉信号，对采集的干涉信号截取出干涉峰，将干涉峰与参考干涉峰进行比较，计算获得由级联的马赫曾德干涉仪引入的静态时间偏移和色散啁啾，通过信号后处理算法消除级联的马赫曾德干涉仪带来的误差，恢复干涉信号的周期性，实现双光频梳光谱系统相位稳定的重频倍频。本发明无需主动闭环或调谐环路即可实现双光频梳系统的有效倍频，能够进行长期相干平均，提高系统的灵敏度和测量速度，可用于提高双光频梳光谱和测距应用的性能。
Description
技术领域
本发明涉及光频梳领域，特别是涉及一种实现双光频梳系统重频倍频的方法及系统。
背景技术
飞秒激光频率梳，简称光频梳，在时域表现为飞秒超快脉冲序列，在频域则表现为宽光谱范围内等距分布的频率纵模。
光频梳诞生于21世纪初，经过20多年的不断发展和探索，其应用领域从早期的光学原子钟拓展到光谱分析、绝对距离测量、成像、系外行星探索等领域。其中，双光频梳光谱技术由于能够实现自动异步光采样的特性，只需单点探测器无需机械扫描部件即可实现高分辨率、高灵敏度、高准确度、宽光谱范围和快速的光谱测量，成为了影响基础和应用物理、分析化学和生物医学领域的一项重要技术，在环境监测、先进制造、国防军工、航空航天和科学研究等领域得到广泛应用。然而，在很多实际光谱测量应用场景下，光频梳的重频远小于光谱探测所需的分辨率，造成测量速度和灵敏度的无谓损失。为了提高双光频梳系统的测量性能，使其能够更适应大规模工业应用推广的需要，需要开发新的具有低复杂度、低成本、易于集成、鲁棒性强、易于实现的双光频梳重频倍频方法。
需要说明的是，在上述背景技术部分公开的信息仅用于对本申请的背景的理解，因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。
发明内容
本发明的主要目的在于弥补上述背景技术存在的不足，提供一种实现双光频梳系统重频倍频的方法及系统。
为实现上述目的，本发明采用以下技术方案：
本发明的第一方面提供一种实现双光频梳系统重频倍频的方法，所述双光频梳包括重复频率不同且拥有重叠的光谱范围的信号光频梳光源和本振光频梳光源，其特征在于，所述方法包括：在信号光频梳光源和本振光频梳光源的输出端接入n阶级联的马赫曾德干涉仪，实时采集信号光频梳和本振光频梳产生的干涉信号，对采集的干涉信号截取出干涉峰，将所述干涉峰与参考干涉信号的干涉峰进行比较，计算获得由级联的马赫曾德干涉仪引入的静态时间偏移和色散啁啾，通过信号后处理算法消除级联的马赫曾德干涉仪带来的误差，恢复干涉信号的周期性，实现双光频梳光谱系统相位稳定的重频倍频。
进一步地：
所述马赫曾德干涉仪通过光纤耦合器、光纤环形器、空间光路、片上光波导中的任一种形式实现。
在所述马赫曾德干涉仪的干涉臂上设置光开关，以调节双光频梳系统的重频倍频系数。
第N阶马赫曾德干涉仪引入的延迟大约设定为：
其中，fr为待插入马赫曾德干涉仪光频梳的重复频率；k为一个任意的整数，n阶马赫曾德干涉仪可将双光频梳干涉信号的重频提高2n倍。
各阶马赫曾德干涉仪同时接入信号光频梳光源和本振光频梳光源的输出端，或者，只接入其中一个光频梳光源的输出端；可选地，每一阶马赫曾德干涉仪的干涉臂数目可大于2。
还包括：在使用马赫曾德干涉仪获得多脉冲模式的两台光频梳之后，将两台光频梳的输出合束，经过待测样品和参考光路，通过光电探测器对待测样品的光学频率响应进行探测，去除光源背景光谱影响后获得物质的吸收光谱信息。
还包括：将信号光频梳和本振光频梳在合光之前分别分出一部分进行探测，获得重频信号，再分别通过锁相环溯源到射频基准上。
还包括：采用两台波长不同但都位于双光频梳重叠光谱范围内的连续激光器，所述连续激光器分别和双光频梳的与自己邻近的纵模做拍，再对拍信号进行混频，获得双光频梳独立纵模之间的拍频，分别提取出双光频梳干涉信号中两根纵模作为误差信号，获得两台光频梳在连续激光器频率处的相对频率噪声；通过数字误差校正，消除相对频率噪声对双光频梳干涉信号产生的畸变干扰。
对实时采集的干涉信号进行如下处理：将实时采集的干涉信号帧进行傅里叶变换，获得相频谱，与参考干涉信号帧的相频谱做差后进行多项式拟合，代表马赫曾德干涉仪引入的重频倍频误差；将实时干涉信号频谱移除拟合得到的相频响应后，反傅里叶变换回时域，得到校正后的干涉信号帧。
本发明的第二方面提供一种实现双光频梳系统重频倍频的系统，包括光频梳光源模块、重频倍频与外差干涉模块以及信号采集与数据处理模块；所述光频梳光源模块包括重复频率不同且拥有重叠的光谱范围的信号光频梳光源和本振光频梳光源；所述重频倍频与外差干涉模块在信号光频梳光源和本振光频梳光源的输出端接入级联的马赫曾德干涉仪，经双光频梳外差干涉之后获得多脉冲干涉信号，实现干涉信号帧数的指数倍增；所述信号采集与数据处理模块对干涉信号进行实时采集，并将其与参考值进行比较，获得各个干涉峰的出现时间和色散啁啾信息，通过信号后处理算法消除级联马赫曾德干涉仪引入的误差，恢复干涉信号的周期性，实现相位稳定的重频倍频。
在本发明的一些实施例中，将n阶级联的马赫曾德干涉仪插入两台光频梳的输出端，由于双光频梳异步光采样原理，单位时间内的干涉峰数目会增加2n倍。n阶马赫曾德干涉仪可将双光频梳干涉信号的重频提高2n倍。对于马赫曾德干涉仪引入的长度偏差和色散啁啾，通过实时采集信号光频梳和本振光频梳产生的干涉信号，对采集的干涉信号截取出干涉峰，将其与任意选定的参考干涉峰进行比较，计算获得马赫曾德干涉仪产生的静态时间偏移和啁啾，通过数字后处理算法补偿这一不理想的相频传递函数，最终获得相位稳定的重频倍频双光频梳光谱系统。
校正后的双光频梳干涉信号恢复了完全的相位稳定性，可进行逐帧的时域相干平均，在相同时间内相比倍频之前帧数增加2n倍，平均后可获得更高信噪比的干涉信号。最终实现的双光频梳光谱系统刷新速率提高2n倍，可用于快速时间分辨光谱测量。
在一些实施例中，各阶马赫曾德干涉仪可同时插入两台光频梳的输出端，也可以只插入其中一个光频梳的输出端，各阶之间的前后次序可任意调节。每一阶马赫曾德干涉仪的干涉臂数目可大于2，各阶干涉臂的数目相乘，能够实现任意整数倍的重频倍频。
在一些实施例中，马赫曾德干涉仪通过1×2和2×2光纤耦合器连接而成，同样也可通过光纤环形器、空间光路、片上光波导等形式实现。在马赫曾德干涉仪的干涉臂上可加入光开关，灵活调节双光频梳系统的重频倍频系数。
在一些实施例中，在使用马赫曾德干涉仪获得多脉冲模式的两台光频梳之后，将两台光频梳的输出合束，经过待测样品和参考光路，使用光电探测器对样品的光学频率响应进行探测，去除光源背景光谱影响后获得物质的吸收光谱信息。
在一些实施例中，两台光频梳光源的输出光分别分出一小部分打在光电探测器上，获得重频信号。通过锁相环将其锁定到射频信号发生器上，从而实现光频梳种子源的重复频率溯源和灵活调谐。锁相环无需高的带宽和锁定精度，只需要保证双光频梳的干涉信号频谱能够无混叠的落在采样定理限制的奈奎斯特频率范围即可。
在一些实施例中，采用数字误差校正方法消除双光频梳之间的相对频率噪声。引入两台自由运行的连续激光器作为光学媒介，每一台连续激光器与两台光频梳的做拍后相互混频，可获得两台光频梳在连续激光器频率处的相对频率噪声。通过数字误差校正，可消除相对频率噪声对双光频梳干涉信号产生的畸变干扰。
在一些实施例中，马赫曾德干涉仪引入的时间延时并不需要特别的精密长度控制，因为残余的光程误差只是会引起时间的静态偏移，可以通过信号后处理算法进行消除。此外，光纤引入的啁啾也无需进行色散补偿，同样可以通过算法进行校正。通过对所述两台光频梳合光产生的干涉信号进行实时采集，对采集的干涉信号逐帧截取出干涉峰，获得实时干涉信号帧，将其与参考信号帧比对，获得马赫曾德干涉仪引入的静态时间偏移和附加啁啾。
在一些实施例中，将实时干涉信号帧进行傅里叶变换，获得相频谱，与参考干涉信号帧的相频谱做差后进行多项式拟合，代表马赫曾德干涉仪引入的重频倍频误差。将实时干涉信号频谱移除拟合得到的相频响应后，反傅里叶变换回时域，即可得到校正后的干涉信号帧，消除静态时间偏移和附加啁啾的影响。
在一些实施例中，在两台光频梳光源的输出端接入级联的保偏光纤马赫曾德干涉仪，调节光频梳光源的偏振态，使其与保偏光纤光轴方向一致，提高光的耦合效率。单个马赫曾德干涉仪通过两个光纤耦合器实现，光频梳通过1×2光纤耦合器之后分成两束，经过不同的传播距离后再次通过耦合器合束。连续通过2×2的耦合器之后即可实现马赫曾德干涉仪的级联。
两台光频梳经过马赫曾德干涉仪之后的光通过2×2耦合器合束，分别经过待测样品和无样品的参考光路，最终获得重频倍频的双光频梳干涉信号。当马赫曾德干涉仪的级联阶数为n时，干涉信号的重复频率值从最初的双光频梳的重频差提高2n倍。
在一些实施例中，为了消除双光频梳相对频率噪声对干涉信号产生的畸变，引入两台自由运行的连续激光器作为光学媒介，分别提取出两台光频梳在不同频率位置的相对纵模拍频信号，希尔伯特变换求取瞬时频率和相位后获得对应干涉信号两根纵模的频率fp1、fp2和相位抖动
δfp2(t)，计算出双光频梳干涉信号的时间抖动
以及载波包络相位抖动
对采集到的双光频梳干涉信号I0(t)，首先消除载波包络相位抖动
再利用时域重采样或插值消除时间抖动
I2(t)＝I1(t-δτ0(t))
I2(t)即为完全校正双光频梳相对频率抖动之后恢复出的干涉信号，此时干涉信号恢复了完全的周期性，只残余激光相对强度噪声、探测器噪声、采样噪声等类型的幅值噪声干扰，可通过相干平均进行信噪比的提高。通过把连续的干涉信号I2(t)按照双光频梳的重频差作为重复周期进行分割，每一帧的信号可表示成
采用如下公式计算m次相干平均后的干涉信号
由于n阶马赫曾德干涉仪的倍频作用，相干平均后的基础帧内包含的干涉峰数量为2n，因此基础帧可再次分割为2n个子帧。由于光纤马赫曾德干涉仪产生的实际延时和理论延时之间的差异，不同的干涉信号子帧之间存在静态时间偏移。此外，干涉仪不同臂引入的光纤色散不同，产生不同的啁啾。为此，将
分割出的2n个子帧
进行傅里叶变换，得到它们的相频谱Pi(f)，选取其中任意一帧的相频谱Pk(f)作为参考，得到其它帧的相对相频谱δPi(f)＝Pi(f)-Pk(f)，对δPi(f)进行多项式拟合得到啁啾和静态时间偏移引入的相频传递函数
采用如下公式进行数字的误差校正
再次相干平均后，可得到最终的干涉信号帧
傅里叶变换后将频率轴从射频域变换到光频域，即可得到双光频梳光源承载的光谱信息。通过比较测量臂和参考臂获得的光谱，即可得到物质的透射光谱和吸光度，实现宽光谱范围内物质组分、浓度、压强、温度等信息的测量。
与现有技术相比，本发明具有如下有益效果：
本发明提供一种实现双光频梳光谱系统重频倍频方法及系统，在两台光频梳的输出端接入级联的马赫曾德干涉仪，双光频梳外差干涉之后获得多脉冲干涉信号，实现干涉信号帧数的指数倍增；信号采集与数据处理过程中，先对干涉信号进行实时采集，获得各个干涉峰的出现时间和色散啁啾信息，并将其与参考值进行比较，通过信号后处理算法消除级联马赫曾德干涉仪引入的误差，恢复干涉信号的周期性，实现相位稳定的重频倍频双光频梳光谱系统。本发明仅使用被动的级联马赫曾德干涉仪和信号后处理算法，无需任何主动闭环或调谐环路即可实现双光频梳系统的有效倍频，能够进行长期相干平均，提高双光频梳系统的灵敏度和测量速度，该方法可用于提高双光频梳光谱和测距应用的性能。
本发明只需使用成熟的光纤通信器件组成级联的马赫曾德干涉仪，结合信号后处理算法能够实现高紧凑度、高鲁棒性、低系统复杂度、高相位稳定性的重频倍频双光频梳光谱系统。相比于基于法布里-珀罗谐振腔的重频倍频方法具有更高的能量利用率，无需闭环反馈控制，易于实现宽光谱范围。随着马赫曾德干涉仪级联阶数的增加，双光频梳干涉信号的获取速度和灵敏度可以以指数倍得到提升，实现更快速的时间分辨光谱测量。对需要高灵敏度的静态光谱测量，能够加速相干平均过程，在相同时间下获得更高的信噪比。为提高光频梳绝对距离测量、光谱分析和精密传感系统的性能提供了一种简单有效的方案。由于马赫曾德干涉仪已经成为了成熟的集成光学组件，本方法能够应用于芯片级集成光频梳的重频倍频，为提升片上光频梳光谱系统性能奠定了基础。
附图说明
图1是本发明一种实施例的双光频梳光谱系统的结构示意图。
图2是本发明一种实施例的信号采集与数据处理的流程图。
图3是采用连续激光器作为媒介消除双光频梳光源噪声的原理示意图。
图4是基于马赫曾德干涉仪的双光频梳异步光采样原理示意图。
具体实施方式
以下对本发明的实施方式作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。
需要说明的是，当元件被称为“固定于”或“设置于”另一个元件，它可以直接在另一个元件上或者间接在该另一个元件上。当一个元件被称为是“连接于”另一个元件，它可以是直接连接到另一个元件或间接连接至该另一个元件上。另外，连接既可以是用于固定作用也可以是用于耦合或连通作用。
需要理解的是，术语“长度”、“宽度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明实施例和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多该特征。在本发明实施例的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
参阅图1和图2，本发明实施例提供一种灵活的重频倍频双光频梳光谱测量方法，包括光频梳光源模块、重频倍频与外差干涉模块以及信号采集与数据处理模块。所述光频梳光源模块包括重复频率不同且拥有重叠的光谱范围的两台光频梳光源。两台光频梳分别通过锁相环将重复频率锁定到频率基准，保证光频梳光源重频fr1和fr2的长期稳定性和可调谐性。双光频梳的重叠光谱范围Δv和重频差Δfr需满足带通采样定理要求。
由于重频锁定只能保证光频梳的长期稳定性，在短的时间尺度上两台光频梳依然处于自由运行状态，优选实施例对双光频梳的相对频率噪声引起的干涉信号的时间抖动和载波包络相位抖动进行实时的提取和补偿。为此，优选实施例中采用两台波长不同但都位于双光频梳重叠光谱范围内的连续激光器。如图3所示，连续激光器分别和两台光频梳与自己邻近的纵模做拍，再对拍信号进行混频，即可获得两台光频梳独立纵模之间的拍频。由于连续激光器只是起到一个中间媒介的作用，其频率噪声不会影响到抖动提取，因此采用自由运行的连续激光器即可满足要求。两台连续激光器能够分别提取出双光频梳干涉信号中两根纵模作为误差信号E1(t)和E2(t)。
为了实现鲁棒、高效、高能量利用率的双光频梳重频倍频，提高干涉信号峰的刷新频率Δfr，可将n阶级联的光纤马赫曾德干涉仪分别插入两台光频梳的输出端。为了避免不同偏振态导致干涉对比度下降，马赫曾德干涉仪全部采用保偏光纤。当光频梳光源的输出采用单模光纤传输时，为了提高的耦合效率，优选实施例在进入干涉仪之前使用偏振控制器调节入射光的偏振态，使其与保偏光纤光轴方向一致。单个马赫曾德干涉仪由两个光纤耦合器组成，光频梳通过1×2光纤耦合器之后分成两束，经过不同的传播距离后再次通过耦合器合束。连续通过2×2的耦合器之后即可实现马赫曾德干涉仪的级联。如图4所示，在光频梳1的输出光经过延时为Δτ1的一阶马赫曾德干涉仪后，与光频梳2合光后，根据异步光采样过程，干涉信号产生的延时Δt为
为了实现干涉信号重频的二倍频，满足Δτ1＝1/2fr1。同理，如果向光频梳2的输出光引入Δτ2的延迟，当Δτ2＝1/2fr2时，也能产生相同的重频倍频效果。因此，为了实现2n倍的重频倍频，将马赫曾德干涉仪进行级联，第N阶马赫曾德干涉仪引入的延迟需大约设定为：
各阶马赫曾德干涉仪可以同时插入两台光频梳的输出端，也可以只插入任意光频梳的输出端，各阶之间的次序可任意调节。当插入光频梳1的输出端时，式中的fr取fr1，反之fr取fr2。各阶马赫曾德干涉仪中均加入光开关，可灵活控制干涉仪的阶数和延迟数目，实现重频倍频系数的灵活调节，使得光谱测量系统能够适应不同应用的需要。
除了使用光纤耦合器，同样也可使用光纤环形器、空间光路、片上光波导等形式实现级联的马赫曾德干涉仪。耦合器的输入输出数目也可以灵活调节，使单阶马赫曾德干涉仪的干涉臂数目为任意整数。此时，级联马赫曾德干涉仪的总体重频倍频系数为各阶干涉臂数目相乘，从而能够实现任意整数倍的重频倍频。
在经过级联的马赫曾德干涉仪之后，两台光频梳的输出脉冲产生了相应的时间延迟，利用2×2耦合器或空间光路合束之后分别经过待测样品和无样品的参考光路，最后入射到探测器，发生异步光采样之后获得多脉冲的干涉信号。其中，测量干涉信号包含了被测样品产生的光谱调制，参考干涉信号既可以通过单独的参考光路获得，也可以通过对测量光谱进行多项式拟合、非线性函数拟合等方式获得双光频梳背景光谱，从而获得样品的光学响应函数。利用低通滤波器滤除干涉信号中高于fr/2的频率成分，进入信号采集与数据处理模块。
图2示出了数据采集、处理与偏频控制的算法流程图，首先，可使用信号采集卡、示波器或FPGA对双光频梳干涉信号I0(t)和误差信号E1(t)和 E2(t)进行ADC采样。在信号调理阶段，对误差信号进行滤波，提取出需要的双光频梳纵模信号，并将其复数化为解析信号，方便后续的校正运算。对干涉信号寻找干涉峰，每2n个截取为一个基础帧，时间长度为双光频梳重频差的倒数1/Δfr。在误差校正与相干平均阶段，首先对双光频梳光源相对频率噪声产生的干涉信号畸变进行误差计算和补偿，对采样到的两个通道误差信号进行希尔伯特变换获得射频纵模频率fp1、fp2和相位抖动
δfp2(t)，采用如下公式可提取出时间抖动δτ0(t)和载波包络相位抖动
对双光频梳干涉信号I0(t)，首先消除载波包络相位抖动
再利用时域重采样或插值消除时间抖动
I2(t)＝I1(t-δτ0(t))
I2(t)即为完全校正双光频梳相对频率抖动之后恢复出的干涉信号，此时干涉信号恢复了完全的周期性，只残余激光相对强度噪声、探测器噪声、采样噪声等类型的幅值噪声干扰，可通过对时间长度为1/Δfr的基础帧
进行相干平均进行信噪比的提高。采用如下公式计算m次相干平均后的干涉信号
n阶级联的马赫曾德干涉仪会在基础帧内产生2n个子帧，但由于光纤马赫曾德干涉仪不同光程的延迟臂引入的啁啾不同，以及实际长度与理论长度之间的偏差导致不同子帧之间的静态时间偏移。
为此，将
分割出的2n个子帧
进行傅里叶变换，得到它们的相频谱Pi(f)，选取其中任意一帧的相频谱Pk(f)作为参考，得到其它帧的相对相频谱δPi(f)＝Pi(f)-Pk(f)，对δPi(f)进行多项式拟合得到啁啾和静态时间偏移引入的相频传递函数
采用如下公式进行数字的误差校正
再次相干平均后，可得到最终的干涉信号帧
对
进行傅里叶变换，可得到干涉信号频谱Sf(f)，除以参考臂干涉信号计算得到的频谱Sref(f)后，消除双光频梳光源背景光谱，使用多项式拟合消除残余基线，获得被测样品的透射光谱T(f)，采用如下公式将频率轴从射频域变换到光频域
从而获得物质的吸收光谱T(υ)。根据朗伯比尔定律
T(υ)＝e-α(v)L
可计算出被测样品的宽带吸光度信息α(v)，通过非线性Voigt线型拟合等分析方法，可获得物质组分、浓度、压强、温度等信息。
本发明的背景部分可以包含关于本发明的问题或环境的背景信息，而不一定是描述现有技术。因此，在背景技术部分中包含的内容并不是申请人对现有技术的承认。
以上内容是结合具体/优选的实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，其还可以对这些已描述的实施方式做出若干替代或变型，而这些替代或变型方式都应当视为属于本发明的保护范围。在本说明书的描述中，参考术语“一种实施例”、“一些实施例”、“优选实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。尽管已经详细描述了本发明的实施例及其优点，但应当理解，在不脱离专利申请的保护范围的情况下，可以在本文中进行各种改变、替换和变更。
Claims (10)
1.一种实现双光频梳系统重频倍频的方法，所述双光频梳包括重复频率不同且拥有重叠的光谱范围的信号光频梳光源和本振光频梳光源，其特征在于，所述方法包括：在信号光频梳光源和本振光频梳光源的输出端接入n阶级联的马赫曾德干涉仪，实时采集信号光频梳和本振光频梳产生的干涉信号，对采集的干涉信号截取出干涉峰，将所述干涉峰与参考干涉信号的干涉峰进行比较，计算获得由级联的马赫曾德干涉仪引入的静态时间偏移和色散啁啾，通过信号后处理算法消除级联的马赫曾德干涉仪带来的误差，恢复干涉信号的周期性，实现双光频梳光谱系统相位稳定的重频倍频。
2.如权利要求1所述的方法，其特征在于，所述马赫曾德干涉仪通过光纤耦合器、光纤环形器、空间光路、片上光波导中的任一种形式实现。
3.如权利要求1所述的方法，其特征在于，在所述马赫曾德干涉仪的干涉臂上设置光开关，以调节双光频梳系统的重频倍频系数。
4.如权利要求1至3任一项所述的方法，其特征在于，第N阶马赫曾德干涉仪引入的延迟大约设定为：
其中，fr为待插入马赫曾德干涉仪光频梳的重复频率；k为一个任意的整数，n阶马赫曾德干涉仪可将双光频梳干涉信号的重频提高2n倍。
5.如权利要求1所述的方法，其特征在于，各阶马赫曾德干涉仪同时接入信号光频梳光源和本振光频梳光源的输出端，或者，只接入其中一个光频梳光源的输出端；可选地，每一阶马赫曾德干涉仪的干涉臂数目可大于2。
6.如权利要求1所述的方法，其特征在于，还包括：在使用马赫曾德干涉仪获得多脉冲模式的两台光频梳之后，将两台光频梳的输出合束，经过待测样品和参考光路，通过光电探测器对待测样品的光学频率响应进行探测，去除光源背景光谱影响后获得物质的吸收光谱信息。
7.如权利要求1所述的方法，其特征在于，还包括：将信号光频梳和本振光频梳在合光之前分别分出一部分进行探测，获得重频信号，再分别通过锁相环溯源到射频基准上。
8.如权利要求1所述的方法，其特征在于，还包括：采用两台波长不同但都位于双光频梳重叠光谱范围内的连续激光器，所述连续激光器分别和双光频梳的与自己邻近的纵模做拍，再对拍信号进行混频，获得双光频梳独立纵模之间的拍频，分别提取出双光频梳干涉信号中两根纵模作为误差信号，获得两台光频梳在连续激光器频率处的相对频率噪声；通过数字误差校正，消除相对频率噪声对双光频梳干涉信号产生的畸变干扰。
9.如权利要求1至8任一项所述的方法，其特征在于，对实时采集的干涉信号进行如下处理：将实时采集的干涉信号帧进行傅里叶变换，获得相频谱，与参考干涉信号帧的相频谱做差后进行多项式拟合，代表马赫曾德干涉仪引入的重频倍频误差；将实时干涉信号频谱移除拟合得到的相频响应后，反傅里叶变换回时域，得到校正后的干涉信号帧。
10.一种实现双光频梳系统重频倍频的系统，其特征在于，包括光频梳光源模块、重频倍频与外差干涉模块以及信号采集与数据处理模块；所述光频梳光源模块包括重复频率不同且拥有重叠的光谱范围的信号光频梳光源和本振光频梳光源；所述重频倍频与外差干涉模块在信号光频梳光源和本振光频梳光源的输出端接入级联的马赫曾德干涉仪，经双光频梳外差干涉之后获得多脉冲干涉信号，实现干涉信号帧数的指数倍增；所述信号采集与数据处理模块对干涉信号进行实时采集，获得各个干涉峰的出现时间和色散啁啾信息，并将其与参考值进行比较，通过信号后处理算法消除级联马赫曾德干涉仪引入的误差，恢复干涉信号的周期性，实现相位稳定的重频倍频。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110635935.7A CN113406838B (zh) | 2021-06-08 | 2021-06-08 | 实现双光频梳系统重频倍频的方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110635935.7A CN113406838B (zh) | 2021-06-08 | 2021-06-08 | 实现双光频梳系统重频倍频的方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113406838A true CN113406838A (zh) | 2021-09-17 |
| CN113406838B CN113406838B (zh) | 2022-09-16 |
Family
ID=77676902
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110635935.7A Active CN113406838B (zh) | 2021-06-08 | 2021-06-08 | 实现双光频梳系统重频倍频的方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113406838B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114337811A (zh) * | 2021-12-24 | 2022-04-12 | 中国人民解放军国防科技大学 | 基于双光频梳和法布里-帕罗滤波器的微波频率测量方法 |
| CN114895502A (zh) * | 2022-04-26 | 2022-08-12 | 东南大学 | 一种基于光开关切换网络的重频大范围可调氮化硅微光梳阵列 |
| CN116930995A (zh) * | 2023-09-19 | 2023-10-24 | 天津大学四川创新研究院 | 调频连续波激光高速目标的速度和距离测量系统及方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080131053A1 (en) * | 2006-08-31 | 2008-06-05 | The Furukawa Electric Co., Ltd. | Arrayed waveguide grating optical multiplexer/demultiplexer |
| CN109060767A (zh) * | 2018-10-19 | 2018-12-21 | 清华大学 | 一种双光频梳光谱聚焦相干反斯托克斯拉曼光谱探测系统 |
| CN110703533A (zh) * | 2019-10-22 | 2020-01-17 | 成都天奥电子股份有限公司 | 一种基于光纤环形干涉仪的光学频率梳重复频率倍频器 |
| CN111952831A (zh) * | 2020-07-29 | 2020-11-17 | 清华大学深圳国际研究生院 | 双光频梳系统及其偏频控制方法 |
| CN112710612A (zh) * | 2021-01-14 | 2021-04-27 | 天津博霆光电技术有限公司 | 一种体积光学相干层析测振系统 |
-
2021
- 2021-06-08 CN CN202110635935.7A patent/CN113406838B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080131053A1 (en) * | 2006-08-31 | 2008-06-05 | The Furukawa Electric Co., Ltd. | Arrayed waveguide grating optical multiplexer/demultiplexer |
| CN109060767A (zh) * | 2018-10-19 | 2018-12-21 | 清华大学 | 一种双光频梳光谱聚焦相干反斯托克斯拉曼光谱探测系统 |
| CN110703533A (zh) * | 2019-10-22 | 2020-01-17 | 成都天奥电子股份有限公司 | 一种基于光纤环形干涉仪的光学频率梳重复频率倍频器 |
| CN111952831A (zh) * | 2020-07-29 | 2020-11-17 | 清华大学深圳国际研究生院 | 双光频梳系统及其偏频控制方法 |
| CN112710612A (zh) * | 2021-01-14 | 2021-04-27 | 天津博霆光电技术有限公司 | 一种体积光学相干层析测振系统 |
Non-Patent Citations (1)
| Title |
|---|
| 张继涛 等: "利用光频梳提高台阶高度测量准确度的方法", 《物理学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114337811A (zh) * | 2021-12-24 | 2022-04-12 | 中国人民解放军国防科技大学 | 基于双光频梳和法布里-帕罗滤波器的微波频率测量方法 |
| CN114337811B (zh) * | 2021-12-24 | 2023-08-29 | 中国人民解放军国防科技大学 | 基于双光频梳和法布里-帕罗滤波器的微波频率测量方法 |
| CN114895502A (zh) * | 2022-04-26 | 2022-08-12 | 东南大学 | 一种基于光开关切换网络的重频大范围可调氮化硅微光梳阵列 |
| CN116930995A (zh) * | 2023-09-19 | 2023-10-24 | 天津大学四川创新研究院 | 调频连续波激光高速目标的速度和距离测量系统及方法 |
| CN116930995B (zh) * | 2023-09-19 | 2023-11-28 | 天津大学四川创新研究院 | 调频连续波激光高速目标的速度和距离测量系统及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113406838B (zh) | 2022-09-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113406838B (zh) | 实现双光频梳系统重频倍频的方法及系统 | |
| CN109357763B (zh) | 一种基于时间分辨光频梳的大气吸收光谱测量系统及方法 | |
| CN103344623B (zh) | 一种提高精度的相干反斯托克斯拉曼散射光梳光谱探测方法 | |
| US8917396B2 (en) | Fourier transform spectrometer with a frequency comb light source | |
| Morozov et al. | Metrological aspects of symmetric double frequency and multi frequency reflectometry for fiber Bragg structures | |
| US10197442B2 (en) | Dual-comb spectroscopy with a free-running bidirectionally mode-locked fiber laser | |
| CN107219002B (zh) | 一种超高分辨率光谱测量方法及系统 | |
| DE112010003904T5 (de) | Optische Signalverarbeitung mit modengekoppelten Lasern | |
| CN107872274B (zh) | 一种光纤色散系数的测量方法 | |
| CA2690910A1 (en) | Referencing of the beating spectra of frequency combs | |
| CN104316186A (zh) | 一种基于光学频率梳的光谱测量方法 | |
| CN112129491B (zh) | 基于单光频梳干涉的光纤时延测量方法及装置 | |
| Hoghooghi et al. | 11-μs time-resolved, continuous dual-comb spectroscopy with spectrally filtered mode-locked frequency combs | |
| CN106197691B (zh) | 一种超短脉冲载波包络相位的探测装置及探测方法 | |
| Soriano-Amat et al. | Common-path dual-comb spectroscopy using a single electro-optic modulator | |
| CN109883348B (zh) | 一种使用伪随机码码分复用的pdh多传感器应变测量装置 | |
| CN215339483U (zh) | 气体分子吸收信号增强系统 | |
| CN110553992A (zh) | 一种红外光谱高速测量系统及方法 | |
| CN113281278A (zh) | 一种快速超高分辨瞬态吸收光谱测量装置及测量方法 | |
| CN108845333B (zh) | 一种抑制振动效应的调频连续波激光测距方法 | |
| CN214893682U (zh) | 一种快速超高分辨瞬态吸收光谱测量装置 | |
| US20230417810A1 (en) | Device for wide-band spectral analysis of a signal of interest | |
| CN115773816A (zh) | 一种可调谐可溯源的光谱定标装置 | |
| CN111397732B (zh) | 非定域波长-时间映射的双光子联合光谱测量方法及系统 | |
| CN110601754B (zh) | 一种基于微波光子下变频的光器件光谱响应测试装置及方法 |
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 |