浅论真振幅权函数的数值计算

真振幅权函数是Kirchhoff 偏移的核函数,数值计算方案与其表达形式有关。为探讨Kirchhoff 型真振幅加权函数的最优计算方案,对出现在文献中典型的3 种表达形式的权函数的数值实现流程进行了对比分析。结果表明: 基于格林函数振幅的计算方案需要进行动力学射线追踪,最为耗时; 基于旅行时和基于坐标与慢度的计算方案仅需运动学射线追踪; 基于旅行时的计算方案把权函数完全表述成旅行时导数的形式,方便易行,但算法的稳定性对旅行时的依赖度很高; 而基于坐标与慢度的计算方案对复杂模型有很好地适应性。此外,文中还特别讨论了2. 5D 介质下基于旅行时的权函数数值计算方案并进行了实现。     

我国全新世泥炭α纤维素稳定碳同位素记录的对比与分析

泥炭是古气候研究的宝贵材料,α纤维素稳定碳同位素（δ¹³C_cell）是泥炭古气候研究的常用代用指标。我国已有数个地点报道了全新世泥炭全样或单种植物δ¹³C_cell记录,对这些记录的古气候解译存在明显冲突。结合相应的孢粉和植物残体证据,对这些记录进行综合对比,结果发现：高山草地植被带和内陆干旱区荒漠草原植被带当中的泥炭地,全新世期间的原地植物种属构成相对稳定且未发生显著变化,因此其δ¹³C_cell记录具有较为明显的长期变化趋势,其古气候解译也更为可靠;森林植被带和林草交互带当中的泥炭地,全新世期间的原地植物种属构成发生过显著变化（如木本植物入侵）,这可能是造成其全新世δ¹³C_cell记录长期趋势不明显或其古气候指示意义复杂难解的主要原因。这些结果显示了泥炭地的原地植物种属构成,对相应δ¹³C_cell记录的古气候解译的重要影响。同时也表明,无论是泥炭δ¹³C_cell研究地点的选择,还是泥炭δ¹³C_cell记录的古气候解译,都需重视其来源植物种属构成和变化。对于全新世泥炭δ¹³C_cell古气候研究,建议开展针对性的现代过程研究,尤其是优势泥炭植物种属的δ¹³C_cell在时间序列上的现代过程研究,以提供可靠的古气候解译基础。

     

均匀介质中F-K反偏移法正演模拟

F-K反偏移法正演模拟是在给定地质模型后通过构建人工偏移场,再对人工偏移场进行F-K反偏移来实现。人工偏移场的构建是在地质模型上直接放上经正确拉伸和幅度校正后的地震脉冲子波。F-K反偏移中插值映射采用两点sinc插值,其能避免假同相轴的出现且具有较快的速度。奇异点处理使用泰勒近似展开。最后对凹陷模型进行了F-K反偏移法正演模拟实现,得到了预期的效果。     

Kirchhoff型真振幅反偏移加权函数分析

Kirchhoff型真振幅反偏移的基本目的是通过等时线叠加得到偏移的输入场,是将深度偏移的结果再反映射到时间域去,而真振幅加权函数是Kirchhoff型真振幅反偏移的核心内容,它没有唯一的数学形式,目前发表在文献当中的加权函数公式是在物理直观下的最优表述.本文在前人研究的基础上推导了真振幅反偏移加权函数公式,对影响加权...     

新疆137Cs定年时标解析

¹³⁷Cs在年代测定和土壤流失等领域得到了广泛的应用。在认定¹³⁷Cs时标前,需要消除粘土和有机质（TOC）对¹³⁷Cs 的富集作用。通过对比新疆多个湖芯钻孔¹³⁷Cs 序列,发现它们异于北半球¹³⁷Cs 沉降模式,即出现了1954、1963、1975 和1986 这四种时标且1986 年的¹³⁷Cs 比活度可能异常的高。通过分析艾里克湖钻孔的¹³⁷Cs序列,并结合新疆其他地区多个钻孔的¹³⁷Cs序列,与北半球¹³⁷Cs沉降特征以及靠近北疆的Belukha冰芯²⁴⁰Pu²³⁹Pu沉降记录进行对比,分析了新疆¹³⁷Cs时标的可靠性以及1986年¹³⁷Cs比活度异常高的原因。在新疆¹³⁷Cs多个时标中,仅1963年的时标是高信度的（即北半球¹³⁷Cs最大沉积年）。1954年¹³⁷Cs蓄积峰代表其全球初始显著沉降,距今已有2个¹³⁷Cs半衰期,当年沉积通量低且可能发生¹³⁷Cs垂直迁移,仅能作为次要时标。1975年的时标反映了1960s- 1970s的中国核试验,但距离罗布泊～300 km的博斯腾湖对中国核试验有截然相反的记录,表明中国核试验不是稳定信号,其核沉降范围以及对湖芯的影响作用还需进一步评估,不能作为可信的时标。1986年蓄积峰对应切尔诺贝利核事故,因核尘埃未大量进入平流层且距离新疆较远（～4 000 km）,但同时考虑到新疆位于核事故的西风环流下风区,所以该事故对新疆有影响但影响不大。1986年的¹³⁷Cs时标可作为次要时标。新疆湖芯1986年的异常高蓄积峰可能是由新疆最近几十年的人类活动（如过牧、耕种等）加上1987/1988年的高降水导致的地表侵蚀增强引起的。增强的地表侵蚀导致含有高¹³⁷Cs含量的表层土壤颗粒被冲刷进入湖泊,叠加在切尔诺贝利核事故¹³⁷Cs沉降之上,进而导致了湖芯沉积物记录到显著的1986年¹³⁷Cs蓄积峰。     

F-K反偏移中的插值映射

均匀介质中的F-K反偏移是实现一般介质下F-K反偏移的基础,而均匀介质中的F-K反偏移的实质是一个插值映射过程。根据插值对象为复函数的特点,采用了实部虚部分别插值、去振荡因子法、模和相位分别插值、直接频率域插值4种插值方法对点脉冲模型进行了试算。结果表明:截断的直接频率域插值在适应性和精度综合衡量下为最佳;去振荡因子法精度虽高,但可应用的模型太具理想,不适合一般模型中的应用;实部虚部分别插值与模和相位分别插值两者精度较低。     

新疆137Cs定年时标解析北大核心CSCD

^(137)Cs在年代测定和土壤流失等领域得到了广泛的应用。在认定^(137)Cs时标前,需要消除粘土和有机质(TOC)对^(137)Cs的富集作用。通过对比新疆多个湖芯钻孔^(137)Cs序列,发现它们异于北半球^(137)Cs沉降模式,即出现了1954、1963、1975和1986这四种时标且1986年的^(137)Cs比活度可能异常的高。通过分析艾里克湖钻孔的^(137)Cs序列,并结合新疆其他地区多个钻孔的^(137)Cs序列,与北半球^(137)Cs沉降特征以及靠近北疆的Belukha冰芯^(240)Pu/^(239)Pu沉降记录进行对比,分析了新疆^(137)Cs时标的可靠性以及1986年^(137)Cs比活度异常高的原因。在新疆^(137)Cs多个时标中,仅1963年的时标是高信度的(即北半球^(137)Cs最大沉积年)。1954年^(137)Cs蓄积峰代表其全球初始显著沉降,距今已有2个^(137)Cs半衰期,当年沉积通量低且可能发生^(137)Cs垂直迁移,仅能作为次要时标。1975年的时标反映了1960s-1970s的中国核试验,但距离罗布泊~300 km的博斯腾湖对中国核试验有截然相反的记录,表明中国核试验不是稳定信号,其核沉降范围以及对湖芯的影响作用还需进一步评估,不能作为可信的时标。1986年蓄积峰对应切尔诺贝利核事故,因核尘埃未大量进入平流层且距离新疆较远(~4 000 km),但同时考虑到新疆位于核事故的西风环流下风区,所以该事故对新疆有影响但影响不大。1986年的^(137)Cs时标可作为次要时标。新疆湖芯1986年的异常高蓄积峰可能是由新疆最近几十年的人类活动(如过牧、耕种等)加上1987/1988年的高降水导致的地表侵蚀增强引起的。增强的地表侵蚀导致含有高^(137)Cs含量的表层土壤颗粒被冲刷进入湖泊,叠加在切尔诺贝利核事故^(137)Cs沉降之上,进而导致了湖芯沉积物记录到显著的1986年^(137)Cs蓄积峰。     

2维和2.5维起伏地表直流电法有限差分数值模拟

起伏地表直流电场数值模拟现多采用有限元法,主要是因为其有灵活的处理曲边界的能力,然而有限元法比有限差分法要复杂,如果让有限差分法也同样具有较好的处理曲边界的能力,那数值模拟将变得更为简单.本文通过在非正则内点处采用不等距差分,在起伏地表点处直接实现边界条件,克服转移法人为改变地表形状的弊端,使得基于笛卡尔网格的有限差分...     

二维起伏地表条件下坐标变换法直流电场数值模拟

针对有限差分方法处理不规则边界难的问题,引入曲化平的思想,实现了一种新的基于坐标变换法的起伏地表条件下直流电场数值模拟方法.通过坐标变换,将物理空间的直流电场满足的基本方程和边界条件变换到计算空间;在计算空间内,利用有限差分方法对变换后的方程和边界条件进行数值分析;最后,将计算结果由计算空间映射回物理空间.精度分析及计算实例表明:本文提出的算法精度高,计算误差主要集中在源附近,且当源离地表较近时平均相对误差在1.0%以下;计算结果满足起伏地表条件下电位及视电阻率的分布规律.     

二维起伏地表直流电场插值法数值模拟

本文提出了一种采用插值法实现边界条件的有限差分数值模拟新方法, 地表函数通过三次样条插值能很好的拟合真实地表形状。其地表边界条件的实现思想是在地表法线方向上寻找一个与地表点位于同一等位线上的点, 该点通过周围点插值来建立离散方程。采用不规则四边形距离加权插值法使插值形式统一化。在水平地表情况下做的数值解与解析解的对比分析表明: 该方法的计算精度适合要求, 对山谷模型的实例计算也得到了视电阻率与地形成镜像对称的结果。