非精确搜索一维全局最优化方法

研究一维全局最优化问题的确定性求解方法。运用逐次建立目标函数的线性下界函数,将不含全局最优解的子区域删除,并基于非精确搜索结合下降算法而得出非精确搜索一维全局最优化方法,使计算量减少且使迭代收敛加快。迭代结束时该算法得到一维全局最优化问题的ε－全局最优解。该方法具有有限收敛性且不需精确的局部优化过程。文中的数值实例表明该算法的有效性。     

基于改进LM-GN算法的磁性目标定位方法研究

利用三轴磁传感器阵列对水下磁性目标进行定位是典型的非线性最小二乘优化问题,传统高斯– 牛顿（GN）法和列文伯格–马夸尔特（LM）算法在求解该问题时具有初值敏感性问题。通过引入信赖域搜索技术,对 LM 算法进行改进,并基于改进的 LM 算法实现磁性目标定位,又通过设置判断阈值来评估迭代点与最优解的距离,提出一种结合改进 LM 算法和高斯–牛顿法特点的改进 LM-GN 算法,既降低了算法对初始值的依赖性,又提高了运行效率。仿真实验结果表明：该方法可以克服现有方法中受初始值影响较大的问题,对目标特征参数的估计更精确,且收敛速度快,具有一定的实际应用价值。     

基于数据分布域变换与贝叶斯神经网络的渗透率预测及不确定性估计

渗透率是储层评价和油气藏开发的关键参数.传统测井方法与常规机器学习方法估算的渗透率都是固定值.但由于测井数据本身存在噪声,渗透率的预测结果可能受到噪声的影响出现测量性的随机误差(即任意不确定性);同时,当测试数据与训练数据存在差异时,机器学习模型在预测渗透率时可能出现模型参数的不确定性(即认知不确定性).为实现渗透率的准确预测并量化两种不确定性对结果的影响,本文提出基于数据分布域变换和贝叶斯神经网络同时实现渗透率预测及其不确定性的估计.提出方法主要包括两个部分：一部分是不同域数据分布的相互转换,另一部分是基于贝叶斯理论的神经网络渗透率建模预测和不确定性估计.由于贝叶斯神经网络存在数据分布的假设,当标签的概率分布与网络的分布保持一致时,贝叶斯神经网络可以更好的学习到数据之间的关系.因此通过寻找一个函数将一个原始域的渗透率标签转换为目标域的与渗透率有关的变量(我们称为目标域渗透率),使得该变量符合贝叶斯神经网络的分布假设.我们使用贝叶斯神经网络预测目标域渗透率以及任意不确定性和认知不确定性.随后,通过分布域的逆变换,我们将目标域渗透率还原回原始域渗透率.应用本文方法到某油田的18口井的测井...     

基于经验格林函数法的2021年MS6.4漾濞地震近场强震动模拟

由于震中附近强震台分布不均,文中联合使用近场2个国家强震动观测台网的强震台和7个云南地震预警台网烈度台的强震记录,采用经验格林函数法构建了特性化震源模型,并利用此模型对近场强震动进行了模拟。结果表明：烈度台的记录可与强震台的记录联合,共同作为强震动模拟的对象,但应注意区分二者的有效频带;在0.20～30.0Hz频带,在53YBX台处的NS分量上,模拟结果的伪加速度反应谱较好地再现了0.1s处的峰值,在EW分量上,合成的速度波形虽然幅值较低,但较好地再现了速度脉冲波段;在53DLY台处的合成波形较好地再现了约2s的长周期地震动;在0.50～30.0Hz频带,合成波形和反应谱与信噪比较高的烈度台的记录较为一致。文中确定的用于模拟强震动的特性化震源模型由一个强震动生成域构成,其面积和相应的短周期范围内加速度震源谱的水平段幅值与地震矩的关系均遵循经验标度律。     

声波方程频率域高精度正演的17点格式及数值实现

频率域正演计算是频率域全波形反演的基础.传统的最优9点格式只具有二阶精度,不能满足高精度地震成像的需要.本文考虑两个四阶精度的格式,即经典的四阶9点格式和优化的17点格式.17点格式可将最小波长内所需网格点数减小到2.56.通过在简单模型和Overthrust模型上的数值实验,比较分析了三种格式的正演效果;简单模型数值实验显示了17点格式克服频散误差的能力优于四阶9点格式和最优9点格式;复杂模型数值实验则进一步承认了算法的可行性.     

香溪河流域白家堡滑坡变形监测初步分析

白家堡滑坡受三峡大坝蓄水及强降雨影响变形较为明显,且仍在变形过程中。利用GDM600S型全站仪及cx-03D型钻孔测斜仪分别对滑坡地表及滑坡深部位移进行监测。监测结果表明:白家堡滑坡体呈"D"型整体运动;滑坡后缘滑动面(带)和中部滑动面(带)分别在地下11m和28.5m左右;滑坡整体滑动方向大致为NE55°,且滑坡变形具有自后缘向前缘递减的趋势。     

基于广义旋转法双九点格式标量波方程数值模拟

前人虽然基于传统旋转法提出了四阶精度最优17点差分格式,用于提高频率域地震波场数值模拟精度,但其仅适用于等网格间距的情形.这大大限制了该格式的使用范围.为了进一步提高17点有限差分格式的数值精度并将其推广到矩形网格,本文基于二阶精度有限差分算子利用广义旋转法提出了双九点格式,由于其差分格点与前人17点格式在分布上一样,所以也可称为二阶精度最优17点格式,但由差分格式构造原理上来讲,称其为双九点格式更妥.前人四阶精度17点格式仅为本文等网格间距情形时的特殊情况,本文方法单位波长网格点数仅需要2.2个即可.本文格式在继承传统旋转法良好几何旋转性质的同时,拥有平均导数方法适用于矩形网格的特点,和平均导数法所得的广义17点格式相比,本文格式数值精度更高,数值频散抑制性能和差分算子对称性更好.同时,本文双九点格式方法和思想对于后人借助传统九点格式的构造方法将其扩展到17点格式求解各类波动方程具有十分重要的意义.     

多通道瞬变电磁法频率域比值算法及视电阻率计算研究

视电阻率在电磁法勘探中起重要作用,可以作为初步解释结果,多通道瞬变电磁法通过大地脉冲响应的峰值时刻定义视电阻率.鉴于传统方法计算峰值时刻视电阻率需要恢复完整大地脉冲响应,本文通过研究大地脉冲响应及其频谱特征,提出一种新的频率域比值法.利用参考谱及其与一般谱的伸缩因子计算一般谱的峰值时刻,进而得到视电阻率.该方法只需大地脉冲响应单个频点幅度谱,从而不用估计整个大地脉冲响应.实际计算时,可以用大地脉冲响应多个频点幅度谱得到平均视电阻率.模拟数据与实测数据结果表明,在一维介质中,频率域比值法可以得到合理的视电阻率.     

GA辅助NLS的GNSS-IR土壤湿度反演方法北大核心CSCD

采用具有阻尼因子的函数模型,使用遗传算法(genetic algorithm,GA)辅助非线性最小二乘(nonlinear least squares,NLS)方法对相位参数进行求解。结果表明:1)相较于标准余弦函数模型,该方法的反演相位与土壤湿度的相关系数有较为明显的提升,反演结果也更加稳定,在5°~15°、5°~20°、5°~25°三个高度角范围内的相关系数均大于0.68,不同高度角之间的相关系数差值小于0.07;2)反演精度有不同程度提高,R 2提高5.72%~76.06%,RMSE减小6.12%~24.24%,MAE减小2.7%~28.3%,将该方案所求相位用于多星线性回归模型后平均RMSE减小10%。     

天津港“8·12”爆炸事件记录波形分析

利用天津地震台网23个短周期地震仪记录到的2015年8月12日天津港爆炸事件波形数据,对事件进行初步分析和研究。分别使用能量包络线法和波形匹配法,共识别出4次爆炸事件,并结合波形互相关的主事件定位法对爆炸事件进行定位。分析结果表明:第一次爆炸发生的时间为2015年8月12日23时34分4秒,1.77s后发生了规模更大的二号爆炸事件;34.11s后又发生了最强的三号爆炸事件;86.33s再一次发生了能量较弱的四号爆炸事件,其中二号与三号爆炸事件相距65~70m,二号与四号事件相距160~170m。此外,通过选用爆破、矿塌、核爆和天然地震等不同类型事件的波形数据,从时间域和频率域分别对这些事件进行对比分析,发现其震相特征表现出明显差异。