黄河三角洲油井废弃后滨海湿地生态恢复轨迹

《Frontiers in Plant Science》：Ecological recovery trajectories of coastal wetlands following oil-well abandonment in the Yellow River Delta

　　

摘要 引言：油井退役虽终止生产，但可能在滨海湿地留下物理、化学和生态遗留效应。 方法：研究人员于2021年在中国黄河三角洲评估了不同退役阶段下的土壤条件、重金属风险及植被状况，共涉及15个油田相关样地和6个自然参考湿地。样地被分为活跃期（Active）、废弃（Abandoned）、修复（Restored，即拆除基础设施并实施生态修复）和参考（Reference）四类，退役样地按<3年、3–5年和5–10年三个年龄等级划分。 结果：油田样地的平均植被盖度低于参考湿地（6月：50.29% vs. 69.00%；12月：38.00% vs. 56.23%），同时pH值较高而盐度普遍较低。镉（Cd）占潜在生态风险的54.3%–65.7%，其单因子污染指数在三个样地超过1。修复样地的平均盖度和物种丰富度随年龄等级增加而上升（从<3年的15.8%和1.5种增至5–10年的64.2%和7.0种），而多样性在中龄等级达到峰值。Cd指数、污染负荷指数（PLI）和潜在生态风险指数（RI）并未随退役年龄单调递减。在5–10年组中，修复样地的平均金属风险指数低于废弃样地，但在3–5年组中则更高。土壤pH值与盖度和丰富度呈正相关；Cd指数和PLI与植被指标无显著关联。 讨论：这些描述性时序序列模式表明，植被重组与化学风险衰减并不遵循共同的轨迹，且不能仅凭退役年龄推断恢复情况。湿地退役管理应将样地特异性污染评估与植被和栖息地评价相结合。

论文解读文章

一、研究背景与科学问题

滨海湿地处于海陆动态交界地带，提供生物多样性维持、碳储存、海岸防护和水质净化等至关重要的生态系统服务。许多石油产区在地理空间上与三角洲、河口和滨海平原重叠，钻井平台、道路、堤坝、排水渠及意外泄漏事件会改变维系湿地功能的地貌、水文、土壤和生物过程，并在生产活动停止后以生态遗留效应的形式持续存在。以往关于油田湿地的研究主要聚焦于生产活跃期的生态影响，包括污染、植被退化、栖息地破碎化及水文连通性降低。然而，生产停止并不能自动消除场地、道路、站场或堤坝的物理足迹——土壤压实、微地形改变、水文阻隔及残留污染物可能继续制约植物的定植与演替。

生态恢复并非单一指标的定向变化，而是生态系统各组分的重组过程，并可能向参照状态趋近。恢复路径因胁迫类型和沿海生态系统的物理环境而异。重金属作为化学维度的重要评估对象，其迁移性和有效性取决于吸附、沉积物氧化还原条件、盐度及有机质相互作用，因此区域平均浓度可能掩盖局部富集，而较低的污染指数本身也不足以证明植物群落已经恢复。此外，退役后的管理方式（单纯废弃vs.拆除设施并进行生态修复）可能进一步影响恢复路径。针对上述问题，研究人员以黄河三角洲这一年轻河口湿地复合体为研究对象，利用空间时序序列设计，系统比较了活跃期、废弃、修复及参考湿地四类样地在不同退役年龄阶段下土壤条件、重金属风险指标与植被属性的差异，旨在厘清生产停止与生态恢复之间的非等价关系。

二、关键技术方法

研究人员采用空间时序序列（chronosequence）替代时间追踪观测，以不同退役年龄的样地代表恢复进程的不同阶段。监测网络涵盖黄河三角洲（山东省）15个油田样地和6个天然参考湿地，样地经纬度范围为118.7845°–119.2208°E、37.7248°–37.8716°N。土壤样品分别于2021年6月和12月采集，经风干、均质、过筛处理后，测定pH值和盐度（1:5土水比提取法），并使用原子吸收光谱法和原子荧光光谱法测定砷（As）、镉（Cd）、总铬（Cr）、铜（Cu）、铅（Pb）和汞（Hg）含量。基于上述数据计算单因子污染指数、污染负荷指数（PLI）和潜在生态风险指数（RI）。植被调查在每个油田样地设置3个样方，记录盖度、物种数、密度、高度和频度，进而计算物种丰富度、Gini-Simpson多样性和Shannon多样性。数据分析采用Pearson相关系数考察金属指数间的协变关系，采用Spearman秩相关系数检验土壤条件、污染指数与植被指标之间的关联。

三、研究结果

3.1 油田样地与参考湿地在土壤条件和植被结构上存在差异

油田样地的平均pH值高于参考湿地（9.28±0.31 vs. 8.03±0.05），平均盐度更低（0.46±0.53 vs. 1.90±0.52 ppt）。两次调查中油田样地的植被盖度均低于参考湿地：6月差值18.71个百分点，12月差值18.23个百分点。多样性对比方向随季节反转：6月油田样地的平均丰富度（5.07种）、Gini-Simpson多样性（0.411）和Shannon多样性（0.819）均高于参考湿地（分别为3.33种、0.315和0.633）；12月则相反，油田样地各项指标均低于参考湿地。

3.2 Cd主导样地级重金属风险信号

在所有油田样地中，Cd的单因子污染指数平均值为0.92，其中O02、O03和O10三个样地超过1（分别为1.37、1.21和1.40），其余四种元素均未超过1。PLI均值1.02，六个样地高于1（含全部五个活跃期样地）。RI均值52.58，最大值出现在O10（76.49）。Cd贡献占各样地总RI的54.3%–65.7%，平均贡献率达59.4%。所有金属指数间均为正相关，Pearson系数介于0.658至0.972之间。

3.3 金属风险指数在不同管理状态和退役年龄等级间的分布

活跃期样地的Cd指数、PLI和RI均值最高（分别为1.10、1.18和61.70）。修复样地中，三项指数均在<3年组最低，3–5年组升高，5–10年组略有回落，未呈现单调递减趋势，且无任何修复样地Cd指数超过1。在3–5年组中，废弃样地的三项指数均值均低于修复样地；在5–10年组中则相反，修复样地均值更低。两个较老废弃样地之间存在较大差异：O03的Cd指数为1.21、RI为68.09，而O07仅为0.76和44.67。

3.4 植被指标在不同管理状态和退役年龄等级间的变化

修复样地的6月平均盖度从<3年组的15.8%升至3–5年组的59.2%和5–10年组的64.2%，最老组与参考湿地均值（69.00%）的差距仅4.8个百分点。平均丰富度同样递增（1.5→5.5→7.0种），最老组每个样地6–8种，超出参考湿地范围（2–5种）。但Gini-Simpson和Shannon多样性均在3–5年组达到峰值后下降。废弃样地方面，3–5年组盖度均值（63.3%）略高于同期修复样地（59.2%），但多样性较低；5–10年组则相反。活跃期样地盖度范围跨度大（28.33%–80.00%）。此外，全油田样地合并物种组成显示，6月记录23种，12月仅6种；芦苇（Phragmites australis）在两次调查中分别出现于12/15和13/15个样地，其密度占比从28.4%升至87.4%。

3.5 土壤pH值与植被盖度和丰富度呈正相关

在14个完整配对数据的油田样地中，pH值与6月盖度（Spearman r=0.548, P=0.043）和丰富度（r=0.561, P=0.037）呈显著正相关。盐度与四项植被指标均为负相关（r介于-0.477至-0.406之间），但均未达显著水平。Cd指数和PLI与植被指标无显著相关性（P值均>0.05）。

四、讨论与结论总结

讨论部分指出，生产停止不应被视为湿地恢复的代用指标：所有退役组的金属风险均值均低于活跃期组，但退役组内部并未形成一致的风险递降序列；植被恢复则表现为组分特异性模式，而非简单的“退化—恢复”单向排序。Cd是潜在生态风险的主要贡献者，其主导地位既源于相对本地参比的浓度水平，也因其在Hakanson框架中的高毒性响应系数。修复样地中最老组的盖度接近参考均值，但丰富度超出参考范围，多样性则在中龄组达到峰值后回落，说明植被恢复涉及物种更替与优势度变化，而非单调趋向参照状态。土壤pH与植被的正相关反映的是扰动景观内的协变关系，而非提高pH即可恢复参照条件的因果关系。研究人员强调，浓度基指数不直接代表金属生物有效性，后者受沉积物化学和水文条件调控。

结论部分翻译如下：黄河三角洲油井退役伴随的是金属风险与植被模式的非协同变化。退役组的平均金属风险指数低于活跃期组，但这些指数并未随退役年龄单调递减；Cd始终是潜在生态风险的主导贡献者。较老修复样地的盖度更接近参考均值，而丰富度和多样性呈现不同的变化模式。废弃与修复样地之间的比较也因年龄和指标而异。因此，恢复评估应将样地特异性污染评价与植被和栖息地评估相结合，仅凭时间推移或基础设施拆除不足以证明生态恢复已经实现。