南黄海中部泥质区YE-2孔8.4 kaBP来的沉积特征

王利波, 杨作升, 赵晓辉, 邢磊, 赵美训, Yoshiki Saito, 范德江

王利波, 杨作升, 赵晓辉, 邢磊, 赵美训, Yoshiki Saito, 范德江. 南黄海中部泥质区YE-2孔8.4 kaBP来的沉积特征[J]. 海洋地质与第四纪地质, 2009, 29(5): 1-11. DOI: 10.3724/SP.J.1140.2009.05001
引用本文: 王利波, 杨作升, 赵晓辉, 邢磊, 赵美训, Yoshiki Saito, 范德江. 南黄海中部泥质区YE-2孔8.4 kaBP来的沉积特征[J]. 海洋地质与第四纪地质, 2009, 29(5): 1-11. DOI: 10.3724/SP.J.1140.2009.05001
WANG Libo, YANG Zuosheng, ZHAO Xiaohui, XING Lei, ZHAO Meixun, Yoshiki Saito, FAN Dejiang. SEDIMENTARY CHARACTERISTICS OF CORE YE-2 FROM THE CENTRAL MUD AREA IN THE SOUTH YELLOW SEA DURING LAST 8 400 YEARS AND ITS INTERSPACE COARSE LAYERS[J]. Marine Geology & Quaternary Geology, 2009, 29(5): 1-11. DOI: 10.3724/SP.J.1140.2009.05001
Citation: WANG Libo, YANG Zuosheng, ZHAO Xiaohui, XING Lei, ZHAO Meixun, Yoshiki Saito, FAN Dejiang. SEDIMENTARY CHARACTERISTICS OF CORE YE-2 FROM THE CENTRAL MUD AREA IN THE SOUTH YELLOW SEA DURING LAST 8 400 YEARS AND ITS INTERSPACE COARSE LAYERS[J]. Marine Geology & Quaternary Geology, 2009, 29(5): 1-11. DOI: 10.3724/SP.J.1140.2009.05001

南黄海中部泥质区YE-2孔8.4 kaBP来的沉积特征

基金项目: 

国家重点基础研究发展规划项目(2005CB422304)

国家自然科学基金项目(90211022,40676032)

详细信息
    作者简介:

    王利波(1983-),男,硕士,从事海洋沉积学研究,E-mail:wanglibowlb@gmail.com

  • 中图分类号: P736.21

