Visual Fortran 5.0下的重磁数据处理及图像显示

重磁数据处理在地球物理中占有很重要的地位.利用Fortran 90语言对现有的Fortran 77程序进行了改写,克服了受内存的限制,并且实现了在Visual Fortran环境下的图形图像显示.     

甲玛- 驱龙矿集区地壳结构及成矿指示意义

南拉萨地体位于印度大陆与亚欧板块陆陆碰撞的前缘,发育有大型的铜、钼、金等多金属斑岩型矿床,其地壳深部结构对于研究成矿构造背景和动力学模式有重要的约束作用。基于甲玛- 驱龙矿集区及周边布置的33个宽频带地震台站资料,使用远震P波接收函数方法获得了矿集区的地壳结构和属性信息。地壳厚度和波速比扫描(H- κ扫描)发现矿集区平均地壳厚度为66 km,地壳平均波速比与地壳厚度呈负相关的关系,指示偏长英质物质的相对增加是地壳增厚的主要原因;共转换点叠加(CCP)结果发现矿集区下方没有双层Moho现象,推测矿集区缺失高速的基性下地壳,这也导致了地壳平均属性更偏向于长英质,推测矿集区地壳底部的基性下地壳发生过拆沉作用。新生下地壳是碰撞造山带斑岩铜钼金成矿的重要条件。本研究结果指示新生下地壳的重熔不仅破坏了陆陆碰撞形成的双莫霍面,还使得南拉萨地体地壳进一步分异,下部地壳发生镁铁质矿物堆晶而更偏基性,莫霍面增厚     

青藏高原东北缘岩石圈缩短变形——深地震反射剖面再处理提供的证据

青藏高原是由印度板块和亚洲板块于50～60 Ma碰撞而形成的全球最高最大的高原,已成为多数国内外学者的共识.然而,关于它的岩石圈变形机制却是长期争论的问题.深地震反射剖面是精细揭示岩石圈结构、分辨变形样式的有效技术.重新处理的松潘地块一西秦岭造山带深地震反射剖面揭示出岩石圈变形的细节,以地壳上部的双重逆冲构造、地壳中部...     

深地震探测揭示的华北及东北地区莫霍面深度

从20世纪70年代以来,在我国华北及东北地区进行了大量的深地震探测研究.本文通过对该地区的深地震探测研究的总结和梳理,探讨了该区的莫霍面深度与变化及其地球动力学意义.结果表明:华北地区最深的莫霍面出现在内蒙褶皱带内,最浅的莫霍面出现在渤海湾盆地.东北地区最深的莫霍面出现在大兴安岭地区,最浅的莫霍面出现在依兰-伊通断裂带...     

青藏高原东北缘宽频带地震台阵远震记录波形及背景噪声分析

本文利用远震极性分析和概率密度函数统计法对国家深部探测专项Sino Probe"青藏高原东北缘宽频带地震台阵"40个台站的远震记录波形和背景噪声进行了评估和影响因素分析。结果表明:(1)该台阵对远震有较好的记录和识别能力,且单台定位结果较好;(2)被检测台站背景噪声达到或接近同类地区固定台站的噪声水平标准,观测数据品质较高;(3)根据实际观测资料计算的各台站概率密度函数PDF分析结果建议,在青藏高原东北缘黄土层较厚地区开展宽频带流动观测,应首选基岩台基以保证三分量观测效果,在不得已选择非基岩台基情况下,应深埋地震计以达到最大程度降噪。     

欧拉反褶积与解析信号相结合的位场反演方法

由于解析信号具有不受(二维)或少受磁化方向影响,能够较好反映磁性体边界的特性,因此受到人们的重视.欧拉反褶积法可以确定场源的位置和深度以及形状因子,具有较强的适应性.因此前人提出将二者相结合的方法.针对前人提出的方法中存在受高频干扰严重的问题,本文提出低阶的欧拉反褶积与解析信号相结合的位场反演方法.本方法在反演中只需计...     

广义S变换及其在大地电磁测深数据处理中的应用

广义S变换是一种优于短时窗傅立叶变换和小波变换的时频分析方法,利用广义S变换能够准确定位大地电磁资料中存在的噪声,通过定义时频窗对噪声进行滤除,从而明显提高阻抗视电阻率与相位的估算质量.本文基于S变换和大地电磁测深资料处理的基本原理,研究了基于广义S变换的大地电磁测深资料的处理流程和方法.对理论模拟信号及实测大地电磁场时间序列数据的处理,证实了方法的有效性.     

