重磁位场转换计算中迭代法的综合分析与研究

处理转换计算在重磁资料解释中发挥着重要的作用,但一些计算如向下延拓、化极等有时是很不稳定的,在频率域中则表现为其转换因子具有明显的放大作用,所以其FFT理论计算结果是不稳定的.因此,很多研究工作都是围绕增加计算的稳定性、提高计算效果进行的,其中迭代法是近来在研究中受到普遍重视的方法技术,并取得了较好的成果.但也存在对迭代法研究还不够深入,对其存在的缺点认识不够充分、客观等问题,例如,迭代法进行延拓及化极等计算时,对一些具体应用虽能在一定程度上获得较好的计算结果,但却存在计算结果并不会随着迭代次数的增加而得到持续改善的问题,对于原本不稳定的计算,迭代法在迭代次数比较大时,所得的结果依然是不稳定的.为此,本文在对迭代法进行分析研究的基础上,进一步推导了迭代法的通式,并分析了对迭代法收敛性影响的各种因素.分析结果表明:迭代法收敛到FFT理论直接计算结果的决定因素是计算过程中如何选择原始数据到目标数据的映射函数;在选择了合适的映射函数的情况下,迭代次数不仅仅是决定计算成本,而是决定结果好坏的关键因素;增加迭代次数虽然能够使计算收敛到FFT直接计算理论结果,但如果该理论结果本身就是不稳定的,则迭代法计算如果收敛,也是收敛到一个不稳定的结果.所以针对位场处理转换中一些不稳定计算采用迭代法,并没有从根本上解决计算的不稳定性问题.     

物理过程参数化方案对中国夏季降水日变化模拟的影响

国际耦合模式比较计划第六阶段(CMIP6)新增的高分辨率模式比较计划(HighResMIP)首次提供全球高分辨率(25—50 km)多模式集合的气候模拟试验结果。利用8个CMIP6 HighResMIP模式评估了高分辨率全球气候模式对青藏高原夏季小时降水与极端降水的模拟能力,结果表明：CMIP6高分辨率模式高(低)估了青藏高原地区的降水量和频率(强度),过多的降水量主要来自模式对降水频率的高估,尤其是弱降水(< 2 mm·h^-1)的发生频率。模拟偏差与地形海拔密切相关,偏差大值区主要位于高原南坡和东坡陡峭地形区。模式不能准确再现降水量与海拔之间的关系,高(低)估了高(低)海拔地区的降水量。模式低估了降水强度随海拔升高而降低的变化速率。在日变化方面,模式能够模拟出青藏高原降水傍晚至午夜的峰值特征,但明显低估了降水的日变化振幅。在小时极端降水方面,模式低估了高原区域平均极端降水第95百分位数阈值,仅为观测值的57%。

     

补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

陇中黄土高原丘陵沟壑区不同植被恢复模式下次降雨产流产沙特征

基于甘肃省清水县汤峪河径流小区2015—2017年的观测数据,研究不同植被恢复模式条件下坡面次降雨入渗、产流产沙特征。结果表明:不同植被恢复模式条件下的土壤入渗量与降雨强度呈二次函数关系,存在入渗量达到最大值的临界降雨强度。入渗速率与降雨历时可以用幂函数关系表达,符合考斯恰可夫入渗模型。不同植被恢复模式条件下的产流率在0.003 3～0.003 6 mm·min-1之间,相对裸地的减流率为54%～58%。产流率与降雨强度之间呈二次函数关系(R2>0.88),产流率的主要影响因素是降雨强度。径流含沙量平均值乔灌混合区(3.13 g·L-1)>灌木林(2.95 g·L-1)>乔木林(2.79 g·L-1)>草地(2.58 g·L-1),径流含沙量与降雨强度呈线性递增函数关系。裸地的产沙量显著高于各植被小区(P<0.05),是各植被小区的43～57倍,各植被小区的减沙率在93%～94%之间,减沙效益高于其减流效益。各植被坡面土壤流失量与降雨侵蚀力呈线性递增函数关系;产流率与侵蚀产沙率之间呈极显著正相关关系(P<0.01),二者间可采用二次函数关系表达。本研...     

位场梯度的归一化标准差方法在地质体边界定位问题中的应用

在综合分析前人研究成果基础上,采用位场梯度的归一化标准差(NSTD)方法来识别地质体边界.通过理论模型实验讨论了窗口大小的选择对于计算结果的影响,并验证了NSTD方法不仅可以识别单一地质体边界,而且还可以识别不同埋深的、形态各异的相邻地质体边界.利用此方法对大巴山及其前缘地区的重力数据进行了处理,并与水平总梯度模(hgm)、总梯度模(AS)、局部波速(LW)以及Tilt水平总梯度模(Tilt-hgm)方法的应用效果进行了对比分析.研究结果提供了大巴山及其前缘地区的构造分区信息,对于油气远景预测具有重要的指导意义.     

