苏正1,2,曹运诚1,吴能友1,22,Lawrence ,陈多福1,2
苏正,(1980—),博士,助理研究员,主要从事天然气水合物及盆地流体活动的数值模拟研究,E-mail:。
注:本文曾发表于《地球物理学报》,2009,12:3124-3131,本次出版有修改。
1.中国科学院边缘海地质重点实验室/广州地球化学研究所,广州510640
2.中国科学院广州天然气水合物研究中心/可再生能源与天然气水合物重点实验室/广州能源研究所,广州510640
of Earth&Atmospheric Sciences,Cornell University,Ithaca,New York 14853-1504,USA
摘要:海洋环境中天然气水合物层是理想的毛细管封闭层,游离气被抑制在水合物层下,游离气层的气体压力随气体聚集和气层厚度的增加而升高,当气压超过封闭层的毛细管力时,游离气会克服毛细管进入压力、刺入上覆封闭层孔隙空间,毛细管封闭作用随之消失,从而形成水合物下伏游离气向海底的渗漏。通过对该过程进行的数值模拟计算表明:渗漏气体是以活塞式驱动上覆沉积层中的孔隙水向海底排出,水合物稳定带内流体渗漏速度随水流柱高度的减小而增加,当水流阻抗大于相应沉积层段的静岩压力时,沉积层将转变为流沙,流沙沉积被海流移除后便在海底留下凹陷麻坑。麻坑形成后流体运移通道演化为气体通道,气体快速排放。麻坑深度主要取决于游离气层的厚度和水合物封闭层(底界)的深度,而与沉积层的渗透率无关。麻坑深度一定程度上指示了渗漏前水合物层下伏游离气层的资源量。对布莱克海台海底麻坑的深度数值模拟计算表明,形成4 m深的海底麻坑需要至少22 m厚的游离气层。
关键词:天然气水合物;毛细管封闭;游离气渗漏;麻坑;布莱克海台
Numerical Computation and Case Analysis of the Venting Process of Free Gas Beneath Hydrate Layer
Su Zheng1,2,Cao Yuncheng1,2,Wu Nengyou1,2,Lawrence Duofu1,2
Key Laboratory of Marginal Sea Geology/Guangzhou Institute of Geochemistry,Chinese Academy of Sciences,Guangzhou 510640,China
Guangzhou Centerfor Gas Hydrate Research/CAS Key Laboratory of Renewable Energy and Gas Hydrate/Guangzhou Institute of Energy Conversion,Chinese Academy of Sciences,Guangzhou 510640,China
of Earth&Atmospheric Sciences,Cornell University,Ithaca,New York 14853-1504,USA
Abstract:A hydrate layer is an ideal capillary seal,beneath which free gas is overpressure increases as gas accumulates and gas column seals have the property that they fail completely when gas pressure reaches the point that they are invaded by gas,and thereafter they offer little resistance to gas the seepage is triggered,the venting gas will push the overlying water upward at increasingly higher velocities as the gas “piston”approaches the model shows that as the water velocity increases,the near surface sediments will become quick at a depth that the resistance of water flow exceeds the hydrostatic pressure of the sediment hosting the water quick sediments can then be removed by bottom ocean currents,leaving a hollow pockmark on the gas pathway isformed below the pockmarks and the reservoir gas drains pockmark depth is afunction of thickness of free gas column beneath the hydrate and depth of the hydrate seal (bottom of hydrate layer).Interestingly,pockmark depth does not depend on sediment depth implies the resource amount offree gas beneath hydrate model shows that a 22-m-thick free gas layer at least is needed toform a 4-m-deep pockmark on the rise of Blake ridge.
Key words:gas hydrate;capillary seal;gas seepage;pockrnarks; Blake ridge
0 引言
在海洋环境水合物稳定带内孔隙水溶解甲烷浓度超过甲烷水合物形成的溶解度时,溶解甲烷会结晶形成水合物,随着水合物含量的增加,形成水合物层圈闭,并在其之下发育游离气层[1-4]。在特定的条件水合物层之下的游离气沿通道向上渗漏进入海底,并在海底形成麻坑、自生碳酸盐岩、生物群落、气泡羽状体,如俄勒冈外海水合物脊[5]、布莱克海台等[6]、北刚果陆坡[7-8]、挪威外海[9]以及中国南海[10]。虽然水合物层下伏游离气向上渗漏活动在水合物发育区比较普遍,但是水合物层下伏游离气向上渗漏的机制和泄漏过程中的流体动力学特征,及流体渗漏对海底沉积地层的破坏(形成麻坑)过程并不清楚。
水合物层下伏游离气受到水合物层毛细管作用的封闭,随气体聚集和气层厚度增长,水合物下伏游离气的压力持续增加,当气体超压克服毛细管封闭作用后气体渗漏被激发,超压气体推动孔隙水向上排出,在海底形成麻坑,麻坑深度反映了流体的破坏强度和游离气层的超压幅度。因此,本文将应用水合物层毛细管封闭机理和沉积孔隙流体渗漏动力学,研究水合物稳定带之下游离气如何向上突破的动力学过程,建立游离气层压力状态与麻坑深度之间的数值模型,通过海底麻坑特征揭示水合物系统游离气层的演化规律。
1 毛细管封闭及游离气渗漏机理
海底沉积层中存在2种毛细管力封闭作用。第一类毛细管力封闭作用是存在于小型的气藏顶部的毛细管封闭作用,属于低渗透率的气体捕集封闭。封闭层的孔隙度和渗透率较低,而水更倾向存在于较小的孔隙空间,因此封闭层的孔隙空间完全被水占有,而封闭层之下含气层的孔隙度和渗透率相对较高[11]。碎屑沉积物孔隙介质一般为水润湿相,气液界面处的毛细管力阻止天然气进一步向上运移,使气体处于孔隙较大的沉积层段,但当气体压力超过相应孔隙的气体的毛细管进入压力时,超压气体将刺入封闭层的小孔隙,气藏开始排气,并在上覆沉积层中产生气体的渗漏通道。侵入毛细管压力由拉普拉斯方程给出[12]:
南海天然气水合物富集规律与开采基础研究专集
其中:γ为界面张力,取值 N/m[13],rf和rc分别代表小孔隙和大孔隙的有效孔隙半径。
第二类毛细管封闭作用存在于气-液二相共存的沉积孔隙中,气液二相均可流动[14-15]。由于整个沉积体是由沉积颗粒构成的孔隙介质,孔隙水优先占据并被吸附在孔隙的喉道位置,具有小孔径的孔喉部位产生的毛细管力抑制了孔隙腔中气体的流动。此类毛细管封闭条件是孔隙内2种流体共存,且二者均可流动。在渗漏活动初期这种情况出现在气流柱顶部和气柱周围的气-水混合的部位,沉积层中毛细管封闭线的位置随气柱的发育而变迁,这种毛细管封闭作用约束了气流柱的形状和发育,并使气流柱有一个相对平坦的顶部;同时也会形成一个相对稳定的通道直径,这意味着渗漏气柱顶部的气-水界面相对平坦,在理想均质介质中渗漏气体以“活塞”式向上推进。但是当渗漏气柱遇到渗透率在横向上不均匀或不连续(如断层)的沉积介质时会出现分支或扭曲的气体通道。
海洋环境扩散型水合物稳定带与下伏游离气之间属于第一类毛细管力的封闭,在水合物稳定带底部水合物含量最高[3,16],水合物的形成降低了孔隙介质的有效孔隙度和渗透率,使水合物层的孔隙度低于下伏游离气层的孔隙度,水合物层的有效孔隙半径小于游离气层的有效孔隙半径。亲水性的水合物沉积层内除水合物外的其余孔隙空间被水占据,而下伏沉积体的孔隙空间完全被气体充填,水合物层与游离气层之间就存在一个上覆孔隙水与下伏游离气的界面。因此在水合物层与游离气层界面(大孔隙与小孔隙之间)上产生毛细管力,其方向指向孔隙半径较大的含气层,阻止下伏气体进入上覆含水层(水合物层),抑制气体向上运移。但是当下伏游离气层中的气体压力超过上覆水合物封闭层的毛细管力时,超压气体将刺入水合物封闭层,使水合物层的毛细管封闭作用完全失效或仅剩很小的封闭作用,气体泄漏开始。超压的气体渗漏进入水合物稳定带后,随着气柱的增长气体逐渐侵占原有孔隙水所占的孔隙空间,驱使孔隙水向上排出,并最终泄漏进入海底。水合物稳定带内气柱的增长过程受第二类毛细管封闭作用的控制,使气流柱以“活塞”式增长,而没有出现气流弯曲和分支,这与地球物理资料显示的近于垂直的流体渗漏通道(气囱)特征一致[8-9,17-19]。
图1给出了海洋水合物层下伏游离气渗漏过程。游离气在水合物层底界之下聚集,气层厚度和气体超压逐渐增加(A),当气体压力超过水合物封闭层的毛细管力时,高压气体会在封闭薄弱点或气层最顶端刺穿封闭,使水合物毛细管封闭失效(B)。气流柱在高压作用下向上推进,并驱使上覆沉积孔隙水向外排出。气流柱高度(hg)逐渐增长,而水流柱高度(hw)相应缩短(B到C过程)。如果气压驱动力保持相对恒定,由于岩层对水的黏滞力(或水流阻抗)远大于其对气的黏滞力(或气流阻抗),随水流柱高度hw减小,流体渗漏速度将越来越快,在单位长度水流柱上的压降(等于岩层对水流的黏滞力)随流体速度的增长而增加。在气流接近海底时流体速度明显增强,浅层水流阻抗(即水流对地层的作用力)超过相应沉积体的静岩压力,浅层含水沉积将被流沙化,当流沙化的沉积物被海底底流搬运后,便在海底形成“新鲜的”麻坑,此时麻坑下形成单一的气体运移通道(D)。由于气体黏度远小于水的黏度(约为1/60),气体排放异常迅速,游离气藏中气体会很快排干,流体渗漏通道中的气流逐渐退化(E),孔隙流体压力回归静水压力,孔隙水重新占据水合物封闭层和流体渗漏通道的孔隙空间,在气量通量减小体系温度降低的过程中伴随者水合物的生成(此文中不做详细论述),并因此减小了流体流动速度,少量气体仍可滞留在流体渗漏通道内,在地震记录上显示为气烟囱,水合物层底部的毛细管封闭作用恢复,水合物层之下游离气的聚集过程再次启动(F)。
图1 水合物下伏游离气渗漏概念模型示意图[11]
Z为海底以下深度,h为水合物稳定带厚度(或水合物封闭层深度)。黑色带表示毛细管封闭层,浅灰色表示气体所占据孔隙沉积层。A.气体被封闭在水合物层之下;B.气体刺穿封闭层开始泄漏C.气柱高度增加,推动水流向外排出,水流柱高度相应缩短,流体运移速度不断增加;D.含水流沉积中孔隙压力超过静岩压力,在海底出现麻坑,形成单一的气流通道;E.游离气藏中的天然气被逐渐排空,孔隙超压消失,流体通道中的气流柱逐渐退化;F.气流柱完全消失,在海底留下气烟囱,并有水合物生成,水合物封闭作用恢复,并开始新的气体聚集
2 游离气渗漏过程的数学模型
气体渗漏过程中(图1)气柱和水柱都是在游离气超压的驱动下流动,流体运移的总驱动力等于气体超压(ρw-ρg)gd。气流柱不断增大,并且以同一速度推动渗漏通道内的上覆孔隙水向上流动。假定水合物稳定带为一种均质孔隙介质,渗漏通道内流体(水和气)的渗漏速率相同,孔隙介质内流体渗漏模型可用达西定律描述为
南海天然气水合物富集规律与开采基础研究专集
其中:Δp为流体运移总推动力,是施加在气流柱和水流柱上的压降之和(Δpg+Δpw),或者是气流阻抗与水流阻抗之和,等于气层底部的超压(ρw-ρg)gd;ρ为流体密度;d为游离气层的厚度;μ为流体黏度;V为流体速度;k为沉积体的渗透率;krg和krw分别为沉积体孔隙气和水的相对渗透率;hg和hw分别为气流柱和水流柱的高度。
假定气流柱中气的饱和度和水流柱中水的饱和度均为1,气和水的相对渗透率为1。由方程(1),流体(气体和水)的运移速度表示为
南海天然气水合物富集规律与开采基础研究专集
在方程(2)中,若 可知流体运移速度随气流柱高度(hg=h-hw)的增长而增加。对方程(2)进行积分得到气柱增长方程:
南海天然气水合物富集规律与开采基础研究专集
利用方程(3)既可以计算渗漏气流柱增长到某一高度所需要的时间,也可以计算某一时间点水合物稳定带内气流柱的高度。
由方程(1)和方程(2)可知,孔隙介质中单位长度流体柱所受阻抗随气流柱高度的增加(或水流柱高度的减小)而增加,也就是说沉积物格架所受流体的反作用力(流体阻抗)逐渐增加,当流体阻抗超过相应沉积体的静岩压力时,相应沉积层将被流体化而成为流沙[20],渗漏流体速度须满足 。流沙沉积被海流移除后在海底形成麻坑,被流沙化沉积体的底界确定了麻坑深度。用 替换方程(2)中流体速度V,麻坑深度hpm替换水流柱高度hw,即可得到麻坑深度方程:
南海天然气水合物富集规律与开采基础研究专集
方程(4)中,若μw≌60μg、krw≌krg≌1(假定水流柱中水的饱和度和气流柱中气的饱和度近似为1),方程(4)可简化为
南海天然气水合物富集规律与开采基础研究专集
在一定的温压条件下流体密度和黏度为常数[12]。因此,方程(5)中麻坑深度可近似为水合物下伏的游离气层厚度(d)和水合物封闭层深度(h)的函数,与沉积体的渗透率无关。模型计算中所有参数取国际标准单位。
3 模型应用及讨论
美国卡罗莱纳外海的布莱克海台区是典型的水合物发育区,既有完美的BSR显示,又有游离气的渗漏活动及在海底形成的麻坑[6,21-22]。大洋钻探计划(ocean drilling program)1 64航次对布莱克海台进行了钻探取心研究,其中997站位钻至海底之下750 m,穿过了BSR(海底之下450 m),其中180~462 m 层段含水合物,水合物平均饱和度为6%,位于水合物稳定带底部(462 m)的水合物体积分数最高为24%[4]。996站位于布莱克底辟链的最南端,处于997站位西北98 km,最大钻孔深度为63 m,刚好位于麻坑之中,地震剖面显示该区BSR深度为440 m,深部底辟作用使上覆地层变形、形成小型断层,成为有利的流体渗漏通道,在海底发育有深4 m、直径50 m的麻坑,并且正在发生气体渗漏(图2),钻探获得的水合物体积分数高,最高达沉积孔隙的99%[6,21-23]。
驱动流体运移的气体超压取决于游离气层的厚度。如果下伏游离气层厚度达100 m(图1),其总的流体驱动力(等于气体超压)可达到 MPa;如果游离气层厚度为22 m,流体超压驱动力为 MPa(图3最左端A点)。渗漏开始时水流柱高度分数(等于hw/h)为1,总水流阻抗等于气体总超压,整个气流柱高度增加而降低。但是由于水流速度增加,施加在单位长度水流柱上的驱动力和相应的黏滞力增加,水流阻抗逐渐趋近海底相应深度沉积层静岩压力,且在水流柱高度分别小于40 m(对于游离气层厚度为100 m)和4 m(对于游离气层厚度为22 m)时水流阻抗超过沉积介质的质量(图3D点)。该位置以上的沉积物被流沙化[20],转变成颗粒悬浮的液状混合体,这种流沙化沉积被海流搬运后在海底形成麻坑。利用方程(3)可以计算游离气从水合物稳定带底部渗漏到达海底所需的时间。假定渗漏率为10-12m2时, 100 m厚的游离气层泄漏到海底的时间大约为5 a。
图2 布莱克海台地震反射强度剖面揭示的BSR、底辟构造、海底麻坑及与ODP977站位揭示的BSR深度比较
a.地震反射强度显示布莱克海台水合物发育、气体聚集以及底辟构造顶端的流体渗漏[22];b.为ODP997站位BSR揭示的水合物封闭层深度[21]
图3 渗漏通道中的流体阻抗和含水沉积层的静岩压力曲线交点指示麻坑深度
水合物稳定带中气流柱高度增加(顶部坐标向右),水流柱高度减小(底部坐标向右),水流阻抗和静岩压力随之减小,水流阻抗大于静岩压力时发生流沙破坏,曲线交点位置指示麻坑深度(D点)。布莱克海台100 m的游离气层发生渗漏时在海底可形成40 m深的麻坑,而22 m厚的气层泄漏时可形成4 m深的海底麻坑(最右边灰色阴影)
方程(2)中流体渗漏速率与渗透率成正比,但方程(4)中麻坑深度不依赖于沉积体渗透率,只是水与气体相对渗透率比的函数,而相对渗透率决定于孔隙流体的饱和度[12],因此沉积体渗透率控制流体渗漏速率,但不控制麻坑形成。实际上,渗透率越大,气体渗漏越快,麻坑形成越快;气体超压在水流柱和气流柱之间的分配不依赖于渗透率,而是决定于气体的超压幅度,以及流体黏度和气流柱高度(或水流柱高度)。
利用方程(5)可以简单计算海底麻坑深度,同时在已知水合物底界(封闭层)深度和麻坑深度,也可以通过方程(5)计算游离气层的厚度。图4显示麻坑深度与游离气层厚度和封闭层深度的关系。在给定封闭层深度,麻坑深度随游离气层厚度的增加而增大,相反较深的沉积层厚度削弱了渗漏流体对麻坑的挖掘作用,水合物封闭层越浅,形成一定深度的麻坑所需的游离气层厚度越小。
图4 水合物封闭层深度和麻坑深度与游离气层厚度的关系
麻坑深度主要决定于游离气层厚度和水合物封闭层埋深,与游离气层厚度呈正比,与水合物层埋深呈反比。如果水合物封闭层深700m,形成4m深的麻坑需要27m的游离气层,如果水合物封闭深度为440 m,则需要22 m的游离气层,如果水合物封闭层深100m,仅需要1l m厚的游离气层
地球物理显示布莱克海台ODP996站位周围的BSR深度为440 m,而在ODP996站位正下方游离气藏气体沿底辟构造上升至大约220 m(图2)处,在沿小断层渗漏至海底,由方程(5)可知麻坑深度与渗透率无关,取决于游离气藏的埋深和游离气层的厚度。对于海底4 m深的麻坑,计算表明在水合物层之下至少需要有22 m厚的游离气层。苏正和陈多福[4]计算了布莱克海台997站位的水合物和游离气体积分数分布,在水合物稳定带底界之下26 m处的气体饱和度为28%,底界之下74 m处气体饱和度为%,其中水合物体积分数分布与同一区域的ODP995站位是相近的[24]。28%的气体饱和度大于气体流动所需20%的饱和度,而底界之下74 m处%的气体饱和度不能流动,也不能传递孔隙气体压力。如果20%的饱和度指示可传递气层的底界,则气层的有效压力传递厚度约为30 m,这与笔者22 m厚的游离气层模型计算结果相近(图5)。实际上,该钻位水合物平均体积分数约为6%[4],可封闭气层厚度为24 m(三角点所示),接近模型估计的22 m。此外,在水合物稳定带底部的水合物饱和度达24%[4],其毛细管作用可封闭约33 m的游离气层(菱形点所示),与Flemings等[25]估计的极限破坏厚度29 m相似(虚线所示位置),接近但略小于30 m的参考厚度。然而,在996站位游离气发生泄漏后, 997站位扩散型水合物的体积分数仍在持续增加[26],水合物层的封闭能力也相应增强,游离气层厚度不断增长,因此,997站位游离气厚度(30 m)大于996站位游离气发生泄漏时的22 m气层厚度是合理的。
图5 布莱克海台的水合物饱和度和所能封闭的游离气层厚度
气层厚度随水合物饱和度增加而增高,水平虚线与气层厚度曲线的交点(29 m)为Flemings等预测的997站位气层的临界水力压裂厚度[25],圆形点标示约30 m的实际气层厚度,三角形点显示平均饱和度6%的水合物能封闭24m的气层,而饱和度24%的水合物可封闭33 m的游离气层(菱形点)
4 结语
本文构建了水合物层下伏游离气渗漏动力学过程的数学模型,游离气被水合物层的毛细管作用所圈闭,下伏游离气的超压随游离气层的增长而增加;当气体超压超过作用于水合物与游离气层界面的毛细管阻力时,游离气渗漏进入上覆水合物稳定带,并以“活塞式”驱动上覆孔隙水向外排出,渗漏速度随水流柱高度的减小而增加;当水流阻抗超过相应层段的静岩压力时沉积体变为流沙,流沙沉积被海流带走便在海底留下麻坑。模型显示麻坑深度为游离气层厚度和水合物封闭层埋深的函数,而与沉积介质的渗透率无关。游离气渗漏形成的海底麻坑对水合物下伏游离气层的厚度具有指示作用,在已知水合物封闭层深度和海底麻坑深度条件下,模型可以计算水合物层下伏游离气藏发生渗漏时的气层厚度,在布莱克海台海底发育有4 m深的麻坑,它的形成需要至少22 m厚的游离气层。
致谢:挪威国家石油公司Martin Hovland教授提供了全球麻坑基础资料和最新信息,表示感谢。
参考文献
[1]Xu W ,Ruppel the Occurrence,Distribution,and Evolution of Methane Gas Hydrate in Porous Marine Sediments[J].Journal of Geophysical Research,1999,104:5081-5095.
[2]Davie M K,Buffett B Steady State Model for Marine Hydrate Formation:Constraints on Methane Supply from Pore Water Sulfate Profiles[J].Journal of Geophysical Research,2003,108(B10): 2495,doi:.
[3]Chen Duo-Fu,Su Zheng,Cathles L of Gas Hydrates in Marine Environments and Their Thermodynamic Characteristics[J].Terrestrial Atmospheric and Oceanic Sciences,2006,17(4) :723-737.
[4]苏正,陈多福.海洋环境甲烷水合物溶解度及其对水合物发育特征的控制[J].地球物理学报,2007,50(5): 1518-1526.
[5]Trehu A M,Bohrmann G,Rack F R,et ofthe Ocean Drilling Program,Initial Reports Volume 204[M].TX:Ocean Drilling Program,2003.
[6]Paull C K,Spiess F N,Ussler W Ⅲ,et Plumes on the Carolina Continental rise: Associations with Gas Hydrates[J].Geology,1995,23: 89-92.
[7]Sahling H,Bohrmann G,Spiess V,et in the Northern Congo Fan area,SW Africa: Complex seafloor features shaped by Fluid Flow[J].Marine Geology,2008,249 : 206-225.
[8]Gay A,Lopez M,Berndt C,et Controls on Focused Fluid Flow Associated with Seafloor Seeps in the Lower Congo Basin[J].Marine Geology,2007,244 (1/2/3/4):68-92.
[9]Hovland M,Svensen H,Forsberg C F,et Pockmarks with Carbonate-Ridges off Mid-Norway:Products of Sediment Degassing[J].Marine Geology,2005,218:191-206.
[10]陈多福,李绪宣,夏斌.南海琼东南盆地天然气水合物稳定域分布特征及资源预测[J].地球物理学报,2004,47(3):483-489.
[11]Cathles L in Sub-Water Table Fluid Flow at the End of the Proterozoic and Its Implicationsfor Gas Pulsars and MVT Leadzinc Deposits[J].Geofluids,2007,7(2): 209-226.
[12]Bear of Fluids in Porous Media[M].New York:Elsevier,1972.
[13]Vigil G,Xu Z,Steinberg S,et of Silica Surfaces[J].J Colloid Interface Sci,1994,165:367.
[14]Cathles L Seals as a Cause of Pressure Com-partmentation in Sedimentary Basins:Presented at the Gulf Coast Section SEPM Foundation 21 st Annual Research Conference on Petroleum Systems of Deep-Water Basins,2001:561-572.
[15]Shosa J D,Cathles L Investigation of Capillary Blockage of Two-Phase Flow in Layered Porous Media,in Petroleum Systems of Deep-Water Basins: Global and Gulfof Mexico Experience: Proceedings ofthe GCSSEPM st Annual Bob Research Conference,2001:725-739.
[16]苏正,陈多福.海洋天然气水合物的类型及特征[J].大地构造与成矿学,2006,30(2).
[17]Hovland M,Judd A Pockmarks and on Geology,Biology and the Marine Ltd.,1988.
[18]Hovland M,Svensen Pingoes:Indicators of Shallow Gas Hydrates in a Pockmark at Nyegga,Norwegian Sea[J].Marine Geology,2006,228:15-23.
[19]Gay A,Lopez M,Cochonat P,et so1ated Seafloor Pockmarks Linked to BSRs,Fluid Chimneys,Polygonal Faults and Stacked Oligocene-Miocene Turbiditic Palaeochannels in the Lower Congo Basin[J].Marine Geology,2006,226(1/2):25-40.
[20]Nicholl M J,Karnowski Apparatus for the Demonstration of Quicksand[J].Journal of Geoscience Education,2006,54(5): 578-583.
[21]Matsumoto R,Paull C,Wallace Leg 164 Scientific party[C]//Gas hydrate sampling on the Blake Ridge and Carolina Rise: ODP,Leg 164 Preliminary Report,1996.
[22]Taylor M H,Dillon W P,Pecher I and Migration of Methane Associated with the Gas Hydrate Stability Zone at the Blake Ridge Diapir:New Insights from Seismic Data[J].Marine Geology,2000,164:79-89.
[23]Paull C K,Matsumoto R,Wallace P J,et [C].164:TX:Ocean Drilling Program,2000.
[24]王秀娟,吴时国,刘学伟.天然气水合物和游离气饱和度估算的影响因素[M].地球物理学报,2006,49(2):504-511.
[25]Flemings P B,Liu X,Winters W Pressure and Multiphase Flow in Blake Ridge Gas Hydrates[J].Geology,2003,31: 1057-1060.
[26]Liu X,Flemings P Multiphase Flow Model of Hydrate Formation in Marine Sediments[J].Journal of Geophy sical Resea rch,2007,112,B03101,doi:.
陈芳1,刘广虎1,苏新2,周洋1,陆红锋1,刘坚1,王金莲1
陈芳(1966—),女,教授级高级工程师,主要从事微体古生物学和海洋地质研究。
1.广州海洋地质调查局,广州510760
2.中国地质大学海洋学院,北京100083
摘要:为探讨沉积物粒度与水合物饱和度的关系,对南海神狐海域水合物钻探区的2个获取水合物的钻孔岩心沉积物进行了粒度分析及粒度与水合物饱和度对比分析。结果表明:水合物主要赋存于粉砂中,含水合物沉积物具有φ(粉砂)为>70%,φ(黏土)小于介于15%~30%之间,φ(砂)一般小于10%的基本特征,其中粉砂中以8~32μm和32~63μm粒级的中细-粗粉砂占优势;含水合物层中砂、粗粉砂含量高的层位与水合物饱和度高的层位呈良好的对应关系,说明粗粒沉积物更有利于水合物的形成和发育。
关键词:沉积物粒度;含水合物层;神狐海域;南海
Variation in Grain Size of Gas Hydrate-Bearing Sediments and the Correlation of Gas Hydrate Saturation from Shenhu Area in the Northern South China Sea
Chen Fang1,Liu Guanghu1,Su Xin2,Zhou Yang1,Lu Hongfeng1,Liu Jian1,Wang Jinlian1
Marine Geologic Survey,Guangzhou 510760,China
of Ocean sciences,China University of Geosciences,Beijing 100083,China
Abstract:Grain size of sediments from two drill holes of Shenhu gas hydrate drilling-area from northern South China Sea are analyzed to understand the relationship between the variation in grain-size of sediments and the saturations of gas results suggest that the siltfraction is the dominant component of gas-hydrate-bearing sediments,fealuned by with> 70‰f second component is clay,in a range of 15%~30%,whereas the sandfraction is less than 10%.The Grain size of gas-hydrate-bearing sediments is similar to the one of no gas-hydrate-bearing contents of coarse-grains silt(~) and sand(~) can be correlated with sediment layers where gas hydrates saturations were further suggests that the occurrence of gas hydrates in Shenhu aera is mainly correlated with coarse-grams sized sediments.
Key words:sedimentary grain size;gas-hydrate-bearing sediments; Shenhu Area;South China Sea
0 引言
水合物的形成与分布除了需要特定的温压条件外,更需要合适的沉积条件,以提供充足的气体来源和良好的储集条件。作为水合物存在载体之一的海洋沉积物,其岩性是除温压条件外控制水合物成藏的重要因素。岩性的差异影响着水合物的产状与饱和度,一般地,沉积物越粗,饱和度越高。各海域已发现水合物的水合物稳定带沉积物的岩性各有差异,但相同的是水合物稳定带内沉积物粒度总体相对较粗[1-5]。笔者以我国在南海神狐海域成功钻取水合物钻孔岩心为材料,分析研究南海北部水合物稳定带沉积物的粒度特征及其与水合物饱和度的关系,探讨沉积物粒度对水合物的制约机制。
1 取样与方法
所研究的2个含水合物钻孔SH2B和SH7B沉积物样品,由广州海洋地质调查局于2007年在南海北部陆坡神狐海域实施“我国海域天然气水合物钻探”航次调查所获得。含水合物层岩心取样主要采用非保压、保压的FC、FRPC和FPC方式。鉴于该钻探航次的实际需要等原因,钻探只在关键和部分控制层段钻取岩心。因此,粒度分析样品取样间隔差异较大,介于20~68 m之间,两钻孔共取样176个。钻孔取样位置见参考文献[6]。沉积物粒度分析方法遵照中华人民共和国国家标准《GB/T 海洋调查规范第八部分:海洋地质地球物理调查》执行,粒级标准采用尤登-温德华氏等比制φ值粒级标准,粒度参数计算采用福克和沃德公式。分析方法采用Mastersizer2000型激光粒度仪:取沉积物样品数克置于玻璃杯中,加纯净水适量使样品充分浸泡,浸泡12 h使样品充分分散;加5 m L的 mol/dm3的六偏磷酸钠(〔Na PO3〕),再浸泡12 h,将浸泡充分分散的样品搅拌均匀,取适量加入激光样品槽中,加超声振动和高速离心,使样品再次充分分散,测定各级粒级质量分数。激光粒度分析误差相同粒级差小于3%,均符合国家标准要求,可以满足本次研究的需要。沉积物粒度分类和命名采用谢帕德的沉积物三角图解法分类方案。样品的处理和测试工作是在广州海洋地质调查局测试所完成。
沉积物原位结构扫描电镜分析在广州海洋地质调查局测试所完成,仪器型号为捷克产的VegaⅡ LUM。
2 结果
水合物储层的地层分布
南海水合物钻探航次应用目前世界水合物勘查中多种新的测试方法和手段来探查水合物的存在和分布。如测井获得的异常电阻率记录、红外线图像(IR images)分析得到的温度差值记录、孔隙水低氯离子浓度的记录等。对沉积物岩心的IR扫描、样品在水中的分解过程观察及X射线扫描均提供了SH2B孔和SH7B孔地层中水合物存在的直接证据。取心后证实水合物呈分散状分布在黏土质粉砂和粉砂的孔隙中,肉眼难以观察到,水合物分解后沉积物呈粥状。通过对测井资料、取心资料以及地震资料的详细分析,确定出水合物层的分布区间;其中SH2B孔水合物分布的区间约在海底以下191~225 m之间,SH7B孔水合物分布的区间大体在海底以下155~177 m之间。根据生物地层的分析,水合物分布在上中新统—下上新统含钙质生物的黏土质粉砂和含钙质粉砂中,硅质生物放射虫和硅藻缺失(图1)[6]。
图1 神狐海域水合物储层的地层分布[6]
含水合物层粒组类型含量变化
沉积物粒度按照大小划分为3个粒组类型:> mm的颗粒统称为砂,~ mm的颗粒称为粉砂,< mm的颗粒则统称为黏土。两钻孔沉积物中最主要的粒组类型均是粉砂,粉砂平均体积分数介于%~%之间。含水合物层沉积物粒组类型平均体积分数与其上下层位沉积物的差别不大,但粒组类型体积分数的范围值有差异(表1)。以SH2B孔粉砂体积分数为例,该孔不含水合物层的粉砂体积分数介于%~%之间,而含水合物层的介于%~%之间。含水合物层沉积物具有粉砂大于70%,黏土量介于15%~30%之间,砂一般小于10%的基本特征。
表1 神狐海域含水合物层及其相邻层位沉积物粒组类型体积分数φB/%
含水合物层粒级组分分布
采用尤登-温德华氏等比制φ值粒级标准细分法将两钻孔沉积物粒级划分为以下10个粒级:粗砂> mm,中砂~ mm,细砂~ mm,极细砂~ mm,粗粉砂~ mm,中粉砂~ mm,细粉砂~ mm,极细粉砂~ mm,粗黏土~ mm和细黏土< mm,以便进行更细致地讨论。
SH2B孔含水合物层粒级组分以中粉砂、细粉砂和极细粉砂粒级为主,平均体积分数分别为%、%和%;粗粉砂占%;砂含量偏低,细砂—粗砂未见,仅见极细砂,平均体积分数%;黏土以粗黏土为主,为%。SH2B孔水合物主要分布在中粉砂、极细粉砂和细粉砂沉积物中,总体上,与含水合物层上下层位相比粒级组分体积分数差别不大(图2)。
SH7B孔含水合物层的沉积物较SH2B孔含水合物沉积物要粗,以粗粉砂、中粉砂和细粉砂粒级为主,体积分数分别为%、%和%;各粒级的砂均有出现,以中砂和极细砂为主,分别为%和%;黏土以粗黏土为主,占%。与含水合物层上下层位相比,粗粉砂、中粉砂粒级沉积物明显增加,粗粉砂、中粉砂和细粉砂是SH7B钻孔含水合物沉积物组成的最主要颗粒组分(图3)。但总的来说,在含水合物层的粉砂粒级中,以8~32 mm和32~63 mm这两个粒级占优势。
图2 SH2B孔各粒级组分体积分数分布(%)(阴影部分为含水合物层,下同)
图3 SH7B孔各粒级组分体积分数分布
含水合物层粒度与水合物饱和度的关系
图4 SH2B孔含水合物层粒度体积分数与水合物饱和度的对比
根据测井参数(LWD)推算出来的水合物饱和度值,理论上代表了沉积物中孔隙被水合物充填的程度。利用测井电阻率(LWD-RAB)和测井沉积物孔隙度值,根据Archie方程可以从理论上推算水合物饱和度值(Sh)[7-8]。这一方法在ODP 164航次和ODP 204航次得到充分实践和验证[9-10]。同样的,利用该方法计算出神狐海域水合物饱和度值。神狐海域含水合物层饱和度值变化范围较大,介于%~%之间,不同层位水合物饱和度差异明显。这种差异与沉积物粒度差异相关性较强。SH2B、SH7B孔水合物层沉积物砂、粗粉砂含量高的层位与水合物饱和度高层位有良好地对应关系(图4,图5A,B),即沉积物中砂、粗粉砂含量高,水合物饱和度也高,反之亦然。这种特征在SH7B孔表现得尤其明显,如在1 594~16 663 cm层段,砂含量%~%,饱和度为20%~44%,平均32%; 16 840~17 120 cm层段,砂含量下降,为%~%,饱和度随之下降,为2%~7%,平均4%; 17 160~17 585 cm层段,砂含量上升为%~%,饱和度随之也上升,为8%~23%,平均17%;粗粉砂与水合物饱和度的关系与砂和饱和度的关系相似(图5A,B)。这种粗粉砂粒径与水合物饱和度关系在3个单层中显示更为明显(图5B)。
图5 SH7B孔含水合物层粒度体积分数与水合物饱和度的对比
3 讨论
对南海神狐水合物钻探区钻井含水合物层岩性特征的研究发现,含水合物层岩性与上下不含水合物层位的差异不大。因此,只要温压、气源等条件满足,在南海海底以下任何深度都可能形成水合物。但就南海神狐海域含水合物层而言,粗粒沉积物砂、粗粉砂含量高的层位与水合物饱和度高的层位呈良好地对应关系,说明沉积物的颗粒粒径是影响水合物形成的控制因素。推测粗粒沉积物可以增加沉积物的孔隙度,为水合物的形成提供更大的孔隙空间,这一点在对含水合物层沉积物原位结构研究中得到证实。根据扫描电镜的观察,沉积物中的砂、粉砂和黏土颗粒随机分布,黏土充填在砂、粉砂颗粒间;砂主要由有孔虫和条状形黄铁矿组成,而粉砂主要由不规则的石英和长石组成,黏土主要充填于颗粒间孔隙中(图6),其含量的高低影响孔隙度的发育程度。换言之,粗粒沉积物含量高,黏土含量相对降低,有利于孔隙的发育,黏土由于其黏性和密实性不利于孔隙的发育。而且粗粒沉积物渗透性好,有利于气体的运移和储存。但水合物饱和度与沉积物的孔隙度并未完全呈正相关关系,含水合物层沉积物为松散未固结沉积物,实测的沉积物孔隙度自下而上逐渐增加,而水合物饱和度表现为时高时低,说明孔隙度只是控制水合物饱和度的因素之一;水合物饱和度还受气体通量、孔隙类型和大小、沉积物渗漏性等因素的影响。初步研究发现,南海神狐水合物钻探区钻井含水合物层沉积物的孔隙主要有粒间孔隙和粒中孔隙2种类型,实测的沉积物孔隙度主要由粒间孔隙组成,而粒中孔隙主要存在于沉积物中的有孔虫房室中。由于实测的沉积物孔隙度无法测得粒中孔隙,因此,水合物饱和度表现为时高时低,除了受粒间孔隙影响外,很大程度受粒中孔隙的影响。关于这方面的深入研究结果另有文章介绍。
图6 SH7B孔16835-16860cm原位沉积物颗粒的分布结构
4 结论
对南海神狐海域钻取的含水合物的2个钻孔岩心进行沉积物粒度分析及其与水合物饱和度的对比分析,结果表明:
1)含水合物层沉积物相对较粗,其基本特征为以(含)钙质生物粉砂为主,粉砂含量为>70%,黏土含量介于15%~30%之间,砂一般小于10%。
2)含水合物层沉积物的粉砂粒级,以8~32 μm和32~63μm粒级的中、细—粗粉砂占优势。粗粒沉积物砂、粗粉砂含量高的层位与水合物饱和度高的层位呈良好地对应关系,说明沉积物的粒度是水合物形成的重要控制因素之一。粗粒沉积物有利于孔隙的发育和水合物的形成。
参考文献
[1]Ginsburg G,Soloviev V,Matveeva T,et Grain-Size Control on Hydrate Presence,Sites994,995 and 997[C]//Paul1C K,Matsumotor R,Wallace P of ODP Initial Reports Station,2000:237-245.
[2]Shipboard Scientific 204 summary[C]//TR Hua M,Bohrmann G,Rack F R,et ODP Init :College Station,2003:1-75.
[3]Kraemer LM,Owen R M,Dickens G of the Upper Gas Hydrate Zone,Blake Outer Ridge,a Link Between Diatoms,Porosity,and Gas Hydrate[C]//Paul IC K,Matsumotor R,Wallace P of ODP Inital Reports Station,2000:229-236.
[4]苏新,宋成兵,方念乔.东太平洋水合物海岭BSR以上沉积物粒度变化与气体水合物分布[J].地学前缘,2005,12(1):234-242.
[5]王家生,高钰涯,李清,等.沉积物粒度对水合物形成的制约:来自IODP航次证据[J].地球科学进展,2007,22(7):659-665.
[6]陈芳,苏新,周洋,等.南海北部陆坡神狐海域晚中新世以来沉积物中生物组分变化特征及意义[J].海洋地质与第四纪地质,2009,29,(2):1-8.
[7]Archie G Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics[J].American Insititute of Mining,Metallurgica and Petroleum Engineers,1942,146:54-67.
[8]Collett T Log Evaluation of Gas Hydrate Saturations[C]//SPWLA 39th Annual Logging SPWLA,1998.
[9]CollettT S,Ladd of Gas Hydrate with Downhole Logs and Assessment of Gas Hydrate Concentrations (saturations) and Gas Volumes on the Blake Ridgewith Electrical Resistivity Log Data[C]//Paull C K,Matsumoto R,Wallace P of ODP,Science Results, Station,2000:179-191.
[10]Trehu AM,Long P E,Torres M,et Distribution of Gas Hydrate Beneath Southern Hydrate Ridge:Constraints from ODP Leg 204[J].Earth Planet Science Letters,2004,222:.
梁劲1王宏斌1,2梁金强1
(1.广州海洋地质调查局 广州 510760;2.中国地质大学(北京)北京 100083)
第一作者简介:梁劲,男,1971年生,高级工程师,1995年毕业于成都理工学院信息工程与地球物理系应用地球物理专业,主要从事天然气水合物调查与研究工作。
摘要 本文采用Jason 反演技术对南海北部陆坡A 测线纵波速度进行计算,结合BSR、振幅空白带以及波形极性反转等多种水合物赋存信息的分析,对水合物成矿带的速度特征进行了综合研究,结果表明:低速背景中的高速异常,是天然气水合物赋存的重要特征;高速异常体一般呈平行于海底的带状分布;在高速异常的内部,速度也是不断变化的。一般在异常体的中心速度最高,由中心到边缘速度逐渐降低,反映在水合物矿带内部,水合物饱和度由矿体中心向边缘逐渐降低的特征。本文的研究成果进一步表明高精度速度分析不仅可以帮助寻找水合物矿点,还可以进一步判定水合物的富集层位。
关键词 Jason 反演技术 天然气水合物 速度分析
1 前言
天然气水合物是在低温、高压环境下,由水的冰晶格架及其间吸附的天然气分子组成的笼状结构化合物,广泛分布于海底和永久冻土带。温度和压力是天然气水合物形成和保存最重要的因素(王宏斌等,2004)。针对天然气水合物的野外调查及研究表明:高分辨率的地震勘探方法是天然气水合物调查评价中行之有效的方法。地震反演技术一直是地震勘探中的一项核心技术,其目的是用地震反射资料反推地下的波阻抗、速度、孔隙度等参数的分布,从而估算含天然气水合物层参数,预测天然气水合物分布状况,为天然气水合物勘探提供可靠的基础资料。常用的地震反演技术有Jason、Strata、Seislog和ISIS等,其中Jason反演技术在含天然气水合物层预测中因其分辨率高而得到广泛推崇,它主要由有井约束和无井约束两种方法组成(廖曦等,2002)。
速度异常是判断天然气水合物是否赋存的重要条件之一。结合BSR(Bottom Simulating Reflector)特征、波形极性特征、振幅特征以及AVO特征等目前已成为判断是否存在天然气水合物层主要手段(史斗等,1999)。大量的测试数据显示:水合物的速度与冰的速度较为接近,而比水高。与含水或含游离气沉积层相比,含水合物沉积层的密度降低,声波速率增大,含水合物层的地层速度往往比一般的地层速度高,含水合物沉积层的下部由于充填了水或气,而使水合物底界面出现速度负异常。因此,地层中速度反转是水合物赋存的一个地球物理标志。含水合物地层的声波速度与水合物的含量有关,水合物含量越高,其声波速度越高。从速度方面看,BSR是上覆高速的含水合物地层与下伏较低速的含水层或含气层之间的分界面。通常,海洋中浅层沉积层的地震纵波速度为1600~1800m/s,如果存在水合物,地震波速度将大幅提高,可达1850~2500m/s,如果水合物层下面为游离气层,则地震波速度可以骤减200~500m/s。因此,在速度剖面上,水合物层的层速度变化趋势呈典型的三段式,即上下小、中间大的异常特征(张光学等,2000)。西伯利亚麦索雅哈气田的资料表明,在原为含水砂层内形成水合物之后,其纵波的传播速度会从1850m/s提高到2700m/s;而在胶结砂岩层,这种速度会从3000m/s提高到3500m/s。深海钻探计划的570站位的测井结果表明,由含水砂岩层进入含水合物砂岩层时,密度由降低到,声波传播速度从1700m/s提高到3600m/s,且电导率剧烈下降。
Cascadia海域ODP889站位的VSP测井资料反映水合物底界为强烈的负速度界面,速度从水合物沉积物层的1900m/s陡降到含游离气层的1580m/s,由于VSP测井为地震测井,受钻井因素的影响较少,因此认为VSP测井真实地反映了水合物沉积层底界的速度变化(陈建文等,2004)。
国土资源部广州海洋地质调查局在2001~2004年在南海北部陆坡进行10000多公里的天然气水合物高分辨地震调查。本研究利用Jason反演技术,通过对南海北部陆坡区的地震速度资料的精细分析,在已圈定BSR分布范围的基础上研究陆坡区各沉积层的速度特征,最后对速度值与水合物的关系进行了分析和探讨。
2 方法原理
纯天然气水合物的密度()和海水密度相近,而游离气的含量又十分有限,这就决定了产生BSR的波阻抗差主要由速度造成。速度反演技术的特点是在无井约束时,以地震解释的层位为控制,对所有的地震同相轴来进行外推内插来完成波阻抗反演,这样就克服了地震分辨率的限制,最佳的逼近了测井分辨率,同时又使反演结果保持了较好的横向连续性。速度反演技术的主要原理是:①通过最大的似然反褶积求得一个具有稀疏特性的反射系数系列;②通过最大的似然反演导出波阻抗;③通过波阻抗计算速度。该方法的主要优点是能获得宽频带的反射系数,是一种基于模型的反演,具有多种建模方法,对所建模型进行比较分析,并使地质模型更趋合理,反演结果更加真实可靠(郝银全等,2004)。
波阻抗反演方法的出发点是认为地下的反射系数是稀疏分布的,即地层反射系数由一系列叠加于高斯背景上的强轴组成。具体反演是从地震道中,根据稀疏的原则抽取反射系数,与子波褶积生成合成地震记录,利用合成地震记录与原始地震道的残差修改反射系数,得到新的反射系数序列,然后再求得波阻抗。其具体步骤是:
假设地层的反射系数是较大的反射界面的反射和具有高斯背景的小反射叠加组合而成的,根据这种假设导出一个最小的目标函数(安鸿伟等,2002):
南海地质研究.2006
式中:R(K)为第一个采样点的反射系数,M为反射层数,L为采样总数,N为噪音变量的平方根,λ为给定反射系数的似然值。
最大的似然反演就是通过转换反射系数导出宽带波阻抗的过程。如果从最大的似然反褶积中求得的反射系数式R(t),则波阻抗:
Z(i)=z(i-1)×(1+R(i))/R(1-i) (2)
利用波阻抗和速度的关系式:
v=Z(i)/ρ (3)
即可得到速度值。其中,ρ为地层密度,可从区域测井资料结合该测线重力资料反演求取。
在上述过程中为了得到可靠的反射系数估算值,可以单独输入波阻抗信息作为约束条件,以求得最合理的速度模型。一方面,速度反演结果是一个宽频带的反射序列和波阻抗及速度数据,同时加入了低频分量,使反演结果更能正确反映速度变化规律;另一方面,它有多种质量控制方法,具体表现为监控子波的选取、同相轴的连续追踪、反演结果准确性的判断和提供多种交汇显示的相关性分析。所以利用速度反演可对地震剖面上任一相位进行速度反演,在每一个CDP点都可得到任一个同相轴速度数据,并利用二维的反射波的速度层析成像反演方法得到高度连续的速度剖面,如果地震测线足够密,还可利用三维速度反演得到速度体图像。
3 实现过程
初始模型的确立
在地质规律的指导下,利用地震和测井资料开展沉积特征分析和沉积旋回划分;建立岩石-电性关系,进行砂层组和单砂层对比;在地震剖面上提取各含油砂层组反射波属性,建立地震属与矿体的关系,实现地震-测井综合预测矿体平面分布厚度,开展层间矿体组外推预测;建立初始速度场;在地震属性约束下开展地震反演,反演层间小层矿体厚度。细分层反演层位的标定正确与否直接影响反演结果的精度。因此,在反演过程中对子波提取、能谱特点、信噪比、频谱及反射系数的研究至关重要(闫奎邦等,2004)。技术路线流程如图1所示:
初始速度场的获得
初始速度场的获得首先要对速度谱进行解释,速度谱的解释和取值是否合理,将直接影响均方根速度的计算精度。具体步骤如下:
1)速度谱的解释先从地质条件简单、反射层质量好、能量团强、干扰少的剖面段开始,绘制叠加速度-反射时间曲线,并逐渐向外扩展;
2)结合地震剖面的反射特征,判断速度极值点是否正确,并选择读取能量团最大的极值点。排除干扰波能量团,从而求得有效波的叠加速度;
3)对相邻速度谱进行比较,通过比较速度谱曲线的形状、相同反射层的速度极值等方法予以检查和修改。
4)每隔40个CDP拾取一组数据,利用地震剖面上的反射倾角数据对它们进行校正,便可得到均方根速度(梁劲等,2006)。
图1 速度反演技术线路流程图
The flow chart of the velocity inversion of technical route
子波的提取
子波提取时,要使能量集中于子波的主瓣,与地震子波形态吻合。如果所提子波近于零相位,则从波峰向两侧能量衰减较快,波峰两侧波形对称;在子波的能谱特征分析,要使能量都集中在地震波的主频范围内;有井资料时,要对井资料都作了子波与地震波自动关联质量控制。保证子波能谱与地震波能谱相吻合,是反演中较为重要的一方面,子波能谱的峰值与地震波主频的能谱峰值相吻合。首先了解合成记录与地震记录之间的偏差。通过合成记录与地震记录之间的偏差分析,对Jason反射系数偏差、能谱偏差进行进一步的校正,使合成记录与地震记录之间的偏差减小。然后通过反射系数与地震资料之间偏差分析,采取相应的手段校正,使地层与合成记录反射系数相吻合。再进行信噪比分析,使反演处理后的信噪比得到最大限度的提高。通过一系列质量控制手段,使各油层合成记录与地震记录的标定精度得到了较大的提高。
关于速度反演可信程度,不能完全由反演方法确定,关键在于获取地震记录的质量和反演前处理流程的振幅保真度。另一个影响因素是数值模拟结果应当是比较准确的,这与计算方法有关,也与子波拾取和地质构造模型有关。至于反演结果的灵敏度,主要由拟合误差值和收敛速度来判断。如果给定的初始模型正确,即与实际地质结构一致,则拟合的误差较小且收敛速度快。本文工作由于受实际情况限制,没有实际的测井资料验证,因此反演所得速度的准确性和精度会受到一定程度的影响。
4 速度剖面特征
运用多种特殊地震成像综合分析,是天然气水合物地震资料解释的关键技术。目前一般采用识别BSR、振幅空白带、波形极性反转、速度异常、波阻抗面貌和AVO等天然气水合物地震相应特征来综合分析沉积物中是否含有水合物。高精度的层速度分析可帮助判定水合物的富集层位,速度及振幅异常结构是水合物与下伏游离气共同作用形成的特殊影像,剖面上表现为“上隆下坳”结构,多层叠合构成一明显的垂向“亮斑”这一特殊成像结构在未变形的水合物盆地内较适用于寻找水合物矿点,并可据此定量估算水合物盆地内水合物的数量,分析BSR上下的详细速度结构,是水合物地震资料综合解释的重要手段(张光学等,2003)。
图2 南海北部陆坡测线A道积分剖面
Trace integration profile of the line A in north slope of the South China Sea
图2是南海北部陆坡测线A的地震反射道积分剖面,从图中可以看出,该剖面中部及右下角距海底大约350ms处出现一强振幅反射波,大致与海底反射波平行,与地层斜交,BSR特征明显。在波形极性方面,海底反射波和BSR都表现为成对出现的强振幅双峰波形特征,海底反射波表现为蓝红蓝特征,而BSR表现为红蓝红特征,这表明相对于海底,BSR显示出负极性反射同相轴,即所谓的极性反转(与海底反射相反)。反射波的极性是由反射界面的反射系数决定的,而反射系数则与界面两侧的波阻抗差有关。实际上,海底和BSR都是一个强波阻抗面,海底是海水和表层沉积物的分界面,上部为低速层,下部为相对高速层,反射系数为正值;BSR是含水合物层与下部地层(或含气层)的分界面,上部为高速层(水合物成矿带是相对高速体),下部为相对低速层(如含游离气,则速度更低),反射系数为负值,因此造成了BSR和海底反射波的极性相反现象(沙志彬等,2003)。图3是用速度反演法反演出来的纵波速度剖面,该速度剖面明显显示出一近似平行于海底的相对高速地质体,其位置恰好在BSR上方。高速地质体的纵波速度大约在2000~2400m/s,其上面的低速层的纵波速度大约在1500~1800m/s,而下面的低速层的纵波速度大约在1500~1900m/s,没有明显的游离气存在特征,但根据其高速地质体特征、BSR以及波形极性反转分析,可以认为南海北部陆坡测线A的相对高速地质体极可能是水合物成矿带。
图3 用速度反演法计算的南海北部陆坡测线A纵波速度剖面
P velocity profile of the line A in north slope of the South China Sea computed by velocity inversion
由图3可见,水合物成矿带内部速度是变化的,表明水合物分布不均匀,呈平行于海底的带状分布,中心速度最高,由中心到边缘速度逐渐降低。海底以下有3个近似平行海底的低速和高速带:①海底与高速体之间的相对低速带,为水饱和带;②水合物成矿带;③水合物成矿带下的低速带。水合物成矿带下面的低速带在速度剖面上没有明显的低速特征,由此推断水合物成矿带下可能不含游离气,或者是气体的饱和度很低。
5 结论
水合物的生成除了需要一定的温度和压力条件外,还需要大量的碳氢气体和充足的水。这就需要地层具有较高的孔隙度和渗透率。未固结沉积岩的孔隙度很高,渗透率大,具备水合物生成的物理条件。具备这种特征的未固结沉积岩的地震波速度较低,而含水合物地层的地震波速度增大。这就形成了水合物成矿带作为低速背景中的高速地质体特征。另外,水合物的生成受温度和压力控制,一般情况,等温面和等压面近似平行于海底,因此低速背景中近似平行于海底的相对高速地质体是水合物成矿带的特征(刘学伟等,2003)。
通过对南海北部陆坡A测线纵波速度的计算,并且结合BSR和振幅空白带识别以及波形极性反转等多种特殊地震成像进行综合分析,我们可以进一步了解水合物成矿带的速度特征:揭示水合物成矿带的高速异常一般呈平行于海底的带状分布,在高速异常的内部,速度也是不断变化的,一般在异常体的中心速度最高,由中心到边缘速度逐渐降低,该现象反映在水合物矿带内部,水合物分布并不均匀,水合物饱和度由矿体中心向边缘逐渐降低。分析BSR上下的详细速度结构,是水合物地震资料综合解释的重要手段。高精度速度分析可帮助判定水合物的富集层位,较适用于寻找水合物矿点,并可据此估算水合物资源量。
参考文献
安鸿伟,李正文,李仁甫,等.2002.稀疏脉冲波阻抗反演在XY油田开发中的应用.石油物探,41(1):56~60
陈建文,闫桂京,吴志强,等.2004.天然气水合物的地球物理识别标志.海洋地质动态,6:9~12
郝银全,潘懋,李忠权.多井约束反演技术在油气储层预测中的应用.成都理工大学学报,31(3):297~300
梁劲,王宏斌,郭依群.2006.南海北部陆坡天然气水合物的地震速度研究[J].现代地质,20(1):123~129
廖曦,马波,沈浩,等.2002.应用Jason软件进行砂体及含气性预测.天然气勘探与开发,25(3):34~42
刘学伟,李敏锋,张聿文,等.2005.天然气水合物地震响应研究——中国南海HD152测线应用实例.现代地质,19(1):33~38
沙志彬,杨木壮,梁金强,等.的反射波特征及其对天然气水合物识别的应用.南海地质研究,15(1):55~61
史斗,郑军卫.1999.世界天然气水合物研究开发现状和前景.地球科学进展,14:330~339
王宏斌,梁劲,龚跃华,等.2005.基于天然气水合物地震数据计算南海北部陆坡海底热流.现代地质,19(1):67~73
闫奎邦,李冬梅,吴小泉.反演技术在岩性识别中的应用.石油物探,43(1):54~58
张光学,黄永样,陈邦彦,主编.2003.海域天然气水合物地震学.北京:海洋出版社
张光学,文鹏飞.2000.南海甲烷水合物的地震特征研究,首届广东青年科学家论坛论文集,中国科学技术出版社
The Application of Jason Inversion Technology in Velocity Analysis of Gas hydrate
Liang Jin1Wang Hongbin1,2Liang Jinqiang1
( Marine Geological Survey,Guangzhou, University of Geosciences(Beijing),Beijing,100083)
Abstract:The P velocity of A seismic profile in the north slope of the South China Sea were calculated by Jason inversion velocity characterostic of the gas hydrate bed was researched in detail based on the calculated result and the information of gas hydrate existing including BSR,amplitude blanking and polarity reversion of the shows that:The abnormity of higher velocity in the background of lower velocity is an important characteristic of gas hydrate existing;The abnormity of higher velocity which distribute as a belt usually parallel to the seafloor;The velocity changes gradually at the inner of the abnormity of higher velocity with the highest velocity at the center of the abnormity whereas the lowest velocity at the margin of it,which suggests that the saturation of gas hydrate decreases gradually from the center to the result that mentioned above suggest that high resolution velocity analysis not only help to search the hydrate spot but also help to estimate the rich layer of gas hydrate.
Key Words:Jason Inversion Technology Gas hydrate Velocity Analysis
173 浏览 1 回答
120 浏览 4 回答
115 浏览 2 回答
105 浏览 5 回答
188 浏览 3 回答
253 浏览 3 回答
174 浏览 3 回答
131 浏览 3 回答
104 浏览 3 回答
220 浏览 6 回答
300 浏览 6 回答
182 浏览 3 回答
207 浏览 2 回答
181 浏览 3 回答
118 浏览 4 回答