Method for measuring coherent distance, and apparatus for measuring coherent distance

Abstract

Measuring the coherence distance of a coherence light source operating in a very high frequency band. [Solution] The light to be measured is output from a coherence wavelength-swept light source, separated into signal light and reference light by a first optical coupler, and combined by a second optical coupler. These are then output as a first interference signal, a second interference signal, and a third interference signal, each with a phase difference of 120°. Each of these is detected by a photodetector. A correction matrix is ​​calculated in which factors related to the amplitude of the first and second optical couplers, the photoelectric efficiency of the photodetector, and the phase of the second optical coupler change over time. The phase difference between the signal light and the reference light is calculated using the first, second, and third interference signals and the correction matrix. The spectral width of the light is calculated using the phase difference, and the coherence distance of the light is calculated using the spectral width.

Description

[Technical Field]

[0001] This disclosure relates to a coherence distance measurement method and a coherence distance measurement apparatus. [Background technology]

[0002] To easily measure the oscillation linewidth of a coherent light source during wavelength sweeping, it is known that the light to be measured, which is coherent light whose oscillation wavelength changes periodically, is generated from the wavelength-swept light source, is split into a measurement unit and a reference unit with different optical path lengths, and then combined in an interference unit to obtain interference light, which is detected by a photodetector unit to obtain an interference signal, and the analysis unit performs a Fourier transform on the interference signal obtained by changing the optical path length difference between the measurement unit and the reference unit to obtain the measured value of the noise floor value of the point image distribution function for each optical path length difference, calculates the measured amplitude value of the noise floor value which changes according to the optical path length difference, determines the estimated coherence time so that the measured amplitude value and the calculated amplitude value match, and measures the oscillation linewidth of the coherent light source based on the estimated coherence time (see Patent Document 1 below). [Prior art documents] [Patent Documents]

[0003] [Patent Document 1] Japanese Patent Publication No. 2022-138828 [Non-patent literature]

[0004] [Non-Patent Document 1] Xi Zhang, Fei Yang, Zitong Feng, Fang Wei, Haiwen Cai, and Ronghui Qu, "Narrow-linewidth swept laser phase reconstruction and noise measurement technology and its applications," Opt. Express 26, 32958-32970 (2018) [Overview of the project] [Problems that the invention aims to solve]

[0005] For example, light sources used for fundus measurement require high coherence to widen the measurement range. Tunable VCSELs have recently attracted attention as light sources with long coherence lengths. However, in conventional coherent distance measurement methods, the method of observing the roll-off of the PSF (Point Spread Function) (a method different from the method described in Patent Document 1 above) has a high repetition frequency of wavelength sweep, making it difficult to evaluate Tunable VCSELs that operate in a very high frequency band, partly due to limitations in the frequency band of the detector in the measurement system.

[0006] The method described in Patent Document 1 (referred to as the "noise floor method" for convenience) is a coherence length measurement method that utilizes the behavior of the noise floor of the PSF of an interference signal, which fluctuates with respect to the difference in optical path length. In this method, since only signals in the low-frequency range are used for measurement, the measurement distance and sweep speed are not electrically limited, making it suitable for measuring the coherence length of high-coherence light sources. However, there is also a demand for a method that directly observes the optical electric field output from the light source rather than from noise.

[0007] Therefore, this disclosure aims to measure the coherence distance of a coherence light source operating in a very high frequency band. [Means for solving the problem]

[0008] The coherence distance measurement method of this disclosure outputs light to be measured from a coherence wavelength sweep light source, separates the light to be measured into signal light and reference light using a first optical coupler, combines the signal light and reference light using a second optical coupler, and outputs them as a first interference signal, a second interference signal, and a third interference signal, each with a phase difference of 120°, respectively, which are detected by a photodetector, and calculates a correction matrix using the first, second, and third interference signals to correct the first, second, and third interference signals, in which factors related to amplitude changes due to the first and second optical couplers, the photoelectric efficiency of the photodetector, and the phase change due to the second optical coupler change change over time, calculates the phase difference between the signal light and the reference light using the first, second, and third interference signals and the correction matrix, calculates the spectral width of the light using the phase difference, and calculates the coherence distance of the light using the spectral width.

