JP2011220911A - 空中電磁探査法による地盤の解析方法 - Google Patents

空中電磁探査法による地盤の解析方法 Download PDF

Info

Publication number
JP2011220911A
JP2011220911A JP2010091934A JP2010091934A JP2011220911A JP 2011220911 A JP2011220911 A JP 2011220911A JP 2010091934 A JP2010091934 A JP 2010091934A JP 2010091934 A JP2010091934 A JP 2010091934A JP 2011220911 A JP2011220911 A JP 2011220911A
Authority
JP
Japan
Prior art keywords
data
grid
resistivity
specific resistance
flight
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.)
Granted
Application number
JP2010091934A
Other languages
English (en)
Other versions
JP5243476B2 (ja
Inventor
Katsushi Kawabe
克志 河戸
Minoru Okumura
稔 奥村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Eng Consultants Co Ltd
Nippon Engineering Consultants Co Ltd
Original Assignee
Nippon Eng Consultants Co Ltd
Nippon Engineering Consultants Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Eng Consultants Co Ltd, Nippon Engineering Consultants Co Ltd filed Critical Nippon Eng Consultants Co Ltd
Priority to JP2010091934A priority Critical patent/JP5243476B2/ja
Publication of JP2011220911A publication Critical patent/JP2011220911A/ja
Application granted granted Critical
Publication of JP5243476B2 publication Critical patent/JP5243476B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

【課題】色分け区分の作業を人為的に行うことなく、地質及び地下水、粘土鉱物の含有状況を可視化して、高い精度で評価を行うことができる地盤の解析方法を提供する。
【解決手段】空中電磁法による電磁探査を行って取得した測定データに対し各種の補正、レベリングを行い、各測定点における比抵抗値を計算し、補間処理を行って、周波数毎にグリッド形式の比抵抗データを生成し、これを数値標高モデルのデータと組み合わせることにより、３次元比抵抗モデルを作成し、任意の鉛直断面についてのグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データに対し、比抵抗差分解析を行い、更に、ラプラシアン解析を行って、比抵抗値を評価値に変換して出力する。
【選択図】図１

Description

本発明は、電磁誘導作用を利用して空中から探査を行い、１次磁場に対する２次磁場の強さの割合を測定して地盤の比抵抗値を求め、そのデータを解析することによって、調査対象地盤の地質及び地下水、粘土鉱物の含有状況を高い精度で評価することができる地盤の解析方法、及び、その解析のためのコンピュータプログラムに関する。
地山にトンネルを構築しようとする際の事前調査、或いは、斜面の防災維持管理のための調査（例えば、地すべり、土石流などの斜面災害や道路に隣接した斜面で斜面災害の原因となる地質・地下水の調査など）を行う場合、調査対象地盤の地質及び地下水や粘土鉱物の含有状況を把握する必要があり、そのための有効な方法の一つとして、ＨＥＭ（Helicopter-borne ElectroMagnetic method、ヘリコプタを用いた空中電磁法）と呼ばれる電磁探査法が知られている。
ＨＥＭは、図１に示すように、ヘリコプタ１の機体下方に電磁探査機器２を吊り下げ、この状態で調査対象域の上空を飛行させて探査を実行し、地盤の比抵抗（単位断面積を通る電流に対する単位長さ当たりの電気抵抗のことで、断面積Ｓ（ｍ２）、長さＬ（ｍ）の均質な導体の電気抵抗Ｒ（Ω）は、Ｒ＝ρＬ／Ｓで表され、このときの比抵抗定数ρを比抵抗と呼ぶ。）を計算するためのデータを３次元的に取得し、それらのデータに基づいて算出した地盤の比抵抗値から地質及び地下水、粘土鉱物の含有状況を把握しようとするものである。
より具体的には、電磁探査機器２に備えられている送信コイル３（ループコイル）に空中で交流電流を流し、交流磁場（１次磁場）を発生させ、この交流磁場による電磁波を調査対象地盤４中に透入させる。そうすると、電磁誘導現象により地盤中において渦電流５が誘起され、別の交流磁場（２次磁場）が発生する。この２次磁場の強さは地盤の比抵抗と負の相関があるため、電磁探査機器２に備えられている受信コイル６で、１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定し、それらの測定データ（同相成分比Ｉ、離相成分比Ｑ）を解析することにより、磁場が透入した深度までの地盤の平均的な比抵抗値を求めることができる。
尚、１次磁場の透入深度は、送信する交流電流の周波数に応じて異なる（周波数が高いと透入深度は浅くなり、周波数が低いと透入深度は深くなる）ので、複数の周波数の交流電流を同時に流すことにより、それらの周波数に応じた深度毎のデータ（比抵抗値を計算するためのデータ）を同時に取得することができる。調査対象地では、５０ｍ或いは１００ｍなどの間隔をおいて複数の測線で同じ測定を行い、それらを解析することにより、調査対象地盤の比抵抗分布を３次元的に把握することができる。
特開２００３−４３１５７
従来、ＨＥＭなどの電磁探査法によって取得した測定データに基づいて、調査対象地盤の地質及び地下水、粘土鉱物の含有状況を評価しようとする場合、測定データを解析することにより調査対象地盤の比抵抗値を３次元的に求め、各種の手法（レベリング、マッピング等）を用いて適宜データを加工し、比抵抗分布状況を図化することによって行っていた。
例えば、飛行測線間のデータを補間して比抵抗分布平面図を作成したり、調査対象域のＤＥＭ（Digital Elevation Model、数値標高モデル）と組み合わせて、仮想の３次元比抵抗モデルを作成し、更に、そこから任意の鉛直断面或いは水平断面の比抵抗分布断面図を作成し、図６に示すように、比抵抗の値に着目して分布状況の特性に応じた色分け区分を行って表示することにより、比抵抗分布状況、ひいては地質及び地下水、粘土鉱物の含有状況を視覚的に表現し、評価するという方法が行われていた。
このように、従来は、電磁探査法によって取得した測定データを加工して、比抵抗の値に着目して、その差異を表現するために人為的な色分け区分を行っていたため、資料を作成する作業者或いは作成された資料に基づいて判断を行う担当者において相応のキャリアや知識レベルが要求されることになり、また、それによって得られる評価も相対的なものに止まり、高い精度で評価を行うことが難しいという問題がある。
本発明は、比抵抗値に着目するだけでなく、隣接するグリッド間での比抵抗値の変化の割合に着目することにより、色分け区分の作業を人為的に行うことなく、地質及び地下水、粘土鉱物の含有状況を可視化して、高い精度で評価を行うことができる地盤の解析方法を提供することを目的とする。
本発明の請求項１に係る地盤の解析方法は、ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データに分割し、それぞれレベリングを行うステップ３と、ステップ３においてレベリングを行った飛行測線毎の測定データから、各測定点における比抵抗の値を計算し、補間処理を行い、周波数毎にグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データと、ステップ３においてレベリングを行った飛行測線毎の測定データとからなる入力格子データを生成し、この入力格子データに対しグリッドレベリングを行って、２次元的にエラーを除外してレベリング済み格子データを求め、このレベリング済み格子データから、飛行測線毎の比抵抗データを生成するステップ４と、ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成するステップ５と、ステップ５において生成した周波数毎のグリッド形式の比抵抗データと、数値標高モデルのデータとを組み合わせることにより、３次元比抵抗モデルを作成するステップ６と、ステップ６において作成した３次元比抵抗モデルから、任意の鉛直断面についてのグリッド形式の比抵抗データを生成するステップ７と、ステップ７において生成した任意の鉛直断面についてのグリッド形式の比抵抗データに対し、深度方向の累積誤差を取り除くために比抵抗差分解析を行うステップ８と、ステップ８において比抵抗差分解析を行ったグリッド形式の比抵抗データに対し、スムージング処理を行い、このスムージング処理を行ったグリッド形式の比抵抗データに対し、ラプラシアンフィルターを適用することにより、当該グリッド形式の比抵抗データを構成するすべての格子点の比抵抗値を、評価値に変換して出力するステップ９とを順次実行することにより、調査対象地盤において着目したい任意の鉛直断面についての比抵抗構造を可視化して出力することを特徴としている。
本発明の請求項２に係る地盤の解析方法は、ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、ステップ２において全体の補正を行った測定データについてレベリングを行うステップ３と、ステップ３においてレベリングを行った測定データから、飛行測線直下の鉛直断面についてのグリッド形式の比抵抗データを生成するステップ４と、ステップ４において生成した飛行測線直下の鉛直断面についてのグリッド形式の比抵抗データに対し、深度方向の累積誤差を取り除くために比抵抗差分解析を行うステップ５と、ステップ５において比抵抗差分解析を行ったグリッド形式の比抵抗データに対し、スムージング処理を行い、このスムージング処理を行ったグリッド形式の比抵抗データに対し、ラプラシアンフィルターを適用することにより、当該グリッド形式の比抵抗データを構成するすべての格子点の比抵抗値を、評価値に変換して出力するステップ６とを順次実行することにより、飛行測線直下の鉛直断面についての比抵抗構造を可視化して出力することを特徴としている。
本発明の請求項３に係る地盤の解析方法は、ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データに分割し、それぞれレベリングを行うステップ３と、ステップ３においてレベリングを行った飛行測線毎の測定データから、各測定点における比抵抗の値を計算し、補間処理を行い、周波数毎にグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データと、前記ステップ３においてレベリングを行った飛行測線毎の測定データとからなる入力格子データを生成し、この入力格子データに対しグリッドレベリングを行って、２次元的にエラーを除外してレベリング済み格子データを求め、このレベリング済み格子データから、飛行測線毎の比抵抗データを生成するステップ４と、ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成するステップ５と、ステップ５において生成した周波数毎のグリッド形式の比抵抗データに対し、スムージング処理を行い、このスムージング処理を行ったグリッド形式の比抵抗データに対し、ラプラシアンフィルターを適用することにより、当該グリッド形式の比抵抗データを構成するすべての格子点の比抵抗値を、評価値に変換して出力するステップ９とを順次実行することにより、調査対象地盤の周波数毎の水平断面についての比抵抗構造を可視化して出力することを特徴としている。
本発明の請求項４に係る地盤の解析方法は、ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データに分割し、それぞれレベリングを行うステップ３と、ステップ３においてレベリングを行った飛行測線毎の測定データから、各測定点における比抵抗の値を計算し、補間処理を行い、周波数毎にグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データと、前記ステップ３においてレベリングを行った飛行測線毎の測定データとからなる入力格子データを生成し、この入力格子データに対しグリッドレベリングを行って、２次元的にエラーを除外してレベリング済み格子データを求め、このレベリング済み格子データから、飛行測線毎の比抵抗データを生成するステップ４と、ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成するステップ５と、ステップ５において生成した周波数毎のグリッド形式の比抵抗データと、数値標高モデルのデータとを組み合わせることにより、３次元比抵抗モデルを作成するステップ６と、ステップ６において作成した３次元比抵抗モデルから、等間隔の複数の仮想水平断面、或いは、仮想等深度平面についてのグリッド形式の比抵抗データを抽出し、それぞれのグリッド形式の比抵抗データについて、ラプラシアン解析を実行し、それらの結果を再構築するステップ７とを順次実行することにより、調査対象地盤において着目したい任意の水平断面、或いは、等深度平面についての比抵抗構造を可視化して出力することを特徴としている。
本発明の請求項５に係る地盤の解析方法は、ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データに分割し、それぞれレベリングを行うステップ３と、ステップ３においてレベリングを行った飛行測線毎の測定データから、各測定点における比抵抗の値を計算し、補間処理を行い、周波数毎にグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データと、前記ステップ３においてレベリングを行った飛行測線毎の測定データとからなる入力格子データを生成し、この入力格子データに対しグリッドレベリングを行って、２次元的にエラーを除外してレベリング済み格子データを求め、このレベリング済み格子データから、飛行測線毎の比抵抗データを生成するステップ４と、ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成するステップ５と、ステップ５において生成した周波数毎のグリッド形式の比抵抗データと、数値標高モデルのデータとを組み合わせることにより、３次元比抵抗モデルを作成するステップ６と、ステップ６において作成した３次元比抵抗モデルから、等間隔の複数の仮想鉛直断面、及び／又は、仮想水平断面についてのグリッド形式の比抵抗データを抽出し、それぞれのグリッド形式の比抵抗データについて、ラプラシアン解析を実行し、それらの結果を再構築するステップ７とを順次実行することにより、調査対象地盤において着目したい複数の鉛直断面、及び／又は、水平断面についての複合的な比抵抗構造を容易に把握できる鳥瞰図を出力することを特徴としている。
本発明の地盤の解析方法は、比抵抗の値そのものを評価するのではなく、比抵抗値が、その周辺の比抵抗値に対しどのような相対関係になっているかという点に着目して、評価することとしており、調査対象地盤中における比抵抗の相対変化を抽出し、強調して出力することができるので、変化が微細であったとしても、調査対象地盤中に地下水或いは断層破砕帯が存在している場合には、その存在をより明確に、より客観的に表現することができる。従って、色分け区分の作業を人為的に行うことなく、地質及び地下水、粘土鉱物の含有状況を可視化して、高い精度で評価を行うことができる。
図１は、ヘリコプタを用いた空中電磁法の説明図である。 図２は、本発明に係る地盤の解析方法において用いられるデータアクセスプログラムの説明図である。 図３は、本発明に係る地盤の解析方法において行う比抵抗差分解析の説明図である。 図４は、本発明に係る地盤の解析方法において用いられる比抵抗構造グリッドデータの可視化プログラムの説明図である。 図５は、本発明に係る地盤の解析方法において、ラプラシアン鉛直断面解析を行って出力した任意の鉛直断面についての比抵抗構造を示す断面図の一例を示す図である。 図６は、空中電磁法による測定データに基づき、従来の方法によって作成した比抵抗断面図の一例を示す図である。
以下、本発明「地盤の解析方法」の実施形態について説明する。まず、本発明の第１の実施形態に係る「地盤の解析方法」について説明する。
第１の実施形態の「地盤の解析方法」は、調査対象地盤において着目したい任意の鉛直断面についての比抵抗構造を可視化して出力するものであり、下記のステップ１〜９を順次実行することによって実施される。
ステップ１：測定
ステップ２：測定データ全体のレベリング
ステップ３：飛行測線毎の測定データのレベリング
ステップ４：グリッディング及びグリッドレベリング
ステップ５：周波数毎のマッピング
ステップ６：ＤＥＭを利用した３次元比抵抗モデルの作成
ステップ７：任意の比抵抗鉛直断面図の作成
ステップ８：比抵抗差分解析
ステップ９：ラプラシアン鉛直断面解析
以下、上記各ステップについてそれぞれ詳細に説明する。
（ステップ１：測定）
ＨＥＭ（ヘリコプタを用いた空中電磁法）による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定する。より具体的には、調査対象域に対して予め設定した飛行ルート（調査対象域を、複数本の平行な直線によって均一な幅（例えば、５０ｍ）の細長い帯状の領域に分割し、それらの分割された各領域の中央線を順番に縦断していくようなルート）に沿ってヘリコプタを飛行させながら連続的に電磁波の送受信を行い、予め設定した時間間隔（例えば、０．１秒間隔）をおいて測定データ（同相成分比Ｉ、離相成分比Ｑ）を取得する。また、測定中に実際にヘリコプタが飛行したルート（飛行測線）、及び、飛行測線上の測定点（送受信を行った地点）の位置データ（緯度、経度情報）及び高度を、ＧＰＳ及び高度計により記録する。尚、電磁波の送受信は、本実施形態においては、異なる六つの周波数（３４０Ｈｚ、１，５００Ｈｚ、３，３００Ｈｚ、６，９００Ｈｚ、３１，０００Ｈｚ、１４０，０００Ｈｚ）について同時に行う。
測定には、低ノイズのデジタルセンサーを搭載し、送信コイルの近くに小さなバッキングコイル（受信コイルと逆巻き直列接続）を配置した電磁探査機器を用いる。地盤中において発生する２次磁場よりも、送信コイルへの通電によって発生する１次磁場の方がはるかに大きいため、バッキングコイルを配置することにより、受信コイルに及ぼされる１次磁場の影響を相殺することができ、２次磁場を好適に測定することができる。
電磁探査による測定を行う場合、測定前に、電磁探査機器について空中キャリブレーションを行う。電磁探査機器のキャリブレーションは、従来、地上で行われていたが、地盤の影響を受けてしまうため不可避的な誤差が生じていた。本実施形態においては、地盤の影響を受けない高高度（例えば、高度１５００フィート以上）でキャリブレーションを行うことにより、誤差を解消することができ、測定データを高い精度で取得することができる。
また、測定の途中、電磁探査機器をリセットして、後に行う補正（レベリング）のための基準ポイントを作成する。空中電磁法においては、測定値にある程度の「ずれ」が生じてしまうことがある。例えば、測定の途中や終了時に、地盤の影響を受けない高高度（地盤中において発生させた２次磁場による電磁波を受信できない高度）に上昇した場合、電磁波は記録されないはずであるところ（つまり、測定値は０レベルとなるべきところ）、あるレベルで電磁波が記録されてしまうという問題がある（ドリフト）。この問題を解決して精度の高い測定データを得るためには、測定データを構成する連続した実測値のそれぞれに含まれる「ずれ」の大きさを把握し、それらを加減して測定データを補正する必要がある。
但し、２次磁場による電磁波の受信中（２次磁場の電磁波を受信できる高度を飛行中）は、「ずれ」の大きさを把握することができない。一方、地盤中において発生させた２次磁場による電磁波を受信できない高高度に上昇した場合には、測定値は０レベルとなるはずであるため、「ずれ」が生じているとすれば、その大きさを把握することができる。従って、測定データを構成する連続した測定値の中に、本来０レベルとなるべき基準ポイントを作成しておけば、当該基準ポイントにおける「ずれ」の大きさの値に基づいて、測定データ全体を補正することができる。
本実施形態においては、測定の開始前、測定の途中（例えば、１５〜２０分に１回）と測定終了後に、２次磁場による電磁波の測定値が本来０レベルとなるべき高高度に上昇し、電磁探査機器をリセットするとともにその位置（ＧＰＳで取得した位置データ）を記録して基準ポイントを作成し、当該基準ポイントにおける「ずれ」の大きさを把握することにより、後に行う補正に備えるようになっている。尚、基準ポイント作成のための高高度飛行を行う間、調査対象域の内外を問わず測定を継続的に行い、測定データを連続取得する。
（ステップ２：測定データ全体のレベリング）
ステップ１の電磁探査によって取得した未加工の測定データは相対値であり、しかも、上述したようなドリフトによる「ずれ」を含んでいるため、測定中に高高度で電磁探査機器をリセットすることによって作成した基準ポイント（及び当該基準ポイントにおける「ずれ」の大きさの値）を利用して、測定データ全体の補正（レベリング）を行う。尚、この作業は、電磁探査（フライト）の終了後、地上にて行う。
（ステップ３：飛行測線毎の測定データのレベリング）
ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データ（測線データ）に分割し、測線データ毎に、それぞれレベリングを行う。具体的には、測線データのうち、周波数間の値に矛盾（ずれている点）が無いかどうか、また、飛行測線間の連続性から、エラーが無いかどうか目視によりチェックし、適宜測定データに修正量を手入力して修正する。
（ステップ４：グリッディング及びグリッドレベリング）
ステップ３における測定データのレベリングは、飛行測線単位のデータ調整（上下オフセット）であり、単純な平滑化のプロセスである。この段階では、測線単位のデータ内部には、有用な情報とエラーが混在している。そこで、ステップ３においてレベリングを行った測定データに対し、グリッディング及びグリッドレベリングを行う。グリッドレベリングは、飛行測線単位の測定データ内部からノイズを除去する作業であり、単純な平滑化のプロセスではない。
ａ）グリッディング（入力格子データの生成）
測定中にＧＰＳによって取得した位置データを用いて、ステップ３においてレベリングを行った測定データから、調査対象域外の測定データを除外し、調査対象域内の測定データを抽出する。次に、抽出した測定データを構成する同相成分比Ｉ及び離相成分比Ｑの各値から、各測定点における比抵抗の値を計算する。尚、測定データを構成する飛行測線毎の測定値は、前述の通り、予め設定された時間間隔（０．１秒間隔）で取得されるが、飛行速度は必ずしも一定ではないため、飛行測線上の各測定点間の距離は均一にはならない。また、飛行測線についても、必ずしも直線であるとは限らず、更に、各飛行測線の間隔も均等ではない。
従って、すべての測定点を調査対象域の平面座標上に表現すると、それらは不規則に分布することになる。そこで、本実施形態においては、それらの不規則に分布する各測定点間における比抵抗の値について補間処理を行うとともに、飛行測線間における比抵抗の値についても補間処理を行い、更に、調査対象域を縦方向及び横方向へ等間隔（例えば、１０ｍ間隔）に区分する格子状の座標をプロットし、周波数毎に、格子の各交点（格子点）の座標における比抵抗の値を計算することにより、グリッド形式の比抵抗データを生成する。そして、このグリッド形式の比抵抗データと、測線データとからなる入力格子データを生成する。
ｂ）グリッドレベリング
入力格子データには、測線方向に延びる縞状構造のエラーが含まれている可能性があるので、２次元的にこのエラーを除外する。より具体的には、同じ周波数の比抵抗値は平面的に急変しないため、一定の基準値を設けて周辺の格子点における比抵抗の値との比較を行い、エラーを除外する。
ｂ−１）２次元長方形領域フィルター処理
まず、入力格子データからバックグラウンド格子データを生成する。バックグラウンドは、ステップ３においてレベリングを行った測定データの長周期の傾向を示すものであり、バックグラウンド格子データは、入力格子データ上で、測線に直交する方向を長辺に、測線方向を短辺とする長方形領域を基本単位とするメディアンフィルター処理（２次元長方形領域フィルター処理）で得られる。これは、短辺も解析グリッドを複数含むことから２次元フィルター処理となる。メディアンフィルターの処理単位の形状は、測線データ間でも共通する長周期成分（バックグラウンド）と測線データに包含される有用な短周期成分の大きさから決定する。
ｂ−２）差分処理Ａ
次に、入力格子データから、バックグラウンド格子データを差し引くこと（差分処理）により、測線データの内部に包含される長周期成分を除去する。この処理によって得られたデータを、便宜上「仮出力格子データ」と呼ぶこととする。この仮出力格子データは、入力格子データの内部に包含される長周期成分を除去したに過ぎず、未だ、有用な情報とエラーが混在している。
ｂ−３）１次元フィルター処理
ここで、エラー成分を分離するため、あえて有用な情報である短周期成分を除去する。短周期成分の除去は、仮出力格子データの測線方向を長辺、解析グリッドサイズを短辺とする処理単位（解析グリッドサイズは基本単位であることから１次元フィルターとなる）で、メディアンフィルター処理（１次元フィルター処理）を行うことで達成できる。尚、この処理の結果得られる格子データを、便宜上「エラー値格子データ」と呼ぶこととする。
ｂ−４）差分処理Ｂ
そして、入力格子データからエラー値格子データを差し引くこと（差分処理）により、有用な情報のみを反映したレベリング済み格子データを生成する。
以上のように、２次元の入力格子データから、バックグラウンド格子データを仲介して作成したエラー値格子データを除去することにより、有用な情報のみを反映したレベリング済み格子データを求めるプロセスが「グリッドレベリング」である。
そして、レベリング済み格子データから、飛行測線毎の比抵抗データを生成する。ここで生成した飛行測線毎の比抵抗データは、ステップ３においてレベリングを行うことによってエラーを除外した飛行測線毎の測定データよりも、更に、調査対象地盤の比抵抗値を正確に反映したものとなる。
（ステップ５：周波数毎のマッピング）
ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成し、周波数毎に比抵抗平面図を作成する。補間は、例えば、クリギング法によって行う。
（ステップ６：ＤＥＭを利用した３次元比抵抗モデルの作成）
空中電磁法によって取得した測定データには、各測定点の地盤標高の情報は含まれていない。そのため、調査対象地盤の３次元的な比抵抗モデルを作成するためには、地表面の数値標高モデル（ＤＥＭ）が必要となる。本実施形態においては、ステップ５において生成した周波数毎のグリッド形式の比抵抗データと、別途準備したＤＥＭデータとを組み合わせることにより、任意の３次元位置の比抵抗値を取り出せる比抵抗モデルを作成する。尚、この３次元比抵抗モデルは、２次元データのセット（ステップ５において生成した周波数毎のグリッド形式の比抵抗データ（比抵抗平面図）、及び、ＤＥＭ）と、それにアクセスするプログラム（データアクセスプログラム）によって作成する。
このデータアクセスプログラムは、図２（データアクセスプログラムの概念図）に示すように、鉛直断面の位置、標高、或いは、深度を入力することにより、２次元データのセットから、任意の鉛直断面、等標高断面、或いは、等深度断面を出力するように構成されており、このプログラムを用いることにより、空中電磁法によって取得した比抵抗値の測定データを３次元情報として活用することができる。尚、このデータアクセスプログラムは、２次元データのセットの平面座標を一元的に管理する機能と、入力データの処理に必要なデータ（セット）を任意の平面分解能で２次元データのセットから取得する機能と、ある平面座標位置に対して取得した各周波数の比抵抗値に対して表示深度を計算する機能と、ある平面座標位置に対して取得した全周波数の比抵抗値と、計算した表示深度の組合せで、任意の鉛直間隔の１次元補間を行う機能と、ある平面座標位置に対して取得した標高との組み合わせで深度と標高の表示変換を行う機能と、個別のアクセス、変換結果を統合して、グリッド形式のデータとして出力する機能とを有している。
この段階で、例えばボクセルのような３次元データ構造を作り、画一化したモデル化（例えば１０ｍ×１０ｍ×１０ｍのサイコロ化）を行うと、データの品質が落ちしてしまうという問題があるが、本実施形態においては、上述の通り、必要に応じて必要箇所（ＸＹＺ）の値を、必要な処理（補間やフィルター）を行って取り出すようにしているため、データの品質が落ちてしまうという問題を好適に回避することができる。
（ステップ７：任意の比抵抗鉛直断面図の作成）
ステップ６において作成した３次元比抵抗モデルから、調査対象地盤において着目したい任意の鉛直断面についての比抵抗分布を図化して出力する。具体的には、当該断面における地盤の輪郭形状、及び、グリッド間隔を１ｍに設定した格子座標をプロットするとともに、各格子点の座標に係るグリッド形式の比抵抗データを生成し、これを図化して比抵抗鉛直断面図を作成する。尚、ここでは、着目したい任意の鉛直断面についてのグリッド形式の比抵抗データから、比抵抗鉛直断面図を作成しているが、次のステップを実行するためには、グリッド形式の比抵抗データが得られればよく、必ずしも図化を行って比抵抗鉛直断面図を作成しなければならない訳ではない。
また、３次元比抵抗モデルは、データ密度に異方性があるが、入手できるＤＥＭがレーザープロファイラーデータの１ｍメッシュが主流であること、空中電磁法のデータが測線直下で概ね１ｍ間隔で取得することから、密な方に合わせて１ｍ間隔の格子を標準とする（測線と測線の間は飛行計画に依存する。また、深度方向の分解能はコイル数に依存する）。
（ステップ８：比抵抗差分解析）
空中電磁法による電磁探査においては、送信される電磁波の周波数が低いほど、地盤中への透入深度は深くなる。従って、低い周波数の電磁波によって取得された比抵抗データは、高い周波数の電磁波によって取得された比抵抗データよりも深い位置の地盤情報を反映していることになる。但し、送信される電磁波は地表から投入されるため、低い周波数の電磁波によって取得された比抵抗データは、より浅い地盤（より高い周波数の透入深度の地盤）の影響を受けることになり、その結果、ステップ７において作成した比抵抗鉛直断面図（当該断面図の元となる鉛直断面についてのグリッド形式の比抵抗データ）においては、深部ほど、誤差の累積が大きくなる。そこで、この深度方向の累積誤差を取り除くために、ステップ７において作成した比抵抗鉛直断面図の元となる鉛直断面についてのグリッド形式の比抵抗データに対し、比抵抗差分解析を行う。
ここで行う比抵抗差分解析は、ある周波数の深度（差分深度）と比抵抗（差分比抵抗）を、その周波数よりも１ランク高い周波数についての所定の項目値との差から求めるというものであり、より具体的には、本実施形態において使用される六つの電磁波の周波数のうち、最も高い１４０，０００Ｈｚを「周波数Ａ」とし、以降、周波数の高い順に、３１，０００Ｈｚを「周波数Ｂ」、６，９００Ｈｚを「周波数Ｃ」、３，３００Ｈｚを「周波数Ｄ」、１，５００Ｈｚを「周波数Ｅ」、３４０Ｈｚを「周波数Ｆ」としたとき、周波数Ｂの差分深度Ｂと差分比抵抗Ｂを、周波数Ａとの差から求め、周波数Ｃの差分深度Ｃと差分比抵抗Ｃを、周波数Ｂとの差から求め、更に同様に、周波数Ｄ、Ｅ、Ｆの差分深度と差分比抵抗を、それぞれ周波数Ｃ、Ｄ、Ｅとの差から求める。
例えば、周波数Ｂの差分深度Ｂと差分比抵抗Ｂを、周波数Ａとの差から求める場合、図３に示すように、まず、周波数Ａについて、見掛比抵抗Ａ（解析前の比抵抗）、及び、見掛深度Ａ（解析前の深度）を計算する。次に、下記の計算式（数１〜３）を用いて、表皮深度Ａ、有効深度Ａ、及び、コンダクタンスＡを求める。
一方、周波数Ｂについても、見掛比抵抗Ｂ、及び、見掛深度Ｂを計算し、上記計算式と同様の式により、見掛比抵抗Ｂと周波数Ｂから表皮深度Ｂを求め、見掛深度Ｂと表皮深度Ｂから有効深度Ｂを求め、見掛比抵抗Ｂと表皮深度ＢからコンダクタンスＢを求める。
次に、有効深度Ｂから有効深度Ａを減じることによって「有効深度の差」を求め（有効深度の差＝有効深度Ｂ−有効深度Ａ）、更に、コンダクタンスＢからコンダクタンスＡを減じることによって「コンダクタンスの差」を求める（コンダクタンスの差＝コンダクタンスＢ−コンダクタンスＡ）。そして、下記の計算式（数４、数５）により、差分深度Ｂ（解析後の深度）、及び、差分比抵抗（解析後の比抵抗）を求める。
上記のような、周波数Ｂについての差分深度Ｂ及び差分比抵抗Ｂの計算（比抵抗差分解析）を、同様の要領により、周波数Ｃ〜Ｆについてもそれぞれ実行する。
（ステップ９：ラプラシアン鉛直断面解析）
調査対象地盤において着目したい任意の鉛直断面について、ステップ７により生成したグリッド形式の比抵抗データを用いてコンター図を作成し、比抵抗分布を可視化すれば、その断面における比抵抗構造が把握できる。しかし、コンター図によって比抵抗分布を可視化しただけでは、探査の目的となる地盤の状況を明確に表現できない場合がある。より具体的には、調査対象地盤中における地下水或いは断層破砕帯の検出を目的として探査を行う場合、次のような問題がある。
調査対象地盤中に地下水が存在する場合には、深度方向に比抵抗の値が変化するため、理論的には、そのような比抵抗値の変化を読み取ることによって地下水の有無を検出することができると考えられる。しかしながら現実には、標準的な空中電磁法における深度方向の分解能が十分ではないため（粗いため）、コンター図では、比抵抗値の変化を明瞭に読み取ることができない場合が多い。また、地下水が飽和状態であるか、不飽和状態であるかによって比抵抗の値が異なることを利用して地下水の有無を判断しようとする場合、地質の違い（例えば、砂であるか、粘土であるか）によって、地下水の有無以上の比抵抗値の変化を示すことがあり、この場合もコンター図での判断は難しい。
また、断層破砕帯を検出しようとする場合、層厚が５〜１０ｍ程度で、縦（垂直または斜め）方向に延在する断層破砕帯は、その周辺の領域と比べて低比抵抗となる構造が想定される。このような断層破砕帯が存在している場合、比抵抗断面図において比抵抗コンターの曲線の一部が細長く下方へ突出し、垂れ下がっているような形状として検出されることが多いが、その判読は客観性に欠ける。
調査対象地盤中に地下水或いは断層破砕帯が存在する場合、その領域は、その周辺の領域と比べ、相対的に低比抵抗域となることが多い。従って、そのような比抵抗の相対変化を抽出し、強調して出力することができれば、変化が微細であったとしても、地下水或いは断層破砕帯が地盤中に存在している場合には、その存在をより明確に、より客観的に表現することができる。
ところで、画像処理の技術分野では、周辺画素データとの相対変化を抽出し、強調表示するための手法として、ラプラシアンフィルターを使用することが知られている。これは、例えばデジタルカメラで取得した画像データに対し、エッジ（色分け模様の輪郭）を抽出する処理を行う際に用いられるものであり、格子状に配列された画素データを二次微分するものである。
画像データに対し、ラプラシアンフィルターを用いてエッジの抽出処理を行うという手法を、空中電磁法の測定データに基づいて算出した比抵抗データの処理に応用しようとする場合、データ密度の異方性と、そのために生じる補間領域の処理が問題となる。この点について具体的に説明すると、デジタルカメラの場合は、データ密度が均質で、基本的にエラーを含まない格子状の画素データが得られるのに対し、空中電磁法による比抵抗データは、必要な情報以上にノイズが強調されてしまう可能性がある。従って、空中電磁法による比抵抗データの処理に際し、ラプラシアンフィルターを用いた手法を単純に適用しただけでは、調査対象地盤の特質を的確に表現することはできない。
本実施形態においては、次のような要領で、ラプラシアンフィルターを用いた比抵抗データの処理を実施する。これにより、上述のような問題を回避し、調査対象地盤の特質を的確に表現することができる。
まず、ステップ８において比抵抗差分解析を行ったグリッド形式の比抵抗データ（調査対象地盤において着目したい任意の鉛直断面のグリッド形式の比抵抗データ）に対し、スムージング処理を行う。この処理は、ステップ８において解析を行ったグリッド形式の比抵抗データに含まれるすべての格子点の比抵抗値を、その周辺において一定の距離にある複数（例えば、４４１個）の格子点の比抵抗値の平均値に置き換える、というものである。
次に、スムージング処理を行ったグリッド形式の比抵抗データに対し、ラプラシアンフィルターを適用することにより、当該グリッド形式の比抵抗データを構成するすべての格子点の比抵抗値を、評価値に変換して出力する。具体的には、グリッド形式の比抵抗データを構成する各格子点（評価対象となる格子点）の比抵抗値を、その周辺（上方、下方、右方、及び、左方）に位置する他の四つの格子点（比較対象となる格子点）の比抵抗値と比較（差を計算）して比較値（差）を出力し、その比較値（差）の総和を、評価対象となる格子点についての比抵抗値と置き換え、評価値として出力する。
このように、ある格子点（評価対象となる格子点）の比抵抗値を、その周辺の格子点（比較対象となる格子点）の比抵抗値と比較する意図乃至目的は、評価対象となる格子点の比抵抗値が、その周辺の格子点の比抵抗値に対しどのような相対関係になっているかという点に着目して、その格子点の比抵抗値を評価しようとするところにある。また、ステップ８において解析を行ったグリッド形式の比抵抗データに対し、そのままラプラシアンフィルターを適用した場合、比抵抗値の相対関係を評価しようとするうえで本質に関係のないノイズが顕著となる場合が多く、顕在化させたい構造をうまく出現させることができないという問題があり、本実施形態においてはこのような問題を回避すべく、ラプラシアンフィルターの適用前に、上述のようなスムージング処理を行うこととしている。
尚、画像処理の技術分野においてラプラシアンフィルターを用いる場合、比較対象となるのは、処理対象となる画素データの「周囲に隣接する」画素データであるのに対し、本実施形態において比較対象となるのは、評価対象となる格子点の周辺において、上方、下方、右方、及び、左方へ一定の距離Ｄをおいた位置にある四つの格子点である。尚、距離Ｄの値は、測定データの密度や、検出しようとする対象（地下水、断層破砕帯など）の深度や大きさに応じて、比較を行う際に、適正な値をその都度設定する。例えば、比較的浅部の地下水の検出を対象とする場合は５ｍ（上方、下方、右方、及び、左方へそれぞれ５ｍ）、深部のトンネル付近の断層破砕帯の検出を対象とする場合は１０ｍ（上方、下方、右方、及び、左方へそれぞれ１０ｍ）が標準となる。
評価対象となる格子点の比抵抗値と、比較対象となる四つの格子点の比抵抗値との比較（差の計算）、及び、評価値の出力は、次のようにして行う。評価対象となる格子点の比抵抗値をＭとし、比較対象となる四つの格子点の比抵抗値をＡ１〜Ａ４とした場合、Ｍと、Ａ１〜Ａ４それぞれとの差を計算すると、Ｍ−Ａ１、Ｍ−Ａ２、Ｍ−Ａ３、Ｍ−Ａ４となる。ここで、Ｍが比較対象より小さい場合には、比較値は負の値に、Ｍが比較対象より大きい場合には、比較値は正の値となる。これら四つの比較値の総和（４Ｍ−Ａ１−Ａ２−Ａ３−Ａ４）を計算し、これを評価対象となる格子点における評価値として出力する。例えば、四つの比較値がすべて負の場合は、大きな負の値として累積することになり、正負が入り乱れた場合や、比抵抗値に殆ど差が無い場合には、０に近い値となる。
このようにして、グリッド形式の比抵抗データにラプラシアンフィルターを適用することにより、すべての格子点の比抵抗値について、その周辺の格子点の比抵抗値との相対比較を行い、評価値に置き換えて出力する。そして、出力されたグリッド形式のデータ（比抵抗構造グリッドデータ）を、コンピュータプログラム（比抵抗構造グリッドデータの可視化プログラム）によって図化（可視化）し、コンピュータの画面或いは紙面に出力する。
ここで用いる比抵抗構造グリッドデータの可視化プログラムは、図４に示すように、比抵抗構造鉛直断面図や、等標高比抵抗構造断面図（水平断面図）、等深度比抵抗構造断面図等、比抵抗構造を示すグリッド形式のデータを活用するために、色調や、鳥瞰図における視点を調整して出力することにより可視化するものであり、入力値の正負、中間（ゼロ近傍の正負反転が過敏な領域）に対して独立して色調を設定できる機能、入力値の正負に対してそれぞれヒストグラムを作成し、コントラストを調整する機能、着色した結果をコンピュータのディスプレイに表示するとともに、必要に応じて紙面に印刷する機能、複数の比抵抗構造鉛直断面図（グリッド形式）（或いは、複数の比抵抗構造水平断面図）を統合し、比抵抗構造鳥瞰図を作成し、視点等を調整してコンピュータの画面に表示するとともに、必要に応じて紙面に印刷する機能を有している。
図５は、ステップ９のラプラシアン鉛直断面解析を行って出力した任意の鉛直断面についての比抵抗構造を示す断面図の例である。この図において、地表面を示す線の下側に表示されているグラデーション模様のうち、「青色表示」と示されている部分（コンピュータのディスプレイ上或いは紙面上においては青色に表示される）は、周囲より相対的に比抵抗が高い（色が濃い部分は高比抵抗の度合いが大きく、色が薄い部分は高比抵抗の度合いが小さい）ところを表し、それ以外のグラデーション模様の部分（コンピュータのディスプレイ上或いは紙面上においては赤色に表示される）は、周囲より相対的に比抵抗が低い（色が濃い部分は低比抵抗の度合いが大きく、色が薄い部分は低比抵抗の度合いが小さい）ところを表している。また、グラデーション模様が表示されていない白い部分は、比抵抗が周囲と相対的に変化していないところを表している。
本実施形態の地盤の解析方法は、図５の例に示すように、評価対象となる格子点の比抵抗の値そのものを評価するのではなく、隣接する格子点間での比抵抗値の変化の割合に着目することにより、つまり、評価対象となる格子点の比抵抗値が、その周辺の格子点の比抵抗値に対しどのような相対関係になっているかという点に着目して、その格子点の比抵抗値を評価することとしているため、色分け区分の作業を人為的に行うことなく、地質及び地下水、粘土鉱物の含有状況を可視化して、高い精度で評価を行うことができる。上述の通り、調査対象地盤中に地下水或いは断層破砕帯が存在する場合、その領域は、その周辺の領域と比べ、相対的に低比抵抗域となることが多いが、本実施形態の地盤の解析方法によれば、そのような比抵抗の相対変化を抽出し、強調して出力することができるので、変化が微細であったとしても、地下水或いは断層破砕帯が地盤中に存在している場合には、その存在をより明確に、より客観的に表現することができる。
第２の実施形態に係る「地盤の解析方法」は、飛行測線直下の鉛直断面についての比抵抗構造を可視化して出力するものであり、下記のステップ１〜６を順次実行することによって実施される。
ステップ１：測定
ステップ２：測定データ全体のレベリング
ステップ３：飛行測線の測定データのレベリング
ステップ４：飛行測線直下の比抵抗鉛直断面図の作成
ステップ５：比抵抗差分解析
ステップ６：ラプラシアン鉛直断面解析
本実施形態における地盤の解析方法は、地山にトンネルを構築しようとする際の事前調査等において、計画路線直上のルートに沿ってヘリコプタを飛行させることによって測定データを得ることができたような場合を想定している。このような場合には、第１の実施形態におけるステップ４〜６を実行せずに、ステップ３において飛行測線（計画路線直上の飛行測線）の測定データについてレベリングを行った後、ステップ４として、そのレベリングを行った測定データから、直接的に比抵抗鉛直断面図（第１の実施形態におけるステップ７において作成する比抵抗鉛直断面図に相当）を作成した方が、マッピングによる品質低下を避ける意味で有効な場合があるからである。
従って、本実施形態においては、ステップ１〜３（第１の実施形態におけるステップ１〜３と同一）を実行した後、ステップ４として、第１の実施形態におけるステップ７に相当する断面図の作成を行い、ステップ５として、比抵抗差分解析（第１の実施形態におけるステップ８と同一）を行い、ステップ６として、ラプラシアン鉛直断面解析（第１の実施形態におけるステップ９と同一）を行うこととしている。これらのステップを順次実行することにより、第１の実施形態と同様の効果を期待することができる。
第３の実施形態に係る「地盤の解析方法」は、調査対象地盤の水平断面（周波数毎）についての比抵抗構造を可視化して出力するものであり、下記のステップ１〜６を順次実行することによって実施される。
ステップ１：測定
ステップ２：測定データ全体のレベリング
ステップ３：飛行測線毎の測定データのレベリング
ステップ４：グリッディング及びグリッドレベリング
ステップ５：周波数毎のマッピング
ステップ６：周波数毎のラプラシアン平面解析
本実施形態においては、周波数毎の水平断面の比抵抗構造を解析対象としているため、第１の実施形態におけるステップ６〜８は不要である。従って、本実施形態においては、ステップ１〜５（第１の実施形態におけるステップ１〜５と同一）を実行した後、ステップ６として、第１の実施形態におけるステップ９に相当するラプラシアン解析を行うこととしている。より詳細には、ステップ６として、ステップ５において生成した周波数毎のグリッド形式の比抵抗データ（周波数毎の比抵抗平面図）に対し、ラプラシアン解析（第１の実施形態のステップ９のラプラシアン解析と同一）を実行する。これにより、調査対象地盤の水平断面（周波数毎）についての比抵抗構造を可視化して出力することができる。
第４の実施形態に係る「地盤の解析方法」は、調査対象地盤において着目したい任意の水平断面、或いは、等深度平面についての比抵抗構造を可視化して出力するものであり、下記のステップ１〜７を順次実行することによって実施される。
ステップ１：測定
ステップ２：測定データ全体のレベリング
ステップ３：飛行測線毎の測定データのレベリング
ステップ４：グリッディング及びグリッドレベリング
ステップ５：周波数毎のマッピング
ステップ６：ＤＥＭを利用した３次元比抵抗モデルの作成
ステップ７：ラプラシアン水平断面解析・等深度平面解析
本実施形態における地盤の解析方法は、着目したい任意の水平断面、或いは、等深度平面についての比抵抗構造を可視化しようとするものであるため、第１の実施形態におけるステップ７（任意の比抵抗鉛直断面図の作成）は不要である。
従って、本実施形態においては、ステップ１〜６（第１の実施形態におけるステップ１〜６と同一）を実行した後、ステップ７として、第１の実施形態におけるステップ９に相当するラプラシアン解析を行うこととしている。より詳細には、ステップ７として、ステップ６において生成した３次元比抵抗モデルから、等間隔（例えば、５ｍ間隔）の複数の仮想水平断面、或いは、仮想等深度平面についてのグリッド形式の比抵抗データを抽出し、それぞれのグリッド形式の比抵抗データについて、ラプラシアン解析（第１の実施形態のステップ９のラプラシアン解析と同一）を実行し、それらの結果を再構築する。これにより、調査対象地盤の任意の水平断面、或いは、等深度平面についての比抵抗構造を可視化して出力することができる。
尚、ステップ６において生成した３次元比抵抗モデルから、着目したい任意の比抵抗水平断面、或いは、等深度平面を抽出し、その抽出された断面、或いは、平面についてのみラプラシアン解析を行うと、精度が低下してしまう可能性があり、好ましくない。このため、本実施形態においては、上述の通り、３次元比抵抗モデルから、等間隔の複数の仮想水平断面等についての比抵抗データを抽出して、それぞれについてラプラシアン解析を実行するとともに、それらの結果を再構築することとしている。これにより、精度の低下を好適に防止することができる。
第５の実施形態に係る「地盤の解析方法」は、調査対象地盤において着目したい複数の鉛直断面、及び／又は、水平断面等についての複合的（立体的）な比抵抗構造を容易に把握できる鳥瞰図を出力するものであり、下記のステップ１〜７を順次実行することによって実施される。
ステップ１：測定
ステップ２：測定データ全体のレベリング
ステップ３：飛行測線毎の測定データのレベリング
ステップ４：グリッディング及びグリッドレベリング
ステップ５：周波数毎のマッピング
ステップ６：ＤＥＭを利用した３次元比抵抗モデルの作成
ステップ７：ラプラシアン鳥瞰解析
本実施形態においては、ステップ１〜６（第１の実施形態におけるステップ１〜６と同一）を実行した後、ステップ７として、第１の実施形態におけるステップ９に相当するラプラシアン解析を行うこととしている。より詳細には、ステップ７として、ステップ６において生成した３次元比抵抗モデルから、等間隔（例えば、５ｍ間隔）の複数の仮想鉛直断面、及び／又は、仮想水平断面についてのグリッド形式の比抵抗データを抽出し、それぞれのグリッド形式の比抵抗データについて、ラプラシアン解析（第１の実施形態のステップ９のラプラシアン解析と同一）を実行し、それらの結果を再構築する。これにより、着目したい複数の鉛直断面、及び／又は、水平断面等についての複合的、立体的な比抵抗構造を容易に把握できる鳥瞰図を出力することができる。
尚、ステップ６において生成した３次元比抵抗モデルから、着目したい複数の鉛直断面、及び／又は、水平断面を抽出し、その抽出された断面についてのみラプラシアン解析を行うと、精度が低下してしまう可能性があり、好ましくない。このため、本実施形態においては、上述の通り、３次元比抵抗モデルから、等間隔の複数の鉛直断面、及び／又は、水平断面についてのグリッド形式の比抵抗データを抽出して、それぞれについてラプラシアン解析を実行するとともに、それらの結果を再構築することとしている。これにより、精度の低下を好適に防止することができる。
１：ヘリコプタ、
２：電磁探査機器、
３：送信コイル、
４：地盤、
５：渦電流、
６：受信コイル

Claims (5)

  1. ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、
    前記ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、
    前記ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データに分割し、それぞれレベリングを行うステップ３と、
    前記ステップ３においてレベリングを行った飛行測線毎の測定データから、各測定点における比抵抗の値を計算し、補間処理を行い、周波数毎にグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データと、前記ステップ３においてレベリングを行った飛行測線毎の測定データとからなる入力格子データを生成し、この入力格子データに対しグリッドレベリングを行って、２次元的にエラーを除外してレベリング済み格子データを求め、このレベリング済み格子データから、飛行測線毎の比抵抗データを生成するステップ４と、
    前記ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成するステップ５と、
    前記ステップ５において生成した周波数毎のグリッド形式の比抵抗データと、数値標高モデルのデータとを組み合わせることにより、３次元比抵抗モデルを作成するステップ６と、
    前記ステップ６において作成した３次元比抵抗モデルから、任意の鉛直断面についてのグリッド形式の比抵抗データを生成するステップ７と、
    前記ステップ７において生成した任意の鉛直断面についてのグリッド形式の比抵抗データに対し、深度方向の累積誤差を取り除くために比抵抗差分解析を行うステップ８と、
    前記ステップ８において比抵抗差分解析を行ったグリッド形式の比抵抗データに対し、スムージング処理を行い、このスムージング処理を行ったグリッド形式の比抵抗データに対し、ラプラシアンフィルターを適用することにより、当該グリッド形式の比抵抗データを構成するすべての格子点の比抵抗値を、評価値に変換して出力するステップ９とを順次実行することにより、調査対象地盤において着目したい任意の鉛直断面についての比抵抗構造を可視化して出力することを特徴とする地盤の解析方法。
  2. ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、
    前記ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、
    前記ステップ２において全体の補正を行った測定データについてレベリングを行うステップ３と、
    前記ステップ３においてレベリングを行った測定データから、飛行測線直下の鉛直断面についてのグリッド形式の比抵抗データを生成するステップ４と、
    前記ステップ４において生成した飛行測線直下の鉛直断面についてのグリッド形式の比抵抗データに対し、深度方向の累積誤差を取り除くために比抵抗差分解析を行うステップ５と、
    前記ステップ５において比抵抗差分解析を行ったグリッド形式の比抵抗データに対し、スムージング処理を行い、このスムージング処理を行ったグリッド形式の比抵抗データに対し、ラプラシアンフィルターを適用することにより、当該グリッド形式の比抵抗データを構成するすべての格子点の比抵抗値を、評価値に変換して出力するステップ６とを順次実行することにより、飛行測線直下の鉛直断面についての比抵抗構造を可視化して出力することを特徴とする地盤の解析方法。
  3. ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、
    前記ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、
    前記ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データに分割し、それぞれレベリングを行うステップ３と、
    前記ステップ３においてレベリングを行った飛行測線毎の測定データから、各測定点における比抵抗の値を計算し、補間処理を行い、周波数毎にグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データと、前記ステップ３においてレベリングを行った飛行測線毎の測定データとからなる入力格子データを生成し、この入力格子データに対しグリッドレベリングを行って、２次元的にエラーを除外してレベリング済み格子データを求め、このレベリング済み格子データから、飛行測線毎の比抵抗データを生成するステップ４と、
    前記ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成するステップ５と、
    前記ステップ５において生成した周波数毎のグリッド形式の比抵抗データに対し、スムージング処理を行い、このスムージング処理を行ったグリッド形式の比抵抗データに対し、ラプラシアンフィルターを適用することにより、当該グリッド形式の比抵抗データを構成するすべての格子点の比抵抗値を、評価値に変換して出力するステップ９とを順次実行することにより、調査対象地盤の周波数毎の水平断面についての比抵抗構造を可視化して出力することを特徴とする地盤の解析方法。
  4. ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、
    前記ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、
    前記ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データに分割し、それぞれレベリングを行うステップ３と、
    前記ステップ３においてレベリングを行った飛行測線毎の測定データから、各測定点における比抵抗の値を計算し、補間処理を行い、周波数毎にグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データと、前記ステップ３においてレベリングを行った飛行測線毎の測定データとからなる入力格子データを生成し、この入力格子データに対しグリッドレベリングを行って、２次元的にエラーを除外してレベリング済み格子データを求め、このレベリング済み格子データから、飛行測線毎の比抵抗データを生成するステップ４と、
    前記ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成するステップ５と、
    前記ステップ５において生成した周波数毎のグリッド形式の比抵抗データと、数値標高モデルのデータとを組み合わせることにより、３次元比抵抗モデルを作成するステップ６と、
    前記ステップ６において作成した３次元比抵抗モデルから、等間隔の複数の仮想水平断面、或いは、仮想等深度平面についてのグリッド形式の比抵抗データを抽出し、それぞれのグリッド形式の比抵抗データについて、ラプラシアン解析を実行し、それらの結果を再構築するステップ７とを順次実行することにより、調査対象地盤において着目したい任意の水平断面、或いは、等深度平面についての比抵抗構造を可視化して出力することを特徴とする地盤の解析方法。
  5. ヘリコプタを用いた空中電磁法による電磁探査を行い、調査対象地盤における１次磁場に対する２次磁場の強さの割合を同相成分と離相成分とに分離して測定するステップ１と、
    前記ステップ１の電磁探査によって取得した測定データ全体について、ドリフトによるずれを補正するステップ２と、
    前記ステップ２において全体の補正を行った測定データを、飛行測線毎の測定データに分割し、それぞれレベリングを行うステップ３と、
    前記ステップ３においてレベリングを行った飛行測線毎の測定データから、各測定点における比抵抗の値を計算し、補間処理を行い、周波数毎にグリッド形式の比抵抗データを生成し、このグリッド形式の比抵抗データと、前記ステップ３においてレベリングを行った飛行測線毎の測定データとからなる入力格子データを生成し、この入力格子データに対しグリッドレベリングを行って、２次元的にエラーを除外してレベリング済み格子データを求め、このレベリング済み格子データから、飛行測線毎の比抵抗データを生成するステップ４と、
    前記ステップ４において生成した飛行測線毎の比抵抗データを基に、飛行測線間の比抵抗データを補間するとともに、再びグリッディングを行い、周波数毎にグリッド形式の比抵抗データを生成するステップ５と、
    前記ステップ５において生成した周波数毎のグリッド形式の比抵抗データと、数値標高モデルのデータとを組み合わせることにより、３次元比抵抗モデルを作成するステップ６と、
    前記ステップ６において作成した３次元比抵抗モデルから、等間隔の複数の仮想鉛直断面、及び／又は、仮想水平断面についてのグリッド形式の比抵抗データを抽出し、それぞれのグリッド形式の比抵抗データについて、ラプラシアン解析を実行し、それらの結果を再構築するステップ７とを順次実行することにより、調査対象地盤において着目したい複数の鉛直断面、及び／又は、水平断面についての複合的な比抵抗構造を容易に把握できる鳥瞰図を出力することを特徴とする地盤の解析方法。
JP2010091934A 2010-04-13 2010-04-13 空中電磁探査法による地盤の解析方法 Active JP5243476B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010091934A JP5243476B2 (ja) 2010-04-13 2010-04-13 空中電磁探査法による地盤の解析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010091934A JP5243476B2 (ja) 2010-04-13 2010-04-13 空中電磁探査法による地盤の解析方法

Publications (2)

Publication Number Publication Date
JP2011220911A true JP2011220911A (ja) 2011-11-04
JP5243476B2 JP5243476B2 (ja) 2013-07-24

Family

ID=45038072

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010091934A Active JP5243476B2 (ja) 2010-04-13 2010-04-13 空中電磁探査法による地盤の解析方法

Country Status (1)

Country Link
JP (1) JP5243476B2 (ja)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012053191A (ja) * 2010-08-31 2012-03-15 Asia Air Survey Co Ltd 地盤立体表示システム
CN106226828A (zh) * 2016-07-18 2016-12-14 吉林大学 时间域航空电磁数据约束多项式拟合调平方法
KR101979936B1 (ko) * 2018-06-29 2019-05-21 한국지질자원연구원 실시간 3d 지화학 이상대 탐지 모델 생성 시스템, 모델 생성 방법 및 이를 실행시키기 위한 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체
EP3745546A1 (en) 2011-10-05 2020-12-02 Gigaphoton Inc. Method for irradiating a target and system for the same
JP2020193933A (ja) * 2019-05-30 2020-12-03 株式会社ジオ・ヴェスト 探査方法、プログラム、プログラムを格納したコンピュータ読み取り可能な記録媒体およびコンピュータ
CN113960677A (zh) * 2021-10-26 2022-01-21 北京卫星环境工程研究所 一种快速判断目标体倾斜方向的方法
JP7346757B1 (ja) 2022-09-05 2023-09-19 中国科学院地質與地球物理研究所 複雑な山岳地帯における長大な線形工事のための工学地質ゾーニング方法及びシステム
JP7491579B2 (ja) 2021-04-06 2024-05-28 有限会社ネオサイエンス 湿潤位置の推定方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105572753B (zh) * 2015-12-11 2017-12-05 云南电网有限责任公司电力科学研究院 一种瞬变电磁探测数据等值线成图与高阻区域识别方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003043157A (ja) * 2001-07-31 2003-02-13 Hisatoshi Konishi 空中電磁探査法のドリフト補正方法
JP2003294853A (ja) * 2002-04-04 2003-10-15 Oyo Corp 空中地下構造探査方法
JP2004069384A (ja) * 2002-08-02 2004-03-04 Anritsu Corp Ｘ線異物検出装置，ｘ線異物検出方法及びｘ線異物検出プログラム
JP2009074953A (ja) * 2007-09-21 2009-04-09 Kyushu Univ 電磁探査装置、電磁探査方法及びプログラム
JP2009168746A (ja) * 2008-01-18 2009-07-30 Sumitomo Electric Ind Ltd 検査方法および検査装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003043157A (ja) * 2001-07-31 2003-02-13 Hisatoshi Konishi 空中電磁探査法のドリフト補正方法
JP2003294853A (ja) * 2002-04-04 2003-10-15 Oyo Corp 空中地下構造探査方法
JP2004069384A (ja) * 2002-08-02 2004-03-04 Anritsu Corp Ｘ線異物検出装置，ｘ線異物検出方法及びｘ線異物検出プログラム
JP2009074953A (ja) * 2007-09-21 2009-04-09 Kyushu Univ 電磁探査装置、電磁探査方法及びプログラム
JP2009168746A (ja) * 2008-01-18 2009-07-30 Sumitomo Electric Ind Ltd 検査方法および検査装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN6013010647; 千田敬二　他: '空中電磁法と応用地質分野への最近の適用展開' 日本応用地質学会北海道支部・北海道応用地質研究会研究発表会講演予稿集 号：２９, 20090626, ページ：１７-２０ *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012053191A (ja) * 2010-08-31 2012-03-15 Asia Air Survey Co Ltd 地盤立体表示システム
EP3745546A1 (en) 2011-10-05 2020-12-02 Gigaphoton Inc. Method for irradiating a target and system for the same
CN106226828A (zh) * 2016-07-18 2016-12-14 吉林大学 时间域航空电磁数据约束多项式拟合调平方法
KR101979936B1 (ko) * 2018-06-29 2019-05-21 한국지질자원연구원 실시간 3d 지화학 이상대 탐지 모델 생성 시스템, 모델 생성 방법 및 이를 실행시키기 위한 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체
JP2020193933A (ja) * 2019-05-30 2020-12-03 株式会社ジオ・ヴェスト 探査方法、プログラム、プログラムを格納したコンピュータ読み取り可能な記録媒体およびコンピュータ
JP7240574B2 (ja) 2019-05-30 2023-03-16 株式会社ジオ・ヴェスト 探査方法、プログラム、プログラムを格納したコンピュータ読み取り可能な記録媒体およびコンピュータ
JP7491579B2 (ja) 2021-04-06 2024-05-28 有限会社ネオサイエンス 湿潤位置の推定方法
CN113960677A (zh) * 2021-10-26 2022-01-21 北京卫星环境工程研究所 一种快速判断目标体倾斜方向的方法
CN113960677B (zh) * 2021-10-26 2023-09-15 北京卫星环境工程研究所 一种快速判断目标体倾斜方向的方法
JP7346757B1 (ja) 2022-09-05 2023-09-19 中国科学院地質與地球物理研究所 複雑な山岳地帯における長大な線形工事のための工学地質ゾーニング方法及びシステム

Also Published As

Publication number Publication date
JP5243476B2 (ja) 2013-07-24

Similar Documents

Publication Publication Date Title
JP5243476B2 (ja) 空中電磁探査法による地盤の解析方法
Salvini et al. The use of an unmanned aerial vehicle for fracture mapping within a marble quarry (Carrara, Italy): photogrammetry and discrete fracture network modelling
Watkins et al. Appraisal of fracture sampling methods and a new workflow to characterise heterogeneous fracture networks at outcrop
Bisdom et al. Calibrating discrete fracture-network models with a carbonate three-dimensional outcrop fracture network: Implications for naturally fractured reservoir modeling
RU2442193C2 (ru) Системы поправки на рельеф местности
Haddad et al. Applications of airborne and terrestrial laser scanning to paleoseismology
DE102015211801A1 (de) Ortungssystem mit handgeführter Ortungseinheit
CA2876632C (en) System and method for migration velocity modeling
Minisini et al. Geological data extraction from lidar 3-D photorealistic models: A case study in an organic-rich mudstone, Eagle Ford Formation, Texas
Schmidt et al. Raster was yesterday: using vector engines to process geophysical data
Wiegand et al. Regolith structure analysis—a contribution to understanding the local occurrence of shallow landslides (Austrian Tyrol)
CN105241428A (zh) 一种利用高光谱进行水深反演的方法
Deparis et al. Combined use of geophysical methods and remote techniques for characterizing the fracture network of a potentially unstable cliff site (the ‘Roche du Midi’, Vercors massif, France)
JP2003294853A (ja) 空中地下構造探査方法
Colica et al. Evaluating characteristics of an active coastal spreading area combining geophysical data with satellite, aerial, and unmanned aerial vehicles images
Tronicke et al. Vertical radar profiling: Combined analysis of traveltimes, amplitudes, and reflections
CN108761540A (zh) 一种频率域天然电场三维勘探方法
Romano et al. An Integrated Geophysical and Unmanned Aerial Systems Surveys for Multi-Sensory, Multi-Scale and Multi-Resolution Cave Detection: The Gravaglione Site (Canale di Pirro Polje, Apulia)
Holland et al. Evolution of fractures in a highly dynamic thermal, hydraulic, and mechanical system-(II) Remote sensing fracture analysis, Jabal Shams, Oman Mountains
JP5671184B1 (ja) 断層図作成装置及び断層図作成方法並びに断層図作成プログラム
Peña et al. Cover surfaces as a new technique for 3D GPR image enhancement. Archaeological applications
Kr et al. Application of Multi-Attributes and Spectral Decomposition with RGB blending for understanding the strati-structural features: A Case study
WO2017155526A1 (en) Visualization of multi-pipe inspection results
Kamiński et al. Geophysical model and displacement of active landslide—an example from Jastrzębia Góra Cliff (Northern Poland)
Ji et al. Exploiting the Potential of Urban DAS Grids: Ambient‐Noise Subsurface Imaging Using Joint Rayleigh and Love Waves

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20111201

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130111

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20130312

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130404

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160412

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 5243476

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350