对A·H·扎瓦里茨基岩石化学计算法的探讨

苏联著名岩石学家A·H·扎瓦里茨基提出的岩石化学计算法,在苏联地质界得到了广泛运用和好评。随着A·H·扎瓦里茨基著作的翻译出版,扎氏法在我国也得到了广泛运用。多年来,在大量的岩石学文献中,在我国一些高等院校编写的《岩浆岩岩石学》教材中,都介绍和运用了这个方法。在近几年的岩石学论文中,也大都附有扎氏数值特征,所以,对扎氏法作一些探讨,还是有必要的。     

用多边形截面法评价罗河矿区重力异常

一、引言对已知勘探矿区,在物探工作已完成的情况下,作任意形体重磁异常的正演计算,目前常采用的电算方法有直立线元法、正方形体元法、表而积分法及多边形截面法等。每一种方法在适当勘探条件下,都有其一定的应用效果。但我们认为,中国地质大学陈钟琦老师编制的多边形截面法电算程序,在形体逼近、勘探资料充分利用及对重磁异常正演计算精度上都较其它方法更适合于安微庐江矿区进一步的评价工作。为计算重磁的剩余异常,解决深部盲矿体存在与否,我们采用了多边形截面法程序Ⅰ、Ⅱ进行了资料处理。计算结果表明:只要合理选用野外资料,恰当地对数据进行预处理,采用多边形     

格里戈良分带指数法的改良

以烟筒砬子铜金矿区ZK001元素垂向分带序列的计算为例,对格里戈良分带指数法进行了两个方面的改良:在原始数据换算中,采用正规化法将原始数据变换为[0,1]区间的数值,避免了传统方法的人为因素,减少了计算量;在计算分带指数的变化梯度时采用原格氏法变化度的倒数,解决了原格氏法计算中可能出现的某一层位的分带指数为零则无法计算的问题。新旧两种格氏法计算结果相比,新方法更符合实际地质情况。     

降水自记曲线坐标读图方法及程序框图

现行自记降水电算整编程序的数据整理,均为时段降水量。此法沿袭了人工整编的习惯,对电算而言增加了脑力劳动的强度,工效低,易出错,校对不方便。我们通过实践采用坐标读图法,获得了电算数据整理省事、省时、直观便于校对的效果,提高了数据整理的质量。坐标读图法是在降水自记记录纸上按降水强度变化且不跨8时的原则摘录转折点,摘点以纵横坐标定点。即横坐标的读数是时间,纵坐标的读数是累积降     

岩石化学换算和编码的电算方法及其应用

一、前言岩石化学换算及编码是区调、科研方面应用很广的一项基础工作。在标准矿物换算上,国内、外比较流行的算法和电算程式多采用C·I·P·W法,但里特曼法建立在统计或然性处理基础上,其标准矿物与岩石实际矿物可比性强,是目前岩石化学换算法中最好的一种,但计算步骤繁杂而难以推广。本文的电算方法系笔者将里特曼法(A·Rittmann,1973)中大量图表进行数学处理后提出的一种程序,经与原著12个算例138个标准矿物含量比较,平均相对误差为4.73%,这种误差包括中间值取整数及运算位数不同在内。一般说本程序更为精确     

金属矿电法电算技术方法研究选问论证会报导

为加强开展金属电法电算工作,考虑安排81年电算科研任务,部科技局、物化探局于80年11月28—29日在京主持召开了“金属矿电法电算技术方法”选题论证会。参加会议的有武汉、长春、成都地院、计算所、物探所、矿床所,内蒙、吉林、安徽和湖北地质局物探队的同志,在会     

拉格朗日插值公式病态现象的研究

拉格朗日插值公式病态现象的研究郭永杰（长春水文水资源勘测局）拉格朗日插值公式（以下简称拉氏插值公式）是比较常用的数值插算公式。在水文资料电算整编中，水位流量关系线法推流部分大多采用拉氏插值公式的二次式形式，即一元三点插值法。经过近十年的实际应用，采用...     

球形矿体上航电异常的近似计算方法

引言 航空电法工作之前或在资料解释过程中,往往需要估计待测矿体可能引起的异常幅度、宽度以及它们与飞行高度等因素的关系。板状导体已有模型实验结果可供估计参考。均匀大地也有电算结果可以利用。关于球形导体,已有了用精确公式电算的直立线圈共轴硬架式频域连续波航空电磁响应剖面曲线。但我们可以看到,既使对不导电围岩中的球形导体而言,磁偶极发射源航电异常的精确计算结果也必须采用大型电算机才能     

灰色系统GM(n,m)模型参数的递推估计和时变估计

本文给出灰色系统GM(n,m)模型参数的递推最小二乘估计和时变参数估计,前者的优点是可以避免反复求逆矩阵,因而使计算简便,同时克服了文献1中给出的估计方法当观察数据多、矩阵线数高、分配单元增加时,微型机容量不够的弱点,而且随着新数据的到来,通过递推算法自动对系统参数加以修正,提高了模型的自适应性;时变参数估计可以缩小系统的时变性与模型参数非时变性所产生的差异。     

《电测深资料的一种快速电算解释方法》通过部级评审

—九八二年十二月十二日至十五日,在地质矿产部物化探局主持下,于四川省成都市对我院物探系承担的部属重点科研项目《金属矿电法勘探电算技术研究》的阶段研究成果《电测深资料的一种快速电算解释方法》等进行了评审。该课题研究以确定基岩深度为主要目的的电测深快速电算解释方法和程序。解释的基本方案是平均电阻率法。在研究中发现和利用了G型曲线的新性质,提出了新的解释方法,将传统的平均电阻率法的应用范围扩大到尾段     

有限单元法地下水水质模型的微型机求解

近几年,我们在山东济宁地区开展了区域的溶质运移问题的试验研究工作。本文根据已建立的济宁地区的二维地下水水质模型,利用有限元法的水质预测的微型机(VICTOR-9000)程序,处理了该区东南角面积为3.7平方公里的计算区的实际资料,选定铬为污染物质,在微型机上模拟计算了地下水铬浓度,计算结果基本达到精度标准,计算的铬浓度分布基本上符合实际铬浓度的变化趋势。并在微机上优选确定了计算区的弥散参数——弥散度,这将为该区地下水污染的预测和控制提供了可靠的依据。     

水文数据分散录入的联机系统

全国水文部门每年要刊印10卷74册约4000万个数据的水文年鉴,这些数据是从约2亿多个原始观测数据中经整理计算而得的。近年来,由于采用了电算整编,大大地减轻了工作强度,提高了计算精度和工效。但电算整编的原始数据,一般都要集中到水文总站,人工从终端键盘录入计算机,工作量很大。尤其目前正在建立水文基本数据库,需要集中录入历年的数据,工作量更大。在目前未普遍采用固态存贮器和     

微型机EG—3003在X射线荧光光谱测试化探样品中的应用——八个元素测定的数据处理

在X光射线荧光谱分析化探样品中,我们运用微型机EG—3003对八个元素测定的数据进行了处理。在本程序中对分析样品的基体校正、背景的扣除、谱线干扰的校正、归系数的求取、统计参数的给出、打印分析结果等,这一系列数据处理都是由微型机执行编制好的同一程序完成的。本程序由BASIC语言编写的。通过二千多个化探样品测试数据的处理说明,运用微型计算机进行数据处理大大加快了化探样品的测试速度,摆脱了人工的烦锁的数据计算,提高了样品测试的精确度、准确度,检出限。程序通过打印机给出统计参数,使分析者对每批样品测定质量都能进行监查。数据的输入以及荧光强度的采集,本程序运用“READ—DATA”语句,使计算     

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

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

洪水期水位流量关系绳套曲线的直接拟合

