探地雷达信号的高阶时间域有限差分模拟

探地雷达信号时域有限差分法模拟,多采用二阶精度的中心差分法近似（FDTD（2．2））,虽然计算简单,但数值色散误差较大,影响了模拟精度。在解决较复杂介质分布的探地雷达信号时,不能很好地反映信号的精细变化。而高阶时域有限差分模拟能减少数值色散带来的误差,提高了模拟的精度。采用各向异性完全匹配层（UPML）作为吸收边界条件,实现了高阶计算,有效地吸收边界电磁波反射,而且提高了计算效率。通过模拟的结果分析可知,高阶时域有限差分法能很好地提高模拟精度。     

青海省贵德扎仓寺地热成因机理及开发利用建议

本文通过卫星热遥感影像及航磁、重力异常的宏观性解译热异常空间分布特征,野外实地调查温泉的出露条件和水化学特征,分析研究青海省贵德扎仓寺温泉的热源通道、入渗通道,从而讨论其成因机理和开发利用前景。研究结果表明:扎仓寺温泉泉口温度为93.5℃,与当地气压有关;大气降水和热水沟上游沟谷水为补给源,近东西向张裂为入渗通道,中三叠统砂板岩与印支期花岗闪长岩(γδ₅²)接触变质带为径流通道(热储层),中三叠统砂板岩为盖层,北北西向热光断裂(F₁)为热源通道,属断裂-深循环地热类系统类型;径流通道(热储层)温度可达136.0℃,具备蒸汽发电示范背景条件,而F₁断裂带附近2000 m深度温度可达195.0℃,具备干热岩开发利用前景。     

三维探地雷达数值模拟中UPML边界研究

应用时域有限差分进行探地雷达正演计算时,边界条件的吸收效果是提高雷达波响应模拟精度的关键因素之一。这里推导出了三维各向异性完全匹配空间UPML吸收边界的计算公式,并编程实现。模拟结果表明:UPML吸收边界在吸收效果,反射误差及计算效率等方面,都明显优于常用的Mur、PML吸收边界,具有较好的模拟效果。采用三维UPML边界的时域有限差法,可以提高探地雷达正演模拟的可靠性、准确度。     

地热地球物理勘探新进展

地热作为一种清洁能源具有巨大的开发潜力,将在我国的经济发展中起到巨大的作用.在传统的热水型地热开发的基础上,国际上非常重视热干岩(Enhanced Geothermal)型地热的勘探开发.在这两类地热勘查和开发中,地球物理方法具有非常重要的作用.本文从地热系统的目标体岩石的地球物理性质出发,分析岩石的地球物理性质与温度、压力和含水量等影响因素的关系.例如随着温度的升高,岩石会出现去磁、电阻率降低、密度降低、弹性波速度也现明显降低等现象.进而分析地球物理方法应用到具体的地热勘查地质-地球物理异常模型.结合国际上21世纪以后的新方法技术,分析了重磁、电、地震方法在利用由于岩石温度的升高而出现的特殊地球物理现象,并应用于地热勘探.通过国内外实例介绍了各种地球物理勘探方法在地热勘查中成功应用,为进一步提高我国地热勘查水平,提供一些参考.     

金属矿钻孔雷达探测的数值模拟

为了研究钻孔雷达对金属矿的响应特征,利用时域有限差分法（FDTD）对钻孔雷达探测进行了数值模拟。通过对不同形态的金属矿体,包括等球状体、板状体和一个实际矿体的模拟,得到大量模拟结果。球状目标体在剖面上表现为双曲线形态,通过偏移处理可以精确确定目标体的位置。板状矿体被钻孔穿过时,直达波位置发生错动并变弱,反射同相轴向上下延伸;板状矿体未被钻孔穿过时,反射同相轴呈线状。对于实际矿体,不同的钻孔表现出不同的特征,结合左、中、右3个钻孔的探测结果可以判断出目标体的形态,并且可以判断出此目标矿体具有低阻、高介电常数的特性。     

探地雷达三维高阶时域有限差分法模拟研究

探地雷达数值模拟中,时域有限差分法在时间和空间上一般采用二阶精度的中心差分近似(FDTD(2,2)),其形式简单,但数值色散误差较大,在复杂模型模拟时不能很好地反映模型的精细变化.高阶时域有限差分法能很好地改善数值色散带来的误差,提高模拟精度.本文基于三维高阶时域有限差分法的基本原理实现了探地雷达正演模拟,采用单轴各向异性完全匹配层(UPML)作为吸收边界条件,可以有效地吸收外向传播的电磁波,在大大地提高计算效率的同时,也能很好地改善边界的吸收效果.分析对比正演模拟结果,通过三维高阶时域有限差分正演能获得目标体准确电磁响应信息,并能很好的提高模拟精度.     

三维频散介质中地质雷达信号的FDTD数值模拟

考虑频散介质的电磁波传播,引入随频率变化的电位移矢量D,并对电场强度E和电位移D进行标准化。开发出一种计算介电常数随频率变化满足Debye关系的频散介质中麦克斯韦方程的时间域有限差分（FDTD）解法,并编写了相应的计算原代码。在边界部分引入假象的介电常数和磁导率,实现了边界的无反射吸收,避免了Berenger完全匹配层中场分裂过程,提高了计算效率。计算实例可以看出频散介质中波的传播情况。     

用于航空电磁数据调平处理的导数突变点统计法

航空电磁数据中经常存在条带状误差,为消除这种误差,通常采用窗口滤波。但是,窗口参数及滤波参数要求极为严格,且不能区分调平误差和其类似大小的地质异常,需要操作者多次尝试以寻求合适的处理参数。这里提出了导数突变点统计法,该方法首先在沿伪切割线方向（垂直飞行方向）,对网格化的航空电磁数据求取伪切割线的导数,采用高通滤波方法滤掉幅值比较小的导数突变点;然后统计每条航线上的突变点个数,根据预设的阈值,确定出存在条带误差的航线;再针对这条航线上的数据,采用均值滤波得到调平结果。导数突变点统计法能自动找出存在误差的航线,并只针对误差线进行滤波,能有效消除条带状误差,极大地保留真实地质异常。导数突变点统计方法算法简单,输入参数少,数据处理速度快,为航空电磁数据调平提供了一种新思路。     

公路勘察中滑坡体的地球物理特征与分析——以张榆线公路勘察为例

滑坡是影响路桥安全的重要隐伏地质灾害。地球物理方法作为一种无损高效的勘察方法,可有效查明滑坡体内部的地质构造及滑面等特征。本文以张榆线崇礼隧道出口段滑坡为例,采用高密度电阻率法和剪切波速测试,并结合钻孔资料对滑坡体的地球物理特征进行了分析。该滑坡体视电阻率值和剪切波速值具有明显的垂向差异性,一定程度反映了该滑坡体岩土体性质、密实度、含水率的垂向差异性,综合分析认为视电阻率值和剪切波速值较高的强风化流纹岩顶面为潜在滑面。在此基础之上,基本确定了滑坡体由冲洪积体和坡积物组成,附着于流纹岩体之上,横向长度达300 m,纵向长度260 m,滑坡体厚度达20~30 m。同时,根据高密度电阻率三维剖面分析,认为滑坡体底部滑趾处为剪出口位置。最后,基于滑坡体地球物理特征建立了滑坡体地质模型。     

探地雷达在古河道砂岩体沉积储层刻画中的应用

利用探地雷达勘探技术对鄂尔多斯盆地曲流河三角洲前缘及三角洲平原分支河道露头砂岩体的分布特征进行实测并描述。研究区岩性主要为泥岩、泥岩与粉-细砂岩互层,夹薄层凝灰岩,且多数地段表层覆盖1~5 m厚的黄土。对雷达剖面进行一系列数据处理并定量给出7个主要砂岩体的厚度依次为54 m、83 m、54 m、38 m、104 m、160 m、191 m;高度依次为6 m、7 m、5 m、8 m、6 m、10m、9 m,等参数;砂岩体长厚比值依次为9、11.8、10.8、4.8、17、16、21.2。河道砂岩主体宽与高的比值总体上15,整体呈现砂岩泥岩复合式沉积构造。砂岩体多呈现孤立状顶部平缓底部突起透镜体形态,河道纵向多期沉积明显且有砂岩体横向迁移现象,砂岩体在探地雷达剖面的反射特征与观测结果一致。