天水地震区深部地温场的有限元计算

用有限元方法计算了天水地震区上地壳地温场。结果表明，影响该区浅部（层）地温场的主要因素是基底构造和断裂构造。祁连褶皱系为高热异常区，秦岭褶皱系为低热异常区，西秦岭北缘断裂带为地温梯度带。认为该带历史地震的发生与热应力差有关     

新疆北山地区区域重、磁场特征及其地质意义

本文在探讨新疆北山地区区域重、磁场特征的基础上,系统地叙述了北山裂谷内发育的白地洼、磁海重力高异常特征和基性、超基性岩体群重、磁场特征。指出坡北、笔架山岩带处于EW、NEE、NNW三组断裂所控制的三角形重力高的梯度带和负磁异常带上。提出通过这一特点,可以在北山地区寻找新的基性、超基性岩体群和拔矿靶区圈定。     

水平气压梯度力差分格式误差的定量分析

通过分析发现地形坐标系中气压梯度力差分格式的计算误差应分为二类，其中，第二类误差在水平坐标面倾斜时出现。这种误差可理解为将水平坐标面上的气压(或位势高度)插回到等高面(或等压面)上的插值误差。 构造气压梯度力差分格式的目的正是为了减小这种第二类误差，而不是其他误差。 误差分析表明静力扣除法和Corby 格式的第二类误差都很小，比一般方法的误差小约一个量级。分析还表明第二类误差产生的主要原因是由于气压随高度非线性变化。依据这一分析思路，本文改进了#As(z)#a坐标系中的静力扣除法。 改进后的方法比改进前的误差小约一个量级。     

航磁ΔT平面归一化总梯度方法及其在内蒙陆西凹陷的应用

本文首次提出将剖面归一化总梯度方法推广到平面加以应用,得到平面归一化总梯度方法,并根据剖面归一化总梯度的计算公式相应地推导了平面归一化总梯度的计算公式,将其应用到内蒙陆西凹陷探测油气田,取得了一定的地质效果。     

应用地形信息精化山区大地水准面

在介绍由高程异常确定大地水准面有关理论的同时,提出了用地形质量计算重力异常垂直梯度的方法,结合高程异常资料,可以使山区大地水准面得到精化。将本方法应用于珠穆朗玛峰,获得其大地水准面高为－ ３０．３６ ｍ 。     

A field study of the coupled effects of aquifer stratification,fluid density,and groundwater fluctuations on dispersivity assessments

A number of experimental studies have tackled the issue of solute transport parameter assessments either in the laboratory or in the field. But yet, the behavior of a plume in the field under density driven forces, is not well known due to possible development of instabilities. Some field tracer tests on the fate of plumes denser than native groundwater such as those encountered under waste disposal facilities, have pointed out the processes of sinking and splitting at the early stage of migration. The process of dispersion was widely investigated, but the range of dispersivity values obtained from either experimental tests, or numerical and theoretical calculations is still very large, even for the same type of aquifers. These discrepancies were considered to be essentially caused by soil heterogeneities and scale effects. In the meantime, studies on the influence of sinking and fingering have remained more scarce. The objective of the work is to analyze how transport parameters such as dispersivities can be affected by unstable conditions, which lead to plume sinking and fingering. A series of tracer tests were carried out to study under natural conditions, the transport of a dense chloride solution injected in a shallow two-layered aquifer. Two types of experiments were performed: in the first type, source injection was such that the plume could travel downward from one layer to the other of higher pore velocity, and in the second one, the migration took place only in the faster layer. The results suggest some new insights in the processes occurring at the early stages of a dense plume migration moving in a stratified aquifer under groundwater fluctuations, which can be summarized through the following points: (i) Above a stability criterion threshold, a fingering process and a multi modal plume transport take place, but local dispersivities can be cautiously derived, using breakthrough curves matching. (ii) When water table is subject to some cycling or rising, the plume can be significantly distorted in the transverse direction, leading to unusual values of the ratio between longitudinal and transverse dispersivities. (iii) Under stable conditions, for example in the case of straightforward injection in the faster aquifer layer, longitudinal dispersivity is greater than the transverse component as usually encountered, and the obtained transport parameters are closed to macro dispersivity values, which reach their asymptotic limit at very short distances. (iv) The classical scale effect about the varying dispersivity at short distances could be a process mainly due to the distance required for a plume stabilization.     

井地直流电法三维数值模拟中若干问题研究

讨论了地下垂直线源分段计算和场叠加的方法,并实现了在套管上供直流电的三维数值模拟。讨论了大型容量矩阵的压缩存储方式,采用数组和结构体相结合的方法实现容量矩阵的一维链表式压缩存储。在求解超大型稀疏线形方程组时引入不完全Cholesky分解稳定化的双共轭梯度算法(ICBG),通过与均匀半空间垂直线源解析解的对比,证明了该算法是准确可靠的。     

引力向量及梯度张量在坐标系之间的变换

引力向量和引力梯度张量在局部和全局坐标系间的转换过程是对卫星重力梯度测量进行模拟研究的主要手段和重要步骤.提出了利用计算机符号运算辅助推导该转换过程的新思想.首先通过对该过程进行分析,指出其运算繁琐复杂;然后通过对理论推导过程的分解,确定了引力向量和引力梯度张量在局部和全局坐标系之间转换的各关键步骤;并借助于计算机符号运算实现了这一转换过程的计算机推演,成功地得到了完整的符号公式结果.所编写的符号推演程序在多台计算机的多种平台上获得了一致的结果.实验结果表明,计算机符号运算对于公式体系繁杂的卫星重力学理论研究来说具有极大的增效作用.     

直流激电测深二维反演的若干问题研究

目前,激电数据的二维反演主要针对高密度激电数据,而对于常规的激电测深则很少涉及.据此,针对常规激电测深二维反演所涉及到的一些技术难点进行探讨.在正演过程中,根据实测数据实现有限元网格自动剖分,并采用双网格系统,用相邻节点电导率双线性连续进行有限元模拟,并对模拟的误差进行校正.在反演过程中,利用互换原理和Broyden的方法计算偏导数矩阵;采用与奇异值分解算法精度相当的共轭梯度算法求解偏导数矩阵方程.最后通过对模拟数据和实测数据进行反演试算,验证本计算方法是可行的.     

地面-坑道激电测深方法在山东玲珑断裂带深部找矿中的应用

在山东玲珑断裂带深部找矿中,考虑断裂带处在地电干扰和地形变化复杂地段,且深部地质体产状变化较大,仅依据半空间地面物探方法很难摸清地质体的展布形态,盲目布设钻孔会造成找矿成本的提高。为精确反演深部矿体的边界范围、走向及倾向,提高钻孔见矿率,降低找矿成本,在山东玲珑断裂带地面、坑道采用立体物探技术,试验通过布设多条方位不同的地-坑、坑-坑和地-地纵横激电中梯剖面和双向三极测深剖面,进行了二维电法剖面反演。在地面-坑道激电测深方法中,地面-坑道供电,坑道测量探测深度与坑道深度呈正相关,反映了深部巷道底板以下电阻率和极化率的分布状况,可以实现深地探测目的。由于地面-坑道供电,坑道测量装置处在半空间和全空间两个部分,需要采用"镜像法"和地下倾斜偶极源电流场特征对其装置系数k值重新推导,以反映坑道空间岩矿体真实电阻率的变化。坑道中梯剖面异常与单向三极剖面异常、双向三极电测深反演结果对应较好。经钻探验证,矿体分布与物探推断高极化体基本一致,为快速准确布设深部钻孔提供资料。