CN115235618B - 一种超快激光场三维时空分布的单帧测量系统及方法 - Google Patents

一种超快激光场三维时空分布的单帧测量系统及方法 Download PDF

Info

Publication number
CN115235618B
CN115235618B CN202211092884.9A CN202211092884A CN115235618B CN 115235618 B CN115235618 B CN 115235618B CN 202211092884 A CN202211092884 A CN 202211092884A CN 115235618 B CN115235618 B CN 115235618B
Authority
CN
China
Prior art keywords
light
measured
laser
frequency spectrum
frequency
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
Application number
CN202211092884.9A
Other languages
English (en)
Other versions
CN115235618A (zh
Inventor
李伟
王逍
曾小明
母杰
胡必龙
左言磊
吴朝辉
王晓东
李钊历
粟敬钦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Laser Fusion Research Center China Academy of Engineering Physics
Original Assignee
Laser Fusion Research Center China Academy of Engineering Physics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Laser Fusion Research Center China Academy of Engineering Physics filed Critical Laser Fusion Research Center China Academy of Engineering Physics
Priority to CN202211092884.9A priority Critical patent/CN115235618B/zh
Publication of CN115235618A publication Critical patent/CN115235618A/zh
Application granted granted Critical
Publication of CN115235618B publication Critical patent/CN115235618B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G01J1/00—Photometry, e.g. photographic exposure meter
    • G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
    • G01J1/4257—Photometry, e.g. photographic exposure meter using electric radiation detectors applied to monitoring the characteristics of a beam, e.g. laser beam, headlamp beam
    • 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
    • G01J1/00—Photometry, e.g. photographic exposure meter
    • G01J1/02—Details
    • G01J1/04—Optical or mechanical part supplementary adjustable parts
    • G01J1/0407—Optical elements not provided otherwise, e.g. manifolds, windows, holograms, gratings
    • G01J1/0411—Optical elements not provided otherwise, e.g. manifolds, windows, holograms, gratings using focussing or collimating elements, i.e. lenses or mirrors; Aberration correction
    • 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/0208—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using focussing or collimating elements, e.g. lenses or mirrors; performing aberration correction
    • 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
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/10—Beam splitting or combining systems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/10—Image enhancement or restoration using non-spatial domain filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • G06T7/85—Stereo camera calibration
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Spectrometry And Color Measurement (AREA)

Abstract

本发明公开了一种超快激光场三维时空分布的单帧测量系统及方法，该方法利用F‑P标准具对超快激光的连续频谱进行等间隔频谱采样，利用光栅将不同频谱成分向不同空间方向衍射，使参考光和待测光的不同频谱成分形成不同调制频率的波前干涉图，利用相机采集频谱混叠的波前干涉图，最后在单帧测量中同时提取多个频谱成分的波前干涉信息，进行光场三维分布重建，实现了超快激光三维光场分布的单帧测量。本发明能够实现超快激光场三维时空分布的单帧波前测量，该方法原理简单可靠，定标和测量操作简便，抗干扰能力强，光场重建算法计算量小，特别适用于以低重频或单发次模式运行的超快超强激光三维光场分布的测量。

Description

一种超快激光场三维时空分布的单帧测量系统及方法
技术领域
本发明属于超快激光场测量领域，尤其涉及一种超快激光场三维时空分布的单帧测量系统及方法。
背景技术
超短超强激光在激光加速、激光聚变、二次源产生、实验室天体物理等领域具有重要的应用价值。世界各地的研究者们都致力于实现超短超强激光更高的远场峰值功率密度，以研究全新的强场物理现象。然而，超短超强激光产生系统由于光路调节误差、材料色散和元件面型等缺陷的存在，导致激光脉冲的时域波形随空间位置而改变，这种现象被称为时空耦合畸变，时空耦合畸变是限制远场功率密度提升的主要因素。准确测量超短超强激光的时空耦合畸变（即光场三维时空分布）具有重要的意义：一方面三维光场分布的精确测量是改善并补偿时空耦合畸变的前提；另一方面可以根据三维光场分布的测量结果预测远场光场分布，为强场物理实验的理论研究提供更准确和完备的激光参数。
然而，由于超快激光场三维时空分布包含了两个空间维度以及光谱/时间维度的振幅和相位，想要在单帧测量中获取所有这些信息十分困难，因此现阶段依然缺乏成熟可靠的超快激光三维光场分布的单帧测量手段。目前的大多数超快激光脉冲测量技术，比如频率分辨光学门（Frequency Resolved Optical Grating, FROG），谱相位相干直接电场重建（Spectral Phase Interferometry for Direct Electric-field Reconstruction,SPIDER）以及自参考光谱干涉（Self-Reference Spectral Interferometry， SRSI），利用这些技术进行超快激光脉冲单帧测量仅能得到激光束横截面内的单个空间点或平均时域特征，这些信息无法用来表征光场的三维时空分布；而且现有的大多数超快激光三维光场测量方法需要对空间域、频率域或时间域进行多次扫描，因此其无法实现单帧测量，而且并不适用于低重频的或各发次间有差异的超快超强激光场表征。也有学者提出了一些三维光场分布的单帧测量方法，如Pablo Gabolde和Rick Trebino等提出的波分复用数字全息法（Gabolde P , Trebino R . Single-frame measurement of the completespatiotemporal intensity and phase of ultrashort laser pulses usingwavelength-multiplexed digital holography[J]. Journal of the Optical Societyof America B, 2008, 25(25):A25-A33.），该方法用二维光栅和带通滤波器将各频谱成分的全息图成像在相机感光面不同空间位置上，实现了三维光场的单帧测量，但是该方法相机像素利用率不高，对其中光学元件位姿调节精度要求高，并且需要复杂的标定过程；刘军等人也提出了一种超快激光完全时间空间耦合特性单发测量方法（刘军、申雄、王鹏。飞秒激光脉冲完全时间空间耦合特性单发测量系统及测量方法，CN111879421A[P]. 2020.），该方法用光纤对激光束横截面进行阵列化取样，再将光纤重新排列在成像光谱仪入口狭缝处，最后利用传统的谱干涉方法恢复每根光纤采样区光场的时域特征，受光谱仪入口狭缝空间尺寸限制，该方法的空间分辨能力无法做到很高。
发明内容
有鉴于此，本发明提供一种超快激光场三维时空分布的单帧测量系统及方法，该方法主要针对单发次激光进行波前测量。该方法首先利用F-P标准具对超快激光的连续频谱进行等间隔采样，再用光栅将不同频谱成分向不同空间方向衍射，使参考光和待测光的不同频谱成分形成不同调制频率的波前干涉图，利用CCD相机采集频谱混叠的波前干涉图，最后再结合频分复用的思想将不同频谱成分的波前干涉图从混叠的波前干涉图中解调出来，从而在单帧测量中可同时提取多个频谱成分的波前干涉信息，实现了超快激光三维光场分布的单帧波前测量。
为达此目的，本发明采用以下技术方案：一种超快激光场三维时空分布的单帧测量系统，所述系统包括：反射取样镜、时域测量设备、F-P标准具、分光棱镜Ⅰ、空间滤波器、延迟调节单元、光栅、分光棱镜Ⅱ、消色差透镜和相机；
所述反射取样镜将入射的待测激光分为透射光和反射光；
所述时域测量设备设置在反射取样镜的反射光光路上；
所述F-P标准具和分光棱镜Ⅰ按顺序依次设置在反射取样镜的透射光光路上；
所述F-P标准具对透射光进行频谱等间隔取样；
所述分光棱镜Ⅰ将F-P标准具的取样光分为透射光和反射光，其中分光棱镜Ⅰ的透射光命名为待测光，分光棱镜Ⅰ的反射光命名为参考光；
所述空间滤波器设置在参考光光路上；
所述延迟调节单元和光栅按顺序依次设置在待测光光路上；
所述分光棱镜Ⅱ设置在参考光和待测光交汇处，用于合成参考光和待测光；
所述消色差透镜和相机依次设置在分光棱镜Ⅱ之后的光传输光路上。
优选的，所述时域测量设备为SPIDER、FROG或者SRSI。
优选的，所述空间滤波器包括沿光轴依次设置的离轴抛面镜Ⅰ、针孔和离轴抛面镜Ⅱ。
优选的，所述延迟调节单元包括多个沿光传输方向设置的反射镜。
一种超快激光场三维时空分布的单帧测量方法，所述方法包括：
S1：对待测激光光谱以及参考光光强分布进行标定，得到参考光各频谱成分的光强分布；
S2：在上述超快激光场三维时空分布的单帧测量系统中输入待测激光，利用时域测量设备测量获取入射待测激光的谱相位，利用相机获得待测激光的待测光和参考光的频谱混叠的波前干涉图；
S3：使用频分复用方法从频谱混叠波前干涉图中恢复出待测激光各频谱成分的光强和波前分布，并结合S2得到的谱相位进行三维光场重建。
优选的，所述S1中，对待测激光光谱以及参考光强分布进行标定的过程如下：首先 在F-P标准具后测量待测激光脉冲光谱，获得透过F-P标准具的激光频率
Figure 952972DEST_PATH_IMAGE001
以及各频率成 分的光强度比；接着利用相机测量经过空间滤波器的参考光的光强分布，并根据光谱标定 得到参考光各频率成分的光强度比，计算得到参考光各频谱成分的光强分布。
优选的，所述S2包括：
S21：输入待测激光，利用反射取样镜对输入待测激光进行分光，其中反射取样镜的反射光进入时域测量设备进行谱相位测量，反射取样镜的透射光传输至F-P标准具；
S22：F-P标准具对反射取样镜的透射光进行频谱等间隔取样，F-P标准具的取样光传输至分光棱镜Ⅰ；
S23：分光棱镜Ⅰ将F-P标准具的取样光进一步分为透射光和反射光，将分光棱镜Ⅰ的透射光命名为待测光，将分光棱镜Ⅰ的反射光命名为参考光；
S24：参考光经过空间滤波器进行波前畸变过滤和空间强度调制以及扩束；待测光经过延迟调节单元与参考光同步，再经光栅，将不同频率成分的激光向不同角度衍射；
S25：经过S24处理的待测光和参考光经分光棱镜Ⅱ合束，再通过消色差透镜将频谱混叠的波前干涉图成像在相机的感光平面上。
优选的，所述S3包括：
S31：对S2中得到的频谱混叠的波前干涉图进行二维傅里叶变换，得到频域图；
S32：计算得到频域图中各频谱对应一级谐波的位置，接着在频域图中使用圆形窗口滤波得到单独的一级谐波，将滤波得到的一级谐波进行傅里叶逆变换，得到该频谱成分的空间域二维复数矩阵，所述二维复数矩阵为一级谐波对应频谱成分的子干涉图，依次对频域图中每个一级谐波都进行傅里叶逆变换，得到每个频谱成分对应的空间域二维复数矩阵；
S33：利用子干涉图复数矩阵以及S1中标定得到参考光各频谱成分的光强分布，计算得到待测光各频谱成分的波前和光强分布；
S34：利用时域测量设备测量得到的谱相位对S33中待测光各频谱成分的波前和光 强分布进行拼接，再沿代表频谱
Figure 512129DEST_PATH_IMAGE002
的坐标进行傅里叶逆变换，得到待测光场三维时空分 布。
本发明的有益效果是：本发明提供的超快激光场三维时空分布的单帧测量系统及方法，该系统及该方法采用频分复用+波前干涉的方法，在单帧测量中可同时提取多个频谱成分的波前干涉信息，将传统波前干涉方法获取信息的能力从二维（两个空间维度）提升到三维（两个空间维度+频谱维度），基于“频分复用+波前干涉”，可以利用多个波长成分波前干涉信息实现超快激光场三维时空分布的单帧测量，本发明提出的基于超快激光场三维时空分布的单帧测量方法是一种全新的测量方法，该方法实现了宽带光、多频谱的激光单发次波前测量，该方法原理简单可靠，定标和测量操作简便，抗干扰能力强，光场重建算法计算量小，特别适用于以低重频或单发次模式运行的超快超强激光三维光场分布的测量。
附图说明
图1为本发明的超快激光场三维时空分布的单帧测量系统的结构示意图；
图2为实施例中对频谱混叠的波前干涉图进行三维光场重建的过程示意图；
图中：1.F-P标准具 2.分光棱镜Ⅰ 3.延迟调节单元 4.光栅 5.分光棱镜Ⅱ 6.消色差透镜 7.相机 8.离轴抛面镜Ⅰ 9.离轴抛面镜Ⅱ 10.针孔 11.反射取样镜 12.时域测量设备 13.成像物面 14.成像像面。
具体实施方式
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。
下面结合附图和具体实施例对本发明进行详细说明。
一种如图1所示的超快激光场三维时空分布的单帧测量系统，该述系统包括：反射取样镜11、时域测量设备12、F-P标准具1、分光棱镜Ⅰ2、空间滤波器、延迟调节单元3、光栅4、分光棱镜Ⅱ5、消色差透镜6和相机7；
其中反射取样镜11用于将入射待测激光分为透射光和反射光，经过反射取样镜11反射的取样光进入设置在反射光光路上的时域测量设备12进行谱相位测量，进而获得待测激光的谱相位信息，作为实施例，时域测量设备12可以是SPIDER、FROG或者SRSI（如Wizzler）等测量装置；
经过反射取样镜11透射的透射光用于待测激光各频谱的波前和光强分布测量。设 置在透射光光路上的F-P标准具1对待测激光进行频谱等间隔取样，F-P标准具1在垂直入射 情况下的取样间隔为
Figure 428132DEST_PATH_IMAGE003
，单模线宽为
Figure 24330DEST_PATH_IMAGE004
，式中，c为真空中光速，n 为F-P标准具1的腔内折射率，h为F-P标准具1的腔长，R为F-P标准具1的腔镜反射率，这些参 数都可以根据实际需求选取。
F-P标准具1的取样光，经过设置在光路上的分光棱镜Ⅰ2后分为透射光和反射光，其中将分光棱镜Ⅰ2的透射光命名为待测光，将分光棱镜Ⅰ2的反射光命名为参考光；参考光经过由离轴抛面镜Ⅰ8、针孔10以及离轴抛面镜Ⅱ9构成的空间滤波器，过滤掉参考光的波前畸变与空间强度调制，同时空间滤波器也具有扩束作用，使参考光与待测光空间重叠的部分具有足够平坦的波前和均匀的光强分布。
待测光经过延迟调节单元3进行延时使其与参考光同步，再经光栅4，使不同频率成分向不同角度衍射；
最后待测光和参考光经设置在两光路交汇处的分光棱镜Ⅱ5合束，再通过消色差透镜6将频谱混叠的波前干涉图成像在相机7的感光平面上，成像物面13和成像像面14分别位于图中光栅4所处位置的虚线处和相机7的感光平面所处位置处。
一种超快激光场三维时空分布的单帧测量方法，该方法基于本发明公开的超快激光场三维时空分布的单帧测量系统进行，具体测量过程如下：
第一步：对待测激光光谱以及参考光强分布进行标定，得到参考光各频谱成分的 光强分布，具体为：首先在F-P标准具后测量待测激光脉冲的光谱，获得透过F-P标准具的激 光频率
Figure 455311DEST_PATH_IMAGE001
以及各频率成分的光强度比；接着利用相机测量参考光光强分布
Figure 942793DEST_PATH_IMAGE005
，并 根据光谱标定得到的各激光频率成分的光强度比，计算得到参考光各频谱成分的光强分布
Figure 775620DEST_PATH_IMAGE006
；
第二步：在上述的超快激光场三维时空分布的单帧测量系统中输入待测激光，根据时域测量设备的测量结果获取入射待测激光谱相位，并用相机测量得到待测光和参考光的频谱混叠的波前干涉图，具体如下：
（1）输入待测激光，利用反射取样镜对待测激光进行取样，其中反射取样镜的反射光进入时域测量设备进行待测激光谱相位测量，反射取样镜的透射光传输至F-P标准具；
（2）F-P标准具对反射取样镜的透射光进行频谱等间隔取样，F-P标准具的取样光传输至分光棱镜Ⅰ；
（3）分光棱镜Ⅰ将F-P标准具的取样光进一步分为透射光和反射光，将分光棱镜Ⅰ的透射光命名为待测光，将分光棱镜Ⅰ的反射光命名为参考光；
（4）参考光经过空间滤波器进行波前畸变过滤和空间强度调制以及扩束；待测光经过延迟调节单元与参考光同步，再经光栅，将不同频率成分向不同角度衍射；
（5）经过上述处理的待测光和参考光经分光棱镜Ⅱ合束，再通过消色差透镜将频谱混叠的波前干涉图成像在相机的感光平面上；
建立如图1所示的坐标系，以参考光经过分光棱镜Ⅱ5反射后的传输方向为z轴，相机7感光平面为x-y平面，待测光与参考光在相机7的感光平面上形成的频谱混叠的波前干涉图可以用如下公式描述：
Figure 542719DEST_PATH_IMAGE007
式中，N表示通过F-P标准具1取样的总的频谱数，
Figure 664259DEST_PATH_IMAGE008
表示通过F-P标准具1取样的 第i个频谱成分的电场角频率，
Figure 565219DEST_PATH_IMAGE009
和
Figure 377186DEST_PATH_IMAGE010
分别表示第i个频谱成分的 待测测光与参考光的强度，
Figure 705399DEST_PATH_IMAGE011
是第i个频谱成分的波数，
Figure 720759DEST_PATH_IMAGE012
和
Figure 425410DEST_PATH_IMAGE013
分别表示待测脉 冲第i个频谱成分的波矢
Figure 91884DEST_PATH_IMAGE014
与x轴和y轴方向的夹角；
上式公式描述了激光脉冲中N个频谱成分的参考光与待测光形成的波前干涉图的线性叠加，由于光栅的存在，使得不同频谱成分的待测光和参考光波前夹角不同，从而形成的波前干涉图的空间调制频率不同，因此可以通过对空间频率域滤波单独提取每一个频谱成分的波前干涉图，进而得到各频谱成分的波前和强度的空间分布信息。
第三步：对待测激光脉冲进行三维光场重建，具体过程如图2所示，包括以下步骤：
首先对第二步中得到的频谱混叠的波前干涉图进行二维傅里叶变换，得到频域图，频域图中心处亮斑为0级谐波，包含干涉图的直流成分，上方一系列亮斑和下方一系列亮斑分别为各频谱成分对应的1级谐波和-1级谐波，待测光的波前和强度的空间分布信息包含在各频谱成分对应的1级谐波和-1级谐波之中，频域图中各频谱对应一级谐波的位置如下：
Figure 794261DEST_PATH_IMAGE015
其中，
Figure 421551DEST_PATH_IMAGE016
，分别为x和y方向的空间频率坐标，-1级谐波位置为对 上式取负号即可；
然后对频域图进行圆形窗口滤波，选取频域图中每个一级谐波，选取的一级谐波 经傅里叶逆变换，得到空间域的二维复数矩阵，以
Figure 805259DEST_PATH_IMAGE017
表示，待测光各频谱成分的 光强分布和波前信息包含在二维复数矩阵的模和辐角中：
Figure 936026DEST_PATH_IMAGE018
Figure 996255DEST_PATH_IMAGE019
其中，
Figure 110841DEST_PATH_IMAGE020
为待测激光的第i个频谱成分的波前。参考光光强分布
Figure 94978DEST_PATH_IMAGE006
可通过预先标定已经确定，在已知参考光的光强分布
Figure 486776DEST_PATH_IMAGE006
，可以根据 模公式求得待测激光光强分布
Figure 327693DEST_PATH_IMAGE021
；在已知
Figure 319789DEST_PATH_IMAGE022
和
Figure 904354DEST_PATH_IMAGE023
的情况下，可以根据辐角公 式得到待测激光波前
Figure 150659DEST_PATH_IMAGE024
，此处假定参考光具有平坦的波前，测量光路中离轴抛 面镜Ⅰ8、离轴抛面镜Ⅱ9以及针孔10构成的空间滤波器结构可以使参考光满足这一假设，
Figure 162477DEST_PATH_IMAGE022
和
Figure 189339DEST_PATH_IMAGE023
可以通过公式
Figure 702229DEST_PATH_IMAGE025
以及1级谐波在频域图中 的位置确定。
最后，利用时域测量设备12得到的谱相位以及计算得到的待测光各频谱成分的波 前和光强分布，沿代表频谱
Figure 193253DEST_PATH_IMAGE026
的z坐标进行傅里叶逆变换，得到待测光场三维时空分布
Figure 985759DEST_PATH_IMAGE027
：
Figure 827813DEST_PATH_IMAGE028
式中，
Figure 878815DEST_PATH_IMAGE029
为待测激光光强，
Figure 958766DEST_PATH_IMAGE030
为待测激光波前，
Figure 781229DEST_PATH_IMAGE031
为 待测脉冲的绝对谱相位，其中不包含随空间位置变化的信息，由时域测量设备12测量得到。
实施例
对中心波长800nm，带宽200nm，光束口径8mm，单脉冲能量10
Figure 454787DEST_PATH_IMAGE032
的超快激光场时空 分布进行测量：
按照图1所示光路搭建测量装置，其中F-P标准具1腔长为h=44
Figure 184845DEST_PATH_IMAGE033
，腔镜反射率R= 95%，在700nm-900nm范围内可取样29个频谱成分，每个频谱成分线宽约为0.5nm；光栅4刻线 密度600线/mm，待测脉冲在光栅上的入射角为13.7°，光栅前待测脉冲传输方向位于y-z平 面内且与z轴之间夹角为1.5°，收集经过光栅的-1级衍射光用于波前干涉；消色差透镜6将 光栅处的光场1:1成像在相机7的感光平面上，相机分辨率为2048*2048，像素间距为5.5μm； 离轴抛面镜8和9的焦距分别为200mm和500mm，滤波小孔10的直径为20
Figure 243937DEST_PATH_IMAGE033
；时域测量设备 12采用商用Wizzler仪器，用于测量激光脉冲的时域特性或绝对谱相位；分光镜11的反射/ 透射比为9:1，将大部分脉冲能量反射，用于时域特性或绝对谱相位测量。
（1）对光谱和参考光强度分布进行标定。在F-P标准具后测量激光脉冲光谱，获得 透过标准具的激光频率
Figure 299618DEST_PATH_IMAGE008
以及各频率成分的强度比；用相机测量参考光光强分布
Figure 460472DEST_PATH_IMAGE034
，并根据光谱标定得到的各频率成分的强度比，计算参考光各频谱成分的光强分 布
Figure 728642DEST_PATH_IMAGE035
；
（2）输入激光脉冲，用相机7采集频谱混叠的波前干涉图，同时用时域测量设备12 测量激光脉冲的绝对谱相位
Figure 986448DEST_PATH_IMAGE036
；
（3）对（2）采集到的频谱混叠的波前干涉图执行二维傅里叶变换，得到频域图；
（4）识别频域图中各频谱对应一级谐波的位置，根据公式计算
Figure 72085DEST_PATH_IMAGE022
和
Figure 579289DEST_PATH_IMAGE023
；
（5）在频域图中使用圆形窗口滤波，单独选取每一个频谱成分对应的一级谐波，并 对其执行二维傅里叶逆变换，得到一系列二维复数矩阵
Figure 792096DEST_PATH_IMAGE037
，即子干涉图；
（6）结合（5）得到的二维复数矩阵
Figure 966726DEST_PATH_IMAGE037
，根据标定得到的参考光的光强分 布
Figure 957684DEST_PATH_IMAGE035
计算待测激光的光强分布
Figure 421027DEST_PATH_IMAGE038
，并根据（4）得到的
Figure 30999DEST_PATH_IMAGE022
和
Figure 935502DEST_PATH_IMAGE023
计算 待测激光波前
Figure 972728DEST_PATH_IMAGE039
。
（7）根据得到的绝对谱相位
Figure 110317DEST_PATH_IMAGE040
，和待测激光的光强分布
Figure 523981DEST_PATH_IMAGE038
和波前
Figure 282989DEST_PATH_IMAGE039
，计算激光脉冲的三维时空场分布
Figure 959958DEST_PATH_IMAGE041
。

Claims (7)

1.一种超快激光场三维时空分布的单帧测量系统，其特征在于，所述系统包括：反射取样镜、时域测量设备、F-P标准具、分光棱镜Ⅰ、空间滤波器、延迟调节单元、光栅、分光棱镜Ⅱ、消色差透镜和相机；
所述反射取样镜将入射的待测激光分为透射光和反射光；
所述时域测量设备设置在反射取样镜的反射光光路上；
所述F-P标准具和分光棱镜Ⅰ按顺序依次设置在反射取样镜的透射光光路上；
所述F-P标准具对透射光进行频谱等间隔取样；
所述分光棱镜Ⅰ将F-P标准具的取样光分为透射光和反射光，其中分光棱镜Ⅰ的透射光命名为待测光，分光棱镜Ⅰ的反射光命名为参考光；
所述空间滤波器设置在参考光光路上；
所述延迟调节单元和光栅按顺序依次设置在待测光光路上；
所述分光棱镜Ⅱ设置在参考光和待测光交汇处，用于合成参考光和待测光；
所述消色差透镜和相机依次设置在分光棱镜Ⅱ之后的光传输光路上。
2.根据权利要求1所述的超快激光场三维时空分布的单帧测量系统，其特征在于，所述时域测量设备为SPIDER、FROG或者SRSI。
3.根据权利要求1所述的超快激光场三维时空分布的单帧测量系统，其特征在于，所述空间滤波器包括沿光轴依次设置的离轴抛面镜Ⅰ、针孔和离轴抛面镜Ⅱ。
4.根据权利要求1所述的超快激光场三维时空分布的单帧测量系统，其特征在于，所述延迟调节单元包括多个沿光传输方向设置的反射镜。
5.一种超快激光场三维时空分布的单帧测量方法，其特征在于，所述方法包括：
S1：对待测激光光谱以及参考光强度分布进行标定，得到参考光各频谱成分的光强度分布；
S2：在权利要求1~4中任一项所述的超快激光场三维时空分布的单帧测量系统中输入待测激光，利用时域测量设备测量获取待测激光的谱相位，利用相机获得待测激光的待测光和参考光的频谱混叠波前干涉图；
S3：使用频分复用方法从频谱混叠波前干涉图中恢复出待测激光各频谱成分的强度和波前分布，并结合S2得到的谱相位进行三维光场重建，其过程如下：
S31：对S2中得到的频谱混叠的波前干涉图进行二维傅里叶变换，得到频域图；
S32：计算得到频域图中各频谱对应一级谐波的位置，接着在频域图中使用圆形窗口滤波得到单独的一级谐波，将滤波得到的一级谐波进行傅里叶逆变换，得到该频谱成分的空间域二维复数矩阵，所述二维复数矩阵为一级谐波对应频谱成分的子干涉图，依次对频域图中每个一级谐波都进行傅里叶逆变换，得到每个频谱成分对应的空间域二维复数矩阵；
S33：利用子干涉图复数矩阵以及S1中标定得到参考光各频谱成分的强度分布，计算得到待测光各频谱成分的波前和光强分布；
S34：利用时域测量设备测量得到的谱相位对S33中待测光各频谱成分的波前和光强分布进行拼接，再沿代表频谱ω的坐标进行傅里叶逆变换，得到待测光场三维时空分布。
6.根据权利要求5所述的超快激光场三维时空分布的单帧测量方法，其特征在于，所述S1中，对待测激光光谱以及参考光强度分布进行标定的过程如下：首先在F-P标准具后测量待测激光脉冲光谱，获得透过F-P标准具的激光频率ω i 以及各频率成分的强度比；接着利用相机测量经过空间滤波器的参考光强度空间分布，并根据光谱标定得到参考光各频率成分的光强度比，计算得到参考光各频谱成分的光强分布。
7.根据权利要求5所述的超快激光场三维时空分布的单帧测量方法，其特征在于，所述S2包括：
S21：输入待测激光，利用反射取样镜对输入待测激光进行分光，其中反射取样镜的反射光进入时域测量设备进行谱相位测量，反射取样镜的透射光传输至F-P标准具；
S22：F-P标准具对反射取样镜的透射光进行频谱等间隔取样，F-P标准具的取样光传输至分光棱镜Ⅰ；
S23：分光棱镜Ⅰ将F-P标准具的取样光进一步分为透射光和反射光，将分光棱镜Ⅰ的透射光命名为待测光，将分光棱镜Ⅰ的反射光命名为参考光；
S24：参考光经过空间滤波器进行波前畸变过滤和空间强度调制以及扩束；待测光经过延迟调节单元与参考光同步，再经光栅，将不同频率成分的激光向不同角度衍射；
S25：经过S24处理的待测光和参考光经分光棱镜Ⅱ合束，再通过消色差透镜将频谱混叠的波前干涉图成像在相机的感光平面上。
CN202211092884.9A 2022-09-08 2022-09-08 一种超快激光场三维时空分布的单帧测量系统及方法 Active CN115235618B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211092884.9A CN115235618B (zh) 2022-09-08 2022-09-08 一种超快激光场三维时空分布的单帧测量系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211092884.9A CN115235618B (zh) 2022-09-08 2022-09-08 一种超快激光场三维时空分布的单帧测量系统及方法

Publications (2)

Publication Number Publication Date
CN115235618A CN115235618A (zh) 2022-10-25
CN115235618B true CN115235618B (zh) 2022-12-13

Family

ID=83682282

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211092884.9A Active CN115235618B (zh) 2022-09-08 2022-09-08 一种超快激光场三维时空分布的单帧测量系统及方法

Country Status (1)

Country Link
CN (1) CN115235618B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116448243B (zh) * 2023-06-19 2023-09-22 中国工程物理研究院激光聚变研究中心 一种基于交叉偏振波的三维光场自参考测量装置及方法
CN116907804B (zh) * 2023-07-24 2024-03-22 中国工程物理研究院激光聚变研究中心 一种单脉冲超快激光时空耦合畸变测量仪

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6639683B1 (en) * 2000-10-17 2003-10-28 Remy Tumbar Interferometric sensor and method to detect optical fields
CN106989834A (zh) * 2017-03-28 2017-07-28 中国工程物理研究院激光聚变研究中心 一种能同时诊断超短脉冲激光的啁啾特性与时空分布特性的方法
CN109100029A (zh) * 2018-09-05 2018-12-28 中国科学院上海光学精密机械研究所 飞秒激光脉冲时空参数单发测量装置
CN109411312A (zh) * 2018-10-31 2019-03-01 中国工程物理研究院流体物理研究所 基于飞秒激光调制的超快电子枪及其验证方法
CN111536868A (zh) * 2020-04-07 2020-08-14 华东师范大学 一种压缩超快的成像型任意反射面速度干涉仪
CN111879421A (zh) * 2020-07-27 2020-11-03 中国科学院上海光学精密机械研究所 飞秒激光脉冲完全时间空间耦合特性单发测量系统及测量方法
CN114739521A (zh) * 2022-03-24 2022-07-12 中国科学院上海光学精密机械研究所 基于空频复用的单次曝光超短脉冲时空测量装置和方法
CN114894325A (zh) * 2022-05-06 2022-08-12 华中科技大学 一种阿秒脉冲时空结构的测量方法和测量装置
CN114894718A (zh) * 2022-04-14 2022-08-12 北京理工大学 基于各向异性超表面的相位复原方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050259249A1 (en) * 2004-01-27 2005-11-24 Dombeck Daniel A Nonlinear optical detection of fast cellular electrical activity
EP2870442B1 (en) * 2012-08-01 2017-03-22 Institut National De La Recherche Scientifique Spectral-domain interferometric method and system for characterizing terahertz radiation

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6639683B1 (en) * 2000-10-17 2003-10-28 Remy Tumbar Interferometric sensor and method to detect optical fields
CN106989834A (zh) * 2017-03-28 2017-07-28 中国工程物理研究院激光聚变研究中心 一种能同时诊断超短脉冲激光的啁啾特性与时空分布特性的方法
CN109100029A (zh) * 2018-09-05 2018-12-28 中国科学院上海光学精密机械研究所 飞秒激光脉冲时空参数单发测量装置
CN109411312A (zh) * 2018-10-31 2019-03-01 中国工程物理研究院流体物理研究所 基于飞秒激光调制的超快电子枪及其验证方法
CN111536868A (zh) * 2020-04-07 2020-08-14 华东师范大学 一种压缩超快的成像型任意反射面速度干涉仪
CN111879421A (zh) * 2020-07-27 2020-11-03 中国科学院上海光学精密机械研究所 飞秒激光脉冲完全时间空间耦合特性单发测量系统及测量方法
CN114739521A (zh) * 2022-03-24 2022-07-12 中国科学院上海光学精密机械研究所 基于空频复用的单次曝光超短脉冲时空测量装置和方法
CN114894718A (zh) * 2022-04-14 2022-08-12 北京理工大学 基于各向异性超表面的相位复原方法
CN114894325A (zh) * 2022-05-06 2022-08-12 华中科技大学 一种阿秒脉冲时空结构的测量方法和测量装置

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Single-frame measurement of the complete spatiotemporal intensity and phase of ultrashort laser pulses using wavelength-multiplexed digital holography;Pablo Gabolde等;《J.Opt.Soc.Am.B》;20080407;第25卷(第6期);A25-A33 *
Single-shot measurement of the complete spatiotemporal field of an ultrashort laser pulse;Pablo Gabolde等;《Proceedings of SPIE》;20060908;第6343卷;63433Z *
Single-shot spatiotemporal characterization of a multi-PW laser using a multispectral wavefront sensing method;Yeong Gyu Kim等;《Optics Express》;20210621;第29卷(第13期);19506-19514 *
基于空谱干涉和频域分割的超快激光时空耦合特性的单次测量方法;李伟;《物理学报》;20220215;第71卷(第3期);034203 *
基于空谱干涉扫描法测量超宽带激光时空耦合特性;李伟 等;《物理学报》;20211214;第70卷(第23期);234201 *
时空分辨的飞秒激光光场测量研究进展;李政言 等;《中国激光》;20220630;第49卷(第12期);1201003 *
超短脉冲时空畸变测量及补偿技术研究;周松 等;《物理学报》;20171231;第66卷(第4期);044204 *

Also Published As

Publication number Publication date
CN115235618A (zh) 2022-10-25

Similar Documents

Publication Publication Date Title
CN115235618B (zh) 一种超快激光场三维时空分布的单帧测量系统及方法
Millerd et al. Pixelated phase-mask dynamic interferometers
US8786864B2 (en) Circular common-path point diffraction interference wavefront sensor
US20030025911A1 (en) Optical pulse measurement
CN104655290A (zh) 斐索型双波长激光调谐移相干涉测试装置及其测试方法
CN113155058A (zh) 一种基于空域载频相移的多波长动态相位测量装置及方法
CN109000781B (zh) 一种结构微振动线域测量装置及方法
Li et al. Three-dimensional laser pulse intensity diagnostic for photoinjectors
CN113418469B (zh) 光谱共焦扫描共光路数字全息测量系统及测量方法
CN110160624B (zh) 一种用于三维振动测量的光纤点衍射装置及测量方法
CN108279068B (zh) 基于四波横向剪切干涉的激光光束质量动态测量装置
CN105277136B (zh) 基于双波长数字全息技术的透射式显微成像装置及其方法
CN112556842B (zh) 一种双波段高光谱分辨率闪电高速成像仪
CN102967379A (zh) 一种用于太阳自适应光学系统的波前传感器
CN116907804B (zh) 一种单脉冲超快激光时空耦合畸变测量仪
CN112595425B (zh) 一种超短激光脉冲测量方法及测量系统
CN109343321A (zh) X射线单次曝光相移径向剪切数字全息成像方法
CN109324023B (zh) 一种紧凑型差分干涉成像光谱仪及其成像方法
CN116379961B (zh) 一种相位测量系统及方法
CN104280139A (zh) 一种动态相位测量装置及方法
US20140098367A1 (en) Device And Method For Characterizing A Light Beam
CN104655029B (zh) 一种位相增强型薄膜厚度测量方法和系统
CN214095898U (zh) 一种基于偏振同步相移的轻离轴翻转干涉数字全息检测装置
CN113218517B (zh) 一种计算成像同轴全息光纤模式分解方法及装置
CN101782433A (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