[go: up one dir, main page]

JP7352962B2 - Brillouin frequency shift measurement device and Brillouin frequency shift measurement method - Google Patents

Brillouin frequency shift measurement device and Brillouin frequency shift measurement method Download PDF

Info

Publication number
JP7352962B2
JP7352962B2 JP2020024536A JP2020024536A JP7352962B2 JP 7352962 B2 JP7352962 B2 JP 7352962B2 JP 2020024536 A JP2020024536 A JP 2020024536A JP 2020024536 A JP2020024536 A JP 2020024536A JP 7352962 B2 JP7352962 B2 JP 7352962B2
Authority
JP
Japan
Prior art keywords
optical
light
frequency comb
optical frequency
brillouin
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
JP2020024536A
Other languages
Japanese (ja)
Other versions
JP2021128131A (en
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.)
Tokyo University of Agriculture and Technology NUC
Original Assignee
Tokyo University of Agriculture and Technology NUC
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 Tokyo University of Agriculture and Technology NUC filed Critical Tokyo University of Agriculture and Technology NUC
Priority to JP2020024536A priority Critical patent/JP7352962B2/en
Publication of JP2021128131A publication Critical patent/JP2021128131A/en
Application granted granted Critical
Publication of JP7352962B2 publication Critical patent/JP7352962B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Optical Devices Or Fibers (AREA)
  • Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)

Description

本発明は、ブリルアン周波数シフト測定装置及びブリルアン周波数シフト測定方法に関する。 The present invention relates to a Brillouin frequency shift measuring device and a Brillouin frequency shift measuring method.

従来から、光ファイバ内のブリルアン散乱を利用した分布型センサの研究が数多くなされている。初期の研究としては、単一モード光源を用いるものが提案されている(非特許文献1)。また、これまでの研究は、非特許文献2にまとめられている。最近、光周波数コムを利用した手法が、周波数掃引なしでブリルアン利得スペクトルを検出できる手法として提案されている(非特許文献3~6)。一方、単一モード光源を用いた手法であるが、ブリルアン利得スペクトルの変化が線形とみなせる領域を利用して、散乱光の強度変化から、温度や歪みを測定する手法が提案されている(非特許文献7)。また、単一モードのポンプ光に変調をかけることで多モードとし、各モードに対応するブリルアン利得スペクトルが重なるようにすることでブリルアン利得スペクトルを整形し、ブリルアン利得スペクトルの線形とみなせる領域を広くするものが提案されている(非特許文献8、9)。 Conventionally, many studies have been conducted on distributed sensors that utilize Brillouin scattering within optical fibers. In early research, a method using a single mode light source was proposed (Non-Patent Document 1). Further, the research conducted so far is summarized in Non-Patent Document 2. Recently, a method using an optical frequency comb has been proposed as a method that can detect the Brillouin gain spectrum without frequency sweeping (Non-Patent Documents 3 to 6). On the other hand, although this method uses a single mode light source, a method has been proposed that uses a region where changes in the Brillouin gain spectrum can be considered linear to measure temperature and strain from changes in the intensity of scattered light (non-linear). Patent Document 7). In addition, by modulating the single mode pump light, it becomes multi-mode, and by making the Brillouin gain spectra corresponding to each mode overlap, the Brillouin gain spectrum is shaped, widening the region where the Brillouin gain spectrum can be considered linear. have been proposed (Non-patent Documents 8 and 9).

T. Horiguchi, K. Shimizu, T. Kurashima, M. Tateda, and Y. Koyamada, “Development of a distributed sensing technique using Brillouin scattering,” J. Lightwave Technol. 13, 1296-1302 (1995).T. Horiguchi, K. Shimizu, T. Kurashima, M. Tateda, and Y. Koyamada, “Development of a distributed sensing technique using Brillouin scattering,” J. Lightwave Technol. 13, 1296-1302 (1995). K. Y. Song, K. Hotate, W. Zou, and Z. He, “Applications of Brillouin dynamic grating to distributed fiber sensors,” J. Lightwave Technol. 35, 3268-3280 (2017).K. Y. Song, K. Hotate, W. Zou, and Z. He, “Applications of Brillouin dynamic grating to distributed fiber sensors,” J. Lightwave Technol. 35, 3268-3280 (2017). A. Voskoboinik, J.Wang, B. Shamee, S. R. Nuccio, L. Zhang, M Chitgarha, A. E. Willner, and M. Tur, “SBS-based fiber optical sensing using frequency-domain simultaneous tone interrogation,” J. Lightwave Technol. 29, 1729-1735 (2011).A. Voskoboinik, J.Wang, B. Shamee, S. R. Nuccio, L. Zhang, M Chitgarha, A. E. Willner, and M. Tur, “SBS-based fiber optical sensing using frequency-domain simultaneous tone interrogation,” J. Lightwave Technol 29, 1729-1735 (2011). C. Jin, L. Wang, Y. Chen, N. Guo, W. Chung, H. Au, Z. Li, H-Y. Tam, and C. Lu, “Single-measurement digital optical frequency comb based phase-detection Brillouin optical time domain analyzer,” Opt. Exp. 25, 9213-9224 (2017).C. Jin, L. Wang, Y. Chen, N. Guo, W. Chung, H. Au, Z. Li, H-Y. Tam, and C. Lu, “Single-measurement digital optical frequency comb based phase-detection Brillouin optical time domain analyzer,” Opt. Exp. 25, 9213-9224 (2017). Y. Tanaka and Y. Ozaki, “Brillouin frequency shift measurement with virtually controlled sensitivity,” Appl. Phys. Exp. 10, 062504 (2017).Y. Tanaka and Y. Ozaki, “Brillouin frequency shift measurement with virtually controlled sensitivity,” Appl. Phys. Exp. 10, 062504 (2017). Y. Tanaka, Y. Ozaki, and Y. So, “Scanless Brillouin gain spectrum measurement based on multi-heterodyne detection,” in Tech. Digest of International Conf. on Optical Fiber Sensors 2018, TuE88 (2018).Y. Tanaka, Y. Ozaki, and Y. So, “Scanless Brillouin gain spectrum measurement based on multi-heterodyne detection,” in Tech. Digest of International Conf. on Optical Fiber Sensors 2018, TuE88 (2018). H. Lee, N. Hayashi, Y. Mizuno, and K. Nakamura, “Slope-assisted Brillouin optical correlation-domain reflectometry: proof of concept,” Photon. Jour. 8, 6802807 (2016).H. Lee, N. Hayashi, Y. Mizuno, and K. Nakamura, “Slope-assisted Brillouin optical correlation-domain reflectometry: proof of concept,” Photon. Jour. 8, 6802807 (2016). Guangyao Yang, Xinyu Fan, Bin Wang, and Zuyuan He, "Enhancing strain dynamic range of slope-assisted BOTDA by manipulating Brillouin gain spectrum shape," Opt. Express 26, 32599-32607 (2018)Guangyao Yang, Xinyu Fan, Bin Wang, and Zuyuan He, "Enhancing strain dynamic range of slope-assisted BOTDA by manipulating Brillouin gain spectrum shape," Opt. Express 26, 32599-32607 (2018) J. Marinelarena, J. Urricelqui, and A. Loayssa, IEEE Photonics J. 9,6802710 (2017).J. Marinelarena, J. Urricelqui, and A. Loayssa, IEEE Photonics J. 9,6802710 (2017).

従来の単一モード光源を利用する手法において、プローブ光の周波数掃引を用いる手法では、本質的に測定時間がかかり、また、光源の周波数掃引には複雑な制御が必要となる。また、散乱光とポンプ光(基準となる光)との周波数差を検出する場合、散乱光が微弱で周波数にも揺らぎがある。そのため、大量のデータをとって平均化を行う必要があり、やはり測定時間を要する。 In the conventional method using a single mode light source, the method using the frequency sweep of the probe light inherently requires a long measurement time, and the frequency sweep of the light source requires complicated control. Furthermore, when detecting the frequency difference between the scattered light and the pump light (reference light), the scattered light is weak and has fluctuations in frequency. Therefore, it is necessary to collect and average a large amount of data, which also requires measurement time.

光周波数コムを用いてブリルアン利得スペクトルの一括検出を行う従来の手法では、周波数掃引が不要となるものの、スペクトルを再現する際の信号処理が複雑になる。 Although the conventional method of detecting the Brillouin gain spectrum all at once using an optical frequency comb eliminates the need for frequency sweeping, it complicates signal processing when reproducing the spectrum.

単一モード光源を使用し、光源の変調等によって変形したブリルアン利得スペクトルについて、その変化が周波数軸上で線形とみなせる領域を利用する手法は、測定の短時間化には有効であるが、強度変化と温度・歪み等の変化とを対応付ける事前のキャリブレーションが複雑になる。また、ポンプ光の各モードが作るブリルアン利得スペクトルを重ねることによってブリルアン利得スペクトルを整形する手法は、測定の短時間化には有効であるが、変調の調整は複雑になる。 A method that uses a single mode light source and uses a region where the change in the Brillouin gain spectrum can be considered linear on the frequency axis, which is modified by light source modulation, etc., is effective for shortening the measurement time, but the intensity Preliminary calibration that associates changes with changes in temperature, strain, etc. becomes complicated. Further, a method of shaping the Brillouin gain spectrum by overlapping the Brillouin gain spectra created by each mode of the pump light is effective for shortening the measurement time, but it complicates modulation adjustment.

本発明は、以上のような課題に鑑みてなされたものであり、その目的とするところは、測定時間を短縮化しつつ信号処理を簡便化することが可能なブリルアン周波数シフト測定装置等を提供することにある。 The present invention has been made in view of the above-mentioned problems, and its purpose is to provide a Brillouin frequency shift measuring device etc. that can shorten measurement time and simplify signal processing. There is a particular thing.

(1)本発明は、レーザー光源からの光を2つに分岐する分岐部と、分岐された2つの光のいずれかの光の周波数をシフトする光周波数シフト部と、分岐された一方の光から第1の光周波数コムを発生する第1光周波数コム発生部と、分岐された他方の光から、前記第1の光周波数コムと周波数間隔が異なる第2の光周波数コムを発生する第2光周波数コム発生部と、前記第1の光周波数コム及び/又は前記第2の光周波数コムのスペクトル強度を整形するスペクトル整形部と、前記第1の光周波数コムをパルス化してプローブ光を生成する第1光変調部と、前記第2の光周波数コムをパルス化してポンプ光を生成する第2光変調部と、測定対象である光ファイバの一端側から前記プローブ光を入射し前記光ファイバの他端側から前記ポンプ光を入射したときに前記光ファイバの他端側から出射される光を検出する光検出部と、前記光検出部で検出された光強度に基づいて、前記光ファイバ中のブリルアン周波数シフトを求める制御部とを含むブリルアン周波数シフト測定装置に関する。 (1) The present invention includes a branching section that branches light from a laser light source into two, an optical frequency shift section that shifts the frequency of one of the two branched lights, and an optical frequency shift section that shifts the frequency of one of the two branched lights. a first optical frequency comb generating section that generates a first optical frequency comb from the optical frequency comb; and a second optical frequency comb generating section that generates a second optical frequency comb having a frequency interval different from that of the first optical frequency comb from the other branched light. an optical frequency comb generator, a spectrum shaping unit that shapes the spectral intensity of the first optical frequency comb and/or the second optical frequency comb, and generates a probe light by pulsing the first optical frequency comb. a first optical modulator that pulses the second optical frequency comb to generate pump light; and a second optical modulator that pulses the second optical frequency comb to generate pump light; a photodetector that detects light emitted from the other end of the optical fiber when the pump light is input from the other end; The present invention relates to a Brillouin frequency shift measuring device including a control section for determining a Brillouin frequency shift in a Brillouin frequency shift.

また本発明は、レーザー光源からの光を2つに分岐する分岐ステップと、分岐された2つの光のいずれかの光の周波数をシフトする光周波数シフトステップと、分岐された一方の光から第1の光周波数コムを発生する第1光周波数コム発生ステップと、分岐された他方の光から、前記第1の光周波数コムと周波数間隔が異なる第2の光周波数コムを発生する第2光周波数コム発生ステップと、前記第1の光周波数コム及び/又は前記第2の光周波数コムのスペクトル強度を整形するスペクトル整形ステップと、前記第1の光周波数コムをパルス化してプローブ光を生成する第1光変調ステップと、前記第2の光周波数コムをパルス化してポンプ光を生成する第2光変調ステップと、測定対象である光ファイバの一端側から前記プローブ光を入射し前記光ファイバの他端側から前記ポンプ光を入射したときに前記光ファイバの他端側から出射される光を検出する光検出ステップと、前記光検出ステップで検出された光強度に基づいて、前記光ファイバ中のブリルアン周波数シフトを求める制御ステップとを含むブリルアン周波数シフト測定方法に関する。 Further, the present invention includes a branching step of branching light from a laser light source into two, an optical frequency shifting step of shifting the frequency of one of the two branched lights, and a step of shifting the frequency of one of the two branched lights. a first optical frequency comb generating step of generating one optical frequency comb; and a second optical frequency comb generating step of generating a second optical frequency comb having a frequency interval different from that of the first optical frequency comb from the other branched light. a comb generation step, a spectral shaping step of shaping the spectral intensity of the first optical frequency comb and/or the second optical frequency comb, and a step of pulsing the first optical frequency comb to generate probe light. a second optical modulation step of pulsing the second optical frequency comb to generate pump light; and a second optical modulation step of pulsing the second optical frequency comb to generate pump light, and inputting the probe light from one end side of the optical fiber to be measured, a light detection step of detecting light emitted from the other end of the optical fiber when the pump light is input from the end; and a control step for determining a Brillouin frequency shift.

本発明によれば、光周波数コムを用いることで、測定精度を劣化させずに測定時間を短縮化することができる。また、本発明では、光ファイバからの出射光から周波数情報を取
り出すのではなく、出射光の強度からブリルアン周波数シフトを求めるため、信号処理を簡便化することができる。
According to the present invention, by using an optical frequency comb, measurement time can be shortened without deteriorating measurement accuracy. Further, in the present invention, the Brillouin frequency shift is determined from the intensity of the emitted light rather than extracting frequency information from the emitted light from the optical fiber, so signal processing can be simplified.

(2)また本発明に係るブリルアン周波数シフト測定装置及びブリルアン周波数シフト測定方法では、前記スペクトル整形部は(前記スペクトル整形ステップでは)、前記ブリルアン周波数シフトと前記光検出部で検出される光強度との関係が所定の関数に当てはまるように、前記第1の光周波数コム及び/又は前記第2の光周波数コムのスペクトル強度を整形してもよい。 (2) Furthermore, in the Brillouin frequency shift measuring device and the Brillouin frequency shift measuring method according to the present invention, the spectrum shaping section (in the spectrum shaping step) adjusts the Brillouin frequency shift and the light intensity detected by the light detection section. The spectral intensities of the first optical frequency comb and/or the second optical frequency comb may be shaped so that the relationship of the above applies to a predetermined function.

