MANTLE PLUME(HOTSPOT)-RIDGE INTERACTION
-
摘要: 地表热点是认识地幔柱假说以及地幔柱动力学的窗口,发生在洋脊与热点之间的相互作用则是了解地球上两大动力系统(板块构造和地幔柱)的直接联系的最有利地区。研究洋脊-热点之间的相互作用对于揭示地幔动力学、热点附近洋壳构造的演变以及与热点密切相关的洋中脊处的岩浆热液活动具有重要的意义。在肯定地幔柱假说的基础上,对洋脊-地幔柱(热点)的模拟实验以及三大洋中不同扩张脊与热点相互作用的最新研究作了系统的介绍和评述,指出室内模拟实验以及地质学、岩石学、地球化学和地球物理学(特别是高分辨率的地震技术)的结合研究将是本领域研究的发展趋势。
-
关键词:
- 地幔柱 /
- 洋中脊 /
- 热点 /
- 热点-洋中脊相互作用
Abstract: Hotspot is a window to understand mantle plume hypothesis and dynamics of mantle plume, and the area where ridge interactions with hotspot is the best place to directly find out relationship between plate tectonics and mantle plume. Based on affirming mantle plume hypothesis, the authors introduce several 2D or 3D simulation experiments about ridge-plume (hotspot) interaction and some examples of hotspot-ridge interactions existing in the three oceans. It is further pointed out that simulation experiments combined with geology, petrology, geochemistry and geophysics(especially for high resolution seismic technique) in studying mantle(hotspot)-ridge interaction will play an important role in such reseaches as plume-ridge interactions.-
Keywords:
- mantle plume /
- mid-oceanic-ridge /
- hotspot /
- plume-ridge interaction
-
天然气水合物是一种由充足的气源、水源在高压(一般大于10 MPa)和低温(一般小于10 ℃)条件下络合而成的笼形结晶化合物[1-2],主要分布在陆地的永久冻土带和大陆边缘水深超过300 m的陆坡带[3-4]。水合物存在的区域往往会伴随着游离气储层,水合物作为游离气体的盖层,二者共同成藏[5]。天然气水合物具有能量密度高、分布广、规模大、埋藏浅等特点,被认为是21世纪的重要新能源[6],地质历史演变中的天然气水合物演化对全球变化与全球碳循环起着重要作用[7]。
中国天然气水合物主要分布在高纬度地区的冻土层和南海海域的浅部地层中。当前对南海水合物资源量的研究表明,其北部陆坡区有相当可观的水合物资源[8]。中国陆地和海域的水合物赋存类型有所不同,祁连山冻土区水合物主要以裂隙型和孔隙型两种状态产出,而海域生物成因气多发育为薄层状或分散状水合物,深部热解气则常依附断层形成厚层状水合物[9]。受成岩影响,中国冻土区水合物储层孔隙、裂隙发育相对海域来说整体性要差;而海洋水合物储层孔隙、裂隙普遍发育,整体表现为高孔、低渗和强非均质性的典型特点[10]。Dai等[11]在研究墨西哥湾深水区天然气水合物过程中提出了颗粒间胶结、颗粒包裹、颗粒支撑、孔隙填充、基质和填充物共存、裂缝充填物等6种微观分布模式。范金等[12]通过分析这6种主要赋存状态的水合物,总结了每种赋存状态之间存在一定的相互转化关系,研究了各种水合物赋存地层的地震响应特征。陈鸣等[13]通过高分辨率随钻成像技术,指出琼东南盆地天然气水合物具有两种宏观赋存状态,即分散状水合物及薄层状水合物,且这两种水合物均具有一定非均质性;水合物的微观赋存状态主要包括孔隙填充型及骨架支撑型;水合物储层岩性主要为泥质粉砂岩及粉砂质泥岩。樊奇等[14]通过水合物的岩芯分析,揭示了琼东南盆地水合物富集区储层具有“细粒性、未成岩、易溃散、高电阻、高声速、强非均质性”等六大特征;同时,识别出泥质粉砂孔隙型储层和南海新发现的水合物储层类型泥质粉砂微裂缝型储层两类主要储层。
地震方法是水合物储层识别中最重要也最常见的方法之一,如潘冬阳等[5]采用叠后地震属性分析方法对南海北部陆坡神狐海域地震数据进行了天然气水合物储层的识别,研究表明地震“三瞬”属性在BSR和天然气水合物储层的识别上具有明显的效果,可以为圈定水合物成藏的目标区域提供依据。孟大江等[15]利用叠前反演技术对南海北部陆坡神狐海域天然气水合物以及气态烃储层进行了预测,研究表明纵横波速度交会是研究区识别储层最敏感的弹性参数,纵横波速度识别因子识别储层效果较好,能够清晰刻画水合物和气态烃储层,其预测的天然气水合物矿体与钻井情况吻合较好。
虽然国际上对天然气水合物产状、分布和特征的认识已取得显著进展,并开展了一系列陆地多年冻土区和海域天然气水合物试采,但天然气水合物开采仍面临诸多问题和科学挑战[16],如天然气水合物成因成藏类型与气源供给及产出赋存特点、天然气水合物成藏机理及主控因素等[17]。所以,需要深入研究水合物的成藏特点及赋存特征,探寻新的水合物识别方法,进而为水合物的勘探和开发提供有意义的方法指导。
Ikelle 等[18]于1993年根据统计学理论,利用指数型自相关函数模拟了二维随机扰动,最早建立了二维随机介质模型,由此开启了有关随机介质研究的先河。Ergintav 和 Canitez[19]于1997年在Ikelle 等研究的基础上,采用高斯型自相关函数建立了二维随机介质模型,同时进行了正演模拟。Müller和Shapiro [20]于1999年根据散射理论实现了地震波在二维及三维随机介质中传播的格林函数构造。国内对随机介质理论的研究从2001 年开始,由奚先和姚姚等[21-22]率先利用椭圆自相关函数模拟了二维平稳随机介质模型,且对模型进行了正演模拟及波场分析,后续进一步对非平稳随机介质理论、复杂随机介质模型模拟以及弹性介质中的波场特征分析等开展了一系列研究[23-25]。李灿苹等[26]将随机介质理论应用到冷泉羽状流模型的构建中,并对其进行了地震波正演模拟,探讨了羽状流的地震波场特征,对羽状流进行叠前偏移成像,成像效果良好;后期对模型进行了改进[27],最终获得的模型更接近实际羽状流赋存特征。目前,国内还未见将随机介质理论应用到水合物储层中的报道。
中国海域水合物储层具有明显的强非均质性特点,因此可以借助随机介质理论对水合物储层进行数值模拟研究,储层内水合物团块区域随机分布在地层中,将其视为大尺度非均匀性;而团块区域内部水合物随机充填在岩石孔隙、裂隙和缝隙之间,将其视为小尺度非均匀性。由此构建的水合物储层模型更加贴合实际赋存状态。因此,本文基于水合物储层和游离气储层的实际赋存特征,根据随机介质理论,建立了上覆均匀地层、中间水合物储层和游离气储层、下覆均匀地层的四层介质起伏界面的速度模型。通过该模型与实际水合物和游离气储层的对比,四层模型可以较好地模拟水合物储层和游离气储层内水合物和游离气分布的复杂非均匀性。该研究为下一步更好地了解水合物成藏特点和赋存特征,研究水合物储层的地震响应特征奠定基础,从而为水合物识别提供理论方法指导。
1. 随机介质的基本概念及其算法实现
将具有大量的、随机分布的、小尺度异常的复杂非均匀介质看成随机介质,介质中的非均匀体可看成一个空间随机过程,用统计学方法描述介质非均匀性所形成的非均匀介质模型就是随机介质模型[28]。
随机介质模型由非均匀性大、小尺度两部分所组成。大尺度非均匀性描述介质的平均特性,是传统意义上的地质模型;小尺度非均匀性是加在地质模型上的随机扰动,通常使用一个均值为零的二阶平稳的空间随机过程来表示,在小尺度上的空间扰动用几个统计量(空间自相关函数、相关长度、均值及标准方差)来描述[29]。
在二维声波随机介质中,空间坐标点$(x,z) $处的速度$ v(x,z) $[30]可分解为:
$$ v(x,z) = {v_0} + \delta v(x,z) $$ (1) 其中,$ {v_0} $是大尺度非均匀介质参数,可设为常数或缓慢变化;$ \delta v(x,z) $是小尺度非均匀扰动量。速度的相对扰动会呈现出线性特点,小尺度上的非均匀性可以用一个速度的扰动量来代表[30]。设
$$ \sigma (x,z) = \delta v/{v_0} $$ (2) 式中,$ \sigma (x,z) $为介质空间的相对扰动,$ \sigma = \sigma (x,z) $假设是具有零均值及自相关函数、方差的空间二阶平稳随机过程,由式(1)和(2)可得:
$$ v(x,z) = {v_0} + \delta v(x,z) = {v_0}(1 + \sigma ) $$ (3) 根据随机过程理论可知[29],随机过程$ \sigma (x,z) $的功率谱就是由自相关函数$ \varphi (x,z) $进行傅里叶变换得到的频谱函数$ \phi ({k_x},{k_z}) $;根据随机过程的谱展开理论,可对功率谱函数$ \phi ({k_x},{k_z}) $进行二维随机运算,得到随机功率谱函数$ \Psi ({k_x},{k_z}) $;对随机功率谱函数$ \Psi ({k_x},{k_z}) $做傅里叶逆变换便可获得由它描述的随机扰动$ \sigma (x,z) $;最后规范化$ \sigma (x,z) $的均值和方差,产生以$ \varphi (x,z) $为自相关函数、具有指定均值及方差的随机扰动[31],代入式(3)便可产生随机介质模型。
本文选用指数型自相关函数产生随机扰动,其具体算法流程参见文献[31],本文不再赘述。
2. 水合物和游离气储层特点及建模思路
2.1 水合物和游离气储层实际赋存状态
地层中水合物成藏需要满足诸多条件,如水合物成藏的环境需要低温高压,水合物储层要有充足的气源以及大量的水,且拥有流体运移和汇聚的通道和空间[10]。在水合物稳定区下部,水合物不能稳定存在,游离气和水存在水合物层下部,一般被称为游离气区[32],所以游离气储层一般是伴随水合物共同成藏的。
水合物在地层中有多种不同的赋存模式[11-12],如图1所示,有像固体颗粒状镶嵌在地层中的,有像流体一样填充在地层孔隙中,这两种情况下探测到的水合物区域速度相对于不含水合物地层的速度,前者速度增加幅度较大,后者速度增加幅度较小。因此,水合物的饱和程度直接关系到在水合物储层中水合物的形态及其高低速区域分布特点。
图2是南海北部陆坡神狐海域由地震属性分析得到的天然气水合物储层[5],图中黑色的BSR(Bottom Simulating Reflector,似海底反射)界面上层有明显的三个高饱和度的水合物区域(黄色椭圆区域),BSR下方伴随着大量的游离气,即水合物储层作为游离气储层的盖层赋存于游离气储层的上面,而水合物储层上层和游离气储层下层为较均匀的地层;水合物和游离气赋存区域为团块状,且团块大小不一(几十米至几百米),其中大块水合物为长条不规则的形状;水合物储层与游离气储层、水合物储层与上层均匀介质的分界面是上下起伏的。图3是中国海域实际钻探获得的典型水合物样品,能够明显看出乳白色的水合物以不同赋存状态随机分布在黑色的岩层中[10]。中国海域生物成因气多发育为薄层状或分散状水合物,深部热解气则常依附断层形成厚层状水合物[9],其中神狐海域水合物主体呈分散状(图3d、f)和厚层状交替发育[10]。
图 3 中国海域典型水合物样品[10]a:块状水合物, b:瘤状水合物, c:脉状水合物, d:粉砂质泥岩中分散状水合物, e:层状水合物, f:泥质粉砂岩中分散状水合物。Figure 3. Typical hydrate samples from Chinese seas[10]a: Massive hydrates, b: nodular hydrates, c: vein hydrates, d: dispersed hydrates in silty mudstone, e: layered hydrates, f: dispersed hydrates in muddy siltstone.综合海底水合物赋存状态(图1)、水合物储层(图2)和水合物样品(图3),都显示出水合物是随机分布在海底沉积物中,其形态也是不规则的,有层状,有分散状态;而水合物储层下方存在游离气,游离气亦随机充填在岩石孔隙和裂缝中。该描述是针对水合物和游离气的赋存区域内部,即小尺度范围。
2.2 水合物储层和游离气储层速度模型构建
根据图2中神狐海域水合物和游离气储层分布特点,即水合物上层为常规地层,水合物下层常富含游离气,游离气下层为常规地层,由此,构建符合实际水合物赋存的四层地层架构图(图4),图中储层内水合物团块区域随机分布在地层中,此为大尺度;同理,游离气储层内游离气团块区域也随机分布在地层中,亦为大尺度。
实际测井资料及地震探测结果显示[33-34],海底浅层地震纵波速度为
1600 ~1800 m/s,含水合物沉积层的纵波速度约为2300 m/s,游离气地层速度约为1400 m/s。由此设定了图4中每一层的地层速度,以及水合物团块区域、游离气团块区域的背景速度。图1中水合物赋存的6种状态显示,水合物在地层中的微观分布具有非均匀性(此即小尺度范围),水合物存在区域相当于随机介质,即水合物随机分布在固体颗粒、基质、结核/裂缝之间。所以,当地震波传播到含有水合物的地层时,水合物区域将产生散射,该区域地层速度应不是纯水合物速度,而是处于水合物和地层速度之间的一个波动速度,从含水合物地层的实际测井曲线[33]可以得到证实。因此,水合物团块区域内部速度的实现可借助于随机介质理论的速度构建方法。同理,由于游离气是充填在岩石孔隙或地层裂隙中,而孔隙或裂隙分布具有随机性或非均匀性,所以,其游离气赋存区域的速度实现也可借助于随机介质理论的速度构建方法。
3. 利用随机介质理论模拟水合物储层
根据以上水合物储层赋存特点及模型构建思路,可以建立一个水合物赋存区随机分布的速度模型。成藏在海底沉积物中的水合物区域,其形态可以模拟为横向不规则的长条状,使其更接近海底沉积物中水合物的赋存状态。为此,先构建水合物储层赋存区轮廓模板(简称储层模版),再结合以水合物背景速度建立的二维随机介质模型,从而实现水合物储层速度模型。
3.1 水合物储层模板
储层模板既要使水合物的赋存状态合理,还需明确水合物的大小。参考中国海域水合物的赋存形态及水合物稳定带厚度(110~350 m)及水合物层厚度(10~69 m)[10],设定水合物团块区域在横向长度上最大不超过200 m,在纵向深度上最大不超过50 m。根据中国海域水合物储层具有明显的强非均质性特点,并参考水合物钻探岩芯样品[10],设置大团块区域较分散,大都呈长条状水平分布,也有倾斜分布,与层状和脉状水合物赋存区域相对应;较大团块中间还存在小团块,与块状和分散状水合物赋存区相对应。构建模板的过程为,先将水合物模板和随机介质模型设置成大小相同的数组,然后在模板中构建出水合物的形态及大小,初步实现水合物模板的大概轮廓;在轮廓图的基础上,将轮廓图转换为二维数组,并对图片进行截取及调整大小,然后读入修改后的图片,识别出水合物团块;最后对水合物区域和非水合物区域分别赋0和1,由此构建出水合物储层模板,结果如图5所示,图中黑色区域即为富含水合物区域。
3.2 水合物储层二维随机介质速度模型
根据2.2节水合物储层模型构建思路,水合物团块内部速度为随机介质的速度,所以,图4中水合物团块区域是二维随机介质。为此,首先建立该随机介质,然后再将其与水合物储层模板匹配,从而实现真正的水合物储层。
在构建水合物层的二维随机介质模型时,模型的规格与水合物储层的模板大小必须相同,以保证水合物储层的模板可以运用在二维随机介质模型上。根据第1节中建立随机介质的方法,给定模型参数:自相关长度a=b=1 m,标准差0.05,步长1 m,背景速度v0=
2300 m/s(即水合物储层中水合物团块区域的背景速度),最后得到在水合物背景速度下扰动的二维随机介质模型(图6)。从图6中可以看出,介质中速度围绕背景速度上下扰动,非均匀性特点明显。3.3 水合物储层速度模型
根据前文所述水合物的赋存特征,水合物为固体颗粒时,检测到的速度相对地层的速度增加幅度较大,水合物像流体一样填充在地层时,检测到的速度相对地层的速度增加幅度较小,所以,水合物团块区域的速度也是或大或小,具有波动变化的特征。因此,水合物储层要赋予一个较小的地层速度v01=
1900 m/s,则水合物团块的背景速度赋予v0=2300 m/s,使围绕水合物团块背景速度v0产生随机扰动的二维随机介质模型在速度上大于地层速度v01,以符合实际探测结果[34]。结合图5水合物储层模板和图6二维随机介质模型建立水合物储层速度模型,两个数组满足水合物储层内及水合物团块的参数要求。水合物储层模板的数组元素中只存在两种数值,即水合物团块区域内元素数值为0,水合物团块区域外元素数值为1。因此,将二维随机介质模型加载到水合物储层的模板中,对模板中水合物团块区域内元素的数值赋予同行同列中二维随机介质模型元素的值,而水合物团块区域外元素的数值则赋予地层速度v01,最终得到的模型即是水合物储层速度模型,结果如图7所示。
从图7可以看出,不同规格不同形态的水合物团块区域随机分布在储层中,水合物团块内部模拟出了水合物在地层中处于不同赋存状态时的或高或低的速度[33],速度变化范围与实际探测的富含水合物地层相符[34]。整体呈现出水合物赋存区域在地层中分布的不均匀性特征,与实际的赋存状态比较接近。
4. 利用随机介质理论模拟游离气储层
由2.1和2.2节可知,水合物储层底界面下层往往有大量的游离气体,而游离气是充填在岩石孔隙或地层裂隙中,所以处于超高压底界面下的游离气赋存区可以模拟为大小不一的块状区域,形似充满孔隙的岩石或裂隙结构,使其更接近游离气的赋存状态。为此,同样先构建游离气储层赋存区轮廓模板(简称储层模版),再结合以游离气背景速度建立的二维随机介质模型,实现游离气储层速度模型。
4.1 游离气储层模板
与水合物储层模版类似,构建游离气储层模板时,首先将游离气模板和随机介质模型设置成相同规格的数组;其次要明确游离气在储层中的赋存形态,运用圆和椭圆函数在数组中进行块状游离气区域的划分。具体实现过程为:先构建一个零值的200×600数组,游离气团块区域赋予一个不同于游离气地层的值以使游离气的区域在图形中显示出来;参考游离气的赋存特征,游离气赋存区团块大小要小于大型水合物赋存区[10],因此,游离气团块区域范围不超过50 m×50 m。游离气团块之间分布较密集,随着深度变浅,压力变小,较大的游离气团块聚集在较浅位置,较小的游离气团块分布于较深位置。最后将游离气区域赋值为0,非游离气区域赋值为1,由此完成游离气储层模板的构建,结果如图8所示,图中黑色区域即为富含游离气区域。
4.2 游离气储层二维随机介质模型
根据2.2节游离气储层模型构建思路,游离气团块内部速度为随机介质的速度,所以,图4中游离气团块区域也是二维随机介质,并且以含游离气地层速度v1=
1400 m/s为背景速度。为此,首先建立该随机介质,然后再将其与游离气储层模板匹配,从而实现真正的游离气储层。游离气储层二维随机介质模型的构建方法与水合物储层相同,模型参数相同,只是背景速度不同,按图4设定背景速度,最后实现的随机介质模型如图9所示。同样,从图9中可以看出,介质中的速度围绕背景速度上下扰动,非均匀性特点明显。
4.3 游离气储层速度模型
纯水合物和游离气速度与地层速度相比较时,在含水合物地层中,水合物相比其地层速度是增加的;而在含游离气地层中,游离气相比其地层速度是减少的,所以块状游离气区域的背景速度v1要小于游离气储层地层的速度v11。又因游离气储层位于水合物储层的下层,所以块状游离气周围的地层速度v11要大于水合物的地层速度v01=
1900 m/s,因此游离气储层中地层速度设定为v11=2100 m/s。参照水合物储层速度模型的构建方法,构建游离气储层速度模型。结合图8游离气储层模板和图9二维随机介质模型,对应的两数组满足游离气储层层内及块状游离气的参数要求。游离气储层模板数组中只存在两种数值,即数值为0的块状游离气区域和数值为1的非块状游离气区域。因此,将二维随机介质模型加载到游离气储层的模板中,对模板中块状游离气区域内的数值0赋予同行同列中二维随机介质模型的值,而非块状游离气区域的数值1则赋予地层速度v11,最终得到的模型即是游离气储层速度模型,结果如图10所示。
从图10可看出,不同规格不同形态的游离气团块随机分布在储层中,游离气团块内部模拟出了游离气在地层中充填在岩石孔隙或地层裂隙中或高或低的速度[33](气充填的多,则速度降的多),速度变化范围与实际探测的富含游离气地层相符[34];整体上呈现出游离气团块区域在地层中分布的不均匀性特征,与实际的赋存状态比较接近。
5. 富含水合物和游离气储层速度模型
实际富含水合物地层,往往与游离气共同成藏,而且水合物作为游离气的盖层,处于游离气之上。因此,根据以上建立的水合物储层模型和游离气储层模型以及实际地层赋存特点,建立富含水合物和游离气储层的地层模型,并与实际赋存的水合物储层进行比对,讨论所建立模型的合理性。
5.1 四层介质速度模型
在地层中,水合物储层和游离气储层之间的地层接触不是简单的直线型,温度和压强的微小变化,都会促使水合物的分解,引发水合物发育区的上升和下降,形成新的底界面,所以,水合物储层与游离气储层之间的地层接触呈上下起伏状。水合物一般埋藏较浅,所以水合物储层深度设置为海底100 m以下;水合物稳定带厚度为200 m左右[10],所以水合物储层模型规格为深度200 m、长度600 m,因为大型水合物富集区需要有充足的气源,因此可以使构建的游离气储层模型与水合物储层模型大小相同。为更接近实际水合物所处的地层环境,在水合物储层之上和游离气储层之下分别设置一层简单的均匀速度地层,两层速度见图4所示。两地层均设置为深度100 m、长度600 m,同时两地层与水合物储层、游离气储层之间的接触关系为起伏曲线型,由此所构建的模型更符合实际地层分布。
按以上要求以及图4架构图,利用第3节和第4节构建的水合物储层速度模型和游离气储层速度模型,最终实现的四层介质起伏界面速度模型结果如图11所示。
5.2 四层介质速度模型与实际水合物和游离气储层的比对
将四层起伏界面速度模型(图11)与水合物在海底的赋存状态(图1)、由实际地震数据反演的水合物储层(图2)以及实际开采出的水合物岩芯样品(图3)进行对比,结果如下:
(1)模型中水合物储层和游离气储层内水合物和游离气赋存区团块随机分布,符合图2中水合物和游离气大尺度宏观范围的复杂非均匀分布的赋存状态。
(2)模型中两层之间有非常明显的起伏分界面,且水合物储层中只存在水合物,游离气储层中只存在游离气,与图2中明显的水合物层及游离气层的分界面,即图中标识的BSR相吻合。
(3)模型中水合物团块区域内部速度围绕纯水合物背景速度上下波动,这与实际水合物在地层中处于不同赋存状态时或高或低的速度相符;游离气团块区域内部速度也是上下波动,与游离气在地层中充填在岩石孔隙或地层裂隙中或高或低的速度相符[33]。此即小尺度范围的非均匀,与图1和图3中水合物和游离气的微观分布模式相对应。
(4)水合物团块区域内部的速度变化范围与实际探测的富含水合物地层相符;游离气团块区域内部的速度变化范围与实际探测的富含游离气地层相符[34]。
(5)模型中水合物赋存区域整体速度比游离气赋存区域速度高,符合实际探测情况[33-34]。
(6)通过设置水合物和游离气储层模型中随机介质的自相关长度,便可以模拟出符合实际开采的水合物尺度。
基于以上,图11较好地模拟了发生变化后水合物区域下形成的新游离气区域和高压底界面,且水合物和游离气赋存区域满足横向不规则长条的形状及块状游离气的条件。因此,四层速度模型在大尺度上很好地模拟了富含水合物地层中实际水合物和游离气赋存区域,在小尺度上实现了富含水合物地层中水合物和游离气分布的复杂非均匀性(即赋存区内部介质速度的非均匀性特点)。
6. 结论
(1)水合物储层和游离气储层速度模型:根据实际赋存的水合物和游离气两层结构特点及水合物和游离气的分布特点和变化特征,描绘出水合物和游离气轮廓图,由此构建出具有长条状水合物赋存区的水合物储层模板及小块状游离气赋存区的游离气储层模板。根据水合物区域的速度为正异常指示,游离气区域的速度为负异常指示,设定水合物储层和游离气储层速度参数;利用随机介质理论建立了水合物和游离气储层的随机介质模型。最后将二维随机介质模型加载到对应的模板中,由此建立了水合物储层速度模型及游离气储层速度模型。
(2)四层介质起伏界面速度模型:根据水合物和游离气之间地层界面是起伏变化的赋存特点,利用已建立的水合物储层模型和游离气储层模型,将二者之间设置成上下起伏的界面。为更接近实际水合物所处地层环境,在水合物层之上及游离气层之下分别设置了均匀速度地层,且与水合物储层及游离气储层接触关系为幅度较大的曲线型,从而建立了四层介质起伏界面速度模型。
(3)四层速度模型与实际水合物储层和游离气储层的比对:通过对比发现,模拟的水合物储层与实际水合物储层中水合物赋存状态大致相同,都为横向长条不规则的复杂非均匀体;两者的游离气赋存区,也都为复杂非均匀体,但模拟的游离气为块状游离气且只存在于水合物储层之下,而实际水合物储层中水合物的周围存在游离气;水合物和游离气储层的速度变化范围与实际探测的地层速度相符,且水合物储层速度高于游离气储层速度。因此,四层速度模型符合水合物区域赋存特点。
综上,基于随机介质理论建立的四层介质起伏界面速度模型可以方便灵活地描述富含水合物地层中水合物及游离气赋存区域以及内部的复杂非均匀性,能够较好地模拟实际水合物存在于固体颗粒、基质、结核/裂缝之间的状态以及游离气充填在岩石孔隙或地层裂隙中的特点,从而可以为研究水合物和游离气储层的地震响应特征提供地质模型。
-
[1] Wilson J T. Evidence from oceanic islands suggesting movement in the Earth[J]. Philos. Trans. R. Soc.London Ser. A, 1965,258:145-167.
[2] Morgan W J. Plate motions and deep convection[J]. Geol. Soc. Am. Mem., 1972,132:7-22.
[3] Morgan W J. Rodriguez, Darwin, Amsterdam:a second type of hotspot island[J]. J. Geophys. Res.,1978, 85:5355-5360.
[4] Anderson D L. The thermal state of the upper mantle:no role for mantle plume[J]. Geophys. Res. Lett., 2000,27:3623-3626.
[6] Chen Y, Morgan W J. A nonlinear-rheology model for mid-ocean ridge axis topography[J]. J. Geophys. Res., 1990,95:17571-17581.
[7] 陈永顺.海底扩张和大洋中脊动力学[M]//地球的结构、演化和动力学.高等教育出版社,2002:283-329.[CHEN Yong-shun. Seafloor spreading and dynamics of Mid-Ocean ridges[M]//ZHANG You-xue, YIN An. Structure Evolution and Dynamics of the Earth. Higher Education Press, 2002:283 -329.]
[8] Wilson J T. A possible origin of the Hawaiian Islands[J]. Canadian J.Phys.,1963,41:863-870.
[9] Morgan W J. Convective plumes in the lower mantle[J]. Nature, 1971, 230:42-43.
[10] Burke K C, Wilson J T. Hot spots on the Earth's surface[J].J. Geophys. Res.,1976,93:7690-7708.
[11] Crough S T, Jurdy D M. Subducted lithosphere, hotspots, and the geoid[J]. EPSL, l980, 48:15-22.
[12] Courtillot V, Davaille A, Besse J, et al. Three distinct types of hotspots in the Earth mantle[J]. EPSL, 2003,205:295-308.
[13] Bourdon B, Langmuir C H, Zindler A. Ridge-hotspot interaction along the Mid-Atlantic Ridge between 37°30' and 40°30'N:the U-Th disequilibrium evidence[J]. EPSL, 1996,142(1-2):175-189.
[14] Mara M Yale, Jason P M. Asthenosphere flow model of hotspot-ridge interactions:a comparison of Iceland and Kerguelen[J]. EPSL,1998,161:45-56.
[15] Zhang Y S,Tanimoto T. Ridges, hotspots and their interaction as observed in seismic velocity maps[J]. Nature, 1992,355:44-49.
[16] Ito G, Lin J. Oceanic spreading center-hotspot interactions:constraints from along-isochron bathymetric and gravity anomalies[J]. Geology,1995,23(7):657-660.
[17] Small C. Observations of ridge-hotspot interactions in the Southern Ocean[J].J. Geophys. Res., 1995,100:17931-17946.
[18] Mittelstaedt E, Ito G. Hotspot-ridge interaction:shaping the geometry of Mid-Ocean Ridges[J]. AGU, Fall Meeting 2001,abstract.
[19] Freighner M A, Mark A R. The fluid dynamics of plume-ridge and plume-plate interactions:An experimental investigation[J]. EPSL,1995a,129:171-182.
[20] Ribe N M, Christensen U R, Thei βing J. The dynamics of plume-ridge interaction, 1:Ridge-centered plumes[J]. EPSL,1995,134:155-168.
[21] Ito G, Lin J, Gable C W. Dynamics of mantle flow and melting at a ridge-centered hotspot:Iceland and the Mid-Atlantic Ridge[J]. EPSL, 1996,144:53-74.
[22] Kincaid C, Ito G, Gable C. Laboratory investigation of the interaction of off-axis mantle plumes and spreading centers[J]. Nature, 1995,376:758-761.
[23] Ribe N M. The dynamics of plume-ridge interaction, 2, Off-ridge plumes[J]. J. Geophys. Res., 1996, 101:16195-16204.
[24] Ito G, Lin J,Gable C W. Interaction of mantle plumes and migrating mid-ocean ridges:Implications for the Galapagos pume-ridge system[J]. J. Geophys. Res.,1997,102(B7):15403-15417.
[25] Mittelstaedt E, Ito G. Plume-ridge interaction, lithospheric stresses, and the origin of near-ridge volcanic lineaments[J]. Geochem. Geophys. Geosyst., 2005, 6(6):1-20.
[26] Rhodes J M,Morgan C,Liias R A. Geochemistry of axial seamount lavas:Magmatic relationship between the Cobb hotspot and the Juan de Fuca ridge[J].J. Geophys. Res., 1990,95(B8):12713-12733.
[27] West M, Menke William, Tolstoy M. Focused magma supply at the intersection of the Cobb hotspot and Juan de Fuca ridge[J].Geophys. Res. Lett., 2003,30(14):1724-1727.
[28] Devey C W, Hekinian R, Ackermand D. The foundation seamount chain:a first survey and sampling[J]. Marine Geology, 1997,137:191-200.
[29] Hekinian R, Stoffers P,Ackermand D,et al.Ridge hotspot interaction:the Pacific Antarctic Ridge and the foundation seamounts[J]. Marine Geology, 1999,160:199-223.
[30] Maia M, Ackermandb D, Dehghanic G A,et al.The Pacific-Antarctic Ridge Foundation hotspot interaction:a case study of a ridge approaching a hotspot[J]. Marine Geology,2000,167:61-84.
[31] Canales J P, Ito G, Detrick R S,et al. Crustal thickness along the western Galapagos spreading center and the compensation of the Galapagos hotspot swell[J].EPSL, 2002, 203:311-327.
[32] Hoernle K, Werner R, Morgan J P. Existence of complex spatial zonation in the Gal pagos plume[J]. Geology, 2000,28(5):435-438.
[33] Verma S P, Schilling J G, Waggoner D G. Neodymium isotopic evidence for Galapagos hotspot-spreading centre system evolution[J]. Nature,1983,306:654-657.
[34] Cushman B, Sinton J, Ito G, et al. Glass compositions, plume-ridge interaction, and hydrous melting along the Gal pagos spreading center, 90.5°W to 98°W[J]. Geochem. Geophys. Geosyst., 2004, 5:Q089E17, doi: 10.1029/2004GC000709.
[35] Haase K M. Geochemical constraints on magma sources and mixing processes in Easter Microplate MORB (SE Pacific):a case study of plume ridge interaction[J]. Chemical Geology, 2002, 182:335-355.
[36] Richard M A, Guust N, Morgan W J. Plume-driven plumbing and crustal formation in Iceland[J]. J. Geophys. Res., 2002,107(B8)abstract.
[37] Schilling J G,Meyer P S, Kingsley R H. Evolution of the Iceland hotspot[J].Nature,1982,296:313-320.
[38] Thirlwall M F. Generation of the Pb isotopic characteristics of the Iceland plume[J]. J. Geological Soc., 1995, 152(6):991-996.
[39] Cannat M, Briais A, Deplus C. Mid -Atlantic Ridge-Azores hotspot interactions:along-axis migration of a hotspot-derived event of enhanced magmatism 10 to 4 Ma ago[J]. EPSL, 1999, 173:257-269.
[40] Schilling J G. Azores mantle blob:the rare earth evidence[J]. EPSL,1975, 25:103-105.
[41] Yu D Y, Fontignie D, Schilling J S. Mantle plume-ridge interactions in the Central North Atlantic——A Nd isotope study of Mid-Atlantic Ridge basalts from 30°N to 50°N[J]. EPSL, 1997, 146:259-272.
[42] Schilling J G, Thompson G, Kingsley R,et al.Hotspot-migrating ridge interaction in the South Atlantic[J]. Nature, 1985,313:187-191.
[43] Hanan B B, Kingsley R H,Schilling J G. Pb isotope evidence in the South Atlantic for migrating ridge-hotspot interactions[J]. Nature,1986,322:137-144.
[44] Fontignie D, Schilling J G. Mantle heterogeneities beneath the South Atlantic:A Nd-Sr-Pb isotope study along the Mid-Atlantic Ridge (3°S~46°S)[J]. EPSL, 1996,142:209-221.
[45] Dyment J, Hemond C,Guillou H. Central Indian Ridge and Reunion Hotspot in Rodrigues Area:Another Type of Ridge-Hotspot Interaction?[J]. AGU, Fall Meeting 2001, abstract.
[46] Nauret F, Hemond C, Hofmann A W, et al. Hotspot-ridge interaction between 18°30'S and 20°S along the Central Indian Ridge[J]. AGU, Fall Meeting 2001, abstract.
[47] Georgen J E, Lin J, Dick H J B. Evidence from gravity anomalies for interactions of the Marion and Bouvet hotspots with the Southwest Indian Ridge:effects of transform offsets[J]. EPSL,2001,187:283-300.
[48] Graham D W, K Johnson T M, Priebe L D. Hotspot-ridge interaction along the Southeast India Ridge near Amsterdam and St. Paul islands:helium isotope evidence[J]. EPSL, 1999,167:297-310.
[49] Nakanishi M, Tamaki K, Fujimoto H,et al.Topographic expression of the Gulf of Aden spreading system and its tectonics:hotspots-ridge interaction in Gulf of Aden[J]. AGU, Fall Meeting 2001,abstract.
[50] Tamaki K,Fujimoto H, Orihashi Y, et al. Hotspot-ridge interaction in the Gulf of Aden, East Africa[J]. AGU, Fall Meeting 2001, abstract.
[51] Niu Y L.On the great plume debate[J]. Chinese Science Bulletin, 2005,50(15):1537-1540.
计量
- 文章访问数: 1968
- HTML全文浏览量: 177
- PDF下载量: 26