Rock physics calculation and pre-stack inversion of shear modulus parameters for natural gas hydrate-bearing sediments
-
摘要:
利用地震资料进行天然气水合物识别,需要解决两个重要问题,一是建立含天然气水合物地层的弹性参数与水合物物性参数之间的定量关系;二是建立地震叠前数据与弹性参数之间的定量关系。针对上述问题,本文从神狐海域天然气水合物的微观赋存形态出发,重点考虑水合物的非零剪切模量特征,利用非耦合微分有效介质(DEM)理论、斑块饱和理论进行岩石物理建模,建立了天然气水合物饱和度与弹性参数之间的非线性关系,并通过弹性参数的敏感性分析以及交会分析,确定了剪切模量对水合物含量有较好的指示作用;使用叠前贝叶斯地震反演方法,从叠前地震数据中提取水合物剪切模量参数,结果与SH2站位测井数据有较好的一致性,有效地进行了天然气水合物识别。本方法可为海域天然气水合物识别和预测提供有力的技术支持。
Abstract:To accurately estimate gas hydrate content using seismic data, two key problems are needed to be solved: the first is to establish the quantitative relationship between the elastic parameters and physical parameters of hydrate-bearing sediments; the second is to establish the quantitative relationship between seismic pre-stack data and the elastic parameters. In view of the above problems, we combined the microscopic structure of gas hydrate-bearing sediments in Shenhu area on the focus of its characteristics of non-zero shear modulus, used the uncoupled differential effective medium (DEM) theory and patchy saturation theory to carry out the modeling of rock physics, and established the nonlinear relationship between gas hydrate saturation and elastic parameters. Through the analysis of sensitivity and crossplots of elastic parameters, we determined that the shear modulus is a good indicator to the hydrate content. Using the pre-stack Bayesian seismic inversion method, the shear modulus of hydrate-bearing sediments were extracted from the pre-stack seismic data, and the natural gas hydrate was effectively identified. The method in this study can provide a technical support for the identification and prediction of natural gas hydrate in sea areas.
-
Keywords:
- gas hydrate /
- rock physics /
- shear modulus /
- seismic inversion
-
西湖凹陷是东海陆架盆地中最大的富烃凹陷,凹陷内部充填新生代地层,最大厚度可达17 km左右,已成为中国近海重要的富含油气单元[1-3]。自勘探以来,西湖凹陷已有8个油气田与4个含油气构造发现,主要分布于西部斜坡带与中央反转构造带。孔雀亭地区位于西湖凹陷西部斜坡带中部,该区域7口探井中,有6口可见工业性油流,其中天然气探明储量约为82亿m3,凝析油探明储量约247万m3,多种迹象表明,孔雀亭地区为西湖凹陷内部具有巨大油气勘探潜力的区块。
经过近50年的勘探,前人已通过地震资料、测井资料、岩屑录井资料、取芯资料等对西湖凹陷的沉积环境进行了诸多解释,认为始新统宝石组、平湖组为三角洲—半封闭海湾环境,渐新统花港组与中新统龙井组为河流—湖泊环境。目前关于西湖凹陷的沉积相刻画多集中于凹陷西南部与南部[4-7],然而针对孔雀亭地区的沉积相解释与控制因素的相关研究仍然很少。与此同时,对西湖凹陷已开展的古地貌研究主要是基于地层印模法,尚未进行去压实、古水深校正,因此西湖凹陷古地貌研究也不够深入。基于上述背景,本文以西湖凹陷孔雀亭地区平湖组沉积为研究对象,精细解剖平湖组沉积微相,并系统分析全球海平面、区域构造活动、古气候、古地貌对沉积微相的控制作用,其中古地貌重建过程考虑了去压实与古水深校正。
1. 区域地质概况
东海陆架盆地为中国近海重要的弧后扩张型盆地,主要受始新世—渐新世时期菲律宾板块向欧亚板块俯冲过程影响(图1)。古新世时期,弧后裂陷作用使得东海陆架盆地内产生了隆洼相间的地形格局,而西湖凹陷则是最大的负地形单元,所占面积约为45000 km2[7]。自始新世以来,西湖凹陷经历了裂陷期、断拗转换期、构造反转期、区域热沉降期[8]。平面上,西湖凹陷以西为海礁隆起与渔山东低隆起,以东为陆架外缘。西湖凹陷内部具有“东西分带、南北分块”的构造格局,沿着东西方向可以将西湖凹陷划分出西侧斜坡带、中部反转构造带、东部断阶带[9-12]。平湖斜坡带位于西湖凹陷西侧,是极富油气的构造单元之一[13-14]。
西湖凹陷内部充填了巨厚古近系、新近系、第四系沉积,包括古新统、始新统(宝石组、平湖组)、渐新统(花港组)、中新统(龙井组、玉泉组、柳浪组)、上新统(三潭组)、更新统(东海组),沉积环境由早期的滨浅海,经历半封闭海湾、三角洲-潮坪、河流、海陆过渡,最终演变为现今的浅海环境[11]。针对平湖组地层,多数学者认为西湖凹陷为受东侧钓鱼岛隆起带所影响的半封闭海湾[4,15-16]。然而,对西湖凹陷不同位置平湖组地层开展的研究结果表明,半封闭海湾内的潮坪与三角洲沉积在不同区域发育的规模与期次存在显著差异性[2, 7, 15, 17-19],这也反映了西湖凹陷内复杂的构造演变史与古地貌特征对局部的沉积相发育及其展布特征具有明显的影响。
2. 平湖组沉积相类型
综合岩芯观察、测井相分析,对孔雀亭地区平湖组沉积相进行了精细刻画与解释,识别出辫状河三角洲相与潮坪相,其中辫状河三角洲相主要包含辫状河三角洲前缘亚相(水下分流河道、分流间湾、河口坝、水下天然堤、席状砂),潮坪相主要包含潮间坪亚相(沙坪、泥坪、混合坪),详细信息见表1。
2.1 辫状河三角洲相
孔雀亭地区平湖组发育的辫状河三角洲相主要包括辫状河三角洲前缘亚相,其可进一步细分为水下分流河道、分流间湾、天然堤、河口坝、席状砂等微相。
(1)水下分流河道微相。岩芯上常见底部冲刷面,其上发育含砾滞留沉积,向上发育含砂砾岩、细砂岩、粉砂岩、粉砂质泥岩以及泥质粉砂岩,整体呈现出向上变细的正韵律特征。楔状交错层理、板状交错层理常见。测井曲线上,伽马曲线与电阻率曲线呈现出钟形特征或者箱型特征。钟形特征体现了河道内沉积物粒度逐渐变小的规律,而箱型特征则体现了河道内以相对较粗的沉积物为主,其为晚期河道切割早期水道从而产生多期水道叠置的结果。该微相主要发育于平中段与平上段,单个水道厚度范围为10~30 m,多期水道厚度累计可达85 m。
(2)分流间湾微相。是发育于分流河道之间以及三角洲朵体之间的沉积微相。岩芯上表现为灰色粉砂质泥岩、泥质粉砂岩,水平层理常见。此外,可见生物钻孔构造。测井曲线上,自然电位曲线整体很低,而自然伽马整体很高。该微相主要发育于平中段与平上段。
(3)水下天然堤微相。是发育于水下分流河道两侧的微相类型。岩芯上可见典型的爬升层理,主要由粉砂岩组成。自然伽马值整体较高,电阻率曲线较为平直。该微相主要与水下分流河道伴生,岩芯内常常位于河道沉积之上,代表河道侧向迁移后晚期天然堤覆盖早期河道砂体。岩芯上该微相厚度常为3~5 m。
(4)河口坝微相。岩芯上主要以灰色细砂岩、粉砂岩为主,分选较好,重要的特征为下细上粗的反韵律特征。交错层理与平行层理常见。自然伽马曲线与电阻率曲线均表现为漏斗形特征,同样反映了向上变粗的沉积递变序列。该微相发育于平中段与平上段,砂体厚度常在10 m左右。
(5)席状砂微相。该微相是由河口坝经过海水波浪作用改造后再次发生的沉积。研究区岩芯中席状砂主要由粉砂岩与泥质粉砂岩组成,其厚度较薄,为3~5 m。沉积构造可见浪成交错层理。自然伽马曲线与电阻率曲线表现为尖峰凸起形态。该微相发育于平中段与平上段研究区靠近东部区域。
2.2 潮坪相
孔雀亭地区发育的潮坪相主要包括潮间坪亚相,其可进一步细分为沙坪、泥坪、混合坪等微相。
(1)沙坪。岩性主要以粉砂岩与细砂岩为主。沉积构造常见交错层理、脉状层理。粒度分析显示,跳跃部分包括两个次总体,反映了沉积物经历了冲刷回流过程。自然伽马曲线整体呈现中间低、两端高的特征,自然电位曲线则呈现中间高、两端低的特征。该微相主要发育于平下下段与平下上段研究区向海一侧。
(2)泥坪。岩性主要为泥岩或者粉砂质泥岩。可见透镜状层理。自然伽马呈现高值,电阻率呈现低值。该微相发育相较于沙坪微相更少,分布面积更小。
(3)混合坪。岩性上主要为薄层粉砂岩与薄层泥岩频繁互层。沉积构造可见波状层理、双黏土层构造。自然伽马与电阻率曲线上呈现出剧烈变化的锯齿状特征。混合坪是潮坪中发育最广泛的一种类型,是平下下段与平下上段常见的沉积微相。
3. 平湖组沉积相空间展布特征
根据研究区内各单井沉积微相精细划分结果建立连井剖面,从而观察沉积微相在井间的变化样式。顺着物源方向的连井剖面(即北西-南东方向,图2a)显示,辫状河三角洲在隆起周缘最为发育,在C井附近的整个平湖组均为辫状河三角洲沉积。向海方向(即B井与E井附近),平下下段与平下上段的辫状河三角洲沉积快速减薄,呈现出潮坪沉积特征,以混合坪沉积为主。平中段与平上段在B井与E井均显示为厚层三角洲沉积,显示出该阶段的辫状河三角洲顺着物源向前快速进积。此外,B井与E井位置的辫状河三角洲沉积比C井位置的辫状河三角洲沉积更厚,反映了洼陷区比隆起区具有更大的可容纳空间,能够堆积更多的三角洲沉积物。垂直于物源方向的连井剖面(即南西-北东向,图2b)显示,平下下段与平下上段发育厚层潮坪沉积,主要为混合坪微相(岩性上表现为薄层粉砂岩与泥岩互层),沙坪微相次之,泥坪微相最少。平中段与平上段发育厚层辫状河三角洲沉积。三角洲前缘的水下分流河道频繁出现,局部可见多期水道垂向叠置。席状砂沉积在不同钻井之上均发育较多。相比之下,河口坝沉积相对较少,可能反映了平中段与平上段时期波浪作用对河口坝沉积体具有明显的改造作用,对其砂体进行进一步淘洗,发生再沉积作用,形成了席状砂。此外,平上段水下分流河道所占比例比平中段比例更大,反映了三角洲向前进积、发育规模逐渐增大的趋势。综上,顺物源方向的剖面与垂直于物源方向的剖面均显示,平湖组由早期的“以潮坪为主”的沉积相向上演变为“以辫状河三角洲为主”的沉积相。
根据地震资料的分析,可以更加全面地识别出沉积相在剖面以及平面上的展布特征(图3)。经过井震资料对比以及地震轴连续追踪,识别出研究区平湖组重要的层段界面,包括T34、T33、T32、T31、T30,其中T34与T33之间为平下下段,T33与T32之间为平下上段,T32与T31之间为平湖组中段,T31与T30之间为平湖组上段。通过分析地震同相轴的终止关系,识别出平湖组发育朝着东南方向的低角度“S”型前积结构与高角度斜交型前积结构,其中以低角度“S”型前积结构为主,地震同相轴下超于地层界面之上,其反映的是在水体相对较浅、沉积物供给充足的情况下发育的三角洲进积作用。结合钻井岩芯与测井资料分析,认为在研究区西北侧发育大规模的辫状河三角洲沉积(图3)。向海方向,地震反射同相轴开始表现为平行结构、亚平行结构,与此同时,地震同相轴的连续性有所增强。结合钻井岩芯与测井资料分析,认为研究区东南侧发育的是潮坪沉积(图3)。在平下下段时期,前积结构主要出现在隆起单元附近,向海方向推进距离有限,反映了该阶段辫状河三角洲沉积主要在距离源区不远的范围内有所发育,而潮坪沉积才是该时期主要的沉积相类型。进入平下上段时期后,前积结构向研究区东南方向有一定程度扩张,反映了辫状河三角洲沉积面积向海方向发生一定程度的增大。在平中段与平上段,前积结构向海方向快速推进,反映了该阶段辫状河三角洲沉积向海推进较远,占据研究区大部分范围,而潮坪沉积则是向海方向大范围退缩。因此,地震剖面也显示出研究区内由平下下段、平下上段至平中段与平上段辫状河三角洲沉积规模逐渐增大,而潮坪沉积规模逐渐减小(图3b)。
为了查明沉积微相在研究区平面上的展布特征,本文采用层间均方根振幅进行属性提取(图4)。平下下段在靠近隆起单元西部与中部的区域呈现出高振幅区,结合位于相同区域的钻井资料(A井、D3井、D2井),认为该区主要发育辫状河三角洲沉积。远离隆起单元往东南方向,振幅呈现出中-低值,结合位于同一区域的钻井资料(如D4井、G井),认为靠近海洋一侧主要发育潮坪沉积,其中振幅相对较高值区(亮点区域)为沙坪沉积,振幅相对较低值区为泥坪沉积,介于两者之间则为混合坪沉积。在研究区东侧与南侧,振幅属性呈现出最低值,其反映了泥坪沉积。与平下下段相比,平下上段中高振幅区具有更大的分布范围,反映了辫状河三角洲沉积具有更大的展布面积。与此同时,研究区南部出现条带状中等振幅,从而反映了平下上段研究区南部出现了较多的沙坪沉积。该段弱振幅比平下下段区域更小,反映了泥坪沉积的减少。平中段时期,均方根振幅高值区进一步南移,代表富砂的辫状河三角洲沉积进一步扩大,潮坪相有所退却。研究区正南部存在中等—强振幅,代表沙坪沉积,而在东部存在弱振幅,代表泥坪沉积。进入平上段之后,研究区的均方根振幅值有所降低。该层段的钻井资料显示发育厚层砂质沉积,是由于砂质沉积之间的非均质性降低导致了均方根振幅有所下降。结合地震属性与钻井资料,平上段三角洲比平中段三角洲沉积厚度更大,面积更广,特别是研究区西侧辫状河三角洲具有比平中段更大的朵体面积。
综合分析钻井岩心、测井曲线、地震剖面、地震属性,获得不同层段沉积相类型分布图(图5)。平下下段辫状河三角洲沉积局限于研究区西北部,面积较小,水下分流河道向前推进程度有限。该阶段研究区中部及东部发育大片潮坪沉积,以混合坪沉积为主,沙坪沉积分布在研究区南部与中部,而泥坪沉积在研究区东部、北部、南部均有发育。平下上段辫状河三角洲沉积范围有所增大,水下分流河道向前推进程度增加。混合坪沉积仍然是主要的潮坪沉积类型,沙坪沉积相对集中分布于研究区南部,而泥坪沉积分布范围有所减小。平中段时期辫状河水下分流河道进一步向东南方向延展,三角洲分布面积进一步扩大。混合坪、沙坪、泥坪沉积均有所较少。平上段时期,辫状河三角洲分布面积达到最大,而潮坪沉积只在研究区东南侧有所发育。沙坪沉积分布于南部与东部,而泥坪沉积主要分布于东部。
4. 沉积相控制因素
西湖凹陷孔雀亭地区平湖组沉积相在平下下段时期以潮坪为主、辫状河三角洲为辅,随着时间的推移,潮坪沉积分布面积逐渐降低,而辫状河三角洲沉积分布面积逐渐增高。沉积相的迁移往往和全球海平面与区域构造、古气候、古地貌有关[20-29]。
4.1 全球海平面与区域构造
全球海平面上升和区域构造沉降均会使得相对海平面上升,从而造成滨线向陆方向迁移,影响沉积相带的分布与发育[30-32]。根据前人已有的年代框架,认为平湖组发育在始新世,其中平湖组下段的底界面为41 Ma,平湖组中段的底界面为39 Ma,平湖组上段的底界面为37.2 Ma,平湖组的顶为33.9 Ma[19, 33]。根据Miller重建的全球海平面变化曲线,全球海平面在41~40 Ma间逐渐升高,40~37 Ma全球海平面逐渐降低,37~33.9 Ma全球海平面重新出现上升趋势[34](图6)。区域构造方面,由于平湖组时期西湖凹陷处于断陷沉降阶段,因而可以通过计算、统计研究区的沉降速率,来评价区域内构造下沉快慢程度。通过计算研究区内391条断层活动速率,获得了平下段、平中段、平上段各时期内沉降速率值。平下段时期,区域沉降速率最大值为100.5 m/Ma,平均值为19 m/Ma;平中段时期,区域沉降速率最大值为66.7 m/Ma,平均值为7.5 m/Ma;平上段时期,区域沉降速率最大值为42.7 m/Ma,平均值为3.7 m/Ma(图6)。因此,由平下段-平中段至平上段,构造沉降速率快速降低,这也符合西湖凹陷在平湖组时期处于断陷期至断陷晚期的背景[11]。由于目前尚未对平下下段与平下上段界面进行定年工作,因而无法对平下段内部构造沉降历史进行进一步细分,但根据区域背景演化趋势,有理由认为平下下段时期的沉降速率大于平下上段的沉降速率。结合全球海平面与区域构造信息,研究区内的相对海平面在平下段时期最高,在平中段时期有所下降,至平上段时期再次出现上升。研究区内,平下段辫状河三角洲只出现在隆起单元附近,潮坪沉积呈现主体地位,反映了相对海平面对沉积相发育过程的控制。高的相对海平面使得辫状河三角洲只在入水口附近有所发育,向海方向推进程度有限,而在该环境背景下,潮坪环境可大量发育。平下段内部,由平下下段至平下上段,辫状河三角洲出现轻度扩张,而潮坪沉积存在轻微向海方向退却,反映了从平下下段过渡至平下上段时期,构造沉降速率逐渐减小,导致相对海平面缓慢下降。在平中段时期,由于相对海平面快速下降,为辫状河三角洲的向海推进提供了极佳条件,因此平中段时期的辫状河三角洲比平下段时期的辫状河三角洲存在明显的扩张。由平中段至平上段,相对海平面再次上升,在该背景下,理应出现辫状河三角洲向陆方向退积,而事实上,平上段的辫状河三角洲比平中段的辫状河三角洲有所扩张,而潮坪沉积向海方向退却,因此,这反映了除全球海平面与区域构造之外,还存在其他因素对平湖组的沉积相存在控制作用。
4.2 古气候
古代气候条件可以通过温度、降雨以及其带来的风化作用来对大陆沉积进行改造,进而通过河流搬运输入海洋[35-36]。新生代CO2浓度指标重建结果显示,37 Ma左右是新生代CO2浓度最高的时间[34](图6)。前人已有的研究显示,高的CO2浓度会使得气候温度升高,降雨增多,从而加速大陆隆起区的物理风化,产生大量的硅质陆源碎屑,并在河流输送量充足的条件下,将大量陆源碎屑沉积物搬运至海洋沉积下来。因此,平中段与平上段时期,正是处于CO2浓度最高值范围,从而当时高的温度、降雨过程使得物源区范围增大,剥蚀强度增大,大量的沉积物经河流搬运作用输入至西湖凹陷,使得平中段与平上段时期的辫状河三角洲得到大范围扩张,向海推进,而潮坪沉积则向海方向发生退却。
4.3 古地貌
古地貌(即古代沉积环境中的隆起与洼陷格局)对沉积相的发育与分布具有重要的控制作用[37-39]。本文利用BacKang软件,基于逐层回剥的思想,通过输入层位深度、地层物性组成、原始孔隙度、水体密度、沉积物密度、压实系数等参数,进行去压实处理与古水深校正,进而获取研究区高精度的微地貌特征(图7)。由于研究区钻井较少,不足以完成古地貌重建,因此本文首先将三维地震数据体从时间域转化为深度域,进而将三维数据设为多口虚拟井,对每口虚拟井进行一系列的数据处理。文章中共使用约170万口虚拟井进行了计算,时间-深度转换公式采用的是平北地区构建的三维速度模型。
古地貌重建结果表明,平下下段古地貌以断陷与低凸起间互出现为特征。研究区西北部呈现带状隆起单元。中部主要为小型断陷带,呈东西方向分布,断陷中间出现局部低凸起,将两侧断陷分隔开。东南部主要为水下低凸起,正南部与正东部为深洼带。平下下段的隆起带周缘为辫状河三角洲分布区,反映了在水体较浅的环境中,河流输入的砂体由于卸载作用而逐渐沉积。在三角洲朵体之间,可见块状分布出现的分流间湾沉积。值得注意的是,该阶段的辫状河三角洲分布面积较小,反映了隆凹相间的古地貌特征限制了辫状河三角洲向海迁移。在平下下段时期的研究区中部以及东南部,由于水深较大以及距离西北部源区较远,潮坪沉积广泛发育,以混合坪沉积为主,沙坪与泥坪沉积相对较少。
平下上段古地貌与平下下段古地貌相比,隆起单元附近的浅水区更加平缓,反映了在经历平下下段辫状河三角洲充填之后,该区域坡度缓、延伸远,为平下上段辫状河三角洲沉积的扩张提供了优势条件。深洼带位于研究区正中央与正南方。沙坪沉积主要出现于DA井南部,从洼陷带至低凸起斜坡之上均有发育。
平中段时期,地貌低势区主要位于研究区东部,而西部区域除了部分断陷出现,大部分以较平缓的斜坡带为特征。平中段时期平缓的地势非常有利于辫状河三角洲的发育,结合该时期有所增强的大陆风化后输入的大量沉积物,辫状河三角洲分布面积快速增加。与此同时,局部洼陷带对砂体具有阻挡作用。DA井西南侧的断层阻碍了辫状河三角洲进一步向南扩张。G井东部的断裂使得其周缘出现深洼带,该深洼带的出现有效捕获了大量沉积砂体,也使得三角洲沉积不再往东南方向推进。远离西北部隆起单元的低凸起区之上主要分布潮坪沉积,以混合坪与沙坪沉积为主。
平上段时期,研究区地势进一步变得平缓,主要是由于平中段时期大量辫状河沉积物的输入,填充了平中段时期的洼陷带。平上段时期平缓的地势与充足的沉积物供给使得辫状河三角洲进一步向前推进,分布面积比平中段时期的辫状河三角洲更大,潮坪沉积只在研究区东南侧小范围发育。
5. 结论
(1)西湖凹陷孔雀亭地区平湖组主要发育辫状河三角洲沉积与潮坪沉积,其中辫状河沉积主要发育辫状河前缘亚相,包括水下分流河道、水下天然堤、河口坝、席状砂、分流间湾等微相,潮坪沉积主要发育潮间坪亚相,主要包括混合坪、沙坪、泥坪等微相。
(2)由平下下段、平下上段至平中段、平上段,辫状河三角洲沉积规模逐渐增大,往东南方向有所扩张,而潮坪沉积规模逐渐减小,往海洋方向退却收缩。
(3)平湖组沉积受相对海平面(全球海平面与区域构造)、古气候、古地貌的联合控制。由平下段至平中段相对海平面下降,使得辫状河三角洲沉积向前进积,潮坪沉积向海后退。在平中段、平上段时期,由于CO2浓度极高,使得大陆风化快速增强,大量的沉积物由河流输送至海洋发生堆积。平下段时期古地貌断陷发育较多,对辫状河三角洲的发育存在限制作用,而在经历了平中段大量沉积物填平补齐后,研究区在平上段时期已经较为平缓,从而有利于三角洲往海方向推进。
-
图 3 实验室获取的水合物CT图像[42]及水合物等效介质模型
a:原始图像,b:对a图红色框线内区域进行灰度分析、着色后放大图,c:等效介质模型。
Figure 3. CT images of gas hydrate obtained in laboratory and equivalent medium model of gas hydrate-bearing sediments
a: Original image, b: enlarged image of the area in the red frame after grey analysis and coloring, c: equivalent medium model of gas hydrate-bearing sediments.
图 8 SH2站位水合物地层的水合物饱和度及弹性参数计算结果
a:根据实测电阻率曲线计算的水合物饱和度,b:纵波速度(红色曲线表示实测值,黑色曲线表示计算值),c:横波速度计算值;d:剪切模量计算值。
Figure 8. Calculation results of hydrate saturation and elastic parameters of hydrate strata at SH2 site
a: hydrate saturation calculated according to the measured resistivity, b: P-wave velocity (red curve represents the measured value, black curve represents the calculated value), c: calculated shear wave velocity; d: calculated shear modulus.
图 10 SH2站位水合物地层弹性参数交会图
蓝色曲线为SH2井的实测纵波速度曲线。a:中角度叠加地震剖面,中心角度为20°;b:中角度弹性阻抗剖面,中心角度为20°;c:反演得到的剪切模量剖面。
Figure 10. Seismic, elastic impedance, and inverted shear modulus profiles of SH2 site
The blue curve is the measured P-wave velocity of SH2 site. a: Partial stack seismic profile at centric angle of 20°; b: elastic impedance profile at centric angle of 20°; c: inverted shear modulus profile.
表 1 数值计算所选用的弹性常数
Table 1 The elastic constants used for numerical calculation
组分 体积模量/GPa 剪切模量/GPa 密度/(g/cm3) 来源 石英 37.0 44.0 2.65 Carmichale, 1989 黏土 21.0 7.0 2.60 Tosaya and Nur, 1982 水合物 7.7 3.2 0.91 Waite et al, 2000 水 2.25 1000 -
[1] Miller J J, Lee M W, von Huene R. An analysis of a seismic reflection from the base of a gas hydrate zone, offshore Peru[J]. AAPG Bulletin, 1991, 75(5):910-924.
[2] 宋海斌, 吴时国, 江为为. 南海东北部973剖面BSR及其热流特征[J]. 地球物理学报, 2007, 50(5):1508-1517 doi: 10.3321/j.issn:0001-5733.2007.05.027 SONG Haibin, WU Shiguo, JIANG Weiwei. The characteristics of BSRs and their derived heat flow on the profile 973 in the northeastern South China Sea[J]. Chinese Journal of Geophysics, 2007, 50(5):1508-1517.] doi: 10.3321/j.issn:0001-5733.2007.05.027
[3] 徐华宁, 陆敬安, 梁金强. 珠江口盆地东部海域近海底天然气水合物地震识别及地质成因[J]. 地学前缘, 2017, 24(4):57-65 XU Huaning, LU Jing’an, LIANG Jinqiang. Seismic identification and geological origin of gas hydrate in near seafloor sediments in the eastern part of the Pearl River Mouth Basin[J]. Earth Science Frontiers, 2017, 24(4):57-65.]
[4] 杨胜雄, 梁金强, 陆敬安, 等. 南海北部神狐海域天然气水合物成藏特征及主控因素新认识[J]. 地学前缘, 2017, 24(4):1-14 YANG Shengxiong, LIANG Jinqiang, LU Jing’an, et al. New understandings on the characteristics and controlling factors of gas hydrate reservoirs in the Shenhu area on the northern slope of the South China Sea[J]. Earth Science Frontiers, 2017, 24(4):1-14.]
[5] 梁劲, 王静丽, 杨承志, 等. 珠江口盆地东部海域含天然气水合物沉积层的地球物理特征[J]. 天然气工业, 2017, 37(2):126-133 doi: 10.3787/j.issn.1000-0976.2017.02.017 LIANG Jin, WANG Jingli, YANG Chengzhi, et al. Geophysical characteristics of gas hydrate bearing sediments in the eastern sea area of the Pearl River Mouth Basin[J]. Natural Gas Industry, 2017, 37(2):126-133.] doi: 10.3787/j.issn.1000-0976.2017.02.017
[6] 张伟, 梁金强, 苏丕波, 等. 南海北部陆坡高饱和度天然气水合物气源运聚通道控藏作用[J]. 中国地质, 2018, 45(1):1-14 doi: 10.12029/gc20180101 ZHANG Wei, LIANG Jinqiang, SU Pibo, et al. Migrating pathways of hydrocarbons and their controlling effects associated with high saturation gas hydrate in Shenhu area, northern South China Sea[J]. Geology in China, 2018, 45(1):1-14.] doi: 10.12029/gc20180101
[7] Wang X J, Collett T S, Lee M W, et al. Geological controls on the occurrence of gas hydrate from core, downhole log, and seismic data in the Shenhu Area, South China Sea[J]. Marine Geology, 2014, 357:272-292. doi: 10.1016/j.margeo.2014.09.040
[8] 刘昌岭, 孙运宝. 海洋天然气水合物储层特性及其资源量评价方法[J]. 海洋地质与第四纪地质, 2021, 41(5):44-57 LIU Changling, SUN Yunbao. Characteristics of marine gas hydrate reservoir and its resource evaluation methods[J]. Marine Geology & Quaternary Geology, 2021, 41(5):44-57.]
[9] Dvorkin J, Nur A. Elasticity of high-porosity sandstones: theory for two North Sea data sets[J]. Geophysics, 1996, 61(5):1363-1370. doi: 10.1190/1.1444059
[10] Ecker C, Dvorkin J, Nur A. Sediments with gas hydrates: internal structure from seismic AVO[J]. Geophysics, 1998, 63(5):1659-1669. doi: 10.1190/1.1444462
[11] Helgerud M B, Dvorkin J, Nur A, et al. Elastic-wave velocity in marine sediments with gas hydrates: effective medium modeling[J]. Geophysical Research Letters, 1999, 26(13):2021-2024. doi: 10.1029/1999GL900421
[12] Ecker C, Dvorkin J, Nur A M. Estimating the amount of gas hydrate and free gas from marine seismic data[J]. Geophysics, 2000, 65(2):565-573. doi: 10.1190/1.1444752
[13] Lee M W, Collett T S. Pore-and fracture-filling gas hydrate reservoirs in the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II Green Canyon 955 H well[J]. Marine and Petroleum Geology, 2012, 34(1):62-71. doi: 10.1016/j.marpetgeo.2011.08.002
[14] Qian J, Wang X J, Collett T S, et al. Gas hydrate accumulation and saturations estimated from effective medium theory in the eastern Pearl River Mouth Basin, South China Sea[J]. Interpretation, 2017, 5(3):SM33-SM48. doi: 10.1190/INT-2016-0217.1
[15] Terry D A, Knapp C C. A unified effective medium model for gas hydrates in sediments[J]. Geophysics, 2018, 83(6):MR317-MR332. doi: 10.1190/geo2017-0513.1
[16] Winters W J, Waite W F, Mason D H, et al. Methane gas hydrate effect on sediment acoustic and strength properties[J]. Journal of Petroleum Science and Engineering, 2007, 56(1-3):127-135. doi: 10.1016/j.petrol.2006.02.003
[17] Mallick S, Huang X R, Lauve J, et al. Hybrid seismic inversion: a reconnaissance tool for deepwater exploration[J]. The Leading Edge, 2000, 19(11):1230-1237. doi: 10.1190/1.1438512
[18] Lu S M, McMechan G A. Estimation of gas hydrate and free gas saturation, concentration, and distribution from seismic data[J]. Geophysics, 2002, 67(2):582-593. doi: 10.1190/1.1468619
[19] Lu S M, McMechan G A. Elastic impedance inversion of multichannel seismic data from unconsolidated sediments containing gas hydrate and free gas[J]. Geophysics, 2004, 69(1):164-179. doi: 10.1190/1.1649385
[20] Wang X J, Wu S G, Xu N, et al. Estimation of gas hydrate saturation using constrained sparse spike inversion: case study from the Northern South China Sea[J]. Terrestrial Atmospheric and Oceanic Sciences, 2006, 17(4):799-813. doi: 10.3319/TAO.2006.17.4.799(GH)
[21] 王秀娟, 吴时国, 余鹏, 等. 混合地震反演技术及其在水合物勘探中的应用[J]. 石油物探, 2007, 46(3):278-282,301 doi: 10.3969/j.issn.1000-1441.2007.03.009 WANG Xiujuan, WU Shiguo, YU Peng, et al. Hybrid seismic inversion and its application in hydrate exploration[J]. Geophysical Prospecting for Petroleum, 2007, 46(3):278-282,301.] doi: 10.3969/j.issn.1000-1441.2007.03.009
[22] 杨睿, 吴能友, 雷新华, 等. 波阻抗反演在南海北部神狐海域天然气水合物勘探中的应用[J]. 现代地质, 2010, 24(3):495-500 doi: 10.3969/j.issn.1000-8527.2010.03.012 YANG Rui, WU Nengyou, LEI Xinhua, et al. Impedance inversion and its application in gas hydrate exploration in Shenhu Area, Northern South China Sea[J]. Geoscience, 2010, 24(3):495-500.] doi: 10.3969/j.issn.1000-8527.2010.03.012
[23] Wang X J, Wu S G, Guo Y Q, et al. Geophysical indicators of gas hydrate in the Northern Continental Margin, South China Sea[J]. Journal of Geological Research, 2011, 2011:359597.
[24] Chen M A P, Riedel M, Hyndman R D, et al. AVO inversion of BSRs in marine gas hydrate studies[J]. Geophysics, 2007, 72(2):C31-C43. doi: 10.1190/1.2435604
[25] Zhang R W, Li H Q, Zhang B J, et al. Detection of gas hydrate sediments using prestack seismic AVA inversion[J]. Applied Geophysics, 2015, 12(3):453-464. doi: 10.1007/s11770-015-0503-3
[26] Liu T, Liu X W. Identification of morphologies of gas hydrate distribution based on amplitude variation with angle analysis[J]. Geophysics, 2018, 83(3):B143-B154. doi: 10.1190/geo2017-0072.1
[27] 刘斌, 张衡. 海底水合物OBS数据全波形反演数值例子[J]. 地球物理学进展, 2018, 33(1):379-384 doi: 10.6038/pg2018AA0580 LIU Bin, ZHANG Heng. Full waveform inversion of OBS data from ocean gas hydrate: synthetic example[J]. Progress in Geophysics, 2018, 33(1):379-384.] doi: 10.6038/pg2018AA0580
[28] 霍元媛, 杨睿, 潘纪顺, 等. 波形反演在天然气水合物中的应用研究进展[J]. 海洋地质与第四纪地质, 2022, 42(4):207-221 HUO Yuanyuan, YANG Rui, PAN Jishun, et al. Application of full waveform inversion to gas hydrate research[J]. Marine Geology & Quaternary Geology, 2022, 42(4):207-221.]
[29] Holbrook W S, Hoskins H, Wood W T, et al. Methane hydrate and free gas on the Blake Ridge from vertical seismic profiling[J]. Science, 1996, 273(5283):1840-1843. doi: 10.1126/science.273.5283.1840
[30] 刘学伟, 李敏锋, 张聿文, 等. 天然气水合物地震响应研究: 中国南海HD152测线应用实例[J]. 现代地质, 2005, 19(1):33-38 doi: 10.3969/j.issn.1000-8527.2005.01.005 LIU Xuewei, LI Minfeng, ZHANG Yuwen, et al. Studies of seismic characteristics about gas hydrate: a case study of Line HD152 in the South China Sea[J]. Geoscience, 2005, 19(1):33-38.] doi: 10.3969/j.issn.1000-8527.2005.01.005
[31] 印兴耀, 宗兆云, 吴国忱. 岩石物理驱动下地震流体识别研究[J]. 中国科学(地球科学), 2015, 45(1): 8-21 YIN Xingyao, ZONG Zhaoyun, WU Guochen. Research on seismic fluid identification driven by rock physics[J]. Science China Earth Sciences, 2015, 58(2): 159-171.]
[32] 吴能友, 杨胜雄, 王宏斌, 等. 南海北部陆坡神狐海域天然气水合物成藏的流体运移体系[J]. 地球物理学报, 2009, 52(6):1641-1650 doi: 10.3969/j.issn.0001-5733.2009.06.027 WU Nengyou, YANG Shengxiong, WANG Hongbin, et al. Gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, Northern South China Sea[J]. Chinese Journal of Geophysics, 2009, 52(6):1641-1650.] doi: 10.3969/j.issn.0001-5733.2009.06.027
[33] 徐华宁, 杨胜雄, 郑晓东, 等. 南中国海神狐海域天然气水合物地震识别及分布特征[J]. 地球物理学报, 2010, 53(7):1691-1698 doi: 10.3969/j.issn.0001-5733.2010.07.020 XU Huaning, YANG Shengxiong, ZHENG Xiaodong, et al. Seismic identification of gas hydrate and its distribution in Shenhu Area, South China Sea[J]. Chinese Journal of Geophysics, 2010, 53(7):1691-1698.] doi: 10.3969/j.issn.0001-5733.2010.07.020
[34] 吴能友, 张海啟, 杨胜雄, 等. 南海神狐海域天然气水合物成藏系统初探[J]. 天然气工业, 2007, 27(9):1-6 doi: 10.3321/j.issn:1000-0976.2007.09.001 WU Nengyou, ZHANG Haiqi, YANG Shengxiong, et al. Preliminary discussion on natural gas hydrate (NGH) reservoir system of Shenhu Area, North Slope of South China Sea[J]. Natural Gas Industry, 2007, 27(9):1-6.] doi: 10.3321/j.issn:1000-0976.2007.09.001
[35] 陆红锋, 陈弘, 陈芳, 等. 南海神狐海域天然气水合物钻孔沉积物矿物学特征[J]. 南海地质研究, 2009, 20:28-39 LU Hongfeng, CHEN Hong, CHEN Fang, et al. Mineralogy of the sediments from gas-hydrate drilling sites, Shenhu Area, South China Sea[J]. Geological Research of South China Sea, 2009, 20:28-39.]
[36] Wang X J, Hutchinson D R, Wu S G, et al. Elevated gas hydrate saturation within silt and silty clay sediments in the Shenhu area, South China Sea[J]. Journal of Geophysical Research: Solid Earth, 2011, 116(B5):B05102.
[37] Schultheiss P, Holland M, Humphrey G. Wireline coring and analysis under pressure: recent use and future developments of the HYACINTH system[J]. Scientific Drilling, 2009, 7:44-50. doi: 10.5194/sd-7-44-2009
[38] 李承峰, 叶旺全, 陈亮, 等. 南海神狐海域天然气水合物微观赋存特征的超分辨率CT图像识别[J]. 海洋地质与第四纪地质, 2024, 44(3):149-159 LI Chengfeng, YE Wangquan, CHEN Liang, et al. Super-resolution CT image recognition of micro-occurrence characteristics of natural gas hydrates from Shenhu area in northern South China Sea[J]. Marine Geology & Quaternary Geology, 2024, 44(3):149-159.]
[39] Lee M W, Waite W F. Estimating pore-space gas hydrate saturations from well log acoustic data[J]. Geochemistry, Geophysics, Geosystems, 2008, 9(7):Q07008.
[40] Guerin G, Goldberg D, Meltser A. Characterization of in situ elastic properties of gas hydrate-bearing sediments on the Blake Ridge[J]. Journal of Geophysical Research: Solid Earth, 1999, 104(B8):17781-17795. doi: 10.1029/1999JB900127
[41] 刘昌岭, 业渝光, 孟庆国, 等. 南海神狐海域天然气水合物样品的基本特征[J]. 热带海洋学报, 2012, 31(5):1-5 doi: 10.3969/j.issn.1009-5470.2012.05.001 LIU Changling, YE Yuguang, MENG Qingguo, et al. Characteristics of gas hydrate samples recovered from Shenhu Area in the South China Sea[J]. Journal of Tropical Oceanography, 2012, 31(5):1-5.] doi: 10.3969/j.issn.1009-5470.2012.05.001
[42] 胡高伟, 李承峰, 业渝光, 等. 沉积物孔隙空间天然气水合物微观分布观测[J]. 地球物理学报, 2014, 57(5):1675-1682 doi: 10.6038/cjg20140530 HU Gaowei, LI Chengfeng, YE Yuguang, et al. Observation of gas hydrate distribution in sediment pore space[J]. Chinese Journal of Geophysics, 2014, 57(5):1675-1682.] doi: 10.6038/cjg20140530
[43] Berryman J G. Mixture theories for rock properties[M]//Ahrens T J. Rock Physics & Phase Relations: A Handbook of Physical Constants. Washington: American Geophysical Union, 1995: 205-228.
[44] 印兴耀, 化世榜, 宗兆云. 基于线性近似的微分等效介质方程解耦方法[J]. 石油地球物理勘探, 2016, 51(2):281-287 YIN Xingyao, HUA Shibang, ZONG Zhaoyun. A decoupling approach for differential equivalent equations based on linear approximation[J]. Oil Geophysical Prospecting, 2016, 51(2):281-287.]
[45] Xu S Y, White R E. A new velocity model for clay-sand mixtures[J]. Geophysical Prospecting, 1995, 43(1):91-118. doi: 10.1111/j.1365-2478.1995.tb00126.x
[46] Ciz R, Shapiro S A. Generalization of Gassmann equations for porous media saturated with a solid material[J]. Geophysics, 2007, 72(6):A75-A79. doi: 10.1190/1.2772400
[47] Hill R. Elastic properties of reinforced solids: some theoretical principles[J]. Journal of the Mechanics and Physics of Solids, 1963, 11(5):357-372. doi: 10.1016/0022-5096(63)90036-X
[48] Reuss A. Calculation of the flow limits of mixed crystals on the basis of the plasticity of monocrystals[J]. Zeitschrift fur Angewandte Mathematik und Mechanik, 1929, 9(1):49-58. doi: 10.1002/zamm.19290090104
[49] 杨晓璐, 钟思玲, 万志峰. 泥底辟/泥火山流体热效应及其对天然气水合物赋存的影响[J]. 海洋地质前沿, 2018, 34(7):15-23 YANG Xiaolu, ZHONG Siling, WAN Zhifeng. The thermodynamics of mud diapir/volcano fluid and its influence on gas hydrate occurrence[J]. Marine Geology Frontiers, 2018, 34(7):15-23.]
[50] 宗兆云, 印兴耀, 吴国忱. 基于叠前地震纵横波模量直接反演的流体检测方法[J]. 地球物理学报, 2012, 55(1):284-292 doi: 10.6038/j.issn.0001-5733.2012.01.028 ZONG Zhaoyun, YIN Xingyao, WU Guochen. Fluid identification method based on compressional and shear modulus direct inversion[J]. Chinese Journal of Geophysics, 2012, 55(1):284-292.] doi: 10.6038/j.issn.0001-5733.2012.01.028
-
期刊类型引用(9)
1. 陈生华,王健伟,刘舒,严曙梅,韩建辉,傅恒,谢才铸,孙莉. 东海西湖凹陷孔雀亭地区始新统平湖组中段沉积特征. 岩性油气藏. 2025(02): 103-114 . 百度学术
2. 郭刚,苏圣民,徐建永,刘志峰,廖计华,张晓庆. 东海盆地西湖凹陷平湖斜坡K构造带油气沿断层走向运聚模式及控制因素. 天然气地球科学. 2024(03): 393-404 . 百度学术
3. 何贤科,李文俊,段冬平,荣乘锐,夏振通. 西湖凹陷平湖斜坡北段平湖组薄煤层与泥岩的微观岩石学特征及其沉积学意义. 海洋地质与第四纪地质. 2024(02): 210-222 . 本站查看
4. 李帅,俞伟哲,秦兰芝,张粲. 西湖凹陷平湖斜坡带物源-坡折耦合控砂模式. 海洋地质前沿. 2024(07): 36-44 . 百度学术
5. 张尚虎,黄建军,李昆,万丽芬,庄建建,王丹萍,王修平,蒋涔. 西湖凹陷孔雀亭地区复合圈闭发育模式与油气富集差异控制因素. 海洋地质与第四纪地质. 2023(01): 128-137 . 本站查看
6. 李盛谦,曾溅辉,刘亚洲,李淼,焦盼盼. 东海盆地西湖凹陷孔雀亭地区古近系平湖组储层成岩作用及孔隙演化. 岩性油气藏. 2023(05): 49-61 . 百度学术
7. 吕鹏,雷蕾,孙莉,王健伟,闫华,林立新,刘世鹏,葛家旺,吴文雯. 东海盆地西湖凹陷W气田平湖组煤系潮控砂体储层地震预测. 海洋地质前沿. 2023(10): 66-76 . 百度学术
8. 王安龙,翁冬子,任培罡,吕鹏,潘潞,吕艾新,魏锋. 西湖凹陷X气田致密气储层测井评价及产能预测. 海洋石油. 2023(04): 77-82 . 百度学术
9. 马耑月,宗廷博,巴魏魏,杨川. 环江地区长3段沉积微相特征研究. 石油化工应用. 2022(09): 94-100 . 百度学术
其他类型引用(3)