本実施形態に係るブリルアン周波数シフト測定装置の構成を示す図。FIG. 1 is a diagram showing the configuration of a Brillouin frequency shift measuring device according to the present embodiment. プローブ光、ポンプ光、ブリルアン利得スペクトル及び増幅されたプローブ光の一例を示す図。FIG. 3 is a diagram showing an example of probe light, pump light, Brillouin gain spectrum, and amplified probe light. ブリルアン周波数シフトと増幅されたプローブ光の強度との関係を示す図。FIG. 3 is a diagram showing the relationship between the Brillouin frequency shift and the intensity of amplified probe light. 合成利得係数が一定の周波数範囲で定数になるときの例を示す図。FIG. 7 is a diagram illustrating an example when a composite gain coefficient becomes a constant in a certain frequency range. 任意の関数を直線とした場合のブリルアン周波数シフトと増幅されたプローブ光の強度との関係をシミュレーションにより求めた結果を示す図。The figure which shows the result calculated|required by the simulation of the relationship between Brillouin frequency shift and the intensity of amplified probe light when an arbitrary function is made into a straight line. 本実施形態の手法による効果を確認する確認実験の実験系を示す図。FIG. 3 is a diagram showing an experimental system for a confirmation experiment to confirm the effect of the method of the present embodiment. 30mの被測定ファイバの端5mの部分を暖めたときの、被測定ファイバにおける温度分布の変化を観測した結果を示す図。The figure which shows the result of observing the change of the temperature distribution in the measured fiber when the end 5 m of the 30 m long measured fiber is heated. ポンプ光とプローブ光との相対周波数差を変化させ、増幅されたプローブ光の強度を測定した結果を示す図。The figure which shows the result of changing the relative frequency difference of pump light and probe light, and measuring the intensity of amplified probe light.

以下、本実施形態について説明する。なお、以下に説明する本実施形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また本実施形態で説明される構成の全てが、本発明の必須構成要件であるとは限らない。 This embodiment will be described below. Note that this embodiment described below does not unduly limit the content of the present invention described in the claims. Furthermore, not all of the configurations described in this embodiment are essential components of the present invention.

図1は、本発明の実施形態に係る測定装置(ブリルアン周波数シフト測定装置)の構成を示す図である。測定装置1は、光ファイバ2のセンシング領域SA(測定対象)でのブリルアン周波数シフトを測定する装置であり、単一波長のレーザー光源10と、光分岐器11(分岐部)と、光周波数シフタ20(光周波数シフト部)と、光周波数コム発生器30,31(第1光周波数コム発生部、第2光周波数コム発生部)と、波長フィルタ40,41(スペクトル整形部)と、光変調器50,51(第1光変調部、第2光変調部)と、信号発生器52と、光サーキュレータ60と、光検出器70(光検出部)と、制御装置80(制御部)とを含む。 FIG. 1 is a diagram showing the configuration of a measuring device (Brillouin frequency shift measuring device) according to an embodiment of the present invention. The measuring device 1 is a device that measures the Brillouin frequency shift in the sensing area SA (measurement target) of the optical fiber 2, and includes a single wavelength laser light source 10, an optical splitter 11 (branching part), and an optical frequency shifter. 20 (optical frequency shift section), optical frequency comb generators 30 and 31 (first optical frequency comb generation section, second optical frequency comb generation section), wavelength filters 40 and 41 (spectrum shaping section), and optical modulation devices 50 and 51 (first optical modulation section, second optical modulation section), signal generator 52, optical circulator 60, photodetector 70 (photodetection section), and control device 80 (control section). include.

光分岐器11は、レーザー光源10からの光(周波数ν)を、所定の分岐比(例えば、1対1)となるように2つに分岐する。 The optical splitter 11 branches the light (frequency ν 0 ) from the laser light source 10 into two at a predetermined splitting ratio (for example, 1:1).

光周波数シフタ20は、光分岐器11で分岐された2つの光のいずれかの光の周波数をシフトする。図1に示す例では、分岐された2つの光のうち光周波数コム発生器31に入射する光の周波数をシフトしている。光周波数シフタ20による周波数シフト量をνとする。 The optical frequency shifter 20 shifts the frequency of one of the two lights split by the optical splitter 11. In the example shown in FIG. 1, the frequency of the light incident on the optical frequency comb generator 31 of the two split lights is shifted. Let the amount of frequency shift by the optical frequency shifter 20 be ν s .

光周波数コム発生器30は、分岐された一方の光から光周波数コム(第1の光周波数コム)を発生する。光周波数コムは、周波数軸上に等間隔に並んだ成分(モード)からなる櫛形のスペクトルを持つ光である。光周波数コム発生器31は、分岐された他方の光(光周波数シフタ20から入射する光)から、第1の光周波数コムと周波数間隔(スペクトルにおいて隣り合うモードの周波数間隔)が僅かに異なる光周波数コム(第2の光周波数コ
ム)を発生する。
The optical frequency comb generator 30 generates an optical frequency comb (first optical frequency comb) from one of the branched lights. An optical frequency comb is light with a comb-shaped spectrum consisting of components (modes) arranged at equal intervals on the frequency axis. The optical frequency comb generator 31 generates light whose frequency interval (frequency interval between adjacent modes in the spectrum) is slightly different from that of the first optical frequency comb from the other branched light (the light incident from the optical frequency shifter 20). A frequency comb (second optical frequency comb) is generated.

波長フィルタ40は、第1の光周波数コムのスペクトル強度(各モードの強度)を整形する。波長フィルタ41は、第2の光周波数コムのスペクトル強度を整形する。図1に示す例では、波長フィルタ40は、第1の光周波数コムのスペクトルにおける各モードの強度が互いに異なるように第1の光周波数コムのスペクトル強度を整形し、波長フィルタ41は、第2の光周波数コムのスペクトルにおける各モードの強度が等しくなるように第2の光周波数コムのスペクトル強度を整形している。 The wavelength filter 40 shapes the spectral intensity (intensity of each mode) of the first optical frequency comb. The wavelength filter 41 shapes the spectral intensity of the second optical frequency comb. In the example shown in FIG. 1, the wavelength filter 40 shapes the spectral intensity of the first optical frequency comb so that the intensities of each mode in the spectrum of the first optical frequency comb are different from each other, and the wavelength filter 41 shapes the spectral intensity of the second optical frequency comb. The spectral intensity of the second optical frequency comb is shaped so that the intensity of each mode in the spectrum of the second optical frequency comb is equal.

光変調器50は、整形された第1の光周波数コムに一定のタイミングで時間ゲートをかけて(周期的なゲートパルスによる強度変調をかけて)パルス化してプローブ光Prを生成する。プローブ光Prは、光ファイバ2のセンシング領域SAの一端側から入射する。光変調器51は、整形された第2の光周波数コムに一定のタイミングで時間ゲートをかけてパルス化してポンプ光Pmを生成する。信号発生器52は、光変調器50,51を駆動するための信号(強度変調信号)を発生する。 The optical modulator 50 generates probe light Pr by applying a time gate to the shaped first optical frequency comb at a constant timing (applying intensity modulation using a periodic gate pulse). The probe light Pr enters the sensing area SA of the optical fiber 2 from one end side. The optical modulator 51 applies a time gate to the shaped second optical frequency comb at a constant timing to generate a pulse, thereby generating pump light Pm. The signal generator 52 generates a signal (intensity modulation signal) for driving the optical modulators 50 and 51.

光サーキュレータ60は、第1のポートに入力するポンプ光Pmを、第2のポートから出力してセンシング領域SAの他端側に導くとともに、センシング領域SAの他端側から出射され第2のポートに入力する光(増幅されたプローブ光AP)を、第3のポートから出力して光検出器70に導く。光検出器70は、増幅されたプローブ光APを検出し、検出した光強度を電気信号として出力する。 The optical circulator 60 outputs the pump light Pm input to the first port from the second port and guides it to the other end of the sensing area SA, and also outputs the pump light Pm from the other end of the sensing area SA to the second port. The light (amplified probe light AP) input to the probe is output from the third port and guided to the photodetector 70. The photodetector 70 detects the amplified probe light AP and outputs the detected light intensity as an electrical signal.