[0009] Furthermore, the coherence distance measurement method of the present disclosure is a coherence distance measurement device comprising: a coherence wavelength sweep light source that outputs light to be measured; a first optical coupler that separates light to be measured into signal light and reference light; a second optical coupler that, after combining the signal light and reference light, outputs a first interference signal, a second interference signal, and a third interference signal, each with a phase difference of 120°; and a photodetector that detects the first interference signal, the second interference signal, and the third interference signal, respectively, wherein the analysis unit analyzes the first interference signal, the second interference signal, Using the third interference signal, a correction matrix is ​​calculated in which factors related to amplitude changes due to the first and second optical couplers, the photoelectric efficiency of the photodetector, and the phase change due to the second optical coupler change over time, correcting the first, second, and third interference signals. Using the first, second, and third interference signals and the correction matrix, the phase difference between the signal light and the reference light is calculated. Using the phase difference, the spectral width of the light is calculated, and using the spectral width, the coherence distance of the light is calculated. [Effects of the Invention]

[0010] According to the present disclosure using the above solution, there is a benefit in measuring the coherence distance of a coherence light source operating in a very high frequency band. [Brief explanation of the drawing]

[0011] [Figure 1] It is a diagram showing a conventional measurement system for measuring the coherence distance. [Figure 2] It is a diagram showing a system configured to measure the phase difference using a 3x3 coupler. [Figure 3] It is a diagram showing a static correction matrix and mathematical formulas for each factor to calculate the phase difference using three interference signals. [Figure 4] It is a diagram showing a correction matrix and mathematical formulas in which each factor changes with time to calculate the phase difference using three interference signals. [Figure 5] It shows a diagram of Hilbert-transforming three interference signals. [Figure 6] It shows a diagram of calculating the phase shift between three interference signals. [Figure 7] It is a flowchart showing the processing flow of the coherence distance measurement method of the present disclosure. [Figure 8] It is a diagram showing the intermediate processing results when calculated using a static correction matrix and when calculated using a time-varying correction matrix. [Figure 9] It is a table showing the calculated spectral width.

Mode for Carrying Out the Invention

[0012] In the present disclosure, the coherence distance (hereinafter, also appropriately referred to as the coherence length) is measured using a 3x3 coupler. For the smooth progress of the understanding of the present disclosure, a method of observing the roll-off of the PSF using a conventionally known 1x2 coupler and 2x2 coupler will also be briefly described.

[0013] Figure 1 shows a conventional measurement system for measuring coherence distance. The measurement device 1 comprises a wavelength-swept light source 10, an optical branching unit 20, a reference unit 70, a measurement unit 30, an interference unit 40, a photodetector 50, and an analysis unit 60. Light emitted from the wavelength-swept light source 10 propagates through an optical guide path 11a, which is made of an optical fiber or the like, and reaches the optical branching unit 20. The optical branching unit 20 splits the light emitted from the wavelength-swept light source 10 into light to be measured and reference light, and is, for example, a 1x2 optical coupler for optical path splitting. The light branched by the optical branching unit 20 propagates to the reference unit 70 via the optical guide path 11b and to the measurement unit 30 via the optical guide path 11c, respectively.

[0014] The reference unit 70 is a system that serves as a reference for the light being measured in the measurement unit 30. Light propagating through the reference unit 70 is emitted from the light guide path 11f via the light circulator 71, shaped into parallel light by an optical element 72 such as a collimating lens, and reflected by a reference mirror 73, for example, a metal mirror surface. The reflected reference light (return light of the reference light) is guided to the interference unit 40 via the light guide path 11d through the optical element 72, which is also the light receiving system of the reference unit 70.

[0015] In the measurement unit 30, the light to be measured is emitted from the light guide path 11g via the light circulator 31, shaped into parallel light by an optical element 32 such as a collimating lens, and, similar to the reference unit 70, the light to be measured is directed toward the moving mirror M by the optical element 32 such as a collimating lens. The moving mirror M can move in one axis direction on the stage 34, thereby creating an optical path length difference of distance Z from the reference unit 70. The light returned from the moving mirror M (light including reflected and scattered light from the moving mirror M) is guided to the interference unit 40 via the optical element 32, which is also the light receiving system of the measurement unit 30, through the light guide path 11e.

[0016] The interference unit 40 is a system that combines the light to be measured, which has been reflected or scattered by each mirror, with the reference light to generate interference light. The reflected light from the measurement unit 30 and the reference unit 70 is coupled by an optical coupling unit 41, which is, for example, a 2x2 optical coupler for optical path combination, and interferes with each other to become interference light, which is then guided to the photodetector unit 50 by optical guide paths 11h and 11i.

[0017] The photodetector 50 is for detecting interference light and is a balanced photodetector (BPD) that outputs a differential using, for example, two photoelectric conversion elements such as two photodiodes. The photodetector 50 outputs an interference signal, which is input to the analysis unit 60 via the filter 51. As shown in Figure 1, the measurement light (signal light) and the reference light arrive at the optical coupling unit 41 at different times τ0 (delay time) due to the difference in optical path length, and the frequency difference is the beat frequency f b This is how it appears. The interference signal is processed by FFT (Fast Fourier Transform) in the analysis unit 60 and output as a calculation result. The FFT converts the interference signal from the time domain to the frequency domain. Note that since a wavelength-swept light source whose wavelength (frequency) changes with time is used as the light source, it is essentially a conversion from the frequency domain to the time domain.

[0018] A window function is applied to the interference signal to obtain a Point Spread Function (PSF). The coherence distance is then measured based on how much the peak value in the PSF (vertical axis: dB) decreases (rolls off) relative to the Depth (horizontal axis: mm).

[0019] However, due to the frequency bandwidth limitations of balanced photodetectors (BPDs) that output differential signals, it is difficult to acquire differential signals in very high frequency ranges, making it challenging to measure coherence distance using the roll-off method in high bandwidths. For example, in tunable VCSELs, the sweep speed can reach 100-400 kHz. For instance, when measuring the coherence length with a sweep wavelength of 90 nm, a sweep speed of 400 kHz, and a depth of 100 mm, the frequency of the interference signal can reach 15 GHz.

[0020] Therefore, the inventors considered a different method for measuring coherence distance, while also referring to Non-Patent Document 1 mentioned above. Non-Patent Document 1 discloses a method for measuring coherence distance by acquiring three interference signals using a 3x3 coupler. Figure 2 shows a system assembled to measure phase difference using a 3x3 coupler. In this figure, for example, there is a light source 100 which is a coherence wavelength sweep light source using a tunable VCSEL, a first optical coupler 110 which is a 2x2 coupler, a signal coil 120, a reference coil 130, and a second optical coupler 140 which is a 3x3 coupler. The optical path is made of optical fiber, and the signal coil 120 and the reference coil 130 are made of optical fiber, not metal, just to be clear. Light emitted from the light source 100 is separated by the first optical coupler 110, passes through the optical path (signal arm) of the signal coil 120 which is the object of measurement, and the optical path (reference arm) of the reference coil 130 with a different optical path length, and is then combined and interfered with by the second optical coupler 140, and each signal is output. Hereinafter, these three signals will also be referred to as three output signals or three interference signals (interference signal 1, interference signal 2, and interference signal 3).

[0021] The second optical coupler 140 is designed so that its three output signals have a phase difference of 120 degrees. Therefore, the three output signals are independent interference signals, each represented as a sinusoidal function of the optical path length difference between the signal arm and the reference arm. The three output signals contain complex electric field vector components representing the amplitude and phase of the electric field as light propagates. In other words, such a system can be described as a 120-degree phase difference interferometer, where the three output signals have a phase difference of 120 degrees. The analysis unit 160 may be a computing terminal or dedicated processing unit equipped with one or more processors.

[0022] At the signal output, photodetectors 150A, 150B, and 150C, which use photodiodes (PDs) or the like, measure the intensity (square of the amplitude) of the electric field component E, and this is detected as the light intensity I. Since the three output signals are each detected as a single interference signal, detection can be performed without bandwidth limitations like in a BPD.