受洪水涨落影响的水痊流量关系通常为时序型绳套曲线,现有的校正因数法在电算整编中存在一定的困难,通过对校正因数法公式的幂级数展开,提出了水位流量关系绳套曲线的计算机直接拟合数学方法,并对该法原理、步骤作了较详细的介绍。实例验证结果表明,该法能有效地模拟受洪水涨落影响的水位流量关系绳套曲线、且具有计算速度快,推流精度高等优点。     

直立共面航电系统在球体上的电磁响应

本文讨论直立共面航电系统(HDY-202)在球体上的电磁响应。它可表为一个级数,其中每一项表示一种复极的二次磁场。通过电算给出磁偶极场中球体的响应函数值及实际解释用的实、虚分量复平面图。当围岩为高阻时,这种图是有实际意义的。     

在重磁异常解释中应用电算资料的一些问题

随着电算技术的普及,目前都把重磁资料的电算处理(如化极,求导,上下延拓等)作为异常解释的先行步骤.这样作的结果,一方面深化了对重磁场特征的认识,另一方面也为定性、定量解释提供了一些有效手段,对提高重磁异常推断水平也是有益的.但在近几年的实践中,通过电算与实测资料的对比,也发现了一些具体问题值得引起重视,现在提出来供同志们参考,不当之处请予指正.     

压实黄土非线性压缩应力-应变关系在地基沉降计算中的应用

目前以压缩试验得到的e-p曲线为基础的分层总和法仍然是地基沉降分析中最常用的方法。割线模量法是近年来为人们所熟知的计算地基沉降的新方法,具有不受初始孔隙比影响、便于电算等优点。传统的割线模量法将土体的压缩应力-应变关系视为双曲线形式。通过对压实黄土压缩应力-应变关系曲线的研究和对传统割线模量法计算的割线模量与压缩模量的比较表明,导致割线模量法和e-p曲线法计算结果差异的主要原因是土体压缩应力-应变关系符合双曲线的假定,将压实黄土的应力-应变关系用更符合实际情况的幂函数进行拟合,得到了基于幂函数形式的压缩应力-应变关系曲线的沉降量计算公式,并结合具体工程进行了沉降计算对比分析,将割线模量法进行了推广应用。     

频率-波数领域震波模型研究

虽然傅氏转换能够筒化波动方程式的解法，但在频率-波数领域中，求解波动方程式，计算震波模型，该方法不能适应传波介质之界面弯斜或传波速度横向变化。把界面视为一组不连续的绕射点，用各点反射系数作为讯号源强度，各绕射点分别在频率-波数领域中实施波场向上外插，再把各点外插波场相加，由此得到震波模型。利用相位转移函数，取代每一绕射点的傅氏转换，能节省大量电算时间。这种频率-波数领域震波模型法，比时间-空间领域法更容易处理复杂介质模型。用新的方法对几种典型之地质模型，包括断层、背斜、尖灭、多层簿地层岩相突变等进行了计算，所得震波特性相当合理，符合震波传播理论。实例表明，即使对于岩相突变，计算的震波模型与实际震波记录，亦相当吻合。频率-波数领域震波模型法，对于复杂界面和速度变化之地质模型，其数学处理上远比射线理论和克希霍夫积分方程来得可行且单纯。且比定差法和定元法显著节省网格点和容易处理边界条件。唯这种频率-波数领域绕射点法，随绕射点数之增加，其电算时间之增加幅度似乎较高。关键词震波模型傅氏转换波动方程式频率-波数领域波场外插相位转移函数     

应用格氏法确定元素分带序列具体问题的探讨

格氏法的实质是利用元素的最大聚集重心的相对空间位置确定元素分带序列的,是通过计算某一元素的线金属量与该剖面所研究的几种指示元素线金属量之和(总金属量)的比值得以实现的.该方法在向定量化判断矿体剥蚀深度、寻找盲矿体及抑制偶然误差等方面有独到处,可在生产研究中试用.因该法已在几种书刊上作过专门介绍,在此从略.下边谨以确定元素垂直分带为例进行具体方法探讨: 一、几个元素的分带指数最大值同时位于相对最上部或最下部中段时,采用的公式完全相同,都是计算变化指数: