ρ（L_（r。w））的界和迭代法收敛性

ρ（Ｌ_（ｒ。ｗ））的界和迭代法收敛性刘光清，辛迪远（长春地质学院基础科学系．长春１３００２６）（四平电视大学．四平１３６０００）本文给出了广义ＡＯＲ法（广义ＳＳＯＲ，广义ＳＯＲ，广义Ｇ－Ｓ及广义西法）迭代矩阵谱半径较理想的上界及其收敛准则，推广了文?..     

并行布尔模拟算法

布尔模拟是一种基于目标的模拟方法,主要用于模拟河道、泥岩夹层等具有相对简单几何形状对象的分布。地质统计学软件库（GSLIB）中的Ellipsim程序是布尔模拟的经典程序,采用Fortran编写。笔者主要介绍了如何利用VisualC＃4．0和．NETFramework4．0的并行扩展把Ellipsim从串行化模式改写为并行化模式。通过自动测试模块和耗时统计模块统计了并行化的加速比,结果表明并行后的程序在很大程度上提高了计算速度。布尔模拟的并行化方法对于其他随机模拟程序的改写具有很好的借鉴意义。     

加速型改进迭代法在非饱和土渗流中的应用研究

Richards方程常用于非饱和土渗流问题,并且应用广泛。在数值求解中,对Richards方程线性化,进而采用有限差分法进行数值离散以及迭代计算。其中传统的迭代法比如Jacobi迭代、Gauss-Seidel迭代法（GS）和连续超松驰迭代法（successive over-relaxation method,简称SOR）迭代收敛率较慢,尤其在离散空间步长较小以及离散时间步长较大时。因此,采用整体校正法以及多步预处理法对传统迭代法进行改进,提出一种基于整体校正法的多步预处理Gauss-Seidel迭代法（improved Gauss-Seidel iterative method with multistep preconditioner based on the integral correction method,简称ICMP(m)-GS）求解Richards方程导出的线性方程组。通过非饱和渗流算例,并与传统迭代法和解析解对比,对改进算法的收敛率和加速效果进行了验证。结果表明,提出的ICMP(m)-GS可以很大程度地改善线性方程组的病态性,相较于常规方法GS,SOR以及单一改进方法,ICMP(m)-GS具有更快的收敛率,更高的计算效率和计算精度。该方法可以为非饱和土渗流的数值模拟提供一定参考。     

基于OpenMP的非连续变形分析并行计算方法

非连续变形分析（DDA）方法严格满足平衡要求和能量守恒,具有完全的运动学及数值可靠性,但对大规模岩土工程问题的数值模拟耗时太长,尤其是线性方程组求解,并行计算可以很好地解决该问题。首先基于DDA方法的基本理论,阐述了适用于DDA方法中的基于块的行压缩法和基于“试验-误差”迭代格式的非零位置记录;其次,引入块雅可比迭代法并行求解DDA方法的线性方程组,并改进了相应的非零存储方法;最后,基于OpenMP实现了DDA线性方程组求解并行计算,并将其应用于地下洞室群的破坏过程分析,以加速比为并行效率的指标评价,结果表明,该并行计算策略可以极大提高DDA的计算效率,而且适合各种规模的问题。     

基于GM(1,1)优化模型的岩石边坡变形预测

岩石边坡系统是一典型灰色系统。其变形发展过程可以用灰色预测模型完成。尽管传统GM（1，1）模型预测有很多成功的实例，但是也存在一些预测偏差过大的情况，必须对其进行优化。逐步迭代法GM（1，1）模型不仅收敛速度快，而且与原始数据序列的凹凸性保持一致。利用自编的计算程序对马步坎边坡观测点G1沉降和开裂进行预测分析，结果表明逐步迭代法GM（1，1）优化模型计算精度较传统GM（1，1）模型和背景构造法GM（1，1）优化模型高，较好地反映了岩石边坡的变形趋势。     

基于GM(1,1)优化模型的岩石边坡变形预报

岩石边坡系统是一典型灰色系统，其变形发展过程可以用灰色预测模型完成。尽管传统GM（1，1）模型预测有很多成功的实例，但是也存在一些预测偏差过大的情况，必须对其进行优化。逐步迭代法GM（1，1）模型不仅收敛速度快，而且与原始数据库序列的凹凸性保护一致。利用自编的计算程序对马步坎边坡预测测点G1沉降和开裂进行预测分析，结果表明逐步迭代法GM（1，1）优化模型计算精度较传统GM（1，1）模型和背景构造法GM（1，1）优化模型高，较好地反映了岩石边坡的变形趋势。     

一种高频面波频散函数的快速算法—改进的Abo-Zena法

