买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网获悉中国电波传播研究所(中国电子科技集团公司第二十二研究所)申请的专利一种实验室内阵列声波慢度计算方法获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN115407405B 。

龙图腾网通过国家知识产权局官网在2025-08-22发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202211027546.7，技术领域涉及：G01V1/40；该发明授权一种实验室内阵列声波慢度计算方法是由董兴蒙;夏济根;翟宇文;茹聪;路艳齐设计研发完成，并于2022-08-25向国家知识产权局提交的专利申请。

本一种实验室内阵列声波慢度计算方法在说明书摘要公布了：本发明公开了一种实验室内阵列声波慢度计算方法，包括如下步骤：步骤1，将采集的单极子阵列波形通过二维傅立叶变换，获得时域地层纵波信息：步骤2，利用条件模糊C均值聚类算法确定每一列波形初至起跳的位置：步骤3，应用贝叶斯信息准则确定初至的准确位置：步骤4，应用主成分分析PCA方法对初至点进行异常检测，并通过范数将异常点拟合到真实初至点；步骤5，应用慢度求解公式计算高分辨率纵波慢度。本发明所公开的方法，人为干涉少，参数设置简单，泛化能力强，同时适用于常规的数字声波仪器和阵列声波仪器的纵波慢度提取；使用条件模糊聚类和BIC信息准则的无监督机器学习方法计算效率高。

本发明授权一种实验室内阵列声波慢度计算方法在权利要求书中公布了：1.一种实验室内阵列声波慢度计算方法，其特征在于，包括如下步骤： 步骤1，将采集的单极子阵列波形通过二维傅立叶变换，将时间—空间域波形X0z,t转换至频率—波形域Xk,ω进行滤波，而后通过二维傅立叶逆变换，获得时域地层纵波信息Xz,t： Xz,t＝∫∫Xk,ω·Qk,ωeikz-ωtdkdω 上式中，t表示时间变量，z表示空间变量，k表示波数，i是虚数，ω表示角向量，Qk,ω表示滤波因子； 步骤2，利用条件模糊C均值聚类算法确定每一列波形初至起跳的位置p0： 条件模糊C均值聚类算法的目标函数为： 上式中，vg为第g个聚类中心，xhz,t为第h个样本点，ugh为第g个聚类中心第h个样本点的隶属度值，m是模糊系数，C为聚类个数，N为样本个数； 然后通过迭代更新隶属度值ugh和聚类中心vg来求解目标函数： 上式中，fh表示模糊C均值聚类算法的条件值，其定义为： 上式中，σg表示聚类g的方差，σmax表示所有聚类中σg的最大值，Mh、Eh和Rh分别是Xz,t波形中一组序列dh的绝对平均值、峰值功率谱密度和短期与长期平均比率： Eh＝max|Dh,ω|2 上式中，dh表示Xz,t波形中的第h组序列dh中的中心元素，常数w表示h附近窗口长度的一半，Dh,ω是dh的二维傅立叶变换的模，SW和LW分别是短期和长期窗口的长度，dj表示第h组序列中心元素dh构建的短期或长期窗口中的第j个元素值。 通过目标函数得到能够表征波列数据特征的隶属度值，当隶属度值大于预先设定的阈值时，其值对应的样本点被划分为波形信号类，并选择波形信号类的第一个分量为波列数据的初至p0； 步骤3，应用贝叶斯信息准则确定初至的准确位置p1： 贝叶斯信息准则BIC函数定义为： BICp0＝p0lnvar{X1,p0}+N-p0-1lnvar{Xp0+1,N}-p0lnN 上式中，X1,p0表示阵列波形X前p0个数据点组成的向量，Xp0+1,N表示阵列波形X剩余数据点组成的向量； 将BIC函数应用于p0附近的区间[p0-e，p0+e]，e表示一个常数，BIC值最小的点即为初至的准确位置p1； 步骤4，应用主成分分析PCA方法对初至点进行异常检测，并通过L1范数将异常点拟合到真实初至点； 上式中，Z1对应的是原始数据方差最小的方向，Z2对应的是原始数据方差最大的方向，px和py分别表示初至点的横坐标和纵坐标，pxo和pyo分别表示Z1和Z2相交点的横坐标和纵坐标，θ表示Z1方向与x方向的夹角；当异常数据投影在残差子空间的投影值大于预设阈值时，则判断为异常数据q； 剔除掉异常点后，使用L1范数通过直线拟合获得异常点位置处的初至数据，用于拟合直线的L1范数描述如下： 上式中，yl与xl之间为线性关系： yx＝a+bx 上式中，a和b表示变量x的系数； 步骤5，应用慢度求解公式计算高分辨率纵波慢度： 通过初至点pg获得对应的首波达到时间tg，利用慢度公式计算高分辨率纵波慢度： 上式中，s表示纵波慢度，RR表示接收器间距。

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人中国电波传播研究所(中国电子科技集团公司第二十二研究所)，其通讯地址为：266107 山东省青岛市城阳区仙山东路36号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。