基性、超基性岩地质数据库及库内数据计算程序(二)

三、岩石化学成分计算程序在研究岩石化学数据时,许多研究者从不同角度提出一些计算方法。我们根据基性、超基性岩石的特点,在数据库内配有A·H·查瓦里茨基、王恒升-白文吉、P·尼格里、A·里特曼、王一圣、谌义睿以及西安地质矿产研究所四室祁连山组等方法的计算程序。上述各种计算程序,每计算300个样品,仅需5分钟(包括计算结果打印时间)。下面仅介绍前两种方法的计算程序。     

TI—59计算器程序运算实例——褶积、滤波、圆滑、延拓

我们利用TI—59计算器编了一个褶积运算程序,如作单点运算,褶积因子点数最大为25点,在5点圆滑情况下也可启动,一次可连续算出多达41点的圆滑值。在生产中,单点运算使用起来比较方便,除计算第1点时要输入全部因子和与因子点数相同的信号数据外,以后每计算一点(不问因子长短),只需输入后继的一个信号数据,按一个标志键。由于滤波、圆滑,延拓等运算与褶积运算在数学上有着相同的形式,所以本程序可用于滤波、圆滑,延拓等多种数据处理。     

Rosenblueth法在土坡稳定可靠度分析中的应用

Rosenblueth方法是一种简单实用的可靠度指标计算方法 ,不需要了解各种状态变量的概率分布 ,只要利用它们的均值和方差 ,就可以求得边坡的可靠度指标 ,其计算满足精度要求。Rosenblueth方法计算所需的数据较少 ,利于计算机处理。论述了Rosenblueth方法的基本原理 ,用FORTRAN语言编制了计算程序 ,并通过算例比较证明了Rosenblueth方法的优点。在工程实际分析运用中 ,是一种值得推广应用的方法。     

Rietveld方法的R.A.Young程序的改进及其在矿物学中的应用

Ｒｉｅｔｅｖｌｄ方法是采用粉末衍射进行晶体结构修正及解析的一种方法。Ｒｉｅｔｖｅｌｄ方法的代表软件是Ｒ．Ａ．Ｙｏｕｎｇ程序。该程序给用户使用的主要困难是由于其输入控制文件（ＩＣＦ），是一个严格格式化的文本文件，其次是在运行时受到常规内存的限制。本文介绍了用ＭｉｃｒｏｓｏｆｔＣ７．０及ＭｉｃｒｏｓｏｆｔＳＤＫ３．１对Ｒ．Ａ．Ｙｏｕｎｇ程序的改进，（１）用对话框输入各种控制参数，形成初始化文件（ＩＮＩ），然后再将其转换为输入控制文件（ＩＣＦ）；（２）用Ｗｉｎｄｏｗｓ多文档界面（ＭＤＩ）实现图形显示，编辑和存储；（３）设计一软件接口，使在Ｗｉｎｄｏｗｓ下启动Ｒ．Ａ．Ｙｏｕｎｇ计算程序；（４）用ＦＯＲＴＲＡＮＰｏｗｅｒＳｔａｔｉｏｎ编译器对Ｒ．Ａ．Ｙｏｕｎｇ程序源代码进行重新编译，使其不受常规内存限制。用改进后的Ｒ．Ａ．Ｙｏｕｎｇ程序，对绿辉石的结构进行了修正     

一种定量确定断层位移性质和位移距离的新方法

介绍一种利用构造面和构造线的空间几何关系,运用矢量代数理论,通过数学计算定量地确定断层位移性质和位移距离的新方法。该方法克服了赤平投影和几何作图过程中的繁琐和误差,只要采用文中所附的电算程序,按要求输入数据便可迅速、准确地得到计算结果。     

York同位素等时线拟合方法的Matlab程序设计

等时线拟合方法的选择是提高同位素测年精度的关键因素之一。这里利用Maflab语言编写了基于York双误差回归模型的等时线拟合程序,包括isoacc．m、samw．m、errs．m和resi．m等模块,实现了读入数据,输入衰变常数和斜率迭代初始值,选择权重方法,以及计算截距、年龄、MSWD和绘图的功能。选用Brooks1972的数据,分别作三种权重等时线拟合,结果与Faure 1977的计算结果一致。程序简练,易使用,运算速度快,运行稳定。     

基于MATLAB的面积计算方法

由于在工程上只能得到一些离散数据，所以在面积的计算上就显得十分复杂，尤其当边界是不规则曲线时。这里提出一种基于MATLAB语言的简单、可行的计算方法，并对该方法的原理、算法、精度分析作了详细的阐述。通过模拟实验数据的计算和分析，说明该方法计算面积的精度较高。另外，所给程序还能绘出工程上需要的边界图形。     

计算机在铝土矿地质中的应用：地质图形量算微机数据处理系统

地质图形微机处理系统建立原始地质数据库是关键,各种计算方法都与数据程序相连。该系统对提高计算质量和效率起到保证作用。 地质数据库可用于各种储量计算以及地质数理统计。 程序采用BASIC语言编制,具有直观、通用、易学和会话性,有利于铝土矿地质工作者掌握。     

EL-5100型计算器用于野外磁法计算

EL-5100型袖珍式电子计算器,比一般计算器多十个记忆存贮器,可进行手控程序计算.计算过程中不必用任何算法语言,只要按照键盘上按键的功能,直接输入数学表达式和数据即可.输入的程序和输出结果都直观地显示在液     

微型计算机的岩石和矿物数据处理系统

中国科学院地质所岩石研究室经过多年实际运用,有下列主要岩石学和矿物学数据处理系统程序: 1.标准矿物分子(CIPW)计算程序:输入岩石化学分析数据,可计算28个标准矿物分子和39个常用的岩石化学特征参数,在程序中设计了几种Fe~(3+)和Fe~(2+)修正可选择开关, 2.矿物晶体化学式和某些矿物端员分子数计算程序:包括氧原子法和阳离子法,适用于硅酸盐矿物(含水矿物除外,但可计算含OH矿物的分子式)和氧化物。使用者可自由地进     

CS-3微型计算机配接光电机的几个问题

美国CROMEMCO公司生产的CS—3型微型计算机,不带光电纸带输入机,因此程序和数据都得从键盘输入。这样不但占用大量的机器时间,而且在程序比较长、数据量又大的情况下很不方便。为适应物化探数据处理工作的需要,提高机器的使用效率,就很有必要配接光电纸带输入机,用来输入程序和数据。我们在学习了清华大学配接小型光电机的经验后,于去年配接了一台大型快启停式光电纸带输入机;并在软件上实现了在BASIC和FORTRAN编译状态下,用程序中的输入语句由纸带输入数据。现把我们做的工作介绍如下:     

位置和尺度估计的计算机程序

近年来,稳健统计学在地质科学应用上引人注目。其方法在处理地质、化探观测数据时有着一定的发展前景。为了提高其应用的普遍性,笔者采用BASIC语言编写了位置和尺度估计的计算程序,供读者借鉴。一、计算公式位置和尺度估计的计算,据有关资料介绍方法很多,其计算结果大同小异。所以本文只介绍一种常用的极大似然估计法,简称M估计法,其计算公式如下:     

多种结构计算应用程序的比较

在结构设计的过程中，我们经常运用各种各样的结构计算程序，各种程序的结构分析模型的简化原则不同，模型的输入和输出结果的形式不同，计算的结果也不尽不同。我们所采用的程序是否能真正地反映出建筑物的实际受力状况呢?本文将对此进行浅要的比较。     

尼格里岩石化学计算程序(用FX-502P程序型计算器)

