CN115437008B 一种基于地质统计学的瑞雷波频散曲线反演方法及系统 （青岛地质工程勘察院(青岛地质勘查开发局)）

一种基于地质统计学的瑞雷波频散曲线反本发明涉及一种基于地质统计学的瑞雷波震横波速度结构和位置确定每一台间待估点的用第二目标函数确定最优的第二瑞雷波频散曲2对于每一台站，获取所述台站的第一初始地震横波速度结构和第一瑞雷波频散曲线；根据所有所述台站的第一地震横波速度结构和位置确定每一台间待估点的多个第二对于每一所述台间待估点的每一所述第二初始地震横波速度结对于每一所述台间待估点，以所述台间待估点对应的所有所述第二瑞线对应的第二初始地震横波速度结构即为所述台间待估点的第二地震横波速度个台间待估点的第j个第二瑞雷波频散曲线的瑞雷所述第i个台间待估点的瑞雷波相速度的估计值的计算对于每一所述层位，根据所述层位的岩石类型确定所述层位的初始s为3相邻两次迭代中第一地震横波速度结构的差值；ΔCR为相邻两次迭代中瑞雷波相速5.根据权利要求4所述的反演方法，其特征波速度结构和位置确定每一台间待估点的多个第二初始地震横波速对于每一所述随机模拟次序，根据所有所述台站的第一地震横波速根据所有所述台站的第一地震横波速度结构和位置以及所有速度已知的所述台间待确定下一个所述台间待估点的第二初始地震横若是，则结束迭代；所有所述随机模拟次序下的每一所述台间横波速度结构组成每一所述台间待估点的多个第二初始地震横波速度若否，则返回“根据所有所述台站的第一地震横波速度结构和位置正演模块，用于对于每一所述台间待估点的每一所述第二所述第二初始地震横波速度结构为输入，利用所述频散方程计算得到第二瑞雷波频散曲最优的第二瑞雷波频散曲线对应的第二初始地震横波速度结构即为所述台间待估点的第4所述第i个台间待估点的瑞雷波相速度的估计值的计算5所述瑞雷波频散曲线为瑞雷波相速度随频率变[0008]根据所有所述台站的第一地震横波速度结构和位置确定每一台间待估点的多个[0014]正演模块，用于对于每一所述台间待估点的每一所述第二初始地震横波速度结6有所述第二瑞雷波频散曲线为输入，利用第二目标函数确定最优的第二瑞雷波频散曲线，所述最优的第二瑞雷波频散曲线对应的第二初始地震横波速度结构即为所述台间待估点根据所有台站的第一地震横波速度结构和位置确定每一台间待估点的多个第二初始地震78TL)n+(ΔVs)n+1；[0053]由于第一目标函数中含有正则化参数λ,本实施例可利用L曲线法量化正则化参数的地震横波速度结构即为S1输出的台站的第一地震横波速9[0058]将已经获得的全部地震台站地下介质的第一地震横波速度结构放置在其对应的[0060](2)对于每一随机模拟次序，根据所有台站的第一地震横波速度结构和位置拟合[0061](3)根据所有台站的第一地震横波速度结构和位置以及所有速度已知的台间待估意h对应的函数v(h)，本实施例的变差函数的类型为球状模型、指数函数模型或幂函数模波速度结构组成每一台间待估点的多个第二初始地震横波波频散曲线对应的第二初始地震横波速度结构即为所述台间待估点的第二地震横波速度PVi为第i个台间待估点的瑞雷波相速度的估计值；COVN1为距离第i个台间待估点最近的N1个节点，N1的大小可自s为第s个节点到第i个台间待估点的距离的倒数；PVs为第s的随机模拟方法可以在充分考虑空间统计相关性的前提下，增加数据的随机性和多样性，站的第一地震横波速度结构；所述瑞雷波频散曲线为瑞雷波相速度随频率变化的变化曲所有所述第二