基于Monte Carlo方法的由GPS观测计算地应变率的误差分析

在利用GPS观测计算地应变率中,应变率的误差分布在地球动力学研究中非常重要. 为解决这一问题,文中提出基于Monte Carlo技术,利用GPS观测的速度计算地应变率进行误差分析的方法;以青藏高原地区为例,利用带有误差的GPS观测数据进行了大量、独立的应变率重复计算,然后运用统计方法给出了应变率分量及其误差的分布. 结果显示青藏高原的喜马拉雅及中东部地区地应变率的东西向、南北向及剪应变率分量的误差相对较大,但整个青藏高原及附近地区应变率结果稳定可靠;研究结果进一步表明,由GPS观测计算应变率的误差主要决定于GPS观测速度的误差, 在变形小的区域（如中国东部地区）,目前GPS的观测精度还不能完全满足地应变率计算的要求.     

全新世初期气候的数值模拟研究

本文尝试用一个完全耦合的大气-海洋-植被模式(Atmosphere-Ocean-Vegetation General Circulation Model, AOVGCM)模拟全新世初期的地球气候.模拟结果表明,耦合模式成功地再现了全新世初期复杂气候条件下的基本特征,为研究这一时期的气候状况提供了重要的模型支持.就全球平均而言,11ka BP冬季地表气温比现代约低1.6 K,夏季比现代低约0.3 K;大气温度从低层到高层有不同的表现,高层反映了太阳辐射的重要作用,而低层气候对下垫面（如冰川、植被和海洋等）的影响比较敏感.从区域分布来看,11 ka BP冬季大部地区比现在偏干,但热带太平洋和南半球少部地区降水偏多;夏季大部地区比现在偏湿,亚洲和非洲季风偏强,主要的季风区降水偏多.     

大气掩星反演误差特性初步分析

GPS大气掩星探测技术可以获得全球大气折射率、气压、密度、温度和水汽压等气象参数,该技术基本原理是基于几何光学近似的Abel积分反演.地球扁率、电离层传播时间延迟、大气大尺度水平梯度、多路径传播现象等因素在某些高度范围影响大气反演的精度.本文采用模拟的方法,分析其中地球扁率及电离层对反演结果的影响,并讨论局部圆弧修正及电离层修正的效果.利用CHAMP掩星实测轨道数据和有关电离层和大气经验模式、采用三维射线追踪方法模拟计算几种情形下的GPS掩星观测附加相位数据,对模拟数据进行反演,将反演气象参量剖面与模拟时给定模式剖面进行比较,得到了0～60 km高度范围内的反演误差.误差统计分析结果表明,局部圆弧中心的修正以及电离层修正,对于高精度的GPS掩星反演是非常重要的;电离层修正残差仍是制约30～60 km高度范围内反演精度的重要因素,进一步完善和优化大气掩星反演需要发展新的电离层修正算法.     

井地电阻率法歧离率确定高阻油气藏边界

根据油气藏的特征,抽象出地球物理模型,利用有限元软件FEPG对三维油气藏模型的稳恒电流场进行了数值模拟.根据井地电阻率法计算结果,定义了一种新的参数歧离率,分析了不同深度油气藏和油气藏偏离钻井位置时歧离率的变化规律,并利用该参数及其梯度确定油气藏边界范围.理论分析表明,参数歧离率可很好地突出油气藏异常并确定其横向边界范围.     

灌注桩钢筋笼的充电电场特征研究

充电法是用于良导电矿体勘探的传统方法.利用良导椭球体的充电电场特征,结合模型桩试验,分析了灌注桩中钢筋笼周围的充电电场特征.结果表明在钢筋笼底边缘附近,电位曲线急剧下降,出现明显拐点,梯度曲线相应出现极小值.研究成果为检测灌注桩钢筋笼长度提供了理论依据.     

瞬变电磁法拟地震成像研究进展

瞬变电磁法是工程地球物理勘察中应用较多的一种勘探方法之一,瞬变电磁法拟地震成像方法研究是当前电磁探测理论与应用研究的热门.文中详细叙述了瞬变电磁法拟地震成像的两种不同思路,一种是建立时间—频率转换经验公式,把TEM数据转换成平面波数据,借鉴大地电磁法拟地震解释思路进行成像,另一种是通过积分变换,把扩散场变成虚拟波场,把电磁衰减信号变成拟地震子波,然后按照地震解释法进行解释,并分析了各自存在的的问题.指出下一步研究的方向.     

Born近似快速三维反演井地电法数据

本篇研究了井中电偶极激发地面接收的井地电法的快速反演成像问题.我们采用了Born近似方法和重加权正则化共轭梯度法(RRCG)算法.数值计算的结果表明Born近似是一种有效的井地电法三维快速反演方法,同时也说明井地电法监测油水前驱和储层边界预测的观测数据可以用该方法进行快速三维反演成像.     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.     

有限元方法在地震孕育过程与前兆机理研究中的应用进展

有限元方法的基本思想是将连续的求解区域离散为一组有限个、 且按一定方式相互联结在一起的单元的组合体。 有限元方法是目前应用最为广泛的数值模拟方法, 由于能够将复杂介质的力学本构关系及边界条件问题转化成为常规问题的计算程序, 所以越来越多的地震学者将其应用于地震成因与前兆机理研究中。 该文着重介绍了有限元数值分析方法及其模拟在地震孕育过程与前兆机理研究中的应用进展, 并提出开展小尺度规模如岩石中矿物颗粒之间的力学行为的有限元模拟研究等几点建议。     

地球物理学中的电磁场正演与反演

本文在近年来众多的地球物理研究者的研究基础上,总结了当前地电磁模型正反演已有成果,定量分析了各种主要正反演的性能测试,指出不同正反演法的优点、缺点以及应用范围局限,提出了各种方法的发展趋势以及当前计算地球物理领域的核心内容,指出了计算地球物理领域的数值模拟发展方向.