Question

矿物化学式与端员组分计算

Answer 1

1.程序功能

程序MIFORM的功能包括：①常见造岩矿物晶体化学式计算；②硫化物族矿物晶体化学式计算；③常见铁镁矿物电子探针分析结果的Fe³⁺、Fe²⁺含量计算；④常见造岩矿物端员组分计算；⑤矿物摩尔质量的计算；⑥辉石族矿物八面体阳离子占位的计算；⑦辉石族矿物晶胞体积、密度和平均折射率［n＝（N_g+n_m+n_p）/3］计算；⑧矿物的比折射能（即格拉斯顿-代尔常数）计算。根据格拉斯顿-代尔定律，由矿物的比折射能和密度，可进一步求出其平均折射率（叶大年，1988）。

2.方法原理

本程序中计算矿物晶体化学式的方法采用Jackson等（1976）提出的以阴离子为基准计算矿物化学式的氢当量法（马鸿文，1993a）。硫化物晶体化学式的计算则采用以阴离子为基准的硫原子法（含类质同像替代元素Se、Te、As、Sb）。

常见铁镁矿物中Fe³⁺、Fe²⁺含量的计算采用电价平衡原理（Droop,1987）。其中闪石族矿物中的Fe³⁺、Fe²⁺含量计算采用Holland和Blundy（1994）的方法。

常见造岩矿物端员组分的计算见马鸿文（1993a）归纳的方法。辉石族矿物八面体阳离子占位的计算据Brey等（1990）。

辉石族的晶胞体积、密度和平均折射率计算据马鸿文等（1998，未发表资料）。矿物比折射能和平均折射率的计算据叶大年（1988）。

3.程序结构

程序框图见下页。

4.使用说明

（1）输入格式

程序运行过程中，按照屏幕提示，依次提供下列参数：

Minerl　矿物族名

IFN/OFN　输入/输出文件名（选择硫化物时IFN缺省）

JFe3　Fe³⁺/Fe²⁺标定选择

Totl　分析总量下限

Totu　分析总量上限

输入文件格式：A7,13F6.2。

各变量排列顺序依次为：Sample（样品号）、SiO₂、TiO₂、Al₂O₃、Cr₂O₃、Fe₂O₃、FeO、MnO、NiO、MgO、CaO、Na₂O、K₂O、Li₂O。

每次计算的样品个数和矿物数不限。

若实际应用中需增加其他变量，只需改变parameter语句中的参数m和data语句中的OXT参数，不影响程序的正常运行。

（2）输出格式

全部计算结果输出到文件OFN中。输出结果除矿物阳离子系数和端员组分外，还包括下列内容：

GFW　摩尔质量（即分子量）

REK　矿物比折射能K

Vcell　晶胞体积（单斜辉石）或其值的1/2（斜方辉石）

DST　密度（辉石族）

RI　平均折射率（辉石族）

结晶岩热力学软件

5.程序文本

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

结晶岩热力学软件

6.计算实例

新疆西准噶尔唐巴勒变质橄榄岩中单斜辉石的晶体化学式计算（马鸿文，1990，博士后基金和国家305项目研究报告）。

输入文件：examll.dat

结晶岩热力学软件

结晶岩热力学软件

输出文件：exam12.dat

结晶岩热力学软件

结晶岩热力学软件