Application of geostatistical grain size trend analysis in the Luanhe River Subaqueous Delta
-
摘要: 利用均匀分布于滦河水下三角洲的85个表层沉积物粒度数据，基于地统计学的粒径输运趋势方法对滦河水下三角洲沉积物的输运趋势及影响因素进行了研究。半方差函数分析表明，以0.015°（地球大圆弧度）为插值半径加密后（规则分布）的平均粒径半方差值具有随变程增加而增加的趋势，并在变程值达到0.09°后趋于稳定，而原始数据的半方差值与变程相关性较差。以此（0.09°）为特征距离计算的粒径输运结果显示，滦河水下三角洲南部水深12 m以浅区域粒径输运趋势以向南为主，北部具有向西和向北的趋势；水深12～15 m范围内粒径输运趋势近似西南方向，且趋势较大；南部水深15 m以深地区粒径输运趋势方向为西北向，北部为东南向。粒径输运趋势整体上与实测潮流、余流和泥沙输运方向一致，显示潮流对研究区沉积物的输运具有控制作用。滦河水下三角洲地区沉积环境和输运趋势受河流、波浪、潮流、残留沉积和地形等多种因素影响，粒径输运趋势仅能够较好地解释潮流影响，对于残留沉积等复杂沉积过程的解释则相对有限。Abstract: With the grain size data of 85 surface sediment samples collected from the Luanhe River Subaqueous Delta (LRSD), we studied the sediment transportation trend using the geostatistical grain size trend analysis method (GSTA) in this paper. It is revealed that the semi-variance value of mean grain size in a regular grid increases with distance until the distance reaches 0.09 decimal degree, while the irregular data shows no such a correlation with distance. Thus, the decimal degree of 0.09 is chosen as the character distance for grain size trend analysis. It is further revealed that the direction of the semi-variance is 55°, which is parallel to the tidal current, suggesting that the character distance is tide related. Geostatistical GSTA result further indicates that the net sediment transport trend direction in the south area, where the water depth is less than 12 m, is in southwest, whereas that in the north area, where the depth is less than 12 m, is in west and southwest directions. In the area with water depth between 12 and 15 m, the net sediment transport trend is in southwestern direction, roughly parallel to the coast line. In the area under water depth more than 15 m, the net sediment transport trend is northwest in direction in the southern part, but southeast direction in the northern part. The grain size net sediment transport trends mentioned above are similar to the field measured tidal current, residual current, and sediment transport directions suggesting a tide predominated sediment transportation. However, relict sediments, wave, river and bedform are more complex, which may give influences the grain size net sediment transport trend. The research results of the net sediment transport trend in LRSD may reveal the influence of tidal current. However, for an explanation to relict sediments, further researches are required.
-
琼东南盆地是中国近海最重要的天然气勘探领域，中新统浊积体是其主要的勘探开发目标，探明储量主要集中在黄流组中央峡谷浊积砂体，而梅山组和三亚组浊积扇研究薄弱、勘探程度低[1-3]。已有研究表明，琼东南盆地梅山组亦广泛发育浊积扇，截至目前，已有20余口探井钻遇浊积扇，且钻探发现有逐年增加的趋势，显示了良好的勘探潜力。然而，梅山组浊积扇存在储层致密、泥质含量高、砂层薄、顶油底水、油藏封闭性差等问题，储层预测成功率低，储层品质及有效性成为制约其勘探成功的瓶颈问题。
本文综合利用钻遇梅山组浊积扇探井的岩芯、井壁取芯及镜下薄片、压汞等资料，分析了梅山组浊积扇储层主控因素，并从岩性侧封、扇体上倾尖灭角度分析了储层成藏有效性。研究成果对梅山组浊积扇进一步勘探具有重要的指导意义，对三亚组浊积扇勘探亦具有重要参考价值。
1. 梅山组浊积扇储层特征
琼东南盆地在中中新世梅山组沉积期处于裂后热沉降演化阶段，沉积背景总体以滨浅海-半深海为特点，发育了2个三级层序（图1）。研究认为每个三级层序的不同体系域均沉积了一系列规模不同、类型多样的海底浊积扇体[2,4-8]，诸如在斜坡脚和盆底深洼区发育的各类盆底扇，斜坡上堆积的斜坡扇和高位期陆棚边缘发育的陆棚扇等[9-12]。综合海底地形地貌、断裂活动、海平面变化等因素，李建平等认为盆底扇可划分为断控缓坡根部盆底扇、断控陡坡根部盆底扇、深洼盆底扇、峡谷切割深洼盆底扇等4种类型[3,13]。钻井资料揭示，这些海底扇沉积物粒度总体较细，以细-粉砂岩为主，泥质粉砂岩次之，少量为中-粗砂岩，局部可见砾岩；埋深总体较大，绝大部分超过了3 000 m，胶结程度较高，构成了琼东南盆地中新统重要的油气储层[2,8]。
![]()
图 1 琼东南盆地梅山组浊积扇平面分布及新近系综合地层柱状图Figure 1. Distribution of submarine fan of the Meishan Formation and stratigraphic sequence of the Neogene in the Qiongdongnan Basin梅山组浊积扇储层岩石类型以长石岩屑砂岩、长石岩屑石英砂岩、岩屑石英砂岩和岩屑砂岩为主，砂岩中偶见砾石颗粒，砾石最大可达5 cm以上。砾石成分复杂，可见花岗岩、灰岩、火山岩、碎屑岩。砂岩沉积结构特征显示，靠近盆地西北段的A-7、A-12井区距离物源区近，沉积物粒度粗（砾石），分选差，磨圆度变化大。靠近盆地东北段的A-1、A-5井区距离物源区相对远，沉积物磨圆度变好，由棱角向圆状过渡（图2）。
![]()
图 2 琼东南盆地梅山组浊积体储层岩石薄片照片Figure 2. Pictures of turbidite thin-sections in the Meishan Formation, Qiongdongnan Basin琼东南盆地梅山组已钻遇海底扇储层物性资料表明，整体上随埋深增加，孔隙度和渗透率降低，局部分布异常高孔、高渗段。孔隙度主要分布在10%～20%，平均值为12.4%，渗透率主要为（0.05～100）×10−3 μm2，平均值为12.6×10−3 μm2，以低孔-中渗储层为主（图3），孔隙度和渗透率相关性差，孔渗测试数据离散，说明影响孔渗的因素多而复杂（图4）。
![]()
图 3 琼东南盆地梅山组浊积体储层物性与埋深关系图Figure 3. Relationship between reservoir physical properties and burial depth of turbidite in the Meishan Formation, Qiongdongnan Basin![]()
图 4 琼东南盆地梅山组浊积体岩心、旋转井壁取心孔渗交会图Figure 4. Cross plot of porosity and permeability of turbidite sandstones from core and rotary sidewall coring hole of the Meishan Formation in the Qiongdongnan Basin琼东南盆地梅山组浊积扇以产出天然气为主，本次研究参考莺琼盆地天然气储层分类标准，从生产实际出发（海上石油天然气储量估算规范，2020），提出了梅山组浊积体天然气储层分类方案，结合中国近海海域天然气单井日产量下限、产能分类，分高产（15×104 m3/d）、具有商业产能（2.5×104 m3/d）、有一定产能、见地层流体、见泥浆滤液，分别对应于优质、好、有效、致密储层（表1）。
表 1 琼东南盆地梅山组天然气储层分类标准Table 1. Classification standard of natural gas reservoir in the Meishan Formation of the Qiongdongnan Basin中石油分类 本次研究分类 储层分类 孔隙度/% 渗透率/(10−3 μm2) 评价 分类 评价 I ＞20 ＞1000 最好 优质储层 高产 II 15～20 100～1000 好 III 10～15 10～100 较好 IV 5～10 1～10 较差 好储层 商业产能 V ＜5 0.1～1 差 有效储层 少量地层流体 ＜0.1 致密层 泥浆滤液 2. 储层主控因素分析
琼东南盆地不同凹陷梅山组浊积扇埋深有较大差异，其浊积扇储层物性除受到埋深影响外，还受到物源、海平面变化、扇体类型及内部构成以及底流改造的影响。
2.1 沉积物来源
沉积物源对储层发育特征的影响主要体现在不同来源沉积物其碎屑成分、矿物组合等可能存在较大差异，进而导致其储层物性存在一定差异[14-16]。琼东南盆地梅山组沉积期发育北部海南岛物源、西部越南物源以及东部神狐物源，利用岩屑成分、重矿物资料以及沉积相平面图可识别出物源差异。
岩屑成分特征显示，A-2井和A-11井岩屑包括了岩浆岩、变质岩、沉积岩。A-2井区岩屑含量呈跳跃分布，A-11井岩屑含量相对稳定，其中喷出岩岩屑较多。反映前者物源环境有不稳定性。A-2井、A-14井和A-11井花岗岩、变质岩和沉积岩屑分布与A-7井有一定相似性，说明可能受同一物源影响（图5）。
![]()
图 5 琼东南盆地梅山组岩石碎屑成分分布图Figure 5. Distribution of rock clastic composition of the Meishan Formation, Qiongdongnan Basin重矿物资料显示，琼东南盆地梅山组沉积期不同构造带的源-汇供给存在差异。A-7井区及以西区域，重矿物组合以锆石+电气石+白钛矿为主，ZTR指数值高，表明A-7井区远离物源，不仅受海南岛物源控制，还受红河物源注入的影响。A-12井区和A-1井区及西南和松东凹陷东南向区域，重矿物组合以锆石+电气石+白钛矿为主，白钛矿含量占据绝对优势。靠近海南隆起区域ZTR数值较A-1井区高，相对靠近物源，与海南岛河流砂岩重矿物特征一致。A-5井区重矿物组合以锆石+电气石+白钛矿为主，ZTR数值占据绝对优势，应属于海南岛物源独立分支，且远离海南岛母岩物源方向（图6）。
![]()
图 6 琼东南盆地中央拗陷西北段梅山组储层重矿物平面分布Figure 6. Plane distribution of heavy minerals in the Meishan Formation reservoir in the northwest of central depression of the Qiongdongnan Basin岩石薄片、岩矿组分和轻重矿物组合特征，以及宏观构造演化和河道发育分析结果表明，A-7井区是海南岛、越南中部和红河多物源供给，沿红河断裂带1号和13号断裂和多期河道，NW至SE向搬运和沉积。A-12和A-1井区是海南岛物源供给，经5号、3号和2号断裂和多期河道，NNW至SSE向搬运和沉积。A-5井区为海南岛不同分支物源，经5号、6号和2号断裂和多期河道，NE至SW向搬运和沉积（图6）。琼东南盆地梅山组沉积物总体受海南岛物源的绝对控制，表现为海南岛、红河和越南中部等多物源供给格局、沿继承性活动断裂和多期发育的河道汇聚、沉积在深洼中央。
物源条件决定了储层的原始碎屑成分及其分选情况，进而从根本上影响和控制储层物性特征。总体上，海南岛物源供给的浊积扇储层的碎屑成分成熟度和结构成熟度都均高于西部越南物源为主供给的浊积扇储层；砂岩中值粒径方面，海南岛物源供给浊积扇储层低于西部越南物源供给浊积扇储层。物源条件决定了梅山组浊积扇储层原始成分和结构的差异，进而影响其物性特征的差异。
2.2 海平面变化
海平面变化对储层发育特征的影响主要反映在海平面变化不同阶段发育的海底扇其储层物性特征存在较大差异。A-7井连续超180 m的取芯段上，以正常深海泥岩为界，可识别出3期海底扇浊积体，钻遇的梅山组海底扇为高位体系域沉积，黄流组海底扇为低位体系域沉积（图7）。
![]()
图 7 A-7井系统取芯段黄流组、梅山组浊积扇物性差异变化Figure 7. Variations in porosity and permeability of turbidite fan in the Huangliu and Meishan Formations in coring section of the well A-7浊积体碎屑颗粒组成复杂，岩性从中细砂岩到砾岩，砾石最大可达5 cm以上；砾石成分上，花岗岩、灰岩、火山岩、碎屑岩均可见，且分选磨圆都很差；间夹异重流沉积、见有片状泥砾、假团块。这些沉积特点增强了储层内部非均质性，是储层变差的重要因素。颗粒粗使其保存有一定数量的大孔隙，因此仍有好储层存在，但储层物性变化很大。其中异重流沉积属于洪水事件重力流沉积，一次性搬运，分选磨圆差。从图7可以看出，低位体系域海底扇储层物性明显优于高位体系域。
低水位晚期，陆坡浊积沟道会有溢岸沉积，颗粒细，泥质含量高，斜坡扇好-优质储层仅占7%；高水位时期，陆架三角洲前端可能发育陆棚扇。由于可能夹有洪水期的异重流沉积，分选差，陆棚扇好-优质储层仅占13%。
2.3 浊积扇类型及内部构成
研究认为不同类型的海底扇其储层物性有较大差别[4,14,17]。盆底扇储层中好-优质储层占68%，但不同类型间差别很大。断控缓坡根部盆底扇好-优质储层在76%以上；而断控陡坡根部盆底扇储层品质总体不好；深洼盆底扇优质储层占比高达85%，但上倾不尖灭是关键风险；峡谷切割深洼盆底扇与之类似，但侧封是最大风险（表2）。
表 2 琼东南盆地梅山组已钻盆底扇储层厚度统计Table 2. Statistics in the thickness that drilled into the basin-floor fan reservoir in the Meishan Formation of Qiongdongnan Basin浊积扇类型 各种类型储层厚度占比/% 大类 亚类 优质+好 有效 致密 盆
底
扇深洼盆底扇
峡谷切割深洼盆底扇85 5 10 断控缓坡根部盆底扇 76 5 19 断控陡坡根部盆底扇 10 87.5 2.5 斜坡扇 5 15 80 陆棚扇 0 13 87 此外，每期海底扇不同相带的储层物性亦存在较大差别[18-20]。A-7井（图7）梅山组海底扇中扇为均匀的中细砂岩，储层物性好，油迹清楚可见，以优质储层为主；内扇以中粗砂岩为主，夹顺层分布的片状泥砾，储层物性差，非均质性增强；外扇为砂泥互层，砂岩颗粒变细、泥质含量增加，储层物性有所降低。
2.4 底流改造
琼东南盆地陵水凹陷北部陆坡梅山组可见底流对重力流沉积物改造的标志[17]（图8）。当流速超过一定值时，底流可侵蚀和携带一定粒级的碎屑物质，尤其是细粒泥质沉积物，进而导致浊积砂体中泥质含量显著降低；使片状泥质碎屑顺层分布；清除孔隙内杂质；保持大孔喉半径并使喉道壁更加光滑。比较而言，底流对敞流环境沉积物，如海底扇中扇、外扇改造强，对限制环境（浊积沟道）的内扇改造弱[21-23]。
![]()
图 8 底流改造部分识别标志（A-1井）A. 3795.66m，压扁层理，泥岩分叉和颗粒定向排列；B. 3795.7m，生物潜穴、生物食遗被底流改造。Figure 8. Identification marks of turbidite fan reworked by bottom currents in the Well A-1A: 3795.66m，flattening bedding，mudstone bifurcation and grain orientation; B: 3795.7m，animal burrow and feeding tracks that reworked by bottom flow.表3显示了A-8井不同深度段壁芯样品的储层物性特征。如表所示，3个壁芯样品均为海底扇沉积物，岩性为粉砂岩，但因其所处相带不同，受底流改造程度不同，导致其储层物性存在较大差异。壁心1属内扇浊积沟道底部沉积，泥质含量高达25%，且含灰质，渗透率小于0.1×10−3 μm2，以细小孔喉为主，属致密层；壁心2属内扇上部沉积，泥质含量12%，渗透率0.711×10−3 μm2，存在一定的中大孔喉，但占比不高，属有效储层；壁心3为外扇沉积，泥质含量9%，渗透率3.72×10−3 μm2，大中孔喉占比很高，属于好储层（表3）。
表 3 A-8井梅山组海底扇储层实测物性特征Table 3. Measured physical property data of submarine fan reservoir of Well A-8序号 井深
/m沉积相 孔隙度
/%渗透率
/(10−3 μm2)泥质含量
/%描述 1 3 639 内扇浊积沟道 11.6 ﹤0.05 25 致密储层
小孔喉占比高2 3 615 内扇上部沉积 16.8 0.52 12 有效储层
中大孔喉普遍存在，但占比不高3 3 550.5 外扇沉积 18.2 3.33 9 好储层
中等孔喉占比很高综上所述，浊积扇储层以盆底扇为最有利，断控陡坡根部盆底扇、深洼盆底扇和峡谷切割深洼盆底扇等3种属于有利于储层发育的类型。海底扇相带中以中扇和外扇发育好储层、内扇物性相对较差。敞流环境底流改造强、限制性环境底流改造弱。底流活动强会显著降低泥质含量，制约碳酸盐的沉淀及其自生矿物的产生，有利于大孔隙的保存[4,14,17]。
3. 储层有效性分析
前已述及，琼东南盆地梅山组浊积扇发育好甚至优质储层。然而，好-优质储层并不等于成藏有效，只有圈闭条件好，与油气输导通道连通，才有可能成藏。探井钻前预测与实钻结果对比显示，仅18%的钻井基本符合；26%的钻井符合性较差；高达56%的钻井基本不符合。鉴于琼东南盆地裂后期总体处于伸展背景，构造圈闭不发育，岩性圈闭成为该区主要勘探目标。造成梅山组浊积扇储层成藏有效性成功率低的原因，除储层不发育外，沉积因素造成的岩性侧封不好、储层上倾不尖灭也是主要原因，前者主要见于峡谷切割型盆底扇、后主要见于深洼盆底扇。因此，两者构成了制约储层有效性的关键因素。
3.1 富泥峡谷充填侧封
琼东南盆地中央坳陷带发育了一条自西向东延伸的中央峡谷，该峡谷切割了梅山组甚至更深地层[18,24-25]。若梅山组富砂海底浊积扇被中央峡谷切割，峡谷充填过程及其内部组成、断控峡谷壁等对浊积扇储层成藏是否有效非常重要。中央峡谷的富泥充填形成的削截尖灭可以在浊积体上倾端或侧翼形成或协助形成圈闭。富泥沟道、沉积尖灭、局部断层共同控制了浊积扇油藏的形成和分布。
A-2井钻遇的梅山组海底扇被中央峡谷切割后与黄流组峡谷充填物形成侧向对接，峡谷下部充填为块体流沉积，空白弱反射，泥质为主，上部为平行连续强反射，属于黄流组II砂组（图9A）。富泥块体流沉积与梅山组海底扇T40A砂体侧向对接，形成了侧向封挡型圈闭，测试获得13层累计厚度18.7 m的气层。与A-2井截然不同的是，A-21井梅山组砂体被中央峡谷切割后，与峡谷充填的黄流组浊积砂体对接，侧向封挡失败（图9B），导致梅山组浊积扇储层含水。
![]()
图 9 梅山组峡谷切割海底扇过井地震剖面图Figure 9. Crossing well seismic sections of basin floor fan cut by canyon in the Meishan Formation3.2 储层上倾尖灭
钻探实践表明，琼东南盆地裂后期以发育岩性圈闭为主，且超过50%的岩性圈闭属于砂岩上倾尖灭型圈闭。梅山组海底浊积扇体因显著的地形高差和过路不沉积，虽然多数上倾尖灭风险不大，但若物源区富砂，碎屑颗粒粗，物源充足，则有可能存在储层上倾尖灭差的风险。
琼东南盆地乐东凹陷东北部梅山组B盆底扇紧邻物源，均方根振幅属性展示砂体的东南部发育很好的地层尖灭，但在NE、NW方向存在一定上倾尖灭风险（图10）。B-1井位于该圈闭 B 砂体较高部位，钻井揭示属于海底扇内扇，碎屑颗粒成分复杂、粒度粗，分选磨圆差、碳酸盐胶结严重。因中粗颗粒占比高，仍保留了不少原生孔隙、普遍发育好储层，渗透率在1×10−3 μm2以上，孔喉半径以中大者为主，且占一定比例（表4），测试结果为水层。造成储层含水的主要原因是扇体西侧浊积沟道含砂率高，保留了部分渗透性好的储层，导致B-1井所钻圈闭的砂体上倾尖灭不彻底，增加了圈闭成藏的风险。
![]()
图 10 B盆底扇均方根振幅地震属性图Figure 10. Seismic attribute of the root mean square amplitude of the basin floor fan B表 4 B-1井梅山组海底扇储层实测物性特征Table 4. Measured physical property data of submarine fan reservoir of Well B-1序号 井壁取芯 压汞资料 深度/m 岩性 照片 孔隙度/% 渗透率/(10−3μm2) 储层
类型孔隙度/% 渗透率
/(10−3μm2)压汞曲线及孔喉半径 1 3899 中砂岩 
6.7
6.1＜0.05
8.17致密 6 0.0438 
2 3950 中砂岩 
18.6
18.44.61
13好-优质 17.9 4.85 3 3968 粗砂岩 
16.1
15.62.08
4.85好-优质 15.7 2.46 4 3980.5 中砂岩 
14.7 0.88 好-优质 15.2 2.79 
5 3994 含砾
粗砂岩
13.6 1.89 好-优质 12.8 2.19 6 4002 中砂岩 
16.1 1.4 好-优质 15.7 1.77 7 4017 粗砂岩 
14.7 0.5 有效 14.7 0.75 4. 结论与建议
（1）琼东南盆地梅山组浊积扇埋深有较大差异，其浊积扇储层物性除受到埋深影响外，还受到物源、海平面变化、扇体类型、内部构成以及底流改造的影响。浊积扇类型及其内部组成、底流改造等是影响和控制浊积体储层品质的主要因素。低位时期的断控陡坡根部盆底扇、深洼盆底扇和峡谷切割深洼盆底扇是好-优质储层相对发育的浊积体类型；从相带来看，中扇和外扇发育好储层，内扇储层非均质性强；此外，在底流强烈活动区，底流改造带走泥质、改善孔隙内壁，明显提高了储层物性。
（2）梅山组浊积扇储层成藏有效性主要取决于岩性侧封和储层上倾尖灭两方面。切割浊积扇的富泥峡谷充填和块体流沉积可以形成有效的顶部和侧向封堵，岩性及断裂侧向封堵是海底扇岩性圈闭形成的关键，特别是精细的断层分析与成图是浊积扇储层成藏有效性评价的关注的重点。
-
![]()
图 1 现代滦河三角洲位置及表层沉积物取样站位
Figure 1. Locations of the LRSD and surface sediment sampling stations in the study area
![]()
图 2 滦河年内径流率、输沙率（a）与年际径流量、输沙量（b）
Figure 2. Monthly water discharge rate and sediment load rate （a） and annual water discharge and sediment load （b） of the Luanhe River
![]()
图 3 滦河三角洲表层沉积物粒径平面分布
Figure 3. Spatial distribution of grain size parameters of the surface sediment samples collected from the LRSD
![]()
图 4 原始（a、b）、规则（c、d）沉积物站位分布与半方差函数
Figure 4. Sampling locations and semi-variances for the irregular （a, b） and regular grids （c, d）
![]()
图 5 滦河水下三角洲粒径净输运趋势
a. 插值后CB+模式，特征距离0.09°；b. 插值后FB-模式，特征距离0.09°；c. 未插值FB-模式，特征距离0.042°；d. 未插值FB-模式，特征距离0.059°。
Figure 5. Sediment transport trend in the modern LRSD
a. CB+ case of regular data with characteristic distance of 0.09 decimal degree; b. FB- case of regular data with characteristic distance of 0.09 decimal degree; c.FB- case of irregular data with characteristic distance of 0.042 decimal degree; d.FB- case of irregular data with characteristic distance of 0.059 decimal degree.
![]()
图 8 滦河水下三角洲沉积物输运模式
红色实线为枯季余流，红色虚线表示枯季输沙，绿色实线为洪季余流，绿色虚线为洪季输沙。
Figure 8. Model of the sediment transportation in the LRSD
Red solid lines denote the residual current in dry season, red dotted lines denote the sediment transport in dry season, green solid lines denote the residual current in flood season, and green dotted lines denote the sediment transport in flood season.
![]()
图 9 水深-离岸距离与平均粒径-离岸距离图 （位置见图1）
Figure 9. Plots of depth versus offshore distance and mean grain size versus offshore distance （see fig.1 for locations）
表 1 粒径参数的Wilcoxon符号秩检验（α=0.99）
Table 1 Wilcoxon non-parameter test for the grain size parameters（α=0.99）
插值半径/（°） 粒径参数 P 值 h 0.016 平均粒径 1 0 分选系数 1 0 偏态 1 0 0.015 平均粒径 1 0 分选系数 1 0 偏态 1 0 0.014 平均粒径 0.801 1 分选系数 0.801 1 偏态 0.801 1 
下载: 导出CSV
-
[1] Yamashita S, Naruse H, Nakajo T. Reconstruction of sediment-transport pathways on a modern microtidal coast by a new grain-size trend analysis method [J]. Progress in Earth and Planetary Science, 2018, 5: 7. doi: 10.1186/s40645-018-0166-9
[2] McCave I N. Grain-size trends and transport along beaches: example from eastern England [J]. Marine Geology, 1978, 28(1-2): M43-M51. doi: 10.1016/0025-3227(78)90092-0
[3] McLaren P. An interpretation of trends in grain size measures [J]. Journal of Sedimentary Research, 1981, 51(2): 611-624.
[4] McClaren P, Ren P. Sediment transport and its environmental implications in the lower Fraser River and Fraser delta[M]. Environment Canada: Environmental Conservation, 1995.
[5] Gao S, Collins M. Net sediment transport patterns inferred from grain-size trends, based upon definition of “transport vectors” [J]. Sedimentary Geology, 1992, 81(1-2): 47-60. doi: 10.1016/0037-0738(92)90055-V
[6] Ma F, Wang Y P, Li Y, et al. The application of geostatistics in grain size trend analysis: A case study of eastern Beibu Gulf [J]. Journal of Geographical Sciences, 2010, 20(1): 77-90. doi: 10.1007/s11442-010-0077-1
[7] Sánchez A, Choumiline E, López-Ortiz B E, et al. Sediment transport patterns in Magdalena Bay, Baja California Sur, Mexico, inferred from grain-size trends [J]. Latin American Journal of Aquatic Research, 2010, 38(2): 167-177. doi: 10.3856/vol38-issue2-fulltext-1
[8] Choi T J, Choi J Y, Park J Y, et al. The effects of nourishments using the grain-size trend analysis on the intertidal zone at a sandy macrotidal beach [J]. Journal of Coastal Research, 2018, 85(S1): 426-430.
[9] De Falco G, Molinaroli E, Baroli M, et al. Grain size and compositional trends of sediments from <italic>Posidonia oceanica</italic> meadows to beach shore, Sardinia, western Mediterranean [J]. Estuarine, Coastal and Shelf Science, 2003, 58(2): 299-309. doi: 10.1016/S0272-7714(03)00082-9
[10] Nugroho S H, Putra P S. Spatial distribution of grain size and depositional process in tidal area along Waikelo Beach, Sumba [J]. Marine Georesources & Geotechnology, 2017, 36(3): 299-307.
[11] Pedreros R, Howa H, Michel D. Application of grain size trend analysis for the determination of sediment transport pathways in intertidal areas [J]. Marine geology, 1996, 135(1-4): 35-49. doi: 10.1016/S0025-3227(96)00042-4
[12] Clarke D W, Boyle J F, Chiverrell R C, et al. A sediment record of barrier estuary behaviour at the mesoscale: interpreting high-resolution particle size analysis [J]. Geomorphology, 2014, 221: 51-68. doi: 10.1016/j.geomorph.2014.05.029
[13] Hegde V S, Nayak S, Krishnaprasad P, et al. Evolution of diverging spits across the tropical river mouths, central west coast of India [J]. Journal of Coastal Zone Management, 2015, 18(2): 1000402.
[14] Li T, Li T J. Sediment transport processes in the Pearl River Estuary as revealed by grain-size end-member modeling and sediment trend analysis [J]. Geo-Marine Letters, 2018, 38(2): 167-178. doi: 10.1007/s00367-017-0518-2
[15] Cheng P, Gao S, Bokuniewicz H. Net sediment transport patterns over the Bohai Strait based on grain size trend analysis [J]. Estuarine, Coastal and Shelf Science, 2004, 60(2): 203-212. doi: 10.1016/j.ecss.2003.12.009
[16] 程鹏, 高抒. 北黄海西部海底沉积物的粒度特征和净输运趋势[J]. 海洋与湖沼, 2000, 31(6):604-615. [CHENG Peng, GAO Shu. Net sediment transport patterns over the northwestern Yellow Sea, based upon grain size trend analysis [J]. Oceanologia et Limnologia Sinica, 2000, 31(6): 604-615. doi: 10.3321/j.issn:0029-814X.2000.06.004 [17] Poizot E, Méar Y, Biscara L. Sediment Trend Analysis through the variation of granulometric parameters: A review of theories and applications [J]. Earth Science Reviews, 2008, 86(1-4): 15-41. doi: 10.1016/j.earscirev.2007.07.004
[18] Gao S, Collins M, Mclaren P, et al. A critique of the "McLaren Method" for defining sediment transport paths; discussion and reply [J]. Journal of Sedimentary Research, 1991, 61(1): 143-147. doi: 10.1306/D42676A9-2B26-11D7-8648000102C1865D
[19] Gao S, Xie Q C, Feng Y J. Fine-grained sediment transport and sorting by tidal exchange in Xiangshan Bay, Zhejiang, China [J]. Estuarine, Coastal and Shelf Science, 1990, 31(4): 397-409. doi: 10.1016/0272-7714(90)90034-O
[20] Gao S, Collins M. Tidal Inlet equilibrium, in relation to cross-sectional area and sediment transport patterns [J]. Estuarine, Coastal and Shelf Science, 1994, 38(2): 157-172. doi: 10.1006/ecss.1994.1010
[21] Poizot E, Mear Y, Thomas M, et al. The application of geostatistics in defining the characteristic distance for grain size trend analysis [J]. Computers & Geosciences, 2006, 32(3): 360-370.
[22] 高抒. 沉积物粒径趋势分析: 原理与应用条件[J]. 沉积学报, 2009, 27(5):826-836. [GAO Shu. Grain size trend analysis: principle and applicability [J]. Acta Sedimentologica Sinica, 2009, 27(5): 826-836. [23] 高善明, 李元芳, 安凤桐, 等. 滦河三角洲滨岸沙体的形成和海岸线变迁[J]. 海洋学报, 1980, 2(4):102-114. [GAO Shanming, LI Yuanfang, AN Fengtong, et al. The formation of sand bars on the Luanhe River Delta and the change of the coast line [J]. Acta Oceanologica Sinica, 1980, 2(4): 102-114. [24] 李从先, 陈刚, 王传广, 等. 论滦河冲积扇-三角洲沉积体系[J]. 石油学报, 1984, 5(4):27-36. [LI Congxian, CHEN Gang, WANG Guangchuan, et al. On the Luanhe River Alluvial Fan-delta complex [J]. Acta Petrolei Sinica, 1984, 5(4): 27-36. doi: 10.7623/syxb198404004 [25] 薛春汀. 滦河冲积扇-三角洲的范围和类型及其演化[J]. 海洋地质与第四纪地质, 2016, 36(6):13-22. [XUE Chunting. Extents, type and evolution of Luanhe River Fan-delta system, China [J]. Marine Geology & Quaternary Geology, 2016, 36(6): 13-22. [26] 王保民. 河北省海洋资源调查与评价综合报告[M]. 北京: 海洋出版社, 2007. WANG Baomin. Comprehensive Report on Investigation and Evaluation of Marine Resources in Hebei Province[M]. Beijing: Ocean Presee, 2007.
[27] 王保民. 河北省海洋资源调查与评价专题报告(上册、下册)[M]. 北京: 海洋出版社, 2007. [ WANG Baomin. Thematic Reports (Volume I, Volume II) on Investigation and Evaluation of Marine Resources in Hebei Province[M]. Beijing: Ocean Presee, 2007. ]
[28] 康雪宁, 印萍, 刘金庆. 我国中小河流入海水沙变化对人类活动响应——以滦河为例[J]. 海洋地质与第四纪地质, 2016, 36(6):1-6. [KANG Xuening, YIN Ping, LIU Jinqing. Variations in water and sediment discharges of mediam and small rivers and their response to human activities: a case study on the Luan River [J]. Marine Geology & Quaternary Geology, 2016, 36(6): 1-6. [29] 黎刚, 殷勇. 滦河下游河道及三角洲地貌的近期演化[J]. 地理研究, 2010, 29(9):1606-1615. [LI Gang, YIN Yong. Recent geomorphological evolution of downstream channel and delta of Luanhe River [J]. Geographical Research, 2010, 29(9): 1606-1615. [30] Mcmanus J. Grain size determination and interpretation[M]//Tucker M. Techniques in Sedimentology. Oxford, UK: Blackwell Scientific Publications, 1988: 63-85.
[31] Poizot E, Méar Y. eCSedtrend: A new software to improve sediment trend analysis [J]. Computers & Geosciences, 2008, 34(7): 827-837.
[32] Yu X X, Li T G, Gu D Q, et al. Sediment transport in the Luanhe River delta: grain size trend analysis [J]. Journal of Oceanology and Limnology, 2019, 37(3): 982-997. doi: 10.1007/s00343-019-8156-3
[33] Ashley G M. Interpretation of Polymodal Sediments [J]. The Journal of Geology, 1978, 86(4): 411-421. doi: 10.1086/649710
[34] 孙东怀, 安芷生, 苏瑞侠, 等. 古环境中沉积物粒度组分分离的数学方法及其应用[J]. 自然科学进展, 2001, 11(3):269-276. [SUN Donghuai, AN Zhisheng, SU Ruixia, et al. Method and application of numerical partitioning of the sediment components in paleoclimate environments [J]. Progress in Natural Science, 2001, 11(3): 269-276. doi: 10.3321/j.issn:1002-008X.2001.03.008 [35] 孙东怀, 鹿化煜, Rea D, et al. 中国黄土粒度的双峰分布及其古气候意义[J]. 沉积学报, 2000, 18(3):327-335. [SUN Donghuai, LU Huayu, Rea D, et al. Bimode Grain-size distribution of Chinese loess and its paleoclimate implication [J]. Acta Sedimentologica Sinica, 2000, 18(3): 327-335. doi: 10.3969/j.issn.1000-0550.2000.03.001 [36] 段晓勇, 印萍, 刘金庆, 等. 滦河口表层沉积物中重金属和多环芳烃的分布、来源及风险评估[J]. 中国环境科学, 2016, 36(4):1198-1206. [DUAN Xiaoyong, YIN Ping, LIU Jinqing, et al. Heavy metals and polycyclic aromatic hydrocarbons in surface sediments of Luan River estuary: distributions, sources and ecological risk assessments [J]. China Environmental Science, 2016, 36(4): 1198-1206. doi: 10.3969/j.issn.1000-6923.2016.04.036 [37] 刘金庆, 印萍, 张勇, 等. 滦河口沉积物重金属分布及生态风险评价[J]. 海洋地质与第四纪地质, 2016, 36(5):43-52. [LIU Jinqing, YIN Ping, ZHANG Yong, et al. Distribution of heavy metals in surface sediments of the Luanhe River Estuary and ecological risk assessment [J]. Marine Geology & Quaternary Geology, 2016, 36(5): 43-52. [38] 乔淑卿, 石学法, 王国庆, 等. 渤海底质沉积物粒度特征及输运趋势探讨[J]. 海洋学报, 2010, 32(4):139-147. [QIAO Shuqing, SHI Xuefa, WANG Guoqing, et al. Discussion on grain-size characteristics of seafloor sediment and transport pattern in the Bohai Sea [J]. Acta Oceanologica Sinica, 2010, 32(4): 139-147. [39] Xue Z, Feng A P, Yin P, et al. Coastal erosion induced by human activities: a northwest Bohai Sea case study [J]. Journal of Coastal Research, 2009, 25(3): 723-733.
[40] 刘振夏. 现代滦河三角洲的影响因素和沉积物分区[J]. 海洋科学进展, 1989, 7(4):55-64. [LIU Zhenxia. The influential factors and zoning of the modern Luanhe River Delta [J]. Journal of Oceanography of Huanghai & Bohai Seas, 1989, 7(4): 55-64. [41] 姜太良, 房宪英, 徐洪达, 等. 滦河口径流量和滦河口输沙分析[J]. 海洋科学进展, 1986, 4(4):100-113. [JIANG Tailiang, FANG Xianying, XU Hongda, et al. Analysis of the water discharge and sediment load in the Luanhe River Estuary [J]. Journal of Oceanography of Huanghai & Bohai Seas, 1986, 4(4): 100-113. -
期刊类型引用(2)
1. 冯永财，郝连成，胡延斌，褚宏宪，李佳林，远继东，黄杏，陈晓日，王超. 渤海湾曹妃甸海域沉积物粒度分布特征及运移趋势分析. 海洋地质前沿. 2024(11): 57-69 . 
百度学术
2. 徐利，郝桂珍，李思敏，段雯瑜，张亚婧. 滦河上游流域水质特征及污染源解析. 科学技术与工程. 2023(16): 7136-7144 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 2413
- HTML全文浏览量: 553
- PDF下载量: 43
- 被引次数: 4
