点电流源直流场的三维正演模拟及可视化的初步研究

为快速获取直流电法不同三维模型参数下的视电阻率的空间分布及其变化规律,基于异常电位法实现了二极装置下的直流电阻率的三维有限差分正演算法,并基于MATLAB的GUI(Graphical User Interfaces)平台开发实现了人机交互的直流电法三维数值模拟系统,通过不同的三维模型和参数测试验证了算法的正确性和系统的...     

壳—幔侧面波（Pn）速度的横向变化与界面温度横向差异的讨论

本文得到壳—幔界面侧面地震波速度的横向变化与界面温度的横向差异密切相关。侧面波（Pn）速度VPn与地面大地热流值q存在统计关系:     

关于辽宁地区b值方法定量预测的研究

本文利用地震学分析预报方法程式指南和“八五”攻关成果中提出的地震b值定量预测方法,对辽宁及邻区1969年以来发生的8次地震震例作了回顾性检验研究。在此基础上,根据震前b值的变化特征,总结出了可能适用于辽宁地区半年左右尺度的b值预测方法的异常羊别指标及预报规则。     

微分算子展开与耦合势三维有限体积法快速计算随钻超深前视电阻率测井响应非线性Born逼近

本文将非均质各向异性地层中电磁场耦合势Helmhotz方程分解为各向异性背景场方程与散射场方程,并应用算子展开技术推导出耦合势Born级数解,建立一套三维非均质各向异性地层中随钻超深前视电阻率测井响应各阶散射电磁场与非线性Born逼近的有效计算方法.首先,利用传输线法与二维插值技术确定背景电磁场空间分布,得到零阶Born逼近解,然后根据电磁场耦合势的Born级数解推导出各阶散射电磁场方程与递推关系,通过三维有限体积法实现各阶散射场耦合势方程的离散,在此基础上结合MKL PARDISO并行技术,通过递推方式逐步计算各阶散射电磁场,并根据Born级数部分和确定不同阶次的非线性Born逼近.最后利用数值模拟结果对正演算法加以检验,并深入研究分析各阶散射场的变化特征与Born级数的收敛性质.

     

地面瞬变电磁勘探中源外地形的影响

起伏地形为中国矿产资源的发现和勘查造成了极大困难,为了避免复杂的地形效应,在地面瞬变电磁勘探中尽量将发射源与接收点布置在地形平缓处,但源外起伏地形的影响仍然存在。目前,对地形影响的分析大多从测道曲线图或者视电阻率断面图进行定性判断,较少从整个观测时间段内量化评估地形影响的程度。基于有限元法实现了带地形条件下矩形回线源的瞬变电磁三维正演,在源外构造了山脊、山谷地形,将地形与发射源的距离、大地电导率、异常体电导率作为变化参数,给出了整个观测时间段内在不同接收点处地形造成的相对误差分布,从而分析了源外地形对电磁响应的影响规律。结果表明：源外地形对电磁响应的影响主要集中在某个时间段内,随着与发射源距离增大,地形作用减弱,受影响区域也推移到较晚期时间段上;随着大地电导率增大,地形主要影响区域会推移到较晚期时间上,但受影响程度基本一致;相同几何尺寸的山谷地形比山脊地形对电磁响应的影响更大;地形、异常体对电磁响应的影响可以分离,如果异常体产生的二次场比较强,地形效应的影响相对较小,可以忽略地形影响。研究结果为地面瞬变电磁勘探中地形效应的识别与后期处理解释提供了理论依据。     

极低频电磁法超高速振荡Bessel积分数值计算的解决方案——积分外推法

极低频电磁法需要考虑电离层和位移电流的影响,Bessel积分核函数超高速振荡,超过了现有滤波系数的截止频率,收敛慢,甚至不收敛,导致数值计算极为困难,即使采用现有的、最高采样密度的Hankel滤波系数进行数值计算,高频时的计算误差大,不能满足需求.针对以上问题,本文将积分外推算法应用到极低频电磁法的数值计算中.首先通过计算CSAMT在均匀半空间的数值解与解析解验证算法的准确性;然后研究了极低频电磁法三层模型中给定收发距的波场频率响应曲线及给定频率的波场结构特征,积分外推法的低频结果与高采样密度Hankel滤波计算结果一致,高频结果优于Hankel滤波计算结果.本文研究表明,积分外推算法可有效解决超高速振荡Bessel积分数值计算难题.     

爆破加速度测量及标定误差分析

本文首先介绍了在频率域上对爆破加速度测量系统的低频响应误差和高频响应误差的详细分析过程,并对输入、输出波形采用频谱分析及计算幅频曲线,用幅频曲线所包围的不同面积来计算经过测量系统后所引起的误差.其次介绍了采用傅里叶变换法如何在频率域上标定冲击加速度计,给出了确定加速度灵敏度的方法.