磁日变可视化编辑软件的开发

通过对航磁生产中磁日变观测的人文干扰产生的原因及特点的分析,利用计算机技术,采取了适当数学方法,设计可视化磁日变人文干扰编辑软件,提高磁日变校正精度,解决生产中实际问题。     

在波数域计算一维重磁异常导数的Matlab语言算法

利用Matlab内建的快速傅氏变换函数可以方便地在波数域计算重磁异常导数。介绍了基于Matlab语言的波数域求导的算法,给出了程序源代码,讨论了一些有助于提高计算精度的编程技巧。通过模型试验和数据分析,发现在计算垂向导数时波数域求导算法的精度比傅氏级数的精度有明显改善;而水平导数的计算,2种方法的精度相当。在某区钾盐勘探中,用该方法处理高精度重力剖面数据,取得了较好效果。     

利用重力场研究东北地区断裂分布及构造分区

为了深入研究东北地区断裂分布及构造分区,利用全国1:2 500 000布格重力异常图,分析了泛东北地区(包括东北3省及内蒙古自治区的东北部分地区)的区域重力场特征.利用布格异常变化形态,并参考地质资料,在研究区内划分出岩石圈断裂14条,壳内断裂11条.根据大兴安岭和依兰-伊通两条明显的重力异常梯级带把研究区划分出兴安、松嫩以及张广才岭等三大重力异常区,在划分断裂和分析重力场特征的基础上进行了构造分区.其研究成果为认识和研究泛东北地区油气与矿产分布提供了重力场方面的依据.     

GPS高程拟合转换正常高的研究

通过对GPS测高与水准测量理论及其异同分析的基础上,阐述了确定似大地水准面的原理与方法,分析了用数学模型法,用少量GPS与水准重合点将GPS大地高直接转换为具有cm量级正常高的实现方法。实验结合黑龙江省虎林地区的地形特点,提出了用线性内插法、平面模型法和二次曲面模型法等来转GPS高程,证明在黑龙江省虎林地区可以通过少量且分布合理的水准点来直接求出该地区具有cm量级的正常高,且精度可以达到四等水准测量的精度要求,满足一般工程的需要。     

基于信号重构的地震道插值

在简要回顾几种插值方法的基础上 ,提出了基于信号重构的地震道插值方法。这种方法无需任何地质和地球物理假设 ,主要用于恢复缺失的地震道。本方法的基础是傅立叶变换理论和最小平方反演技术。文中给出了原理、算法、实现步骤及简要的流程框图。如果对变换域的参数进行规则采样 ,则该方法可以用高效、快速的算法来实现。本插值方法适用于规则和非规则采样的数据。文中给出了理论合成数据和实际数据的试验结果 ,并对结果进行了讨论。同时 ,本方法对于 AVO分析、DMO处理、三维叠前偏移成像等研究也将是有益的     

维纳滤波技术及其在新疆地区重力场划分中的应用

基于维纳滤波原理,应用多层格林等效层模型模拟实测重力场的径向对数功率谱,采用改进的格拉姆-施密特正交化方法计算功率谱密度函数,确保了数值的稳定性,得到了较为满意的结果。通过理论模型实验,证实了利用维纳滤波方法提取深部场时可以保证低频成分不衰减的功能,能够较为精确地提取区域异常,从而实现频率域中位场的分离。运用维纳滤波方法对新疆地区布格重力异常进行了分离处理,获得了基于最佳滤波原理计算的新疆地区区域重力异常,并根据新疆地区布格异常和区域场特征,划分了9条深大断裂和5个一级构造单元。     

龙门山中、南段中-新生代隆升史:来自裂变径迹的证据

运用裂变径迹年代学方法对龙门山中、南段中生代地层内的磷灰石、锆石进行分析,并进行了热历史模拟,结果表明:①龙门山中、南段天全、怀远一带中生代地层最主要的强烈隆升阶段为喜山晚期,与该区的强烈构造运动时期相一致;②怀远镇西南侧及西侧剖面,喜山晚期强烈隆升阶段的时间自西侧的20Ma左右开始至东侧的7～5Ma开始并延续至今,开始强烈隆升的时间自西向东逐步变新;与龙门山造山带逆冲推覆作用在时空上具向东前展式渐进推覆的特点相对应;③部分地区存在163～155Ma、140～77 Ma、55～22Ma的隆升事件.     

拉萨地体内松多榴辉岩的同碰撞折返:来自构造变形和40Ar-39Ar年代学的证据

松多地区的区域构造变形与糜棱质白云母石英片岩和绿片岩的白云母单矿物40Ar-39Ar年代学测试表明拉萨地体内的松多地区于220～240 Ma经历过印支期碰撞造山事件.这次造山事件为晚二叠世松多榴辉岩带代表的古特提斯洋盆消失闭合之后北拉萨地体与南冈瓦那大陆碰撞的结果.该区榴辉岩与退变榴辉岩白云母和角闪石的40Ar-39A...     

利用重力归一化总梯度及相位法研究断裂构造

用波数域的位场转换和快速傅立叶变换理论研究重力归一化总梯度和归一化相位;在波数域向下延拓和导数计算的滤波算子中,分别引入圆滑滤波因子,抑制了向下延拓计算中对噪声干扰的放大作用,增加了计算的稳定性;利用改进的方法,计算并阐述了4种断裂构造GH场等值线和相位曲线的变化特征;别列兹金重力归一化总梯度法中极大值理论在研究断裂构造中具有局限性,对于无限延伸断裂构造,应根据GH场等值线的走向、形状的特点,结合归一化相位曲线的转折点确定断裂空间位置.