1、南海北部神狐海域天然气水合物成藏动力学模拟第32卷 第2期2011年3月石油学报A CT A PETROLEI SINICAV o l. 32N o. 2M ar. 2011文章编号:0253O 2697(2011 02O 0226O 08南海北部神狐海域天然气水合物成藏动力学模拟苏丕波1 梁金强2 沙志彬2 付少英2 雷怀彦1, 3, 4 龚跃华2(11厦门大学海洋与环境学院 福建厦门 361005; 21广州海洋地质调查局 广东广州 510760;31厦门大学近海海洋环境国家重点实验室 福建厦门 361005; 41中国科学院地质与地球物理研究所兰州油气中心 甘肃兰州 730000摘要:为
2、进一步了解南海北部神狐海域天然气水合物的成藏匹配条件, 利用典型二维地震剖面, 构建了该区的地质模型, 并对其进行了天然气水合物成藏动力学的模拟。研究结果表明:神狐海域具备天然气水合物成藏的温度、压力条件; 生物气和热解气的资源潜力巨大, 满足水合物形成的气源条件; 运移条件优越, 有利于天然气水合物的聚集成藏。并提出了该区天然气水合物的成藏模式。关键词:南海; 神狐海域; 天然气水合物; 成藏模式; 生物气; 热解气中图分类号:T E 112 文献标识码:ADynamic simulation of gas hydrate reservoirs in the Shenhu area,the
3、northern South China SeaSU Pibo LIANG Jinqiang SH A Zhibin FU Shaoying LEI H uaiyan12221, 3, 4GON G Yuehua2(1. Colleg e of Oceanogr ap hy &Envir onmental S cience, X iamen Univer sity , X iamen 361005, China; 2. Guangz ho u M ar ine Geological S ur v ey , Guangz hou 510760, China; 3. S tate K ey L a
4、bor ator y of M ar ine Env ir onmental Science, X iamen Univ ers ity , X iamen 361005, China; 4. L anz hou Resear ch Center of Petr oleumRes our ces , I nstitute of Geology &Geop hy s ics , Chinese A cademy o f Sciences , L anz hou 730000, ChinaAbstract :T he no rth slo pe of the South China Sea is
5、str uctur ally char act er ized by a passive co nt inental mar gin and has var ious Cenozo -ic o i-l and gas -bearing deposit ional basins, w her e hy dr ocarbo n r eso ur ces ar e abundant, and the g eolog ical and tect onic settings, the temper atur e and pressur e reg imes as w ell as the methane
6、 -generat ive po tential o f t hick or ganic -r ich sediment s ar e all favo rable fo r the format ion of g as hy dr ate. A s a result, this reg ion is also an ideal area for study ing gas hydrate accumulations. I n or der to under stand natura l conditions of the g as hydrate format ion, based on t
7、he ty pical 2-dimensio nal seismic image of the Shenhu area, a geo log ical model of g as hy dr ate reserv oir s w as est ablished and a dynamic simulatio n o f g as hydrate accumulations w as car ried out as well. T he results indicated that the temper ature and pressur e of t he Shenhu ar ea are a
8、ppro priate fo r gas hydr ate accumulations, source rocks ther e have a hug e gas -g enerat ive potential and hy dr ocarbon mig ratio n conditio ns are favo rable, pro viding an abundant gas source and go od pathw ays for the g as hydrate for matio n. A forecasting mo del of t he g as hydrate fo rma
9、tion was g iv en based on basin analysis, indicating that the Shenhu ar ea is a better pr ospect a rea of g as hy dr ate because o f its hug e potential o f gas hydrat e resources. Key words :So uth China Sea; Shenhu a rea; g as hy dr ate; r eser vo ir fo rming model; biogas; thermog enic gas天然气水合物是
10、在低温、高压环境下由水和天然气组成的类冰结晶化合物, 主要赋存在陆地永久冻土带和水深超过300m 的海洋沉积物中。目前发现的海底天然气水合物主要分布于世界各大洋边缘海域的大陆斜坡、陆隆海台和盆地, 以及一些内陆海区的大洋沉积物中, 水深一般为3004000m, 赋存沉积物一般分布于海底以下01500m 深度1。控制海洋天然气水合物成藏的关键因素包括温度、压力、气体组分、气体饱和度及孔隙水组成, 水合物的结晶和生长还取决于沉积物颗粒大小、形状和组成, 但是这些因素受到海洋中一系列构造和沉积作用的影响, 在不同的时间尺度上可能导致多种天然气水合物成藏的动力学反映3-52。目前, 国内外对天然气水合
11、物赋存及分布主控因素的研究仍局限于对影响水合物成藏的个别因素探讨上, 如全球气温变化、构造活动与地热史、沉积作用效应、地温梯度和冰川性海平面相对移位等6, 这些因素均可改变天然气水合物形成所需要的温压条件与沉积物的物性特征, 从而影响基金项目:国土资源部公益性行业科研专项(200811014 /天然气水合物成藏的控制因素研究0、国家自然科学基金项目(No. 40930845 /南海北部深水盆地油气渗漏系统与天然气水合物富集机制研究0和国家重点基础研究发展规划(973 项目(2009C B 219501 /南海北部天然气水合物成藏的气源条件研究0联合资助。第一作者及通讯作者:苏丕波, 男, 19
12、81年3月生, 2004年毕业于长江大学, 现为厦门大学在读博士研究生, 主要从事天然气水合物的气源条件与成藏模拟研究。E -mail:spb _525s ina. com天然气水合物系统的稳定性。除温压条件外, 充足的气体供应是天然气水合物形成的另一个重要控制因素, 从动态过程来考虑, 除了烃类气体的供应外, 还涉及烃类气体到达天然气水合物稳定带的运移通道、天然气水合物形成的构造环境等。南海北部陆坡含油气盆地发育, 气源丰富、类型众多, 深部热解气、浅层生物气均有可能形成天然气水合物, 虽然部分学者分别就烃类气体供应问题、烃类运移条件、岩层和构造对天然气水合物产状与分布影响或控制做过单方面的
13、研究7-9, 但还没有将它们作为一个有机整体在时空尺度上开展水合物的成藏系统研究。笔者选取南海北部神狐海域的典型地震剖面, 通过水合物成藏动力学模拟, 对其天然气水合物成藏模式进行了初步探讨。条件, 而褶皱构造易于捕获天然气, 促使水合物的形成。同时, 神狐海域海底滑塌作用非常强烈, 有分析认为可能与水合物的形成和分解有关15。此外, 根据沉积相分析, 南海北部陆坡自晚渐新世以来处于坳陷沉降期, 以滨、浅海 半深海沉积环境为主, 陆源碎屑供给充足, 沉积速率大, 厚度大, 粒度总体上中等偏细。其中, 晚中新世以来神狐海域研究区以三角洲、扇三角洲、滑塌扇、浊积扇沉积为主, 重力流非常发育, 特别
14、是第四纪, 广泛发育滑塌沉积。这些沉积体普遍具有较高的沉积速率, 沉积厚度相对较大, 含有大量的有机质, 并能得以有效地保存, 能为天然气水合物的形成提供充足的气源。总之, 神狐海域研究区具备良好的天然气水合物成藏地质条件。16-17O2 模型选择及参数选取由于神狐海域探井缺乏, 因此, 根据既要在神狐海域水合物研究区具有代表性, 还要兼顾该区及邻区是否有可以借鉴的模拟参数资料的原则, 选取神狐海域水合物研究区的二维地震测线line A 进行模拟。该测线处水深为4001700m, 地层自下而上发育有始新统文昌组, 渐新统恩平组、珠海组, 中新统珠江组、韩江组、粤海组, 上新统万山组和第四系,
15、并且在水深5001000m 的万山组内识别出了指示水合物存在的BSR(Bottom Simulating Reflection, 似海底反射 特征标志(图2 。1 研究区地质概况神狐海域水合物研究区地理上位于南海北部陆缘的陆坡区中段的神狐暗沙东南海域, 即西沙海槽与东沙全岛之间海域, 构造上位于珠江口盆地珠二坳陷白云凹陷(图1 。白云凹陷水深2002000m, 面积约为20000km 2, 新生代最大沉积厚度约为12000m, 地史上经历多次地壳运动和多阶段的构造演化, 地质构造复杂, 断裂-褶皱体系非常发育10-13 。图1 研究区位置及范围Fig. 1 Location and range
16、 of the studied area神狐海域研究区晚期断裂极其发育14, 新近纪断层主要形成于晚中新世和上新世以来两个时期。晚中新世断层以NW 向为主, 断层大部分切割上中新统, 部分切割上新统, 是研究区最主要的断层; 上新世以来活动断层以NEE 向为主, 断层活动较弱, 但数目众多, 部分断层切穿较新的沉积层延伸至海底附近。深部断层为天然气向浅部水合物稳定带运移创造了有利图2 神狐海域line A 地震解释剖面及模拟地质模型Fig. 2 Seismic profile and its geological interpretationof line A in Shenhu areaO
17、于兴河, 苏新, 陈芳, 等. 南海天然气水合物成矿的沉积条件初步研究. 中国地质大学(北京 、广州海洋地质调查局, 2002.采用IES 软件中的PetroM od 2D 模块, 主要对研究区新生界地层的温压场、有机质热演化指数R o 和流体运移进行了模拟。地层压力的演化基于2个假设, 应用有限元模拟方法来模拟:首先假设岩石和孔隙流体在压缩和变形过程中保持质量平衡; 其次压实过程中, 流体排出极其缓慢, 能够以达西定律来描述。热史恢复则采用地球热力学和地球化学相结合方法, 即将正演技术与反演技术、地史恢复与热史恢复结合起来, 利用已知的地层信息和古温标资料作为约束条件, 对研究区的热演化史进
18、行模拟。有机质成熟度的计算采用Sw eeney 和Burnham 提出的EA S -Y %R o 模型18-19, 它是目前用于成熟度计算最为完善的一种模型, 不仅考虑了众多一级平行化学反应及其相应反应的活化能, 而且还考虑了加热速率, 适用范围广, 能比较精确地模拟地质过程中有机质成熟度演化。模拟中主要需要岩石性质、地质界面、烃源岩地球化学和断层活动性等参数。对这些模拟参数的选取, 综合借鉴了研究区各方面的研究成果。岩性参数来源于中海油钻探资料20。地质界面参数中, 古水深来源于高红芳等21的研究结果; 热流值来源于ODP 184航次调查成果22-23; 古地温由IES 系统根据剖面所在的全
19、球位置, 利用全球平均地表温度窗口以及古水深变化来计算。综合目前研究资料及地质分析, 认为该区主要烃源岩层为文昌组和恩平组, 对于其烃源岩地球化学参数, 采用傅宁等24的研究成果, 其中恩平组TOC 平均值为2119%, H I 平均值为15714mg/g; 文昌组T OC 平均值为2194%, H I 平均值为48314mg/g 。而断层活动性的分析主要是基于断层在地震剖面上断穿的层位以及研究区构造活动的时间来判断和估算, 本次模拟研究中, 根据断层活动期次将其划分为始新世及之前形成的活动断层、晚中新世形成的活动断层以及上新世以来的活动断层, 对剖面经过的每一条断层均进行了属性定义; 在模拟
20、过程中, 各断层活动性自构造活动时间开始均设为完全开启状态。比(图3 可以看出, 测线点模拟曲线与井测试值趋势比较一致, 说明模拟结果比较准确, 可以用于相关 解释。图3 研究区井B 地温、R o 实测值与模拟值对比Fig. 3 Relations of modeling temperature and R o trend withmeasured values in Well B311 温压场模拟天然气水合物的形成与成藏需要特定的温压条件, 低温和高压有利于水合物的形成和稳定赋存。通过模拟得到该区现今的温度场(图4 与压力场(图5 , 在剖面上BSR 所处温度约为16e , 压力约为15M
21、Pa 。对比世界上已知天然气水合物区, 结合甲烷在海水中形成水合物的相平衡曲线表明27, 该测线剖面BSR 区域处于天然气水合物稳定存在的温压场范围内, 符合天然气水合物的成藏要求。312 有机质成熟度模拟对神狐海域地质调查站位资料28O 的分析表明:研究区浅表层沉积物中普遍存在游离气, 其甲烷碳同位素D C 1值为-4612j -7413j (PDB , 平均为-6019j (PDB , 除2个样品的D C 1值为-4612j 和-51j (PDB 外, 大多数样品的D C 1值小于-57j (PDB , 证实神狐海域浅表层沉积物顶空气主要来源于生物气。同时, 许多调查站位顶空气甲烷的含量在
22、垂向上保持了相对较高的丰度, 特别是在调查区北部白云凹陷内的两个调查站位, 其顶空气甲烷的含量分别约为120L L/kg 和200L L/kg , 暗示其深部可能有持续稳定的游离甲烷供应(来源于深部的热解气 。王建桥等29131313263 模拟结果分析模拟结果是否可靠需要通过其与钻井实测值的对比来进行检验。研究区番禺低隆起有部分探井, 其中井B 有实测的地温和镜质体反射率25, 且该井与测线line A 较近, 两者的演化环境与受热历史相差不大, 因此, 可以利用该井的实测值对模拟结果进行检验。从距该井最近的剖面点模拟结果与实际井资料的对对研究区东部ODP 1146站位顶空气样品的分析结果显
23、示为混合气体的特征。由此推测, 研究区浅部地层中的天然气可能兼有生物气和热解气两种来源。镜质体反射率R o 值是反映烃源岩成熟度的重要指标。通常, 生物气的烃源岩应处于未熟 低成熟阶O 郭依群, 梁劲, 龚跃华, 等. 南海北部神狐海区天然气水合物资源概查报告. 广州海洋地质调查局, 2004. 图4 神狐海域line A 现今温度场模拟F ig. 4 Temperature field sim ulation of line A in Shenhu area图5 神狐海域line A 现今压力场模拟Fig. 5 Pressure field simulation of line A in S
24、henhu area段, 其R o 2%, 最大超过3%, 处于过成熟生干气阶段, 已产生大量热解气; 恩平组R o 为113%216%, 处于高演化阶段, 以生气为主。高分辨率地震资料解释结果显示O , 文昌组在白云凹陷中面积达1900km 2, 厚度为17003000m, T OC 平均值为2194%, 氯仿沥青/A 0平均值为01225%, 有机质类型以1型为主, H I 平均为48314mg /g ; 恩平组在白云凹陷中分布面O 梁金强, 郭依群, 沙志彬, 等. 天然气水合物资源量评价方法及成矿远景研究. 广州海洋地质调查局, 2002. 图6 神狐海域line A 有机质成熟度模拟
25、Fig. 6 Organic maturity simulation of line A in Shenhu area积为2860km 2, 厚度为11002300m , T OC 平均值为2119%, 氯仿沥青/A 0平均值为011976%, 有机质类型以2 型为主。岩石热解分析测定恩平组烃源岩S 1+S 2为012234136mg/g, 平均为3111mg/g, H I 为411640010m g/g , 平均为15714mg /g 。研究区热解生气潜力同样巨大。313 流体运移模拟从测线剖面所在区域的文昌组和恩平组烃源岩产生的油气流体运移模拟结果可以看到(图7 , 深部的文昌组和恩平组烃
26、源岩已经开始产生大量的热解气, 并且所产热解气通过断层或上部渗透率高的岩层, 可以运移至浅部水合物稳定带, 为水合物成藏提供一定的热解气。同时也应注意到, 虽然深部烃源岩层能够大量产气, 但是大部分气体运移至珠海组和珠江组的有利构造部位集聚成藏, 然后, 这些成藏的气体部分以断裂为主要运移通道向上运移至浅部水合物稳定带。而当断层断至海底时, 气体将沿着断层逸散至海面, 造成气体的散失, 不利于水合物的成藏。另外, 深部热解图7 神狐海域line A 油气运移模拟7 simulation of oil gas of area第2期 苏丕波等: 南海北部神狐海域天然气水合物成藏动力学模拟 231
27、气也可以随超压孔隙流体向上运移, 与浅部生物气混 合形成水合物。而在浅部, 由于断裂构造不发育, 受流 体势控制, 浅部生物气主要以侧向运移的方式运移至 水合物稳定带区域。 远处, 有中国第一口深水钻井 L W3 1 1 井, 该井在上 -渐新统珠 海组 和下 中 新统 珠 江组 钻 遇了 大 量天 然 气, 累计探 明天 然气 地质储 量约 为( 800 1 100 10 m 。据此推测, 该区域深部烃源岩在一定程 度上可以产生大 量热解 气, 这些热 解气通 过合适 的 断层和底辟为天然气水合物的成藏提供一定的热解 气源。 同时, 近海油气勘探表明 , 南海北部边缘盆地 生物气的烃源岩分布
28、相当广泛, 纵向上从上中新统至 第四系, 甚至在局部区域的中中新统的不同层段也有 分布; 区域上盆地内均有大套浅海相 半深海相的泥 质烃源岩展布, 其有机质丰度相对较高, 已达到了作为 生物气烃源岩的标准, 具有一定的生烃潜力。并且已 在白云凹陷北斜坡 PY34 1 和 PY30 1 构造的浅层发 现了生物气气藏。 在上述分析基础上, 构建了该区的水合物成藏模 式( 图 8 : 神狐海域水合物气源由深海平原生物气和 深部热解气两部分组成。在硫酸盐 甲烷界面( SMI 之上是富硫酸盐的沉积层, 由于硫酸盐对产甲烷菌有 明显的抑制作用, 不可能存在甲烷水合物, 在 SMI 之 下, 由于硫酸盐被消
29、耗掉, 甲烷菌不受抑制作用, 在一 定条件下开始产生大量甲烷生物气并主要通过横向迁 移至有利部位形成水合物; 深部热解烃源岩具有良好 的生烃能力, 生成的大量气体以活动断裂为主要运移 通道向上运移, 并在合适的条件下在源岩上部有利构 造部位形成一定规模的天然气气藏。同时, 这些深源 高成熟气体持续以断裂为主要运移通道或者随超压孔 隙流体向上运移, 至浅部与浅部生物成因气混合在一 起, 在合适的温压域内形成水合物藏。 34 8 3 32 33 - 4 天然气水合物成藏模式的构建 天然气水合物成藏是一个复杂的过程。其成藏系 统包括烃类生成体系、 流体运移体系和成藏富集体系, 它们在时间和空间上的有
30、效匹配共同决定着天然气水 合物的成藏特征。始新世 早渐新世, 白云凹陷在潮 湿的气候环境、 全封闭的深洼陷及高的沉积速率下形 成了巨厚的文昌组、 恩平组烃源岩, 随后, 这两套烃源 岩在裂后相对构造平静期大量生烃, 而以高沉积速率 的深水细粒为主的充填作用导致白云凹陷形成超压, 随后的东沙运动使白云凹陷发育大型底辟构造和大量 NW 向张扭断裂, 压力随之得到释放, 逐步形成现今趋 于正常地层压力的状态 30 。超压存在说明油气运移 曾经不畅, 现今白云凹陷趋于正常压力, 则表明超压得 到了有效释放, 油气运移通畅, 大量油气 已经运移出 来。因此, 可以认为晚期底辟和断裂产生的垂向通道 为油气
31、垂向输导的有效通道。油气勘探也显示白云凹 陷北坡天然气藏具有晚期断裂控制成藏的特点, 同时 白云凹陷深水区存在的大量与底辟构造和断裂相关的 浅层亮点气异常反射, 也证明了凹陷深部的油气被垂 直输导到浅部地层, 即白云凹陷存在晚期活动的断裂 和底辟带的垂向输导系统, 可以大大改善天然气的垂 向运移条件。代一 丁等 31 通 过盆 地模 拟表明: 文 昌 组和恩平组两套烃源岩层在白云凹陷已处于产生裂 解气的阶段, 这与本次模拟吻合。另外, 离该测线不 图 8 神狐海域天然气水合物 成藏模式 Fig. 8 Reservoir f orming model for gas hydrate in Shenhu area 232 石 油 学 报 2011 年 第 32 卷 5 结 论 ( 1 神狐海域具备有利于天然气水合物成藏的水
