• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to previous page

Development of multiscale simulator for revealing actual behavior of long-term permeability evolution after fracture generation within rocks

Research Project

Project/Area Number 23K22859
Project/Area Number (Other) 22H01589 (2022-2023)
Research Category

Grant-in-Aid for Scientific Research (B)

Allocation TypeMulti-year Fund (2024)
Single-year Grants (2022-2023)
Section一般
Review Section Basic Section 22030:Geotechnical engineering-related
Research InstitutionOsaka University

Principal Investigator

緒方 奨  大阪大学, 大学院工学研究科, 助教 (50868388)

Co-Investigator(Kenkyū-buntansha) 清水 麻由子  国立研究開発法人日本原子力研究開発機構, 核燃料・バックエンド研究開発部門 地層処分研究開発推進部, 研究職 (10751191)
岸田 潔  京都大学, 工学研究科, 教授 (20243066)
中島 伸一郎  山口大学, 大学院創成科学研究科, 准教授 (70346089)
安原 英明  京都大学, 工学研究科, 教授 (70432797)
福田 大祐  北海道大学, 工学研究院, 准教授 (80647181)
下茂 道人  公益財団法人深田地質研究所, その他部局等, 主席研究員 (90393694)
Project Period (FY) 2022-04-01 – 2026-03-31
Project Status Granted (Fiscal Year 2024)
Budget Amount *help
¥17,550,000 (Direct Cost: ¥13,500,000、Indirect Cost: ¥4,050,000)
Fiscal Year 2025: ¥3,250,000 (Direct Cost: ¥2,500,000、Indirect Cost: ¥750,000)
Fiscal Year 2024: ¥3,250,000 (Direct Cost: ¥2,500,000、Indirect Cost: ¥750,000)
Fiscal Year 2023: ¥3,120,000 (Direct Cost: ¥2,400,000、Indirect Cost: ¥720,000)
Fiscal Year 2022: ¥7,930,000 (Direct Cost: ¥6,100,000、Indirect Cost: ¥1,830,000)
Keywords岩盤内亀裂 / 透水特性 / 熱・水・力学・化学連成現象 / マルチスケール・マルチフィジックスシミュレータ / 岩盤内亀裂のミクロ構造 / ミクロスケールの岩石亀裂生成解析 / マルチGPU並列計算 / 岩石亀裂生成解析の高速化 / 鉱物分布モデル / 亀裂のミクロ構造 / 亀裂内マルチフィジックス現象 / 亀裂の透水性変化 / GPGPU並列計算 / 岩石鉱物の溶解・沈殿
Outline of Research at the Start

高レベル放射性廃棄物を岩盤内に隔離し、その長期閉じ込め性能を保証するためには、岩盤内亀裂の透水性を長期に渡り予察する必要がある。そのためには、地層処分時に岩盤中に亀裂が生成した後、熱・水・力学・化学連成現象の影響で亀裂の透水性が長期変化していくまでを予測評価する必要がある。しかし、従来の岩盤連成解析では、亀裂での連成現象に強く寄与する亀裂内ミクロ構造の影響を考慮できておらず、透水性変化の実像を捉えきれない。本研究では、亀裂生成から連成現象に伴う長期透水性変化まで、亀裂のミクロ構造特性を考慮した上で精密に再現するマルチスケール岩盤シミュレータを創成し、地層処分時の岩盤の透水特性変化の実像を導く。

Outline of Annual Research Achievements

本年度に取り組んだ各課題の成果概要を記す。
①ミクロ構造レベルで亀裂生成形態を予測する解析モデル構築に関する成果:昨年度に構築したミクロスケール亀裂生成解析モデルに岩石表面の化学元素分析のデータも取り入れることで、より正確な鉱物分布を反映可能となった。この改良版モデルを用いて花崗岩圧裂試験の解析を行った。その結果、実験より得られた亀裂の凹凸度(Z2、JRC)を解析結果は定量的によく再現しており、亀裂の造成形態をミクロ構造レベルで予測し得ることが示された。さらに、昨年度に亀裂生成解析に実装していたシングルGPUを用いたGPGPU並列計算をマルチGPUに対応可能な形へと拡張し、GPU枚数に概ね比例した計算高速化を達成した。この拡張により、スパコンでの超多数のGPUを駆使した大規模解析が実行可能となった。
②ミクロ構造に依存した亀裂内マルチフィジックス現象を記述するメゾスケール構成モデル構築に関する成果:昨年度構築した実験機構を用いて、地層処分時に想定される温度・拘束圧条件(拘束圧~3MPa、温度25度~60度)下での単一亀裂を有する花崗岩への透水試験及びX線CTによるその場観察を100日以上の長期スパンで実施した。一連の実験・観察より、亀裂全体としての透水性の経時低下挙動及びそれをもたらしている亀裂内の時々刻々の構造変化(開口部が時間と共に閉塞し接触していく)を定量評価することができた。これらの結果をふまえ、亀裂全体としての透水性の経時低下挙動を表現可能な構成モデルを確立できた。
③亀裂のミクロ構造を考慮可能なマルチスケール・マルチフィジックスシミュレータの構築に関する成果:構築予定のシミュレータのベースとなるメゾ~マクロスケールの岩盤連成解析コードを高度化した。具体的には、空洞掘削解析における破壊判定を有効応力ベースに変更し実際の空洞損傷領域の分布傾向をよく再現できるように改良した。

Current Status of Research Progress
Current Status of Research Progress

2: Research has progressed on the whole more than it was originally planned.

Reason

当初計画では本年度は、(1)「ミクロ構造レベルで亀裂生成形態を予測する解析モデルの構築」、(2)「ミクロ構造に依存した亀裂内マルチフィジックス現象を記述する亀裂の構成モデルの構築」、(3) 「亀裂のミクロ構造を考慮可能なマルチスケール・マルチフィジックスシミュレータの構築」という3課題に取り組む予定であった。
(1)の課題について当初計画通り、今年度終了の時点までに予定していた一連の内容を遂行することができ、亀裂造成形態をミクロ構造レベルで予測し得ることが示された。さらに、想定していなかった成果として、昨年度中に亀裂生成解析に搭載済みであったシングルGPU型の並列計算をマルチGPU型に拡張し更なる計算高速化を達成できている。亀裂を陽に表現可能な亀裂生成解析に対しマルチGPU型の並列計算まで搭載できた事例は、世界的に見ても極めて稀有である。
(2)の課題について当初計画では、今年度は(2-1)【温度・圧力・化学条件を制御した岩石亀裂透水試験とその場観察】のサブ課題を長期に渡り継続実施するとともに、得られた計測・観察結果を活用し、(2-2)【ミクロ構造に依存した亀裂内マルチフィジックス現象を記述する亀裂の構成モデルの定式化】のサブ課題を遂行する予定であった。結果として今年度は、(2-1)の実験を100日以上長期に渡り実施するとともに、実験結果の傾向を表現可能な(2-2)の構成モデルを提案することができている。
(3)の課題について当初計画では、今年度と次年度でシミュレータを構築予定であり、今年度は構築予定のシミュレータのベースとなるメゾ~マクロスケールの岩盤連成解析コードの整備・性能検証まで実施できた。
上記の通り、(1)と(2)の研究課題で順調に成果があがっており、(3)の課題についても想定の範囲内で進捗しており、総合的には研究は概ね順調に進捗していると判断した。

Strategy for Future Research Activity

本研究で実施予定としている(1)~(4)の研究課題の内、(1) 「ミクロ構造レベルで亀裂生成形態を予測する解析モデルの構築」の課題が基本的に完遂できたので、今後は残りの課題に注力していく。(2)「ミクロ構造に依存した亀裂内マルチフィジックス現象を記述する亀裂の構成モデルの構築」の課題については、岩石亀裂透水試験とその場観察を引き続き実施し実験デーを拡充していき、今年度提案した構成モデルの更なる高精度化を図っていく予定である。また、(3) 「亀裂のミクロ構造を考慮可能なマルチスケール・マルチフィジックスシミュレータの構築」の課題については、今年度整備が完了した連成解析コードに(1)と(2)の課題で得られたミクロスケールモデルとメゾスケールモデルを次年度中に組み込み、ミクロ~メゾ~マクロを繋ぎ合わせたマルチスケール解析へと昇華させていきたい。また、JAEAの瑞浪超深地層研究所・幌延深地層研究センターで行われた現場実験の再現解析を行い、シミュレータの実現場レベルの妥当性を検証していく予定である。最後に、(4) 「亀裂のミクロ構造を考慮可能なマルチスケール・マルチフィジックスシミュレータの構築」の課題については、(3)にて構築・検証されたシミュレータを用いて地層処分環境を想定した条件下で岩盤の長期透水特性を数万年オーダーで予測解析を行っていく予定である。

