JP4311745B2 - 地中探査装置、システム及び方法 - Google Patents
地中探査装置、システム及び方法 Download PDFInfo
- Publication number
- JP4311745B2 JP4311745B2 JP2005512913A JP2005512913A JP4311745B2 JP 4311745 B2 JP4311745 B2 JP 4311745B2 JP 2005512913 A JP2005512913 A JP 2005512913A JP 2005512913 A JP2005512913 A JP 2005512913A JP 4311745 B2 JP4311745 B2 JP 4311745B2
- Authority
- JP
- Japan
- Prior art keywords
- voltage
- impedance
- electrodes
- current
- unit
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/02—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current
- G01V3/06—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current using ac
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/043—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a granular material
Description
この発明は地中の漏水状態、地下水脈や埋蔵資源あるいは土壌汚染の位置・形状・深度・収蔵量、土壌改善の進行状況を低侵襲性の交流電気試験手段により探査あるいは観測する計測装置に関する。さらに詳しくは、この発明は、大地に電流を通流することにより出現する電圧の多点計測を行い、２以上の周波数における計測結果を得、有限要素法または境界要素法またはインピーダンス網等を用いた大地モデルを使用して、大地モデルの局所領域の導電率と周波数別の誘電率、位置情報等を未知数量として実測値と計算値の誤差が小さくなるように未知数量を変更することにより地中物質の同定や地質の変化の検出を行い、大地モデルの入力情報と推定処理結果を２次元あるいは３次元表示することにより実現する地中探査装置に関する。
地中を探査する従来技術としては、地中を部分的に掘り返す侵襲的なボーリング方法以外に、非侵襲的な地中探査方法としては１）人工的に振動を加えて局所地震を起こす音響探査方法、２）軟Ｘ線を照射して高導電率の物質からの後方散乱光を撮影するＸ線バック・スキャッタ法、３）マイクロ波照射による地中レーダー探査法、４）上空から近赤外線領域を撮影する赤外線撮影法、５）上空からミリ波散乱をモニターする探査方法等が知られている。
人工的に振動を加えて局所地震を起こす音響探査方法は、同一地質の領域内では物質密度がほぼ同じで音響インピーダンスが等しいので振動波の反射は殆ど発生せず、地質が異なる境界面においては物質密度が異なり音響インピーダンスが不連続であるために振動波が反射する現象を利用しているので、地質不連続面からの反射波を捉えることにより各地層の地質に対応する音速が既知であれば地層の厚みや深度を３次元的に精確に計測して可視化表示することができる。しかし、未知の物質が存在する地層は密度がほぼ同一である場合には振動波の反射が発生しないために発見することはできない欠点がある。即ち、音響探査方法では地層の深さや厚みを知ることはできるが、その地層がどのような物質で構成されているかを適切なボーリング調査を行わずに計測することはできない。例えば、地下水脈の深さや埋蔵資源の存在検知、あるいは土壌が化学物質で汚染されているか否かを計測することはできないのである。
近年実現された地中探査方法として軟Ｘ線を照射して高導電率の物質からの後方散乱光を撮影するＸ線バック・スキャッタ法がある。これは軟Ｘ線が電磁波と粒子の性質を合わせ持ったものであり、電磁波の表皮効果のために導電性や透磁率が低い値を有する物質あるいは誘電率の高い物質あるいは軽い原子核が含まれる物質に対しては透過性に優れるが、導電性や透磁率が高い値を有する物質あるいは誘電率の低い物質あるいは重い原子核が含まれる物質に対しては透過性が低く、散乱あるいは反射する性質を利用した探査方法である。深さ数１０ｃｍ程度の浅い地中の埋設物の有無を調べるのには適した方法であるが、含水量の多い土壌ではほぼ表面付近しか探査できない欠点や地質判定が行えない欠点がある。
マイクロ波照射による地中レーダー探査法も電磁波の表皮効果のために含水量が多い土壌には電磁波が入り込めず減衰してしまうので比較的浅深度の探査にしか適用することができない。また、地質を判定することができない。
上空から近赤外線領域を撮影する赤外線撮影法では大地から輻射される近赤外光を捉えるので、地中に周囲に比して高温あるいは低温の領域が存在する等の温度に関連した情報収集には適しているが、物質の種類や深さを同定することはできない。
上空からミリ波散乱をモニターする探査方法は、太陽活動に伴って地球を取り巻くバンアレン帯から入射するミリ波が地上で散乱しており、鉄分を多く含んだ土壌はミリ波を吸収する性質を利用した探査方法である。しかし、ほぼ地表付近の地質探査にしか適用することができない欠点がある。
この他に、非侵襲あるいは低侵襲性の地中探査方法として、大地に電流を通流して、電位計測あるいは導電率計測を行う方法が知られている。
特許庁の特許電子図書館の公報テキスト検索により検索式として［請求の範囲＝ｏｒ（地中 地質 探査）］ＡＮＤ［請求の範囲＝ａｎｄ（電流 電極）］を用いた結果では以下の文献が検索された。
埋設管の埋設位置あるいは腐食位置あるいは塗膜損傷部位の検出技術としては以下の文献が存在した。
１）特開２００２−０７１８２２ Ｇ０１Ｖ １／００ 大山 順彦 埋設管の位置検出方法
２）特開２００３−００４６８６ Ｇ０１Ｎ ２７／２０ 新日本製鐵（株） 積算手段を用いた埋設金属管類の防食被覆損傷位置検出方法
３）特開２００３−００４６８７ Ｇ０１Ｎ ２７／２０ 新日本製鐵（株） ２種類の周波数信号を用いた埋設金属管類の防食被覆損傷位置検出方法
４）特開２０００−２５８２８０ Ｇ０１Ｍ ３／１８ 富士電機（株）；高嶋 武 水道管の漏水検知方法および装置
５）特開２００１−０１３１１１ Ｇ０１Ｎ ２７／８２ 日本鋼管（株） 埋設塗覆装鋼管の塗膜損傷位置検出方法及び装置
６）特開２０００−０１９１５９ Ｇ０１Ｎ ２７／８２ 日本鋼管（株） 埋設鋼管の塗膜損傷位置検出方法及び装置
７）特開２０００−０１９１５８ Ｇ０１Ｎ ２７／８２ 日本鋼管（株） 埋設鋼管の塗膜損傷位置検出方法及び装置
８）特開２０００−０８８６９１ Ｇ０１Ｍ ３／１６ 西松建設（株） 漏水測定装置及び漏水測定方法
９）特開平０６−０９４５６８ Ｇ０１Ｍ ３／４０ ７３２４−２Ｇ 丸山製作所：（有） 漏水検出方法および漏水検出装置
文献１）は２電極を用いて直流による導電率計測を行い、測定位置を変えながら間接的に計測する方法で、埋設管を挟んだ位置で測定する場合には通流分布が変化することを利用したものであるが、含水量、埋設管径、埋設方向、埋設管種、深さ、他の埋設物の影響を受け易い欠点がある。文献２）、３）、４）、５）、６）、及び７）は電流印加電極の一方を埋設金属管にする方法である。文献８）は予め配置した多数の電極は配列上に漏水が発生した場合に電極間インピーダンスが変化することを利用した直接計測法である。文献９）は導水経路を挟んだ大地の電気抵抗が閉栓時と開栓時で変化することを利用した間接的計測方法である。
交流電気試験法を用いた地質調査の先行技術としては以下の文献が存在した。
１０）特開平１０−２９３１８１ Ｇ０１Ｖ ３／０６ 応用地質（株）キャパシタ法による多チャンネル電気探査方法
１２）特開２００１−０７４８５０ Ｇ０１Ｖ ３／０８ 応用地質（株） 無分極電極を用いる電気探査方法
１１）特開平１０−２９３１８１ Ｇ０１Ｖ ３／０６ 応用地質（株） キャパシタ法による多チャンネル電気探査方法
１３）特開平０９−１２７２５３ Ｇ０１Ｖ ３／０８ 応用地質（株） キャパシタ法によるマルチチャンネル電気探査システム
１４）特開２００２−１５６４６０ Ｇ０１Ｖ ３／０２ 産学連携機構九州：（株） 電気探査方法とそれを用いた電気探査装置及び地雷探知装置
１５）特開２０００−０２８７４３ Ｇ０１Ｖ ３／１２ 三井造船（株） 埋設物探査方法および空洞検知用非接触高周波電流センサ
１６）特開平０７−０１２７６６ Ｇ０１Ｎ ２７／００ Ｄ ９１１５−２Ｊ 鹿島建設（株） 比抵抗利用の地下流れ特性測定方法及び装置
文献１０）、１１）は一般的な交流電気試験法により周波数と距離減衰の影響を補正して直流比抵抗を精度よく測定する方法について触れた文献である。文献１２）は一般的な交流電気試験法により電極の充電率と直流比抵抗を測定する方法である。文献１３）は交流電気試験により比抵抗を計測するシステムの実施事例である。文献１４）は無限遠方電流電極を設けることにより電流経路に制約を与え、測定用電流電極の位置を変更しながら測定用電流電極周囲の電位を測定することで測定用電流電極の直下の埋設物の有無を検出する装置である。文献１５）は交流電気試験により導電率を計測する方法である。文献１６）はトレース流体注入と交流電気試験を併用する方法である。
また、交流電気試験法ではあるが、電位を計測するのではなく印加周波数と同じ共振周波数の検出器を使う信号追尾法（シグナル・トレーシング法）の文献としては次のものが存在した。
１７）特開平０６−１３０１５６ Ｇ０１Ｖ ３／０８ Ｂ ７２５６−２Ｇ 日本電信電話（株） 金属埋設物探知方法及び装置
電流の通流経路を工夫して比抵抗を数値演算で求めてトモグラフ表示させる装置の文献として次のものが存在した。
特開平０５−１０００４４ Ｇ０１Ｖ ３／０２ Ｃ ７２５６−２Ｇ ダイヤコンサルタント：（株） 自動電気探査方法
上記の先行技術においては、地中に埋設された水道管からの漏水状態に関しては工夫次第で検出することはできるのであるが、例えばガス管の埋設されている場所に近接した場所に埋設された水道管からガス管方向に対して発生した漏水や急激な漏水に伴って発生するサンドブラスト現象（砂や土砂が水流とともに噴きつける現象）の発生を検出することはできない欠点が存在している。
また、上記の先行技術においては、比抵抗（導電率）の分布状態は検出することが可能であるが、導電率は周波数によっても殆ど変化しないために、地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を計測あるいはモニターすることはできない欠点が存在している。
本発明が解決しようとする課題は、低侵襲性の交流電気試験法により地中探査を行う技術を実現することである。さらに詳細に言えば、地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を計測あるいはモニターすることが可能な技術的手段を実現することである。
さらに別の解決しようとする課題は、送電線や電車架線を流れる電流が大地に帰還あるいは漏洩することによって生じる交流雑音障害からの影響を受けることなく、交流電気試験法を行う方法である。
そしてさらに別の解決しようとする課題は、ガス管近傍に埋設された水道管からの漏水により発生したサンドブラスト現象を検出することである。
人工的に振動を加えて局所地震を起こす音響探査方法は、同一地質の領域内では物質密度がほぼ同じで音響インピーダンスが等しいので振動波の反射は殆ど発生せず、地質が異なる境界面においては物質密度が異なり音響インピーダンスが不連続であるために振動波が反射する現象を利用しているので、地質不連続面からの反射波を捉えることにより各地層の地質に対応する音速が既知であれば地層の厚みや深度を３次元的に精確に計測して可視化表示することができる。しかし、未知の物質が存在する地層は密度がほぼ同一である場合には振動波の反射が発生しないために発見することはできない欠点がある。即ち、音響探査方法では地層の深さや厚みを知ることはできるが、その地層がどのような物質で構成されているかを適切なボーリング調査を行わずに計測することはできない。例えば、地下水脈の深さや埋蔵資源の存在検知、あるいは土壌が化学物質で汚染されているか否かを計測することはできないのである。
近年実現された地中探査方法として軟Ｘ線を照射して高導電率の物質からの後方散乱光を撮影するＸ線バック・スキャッタ法がある。これは軟Ｘ線が電磁波と粒子の性質を合わせ持ったものであり、電磁波の表皮効果のために導電性や透磁率が低い値を有する物質あるいは誘電率の高い物質あるいは軽い原子核が含まれる物質に対しては透過性に優れるが、導電性や透磁率が高い値を有する物質あるいは誘電率の低い物質あるいは重い原子核が含まれる物質に対しては透過性が低く、散乱あるいは反射する性質を利用した探査方法である。深さ数１０ｃｍ程度の浅い地中の埋設物の有無を調べるのには適した方法であるが、含水量の多い土壌ではほぼ表面付近しか探査できない欠点や地質判定が行えない欠点がある。
マイクロ波照射による地中レーダー探査法も電磁波の表皮効果のために含水量が多い土壌には電磁波が入り込めず減衰してしまうので比較的浅深度の探査にしか適用することができない。また、地質を判定することができない。
上空から近赤外線領域を撮影する赤外線撮影法では大地から輻射される近赤外光を捉えるので、地中に周囲に比して高温あるいは低温の領域が存在する等の温度に関連した情報収集には適しているが、物質の種類や深さを同定することはできない。
上空からミリ波散乱をモニターする探査方法は、太陽活動に伴って地球を取り巻くバンアレン帯から入射するミリ波が地上で散乱しており、鉄分を多く含んだ土壌はミリ波を吸収する性質を利用した探査方法である。しかし、ほぼ地表付近の地質探査にしか適用することができない欠点がある。
この他に、非侵襲あるいは低侵襲性の地中探査方法として、大地に電流を通流して、電位計測あるいは導電率計測を行う方法が知られている。
特許庁の特許電子図書館の公報テキスト検索により検索式として［請求の範囲＝ｏｒ（地中 地質 探査）］ＡＮＤ［請求の範囲＝ａｎｄ（電流 電極）］を用いた結果では以下の文献が検索された。
埋設管の埋設位置あるいは腐食位置あるいは塗膜損傷部位の検出技術としては以下の文献が存在した。
１）特開２００２−０７１８２２ Ｇ０１Ｖ １／００ 大山 順彦 埋設管の位置検出方法
２）特開２００３−００４６８６ Ｇ０１Ｎ ２７／２０ 新日本製鐵（株） 積算手段を用いた埋設金属管類の防食被覆損傷位置検出方法
３）特開２００３−００４６８７ Ｇ０１Ｎ ２７／２０ 新日本製鐵（株） ２種類の周波数信号を用いた埋設金属管類の防食被覆損傷位置検出方法
４）特開２０００−２５８２８０ Ｇ０１Ｍ ３／１８ 富士電機（株）；高嶋 武 水道管の漏水検知方法および装置
５）特開２００１−０１３１１１ Ｇ０１Ｎ ２７／８２ 日本鋼管（株） 埋設塗覆装鋼管の塗膜損傷位置検出方法及び装置
６）特開２０００−０１９１５９ Ｇ０１Ｎ ２７／８２ 日本鋼管（株） 埋設鋼管の塗膜損傷位置検出方法及び装置
７）特開２０００−０１９１５８ Ｇ０１Ｎ ２７／８２ 日本鋼管（株） 埋設鋼管の塗膜損傷位置検出方法及び装置
８）特開２０００−０８８６９１ Ｇ０１Ｍ ３／１６ 西松建設（株） 漏水測定装置及び漏水測定方法
９）特開平０６−０９４５６８ Ｇ０１Ｍ ３／４０ ７３２４−２Ｇ 丸山製作所：（有） 漏水検出方法および漏水検出装置
文献１）は２電極を用いて直流による導電率計測を行い、測定位置を変えながら間接的に計測する方法で、埋設管を挟んだ位置で測定する場合には通流分布が変化することを利用したものであるが、含水量、埋設管径、埋設方向、埋設管種、深さ、他の埋設物の影響を受け易い欠点がある。文献２）、３）、４）、５）、６）、及び７）は電流印加電極の一方を埋設金属管にする方法である。文献８）は予め配置した多数の電極は配列上に漏水が発生した場合に電極間インピーダンスが変化することを利用した直接計測法である。文献９）は導水経路を挟んだ大地の電気抵抗が閉栓時と開栓時で変化することを利用した間接的計測方法である。
交流電気試験法を用いた地質調査の先行技術としては以下の文献が存在した。
１０）特開平１０−２９３１８１ Ｇ０１Ｖ ３／０６ 応用地質（株）キャパシタ法による多チャンネル電気探査方法
１２）特開２００１−０７４８５０ Ｇ０１Ｖ ３／０８ 応用地質（株） 無分極電極を用いる電気探査方法
１１）特開平１０−２９３１８１ Ｇ０１Ｖ ３／０６ 応用地質（株） キャパシタ法による多チャンネル電気探査方法
１３）特開平０９−１２７２５３ Ｇ０１Ｖ ３／０８ 応用地質（株） キャパシタ法によるマルチチャンネル電気探査システム
１４）特開２００２−１５６４６０ Ｇ０１Ｖ ３／０２ 産学連携機構九州：（株） 電気探査方法とそれを用いた電気探査装置及び地雷探知装置
１５）特開２０００−０２８７４３ Ｇ０１Ｖ ３／１２ 三井造船（株） 埋設物探査方法および空洞検知用非接触高周波電流センサ
１６）特開平０７−０１２７６６ Ｇ０１Ｎ ２７／００ Ｄ ９１１５−２Ｊ 鹿島建設（株） 比抵抗利用の地下流れ特性測定方法及び装置
文献１０）、１１）は一般的な交流電気試験法により周波数と距離減衰の影響を補正して直流比抵抗を精度よく測定する方法について触れた文献である。文献１２）は一般的な交流電気試験法により電極の充電率と直流比抵抗を測定する方法である。文献１３）は交流電気試験により比抵抗を計測するシステムの実施事例である。文献１４）は無限遠方電流電極を設けることにより電流経路に制約を与え、測定用電流電極の位置を変更しながら測定用電流電極周囲の電位を測定することで測定用電流電極の直下の埋設物の有無を検出する装置である。文献１５）は交流電気試験により導電率を計測する方法である。文献１６）はトレース流体注入と交流電気試験を併用する方法である。
また、交流電気試験法ではあるが、電位を計測するのではなく印加周波数と同じ共振周波数の検出器を使う信号追尾法（シグナル・トレーシング法）の文献としては次のものが存在した。
１７）特開平０６−１３０１５６ Ｇ０１Ｖ ３／０８ Ｂ ７２５６−２Ｇ 日本電信電話（株） 金属埋設物探知方法及び装置
電流の通流経路を工夫して比抵抗を数値演算で求めてトモグラフ表示させる装置の文献として次のものが存在した。
特開平０５−１０００４４ Ｇ０１Ｖ ３／０２ Ｃ ７２５６−２Ｇ ダイヤコンサルタント：（株） 自動電気探査方法
上記の先行技術においては、地中に埋設された水道管からの漏水状態に関しては工夫次第で検出することはできるのであるが、例えばガス管の埋設されている場所に近接した場所に埋設された水道管からガス管方向に対して発生した漏水や急激な漏水に伴って発生するサンドブラスト現象（砂や土砂が水流とともに噴きつける現象）の発生を検出することはできない欠点が存在している。
また、上記の先行技術においては、比抵抗（導電率）の分布状態は検出することが可能であるが、導電率は周波数によっても殆ど変化しないために、地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を計測あるいはモニターすることはできない欠点が存在している。
本発明が解決しようとする課題は、低侵襲性の交流電気試験法により地中探査を行う技術を実現することである。さらに詳細に言えば、地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を計測あるいはモニターすることが可能な技術的手段を実現することである。
さらに別の解決しようとする課題は、送電線や電車架線を流れる電流が大地に帰還あるいは漏洩することによって生じる交流雑音障害からの影響を受けることなく、交流電気試験法を行う方法である。
そしてさらに別の解決しようとする課題は、ガス管近傍に埋設された水道管からの漏水により発生したサンドブラスト現象を検出することである。
本発明の第一の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電極３と、電極３の両端に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンスの差を演算して発振周波数の違いによるインピーダンス差を計算する演算部６と、演算部６が演算するインピーダンス差に基づいて短時間に電極３の２点間のインピーダンスが少なくとも一方の周波数で大きく変化したことを検出判定する判定部７により構成される地中探査装置である。
また、本発明の第二の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される１対の電圧測定電極８１、８２と、電圧測定電極８１、８２の２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置である。
また、本発明の第三の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置である。
また、本発明の第四の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した少なくとも２以上の総数Ｍ個の電流印加電極３１、３２、・・・３Ｍと、電圧電流変換器２が出力する電流を電流印加電極３１、３２、・・・３Ｍから選択した任意の１対の電極に供給する電流供給選択部１２と、電流印加電極３１、３２、・・・３Ｍの設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置である。
前記の第一から第四の発明の地中探査装置において、発振器１が発生する電圧信号が正弦波、三角波、のこぎり波、あるいは矩形波であるか、前記正弦波、三角波、のこぎり波、あるいは矩形波を基本発振周波数の整数分の１の周波数の２進最長符号系列で２相変調を施したものであることを特徴とする。
前記の第二の発明の地中探査装置において、１対の電圧測定電極８１、８２を電極間距離を一定値に固定して地表面を移動することが可能な機構を具備し、地表面上を移動しながら測定することを特徴とする。
前記の第三および第四の発明の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元マップを表示する表示部１１を備えることを特徴とする。
前記第三および第四の発明の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元補間処理を施し、さらに得られた２次元補間値に対して前記２次元を平面とみなして互いに直交する平面上の２軸ｘとｙに対してそれぞれ軸方向の微分処理を行い、ｘ方向微分値とｙ方向微分値をベクトルの要素とみなした２次元ベクトル表示またはアローマップ表示を行う表示部１１を備えることを特徴とする。
前記の地中探査装置において、電流印加電極３１、３２、・・・３Ｍ、電圧測定電極８１、８２、・・・８Ｎの３次元座標と地形データに基づいて印加電流順方向問題を有限要素法または境界要素法またはインピーダンス網で解いて大地モデルの電圧測定電極に出現する高周波電圧振幅値を求め、実測値と計算値の誤差が小さくなるように大地モデルの局所誘電率、導電率を変更することにより地中物質の推定を行う処理部１０と、処理部１０より入力情報と推定処理結果の供給を受けて２次元あるいは３次元表示処理を行う表示部１１を備えることを特徴とする。
本発明の第五の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電極３と、電極３の両端に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンスの差を演算して発振周波数の違いによるインピーダンス差を計算する演算部６と、演算部６が演算するインピーダンス差に基づいて電極３の２点間のインピーダンスが少なくとも一方の周波数で変化したことを検出判定する判定部７により構成される地中探査システムである。
本発明の第六の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される１対の電圧測定電極８１、８２と、電圧測定電極８１、８２の２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システムである。
本発明の第七の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システムである。
本発明の第八の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した少なくとも２以上の総数Ｍ個の電流印加電極３１、３２、・・・３Ｍと、電圧電流変換器２が出力する電流を電流印加電極３１、３２、・・・３Ｍから選択した任意の１対の電極に供給する電流供給選択部１２と、電流印加電極３１、３２、・・・３Ｍの設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システムである。
本発明の第九の発明は、大地の離れた位置に設置あるいは埋設した１対の電流印加電極または大地の離れた位置に設置あるいは埋設した複数の電流印加電極から選択した任意の１対の電極に異なる２以上の発振周波数の電流を通流する工程と、
大地の離れた位置に設置した１対の電圧測定電極または大地の離れた位置に設置した複数の電圧測定電極から選択した任意の１対の電極該電極間に発生する電圧を測定して２点間のインピーダンスを検出する工程と、
該インピーダンス値を上記発振周波数毎に保持する工程と、
該インピーダンス値の差を演算して発振周波数の違いによるインピーダンス差を演算する工程を有する地中探査方法である。
前記の第九の発明の地中探査方法において、前記インピーダンス差が変化したことを検出判定する工程を有することを特徴とする。
また、本発明の第二の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される１対の電圧測定電極８１、８２と、電圧測定電極８１、８２の２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置である。
また、本発明の第三の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置である。
また、本発明の第四の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した少なくとも２以上の総数Ｍ個の電流印加電極３１、３２、・・・３Ｍと、電圧電流変換器２が出力する電流を電流印加電極３１、３２、・・・３Ｍから選択した任意の１対の電極に供給する電流供給選択部１２と、電流印加電極３１、３２、・・・３Ｍの設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置である。
前記の第一から第四の発明の地中探査装置において、発振器１が発生する電圧信号が正弦波、三角波、のこぎり波、あるいは矩形波であるか、前記正弦波、三角波、のこぎり波、あるいは矩形波を基本発振周波数の整数分の１の周波数の２進最長符号系列で２相変調を施したものであることを特徴とする。
前記の第二の発明の地中探査装置において、１対の電圧測定電極８１、８２を電極間距離を一定値に固定して地表面を移動することが可能な機構を具備し、地表面上を移動しながら測定することを特徴とする。
前記の第三および第四の発明の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元マップを表示する表示部１１を備えることを特徴とする。
前記第三および第四の発明の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元補間処理を施し、さらに得られた２次元補間値に対して前記２次元を平面とみなして互いに直交する平面上の２軸ｘとｙに対してそれぞれ軸方向の微分処理を行い、ｘ方向微分値とｙ方向微分値をベクトルの要素とみなした２次元ベクトル表示またはアローマップ表示を行う表示部１１を備えることを特徴とする。
前記の地中探査装置において、電流印加電極３１、３２、・・・３Ｍ、電圧測定電極８１、８２、・・・８Ｎの３次元座標と地形データに基づいて印加電流順方向問題を有限要素法または境界要素法またはインピーダンス網で解いて大地モデルの電圧測定電極に出現する高周波電圧振幅値を求め、実測値と計算値の誤差が小さくなるように大地モデルの局所誘電率、導電率を変更することにより地中物質の推定を行う処理部１０と、処理部１０より入力情報と推定処理結果の供給を受けて２次元あるいは３次元表示処理を行う表示部１１を備えることを特徴とする。
本発明の第五の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電極３と、電極３の両端に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンスの差を演算して発振周波数の違いによるインピーダンス差を計算する演算部６と、演算部６が演算するインピーダンス差に基づいて電極３の２点間のインピーダンスが少なくとも一方の周波数で変化したことを検出判定する判定部７により構成される地中探査システムである。
本発明の第六の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される１対の電圧測定電極８１、８２と、電圧測定電極８１、８２の２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システムである。
本発明の第七の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システムである。
本発明の第八の発明は、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した少なくとも２以上の総数Ｍ個の電流印加電極３１、３２、・・・３Ｍと、電圧電流変換器２が出力する電流を電流印加電極３１、３２、・・・３Ｍから選択した任意の１対の電極に供給する電流供給選択部１２と、電流印加電極３１、３２、・・・３Ｍの設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システムである。
本発明の第九の発明は、大地の離れた位置に設置あるいは埋設した１対の電流印加電極または大地の離れた位置に設置あるいは埋設した複数の電流印加電極から選択した任意の１対の電極に異なる２以上の発振周波数の電流を通流する工程と、
大地の離れた位置に設置した１対の電圧測定電極または大地の離れた位置に設置した複数の電圧測定電極から選択した任意の１対の電極該電極間に発生する電圧を測定して２点間のインピーダンスを検出する工程と、
該インピーダンス値を上記発振周波数毎に保持する工程と、
該インピーダンス値の差を演算して発振周波数の違いによるインピーダンス差を演算する工程を有する地中探査方法である。
前記の第九の発明の地中探査方法において、前記インピーダンス差が変化したことを検出判定する工程を有することを特徴とする。
図１は、この発明の実施の形態による地中探査装置の構造を説明する図である。（実施例３）
図２は、土壌物質や様々な化学物質ごとに異なった比誘電率を持っていることを説明する図表である。
図３は、誘電率の周波数依存性とその要因を説明する図である。
図４Ａ及び４Ｂは、界面分極と配向分極について説明する図である。
図５Ａ、５Ｂ、５Ｃ及び５Ｄは、印加電流の周波数による誘電率の変化に起因して大地の深部地質の部分的差異によって通流経路が変わることを説明する図である。
図６Ａ、６Ｂ、６Ｃ及び６Ｄは、この発明の実施の形態による地中探査装置の処理の仕組みを説明する図である。（実施例３）
図７は、この発明の実施の形態による地中探査装置の構造を説明する図である。（実施例１）
図８Ａ、８Ｂ及び８Ｃは、この発明の実施の形態による漏水に伴うサンドブラスト現象を検出する動作原理を説明する図である。（実施例１）
図９は、この発明の実施の形態による地中探査装置の発振器１の構造を説明する図である。（実施例２）
図１０は、この発明の実施の形態による地中探査装置の構造を説明する図である。（実施例４）
図１１Ａ及び１１Ｂは、この発明の実施の形態による地中探査装置の大地モデルの例を説明する図である。（実施例３）
図１２Ａ及び１２Ｂは、この発明の実施の形態による地中探査装置の大地モデルの例を説明する図である。（実施例３）
図１３Ａ、１３Ｂ及び１３Ｃは、この発明の一実施の形態による地中探査装置の構造を説明する図である。
図２は、土壌物質や様々な化学物質ごとに異なった比誘電率を持っていることを説明する図表である。
図３は、誘電率の周波数依存性とその要因を説明する図である。
図４Ａ及び４Ｂは、界面分極と配向分極について説明する図である。
図５Ａ、５Ｂ、５Ｃ及び５Ｄは、印加電流の周波数による誘電率の変化に起因して大地の深部地質の部分的差異によって通流経路が変わることを説明する図である。
図６Ａ、６Ｂ、６Ｃ及び６Ｄは、この発明の実施の形態による地中探査装置の処理の仕組みを説明する図である。（実施例３）
図７は、この発明の実施の形態による地中探査装置の構造を説明する図である。（実施例１）
図８Ａ、８Ｂ及び８Ｃは、この発明の実施の形態による漏水に伴うサンドブラスト現象を検出する動作原理を説明する図である。（実施例１）
図９は、この発明の実施の形態による地中探査装置の発振器１の構造を説明する図である。（実施例２）
図１０は、この発明の実施の形態による地中探査装置の構造を説明する図である。（実施例４）
図１１Ａ及び１１Ｂは、この発明の実施の形態による地中探査装置の大地モデルの例を説明する図である。（実施例３）
図１２Ａ及び１２Ｂは、この発明の実施の形態による地中探査装置の大地モデルの例を説明する図である。（実施例３）
図１３Ａ、１３Ｂ及び１３Ｃは、この発明の一実施の形態による地中探査装置の構造を説明する図である。
本発明の発明者の一人である上田智章は、鋭意研究を重ねた結果、大地を構成する各地質構成物質や汚染物質あるいは地下水の電気導電率σは印加する電流の周波数によって殆ど変化しないが、誘電率εは大きく変化し、分布的に存在する静電容量Ｃも大きな変化をすることを見出した。
図２に示すように物質ごとに比誘電率は異なることが知られている。比誘電率とは真空の誘電率に対する比を表している数値である。この誘電率は図３に示すように周波数依存性があることが知られている。分極しやすい物質ほど誘電率は大きい。この分極の原因には、１）空間電荷の移動に伴う界面分極と、２）物質を構成する分子や粒子の分極に関係する双極子分極と、３）原子核の移動によるイオン分極と、４）電子の移動による電子分極があることが知られている。
上記の分極のうち、界面分極と双極子分極を説明する図を図４Ａ及び４Ｂに示す。界面分極とは外部から印加された電界によって粒子と粒子の隙間を荷電粒子が移動する現象である。よって界面分極は地質を構成する土砂や石の粒子径や性質に依存するところが大きいと考えられ、界面分極に依存する誘電率の変化は直流（ＤＣ）から数ｋＨｚ付近の周波数において顕著に現れることが知られている。また、双極子分極は外部から印加された電界によって分子配向が揃うことによって発生する現象であり、地層に含まれる物質の分子や粒子が印加された電界によって受けるモーメント力による回転に依存するところが大きいと考えられる。
それ故、印加する電流の周波数を適切に選択して、２以上の周波数において交流電気試験を行えば、誘電率が地質を構成する物質や含有される物質に依存して変化する周波数が異なるので、物質を特定することが可能なのである。本発明の発明者の一人である上田智章はこの点に着目した。
この発明による地中探査装置によれば、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電極３と、電極３の両端に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンスの差を演算して発振周波数の違いによるインピーダンス差を計算する演算部６と、演算部６が演算するインピーダンス差に基づいて短時間に電極３の２点間のインピーダンスが少なくとも一方の周波数で大きく変化したことを検出判定する判定部７により構成されており、電極の導電率は大地を構成する物質の導電率に比して非常に低いので、２端子法によってもインピーダンス計測が正確に行え、計測結果である複数の周波数によるインピーダンスの差から誘電率の変化を的確に捕捉することができ、ひいては１対の電極３間の最短電流経路近傍に起こる地中の変化を検出することができる。
好ましくは、１対の電極３を埋設ガス管の近傍に配置することである。これによって電極３間の最短電流経路近傍に埋設された水道管からの漏水に伴うサンドブラスト現象がガス管の保全に対して重大な影響を及ぼす恐れがある場合の兆候を漏水部の土壌の誘電率の変化により確実に検出することができる。これは一般の土壌の比誘電率は５から４０であるのに対して、水の比誘電率は８０程度であり、交流電流が漏水部に通流しやすい性質に依る。
また、好ましくは、印加する周波数を１ｋＨｚから１０ＭＨｚまでに設定し、直流を印加しないことである。これによってガス管の塗膜に欠陥がある場合においても腐食の原因となる腐食電流が流れないよう防止することができ、微弱な電流においても大きなインピーダンス変化として捉えることができる。
この発明の他の局面に従うと、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される１対の電圧測定電極８１、８２と、電圧測定電極８１、８２の２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成されている地中探査装置であれば、電圧測定電極８１、８２は地表上の測定点を移動しながら測定することが可能であり、大地表面上でインピーダンスが急激に変化する点を探索することができる。
さらにこの発明の他の局面に従うと、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置であれば、大地表面における電位観測により地表付近だけでなく深部の地中探査も行うことができる。即ち、地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を計測あるいはモニターすることができる。
好ましくは、電圧測定電極８１、８２、・・・８Ｎが２次元的な網状シート上に配置されたものであり、網状シートの敷設によって大地上に電位観測点を容易に展開できるものである。また、電圧測定電極は図１３ＡまたはＢに示すように、電流印加電極と一体的な構造のものであってもよいし、図１３Ｃのように別個に設けてもよい。
さらにこの発明の他の局面に従うと、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した少なくとも２以上の総数Ｍ個の電流印加電極３１、３２、・・・３Ｍと、電圧電流変換器２が出力する電流を電流印加電極３１、３２、・・・３Ｍから選択した任意の１対の電極に供給する電流供給選択部１２と、電流印加電極３１、３２、・・・３Ｍの設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置であれば、通流経路を変更することが容易であるので、地中埋設物や土壌汚染の方向性や広がり方に偏りがある場合においても適切に地中探査を行うことができる。即ち、方向性を持った地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を精確に計測あるいはモニターすることができる。
さらにこの発明の他の局面に従うと、請求項１または請求項２または請求項３または請求項４に記載の地中探査装置において、発振器１が発生する電圧信号が正弦波、三角波、のこぎり波、あるいは矩形波であるか、前記正弦波、三角波、のこぎり波、あるいは矩形波を基本発振周波数の整数分の１の周波数の２進最長符号系列で２相変調を施したものであることを特徴とする地中探査装置であれば、送電線や電車の架線を流れる電流の大地帰還成分や電力機器からの漏えい電流の影響を受けることなく、高いＳ／Ｎ比でノイズに埋もれた微弱な電位を容易に捕捉することができ、強いノイズ耐性の地中探査を行うことが可能である。
好ましくは、矩形波を基本発振周波数の整数分の１の周波数の２進最長符号系列として、長周期のＭ系列符号を用いたものであることである。Ｍ系列符号の有する自己相関特性によって、漏えい電流や大地に蓄積した静電気（自由電荷）の影響を抑制することができる。
さらにこの発明の他の局面に従うと、請求項２に記載の地中探査装置において、１対の電圧測定電極８１、８２を電極間距離を一定値に固定して地表面を移動することが可能な機構を具備し、地表面上を移動しながら測定することを特徴とする地中探査装置であれば、常に電極間距離を一定値に固定した計測を行うことができるので、さらに測定精度を向上することができるだけでなく、測定時間を著しく短縮することができる。
さらにこの発明の他の局面に従うと、請求項３または請求項４に記載の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元マップを表示する表示部１１を備えることを特徴とする地中探査装置であるので、印加周波数が変わったことによって通流経路が変わった部分のみに電流が流れた場合に電位測定面に誘起する高周波電圧の振幅を観測するのと同等の効果を得ることができる。
図５Ａに示すように、導電率がσ１で誘電率がε１の特性を持った大地の地中に、導電率がσ２で誘電率がε２の特性を持った領域が存在している場合、図５Ｂに示すような等価回路モデルが成立する。このとき、大地を抵抗Ｒ１とコンデンサＣ１で表現し、領域を抵抗Ｒ２とコンデンサＣ２で表現することができる。いま、周波数ｆ１においては誘電率ε１より誘電率ε２の方が大きいとすると、周波数ｆ１においては図５Ｃに示すように印加電流は主として領域の方に流れている。周波数ｆ２において領域の方の誘電率ε２が分極要因の１つの追従速度を超えたことにより著しく低下した場合には容量Ｃ２は低下するので領域のインピーダンスは高くなり、図５Ｄに示すように領域内には電流が集中しなくなる。従って、両者のインピーダンス測定値の差を取るということは、領域内に流れていた電流が分散したことによって生じた差異を得ることを意味し、ひいては領域内に集中的に流れていた電流に起因する電位分布パターンを得ることを意味する。この２次元電位分布パターンはダイポール性が高いので、インピーダンスの差のデータに基づいて２次元マップを表示すれば、領域直上付近に図６Ｃに示すようなダイポール性の強い等高線図を得ることができる。即ち、地表表面の観測から地中の深部の領域の検出が行えるという特有の効果を奏することができる。
さらにこの発明の他の局面に従うと、請求項３または請求項４に記載の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元補間処理を施し、さらに得られた２次元補間値に対して前記２次元を平面とみなして互いに直交する平面上の２軸ｘとｙに対してそれぞれ軸方向の微分処理を行い、ｘ方向微分値とｙ方向微分値をベクトルの要素とみなした２次元ベクトル表示またはアローマップ表示を行う表示部１１を備えることを特徴とする地中探査装置であるので、片方の周波数で地中に集中して電流が流れた領域が存在する場合に、インピーダンスの差のデータに基づいて２次元補間処理を施して得られた大地表面の２次元電位マップに基づいて、電流集中部位の直上は変化が激しい性質を有したマップであるために、得られた２次元補間値に対して前記２次元を平面とみなして互いに直交する平面上の２軸ｘとｙに対してそれぞれ軸方向の微分処理を行い、ｘ方向微分値とｙ方向微分値をベクトルの要素とみなした２次元ベクトル表示またはアローマップ表示を行えば、領域内を通流する電流とほぼ同じ方向ベクトルを持った図を得ることができる。即ち、通流電流の密度変化を起こした部位を２次元ベクトル表示またはアローマップ表示という形で可視化することができる。
好ましくは、数式１を用いてベクトルを得るものであることである。
さらにこの発明の他の局面に従うと、請求項７または請求項８に記載の地中探査装置において、電流印加電極３１、３２、・・・３Ｍ、電圧測定電極８１、８２、・・・８Ｎの３次元座標と地形データに基づいて印加電流順方向問題を有限要素法または境界要素法またはインピーダンス網で解いて大地モデルの電圧測定電極に出現する高周波電圧振幅値を求め、実測値と計算値の誤差が小さくなるように大地モデルの局所誘電率、導電率を変更することにより地中物質の推定を行う処理部１０と、処理部１０より入力情報と推定処理結果の供給を受けて２次元あるいは３次元表示処理を行う表示部１１を備えることを特徴とする地中探査装置であれば、大地モデルによる逐次法を用いた逆問題を解くことができるので、大地表面上のおよその領域の位置だけでなく、深さ、領域の大きさ、形状、導電率、誘電率を推定することができる。
導電率および誘電率を推定することができるので、その領域の物質を特定することができるので、漏水状態、地下水脈や埋蔵資源あるいは土壌汚染の位置・形状・深度・収蔵量、土壌改善の進行状況等を観測あるいはモニターすることができる特徴を有する。
好ましくは大地モデルとして、図１１Ａ及び１１Ｂに示すような等価回路網を用いたインピーダンス網モデルか図１２Ａに示すような有限要素法モデル、あるいは図１２Ｂに示すような境界要素法モデルを用いることである。これによって少ないパラメータ数で領域のモデル化が行え、ひいては逆問題を解くための演算時間を著しく短縮することができるのである。
図２に示すように物質ごとに比誘電率は異なることが知られている。比誘電率とは真空の誘電率に対する比を表している数値である。この誘電率は図３に示すように周波数依存性があることが知られている。分極しやすい物質ほど誘電率は大きい。この分極の原因には、１）空間電荷の移動に伴う界面分極と、２）物質を構成する分子や粒子の分極に関係する双極子分極と、３）原子核の移動によるイオン分極と、４）電子の移動による電子分極があることが知られている。
上記の分極のうち、界面分極と双極子分極を説明する図を図４Ａ及び４Ｂに示す。界面分極とは外部から印加された電界によって粒子と粒子の隙間を荷電粒子が移動する現象である。よって界面分極は地質を構成する土砂や石の粒子径や性質に依存するところが大きいと考えられ、界面分極に依存する誘電率の変化は直流（ＤＣ）から数ｋＨｚ付近の周波数において顕著に現れることが知られている。また、双極子分極は外部から印加された電界によって分子配向が揃うことによって発生する現象であり、地層に含まれる物質の分子や粒子が印加された電界によって受けるモーメント力による回転に依存するところが大きいと考えられる。
それ故、印加する電流の周波数を適切に選択して、２以上の周波数において交流電気試験を行えば、誘電率が地質を構成する物質や含有される物質に依存して変化する周波数が異なるので、物質を特定することが可能なのである。本発明の発明者の一人である上田智章はこの点に着目した。
この発明による地中探査装置によれば、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電極３と、電極３の両端に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンスの差を演算して発振周波数の違いによるインピーダンス差を計算する演算部６と、演算部６が演算するインピーダンス差に基づいて短時間に電極３の２点間のインピーダンスが少なくとも一方の周波数で大きく変化したことを検出判定する判定部７により構成されており、電極の導電率は大地を構成する物質の導電率に比して非常に低いので、２端子法によってもインピーダンス計測が正確に行え、計測結果である複数の周波数によるインピーダンスの差から誘電率の変化を的確に捕捉することができ、ひいては１対の電極３間の最短電流経路近傍に起こる地中の変化を検出することができる。
好ましくは、１対の電極３を埋設ガス管の近傍に配置することである。これによって電極３間の最短電流経路近傍に埋設された水道管からの漏水に伴うサンドブラスト現象がガス管の保全に対して重大な影響を及ぼす恐れがある場合の兆候を漏水部の土壌の誘電率の変化により確実に検出することができる。これは一般の土壌の比誘電率は５から４０であるのに対して、水の比誘電率は８０程度であり、交流電流が漏水部に通流しやすい性質に依る。
また、好ましくは、印加する周波数を１ｋＨｚから１０ＭＨｚまでに設定し、直流を印加しないことである。これによってガス管の塗膜に欠陥がある場合においても腐食の原因となる腐食電流が流れないよう防止することができ、微弱な電流においても大きなインピーダンス変化として捉えることができる。
この発明の他の局面に従うと、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される１対の電圧測定電極８１、８２と、電圧測定電極８１、８２の２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成されている地中探査装置であれば、電圧測定電極８１、８２は地表上の測定点を移動しながら測定することが可能であり、大地表面上でインピーダンスが急激に変化する点を探索することができる。
さらにこの発明の他の局面に従うと、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置であれば、大地表面における電位観測により地表付近だけでなく深部の地中探査も行うことができる。即ち、地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を計測あるいはモニターすることができる。
好ましくは、電圧測定電極８１、８２、・・・８Ｎが２次元的な網状シート上に配置されたものであり、網状シートの敷設によって大地上に電位観測点を容易に展開できるものである。また、電圧測定電極は図１３ＡまたはＢに示すように、電流印加電極と一体的な構造のものであってもよいし、図１３Ｃのように別個に設けてもよい。
さらにこの発明の他の局面に従うと、２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２、高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した少なくとも２以上の総数Ｍ個の電流印加電極３１、３２、・・・３Ｍと、電圧電流変換器２が出力する電流を電流印加電極３１、３２、・・・３Ｍから選択した任意の１対の電極に供給する電流供給選択部１２と、電流印加電極３１、３２、・・・３Ｍの設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置であれば、通流経路を変更することが容易であるので、地中埋設物や土壌汚染の方向性や広がり方に偏りがある場合においても適切に地中探査を行うことができる。即ち、方向性を持った地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を精確に計測あるいはモニターすることができる。
さらにこの発明の他の局面に従うと、請求項１または請求項２または請求項３または請求項４に記載の地中探査装置において、発振器１が発生する電圧信号が正弦波、三角波、のこぎり波、あるいは矩形波であるか、前記正弦波、三角波、のこぎり波、あるいは矩形波を基本発振周波数の整数分の１の周波数の２進最長符号系列で２相変調を施したものであることを特徴とする地中探査装置であれば、送電線や電車の架線を流れる電流の大地帰還成分や電力機器からの漏えい電流の影響を受けることなく、高いＳ／Ｎ比でノイズに埋もれた微弱な電位を容易に捕捉することができ、強いノイズ耐性の地中探査を行うことが可能である。
好ましくは、矩形波を基本発振周波数の整数分の１の周波数の２進最長符号系列として、長周期のＭ系列符号を用いたものであることである。Ｍ系列符号の有する自己相関特性によって、漏えい電流や大地に蓄積した静電気（自由電荷）の影響を抑制することができる。
さらにこの発明の他の局面に従うと、請求項２に記載の地中探査装置において、１対の電圧測定電極８１、８２を電極間距離を一定値に固定して地表面を移動することが可能な機構を具備し、地表面上を移動しながら測定することを特徴とする地中探査装置であれば、常に電極間距離を一定値に固定した計測を行うことができるので、さらに測定精度を向上することができるだけでなく、測定時間を著しく短縮することができる。
さらにこの発明の他の局面に従うと、請求項３または請求項４に記載の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元マップを表示する表示部１１を備えることを特徴とする地中探査装置であるので、印加周波数が変わったことによって通流経路が変わった部分のみに電流が流れた場合に電位測定面に誘起する高周波電圧の振幅を観測するのと同等の効果を得ることができる。
図５Ａに示すように、導電率がσ１で誘電率がε１の特性を持った大地の地中に、導電率がσ２で誘電率がε２の特性を持った領域が存在している場合、図５Ｂに示すような等価回路モデルが成立する。このとき、大地を抵抗Ｒ１とコンデンサＣ１で表現し、領域を抵抗Ｒ２とコンデンサＣ２で表現することができる。いま、周波数ｆ１においては誘電率ε１より誘電率ε２の方が大きいとすると、周波数ｆ１においては図５Ｃに示すように印加電流は主として領域の方に流れている。周波数ｆ２において領域の方の誘電率ε２が分極要因の１つの追従速度を超えたことにより著しく低下した場合には容量Ｃ２は低下するので領域のインピーダンスは高くなり、図５Ｄに示すように領域内には電流が集中しなくなる。従って、両者のインピーダンス測定値の差を取るということは、領域内に流れていた電流が分散したことによって生じた差異を得ることを意味し、ひいては領域内に集中的に流れていた電流に起因する電位分布パターンを得ることを意味する。この２次元電位分布パターンはダイポール性が高いので、インピーダンスの差のデータに基づいて２次元マップを表示すれば、領域直上付近に図６Ｃに示すようなダイポール性の強い等高線図を得ることができる。即ち、地表表面の観測から地中の深部の領域の検出が行えるという特有の効果を奏することができる。
さらにこの発明の他の局面に従うと、請求項３または請求項４に記載の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元補間処理を施し、さらに得られた２次元補間値に対して前記２次元を平面とみなして互いに直交する平面上の２軸ｘとｙに対してそれぞれ軸方向の微分処理を行い、ｘ方向微分値とｙ方向微分値をベクトルの要素とみなした２次元ベクトル表示またはアローマップ表示を行う表示部１１を備えることを特徴とする地中探査装置であるので、片方の周波数で地中に集中して電流が流れた領域が存在する場合に、インピーダンスの差のデータに基づいて２次元補間処理を施して得られた大地表面の２次元電位マップに基づいて、電流集中部位の直上は変化が激しい性質を有したマップであるために、得られた２次元補間値に対して前記２次元を平面とみなして互いに直交する平面上の２軸ｘとｙに対してそれぞれ軸方向の微分処理を行い、ｘ方向微分値とｙ方向微分値をベクトルの要素とみなした２次元ベクトル表示またはアローマップ表示を行えば、領域内を通流する電流とほぼ同じ方向ベクトルを持った図を得ることができる。即ち、通流電流の密度変化を起こした部位を２次元ベクトル表示またはアローマップ表示という形で可視化することができる。
好ましくは、数式１を用いてベクトルを得るものであることである。
導電率および誘電率を推定することができるので、その領域の物質を特定することができるので、漏水状態、地下水脈や埋蔵資源あるいは土壌汚染の位置・形状・深度・収蔵量、土壌改善の進行状況等を観測あるいはモニターすることができる特徴を有する。
好ましくは大地モデルとして、図１１Ａ及び１１Ｂに示すような等価回路網を用いたインピーダンス網モデルか図１２Ａに示すような有限要素法モデル、あるいは図１２Ｂに示すような境界要素法モデルを用いることである。これによって少ないパラメータ数で領域のモデル化が行え、ひいては逆問題を解くための演算時間を著しく短縮することができるのである。
以下、この発明の実施の形態を図面を参照して詳しく説明する。なお、図中同一または相当部分には同一符号を付してその説明は繰返さない。
［実施例１］
図２は、この発明の測定原理の基となっている大地を構成する物質や化学物質等の比誘電率が異なる値を有することを説明する図である。比誘電率とは真空誘電率に対して対応する物質の分極しやすさを示す係数である。実際の誘電率εは真空誘電率にその物質の比誘電率を乗じて求めることができる。一般に誘電率は測定する周波数や温度によって異なる。
図３は誘電率の周波数依存性とその変動要因を説明する図である。誘電率は分極し易さを表す係数である。物質の分極要因としては、界面分極、双極子分極、イオン分極、電子分極の４要因が知られている。図４Ａ及び４Ｂは界面分極と双極子分極を説明する図である。
界面分極は図４Ａに示すように粒界の隙間に存在する荷電粒子が外部からの電界によって移動することによって発生する。界面分極は大きな質量の物質が粒界を時間をかけて移動することによって起きる現象であるので、最も低い周波数で発生する分極である。通常、直流（ＤＣ）から数ｋＨｚ程度の周波数で発生し、印加周波数が数ｋＨｚを超えると移動速度が追従しなくなり、この形態での分極は発生しなくなり、誘電率が低下する。粒界の隙間の大きさや移動する荷電粒子の粘性によっても左右され、温度依存性も存在する。従って、界面分極は土壌の性質や含有する荷電粒子の種類に依存するので地質と含有物質を判断する材料となっている。
双極子分極は図４Ｂに示すように荷電粒子を構成する分極した分子が外部の電界によって配向を変えることによって発生する。従って、荷電粒子自体が位置を移動するのに比して短い時間で配向を変えることができるため、追従限界となる周波数は界面分極に比して高く、１ＭＨｚから１００ＭＨｚとなる。この限界周波数は分子量や分子分極の大きさに応じて決まるので、分子を特定する手がかりとなるのである。
図５Ａ乃至５Ｄは本件発明に最も関係が深い測定原理を説明する図である。図５Ａに示すように大地の２点に電流印加電極を設け、２極間に高周波電流を印加するものとする。大地の導電率をσ１とし誘電率をε１とする。また、地中の一部の領域内は周囲の地質とは電気的性質が異なり、導電率をσ２とし誘電率をε２とする。一般に導電率の方は周波数によって殆ど変動しない。図５Ｂに示すようにこの大地と領域はそれぞれ抵抗とコンデンサを並列に接続した等価回路で表すことができる。大地の等価回路の抵抗値をＲ１、コンデンサをＣ１とし、領域の等価回路の抵抗値をＲ２、コンデンサをＣ２とする。
ここで大地を土とし、領域を水がしみこんだ土とする。図２の表に見られるように土の比誘電率は５から４０であり、水の比誘電率は８０程度である。水分子は酸素原子１個に水素原子２個が偏って共有結合したものであるために弱い分極性を有しており、摂氏２０度前後では６個から７個程度の水分子が弱い分極性で紐状に連結した構造を持っている。従って、この短い連結鎖自体も分極しており界面分極の性質がある。低い周波数ｆ１では図５Ｃに示すように水分を含んだ領域の方が周囲の大地に比べて高い誘電率をもっており、それだけこの領域のインピーダンスは低くなるので、印加電流が周囲に比べて高密度に通流する。しかし、界面分極の追従速度を越えた高い周波数ｆ２では界面分極が無くなるために誘電率は低くなり、図５Ｄに示すようにほぼ均質に通流するようになる。
図７に本件発明の１実施形態を示す。１は印加電流の基になる信号を発生する発振器、２は電圧信号を電流に変換する電圧電流変換器、３は電流印加電極、３１、３２、・・・３Ｍは電流印加電極、４は発振器が発生する信号に同期して検波する同期検波によりインピーダンスを測定する測定部、５は印加する周波数ごとに測定データを保持する保持部、６は周波数を変更したことによるインピーダンスの差を演算する演算部、７はインピーダンスの差によってサンドブラストか否かを判定する判定部、８１、８２、・・・８Ｎは電圧測定電極、９は任意の電圧測定電極間の電圧を選択する電圧切換部、１０は大地モデルによる逆問題解析や補間処理等を行う処理部、１１は表示部、１２は電流印加電極を選択して電流供給経路を変更するための電流供給選択部を各々示す。ガス管に近接した位置に埋設された水道管からガス管の方向に急激な漏水が発生した場合に水流とともに土砂が吹きつけガス管を損傷させ、ガス管内に水が入り込む被害が発生する。土砂を巻き込んだ急激な漏水はサンドブラスト現象と呼ばれている。図７は地表からでは発見が遅れるサンドブラスト現象を早期に検知する検出システムの例を示している。
図７において、発振器１は２以上の発振周波数を選択可能な高周波電圧信号発振器であり、電圧信号を電流に変換する電圧電流変換器２に信号を供給している。電圧電流変換器は、サンドブラスト早期警戒対象となるガス管の近傍に埋設された２個の電極３に対して高周波電流を供給している。漏水が発生していない状態においては図７や図８Ａに示すように地中の電極３間を電流がほぼ均質にガス管に沿って流れることになる。電極３間に流れる電流と電極間の土壌のインピーダンスに応じて電極には高周波電圧が発生する。図７に詳細を図示しているように測定部４は発振器１から同期検波用の高周波電圧信号の供給を受けて、電極３間に誘起した電圧を増幅した後で乗算処理を行い、高周波リップル分を低域ろ過器を通すことにより、同期検波処理を行っている。即ち、固定振幅で電流を印加した結果として電極間に発生する高周波電圧の振幅値を得ているので、電極間インピーダンスを測定していることになる。測定部４からのインピーダンス値の供給を受け、保持部５は内部のメモリに印加周波数ごとのインピーダンス値を保持している。保持部５からの測定データの供給を受け、演算部６では印加周波数を変更したことによって発生したインピーダンスの差（変化値）を演算子、判定部７に供給している。判定部７は予め設定された以上のインピーダンスの変化が短時間に起こったか否かを判定することにより、ガス管近傍に近くに埋設されている水道管からの漏水が発生し、サンドブラスト現象を発生する直前状態であることを早期検知することができる。
本実施例の動作をさらに詳細に説明する。図８Ｂはガス管に近接した位置に埋設された水道管が急激な漏水を引き起こしている状態を説明する図である。図８Ｂに図示されているように、電極３に印加された電流はよりインピーダンスが低い領域である漏水箇所に引き込まれやすい。印加電流はガス管に沿って通流されているので、漏水箇所が大きくなり、ガス管に近づくほど電流の引き込みは顕著になるのである。図８Ｃはガス管の切断面から見た電流密度を等高線で示した図である。低い周波数ｆ１において急激にインピーダンスの低下が生じた場合に高い周波数ｆ２においても測定を行い、インピーダンスの差を求めることによりサンドブラスト現象の発生を検知することができるのである。上記の実施例であれば、ガス管に近接して発生した漏水のうち、サンドブラスト現象に移行し、ガス管を損傷する恐れのある場合のみを精確に検出することが可能である。
［実施例２］
図９はこの発明の別の実施の形態による発振器１の詳細構造を説明する図である。地中には交流電気試験にとって障害となる様々な電流が流れている。大地に帯電した静電気の移動に伴う電流、送電線や電車の架線を流れる電流の大地を帰還して流れる電流、電力機器から漏洩する電流等である。これらの影響を受けることなく、微弱なインピーダンスの変化を検出する場合には図９に示すような構造の発振器を用いることが好ましい。
図９に記載されている正弦波発生器は正弦波を発生し、信号を乗算器に供給するとともに、基本周波数に同期してクロック信号をＭ系列発生器に供給している。Ｍ系列発生器とはＬ段のレジスタを用いて構成できる２進符号発生器において最も長周期の２進符号である最長符号系列（Ｍ系列符号）を発生する電子回路である。Ｌ段のレジスタを用いて発生できる最長の符号長は２のＬ乗−１ビット長である。Ｍ系列符号は別名擬似雑音とも呼ばれており、強い自己相関特性を有している。２進論理の“Ｈ”状態に対応して＋１を、“Ｌ”状態に対応して−１を発生する場合には同一符号の０位相で乗算すると、（＋１）×（＋１）＝１となり、（−１）×（−１）＝１となる性質がある。
従って、正弦波にＭ系列を乗じた２相変調された電流信号の印加の結果得られる高周波電圧も印加したＭ系列で２相変調された正弦波と同相か逆相の信号であるため、もう一度送信波を乗算することにより同期検波と同様の復号を行うことができる。しかしながら、受信された高周波信号に直流または交流の障害電流が重畳していても擬似雑音で変調された送信波とは無相関な信号であるために、障害電流成分は擬似雑音により拡散され、Ｍ系列１周期分の相関値を取った場合には素の信号強度の符号長分の１の大きさになるのである。即ち、十分に長い符号長のＭ系列符号を用いることにより、計測に無関係な電流雑音を低減することができる。
上記の構成であれば、非常に大きな環境電流ノイズが存在している大地の地中探査においても障害電流の影響を受けることなく、極めて精確なインピーダンス計測を行うことが可能であり、ひいては微弱な誘電率の変化に伴うインピーダンスの差を求めることができるのである。
［実施例３］
図１はこの発明の別の実施の形態による地中探査装置の構造を示す図である。図１において、発振器１で発生された送信波形は電圧電流変換器２で電流に変換された後で１対の電流印加電極３の両端を介して大地に通流されている。通流の結果、地表に誘起する高周波電圧を測定する目的で設置された総数Ｎ個の電圧測定電極８１、８２，８３，８４、・・・・、８Ｎのうち選択された任意の２電極間に発生した高周波電圧が電圧切換部９を介して増幅器に供給されており、送信波に同期して乗算され、低域ろ過器を通されることにより同期検波処理が施されている。同期検波回路は増幅器と乗算器と低域ろ過きにより構成された測定部４により構成されている。この測定部の出力信号は、送信波の電流振幅が一定値であることによりインピーダンスに比例した値である。測定部４の出力は保持部５において測定周波数ごとに保持され、演算部６で印加周波数を変更したことによって生じたインピーダンスの差が演算され、処理部１０に供給される。処理部１０においては測定データであるインピーダンスの差値に基づいて複数の情報処理を行った後、必要に応じて表示部１１で２次元または３次元表示処理を施すことができる構成となっている。
処理部１０において施される処理についての説明を行う。２以上の周波数でインピーダンスを測定して求めた差に基づいて、電位測定電極を配置した面を平面と考えて、測定値の２次元補間値を演算して得られる図６Ａと図６Ｂの差を演算することには次の意味がある。即ち、周波数ｆ１で測定したインピーダンス分布と周波数ｆ２で測定したインピーダンス分布の差が上記の処理結果と同等であり、周波数の変化に伴って電流が変化した部分のみに変化した電流密度と同値の電流を通流した際に発生する電位パターンを得ることになる。図６Ｃの等高線図は上記の差電流に起因して発生するはずの電位パターンであるため、周波数を変えることで誘電率が大きく変化した領域の直上にダイポールパターンとなって出現する。
この２次元補間値に対して数１による情報処理を施せば、図６Ｄに示すように誘電率が大きく変化した領域の直上に電流の通流方向にほぼ一致するベクトルを持った２次元図を得ることができる。表示部１１において表示することができる。
上記の処理であれば、誘電率が異なる領域の直上位置の同定が容易に行える。
さらに、処理部１０において、大地モデルとして、図１１Ａ及び１１Ｂに図示のインピーダンス網か、図１２Ａに図示の有限要素モデルまたは図１２Ｂに図示の境界要素モデルを用いて、実測定値と大地モデルによる計算値との誤差が小さくなるように逐次法を用いて未知数である導電率、誘電率、３次元座標、大きさ、形状を推定する逆問題を解くことが可能である。
本件発明の発明者の一人である上田智章は、１９９５年に発行された日本生体磁気学会論文誌Ｎｏ．２の７頁から１３頁の「逆問題解析のための生体組織導電率推定の試み」において境界要素法モデルによる逐次法による未知数の推定処理方法について述べている。
上記の構成であれば、大地モデルを用いた地中深部の未知数量を推定することができるため、測定に用いた周波数情報と推定された導電率や誘電率情報により、地質を構成する物質や含有物質の同定が行うことができる。
［実施例４］
図１０は、この発明の別の実施の形態による地中探査装置の構造を示す図である。図１との違いは電流印加電極を３個以上設置した点と、電流を供給するための電流印加電極を選択するための電流供給選択部１２を設けた点のみである。
上記の構成であれば、通流方向を自由に変更することが可能であるために、地質が異なる領域が方向性を有した広がりを持つ場合においても精確にその領域を検出し、その領域の未知数を推定することができ、その領域の地質や含有物質を特定することが可能である。
開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上記した説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
［実施例１］
図２は、この発明の測定原理の基となっている大地を構成する物質や化学物質等の比誘電率が異なる値を有することを説明する図である。比誘電率とは真空誘電率に対して対応する物質の分極しやすさを示す係数である。実際の誘電率εは真空誘電率にその物質の比誘電率を乗じて求めることができる。一般に誘電率は測定する周波数や温度によって異なる。
図３は誘電率の周波数依存性とその変動要因を説明する図である。誘電率は分極し易さを表す係数である。物質の分極要因としては、界面分極、双極子分極、イオン分極、電子分極の４要因が知られている。図４Ａ及び４Ｂは界面分極と双極子分極を説明する図である。
界面分極は図４Ａに示すように粒界の隙間に存在する荷電粒子が外部からの電界によって移動することによって発生する。界面分極は大きな質量の物質が粒界を時間をかけて移動することによって起きる現象であるので、最も低い周波数で発生する分極である。通常、直流（ＤＣ）から数ｋＨｚ程度の周波数で発生し、印加周波数が数ｋＨｚを超えると移動速度が追従しなくなり、この形態での分極は発生しなくなり、誘電率が低下する。粒界の隙間の大きさや移動する荷電粒子の粘性によっても左右され、温度依存性も存在する。従って、界面分極は土壌の性質や含有する荷電粒子の種類に依存するので地質と含有物質を判断する材料となっている。
双極子分極は図４Ｂに示すように荷電粒子を構成する分極した分子が外部の電界によって配向を変えることによって発生する。従って、荷電粒子自体が位置を移動するのに比して短い時間で配向を変えることができるため、追従限界となる周波数は界面分極に比して高く、１ＭＨｚから１００ＭＨｚとなる。この限界周波数は分子量や分子分極の大きさに応じて決まるので、分子を特定する手がかりとなるのである。
図５Ａ乃至５Ｄは本件発明に最も関係が深い測定原理を説明する図である。図５Ａに示すように大地の２点に電流印加電極を設け、２極間に高周波電流を印加するものとする。大地の導電率をσ１とし誘電率をε１とする。また、地中の一部の領域内は周囲の地質とは電気的性質が異なり、導電率をσ２とし誘電率をε２とする。一般に導電率の方は周波数によって殆ど変動しない。図５Ｂに示すようにこの大地と領域はそれぞれ抵抗とコンデンサを並列に接続した等価回路で表すことができる。大地の等価回路の抵抗値をＲ１、コンデンサをＣ１とし、領域の等価回路の抵抗値をＲ２、コンデンサをＣ２とする。
ここで大地を土とし、領域を水がしみこんだ土とする。図２の表に見られるように土の比誘電率は５から４０であり、水の比誘電率は８０程度である。水分子は酸素原子１個に水素原子２個が偏って共有結合したものであるために弱い分極性を有しており、摂氏２０度前後では６個から７個程度の水分子が弱い分極性で紐状に連結した構造を持っている。従って、この短い連結鎖自体も分極しており界面分極の性質がある。低い周波数ｆ１では図５Ｃに示すように水分を含んだ領域の方が周囲の大地に比べて高い誘電率をもっており、それだけこの領域のインピーダンスは低くなるので、印加電流が周囲に比べて高密度に通流する。しかし、界面分極の追従速度を越えた高い周波数ｆ２では界面分極が無くなるために誘電率は低くなり、図５Ｄに示すようにほぼ均質に通流するようになる。
図７に本件発明の１実施形態を示す。１は印加電流の基になる信号を発生する発振器、２は電圧信号を電流に変換する電圧電流変換器、３は電流印加電極、３１、３２、・・・３Ｍは電流印加電極、４は発振器が発生する信号に同期して検波する同期検波によりインピーダンスを測定する測定部、５は印加する周波数ごとに測定データを保持する保持部、６は周波数を変更したことによるインピーダンスの差を演算する演算部、７はインピーダンスの差によってサンドブラストか否かを判定する判定部、８１、８２、・・・８Ｎは電圧測定電極、９は任意の電圧測定電極間の電圧を選択する電圧切換部、１０は大地モデルによる逆問題解析や補間処理等を行う処理部、１１は表示部、１２は電流印加電極を選択して電流供給経路を変更するための電流供給選択部を各々示す。ガス管に近接した位置に埋設された水道管からガス管の方向に急激な漏水が発生した場合に水流とともに土砂が吹きつけガス管を損傷させ、ガス管内に水が入り込む被害が発生する。土砂を巻き込んだ急激な漏水はサンドブラスト現象と呼ばれている。図７は地表からでは発見が遅れるサンドブラスト現象を早期に検知する検出システムの例を示している。
図７において、発振器１は２以上の発振周波数を選択可能な高周波電圧信号発振器であり、電圧信号を電流に変換する電圧電流変換器２に信号を供給している。電圧電流変換器は、サンドブラスト早期警戒対象となるガス管の近傍に埋設された２個の電極３に対して高周波電流を供給している。漏水が発生していない状態においては図７や図８Ａに示すように地中の電極３間を電流がほぼ均質にガス管に沿って流れることになる。電極３間に流れる電流と電極間の土壌のインピーダンスに応じて電極には高周波電圧が発生する。図７に詳細を図示しているように測定部４は発振器１から同期検波用の高周波電圧信号の供給を受けて、電極３間に誘起した電圧を増幅した後で乗算処理を行い、高周波リップル分を低域ろ過器を通すことにより、同期検波処理を行っている。即ち、固定振幅で電流を印加した結果として電極間に発生する高周波電圧の振幅値を得ているので、電極間インピーダンスを測定していることになる。測定部４からのインピーダンス値の供給を受け、保持部５は内部のメモリに印加周波数ごとのインピーダンス値を保持している。保持部５からの測定データの供給を受け、演算部６では印加周波数を変更したことによって発生したインピーダンスの差（変化値）を演算子、判定部７に供給している。判定部７は予め設定された以上のインピーダンスの変化が短時間に起こったか否かを判定することにより、ガス管近傍に近くに埋設されている水道管からの漏水が発生し、サンドブラスト現象を発生する直前状態であることを早期検知することができる。
本実施例の動作をさらに詳細に説明する。図８Ｂはガス管に近接した位置に埋設された水道管が急激な漏水を引き起こしている状態を説明する図である。図８Ｂに図示されているように、電極３に印加された電流はよりインピーダンスが低い領域である漏水箇所に引き込まれやすい。印加電流はガス管に沿って通流されているので、漏水箇所が大きくなり、ガス管に近づくほど電流の引き込みは顕著になるのである。図８Ｃはガス管の切断面から見た電流密度を等高線で示した図である。低い周波数ｆ１において急激にインピーダンスの低下が生じた場合に高い周波数ｆ２においても測定を行い、インピーダンスの差を求めることによりサンドブラスト現象の発生を検知することができるのである。上記の実施例であれば、ガス管に近接して発生した漏水のうち、サンドブラスト現象に移行し、ガス管を損傷する恐れのある場合のみを精確に検出することが可能である。
［実施例２］
図９はこの発明の別の実施の形態による発振器１の詳細構造を説明する図である。地中には交流電気試験にとって障害となる様々な電流が流れている。大地に帯電した静電気の移動に伴う電流、送電線や電車の架線を流れる電流の大地を帰還して流れる電流、電力機器から漏洩する電流等である。これらの影響を受けることなく、微弱なインピーダンスの変化を検出する場合には図９に示すような構造の発振器を用いることが好ましい。
図９に記載されている正弦波発生器は正弦波を発生し、信号を乗算器に供給するとともに、基本周波数に同期してクロック信号をＭ系列発生器に供給している。Ｍ系列発生器とはＬ段のレジスタを用いて構成できる２進符号発生器において最も長周期の２進符号である最長符号系列（Ｍ系列符号）を発生する電子回路である。Ｌ段のレジスタを用いて発生できる最長の符号長は２のＬ乗−１ビット長である。Ｍ系列符号は別名擬似雑音とも呼ばれており、強い自己相関特性を有している。２進論理の“Ｈ”状態に対応して＋１を、“Ｌ”状態に対応して−１を発生する場合には同一符号の０位相で乗算すると、（＋１）×（＋１）＝１となり、（−１）×（−１）＝１となる性質がある。
従って、正弦波にＭ系列を乗じた２相変調された電流信号の印加の結果得られる高周波電圧も印加したＭ系列で２相変調された正弦波と同相か逆相の信号であるため、もう一度送信波を乗算することにより同期検波と同様の復号を行うことができる。しかしながら、受信された高周波信号に直流または交流の障害電流が重畳していても擬似雑音で変調された送信波とは無相関な信号であるために、障害電流成分は擬似雑音により拡散され、Ｍ系列１周期分の相関値を取った場合には素の信号強度の符号長分の１の大きさになるのである。即ち、十分に長い符号長のＭ系列符号を用いることにより、計測に無関係な電流雑音を低減することができる。
上記の構成であれば、非常に大きな環境電流ノイズが存在している大地の地中探査においても障害電流の影響を受けることなく、極めて精確なインピーダンス計測を行うことが可能であり、ひいては微弱な誘電率の変化に伴うインピーダンスの差を求めることができるのである。
［実施例３］
図１はこの発明の別の実施の形態による地中探査装置の構造を示す図である。図１において、発振器１で発生された送信波形は電圧電流変換器２で電流に変換された後で１対の電流印加電極３の両端を介して大地に通流されている。通流の結果、地表に誘起する高周波電圧を測定する目的で設置された総数Ｎ個の電圧測定電極８１、８２，８３，８４、・・・・、８Ｎのうち選択された任意の２電極間に発生した高周波電圧が電圧切換部９を介して増幅器に供給されており、送信波に同期して乗算され、低域ろ過器を通されることにより同期検波処理が施されている。同期検波回路は増幅器と乗算器と低域ろ過きにより構成された測定部４により構成されている。この測定部の出力信号は、送信波の電流振幅が一定値であることによりインピーダンスに比例した値である。測定部４の出力は保持部５において測定周波数ごとに保持され、演算部６で印加周波数を変更したことによって生じたインピーダンスの差が演算され、処理部１０に供給される。処理部１０においては測定データであるインピーダンスの差値に基づいて複数の情報処理を行った後、必要に応じて表示部１１で２次元または３次元表示処理を施すことができる構成となっている。
処理部１０において施される処理についての説明を行う。２以上の周波数でインピーダンスを測定して求めた差に基づいて、電位測定電極を配置した面を平面と考えて、測定値の２次元補間値を演算して得られる図６Ａと図６Ｂの差を演算することには次の意味がある。即ち、周波数ｆ１で測定したインピーダンス分布と周波数ｆ２で測定したインピーダンス分布の差が上記の処理結果と同等であり、周波数の変化に伴って電流が変化した部分のみに変化した電流密度と同値の電流を通流した際に発生する電位パターンを得ることになる。図６Ｃの等高線図は上記の差電流に起因して発生するはずの電位パターンであるため、周波数を変えることで誘電率が大きく変化した領域の直上にダイポールパターンとなって出現する。
この２次元補間値に対して数１による情報処理を施せば、図６Ｄに示すように誘電率が大きく変化した領域の直上に電流の通流方向にほぼ一致するベクトルを持った２次元図を得ることができる。表示部１１において表示することができる。
上記の処理であれば、誘電率が異なる領域の直上位置の同定が容易に行える。
さらに、処理部１０において、大地モデルとして、図１１Ａ及び１１Ｂに図示のインピーダンス網か、図１２Ａに図示の有限要素モデルまたは図１２Ｂに図示の境界要素モデルを用いて、実測定値と大地モデルによる計算値との誤差が小さくなるように逐次法を用いて未知数である導電率、誘電率、３次元座標、大きさ、形状を推定する逆問題を解くことが可能である。
本件発明の発明者の一人である上田智章は、１９９５年に発行された日本生体磁気学会論文誌Ｎｏ．２の７頁から１３頁の「逆問題解析のための生体組織導電率推定の試み」において境界要素法モデルによる逐次法による未知数の推定処理方法について述べている。
上記の構成であれば、大地モデルを用いた地中深部の未知数量を推定することができるため、測定に用いた周波数情報と推定された導電率や誘電率情報により、地質を構成する物質や含有物質の同定が行うことができる。
［実施例４］
図１０は、この発明の別の実施の形態による地中探査装置の構造を示す図である。図１との違いは電流印加電極を３個以上設置した点と、電流を供給するための電流印加電極を選択するための電流供給選択部１２を設けた点のみである。
上記の構成であれば、通流方向を自由に変更することが可能であるために、地質が異なる領域が方向性を有した広がりを持つ場合においても精確にその領域を検出し、その領域の未知数を推定することができ、その領域の地質や含有物質を特定することが可能である。
開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上記した説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
以上のように、この発明によれば、低侵襲性の交流電気試験によって地質やその含有物質を調べる地中探査を環境電流ノイズの影響を受けることなく精確に行うことができるという特有の効果を奏する。
また、この発明によれば、ガス管に重大な損傷を与える恐れのある近接して埋設された水道管の漏水に起因して発生するサンドブラスト現象を早期に検出することができるという特有の効果を奏する。
さらに、この発明によれば、大地モデルによる未知数量の推定処理により、地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を計測あるいはモニターすることが可能であるという特有の効果を奏する。
また、汚染物質や資源の埋蔵位置のおよその位置を地表面上に配置した多数の電位測定電極の配置面を平面と仮定して補間処理を施した２次元分布図に基づいて作成されるベクトル図あるいはアローマップを使用して容易に同定することができるという特有の効果を奏する。
また、この発明によれば、ガス管に重大な損傷を与える恐れのある近接して埋設された水道管の漏水に起因して発生するサンドブラスト現象を早期に検出することができるという特有の効果を奏する。
さらに、この発明によれば、大地モデルによる未知数量の推定処理により、地下水脈や埋蔵資源あるいは土壌汚染の位置・深度・収蔵量、地層に含まれる物質の特定、あるいは土壌改善の進行状況を計測あるいはモニターすることが可能であるという特有の効果を奏する。
また、汚染物質や資源の埋蔵位置のおよその位置を地表面上に配置した多数の電位測定電極の配置面を平面と仮定して補間処理を施した２次元分布図に基づいて作成されるベクトル図あるいはアローマップを使用して容易に同定することができるという特有の効果を奏する。
Claims (15)
- ２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電極３と、電極３の両端に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンスの差を演算して発振周波数の違いによるインピーダンス差を計算する演算部６と、演算部６が演算するインピーダンス差に基づいて電極３の２点間のインピーダンスが少なくとも一方の周波数で変化したことを検出判定する判定部７により構成される地中探査装置。
- ２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される１対の電圧測定電極８１、８２と、電圧測定電極８１、８２の２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置。
- ２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置。
- ２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、高周波電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した少なくとも２以上の総数Ｍ個の電流印加電極３１、３２、・・・３Ｍと、電圧電流変換器２が出力する電流を電流印加電極３１、３２、・・・３Ｍから選択した任意の１対の電極に供給する電流供給選択部１２と、電流印加電極３１、３２、・・・３Ｍの設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する高周波電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査装置。
- 請求項１〜４のいずれか一項に記載の地中探査装置において、発振器１が発生する電圧信号が正弦波、三角波、のこぎり波、あるいは矩形波であるか、前記正弦波、三角波、のこぎり波、あるいは矩形波を基本発振周波数の整数分の１の周波数の２進最長符号系列で２相変調を施したものであることを特徴とする地中探査装置。
- 請求項２に記載の地中探査装置において、１対の電圧測定電極８１、８２を電極間距離を一定値に固定して地表面を移動することが可能な機構を具備し、地表面上を移動しながら測定することを特徴とする地中探査装置。
- 請求項３または４に記載の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元マップを表示する表示部１１を備えることを特徴とする地中探査装置。
- 請求項３または４に記載の地中探査装置において、総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎは地表面上に２次元的に配置されたものであり、演算部６から供給される印加発振周波数の違いによるインピーダンスの差のデータに基づいて２次元補間処理を施し、さらに得られた２次元補間値に対して前記２次元を平面とみなして互いに直交する平面上の２軸ｘとｙに対してそれぞれ軸方向の微分処理を行い、ｘ方向微分値とｙ方向微分値をベクトルの要素とみなした２次元ベクトル表示またはアローマップ表示を行う表示部１１を備えることを特徴とする地中探査装置。
- 請求項７または８に記載の地中探査装置において、電流印加電極３１、３２、・・・３Ｍ、電圧測定電極８１、８２、・・・８Ｎの３次元座標と地形データに基づいて印加電流順方向問題を有限要素法または境界要素法またはインピーダンス網で解いて大地モデルの電圧測定電極に出現する高周波電圧振幅値を求め、実測値と計算値の誤差が小さくなるように大地モデルの局所誘電率、導電率を変更することにより地中物質の推定を行う処理部１０と、処理部１０より入力情報と推定処理結果の供給を受けて２次元あるいは３次元表示処理を行う表示部１１を備えることを特徴とする地中探査装置。
- ２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電極３と、電極３の両端に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンスの差を演算して発振周波数の違いによるインピーダンス差を計算する演算部６と、演算部６が演算するインピーダンス差に基づいて電極３の２点間のインピーダンスが少なくとも一方の周波数で変化したことを検出判定する判定部７により構成される地中探査システム。
- ２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される１対の電圧測定電極８１、８２と、電圧測定電極８１、８２の２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システム。
- ２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電圧電流変換器２が供給する電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した１対の電流印加電極３と、電流印加電極３の設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システム。
- ２以上の基本発振周波数を選択可能な発振器１と、発振器１が発生した電圧信号を電流に変換する電圧電流変換器２と、電流を大地に通流するために用い大地の離れた位置に設置あるいは埋設した少なくとも２以上の総数Ｍ個の電流印加電極３１、３２、・・・３Ｍと、電圧電流変換器２が出力する電流を電流印加電極３１、３２、・・・３Ｍから選択した任意の１対の電極に供給する電流供給選択部１２と、電流印加電極３１、３２、・・・３Ｍの設置場所とは異なり互いに離れた位置に設置される２以上の総数Ｎ個の電圧測定電極８１、８２、・・・８Ｎと、Ｎ個の電圧測定電極８１、８２、・・・８Ｎから任意の１対の電極を選択して電圧を供給する電極切換部９と、電圧切換部９が選択した２電極間に発生する電圧を測定して２点間のインピーダンスを検出する測定部４と、測定部４が出力するインピーダンス値を発振周波数毎に保持する保持部５と、保持部５が保持するインピーダンス値に基づいて印加発振周波数の違いによるインピーダンスの差を演算する演算部６により構成される地中探査システム。
- 大地の離れた位置に設置あるいは埋設した１対の電流印加電極または大地の離れた位置に設置あるいは埋設した複数の電流印加電極から選択した任意の１対の電極に異なる２以上の発振周波数の電流を通流する工程と、
大地の離れた位置に設置した１対の電圧測定電極または大地の離れた位置に設置した複数の電圧測定電極から選択した任意の１対の電極該電極間に発生する電圧を測定して２点間のインピーダンスを検出する工程と、
該インピーダンス値を上記発振周波数毎に保持する工程と、
該インピーダンス値の差を演算して発振周波数の違いによるインピーダンス差を演算する工程を有する地中探査方法。 - 前記インピーダンス差が変化したことを検出判定する工程を有する請求項１４に記載の地中探査方法。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003197537 | 2003-07-16 | ||
| JP2003197537 | 2003-07-16 | ||
| PCT/JP2004/010446 WO2005015262A1 (ja) | 2003-07-16 | 2004-07-15 | 地中探査装置、システム及び方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2005015262A1 JPWO2005015262A1 (ja) | 2006-10-05 |
| JP4311745B2 true JP4311745B2 (ja) | 2009-08-12 |
Family
ID=34131357
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005512913A Expired - Fee Related JP4311745B2 (ja) | 2003-07-16 | 2004-07-15 | 地中探査装置、システム及び方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US7034539B2 (ja) |
| EP (1) | EP1655621A1 (ja) |
| JP (1) | JP4311745B2 (ja) |
| WO (1) | WO2005015262A1 (ja) |
Families Citing this family (34)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4311745B2 (ja) | 2003-07-16 | 2009-08-12 | キヤノン株式会社 | 地中探査装置、システム及び方法 |
| EP1718996A1 (en) * | 2004-02-26 | 2006-11-08 | ExxonMobil Upstream Research Company | Electrode configurations for suppression of electroseismic source noise |
| JP4817742B2 (ja) * | 2004-10-29 | 2011-11-16 | キヤノン株式会社 | 土壌中の汚染物質の高濃度域方向検知方法 |
| JP2006126120A (ja) * | 2004-11-01 | 2006-05-18 | Canon Inc | 地中探査装置 |
| US7109716B2 (en) * | 2004-11-24 | 2006-09-19 | Operations Technology Development, Nfp | Soil and time modeling for improved electromagnetic locators for underground utilities |
| JP4710061B2 (ja) * | 2005-02-04 | 2011-06-29 | 株式会社竹中工務店 | コンクリート含有成分測定装置および測定方法 |
| GB0516153D0 (en) * | 2005-08-05 | 2005-09-14 | Mtem Ltd | Multi-transient dc resistivity measurements |
| US7775146B1 (en) | 2006-08-02 | 2010-08-17 | Xtreme Ads Limited | System and method for neutralizing explosives and electronics |
| EP1947480B1 (en) * | 2007-01-22 | 2015-10-07 | Services Pétroliers Schlumberger | A method and apparatus for electrical investigation of a borehole |
| JP4989310B2 (ja) * | 2007-05-22 | 2012-08-01 | 日本電信電話株式会社 | 鉄筋診断装置 |
| JP4206421B2 (ja) * | 2007-06-08 | 2009-01-14 | 有限会社地圏探査技術研究所 | 擬似ノイズ信号を通電波形とした多点同時通電方式による効率的比抵抗計測法 |
| US7863902B2 (en) * | 2008-06-30 | 2011-01-04 | Baker Hughes Incorporated | Coherent in-phase demodulation for resistivity image tools with spread spectrum continuous phase modulation stimulus |
| DE102009043972A1 (de) * | 2009-09-10 | 2011-03-17 | Bucyrus Europe Gmbh | Sensoreinrichtung und Verfahren zur geoelektrischen Erkundung von mineralischen Rohstofflagerstätten |
| JP2011247678A (ja) * | 2010-05-25 | 2011-12-08 | Railway Technical Research Institute | 地下水位観測システム、地下水位観測方法 |
| US9151861B2 (en) | 2011-03-02 | 2015-10-06 | Multi-Phase Technologies, Llc | Method and apparatus for measuring the electrical impedance properties of geological formations using multiple simultaneous current sources |
| JP6046927B2 (ja) | 2011-08-09 | 2016-12-21 | キヤノン株式会社 | 画像処理装置及びその制御方法 |
| US8683907B1 (en) | 2011-09-07 | 2014-04-01 | Xtreme Ads Limited | Electrical discharge system and method for neutralizing explosive devices and electronics |
| US9243874B1 (en) | 2011-09-07 | 2016-01-26 | Xtreme Ads Limited | Electrical discharge system and method for neutralizing explosive devices and electronics |
| US9256003B2 (en) * | 2013-01-15 | 2016-02-09 | Shan Dong University | Three-dimensional focusing induced polarization equipment for advanced geological prediction of water inrush disaster source in underground engineering |
| WO2014116673A1 (en) * | 2013-01-22 | 2014-07-31 | Cidra Corporate Services Inc. | Tomographic imaging of multiphase flows and liquid/froth interfaces using code division multiplexed electrical impedance sensing |
| CN104122593B (zh) * | 2013-04-26 | 2016-10-19 | 邓友明 | 一种对勘探测井的电容率频散测量方法和应用方法 |
| KR101521473B1 (ko) * | 2014-06-02 | 2015-05-21 | 한국지질자원연구원 | 수중탐지장치 및 수중탐지방법 |
| CN104215669B (zh) * | 2014-09-15 | 2016-08-24 | 国家电网公司 | 变压器绝缘系统的含水量测量方法 |
| US9945942B2 (en) | 2015-03-24 | 2018-04-17 | Utilis Israel Ltd. | System and method of underground water detection |
| US10884128B2 (en) | 2015-03-24 | 2021-01-05 | Utilis Israel Ltd. | System and method of underground water detection |
| US10514341B2 (en) * | 2015-03-24 | 2019-12-24 | Utilis Israel Ltd. | System and method of detecting underground gas leakage |
| RU2632998C1 (ru) * | 2016-07-29 | 2017-10-11 | Общество с ограниченной ответственностью "Геологический центр СПбГУ" | Способ выявления загрязнений в почвах и грунтовых водах |
| JP6834593B2 (ja) * | 2017-02-28 | 2021-02-24 | 株式会社Ihi | 計測装置および計測方法 |
| CN107102363A (zh) * | 2017-05-02 | 2017-08-29 | 电子科技大学 | 一种金属物体探测装置及金属探测方法 |
| CN110596623A (zh) * | 2019-09-05 | 2019-12-20 | 国网内蒙古东部电力有限公司检修分公司 | 一种基于混合土壤模型的接地极环境与入地电流测量平台 |
| CN111123372B (zh) * | 2019-12-24 | 2020-11-17 | 中国矿业大学(北京) | 一种煤矿底板钻孔激发极化法突水监测预报系统及方法 |
| CN111487618B (zh) * | 2020-06-11 | 2022-11-11 | 中国地质大学（北京） | 应用于探地雷达的地表可重构阻抗匹配方法及装置 |
| WO2023149579A1 (ja) * | 2022-02-07 | 2023-08-10 | 坂田電機株式会社 | 内部構造推定システム、内部構造推定装置、内部構造推定方法及びプログラム |
| CN116717242B (zh) * | 2023-05-31 | 2024-02-02 | 中国地质大学（武汉） | 电容梯级投切型油井油水界面实时测量系统及方法 |
Family Cites Families (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US1808397A (en) * | 1927-12-02 | 1931-06-02 | Oil Finding Corp | Method of and apparatus for locating deposits of oil, gas and other dielectric subterranean bodies |
| JPS52801B2 (ja) * | 1974-01-25 | 1977-01-11 | ||
| US4714889A (en) * | 1985-07-31 | 1987-12-22 | Chevron Research Company | Method of interpreting impedance distribution of an earth formation using precursor logging data provided by a multi-electrode logging array stationary within a borehole |
| JP2883087B2 (ja) * | 1988-09-01 | 1999-04-19 | 日本電信電話株式会社 | 大地抵抗率測定装置 |
| JPH079465B2 (ja) | 1991-10-11 | 1995-02-01 | 株式会社ダイヤコンサルタント | 自動電気探査方法 |
| JPH0694568A (ja) | 1992-09-11 | 1994-04-05 | Maruyama Seisakusho:Yugen | 漏水検出方法および漏水検出装置 |
| JPH06130156A (ja) | 1992-10-19 | 1994-05-13 | Nippon Telegr & Teleph Corp <Ntt> | 金属埋設物探知方法及び装置 |
| JP2790597B2 (ja) | 1993-06-28 | 1998-08-27 | 鹿島建設株式会社 | 比抵抗利用の地下流れ特性測定方法及び装置 |
| JP3488559B2 (ja) * | 1995-11-01 | 2004-01-19 | 応用地質株式会社 | キャパシタ法によるマルチチャンネル電気探査システム |
| JP3695676B2 (ja) | 1997-04-17 | 2005-09-14 | 応用地質株式会社 | キャパシタ法による多チャンネル電気探査方法 |
| JP2935229B1 (ja) * | 1998-03-04 | 1999-08-16 | 日本電信電話株式会社 | 接地抵抗測定装置及びこれを用いた接地抵抗測定方法 |
| US6373254B1 (en) * | 1998-06-05 | 2002-04-16 | Schlumberger Technology Corporation | Method and apparatus for controlling the effect of contact impedance on a galvanic tool in a logging-while-drilling application |
| JP2000019158A (ja) | 1998-06-30 | 2000-01-21 | Nkk Corp | 埋設鋼管の塗膜損傷位置検出方法及び装置 |
| JP2000019159A (ja) | 1998-06-30 | 2000-01-21 | Nkk Corp | 埋設鋼管の塗膜損傷位置検出方法及び装置 |
| JP4130010B2 (ja) | 1998-07-09 | 2008-08-06 | 三井造船株式会社 | 埋設物探査方法および空洞検知用非接触高周波電流センサ |
| JP2000088691A (ja) | 1998-09-09 | 2000-03-31 | Nishimatsu Constr Co Ltd | 漏水測定装置及び漏水測定方法 |
| JP3625257B2 (ja) | 1999-03-09 | 2005-03-02 | 富士電機ホールディングス株式会社 | 水道管の漏水検知方法および装置 |
| JP3611019B2 (ja) | 1999-07-01 | 2005-01-19 | Ｊｆｅエンジニアリング株式会社 | 埋設塗覆装鋼管の塗膜損傷位置検出方法 |
| JP2001074850A (ja) | 1999-09-09 | 2001-03-23 | Oyo Corp | 無分極電極を用いる電気探査方法 |
| JP2002071822A (ja) | 2000-09-04 | 2002-03-12 | Nobuhiko Oyama | 埋設管の位置検出方法 |
| JP2002156460A (ja) | 2000-11-20 | 2002-05-31 | Sangaku Renkei Kiko Kyushu:Kk | 電気探査方法とそれを用いた電気探査装置及び地雷探知装置 |
| JP4106202B2 (ja) | 2001-06-22 | 2008-06-25 | 新日鉄エンジニアリング株式会社 | 積算手段を用いた埋設金属管類の防食被覆損傷位置検出方法 |
| JP4044303B2 (ja) | 2001-06-25 | 2008-02-06 | 新日鉄エンジニアリング株式会社 | ２種類の周波数信号を用いた埋設金属管類の防食被覆損傷位置検出方法 |
| US6765386B2 (en) * | 2002-04-10 | 2004-07-20 | Halliburton Energy Services, Inc. | Galvanic method of measuring electrical anisotropy |
| JP4311745B2 (ja) | 2003-07-16 | 2009-08-12 | キヤノン株式会社 | 地中探査装置、システム及び方法 |
-
2004
- 2004-07-15 JP JP2005512913A patent/JP4311745B2/ja not_active Expired - Fee Related
- 2004-07-15 EP EP04747835A patent/EP1655621A1/en not_active Withdrawn
- 2004-07-15 WO PCT/JP2004/010446 patent/WO2005015262A1/ja active Application Filing
- 2004-11-19 US US10/991,404 patent/US7034539B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPWO2005015262A1 (ja) | 2006-10-05 |
| EP1655621A1 (en) | 2006-05-10 |
| US20050093548A1 (en) | 2005-05-05 |
| US7034539B2 (en) | 2006-04-25 |
| WO2005015262A1 (ja) | 2005-02-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4311745B2 (ja) | 地中探査装置、システム及び方法 | |
| US9891339B2 (en) | Method and apparatus for detecting and mapping subsurface anomalies | |
| CN110940725A (zh) | 一种堤防内部渗漏通道的探测装置及方法 | |
| Huang et al. | Conductivity-depth imaging of helicopter-borne TEM data based on a pseudolayer half-space model | |
| Bechtel et al. | Geophysical methods | |
| MX2011009449A (es) | Topografia electromagnetica usando campos electromagneticos que se presentan de manera natural como una fuente. | |
| Chang et al. | An investigation into the debris flow induced by Typhoon Morakot in the Siaolin Area, Southern Taiwan, using the electrical resistivity imaging method | |
| RU2545472C2 (ru) | Усовершенствованный способ электромагнитного исследования с управляемым источником | |
| Xue et al. | Discovery of a major coal deposit in China with the use of a modified CSAMT method | |
| Jessop et al. | Magnetometric resistivity: a new approach and its application to the detection of preferential flow paths in mine waste rock dumps | |
| Schneider et al. | Inversion of geo-magnetic full-tensor gradiometer data | |
| Won | Small frequency-domain electromagnetic induction sensors: How in the world does a small broadband EMI sensor with little or no source-receiver separation work? | |
| Braun et al. | Evaluation of the influence of 2-D electrical resistivity on magnetic resonance sounding | |
| Guillemoteau et al. | Evaluation of a rapid hybrid spectral-spatial domain 3D forward-modeling approach for loop-loop electromagnetic induction quadrature data acquired in low-induction-number environments | |
| Tezkan | Radiomagnetotellurics | |
| Huang et al. | 2D Magnetotelluric forward modelling for deep buried water-rich fault and its application | |
| Liu et al. | Mapping water pipeline leakage by ground-penetrating radar diffraction imaging | |
| RU2328020C2 (ru) | Комплексный способ обнаружения неметаллических трубопроводов и повреждений на них | |
| Galagedara et al. | Assessment of soil moisture content measured by borehole GPR and TDR under transient irrigation and drainage | |
| CN106770665B (zh) | 一种基于瞬变电磁法的接地网成像方法 | |
| Sadatcharam et al. | Depth sensitivity of apparent magnetic susceptibility measurements using multi-coil and multi-frequency electromagnetic induction | |
| Goddard et al. | Detection and location of underground cables using magnetic field measurements | |
| Shubitidze et al. | Linear Current Sensing for Detecting and Locating Underground Structures | |
| Gelius | Multi-component processing of sea bed logging data | |
| He et al. | A fusion approach for suppression of environmental noise in spread spectrum induced polarization data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070713 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070713 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070713 |
|
| 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: 20090508 |
|
| 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: 20090511 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120522 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |