天然气水合物沉积层渗流特性的模拟.pdf

15 7 6 化 工 进 展 C HE MI C A L I ND US T R Y A ND E NG INE E R I N G P R OG R E S S 2 0 1 5年第 3 4卷第 6期 天然气水合物沉积层渗流特性的模拟 车雯，梁海峰，孙 国庆 ，苏向东 ，吕亮国 太原理工大学化学化工学院，山西 太原 0 3 0 0 2 4 摘要天然气水合物开采过程中水合物饱和度的变化会引起储层渗透率的相应变化，对开采过程造成影响。为 研究天然气水合物对多孔介质渗流特性的影响，本文基于孔隙网络模型模拟研究了水合物生成于壁面与中心两 种方式下，多孔介质渗流特性变化，并与相关模型进行比较。结果表明，水合物生成于中心时绝对渗透率小于 生成于壁面时；水合物饱和度相同时多孔介质孔径越大，渗透率越大；水合物生成于中心时两相相对渗透率等 渗点小于生成于壁面时；当水合物饱和度变化时两相相对渗透率几乎不变。说明了储层渗透率与水合物饱和度 之 间有相对应的关系。 关键词天然气水合物；孔隙网络模型；绝对渗透率；相对渗透率 中图分类号T E 3 1 2 文献标志码A 文章编号1 0 0 06 6 1 3 2 0 1 5 0 61 5 7 60 6 DO I 1 0 ． 1 6 0 8 5 ． i s s n ． 1 0 0 0 - 6 6 1 3 ． 2 0 1 5 ． 0 6 ． 0 1 2 S i mu l a t i o n s t u d y o n t h e s e e pa g e c h a r a c t e r i s t i c s o f n a t u r a l g a s h y d r a t e s e di m e n t C HE W e n ，L I A NGHa if e n g，S U NGu o q i n g ，S UXi a n g d o n g，L UL i a n g g u o T a i y u a n Un i v e r s i t y o f T e c h n o l o g y，T a i y u a n 0 3 0 0 2 4，S h a n x i ，C h i n a Abs t r a c t T h e c h a n g e o f h y d r a t e s a t u r i o n i n t h e p r o c e s s o f n a t u r a l g a s h y d r a t e e x p l o i t a t i o n wi l l c a u s e t h e c o r r e s p o n d i n g c h a n g e o f r e s e r v o i r p e r me a b i l i t y ． An d i t wi l l a f f e c t the p r o c e s s o f n a t u r a l g a s h y d r a t e e x p l o i t a t i o n ． To s tud y the e f f e c t o f n a t u r a l g a s h y d r a t e o n p o r o u s me d i a s e e p a g e c h a r a c t e r i s t i c s ， the p o r e n e t wo r k mo d e l i s a p p l i e d t o s tud y t h e c h a n g e o f s e e p a g e c h a r a c t e r i s t i c s o f p o r o u s me d i a ． Hy d r a t e f o r ms i n t wo wa y s ，wh i c h a r e wa l l s u r f a c e a n d p o r e c e n t e r ，c o mp a r e d wi t h r e l e v a n t mo d e l s ． I t i s f o u n d t h a t t h e a b s o l u t e p e rm e a b i l i t y o f h y d r a t e f o r ma t i o n t h e p o r e c e n t e r i s l e s s tha n tha t g e n e r a t e d o n wa l l s u r f a c e ．Wh e n the h y dra t e s a t u r a t i o n i s the s a me ，the r a d i u s o f p o r o u s me d i a i s l a r g e r ，t h e n t h e p e rm e a b i l i t y i s b i g g e r ． T h e i s o t o n i c p o i n t o f t wo p h a s e r e l a t i v e p e rm e a b i l i t y o f h y dra t e f o rm a t i o n i n p o r e c e n t e r i s l e s s t h a n tha t g e n e r a t e d o n wa l l s u r f a c e ． Wh e n h y dra t e s a t u r a t i o n c h a n g e s ， tw o p h a s e r e l a t i v e p e rm e a bi l i t y i s a l mo s t c o ns t a nt ．I t a l s o g i v e s t he c o r r e s po n di ng r e l a t i o ns h i p be tw e e n r e s e r vo i r p e rm e a b i l i t y a n d h y dra t e s a t u r a t i o n ． Ke y wo r d s n a t u r a l g a s h y d r a t e ；p o r e n e tw o r k mo d e l ；a b s o l u t e p e rm e a b i l i ty；r e l a t i v e p e rm e a b i l i t y 天然气水合物作为储量丰富、燃烧清洁、能量 密度高的优质能源得到越来越多的关注， 是 2 l世纪 最具商业开发前景的清洁能源之一L 1 J 。由于天然气 水合物存在于海底沉积层及永久冻土带 ，所 以其所 处的地质环境复杂 ，开采过程 中影响因素多，以至 于其开采过程复杂 ， 需要考虑水合物藏 的地质条件、 渗流特性、化学性质及开采后的处理等问题。在水 合物开采过程 中，随着水合物饱和度的变化，会 弓 收稿 日期2 0 1 4 ． 1 1 ． 0 5 ；修改稿日期2 0 1 4 1 2 ． 1 2 。 基金项 目 山西省青年基金 2 0 1 2 0 2 1 0 0 2 2 ． 5 、国家 自然科学基 金 5 1 0 7 4 1 1 1 及国家自然科学基金青年科学基金 5 1 1 0 6 1 0 4 项目。 第一作者车雯 1 9 8 9 一 ，女，硕士研究生，主要从事天然气水合物 渗流特性模拟研究 E - ma i l s h a n x i c h e w e n 1 6 3 ．c o rn 联系人梁海峰， 博士，讲师，研究方向为多孔介质内传热传质及天然气水合物开采技 术。E ． ma i l l i a n g h a i f e n e t y u t ． e d u ．c rl 。 学兔兔 w w w .x u e t u t u .c o m 第 6期 车雯等天然气水合物沉积层渗流特性的模拟 1 5 7 7 起沉积层渗透率、 孔 隙度等储层参数的变化 。因此 ， 对于水合物的实际开采 ，进行储层渗流特性变化研 究是很有必要的。 在水合物开采过程 中，随着水合物饱和度 的变 化，储层渗透率也会随之变化。 目前国内外许多学 者 对 含 水 合 物 多 孔 介 质 渗 透 率 进 行 了研 究 。 Y o s h i h i r o等【 2 J 用一套创新的设备测量 了有 甲烷水合 物存在的沉积层有效的气一 水渗透率 。An j a n i 等【 3 J 用 实验方法研究了不同水合物饱和度下多孔介质 中气 体 的渗透率，并将实验结果与经验公式计算结果对 比得到了当初始含水饱和度小于 3 5 ％时水合物形成 于颗粒的表面， 但是 当初始含水饱和度大于 3 5 ％时， 水合物主要生成于颗粒间隙。 An d r e w等【4 J 用实验方 法研究了含不同饱和度水合物 的松散岩心样 品的渗 透特性，实验结果表 明，水合物饱和度为 1 ． 5 ％～ 3 6 ％时可以显著降低多孔介质 的渗透率 。R u a n等I S 1 用格子 Bo l t z ma n n模型对多孔介质 中水合物生成、 分解过程中饱和度变化对多孔介质渗透特性影响进 行模拟，结果表明多孔介质中单相流的流场分布与 孔隙直径和渗透率有关，沉积物 中水合物 的生成使 多孔介质渗透率大大降低 。 L i a n g等【 6 ] 用三维孔隙网 络模型研究了水合物的生成与生长习惯对渗透率的 影响，该模型表明随着水合物饱和度 的增加，渗透 率呈指数递减 。L i 等[ 7 j 用实验方法研究了多孔介质 的流动特性与水的渗透率在有水合物沉积和没有水 合物沉积的疏松多孔介质中的关系，结果表明残余 气体在水 的流动中起着 阻碍作用 ，并且会略微降低 水在多孔介质 中的渗透率，这种影响在高压条件下 表现 的尤为 明显。宋永臣等 J 开发 了一套测量渗透 性的实验装置 ，研究了多孔介质孔隙度与渗透率的 关系 ，不 同水合物饱和度下渗透率的变化情况 ，结 果表明多孔介质中甲烷水合物的存在会导致其渗透 率急剧下降，饱和度． 渗透率曲线呈指数分布 。国内 外一些学者【 9 2 J 已经开始采用孔隙网络模型研究水 合物的分布和生成 ，其 中 T s i mp a n o g i a n n i 等用二维 孔隙网络模型研究了天然气水合物分解时 甲烷气体 释放模式。但是他们并没有考虑在水合物的存在下 多孔介质中流体 的流动特性，本文将采用孔隙网络 模型来研究水合物存在下多孔介质渗透率的变化 。 1 模型建立 运用孔 隙网络模 型研究多孔介质 中的流体流 动，首先要建立二维或是三维孔隙网络结构，同时 要确定孔隙网络的空间参数。描述孔隙网络的空间 参数主要有孔喉半径、配位数、网络的维数、孔 喉长度 、孔喉 比、形状因子、孔喉分布规律 以及其 他性质参数 。孔隙和吼道常见的分布函数有均匀随 机分布 、正态分布、对数正态分布、伽玛分布、截 断正态分布和截断威布尔分布。 采用孔隙网络模型计算多孔介质的渗透率，还 需要对吼道 中流体流量进行计算 ，图 1 为孔隙 i 、 J 之间的流体流动示意图。认为节点 i 、 J之间发生层 流 ，流体的体积流量计算公式如式 1 。 _ g／j -一 1 ’ 式 中， 、 是节点 i 、 J的压力； 为节点 i 、 的孔隙中心间距 ；g u是节点 、 之间流体的导 流率 。 孔隙之间的导流率可以认为是各个导流率 的调 和平均数，其计算公式如式 2 。 土l_ L f 土 土1 2 g gt 2 gi g J 式中， 孔隙体长度 ， f 、 ， 『 分别是孔隙 、 的长度 ； 下角标 t 指连通的喉道 。 图 1 孔 隙 i 、J之间的导流率示意 图 对于孔隙 i 流量守恒 ，由此可 以得式 3 。 0 3 7 学者 P a t z e k等 。 通过求解 P o i s s o n方程得到单 相流体流动的导流率方程如式 4 。 gc G 2 o 4 式 中， 为流体黏度；G 为柱体形状 因子 ；A 为柱体横截面积；c为形状系数，正方形、三角形 和圆形界面孔喉，c值分别为 0 ． 5 6 2 3 、0 ． 6 、0 ． 5 。 在建立的孔隙网络结构中，认为孔隙是 由喉道 产生 的，从而不考虑孔隙的体积，网络结构的配位 数为 6 ，图 2为三维孔隙网络结构示意图。孔隙之 间喉道为圆柱形毛细管，本研究中，认为管径大小 服从对数正态分布，其表达式如式 5 。 e ㈤ 学兔兔 w w w .x u e t u t u .c o m 化 工 进 展 2 0 1 5 年第 3 4卷 式 中，D为喉道直径； 为平均半径； 为标准 方差。 假设孔隙之间连接的喉道长度是相等的，其长 度大小确定多孔介质 的孔隙率，孔 隙率 小可以 通过式 6 计算。 ，⋯， ㈤ 式中， 是多孔介质的总体积 ； A 是第 i 连接 喉道的横截面积；J是网络模型的总连接喉道数 。 图 2 三维孔隙网络结构示意图 2 多孔介质渗流特性 2 ． 1 绝对渗透率 多孔介质 中流体的运动方程实质是动量守恒定 律的数学表达式，即控制单元体中运动流体的动量 变化率等于所有有效的作用外力之和，如式 7 。 口 F ． 7 ‘ 其 中，F为各种质量重力；盯为各种表面力； 为是流 体质 点的运动速度 ，式 7 也就 是 Na v i e r - S t o k e s 方程。 1 9 5 6年学者 H u b b e r t [ 对 Na v i e r - S t o k e s 方程简化得到 了 D a r c y线性渗流定律。绝对渗透率 计算公式如式 8 。 K垒 8 式中，Q为流体的体积流量；A为岩心的渗流 横截面积； 为流体的黏度； △ P为岩心两端的压差； 为岩心的长度 。 采 用孔 隙 网络 模 型计 算 绝对 渗 透率 的步骤 如 下 。 整个网络完全饱和一种流体；计算网络模型中 孔隙和吼道的水力传导率；网络模型的所有节点体 积守恒；对进 出口压力进行赋值；计算每个节点的 压力；计算入 口流量和 出口流量；应用达西定律计 算绝对渗透率。 2 ． 2 两相流相对渗透率 在水合物实际开采过程 中，会存在气、水两相 共存的状态 ，这就需要考虑气、水两相的相对渗透 率，所 以本研究拟采用孔隙网络模型来计算气 、水 的相对渗透率，且计算两相相对渗透率时假设水合 物不分解和整个多孔介质 中温度是均匀分布的。其 中气相、 水相 的有效渗透率计算如式 9 、 式 1 0 。 9 K 1 0 A P l P 2 式中，q g 是气体流量；q w 是水流量 ； 是气体 黏度； 是水黏度 ；L为岩心的长度； 为岩心的 渗流横截面积； P 1 、 P 2 分别是进出口压力； P 是大 气压力。 气相、水相的相对渗透率计算方法如式 1 1 、 式 1 2 。 1 1 ’ 惫 式中， 为束缚水状态下气相渗透率； 为残 余气状态下水相渗透率。 用网络模型计算两相相对渗透率时，将气相与 水相看作一个整体， 是水相饱和度 ，是气相饱 和度 ，其两相饱和度 的关系如式 1 3 。 1 1 3 采用孔隙网络模型计算两相相对渗透率的步骤 为整个网络完全饱和一种流体 ；另一种流体由一 恒定流速进入网络；对进出口压力进行赋值；计算 网络模型中孔隙和吼道的流体传导率；模型的网络 所有节点体积守恒；计算每个节点的压力；计算入 口流量和出 口流量；应用式 9 ～式 1 2 计算气、 水两相相对渗透率。 3 模拟结果及讨论 3 ． 1 单相流模拟 研究水合物存在下多孔介质绝对渗透率，本文 选用水作为流体来进 行模拟 。网络模 型大小选取 N 1 0 ，2 2 0 ． 5 1． t in， 0 ． 5 1a m，并且假设整个多孔介 学兔兔 w w w .x u e t u t u .c o m 第 6期 车雯等天然气水合物沉积层渗流特性的模拟 质中温度是均匀分布的；水合物分别生成于孔隙中 心和孔隙壁面。 目前 已有许多模型用来估算水合物 存在 下多孔介质 的渗透率 ，本研究将网络模型与平 行毛细管模型[ 1 5 J 、Ma s u d a 模型L l 6 J 进行 了比较，图 3 为几种模型的比较。从图 3中可以看出，不论水合 物是生成于孔隙中心还是孔隙壁面，随着饱和度的 增加 ， 渗透率都是减小的； 并且对应于相 同饱和度 ， 水合物生成于中心时的渗透率小于生成于壁面时的 渗透率 ；生成 于中心时网络模 型与 Ma s u d a模型 N 1 0时比较接近 ， 生成于壁面时网络模型与平行毛 细管模型基本一致。说明应用网络模型对多孔介质 渗透率进行估算也是可行的，且还可 以考虑孔径大 小与分布。 图 3 网络模型与平行毛细管模型、Ma s u d a 模型的比较 3 ． 2 孔径大小对单相流的影响 当多孔介质有水合物存在时，其渗透率不仅与 水合物饱和度有关系，也与多孔介质的孔径大小有 关系 。图 4 、图 5分别为水合物生成于中心和壁面 时不 同饱和度下，孔径大小与渗透率 的关系。从图 4 、图 5中可 以看到，对于两种生成方式，当孔径较 小时渗透率变化幅度小 ，反之则较大；并且对应于 相 同饱和度 ，孔径大了渗透率大 。可见当饱和度相 同时，由于流体通过孔径大的多孔介质其流通面积 霸 图 4 水合物生成于中心孔径大小与渗透率的关系 图 5 水合物生成于壁面孔径大小与渗透率的关系 相对较大，故其渗透率也就随之变大。 3 ． 3 两相流模拟 水合物存在下气、 水两相相对渗透率模拟研究， 网络结构 的参数与单相时一样 ，水合物生成于孔隙 中心与孔隙壁面。 其中束缚水饱和度为 0 ． 3 ， 残余气 饱和度为 0 ． O 5 。在水合物 降压分解研究中，对相对 渗透率的计算，有时会采用修正的 B r o o k s C o r e y模 型[ 1 7 ] ，图 6为 B r o o k s ． C o r e y模型与网络模型所计算 的两相相对渗透率。从图 6中可以看到，水合物生 成于 中心的等渗点低于生成于壁面 ，这是由于生成 于中心使流体流动方式为环流，流动阻力变大 ，渗 透率变小；B r o o k s C o r e y模型的等渗点介于网络模 型之间， 表 明网络模型能更好模拟两相相对渗透率。 3 ． 4 水合物饱和度对两相流的影响 在水合物开采过程中，随着开采 的进行，水合 物饱和度会发生变化 ，致使多孔介质中流体的流动 状态发生改变。图 7 、图 8为不 同水合物饱和度下 气 、水两相相对渗透率的变化曲线。由图 8中可以 看到，对于气、水两相相对渗透率，当水合物饱和 度发生变化时，渗透率几乎不变 。这是由于当水合 物饱和度发生变化时，多孔介质 的孔隙随之变化 ， 相应的气、水两相有效渗透率和绝对渗透率也会变 - ， C o f e y 模型、 _． c o r e y 模型 ， 网络模型壁面 ， ／ e L网 络模 型壁 面 ， 网络模型中心 - ， 网络模型中心 0 0 ． 2 0 ．4 0 ． 6 0 ． 8 1 ． 0 ＆ 图6 B r o o k s - C o r e y模型与网络模型比较 O 8 6 4 2 1 O O O O 学兔兔 w w w .x u e t u t u .c o m 1 5 8 0 化 工 进 展 2 0 1 5年第 3 4卷 图 7 生成于中心两相渗透率与水合物饱和度的关系 图9 毛细管压力随水饱和度的变化 图 8 生成于壁面两相渗透率与水合物饱和度的关系 化； 相对渗透率是有效渗透率与绝对渗透率 的比值， 两者都在变化，其比值几乎不变 ，故两相相对渗透 率几乎不变。 4 毛细压力 多孔介质多相流时，气、水之 间会存在毛细管 压力。当气、水之间状态发生变化时，毛细管压力 也会随之变化。本网络模型根据 V a n Ge n u c h t e n t 】 修正模 型计 算毛细 管压力 ，计算 公式如 1 4 、 式 1 5 。 。 w - - 1 r 二 1 5 1 一S 一s 式 中 ，P c 。 l k P a ，m 0 ． 4 5 ， 0 _ 3 ， 0 ． O 5 。 毛细管压力与气、水饱和度有关系，图 9为毛 细管压力随水饱和度变化图。由图 9可以看到，在 气、水两相状态下，水饱和度变大，毛细管压力变 小；饱和度在 0 ． 3 ～0 ． 4 ，毛细管压力变化幅度较大； 饱和度大于 0 ．4时毛细管压力变化缓慢，直至降低 为 0 。这是因为随着水饱和度的增加，气体压力在 逐渐减小，气、水两相的压力差也随之变小。 5 结 论 1 两种生成方式下， 水合物生成于壁面时多 孔介质绝对渗透率较大 ；网络模型计算的渗透率介 于平行毛细管模型和 Ma s u d a模型之间。 2 多孔介质孔径越小， 渗透率变化幅度越小； 且当水合物饱和度相同时， 孔径越大 ， 渗透率越大 。 3 水合物生成于中心时两相相对渗透率的等 渗点小于生成于壁面时；网络模型与 B r o o k s C o r e y 模型相比较，B r o o k s C o r e y模型计算的等渗点介于 网络模型之间，且水合物饱和度发生变化时，两相 相对渗透率基本不变 。 4 在气、水两相状态下，毛细管压力随水饱 和度变化 ，且饱和度越小，毛细管压力越大。 5 通过储层渗透率的模拟研究， 得到了水合 物饱和度与储层渗透率之间的关系，从而为水合物 实际开采提供理论支持 。 参考文献 [ 1 ] S l o a n E D， Ko h A H ． C l a t h r a t e H y d r a t e s o f N a t u r a l G a s e s [ M] ． 3 r c l e d ． Bo c a Ra t o n，FLC RC P r e s s ，T a y l o r＆ F r a n c i s Gr o u p， 2 0 0 8 ． [ 2 】 Y o s h i h i r o Ko n n o， Y u s u k e J in， T a k a s h i Ue h i u mi ， e t a 1 ． M u n l e - p r e s s u r e t a p p e d c o r e h o l d e r c o mb i n e d wi t h X r a y c o mp u t e d T o mo g r a p h y s c a n n i n g f o r g a s - wa t e r p e r me a b i l i t y me a s e me m s o f me thane h y dra t e - b e a r i n g s e d i me n t s [ J ] ． R e v i e w o f S c i e n t ifi c I n s t r u me n t s ，2 0 1 3，8 4 6 0 6 45 0 1 ． 5 ． [ 3 】 A n j ani K u mar，B r ij Ma i n i ，B i s h n o i P R，e t a 1 ．E x p e r i me n t a l d e t e r mi n a t i o n o f p e r me a b i l i t y i n t h e p r e s e n c e o f h y d r a t e s an d i t s e ffe c t o n t h e d i s s o c i a t io n c h a r a c t e ris t i c s o f g a s h y dr a t e s in p o r o u s me d i a [ J ] ． J o u r n a l o f P e tr o l e u m S c i e n c e a n d E n g i n e e r i n g ，2 0 1 0 ，7 0 1 1 4 - 1 2 2 ． [ 4 】 An d r e w J o h n s o n ，S h ir i s h P a t i l ，A b h i j i t D and e k a r ． E x p e r i me n t a l i n v e s t i g a t i o n o f g a s wa t e r r e l ati v e p e r me a b i l i ty f o r g a s h y dra t e b e a r i n g s e d i me n t sf r o m t h eM o u n t El b e r tGa sHy d r a t e S t r a t i g r a p h i c T e s t We l l ，Al a s k a No r t h S l o p e [ J ] ． Ma r i n e a n d P e t r o l e u m G e o l o g y ， 学兔兔 w w w .x u e t u t u .c o m 第 6期 车雯等天然气水合物沉积层渗流特性的模拟 1 5 8 1 2 01 1 ， 2 8 2 4 1 9 4 2 6 ． 【 5 】 R u a n X u k e ，S o n g Y o n g c h e n ，L i a n g Ha i f e n g ，e t a 1 ． Mo d e l i n g t h e e ffe c t o f p e r me a b i l i t y o n me tha n e g a s p r o d u c t i o n f r o m h y d r a t e s i n p o r o u s m e d ia l J ] ． A p p l i e d Me e h a n s a n dMa t e r i a l s ， 2 0 1 0 ， 2 9 3 2 l 7 6 2 ． 1 7 6 7 ． 【 6 ] L i ang Ha i f e n g，S o n g Y o n g c h e n，L Y u ，e t a 1 ．S t u d y o f t h e p e r me a b i l i t y c h a r a c t e r i s t i c s o f p o r o u s me d i a wi th me tha n e h y dra t e b y p o r e n e t wo r km o d e l [ J ] ． J o u r n a l o f N a t u r a l G a s C h e m i s t r y ，2 0 1 0 ，1 9 3 2 5 5 - 2 6 0 ． [ 7 ] L i B o，L i X i a o s e n，L i G a n g ，e t a 1 ．Me asu r e me n t s o f wa t e r p e r me a b i l i ty in u n c o n s o l id a t e d p o r o u s me d i a wi m me t h ane h y dra t e f o r ma ti o n l [ J ] ． E n e r g i e s ，2 0 1 3 ，6 3 6 2 2 3 6 3 6 ． [ 8 ] 宋永臣，黄兴，刘瑜 ，等．含甲烷水合物多孔介质渗透性的实验 研究[ J ] ． 热科学与技术，2 0 1 0 ，9 1 5 1 5 7 ． [ 9 】 S i m p ano g i a uni s I N，L i c h t n e r P C ． P o r e n e two r k s t u d y o f me thane c l a t h r a t e h y dra te dis s o c i a t i o n [ J ] ． T r a n s a c t i o n s ， A m e r i c a n G e o p h y s i c a l Un i o n，2 0 0 3 1 4 6 8 4． [ 1 0 ] T s i m p a n o g i a n n i s I N，L i c h t n e r P C ． R o l e o f c ri ti c a l g as s a t u r a t i o n in me t h a n e p r o d u c tio n f r o m h y dr a t e d i s s o c i a t i o n a t the p o r e - n e t wo r k s c a l e [ C ] ／ ／ S P E A n n u a l T e c h n i c a l C o n f e r e n c e a n d E x h i b i t i o n， Ho u s t o n ， TX，2 0 0 4 ． [ 1 1 ] K eng Q，T s i mp ano g i a n n i s I N，Z h ang D，e t a 1 ． Nu me r i c a l m o d e li n g [ 1 2 】 【 1 3 ] [ 1 4 ] [ 1 5 ] [ 1 6 ] [ 1 7 】 [ 1 8 ] o f p o r e s c a l e ph e n o me n a d u r i n g C02 s e q u e s t r a ti o n i n o c e an i c s e d ime n t s [ J ] ． F u e l P r o c e s s T e c h n o l o g y ， 2 0 0 5 ， 8 4 1 4 - 1 5 1 6 4 7 - 1 6 6 5 ． T s i mp an o g i a n n i s I N ， L i c h t n e r P C． P o r e n e tw o r k s t u d y o f me t h an e h y dra t e d i s s o c i a t i o r I[ J 】 ． P h y s i c a l R e v i e wE，2 0 0 7 ，7 4 5 I - 1 3 ． P a t z e k T W ，S i l in D B． S h a p e f a c t o r an d h y dra u l i c c o n d u c t a n c e i n n o n c i r c u l a r c a p i l l a r i e s IO n e - p h a s e c r e e p i n g fl o w[ J ] ． J o u rna l o f c 0 ， f 0 a n d l n t e rf a c e S c i e n c e ，2 0 0 1 ，2 3 6 2 2 9 5 3 0 4 ． Hu b b e M K．Da r c y ’ s l a w an d t h e fie l d e q u a t i o n s o f fl o w o f u n d e r g r o u n d fl u i d s [ J ] ． T r a n s A I ME，1 9 5 6 ，2 0 7 2 2 2 2 3 9 ． Kl e i n b e r g L， Br e we r G， Ye s i n o ws k i J R De e p s e a NM R M e t h a n e h y dra t e gro wt h h a b i t in p o r o u s me dia an d i t s r e l a t i o n s h i p t o h y d r a u l i c p e r me a b i l i t y ， d e p o s i t a c c u mu l a t i o n ，a n d s u b m arine s l o p e s t a b i l i ty【 J ] ． J o u r n a l o fG e o p h y s ic s R e s e a r c h ，2 0 0 3 ，1 0 8 B 4 2 5 0 8 2 5 7 7 ． M a s u d a Y S，Na g a n a wa S， S a t o K． Nume ri c a l c a l c u l a t i o n o f gas h y d r a t e p r o d u c t i o n p e rfo r ma n c e fro m r e s e r v o i r s c o n t a i n i n g n a t u r a l g as h y d r a t e s [ C ] ／ ／ S P E As i a P a c i fi c Oi l a n d G a s C o n f erenc e ，Ku a l a Lum p u r ，M a l a y s i a ， 1 9 9 7 ． S u n X， Nan c h a r y N ， M o h an t y K K．1 - D mo d e l i n g o f h y dra t e d epr e s s u r i z a t i o n i n p o r o us me d i a [ J ] ． T r a n s p o r t in P o rou s Me d i a ， 2 0 05 ，5 8 3 31 5 ． 3 3 8 ． E r tk i n T， Ab o u - Ka s s e m J H， Ki n g G R．Basi c Ap p l i e d Re s e r