3-448

土木学会第63回年次学術講演会(平成20年9月)

メタンハイドレート生産時の地盤変形に関する室内模型実験
エンジニアリング振興協会（清水建設） フェロー会員 ○荻迫栄治

１．はじめに同上

正会員

西尾伸也

正会員

傳田 篤

京都大学フェロー会員

岡二三生

正会員

木元小百合

エンジニアリング振興協会（伊藤忠テクノソリューションズ）

正会員

岡部直司

メタンハイドレートは次世代資源として注目されており、日本周辺海底にも日本が消費している天然ガスの約 100 年分の量が存在すると推定されている。経済産業省が策定した「我が国におけるメタンハイドレート開発計画」に従い、「メタンハイドレート資源開発研究コンソーシアム」が組織され、（財）エンジニアリング振興協会が環境影響評価に関する研究開発を行っており、筆者らは、この中で、メタンハイドレートからのメタンガス生産に伴う海底地盤の変形の可能性について検討を進めている。今回、室内模型実験によりメタンハイドレート生産時の地盤変形挙動について検討を行うとともに、実験を対象に弾粘塑性構成式を用いた数値解析を行った。
２．実験概要実験は、できるだけ実際の海底地盤のスケールに近づけるために遠心載荷装置で行った。

遠心載荷装置に設置する土槽の寸法は、内寸で高さ 750mm・幅 800mm・奥行き 300mm である。土槽前面は変形が観察可能なようにアクリル版とした。土槽内にメタンハイドレートの生産を模擬するためのラバーバルーンを設置し、バルーンの圧力を下げることで生産を模擬する。土槽の前面および背面にはシリコンシートを貼り、その上からシリコンオイルを塗布し、ラバーメンブレンを貼付して壁面摩擦の軽減を図った。地
表面の変位は変位計により、また、地盤内の変位は上記ラバーメンブレンに記した標点の変位を前方からデジタルカメラで撮影することにより、地盤の変形を観察した。図-1 に実験土槽の概要を示す。実験は以下のような手順で行った。(1)ラバーバルーン下面の高さまで土槽に試料土をセットする。(2)土槽にラバーバルーンを設置し、バルーン内を脱気水で満たし、給水用バルブを閉じる。バルーンの初期体積は 2300cm3 である。
(3)所定の高さまで土槽に試料土をセットする。(4)飽和土の場合は、給水タンクを用い土槽底面の給水用穴から通水し、飽和させる。(5)地表面に計測用の変位計をセットする。(6)

変位計

土槽を遠心載荷装置にセットし、遠心加速度を作用させる。(7)コント
ロールシリンダーにてバルーンの体積を減少させることによってメタ50cm）、飽和砂（地盤厚さ 60cm）、飽和砂（地盤厚さ 50cm）の 4 ケー

750
100は 30G である。実験は、乾燥砂（地盤厚さ 60cm）、乾燥砂（地盤厚さ200

落下法により乾燥密度約 1.5g/cm3 の模型地盤を作成した。遠心加速度

600または500ンハイドレート生産を模擬する。実験に用いた土は 7 号硅砂で、空中ス行った。
３．実験結果図-2 は、乾燥砂（地盤厚さ 60cm）のケースについ800

ラバーバルーン
外径100mm

土側面図て、バルーン体積減少開始時点を時刻０としたときのバルーン体積変ラバーバルーンある。バルーン体積変化量はコントロールシリンダーのストロークの300化量ならびに変位計で計測した地表面変位の経時変化を示したもので変位量より算出した。変位 1 は土槽中央部の変位、変位 2〜5 はそれぞ
100

れ土槽中央から 80, 160, 240, 320mm 離れた位置における変位である。
地表面変位はバルーン体積の減少につれて増加するが、その増加傾向は一定ではない。すなわち、当初は緩やかに増加するが、時刻 10min

800
平面図

単位：mm

図-1 実験土槽の概要

キーワード：メタンハイドレート，海底地盤，変形
連絡先：〒135-8530 東京都江東区越中島 3-4-17 清水建設(株)技術研究所 Tel:03-3820-6476 Fax:03-3820-5955-895-土木学会第63回年次学術講演会(平成20年9月)

00‑500‑5

地表面変位(mm)

バルーン体積変化量(cm3)3-448‑1000‑1500‑2000‑2500
051001520時間(min)2530

バルーン体積変化量‑10‑400cm3‑800cm3‑1000cm3‑1200cm3‑1400cm3‑1800cm3‑15‑2035‑25‑40‑30‑20‑10010203040地表面変位(mm)

土槽中心からの距離(cm)

図-3 地表面変位分布‑10変位1変位2変位3変位4変位5‑20‑3005101520253035時間(min)

図-2 バルーン体積変化量および地表面変位の経時変化付近から増加の割合が大きくなっている。そして、バルーン体積の減少が止まって以降はほぼ一定の値

写真-1 地盤の変形状況（バルーン体積変化-1800cm3 時）に収束している。図-3 は、バルーン体積変化量ごとにプロットした水平方向の地表面変位分布である。地表面変位はバルーン直上の土槽中央部で最も大きく、中央から離れるにつれて次第に小さくなっている。バルーンの体積減少量が大きくなるにつれて土槽中央部の変位の進展が大きくなっていくことがわかる。写真-1 は、バルーン体積変化-1800cm3、ほぼ 23 分後の土槽前面から撮影した変形状況の写真である。標点の動きからバルーン直上の地盤が変形している様子がうかがえる。
４．実験の数値解析

上記の室内模型実験を対象に有限要素法による数値解析を実施した。解析手法としては 2 次元平面ひずみ弾粘塑性解析である。なお、乾燥砂のケースのため間隙水はないものとした。土の構成式としては岡ら 1) の提案する凍結砂のひずみ軟化型弾粘塑性構成式を用いた。
なお、解析に際しては弾性係数の拘束圧依存性とひずみ依存性を考慮し、{ ( )}
β′ /σ m′ 0 ， G / G0 = 1 / 1 + αγ p
E / E0 = σ m′ ，σ m′ 0 ：平均有効応力， G ，G ：せん断
ここで、 E ，E0 ：ヤング係数，σ m0
次式を仮定した。

弾性係数， γ p = ∫α ，β ：定数で添え字

deijp deijp ：蓄積塑性せん断ひずみ，の 0 は初期値を示す。解析手順としては、まず、バルーン部の節点に実験

表-1 解析に用いた物性パラメータ
ヤング係数 E (kPa)21200

体積弾性係数 K (kPa)23600

ひずみ硬化－軟化パラメータ G'8.50

ひずみ硬化－軟化パラメータ Mf*

1.22

塑性ポテンシャルパラメータ σmb (kPa)90塑性ポテンシャルパラメータ b (kPa)
過圧密境界面パラメータ Mm*

01.03

応力履歴パラメータ τ0 (sec)3000

応力履歴パラメータ a0.97

ひずみ依存パラメータ α100

ひずみ依存パラメータ β4.0(m)0.
-0.1-0.2時の初期水圧値に相当する圧力を載荷した状態で土の自重解析を行い、次-0.3-0.4に、実験時のバルーンの体積変化量から算出した変位量をバルーン部の節-0.5

点に簡単のため鉛直方向の強制変位として作用させた。解析に用いた物性-0.6

パラメータは、実験に用いた試料の三軸圧縮試験結果を基に表-1 のように

-0.8-0.7-0.9設定した。図-4 に、時刻 500min における鉛直変位コンター図の解析結果-1.
-1.1を示す。バルーン直上部の地盤の鉛直変位が増加している。しかし、鉛直-1.2変位の発生する領域はあまり大きくなく、バルーン直上付近の地盤にある-1.3程度限られており、これは写真-1 に示した実験結果の傾向と一致している。

-1.5-1.4-1.6５．おわりに室内模型実験によりメタンハイドレート生産時の地盤変

図-4 鉛直変位コンター図（500min）

形挙動について検討を行った。その結果、生産を模擬したバルーンの直上付近の地盤に主要な変形が発生していることが確認された。また、数値解析の結果も同様の傾向を示していた。
参考文献

1) 足立，岡：凍結砂のひずみ軟化型弾粘塑性構成式，土木学会論文集，No.454/III-20，p.75-81, 1992.

-896-