JP3763400B2

JP3763400B2 - 連続電磁波型トモグラフィ及びその装置

Info

Publication number: JP3763400B2
Application number: JP2001045531A
Authority: JP
Inventors: 克戸井田; 一彦升元; 啓介日比谷; 真弓田中
Original assignee: Kajima Corp
Current assignee: Kajima Corp
Priority date: 2001-02-21
Filing date: 2001-02-21
Publication date: 2006-04-05
Anticipated expiration: 2021-02-21
Also published as: JP2002243850A

Description

【０００１】
【発明の属する技術分野】
本発明は連続電磁波型トモグラフィ及びその装置に関し、とくに周波数が段階状に切替わる多周波数連続電磁波の透過によるトモグラフィ及びその装置に関する。本発明において「トモグラフィ」とは、多方向からの電磁波の透過量の計測値から計算によって断面像、即ちトモグラムを得る方法をいう。
【０００２】
【従来の技術】
岩盤構造物の設計・施工に際しては、岩盤の構造や各種物性値を精度良く把握することが必要である。特に、高レベル放射性廃棄物の地層処分のためには、事前調査段階において数少ないボーリング孔を利用して地下深部の岩盤を評価することが望まれるが、このような場合には、使用する探査手法の探査範囲・精度を十分把握した上で調査計画を立案すべきである。
【０００３】
三次元的な岩盤の水理地質構造評価に関しては、電磁波が含水状態に対して鋭敏であって、弾性波や比抵抗等よりも高い分解能を与えるため、電磁波利用のトモグラフィが適当であると考えられる。
【０００４】
図２を参照するに、地盤１の地質分布などを計測する従来のトモグラフィでは、調査すべき断面を選定し、その断面の両端に送信側の計測孔２と受信側の計測孔３とを穿つ。計測孔２内の送信機４からパルス波を送出し、これを計測孔３内の受信機５により受信すると共に、移動手段17によって受信機５を移動し複数の受信位置で受信する。送信機４と各受信位置の受信機５との間に図３に示すようなレイパス18を想定し、各レイパスについて両計測孔２、３間の地盤を透過する際の電磁波の時間遅れ、振幅減衰などを計測する。
【０００５】
更に送信機４の計測孔２内における位置を移動手段17によって順次移動させ、複数の送信位置から電磁波を送信する。各送信位置毎に、送信電磁波を複数の受信位置で受信する。こうして調査すべき断面を覆う図３のレイパス網19を、前記レイパス18の集合として得る。各レイパス18には、例えば両計測孔２、３間における電磁波信号の時間遅れ又は振幅減衰の計測値が付随する。レイパス網19の送・受信位置情報及び前記計測値の情報に対しトモグラフィ処理を施すことによって、図２の速度分布トモグラム22s又は振幅減衰分布トモグラム22aを作成し、地盤の分析を行う。
【０００６】
【発明が解決しようとする課題】
しかし従来の電磁波トモグラフィでは、電磁波としてパルス波を用いていた。このため、送信した電磁波の伝播距離が短く、両計測孔２、３間の距離を大きく取れないという問題点があった。
【０００７】
更にパルス波による信号は、信号の伝播速度や振幅減衰を当該パルス波に含まれる周波数成分（図２のfm〜fn）の全体に対して１つ与えるに過ぎず、いわば全周波数成分に対する平均値として与えるに過ぎない。送信位置から受信位置へ進行する電磁波信号の特性のうち、例えば伝播速度や振幅減衰等の周波数依存性のあるものについて、個別の周波数での計測値ではなく、信号に含まれる全周波数成分に対する平均値として与える場合には、周波数毎の特性が明確に把握できない欠点を生じる。例えば周波数f1の伝播速度が開口亀裂の検出に対して鋭敏であり、周波数f2の伝播速度が亀裂密集帯の検出に対して鋭敏である場合に、周波数f1及びf2に対する伝播速度の平均値のみが与えられ、各周波数f1、f2の個別伝播速度の計測値が与えられない時は、正確な評価が困難になる。
【０００８】
平均値のみが与えられた場合の困難を解決するため、図１１に示すようにスペクトル解析を行い、図１２の振幅スペクトル及び位相スペクトルを得て、個別周波数の振幅及び位相を推定することが理論的には可能である。しかし、実際には計測値における誤差と分析における誤差とが重畳される問題点があり、実用されてはいない。
【０００９】
そこで本発明の目的は、高精度の周波数別トモグラム作像が可能な連続電磁波型トモグラフィ及びその装置を提供することにある。
【００１０】
【課題を解決するための手段】
前記目的を達成するため、本発明者はトモグラフィにおいて、作像すべき断面に透過させる信号を複数周波数の電磁波が段階状に切替わる連続電磁波とし且つ連続電磁波に含まれる周波数毎に送・受信電磁波間の時間遅れなどの特性を計測すれば、周波数別の断面像を正確に作像できることに注目した。
【００１１】
図１及び２を参照するに、本発明の連続電磁波型トモグラフィは、地盤１中の送信位置列内の各送信位置で送信機４から周波数f1〜f4が段階状に切替わる多周波数連続電磁波16（図６参照）を送信し、連続電磁波16が受信可能な所定距離だけ送信位置列から隔てた地盤１中の受信位置列内の受信位置において受信機５により各送信位置から地盤１透過後の連続電磁波16を受信し、送信機４及び受信機５に接続した計測装置９により各受信位置における連続電磁波16の各周波数f1〜f4毎に送・受信電磁波間の時間遅れｔ（図９参照）又は振幅減衰A／A₀（図９参照）を検出し、検出された時間遅れｔ又は振幅減衰A／A₀に基づき送信位置列と受信位置列との間における周波数別の透過型電磁波速度トモグラム22_sf1〜22_sf4又は透過型電磁波振幅減衰トモグラム22_af1〜22_af4を作成し、比較的低い周波数帯の速度トモグラム 22 _sf1 〜 22 _sf2 又は振幅減衰トモグラム 22 _af1 〜 22 _af2 により地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の速度トモグラム 22 _sf3 〜 22 _sf4 又は振幅減衰トモグラム 22 _af3 〜 22 _af4 により地盤中の小さな亀裂の存在可能性を検出してなるものである。
【００１２】
好ましくは、連続電磁波型トモグラフィにおいて透過型電磁波速度トモグラム22_sf1〜22_sf4中の電磁波速度の分布に基づいて、地盤１中の亀裂の大小、粗密、地盤の硬軟その他の性状を検出する。また、前記速度トモグラム22 _sf1 〜 22 _sf4 又は振幅減衰トモグラム22 _af1 〜 22 _af4 に基づき送信位置と受信位置との間の周波数別電気伝導度分布図（図示せず）を作成し、比較的低い周波数帯の電気伝導度分布図により地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の電気伝導度分布図により地盤中の小さな亀裂の存在可能性を検出することもできる。
【００１３】
再び図１及び２を参照するに、本発明の連続電磁波型トモグラム作像装置は、周波数が段階状に切替わる多周波数連続信号16a（図６参照）を発生する信号発生器15；該連続信号16aを連続電磁波16として送信する送信機４；該送信機４を被計測地盤内の所定送信位置列内の各送信位置へ順次移動する移動手段17；送信機４からの連続電磁波16を受信する受信機５；該受信機５を被計測地盤内の所定受信位置列内の各受信位置へ順次移動する移動手段17；送信機４の送信電磁波と受信機５の受信電磁波との間の時間遅れｔ（図９参照）又は振幅減衰A／A₀（図９参照）を連続信号中の周波数別に計測する計測装置９；並びに各送信位置と各受信位置との間に地質分布を想定し、各想定地質の周波数別の電磁波伝播速度又は電磁波減衰率に基づき送・受信電磁波間の時間遅れ又は振幅減衰の周波数別試算値を算出し、周波数別の計測値の該周波数別試算値に対する残差を求め、地質分布の修正の反復により周波数別の残差を所定許容値以下に収束させることにより周波数別の速度トモグラム22 _sf1 〜 22 _sf4又は振幅減衰トモグラム22 _af1 〜 22 _af4を作成し、比較的低い周波数帯の速度トモグラム 22 _sf1 〜 22 _sf2 又は振幅減衰トモグラム 22 _af1 〜 22 _af2 により地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の速度トモグラム 22 _sf3 〜 22 _sf4 又は振幅減衰トモグラム 22 _af3 〜 22 _af4 により地盤中の小さな亀裂の存在可能性を検出する演算手段20を備えてなるものである。
【００１４】
また、本発明の連続電磁波型トモグラム作像プログラムは、周波数が段階状に切替わる多周波数連続電磁波（図６参照）をそれぞれ送信する複数の送信位置の座標と、送信位置からの連続電磁波をそれぞれ受信する複数の受信位置の座標と、送信位置の送信電磁波と受信位置の受信電磁波との間の周波数別の時間遅れｔ（図９参照）又は振幅減衰A／A₀（図９参照）の計測値と、想定される地質の周波数別の電磁波伝播速度及び電磁波減衰率とを入力し；各送信位置と各受信位置との間に地質分布を想定し；想定地質分布と前記電磁波伝播速度及び電磁波減衰率とに基づき前記送・受信電磁波間の周波数別の時間遅れ又は振幅減衰の試算値を算出し；周波数別の前記計測値と該試算値との間の残差を求め；前記地質分布の修正の反復により周波数別の前記残差を所定許容値以下に収束させることにより周波数別の速度トモグラム22 _sf1 〜 22 _sf4又は振幅減衰トモグラム22 _af1 〜 22 _af4を作成し、比較的低い周波数帯の速度トモグラム 22 _sf1 〜 22 _sf2 又は振幅減衰トモグラム 22 _af1 〜 22 _af2 により地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の速度トモグラム 22 _sf3 〜 22 _sf4 又は振幅減衰トモグラム 22 _af3 〜 22 _af4 により地盤中の小さな亀裂の存在可能性を検出してなるものである。
【００１５】
好ましくは、前記周波数別の伝播速度トモグラム22 _sf1 〜 22 _sf4又は振幅減衰トモグラム22 _af1 〜 22 _af4に基づき前記送信位置と受信位置との間の周波数別電気伝導度分布図を作成し、比較的低い周波数帯の電気伝導度分布図により地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の電気伝導度分布図により地盤中の小さな亀裂の存在可能性を検出する（図１のステップ104参照）。
【００１６】
【発明の実施の形態】
図１〜９に示す実施例により本発明の作用を説明する。先ず図１の流れ図のステップ101に示すように、送信機４から送出すべき計測電磁波である多周波数連続信号16aの波形を選択する。この実施例では、信号発生器15によって、図６に例示するような多周波数連続信号16aを発生させる。この信号16aは、異なった多数の周波数部分を含む連続信号であって、各周波数部分においては一定周波数と一定振幅を保つ。連続信号に多数の周波数部分を含めるのは、送受信位置間の電磁波信号の時間遅れ、振幅減衰等のデータを異なる周波数に対して個別に計測して取得するためである。ただし、異なる周波数部分の振幅を等しくする必要はない。
【００１７】
図６（Ａ）の場合には、振幅が常時一定であり、周波数が時間部分T₀〜T₁の部分では一定角周波数ω₀であり、以下時間部分T₁〜T₂、T₂〜T₃、……において、それぞれ角周波数がω₁、ω₂……で一定である。
【００１８】
図６（Ｂ）の例では、時間部分T₀〜T₁の部分で周波数F1が約0.8MHz（≒（1／1.2）×10⁶）であり、次の時間部分T₁〜T₂では周波数が2F1となっており、以下同様に周波数が倍増しているが、振幅は一定である。
【００１９】
図４は、地盤１に深さ30ｍの計測孔２、３を間隔25ｍを隔てて掘削し、計測孔２に送信機４による送信位置を15点想定し、計測孔３に受信機５による受信位置を15点想定した例を示す。
【００２０】
トモグラフィ計測に先立ち本装置の伝播特性を調べるために、図５に示すように、トンネル７の底地盤８に穿った計測孔２内に送信機４を固定し、地表受信機６を底地盤８の表面沿いに移動して計測を行った。計測結果を同図に併記する。従来のパルス波では水平距離20ｍ程度で受信信号がノイズレベル以下になるのに対し、図６の連続電磁波を用いればパルス波による場合の２〜３倍程度の距離まで信号が検出可能であることを確認した。
【００２１】
再び図２を参照するに、計測孔２内に配置した送信機４へ信号発生器15の連続信号16aを送り、連続電磁波16として送信する。計測孔３内に配置した受信機５によって、送信機４からの連続電磁波16を受信し、計測装置９によって送信電磁波と受信電磁波とを比較する。図示例では、この比較のために信号発生器15の連続信号16aを計測装置９経由で送信機４へ送る。図示例では、図１の流れ図のステップ102及び103に示すように送信電磁波と受信電磁波との間の時間遅れｔ（図９参照）及び振幅減衰A／A₀（図９参照）を計測する。計測項目は、時間遅れと振幅減衰に限定されず、他の物理量、例えば電磁波伝播特性などを対象にすることができる。
【００２２】
本発明において重要なことは、送信電磁波と受信電磁波との間の時間遅れ及び振幅減衰の計測を、段階状の連続信号16aに含まれる全ての周波数について個別に行う点にある。即ち、各周波数の時間遅れ及び振幅減衰を、単一のパルス信号による計測値からの計算による誘導値ではなく、直接の計測値として求める。
【００２３】
一つの送信位置からの電磁波を、計測孔３内に想定した受信位置、例えば図４の15位置で計測する。このため、受信機５の位置を移動手段17によって順次調節する。次に、計測孔２内の送信機４を移動手段17により次回の送信位置へ順次移動させ、移動後の送信位置から送出される電磁波を想定した受信位置で受信して前記計測を繰り返す。このような送信機４及び受信機５の移動の反復により、図３に示すようなレイパス網19を作成する。
【００２４】
本発明によればレイパス網19の各レイパス18に対して、複数の周波数における計測値が付随しているので、それらの計測値をデータとして演算手段20に与えれば、周波数別の複数枚の透過型電磁波速度トモグラム及び透過型振幅減衰トモグラムを得ることができる。図６（Ｂ）の例において、時間部分T₀〜T₁、T_１〜T_２、T_２〜T_３、T_３〜T_４の周波数が夫々F1、F2、F3、F4であって、これらの４つの周波数計測が行われたとすると、図２に示すように、周波数別に、４枚の透過型電磁波速度トモグラム22_sf1〜22_sf4及び４枚の透過型電磁波振幅減衰トモグラム22_af1〜22_af4を作成できる。更に、図６（Ａ）のように（ｎ＋１）個の時間部分T₀〜T₁、T_１〜T_２、……、T_ｉ〜T_i _＋１、……、T_n〜T_n _＋１に（ｎ＋１）個の角周波数ω_０、ω₁、……、ω_i、……、ω_nが夫々割当てられた段階状多周波数連続電磁波16を用いる場合は、（ｎ＋１）枚の透過型電磁波速度トモグラム22_s _ω ₀〜22_s _ω _n及び（ｎ＋１）枚の透過型電磁波振幅減衰トモグラム22_a _ω ₀〜22_a _ω _nを作成することができる。レイパス網19及び各レイパス18に対する個別周波数の計測値が与えられた時に、周波数別トモグラム22を作成するには種々の方法を用いることができる。
【００２５】
一例として、物理探査学会の「物理探査ハンドブック、手法編、物理探査学会1998」の778〜779頁に記載の繰返し近似演算による方法を、本発明の理解に必要な限度において図７の流れ図により説明する。ステップ701において、例えば利用できる初期モデルがある場合のデータを演算手段20に入力する。
【００２６】
レイパス網19中の最初に検討するレイパス18のデータをステップ702で読み取る。電磁波の伝播速度Vと電気伝導度Ecとの間には次の式（1）の関係があるので、地質分布を想定すれば、レイパス18に沿う各周波数成分毎の電磁波伝播の遅れを算出することができる（物理探査学会、第６回物理探査基礎講座、講習会テキスト、1990年6月、第254頁参照）。ステップ703で最初のレイパス18に沿って想定した特定周波数の速度分布を、前記地質分布から割出す。割出した値を、当該レイパス18の速度分布仮定値として演算手段20上に書込む。ステップ703において、そのレイパス18上の速度分布を演算手段20によってトレイスしながら送信信号と受信信号との間の時間遅れを試算する。ステップ704において、その試算値と、ステップ702で読取った伝播時間、即ち時間遅れｔの実測値との間の差を、残差として計算する。残差のこの計算は、前記繰返し演算の一環として行う。次いでステップ705において、他のレイパス18について前記残差を計算して中間モデルを作成する。
【００２７】
【数１】
Ec＝（C₀／60π）・（A／A₀）・（1／V） ………………（1）
ここに、Ec：電気伝導度
A／A₀：振幅減衰
V：電磁波の伝播速度
C₀：空気中の電磁波の伝播速度
【００２８】
前記中間モデルの速度分布仮定値が収束しているか否かを、ステップ706において例えば今回の中間モデルの前記残差と前回の残差との違いの大きさで判定する。否と判定された場合は、当該レイパス18に沿う速度分布の想定値を更新し、前記の収束判定を繰り返す。収束と判定された時には、ステップ707において最終速度モデル、即ち速度トモグラム22sが得られる。他の周波数に対する速度トモグラムも同様な手法により得られる。更に、振幅減衰トモグラム22aも、伝播時間の替わりに振幅減衰を用い、同様な手法により作像することができる。
【００２９】
こうして連続電磁波16に含まれる全ての周波数別のトモグラム22を作像することができる。従って、本発明の目的である高精度の周波数別トモグラム作像が可能な連続電磁波型トモグラフィ及びその装置の提供が達成される。
【００３０】
図１のステップ104は、例えば、前記式（1）の関係を用い、速度トモグラム22s又は振幅減衰トモグラム22aから電気伝導度分布図を作る過程を示す。式（1）が利用可能である場合には、この過程を比例計算によって行うことができる。
【００３１】
【実施例】
図８は、三つの周波数帯、即ち2.554〜3.781MHz、3.781〜5.008MHz及び6.234〜7.4MHzの各々に対する周波数帯別速度トモグラムを、図５のトンネル７の底地盤８に対して実測により作成した結果の例示である。比較的低い周波数帯を用いることにより、地盤中の亀裂密集域の明確な検出が、図８の2.554〜3.781MHz帯速度トモグラムにおける電磁波伝播速度が特に低い領域10として可能になる。亀裂密集域は、電磁波伝播速度の低いことが知られているからである。図１の流れ図におけるステップ105の地質構造の把握は、このような地質構造中の特定性質部分の検出把握を示す。図示例の亀裂密集域の領域10は、本発明のトモグラフィによリ作成したトモグラムにおいて、電磁波速度が、予め選定された最低電磁波速度40m/nsec以下である部分の中に他の関連要因をも考慮して推定されたものである。
【００３２】
図１０の「トンネル壁面の亀裂」の欄は、トンネル底面に亀裂が機械的な手法により観察されることを示し、図８における亀裂密集帯10、小さな亀裂11、及び開口亀裂12の位置に関する判定結果の妥当性を裏付けている。
【００３３】
図８の3.781〜5.008MHz帯域の速度トモグラムは開口亀裂12を、電磁波伝播速度が特に低い領域として示している。更に6.234〜7.4MHz帯域の速度トモグラムは小さな亀裂11を、電磁波伝播速度がやや低い領域として示している。図１０の「トンネル壁面の亀裂」の欄は、トンネルの底地盤８及び計測孔２、３の内面におけるこれらの開口亀裂12及び小さな亀裂11の端面が、テレビカメラによる撮像と画像観察などの手法により検出できることを示す。即ち、本発明のトモグラフィ計測による開口亀裂12及び小さな亀裂11の判定結果の妥当性を図１０の観察結果が裏付けている。図示例の開口亀裂12の存在は、前記最低電磁波速度40m/nsec以下である部分の中に関連要因をも勘案して推定されたものである。また、小さな亀裂11の存在は、予め選定された限界電磁波速度60m/nsec以下であるが前記最低電磁波速度40m/nsec以上である部分の中に他の関連要因をも考慮して推定されたものである。
【００３４】
以上の説明において、地盤構造の断面像への本発明のトモグラフィの適用に言及したが、本発明の利用は地盤構造のみに限定されるものではなく、地下埋蔵物の非接触検出などにも応用できる。また、トモグラムは二次元に限定されるものではなく、電磁波透過面を複数にし各透過面のトモグラムを立体的に重畳して三次元トモグラムを作像することも可能である。
【００３５】
【発明の効果】
以上説明したように本発明の連続電磁波型トモグラフィは、複数周波数が含まれる連続電磁波を使用して周波数別のトモグラムを作成するので、次の顕著な効果を奏する。
【００３６】
（イ）周波別の出力を発生するので、広い周波数帯域が含まれるが持続時間の短いパルスによる従来のトモグラフィでは、出力が全周波数帯域に平均的となるために検出できなかった現象を、高精度で確実に評価できるようになる。
（ロ）持続時間の短いパルスではエネルギー不足で評価できなかった事柄を、多周波数の連続電磁波の大きなエネルギーによって確実に評価できる。
（ハ）亀裂の大小や亀裂密度の高低等の地盤性状を、異なる周波数のトモグラムによって識別することが期待できる。
（ニ）地盤構造の三次元計測が理論的に可能になる。
（ホ）岩盤の詳細な含水状態の検出が可能になる。
（ヘ）送信機から送出する電磁波の到達距離を長くすることにより、トモグラフィによって評価できる面積を広げることができる。
【図面の簡単な説明】
【図１】は、本発明によるトモグラフィの流れ図である。
【図２】は、本発明のトモグラフィによる計測及び演算過程の図式的説明図である。
【図３】は、電磁波信号が伝播するレイパス網を示す説明図である。
【図４】は、送信機及び受信機が配置される計測孔の一例の説明図である。
【図５】は、送受信機の配置の一例の説明図である。
【図６】は、本発明で使う段階状多周波数連続信号の波形の図式的説明図である。
【図７】は、トモグラム作成演算の一例の流れ図である。
【図８】は、トモグラムの一例の説明図である。
【図９】は、送信電磁波及び受信電磁波の波形を示す図である。
【図１０】は、本発明のトモグラフィによる計測結果の一例を示す図である。
【図１１】は、従来のスペクトル分析の一例の流れ図である。
【図１２】は、従来のスペクトル分析の結果の一例を示すグラフである。

