地すべり安全率 Bishop 簡便法シミュレーター 戻る
地盤工学・防災

地すべり安全率 Bishop 簡便法シミュレーター

円弧すべり面を仮定した斜面の安全率 Fs を、Bishop 簡便法(1955)の反復計算で求めます。斜面高さ・斜面角度・粘着力・内部摩擦角・間隙水圧比・地震係数を変えて、降雨や地震が地すべりに与える影響を体感できます。

パラメータ設定
斜面高さ H
m
斜面角度 β
°
粘着力 c'
kPa
有効応力法での粘着力
内部摩擦角 φ'
°
単位体積重量 γ
kN/m³
間隙水圧比 r_u
降雨・地下水の影響。0.5 で水位が斜面表面まで
地震係数 k_h
擬似静的解析の水平震度。レベル1で0.10〜0.15
破壊形態
本ツールは Bishop 円弧すべりを中心に解析
計算結果
安全率 Fs (Bishop)
限界斜面角 (Fs=1) (°)
臨界 r_u (Fs=1)
抵抗モーメント
駆動モーメント
リスク評価
斜面断面・円弧すべり面・スライス分割

青線=円弧すべり面、点線=スライス境界、水色帯=地下水位、矢印=駆動/抵抗の力。色は Fs(緑=安定/橙=要注意/赤=破壊危険)。

降雨影響 — Fs vs 間隙水圧比 r_u
地震影響 — Fs vs 地震係数 k_h
理論・主要公式

$$F_s = \frac{\sum \left[c'b + (W - ub)\tan\phi'\right]/m_\alpha}{\sum W\sin\alpha + k_h W\cos\alpha},\quad m_\alpha = \cos\alpha + \frac{\sin\alpha \tan\phi'}{F_s}$$

Bishop 簡便法の反復解。c'=有効粘着力 [kPa]、φ'=有効内部摩擦角 [°]、W=スライス重量 [kN/m]、u=スライス底面の間隙水圧 [kPa]、b=スライス幅 [m]、α=スライス底面の傾斜角、k_h=水平震度。Fs を未知数として右辺に含むため、3〜5 回の反復で収束する。

$$u = r_u \cdot \gamma \cdot h, \qquad r_{u,\text{crit}} \approx \frac{c'/(\gamma H) + \tan\phi'\cos^2\beta - \sin\beta\cos\beta}{\tan\phi'\cos^2\beta}$$

間隙水圧比 r_u の定義と、Fs=1 を与える臨界 r_u(無限斜面近似)。降雨で地下水位が上がると u が増え、有効応力が減って摩擦抵抗が落ちる。

地すべり安全率 Fs — Bishop 簡便法による斜面安定解析

🙋
「地すべり安全率」って、ニュースで「大雨で斜面の安全率が下がった」とか聞きますけど、具体的にはどうやって計算するんですか?
🎓
ざっくり言うと、斜面を「すべり落ちようとする力(駆動力)」と「すべりに抵抗する力(抵抗力)」の比、それが Fs = 抵抗 / 駆動 だよ。Fs > 1 ならまだ抵抗が勝っていて安全、Fs < 1 なら駆動が勝っていて、いつ崩れてもおかしくない状態。実務では円形の「すべり面」を仮定して、その面に沿って斜面が一塊で滑り落ちるモデルで考えるんだ。これを「円弧すべり」、計算法では Bishop 簡便法(1955)が世界中で一番使われている。
🙋
Bishop さんって人の方法ってことですね。普通の力のつり合いと何が違うんですか?
🎓
斜面を「スライス」って呼ばれる細長いブロックに切って(このツールでは 10 個)、各スライスごとに重量・水圧・摩擦を計算するんだ。Bishop の特徴は、スライス側面の水平力を考慮しつつ、安全率 Fs 自体を未知数として反復計算するところ。式の右辺に Fs が m_α として入るから、最初 Fs=1.5 と仮定 → 計算 → 新しい Fs → また計算…と 3〜5 回回すと収束する。本ツールはこの反復を自動で 10 回やっているよ。
🙋
スライダーの「間隙水圧比 r_u」を 0.3 から 0.6 にしたら、Fs がガクッと 1.5 → 1.0 くらいまで下がりました!これが「雨で地すべりが起きる」ってやつですか?
🎓
まさにそれだ。降雨で地下水位が上がると、すべり面の間隙水圧 u が増える。Bishop の式の (W - ub)·tan φ' という項が「有効応力による摩擦抵抗」で、u が増えるとここが小さくなって Fs が落ちる。日本で地すべりの 7〜8 割が梅雨や台風シーズンに集中するのも、これが理由なんだ。r_u = 0.5 はだいたい「水位が斜面表面まで上がった」状態で、ここまで来ると本来 Fs=1.5 の斜面でも Fs=1.0 を割ることがある。だから雨量規制や避難勧告のしきい値はこの r_u をモニタリングして決められる。
🙋
地震係数 k_h ってのもありますね。これを 0.15 にしたら駆動モーメントが急に増えて Fs が落ちました。地震でも地すべり起きますもんね。
🎓
それが「擬似静的解析」と呼ばれる方法だよ。本当は地震動は時々刻々変わる波だけど、設計実務では「水平方向に重力の k_h 倍の慣性力が定常的にかかる」と単純化する。日本の道路土工指針だとレベル 1 地震動で k_h=0.10〜0.15、レベル 2 で 0.15〜0.20 が目安。Vajont ダム滑落(1963年、イタリア、約2000人犠牲)や台湾の 八八水災(2009年、小林村が壊滅)も、貯水位上昇+豪雨や台風という複合トリガーで起きた。詳細解析には FLAC や PLAXIS Dynamic を使うけど、まずこのツールのような擬似静的解析で「だいたいどれくらい Fs が下がるか」を当たり付けする、というのが標準的な手順だね。
🙋
Fs=1.5 なら安全と思っていいですか?安全率の目安はどれくらいですか?
🎓
日本の宅地造成等規制法や道路土工指針では、常時(無地震・通常水位)で Fs ≥ 1.5、地震時で Fs ≥ 1.0、降雨時で Fs ≥ 1.2 が一般的な目安だよ。ただし注意したいのは、土質パラメータ(c'、φ'、γ)にはバラツキが大きいこと。室内試験で c'=30kPa と出ても、実斜面では場所ごとに 20〜50kPa とばらつく。だから Fs=1.5 でも、パラメータが 20% 悪い側にずれただけで Fs=1.0 を割ることがある。実務では信頼性解析(Reliability-Based Design、破壊確率で評価)も併用するし、対策工も「Fs を 1.2 から 1.5 に上げる」のように余裕を持って設計するんだ。

よくある質問

Bishop 簡便法(1955)は、円弧すべり面の安全率 Fs を求める Limit Equilibrium 法のひとつです。Fellenius 法(簡易分割法)が各スライスのモーメントだけを取るのに対し、Bishop はスライス側面の水平力を考慮し、Fs を未知数として反復的に解きます。一般に Bishop のほうが Fellenius より 5〜20% 大きな Fs を与え、円弧すべりに対しては Spencer 法や Morgenstern-Price 法とほぼ同じ結果になります。本ツールは 10 スライスで 10 回反復しており、典型条件で 3〜5 反復で収束します。
間隙水圧比 r_u は、スライス底面での間隙水圧 u を、上載土の鉛直全応力 γh で割った無次元量です(r_u = u / (γh))。地下水位が斜面表面まで上がり、流れがない静水圧で満たされた極限で r_u ≈ 0.5、流れがある場合は 0.5〜0.7 に達します。本ツールの既定条件(H=20m, β=30°, c'=30kPa, φ'=25°)では、r_u = 0.3 で Fs ≈ 1.5、r_u = 0.6(豪雨直後)で Fs ≈ 1.0 近くまで低下します。これが「雨で地すべりが起きる」物理的理由です。
擬似静的解析では、地震動による水平加速度 a を重力 g で割った値 k_h = a/g を、各スライスに水平慣性力として加えます。日本の道路土工指針(NEXCO)では、レベル 1 地震動で k_h = 0.10〜0.15、レベル 2 で 0.15〜0.20 が一般的な目安です。米国 USACE は 0.10〜0.15 を推奨。本ツールで k_h = 0.15 を入れると、駆動モーメントが約 30〜40% 増え、Fs は静的値の 70% 程度に低下します。詳細な動的応答解析(FLAC、PLAXIS Dynamic)には及びませんが、設計初期の擬似静的判定には十分使えます。
Fs = 1.0 は「ぎりぎり釣り合っている」状態で、実務では決して安全とは言いません。日本の宅地造成等規制法・斜面安定基準では、常時(無地震・通常水位)で Fs ≥ 1.5、地震時で Fs ≥ 1.0、降雨時で Fs ≥ 1.2 が一般的な目安です。本ツールの「リスク評価」も Fs > 1.5 を安定、1.2〜1.5 を要注意、1.0〜1.2 を警戒、< 1.0 を破壊危険としています。さらに、土質パラメータには大きな不確実性があるため、Fs = 1.5 でも、c' や φ' が実測値の 20% 低めだとあっという間に 1.0 を割ることに注意してください。

実世界での応用

急傾斜地の防災・ハザードマップ:日本では土砂災害防止法に基づいて、全国の急傾斜地が「土砂災害警戒区域(イエローゾーン)」「特別警戒区域(レッドゾーン)」に指定されています。これらの区域指定や住宅造成の許認可では、本ツールのような Bishop 簡便法による Fs 計算が基本評価指標となります。国土地理院の土砂災害ハザードマップや、都道府県が公開する Web ハザードマップの背後では、地形・土質データと組み合わせた数千〜数万箇所の Fs 評価が走っているのです。

道路・鉄道の切土・盛土設計:NEXCO(高速道路)・JR各社の土構造物では、切土斜面の Fs を常時 ≥ 1.5、地震時 ≥ 1.0、降雨時 ≥ 1.2 で照査します。Fs が不足する場合の対策工は、軽い順から、(1) 排水工(横ボーリング・暗渠)で r_u を下げる、(2) 表面工(モルタル吹付・植生)、(3) アンカー工・抑止杭、(4) 抑え盛土。本ツールで「r_u を 0.5 から 0.3 に下げる」シミュレーションをすると、排水工による Fs 回復が定量的に見えます。

ダム・貯水池の斜面安定:貯水位の急上昇・急低下(Rapid Drawdown)は典型的な地すべりトリガーです。1963年のイタリア・Vajont ダム滑落(約2億 m³ の岩盤が貯水池に崩落、約 2000 人犠牲)はまさにこの現象でした。ダム上下流の斜面では、定常時・湛水時・地震時・急低下時の 4 条件で Fs 照査が義務付けられ、特に急低下時は間隙水圧の消散が追いつかず Fs が最小になります。本ツールで r_u を高くしたまま β を上げる感覚は、急低下時の不安定化と対応しています。

地すべりモニタリングと早期警報:近年は IoT 傾斜計・GNSS・伸縮計を斜面に設置し、リアルタイムで動きを監視するシステムが普及しています。観測された変位速度や雨量と、本ツールのような Fs 計算とを組み合わせ、「累積雨量 X mm を超え、変位 Y mm/h を超えたら避難勧告」というしきい値を設計します。広島県の土砂災害(2014年)、熱海伊豆山土石流(2021年)以降、こうした AI を含む統合警報システムへの投資が急増しています。

よくある誤解と注意点

まず最大の落とし穴が、「土質パラメータの不確実性を Fs 計算が隠してしまう」こと。本ツールに c' = 30 kPa と入れると Fs が一意に求まりますが、実斜面では場所・深さ・季節によって c' は 10〜60 kPa とばらつき、室内三軸試験と原位置試験でも 2〜3 倍違うのが普通です。Fs = 1.5 と計算されても、パラメータが下方 20% にずれただけで Fs = 1.0 を割ります。実務では信頼性解析(Reliability-Based Design、Monte Carlo シミュレーション)で破壊確率 P_f を併記したり、感度解析でパラメータごとの影響度を可視化するのが望ましい。決定論的 Fs 一つで「OK」と判定するのは危険です。

次に、「Bishop 簡便法ですべての斜面を解けると思い込む」こと。Bishop は円弧すべり面が前提なので、層構造が傾いている岩盤の平面すべり、断層や弱面に沿うウェッジ崩壊、流動性の高い土石流などには適用できません。地層が斜面と平行に堆積している風化岩・崩積土の斜面では、平面すべりや Janbu の一般化分割法のほうが適切です。複雑な地質ではFEM/FDM(PLAXIS、FLAC)による弾塑性解析が必要になります。本ツールのモード選択で「平面すべり」「ウェッジ」を選んでも、解析自体は Bishop モデルで実行している点に注意してください(教育用途では十分な近似)。

最後に、「擬似静的解析の k_h を絶対視する」こと。k_h = 0.15 で Fs ≥ 1.0 を満たせば耐震 OK とする設計は便利ですが、実際の地震は数十秒の繰り返し動荷重で、累積変位・液状化・残留強度の低下などが効きます。1995年阪神・淡路、2011年東日本大震災、2018年北海道胆振東部地震では、擬似静的解析で Fs > 1 だった斜面でも、累積変位による破壊が観測されました。重要構造物近傍や大規模斜面では Newmark 法(変位ベース)や動的 FEM が併用されます。本ツールの k_h はあくまで「初期スクリーニング」と位置づけてください。

使い方ガイド

  1. 斜面高さ(m)、斜面角度(°)を入力。例:高さ15m、角度35°
  2. 土質定数を設定。粘着力c(kPa)と内部摩擦角φ(°)を入力。例:砂質土でc=5kPa、φ=32°
  3. 「計算実行」ボタンでBishop簡便法により安全率Fsを自動算出
  4. 降雨時の間隙水圧比r_uを増加させてケーススタディ実施。Fs=1となる臨界r_uを確認
  5. 抵抗モーメント・駆動モーメント比からすべりリスクを評価

具体的な計算例

火山灰質粘性土の盛土斜面:高さH=12m、斜面角β=32°、粘着力c=15kPa、内部摩擦角φ=28°、湿潤単位体積重量γt=18.5kN/m³を入力。Bishop簡便法による円弧すべり面の解析結果、安全率Fs=1.35が得られる。同一斜面で連続降雨によりr_u=0.35に達すると、Fs≒1.02となり限界状態に接近。このとき臨界斜面角は27°に低下する。抵抗モーメント850kN·m、駆動モーメント790kN·mの比率1.08で安定性の余裕が縮小することが判定される。

実務での注意点

  1. Bishop簡便法は円弧すべり面を仮定するため、層理面や断層に沿う直線すべりが支配的な場合はFellenius法や条分法の併用を推奨
  2. 間隙水圧比r_u=0.2~0.3で既成斜面、r_u≧0.4で大規模降雨時として実務で使い分け。瓦礫質斜面ではφ値の低下を考慮し5°程度割引くこと
  3. 地震時加速度0.2gを入力する場合、見かけ粘着力を10kPa低減する現地慣行あり。安全率1.2以上を一般的な最小許容値とする
  4. 間隙水位が斜面高さの50%以上に達する場合、簡便法の精度限界を認識し3次元有限要素法(FEM)での詳細検討を検討