Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP3065833B2 - Temperature distribution detector - Google Patents
[go: Go Back, main page]

JP3065833B2 - Temperature distribution detector - Google Patents

Temperature distribution detector

Info

Publication number
JP3065833B2
JP3065833B2 JP5000646A JP64693A JP3065833B2 JP 3065833 B2 JP3065833 B2 JP 3065833B2 JP 5000646 A JP5000646 A JP 5000646A JP 64693 A JP64693 A JP 64693A JP 3065833 B2 JP3065833 B2 JP 3065833B2
Authority
JP
Japan
Prior art keywords
light
optical fiber
signal
scattered light
temperature distribution
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.)
Expired - Fee Related
Application number
JP5000646A
Other languages
Japanese (ja)
Other versions
JPH06201488A (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.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP5000646A priority Critical patent/JP3065833B2/en
Priority to KR94000146A priority patent/KR0133488B1/en
Publication of JPH06201488A publication Critical patent/JPH06201488A/en
Priority to US08/528,594 priority patent/US5639162A/en
Priority to US08/813,980 priority patent/US5825804A/en
Application granted granted Critical
Publication of JP3065833B2 publication Critical patent/JP3065833B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Measuring Temperature Or Quantity Of Heat (AREA)

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【産業上の利用分野】本発明は光ファイバ内の後方散乱
光を利用した温度分布検出装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a temperature distribution detecting device utilizing backscattered light in an optical fiber.

【0002】[0002]

【従来の技術】光ファイバの後方散乱光を利用したセン
サでは、図8に示す如く光ファイバ100に入射したレ
ーザパルスによって光ファイバ100内の散乱光を励起
し、励起散乱光のスペクトルや強度、偏光等の状態に被
散乱場所の温度情報を含み、光ファイバ100内の後方
となる励起レーザパルスの入射側に伝搬する後方散乱光
を時系列信号として検出処理することにより、光ファイ
バ100に沿った温度の一次元分布を測定する。
2. Description of the Related Art In a sensor using backscattered light from an optical fiber, as shown in FIG. 8, a laser pulse incident on the optical fiber 100 excites the scattered light in the optical fiber 100, and the spectrum and intensity of the excitation scattered light are reduced. The state of the polarized light or the like includes the temperature information of the scattered place, and the backscattered light propagating to the incident side of the excitation laser pulse that is located rearward in the optical fiber 100 is detected and processed as a time-series signal. Measure the one-dimensional distribution of temperature.

【0003】この場合、光ファイバ100内にレーザパ
ルスを入射したとき発生する散乱光としては、密度揺ら
ぎによるレーリー散乱光、伝搬性の揺らぎによるブリル
アン散乱光、分子の振動回転によるラマン散乱光があ
る。
In this case, the scattered light generated when a laser pulse enters the optical fiber 100 includes Rayleigh scattered light due to density fluctuation, Brillouin scattered light due to propagation fluctuation, and Raman scattered light due to molecular rotation. .

【0004】このうち、ブリルアン散乱光と、ラマン散
乱光は非弾性散乱光であり、励起光とスペクトルの異な
った散乱光となる。温度情報は3つの散乱光、すべてに
含まれているが、温度に対する感度が最も高いのはラマ
ン散乱光で、その強度が温度に依存して変化する。
[0004] Among them, Brillouin scattered light and Raman scattered light are inelastic scattered light, and become scattered light having a spectrum different from that of excitation light. The temperature information is contained in all three scattered lights, but Raman scattered light has the highest sensitivity to temperature, and its intensity changes depending on the temperature.

【0005】そして、ラマン散乱光を使用して温度を計
測する場合、波長が励起光より長波長側にシフトするス
トークスラマン散乱光と、短波長側にシフトするアンチ
ストークスラマン散乱光とをフィルタで選別し、2つの
散乱光の比を基本とした値から温度分布を算出してい
る。また、比を用いないまでも、フィルタで一方の散乱
光を選別する必要がある。
When measuring temperature using Raman scattered light, Stokes Raman scattered light whose wavelength shifts to a longer wavelength side than excitation light and anti-Stokes Raman scattered light whose wavelength shifts to a shorter wavelength side are filtered. The temperature distribution is calculated from the values based on the ratio of the two scattered lights. Further, even if the ratio is not used, it is necessary to select one scattered light with a filter.

【0006】[0006]

【発明が解決しようとする課題】ところで、このような
分布型のセンサでは、被測定物理量の精度とともに、そ
の位置分解能も重要な量になる。
In such a distributed sensor, not only the accuracy of the physical quantity to be measured, but also its positional resolution becomes an important quantity.

【0007】この場合、分解能は一般的に入射パルス幅
によって決まり、さらに後方散乱信号が離散値系となっ
ているときには、さらにサンプリング周波数によっても
左右される。
In this case, the resolution is generally determined by the incident pulse width, and further depends on the sampling frequency when the backscattered signal is a discrete value system.

【0008】但し、サンプリング周期を高くする方法で
は、サンプリング周期が半値幅に対し、十分に短いとす
るなら、ある時間に計測される後方散乱光は入射パルス
の半値幅にわたる散乱光の合計となるため、測定される
温度は半値幅にわたる平均温度となるため、入射パルス
の幅以上に位置分解能を高くすることはできない。
However, in the method of increasing the sampling period, if the sampling period is sufficiently short with respect to the half width, the backscattered light measured at a certain time is the sum of the scattered light over the half width of the incident pulse. Therefore, the measured temperature is the average temperature over the half-value width, so that the position resolution cannot be increased beyond the width of the incident pulse.

【0009】そこで、位置分解能を高くするために、入
射パルスの幅を狭くする方法が有効であるが、この方法
を実現する方法として、入射パルスの自体の幅を非常に
狭くする方法と、有限幅の光パルスに対する応答を1つ
の変換過程と考えて逆変換を行なってインパルス応答を
求める方法とが考えられる。
In order to increase the position resolution, it is effective to reduce the width of the incident pulse. However, as a method for realizing this method, a method of extremely narrowing the width of the incident pulse itself and a method of finite reduction A method is conceivable in which a response to an optical pulse having a width is considered as one conversion process and an inverse conversion is performed to obtain an impulse response.

【0010】しかしながら、入射パルス自体の幅を狭く
する方法は、発光系の立ち上がり特性、立ち下がり特性
によって決まるパルス幅以上に、光パルスの幅を狭くす
ることができないため、現在の開発されている素子の性
能以上に入射パルスの幅を狭くすることが難しく、位置
分解能を高める方法としては不適切である。
However, the method of reducing the width of the incident pulse itself is currently being developed because the width of the light pulse cannot be reduced beyond the pulse width determined by the rising and falling characteristics of the light emitting system. It is difficult to make the width of the incident pulse narrower than the performance of the element, and it is not suitable as a method for increasing the position resolution.

【0011】これに対し、インパルス応答を求める方法
は位置分解能を高める方法として有効であると考えられ
る。
On the other hand, a method for obtaining an impulse response is considered to be effective as a method for increasing the position resolution.

【0012】そこで、このインパルス応答を求める方法
について、ここで、詳細に検討してみる。
Therefore, the method for obtaining the impulse response will be examined in detail here.

【0013】まず、インパルス応答を求める方法で使用
される変換は線形近似の範囲でコンボルーション過程で
あることから、インパルス応答をh(t) 、入力光パルス
をP(t) とすると、後方散乱光信号g(t) は、次式で表
わすことができる。
First, since the transform used in the method for obtaining the impulse response is a convolution process within the range of linear approximation, if the impulse response is h (t) and the input light pulse is P (t), the backscatter The optical signal g (t) can be expressed by the following equation.

【0014】[0014]

【数1】 そして、インパルス応答h(t) は入力信号となる入力光
パルスP(t) を予め測定しておき、この測定結果に対し
てデコンボルーションを行なうことによって求めること
ができる。
(Equation 1) The impulse response h (t) can be obtained by measuring an input optical pulse P (t) as an input signal in advance and performing deconvolution on the measurement result.

【0015】しかしながら、実際には、実測データに対
し、この逆変換過程を適用すると、測定信号中に含まれ
る雑音成分のため、正しいインパルス応答を求めること
ができないことが多い。この理由としては、次に述べる
ようにすると理解しやすい。
However, in practice, when this inverse transformation process is applied to actually measured data, a correct impulse response cannot often be obtained due to noise components included in the measurement signal. The reason is easy to understand as described below.

【0016】まず、前記(1)式をフーリェ変換して後
方散乱光に雑音成分N( ω) を追加すると、 G( ω) +N( ω) =H( ω) ・P( ω) …(2) となり、逆変換のために、この(2)式の両辺をP(
ω) で除算すると、 H( ω) =G( ω) /P( ω) +N( ω) /P( ω) …(3) となる。
First, when the noise component N (ω) is added to the backscattered light by Fourier-transforming the equation (1), G (ω) + N (ω) = H (ω) · P (ω) (2) ) And both sides of this equation (2) are P (
When divided by (ω), H (ω) = G (ω) / P (ω) + N (ω) / P (ω) (3)

【0017】ここで、問題となるのは、右辺の第2項目
の部分であり、入力光パルス成分P( ω) のスペクトル
は有限であるのに対し、雑音成分N( ω) のスペクトル
がかなりの広域まで伸びているため、周波数の高い領域
(角速度ωが大きな領域)で除算結果が発散してしま
う。
The problem here is that of the second item on the right side. The spectrum of the input light pulse component P (ω) is finite, while the spectrum of the noise component N (ω) is quite large. , The division result diverges in a high frequency region (a region where the angular velocity ω is large).

【0018】このため、この発散によって総合的な位置
測定精度が低下してしまうという問題があった。
For this reason, there has been a problem that the divergence lowers the overall position measurement accuracy.

【0019】本発明は上記の事情に鑑み、逆変換過程に
おいて雑音の影響による解の発散を防止し、これによっ
て光パルスの幅を狭くすることなく総合精度を最も改善
させることができる最適な近似を与えることができると
ともに、光パルス幅を狭くした場合に生じるS/N比の
低下やこれに伴う計測時間の増大、半導体レーザの駆動
回路の困難さを生じることなく、位置分解能を向上させ
てホットスポットなどの異常加熱部などの検出を容易に
することができる温度分布検出装置を提供することを目
的としている。
The present invention has been made in view of the above circumstances, and has an optimal approximation capable of preventing the divergence of a solution due to the influence of noise in an inverse transformation process, thereby maximizing the overall accuracy without reducing the width of an optical pulse. And the position resolution can be improved without lowering the S / N ratio caused by narrowing the optical pulse width, increasing the measurement time associated therewith, and making the semiconductor laser drive circuit difficult. It is an object of the present invention to provide a temperature distribution detecting device capable of easily detecting an abnormally heated portion such as a hot spot.

【0020】[0020]

【課題を解決するための手段】上記の目的を達成するた
めに本発明は、光パルスを発生する光源と、前記光パル
スを光ファイバに入射し、その後方散乱光から2種類の
ラマン散乱光を選択する光学系と、前記ラマン散乱光を
電気信号に変換増幅する光電変換部と、この光電変換部
の出力信号をデジタル信号に変換してアベレージングを
行なうデジタルアベレージャ部と、このデジタルアベレ
ージャ部の処理結果に基づいて前記光ファイバに沿った
温度分布を計算する演算部とを有する温度分布検出装置
において、光ファイバの端面からのフレネル反射光の波
形データに対し、デコンボルーション処理を行なって前
記波形データの波形が所定の形状になるまでの反復回数
を求め、温度測定を行なうとき、2種類のラマン散乱光
の波形データに対し、前記反復回数だけデコンボルーシ
ョン処理を繰り返して光ファイバに沿った温度分布を計
算することを特徴としている。
SUMMARY OF THE INVENTION In order to achieve the above object, the present invention provides a light source for generating an optical pulse, two types of Raman scattered light incident on an optical fiber, An optical system for converting the Raman scattered light into an electric signal, a digital averaging unit for converting an output signal of the photoelectric conversion unit into a digital signal and performing averaging, and a digital averaging unit. And a calculation unit for calculating a temperature distribution along the optical fiber based on the processing result of the optical fiber unit, wherein a deconvolution process is performed on the waveform data of the Fresnel reflected light from the end face of the optical fiber. The number of repetitions until the waveform of the waveform data has a predetermined shape is obtained and the temperature measurement is performed. It is characterized by calculating the temperature distribution along the optical fiber by repeating the repetition number of times deconvolution process.

【0021】[0021]

【作用】上記の構成において、光ファイバの端面からの
フレネル反射光の波形データに対し、デコンボルーショ
ン処理を行なって前記波形データの波形が所定の形状に
なるまでの反復回数を求め、温度測定を行なうとき、2
種類のラマン散乱光の波形データに対し、前記反復回数
だけデコンボルーション処理を繰り返して光ファイバに
沿った温度分布を計算することにより、逆変換過程にお
いて雑音の影響で解が発散してしまうのを防止し、これ
によって光パルスの幅を狭くすることなく総合精度を最
も改善させることができる最適な近似を与え、光パルス
幅を狭くした場合に生じるS/N比の低下やこれに伴う
計測時間の増大、半導体レーザの駆動回路の困難さを生
じることなく、位置分解能を向上させてホットスポット
などの異常加熱部などの検出を容易にする。
In the above arrangement, deconvolution processing is performed on the waveform data of the Fresnel reflected light from the end face of the optical fiber to determine the number of repetitions until the waveform of the waveform data has a predetermined shape, and the temperature is measured. When performing
By calculating the temperature distribution along the optical fiber by repeating the deconvolution process for the number of iterations for the waveform data of the Raman scattered light of the type, the solution diverges due to the influence of noise in the inverse conversion process. , Thereby providing the optimum approximation that can improve the overall accuracy without reducing the width of the light pulse, and reducing the S / N ratio caused by narrowing the light pulse width and the measurement accompanying this. The position resolution is improved and the detection of an abnormally heated portion such as a hot spot is facilitated without increasing the time and making the driving circuit of the semiconductor laser difficult.

【0022】[0022]

【実施例】まず、本発明による温度分布検出装置の具体
的な説明に先だって本発明による温度分布検出装置の測
定原理を説明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Prior to the detailed description of the temperature distribution detecting device according to the present invention, the measurement principle of the temperature distribution detecting device according to the present invention will be described.

【0023】今、温度計測用の光ファイバにおいて、単
位長さ当たりの後方散乱係数をS(z) 、入射光パルスが
受ける減衰率をRf(z) 、後方散乱光が入射端までに受
ける減衰率をRb(z) 、光ファイバ中における光パルス
の群速度をVg とすると、この光ファイバに光パルスを
入射させて得られる後方散乱光のインパルス応答の電力
h(t) は次式で表わすことができる。
Now, in an optical fiber for temperature measurement, the backscattering coefficient per unit length is S (z), the attenuation rate received by the incident light pulse is Rf (z), and the attenuation received by the backscattered light up to the incident end. Assuming that the ratio is Rb (z) and the group velocity of the light pulse in the optical fiber is Vg, the power h (t) of the impulse response of the backscattered light obtained by making the light pulse incident on the optical fiber is expressed by the following equation. be able to.

【0024】 h(t) =Rf(z) ・S(z) ・Rb(z) ・Vg/2 …(4) そして、入射される光パルスPとのコンボルーションと
なるのは、アンチストークスラマン散乱光と、ストーク
スラマン散乱光の各係数Sであるため、逆変換はこれら
アンチストークスラマン散乱光と、ストークスラマン散
乱光の各後方散乱信号に対して行なった後、光ファイバ
の減衰効果を補償するため、これらの各後方散乱信号の
除算を行ない、温度分布を計算する手順となる。
H (t) = Rf (z) · S (z) · Rb (z) · Vg / 2 (4) Then, the convolution with the incident light pulse P is anti-Stokes Raman. Since the coefficients are S of the scattered light and Stokes Raman scattered light, inverse conversion is performed on each of the anti-Stokes Raman scattered light and the backscattered signal of the Stokes Raman scattered light, and then the attenuation effect of the optical fiber is compensated. Therefore, the procedure is to calculate the temperature distribution by dividing each of these backscattered signals.

【0025】しかしながら、実際の装置では、後方散乱
光はフォトダイオードなどで電気信号が変換された後、
A/D変換されて離散値系デジタル信号として扱われる
ため、積分形で表現されたコンボルーションは図2に示
す如く行列形式となる。
However, in an actual device, after the backscattered light is converted into an electric signal by a photodiode or the like,
Since the signal is A / D converted and treated as a discrete-valued digital signal, the convolution expressed in an integral form has a matrix form as shown in FIG.

【0026】そして、光パルス入力はn×nの行列とな
り、時間の推移とともに、入射パルスの位置が移動して
いく対角近傍だけを持つ行列となる。そして、各行の要
素は光パルス波形の実測値が移動したものが代入され
る。
Then, the optical pulse input is an n × n matrix, which is a matrix having only the diagonal vicinity where the position of the incident pulse moves with time. Then, as the element of each row, the one in which the measured value of the optical pulse waveform has moved is substituted.

【0027】この後、インパルス応答hを求めるため、
逆行列を両辺から掛け合わせることになるが、雑音の影
響を監視しながら、インパルス応答hを求めるためには
反復法(Jacobi Gauss-Seidel 法)により解を収束させ
ながら求める方法が適している。
Thereafter, in order to obtain an impulse response h,
The inverse matrix is multiplied from both sides. To find the impulse response h while monitoring the influence of noise, a method of finding the solution while converging the solution by an iterative method (Jacobi Gauss-Seidel method) is suitable.

【0028】この反復法では、反復回数をkにすると、
k+1回目のi番目の解hi は次式で与えられる。
In this iterative method, if the number of repetitions is k,
The (k + 1) th i-th solution hi is given by the following equation.

【0029】[0029]

【数2】 一方、図3に示す如く信号処理装置101と、80℃の
恒温槽102と、1Kmの光ファイバ103とを用意
し、この光ファイバ103の所定部分を2m、5m、1
0m、20m、50mの各長さだけ恒温槽102に入れ
て、その温度を80℃に保ちがながら、信号処理装置1
01によってOTDR法でアンチストークスラマン散乱
光、ストークスラマン散乱光を計測し、これらの比を求
めると、図4に示す波形となる。
(Equation 2) On the other hand, as shown in FIG. 3, a signal processing device 101, a constant temperature bath 102 at 80 ° C., and an optical fiber 103 of 1 km are prepared.
0 m, 20 m, and 50 m are put into the constant temperature bath 102, and the temperature is kept at 80 ° C. while the signal processing device 1
When the anti-Stokes Raman scattered light and the Stokes Raman scattered light are measured by the OTDR method using 01, and their ratios are obtained, the waveform shown in FIG. 4 is obtained.

【0030】この場合、測定に用いた光パルスの半値幅
は100nsec当たり、光ファイバの長さで10mに
相当することから、光ファイバ103のうち、恒温槽1
02に入れられている部分の演算結果が全て同じ値にな
るべきであるが、図4から明らかなように、コンボルー
ションのため、実際には10m以下の5mや2mの部分
で十分な振幅を得ることができない。
In this case, the half width of the light pulse used for the measurement is 100 nsec, which is equivalent to 10 m in the length of the optical fiber.
02 should have the same value, but as is clear from FIG. 4, due to the convolution, in practice, a sufficient amplitude can be obtained at 5 m or 2 m below 10 m. I can't get it.

【0031】そこで、このようなコンボルーションの影
響を除くために、上述した測定動作によって得られたデ
ータに対して、上述した反復法を適用してデコンボルー
ションを行なうことにより、デコンボルーションを行な
う前におけるデータの波形が図5(a)に示す形状のと
き、このデータに対し、反復回数が5回のデコンボルー
ションを行なったとき、図5(b)に示す波形が得ら
れ、反復回数が10回のデコンボルーションを行なった
とき、図5(c)に示す波形が得られる。
Therefore, in order to eliminate the influence of such convolution, deconvolution is performed on the data obtained by the above-described measurement operation by applying the above-mentioned iterative method to deconvolution. When the waveform of the data before the execution is the shape shown in FIG. 5A, when the data is subjected to deconvolution with the number of repetitions of 5 times, the waveform shown in FIG. 5B is obtained. When the deconvolution is performed ten times, a waveform shown in FIG. 5C is obtained.

【0032】この図から明らかなように、デコンボルー
ションにより振幅が不足していた箇所の改善されるとと
もに、立ち上がりの改善されているのが分かる。
As is apparent from this figure, the portion where the amplitude is insufficient is improved by the deconvolution, and the rise is improved.

【0033】しかし、このような改善ととともに、反復
回数の増加とともに、雑音の影響が顕著になってくるの
も分かる。
However, it can be seen that the effect of noise becomes remarkable as the number of repetitions increases with the improvement.

【0034】そこで、総合精度としてコンボルーション
のために振幅が所定値に達しない誤差を系統誤差、雑音
の影響によるバラツキをランダム誤差とし、デコンボル
ーションの反復回数と、系統誤差と、ランダム誤差との
関係を調べるために、加熱長5mの部分に関し、ランダ
ム雑音として加熱直前部のデータを使用して上述したデ
コンボルーション処理を行なった。このとき、図6に示
す表が得られた。
Therefore, an error in which the amplitude does not reach a predetermined value due to convolution is regarded as a systematic error, and a variation due to the influence of noise is regarded as a random error, so that the number of deconvolution repetitions, a systematic error, a random error and In order to investigate the relationship, the above-mentioned deconvolution processing was performed on a portion having a heating length of 5 m using the data immediately before the heating as random noise. At this time, the table shown in FIG. 6 was obtained.

【0035】この図6から明らかなように、反復回数を
増加させていくと、ランダム誤差は一定の比率で増加す
るが、系統誤差は最初、一定の比率で低減し、ある回数
を境に増加に転じ、以後一定の比率で増加する。
As is apparent from FIG. 6, when the number of repetitions is increased, the random error increases at a constant rate, but the systematic error decreases at a constant rate at first, and increases after a certain number of times. And then increase at a constant rate.

【0036】すなわち、系統誤差と、ランダム誤差とを
加算した総合誤差を最も小さくする反復回数が存在す
る。
That is, there is the number of iterations that minimizes the total error obtained by adding the systematic error and the random error.

【0037】そして、最適な反復回数はここで与えた逆
変換基準温度部のように、予め温度分布が正確に分かっ
ている場合には、反復計算毎に総合誤差を計算し、その
値が極値となる回数を求めることができる。
If the temperature distribution is known in advance as in the case of the inverse transformation reference temperature section given here, the optimum number of repetitions is calculated as a total error for each repetition calculation. The number of times that the value is obtained can be obtained.

【0038】但し、このような逆変換方式では、計算回
数とともに、ランダム雑音が増加するため、系統誤差が
ほとんどない場合、例えば温度分布がほぼ一様な部分で
は、ランダム雑音のみが増加することになる。
However, in such an inverse conversion method, since random noise increases with the number of calculations, when there is almost no systematic error, for example, in a portion where the temperature distribution is almost uniform, only random noise increases. Become.

【0039】つまり、この逆変換方式が有利になるの
は、温度変化が空間的に激しく生じる場合、特に局所的
に高温、低温となるような場合の温度分布を計測しよう
とする場合である。
In other words, the advantage of this inverse conversion method is when the temperature change occurs spatially intensely, particularly when the temperature distribution is to be measured when the temperature is locally high or low.

【0040】したがって、このような逆変換方式の使用
対象となるプラントシステムとしては、システムの異常
が局所的な温度上昇として顕在化するようなシステムで
あり、このようなシステムで局所的に温度が上昇したと
き、最適な反復回数を使用することにより、これを効率
良く検出することができる。
Accordingly, a plant system to which such an inverse conversion method is used is a system in which an abnormality of the system is manifested as a local temperature rise, and the temperature is locally increased in such a system. When it rises, it can be detected efficiently by using the optimal number of iterations.

【0041】そして、本発明では、逆変換過程におい
て、最適な反復回数を得るために、図7(a)に示す如
く光ファイバの入射面や終端面からのフレネル反射光を
取り込み、このフレネル反射光に対するデコンボルーシ
ョン処理を繰り返し行なって、図7(b)に示す如くそ
のパルス波形が線状になる反復回数を求めて最適な近似
解を求める最適計算過程を選択し、この最適計算過程を
光ファイバの全領域にわたって適用して逆変換を施し、
これによって光ファイバによる温度分布測定を行なうと
き、雑音の影響を最小限に食い止めながら位置分解能を
大幅に改善させるようにしている。
In the present invention, in order to obtain the optimum number of repetitions in the inverse conversion process, the Fresnel reflected light from the entrance surface or the end surface of the optical fiber is taken in as shown in FIG. By repeatedly performing the deconvolution process on the light, the number of repetitions at which the pulse waveform becomes linear as shown in FIG. 7B is selected, and an optimal calculation process for obtaining an optimal approximate solution is selected. Applied over the entire area of the optical fiber to perform the inverse transformation,
In this way, when measuring the temperature distribution using an optical fiber, the positional resolution is greatly improved while minimizing the influence of noise.

【0042】図1は本発明による温度分布検出装置の一
実施例を示す構成図である。
FIG. 1 is a block diagram showing an embodiment of a temperature distribution detecting device according to the present invention.

【0043】この図に示す温度分布検出装置は、光ファ
イバ1と、信号処理装置2とを備えており、初期化時お
よび定期的に、信号処理装置2によって非常に幅の狭い
光パルスを生成してこれを光ファイバ1に入射させ、こ
の光ファイバ1の終端面からのフレネル反射光を取り込
むとともに、このフレネル反射光に基づき反復法によっ
てデコンボルーションを行なってこのフレネル反射光に
対応する信号の波形が線状となる最適な反復回数を求め
る。そして、通常の測定を行なうとき、信号処理装置2
によって非常に幅の狭い光パルスを生成してこれを光フ
ァイバ1に入射させ、この後方散乱光を取込んで、その
アンチストークスラマン散乱光と、ストークスラマン散
乱光とに基づき前記最適な反復回数に基づいてデコンボ
ルーションを行なって光ファイバ1が敷設されている各
部の温度分布を計測する。
The temperature distribution detecting device shown in this figure includes an optical fiber 1 and a signal processing device 2, and generates a very narrow light pulse by the signal processing device 2 during initialization and periodically. Then, the light is made incident on the optical fiber 1, and the Fresnel reflected light from the end face of the optical fiber 1 is taken in. The deconvolution is performed based on the Fresnel reflected light by an iterative method to obtain a signal corresponding to the Fresnel reflected light. The optimum number of repetitions at which the waveform becomes linear is determined. When performing a normal measurement, the signal processing device 2
A very narrow optical pulse is generated, and the pulse is incident on the optical fiber 1. The backscattered light is taken in, and the optimum number of repetitions is determined based on the anti-Stokes Raman scattered light and the Stokes Raman scattered light. The temperature distribution of each part where the optical fiber 1 is laid is measured by performing deconvolution based on.

【0044】光ファイバ1は、温度の測定対象となる各
部分を通るように敷設されており、入射端から光パルス
が入射されたとき、これを取り込むとともに、各部の温
度に応じた後方散乱光および終端面でのフレネル反射光
を前記入射端に戻す。
The optical fiber 1 is laid so as to pass through each portion to be measured for temperature. When an optical pulse is incident from the incident end, the optical fiber 1 takes in the optical pulse, and the backscattered light corresponding to the temperature of each portion. And returning the Fresnel reflected light at the end face to the incident end.

【0045】また、信号処理装置2は、タイミング制御
回路4と、レーザ駆動回路5と、レーザ光源6と、カプ
ラ7と、光学フィルタ8と、光スイッチ9と、バイアス
制御回路10と、フォトダイオード11と、光電流増幅
回路12と、デジタルアベレージャ回路13と、制御記
憶演算回路14とを備えており、初期化時および定期的
に、非常に幅の狭い光パルスを生成してこれを光ファイ
バ1に入射させ、この光ファイバ1の終端面からのフレ
ネル反射光を取り込むとともに、このフレネル反射光に
基づき反復法によってデコンボルーションを行なってこ
のフレネル反射光に対応する信号の波形が線状となる最
適な反復回数を求める。そして、通常の測定を行なうと
き、非常に幅の狭い光パルスを生成してこれを光ファイ
バ1に入射させるとともに、この後方散乱光を取込ん
で、そのアンチストークスラマン散乱光と、ストークス
ラマン散乱光とに基づき前記最適な反復回数に基づいて
デコンボルーションを行なって光ファイバ1が敷設され
ている各部の温度分布を計測する。
The signal processing device 2 includes a timing control circuit 4, a laser driving circuit 5, a laser light source 6, a coupler 7, an optical filter 8, an optical switch 9, a bias control circuit 10, a photodiode 11, a photocurrent amplifier circuit 12, a digital averager circuit 13, and a control storage operation circuit 14. At initialization and at regular intervals, a very narrow optical pulse is generated to generate an optical pulse. The incident light is incident on the fiber 1, and the Fresnel reflected light from the end face of the optical fiber 1 is taken in. The deconvolution is performed based on the Fresnel reflected light by an iterative method, and the waveform of the signal corresponding to the Fresnel reflected light is linear. Find the optimal number of iterations. When a normal measurement is performed, a very narrow light pulse is generated and incident on the optical fiber 1, and the backscattered light is taken in. The anti-Stokes Raman scattering light and the Stokes Raman scattering Based on the light, deconvolution is performed based on the optimum number of repetitions, and the temperature distribution of each part where the optical fiber 1 is laid is measured.

【0046】タイミング制御回路4は、前記制御記憶演
算回路14と信号の授受を行ないながら、発光タイミン
グ信号や同期加算タイミング信号などを生成する部分で
あり、生成動作によって得られた発光タイミング信号を
レーザ駆動回路5に供給してこれを制御するとともに、
同期加算タイミング信号をデジタルアベレージャ回路1
3に供給してこれを制御する。
The timing control circuit 4 generates a light emission timing signal and a synchronous addition timing signal while transmitting and receiving signals to and from the control storage arithmetic circuit 14. The timing control circuit 4 converts the light emission timing signal obtained by the generation operation into a laser. This is supplied to the drive circuit 5 to control it,
Synchronous addition timing signal to digital averager circuit 1
3 to control this.

【0047】レーザ駆動回路5は、前記タイミング制御
回路4から出力される発光タイミング信号に基づいて非
常に幅の狭い駆動信号を生成する部分であり、この駆動
信号を前記レーザ光源6に供給してこれを発光させる。
The laser drive circuit 5 is a portion for generating a drive signal having a very narrow width based on the light emission timing signal output from the timing control circuit 4. The drive signal is supplied to the laser light source 6. This emits light.

【0048】レーザ光源6は、前記レーザ駆動回路5か
ら駆動信号が出力されている間、発光して予め設定され
ている波長の光パルスを生成する部分であり、生成した
光パルスをカプラ7に供給する。
The laser light source 6 emits light while the driving signal is being output from the laser driving circuit 5 and generates an optical pulse having a predetermined wavelength. Supply.

【0049】カプラ7は、前記レーザ光源6から出力さ
れる光パルスを前記光ファイバ1の入射端に導いたり、
この入射端から出射される後方散乱光および終端面での
フレネル反射光を前記光学フィルタ8に導いたりする部
分であり、前記レーザ光源6から出力される光パルスを
取り込んでこれを前記光ファイバ1の入射端に導いて光
ファイバ1中に入射させるとともに、前記光パルスが前
記光ファイバ1中を伝搬する際に発生し、前記入射端か
ら出射される後方散乱光や終端面で生じるフレネル反射
光を取り込んでこれを前記光学フィルタ8に導く。
The coupler 7 guides an optical pulse output from the laser light source 6 to the incident end of the optical fiber 1,
This portion guides the backscattered light emitted from the incident end and the Fresnel reflected light at the terminal surface to the optical filter 8. The optical pulse output from the laser light source 6 is taken in and the optical pulse is transmitted to the optical fiber 1. And incident on the optical fiber 1 while the optical pulse propagates through the optical fiber 1, and the backscattered light emitted from the incident end and the Fresnel reflected light generated at the terminal surface And guide it to the optical filter 8.

【0050】光学フィルタ8は、前記カプラ7から出射
される後方散乱光やフレネル反射光を取り込むととも
に、この後方散乱光やフレネル反射光を波長弁別してア
ンチストークスラマン散乱光と、ストークスラマン散乱
光とに分離する部分であり、この分離動作によって得ら
れたアンチストークスラマン散乱光と、ストークスラマ
ン散乱光とを光スイッチ9に導く。
The optical filter 8 takes in the backscattered light and Fresnel reflected light emitted from the coupler 7 and discriminates the wavelength of the backscattered light and Fresnel reflected light to generate anti-Stokes Raman scattered light and Stokes Raman scattered light. The anti-Stokes Raman scattered light and the Stokes Raman scattered light obtained by this separation operation are guided to the optical switch 9.

【0051】光スイッチ9は、前記制御記憶演算回路1
4からの選択信号に基づいて前記光学フィルタ8から出
射される前記アンチストークスラマン散乱光またはスト
ークスラマン散乱光のいずれか一方を選択する部分であ
り、この選択処理によって得られたアンチストークスラ
マン散乱光またはストークスラマン散乱光のいずれか一
方をフォトダイオード11に導く。
The optical switch 9 is connected to the control storage arithmetic circuit 1
4 is a portion for selecting either the anti-Stokes Raman scattered light or the Stokes Raman scattered light emitted from the optical filter 8 based on the selection signal from the optical filter 8, and the anti-Stokes Raman scattered light obtained by the selection processing. Alternatively, one of the Stokes Raman scattered light is guided to the photodiode 11.

【0052】また、バイアス制御回路10は、前記制御
記憶演算回路14から出力されるバイアス指令信号に基
づいてバイアス電圧を生成する部分であり、前記バイア
ス指令信号が高増幅率指令であるときには、予め設定さ
れている高い電圧を生成してこれを前記フォトダイオー
ド11に供給し、また前記バイアス指令信号が低増幅率
指令であるときには、予め設定されている低い電圧を生
成してこれを前記フォトダイオード11に供給する。
The bias control circuit 10 generates a bias voltage based on the bias command signal output from the control storage arithmetic circuit 14. When the bias command signal is a high amplification rate command, A preset high voltage is generated and supplied to the photodiode 11, and when the bias command signal is a low amplification rate command, a low voltage set in advance is generated and is supplied to the photodiode 11. 11

【0053】フォトダイオード11は、シリコン等によ
って構成されるアバランシェダイオードであり、前記バ
イアス制御回路10から出力されるバイアス電圧が高い
電圧であるとき、前記光スイッチ9から出射されるアン
チストークスラマン散乱光またはストークスラマン散乱
光を高い増幅率で光電流に変換してこれを光電流増幅回
路12に供給し、また前記バイアス制御回路10から出
力されるバイアス電圧が低い電圧であるとき、前記光ス
イッチ9から出射されるアンチストークスラマン散乱光
またはストークスラマン散乱光を低い増幅率(例えば、
増幅率“1”)で光電流に変換してこれを光電流増幅回
路12に供給する。
The photodiode 11 is an avalanche diode made of silicon or the like. When the bias voltage output from the bias control circuit 10 is high, the anti-Stokes Raman scattering light emitted from the optical switch 9 is used. Alternatively, the Stokes Raman scattered light is converted into a photocurrent at a high amplification rate and supplied to the photocurrent amplifier circuit 12. When the bias voltage output from the bias control circuit 10 is a low voltage, the optical switch 9 Anti-Stokes Raman scattered light or Stokes Raman scattered light emitted from a low amplification factor (for example,
The photocurrent is converted at the amplification factor “1”) and supplied to the photocurrent amplifier circuit 12.

【0054】光電流増幅回路12は、前記制御記憶演算
回路14から出力される増幅率指定信号に応じた増幅率
で前記フォトダイオード11の光電流を増幅して散乱信
号やフレネル反射信号を生成する部分であり、増幅動作
によって得られた散乱信号やフレネル反射信号をデジタ
ルアベレージャ回路13に供給する。
The photocurrent amplifier circuit 12 amplifies the photocurrent of the photodiode 11 at an amplification factor according to the amplification factor designation signal output from the control storage operation circuit 14 to generate a scattered signal and a Fresnel reflection signal. And a scattered signal and a Fresnel reflected signal obtained by the amplification operation are supplied to the digital averager circuit 13.

【0055】デジタルアベレージャ回路13は、前記タ
イミング制御回路4から出力される同期加算タイミング
信号に基づいて前記光電流増幅回路12から出力される
散乱信号やフレネル反射信号を取り込んでこれを同期加
算してアンチストークスラマン散乱信号のアベレージ出
力またはストークスラマン散乱信号のアベレージ出力、
フレネル反射信号のアベレージ出力を生成してこれを制
御記憶演算回路14に供給する。
The digital averager circuit 13 takes in the scattered signal and the Fresnel reflection signal output from the photocurrent amplifier circuit 12 based on the synchronous addition timing signal output from the timing control circuit 4 and synchronously adds them. Average output of anti-Stokes Raman scattering signal or average output of Stokes Raman scattering signal,
An average output of the Fresnel reflection signal is generated and supplied to the control storage operation circuit 14.

【0056】制御記憶演算回路14は、この温度分布検
出装置全体の制御を行なう部分であり、前記タイミング
制御回路4やバイアス制御回路10を制御して発光タイ
ミングや同期加算タイミングを調整する処理や前記デジ
タルアベレージャ回路13から出力されるアンチストー
クスラマン散乱信号のアベレージ出力、またはストーク
スラマン散乱信号のアベレージ出力、フレネル反射光の
アベレージ出力を取り込んでこれを記憶する処理、初期
化時や予め設定されている周期毎に、記憶しているフレ
ネル反射光のアベレージ出力に対し、上述したデコンボ
ルーションを行なって最適な反復回数(最適反復回数)
を求める処理、測定指令が設定されているとき、記憶し
ているアンチストークスラマン散乱信号のアベレージ出
力およびストークスラマン散乱信号のアベレージ出力に
対し、前記最適反復回数だけ上述したデコンボルーショ
ンを行なって温度信号を生成する処理、この処理によっ
て得られた温度信号を温度分布出力インタフェースを介
して他の装置に出力する処理などを行なう。
The control storage operation circuit 14 controls the entire temperature distribution detecting device, and controls the timing control circuit 4 and the bias control circuit 10 to adjust the light emission timing and the synchronous addition timing. An average output of the anti-Stokes Raman scattering signal output from the digital averager circuit 13, an average output of the Stokes Raman scattering signal, an average output of the Fresnel reflected light, and a process of storing the average output. For each period, the above-mentioned deconvolution is performed on the stored average output of the Fresnel reflected light, and the optimum number of repetitions (optimal number of repetitions) is obtained.
When the measurement command is set, the average output of the stored anti-Stokes Raman scattering signal and the average output of the Stokes Raman scattering signal are subjected to the above-described deconvolution for the optimal number of repetitions, and the temperature is calculated. A process of generating a signal, a process of outputting a temperature signal obtained by this process to another device via a temperature distribution output interface, and the like are performed.

【0057】次に、図1に示すブロック図を参照しなが
ら、この実施例の最適反復回数検出動作、温度測定動作
を順次、説明する。
Next, the operation of detecting the optimum number of repetitions and the operation of measuring the temperature of this embodiment will be sequentially described with reference to the block diagram shown in FIG.

【0058】《最適反復回数検出動作》まず、この温度
分布検出装置の電源が投入されて回路の初期化を行なう
ときや予め設定されている較正タイミングになる毎に、
信号処理装置の制御記憶演算回路14はバイアス制御回
路10を制御して低い電圧を生成させ、フォトダイオー
ド11の増幅率を低い値(例えば、増幅率“1”)にさ
せた後、タイミング制御回路4を制御してレーザ駆動回
路5を駆動させ、このレーザ駆動回路5に接続されたレ
ーザ光源6に非常に幅の狭い光パルスを所定のタイミン
グで連続的に生成させて、これをカプラ7に供給させ
る。
<< Optimal Number of Repetitions Detecting Operation >> First, when the temperature distribution detecting device is powered on to initialize the circuit or every time a preset calibration timing is reached,
The control storage operation circuit 14 of the signal processing device controls the bias control circuit 10 to generate a low voltage and set the amplification factor of the photodiode 11 to a low value (for example, the amplification factor “1”), and then the timing control circuit 4 to drive a laser drive circuit 5 to cause a laser light source 6 connected to the laser drive circuit 5 to continuously generate a very narrow light pulse at a predetermined timing. Let it be supplied.

【0059】そして、このカプラ7によって前記光パル
スが光ファイバ1に導かれ、この光ファイバ1からフレ
ネル反射光が出射される毎に、これが光学フィルタ8お
よび光スイッチ9を介してフォトダイオード11に導か
れて光電流に変換されるとともに、光電流増幅回路12
によって増幅されてフレネル反射信号にされた後、デジ
タルアベレージャ回路13によって指定された回数だけ
同期加算されてフレネル反射信号のアベレージ出力が生
成され、これが制御記憶演算回路14に供給されて記憶
される。
The optical pulse is guided to the optical fiber 1 by the coupler 7, and every time the Fresnel reflected light is emitted from the optical fiber 1, this light is transmitted to the photodiode 11 via the optical filter 8 and the optical switch 9. While being guided and converted into a photocurrent, the photocurrent amplifier circuit 12
The signal is amplified by the digital averager circuit 13 to be a Fresnel reflection signal, and is synchronously added a specified number of times by the digital averager circuit 13 to generate an average output of the Fresnel reflection signal, which is supplied to the control storage arithmetic circuit 14 and stored. .

【0060】そして、この処理が終了すれば、制御記憶
演算回路14は記憶しているフレネル反射信号のアベレ
ージ出力に対し、デコンボルーション処理を繰り返しな
がら、波形が線状になる回数を最適反復回数として記憶
する。
When this processing is completed, the control storage arithmetic circuit 14 repeats the deconvolution processing on the stored average output of the Fresnel reflection signal, and determines the number of times the waveform becomes linear to the optimum number of repetitions. To be stored.

【0061】《温度測定動作》そして、上述した最適反
復回数検出動作が終了すると、制御記憶演算回路14は
光スイッチ9を制御してアンチストークスラマン散乱光
を選択させた後、タイミング制御回路4を制御してレー
ザ駆動回路5を駆動させ、このレーザ駆動回路5に接続
されたレーザ光源6に非常に幅の狭い光パルスを所定の
タイミングで連続的に生成させて、これをカプラ7に供
給させる。
<< Temperature Measuring Operation >> When the above-described operation for detecting the optimum number of repetitions is completed, the control storage arithmetic circuit 14 controls the optical switch 9 to select the anti-Stokes Raman scattered light, and then switches the timing control circuit 4 The laser drive circuit 5 is controlled to drive the laser drive circuit 5, and the laser light source 6 connected to the laser drive circuit 5 continuously generates a very narrow light pulse at a predetermined timing and supplies the light pulse to the coupler 7. .

【0062】そして、このカプラ7によって前記光パル
スが光ファイバ1に導かれて、この光ファイバ1から後
方散乱光が出射される毎に、これが光学フィルタ8によ
って波長弁別されてアンチストークスラマン散乱光と、
ストークスラマン散乱光とに分離される。
The optical pulse is guided to the optical fiber 1 by the coupler 7, and each time the backscattered light is emitted from the optical fiber 1, the backscattered light is discriminated in wavelength by the optical filter 8 and the anti-Stokes Raman scattered light is separated. When,
It is separated into Stokes Raman scattered light.

【0063】次いで、光スイッチ9によって前記光学フ
ィルタ8の波長弁別処理で得られたアンチストークスラ
マン散乱光が選択されてこれがフォトダイオード11に
導かれて光電流に変換されるとともに、光電流増幅回路
12によって増幅されてアンチストークスラマン散乱信
号にされた後、デジタルアベレージャ回路13によって
指定された回数だけ同期加算されてアンチストークスラ
マン散乱信号のアベレージ出力が生成され、これが制御
記憶演算回路14に供給されて記憶される。
Next, the anti-Stokes Raman scattered light obtained by the wavelength discrimination processing of the optical filter 8 is selected by the optical switch 9 and guided to the photodiode 11 to be converted into a photocurrent. The amplified signal is amplified by an anti-Stokes Raman scattering signal 12 and is synchronously added by a digital averager circuit 13 a specified number of times to generate an average output of the anti-Stokes Raman scattering signal, which is supplied to a control storage arithmetic circuit 14. Is stored.

【0064】この後、この処理が終了すると、制御記憶
演算回路14は光スイッチ9を制御してストークスラマ
ン散乱光を選択させた後、タイミング制御回路4を制御
してレーザ駆動回路5を駆動させ、このレーザ駆動回路
5に接続されたレーザ光源6に非常に幅の狭い光パルス
を所定のタイミングで連続的に生成させて、これをカプ
ラ7に供給させる。
Thereafter, when this processing is completed, the control storage arithmetic circuit 14 controls the optical switch 9 to select Stokes Raman scattered light, and then controls the timing control circuit 4 to drive the laser drive circuit 5. The laser light source 6 connected to the laser drive circuit 5 continuously generates a very narrow light pulse at a predetermined timing and supplies the light pulse to the coupler 7.

【0065】そして、このカプラ7によって前記光パル
スが光ファイバ1に導かれて、この光ファイバ1から後
方散乱光が出射される毎に、これが光学フィルタ8によ
って波長弁別されてアンチストークスラマン散乱光と、
ストークスラマン散乱光とに分離される。
The optical pulse is guided to the optical fiber 1 by the coupler 7, and each time the backscattered light is emitted from the optical fiber 1, the backscattered light is discriminated in wavelength by the optical filter 8, and the anti-Stokes Raman scattered light is separated. When,
It is separated into Stokes Raman scattered light.

【0066】次いで、光スイッチ9によって前記光学フ
ィルタ8の波長弁別処理で得られたストークスラマン散
乱光が選択されてこれがフォトダイオード11に導かれ
て光電流に変換されるとともに、光電流増幅回路12に
よって増幅されてストークスラマン散乱信号にされた
後、デジタルアベレージャ回路13によって指定された
回数だけ同期加算されてストークスラマン散乱信号のア
ベレージ出力が生成され、これが制御記憶演算回路14
に供給されて記憶される。
Next, the Stokes Raman scattered light obtained by the wavelength discrimination processing of the optical filter 8 is selected by the optical switch 9 and guided to the photodiode 11 to be converted into a photocurrent. The signal is amplified and converted into a Stokes Raman scattering signal, and is synchronously added a number of times designated by the digital averager circuit 13 to generate an average output of the Stokes Raman scattering signal.
And stored.

【0067】そして、この処理が終了すると、制御記憶
演算回路14は記憶しているアンチストークスラマン散
乱信号のアベレージ出力と、ストークスラマン散乱信号
のアベレージ出力とに対し、上述した最適反復回数検出
動作によって得られた最適反復回数だけデコンボルーシ
ョン処理を繰り返しながら、両者の比を計算して最終的
に得られた温度信号を温度分布出力インタフェースを介
して他の装置に出力する。
When this processing is completed, the control storage arithmetic circuit 14 performs the above-described optimum repetition number detection operation on the stored average output of the anti-Stokes Raman scattering signal and the average output of the Stokes Raman scattering signal. While repeating the deconvolution process by the obtained optimum number of repetitions, the ratio between the two is calculated, and the finally obtained temperature signal is output to another device via the temperature distribution output interface.

【0068】このようにこの実施例においては、初期化
時および定期的に、信号処理装置2によって非常に幅の
狭い光パルスを生成してこれを光ファイバ1に入射さ
せ、この光ファイバ1の終端面からのフレネル反射光を
取り込むとともに、このフレネル反射光に基づき反復法
によってデコンボルーションを行なってこのフレネル反
射光に対応する信号の波形が線状となる最適な反復回数
を求め、通常の測定を行なうとき、信号処理装置2によ
って非常に幅の狭い光パルスを生成してこれを光ファイ
バ1に入射させ、この後方散乱光を取込んで、そのアン
チストークスラマン散乱光と、ストークスラマン散乱光
とに基づき前記最適な反復回数に基づいてデコンボルー
ションを行なって光ファイバ1が敷設されている各部の
温度分布を計測するようにしたので、逆変換過程におい
て雑音の影響で解が発散してしまうのを防止し、これに
よって光パルスの幅を狭くすることなく総合精度を最も
改善させることができる最適な近似を与えることができ
るとともに、光パルス幅を狭くした場合に生じるS/N
比の低下やこれに伴う計測時間の増大、半導体レーザの
駆動回路の困難さを生じることなく、位置分解能を向上
させてホットスポットなどの異常加熱部などの検出を容
易にすることができる。
As described above, in this embodiment, at the time of initialization and at regular intervals, a very narrow optical pulse is generated by the signal processing device 2 and is incident on the optical fiber 1. While taking in the Fresnel reflected light from the terminal surface, deconvolution is performed based on the Fresnel reflected light by an iterative method, and the optimum number of repetitions at which the waveform of the signal corresponding to the Fresnel reflected light becomes linear is obtained. When performing the measurement, a very narrow light pulse is generated by the signal processing device 2 and made incident on the optical fiber 1, and the backscattered light is taken in, and the anti-Stokes Raman scattered light and the Stokes Raman scattered light are taken. The deconvolution is performed based on the optimal number of repetitions based on the light and the temperature distribution of each part where the optical fiber 1 is laid is measured. Therefore, it is necessary to prevent the solution from diverging due to the influence of noise in the inverse conversion process, and thereby to give an optimal approximation that can improve the overall accuracy most without reducing the width of the light pulse. S / N generated when the optical pulse width is reduced
It is possible to improve the positional resolution and to easily detect an abnormally heated portion such as a hot spot without lowering the ratio, increasing the measuring time and making the driving circuit of the semiconductor laser difficult.

【0069】また、上述した実施例においては、フォト
ダイオード11に印加するバイアス電圧の値を切り替え
て増幅率を小さくすることにより、フレネル反射光を光
電流に変換するようにしているが、光ファイバ1の入射
面や終端面にコーティング等の処理を施してフレネル反
射光の大きさをアンチストークスラマン散乱光やストー
クスラマン散乱光の大きさと同程度にしてフォトダイオ
ード11のバイアス電圧を変更することなく、このフォ
トダイオード11の飽和を防止しながら、フレネル反射
光を光電流に変換するようにしても、上述した実施例と
同様な効果を得ることができる。
In the above embodiment, the Fresnel reflected light is converted into a photocurrent by switching the value of the bias voltage applied to the photodiode 11 to reduce the amplification factor. No. 1 is subjected to a treatment such as coating on the incident surface or the terminal surface so that the magnitude of the Fresnel reflected light is substantially equal to the magnitude of the anti-Stokes Raman scattered light or the Stokes Raman scattered light without changing the bias voltage of the photodiode 11. Also, even if the Fresnel reflected light is converted into a photocurrent while preventing the photodiode 11 from being saturated, the same effect as in the above-described embodiment can be obtained.

【0070】[0070]

【発明の効果】以上説明したように本発明によれば、逆
変換過程において雑音の影響による解の発散を防止し、
これによって光パルスの幅を狭くすることなく総合精度
を最も改善させることができる最適な近似を与えること
ができるとともに、光パルス幅を狭くした場合に生じる
S/N比の低下やこれに伴う計測時間の増大、半導体レ
ーザの駆動回路の困難さを生じることなく、位置分解能
を向上させてホットスポットなどの異常加熱部などの検
出を容易にすることができる。
As described above, according to the present invention, it is possible to prevent the divergence of the solution due to the influence of noise in the inverse transformation process.
As a result, an optimum approximation that can improve the overall accuracy without reducing the width of the light pulse can be provided, and the S / N ratio decreases when the light pulse width is reduced and the measurement accompanying the reduction can be performed. The position resolution can be improved and the detection of an abnormally heated portion such as a hot spot can be facilitated without increasing the time and making the driving circuit of the semiconductor laser difficult.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明による温度分布検出装置の一実施例を示
す構成図である。
FIG. 1 is a configuration diagram showing one embodiment of a temperature distribution detecting device according to the present invention.

【図2】本発明による温度分布検出装置の基本原理で使
用されるデコンボルーション処理の手順例を示す模式図
である。
FIG. 2 is a schematic diagram showing an example of a procedure of a deconvolution process used in the basic principle of the temperature distribution detecting device according to the present invention.

【図3】本発明による温度分布検出装置の基本原理を説
明するための装置構成例を示す構成図である。
FIG. 3 is a configuration diagram showing an example of a device configuration for explaining a basic principle of a temperature distribution detection device according to the present invention.

【図4】図3に示す温度分布検出装置によって得られる
温度信号例を示す波形図である。
FIG. 4 is a waveform chart showing an example of a temperature signal obtained by the temperature distribution detecting device shown in FIG.

【図5】図4に示す温度信号に対するデコンボルーショ
ン処理例を示す波形図である。
5 is a waveform chart showing an example of deconvolution processing for the temperature signal shown in FIG.

【図6】本発明による温度分布検出装置の基本原理で使
用されるデコンボルーション処理の系統誤差と、ランダ
ム誤差との関係例を示す表図である。
FIG. 6 is a table showing an example of a relationship between a systematic error of a deconvolution process used in the basic principle of the temperature distribution detecting device according to the present invention and a random error.

【図7】本発明による温度分布検出装置の基本原理で使
用されるデコンボルーション処理の反復回数と、フレネ
ル反射光の信号との関係を説明するための波形図であ
る。
FIG. 7 is a waveform diagram for explaining the relationship between the number of repetitions of the deconvolution process used in the basic principle of the temperature distribution detecting device according to the present invention and the Fresnel reflected light signal.

【図8】光ファイバの後方散乱光を利用したセンサの温
度計測原理を示す模式図である。
FIG. 8 is a schematic diagram illustrating a principle of measuring a temperature of a sensor using backscattered light of an optical fiber.

【符号の説明】[Explanation of symbols]

1 光ファイバ 2 信号処理装置 4 タイミング制御回路 5 レーザ駆動回路 6 レーザ光源(光源) 7 カプラ 8 光学フィルタ(光学系) 9 光スイッチ(光学系) 10 バイアス制御回路 11 フォトダイオード(光電変換部) 12 光電流増幅回路(光電変換部) 13 デジタルアベレージャ回路(デジタルアベレージ
ャ部) 14 制御記憶演算回路(演算部)
Reference Signs List 1 optical fiber 2 signal processing device 4 timing control circuit 5 laser drive circuit 6 laser light source (light source) 7 coupler 8 optical filter (optical system) 9 optical switch (optical system) 10 bias control circuit 11 photodiode (photoelectric conversion unit) 12 Photocurrent amplifier circuit (photoelectric conversion unit) 13 Digital averager circuit (digital averager unit) 14 Control storage operation circuit (operation unit)

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G01K 11/12 G01K 3/00 ──────────────────────────────────────────────────続 き Continued on the front page (58) Field surveyed (Int.Cl. 7 , DB name) G01K 11/12 G01K 3/00

Claims (1)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 光パルスを発生する光源と、前記光パル
スを光ファイバに入射し、その後方散乱光から2種類の
ラマン散乱光を選択する光学系と、前記ラマン散乱光を
電気信号に変換増幅する光電変換部と、この光電変換部
の出力信号をデジタル信号に変換してアベレージングを
行なうデジタルアベレージャ部と、このデジタルアベレ
ージャ部の処理結果に基づいて前記光ファイバに沿った
温度分布を計算する演算部とを有する温度分布検出装置
において、 光ファイバの端面からのフレネル反射光の波形データに
対し、デコンボルーション処理を行なって前記波形デー
タの波形が所定の形状になるまでの反復回数を求め、温
度測定を行なうとき、2種類のラマン散乱光の波形デー
タに対し、前記反復回数だけデコンボルーション処理を
繰り返して光ファイバに沿った温度分布を計算する、 ことを特徴とする温度分布検出装置。
1. A light source for generating an optical pulse, an optical system for injecting the optical pulse into an optical fiber and selecting two types of Raman scattered light from back scattered light, and converting the Raman scattered light into an electric signal A photoelectric conversion unit that amplifies, a digital averager unit that converts an output signal of the photoelectric conversion unit into a digital signal and performs averaging, and a temperature distribution along the optical fiber based on a processing result of the digital averager unit. A temperature distribution detecting device having a calculation unit for calculating the waveform data of the Fresnel reflected light from the end face of the optical fiber. When the number of times is determined and the temperature is measured, the deconvolution process is repeated for the waveform data of two types of Raman scattered light by the number of repetitions. Calculating the temperature distribution along the optical fiber Te, the temperature distribution detection device according to claim.
JP5000646A 1993-01-06 1993-01-06 Temperature distribution detector Expired - Fee Related JP3065833B2 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP5000646A JP3065833B2 (en) 1993-01-06 1993-01-06 Temperature distribution detector
KR94000146A KR0133488B1 (en) 1993-01-06 1994-01-06 Temperature distribution detector using optical fiber
US08/528,594 US5639162A (en) 1993-01-06 1995-09-15 Temperature distribution measuring apparatus using an optical fiber
US08/813,980 US5825804A (en) 1993-01-06 1997-03-10 Temperature distribution measuring apparatus using an optical fiber

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP5000646A JP3065833B2 (en) 1993-01-06 1993-01-06 Temperature distribution detector

Publications (2)

Publication Number Publication Date
JPH06201488A JPH06201488A (en) 1994-07-19
JP3065833B2 true JP3065833B2 (en) 2000-07-17

Family

ID=11479476

Family Applications (1)

Application Number Title Priority Date Filing Date
JP5000646A Expired - Fee Related JP3065833B2 (en) 1993-01-06 1993-01-06 Temperature distribution detector

Country Status (1)

Country Link
JP (1) JP3065833B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4443988B2 (en) * 2004-04-16 2010-03-31 三菱重工業株式会社 Status detector for radioactive material containers
CN102230831B (en) * 2011-03-30 2013-04-10 中国科学院苏州纳米技术与纳米仿生研究所 Method for measuring substrate temperature and device layer temperature in LED (light-emitting diode) simultaneously
WO2014064770A1 (en) * 2012-10-23 2014-05-01 富士通株式会社 Abnormality detecting system and abnormality detecting method
CN114486179B (en) * 2022-02-21 2023-06-13 电子科技大学中山学院 High-precision optical fiber quality detection method and system for deconvolution filtering

Also Published As

Publication number Publication date
JPH06201488A (en) 1994-07-19

Similar Documents

Publication Publication Date Title
US5825804A (en) Temperature distribution measuring apparatus using an optical fiber
US9322721B2 (en) Optic fiber distributed temperature sensor system with self-correction function and temperature measuring method using thereof
US5054935A (en) Temperature-measuring method and distributed optical fiber temperature sensor
EP0502283B1 (en) Temperature distribution analyzer using optical fiber
US7873273B2 (en) Apparatus for measuring the characteristics of an optical fiber
US20210116311A1 (en) Self-calibration detection device and temperature demodulation method oriented to fiber raman temperature sensing system
CN104880209B (en) Distributed optical sensing with two-step assessment
EP3141880B1 (en) Fiber measurement with pulse shaping
JP2009156718A (en) Optical pulse test equipment
JP3065833B2 (en) Temperature distribution detector
JP3065832B2 (en) Temperature distribution detector
CN212320747U (en) Linear frequency modulation&#39;s distributed optical fiber sensing device
JP2584478B2 (en) Method for processing received signal of optical fiber backscattered light
CN109444895B (en) Vibration information positioning method for eliminating interference fading of distributed vibration sensor
JP5515199B2 (en) Optical pulse test apparatus and adjustment method thereof
JP2747565B2 (en) Method and apparatus for measuring curvature distribution of optical fiber
JPS58165031A (en) Method and device for measuring light transmission characteristic of optical fiber
JP2009281813A (en) Optical fiber measuring device
US20040227931A1 (en) Optic communication or transmission media sensing
CN116907640A (en) Method for determining ROTDR optimal optical pulse power based on fibonacci number
JP3167202B2 (en) Temperature measurement method and device
CN112665753B (en) Encoding device and method for distributed optical fiber temperature sensing under optical amplification
JP2885979B2 (en) Temperature distribution detector
JP4803083B2 (en) Optical fiber measuring apparatus and sampling method
EP0293441B1 (en) Improvements relating to optical sensing systems

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees