[发明专利]一种淤泥质海岸浮泥形成条件的定量计算方法

权利要求书

1.一种淤泥质海岸浮泥形成条件的定量计算方法，其特征在于，包括如下

步骤：

步骤A，根据港口工程的位置，建立工程海域的二维潮流泥沙数学模型，采用有限体积法求解二维潮流泥沙数学模型，得到浮泥形成条件定量计算所需的水动力及泥沙参数，所述二维潮流泥沙数学模型的控制方程见式(1)～(6)：

①连续方程：

∂h∂t+∂hu‾∂x+∂hv‾∂y=hS---(1)]]>

②运动方程：

∂hu‾∂t+∂hu‾2∂x+∂hu‾v‾∂y=fv‾h-gh∂η∂x-hρ0∂pa∂x-]]>

                                   (2)

gh22ρ0∂ρ∂x+τsxρ0-τbxρ0+∂∂x(hTxx)+∂∂y(hTxy)+husS]]>

∂hv‾∂t+∂hu‾v‾∂x+∂hv‾2∂y=-fu‾h-gh∂η∂y-hρ0∂pa∂y-]]>

                                    (3)

gh22ρ0∂ρ∂y+τsyρ0-τbyρ0+∂∂x(hTxy)+∂∂y(hTyy)+hvsS]]>

Txx=2A∂u‾∂x,]]>Txy=A(∂u‾∂y+∂v‾∂x),]]>Tyy=2A∂v‾∂y---(4)]]>

A=cs2l22SijSij,]]>Sij=12(∂ui∂xj+∂uj∂xi)]]>(i，j＝1，2)(5)

③悬沙运动方程：

∂c‾∂t+u∂c‾∂x+v∂c‾∂y=]]>

                                (6)

1h∂∂x(hDx∂c‾∂x)+1h∂∂y(hDy∂c‾∂y)+QLCL1h+E]]>

上式中：x、y-直角坐标系坐标；t-时间变量；η-波面高程；d-静水水深；h＝η+d-总水深；u、v-x、y方向上的流速；-x、y方向上的垂线平均速度；P_a-大气压强；ρ-水流密度；ρ₀-水密度的参考；S-点源处的流量的大小；g-重力加速度；u_s、v_s-点源的流速；τ_sx、τ_sy、τ_bx、τ_bx-x、y方向上的表面风速和底部应力；f-科氏参数，f＝2ωsinφ，ω为地球旋转角速度，φ为纬度；T_xx、T_xy、T_yy一横向应力；A为水平涡粘系数；c_s为常数；l为网格特征长度；为垂线平均含沙量；D_x、D_y分别为x、y方向上的扩散系数；Q_L为单位面积上的点源流量；C_L位点源含沙量；E为底部泥沙冲、淤项；

步骤B，建立浮泥形成条件定量研究模型，具体分步骤如下：

B-1，采用下式确定含沙水体密度与含沙量的关系：

ρ′=ρ+(1-ρρs)S---(7)]]>

上式中：ρ′为含沙水体的密度(kg/m³)；ρ为水的密度(kg/m³)，ρ_s为泥沙的密度(kg/m³)，S为含沙量(kg/m³)；

B-2，采用异重流发生判别公式来定量研究淤泥质海岸浮泥的形成条件，所述异重流发生判别公式如下：

U0ΔρgH/ρ′≤α---(8)]]>

上式中：U₀为垂线平均流速(m/s)，H为水深(m)；Δρ为水体的密度差(kg/m³)；α为浑水的弗劳德系数，α的取值范围是0.65～1.0；g为重力加速度(m/s²)；

B-3，定义淤泥质海岸异重流发生时刻的水体垂线平均含沙量作为浮泥即将形成的临界含沙量，将式(7)代入式(8)，得到临界含沙量计算公式(9)，如下：

S≥(U02α2gHρ+ρs-ρρsS0)/(1-U02α2gH)ρs-ρρs---(9)]]>

式中：S₀为背景含沙量；

步骤C，计算分析步骤：

采用公式(3)，根据步骤A建立的二维潮流泥沙数学模型，获取水动力参数，并选取合适的背景含沙量条件，根据步骤B-3得到的临界含沙量计算公式，分析计算不同航道设计尺度下浮泥的形成条件。