热带海洋学报

• • 上一篇    下一篇

基于Bekker模型的砂质海滩车辆陷车风险蒙特卡洛评估

李帅磊1,2,3,唐诚1,3,郑向阳1,3,李东1,李艳芳1,3,赵明亮1   

  1. 1. 中国科学院烟台海岸带研究所,山东 烟台 264003;

    2. 中国科学院大学,北京 100049;

    3. 山东省海岸带环境过程与生态安全重点实验室,山东 烟台 264003



  • 收稿日期:2026-01-16 修回日期:2026-03-30 接受日期:2026-04-08
  • 通讯作者: 唐诚
  • 基金资助:
    国家自然科学基金(41006055)

Monte Carlo Assessment of Vehicle Sinking Risk on Sandy Beaches Based on the Bekker Model

LI Shuailei1,2,3,TANG Cheng1,3,ZHENG Xiangyang1,3,LI Dong1,LI Yanfang1,3,ZHAO Mingliang1   


  • Received:2026-01-16 Revised:2026-03-30 Accepted:2026-04-08
  • Supported by:

    National Natural Science Foundation of China(41006055)

摘要: 车辆在滨海作业和两栖登陆时频发陷车问题,其根源在于砂质滩面力学性质受潮汐驱动而表现出的时变特性。由于地下水系统渗流对外部潮汐的响应存在相位滞后,浅层砂体在涨潮中后期容易出现短时间内的高孔隙水压力与低有效应力状态,造成抗剪强度和地基承载能力骤减。为量化水力耦合引发的承载弱化效应,本文构建了基于“时变Bekker沉陷模型”和Monte Carlo评估框架的车辆陷车风险评估框架。通过把有效应力的时态演化映射成地基承载参数的非线性退化,并引入渗透扩散和孔压滞后等水文动力学参数的不确定性。计算结果表明,陷车概率随潮位的上涨而明显增加,但风险峰值通常滞后于潮位峰值约0.5到1.0小时,揭示滩面失稳受“孔压累积-承载退化-沉陷突增”的滞后跃迁机制控制。此外,敏感性分析结果表明,地基最小承载参数和临界沉陷判据决定着陷车风险水平;由于在承载下降和液化软化之间多物理参数的共同作用下,陷车概率的放大效应可达单因素独立演化的1.65倍。该模型评估从机理层面解释了潮间带车辆失稳的迟滞现象,可为潮滩车辆通行安全预警、两栖登陆窗口期选择及滨海作业风险评价提供定量参考。

关键词: 潮-地下水相互作用, 砂质海滩, 车辆陷车风险, Bekker沉陷模型, 蒙特卡洛模拟

Abstract: Vehicles frequently encounter immobilization problems during coastal operations and amphibious landings, which is fundamentally rooted in the time-varying characteristics of the mechanical properties of sandy beaches driven by tidal forcing. Due to the phase-lagged response of groundwater seepage to external tides, shallow sand bodies are highly susceptible to transient states of high pore water pressure and low effective stress during the mid-to-late stages of the rising tide, leading to a sharp decline in shear strength and ground bearing capacity. To quantify the bearing weakening effect induced by hydro-mechanical coupling, this paper constructs a vehicle immobilization risk assessment framework based on a "time-varying Bekker sinkage model" and Monte Carlo simulations. This framework maps the temporal evolution of effective stress onto the non-linear degradation of ground bearing parameters, incorporating the uncertainties of hydrodynamic parameters such as hydraulic diffusivity and pore pressure lag. Computational results demonstrate that the immobilization probability increases significantly with the rising tide; however, the risk peak typically lags behind the tidal peak by approximately 0.5 to 1.0 hours. This reveals that beach surface instability is governed by a delayed transition mechanism characterized by "pore pressure accumulation-bearing degradation-sinkage surge". Furthermore, sensitivity analysis indicates that the minimum ground bearing parameter and the critical sinkage criterion determine the overall level of immobilization risk. Under the synergistic action of multiple physical parameters between bearing capacity reduction and liquefaction softening, the amplification effect on the immobilization probability can reach 1.65 times that of a single-factor independent evolution. This assessment model mechanistically explains the hysteresis phenomenon of vehicle instability in the intertidal zone, providing a quantitative reference for the early safety warning of tidal-flat vehicle traffic, the selection of amphibious landing windows, and the risk assessment of coastal operations.

Key words: Tide-groundwater interaction, Sandy beach, Vehicle immobilization risk, Bekker pressure-sinkage model, Monte Carlo simulation