地球化学图相关计算
14-微量元素地球化学
5.微量元素的岩石成因
硅质岩成因
5.微量元素的岩石成因
太原组硅质岩的热液成因解释
5.微量元素的岩石成因
变质岩原岩恢复
地壳不同变 质原岩的 REE与 La/Yb比值判 别图，可用 于区分不同 类型的玄武 岩、花岗岩 和碳酸盐岩
5.微量元素的岩石成因
重要元素对岩石成因的指示意义
元素 特征解释
高度相容元素，N i 和C o 赋存在Ol 中，C r在Sp和C px中，这些元素的高度富集（例如N i ＝250-300 ppm, C r=500Ni, Co, Cr 600 ppm）暗示着岩石母岩为地幔橄榄岩，如果岩石系列显示N i 的逐渐降低 （C o 也可以显示同样规律） 则预示着 Ol 的分离结晶作用。C r的逐渐降低代表尖晶石或者C px的分离结晶作用。 它们在部分熔融和分离结晶过程中显示相似的特征。都倾向于进入Fe-T i 氧化物（钛铁矿和钛磁铁矿)， 是钛铁矿和 钛磁铁矿结晶分异的示踪剂。如果V和T i 显示差异性质， 则T i 可以类质同象进入一些副矿 物相，例如榍石和金红石。 极不相容元素，基本不进入主要的地幔矿物相，有时可以与T i 类质同象进入副矿物相，例如榍石和金红石。 不相容元素，在钾长石。云母或者角闪石中可以替换K。Rb在角闪石中类质同象替换能力弱于在钾长石和云母中，因 此K/Ba比值可以用来鉴别这些矿物相。 在Pl 中容易类质同象替换C a（但是在Py中不取代C a），在钾长石中替换K，在浅部低压条件下当Pl 作为早期结晶相 的时候，显示相容元素特征，因此Sr或者C a/Sr比值是鉴别Pl 参与的有力指示剂。但是Sr在高压的地幔条件下，当Pl 不 稳定的时候，显示不相容元素特征。 石榴石(Opx和角闪石稍弱)容纳重稀土元素，因此会导致轻稀土的分异。榍石和Pl 倾向于吸纳轻稀土元素。C px仅导 致REE轻度分异。Eu强烈倾向进入Pl 中，因此Eu异常可以鉴别是否有Pl 的参与。 常类似于HREE，显示不相容元素特征。强烈倾向进入石榴石和角闪石中，辉石次之。榍石和磷灰石也富集Y，因此如 果岩石中存在这些副矿物，将明显影响Y的分异。
岩石地球化学计算
岩石地球化学计算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异常。
地球化学-稀土元素标准化计算
表中数据为山东济南辉长岩、沂南花岗岩7件样品的REE组成(PPm)1, 用球粒陨石值对样品的REE组成进行标准化，作其分配模式图，对图件中表达的地球化学特征进行说明；2, 计算各样品的Eu/Eu* ,并对其地球化学意义进行说明；,3,假设辉长岩中造岩矿物的组成为：CPX45%, PL35%, OL20%。
结合课件中提供的REE在矿物和熔体间的分配系数，计算与辉长岩平衡的熔体的REE组成, 并作REE配分模式图。
解答:1如下表1-1为常用球粒陨石和原始地幔稀土元素组成，我采用C1球粒陨石数据(Sun & MCDOnough,1989)对样品的REE进行标准化，得到了下表1-2,再根据对样品REE标准化的数据进行作样品的分配模式图，得到了图1-1表1-1MJN0608MJNO607MJN0609MJN0606MYN0625MYN0625MYNO607La21. 055 25. 570 27. 38476 371489. 451 337.975464,135 Ce20, 261 24. 837 124.673 66 340 369* 281 254, 902 341. 503 Pr18. 421 22. 421 ΞL E7952； 421 235. 789 164, 211 205. 263 Nd17. 880 32. 270 21. 370 46. 467 165. 739 114. 347 131.692Sm14. 96717. 320 16. 60132.026 75. 163 50. 327 46. 993EU13. 793 14. 4S3 19. 130 27÷41424.138 16.897 21, 379Gd9. 732 11.290 IL 33S Ξ0.00035,961 25, 937 Γ 18,735Tb8. 824 9. 626 9. 626 17. 112 Ξ8, 34219, 786 11. 230 Dy7. 953 8. 7019. 094 16. 024 24. 291 16. 811 7. 795HO7. 067 7. 774 8. 12714. 311ΞC. S4S14. 841 6L 007Er 5. 921 6. 6477. 130 12, 085 181852 12. 50S 5. 498Tm 5. 098 5. 4&0 1 6. 275 10, &80 18, 039 IL 765 Γ 0. 490Yb 5. 588 6. 059 6. 706 12. 294 18, 647 12. 471 6. 176LU 5. 118 5. 118 5. 906 11. 024 16. 929IL 811 5. 906表1-2→-IJN06Q8 →-IJN0607IJN0609IJNO 606 T^iYNo &药 →-lYN0625 -^lYNO 607图1-1通过对样品配分模式图进行分析可知道，沂南花岗岩样品中富集轻稀土元素 而亏损重稀土元素，这与花岗岩的成分岩性有一定关系， 花岗岩为酸性岩，主要 矿物为长石、石英和云母，而这矿物主要富集轻稀土元素，并且从图中可以看出 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等气体元素，表明由宇宙物质形成地球的演化过程必然伴随气态元素的散失。
地球化学图制作步骤
打开数字地质调查综合信息平台。打开工程文件原始数据点，综合数据处理 →等值线图→出现对话框（图 11-1），依次选择点图层→字段名称（各元素名， 如 Au）→单击计算数字特征→单击下一步。出现对话框（图 11-2 设置），选择 Tin 方法→自定义 文件名改为 Au 极值→单击下一步→出现设置等值线参数对 话框，将步骤 7 中保存的色阶文件导入（单击装设置...找到色阶文件打开）→确 定→上一步→上一步，回到了图 11-1 对话框→字段名称改为其它元素（如 Ag） →单击计算数字特征→单击下一步。出现对话框（图 11-2 设置）→自定义 文件 名改为 Ag 极值→单击下一步→出现设置等值线参数对话框，将步骤 7 中保存的 色阶文件导入→单击确定→上一步→如此往返直到所有元素都做完为止。
10
广西浦北
往返，直到添加完所有元素为止，关闭工程文件。 9、裁剪地球化学图：
打开元素地球化学图的工程文件（如 Ag），观察图件边框外的区，如出现图 9-1 情形（红色箭头所指，应该叫窜色），就得加以修改，特别是图 9-1 所示在边
图 9-1
图 9-2
图 9-3（左侧为修改前、右侧为修改后） 框内就已经窜色的，以备裁剪后易于修改。修改后效果见 9-2。出现图 9-3 左侧
3.2
1.6
4.0
1.7
5.0
1.8
6.3
1.9
7.9
2.0
10
…
图上标注的真值 µg/g (或 ng/g)
12.5 16.0 20.0 25.0 32.0 40.0 50.0 63.0 79.0 100.0 …
图 7-9
9
广西浦北
新规范地球化学图色区设置见表 7-2，建议采用表 7-2 色区设置。
地球化学图制作步骤
图 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
地球化学+地质填图教程总结
地球化学+地质填图常用教程总结一、1:5万土壤采样点位设计(水产沟幅为例)1、标准图框生成首先利用Mapgis—实用服务—投影变换—系列标准图框—生成1:5万图框.图框模式：高斯坐标实线公里网.投影参数：起点纬度491000. 起点经度1210000. 网间间距1KM.图框文件名：水产沟图框.图框参数输入：图框内容：图幅名称：水产沟幅.图框参数选择：输入并绘制接图表+绘制图幅比例尺+图框外图廓线.2、图形校正图片JPG转换成msi文件： Mapgis—图像处理—图像分析—文件—数据输入—添加文件(图片JPG)—转换—转化成msi文件.校正到标准图框内：Mapgis—图像处理—图像分析—文件—打开影像—msi文件—镶嵌融合--打开参照文件—参照线文件(标准框). 镶嵌融合—添加控制点—(先点击带校正图片某点，然后点击对应的标准框点，最后按两次空格键，选择是，以此类推，尽量多的找几处控制点，通常不少于13个控制点，要求残差<1为准)—镶嵌融合—校正预览—透明度调节--回车--镶嵌融合—影像几何校正—几何校正--重命名保存文件.3、图片数字化对地形线、河流、公路、村庄等地形要素进行数字化过程. 同时，准备大格号、图例、网格线、采样位置、重复样和监控样分析表格等文件.4、采样点位设计四幅图框：据前述方法，生成5万四幅合体的图框（注意修饰）。
将四幅图之间的界线用红线标出，根据阵列复制线，画出网格线，在网格线的交汇处定出大格号标注，然后建立采样点位文件，在上面点击采样点位置，辅助工具1—导入导出功能—点位置转属性. 辅助工具2，水系沉积物自动编号，0001大格号的左上角坐标为5488000,355000，一行上的大格号从0001—0004，点击确定，重新显示文件，出来小格号标注文件.5、重复样、监控样在每个单幅的采样点位图中，设置重复样和监控样，按从上到下，从左到右的原则，即a—c—b —d，查出每50个样品，其中1个是重复样，4个是监控样，45个是普通样品，导出excel表格之后排序，对其重复样和监控样设置，然后在单幅采样点位图中标出来.二、1:1万土壤采样点位设计1、基点基线准备在地形图基础上，可以在图框文件里输入一个点，作为基点（待设计完毕可以删除，例如基点为200/100，点号为200，线号为100），图框外部边缘作为基线。
实习一地球化学剖面图绘制
• 采用目估法确定异常下限； • 按等差法对异常进行分级；
• 分析数据，采用内插法进行 等值线勾绘。
三、方法
• 一、异常下限确定方法 1、经验法 2、长剖面法 3、计算法
三、方法
• 一、异常分级方法 1、0.1lg分级 2、等差、等分级 3、累积频率分级
三、方法
• 三、等值线内插法
（板书演示）
四、坐标模版
五、思考题
• 目估法确定背景值和异常下限需要哪些 前提？ 植物中元素指示意义受哪些因素影响？ 完整的地球化学剖面图需含哪些方面的 内容？
• •
实习二 地球化学平面Biblioteka 绘制 （4学时）一、实习目的
1、学会手工绘制地球化学 平面图方法； 2、学会确定异常下限和异 常分级。
二、实习要求
三、原始数据
福建峰岩铅锌矿区植物地球化学测量（蕨菜）结果
距起点(m) Pb（μg/g） Zn（μg/g） As（μg/g） Hg（ng/g） 距起离(m) Pb（μg/g） Zn（μg/g） As（μg/g） Hg（ng/g） 0 17.3 21.8 0.18 21.2 2050 22 96.9 0.45 17. 2 100 83.5 25.6 0.48 25.2 2150 74.2 132.5 0.73 22.2 200 28.5 25.2 0.33 19.17 2450 31.4 138.6 0.43 13.1 300 44.6 132.4 0.78 16.1 2550 155.7 458.3 0.66 23.2 400 72.4 76.8 0.453 18.2 2650 114 504.2 0.60 16.1 500 33.4 36.5 0.58 19.2 2750 50.2 206.8 0.55 19.2 700 10 32.9 0.35 18.2 2850 28.4 210.8 0.73 18.2 800 36.5 31.3 0.38 17. 2 2950 23.9 94.6 0.48 15.1 900 13.9 48.3 0.40 20.2 3050 27 83.1 0.55 16.1 1000 9 32.5 0.48 31.3 3150 33 99.4 0.55 24.2 1050 10.6 72.2 0.45 16.1 3250 29.4 67.8 0.68 15.2 1250 12 41.7 0.40 21.2 3350 24.4 83.1 0.45 25.2 1750 10.1 52.4 0.40 13.1 3450 15.3 46.1 0.50 19.2 1850 10.4 64.2 0.33 16.1 3550 22.3 48.1 0.55 15.1 1950 10.7 57.6 0.76 12.1
地球化学图相关计算
地球化学图相关计算一、计算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.02．列出所测样品的Al2O3、CaO、Na2O和K2O含量岩性样品号Al2O3含量CaO含量Na2O含量K2O含量(xηγT1)P4XT0614.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.303．求出分子(摩尔)数(样品中Al2O3、CaO、Na2O和K2O含量除以各自氧化物的分子量)岩性样品号Al2O3分子数CaO分子数Na2O分子数K2O分子数(xηγT1)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ηγT1)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ηγT1)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ηγT1)P4XT10（13.82/102.0=）0.135（0.56/56.0=）0.010（4.23/62.0=）0.068（4.30/94.0=）0.0464．分别求出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.1895．投图以Al2O3/(CaO+Na2O+K2O) 或A/CNK分子（摩尔）比X为轴和以Al2O3/(Na2O+K2O)或A/NK 为Y轴，投点即可。
地球化学图相关计算
地球化学图相关计算 The manuscript was revised on the evening of 2021地球化学图相关计算一、计算Al2O3/(CaO+Na2O+K2O)和Al2O3/(Na2O+K2O)或A/CNK-A/NK分子(摩尔)比步骤与投图1．列出Al2O3、CaO、Na2O和K2O的分子量Al2O3分子量CaO分子量Na2O分子量K2O分子量主要氧化物分子量氧化物Al2O3 CaO 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 根据工作比例尺，并在相应的比例尺下的图纸上，用实心小圆点标注出测区内所有的采样点位置及其相应的位置编号（图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轴，投点即可。
地球化学-稀土元素标准化计算
表中数据为辉长岩、沂南花岗岩7件样品的REE组成(ppm)1，用球粒陨石值对样品的REE组成进行标准化，作其分配模式图，对图件中表达的地球化学特征进行说明；2，计算各样品的Eu/Eu*，并对其地球化学意义进行说明；,3，假设辉长岩中造岩矿物的组成为：CPX45%，PL35%，OL20%。
结合课件中提供的REE在矿物和熔体间的分配系数，计算与辉长岩平衡的熔体的REE组成，并作REE配分模式图。
解答：1，如下表1-1为常用球粒陨石和原始地幔稀土元素组成，我采用C1 球粒陨石数据（Sun & McDonough,1989）对样品的REE进行标准化，得到了下表1-2，再根据对样品REE标准化的数据进行作样品的分配模式图，得到了图1-1表1-1表1-2图1-1通过对样品配分模式图进行分析可知道，沂南花岗岩样品中富集轻稀土元素而亏损重稀土元素，这与花岗岩的成分岩性有一定关系，花岗岩为酸性岩，主要矿物为长石、石英和云母，而这矿物主要富集轻稀土元素，并且从图中可以看出Eu的负异常，说明在岩浆结晶形成花岗岩之前就有长石结晶出来，使岩浆呈Eu 的负异常。
辉长岩的样品配分模式图表现出来的富集轻稀土元素没有沂南花岗岩样品那么显著，富集程度较低，这也与辉长岩的岩性成分有关，辉长岩中主要矿物为辉石和长石，长石富集轻稀土元素较为显著，而辉石相对较富集重稀土元素，但程度不是很显著，所以岩石总体表现较为富集轻稀土元素，但程度不是那么显著。
并且从图中可以看出Eu的正异常，只是不是很显著，说明长石结晶出来使岩石呈Eu的正异常。
2，Eu/Eu*=2×Eu/(Sm+Gd)(其中Eu、Sm、Gd都是为球粒陨石标准化值)，根据这个求出各样品中的Eu/Eu*，如下表1-3：表1-3由上表中的Eu/Eu*值可知的辉长岩为Eu的正异常，说明在岩浆结晶时，长石和辉石先结晶出去形成辉长岩，而长石中富集Eu元素，所以在辉长岩中Eu 为正异常，而后期岩浆因长石的结晶分异而呈Eu的负异常，并且逐渐向酸性过渡，结晶形成酸性岩。
地球化学图和原始图件的编制规范(正式)
a.采样点要分布均匀、合理。要使1：50000地形图上按1km2划出的采样格子（称小格）中的绝大多数（95%以上）小格内都有采样点分布，应做到基本不出现或很少出现连续3个以上空白小格。采样密度平均为1~2个点/km2时，小格内样品数一般不要超过2个。采样点分布均匀并不意味着将采样点布置在采样格子中央，而应将采样点布置在每一个格子中能最大限度控制汇水面积处（见图2）；
GB/T14839地球化学勘查技术符号
3总则
3.1区域化探是基础性、公益性、战略性基础地质调查工作之一。工作范围应覆盖我国全部可工作的国土面积。其调查成果，主要为区域矿产勘查中成矿远景区和中大型找矿靶区的圈定、资源潜力预测、地质填图与基础地质研究、环境与土地质量评价以及社会其它行业和领域的应用提供数十种元素含量空间分布的地球化学资料。此项工作按国际分幅部署测区。基本工作比例尺为1：250000，交通很不发达、工作条件困难的边远地区也可采用1：500000工作比例尺。它们的基本采样密度（以水系沉积物测量为例）分别为1点/1~4km2和1点/16km2。成果图件视不同行业、不同领域和不同群体对地球化学资料的需求，可制作任意比例尺和不同图示表达形式的图件（册）提供应用。
d.重要异常的查证：包括拟采用的查证方法、使用的现场分析手段和设备、查证程度、提交异常查证的成果图件、简报等。
e.技术、经济指标：包括设备、材料计划、经济预算、组织编制和工作进度安排等；
f.提交成果的内容与时间：包括计划上交的报告、图件等资料，提交时间等。
4.5区域化探项目都应单独编写设计书。设计书由项目负责人负责编写。
3.2按地理特点和区域化探工作条件，可把我国大致分为以下两类地区（见图1）：
a.内地及沿海中低山、丘陵地区；
13-微量元素地球化学
1000
碳酸盐岩
La/Yb
100
沉积岩 钙质泥岩
N
EuN，SmN和GdN均为相应元素实测值的球粒陨石标
准化值。δEu(或Eu/Eu*)>1为正异常，δEu<1为负
异常，δEu＝1无异常。
4.稀土元素地球化学
A negative Eu anomaly is typical of many
continental rocks, as well as most sediments and seawater.
上次课回顾
3.岩浆过程中微量元素定量模型
两花岗岩体，经采样，测得La、Sm含量（ppm）分别为：
样 品
花岗岩A
花岗岩B
La 7.5 11.0 33.1 38.2 20.5 42.3 50.2 30.8 38.4 68.5
Sm 6.0 4.5 5.8 6.2 6.1 8.4 10.1 7.9 5.97 13.9
REE球粒陨石标准化图 解，表示Eu异常的计算
4.稀土元素地球化学
δEu(或Eu/Eu*)计算以曾田彰正-科里尔图解为根据，
无Eu异常时，Eu的应有含量值为标准化曲线上旁侧
两个元素Sm和Gd的丰度值以内差法求得。δEu(或
Eu/Eu*)按下式得出:
δEu ＝ Eu/Eu*＝
EuN
(
Sm
2
Gd
)
4.稀土元素地球化学
中国西南地区76种元素地球化学填图
第29卷 增刊物探化探计算技术2007年10月收稿日期:2007-06-30文章编号:1001)1749(2007)增刊(1))0174)06中国西南地区76种元素地球化学填图程志中,谢学锦(中国地质科学院 地球物理地球化学勘查研究所,河北廊坊065000)摘 要:人类生活的地球是由元素周期表中的所有元素构成的,了解所有元素在地球的分布规律对于人类认识自然以及资源与环境研究,具有十分重要的作用。
利用西南地区区域化探扫面的副样,采用1B 50000图幅一个组合成一个分析样的方式,收集水系沉积物样品,研制开发了76种元素分析方案,特别是对一些困难元素进行重点研究,使得所有元素的分析检出限低于地壳丰度值。
编制完成76种元素的地球化学图,利用地球化学资料进行资源远景的初步评价,区域成矿的划分等基础性研究工作。
关键词:76种元素;地球化学填图;分析方法;成矿区带中图分类号:P 632 文献标识码:A0 前言人类生存所依赖的地球环境,是由元素周期表中所有元素组成,了解周期表中所有元素的空间分布情况,分析它们的空间分布规律将提供前所未有的新信息:找寻大型矿床;探索新型资源利用的可能性;为区域成矿学提供新的研究方向。
1962年H.E .H a w kes 和J .S .W ebb 在5矿产勘查的地球化学6[1]一书中首次指出地球化学填图的重要意义。
指出:/将元素的变化来制图也是研究地球的一种办法,与制作土壤图,岩石类型图,构造图和变质作用图件的作用一样。
01973年J .S .W ebb 等出版了第一部地球化学图集5北爱尔兰地球化学实验图集6[2]。
从1975年W ebb [3]及他在伦敦帝国学院的同事们开始在赞比亚、塞拉利昂、马来西亚、爱尔兰等地用稀密度的水系沉积物采样及多元素分析进行区域性地球化学填图。
从1973年至今,世界上共有50余项区域性及国家性地球化学填图计划已完成、或正在进行[4~16],分析元素的数目从数种至最多40种,但由于测定的元素种类较少,且很不一致,痕量或次痕量元素信息量不够,特别是对一些痕量及超痕量元素分析检出限偏高或根本未分析,取得结果不能全球对比,使得绝大部分地球化学图远远达不到地球化学图对资源与环境问题本应做出的贡献。
地球化学图相似性定量评价方法
ｓｍｉ ｒ ｙ ｏ ｈ ｐ ｏ ｌ ｅ ｅ ａ ｕ ｔ ｄ ｑ ａ ｔ ａｉ ｅｙ Ｔ ｉ ｍｅｈ ｄ ｃ ｎ ｂ ｓ ｐ ｌ ｄ ｔ ｈ ｅ ｅ ｒ ｈ ｏ ｉ ｌ ｉ ｆ ｔ ｅ ｍａ ｓ ｃ ｕ ｄ ｂ ｖ ｌａ ｅ ｕ ｎ ｉ ｔ ｌ ． ｈ ｓ ａｔ ｔ ｖ ｔ ｏ ａ ｅ ａ ｏ ａ ｐ ｉ ｏ ｔ ｅ ｒｓ ａ ｃ ｎ ｌ ｅ
Ｑｕ ｎｉ ｔ ｅＡｓｅｓ ｎ ｅｈ ｄ ｆ ｒｔｅＳｍｉ ｒｔ ｆＧｅ ｃ ｅ ｃ ｌ ａ ｓ ａ ｔ ａｉ ｓ ｓｍｅ ｔ ｔ ｏ ｏ ｈ ｉ ｌ ｉ ｏ ｏ ｈ ｍｉａ ｐ ｔ ｖ Ｍ ａ ｙ Ｍ
ＪＡ ｎ Ｉ ＮＧ Ｈｏ ｇ
（ ｉ ｎＺｊ ｎｎ ｎ ｔｌ ｒｙＴ ｃ ｎ ｌｇ ．Ｌｄ ，Ｘａ ｎ ３ １ ０ ， ｉａ Ｘａ ｍｅ ｉ Ｍｉｉｇａ ｄ Ｍｅａｌ ｇ ｅ ｈ ｏ ｙＣｏ ｔ． ｉ ｉ ｎ ｕ ｏ ｍｅ ６ １ Ｃｈｎ ） １
ｅｅ ｎ ｉｔ ｂ ｔｏ ｎ ｅ ｏ ａ ｎ ｐ ｔ ｌｖ ｒａｉｎ ｏ ｎ ｉｏ ｍ ｅ ｔ ｌｐｏｌ ｉｎ ｒ ｇｏ ｓ ｌ ｍｅ ｔｄ ｓｒ ｕ ｉｎ ａ ｄ ｔｍｐ ｒ ｌａ ｄ ｓ ａｉ ａ ｉｔｏ ｆｅ ｖｒ ｎ ｎ ａ ｌ ｏ ｅ ｉ ｎ ． ｉ ａ ｕｔ
ｑ ａ ｉ ｔ ｅ ａ ｄ ｕ ｚ ｃ ｎｃｕ ｉｎ ． Ｂｙ ａ ｃ ｌｔｎ ｔ ｅ ｖ ｄａ ａ ｅ ｓ ｆ ｈ ｍａ ａ ｄ ｎｒ ｄ ｃ ｎ ｔ ｅ ＦＰ ｕ ｔ ｉ ｎ ｆ ｚ ｙ ｏ ｌ ｓｏ ｓ ｌ ａｖ ｃ ｌｕ ａｉ ｇ ｈ ｏ ｅ ｐ ｒ ａ ｏ ｔ ｅ ｐｓ ｎ ｉ ｔｏ ｕ ｉ ｇ ｈ Ｆ
Ａｂ ｔ ａ ｔ：Ｃｏ ａ ｉ ｇ ｔ ｅ ｓｍｉ ｒｔ ｆ ｖｒｕ ｌｇ ｏ ｈ ｍｉ ａ ｐｓ ｍａ ｅ ｆｏ ｔ ｅ ｄ ｔ ｆ ｄ ｆ ｒ ｎ ａ ｓ ｉ ｅ ｓｒ ｃ ｍｐ ｒｎ ｈ ｉ ｌ ｉｙ ｏ ｉｔａ ｅ ｃ ｅ ｃ ｌｍａ ｄ ｒ ｍ ｈ ａ ａ ｏ ｉ ｅ ｅ ｔｌ ｂ ｓ ａ ｎ ｗ ａ ｆ ｍｅｈ ｄ ｆｒ ｇ ｏ ｈ ｍｉａ ｎ ｌ ｔ ａ ｕ ｌ ｙ ｃ ｎｒ ｌａ ｄ ｅ ａｕ ｔ ｎ．Ｂｕ ｈ ｏ ｔ ｏ ｏ ｅ ｃ ｅ ｃ ｌａ ａｙ ｉ ｌｑ ａ ｉ ｏ ｔｏ ｎ ｖ ｌ ａｉ ｃ ｔ ｏ ｔｔ ｅ ｃ ｍｐａ ｉｏ ｙ ｅ ｅ ｉｈｔｃ ｕｄ ｇ ｔｏ ｌ ｒｓ ｎ ｂ ｙ ｓｇ ｏ ｌ ｅ ｎ ｙ
地球化学异常下限确定方法
地球化学异常下限确定方法一、地球化学数据处理基础数据处理的意义是获得较为准确的平均值(背景)和异常下限。
1、地球化学数据处理归根结底仍属于统计学的范畴，所以要求数据应是正态分布的，不是拿来数据就能应用的，特别是用公式计算时更要注意这一点。
正态(μ =0, δ =1)----(偏态)。
大数定理:又称大数法则、大数率。
在一个随机事件中，随着试验次数的增加，事件发生的频率趋于一个稳定值;同时，在对物理量的测量实践中，测定值的算术平均也具有稳定性。
所以如果在计算时，数据中包含较多的野值时，实际获得的是一个不具稳定性的算术平均，它实际不能替代背景值。
2、异常是一个相对概念，有不同尺度上的要求，所以不要将其看作一个定值。
在悉尼国际化探会议上(1976)，对异常下限定义:异常下限是地球化学工作者根据某种分析测试结果对样品所取定的一个数值，据此可以圈定能够识别出与矿化有关的异常。
并对异常下限提出了一个笼统的定义:凡能够划分出异常和非异常数据的数值即为异常下限。
据此，异常下限不能简单的理解为背景上限。
二、异常下限确定方法具体异常下限确定方法较多:地化剖面法、概率格纸法、直方图法、马氏距离法、单元素计算法、数据排序法、累积频率法……下面逐一介绍:1、地化剖面法:(可以不考虑野值)在已知区做地化剖面:要求剖面较长，穿过矿化区(含蚀变区)和正常地层(背景)，能区分含矿区和非矿区就可确定为下限。
2、概率格纸法:(可以不考虑野值)以含量和频率作图15%--负异常50%--背景值85%--X+δ(高背景)98%-- ( X+2δ)异常下限3、直方图法:(可以不考虑野值)能分解出后期叠加的值就为异常下限4、马氏距离法:(在计算时已考虑野值)针对样本，实际为建立在多元素正态分布基础之上—多重样本的正态分布，超出椭球体时—异常样(如P3点)。
相似于因子得分的计算，最后为一个剔除异常样本时的计算值，实际计算出综合异常边界线。
当令m=1时，上式化解为Xa=Xo?KS，这是我们较为熟悉的单元素(一维)计算异常下限常用公式。
地球化学图解应用
标准矿物岩石命名
CIPW标准矿物（Cross、Iddings、Pirrson、Washingdon，1903） 阳离子标准矿物（Barth-Niggli）
5
R1-R2岩石命名
另外，沉积岩、变质岩中矿物成分和化学 成分之间的简单关系十分困难，不能用简单 的图解来进行分类。
6
（2）协变图解
三变量图解
二变量图解
Harker图解：反映岩浆分离
结晶、部分熔融。沉积岩中
不同矿物组分混合的结果； 变质岩中的混合作用等等。
7
（3）成因系列判别
花岗岩K2O-SiO2图解
花岗岩SiO2-AR图解
花岗岩A/NK-A/KNC判别图
花岗岩SiO2-AR图解 8
2、微量元素图解
当地幔发生部分熔融作用时，微量元素优 先进入矿物相的元素称为相容元素；择优 进入熔体相的微量元素叫做不相容元素 （亲岩浆元素）。
酸性熔岩及玄武岩； （3）识别火山弧玄武岩
效果特别好。 注意：
样品不能含有大量的 蚀变玻璃及磁铁矿。
20
火山弧玄武岩优先使用的图解
La/10-Y/15Nb/8图解
Cr-Y图解
Cr-Ce/Sr图解
21
F1-F2-F3图解
成分范围：20wt%>CaO+MgO>12wt%
22
MgO-FeO-Al2O3图解
31
其他补充
U-Pb同位素：中生代及其以前的岩浆岩、变质岩、 沉积岩的沉积岩年龄、变质年龄、热事件年龄。
氢、氧、硫、碳同位素：计算成岩、成矿温度等物 理化学条件。
铅同位素：计算模式年龄，判别成因。
关于岩石成因系列：涉及到成因岩石学，不同的划 分方法和种类，内容繁琐。
地球化学图相关计算完整版
地球化学图相关计算 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】地球化学图相关计算一、计算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含量岩性样品号Al2O3含量CaO含量Na2O含量K2O含量(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/(CaO+Na2O+K2O)分子比（摩尔）Al2O3/(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)。
地球化学图相关计算集团标准化工作小组 [Q8QX9QT-X8QQB8Q8-NQ8QJ8-M8QMN]
地球化学图相关计算
一、计算Al2O3/(CaO+Na2O+K2O)和Al2O3/(Na2O+K2O)或A/CNK-A/NK分子(摩
尔)比步骤与投图
1．列出Al
2O
3
、CaO、Na
2
O和K
2
O的分子量
Al2O3分子量CaO分子量Na2O分子量K2O分子量主要氧化物分子量
氧化物Al
2O
3
CaO Na
2
O K
2
O Si
2
O TiO
2
Fe
2
O
3
FeO MnO MgO P
2
O
5
分子量
2．列出所测样品的Al
2O
3
、CaO、Na
2
O和K
2
O含量
岩性样品号Al
2O
3
含量CaO含量Na
2
O含量K
2
O含量
(xηγT
1)
P4XT06
(xηγT
1)
P15XT02
(zηγT
1)
P15XT03-
1
(zηγT
1)
P4XT10
3．求出分子(摩尔)数(样品中Al
2O
3
、CaO、Na
2
O和K
2
O含量除以各自氧化物的分
子量)
岩性样品号Al2O3分子数CaO分子数Na2O分子数K2O分子数
(xηγT
1)P4XT06
（=）（=）（=）（=）
(xηγT
1)P15XT02
（=）（=）（=）（=）
(zηγT
1)P15XT03-1
（=）（=）（=）（=）
(zηγT
1)P4XT10
（=）（=）（=）（=）
4．分别求出Al
2O
3
/(CaO+Na
2
O+K
2
O)和Al
2
O
3
/(Na
2
O+K
2
O)或A/CNK-A/NK分子 (摩尔)
比
岩性样品号Al2O3/(CaO+Na2O+K2O)分子比（摩尔）Al2O3/(Na2O+K2O)分子比（摩尔）(xηγT1)P4XT06++= +=
(xηγT1)P15XT02++= +==
(zηγT1)P15XT03-1++= +=
(zηγT1)P4XT10++= +=
5．投图
以Al
2O
3
/(CaO+Na
2
O+K
2
O) 或A/CNK分子（摩尔）比X为轴和以Al
2
O
3
/(Na
2
O+K
2
O)或
A/NK 为Y轴，投点即可。
换算公式 n = m / M，n是摩尔数，m是质量或含量，M是分子量。
摩尔比和分子比可以看做是一个意思。
二、微量元素蛛网图K、P、Ti等的计算
1. K=K
2
O××10000
2. P=P
2O
5
××10000
3. Ti=TiO
2
××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=（样品中SiO
2含量/SiO
2
分子重量）×阳离子个数=（）×=
Na=（样品中Na
2O含量/Na
2
O分子重量）×阳离子个数=（）×=
K=（样品中K
2O含量/K
2
O分子重量）×阳离子个数=（）×=
Fe=（样品中FeO含量/FeO分子重量）×阳离子个数=（）×=或（FeO+ Fe2O3）阳离子之和。
Ti=（样品中TiO
2含量/ TiO
2
分子重量）×阳离子个数=（）×=
Ca =（样品中CaO含量/ CaO分子重量）×阳离子个数=（）×= Mg =（样品中MgO含量/ MgO分子重量）×阳离子个数=（）×=
Al=（样品中Al
2O
3
含量/ Al
2
O
3
分子重量）×阳离子个数=（）×=
（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．求出R
2和R
1
阳离子数或原子数
将Si、Na、K、Fe、Ti各千阳离子数代入公式R1=4Si-11(Na+K)-2(Fe+Ti)，求
出R
1
阳离子数或原子数；同理，将Ca、Mg、Al千阳离子数代入公式R2=6Ca+2Mg+Al
求出R
2
阳离子数或原子数。
四、其它相关计算
1. Al=Al
2O
3×%
2. FeOt= FeO+ Fe
2O 3
3. 氧化指数 FeOt/MgO
4. Mg#（镁指数）
计算方式分为如下两步：
（1）FeOT (wt.%) = FeO (wt.%) + Fe2O3 (wt.%) * = FeO (wt.%) + Fe2O3 (wt.%) * 2))
（2）Mg#=100*(MgO//(MgO/+FeOT/
上述是分别测试分析了FeO和Fe2O3的计算方法，如果是测试的全铁，也可以近似计算。
通常说的高Mg，是指岩石具有较高的MgO含量，如火山岩中的高镁安山岩（通常情况下，异常高的MgO含量指示着可能有地幔物质参与，如俯冲带地幔楔或者软流圈熔体上涌等等）。
Mg#（镁指数）也可以定量的表示岩石中的Mg含量高低。
Mg#通常用于镁铁质岩石，可以粗略指示地幔岩石的部分熔融程度，高Mg#的地幔橄榄岩可能经历了更高程度的部分熔融，常在92-93左右，而原始地幔会相对富集，Mg#较低，在88-89左右。