論文報告

公益社団法人地盤工学会中国支部論文報告集
地盤と建設Vol.37,No.1,2019

単孔を利用した透水試験方法の非定常法における直線勾配法の解釈に関する検討
Interpretation of the Straight-line Technique of the Transient Method for Hydraulic Conductivityof Aquifer by Single Borehole Testing小松満
中本萌

Mitsuru KOMATSU （岡山大学大学院環境生命科学研究科）
Moe NAKAMOTO （岡山大学大学院環境生命科学研究科）

単孔を利用した透水試験方法の非定常法における結果の整理方法のうち，直線勾配法により透水係数を求める理論式である Hvorslev 式は，その導入過程から初期定常状態からの試験データを適用するのが本来の形である．しかし，現在の地盤工学会基準の試験方法は，測定用パイプ内の水位を瞬時に変化させるスラグ試験が示されており，理論式の誘導過程の条件とは必ずしも一致しない．本研究では，被圧帯水層を対象とした数値解析により，スラグ試験結果を Hvorslev 式で評価した場合に求まる透水係数に含まれる誤差について，定常試験後に連続して非定常試験を行う場合と比較する形で検討した．

キーワード：原位置透水試験，透水係数，直線勾配法，スラグ試験 （IGC：D-04)
1．はじめに

水試験（以下，単孔式透水試験と称する）があり，揚水試験は揚水井と複数の観測井を用いた一定揚水時の観測

地下水に関する諸問題は以前より存在しているが，近井における水位低下量から透水係数を求める試験である．

年それらはますます複雑化している．地盤の浸透特性に

帯水層全体の平均的な透水特性を評価することができる

基づき，実際の地下水の挙動を正確に把握するためのシ

上に貯留性も得られ，信頼性が高い試験法であるものの，ミュレーション技術は電気計算機の汎用化及び高精度化

経済性の理由から比較的大規模な工事の際に実施されて1)．数値解析手法による地下水

いる．一方，単孔式透水試験は 1 本のボーリング孔や井シミュレーションは，モデルの選定・設定を適切に行な戸を用い，測定用パイプ内で揚水あるいは注水した際のえば，理論解で解くことができない複雑な地盤条件や施地下水位の回復または低下挙動から透水係数を求める試工条件等での設計が可能となる有効なツールであるが，験である．ボーリング孔先端部，スクリーン設置部など

実際にはこれらの条件に即した解析モデルが使われていの局所的な地盤の透水性を評価することができ，深さやなかったり，入力すべき情報が十分でないにも関わらず位置を変えて複数回試験を行うことで三次元的評価が可

高度な解析手法が使われていたりすることから，不適切能になる．しかし，透水係数の算定式の前提条件や境界な解析結果であっても定量的な評価は信頼されやすく，条件が満たされていない場合，信頼性が低下する懸念が設計と施工の乖離が生じることが懸念される 2)．

ある 5）．つまり，試験方法や試験孔仕様の違いによる点，に伴って高度化している

地下水シミュレーションを実施するにあたり，その入

さらにはこれらに基づく理論解を必ずしも正しく適用し力情報の一つである透水係数は，地盤の特性を示す重要ていない点などが信頼性が低下する要因として挙げられなパラメータである

3)．最近は，理工学のさまざまな分

る．

野でデータ同化や逆問題の研究が活発に行われており，

単孔式透水試験の試験方法には定常法と非定常法があ

不均質地盤における透水係数の分布の評価に適用されたるが，このうち非定常法における結果の整理方法は曲線事例が報告されている 4）．しかし，これらの技術の実用一致法（Cooper 法）と直線勾配法（Hvorslev 法）が用い化には，信頼性と計算速度を高めることが必要であり，られる．曲線一致法が試験開始から平衡水位に達するま

現状での実務において，透水係数は室内透水試験及び原でのデータを用いるのに対し，直線勾配法は試験開始直位置透水試験により直接測定されている 5)．このうち，後のデータを用いることから，平衡水位までのデータを室内透水試験ではサンプリングした試料に対して諸条件用いないため，試験が早く終了し，結果の整理方法も簡の影響を除外した正確な透水係数を求めることができる

易である等のメリットがある．しかし，水位差と時間のが，実際の地盤の不均一性を考慮することができない．

関係に直線が見られない場合には適用できない．元々，

一方，原位置透水試験では気象条件や地下水変動等の影

この直線勾配の理論式である Hvorslev 式は，定常式に非響を受けるものの，より実際の地盤状況に即した透水係

定常項を導入して誘導されており，初期定常状態からの数を求めることが可能である．

試験データを適用するのが本来の形であると言える 7)．

原位置透水試験

6）には，揚水試験と単孔を利用した透

しかし，現在の学会基準（JGS1314-2012） 6) では，測定－ 53 －

小松・中本用パイプ内の水位を瞬時に変化させるスラグ試験が示さ

非定常法の試験結果の整理のうち，直線勾配法

れており，理論式の誘導過程の条件とは必ずしも一致し（Hvorslev 法）は，図-2 に示すように片対数グラフに平ない．そこで本研究では，被圧帯水層を対象とした数値

衡状態からの水位差 s (cm)と時間 t (s)をとって log10s‒t

解析により，スラグ試験結果を Hvorslev 式で評価した場曲線を作成し，得られた直線の勾配から透水係数を算定合に求まる透水係数に含まれる誤差について，定常試験する方法であり，試験結果が地盤の貯留性の影響を受け後に連続して非定常試験を行う場合と比較して検討した． ていないと判断される場合に適用することができる．

2．現行の地盤工学会基準
2.1

単孔を利用した透水試験方法（JGS1314-2012）6)単孔を利用した透水試験方法には前述のように，定常
法と非定常法の 2 種類がある．図-1 にこれらの試験方法の概要を示す．定常法は揚水または注水を続け，測定用パイプ内の水位が一定になったときの流量を測定して地盤の透水係数を求める方法である．非定常法は，測定用パイプ内の水位を一時的に低下または上昇させ，平衡状態に戻る時の水位変化を経時的に測定する方法であり，

図-2log10 s-t 曲線の例 6)

スラグ試験と呼ばれる．非定常法は，透水係数が 10-2cm/s程度以上と予想される砂質・礫質地盤においては，水位2.2の経時変化が速く，計測しにくいため注意を要する．ま非定常試験の理論的背景

非定常試験の 2 種類の方法は，透水係数 k を求める理

た，試験時間内に初期水位の 90％程度まで水位が回復し， 論式の誘導における基礎式の元々の形が異なる．このう
この時間に少なくとも 10 点以上の測定データが得られ
ち，曲線一致法の Cooper 式は非定常試験時の地下水挙動る際に適している．

を無限の広がりを有する被圧帯水層内の二次元放射状流で表したときの支配方程式に基づいている．
� � � 1 �� � ������ � � �� � ��(1)

ここで，S：貯留係数

T：透水量係数である．

一方，直線勾配法の Hvorslev 式は，下式のように定常
式に非定常項を導入して誘導されている 7)．
��� � ��・����

ここで，K：透水係数(2)
H：孔内外水位差

F：形状係数，

A：試験孔内外水位変動域のパイプ断面積である．
式(2)は試験孔に流入する流量によって孔内水位が上昇した時，帯水層は孔内水位に即座に呼応して定常になる，すなわち時間遅れが無い状態を表している．帯水層(a) 定常法の貯留性が十分に小さい，または孔内貯留が十分に大きいと，孔内水位変動の影響が瞬時に帯水層に伝わるので，時間遅れが無いとみなすことができ，Hvorslev 式が成立する．一方，帯水層の貯留性が無視できないほど大きい，またはスラグ試験のように初期水位が定常状態でない場
合，孔内水位変動の影響が瞬時に帯水層に伝わり難いので，Hvorslev 式では不整合が生じる．よって，帯水層内では，少なくとも孔内水位に呼応した水圧分布になっていることが必要条件であることから，図-3 に示すような
初期定常状態から非定常試験によるデータを適用するのが適切であると考えることができる．しかし，実務で実施されている手法は，図-4 に示すように，測定用パイプ

図-1(b) 非定常法

内の水位を一時的に低下させるか上昇させるスラグ試験

試験方法の概要 6)について整理している事例が多く，理論式の誘導過程の条件とは必ずしも一致しない．
－ 54 －

単孔を利用した透水試験方法の非定常法における直線勾配法の解釈に関する検討

図-3 定常後の非定常試験概念図

図-5

解析モデル図表-1
帯水層透水係数
影響圏半径
初期水頭
初期水位差
井戸半径
試験区間長
層厚

解析条件
k=10-2cm/s
R=20m
H0=700cms=200cm
r=4.3cm
L=206.4cm
D=412.8cm
Q

Q

水面

圧力水頭

rc

図-4 スラグ試験概念図

rw

圧力水頭

k’

Cw

圧力水頭

Ss’

3．解析モデルおよび検討結果(a) 井戸貯留効果

