Zengmu Xi fault on the western-border of South China Sea and its control over the sedimentation of the Basin
-
摘要: 通过解释重新处理的曾母盆地多道地震数据，结合重、磁异常资料，对前人提出南海西缘最南部发育近S-N向“曾母西断裂”的分布区域进行了重新厘定，详细解释了曾母西断裂的构造特征，具体分析了其结构、组成和空间展布特征。确定曾母西断裂由①号断层和②号断层所构成，曾母西①号断层位于4°30′～7°N、109°30′～110°E，主断层面倾向ESE，走向NEN，其为康西坳陷的大型控边断层；曾母西②号断层分布于3°～5°N、109°30′～110°30′E，主断层面倾向NE，走向NNW，其为塔陶垒堑的大型控边断层，曾母西①号断裂往南的延伸终止于曾母西②号断裂。同时分析讨论了曾母西断裂的构造活动特征以及对曾母盆地沉积发育的控制作用。Abstract: Reprocessed multi-channel seismic data collected from the Zengmu Basin is interpreted in this paper combined with gravity and magnetic anomaly data. It is confirmed that the “Zengmu Xi Fault” is located in the southernmost part of the western-margin of the South China Sea in a nearly S-N direction. The structural characteristics of the fault, especially its composition and spatial framework, are carefully described and explained. It is found that the Zengmu Xi fault is composed of two faults, the Fault 1 and Fault 2. The Fault 1 is located at 4°30′N～7°00′N, 109°30′～110°E, with a main tendency to ESE, and a strike towards NEN. It is a large boundary controlling fault of the Kangxi Depression, while the fault 2 is distributed between 3° and 5°N, 109°30′E and 110°30′E, with a main fault tendency to NE and a strike towards NNW. It is a large boundary controlling fault of the Tatao horst graben. The southward extension of Zengmu Xi fault 1 ends when it meets the fault 2. Upon the basis, the characteristics of the tectonic activity of the Zengmu Xi fault and its control over the depositional evolution of the Zengmu basin are analyzed and discussed.
-
天然气水合物是在低温、高压条件下，甲烷等气体与水作用形成的笼形晶体化合物[1]，分布广、储量大、能量密度高，是最为重要的替代能源之一。我国陆地和海洋中的水合物资源相当丰富[2-3]。根据地质构造和储层条件，水合物藏可分为4类（Ⅰ—Ⅳ类）[4-5]。如图1所示，其中，Ⅱ类水合物藏由水合物层和下伏水层组成，顶底为非渗透层。从地质学、地球化学及热力学等角度分析，Ⅱ类水合物藏是分布最为广泛的一种类型，最有希望得到大规模开发利用[6]。
从经济和技术角度看，降压法[8]和热激法[9]是实际水合物藏开采最为可行的方式，常用于现场的开采试验[10-14]。但与常规油气藏不同，水合物藏在开采过程中会发生相变，即固态的水合物吸热分解为可动流体（气和水），单一的降压法或热激法的作用效果往往十分有限。Moridis等[15]研究了定流量抽取地层水，以实现Ⅱ类水合物藏储层降压的方法，但产气效果并不理想。Gao等[16]开展的降压分解实验也表明，在单纯降压条件下，Ⅱ类水合物岩心的生产指标欠佳。杨圣文[17]研究了定井底流压结合井筒加热进行Ⅱ类水合物藏开采的方法，分析了不同井底流压下的开采效果，结果表明，降低井底流压能改善水合物藏的开采效果。Moridis等[18]提出了单井条件下，分阶段进行井筒加热和注热水开采Ⅱ类水合物藏的方法，但施工程序较为复杂，实际应用存在一定困难。Reagan等[19]对Moridis等提出的方法进行了敏感性分析，结果表明孔隙度、储层非均质性等对开采效果有明显影响。
实验研究[20]和数值模拟研究[21]都表明，热水驱替法能够结合降压法和热激法的优点，是水合物藏开采的一种有效方法。对Ⅱ类水合物藏而言，将热水驱替和井组[22]结合具有明显优势。鉴于Ⅱ类水合物藏分布的广泛性，且相关开采模拟仍有待加强，本文结合实际水合物藏参数，使用数值模拟方法研究了热水驱替方式下Ⅱ类水合物藏的开采规律，并与降压法的开采效果进行了对比，以加强对水合物资源开发利用的认识。
1. 研究方法
1.1 开采方法
采用五点井网，热水驱替开采Ⅱ类水合物藏的示意图如图2所示。由于水合物层的渗透率很低、注入性较差，中心注入井在下伏水层的上部区域射孔，向储层注入热水；4口生产井定井底流压生产。在生产过程中，生产井近井地带的低压环境和注入井持续的热水注入使得储层中的水合物大量分解，分解得到的气体运移至生产井采出。
1.2 数值模拟方法
由于水合物藏开采试验不多，数值模拟是目前研究水合物藏开采动态规律的主要手段。HydrateResSim（HRS）是专门用于水合物藏开采的开源学术代码[23]，它考虑了水合物藏开采过程中的相变、传热、多相渗流等机理，能够对降压法、热激法等方式下的水合物藏开采进行有效模拟[24]。在HRS中，相变过程采用主变量变换法处理，数学模型的离散采用积分有限差分法，离散得到的非线性方程组使用Newton-Raphson方法迭代求解。
本文在HRS基础上开展模拟研究，研究的水合物假设为单一的甲烷水合物，体系考虑4相3组分，相包括气相（G）、水相（A）、水合物相（H）和冰相（I），这4相均为储层孔隙中的一部分。其中冰相和水合物相为不可流动相，水相与气相为可流动相，其流动遵循达西定律。组分包括甲烷组分（m）、水组分（w）和水合物组分（由甲烷组分和水组分根据水合数表示，本文水合数取值为6）。其中，水相、气相和水合物相中都存在水组分和甲烷组分，冰相中只存在水组分。
组分κ的质量守恒方程如式（1）所示：
∂∂t(∑β=A,G,H,IφSβρβXκβ)+∇⋅(∑β=A,GXκβ→Fβ)=qκ (1) 式中，φ为储层孔隙度；Sβ为相β（=G, A, I, H）的饱和度；ρβ为相β的密度，kg·m−3；
Xκβ 为组分κ（=m, w）在相β中的质量分数；→Fβ 为相β的质量流速，kg·m−2·s−1；qκ为组分κ的源汇项，kg·m−3·s−1。体系的能量守恒方程如式（2）所示：
∂∂t((1−φ)ρRCRT+∑β=A,G,H,IφSβρβHβ+Hd)+∇⋅{−((1−φ)KR+∑β=A,G,H,IφSβKβ))∇T+∑β=A,GHβ→Fβ}=qe (2) 式中，ρR为岩石的密度，kg·m−3；CR为岩石的比热容，J·kg−1·K−1；Hβ为相β的焓，J·kg−1；Hd为水合物的形成/分解焓，J·m−3；KR为岩石的导热系数，W·m−1·K−1；Kβ为相β的导热系数，W·m−1·K−1；qe为热量的源汇项，J·m−3·s−1。
模拟研究时，采用平衡模型描述水合物的分解与形成，采用的相对渗透率模型[23]和毛管力模型[23]分别如式（3）、式（4）所示：
krA=min \begin{aligned}{k_{rG}} = \min \left\{ {{{\left[ {\frac{{{S_G} - {S_{irG}}}}{{1 - {S_{irA}}}}} \right]}^{{n_G}}},1} \right\}\end{aligned} (3) 式中，krA为水相相对渗透率；SA为水相饱和度；SirA为束缚水饱和度；nA为水相相对渗透率递减指数，本文取3.572[25]；krG为气相相对渗透率；SG为气相饱和度；SirG为束缚气饱和度；nG为气相相对渗透率递减指数，本文取3.572[25]。
{p_c} = - {p_{c0}}{\left[ {{{\left( {\frac{{{S_A} - {S_{irA}}}}{{1 - {S_{irA}}}}} \right)}^{ - 1/\lambda }} - 1} \right]^{1 - \lambda }} (4) 式中，pc为气水间毛管力，Pa；pc0为模型参数，本文取2×103 Pa[25]；λ为模型参数，本文取0.45[25]。
2. 研究实例
2.1 数值模拟模型
参考Mallik地区Ⅱ类水合物藏的参数[26]建立基础模型（如图3所示。模型大小为165 m×165 m×90 m，从上到下，依次由顶部非渗透层、水合物层、下伏水层和底部非渗透层组成。采用直角网格系统，x，y，z方向对应的网格数分别为33×33×22。模型参数如表1所示。采用五点井网热水驱替方式进行模拟研究，生产时间为3年。其中，注入井和生产井的井径均为0.1 m。各生产井的射孔范围为[−1 m，15 m]，以3 MPa进行定井底流压生产；注入井的射孔范围为[−1 m，0 m]，以200 t/d的速度注入50 ℃的热水。
表 1 基础模型参数Table 1. Basic model parameters of the Class Ⅱ hydrate reservoir参数 水合物层 下伏水层 顶底非渗透层 厚度/m 15 15 30 绝对渗透率/10−3 μm2 1 000 1 000 0 孔隙度 0.35 0.35 0 水合物饱和度 0.7 0 − 含水饱和度 0.3 1 − 底部初始压力/MPa 10.67 − − 底部初始温度/℃ 13.3 − − 温度梯度/（℃/100 m） 3.0 束缚水饱和度 0.20 − 束缚气饱和度 0.02 − 由于相变的存在，水合物在开采过程中会与周围环境进行大量的热交换，热效应显著。而水合物藏的顶底非渗透层虽然几乎没有渗透性，但通常赋存着大量热量，在开采过程中能够以热传导的方式为水合物的分解提供能量[27]，因而对水合物藏的开采效果有直接影响。在本文的模型中，顶底非渗透层均为30 m，这个厚度足以刻画开采过程中的热效应[28]。
2.2 开采规律
2.2.1 产气动态
产气速率和累产气、分解气速率和累分解气曲线分别如图4和图5所示。从图4可以看出，产气速率可大致分为两个阶段（阶段①和阶段②），呈现先快速上升，然后以较快速度下降至趋于相对稳定的变化规律。在阶段①，经50 d达到峰值产气速率约为7.0×104 m3/d。在阶段②，产气速率先以较快速度降低，然后逐渐趋于稳定，在生产末期，产气速率约为5×103 m3/d。分解气速率的变化与产气速率变化相同（图5）。对应产气速率和分解气速率的变化，累产气和累分解气前期快速上升，而后近似线性增加。截至生产结束，累产气1.31×107 m3，累分解气1.33×107 m3。除极少量溶于地层水中的气体外，Ⅱ类水合藏中的采出气全部来源于水合物的分解，因此，分解气几乎被全部采出（＞99%），地层中的残余极少。
出现上述变化的原因主要是：在阶段①，由于定压生产以及水层中水的产出，储层压力下降较快，生产井近井地带的水合物快速分解。同时，热水的注入也导致储层水合物的大量分解，但产出的气体主要来自生产井近井地带水合物的分解。在阶段②，随生产时间的增加，在生产井的降压作用和注入井的热水作用（图6）下，水合物持续分解，但由于水合物饱和度持续下降，总体产气速率逐渐降低。
储层温度场的变化如图6所示（剖面图）。可以看出，随生产时间的增加，水合物的分解导致温度降低，模型的温度整体呈逐渐下降趋势，但注入井近井地带由于热水的注入，高温区域不断扩大，能够持续为水合物的分解提供能量。
2.2.2 水合物饱和度
模型水合物饱和度场的变化如图7所示（剖面图）。可以看出，随生产时间的增加，在生产井的降压作用和注入井的热水作用下，生产井周围、水合物层与下伏水层接触面、注入井周围的水合物饱和度逐渐降低；至生产末期，相当一部分水合物已经分解，剩余的水合物主要存在于井间地带。顶部非渗透层携带的热量也促进了水合物的分解，因而水合物层顶部（即顶部非渗透层与水合物层的交界处）的水合物分解明显。当生产结束时，水合物层的平均水合物饱和度降至0.176。
![图 7 水合物饱和度场]()
图 7 水合物饱和度场注：为展示水合物层水合物饱和度的变化，示意图中只显示了部分网格。Figure 7. Hydrate saturation field evolution2.3 开采效果对比
鉴于降压法在水合物藏开采方法中具有重要地位，为对比热水驱替开采水合物藏的效果，建立了降压法开采的对比模型。在降压法开采模型中，将中心注入井关闭，生产井以定井底流压方式生产（3 MPa），其他储层参数及生产参数不变。
在热水驱替和降压开采条件下，水合物藏的气体采出程度、水合物分解程度及累积气水比如图8所示。其中，气体采出程度为采出气量与气体总储量的比值，如式（5）所示：
![图 8 气体采出程度、水合物分解程度和累积气水比]()
图 8 气体采出程度、水合物分解程度和累积气水比Figure 8. Gas recovery percent, hydrate dissociation percent and gas-water ratio\eta = \frac{{{N_G}}}{{{R_G}}} (5) 式中，η为气体采出程度；NG为采出的甲烷量，m3；RG为水合物藏中甲烷初始总储量（包括水相中的溶解气以及水合物中蕴藏的甲烷总量），m3。
水合物分解程度为分解的水合物量与水合物总储量的比值，如式（6）所示：
\beta = \frac{{{D_H}}}{{{R_H}}} (6) 式中，β为水合物分解程度；DH为分解的水合物量，kg；RH为水合物藏中的水合物总储量，kg。
累积气水比为累积产气量与累积产水量的比值，如式（7）所示：
r = \frac{{{V_G}}}{{{m_A}}} (7) 式中，r为累积气水比；VG为累积产气量，m3；mA为累积产水量，t。
从图8可以看出，截至生产结束，相对于降压法，热水驱替的气体采出程度和水合物分解程度更高（分别高出约40%和24%），且热水驱替条件下的气体采出程度和水合物分解程度均大于60%，处于较高水平。两种模型的差异来源于外源热水的注入。相对于降压模型，热水驱替模型中热水的注入扩大了水合物的分解范围，携带的热量提高了水合物的分解程度，相应地气体采出程度也较高；在降压和热水的综合作用下，开采指标优于降压模型。因此，热水驱替对Ⅱ类水合物藏的开采具有一定的适应性，开采指标较好。但热水驱替的产水量较高，其累积气水比仅为降压开采的1/3，这是由外源热水的注入导致的不可避免的结果，因而开采过程中的水处理是十分重要的问题。
3. 结论
（1）热水驱替开采Ⅱ类水合物藏时，水合物藏的产气速率和分解气速率首先快速上升，然后以较快速度下降至趋于相对稳定；累产气和累分解气先快速上升，然后近似线性增加。气体采出程度和水合物分解程度处于较高水平（＞60%）。
（2）热水驱替对Ⅱ类水合物藏的开采具有一定的适应性。相对于降压开采，在热水驱替条件下，水合物藏的气体采出程度更高，水合物饱和度的降低更显著；但累积气水比较低，伴随较大的产水量。
-
![]()
图 3 曾母盆地地震反射界面解释的地层结构
地震测线位置见图2。
Figure 3. A selected seismic profile showing the Cenozoic structure of the Zengmu Basin
![]()
图 5 曾母西断裂布格重力异常特征
Figure 5. Characteristics of Bouguer gravity anomaly of the Zengmu Xi fault
![]()
图 6 曾母西断裂磁力异常特征
Figure 6. Characteristics of Magnetic gravity anomaly of the Zengmu Xi fault
-
[1] 邱燕, 王立飞, 黄文凯, 等. 中国海域中新生代沉积盆地[M]. 北京: 中国地质出版社, 2016. QIU Yan, WANG Lifei, HUANG Wenkai, et al. Sedimentary Basins in Mesozoic and Cenozoic in the China Sea[M]. Beijing: Geological Publishing House, 2016.
[2] Hall R. Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations [J]. Journal of Asian Earth Sciences, 2002, 20(4): 353-431. doi: 10.1016/S1367-9120(01)00069-4
[3] 李文勇, 李东旭. 中国南海不同板块边缘沉积盆地构造特征[J]. 现代地质, 2006, 20(1):19-29. [LI Wenyong, LI Dongxu. Tectonic characteristics on the sedimentary Basins with different plate margins in the South China Sea [J]. Geoscience, 2006, 20(1): 19-29. doi: 10.3969/j.issn.1000-8527.2006.01.002 [4] Golonka J, Krobicki M, Paja K J, et al. Global Plate Tectonics and Paleogeography of Southeast Asia[M]. Arkadia, Kraków: Ministry of the Environment, 2006: 119-127.
[5] 刘宝明, 夏斌, 李绪宣, 等. 红河断裂带东南的延伸及其构造演化意义[J]. 中国科学 D辑: 地球科学, 2006, 36(10):914-924. [LIU Baoming, XIA Bin, LI Xuxuan, et al. Southeastward extension of the Red River Fault Zone and implication for its tectonic evolution significance [J]. Science in China (Series D) Earth Sciences, 2006, 36(10): 914-924. [6] 高红芳. 南海西缘断裂带走滑特征及其形成机理初步研究[J]. 中国地质, 2011, 38(3):537-543. [GAO Hongfang. A tentative discussion on strike-slipping character and formation mechanism of western-edge fault belt in South China Sea [J]. Geology in China, 2011, 38(3): 537-543. doi: 10.3969/j.issn.1000-3657.2011.03.003 [7] 陈强, 金庆焕. 南海西部万安断裂活动特征及其对万安盆地的控制作用[J]. 海洋地质前沿, 2017, 33(10):1-8. [CHEN Qiang, JIN Qinghuan. Activity features of Wan’an fault and its constraints on Wan’an basin, western south China Sea [J]. Marine Geology Frontiers, 2017, 33(10): 1-8. [8] 邱燕. 南沙海域曾母盆地构造样式特征及其与油气聚集和圈闭的关系[J]. 地质通报, 2005, 24(1):16-22. [QIU Yan. Typical structural styles and their relationships with hydrocarbon accumulation and traps in the Zengmu Basin, Nansha sea area [J]. Geological Bulletin of China, 2005, 24(1): 16-22. doi: 10.3969/j.issn.1671-2552.2005.01.003 [9] 邱燕, 陈国能, 解习农, 等. 南海西南海域曾母盆地新生界沉积充填演化研究[J]. 热带海洋学报, 2005, 24(5):43-52. [QIU Yan, CHEN Guoneng, XIE Xinong, et al. Sedimentary filling evolution of Cenozoic strata in Zengmu Basin, Southwestern South China Sea [J]. Journal of Tropical Oceanography, 2005, 24(5): 43-52. doi: 10.3969/j.issn.1009-5470.2005.05.005 [10] 陈玲. 南沙海域曾母盆地西部地质构造特征[J]. 石油地球物理勘探, 2002, 37(4):354-362. [CHEN Ling. Geologic structural feature in the west of Zengmu Basin, Nansha sea area [J]. Oil Geophysical Prospecting, 2002, 37(4): 354-362. doi: 10.3321/j.issn:1000-7210.2002.04.010 [11] 郭佳, 谢晓军, 刘世翔, 等. 南海曾母盆地新生代沉积体系特征[J]. 中国石油勘探, 2016, 21(4):99-107. [GUO Jia, XIE Xiaojun, LIU Shixiang, et al. Cenozoic sedimentary systems in Zengmu Basin, South China Sea [J]. China Petroleum Exploration, 2016, 21(4): 99-107. doi: 10.3969/j.issn.1672-7703.2016.04.0011 [12] 姜绍仁, 周效中, 黄慈流, 等. 南海曾母盆地地质构造特征及演化历史[J]. 海洋地质与第四纪地质, 1990, 10(1):25-36. [JIANG Shaoren, ZHOU Xiaozhong, HUANG Ciliu, et al. Zengmu Basin: geologic structure and evolution [J]. Marine Geology and Quaternary Geology, 1990, 10(1): 25-36. [13] Hutchison C S. Marginal basin evolution: the southern South China Sea [J]. Marine and Petroleum Geology, 2004, 21(9): 1129-1148. doi: 10.1016/j.marpetgeo.2004.07.002
[14] 钟广见, 王嘹亮. 南海西南部新生代盆地类型及演化历史[J]. 海洋湖沼通报, 1996(2):24-32. [ZHONG Guangjian, WANG Liaoliang. The types of Cenozoic basin in the southwest of South China Sea and their evolution [J]. Transactions of Oceanology and Limnology, 1996(2): 24-32. [15] 林珍. 南沙中部海域北康·曾母盆地重磁异常特征及解释[J]. 物探与化探, 2003, 27(4):263-268. [LIN Zhen. Characteristics and interpretation of gravity and magnetic anomalies in the Beikang Zengmu Basin of central Nansha islands [J]. Geophysical and Geochemical Exploration, 2003, 27(4): 263-268. doi: 10.3969/j.issn.1000-8918.2003.04.005 [16] 安慧婷, 李三忠, 索艳慧, 等. 南海西部新生代控盆断裂及盆地群成因[J]. 海洋地质与第四纪地质, 2012, 32(6):95-111. [AN Huiting, LI Sanzhong, SUO Yanhui, et al. Basin-controlling faults and formation mechanism of the Cenozoic basin groups in the western south China Sea [J]. Marine Geology & Quaternary Geology, 2012, 32(6): 95-111. [17] 张功成, 谢晓军, 王万银, 等. 中国南海含油气盆地构造类型及勘探潜力[J]. 石油学报, 2013, 34(4):611-627. [ZHANG Gongcheng, XIE Xiaojun, WANG Wanyin, et al. Tectonic types of petroliferous basins and its exploration potential in the South China Sea [J]. Acta Petrolei Sinica, 2013, 34(4): 611-627. doi: 10.7623/syxb201304001 [18] 姚伯初, 万玲, 刘振湖, 等. 南海南部海域新生代万安运动的构造意义及其油气资源效应[J]. 海洋地质与第四纪地质, 2004, 24(1):69-77. [YAO Bochu, WAN Ling, LIU Zhenhu, et al. Tectonic significance and its petroleum effect of the Wan’an tectonic movement in the south of the South China Sea [J]. Marine Geology & Quaternary Geology, 2004, 24(1): 69-77. [19] 蔡利飘, 周娟. 银川盆地断裂体系发育特征及其对盆地的控制作用[J]. 油气地球物理, 2018, 16(2):53-59. [CAI Lipiao, ZHOU Juan. Development characteristics of the fault system in Yinchuan Basin and its control action on the Basin [J]. Petroleum Geophysics, 2018, 16(2): 53-59. [20] 朱晨曦. 南海南部断裂及前新生代基底岩性分布特征[D]. 吉林大学硕士学位论文, 2019. ZHU Chenxi. Characteristics of fault and distribution of Pre-Cenozoic basement lithology in Southern South China Sea[D]. Master Dissertation of Jilin University, 2019.
[21] 刘海龄, 杨树康, 周蒂, 等. 南沙北部伸展构造的基本特征及其动力学意义[J]. 高校地质学报, 1998, 4(1):64-72. [LIU Hailing, YANG Shukang, ZHOU Di, et al. Basic characteristiecs of extension structure in northern Nansha islands, China, and its dynamical implications [J]. Geological Journal of China Universities, 1998, 4(1): 64-72. [22] 杨木壮, 吴进民. 曾母盆地地层发育特征[J]. 海洋地质与第四纪地质, 1996, 17(1):18-27. [YANG Muzhuang, WU Jinmin. The stratigraphy of Zengmu Basin [J]. Marine Geology and Quaternary Geology, 1996, 17(1): 18-27. -
期刊类型引用(4)
1. 邱燕，杜文波，鞠东，韩冰，马浩明，胡小三. 南海主要断裂（带）与构造分区. 海洋地质与第四纪地质. 2025(01): 136-153 . 
本站查看
2. 鞠东，邱燕，黄文凯，杜文波. 南海主要控凹断层活动特征与新生代地壳初始张裂. 海洋地质前沿. 2022(09): 59-70 . 百度学术
3. 杜文波，朱荣伟，熊量莉，朱红涛，舒誉，黄众. 珠江口盆地恩平凹陷恩平组物源体系分析及未来大洋钻探建议. 地质学报. 2022(08): 2788-2799 . 百度学术
4. 李法坤，戴黎明，杜晓东，蔡国富，李三忠，董昊，王宇. 构造-沉积耦合过程的数值模拟：以南海北部阳江凹陷为例. 海洋地质与第四纪地质. 2021(05): 139-150 . 本站查看
其他类型引用(0)
![图 1 水合物藏类型示意图（据文献[7]修改）](/fileHYDZYDSJDZ/journal/article/hydzydsjdz/2021/2/2018102603-1.jpg)
下载:
