(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202111305987.4 (22)申请日 2021.11.05 (65)同一申请的已公布的文献号 申请公布号 CN 113987972 A (43)申请公布日 2022.01.28 (73)专利权人 中国石油大 学(北京) 地址 102249 北京市昌平区府学路18号 (72)发明人 胡晓东　周福建　李卓龙　丘阳　 罗英浩　易普康　李宇娇　左洁　 杨凯　梁天博　李奔　曲鸿雁　 姚二冬　王博　刘雄飞　 (74)专利代理 机构 北京三友知识产权代理有限 公司 11127 专利代理师 徐焕　童磊(51)Int.Cl. G06F 30/28(2020.01) G06F 119/14(2020.01) G06F 113/08(2020.01) G06F 113/14(2020.01) (56)对比文件 CN 103822097 A,2014.0 5.28 CN 112861740 A,2021.0 5.28 US 20191281 12 A1,2019.0 5.02 CN 113420710 A,2021.09.21 CN 111550230 A,2020.08.18 刘明亮 等.基 于小波的输油管道泄漏信号 去噪处理. 《现代电子技 术》 .2009,(第310期), 审查员 范玉霞 (54)发明名称 水击压力波波速确定方法、 装置及电子设备 (57)摘要 本说明书公开了水击压力波波速确定方法、 装置及电子设备， 其中所述方法包括： 获取填充 目标压裂液的管道中产生水击压力波的至少一 个位置处的压力值序列； 根据所述压力值序列的 频率确定分解所述压力值序列的最优小波基函 数； 采用所述最优小波基函数对 所述压力值序列 进行处理得到数据序列； 根据所述数据序列计算 水击压力波在所述目标压裂液中的波速。 该方案 较为简单， 所确定的目标压裂液中水击压力波波 速较为准确， 从而进一步确定的裂缝深度的准确 性较高。 权利要求书3页 说明书10页 附图3页 CN 113987972 B 2022.09.27 CN 113987972 B 1.一种水击 压力波波速确定方法， 其特 征在于， 包括： 获取填充目标压裂液的管道中产生水击 压力波的至少一个位置处的压力值序列； 根据所述压力值序列的频率确定分解所述压力值序列的最优小 波基函数； 采用所述 最优小波基函数对所述压力值序列进行处 理得到数据 序列； 根据所述数据 序列计算水击 压力波在所述目标压裂液中的波速； 获取填充目标压裂液的管道中产生水击压力波的至少一个位置处的压力值序列， 包 括： 获取填充目标压裂液的管道产生水击压力波的第一位置处的第一压力值序列和第二 位置处的第二压力值序列； 其中， 所述第一位置和所述第二位置沿管道延伸方向间隔预定 距离； 相应地， 采用所述最优小波基函数对所述压力值序列进行处理得到数据序列， 包括： 采 用所述最优小波基函数对所述第一压力值序列进 行处理得到第一数据序列， 采用所述最优 小波基函数对所述第二压力值序列进行处 理得到第二数据 序列； 根据所述数据序列计算水击压力波在所述目标压裂液中的波速， 包括： 根据所述第一 数据序列和所述第二数据 序列计算水击 压力波在所述目标压裂液中的波速； 根据所述第一数据序列和所述第二数据序列计算水击压力波在所述目标压裂液中的 波速， 包括： 计算所述第一数据 序列和所述第二数据 序列之间的延迟时间； 计算所述预定距离除以所述延迟时间得到的商， 将所述商作为水击压力波在所述目标 压裂液中的波速； 计算所述第一数据 序列和所述第二数据 序列之间的延迟时间， 包括： 取k为[‑N,N]的多个整数， 分别按照以下公式计算所述第一数据序列和所述第二数据 序列的相关程度： 其中， N为所述第一数据序列和所述第二数据序列中数据点的个数， P1(m)为第一数据序 列中第m个数据点， P2(m+k)为第二数据序列中第m+k个数据点， R12(k)为在第二数据序列相 对第一数据序列滞后k个数据点的情况下所述第一数据序列和所述第二数据序列的相关程 度； 从所述多个整数中确定相关程度取值 最大时的k 值； 按照以下公式计算所述第一数据 序列和所述第二数据 序列之间的延迟时间： k*为相关程度取值 最大时的k取值， f为压力传感器的采集频率， t为延迟时间。 2.根据权利要求1所述的方法， 其特征在于， 根据 所述数据序列计算水击压力波在所述 目标压裂液中的波速之后， 还 包括： 根据水击压力波在所述目标压裂液中的波速确定填充所述目标压裂液的井筒周围裂 缝的深度。 3.根据权利要求1所述的方法， 其特征在于， 根据 所述压力值序列的频率确定分解所述权　利　要　求　书 1/3 页 2 CN 113987972 B 2压力值序列的最优小 波基函数， 包括： 根据所述压力值序列确定水击频率； 采用目标小波基函数对所述压力值序列进行逐级小波分解得到低频段和高频段， 直至 分解出的低频 段最小且 包括所述水击频率后， 确定分解级数； 根据确定的分解级数采用多种小 波基函数对所述压力值序列进行分解； 根据分解结果从所述多种小波基函数中确定出分解所述压力值序列的最优小波基函 数。 4.一种水击 压力波波速确定装置， 其特 征在于， 包括： 获取模块， 用于获取填充目标压裂液的管道中产生水击压力波的至少一个位置处 的压 力值序列； 第一确定模块， 用于根据所述压力值序列的频率确定分解所述压力值序列的最优小波 基函数； 处理模块， 用于采用所述 最优小波基函数对所述压力值序列进行处 理得到数据 序列； 计算模块， 用于根据所述数据 序列计算水击 压力波在所述目标压裂液中的波速； 所述获取模块包括： 第一获取子模块， 用于获取填充目标压裂液的管道产生水击压力波的第 一位置处 的第 一压力值序列和第二位置处的第二压力值序列； 其中， 所述第一位置和所述第二位置沿管 道延伸方向间隔预定距离； 相应地， 所述处 理模块包括： 第一处理子模块， 用于采用所述最优小波基函数对所述第 一压力值序列进行处理得到 第一数据序列， 采用所述最优小波基函数对所述第二压力值序列进行 处理得到第二数据序 列； 所述计算模块包括： 第一计算子模块， 用于根据所述第一数据序列和所述第二数据序 列计算水击 压力波在所述目标压裂液中的波速； 其中， 所述第一计算子模块包括： 第二计算子模块， 用于计算所述第一数据 序列和所述第二数据 序列之间的延迟时间； 第三计算子模块， 用于计算所述预定距离除以所述延迟时间得到的商， 将所述商作为 水击压力波在所述目标压裂液中的波速； 其中， 所述第二计算子模块包括： 第四计算子模块， 用于取k为[ ‑N,N]的多个整数， 分别 按照以下公式计算所述第一数据 序列和所述第二数据序列的相关程度： 其中， N为所述第一数 据序列和所述第 二数据序列中数据点的个数， P1(m)为第一数据序列中第m个数据点， P2(m+ k)为第二数据序列中第m+k个数据点， R12(k)为在第二数据序列相对第一数据序列滞后k个 数据点的情况 下所述第一数据 序列和所述第二数据 序列的相关程度； 第七确定 子模块， 用于从所述多个整数中确定相关程度取值 最大时的k 值； 第五计算子模块， 用于按照以下公式计算所述第 一数据序列和所述第 二数据序列之间 的延迟时间： k*为相关程度取值最大时的k取值， f为压力传感器的采集频率， t为延权　利　要　求　书 2/3 页 3 CN 113987972 B 3