Claims

地盤中の送信位置列内の各送信位置で送信機から周波数が段階状に切替わる多周波数連続電磁波を送信し、前記連続電磁波が受信可能な所定距離だけ送信位置列から隔てた地盤中の受信位置列内の受信位置において受信機により前記各送信位置からの地盤透過後の連続電磁波を受信し、前記送信機及び受信機に接続した計測装置により各受信位置における前記連続電磁波の各周波数毎の送・受信電磁波間の時間遅れ又は振幅減衰を検出し、検出した時間遅れ又は振幅減衰に基づいて前記送信位置列と受信位置列との間における周波数別の透過型電磁波速度トモグラム又は透過型電磁波振幅減衰トモグラムを作成し、比較的低い周波数帯の前記トモグラムにより地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の前記トモグラムにより地盤中の小さな亀裂の存在可能性を検出してなる連続電磁波型トモグラフィ。
請求項１のトモグラフィにおいて、前記周波数別の速度トモグラム又は振幅減衰トモグラムに基づき前記送信位置と受信位置との間の周波数別電気伝導度分布図を作成し、比較的低い周波数帯の電気伝導度分布図により地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の電気伝導度分布図により地盤中の小さな亀裂の存在可能性を検出してなる連続電磁波型トモグラフィ。
請求項１又は２のトモグラフィにおいて、前記連続電磁波の周波数を最低周波数から最高周波数まで所定継続時間毎に段階的に切替え、前記所定継続時間中は一定周波数としてなる連続電磁波型トモグラフィ。
請求項１から３の何れかのトモグラフィにおいて、前記透過型電磁波速度トモグラム中の所定最低電磁波速度以下の領域により大きな亀裂又は亀裂密集域の存在可能性を検出し、前記所定最低電磁波速度より高く所定の限界電磁波速度以下の領域により小さな亀裂の存在可能性を検出してなる連続電磁波型トモグラフィ。
周波数が段階状に切替わる多周波数連続信号を発生する信号発生器；該連続信号を連続電磁波として送信する送信機；該送信機を被計測地盤内の所定送信位置列内の各送信位置へ順次移動する移動手段；前記送信機からの連続電磁波を受信する受信機；該受信機を被計測地盤内の所定受信位置列内の各受信位置へ順次移動する移動手段；前記送信機の送信電磁波と受信機の受信電磁波との間の時間遅れ又は振幅減衰を前記連続信号中の周波数別に計測する計測装置；及び前記各送信位置と各受信位置との間に地質分布を想定し、各想定地質の周波数別の送・受信電磁波間の時間遅れ又は振幅減衰の周波数別試算値を算出し、周波数別の前記計測値に対する該周波数別試算値の残差を求め、前記地質分布の修正の反復により周波数別の前記残差を所定許容値以下に収束させることにより周波数別の速度トモグラム又は振幅減衰トモグラムを作成し、比較的低い周波数帯の前記トモグラムにより地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の前記トモグラムにより地盤中の小さな亀裂の存在可能性を検出する演算手段を備えてなる連続電磁波型トモグラム作像装置。
請求項５の作像装置において、前記演算手段により、前記周波数別の速度トモグラム又は振幅減衰トモグラムに基づき前記送信位置と受信位置との間の周波数別電気伝導度分布図を作成し、比較的低い周波数帯の電気伝導度分布図により地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の電気伝導度分布図により地盤中の小さな亀裂の存在可能性を検出してなる連続電磁波型トモグラム作像装置。
周波数が段階状に切替わる多周波数連続電磁波をそれぞれ送信する複数の送信位置の座標と、前記送信位置からの連続電磁波をそれぞれ受信する複数の受信位置の座標と、前記送信位置の送信電磁波と受信位置の受信電磁波との間の周波数別の時間遅れ又は振幅減衰の計測値と、想定される地質の周波数別の電気伝導度及び電磁波減衰率とを入力し；前記各送信位置と各受信位置との間に地質分布を想定し；想定地質分布と前記電気伝導度及び電磁波減衰率とに基づき前記送・受信電磁波間の周波数別の時間遅れ又は振幅減衰の試算値を算出し；周波数別の前記計測値と該試算値との間の残差を求め；前記地質分布の修正の反復により周波数別の前記残差を所定許容値以下に収束させることにより周波数別の速度トモグラム又は振幅減衰トモグラムを作成し、比較的低い周波数帯の前記トモグラムにより地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の前記トモグラムにより地盤中の小さな亀裂の存在可能性を検出してなる連続電磁波型トモグラム作像プログラム。
請求項７のプログラムにおいて、前記周波数別の速度トモグラム又は振幅減衰トモグラムに基づき前記送信位置と受信位置との間の周波数別電気伝導度分布図を作成し、比較的低い周波数帯の電気伝導度分布図により地盤中の大きな亀裂又は亀裂密集域の存在可能性を検出し、比較的高い周波数帯の電気伝導度分布図により地盤中の小さな亀裂の存在可能性を検出してなる連続電磁波型トモグラム作像プログラム。