特許 7511745 飽和によるＤＡＳの位相ジャンプの回避

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-06-27

(45)【発行日】2024-07-05

(54)【発明の名称】飽和によるＤＡＳの位相ジャンプの回避

(51)【国際特許分類】

   G01D 5/353 20060101AFI20240628BHJP

【ＦＩ】

G01D5/353 B

【請求項の数】 5

(21)【出願番号】P 2023512675

(86)(22)【出願日】2021-10-27

(65)【公表番号】

(43)【公表日】2023-10-05

(86)【国際出願番号】 US2021056797

(87)【国際公開番号】W WO2022093920

(87)【国際公開日】2022-05-05

【審査請求日】2023-03-17

(31)【優先権主張番号】63/105,984

(32)【優先日】2020-10-27

(33)【優先権主張国・地域又は機関】US

(31)【優先権主張番号】17/511,381

(32)【優先日】2021-10-26

(33)【優先権主張国・地域又は機関】US

(73)【特許権者】

【識別番号】504080663

【氏名又は名称】エヌイーシー ラボラトリーズ アメリカ インク

【氏名又は名称原語表記】ＮＥＣ Ｌａｂｏｒａｔｏｒｉｅｓ Ａｍｅｒｉｃａ， Ｉｎｃ．

(74)【代理人】

【識別番号】100123788

【弁理士】

【氏名又は名称】宮崎 昭夫

(74)【代理人】

【識別番号】100127454

【弁理士】

【氏名又は名称】緒方 雅昭

(72)【発明者】

【氏名】ホ、 ジュンチャン

【審査官】吉田 久

(56)【参考文献】

【文献】米国特許出願公開第２０１２／０１３７７８１（ＵＳ，Ａ１）

【文献】国際公開第２０２０／０７０２２９（ＷＯ，Ａ１）

【文献】特開２００８－１７５７４６（ＪＰ，Ａ）

【文献】特開昭６１－４１９３２（ＪＰ，Ａ）

【文献】特表２０１５－５２３５５１（ＪＰ，Ａ）

【文献】米国特許出願公開第２００８／００８８８４６（ＵＳ，Ａ１）

【文献】特表２０１１－５０９０５８（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０１Ｄ ５／２６－５／３８

Ｇ０１Ｂ １１／００－１１／３０

Ｇ０２Ｂ ６／００

Ｇ０１Ｍ １１／００－１１／０２

Ｇ０１Ｊ ９／００－９／０４

Ｇ０１Ｖ １／００－１／５２

Ｇ０１Ｈ ９／００

(57)【特許請求の範囲】

【請求項1】

光ファイバと、
光パルスを生成し、それら光パルスを前記光ファイバに取り込み、前記光ファイバから一連のレイリー反射信号を受信する光インタロゲータユニットと、
前記レイリー反射信号から情報を抽出するように構成されたコヒーレント受信ユニットと、を有するＤＯＦＳ／ＤＡＳシステムのための分散型光ファイバセンシング（ＤＯＦＳ）／分散型音響センシング（ＤＡＳ）方法であって、
前記ＤＯＦＳ／ＤＡＳシステムを動作させて、前記光ファイバに沿った複数の位置について、ビート積のシーケンスを取得（サンプリング)することと、
前記シーケンス内の２つのサンプルごとの位相差がオーバーフロー状態になるかどうかを決定することと、
決定されたオーバーフローの位相差を、固定値だけ増減して正常値になるように調整することと、
前記オーバーフロー状態と該オーバーフローの方向とを示すフラグを出力することと、
前記調整された位相差およびフラグをハイパスフィルタに提供することと、を含み、
前記ハイパスフィルタにおいて、バッファされたｘ（ｎ－ｍ）個のサンプル（ｍ＝０，１,２，...）をｘ’（ｎ－ｍ）＝ｘ（ｘ－ｍ）＋Ｎ・２πに置き換え、ここでＮは整数である、方法。

【請求項2】

