JP4229371B2 - 地中空洞探査方法 - Google Patents
地中空洞探査方法 Download PDFInfo
- Publication number
- JP4229371B2 JP4229371B2 JP2003147120A JP2003147120A JP4229371B2 JP 4229371 B2 JP4229371 B2 JP 4229371B2 JP 2003147120 A JP2003147120 A JP 2003147120A JP 2003147120 A JP2003147120 A JP 2003147120A JP 4229371 B2 JP4229371 B2 JP 4229371B2
- Authority
- JP
- Japan
- Prior art keywords
- cavity
- exploration
- underground
- reaction
- frequency
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 96
- 238000005259 measurement Methods 0.000 claims description 93
- 238000006243 chemical reaction Methods 0.000 claims description 86
- 230000000694 effects Effects 0.000 claims description 54
- 102100021809 Chorionic somatomammotropin hormone 1 Human genes 0.000 claims description 19
- 101000895818 Homo sapiens Chorionic somatomammotropin hormone 1 Proteins 0.000 claims description 19
- 230000006698 induction Effects 0.000 claims description 10
- 230000005856 abnormality Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 55
- 230000002093 peripheral effect Effects 0.000 description 34
- 238000004364 calculation method Methods 0.000 description 17
- 230000005540 biological transmission Effects 0.000 description 15
- 230000005669 field effect Effects 0.000 description 14
- 230000006870 function Effects 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 9
- 238000004141 dimensional analysis Methods 0.000 description 8
- 238000003199 nucleic acid amplification method Methods 0.000 description 8
- 230000003321 amplification Effects 0.000 description 7
- 239000002689 soil Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000005684 electric field Effects 0.000 description 5
- 238000011835 investigation Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000010587 phase diagram Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000005672 electromagnetic field Effects 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 239000003990 capacitor Substances 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 239000003673 groundwater Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000009933 burial Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 229910000679 solder Inorganic materials 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、高周波ＣＳＭＴ法により地中空洞の位置を高精度に推定する地中空洞探査方法に関する。
【０００２】
【従来の技術】
従来、地表から２ｍ以上の浅部の探査においては地中レーダーが有効とされている。しかし、この探査方法は信号の減衰が大きいので、深度１０ｍを超えると事実上探査ができなくなるものであった。地下の空洞を探査するときにも２ｍ〜１０ｍが限界であった。
【０００３】
電気探査方法にはこのような制限がなく、これを使用して探査することも考えられる。しかし、電気探査方法は深度が深くなると水平分解能が低下する傾向がある。例えば、地下の状態を推定する有力な手法として、電気比抵抗を測定してそのコントラストを測定する比抵抗調査法がある。中でも、直流比抵抗法を使用した高密度電気探査法（比抵抗影像法）がもっとも一般的で普及している。この直流比抵抗法は、探査深度を深くするとき電極間隔を広げる必要があり、これによって水平分解能が低下してしまう。しかも、測定時の作業負担が大きく、探査費用も嵩むものであった。
【０００４】
また、このような比抵抗調査法の一種である電磁探査法にＭＴ法（地磁気地電流探査）と可聴領域（Ａ）でコントロール（Ｃ）された電磁波を発信源（Ｓ；ソース）から発信し測定するＣＳＡＭＴ法があるが、これらは周波数を変化させて深部を探査するため電極間隔を狭く保て、測定点が信号源から十分離れている場合に電磁場を平面波とみなして測定するため、条件を満たせば比較的高い精度が期待できる。電磁波を利用するため作業負担、費用に関しても優れている。
【０００５】
このうち、ＭＴ法は地球磁気圏や雷放電により発生する自然電磁場の１０−３Ｈｚ〜２５０Ｈｚの周波数を測定するものであり、自然電磁場を利用して地下の比抵抗構造を分析するため、データが得られなかったり、ノイズに弱く、このような意味で精度が高いとはいえなかった。
【０００６】
これに対して、人工信号源からの電磁波を測定するＣＳＡＭＴ法はノイズに強く、測定費用も安価であり、周波数を変化させて探査するので電気探査、直流比抵抗法と比較して深部における水平分解能が高いという利点を有している。しかし、低周波数域ではニアフィールド効果を受け易いという欠点も有している。このニアフィールド効果は、このような低周波域においては平面波の仮定が崩れて見掛比抵抗が増加し、直流電界が強く、測定点直下の正確な比抵抗や位相を得るのが難しくなるものである。
【０００７】
このニアフィールド効果を減らす１つの手段として、従来、水平磁気双極子アンテナが採用されてきた。しかし、これによってＣＳＡＭＴ法に伴うニアフィールド効果が根本的に解消されるようなものではない。また、深度１０ｍ〜１００ｍの浅部に存在する高比抵抗の空洞を探査するのは、高周波情報を利用していない従来のＣＳＡＭＴ法では事実上困難と考えられてきた。
【０００８】
そこでニアフィールドの問題に対しては、ＣＳＭＴ法の特徴を活かしてニアフィールド効果を除く方法が提案された（特許文献１参照）。これは、ＣＳＭＴ法による探査とＭＴ法による探査を重複して行い、ＣＳＭＴ法によって得たデータのうち、ニアフィールド効果が発生する周波数領域を、ＭＴ法によるデータによって置き換えて地下構造を解析するものである。しかしながら、この方法ではＭＴ法のノイズに弱いという欠点はそのまま残るものであり、また、数百〜数千ｍの大深度の調査に適応されるものであるため、深度１０ｍ〜１００ｍの浅部の調査にＣＳＡＭＴ法は利用されていなかった。
【０００９】
そこで本発明者らは、ＣＳＭＴ法の水平分解能が良いという特徴を活かして深度１００ｍまでの浅部を探査するため、高周波ＣＳＭＴ法による探査装置を提案した（非特許文献１参照）。１ｋＨｚ〜１００ｋＨｚ（１６測点）で測定を行うものである。アンテナは１〜１０ｋＨｚの低周波と、１０ｋＨｚ〜１００ｋＨｚの高周波の２帯域で異なったループアンテナを使っている。高周波を用いることにより、従来のＣＳＭＴ法と比較して、（１）装置、とくに送信部が小型、軽量化され、（２）探査時間が１周波数について１０秒、１測点について１６０秒と、大幅な短縮が可能であり、（３）水平方向の分解能が高く、操作性もよく、地形の制約を受けることが少なく、探査費用が低減できるものであった。そして、非特許文献１で提案した探査装置で分析した比抵抗構造は、高密度電気探査法とほぼ同じになるという結果を得ている。
【００１０】
【特許文献１】
特開平５−３２３０３８号公報
【非特許文献１】
半田，外３名，「高周波ＣＳＭＴ装置の開発」，物理探査，物理探査学会，５５巻，３号，ｐ．１９９−２０６
【００１１】
【発明が解決しようとする課題】
以上説明したように、従来のＣＳＭＴ法は低周波数域ではニアフィールド効果を受け易いという問題を有していた。ニアフィールド領域では、直流電界が強くなり、その結果、見掛比抵抗（以下、比抵抗）が増加し、解析による測定点直下の「真の比抵抗」の算出ができなくなる。このため、従来は如何にニアフィールド効果を低減させるかが、探査性能の向上に重要と考えられていた。
【００１２】
図１２は従来から、ニアフィールド効果が少ないので多く使用されてきた水平磁気双極子アンテナを使用したＣＳＭＴ法による比抵抗と位相の測定図である。水平磁気双極子アンテナを使用したＣＳＭＴ法によって空洞探査した場合の事例を示す。この事例では、同じ発信地点で水平磁気双極子アンテナを使用していることから、ニアフィールド効果は垂直磁場双極子アンテナに比べ少ないと考えてよい。比抵抗に関しては図１２に示すように、周波数１．４ｋＨｚ〜９７ｋＨｚで１０２Ω・ｍ以下で変化に乏しく、１３ｋＨｚで僅かに空洞部で高比抵抗となっているに過ぎない。同様に、位相に関しても図１２に示すように、周波数１．４ｋＨｚ〜４５ｋＨｚでは変化の乏しい位相６０°〜８０°を示すだけであり、周波数４５ｋＨｚ〜９７ｋＨｚの範囲でも空洞を明確に示唆するような情報を含まないものである。このように、比抵抗および位相の地盤応答だけでは、高周波ＣＳＭＴ探査機を用いても空洞を明確に識別できない。
【００１３】
また、非特許文献１で説明した本発明者らの高周波ＣＳＭＴ法による探査装置は、深さ１００ｍ以下の浅部探査であるためにループアンテナを採用し、ＰＬＬ回路を使用してフーリエ変換するが、参照信号と送信信号との差で出力が変動してしまい、高周波分解能が悪く、高Ｓ／Ｎ比を実現するのが難しいものであった。そして高密度電気探査法と同程度の探査が可能になるという比較結果は得られたが、上述の出力変動でＳ／Ｎ比が上がらず、分析した比抵抗構造は粗くて、空洞位置を正確に予測することなど期待できないものであった。
【００１４】
空気が充満する空洞においては地盤より際立って比抵抗が大きくなるが、非特許文献１は、こうした空洞の場合について未研究で、空洞の比抵抗が示す挙動、位相が示す挙動のいずれもがまったく未知であった。そして、こうした空洞様反応が得られたとしても、空洞位置をどのように分析し、判定するのか、その手法はまったく未解決であった。
【００１６】
そこで本発明は、浅部の空洞の位置を精度よく安価に探査することができる地中空洞探査方法を提供することを目的とする。
【００１７】
【課題を解決するための手段】
以上説明した課題を解決するために本発明の地中空洞探査方法は、地中探査のため高周波信号を電磁波として送信するとともに、一対の地電流電極で測線上の測点の電位差を検出するとともに該測線と直交する方向に配設されたインダクションコイルにより該測点の磁場を測定し、該測点における比抵抗と位相を演算する高周波ＣＳＭＴ法による地中空洞探査方法であって、空洞の所在を示す空洞様反応をニアフィールド効果によって増幅するための垂直磁気双極子アンテナを空洞から所定距離おいた位置に配置して高周波信号を送信し、探査領域における周波数と測点に関する等比抵抗データを取得して、該等比抵抗データの中からニアフィールド効果により増幅された空洞様反応領域を抽出し、該空洞様反応領域の測線方向幅の中央を空洞中心として推定することを特徴とする。
【００１８】
本発明は、この構成により、ニアフィールド効果を避けるのではなく、これを利用し、ニアフィールド効果による空洞様反応を増幅させて可視化することによって、浅部の空洞の位置を精度よく安価に且つ自動的に探査することができる。
【００１９】
【発明の実施の形態】
本発明の第１の発明は、地中探査のため高周波信号を電磁波として送信するとともに、一対の地電流電極で測線上の測点の電位差を検出するとともに該測線と直交する方向に配設されたインダクションコイルにより該測点の磁場を測定し、該測点における比抵抗と位相を演算する高周波ＣＳＭＴ法による地中空洞探査方法であって、空洞の所在を示す空洞様反応をニアフィールド効果によって増幅するための垂直磁気双極子アンテナを空洞から所定距離おいた位置に配置して高周波信号を送信し、探査領域における周波数と測点に関する等比抵抗データを取得して、該等比抵抗データの中からニアフィールド効果により増幅された空洞様反応領域を抽出し、該空洞様反応領域の測線方向幅の中央を空洞中心として推定することを特徴とする地中空洞探査方法であり、高周波の電磁波を送信する高周波ＣＳＭＴ法による探査を基本とし、これまでの探査手法では測定精度が悪かった地下１０ｍ以深（１０ｍ〜１００ｍ）を含む地下１ｍ〜１００ｍ程度の浅部から「大深度地下の公共的使用に関する特別措置法」の定める大深度地下の空洞位置を検出でき、垂直磁気双極子アンテナを使用してニアフィールド効果により空洞の所在を示す空洞様反応を増幅することによって比抵抗における空洞様反応を強調でき、高Ｓ／Ｎ比で、精度よく安価に探査することができる。空洞様反応領域の中から、ニアフィールド効果でコントラストが強調された等比抵抗データの空洞様反応を抽出するから、空洞中心の判別が正確になる。
【００２８】
本発明の第２の発明は、第１の発明に従属する発明であって、探査領域における周波数と測点に関する等比抵抗データを取得して、該空洞様反応領域の測線方向幅の中央を空洞中心として推定するのに代えて、探査領域における周波数と測点に関する等位相データの中から空洞様反応のピークを含む低位相の空洞様反応領域を抽出し、この空洞様反応領域に対応する等比抵抗データの中から比抵抗のピークを抽出して空洞中心を推定することを特徴とする地中空洞探査方法であり、等位相データからピークを抽出して空洞様反応領域の存在を確認し、さらにニアフィールド効果でコントラストが強調された等比抵抗データの周縁効果と中心反応示す空洞様反応のピークを抽出することにより、空洞中心の判別がさらに正確になる。
【００３２】
本発明の第３の発明は、第１又は第２の発明に従属する発明であって、測線の両側から電磁波を送信し、該測線の直下にある空洞と測線の直下に無い空洞とを、シャドウ効果による空洞様反応に近似した比抵抗または位相異常の有無から判別することを特徴とする地中空洞探査方法であり、側線の直下の空洞と、直下にない空洞とを判別できる。
【００３７】
（実施の形態１）
以下、本発明の実施の形態１における地中空洞探査方法について説明する。実施の形態１は、空洞上に測線が空洞に対し、ほぼ直交方向に設置されている場合に空洞位置を把握するものである。なお、実施の形態１の地中空洞探査装置は、測点間距離等の測定条件を変更し、また送受信方向を選択して、繰り返し測定を行うことにより、この測定条件を満たす状態を実現し、実際の測定を行うものである。図１は本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置の構成図、図２（ａ）は本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置の測線と送信部のＸ−Ｙ平面配置図、図２（ｂ）は（ａ）の地中空洞探査方法を行うための地中空洞探査装置の測線とセンサと地下空洞のＸ−Ｙ−Ｚ配置図、図２（ｃ）は（ａ）の地中空洞探査装置における空洞様反応の出現位置の説明図、図３（ａ）は比抵抗モデルの説明図、図３（ｂ）は（ａ）の比抵抗の説明図、図３（ｃ）は（ａ）の位相の説明図、図４は本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置が空洞位置に対して測定した見掛比抵抗と位相の測定図である。
【００３８】
図１において、１は地中空洞探査装置から信号波を送信するために発振を行う発振部、２は発振部１で発振した１００Ｈｚの基本周波数が入力されると１ｋＨｚ〜１００ｋＨｚの周波数帯を対数的に１６分割した周波数を得るため基本周波数の逓倍を行う周波数逓倍部、２ａは周波数逓倍部２からの信号で±１３Ｖの電圧で矩形波信号を形成するスイッチング部、３は矩形波信号を増幅するドライバ部、４ａはループコイルとコンデンサで共振回路を構成した後述するニアフィールド効果を増幅するための高周波アンテナ、４ｂは同じく、ループコイルとコンデンサで共振回路を構成したニアフィールド効果を増幅するための低周波アンテナである。高周波アンテナ４ａは９７ｋＨｚ，７３ｋＨｚ，５５ｋＨｚ，４２ｋＨｚ，３１ｋＨｚ，２４ｋＨｚ，１８ｋＨｚ，１３ｋＨｚの８つの周波数を送信し、低周波アンテナ４ｂは１０ｋＨｚ，７．６ｋＨｚ，５．８ｋＨｚ，４．３ｋＨｚ，３．３ｋＨｚ，２．５ｋＨｚ，１．８ｋＨｚ，１．４ｋＨｚの８つの周波数を送信する。そして、この高周波アンテナ４ａと低周波アンテナ４ｂはニアフィールド効果が強く現われる特性をもつ垂直磁気双極子アンテナを構成するものである。なお、周波数がこれらの周波数に限られないのは当然である。
【００３９】
５は中央処理装置（以下、ＣＰＵ）等から構成され、制御プログラムを読み込んで機能実現手段として機能を実現する地中空洞探査装置の送信側制御部である。６は１６の周波数の信号波を１０秒ごとに切り替えるための計時手段、７は制御プログラム等を記憶した記憶部、８は発振部１が発振した１００Ｈｚの基本周波数の信号を地中空洞探査装置の受信側に送って同期をとるための光ファイバ、８ａは発振部１の電気信号を光信号に変換する信号変換部、８ｂは光信号を電気信号に変化する信号変換部、９は発振部１と信号変換部８ａから構成される地中空洞探査装置の同期部である。
【００４０】
１１ａは共振周波数９７ｋＨｚのインダクションコイルＨから構成される高周波用磁気センサ、１１ｂは共振周波数９．７ｋＨｚのインダクションコイルＬから構成される低周波用磁気センサ、１２ａ，１２ｂはステンレス製の地電流電極、１３ａ，１３ｂはプリアンプ、１４ａ，１４ｂは入力インピーダンスを調整するバッファアンプである。１５，１６はバンドパスフィルタ、１７，１８は増幅部、１９はフーリエ変換部、１９ａはＰＬＬ回路を使って余弦波成分と正弦波成分でフーリエ変換する高周波変換部、１９ａはＰＬＬ回路を使って余弦波成分と正弦波成分でフーリエ変換する低周波変換部、２０はＡ／Ｄ変換部である。
【００４１】
２１はＣＰＵ等から構成され、制御プログラムを読み込んで機能実現手段として演算機能，システム制御機能等を実現する地中空洞探査装置の制御部、２２はデジタル化された測定値から各演算を行う機能実現手段としての演算手段、２３は地中空洞探査装置のシステム制御を行う機能実現手段としてのシステム制御手段である。システム制御手段２３は高周波用磁気センサ１１ａ，低周波用磁気センサ１１ｂの切り替え、増幅部１７，１８のゲインやバンドパスフィルタ１５，１６の中心周波数の切り替え等を行う。
【００４２】
２４はフーリエ変換部１９からの出力を使って比抵抗を演算する比抵抗算出手段、２５はフーリエ変換部１９の出力から位相を演算する位相算出手段、２６は比抵抗構造を解析するためボスティック法等による１次元解析断面の演算を行う解析手段である。後述するニアフィールド効果を避けて解析するため、高周波側のデータだけで解析するのも好適である。1次元インバージョン解析では深度誤差が大きい場合、空洞上部の地盤の緩みや崩壊、あるいは空洞とは無関係の地質要因による異常値を示す可能性があるが、解析手段２６でさらに２次元インバージョンを行うことにより、空洞位置の解析精度を上げることができる。
【００４３】
２７は位相算出手段２５が算出した周波数と測点の等位相図（同一位相の周波数と測点を示すもので、本発明の等位相データ）に基づいて、低位相の空洞様反応領域から中心反応と周縁効果を示す一対の空洞様反応のピークを抽出し、また、比抵抗算出手段２４によって作成された比抵抗の等比抵抗図（同一比抵抗の周波数と測点を示すもので、本発明の等比抵抗データ）に基づいて、高比抵抗の空洞様反応領域から中心反応と周縁効果を示す空洞様反応のピークを抽出し、このピークから地下空洞の位置を算出する空洞位置算出手段である。この空洞様反応領域と、それに含まれる中心反応と周縁効果の詳細については後述する。
【００４４】
実施の形態１の空洞位置算出手段２７は、位相に関する周縁効果のピークをチェックし、さらに高比抵抗の空洞様反応領域から中心反応と周縁効果を示す空洞様反応のピークを抽出して空洞中心を推定する。このとき、等位相図で周縁効果を示す一対のピークからその中央位置を計算し、この中央位置を簡易に空洞の中心位置として推定することも可能である。なお、中心反応と周縁効果を示す空洞様反応のピークは最適な測定条件下において整った形態で出現するが（後述の図２（ｃ）参照）、状況によっては顕著には出現しない場合もある。このため実施の形態１の地中空洞探査装置においては、地電流電極間の距離を大きい間隔から徐々に小さくして測定し、また受信方向、その距離を変更し、あるいはその他の測定条件を変更して、繰り返し測定を行うことにより最適な測定条件を実現している。最適条件を満たしているか否かは、中心反応と周縁効果のピークが整ったパターンで得られたことで判断する。
【００４５】
以上説明した実施の形態１の地中空洞探査装置は、基本的に、高周波アンテナ４ａ、低周波アンテナ４ｂからそれぞれ送信信号を送信する送信機部と、高周波用磁気センサ１１ａと低周波用磁気センサ１１ｂで受信する受信機部とから構成される。なお、実施の形態２で説明するように送信機部は受信機部と同期をとって２台、あるいは複数設けることができる。
【００４６】
次に、実施の形態１の地中空洞探査装置と地中空洞探査方法で空洞位置を決定できる原理を説明する。まず、実施の形態１の地中空洞探査装置で探査を行うときの測線の設け方、測定の方法について説明する。本発明が探査する典型的な場合は、地下に空洞が存在することが予測される調査地で、どこに空洞があるかを探査する場合等であり、測定条件を変更しての測定はそれほど面倒になることはない。なお、空気が充満した空洞を探査する事例で説明するが、地下水が溜まった空洞や土の陥没で埋まった空洞は、空洞に溜まった地下水や充填土と周辺地質の比抵抗の違いによって探査するため、空洞様反応の感度が低くなる。
【００４７】
実施の形態１では、地下空洞に対して図２（ａ）（ｂ）に示すように測線と送信機Ａ、地電流電極１２ａ，１２ｂ、高周波用磁気センサ１１ａ，低周波用磁気センサ１１ｂのインダクションコイルＨ，Ｌを配置して測定を行う。想定される地下の空洞に直交すると思われる方向（Ｙ方向）に測線を設置し、測線の側方（Ｘ方向）に５０ｍ〜１００ｍ測線から離した地点に送信機Ａを置く。また、受信機の地電流電極１２ａ，１２ｂを測線方向に，インダクションコイルＨ，Ｌを測線の直交方向に設置する。地電流電極１２ａ，１２ｂの間隔は、想定される空洞幅の１／２〜３／４程度とする。測定間隔は、地電流電極１２ａ，１２ｂの間隔と同一か、短くする。なお、上述したように、中心反応と周縁効果のピークが整った形で得られたときが最も適切な測定条件のときであるが、この状態は地電流電極１２ａ，１２ｂと空洞間隔を上記の通り調整することで実現できる。1測点に付き送信機Ａからの信号を測定する。
【００４８】
上述のように測線とインダクションコイルを配置したときは、図２（ａ）（ｂ）のように空洞は電場方向に対して直交する方向（磁場方向）となる。このような電場方向に対して直交する方向または高角度で斜交する状態（以下、ＴＭモードという）においては、空洞と地盤との境界では比抵抗が不連続となり、この不連続点で比抵抗が極端な変化を示す。なお、これが逆のものがＴＥモードである。図３（ａ）に示した比抵抗モデルは、比抵抗ρ１＝１０Ωｍの媒体と比抵抗ρ２＝１００Ωｍの媒体が接しているものである。図３（ｂ）（ｃ）の×点は、比抵抗ρ１＝１０Ωｍの媒体と比抵抗ρ２＝１００Ωｍの媒体が接しているときの、ＴＭモードにおける境界での比抵抗、位相の変化を示している。○点は、空洞を電場方向と平行に向けた場合の比抵抗、位相の変化である。ＴＭモードにおいては２つの媒体の境界面で比抵抗が鋭敏な変化をし、位相は高比抵抗領域に入ると低下している。
【００４９】
これを地中の空洞の場合で説明すると、図２（ａ）（ｃ）に示すように、ＴＭモードの空洞と地盤の境界において比抵抗は地盤側では非常に小さく空洞側で非常に大きくなり、コントラストが際立つものとなる。同様に位相においても、地盤側境界で大きなピーク値をとって空洞では低位相を示すことを示す。
【００５０】
しかし、実際の空洞においては少し複雑になる。すなわち、図３の単純な２次元構造の場合と異なり、図２（ｃ）に示すように空洞の直上の測点で比抵抗は高く、位相は低くなるが、空洞側壁を跨がない最も空洞に近い測点でも同様に、比抵抗は高く、位相は低くなる。そして、この高比抵抗および低位相の異常値は空洞から遠ざかるほど小さくなる。この現象は、現場試験と数値解析によるモデル実験で確認することができ、ＴＭモードにおける空洞に対するＭＴ応答の特徴として定義できるものである。そして図２（ｃ）の状態は測線が空洞上を直交して交差した場合の測点、地電流電極間隔および空洞位置に関係する空洞様反応のピークパターンで、これらの位置関係が悪いとピークが隠れてしまうことがある。
【００５１】
このＭＴ応答の特徴を挙げると、（１）空洞直上は高比抵抗と低位相のピークとなる。本明細書においてはこれを「中心反応」という。また、（２）空洞側壁を跨がない最も空洞に近い測点は高比抵抗と低位相のピークとなる。本明細書においてはこれを「周縁効果」という。（３）周縁効果は空洞から離れるほど小さくなる。本明細書においては「中心反応」「周縁効果」を一括して「空洞様反応」という。とくに、中心反応と空洞両翼に生じる一対の周縁効果が生じている範囲を本明細書においては「空洞様反応域」と定義する。ここで、「様」を付加したのは実際の空洞調査において，この特徴を持ったＭＴ応答が現れたとしても、直ちに空洞と断定できないので、「空洞あるいは空洞に似た反応」のニュアンスを残すためである。
【００５２】
そこで、本発明の地中空洞探査装置と地中空洞探査方法の基本となるこの空洞様反応の出現の仕方を、実測した等位相図と等比抵抗図で説明する。図４は、地中空洞Ａが測点５、地中空洞Ｂが測点１１に存在するもので、図１２と同じ地中空洞Ａ，地中空洞Ｂを図２（ａ）（ｂ）（ｃ）に示すような測点で同一の条件で測定したものである。後述の実施例１，２と同一の測定対象で測定したものである。空洞の天井までの深度は１３．７ｍ、空洞の高さは２．２ｍで、この比を意味する土被比は約６．３である。
【００５３】
図４において、測点Ｎｏ.３〜Ｎｏ.７が空洞Ａに対する空洞様反応領域を示す。測点Ｎｏ.１０〜Ｎｏ.１３が空洞Ｂに対する空洞様反応を示している。空洞様反応領域は、一対の低位相のピーク(空洞Ａでは測点Ｎｏ.３とＮｏ.７、空洞Bでは測点Ｎｏ.１０とＮｏ.１３)と、それに挟まれた低位相の測点(空洞Ａでは測点Ｎｏ.５、空洞Ｂでは測点Ｎｏ.１１)の領域として現れる。この一対の低位相のピークが図２（ｃ）の周縁効果の位置を示し、周縁効果に挟まれた低位相のピークが中心反応である。なお、空洞による空洞様反応では、中心反応、周縁効果のピークはほぼ等間隔に位置している。Ｎｏ.３、Ｎｏ.５、Ｎｏ.７の３測点間はほぼ等間隔であり、中央部の測点５、測点１１に空洞が存在すると推定できる。
【００５４】
また、空洞様反応の高位相のピークは比抵抗図では不明瞭で確認し辛いので、空洞別の空洞様反応領域を見定めるためには、これが明瞭に現われる等位相図を利用するのが好適である。
【００５５】
次に、空洞Ａ，Ｂの比抵抗の測定を行うと図４のような高比抵抗域が現れ、これが比抵抗における空洞様領域である。ＴＭモードでは、通常は高比抵抗のピークは空洞Ａ，Ｂの直上及び空洞Ａ，Ｂの両側方に帯状に現れる。各ピークが、比抵抗における中心反応と周縁効果を示す。
【００５６】
ところで、上述したように位相では低位相、および比抵抗では高比抵抗の中でピークを示す空洞様反応は、地電流電極と空洞の位置関係や、間隔及び幅寸法によって出現しない場合がある。例えば図４の空洞Ａ，Ｂはいずれも空洞直上の測点Ｎｏ.５とＮｏ.１１の高比抵抗のピーク値が欠落している。同様に、空洞両側方の周縁効果である低位相の一対のピークや、高比抵抗域の一対のピークの両者またはその一方が欠落する場合もある。
【００５７】
こうした場合、等位相図および等比抵抗図で空洞様反応のピークの出方が不明瞭であり、空洞の位置を正確に把握する必要がある場合には、ＴＭモード実現のため、地電流電極間隔，測点間隔，測点位置，送受信距離，受信方向等の測定条件を適宜調整して、空洞様反応における中心反応と周縁効果を明確にすればよい。両側に周縁効果が現れ、中央に中心反応が現れたときは、側線が空洞の直上で概ね直交して測定している状態であり、周縁効果の内側の低位相領域との境目付近が空洞内表面の位置を示すものと推定される。また、このとき空洞幅は、経験的に地電流電極１２ａ，１２ｂの間隔の１／２〜３／４程度であり、地電流電極１２ａ，１２ｂの間隔に基づいて空洞幅の大きさも推定できる。そして、空洞中心を簡易に判断する場合においては、周縁効果を示す一対のピーク間の中央位置を直ちに空洞中心とすればよい。
【００５８】
このように本発明の地中空洞探査装置と地中空洞探査方法は、空洞様反応を利用し、測定条件を繰り返し変えて測定することで図２（ｃ）に示すようなパターンを実現し、この状態で測定をすることにより空洞位置を推定するものである。なお、浅部の地下の空洞探査においては、目的の空洞の大きさ等の情報が別途得られることが多く、この繰り返し調査が面倒になる場合は多くない。
【００５９】
ところで実施の形態１においては、空洞中心の位置を精度よく推定するために、等位相図と等比抵抗図の両方を組み合わせて利用する。すなわち、まず、等位相図に現れた周縁効果を示す一対のピーク等を見つけて、これで空洞様反応の存在と位置に対する第１のチェックを行い、これを第２のチェックで確認するため、このピークと対応する測点とその内側領域の比抵抗のピークを抽出し、そのピークが３個であれば周縁効果を示すピークのほかに空洞直上の中心反応のピークが出現しているとして、中心反応のピークの位置を空洞中心と予測する。
【００６０】
なお、測定条件を繰り返し変えて測定しているときには、等位相図で１対のピークが得られなかったり、等位相図での一対のピークの内側に属する測点領域に、比抵抗のピーク値が１個しか現れず、図２（ｃ）のパターンを実現できない場合もある。このような場合は適宜地電流電極間隔を調整する等、測定条件を変えて再度測定する。測定条件を変えても図２（ｃ）のような適切な状態が得られない場合、例えばピークが４個以上現れたような場合には測定不能と判断する。
【００６１】
以上説明したように、本発明の実施の形態１の地中空洞探査装置と地中空洞探査方法によれば空洞位置を推定できるが、それを支える構成について説明する。高周波ＣＳＭＴ法による地中空洞探査方法は、ニアフィールド効果を受け易いため、従来から水平磁気双極子アンテナを使用することが多かった。その理由は、ニアフィールド効果は電磁波の透過深度と送受信間距離の関数であって、垂直磁気双極子アンテナではニアフィールド効果が、透過深度が送受信間距離の約１／５で現れ、水平磁気双極子アンテナでは透過深度が送受信間距離の約１／３で現れるからである。なお、水平磁気双極子アンテナはループアンテナを地表に垂直に設置し、磁場が地表に対して水平方向に形成されるアンテナであり、垂直磁気双極子アンテナは、ループアンテナを地表に水平に設置し、磁場が地表に対して垂直方向に形成されるアンテナである。
【００６２】
しかし、実施の形態１の地中空洞探査装置では、むしろニアフィールド効果を利用するために、敢えてニアフィールド効果が強く現われる垂直磁気双極子アンテナ等のアンテナを採用している。そしてアンテナはニアフィールド効果を発生し易いそれぞれ固有の範囲内に設置される。垂直磁気双極子アンテナであれば、上述したように、空洞の位置からおおむね想定される空洞の深度の３倍以上、５倍未満の距離に設置すればよい。
【００６３】
このニアフィールド効果は、低周波域で比抵抗や位相の測定値を増加させるという欠点があるため、従来から測定値を不明確にするものとして忌み嫌われているが、本発明でこれを利用するのは以下の理由による。すなわち、空洞部の比抵抗は無限大であり、空洞が無い所に比べニアフィールド効果が強いため、ニアフィールド領域で測定を行うと、空洞の有無で、ニアフィールド効果の強弱が生じ、その結果、比抵抗コントラストや位相コントラストが増幅すると考えられる。このため、通常とは逆にニアフィールド効果による比抵抗コントラストや位相コントラストの増幅を利用すると、空洞様反応のピークのコントラストを際立たせ、空洞位置の抽出をより明確に行うことが可能になる。また、垂直磁場双極子アンテナは、志向性が無いことや、設置が容易で取り扱い安いという利点があることから、この探査手法のアンテナとして最適である。
【００６４】
ところで、垂直磁気双極子アンテナは発信磁場の方向が受信磁場の方向に対し９０°のズレがある。このため垂直磁気双極子アンテナは水平磁気双極子アンテナと比較し、電磁波強度が弱く、高Ｓ／Ｎ比が得られ難い。従って垂直磁気双極子アンテナを採用しただけでは高周波分解能、Ｓ／Ｎ比が低く、実用性のある地中空洞探査装置は得られない。そこで、本実施の形態１においては、光ファイバ８を使って同期検波を行い、送信機部と受信機部の同期をとって高周波分解能を高め、高Ｓ／Ｎ比を実現したものである。なお、垂直磁気双極子アンテナのほかにも、垂直磁気双極子アンテナと同様にニアフィールド効果を増幅できるアンテナおよびその配置であれば使用することができる。
【００６５】
このように本実施の形態１においては、等比抵抗図でピークをより鋭敏に抽出するために、従来とは逆にニアフィールド効果を積極的に利用できる垂直磁気双極子アンテナを採用し、光ファイバ８を使って同期検波を行って高周波分解能を高め、高精度に空洞位置を決定するものである。
【００６６】
続いて、実施の形態１の地中探査方法において、空洞中心を推定して表示する一連の判断手順について説明する。図５は本発明の実施の形態１の地中探査方法で空洞位置を推定する判断手順のフローチャート、図６は本発明の実施の形態１の地中空洞探査装置と地中探査方法で空洞位置を推定して表示する動作のフローチャートである。図５において、図２（ａ）に示すように測線と送受信機部を配置して各測点の測定を行う。フーリエ変換したデータを入力して、１ｋＨｚ〜１００ｋＨｚの複数の周波数で、測点と比抵抗を計算するとともに測点と位相を計算し、これによって所定幅ごとに比抵抗値を分けて等高線を描いた等比抵抗図、所定幅ごとに位相を分けて等高線を描いた等位相図、さらに、１次元解析断面の演算を行い、比抵抗断面図を作成する（ｓｔｅｐ１）。
【００６７】
次いで、等位相図に現れる周縁効果を示す一対のピークを抽出できるか否かをチェックする（ｓｔｅｐ２）。一対のピークを抽出できない場合は、等位相図の低位相領域がほぼ対称か否かをチェックする（ｓｔｅｐ３）。ｓｔｅｐ３において低位相領域が対称の場合、おおむね対称軸が空洞の中心位置となり（ｓｔｅｐ４）、対称でない場合、ピークの位置が周縁効果の一方を示している（ｓｔｅｐ５）。ｓｔｅｐ４，５の後、ｓｔｅｐ２で一対のピークを抽出した場合とともに、一対のピークと対応した測点間の中心位置を算出すればよい（ｓｔｅｐ７）。ｓｔｅｐ１〜７において、等位相図だけでもデータが揃えば空洞中心を推定できることが分かる。
【００６８】
実施の形態１においては、さらに比抵抗図によってより正確な空洞位置を算出するため、位相図における一対のピークを示す測点間の領域で比抵抗のピークを抽出する（ｓｔｅｐ８）。ｓｔｅｐ８において、ピークが存在して、その数が３個以下であるか否かをチェックし（ｓｔｅｐ９）、ピークがないか４個以上のときはｓｔｅｐ２１へ進み、ピークが２個かどうかをチェックする（ｓｔｅｐ１０）。ピーク間の内側の地電流電極の電極間隔が空洞幅より大きいかどうかを判断し（ｓｔｅｐ１１）、大きい場合は２個のピークの中間の位置を空洞中心とし（ｓｔｅｐ１２）、電極間隔が空洞幅より小さい場合、２個のうち右側のピークが位相図における一対のピークの中央付近に存在するかどうかをチェックする（ｓｔｅｐ１３）。
【００６９】
すなわち、電極間隔と空洞幅の関係で、２個のピークが出現した場合は、周縁効果を示す一対のピークだけが現れた場合と、比抵抗のピークが１本しか現れない場合の変形で、空洞直上の比抵抗のピークと周縁効果の一方が現れた場合の２つがあり、ｓｔｅｐ１２で前者の場合を判定し、ｓｔｅｐ１３で後者の場合を判定する。ｓｔｅｐ１３において、右側のピークの位置が等位相図の一対のピークの中央付近に存在する場合、右側のピークの位置を空洞中心とし（ｓｔｅｐ１４）、右側のピークの位置が等位相図の一対のピークの中央付近に存在しない場合、左側のピークの位置を空洞中心とする（ｓｔｅｐ１５）。
【００７０】
ｓｔｅｐ１０において、２個のピークでない場合、３個のピークが存在するか否かをチェックし（ｓｔｅｐ１６）、ピークが１個の場合には、ｓｔｅｐ７の中央位置と一致しているかを確認する（ｓｔｅｐ１７）。一致しない場合はｓｔｅｐ２１へ進み、ｓｔｅｐ７の中央位置を空洞中心とするか、エラーとして終了する。ｓｔｅｐ１７において、ｓｔｅｐ７の中央位置と一致している場合は、このピーク位置が中心効果を示すものとして、これを空洞中心とする（ｓｔｅｐ１９）。また、ｓｔｅｐ１６において、ピークが３個の場合は、中央のピークの位置を空洞中心とする（ｓｔｅｐ１８）。
【００７１】
ｓｔｅｐ１２，１４，１５，１８，１９で空洞中心が定まると、この空洞中心が同一位置または大きくずれても隣接する測点にない場合（所定の範囲内の場合）は、１次元解析を行って（ｓｔｅｐ２２）、空洞の想定される深度と一致するか否かをチェックする（ｓｔｅｐ２３）。空洞の想定深度と一致しているときは空洞との推定が確定でき（ｓｔｅｐ２４）、終了する。ｓｔｅｐ２３において、深度が一致しない場合は別要因による反応があったものとして（ｓｔｅｐ２５）終了する。
【００７２】
続いて、図６に基づいて実施の形態１の地中空洞探査装置で空洞位置を推定して表示する動作について説明する。地電流電極間の距離を大きい間隔から徐々に小さくするなどし、測定条件を変えては、以下のｓｔｅｐ３１〜３３を繰返し、等位相図で適切なピークが得られたかをチェックし、等比抵抗図で３個のピークを求める。
【００７３】
ｓｔｅｐ３１において、フーリエ変換したデータを入力されると、演算手段２２の比抵抗算出手段２４と位相算出手段２５が、１ｋＨｚ〜１００ｋＨｚの複数の周波数で、測点と比抵抗を計算するとともに測点と位相を計算し、これによって所定幅ごとに比抵抗値を分けて等高線を描いた等比抵抗図、所定幅ごとに位相を分けて等高線を描いた等位相図を作成する。さらに、解析手段２６によって１次元解析断面の演算を行い、比抵抗断面図を作成する（ｓｔｅｐ３１）。
【００７４】
次いで、空洞位置算出手段２７が等位相図のデータに基づき、空洞様反応領域を示す周縁効果の一対のピークが得られたかチェックを行う（ｓｔｅｐ３２）。ここで、このピークの位置の計算を容易にするため、例えば各測点ごとにすべての周波数に関する位相の総和を計算し、「測点−位相総和グラフ」を作成し、このグラフにおける低位相のピークを抽出すればよい。
【００７５】
測定条件を変えたとき、地電流電極間隔が空洞幅に対して適切な場合（空洞幅が電極間隔の1／２〜３／４の場合）、空洞におけるピークは、図２（ｃ）における左右の図のように２パターンで示され、いずれもピークが出現する。この一対の周縁効果の現れた測点の内側の地電流電極の間が空洞の存在する範囲となり、その中心が空洞の中心となる。このうち、図２（ｃ）の左図の場合空洞中心に整った形で中心反応も現れるので、測定条件を変えてこの状態が得られるまで繰り返せばよい。このときより正確に空洞位置が推定できる。なお、実施の形態１の地中空洞探査装置では等位相図で一対のパターンを得て第１のチェックを行い、ニアフィールド効果で増幅された等比抵抗図のピークを使って第２の確認を行うから、等位相図でのピークは一対のピークを検出することで足りる。
【００７６】
そして、本地中空洞探査装置においては、この「測点−位相総和グラフ」による総和とピーク値の手法が自動化に馴染むため、これを採用しているが、これに限られるものではない。なお、ピーク値は土被比と地盤の比抵抗およびニアフィールド効果による増幅率の関数で決まるので、一律に閾値を設定するのは適当とはいえない。
【００７７】
さらに、位相の総和の測点に関する分布を対象として、空洞が存在する範囲となる分布形状の左右対称性をチェックしてもよい。これはピーク間の距離で判定できる。なお計算に際して、周波数として周縁効果を十分明確に示す高周波側の１０ｋＨｚ〜１００ｋＨｚを選び、この周波数だけを対象としてピークの位置を抽出するのも好適である。
【００７８】
ｓｔｅｐ３２において周縁効果を示す一対のピークが得られた場合、等比抵抗図のデータに基づき、３個の比抵抗のピークが適切な配置で存在するか確認する（ｓｔｅｐ３３）。計算を容易にするため、上述の位相と同様に、各測点ごとにすべての周波数における比抵抗の総和を計算し、ピーク値の間隔が図２（ｃ）の２パターンの何れかに対応すればよい。なお計算に際して、ニアフィールド効果により比抵抗コントラストが増幅された低周波領域となり、周縁効果を明確に示す１ｋＨｚ〜１０ｋＨｚを選び、この周波数を対象としてピークの位置を抽出するのも好適である。なお、比抵抗におけるピーク値も土被比と地盤の比抵抗およびニアフィールド効果による増幅率の関数で決まるので、一律に閾値を設定するのは適当とはいえない。
【００７９】
ｓｔｅｐ３３で抽出したピークが３個の場合、中央のピークの位置を空洞中心とし（ｓｔｅｐ３４）、終了する。ｓｔｅｐ３３において、３個以外の場合、電極間距離等の測定条件を変更し、測定条件の変更がすべて終わったか確認する（ｓｔｅｐ３５）。例えば、地電流電極間の距離を大きい間隔から徐々に小さくして測定している場合には、予定された測定条件の変更が終了したかを確認し、終了している場合は、測定不能である旨の表示を行ない（ｓｔｅｐ３６）、終了していない場合はｓｔｅｐ３１に戻って、測定を繰り返す。
【００８０】
以上説明したように本発明の実施の形態１の地中空洞探査装置と地中探査方法は、空洞の水平位置を推定するために、従来の電気探査のように数値解析を行って探査するのではなく、位相と比抵抗を相対的に評価することで、簡単な構成で明確に空洞の位置を推定することができる。そして、電気比抵抗が均一な地盤や、電線や電話線等から生じる電磁波ノイズが小さいところ、ノイズとなる地下埋設物が無い測定地においては、とくに顕著な結果を得ることができる。また、空洞位置の捕捉精度、水平分解能は、測点間隔を密にすることで容易に向上させることができる。さらに、比抵抗および位相のピーク値は土被比と地盤の比抵抗およびニアフィールド効果による増幅率の関数で決まるので、深く土被り比が大きい空洞におけるピーク値の低下に対して、送受信間距離を短くすることにより、ニアフィールド効果による増幅率を上げることができる。なお、実施の形態１の地中空洞探査装置と地中探査方法は、空洞そのものを確認しているものではないので、最終確認のためのボーリング等を行うことが望まれる。
【００８１】
（実施の形態２）
本発明の実施の形態２の地中空洞探査装置と地中探査方法は、調査地内において空洞位置の推定が困難な場合、闇雲に測線を設けるのではなく、簡便に測線を空洞上に空洞に対し、ほぼ直交方向に設置するまでの形態である。従って、実施の形態１の地中空洞探査装置と地中探査方法のいわば前提となるものである。
【００８２】
ＣＳＡＭＴ探査には、シャドウ効果という現象が知られている。これは、送受信機間に比抵抗異常がある場合に、その影が受信点の測定に影響を与える現象であり、発信位置を移動させるというシャドウ効果の排除方法が既存の研究で提案されている。本発明の実施の形態２は、シャドウ効果を排除するのでは無くこれを利用し、測線直下の空洞の有無、測線に対する未知の空洞の方向、測線と未知の空洞との距離、未知の空洞の向きを推定することで、測線を空洞上に空洞に対してほぼ直交方向に設置することができるものである。図７は本発明の実施の形態２における地中空洞探査方法を行うための地中空洞探査装置のシャドウ効果の説明図、図８は本発明の実施の形態２における地中空洞探査方法を行うための地中空洞探査装置の測線と送信部のＸ−Ｙ平面配置図である。
【００８３】
図７に示すように空洞調査のときのシャドウ効果は、探査測線の地下に空洞がなくても、空洞の延長上にシャドウ効果により空洞があるときと似た空洞反応が生じる。このシャドウ効果を除くために実施の形態１の地中空洞探査装置では、送信機部を測線の両翼に２箇所設けて両側から送信される電磁波を測定する。
【００８４】
実施の形態２の地中空洞探査装置は、測線と２つの送信機部を図８に示すように配置する。送信機部が２箇所になった点を除き、基本的に実施の形態１と同一であるから、図１〜図５を参照し、同一符号の説明は実施の形態１に譲ってここでは省略する。
【００８５】
実施の形態２においては、地下空洞に対して図８に示すように（併せて図１、図２も参照）、測線と送信機部Ａ，Ｂ、地電流電極、高周波用磁気センサ，低周波用磁気センサのインダクションコイルを配置する。想定される地下の空洞に直交すると思われる方向（Ｙ方向）に測線を設置し、測線の両翼（Ｘ方向）に５０ｍ〜１００ｍ測線から離した地点に送信機部Ａ，Ｂを置く。また、受信機部の地電流電極を測線方向に，インダクションコイルを測線の直交方向に設置する。地電流電極の間隔は、想定される空洞幅の１／２〜３／４程度とする。測定間隔は、地電流電極の間隔と同じかより短くする。1測点に付き送信機Ａ，Ｂからの信号を個別に測定する。
【００８６】
実施の形態２の地中空洞探査装置は、探査測線上において空洞の延長部に高比抵抗異常が検出されるため、シャドウ効果が出た測点で、図７の周縁効果の幅方向に測定方向を回転させるものである。この回転をさせることで、最大の比抵抗値と最小の位相値が得られた方向、あるいは最小の比抵抗値と最大の位相値が得られた方向の直交方向に空洞が存在することが分かる。この未知の空洞までの距離によってシャドウ効果の広がりが異なるため、この広がりを測定することでこの空洞までの距離を推定することができる。
【００８７】
ところで、測線の直下に空洞が存在することが未確定の場合、まず、送信機部Ａ，Ｂごとに図４に示すような比抵抗と等位相図を作成する。送信機部Ａ，Ｂの比抵抗と等位相図の両者で、空洞様反応が出た場合は測線の直下に空洞があると判断される。片方の発信機部で空洞様反応が生じた場合、シャドウ効果により空洞様反応が出た送信機と測線の間に空洞が存在する。
【００８８】
次に、強いシャドウ効果が現れた測点で、測定方向を少しずつ回転させ、最大の比抵抗値と最小の位相値が得られた方向と最小の比抵抗値と最大の位相値が得られた方向を記録する。理論的に前者の方向および後者の直交方向が未知の空洞がある方向を示す。このシャドウ効果の角度測定が周縁効果に対応するピークの２箇所で行えた場合、未知空洞の横幅を推定することにより、空洞までの距離が作図より求められる。
【００８９】
実施の形態２の地中空洞探査装置を使うことによって、送受信機部間の空洞の存在を推定することができ、今後の空洞調査の経費削減に大きく寄与することができる。測線をグリッド状に張り巡らせるのではなく、地中空洞の位置を推定することができる。
【００９０】
【実施例】
（実施例１）
鹿児島県鹿屋市周辺のシラス台地を探査地点として選んだ。この探査地においては、防空壕が多数存在する。図９は探査地における測線と測定番号の配置図、図１０（ａ）はＮｏ.１測線の等比抵抗図、図１０（ｂ）はＮｏ.１測線の等位相図、図１０（ｃ）は実測空洞の位置を表示したＮｏ.１測線の比抵抗２次元解析断面図である。
【００９１】
実施例１のＮｏ.１測線において、空洞Ａは図１０（ａ）（ｂ）から測点３と推定され、空洞Ｂは、図１０（ｂ）から選んだ測点８〜測点１１の範囲のうち、図１０（ａ）の測点９と推定される。空洞Ｃは、図１０（ｂ）において周縁反応のピーク位置は測点１３と測点１７であり、その中心を取って測点１５が空洞の中心と推定されるが、図１０（ａ）を見ると、その対称性から測点１５が空洞に関係していると見ることができる。空洞の実測位置は測点１４〜１５でこの推定と一致する。空洞Ｄは、図１０（ｂ）からも分かるように一対のピーク値の一方しかないため、断定は難しい。
【００９２】
しかし、実施例１の探査結果はきわめて正確に空洞位置を測定しているのが分かる。
【００９３】
（実施例２）
図１１（ａ）はＮｏ.２測線の等比抵抗図、図１１（ｂ）はＮｏ.２測線の等位相図、図１１（ｃ）は実測空洞の位置を表示したＮｏ.２測線の比抵抗２次元解析断面図である。
【００９４】
実施例２のＮｏ.２測線において、空洞Ａは、図１１（ｂ）から選んだ測点２〜測点５の範囲のうち中央の測点３と推定され、空洞Ｂは、図１１（ｂ）から選んだ測点９〜測点１０の範囲のうち、図１１（ａ）の測点１０と推定される。空洞Ｃは、図１１（ｂ）から選んだ測点１４〜測点１６の範囲のうち中央の測点１５と推定される。空洞の本坑と測線が低角度で交差しているためか、図１１（ａ）は測点１８まで高比抵抗が広がっている。空洞Ｄは、図１１（ｂ）から測点２０が周縁効果の一方で、空洞は終点側にあると見ることができる。このとき空洞が他のＡ，Ｂ，Ｃ孔とほぼ同じ大きさだと仮定すれば、空洞位置は測点２２と推定される。しかし、この測線における測点２０〜終点（測点２２）の空洞様反応は、空洞の延長部で生じたシャドウ効果である。すなわち、実施例２のように発信が１箇所だけで測定するとシャドウ効果が現れた場合に、殆ど測線直下にある空洞と判別できない。
【００９５】
このように、実施例２の探査結果は、空洞Ａ，Ｂで正確に推定できているが、空洞Ｃでは斜交部分で等比抵抗図（図１１（ａ））に異常が表れ、空洞Ｄは直上に測線が設けられていないにも係わらず、シャドウ効果により空洞様反応が現れていることが分かる。
【００９６】
【発明の効果】
本発明の地中空洞探査方法によれば、高周波の電磁波を送信する高周波ＣＳＭＴ法による探査を基本とし、これまでの探査手法では測定精度が悪かった地下１０ｍ以深（１０ｍ〜１００ｍ）を含む地下１ｍ〜１００ｍ程度の浅部から「大深度地下の公共的使用に関する特別措置法」の定める大深度地下の空洞位置を検出でき、垂直磁気双極子アンテナ等のアンテナを使用してニアフィールド効果により空洞の所在を示す空洞様反応を増幅することによって比抵抗や位相における空洞様反応を強調でき、送信機部と受信機部とで同期検波を行うためＳ／Ｎ比が上がり、精度よく安価に且つ自動的に探査することができる。また、シャドウ効果による比抵抗異常を判別することができ、さらには測線直下の空洞の有無，測線直下に無い空洞がある方向，測線直下に無い空洞の向き，およびその空洞までの距離を推定することができる。
【図面の簡単な説明】
【図１】 本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置の構成図
【図２】 （ａ）本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置の測線と送信部のＸ−Ｙ平面配置図、
（ｂ）（ａ）の地中空洞探査装置の測線とセンサと地下空洞のＸ−Ｙ−Ｚ配置
（ｃ）（ａ）の地中空洞探査装置の周縁効果の出現位置の説明図
【図３】 （ａ）比抵抗モデルの説明図
（ｂ）（ａ）の比抵抗の説明図
（ｃ）（ａ）の位相の説明図
【図４】 本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置が空洞位置に対して測定した等比抵抗と等位相の測定図
【図５】 本発明の実施の形態１の地中探査方法で空洞位置を推定する判断手順のフローチャート
【図６】 本発明の実施の形態１の地中空洞探査装置と地中探査方法で空洞位置を推定して表示する動作のフローチャート
【図７】 地中空洞探査装置のシャドウ効果の説明図
【図８】 本発明の実施の形態２における地中空洞探査装置の測線と送信部のＸ−Ｙ平面配置図
【図９】 探査地における測線と測定番号の配置図
【図１０】 （ａ）Ｎｏ.１測線の等比抵抗図
（ｂ）Ｎｏ.１測線の等位相図
（ｃ）実測空洞の位置を表示したＮｏ.１測線の比抵抗２次元解析断面図
【図１１】 （ａ）Ｎｏ.１測線の等比抵抗図
（ｂ）Ｎｏ.２測線の等位相図
（ｃ）実測空洞の位置を表示したＮｏ.２測線の比抵抗２次元解析断面図
【図１２】 図４と同測線における水平磁気双極子アンテナによる等比抵抗と等位相の測定図
【符号の説明】
１ 発振部
２ 周波数逓倍部
２ａ スイッチング部
３ ドライバ部
４ａ 高周波アンテナ
４ｂ 低周波アンテナ
５ 送信側制御部
６ 計時手段
７ 記憶部
８ 光ファイバ
８ａ，８ｂ 信号変換部
９ 同期部
１１ａ 高周波用磁気センサ
１１ｂ 低周波用磁気センサ
１２ａ，１２ｂ 地電流電極
１３ａ，１３ｂ プリアンプ
１４ａ，１４ｂ バッファアンプ
１５，１６ バンドパスフィルタ
１７，１８ 増幅部
１９ フーリエ変換部
１９ａ 高周波変換部
【発明の属する技術分野】
本発明は、高周波ＣＳＭＴ法により地中空洞の位置を高精度に推定する地中空洞探査方法に関する。
【０００２】
【従来の技術】
従来、地表から２ｍ以上の浅部の探査においては地中レーダーが有効とされている。しかし、この探査方法は信号の減衰が大きいので、深度１０ｍを超えると事実上探査ができなくなるものであった。地下の空洞を探査するときにも２ｍ〜１０ｍが限界であった。
【０００３】
電気探査方法にはこのような制限がなく、これを使用して探査することも考えられる。しかし、電気探査方法は深度が深くなると水平分解能が低下する傾向がある。例えば、地下の状態を推定する有力な手法として、電気比抵抗を測定してそのコントラストを測定する比抵抗調査法がある。中でも、直流比抵抗法を使用した高密度電気探査法（比抵抗影像法）がもっとも一般的で普及している。この直流比抵抗法は、探査深度を深くするとき電極間隔を広げる必要があり、これによって水平分解能が低下してしまう。しかも、測定時の作業負担が大きく、探査費用も嵩むものであった。
【０００４】
また、このような比抵抗調査法の一種である電磁探査法にＭＴ法（地磁気地電流探査）と可聴領域（Ａ）でコントロール（Ｃ）された電磁波を発信源（Ｓ；ソース）から発信し測定するＣＳＡＭＴ法があるが、これらは周波数を変化させて深部を探査するため電極間隔を狭く保て、測定点が信号源から十分離れている場合に電磁場を平面波とみなして測定するため、条件を満たせば比較的高い精度が期待できる。電磁波を利用するため作業負担、費用に関しても優れている。
【０００５】
このうち、ＭＴ法は地球磁気圏や雷放電により発生する自然電磁場の１０−３Ｈｚ〜２５０Ｈｚの周波数を測定するものであり、自然電磁場を利用して地下の比抵抗構造を分析するため、データが得られなかったり、ノイズに弱く、このような意味で精度が高いとはいえなかった。
【０００６】
これに対して、人工信号源からの電磁波を測定するＣＳＡＭＴ法はノイズに強く、測定費用も安価であり、周波数を変化させて探査するので電気探査、直流比抵抗法と比較して深部における水平分解能が高いという利点を有している。しかし、低周波数域ではニアフィールド効果を受け易いという欠点も有している。このニアフィールド効果は、このような低周波域においては平面波の仮定が崩れて見掛比抵抗が増加し、直流電界が強く、測定点直下の正確な比抵抗や位相を得るのが難しくなるものである。
【０００７】
このニアフィールド効果を減らす１つの手段として、従来、水平磁気双極子アンテナが採用されてきた。しかし、これによってＣＳＡＭＴ法に伴うニアフィールド効果が根本的に解消されるようなものではない。また、深度１０ｍ〜１００ｍの浅部に存在する高比抵抗の空洞を探査するのは、高周波情報を利用していない従来のＣＳＡＭＴ法では事実上困難と考えられてきた。
【０００８】
そこでニアフィールドの問題に対しては、ＣＳＭＴ法の特徴を活かしてニアフィールド効果を除く方法が提案された（特許文献１参照）。これは、ＣＳＭＴ法による探査とＭＴ法による探査を重複して行い、ＣＳＭＴ法によって得たデータのうち、ニアフィールド効果が発生する周波数領域を、ＭＴ法によるデータによって置き換えて地下構造を解析するものである。しかしながら、この方法ではＭＴ法のノイズに弱いという欠点はそのまま残るものであり、また、数百〜数千ｍの大深度の調査に適応されるものであるため、深度１０ｍ〜１００ｍの浅部の調査にＣＳＡＭＴ法は利用されていなかった。
【０００９】
そこで本発明者らは、ＣＳＭＴ法の水平分解能が良いという特徴を活かして深度１００ｍまでの浅部を探査するため、高周波ＣＳＭＴ法による探査装置を提案した（非特許文献１参照）。１ｋＨｚ〜１００ｋＨｚ（１６測点）で測定を行うものである。アンテナは１〜１０ｋＨｚの低周波と、１０ｋＨｚ〜１００ｋＨｚの高周波の２帯域で異なったループアンテナを使っている。高周波を用いることにより、従来のＣＳＭＴ法と比較して、（１）装置、とくに送信部が小型、軽量化され、（２）探査時間が１周波数について１０秒、１測点について１６０秒と、大幅な短縮が可能であり、（３）水平方向の分解能が高く、操作性もよく、地形の制約を受けることが少なく、探査費用が低減できるものであった。そして、非特許文献１で提案した探査装置で分析した比抵抗構造は、高密度電気探査法とほぼ同じになるという結果を得ている。
【００１０】
【特許文献１】
特開平５−３２３０３８号公報
【非特許文献１】
半田，外３名，「高周波ＣＳＭＴ装置の開発」，物理探査，物理探査学会，５５巻，３号，ｐ．１９９−２０６
【００１１】
【発明が解決しようとする課題】
以上説明したように、従来のＣＳＭＴ法は低周波数域ではニアフィールド効果を受け易いという問題を有していた。ニアフィールド領域では、直流電界が強くなり、その結果、見掛比抵抗（以下、比抵抗）が増加し、解析による測定点直下の「真の比抵抗」の算出ができなくなる。このため、従来は如何にニアフィールド効果を低減させるかが、探査性能の向上に重要と考えられていた。
【００１２】
図１２は従来から、ニアフィールド効果が少ないので多く使用されてきた水平磁気双極子アンテナを使用したＣＳＭＴ法による比抵抗と位相の測定図である。水平磁気双極子アンテナを使用したＣＳＭＴ法によって空洞探査した場合の事例を示す。この事例では、同じ発信地点で水平磁気双極子アンテナを使用していることから、ニアフィールド効果は垂直磁場双極子アンテナに比べ少ないと考えてよい。比抵抗に関しては図１２に示すように、周波数１．４ｋＨｚ〜９７ｋＨｚで１０２Ω・ｍ以下で変化に乏しく、１３ｋＨｚで僅かに空洞部で高比抵抗となっているに過ぎない。同様に、位相に関しても図１２に示すように、周波数１．４ｋＨｚ〜４５ｋＨｚでは変化の乏しい位相６０°〜８０°を示すだけであり、周波数４５ｋＨｚ〜９７ｋＨｚの範囲でも空洞を明確に示唆するような情報を含まないものである。このように、比抵抗および位相の地盤応答だけでは、高周波ＣＳＭＴ探査機を用いても空洞を明確に識別できない。
【００１３】
また、非特許文献１で説明した本発明者らの高周波ＣＳＭＴ法による探査装置は、深さ１００ｍ以下の浅部探査であるためにループアンテナを採用し、ＰＬＬ回路を使用してフーリエ変換するが、参照信号と送信信号との差で出力が変動してしまい、高周波分解能が悪く、高Ｓ／Ｎ比を実現するのが難しいものであった。そして高密度電気探査法と同程度の探査が可能になるという比較結果は得られたが、上述の出力変動でＳ／Ｎ比が上がらず、分析した比抵抗構造は粗くて、空洞位置を正確に予測することなど期待できないものであった。
【００１４】
空気が充満する空洞においては地盤より際立って比抵抗が大きくなるが、非特許文献１は、こうした空洞の場合について未研究で、空洞の比抵抗が示す挙動、位相が示す挙動のいずれもがまったく未知であった。そして、こうした空洞様反応が得られたとしても、空洞位置をどのように分析し、判定するのか、その手法はまったく未解決であった。
【００１６】
そこで本発明は、浅部の空洞の位置を精度よく安価に探査することができる地中空洞探査方法を提供することを目的とする。
【００１７】
【課題を解決するための手段】
以上説明した課題を解決するために本発明の地中空洞探査方法は、地中探査のため高周波信号を電磁波として送信するとともに、一対の地電流電極で測線上の測点の電位差を検出するとともに該測線と直交する方向に配設されたインダクションコイルにより該測点の磁場を測定し、該測点における比抵抗と位相を演算する高周波ＣＳＭＴ法による地中空洞探査方法であって、空洞の所在を示す空洞様反応をニアフィールド効果によって増幅するための垂直磁気双極子アンテナを空洞から所定距離おいた位置に配置して高周波信号を送信し、探査領域における周波数と測点に関する等比抵抗データを取得して、該等比抵抗データの中からニアフィールド効果により増幅された空洞様反応領域を抽出し、該空洞様反応領域の測線方向幅の中央を空洞中心として推定することを特徴とする。
【００１８】
本発明は、この構成により、ニアフィールド効果を避けるのではなく、これを利用し、ニアフィールド効果による空洞様反応を増幅させて可視化することによって、浅部の空洞の位置を精度よく安価に且つ自動的に探査することができる。
【００１９】
【発明の実施の形態】
本発明の第１の発明は、地中探査のため高周波信号を電磁波として送信するとともに、一対の地電流電極で測線上の測点の電位差を検出するとともに該測線と直交する方向に配設されたインダクションコイルにより該測点の磁場を測定し、該測点における比抵抗と位相を演算する高周波ＣＳＭＴ法による地中空洞探査方法であって、空洞の所在を示す空洞様反応をニアフィールド効果によって増幅するための垂直磁気双極子アンテナを空洞から所定距離おいた位置に配置して高周波信号を送信し、探査領域における周波数と測点に関する等比抵抗データを取得して、該等比抵抗データの中からニアフィールド効果により増幅された空洞様反応領域を抽出し、該空洞様反応領域の測線方向幅の中央を空洞中心として推定することを特徴とする地中空洞探査方法であり、高周波の電磁波を送信する高周波ＣＳＭＴ法による探査を基本とし、これまでの探査手法では測定精度が悪かった地下１０ｍ以深（１０ｍ〜１００ｍ）を含む地下１ｍ〜１００ｍ程度の浅部から「大深度地下の公共的使用に関する特別措置法」の定める大深度地下の空洞位置を検出でき、垂直磁気双極子アンテナを使用してニアフィールド効果により空洞の所在を示す空洞様反応を増幅することによって比抵抗における空洞様反応を強調でき、高Ｓ／Ｎ比で、精度よく安価に探査することができる。空洞様反応領域の中から、ニアフィールド効果でコントラストが強調された等比抵抗データの空洞様反応を抽出するから、空洞中心の判別が正確になる。
【００２８】
本発明の第２の発明は、第１の発明に従属する発明であって、探査領域における周波数と測点に関する等比抵抗データを取得して、該空洞様反応領域の測線方向幅の中央を空洞中心として推定するのに代えて、探査領域における周波数と測点に関する等位相データの中から空洞様反応のピークを含む低位相の空洞様反応領域を抽出し、この空洞様反応領域に対応する等比抵抗データの中から比抵抗のピークを抽出して空洞中心を推定することを特徴とする地中空洞探査方法であり、等位相データからピークを抽出して空洞様反応領域の存在を確認し、さらにニアフィールド効果でコントラストが強調された等比抵抗データの周縁効果と中心反応示す空洞様反応のピークを抽出することにより、空洞中心の判別がさらに正確になる。
【００３２】
本発明の第３の発明は、第１又は第２の発明に従属する発明であって、測線の両側から電磁波を送信し、該測線の直下にある空洞と測線の直下に無い空洞とを、シャドウ効果による空洞様反応に近似した比抵抗または位相異常の有無から判別することを特徴とする地中空洞探査方法であり、側線の直下の空洞と、直下にない空洞とを判別できる。
【００３７】
（実施の形態１）
以下、本発明の実施の形態１における地中空洞探査方法について説明する。実施の形態１は、空洞上に測線が空洞に対し、ほぼ直交方向に設置されている場合に空洞位置を把握するものである。なお、実施の形態１の地中空洞探査装置は、測点間距離等の測定条件を変更し、また送受信方向を選択して、繰り返し測定を行うことにより、この測定条件を満たす状態を実現し、実際の測定を行うものである。図１は本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置の構成図、図２（ａ）は本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置の測線と送信部のＸ−Ｙ平面配置図、図２（ｂ）は（ａ）の地中空洞探査方法を行うための地中空洞探査装置の測線とセンサと地下空洞のＸ−Ｙ−Ｚ配置図、図２（ｃ）は（ａ）の地中空洞探査装置における空洞様反応の出現位置の説明図、図３（ａ）は比抵抗モデルの説明図、図３（ｂ）は（ａ）の比抵抗の説明図、図３（ｃ）は（ａ）の位相の説明図、図４は本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置が空洞位置に対して測定した見掛比抵抗と位相の測定図である。
【００３８】
図１において、１は地中空洞探査装置から信号波を送信するために発振を行う発振部、２は発振部１で発振した１００Ｈｚの基本周波数が入力されると１ｋＨｚ〜１００ｋＨｚの周波数帯を対数的に１６分割した周波数を得るため基本周波数の逓倍を行う周波数逓倍部、２ａは周波数逓倍部２からの信号で±１３Ｖの電圧で矩形波信号を形成するスイッチング部、３は矩形波信号を増幅するドライバ部、４ａはループコイルとコンデンサで共振回路を構成した後述するニアフィールド効果を増幅するための高周波アンテナ、４ｂは同じく、ループコイルとコンデンサで共振回路を構成したニアフィールド効果を増幅するための低周波アンテナである。高周波アンテナ４ａは９７ｋＨｚ，７３ｋＨｚ，５５ｋＨｚ，４２ｋＨｚ，３１ｋＨｚ，２４ｋＨｚ，１８ｋＨｚ，１３ｋＨｚの８つの周波数を送信し、低周波アンテナ４ｂは１０ｋＨｚ，７．６ｋＨｚ，５．８ｋＨｚ，４．３ｋＨｚ，３．３ｋＨｚ，２．５ｋＨｚ，１．８ｋＨｚ，１．４ｋＨｚの８つの周波数を送信する。そして、この高周波アンテナ４ａと低周波アンテナ４ｂはニアフィールド効果が強く現われる特性をもつ垂直磁気双極子アンテナを構成するものである。なお、周波数がこれらの周波数に限られないのは当然である。
【００３９】
５は中央処理装置（以下、ＣＰＵ）等から構成され、制御プログラムを読み込んで機能実現手段として機能を実現する地中空洞探査装置の送信側制御部である。６は１６の周波数の信号波を１０秒ごとに切り替えるための計時手段、７は制御プログラム等を記憶した記憶部、８は発振部１が発振した１００Ｈｚの基本周波数の信号を地中空洞探査装置の受信側に送って同期をとるための光ファイバ、８ａは発振部１の電気信号を光信号に変換する信号変換部、８ｂは光信号を電気信号に変化する信号変換部、９は発振部１と信号変換部８ａから構成される地中空洞探査装置の同期部である。
【００４０】
１１ａは共振周波数９７ｋＨｚのインダクションコイルＨから構成される高周波用磁気センサ、１１ｂは共振周波数９．７ｋＨｚのインダクションコイルＬから構成される低周波用磁気センサ、１２ａ，１２ｂはステンレス製の地電流電極、１３ａ，１３ｂはプリアンプ、１４ａ，１４ｂは入力インピーダンスを調整するバッファアンプである。１５，１６はバンドパスフィルタ、１７，１８は増幅部、１９はフーリエ変換部、１９ａはＰＬＬ回路を使って余弦波成分と正弦波成分でフーリエ変換する高周波変換部、１９ａはＰＬＬ回路を使って余弦波成分と正弦波成分でフーリエ変換する低周波変換部、２０はＡ／Ｄ変換部である。
【００４１】
２１はＣＰＵ等から構成され、制御プログラムを読み込んで機能実現手段として演算機能，システム制御機能等を実現する地中空洞探査装置の制御部、２２はデジタル化された測定値から各演算を行う機能実現手段としての演算手段、２３は地中空洞探査装置のシステム制御を行う機能実現手段としてのシステム制御手段である。システム制御手段２３は高周波用磁気センサ１１ａ，低周波用磁気センサ１１ｂの切り替え、増幅部１７，１８のゲインやバンドパスフィルタ１５，１６の中心周波数の切り替え等を行う。
【００４２】
２４はフーリエ変換部１９からの出力を使って比抵抗を演算する比抵抗算出手段、２５はフーリエ変換部１９の出力から位相を演算する位相算出手段、２６は比抵抗構造を解析するためボスティック法等による１次元解析断面の演算を行う解析手段である。後述するニアフィールド効果を避けて解析するため、高周波側のデータだけで解析するのも好適である。1次元インバージョン解析では深度誤差が大きい場合、空洞上部の地盤の緩みや崩壊、あるいは空洞とは無関係の地質要因による異常値を示す可能性があるが、解析手段２６でさらに２次元インバージョンを行うことにより、空洞位置の解析精度を上げることができる。
【００４３】
２７は位相算出手段２５が算出した周波数と測点の等位相図（同一位相の周波数と測点を示すもので、本発明の等位相データ）に基づいて、低位相の空洞様反応領域から中心反応と周縁効果を示す一対の空洞様反応のピークを抽出し、また、比抵抗算出手段２４によって作成された比抵抗の等比抵抗図（同一比抵抗の周波数と測点を示すもので、本発明の等比抵抗データ）に基づいて、高比抵抗の空洞様反応領域から中心反応と周縁効果を示す空洞様反応のピークを抽出し、このピークから地下空洞の位置を算出する空洞位置算出手段である。この空洞様反応領域と、それに含まれる中心反応と周縁効果の詳細については後述する。
【００４４】
実施の形態１の空洞位置算出手段２７は、位相に関する周縁効果のピークをチェックし、さらに高比抵抗の空洞様反応領域から中心反応と周縁効果を示す空洞様反応のピークを抽出して空洞中心を推定する。このとき、等位相図で周縁効果を示す一対のピークからその中央位置を計算し、この中央位置を簡易に空洞の中心位置として推定することも可能である。なお、中心反応と周縁効果を示す空洞様反応のピークは最適な測定条件下において整った形態で出現するが（後述の図２（ｃ）参照）、状況によっては顕著には出現しない場合もある。このため実施の形態１の地中空洞探査装置においては、地電流電極間の距離を大きい間隔から徐々に小さくして測定し、また受信方向、その距離を変更し、あるいはその他の測定条件を変更して、繰り返し測定を行うことにより最適な測定条件を実現している。最適条件を満たしているか否かは、中心反応と周縁効果のピークが整ったパターンで得られたことで判断する。
【００４５】
以上説明した実施の形態１の地中空洞探査装置は、基本的に、高周波アンテナ４ａ、低周波アンテナ４ｂからそれぞれ送信信号を送信する送信機部と、高周波用磁気センサ１１ａと低周波用磁気センサ１１ｂで受信する受信機部とから構成される。なお、実施の形態２で説明するように送信機部は受信機部と同期をとって２台、あるいは複数設けることができる。
【００４６】
次に、実施の形態１の地中空洞探査装置と地中空洞探査方法で空洞位置を決定できる原理を説明する。まず、実施の形態１の地中空洞探査装置で探査を行うときの測線の設け方、測定の方法について説明する。本発明が探査する典型的な場合は、地下に空洞が存在することが予測される調査地で、どこに空洞があるかを探査する場合等であり、測定条件を変更しての測定はそれほど面倒になることはない。なお、空気が充満した空洞を探査する事例で説明するが、地下水が溜まった空洞や土の陥没で埋まった空洞は、空洞に溜まった地下水や充填土と周辺地質の比抵抗の違いによって探査するため、空洞様反応の感度が低くなる。
【００４７】
実施の形態１では、地下空洞に対して図２（ａ）（ｂ）に示すように測線と送信機Ａ、地電流電極１２ａ，１２ｂ、高周波用磁気センサ１１ａ，低周波用磁気センサ１１ｂのインダクションコイルＨ，Ｌを配置して測定を行う。想定される地下の空洞に直交すると思われる方向（Ｙ方向）に測線を設置し、測線の側方（Ｘ方向）に５０ｍ〜１００ｍ測線から離した地点に送信機Ａを置く。また、受信機の地電流電極１２ａ，１２ｂを測線方向に，インダクションコイルＨ，Ｌを測線の直交方向に設置する。地電流電極１２ａ，１２ｂの間隔は、想定される空洞幅の１／２〜３／４程度とする。測定間隔は、地電流電極１２ａ，１２ｂの間隔と同一か、短くする。なお、上述したように、中心反応と周縁効果のピークが整った形で得られたときが最も適切な測定条件のときであるが、この状態は地電流電極１２ａ，１２ｂと空洞間隔を上記の通り調整することで実現できる。1測点に付き送信機Ａからの信号を測定する。
【００４８】
上述のように測線とインダクションコイルを配置したときは、図２（ａ）（ｂ）のように空洞は電場方向に対して直交する方向（磁場方向）となる。このような電場方向に対して直交する方向または高角度で斜交する状態（以下、ＴＭモードという）においては、空洞と地盤との境界では比抵抗が不連続となり、この不連続点で比抵抗が極端な変化を示す。なお、これが逆のものがＴＥモードである。図３（ａ）に示した比抵抗モデルは、比抵抗ρ１＝１０Ωｍの媒体と比抵抗ρ２＝１００Ωｍの媒体が接しているものである。図３（ｂ）（ｃ）の×点は、比抵抗ρ１＝１０Ωｍの媒体と比抵抗ρ２＝１００Ωｍの媒体が接しているときの、ＴＭモードにおける境界での比抵抗、位相の変化を示している。○点は、空洞を電場方向と平行に向けた場合の比抵抗、位相の変化である。ＴＭモードにおいては２つの媒体の境界面で比抵抗が鋭敏な変化をし、位相は高比抵抗領域に入ると低下している。
【００４９】
これを地中の空洞の場合で説明すると、図２（ａ）（ｃ）に示すように、ＴＭモードの空洞と地盤の境界において比抵抗は地盤側では非常に小さく空洞側で非常に大きくなり、コントラストが際立つものとなる。同様に位相においても、地盤側境界で大きなピーク値をとって空洞では低位相を示すことを示す。
【００５０】
しかし、実際の空洞においては少し複雑になる。すなわち、図３の単純な２次元構造の場合と異なり、図２（ｃ）に示すように空洞の直上の測点で比抵抗は高く、位相は低くなるが、空洞側壁を跨がない最も空洞に近い測点でも同様に、比抵抗は高く、位相は低くなる。そして、この高比抵抗および低位相の異常値は空洞から遠ざかるほど小さくなる。この現象は、現場試験と数値解析によるモデル実験で確認することができ、ＴＭモードにおける空洞に対するＭＴ応答の特徴として定義できるものである。そして図２（ｃ）の状態は測線が空洞上を直交して交差した場合の測点、地電流電極間隔および空洞位置に関係する空洞様反応のピークパターンで、これらの位置関係が悪いとピークが隠れてしまうことがある。
【００５１】
このＭＴ応答の特徴を挙げると、（１）空洞直上は高比抵抗と低位相のピークとなる。本明細書においてはこれを「中心反応」という。また、（２）空洞側壁を跨がない最も空洞に近い測点は高比抵抗と低位相のピークとなる。本明細書においてはこれを「周縁効果」という。（３）周縁効果は空洞から離れるほど小さくなる。本明細書においては「中心反応」「周縁効果」を一括して「空洞様反応」という。とくに、中心反応と空洞両翼に生じる一対の周縁効果が生じている範囲を本明細書においては「空洞様反応域」と定義する。ここで、「様」を付加したのは実際の空洞調査において，この特徴を持ったＭＴ応答が現れたとしても、直ちに空洞と断定できないので、「空洞あるいは空洞に似た反応」のニュアンスを残すためである。
【００５２】
そこで、本発明の地中空洞探査装置と地中空洞探査方法の基本となるこの空洞様反応の出現の仕方を、実測した等位相図と等比抵抗図で説明する。図４は、地中空洞Ａが測点５、地中空洞Ｂが測点１１に存在するもので、図１２と同じ地中空洞Ａ，地中空洞Ｂを図２（ａ）（ｂ）（ｃ）に示すような測点で同一の条件で測定したものである。後述の実施例１，２と同一の測定対象で測定したものである。空洞の天井までの深度は１３．７ｍ、空洞の高さは２．２ｍで、この比を意味する土被比は約６．３である。
【００５３】
図４において、測点Ｎｏ.３〜Ｎｏ.７が空洞Ａに対する空洞様反応領域を示す。測点Ｎｏ.１０〜Ｎｏ.１３が空洞Ｂに対する空洞様反応を示している。空洞様反応領域は、一対の低位相のピーク(空洞Ａでは測点Ｎｏ.３とＮｏ.７、空洞Bでは測点Ｎｏ.１０とＮｏ.１３)と、それに挟まれた低位相の測点(空洞Ａでは測点Ｎｏ.５、空洞Ｂでは測点Ｎｏ.１１)の領域として現れる。この一対の低位相のピークが図２（ｃ）の周縁効果の位置を示し、周縁効果に挟まれた低位相のピークが中心反応である。なお、空洞による空洞様反応では、中心反応、周縁効果のピークはほぼ等間隔に位置している。Ｎｏ.３、Ｎｏ.５、Ｎｏ.７の３測点間はほぼ等間隔であり、中央部の測点５、測点１１に空洞が存在すると推定できる。
【００５４】
また、空洞様反応の高位相のピークは比抵抗図では不明瞭で確認し辛いので、空洞別の空洞様反応領域を見定めるためには、これが明瞭に現われる等位相図を利用するのが好適である。
【００５５】
次に、空洞Ａ，Ｂの比抵抗の測定を行うと図４のような高比抵抗域が現れ、これが比抵抗における空洞様領域である。ＴＭモードでは、通常は高比抵抗のピークは空洞Ａ，Ｂの直上及び空洞Ａ，Ｂの両側方に帯状に現れる。各ピークが、比抵抗における中心反応と周縁効果を示す。
【００５６】
ところで、上述したように位相では低位相、および比抵抗では高比抵抗の中でピークを示す空洞様反応は、地電流電極と空洞の位置関係や、間隔及び幅寸法によって出現しない場合がある。例えば図４の空洞Ａ，Ｂはいずれも空洞直上の測点Ｎｏ.５とＮｏ.１１の高比抵抗のピーク値が欠落している。同様に、空洞両側方の周縁効果である低位相の一対のピークや、高比抵抗域の一対のピークの両者またはその一方が欠落する場合もある。
【００５７】
こうした場合、等位相図および等比抵抗図で空洞様反応のピークの出方が不明瞭であり、空洞の位置を正確に把握する必要がある場合には、ＴＭモード実現のため、地電流電極間隔，測点間隔，測点位置，送受信距離，受信方向等の測定条件を適宜調整して、空洞様反応における中心反応と周縁効果を明確にすればよい。両側に周縁効果が現れ、中央に中心反応が現れたときは、側線が空洞の直上で概ね直交して測定している状態であり、周縁効果の内側の低位相領域との境目付近が空洞内表面の位置を示すものと推定される。また、このとき空洞幅は、経験的に地電流電極１２ａ，１２ｂの間隔の１／２〜３／４程度であり、地電流電極１２ａ，１２ｂの間隔に基づいて空洞幅の大きさも推定できる。そして、空洞中心を簡易に判断する場合においては、周縁効果を示す一対のピーク間の中央位置を直ちに空洞中心とすればよい。
【００５８】
このように本発明の地中空洞探査装置と地中空洞探査方法は、空洞様反応を利用し、測定条件を繰り返し変えて測定することで図２（ｃ）に示すようなパターンを実現し、この状態で測定をすることにより空洞位置を推定するものである。なお、浅部の地下の空洞探査においては、目的の空洞の大きさ等の情報が別途得られることが多く、この繰り返し調査が面倒になる場合は多くない。
【００５９】
ところで実施の形態１においては、空洞中心の位置を精度よく推定するために、等位相図と等比抵抗図の両方を組み合わせて利用する。すなわち、まず、等位相図に現れた周縁効果を示す一対のピーク等を見つけて、これで空洞様反応の存在と位置に対する第１のチェックを行い、これを第２のチェックで確認するため、このピークと対応する測点とその内側領域の比抵抗のピークを抽出し、そのピークが３個であれば周縁効果を示すピークのほかに空洞直上の中心反応のピークが出現しているとして、中心反応のピークの位置を空洞中心と予測する。
【００６０】
なお、測定条件を繰り返し変えて測定しているときには、等位相図で１対のピークが得られなかったり、等位相図での一対のピークの内側に属する測点領域に、比抵抗のピーク値が１個しか現れず、図２（ｃ）のパターンを実現できない場合もある。このような場合は適宜地電流電極間隔を調整する等、測定条件を変えて再度測定する。測定条件を変えても図２（ｃ）のような適切な状態が得られない場合、例えばピークが４個以上現れたような場合には測定不能と判断する。
【００６１】
以上説明したように、本発明の実施の形態１の地中空洞探査装置と地中空洞探査方法によれば空洞位置を推定できるが、それを支える構成について説明する。高周波ＣＳＭＴ法による地中空洞探査方法は、ニアフィールド効果を受け易いため、従来から水平磁気双極子アンテナを使用することが多かった。その理由は、ニアフィールド効果は電磁波の透過深度と送受信間距離の関数であって、垂直磁気双極子アンテナではニアフィールド効果が、透過深度が送受信間距離の約１／５で現れ、水平磁気双極子アンテナでは透過深度が送受信間距離の約１／３で現れるからである。なお、水平磁気双極子アンテナはループアンテナを地表に垂直に設置し、磁場が地表に対して水平方向に形成されるアンテナであり、垂直磁気双極子アンテナは、ループアンテナを地表に水平に設置し、磁場が地表に対して垂直方向に形成されるアンテナである。
【００６２】
しかし、実施の形態１の地中空洞探査装置では、むしろニアフィールド効果を利用するために、敢えてニアフィールド効果が強く現われる垂直磁気双極子アンテナ等のアンテナを採用している。そしてアンテナはニアフィールド効果を発生し易いそれぞれ固有の範囲内に設置される。垂直磁気双極子アンテナであれば、上述したように、空洞の位置からおおむね想定される空洞の深度の３倍以上、５倍未満の距離に設置すればよい。
【００６３】
このニアフィールド効果は、低周波域で比抵抗や位相の測定値を増加させるという欠点があるため、従来から測定値を不明確にするものとして忌み嫌われているが、本発明でこれを利用するのは以下の理由による。すなわち、空洞部の比抵抗は無限大であり、空洞が無い所に比べニアフィールド効果が強いため、ニアフィールド領域で測定を行うと、空洞の有無で、ニアフィールド効果の強弱が生じ、その結果、比抵抗コントラストや位相コントラストが増幅すると考えられる。このため、通常とは逆にニアフィールド効果による比抵抗コントラストや位相コントラストの増幅を利用すると、空洞様反応のピークのコントラストを際立たせ、空洞位置の抽出をより明確に行うことが可能になる。また、垂直磁場双極子アンテナは、志向性が無いことや、設置が容易で取り扱い安いという利点があることから、この探査手法のアンテナとして最適である。
【００６４】
ところで、垂直磁気双極子アンテナは発信磁場の方向が受信磁場の方向に対し９０°のズレがある。このため垂直磁気双極子アンテナは水平磁気双極子アンテナと比較し、電磁波強度が弱く、高Ｓ／Ｎ比が得られ難い。従って垂直磁気双極子アンテナを採用しただけでは高周波分解能、Ｓ／Ｎ比が低く、実用性のある地中空洞探査装置は得られない。そこで、本実施の形態１においては、光ファイバ８を使って同期検波を行い、送信機部と受信機部の同期をとって高周波分解能を高め、高Ｓ／Ｎ比を実現したものである。なお、垂直磁気双極子アンテナのほかにも、垂直磁気双極子アンテナと同様にニアフィールド効果を増幅できるアンテナおよびその配置であれば使用することができる。
【００６５】
このように本実施の形態１においては、等比抵抗図でピークをより鋭敏に抽出するために、従来とは逆にニアフィールド効果を積極的に利用できる垂直磁気双極子アンテナを採用し、光ファイバ８を使って同期検波を行って高周波分解能を高め、高精度に空洞位置を決定するものである。
【００６６】
続いて、実施の形態１の地中探査方法において、空洞中心を推定して表示する一連の判断手順について説明する。図５は本発明の実施の形態１の地中探査方法で空洞位置を推定する判断手順のフローチャート、図６は本発明の実施の形態１の地中空洞探査装置と地中探査方法で空洞位置を推定して表示する動作のフローチャートである。図５において、図２（ａ）に示すように測線と送受信機部を配置して各測点の測定を行う。フーリエ変換したデータを入力して、１ｋＨｚ〜１００ｋＨｚの複数の周波数で、測点と比抵抗を計算するとともに測点と位相を計算し、これによって所定幅ごとに比抵抗値を分けて等高線を描いた等比抵抗図、所定幅ごとに位相を分けて等高線を描いた等位相図、さらに、１次元解析断面の演算を行い、比抵抗断面図を作成する（ｓｔｅｐ１）。
【００６７】
次いで、等位相図に現れる周縁効果を示す一対のピークを抽出できるか否かをチェックする（ｓｔｅｐ２）。一対のピークを抽出できない場合は、等位相図の低位相領域がほぼ対称か否かをチェックする（ｓｔｅｐ３）。ｓｔｅｐ３において低位相領域が対称の場合、おおむね対称軸が空洞の中心位置となり（ｓｔｅｐ４）、対称でない場合、ピークの位置が周縁効果の一方を示している（ｓｔｅｐ５）。ｓｔｅｐ４，５の後、ｓｔｅｐ２で一対のピークを抽出した場合とともに、一対のピークと対応した測点間の中心位置を算出すればよい（ｓｔｅｐ７）。ｓｔｅｐ１〜７において、等位相図だけでもデータが揃えば空洞中心を推定できることが分かる。
【００６８】
実施の形態１においては、さらに比抵抗図によってより正確な空洞位置を算出するため、位相図における一対のピークを示す測点間の領域で比抵抗のピークを抽出する（ｓｔｅｐ８）。ｓｔｅｐ８において、ピークが存在して、その数が３個以下であるか否かをチェックし（ｓｔｅｐ９）、ピークがないか４個以上のときはｓｔｅｐ２１へ進み、ピークが２個かどうかをチェックする（ｓｔｅｐ１０）。ピーク間の内側の地電流電極の電極間隔が空洞幅より大きいかどうかを判断し（ｓｔｅｐ１１）、大きい場合は２個のピークの中間の位置を空洞中心とし（ｓｔｅｐ１２）、電極間隔が空洞幅より小さい場合、２個のうち右側のピークが位相図における一対のピークの中央付近に存在するかどうかをチェックする（ｓｔｅｐ１３）。
【００６９】
すなわち、電極間隔と空洞幅の関係で、２個のピークが出現した場合は、周縁効果を示す一対のピークだけが現れた場合と、比抵抗のピークが１本しか現れない場合の変形で、空洞直上の比抵抗のピークと周縁効果の一方が現れた場合の２つがあり、ｓｔｅｐ１２で前者の場合を判定し、ｓｔｅｐ１３で後者の場合を判定する。ｓｔｅｐ１３において、右側のピークの位置が等位相図の一対のピークの中央付近に存在する場合、右側のピークの位置を空洞中心とし（ｓｔｅｐ１４）、右側のピークの位置が等位相図の一対のピークの中央付近に存在しない場合、左側のピークの位置を空洞中心とする（ｓｔｅｐ１５）。
【００７０】
ｓｔｅｐ１０において、２個のピークでない場合、３個のピークが存在するか否かをチェックし（ｓｔｅｐ１６）、ピークが１個の場合には、ｓｔｅｐ７の中央位置と一致しているかを確認する（ｓｔｅｐ１７）。一致しない場合はｓｔｅｐ２１へ進み、ｓｔｅｐ７の中央位置を空洞中心とするか、エラーとして終了する。ｓｔｅｐ１７において、ｓｔｅｐ７の中央位置と一致している場合は、このピーク位置が中心効果を示すものとして、これを空洞中心とする（ｓｔｅｐ１９）。また、ｓｔｅｐ１６において、ピークが３個の場合は、中央のピークの位置を空洞中心とする（ｓｔｅｐ１８）。
【００７１】
ｓｔｅｐ１２，１４，１５，１８，１９で空洞中心が定まると、この空洞中心が同一位置または大きくずれても隣接する測点にない場合（所定の範囲内の場合）は、１次元解析を行って（ｓｔｅｐ２２）、空洞の想定される深度と一致するか否かをチェックする（ｓｔｅｐ２３）。空洞の想定深度と一致しているときは空洞との推定が確定でき（ｓｔｅｐ２４）、終了する。ｓｔｅｐ２３において、深度が一致しない場合は別要因による反応があったものとして（ｓｔｅｐ２５）終了する。
【００７２】
続いて、図６に基づいて実施の形態１の地中空洞探査装置で空洞位置を推定して表示する動作について説明する。地電流電極間の距離を大きい間隔から徐々に小さくするなどし、測定条件を変えては、以下のｓｔｅｐ３１〜３３を繰返し、等位相図で適切なピークが得られたかをチェックし、等比抵抗図で３個のピークを求める。
【００７３】
ｓｔｅｐ３１において、フーリエ変換したデータを入力されると、演算手段２２の比抵抗算出手段２４と位相算出手段２５が、１ｋＨｚ〜１００ｋＨｚの複数の周波数で、測点と比抵抗を計算するとともに測点と位相を計算し、これによって所定幅ごとに比抵抗値を分けて等高線を描いた等比抵抗図、所定幅ごとに位相を分けて等高線を描いた等位相図を作成する。さらに、解析手段２６によって１次元解析断面の演算を行い、比抵抗断面図を作成する（ｓｔｅｐ３１）。
【００７４】
次いで、空洞位置算出手段２７が等位相図のデータに基づき、空洞様反応領域を示す周縁効果の一対のピークが得られたかチェックを行う（ｓｔｅｐ３２）。ここで、このピークの位置の計算を容易にするため、例えば各測点ごとにすべての周波数に関する位相の総和を計算し、「測点−位相総和グラフ」を作成し、このグラフにおける低位相のピークを抽出すればよい。
【００７５】
測定条件を変えたとき、地電流電極間隔が空洞幅に対して適切な場合（空洞幅が電極間隔の1／２〜３／４の場合）、空洞におけるピークは、図２（ｃ）における左右の図のように２パターンで示され、いずれもピークが出現する。この一対の周縁効果の現れた測点の内側の地電流電極の間が空洞の存在する範囲となり、その中心が空洞の中心となる。このうち、図２（ｃ）の左図の場合空洞中心に整った形で中心反応も現れるので、測定条件を変えてこの状態が得られるまで繰り返せばよい。このときより正確に空洞位置が推定できる。なお、実施の形態１の地中空洞探査装置では等位相図で一対のパターンを得て第１のチェックを行い、ニアフィールド効果で増幅された等比抵抗図のピークを使って第２の確認を行うから、等位相図でのピークは一対のピークを検出することで足りる。
【００７６】
そして、本地中空洞探査装置においては、この「測点−位相総和グラフ」による総和とピーク値の手法が自動化に馴染むため、これを採用しているが、これに限られるものではない。なお、ピーク値は土被比と地盤の比抵抗およびニアフィールド効果による増幅率の関数で決まるので、一律に閾値を設定するのは適当とはいえない。
【００７７】
さらに、位相の総和の測点に関する分布を対象として、空洞が存在する範囲となる分布形状の左右対称性をチェックしてもよい。これはピーク間の距離で判定できる。なお計算に際して、周波数として周縁効果を十分明確に示す高周波側の１０ｋＨｚ〜１００ｋＨｚを選び、この周波数だけを対象としてピークの位置を抽出するのも好適である。
【００７８】
ｓｔｅｐ３２において周縁効果を示す一対のピークが得られた場合、等比抵抗図のデータに基づき、３個の比抵抗のピークが適切な配置で存在するか確認する（ｓｔｅｐ３３）。計算を容易にするため、上述の位相と同様に、各測点ごとにすべての周波数における比抵抗の総和を計算し、ピーク値の間隔が図２（ｃ）の２パターンの何れかに対応すればよい。なお計算に際して、ニアフィールド効果により比抵抗コントラストが増幅された低周波領域となり、周縁効果を明確に示す１ｋＨｚ〜１０ｋＨｚを選び、この周波数を対象としてピークの位置を抽出するのも好適である。なお、比抵抗におけるピーク値も土被比と地盤の比抵抗およびニアフィールド効果による増幅率の関数で決まるので、一律に閾値を設定するのは適当とはいえない。
【００７９】
ｓｔｅｐ３３で抽出したピークが３個の場合、中央のピークの位置を空洞中心とし（ｓｔｅｐ３４）、終了する。ｓｔｅｐ３３において、３個以外の場合、電極間距離等の測定条件を変更し、測定条件の変更がすべて終わったか確認する（ｓｔｅｐ３５）。例えば、地電流電極間の距離を大きい間隔から徐々に小さくして測定している場合には、予定された測定条件の変更が終了したかを確認し、終了している場合は、測定不能である旨の表示を行ない（ｓｔｅｐ３６）、終了していない場合はｓｔｅｐ３１に戻って、測定を繰り返す。
【００８０】
以上説明したように本発明の実施の形態１の地中空洞探査装置と地中探査方法は、空洞の水平位置を推定するために、従来の電気探査のように数値解析を行って探査するのではなく、位相と比抵抗を相対的に評価することで、簡単な構成で明確に空洞の位置を推定することができる。そして、電気比抵抗が均一な地盤や、電線や電話線等から生じる電磁波ノイズが小さいところ、ノイズとなる地下埋設物が無い測定地においては、とくに顕著な結果を得ることができる。また、空洞位置の捕捉精度、水平分解能は、測点間隔を密にすることで容易に向上させることができる。さらに、比抵抗および位相のピーク値は土被比と地盤の比抵抗およびニアフィールド効果による増幅率の関数で決まるので、深く土被り比が大きい空洞におけるピーク値の低下に対して、送受信間距離を短くすることにより、ニアフィールド効果による増幅率を上げることができる。なお、実施の形態１の地中空洞探査装置と地中探査方法は、空洞そのものを確認しているものではないので、最終確認のためのボーリング等を行うことが望まれる。
【００８１】
（実施の形態２）
本発明の実施の形態２の地中空洞探査装置と地中探査方法は、調査地内において空洞位置の推定が困難な場合、闇雲に測線を設けるのではなく、簡便に測線を空洞上に空洞に対し、ほぼ直交方向に設置するまでの形態である。従って、実施の形態１の地中空洞探査装置と地中探査方法のいわば前提となるものである。
【００８２】
ＣＳＡＭＴ探査には、シャドウ効果という現象が知られている。これは、送受信機間に比抵抗異常がある場合に、その影が受信点の測定に影響を与える現象であり、発信位置を移動させるというシャドウ効果の排除方法が既存の研究で提案されている。本発明の実施の形態２は、シャドウ効果を排除するのでは無くこれを利用し、測線直下の空洞の有無、測線に対する未知の空洞の方向、測線と未知の空洞との距離、未知の空洞の向きを推定することで、測線を空洞上に空洞に対してほぼ直交方向に設置することができるものである。図７は本発明の実施の形態２における地中空洞探査方法を行うための地中空洞探査装置のシャドウ効果の説明図、図８は本発明の実施の形態２における地中空洞探査方法を行うための地中空洞探査装置の測線と送信部のＸ−Ｙ平面配置図である。
【００８３】
図７に示すように空洞調査のときのシャドウ効果は、探査測線の地下に空洞がなくても、空洞の延長上にシャドウ効果により空洞があるときと似た空洞反応が生じる。このシャドウ効果を除くために実施の形態１の地中空洞探査装置では、送信機部を測線の両翼に２箇所設けて両側から送信される電磁波を測定する。
【００８４】
実施の形態２の地中空洞探査装置は、測線と２つの送信機部を図８に示すように配置する。送信機部が２箇所になった点を除き、基本的に実施の形態１と同一であるから、図１〜図５を参照し、同一符号の説明は実施の形態１に譲ってここでは省略する。
【００８５】
実施の形態２においては、地下空洞に対して図８に示すように（併せて図１、図２も参照）、測線と送信機部Ａ，Ｂ、地電流電極、高周波用磁気センサ，低周波用磁気センサのインダクションコイルを配置する。想定される地下の空洞に直交すると思われる方向（Ｙ方向）に測線を設置し、測線の両翼（Ｘ方向）に５０ｍ〜１００ｍ測線から離した地点に送信機部Ａ，Ｂを置く。また、受信機部の地電流電極を測線方向に，インダクションコイルを測線の直交方向に設置する。地電流電極の間隔は、想定される空洞幅の１／２〜３／４程度とする。測定間隔は、地電流電極の間隔と同じかより短くする。1測点に付き送信機Ａ，Ｂからの信号を個別に測定する。
【００８６】
実施の形態２の地中空洞探査装置は、探査測線上において空洞の延長部に高比抵抗異常が検出されるため、シャドウ効果が出た測点で、図７の周縁効果の幅方向に測定方向を回転させるものである。この回転をさせることで、最大の比抵抗値と最小の位相値が得られた方向、あるいは最小の比抵抗値と最大の位相値が得られた方向の直交方向に空洞が存在することが分かる。この未知の空洞までの距離によってシャドウ効果の広がりが異なるため、この広がりを測定することでこの空洞までの距離を推定することができる。
【００８７】
ところで、測線の直下に空洞が存在することが未確定の場合、まず、送信機部Ａ，Ｂごとに図４に示すような比抵抗と等位相図を作成する。送信機部Ａ，Ｂの比抵抗と等位相図の両者で、空洞様反応が出た場合は測線の直下に空洞があると判断される。片方の発信機部で空洞様反応が生じた場合、シャドウ効果により空洞様反応が出た送信機と測線の間に空洞が存在する。
【００８８】
次に、強いシャドウ効果が現れた測点で、測定方向を少しずつ回転させ、最大の比抵抗値と最小の位相値が得られた方向と最小の比抵抗値と最大の位相値が得られた方向を記録する。理論的に前者の方向および後者の直交方向が未知の空洞がある方向を示す。このシャドウ効果の角度測定が周縁効果に対応するピークの２箇所で行えた場合、未知空洞の横幅を推定することにより、空洞までの距離が作図より求められる。
【００８９】
実施の形態２の地中空洞探査装置を使うことによって、送受信機部間の空洞の存在を推定することができ、今後の空洞調査の経費削減に大きく寄与することができる。測線をグリッド状に張り巡らせるのではなく、地中空洞の位置を推定することができる。
【００９０】
【実施例】
（実施例１）
鹿児島県鹿屋市周辺のシラス台地を探査地点として選んだ。この探査地においては、防空壕が多数存在する。図９は探査地における測線と測定番号の配置図、図１０（ａ）はＮｏ.１測線の等比抵抗図、図１０（ｂ）はＮｏ.１測線の等位相図、図１０（ｃ）は実測空洞の位置を表示したＮｏ.１測線の比抵抗２次元解析断面図である。
【００９１】
実施例１のＮｏ.１測線において、空洞Ａは図１０（ａ）（ｂ）から測点３と推定され、空洞Ｂは、図１０（ｂ）から選んだ測点８〜測点１１の範囲のうち、図１０（ａ）の測点９と推定される。空洞Ｃは、図１０（ｂ）において周縁反応のピーク位置は測点１３と測点１７であり、その中心を取って測点１５が空洞の中心と推定されるが、図１０（ａ）を見ると、その対称性から測点１５が空洞に関係していると見ることができる。空洞の実測位置は測点１４〜１５でこの推定と一致する。空洞Ｄは、図１０（ｂ）からも分かるように一対のピーク値の一方しかないため、断定は難しい。
【００９２】
しかし、実施例１の探査結果はきわめて正確に空洞位置を測定しているのが分かる。
【００９３】
（実施例２）
図１１（ａ）はＮｏ.２測線の等比抵抗図、図１１（ｂ）はＮｏ.２測線の等位相図、図１１（ｃ）は実測空洞の位置を表示したＮｏ.２測線の比抵抗２次元解析断面図である。
【００９４】
実施例２のＮｏ.２測線において、空洞Ａは、図１１（ｂ）から選んだ測点２〜測点５の範囲のうち中央の測点３と推定され、空洞Ｂは、図１１（ｂ）から選んだ測点９〜測点１０の範囲のうち、図１１（ａ）の測点１０と推定される。空洞Ｃは、図１１（ｂ）から選んだ測点１４〜測点１６の範囲のうち中央の測点１５と推定される。空洞の本坑と測線が低角度で交差しているためか、図１１（ａ）は測点１８まで高比抵抗が広がっている。空洞Ｄは、図１１（ｂ）から測点２０が周縁効果の一方で、空洞は終点側にあると見ることができる。このとき空洞が他のＡ，Ｂ，Ｃ孔とほぼ同じ大きさだと仮定すれば、空洞位置は測点２２と推定される。しかし、この測線における測点２０〜終点（測点２２）の空洞様反応は、空洞の延長部で生じたシャドウ効果である。すなわち、実施例２のように発信が１箇所だけで測定するとシャドウ効果が現れた場合に、殆ど測線直下にある空洞と判別できない。
【００９５】
このように、実施例２の探査結果は、空洞Ａ，Ｂで正確に推定できているが、空洞Ｃでは斜交部分で等比抵抗図（図１１（ａ））に異常が表れ、空洞Ｄは直上に測線が設けられていないにも係わらず、シャドウ効果により空洞様反応が現れていることが分かる。
【００９６】
【発明の効果】
本発明の地中空洞探査方法によれば、高周波の電磁波を送信する高周波ＣＳＭＴ法による探査を基本とし、これまでの探査手法では測定精度が悪かった地下１０ｍ以深（１０ｍ〜１００ｍ）を含む地下１ｍ〜１００ｍ程度の浅部から「大深度地下の公共的使用に関する特別措置法」の定める大深度地下の空洞位置を検出でき、垂直磁気双極子アンテナ等のアンテナを使用してニアフィールド効果により空洞の所在を示す空洞様反応を増幅することによって比抵抗や位相における空洞様反応を強調でき、送信機部と受信機部とで同期検波を行うためＳ／Ｎ比が上がり、精度よく安価に且つ自動的に探査することができる。また、シャドウ効果による比抵抗異常を判別することができ、さらには測線直下の空洞の有無，測線直下に無い空洞がある方向，測線直下に無い空洞の向き，およびその空洞までの距離を推定することができる。
【図面の簡単な説明】
【図１】 本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置の構成図
【図２】 （ａ）本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置の測線と送信部のＸ−Ｙ平面配置図、
（ｂ）（ａ）の地中空洞探査装置の測線とセンサと地下空洞のＸ−Ｙ−Ｚ配置
（ｃ）（ａ）の地中空洞探査装置の周縁効果の出現位置の説明図
【図３】 （ａ）比抵抗モデルの説明図
（ｂ）（ａ）の比抵抗の説明図
（ｃ）（ａ）の位相の説明図
【図４】 本発明の実施の形態１における地中空洞探査方法を行うための地中空洞探査装置が空洞位置に対して測定した等比抵抗と等位相の測定図
【図５】 本発明の実施の形態１の地中探査方法で空洞位置を推定する判断手順のフローチャート
【図６】 本発明の実施の形態１の地中空洞探査装置と地中探査方法で空洞位置を推定して表示する動作のフローチャート
【図７】 地中空洞探査装置のシャドウ効果の説明図
【図８】 本発明の実施の形態２における地中空洞探査装置の測線と送信部のＸ−Ｙ平面配置図
【図９】 探査地における測線と測定番号の配置図
【図１０】 （ａ）Ｎｏ.１測線の等比抵抗図
（ｂ）Ｎｏ.１測線の等位相図
（ｃ）実測空洞の位置を表示したＮｏ.１測線の比抵抗２次元解析断面図
【図１１】 （ａ）Ｎｏ.１測線の等比抵抗図
（ｂ）Ｎｏ.２測線の等位相図
（ｃ）実測空洞の位置を表示したＮｏ.２測線の比抵抗２次元解析断面図
【図１２】 図４と同測線における水平磁気双極子アンテナによる等比抵抗と等位相の測定図
【符号の説明】
１ 発振部
２ 周波数逓倍部
２ａ スイッチング部
３ ドライバ部
４ａ 高周波アンテナ
４ｂ 低周波アンテナ
５ 送信側制御部
６ 計時手段
７ 記憶部
８ 光ファイバ
８ａ，８ｂ 信号変換部
９ 同期部
１１ａ 高周波用磁気センサ
１１ｂ 低周波用磁気センサ
１２ａ，１２ｂ 地電流電極
１３ａ，１３ｂ プリアンプ
１４ａ，１４ｂ バッファアンプ
１５，１６ バンドパスフィルタ
１７，１８ 増幅部
１９ フーリエ変換部
１９ａ 高周波変換部
Claims (3)
- 地中探査のため高周波信号を電磁波として送信するとともに、一対の地電流電極で測線上の測点の電位差を検出するとともに該測線と直交する方向に配設されたインダクションコイルにより該測点の磁場を測定し、該測点における比抵抗と位相を演算する高周波ＣＳＭＴ法による地中空洞探査方法であって、
空洞の所在を示す空洞様反応をニアフィールド効果によって増幅するための垂直磁気双極子アンテナを前記空洞から所定距離おいた位置に配置して高周波信号を送信し、探査領域における周波数と測点に関する等比抵抗データを取得して、該等比抵抗データの中からニアフィールド効果により増幅された空洞様反応領域を抽出し、該空洞様反応領域の測線方向幅の中央を空洞中心として推定することを特徴とする地中空洞探査方法。 - 請求項１の地中空洞探査方法において、探査領域における周波数と測点に関する等比抵抗データを取得して、該空洞様反応領域の測線方向幅の中央を空洞中心として推定するのに代えて、探査領域における周波数と測点に関する等位相データの中から空洞様反応のピークを含む低位相の空洞様反応領域を抽出し、この空洞様反応領域に対応する等比抵抗データの中から比抵抗のピークを抽出して空洞中心を推定することを特徴とする地中空洞探査方法。
- 前記測線の両側から電磁波を送信し、該測線の直下にある空洞と測線の直下に無い空洞とを、シャドウ効果による空洞様反応に近似した比抵抗または位相異常の有無から判別することを特徴とする請求項１又は２に記載の地中空洞探査方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003147120A JP4229371B2 (ja) | 2003-05-26 | 2003-05-26 | 地中空洞探査方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003147120A JP4229371B2 (ja) | 2003-05-26 | 2003-05-26 | 地中空洞探査方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004347541A JP2004347541A (ja) | 2004-12-09 |
| JP4229371B2 true JP4229371B2 (ja) | 2009-02-25 |
Family
ID=33533741
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003147120A Expired - Fee Related JP4229371B2 (ja) | 2003-05-26 | 2003-05-26 | 地中空洞探査方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4229371B2 (ja) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO323889B1 (no) * | 2005-11-03 | 2007-07-16 | Advanced Hydrocarbon Mapping A | Framgangsmate for kartlegging av hydrokarbonreservoarer samt apparat for anvendelse ved gjennomforing av framgangsmaten |
| KR100945362B1 (ko) | 2007-10-09 | 2010-03-05 | 한국지질자원연구원 | 전자탐사장치, 전자탐사방법 및 프로그램을 기록한 컴퓨터로 읽을 수 있는 매체 |
| JP4512790B2 (ja) * | 2007-11-22 | 2010-07-28 | 独立行政法人 日本原子力研究開発機構 | 地磁気地電流法における観測データの処理方法及び装置 |
| RU2479858C1 (ru) * | 2011-12-05 | 2013-04-20 | Федеральное Государственное Автономное Образовательное Учреждение Высшего Профессионального Образования "Сибирский Федеральный Университет" (Сфу) | Электроразведочное устройство |
| KR101348787B1 (ko) | 2013-10-16 | 2014-01-07 | 한국광물자원공사 | 2차원 수치모델링 및 역산을 이용하여 사전정보를 반영하는 3차원 자력역산 방법 |
| US20180209787A1 (en) * | 2015-07-14 | 2018-07-26 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Method for detecting change in underground environment by using magnetic induction, detection sensor and detection system |
| JP7030496B2 (ja) * | 2017-12-04 | 2022-03-07 | 三菱電機株式会社 | 地中空洞検出装置および地中空洞検出方法 |
| WO2020194371A1 (ja) * | 2019-03-22 | 2020-10-01 | 日本電気株式会社 | 地下空洞検査システム、及び地下空洞検査方法 |
| CN112462432B (zh) * | 2020-11-16 | 2024-02-02 | 新疆维吾尔自治区煤田灭火工程局 | 煤田火区空洞高精度探测方法及装置 |
| JP7508420B2 (ja) | 2021-09-01 | 2024-07-01 | 三郎 上森 | 空洞探査器 |
-
2003
- 2003-05-26 JP JP2003147120A patent/JP4229371B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2004347541A (ja) | 2004-12-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7030617B2 (en) | System, apparatus, and method for conducting electromagnetic induction surveys | |
| US10753722B1 (en) | Systems and methods for locating buried or hidden objects using sheet current flow models | |
| EP3341766B1 (en) | Locating devices, systems, and methods using frequency suites for utility detection | |
| CN103869371B (zh) | 人工场源频率域全梯度电磁测量方法 | |
| US20110066379A1 (en) | survey system for locating geophysical anomalies | |
| KR101267017B1 (ko) | 지피알 탐사장비의 신호처리에 의한 지하시설물 탐측시스템 | |
| JP4229371B2 (ja) | 地中空洞探査方法 | |
| CN102062876A (zh) | 一种全区电偶源频率域电测深方法 | |
| US6525539B2 (en) | Apparatus and method for locating subsurface objects in conductive soils by measurements of magnetic fields by induced currents with an array of multiple receivers | |
| KR100399984B1 (ko) | 전자 지하탐사방법 및 지하탐사시스템 | |
| CN104749643A (zh) | 一种用于电磁干扰地区的csamt测量方法和装置 | |
| WO2021098667A1 (zh) | 一种全覆盖抗干扰人工源电磁探测方法 | |
| Sadatcharam et al. | Depth sensitivity of apparent magnetic susceptibility measurements using multi-coil and multi-frequency electromagnetic induction | |
| Sogade et al. | Electromagnetic cave-to-surface mapping system | |
| CN114280673B (zh) | 基于静态效应的物探仪及数据采集方法 | |
| CN116430464A (zh) | 提高时频电磁法频率域资料处理质量的方法 | |
| US10156134B2 (en) | Method and apparatus for passive detection of near-surface human-scale underground anomalies using earth field measurements | |
| CN114879267B (zh) | 一种基于mt静态效应的天然磁场选频勘探方法 | |
| CN114897920B (zh) | 一种基于索贝尔算法的道路空洞边缘分割方法 | |
| WO2014113516A1 (en) | Phase-based electromagnetic surveys for geological formations | |
| KR200311086Y1 (ko) | 전자 지하탐사장치 | |
| Zhang et al. | Efficient and High-Resolution Detection for Urban Underground Space Using a Grounded Wire Source Frequency-Domain Electromagnetic Gradient Method | |
| JP3041415B2 (ja) | 楕円体パラメータを測定する電磁探査法 | |
| CN114153006A (zh) | 一种全向视电导率示踪勘探方法 | |
| Waff | Electrical Resistivity Imaging of Groundwater Systems using Natural and Controlled Source Magnetotellurics |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20060515 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20060517 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060518 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20060515 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20060517 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070810 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080902 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20081024 |
|
| 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: 20081117 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20081128 |
|
| R150 | Certificate of patent (=grant) or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111212 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121212 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121212 Year of fee payment: 4 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: R3D02 |
|
| LAPS | Cancellation because of no payment of annual fees |