作者在本文中首先介绍了改进的Anas Abo-Zena传递矩阵法面波频散函 数的计算问题。这是一种高频时稳定的算法，适用于一般地基勘查和无损检测的面波频散函数计算。其次，还讨论了用Monte Carlo法求解面波频散函数的问题。并用这种算法和Newton迭代法进行了对比，显示了对于Newton迭代法不能计算的速度逆转剖面（即层速度自上而下逐层减递），Monte Carlo法也适用。     

全吸力范围内土-水特征曲线的快速测定方法

传统的宽（全）吸力范围内的土-水特征曲线（SWCC）测量方法动辄耗时数月。为了更好地使非饱和土力学理论指导工程实践,SWCC的快速测定显得尤为关键。为开发宽（全）吸力范围内SWCC的快速测定方法,联合压力板法、并行滤纸法、露点水势仪法进行了两种不同土体的宽吸力范围内SWCC的测定试验。试验结果表明：（1）并行滤纸法与串行滤纸法的试验结果能较好地吻合,故可采用并行滤纸法代替串行滤纸法以缩短测量时间;（2）当测量的试验数据满足SWCC的大致形状时（或满足工程设计需求时）,可适当减少试验数量,以提高测定效率,联合测定方法测量一条完整的SWCC仅需要10～12个点即可完成;（3）联合测定方法可将SWCC测量时间由数月缩短至7～10 d。所提出的联合测定方法实现了宽（全）吸力范围内SWCC的快速、精准测定,有望使得SWCC的测量成为常规的土工试验。     

黑龙江省太平沟幅1：20万金矿矿化集中区地球化学预测

黑龙江省太平沟幅１：２０万金矿矿化集中区地球化学预测孟宪伟（中国科学院地球化学研究所贵州贵阳５５０００２）吴锡生，余先川（长春地质学院吉林长春１３００２６）本文根据１：ＺＯ万地球化学数据结构的特点，提出应在地质单元划分的基础上，运用以迭代法为核心的稳...     

利用稀土分布分形结构因子对碱性花岗岩稀土组成分类——以东准噶尔花岗岩为例

尝试用稀土分布的分形结构因子（FSF）值对我国碱性花岗岩的稀土组成进行了分类。可将我国碱性花岗岩的稀土分布分形结构因子划分为4种类型：（1）弱变异型（0．98－0．95）；（2）中变异型（0．95－0．92）；（3）强变异型（0．92－0．89）；（4）超强变异型（0．89－0．86）。其中前两种类型可很好地对应重稀土富集型和中稀土亏损型花岗岩，而后两种类型可对应轻稀土相对富集型和轻稀土特富型花岗岩。利用FSF值对新疆东准噶尔地区三带富碱花岗岩按稀土组成进行了归类。     

探测管线的倾斜压线法和垂直压线法

通过分析发射线圈与管线间耦合（感应激发方向）的几个特点,提出管线探测的倾压线法和垂直压线法。通过实例说明,倾斜压线法适用于近间距的并行管线,但近于上下并行的不适用;垂直压线法适用于上下并行的管线。当两线间距较小时,用垂直压线法无效,而用倾斜压线法就可探测到明显的目标管导演。在实践中,两种方法灵活运用,即可取得较好的应用效果。     

辽西热河生物群箭蜒科昆虫新类群

描述了箭蜓科（Gomphidae）一新属种一宽翅辽宁澜箭蜓（Liaoninglanthus latus gen．et sp．nov．）及新属一未定种（Liaoninglanthus．sp．），还对义县光箭蜓（Liogomphus yixianensis Ren et Gao，1996）进行了补充描述。化石全部采自辽西北票黄半吉沟义县组下部（晚侏罗世晚期一早白垩世）。     

重力向下延拓的迭代法对比分析研究

向下延拓是位场数据处理和解释的基本方法,但其不适定性导致计算结果稳定性差。近年来,迭代法在提高向下延拓计算稳定性和计算精度方面上的优势备受重视。为比较各种迭代法的优缺点,对积分迭代法、泰勒级数迭代法和补偿法进行了对比分析研究。在模型试算的基础上,将三种方法应用到了江西相山铀矿田重力数据处理之中。研究表明,三种迭代法的迭代通式具有较强的相似性,只是选取的初始滤波算子不同,但初值滤波算子选取的合理性直接决定着迭代法对低频成分的转换能力和高频成分的压制能力。补偿法在向下延拓计算中具有更高的计算精度和更强的计算稳定性,因此选择补偿法具有更好的实际应用价值。     

基于区域分解法的地下水有限元与边界元耦合模型—淄博市王旺庄水源地地下水数值模拟

区域分解法（DDM）是20世纪90年代兴起的一种求解偏微分方程的新方法，方法本身独到的耦合思想和高效的并行计算机理，对于求解复杂的、大型的地下水问题具有相当的优势和广阔的应用前景。本文以淄博市王旺庄水源地地下水流模型为例，应用重叠型区域分解法（DDM）构造了边界单元法（BEM）与有限单元法（FEM）耦合模型，在两种数值方法各自优点的基础上，更形象地再现了实际水文地质原型，有效地消除了人为边界造成的流场失真。     