FX-502P是一种程序型计算器,进口价格约为人民币300元.它轻便,可随身携带;断电后能保留程序存贮;亦可利用程序输送器将程序装入录音磁带或从磁带调回计算器中.最近,我们编制了一个适用于该计算器的尼格里岩石化学计算程序.我们以这个程序利用FX-502P验算了《尼格里岩石化学计算的电算程序》一文(刊载于《地质与勘探》1981年第六期28～30页)提供的全部数据.平均每个样品的参数计算只用了三、五分钟.这样就免去了使用国产108机必须的纸带穿孔、检查、校正,程序上机编译,数据光电输入等等全部过程.节省了使用电算机的费用,更重     

计算岩层真厚度的最简单程序

当计算岩层真厚度时,人们力图摆脱某些条件的人工判别,以简化计算. 本文正是为这一目的,为TI-59型计算器编制的一个十分简便的程序.运算时只要输入有关变量,它就会自动运算出正确结果.这一程序仅用了58步,每次计算只需25秒钟.运算过程中,各数据适用性的判别及选择,由计算器自行完成.计算操作很简单,一般人都能掌握.     

硅灰石矿矿物成分及其含量计算的电子计算机方法

硅灰石矿床的质量评价,主要以矿石中有益有害矿物含量为依据。由矿石化学分析结果,计算硅灰石矿的矿物成分及其含量是一项极为繁杂的工作。为提高计算速度和精度,笔者根据张大山编写的河北省迁西小西营硅灰石矿的矿物含量计算方法,编制了 BD—200语言电算程序,并对该矿小批量的化学分析数据在6912—5机上计算结果证明,输出的各种数据可靠,精度大大提高。为推广电算技术,现将该程序及其使用方法做一介绍,供从事硅灰石矿物研究者参考。     

ClPW标准矿物电算程序

一前言岩石化学分析数据的换算、整理和对比的方法很多,换算工作量也较大.笔者就国内外常用的CIPW标准矿物计算方法,用ALGOL60语言编制了电算程序.该程序除计算并输出经典的CIPW标准矿物含量外,还同时计算并输出近十年来岩石化学领域常用到的其他一些化学成分特征参数,供岩系划分、岩石分类命名时参考.本程序尽量做到一次计算后能获得更多的数据.     

水准测量平差计算及程序设计

为了对地面沉降和地壳垂直形变进行研究,上海市地质处测量队曾对上海市一、二等水准测量的精度进行了估算。由于工作量比较大,因此后来改成用电子计算机进行计算,关于各种计算的结果已经作了总结。本文仅介绍在用电子机算机进行水准测量平差计算时的一些方法,对从事测量计算工作的同志,可能有所帮助。由于该项工作是在1974年完成的,而当时计算机的使用还不大普遍,加上作者对测量技术的无知,所以在程序设计时一些方法还是非常粗糙的。根据近年来发展起来的一些数学方法,对于大型测量网的数据信息安排还有更科学的方法以解决计算存量和速度方面的一些问题。     

CS—3机和108机数据纸带共享

为了使CS—3机和108机共同享用同一种数据纸带,我们设计了输入程序AGRD。一、程序功能能把108机上使用的ALGOL—60语言的数据纸带读到CS—3机中,转换成适于FORTRAN语言所使用的数据格式,存放在磁盘上,提供FORTRAN程序使用。本程序对纸带的穿孔格式没有任何特殊要求,凡是符合ALGOL—60语言所要求的五单位纸带,都可以使用。     

应用图形处理器快速计算逆时偏移

逆时偏移是目前精度最高的地震数据叠前深度偏移方法,但高强度的计算需求限制了其在工业生产领域的大规模应用。可编程图形处理器的发展为逆时偏移的快速计算提供了一种新的计算选择。围绕如何在图形处理器上开展逆时偏移计算展开,总结了图形处理器计算的优化关键,并根据逆时偏移的特点着重介绍了两个优化环节：一个是应用随机边界条件,以计算换存储,减少数据在主机和图形处理器间的传输;二是应用共享存储器来存储正演计算的波场,相比全局存储器,提高了数据读取的带宽。应用Marmousi模型数据对经过上述优化后的程序进行了测试,结果表明,图形处理器逆时偏移程序得到了很好的优化,提高了计算效率。