制御装置80は、プロセッサ(CPU等)と記憶部(RAM、ハードディスク等)とを備えた情報処理装置(例えば、パーソナルコンピュータ)であり、光検出器70からの信号の強度(光検出器70で検出された光強度)に基づいて、光ファイバ2中(センシング領域SA)のブリルアン周波数シフト(ブリルアン利得スペクトルの周波数シフト)を算出する。また、制御装置80は、信号発生器52に制御信号を出力して信号発生器52を制御する。 The control device 80 is an information processing device (for example, a personal computer) equipped with a processor (such as a CPU) and a storage unit (RAM, hard disk, etc.), and controls the intensity of the signal from the photodetector 70 (the intensity of the signal from the photodetector 70). Based on the detected light intensity), the Brillouin frequency shift (frequency shift of the Brillouin gain spectrum) in the optical fiber 2 (sensing area SA) is calculated. Further, the control device 80 outputs a control signal to the signal generator 52 to control the signal generator 52.

光ファイバ2中では、熱的に誘起された音響波により屈折率分布(動く回折格子)が生じる。光ファイバ2(センシング領域SA)の他端側に入射したポンプ光Pmがこれにより反射され、周波数がダウンシフトしたストローク光が発生し後方に(光ファイバの他端側に)伝搬する。更に、ストローク光とポンプ光Pmのビートにより電歪が生じ、屈折率分布が生じる。この一連の過程が繰り返され、ブリルアン利得スペクトル(BGS:Brillouin Gain Spectrum)が生じる。プローブ光Prがポンプ光Pmと対向して(光ファイバの一端側に)入射すると、ブリルアン利得スペクトルによりプローブ光Prが増幅され、光ファイバの他端側から増幅されたプローブ光APが出射される。 In the optical fiber 2, a refractive index distribution (moving diffraction grating) is generated by thermally induced acoustic waves. The pump light Pm incident on the other end of the optical fiber 2 (sensing area SA) is reflected thereby, and a stroke light whose frequency is downshifted is generated and propagates backward (toward the other end of the optical fiber). Furthermore, electrostriction occurs due to the beats of the stroke light and the pump light Pm, resulting in a refractive index distribution. This series of processes is repeated to generate a Brillouin Gain Spectrum (BGS). When the probe light Pr enters opposite the pump light Pm (to one end of the optical fiber), the probe light Pr is amplified by the Brillouin gain spectrum, and the amplified probe light AP is emitted from the other end of the optical fiber. .

図2に、プローブ光Pr、ポンプ光Pm、ブリルアン利得スペクトル(BGS)及び増幅されたプローブ光APの一例を示す。図中の上向き矢印は、光周波数コムの各モードを示している。また、図中のBGS(n=1~4)は、ポンプ光Pmのn番目のモードにより発生したブリルアン利得スペクトルを示している。 FIG. 2 shows an example of the probe light Pr, the pump light Pm, the Brillouin gain spectrum (BGS), and the amplified probe light AP. The upward arrows in the figure indicate each mode of the optical frequency comb. Further, BGS n (n=1 to 4) in the figure indicates the Brillouin gain spectrum generated by the n-th mode of the pump light Pm.

本実施形態の測定装置1では、プローブ光Prである光周波数コムの周波数間隔Δνprobeとポンプ光Pmである光周波数コムの周波数間隔Δνpumpとが、僅かに異なるようにする。光周波数コムのモード数が多いほど、周波数間隔Δνprobeと周波数間隔Δνpumpとの差が小さくなるようにする。また、ポンプ光Pmの周波数をプローブ光Prの周波数に対してシフトする。光ファイバ2中では、ポンプ光Pmの光強度が一定以上のとき、ポンプ光Pmの周波数よりも低い一定の周波数範囲にブリルアン利得スペクトルが生じ、その周波数範囲において、ポンプ光Pmと対向伝搬する光(プローブ光P
r)を増幅する。例えば、ポンプ光Pmの周波数が波長1.5μmのとき、約11GHz低い周波数にブリルアン利得スペクトルが生じる。ポンプ光Pmの1つのモードの周波数(スペクトルピークの中心周波数)と、当該モードによって発生したブリルアン利得スペクトルの中心周波数との周波数差を、ブリルアン周波数シフト(BFS:Brillouin Frequency Shift)と呼ぶ。図2に示すように、プローブ光Prの各モードは、ポンプ光Pmの各モードの周波数よりも低い周波数に生じた各ブリルアン利得スペクトル(BGS~BGS)により増幅されている。なお、光周波数コムの周波数間隔Δνpump,Δνprobeは、ブリルアン周波数シフトよりも広くなるようにしている。
In the measuring device 1 of this embodiment, the frequency interval Δν probe of the optical frequency comb that is the probe light Pr is made to be slightly different from the frequency interval Δν pump of the optical frequency comb that is the pump light Pm. The larger the number of modes of the optical frequency comb, the smaller the difference between the frequency interval Δν probe and the frequency interval Δν pump . Furthermore, the frequency of the pump light Pm is shifted relative to the frequency of the probe light Pr. In the optical fiber 2, when the light intensity of the pump light Pm is above a certain level, a Brillouin gain spectrum occurs in a certain frequency range lower than the frequency of the pump light Pm, and in that frequency range, light propagating counter-propagating with the pump light Pm (Probe light P
r) amplify. For example, when the frequency of the pump light Pm has a wavelength of 1.5 μm, a Brillouin gain spectrum occurs at a frequency approximately 11 GHz lower. The frequency difference between the frequency of one mode of the pump light Pm (the center frequency of the spectrum peak) and the center frequency of the Brillouin gain spectrum generated by the mode is called a Brillouin frequency shift (BFS). As shown in FIG. 2, each mode of the probe light Pr is amplified by each Brillouin gain spectrum (BGS 1 to BGS 4 ) generated at a frequency lower than the frequency of each mode of the pump light Pm. Note that the frequency interval Δν pump and Δν probe of the optical frequency comb is set to be wider than the Brillouin frequency shift.

ブリルアン利得スペクトルの中心周波数は、光ファイバ2に加わる熱や歪みに比例してシフトする。従って、ブリルアン周波数シフトを求めることで、光ファイバ2の温度や歪みを測定することができる。本実施形態の測定装置1では、増幅されたプローブ光APの強度(全パワー)を光検出器70で検出する。プローブ光Prとなる光周波数コムとポンプ光Pmとなる光周波数コムのスペクトル強度を予め適切に設定してスペクトル整形部(波長フィルタ40,41)により整形することで、例えば図3に示すように、広い周波数範囲にわたって、ブリルアン周波数シフト(BFS)に対して増幅されたプローブ光APの検出強度が比例するようになる。すなわち、スペクトル整形部(波長フィルタ40,41)は、ブリルアン周波数シフトと光検出器70で検出される光強度との関係が、所定の関数に当てはまる(例えば、単調増加や単調減少の関係となる)ように、光周波数コムのスペクトル強度を整形する。これにより、光検出器70で検出された光強度から簡単な計算で一意にブリルアン周波数シフトを求めることができる。 The center frequency of the Brillouin gain spectrum shifts in proportion to the heat and strain applied to the optical fiber 2. Therefore, by determining the Brillouin frequency shift, the temperature and strain of the optical fiber 2 can be measured. In the measuring device 1 of this embodiment, the intensity (total power) of the amplified probe light AP is detected by the photodetector 70. By appropriately setting the spectral intensities of the optical frequency comb that becomes the probe light Pr and the optical frequency comb that becomes the pump light Pm in advance and shaping them by the spectrum shaping section (wavelength filters 40 and 41), for example, as shown in FIG. , the detection intensity of the amplified probe light AP becomes proportional to the Brillouin frequency shift (BFS) over a wide frequency range. That is, the spectrum shaping section (wavelength filters 40, 41) determines whether the relationship between the Brillouin frequency shift and the light intensity detected by the photodetector 70 corresponds to a predetermined function (for example, a monotonically increasing or monotonically decreasing relationship). ) to shape the spectral intensity of the optical frequency comb. Thereby, the Brillouin frequency shift can be uniquely determined from the light intensity detected by the photodetector 70 by simple calculation.

