改进麻雀搜索算法的自然电位数据反演

自然电位法是一种成本低廉,野外观测操作简便的天然源地球物理勘探方法。自然电位数据反演具有病态和非线性的特征。常规反演算法分为局部优化类和全局搜索类,但基于梯度运算的局部优化算法难以求得全局最优解且反演效果依赖于初始模型的构建,而传统全局搜索算法又存在收敛速度慢,易陷入局部极值和不稳定的缺点。基于此,对一种新的全局优化策略(麻雀搜索算法)进行改进,通过混沌映射叠加反向学习策略初始化麻雀种群,再依随机概率使用Levy飞行策略更新麻雀个体位置来进一步提升算法对于解空间的探索能力和增强算法跳出局部极值的可能性。将改进前后的麻雀搜索算法分别应用于合成自然电位数据(不含与含10%、30%的随机噪声)与来自印度和法国的实测数据的反演解释中以对比检验改进算法的反演效果。理论测试结果表明：麻雀搜索算法(SSA)在无噪声干扰下的垂直圆柱和倾斜板模型数据反演误差为0.42%和0.25%,相同情况下改进麻雀搜索算法(ISSA)的反演误差为0.06%和0.07%,改进后算法拟合精度提高到3～7倍,对比目标函数收敛曲线图中ISSA的收敛速度与收敛精度都要明显优于SSA;SSA反演参数的稳定性、精度和异常响应曲线拟...     

地球物理非线性反演模拟退火法的改进

在地球物理非线性反演方法中，模拟退火法是一种较先进的启发式蒙特卡洛（ＭｏｎｔｅＣａｒｌｏ）方法．但是，在处理实际资料时，该方法存在着计算效率不够高的缺点，有时还会失效．为此，从模拟退少法的关键问题──最低温度的选择入手，根据模拟退火法与统计力学的吉布斯－马尔柯夫（Ｇｉｂｂｓ－Ｍａｒｋｏｖ）模型之间的关系导出临界温度的近似表达式；用此式分析目标函数超曲面形状对模拟退火法计算的影响；提出利用模糊先验信息确定最低温度、改造目标函数等改进措施．     

中国及其邻区周围板块作用力的研究

给出了如何根据观测的应力方向,利用平滑方法处理为在规则格点上均匀分布的数据的过程,并利用新编的中国及其邻区地震震源机制解推断的观测应力方向,反演了周围板块作用力的相对大小．结果表明,在周围三大板块中,以印度板块的作用力最大;其次是太平洋板块;而在菲律宾海板块边界上,作用力大小差异较大,其中又以与台湾相邻的边界段受力最大,琉球岛弧受力最小．     

改进的约束变密度界面反演策略及其应用

计算密度分界面的起伏变化在区域地质构造研究和石油矿产资源勘探中具有重要意义.已有密度界面反演方法更多侧重约束变密度界面反演算法,而对约束信息的准确性、研究区横向密度变化往往考虑不足,影响了最终反演结果的可信度.本文在变密度界面正反演算法基础之上,结合实际需求,提出已知深度信息约束下的变密度界面反演策略.该策略主要包括变密度约束反演算法、已知深度约束信息校验和分区变密度模型三个方面.其中反演算法提供了带已知深度约束信息的密度界面迭代反演方法;约束信息校验用于评估约束信息精度,通过调差降低约束信息的系统误差;在反演过程中引入水平密度分区以应对不同构造背景密度界面模型,提高反演结果的可信度.最后将本策略应用于南海莫霍面深度反演计算中,结果显示借助已知约束信息,利用分区密度模型能够获得更为可信的深度反演结果,验证了该策略的正确性.     

2013年芦山Ms7.0地震前重力场源变化分析

陆地重力观测获取的时变重力信号是研究地壳内部孕震区介质变化的重要依据。利用贝叶斯平差方法处理了川西地区2010—2012年的流动重力观测资料,发现相对重力仪的格值系数在两年内变化可达2×10^-4,而且仪器存在较小的非线性漂移率;重力场累积动态变化图像显示,2013年芦山MS7.0地震发生在重力梯度带零值线与龙门山断裂带交会部位附近。进一步采用等效源方法划定了测网内以雅安、康定和石棉为中心的场源分辨能力较高的区域,并反演了地震前的地壳视密度变化。场源视密度动态变化结果表明:2013年芦山MS7.0地震前2年,在龙门山断裂带的南侧和鲜水河断裂南段的东侧区域形成了范围较大的等效视密度增加区,而在小金至康定一带和都江堰周边地区也存在小规模的视密度增加。本文基于陆地时变重力测量结果反演得到芦山地震前的地壳内部等效视密度变化,可为研究壳内深部孕震环境变化特征提供参考。     

川滇地区中小地震震源机制解及构造应力场的研究[

利用目前国际上较为流行和普遍接受的基于P波初动符号和振幅比HASH方法,研究给出川滇地区2003年1月1日-2012年12月31日1 893次M3.0以上地震震源机制解。选择1 651个可靠的中小地震震源机制解,采用基于中小地震震源机制解的Hardebeck和Michael的阻尼区域应力场反演方法,研究区域的水平最大主压应力方向。结果表明：川滇菱形块体南段红河断裂带尾部左右两侧应力方向相同,水平最大主压应力NW-SE向;川滇菱形块体内部以丽江—小金河为界呈现出不同的应力状态,位于分界线以南的滇中块体水平最大主压应力NW-SE向,而在分界线以北取向基本上沿NS方向;从青藏高原内部到川滇菱形块体东边界应力方向整体有一定的顺时针旋转趋势;块体东边界从北部的NNE-SWW向逐渐过度到南部的NW-SE向。滇西及滇西南水平最大主压应力方向与构造方向基本平行,为近SWW-NEE向。     

青藏块体东北缘弹性块体边界负位错反演与强震地点预测

利用青藏块体东北缘1993——1999年、1999——2001年和2001——2003年3个时段的GPS水平运动速度场成果,通过研究建立的弹性块体边界负位错模型反演,给出该区构造块体及其边界断裂运动与变形差异、构造应变场时空分布、高应变能积累部位以及强度的定性分析和定量描述. 在充分考虑区域构造及块体应变的前提下,研究与该区6级左右及6级以上地震地点预测相关的背景性前兆的共性特征.      

南海珊瑚礁光学浅水海域多光谱遥感水深反演

海洋水深信息对研究珊瑚礁海域资源与环境具有重要作用。南海珊瑚礁海域测深数据受多种条件限制施测困难,在时间与空间方面数量非常有限。文章针对南海岛礁海域以I类水体为主导的海水光学特性,以南沙群岛库归沙洲海域为例,使用Sentinel-2多光谱卫星遥感影像和同期过境的MODIS卫星数据,构建底质光谱,采用半分析半经验模型计算海水表面遥感反射率与海水叶绿素浓度,通过对数比值模型进行该地区光学浅水海域遥感水深反演分析,并进一步通过多时相反演水深融合提升精度。经与多波束实测水深数据验证,研究区域反演水深总体均方根误差和平均相对误差分别为2.68 m 和9.99%。该方法通过叶绿素浓度推演部分海水光学特性,可以从多光谱卫星影像中快速获取南海岛礁光学浅水海域初步水深信息,供相关海洋领域分析与应用。     

重磁三维反演伊犁盆地中部密度和磁性结构

伊犁盆地中部区域地处多个构造单元的交汇区,构造复杂,缺乏高精度地质资料,制约了对盆地构造演化的进一步认识.本文联合使用多种先验信息作为约束条件,利用高精度重磁数据反演获得了研究区地下海拔-10 km以浅的三维密度和磁性结构,增强了反演结果的可靠性.结果揭示以白石墩次凸为中心发育了一个"北断南超"的南西-北东向不对称型凹陷,北部沉积厚度大于南部;区内主要断裂具有高角度特征,生烃中心受断裂控制,阿吾拉勒山前凹陷浅部可能发育逆掩断裂;火成岩以侵入的基性火成岩为主,伴有部分喷出火成岩和中酸性火成岩.     

利用地脉动台阵记录反演场地浅层Vs结构

比较了影响Love波相速度频散曲线的主要因素,通过虚拟反演,说明了利用Love波相速度反演浅层剪切波速度结构的可行性,以地脉动台阵3分量观测记录提取的Love波相速度频散曲线为目标函数,利用遗传算法与单纯形法混合算法,分两种情况反演了某实际工程场地的浅层剪切波速度结构。结果表明,利用从地脉动台阵3分量观测资料中提取的Love波相速度反演工程场地浅层的剪切波速度Vs结构是一种很有前景的技术方法。