(57)【要約】【課題】不均質な透水係数の分布を考慮した地下水流動解析を行う際の解析データ作成作業の簡易化を図り、効率的に解析を行うことのできる地下水流動解析方法及びそのプログラムを提供する。【解決手段】二重相反境界要素法を用いて、解析領域における透水係数を与えて境界積分方程式を解くことにより解析領域内の地下水の全水頭を求める方法において、解析領域における点ｘm（ｍ＝１〜Ｍ）に透水係数Ｋ（ｘm）を与え、任意の点ｘの透水係数Ｋ（ｘ）の分布を、ＲＢＦ又はＣＳ−ＲＢＦｓを用いた近似関数ｆ（ｘ，ｘm）を用いた式で与え、この透水係数の分布を表すＫ（ｘ）と、その微分形である∂Ｋ（ｘ）／∂ｘを、境界積分方程式に与える。【選択図】        図１
2010014619.tif 000002
【特許請求の範囲】
【請求項１】
二重相反境界要素法を用いて、解析領域における透水係数を与えて境界積分方程式を解くことにより前記解析領域内の地下水の全水頭を求める地下水流動解析方法において、
前記解析領域における点ｘm（ｍ＝１〜Ｍ）に透水係数Ｋ（ｘm）を与え、前記解析領域における任意の点ｘの透水係数Ｋ（ｘ）の分布を、近似関数ｆ（ｘ，ｘm）を用いて次式
【数１】
2010014619.tif 000042
で与え、
前記近似関数ｆ（ｘ，ｘm）として、ＲＢＦ又はＣＳ−ＲＢＦｓを適用し、
前記透水係数Ｋ（ｘ）の微分形を、次式
【数２】
2010014619.tif 000043
で与えることを特徴とする地下水流動解析方法。
【請求項２】
前記近似関数ｆ（ｘ，ｘm）として、次式
【数３】
2010014619.tif 000044
を適用することを特徴とする請求項１に記載の地下水流動解析方法。
【請求項３】
請求項１又は２に記載の地下水流動解析方法を、コンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、二重相反境界要素法を用いた地下水流動解析方法及びそのプログラムに関するものである。
【背景技術】
【０００２】
原子力発電所等から発生する高レベル放射性廃棄物の処分方法として、図１０に示すように、地下３００ｍ以深にトンネル群を建設してその中に放射性廃棄物を埋設する「地層処分」が検討されている。地層処分では、溶解したガラスに放射性廃棄物を溶け込ませてガラス固化体とした後に、これをオーバーパックと称される炭素鋼製の容器に収容し、処分トンネル内に埋設する。
【０００３】
この地層処分の安全性を評価する項目（シナリオ）の一つとして、地下水シナリオと呼ばれるものがある。この検討では、ガラス固化体から溶出した放射線核種が、オーバーパック周辺のベントナイトからなる人工バリアを通って、周辺岩盤からなる天然バリアの地下水に溶出し、地表近傍の生物圏へ輸送されて人間が被爆するまでのプロセスが対象となる。そのため、想定される放射性廃棄物の埋設箇所を含む数十キロの領域を対象とした地下水の流動状況を把握することが必要となる。そこで、近年、有限要素法（ＦＥＭ）や有限差分法（ＦＤＭ）に代表される数値解析手法を用いた地下水流動解析（浸透流解析）が行われている。この地下水流動解析では、地層の物性のひとつである透水係数をもとに、解析領域内の地下水の全水頭及び流速ベクトル等を求める（例えば特許文献１を参照）。
【０００４】
ところで、わが国は環太平洋の地質擾乱帯に位置することから、広域の地下水流動解析モデルの地質構成は複雑になる。このため、現状の有限要素法等を用いた解析では、解析モデルはかなり複雑なものとなることが想定される。
【０００５】
さらに、同じ地層内であっても、実際には透水係数は一定ではなくばらつきがある。すなわち、透水係数の分布が不均質である。有限要素法等の領域離散型の解析手法を用いて、透水係数の不均質性を考慮したモデルを作成する場合、解析領域を有限要素に分割して個々の要素に透水係数を与える必要があるため、解析モデルの作成作業に労力を要するという問題がある。
【０００６】
一方、解析対象領域内部の要素分割を必要としない境界要素法（ＢＥＭ）を用いた解析手法も検討されている。境界要素法は、一般には、解析領域の物性が均質な場合の基本解を用いている。このため、解析領域の物性（透水係数）が不均質な場合には、境界積分方程式への定式化において、有限要素法や有限差分法と同様に離散化した内部領域に対する領域積分を行う必要がある。領域積分をするためには解析領域を要素分割する必要があるため、結局、上述した有限要素法と同様の労力を要することとなる。
【０００７】
これに対し、二重相反境界要素法（DRBEM:Dual Reciprocity Boundary Element Method）を用いて、上記の領域積分を境界積分に置換する研究が報告されている。すなわち、この手法は、領域積分を境界積分に置換することで、解析対象領域に対して内部要素分割を不要としたものである。そこで、本発明者は、非特許文献１において、この解析領域の要素分割を必要としない二重相反境界要素法を、上述した地下水の流動解析に活用することを検討した。
【０００８】
【特許文献１】特開２００６−２４２８９３号公報
【非特許文献１】第３７回岩盤力学に関するシンポジウム講演集、（社）土木学会、２００８年１月１０日〜１１日、講演番号４１、pp.233-238
【発明の開示】
【発明が解決しようとする課題】
【０００９】
上述した二重相反境界要素法を適用した境界積分方程式の定式化においては、解析領域内の透水係数の不均質性を、透水係数Ｋ（ｘ）の関数およびその微分形∂Ｋ（ｘ）／∂ｘで表す方法が用いられ、これらの関数を境界積分方程式に与えて計算を行っている。この場合、解析領域の透水係数Ｋ（ｘ）の分布が三角関数や多項式のように規則的な分布の場合には有効である。
【００１０】
しかしながら、一般的には地層中の透水係数分布は複雑であるため、三角関数を用いたフーリエ級数などで複雑な分布を表現することには限界がある。さらに、この方法では、対象とする地層ごとに、透水係数Ｋ（ｘ）の関数とその微分した関数の両方をプログラム中に組み込む必要があるため、効率が悪いという問題がある。
【００１１】
本発明は、上記の点に鑑み、不均質な透水係数の分布を考慮した地下水流動解析を行う際の解析データ作成作業の簡易化を図り、効率的に解析を行うことのできる地下水流動解析方法及びそのプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１２】
本発明の請求項１に係る地下水流動解析方法は、二重相反境界要素法を用いて、解析領域における透水係数を与えて境界積分方程式を解くことにより前記解析領域内の地下水の全水頭を求める地下水流動解析方法において、前記解析領域における点ｘm（ｍ＝１〜Ｍ）に透水係数Ｋ（ｘm）を与え、前記解析領域における任意の点ｘの透水係数Ｋ（ｘ）の分布を、近似関数ｆ（ｘ，ｘm）を用いて次式
【数４】
2010014619.tif 000003
で与え、前記近似関数ｆ（ｘ，ｘm）としてＲＢＦ又はＣＳ−ＲＢＦを適用し、前記透水係数Ｋ（ｘ）の微分形を、次式
【数５】
2010014619.tif 000004
で与えることを特徴とする。なお、これらの定数γmは、点ｘm（ｍ＝１〜Ｍ）のＫ（ｘm）が既知であることより、ベクトル表示で｛γ｝＝［ｆ]-1{Ｋ}で与えられる。
【００１３】
また、本発明の請求項２に係る地下水流動解析方法は、上記請求項１における前記近似関数ｆ（ｘ，ｘm）として、次式
【数６】
2010014619.tif 000005
を適用することを特徴とする。
【００１４】
また、本発明の請求項３に係るプログラムは、請求項１又は２に記載の地下水流動解析方法をコンピュータに実行させるものである。
【発明の効果】
【００１５】
本発明の地下水流動解析方法及びそのプログラムによれば、二重相反境界要素法を用いて、不均質な透水係数の分布を考慮した解析を行う際に、解析領域内の特定の点ｘmに透水係数のデータＫ（ｘm）を与えるだけで解析を行うことが可能になるため、従来のように透水係数の複雑な分布をフーリエ級数等の関数としてプログラムする必要がなくなる。その結果、解析データの作成作業の簡易化を図ることができ、効率的に解析を行うことができるようになる。
【発明を実施するための最良の形態】
【００１６】
以下に、添付図面を参照して、本発明に係る地下水流動解析方法及びそのプログラムの好適な実施の形態について詳細に説明する。
【００１７】
まず、二重相反境界要素法（DRBEM:Dual Reciprocity Boundary Element Method）を用いた地下水流動解析の概要について説明する。なお、本実施の形態で対象となる地下水流動解析とは、地下水の時間的変動を考慮しない定常浸透流解析である。
【００１８】
[境界積分方程式の導出]
境界要素法は、一般には、解析領域の物性が均質な場合の基本解を用いている。このため、解析領域内の物性の不均質性を考慮する場合、境界積分方程式への定式化において、有限要素法（ＦＥＭ）や有限差分法（ＦＤＭ）と同様に離散化した内部領域に対する領域積分を行う必要があった。これに対し、Nadiniらは、二重相反境界要素法を用いてこの領域積分を境界積分に置換する研究を報告している。本実施の形態では、この手法を活用して、透水係数の分布が不均質な解析領域の地下水流動解析を行う。以下、境界積分方程式の導出について説明する。
【００１９】
地下浸透流問題において透水係数の空間変動（不均質さ）を考慮した支配方程式は次式で表される。
【数７】
2010014619.tif 000006
ここに、Ｋ＝Ｋ（ｘ）は透水係数、φ＝φ（ｘ）は全水頭（ポテンシャル）、∇はナブラである。
【００２０】
（数７）式を展開して整理すると次式が得られる。
【数８】
2010014619.tif 000007
【００２１】
上式の右辺をｂ（ｘ）とおけば、次式のポアソン方程式が得られる。
【数９】
2010014619.tif 000008
【数１０】
2010014619.tif 000009
このポアソン方程式を満たす特解をφ′とすると、（数９）式は次式と表すことができる。
【数１１】
2010014619.tif 000010
すなわち次式となる。
【数１２】
2010014619.tif 000011
よって、φnew＝φ−φ′と置くことで、φnewに対する新たなラプラス方程式が得られる。従って、従来の境界要素法と同様の定式化により境界積分方程式を導出し、解くことができる。
【００２２】
しかしながら、一般的な透水係数分布で与えられる関数ｂ（ｘ）に対して、この特解φ′を具体的に求めることは難しい。このため、本手法では、まず任意の点ｘのｂ（ｘ）を、解析対象領域を離散化した境界節点（１〜Ｎ）と領域内節点（１〜Ｌ）の各節点ｘjにおいて、係数αjとスカラーの近似関数ｆ（ｘ，ｘj）を用いた線形和として次式で表す。
【数１３】
2010014619.tif 000012
【００２３】
ここで、上式の近似関数ｆ（ｘ，ｘj）の特解を利用することを考える。すなわち、（数１１）式を満たす特解φ′も（数１３）式と同様に係数αjと非同次項を近似関数ｆ（ｘ，ｘj）とする特解ψj＝ψ（ｘ，ｘj）の線形和で表すことができる。ここで、ψ（ｘ，ｘj）は非同次項をｆ（ｘ，ｘj）とするポアソン方程式の特解である。
【数１４】
2010014619.tif 000013
【００２４】
上式より次式が得られる。
【数１５】
2010014619.tif 000014
【００２５】
よって、（数９）式を（数１２）式と表現したのと同様に次式を得ることができる。
【数１６】
2010014619.tif 000015
【００２６】
上式から、境界積分方程式は次式と表される。
【数１７】
2010014619.tif 000016
ここに、ciは、φの境界要素及び内部節点の幾何学的条件により定められる定数、Γは領域境界、ｎは境界の法線方向、φ*はラプラス方程式の基本解であり、二次元の基本解は次式で与えられる。
【数１８】
2010014619.tif 000017
【００２７】
（数１７）式の各項をベクトルとマトリクスで表現すると次式のようになる。
【数１９】
2010014619.tif 000018
ここに、φ，α，ｑは各節点ｉのφiとαi及び境界点ｉの∂φ／∂ｎによるベクトルであり、マトリクスＨ，Ｇ，Ψ，Ｑの各成分は次式のようになる。
【数２０】
2010014619.tif 000019
【数２１】
2010014619.tif 000020
【数２２】
2010014619.tif 000021
【数２３】
2010014619.tif 000022
ここに、δijは、クロネッカのデルタである。
【００２８】
（数１９）式の係数ベクトルαは、（数１３）式をマトリックスとベクトルで表示して用いれば、次式のように表される。
【数２４】
2010014619.tif 000023
【００２９】
ベクトルｂは領域内のＫ（ｘ）とφ（ｘ）で表される。そこで、φ（ｘ）を（数１３）式と同様に各点の係数βjと近似関数ｆ（ｘ，ｘj）を用いた線形和として次式で表現する。
【数２５】
2010014619.tif 000024
【００３０】
従って、（数１０）式に代入し次式を得る。
【数２６】
2010014619.tif 000025
【００３１】
上式をベクトルとマトリクスで表示すると次式となる。
【数２７】
2010014619.tif 000026
ここに、ＦKは次式のマトリクスとなる。
【数２８】
2010014619.tif 000027
【００３２】
従って、（数２４）式は、次式に書き直すことができる。
【数２９】
2010014619.tif 000028
【００３３】
（数１９）式に代入し次式が得られる。
【数３０】
2010014619.tif 000029
上式は、境界要素と内部節点のφ及び境界要素の∂φ／∂ｎによってつくられる連立方程式となる。この連立方程式を解くことにより未知境界量のφ、∂φ／∂ｎを求めることができる。
【００３４】
[近似関数ｆ（ｘ，ｘj）と特解ψ]
ところで、上述したｂ（ｘ）の近似関数ｆ（ｘ，ｘj）を一意的に定めることは難しく、従来、ｒの多項式やlog関数などが用いられてきたが、近年、ChenらによりWendlandのＲＢＦ（radial basis functions）やＣＳ−ＲＢＦｓ（compactly suppo-rted radial basis functions）を用いることが提案されている。そこで、本実施の形態においても、Ｃ２クラスの関数として以下のＲＢＦを用いている。
【数３１】
2010014619.tif 000030
ここに、ｒはｒ＝｜ｘ−ｘj｜，ａは定数である。
【００３５】
（数３１）式より（数１４）式の特解ψ，∂ψ／∂ｎは、次式で与えられる。
【数３２】
2010014619.tif 000031
【数３３】
2010014619.tif 000032
【００３６】
[透水係数Ｋ（ｘ）の表現方法]
上述した二重相反境界要素法を用いた境界積分方程式の定式化においては、（数２６）式に示されるように境界節点と内部節点において透水係数Ｋ（ｘ）及びその微分∂Ｋ（ｘ）／∂ｘを与える必要がある。上述した非特許文献１では、透水係数Ｋ（ｘ）とその微分形を直接関数として与えて、（数２８）式のマトリクスＦKを計算している。
【００３７】
しかしながら、実際の地層内の透水係数Ｋ（ｘ）の分布は複雑であるため、三角関数を用いたフーリエ級数などで複雑な分布を表現することには限界があると考えられる。また、対象とする地層ごとに関数とその微分形を作成してプログラムの中に導入するのは効率的ではない。
【００３８】
そこで、本発明者が鋭意検討した結果、透水係数Ｋ（ｘ）とその微分形∂Ｋ（ｘ）／∂ｘを以下のように表現することで、解析データの作成作業の簡易化を図り、解析効率を向上させることが可能になるとの知見を得た。
【００３９】
図２は、本実施の形態で適用する透水係数Ｋ（ｘ）の空間分布を模式的に示した図である。解析領域（解析モデル）において、白抜きの丸は境界節点（１〜Ｎ）であり、白抜きの四角は領域内節点（１〜Ｌ）である。図２に示すように解析領域を覆う透水係数Ｋの分布を考える。そして、解析領域におけるＮ個の境界節点やＬ個の内部節点とは必ずしも一致しない複数の点ｘm（ｍ＝１〜Ｍ：図２に示す白抜きの三角）に、それぞれ透水係数Ｋ（ｘm）を与える。この点ｘmに与える透水係数Ｋ（ｘm）は、解析対象となる実際の地層を測定して得たデータ、又は、これらの測定値や地質情報から推定されるデータである。これを用いて解析領域内の透水係数分布を（数１３）式又は（数２５）式と同様に次式で表す。
【数３４】
2010014619.tif 000033
なお、これらの定数γmは、点ｘm（ｍ＝１〜Ｍ）のＫ（ｘm）が既知であることより、ベクトル表示で｛γ｝＝［ｆ]-1{Ｋ}で与えられる。
【００４０】
（数３４）式において、ｆ（ｘ，ｘm）は近似関数である。すなわち、（数３４式）は、解析領域における任意の点ｘの透水係数Ｋ（ｘ）の近似分布を示す式である。近似関数ｆ（ｘ，ｘm）としては様々な関数が考えられるが、上述したｂ（ｘ）の近似関数ｆ（ｘ，ｘj）と同様に、ＲＢＦ（radial basis functions）又はＣＳ−ＲＢＦｓ（compactly supported radial basis functions）を用いる。本実施の形態では、上述したｂ（ｘ）の近似関数ｆ（ｘ，ｘj）と同様に、以下に示すＣ２クラスのＣＳ−ＲＢＦｓを用いる。
【数３５】
2010014619.tif 000034
【００４１】
また、（数３４）式の微分形は次式となる。
【数３６】
2010014619.tif 000035
【００４２】
なお、近似関数ｆ（ｘ，ｘm）としては、上記の（数３５）式以外に、３次元の場合には以下に示すＣ０クラスのＲＢＦである（数３７）式、Ｃ２クラスのＲＢＦである（数３８）式、Ｃ４クラスのＲＢＦである（数３９）式、Ｃ６クラスのＲＢＦである（数４０）式も用いることが可能である。
【数３７】
2010014619.tif 000036
【数３８】
2010014619.tif 000037
【数３９】
2010014619.tif 000038
【数４０】
2010014619.tif 000039
【００４３】
[解析手順]
図１は、本実施の形態である地下水流動解析方法の解析手順を示す概略フローチャートである。以下、フローチャートを参照しながら解析手順について説明する。
【００４４】
まず、解析領域の境界条件、境界値データ、境界節点データ、内部節点データ等の各種データの入力、及び、透水係数の空間分布を定めるための点ｘm（ｍ＝１〜Ｍ）の透水係数Ｋ（ｘm）を入力することにより、解析モデルを作成する（ステップＳ１）。次いで、透水係数の分布を表す（数３４）式におけるγの計算を行う（ステップＳ２）。次いで、透水係数の分布を表す（数３４）式と、その微分形である（数３６）式を、積分境界方程式の（数２６）式に代入し、（数３０）式の各マトリックスの計算を行う（ステップＳ３）。次いで、（数３０）式を解いて、未知境界のφ、∂φ／∂ｎを求める（ステップＳ４）。ステップＳ４で求めたφ、∂φ／∂ｎから、解析領域内の各境界節点及び各内部節点での流速ベクトルを求める。なお、ステップＳ４で求めたポテンシャルφ及びステップＳ５で求めた流速ベクトルは、ディスプレイやプリンタ等の出力手段を通じて出力を行うことができる。
【００４５】
以上のステップＳ１〜Ｓ５は、コンピュータと、そのコンピュータに実行させるプログラムによって実現することができ、そのプログラムは、コンピュータが読み取り可能な磁気ディスク、光ディスク、半導体メモリ等の記録媒体に格納することができる。
【００４６】
[計算例]
以下、上述した地下水流動解析方法の計算例について説明する。
【００４７】
（解析モデル）
図１０に示すような地層処分施設において放射性廃棄物の埋設箇所周辺を想定した比較的広い範囲の解析領域を対象とする。図３は、解析領域の寸法と境界条件を示す図である。解析領域の大きさは、図３に示すように、左側方１０００ｍ、右側方１００１ｍ、底盤１０００ｍの台形とした。境界条件は、左側方を静水圧条件（ポテンシャルφ＝１０００）、右側方を静水圧条件（ポテンシャルφ＝１００１）、底盤を不透水境界（∂φ／∂ｎ＝０）、上部境界面を水面（ポテンシャルφ＝１０００＋０．００１ｘ）とした。
【００４８】
図４は、以下で説明する比較例１と実施例１の二重相反境界要素法で用いる節点の配置を示す図である。図４に示すように、解析領域の境界と内部に縦４０個×横４０個の節点を配置した。図４において、白抜きの丸が境界節点であり、白抜きの四角が内部節点である。また、図６は、以下で説明する比較例２の有限差分法で用いる要素（メッシュ）を示す図である。図６に示すように解析領域を４０×４０の要素（メッシュ）に分割した。
【００４９】
上記の解析モデルを用いて、透水係数の分布が周期性を有する場合について、以下の３種類の計算を行った。
【００５０】
（比較例１）
上記解析モデルを用いて二重相反境界要素法（ＤＲＢＥＭ）による計算を行った。透水係数Ｋ（ｘ）の分布が周期性を有する場合として次式を用いた。
【数４１】
2010014619.tif 000040
【数４２】
2010014619.tif 000041
（数３７）式及び（数３８）式と、（数３７）式の微分形を（数２６）式に代入し、（数３０）式を解いて、図４に示す各境界節点及び各内部節点におけるポテンシャルφを求めた。また、各ポテンシャルφから各境界節点及び各内部節点における流速ベクトルを求めた。
【００５１】
（比較例２）
上記解析モデルを用いて有限差分法（ＦＤＭ）による計算を行った。透水係数の変化は、（数３７）式及び（数３８）式を用いて、図６に示す個々の要素の中心の座標に対して与えた。図６に示す各要素中心におけるポテンシャルφ及び流速ベクトルを求めた。
【００５２】
（実施例１）
上記解析モデルを用いて、二重相反境界要素法（ＤＲＢＥＭ）による計算を行った。図５は、解析領域に透水係数を与える複数の点ｘm（白抜き三角）を示す図である。（数３７）式と実質的に同じ周期をもつ透水係数Ｋ（ｘm）を、図５に示す各点に与え、（数３４）式の透水係数分布Ｋ（ｘ）を求めた。（数３４）式における近似関数は（数３５）式を用いた。（数３４）式とその微分形である（数３６）式を（数２６）式に代入し、（数３０）式を解いて、図４に示す各境界節点及び各内部節点におけるポテンシャルφを求めた。また、各ポテンシャルφから各境界節点及び各内部節点における流速ベクトルを求めた。
【００５３】
比較例１、比較例２、実施例１で求めたポテンシャルφのコンター分布を図７に示す。図７に示されるように、右下底盤付近を除き、３つの解析結果がほぼ一致することが分かる。
【００５４】
二重相反境界要素法をそれぞれ適用した実施例１及び比較例１で求めた流速ベクトルを図８に示す。図７に示したポテンシャルφのコンター分布と同様に、右下底盤部分付近で流速ベクトルに若干の相違が見られるが、全体として両者はほとんど差がないことが分かる。
【００５５】
また、二重相反境界要素法を適用した実施例１と有限差分法を適用した比較例２の流速ベクトルを図９に示す。流速ベクトルは、実施例１では内部節点上、比較例２では要素中心で表示しているため、両者の流速ベクトルはずれているが、両者の流速ベクトルの大きさと方向は整合していることが分かる。
【００５６】
上記の解析結果から、本実施の形態である地下水流動解析方法（実施例１）の有効性を確認することができた。
【００５７】
以上説明したように、本実施の形態の地下水流動解析方法及びそのプログラムによれば、二重相反境界要素法を用いて、不均質な透水係数の分布を考慮した解析を行う際に、解析領域内の特定の点に透水係数のデータを与えるだけで解析を行うことが可能になるため、従来のように透水係数の複雑な分布をフーリエ級数等の関数としてプログラムする必要がなくなる。その結果、解析データの作成作業の簡易化を図ることができ、効率的に解析を行うことができるようになる。
【００５８】
また、対象とする地層ごとに、透水係数の複雑な分布をフーリエ級数等の関数としてプログラム必要がないため、複数の地層の解析を短時間で効率よく行うことができる。
【００５９】
また、透水係数Ｋ（ｘ）の近似関数ｆ（ｘ，ｘm）としてＣＳ−ＲＢＦｓを用いたことで、透水係数Ｋ（ｘ）の微分値を精度よく求めることができる。
【００６０】
また、有限要素法や有限差分法のように、解析領域を分割した個々の要素に透水係数を設定するためのデータの加工を行う必要がない。
【産業上の利用可能性】
【００６１】
本発明の地下水流動解析方法は、地層処分の安全性の評価を行うに際し、放射性核種が地下水により輸送されることによって地表近傍の生物圏へ及ぼす影響を評価するための地下水流動状況の把握に有効に適用することができる。
【図面の簡単な説明】
【００６２】
【図１】本実施の形態の地下水流動解析方法における解析手順を示したフローチャートである。
【図２】本実施の形態の地下水流動解析方法で適用する透水係数Ｋ（ｘ）の空間分布を模式的に示した図である。
【図３】計算例で適用する解析領域の寸法と境界条件を示す図である。
【図４】二重相反境界要素法で用いる節点の配置を示す図である。
【図５】解析領域に透水係数Ｋを与える複数の点を示す図である
【図６】有限差分法で用いる要素（メッシュ）を示す図である。
【図７】実施例１、比較例１、比較例２で求めたポテンシャルφのコンター分布を示す図である。
【図８】二重相反境界要素法を適用した実施例１及び比較例１で求めた流速ベクトルを示す図である。
【図９】二重相反境界要素法を適用した実施例１及び有限差分法を適用した比較例２で求めた流速ベクトルを示す図である。
【図１０】放射性廃棄物を埋設する地層処分施設の概略図である。
【図１】
2010014619.tif 000045
【図２】
2010014619.tif 000046
【図３】
2010014619.tif 000047
【図４】
2010014619.tif 000048
【図５】
2010014619.tif 000049
【図６】
2010014619.tif 000050
【図７】
2010014619.tif 000051
【図８】
2010014619.tif 000052
【図９】
2010014619.tif 000053
【図１０】
2010014619.tif 000054