また、プローブ光Prの増幅は、ポンプ光Pmとプローブ光Prが衝突する点においてのみ起こる。そして、パルス光であるポンプ光Pmとプローブ光Prの衝突点は、光変調器50,51によるゲートのタイミング(パルスの出射タイミング)で制御することができる。従って、信号発生器52を制御して、センシング領域SAにおける任意の位置に衝突点を設定することで、任意の位置での増幅されたプローブ光APの強度を取得することができる。そして、センシング領域SAにおいて複数の位置に衝突点を設定し、設定した衝突点ごとの増幅されたプローブ光APの強度を取得し、設定した衝突点ごとのブリルアン周波数シフトを求めることで、センシング領域SAにおける熱や歪みの分布を測定することができる。 Further, amplification of the probe light Pr occurs only at the point where the pump light Pm and the probe light Pr collide. The collision point of the pump light Pm, which is pulsed light, and the probe light Pr can be controlled by the gate timing (pulse emission timing) by the optical modulators 50 and 51. Therefore, by controlling the signal generator 52 and setting the collision point at an arbitrary position in the sensing area SA, it is possible to obtain the intensity of the amplified probe light AP at an arbitrary position. Then, by setting collision points at multiple positions in the sensing area SA, obtaining the intensity of the amplified probe light AP for each set collision point, and finding the Brillouin frequency shift for each set collision point, the sensing area The distribution of heat and strain in SA can be measured.

次に、光周波数コムのスペクトル強度の整形手法の一例について説明する。N本のモードからなるポンプ光Pmについては、波長フィルタ41により各モードの強度を等しくする。また、プローブ光PrもN本のモードからなるとして、プローブ光Prの周波数の低い方からk番目のモードの強度をpとする。ここで、ブリルアン利得スペクトルにより増幅されたプローブ光APの強度Pap(νBFS)は、以下の式(1)により表される。 Next, an example of a method for shaping the spectral intensity of an optical frequency comb will be described. Regarding the pump light Pm consisting of N modes, the wavelength filter 41 equalizes the intensity of each mode. Further, assuming that the probe light Pr is also composed of N modes, the intensity of the k-th mode from the lowest frequency of the probe light Pr is defined as pk . Here, the intensity P apBFS ) of the probe light AP amplified by the Brillouin gain spectrum is expressed by the following equation (1).

Figure 0007352962000001
ここで、G(νBFS,νk)は、k番目のモードに対する利得係数であり、ブリルアン周波数シフトνBFSと、k番目のモードの強度pの関数となっている。ここで、合成利得係数G(νBFS)を、以下の式(2)で定義する。
Figure 0007352962000001
Here, G 0BFS , νk) is a gain coefficient for the kth mode and is a function of the Brillouin frequency shift ν BFS and the intensity p k of the kth mode. Here, the composite gain coefficient G SBFS ) is defined by the following equation (2).

Figure 0007352962000002
合成利得係数G(νBFS)は、増幅されたプローブ光APの強度Pap(νBFS)のブリルアン周波数シフトνBFSに対する比である。k番目のモードに対する利得係数G(νBFS,νk)が予め測定されており、且つプローブ光Prの強度が大きくなくプローブ光Prの影響を受けないとすれば、合成利得係数G(νBFS)が一定の周波数範囲で定数になる(すなわち、νBFSとPap(νBFS)との関係が一定の周波数範囲で比例関係となる)ような各モードの強度Pは、数値計算で容易に求めることができる。図4は、合成利得係数G(νBFS)が一定の周波数範囲で定数になるときの例であり、図中点線の各曲線は、増幅されたプローブ光APの各モードの強度とブリルアン周波数シフト(BFS)との関係を示し、図中実線の直線は、増幅されたプローブ光APの強度(全パワー)とブリルアン周波数シフトとの関係を示している。
Figure 0007352962000002
The composite gain coefficient G SBFS ) is the ratio of the intensity P apBFS ) of the amplified probe light AP to the Brillouin frequency shift ν BFS . If the gain coefficient G 0BFS , νk) for the k-th mode has been measured in advance and the intensity of the probe light Pr is not large and is not affected by the probe light Pr, then the combined gain coefficient G S (ν The strength P k of each mode such that BFS ) becomes constant in a certain frequency range (that is, the relationship between ν BFS and P apBFS ) is proportional in a certain frequency range) can be calculated by numerical calculation. can be easily determined. FIG. 4 is an example when the composite gain coefficient G SBFS ) becomes a constant in a certain frequency range, and each dotted curve in the figure shows the intensity of each mode of the amplified probe light AP and the Brillouin frequency. The solid straight line in the figure shows the relationship between the intensity (total power) of the amplified probe light AP and the Brillouin frequency shift.

N本のモードからなるポンプ光Pmによって生じた各ブリルアン利得スペクトルは、例えばローレンツ関数で与えられる。そのとき、ブリルアン利得スペクトル全体GBFS(ν)は、以下の式(3)により表される。 Each Brillouin gain spectrum generated by the pump light Pm consisting of N modes is given by a Lorentz function, for example. At that time, the entire Brillouin gain spectrum G BFS (v) is expressed by the following equation (3).

Figure 0007352962000003
ここで、gは、ポンプ光Pmのk番目のモードによって生成したブリルアン利得スペクトルの最大値を決める定数であり、特に、ポンプ光Pmの各モードの強度が等しいとき、g=g(一定)である。また、νpump,kは、ポンプ光Pmのk番目のモードの周波数であり、Δνpumpは、ポンプ光Pmのスペクトル間隔(隣り合うモードの周波数間隔)であり、Δνは、定数である。
Figure 0007352962000003
Here, g k is a constant that determines the maximum value of the Brillouin gain spectrum generated by the k-th mode of pump light Pm. In particular, when the intensity of each mode of pump light Pm is equal, g k = g 0 ( constant). Further, ν pump,k is the frequency of the k-th mode of the pump light Pm, Δν pump is the spectral interval (frequency interval between adjacent modes) of the pump light Pm, and Δν is a constant.

また、プローブ光PrのスペクトルPprobe(ν)は、以下の式(4)により表される。 Further, the spectrum P probe (ν) of the probe light Pr is expressed by the following equation (4).

Figure 0007352962000004
ここで、νprobe,kは、プローブ光Prのk番目のモードの周波数であり、Δνprobeは、プローブ光Prのスペクトル間隔である。また、fは、プローブ光Prの各モードのスペクトル形状を与える関数であり、f(νprobe,k)=1とみなしてよい。
Figure 0007352962000004
Here, ν probe,k is the frequency of the k-th mode of the probe light Pr, and Δν probe is the spectral interval of the probe light Pr. Further, f is a function that gives the spectral shape of each mode of the probe light Pr, and may be regarded as f(ν probe ,k)=1.

ここでは、ポンプ光Pmの各モードの強度が等しいとする。また、ポンプ光Pm、プローブ光Prのスペクトル線幅は十分に狭いものとする(数学的には、近似的にデルタ関数列と考える)。ブリルアン利得スペクトルにより増幅されるプローブ光APの各モードの強度は、pBFS(νprobe,k)に比例する。従って、ブリルアン利得スペク
トルにより増幅されるプローブ光APの強度Pap(νBFS)は、以下の式(5)で与えられる。
Here, it is assumed that the intensity of each mode of the pump light Pm is equal. Further, the spectral line widths of the pump light Pm and the probe light Pr are assumed to be sufficiently narrow (mathematically, they are approximately considered as a delta function sequence). The intensity of each mode of the probe light AP amplified by the Brillouin gain spectrum is proportional to p k G BFS (v probe,k ). Therefore, the intensity P apBFS ) of the probe light AP amplified by the Brillouin gain spectrum is given by the following equation (5).

