(19)中华 人民共和国 国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202110934284.1 (22)申请日 2021.08.13 (65)同一申请的已公布的文献号 申请公布号 CN 113552630 A (43)申请公布日 2021.10.26 (73)专利权人 广州海洋地质调查局 地址 510000 广东省广州市南沙区环市大 道南25号南沙科技创新中心A4栋2楼 208-218房 专利权人 南方海洋科学与工程广东省实验 室 （广州） (72)发明人 邓炜　梁金强　匡增桂　钟桐　 何玉林　张伟　谢莹峰　孟苗苗　 (74)专利代理 机构 广州君咨知识产权代理有限 公司 44437 代理人 张华华 (51)Int.Cl. G01V 1/30(2006.01) G01V 1/40(2006.01)G01N 15/08(2006.01) G06F 30/28(2020.01) G06F 30/27(2020.01) G06N 3/08(2006.01) G06F 111/10(2020.01) G06F 113/08(2020.01) G06F 119/14(2020.01) (56)对比文件 CN 111859632 A,2020.10.3 0 CN 109882134 A,2019.0 6.14 CN 111601946 A,2020.08.28 CN 112946783 A,2021.0 6.11 CN 112213780 A,2021.01.12 US 8301383 B2,2012.10.3 0 US 201417217 7 A1,2014.0 6.19 WO 2008066389 A1,20 08.06.05 邓炜等.高渗水合物储层的孔隙可压缩性反 演. 《地球物理学报》 .2021,第64卷(第1期),第 298-313页. 审查员 王叶 (54)发明名称 基于弹性 阻抗的未固结地层渗透率预测方 法及处理终端 (57)摘要 本发明涉及一种基于弹性 阻抗的未固结地 层渗透率预测方法及处理终端， 所述方法包括： 计算出水合物储层的流体模量下限和流体模量 上限； 根据流体模量下限和流体模量上限计算水 合物储层的流体模量； 根据水合物储层的流体模 量计算出水合物储层的双状态模量差， 通过中子 孔隙度反演计算得到中子孔隙度； 将双状态模量 差和中子孔隙度作为输入参数， 输入 预测神经网 络进行处理， 预测神经网络的输出即为水合物储 层的渗透率。 本发明能够很好地对作为未固结地 层的水合物储层的渗透率进行预测， 且能够进行 定量计算， 预测结果稳定可靠， 实用性强。 权利要求书2页 说明书7页 附图3页 CN 113552630 B 2022.03.04 CN 113552630 B 1.一种基于弹性阻抗的未固结地层渗透率预测方法， 其特在于， 包括如下步骤： 步骤1： 按公式 ①计算出水合物储层的流体模量 式中， 公式计算得到的Kf作为 其中， EI表示水合物储层的弹性阻抗， EI0表示入射角度 θ为 0时的弹性阻抗， θ表示入射角度， Kf0表示Kf的平均值， φ表示中子孔隙度， φ0表示φ的平均 值， fm表示固体刚度系数， fm0表示fm的平均值， ρ 表示密度， ρ0表示ρ 的平均值， vp表示纵波速 度， vs表示横波速度， 基于碎屑岩饱和岩石的固解液解耦方程计算出 水合物储层的流体模量上限 步骤2： 按公式 ②计算水合物储层的流体模量K ′f： 式中， C1和C2均为正数的权 重系数， 且C1+C2＝1； 步骤3： 按公式 ③计算出水合物储层的双状态模量差PCI: 其中， c为常数， μ表示剪切模 量， φc表示临界孔隙度， 双状态模 量差PCI表征水合物储层 为具有可压缩性且处于未固结状态的岩石， 通过中子孔隙度反演计算得到中子孔隙度； 步骤4： 将双状态模量差PCI和中子孔隙度作为输入参数， 输入预测神经网络进行处理， 预测神经网络的输出即为水合物储层的渗透率。 2.根据权利要求1所述的基于弹性阻抗的未固结地层渗透率预测方法， 其特在于， 水合 物储层的弹性阻抗EI 通过AVO地震反射特 征方程推演得到， A VO地震反射特 征方程如下：权　利　要　求　书 1/2 页 2 CN 113552630 B 2式中， ΔKf表示Kf的改变量， Δφ表示φ的改变量， Δ fm表示fm的改变量， Δρ 表示ρ 的改 变量。 3.一种处 理终端， 其特 征在于， 其包括： 存储器， 用于存 储程序指令； 处理器， 用于运行所述程序指令， 以执行如权利要求1或2所述的基于弹性阻抗的未固 结地层渗透率预测方法的步骤。权　利　要　求　书 2/2 页 3 CN 113552630 B 3