変形計測システムおよび変形計測方法
【課題】導電率の影響を考慮して２以上の三軸センサの相対位置を計測することにより、盛立構造物の変形を高精度に計測することが可能な変形計測システム等を提供することを目的とする。
【解決手段】本発明にかかる変形計測システムの代表的な構成は、盛立構造物１０２に埋設される移動用三軸センサ１１０と、盛立構造物の表面近傍に配置される基準用三軸センサ１１２と、較正用三軸センサ１１４と、三軸センサから出力された信号を受信する信号受信部１１６と、導電率を入力する入力部１１８と、第１の磁場解析を行い磁場強度Ｈ１を導出する較正用磁場導出部１２２と、換算比率を導出する換算比率導出部１２４と、三軸センサが検出した誘起電圧Ｅ２と前記換算比率とから磁場強度Ｈ２を求める磁場強度換算部１２６と、第２の磁場解析を行い前記誘起電圧Ｅ２を誘起させる磁場を励磁させた三軸センサの位置を同定する位置同定部１２８とを備えることを特徴とする。
【解決手段】本発明にかかる変形計測システムの代表的な構成は、盛立構造物１０２に埋設される移動用三軸センサ１１０と、盛立構造物の表面近傍に配置される基準用三軸センサ１１２と、較正用三軸センサ１１４と、三軸センサから出力された信号を受信する信号受信部１１６と、導電率を入力する入力部１１８と、第１の磁場解析を行い磁場強度Ｈ１を導出する較正用磁場導出部１２２と、換算比率を導出する換算比率導出部１２４と、三軸センサが検出した誘起電圧Ｅ２と前記換算比率とから磁場強度Ｈ２を求める磁場強度換算部１２６と、第２の磁場解析を行い前記誘起電圧Ｅ２を誘起させる磁場を励磁させた三軸センサの位置を同定する位置同定部１２８とを備えることを特徴とする。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、２以上の三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測システムおよび変形計測方法に関するものである。
【背景技術】
【０００２】
フィルダムや堤防の堤体に代表される盛立構造物においては、安全上の理由から（崩壊の危険性を把握するために）施工中および完成後その変形を監視する必要がある。そのため、従来、基礎岩盤に固定された支柱やワイヤ等によって連結された沈下計などの埋設型センサを埋設して、それらの相対位置を計測することによって盛立構造物の変形を計測していた。
【０００３】
しかし、支柱やワイヤは盛立構造物にとっては異物であり、強度等の性質低下を招きかねない。また、かかる手法では、支柱やワイヤと埋設型センサとの相対変位、すなわち１次元の変位しか計測することができない。その上、従来の埋設型センサは、信号を入出力するためにケーブルを要するので、工事中におけるケーブルの養生やメンテナンスの労力といった課題を生じることとなる。さらに、これらのケーブルが水路（みずみち：水の通り路）を形成するおそれもある。
【０００４】
そこで、このような課題を解決する技術が求められている。これに関連する従来技術としては、特許文献１に、三軸センサを２点に配置して一方の三軸センサを励磁した場合に他方の三軸センサに誘起される電圧（Ｅ１）の真空中における距離減衰特性を利用して、２点間の相対的な距離（ｒ１）を以下の式（式１）より算出する技術が開示されている。かかる技術は、地中や水中であっても適用可能とされている。
ｒ１＝｛（Ｅ０／Ｅ１）１／３｝・ｒ０ （式１）
ここで、ｒ０：既知の距離
Ｅ０：２点間の距離がｒ０の場合に三軸センサに誘起される電圧
ｒ１：２点間の距離（算出対象）
Ｅ１：算出対象時において三軸センサに誘起される電圧
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００３−４４０９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１の技術は、真空中すなわち導電率の影響を排除できる環境を対象としている。現実的には、盛立構造物の地中においては、その構成要素である岩石、砂利、水等の導電率（例えばロックフィルダム堤体の導電率は０．００１から０．１Ｓ/ｍ）により磁場の分布に歪みを生じるので、当然ながらその影響を考慮する必要がある。
【０００７】
よって、導電率の影響を排除し得ない盛立構造物に対し、特許文献１の技術を適用した場合には、測定に無視できない誤差を生じることとなる。すなわち、この技術は精度の点に課題がある。その上、特にフィルダムでは中央から外側に行くにつれて土質（岩石、砂利、土等）を異ならせるため、箇所によって導電率が異なり、また三軸センサの送信側と受信側の間に土質の境界が存在する場合もあるため、既定のデータ等から一律に導電率を定めることに対して難がある。
【０００８】
そこで、本発明は、盛立構造物の導電率の影響を考慮して２以上の三軸センサの相対位置を計測することにより、盛立構造物の変形を高精度に計測することが可能な変形計測システムおよび変形計測方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するために本発明の代表的な構成は、軸方向が異なる３つのコイルを備え、コイルを励磁することによって磁場を生じ、またはコイルが磁場を受けて発生した誘起電圧を検出して信号として出力する三軸センサを用いて、２以上の三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測システムであって、盛立構造物に埋設される移動用三軸センサと、盛立構造物の表面近傍に配置される基準用三軸センサと、移動用三軸センサから所定の基準距離に配置される較正用三軸センサと、移動用三軸センサまたは基準用三軸センサ、もしくは較正用三軸センサから出力された信号を受信する信号受信部と、盛立構造物の導電率を入力する入力部と、移動用三軸センサから出力された磁場について基準距離および導電率を用いて第１の磁場解析を行い、較正用三軸センサの位置における磁場強度Ｈ１を導出する較正用磁場導出部と、較正用三軸センサの位置における磁場強度Ｈ１とこの較正用三軸センサが検出した誘起電圧Ｅ１とから、盛立構造物において磁場強度と誘起電圧とを換算する換算比率を導出する換算比率導出部と、移動用三軸センサまたは基準用三軸センサの一方が検出した誘起電圧Ｅ２と換算比率とから、誘起電圧Ｅ２を検出した三軸センサの位置における磁場強度Ｈ２を求める磁場強度換算部と、導電率を反映した第２の磁場解析を行い、移動用三軸センサまたは基準用三軸センサのうち、誘起電圧Ｅ２を検出した一方の三軸センサの位置における磁場強度が磁場強度換算部が求めた磁場強度Ｈ２となるように、他方の三軸センサの現在位置を同定する位置同定部と、を備えることを特徴とする。
【００１０】
かかる構成によれば、盛立構造物の導電率の影響を考慮して２以上の三軸センサの相対位置を計測することが可能となり、盛立構造物の変形を高精度に計測することができる。特に、磁場解析による磁場強度と較正用に実測した誘起電圧によって求めた換算比率を用いることにより、その盛立構造物における磁場の強度と誘起電圧の関係を求めることができ、埋設された三軸センサの位置を同定することができる。
【００１１】
上記導電率は、盛立構造物を測定した実測値であるとよい。上述したように、盛立構造物は、地点によって構成要素（岩石、砂利、水等）が異なるため、既定のデータ等から一律に導電率を定めることが困難である。一方、盛立構造物の実測によれば、現場に則した値を取得できるため最適な導電率を設定することができる。
【００１２】
上記移動用三軸センサは発信機、基準用三軸センサは受信機であって、位置同定部は、移動用三軸センサが埋設された初期位置の周囲に複数の仮想発信点を設定し、それぞれの仮想発信点に移動用三軸センサがあると仮定した場合に、この移動用三軸センサから出力された磁場について第２の磁場解析を行うことにより磁場強度換算部が求めた磁場強度Ｈ２となる複数の仮想受信点を求め、複数の仮想受信点と基準用三軸センサの位置関係と、複数の仮想発信点の設定位置とから、比例配分により移動用三軸センサの現在位置を同定するとよい。これにより、簡潔かつ容易に移動用三軸センサの現在位置を同定することができる。
【００１３】
上記移動用三軸センサは発信機、基準用三軸センサは受信機であって、位置同定部は、移動用三軸センサが埋設された初期位置の周囲に複数の仮想発信点を設定し、それぞれの仮想発信点に移動用三軸センサがあると仮定した場合に、この移動用三軸センサから出力された磁場について第２の磁場解析を行うことにより基準用三軸センサの位置における複数の磁場強度を求め、複数の磁場強度および磁場強度換算部が求めた磁場強度Ｈ２との大小関係と、複数の仮想発信点の設定位置とから、比例配分により移動用三軸センサの現在位置を同定してもよい。これにより、簡潔かつ容易に移動用三軸センサの現在位置を同定することができる。
【００１４】
上記位置同定部は、空間をメッシュ状に分割する３次元解析によって第２の磁場解析を行い、メッシュの１辺の長さは、移動用三軸センサが備えるコイルの直径の１／４以上１／２以下に設定するとともに、初期位置から少なくとも２０ｍの範囲内の空間は他の範囲よりも細かく分割するとよい。これにより、高精度に移動用三軸センサの現在位置の同定を行うことができる。
【００１５】
上記盛立構造物は、土砂や岩石を盛り立てて構成されるフィルダムの堤体であるとよい。すなわち、当該変形計測システムは、フィルダムを主な適用対象とする。
【００１６】
上記課題を解決するために本発明の他の代表的な構成は、軸方向が異なる３つのコイルを備え、コイルを励磁することによって磁場を生じ、またはコイルが磁場を受けて発生した誘起電圧を検出して信号として出力する三軸センサを用いて、２以上の三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測方法であって、盛立構造物の初期位置に移動用三軸センサを埋設し、盛立構造物の導電率を測定し、移動用三軸センサから所定の基準距離の位置に較正用三軸センサを配置し、較正用三軸センサを用いて誘起電圧Ｅ１を測定し、移動用三軸センサから出力される磁場について基準距離および導電率を用いて第１の磁場解析を行い、較正用三軸センサの位置における磁場強度Ｈ１を導出し、較正用三軸センサの位置における磁場強度Ｈ１とこの較正用三軸センサが検出した誘起電圧Ｅ１とから、盛立構造物において磁場強度と誘起電圧とを換算する換算比率を導出し、盛立構造物の表面近傍に基準用三軸センサを配置し、移動用三軸センサまたは基準用三軸センサの一方によって誘起電圧Ｅ２を測定し、誘起電圧Ｅ２と換算比率とから、誘起電圧Ｅ２を検出した三軸センサの位置における磁場強度Ｈ２を導出し、導電率を反映した第２の磁場解析を行い、移動用三軸センサまたは基準用三軸センサのうち、誘起電圧Ｅ２を検出した一方の三軸センサの位置における磁場強度が磁場強度換算部が求めた磁場強度Ｈ２となるように、他方の三軸センサの位置を同定するとよい。
【００１７】
上述した変形計測システムにおける技術的思想に対応する構成要素やその説明は、当該変形計測方法にも適用可能である。
【発明の効果】
【００１８】
本発明によれば、盛立構造物の導電率の影響を考慮して２以上の三軸センサの相対位置を計測することにより、盛立構造物の変形を高精度に計測することが可能な変形計測システムおよび変形計測方法を提供することができる。
【図面の簡単な説明】
【００１９】
【図１】三軸センサについて説明する図である。
【図２】本実施形態にかかる変形計測システムの構成を示すブロック図である。
【図３】変形計測システムが盛立構造物の変形を計測する一連の動作を説明するフローチャートである。
【図４】盛立構造物の盛立途中および盛立完了後を示す図である。
【図５】盛立構造物の導電率の実測手法の一例を示す図である。
【図６】移動用三軸センサの現在位置を同定する手法について説明する図である。
【図７】移動用三軸センサの現在位置を同定する手法について説明する図である。
【図８】実施例および比較例の解析モデルを示す図である。
【図９】実施例および比較例における発信点からの距離と磁場強度の解析誤差の相関を示す図である。
【発明を実施するための形態】
【００２０】
以下に添付図面を参照しながら、本発明の好適な実施形態、実施例について詳細に説明する。かかる実施形態、実施例に示す寸法、材料、その他具体的な数値などは、発明の理解を容易とするための例示に過ぎず、特に断る場合を除き、本発明を限定するものではない。なお、本明細書及び図面において、実質的に同一の機能、構成を有する要素については、同一の符号を付することにより重複説明を省略し、また本発明に直接関係のない要素は図示を省略する。
【００２１】
（三軸センサ）
図１は、三軸センサ１０８ａ、１０８ｂについて説明する図である。特に、図１（ａ）は発信機である三軸センサ１０８ａについて説明する図、図１（ｂ）は受信機である三軸センサ１０８ｂについて説明する図である。以下、三軸センサ１０８ａ、１０８ｂについて説明するが、その詳細については特許文献１（特開２００３−４４０９号公報）に開示されているため以下簡略化して説明する。
【００２２】
図１（ａ）に示すように、三軸センサ１０８ａ（発信機）は、発信回路２１２、増幅器２１４、励磁コイル２１６、切替回路２１８を包含する。
【００２３】
発信回路２１２は、所定の周波数の交流信号を発生させて、増幅器２１４へと出力する。このときの周波数は、高くなると盛立構造物に吸収されやすくなってしまい、低すぎると検波しにくくなるため、可聴帯域（２０Ｈｚ〜１５ｋＨｚ程度）を用いることができ、さらに好適には１０ｋＨｚ程度を用いることができる。増幅器２１４は、発信回路２１２から出力された交流信号の電流量を増幅して、励磁コイル２１６へと出力する。
【００２４】
励磁コイル２１６は、軸方向が互いに直交する三軸座標系のそれぞれの軸に配置された３つコイル（ｘ,ｙ,ｚ）を備えており、それぞれのコイルは増幅器２１４から与えられた交流信号により励磁され磁場を生じる。切替回路２１８は、交流信号の出力対象とされるコイル（ｘ,ｙ,ｚ）を順次切り替える。
【００２５】
図１（ｂ）に示すように、三軸センサ１０８ｂ（受信機）は、検出コイル２２２、切替回路２２４、検波回路２２６、電圧計２２８、誘起電圧演算部２３０、信号発信部２３２を包含する。
【００２６】
検出コイル２２２は、軸方向が互いに直交する三軸座標系のそれぞれの軸に配置された３つのコイル（ｕ,ｖ,ｗ）であって、励磁コイル２１６（ｘ,ｙ,ｚ）が発生する磁場に伴い誘起電圧を発生する。それぞれのコイル（ｕ,ｖ,ｗ）に発生した誘起電圧は、検波回路２２６を介して電圧計２２８により測定される。切替回路２２４は、電圧計２２８の測定対象とされるコイル（ｕ,ｖ,ｗ）を順次切り替える。
【００２７】
よって、発信機の励磁コイル２１６（ｘ,ｙ,ｚ）をそれぞれ励磁した場合、受信機ではコイルｘの励磁に起因するそれぞれのコイル（ｕ,ｖ,ｗ）の誘起電圧Ｅxu、Ｅxv、Ｅxw、コイルｙの励磁に起因するそれぞれのコイル（ｕ,ｖ,ｗ）の誘起電圧Ｅyu、Ｅyv、Ｅyw、コイルｚの励磁に起因するそれぞれのコイル（ｕ,ｖ,ｗ）の誘起電圧Ｅzu、Ｅzv、Ｅzwの合計９種類の電圧が測定される。
【００２８】
誘起電圧演算部２３０は、これら９つのパラメータ（誘起電圧）を用いて、受信機が配置される位置における誘起電圧のｘ成分Ｅx、ｙ成分Ｅy、ｚ成分Ｅzを以下の式（式２から式４）に基づき算出する。
Ｅx＝（Ｅxu＋Ｅxv＋Ｅxw）1/2 （式２）
Ｅｙ＝（Ｅyu＋Ｅyv＋Ｅyw）1/2 （式３）
Ｅｚ＝（Ｅzu＋Ｅzv＋Ｅzw）1/2 （式４）
【００２９】
信号発信部２３２は、誘起電圧演算部２３０が算出した誘起電圧を、後述する信号受信部１１６へと発信する。なお、発信側の三軸センサ１０８ａにも不図示の信号受信部を備えて、算出した誘起電圧を発信側の三軸センサ１０８ａを介して信号受信部１１６へと発信してもよい。本実施形態では、この誘起電圧Ｅ（Ｅx、Ｅy、Ｅz）を用いて、発信機が配置される位置と、受信機が配置される位置の間の距離（相対位置）を磁場解析により求める。
【００３０】
なお上記説明では三軸センサ１０８ａ、１０８ｂの構成を送信機と受信機として説明したが、コイルを共用として両方の構成を備えることにより、送受信のいずれも可能な三軸センサとすることができる。
【００３１】
（変形計測システムおよび変形計測方法）
図２は、本実施形態にかかる変形計測システム１００の構成を示すブロック図である。図２に示すように、変形計測システム１００は、移動用三軸センサ１１０、基準用三軸センサ１１２、較正用三軸センサ１１４、信号受信部１１６、入力部１１８、記憶部１２０、較正用磁場導出部１２２、換算比率導出部１２４、磁場強度換算部１２６、位置同定部１２８、変形計測部１３０、報知部１３２を包含する。移動用三軸センサ１１０、基準用三軸センサ１１２、較正用三軸センサ１１４は、いずれも装置としては図１を用いて説明した三軸センサ１０８ａ（発信機）または三軸センサ１０８ｂ（受信機）であるが、目的が異なるために、説明の便宜上名称を変えて示している。本実施形態では、移動用三軸センサ１１０を発信機、基準用三軸センサ１１２および較正用三軸センサ１１４を受信機とする。以下、各部の構成をその動作と併せて説明する。
【００３２】
図３は、変形計測システム１００が盛立構造物１０２の変形を計測する一連の動作を説明するフローチャートである。また、図４は、盛立構造物１０２の盛立途中および盛立完了後を示す図である。特に、図４（ａ）は盛立構造物１０２の盛立途中を示す図、図４（ｂ）は盛立構造物１０２の盛立完了後を示す図である。なお、図４では、理解を容易にするために移動用三軸センサ１１０、基準用三軸センサ１１２を１つとして図示しているが、実際にはこれらは複数設けられ得る。
【００３３】
変形計測システム１００の導入に当たっては、まず変形を計測する盛立構造物１０２が所定の高さまで盛り立てられる（Ｓ１５０）。本実施形態において、盛立構造物１０２とは、段階的に盛立てられて構築される構造物のことであり、フィルダムや堤防の堤体、盛土等が例示される。さらに、掘削された後に段階的に埋め戻される廃棄物処分場なども包含される。
【００３４】
次に、図４（ａ）に示すように、移動用三軸センサ１１０が初期位置に配置される（Ｓ１５２）。移動用三軸センサ１１０は、盛立構造物１０２に埋設され、完成後にその位置（変位）を検知することによって盛立構造物１０２の変形を検知するためのセンサである。
【００３５】
次に、移動用三軸センサ１１０から所定の基準距離Ｌに較正用三軸センサ１１４が配置される（Ｓ１５４）。較正用三軸センサ１１４は、次に述べる換算比率を求めるために較正（キャリブレーション）を行うためのセンサである。
【００３６】
基準距離Ｌは任意に設定することができるが、盛立完成時における移動用三軸センサ１１０と基準用三軸センサ１１２の距離（予定された距離）とすることが好ましい。このように、運用時の位置関係に近い条件で較正を行うことにより、位置測定精度を向上させることができる。また基準距離Ｌは、盛立構造物１０２の変形によって移動した移動用三軸センサ１１０と基準用三軸センサ１１２の距離（予測される距離）としてもよい。特に大変形時には、発信機の予測される位置を基準とすることで、さらに位置測定精度を向上させることができる。
【００３７】
なお、較正用三軸センサ１１４を設置する位置は、移動用三軸センサ１１０と同じ盛立高さであってもよいが、移動用三軸センサ１１０の位置が明確である限りにおいて、盛立が進んだ位置（高さ）であってもよい。また較正用三軸センサ１１４を設置する位置の土質は、移動用三軸センサ１１０と連続する同じ土質であってもよいが、異なる土質の位置であってもよい。その場合は、それぞれの土質についての導電率を境界条件として入力すればよい。
【００３８】
そして、移動用三軸センサ１１０を励磁して、較正用三軸センサ１１４に誘起電圧Ｅ１（Ｅx1、Ｅy1、Ｅz1）を発生させる。較正用三軸センサ１１４に生じた誘起電圧Ｅ１は、信号発信部２３２から発信された信号を受信する信号受信部１１６が取得する（Ｓ１５６）。
【００３９】
次に、オペレータ（使用者）が、入力部１１８から、盛立構造物１０２の導電率σを入力する（Ｓ１５８）。導電率σは、既定のデータ（設計条件、文献調査等）から設定してもよいが、盛立構造物１０２を測定した実測値であるのが好ましい。これは、現場に則した導電率を設定できるからである。また、三軸センサが配置される周辺の値を取得するのがさらに好ましい。
【００４０】
地盤の電導率を測定する方法としては古くから様々な手法が知られているが、本実施形態では一例としてウェンナ法について説明する。図５は、盛立構造物１０２の導電率の実測手法の一例を示す図である。図５に示すように、盛立構造物１０２の導電率を実測するには、その測定対象域に直線上に４本の電極１４０、１４２、１４４、１４６を距離ａの等間隔で配置し、交流電源１４８によって両端の電極１４０、１４６に交流電圧を印加する。そして、電流計１５０が計測した電流ｉ（電極１４０、１４６間を流れる電流）と、電圧計１５２が計測した電圧ｖ（電極１４２、１４４間に生じる電圧）から、以下の式（式５）より導電率が算出される。
σ＝ｉ／（２×π×ａ×ｖ） （式５）
【００４１】
図３に戻って、次に、較正用磁場導出部１２２が、移動用三軸センサ１１０から出力された磁場について、基準距離Ｌおよび導電率σを用いて第１の磁場解析を行い、較正用三軸センサ１１４の位置における磁場強度Ｈ１（Ｈx1、Ｈy1、Ｈz1）を導出する（Ｓ１６０）。
【００４２】
較正用磁場導出部１２２は、有限要素法解析（ＦＥＭ）や境界要素法解析（ＢＥＭ）などの数値解析を用いて、磁場強度Ｈ１を求めることができる。移動用三軸センサ１１０が出力する磁場強度Ｈ０は励磁コイル２１６をリング状のコイル要素でモデル化して流れる電流を入力することで理論的に求めることができるから、盛立構造物１０２の形状、材質、移動用三軸センサ１１０から較正用三軸センサ１１４までの基準距離、および入力された導電率を境界条件として設定することにより、磁場強度Ｈ１を求めることができる。このように、盛立構造物１０２の導電率σを反映した第１の磁場解析を行うことで、導電率σが大きい場合や変化する場合にも影響を適切に考慮することができる。なお、励磁コイル２１６から生じる磁場の解析においては、渦電流の影響を考慮した定常渦電流解析を行うことが好ましい。
【００４３】
第１の磁場解析における解析モデル（盛立構造物１０２の形状）は、図４（ａ）に示すように盛立途中の形状としてもよく、または図４（ｂ）に示すように盛立完了後の形状としてもよい。盛立途中の形状では、三軸センサ（移動用三軸センサ１１０、較正用三軸センサ１１４）の上側が気中になっている条件を反映して、より実際の条件に近い磁場強度を求めることができる。しかし、完了後の形状で磁場解析しても磁場強度への影響は大差なく、十分な精度で磁場強度が求められることが分かっている。さらに、完了後の形状で解析モデルを作成しておけば、後述の第２の磁場解析にもその解析モデルを流用できる利点も有する。
【００４４】
次に、換算比率導出部１２４が、磁場強度Ｈ１と誘起電圧Ｅ１とから、換算比率Ｋ（Ｋｘ＝Ｈx1／Ｅx1、Ｋｙ＝Ｈy1／Ｅy1、Ｋｚ＝Ｈz1／Ｅz1）を導出する（Ｓ１６２）。換算比率Ｋは、土質の導電率が変化しなければ、距離にかかわらず一定となる。したがって較正によって地盤の導電率σに対応した磁場強度Ｈ１と誘起電圧Ｅ１との換算比率Ｋを求めておけば、運用時に誘起電圧Ｅ２を測定した場合に、これに対応する磁場強度Ｈ２を求めることができる。
【００４５】
なお、導電率σや換算比率Ｋは、ＲＯＭ、ＲＡＭ、ＥＥＰＲＯＭ、不揮発性ＲＡＭ、フラッシュメモリ、ＨＤＤ等からなる記憶部１２０に記憶される。
【００４６】
次に、複数（多数）の移動用三軸センサ１１０が埋設された状態で、盛立構造物１０２の盛立完了まで盛立が実施される（Ｓ１６４）。なお、較正用三軸センサ１１４に関しては、Ｓ１５６で誘起電圧Ｅ１を取得した後に盛立構造物１０２から撤去してよく、基準位置に設置して基準用三軸センサ１１２として利用してもよい。校正時と本計測時とで同じセンサを用いることで、センサ特性を反映した校正ができ、さらに位置測定精度を向上させることができる。また受信機能だけでなく送信機能も備えているのであれば、較正後にそのまま埋設して移動用三軸センサ１１０として利用することもできる。
【００４７】
そして、図４（ｂ）に示すように、盛立構造物１０２の表面近傍（基準位置）に基準用三軸センサ１１２が配置される（Ｓ１６６）。基準用三軸センサ１１２は１つであってもよいが、移動用三軸センサ１１０からの磁場の到達距離と盛立構造物１０２の規模を勘案して、移動用三軸センサ１１０をグルーピングし、それぞれのグループを担当する複数の基準用三軸センサ１１２を設けてもよい。さらに、同じ移動用三軸センサ１１０に対して異なる位置に配置した２以上の基準用三軸センサ１１２から位置の同定を行うことにより、それぞれが同定した位置座標の平均処理を行うことで、測定精度をさらに向上させることができる。
【００４８】
本実施形態では、移動用三軸センサ１１０すなわち発信機を埋設し、基準用三軸センサすなわち受信機を表面近傍に配置している。よって、発信機に電源を内蔵し、リモート制御とすれば盛立構造物１０２内へのケーブル配線の設置を回避できる。一方、受信機を埋設する場合には、信号発信部２３２が誘起電圧の計測値を信号受信部１１６へと発信する必要がある。そのため、電源を内蔵するには、消費電力の点で課題があり、信号受信部１１６へと信号を伝達するケーブル配線を敷設しなければならない可能性が高い。その上、複数の発信機を埋設する場合には１つの受信機で複数の発信機からの磁場強度を測定することができるため、受信機の数が少なくて済み、データ処理も容易となるため、システム構成を簡略化することができる。
【００４９】
基準用三軸センサ１１２を設置する基準位置は、その３次元座標（絶対座標：地球座標とも言う）が明確であればよく、盛立構造物１０２の表面近傍ばかりではなく、ある程度埋設された位置であってもよい。また基準位置は、例えばＧＰＳを用いたり、路面上の印を基準として測量したりして、随時更新（確認）することが好ましい。
【００５０】
運用時においては、信号受信部１１６が、移動用三軸センサ１１０を励磁した際に基準用三軸センサ１１２に生じる誘起電圧Ｅ２（Ｅx2、Ｅy2、Ｅz2）を取得する（Ｓ１６８）。そして、磁場強度換算部１２６が、誘起電圧Ｅ２と記憶部１２０に記憶された換算比率Ｋ（Ｋｘ、Ｋｙ、Ｋｚ）を用いて、以下の式（式６）より基準用三軸センサ１１２の位置における磁場強度Ｈ２（Ｈx2、Ｈy2、Ｈz2）を求める（Ｓ１７０）
Ｈ２＝Ｋ・Ｅ２＝（Ｈ１／Ｅ１）・Ｅ２ （式６）
なお、これを、ｘ成分、ｙ成分、ｚ成分別にすると、以下の式（式７〜式９）に表される。
Ｈx2＝Ｋｘ・Ｅx2＝（Ｈx1／Ｅx1）・Ｅx2 （式７）
Ｈy2＝Ｋｙ・Ｅy2＝（Ｈy1／Ｅy1）・Ｅy2 （式８）
Ｈz2＝Ｋｚ・Ｅz2＝（Ｈz1／Ｅz1）・Ｅz2 （式９）
【００５１】
次に、位置同定部１２８が、導電率σを反映した第２の磁場解析によって、基準用三軸センサ１１２の位置における磁場強度が算出した磁場強度Ｈ２となるように、移動用三軸センサ１１０の現在位置を同定する（Ｓ１７２）。第２の磁場解析における解析モデル（盛立構造物１０２の形状）は、図４（ｂ）に示すように盛立完了後の形状とする。詳細には後述するが、本実施形態では移動用三軸センサ１１０が発信機であるため、位置同定部１２８は、逆解析的手法に基づき移動用三軸センサ１１０の現在位置を同定する。なお、移動用三軸センサ１１０が受信機である場合には、通常の磁場解析によって、その現在位置を同定することができる。そして位置同定部１２８は、基準用三軸センサ１１２の基準位置に対する移動用三軸センサ１１０の現在位置の位置関係と、基準位置の絶対座標とから、移動用三軸センサ１１０の現在位置を絶対位置として同定することができる。
【００５２】
次に、変形計測部１３０が、盛立構造物１０２の変形を計測する。具体的には、移動用三軸センサ１１０の現在位置の初期位置からの変位を計算することにより、移動用三軸センサ１１０の移動量、すなわち盛立構造物１０２のその箇所の変形量を計測することができる。このような計測によって、盛立構造物１０２の変形が所定以上であると判断された場合（Ｓ１７４のＹｅｓ）には、報知部１３２がオペレータに対して警告を発する（Ｓ１７６）。また所定未満であった場合には（Ｓ１７４のＮｏ）、引き続き盛立構造物１０２の変形の計測が継続される。
【００５３】
（移動用三軸センサ１１０の現在位置の同定）
以下、位置同定部１２８による移動用三軸センサ１１０の現在位置の同定手法について説明する。ここで、上述したように、磁場強度換算部１２６が導出した磁場強度Ｈ２（Ｈx2、Ｈy2、Ｈz2）および基準用三軸センサ１１２の基準位置は既知である。
【００５４】
有限要素法解析（ＦＥＭ）や境界要素法解析（ＢＥＭ）などの数値解析は、空間をメッシュ状に細かく分割し、その１つの要素（区画）ごとに順に計算を行う。詳細については、後述する実施例に記載するが、このメッシュの１辺の長さを移動用三軸センサ１１０が備えるコイルの直径の１／４以上１／２以下とし、移動用三軸センサ１１０の初期位置から少なくとも２０ｍの範囲内の空間を他の範囲よりも細かく分割するとよい。これにより、移動用三軸センサ１１０の位置の同定精度を実用的なレベルまで向上できる。
【００５５】
現在位置を同定する１つ目の手法としては、繰り返し計算によって正解に漸近する方法が挙げられる。すなわち、盛立構造物１０２のモデルのなかで移動用三軸センサ１１０を仮想発信点に設定し、仮想発信点の位置を変えながら基準用三軸センサ１１２の基準位置における磁場強度Ｈ’（Ｈx’、Ｈy’、Ｈz’）を算出し、算出した磁場強度Ｈ’と先に導出された磁場強度Ｈ2（Ｈx2、Ｈy2、Ｈz2）との残差が最小となるような発信機（移動用三軸センサ１１０）の位置を現在位置を同定する。残差は、次の式（式１０）によって表される。
（Ｈx2−Ｈx’）２＋（Ｈy2−Ｈy’）２＋（Ｈz2−Ｈz’）２＝ｍｉｎ （式１０）
ここで、あらかじめ複数の仮想発信点について磁場強度Ｈ’を算出して残差が最小となる仮想発信点を選択してもよいが、残差が少なくなる方向に仮想発信点をずらしながら残差が十分に小さくなるまで計算を繰り返すことが好ましい。
【００５６】
図６は、移動用三軸センサ１１０の現在位置を同定する手法について説明する図である。以下、図６を参照しながら、現在位置を同定する２つ目の手法について説明する。
【００５７】
図６に示すように、まず移動用三軸センサ１１０が埋設された初期位置の周囲に複数の仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａを設定する（実際には立体であるため８点をとるべきであるが、図では説明の便宜上４点をとって図示および説明をする。）。この仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａは、移動用三軸センサ１１０の初期位置からの変位を想定した上で移動用三軸センサ１１０の現在位置を囲うように設定される。移動後の現在位置は未知であるが、移動範囲は想定可能である。
【００５８】
そして、それぞれの仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａ（図６の白丸）に移動用三軸センサ１１０があると仮定した場合に、導出された磁場強度Ｈ2となる複数の仮想受信点１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂ（図６の黒丸）を、導電率σの影響を考慮した第２の磁場解析によって導出する。なお、第２の磁場解析を行うにあたり、仮想発信点で囲まれた範囲内の分割幅は細かくすることが好ましく、例えば１ｍ以下とすることができる。
【００５９】
４つの仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａで囲まれた領域の形状は、４つの仮想受信点１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂで囲まれた領域の形状とは一致しない。これは、導電率σによって磁場が歪みを生じるためである。しかし、仮想発信点同士と仮想受信点同士の相対位置関係は比例関係にあると考えることができる。そこで、仮想発信点同士と仮想受信点同士によって囲まれる空間を同じ数で分割する。そして、仮想受信点１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂに対する基準用三軸センサ１１２の位置関係を、仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａに当てはめることによって、比例配分により移動用三軸センサ１１０の現在位置を同定することができる。
【００６０】
すなわち、図６では、仮想受信点１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂで囲まれた領域の右から２つ目、上から３つ目のマス目に基準用三軸センサ１１２が配置されているので、移動用三軸センサ１１０の現在位置も仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａで囲まれた領域の右から２つ目、上から３つ目のマス目であると同定することができる。故に、１つ目の手法のように繰り返し計算を行うよりも、簡潔かつ高速に移動用三軸センサ１１０の現在位置を同定することができる。
【００６１】
図７は、移動用三軸センサ１１０の現在位置を同定する手法について説明する図である。以下、図７を参照しながら、現在位置を同定する３つ目の手法について説明する。
【００６２】
図７に示すように、まず、移動用三軸センサ１１０が埋設された初期位置の周囲に複数の仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａを設定する。上記と同様に、移動用三軸センサの現在位置を囲うように仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａは設定される。そして、それぞれの仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａに移動用三軸センサ１１０があると仮定した場合に、基準用三軸センサ１１２の基準位置での磁場強度を、導電率σの影響を考慮した第２の磁場解析によって導出する。すると、仮想発信点ごとに異なる値の磁場強度が得られる（基本的には、遠いものほど弱くなる）。
【００６３】
ここで、図７に示すように、仮想発信点１８０ａからの磁場強度が１００、仮想発信点１８２ａからの磁場強度が１００、仮想発信点１８４ａからの磁場強度が８０、仮想発信点１８６ａからの磁場強度が８０、そして磁場強度換算部１２６が導出した磁場強度Ｈ２が９０であったとする。このとき、それぞれの仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａに対応する磁場強度と導出された磁場強度Ｈ2の大小関係を定めることができるので、その大小関係に対応する比例配分により移動用三軸センサ１１０の現在位置を同定することができる。この場合においても、１つ目の手法のように繰り返し計算を行うよりも、簡潔かつ高速に移動用三軸センサ１１０の現在位置を同定することができる。
【００６４】
以上、本実施形態にかかる変形計測システム１００および変形計測方法によれば、盛立構造物１０２の導電率の影響を考慮して移動用三軸センサ１１０と基準用三軸センサ１１２の相対位置を計測することにより、導電率が大きい場合や、発信機と受信機の間で導電率が変化する（土質をまたぐ）場合であっても、盛立構造物１０２の変形を高精度に計測することができる。なお、この変形計測システム１００および変形計測方法は、地中、水中問わず適用することができ、光学計測機器では測定不可能な盛立構造物１０２の堤体内部や水中における測定も可能である。そのため、崩壊に至る前兆を確実に高精度に検出することができ、ひずみ発生段階での警報を行うことも可能になる。
【００６５】
なお、本実施形態では、移動用三軸センサ１１０を発信機、基準用三軸センサ１１２、および較正用三軸センサを受信機として説明した。しかし、移動用三軸センサ１１０を受信機、基準用三軸センサ１１２および較正用三軸センサを発信機としてもよい。このような場合には、埋設される移動用三軸センサ１１０の現在位置を逆解析的手法により求める必要がないので、数値解析の回数を削減でき、また同定精度のさらなる向上を図り得る。
【００６６】
（磁場解析方法の詳細）
次に、磁場解析方法の詳細について説明する。
【００６７】
（１）Ｍａｘｗｅｌｌの方程式
磁場解析における電磁気の現象を記述する構成則として、Maxwellの方程式を採用した。この方程式は、次の４式から構成される。
【００６８】
【数１】
【００６９】
また、物質の構成式は、次のように表される。
【００７０】
【数２】
【００７１】
（２）解析手法の検討
送信アンテナから発せられる交流の周波数は約１０ｋＨｚであるとすると、電磁波の進行速度は概ね光速（約３０万ｋｍ／ｓ）であることから、伝搬する電磁波の１波長は約３０ｋｍとなる。これに対して、計測対象範囲は送受信間距離が高々数十ｍであることから、電磁波の一波長に比べると非常に短い区間が対象となる。このため、送信アンテナから発せられる磁場の振幅と受信側のアンテナで計測される磁場の振幅の低下(振動としての振幅の減少)を捉えることが容易である。このような視点から見た電磁波は、波動として捉える動的問題から、ポテンシャルが低下する静的問題に置き換えることが可能である。
【００７２】
電磁波の問題を動的問題から静的問題へ転換した場合、電場の変化に伴う変位電流（過渡的な電流）を考慮しないことから、式１１〜式１４に示したＭａｘｗｅｌｌの方程式の内、式１１から式１３を連立して解くこととなる。また、同様の理由により、式１２の右辺第三項は恒等的に０として取り扱うこととなる。所謂、定常渦電流解析を実施することとなる。
【００７３】
（３）定常渦電流解析による導電率σを考慮した真の磁場強さＨの算出方法
磁場が時間的に変化すると電磁誘導により電界が生じる。電界の中に導体が置かれていると、その中には渦電流が発生する、と同時に渦電流により磁場が発生する。本解析では、この渦電流の影響を考慮した。
【００７４】
任意のベクトルをＡと書き、Ａの回転を取った式１６のベクトル▽×Ａを考えると、任意のベクトルに回転を取ったベクトル▽×Ａは、常に式１３を満たすことが分かっている。
【００７５】
【数３】
【００７６】
電磁気学の世界では、このベクトルＡをベクトルポテンシャルと呼ぶ。▽×Ａが式１３を満足することから、逆に、式１７に示すように、これを磁束密度Ｂとおくことができる。
Ｂ＝▽×Ａ （式１７）
【００７７】
電界の強さとベクトルポテンシャルＡ、電位φ（一般に言われる電圧）との間には、式１８に示す関係があることが分かっている。
【００７８】
【数４】
【００７９】
式１７と式１８を式１２に代入すると、定常渦電流の基礎式１９を得る。
【００８０】
【数５】
【００８１】
式１７、式１８を仮定することにより、式１１は数学的に自動的に満足されることから、具体的には、式１９を解くことによりＭａｘｗｅｌｌの方程式の解に帰着することとなる。
【００８２】
式１９はφとＡを未知数とする方程式であり、解を得るためには、Ａ、φに関する他の関係式が必要となる。これについては、式２０に示す電流に対するガウスの法則（定常電流では、導体の閉曲面に出入りする電流の総和は常に０となる）を用いる。
【００８３】
【数６】
【００８４】
上記２式が渦電流解析で満足すべき方程式となり、解析的には両式を連立させて、境界を貫く磁束密度が０となる境界条件の下で、２つの未知数Ａとφを求めることとなる。μ、σ、Ｊ０は与条件であるため、Ａが得られれば、導電率σに応じた任意点の磁場の強度Ｈが、式２１として解析から求められることとなる。
【００８５】
【数７】
【００８６】
発信器は互いに直交する三方向のコイルに順に励磁するため、解析も発信器と同じ三方向のコイル（仮にコイルｘ、コイルｙ、コイルｚとする）に対して実施すれば、コイルｘによる磁場の強さＨｘ、コイルｙによる磁場の強さＨｙ、コイルｚによる磁場の強さＨｚがそれぞれ求められることになる。
【実施例】
【００８７】
（コイル直径Ｄとメッシュの１辺の長さΔｘとの比に起因する解析精度）
数値解析において移動用三軸センサ１１０の現在位置を同定する際のメッシュの大きさとコイル寸法の関係を把握するため、コイル近傍のコイル理論値と解析結果を比較した。
【００８８】
表１に示すように、コイル直径Ｄ＝１．０ｍ（半径ｒ＝０．５ｍ）として、Ｄ／Δｘ（メッシュの１辺の長さ）＝１．０〜４．２となるようにΔｘを変化させた。解析モデルとしては、一辺５ｍの立方体を用い、コイル中心位置を解析領域の中心とした。コイル電流は１０００Ａとした。
【００８９】
表１に示すように、Ｄ／Δｘが大きくなるほど誤差が小さくなり、磁場強度の解析値（表１の値は解析モデル中の代表点）が理論値に近づくことが分かる。また、Ｄ／Δｘが２以上であれば十分な解析精度は得られると判断し得る。Ｄ／Δｘが４以下であれば、実用的な解析時間に収まることから、Ｄ／Δｘは２〜４とするのが好ましいと判断される。なお、ここでの誤差は以下の式（式２２）より算出される。
誤差＝｛解析値−理論値｝／理論値 （式２２）
【００９０】
【表１】
【００９１】
（メッシュの細分化領域）
メッシュの細分化領域（１〜２ｍ間隔でメッシュを細かく分割する領域）の大きさと解析精度との関係を把握するため、細分化領域の大きさを変化させて解析結果と理論値とを比較した。
【００９２】
図８は、実施例および比較例の解析モデルを示すである。図８に示すように、実施例では、送信コイルから１０ｍまでは１ｍ間隔でメッシュを細分割し、送信コイルから１０ｍ〜２０ｍ間を２ｍ間隔でメッシュを細分割し、それ以降は分割領域を増加させた。比較例では送信コイルから１０ｍまでは１ｍ間隔でメッシュを細分割したが、それ以降は分割間隔を増加させた。コイル直径Ｄ＝３．０ｍとし、コイル軸と解析モデルのｘ軸の角度を４５°とした。
【００９３】
図９は、実施例および比較例における発信点からの距離と磁場強度の解析誤差の相関を示す図である。図９に示すように、実施例では発信点からの距離に関わらず、磁場強度の解析誤差を殆ど０％とすることができる。しかし、比較例では、メッシュ分割が粗くなる発信点から１０ｍを境界として解析誤差が急増しており、その誤差は６０％程度にも達し得る。よって、解析精度を維持するためには、発信点（励磁コイル２１６）の位置から少なくとも２０ｍの範囲は空間を細かく分割する必要があると判断できる。
【００９４】
以上、本発明の好適な実施形態、実施例について詳細に説明したが、本発明は係る実施形態、実施例に限定されないことは言うまでもない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【産業上の利用可能性】
【００９５】
本発明は、２以上の三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測システムおよび変形計測方法に利用することができる。
【符号の説明】
【００９６】
１００…変形計測システム、１０２…盛立構造物、１０８ａ…三軸センサ（発信機）、１０８ｂ…三軸センサ（受信機）、１１０…移動用三軸センサ、１１２…基準用三軸センサ、１１４…較正用三軸センサ、１１６…信号受信部、１１８…入力部、１２０…記憶部、１２２…較正用磁場導出部、１２４…換算比率導出部、１２６…磁場強度換算部、１２８…位置同定部、１３０…変形計測部、１３２…報知部、１４０、１４２、１４４、１４６…電極、１４８…交流電源、１５０…電流計、１５２…電圧計、１８０ａ、１８２ａ、１８４ａ、１８６ａ…仮想発信点、１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂ…仮想受信点、２１２…発信回路、２１４…増幅器、２１６…励磁コイル、２１８…切替回路、２２２…検出コイル、２２４…切替回路、２２６…検波回路、２２８…電圧計、２３０…誘起電圧演算部、２３２…信号発信部
【技術分野】
【０００１】
本発明は、２以上の三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測システムおよび変形計測方法に関するものである。
【背景技術】
【０００２】
フィルダムや堤防の堤体に代表される盛立構造物においては、安全上の理由から（崩壊の危険性を把握するために）施工中および完成後その変形を監視する必要がある。そのため、従来、基礎岩盤に固定された支柱やワイヤ等によって連結された沈下計などの埋設型センサを埋設して、それらの相対位置を計測することによって盛立構造物の変形を計測していた。
【０００３】
しかし、支柱やワイヤは盛立構造物にとっては異物であり、強度等の性質低下を招きかねない。また、かかる手法では、支柱やワイヤと埋設型センサとの相対変位、すなわち１次元の変位しか計測することができない。その上、従来の埋設型センサは、信号を入出力するためにケーブルを要するので、工事中におけるケーブルの養生やメンテナンスの労力といった課題を生じることとなる。さらに、これらのケーブルが水路（みずみち：水の通り路）を形成するおそれもある。
【０００４】
そこで、このような課題を解決する技術が求められている。これに関連する従来技術としては、特許文献１に、三軸センサを２点に配置して一方の三軸センサを励磁した場合に他方の三軸センサに誘起される電圧（Ｅ１）の真空中における距離減衰特性を利用して、２点間の相対的な距離（ｒ１）を以下の式（式１）より算出する技術が開示されている。かかる技術は、地中や水中であっても適用可能とされている。
ｒ１＝｛（Ｅ０／Ｅ１）１／３｝・ｒ０ （式１）
ここで、ｒ０：既知の距離
Ｅ０：２点間の距離がｒ０の場合に三軸センサに誘起される電圧
ｒ１：２点間の距離（算出対象）
Ｅ１：算出対象時において三軸センサに誘起される電圧
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００３−４４０９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１の技術は、真空中すなわち導電率の影響を排除できる環境を対象としている。現実的には、盛立構造物の地中においては、その構成要素である岩石、砂利、水等の導電率（例えばロックフィルダム堤体の導電率は０．００１から０．１Ｓ/ｍ）により磁場の分布に歪みを生じるので、当然ながらその影響を考慮する必要がある。
【０００７】
よって、導電率の影響を排除し得ない盛立構造物に対し、特許文献１の技術を適用した場合には、測定に無視できない誤差を生じることとなる。すなわち、この技術は精度の点に課題がある。その上、特にフィルダムでは中央から外側に行くにつれて土質（岩石、砂利、土等）を異ならせるため、箇所によって導電率が異なり、また三軸センサの送信側と受信側の間に土質の境界が存在する場合もあるため、既定のデータ等から一律に導電率を定めることに対して難がある。
【０００８】
そこで、本発明は、盛立構造物の導電率の影響を考慮して２以上の三軸センサの相対位置を計測することにより、盛立構造物の変形を高精度に計測することが可能な変形計測システムおよび変形計測方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するために本発明の代表的な構成は、軸方向が異なる３つのコイルを備え、コイルを励磁することによって磁場を生じ、またはコイルが磁場を受けて発生した誘起電圧を検出して信号として出力する三軸センサを用いて、２以上の三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測システムであって、盛立構造物に埋設される移動用三軸センサと、盛立構造物の表面近傍に配置される基準用三軸センサと、移動用三軸センサから所定の基準距離に配置される較正用三軸センサと、移動用三軸センサまたは基準用三軸センサ、もしくは較正用三軸センサから出力された信号を受信する信号受信部と、盛立構造物の導電率を入力する入力部と、移動用三軸センサから出力された磁場について基準距離および導電率を用いて第１の磁場解析を行い、較正用三軸センサの位置における磁場強度Ｈ１を導出する較正用磁場導出部と、較正用三軸センサの位置における磁場強度Ｈ１とこの較正用三軸センサが検出した誘起電圧Ｅ１とから、盛立構造物において磁場強度と誘起電圧とを換算する換算比率を導出する換算比率導出部と、移動用三軸センサまたは基準用三軸センサの一方が検出した誘起電圧Ｅ２と換算比率とから、誘起電圧Ｅ２を検出した三軸センサの位置における磁場強度Ｈ２を求める磁場強度換算部と、導電率を反映した第２の磁場解析を行い、移動用三軸センサまたは基準用三軸センサのうち、誘起電圧Ｅ２を検出した一方の三軸センサの位置における磁場強度が磁場強度換算部が求めた磁場強度Ｈ２となるように、他方の三軸センサの現在位置を同定する位置同定部と、を備えることを特徴とする。
【００１０】
かかる構成によれば、盛立構造物の導電率の影響を考慮して２以上の三軸センサの相対位置を計測することが可能となり、盛立構造物の変形を高精度に計測することができる。特に、磁場解析による磁場強度と較正用に実測した誘起電圧によって求めた換算比率を用いることにより、その盛立構造物における磁場の強度と誘起電圧の関係を求めることができ、埋設された三軸センサの位置を同定することができる。
【００１１】
上記導電率は、盛立構造物を測定した実測値であるとよい。上述したように、盛立構造物は、地点によって構成要素（岩石、砂利、水等）が異なるため、既定のデータ等から一律に導電率を定めることが困難である。一方、盛立構造物の実測によれば、現場に則した値を取得できるため最適な導電率を設定することができる。
【００１２】
上記移動用三軸センサは発信機、基準用三軸センサは受信機であって、位置同定部は、移動用三軸センサが埋設された初期位置の周囲に複数の仮想発信点を設定し、それぞれの仮想発信点に移動用三軸センサがあると仮定した場合に、この移動用三軸センサから出力された磁場について第２の磁場解析を行うことにより磁場強度換算部が求めた磁場強度Ｈ２となる複数の仮想受信点を求め、複数の仮想受信点と基準用三軸センサの位置関係と、複数の仮想発信点の設定位置とから、比例配分により移動用三軸センサの現在位置を同定するとよい。これにより、簡潔かつ容易に移動用三軸センサの現在位置を同定することができる。
【００１３】
上記移動用三軸センサは発信機、基準用三軸センサは受信機であって、位置同定部は、移動用三軸センサが埋設された初期位置の周囲に複数の仮想発信点を設定し、それぞれの仮想発信点に移動用三軸センサがあると仮定した場合に、この移動用三軸センサから出力された磁場について第２の磁場解析を行うことにより基準用三軸センサの位置における複数の磁場強度を求め、複数の磁場強度および磁場強度換算部が求めた磁場強度Ｈ２との大小関係と、複数の仮想発信点の設定位置とから、比例配分により移動用三軸センサの現在位置を同定してもよい。これにより、簡潔かつ容易に移動用三軸センサの現在位置を同定することができる。
【００１４】
上記位置同定部は、空間をメッシュ状に分割する３次元解析によって第２の磁場解析を行い、メッシュの１辺の長さは、移動用三軸センサが備えるコイルの直径の１／４以上１／２以下に設定するとともに、初期位置から少なくとも２０ｍの範囲内の空間は他の範囲よりも細かく分割するとよい。これにより、高精度に移動用三軸センサの現在位置の同定を行うことができる。
【００１５】
上記盛立構造物は、土砂や岩石を盛り立てて構成されるフィルダムの堤体であるとよい。すなわち、当該変形計測システムは、フィルダムを主な適用対象とする。
【００１６】
上記課題を解決するために本発明の他の代表的な構成は、軸方向が異なる３つのコイルを備え、コイルを励磁することによって磁場を生じ、またはコイルが磁場を受けて発生した誘起電圧を検出して信号として出力する三軸センサを用いて、２以上の三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測方法であって、盛立構造物の初期位置に移動用三軸センサを埋設し、盛立構造物の導電率を測定し、移動用三軸センサから所定の基準距離の位置に較正用三軸センサを配置し、較正用三軸センサを用いて誘起電圧Ｅ１を測定し、移動用三軸センサから出力される磁場について基準距離および導電率を用いて第１の磁場解析を行い、較正用三軸センサの位置における磁場強度Ｈ１を導出し、較正用三軸センサの位置における磁場強度Ｈ１とこの較正用三軸センサが検出した誘起電圧Ｅ１とから、盛立構造物において磁場強度と誘起電圧とを換算する換算比率を導出し、盛立構造物の表面近傍に基準用三軸センサを配置し、移動用三軸センサまたは基準用三軸センサの一方によって誘起電圧Ｅ２を測定し、誘起電圧Ｅ２と換算比率とから、誘起電圧Ｅ２を検出した三軸センサの位置における磁場強度Ｈ２を導出し、導電率を反映した第２の磁場解析を行い、移動用三軸センサまたは基準用三軸センサのうち、誘起電圧Ｅ２を検出した一方の三軸センサの位置における磁場強度が磁場強度換算部が求めた磁場強度Ｈ２となるように、他方の三軸センサの位置を同定するとよい。
【００１７】
上述した変形計測システムにおける技術的思想に対応する構成要素やその説明は、当該変形計測方法にも適用可能である。
【発明の効果】
【００１８】
本発明によれば、盛立構造物の導電率の影響を考慮して２以上の三軸センサの相対位置を計測することにより、盛立構造物の変形を高精度に計測することが可能な変形計測システムおよび変形計測方法を提供することができる。
【図面の簡単な説明】
【００１９】
【図１】三軸センサについて説明する図である。
【図２】本実施形態にかかる変形計測システムの構成を示すブロック図である。
【図３】変形計測システムが盛立構造物の変形を計測する一連の動作を説明するフローチャートである。
【図４】盛立構造物の盛立途中および盛立完了後を示す図である。
【図５】盛立構造物の導電率の実測手法の一例を示す図である。
【図６】移動用三軸センサの現在位置を同定する手法について説明する図である。
【図７】移動用三軸センサの現在位置を同定する手法について説明する図である。
【図８】実施例および比較例の解析モデルを示す図である。
【図９】実施例および比較例における発信点からの距離と磁場強度の解析誤差の相関を示す図である。
【発明を実施するための形態】
【００２０】
以下に添付図面を参照しながら、本発明の好適な実施形態、実施例について詳細に説明する。かかる実施形態、実施例に示す寸法、材料、その他具体的な数値などは、発明の理解を容易とするための例示に過ぎず、特に断る場合を除き、本発明を限定するものではない。なお、本明細書及び図面において、実質的に同一の機能、構成を有する要素については、同一の符号を付することにより重複説明を省略し、また本発明に直接関係のない要素は図示を省略する。
【００２１】
（三軸センサ）
図１は、三軸センサ１０８ａ、１０８ｂについて説明する図である。特に、図１（ａ）は発信機である三軸センサ１０８ａについて説明する図、図１（ｂ）は受信機である三軸センサ１０８ｂについて説明する図である。以下、三軸センサ１０８ａ、１０８ｂについて説明するが、その詳細については特許文献１（特開２００３−４４０９号公報）に開示されているため以下簡略化して説明する。
【００２２】
図１（ａ）に示すように、三軸センサ１０８ａ（発信機）は、発信回路２１２、増幅器２１４、励磁コイル２１６、切替回路２１８を包含する。
【００２３】
発信回路２１２は、所定の周波数の交流信号を発生させて、増幅器２１４へと出力する。このときの周波数は、高くなると盛立構造物に吸収されやすくなってしまい、低すぎると検波しにくくなるため、可聴帯域（２０Ｈｚ〜１５ｋＨｚ程度）を用いることができ、さらに好適には１０ｋＨｚ程度を用いることができる。増幅器２１４は、発信回路２１２から出力された交流信号の電流量を増幅して、励磁コイル２１６へと出力する。
【００２４】
励磁コイル２１６は、軸方向が互いに直交する三軸座標系のそれぞれの軸に配置された３つコイル（ｘ,ｙ,ｚ）を備えており、それぞれのコイルは増幅器２１４から与えられた交流信号により励磁され磁場を生じる。切替回路２１８は、交流信号の出力対象とされるコイル（ｘ,ｙ,ｚ）を順次切り替える。
【００２５】
図１（ｂ）に示すように、三軸センサ１０８ｂ（受信機）は、検出コイル２２２、切替回路２２４、検波回路２２６、電圧計２２８、誘起電圧演算部２３０、信号発信部２３２を包含する。
【００２６】
検出コイル２２２は、軸方向が互いに直交する三軸座標系のそれぞれの軸に配置された３つのコイル（ｕ,ｖ,ｗ）であって、励磁コイル２１６（ｘ,ｙ,ｚ）が発生する磁場に伴い誘起電圧を発生する。それぞれのコイル（ｕ,ｖ,ｗ）に発生した誘起電圧は、検波回路２２６を介して電圧計２２８により測定される。切替回路２２４は、電圧計２２８の測定対象とされるコイル（ｕ,ｖ,ｗ）を順次切り替える。
【００２７】
よって、発信機の励磁コイル２１６（ｘ,ｙ,ｚ）をそれぞれ励磁した場合、受信機ではコイルｘの励磁に起因するそれぞれのコイル（ｕ,ｖ,ｗ）の誘起電圧Ｅxu、Ｅxv、Ｅxw、コイルｙの励磁に起因するそれぞれのコイル（ｕ,ｖ,ｗ）の誘起電圧Ｅyu、Ｅyv、Ｅyw、コイルｚの励磁に起因するそれぞれのコイル（ｕ,ｖ,ｗ）の誘起電圧Ｅzu、Ｅzv、Ｅzwの合計９種類の電圧が測定される。
【００２８】
誘起電圧演算部２３０は、これら９つのパラメータ（誘起電圧）を用いて、受信機が配置される位置における誘起電圧のｘ成分Ｅx、ｙ成分Ｅy、ｚ成分Ｅzを以下の式（式２から式４）に基づき算出する。
Ｅx＝（Ｅxu＋Ｅxv＋Ｅxw）1/2 （式２）
Ｅｙ＝（Ｅyu＋Ｅyv＋Ｅyw）1/2 （式３）
Ｅｚ＝（Ｅzu＋Ｅzv＋Ｅzw）1/2 （式４）
【００２９】
信号発信部２３２は、誘起電圧演算部２３０が算出した誘起電圧を、後述する信号受信部１１６へと発信する。なお、発信側の三軸センサ１０８ａにも不図示の信号受信部を備えて、算出した誘起電圧を発信側の三軸センサ１０８ａを介して信号受信部１１６へと発信してもよい。本実施形態では、この誘起電圧Ｅ（Ｅx、Ｅy、Ｅz）を用いて、発信機が配置される位置と、受信機が配置される位置の間の距離（相対位置）を磁場解析により求める。
【００３０】
なお上記説明では三軸センサ１０８ａ、１０８ｂの構成を送信機と受信機として説明したが、コイルを共用として両方の構成を備えることにより、送受信のいずれも可能な三軸センサとすることができる。
【００３１】
（変形計測システムおよび変形計測方法）
図２は、本実施形態にかかる変形計測システム１００の構成を示すブロック図である。図２に示すように、変形計測システム１００は、移動用三軸センサ１１０、基準用三軸センサ１１２、較正用三軸センサ１１４、信号受信部１１６、入力部１１８、記憶部１２０、較正用磁場導出部１２２、換算比率導出部１２４、磁場強度換算部１２６、位置同定部１２８、変形計測部１３０、報知部１３２を包含する。移動用三軸センサ１１０、基準用三軸センサ１１２、較正用三軸センサ１１４は、いずれも装置としては図１を用いて説明した三軸センサ１０８ａ（発信機）または三軸センサ１０８ｂ（受信機）であるが、目的が異なるために、説明の便宜上名称を変えて示している。本実施形態では、移動用三軸センサ１１０を発信機、基準用三軸センサ１１２および較正用三軸センサ１１４を受信機とする。以下、各部の構成をその動作と併せて説明する。
【００３２】
図３は、変形計測システム１００が盛立構造物１０２の変形を計測する一連の動作を説明するフローチャートである。また、図４は、盛立構造物１０２の盛立途中および盛立完了後を示す図である。特に、図４（ａ）は盛立構造物１０２の盛立途中を示す図、図４（ｂ）は盛立構造物１０２の盛立完了後を示す図である。なお、図４では、理解を容易にするために移動用三軸センサ１１０、基準用三軸センサ１１２を１つとして図示しているが、実際にはこれらは複数設けられ得る。
【００３３】
変形計測システム１００の導入に当たっては、まず変形を計測する盛立構造物１０２が所定の高さまで盛り立てられる（Ｓ１５０）。本実施形態において、盛立構造物１０２とは、段階的に盛立てられて構築される構造物のことであり、フィルダムや堤防の堤体、盛土等が例示される。さらに、掘削された後に段階的に埋め戻される廃棄物処分場なども包含される。
【００３４】
次に、図４（ａ）に示すように、移動用三軸センサ１１０が初期位置に配置される（Ｓ１５２）。移動用三軸センサ１１０は、盛立構造物１０２に埋設され、完成後にその位置（変位）を検知することによって盛立構造物１０２の変形を検知するためのセンサである。
【００３５】
次に、移動用三軸センサ１１０から所定の基準距離Ｌに較正用三軸センサ１１４が配置される（Ｓ１５４）。較正用三軸センサ１１４は、次に述べる換算比率を求めるために較正（キャリブレーション）を行うためのセンサである。
【００３６】
基準距離Ｌは任意に設定することができるが、盛立完成時における移動用三軸センサ１１０と基準用三軸センサ１１２の距離（予定された距離）とすることが好ましい。このように、運用時の位置関係に近い条件で較正を行うことにより、位置測定精度を向上させることができる。また基準距離Ｌは、盛立構造物１０２の変形によって移動した移動用三軸センサ１１０と基準用三軸センサ１１２の距離（予測される距離）としてもよい。特に大変形時には、発信機の予測される位置を基準とすることで、さらに位置測定精度を向上させることができる。
【００３７】
なお、較正用三軸センサ１１４を設置する位置は、移動用三軸センサ１１０と同じ盛立高さであってもよいが、移動用三軸センサ１１０の位置が明確である限りにおいて、盛立が進んだ位置（高さ）であってもよい。また較正用三軸センサ１１４を設置する位置の土質は、移動用三軸センサ１１０と連続する同じ土質であってもよいが、異なる土質の位置であってもよい。その場合は、それぞれの土質についての導電率を境界条件として入力すればよい。
【００３８】
そして、移動用三軸センサ１１０を励磁して、較正用三軸センサ１１４に誘起電圧Ｅ１（Ｅx1、Ｅy1、Ｅz1）を発生させる。較正用三軸センサ１１４に生じた誘起電圧Ｅ１は、信号発信部２３２から発信された信号を受信する信号受信部１１６が取得する（Ｓ１５６）。
【００３９】
次に、オペレータ（使用者）が、入力部１１８から、盛立構造物１０２の導電率σを入力する（Ｓ１５８）。導電率σは、既定のデータ（設計条件、文献調査等）から設定してもよいが、盛立構造物１０２を測定した実測値であるのが好ましい。これは、現場に則した導電率を設定できるからである。また、三軸センサが配置される周辺の値を取得するのがさらに好ましい。
【００４０】
地盤の電導率を測定する方法としては古くから様々な手法が知られているが、本実施形態では一例としてウェンナ法について説明する。図５は、盛立構造物１０２の導電率の実測手法の一例を示す図である。図５に示すように、盛立構造物１０２の導電率を実測するには、その測定対象域に直線上に４本の電極１４０、１４２、１４４、１４６を距離ａの等間隔で配置し、交流電源１４８によって両端の電極１４０、１４６に交流電圧を印加する。そして、電流計１５０が計測した電流ｉ（電極１４０、１４６間を流れる電流）と、電圧計１５２が計測した電圧ｖ（電極１４２、１４４間に生じる電圧）から、以下の式（式５）より導電率が算出される。
σ＝ｉ／（２×π×ａ×ｖ） （式５）
【００４１】
図３に戻って、次に、較正用磁場導出部１２２が、移動用三軸センサ１１０から出力された磁場について、基準距離Ｌおよび導電率σを用いて第１の磁場解析を行い、較正用三軸センサ１１４の位置における磁場強度Ｈ１（Ｈx1、Ｈy1、Ｈz1）を導出する（Ｓ１６０）。
【００４２】
較正用磁場導出部１２２は、有限要素法解析（ＦＥＭ）や境界要素法解析（ＢＥＭ）などの数値解析を用いて、磁場強度Ｈ１を求めることができる。移動用三軸センサ１１０が出力する磁場強度Ｈ０は励磁コイル２１６をリング状のコイル要素でモデル化して流れる電流を入力することで理論的に求めることができるから、盛立構造物１０２の形状、材質、移動用三軸センサ１１０から較正用三軸センサ１１４までの基準距離、および入力された導電率を境界条件として設定することにより、磁場強度Ｈ１を求めることができる。このように、盛立構造物１０２の導電率σを反映した第１の磁場解析を行うことで、導電率σが大きい場合や変化する場合にも影響を適切に考慮することができる。なお、励磁コイル２１６から生じる磁場の解析においては、渦電流の影響を考慮した定常渦電流解析を行うことが好ましい。
【００４３】
第１の磁場解析における解析モデル（盛立構造物１０２の形状）は、図４（ａ）に示すように盛立途中の形状としてもよく、または図４（ｂ）に示すように盛立完了後の形状としてもよい。盛立途中の形状では、三軸センサ（移動用三軸センサ１１０、較正用三軸センサ１１４）の上側が気中になっている条件を反映して、より実際の条件に近い磁場強度を求めることができる。しかし、完了後の形状で磁場解析しても磁場強度への影響は大差なく、十分な精度で磁場強度が求められることが分かっている。さらに、完了後の形状で解析モデルを作成しておけば、後述の第２の磁場解析にもその解析モデルを流用できる利点も有する。
【００４４】
次に、換算比率導出部１２４が、磁場強度Ｈ１と誘起電圧Ｅ１とから、換算比率Ｋ（Ｋｘ＝Ｈx1／Ｅx1、Ｋｙ＝Ｈy1／Ｅy1、Ｋｚ＝Ｈz1／Ｅz1）を導出する（Ｓ１６２）。換算比率Ｋは、土質の導電率が変化しなければ、距離にかかわらず一定となる。したがって較正によって地盤の導電率σに対応した磁場強度Ｈ１と誘起電圧Ｅ１との換算比率Ｋを求めておけば、運用時に誘起電圧Ｅ２を測定した場合に、これに対応する磁場強度Ｈ２を求めることができる。
【００４５】
なお、導電率σや換算比率Ｋは、ＲＯＭ、ＲＡＭ、ＥＥＰＲＯＭ、不揮発性ＲＡＭ、フラッシュメモリ、ＨＤＤ等からなる記憶部１２０に記憶される。
【００４６】
次に、複数（多数）の移動用三軸センサ１１０が埋設された状態で、盛立構造物１０２の盛立完了まで盛立が実施される（Ｓ１６４）。なお、較正用三軸センサ１１４に関しては、Ｓ１５６で誘起電圧Ｅ１を取得した後に盛立構造物１０２から撤去してよく、基準位置に設置して基準用三軸センサ１１２として利用してもよい。校正時と本計測時とで同じセンサを用いることで、センサ特性を反映した校正ができ、さらに位置測定精度を向上させることができる。また受信機能だけでなく送信機能も備えているのであれば、較正後にそのまま埋設して移動用三軸センサ１１０として利用することもできる。
【００４７】
そして、図４（ｂ）に示すように、盛立構造物１０２の表面近傍（基準位置）に基準用三軸センサ１１２が配置される（Ｓ１６６）。基準用三軸センサ１１２は１つであってもよいが、移動用三軸センサ１１０からの磁場の到達距離と盛立構造物１０２の規模を勘案して、移動用三軸センサ１１０をグルーピングし、それぞれのグループを担当する複数の基準用三軸センサ１１２を設けてもよい。さらに、同じ移動用三軸センサ１１０に対して異なる位置に配置した２以上の基準用三軸センサ１１２から位置の同定を行うことにより、それぞれが同定した位置座標の平均処理を行うことで、測定精度をさらに向上させることができる。
【００４８】
本実施形態では、移動用三軸センサ１１０すなわち発信機を埋設し、基準用三軸センサすなわち受信機を表面近傍に配置している。よって、発信機に電源を内蔵し、リモート制御とすれば盛立構造物１０２内へのケーブル配線の設置を回避できる。一方、受信機を埋設する場合には、信号発信部２３２が誘起電圧の計測値を信号受信部１１６へと発信する必要がある。そのため、電源を内蔵するには、消費電力の点で課題があり、信号受信部１１６へと信号を伝達するケーブル配線を敷設しなければならない可能性が高い。その上、複数の発信機を埋設する場合には１つの受信機で複数の発信機からの磁場強度を測定することができるため、受信機の数が少なくて済み、データ処理も容易となるため、システム構成を簡略化することができる。
【００４９】
基準用三軸センサ１１２を設置する基準位置は、その３次元座標（絶対座標：地球座標とも言う）が明確であればよく、盛立構造物１０２の表面近傍ばかりではなく、ある程度埋設された位置であってもよい。また基準位置は、例えばＧＰＳを用いたり、路面上の印を基準として測量したりして、随時更新（確認）することが好ましい。
【００５０】
運用時においては、信号受信部１１６が、移動用三軸センサ１１０を励磁した際に基準用三軸センサ１１２に生じる誘起電圧Ｅ２（Ｅx2、Ｅy2、Ｅz2）を取得する（Ｓ１６８）。そして、磁場強度換算部１２６が、誘起電圧Ｅ２と記憶部１２０に記憶された換算比率Ｋ（Ｋｘ、Ｋｙ、Ｋｚ）を用いて、以下の式（式６）より基準用三軸センサ１１２の位置における磁場強度Ｈ２（Ｈx2、Ｈy2、Ｈz2）を求める（Ｓ１７０）
Ｈ２＝Ｋ・Ｅ２＝（Ｈ１／Ｅ１）・Ｅ２ （式６）
なお、これを、ｘ成分、ｙ成分、ｚ成分別にすると、以下の式（式７〜式９）に表される。
Ｈx2＝Ｋｘ・Ｅx2＝（Ｈx1／Ｅx1）・Ｅx2 （式７）
Ｈy2＝Ｋｙ・Ｅy2＝（Ｈy1／Ｅy1）・Ｅy2 （式８）
Ｈz2＝Ｋｚ・Ｅz2＝（Ｈz1／Ｅz1）・Ｅz2 （式９）
【００５１】
次に、位置同定部１２８が、導電率σを反映した第２の磁場解析によって、基準用三軸センサ１１２の位置における磁場強度が算出した磁場強度Ｈ２となるように、移動用三軸センサ１１０の現在位置を同定する（Ｓ１７２）。第２の磁場解析における解析モデル（盛立構造物１０２の形状）は、図４（ｂ）に示すように盛立完了後の形状とする。詳細には後述するが、本実施形態では移動用三軸センサ１１０が発信機であるため、位置同定部１２８は、逆解析的手法に基づき移動用三軸センサ１１０の現在位置を同定する。なお、移動用三軸センサ１１０が受信機である場合には、通常の磁場解析によって、その現在位置を同定することができる。そして位置同定部１２８は、基準用三軸センサ１１２の基準位置に対する移動用三軸センサ１１０の現在位置の位置関係と、基準位置の絶対座標とから、移動用三軸センサ１１０の現在位置を絶対位置として同定することができる。
【００５２】
次に、変形計測部１３０が、盛立構造物１０２の変形を計測する。具体的には、移動用三軸センサ１１０の現在位置の初期位置からの変位を計算することにより、移動用三軸センサ１１０の移動量、すなわち盛立構造物１０２のその箇所の変形量を計測することができる。このような計測によって、盛立構造物１０２の変形が所定以上であると判断された場合（Ｓ１７４のＹｅｓ）には、報知部１３２がオペレータに対して警告を発する（Ｓ１７６）。また所定未満であった場合には（Ｓ１７４のＮｏ）、引き続き盛立構造物１０２の変形の計測が継続される。
【００５３】
（移動用三軸センサ１１０の現在位置の同定）
以下、位置同定部１２８による移動用三軸センサ１１０の現在位置の同定手法について説明する。ここで、上述したように、磁場強度換算部１２６が導出した磁場強度Ｈ２（Ｈx2、Ｈy2、Ｈz2）および基準用三軸センサ１１２の基準位置は既知である。
【００５４】
有限要素法解析（ＦＥＭ）や境界要素法解析（ＢＥＭ）などの数値解析は、空間をメッシュ状に細かく分割し、その１つの要素（区画）ごとに順に計算を行う。詳細については、後述する実施例に記載するが、このメッシュの１辺の長さを移動用三軸センサ１１０が備えるコイルの直径の１／４以上１／２以下とし、移動用三軸センサ１１０の初期位置から少なくとも２０ｍの範囲内の空間を他の範囲よりも細かく分割するとよい。これにより、移動用三軸センサ１１０の位置の同定精度を実用的なレベルまで向上できる。
【００５５】
現在位置を同定する１つ目の手法としては、繰り返し計算によって正解に漸近する方法が挙げられる。すなわち、盛立構造物１０２のモデルのなかで移動用三軸センサ１１０を仮想発信点に設定し、仮想発信点の位置を変えながら基準用三軸センサ１１２の基準位置における磁場強度Ｈ’（Ｈx’、Ｈy’、Ｈz’）を算出し、算出した磁場強度Ｈ’と先に導出された磁場強度Ｈ2（Ｈx2、Ｈy2、Ｈz2）との残差が最小となるような発信機（移動用三軸センサ１１０）の位置を現在位置を同定する。残差は、次の式（式１０）によって表される。
（Ｈx2−Ｈx’）２＋（Ｈy2−Ｈy’）２＋（Ｈz2−Ｈz’）２＝ｍｉｎ （式１０）
ここで、あらかじめ複数の仮想発信点について磁場強度Ｈ’を算出して残差が最小となる仮想発信点を選択してもよいが、残差が少なくなる方向に仮想発信点をずらしながら残差が十分に小さくなるまで計算を繰り返すことが好ましい。
【００５６】
図６は、移動用三軸センサ１１０の現在位置を同定する手法について説明する図である。以下、図６を参照しながら、現在位置を同定する２つ目の手法について説明する。
【００５７】
図６に示すように、まず移動用三軸センサ１１０が埋設された初期位置の周囲に複数の仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａを設定する（実際には立体であるため８点をとるべきであるが、図では説明の便宜上４点をとって図示および説明をする。）。この仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａは、移動用三軸センサ１１０の初期位置からの変位を想定した上で移動用三軸センサ１１０の現在位置を囲うように設定される。移動後の現在位置は未知であるが、移動範囲は想定可能である。
【００５８】
そして、それぞれの仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａ（図６の白丸）に移動用三軸センサ１１０があると仮定した場合に、導出された磁場強度Ｈ2となる複数の仮想受信点１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂ（図６の黒丸）を、導電率σの影響を考慮した第２の磁場解析によって導出する。なお、第２の磁場解析を行うにあたり、仮想発信点で囲まれた範囲内の分割幅は細かくすることが好ましく、例えば１ｍ以下とすることができる。
【００５９】
４つの仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａで囲まれた領域の形状は、４つの仮想受信点１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂで囲まれた領域の形状とは一致しない。これは、導電率σによって磁場が歪みを生じるためである。しかし、仮想発信点同士と仮想受信点同士の相対位置関係は比例関係にあると考えることができる。そこで、仮想発信点同士と仮想受信点同士によって囲まれる空間を同じ数で分割する。そして、仮想受信点１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂに対する基準用三軸センサ１１２の位置関係を、仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａに当てはめることによって、比例配分により移動用三軸センサ１１０の現在位置を同定することができる。
【００６０】
すなわち、図６では、仮想受信点１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂで囲まれた領域の右から２つ目、上から３つ目のマス目に基準用三軸センサ１１２が配置されているので、移動用三軸センサ１１０の現在位置も仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａで囲まれた領域の右から２つ目、上から３つ目のマス目であると同定することができる。故に、１つ目の手法のように繰り返し計算を行うよりも、簡潔かつ高速に移動用三軸センサ１１０の現在位置を同定することができる。
【００６１】
図７は、移動用三軸センサ１１０の現在位置を同定する手法について説明する図である。以下、図７を参照しながら、現在位置を同定する３つ目の手法について説明する。
【００６２】
図７に示すように、まず、移動用三軸センサ１１０が埋設された初期位置の周囲に複数の仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａを設定する。上記と同様に、移動用三軸センサの現在位置を囲うように仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａは設定される。そして、それぞれの仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａに移動用三軸センサ１１０があると仮定した場合に、基準用三軸センサ１１２の基準位置での磁場強度を、導電率σの影響を考慮した第２の磁場解析によって導出する。すると、仮想発信点ごとに異なる値の磁場強度が得られる（基本的には、遠いものほど弱くなる）。
【００６３】
ここで、図７に示すように、仮想発信点１８０ａからの磁場強度が１００、仮想発信点１８２ａからの磁場強度が１００、仮想発信点１８４ａからの磁場強度が８０、仮想発信点１８６ａからの磁場強度が８０、そして磁場強度換算部１２６が導出した磁場強度Ｈ２が９０であったとする。このとき、それぞれの仮想発信点１８０ａ、１８２ａ、１８４ａ、１８６ａに対応する磁場強度と導出された磁場強度Ｈ2の大小関係を定めることができるので、その大小関係に対応する比例配分により移動用三軸センサ１１０の現在位置を同定することができる。この場合においても、１つ目の手法のように繰り返し計算を行うよりも、簡潔かつ高速に移動用三軸センサ１１０の現在位置を同定することができる。
【００６４】
以上、本実施形態にかかる変形計測システム１００および変形計測方法によれば、盛立構造物１０２の導電率の影響を考慮して移動用三軸センサ１１０と基準用三軸センサ１１２の相対位置を計測することにより、導電率が大きい場合や、発信機と受信機の間で導電率が変化する（土質をまたぐ）場合であっても、盛立構造物１０２の変形を高精度に計測することができる。なお、この変形計測システム１００および変形計測方法は、地中、水中問わず適用することができ、光学計測機器では測定不可能な盛立構造物１０２の堤体内部や水中における測定も可能である。そのため、崩壊に至る前兆を確実に高精度に検出することができ、ひずみ発生段階での警報を行うことも可能になる。
【００６５】
なお、本実施形態では、移動用三軸センサ１１０を発信機、基準用三軸センサ１１２、および較正用三軸センサを受信機として説明した。しかし、移動用三軸センサ１１０を受信機、基準用三軸センサ１１２および較正用三軸センサを発信機としてもよい。このような場合には、埋設される移動用三軸センサ１１０の現在位置を逆解析的手法により求める必要がないので、数値解析の回数を削減でき、また同定精度のさらなる向上を図り得る。
【００６６】
（磁場解析方法の詳細）
次に、磁場解析方法の詳細について説明する。
【００６７】
（１）Ｍａｘｗｅｌｌの方程式
磁場解析における電磁気の現象を記述する構成則として、Maxwellの方程式を採用した。この方程式は、次の４式から構成される。
【００６８】
【数１】
【００６９】
また、物質の構成式は、次のように表される。
【００７０】
【数２】
【００７１】
（２）解析手法の検討
送信アンテナから発せられる交流の周波数は約１０ｋＨｚであるとすると、電磁波の進行速度は概ね光速（約３０万ｋｍ／ｓ）であることから、伝搬する電磁波の１波長は約３０ｋｍとなる。これに対して、計測対象範囲は送受信間距離が高々数十ｍであることから、電磁波の一波長に比べると非常に短い区間が対象となる。このため、送信アンテナから発せられる磁場の振幅と受信側のアンテナで計測される磁場の振幅の低下(振動としての振幅の減少)を捉えることが容易である。このような視点から見た電磁波は、波動として捉える動的問題から、ポテンシャルが低下する静的問題に置き換えることが可能である。
【００７２】
電磁波の問題を動的問題から静的問題へ転換した場合、電場の変化に伴う変位電流（過渡的な電流）を考慮しないことから、式１１〜式１４に示したＭａｘｗｅｌｌの方程式の内、式１１から式１３を連立して解くこととなる。また、同様の理由により、式１２の右辺第三項は恒等的に０として取り扱うこととなる。所謂、定常渦電流解析を実施することとなる。
【００７３】
（３）定常渦電流解析による導電率σを考慮した真の磁場強さＨの算出方法
磁場が時間的に変化すると電磁誘導により電界が生じる。電界の中に導体が置かれていると、その中には渦電流が発生する、と同時に渦電流により磁場が発生する。本解析では、この渦電流の影響を考慮した。
【００７４】
任意のベクトルをＡと書き、Ａの回転を取った式１６のベクトル▽×Ａを考えると、任意のベクトルに回転を取ったベクトル▽×Ａは、常に式１３を満たすことが分かっている。
【００７５】
【数３】
【００７６】
電磁気学の世界では、このベクトルＡをベクトルポテンシャルと呼ぶ。▽×Ａが式１３を満足することから、逆に、式１７に示すように、これを磁束密度Ｂとおくことができる。
Ｂ＝▽×Ａ （式１７）
【００７７】
電界の強さとベクトルポテンシャルＡ、電位φ（一般に言われる電圧）との間には、式１８に示す関係があることが分かっている。
【００７８】
【数４】
【００７９】
式１７と式１８を式１２に代入すると、定常渦電流の基礎式１９を得る。
【００８０】
【数５】
【００８１】
式１７、式１８を仮定することにより、式１１は数学的に自動的に満足されることから、具体的には、式１９を解くことによりＭａｘｗｅｌｌの方程式の解に帰着することとなる。
【００８２】
式１９はφとＡを未知数とする方程式であり、解を得るためには、Ａ、φに関する他の関係式が必要となる。これについては、式２０に示す電流に対するガウスの法則（定常電流では、導体の閉曲面に出入りする電流の総和は常に０となる）を用いる。
【００８３】
【数６】
【００８４】
上記２式が渦電流解析で満足すべき方程式となり、解析的には両式を連立させて、境界を貫く磁束密度が０となる境界条件の下で、２つの未知数Ａとφを求めることとなる。μ、σ、Ｊ０は与条件であるため、Ａが得られれば、導電率σに応じた任意点の磁場の強度Ｈが、式２１として解析から求められることとなる。
【００８５】
【数７】
【００８６】
発信器は互いに直交する三方向のコイルに順に励磁するため、解析も発信器と同じ三方向のコイル（仮にコイルｘ、コイルｙ、コイルｚとする）に対して実施すれば、コイルｘによる磁場の強さＨｘ、コイルｙによる磁場の強さＨｙ、コイルｚによる磁場の強さＨｚがそれぞれ求められることになる。
【実施例】
【００８７】
（コイル直径Ｄとメッシュの１辺の長さΔｘとの比に起因する解析精度）
数値解析において移動用三軸センサ１１０の現在位置を同定する際のメッシュの大きさとコイル寸法の関係を把握するため、コイル近傍のコイル理論値と解析結果を比較した。
【００８８】
表１に示すように、コイル直径Ｄ＝１．０ｍ（半径ｒ＝０．５ｍ）として、Ｄ／Δｘ（メッシュの１辺の長さ）＝１．０〜４．２となるようにΔｘを変化させた。解析モデルとしては、一辺５ｍの立方体を用い、コイル中心位置を解析領域の中心とした。コイル電流は１０００Ａとした。
【００８９】
表１に示すように、Ｄ／Δｘが大きくなるほど誤差が小さくなり、磁場強度の解析値（表１の値は解析モデル中の代表点）が理論値に近づくことが分かる。また、Ｄ／Δｘが２以上であれば十分な解析精度は得られると判断し得る。Ｄ／Δｘが４以下であれば、実用的な解析時間に収まることから、Ｄ／Δｘは２〜４とするのが好ましいと判断される。なお、ここでの誤差は以下の式（式２２）より算出される。
誤差＝｛解析値−理論値｝／理論値 （式２２）
【００９０】
【表１】
【００９１】
（メッシュの細分化領域）
メッシュの細分化領域（１〜２ｍ間隔でメッシュを細かく分割する領域）の大きさと解析精度との関係を把握するため、細分化領域の大きさを変化させて解析結果と理論値とを比較した。
【００９２】
図８は、実施例および比較例の解析モデルを示すである。図８に示すように、実施例では、送信コイルから１０ｍまでは１ｍ間隔でメッシュを細分割し、送信コイルから１０ｍ〜２０ｍ間を２ｍ間隔でメッシュを細分割し、それ以降は分割領域を増加させた。比較例では送信コイルから１０ｍまでは１ｍ間隔でメッシュを細分割したが、それ以降は分割間隔を増加させた。コイル直径Ｄ＝３．０ｍとし、コイル軸と解析モデルのｘ軸の角度を４５°とした。
【００９３】
図９は、実施例および比較例における発信点からの距離と磁場強度の解析誤差の相関を示す図である。図９に示すように、実施例では発信点からの距離に関わらず、磁場強度の解析誤差を殆ど０％とすることができる。しかし、比較例では、メッシュ分割が粗くなる発信点から１０ｍを境界として解析誤差が急増しており、その誤差は６０％程度にも達し得る。よって、解析精度を維持するためには、発信点（励磁コイル２１６）の位置から少なくとも２０ｍの範囲は空間を細かく分割する必要があると判断できる。
【００９４】
以上、本発明の好適な実施形態、実施例について詳細に説明したが、本発明は係る実施形態、実施例に限定されないことは言うまでもない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【産業上の利用可能性】
【００９５】
本発明は、２以上の三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測システムおよび変形計測方法に利用することができる。
【符号の説明】
【００９６】
１００…変形計測システム、１０２…盛立構造物、１０８ａ…三軸センサ（発信機）、１０８ｂ…三軸センサ（受信機）、１１０…移動用三軸センサ、１１２…基準用三軸センサ、１１４…較正用三軸センサ、１１６…信号受信部、１１８…入力部、１２０…記憶部、１２２…較正用磁場導出部、１２４…換算比率導出部、１２６…磁場強度換算部、１２８…位置同定部、１３０…変形計測部、１３２…報知部、１４０、１４２、１４４、１４６…電極、１４８…交流電源、１５０…電流計、１５２…電圧計、１８０ａ、１８２ａ、１８４ａ、１８６ａ…仮想発信点、１８０ｂ、１８２ｂ、１８４ｂ、１８６ｂ…仮想受信点、２１２…発信回路、２１４…増幅器、２１６…励磁コイル、２１８…切替回路、２２２…検出コイル、２２４…切替回路、２２６…検波回路、２２８…電圧計、２３０…誘起電圧演算部、２３２…信号発信部
【特許請求の範囲】
【請求項１】
軸方向が異なる３つのコイルを備え、前記コイルを励磁することによって磁場を生じ、または前記コイルが磁場を受けて発生した誘起電圧を検出して信号として出力する三軸センサを用いて、２以上の前記三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測システムであって、
盛立構造物に埋設される移動用三軸センサと、
前記盛立構造物の表面近傍に配置される基準用三軸センサと、
前記移動用三軸センサから所定の基準距離に配置される較正用三軸センサと、
前記移動用三軸センサまたは前記基準用三軸センサ、もしくは前記較正用三軸センサから出力された信号を受信する信号受信部と、
前記盛立構造物の導電率を入力する入力部と、
前記移動用三軸センサから出力された磁場について前記基準距離および前記導電率を用いて第１の磁場解析を行い、前記較正用三軸センサの位置における磁場強度Ｈ１を導出する較正用磁場導出部と、
前記較正用三軸センサの位置における磁場強度Ｈ１と該較正用三軸センサが検出した誘起電圧Ｅ１とから、前記盛立構造物において磁場強度と誘起電圧とを換算する換算比率を導出する換算比率導出部と、
前記移動用三軸センサまたは基準用三軸センサの一方が検出した誘起電圧Ｅ２と前記換算比率とから、前記誘起電圧Ｅ２を検出した三軸センサの位置における磁場強度Ｈ２を求める磁場強度換算部と、
前記導電率を反映した第２の磁場解析を行い、前記移動用三軸センサまたは基準用三軸センサのうち、前記誘起電圧Ｅ２を検出した一方の三軸センサの位置における磁場強度が前記磁場強度換算部が求めた磁場強度Ｈ２となるように、他方の三軸センサの現在位置を同定する位置同定部と、
を備えることを特徴とする変形計測システム。
【請求項２】
前記導電率は、前記盛立構造物を測定した実測値であることを特徴とする請求項１に記載の変形計測システム。
【請求項３】
前記移動用三軸センサは発信機、前記基準用三軸センサは受信機であって、
前記位置同定部は、
前記移動用三軸センサが埋設された初期位置の周囲に複数の仮想発信点を設定し、
それぞれの前記仮想発信点に前記移動用三軸センサがあると仮定した場合に、該移動用三軸センサから出力された磁場について前記第２の磁場解析を行うことにより前記磁場強度換算部が求めた磁場強度Ｈ２となる複数の仮想受信点を求め、
前記複数の仮想受信点と前記基準用三軸センサの位置関係と、前記複数の仮想発信点の設定位置とから、比例配分により前記移動用三軸センサの現在位置を同定することを特徴とする請求項１に記載の変形計測システム。
【請求項４】
前記移動用三軸センサは発信機、前記基準用三軸センサは受信機であって、
前記位置同定部は、
前記移動用三軸センサが埋設された初期位置の周囲に複数の仮想発信点を設定し、
それぞれの前記仮想発信点に前記移動用三軸センサがあると仮定した場合に、該移動用三軸センサから出力された磁場について前記第２の磁場解析を行うことにより前記基準用三軸センサの位置における複数の磁場強度を求め、
前記複数の磁場強度および前記磁場強度換算部が求めた磁場強度Ｈ２との大小関係と、前記複数の仮想発信点の設定位置とから、比例配分により前記移動用三軸センサの現在位置を同定することを特徴とする請求項１に記載の変形計測システム。
【請求項５】
前記位置同定部は、空間をメッシュ状に分割する３次元解析によって第２の磁場解析を行い、
前記メッシュの１辺の長さは、前記移動用三軸センサが備えるコイルの直径の１／４以上１／２以下に設定するとともに、前記初期位置から少なくとも２０ｍの範囲内の空間は他の範囲よりも細かく分割することを特徴とする請求項３または請求項４に記載の変形計測システム。
【請求項６】
前記盛立構造物は、土砂や岩石を盛り立てて構成されるフィルダムの堤体であることを特徴とする請求項１から５のいずれか１項に記載の変形計測システム。
【請求項７】
軸方向が異なる３つのコイルを備え、前記コイルを励磁することによって磁場を生じ、または前記コイルが磁場を受けて発生した誘起電圧を検出して信号として出力する三軸センサを用いて、２以上の前記三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測方法であって、
盛立構造物の初期位置に移動用三軸センサを埋設し、
前記盛立構造物の導電率を測定し、
前記移動用三軸センサから所定の基準距離の位置に較正用三軸センサを配置し、
前記較正用三軸センサを用いて誘起電圧Ｅ１を測定し、
前記移動用三軸センサから出力される磁場について前記基準距離および前記導電率を用いて第１の磁場解析を行い、前記較正用三軸センサの位置における磁場強度Ｈ１を導出し、
前記較正用三軸センサの位置における磁場強度Ｈ１と該較正用三軸センサが検出した誘起電圧Ｅ１とから、前記盛立構造物において磁場強度と誘起電圧とを換算する換算比率を導出し、
前記盛立構造物の表面近傍に基準用三軸センサを配置し、
前記移動用三軸センサまたは基準用三軸センサの一方によって誘起電圧Ｅ２を測定し、
前記誘起電圧Ｅ２と前記換算比率とから、前記誘起電圧Ｅ２を検出した三軸センサの位置における磁場強度Ｈ２を導出し、
前記導電率を反映した第２の磁場解析を行い、前記移動用三軸センサまたは基準用三軸センサのうち、前記誘起電圧Ｅ２を検出した一方の三軸センサの位置における磁場強度が前記磁場強度換算部が求めた磁場強度Ｈ２となるように、他方の三軸センサの位置を同定することを特徴とする変形計測方法。
【請求項１】
軸方向が異なる３つのコイルを備え、前記コイルを励磁することによって磁場を生じ、または前記コイルが磁場を受けて発生した誘起電圧を検出して信号として出力する三軸センサを用いて、２以上の前記三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測システムであって、
盛立構造物に埋設される移動用三軸センサと、
前記盛立構造物の表面近傍に配置される基準用三軸センサと、
前記移動用三軸センサから所定の基準距離に配置される較正用三軸センサと、
前記移動用三軸センサまたは前記基準用三軸センサ、もしくは前記較正用三軸センサから出力された信号を受信する信号受信部と、
前記盛立構造物の導電率を入力する入力部と、
前記移動用三軸センサから出力された磁場について前記基準距離および前記導電率を用いて第１の磁場解析を行い、前記較正用三軸センサの位置における磁場強度Ｈ１を導出する較正用磁場導出部と、
前記較正用三軸センサの位置における磁場強度Ｈ１と該較正用三軸センサが検出した誘起電圧Ｅ１とから、前記盛立構造物において磁場強度と誘起電圧とを換算する換算比率を導出する換算比率導出部と、
前記移動用三軸センサまたは基準用三軸センサの一方が検出した誘起電圧Ｅ２と前記換算比率とから、前記誘起電圧Ｅ２を検出した三軸センサの位置における磁場強度Ｈ２を求める磁場強度換算部と、
前記導電率を反映した第２の磁場解析を行い、前記移動用三軸センサまたは基準用三軸センサのうち、前記誘起電圧Ｅ２を検出した一方の三軸センサの位置における磁場強度が前記磁場強度換算部が求めた磁場強度Ｈ２となるように、他方の三軸センサの現在位置を同定する位置同定部と、
を備えることを特徴とする変形計測システム。
【請求項２】
前記導電率は、前記盛立構造物を測定した実測値であることを特徴とする請求項１に記載の変形計測システム。
【請求項３】
前記移動用三軸センサは発信機、前記基準用三軸センサは受信機であって、
前記位置同定部は、
前記移動用三軸センサが埋設された初期位置の周囲に複数の仮想発信点を設定し、
それぞれの前記仮想発信点に前記移動用三軸センサがあると仮定した場合に、該移動用三軸センサから出力された磁場について前記第２の磁場解析を行うことにより前記磁場強度換算部が求めた磁場強度Ｈ２となる複数の仮想受信点を求め、
前記複数の仮想受信点と前記基準用三軸センサの位置関係と、前記複数の仮想発信点の設定位置とから、比例配分により前記移動用三軸センサの現在位置を同定することを特徴とする請求項１に記載の変形計測システム。
【請求項４】
前記移動用三軸センサは発信機、前記基準用三軸センサは受信機であって、
前記位置同定部は、
前記移動用三軸センサが埋設された初期位置の周囲に複数の仮想発信点を設定し、
それぞれの前記仮想発信点に前記移動用三軸センサがあると仮定した場合に、該移動用三軸センサから出力された磁場について前記第２の磁場解析を行うことにより前記基準用三軸センサの位置における複数の磁場強度を求め、
前記複数の磁場強度および前記磁場強度換算部が求めた磁場強度Ｈ２との大小関係と、前記複数の仮想発信点の設定位置とから、比例配分により前記移動用三軸センサの現在位置を同定することを特徴とする請求項１に記載の変形計測システム。
【請求項５】
前記位置同定部は、空間をメッシュ状に分割する３次元解析によって第２の磁場解析を行い、
前記メッシュの１辺の長さは、前記移動用三軸センサが備えるコイルの直径の１／４以上１／２以下に設定するとともに、前記初期位置から少なくとも２０ｍの範囲内の空間は他の範囲よりも細かく分割することを特徴とする請求項３または請求項４に記載の変形計測システム。
【請求項６】
前記盛立構造物は、土砂や岩石を盛り立てて構成されるフィルダムの堤体であることを特徴とする請求項１から５のいずれか１項に記載の変形計測システム。
【請求項７】
軸方向が異なる３つのコイルを備え、前記コイルを励磁することによって磁場を生じ、または前記コイルが磁場を受けて発生した誘起電圧を検出して信号として出力する三軸センサを用いて、２以上の前記三軸センサの相対位置を計測することによって盛立構造物の変形を計測する変形計測方法であって、
盛立構造物の初期位置に移動用三軸センサを埋設し、
前記盛立構造物の導電率を測定し、
前記移動用三軸センサから所定の基準距離の位置に較正用三軸センサを配置し、
前記較正用三軸センサを用いて誘起電圧Ｅ１を測定し、
前記移動用三軸センサから出力される磁場について前記基準距離および前記導電率を用いて第１の磁場解析を行い、前記較正用三軸センサの位置における磁場強度Ｈ１を導出し、
前記較正用三軸センサの位置における磁場強度Ｈ１と該較正用三軸センサが検出した誘起電圧Ｅ１とから、前記盛立構造物において磁場強度と誘起電圧とを換算する換算比率を導出し、
前記盛立構造物の表面近傍に基準用三軸センサを配置し、
前記移動用三軸センサまたは基準用三軸センサの一方によって誘起電圧Ｅ２を測定し、
前記誘起電圧Ｅ２と前記換算比率とから、前記誘起電圧Ｅ２を検出した三軸センサの位置における磁場強度Ｈ２を導出し、
前記導電率を反映した第２の磁場解析を行い、前記移動用三軸センサまたは基準用三軸センサのうち、前記誘起電圧Ｅ２を検出した一方の三軸センサの位置における磁場強度が前記磁場強度換算部が求めた磁場強度Ｈ２となるように、他方の三軸センサの位置を同定することを特徴とする変形計測方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【公開番号】特開２０１１−３９０１０（Ｐ２０１１−３９０１０Ａ）
【公開日】平成２３年２月２４日（２０１１．２．２４）
【国際特許分類】
【出願番号】特願２００９−１８９５６３（Ｐ２００９−１８９５６３）
【出願日】平成２１年８月１８日（２００９．８．１８）
【出願人】（３０１０３１３９２）独立行政法人土木研究所 (107)
【出願人】（０００００３６８７）東京電力株式会社 (2,580)
【出願人】（３９００２７１７７）坂田電機株式会社 (16)
【Ｆターム（参考）】
【公開日】平成２３年２月２４日（２０１１．２．２４）
【国際特許分類】
【出願日】平成２１年８月１８日（２００９．８．１８）
【出願人】（３０１０３１３９２）独立行政法人土木研究所 (107)
【出願人】（０００００３６８７）東京電力株式会社 (2,580)
【出願人】（３９００２７１７７）坂田電機株式会社 (16)
【Ｆターム（参考）】
[ Back to top ]