SEDIMENTARY CHARACTERISTICS OF CORE YE-2 FROM THE CENTRAL MUD AREA IN THE SOUTH YELLOW SEA DURING LAST 8 400 YEARS AND ITS INTERSPACE COARSE LAYERS

  • 摘要: 对取自南黄海中部泥质区YE-2钻孔的6.1 m岩心进行了年代学、沉积学、微体古生物学和U37K'海表温度分析,结果显示,该孔沉积物基本类型为黏土质粉砂,8 400 aBP来的沉积环境可分为3段:8.4~6.9 kaBP为砂含量达13%的滨岸浅水环境;6.9~6.1 kaBP为砂含量8%的滨岸浅水向稳定浅海过渡的环境;6.1~0 kaBP为基本不含砂的黏土质粉砂稳定浅海环境。表层海水温度(SST)在滨岸浅水阶段较低,平均约14.2℃;在第2阶段从6.9 kaBP开始增加,最高达18.1℃,指示了黄海暖流的入侵;6.4 kaBP增温结束后,浅海阶段SST稳定在15~16℃之间,同时有3个SST平均为13.9℃的短暂低温期,大致为5.7~5.2、3.7~3.2和2.3~1.8 kaBP。在稳定浅海环境的黏土质粉砂地层中发现了罕见的3个砂质粉砂质粗粒级夹层,大致形成于6.0、3.9和2.5 kaBP,其粒度具有风暴沉积特征,重矿物表现出变质岩特点,U37K'-SST和黏土矿物指示其沉积环境的温度偏高。上述特征与岩心其他所有层位差异较大,显示了异地搬运的特征,可能为风暴沉积。根据AMS 14C测年,不考虑风暴沉积,YE-2岩心的沉积速率在50.7~75.6 cm/ka之间。
    Abstract: Based on chronological, sedimentological, micropaleontological and U37K'-sea surface temperature data of core YE-2 from the central mud area of the South Yellow Sea, evolution of the marine environment in the past 8 400 years and three interspace sandy layers in the core were studied. The core is mainly composed of clayey silt, and evolution of marine environment during the last 8 400 years experienced three phases:(1) Near-shore shallow water environment with sand content of 13% at about 8.4~6.9 kaBP; (2) Transitional marine environment from near-shore shallow water to stable shallow sea with sand content of 8% at about 6.9~6.1 kaBP; (3) Stable shallow sea environment with sediment without sand at about 6.1~0 kaBP. Sea surface temperature (SST) was low in the near-shore shallow water phase, about 14.2℃ on average. And the SST began to increase at the beginning of the second phase, about 6.9 kaBP, indicating intrusion of the Yellow Sea Warm Current with the highest SST up to 18.1℃. The SST stopped increasing at 6.4 kaBP. After that the SST was stable from 15℃ to 16℃ with some fluctuations in the shallow sea phase, excluding three low temperature periods at about 5.7~5.2, 3.7~3.2 and 2.3~1.8 kaBP respectively, whose temperature was about 13.9℃ on average. Three abnormal sandy silt coarse interlayers deposited in stable shallow sea environment,were speculated to have formed in 6.0, 3.9 and 2.5 kaBP respectively. Grain-sizes of the sandy layers show characteristics of storm deposition.Heavy mineral compositions of the sandy layers show characteristics of the origin of metamorphic rocks. U37K'-SST and clay mineral compositions of the sandy layers indicate a sedimentary environment with higher temperature than that in shallow sea. These characteristics of sandy layers are different from other sediment in the core, indicating great differentce between the sandy layer and other strata in provenance.Therefore, the three coarse layers are speculated to be the storm deposition. Without consideration of three storm depositions, the sedimentation rates in core YE-2 range from 50.7 to 75.6 cm/ka, based on AMS 14C data.
  • 天然气水合物主要赋存于大陆边缘和大陆永久冻土区[1-2],以层状、块状、分散状、胶结状、瘤状等形式赋存于地层中[3],可归结为两大类,即孔隙充填型和裂隙充填型。孔隙充填型水合物替代了孔隙中的流体,赋存于沉积物孔隙中[4];裂隙充填型水合物则以纯水合物形式与沉积物互相接触,一般存在于沉积物裂隙、地层裂隙、断层或泥火山等特殊地质体中。因此,由于赋存状态和岩石物理性质不同,针对孔隙充填型、裂隙充填型水合物往往采用不同的方法进行资源评估。然而,在野外勘探过程中,如何通过简便的方法识别水合物类型,尚无较为成熟的识别方法。

    已有钻探结果为认识孔隙充填型和裂隙充填型水合物打下了良好基础[5-6]。中国、美国、印度、韩国等钻探表明,中国南海神狐海域天然气水合物以孔隙充填为主[7],而美国墨西哥湾[8]、韩国Ulleung盆地[9]和印度K-G盆地[10-11]的细粒沉积物中均发现了裂隙充填型水合物。2006年,印度国家天然气水合物计划(NGHP01)在K-G盆地通过测井和钻探发现大量裂隙充填型水合物富集,测井和红外成像显示水合物多呈结核状、脉状、瘤状,充填于高角度(60°~120°)裂隙以及近水平层理中,其中充填水合物的裂隙层位厚度最大可达135 m,长度达数米以上[1, 5, 11]。孔隙充填型水合物和裂隙充填型水合物最主要的区别在于非微观尺度研究中,孔隙充填型水合物通常被看作是各向同性介质,而裂隙充填型水合物则是各向异性介质[12]。中国地质调查局广州海洋地质调查局2013年在中国南海北部海域开展的第二次水合物勘探GMGS2,在13个站位布设了21个钻井,并取得了保压岩心[13],其中5个站位发现了水合物,这一航次中同时发现了孔隙充填型和裂隙充填型水合物。

    由于孔隙充填型和裂隙充填型水合物的形成机理和地球物理响应都存在差异,因此,在无钻探井位的区域,准确识别两类水合物在沉积层中的分布形态,对于水合物资源评价、钻采安全以及环境等方面都至关重要[14]。目前为止,关于孔隙充填型和裂隙充填型水合物识别的文献相对较少。戴丹青等[15]利用中观衰减理论,用斑块饱和模型和双孔隙模型针对水合物在地层中的赋存状态,研究了地震频带内BSR界面反射系数的变化特征,并将研究结果应用到存在孔隙充填型水合物的南海神狐海域和存在裂隙充填型水合物的布莱克海台,模拟结果与实际情况吻合较好。王吉亮等[3]为识别印度东海岸NGHP01-10D井的水合物赋存状态,利用等效介质理论(EMT)和层状介质横向各向异性模型分别计算了水合物的饱和度,发现层状介质横向各向异性模型计算的结果与实际情况吻合较好,且存在的裂隙多为高角度裂隙。Lee等[16]利用各向同性模型和横向各向同性理论计算了韩国东海Ulleung 盆地UBGH2-3B井中水合物层的电阻率、纵波速度和水合物饱和度,研究结果表明井中局部范围水合物的赋存类型为裂隙充填型,气烟囱上半部分(5~14 mbsf)地层中主要赋存裂隙角度约7°的裂隙充填型水合物,下半部分(45~103 mbsf)中赋存高角度裂隙充填型水合物。Liu等[17]提出利用纵波速度和电阻率交汇图的斜率来识别两类水合物。Sriram等[18]对含孔隙充填型/裂隙充填型水合物地层的BSR进行了AVA(振幅随入射角的变化)模拟,通过模拟将两者的AVA类别进行归类,从而达到水合物在储层中赋存形态的识别。Liu等[19]通过岩石物理模拟分别计算了孔隙充填型和裂隙充填型水合物的密度和速度,提出用密度速度属性识别孔隙充填型及裂隙充填型水合物的方法。

    迄今为止,对水合物赋存形态识别的研究多以模拟手段为主来计算两种类型水合物储层的电阻率、纵波速度、地震衰减/反射系数及AVA/AVO属性。截至目前,未发现将模拟结果与实验实测数据进行对比分析的案例。本文针对孔隙充填型水合物选取了EMT(A、B模式)、BGTL和STPE模型,针对裂隙充填型水合物(本文仅以裂隙空间发育的纯水合物为研究对象)(图1a)选取了横向各向同性理论(TIT),模拟了两种类型水合物的密度和声波速度与水合物体积分数间的变化特征,并通过声学实验获取了孔隙充填型和裂隙充填型水合物的纵波速度特征。

    图  1  裂隙充填型水合物的类型[20]
    a. 裂隙空间发育的纯水合物,b. 黏土质沉积物夹杂含水合物的砂质沉积介质, c. 黏土沉积物中发育连续微裂隙中的水合物,d. 黏土沉积物中不连续微裂隙中的水合物。
    Figure  1.  Types of fracture-filling hydrates[20]
    a. Pure hydrate in fracture space, b. sandy sedimentary medium with hydrate in clay sediment,c. hydrate in continuous microfractures developed in clay sediments,d. hydrate in discontinuous microfractures in clay sediments.

    Liu等[18]、胡高伟[21]、王吉亮等[3]等关于孔隙充填型水合物岩石物理模拟的研究结果表明,等效介质理论(EMT)、BGTL理论和简化三相方程(STPE)的计算结果与实验测试波速的吻合程度相比其他岩石物理模型更具优势。因此,本文针对孔隙充填型水合物选取了EMT-A/B、BGTL和STPE岩石物理模型,计算了含孔隙充填型水合物的南海沉积物(孔隙度为0.4)的密度和纵波速度。表1为模拟计算所用的南海沉积物的弹性参数及含量表,模型计算公式见附录。

    表  1  南海泥质黏土的体积模量(K)、剪切模量(G)、密度(ρ)和含量(C[22]
    Table  1.  The bulk modulus (K), shear modulus (G), density (ρ) and content (C) of the silty clay in the South China Sea [22]
    沉积物组分K/GPaG/GPaρ/(g/cm3C/%
    方解石76.8322.7114
    石英36.6452.6528
    长石76262.7112
    云母62412.6826
    黏土20.96.852.5820
    水合物5.62.40.9
    2.501.03
    下载: 导出CSV 
    | 显示表格

    由于EMT-A和BGTL的密度计算公式相同,因此,两者模拟的密度曲线重合(图2),而EMT-B模式中,水合物形成后充当骨架的一部分,同时使得沉积物孔隙度减小,该模式计算密度时对沉积物孔隙度进行了修正,但最终计算结果与STPE相同。各模式计算的孔隙充填型水合物介质的密度尽管有差异,但最大差值仅为1.25%。

    图  2  EMT、BGTL、STPE模型计算的含孔隙充填型水合物的密度
    Figure  2.  The density of pore-filling hydrate calculated by EMT、BGTL and STPE models

    3类岩石物理模型计算的纵波速度均随水合物体积分数的增加而增大(图3),但不同模型的波速计算方式不同,波速初始起跳值也不同。当沉积物孔隙中未生成水合物(体积分数为0)时,EMT-A/B模式计算的沉积物纵波速度大于BGTL和STPE模型计算的波速,3类模型计算的纵波速度大小分两种情况:①当水合物体积分数小于17%时,波速大小为:EMT-B>EMT-A>STPE>BGTL;②水合物体积分数大于17%时,EMT-B>STPE>EMT-A>BGTL。EMT-B模式计算的波速值最大,BGTL的计算值最小,图3中横坐标体积分数40%对应水合物饱和度为100%。

    图  3  EMT、BGTL、STPE模型计算的含孔隙充填型水合物的纵波速度
    Figure  3.  The P-wave velocity of pore-filling hydrate calculated by EMT、BGTL and STPE models

    裂隙充填型水合物模型用横向各向同性理论(TIT)进行计算。模型由两个端元组成(图4),端元Ⅰ为裂隙,100%由纯水合物充填,端元Ⅱ为饱和水的南海沉积物(泥质黏土,孔隙度为0.4),其弹性参数及含量见表1。端元Ⅱ的体积模量和剪切模量基于Hill平均方程(B-5)和STPE模型计算。

    图  4  裂隙充填型水合物端元模型[23]
    Ⅰ为100%纯水合物充填的裂隙,Ⅱ为饱和水泥质黏土。
    Figure  4.  Fracture-filling hydrate end member model[23]
    Ⅰ: Fracture filled by pure hydrate; Ⅱ: Silty clay with saturated water.

    由横向各向同性理论计算的含裂隙充填型水合物的南海沉积介质的密度和纵波速度结果(图5)发现,随着水合物体积分数增大,沉积介质的总密度减小,且当模型中端元Ⅰ占比为100%时,裂隙中为纯水合物,此时介质密度最小,水合物的密度为0.93 g/cm3;纵波速度曲线最小值代表南海沉积物的模型理论速度(1.84 km/s),最大值代表纯水合物的速度(3.08 km/s),纵波速度随着水合物体积分数增大而增大。

    图  5  横向各向同性理论计算的含裂隙充填型水合物的密度和纵波速度
    Figure  5.  The density and P-wave velocity of fracture-filling hydrate calculated by layered media model

    孔隙充填型水合物声学特性模拟实验装置如图6。实验之前需要先对弯曲元超声测速探头的固有声波走时和时域反射技术(TDR)测定沉积物中含水量的误差进行标定。速度测试时需要减去超声探头的固有声波走时,标定结果为$y = 2.242x \times $$ {10^{ - 3}} - 0.01246$$y = 9.291x \times {10^{ - 4}} - 9.812 \times {10^{ - 3}}$,弯曲元换能器中纵波的固有传播时间为5.56 μs,横波固有传播时间为10.56 μs;TDR实测值与理论值间的误差范围为0.17%~6.75%,由于TDR测量的结果为一个场平均值,因此,TDR系统测试精度满足实验要求。

    图  6  天然气水合物声学模拟实验装置
    Figure  6.  Acoustic experimental apparatus for natural gas hydrate

    实验模拟过程为:将饱和度水南海沉积物装入高压反应釜,并加入SDS(十二烷基硫酸钠)溶液(加快水合物的生成),通入纯甲烷气体(纯度99.99%)使釜内压力达到目标压力(约7 MPa),并放置24 h使甲烷气溶解到SDS溶液中,随后开启制冷系统使反应釜内部温度到约1 ℃。随着水合物生成,由TDR系统测试水合物的饱和度,弯曲元声波系统测试波速。

    实验中,孔隙充填型水合物的密度由沉积物中矿物及生成的水合物的密度及其含量计算得来。纵波速度为实验测试结果(图7)。由于南海沉积物中短期很难生成高饱和度水合物,因此,反应釜中生成水合物的最大体积分数为9.6%(对应饱和度为24%)。基于实验波速数据发现,纵波速度随水合物体积分数呈正比趋势,当水合物体积分数为9.6%时,沉积介质纵波速度达到最大值2.06 km/s。

    图  7  南海沉积物中孔隙充填型水合物的纵波速度和密度
    Figure  7.  P-wave velocity and density of pore-filling hydrate in deposits of the South China Sea

    先利用实验装置(图8)进行裂隙充填型水合物实验合成,其过程为:利用雾化喷头将纯水雾化后喷洒到液氮中,合成所需量的冰粉,将冰粉与沉积物按照比例装入多孔塑胶管中,反应釜密封后通入甲烷气使釜内压力达到水合物合成所需目标压力,再将反应釜置于循环水域控温装置中进行降温(约1 ℃),直至水合物合成完毕。最后开釜取出样品,连同塑胶管一同置于岩心夹持器(图9)上,利用超声探头发射/接收声波,进行波速测试。

    图  8  裂隙充填型水合物合成实验装置
    Figure  8.  Experimental apparatus for fracture-filling hydrate synthesis
    图  9  岩心夹持测试波速装置
    Figure  9.  Rock holding device for testing wave velocity

    不加入水合物催化溶液时,南海沉积物中短期很难生成孔隙充填型水合物,且长时间内合成的水合物饱和度也很低[21]。该次裂隙充填型水合物合成实验中未向沉积物中加入SDS溶液,一个实验周期(约7天)内沉积物孔隙中几乎不会有水合物生成,实验介质可视为纯裂隙充填型水合物。

    含裂隙充填型水合物(裂隙为水平裂隙)的南海沉积介质的密度与纵波速度结果中(图10),密度计算公式为$\rho = {\eta _1}{\rho _{\rm{h}}} + {\eta _2}\left[ {{\rho _{\rm{s}}}(1 - \phi ) + {\rho _{\rm{w}}}\phi } \right]$,其中,${{\rm{\eta }}_{\rm{1}}}$${{\rm{\eta }}_{\rm{2}}}$分别为裂隙和沉积物的体积分数,${\rho _{\rm{h}}}$为纯水合物的密度,计算取值为0.93 g/cm3${\rho _{\rm{s}}}$为南海泥质黏土沉积物的骨架密度(2.66 g/cm3),${\rm{\phi }}$为沉积物孔隙度,${\rho _{\rm{w}}}$为水的密度(1 g/cm3);纵波速度为实验超声测试结果。当水合物体积分数为0时,整个沉积介质的密度即为沉积物的骨架密度(2.66 g/cm3),随着水合物体积分数增大,整个沉积介质(水合物+沉积物)的密度呈减小趋势,而当水合物体积分数为100%时,介质的密度即为纯水合物的密度值(0.93 g/cm3)。水合物体积分数为0和100%分别对应100%南海沉积物和100%纯水合物的波速,当水合物体积分数增大时,纵波速度也随之增大。

    图  10  含裂隙充填型水合物的南海沉积介质的密度与纵波速度
    Figure  10.  Density and P-wave velocity of sediment medium of South China Sea with fracture-filling hydrate

    两种类型水合物的岩石物理模拟和实验纵波速度结果中(图11),波速均随水合物体积分数增大而增大。实验中生成的孔隙充填型水合物的最大体积分数为9.6%(对应饱和度为24%)。孔隙充填型水合物的纵波速度在体积分数小于4%(对应饱和度为10%)时低于岩石物理模拟的波速,推断造成该现象的主要原因是当南海沉积物孔隙空间生成的水合物饱和度较低时,水合物呈悬浮状分散于沉积物孔隙内,导致声波信号发生衰减,且这种悬浮状态对沉积物骨架支撑几乎不起作用,故当水合物体积分数(饱和度)较低时,实验模拟的波速低于岩石物理模拟值。而当体积分数大于4%时,实验波速趋向于BGTL模型计算的理论波速。横向各向同性理论(TIT)计算的裂隙充填型(水平裂隙)水合物理论波速在水合物体积分数较大时(>40%)与实验测试的纵波速度吻合较好,但当水合物体积分数较小时,两者间的最大差为0.27 km/s,整体吻合较好。

    图  11  岩石物理模拟和实验测试的纵波速度
    Figure  11.  The P- wave velocity of petrophysical simulation and experimental test

    通过将密度与速度参数结合,计算了两种类型水合物的阻抗(ρVp)(图12)和$\rho {\sqrt V _{\rm{p}}}$属性(图13)。由于裂隙填充型水合物的实验阻抗斜率非严格负值,因此,基于阻抗无法进行两种类型水合物的判别。岩石物理模拟和实验测试的$\rho {\sqrt V _{\rm{p}}}$属性中,孔隙充填型水合物和裂隙充填型水合物曲线的斜率符号相反:孔隙充填型水合物的$\rho {\sqrt V _{\rm{p}}}$属性曲线斜率为正,而裂隙充填型水合物为负。孔隙充填型和裂隙充填型水合物的这一特征为两种类型水合物的判别提供了有效手段。

    图  12  岩石物理模拟和实验测试的孔隙充填型和裂隙充填型水合物的阻抗
    Figure  12.  The impedance of pore-filling and fracture-filling hydrates obtained by petrophysical simulation and experimental test
    图  13  岩石物理模拟和实验测试的孔隙充填型和裂隙充填型水合物的$\rho {\sqrt V _{\rm{p}}}$属性
    Figure  13.  The $\rho {\sqrt V _{\rm{p}}}$properties of pore-filling and fracture-filling hydrates obtained by petrophysical simulation and experimental test

    中国第二次天然气水合物钻探航次(GMGS2)在中国南海珠江口盆地东部第16井中同时发现了孔隙充填和裂隙充填型水合物,利用该井的密度和测井资料对$\rho {\sqrt V _{\rm{p}}}$属性识别水合物类型的方法进行验证。该16井中孔隙充填和裂隙充填型水合物平均饱和度分别为42%和50%,其中孔隙充填和裂隙充填型水合物分别位于井深(泥线以下:BML)189~198和10~20 m[18, 24]

    该井上部10~20 m验证结果(图14)表明该范围水合物类型以裂隙充填型水合物为主,夹杂部分孔隙充填型水合物。井深10~11.6、14~18.5 m深度含水合物层与背景层的$\rho {\sqrt V _{\rm{p}}}$属性差值小于零,指示该范围内存在裂隙充填型水合物,12~13.6、18.5~20 m深度$\rho {\sqrt V _{\rm{p}}}$属性差值大于零,代表赋存的水合物类型为孔隙充填型。井底部189~198 m赋存的水合物类型以孔隙充填型为主(图15),在194~196 m深度夹杂少量裂隙充填型水合物。从验证结果来看,基本与野外观测结果一致。

    图  14  GMGS2-16井上部含裂隙充填型水合物储层$\rho {\sqrt V _{\rm{p}}}$属性差值结果
    Figure  14.  Result of $\rho {\sqrt V _{\rm{p}}}$ attribute difference of fracture-filling hydrate reservoir in the upper part range of MGMS2-16 Well
    图  15  GMGS2-16井下部含孔隙充填型水合物储层$\rho {\sqrt V _{\rm{p}}}$属性差值结果
    Figure  15.  Result of $\rho {\sqrt V _{\rm{p}}}$ attribute difference of pore-filling hydrate reservoir in the lower part of GMGS2-16 Well

    密度和纵波速度间的差异提供了识别孔隙充填型水合物和裂隙充填型水合物在沉积物中赋存形态的基础。通过结合岩石物理模型和实验模拟,获得了南海沉积物中孔隙充填型和裂隙充填型水合物的密度和纵波速度参数,测试了两种类型水合物的阻抗和$\rho {\sqrt V _{\rm{p}}}$属性。研究结果表明基于$\rho {\sqrt V _{\rm{p}}}$属性可以有效识别孔隙型充填和裂隙充填型水合物:对于含孔隙充填型水合物的南海沉积物,岩石物理模型的$\rho {\sqrt V _{\rm{p}}}$属性与实验测试的$\rho {\sqrt V _{\rm{p}}}$属性曲线都是正斜率;对于含裂隙充填型水合物的南海沉积物,岩石物理模型和实验测试的$\rho {\sqrt V _{\rm{p}}}$属性都为负斜率。对于含孔隙充填型水合物的南海沉积介质,当水合物体积分数大于4%时,岩石物理模型计算的$\rho {\sqrt V _{\rm{p}}}$属性与实验数据吻合较好,尤其是STPE模型与实验结果吻合最佳,因此,STPE模型在计算南海沉积物中水合物的饱和度时相比其他模型更具优势。对于裂隙充填型水合物,当水合物体积分数较低(<40%)时,横向各向同性理论计算的$\rho {\sqrt V _{\rm{p}}}$属性和实验测试结果存在明显差异。因此,当水合物体积分数低于40%时,计算水合物的$\rho {\sqrt V _{\rm{p}}}$属性过程中需要对模型进行适当修正。

    最后利用$\rho {\sqrt V _{\rm{p}}}$属性对GMGS2航次16井赋存的孔隙充填和裂隙充填型水合物进行了验证表明结果与实际钻探结果一致,因此,利用$\rho {\sqrt V _{\rm{p}}}$属性识别孔隙充填和裂隙充填型水合物具有理论和实际可行性。

    Helgerud等[25-26]对水合物微观分布形态进行了研究,认为水合物与沉积物颗粒间的共生关系主要是充填于沉积物孔隙中,或胶结了沉积物骨架,与骨架呈支撑形态(图A-1)。对于模式A,水合物的生成替代了沉积物孔隙中的流体并占据了流体空间,因而该模式中,水合物也被认为是孔隙流体的一部分;模式B中,水合物被当作沉积物骨架的一部分,孔隙空间中只有流体,该模式对沉积物骨架的密度、体积模量和剪切模量产生了影响[21];模式C中,孔隙度的降低等同于模式B,另一方面,岩石骨架的体积模量和剪切模量发生了改变。南海沉积物中孔隙充填型水合物模拟选用A(水合物充当流体部分)和B(水合物充当骨架部分)模式进行计算。

      A-1  水合物在沉积物中的3种赋存形态
      A-1.  Three types of hydrate occurrence in sediment

    (1)模式A

    沉积介质的纵波速度及密度计算公式如下:

    $$\tag{A-1}\begin{array}{l} {V_{\rm{p}}} = \sqrt {\frac{{{K_{{\rm{sat}}}} + \displaystyle\frac{4}{3}{G_{{\rm{sat}}}}}}{{{\rho _{\rm{b}}}}}} \\ {\rho _{\rm{b}}} = (1 - \phi ){\rho _{\rm{s}}} + \phi {\rho _{\rm{f}}} \end{array}$$

    式中,KsatGsat分别是等效介质的体积模量和剪切模量;${\rho _{\rm{s}}}$${\rho _{\rm{f}}}$分别为岩石固相和流体相的体积密度。ϕc为临界孔隙度,一般取0.36~0.40 [27]。当沉积物中充填体积模量为Kf的流体时,可以根据Gassmann方程,通过如下公式得到沉积物的体积模量Ksat和剪切模量Gsat

    $$\tag{A-2}{K_{\rm{sat}}} = {K_{\rm{ma}}}\frac{{\phi {K_{\rm{dry}}} - (1 + \phi ){K_{\rm f}}{K_{\rm{dry}}}/{K_{\rm{ma}}} + {K_{\rm f}}}}{{(1 - \phi ){K_{\rm f}} + \phi {K_{\rm{ma}}} - {K_{\rm f}}{K_{\rm{dry}}}/{K_{\rm{ma}}}}}$$
    $$\tag{A-3}{G_{{\rm{sat}}}} = {G_{{\rm{dry}}}}$$

    式中,Kma为岩石固相的体积模量,KdryGdry分别为干岩石的体积模量和剪切模量,Kf为流体的体积模量。模式A中,孔隙中有水合物生成,因此,Kf计算公式为:

    $$\tag{A-4}{K_{\rm f}} = {\left[\frac{{1 - {S_{\rm h}}}}{{{K_{\rm w}}}} + \frac{{{S_{\rm h}}}}{{{K_{\rm h}}}}\right]^{ - 1}}$$

    式中,ShKh分别代表水合物占孔隙的体积分数和水合物的体积模量,Kw为水的体积模量。

    式(A-2)中,KdryGdry的计算公式为:

    $$\tag{A-5}{K_{\rm{dry}}} \!= \left\{ {\begin{array}{*{20}{c}} {{{\left[\frac{{\phi /{\phi _{\rm c}}}}{{{K_{\rm{HM}}} + \displaystyle\frac{4}{3}{G_{\rm{HM}}}}} + \frac{{1 - \phi /{\phi _{\rm c}}}}{{{K_{\rm{ma}}} + \displaystyle\frac{4}{3}{G_{\rm{HM}}}}}\right]}^{ - 1}} \!-\! \displaystyle\frac{4}{3}{G_{\rm{HM}}},\phi < {\phi _{\rm c}}} \\ {{{\left[\frac{{(1 - \phi )/(1 - {\phi _{\rm c}})}}{{{K_{\rm{HM}}} + \displaystyle\frac{4}{3}{G_{\rm{HM}}}}} + \frac{{(\phi - {\phi _{\rm c}})/(1 - {\phi _{\rm c}})}}{{\displaystyle\frac{4}{3}{G_{\rm{HM}}}}}\right]}^{ - 1}} \!-\! \displaystyle\frac{4}{3}{G_{\rm{HM}}},\phi \geqslant {\phi _{\rm c}}} \end{array}} \right.$$
    $$\tag{A-6}{G_{\rm{dry}}} \!= \left\{ {\begin{array}{*{20}{l}} {{{\left[\displaystyle\frac{{\phi /{\phi _{\rm c}}}}{{{G_{\rm{HM}}} + Z}} \!+\! \frac{{1 - \phi /{\phi _{\rm c}}}}{{G + Z}}\right]}^{ - 1}} - Z,\phi < {\phi _{\rm c}}} \\ {{{\left[\displaystyle\frac{{(1 \!-\! \phi )/(1 - {\phi _{\rm c}})}}{{{G_{\rm{HM}}} + Z}} \!+\! \frac{{(\phi \!-\! {\phi _{\rm c}})/(1 \!-\! {\phi _{\rm c}})}}{Z}\right]}^{ - 1}} \!- \!Z,\phi \geqslant {\phi _{\rm c}}} \end{array}} \right.$$
    $$\tag{A-7}Z = \frac{{{G_{\rm{HM}}}}}{6}\left(\frac{{9{K_{\rm{HM}}} + 8{G_{\rm{HM}}}}}{{{K_{\rm{HM}}} + 2{G_{\rm{HM}}}}}\right).$$

    其中,

    $$\tag{A-8}\begin{split} {K_{\rm{HM}}} =& {\left[\displaystyle\frac{{{n^2}{{(1 - {\phi _{\rm c}})}^2}{G_{\rm{ma}}}^2}}{{18{\pi ^2}{{(1 - \nu )}^2}}}P\right]^{\frac{1}{3}}}\\ {G_{\rm{HM}}} =&\displaystyle\frac{{5 - 4\nu }}{{5(2 - \nu )}}{\left[\frac{{3{n^2}{{(1 - {\phi _{\rm c}})}^2}{G_{\rm{ma}}}^2}}{{2{\pi ^2}{{(1 - \nu )}^2}}}P\right]^{\frac{1}{3}}} \end{split} $$

    式(A-8)中,P为有效压力;KmaGma分别是岩石骨架的体积模量、剪切模量,由公式(A-13)计算;ν是岩石骨架的泊松比,且$\nu = 0.5({K_{\rm{ma}}} - \displaystyle\frac{2}{3}{G_{\rm{ma}}})/({K_{\rm{ma}}} +$$ {G_{\rm{ma}}}/3) $n为临界孔隙度时单位体积内颗粒平均接触的数目,一般取8~9.5[28-29]

    (2)模式B

    模式B中水合物被认为是岩石骨架的一部分,产生了两个效应:一个是使孔隙度减小,另一个是改变了骨架的体积模量和剪切模量。因此,在模式A的基础上,需对沉积物孔隙度进行修正,即ϕr=ϕ(1-Sh)。同时,应将水合物作为矿物组分代入公式(A-9)来计算岩石的KmaGma。此外,沉积物孔隙中只有水,孔隙流体密度和体积模量等直接用水的替代。

    BGTL理论(Biot-Gassmann Theory modified by Lee的简称)建立在经典的BGT理论(Biot-Gassmann Theory)上,在预测速度时不仅考虑了分压的影响,而且还考虑了岩石的孔隙度、固结度等因素的影响,其公式为:

    $$\tag{A-9}{V_{\rm{p}}} = \sqrt {\frac{{k + 4\mu /3}}{\rho }} ,{V_{\rm s}} = \sqrt {\frac{\mu }{\rho }} $$

    式中,$k$$\mu $分别为沉积介质的体积模量和剪切模量。

    $$\tag{A-10}k = {k_{\rm{ma}}}(1 - \beta ) + {\beta ^2}M,\frac{1}{M} = \frac{{\beta - \phi }}{{{k_{\rm{ma}}}}} + \frac{\phi }{{{k_{\rm{fl}}}}}$$

    式中,${k_{\rm{ma}}}$为岩石骨架的体积模量,${k_{\rm{fl}}}$为孔隙中流体的体积模量;$\beta $是Biot系数,表征了流体体积变化与岩石体积变化的比值,与沉积介质的孔隙度有关;M为一模量,表征了沉积介质体积不变的情况下,将一定量的水压入沉积介质所需要的静水压力增量[30]。BGTL理论假设沉积介质速度比率与沉积物基质速度比率之间存在如下关系:

    $$\tag{A-11}{V_s} = {V_{\rm{p}}}G\alpha {(1 - \phi )^n}$$

    式中,$\alpha $为基质部分的Vs/VpG是与沉积物中黏土含量有关的常数,n取决于分压大小及岩石的固结程度。

    由公式(B-2)和(B-3)得出沉积介质的剪切模量为:

    $$\tag{A-12}\mu = \frac{{{\mu _{\rm{ma}}}{G^2}{{(1 - \phi )}^{2n}}k}}{{{k_{\rm{ma}}} + 4{\mu _{\rm{ma}}}[1 - {G^2}{{(1 - \phi )}^{2n}}]/3}}$$

    式中,${\mu _{\rm{ma}}}$为岩石骨架的剪切模量。${k_{\rm{ma}}}$${\mu _{\rm{ma}}}$由Hill平均方程计算[31]:

    $$\tag{A-13}\begin{array}{l}{k_{\rm{ma}}} = \displaystyle\frac{1}{2}\left[\sum\limits_{i = 1}^m {{f_i}} {k_i} + {\left(\sum\limits_{i = 1}^m {{f_i}/{k_i}} \right)^{ - 1}}\right],\\ {\mu _{\rm{ma}}} = \displaystyle\frac{1}{2}\left[\sum\limits_{i = 1}^m {{f_i}{\mu _i} + {{\left(\sum\limits_{i = 1}^m {{f_i}/{\mu _i}} \right)}^{ - 1}}} \right] \end{array}$$

    式中,m为岩石固相部分中矿物的种数,${f_i}$为第i种矿物占固相部分的体积分数;${k_i}$${\mu _i}$分别为第i种矿物的体积模量和剪切模量。

    利用Lee[32]的简化三相方程(STPE)可以计算各向同性的气体水合物填充沉积物孔隙空间的速度。Lee推导出低频下气体水合物填充沉积物孔隙空间的体积模量和剪切模量,用于测井和地震数据:

    $$\tag{A-14}k = {K_{\rm{ma}}}\left( {1 - {{\rm{\beta }}_{\rm{p}}}} \right){\rm{ + }}{{\rm{\beta }}_{\rm{p}}}^{\rm{2}}{K_{{\rm{av}}}},{\rm{\mu }} = {{\rm{\mu }}_{\rm{ma}}}\left( {{\rm{1 - }}{{\rm{\beta }}_{\rm{s}}}} \right)$$
    $$\tag{A-15}\frac{1}{{{K_{\rm{av}}}}}{\rm{ = }}\frac{{{{\rm{\beta }}_{\rm{p}}}{\rm{ - }}{\rm{\phi }}}}{{{K_{\rm{ma}}}}} + \frac{{{{\rm{\phi }}_{\rm w}}}}{{{k_{\rm w}}}} + \frac{{{{\rm{\phi }}_{\rm h}}}}{{{K_{\rm h}}}},{{\rm{\beta }}_{\rm p}} = \frac{{{{\rm{\phi }}_{\rm{as}}}\left( {1 + {\rm{\alpha }}} \right)}}{{1 + {\rm{\alpha }}{{\rm{\phi }}_{\rm{as}}}}},{{\rm{\beta }}_{\rm s}} = \frac{{{{\rm{\phi }}_{\rm {as}}}\left( {1 + {\rm{\gamma \alpha }}} \right)}}{{1 + {\rm{\gamma \alpha }}{{\rm{\phi }}_{\rm {as}}}}}$$
    $$\tag{A-16}{{\rm{\phi }}_{{\rm{as}}}} = {{\rm{\phi }}_{\rm w}} + {\rm{\varepsilon }}{{\rm{\phi }}_{\rm h}},{{\rm{\phi }}_{\rm w}} = (1 - {S_{\rm h}}){\rm{\phi }},{\rm{\phi }}h = {S_{\rm h}}{\rm{\phi }}$$

    式中,α为固结参数,γ是与剪切模量相关参数,可由公式γ=(1+2α)/(1+α)得到,KmaKwKh是骨架、水和水合物的体积模量,μma是骨架的剪切模量,ϕ为孔隙度。参数ε是水合物形成使沉积物骨架发生硬化的降低量,Lee推荐使用ε=0.12为建模数值。气体水合物填充沉积物孔隙空间的纵横波速度可由下式得:

    $$\tag{A-17}{V_{\rm{p}}} = \sqrt {\frac{{k + 4\mu /3}}{\rho }} ,{V_{\rm{s}}} = \sqrt {\frac{\mu }{\rho }} $$
    $$\tag{A-18}{\rho _{\rm{b}}} = {\rho _{\rm{s}}}\left( {1 - {\rm{\phi }}} \right) + {\rho _{\rm{w}}}{\rm{\phi }}\left( {1 - {S_{\rm{h}}}} \right){\rm{ + }}{\rho _{\rm{h}}}{\rm{\phi }}{S_{\rm{h}}}$$

    式中,$\rho $是水合物填充沉积物孔隙空间的体积密度,固结参数α取决于沉积物的有效压力和固结程度。Mindlin[33]表明体积模量和剪切模量取决于有效压力的1/3次方。基于Mindlin的这一理论,深度或有效压力相关的α可由下式得:

    $$\tag{A-19}{{\rm{\alpha }}_{\rm{i}}} = {{\rm{\alpha }}_0}{({P_0}/{P_i})^n} \approx {{\rm{\alpha }}_0}{({d_0}/{d_i})^n}$$

    式中,α0是有效压力P0或深度d0处的固结参数,αi是有效压力Pi或深度di处的固结参数。

    层状介质的两端元模型(图7)表示为:

    $$\tag{B-1}\left\langle G \right\rangle = {\eta _1}{G_1} + {\eta _2}{G_2},{\left\langle {\frac{1}{G}} \right\rangle ^{ - 1}} = {\left( {\frac{{{\eta _1}}}{{{G_1}}} + \frac{{{\eta _2}}}{{{G_2}}}} \right)^{ - 1}}$$

    式中,${G_1}$${G_{\rm{2}}}$分别代表模型中组分1和组分2的任意弹性参数或参数组合,${\eta _1}$${\eta _{\rm{2}}}$分别代表组分1和组分2的体积分数。模型相速度可以根据White由拉梅常数[34]求得:

    $$\tag{B-2}{V_p} = {\left( {\frac{{A{{\sin }^2}\varphi + C{{\cos }^2}\varphi + L + Q}}{{2\rho }}} \right)^{1/2}}$$
    $$\tag{B-3}{V_s}^H = {\left( {\frac{{N{{\sin }^2}\varphi + L{{\cos }^2}\varphi }}{\rho }} \right)^{1/2}}$$

    其中:

    $$\tag{B-4}A = \left\langle {\frac{{4\mu (\lambda + \mu )}}{{(\lambda + 2\mu )}}} \right\rangle + {\left\langle {\frac{1}{{(\lambda + 2\mu )}}} \right\rangle ^{ - 1}}{\left\langle {\frac{\lambda }{{(\lambda + 2\mu )}}} \right\rangle ^2}$$
    $$\tag{B-5}C = {\left\langle {\frac{1}{{(\lambda + 2\mu )}}} \right\rangle ^{ - 1}}$$
    $$\tag{B-6}F = {\left\langle {\frac{1}{{(\lambda + 2\mu )}}} \right\rangle ^{ - 1}}\left\langle {\frac{\lambda }{{(\lambda + 2\mu )}}} \right\rangle $$
    $$\tag{B-7}L = {\left\langle {\frac{1}{\mu }} \right\rangle ^{ - 1}}$$
    $$\tag{B-8}N = \left\langle \mu \right\rangle $$
    $$\tag{B-9}\rho = \left\langle \rho \right\rangle $$
    $$\tag{B-10}Q \!=\! \sqrt {{{[(A \!-\! L){{\sin }^2}\varphi \!-\! (C \!-\! L){{\cos }^2}\varphi ]}^2} \!+\! 4{{(F \!+\! L)}^2}{{\sin }^2}\varphi {{\cos }^2}\varphi } $$

    其中,VpVsH分别代表纵波速度和水平极化横波速度,参数$ \varphi $代表入射波传播方向相对于裂隙对称轴的夹角。

    根据Thomson[35],群速度与相速度间的关系可表述为:

    $$\tag{B-11}{V_{\rm{p}}}\left( {\rm{\varphi }} \right){\rm{ = }}G{V_{\rm{p}}}\left( {{{\rm{\varphi }}_{\rm{g}}}} \right),{V_{\rm{s}}}^H\left( {\rm{\varphi }} \right) = G{V_{\rm{s}}}^H\left( {{{\rm{\varphi }}_{\rm{g}}}} \right)$$

    式中:GVp和GVsH分别代表VpVsH的群速度,$ {\varphi }_{g} $是射线方向和裂隙同相轴间的夹角,可由下式求得:

    纵波:

    $$\tag{B-12}\tan {{\rm{\varphi }}_{\rm{g}}} = {\rm{tan}}{\rm{\varphi }}[1 + 2\delta + {\rm{4}}(\varepsilon - \delta ){\rm{si}}{{\rm{n}}^2}{\rm{\varphi }}]$$

    SH波:

    $$\tag{B-13}\tan {{\rm{\varphi }}_{\rm g}} = \tan {\rm{\varphi }}(1 + 2\gamma )$$

    其中:${\rm{\gamma }} = \displaystyle\frac{{N - L}}{{2L}},{\rm{\delta }} = \frac{{{{(F + L)}^2} - {{(C - L)}^2}}}{{2C(C - L)}},{\rm{\varepsilon }} = \frac{{A - C}}{{2C}}$,γ,$ \mathrm{\delta } $$ \mathrm{\varepsilon } $是表征横向各向同性介质属性的参数。

  • [1]

    Kim J M, Kennett J P. Paleoenvironmental changes associated with the Holocene marine transgression, Yellow Sea (Hwanghae)[J]. Marine Micropaleontology, 1998, 34(1-2):71-89.

    [2]

    Kim J M, Kucera Michal. Benthic foraminifer record of environmental changes in the Yellow Sea (Hwanghae) during the last 15000 years[J]. Quaternary Science Reviews, 2000, 19(11):1067-1085.

    [3] 刘健, 李绍全, 王圣洁, 等. 末次冰消期以来黄海海平面变化与黄海暖流的形成[J]. 海洋地质与第四纪地质, 1999, 19(1):13-24.

    [LIU Jian, LI Shaoquan, WANG Shengjie, et al. Sea level changes of the Yellow Sea and formation of the Yellow Sea warm current since the last deglaciation[J]. Marine Geology and Quaternary Geology, 1999, 19(1):13-24.]

    [4]

    Liu Jian, Saito Y, Wang Hong, et al. Sedimentary evolution of the Holocene subaqueous clinoform off the Shandong Peninsula in the Yellow Sea[J]. Marine Geology, 2007, 236(3-4):165-187.

    [5] 李铁刚, 李绍全, 苍树溪, 等. YSDP102钻孔有孔虫动物群与南黄海东南部古水文重建[J]. 海洋与湖沼, 2000, 31(6):588-595.

    [LI Tiegang, LI Shaoquan, CANG Shuxi, et al. Paleo-hydrological reconstruction of the southern Yellow Sea inferred from foraminiferal fauna in core YSDP102[J]. Oceanologia et Limnologia Sinica, 2000, 31(6):588-595.]

    [6] 李铁刚, 江波, 孙荣涛, 等. 末次冰消期以来东黄海暖流系统的演化[J]. 第四纪研究, 2007, 27(6):945-954.

    [LI Tiegang, JIANG Bo, SUN Rongtao, et al. Evolution pattern of warm current system of the East China Sea and the Yellow Sea since the last deglaciaion[J]. Quaternary Sciences, 2007, 27(6):945-954.]

    [7]

    Xiang Rong, Yang Zuosheng, Saito Y, et al. Paleoenvironmental changes during the last 8400 years in the southern Yellow Sea:Benthic foraminiferal and stable isotopic evidence[J]. Marine Micropaleontology, 2008, 67(1-2):104-119.

    [8] 申顺喜, 李安春, 袁巍. 南黄海中部的低能沉积环境[J]. 海洋与湖沼, 1996, 27(5):518-523.

    [SHEN Shunxi, LI Anchun, YUAN Wei. Low energy environment of the central south Yellow Sea[J]. Oceanologia et Limnologia Sinica, 1996, 27(5):518-523.]

    [9] 秦蕴珊, 赵一阳, 陈丽蓉, 等. 黄海地质[M]. 北京:海洋出版社, 1989.[QIN Yunshan, ZHAO Yiyang, CHEN Lirong, et al. Geology of the Yellow Sea[M]. Beijing:China Ocean Press, 1989.]
    [10] 赵一阳, 朴龙安, 秦蕴珊, 等. 南黄海沉积学研究新进展——中韩联合调查[J]. 海洋科学, 1998(1):34-37.[ZHAO Yiyang, Park Y, QIN Yunshan, et al. Recent development in the southern Yellow Sea sedimentology-The China-Korea Joint Investigation[R]. Marine Sciences, 1998

    (1):34-37.]

    [11] 赵一阳, 李凤业, DeMaster D J, 等. 南黄海沉积速率和沉积通量的初步研究[J]. 海洋与湖沼, 1991, 22(1):38-43.

    [ZHAO Yiyang, LI Fengye, DeMaster D J, et al. Preliminary studies on sedimentation rate and sediment flux of the south Huanghai Sea[J]. Oceanologia et Limnologia Sinica, 1991, 22(1):38-43.]

    [12]

    Alexander C R, DeMaster D J, Nittrouer C A. Sediment accumulation in a modern epicontinental-shelf setting:The Yellow Sea[J]. Marine Geology, 1991, 98(1):51-72.

    [13] 李凤业, 高抒, 贾建军, 等. 黄渤海泥质沉积区现代沉积速率[J]. 海洋与湖沼, 2002, 33(4):364-369.

    [LI Fengye, GAO Shu, JIA Jianjun, et al. Contemporary deposition rates of fine-grained sediment in the Bohai and Yellow Seas[J]. Oceanologia et Limnologia Sinica, 2002, 33(4):364-369.]

    [14]

    Cheng Peng, Gao Shu, Park Y A, et al. Grain size characteristics of the central Yellow Sea mud deposits[J]. The Yellow Sea, 1999, 5:1-7.

    [15] 王琦, 杨作升. 黄海南部表层沉积中的自生黄铁矿[J]. 海洋与湖沼, 1981, 12(1):25-32.

    [WANG Qi, YANG Zuosheng. Authigenic pyrite in the surface sediments of the southern Huanghai Sea[J]. Oceanologia et Limnologia Sinica, 1981, 12(1):25-32.]

    [16] 陈庆. 南海沉积物中自生黄铁矿研究[J]. 地质学报, 1981, 55(3):56-61.

    [CHEN Qing. Study on authigenic pyrites in sediments of the south Huanghai Sea[J]. Acta Geologica Sinica, 1981, 55(3):56-61.]

    [17] 李安春, 陈丽蓉, 申顺喜. 南黄海H-106岩柱中自生黄铁矿的硫同位素研究[J]. 科学通报, 1991, 36(12):928-930.

    [LI Anchun, CHEN Lirong, SHEN Shunxi. The sulfur isotopes in authigenic pyrite from the south Huanghai Sea[J]. Chinese Science Bulletin, 1991, 36(18):1549-1552.]

    [18] 庄丽华, 常凤鸣, 李铁刚, 等. 南黄海EY02-2孔底栖有孔虫群落特征与全新世沉积速率[J]. 海洋地质与第四纪地质, 2002, 22(4):7-14.

    [ZHUANG Lihua, CHANG Fengming, LI Tiegang, et al. Foraminiferal faunas and Holocene sedimentation rates of core EY02-2 in the South Yellow Sea[J]. Marine Geology and Quaternary Geology, 2002, 22(4):7-14.]

    [19]

    Kong G S, Park S C, Han H C, et al. Late Quaternary paleoenvironmental changes in the southeastern Yellow Sea, Korea[J]. Quaternary International, 2006, 144(1):38-52.

    [20] 李广雪, 杨子赓, 刘勇. 中国东部海域海底沉积物类型图[Z]. 北京:科学出版社, 2005.[LI Guangxue, YANG Zigeng, LIU Yong. Sediment Distribution Map of the East China Seas[Z]. Beijing:Science Press, 2005.]
    [21]

    Biscaye P E. Mineralogy and sedimentation of recent deep-sea clay in the Atlantic Ocean and adjacent seas and oceans[J]. Bulletin of the Geological Society of America, 1965, 76(7):803-831.

    [22]

    Petschick R, Kuhn G, Gingele F. Clay mineral distribution in surface sediments of the South Atlantic:sources, transport, and relation to oceanography[J]. Marine Geology, 1996, 130(3):203-229.

    [23]

    Stuiver M, Reimer P J. Extended 14C database and revised CALIB radiocarbon calibration program[J]. Radiocarbon, 1993, 35(1):215-230.

    [24]

    Hughen K A, Baillie M G L, Bard E, et al. Marine radiocarbon age calibration, 0~26 cal.kyrBP[J]. Radiocarbon, 2004, 46(3):1059-1086.

    [25]

    Schöe B R, Oschmann W, Tanabe K, et al. Holocene seasonal environmental trends at Tokyo Bay, Japan, reconstructed from bivalve mollusk shells:implications for changes in the East Asian monsoon and latitudinal shifts of the Polar Front[J]. Quaternary Science Reviews, 2004, 23(9-10):1137-1150.

    [26]

    LIU Jingpu, Milliman J D, GAO Shu, et al. Holocene development of the Yellow River's subaqueous delta, North Yellow Sea[J]. Marine Geology, 2004, 209(1-4):45-67.

    [27]

    Füchtbauer H, Müller G. Sedimente und Sedimentgesteine, Sediment-Petrologie Teil Ⅱ[M]. Stuttgart:Schweizerbart'sche Verlagsbuch-handlung, 1988:129-139.

    [28] 徐海, 洪业汤, 林庆华, 等. 红原泥炭纤维素氧同位素指示的距今6 ka温度变化[J]. 科学通报, 2002, 47(15):1181-1186.

    [XU Hai, HONG Yetang, LIN Qinghua, et al. Temperature variations in the past 6000 years inferred from δ18O of peat cellulose from Hongyuan, China[J]. Chinese Science Bulletin, 2002, 47(18):1578-1584.]

    [29] 洪业汤, 姜洪波, 陶发祥, 等. 近5 ka温度的金川泥炭δ18O记录[J]. 中国科学D辑, 1997, 27(6):525-530.

    [HONG Yetang, JIANG Hongbo, TAO Faxiang, et al. Temperature variations in the past 5000 years inferred from δ18O of peat cellulose from Jinchuan, China[J]. Science in China (Ser. D), 1997, 27(6):520-524.]

    [30] 姚檀栋, Thompson L G. 敦德冰心记录与过去5 ka温度变化[J]. 中国科学B辑, 1992, 10:1089-1093.[YAO Tandong, Thompson L G. The past 5

    ka years temperature variation in Dunde ice core[J]. Science in China (Ser. B), 1992, 10:1089-1093.]

    [31]

    Wang Yongji, Cheng Hai, Edwards R L, et al. The Holocene Asian monsoon:links to solar changes and North Atlantic climate[J]. Science, 2005, 308:854-857.

    [32]

    Kreisa R D. Storm-generated sedimentary structures in subtidal marine facies with examples from the Middle and Upper Ordovician of southwestern Virginia[J]. Journal of Sedimentary Research, 1981, 51(3):823-848.

    [33] 严钦尚. 论滨岸浅海的风暴沉积[J]. 海洋与湖沼, 1984, 15(1):14-19.

    [YAN Qinshang. Overview of the storm-generated deposits on near shore zone and open shelf[J]. Oceanologia et Limnologia Sinica, 1984, 15(1):14-19.]

    [34]

    Schulz H M, Schoner A, Emeis K C. Long-chain alkenone patterns in the Baltic Sea:an ocean-freshwater transition[J]. Geochimica et Cosmochimica Acta, 2000, 64(3):469-477.

    [35]

    CHU Guoqiang, SUN Qing, LI Shengqiang, et al. Long-chain alkenone distributions and temperature dependence in lacustrine surface sediments from China[J]. Geochimica et Cosmochimica Acta, 2005, 69(21):4985-5003.

    [36] 杨作升. 黄河、长江、珠江沉积物中黏土的矿物组合、化学特征及其与物源区气候环境的关系[J]. 海洋与湖沼, 1988, 19(4):336-346.

    [YANG Zuosheng. Mineralogical assemblages and chemical characteristics of clays from sediments of the Huanghe, Changjiang, Zhujiang rivers and their relationship to the climate environment in their sediment source areas[J]. Oceanologia et Limnologia Sinica, 1988, 19(4):336-346.]

    [37] 王昆山, 石学法, 林振宏. 南黄海和东海北部陆架重矿物组合分区及来源[J]. 海洋科学进展, 2003, 21(1):31-40.

    [WANG Kunshan, SHI Xuefa, LIN Zhenhong. Assemblages, provinces and provenances of heavy minerals on the shelf of the southern Yellow Sea and northern East China Sea[J]. Advances in Marine Science, 2003, 21(1):31-40.]

    [38] 向荣, 杨作升, Saito Y, 等. 济州岛西南泥质区近2300 a来环境敏感粒度组分记录的东亚冬季风变化[J]. 中国科学D辑, 2006, 36(7):654-662.

    [XIANG Rong, YANG Zuosheng, Saito Y, et al. East Asia Winter Monsoon changes inferred from environmentally sensitive grain-size component records during the last 2300 years in mud area southwest off Cheju Island, ECS[J]. Science in China (Ser. D), 2006, 49(6):604-614.]

  • 期刊类型引用(3)

    1. 许晓燕,窦衍光,张勇,陈晓辉,孙治雷,白凤龙,翟滨,王利波,邹亮. 冲绳海槽中南部柱状沉积物氧化还原敏感元素赋存机理与环境指示意义. 海洋学报. 2022(10): 63-79 . 百度学术
    2. 王飞,丛柏林,张朝晖,杨黄浩. 南极罗斯海沉积物中可培养菌株的分离及胞外水解酶活性检测. 极地研究. 2020(01): 68-81 . 百度学术
    3. XIU Chun,DU Ming,ZHANG Xu,HUO Suxia. Changes of marine productivity and sedimentary environment recorded by biogenic components in the Antarctica Ross Sea since the last deglaciation. Journal of Oceanology and Limnology. 2020(06): 1746-1754 . 必应学术

    其他类型引用(2)

计量
  • 文章访问数:  1719
  • HTML全文浏览量:  216
  • PDF下载量:  22
  • 被引次数: 5
出版历程
  • 收稿日期:  2009-02-25
  • 修回日期:  2009-05-03

目录

/

返回文章
返回