新疆昆其布拉克地区床板珊瑚、日射珊瑚新资料

描述了由施恩全等于1987年采自新疆且末县昆其布拉克地区上奥陶统孔其布拉克组中上部的晚奥陶世床板珊瑚，日射珊瑚化石，计有8属，13种（其中新种6个）：AgetolitesaltunensisWang，A．KingqiensisWang（Sp．nov．），Cateniporaaqiangensiswang（sp．nov．），Taeniolitescf．zhejiangensisLinetChow，Proheliolitessp．,Priporacf．multispinosaLinetChow，PlasmoporellaqiemoensisWang（sp．nov．），Pl．altunensiWang（sp．nov．），Pl．of．columellaBondarenko，Pl．com-plexawang（sp．nov），Pl．sp，Sibiriolitescf．sibiricusmongolicusBond.etMin.,Diploeporaaltune-sisWnagsp．nov。笔者将其归入Agetolites—Plasmoporella组合。这些化石首次被发现，对于研究阿尔金山地区上奥陶统生物群和地层对比具有十分重要的理论意义和实践意义。     

一种新的并行整体积分和局部微分方程分解（GILD）非线性反演的…

在本文中介绍了基于整体积分和局部微分方程分解（ＧＩＬＤ）的一种新的有效的３－Ｄ电磁模拟和非线性反演算法。这种整体积分和局部微分方程分解并行非线性反算法包括三部分;（１）域被分解成两个子域;子域１（ｓｄ－１）和子域Ⅱ（ｓｄ－Ⅱ）。（２）子域Ⅰ中整体的新的磁积分方程和子域Ⅱ中局部的磁微分方程将共同用来在模拟步骤获得磁场。（３）在子域１中整体的新的磁积分方程和子域Ⅱ中局部的磁微分方法将共同用于反演中不断     

裂变径迹定年法Zeta常数的测定

用第６届国际裂变径迹测定地质年代学术讨论会推荐的三种锆石、两种磷灰石和一种榍石国际年龄标样对Ｚｅｔａ值进行了系统刻度，对ＮＢＳ－ＳＲＭ６１２铀标准玻璃测定的锆石Ｚｅｔａ值为（３３５．６±８．２）ａｃｍ￣２，磷灰石（３４３．３±１４．６）ａｃｍ￣２，榍石为（３４７．７±１５．８）ａｃｍ￣２，三种矿物之间的Ｚｅｔａ值无很大差异，统一的加权平均Ｚｅｔａ值为（３３９．２±６．６）ａｃｍ￣２。     

刚性基础-桩-土非线性共同作用的近似分析方法

提出了一种刚性基础-桩-土共同作用非线性分析更为简便的近似方法。刚性基础基底单元（桩或土单元）本身受力变形的非线性关系采用静载试验得到的Q-s曲线,这种单元自身非线性及单元之间线弹性共同作用,可同时用迭代法计算。算例表明,这种方法方便、实用,能较好地反映桩、土非线性特性对基础受力变形的影响。     

中国新矿物种名录

由Dr．M．Fleischer（美国）和Dr‘J．Mandarlno（加拿大）合著的“GlossaryofMineralSpecies,1995”（矿物种词汇,1995）一书,共计239页（16开本）。该书是一本工具书,收录了截至1994年8月底以前经IMA（国际矿物学协会）新矿物与矿物命名委员会审查通过的独立矿物种共计3519个（变种,同义词除外）。该“矿物种词汇”从1971年开始到1995年为止,共出版了7种版本（1971年,1975年,198O年,1983年,1987年,1991年,1995年）,每一次再版均又增加了一些新的矿物种（也废除了一些矿物种）,这充分说明,当今各国矿物学家对矿物研究…     

三峡库区中侏罗世植物化石及所指示的气候环境

报道了陈家湾组的植物化石８属１１种（包括１新种）即Ｎｅｏｃａｌａｍｉｔｅｓｓｐ．１,Ｎ．ｓｐ２．Ｐｈｌｅｂｏｔｅｎｓｏｆ．ｐｏｌｙｐｏｄｉｏｄｉｅｓＢｒｏｎｇｎｉａｒｔ,Ｔｈａｕｍａｔｏｐｔｅｒｉｓ,ｓｐ．,Ｃｏｎｉｏｐｔｅｒｉｓｓｉｎｐｌｅｘ（Ｌ．ｅｔＨ．）,ＣｌａｄｏｐｈｌｅｂｉｓｓｈａｎｓｉｅｎｓｉｓＳｚｅ,Ｃ．ｘｉｅｔａｎｅｎｓｉｓＭｅｎｇ（ｓｐ．ｎｏｖ．）Ａｎｏｍｏｚａｍｉｔｅｓｇｒａｃ