介绍两个用 Basic 语言编制的古地磁绘图程序的功能

我们用Basic语言编制了用于古地磁研究的古纬度绘图程序和球面上两个或几个板块绕一欧拉极发生相对转动的绘图程序。其中古纬度绘图程序采用了赤平极射和墨卡托两种投影方式,它能依据古地磁数据绘出板块的古纬度和古方位。和手工绘图相比,它具有直观、快速、准确,适应性强的优点,只要改变计算参数就可以得到不同角度的投影,并可采用不同比例尺,对局部图形加以放大或缩小。使用以上两种程序可帮助人们利用古地磁资料来研究板块的水平运动历史,重建大陆的古地理位置。赤平极射投影绘图程序也可以用于地质学的其它方面(如结晶学等)的绘图工作。古地理图还能直观显示某些沉积矿床生成的古纬度条件,并对某些矿床的现代地理分布作出预测     

Euler Deconvolution with Improved Accuracy and Multiple Different Structural Indices

Euler deconvolution is a semi-automatic interpretation method that is frequently used with magnetic and gravity data. For a given source type, which is specified by its structural index (SI), it provides an estimate of the source location. It is demonstrated here that by computing the solution space of individual data points and selecting common source locations the accuracy of the result can be improved. Furthermore, only a slight modification of the method is necessary to allow solutions for any number of different Sis to be obtained simultaneously. The method is applicable to both evenly and unevenly sampled geophysical data and is demonstrated on gravity and magnetic data. Source code (in Matlab format) is available from www.iamg.org.     

阿穆尔板块相对欧亚板块的欧拉运动

利用中国东北地区1998～2005年的GPS数据,结合俄罗斯、韩国和蒙古的部分GPS观测资料,用加权最小二乘方法以迭代的方式解算了阿穆尔板块相对欧亚板块的欧拉运动矢量。最终的33个测站的位移速率加权均方根残差(WRMS)为0.98mm/a。计算结果表明阿穆尔相对欧亚板块运动的欧拉极位于阿穆尔板块的北部边界外兴安岭东端(54.06±1.13°N,135.87±7.16°E),欧拉旋转量为0.099±0.00735°/Ma。     

最小熵反褶积的微机程序

讨论了最小熵反褶积的基本原理,并给出BASIC微机程序,利用该程序处理了合成地震记录。      

利用微型计算机进行世界古大陆再造成图——方法与实例

进行世界占大陆再造需要3个步骤:（1）确定各古大陆块体的边界;（2）进行必要的陆块拼合处理;（3）恢复各个占大陆块体正当时的方位和位置。第一步是地质学上的判断问题,第二步和第三步是应用微型计算机进行处理,涉及大陆地块在地球表面运动时球面几何学中的欧拉极和欧拉旋转问题。本文介绍在IBM-PC/AT微型电子汁算机上进行占大陆再造的方法,同时讨论了古生代古大陆再造的特殊性。最后给出了震旦纪、泥盆纪5张不同投影实例图作为说明。     

海上地震资料子波提取鬼波压制技术

海洋地震勘探由于受到海面强反射界面的影响,记录数据中存在鬼波,致使数据频谱中产生陷波效应,限制数据的频带宽度,影响数据的分辨率,干扰地层反射的识别.近年来发展的鬼波压制方法,能够较好的压制地震数据中的鬼波,然而现有鬼波压制技术容易产生延续相位,且计算效率较低.为进一步提高鬼波压制效果和工作效率,本文对海上地震资料特性进行分析,研究认为,海上地震资料中的海底反射信号真实记录了海底有效反射、震源鬼波和检波点鬼波的形态,可以利用提取的海底反射子波对数据进行确定性子波反褶积处理,能够有效的压制震源鬼波和检波点鬼波.模型海上实际数据测试结果表明,该方法具有鬼波压制效果好、计算效率高的优势,处理之后的数据频带得到有效拓宽,低频和高频信息均得到加强.     

非调制大容量气枪信号处理中反褶积方法的对比研究

利用云南宾川主动源资料,对比分析了阻尼系数法、水准法、最小二乘法、时间域迭代法以及维纳法去除气枪信号震源响应的效果。结果表明:①在参数选择正确时,以上方法得到的波形差异不大。其中阻尼系数法、水准法参数选取难度最大,最小二乘法虽然能够自动找到最优解,但步长设置过小时,计算效率低,步长设置过大时,易出现假极值点。时间域迭代法及维纳法的参数对结果的影响相对较小。②对于高信噪比叠加信号,时间域迭代法得到的信号信噪比最高,计算时间对比其他方法增加不大,是一种较好的反褶积方法。对于单枪低信噪比信号,当只需要到时信息时,可使用时间域迭代法,但需波形信息时,从波形的可信度及计算效率上考虑,维纳法是更加折中的方法。     

气枪震源波形走时变化分析的反褶积处理

利用云南省宾川主动源基地周边10km范围内7个流动台站2017年2月气枪震源集中激发时段记录的气枪信号,使用干涉法计算了直接使用波形和将波形与参考台反褶积后的格林函数两种处理方法下气枪信号的走时变化。使用了气枪震源附近不同方位角的两个参考台的数据,以分析不同参考台对走时变化结果的影响。结果表明:反褶积能较为有效的压制气枪震源的时间、空间和频谱轻微不重复性导致的走时变化;反褶积能降低一些来源不明的伪走时变化;使用CKT2的数据反褶积的结果略优于CKT0,其可能原因是气枪阵列沿方位角的非均匀辐射。     

基于相似性重构低频数据的金属矿频域全波形反演

由于深部金属矿埋深和自身的复杂性,利用重、磁、电方法和一般的地震方法很难有效地对其进行高精度定位.全波形反演通过最小化模拟数据与观测数据的差异使深部金属矿的高精度探测成为可能,但全波形反演是一个局部优化过程,需要准确的低频数据作为起始,而这在一般的地震数据采集中难以做到.本文先在频域中使用伴随状态震源函数反演方法,通过震源附近的直达波能精确地反演出震源函数的形态.然后利用得到的高精度震源子波结合褶积与反褶积思想及相似性现象重构含有低频成分的自激自收数据.将该数据应用到全波形反演中,有效缓减了反演过程中出现的周波跳跃现象,并提高了模型反演的正确性.Marmousi模型和金属矿模型的数值模拟实验证明了新方法改善了在没有低频数据时的全波形反演结果,并有较好抗噪性.