Figure 0007352962000005
ここで、Aは、比例係数である。Pap(νBFS)は、ブリルアン周波数シフトνBFSの関数となっており、プローブ光Prの各モードの強度pを適切に設定することで、任意の関数に近似した関数を合成する(νBFSとPap(νBFS)との関係が所定の関数に当てはまるようにする)ことができる。なお、比例係数Aは、事前に一定のポンプ光(単一波長)に対して、プローブ光(単一波長)の強度とブリルアン利得スペクトルにより増幅されるプローブ光の強度との関係を測定しておくことで求めることができる。
Figure 0007352962000005
Here, A is a proportionality coefficient. P apBFS ) is a function of the Brillouin frequency shift ν BFS , and by appropriately setting the intensity p k of each mode of the probe light Pr, a function that approximates an arbitrary function can be synthesized (ν The relationship between BFS and P ap (v BFS ) can be made to fit a predetermined function). Note that the proportionality coefficient A is determined by measuring the relationship between the intensity of the probe light (single wavelength) and the intensity of the probe light amplified by the Brillouin gain spectrum with respect to a fixed pump light (single wavelength) in advance. You can find it by setting it.

例えば、任意の関数として直線を合成すれば、ブリルアン利得スペクトルにより増幅されるプローブ光の強度Pap(νBFS)が、光ファイバの温度や歪みによるブリルアン周波数シフトνBFSの変化に比例した出力となるようにすることができる。任意の関数を直線とする場合、所望の直線(Pap(νBFS)=aνBFS+b)に対して、誤差関数e(νBFS)を以下の式(6)で定義する。 For example, if straight lines are synthesized as arbitrary functions, the intensity P apBFS ) of the probe light amplified by the Brillouin gain spectrum becomes an output proportional to the change in the Brillouin frequency shift ν BFS due to the temperature and strain of the optical fiber. You can make it happen. When an arbitrary function is a straight line, an error function e(ν BFS ) is defined by the following equation (6) for a desired straight line (P apBFS )=aν BFS +b).

Figure 0007352962000006
そして、e(νBFS)=0とおき、この関数に対して最小二乗法によりプローブ光Prの各モードの強度pを決定する。図5に、任意の関数を直線とした場合のブリルアン周波数シフトνBFSと増幅されるプローブ光の強度Pap(νBFS)との関係をシミュレーションにより求めた結果を示す。
Figure 0007352962000006
Then, e(ν BFS )=0 is set, and the intensity p k of each mode of the probe light Pr is determined using the least squares method for this function. FIG. 5 shows the results obtained by simulation of the relationship between the Brillouin frequency shift ν BFS and the intensity P apBFS ) of the amplified probe light when an arbitrary function is a straight line.

図6は、本実施形態の手法による効果を確認する確認実験の実験系を示す図である。図6において、LDは、単一波長レーザー光源であり、SSは、スペクトル整形器であり、IMは、光強度変調器であり、PSは、偏波スクランブラであり、BPFは、光バンドパスフィルタであり、PDは、光検出器であり、LIAは、ロックインアンプであり、FUTは、被測定ファイバ(Fiber Under Test、センシング領域SAに相当)である。 FIG. 6 is a diagram showing an experimental system for a confirmation experiment to confirm the effect of the method of this embodiment. In FIG. 6, LD is a single wavelength laser light source, SS is a spectrum shaper, IM is an optical intensity modulator, PS is a polarization scrambler, and BPF is an optical bandpass PD is a photodetector, LIA is a lock-in amplifier, and FUT is a fiber under test (corresponding to sensing area SA).

単一波長レーザー光源LDの波長は、1553nmとした。実験系の上側の光路は、プローブ光Prを生成する光路であり、光源からの光に、光強度変調器IMにより12.02GHzで深い変調をかけ、周波数間隔(Δνprobe)が12.02GHzの5モードからなる光周波数コム(第1の光周波数コム)を生成した。すなわち、この光強度変調器IMは、第1光周波数コム発生部として機能する。 The wavelength of the single wavelength laser light source LD was 1553 nm. The upper optical path of the experimental system is the optical path that generates the probe light Pr, in which the light from the light source is deeply modulated at 12.02 GHz by the optical intensity modulator IM, and the frequency interval (Δν probe ) is 12.02 GHz. An optical frequency comb (first optical frequency comb) consisting of five modes was generated. That is, this optical intensity modulator IM functions as a first optical frequency comb generator.

実験系の下側の光路は、ポンプ光Pmを生成する光路であり、光源からの光に、最初の光強度変調器IMにより10.86GHzで浅い変調をかけ、周波数軸上で元の光周波数スペクトルの両側に生じる側波帯のうち、高周波側のスペクトルを光バンドパスフィルタBPFで切り出した。これに2番目の光強度変調器IMにより12.00GHzで深い変調をかけ、周波数間隔(Δνpump)が12.00GHzの5モードからなる光周波数
コム(第2の光周波数コム)を生成した。すなわち、最初の光強度変調器IM、光バンドパスフィルタBPF及び2番目の光強度変調器IMは、光周波数シフト部及び第2光周波数コム発生部として機能する。
The lower optical path of the experimental system is the optical path that generates the pump light Pm, and the light from the light source is shallowly modulated at 10.86 GHz by the first optical intensity modulator IM, so that the original optical frequency on the frequency axis is Among the sidebands occurring on both sides of the spectrum, the spectrum on the high frequency side was cut out using an optical bandpass filter BPF. This was deeply modulated at 12.00 GHz by a second optical intensity modulator IM to generate an optical frequency comb (second optical frequency comb) consisting of five modes with a frequency interval (Δν pump ) of 12.00 GHz. That is, the first optical intensity modulator IM, the optical bandpass filter BPF, and the second optical intensity modulator IM function as an optical frequency shifter and a second optical frequency comb generator.

また、実験系の上側の光路及び下側の光路のそれぞれにおいて、スペクトル整形器SSにより、増幅されるプローブ光APの強度がブリルアン周波数シフトに比例するように光周波数コムを整形した。ここでは、スペクトル整形器SSとして、入射光を分光して、異なる波長の光を空間上に並べ、それぞれの光の強度を可変減衰器で調整したのち、再び合波するもの(例えば、導波路型のスペクトル整形器)を使用した。なお、スペクトル整形器SSとして、空間型液晶光変調器を用いたものを使用することもできる。また、光強度変調器IMによって、周期的なゲートパルスによりそれぞれの光周波数コムに強度変調をかけてパルス化して、プローブ光Pr、ポンプ光Pmを生成した。このゲートパルスは、空間に変換すると幅3m(15ns)、周期30m(150ns)であり、空間分解能3m、測定レンジ30mを与える。 Furthermore, in each of the upper and lower optical paths of the experimental system, the optical frequency combs were shaped by a spectrum shaper SS so that the intensity of the amplified probe light AP was proportional to the Brillouin frequency shift. Here, the spectrum shaper SS is one that splits the incident light, arranges the lights of different wavelengths in space, adjusts the intensity of each light with a variable attenuator, and then combines them again (for example, a waveguide A type of spectral shaper) was used. Note that as the spectrum shaper SS, one using a spatial liquid crystal light modulator can also be used. In addition, each optical frequency comb was intensity-modulated by a periodic gate pulse using an optical intensity modulator IM and pulsed to generate probe light Pr and pump light Pm. This gate pulse has a width of 3 m (15 ns) and a period of 30 m (150 ns) when converted into space, giving a spatial resolution of 3 m and a measurement range of 30 m.

30mの被測定ファイバFUTの端5mの部分(図中(b)で示す部分)を暖め、被測定ファイバFUTの各測定点(パルスの衝突点)での増幅されたプローブ光APの強度を時系列で取得して各測定点でのブリルアン周波数シフトの変化を測定し、被測定ファイバFUTにおける温度分布の変化を観測した。結果を図7に示す。図7のグラフは、被測定ファイバFUTの端から3mずつ温度(ブリルアン周波数シフト(BFS))の変化を測定した結果であり、上から、0~3m区間の結果、3~6m区間の結果、6~30m区間の結果を示す。被測定ファイバFUTの暖めた部分に近い区間ほどブリルアン周波数シフトが大きく変化しており、暖めた部分から遠い区間ではブリルアン周波数シフトが変化していない。この結果から、光ファイバの温度分布の変化を精度良く測定できることが確認された。 The 5 m end portion (the part shown in (b) in the figure) of the 30 m long fiber under test FUT is heated, and the intensity of the amplified probe light AP at each measurement point (pulse collision point) of the fiber under test FUT is measured over time. A change in the Brillouin frequency shift at each measurement point was measured by acquiring the data in series, and a change in the temperature distribution in the fiber FUT to be measured was observed. The results are shown in FIG. The graph in FIG. 7 shows the results of measuring the change in temperature (Brillouin frequency shift (BFS)) from the end of the fiber under test FUT by 3 m. From the top, the results for the 0 to 3 m section, the results for the 3 to 6 m section, The results for the 6-30m section are shown. The Brillouin frequency shift changes more greatly in the section closer to the heated part of the fiber under test FUT, and the Brillouin frequency shift does not change in the section farther from the heated part. From this result, it was confirmed that changes in the temperature distribution of the optical fiber can be measured with high accuracy.

また、本実験系において、ポンプ光Pmとプローブ光Prとの相対周波数差(ブリルアン周波数シフトと等価)を変化させ、増幅されたプローブ光APの強度を測定した。結果を図8に示す(図中黒塗り点)。この結果から、100MHz(100度の温度変化に相当)の広い周波数範囲にわたり、相対周波数差の変化に対して増幅されたプローブ光APの強度が線形に変化していることが確認された。 Furthermore, in this experimental system, the relative frequency difference (equivalent to Brillouin frequency shift) between the pump light Pm and the probe light Pr was varied, and the intensity of the amplified probe light AP was measured. The results are shown in Figure 8 (black dots in the figure). From this result, it was confirmed that over a wide frequency range of 100 MHz (corresponding to a temperature change of 100 degrees), the intensity of the amplified probe light AP changes linearly with respect to changes in the relative frequency difference.

本実施形態のブリルアン周波数シフト測定装置によれば、光周波数コムを用いることで、測定精度を劣化させずに測定時間を短縮化することができる。また、光ファイバからの出射光から周波数情報を取り出すのではなく、出射光の強度からブリルアン周波数シフトを求めるため、信号処理を簡便化することができる。また、光周波数コムのスペクトル強度(各モードの強度)を予め設計したフィルタで整形すればよいため、キャリブレーションが容易である。 According to the Brillouin frequency shift measurement device of this embodiment, by using an optical frequency comb, measurement time can be shortened without deteriorating measurement accuracy. Further, since the Brillouin frequency shift is determined from the intensity of the emitted light rather than extracting frequency information from the emitted light from the optical fiber, signal processing can be simplified. Further, since the spectral intensity (intensity of each mode) of the optical frequency comb can be shaped using a pre-designed filter, calibration is easy.

なお、本発明は、上述の実施の形態に限定されるものではなく、種々の変更が可能である。本発明は、実施の形態で説明した構成と実質的に同一の構成(例えば、機能、方法及び結果が同一の構成、あるいは目的及び効果が同一の構成)を含む。また、本発明は、実施の形態で説明した構成の本質的でない部分を置き換えた構成を含む。また、本発明は、実施の形態で説明した構成と同一の作用効果を奏する構成又は同一の目的を達成することができる構成を含む。また、本発明は、実施の形態で説明した構成に公知技術を付加した構成を含む。 Note that the present invention is not limited to the above-described embodiments, and various modifications are possible. The present invention includes configurations that are substantially the same as those described in the embodiments (for example, configurations that have the same functions, methods, and results, or configurations that have the same objectives and effects). Further, the present invention includes a configuration in which non-essential parts of the configuration described in the embodiments are replaced. Further, the present invention includes a configuration that has the same effects or a configuration that can achieve the same objective as the configuration described in the embodiment. Further, the present invention includes a configuration in which a known technique is added to the configuration described in the embodiment.

例えば、上記実施形態では、プローブ光となる光周波数コムのスペクトル強度と、ポンプ光となる光周波数コムのスペクトル強度とを整形する場合について説明したが、プローブ光となる光周波数コムのスペクトル強度のみを整形するようにしてもよいし、ポンプ光
となる光周波数コムのスペクトル強度のみを整形するようにしてもよい。
For example, in the above embodiment, a case has been described in which the spectral intensity of the optical frequency comb that becomes the probe light and the spectral intensity of the optical frequency comb that becomes the pump light are shaped, but only the spectral intensity of the optical frequency comb that becomes the probe light Alternatively, only the spectral intensity of the optical frequency comb serving as the pump light may be shaped.

1…測定装置(ブリルアン周波数シフト測定装置)、2…光ファイバ、10…レーザー光源、11…光分岐器、20…光周波数シフタ、30,31…光周波数コム発生器、40,41…波長フィルタ、50,51…光変調器、52…信号発生器、60…光サーキュレータ、70…光検出器、80…制御装置 1... Measuring device (Brillouin frequency shift measuring device), 2... Optical fiber, 10... Laser light source, 11... Optical splitter, 20... Optical frequency shifter, 30, 31... Optical frequency comb generator, 40, 41... Wavelength filter , 50, 51... Optical modulator, 52... Signal generator, 60... Optical circulator, 70... Photodetector, 80... Control device

Claims (2)

レーザー光源からの光を2つに分岐する分岐部と、
分岐された2つの光のいずれかの光の周波数をシフトする光周波数シフト部と、
分岐された一方の光から第1の光周波数コムを発生する第1光周波数コム発生部と、
分岐された他方の光から、前記第1の光周波数コムと周波数間隔が異なる第2の光周波数コムを発生する第2光周波数コム発生部と、
前記第1の光周波数コム及び/又は前記第2の光周波数コムのスペクトル強度を整形するスペクトル整形部と、
前記第1の光周波数コムをパルス化してプローブ光を生成する第1光変調部と、
前記第2の光周波数コムをパルス化してポンプ光を生成する第2光変調部と、
測定対象である光ファイバの一端側から前記プローブ光を入射し前記光ファイバの他端側から前記ポンプ光を入射したときに前記光ファイバの他端側から出射される光を検出する光検出部と、
前記光検出部で検出された光強度に基づいて、前記光ファイバ中のブリルアン周波数シフトを求める制御部とを含み、
前記スペクトル整形部は、
前記ブリルアン周波数シフトと前記光検出部で検出される光強度との関係が所定の関数に当てはまるように、前記第1の光周波数コム及び/又は前記第2の光周波数コムのスペクトル強度を整形する、ブリルアン周波数シフト測定装置。
A branching part that branches the light from the laser light source into two,
an optical frequency shifter that shifts the frequency of one of the two branched lights;
a first optical frequency comb generator that generates a first optical frequency comb from one of the branched lights;
a second optical frequency comb generator that generates a second optical frequency comb having a frequency interval different from that of the first optical frequency comb from the other branched light;
a spectrum shaping unit that shapes the spectral intensity of the first optical frequency comb and/or the second optical frequency comb;
a first optical modulator that pulses the first optical frequency comb to generate probe light;
a second optical modulator that pulses the second optical frequency comb to generate pump light;
a light detection unit that detects light emitted from the other end of the optical fiber when the probe light is input from one end of the optical fiber to be measured and the pump light is input from the other end of the optical fiber; and,
a control unit that determines a Brillouin frequency shift in the optical fiber based on the light intensity detected by the light detection unit ;
The spectrum shaping section is
Shaping the spectral intensity of the first optical frequency comb and/or the second optical frequency comb so that the relationship between the Brillouin frequency shift and the light intensity detected by the photodetector fits a predetermined function. , Brillouin frequency shift measurement device.
レーザー光源からの光を2つに分岐する分岐ステップと、
分岐された2つの光のいずれかの光の周波数をシフトする光周波数シフトステップと、
分岐された一方の光から第1の光周波数コムを発生する第1光周波数コム発生ステップと、
分岐された他方の光から、前記第1の光周波数コムと周波数間隔が異なる第2の光周波数コムを発生する第2光周波数コム発生ステップと、
前記第1の光周波数コム及び/又は前記第2の光周波数コムのスペクトル強度を整形するスペクトル整形ステップと、
前記第1の光周波数コムをパルス化してプローブ光を生成する第1光変調ステップと、
前記第2の光周波数コムをパルス化してポンプ光を生成する第2光変調ステップと、
測定対象である光ファイバの一端側から前記プローブ光を入射し前記光ファイバの他端側から前記ポンプ光を入射したときに前記光ファイバの他端側から出射される光を検出する光検出ステップと、
前記光検出ステップで検出された光強度に基づいて、前記光ファイバ中のブリルアン周波数シフトを求める制御ステップとを含み、
前記スペクトル整形ステップでは、
前記ブリルアン周波数シフトと前記光検出ステップで検出される光強度との関係が所定の関数に当てはまるように、前記第1の光周波数コム及び/又は前記第2の光周波数コムのスペクトル強度を整形する、ブリルアン周波数シフト測定方法。
a branching step of branching the light from the laser light source into two;
an optical frequency shifting step of shifting the frequency of one of the two branched lights;
a first optical frequency comb generation step of generating a first optical frequency comb from one of the branched lights;
a second optical frequency comb generation step of generating a second optical frequency comb having a different frequency interval from the first optical frequency comb from the other branched light;
a spectral shaping step of shaping the spectral intensity of the first optical frequency comb and/or the second optical frequency comb;
a first optical modulation step of pulsing the first optical frequency comb to generate probe light;
a second optical modulation step of pulsing the second optical frequency comb to generate pump light;
a photodetection step of detecting light emitted from the other end of the optical fiber when the probe light is input from one end of the optical fiber to be measured and the pump light is input from the other end of the optical fiber; and,
a control step of determining a Brillouin frequency shift in the optical fiber based on the light intensity detected in the light detection step ,
In the spectrum shaping step,
Shaping the spectral intensity of the first optical frequency comb and/or the second optical frequency comb so that the relationship between the Brillouin frequency shift and the light intensity detected in the photodetection step fits a predetermined function. , Brillouin frequency shift measurement method.
JP2020024536A 2020-02-17 2020-02-17 Brillouin frequency shift measurement device and Brillouin frequency shift measurement method Active JP7352962B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020024536A JP7352962B2 (en) 2020-02-17 2020-02-17 Brillouin frequency shift measurement device and Brillouin frequency shift measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020024536A JP7352962B2 (en) 2020-02-17 2020-02-17 Brillouin frequency shift measurement device and Brillouin frequency shift measurement method

Publications (2)

Publication Number Publication Date
JP2021128131A JP2021128131A (en) 2021-09-02
JP7352962B2 true JP7352962B2 (en) 2023-09-29

Family

ID=77488669

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020024536A Active JP7352962B2 (en) 2020-02-17 2020-02-17 Brillouin frequency shift measurement device and Brillouin frequency shift measurement method

Country Status (1)

Country Link
JP (1) JP7352962B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114353685B (en) * 2022-01-12 2024-05-28 太原理工大学 High-frequency dynamic strain measurement device and method based on chaotic Brillouin phase spectrum
CN115265613B (en) * 2022-07-28 2024-07-02 浙江大学绍兴研究院 Multi-frequency interval Brillouin frequency shift extraction method and device

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130025374A1 (en) 2011-03-03 2013-01-31 University Of Southern California Sweep-Free Stimulated Brillouin Scattering-Based Fiber Optical Sensing
WO2016080415A1 (en) 2014-11-17 2016-05-26 国立大学法人東京農工大学 Measurement device and sensor system
JP2019203859A (en) 2018-05-25 2019-11-28 日本電信電話株式会社 Device and method for measuring brillouin frequency shift
JP2019203860A (en) 2018-05-25 2019-11-28 日本電信電話株式会社 Device and method for measuring brillouin gain spectrum

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130025374A1 (en) 2011-03-03 2013-01-31 University Of Southern California Sweep-Free Stimulated Brillouin Scattering-Based Fiber Optical Sensing
WO2016080415A1 (en) 2014-11-17 2016-05-26 国立大学法人東京農工大学 Measurement device and sensor system
JP2019203859A (en) 2018-05-25 2019-11-28 日本電信電話株式会社 Device and method for measuring brillouin frequency shift
JP2019203860A (en) 2018-05-25 2019-11-28 日本電信電話株式会社 Device and method for measuring brillouin gain spectrum

Also Published As

Publication number Publication date
JP2021128131A (en) 2021-09-02

Similar Documents

Publication Publication Date Title
JP7115387B2 (en) Light source device for light measurement, spectroscopic measurement device and spectroscopic measurement method
JP5021221B2 (en) Distributed optical fiber sensor
JP7086225B2 (en) Fiber optic BOCDA sensor using phase code modulation of pump light and probe light with time lag
JP7215060B2 (en) Spectroscopic analysis light source, spectroscopic analysis device and spectroscopic analysis method
CN109556527B (en) Optical fiber strain measuring device and optical fiber strain measuring method
JP6868246B2 (en) Equipment and methods for measuring Brillouin frequency shift
JP7352962B2 (en) Brillouin frequency shift measurement device and Brillouin frequency shift measurement method
KR101310783B1 (en) Distributed optical fiber sensor and sensing method using simultaneous sensing of brillouin gain and loss
KR101889351B1 (en) Spatially-selective brillouin distributed optical fiber sensor with increased effective sensing points and sensing method using brillouin scattering
JP3883458B2 (en) Reflective Brillouin spectral distribution measuring method and apparatus
CN103148940A (en) Light asynchronous sampling signal measurement method and system
CN105806374A (en) Fiber bragg grating wavelength demodulation method
CA2392308A1 (en) System and method for resolving polarization mode dispersion in optical fibers
JP7493772B2 (en) Distance measuring device and distance measuring method
JP7396382B2 (en) Optical fiber sensor and Brillouin frequency shift measurement method
JP6602689B2 (en) Optical line characteristic analyzer and signal processing method
JP3905780B2 (en) Brillouin spectral distribution measuring method and apparatus
US6417926B1 (en) Wavelength measuring system
KR20160122319A (en) Fiber optic botda sensor using multiple light sources and method for sensing thereof
Cooper et al. Simple and highly sensitive method for wavelength measurement of low-power time-multiplexed signals using optical amplifiers
CN111537069A (en) Terahertz time-domain spectroscopy and imaging system with distributed temperature measurement function
JP2022039684A (en) Optical fiber strain measuring device and optical fiber strain measurement method
JP7521612B2 (en) Brillouin gain spectrum distribution measuring method and apparatus
CN212340434U (en) Terahertz time-domain spectroscopy and imaging system with distributed temperature measurement function
WO2022181373A1 (en) Measuring device, method for adjusting measuring device, and measuring method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221115

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230731

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230808

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230822

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20230905

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230911

R150 Certificate of patent or registration of utility model

Ref document number: 7352962

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150