基于差分干涉雷达的汶川地震同震形变特点

地震同震形变场是认识震源机制、确定发震断裂、分析发震断裂活动方式、评估震害损失及揭示未来发震趋势的重要依据。但在龙门山地区复杂的自然环境中, 现有的地震同震形变场获取方法(GPS测量、构造形迹分析、震源机制解结合数值模拟计算)难以快速、全面的捕捉到这种信息, 这突显出差分干涉雷达技术(D-InSAR, Differential Interferometry SAR)在同震形变场研究中的优势。本文介绍了其观测原理, 在地震同震形变监测中的研究现状, 并重点分析了已获取的D-InSAR汶川地震同震形变场的观测结果。基于差分干涉雷达观测及相关资料, 对龙门山断裂在本次地震中活动性进行了初步分析, 结果表明干涉雷达形变图像与地表破裂范围、逆冲角度的变化、上下盘升降关系及大型余震的展布存在很强的相关性。      

甚高频时岩石介电常数的测试方法和特性

本文讨论20—200MHz频率范围内平行板电容法测量岩石介电常数和电导率的方法,指出常用的静电场公式的局限性,给出较精确的似稳场动态公式,并讨论20—50MHz频率范围内的岩石介电特性。     

微波消解-电感耦合等离子体质谱法测定土壤和沉积物中痕量稀土元素

建立了微波消解-电感耦合等离子体质谱法测定土壤和沉积物中15种痕量稀土元素的分析方法。研究了溶样试剂、微波消解程序、标准溶液配制、质谱干扰与内标元素对稀土元素测定的影响。加入氢氟酸能有效打开样品,以HNO₃-HF-H₂O₂混合酸溶解样品,稀土元素的溶出率较高。采用模拟土壤、沉积物中稀土元素天然组成比值的校正溶液,对稀土元素间的干扰具有明显的抑制作用。通过测定单个La、Ce、Pr、Nd和Ba的氧化物及氢氧化物产率,计算出等效干扰浓度,进而校正多原子离子干扰。利用¹⁰³Rh内标校正系统,有效地抑制了分析信号的动态漂移。方法检出限为1.2~7.1 ng/g,精密度(RSD)≤5.3%(n=6),加标回收率为86.1%~110.1%。使用土壤、沉积物标准物质进行验证,测定结果与标准值相符。建立的方法样品处理程序简单快速、线性范围宽、分析重现性好、结果准确,适用于大批量地质样品的分析。     

真实基流中非线性Rossby波演变特征(二):能量、结构演变及初始场影响

针对非线性的准地转正压位涡方程,利用自行设计的差分格式和高斯函数拟合得到的真实基流分布,数值研究了线性和非线性Rossby波流场结构和总能量的演变以及初值对总能量演变的影响.发现在非线性的真实基流中,线性和非线性Rossby波的相对总能量出现振荡型增长或衰减,非线性波动的振荡周期明显小于线性波动,非线性项不仅抑制能量的快速发展,也抑制能量的快速衰减;线性Rossby波的流场结构一直是以标准的正弦(或余弦)方式演变,而非线性Rossby波则是以类似椭圆余弦波的方式演变,并会出现长波调整现象;非线性Rossby波对初始场和基流结构具有强烈的依赖性,在其总能量演变过程中,可出现类似混沌的现象.     

利用远震研究四川盆地及其周缘瑞利面波相速度和方位各向异性

基于四川盆地及周边的245个宽频带台站2010年9月-2014年9月期间的远震记录,提取双台路径瑞利面波相速度频散资料,反演得到四川盆地20~120 s的高分辨率瑞利面波相速度及各向异性空间分布.在丰富区域地球物理基础数据的同时,结合已有研究成果对地壳上地幔变形耦合进行探讨,结果表明短周期（20~30 s）的相速度分布与四川盆地的地质构造特征相吻合,作为川滇地块、松潘-甘孜地块和四川盆地之间的边界——龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用;松潘-甘孜地块,特别是川滇地块中下地壳普遍表现为明显的低速异常,表明中下地壳相对软弱;而四川盆地的中下地壳整体呈现相对高速,表明四川盆地具有相对坚硬的中下地壳.研究区域东南角接近北扬子地块与南扬子地块的缝合部位,呈现高速异常.四川盆地南部和东南邻区不同周期均具有较强的各向异性,且快波方向较为一致,反映这些地区不同深度变形耦合较好.四川盆地西部、北部及东北部邻区,不同周期的各向异性快波方向变化较大,不同深度变形耦合较差.这些特征与绕喜马拉雅东构造结的物质流动被扬子地块的高速地壳阻挡的宏观认识基本一致.