海底水合物稳定带甲烷溶解度计算模型.pdfVIP

海底水合物稳定带甲烷溶解度计算模型 官宝聪1 雷怀彦1，2龚承林1黄磊1 1.厦门大学海洋与环境学院、近海海洋环境科学国家重点实验室，福建厦门3610052.中国科学院地质与地球物理研究所兰州 油气资源研究中心，兰州,730000 摘要： 海底甲烷水合物生成模型的建立，需要对沉积物中甲烷的溶解度进行正确测算。虽然甲烷溶解 性的相关问题已经有了较为深入的研究，但是对于海底环境下甲烷溶解度测算问题的研究才刚刚起步.本 文通过测定海水深度，海底温度和地温梯度，结合实验获得的相平衡计算式，建立了水合物稳定带(HSZ) 底部及其内部甲烷溶解度的测算方法，同时运用该方法对四个典型的甲烷水合物研究区域进行评估，求得 了各区域甲烷气溶解度与深度的相关性，研究结果表明该方法简单实用，可以作为海底甲烷水合物生成模 型中甲烷气源参数的计算方法. 关键词: 甲烷水合物；水合物稳定带； 溶解度； 建模 海底甲烷水合物生成模型的建立需要明确甲烷溶解度与水合物稳定带甲烷气之间的相 互关系，这可反映甲烷水合物在沉积物中的垂向分布特征12l，而且通过测算水合物稳定带区 域的甲烷溶解度，还可以推算甲烷水合物气的资源量13l。本文依据水，水合物和游离气三相 平衡温度Ts(P)与地温梯度的交点来确定稳定带底面，利用实验所得的相平衡计算式来测 算三相平衡条件F甲烷的溶解度，并通过修正参数后的数值模型，推算出水合物稳定带内及 其以下部分甲烷的溶解度。 2 甲烷水合物稳定带底面的确定 海洋沉积环境中。水合物稳定带之下通常会形成游离气带，从而使P波速度产生巨大的 变化，这就是海底地震剖面上显示的似海底反射面(BSR)115,161，也就是水合物稳定带的底 面。根据水合物稳定带温度和压力随深度的变化关系，可得到水，水合物和游离气的三相平 衡条件。如图1中所示，纯水【12】和盐水【13】体系中三相平衡实验数据用圆点标注，Brown和 Bangs[14J经验值用实线表示。盐水体系的三相平衡温度T3(P)(左边曲线)要比纯水体系 的三相平衡温度(右边曲线)低大约1．5℃。通过水．水合物．游离气三相平衡点温度T3(P) 与地温梯度线的交点(图中红圈标注)所得平行于海底的面来确定稳定带底面(如图l中虚 线所示)． 假设存在流体静力学关系式(1)，表示海底不同深度的压力值： P(z)=励留(日+z) (1) 式中，Z表示海底以下深度，∥表示海水的密度(1035kg／m3)，g表示重力加速度(9．8m／s2)， 71。 H表示水深【1 ，’ 薯 兰 型 器 图2水．水合物．游离气体系中的三相平衡温度(T3)和压力(P3)之间的比较 of and the the Fig．2Comparisonthree·phaseequilibriumpressure(P9temperature(T3)inwater-hydrate-freegas system． 假设海底以F温度与深度的线性关系式为： r(z)=7’(O)+Gz (2) 式中，T(O)表示海底表面温度，G表示海底地温梯度。理论上，随着深度Z的逐渐递增， 温度T(z)也随之增加。当T(z)=T3(P)时的Z值就是水合物稳定带底面所处的深度。 T3(P)可以由经验公式或者插值法求得。若T(z)<T3(P)，则不断增大Z值代入关系式 fz)和P(z)即是水合物稳定带底面 (2)计算矗到T(z)=T3(P)为止，最终得到的T 的温度和压力大小。如图2