3.1 解析モデル

図-6

図-5 および表-1 に本研究で用いた被圧帯水層の解析
モデル図と解析条件を示す．ここで，解析条件は任意の(b) 等価井戸材質

井戸貯留項と等価な井戸材質モデル 8)

3.2.2 微小要素の有無に対する検討値を用いた．また，井戸部分には図-6 に示すように井戸井戸境界での極端に不均質な透水性に対応するため，貯留と等価な特性を有する材質を与えた 8)．この理由と井戸底および井戸壁の内外に微小要素として縦横 1mm 幅

して，数値解析で井戸貯留を考慮した場合，井戸内水位の要素を入れた場合と入れていない場合の計算結果の違が非線形境界条件として働き，解の収束性が悪く，繰りいを検討した．さらに，井戸内の最も細かい要素の縦幅返し計算を行うことによる解析精度が低下することが懸を井戸内の最小要素縦幅と定義し，17.2，8.6，4.3，2.15，

念されるためである．なお，用いた解析コードは二次元

1.075cm の 5 段階に変化させて検討を行った．図-7 に井軸対称浸透流解析の AC-UNSAF2D 3)である．

戸内の最小要素縦幅と定常流量の変化を示す．井戸底，
井戸壁内外に微小要素を入れることで，流量の低下が見

3.2 定常解析

られ，理論解の流量 Q =568.25cm3/s に近づいた． した

3.2.1 定常時の流量計算

がって，井戸底および井戸壁内外に縦横 1mm 幅の微小要次式(3)で表される理論式から図-5 の解析モデルにおける定常流量を計算した結果，Q =568.25 cm3/s となった．
��

2����
4��
�
���� �
�素を入れたモデルとした．
3.2.3 マトリックス解法の選択
次に，マトリックス解法として PCG 法（間接解法）と(3)

ガウスの消去法（直接解法）の 2 つの解析手法を比較した．一般的に，ガウスの消去法の方が計算精度が良好と

ここで，�：水位低下量(cm)，d：孔径(cm)，�：試験区間

されているが，要素数が多くなるとマトリックスが大き長(cm)，�：流量(cm3/s)，λ：異方性度である．

くなり，計算時間を要する．そこで，精度に加え，計算に要する時間も比較した．図-8 に 2 つの方法による井戸－ 55 －

小松・中本内の最小要素縦幅と定常流量の変化を示す．結果として

600.0

定常流量に違いは認められなかった．なお，計算に要し599.5

た時間は，PCG 法が 1.1 秒，ガウスの消去法が 2.0 秒となった．計算時間もさほど変わりが無いため，以後の計流量(cm3/s)

599.0

算は精度が高いガウスの消去法とした．
3.2.4 井戸内の透水係数に関する検討

598.5
598.0

定常状態で井戸境界で同じ流量になるためには，透水係数が非常に高い井戸内では動水勾配が小さくなる．し

流入

597.5
597.0

低下する．そこで，井戸内透水係数 k’を 1.0×105cm/s から 1.0×104cm/s に変更させた場合の差異を検討した結果

流出

597.44

たがって，水頭差が非常に小さくなるために計算精度が
0

2

4

6

を図-9 に示す．透水係数を低下させると僅かではあるが8 10 12 14 16 18 20
最小要素縦幅(cm)(a) PCG法

流量の低下があり，理論解の流量に近づいた．また，井600.0

戸内節点の全水頭分布を確認し，井戸内が静水圧分布に599.5なっていることを確認した．ここで，井戸内透水係数が1.0×104cm/s で井戸内の最小要素縦幅が 1.075cm であると流量(cm3/s)

599.0

き，流量の解析値と理論値の誤差は約 5％となった．
3.2.5 解析モデルと解析条件の決定
これらの検討結果を踏まえ，解析モデルと解析条件を，

598.5
598.0
流入

①縦横 1mm 幅の微小要素を井戸底および井戸壁に作成す

597.5

る，②マトリックス解法にはガウスの消去法を使う，③597.0

井戸内透水係数は k’=1.0×104cm/s とする，④井戸内の最流出597.4402

4

小要素縦幅は 1.075cm とする，の各条件を踏まえて定常流量を解析した結果，Q=597.23cm3/s となり，倍精度で6 8 10 12 14 16 18 20
最小要素縦幅(cm)(b) ガウスの消去法の数値解析結果と理論値との誤差は約 5％であった．

図-8

解析手法による定常流量の変化

600.0
600.0

599.5

599.0

599.0

流量(cm3/s)

流量(cm3/s)

599.5
599.22

598.5

598.5
598.0

598.0
597.5
597.0

597.0

流出
0

2

4

6

8 10 12 14 16 18 20
最小要素縦幅(cm)2

4

6 8 10 12 14 16 18 20
最小要素縦幅(cm)

599.5

599.5

599.0
流量(cm3/s)

流量(cm3/s)

0

600.0

600.0

599.0
598.5
598.0

598.5
598.0
流入

597.5

0

2

597.0

流出

597.444流出597.23流入

597.5

6 8 10 12 14 16 18 20
最小要素縦幅(cm)

0

2

4

6

8 10 12 14 16 18 20
最小要素縦幅(cm)(b) 井戸内透水係数 k’=1.0×104cm/s
図-9(b) 微小要素あり
図-7

流出

597.44(a) 井戸内透水係数 k’=1.0×105cm/s(a) 微小要素なし

597.0

流入

597.5
流入

微小要素の有無による定常流量の変化－ 56 －

井戸内透水係数の変化による定常流量の変化

単孔を利用した透水試験方法の非定常法における直線勾配法の解釈に関する検討
3.3 非定常解析(2) 帯水層の比貯留係数Ssによる変化

3.3.1 結果の整理方法

図-11(a)に定常からの非定常試験を対象に帯水層比貯地盤工学会基準では，試験時間内に初期水位差の90％留係数を 10-7～10-6cm-1 まで変化させた際の log10s-t グラ

程度まで水位が回復することが望ましい 6)とされている． フの挙動の初期変化部分を示す．帯水層の比貯留係数 Ss
ここで，定常解析と同条件で式(4)に示すHvorslev法の非が大きくなると，Hvorslev 理論式の直線から外れる傾向

定常式（以下，Hvorslev理論式）から90％回復するまでが見られた．ここで 50%回復までに限定すると，Ss=10-6の時間を求めると46.97秒となった．したがって，2.5秒

cm-1 の場合でもほぼ一致しているとみることができる．

間隔で60秒まで計算を行い，その結果を基にlog10s-tグラ

また，Ss をさらに大きな値である 10-6～10-4cm-1 の間で変フを作成した．そして，種々の条件を変化させながら試化させた際は，同様に帯水層の比貯留係数 Ss が大きくな行を繰り返し，Hvorslev理論式との比較を行った．

ると直線から外れることが分かった．また，50%回復ま��

�2.3�� ��8�

2� ������ �� ������ � ��� � ��
�(4)での範囲でも一致しないことが判明した．一方，図-11(b)より，スラグ試験でも定常からの非定常試験と同様，帯

ただし，�/�≥4．ここで，de：有効直径（cm）（測定用パ

水層の比貯留係数 Ss が大きくなると直線から外れる傾

イプ内の断面積d（cm）から水位測定ケーブルの断面積

向が見られた．しかし，いずれも試験開始直後から円弧を差し引いて求めた有効断面積と等価な面積を有する円を描き，50%回復までの範囲に限定しても直線に一致しの直径），D：試験区間の孔径あるいは測定用パイプのスているとは言い難い結果を示した．また，Ss をさらに小クリーン外径（cm），L：試験区間の長さ（cm）である．

さい値である 10-7～10-10cm-1 の間で変化させた際は，直3.3.2 検討項目

線性が見られることが確認できた．しかし，これらは実(1) 試験方法による違い

際の地盤では取り得ないような小さな値である．

定常後の非定常試験およびスラグ試験での水位変化の違いとHvorslev理論解との差異を確認した．
(2) 帯水層の比貯留係数Ssによる違い

1

0

10

(3) 井戸内の比貯留係数Ss’による違い

60

定常からの非定常試験
定常→非定常水位差logs(cm)

復の様子の変化を確認した．

50

90％回復（s=20cm）

生じる．そのため，それぞれの試験方法で，一般的な砂3.0×10-7，5.0×10-7，1.0×10-6cm-1の1オーダー間での水位回40理論

帯水層の比貯留係数の影響で，水位回復に時間遅れが礫 層 の 比 貯 留 係 数 の 範 囲 で あ る 1.0×10-7 ， 2.0×10-7 ，時間t(s)
30

20

瞬時低下→非定常
スラグ試験
10
20

井戸内比貯留係数は，以下の式(5)で表される8)．
100

(5)

200

ここで，Cw：井戸貯留項，V：井戸内水面下体積，rc：ケ

図-10

ーシング半径，rw：井戸半径，L：ストレーナ長である．
90
90

井 戸 内 比 貯 留 係 数 を 小 さ く す る 試 行 （ rw ＞ r c ） で の

100
100

Hvorslev理論解との合致具合を確認する．なお，岩盤のような低い透水性の地盤では回復に時間を要すため，rw＞rc（d＞de）になるような試験装置を使う場合がある．

00

10
10

時間t(s)
時間t(s)
20
20-6 ～1×10-4cm/s程度が想定されている 6) ．そこで，

90
90

00

10
10

時間t(s)
時間t(s)
2030
30

水位差logs(cm)
水位差logs(cm)

図-10 に定常からの非定常試験とスラグ試験の水位変20
20

時間t(s)
時間t(s)

40
40

60
60

理論
理論
Ss=1×10(-6)(/cm)
Ss=1E-6(cm-1)
Ss=2×10(-6)(/cm)
Ss=2E-6(cm-1)
Ss=5×10(-6)(/cm)
Ss=5E-6(cm-1)
Ss=1×10(-5)(/cm)
Ss=1E-5(cm-1)
Ss=2E-5(cm-1)
Ss=2×10(-5)(/cm)
Ss=5E-5(cm-1)
Ss=5×10(-5)(/cm)
Ss=1E-4(cm-1)
Ss=1×10(-4)(/cm)
90％回復
90％回復
50％回復
50%回復

理論
理論
Ss=1×10(-7)(/cm)
Ss=1E-7(cm-1)
Ss=2×10(-7)(/cm)
Ss=2E-7(cm-1)
Ss=3×10(-7)(/cm)
Ss=3E-7(cm-1)
Ss=5×10(-7)(/cm)
Ss=5E-7(cm-1)
Ss=1×10(-6)(/cm)
Ss=1E-6(cm-1)
50％回復
50％回復

化の解析結果と Hvorslev 理論式の結果を示す．定常からの非定常試験では，初期段階は直線になり，時間が経過

00

20
20

時間t(s)
時間t(s)

40
40

理論
理論
Ss=1×10(-10)(/cm)
Ss=1E-10(cm-1)
Ss=1×10(-9)(/cm)
Ss=1E-9(cm-1)
Ss=1×10(-8)(/cm)
Ss=1E-8(cm-1)
Ss=1×10(-7)(/cm)
Ss=1E-7(cm-1)
90％回復
90％回復
50％回復
50%回復

(b) スラグ試験
図-11－ 57 －

90
90

200

200

するにつれ水位回復に遅れが生じた．一方，スラグ試験では初期段階から円弧を描くような曲線となった．

40
40

100
100

100

帯水層の透水係数がこの範囲での検討も実施した．
(1) 試験方法による違い

00(a) 定常からの非定常試験

細粒分を多く含んだ地盤を対象とし，透水係数が3.3.3 解析結果

90
90

200
200

200
200

式は，本来，試験区間の孔壁が自立するような，比較的
k=1×10

40
40

理論
理論
Ss=1×10(-7)(/cm)
Ss=1E-7(cm-1)
Ss=2×10(-7)(/cm)
Ss=2E-7(cm-1)
Ss=3×10(-7)(/cm)
Ss=3E-7(cm-1)
Ss=5×10(-7)(/cm)
Ss=5E-7(cm-1)
Ss=1×10(-6)(/cm)
Ss=1E-6(cm-1)
50％回復
50％回復

(4) 帯水層透水係数に対する検討
JGS1314-2012に示されている定常および非定常の理論

30
30

100
100

水位差logs(cm)
水位差logs(cm)rw＝rcの時，Ss’＝1/L=4.845×10-3cm-1となる．この値から，試験方法による違い(Ss=10-6cm-1)

水位差logs(cm)
水位差logs(cm)

��
��� �
�
�
��� � �水位差logs(cm)
水位差logs(cm)��� �

帯水層比貯留係数 Ss による変化

60
60

小松・中本
0

1

10

20

時間t(s)
30

40

50

rc を小さくし，早く水位を回復させる試行を行った．図60-12に式(5)の計算値であるSs’=4.845×10-3cm-1 とSs’=1.0×

水位差logs(cm)理論
定常からの非定常試験
定常→非定常10-3cm-1の結果をそれぞれ示す．井戸貯留項Cwを小さくすると理論式に一致しなくなることが判明した．

瞬時低下→非定常
スラグ試験
90％回復

(4) 帯水層の透水係数に対する検討

10

帯水層の透水係数を k=1.0×10-6，k=1×10-4cm/s として同様の検討を行った結果を図-13 に示す．透水係数を変化させても，同様の傾向が確認できた．

100
200

1

4．考察

(a) Ss’=4.845×10-3cm-1
0

10

20

時間t(s)
30

40

50

60

4.1 試験結果の妥当性の検討
通常，Hvorslev 法で試験結果を整理する際は初期水位

水位差logs(cm)

差の 90％程度まで水位が回復するまで試験を続ける 6)が，実際の現場では時間の制約上，試験を短時間で終了
10

理論
定常からの非定常試験
定常→非定常瞬時低下→非定常
スラグ試験
90％回復

させたい場合があることが想定される．そのため，これまでの解析結果に基づき，初期段階に限定した場合の試験結果の妥当性について検討した．まず，直線とみなす
基準として，Ss=1.0×10-6cm-1 の場合において定常後の非100定常試験で初期水位差からの回復率が 50％を超えた時200間までのデータに設定した．そして，この時の残差平方(b) Ss’=1.0×10-3cm-1
図-12

水位差logs(cm)

0E+0
1

1E+3

平均（MSE）を直線とみなす許容範囲とし，この MSE

井戸貯留項による変化
2E+3

時間t(s)
3E+3

4E+3の値を基準として定常後の非定常試験およびスラグ試験
5E+3

6E+3

結果の妥当性は，透水係数の誤差（kE）で評価し，許容理論解（Hvorslev式）
定常→非定常(Ss=1E-6)
スラグ(Ss=1E-6)
定常→非定常(Ss=1E-5)
スラグ(Ss=1E-5)
90％回復
50％回復

10

の直線とみなせる区間をそれぞれ決定した．また，試験範囲は理論値の 10％以内とした．
4.4.1 試験方法の違いによる試験結果の比較
表-2 に直線とみなせる限界の回復率とその回復率での透水係数の誤差の関係を試験方法ごとにまとめた．定20

常後の非定常試験は，許容範囲内の回復率では透水係数の誤差もほぼ許容範囲内であった．スラグ試験は，許容範囲内の回復率であっても，透水係数の誤差が許容範囲100外であることが判明した．また，同じ比貯留係数の試験

200

結果を比較すると，スラグ試験は定常後の非定常試験に(a) k=1.0×10-4cm/s

水位差logs(cm)

0E+0
1

1E+5

2E+5

時間t(s)
3E+5

4E+5

比べ直線視できる区間が小さくなり，求まる透水係数の
5E+5

6E+5

理論解（Hvorslev式）
定常→非定常(Ss=1E-6)
スラグ(Ss=1E-6)
定常→非定常(Ss=1E-5)
スラグ(Ss=1E-5)
90％回復
50％回復

10

誤差も大きくなる傾向が認められた．
表-2 限界回復率と透水係数の誤差の関係(a) 定常後の非定常試験(b) スラグ試験
Ss(cm-1)

20

1.0E-6
5.0E-7
3.0E-72.0E-7
1.0E-7100200(b) k=1.0×10-6cm/s
図-13

帯水層透水係数に対する検討

限界
回復率
(％)
53.1
58.7
63.7
71.4
100.0

kE
(％)

Ss
(cm-1)

11.7
9.7
8.1
7.1
4.1

1.0E-6
5.0E-7
3.0E-72.0E-7
1.0E-7

限界
回復率
(％)
38.5
45.1
44.1
50.1
54.7

kE
(％)
32.4
22.4
18.4
13.3
7.6

4.4.2 定常からの非定常試験結果の妥当性の検討
図-14に比貯留係数Ssと残差平方平均（MSE）が許容範(3) 井戸貯留項Cwによる変化

囲内に入る限界の回復率の関係，図-15に直線とみなせる

井戸貯留項Cwを理論値より小さく，つまり，有効半径範囲内（MSE≦1.7×10-5）での最大の透水係数の誤差kE－ 58 －

単孔を利用した透水試験方法の非定常法における直線勾配法の解釈に関する検討
および試験時間tの関係を示す．帯水層の比貯留係数Ssが囲内に入る限界の回復率の関係，図-18に直線とみなせる

大きくなると，直線視できる範囲は小さくなる．また，範囲内（MSE≦1.7×10-5）での最大の透水係数の誤差kE，

砂礫層が対象となる比貯留係数の範囲では，透水係数の試験時間tの関係を示す．定常からの非定常試験と同様，

誤差もほぼ許容範囲内で，試験時間も短くなっている．