前記固定値は、２πおよび２Ｎπからなる群から選択されるものであり、ここで、２Ｎπ≦（正常値のサポート範囲）であり、２（Ｎ＋１）π≧（正常値のサポート範囲）である、請求項１に記載の方法。

【請求項3】

前記ハイパスフィルタは、カスケード接続された２次ＩＩＲサブフィルタを使用して実装される無限インパルス応答（ＩＩＲ）を有するハイパスフィルタである、請求項２に記載の方法。

【請求項4】

第１の２次ＩＩＲフィルターは、

【数1】

で表される直接型であり、ｘ（ｎ）がオーバフローした場合、

【数2】

について、ｘ（ｎ＋ｍ）はｘ’（ｎ＋ｍ）＝ｘ（ｎ＋ｍ）－Δ（ｎ）で置き換えられ、

【数3】

となる、請求項３に記載の方法。

【請求項5】

前記ハイパスフィルタは、有限インパルス応答（ＦＩＲ）を有するハイパスフィルタであり、前記シーケンス内のｘ（ｎ）より前のサンプルは全て、位相アンラップと同じ値である、請求項２に記載の方法。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、一般に、分散型光ファイバセンシング（ＤＯＦＳ）システム、方法、および構造に関する。より詳細には、本開示は、コヒーレント分散型音響センシング（ＤＡＳ）における飽和による位相ジャンプの回避を提供するシステム、方法、および構造に関する。

【背景技術】

【0002】

光ファイバ内のコヒーレントレイリー後方散乱を用いた分散型音響センシング（ＤＡＳ）は、多くの安全性、セキュリティ、および完全性監視システムにおけるその有用な応用により、リアルタイムで複数の動的事象を監視するユビキタスな技術となっている。そのため、本技術の改善は、歓迎すべきものであろう。

【発明の概要】

【0003】

複素シーケンスの位相から情報を抽出し、アンラップ中にハイパスフィルタの出力にオーバーフローまたは「スパイク」を示すＤＡＳシステムにおいて、飽和による位相ジャンプを回避するためのシステム、方法、および構造を対象とする本開示の態様によれば、当技術分野における進歩がもたらされる。

【0004】

本開示の一態様によれば、従来技術とは対照的に、位相アンラップ中（出力信号がサポート範囲を超える場合）に、調整されるべき方向に応じて、Ｎが負または正の整数である位相にＮ・２πを加算することによって、範囲内に戻るように調整される。

【0005】

調整された位相は、フラグ信号と共に、ハイパスフィルタ（ＨＰＦ）に出力され、ＨＰＦにバッファされたｘ（ｎ－ｍ）個のサンプル（ｍ＝０，１,２，...）は、ｘ’（ｎ－ｍ）＝ｘ（ｘ－ｍ）＋Ｎ・２πに置き換えられ、出力のスパイクが除去される。

【図面の簡単な説明】

【0006】

本開示のより完全な理解は、添付図面を参照することによって実現され得る。

【0007】

【図1】本開示の態様による、例示的なコヒーレントＤＡＳ処理手順の概略フローチャート図である。

【0008】

【図2】本開示の態様による、アンラップ位相飽和調整およびハイパスフィルタ（ＨＰＦ）出力における対応するスパイクを示すプロットである。

【0009】

【図3(A)】本開示の態様による方法の例示的な処理ステップを示す一対のフローチャートである。

【図3(B)】本開示の態様による方法の例示的な処理ステップを示す一対のフローチャートである。

【0010】

【図4】従来技術で実行される位相アンラップを示す概略図である。

【0011】

【図5】本開示の態様による、オーバーフローが発生したときのアンラップ位相処理を示すプロットである。

【0012】

【図6】従来技術で実行される直接型ＩＩＲフィルタを示す概略図である。

【0013】

【図7】本発明の態様に係る第１の２次ＤＦ－Ｉ型ＩＩＲフィルタの実装例を示す概略図である。

【0014】

【図8】本開示の態様による改善されたＨＰＦ出力を示すプロットである。

【0015】

【図9】本開示の態様による、インタロゲータとコヒーレント受信機と解析システムとを含む例示的なＤＦＳシステムの概略ブロック図である。

