CIPW标准矿物及岩石化学参数计算
课题研究程序：包括研究方案、研究开题、实施研究和课题总结
1.岩石学特征
岩墙群的岩石结构和矿物成分较单一，除粒度由边部向中间由细变粗之外，其岩石结构和矿物成分没有什么变化，岩石类型为辉绿岩和辉绿玢岩，其中早侏罗世为以辉绿岩为主，中侏罗世和早白垩世以辉绿玢岩为主。
辉绿岩：岩石呈灰绿色、黑绿色，块状构造、辉绿结构，主要矿物成分为单斜辉石和斜长石，副矿物为磁铁矿、钛磁铁矿和磷灰石。岩墙边部的矿物粒度细，为0.1mm左右，中间部分粒度粗，为0.8～0.6mm。
辉绿玢岩：岩石呈灰绿色、墨绿色，斑状构造，基质为辉绿结构，斑晶为斜长石和单斜辉石。镜下见单斜辉石和斜长石斑晶穿插现象，说明它们几乎是同时结晶的。基质也主要由斜长石和单斜辉石组成，副矿物有磁铁矿、钛磁铁和磷灰石。斜长石斑晶粒度一般为2mm左右，最大可达5mm；单斜辉石斑晶粒度一般小于斜长石斑晶，基质斜长石粒度一般为0.1～0.4mm，单斜辉石粒度与基质斜长石相近。斑晶斜长石含量为5%～25%，斑晶单斜辉石含量为0～10%；基质斜长石含量35%～55%，基质单斜辉石含量35%～50%，其他少量。
闪长玢岩：岩石呈灰黑色，具斑状结构，斑晶为斜长石及少量角闪石。
2.矿物学特征
通过电子探针分析对矿物学组成特征描述如下。
（1）辉石
辉石是辉绿岩中最主要的造岩矿物，我们所研究的岩墙群中只有单斜辉石，未见斜方辉石。单斜辉石斑晶为半自形—自形，基质中为他形，常见简单双晶和聚片双晶。其电子探针分析结果及端元组分见表8-2。由表8-2可知，所研究的基性岩墙群的单斜辉石以含较高的MgO、CaO、P2O5为特征，多数样品质量分数 MgO ＞13%、CaO ＞18.5%、P2O5＞0.20%。另外，斑晶辉石比基质辉石富含Al2O3，可能表明斑晶结晶时的压力大于基质结晶。早白垩世岩墙中辉石MgO含量（14.57%～17.13%）高于早、中侏罗世岩墙的辉石（12.09%～15.99%）。单斜辉石的Wo-En-Fs关系图（图8-3）中亦反映出这个特点，即早、中侏罗世岩墙中的单斜辉石全部落在普通辉石区，而早白垩世岩墙的单斜辉石大部分落在顽（镁）透辉石区。早侏罗世辉绿岩墙的单斜辉石含 TiO2较高（1.99%～3.13%）。辉石含Ti高反映早期岩浆源区压力比较大，晚期辉石含Mg高反映源区岩浆富镁，来源深度大，表示地壳伸展程度和渗透程度增加。
表8-2 单斜辉石的电子探针分析结果
注：表8-2和表8-3的探针分析数据均由北京大学地球与空间科学学院舒桂明测定，*为辉石斑晶。
图8-3 单斜辉石的Wo-En-Fs关系图
Ⅰ—次透辉石；Ⅱ—顽（镁）透辉石；Ⅲ—普通辉石
1—早侏罗世辉石；2—中侏罗世辉石；3—早白垩世辉石；4—后期中基性岩脉中辉石
（2）斜长石
斜长石是辉绿岩/辉绿玢岩中最主要的造岩矿物，斑晶斜长为板状，长条状，自形—半自形晶，不论是斑晶还是基质，都发育有各类双晶，常见为聚片双晶，次为卡纳复合双晶，部分斜长石电子探针分析结果及其牌号见表8-3。从表中可知，早侏罗世岩墙之An为63-66之间。
通常岩墙代表一次岩浆快速侵位的产物，岩墙边部和中心的成分差异不甚明显，但是岩墙两侧率先结晶的冷凝边往往代表初始未演化的岩浆组分，如果是基性岩浆，其斜长石牌号有可能比岩墙中心部分略微偏大。斜长石成分的细微变化往往有规律地反映出岩石形成过程物质成分演变的特点。
表8-3中G1、G2、G3、G4 是从同一条1.34m宽的早侏罗世岩墙中依次采集的样品，G3位于中间，粒度最粗，G2、G4分别在两边，粒度中细，G1为最东侧的冷凝边最细（西侧冷凝边因被后期碱性岩脉侵入蚀变，未取样）。虽然它们对应的斜长石牌号相差甚小，但还是表现出从中间向两边变小的趋势，这与一般单脉岩墙正好相反（陈孝德等，1983），证实它可能不是一次侵位的单脉岩墙，后期有基性程度更高的岩浆补充。
中侏罗世辉绿玢岩的斜长石斑晶中出现环带构造，环带中心的斜长石偏基性，An为85，外围部分An为66，差异之大可能反映早期在偏基性的岩浆中结晶的长石被后期岩浆捕获。岳永君（1991）发现同一样品的斜长石斑晶和基质的牌号An分别为76和44，也反映后期岩浆偏中性。上述现象反映在岩浆演化过程中（或许在岩浆房内）有过多次基性程度不同的岩浆补给及相互作用。
早白垩世辉绿玢岩中斜长石钠质程度明显提高（4.52%），基质中斜长石牌号（An56）也提高，反映早白垩世岩浆更偏基性和碱性，表明岩浆来源不断加深。
表8-3 斜长石电子探针分析结果
3.岩石化学成分特征
本区的基性岩墙群SiO2含量变化于48.53%～52.64%，MgO变化于3.82%～8.45%，CaO变化于 6.34%～11.41%，Na2O变化于 1.58%～3.60%，K2O变化于 0.87%～1.50%，Na2O＞K2O。从CIPW标准矿物计算看，早、中侏罗世含标准矿物石英（Q），早白垩世含标准矿物橄榄石（Ol），从岩墙群的Mg′值［100×Mg/（Mg+Fe）］看，早白垩世Mg′值从早、中侏罗世的51～61增到67（表8-4）。在TAS图解（图8-4）中除一个样品投影于碱性系列玄武岩区以外，其余均位于拉斑系列玄武岩区。随着时间的推移，SiO2、Na2O+K2O含量总体降低。
表8-4 辉绿岩墙群的岩石化学分析（wB/%）结果及CIPW标准矿物计算结果
注：F3、F9、F10、Sh1由北京有色冶金设计研究总院中心化验室分析；G编号为中国科学院矿物资源中心X荧光法分析；C56和C12数据引自佘宏泉。
在R1-R2图解（图8-5）中，除了早侏罗世外，从晚三叠世、中侏罗世到早白垩世的岩墙群投点连续变化，从二长闪长岩经辉长岩演变到辉长苏长岩和橄榄辉长岩，岩墙基性程度逐渐提高，早侏罗世和早白垩世两期岩墙基性程度高，与岩浆来源加深有关。
与上述斜长石部分阐述的早侏罗世的复式岩墙的斜长石牌号变化规律一致，G1至G5为依次的全岩成分，具有对称分布的规律（表8-5），两侧的G1和G5 相近，G2和G4相近，G3显示了MgO、Na2O和Mg′值最高的特点，结合前述中间带G3号样品的斜长石牌号增高的特点，可以认为这种变化不是岩浆自身演化造成的，可能来自岩浆的多次补给与充填，而且这种再次充填和补给时间相隔不长，在接近的温度条件下，使得新生岩墙两侧不具冷凝边，只保留了最初侵位时形成的冷凝边。反之如果是新上涌的岩浆从两侧边部充填，则应有不止一对冷凝边存在。因此作者将其定为复式岩墙。这种从中心向两侧对称扩展的填充方式类似大洋中脊的扩张。早侏罗世和早白垩世复式岩墙的存在也表明这两个时期伸展作用最强。这种对称扩展的岩墙显示了深部岩浆上涌造成水平方向的伸展作用。
图8-4 辉绿岩墙TAS图解
1—早白垩世；2—中侏罗世；3—早侏罗世；4—晚三叠世
图8-5 R1-R2图解
图例同图8-4
表8-5 林西富林构早侏罗世复式岩墙的组成（wB/%）
注：由中国科学院地质地球物理所矿物资源中心测试，Mg′=100×Mg/（Mg+Fe）。
4.微量元素特征
表8-6 中列出了研究区辉绿岩岩墙群的微量元素分析结果。和前述 R1-R2图解（图8-5）中显示的结果类似，除了早侏罗世外，从晚三叠世、中侏罗世到早白垩世岩墙群的稀土元素总量（ΣREE）不断降低，早白垩世的ΣREE只有42.5×10-6，近于晚三叠世ΣREE的1/3。图8-6中也表现出晚三叠世、中侏罗世到早白垩世的稀土配分曲线斜率依次变小，早白垩世呈现近水平状，具有轻微的Eu负异常，早侏罗世曲线与之平行，只是呈平坦型。轻、重稀土的富集与否及其程度与岩石成因有关，晚三叠世稀土配分型式反映岩浆来源以地壳的局部熔融作用为主，因为局部熔融时轻稀土优先进入熔体，重稀土在部分熔融时易于保留在残余固相中。而早白垩世的配分型式表示有未分异的原始地幔岩浆的加入，受地壳混染程度较低。
表8-6 各时代岩墙群稀土元素分析（wB/10-6）
图8-6 各时代岩墙群的稀土配分曲线
（样号同表8-6）
从岩墙群岩石的微量元素（表8-6）和蛛网图（图8-7）看，晚三叠世与早、中侏罗世和早白垩世的配分型式明显不同，前者Th富集，Nb、Zr、Hf亏损，而后者Nb、Ta富集，Th相对亏损。Th的富集和Nb的亏损表明岩浆局部有陆壳特征或受到地壳混染。与此相反，Nb的富集和Th的亏损说明幔源岩浆受地壳混染程度很低或基本未受混染（Oliveira et al.，1990），裂谷初期玄武岩以Th的亏损为典型（李昌年，1992）。由此看来，晚三叠世岩墙群具有大陆壳的特征，而早、中侏罗世和早白垩世辉绿岩或辉绿玢岩来自地幔，受混染程度低，具有裂谷初期玄武岩的特征。
图8-7 各时代岩墙群的微量元素蛛网图
（样号同表8-6）
Thompson等认为La/Nb比是地壳混染的一个有用指数（转引自Wilson，1989），洋岛玄武岩、大陆碱性玄武岩和金伯利岩La/Nb＜1，而大陆溢流玄武岩La/Nb为0.5～7，从表8-7中可见除了晚三叠世C12的La/Nb＞1外，其他样品La/Nb＜1。结合部分微量元素比值（Hart et al.，1992）来看，赤峰地区晚三叠世闪长质岩墙群的比值与大陆地壳十分接近，而北部林西地区的侏罗-白垩纪的辉绿岩或辉绿玢岩岩墙群大部分与原始地幔的微量元素比值接近，个别元素之比，如Th/Nb显示了亏损地幔的特征，这与其他手段的分析结果一致。
表8-7 岩墙群的部分微量元素比值和原始地幔部分微量元素丰度
（据Hart et al.，1992）
5.Sr-Nd同位素特征
林西地区的两个采样地点位于大兴安岭主峰带，岩墙群侵位于古生代增生陆壳基底上，早-中侏罗世和早白垩世辉绿岩岩墙群的ISr比值分别为0.7036、0.7044和0.7048（表8-8），初始比值低，显示岩浆来源于上地幔，未受到明显的地壳混染，TDM模式年龄表明基底不存在古老陆壳。而位于华北克拉通北缘的赤峰柴胡栏子晚三叠世的辉长闪长岩ISr值为0.7072，显示了岩浆来源与地壳部分熔融关系密切。
表8-8 岩墙群的Sr-Nd同位素组成
注：均由中国科学院地质地球物理研究所质谱实验室乔广生、张任祜测试。
1.程序功能
程序CHECAL具有以下功能：①全岩化学分析数据的检查和修正，包括检查H2O+、CO2的含量和分析总量，检查和修正Fe2O3/FeO*比值，干成分换算等；②CIPW标准矿物计算；③常用的40余种岩石化学参数计算。
2.方法原理
对全岩化学分析数据的检查和修正，程序中未设定具体的标准，可由用户自行选择H2O+、CO2的上限和分析总量的上、下限。Le Bas等（1986）以w（H2O+）＜2%、w（CO2）＜0.5%作为判断新鲜火山岩的主要标准，用户在实际应用中可以此作为参考。全岩化学分析的总量，应限定在99.00%～101.00%之间。对于Fe2O3/FeO*比值的检查和Fe2O3、FeO含量的修正，程序中提供了8种可供选择的方法（马鸿文，1988，未发表资料：Irvine et al.，1971；Hughes，1982；Le Maitre，1976）。若化学分析结果中的挥发分含量过高（一般＞3%），则在进行CI'PW标准矿物和岩石化学参数计算前，应首先将其换算为干成分。
CIPW，际准矿物计算方案主要采用Holmes（1921）的算法。常用的岩石化学参数参考《岩浆岩岩石学》教材（邱家骧，1984）或其它岩石化学方面的参考教材或专著。
3.程序结构
结晶岩热力学软件
4.使用说明
（1）输入格式
程序运行过程中，按照屏幕提示，依次提供以下参数：
IFN 输入文件名
CO2 CO2的上限
H2O 2的上限
ITA 分析总量选择
Option Fe2O3/Fe*比值修正方法选择
Vfree 干成分换算选择
OFNl 输出文件名，保存CIPW标准矿物计算结果
OFN2 输出文件名，保存岩石化学参数计算结果
样品名称和氧化物原始分析结果由输入文件读入。
原始数据输入格式为：A6,15F6.2。
各变量排列顺序依次为：Sample（样品号）、SiO2、TiO2、Al2O3、Cr2O3、Fe2O3、FeO、MnO、MgO、CaO、Na2O、K2O、P2O5、CO2、H2O+、H2O-。
每次计算的样品个数不限。
（2）输出格式
原始数据的检查和修正结果输出到文件clboxd.dat中。若样品的CO2、H2O+或分析总量超过所设定的上、下限，则样品被自动剔除。
CIPW标准矿物计算结果以表格形式保存在输出文件中。CIPW标准矿物符号可参阅马鸿文（1993a）编著的《结晶岩热力学概论》附录三。除CIPW标准矿物含量外，输出结果中还包括分异指数（DI）、Di-Ab-An、An-Ab-Or、Q-Ab-Or三元系相图的端员组分含量、用于岩石化学分类命名的Q’-ANOR-F’参数（Streckeisen et al.,1979）及岩石的镁值（Mg#）。
常用岩石化学参数计算结果同样以表格形式保存在输出文件中。各参数的含义如下：
Totalk 总碱含量
IFels 长英指数
IAlk 碱质指数
IK 钾质指数
FeOt 全铁含量
Fe/Mg 铁镁比值
Cox 氧化率
IFe 铁质指数
IMg 镁质指数
SI 固结指数
ILars 拉森指数
Sigma 里特曼指数σ（Rittmann,1970）
Lgsgm 里特曼指数的对数lgσ
LgTau 戈梯尼指数对数lgτ
Lgsg25 Loffler指数的对数lgσ25（Loffler,1979）
F1,F2,F3 玄武岩构造环境判别函数（Pearce,1976）
Na-K-Ca 判别花岗岩成因的Na-K-Ca三角图（Raju et al.,1972）的端员组分含量
Fe-Alk-Mg FAM图解的端员组分含量
Fe-Mg-Al 判别玄武岩类构造环境的Fe*-MgO-Al2O3三角图（Pearce et al.,1977）的端员组分含量
Ti-Mn-P 判别玄武岩类构造环境的TiO2-MnO-P2O5三角图（Mullen,1983）的端员组分含量
Ti-K-P 判别玄武岩类构造环境的TiO2-K2O-P2O5三角图（Pearce et al.,1975）的端员组分含量
Ti-Al-K 判别玄武岩类构造环境的TiO2-Al2O3-K2O三角图（赵崇贺，1989）的端员组分含量
Mg# 镁值
S/Iidx S/I指数（Griffiths et al.,1983）
（Al-NaK）-Ca-（Fe2+Mg） 区分花岗岩成因类型的（Al-Na-K）-Ca-（Fe2++Mg）三角图（Hine et al.,1978）的端员组分含量
R1,R2 火山岩化学分类命名的R1-R2参数（De La Roche et al.,1980）
上述CIPW标准矿物含量和常用岩石化学参数的计算结果，按照每10个样品以相同的列表形式，顺序输出到相应的文件中。
5.程序文本
结晶岩热力学软件
结晶岩热力学软件
结晶岩热力学软件
结晶岩热力学软件
结晶岩热力学软件
结晶岩热力学软件
结晶岩热力学软件
结晶岩热力学软件
end if
****calculate total silica used and perform necessary desilications****
TSi=cnm（i,1）+cnm（i,1 2）+6*（cnm（i,18）+cnm（i,19）+2*（cnm（i,23）
$ +cnm（i,29）+cnm（i,24）+cnm（i,25）+4*cnm（i,27）+cnm（i,28）+cnm（i,30）
if（TSi.le.oxd（1）then
cnm（i,17）=oxd（1）-TSi
goto 300
else
DSi=TSi-oxd（1）
end if
if（DSi.le.cnm（i,30）/2）then
cnm（i,31）=DSi
cnm（i,30）=cnm（i,30）-2*DSi
goto 300
else
cnm（i,31）=cnm（i,30）/2
cnm（i,30）=0
DSi=DSi-cnm（i,31）
end if
if（DSi.le.cnm（i,11））then
cnm（i,11）=cnm（i,11）-DSi
cnm（i,12）=cnm（i,12）+DSi
goto 300
else
DSi=DSi-cnm（i,11）
cnm（i,12）=cnm（i,12）+cnm（i,11）
cnm（i,11）=0
end if
if（DSi/4.le.cnm（i,19）then
cnm（i,21）=DSi/4
cnm（i,19）=cnm（i,19）-DSi/4
goto 300
else
cnm（i,21）=cnm（i,19）
cnm（i,19）=0
DSi=DSi-4*cnm（i,21）
end if
if（DSi/2.le.cnm（i,18）then
cnm（i,20）=DSi/2
cnm（i,18）=cnm（i,18）-DSi/2
goto 300
else
cnm（i,20）=cnm（i,18）
cnm（i,18）=0
DSi=DSi-2*cnm（i,20）
end if
if（DSi.le.cnm（i,28）/2）then
cnm（i,26）=DSi
cnm（i,28）=cnm（i,28）-2*DSi
goto 300
else
cnm（i,26）=cnm（i,28）/2
cnm（i,28）=0
DSi=DSi-cnm（i,26）
end if
if（DSi.le.cnm（i,29）then
cnm（i,26）=cnm（i,26）+DSi/2
cnm（i,31）=cnm（i,31）+DSi/2
cnm（i,29）=cnm（i,29）-DSi
else
cnm（i,26）=cnm（i,26）+cnm（i,29）/2
cnm（i,31）=cnm（i,31）+cnm（i,29）/2
DSi=DSi-cnm（i,29）
cnm（i,29）=0
cnm（i,22）=DSi/2
cnm（i,20）=cnm（i,20）-DSi/2
end if
if（cnm（i,20）.lt.0）then
write（*,280）Sample（i）
280 format（/5X，’-lc norm for sample’，A6,
$ ’；poor analyses or sampling！’）
end if
**** calCulate formula weight of mafic norms ****************************
300 RMF=oxd（8）/TMF
cwt（29）=248.0944-RMF*31.542
cwt（30）=131.9307-RMF*31.542
cwt（31）=203.7771-RMF*63.084
sum=0
do j=1,31
cnm（i,j）=cwt（j）*cnm（i,j）
sum=sum+cnm（i,j）
end do
cnm（i,32）=sum
di=0
do j=17,22
di=di+cnm（i,j）
end dO
cnm（i,33）=di
sum=0.01*（cnm（i,29）+cnm（i,19）+cnm（i,23）+1.0E-8
cnm（i,34）=cnm（i,29）/sum
cnm（i,35）=cnm（i,19）/sum
cnm（i,36）=cnm（i,23）/sum
sum=0.01*（cnm（i,23）+cnm（i,19）+cnm（i,18）+1.OE-8
cnm（i,37）=cnm（i,23）/sum
cnm（i,38）=cnm（i,19）/sum
cnm（i,39）=cnm（i,18）/sum
Fld=0.01*（cnm（i,20）+cnm（i,21）+cnm（i,22）+1.0E-8
cnm（i,45）=Fld/（Fld+sum）
sum=0.01*（cnm（i,17）+cnm（i,19）+cnm（i,18）+1.0E-8
cnm（i,40）=cnm（i,17）/sum
cnm（i,41）=cnm（i,19）/sum
cnm（i,42）=cnm（i,18）/sum
cnm（i,43）=cnm（i,17）/（sum+0.01*cnm（i,23）
cnm（i,44）=100*cnm（i,23）/（cnm（i,18）+cnm（i,23）+1.0E-8）
if（i.lt.n）goto 15
500 If（ns.lt.n）i=i-1
write（3,530）（Sample（k）,k=1，i）.
530 format（/1X，’Title’，10A7）
531 do j=1,m
do k=1,i
if（OXDS（k,j）.gt.0.005）then
write（3,535）OXT（j）,（OXDS（kk,j）,kk=1,i）
exit
end if
end do
end do
write（3,534）（Tot（k）,k=1,i）
结晶岩热力学软件
45 format（A6,15F6.2）
ns=ns+1
SiO2O=OXD（1）
Ti02O=OXD（2）
Al2O3=OXD（3）
Fe2O3=OXD（5）
FeO=OXD（6）
MnO=OXD（7）
MgO=OXD（8）
CaO=OXD（9）
Na2O=OXD（10）
K2O=OXD（11）
P205=OXD（12）
if（MgO.eq.0）MgO=1.0E-8
TA=Na2O+K2O
if（TA.eq.0）TA=1.0E-8
PRM（i, 1）=TA
CTA=CaO+TA
PRM（i,2）=100*TA/CTA
TFeO=Fe2O3+FeO
TFM=TFeO+MgO
PRM（i,3）=100*TA/TFM
PRM（i,4）=100*K2O/TA
FeOt=FeO+0.8998*Fe2O3
PRM（i,5）=FeOt
PRM（i,6）=FeOt/MgO
PRM（i,7）=Fe2O3/TFeO
PRM（i,8）=100*TFeO/TFM
PRM（i,9）=100*MgO/TFM
PRM（i,10）=100*MgO/（TFM+TA）
PRM（i,11）=SiO2O/3+K2O-（CaO+MgO+FeO）
if（SiO2O.gt.43.0）then
Sigma=TA**2/（SiO2O-43）
PRM（i,12）=Sigma
else
PRM（i,12）=-1
end if
PRM（i,13）=Sigma
PRM（i,14）=（Al2O3-Na2O）/（Ti02O+1.0E-8）
PRM（i,15）=100*TA**2/（SiO2O-25）
do j=13,15
if（PRM（i,j）.gt.0）then
PRM（i,j）=-1
else
PRM（i,j）=Log10（PRM（i,j）
end if
end do
PRM（i,16）=0.0088*SiO2O-0.0774*Ti02O+0.0102*Al2O3+0.0066*FeO
$ -0.0017*MgO-0.0143*CaO-0.0155*Na2O-0.0007*K2O
PRM（i,17）=-0.0130*SiO2O-0.0185*Ti02O-0.0129*Al2O3-0.0134*FeO
$ -0.0300*MgO-0.0204*CaO-0.0481*Na2O+0.0715*K2O
PRM（i,18）=-0.0221*SiO2O-0.0532*Ti02O-0.0361*Al2O3-0.0016*FeO
$ -0.0310*MgO-0.0237*CaO-0.0614*Na2O-0.0289*K2O
PRM（i,19）=100*Na2O/CTA
PRM（i,20）=100*K2O/CTA
PRM（i,21）=100*CaO/CTA
TFMA=TFM+TA
PRM（i,22）=100*TFeO/TFMA
PRM（i,23）=100*TA/TFMA
PRM（i,24）=100*MgO/TFMA
TAFM=FeOt+MgO+Al2O3
PRM（i,25）=100*FeOt/TAFM
PRM（i,26）=100*MgO/TAFM
PRM（i,27）=100*Al2O3/TAFM
TMP=Ti02O/10+MnO+P205
if（TMP.eq.0）TMP=l.0E-8
PRM（i,28）=10*Ti02O/TMP
PRM（i,29）=100*MnO/TMP
PRM（i,30）=100*P205/TMP
TKP=Ti02O+K2O+P205
if（TKP.eq.0）TKP=1.0E-8
PRM（i,31）=100*Ti02O/TKP
PRM（i,32）=100*K2O/TKP
PRM（i,33）=100*P205/TKP
TAK=Ti02O+Al2O3/10+K2O
PRM（i,34）=100*Ti02O/TAK
PRM（i,35）=10*Al2O3/TAK
PRM（i,36）=100*K2O/TAK
do j=1,m
OXD（j）=OXD（j）/GFW（j）
end do
PRM（i,37）=OXD（8）/（OXD（6）+OXD（8））
PRM（i,38）=OXD（3）/（OXD（9）+OXD（10）+OXD（11））
AlNaK=2*（OXD（3）-OXD（10）-OXD（11））
Sum=AlNaK+OXD（9）+（OXD（6）+OXD（8））
PRM（i,39）=100*AlNaK/Sum
PRM（i,40）=100*OXD（9）/Sum
PRM（i,41）=100*（OXD（6）+OXD（8））/Sum
PRM（i,42）=1000*（4*OXD（1）-22*（OXD（10）+OXD（11））-2*（OXD（2）+2*OXD（5）
$ +OXD（6）））
PRM（i,43）=1000*（6*OXD（9）+2*（OXD（8）+OXD（3））
if（i.lt.n）goto 40
200 if（ns.lt.n）i=i-1
write（4,215）（Sample（k）,k=1,i）
215 format（/lX，'paramt'，10A7）
do j=1,npr-2
write（4,218）Index（j）,（PRM（k,j）,k=1,i）
218 format（1X,A6,10F7.2）
end do
do j=npr-1,npr
write（4,220）Index（j）,（PRM（k,j）,k=1,i）
220 format（1X,A6,10F7.0）
end do
if（i.eq.n）goto 35
goto 236
230 write（*,*）'File read error,data skipped！’
236 write（*,240）OFN
240 format（/'Edit'，A，'to look over the results！'）
end
*************************************************************************
6.计算实例
甘肃北祁连九个泉地区，早古生代拉斑玄武质海相火山岩的CIPW标准矿物含量和常用岩石化学参数计算（马鸿文等，1994，未发表资料）。
输入文件：exam83.dat
83105 48.87 1.97 15.72 4.50 7.17 0.22 6.61 4.41 5.10 0.46 0.33 0.31 3.72 0.37
83106 50.45 1.67 14.98 7.51 4.63 0.19 4.83 6.58 5.21 0.10 0.27 0.26 2.98 0.29
90410 49.95 1.10 15.93 1.83 7.33 0.18 6.57 8.68 3.66 0.42 0.12 0.35 3.56 0.10
72404 50.36 2.19 12.91 3.47 10.39 0.21 6.45 5.46 3.76 0.64 0.37 0.00 3.24 0.49
71913 50.97 1.46 15.34 4.13 5.41 0.15 6.42 6.54 4.82 0.09 0.14 0.91 2.83 0.56
71607 48.53 1.19 15.51 2.51 7.00 0.22 7.12 8.62 3.86 0.20 0.07 0.00 4.34 0.50
71902 48.82 1.58 15.40 3.40 5.00 0.16 6.63 7.17 4.10 0.84 0.18 2.67 3.39 0.53
72409 45.83 2.60 13.06 7.62 8.38 0.20 7.36 7.37 3.32 0.26 0.31 0.00 3.57 0.10
71509 51.60 1.06 15.36 2.68 6.15 0.34 7.03 6.26 4.83 0.04 0.06 0.00 3.74 0.49
71604 47.24 0.94 18.73 2.02 5.51 0.15 7.05 10.31 2.15 1.63 0.09 1.12 3.06 0.29
输出文件：exam84.dat
结晶岩热力学软件
结晶岩热力学软件
结晶岩热力学软件
龙先醋酸： CIPW标准矿物计算1. 功能简介本程序的CIPW标准矿物计算主要用于火成岩 CIPW 标准矿物计算,适用于基性岩到酸性岩各类岩石.为方便作图(与散点图程序配套),该程序部分还附带有两组岩石化学参数计算,即 Batchelor 的多阳离子R1-R...
子长县15064932089： 怎么用软件计算岩石的CIPW标准矿物?geokit之类的软件处理的,非常感谢 - ？
龙先醋酸： 安装geokit最新版本,然后按照软件中各元素的顺序将自己取样的数据复制进去,点计算就行了,非常方便.
子长县15064932089： (二)拉斑玄武岩 - ？
龙先醋酸： 拉斑玄武岩(ic basalt)并非具有拉玄结构,其主要特点是化学成分中w(SiO2)相对较高,而碱含量较低,w(Na2O+K2O)多为2%~4%,为亚碱性玄武岩的典型代表.主要矿物成分为斜长石(斑晶多为拉-倍长石,基质为拉-中长石)、普通辉石和...
子长县15064932089： 请教:化学全分析结果转换为矿物含量？
龙先醋酸： 你是要求什么矿物含量,有CIPW,尼格里矿物含量,或者是花岗岩,变质岩等,要有具体的说明才能用不同的方法计算呀.
子长县15064932089： 火山岩的分类和命名 - ？
龙先醋酸： 按火山作用方式的不同可分为熔岩和火山碎屑岩两大类. 7.3.1.1 熔岩的分类和命名 熔岩的分类有矿物定量分类、化学定量分类和定性分类等方案,三者常互相配合使用.对于全晶质火山岩,可使用矿物定量分类;对于不能准确识别矿物成分及...
子长县15064932089： 钻井工程技术手册 钻井手册 哪个好 - ？
龙先醋酸： 钻井工程技术手册 本手册是针对石油钻井工程需要而编写的,系统介绍了涉及石油钻井工程的工艺技术、钻井工具、技术要点等内容.主要章节包括钻井基础数据及有关计算、钻井设备、钻井管材、钻井工具、随钻震击器与减震器、井下动力钻...
子长县15064932089： 玄武岩的主要岩石类型有哪些? - ？
龙先醋酸： 玄武岩类的岩石类型较多,但自然界中常见的类型主要为玄武岩、拉斑玄武岩、碱性玄武岩、粗玄岩、玻基玄武岩和细碧岩.其中玄武岩可以认为是这一类火山熔岩的代表,特别是在没有化学分析资料、又不具备碱性玄武岩、拉斑玄武岩等的显...
子长县15064932089： 岩石化学计算方法？
龙先醋酸： ① 火成岩 .来自地球内部的熔融物质,在不同地质条件下冷凝固结而成的岩石.当熔浆由火山通道喷溢出地表凝固形成的岩石,称喷出岩或称火山岩.常见的火山岩有 玄武岩 、 安山岩 和 流纹岩 等.当熔岩上升未达地表而在地壳一定深度凝结...
子长县15064932089： 铁矿石的抗压等级是多少?如何计算矿石的岩石硬度等级 - ？
龙先醋酸： 1.与盐酸反应,看他生成多少氢气,2.通过Fe+2HCl=H2+FeCl2 算出铁的质量,3.再除以铁矿石的质量.对于矿石的岩石硬度等级计算一般用;普氏分级法
子长县15064932089： 查氏系数的概念是什么?有什么意义?？
龙先醋酸： 不好意思,从来没有听说过查氏系数,各种文献里也找不到,请问是不是打字有误? 若指硬度,那当然就是普氏系数. 普氏硬度表示矿岩的坚固性的量化指标. 人们在长期的实践中认识到,有些岩石不容易破坏,有一些则难于破碎.难于破碎的...