帯水層比貯留係数Ssが大きくなると，直線視できる範囲-5，透水係数のは小さくなっている．また，定常からの非定常試験に比次に，図-16に残差平方平均MSE＝1.7×10

誤差kE＝10％以内に入る限界回復率をそれぞれ計算し，べ，直線性を満たす限界回復率の値は小さい．砂礫層の比貯留係数Ssとの関係を整理した．比貯留係数が大きく範囲で，直線視できる範囲でも透水係数の誤差が許容範なると，直線性の限界よりも先に結果の妥当性を示す透囲外となることが示唆された．

水係数の誤差が許容範囲内から外れる傾向にあることが判明した．

限界回復率(％)
100
90

限界回復率(％)

80

100

70

90

60

80

限界回復率
（MSE≦1.7E-5）

70

限界回復率
≦1.7E-5)（MSE（MSE≦1.7E-5）

50
40

60

30

50

20

40

10

30

0
1.0E-1020

1.0E-910

1.0E-81.0E-7砂礫

1.0E-6

帯水層比貯留係数Ss(/cm)

01.0E-7図-14

1.0E-6

砂礫

1.0E-5図-17

1.0E-4帯水層の比貯留係数と限界回復率の関係

帯水層比貯留係数Ss（/cm）（スラグ試験）帯水層の比貯留係数と限界回復率の関係

最大透水係数誤差kE（％）
100（定常からの非定常試験）

最大kE
最大
kE

90
試験時間（s）
60

50

試験時間
試験時間

80
最大透水係数誤差kE（％）
60

試験時間t(s)
60

70

40

60
50

50

50

最大kE

40

40

30
20

20

10

10

20
0
1.0E-10

砂礫

図-15

図-18

となった．つまり，試験開始直後は帯水層Ssの影響が少限界回復率（MSE≦1.7E-5）
結果の妥当性が先に許容範囲を超えるなく，試験後半になると貯留性の影響が大きくなる．定限界回復率（kE≦10％）常状態では，井戸付近の水位は下がっているので，初期の井戸内外の水位差は小さい．そのため，井戸内の水位

直線性が先に許容範囲を超えるが回復するのと同時に井戸付近の帯水層の水位も上昇す
る．この水位変化が貯留性の影響で遅れるため，井戸内

20外の水位差，すなわち動水勾配は比較的小さい状況が続

100
1.0E-7図-16帯水層の比貯留係数と透水係数の誤差および定常からの非定常試験では，試験開始直後はHvorslev

60

30

砂礫

理論の直線に概ね合致し，試験後半になると外れる結果

90

40

0
1.0E-6

1.0E-74.2 試験方法の違いによる解析結果に対する考察

限界回復率(%)
10050

1.0E-8

試験時間との関係（スラグ試験）

試験時間との関係（定常からの非定常試験）

70

1.0E-9帯水層比貯留係数Ss（/cm）

帯水層の比貯留係数と透水係数の誤差および80

透水係数誤差許容範囲

帯水層比貯留係数Ss(/cm)

01.0E-4

1.0E-51010透水係数誤差許容範囲
1.0E-6

20

30

30

0
1.0E-7

30

40

試験時間

く．そのため，井戸内への流入流量が少なくなることで，
1.0E-6
1.0E-5
帯水層比貯留係数Ss（/cm）

1.0E-4

水位回復が遅くなっていると考えられる．また，試験開帯水層の比貯留係数と限界回復率との関係

始直後にHvorslev理論の直線に合致するのは，初期段階
では流入流量が帯水層の貯留を回復させるより，孔に流4.4.3 スラグ試験結果の妥当性の検討入させる方が優先的に働いていることによる．その後，

図-17に比貯留係数Ssと残差平方平均（MSE）が許容範

孔に入る水の動水勾配が小さくなると，流入流量が帯水－ 59 －

小松・中本層の貯留を回復させるため，貯留性の影響が出てくるも非定常試験よりも水位が早く回復することから，算のと考えることができる．

定される透水係数は大きくなることが想定される．

一方，スラグ試験では，試験開始時に井戸内外の水位(5) これらの結果を踏まえると，Hvorslev式の適用は，差が大きい．試験開始直後に井戸内の水位は回復するが，定常からの非定常試験が望ましいが，定常状態を判

井戸付近の帯水層の水位は低下しようとする．この水位断する必要があることから場合によっては長時間低下が貯留性の影響で遅れるため，井戸内外の水位差やの試験となる可能性がある．そのため，実務での実動水勾配は比較的大きい状態が続く．そのため，井戸内施を鑑みると，必ずしも定常からの非定常試験が推への流入量が多い状態が続き，水位回復も早いままであ