[0023] First, let's briefly explain how the analysis unit 160 calculates the spectral width from these three output signals. From these three output signals, the phase difference of the interference signals is obtained using a correction matrix, which will be described later. Unwrapping is performed on this phase difference. From the phase difference of the unwrapped interference signals, the time derivative of the instantaneous phase of the light is reconstructed. The reconstructed instantaneous phase is differentiated with respect to time and divided by 2π. This gives the instantaneous frequency. Furthermore, a quasi-linear function of the frequency sweep is obtained using polynomial fitting. Then, the residual frequency error is obtained by subtracting the quasi-linear function from the instantaneous frequency. Then, the power spectral density (PSD) of the frequency fluctuation within a specific observation time window is estimated using the residual frequency error. Here, the PSD is not treated as the absolute value of the spectrum, but rather as information representing the frequency fluctuation (shift) centered on the sweep frequency. The spectral width is calculated using the obtained PSD.

[0024] If the spectral width can be calculated, the coherence distance can be calculated. It is known that the coherence distance Lc can be approximately calculated using the following formula. Lc~λ 2 / Δλ λ: center wavelength of the light source, Δλ: spectral width (full width at half maximum FWHM)

[0025] Figure 3 shows the static correction matrices and formulas for each factor used to calculate the phase difference using three interference signals. Here, the phase difference between the signal light and the reference light is as shown in the formula in Figure 3, η(t)=φ(tT)=arctan(X′2(t) / X′1(t)) This is expressed as follows: T is the delay time of the signal light relative to the reference light. X'2(t) and X'1(t) are derived from the three output signals of the 120-degree phase difference interferometer shown in Figure 2.

[0026] I1(t), I2(t), and I2(t) each represent three interference signals, and X'1(t), X'2(t), and X'3(t) represent their orthogonal components. Here, to calculate the phase difference, it is sufficient to obtain two values, X'1(t) and X'2(t).

[0027] The values ​​of η, ζ, and ξ are parameters that constitute a 3x3 matrix, which is a correction matrix used to compensate for various factors in a real 120-degree phase difference interferometer. Specifically, it includes a combination of elements such as amplitude changes due to 2x2 and 3x3 couplers, photoelectric efficiency of the photodetector, and phase changes due to the 3x3 coupler.

[0028] Therefore, in actual measurements, by appropriately setting the above parameters and applying corrections to the three output signals, X'1(t) and X'2(t) can be obtained, and the phase difference between the signal light and the reference light can be calculated from these.

[0029] However, the correction matrix shown in Figure 3, also described in Non-Patent Document 1, has static values ​​for each of η, ζ, and ξ. In other words, the first and second couplers are wavelength-dependent, but they do not correspond to the change in wavelength in the wavelength-swept light source. Therefore, the amplitude due to the first and second optical couplers, the photoelectric efficiency of the photodetector, and the phase due to the second optical coupler, as described above, have properties that change over time, and accuracy cannot be high unless they change over time.

[0030] Furthermore, while photoelectric efficiency includes factors related to the coupling efficiency between, for example, an optical fiber and an optical coupler, this coupling efficiency is an extremely unreproducible value that changes each time operations such as reconnecting the coupler and fiber are performed.

[0031] For the reasons stated above, in order to improve the accuracy of the final calculation of spectral linewidth and coherence distance, the correction matrix needs to be one in which factors related to amplitude changes due to the first and second optical couplers, the photoelectric efficiency of the photodetector, and the phase change due to the second optical coupler change over time.

[0032] FIG. 4 shows equations and the like for explaining a correction matrix (hereinafter, simply referred to as a time-varying correction matrix or a time-varying matrix for convenience) in which factors related to the amplitude change by the first optical coupler and the second optical coupler, the photoelectric efficiency of the photodetector, and the phase change by the second optical coupler change with time using three interference signals.

[0033] Specifically explain how each parameter included in the time-varying correction matrix relates to the above factors. As shown in this FIG. 4, η n (t)=2c 1n (t)c 2n (t)r n (t)cosφ n (t) ζ n (t)=2c 1n (t)c 2n (t)r n (t)sinφ n (t) ξ n (t)=r n (t)(c 1n (t) 2 +c 2n (t) 2 ) where c 1n (t) is a value obtained by multiplying the amplitude change by outputting from the first optical coupler 110 to the signal coil 120 (output tap on the side) and the amplitude change of the signal output from the nth (from 1 to 3) output tap of the second optical coupler 140 input from the signal coil 120, c 2n (t) is a value obtained by multiplying the amplitude change by outputting to the reference coil 130 (output tap on the side) of the first coupler 110 and the amplitude change of the signal output from the nth (from 1 to 3) output tap of the second optical coupler 140 input from the reference coil 130, r n (t) is the photoelectric efficiency, cosφ n [[ID=,52]](t) and sinφ n (t) are functions related to time of the phase change by the second coupler. Therefore, if these can be obtained from the three output interference signals, the time-varying correction matrix can be obtained.

