黄河口泥沙异重流基本控制参数的数值试验

《黄河口泥沙异重流基本控制参数的数值试验》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、党委成员工作鉴定根据镇党委分工，***同志协助镇长**泥沙异重流基本控制参数对泥沙异重流流动的影响，本文设计了一系列数值试验来展示各种参数与泥沙异重流流动特征之间的定量关系。采用ADI法求解三维方程经垂向积分后的平面二维泥沙异重流方程。所选取的基本控制参数包括水下斜坡坡角总阻力系数，水力挟带系数，泥沙饱和浓度和泥沙沉速。通过数值试验得出水下斜坡坡角0、泥沙饱和浓度Es越大，总阻力系数cd、泥沙沉速o越小，相应地泥沙异重流流速、流动厚度以及悬沙浓度越大。其中斜坡坡角对泥沙异重流影响最大。当泥沙异重流流动处于超临界流状态，水力挟带系数Ew

2、大于零，这时由于有一定水量的卷入，所以流动厚度沿程增加。关键词：黄河口泥沙异重流控制参数数值试验ADI法泥沙异重流及形成的地质体在沉积学和石油地质中有重要意义。在黄河口，泥沙异重流是泥沙搬运的重要方式，是河口区沉积动力的重要过程，是认识河口三角洲发育演化的重要途径。为使研究深入一步，我们进行了泥沙异重流的数值模拟。近年来，国内对异重流数模研究较少，仅见方春明等（1997）关于泥沙异重流潜入时的立面二维数值模拟以及邱晨霞（1995）对盐水异重流的二维两层数值计算。1平面二维泥沙异重流控制方程及数值模拟目前，包含三个基本控制方程(流体质量

3、、动量守恒方程和泥沙质量守恒方程)的泥沙异重流数值模拟被广泛应用，本文把一维方程推广到平面二维的情形，成为四方程模型，如下所示［1］(1)(2)(3)(4)其中u,v分别为X，y方向速度，0x，9y分别为X，y方向坡度，Ewx与Ewy分别为x、y方向的水力挟带系数，TX、Ty为X方向与y方向的剪切摩擦力，kx、ky分别为X方向与y方向的悬沙扩散系数，t为时间，h为异重流厚度，s为深度平均的悬沙含量，ps，pw，Pf分别为泥沙颗粒、水体和泥沙异重流的密度，Es为泥沙饱和浓度，《为泥沙颗粒的沉降速度，a为泥沙运动饱和系数，g为重力加速度，

4、f为柯氏力参量。在方程(1)中，水力挟带系数Ewx与Ewy是流动Richardson数(Rix、Riy)的函数，可用若干经验式表迗［5］。以x方向为例，它们是Ashidaetal.(1975):Ewx=/Rix;Parkeretal.(1987):Ewx=/(1+718Rix)；Fukushimetal.(1985):Ewx=/(+Rix);以及俞维升(1991):Ewx=/Rix等。在方程(2)、(3)中，剪切摩擦力TX、Ty是作用在异重流上下界面摩擦力的总和，它们由下式所示。其中cd=cf(l+a')，被称为总阻力系数，范家骅等(

5、1980)的cd值约，Parkeretal.(1987)的cd值在〜之间，舍维升(1991)的cd值在〜之间，可见cd值有很大的不确定性。在方程(4)中，泥沙饱和浓度Es是描述泥沙异重流流动特征参数的函数，对黄河口可用下式表示(中国水科院(1997)公式[2，3]s<15kg/m3,Es=123(V2/gh)(V/o)-(1-Ah/h)s〉15kg/m3,Es=(V2/gh)(V/o)(1-Ah/h)V为合流速，h为水深，《为泥沙沉速，Ah为潮差。张青玉经验式Es=(V2/h)定解条件：边界条件：陆地边界：Vn=0(n为岸线的法线方向

6、)，s'=0('表示陆地)水边界：h*(t)=h(t)，s*=s(*表示水边界);河口边界:u=(m/s)s=50(kg/m3)初始条件：u=0v=0s=0h=(m)采用有限差分法中的ADI法求解方程组(1)、(2)、(3)、(4)。事实上，我们应当同时求解描述泥沙异重流流动的方程，因为它们是一组相互耦合的方程组。然而同时求解相当复杂，我们便采取一种叫冻结系数法［4］的近似耦合法来求解。计算黄河口泥沙异重流时，泥沙异重流起点在黄河入海处。数值计算区域如图1所示，X、y方向空间步长均取1/3千米，时间步长取24秒；计算区海底地形如图2所

7、示；柯氏参数f=2Qsind），Q取X10-5s-l，（1）取°;泥沙密度Ps取/cm3,水体密度Pw取/cm3,泥沙异重流密度Pf取/cm3。1泥沙异重流数值计算区域Areaofnumericalsimulationofsedimentdensitycurrent2计算区海底地形图Mapofareaofnumericalsimulation2平面二维泥沙异重流的参数试验为了定量地评价基本控制参数对泥沙异重流流动的影响，本文设计了一系列数值试验（或称敏感性试验）来展示各种参数与泥沙异重流流动特征之间的定量关系。所选取的基本控制参数包括

8、水下斜坡坡角0x、0y，总阻力系数cd，水力挟带系数Ew，泥沙饱和浓度Es和泥沙沉速《，进行数值试验时它们各自的输入值见表1。以试验13为标准进行对比，试验13计算结果如图3。表1数值试验中基本控制参数的输入值Input