申请/专利权人：青岛海洋地质研究所

申请日：2023-12-14

公开（公告）日：2024-03-12

公开（公告）号：CN117408186B

主分类号：G06F30/28

分类号：G06F30/28;G06F30/23;G06F113/08;G06F111/10;G06F119/14

优先权：

专利状态码：有效-授权

法律状态：2024.03.12#授权;2024.02.02#实质审查的生效;2024.01.16#公开

摘要：本发明所述的地震剖面中天然气水合物似海底反射层的判断方法，属于海洋地球物理领域，该方法提出在数值模型中求解二维剖面中天然气水合物似海底反射层位置的技术手段，以期通过两矩阵相减结果的临界数值判断依据提高BSR存在解释结果的准确率和可信度，同时达到缩短计算BSR二维深度步骤与时长的设计目的。包括下述步骤：步骤1、确定模型输入参数；步骤2、建立数值模型；步骤3、设置边界条件；步骤4、计算沉积物节点温度；步骤5、计算第一压力矩阵M1；步骤6、计算第二压力矩阵M2；步骤7、获得单个BSR位置；步骤8、获得多个BSR位置；步骤9、对比实际资料进标记。

主权项：1.一种地震剖面中天然气水合物似海底反射层的判断方法，其特征在于：包括下述步骤，步骤1、确定模型输入参数；依据目标二维剖面情景，沿横向和纵向分别建立x轴和z轴，确定数值模型的总长度L和总高度H，获取当前海底x和z坐标位置；获取当前海底实测温度、沉积物地温梯度、海水盐度和密度；步骤2、建立数值模型；根据模型尺寸确定x轴和z轴方向上的节点数量；采用二维热传导稳态方程建立数值模型， 其中，T为温度，kx和kz分别是沉积物在水平和垂直方向上的热导率；步骤3、设置边界条件；将步骤2中的数值模型赋予边界条件：模型左边界海底下方沉积物节点温度为海底温度加上地温梯度增值；模型右边界海底下方沉积物节点温度为海底温度加上地温梯度增值；模型下边界的温度梯度为固定值；海底温度为初始赋值，数值来自步骤1提供的实测值；步骤4、计算沉积物节点温度；将主控方程使用有限差分法进行离散化处理，求解方程；在求解过程中调整矩阵参数排列位置，使之成为对角线元素矩阵，之后使用spdiags函数命令和MLDIVIDE方程进行求解，获得稳态条件下的沉积物节点温度矩阵T1；步骤5、计算第一压力矩阵M1；根据实验室内水合物相平衡方程，将沉积物节点温度矩阵T1中的各数值转化成压力值，按照节点位置形成二维第一压力矩阵M1，数值单位为兆帕；步骤6、计算第二压力矩阵M2；按下述公式计算海底下方沉积物内各节点的静水压力： 其中，Phy是静水压力，h1和h2分别是具体节点上方的海水深度和沉积物厚度，ρswz是海水密度，g是重力加速度；计算后按照节点位置形成二维第二压力矩阵M2，结果单位为兆帕；步骤7、获得单个BSR位置；将第一压力矩阵M1和第二压力矩阵M2计算相减，将差值0.001的位置标记为数值模型节点中天然气水合物底界面的位置；将上述标记的位置在目标二维剖面中连成一条连续的曲线，将符合标记条件的位置坐标x0，y0导出；步骤8、获得多个BSR位置；设置地温梯度数值范围，赋予多个数值，重复以上过程，获得一组BSR位置坐标x1，y1；x2，y2，……，xn，yn，上述位置坐标范围用以约束BSR出现的空间；步骤9、对比实际资料进行结果标记；将上述步骤7和步骤8计算得出的坐标x，y数据导出至二维剖面中，沿导出的坐标数据形成的曲线范围寻找BSR相关特征的反射层；标出符合上述条件的反射层，该反射层即为天然气水合物似海底反射层BSR，流程结束。

全文数据：

权利要求：

百度查询： 青岛海洋地质研究所 地震剖面中天然气水合物似海底反射层的判断方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD