日本地球惑星科学連合2015年大会

講演情報

口頭発表

セッション記号 H (地球人間圏科学) » H-RE 応用地質学・資源エネルギー利用

[H-RE28] 地球温暖化防止と地学(CO2地中貯留・有効利用,地球工学)

2015年5月25日(月) 11:00 〜 12:45 105 (1F)

コンビーナ:*徳永 朋祥(東京大学大学院新領域創成科学研究科環境システム学専攻)、薛 自求(財団法人 地球環境産業技術研究機構)、徂徠 正夫(独立行政法人産業技術総合研究所地圏資源環境研究部門)、座長:徂徠 正夫(独立行政法人産業技術総合研究所地圏資源環境研究部門)

11:15 〜 11:30

[HRE28-08] CO2分布に対する地層不均質性の効果:長岡サイトCO2貯留層中の数値シミュレーション

*中島 崇裕1伊藤 拓馬1薛 自求1千代延 俊2 (1.(公財)地球環境産業技術研究機構、2.秋田大学国際資源学部)

キーワード:CO2地中貯留, 長岡サイト, 多相流, 不均質

This paper discusses heterogeneity of petrophysical properties in the reservoir for geological carbon dioxide (CO2) storage at Nagaoka pilot site in Japan. Detailed site characterization is critical for successful geological storage of CO2. Not only the traditional site characterization techniques, but also the injection and monitoring data can be used for the geological modeling with heterogeneity. In this study, we carried out numerical simulations using reservoir models with heterogeneity, and compared the results to investigate effects of heterogeneity in the reservoir.
Nagaoka pilot-scale CO2 injection site was used in this study. Nagaoka project was undertaken in order to verify an ability of CO2 injection into Japanese formation. CO2 was injected into a thin permeable zone at a depth of 1100m, and the total amount of injected CO2 was 10.4 k-tons. Bottom-hole pressure, time-lapse well loggings and cross-well seismic tomography were conducted as the monitoring of injected CO2.
For the modeling of the reservoir at this site, the method presented by Ito et al. (2015) was used; constructed a profile of the reservoir using seismic results, and generated horizontal distribution of petropysical properties by applying Random function Gaussian Simulation. It is worth to note that the feature along N-S direction had reported at this site (Chiyonobu et al., 2013). The derived models showed that the high porosity and permeability area located in the eastern part of the reservoir.
For the reservoir simulation, hydrological properties for the multiphase flow are also important. We measured capillary pressure function using mercury injection tests. The measurements indicated that the functions in sandy and shaly rocks were different. For the relative permeability function, we referred the results reported by Ohtake (2013), and approximated to van Genchten function. We did not assume hysteretic functions in this study.
We made simulation using the reconstructed models and TOUGH2 simulator. For the history matching, bottom-hole pressure and CO2 saturation data were used. Several absolute permeability models were tested manually in order to find the best match between the monitoring data and the reservoir simulation. We could find a geological model reasonably matched to the monitoring data. The result of CO2 distribution was also consistent with the observed velocity anomalies by the cross-well tomography. The numerical results revealed the migration of CO2 plume to up-dip direction. This explained the behavior of CO2 distribution observed by the wire-line logging very well.
During the history matching, sensitivity analysis of the model was conducted. The results suggested that anisotropy of the permeability was essential to explain the monitored CO2 behaviors. This anisotropy could be created during the depositional process of the reservoir. It should be noted that not only the permeability distribution, but the properties depending on facies were also important to the matching.

This work was supported by Ministry of Economy, Trade and Industry of Japan under the research contract "Development of Safety Assessment Technology for Carbon Dioxide Capture and Storage".