CN116429252A - A Comprehensive Parameter Evaluation Method for Laser Beam Quality - Google Patents
A Comprehensive Parameter Evaluation Method for Laser Beam Quality Download PDFInfo
- Publication number
- CN116429252A CN116429252A CN202310411616.7A CN202310411616A CN116429252A CN 116429252 A CN116429252 A CN 116429252A CN 202310411616 A CN202310411616 A CN 202310411616A CN 116429252 A CN116429252 A CN 116429252A
- Authority
- CN
- China
- Prior art keywords
- wavefront
- laser
- light
- measured
- laser beam
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 8
- 238000005259 measurement Methods 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000003287 optical effect Effects 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 210000001747 pupil Anatomy 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
-
- 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
- G01J2009/002—Wavefront phase distribution
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明涉及一种激光光束质量的综合参数评价方法,属于激光光束特性测量领域。该方法通过多次测量平均来消除气流扰动的影响;考虑多模激光光源的尺寸,采用点扩散函数与光源进行卷积,来计算多模激光的远场光斑图像,解决多模激光测量错误的问题;在缩束或扩束时,采用扣除背景Wbg后的波前来计算远场光斑图像,减少远场光斑的误差;利用夏克哈特曼波前探测器测量波前、PV和RMS等参数,克服了现有的夏克哈特曼波前探测器测量激光光束质量因子M2中测量误差大,只能定性测量的问题,达到了准确测量单模或多模激光的光束质量因子、激光波前、PV和RMS等参数,实现激光光束质量综合参数的测量。
The invention relates to a comprehensive parameter evaluation method of laser beam quality, which belongs to the field of laser beam characteristic measurement. This method eliminates the influence of airflow disturbance by averaging multiple measurements; considering the size of the multi-mode laser light source, the point spread function is used to convolve with the light source to calculate the far-field spot image of the multi-mode laser, which solves the problem of multi-mode laser measurement errors Problem: When reducing or expanding the beam, use the wavefront after deducting the background W bg to calculate the far-field spot image to reduce the error of the far-field spot; use the Shack-Hartmann wavefront detector to measure the wavefront, PV and RMS and other parameters, it overcomes the problem that the existing Shack-Hartmann wavefront detector measures the laser beam quality factor M2 , which has a large measurement error and can only be measured qualitatively, and achieves accurate measurement of the beam quality factor of single-mode or multi-mode lasers , laser wavefront, PV and RMS and other parameters to realize the measurement of comprehensive parameters of laser beam quality.
Description
技术领域technical field
本发明涉及一种激光光束质量的综合参数评价方法,属于激光光束特性测量领域。The invention relates to a comprehensive parameter evaluation method of laser beam quality, which belongs to the field of laser beam characteristic measurement.
背景技术Background technique
1917年爱因斯坦提出“受激发射”理论,为激光的产生提供了重要的理论基础。1960年美国物理学家梅曼发明红宝石激光器,发出第一束可实际应用的激光。随着科学技术的不断进步,激光在军事、工业、医疗等领域中的应用越来越普遍。随着激光技术和应用的快速发展,检测激光光束质量也越来越重要,国际标准化组织将用于评价光束质量的参数——激光光束质量因子M2定义为实际光束的束腰半径(光束传播方向上光束半径最小的地方称为束腰)与其远场发散角的乘积与理想高斯基模光束的相应乘积之比,由此,为了评价光束质量,需要获取实际光束的束腰半径以及其远场发散角,这两个参数都属于聚焦光斑尺寸参数。In 1917, Einstein proposed the theory of "stimulated emission", which provided an important theoretical basis for the generation of laser light. In 1960, the American physicist Maiman invented the ruby laser, which issued the first practically applicable laser. With the continuous advancement of science and technology, the application of lasers in military, industrial, medical and other fields is becoming more and more common. With the rapid development of laser technology and applications, it is more and more important to detect the quality of laser beams. The International Organization for Standardization defines the parameters used to evaluate the quality of beams—the laser beam quality factor M 2 as the beam waist radius of the actual beam (beam propagation The place where the beam radius is the smallest in the direction is called the beam waist) and the product of its far-field divergence angle and the corresponding product of the ideal Gaussian mode beam. Therefore, in order to evaluate the beam quality, it is necessary to obtain the beam waist radius of the actual beam and its far-field Field divergence angle, both of these parameters belong to the focusing spot size parameter.
现有测量聚焦光斑尺寸的方法有阵列探测法和机械扫描法。其中机械扫描法又有狭缝扫描法、移动刀口法以及可变光阑法,在国外,机械扫描技术已经比较成熟,许多国外的公司研发专用的光束质量测量仪。比如赵俊成等人采用一种滚轮狭缝式的测量方法[赵俊成,阴万宏,刘建平,等.正交狭缝扫描式激光光束质量测量方法研究[J].应用光学,2020,41(04):704-710]就是一种机械扫描技术,直接利用这种滚轮扫描狭缝扫描被测激光光束进而计算光束质量。而阵列探测技术(参见[Sheldakova Y V,Cherezova T Y,Kudryashov A V.M2 sensor for the adaptive optical system[C].SeventhInternational Conference on Laser and Laser-InformationTechnologies.International Society for Optics and Photonics,2002,4644:392-399])是基于CCD相机等光电探测器而衍生的一种比较具有发展前景的光束检测技术,可以通过CCD直接测得激光光强分布,并且可以用相机像素值获得每一点的光强度值。根据CCD相机上光强度值分布范围可以计算出激光的光斑尺寸,即从不同位置的光强度衰减值情况就可以得到光斑尺寸。根据激光传输性质,利用不同位置的光斑尺寸测量数据就可以拟合激光传输特性曲线,进而得到激光远场发散角和光束质量因子等光束质量参数,由于使用CCD不需要改变激光光束内部结构,因此不会对激光光束质量产生影响,所以阵列探测法也是目前测量光束质量的最常用一种方法。但是无论是机械扫描法还是阵列探测法都是采用多点法测量光束质量,这种方法不能实时测量光束质量,而且往往会引入一定的测量误差,还会因为误差累计导致最终计算出的M2存在更大的误差。The existing methods for measuring the focus spot size include array detection method and mechanical scanning method. Among them, the mechanical scanning method includes the slit scanning method, the moving knife-edge method and the iris method. In foreign countries, the mechanical scanning technology has been relatively mature, and many foreign companies have developed special beam quality measuring instruments. For example, Zhao Juncheng and others adopted a roller slit measurement method [Zhao Juncheng, Yin Wanhong, Liu Jianping, et al. Research on the quality measurement method of laser beams with orthogonal slit scanning [J]. Applied Optics, 2020, 41( 04):704-710] is a mechanical scanning technology, which directly uses the roller to scan the slit to scan the measured laser beam and then calculate the beam quality. And array detection technology (see [Sheldakova YV, Cherezova TY, Kudryashov A V. M2 sensor for the adaptive optical system [C]. Seventh International Conference on Laser and Laser-Information Technologies. International Society for Optics and Photonics, 2002, 4644: 392- 399]) is a relatively promising beam detection technology derived from photodetectors such as CCD cameras. The laser light intensity distribution can be directly measured through the CCD, and the light intensity value of each point can be obtained by using the camera pixel value. The spot size of the laser can be calculated according to the distribution range of the light intensity value on the CCD camera, that is, the spot size can be obtained from the light intensity attenuation values at different positions. According to the properties of laser transmission, the laser transmission characteristic curve can be fitted by using the measurement data of spot size at different positions, and then the beam quality parameters such as laser far-field divergence angle and beam quality factor can be obtained. Since the use of CCD does not need to change the internal structure of the laser beam, so It will not affect the quality of the laser beam, so the array detection method is also the most commonly used method to measure the quality of the beam. However, whether it is the mechanical scanning method or the array detection method, the multi-point method is used to measure the beam quality. This method cannot measure the beam quality in real time, and often introduces certain measurement errors, and the final calculated M 2 will also be caused by the accumulation of errors. There are larger errors.
许多国内外学者开始研究如何实时测量光束质量,他们发现可以利用波前相位重构复振幅,再利用光束衍射传输公式计算远场发散角和束腰宽度进而测量出光束质量,如基于斜率测量恢复相位的夏克哈特曼波前测量技术[Sheldakova J V,Kudryashov A V,Zavalova V Y.Beam quality measurements with Shack-Hartmann wavefront sensorand M2-sensor:comparison of two methods[C].Laser Resonators and Beam ControlIX,2007,6452:645207]。但它测量的数值不准确,只能定性地进行测量,且只能针对单模激光进行测量,对于多模激光测量结果是错误的,这些限制了它的实际应用。主要是因为在实际光路测试过程中,光路中透镜、气流等影响会引入畸变,光路的畸变即使可以通过非常精密地进行光路装调来减少,但无法完全消除;室内气流扰动引起的畸变不能通过光路的精密装调来消除,当这些畸变参与到远场光斑的计算时,它会影响光斑的大小,使计算M2结果出现误差,因此,常规的夏克哈特曼波前测量技术只能定性地测量激光光束质量因子M2,无法广泛应用。Many scholars at home and abroad have begun to study how to measure the beam quality in real time. They found that the complex amplitude can be reconstructed using the wavefront phase, and then the beam diffraction and transmission formula can be used to calculate the far-field divergence angle and beam waist width to measure the beam quality. For example, based on the slope measurement recovery Shack-Hartmann wavefront measurement technique of phase[Sheldakova JV, Kudryashov AV, Zavalova V Y. Beam quality measurements with Shack-Hartmann wavefront sensor and M2-sensor: comparison of two methods[C].Laser Resonators and Beam ControlIX,2007 ,6452:645207]. But the value it measures is inaccurate, it can only be measured qualitatively, and it can only be measured for single-mode lasers, and the measurement results for multi-mode lasers are wrong, which limits its practical application. The main reason is that in the actual optical path test process, the lens, airflow and other influences in the optical path will introduce distortion. Even if the optical path distortion can be reduced by very precise optical path adjustment, it cannot be completely eliminated; the distortion caused by indoor airflow disturbance cannot be passed. When these distortions are involved in the calculation of the far-field spot, it will affect the size of the spot and cause errors in the calculation of M2 . Therefore, the conventional Shack-Hartmann wavefront measurement technology can only The qualitative measurement of the laser beam quality factor M 2 cannot be widely used.
针对前面的问题,采用夏克哈特曼波前传感器代替传统接收器CCD来进行光束质量评价,同时采用合理的方法来克服不利因素的影响,减少测量误差,夏克哈特曼波前传感器所构成设备操作方便,并且可以有效的缩短光束评价时间。采用新的算法测量激光光束时,我们可以同时获得激光光束的波前分布、波前的PV和RMS值、光强分布、光强峰值、激光光束质量因子M2等关键参数,这种综合参数对激光研究具有重要的意义和广泛的应用价值。In view of the previous problems, the Shack-Hartmann wavefront sensor is used instead of the traditional receiver CCD to evaluate the beam quality, and a reasonable method is used to overcome the influence of unfavorable factors and reduce measurement errors. The Shack-Hartmann wavefront sensor The device is easy to operate and can effectively shorten the beam evaluation time. When using the new algorithm to measure the laser beam, we can simultaneously obtain key parameters such as the wavefront distribution of the laser beam, the PV and RMS values of the wavefront, the light intensity distribution, the peak light intensity, and the laser beam quality factor M2 . It is of great significance and wide application value to laser research.
发明内容Contents of the invention
为了解决目前存在的夏克哈特曼波前测量技术只能定性地测量激光光束质量因子M2问题,本发明提供了一种激光光束质量的综合参数评价方法,所述方法基于夏克哈提曼波前探测器SHWFS实现,所述夏克哈提曼波前探测器由微透镜阵列和光电探测器构成,所述微透镜阵列中的所有微透镜具有相同的直径和焦距;所述方法包括:In order to solve the existing Shaker-Hartman wavefront measurement technology can only qualitatively measure the laser beam quality factor M 2 problem, the invention provides a kind of comprehensive parameter evaluation method of laser beam quality, said method is based on Shaker-Hartman The wavefront detector SHWFS is realized, and the Shaker Hartimann wavefront detector is composed of a microlens array and a photodetector, and all microlenses in the microlens array have the same diameter and focal length; the method includes:
步骤1,以平行光作为光源,利用所述光电探测器测量得到光路的背景畸变Wbg;
步骤2,以待测激光作为光源,利用所述光电探测器多次采集激光光束的光点阵图案,获得畸变波前,并对多次测量得到的畸变波前取平均值,利用平均的畸变波前,计算峰谷值PV和均方根RMS;
步骤3,当待测激光光源为单模激光时,利用扣除背景Wbg后的畸变波前计算远场光斑图像,获得实际远场光斑的图像;当待测激光光源为多模激光时,考虑光路的缩放比例和光源的尺寸,利用扣除背景Wbg后的畸变波前计算的点扩散函数与光源缩放后的图像进行卷积,获得实际远场光斑的图像;
步骤4,根据远场光斑的图像,分别计算待测激光的远场光斑尺寸及理想远场光斑的尺寸,并利用它们的比值计算激光光束质量因子M2参数。Step 4: Calculate the far-field spot size of the laser to be tested and the ideal far-field spot size according to the far-field spot image, and use their ratio to calculate the laser beam quality factor M 2 parameter.
可选的,所述步骤1包括:Optionally, the
以平行光作为光源,利用所述光电探测器获取相应的光点阵图案,所述光点阵图案中每个子孔径内光点作为参考光点的位置,所有子孔径内光点相对于孔径中心位置的偏移称为光路的背景畸变Wbg。Using parallel light as a light source, using the photodetector to obtain a corresponding light lattice pattern, the light point in each sub-aperture in the light lattice pattern is used as the position of the reference light point, and the light points in all sub-apertures are relative to the center of the aperture. The shift in position is called the background distortion W bg of the optical path.
可选的,所述步骤2中利用所述光电探测器多次采集激光光束的光点阵图案,获得畸变波前,包括:Optionally, in the
以待测激光作为光源,利用所述光电探测器多次采集激光光束的光点阵图案;Using the laser to be measured as a light source, using the photodetector to collect the light lattice pattern of the laser beam multiple times;
测量光点阵图案中每个子孔径内光点相对于参考光点位置的偏移量;Measuring the offset of the light spot in each sub-aperture in the light dot matrix pattern relative to the position of the reference light spot;
根据每个子孔径内光点相对于参考光点位置的偏移量计算子孔径中的波前斜率,并基于畸变波前的局部斜率的测量来重构畸变波前,得到整个畸变波前W:Calculate the wavefront slope in each subaperture based on the offset of the spot relative to the reference spot position within each subaperture and based on the local slope of the distorted wavefront The measurement of the distortion wavefront is reconstructed, and the entire distortion wavefront W is obtained:
其中,(Gix,Giy)是SHWFS测得的第i个子透镜测得的沿x方向、沿y方向上的波前斜率值;Zi表示第i项Zernike模式,i取值范围1到n,n为微透镜阵列中包含的透镜个数;[a1... an]T表示Zernike多项式的系数;Among them, (G ix ,G iy ) is the wavefront slope value along the x direction and along the y direction measured by the i-th sub-lens measured by SHWFS; Z i represents the i-th Zernike mode, and the value of i ranges from 1 to n, n is the number of lenses included in the microlens array; [a 1 ... a n ] T represents the coefficient of the Zernike polynomial;
将式(1)表示成矩阵形式:Express formula (1) in matrix form:
G=F·A (2)G=F·A (2)
其中,G是由SHWFS测得的波前斜率向量,F是波前重构矩阵,其维度为2m×n,A是Zernike多项式的系数向量;Among them, G is the wavefront slope vector measured by SHWFS, F is the wavefront reconstruction matrix with a dimension of 2m×n, and A is the coefficient vector of the Zernike polynomial;
根据测量得到的斜率G,则待测波前的Zernike系数向量A的最小二乘解为:According to the measured slope G, the least squares solution of the Zernike coefficient vector A of the wavefront to be measured is:
A=(FTF)-1G (3)A=(F T F) -1 G (3)
求出系数A后,代入到公式(4)即可求得整个畸变波前W的相位信息;After the coefficient A is obtained, it can be substituted into formula (4) to obtain the phase information of the entire distorted wavefront W;
由此,得到待测波前W的分布情况,W是由N×N格点组成的矩阵,其中,有效数据为内截圆的区域中的数据,共M个。Thus, the distribution of the wavefront W to be measured is obtained. W is a matrix composed of N×N grid points, wherein the effective data are the data in the area of the inner truncated circle, and there are M pieces in total.
可选的,假设步骤2中共采集M次激光光束的光点阵图案,每次采集间隔固定时间;所述步骤2中峰谷值PV和均方根RMS计算公式如下:Optionally, it is assumed that the light lattice pattern of the laser beam is collected M times in
PV=Wmax-Wmin (5)PV=W max -W min (5)
其中,Wmax表示单次测量得到的波前W中M个有效数据中的最大值,Wmin表示单次测量得到的波前W中M个有效数据中的最小值;Wj表示根据第j次测量得到的波前。Among them, W max represents the maximum value among the M effective data in the wavefront W obtained by a single measurement, W min represents the minimum value among the M valid data in the wavefront W obtained by a single measurement; W j represents the The measured wavefront.
可选的,所述步骤3中,当待测激光光源为单模激光时,所述实际远场光斑图像的直径直接根据畸变波前W计算光强分布:Optionally, in
其中,I(x,y)表示激光光束的幅度;W(x,y)表示激光光束的相位;x表示x方向上的坐标,y表示y方向上的坐标;λ是激光波长,exp表示自然常数e为底的指数函数;i表示复数的基本单位。in, I(x,y) represents the amplitude of the laser beam; W(x,y) represents the phase of the laser beam; x represents the coordinate in the x direction, y represents the coordinate in the y direction; λ is the laser wavelength, exp represents the natural constant e Base exponential function; i represents the basic unit of complex numbers.
可选的,所述步骤3中,当待测激光光源为多模激光时,计算其点扩散函数PSF:Optionally, in
则实际光斑为光源的理想成像O与点扩散函数PSF的卷积:Then the actual spot is the convolution of the ideal image O of the light source and the point spread function PSF:
I=O*PSF (9)。I=O*PSF (9).
可选的,所述步骤4激光光束质量因子M2参数计算公式如下:Optionally, the formula for calculating the parameters of the laser beam quality factor M in step 4 is as follows:
其中,d0是近场中的光束直径,发散角df是焦点为f的透镜焦平面上的光束直径,因此,where d0 is the beam diameter in the near field and the divergence angle df is the beam diameter at the focal plane of the lens with focal point f, so,
其中,d(Ifar field)表示待测激光的远场光斑尺寸;dgauss表示理想远场光斑的尺寸。Among them, d(I far field ) represents the far-field spot size of the laser to be measured; d gauss represents the ideal far-field spot size.
可选的,若所述待测激光的直径与所述夏克哈特曼波前探测器的口径不匹配,则使用透镜组进行缩束或扩束。Optionally, if the diameter of the laser to be measured does not match the aperture of the Shack-Hartmann wavefront detector, a lens group is used to reduce or expand the beam.
可选的,所述微透镜阵列至少包含20*20个微透镜。Optionally, the microlens array includes at least 20*20 microlenses.
本发明有益效果是:The beneficial effects of the present invention are:
通过多次测量平均来消除气流扰动的影响;考虑多模激光光源的尺寸,采用点扩散函数与光源进行卷积,来计算多模激光的远场光斑图像,解决多模激光测量错误的问题;在缩束或扩束时,采用扣除背景Wbg后的波前来计算远场光斑图像,减少远场光斑的误差;利用夏克哈特曼波前探测器测量波前、PV和RMS等参数,克服了现有的夏克哈特曼波前探测器测量激光光束质量因子M2中测量误差大,只能定性测量的问题,达到了准确测量单模或多模激光的光束质量因子、激光波前、PV和RMS等参数,实现激光光束质量综合参数的测量。Eliminate the influence of airflow disturbance by averaging multiple measurements; consider the size of the multi-mode laser light source, use the point spread function to convolve with the light source to calculate the far-field spot image of the multi-mode laser, and solve the problem of multi-mode laser measurement errors; When shrinking or expanding the beam, use the wavefront after deducting the background W bg to calculate the far-field spot image, reduce the error of the far-field spot; use the Shack-Hartmann wavefront detector to measure parameters such as wavefront, PV and RMS , overcomes the problem that the existing Shaker-Hartmann wavefront detector measures the laser beam quality factor M 2 with large measurement errors and can only be measured qualitatively, and achieves accurate measurement of the beam quality factor of single-mode or multi-mode lasers, laser Wavefront, PV, RMS and other parameters to realize the measurement of comprehensive parameters of laser beam quality.
附图说明Description of drawings
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings that need to be used in the description of the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present invention. For those skilled in the art, other drawings can also be obtained based on these drawings without creative effort.
图1为需要缩束或扩束时测量光路原理示意图;1为激光器,2和3为透镜,4为夏克哈特曼波前探测器。Figure 1 is a schematic diagram of the principle of the measurement optical path when beam shrinkage or beam expansion is required; 1 is a laser, 2 and 3 are lenses, and 4 is a Shack-Hartmann wavefront detector.
图2为直接测量光路原理示意图。1为激光器,4为夏克哈特曼波前探测器。Figure 2 is a schematic diagram of the principle of the direct measurement optical path. 1 is a laser, and 4 is a Shack-Hartmann wavefront detector.
图3为实验光路图;有633nm的激光器,用于衰减光强的衰减片,反射镜用于光路的折轴,第一透镜和第二透镜组成一个扩束器扩大激光光束的直径;最后用夏克哈特曼波前探测器对激光光束进行测量。Fig. 3 is the experimental light path diagram; The laser device of 633nm is arranged, the attenuation plate that is used for attenuating light intensity, and reflective mirror is used for the folding axis of light path, and the first lens and the second lens form a beam expander to expand the diameter of laser beam; Finally use The Shack-Hartmann wavefront detector measures the laser beam.
图4为夏克哈特曼波前探测器参考光斑和畸变光斑的仿真图,其中(a)为平面波形成的光点阵图;(b)为畸变波前像差的光点阵图。Fig. 4 is a simulation diagram of the reference spot and the distorted spot of the Shack-Hartmann wavefront detector, wherein (a) is a dot matrix diagram formed by a plane wave; (b) is a dot matrix diagram of a distorted wavefront aberration.
图5为波前斜率计算原理示意图。Fig. 5 is a schematic diagram of the calculation principle of the wavefront slope.
图6中(a)为畸变波前对应的点扩散函数仿真图,(b)为畸变波前仿真图,(c)为光点阵图。In Fig. 6, (a) is the simulation diagram of the point spread function corresponding to the distorted wavefront, (b) is the simulation diagram of the distorted wavefront, and (c) is the light dot matrix diagram.
图7中(a)为平面波的点扩散函数仿真图,(b)瞳孔函数仿真图。In Fig. 7, (a) is a simulation diagram of a point spread function of a plane wave, and (b) is a simulation diagram of a pupil function.
具体实施方式Detailed ways
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。In order to make the purpose, technical solution and advantages of the present invention clearer, the following will further describe in detail the embodiments of the present invention in conjunction with the accompanying drawings.
实施例一:Embodiment one:
本实施例提供一种激光光束质量的综合参数评价方法,主要步骤如下:This embodiment provides a comprehensive parameter evaluation method for laser beam quality, the main steps are as follows:
1)当激光光束的直径与夏克哈特曼波前探测器的口径不匹配,需要用两个透镜进行缩束或扩束,如图1所示,1为激光器,2和3为透镜,4为夏克哈特曼波前探测器。透镜2和透镜3的焦点重合,激光器1位于透镜2的物面,夏克哈特曼波前探测器4位于透镜3的像面。当激光光束的直径与夏克哈特曼波前探测器的口径匹配时,可直接测量,如图2所示,1为激光器,4为夏克哈特曼波前探测器。1) When the diameter of the laser beam does not match the aperture of the Shack-Hartmann wavefront detector, two lenses are needed to shrink or expand the beam, as shown in Figure 1, 1 is the laser, 2 and 3 are lenses, 4 is the Shack-Hartmann wavefront detector. The focal points of the
2)对于图1的光路,可以利用zygo干涉仪的平行光输出作为光源,替代1光源,此时,夏克哈特曼波前探测器测量的畸变,即为对光路的畸变进行标定,作为光路的背景畸变Wbg。2) For the optical path in Figure 1, the parallel light output of the zygo interferometer can be used as the light source instead of 1 light source. At this time, the distortion measured by the Shack-Hartmann wavefront detector is to calibrate the distortion of the optical path, as Background distortion W bg of the light path.
3)利用夏克哈特曼波前探测器采集光点阵图像,获得畸变波前,间隔1秒重复测量波前,测量次数至少10次,对测量的波前取平均,以消除室内气流扰动的影响;利用平均的畸变波前,计算峰谷值PV和均方根RMS;3) Use the Shack-Hartmann wavefront detector to collect light lattice images to obtain the distorted wavefront, repeat the wavefront measurement at intervals of 1 second, and measure the wavefront at least 10 times, and average the measured wavefront to eliminate indoor airflow disturbance The influence of; use the average distortion wavefront to calculate the peak-to-valley PV and root mean square RMS;
4)当待测激光光源为单模激光时,利用扣除背景Wbg后的畸变波前计算远场光斑图像,获得实际远场光斑的图像;当待测激光光源为多模激光时,需要考虑光路的缩放比例和光源的尺寸,将利用扣除背景Wbg后的畸变波前计算的点扩散函数与光源缩放后的图像进行卷积,获得实际远场光斑的图像;4) When the laser source to be tested is a single-mode laser, the far-field spot image is calculated by using the distorted wavefront after deducting the background Wbg to obtain the actual far-field spot image; when the laser source to be tested is a multi-mode laser, it is necessary to consider For the scaling ratio of the optical path and the size of the light source, the point spread function calculated by using the distorted wavefront after deducting the background W bg is convolved with the image scaled by the light source to obtain the actual far-field spot image;
5)根据远场光斑的图像,分别计算待测激光的远场光斑尺寸及理想远场光斑的尺寸,并利用它们的比值计算激光光束质量因子M2参数;5) Calculate the far-field spot size of the laser to be measured and the ideal far-field spot size according to the far-field spot image, and use their ratio to calculate the laser beam quality factor M parameter;
实施例二:Embodiment two:
本实施例提供一种激光光束质量的综合参数评价方法,所述方法用于同时测量激光光束波前分布、波前的PV和RMS值、光强分布、光强峰值、M2因子等参数,本实施例以激光光束的直径与夏克哈特曼波前探测器的口径不匹配为例进行介绍,搭建图3所示的光路系统,该光路系统包括:激光器、衰减片、反射镜、第一透镜、第二透镜和夏克哈提曼波前探测器SHWFS;所述夏克哈提曼波前探测器由微透镜阵列和位于透镜焦点位置的光电探测器CCD构成,微透镜阵列至少包含20*20个微透镜,且所有微透镜具有相同的直径和焦距。激光器发出的激光束通过衰减片后经反射镜改变光束方向后,依次透过第一透镜和第二透镜后到达夏克哈提曼波前探测器,入射到SHWFS的波前被微透镜分割成近似的平面波然后成像于后焦面上的CCD,并形成相应的光点阵图案。The present embodiment provides a kind of comprehensive parameter evaluation method of laser beam quality, and described method is used for simultaneously measuring parameters such as laser beam wavefront distribution, wavefront PV and RMS value, light intensity distribution, light intensity peak value, M factor , In this embodiment, the diameter of the laser beam does not match the aperture of the Shack-Hartmann wavefront detector as an example. The optical path system shown in Figure 3 is built. The optical path system includes: lasers, attenuators, mirrors, and A lens, a second lens, and a Shaker Hartimann wavefront detector SHWFS; the Shaker Hartimann wavefront detector is composed of a microlens array and a photodetector CCD positioned at the focal point of the lens, and the microlens array includes at least 20* 20 microlenses, and all microlenses have the same diameter and focal length. The laser beam emitted by the laser passes through the attenuation sheet, changes the direction of the beam through the mirror, passes through the first lens and the second lens in turn, and then reaches the Shaker-Hartimann wavefront detector. The wavefront incident on the SHWFS is divided into approximately The plane wave is then imaged on the CCD on the back focal plane, and forms a corresponding light lattice pattern.
基于所述光点阵图案测量波前,包括:Measuring a wavefront based on the light lattice pattern, comprising:
由于微透镜阵列将整个波前分割成多个子孔径,故首先测量光点阵图案中每个子孔径内光点相对于参考光点位置的偏移量,请参考图4;参考光点位置可以利用zygo干涉仪的平行光输出作为光源替代激光器得到。Since the microlens array divides the entire wavefront into multiple sub-apertures, first measure the offset of the light spot in each sub-aperture in the light lattice pattern relative to the reference light spot position, please refer to Figure 4; the reference light spot position can be used The parallel light output of the zygo interferometer is obtained as a light source instead of a laser.
根据每个子孔径内光点相对于参考光点位置的偏移量计算子孔径中的波前斜率(可参考图5),并基于畸变波前的局部斜率的测量来重构畸变波前;Calculate the wavefront slope in the subaperture according to the offset of the light spot in each subaperture relative to the reference spot position (see Figure 5), and based on the local slope of the distorted wavefront measurements to reconstruct the distorted wavefront;
式(1)中(Gix,Giy)是SHWFS测得的第i个子透镜测得的沿x方向和y方向上的波前斜率值;Zi表示第i项Zernike模式,i取1到n的值。In formula (1), (G ix ,G iy ) is the wavefront slope value along the x-direction and y-direction measured by the i-th sub-lens measured by SHWFS; Z i represents the i-th Zernike mode, i ranges from 1 to The value of n.
将式(1)表示成矩阵形式:Express formula (1) in matrix form:
G=F·A (2)G=F·A (2)
其中,G是由SHWFS测得的波前斜率向量,F是波前重构矩阵,其维度为2m×n,A是Zernike多项式的系数向量。因此,根据测量得到的斜率G,则待测波前的Zernike系数向量A的最小二乘解为:Among them, G is the wavefront slope vector measured by SHWFS, F is the wavefront reconstruction matrix with a dimension of 2m×n, and A is the coefficient vector of the Zernike polynomial. Therefore, according to the measured slope G, the least squares solution of the Zernike coefficient vector A of the wavefront to be measured is:
A=(FTF)-1G (3)A=(F T F) -1 G (3)
求出系数A后,代入到公式(4)即可求得整个畸变波前W的相位信息。After the coefficient A is obtained, it can be substituted into formula (4) to obtain the phase information of the entire distorted wavefront W.
由此,得到待测波前的分布情况。Thus, the distribution of the wavefront to be measured is obtained.
在此基础上,根据波前W的最大值减去最小值得到波前的PV值,根据波前W计算出它的均方根RMS:On this basis, the PV value of the wave front is obtained by subtracting the minimum value from the maximum value of the wave front W, and its root mean square RMS is calculated according to the wave front W:
PV=Wmax-Wmin (5)PV=W max -W min (5)
SHWFS的微透镜数目足够多时,例如超过20x20个微透镜,SHWFS的CCD上采集到的光强也反映了垂直光束传播方向上的光强分布。When the number of microlenses in SHWFS is large enough, for example, more than 20x20 microlenses, the light intensity collected by the CCD of SHWFS also reflects the light intensity distribution in the direction of vertical beam propagation.
获得激光光束的相位W(x,y)和幅度I(x,y)后,可以计算出M2光束参数。After obtaining the phase W(x,y) and amplitude I(x,y) of the laser beam, the M2 beam parameters can be calculated.
在已有的文献中计算远场光斑直径是直接根据畸变波前W计算光强分布,但这实际是点扩散函数,因此,对于点光源来说,这种方法可行,但对于有一定大小的光源来说就会出现较大的误差。因此,对于有一定大小的目标,本申请用傅里叶变换算法计算出点扩散函数后,再用点扩散函数对目标进行卷积得到激光光斑的实际图像。In the existing literature, the calculation of the far-field spot diameter is to directly calculate the light intensity distribution based on the distorted wavefront W, but this is actually a point spread function. Therefore, this method is feasible for a point light source, but for a certain size For the light source, there will be a large error. Therefore, for a target with a certain size, the application uses the Fourier transform algorithm to calculate the point spread function, and then uses the point spread function to convolve the target to obtain the actual image of the laser spot.
对于点光源:For point lights:
其中, in,
对于有一定尺寸的光源,首先计算其点扩散函数PSF:For a light source with a certain size, first calculate its point spread function PSF:
其中, in,
则实际光斑为光源的理想成像O与点扩散函数PSF的卷积:Then the actual spot is the convolution of the ideal image O of the light source and the point spread function PSF:
I=O*PSF (9)I=O*PSF (9)
由国际标准ISO/DIS 11146中给出的关系式计算得出激光光束质量因子M2:The laser beam quality factor M 2 is calculated from the relation given in the international standard ISO/DIS 11146:
其中,d0是近场中的光束直径,θ是发散角,λ是激光波长。where d0 is the beam diameter in the near field, θ is the divergence angle, and λ is the laser wavelength.
发散角θ由方程定义,其中df是焦点为f的透镜焦平面上的光束直径。公式(10)也可以近似表示为:The divergence angle θ is given by the equation Definition, where df is the beam diameter at the focal plane of the lens with focal point f. Formula (10) can also be approximated as:
至此,我们可以公式(1)~(4)计算波前W,根据(5)和(6)计算PV和RMS,根据(7)~(9)可计算点光源或有一定尺寸光源的远场光强分布。最后根据(10)或(11)计算出激光光束质量因子M2。So far, we can calculate the wavefront W according to formulas (1)~(4), calculate PV and RMS according to (5) and (6), and calculate the far field of a point light source or a light source with a certain size according to (7)~(9) light intensity distribution. Finally, the laser beam quality factor M 2 is calculated according to (10) or (11).
本实施例中激光器为单模激光,激光波长λ=633nm,第一透镜的焦点f=25nm,第二透镜的焦点f=100nm,夏克哈提曼波前探测器的微透镜阵列包含29*29个微透镜。In the present embodiment, the laser is a single-mode laser, the laser wavelength λ=633nm, the focal point f=25nm of the first lens, the focal point f=100nm of the second lens, and the microlens array of the Shaker-Hartimann wavefront detector includes 29*29 a microlens.
基于上述搭建的光路,所得到的畸变波前对应的点扩散函数如图6中的(a),图6中的(b)为畸变波前,图6中的(c)为光点阵;图7中(a)为平面波的点扩散函数,图7中(b)为瞳孔函数。计算结果M2为10.48,由此可见,本申请方案可以定量地测量激光光束质量因子M2。Based on the optical path built above, the point spread function corresponding to the obtained distorted wavefront is shown in (a) in Figure 6, (b) in Figure 6 is the distorted wavefront, and (c) in Figure 6 is the light lattice; (a) in Fig. 7 is the point spread function of the plane wave, and (b) in Fig. 7 is the pupil function. The calculation result M 2 is 10.48. It can be seen that the scheme of the present application can quantitatively measure the laser beam quality factor M 2 .
本发明实施例中的部分步骤,可以利用软件实现,相应的软件程序可以存储在可读取的存储介质中,如光盘或硬盘等。Part of the steps in the embodiments of the present invention can be realized by software, and the corresponding software program can be stored in a readable storage medium, such as an optical disk or a hard disk.
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included in the protection of the present invention. within range.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310411616.7A CN116429252A (en) | 2023-04-18 | 2023-04-18 | A Comprehensive Parameter Evaluation Method for Laser Beam Quality |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310411616.7A CN116429252A (en) | 2023-04-18 | 2023-04-18 | A Comprehensive Parameter Evaluation Method for Laser Beam Quality |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116429252A true CN116429252A (en) | 2023-07-14 |
Family
ID=87079432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310411616.7A Pending CN116429252A (en) | 2023-04-18 | 2023-04-18 | A Comprehensive Parameter Evaluation Method for Laser Beam Quality |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116429252A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118376157A (en) * | 2024-06-21 | 2024-07-23 | 钛玛科(北京)工业科技有限公司 | Method, device and equipment for detecting size of collimated beam-expanding laser |
CN118518211A (en) * | 2024-07-23 | 2024-08-20 | 中国科学院长春光学精密机械与物理研究所 | High-energy laser atmospheric transmission effect simulation and detection device and method |
-
2023
- 2023-04-18 CN CN202310411616.7A patent/CN116429252A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118376157A (en) * | 2024-06-21 | 2024-07-23 | 钛玛科(北京)工业科技有限公司 | Method, device and equipment for detecting size of collimated beam-expanding laser |
CN118518211A (en) * | 2024-07-23 | 2024-08-20 | 中国科学院长春光学精密机械与物理研究所 | High-energy laser atmospheric transmission effect simulation and detection device and method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110160751B (en) | A broadband wavefront error detection device and detection method based on phase recovery | |
CN104034416B (en) | High-dynamic-range laser far-field focal spot measuring device and method | |
CN116429252A (en) | A Comprehensive Parameter Evaluation Method for Laser Beam Quality | |
CN110375964B (en) | A wavefront error detection device and detection method based on extended Naipor-Zernike mode optimized phase recovery | |
CN102252832B (en) | Large aperture collimation system wavefront quality detection device and method | |
CN102564611B (en) | High-power laser wave front measuring instrument and wave front measuring method | |
CN111551129B (en) | Medium-low-order surface shape detection device and system of large-caliber plane mirror and storage medium | |
CN110726381A (en) | Optical free-form surface full-band aberration detection system and detection method | |
CN114216659A (en) | System and method for measuring parallelism of large-caliber long-focus optical axis | |
CN101936779B (en) | Double-optical-wedge spliced rectangular pyramid wavefront sensor | |
Chen et al. | Transition imaging phase measuring deflectometry for high-precision measurement of optical surfaces | |
CN113834421B (en) | Imaging lens group and interferometer using same | |
CN115127473A (en) | A large-diameter optical element surface shape detection device and detection image distortion correction method | |
CN111076904B (en) | A high-power thin-film laser dynamic wavefront aberration detection device and method | |
CN210294682U (en) | Dual-channel infrared scene simulator device | |
JP2003121300A (en) | Wavefront aberration measurement method and device of optical system | |
CN106770335B (en) | A Phase Defect Detection System and Method Based on Reflective Point Diffraction Interferometer | |
CN117571263A (en) | An optical performance measurement method for free-form afocal optical systems | |
CN110376732A (en) | A kind of binary channels Infrared scene simulation with image device device | |
CN100562723C (en) | Positive branch confocal unstable intracavity aberration detection system | |
CN107589542B (en) | A Method for Selecting the Center Wavelength in Reconstruction of Broadband Phase-difference Image | |
US11333487B2 (en) | Common path mode fiber tip diffraction interferometer for wavefront measurement | |
JP6823334B2 (en) | Outlet wave surface measurement method and outlet wave surface measurement system for high NA condensing elements | |
CN115585894A (en) | Wavefront detection method for realizing extremely small four-wavefront transverse shear rate based on 4f system | |
RU2680615C1 (en) | Method for determining deformations of wavefront of light beam caused by waviness of surfaces of optical system |
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 |