北黄海盆地褶皱分布与成因类型

王后金, 梁世荣, 陈强

王后金, 梁世荣, 陈强. 北黄海盆地褶皱分布与成因类型[J]. 海洋地质与第四纪地质, 2006, 26(6): 65-70.
引用本文: 王后金, 梁世荣, 陈强. 北黄海盆地褶皱分布与成因类型[J]. 海洋地质与第四纪地质, 2006, 26(6): 65-70.
WANG Hou-jin, LIANG Shi-rong, CHEN Qiang. DISTRIBUTION AND ORIGIN TYPES OF THE FOLDS IN THE NORTH YELLOW SEA BASIN[J]. Marine Geology & Quaternary Geology, 2006, 26(6): 65-70.
Citation: WANG Hou-jin, LIANG Shi-rong, CHEN Qiang. DISTRIBUTION AND ORIGIN TYPES OF THE FOLDS IN THE NORTH YELLOW SEA BASIN[J]. Marine Geology & Quaternary Geology, 2006, 26(6): 65-70.

北黄海盆地褶皱分布与成因类型

详细信息
    作者简介:

    王后金(1972-),男,高级工程师,硕士生,从事海洋地质和油气地质研究

  • 中图分类号: P738

DISTRIBUTION AND ORIGIN TYPES OF THE FOLDS IN THE NORTH YELLOW SEA BASIN

  • 摘要: 北黄海盆地为叠置在华北地台上的中新生代断陷盆地,在地史上经受了多期区域构造运动的改造,地层变形强烈,褶皱构造相当发育,分为伸展型、压缩型、重力型和热力型4大成因类型。其中,与张性正断层相关的伸展型褶皱分布范围最广,可进一步划分为牵引褶皱、逆牵引褶皱、滚动褶皱、翘倾褶皱、补偿性挤压背斜和断层-位移褶皱;形成于区域挤压应力场的压缩型褶皱包括非反转挤压褶皱和反转褶皱,在东部和中部坳陷表现明显;重力型褶皱以发育于古潜山之上的披覆背斜比较常见;由岩浆活动引起的热力型褶皱不甚发育,目前仅见于东部坳陷。
    Abstract: North Yellow Sea basin is a Mesozoic/Cenozoic fault basin superposing on North China platform, which has distorted stratum and well-developed fold structures due to the reconstruction of several structure movements. The folds in the North Yellow Sea basin are mainly classified into four origin types, namely extensional folds, compressional folds, gravitational force folds and thermogene folds. The extensional folds related to normal faults are distributed most widely, which can be further divided into drag folds, reverse-drag folds, rolled folds, compensation squeezed anticlines and fault-displacement folds;the compressional folds including non-reverse squeezed folds and reverse folds are caused by territorial extrusion stress and spread both in the eastern sag and the mid-sag apparently. As a kind of gravitational force folds, buried-hill drape anticlines are quite easily observed in the North Yellow Sea basin, while the thermogene folds only appear in the eastern sag owing to inactive magmation.
  • 天然气水合物是一种由水分子和气体分子在低温高压条件下形成的似冰状笼形结晶化合物,广泛赋存于海底沉积物和高原冻土带中,被公认为资源量巨大的潜在能源[1],得到了科技界和各国政府的关注。我国已于2017年5月在南海神狐海域成功进行了天然气水合物的试采[2]

    含天然气水合物沉积介质是赋存有天然气水合物的多孔介质。在自然条件下,外界温度和压力条件的改变会导致地层中的水合物分解。水合物的分解一方面降低了颗粒之间的胶结性,另一方面分解产生的液体会使地层液化,降低含水合物地层的稳定性,进而可能诱发海底地质灾害[3]。在水合物开采过程中,降压、热激等方法的应用以及钻井液的入侵都会改变含水合物沉积介质的力学性质,弱化沉积介质的强度,影响井壁稳定性和储层稳定性,可能诱发生产事故和地质灾害[4]。对含水合物沉积介质的力学性质进行研究,建立能有效描述含水合物沉积介质力学特性的力学模型,是研究含水合物地层井壁稳定性和储层稳定性的前提。

    前期系列试验研究初步揭示了含水合物沉积介质的力学特性[5-16]。研究表明水合物饱和度[7, 9-11, 15]、围压[5, 6]、赋存模式[17]、剪切速率[8]、温度[8, 15]等都会影响含水合物沉积介质的力学性质。因此,在构建力学模型的过程中应该尽可能全面地考虑这些影响因素。在大量试验的基础上,一些学者建立了不同的含水合物沉积介质的力学模型。Klar等[18]将含水合物沉积介质的刚度、粘聚力和剪胀角看作水合物饱和度的函数,利用摩尔—库伦破坏准则,建立了含水合物沉积介质的弹塑性力学模型。该模型参数较少,能较好地描述含水合物沉积介质峰值强度及弹性模量与水合物饱和度之间的关系,但不能反映出应变软化特性。Uchida等[19]在传统临界状态模型的基础上通过修正屈服函数,建立了适用于降压开采过程的含水合物沉积介质的力学模型。该模型虽然可以描述沉积介质应力-应变的全过程,但参数多且部分参数物理意义不明确。吴二林等[20, 21]将含水合物沉积介质假定为线弹性材料,采用几何损伤理论建立了含水合物沉积介质的本构模型。但损伤模型假设材料破坏后完全失去承载力,这与实际情况不符[22]。在众多的岩土力学模型中,邓肯-张模型(Duncan-Chang)因其可以准确地描述含水合物沉积介质应力-应变的非线性弹性关系,且参数物理意义明确,得到了广泛应用。Miyazaki等[23]将水合物饱和度和有效围压引入经典的邓肯-张模型,建立了适用于含天然气水合物的非线性弹性模型。

    虽然有学者[10, 23, 24]借鉴Miyzaki[17]的研究对含水合物沉积介质抗剪强度与水合物饱和度和有效围压之间的数学关系进行了探索,但尚未建立细砂质含水合物沉积介质的非线性弹性模型。笔者以细砂质含水合物沉积介质的三轴力学试验实测数据为基础,结合Miyazaki等[23]提出的含水合物沉积介质邓肯-张模型,建立了适合南海含水合物沉积介质的非线性弹性模型,并对模型的准确性进行了验证。

    本次试验采用的仪器为中国科学院广州能源研究所自建的含天然气水合物三轴力学试验装置。装置的技术指标为:反应釜为圆柱形,尺寸为φ50mm×100mm;三轴腔室的最大围压为30MPa,最大孔隙压力为30MPa,最大反压为20MPa,最大轴向压力为250KN,压力传感器的精度误差优于0.12%F.S。温度的控制范围为—30~50℃。

    本文试验用沉积细砂由广州海洋地质调查局提供,由“海洋四号”地质调查船在南海北部海底钻取。沉积细砂经过标准检验筛(60~100目)筛选[25]

    试验步骤主要包括:(1)检漏;(2)制备沉积试样;(3)原位法生成水合物;(4)三轴剪切;(5)收集水合物分解气反算饱和度。试验装置和试验步骤的具体细节详见前期研究成果[25]

    通过一系列试验,分别选取了单调加载和循环加载的试验结果进行研究,具体的试验结果如下:

    图 1给出了本次试验的应力摩尔圆,并通过应力摩尔圆的公切线得到了强度包络线,强度包络线在应力轴的截距为粘聚力c,强度包络线与应变轴的交角为内摩擦角φ

    图  1  含水合物沉积介质摩尔圆
    Figure  1.  The mohr circles of hydrate-bearing sediments

    图 2给出了不含水合物沉积试样和水合物饱和度为35%时沉积试样的应力-应变关系。从图中可以看出,在有效围压不变的条件下,随着饱和度的增大,含水合物沉积介质的力学强度和刚度也明显增大。这与前人的试验结论[5, 26]一致。

    图  2  含水合物沉积介质的应力-应变曲线计算值与试验值对比
    Figure  2.  Comparison of calculated results(solid lines) and experimental results(broken lines) for hydrate-bearing sediments

    图 3给出了循环加载试验的结果,循环加载主要是分2次降低了沉积介质中水合物饱和度。从图中可知,在水合物饱和度和有效围压同时改变的条件下,沉积介质的应力-应变曲线依然表现出非线性弹性关系。

    图  3  循环加载条件下含水合物沉积介质应力-应变滞回环
    (虚线为完整滞回环,实线为循环加载弹塑性阶段的模型计算值)
    Figure  3.  Stress-strain hysteresis loop of hydrate-bearing sediments
    (the dotted line is a complete hysteresis loop, the solid line is the calculated value of the model)

    Duncan等[27]基于可以用双曲线拟合应力-应变关系的假设,提出了一种目前被广泛应用的增量弹性模型,称为邓肯-张模型。其应力-应变关系如下:

    $$ \sigma_{1}-\sigma_{3}=\frac{\varepsilon_{a}}{a+b \varepsilon_{a}} $$ (1)

    其中,σ1为轴向应力(MPa);σ3为围压(MPa);σ1-σ3为偏差应力(MPa);εa为轴向应变(%);ab为试验常数。通过引入初始变形模量Ei、强度(σ1-σ3)f和破坏比Rf,式(1)可以改写为(2):

    $$ \sigma_{1}-\sigma_{3}^{\prime}=\frac{\varepsilon_{a}}{\frac{1}{E_{i}}+\frac{R_{f} \cdot \varepsilon_{a}}{\left(\sigma_{1}-\sigma_{3}^{\prime}\right)_{f}}} $$ (2)

    需要说明的是:在含水合物沉积力学研究中,认为含水合物沉积介质的力学强度是轴向应力σ1与有效围压σ3的差值,因此在公式(2)中将(1)中的σ3改为σ3

    由2.1的公式推导可知,构建含水合物沉积介质的非线性弹性模型,需要建立沉积介质的破坏强度(σ1-σ3)f、初始弹性模量Ei和破坏比Rf的数学关系。下面给出这3个参数的数学表达。

    含水合物沉积介质的力学强度主要受水合物饱和度和有效围压的影响。沉积介质的粘聚力可视作饱和度的函数,内摩擦角可视为常量[22]。以摩尔-库伦强度准则为基础,考虑水合物饱和度的影响,含水合物沉积介质的力学强度(σ1-σ3)f可以表示为[23]

    $$\left(\sigma_{1}-\sigma_{3}^{\prime}\right)_{f}=\frac{2 \cdot \cos \varphi}{1-\sin \varphi} c+\alpha \cdot S_{h}^{\beta}+\frac{2 \cdot \sin \varphi}{1-\sin \varphi} \sigma_{3}^{\prime} $$ (3)

    其中,φ为内摩擦角(°);cSh=0条件下沉积物的粘聚力;σ3为有效围压;αβ为试验参数。

    根据摩尔-库伦强度理论,通过摩尔圆的公切线可以得到内摩擦角和粘聚力的值(图 1)。利用不同饱和度下沉积介质的强度值(图 2图 3),通过最小二乘法可以回归得到αβ的值。公式(3)中的参数值见表 1

    表  1  细砂质含水合物沉积介质强度的基本参数
    Table  1.  Basic parameter for strength of fine sandy hydrate-bearing sediments
    试样 α(-) β(-) φ/(°) c0/MPa
    含水合物沉积介质 0.003 2.1 9.54 1.414
    下载: 导出CSV 
    | 显示表格

    Miyazaki等[23]研究认为,含水合物沉积介质的初始弹性模量受饱和度和有效围压的影响,与有效围压之间呈幂关系变化。初始弹性模量:

    $$E_{i}=\left(1+\gamma \cdot S_{h}^{\delta}\right) \cdot e_{i 0} \cdot\left(\sigma_{3}^{\prime}\right)^{n} $$ (4)

    其中,ei0为饱和度为零时(Sh=0),沉积介质的初始弹性模量;γδn为试验参数。公式(4)的等号两边物理量纲有差异,在未来的研究中可以进一步改进。

    通过含水合物沉积介质应力-应变关系曲线(图 2图 3),在图 2中可以获得2组初始切线模量,在图 3中可以获得3组不同饱和度下的初始切线模量,利用饱和度为零的初始切线模量ei0和另外4组切线模量,利用最小二乘法回归得到的公式(4)中参数的数值(表 2)。

    表  2  细砂质含水合物沉积介质初始弹性模量公式的基本参数
    Table  2.  Basic parameters for initial tangent elastic modulus of fine sandy hydrate-bearing sediments
    试样 γ(-) δ(-) ei0/MPa n(-)
    含水合物沉积介质 0.0015 2.01 395 0.7
    下载: 导出CSV 
    | 显示表格

    Duncan等[20]在总结大量的土力学试验数据的基础上,给出了破坏比Rf的经验公式,一般取值范围为0.7~0.9。Miyazaki等[17]通过对利用浦丰砂以及7号、8号硅砂合成的水合物沉积介质的三轴试验,总结发现破坏比Rf的取值范围同样为0.7~0.9,但尚未给出破坏比Rf与有效围压和饱和度的关系。本文以Duncan的经验公式为基础,推导了含水合物沉积介质破坏比的计算公式,具体方法如下:

    结合公式(2)和公式(3),可以得到破坏比Rf的经验表达式:

    $$ R_{f}=\frac{\frac{2 \cdot \cos \varphi}{1-\sin \varphi}\left(c_{0}+\sigma_{3}^{\prime}\right)+\alpha \cdot S_{h}^{\beta}}{\left(\sigma_{1}-\sigma_{3}\right)_{u l t}} $$ (5)

    Duncan等[26]给出了求取极限偏差应力(σ1-σ3)ult的经验公式:

    $$\left(\sigma_{1}-\sigma_{3}\right)_{u l t}=\frac{\left(\varepsilon_{a}\right)_{95 \%}-\left(\varepsilon_{a}\right)_{70 \%}}{\left(\frac{\varepsilon_{a}}{\sigma_{1}-\sigma_{3}}\right)_{95 \%}-\left(\frac{\varepsilon_{a}}{\sigma_{1}-\sigma_{3}}\right)_{70 \%}} $$ (6)

    其中,$\left(\varepsilon_{a}\right) _{70 \% / 95 \%}$和$\left(\frac{\varepsilon_{a}}{\sigma_{1}-\sigma_{3}}\right)_{70 \% / 95 \%}$为(σ1-σ3)达到70%/95%(σ1-σ3)f时的试验值。该公式是总结常规土力学试验得到的,考虑到含水合物沉积介质与常规土的不同,公式的准确度还需要大量含水合物沉积介质的力学试验数据的检验和校正。

    因此,由公式(5)和公式(6)可以计算含水合物沉积介质的破坏比Rf。本研究只是通过理论推导得到破坏比Rf的表达式,该公式的适用范围还有待更多的试验数据进一步的验证。

    由前文的推导可知,含水合物沉积介质的应力-应变关系可以表示为:

    $$ \sigma_{1}-\sigma_{3}=\\ \frac{{{\varepsilon _a}}}{{\frac{1}{{\left[ {1 + \gamma \cdot S_h^\delta } \right] \cdot {e_{i0}} \cdot {{\left( {\sigma _3^\prime } \right)}^n}}} + \frac{{{R_f} \cdot {\varepsilon _a}}}{{\frac{{2 \cdot \cos \varphi }}{{1 - \sin \varphi }}{c_0} + \alpha \cdot S_h^\beta + \frac{{2 \cdot \cos \varphi }}{{1 - \sin \varphi }}\sigma _3^\prime }}}} $$ (7)

    其中,破坏比Rf取0.86(取值范围为0.7~0.9),其他参数取值参见表 1表 2

    由公式(7)可知,含水合物沉积介质的应力-应变关系受水合物饱和度和有效围压的影响,但试验结果表明,应力-应变关系还受剪切速率和温度的影响,该模型并未体现出二者的影响。因此模型还有进一步改进的空间。

    为了说明本文给出的细砂质含水合物沉积介质非线性弹性力学模型的有效性,将模型的计算结果与试验结果进行了对比(图 2图 3),发现模型计算值与试验值符合度较好。

    图 2图 3可以看出本文给出的非线性弹性模型不仅能够描述细砂质含水合物沉积介质在常规三轴剪切条件下的非线性弹性应力-应变关系,还能够描述在水合物饱和度和有效围压同时改变的条线下,沉积介质应力-应变的非线性弹性关系。但是本文给出的模型和参数只能描述符合双曲线特征的应力-应变关系,无法描述应变软化特征。

    (1) 以细砂质含水合物沉积介质力学三轴试验数据为基础,结合前人的研究成果,给出适合细砂质含水合物沉积介质的非线性弹性模型及其参数。该模型形式简单,参数物理意义明确。通过将模型的计算结果与试验结果进行对比验证了模型的适用性。

    (2) 通过理论推导给出了计算破坏比的公式,但该公式的适用性还需要更多的试验数据进一步的检验。初始弹性模量的表达式存在量纲不一致的问题,未来可以做进一步修正。

    (3) 给出的模型中含水合物沉积介质应力-应变关系只受饱和度和有效围压的影响,但应力-应变关系还受到剪切速率和温度等因素的影响,后续研究将对模型进行进一步改进,使其能考虑更多的影响因素。

  • [1] 蔡峰.北黄海盆地中新生代沉积环境及特征[J].海洋地质动态,1997(8):1-4.[CAI Feng. Depositional environment and feature of Cenozoic basin in North Yellow Sea[J]. Marine Geology Letters, 1997

    (8):1-4.]

    [2] 杜旭东,陆克政,漆家福.黄骅坳陷褶皱分类与成因分析[J].石油大学学报自然科学版,1999,23(1):1-6.

    [DU Xu-dong, LU Ke-zheng, QI Jia-fu.The causes and classification of fold in Huanghua Depression[J]. Journal of China University of Petroleum(Edition of Natural Science), 1999, 23(1):1-6.]

    [3] 姚超,焦贵浩,王同和,等.含油气构造样式[M].北京:石油工业出版社,2004.[YAO Chao, JIAO Gui-hao, WANG Tong-he, et al. Oil Gas Structure Styles in China[M].Beijing:Petroleum Industry Press, 2004.]
    [4] 陆克政,朱筱敏,漆家福.含气盆地分析[M].石油大学出版社,2001.[LU Ke-zheng, ZHU Xiao-min, QI Jia-fu.Petroliferous Basin Analysis[M]. China University of Petroleum Press, 2001.]
    [5] 沈华,刘俊,刘建党.贝尔凹陷褶皱构造样式及成因分析[J].中国矿业,2005,14(3):79-11.

    [SHEN Hua, LIU Jun, LIU Jian-dang. Structure styles and origin of fold in Bei'er depression[J].China University of Mining and Technology Press, 2005, 14(3):79-11.]

  • 期刊类型引用(3)

    1. 赵亚鹏,刘乐乐,孔亮,刘昌岭,吴能友. 含天然气水合物土微观力学特性研究进展. 力学学报. 2021(08): 2119-2140 . 百度学术
    2. 李淑霞,郭尚平,陈月明,张宁涛,武迪迪. 天然气水合物开发多物理场特征及耦合渗流研究进展与建议. 力学学报. 2020(03): 828-842 . 百度学术
    3. 董林,廖华林,李彦龙,刘昌岭. 天然气水合物沉积物力学性质测试与评价. 海洋地质前沿. 2020(09): 34-43 . 百度学术

    其他类型引用(2)

计量
  • 文章访问数:  1869
  • HTML全文浏览量:  162
  • PDF下载量:  13
  • 被引次数: 5
出版历程
  • 收稿日期:  2006-03-29
  • 修回日期:  2006-09-29

目录

/

返回文章
返回