(57)【要約】【課題】地層の比抵抗が小さい場合に、ケーシングパイプ１０を用いる流電電位法による比抵抗測定方法の探査精度を改善する。【解決手段】ケーシングパイプ１０内に電流源を挿入して、電流源深度を変更しつつ又は／及び地表測線に沿う付加電流源を設けて及び／又はケーシングパイプ１０内に一定間隔で電流源を配設してケーシングパイプ１０に通電することで比抵抗トモグラフィ測定を実施すると共に、そのデータ解析においては、線要素およびジョイント要素を使ってケーシングパイプ１０の挙動を適切に評価しつつ解析を行うか、別途ケーシングパイプ１０内で計測した電流分布のデータを基に解析を行うかいずれか／又は両方を組み合わせた方法によって、地下の２次元又は３次元的な比抵抗分布を精度よく求める。【選択図】        図１
2007285729.tif 000002
【特許請求の範囲】
【請求項１】
導電性のケーシングパイプを有するボーリング孔中に電極を挿入し、ボーリング孔内の電極と地表に設置した電極アレイ間又は、ボーリング孔中の電極と他のボーリング孔中の電極間で比抵抗トモグラフィ測定を実施して、地下の２次元的又は３次元的な比抵抗分布を求めることを特徴とする地層内の比抵抗測定方法。
【請求項２】
前記ケーシングパイプへの通電を電流源深度を変更しつつ実施することを特徴とする請求項１記載の地層内の比抵抗測定方法。
【請求項３】
前記ケーシングパイプをＦＥＭの線要素及び接地抵抗要素を使ってモデル化してケーシングパイプの影響を適切に評価し、ケーシングパイプ上の電流密度を評価することによって、比抵抗の解析精度を向上することを特徴とする請求項１又は２記載の地層内の比抵抗測定方法。
【請求項４】
ケーシングパイプ電流分布測定装置を用いて計測した電流密度分布を使い、測定されたケーシングパイプから流出する電流を点電流源又は微小線電流源の重ね合わせとして計算し、ＦＥＭ又は差分法によって解析することを特徴とする請求項１又は２記載の地層内の比抵抗測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、比抵抗トモグラフィを用いて地層内の比抵抗を測定する方法に関する。さらに詳しくは、鋼製など導電性のケーシングパイプで保孔されたボーリング孔を使って、ケーシングパイプ内の多数の深度に電流源又は電位測定点をおいて、地下の比抵抗分布を測定する方法に関する。
【０００２】
なお、本発明においてケーシングパイプとは、ボーリング孔内壁を被覆するケーシングチューブパイプを言う。
【背景技術】
【０００３】
比抵抗トモグラフィは、ボーリング孔内および地表に電極を設置し、任意の１つの電極に電流を通電することによって励起された電位を他の電極で計測することを繰り返し、多数の電極の組み合わせによって測定されたデータから地下の２次元的又は３次元的な比抵抗の分布を推定する探査手法である。
【０００４】
一般に、ボーリング孔内に電極を設置して比抵抗トモグラフィを実施する場合は、裸孔（ケーシングパイプが入っていないボーリング孔）であることを測定の前提条件とし、塩化ビニル管等の絶縁性のケーシングパイプが挿入されている場合は、ケーシングパイプにストレイナーと呼ばれるスリット又は円孔をあけて、周囲の地盤や岩盤とボーリング孔内を、水を介入して電気的につながった状態にして測定を行う。
【０００５】
しかし、ボーリング孔に鋼製のケーシングパイプが入っている場合は、ケーシングパイプ自体が一つの線電極として振舞うため通常の比抵抗トモグラフィを適用することはできなかった。このような場合には、ケーシングパイプを１つの線電極として使った探査法が開発されており、流電電位法、鉱体流電法、Ｍｉｓｅ−ａ−ｌａ−Ｍａｓｓｅ  Ｍｅｔｈｏｄ等の様々な名称で呼ばれている。
【０００６】
最近、石油分野において、石油の回収率を向上させるためボーリング孔から水中に水や水蒸気、ＣＯ２等を注入して、これによって移動した石油を別の生産井にて回収するＥＯＲ（Ｅｎｈａｎｃｅｄ  Ｏｉｌ  Ｒｅｃｏｖｅｒｙ）と呼ばれる石油生産技術がある。この技術においても、流電電位法は石油貯留槽のモニタリング・ツールとして用いることが研究されている。
【０００７】
従来の流電電位法では、ケーシングパイプを導電率が無限大の理想的な線電極と考え、その応答を均質な半無限媒質中の理論解で近似して解析していた。しかし、石油を対象とする探査では、ケーシングパイプを囲繞する区域の比抵抗（バックグランド比抵抗）が低いことが多く、流電電位法の適用に際してもケーシングパイプ自体の抵抗を無視することができなくなり、特に電流がケーシングパイプの先端まで到達しなくなることが分かってきた。その場合には流電電位法を適用してもケーシングパイプ先端付近の地下深部の情報を得ることは難しい。
【０００８】
またシミュレーション検討によると、バックグランド比抵抗が小さい場合には、ケーシングパイプ内の各深度位置に電流源を配置すると、電流源の深度によって異なった電位分布が測定されることも分かってきた。この性質を利用すると、ケーシングパイプ内の異なった多数の深度に測定点を配置することにより、従来は不可能とされた鋼製のケーシングパイプを有するボーリング孔における比抵抗トモグラフィの適用が可能となり、地下の比抵抗分布を精度よく求めることができる。本発明技術は従来の比抵抗トモグラフィと流電電位法の間を結ぶ技術である。
【０００９】
また、データ解析においては、ケーシングパイプをＦＥＭ（有限要素法）の線要素で近似する他、周辺地盤とケーシングパイプとの接地抵抗やケーシングパイプ径の影響をＦＥＭのジョイント要素を使って考慮することで、ケーシングパイプ周辺の電位応答を正確に計算する方法が有効である（例えば、非特許文献１参照。）。
【非特許文献１】物理探査学会  第９７回  学術講演会論文集（１９９７）  ｐ８４〜８８、        『比抵抗トモグラフィにおける孔内水の影響』
【発明の開示】
【発明が解決しようとする課題】
【００１０】
本発明は導電性のケーシングパイプが挿入された大深度ボーリングを使って比抵抗トモグラフィを適用することで、従来の流電電位法の探査精度、特に鉛直方向の分解能を改善し、例えば効率的な石油回収のためのＥＯＲモニタリングやＣＯ２地下貯留におけるモニタリング・ツールとして比抵抗測定方法の技術的基礎を確立することを目的とする。
【００１１】
なお、ケーシングパイプの内部に電流源をおいた時の電位分布を、シミュレーションする技術は、本探査手法の実現には欠かせない要素であり、これを確立することを課題とする。
【課題を解決するための手段】
【００１２】
本発明は、上記問題点を解決するためになされたもので、次の技術手段を講じたことを特徴とする地層内の比抵抗測定方法である。
【００１３】
すなわち、本発明は、導電性のケーシングパイプを有するボーリング孔中に電極を挿入し、ボーリング孔内の電極と地表に設置した電極アレイ間又は、ボーリング孔中の電極と他のボーリング孔中の電極間で比抵抗トモグラフィ測定を実施して、地下の２次元的又は３次元的な比抵抗分布を求めることを特徴とする地層内の比抵抗測定方法である。
【００１４】
本発明は、比抵抗を測定するに当たり、導電性のケーシングパイプを有するボーリング孔に電極を挿入し所定の深度で通電する。通電によって励起された電位分布を地表や別のボーリング孔内にセットした電極で計測する。通電深度を変えて同様の測定を繰り返し、多数の電流源―電位測定点の組み合わせのデータを取得し、これらのデータから逆解析プログラムを使って比抵抗分布を求める。
【００１５】
この場合にケーシングパイプへの通電は、電流源深度を変更しつつ行うことにすれば測定精度が向上し、好ましい。
【００１６】
解析のアルゴリズムは一般の比抵抗逆解析と同様の手順に従う。まず初期比抵抗モデルを与え、それに基づいて測定データと同じ電極配置でシミュレーションを行う。シミュレーションで求められた電位（抵抗又は見掛け比抵抗の場合もある）と測定された電位（抵抗又は見掛け比抵抗の場合もある）の食い違いが小さくなるように初期比抵抗モデルを修正する。修正された比抵抗モデルを初期モデルとして同様の修正を繰り返し、シミュレーション結果と測定結果が十分一致するような比抵抗モデルを求め、これを最終結果とする。
【００１７】
従来の比抵抗トモグラフィ解析との相違点は、シミュレーションにおいて導電性のケーシングパイプの影響を適切に評価することである。通常の比抵抗トモグラフィの解析プログラムでは電流源を点ソースとしてシミュレーションが行われるが、ケーシングパイプ中に電流源が存在する場合には、ケーシングパイプに沿った線状区間で連続的に電流が流れ出す（線電流源）。その効果を適切に反映させて計算を行うためには２つの方法が考えられる。
【００１８】
一つは、有限要素法（ＦＥＭ）を使った電位シミュレーションにおいて、ケーシングパイプをこれと等価な線要素及び接地抵抗要素（ジョイント要素とも云う）によって近似して、解析モデル中に組み込む方法である。ケーシングパイプをＦＥＭの線要素及び接地抵抗要素を使ってモデル化してケーシングパイプの影響を適切に評価し、ケーシングパイプ上の電流密度を評価することによって、比抵抗の解析精度を向上することができる。
【００１９】
もう一つは、ケーシングパイプ中に電流源を設置したときの電流分布を別途計測したのち、これを基に線電流源を多数の点ソース又は微小線ソースの重ね合わせとして取り扱う方法である。ケーシングパイプ電流分布測定装置を用いて計測した電流密度分布を使い、測定されたケーシングパイプから流出する電流を点電流源又は微小線電流源の重ね合わせとして計算し、ＦＥＭ又は差分法によって解析することができる。
【００２０】
後者の方法による場合には、既往の差分法又は有限要素法による比抵抗トモグラフィ解析プログラムを若干改良するだけで容易に対応することができる。
【発明の効果】
【００２１】
本発明によれば、バックグランド比抵抗が小さく、かつボーリングの深度が数１００ｍ以上の場合には、ボーリング孔に導電性のケーシングパイプが挿入されていても、ボーリング孔内に測定点を設けた比抵抗トモグラフィを実施することができる。これによって、ケーシングパイプの１ヶ所に電流源を設けて測定を行う通常の流電電位法に比べると深度方向の比抵抗分布の情報がより多く得られるとともに、ケーシングパイプ先端付近の地下深部の探査精度も大幅に改善される。
【発明を実施するための最良の形態】
【００２２】
以下図面を参照して本発明の実施の形態を説明する。
【００２３】
図６は、本発明に係る比抵抗測定方法の全体を概念的に示す説明図であって、本発明を実施する対象の地層内を透視した模型図に、本発明の実施に用いる解析装置を配置したものを示した。図６の配置に加えて複数のボーリング孔があってもかまわない。また、この他に地表の測定点を直線上に並べた２次元探査の配置も可能である。
【００２４】
図６に示すように、地層の比抵抗解析対象区域１００に地表からケーシングパイプ１０が沈設され、その頭部又は内部に電極が設置される。トランスミッタ１１は電線ケーブルを介して電極に電圧を負荷する。
【００２５】
地表１０１に、所望の配列で多数の各測点が１０２が配設されており、これらの測点の２点間又はこれらの測点と遠方の基準電極（遠電極）間の２点間のそれぞれの電位差の測定を行う。この２点間の電位差測定値はレコーダ２０に記録される。この記録はコンピュータ２１に送られ、コンピュータ２１はこのデータを解析して地層の比抵抗を算出する。
【００２６】
先に言及したように、ケーシングパイプ１０内に電流源がある場合の比抵抗トモグラフィの解析では、一般にケーシングパイプ１０から流出する電流１２の密度がケーシングパイプ１０の長手方向に沿って一定であるものとして近似計算が行われていた（流電電位法）。しかし、ケーシングパイプ１０から流出する電極上の実際の電流１２の密度分布はケーシングパイプと地山との間の接地抵抗及びケーシングパイプ１０自身の電気抵抗によって変化する。
【００２７】
特に比抵抗解析対象区域１００の地山の比抵抗（バックグランド比抵抗）が低い場合、流電電位探査方法の解析におけるケーシングパイプ１０の抵抗が近似精度に与える影響を評価することが重要な課題であり、このためには、何らかの手段によってケーシングパイプ１０の長手方向に沿って流出する電流分布を測定することが必要である。
【００２８】
また、計測された電流分布を使って流電電位探査方法の解析を行うことができれば、解析の精度はこれまでより以上に向上することを期待することができる。
【００２９】
ケーシングパイプ１０から流出する電流密度に影響を与える第１の要因であるケーシングパイプ１０と地山との間の接地抵抗は、ケーシングパイプ１０と地山とのクリアランス、ケーシングパイプ１０の周囲の腐食、ケーシングパイプ１０の周りに充填したセメントの影響、ケーシングパイプ１０の径、あるいは周辺地山の比抵抗というような様々な要因と関連している。
【００３０】
また、第２の要因であるケーシングパイプ１０自身の電気抵抗は、バックグランド比抵抗が１Ωｍ以下の極端に低い区域において、特に大きな影響を与えることがシミュレーションによって明らかになっている。
【００３１】
ＦＥＭ（有限要素法）による解析プログラムではケーシングパイプ１０と地山との間の接地抵抗及びケーシングパイプ１０自身の電気抵抗を線要素と接地抵抗要素で近似することによって、これらの影響をある程度考慮した電位シミュレーションを行うことができる。しかし、そのためにはケーシングパイプ１０と地山との間の接地抵抗及びケーシングパイプ１０自身の電気抵抗を予め精度良く求めておく必要がある。
【００３２】
図５に電位シミュレーションの例を模式的に示した。例えば、長さＬ＝１０００ｍ、直径２０ｃｍ、肉厚６ｍｍ、導電率１０６Ｓ／ｍのケーシングパイプ１０を対象区域１００内で鉛直方向（深度ｚ方向）に沈下し、ケーシングパイプ１０に直交する地表面１０１を設定し、この地表面１０１に距離ｐをとって多数の測定点を設定する。地表面１０１とケーシングパイプ１０との交点１０３を座標（０、０）とすると、任意の測点１０２は座標（ｐ、ｚ）で表わされる。
【００３３】
電位シミュレーションは３次元ＦＥＭによって実施することとし、図１〜図４に示すような解析用３次元ＦＥＭモデル１を作成した。実際の適用で２次源測線を設定した場合には、３次元ＦＥＭの代わりに２．５次元ＦＥＭを使ってシミュレーションすることも可能である。
【００３４】
図１は比抵抗を測定すべき対象区域のＦＥＭモデル１の斜視図、図２はその平面図、図３はその中心位置の立面断面図、図４は図３の１ａ部の部分拡大図である。図１に示すように、対象区域はケーシングパイプ１０を中心とし、Ｘ，Ｙ方向がそれぞれ１００，０００ｍ、深度５０，０００ｍの直方体状の区域であり、ケーシングパイプ１０が、図４に示すように、深度１０００ｍまで打ち込まれた場合のＦＥＭモデル１である。
【００３５】
ＦＥＭの計算では、半無限媒質を四面体要素で分割するとともに、ケーシングパイプ１０については線要素（長さ５ｍ）および接地抵抗を表す接地抵抗要素を組み合わせてモデルを作成した。線要素の物性値パラメータに関しては、ケーシングパイプ１０の断面積と比抵抗から単位長さあたりの抵抗を算出した。線要素は太さを表すパラメータを本来持たないが、離散化効果によって周辺の三次元領域の要素の形状やサイズに従った等価な電極径に対する電位応答を示すことが知られているのでこれを利用した。
【００３６】
図７に、３次元ＦＥＭ（四面体要素を使った場合）の要素サイズおよび要素形状と等価電極径（直径）との関係を示した。図７は、ＦＥＭモデルの四面体要素のサイズ（ｍ）を横軸に取り、縦軸に等価電極径（ｃｍ）を取って、その関係をプロットしたものである。解析結果次式の関係が得られている。
【００３７】
Ｄ＝０．３６１１×ｘ１．０２４３        ……（１）
ここで、  Ｄ：等価電極径（ｃｍ）
ｘ：線要素サイズ（ｍ）
である。
【００３８】
図１に示した電位シミュレーションのためのＦＥＭモデル１では要素サイズが５ｍであり要素形状は正方形であるから、図７に示す経験式からその等価電極径は４３．８ｃｍと評価される。そのままでは電極径４３．８ｃｍの電極に対する電位分布が計算されることになる。そこで直径２０ｃｍの電極に対する電位シミュレーションを行うため、目標とする電極径と等価電極径の差を接地抵抗に換算して、これに相当する接地抵抗要素を、線要素と３次元要素の間に挿入した。
【００３９】
次に、ケーシングパイプ長１０００ｍの場合の電位シミュレーション結果について説明する。
【００４０】
バックグランド比抵抗を１Ωｍ、１０Ωｍ、１００Ωｍ、１０００Ωｍとし、それぞれについて１Ａ通電したときの計算電位の分布を図９〜図１２に示す。図９〜図１２は等電位線によって分布を示し、等電位線に電位の値を併記した。
【００４１】
図９に示すように、バックグランド比抵抗が１Ωｍの場合、測定範囲内（地表部）における電位が概ね０．１ｍＶ〜８ｍＶであり、測定電位レベルとしてはやや小さいが、電流量を５Ａ〜１０Ａ程度まで増やすことによって十分に測定可能なレベルとなる範囲である。また、図１０〜図１２に示すように、バックグランド比抵抗がそれぞれ１０Ωｍ、１００Ωｍ、１０００Ωｍと大きくなるにつれて測定電位も大きくなり、十分に測定可能な電位レベルが得られることが確認された。
【００４２】
ＦＥＭによって定式化された電流と電位との関係式を使うことによって、解析された電位から要素間でやりとりされる電流を求めることができる。
【００４３】
図８は、このようにして線電極上の電位から求めた、ケーシングパイプとバックグランド領域の間でやりとりされる電流密度分布を示したものである。図８に示す４本の電流密度曲線３１、３２、３３、３４はそれぞれ図９〜図１２に対応するもので、バックグランド比抵抗が１Ωｍ、１０Ωｍ、１００Ωｍ、１０００Ωｍの場合を示している。
【００４４】
バックグランド比抵抗が１０００Ωｍの場合（曲線３４）では、ケーシングパイプ上の電流密度は深度１０００ｍに亘ってほぼ一定となり、理想的な線電極周辺の電位分布を示すが、曲線３３、３２、３１に示すように、バックグランド比抵抗が小さくなるにつれてケーシングパイプ１０自身の電気抵抗を無視することができなくなり、ケーシングパイプ先端まで電流が到達しなくなることが分かる。
【００４５】
図９に示す電位分布でも、バックグランド比抵抗が１Ωｍのケースでは、通常の点電流源に対する電位分布とそれほど変わらない電位分布を示すようになる。
【００４６】
このように、バックグランド比抵抗が低い場合には、ケーシングパイプ１０に通電しても従来考えられてきたように線電極の電位応答を示すわけではない。そのため、流電電位法をそのまま適用しても深部の情報を取得できないということも起こり得る。これは、実際のケーシングパイプが理想的な導電体でないため、ケーシングパイプ自身の電気抵抗の影響が現れるためである。しかし逆に考えると、このことは電流源の位置情報がケーシングパイプによって完全にマスクされないということを意味しており、電流源深度を変えることによって異なる配置のデータを取得することができるので、ケーシングパイプで保孔されたボーリング孔を使った比抵抗トモグラフィの可能性を示すものでもある。
【００４７】
次に、流体流動電位法をＥＯＲに適用することを目的として流動電位法におけるケーシングパイプの抵抗の影響を検討するため、塩水で飽和した低比抵抗の堆積岩地帯に流電電位法を適用したケースを想定してシミュレーション・モデルを作成し、シミュレーション検討を行った。
【００４８】
シミュレーションに用いた比抵抗モデルを図１３、図１４に示す。図１３は平面図でケーシングパイプ１０を中心とするＸ，Ｙ方向１０００ｍの解析対象区域１００を示している。図１４はそのケーシングパイプ１０の位置を通る立面断面図である。ケーシングパイプ１０が深度５００ｍまで挿入されており、その下端から１０ｍ離れた側方に比抵抗１Ωｍの厚さ４０ｍの扇形平板状の比抵抗異常体４０が設けられている。図１３に示すように、その平板状の比抵抗異常体４０の平面形は中心角４５度の扇形である。この扇形内に比抵抗測定部が配置されている。この解析対象区域１００のバックグランド比抵抗は１０Ωｍである。
【００４９】
このシミュレーション・モデルは比抵抗１０Ωｍの媒質中にある深度５００ｍのケーシングパイプに通電することを想定しており、比抵抗異常体４０を検出しようとするものである。
【００５０】
ケーシングパイプ１０は直径２０ｃｍ、肉厚６ｍｍ、比抵抗１０−６Ωｍとした。
【００５１】
図１３に示すように、ケーシングパイプ１０を中心として１２方向に展開した長さ５００ｍの測線４１上に２０ｍ間隔で電位測定点４２（図中に●で示す電位電極）を設定した。また、電流源４３（図中に○で示す電流源）の位置については以下の（１）〜（４）の４ケースを想定してシミュレーション・データを作成し、それぞれを解析した。
【００５２】
また、図１４に示すようにケーシングパイプ１０内に深度５０ｍごとにケーシングパイプ内電流源４４を配置した。
（１）流電電位法（３００データ）
ケーシングパイプの上端部及び下端部に１点通電した。これは通常の流電電位法測定に相当する。
（２）流電電位法に地表電流源を追加（１８，２４０データ）
流電電位法測定に加えて、地表測線に１００ｍ間隔で電流源を設けて測定することとした。
（３）流電電位法にケーシングパイプ内電流源を追加（３，３００データ）
上記（１）の流電電位法に加えて、ケーシングパイプ内に５０ｍ問隔で電流源を設けて測定することとした。
（４）流電電位法に地表電流源及びケーシングパイプ内電流源を追加（２１，２４０データ）
上記の（２）と（３）を併用した測定を行うこととした。
【００５３】
ここで検討したモデルでは、バックグランド比抵抗が１０Ωｍと低いのでケーシングパイプ自身の電気抵抗の影響が現れて完全な導体の挙動は示さない。
【００５４】
図１５、図１６はそれぞれ、電流源４３の深度が異なる場合の電位分布の変化を示す立面図で、図１５は電流源４３の深度が深度０ｍ、図１６は電流源４４が深度５００ｍの位置にあるときの電位分布を示している。電流源４３，４４がケーシングパイプ１０内にあっても電流源４３，４４の深度によってケーシングパイプ１０外の電位分布は若干変化する。
【００５５】
この情報をトモグラフィ・データとして活用することによって、深部のトモグラフィの解析精度を向上することができる。このような測定配置は、地震探査におけるオフセットＶＳＰ（Ｖｅｒｔｉｃａｌ  Ｓｅｉｓｍｉｃ  Ｐｒｏｆｉｌｉｎｇ）や坑井と地表間のトモグラフィに相当するものである。
【００５６】
次に、図１７〜図１８に示すように、ケーシングパイプを中心として同心円筒形にグリッド５１を配置した。図１７は水平断面図、図１８は図１７の中心を通る立面断面図（鉛直断面図）である。各グリッド点５２に比抵抗を与えることとし、各要素の比抵抗はグリッド点５２の比抵抗から内挿されるものとして定義することとした。このように各要素５１の比抵抗を定義すると、要素５１の比抵抗とグリッド点５２の比抵抗は線形関係となる。従って、各要素５１の感度に内挿の際の重み係数（内挿関数の係数）をかけることによって、容易にグリッドの感度を導くことができる。このようにして求められた各グリッドの感度を使って制約付最小二乗法による逐次修正法を適用し、最適な比抵抗モデルを求めた。なお、初期モデルは均一な比抵抗分布であるとし、比抵抗の値は全データの見掛比抵抗のメディアン値を使った。
【００５７】
逆解析した結果を図１９〜図２２に示した。図１９〜図２２はケーシングパイプ１０を含む縦断面図で、測定された比抵抗分布を等高線で描いたもので、等高線に付されている数字は比抵抗値（Ωｍ単位）を示したものである。
【００５８】
図１９は通常の流電電位法による測定に係るもので、ケーシングパイプ１０の孔口に電流源４３を設けて通電した通常の流電電位法の測定結果である。深度５００ｍに存在する低比抵抗異常体４０は全く検出されずに、その上方の深度１００ｍ付近にバックグランド比抵抗の値よりやや低い比抵抗の擬像６１が現れており、低比抵抗異常体４０について、検出精度及び深度方向位置決定精度が低いことが確認された。
【００５９】
これに対して図２０に示すように、地表１０１に多数の電流源４３を設けて測定データを追加すると、精度はやや改善され、深度５００ｍのケーシングパイプ１０の先端付近に右上に伸びる低比抵抗部異常体４０の像６２が検出された。このように、通常の流電電位法にケーシングパイプ１０以外の地表部に電流源４３を迫加することによって、探査精度が向上することが示された。
【００６０】
図２１は、図２０にケーシングパイプ１０の深度方向に沿って５０ｍごとにケーシングパイプ内電流源４４を追加して測定した結果を示す断面図で、ケーシングパイプ１００の先端付近の低比抵抗異常体４０の直上にやや上向きに拡がる像６３が明瞭に検出されている。
【００６１】
さらに、地表１０１及びケーシングパイプ１０内の双方に電流源４３及び４４を追加した結果、図２２に示すように、低比抵抗異常体４０の近傍に、その明確な像６４が現出した。従って、探査精度の大幅な改善が認められた。
【００６２】
ケーシングパイプ１０自身の電気抵抗を無視することができる場合は、坑井内のどの位置に電流源があってもケーシングパイプ外の電位分布は変わらないため、ケーシングパイプ１０内で電流源の位置を変えた測定を実施しても探査精度の向上は望めない。しかし、バックグランド比抵抗が小さな場合には、ケーシングパイプ１０内の各深度の位置に電流源を配設したトモグラフィ的な測定によって、探査精度が大きく向上することが分かった。
【００６３】
以上のように、塩水で飽和した低比抵抗の堆積岩地帯に流電電位法を適用した場合の問題点をミュレーションによって検討した結果、低比抵抗領域で流電電位法を適用した場合には、ケーシングパイプ長の増大とともにケーシングパイプの抵抗の影響によって電流が先端まで流れなくなり、その場合には流電電位法を適用してもケーシングパイプ先端付近の情報が得られないという可能性があったが、ケーシングパイプ内部の複数の深度に電流源をおいて逆ＶＳＰ型の測定を行うことも原理的には可能であり、これによって探査精度を大幅に改善することができた。
【図面の簡単な説明】
【００６４】
【図１】本発明の解析用ＦＥＭモデルの斜視図である。
【図２】図１の平面図である。
【図３】図１の中心位置の立面断面図である。
【図４】図３の部分拡大図である。
【図５】電位シミュレーションを示す模式図である。
【図６】本発明方法の全体説明図である。
【図７】ＦＥＭ要素サイズ及び要素形状と等価電極径との関係を示すグラフである。
【図８】ケーシングパイプの長手方向に沿う流出電流密度分布を示すグラフである。
【図９】バックグランド比抵抗が１Ωｍのときの地中電位分布を示す図である。
【図１０】バックグランド比抵抗が１０Ωｍのときの地中電位分布を示す図である。
【図１１】バックグランド比抵抗が１００Ωｍのときの地中電位分布を示す図である。
【図１２】バックグランド比抵抗が１０００Ωｍのときの地中電位分布を示す図である。
【図１３】シミュレーションに用いた比抵抗モデルの平面図である。
【図１４】図１３の中心を通る立面断面図である。
【図１５】電流源深度が０ｍの場合の電位分布を示す断面図である。
【図１６】電流源深度が５００ｍの場合の電位分布を示す断面図である。
【図１７】同心円筒形にグリッドを配したモデルの平面図である。
【図１８】図１７の中心を通る立面断面図である。
【図１９】逆解析結果を示す低比抵抗部の検出を示す立面断面図である。
【図２０】逆解析結果を示す低比抵抗部の検出を示す立面断面図である。
【図２１】逆解析結果を示す低比抵抗部の検出を示す立面断面図である。
【図２２】逆解析結果を示す低比抵抗部の検出を示す立面断面図である。
【符号の説明】
【００６５】
１    ＦＥＭモデル
１０    ケーシングパイプ
１１    トランスミッタ
１２    電流
１３    流体流
２０    レコーダ
２１    コンピュータ
３１、３２、３３、３４    電流密度曲線
４０    比抵抗異常体
４１    測線
４２    電位測定点
４３、４４    電流源
５１    グリッド（要素）
５２    グリッド点
６１    擬像
６２、６３、６４    像
１００    解析対象区域
１０１    地表
１０２    測点
１０３    交点
【図１】
2007285729.tif 000003
【図２】
2007285729.tif 000004
【図３】
2007285729.tif 000005
【図４】
2007285729.tif 000006
【図５】
2007285729.tif 000007
【図６】
2007285729.tif 000008
【図７】
2007285729.tif 000009
【図８】
2007285729.tif 000010
【図９】
2007285729.tif 000011
【図１０】
2007285729.tif 000012
【図１１】
2007285729.tif 000013
【図１２】
2007285729.tif 000014
【図１３】
2007285729.tif 000015
【図１４】
2007285729.tif 000016
【図１５】
2007285729.tif 000017
【図１６】
2007285729.tif 000018
【図１７】
2007285729.tif 000019
【図１８】
2007285729.tif 000020
【図１９】
2007285729.tif 000021
【図２０】
2007285729.tif 000022
【図２１】
2007285729.tif 000023
【図２２】
2007285729.tif 000024