(57)【特許請求の範囲】
【請求項１】
浸透流方程式を適用した、地下水浸透流数値解析の有限要素法モデル化で岩盤層要素を通過するように薄い地層や割れ目が存在する場合、該薄い地層や割れ目を薄層要素として所定の透水係数を設定して解析要素モデル化を行うようにした地下水浸透流数値解析のモデル化方法において、前記解析要素モデル化の処理は、
前記薄層要素と、該薄層要素が通過する前記岩盤層要素とを、一定厚さｔeの三層構造モデルに換算するステップと、
該三層構造モデルの一定厚さｔeを、該岩盤層要素の体積の総和Ｖeを前記薄層要素の片面積Ａfで除して求めるステップと、
一定厚さｔeが求められた前記三層構造モデルに対して、各層の厚さｔiと各層の透水係数ｋiをもとに、単一層の層構造に沿う方向（ｓ）及びそれに直交する方向（ｎ）の等価透水係数ｋs，ｋnである下式
0004352199.tif 000007
を適用するステップと
を含み、前記岩盤層と、岩盤層を通過する薄い地層や割れ目とを一定厚さの要素モデルとして、該要素モデルに所定の等価透水係数が設定されるようにしたことを特徴とする地下水浸透流数値解析におけるモデル化方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は浸透流方程式を適用した、地下水浸透流数値解析におけるモデル化方法に係り、特に岩盤層の母岩に薄い地層や割れ目が存在している場合の地下水浸透流の数値解析において、薄い地層等の材料のモデル化を簡便に行えるようにしたモデル化方法に関する。
【０００２】
【従来の技術及び発明が解決しようとする課題】
以下の（１）式で支配されるような濃度拡散場や熱伝導場、地下水浸透流等の数値解析において、有限要素法や有限差分法などの数値解析手法が広く利用されている。
【０００３】
［数２］
0004352199.tif 000002
【０００４】
これらの解析手法には、解析領域を有限要素や格子によって分割したモデル（メッシュ）が必要となるが、そのメッシュの作成が一連の解析プロセスの中で最も煩雑な人手を要する作業となっている。特に地下水挙動の数値解析では、実際の地層中に非常に薄い地層や割れ目が存在している場合が多く、それを忠実に反映したモデル（メッシュ）化するのに多大な労力を要している。
【０００５】
有限要素法地下水解析において、この種の薄い地層や割れ目を効率よく考慮できるモデルとしてスメアード割れ目モデルが提案されている。このモデルは、地層内等に存在する割れ目を要素としてモデル化するのではなく、その割れ目が通過する母岩要素の透水係数を、提案式における面積比例配分により、母岩要素の透水係数を、母岩と割れ目の断面積を用いて平均化し、等価な物性値に変換して解析する手法である。
【０００６】
しかしながら、このスメアード割れ目モデルでは、母岩の各要素ごとに断面積の平均化を行う必要があり、三次元的に分布する複雑形状の割れ目（薄い地層）の場合、透水係数の平均化の計算は容易ではない。
【０００７】
そこで、本発明の目的は上述した従来の技術が有する問題点を解消し、薄い層を簡易にかつ精度良くモデル化できるようにした浸透流方程式を適用した数値解析におけるモデル化方法を提供することにある。
【０００８】
【課題を解決するための手段】
上記目的を達成するために、本発明は浸透流方程式を適用した、地下水浸透流数値解析の有限要素法モデル化で岩盤層要素を通過するように薄い地層や割れ目が存在する場合、該薄い地層や割れ目を薄層要素として所定の透水係数を設定して解析要素モデル化を行うようにした地下水浸透流数値解析のモデル化方法において、前記解析要素モデル化の処理は、前記薄層要素と、該薄層要素が通過する前記岩盤層要素とを、一定厚さｔeの三層構造モデルに換算するステップと、該三層構造モデルの一定厚さｔeを、該岩盤層要素の体積の総和Ｖeを前記薄層要素の片面積Ａfで除して求めるステップと、一定厚さｔeが求められた前記三層構造モデルに対して、各層の厚さｔiと各層の透水係数ｋiをもとに、単一層の層構造に沿う方向（ｓ）及びそれに直交する方向（ｎ）の等価透水係数ｋs，ｋnである下式
0004352199.tif 000003
を適用するステップとを含み、前記岩盤層と、岩盤層を通過する薄い地層や割れ目とを一定厚さの要素モデルとして、該要素モデルに所定の等価透水係数が設定されるようにしたことを特徴とする。
【００１３】
【発明の実施の形態】
以下、本発明の浸透流方程式を適用した数値解析におけるモデル化方法の一実施の形態について、添付図面を参照して説明する。
【００１４】
本発明は上記の（１）式で支配される濃度拡散場や熱伝導場、地下水浸透流の解析のモデル化において、その薄い材料を簡便なモデルとして取り扱うための手法であるが、本実施の形態では、説明を具体的、かつ容易にするために、地下水挙動に関する数値解析を例に説明する。
図１（ａ）は、地層内に薄い地層１０が斜行するように存在するモデルを示している。同図に示したように、四角形要素のうち、ハッチをつけた要素２０内を、濃色で示した薄い地層１０が通過している。このとき、図１（ｂ）に示したように、斜めに階段状をなすハッチをつけた要素２０の平均的な厚さｔeを求め、厚さｔfの薄い地層１０が挟在する一定厚さｔeの地層２５に換算する。
【００１５】
このとき薄い地層１０が通過するハッチをつけた要素２０の平均厚さｔeは、薄い地層１０の厚さｔfがほぼ一定であるとすれば、要素２０の体積の総和をＶe、薄い地層の片面積（周辺母岩との接触面積）をＡfとすると次（２）式で概算できる。
【００１６】
［数４］
0004352199.tif 000004
【００１７】
したがって、解析モデルは、一定厚さｔeの母岩中に厚さｔfの薄層が挟在する三層構造と見なすことができる。図１（ｂ）に示した一定幅の地層２５に換算した図１（ａ）の薄い層が通過する要素２０に対しては、Leonardによって示された以下の（３），（４）式による成層構造を成す地層の等価透水係数を適用することが好ましい。
【００１８】
［数５］
0004352199.tif 000005
【００１９】
ここで、本発明の材料特性としてのｋsとｋnは、図２に示したような成層構造の地層を単一層と見なした場合の層構造に沿う方向とそれに直交する方向の等価透水係数に相当する。また、ｔiは各地層の厚さ、ｋiは各地層の透水係数である。このとき各地層は等方、均質モデルとしている。
【００２０】
したがって、図１（ａ）に示した薄い層１０と薄い層が通過する要素２０を図１（ｂ）のような三層構造の地層２５とみなした場合、要素２０の等価透水係数は次式、（５），（６）式により計算できる。
【００２１】
［数６］
0004352199.tif 000006
【００２２】
ここで、ｋoは薄い地層を取り巻く母岩の透水係数であり、ｋfは薄層の透水係数である。ｋsは図１（ａ）に示した薄い層が通過する要素２０の薄層に沿う方向の透水係数、ｋnはそれに直交する方向の透水係数である、ただし、薄層の透水係数ｋfが母岩の透水係数ｋoがより小さい場合は、ｋsの寄与を排除するためｋsをゼロとおく。これは薄層の透水性が周辺母岩より小さい場合は、薄層に直交する方向の透水係数ｋnが支配的になるが、図１左図に示したように、実際には三層構造は凸凹のメッシュでモデル化されているためｋsが解析に大きく寄与してしまうことを排除するためである。
【００２３】
【実施例】
この発明の有効性を示すため、図３に示した三層構造の解析モデルを用い、厚さ５ｍの薄層の材料特性としての透水係数ｋfが周辺母岩の透水係数ｋoと大きく異なる場合の数値解析を実施した。
［解析モデル］
解析モデルのメッシュとしては、一辺が１０ｍの立方体で分割した図４に示したメッシュ（モデル１）、一辺が２０ｍの立方体で分割した図５に示したメッシュ（モデル２）、及び比較のため薄層を厚さ５ｍの要素で忠実に分割した図６のメッシュ（モデル３）を想定した。境界条件は、モデルの両端で全水頭（ｈ）を拘束し、その他の側面は不透水条件とした（図３参照）。
［解析結果］
(1)ｋf＝100ｋoの場合
図７はｋf＝100ｋoにおける解析結果である。同図は、図３に示した評価断面でのｈ（全水頭）の分布を示している。また同図には、モデル１，モデル２，モデル３の解析結果（等値線）が描かれている。３つの解析結果はよく一致しており、本発明における薄層のモデル化方法の有効性が確認できる。
(2)ｋf＝1/100ｋoの場合
図８はｋf＝1/100ｋoの場合の解析結果である。図中、モデル１，モデル２，モデル３の結果はいずれも良好な一致を示している。なお、図９は、図５に示した解析において、ｋs＝０と置かずに計算した結果であり、モデル１とモデル２の解析結果がモデル３の結果から外れていることが確認できる。
【００２４】
なお、本発明では前述した（５），（６）式の等価透水係数の計算において、薄い地層が通過する要素の体積の総和Ｖe、薄い地層の片面積（母岩との接触面積）Ａｆの値が必要となるが、これらの値は、公知の有限要素法プリプロセッサ等により簡単に算出することができる。これらのＶe，Ａfが得られれば、（５），（６）式の計算は非常に容易であり、実用性の高い方法といえる。
【００２５】
また、単一形状の立方体あるいは直方体要素で解析領域を分割するボクセル（Voxel）解析手法では薄層を忠実にモデル化できないので、ボクセル解析に本発明を適用することは特に有効である。
【００２６】
さらに、本発明では作成したメッシュデータを変更することなく、種々の薄層を考慮できるので、薄層の位置を変化させた設計検討も容易になる。
【００２７】
以上の実施の形態、実施例では、地層内に薄層を有する場合の地下水挙動の解析において、透水係数を材料特性として本発明の説明を行ったが、前述の（１）式で支配される濃度拡散場や熱伝導場等の解析対象においては、拡散係数、熱伝導率等が材料特性として適用される。
【００２８】
【発明の効果】
以上に述べたように、上記の（１）式で支配される地下水浸透流の解析のモデル化において、母岩に対する薄層や割れ目の状態を簡易にモデル化することができるという効果を奏する。
【図面の簡単な説明】
【図１】本発明による浸透流方程式を適用した数値解析におけるモデル化方法の一実施の形態を示したモデル図。
【図２】本発明の浸透流方程式を適用した数値解析におけるモデル化方法における成層構造をなす地層の等価透水係数の設定例を示したモデル図。
【図３】本発明において、薄い地層の検討用解析モデル図。
【図４】本発明のモデル化方法により、立方体分割によりモデル化された分割メッシュ（モデル１）
【図５】本発明のモデル化方法により、立方体分割によりモデル化された分割メッシュ（モデル２）
【図６】従来のモデル化方法により、薄層が忠実にモデル化された分割メッシュ（モデル３）
【図７】解析結果図（ｋf＝100ｋoの場合）
【図８】解析結果図（ｋf＝1/100ｋoの場合）
【図９】解析結果図（ｋf＝1/100ｋo，ｋs＝０と置かない場合）
【符号の説明】
１  地層モデル図
１０  薄い地層
２０  薄い地層が通過する要素（ハッチをつけた要素）
２５  一定幅の厚さの地層
ｋs，ｋn  材料特性（等価透水係数）
ｔe  薄層が通過する要素の平均層厚さ
Ｖe  薄層が通過する要素の体積の総和
Ａf  薄層の片面積
【図１】
0004352199.tif 000008
【図２】
0004352199.tif 000009
【図３】
0004352199.tif 000010
【図４】
0004352199.tif 000011
【図５】
0004352199.tif 000012
【図６】
0004352199.tif 000013
【図７】
0004352199.tif 000014
【図８】
0004352199.tif 000015
【図９】
0004352199.tif 000016