广东三水盆地玄武岩源区特征与南海早期演化

张煜, 方念乔

张煜, 方念乔. 广东三水盆地玄武岩源区特征与南海早期演化[J]. 海洋地质与第四纪地质, 2021, 41(3): 95-113. DOI: 10.16562/j.cnki.0256-1492.2020092902
引用本文: 张煜, 方念乔. 广东三水盆地玄武岩源区特征与南海早期演化[J]. 海洋地质与第四纪地质, 2021, 41(3): 95-113. DOI: 10.16562/j.cnki.0256-1492.2020092902
shu
ZHANG Yu, FANG Nianqiao. Source characteristics of basalts in Sanshui Basin and the early tectonic evolution stage of the South China Sea[J]. Marine Geology & Quaternary Geology, 2021, 41(3): 95-113. DOI: 10.16562/j.cnki.0256-1492.2020092902
Citation: ZHANG Yu, FANG Nianqiao. Source characteristics of basalts in Sanshui Basin and the early tectonic evolution stage of the South China Sea[J]. Marine Geology & Quaternary Geology, 2021, 41(3): 95-113. DOI: 10.16562/j.cnki.0256-1492.2020092902
shu

广东三水盆地玄武岩源区特征与南海早期演化

  • 中国地质大学海洋学院，北京 100083

基金项目: 国家自然科学基金面上项目“三水盆地及周边地区古近纪火山-沉积组合与南海早期演化”（41572207）
详细信息
    作者简介:

    张煜（1992—），男，博士研究生，研究方向为岩石大地构造，E-mail：zhangyu_cugb@foxmail.com

  • 中图分类号: P736.1

Source characteristics of basalts in Sanshui Basin and the early tectonic evolution stage of the South China Sea

  • School of Ocean Sciences, China University of Geosciences, Beijing 100083, China

摘要

    摘要
  • 摘要: 南海在扩张前是否经历了陆内裂谷阶段是南海成因研究中一个重要的问题。三水盆地位于南海北部陆缘，其新生代以来喷发的双峰式火山岩具备大陆裂谷的岩石组合特征。通过对其中玄武岩主微量元素分析认为三水盆地玄武岩可以分为亚碱性和碱性玄武岩系列，两者均显示出明显的Nb、Ta正异常，相对于大陆地壳具有较低的Th/Sc、La/Nb和U/Al×1000，陆壳混染程度低；首次对盆地内玄武岩进行40Ar-39Ar测年，结合前人年代学结果表明玄武质岩浆强烈喷发的时段为61～54 Ma，其中亚碱性玄武岩喷发时间（60 Ma）早于碱性玄武岩（56 Ma）；通过熔融柱模型反演得到亚碱性岩浆源区起止熔融温压分别为1 517 ℃（3.03 GPa）和1 471 ℃（2.25 GPa），深度为101～76 km，碱性岩浆源区起止熔融温压分别为1 555 ℃（3.33 GPa）和1 506 ℃（2.48 GPa），深度为110～84 km，整体为石榴石-尖晶石橄榄岩过渡区且呈逐渐变深的趋势。综合岩浆源区特征以及岩石组合特征认为三水盆地在古新世具备大陆裂谷特征。通过对比三水盆地与南海扩张期岩浆活动的分布时段及源区特征，发现三水盆地与南海扩张期岩浆活动时间分布存在较长间隔，深部过程差异较大，三水盆地岩浆活动与南海扩张并无直接因果联系。
    Abstract: It is an important issue in the study of the origin of the South China Sea (SCS) whether the SCS evolved in an intracontinental rift stage before its opening up. The Sanshui Basin is located on the northern margin of the SCS. The Cenozoic eruption of the bimodal volcanic suggest that the basin was a continental rift. The major and trace elements of the basalts indicate that the basalts in Sanshui Basin can be divided into sub-alkaline and alkaline basalts series. Both of them show obvious Nb and Ta positive anomalies. Compared with the continental crust, they have lower Th/Sc, La/Nb and U/Al×1000. 40Ar-39Ar age and previous chronology results suggests that the eruption of basaltic magma was intense in 61～54 Ma, while the eruption of sub-alkaline basalt (60 Ma) is earlier than alkaline basalt (56 Ma). Based on the mantle melting column model, the temperature and pressure of the subalkaline magma source range from 1 517℃ (3.03 GPa) to 1 471℃ (2.25 GPa), in a depth of 101～76 km, while the alkaline magma source range from 1 555℃ (3.33 GPa) to 1 506℃ (2.48 GPa) in a depth of 110～84 km. According to the source of magmatism and rock assemblage, it is concluded that the Sanshui Basin was a Paleocene continental rift. Comparing the temporal distribution and source magmatic activity in the spreading period between the Sanshui Basin and the South China Sea, it is inferred that there is a long time interval between the magmatic activity of the Sanshui Basin and the spreading of SCS, and the deep process is quite different. There is no direct relationship between magmatism in the Sanshui Basin and the spreading of SCS.
    • HTML全文

  • 天然气水合物（以下简称水合物）储量丰富，是最有可能代替化石燃料的清洁能源之一[1-3]。水合物沉积物层存在于低温、高压的自然环境中，陆地与海洋均有分布。中国海域水合物主要分布于南海地区[4]，据估算其储量约为6.4×1013 m3，与中国常规天然气资源总量基本相当[5]，具有巨大的开发潜力。随着中国对于海域水合物研究的不断深入，中国地质调查局分别在2017年和2020年进行了海域水合物的试采[6-7]，实现短期稳定产气。

    国内外学者针对含水合物沉积物的力学特性开展了丰富的力学试验研究[8-11]。Winters等[12-13]通过核磁共振和三轴压缩试验分别对原状样品和人工合成样品进行测量，得到了水合物沉积物的强度和压缩波速。Masui等[14]通过实验证实了原状样品和人工合成样品具有相似的力学强度，但是应力-应变关系随初始孔隙率和颗粒级配的变化而变化。Yun等[15]通过室内试验证明了四氢呋喃水合物（THF）和甲烷水合物具有相似的力学强度，建议使用THF进行力学性质研究。对于孔隙型水合物的研究较为广泛，Ren[16]等通过核磁共振发现水合物首先生长于小孔隙中，黏土会延缓整个生长过程。相较于均匀分布的含水合物沉积物样品，非均匀分布会对沉积物的物理和力学性质产生影响。Xie等[17-18]考虑了水合物沉积物的非均匀性，在人工样品中加入不同形状的水合物块体模拟不同的赋存形式，得到了不同形式下的水合物沉积物的强度和描述模型。含水合物沉积物的渗透率和破裂压力等均与水合物饱和度及分布形式有关[19-20]。然而大规模商业化开采水合物需要中长期稳定产气，因此还需要考虑储层的蠕变特性[21]。

    目前，针对粉黏土蠕变特性的研究已经较为成熟[22-24]。余云燕等[25]针对非饱和盐渍土开展了一系列三轴固结不排水蠕变试验，在Singh-Mitchell模型的基础上进行修正，提出了适合非饱和盐渍土的经验模型。罗庆姿等[26]通过三轴不排水压缩试验对汕头的软土进行了分析。Lai等[27]针对三峡滑土带中的非饱和土进行三轴蠕变试验，引入了考虑基质吸力的相关参数，建立了适用于非饱和土的Singh-Mitchell蠕变模型。刘业科等[28]通过梯级载荷加载方法进行室内三轴固结不排水的蠕变试验，在此基础上对Singh-Mitchell模型进行修正使其能描述零应力零应变水平。关于常规粉黏土的蠕变特性试验已经开展了较多研究，积累了较多经验[29]。区别于常规软黏土的蠕变特性，含水合物粉质黏土在开采过程中存在水合物相变、渗流等过程，水合物分解前后沉积物的孔隙结构也相应发生改变，目前针对含水合物粉黏土的蠕变特性研究和相关模型的研究较少[30]。

    本工作采用固结仪进行蠕变试验，研究了含水合物粉黏土的蠕变特性，并提出了描述模型。首先，通过试验获得了含水合物粉黏土的蠕变特性曲线，分析了含水合物粉黏土分解后的蠕变特性；在Singh-Mitchell模型中引入水合物饱和度参数，构建并验证修正后的Singh-Mitchell模型，得到了可以描述含水合物粉黏土的蠕变特性模型。

    1.   区域概况

    研究区域位于中国南海北部的西北大陆边缘，呈NE向延展，面积达4.5×104 km2，具备稳定的高压低温环境，可以形成天然气水合物，也是我国重要的常规石油富集区域[31-32](图1)。试验采用研究区海域水合物储层的粉黏土沉积物，水合物赋存于海底7～158 m范围内的第四系沉积物中，该区域平均水深约为1 000 m。

    图  1  研究区域位置[31]
    Figure  1.  Location of the research area[31]
    下载: 全尺寸图片 幻灯片

    2   试验样品来源及方法

    2.1   试验材料

    土体的基本物理和力学性质如下：土体颗粒比重Gs=2.65，孔隙度n=0.509，平均粒径d50=0.0067，其中黏质组分（d≤4 μm）约占50%，砂质组分（4 μm＜d≤63 μm）约占45%，颗粒级配曲线如图2所示。根据《土的工程分类标准》，试验用土粒径小于0.075 mm的颗粒质量分数大于50%，属于细粒土；不均匀系数Cu = 13.56，属于不均匀土，曲率系数Cc = 0.57，级配不连续。综上，试验土为级配不良的细粒土。

    图  2  沉积物颗粒级配曲线图
    Figure  2.  Curves of sediment particle size distribution
    下载: 全尺寸图片 幻灯片

    试验所用的四氢呋喃溶液采购于中国国药试剂网，浓度为99.9%，通过加入蒸馏水配置成不同浓度的四氢呋喃溶液。质量分数为19%的四氢呋喃溶液在低温常压条件下可以完全合成四氢呋喃水合物，考虑到四氢呋喃的挥发性，试验中使用四氢呋喃溶液的质量分数为21%。根据试验样品所需的水合物饱和度，计算并配置一定质量分数的四氢呋喃溶液。

    2.2   试验样品制备与试验设置

    试验介质为粉黏土，基于试验区域水合物储层实际工况和土工试验结果，粉黏土比重Gs = 2.65，制备水合物沉积物样品时选取干密度$ {\rho }_{\mathrm{d}}\text{=}\text{1.3}\;\text{}\text{g/}{\text{cm}}^{\text{3}} $，研究在恒定的上覆层荷载下水合物储层分解后对水合物储层的蠕变力学特性的影响。具体制备过程如下：

    土骨架制备：将洗盐后的粉黏土沉积物置于烘箱中在105℃条件下烘干12 h后取出冷却，研磨土体至颗粒完全粉碎（土体颗粒粒径小于0.2 mm）。根据干密度称取一定质量的沉积物平均二等分，将每份土体加入适量水后填入环刀内砸实至所要求的高度，在每层样品制备完成后将试样拉毛，以降低其分层的影响，制备成面积为30 cm2、高2 cm样品，置于烘箱中在105℃条件下烘干12 h后取出冷却。

    水合物合成：根据试验中水合物饱和度的设置配置不同浓度的四氢呋喃溶液，将制备好的土骨架置于密封罐内，密封罐与真空装置和储液装置相连，开启抽真空装置对样品抽取真空两小时，完成后通过储液装置自下而上浇注配置好的四氢呋喃溶液，保持密封状态一天至土样含液量稳定。最后取出浸泡完成的沉积物样品置于0～2 ℃的冷柜中1～2 d合成水合物，如图3所示。

    图  3  试验前后的试样
    Figure  3.  Sample before and after test
    下载: 全尺寸图片 幻灯片

    在试验前取出含水合物沉积物样品置于密封罐内分解，从而得到不同水合物饱和度分解后的沉积物样品。采用分别加载的方式进行，在荷载的作用下，连续观测。加载应力分别为：①号0.5 MPa、②号2 MPa、③号3 MPa条件下进行试验；分解前水合物饱和度分别为 0%、20%、60%、80%。试验仪器使用固结仪如图4，通过采用单杠杆固结仪试验获得水合物沉积物的蠕变力学响应参数。

    图  4  高压固结仪实物图
    Figure  4.  High-pressure consolidation apparatus
    下载: 全尺寸图片 幻灯片

    高压固结仪放置于20℃的恒温室中，在试验过程中保持样品处于饱和状态，当沉降值稳定后样品停止变形，即认为试样的变形不超过0.001 mm/d后达到稳定状态。

    3.   结果和分析

    3.1   试验结果

    通过对原始数据的整理，得到了分别加载方式下，不同水合物饱和度（SH）粉质黏土样品分解后在轴向压力σ下试验时间t和应变ε的蠕变试验曲线（图5）。

    图  5  不同饱和度下蠕变全过程曲线图
    Figure  5.  The whole-process curve of creeping under different saturation degrees
    下载: 全尺寸图片 幻灯片

    由图5中可以看出，在分别加载的条件下，含水合物沉积物分解后样品在各级荷载下的应变-时间关系曲线的变化规律基本相同，只有初始蠕变阶段和稳定蠕变阶段，并未出现加速蠕变阶段。施加荷载的初期，均产生了较大的瞬间变形，曲线呈飞跃式增长。随着时间的推移，曲线的斜率逐渐减小，表示试样的变形速率逐渐减小。最终，试样的变形速率趋于稳定，曲线呈斜率较小的直线，试样的变形随时间不断增大，但该阶段变形值较小。

    含水合物粉质黏土样品的蠕变特性曲线可以分为瞬时变形、固结变形和蠕变变形3个阶段[33]。水合物的赋存改变了粉黏土的初始孔隙结构，在施加荷载的瞬间会产生较大的变形进入瞬时变形阶段；当样品孔隙被压缩，孔隙水被排出，而由于黏土沉积物骨架的渗透率较低因此会产生超静孔隙水压力并需要较长时间消散，处于固结变形阶段；最后固结完成后，在恒定荷载的作用下，由于恒定轴向压力的作用，土体骨架持续变形，应变不断增加但总体应变率较小，最后趋于稳定，处于蠕变变形阶段（图6）。无水合物沉积物应变率和时间的曲线见图7。瞬时变形认为当轴向荷载加载完成后完成，取1 min对应的应变作为瞬时变形。将图7中的转折点作为固结变形和蠕变变形分界点，即将应变率1×10−5/min作为固结变形和蠕变变形分界点。

    图  6  变形过程三阶段示意图
    Figure  6.  Diagram of three-stage deformation process
    下载: 全尺寸图片 幻灯片
    图  7  无水合物沉积物应变率和时间的曲线
    Figure  7.  Curves of strain rate and time for hydrate-free sediments
    下载: 全尺寸图片 幻灯片

    3.2   Singh-Mitchell蠕变模型

    Singh和Mitchell等[34]在总结大量土体蠕变的基础上，使用指数函数及幂函数的形式描述黏土的应力-应变-时间关系，并在此基础上提出了Singh-Mitchell蠕变模型。Singh-Mitchell蠕变模型具有参数少、适用范围广等特点，因此本文先择此模型来描述粉黏土的蠕变特性。Singh-Mitchell蠕变模型可以表示为：

    $$ \dot{{ \varepsilon }}={{A}{\rm e}}^{\text{α}{{D}}_{\text{r}}}{\left(\frac{{t}}{{{t}}_{{1}}}\right)}^{{m}} $$ (1)

    式(1)中：$ \dot{{ \varepsilon }} $表示任意时刻t的应变速率；A、α、m为模型参数，m为$ \text{ln}\text{}{ \varepsilon }\text{}\text{-}\text{}\text{ln}{}{t} $关系图中的直线斜率的绝对值，一般m值为0.75～1.0，并且$ \text{m}\text{}\text{≠}\text{}\text{1} $，$ {{D}}_{\text{r}} $为偏应力水平；$ {{t}}_{\text{1}} $为单位参考时间；t为蠕变加载时间。式中需要确定A、α、m三个参数。将式(1)两边分别积分，可以得到式(2)：

    $$ { \varepsilon }{}={{ \varepsilon }}_{{0}}+{}\frac{{{A}{{t}}_{{1}}}^{{m}}}{{1-}{m}}\text{exp}{{(}\overline{\alpha}}{{D}}_{\text{r}})\left(\frac{t}{t_1}\right)^{{1-}{m}} $$ (2)

    式(2)中$ { \varepsilon } $为轴向应变，$ {{ \varepsilon }}_{\text{0}} $为初始轴向应变；在不考虑初始的轴向应变时，$ {{ \varepsilon }}_{\text{0}}\text{}\text{=}\text{}\text{0} $，令$ \lambda = 1- m $，$ {B}\text{}\text{=}\text{}\dfrac{{{{At}}_{{1}}}^{{m}}}{{1-}{m}} $，$ \beta=\overline{\alpha} $，得到式(3)：

    $$ { \varepsilon }={B}\text{exp(}\beta{{D}}_{\text{r}})\left(\dfrac{{t}}{{{t}}_{{1}}}\right)^{{ \lambda }} $$ (3)

    式(3)即蠕变方程，其中需要确定的参数有3个分别为$ {B}\text{、}\beta\text{、}{ \lambda } $。

    当${t}={t}_{1} $时，有$ { \varepsilon }={ \varepsilon}_{1} $，可以得到方程(4)，两边取对数得到方程(5)：

    $$ { \varepsilon }={B}\text{exp}\left(\beta{{D}}_{\text{r}}\right) $$ (4)
    $$ \text{ln}\text{}{ \varepsilon }=\beta{{D}}_{\text{r}}+\text{ln}{B} $$ (5)

    β和B的值可以直接通过单位参考时间$ {{t}}_{\text{1}} $的$ \text{ln}\text{}{ \varepsilon }\text{}\text{-}\text{}{{D}}_{\text{r}} $直线关系图得到。土体的非线性蠕变是一簇曲线，从上式中可以看出Singh-Mitchell蠕变模型也是一簇指数形式的曲线，形状与现有的粉黏土蠕变试验结果吻合，因此选择Singh-Mitchell蠕变模型来描述粉黏土的蠕变特性是合适的。

    3.3   试验用粉黏土蠕变特性分析

    将蠕变特性曲线绘制在双对数坐标轴中，可以得到$ \text{ln}\text{}{ \varepsilon }\text{}\text{-}\text{}\text{ln}\text{}{t} $关系图（图8），图中可以看出，在不同的剪切应力水平下样品的蠕变试验曲线为一组平行近似的直线，直线的斜率为模型中${ \lambda } $值。将$ {{t}}_{\text{1}}\text{=1}\text{}\text{min} $时的应变带入式(4)和式(5)中得到$ \text{ln}\text{}{{ \varepsilon }}_{{\rm r}}\text{}\text{-}\text{}{{D}}_{\text{r}} $关系曲线，拟合后具有良好的线性关系，通过式(5)可以计算出不含水合物沉积物的β和B值（图9）。

    图  8  无水合物试样$ \text{lg}{ \varepsilon }\;\text{vs lg}{t} $蠕变特性曲线
    Figure  8.  The lgε vs lgt creep characteristic curve of hydrate-free samples
    下载: 全尺寸图片 幻灯片
    图  9  $ \text{ln}{{ \varepsilon }}_{\text{r}}\text{}\;\text{vs}\;\text{}{{D}}_{\text{r}} $关系拟合曲线
    Figure  9.  Fitting curve $ \text{ln}{{ \varepsilon }}_{\text{r}}\;\text{vs}\;\text{}{{D}}_{\text{r}} $ relationship
    下载: 全尺寸图片 幻灯片

    考虑含水合物沉积物的剪切应力，偏应力的计算公式为q$ \text{=}\text{}\left(\text{1-}{{K}}_{\text{0}}\right)\text{*}\text{σ} $，其中$ {{K}}_{\text{0}} $为静止土压力系数，σ为轴向应力。根据Simpson等[35]正常固结黏土$ {K}_{0} $经验公式：

    $$ {{K}}_{\text{0}}\text{}\text=\text{}\frac{{1-}\dfrac{\text{1}}{\sqrt{\text{2}}}{\sin}\varphi'}{\text{1+}\dfrac{\text{1}}{\sqrt{\text{2}}}{\sin}\varphi'} $$ (6)

    将三轴剪切试验内摩擦角$\varphi '=5^{\text °} $带入式(6)后得到K0=0.84，即q=0.16σ，同时注意到不同围压下qf=0.15~0.2 MPa[36]，$ \dfrac{{1-}{{K}}_{{0}}}{{{q}}_{\text{f}}}\text{≈ 1} $，因此后续取值Dr=σ。

    带入上述参数计算Singh-Mitchell蠕变模型中的相关参数值，结果如表1所示。

    表  1  Singh-Mitchell蠕变模型中参数计算值
    Table  1.  Parameter values of the Singh-Mitchell creeping model
    水合物饱和度σ / MPaDrBβλTr / min
    00.50.513.30.2060.01251
    2.02.013.30.2060.01011
    3.03.013.30.2060.00811
    下载: 导出CSV 
    | 显示表格

    将表1中计算的模型参数带入式(4)中，以轴向压力σ为0.5 MPa为例得到Singh-Mitchell蠕变模型的表达式：

    $$ { \varepsilon }\text{}\text=\text{}\text{13.3}{\text{e}}^{\text{0.103}}{{t}}^{\text{0.0}\text{125}} $$ (7)

    将模型计算曲线和试验值进行对比（图10）。可以看出，Singh-Mitchell蠕变模型计算曲线和蠕变试验曲线的发展趋势具有较好的一致性，在加载瞬间应变快速增长，随着时间的增加轴向应变迅速衰减。在较高的剪切应力水平下Singh-Mitchell蠕变模型和试验曲线具有较好的一致性，在较低的剪切应力水平下模型计算值略大于试验值，有一定的误差，总体拟合关系较好。

    图  10  无水合物模型计算曲线和试验值对比
    Figure  10.  Comparison of model calculation curve to the test result of hydrate-free sediments
    下载: 全尺寸图片 幻灯片

    3.4   含水合物试样蠕变特性分析

    上一节中讨论了不含水合物沉积物样品的蠕变特性曲线，通过结果可以看出Singh-Mitchell蠕变模型可以描述其蠕变特性，为了使Singh-Mitchell蠕变模型可以描述含水合物沉积物样品，需要对模型进行修正。引入水合物饱和度SH来描述水合物含量，SH是指水合物与孔隙体积之比[37]：

    $$ {{S}}_{\text{H}}\text{}\text=\text{}\frac{{{V}}_{{{S}}_{\text{H}}}}{{{V}}_{\text{v}}} $$ (8)

    水合物在合成时会膨胀从而引起土体孔隙的变化，在经典模型中，模型参数A反映出土的组成、结构和应力历史等特性，因此考虑通过对参数A进行修正来表征水合物饱和度的影响，在Singh-Mitchell蠕变模型的基础上引入模型如式(9)：

    $$ { \varepsilon }\text={{ \varepsilon }}_{\mathrm{r}}{\text{(1+}{{S}}_{\text{H}}\text{)}}^{\textit{α}} $$ (9)

    式中，α为模型参数，根据试验数据，本文中取0.37。

    将修正后的参数A和不同水合物饱和度经典模型中的参数A进行对比（图11），可以看出其具有良好的对应关系。

    图  11  经典模型和修正模型参数$ \text{A} $对比图
    Figure  11.  Comparison of classical model and modified model parameter A
    下载: 全尺寸图片 幻灯片

    由此建立了一种修正的Singh-Mitchell蠕变模型：

    $$ \dot{{ \varepsilon }\text{}}\text=\text{}{{{A}}_{\text{0}}\text{*}{\left(\text{1+}{{S}}_{\text{H}}\right)}^{\textit{α}}{e}}^{\textit{β}{{D}}_{\text{r}}}{\left(\frac{{t}}{{{t}}_{\text{1}}}\right)}^{{ \lambda }} $$ (10)

    当SH=0时，模型又退化成经典的Singh-Mitchell蠕变模型。

    通过计算，得到不同工况下修正的Singh-Mitchell蠕变模型参数，得到参数（表2）。

    表  2  修正的 Singh-Mitchell 蠕变模型参数值
    Table  2.  Parameter values of modified Singh-Mitchell creeping model
    水合物饱和度/ %σ / MPaDrA(1+SH)αβλ
    200.50.514.230.2000.0123
    2.02.014.230.2000.0160
    3.03.014.230.2000.0114
    600.50.515.830.1920.0061
    2.02.015.830.1920.0117
    3.03.015.830.1920.0074
    800.50.516.530.2260.0076
    2.02.016.530.2260.0092
    3.03.016.530.2260.0064
    下载: 导出CSV 
    | 显示表格

    将表2中的参数带入公式(10)中，得到模型计算值，同时和试验数据对比如图12所示。

    图  12  不同水合物饱和度模型计算曲线和试验值对比
    Figure  12.  Comparison of model calculated curves and the test values at different hydrate saturation degrees
    下载: 全尺寸图片 幻灯片

    从图12可以看出，修正后的Singh-Mitchell蠕变模型在不同水合物饱和度和剪切应力下都能够较好地预测均匀分布的含水合物粉质黏土样品的蠕变特性，因此可以使用此修正的Singh-Mitchell蠕变模型来预测含水合物沉积物样品的蠕变特性曲线。修正的Singh-Mitchell蠕变模型参数计算简便、参数较少并且容易获取，使用较为方便，因此本文建立的修正Singh-Mitchell蠕变模型对描述含水合物沉积物样品的蠕变特性曲线更为合适，对我国南海水合物储层长期开采后的沉降预测具有一定的指导意义。

    4.   讨论

    4.1   蠕变试验结果分析

    南海海域水合物主要赋存于海底以下到200 m范围内的粉黏土等未固结的沉积物内，具有含水高、孔隙度高的特点，水合物饱和度最高可达90%以上[38]。因此本文根据储层的实际赋存条件，以水合物储层的粉黏土为介质，在0.5、2和3 MPa的压缩加载条件下进行了0、20%、60%、80%水合物饱和度分解后的蠕变特性试验。试验表明：加载过程中，含水合物沉积物经历瞬时变形、固结变形和蠕变变形3个阶段；在加载的短时间内，瞬时变形阶段完成，引起的轴向变形约占总变形量的70%～80%，最大值为85.8%，在固结变形阶段，孔隙水在稳定的应力条件下从孔隙中排出，一般数小时内完成，引起的轴向应变占应变量的15%～25%，最后的蠕变变形阶段一般数天内完成，引起的轴向应变小于总应变量的5%。轴向压力和水合物饱和度均对轴向变形有影响。轴向压力越大，蠕变变形越大，轴向压力在较低水平时变化对轴向变形的影响大，随着轴向压力的增加，这种影响程度减小。水合物饱和度主要通过改变样品的微观孔隙结构进而影响蠕变特性，分解前的水合物饱和度越高，蠕变变形越大。水合物饱和度在低轴压条件下对变形影响小，在轴压为0.5 MPa时，从无水合物到水合物饱和度增加到80%时，轴向变形从17.25%增加到18.8%，而在轴压为3 MPa时，轴向变形从26.83%增加到了34.73%。

    4.2   基于Singh-Mitchell的修正模型

    数学模型是定量描述各个试验变量之间相互关系的常用方法。之前的海洋地质调查海域水合物储层粉质黏土进行蠕变特性的研究不足，本文利用蠕变试验结果，通过建立修正的Singh-Mitchell蠕变模型，对研究区域内水合物储层土的蠕变特性进行初步分析。从蠕变模型中可以看出，通过轴向压力和水合物饱和度两个物理指标能够建立出较好的蠕变模型，反映出这两个物理指标是蠕变特性的主要影响因素。模型中参数α和土体初始条件有关，主要根据试验获取，反映出土体的蠕变特性与土体孔隙微观结构、颗粒级配、应力历史等多种因素相关。

    5.   结论

    （1）含水合物沉积物样品蠕变特性曲线可以分为瞬时变形、固结变形和蠕变变形3个阶段。瞬时变形数分钟内完成，引起较大的变形，固结变形数小时内完成，引起较小的变形，蠕变变形数天内完成，引起极小的变形。

    （2）蠕变变形随着加载应力水平、水合物饱和度的增大而增加，水合物饱和度主要通过改变沉积物样品的孔隙结构从而影响蠕变变形，水合物饱和度越高，孔隙结构受影响程度也越大，试样的物理和力学性质越差。在低应力水平时，水合物饱和度对蠕变变形的影响不明显。

    （3）基于应力水平和水合物饱和度建立的修正Singh-Mitchell蠕变模型，可以有效地预测含水合物沉积物的蠕变特性，对水合物储层长期开采后的沉降变形预测具有指导意义，便于实际工程中推广。但由于试验组数少等原因本文未对水合物沉积物蠕变变形机理进行深入研究，针对非均匀分布的含水合物沉积物蠕变特性也需要进一步研究。

  • 图  1   三水盆地位置及地质简图

    a.南海北部陆域火山岩盆地分布（据文献[13]进行修改），b.三水盆地地质简图。

    Figure  1.   Location and geological map of Sanshui Basin

    a. the distribution of volcanic basins to the northern part of South China Sea （modified from reference [13]）, b. the geological map of Sanshui Basin.

    下载: 全尺寸图片 幻灯片

    图  2   三水盆地玄武岩TAS图解

    橙色线为碱性-亚碱性界线，据文献[28]。

    Figure  2.   TAS diagram of basalt in Sanshui Basin

    The orange line is the alkaline-subalkaline boundary, according to reference[28].

    下载: 全尺寸图片 幻灯片

    图  3   三水盆地玄武岩稀土及微量元素分布

    a.三水盆地玄武岩球粒陨石标准化稀土元素配分图，b.原始地幔标准化微量元素蛛网图；原始地幔以及球粒陨石数据来自文献[28]，LCC（下地壳）和UCC（上地壳）数据来自文献[29]，U1500（约33 Ma）为南海扩张初期洋壳样品，数据来自文献[30]；a、b图例相同，ZD、WJG、DQC及STC含义同表1

    Figure  3.   Distribution of rare earth and trace elements in basalts in Sanshui Basin

    a.trace element pattern, b.REE pattern of the basalt; data are normalized to primitive mantle and chondrite of [28], LCC （lower continental crust） and UCC （upper continental crust） is from [29], U1500 （about 33 Ma） are the ocean crust of SCS from [30]; the legends of a and b are the same, and the meanings of ZD, WJG, DQC and STC are same as those in Table 1.

    下载: 全尺寸图片 幻灯片

    图  4   三水盆地玄武岩40Ar/39Ar测年结果

    Figure  4.   40Ar/39Ar dating results of basalt in Sanshui Basin

    下载: 全尺寸图片 幻灯片

    图  5   三水盆地玄武岩Haker图解

    Figure  5.   Haker diagram for basalts in Sanshui Basin

    下载: 全尺寸图片 幻灯片

    图  6   三水盆地亚碱性-碱性玄武岩FeOT、Na2O含量回归校正

    a.FeOT回归校正，b.Na2O回归校正；ZD-紫洞，WJG-王借岗。

    Figure  6.   FeOT and Na2O regression correction of sub-alk （ZD） and alk （WJG） basalts

    a.regression correction of FeOT, b. regression correction of Na2O; ZD-ZiDong, WJG-WangJieGang.

    下载: 全尺寸图片 幻灯片

    图  7   亚碱性（ZD）-碱性（WJG）玄武岩熔融柱模型FeO-Na2O关系图

    P0表示部分熔融作用发生时的压力，Pf表示部分熔融作用结束时的压力；ZD-紫洞，WJG-王借岗。

    Figure  7.   FeO-Na2O diagram of sub-alk （ZD） and alk （WJG） mantle melting column

    P0 represents the pressure when partial melting occurs, and Pf represents the pressure at the end of partial melting; ZD-ZiDong, WJG-WangJieGang.

    下载: 全尺寸图片 幻灯片

    图  8   三水盆地源区性质(La/Yb)N-(Sm/Yb)N图解

    Figure  8.   Mantle source （La/Yb）N-（Sm/Yb）N diagram of Sanshui Basin

    下载: 全尺寸图片 幻灯片

    图  9   三水盆地火山岩年龄分布

    年龄数据来自[16]，[17]及[21]。

    Figure  9.   Age distribution of volcanic rocks in Sanshui Basin

    Age datas are from [16], [17] and [21].

    下载: 全尺寸图片 幻灯片

    图  10   三水盆地演化模式及其与南海演化的关系

    Figure  10.   The evolution model of the Sanshui Basin and its relationship with the evolution of the South China Sea

    下载: 全尺寸图片 幻灯片

    图  11   三水盆地玄武岩与OIB、MORB、玳瑁海山以及U1500样品源区对比

    彩色箭头为不同区域样品熔融柱模型计算结果，箭尾和箭头分别为部分熔融作用开始和停止时的的源区特征；彩色区块为前人对不同区域样品源区范围的计算结果，其中夏威夷OIB、东太平洋海岭及中大西洋洋脊MORB源区据文献[52]；灰色实线（近垂直）为固体地幔在不发生熔融条件下绝热上升线，据文献[52]；灰色虚线为熔融分数等值线，二辉橄榄岩固相线及熔融分数等值线据文献[53]；F-部分熔融分数；U1500 MORB源区计算结果据文献[30]。

    Figure  11.   Comparison of basalt source in Sanshui Basin with OIB, MORB, Daimao seamount and U1500 samples.

    The colored arrows are the calculation results of the melting column model for samples in different regions, the arrow tails and arrows represents the source area characteristics when partial melting starts and stops respectively; the colored blocks are the calculation results of the source area of the samples in different areas, while Hawaii OIB, the Eastern Pacific Rise and the Mid-Atlantic Ridge MORB source areas are based on [52]; the solid gray line （near vertical） represents the solid mantle adiabats, according to [52]; the gray dashed line is the melt fraction isopleths, according to [53]; F-partial melting fraction; the source of U1500 MORB according to [30].

    下载: 全尺寸图片 幻灯片

    表  1   三水盆地玄武岩野外及镜下特征

    Table  1   Field and microscopic characteristics of basalts in Sanshui Basin

    采样地区野外特征镜下结构主要矿物及含量矿物特征
    紫洞（ZD）灰黑色、块状构造，露头呈现明显的柱状节理间粒-间隐结构橄榄石（1%～3%）、单斜辉石（10%～12%）、斜长石（15%～19%）橄榄石无色（0.2～0.3 mm），他形粒状包裹于辉石颗粒中；辉石（2 mm）具筛状熔蚀特征，自形短柱状；斜长石（2～6 mm），发育聚片双晶，可见环带结构。
    王借岗（WJG）黑色、块状构造间隐结构橄榄石（3%～5%）、辉石（10%～15%）、斜长石（20%～36%）橄榄石呈椭圆形粒状（1.5～3 mm）。部分细粒橄榄石被具有环带的辉石包裹在核部。单斜辉石裂隙发育，多见呈八边形的横切面。
    邓群村（DQC）黑色、块状构造间隐结构辉石（5%～10%）、斜长石（约20%）辉石（1 mm）分为两类，一类为未发生变质的新鲜辉石，呈他形，另一类为表面发生滑石化的辉石。
    石头村（STC）黑色、块状构造间粒-间隐结构辉石（约15%）、斜长石（约15%）单斜辉石多见呈八边形的横切面，亦可见发育裂纹，发育较弱的环带。斜长石自形板状，发育清晰、完整的环带结构。部分长石具有熔蚀结构。基质具间隐结构，主要为斜长石微晶。
    下载: 导出CSV

    表  2   三水盆地玄武岩主量元素测试及标准矿物计算结果

    Table  2   Major element and CIPW results of basalts in Sanshui Basin

    主量元素含量/（wt%） 标准矿物计算/%
    样品编号SiO2Al2O3TiO2Fe2O3Fe2O3TFeOCaOMgOK2ONa2OMnOP2O5FeOTQOrAbAnNeDiHyOlMtIlHmAp
    紫洞ZD-201-A51.6615.881.963.43−7.017.866.631.613.310.160.5010.100.309.4928.0023.740.009.6718.950.004.973.720.001.16
    ZD-20151.7716.201.953.39−6.767.966.331.643.350.160.509.820.499.6728.3124.350.009.5717.830.004.923.700.001.16
    14SSZD-N50.9815.842.03−11.43−7.965.851.713.570.150.4910.280.0010.1730.4022.310.0011.4114.410.006.273.890.001.15
    14SSZD-S251.3515.631.91−10.97−7.576.282.043.650.150.459.860.0012.1231.1120.380.0011.5811.192.746.183.650.001.05
    ZD-10152.4615.981.892.97−6.857.516.172.073.500.140.479.520.0012.2429.6221.780.0010.0616.231.114.313.580.001.08
    14SSZD-S152.4416.001.91−10.79−7.465.222.003.560.150.459.712.1511.8730.3521.930.009.9512.940.006.113.660.001.05
    ZD-20251.4416.162.016.31−4.348.016.041.653.350.160.5110.032.6610.0329.2024.870.009.9312.940.005.243.920.001.22
    ZD-20351.7016.052.005.53−4.967.685.811.983.630.150.529.940.7311.9231.3722.090.0010.6213.000.005.183.880.001.22
    18SS082-151.3116.252.020.0011.206.878.245.311.663.330.140.5310.081.309.8628.4224.670.0010.5414.880.005.213.870.001.25
    18SS082-251.8416.122.030.0010.827.948.175.151.803.380.150.539.740.0310.7428.8723.690.0011.2517.370.002.933.880.001.24
    王借岗WJG-20147.2616.782.855.35−4.838.547.171.764.510.150.819.64 0.0010.3926.6620.336.2513.050.008.277.765.420.001.87
    WJG-20347.8017.362.794.70−5.208.476.581.684.450.140.829.060.009.9228.9222.464.7211.130.008.836.825.300.001.91
    WJG-20447.0517.272.764.81−5.538.786.781.694.370.150.819.360.009.9825.7422.556.0712.410.009.176.975.250.001.87
    14SS012c46.5216.922.88−11.43−8.607.002.203.580.150.7210.280.0013.1322.3823.744.4311.490.0011.506.155.510.001.69
    14SS012b46.8616.652.82−11.04−8.977.092.203.510.150.729.930.0013.0822.1623.324.2213.320.0010.915.935.390.001.68
    14SS01246.8916.882.76−10.99−8.677.232.243.460.150.739.890.0013.3122.8024.093.6211.440.0011.855.915.280.001.71
    18SS083-347.7317.022.79−10.578.438.536.531.903.880.150.889.510.0011.3425.0923.624.3810.790.0015.601.765.360.002.07
    石头村14SS004-247.7417.402.79−12.11−9.705.081.712.830.150.4910.89 0.0010.1724.1329.980.0012.337.233.546.145.340.001.15
    18SS079-247.7417.482.81−11.798.049.825.411.372.910.150.5310.610.008.1424.8030.910.0012.046.556.764.185.390.001.23
    17SS060-147.7217.342.85−11.956.8310.005.291.342.850.170.5010.750.008.0024.3030.780.0012.759.751.446.375.460.001.17
    17SS060-247.9817.712.81−11.698.079.955.131.302.720.170.5310.520.007.7523.2332.570.0011.2311.623.003.985.390.001.23
    18SS079-147.9017.282.82−11.758.4110.015.231.362.960.150.5410.570.008.0925.2630.160.0013.445.107.773.515.420.001.26
    18SS081-847.9017.492.84−11.867.609.905.131.472.710.160.5410.670.008.7523.1431.490.0011.7411.082.124.985.450.001.26
    18SS081-748.1317.362.81−11.838.109.945.281.382.560.160.5510.640.008.2421.8432.100.0011.4814.481.074.135.390.001.28
    邓群村17SS063-248.9416.262.73−13.066.048.334.021.803.460.231.1611.75 2.7810.7329.4823.670.008.147.990.009.275.220.002.72
    19SS023-149.2515.832.83−13.676.527.734.121.613.510.261.1912.313.959.5629.9322.850.006.319.830.009.405.410.002.78
    17SS063-149.1116.332.70−12.907.358.553.951.643.370.201.2511.612.149.7928.7424.790.007.9611.580.006.925.180.002.91
    下载: 导出CSV

    表  3   三水盆地玄武岩微量元素测试结果

    Table  3   Trace elements of basalts in Sanshui Basin

    μg/g
    样品编号YLaCePrNdSmEuGdTbDyHoErTmYbLuLiBeScVCrCo
    紫洞ZD-201-A28.7330.1753.386.8828.276.662.176.241.116.051.102.670.462.880.6222.481.8823.08165.40217.6041.94
    ZD-20127.2629.0751.276.5527.066.502.125.951.065.651.062.540.442.770.6021.491.6521.50155.00192.0037.32
    14SSZD-N29.4028.1053.406.6127.406.232.096.231.005.591.062.810.412.420.36−1.6121.20159.00139.0042.80
    14SSZD-S228.2028.7054.006.5726.706.011.985.910.955.411.022.730.402.450.37−1.6720.10148.00184.0038.90
    ZD-10127.4629.5651.976.6127.056.342.095.981.065.571.052.580.432.670.6311.141.7522.47164.90199.3040.25
    14SSZD-S128.5029.3054.706.6927.206.112.016.090.975.441.042.780.412.440.37−1.5620.70151.00196.0040.80
    ZD-20228.8631.2656.356.7827.536.122.155.810.975.821.072.930.442.830.4719.191.6621.52153.10155.7041.29
    ZD-20329.1032.9958.527.0528.696.302.226.091.015.991.123.060.472.930.4810.121.6421.26152.70147.2039.92
    18SS082-130.0027.8051.606.6427.206.442.255.000.995.171.072.400.462.930.3810.501.3522.10177.00165.0062.30
    18SS082-230.6028.1053.606.6927.306.172.175.280.994.901.072.420.472.840.377.571.7121.90185.00173.0060.50
    王借岗WJG-20128.2641.8976.199.0735.987.312.536.721.045.911.072.890.432.660.4228.131.5823.85211.20154.4041.93
    WJG-20329.7042.8677.979.2136.707.432.556.771.045.941.092.980.432.740.4435.391.6822.69206.40116.2043.74
    WJG-20426.4639.9671.278.5233.696.762.336.290.965.470.992.720.402.520.4132.641.5820.69191.00109.7040.45
    14SS012c27.1037.6070.708.6034.707.062.396.691.005.311.002.610.372.070.32−1.5621.00211.00116.0045.60
    14SS012b28.5039.2074.209.0136.407.512.507.001.065.651.052.800.392.280.34−1.5123.40217.00153.0043.50
    14SS01226.1037.2070.008.4834.207.022.386.540.995.230.972.530.352.010.30−1.5020.70197.00119.0038.20
    18SS083-329.5037.1070.209.1035.707.272.425.761.074.991.092.290.382.520.3725.001.5320.60215.00113.0063.10
    石头村14SS004-224.2024.0049.706.4827.406.192.165.960.925.010.912.400.341.980.29−1.2524.30240.0067.3042.00
    18SS079-225.1023.1046.306.2025.906.092.204.630.914.460.952.090.352.290.3020.001.2024.90256.0066.7055.30
    17SS060-129.2026.2051.806.8429.606.472.176.081.146.001.122.900.462.720.3818.501.3825.40284.0069.3040.40
    17SS060-226.4023.9050.006.2628.806.362.445.030.974.921.012.120.362.220.3021.701.2624.40272.0068.5081.10
    18SS079-125.1022.8048.406.0927.306.342.134.650.934.711.002.120.342.350.3116.501.2024.60263.0072.9058.10
    18SS081-824.5023.0047.606.4427.506.222.254.730.894.440.982.160.372.290.3132.601.1825.40255.0069.8060.30
    18SS081-725.8023.8050.506.1728.006.442.324.850.994.540.992.140.372.360.3016.301.2023.90248.0066.1058.20
    邓群村17SS063-251.4038.0086.3010.5047.8010.903.8610.501.9410.301.935.090.814.790.6714.001.4929.50176.0040.2038.90
    19SS023-150.4032.5070.209.4342.509.853.599.421.789.211.774.650.764.430.6426.101.7730.20169.0056.0066.90
    17SS063-145.9035.1076.2010.2045.0011.004.077.811.588.481.743.750.663.890.5814.301.5030.10168.0043.0064.00
    样品编号NiCuZnGaRbSrNbMoCdInCsBaTaWTlPbBiThUZrHf
    紫洞ZD-201-A113.1066.44113.9021.9220.80418.2452.551.860.110.080.49394.903.280.640.134.600.033.721.20251.806.03
    ZD-20195.2861.32106.3020.7919.03412.2648.291.660.110.080.39374.303.040.600.223.560.033.751.14250.706.02
    14SSZD-N108.0051.10107.0022.4024.30447.0071.80−−−0.48402.004.03−−2.41−3.821.02256.005.82
    14SSZD-S2126.0042.80102.0021.6039.80414.0074.20−−−0.36358.004.15−−2.60−4.271.15276.006.18
    ZD-101108.2067.06113.1022.3740.96412.9655.071.600.130.080.37363.003.410.510.163.890.033.631.29250.805.79
    14SSZD-S1132.0042.90106.0022.0020.50388.0073.20−−−0.18370.004.12−−2.91−4.271.14280.006.22
    ZD-20292.6844.73102.3821.6520.92444.7054.922.64−0.080.42398.293.170.530.113.120.0216.471.18263.405.72
    ZD-20386.1040.65105.7021.9134.93410.0054.402.81−0.090.31367.213.130.600.123.290.0110.291.23239.905.35
    18SS082-1101.0064.80156.0022.3020.90451.0048.802.780.250.070.29391.002.84138.000.092.690.023.310.97129.005.40
    18SS082-2108.0062.60161.0022.7024.50437.0050.003.100.220.070.14367.002.91116.000.072.090.013.140.89137.005.68
    下载: 导出CSV
    续表 3
    样品编号NiCuZnGaRbSrNbMoCdInCsBaTaWTlPbBiThUZrHf
    王
    借
    岗    		WJG-20179.5738.12101.9119.0734.53774.8069.841.77−0.080.97693.932.920.440.108.750.037.991.10181.804.13
    WJG-20377.0638.4788.0419.8633.45825.2074.202.90−0.081.01739.073.540.470.103.190.0212.401.23202.704.37
    WJG-20472.9434.1687.3718.7831.68829.5070.361.73−0.070.90742.933.700.450.093.000.0217.891.13192.304.32
    14SS012c93.7036.8083.4019.3046.10895.0090.20−−−0.55572.005.07−−2.51−3.781.09155.003.94
    14SS012b98.3040.3082.3019.7046.90901.0098.00−−−0.56605.005.31−−2.71−3.891.67168.004.34
    14SS01290.7037.2078.4019.0045.80874.0089.70−−−0.50572.004.92−−2.44−3.811.36150.003.94
    18SS083-382.0052.30162.0019.1045.90876.0064.602.830.230.080.71579.003.58116.000.061.780.003.690.83101.004.69
    石
    头
    村    		14SS004-240.7045.60100.0022.6032.40768.0053.10−−−0.32318.003.13−−2.57−2.550.73187.004.56
    18SS079-241.8059.20167.0021.1028.80871.0038.101.990.210.071.90291.002.3096.700.131.380.012.270.6294.504.30
    17SS060-134.7049.00116.0021.4022.50914.0041.101.770.090.071.23699.002.570.350.062.030.012.710.75258.005.55
    17SS060-243.0063.60177.0023.2023.60756.0038.102.070.170.071.01519.002.44201.000.061.600.012.530.6798.804.60
    18SS079-140.0062.20171.0022.0027.20693.0039.902.030.170.071.61329.002.24109.000.104.080.002.410.6398.804.17
    18SS081-841.8066.10172.0022.9025.70610.0038.801.980.210.080.37919.002.24128.000.071.540.002.290.6097.504.48
    18SS081-742.8063.80167.0021.5026.80603.0037.701.900.210.070.58407.002.33107.000.071.480.012.440.6793.304.51
    邓
    群
    村    		17SS063-226.6039.20138.0023.1036.20507.0045.702.100.140.093.07498.002.690.470.323.870.012.830.92322.006.75
    19SS023-131.1035.70143.0022.7023.00447.0041.102.160.270.1211.50411.002.99177.000.366.700.034.170.85266.007.16
    17SS063-135.5050.50200.0025.0030.20490.0044.902.510.220.092.50464.002.56129.000.493.390.012.710.85128.005.56
    下载: 导出CSV

    表  4   三水盆地玄武岩稀土元素特征值

    Table  4   The characteristic values for REE of basalts in Sanshui Basin

    样品编号 ΣREE ΣLREE ΣHREE LREE/HREE δEu δCe (La/Sm)N (La/Yb)N (Sm/Nd)N (Gd/Yb)N
    ZD-201-A148.66127.5321.136.0354951.0137540.8595642.8495357.0626460.7248781.748391
    ZD-201142.64122.5720.076.1071251.0245220.8605972.8132267.0753810.7390981.73334
    14SSZD-N143.71123.8319.886.2288731.0154990.9126722.8372087.8284460.6996072.077396
    14SSZD-S2143.2123.9619.246.4428271.0044460.9127953.0038657.8976960.6925961.946561
    ZD-101143.59123.6219.976.1902851.0228250.8602572.9328387.4641050.7211721.807326
    14SSZD-S1145.55126.0119.546.4488230.9972110.9065553.0164728.0958490.6911762.014067
    ZD-202150.53130.1920.346.4006881.0870790.8917463.2129987.4470990.6840091.656671
    ZD-203156.92135.7721.156.4193851.0821720.8821313.2939327.5910050.6756581.677244
    18SS082-1140.33121.9318.46.626631.1699540.8863262.7153886.3967850.7285071.377047
    18SS082-2142.37124.0318.346.7628141.1349390.9119612.8647986.6707180.6954071.500245
    WJG-201194.11172.9721.148.1821191.0852930.9003013.60467310.617280.6251342.03861
    WJG-203198.15176.7221.438.2463841.0801070.9029013.62857610.545960.622931.993814
    WJG-204182.29162.5319.768.2252021.0754670.887653.71836210.690780.6173942.014172
    14SS012c180.42161.0519.378.3144041.048340.9124943.35008712.246220.6260252.607968
    14SS012b189.39168.8220.578.2070981.0379520.9169233.28336411.59140.6348272.477477
    14SS012178.2159.2818.928.4186051.0573210.91433.33333312.477610.6315792.625598
    18SS083-3180.26161.7918.478.759611.106330.8944493.2100559.9256270.6265891.844457
    14SS004-2133.74115.9317.816.5092641.0734240.9424342.4388978.1720430.6951152.429
    18SS079-2125.77109.7915.986.8704631.2190460.9143752.3859846.8008170.7234931.631519
    17SS060-1143.88123.0820.85.9173081.0422430.9120522.547246.494070.6725571.80377
    17SS060-2134.69117.7616.936.95571.2759290.963822.3638167.258210.6794871.828359
    18SS079-1129.47113.0616.416.8897011.1484690.970322.2621356.5411120.7145671.59673
    18SS081-8129.18113.0116.176.9888681.2205840.9278742.3260046.7713760.6959441.666757
    18SS081-7133.77117.2316.547.0876661.2198530.9815632.3246846.7990710.7076921.65835
    17SS063-2233.39197.3636.035.4776581.0891351.0238092.1929565.3485080.7016411.768888
    19SS023-1200.73168.0732.665.146051.1243460.9540372.0754874.9461150.7131221.715907
    17SS063-1210.06181.5728.496.3731131.279370.9582462.0071856.083340.7521371.620123
    下载: 导出CSV

    表  5   三水盆地岩浆熔融柱计算结果

    Table  5   Calculation results of mantle melting column in Sanshui Basin

    FeO8.0/%Na2O8.0/%FeOprim/%Na2Oprim/%P0 /GPaPf /GPaT0 /℃Tf /℃F/ %Z0 /kmZf /km
    ZD10.033.4210.092.923.032.251 5171 4719.22101.287376.6719
    WJG10.723.2510.682.723.322.481 5551 50610.06110.392584.01884
    　　注：FeO8.0、Na2O8.0分别表示样品中的FeO、Na2O在MgO=8%时的含量；FeOprim、Na2Oprim分别表示样品所代表的母岩浆在未发生橄榄石分离结晶时FeO、Na2O的含量；P0、T0、Z0分别表示部分熔融作用发生时的初始压力、温度及深度；Pf、Tf、Zf分别表示部分熔融作用结束时的终止压力、温度及深度；F表示压力从P0下降至Pf的过程中的总部分熔融程度。
    下载: 导出CSV

    表  6   地幔源区矿物组成及La、Yb、Sm总分配系数

    Table  6   Mineral composition in mantle source and total partition coefficients of La, Yb and Sm

    岩石类型矿物相/% 总分配系数
    橄榄石单斜辉石斜方辉石石榴石尖晶石LaYbSm
    石榴石橄榄岩0.540.090.170.200.005112780.868120.129968
    尖晶石橄榄岩0.460.180.2800.080.009839220.116340.05835
    三水盆地源区0.510.140.240.060.050.007717570.3361550.075927
    下载: 导出CSV
    • 参考文献(53)
  • [1]

    Tapponnier P, Peltzer G, Le Dain A Y, et al. Propagating extrusion tectonics in Asia: new insights from simple experiments with plasticine [J]. Geology, 1982, 10(12): 611-616. doi: 10.1130/0091-7613(1982)10<611:PETIAN>2.0.CO;2
    [2]

    Mai H A, Chan Y L, Yeh M W, et al. Tectonic implications of mesozoic magmatism to initiation of cenozoic basin development within the passive south China Sea margin [J]. International Journal of Earth Sciences, 2018, 107(3): 1153-1174. doi: 10.1007/s00531-017-1537-y
    [3]

    Sun W D. Initiation and evolution of the South China Sea: an overview [J]. Acta Geochimica, 2016, 35(3): 215-225. doi: 10.1007/s11631-016-0110-x
    [4]

    Holloway N H. North palawan block, philippines - its relation to asian mainland and role in evolution of South China Sea [J]. AAPG Bulletin, 1982, 66(9): 1355-1383.
    [5]

    Wu J, Suppe J. Proto-South China Sea plate tectonics using subducted slab constraints from tomography [J]. Journal of Earth Science, 2018, 29(6): 1304-1318. doi: 10.1007/s12583-017-0813-x
    [6]

    Xu Y G, Wei J X, Qiu H N, et al. Opening and evolution of the South China Sea constrained by studies on volcanic rocks: preliminary results and a research design [J]. Chinese Science Bulletin, 2012, 57(24): 3150-3164. doi: 10.1007/s11434-011-4921-1
    [7]

    Chen L, Hu J W, Yang D H, et al. Kinematic models for the opening of the South China Sea: an upwelling divergent flow origin [J]. Journal of Geodynamics, 2017, 107: 20-33. doi: 10.1016/j.jog.2017.03.002
    [8] 林间, 李家彪, 徐义刚, 等. 南海大洋钻探及海洋地质与地球物理前沿研究新突破[J]. 海洋学报, 2019, 41(10):125-140

    LIN Jian, LI Jiabiao, XU Yigang, et al. Ocean drilling and major advances in marine geological and geophysical research of the South China Sea [J]. Acta Oceanologica Sinica, 2019, 41(10): 125-140.
    [9]

    Sun Z, Lin J, Qiu N, et al. The role of magmatism in the thinning and breakup of the south China sea continental margin: special topic: the South China Sea ocean drilling [J]. National Science Review, 2019, 6(5): 871-876. doi: 10.1093/nsr/nwz116
    [10]

    Ebinger C J, Yemane T, Woldegabriel G, et al. Late eocene - recent volcanism and faulting in the Southern Main Ethiopian rift [J]. Journal of the Geological Society, 1993, 150(1): 99-108. doi: 10.1144/gsjgs.150.1.0099
    [11]

    Christiansen R L, Foulger G R, Evans J R. Upper-mantle origin of the yellowstone hotspot [J]. GSA Bulletin, 2002, 114(10): 1245-1256. doi: 10.1130/0016-7606(2002)114<1245:UMOOTY>2.0.CO;2
    [12]

    Putirka K D, Perfit M, Ryerson F J, et al. Ambient and excess mantle temperatures, olivine thermometry, and active vs. passive upwelling [J]. Chemical Geology, 2007, 241(3-4): 177-206. doi: 10.1016/j.chemgeo.2007.01.014
    [13] 阎贫, 刘海龄. 南海及其周缘中新生代火山活动时空特征与南海的形成模式[J]. 热带海洋学报, 2005, 24(2):33-41 doi: 10.3969/j.issn.1009-5470.2005.02.005

    YAN Pin, LIU Hailing. Temporal and spatial distributions of meso-enozoic igneous rocks over south China Sea [J]. Journal of Tropical Oceanography, 2005, 24(2): 33-41. doi: 10.3969/j.issn.1009-5470.2005.02.005
    [14]

    Chung S L, Cheng H, Jahn B M, et al. Major and trace element, and Sr-Nd Isotope constraints on the origin of paleogene volcanism in South China Prior to the South China Sea opening [J]. Lithos, 1997, 40(2-4): 203-220. doi: 10.1016/S0024-4937(97)00028-5
    [15] 董月霞, 肖龙, 周海民, 等. 广东三水盆地双峰式火山岩: 空间展布、岩石学特征及其盆地动力学意义[J]. 大地构造与成矿学, 2006, 30(1):82-92 doi: 10.3969/j.issn.1001-1552.2006.01.010

    DONG Yuexia, XIAO Long, ZHOU Haimin, et al. Spatial distribution and petrological characteristics of the bimodal volcanic rocks from Sanshui Basin, Guangdong Province: implication for basin dynamics [J]. Geotectonica et Metallogenia, 2006, 30(1): 82-92. doi: 10.3969/j.issn.1001-1552.2006.01.010
    [16] 肖龙, 周海民, 董月霞, 等. 广东三水盆地火山岩: 地球化学特征及成因-兼论火山岩性质的时空演化和南海形成的深部过程[J]. 大地构造与成矿学, 2006, 30(1):72-81 doi: 10.3969/j.issn.1001-1552.2006.01.009

    XIAO Long, ZHOU Haimin, DONG Yuexia, et al. Geochemistry and petrogenesis of cenozoic volcanic rocks from Sanshui basin: implications for spatial and temporal variation of rock types and constraints on the formation of South China Sea [J]. Geotectonica et Metallogenia, 2006, 30(1): 72-81. doi: 10.3969/j.issn.1001-1552.2006.01.009
    [17] 袁晓博. 三水盆地新生代岩浆记录与南海早期演化[D]. 中国地质大学(北京)博士学位论文, 2019.

    YUAN Xiaobo. The record of cenozoic magmatism in Sanshui basin and its relationship with the early tectonic evolution stage of the South China Sea[D]. Doctor Dissertation of Chian University of Geoscience (Beijing), 2019.
    [18] 袁晓博, 方念乔. 三水盆地中渐新世火山记录的新建与南海扩张[J]. 地质通报, 2019, 38(4):689-695

    YUAN Xiaobo, FANG Nianqiao. The new volcanics record in sanshui basin and its relationship with the spreading of the South China Sea [J]. Geological Bulletin of China, 2019, 38(4): 689-695.
    [19] 张维, 方念乔. 广东三水盆地始新世火山岩地球化学特征[J]. 地球科学—中国地质大学学报, 2014, 39(1):37-44 doi: 10.3799/dqkx.2014.004

    ZHANG Wei, FANG Nianqiao. Geochemistry characteristics of eocene volcanic rocks in Sanshui basin, Guangdong [J]. Earth Science—Journal of China University of Geosciences, 2014, 39(1): 37-44. doi: 10.3799/dqkx.2014.004
    [20] 杨蜀颖. 南海玳瑁海山与相邻陆域玄武岩的地球化学特征及其构造意义[D]. 中国地质大学(北京)博士学位论文, 2015.

    YANG Shuying. Geochemical characteristics of basalts from the daimao seamount in the South China Sea (SCS) and from the SCS's neighboring lands: implications for the regional tectonic evolution[D]. Doctor Dissertation of Chian University of Geoscience (Beijing), 2015.
    [21] 朱炳泉, 王慧芬, 陈毓蔚, 等. 新生代华夏岩石圈减薄与东亚边缘海盆构造演化的年代学与地球化学制约[J]. 地球化学, 2002, 31(3):213-221 doi: 10.3321/j.issn:0379-1726.2002.03.001

    ZHU Bingquan, WANG Huigfen, CHEN Yuwei, et al. Geochronological and geochemical constraint on the cenozoic extension of cathaysian lithosphere and tectonic evolution of the Border Sea Basins in East Asia [J]. Geochimica, 2002, 31(3): 213-221. doi: 10.3321/j.issn:0379-1726.2002.03.001
    [22] 唐忠驭. 广东三水盆地白垩纪—早第三纪裂谷型火山作用[J]. 广东地质, 1994, 9(1):49-57

    TANG Zhongyu. Cretaceous-eogene rift valley-type volcanism in Sanshui Basin, Guangdong [J]. Guangdong Geology, 1994, 9(1): 49-57.
    [23]

    Ackerman L, Ulrych J, Řanda Z, et al. Geochemical characteristics and petrogenesis of phonolites and trachytic rocks from the ceske stredohori volcanic complex, the Ohre Rift, Bohemian Massif [J]. Lithos, 2015, 224-225: 256-271. doi: 10.1016/j.lithos.2015.03.014
    [24]

    Yoder H S. Citation-Classic Origin of basalt magmas - an experimental-study of natural and synthetic rock systems [J]. Current Contents, 1986(39): 14-20.
    [25]

    Shand S J. The problem of the alkaline rocks [J]. Proceedings of the Geological Society of South Africa, 1922, 25: 19-33.
    [26]

    Middlemost E A K. Naming materials in the magma/igneous rock system [J]. Earth-Science Reviews, 1994, 37(3-4): 215-24. doi: 10.1016/0012-8252(94)90029-9
    [27] 焦守涛, 张旗, 葛粲, 等. 碱性岩及碱性与亚碱性岩系列的界线: 基于全球火山岩数据的探讨[J]. 地质通报, 2019, 38(12):1955-1962

    JIAO Shoutao, ZHANG Qi, GE Can, et al. Alkaline rock and the distinction between alkaline and sub-alkaline: a discussion on data of global volcanic rocks [J]. Geological Bulletin of China, 2019, 38(12): 1955-1962.
    [28]

    McDonough W F, Sun S S, Ringwood A E, et al. Potassium, rubidium, and cesium in the earth and moon and the evolution of the mantle of the earth [J]. Geochimica et Cosmochimica Acta, 1992, 56(3): 1001-1112. doi: 10.1016/0016-7037(92)90043-I
    [29]

    Gaschnig R M, Rudnick R L, McDonough W F, et al. Compositional evolution of the upper continental crust through time, as constrained by ancient glacial diamictites [J]. Geochimica et Cosmochimica Acta, 2016, 186: 316-43. doi: 10.1016/j.gca.2016.03.020
    [30]

    Yu X, Liu Z F. Non-mantle-plume process caused the initial spreading of the South China Sea [J]. Scientific Reports, 2020, 10(1): 10. doi: 10.1038/s41598-019-56089-4
    [31]

    Frey F A, Green D H, Roy S D. Integrated models of basalt petrogenesis: a study of quartz tholeiites to olivine melilitites from south Eastern Australia utilizing geochemical and experimental petrological data [J]. Journal of Petrology, 1978, 19(3): 463-513. doi: 10.1093/petrology/19.3.463
    [32]

    Langmuir C H, Klein E M, Plank T. Petrological systematics of mid-ocean ridge basalts: constraints on melt generation beneath ocean ridges[M]//Morgan J P, Blackman D K, Sinton J M. Mantle Flow and Melt Generation at Mid-Ocean Ridges. Washington, D. C. : Geophysical Monograph Series, 1992: 183-280.
    [33] 梁涛, 罗照华, 李德东, 等. 托云盆地新生代幔源岩浆源区起止深度的限定[J]. 岩石学报, 2008, 24(12):2820-2838

    LIANG Tao, LUO Zhaohua, LI Dedong, et al. Source location identification of cenozoic mantle-derived magma in Tuyon Basin [J]. Acta Petrologica Sinica, 2008, 24(12): 2820-2838.
    [34]

    Wang K, Plank T, Walker J D, et al. A mantle melting profile across the basin and range, SW USA [J]. Journal of Geophysical Research: Solid Earth, 2002, 107(B1): ECV 5-1-ECV 5-21. doi: 10.1029/2001JB000209
    [35]

    Cheng L L, Liang T, Zeng L, et al. Mantle melting column software named Calmantle 1.0 and the preliminary discussion on the thickness variation of the Tengchong Cenozoic lithosphere [J]. Earth Science Frontiers, 2012, 19(4): 126-134.
    [36]

    Jones C H, Wernicke B P, Farmer G L, et al. Variations across and along a major continental rift: an interdisciplinary study of the basin and Range province, Western USA [J]. Tectonophysics, 1992, 213(1-2): 57-96. doi: 10.1016/0040-1951(92)90252-2
    [37] 黄海波, 郭兴伟, 夏少红, 等. 华南沿海地区地壳厚度与泊松比研究[J]. 地球物理学报, 2014, 57(12):3896-3906 doi: 10.6038/cjg20141204

    HUANG Haibo, GUO Xingwei, XIA Shaohong, et al. Crustal thickness and poisson's ratio in the coastal areas of South China [J]. Chinese Journal of Geophysics, 2014, 57(12): 3896-3906. doi: 10.6038/cjg20141204
    [38]

    Kelemen P B, Shimizu N, Dunn T. Relative depletion of niobium in some arc magmas and the continental crust: partitioning of K, Nb, La And Ce during melt/rock reaction in the upper mantle [J]. Earth and Planetary Science Letters, 1993, 120(3-4): 111-134. doi: 10.1016/0012-821X(93)90234-Z
    [39] 徐义刚, 钟孙霖. 峨眉山大火成岩省: 地幔柱活动的证据及其熔融条件[J]. 地球化学杂志, 2001, 30(1):1-9

    XU Yigang, ZHONG Sunlin. The Emeishan large Igneous province: evidence for mantle plume activity and melting conditions [J]. Geochimica, 2001, 30(1): 1-9.
    [40] 陈盼盼. 三水盆地晚白垩世-始新世火山-沉积序列对南海北缘构造演化的响应[D]. 博士学位论文中国地质大学(北京), 2018.

    CHEN Panpan. The response of late cretaceous-eocene epoch volcanic and sedimentary sequence in sanshui basin to the tectonic evolution of the Northern margin of Southe China Sea[J]. Doctor Dissertation of Chian University of Geoscience (Beijing), 2018.
    [41] 邹和平, 李平鲁, 饶春涛. 珠江口盆地新生代火山岩地球化学特征及其动力学意义[J]. 地球化学, 1995, 24(S1):33-45

    ZOU Heping, LI Pinglu, RAO Chuntao. Geochemistry of cenozoic volcanic rocks in zhu jiangkou basin and its geodynamic significance [J]. Geochimica, 1995, 24(S1): 33-45.
    [42] 张斌, 王璞珺, 张功成, 等. 珠—琼盆地新生界火山岩特征及其油气地质意义[J]. 石油勘探与开发, 2013, 40(6):657-665 doi: 10.11698/PED.2013.06.03

    ZHANG Bin, WANG Pujun, ZHANG Gongcheng, et al. Cenozoic volcanic rocks in the pearl river mouth and southeast Hainan Basins of South China Sea and their implications for petroleum geology [J]. Petroleum Exploration and Development, 2013, 40(6): 657-665. doi: 10.11698/PED.2013.06.03
    [43] 李思伟. 珠江口盆地惠州凹陷新生代火山岩: 从岩石成因到火山岩储层[D]. 吉林大学博士学位论文, 2020.

    LI Siwei. Cenozoic volcanic rocks in Huizhou Sag of the Pearl River mouth basin: from petrogenesis to volcanic reservoir[D]. Doctor Dissertation of Jilin University, 2020.
    [44] 汪云亮, 张成江, 修淑芝. 玄武岩类形成的大地构造环境的Th/Hf-Ta/Hf图解判别[J]. 岩石学报, 2001, 17(3):413-421 doi: 10.3969/j.issn.1000-0569.2001.03.009

    WANG Yunliang, ZHANG Chengjiang, XIU Shuzhi. Th/Hf-Ta/Hf identification of tectonic setting of basalts [J]. Acta Petrologica Sinica, 2001, 17(3): 413-421. doi: 10.3969/j.issn.1000-0569.2001.03.009
    [45]

    Fitton J G, Saunders A D, Norry M J, et al. Thermal and Chemical Structure of the Iceland Plume [J]. Earth and Planetary Science Letters, 1997, 153(3-4): 197-208. doi: 10.1016/S0012-821X(97)00170-2
    [46]

    Bevins R E, Kokelaar B P, Dunkley P N. Petrology and geochemistry of lower to middle ordovician igneous rocks in wales: a volcanic arc to marginal basin transition [J]. Proceedings of the Geologists’ Association, 1984, 95(4): 337-347. doi: 10.1016/S0016-7878(84)80064-4
    [47]

    Weaver B L. The origin of ocean island basalt end-member compositions: trace element and isotopic constraints [J]. Earth and Planetary Science Letters, 1991, 104(2-4): 381-391. doi: 10.1016/0012-821X(91)90217-6
    [48]

    Hart S R, Hauri E H, Oschmann L A, et al. Mantle plumes and entrainment: isotopic evidence [J]. Science, 1992, 256(5056): 517-520. doi: 10.1126/science.256.5056.517
    [49] 阙晓铭, 李元森, 陈会霞, 等. 深部地幔在白云凹陷的岩浆记录: 基于BY7火山岩的地球化学研究[J]. 华南地质与矿产, 2013, 29(2):105-115

    QUE Xiaoming, LI Yuansen, CHEN Huixia, et al. Geochemistry research on the deep mantle activity in baiyun sag during the cenozoic from the volcanic rocks of well By7 [J]. Geology and Mineral Resources of South China, 2013, 29(2): 105-115.
    [50]

    Zhang G L, Sun W D, Seward G. Mantle source and magmatic evolution of the dying spreading ridge in the South China Sea [J]. Geochemistry, Geophysics, Geosystems, 2018, 19(11): 4385-4399. doi: 10.1029/2018GC007570
    [51]

    Zou H B, Fan Q C. U-Th isotopes in Hainan basalts: Implications for sub-asthenospheric origin of EM2 mantle endmember and the dynamics of melting beneath Hainan Island [J]. Lithos, 2010, 116(1-2): 145-152. doi: 10.1016/j.lithos.2010.01.010
    [52]

    Lee C T A, Luffi P, Plank T, et al. Constraints on the depths and temperatures of basaltic magma generation on Earth and other terrestrial planets using new thermobarometers for mafic magmas [J]. Earth and Planetary Science Letters, 2009, 279(1-2): 20-33. doi: 10.1016/j.jpgl.2008.12.020
    [53]

    Katz R F, Spiegelman M, Langmuir C H. A new parameterization of hydrous mantle melting [J]. Geochemistry, Geophysics, Geosystems, 2003, 4(9): 1073.
    • 相关文章
    • 施引文献
    • 资源附件(0)
图(11)  /  表(7)
计量
  • 文章访问数:  1629
  • HTML全文浏览量:  335
  • PDF下载量:  39
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-09-28
  • 修回日期:  2020-12-22
  • 网络出版日期:  2021-02-21
  • 刊出日期:  2021-06-27

目录

摘要
HTML全文

  • 1.   区域概况

  • 2   试验样品来源及方法

  • 2.1   试验材料

  • 2.2   试验样品制备与试验设置

  • 3.   结果和分析

  • 3.1   试验结果

  • 3.2   Singh-Mitchell蠕变模型

  • 3.3   试验用粉黏土蠕变特性分析

  • 3.4   含水合物试样蠕变特性分析

  • 4.   讨论

  • 4.1   蠕变试验结果分析

  • 4.2   基于Singh-Mitchell的修正模型

  • 5.   结论

参考文献(53)
相关文章
施引文献
资源附件(0)

/

返回文章
返回
微信公众号 微信公众号
订购码 订购码

版权所有 © 《海洋地质与第四纪地质》编辑部    备案号：鲁ICP备15036216号-7

地址：山东省青岛市即墨区观山路596号    邮编：266237    联系电话：0532-85755823    E-mail：hydzdsjdz@163.com

本系统由北京仁和汇智信息技术有限公司开发 技术支持： info@rhhz.net 百度统计