[0034] Here, within the wavelength band where the accuracy of the 3x3 coupler is guaranteed, the phase difference of the interference signals of each output is guaranteed to be 120°. Therefore, the distribution ratio of all couplers is set to the same value of 1 / 6 (c 1n =c 2n (=1 / 6). However, since the wavelength of a wavelength-swept light source changes with time, in wavelength bands where the accuracy of the coupler is not strictly guaranteed, the distribution ratio and phase difference of the coupler change with time. This point is compensated for by including it in the photoelectric efficiency of the photodetector.

[0035] Figure 5 shows a graph obtained by Hilbert transforming three interference signals to extract factors related to amplitude changes and photoelectric efficiency due to the optical coupler. This allows the amplitude components of each interference signal to be derived. The difference in amplitude of the three interference signals is due to the time dependence of the photoelectric efficiency of each photodetector and the distribution ratio of the coupler, and the temporal change of this amplitude component can be substituted as the photoelectric efficiencies r1(t), r2(t), and r3(t) of each photodetector.

[0036] Figure 6 shows a graph of the phase shift of each interference signal calculated to extract the factors related to the phase change caused by the optical coupler. In other words, each interference signal should have a phase difference of 120° (±2π / 3), but if the difference is greater than that, it is not the designed phase difference, but rather a phase shift. Here, using interference signal 2 as the reference, if we set φ2(t)=0, we can substitute the phase shift between interference signal 2 and interference signal 1, and the phase shift between the phase of interference signal 2 and interference signal 3 shown in Figure 6, as φ1(t) and φ3(t), respectively.

[0037] Based on the above, the time-varying correction matrix shown in Figure 4 can be obtained from the three interference signals. Below, we will summarize the method for measuring the coherence distance using this correction matrix.

[0038] Figure 7 is a flowchart showing the processing flow of the coherence distance measurement method of this disclosure.

[0039] First, in step S101, wavelength-swept light is output from the light source 100 as the light to be measured.

[0040] Next, in step S102, the light to be measured is split into signal light and reference light by the first optical coupler 110.

[0041] Next, in step S103, the signal light and the reference light are combined by the second optical coupler 140 and output as a first interference signal, a second interference signal, and a third interference signal, each with a phase difference of 120°.

[0042] Next, in step S104, the three interference signals are detected and acquired by the photodetectors 150A, 150B, and 150C, respectively.

[0043] Next, in step S105, the analysis unit 160 uses the first interference signal, the second interference signal, and the third interference signal to calculate a correction matrix that corrects the first interference signal, the second interference signal, and the third interference signal, in which factors related to the amplitude change due to the first optical coupler 110 and the second optical coupler 140, the photoelectric efficiency of the photodetectors 150A, 150B, and 150C, and the phase change due to the second optical coupler 140 change over time.

[0044] Next, in step S106, the analysis unit 160 calculates the phase difference between the signal light and the reference light using the first interference signal, the second interference signal, the third interference signal, and the correction matrix.

[0045] Next, in step S107, the analysis unit 160 calculates the spectral width of the wavelength-swept light source using the phase difference. First, an unwrapping process is performed on the phase difference. The time derivative of the instantaneous phase of the light is reconstructed from the phase difference of the unwrapped interference signal. The reconstructed instantaneous phase is differentiated with respect to time and divided by 2π. This gives the instantaneous frequency.

[0046] Furthermore, a quasi-linear function for the frequency sweep is obtained using polynomial fitting. The residual frequency error is then obtained by subtracting the quasi-linear function from the instantaneous frequency. The power spectral density (PSD) of the frequency variation within a specific observation time window is then estimated using the residual frequency error. Here, the PSD is not treated as the absolute value of the spectrum, but rather as information representing the frequency fluctuation (shift) centered on the sweep frequency. The spectral width is calculated using the obtained PSD.

[0047] Figure 8 shows an intermediate processing result that provides a visually easy-to-understand comparison between calculations using a static correction matrix and those using a time-varying correction matrix. In this figure, the electric field, the SFFT (Short-time Fourier Transform) of the electric field, and the PSD are shown for each of the two cases described above. It can be visually understood that the oscillation spectral width is more uniform over the entire wavelength sweep time, i.e., the entire wavelength sweep range, when calculated using the time-varying correction matrix, indicating improved measurement accuracy. Furthermore, it can be seen that the PSD is close to the shape of an ideal Lorentz laser.

[0048] Figure 9 is a table comparing the spectral widths measured and calculated for a coherence light source swept at sweep frequencies of 100 to 400 kHz, using a static correction matrix, a time-varying correction matrix, and the noise floor method, which is also introduced as an invention described in Patent Document 1. In this table, assuming that the accuracy of the spectral width calculation values ​​using the noise floor method is high, it can be seen that, at all sweep frequencies, the values ​​calculated using the time-varying correction matrix have higher measurement accuracy than the values ​​calculated using the static correction matrix.

[0049] Next, in step S108, the coherence distance is calculated from the spectral width.

[0050] As described above, the coherence distance measurement method of this disclosure can be used to measure the coherence distance of a light source. This makes it possible to measure the coherence distance of a coherence light source operating in a very high frequency band.

[0051] This concludes the explanation provided in this disclosure. However, the new technologies described herein can be realized in various other forms, and parts of the content may be omitted, modified, or replaced without departing from the spirit of this disclosure. The embodiments and variations thereof shown in this disclosure are also included in the scope and spirit of this disclosure and shall be treated as equivalent and comparable to the technologies protected under the claims. [Explanation of Symbols]

[0052] 1. Measuring device 10 wavelength swept light source 11a~i Light guide path 20 Optical branching section 30 Measuring part 32 Optical Circumference 40 Interference section 41 Optical coupling part 50 Light detection unit 60 Analysis Department 70 Reference section 71 Light Circulator 72 Optical elements 73 Reference Mirror M Mobile Mirror 2. Coherence Distance Measuring Device 100 light sources 110 First Optical Coupler 120 signal coils 130 Reference coil 140 Second Optical Coupler 150A, 150B, 150C photodetectors 160 Analysis Department

Claims

1. A method for measuring coherence distance, The light to be measured is output from a coherence wavelength-swept light source. The light to be measured is separated into a signal light and a reference light using the first optical coupler. The signal light and the reference light are combined by the second optical coupler and output as a first interference signal, a second interference signal, and a third interference signal, each with a phase difference of 120°, respectively, and each is detected by a photodetector. Using the first interference signal, the second interference signal, and the third interference signal, a correction matrix is ​​calculated in which factors related to the amplitude change due to the first and second optical couplers, the photoelectric efficiency of the photodetector, and the phase change due to the second optical coupler change over time, thereby correcting the first, second, and third interference signals. Using the first interference signal, the second interference signal, the third interference signal, and the correction matrix, the phase difference between the signal light and the reference light is calculated. Using the phase difference, the spectral width of the light is calculated. Using the spectral width, the coherence distance of the light is calculated. Method for measuring coherence distance.

2. A coherent distance measuring device, A coherence wavelength-swept light source that outputs the light to be measured, A first optical coupler that separates the light to be measured into a signal light and a reference light, A second optical coupler combines the signal light and the reference light and outputs them as a first interference signal, a second interference signal, and a third interference signal, each with a phase difference of 120°, respectively. The system comprises a photodetector that detects the first interference signal, the second interference signal, and the third interference signal, respectively. The analysis unit, Using the first interference signal, the second interference signal, and the third interference signal, a correction matrix is ​​calculated in which factors related to the amplitude change due to the first and second optical couplers, the photoelectric efficiency of the photodetector, and the phase change due to the second optical coupler change over time, thereby correcting the first, second, and third interference signals. Using the first interference signal, the second interference signal, the third interference signal, and the correction matrix, the phase difference between the signal light and the reference light is calculated. Using the phase difference, the spectral width of the light is calculated. Using the spectral width, the coherence distance of the light is calculated. Coherence distance measuring device.

Images