[发明专利]盐底劈构造带水合物稳定底界的二维数值模拟方法及系统

权利要求书

1.一种盐底劈构造带水合物稳定底界的二维数值模拟方法，其特征在于，所述盐底劈构造带水合物稳定底界的二维数值模拟方法包括：

综合分析地温梯度、盐和沉积物不同的热导率以及盐底劈构造带的盐度高异常对水合物稳定底界的影响，通过选取合适的水合物相平衡公式，对BSR进行数值模拟；通过设定不同的地温梯度、热导率、盐度参数值模拟不同参数条件下的BSR位置；再将模拟BSR位置与观察到的BSR位置进行拟合程度的比较；

所述盐底劈构造带水合物稳定底界的二维数值模拟方法具体包括：

步骤一：在盐底劈构造带，选择水合物发育且地震资料上BSR特征明显的地震剖面；

步骤二：确定地震剖面上海底界面的双程走时T_sb、观测BSR的双程走时T_bsr和盐底劈构造带顶部边界的双程走时T_salt；

步骤三：根据分析区天然气水合物的气体组分信息，选择合适的水合物相平衡曲线；纯甲烷水合物的相平衡公式如下：

In(P)＝a+bT+cT²+dT³+fT⁴+gT⁵

式中，P和T分别为甲烷水合物的稳定压力条件和稳定温度条件，a、b、c、d、f、g为经验常数，分别为a＝-1.94138504464560×10⁵,b＝3.31018213397926×10³,c＝-2.25540264493806×10¹,d＝7.67559117787059×10^-2,f＝-1.30465829788791×10^-4,g＝8.8606531668757×10^-8；

步骤四：地震剖面上BSR发育深度处的压力值为静水压力，压力通过以下公式求取：

P_bsr＝ρ_sw g H_bsr

式中，ρ_sw为海水密度，为1028kg/m³，g是重力加速度，为9.81m/s²，H_bsr是BSR的深度值；

步骤五：BSR的深度值通过下面公式求得：

H_bsr＝V_sw×T_bsr/2

式中，V_sw为海水速度，为1500m/s；T_bsr为BSR位置的双程走时，单位是s，通过地震剖面读取；

步骤六：通过数据库查询，确定地震剖面上海底各点的温度值T；

步骤七：通过二维稳态热传导方程，计算地震剖面的温度场分布：

式中，T为温度℃，x为横向距离km，z为垂向距离km，k_x为横向导热率W m^-1K^-1，k_z为垂向导热率W m^-1K^-1；一般认为沉积物各向同性且均质性较强，因此k_x＝k_z，但沉积物和盐的热导率差别较大，分别设置为2.5W m^-1K^-1和5.9W m^-1K^-1；

步骤八：将水合物稳定底界移向低的温度条件，

式中，m为盐度对水合物稳定条件的影响参数，计算公式如下：

式中，S_w为盐度值；

步骤九：通过设置不同的地温梯度，模拟对应的BSR深度，将模拟的BSR位置与地震剖面上的观测BSR进行对比。

2.一种实现权利要求1所述盐底劈构造带水合物稳定底界的二维数值模拟方法的信息数据处理终端。

3.一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如权利要求1所述盐底劈构造带水合物稳定底界的二维数值模拟方法。