大地电磁三维数据空间反演并行算法研究

目前大地电磁三维反演实际应用的主要问题是计算效率低.在对大地电磁三维数据空间反演算法进行深入分析的基础上,本文提出了基于频点和矩阵划分的大粒度并行反演方案和具体实现步骤,并在曙光TC5000A高性能计算平台上实现了基于MPI的大地电磁三维数据空间反演并行算法.该算法实现了包括三维正演、灵敏度矩阵、叉积矩阵以及模型改正量的并行执行,不仅计算效率高,而且每个节点机上灵敏度矩阵的存储空间只需原来微机上的2/N(N是参加并行计算的节点机个数),大大地减少了内存开销.通过两个理论模型合成的数据对实现的三维数据空间反演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维数据空间反演并行算法是可行的、高效的,与单机相比,不仅可以提高运行速度,缩短计算时间,而且还可以扩大计算规模,极大地推动了大地电磁三维反演的实用化.     

基于褶积模型的地球物理反演模型增强

地球物理数据在采集和处理过程中,由于存在噪声、模型误差、以及数据离散化误差等系统误差,导致了异常体边界模糊和模型分辨率降低等一些不可避免的不良系统退化效应的产生.本文提出了一种新的地球物理反演模型增强方法,通过消除反演估计模型中的系统误差,压制模型中的不良系统退化效应,增强反演模型的分辨率.文章从理论上分析了数据中存在的系统误差对模型求解的影响,提出了一个新的系统误差褶积退化模型,并根据该模型提出了一种基于混合范数总变分正则化的盲反褶积模型增强算法.最后,文章通过1D线性反演增强试验和2D大地电磁反演增强试验,验证了所提出的地球物理系统退化模型的正确性,以及盲反褶积增强算法的有效性.试验结果表明,方法可以有效地提高反演参数模型的分辨率.     

二维波动方程速度的正则化-同伦-测井约束反演

针对二维波动方程反问题,将大范围收敛的同伦方法引入速度参数的反演过程中,并将其与求解不适定问题的Tikhonov正则化有机结合,提出了一种新的、特别适用于非线性的、不适定的、多极值的地震勘探反演问题的反演策略：正则化－同伦方法. 为了充分利用测井资料和地震资料的互补特征,进一步提高反演分辨率并压制噪声,设计了正则化-同伦－测井约束联合反演方法. 大量数值试验结果表明了这两种方法的有效性.     

基于高斯牛顿法的频率域可控源电磁三维反演研究

三维反演解释是电磁法勘探发展的重要趋势,而如何提高三维反演的可靠性、稳定性和计算效率是算法开发者们目前的研究重点.本文实现了一种频率域可控源电磁(CSEM)三维反演算法.其中正演基于拟态有限体积法离散化,利用直接矩阵分解技术来求解大型线性系统方程,不仅准确、稳定,而且特别有利于含有大量发射场源位置的CSEM勘探情况;对目标函数的最优化采用高斯牛顿法(GN),具有近似二次的收敛性;使用预条件共轭梯度法(PCG)求解每次GN迭代所得到的法方程,避免了显式求解和存储灵敏度矩阵,减小了计算量.以上这些方法的结合应用,使得本文的三维反演算法准确、稳定且高效.通过陆地和海洋CSEM勘探场景中的典型理论模型的反演测试,验证了本文算法的有效性.     

基于拟态有限体积法的频率域可控源三维正演计算

大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法.本文通过求解离散化的三维电场矢量Helmholtz方程,实现了基于有限体积法的频率域可控源电磁(CSEM)三维正演算法.为模拟具有强电性差异的三维电性介质,该算法采用拟态有限体积法(MFV)对Maxwell方程组进行离散化;另外,为获得稳定、高精度的正演数值结果,采用直接矩阵分解技术来求解离散所得到的大型稀疏线性方程组.对于具有多个发射源的CSEM测量来说,一次矩阵分解结果能够用于同频率下所有场源的正演计算.为降低场源奇异性及边界条件对数值精度的影响,采用虚拟场源校正技术,避免了散射场公式中在构建场源项时所需的大量时间.对于具有多个频率的CSEM的模拟计算,采用分频并行策略来加快三维正演计算.最后,通过与一维层状模型及三维模型的数值结果的对比验证了本文所开发的正演算法对频率域CSEM模拟计算的准确性及有效性,表明该正演算法能够有效应用于三维介质的数值计算.另外,对于多频率CSEM的并行测试结果表明基于分频并行策略的并行计算能够显著地降低正演计算时间.     

