河口模型中地貌演变与潮汐能量耗散研究

《Oceans》：On Morphological Change and Tidal Energy Dissipation in a Modelled Estuary

　　

随着潮汐河口的地形（bathymetry）发生变化，随着河道、曲流与浅滩发生调整，整体能量耗散率也会改变。这一过程符合熵产增加原理，因为不断演化的海床特征增加了系统的无序度。在一个收敛型潮汐河口中，研究人员对250年地形演变进行了模型模拟，结果显示能量耗散随时间降低，在演化河道中变得局部化，同时又从初始以潮间带河段为主的状态向下游重新分布。随着地形趋向平衡状态，地貌活动减少，其特征是平均河道宽度与平均能量耗散率沿河口向上均近似线性下降。有无形态加速（morphological acceleration）的模拟都显示出能量耗散降低的相同总体行为，但在采用形态因子（morphological factor, MF）时，耗散下降更慢。这部分原因似乎在于：未采用MF时，内河口预测的浅滩与河道纹理更细，因而比采用MF时的预测结果具有更低的耗散率。该研究确认，随着河口向平衡演化，熵产下降。

研究背景方面，河口地貌演化兼具人类活动与生态环境影响，既关系战略区位的经济价值，也支撑丰富生态系统。理解该演化过程需要兼顾泥沙输运过程与更广泛物理原理。早期研究多聚焦河流系统，提出最小_stream功率与最小功原理；潮汐河口因沿程流量可变，是比河流更一般的情形，已有研究将热力学类比如熵产与最小功条件用于分析地貌变化与能量耗散的关系。既有长期过程模型常借助形态因子（MF，morphological factor，即在每个水动力时间步放大床面高程变化以缩短计算时间的缩放系数）加速模拟，但少有研究细致展示理想化河口中潮汐能量耗散随地貌演化的基础过程，且MF对细尺度河道/浅滩纹理及耗散率的影响尚缺乏详细比对。

研究人员以收敛型理想河口为对象，用TELEMAC系统开展250年二维条形（barotropic）水动力与泥沙模块SISYPHE双向耦合模拟，基线取MF=1，即不加速；对比MF=10、20、40。论文发表于《Oceans》。研究目的在于确定250年内能量耗散下降速率、宽平均地形剖面变化、MF有无是否得到相似河道/浅滩与耗散结果，并从河道宽度所刻画的细纹理解释差异。结论上，河口向平衡演化时形态活动减弱，域平均能量耗散率下降，熵产降低；MF使耗散下降更慢，但因跳过中间过程而损失内河口细纹理，MF1预测更细河道/浅滩体系。

主要关键技术方法方面，研究人员采用TELEMAC水文信息系统，以telemac2d进行二维条形水动力计算，以SISYPHE进行泥沙输运与地貌演化，二者双向耦合并反馈时变床面粗糙度k_s；边界为口外开边界给定M₂分潮、侧边界与上端滑移边界、口外床面固定；利用超级计算避免MF以MF1跑满250年，并以MF10/20/40作加速对照；床沙取单粒径0.2 mm砂，输沙用Bijker公式；河道宽度通过规则网格上掩膜浅滩后最大内切圆拟合估算；宽平均地形用4 km盒式滤波平滑；能量耗散率按床剪应力做功τ_bU时间平均与空间平均获得。

3.1. 初始流与泥沙输运

研究人员由第3天结果识别初始潮间带河段，显示口门速度洪枯不对称弱，向上游X=20 km处速度增大且洪潮不对称显著，X=28 km处出现短时强洪流与长时弱枯流；泥沙输运在X=20 km后洪向占优，潮间带以上以向上游脉冲为主。结论为初始不对称驱动后续泥沙再分布与地貌成形。

3.2. 河道与浅滩发育

由10年、50年、240年床面演化可见，无随机扰动下非结构网格也能触发图纹；河道/浅滩自潮间带向上、向下扩展，主曲流在固定边界处加深。240年时MF1显示特定上边界曲流，MF20则出现不同下边界曲流位置，MF10/40图纹亦异。研究人员结论为MF改变演化事件序列，使河道空间排布与细纹理不同，但各情形均见河道宽度向上游变小。

3.3. 宽平均地形剖面发育

宽平均床面演化显示，前40年上游净淤、中部净冲，后随河道形成把物质送向下游；MF情形减少上游淤积、增强近口沉积。研究人员结论为宽度平均剖面受洪枯不对称与河道化共同控制，MF不给出单调剖面差异。

3.4. 能量耗散

耗散率∝床应力×流速，因k_s时空变化可高于三次律。初始域平均0.347 W/m²，50年降至0.257 W/m²，耗散由均匀转为河道内局部化、浅滩上很低。250年内MF1降至初始52%，MF10/20/40分别为61%、62%、60%；沿程峰值由X=20 km下移，后期5–25 km近似线性下降。研究人员结论为熵产随演化下降，MF使下降更慢，耗散向下游重新分布。

3.5. 河道-浅滩点数比

以宽平均剖面为基准划分河道/浅滩点，MF1在约110年后趋稳，浅滩/河道点数比240年达1.31；MF情形比值更高，通道点更多。研究人员结论为系统向后形态活动更低状态演化，浅滩占比上升但MF略偏河道化。

3.6. 河道宽度

用最大内切圆得到河道宽度，活跃区随时间长大；最大宽度沿上游减小且不受MF明显影响，但平均宽度在MF1上游更小，说明MF1保留更细河道/浅滩纹理。研究人员结论为更细纹理使MF1上游耗散更低，耗散与河道宽均向上游近似线性下降可作平衡态表征。

讨论与结论翻译如下：理想化收敛潮汐河口250年地形演化模拟聚焦能量耗散随时变。已知河口趋向低形态活动与低能耗平衡，但细节未被突出。本文详述河道/浅滩演化与耗散下降关系：从平面可冲床面出发，图纹始于潮间带，向口门与河头扩展，洪枯不对称因河口收敛而增强并放大净冲淤。耗散渐聚于河道、浅滩耗散低，全域耗散下降，证实熵产下降、平衡趋近。本研究借超算免用MF，并以MF10/20/40对照；各MF图纹互异，但均现耗散渐进降低，MF使降速更缓。河道宽度为诊断量，所有情形向上游变窄，活跃区扩大；240年后耗散与河宽在活跃区向上游近似线性降。MF1上游耗散略小，因其浅滩/河道体系比MF10/20/40更细。尽管有差，研究显示用MF亦可较简满意预测系统总体行为。