JP2024510957A - 分散型光ファイバセンシングにより地震情報を提供する装置、システム、および方法 - Google Patents
分散型光ファイバセンシングにより地震情報を提供する装置、システム、および方法 Download PDFInfo
- Publication number
- JP2024510957A JP2024510957A JP2023554858A JP2023554858A JP2024510957A JP 2024510957 A JP2024510957 A JP 2024510957A JP 2023554858 A JP2023554858 A JP 2023554858A JP 2023554858 A JP2023554858 A JP 2023554858A JP 2024510957 A JP2024510957 A JP 2024510957A
- Authority
- JP
- Japan
- Prior art keywords
- dfos
- data
- optical fiber
- earthquake
- further embodiment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 311
- 239000000835 fiber Substances 0.000 title claims description 67
- 239000013307 optical fiber Substances 0.000 claims abstract description 172
- 230000001133 acceleration Effects 0.000 claims abstract description 80
- 101001059930 Drosophila melanogaster Transcription factor kayak, isoforms A/B/F Proteins 0.000 claims abstract 87
- 101001059931 Drosophila melanogaster Transcription factor kayak, isoforms D/sro Proteins 0.000 claims abstract 87
- 238000005259 measurement Methods 0.000 claims description 58
- 230000033001 locomotion Effects 0.000 claims description 51
- 230000003287 optical effect Effects 0.000 claims description 36
- 230000000694 effects Effects 0.000 claims description 19
- 238000001228 spectrum Methods 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 18
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000013461 design Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 7
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 230000007423 decrease Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 230000036541 health Effects 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims 1
- 238000004891 communication Methods 0.000 description 18
- 230000008569 process Effects 0.000 description 15
- 230000010287 polarization Effects 0.000 description 10
- 238000001514 detection method Methods 0.000 description 7
- 230000015556 catabolic process Effects 0.000 description 6
- 238000006731 degradation reaction Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 239000002689 soil Substances 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 3
- 230000001131 transforming effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 230000002498 deadly effect Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000001069 Raman spectroscopy Methods 0.000 description 1
- 101100173586 Schizosaccharomyces pombe (strain 972 / ATCC 24843) fft2 gene Proteins 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000003760 hair shine Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000005527 soil sampling Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/01—Measuring or predicting earthquakes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/22—Transmitting seismic signals to recording or processing apparatus
- G01V1/226—Optoseismic systems
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
地表面最大加速度を決定する方法は、少なくとも1つの光ファイバに接続されたDFOS機器を提供することと、DFOSデータをDFOS機器で記録することと、を含む。DFOSデータは、少なくとも1つの光ファイバに沿った歪みデータを含む。また、本方法は、歪みデータを複素フーリエ係数に変換することと、複素フーリエ係数をスケーリングすることと、スケーリングされたフーリエ係数に逆変換を適用することと、少なくとも1つの光ファイバに沿った位置の地表面最大加速度を特定するために、逆変換の出力から最大値を選択することと、を含む。
Description
関連出願の相互参照
本出願は、2021年3月7日に出願された米国仮特許出願第63/157,789号の優先権の利益を主張するものであり、参照することによって同特許出願の全体が本明細書に組み込まれる。
本出願は、2021年3月7日に出願された米国仮特許出願第63/157,789号の優先権の利益を主張するものであり、参照することによって同特許出願の全体が本明細書に組み込まれる。
2011年3月11日に日本の東北沖で発生したM9.0の巨大地震、2011年にニュージーランドのクライストチャーチで発生した地震から10年が経過したが、地震の事象は依然として、世界の人口密集都市にとって急を要する課題である。地震は人命、インフラ両方に多大なリスクをもたらす。地震による初期震動の間に壊滅的な構造崩壊が起こった後、津波、地滑り、液状化、余震といった二次災害が致命的に連続して起こることが多い。多くの場合、これらの二次災害は初期振動よりも致命的である。インフラの損害に関連する大きな経済的損失は、住居や建造物の所有者だけでなく、地震保険を提供する民間・国営金融機関にも持続不可能な負担を強いる。
従来は、地震計の疎なネットワークを利用して地震による地震動の程度を定量化していた。都市周辺地域や、人間の居住地もしくは勤務地から離れた場所に設置されることの多い少数の地震計を、都市全体について一般化することが通常である。実際には、地震波の方向性だけでなく、地域の地質、堆積物の厚さ、および土壌条件によって、地震動は都市全域で劇的に変化する。疎な地震計ネットワークから得られる地震動に関する現行の情報は、番地(address)単位で望ましい解像度の情報を提供するには不十分である。
大地震後の地震計からの情報は不完全であることから、大地震後には技師や地震学者が派遣されて損害の程度を判断するために番地を調査する。多くの場合、建造物内外の壁のひび割れや他の損害の視覚的証拠を探す。こうした調査は著しく非効率的で、交通システムや重要なインフラの復旧を1日以上遅らせることも少なくない。多くの建造物は調査されず、よって損害の記録の立証責任は建造物の所有者が負うこととなる。
さらに、このシステムは大地震(例えば、マグニチュードM>6.5~7.0)を評価することを意図したものだが、マグニチュードが比較的小さい(M2.0~M6.5)地震が頻繁に発生した結果、軽~中度の振動が反復して発生した場合、建造物の長期的な構造的完全性を把握することが困難となり、また保険の観点から地震の損害賠償請求と損害との因果関係を確立することが困難となる。
本開示の実施形態は、地震動の振幅の程度、危険性が高い領域、津波や液状化などの二次的な自然災害の有無など、地震の影響の解析を可能とする方法、システム、および装置を提供する。
当業者は、地震学に関して、今日、検知は主たる関心事ではないことを理解するであろう。地震が発生したことを検知する方法は数多くある。例えば、モバイル機器から外部の地震サービスを利用して地震が発生したことを確認することができる。損害をもたらす地震が発生した直後、最も価値のある情報は地震の検知ではなく、地震の衝撃である。特に、実施形態は、地震動の種類、類型、程度、および当該地震動が建造物や他のインフラにどのような影響を与えたかを特徴付ける方法、システム、および装置を提供する。この観点から、実施形態では、人口密集地域において、個別の番地(建造物や指定された地域など)の評価を可能とする空間分解能(いわゆる「番地別地表面最大加速度(peak ground acceleration by address)」)で地震リスクの予測を改善する方法、装置、およびシステムを提示する。
実施形態では、方法、装置、およびシステムは、高分解能(例えば、1メートルの分解能)で光ファイバ(通信ファイバインフラなど)に沿って行われる分散型光ファイバセンシング(Distributed Fiber Optic Sensing:DFOS)の測定値を利用し、地震エネルギーによって引き起こされる潜在的な建造物の損害、構造的完全性、液状化、二次災害の可能性、および経済的損害に関する定量的情報を提供する。
以下、添付図面を参照して実施形態を詳細に説明する。同様の参照番号は同様の要素を表す。添付の図面は必ずしも原寸に比例しない。一部の図は、他の基本的な特徴をより明確に示す目的で、選択された特徴を省略することによって簡略化してある。一部の図における要素のこうした省略は、対応する書面の説明に明示的に開示されている場合を除いて、例示的な実施形態のいずれかにおける特定の要素の存在もしくは非存在を必ずしも示さない。
図1を参照すると、DFOS機器100は光ファイバ150に光を照射し、時間の関数として戻ってくるエネルギー110を記録する。光は、単純なパルス、チャープパルス、連続波など、どのようなものでもよい。実施形態において、光はレーザ光である。さらなる実施形態において、レーザ光は赤外または近赤外の周波数範囲にある。光ファイバセンシング経路から戻ってきたエネルギーは、光ファイバの長さに沿った光散乱特性の結果としてレイリー散乱されている場合もあり、あるいは戻ってきたエネルギーは、入射波長からのブリルアン遷移またはラマン遷移によって生じた場合もある。戻ってきた光は光学的、デジタル的、あるいはその両方で解析できる。このプロセスの出力は、ここではDFOSデータまたはDFOS記録と呼ばれるデータセットであり、該データセットはファイバのすべてのセンサ位置における特定の時間サンプルの光ファイバの状態に関する値を含む。
実施形態では、ファイバ長さの時間変化率に関する情報を伝達し得る、偏光特性や一端から他端への飛行時間など、ファイバの状態、あるいは、ある点または長さ全体にわたるファイバの応力もしくは歪みの状態に関連する情報を抽出するため、電気通信レシーバの統計値そのものが、一方の端に専用の測定器がなくても、DFOSの入力日付として収集される。歪みには測地歪みが含まれ、これは測地歪みレスポンスによって測定できることが理解されよう。地震の後、震源地付近の地表は、おそらくは横方向に、あるいは断層を挟んで上下に動く。この動きは測地歪みと呼ばれる。
一実施形態において、DFOS機器は、WO2018/045433A1(参照によりその全体が本明細書に組み込まれる)に記載される分散音響センシングのためのユニットなどのインテロゲータユニット(「IU」)であってもよい。
一部のDFOSの場合、物理的な固定長のセンサ位置を有さない。原理的には、ゲージ長は、1つのDFOS値が感度を有する、光ファイバに沿った距離である。例えば、ある種類のパルスレイリーベースの(pulsed Rayleigh-based)DFOSでは、ゲージ長は10mで、記録されたデータは歪みであってもよい。この場合、10mのゲージ長は「基準長」となり、それより大きい変位は歪みの原因となる。ゲージ長は、測定の品質、測定のゲイン、空間エイリアシングなしに解析可能な最も細かい空間サイズ、記録されるDFOSデータ量など、測定の多くの側面に影響を与え得る。ゲージ長は、ハードウェアおよび/またはソフトウェアで設定できる。ゲージ長がソフトウェアで設定されている場合、複数のゲージ長でDFOSデータを記録できる。
DFOS測定は、光ファイバの動きや変形、ひいては光ファイバの周囲の動きや変形の影響を受ける。例えば、地震の際、地震波が地表近くを伝播する際に土壌が圧縮され粗密化し、この動きが光ファイバに伝達される。
DFOS測定は、連続する光ファイバが存在し、その一端をDFOS機器に接続できる場所であればどこでも行うことができる。DFOS測定の各ゲージ長に複数の運動/変形成分を導入するため、ファイバは任意の向きに敷設してもよく、また中央のシリンダにらせん状に巻き付けてもよい(例えば橋脚に巻き付け、あるいは建造物の基礎に巻き付けてもよい)。別の目的で敷設された既存の光ファイバをDFOS測定に利用してもよい。複数のファイバを直列に接続し、1台の装置でDFOSに使用してもよく、また(同一もしくは別々のDFOS機器で解析される)複数のDFOSチャネルを使用して、同じ付近のDFOSデータを記録してもよい。
ファイバ経路に沿って行われたDFOS測定は、地震の地震動から生じる変位、速度、加速度、および/または歪みに関する定量的な情報を提供し、よって建造物の損害、液状化、構造的完全性に生じる危険性を直接定量化するデータを提供する。
各種実施形態において、DFOS測定は、地震波が特定の場所に到達する前(例えば、数秒前)、地震動の間、および地震動の後に行われる。特定の場所は、人口密集地における場所であってもよく、空間分解能は、特定の建造物を特定できるように、メートルオーダー(例えば、1メートル)であってもよいことが理解されよう。
地震は多くの場合、人間の居住地やインフラから少し離れた場所で発生する。例えば、日本やニュージーランド、カナダのブリティッシュ・コロンビア沖で発生した大地震のほとんどは、沿岸都市から100km以上離れた場所で発生している。実施形態では、通信用に海上や陸上に配備されている長尺(例えば、数百キロメートルの長さ)のファイバも、ファイバ内の光速で進む光信号によってファイバに沿った任意の地点でこれらの遠隔地震を検知するために使用できる。
実施形態では、DFOS機器は陸上の都市データセンターに設置され、長い(例えば100km)沖合のファイバ経路に接続される。地震が100km沖合、ファイバの海側端部で発生した場合、DFOS測定の使用は、DFOSの使用が地震の早期警報を提供するというシナリオを提示し得る。この例では、地震の地点から地球を最も速く伝わるP波(縦波)によって海底ファイバが変形し、この海底ファイバの変形がDFOS測定によって検知され、定量化される。光ファイバのセンシング長を100kmとすると、システムは、P波が都市に到達する前に、16.7秒間の地震警報を伴う警報を発令することができる(P波の速度は6km/sを想定)。最も大きい損害をもたらすレイリー波や表面波の到着は、P波に比べて相対的に速度が遅いことから、わずかに遅れる。そのため、この方法では、人口密集地域に対して最大50秒間の警報が可能となり得る(レイリー波の速度を2km/sと仮定)。よって、地震の震源地まで、あるいは震源地を越えて延在する既存の、あるいは新たに敷設された光ファイバケーブルによるDFOSを採用することで、地震の発生を検知すること、および地震によるP波や表面波が都市に到達する前に、人口集中地区に警報を提供することが可能であることが理解されよう。さらに、震央位置の地盤状態を直接観察すること(通信ファイバケーブル上にDFOSを大規模に配備することで、その期待は向上している)で、巨大地震が起こる前にその可能性を予測するために使用できる高密度の地質工学データを提供できることが理解されよう。
実施形態では、DFOS測定値は、(例えば、地表面最大加速度(PGA)の観点から)地震動の類型、方向、および/または振幅および継続時間に関する情報を提供するために使用される。地震学や地震科学、あるいは信号処理の技術に精通した当業者には、提案された方法で得られる測定振幅は、都市部の地震記録の履歴に基づいて、地震のタイミング中に最大となることが明らかであろう。故に、地震検知スキームがない場合でも、地震動の影響について知らせるための測定方法を提供する本開示によるDFOS測定を利用することが可能である。
このことは、この方法を一般に使用する場合にも、単独地震または連続地震という特定の場合に使用する場合にも当てはまる。例えば、一般的なユースケースが何を意味するかを明確にするために、ニュージーランドの地震が起こりやすい地域にある建造物の所有者や、建造物資産の保険業者を考えてみたい。上記所有者や保険業者は、壊滅的な損害をもたらすような事象はなくとも、低レベルの地震が毎年頻繁に発生していることから、近隣や自らの建造物の地表面最大加速度(PGA)の統計分布に関心があるかもしれない。この場合、DFOS測定値を連続的に記録・解析し、対象建造物近傍のファイバのセグメントの平均PGAレベル、または平均と標準偏差、もしくは別の統計的尺度を算出してもよく、このときには、何年にも及ぶ可能性がある一定期間にわたる統計的見地のデータを確実に計算する必要がある。この方法の特定のユースケースが何を意味するかを明確にするために、同一または異なる建造物の所有者もしくは保険業者が、損害をもたらす1回の本震、あるいは1回の本震・余震における本震とそれに続くすべての余震の際のPGA値をどのように知りたいと考え得るかを考えてみたい。PGAのすべての測定値はDFOSから得ることができる。この説明から明らかなように、計算中のPGAの測定は連続的に(例えば1分毎に)実行され得ることから、一般的なユースケースと特定の単独地震のユースケースの両方を充足できる。
実施形態では、DFOSデータを使用して加速度(A)を推定するため、ファイバアレイの複数の位置におけるDFOSからの歪みが記録される。実施形態では、複数の位置はすべてファイバアレイ内の位置である。PGAは、下記の式1の関係から求められる。
式中、uは光ファイバケーブルに沿った水平方向の地盤(地下媒体)の見かけの速度の測定値であり、Aは加速度(例えば、単位はm/s2)である。
式中、uは光ファイバケーブルに沿った水平方向の地盤(地下媒体)の見かけの速度の測定値であり、Aは加速度(例えば、単位はm/s2)である。
DFOSデータは距離および時間において一様にサンプリングされ得る。実施形態では、DFOSデータは空間的に約1~10mで、毎秒100~1000サンプルの細かい時間増分でサンプリングされる。その結果、2次元フーリエ変換(FFT2)を用いて、到着した地震から任意の特定の地震位相を、加速度の計算が可能な量に変換できる。よって、前式は次のように書き換えられる。
式中、εはDFOS法で記録された歪み値を表す。歪み値は、通常100~2000Hzの範囲のアナログ-デジタルコンバータによって設定されたサンプリングレートで、時間tでサンプリングされたデジタル化された記録である。光ファイバケーブルの各セグメントまたはゲージ長は独立したセンサとして機能し、よって歪み場も光ファイバケーブルに沿って、このゲージ長の選択によって決定される空間サンプリングレートで空間xにサンプリングされる。一般的なゲージ長は10メートルである。ゲージ長は、ハードウェアで予め決めておいてもよく、ソフトウェアで予め決めておいてもよく、あるいは、生の光位相データを保存し、様々なゲージ長の後処理を可能としてもよい。
次いで、上記のアプローチを適用する。
式中、εはDFOS法で記録された歪み値、すなわち空間におけるDFOSチャネルの座標であるx、および記録の時間サンプルであるtの両方に依存する量を表し、Eはεに2次元フーリエ変換を適用した結果であり、Eは変換された周波数成分ωおよびkを有する。ωは時間周波数、kは空間周波数を表す。
式中、εはDFOS法で記録された歪み値を表す。歪み値は、通常100~2000Hzの範囲のアナログ-デジタルコンバータによって設定されたサンプリングレートで、時間tでサンプリングされたデジタル化された記録である。光ファイバケーブルの各セグメントまたはゲージ長は独立したセンサとして機能し、よって歪み場も光ファイバケーブルに沿って、このゲージ長の選択によって決定される空間サンプリングレートで空間xにサンプリングされる。一般的なゲージ長は10メートルである。ゲージ長は、ハードウェアで予め決めておいてもよく、ソフトウェアで予め決めておいてもよく、あるいは、生の光位相データを保存し、様々なゲージ長の後処理を可能としてもよい。
次いで、上記のアプローチを適用する。
式中、εはDFOS法で記録された歪み値、すなわち空間におけるDFOSチャネルの座標であるx、および記録の時間サンプルであるtの両方に依存する量を表し、Eはεに2次元フーリエ変換を適用した結果であり、Eは変換された周波数成分ωおよびkを有する。ωは時間周波数、kは空間周波数を表す。
図2は、本開示の主題の実施形態に従った、地表面最大加速度を計算する方法の非限定的な例を示す。S200で、プロセスは開始する。S210で、DFOSデータが記録され、保存される。DFOSデータは歪みデータを含み、歪みデータは歪み振幅を含んでもよい。次いで、S220で、2次元フーリエ変換を適用して、歪みデータを複素フーリエ係数に変換する。次に、S230で、フーリエ係数に量ω/kを乗ずるが、これには式4のように元のデータを見かけの速度でスケーリングする効果がある。次に、S240で、2次元逆フーリエ変換を用いてデータを時空間領域に戻し、時間微分するとともに-1(否定)を乗ずることで、式5のような加速度Aの単位を持つデータセットが得られる。S250で各チャネルにおけるAの最大値が選択され、光ファイバ内の当該位置の地表面最大加速度とみなされる。
他の実施形態では、地表面最大加速度を計算する代替方法が説明される。この代替方法は、震源のおおよその位置(緯度、経度、DFOS測定位置からの半径距離に応じて5km以内の深さ)を用いてレイトレーシングによりuを推定するために用いられる地域地震速度モデルと、ファイバの領域におけるせん断波速度(Vs)の推定値とを使用する。地震のおおよその位置は、地震計やその他の非DFOSセンサから得ることができるが、DFOSセンシングは、対象地域や対象地域内の特定の物体(建造物、橋、道路など)に対する地震の影響を定量化するために使用される。Vsの値は、地質調査、土壌サンプリング、写真測量、もしくは航空画像から知ることができ、あるいは地質図やその他の地質情報のアーカイブから得られる情報より知ることができる。これらの実施形態では、レイトレーシング中に計算されたDFOSアレイの領域で地震波が水平面と成す代表的な角度(θ)が保存される。これらの代表的な角度を用いて、例えばuは以下の式で推定される。
次いで、この関係を用いて加速度値を計算する。この関係を用いて、式1のように加速度を計算する。
次いで、この関係を用いて加速度値を計算する。この関係を用いて、式1のように加速度を計算する。
図3は、地域地震速度モデルに基づいて地表面最大加速度を計算する方法の非限定的な例を示す。プロセスはS300で開始する。DFOSデータは連続的に記録され、S310で保存される。実施形態では、データは循環バッファに記録され、新しいデータで古いデータを上書きする。これにより、ストレージ要件を削減しながら、関連データを解析することができる。S320で地震の発生が検知される。上述したように、この検知は地震計による測定など、DFOS以外の測定によって行うことができる。この検知は、S310で記録された解析DFOSデータをトリガする。
次に、S330で、地震波線経路がレイトレースされる。これには地震の震源地と地域の地質情報を入力する必要があり、それらは例えば米国地質調査所(U.S.Geological Survey:USGS)によって公開されているようなこの情報の外部アーカイブからロードできる。レイトレーシングは、スネルの法則、すなわち地震波が、1つの物質層から次の物質層へと通過する際に、層の境界面に対する入射角度に応じて屈折する方法を応用したものである。対象のDFOSファイバチャネルが存在すればどこでも、レイトレーシングの結果として表面における入射角度θの推定値が得られる。
次に、S340で、S330で得たθの値を、DFOSファイバチャネル付近の地下のせん断波速度Vsとともに用いて、uを式6のように計算する。実施形態では、DFOSファイバ付近の地下とは、DFOSファイバの5m、10m、100m、または500m以内である。
最後に、S350で、記録された歪みのDFOS測定値に時間微分が適用され、その値uが乗算されて加速度A単位でデータセットが得られる。各チャネルにおけるAの最大値を、S360で選択された光ファイバ内の当該位置の地表面最大加速度とする。
病院や、橋・トンネルなどの重要なインフラが、建造物の振動を生じることが知られている特定の周波数で強い振動を記録した場合、電気の通電停止、ガスの自動遮断弁、エレベーターのドア開放といった追加の安全プロトコルが実行できる。このことは、DFOSに使用されるファイバが存在する場所ならば、建造物毎および住所毎に達成され得る。
図4Aを参照すると、光ファイバ信号強度の劣化に基づいて地震の影響を検知し定量化する方法の非限定的な例が示されている。S410で、DFOSデータが記録され、保存される。保存されたデータは、S415で、地震の優位周波数範囲(局部地震では0.5~150Hz、広域地震および直下型地震では0.005~5Hz)の通過帯域内における信号振幅の定量的測定値を使用して解析される。この解析の結果は、表面波の放射が逆二乗則に似るという理解に基づいて測定された領域を超えた伝播距離における予想される振幅の推定値となり得る。一部の実施形態では、DFOS振幅測定値を用いて所定の場所または都市で予想される地震動レベルを予測するため、層状地盤の外在的および内在的減衰効果について記載した、Anderson and Hart(1978)によって提案されたモデルなど一般または専用減衰モデルが採用される。
S420で定量的測定値が所定の閾値を上回った場合、地震発生と推定する。そうでない場合、プロセスはループバックする。信号振幅の劣化は、車両の通行、工事、掘削、DFOS用にデータ記録を行っている光ファイバの近傍で付加される他のノイズやエネルギーなど、地震以外の物理的作用によって生じることもある。このように、S430で地震発生の検証が行われる。このことは、地震データを保存している1つまたは複数のコンピュータシステムに照会して、他のセンサで地震が検知されたか否かを判断することを含んでもよい。
S430で地震が確認されると、プロセスはS440に進み、光ファイバのセンシング距離内にある建造物や橋などの特定の場所について、地震の影響が定量化される。影響の定量化は、図2に示された、前述したプロセスに続いて行ってもよい。具体的には、プロセスは記録されたDFOSデータを既に含み、処理はS220に続いてもよい。
他の実施形態では、専用のDFOSファイバが上記の例の病院などの対象位置の近傍に配置されていなくても、地震動の測定と算出は、他の既存のファイバ通信インフラに基づいて行うことができる。次に図4Bを参照すると、以下に強調する相違点を除いて、ステップは図4Aと同様である。
一実施形態では、専用のDFOSインテロゲータユニットを持たない通信ファイバは、ファイバを通じて送受信される信号の品質(例えば、信号対雑音比または散乱の測定値)を継続的に監視すべく構成されている。このプロセスは、S416で信号品質の評価が行われる点を除いて、図4Aのプロセスと同様である。この評価は連続的であってもよく、DFOSシステム100によって受信され記録された光信号の信号対雑音比の測定を含んでもよい。品質測定値(例えば、SNR)は、S421で閾値と比較される。品質が所定のレベルまで低下すると、地震の地震エネルギーなどによりファイバの物理的変化が生じたと判断され、DFOSに使用可能なデータが記憶装置に記録される。記録は、信号品質検知の前後に行うことができる。実施形態では、信号品質低下が検知される前に、より少ない量のデータがS410で保存される。信号品質低下の検知後、より多い量のデータが解析用に保存される。これにより、(道路工事やその他の事象ではなく)実際に地震が信号品質の低下を引き起こした場合に記録データの解析を行うことが可能となり、一方でローカルストレージ容量が小さく安価に維持されているため、すべてではないにしても多くの光ファイバ通信機器がこのデータ収集を実行できる。品質の低下が地震により生じたのであれば、光ファイバ通信機器に保存されたデータが解析され、通信ファイバに沿った地震動の定量的測定値が決定される。
実施形態では、ファイバケーブルから100m以内の対象位置の地震動を正確に計算できるように、100mの空間分解能が達成可能である。この距離を超えると、土壌の状態が変化してファイバの測定精度が低下することがある。DFOSによる地震動測定値の密度は、地震時に地震動が最小となる都市の安全な場所を計画するのに利用できる。よって、DFOSによる地震動の測定値を場所毎に索引付けして保存する、地震動計算値のデータベースが作成される。これにより、(自然の地盤特性によるものであれ、特定の人為的緩和策によるものであれ)危険な地震動が起こらない場所を特定することが可能となる。
ここで図5を参照すると、図示されたプロットは、開示された主題の実施形態に従って実行されたDFOS記録からの地震記録セクションを示す。濃い色は強い地震動を示す。白い線は、地震波を発生させる断層運動が発生した震源時である。地震エネルギーのパルスは、少し経つと、X軸の光学距離[km]に沿って表されるDFOSチャネル上で可視的となる。すなわちP波とS波である。最初に到達するエネルギーは、Time=約80秒のP波である。約35~40秒後のTime=120秒の第2のエネルギーパルスがS波または表面波である。ソースとレシーバの幾何学により、これらの波の到着には曲率がある。振幅の変動も大きい。光学距離5~16km付近のチャネルは、その直前もしくは直後のチャネルよりも強い地震動を記録する。これは地盤増幅の結果であり、地表面最大加速度の測定により捕捉される効果である。これは、各々の地域で地表面最大加速度レベルがもたらす予期される危険に応じて設計や建築を行い、あるいはこれらの危険に基づいて建築基準法を制定する仕事をしている都市計画者、建造物所有者、地質工学技術者にとって貴重な情報である。
図6は、2021年9月21日にオーストラリアのメルボルンで発生したM5.9の地震の付近の、複数の方位および距離における例示的な地域展開を示す。ここでラベル付けされたオフセットにある各位置は、この期間に記録されたDFOSデータであり、それぞれのデータは図7A~図7Fに示されている。この地震はメルボルン中心部から東北東に約130kmの地点で発生した。
図7Aから図7Fは、2021年にオーストラリアのメルボルンで発生したM5.9の地震のDFOSデータの測定値を示している。左の縦軸は時間を秒単位で表し、横軸は光学距離をキロメートル単位で表す。各図の右側のスケールは、各図のグレースケールの濃淡に対応するパワースペクトル密度を表す。高振幅の到来エネルギーは、アレイの片側から次の側へと拡散する。ファイバアレイのセグメント間にも明確な差異が存在する。図7Aにおいて、図示されたプロットは、開示された主題の実施形態に従って実行された、オーストラリアのメルボルンにおけるDFOS記録からの地震記録セクションを示す。図7Bにおいて、図示されたプロットは、開示された主題の実施形態に従って実行された、オーストラリアのビクトリア沖におけるDFOS記録からの地震記録セクションを示す。図7Cにおいて、図示されたプロットは、開示された主題の実施形態に従って実行された、オーストラリア沖の別の領域におけるDFOS記録からの地震記録セクションを示す。図7Dにおいて、図示されたプロットは、開示された主題の実施形態に従って実行された、オーストラリアのアデレイドの別の領域におけるDFOS記録からの地震記録セクションを示す。図7E、図示されたプロットは、開示された主題の実施形態に従って実行された、オーストラリアのシドニーの別の領域におけるDFOS記録からの地震記録セクションを示す。図7Fにおいて、図示されたプロットは、開示された主題の実施形態に従って実行された、オーストラリアのシドニーの別の領域におけるDFOS記録からの地震記録セクションを示す。図7A~図7Fの各図において、Time=約20秒からTime=380秒までの間に到来する暗い色の水平信号は、地震による地震エネルギーを示す。
図8は、2021年にオーストラリアのメルボルンで発生したM5.9の地震を記録した、メルボルン中心部のファイバの計算結果である。すべての独立したDFOSセンサは、生の歪み単位から加速度単位に変換される固有の地震波形を記録し、ピークサンプルが選択され、示されている範囲内でプロットされる。これにより、地表面最大加速度を、番地別PGA(PGA by address)と呼ばれる番地の細かい粒度で計算できる。
番地別PGA結果を集計し、対象地域の地図にオーバレイすることができる。図9を参照すると、番地別PGAの値がオーバレイされたメルボルン中心部の地図を含むグラフィカルな地震レポートの例が示されている(相対振幅はグレースケールの濃淡で表され、暗いほど振幅が小さく、明るいほど振幅が大きい)。
実施形態では、振動中のDFOS測定値は、建造物内を通る光ファイバケーブルを使用して、フロア毎に建造物内部に関する情報を提供するために使用される。現代の建造物は、通信用の光ファイバケーブルを有することが多い。これらのケーブルには、地震時の構造物の地震動による変形が生じる。シングルポートまたはマルチポートのDFOS機器は、建造物全体を通る単一または多数の光ファイバケーブルの一端に設けることができ、DFOSデータを連続的に記録することができ、もしくは地震が検知された後に処理されるローリングバッファに記録することができる。建造物からのDFOSデータを解析することで、建造物の個々のサブパート(外壁、支柱、床用根太、天井根太など)が地震に受けた影響を明確に定量化することができる。例えば、低層階の層間変位が高層階よりも小さい場合、この種のDFOSデータを用いて相対的な層間変位の大きさを評価することができる。あるいは、垂直方向と水平方向に走るファイバからDFOSデータを記録することで、水平-垂直方向の変形時系列を算出できる。複数の成分をグリッドに記録することで、記録の各成分における歪みの空間微分の計算を通じて構造物の回転を計算できる。実施形態では、これらのサブパートのうち1つまたは複数の特定の最大加速度および/または変位がDFOSデータに基づいて計算され、次いで、許容閾値(建築家または建築技師によって提供される閾値など)と比較されてもよい。
図10は、地震発生の影響に基づく、建造物や橋などのインフラの自己評価プロセスの非限定的な例を示す。一実施形態では、自動建造物自己評価はDFOSシステムによって実行される。プロセスはS1000で開始する。上述したように、DFOSシステムは、建造物内を通る光ファイバケーブルに接続され、DFOSデータを連続的に記録および解析し、建造物構成要素の変位および/または加速度を計算することができる。この処理は、S1005でTと表示される、10分毎などの定期的なスケジュールで実行され得る。実施形態では、Tは1時間、12時間、または24時間である。実施形態では、プロセスは連続的に実行され、すなわち0秒からT期間が経過するまで実行され、プロセスはS1010に進む。
S1010で、DFOSデータが記録され、保存される。記録されるデータは、期間Tを表しても、別の継続時間を表してもよいことが理解されよう。S1020で、S1010で得られたDFOSデータに基づいて、対象構造物(建造物や橋など)の共振ピークが計算される。共振ピークの履歴値はデータストレージに保存される。これらの値は、予め計算された値を含んでいてもよいし、構造物の有限要素モデルシミュレーションから構造物の設計に基づいて計算された工学的パラメータであってもよいし、一般的構造設計であってもよい。S1030で、アーカイブされた値がデータストレージから取り出され、S1040でDFOSに基づく値と比較される。S1040の出力は、算出値とアーカイブ値との差分値であってもよい。S1050で、差分値が閾値と比較される。閾値は、構造物の有限要素シミュレーションに基づいて算出された工学的パラメータであってもよいし、安全規制によって定められた値であってもよい。実施形態では、閾値は適応的であり、履歴測定値に基づいて減少する。これは、累積するものと考えられる、構造物への構造的影響の蓄積の可能性を考慮したものである。実施形態では、システムは、各種サブコンポーネントの加速度および/または変位の許容閾値を含み、算出値を閾値と比較し、1つまたは複数の閾値を超えた場合にS1060で警報を生成することができる。これにより、建造物は少なくとも地震によって引き起こされ得る初期の損害を自己評価することができ、また建造物の安全性を確保するために構造技師による追加調査が必要な建造物を迅速に特定できることが理解されよう。警報は、有線または無線ネットワークを通じて、1人または複数の指定受信者に送信されてもよい。実施形態では、建造物が監視対象の構造物である場合、1人または複数の指定受信者は建造物の居住者である。
振動の後、DFOS測定値は建造物やインフラの構造的完全性に関する地震後警報を提供するために使用できる。例えば、測定された振動が建造物や橋の安全閾値を超えた場合、その建造物やインフラには自動的にレッドタグを付け、人間が不用意に立ち入ったり使用したりしないようにしてもよい。自然固有モードの共振(natural eigenmode resonances)(構造物の振動や捩れの周波数ピーク)を測定し、地震前の共振と比較することができる。この情報は、建造物の損害に関する警報を提供するために使用できる。
二次災害は、地震によって発生した物理的な地表面のプロセスの結果としての、損害をもたらす地震動の後、短時間に相次いで発生することが多い。実施形態では、土壌の振動に関するDFOS測定値は、液状化調査を促進し、もしくは的を絞った戦略的な方法で液状化調査のスケールを洗練してリソースの利用を改善するために利用される。
開示された主題のさらなる実施形態は以下の通りである。第1のさらなる実施形態によれば、地表面最大加速度を決定する方法が提供され、該方法は、少なくとも1つの光ファイバに接続されたDFOS機器を提供することと、前記少なくとも1つの光ファイバに沿った歪みデータを含むDFOSデータを前記DFOS機器で記録することと、前記歪みデータを複素フーリエ係数に変換することと、前記複素フーリエ係数をスケーリングすることと、前記スケーリングされたフーリエ係数に逆変換を適用することと、前記少なくとも1つの光ファイバに沿った位置の地表面最大加速度を特定するために、前記逆変換の出力から最大値を選択することとを含む。第2のさらなる実施形態によれば、前記歪みデータを複素フーリエ係数に変換することは、前記歪みデータに2次元フーリエ変換を適用することを含む、第1のさらなる実施形態の方法が提供される。第3のさらなる実施形態によれば、前記逆変換を適用することは、前記スケーリングされたフーリエ係数に2次元逆フーリエ変換を適用することを含む、第1のさらなる実施形態の方法が提供される。
第4のさらなる実施形態によれば、地域地震速度モデルに基づいて地表面最大加速度を決定する方法が提供され、該方法は、少なくとも1つの光ファイバに接続されたDFOS機器を提供することと、前記少なくとも1つの光ファイバに沿った歪みデータを含むDFOSデータを前記DFOS機器で記録することと、地震の発生を検知することと、前記光ファイバが設置されている地表面における地震エネルギーの入射角度を決定するために、地震発生地点から前記光ファイバまでの地震波経路をレイトレースすることと、前記決定された入射角度と前記光ファイバ周辺の地面のせん断速度に基づいて、地面の見かけの速度の測定値を計算することと、地表面加速度の測定値を得るために、前記計算された見かけの速度の測定値に時間微分を適用することと、前記少なくとも1つの光ファイバに沿った位置の地表面最大加速度を特定するために、前記適用の出力から最大値を選択することとを含む。第5のさらなる実施形態によれば、前記地面の見かけの速度は水平方向における速度である、第4のさらなる実施形態の方法が提供される。第6のさらなる実施形態によれば、前記地震の発生を検知することは、DFOSデータに基づかない、第4のさらなる実施形態の方法が提供される。
第7のさらなる実施形態によれば、地震によって引き起こされる地震動の特性を決定するための方法が提供され、該方法は、少なくとも1つの光ファイバに沿った歪みデータを含むDFOSデータをDFOS機器で記録することと、前記DFOSデータの信号振幅の定量的測定値を得るために、前記記録されたDFOSデータを解析することと、前記信号振幅の前記定量的測定値を所定の閾値と比較することと、前記定量的測定値が前記所定の閾値を上回った場合に、地震が発生したと予備判断することと、前記地震の影響を定量化することとを含む。第8のさらなる実施形態によれば、前記予備判断した後、かつ前記定量化する前に地震が発生したことを検証することをさらに含む、第7のさらなる実施形態の方法が提供される。第9のさらなる実施形態によれば、前記記録されたDFOSデータを解析することは、前記地震の優位周波数の通過帯域内の信号振幅を測定することを含む、第8のさらなる実施形態の方法が提供される。第10のさらなる実施形態によれば、前記優位周波数は、局部地震では0.5~156Hzの範囲である、第9のさらなる実施形態の方法が提供される。第11のさらなる実施形態によれば、前記優位周波数は、広域地震および直下型地震では0.005~5Hzの範囲である、第9のさらなる実施形態の方法が提供される。第12のさらなる実施形態によれば、前記検証することは、非DFOS測定値に基づいて行われる、第9のさらなる実施形態の方法が提供される。第13のさらなる実施形態によれば、前記検証することは、前記記録されたDFOSデータの解析に基づいて行われる、第9のさらなる実施形態の方法が提供される。第14のさらなる実施形態によれば、前記定量化することは、前記少なくとも1つの光ファイバ近傍の地面の地表面最大加速度を決定することを含む、第9のさらなる実施形態の方法が提供される。第15のさらなる実施形態によれば、前記地表面最大加速度を決定することは、前記DFOSデータの前記歪みデータを複素フーリエ係数に変換することと、前記複素フーリエ係数をスケーリングすることと、前記スケーリングされたフーリエ係数に逆変換を適用することと、前記少なくとも1つの光ファイバに沿った位置の地表面最大加速度を特定するために、前記逆変換の出力から最大値を選択することとを含む、第14のさらなる実施形態の方法が提供される。第16のさらなる実施形態によれば、前記歪みデータを複素フーリエ係数に変換することは、前記歪みデータに2次元フーリエ変換を適用することを含む、第15のさらなる実施形態の方法が提供される。第17のさらなる実施形態によれば、前記逆変換を適用することは、前記スケーリングされたフーリエ係数に2次元逆フーリエ変換を適用することを含む、第15のさらなる実施形態の方法が提供される。
第18のさらなる実施形態によれば、地震によって引き起こされる地震動の特性を決定するための方法が提供され、該方法は、少なくとも1つの光ファイバに沿った歪みデータを含むDFOSデータをDFOS機器で記録することと、前記DFOSデータの信号品質を測定することと、前記測定されたDFOSデータの信号品質を所定の閾値と比較することと、前記測定された信号品質が所定の閾値未満である場合に、地震が発生したと予備判断することと、前記地震の影響を定量化することとを含む。第19のさらなる実施形態によれば、前記予備判断した後、かつ前記定量化する前に地震が発生したことを検証することをさらに含む、第18のさらなる実施形態の方法が提供される。第20のさらなる実施形態によれば、前記信号品質を測定することは、前記DFOSデータの信号対雑音比を決定することを含む、第19のさらなる実施形態の方法が提供される。第21のさらなる実施形態によれば、前記信号品質を測定することは、前記記録されたDFOSデータに適用される、第19のさらなる実施形態の方法が提供される。第22のさらなる実施形態によれば、前記信号品質を測定することは、前記DFOSデータを記録することと同時に前記DFOSデータに適用される、第19のさらなる実施形態の方法が提供される。第23のさらなる実施形態によれば、前記定量化することは、前記少なくとも1つの光ファイバ近傍の地面の地表面最大加速度を決定することを含む、第19のさらなる実施形態の方法が提供される。第24のさらなる実施形態によれば、前記地表面最大加速度を決定することは、前記DFOSデータの前記歪みデータを複素フーリエ係数に変換することと、前記複素フーリエ係数をスケーリングすることと、前記スケーリングされたフーリエ係数に逆変換を適用することと、前記少なくとも1つの光ファイバに沿った位置の地表面最大加速度を特定するために、前記逆変換の出力から最大値を選択することとを含む、第23のさらなる実施形態の方法が提供される。第25のさらなる実施形態によれば、前記歪みデータを複素フーリエ係数に変換することは、前記歪みデータに2次元フーリエ変換を適用することを含む、第24のさらなる実施形態の方法が提供される。第26のさらなる実施形態によれば、前記逆変換を適用することは、前記スケーリングされたフーリエ係数に2次元逆フーリエ変換を適用することを含む、第25のさらなる実施形態の方法が提供される。
第27のさらなる実施形態によれば、構造物に対する地震動の影響を評価する方法が提供され、該方法は、少なくとも1つの光ファイバに沿った歪みデータを含むDFOSデータをDFOS機器で記録することと、前記DFOSデータに基づいて前記構造物の共振ピークを計算することと、前記計算された共振ピークを、前記構造物に関連付けられた共振ピーク値をアーカイブする記憶装置に保存することと、前記記憶装置からアーカイブされた共振ピーク値を取り出すことと、前記取り出されたアーカイブ値と前記DFOSデータに基づいて計算された共振ピークとの間の、変化値によって表される変化を検知することと、前記変化値を閾値と比較することと、前記検知された変化が前記閾値を上回った場合に警報を生成することとを含む。第28のさらなる実施形態によれば、前記少なくとも1つの光ファイバを前記構造物の近傍に設置することをさらに含む、第27のさらなる実施形態の方法が提供される。第29のさらなる実施形態によれば、前記少なくとも1つの光ファイバは、前記構造物の内部または下にある、第28のさらなる実施形態の方法が提供される。第30のさらなる実施形態によれば、前記方法はスケジュール通りに繰り返される、第27のさらなる実施形態の方法が提供される。第31のさらなる実施形態によれば、前記スケジュールの期間は5分から48時間の間である、第30のさらなる実施形態の方法が提供される。第32のさらなる実施形態によれば、前記スケジュールの期間は10分、60分、12時間、または24時間である、第31のさらなる実施形態の方法が提供される。第33のさらなる実施形態によれば、前記スケジュールの期間は適応的であり、前記比較の結果に応じて継続期間が減少する、第30のさらなる実施形態の方法が提供される。第34のさらなる実施形態によれば、前記構造物は、建造物、橋、塔のいずれかである、第27のさらなる実施形態の方法が提供される。第35のさらなる実施形態によれば、前記構造物の設計に基づく工学設計パラメータが記憶装置に保存される、第27のさらなる実施形態の方法が提供される。第36のさらなる実施形態によれば、前記記憶装置は、前記構造物に関連付けられた複数の共振ピーク値を記憶し、前記共振ピーク値の各々は、タイムスタンプおよび継続時間のうちの少なくとも1つに関連付けられる、第35のさらなる実施形態の方法が提供される。第37のさらなる実施形態によれば、最後に計算された共振ピークと前記記憶装置に保存された最も古い共振ピーク値との間の変化を表す累積変化値を特定するために、前記検知された変化値を経時的に累積することをさらに含む、第36のさらなる実施形態の方法が提供される。第38のさらなる実施形態によれば、前記地震動は地震によって引き起こされる、第27乃至第37のさらなる実施形態のいずれかの方法が提供される。第39のさらなる実施形態によれば、前記閾値は適応的であり、経時的に低下する、第27のさらなる実施形態の方法が提供される。第40のさらなる実施形態によれば、前記警報を生成することは、有線または無線ネットワークを通じて信号を送信することを含む、第27のさらなる実施形態の方法が提供される。第41のさらなる実施形態によれば、前記警報を生成することは、前記構造物の内部で人間が聞き取れるメッセージを出力することを含む、第27のさらなる実施形態の方法が提供される。第42のさらなる実施形態によれば、前記人間が聞き取れるメッセージは、前記構造物に対する潜在的損害に関する情報を含む、第41のさらなる実施形態の方法が提供される。
第43のさらなる実施形態によれば、地震の震源である第2の場所から離れた第1の場所で地震警報を生成する方法が提供され、該方法は、前記第1の場所から前記第2の場所を取り囲む領域まで延在する光ファイバの、前記第1の場所でDFOSデータを測定することと、前記第2の場所で地震が発生したことを示す表示を前記DFOSデータから計算することと、前記地震が発生したことを示す警報を前記第1の場所で生成することとを含む。第44のさらなる実施形態によれば、前記第1の場所は人口密集地域外にあり、前記第2の場所は人口密集地域内にある、第43のさらなる実施形態の方法が提供される。第45のさらなる実施形態によれば、前記第1の場所は海上である、第44のさらなる実施形態の方法が提供される。第46のさらなる実施形態によれば、前記第1の場所は地方にある、第44のさらなる実施形態の方法が提供される。第47のさらなる実施形態によれば、前記第2の場所は都市部にある、第44のさらなる実施形態の方法が提供される。第48のさらなる実施形態によれば、前記警報を生成することは、前記地震によるP波および/または損害をもたらす地震動が前記第1の場所に到着する前に行われる、第43のさらなる実施形態の方法が提供される。第49のさらなる実施形態によれば、前記DFOSデータを測定することは、DFOS装置から前記光ファイバに光を送信し、前記光ファイバから屈折した光もしくは反射した光を受信することを含む、第43のさらなる実施形態の方法が提供される。第50のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから前記光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第43のさらなる実施形態の方法が提供される。第51のさらなる実施形態によれば、前記DFOSデータを測定することは、歪みを測定することを含む、第43のさらなる実施形態の方法が提供される。第52のさらなる実施形態によれば、前記DFOSデータを測定することは、地震動を測定することを含む、第43のさらなる実施形態の方法が提供される。第53のさらなる実施形態によれば、前記計算することは、前記地震のマグニチュードを計算することを含む、第43のさらなる実施形態の方法が提供される。第54のさらなる実施形態によれば、前記計算することは、前記第2の場所の地理的位置を計算することを含む、第43のさらなる実施形態の方法が提供される。第55のさらなる実施形態によれば、前記計算することは、前記第2の場所の地表面最大加速度を計算することを含む、第43のさらなる実施形態の方法が提供される。第56のさらなる実施形態によれば、前記計算することは、前記第1の場所の地表面最大加速度を計算することを含む、第43のさらなる実施形態の方法が提供される。第57のさらなる実施形態によれば、前記計算することは、前記光ファイバの所定の距離内の場所の地表面最大加速度を計算することを含む、第43のさらなる実施形態の方法が提供される。第58のさらなる実施形態によれば、前記所定の距離は50メートルから500メートルの間である、第57のさらなる実施形態の方法が提供される。第59のさらなる実施形態によれば、前記所定の距離は100メートルである、第5の-第8のさらなる実施形態の方法が提供される。第60のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第43乃至第59のさらなる実施形態のいずれかの方法が提供される。
第61のさらなる実施形態によれば、第1の場所を震源とする地震の影響を定量化する方法が提供され、該方法は、前記地震の間、前記第1の場所からずれた第2の場所でDFOSシステムを用いてDFOSデータを測定することと、前記測定されたDFOSデータを用いて、地表面最大加速度、地表面最大速度、最大地動変位、最大地動歪み(peak ground strain)、最大地動歪み率、最大スペクトル加速度、および/または地震強度の少なくとも1つを計算することとを含む。第62のさらなる実施形態によれば、前記測定されたDFOSデータを記録することと、地震が発生したという通知を受信することと、前記受信された通知に応答して、前記記録されている測定されたDFOSデータを使用して、発生した地震の地震マグニチュード、モーメント、歪み、歪み率、測地変形、位置、深さ、焦点メカニズム、放射パターン、モーメントテンソル、滑りの有限断層領域、滑り速度、および破壊速度のうちの少なくとも1つを計算することとをさらに含む、第61のさらなる実施形態の方法が提供される。第63のさらなる実施形態によれば、前記算出することに基づいて、前記DFOSシステムから警報を生成することと、前記生成された警報を人口密集地域のレシーバに送信することとをさらに含む、第62のさらなる実施形態の方法が提供される。第64のさらなる実施形態によれば、前記DFOSデータを測定することは、DFOS地震センシング用および電気通信トラフィックの伝送用に同時に使用される単一の電気通信ファイバ上でDFOS信号を送受信することを含む、第61のさらなる実施形態の方法が提供される。第65のさらなる実施形態によれば、前記DFOSデータを測定することは前記地震の発生中に行われ、前記計算することは前記地震の発生後に行われる、第61のさらなる実施形態の方法が提供される。第66のさらなる実施形態によれば、前記DFOSデータを測定することおよび前記計算することは、前記地震の発生中に行われる、第61のさらなる実施形態の方法が提供される。第67のさらなる実施形態によれば、前記計算することは、前記第2の場所における物理的影響を表す、第61のさらなる実施形態の方法が提供される。第68のさらなる実施形態によれば、DFOS装置は前記第2の場所に配置される、第61のさらなる実施形態の方法が提供される。第69のさらなる実施形態によれば、前記計算することは、前記第1の場所と前記第2の場所との間に位置する第3の場所における物理的影響を表す、第61のさらなる実施形態の方法が提供される。第70のさらなる実施形態によれば、前記DFOSデータを測定することは、前記DFOSシステムから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含み、前記計算することは、前記光ファイバから所定の距離内に位置する第3の場所における物理的影響を表す、第61のさらなる実施形態の方法が提供される。第71のさらなる実施形態によれば、前記DFOSデータを測定することは、前記DFOSシステムから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第61のさらなる実施形態の方法が提供される。第72のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第61のさらなる実施形態の方法が提供される。第73のさらなる実施形態によれば、前記所定の距離は50メートルから500メートルの間である、第70のさらなる実施形態の方法が提供される。第74のさらなる実施形態によれば、前記所定距離は100メートルである、第73のさらなる実施形態の方法が提供される。第75のさらなる実施形態によれば、前記第2の場所は、建造物、データセンター、病院、空港、重要インフラ、道路、橋、トンネル、住居、ホテル、他の構造物、または地震動の影響を受ける建築物の少なくとも1つの敷地である、第67のさらなる実施形態の方法が提供される。第76のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第61乃至第75のさらなる実施形態のいずれかの方法が提供される。
第77のさらなる実施形態によれば、地震の影響を定量化する方法が提供され、該方法は、少なくとも前記地震の発生中にDFOSデータを測定することと、前記測定されたDFOSデータを用いて、前記地震の発生中の地震動の方向および前記地震の発生中の前記地震動の継続時間の少なくとも1つを計算することとを含む。第78のさらなる実施形態によれば、前記DFOSデータを測定することは、DFOS装置から光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第77のさらなる実施形態の方法が提供される。第79のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第77のさらなる実施形態の方法が提供される。第80のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第77乃至第79のさらなる実施形態のいずれかの方法が提供される。
第81のさらなる実施形態によれば、建造物に対する地震の影響を解析する方法が提供され、該方法は、少なくとも前記地震の発生中にDFOSデータを測定することと、前記測定されたDFOSデータから前記建造物の固有モード共振値を計算することと、前記固有モードから建造物の振動の強さを表す量を決定することとを含む。第82のさらなる実施形態によれば、前記DFOSデータを測定することは、DFOS装置から光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含み、前記光ファイバは、前記建造物から所定距離離れて延在するか、または前記建造物を通過する、第81のさらなる実施形態の方法が提供される。第83のさらなる実施形態によれば、前記DFOSデータを測定することは、DFOS装置から光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第81のさらなる実施形態の方法が提供される。第84のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第81のさらなる実施形態の方法が提供される。第85のさらなる実施形態によれば、前記DFOSデータを測定することは、前記地震によって生じる振動中に行われる、第81のさらなる実施形態の方法が提供される。第86のさらなる実施形態によれば、前記計算することは、前記地震によって生じる振動中に行われる、第81のさらなる実施形態の方法が提供される。第87のさらなる実施形態によれば、前記計算することは、前記地震により生じる振動が終了した後に行われる、第85のさらなる実施形態の方法が提供される。第88のさらなる実施形態によれば、前記決定することは、前記地震の発生後に計算された前記建造物の固有モードを前記地震前の固有モードと比較することを含む、第81のさらなる実施形態の方法が提供される。第89のさらなる実施形態によれば、前記固有モードの変化を前記建造物の構造的健全性の状態変化に関連付けることをさらに含む、第88のさらなる実施形態の方法が提供される。第90のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第81乃至第89のさらなる実施形態のいずれかの方法が提供される。
第91のさらなる実施形態によれば、地震の余波において構造物の安全警報を生成する方法が提供され、該方法は、少なくとも地震の発生中にDFOSデータを測定することと、前記測定されたDFOSデータに基づいて前記構造物の安全スコアを計算することと、前記安全スコアに基づいて安全警報を生成することとを含む。第92のさらなる実施形態によれば、前記安全スコアは、前記測定されたDFOSデータから計算された地表面最大加速度と、前記構造物の地表面最大加速度率(peak ground acceleration rating)との比較に基づく、第91のさらなる実施形態の方法が提供される。第93のさらなる実施形態によれば、前記安全警報は、前記計算された安全スコアの表現を含む、第91のさらなる実施形態の方法が提供される。第94のさらなる実施形態によれば、前記安全警報は複数の重大度レベルに分類され、前記生成することは前記重大度レベルに基づく、第91のさらなる実施形態の方法が提供される。第95のさらなる実施形態によれば、前記生成された安全警報を、重大度レベルに基づいて選択された受信者に送信することをさらに含む、第94のさらなる実施形態の方法が提供される。第96のさらなる実施形態によれば、前記安全警報は、前記構造物が人間の居住者が立ち入ることについて安全であるか否かを示す、第91のさらなる実施形態の方法が提供される。第97のさらなる実施形態によれば、前記DFOSデータを測定することは、DFOSシステムから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第91のさらなる実施形態の方法が提供される。第98のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第91のさらなる実施形態の方法が提供される。第99のさらなる実施形態によれば、前記光ファイバの物理的変化の検知に応答して、前記測定されたDFOSデータを記憶装置に保存することをさらに含む、第98のさらなる実施形態の方法が提供される。第100のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第91乃至第99のさらなる実施形態のいずれかの方法が提供される。
第101のさらなる実施形態によれば、建造物の損害を検知する方法が提供され、該方法は、前記建造物の内部もしくは近傍の光ファイバでDFOSデータを測定することと、前記測定されたDFOSに基づいて、建造物の層間変位、前記建造物の下部もしくは側部の地表面最大加速度、前記建造物自体の構成要素の地表面最大加速度、および液状化のうち少なくとも1つを表す物理量を計算することと、前記計算された物理量を前記建造物の仕様と比較することと、前記計算された物理量と前記建造物の仕様とに基づいて、潜在的な建造物の損害を示す警報を生成することとを含む。第102のさらなる実施形態によれば、前記DFOSデータを測定することは、DFOSシステムから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第101のさらなる実施形態の方法が提供される。第102のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第101のさらなる実施形態の方法が提供される。第104のさらなる実施形態によれば、前記光ファイバは前記建造物の下方に延在する、第102乃至第103のさらなる実施形態のいずれかの方法が提供される。第105のさらなる実施形態によれば、前記光ファイバは前記建造物を貫通して延在する、第102乃至第103のさらなる実施形態のいずれかの方法が提供される。第106のさらなる実施形態によれば、前記光ファイバは前記建造物の複数の階層を通って延在する、第105のさらなる実施形態の方法が提供される。第107のさらなる実施形態によれば、前記光ファイバは前記建造物の基礎の周囲を囲繞する、第106のさらなる実施形態の方法が提供される。第108のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第101乃至第107のさらなる実施形態のいずれかの方法が提供される。
第109のさらなる実施形態によれば、インフラの損害を検知する方法が提供され、該方法は、前記インフラの内部もしくは近傍の光ファイバでDFOSデータを測定することと、前記測定されたDFOSに基づいて、建造物の層間変位、前記建造物の下部もしくは側部の地表面最大加速度、前記インフラ自体の構成要素の地表面最大加速度、および液状化のうち少なくとも1つを表す物理量を計算することと、前記計算された物理量を前記インフラの仕様と比較することと、前記計算された物理量と前記インフラの仕様とに基づいて、潜在的なインフラの損害を示す警報を生成することとを含む。第110のさらなる実施形態によれば、前記DFOSデータを測定することは、DFOSシステムから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第109のさらなる実施形態の方法が提供される。第111のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第109のさらなる実施形態の方法が提供される。第112のさらなる実施形態によれば、第110または第111のさらなる実施形態の方法が提供される。前記光ファイバは前記インフラの下方に延在する第113のさらなる実施形態によれば、前記光ファイバは前記インフラを貫通して延在する、第110または第111のさらなる実施形態の方法が提供される。第114のさらなる実施形態によれば、前記光ファイバは、前記インフラの複数の水平方向の高さ(horizontal levels)を貫通して延在する、第113のさらなる実施形態の方法が提供される。第115のさらなる実施形態によれば、前記光ファイバは前記インフラの基礎の周囲を囲繞する、第113のさらなる実施形態の方法が提供される。第116のさらなる実施形態によれば、前記インフラは、建造物、橋、道路、記念碑、街灯、建設現場、鉱山、地下構造物、トンネル、大量輸送機関の発着所、空港、鉄道駅、工場、および海港のうちの少なくとも1つを含む、第109乃至第115のさらなる実施形態のいずれか1つの方法が提供される。第117のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第109乃至第116のさらなる実施形態のいずれか1つの方法が提供される。第118のさらなる実施形態によれば、光を光ファイバに送信すべく構成された光送信器と、光を光ファイバから受信すべく構成された受信器と、第43乃至第117のさらなる実施形態のいずれかに定義された方法を実行すべく構成されたコントローラとを含むDFOS装置が提供される。第119のさらなる実施形態によれば、第118のさらなる実施形態による1つまたは複数のDFOS装置と、1つまたは複数のDFOS装置に動作可能に接続された1つまたは複数の光ファイバとを含む、地震データを処理するシステムが提供される。第120のさらなる実施形態によれば、前記地震の震源のセンシング距離内に位置する光ファイバにおいてDFOSデータを測定することと、前記地震が発生したことを検知するために前記DFOSデータを解析することとを含む、地震の発生を検知する方法が提供される。第121のさらなる実施形態によれば、DFOS装置を前記光ファイバに接続することをさらに含み、前記DFOSデータを測定することは、前記DFOS装置から前記光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第120のさらなる実施形態の方法が提供される。第122のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから前記光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第120のさらなる実施形態の方法が提供される。第123のさらなる実施形態によれば、前記DFOSデータの解析から前記地震のマグニチュードを計算することをさらに含む、第120のさらなる実施形態の方法が提供される。第124のさらなる実施形態によれば、前記DFOSデータを測定することは、歪みを測定することを含む、第120のさらなる実施形態の方法が提供される。第125のさらなる実施形態によれば、前記DFOSデータを測定することは、地震動を測定することを含む、第120のさらなる実施形態の方法が提供される。第126のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第120乃至第125のさらなる実施形態のいずれか1つの方法が提供される。
第127のさらなる実施形態によれば、DFOSシステムによって地震の影響を解析する方法が提供され、該方法は、前記DFOSシステムでDFOSデータを捕捉することと、前記捕捉されたDFOSデータを保存することと、前記地震が発生したという外部通知を前記DFOSシステムから受信することと、前記地震の前記影響を定量化するために、前記保存されたDFOSデータを解析することとを含む。第128のさらなる実施形態によれば、前記捕捉することは時間的に連続する、第127のさらなる実施形態の方法が提供される。第129のさらなる実施形態によれば、前記保存することは、以前に保存されたDFOSデータを上書きすることを含み、新しいDFOSデータは最も古いDFOSデータを上書きする、第127乃至第128のさらなる実施形態のいずれか1つの方法が提供される。第130のさらなる実施形態によれば、前記DFOSデータを捕捉することは、光ファイバに光を放出することと、前記光ファイバから光を受信することとを含む、第127乃至第129のさらなる実施形態のいずれか1つの方法が提供される。第131のさらなる実施形態によれば、前記通知を受信することは、地震計から発信される信号を受信することを含む、第127乃至第130のさらなる実施形態のいずれか1つの方法が提供される。第132のさらなる実施形態によれば、前記通知を受信することは、無線ネットワークからの信号を受信することを含む、第127乃至第131のさらなる実施形態のいずれか1つの方法が提供される。第133のさらなる実施形態によれば、前記信号は、受振器の処理された出力である、第127乃至第132のさらなる実施形態のいずれか1つの方法が提供される。第134のさらなる実施形態によれば、前記保存されたDFOSデータを解析することは、地表面最大加速度を計算することを含む、第127乃至第133のさらなる実施形態のいずれか1つの方法が提供される。第135のさらなる実施形態によれば、前記保存されたDFOSデータを解析することは、前記光ファイバのセンシング距離内の領域における測地歪みレスポンスを測定することを含む、第127乃至第133のさらなる実施形態のいずれか1つの方法が提供される。
第136のさらなる実施形態によれば、第1の場所における地表面最大加速度を決定する方法が提供され、該方法は、前記第1の場所から第1の距離内に、第2の距離にわたって単一方向に延在する光ファイバを提供することと、DFOSデータに基づいて、前記光ファイバが水平方向に位置する地面の見かけの速度を計算することと、加速度を得るために、前記見かけの速度を時間で微分することと、を含む。第137のさらなる実施形態によれば、前記第1の距離は100メートル未満であり、前記第2の距離は100メートルより大きい、第136のさらなる実施形態の方法が提供される。
第138のさらなる実施形態によれば、地理的領域内の複数の異なる場所に対する地震の影響を視覚的に表現する方法が提供され、該方法は、前記地理的領域のセンシング距離内に少なくとも1つの光ファイバを提供することと、DFOSデータに基づいて、前記複数の異なる場所のそれぞれについて地表面最大加速度値を計算することと、前記計算された地表面最大加速度値を、前記複数の場所を表す地図上にオーバレイすることとを含む。第139のさらなる実施形態によれば、前記地表面最大加速度値を計算することは、前記少なくとも1つの光ファイバからDFOSシステムによってDFOSデータを捕捉することと、前記DFOSデータに基づいて、前記光ファイバが水平方向に位置する地面の見かけの速度を計算することと、加速度を得るために、前記見かけの速度を時間で微分することとを含む、第138のさらなる実施形態の方法が提供される。
第140のさらなる実施形態によれば、前記地震の震源のセンシング距離内に位置する光ファイバにおいてDFOSデータを測定することと、前記地震が発生したことを検知するために前記DFOSデータを解析することとを含む、地震の発生を検知する方法が提供される。第141のさらなる実施形態によれば、DFOS装置を前記光ファイバに接続することをさらに含み、前記DFOSデータを測定することは、前記DFOS装置から前記光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第140のさらなる実施形態の方法が提供される。第142のさらなる実施形態によれば、前記DFOSデータを測定することは、光通信デバイスから前記光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含む、第140のさらなる実施形態の方法が提供される。第143のさらなる実施形態によれば、前記DFOSデータの解析から前記地震のマグニチュードを計算することをさらに含む、第140のさらなる実施形態の方法が提供される。第144のさらなる実施形態によれば、前記DFOSデータを測定することは、歪みを測定することを含む、第140のさらなる実施形態の方法が提供される。第145のさらなる実施形態によれば、前記DFOSデータを測定することは、地震動を測定することを含む、第140のさらなる実施形態の方法が提供される。第146のさらなる実施形態によれば、前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、第140乃至第145のさらなる実施形態のいずれか1つの方法が提供される。
よって、本開示に従って、分散型光ファイバセンシングシステム、装置、および方法が提供されることは明らかである。本開示により、多くの代替、修正、変形が可能となる。開示された実施形態の特徴は、本発明の範囲内で組み合わせ、再配置、省略などして、追加の実施形態を生み出すことができる。さらに、特定の機能は、他の機能を付随して使用せずに有利に使用されることがある。したがって、出願人は、本発明の精神および範囲内にあるそのような代替物、修正物、等価物、および変形物をすべて包含することを意図している。
Claims (96)
- 地表面最大加速度を決定する方法であって、
少なくとも1つの光ファイバに接続されたDFOS機器を提供することと、
前記少なくとも1つの光ファイバに沿った歪みデータを含むDFOSデータを前記DFOS機器で記録することと、
前記歪みデータを複素フーリエ係数に変換することと、
前記複素フーリエ係数をスケーリングすることと、
前記スケーリングされたフーリエ係数に逆変換を適用することと、
前記少なくとも1つの光ファイバに沿った位置の前記地表面最大加速度を特定するために、前記逆変換の出力から最大値を選択することと
を含む方法。 - 前記歪みデータを複素フーリエ係数に変換することは、前記歪みデータに2次元フーリエ変換を適用することを含む、請求項1に記載の方法。
- 前記逆変換を適用することは、前記スケーリングされたフーリエ係数に2次元逆フーリエ変換を適用することを含む、請求項1に記載の方法。
- 地域地震速度モデルに基づいて地表面最大加速度を決定する方法であって、
少なくとも1つの光ファイバに接続されたDFOS機器を提供することと、
前記少なくとも1つの光ファイバに沿った歪みデータを含むDFOSデータを前記DFOS機器で記録することと、
地震の発生を検知することと、
前記光ファイバが設置されている地表面における地震エネルギーの入射角度を決定するために、地震発生地点から前記光ファイバまでの地震波経路をレイトレースすることと、
前記決定された入射角度と前記光ファイバの周辺の地面のせん断波速度に基づいて、地面の見かけの速度の測定値を計算することと、
前記歪みデータに時間微分を適用することと、
前記歪みデータに対する時間微分と、前記計算された見かけの速度の測定値とに基づいて、地表面加速度の測定値を取得することと、
前記少なくとも1つの光ファイバに沿った位置の前記地表面最大加速度を特定するために、前記地表面加速度の出力から最大値を選択することと
を含む方法。 - 前記地面の見かけの速度は水平方向における速度である、請求項4に記載の方法。
- 前記地震の発生を検知することは、DFOSデータに基づかない、請求項4に記載の方法。
- 地震によって引き起こされる地震動の特性を決定する方法であって、
少なくとも1つの光ファイバに沿ってDFOS機器でDFOSデータを記録することと、
前記DFOSデータの信号振幅の定量的測定値を得るために、前記記録されたDFOSデータを解析することと、
前記信号振幅の前記定量的測定値を所定の閾値と比較することと、
前記定量的測定値が前記所定の閾値を上回った場合に、地震が発生したと予備判断することと、
前記地震の影響を定量化することと
を含む方法。 - 前記予備判断した後、かつ前記定量化する前に地震が発生したことを検証することをさらに含む、請求項7に記載の方法。
- 前記記録されたDFOSデータを解析することは、前記地震の優位周波数の通過帯域内の信号振幅を測定することを含む、請求項7に記載の方法。
- 前記優位周波数は、局部地震では0.5~156Hzの範囲である、請求項9に記載の方法。
- 前記優位周波数は、広域地震および直下型地震では0.005~5Hzの範囲である、請求項9に記載の方法。
- 前記検証することは、非DFOS測定値に基づいて行われる、請求項8に記載の方法。
- 前記検証することは、前記記録されたDFOSデータの解析に基づいて行われる、請求項8に記載の方法。
- 前記DFOSデータは歪みデータを含み、前記定量化することは、前記少なくとも1つの光ファイバ近傍の地面の地表面最大加速度を決定することを含む、請求項7に記載の方法。
- 前記地表面最大加速度を決定することは、
前記DFOSデータの前記歪みデータを複素フーリエ係数に変換することと、
前記複素フーリエ係数をスケーリングすることと、
前記スケーリングされたフーリエ係数に逆変換を適用することと、
前記少なくとも1つの光ファイバに沿った位置の地表面最大加速度を特定するために、前記逆変換の出力から最大値を選択することと
を含む、請求項14に記載の方法。 - 前記歪みデータを複素フーリエ係数に変換することは、前記歪みデータに2次元フーリエ変換を適用することを含む、請求項15に記載の方法。
- 前記逆変換を適用することは、前記スケーリングされたフーリエ係数に2次元逆フーリエ変換を適用することを含む、請求項15に記載の方法。
- 地震によって引き起こされる地震動の特性を決定する方法であって、
少なくとも1つの光ファイバに沿ってDFOS機器でDFOSデータを記録することと、
前記DFOSデータの信号品質を測定することと、
前記測定されたDFOSデータの信号品質を所定の閾値と比較することと、
前記測定された信号品質が所定の閾値未満である場合に、地震が発生したと予備判断することと、
前記地震の影響を定量化することと
を含む方法。 - 前記予備判断した後、かつ前記定量化する前に地震が発生したことを検証することをさらに含む、請求項18に記載の方法。
- 前記信号品質を測定することは、前記DFOSデータの信号対雑音比を決定することを含む、請求項18に記載の方法。
- 前記信号品質を測定することは、前記記録されたDFOSデータに適用される、請求項18に記載の方法。
- 前記信号品質を測定することは、前記DFOSデータを記録することと同時に前記DFOSデータに適用される、請求項18に記載の方法。
- 前記DFOSデータは歪みデータを含み、前記定量化することは、前記少なくとも1つの光ファイバ近傍の地面の地表面最大加速度を決定することを含む、請求項18に記載の方法。
- 前記地表面最大加速度を決定することは、
前記DFOSデータの前記歪みデータを複素フーリエ係数に変換することと、
前記複素フーリエ係数をスケーリングすることと、
前記スケーリングされたフーリエ係数に逆変換を適用することと、
前記少なくとも1つの光ファイバに沿った位置の地表面最大加速度を特定するために、前記逆変換の出力から最大値を選択することと
含む、請求項23に記載の方法。 - 前記歪みデータを複素フーリエ係数に変換することは、前記歪みデータに2次元フーリエ変換を適用することを含む、請求項24に記載の方法。
- 前記逆変換を適用することは、前記スケーリングされたフーリエ係数に2次元逆フーリエ変換を適用することを含む、請求項25に記載の方法。
- 構造物に対する地震動の影響を評価する方法であって、
少なくとも1つの光ファイバに沿った歪みデータを含むDFOSデータをDFOS機器で記録することと、
前記DFOSデータに基づいて前記構造物の共振ピークを計算することと、
前記計算された共振ピークを、前記構造物に関連付けられた共振ピーク値をアーカイブする記憶装置に保存することと、
前記記憶装置からアーカイブされた共振ピーク値を取り出すことと、
前記取り出されたアーカイブ値と前記DFOSデータに基づいて計算された共振ピークとの間の、変化値によって表される変化を検知することと、
前記変化値を閾値と比較することと、
前記検知された変化が前記閾値を上回った場合に警報を生成することと
を含む方法。 - 前記少なくとも1つの光ファイバを前記構造物の近傍に設置することをさらに含む、請求項27に記載の方法。
- 前記少なくとも1つの光ファイバは、前記構造物の内部または下にある、請求項28に記載の方法。
- 前記方法はスケジュールの通りに繰り返される、請求項27に記載の方法。
- 前記スケジュールの期間は5分から48時間の間である、請求項30に記載の方法。
- 前記スケジュールの期間は10分、60分、12時間、または24時間である、請求項31に記載の方法。
- 前記スケジュールの期間は適応的であり、前記比較の結果に応じて継続期間が減少する、請求項30に記載の方法。
- 前記構造物は、建造物、橋、塔のいずれかである、請求項27に記載の方法。
- 前記構造物の設計に基づく工学設計パラメータが記憶装置に保存される、請求項27に記載の方法。
- 前記記憶装置は、前記構造物に関連付けられた複数の共振ピーク値を記憶し、前記共振ピーク値の各々は、タイムスタンプおよび継続時間のうちの少なくとも1つに関連付けられる、請求項35に記載の方法。
- 最後に計算された共振ピークと前記記憶装置に保存された最も古い共振ピーク値との間の変化を表す累積変化値を特定するために、前記検知された変化値を経時的に累積することをさらに含む、請求項36に記載の方法。
- 前記地震動は地震によって引き起こされる、請求項27乃至37のいずれか一項に記載の方法。
- 前記閾値は適応的であり、経時的に低下する、請求項27に記載の方法。
- 前記警報を生成することは、有線または無線ネットワークを通じて信号を送信することを含む、請求項27に記載の方法。
- 前記警報を生成することは、前記構造物の内部で人間が聞き取れるメッセージを出力することを含む、請求項27に記載の方法。
- 前記人間が聞き取れるメッセージは、前記構造物に対する潜在的損害に関する情報を含む、請求項41に記載の方法。
- 第1の場所を震源とする地震の影響を定量化する方法であって、
前記地震の間、前記第1の場所からずれた第2の場所でDFOSシステムを用いてDFOSデータを測定することと、
前記測定されたDFOSデータを用いて、地表面最大加速度、地表面最大速度、最大地動変位、最大地動歪み(peak ground strain)、最大地動歪み率、最大スペクトル加速度、および/または地震強度の少なくとも1つを計算することと
を含む方法。 - 前記測定されたDFOSデータを記録することと、
地震が発生したという通知を受信することと、
前記受信された通知に応答して、前記記録されている測定されたDFOSデータを使用して、発生した地震の地震マグニチュード、モーメント、歪み、歪み率、測地変形、位置、深さ、焦点メカニズム、放射パターン、モーメントテンソル、滑りの有限断層領域、滑り速度、および破壊速度のうちの少なくとも1つを計算することと
をさらに含む、請求項43に記載の方法。 - 前記計算することに基づいて、前記DFOSシステムから警報を生成することと、
前記生成された警報を人口密集地域のレシーバに送信することと
をさらに含む、請求項44に記載の方法。 - 前記DFOSデータを測定することは、DFOS地震センシング用および電気通信トラフィックの伝送用に同時に使用される単一の電気通信ファイバ上でDFOS信号を送受信することを含む、請求項43に記載の方法。
- 前記DFOSデータを測定することは前記地震の発生中に行われ、前記計算することは前記地震の発生後に行われる、請求項43に記載の方法。
- 前記DFOSデータを測定することおよび前記計算することは、前記地震の発生中に行われる、請求項43に記載の方法。
- 前記計算することは、前記第2の場所における物理的影響を表す、請求項43に記載の方法。
- DFOS装置は前記第2の場所に配置される、請求項43に記載の方法。
- 前記計算することは、前記第1の場所と前記第2の場所との間に位置する第3の場所における物理的影響を表す、請求項43に記載の方法。
- 前記DFOSデータを測定することは、前記DFOSシステムから光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含み、前記計算することは、前記光ファイバから所定の距離内に位置する第3の場所における物理的影響を表す、請求項43に記載の方法。
- 前記所定の距離は50メートルから500メートルの間である、請求項52に記載の方法。
- 前記所定の距離は100メートルである、請求項53に記載の方法。
- 前記第2の場所は、建造物、データセンター、病院、空港、重要インフラ、道路、橋、トンネル、住居、ホテル、他の構造物、または地震動の影響を受ける建築物の少なくとも1つの敷地である、請求項49に記載の方法。
- 前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、請求項43乃至55のいずれか一項に記載の方法。
- 地震の影響を定量化する方法であって、
少なくとも前記地震の発生中にDFOSデータを測定することと、
前記測定されたDFOSデータを用いて、前記地震の発生中の地震動の方向および前記地震の発生中の前記地震動の継続時間の少なくとも1つを計算することと
を含む方法。 - 前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、請求項57に記載の方法。
- 建造物に対する地震の影響を解析する方法であって、
少なくとも前記地震の発生中にDFOSデータを測定することと、
前記測定されたDFOSデータから前記建造物の固有モードの共振値を計算することと、
前記固有モードから建造物の振動の強さを表す量を決定することと
を含む方法。 - 前記DFOSデータを測定することは、DFOS装置から光ファイバに光を送信し、前記光ファイバから屈折した光を受信することを含み、前記光ファイバは、前記建造物から所定距離離れて延在するか、または前記建造物を通過する、請求項59に記載の方法。
- 前記DFOSデータを測定することは、前記地震によって生じる振動中に行われる、請求項59に記載の方法。
- 前記計算することは、前記地震によって生じる振動中に行われる、請求項59に記載の方法。
- 前記計算することは、前記地震により生じる振動が終了した後に行われる、請求項61に記載の方法。
- 前記決定することは、前記地震の発生後に計算された前記建造物の固有モードを前記地震前の固有モードと比較することを含む、請求項59に記載の方法。
- 前記固有モードの変化を前記建造物の構造的健全性の状態変化に関連付けることをさらに含む、請求項64に記載の方法。
- 前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、請求項59乃至65のいずれか一項に記載の方法。
- 地震の余波において構造物の安全警報を生成する方法であって、
少なくとも地震の発生中にDFOSデータを測定することと、
前記測定されたDFOSデータに基づいて前記構造物の安全スコアを計算することと、
前記安全スコアに基づいて安全警報を生成することと
を含む方法。 - 前記安全スコアは、前記測定されたDFOSデータから計算された地表面最大加速度と、前記構造物の地表面最大加速度率(peak ground acceleration rating)との比較に基づく、請求項67に記載の方法。
- 前記安全警報は、前記計算された安全スコアの表現を含む、請求項67に記載の方法。
- 前記安全警報は複数の重大度レベルに分類され、前記生成することは前記重大度レベルに基づく、請求項67に記載の方法。
- 前記生成された安全警報を、重大度レベルに基づいて選択された受信者に送信することをさらに含む、請求項70に記載の方法。
- 前記安全警報は、前記構造物が人間の居住者が立ち入ることについて安全であるか否かを示す、請求項67に記載の方法。
- 光ファイバの物理的変化の検知に応答して、前記測定されたDFOSデータを記憶装置に保存することをさらに含む、請求項67に記載の方法。
- 前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、請求項67乃至73のいずれか一項に記載の方法。
- インフラの損害を検知する方法であって、
前記インフラの内部もしくは近傍の光ファイバでDFOSデータを測定することと、
前記測定されたDFOSに基づいて、建造物の層間変位、前記建造物の下部もしくは側部の地表面最大加速度、前記インフラ自体の構成要素の地表面最大加速度、および液状化のうち少なくとも1つを表す物理量を計算することと、
前記計算された物理量を前記インフラの仕様と比較することと、
前記計算された物理量と前記インフラの仕様とに基づいて、潜在的なインフラの損害を示す警報を生成することと
を含む方法。 - 前記光ファイバは前記インフラの下方に延在する、請求項75に記載の方法。
- 前記光ファイバは前記インフラを貫通して延在する、請求項75に記載の方法。
- 前記光ファイバは、前記インフラの複数の水平方向の高さ(horizontal levels)を貫通して延在する、請求項77に記載の方法。
- 前記光ファイバは前記インフラの基礎の周囲を囲繞する、請求項77に記載の方法。
- 前記インフラは、建造物、橋、道路、記念碑、街灯、建設現場、鉱山、地下構造物、トンネル、大量輸送機関の発着所、空港、鉄道駅、工場、および海港のうちの少なくとも1つを含む、請求項75乃至79のいずれか一項に記載の方法。
- 前記DFOSデータを測定することは、そうしなければ点灯している光ファイバから未点灯または未使用の光スペクトルの一部を選択すること、またはそうしなければ点灯している光ファイバから偏波固有モードを選択すること、または確立された専用の電気通信ネットワークからダークファイバ全体を選択することを含む、請求項75乃至80のいずれか一項に記載の方法。
- 光を光ファイバに送信すべく構成された光送信器と、
光を光ファイバから受信すべく構成された受信器と、
請求項1乃至81のいずれか一項に定義された方法を実行すべく構成されたコントローラと
を含む、DFOS装置。 - 地震データを処理するシステムであって、
1つまたは複数の請求項82に記載のDFOS装置と、
1つまたは複数のDFOS装置に動作可能に接続された1つまたは複数の光ファイバと
を含むシステム。 - DFOSシステムによって地震の影響を解析する方法であって、
前記DFOSシステムでDFOSデータを捕捉することと、
前記捕捉されたDFOSデータを保存することと、
前記地震が発生したという外部の通知を前記DFOSシステムから受信することと、
前記地震の前記影響を定量化するために、前記保存されたDFOSデータを解析することと
を含む方法。 - 前記捕捉することは時間的に連続する、請求項84に記載の方法。
- 前記保存することは、以前に保存されたDFOSデータを上書きすることを含み、新しいDFOSデータは最も古いDFOSデータを上書きする、請求項84乃至85のいずれか一項に記載の方法。
- 前記DFOSデータを捕捉することは、光ファイバに光を放出することと、前記光ファイバから光を受信することとを含む、請求項84乃至86のいずれか一項に記載の方法。
- 前記通知を受信することは、地震計から発信される信号を受信することを含む、請求項84乃至87のいずれか一項に記載の方法。
- 前記通知を受信することは、無線ネットワークからの信号を受信することを含む、請求項84乃至88のいずれか一項に記載の方法。
- 前記信号は、受振器の処理された出力である、請求項84乃至89のいずれか一項に記載の方法。
- 前記保存されたDFOSデータを解析することは、地表面最大加速度を計算することを含む、請求項84乃至90のいずれか一項に記載の方法。
- 前記保存されたDFOSデータを解析することは、前記光ファイバのセンシング距離内の領域における測地歪みレスポンスを測定することを含む、請求項84乃至90のいずれか一項に記載の方法。
- 第1の場所における地表面最大加速度を決定する方法であって、
前記第1の場所から第1の距離内に、第2の距離にわたって単一方向に延在する光ファイバを提供することと、
DFOSデータに基づいて、前記光ファイバが水平方向に位置する地面の見かけの速度を計算することと、
加速度を得るために、前記見かけの速度を時間で微分することと
を含む方法。 - 前記第1の距離は100メートル未満であり、前記第2の距離は100メートルより大きい、請求項93に記載の方法。
- 地理的領域内の複数の異なる場所に対する地震の影響を視覚的に表現する方法であって、
前記地理的領域のセンシング距離内に少なくとも1つの光ファイバを提供することと、
DFOSデータに基づいて、前記複数の異なる場所のそれぞれについて地表面最大加速度値を計算することと、
前記計算された地表面最大加速度値を、前記複数の場所を表す地図上にオーバレイすることと
を含む方法。 - 前記地表面最大加速度値を計算することは、前記少なくとも1つの光ファイバからDFOSシステムによってDFOSデータを捕捉することと、
前記DFOSデータに基づいて、前記光ファイバが水平方向に位置する地面の見かけの速度を計算することと、
加速度を得るために、前記見かけの速度を時間で微分することと
を含む、請求項95に記載の方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US202163157789P | 2021-03-07 | 2021-03-07 | |
US63/157,789 | 2021-03-07 | ||
PCT/US2022/019113 WO2022192125A2 (en) | 2021-03-07 | 2022-03-07 | Devices, systems, and methods providing earthquake information with distributed fiber-optic sensing |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2024510957A true JP2024510957A (ja) | 2024-03-12 |
Family
ID=83228501
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2023554858A Pending JP2024510957A (ja) | 2021-03-07 | 2022-03-07 | 分散型光ファイバセンシングにより地震情報を提供する装置、システム、および方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20240142646A1 (ja) |
JP (1) | JP2024510957A (ja) |
WO (1) | WO2022192125A2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115856996B (zh) * | 2022-11-08 | 2023-06-30 | 应急管理部国家自然灾害防治研究院 | 一种地震数值预测预报方法和系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5627637A (en) * | 1995-02-24 | 1997-05-06 | Kapteyn; Kelvin L. | Fully distributed optical fiber strain sensor |
US5850622A (en) * | 1996-11-08 | 1998-12-15 | Amoco Corporation | Time-frequency processing and analysis of seismic data using very short-time fourier transforms |
US7916303B2 (en) * | 2007-11-13 | 2011-03-29 | Optoplan As | Non-uniform sampling to extend dynamic range of interferometric sensors |
US8619497B1 (en) * | 2012-11-15 | 2013-12-31 | Cggveritas Services Sa | Device and method for continuous data acquisition |
US11204434B2 (en) * | 2015-12-16 | 2021-12-21 | Halliburton Energy Services, Inc. | Large area seismic monitoring using fiber optic sensing |
KR101946855B1 (ko) * | 2018-09-07 | 2019-02-15 | 대한민국 | 지진동 감지 센서 및 이를 이용한 지진 대응 시스템 |
-
2022
- 2022-03-07 WO PCT/US2022/019113 patent/WO2022192125A2/en active Application Filing
- 2022-03-07 JP JP2023554858A patent/JP2024510957A/ja active Pending
- 2022-03-07 US US18/280,652 patent/US20240142646A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
WO2022192125A3 (en) | 2022-11-17 |
US20240142646A1 (en) | 2024-05-02 |
WO2022192125A2 (en) | 2022-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9435902B2 (en) | Wide area seismic detection | |
US8131121B2 (en) | Optical fiber pipeline monitoring system and method | |
JP4324126B2 (ja) | 地中内観測システムおよび地中内観測方法 | |
Picozzi et al. | Interferometric analysis of strong ground motion for structural health monitoring: the example of the L’Aquila, Italy, seismic sequence of 2009 | |
GB2467419A (en) | Monitoring soil slope displacement rate by detecting acoustic emissions | |
JP3894494B2 (ja) | 土砂災害予知システム、地域情報提供システム及び土砂災害予知方法 | |
US20110265547A1 (en) | Sensor, methods of calibrating a sensor, methods of operating a sensor | |
Kaya et al. | British Columbia smart infrastructure monitoring system | |
JP2024510957A (ja) | 分散型光ファイバセンシングにより地震情報を提供する装置、システム、および方法 | |
McKenna et al. | Remote structural infrasound: Case studies of real-time infrastructure system monitoring | |
Nayak et al. | Addressing the risk of the tsunami in the Indian Ocean | |
Kapogianni et al. | Structural health monitoring of the Athenian Acropolis' Walls via optical fibre sensors, accelerographs and numerical simulations | |
Lin et al. | Time series analysis of acoustic emissions in prefabricated utility tunnel during the sealing test | |
EP2965124B1 (en) | A system and method for evaluating a dam state | |
Fujihashi et al. | Development of seafloor seismic and tsunami observation system | |
Hoult et al. | The 2012 Moe earthquake and earthquake attenuation in south eastern Australia | |
Moya et al. | Integration of monitoring and inspection systems for geohazard assessment on pipelines that cross Amazonian Jungles and the Andes | |
CN109375266B (zh) | 一种采用斜长分布式光纤的地下水封洞库安全监测系统 | |
RU2809469C1 (ru) | Способ и система сейсмоакустического контроля массива горных пород | |
Anjana et al. | Optical Fibre Sensing and Deep Learning-Based Disaster and Climate Change Risk Assessments of Civil Infrastructure: Current Status and Future Perspective | |
Smith et al. | Listening for landslides: method, measurements and the Peace River case study | |
Biller | Detection of seepage-induced internal instability using acoustic emission | |
RU2377602C1 (ru) | Способ определения параметров сейсмического удара | |
Brigante et al. | Vibration-Based Procedure for the Structural Assessment of Heritage Structures | |
D'Angio | Rheological effects related to neo-fracturing processes in rock masses |