远震P波接收函数揭示的张家口(怀来)—中蒙边境(巴音温多尔)剖面地壳厚度与泊松比

本文使用时间域迭代反褶积算法,从张家口(怀来)—巴音温多尔一线布设的41个宽频地震台站、1年期间记录的连续三分量数据中提取到高质量的P波接收函数1844个.用H-κ扫描方法获得了测线下方Moho深度与波速比值(VP/VS)进而计算出泊松比,用共转换点(CCP)叠加方法获得了沿测线Moho面起伏图像.结果显示:(1)测线下方Moho深度平均40km,仅各块体边界处出现Moho深度小尺度急剧变化.整体上,Moho面产状相对于索伦缝合带大致对称,在缝合带南侧的温都尔庙带和白乃庙带下方呈南倾趋势,在缝合带北侧的宝力道带、贺根山杂岩带下方呈北倾趋势.(2)华北克拉通北缘泊松比总体较高,兴蒙造山带整体较低;各次级块体内部泊松比分布相对稳定,块体分界带附近往往存在泊松比值的升降扰动.(3)整条测线地壳厚度和泊松比之间存在弱的负相关关系,表明存在构造作用的影响.(4)整条测线泊松比呈现以索伦缝合带南缘为对称轴的非线性分布.本文所获得的地壳上地幔结构以及泊松比分布特征,支持古亚洲(索伦)洋(南北)双向俯冲,最终沿林西断裂闭合的动力学模式.     

华北克拉通北缘(怀来-苏尼特右旗)地壳结构

2009年,中国地质科学院地质研究所与美国俄克拉荷马大学合作实施了一条长453 km的深地震反射、宽角反射与折射、三分量反射地震联合探测剖面. 剖面南起怀来盆地,向北依次穿过燕山造山带西缘、内蒙地轴、白乃庙弧带、温都尔庙杂岩带,到达索伦缝合带. 其中,宽角反射与折射剖面采用8个0.5~1.5 t炸药震源激发,使用300套Texan单分量数字检波器接收,获得了高质量的地震资料. 通过资料分析和处理,识别出沉积层及结晶基底的折射波（Pg）、来自上地壳底界面的反射波（Pcp）,中地壳底界面的反射波（Plp）,莫霍界面的反射波（Pmp）及上地幔顶部的折射波（Pn）等5个震相. 分别采用Hole有限差分层析成像和Rayinvr算法对华北克拉通北缘及中亚造山带南部进行了上地壳P波速度结构成像和全地壳二维射线追踪反演成像. 结果显示：（1）中亚造山带地壳厚度~40 km,变化平缓,低于全球平均造山带地壳平均厚度,可能为造山后区域伸展的结果. 阴山-燕山带附近莫霍明显加深,推测其为燕山期造山过程形成的山根,但该山根很可能在后期被改造. （2）测线中部地壳上部速度较高,对应地表大面积花岗岩出露,而下地壳速度较低,速度梯度低,呈通道状,推测其可能曾为古亚洲洋向南俯冲消亡的主动陆缘,并在碰撞后演变为伸展环境下岩浆侵入的通道. （3）华北克拉通北缘与中亚造山带显示出不同速度变化特征,前者变化相对缓而后者则变化剧烈,二者的分界出现在赤峰-白云鄂博断裂附近.     

深反射地震资料的偏移处理

本文在对偏移理论和相关算法进行概括性的梳理和总结之外,主要针对深反射地震资料叠前偏移方法进行介绍,对初始数据的规则化、偏移孔径、去假频因子以及速度模型等关键参数进行详细的测试和深入地探讨.适当的规则化处理能提高深反射地震资料中反射信号的连续性和剖面整体能量的均一性;合理的偏移孔径能让波场能量很好地收敛、陡倾角反射层能够准确归位以及岩石圈中、深部的精细结构能较好地成像,同时又不至于引入严重的偏移噪音;通过选择适当的去假频因子来减少剖面上的高频干扰,使得整个剖面的高、低频能量分布适中;精确的速度模型是影响偏移结果的关键性因素,通过对速度模型的不断迭代优化,能显著地提高深反射地震剖面的中、深部偏移成像的精度和准确性,为后续深反射地震资料的综合解释打下良好基础.