Identification of submarine fan in the Meishan Formation of the Ledong-Lingshui Sag in the Qiongdongnan Basin and its depositional model
-
摘要: 近年来,乐东–陵水凹陷新近系岩性圈闭成为南海海域勘探的热点,发现了多个大气田,且已发现的岩性圈闭均与海底扇、浊积水道密切相关,为了获得更多天然气发现,提高岩性圈闭的勘探成功率,应加强海底扇识别及沉积模式研究。本次研究基于乐东–陵水凹陷的钻井和三维地震资料,以经典层序地层学为指导,从钻井取心入手,综合运用测井、FMI成像测井及分析化验等资料,分析了乐东–陵水凹陷梅山组海底扇的岩相、沉积构造、测井相及地震相特征。梅山组海底扇岩性以粉砂岩为主,发育块状构造、小型交错层理和冲刷面,结合粒度、C-M图分析认为其为浊流沉积,后期受牵引流改造。在梅山组地层内识别出三种类型海底扇:水道化型、盆底扇型、底流改造型海底扇,可进一步划分为3个期次,各期海底扇形成演化、主物源方向受古地貌影响由陵水凹陷向乐东凹陷迁移,扇体规模逐渐变大。综合考虑坡折带、相对海平面变化、三角洲物源供给对沉积的控制,建立了乐东–陵水凹陷梅山组“源–沟–坡–扇”的海底扇沉积模式。Abstract: The Neogene lithologic traps in the Ledong-Lingshui Sag have become hotspots for exploration of oil and gas in the South China Sea in recent years. Some large gas fields are discovered closely related to submarine fans and turbidite deposits. In order to obtain more natural gas discoveries to enhance the exploration success rate in the lithologic traps, it has become an urgent task to researchers to identify the submarine fans and to establish their depositional model. With the support of drilling and 3D seismic data from the Ledong-Lingshui Sag and following the principles of sequence stratigraphy, we started researches on the cores from drilling holes. Lithofacies, logging facies and seismic facies as well as sedimentary structures of the Meishan Formation submarine fans in the Ledong-Lingshui sag are carefully studied using well logging, FMI imaging logging and chemical data as tools. The submarine fans of the Meishan Formation is dominated by massive siltstone, with small cross bedding and scouring surfaces, which is recognized as turbidity current deposits later modified by tractive current based on grain size distribution pattern such as C-M diagram. Three types of submarine fans, say, channel type, basin floor type, and underflow reformed type, are defined. The deposits may be further divided into three periods, and the formation and evolution of the submarine fans in each period as well as the main source direction are affected by paleotopograpic framework. With the fans moving from the Lingshui sag to the Ledong sag, the scale of fans increased gradually. Upon the consideration of the type and distribution of slope breaks, relative sea level changes, deltaic source supply, and paleogeomorphological control over sedimentation, a depositional model including source, channel, slope and fan is established for the Meishan Formation in the Ledong-Lingshui Sag.
-
天然气水合物是在低温、高压条件下,甲烷等气体与水作用形成的笼形晶体化合物[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)所示:
$$\frac{\partial }{{\partial t}}\left(\sum\limits_{\beta = A,G,H,I} {\varphi {S_\beta }{\rho _\beta }X_\beta ^\kappa } \right) + \nabla \cdot \left(\sum\limits_{\beta = A,G} {X_\beta ^\kappa {{\vec F}_\beta }}\right) = {q^\kappa }$$ (1) 式中,φ为储层孔隙度;Sβ为相β(=G, A, I, H)的饱和度;ρβ为相β的密度,kg·m−3;
$X_\beta ^\kappa $ 为组分κ(=m, w)在相β中的质量分数;${\vec F_\beta }$ 为相β的质量流速,kg·m−2·s−1;qκ为组分κ的源汇项,kg·m−3·s−1。体系的能量守恒方程如式(2)所示:
$$\begin{split}&\frac{\partial }{{\partial t}}\left( {(1 - \varphi ){\rho _R}{C_R}T + \sum\limits_{\beta = A,G,H,I} {\varphi {S_\beta }{\rho _\beta }{H_\beta }} + {H_d}} \right) + \nabla \cdot \\ & \left\{ { - \left( {(1 - \varphi ){K_R} + \sum\limits_{\beta = A,G,H,I} {\varphi {S_\beta }{K_\beta })} } \right)\nabla T + \sum\limits_{\beta = A,G} {{H_\beta }{{\vec F}_\beta }} } \right\} = {q^e}\end{split}$$ (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)所示:
$${k_{rA}} = \min \left\{ {{{\left[ {\frac{{{S_A} - {S_{irA}}}}{{1 - {S_{irA}}}}} \right]}^{{n_A}}},1} \right\}$$ $$\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。
2.3 开采效果对比
鉴于降压法在水合物藏开采方法中具有重要地位,为对比热水驱替开采水合物藏的效果,建立了降压法开采的对比模型。在降压法开采模型中,将中心注入井关闭,生产井以定井底流压方式生产(3 MPa),其他储层参数及生产参数不变。
在热水驱替和降压开采条件下,水合物藏的气体采出程度、水合物分解程度及累积气水比如图8所示。其中,气体采出程度为采出气量与气体总储量的比值,如式(5)所示:
$$\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)热水驱替对Ⅱ类水合物藏的开采具有一定的适应性。相对于降压开采,在热水驱替条件下,水合物藏的气体采出程度更高,水合物饱和度的降低更显著;但累积气水比较低,伴随较大的产水量。
-
-
[1] Forel F A. Les ravins sous-lacustres des fleuves glaciaires [J]. Comptes Rendus de l’Académie des Sciences Paris, 1885, 101: 725-728.
[2] Bouma A H. Sedimentology of Some Flysch Deposits: A Graphic Approach to Facies Interpretation[M]. Amsterdam: Elsevier, 1962: 168.
[3] Vail P R, Mitchum R M, Thompson S. Seismic stratigraphy and global changes of sea level[M]//Payton C E. Seismic Stratigraphy—Applications to Hydrocarbon Exploration. Tulsa: AAPG Memoir, 1977: 63-81.
[4] Bouma A H. Key controls on the characteristics of turbidite systems[M]//Lomas S A, Joseph P. Confined Turbidite Systems. London: Geological Society, 2004: 9-22.
[5] Muntingh A, Brown Jr L F. Sequence stratigraphy of petroleum plays, post-rift cretaceous rocks (Lower Aptian to Upper Maastrichtian), Orange Basin, western offshore, South Africa[C]//Weimer P, Posamentier H. Siliciclastic Sequence Stratigraphy: Recent Developments and Applications. Tulsa: AAPG, 1993: 71-98.
[6] 熊波, 苏龙, 张功成, 等. 琼东南盆地深水区崖城组煤系烃源岩定量评价[J]. 天然气地球科学, 2014, 25(4):542-550. [XIONG Bo, SU Long, ZHANG Gongcheng, et al. Quantitative evaluation for coal-bearing source rocks of Yacheng Formation in the deepwater area of the Qiongdongnan Basin [J]. Natural Gas Geoscience, 2014, 25(4): 542-550. [7] 张亚雄, 朱筱敏, 张功成, 等. 中国南海琼东南盆地渐新统陵水组沉积特征[J]. 天然气地球科学, 2013, 24(5):956-964. [ZHANG Yaxiong, ZHU Xiaomin, ZHANG Gongcheng, et al. Sedimentary characteristics of Oligocene Lingshui Formation in Qiongdongnan Basin, South China Sea [J]. Natural Gas Geoscience, 2013, 24(5): 956-964. [8] 吴玉坤, 胡明毅, 刘志峰, 等. 琼东南盆地崖北凹陷崖城组沉积相及煤系烃源岩分布[J]. 天然气地球科学, 2013, 24(3):582-590. [WU Yukun, HU Mingyi, LIU Zhifeng, et al. The sedimentary facies and coal measures source rock distribution of Yacheng Formation in Yabei sag, Qiongdongnan Basin [J]. Natural Gas Geoscience, 2013, 24(3): 582-590. [9] 王子嵩, 刘震, 黄保家, 等. 琼东南盆地深水区中央坳陷带东部渐新统烃源岩分布及评价[J]. 天然气地球科学, 2014, 25(3):360-371. [WANG Zisong, LIU Zhen, HUANG Baojia, et al. Distribution and evaluation of Oligocene source rocks in the east of central depression belt in deep-water area, Qiongdongnan Basin [J]. Natural Gas Geoscience, 2014, 25(3): 360-371. [10] 左倩媚, 张道军, 王亚辉, 等. 琼东南盆地深水区新近系海底扇沉积特征与资源潜力[J]. 海洋学报, 2016, 38(11):105-116. [ZUO Qianmei, ZHANG Daojun, WANG Yahui, et al. Sedimentary characteristics and exploration potential of Neogene submarine fan in the Deepwater area of the Qiongdongnan Basin [J]. Haiyang Xuebao, 2016, 38(11): 105-116. [11] 田姗姗, 苏明, 何云龙, 等. 琼东南盆地新近系坡折带特征及其对沉积体系的控制[J]. 科技导报, 2010, 28(10):73-78. [TIAN Shanshan, SU Ming, HE Yunlong, et al. Patterns of Neogene slope-break zones and the depositional system under their control in Qiongdongnan Basin [J]. Science & Technology Review, 2010, 28(10): 73-78. [12] 龙根元, 吴世敏, 刘兵, 等. 琼东南盆地南部坳陷带12号断层性质及成因[J]. 海洋地质与第四纪地质, 2016, 36(2):93-99. [LONG Genyuan, WU Shimin, LIU Bing, et al. Features of NO.12 fault and its forming mechanism in the southern depression of Qiongdongnan Basin [J]. Marine Geology & Quaternary Geology, 2016, 36(2): 93-99. [13] 胡杨, 张伟, 张景茹, 等. 南海北部大陆边缘琼东南盆地有利油气聚集带及勘探方向[J]. 海洋地质与第四纪地质, 2015, 35(4):73-86. [HU Yang, ZHANG Wei, ZHANG Jingru, et al. The favorable petroleum accumulation belts and exploration targets in Qiongdongnan Basin on the margin of northern South China Sea [J]. Marine Geology & Quaternary Geology, 2015, 35(4): 73-86. [14] 杨军, 施小斌, 王振峰, 等. 琼东南盆地张裂期沉降亏损与裂后期快速沉降成因[J]. 海洋地质与第四纪地质, 2015, 35(1):81-90. [YANG Jun, SHI Xiaobin, WANG Zhenfeng, et al. Origin of syn-rift subsidence deficit and rapid post-rift subsidence in Qiongdongnan Basin [J]. Marine Geology & Quaternary Geology, 2015, 35(1): 81-90. [15] 赵民, 张晓宝, 吉利明, 等. 琼东南盆地构造演化特征及其对油气藏的控制浅析[J]. 天然气地球科学, 2010, 21(3):494-502. [ZHAO Min, ZHANG Xiaobao, JI Liming, et al. Characteristics of tectonic evolution in the Qiongdongnan Basin and brief discussion about its controlling on reservoirs [J]. Natural Gas Geoscience, 2010, 21(3): 494-502. [16] 谢文彦, 张一伟, 孙珍, 等. 琼东南盆地断裂构造与成因机制[J]. 海洋地质与第四纪地质, 2007, 27(1):71-78. [XIE Wenyan, ZHANG Yiwei, SUN Zhen, et al. Characteristics and formation mechanism of faults in Qiongdongnan Basin [J]. Marine Geology & Quaternary Geology, 2007, 27(1): 71-78. [17] 范彩伟, 李绪深, 刘昆, 等. 琼东南盆地乐东、陵水凹陷中新统岩性地层圈闭成藏条件[J]. 中国海上油气, 2016, 28(2):53-59. [FAN Caiwei, LI Xushen, LIU Kun, et al. Hydrocarbon accumulation condition of Miocene litho-stratigraphic trap in Ledong & Lingshui sags, Qiongdongnan Basin [J]. China Offshore Oil and Gas, 2016, 28(2): 53-59. [18] 冯杨伟, 屈红军, 张功成, 等. 南海北部琼东南盆地深水区梅山组一段地震相分析[J]. 矿物岩石, 2016, 36(1):82-95. [FENG Yangwei, QU Hongjun, ZHANG Gongcheng, et al. The seismic facies analyzing of deep water area of Miocene Meishan formation in Qiongdongnan Basin, northern South China Sea [J]. Journal of Mineralogy and Petrology, 2016, 36(1): 82-95. [19] 王元. 莺琼盆地烃源岩生源构成、生烃潜力及其控藏作用研究[D]. 中国矿业大学(北京)博士学位论文, 2018. WANG Yuan. Source composition, hydrocarbon generation potential of source rocks and its control on hydrocarbon accumulation in Yingqiong Basin[D]. Doctor Dissertation of China University of Mining and Technology (Beijing), 2018.
[20] 赵忠泉, 钟广见, 冯常茂, 等. 南海北部西沙海槽盆地新生代层序地层及地震相[J]. 海洋地质与第四纪地质, 2016, 36(1):15-26. [ZHAO Zhongquan, ZHONG Guangjian, FENG Changmao, et al. Cenozoic sequence stratigraphy and seismic facies analysis of Xisha Trough Basin in northern South China Sea [J]. Marine Geology & Quaternary Geology, 2016, 36(1): 15-26. [21] 雷振宇, 张莉, 苏明, 等. 南海南部北康盆地中中新世深水沉积体类型、特征及意义[J]. 海洋地质与第四纪地质, 2017, 37(6):110-118. [LEI Zhenyu, ZHANG Li, SU Ming, et al. Middle Miocene deep-water sediments in the Beikang Basin, southern South China Sea: types, characteristics and implications [J]. Marine Geology & Quaternary Geology, 2017, 37(6): 110-118. [22] 张建新, 党亚云, 何小胡, 等. 莺歌海盆地乐东区峡谷水道成因及沉积特征[J]. 海洋地质与第四纪地质, 2015, 35(5):29-36. [ZHANG Jianxin, DANG Yayun, HE Xiaohu, et al. Origin and sedimentary characteristics of canyon channels in Ledong area of Yinggehai Basin [J]. Marine Geology & Quaternary Geology, 2015, 35(5): 29-36. [23] 王远, 严学新, 王治华, 等. 现代长江水下三角洲浅地层地震相特征[J]. 海洋地质与第四纪地质, 2019, 39(2):114-122. [WANG Yuan, YAN Xuexin, WANG Zhihua, et al. Shallow seismic facies characteristics of the modern underwater delta of the Yangtze River [J]. Marine Geology & Quaternary Geology, 2019, 39(2): 114-122.