解地下水渗流问题的有限差分法(二)

五、隐式差分方程的叠代解法 用隐式差分格式解地下水非稳定流问题时,每个时段都要解代数方程组。当含水层导水系数不随水位变化,或变化小到可以忽略不计时,方程是线性的。线性代数方程组的解法虽很多,但都可以归结为两类:直接解法及叠代解法。传统的直接解法工作量随结点数的增多而上升得很快,所以只适用于规模较小的问题。叠代方法也有很多种,对规则网格来说现在有     

旅行时分解法初至折射波静校正及其应用

在煤田地震勘探中,当第四系地层厚度较薄时,从野外观测到资料处理都较难得到来自第四系地层底界的反射波,因此无法解释第四系地层厚度问题,而这一问题又常常是要解决的地质任务之一。并且第四系地层厚度的变化也直接影响了深层反射波的叠加效果和构造形态。就此,文章介绍一种初至折射波静校正方法──旅行时分解法,及其在资料处理与解释中的应用。     

△T球体异常的三元解释法

在△Ｔ异常解释课题中，一个新的角度参量─—视磁化特征角的拟定，最大限度地简化了球体磁场表达式。从而才有可能通过特殊解法，给出了球体磁场的极值与拐点公式。作者据此建立了四种不受正常场制约、仅需异常的局部持征曲线即可全面解释三度异常的解释方法。屏障分解、叠加判别、三元迭代的设计，为解释三度体叠加异常给出了简单而易行的途径。     

线性方程组的正则化解法及其应用

在重力资料反演解释中,正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子,求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中,要求对所求的解答预先有一个初步的推断解释。为此,要求对所解释的问题先进行一番研究,并对综合利用各种方法和手段所获得的先验资料进行对比,提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面,与正则化法和其它线性方程解法的结果作对比,说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。     

二氧化钛—水体系氧同位素分馏系数的实验研究

对ＴｉＯ２的两种常见同质多像变体金红石和锐钛矿与水之间的氧同位素分馏系数的实验研究进行了系统的总结和评述。水热晶化法和水热氧化法以及低温水解法获得的实验金红石－水体系分馏系数与增量方法理论计算结果相一致。低温水解法获得的够太矿－水体系分馏可能是一种不平衡分馏,其分馏系数的大小决定于锐钛矿的形成机制。同时,低温水解实验还揭示,在低温矿物形成和同质多像转变过程中可能存在氧同位素继承性,这对解释低温环境     

全空间瞬变电磁法资料解释方法中的平移算法

全空间视电阻率公式及全区求解方法是地下瞬变电磁理论的重要组成部分。为探讨全空间全区视电阻率求解方法,笔者将半空间全区视电阻率的平移算法引入至全空间并对其进行了研究。同时,分别采用理论数据、正演模拟数据及实测数据对该方法在全空间资料处理与解释中的有效性和精确性进行了分析。结果显示,平移算法与常规的全区视电阻率解法具有相同的理论基础,但前者避免了传统全区计算方法中的多解性;此外,采用由截距参数绘制的时间剖面图和视电阻率断面图进行联合解释产生了较好的应用效果。     

线性方程组的正则化解法及其应用

在重力资料反演解释中，正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子，求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中，要求对所求的解答预先有一个初步的推断解释。为此，要求对所解释的问题先进行一番研究，并对综合利用各种方法和手段所获得的先验资料进行对比，提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面，与正则化法和其它线性方程解法的结果作对比，说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。     

关于求解碾磨机运动问题的相关讨论

针对关于求解碾磨机上一点运动的速度及加速度的问题展开讨论,从运动的不同角度考虑,给出了n种不同的解法。并将此问题拓展到一般情况下进行讨论,得出正确结论。并就一些相关问题进行了分析和解释。     

富里叶变换解释倾斜厚板状体磁异常的应用

迄今,国内外对磁异常进行富里叶变换求其频谱以解释磁测资料方面的文章发表得很多,V.L.S.伯赫美树卡雷姆等人的《用富里叶变换解释岩脉磁异常》一文即是其中之一。该论文在频率域给出了岩脉顶宽、顶部埋深及磁化方向与岩脉之间夹角等参数的准确解法,但尚存在着如下缺陷:     

关于频率电磁测深几个问题的探讨(七)——解析广义逆矩阵反演法

从目标函数的线性化出发,采用奇异值分解法求广义逆矩阵并适当加阻尼,用迭代法直接解超定方程实现改进广义逆矩阵反演。进一步从理论与实际结合上阐明改进广义逆矩阵反演法的灵活、稳定,它涵盖了多种反演方法的优点,又具有自身特点,在实际应中适应性强,并能提供一些辅助信息,可更好地评价解释结果。      

赞坎磁铁矿两种磁异常解释的应用

主要探讨在用磁法寻找磁铁矿的时候的应用,对异常的地质解释和解决其它地质问题时的解释一样,不可能提出一个对任何条件都适用的标准化方法。下面就新疆赞坎磁铁矿,讨论其解释方法和途径。     

地幔对流的数值模拟方法

地幔对流的数值模拟研究为理解地壳运动、地幔内部的性质和地球的动力学机制提供了手段。笔者回顾了地幔对流数值模拟的发展,总结了各种常用的描述地幔对流的数学模型,边界条件以及数值解法。分析了级数展开方法,有限差分方法,有限元方法和多重网格方法的性质和特点,针对某些方法编写了二维地幔对流的计算机软件,使用了实际的地质参数,无穷大的Ｐｒａｎｄｔｌ数,在１０４到１０７之间的Ｒａｙｌｅｉｇｈ数,并给出了数值结果和地质解释。     

植物化石角质层分析实验技术总结及探索

植物化石角质层的微细构造特征是植物化石重要的鉴定和分类依据之一,同时也是获取地质历史时期气候信息的重要途径之一。对保存有角质层的植物化石,如何更好的获取角质层特征,选择合适的实验技术尤为重要。本文采用三种不同的实验方法,即王水+次氯酸钠浸解法,次氯酸钠浸解法以及硝酸浸解法,对地质时期植物化石角质层的获取方法进行了探索性实验研究。实验结果表明,对于砂质泥岩中的植物化石适合用王水+次氯酸钠浸解法进行处理;含煤层或者聚煤期泥岩中的植物化石适合用次氯酸钠浸解法进行处理;灰白色泥岩或者含有硅藻土泥岩中的植物化石适合用硝酸浸解法进行处理。本文采用的三种角质层浸解法在角质层研究方面可根据化石材料的具体情况进行选择,也可同时进行三种方法的尝试进而获得最佳实验效果。该研究总结了不同地质时期植物化石角质层对三种实验方法的反应,为研究植物化石的角质层特征提供更加有效的实验途径。     

边坡稳定分析的非线性有限元混合解法

在对有限元强度折减法和有限元迭代解法的优缺点分析的基础上,提出了边坡稳定分析的非线性有限元混合解法,该法将有限元强度折减法和有限元迭代解法联合运用于边坡稳定分析,充分利用了两种方法的优点。即由有限元强度折减法搜索边坡可能滑动面,将可能滑动面在网格中画出,由迭代解法逐步迭代求得其滑动面的安全系数。并以算例和工程实例说明了此种方法的正确性和合理性。     

利用重磁场对断裂构造进行解释推断

断裂构造是构造运动的一种重要表现形式,在地壳发展演化的各个不同阶段都有发生。断裂构造解释推断工作是重磁资料研究的主要工作内容。利用重磁场对断裂构造进行解释推断,特别是对区域性隐伏的构造的解释效果要优于其它方法。解释推断方法是以重磁场的相对变化为依据,以数据处理为手段。通常数据处理以场的转换、异常分离、定性、定量解释为主要内容。     

改进的周文德法在解释抽水资料上的应用

该文对非稳定流抽水试验资料整理采用的周文德法提出了改进意见,即用迭代法求解特殊函数方程计算T、a、s,代替周的对数纸绘解法,计算简便,求解精度高。作者用实例计算,并以结果与雅柯布法和配线法相比较,提出了雅法的运用条件,得出含水层的参数是随s的变化而不同。而配线法则由于人为因素造成误差,进而认为改进的周法是解释非稳定流抽水资料的一种精确可靠的方法。     

剩余降深法解线性非齐次系统响应矩阵

本文就线性非齐次地下水系统提出了响应矩阵的直接解法——剩余降深方法,并且给出了较严格的理论证明。它把响应矩阵的求解分为三个步骤:一求附加降深场,二求综合降深场,再求剩余降深场,即得到响应矩阵。该方法不仅为以往的间接解法找到了一种对比解法,而且有可能据此向非线性系统拓展。颇具理论和实际意义。     

二维大地电磁资料的人机联作解释

本文介绍了在PC微机上实现的大地电磁资料的二维人机联作解释软件(简称MTI2D)。该软件由建模、二维正演和绘图三部分组成。建模、绘图部分采用TURBO BASIC编程,二维正演部分采用有限元解法,用FORTRAN77写成。三者相互独立,通过SHELL功能将三者连在主莱单下。此外,还说明了在建模过程中涉及的几个问题。 经理论验证和实测资料解释表明,该软件使用方便,操作简单,运行稳定可靠,具有良好的人机对话功能。     

一维大地电磁资料的直接反演法——修改的COEN近似法

本文对一维大地电磁资料的COEN近似法作了一些重要的修改.其中包括将输入数据KC-1/KC+1改为大地电磁的反射函数KC-1/2;采用广泛应用的奇异值分解法(SVD)解线性方程组,并对模型S§加以合理的物理解释.从而取得较满意的结果.     

应用折射层析方法解决塔里木盆地浅部地震速度的分布

首次提出了用于解决塔里木盆地地震资料数据处理中浅部地震速度分析的折射层析成象方法。该方法可以通过求解任意变速介质中的二维程函方程、射线追踪和ＳＩＲＴ迭代解法实现。通过对塔里木盆地地震资料的层析反演取得了较好的结果,揭示了浅部地震速度具有较大的横向变化。这一变化对建立地下速度模型及反射构造成象有较大影响,钻井资料证实了这一结果。这项研究的成功对研究塔里木盆地地震波速度的变化规律及进一步提高该区整体地震资料处理与解释的质量具有重要的实际意义。