Report

(2 results)
  • 2023 Annual Research Report
  • 2022 Annual Research Report
  • Research Products

    (42 results)

All 2024 2023 2022 Other

All Journal Article (10 results) (of which Int'l Joint Research: 2 results,  Peer Reviewed: 8 results,  Open Access: 2 results) Presentation (31 results) (of which Invited: 1 results) Remarks (1 results)

  • [Journal Article] Modeling Rock Permeability Alteration in NaCl Solution Flow-through Experiments2024

    • Author(s)
      Achmad HAFIDZ, Sho OGATA, Hideaki YASUHARA, Naoki KINOSHITA, Yuri IWABUCHI
    • Journal Title

      International Journal of the JSRM

      Volume: 20 Issue: 2 Pages: 1-5

    • DOI

      10.11187/ijjsrm.240202

    • ISSN
      2189-8405
    • Year and Date
      2024-01-05
    • Related Report
      2023 Annual Research Report
    • Peer Reviewed
  • [Journal Article] Laboratory observation of shear stress recovery in rock joints considering joint compaction and dilatancy2024

    • Author(s)
      Zhang Jintong、Kikumoto Mamoru、Yasuhara Hideaki、Ogata Sho、Kishida Kiyoshi
    • Journal Title

      International Journal of Rock Mechanics and Mining Sciences

      Volume: 174 Pages: 105666-105666

    • DOI

      10.1016/j.ijrmms.2024.105666

    • Related Report
      2023 Annual Research Report
    • Peer Reviewed / Int'l Joint Research
  • [Journal Article] Development of a GPGPU-parallelized FDEM based on extrinsic cohesive zone model with master-slave algorithm2024

    • Author(s)
      Maeda Yutaro、Ogata Sho、Fukuda Daisuke、Liu Hongyuan、Inui Toru
    • Journal Title

      Computers and Geotechnics

      Volume: 166 Pages: 105942-105942

    • DOI

      10.1016/j.compgeo.2023.105942

    • Related Report
      2023 Annual Research Report
    • Peer Reviewed / Int'l Joint Research
  • [Journal Article] A Framework of Coupled Thermal-Hydraulic-Mechanical Analysis for Simulating Hydraulic Fracturing Based on Explicit and Implicit Methods2023

    • Author(s)
      Maeda Y.、Ogata S.、Fukuda D.、Inui T.、Yasuhara H.、Kishida K.
    • Journal Title

      The 57th U.S. Rock Mechanics/Geomechanics Symposium

      Volume: 1

    • DOI

      10.56952/arma-2023-0759

    • Related Report
      2023 Annual Research Report
    • Peer Reviewed
  • [Journal Article] Numerical simulations for describing generation of excavation damaged zone: Important case study at Horonobe underground research laboratory2023

    • Author(s)
      Ogata Sho、Yasuhara Hideaki
    • Journal Title

      Rock Mechanics Bulletin

      Volume: 2 Issue: 3 Pages: 100063-100063

    • DOI

      10.1016/j.rockmb.2023.100063

    • Related Report
      2023 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] Multi-physics numerical analyses for predicting the alterations in permeability and reactive transport behavior within single rock fractures depending on temperature, stress, and fluid pH conditions2022

    • Author(s)
      Ogata Sho、Nishira Eita、Yasuhara Hideaki、Kinoshita Naoki、Inui Toru、Kishida Kiyoshi
    • Journal Title

      Soils and Foundations

      Volume: 62 Issue: 6 Pages: 101207-101207

    • DOI

      10.1016/j.sandf.2022.101207

    • Related Report
      2022 Annual Research Report
    • Peer Reviewed / Open Access
  • [Journal Article] Development of hydraulic fracturing simulator based on explicit and implicit methods2022

    • Author(s)
      Yutaro Maeda, Sho Ogata, Daisuke Fukuda, Hideaki Yasuhara, Toru Inui
    • Journal Title

      The proceedings of CouFrac2022

      Volume: No. 165 Pages: 1-4

    • Related Report
      2022 Annual Research Report
    • Peer Reviewed
  • [Journal Article] Coupled THMC Simulation Based on Explicit Fracture Representation Using Extrinsic Cohesive Zone Model2022

    • Author(s)
      Ogata S.、Maeda Y.、Fukuda D.、Yasuhara H.、Inui T.、Kishida K.
    • Journal Title

      Proceedings of 56TH US ROCK MECHANICS / GEOMECHANICS SYMPOSIUM

      Volume: ARMA-2022-0642

    • DOI

      10.56952/arma-2022-0642

    • Related Report
      2022 Annual Research Report
    • Peer Reviewed
  • [Journal Article] Numerical analyses of coupled thermal-hydraulic-mechanical-chemical processes for estimating permeability change in fractured rock induced by alkaline solution2022

    • Author(s)
      Ogata Sho、Yasuhara Hideaki、Kinoshita Naoki、Inui Toru、Nishira Eita、Kishida Kiyoshi
    • Journal Title

      Geomechanics for Energy and the Environment

      Volume: 31 Pages: 100372-100372

    • DOI

      10.1016/j.gete.2022.100372

    • Related Report
      2022 Annual Research Report
  • [Journal Article] Modeling the shearing behavior of discontinuous rock mass incorporating dilation of joint aperture2022

    • Author(s)
      Zhang Jintong、Kikumoto Mamoru、Yasuhara Hideaki、Ogata Sho、Kishida Kiyoshi
    • Journal Title

      International Journal of Rock Mechanics and Mining Sciences

      Volume: 153 Pages: 105101-105101

    • DOI

      10.1016/j.ijrmms.2022.105101

    • Related Report
      2022 Annual Research Report
  • [Presentation] μX線CTと機械学習を組み合わせた 花崗岩の鉱物分布評価手法の開発とGrain-based modelへの応用2023

    • Author(s)
      尾崎 友星, 緒方 奨, 下茂 道人, 木下 尚樹, 東 辰之介, 乾 徹
    • Organizer
      第 50 回岩盤力学に関するシンポジウム
    • Related Report
      2023 Annual Research Report
  • [Presentation] 岩石のCO2水押破砕シミュレータの開発2023

    • Author(s)
      前田 悠太朗, 緒方 奨, 詫間 康平, 渡邉 則昭, 大里 和己, 寺井 周
    • Organizer
      第 50 回岩盤力学に関するシンポジウム
    • Related Report
      2023 Annual Research Report
  • [Presentation] 亜臨界領域での地化学反応に伴う岩石亀裂の透水性変化を考慮した透水-反応輸送連成 シミュレータの開発と検証2023

    • Author(s)
      永瀬 みちる, 緒方 奨, 岡本 敦, 乾 徹
    • Organizer
      第 50 回岩盤力学に関するシンポジウム
    • Related Report
      2023 Annual Research Report
  • [Presentation] Observation of Fracture Permeability by Flow-through Tests and X-ray CT Tests2023

    • Author(s)
      Zhiqi LI, Sho OGATA, Shinichiro NAKASHIMA, Hideaki YASUHARA, Kiyoshi KISHIDA
    • Organizer
      第 50 回岩盤力学に関するシンポジウム
    • Related Report
      2023 Annual Research Report
  • [Presentation] Phase-Field 法に基づく低粘性を考慮した SC-CO2 破砕シミュレータの開発と検証2023

    • Author(s)
      谷澤 侑樹, 緒方 奨, 乾 徹
    • Organizer
      第 50 回岩盤力学に関するシンポジウム
    • Related Report
      2023 Annual Research Report
  • [Presentation] 花崗岩不連続面の塩水浸透に関する室内実験2023

    • Author(s)
      岩渕 友梨, 木下 尚樹, 安原 英明, 緒方 奨, Achmad Hafidz
    • Organizer
      第 50 回岩盤力学に関するシンポジウム
    • Related Report
      2023 Annual Research Report
  • [Presentation] 亀裂発生・進展を考慮した 3次元坑道掘削シミュレーション2023

    • Author(s)
      末武 航弥, 緒方奨, 安原 英明, 岸田 潔
    • Organizer
      第 50 回岩盤力学に関するシンポジウム
    • Related Report
      2023 Annual Research Report
  • [Presentation] Development of Coupled THMC Model for Predicting Rock Permeability Evolution Under a Seawater Condition2023

    • Author(s)
      Achmad HAFIDZ, Sho OGATA, Naoki KINOSHITA, Hideaki YASUHARA
    • Organizer
      第 50 回岩盤力学に関するシンポジウム
    • Related Report
      2023 Annual Research Report
  • [Presentation] Elastoplastic Model for Rock Joint Considering Matching and Adhesion2023

    • Author(s)
      Y. Matsuoka, M. Kikumoto, S. Ogata, K. Kishida
    • Organizer
      The 34th KKHTCNN Symposium on Civil Engineering
    • Related Report
      2023 Annual Research Report
  • [Presentation] Short-term and long-term observation of fractured granite permeability by flow-through tests and X-ray CT tests2023

    • Author(s)
      LI. Zhiqi, OGATA. Sho, NAKASHIMA. Shinichiro, YASUHARA Hideaki, KISHIDA. Kiyoshi
    • Organizer
      The 34th KKHTCNN Symposium on Civil Engineering
    • Related Report
      2023 Annual Research Report
  • [Presentation] 地熱貯留層条件下でのCO2圧入による火山岩の破砕メカニズム2023

    • Author(s)
      詫間 康平, 渡辺 優斗, 坂口 清敏, 渡邉 則昭, 前田 悠太朗, 緒方 奨, 大里和己, 寺井 周
    • Organizer
      日本地熱学会令和五年学術講演会(岐阜大会)
    • Related Report
      2023 Annual Research Report
  • [Presentation] 火山岩地熱貯留層のCO2フラクチャリングシミュレーション2023

    • Author(s)
      詫間康平, 渡辺優斗, 坂口清敏, 渡邉則昭, 緒方奨, 前田悠太朗, 大里和己, 寺井周
    • Organizer
      資源・素材学会東北支部令和五年春季大会
    • Related Report
      2023 Annual Research Report
  • [Presentation] 亜臨界条件における花崗岩亀裂の透水性経時変化予測解析2023

    • Author(s)
      永瀬みちる、緒方奨、乾徹、岡本敦
    • Organizer
      日本材料学会第72期学術講演会
    • Related Report
      2023 Annual Research Report
  • [Presentation] 三次元T-H-M-C連成数値シミュレータによる天然バリアの長期透水特性評価2023

    • Author(s)
      末武航弥、緒方奨、乾徹
    • Organizer
      日本材料学会第72期学術講演会
    • Related Report
      2023 Annual Research Report
  • [Presentation] X-FEM亀裂進展解析とFEM流体浸透解析を連成した水圧破砕解析の提案2023

    • Author(s)
      谷澤侑樹、緒方奨、乾徹
    • Organizer
      日本材料学会第72期学術講演会
    • Related Report
      2023 Annual Research Report
  • [Presentation] 亀裂表面形状を考慮した透水ー反応輸送解析による模擬海水中での岩石亀裂の透水性変化予測2022

    • Author(s)
      西羅瑛太, 緒方奨, 安原英明, 木下尚樹, 乾徹, 岸田潔
    • Organizer
      第49回岩盤力学に関するシンポジウム
    • Related Report
      2022 Annual Research Report
  • [Presentation] 垂直拘束圧条件下における岩石亀裂面挙動評価手法の検討2022

    • Author(s)
      尾崎友星, 緒方奨, 安原英明, 木下尚樹, 佐古大地, 乾徹
    • Organizer
      第49回岩盤力学に関するシンポジウム
    • Related Report
      2022 Annual Research Report
  • [Presentation] 開口幅と固着の影響を考慮した岩石不連続面せん断モデルー定式化と適用性2022

    • Author(s)
      松岡勇樹, 菊本統, 緒方奨, 岸田潔
    • Organizer
      第49回岩盤力学に関するシンポジウム
    • Related Report
      2022 Annual Research Report
  • [Presentation] Extrinsic Cohesive Zone ModelをベースとしたハイブリッドFEM-DEM解析の並列化2022

    • Author(s)
      前田悠太朗, 緒方奨, 福田大祐, 乾徹
    • Organizer
      第49回岩盤力学に関するシンポジウム
    • Related Report
      2022 Annual Research Report
  • [Presentation] 熱―流体ー力学連成モデルによる三次元地熱流体流動解析2022

    • Author(s)
      前原崇志, 緒方奨, 福田大祐, 乾徹
    • Organizer
      第49回岩盤力学に関するシンポジウム
    • Related Report
      2022 Annual Research Report
  • [Presentation] 稲田花崗岩の室内水圧破砕実験と鉱物分布を考慮した解析モデルの検討2022

    • Author(s)
      佐古大地, 緒方奨, 安原英明, 木下尚樹
    • Organizer
      第49回岩盤力学に関するシンポジウム
    • Related Report
      2022 Annual Research Report
  • [Presentation] 地熱資源開発やエネルギー生成後の副産物地下貯留に関する熱-水理-力学-化学連成数値 シミュレーション2022

    • Author(s)
      緒方奨
    • Organizer
      第6回 COMSOLを用いたカーボンニュートラルを支援するためのCAEセミナー
    • Related Report
      2022 Annual Research Report
    • Invited
  • [Presentation] 陽解法-陰解法連成解析法を用いた水圧破砕解析2022

    • Author(s)
      緒方奨, 前田悠太朗, 福田大祐, 陳友晴
    • Organizer
      日本地熱学会令和4年学術講演会
    • Related Report
      2022 Annual Research Report
  • [Presentation] 亀裂性岩盤における熱-水理-力学-化学連成場を解くマルチフィジックスシミュレータの開発2022

    • Author(s)
      緒方奨
    • Organizer
      日本機械学会M&Mカンファレンス
    • Related Report
      2022 Annual Research Report
  • [Presentation] 損傷理論を用いた大土被り条件下の掘削解析における側圧係数の影響2022

    • Author(s)
      三好 航平, 緒方 奨, 安原 英明, 岸田 潔
    • Organizer
      第57回地盤工学発表会
    • Related Report
      2022 Annual Research Report
  • [Presentation] 岩石不連続面のせん断強度回復のモデル化とパラメータの検討2022

    • Author(s)
      松岡勇樹, 菊本統, 緒方奨, 張晋通, 岸田潔
    • Organizer
      第57回地盤工学発表会
    • Related Report
      2022 Annual Research Report
  • [Presentation] 垂直拘束圧条件下における岩石亀裂構造評価手法の検討2022

    • Author(s)
      尾崎友星, 緒方奨, 安原英明, 木下尚樹, 佐古大地, 乾徹
    • Organizer
      2022年度土木学会関西支部年次学術講演会
    • Related Report
      2022 Annual Research Report
  • [Presentation] 陽解法と陰解法を連成した岩石の水圧破砕シミュレータの開発2022

    • Author(s)
      家永凌冴, 緒方奨, 前田悠太朗, 乾徹
    • Organizer
      2022年度土木学会関西支部年次学術講演会
    • Related Report
      2022 Annual Research Report
  • [Presentation] 割れ目の開口幅を考慮した岩石不連続面のせん断シミュレーションとパラメータの分析2022

    • Author(s)
      松岡勇輝, 菊本統, 緒方奨, 岸田潔
    • Organizer
      2022年度土木学会関西支部年次学術講演会
    • Related Report
      2022 Annual Research Report
  • [Presentation] 陽解法-陰解法連結計算による水圧破砕解析手法の開発2022

    • Author(s)
      家永凌冴,緒方 奨,前田悠太朗,福田大祐,乾 徹
    • Organizer
      日本材料学会第71期学術講演会
    • Related Report
      2022 Annual Research Report
  • [Presentation] 垂直拘束圧条件下での花崗岩亀裂の開口幅・接触状態評価2022

    • Author(s)
      尾崎友星,緒方 奨,安原英明,木下尚樹,佐古大地,乾 徹
    • Organizer
      日本材料学会第71期学術講演会
    • Related Report
      2022 Annual Research Report
  • [Remarks] 掘削損傷領域の予測解析の成果が、EurekAlert!にて英文プレスリリース

    • URL

      https://www.eurekalert.org/news-releases/1002125

    • Related Report
      2023 Annual Research Report

URL: 

Published: 2022-04-19   Modified: 2024-12-25  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi