地面核磁共振找水反演

简要介绍了地面核磁共振找水的正演理论方法,采用高斯求积并结合连分式展开的方法对核磁共振核函数中包含双重贝塞尔函数的积分核进行了数值积分,进而计算出精度较高的核函数值,在此基础之上,讨论了导电性对地面核磁共振信号的影响.基于奇异值分解算法,对核磁共振找水理论模型进行了反演研究,在不导电的情况下,重构出了比较理想的含水率分布.总结并改进了模拟退火算法,用改进的模拟退火算法分别对覆盖层高阻和导电条件下的核磁共振人工合成数据进行了反演,试验结果表明,利用改进的模拟退火算法反演地面核磁共振是可行的,而且反演结果较稳定,收敛速度较快.     

基于RRI反演的高频大地电磁测深在深边部矿产勘探中的试验研究

随着我国国民经济的快速发展,对各种矿产资源的需求量越来越大,经过多年对近地表的矿产开采之后,找矿的目标逐渐向矿区深边部发展,本文阐述采用EH-4电导率成像系统的高频大地电磁法,并且利用RRI反演,能够准确快速地探测大深度范围内的矿化异常,对矿区的深边部矿产的赋存状态进行超前地宏观预测,并将野外的成功实例进行分析.     

煤层气解吸的温度效应

为了探索温度对煤层气解吸效果的影响,通过对沁水南部寺河煤矿3号煤的吸附/解吸实验,经过数学分析,获得该煤层气吸附/解吸的数学式与等量吸附热表达式。分析认为：随着储层温度升高,饱和吸附量减小;解吸过程滞后于吸附作用的原因主要是非物理吸附常数 c值作用,随着储层温度增高,c值增大;降压解吸伴随吸热(温降),同时阻碍了持续解吸过程。通过提高储层温度来促进解吸过程,在理论上可行,工程中应该受到重视。     

利用树轮资料重建天山中段南坡巩乃斯地区公元1777年来平均最高气温变化

利用年轮平均灰度年表重建了天山巩乃斯地区自公元1777年来当年5月至8月平均最高气温,方差解释量达41%(调整自由度后为39%).过去232年来,巩乃斯地区平均最高气温变化存在2个偏冷阶段和1个偏暖阶段。重建序列具有154年、77年、2.7年和2.3年的准周期, 并且在1842年前后、1880年前后、1923年前后发生了突变。与北半球平均最高气温网格点资料和多种涛动指数的相关性分析表明,巩乃斯地区平均最高气温变化不仅反映了区域气候变化,对于大范围气候变化也存在一定程度上的响应。     

大地电磁测深方法在内蒙古西部银根-额济纳旗盆地石炭系—二叠系油气地质调查中的应用

对银-额（银根-额济纳旗）盆地以石炭系—二叠系为主要目的层开展了油气地质调查实施的7条大地电磁测深剖面资料的处理和定性分析,在综合信息反演和成像的基础上建立了重磁电综合反演剖面图。以钻井、地震和区域地质调查资料为约束条件,确定了电性层标定和解释的基本原则。利用大地电磁测深反演电阻率低阻标志层,准确地界定了中生界的埋深和厚度,次高阻层半定量地刻画出了石炭系—二叠系的分布,对剖面的地质层位进行了追踪划分,对研究区地层、电性特征进行了系统的总结。研究成果表明,大地电磁测深法是获取中生界、古生界分布的重要手段。在银-额盆地复杂地质条件和沙漠区恶劣的施工条件下,大地电磁测深法显示出的优势和良好的勘探效果尤为明显,在该区石炭系—二叠系油气地质调查中起到了重要作用。     

地震和潮汐对气候波动变化的影响

2004年12月26日印尼地震海啸后,全球低温冻害和暴雪灾害频繁发生."潮汐调温说"和"深海巨震降温说"是一种合理的解释.根据"潮汐调温说"和"深海巨震降温说"理论,2005年以后全球气温将因为地震海啸和强潮汐南北震荡而降低.2009年11月至2010年1月低温暴雪袭击北半球,西方科学家也承认2000~2010年气候的自然变化减缓全球气候变暖效应这一客观事实.潮汐震荡可以解释全球气温的准60a变化,海洋及其边缘的强震能够将深海冷水翻上表面,使全球气候变冷.潮汐和太阳黑子活动不仅有相同的变化周期,而且都和气温变化有很好的对应关系.研究表明,在太平洋十年涛动冷位相时期,强震与低温冻害频繁发生.     

倾子资料三维共轭梯度反演研究