气溶胶光学厚度估测中的LASSO特征选择方法

气溶胶光学厚度估测中通常利用遥感信息构造的多种特征属性作为输入,然而,这些属性中常常存在数据噪音、相互关联性和缺失值,从而降低了估测精度和估测强健性。针对这个问题,基于最小绝对收缩和选择算子（least absolute shrinkage selection operator,LASSO）方法和气溶胶光学厚度反演的先验知识,提出了一种针对遥感卫星观测的高维数据进行特征选择的方法,利用2009年4月2日至2011年4月1日2 a间与全球197个气溶胶地基自动观测网站点时空同步的MODIS（moderate-resolution imaging spectroradiometer）遥感数据,采用常用的人工神经网络作为估测模型进行实验分析,表明该方法能结合反演先验知识对多种异质遥感属性进行分组,通过组间迭代保留关键特征,去除冗余属性,有效进行特征选择,从而显著提高气溶胶光学厚度的估测精度。     

辫曲转换与共存的主要影响因素及对古代河流沉积环境恢复的启示

河流辫曲转换与共存是自然地理学、水力学、河流沉积地质方面的研究热点问题,也对古代河流沉积环境恢复与储层预测有重要借鉴意义.首先对河流辫曲转换与共存特征的主要影响因素进行了探讨,指出构造变动与地貌单元、坡度（坡降）的陡缓差异、物源远近与水动力条件、气候变化与植被发育状况、海/湖平面变化等主要因素控制着辫曲河型转换过程,其中构造变动与地貌单元是最为关键因素;然后总结了辫曲转换与共存这一理念对古代河流沉积恢复研究的4点启示;再以鄂尔多斯盆地大牛地气田中二叠统下部的下石盒子组为例,结合测井曲线特征、地震属性与砂岩厚度分布规律,把辫曲转换与共存的理念应用于古代河流沉积环境恢复之中,最终再现了气田区H21砂层河流沉积分布格局,明确了该区北辫南曲、辫曲转换与共存的规律,指出在辫曲转换地带与相邻的曲流河发育的区域更易于产生废弃河道.      

衡阳山洪地质灾害气象预警系统的研发与应用

衡阳受气候和地理条件的影响,是湖南省内发生山洪、地质灾害较多的地区之一,严重威胁人民群众生命财产安全,制约全市社会经济发展。该文介绍了衡阳山洪地质灾害气象预警系统(以下简称MGHWS)的研发与应用情况,并对其中的关键技术及其实现方法进行了分析与论述。系统的建成并投入应用,使国家"山洪项目建设"的综合效益得到充分发挥,提高了衡阳应对由强降水所引发的山洪地质灾害的监测、预警及服务能力,社会、经济与生态效益显著。     

三度体磁异常的不定积分变换法研究

为发现和圈定干扰覆盖下的不明显异常，作者从三度磁性体的磁异常公式出发，推导磁异常的不定积分变换结果并由参量图方程绘得不定积分变换的磁异常参量图。本文针对2种规则形态磁性体推导其特征点规律，并进而用不同的方法求取磁性体的深度和磁性参数。在不同的空间几何参数条件下，讨论不定积分变换后参量图的曲线变化规律并总结其异常特征。提取判断产状的要素。该方法能简单有效地发现一些重要的地质信息，所得到的空间参数与变换异常间的特征点关系将有助于磁异常的分析。     

基于流体饱和多孔隙介质波动方程的时间推移地震反演

针对二维流体饱和多孔隙介质波动方程,设计了自适应同伦法和小波-自适应同伦法,在此基础上,将局域化思想引入时间推移地震反演,并使用小波-自适应同伦法和自适应同伦法反演两次监测数据,提出了一种适用于时间推移地震勘探的局域化反演算法.为了验证新方法是可行的,文中把它应用于同一地点不同时间的观测模型,通过数值实验表明了这种方法的有效性.