【0016】

例示的な実施形態は、図面および詳細な説明によってより完全に説明される。しかしながら、本開示による実施形態は様々な形態で実施することができ、図面および詳細な説明に記載された特定のまたは例示的な実施形態に限定されない。

【発明を実施するための形態】

【0017】

以下は、単に本開示の原理を例示するものである。したがって、当業者は本明細書に明示的に記載または図示されていないが、本開示の原理を具現化し、その精神および範囲内に含まれる様々な構成を考案することができることが理解されよう。

【0018】

さらに、本明細書に記載されているすべての実施例および条件付き用語は、本開示の原理および技術を促進するために発明者によって寄与された概念を読者が理解するのを助けるための教育目的のためだけのものであることを意図しており、そのような具体的に列挙された実施例および条件に限定されないと解釈されるべきである。

【0019】

さらに、本開示の原理、態様、および実施形態を記載する本明細書のすべての記述、ならびにその具体例は、その構造的および機能的等価物の両方を包含することを意図している。さらに、そのような等価物は、現在知られている等価物と、将来開発される等価物、すなわち、構造に関係なく同じ機能を実行する開発された要素との両方を含むことが意図されている。

【0020】

したがって、たとえば、本明細書の任意のブロック図が、本開示の原理を実施する例示的な回路の概念図を表すことは、当業者には理解されるであろう。

【0021】

本明細書で特に明記しない限り、図面を構成する図は、縮尺通りに描かれていない。

【0022】

まず、作動位相抽出を用いたレイリー後方散乱信号のコヒーレント検出は、感度、縦方向応力に対する線形位相応答、および無制限の位相範囲の利点により、一部に起因して、分散音響センシング（ＤＡＳ）で一般的に使用されている周知の技術であることを最初に述べておく。

【0023】

当業者であれば、理解および認識されるように、複素信号（complex signal）に関する位相決定および計算は、通常、複数の処理ステップを要する。

【0024】

まず、アークタンジェント（ａｔａｎ（））を使用し、入力の符号を考慮することによって、（ｒｅａｌ（Ｓ（ｚ，ｎ）），ｉｍａｇ（Ｓ（ｚ，ｎ）））から（－π，π］範囲の角度φ（ｚ，ｎ）を計算する。ここで、

【数1】

は、位置ｚ、時間ｎにおける差動位相情報を含む複素信号である。

【0025】

次に、φ（ｚ，ｎ）をアンラップして、ｘ（ｚ，ｎ）－ｘ（ｚ，ｎ－１）（任意のｎについて）が（－π，π］の範囲に入るｘ（ｚ，ｎ）を取得する。

【0026】

ハイパスフィルタ（ＨＰＦ）は、通常、振動検出またはその他の目的で、アンラップされた位相に続いてＤＣを除去する。この手順を図１に概略的に示す。

【0027】

実際にはいくつかの理由で、アンラップされた位相の範囲は限られており、出力はその範囲に達する／それを超えて、オーバーフローを引き起こす可能性があることに留意されたい。その理由１つは、利用可能な処理リソースが原因で、特に固定小数点処理を使用するファームウェアで実装されている場合、出力信号のビット数が限られてしまうことである。

【0028】

別の理由は、アンラップされた位相自体が、ノイズによって引き起こされるランダムウォーク、または超大なファイバ応力レベル、またはレーザ位相ドリフトなどの他の要因のために、一方向に移動し、最終的に最大範囲に達する可能性があることである。

【0029】

アンラップされた位相が最大値（正または負のいずれか）に達すると、２つの処理の選択肢があり、すなわち、最大範囲に留まり（つまり、飽和したまま)、信号が反対方向に向かうまで待つか、または０にリセットするなどの正常な範囲に引き戻すか、またはアンラップされた値を使用して（－π、π］の範囲に変更するかのいずれかである。

【0030】

飽和状態のままでは、信号がフラットになるため検出機能が失われ、アンラップされた信号を正常な範囲に引き戻すと、２つのサンプル間に大きな段差（step）が生じ、ＨＰＦ出力に強いスパイクが発生する。このスパイクは、誤って振動として扱われる可能性がある。

