地球化学图相关计算
地球化学 (14)
REE的较大半径削弱着共价键性和静电的相互作 用， 成为阻止REE形成稳定络合物的主要因素之 一。溶液中三价REE离子能同CO23-，Br+，I-， NO3-和SO42-等组成离子对，形成碳酸盐，硫酸盐， 氯化物和氟化物型络合物. ；
在富CO2溶液中REE极活动。 实验证明HREE在共 存硅酸盐和碳酸盐熔体(岩浆熔离)之间优先富集于 碳酸盐熔体中；在共存富CO2蒸气相中REE更加富 集；
尽管REE化学行为相似， 还是能通过某些成岩和 成矿过程发生彼此分离。这是因电子构型对它们 离子价态和半径施加影响的结果，也与REE在造 岩矿物配位多面体类型众多和大小变化有关。
2. REE价态
REE是强正电性元素， 以离子键为特征， 只含极 少共价成分。电离顺序是先移去6s亚层上两个电 子，然后丢失1个5d或4f电子，因为5d和4f电子在 能量上相对接近于1个6s电子。设想再从4f移去1 个，即第4个电子，但是这个电子电离能太高， 不易移去。因此REE在化学和地球化学上均显示 三价离子状态，只有Eu和Yb可呈2价， Ce和Tb可 呈4价。原因：Eu2＋和Tb4＋具有半充满4f亚层， Yb2＋具有全充满4f亚层，Ce4＋具有贵气体氙(Xe) 电子构型，这些电子构型可以提高该价态离子的 稳定性。
变价离子(Eu，Ce等)不同价态的比例取决于 体系的成分、氧逸度、温度和压力；
3．REE的配位和离子半径
矿物中REE占据多种多样的配位多面体，从六次 到十二次，甚至更高的配位均有。较小的稀土元 素占据六次配位位置，但这种情况在矿物中少见。
一般REE在矿物中的配位要大些，最常见的配位 是七次到十二次，如榍石中为七次，锆石中为八 次，独居石中为九次，褐帘石中为十一次和钙钛 矿中为十二次。
岩石地球化学计算
岩石地球化学计算1. TFe2O3=FeO+0.9Fe2O3FeOT(wt.%)=FeO(wt.%)+Fe2O3(wt.%)*0.8998=FeO(wt.%)+Fe2O3(wt.%)*(71.844/(159.6882/2))2. LOI 烧失量3. Mg#=100*(MgO/40.3044)/(MgO/40.3044+FeOT/71.844)FeOm71.85 ；MgOm40.31上述是分别测试分析了FeO和Fe2O3的计算方法，如果是测试的全铁，也可以近似计算。
通常说的高Mg，是指岩石具有较高的MgO含量，如火山岩中的高镁安山岩（通常情况下，异常高的MgO含量指示着可能有地幔物质参与，如俯冲带地幔楔或者软流圈熔体上涌等等）。
Mg#（镁指数）也可以定量的表示岩石中的Mg含量高低。
Mg#通常用于镁铁质岩石，可以粗略指示地幔岩石的部分熔融程度，高Mg#的地幔橄榄岩可能经历了更高程度的部分熔融，常在92-93左右，而原始地幔会相对富集，Mg#较低，在88-89左右。
4. 里特曼组合指数δ或里特曼指数δ=（K2O+Na2O）2/（SiO2-43）（wt%）δ<3.3 者称为钙碱性岩，δ=3.3-9 者为碱性岩，δ>9 者为过碱性岩。
5.A/NK = Al2O3/102/(Na2O/62+K2O/94)6.A/CNK = Al2O3/102/(CaO/56+Na2O/62+K2O/94)7.全碱ALK = Na2O+K2O8.AKI = (Na2O/62+K2O/94)/Al2O3*1029.AR = (Al2O3+CaO+Na2O+K2O)/(Al2O3+CaO-Na2O-K2O)10.固结指数(SI) =MgO×100/(MgO+FeO+F2O3+Na2O+K2O) (Wt%)11.阳离子R1-R2图（岩石氧化物wt%总量不用换算成100%）R1=（4Si-11(Na+K)-2(Fe+Ti)*1000R2=（6Ca+2Mg+Al）*100012.(La/Sm)N对δEu的双变量斜边图解认识Eu异常。
地球化学
1.3.3 地壳化学成分特征和元素克拉克值的地球化学意义1.3.3.1 地壳元素丰度特征分析1) 地壳中元素相对的平均含量极不均匀。
按维氏（1949）值，丰度最大的元素（O=45.6%）比丰度最小的元素（Rn－氡,7×10-17%）在含量上大1017倍，相差十分悬殊。
按克拉克值递减的顺序排列，含量最多的前3种元素（O、Si、Al）即占地壳总重量的81.3％(图1.12)；含量最多的前9种元素（O、Si、Al、Fe、Ca、Na、K、Mg）占地壳总重量的99.1％，其它元素只占0.9%，而前15种元素的重量占99.6－99.8％，其余77种元素总重量仅占地壳总重量的0.4－0.2％。
微量元素在地壳中的分布也是不均匀的，它们的丰度可以相差达107倍。
2) 元素克拉克值与周期表对比。
克拉克值大体上随原子序数的增大而减小。
周期表中前26种元素（从H至Fe）的丰度占地壳总重量的99.74％。
但Li、B、Be及惰性气体的含量并不符合上述规律。
周期表中原子序数为偶数的元素总分布量（占86.36％）大于奇数元素的总分布量（占13.64％），相邻元素偶数序数的元素分布量大于奇数元素分布量，这一规律称为奥多－哈根斯法则。
这一规律仍粗略地与太阳系元素的分布规律相同。
这一事实再次说明地球、地壳在物质上同太阳系其它部分的统一性。
图1.12 地球地壳中的组成、主要岩石和主要矿物(Krauskopf et al,1995)3) 若按元素丰度排列，太阳系、地球、地幔和地壳中主要的10种元素的分布顺序是：太阳系：H>He>O>Ne>N>C>Si>Mg>Fe>S地球：Fe>O>Mg>Si>Ni>S>Ca>Al>Co>Na地幔：O>Mg>Si>Fe>Ca>Al>Na>Ti>Cr>Mn地壳：O>Si>Al>Fe>Ca>Na>K>Mg>Ti>H与太阳系和宇宙相比，地球和地壳明显贫H、He、Ne和N等气体元素，表明由宇宙物质形成地球的演化过程必然伴随气态元素的散失。
GeoPlot地球化学数据投图软件-使用说明
GeoPlot地球化学数据投图软件-使用说明关于GeoPlot的使用介绍一．软件介绍GeoPlot软件是用VBA开发并在E某cel中使用的，主要用于地球化学数据计算和投图分析的专业宏程序软件；其实用性强、操作简便，可帮助研究人员快速分析地球化学数据。
功能如下：1)投图功能:(1)某Y散点图、三角图;(2)岩石地球化学判别图;(3)蛛网图及REE配分图2)数据计算:(1)参数的公式计算;(2)CIPW标准矿物计算;(3)AFC及二元混合计算3)辅助功能:(1)样品分组投图;(2)数据统计结果;(3)数据排序及小数位数;(4)分子量计算;(5)Chart2Picture;(6)系统设置及数据库管理4)二．安装删除运行自解压安装包程序文件“GeoPlot.e某e”，选择一个安装目录（一般安装目录为：C:\ProgramFile\GeoPlot）就轻松安装成功，如图1。
启动GeoPlot软件的方法有两种：（1）运行安装目录下的宏程序GeoPlot.某la，自动启动E某cel并出现GeoPlot的文件菜单和工具栏；打开数据文件后，就可以使用该软件的各项功能；（2）在E某cel菜单中加载安装目录下的宏程序GeoPlot.某la，该方法可以让每次E某cel运行后，都会自动启动GeoPlot。
如果需要不自动启动GeoPlot则可在E某cel下取消宏GeoPlot.某la的加载。
该软件是绿色宏程序软件，所以反安装GeoPlot软件则可直接删除所安装的目录，比如“C:\ProgramFile\GeoPlot”。
图1安装GeoPlot的目录选择三．软件使用1．数据准备使用GeoPlot软件进行计算和投图，地球化学数据存放的表格需要复合一个简单的、也是最常用的数据排列格式—元素名称排第一列，样品名称排第一行，其余数据对应元素名称和样品名称排列，如图2。
关于GeoPlot的使用介绍元素名称常用到软件中查询，所以元素名称中不能有非可见字符（网上复制的数据可能产生不可见字符，常导致不能正确计算和投图，需要用户手动重新输入元素名称一遍）。
水文地球化学
当 pH ＜9时，上式中的 (H+)、 (OH-) 和(CO32-)的浓度比 (Ca2+)和 (HCO3-) 相对小得多，因此，可忽略不计，则方程可变为：
2(Ca2+) = (HCO-3) ，或 (Ca2+) = (HCO3-)/2
以上为描述碳酸平衡系统的最基本的方程
5、地下水系统中的碳酸平衡
碳酸平衡
在稀溶液中，挥发性溶质的分压以巴为单位，等于溶 质的摩尔数（亨利定律），适用于难溶气体
CO 2 ( g ) H 2O H 2CO 3 (aq)
K CO 2
H 2CO 3
PCO 2
H 2CO3 H HCO
3
K1
H HCO
HCO H CO
3、温度和压力
温度和压力对某些矿物，如石英、玉髓、非晶质 SiO2的溶解度有显著影响。
在结晶岩地区的热泉口，由于温度和压力降低，SiO2 的溶解度也降低，形成硅华；
碳酸盐矿物的溶解度取决于CO2分压，pH值，温度等，
一般而言， CO2分压越高、温度越低，气体在水中的
溶解度越大，如CO2分压降低则发生脱碳酸作用，形
2、pH值
pH值是决定天然水中许多元素溶解性的重要因素
Na、Ca、NO3和Cl等少数几种离子可存在于各种pH条
件下的天然水中；
多数金属元素在酸性条件下以阳离子形式存在，而当
pH值增高时，即以氢氧化物或偏碱性的盐类形式从水
中沉淀出来(如Fe3+)；
一些金属元素（如Cu、Zn等）发生沉淀的pH略高一 些（Cu为5.3, Zn为7）
地球化学图制作步骤
图 2-1 打开用户文件投影变换（图 2-2 红色箭头所指图标）
图 2-2 打开文件（上述所存文本文件），设置投影参数（以黑龙江省 22 度带为例， 投影数据坐标为大地坐标）： 用户投影参数设置如图 2-3，结果投影参数设置如图 2-4，椭球参数都选择 1： 北京五四/克拉索夫斯基[1940 年]椭球。
12
广西浦北
图 9-7
图 9-8 图 9-9
13
广西浦北
保存裁剪框区文件。将裁剪框区文件从工程文件删除即可进行工程裁剪。 工程裁剪：为各元素重新新建文件夹，以保存裁剪后文件。步骤：其它→工
程裁剪→选择裁剪后文件的存放目录（图 9-10）→选择要裁剪的文件（图 9-11 左侧蓝色部分）→添加→选择全部（图 9-11）→裁剪类型选内裁、裁剪方式选拓 扑裁剪（图 9-12）→参数应用→装入裁剪框→开始裁剪→退出，关闭工程。
>25%~≤75%
浅红（高背景区） 红（高值区）
X 0+0.5S0~ X 0+1.5S0
X X 0+1.5S0~ 0+2.5S0
>75%~≤95% >95%~≤98.5%
深红（强高值区）
＞ X 0+2.5S0
>98.5%
设置完色阶可将存设置，以备以后再用时就不用再重复设置色阶了，单击存
设置...并将原色阶文件覆盖即可保存。之后单击确定，地球化学图就出来了，见 图 7-10。
等量线值 lgµg/g (或 ng/g)
… 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
图 7-8
GeoPlot地球化学数据投图软件-使用说明
GeoPlot地球化学数据投图软件-使用说明一．软件介绍GeoPlot软件是用VBA开发并在Excel中使用的，主要用于地球化学数据计算和投图分析的专业宏程序软件；其实用性强、操作简便，可帮助研究人员快速分析地球化学数据。
功能如下：1)投图功能: (1)XY散点图、三角图;(2)岩石地球化学判别图;(3)蛛网图及REE配分图2)数据计算: (1)参数的公式计算;(2)CIPW标准矿物计算;(3)AFC及二元混合计算3)辅助功能: (1)样品分组投图;(2)数据统计结果;(3)数据排序及小数位数;(4)分子量计算;(5)Chart 2Picture;(6)系统设置及数据库管理4)二．安装删除运行自解压安装包程序文件“GeoPlot.exe”，选择一个安装目录（一般安装目录为：C:\Program Files\GeoPlot）就轻松安装成功，如图1。
启动GeoPlot软件的方法有两种：（1）运行安装目录下的宏程序GeoPlot.xla，自动启动Excel并出现GeoPlot的文件菜单和工具栏；打开数据文件后，就可以使用该软件的各项功能；（2）在Excel菜单中加载安装目录下的宏程序GeoPlot.xla，该方法可以让每次Excel运行后，都会自动启动GeoPlot。
如果需要不自动启动GeoPlot则可在Excel下取消宏GeoPlot.xla的加载。
该软件是绿色宏程序软件，所以反安装GeoPlot软件则可直接删除所安装的目录，比如“C:\Program Files\GeoPlot”。
图1 安装GeoPlot的目录选择三．软件使用1．数据准备使用GeoPlot软件进行计算和投图，地球化学数据存放的表格需要复合一个简单的、也是最常用的数据排列格式—元素名称排第一列，样品名称排第一行，其余数据对应元素名称和样品名称排列，如图2。
元素名称常用到软件中查询，所以元素名称中不能有非可见字符（网上复制的数据可能产生不可见字符，常导致不能正确计算和投图，需要用户手动重新输入元素名称一遍）。
Tin法制作地球化学图
首先把EXCEL表Y坐标放第一列，另存为CSV格式。
1：50000地化图绘制：1、用记事本打开Ag.csv格式数据，将x（坐标符号），y（坐标符号），Ag （元素符号）改为notgrid，另存为Ag.det格式。
对于1：5万化探，x、y单位为m，成图比例尺为1：50000万，那么将x、y同除以50，就变成mm。
2、依次打开MAPGIS程序主菜单\空间分析\DTM分析（MAPGIS数字地面模型自系统），打开Ag.det文件。
3、打开Tin模型\快速生成三角剖分网。
4、打开Tin模型\删除三角剖分网边。
5、打开Tin模型\三角网内插网格化：左下角坐标设为最小整公里坐标，右上角坐标设为最大整公里坐标；1：50000地化图网格化间距实验证明设为X：125，Y：125最为适当，另存栅格化结果文件Ag.Grd.6、打开Grd模型\平面等值线图绘制：选择：等值线套区、绘制色阶、保留边界线、等值线光滑处理（高程度）、示坡线、线型。
线参数设定：默认。
注记参数设定：标注字体尺寸：2mm，小数部分位数：与该元素线值位数相同，其它默认。
(.VEL)。
以后重新绘图时，可直接装入该文件。
7、确定成图，保存点、线、面文件。
8、极值点标注：等值线图极值点标注，只能在“三角剖分-追踪剖分等值线后”，在“模型应用-高程点标注制图”子菜单下标注。
在标注输出设置中，选择“等值区内部极值点标注”，然后按自己的要求设网上收集资料一、建立数据表1、用电子表格建立成图元素表（不要表头，建议将平面坐标按Y、X、数据…列次序移到表前1、2…列，Y、X坐标分别为7位、6位整数）；2、另存为文本文件（制表符分隔）类型的新文件并取名(如文件名为茨竹箐原始数据表)。
二、Mapgiss投影转换1、打开Mapgiss程序主菜单→点击″实用服务″→点击″投影变换″→点击″P投影转换″→点击″U用户文件投影转换″→打开文本文件(茨竹箐原始数据表)；2、选″指定分隔符″→指定坐标X、Y列（注意X、Y列可能倒置）；3、点击″设置分隔符″→选″Tab″及″空格″→点击″数据类型″→将所有的字符串点为″双精度″→在″属性名称所在行″选框中选择数据（注意数据不能选″无″）→在″加入″列选定所需的坐标及元素值→确定；4、点击″点图元参数″→选择子图参数（建议选编号35的圆点的子图号，高度宽度一般均为1，颜色为1）→确定；5、选″不需要投影″→点击″数据生成″→确定→关闭；6、另取点文件名（如茨竹箐数据投影)→保存→退出。
实验三 地球化学制图
实验四地球化学制图
一实验目的与要求
通过实验，加深对地球化学场和地球化学异常的理解，掌握元素地球化学图和地球化学剖面图的绘制方法。
二实验内容与实验步骤：
（一）元素地球化学图(单元素含量等值线图)的绘制
作图步骤如下：
1 根据工作比例尺，并在相应的比例尺下的图纸上，用实心小圆点标注出测区内所有的采样点位置及其相应的位置编号（图1）；
2 在每个采样点位置标出该点样品的元素含量分析结果，一张图上一般只标注一种元素，最多不超过两种元素；
3 按照对数或算数等间隔采用内插法勾绘等值线。
通常，如果异常强度很高（例如异常强度达到背景值的数倍、十数倍乃至上百倍）时，采用对数等间距；异常强度不高时，则采用算数等间距。
本实验以4ppm 为等值线间隔进行作图。
（二）地球化学剖面图的绘制
作图步骤如下：
1 元素含量按照一定的比例尺，绘出所在测线长度，并按照一定间距标出相应位置标号；
2 将剖面上采集样品的元素分析结果按照一定的纵向比例尺于剖面相应采集点位处的上方（与侧线水平轴线垂直方向）通过实心小圆点标出；
3 将相邻测点上对应的标志点用线段连接在一起，则会形成一条起伏变化的曲线，成为元素含量剖面图；
表1为内蒙某地Ⅴ号剖面采样点位及W元素的含量，请按照1：10000的横向比例尺，1：20（1cm＝20ppm）的纵向比例尺绘制一张元素地球化学剖面图。
表1 内蒙某地Ⅴ号剖面采样点位及W元素的含量
实验报告要求内容：
一、实验目的与要求；
二、实验方法与步骤；
三、绘制地球化学剖面图；
四、绘制元素地球化学图；
图1 某区Co元素地球化学元素含量数据图。
地球化学图相关计算
地球化学图相关计算 Document number：WTWYT-WYWY-BTGTT-YTTYU-2018GT地球化学图相关计算一、计算Al2O3/(CaO+Na2O+K2O)和Al2O3/(Na2O+K2O)或A/CNK-A/NK分子(摩尔)比步骤与投图1．列出Al2O3、CaO、Na2O和K2O的分子量Al2O3分子量CaO分子量Na2O分子量K2O分子量主要氧化物分子量氧化物Al2O3CaO Na2O K2O Si2O TiO2Fe2O3FeO MnO MgO P2O5分子量2．列出所测样品的Al2O3、CaO、Na2O和K2O含量K2O含量岩性样品号Al2O3含量CaO含量Na2O含量(xηγT1) P4XT06(xηγT1) P15XT02(zηγT1) P15XT03-1(zηγT1) P4XT103．求出分子(摩尔)数(样品中Al2O3、CaO、Na2O和K2O含量除以各自氧化物的分子量)岩性样品号Al2O3分子数CaO分子数Na2O分子数K2O分子数(xηγT1) P4XT06 （=）（=）（=）（=）(xηγT1) P15XT02 （=）（=）（=）（=）(zηγT1) P15XT03-1 （=）（=）（=）（=）(zηγT1) P4XT10 （=）（=）（=）（=）4．分别求出Al2O3/(CaO+Na2O+K2O)和Al2O3/(Na2O+K2O)或A/CNK-A/NK分子 (摩尔) 比Al2O3/(Na2O+K2O)分子比（摩尔）岩性样品号Al2O3/(CaO+Na2O+K2O)分子比（摩尔）(xηγT1) P4XT06 ++= +=(xηγT1) P15XT02 ++= +==(zηγT1) P15XT03-1 ++= +=(zηγT1) P4XT10 ++= +=5．投图以Al2O3/(CaO+Na2O+K2O) 或A/CNK分子（摩尔）比X为轴和以Al2O3/(Na2O+K2O)或A/NK 为Y轴，投点即可。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
地球化学图相关计算
一、计算Al2O3/(CaO+Na2O+K2O)和Al2O3/(Na2O+K2O)或A/CNK-A/NK分
子(摩尔)比步骤与投图
1．列出Al2O3、CaO、Na2O和K2O的分子量
Al2O3分子量CaO分子量Na2O分子量K2O分子量
102.00 56.00 62.00 94.00
主要氧化物分子量
氧化物Al2O3CaO Na2O K2O Si2O TiO2Fe2O3FeO MnO MgO P2O5
分子量102.0 56.0 62.0 94.0 60.0 80.0 160.0 72.0 71.0 40.0 142.0
2．列出所测样品的Al2O3、CaO、Na2O和K2O含量
岩性样品号Al2O3含量CaO含量Na2O含量K2O含量
(xηγT1
)P4XT06
14.44 0.41 2.66 5.52
(xηγT1
)
P15XT0213.31 0.41 4.33 4.47
(zηγT1
)
P15XT03-113.34 0.57 4.32 4.25
(zηγT1
)
P4XT1013.82 0.56 4.23 4.30
3．求出分子(摩尔)数(样品中Al2O3、CaO、Na2O和K2O含量除以各自氧化物的分子量)岩性样品号Al2O3分子数CaO分子数Na2O分子数K2O分子数
(xηγT
1)P4XT06（14.44/102.0=）
0.142
（0.41/56.0=）
0.007
（2.66/62.0=）
0.043
（5.52/94.0=）
0.059
(xηγT
1)P15XT02（13.31/102.0=）
0.130
（0.41/56.0=）
0.007
（4.33/62.0=）
0.070
（4.47/94.0=）
0.048
(zηγT
1)P15XT03-1（13.34/102.0=）
0.131
（0.57/56.0=）
0.010
（4.32/62.0=）
0.070
（4.25/94.0=）
0.045
(zηγT
1)P4XT10（13.82/102.0=）
0.135
（0.56/56.0=）
0.010
（4.23/62.0=）
0.068
（4.30/94.0=）
0.046
4．分别求出Al2O3/(CaO+Na2O+K2O)和Al2O3/(Na2O+K2O)或A/CNK-A/NK分子 (摩尔) 比岩性样品号Al2O3/(CaO+Na2O+K2O)分子比（摩尔）Al2O3/(Na2O+K2O)分子比（摩尔）
(xηγT1)P4XT060.142/(0.007+0.043+0.059)=1.299 0.142/(0.043+0.059)=1.393 (xηγT1)P15XT020.130/(0.007+0.070+0.048)=1.046 0.130/(0.070+0.048)==1.112 (zηγT1)P15XT03-10.131/(0.010+0.070+0.045)=1.046 0.131/(0.070+0.045)=1.138
(zηγT1)P4XT100.135/(0.010+0.068+0.046)=1.093 0.135/(0.068+0.046)=1.189
5．投图
以Al2O3/(CaO+Na2O+K2O) 或A/CNK分子（摩尔）比X为轴和以Al2O3/(Na2O+K2O)或A/NK 为Y轴，投点即可。
换算公式 n = m / M，n是摩尔数，m是质量或含量，M是分子量。
摩尔比和分子比可以看做是一个意思。
二、微量元素蛛网图K、P、Ti等的计算
1. K=K2O×0.83013×10000
2. P=P2O5×0.43646×10000
3. Ti=TiO2×0.5995×10000
三、花岗岩R2-R1构造环境判别图解（阳离子或原子数之比）R2和R1的计算步骤
1．列出R
2和R
1
R1=4Si-11(Na+K)-2(Fe+Ti)，R为阳离子数或原子数。
R2=6Ca+2Mg+Al，R为阳离子数或原子数。
2．列出氧化物中阳离子计算步骤
（1）列出样品中氧化物含量%（表1）
（2）列出氧化物分子量（表1）
（3）列出氧化物阳离子数（表1）
表1 平均英云闪长岩的阳离子计算
（4）对应各氧化物分别求出Si、Na、K、Fe、Ti、Ca、Mg、Al阳离子比例（表1）
计算公式：
Si=（样品中SiO2含量/SiO2分子重量）×阳离子个数=（61.52/60.09）×1.00=1.0238 Na=（样品中Na2O含量/Na2O分子重量）×阳离子个数=（3.63/61.98）×2.00=0.1171 K=（样品中K2O含量/K2O分子重量）×阳离子个数=（2.07/94.20）×2.00=0.0439
Fe=（样品中FeO含量/FeO分子重量）×阳离子个数=（3.82/71.85）×1.00=0.0532或（FeO+ Fe2O3）阳离子之和。
Ti=（样品中TiO2含量/ TiO2分子重量）×阳离子个数=（0.73/79.90）×1.00=0.0091 Ca =（样品中CaO含量/ CaO分子重量）×阳离子个数=（5.42/56.08）×1.00=0.0966 Mg =（样品中MgO含量/ MgO分子重量）×阳离子个数=（0.08/40.30）×1.00=0.0695 Al=（样品中Al2O3含量/ Al2O3分子重量）×阳离子个数=（16.48/101.96）×2.00=0.3233（5）求出Si、Na、K、Fe、Ti、Ca、Mg、Al千阳离子数（表1）
Si、Na、K、Fe、Ti、Ca、Mg、Al阳离子比例×1000，得出Si、Na、K、Fe、Ti、Ca、Mg、Al千阳离子如表1。
3．求出R2和R1阳离子数或原子数
将Si、Na、K、Fe、Ti各千阳离子数代入公式R1=4Si-11(Na+K)-2(Fe+Ti)，求出R1阳离子数或原子数；同理，将Ca、Mg、Al千阳离子数代入公式R2=6Ca+2Mg+Al求出R2阳离子数或原子数。
四、其它相关计算
1. Al=Al2O3×5
2.9%
2. FeOt= FeO+0.8998 Fe2O3
3. 氧化指数 FeOt/MgO
4. Mg#（镁指数）
计算方式分为如下两步：
（1）FeOT (wt.%) = FeO (wt.%) + Fe2O3 (wt.%) *0.8998 = FeO (wt.%) + Fe2O3 (wt.%) * (71.844/ (159.6882/2))
（2）Mg#=100*(MgO/40.3044)/(MgO/40.3044+FeOT/71.844)
上述是分别测试分析了FeO和Fe2O3的计算方法，如果是测试的全铁，也可以近似计算。
通常说的高Mg，是指岩石具有较高的MgO含量，如火山岩中的高镁安山岩（通常情况下，异常高的MgO含量指示着可能有地幔物质参与，如俯冲带地幔楔或者软流圈熔体上涌等等）。
Mg#（镁指数）也可以定量的表示岩石中的Mg含量高低。
Mg#通常用于镁铁质岩石，可以粗略指示地幔岩石的部分熔融程度，高Mg#的地幔橄榄岩可能经历了更高程度的部分熔融，常在92-93左右，而原始地幔会相对富集，Mg#较低，在88-89左右。
（注：文档可能无法思考全面，请浏览后下载，供参考。
）。