奨されるということではなく，現在の地盤工学会基る．やがて帯水層の水位と井戸内の水位が一致すると，準の試験法であるスラグ試験を適用する場合は，算帯水層の水位上昇が貯留性によって遅れるため，井戸内定する透水係数にいくらかの誤差が含まれることへ流入する流量が少なくなり，回復に遅れが見られるとを認識した上で結果を評価することが必要である．

考えられる．ここで，帯水層の比貯留係数 Ss’が大きくなると，試験方法によらず直線から離れる．これは，理論

謝辞

式そのものに貯留性の影響を示す項が入っていないため

本研究を遂行するにあたり，日本原燃株式会社埋設事である．また，井戸内の比貯留係数 Ss を小さくすると，

業部の進士喜英氏に貴重な助言をいただきました．ここ

試験方法によらず直線から離れる．そのため，岩盤のよに記して謝意を表します．

うな低い透水性の地盤で直径 d のまま，孔内にケーブルを入れる等して，水位変動断面積を小さくしたり，試験参考文献

削孔は直径 d で掘り切り，試験区間上端に直径 de となる

1)

西垣誠：小特集「地下水流動解析の高度化手法と検証・確認」（その 2）の掲載にあたって，地下水学

内管を取り付けたりして，水位回復を早くする試験を実施すると，誤差が大きくなることが想定される．

会誌，Vol.57，No.3，pp.317-318，2015.
2)

高坂信章：地下水掘削工事における地下水シミュレーションの適用－その手順と留意点－，基礎工，Vol.46，No.6，pp.29-32，2018.

5．まとめ
3)
本研究では，被圧モデルを対象に，スラグ試験結果をHvorslev 式で評価した場合に求まる透水係数の誤差につ地盤工学会・地下水を知る編集委員会編：地下水を知る，丸善，p.16，2008.

4)

村上章：地盤工学における最先端のデータ同化／逆いて，定常試験後に連続して非定常試験を行う場合と比解析手法，地盤工学会誌，Vol.65，No.10，pp.1-5，

較して検討した．以下に得られた知見を列挙する．

2017.

(1) 定常からの非定常試験において，一般的な砂礫層の5)くの間は理論直線に合致し，透水係数の誤差も許容

河野伊一郎：地下水工学，鹿島出版会，pp.85-110，
1989.

比貯留係数での試験では，試験開始直後からしばら
6)

地盤工学会 地盤調査規格・基準委員会編：地盤調査の方法と解説－第 7 編

範囲内にある．その後は直線から離れていき，透水地下水調査，丸善出版，

pp.471-659，2013.

係数の誤差も大きくなる．
(2) 定常からの非定常試験において，帯水層の比貯留係

7)

源流 Hvorslev（1951），地下水学会誌，Vol.58，No.2，

数が大きくなると，直線性の限界よりも先に透水係

pp.241-246，2016.

数の誤差が許容範囲外になることが判明した．ただし，これらは透水係数を一定とした条件での解析結

進士喜英：地下水学の名著を読む(7)現場透水試験の8)

進士喜英，嶋村貞夫，濱野隆司：井戸貯留項と等価果であり．実際には比貯留係数が高い土砂地盤ではな井戸材質モデル，第 31 回地盤工学研究発表会，

透水係数が小さいことに注意する必要がある．

pp.2155-2156，1996.

(3) スラグ試験では，試験開始直後から理論直線には乗9)西垣誠：有限要素法による飽和不飽和浸透流解析

らない曲線を描く結果となった．ある程度水位が回-AC-UNSAF2D-プログラム解説およびユーザーマ

復すると，理論直線と平行な直線を示す．これは，ニュアル，地盤工学会関西支部，第 41 回実技セミ

あくまでも理論直線からは外れており，理論直線とナーテキスト「数値解析の理論と実務」, 1999.

合致するのは，かなり小さな帯水層の比貯留係数の10) 地盤工学会地盤調査規格・基準委員会：地下水調査

場合のみである．

に用いる井戸理論式の整理及び解説（2017 年度版），(4) 一般的な砂礫層の比貯留係数で透水係数の算定を行うと，初期段階でも誤差は10％を超えることから，

https://www.jiban.or.jp/wp-content/uploads/2017/10/0-i
ntroduction.pdf，（参照 2019.3.15).

スラグ試験をHvorslev法で整理する際は，誤差を含む可能性に注意する必要がある．また，定常からの－ 60 －(2019 年 6 月 17 日

受付)