在对倾子响应和共轭梯度算法深入分析的基础上,我们实现了倾子资料三维共轭梯度反演算法.基于倾子资料的三维共轭梯度反演研究,探讨了利用倾子资料进行三维反演定量解释的方法.通过对理论模型合成数据进行反演试算,验证了所实现的倾子资料三维共轭梯度反演算法的有效性和稳定性.该反演算法可用于对大地电磁测深和地磁测深(地震地磁台站进行的三分量地磁观测)所整理出的倾子资料进行三维定量反演,获得地下三维模型的电阻率结构.     

双参数混合正则化方法及在电导率反演成像中的应用

通过引入带有二阶正则算子的正则化项,建立了一种双参数混合正则化方法.为确定最佳正则化参数,这里主要应用L-曲线法、偏差原理和广义交叉校验准则的优化组合来确定.首先对理论模型进行了数值模拟,通过与截断奇异值分解法、共轭梯度法及标准Tikhonov正则化法的结果比较,表明该方法不仅精度高,而且对于数据的随机扰动具有稳定性.然后将此方法应用于对电导率反演成像的计算中,得到的电导率反演成像精确细致可靠,也符合实际情况.这表明该方法求解地球物理反问题时,在增强反演的稳定性及减少多解性的同时,还能进一步提高反演计算速度.     

2018年9月8日墨江地震及周边地区构造应力场特征分析

2018年9月8日,云南省墨江县发生M_S5.9地震并伴随一系列余震,探究该地震周围的应力场对于理解该地震的发生机制和后续地震的发展趋势具有着重要的参考意义.本研究收集了震源及其邻区中前人研究和Global CMT所给出的震源机制资料,对该地区进行了构造应力场反演,并同时利用反演得到的应力张量模拟墨江地区的震源机制解表现.结果表明：①在应力轴整体分布上,自西向东σ₁轴（压轴）从NNE-SSW向逐渐转向NNW-SSE向,σ₃轴（张轴）从WWN-EES向逐渐转向WWS-EEN向,张轴呈弧形分布,压轴呈放射状分布.②在应力轴倾伏角上,研究区域内的压应力轴和张应力轴倾伏角都比较小,即两轴均接近水平.③R值分布大体是在东南部相较于西北部大,结合当地地质背景分析得到,物质逃逸自西北向东南呈逐渐变缓的趋势.④利用反演得到的应力张量和应力状态计算墨江地震震源区的相对剪应力和相对正应力大小.由此推测,墨江地震恰好发生在相对剪切应力值和相对正应力正值最大的节面上.从而可以确定墨江地震的发震节面的基本参数：走向216.32°,倾角86.91°,滑动角0.27°,相对剪应力值0.9,相对正应力值0.3.本研究为此次墨江地震的发震背景和地震动力学研究提供了基础性资料.

     

2013年南斯科舍海岭MW7.8地震的多点震源机制反演

2013年11月17日,在南极南奥克尼群岛北、南极板块与斯科舍板块之间发生了一次M_W7.8级地震(2013年南斯科舍海岭M_W7.8地震),我们利用全球分布的长周期和宽频带地震记录反演确定了这次地震随时间和空间变化的震源机制,验证了提出的一种多点震源机制反演的新方法.首先利用长周期记录的W震相反演了这次地震的矩心矩张量解并利用体波提取了视震源时间函数,同时利用台阵反投影技术从宽频带记录中获得了这次地震的高频源的时空分布,然后基于矩心矩张量解、视震源时间函数以及高频源的时空分布,实现了采用新方法对2013年南斯科舍海岭M_W7.8地震的多点震源机制反演.矩心矩张量解表明,地震矩心在44.50°W/60.18°S,矩心深度19 km,半持续时间49 s,释放标量地震矩4.71×10²⁰ N·m,发震断层走向104°,倾角54°,滑动角8°.视震源时间函数清楚地揭示了地震矩随时间变化的方位依赖性,总体上可以将时间过程分为前60 s和后50 s两个阶段,但前60 s可细分为两次子事件.根据台阵反投影结果,这次地震为沿海沟从西到东的单侧破裂,破裂长度达311 km,可以分为5次子事件,能量释放的峰值点依次为13 s、30 s、51 s、64 s和84 s,平均破裂速度分别为0.6 km·s^-1、2.6 km·s^-1、2.3 km·s^-1、2.8 km·s^-1和3 km·s^-1.多点震源机制反演显示,5次子事件的矩震级分别为M_W7.57,M_W7.48,M_W6.80,M_W7.53和M_W7.08,半持续时间依次为21 s,17 s,6 s,16 s和8 s,走向分别为95°,105°,81°,98°和98°,倾角依次为57°,49°,86°,46°和64°,滑动角-9°,1°,-17°,13°和-4°.这些在震源机制、能量释放以及持续时间方面的变化都是当地构造和应力环境复杂性的反映.