✅ Phreeqcマニュアル

✅ Water-Calcite-CO2反応 

Reference: Appelo  C.A.J. and Postma D. (2005) GEOCHEMISTRY, GROUNDWATER AND POLLUTION. A.A. Balkema Publishers, Leiden, The Netherlands a member of Taylor & Francis Group plc  pp.634

✔ PHREEQCで、Open system (Solution 1) と Closed system (Solution 2)のコードを作成し、pH変化を確認する。

SOLUTION 1　（Open system）

    temp      25

    pH        7

    pe        4

    redox     pe

    units     mmol/kgw

    density   1

    -water    1 # kg

EQUILIBRIUM_PHASES 1

    CO2(g)    -1.5 10

    Calcite   0 10

END  ＃ENDが１つで上記のOpen systemの化学式のCO2とCalciteが一緒に反応することを表現

SOLUTION 2　（Closed system）

    temp      25

    pH        7

    pe        4

    redox     pe

    units     mmol/kgw

    density   1

    -water    1 # kg

EQUILIBRIUM_PHASES 1

    CO2(g)    -1.5 10

SAVE solution 2

END ＃はじめのENDで、上記のClosed systemの化学式のCO2が水と反応することを表現

USE solution 2 ＃CO2と水が反応した溶液を用いることを表現

EQUILIBRIUM_PHASES 2 ＃CO2と水が反応した溶液がCalciteと反応することを表現

    Calcite   0 10

END ＃simulation終了

☑ 結果

✅ Phreeqcモデリング（CO2ガスなし（C1）CO2ガスあり（C2）における重炭酸イオン濃度の計算） 

Reference: Yang H., Mishima T., Katazakai S., Kagabu M (2023). Analytical approach using a chemical equilibrium formula and geochemical modeling for alkalinity measurements of small natural water samples. Applied Geochemistry    https://doi.org/10.1016/j.apgeochem.2022.105535 

✅ Phreeqcモデリング（反応速度解析（玄武岩-CO2-H2Oバッチ試験）） 

1.　圧力10MPa、温度250℃、CO2分圧4MPaにおける炭酸水を用意する。下記の図は、実際行っている実験をモデル化したプログラム

２.　玄武岩の鉱物（斜長石、単斜輝石、橄欖石）、Kinetics、Rate設定（ただし、Rateは長いためにAlbiteのみを表示）