【0031】

図２は、本開示の態様による、アンラップ位相飽和調整と、ハイパスフィルタ（ＨＰＦ）出力における対応するスパイクとを示すプロットである。

【0032】

図２に示される例では、ｘ軸はサンプルのシーケンスまたは時間ｎであり、ｙ軸は振幅であり、上側の線はアンラップ位相であり、下側の線はＤＣを除去するためにハイパスフィルタを通過する対応する出力である。上側の線の大きな段差（step）はアンラップされた位相が飽和するタイミングで起こり、アルゴリズムは、２πを減算することによって、それを正常な範囲に引き戻す。下側の線から分かるように、この大きな段差によって、ＨＰＦの出力に強いスパイクとそれに続く長いリップルが発生する。

【0033】

ここで示し記載するように、本開示の一態様は、図２に示すような場合に、ＨＰＦ出力からスパイクを除去する方法を説明する。説明を簡単にするために、以下では、特に明記しない限り、固定位置ｚを考慮する。

【0034】

位相アンラップ出力

【0035】

アンラップ位相が飽和するとき、本開示の態様による方法は、２πを減算または加算することによって出力を調整し、これは、飽和方向に応じて、（－π－２π，π－２π］または（－π＋２π，π＋２π］の範囲のｘ（ｎ）－ｘ（ｎ－１）を有する。後続の信号ｘ（ｎ＋ｉ＋１）－ｘ（ｎ＋ｉ）（ｉ＝０，１，２，．．．）は、ｘ（ｎ＋ｉ＋１）が飽和するまで、標準のアンラップ手順に従って範囲（－π，π］になる。アンラップ位相出力と共に、フラグを使用して、時間ｎにおける飽和状態と調整方向を示す。

【0036】

例えば、フラグは２ビット（「ｂ００」は「飽和なし」、「ｂ１０」は「飽和し、－２πによって調整」、「ｂ１１」は「飽和し、＋２πによって調整」）を使用することができる。

【0037】

ＨＰＦ処理

【0038】

ＨＰＦがＩＩＲ（無限インパルス応答）フィルタの場合は、カスケード接続された２次直接型Ｉ（ＤＦ－Ｉ）を使用して実装する。飽和のためにアンラップ位相が調整される場合、ｘ（ｎ－１）およびｘ（ｎ－２）を、ｘ（ｎ）に適用されたものと同じ量だけ増減する。

【0039】

例えば、新しい位相入力ｘ（ｎ）と飽和フラグ「ｂ０」（飽和を避けるためにｘ（ｎ）が２πだけ減少したことを意味する）を使用すると、ＤＦ－Ｉ構造のｘ（ｎ－１）とｘ（ｎ－２）はそれぞれｘ’（ｎ－１）＝ｘ（ｎ－１）－２πとｘ’（ｎ－２）＝ｘ（ｎ－２）－２πに置き換える必要があり、ｘ（ｎ），ｘ’（ｎ－１）が次のサイクルの動作で使用される。

【0040】

ＨＰＦが

【数2】

の形式のＦＩＲ（有限インパルス応答）フィルタである場合、

【数3】

について、ｘ（ｎ＋ｍ）はｘ’（ｎ＋ｍ）＝ｘ（ｎ＋ｍ）＋／－２πに置き換えられ、

【数4】

となる。

【0041】

この時点で、本発明の方法のいくつかの特徴的な機能を強調する。

【0042】

第１に、位相アンラップモジュールでは、出力信号がサポートされている範囲を超えると、変化方向に依存して、Ｎ・２π（Ｎは負または正の整数）を加算することによって、出力信号を範囲内に戻す。

【0043】

第２に、調整された位相は、フラグ信号と共にＨＰＦに出力される。

【0044】

第３に、ＨＰＦにおけるバッファされたｘ（ｎ－ｍ）（ｍ＝０，１，２，．．．）をｘ’（ｎ－ｍ）＝ｘ（ｎ－ｍ）＋Ｎ・２πに置き換えて、追加工程によって生じるスパイクとそれに続くリップルを除去する。

【0045】

図３（Ａ）および図３（Ｂ）は、本開示の態様による方法の例示的な処理ステップを示す一対のフローチャートである。

【0046】

これらの図を参照すると、図３（Ａ）では、差動ビートからの入力信号ｓｉｇ（ｎ）を使用して、関数にａｔａｎ（）および位相アンラップが含まれていることが確認できる。入力ｓｉｇ（ｎ）は複素数型であり、そこから

【数5】

を用いて角度を計算し、（－π，π］の範囲に変換する。

【0047】

角度出力に続いて、従来技術で知られている標準的な位相アンラップを行い、ｘ（ｎ）－ｘ（ｎ－1）≧（－π，π］とする。そのアンラップに加えて、本開示の方法は、出力位相がオーバーフローするかどうかをチェックする。正の最大値を超えた場合は、出力値が２π減算され、負の最大値を超えた場合には、出力値が２π増加される。あるいは、出力位相が（－π，π］の範囲に収まるように、適用される調整をＮ・２πにすることもできる。この更新されたｘ（ｎ）および調整Δ（ｎ）（フラグ信号で表される）は、ハイパスフィルタに出力される。

【0048】

図３（Ｂ）は、ハイパスフィルタの動作を示すフローチャートである。従来技術と比較したその動作の際立った特徴の１つは、本開示による本発明の方法では、オーバーフローが発生すると、フィルタ内の前のサイクルからのバッファされた全ての入力信号（すなわち、ｘ（ｎ－ｉ）、ｉ＝１，２,．．．，Ｎ）がΔ（ｎ）によって調整されることである。

【0049】

ＤＡＳにおけるレイリー後方散乱信号のコヒーレント検出では、まず、Ｘ（ｃ，ｎ）・Ｘ（ｃ－ｌ，ｎ）*を計算することによる差動ビートを使用し、ここで、ｌは「ビートタップ」と呼ばれ、センサファイバの長さに沿った２つの位置間のサンプル数である。このビート積の位相は、２つの位置の間の時間ｎにおける全応力に対して線形である。

【0050】

図４は、従来技術で行われている位相アンラップを示す概略図であり、２つの信号ごとの差がπを超えないという仮定に基づくものである。まず、ａｔａｎ（）を用いて複素積から（－π，π］の範囲の角度φ（ｎ）を計算し、次に、前回の角度（ｚ^-1、シングルタップ遅延線、またはレジスタからの(φ（ｎ－１）)で減算し、必要に応じて＋２πまたは－２πのいずれかで差分を（－π，π］の範囲内に変更する。そして、その差分を前回のアンラップ出力ｘ（ｎ－１）に加算して出力ｘ（ｎ）とする。この処理により、ｘ（ｎ）－ｘ（ｎ－1）が（－π，π］の範囲に収まる。

【0051】

デジタル信号処理ではビット数が限られているため、アンラップ信号ｘ（ｎ）がオーバーフローする可能性がある。このようなオーバーフローは、信号の飽和（すなわち、ハイパスフィルタを通過した後の出力がゼロ）を引き起こすか、または正常な範囲に強制的に戻された場合に強いスパイクおよび長いリップルを発生させる。

【0052】

本開示による方法は、オーバーフローが発生したときに、アンラップ信号を正常な範囲に戻す。

【0053】

理解を深めるために、サポートされる範囲を［ｎｅｇ＿ｍａｘ，ｐｏｓ＿ｍａｘ］とする。ここで、ｎｅｇ＿ｍａｘ＜０、ｐｏｓ＿ｍａｘ＞０であり、アンラップ位相はサンプルｎでオーバーフローする。すなわち、－π＜ｘ（ｎ）－ｘ（ｎ－１）＜π、ｎｅｇ＿ｍａｘ≦ｘ（ｎ－１）≦ｐｏｓ＿ｍａｘであり、ｘ（ｎ）＞ｐｏｓ＿ｍａｘまたはｘ（ｎ）＜ｎｅｇ＿ｍａｘである。

【0054】

本開示による方法は、ｘ’（ｎ）＝ｘ（ｎ）－Δ（ｎ）を使用し、ここで、

【数6】

である。

【0055】

ここでは、Ｎ・２π≦ｐｏｓ＿ｍａｘ＜（Ｎ＋１）・２πである。出力は、ｘ（ｎ）をｘ’（ｎ）で置き換える。

【0056】

図５は、本開示の態様による、オーバーフローが発生したときのアンラップ位相処理を示すプロットである。この図の例のように、ｘ（ｎ）がｐｏｓ＿ｍａｘより大きくなると、出力は、ｘ（ｎ）＝ｘ（ｎ）－Δ（ｎ）を使用する。ここで、２つの選択肢で示されるように、Δ（ｎ）＝２πまたはΔ（ｎ）＝Ｎ・２πである。この調整された信号ｘ（ｎ）は、オーバーフローとオーバーフローの方向（正または負）を示すフラグｆ（ｎ）と共に出力される。例えば、f（ｎ）は２ビットを使用し、１ビット目はオーバーフロー（「ｂ１」）するか否か（「ｂ０」）を示し、２ビット目は正（「ｂ０」）方向または負（「’ｂ１」）方向を示すことができる。ｘ（ｎ）とｆ（ｎ）の両方がＨＰＦブロックに入力される。

【0057】

ＨＰＦは、ＩＩＲまたはＦＩＲのいずれかを用いることができる。ＩＩＲの場合は、複数の２次ＩＩＲフィルタをカスケード接続することによって、任意の高次フィルタを実装できることが知られている。２次ＩＩＲフィルタは、以下のように表すことができる。

【数7】

【0058】

２次ＩＩＲフィルタの実装には、直接型Ｉ（ＤＦ－Ｉ）または直接型ＩＩ（ＤＦ－ＩＩ）の２つの選択肢がある。ＤＦ－Ｉ実装には、最終結果ｙ（ｎ）が範囲内にある限り、内部フィルタのオーバーフローの可能性がないという利点がある。これにより、ハードウェアの実装により適したものとなり、これは本発明で採用される選択肢である。

【0059】

図６は、上記の式を直接写像したＤＦ－Ｉフィルタの実装である。

【0060】

アンラップ位相がオーバーフローに遭遇すると、本発明では、そのバッファリングされたｘ（ｎ－１）およびｘ（ｎ－２）にΔ（ｎ）を減算することによって、第１の２次ＩＩＲのフィードフォワード段における動作を変更する。すなわち、ｘ’（ｎ－１）＝ｘ（ｎ－１）－Δ（ｎ）、ｘ’（ｎ－２）＝ｘ（ｎ－２）－Δ（ｎ）、および

【数8】

である。

【0061】

フィードバック段および後続の２次ＩＩＲフィルタはどちらも変更されない。この方法を、図６の従来技術と比較して、チェーン内の第１の２次ＩＩＲフィルタについて図７に示す。

【0062】

ＦＩＲフィルタについても同様に、

【数9】

であり、ｘ（ｎ）がオーバーフローする場合、

【数10】

について、ｘ（ｎ＋ｍ）は、ｘ’（ｎ＋ｍ）＝ｘ（ｎ＋ｍ）－Δ（ｎ）に置き換えられ、

【数11】

となる。

【0063】

図８は、図２の入力と同じデータを使用する、この開示の上記方法からの本開示の態様による改善されたＨＰＦ出力を示すプロットである。この図から容易に分かるように、スパイクおよびリップルは両方とも除去される。

【0064】

図９は、本開示の態様による、インタロゲータとコヒーレント受信機と解析システムとを含む例示的なＤＦＳシステムの概略ブロック図である。

【0065】

いくつかの具体例を用いて本開示を示したが、当業者であれば、本教示がそのように限定されないことを認識するのであろう。したがって、本開示は、本明細書に添付される特許請求の範囲によってのみ限定されるべきである。

【図1】

【図2】

【図3(A)】

【図3(B)】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】