面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.

     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法（线性方法）与改进遗传算法（非线性方法）相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.

     

基于多轨道卫星观测数据分析尼泊尔地震长波辐射特征

综合已经在红外异常提取中应用的涡度和RST（Robust Satellite Technique）算法优点,提出了红外异常指数算法.并基于长时间尺度的中国静止卫星FY-2D和美国极轨卫星NOAA长波辐射数据,应用RST和异常指数算法,分别对2015年4月25日尼泊尔M_S8.1和5月12日M_S7.5地震前后卫星长波辐射变化特征进行了分析,开展了多轨道、多时空分辨率长波辐射同步地震热红外特征研究.结果表明,运用RST算法,两次地震前后,未能在震中周围发现明显的长波辐射异常.运用异常指数算法：（1）对于NOAA卫星,4月15日在M_S8.1地震震中以西出现热红外异常,到4月24日震中以西约100 km处出现异常最大值,随后逐渐消失.5月10日在M_S7.5地震震中以东约200 km发现异常;（2）在NOAA卫星长波辐射异常发现最大值当日,采用FY-2D卫星每3 h的数据分析可发现红外异常的动态演化过程,弥补NOAA卫星分辨率不足.以上结果为利用多轨道卫星监测地震热辐射变化提供了依据.

     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明：新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.

     

水平层状非均质TI地层中仪器偏心情况下三维感应测井响应高效数值模拟与响应特征分析

本文针对中国石油集团测井有限公司（CPL）三维阵列感应测井仪器（Three Dimensional Induction Tool,TDIT）的真实结构与参数,利用2.5维数值模式匹配技术研究建立水平层状非均质TI地层中仪器偏心情况下三维感应测井响应的高效数值模拟算法.首先,利用2.5维数值模式匹配技术给出水平层状非均质TI模型中偏心磁流源并矢Green函数的半解析解.在此基础上,结合三维阵列感应测井仪器参数以及测量过程中出现的仪器偏心与旋转等实际情况,通过叠加原理推导出该仪器响应的有效计算公式.最后,通过数值模拟结果分析考察仪器偏心、仪器自旋角等参数对测井响应的影响以及变化特征,为多分量感应资料处理与解释提供理论基础.

     

异常结构三维综合系统的开发和应用

地球化学异常结构三维综合系统的开发和应用是为了满足地球化学多维异常体系的研究而逐步展开的。应用于深部找矿预测的多维异常结构方法,从钻孔岩芯直接采样,获得了具有真三维属性的地球化学数据,研究者迫切需要在三维场景下对所获三维地球化学数据进行解读和表达。异常结构三维综合系统(3D-MDAS)是利用三维可视化和地质体三维建模技术,在MapGIS K9平台上进行二次开发获得的。本系统实现了地上、地下数据模型的一体化显示,实现了实体模型与属性模型的耦合一体化显示,在此基础上还可以对异常结构进行空间分析,建立三维地球化学找矿模型。笔者详细阐述了三维异常结构3D综合系统的三维建模思路和技术,并用安徽马头矿区的实际应用成果给以说明。     

油气运聚影响下三维地震精细解译的砂岩型铀矿控矿因素分析——以松辽盆地大庆长垣南端铀矿为例

松辽盆地大庆长垣隆起区砂岩型铀矿资源丰富,但成矿条件复杂。以大庆长垣南端四方台组砂岩型铀矿为研究对象,基于三维地震精细解译成果,对长垣隆起区铀矿的分布特征、富集主控因素及成矿模式开展研究。结果表明：①三维地震资料是刻画断层与砂体分布的有效手段;②砂岩型铀矿主要富集在四方台组一段的河道砂体中,与沟通深部油气储层的输导断层邻近;③河道砂体、油气输导断层及油气圈闭是控制砂岩型铀矿富集的3个主要因素;④长垣隆起区的砂岩型铀矿是含铀地下水与油气运聚协调配合局部富集的结果,建立了输导断层与油气圈闭联合控矿模式。研究成果为含油气盆地砂岩型铀矿勘探提供了新的思路。     

矿产资源潜力评价方法对比及其发展趋势探讨

矿产资源潜力评价是预测一个地区矿产资源找矿潜力的评价方法,先后经历3个阶段：探索和应用阶段、快速发展阶段及信息化阶段。笔者梳理了“三步式”、证据权、预测普查组合、成矿系列、地质异常、综合信息预测、地球化学块体和非线性预测等8种矿产资源潜力评价理论或方法,通过相关应用实例阐述其实用性,分析了国内外应用潜力评价方法获取的研究成果,总结其发展趋势。其中,“三步式”和证据权是目前使用较为广泛的方法,成矿系列则是国内研究的热点。随着科学技术的进步以及数学地质的不断发展,矿产资源潜力评价已步入信息化阶段,三维深部预测将是未来潜力评价重点发展方向。     

城市场景结构感知的网格模型简化算法

针对二次误差测度(QEM)网格简化算法全局几何特征信息损失严重的缺点,提出一种具有结构感知功能面向城市三维模型重建的网格简化算法。该算法顾及城市影像中平面结构特征,以代理平面为全局特征约束条件,使模型在简化过程中全局结构特征更多地被保持,以利于多层次细节模型(LOD技术)、网格优化提速等模型后续操作。以倾斜摄影获取影像生成的初始三角网格模型为试验数据,采用所给算法对其进行网格简化并与QEM算法进行对比。结果表明:所给算法简化精度及简化效率均优于QEM算法。     

多视密集匹配并行传播GPU-PatchMatch算法

针对多视密集匹配的效率较低的问题,提出了GPU-PatchMatch多视密集匹配算法。该算法使用GPU提高PatchMatch的计算效率;同时充分利用稀疏场景信息,对深度信息进行规则初始化;为提高传播效率,使用了金字塔红黑板并行传播深度信息。最后在DTU、Strecha和Vaihigen数据集上进行了试验,并与常用的多视密集匹配算法进行对比。试验结果表明,本文算法在重建效率上有较大提高,与CPU算法(PMVS、MVE、OpenMVS)相比有7倍以上提升,与GPU算法相比也有2.5倍以上提升,表明本文算法的有效性。