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

1、黄河口泥沙异重流基本控制参数-数值试验            摘要：为了定量地评价泥沙异重流基本控制参数对泥沙异重流流动的影响，本文设计了一系列数值试验来展示各种参数与泥沙异重流流动特征之间的定量关系。采用ADI法求解三维方程经垂向积分后的平面二维泥沙异重流方程。所选取的基本控制参数包括水下斜坡坡角，总阻力系数，水力挟带系数，泥沙饱和浓度和泥沙沉速。通过数值试验得出水下斜坡坡角、泥沙饱和浓度Es越大，总阻力系数cd、泥沙沉速越小，相应地泥沙异重流流速、流动厚度以及悬沙浓度越大

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

3、制方程及数值模拟目前，包含三个基本控制方程(流体质量、动量守恒方程和泥沙质量守恒方程)的泥沙异重流数值模拟被广泛应用，本文把一维方程推广到平面二维的情形，成为四方程模型，如下所示1(2)(4)其中u，v分别为x，y方向速度，x，y分别为x，y方向坡度，Ewx与Ewy分别为x、y方向的水力挟带系数，x、y为x方向与y方向的剪切摩擦力，kx、ky分别为x方向与y方向的悬沙扩散系数，t为时间，h为异重流厚度，s为深度平均的悬沙含量，s，w，f分别为泥沙颗粒、水体和泥沙异重流的密度，Es为泥沙饱和浓度，为泥沙颗粒的沉降速度，为泥沙运动饱和系数，g为重力加速度，f为柯氏力参量。在方程(1)中，水力挟带系

4、数Ewx与Ewy是流动Richardson数(Rix、Riy)的函数，可用若干经验式表达5。以x方向为例，它们是Ashida et al.(1975): E    2007-04-23        wxx；Parker et al.(1987): Ewx=0.075/(1+718Rix)；Fukushim et al.(1985):Ewx=0.0015/(0.0204+Rix)；以及俞维升(1991)：Ewxx等。在方程(2)、(3)中，剪切摩擦力x、y是作用在异重流上下界

5、面摩擦力的总和，它们由下式所示。图3 数值试验13的计算结果图(本图的坐标与图1的坐标一致)Calculating results of numerical test No. 132.1 水下斜坡坡角的影响水下斜坡坡角越大，有效重力作用越强，泥沙异重流加速越快，可以运行更远的距离，而且可以一直保持较高的含沙量。通过数值试验13、1、2、3可揭示坡角大小对泥沙异重流流动产生的显著影响。这些数值试验的控制参数除坡角不同外，其余参数均相同,见表1。流速变化图4显示坡角越大,流动速度越大。在河口以东2km处(图4)，试验13的流速为0.80m/s，试验1、2、3的流速分别为1.41m/s, 1.82m

6、/s和2.67m/s。数值试验13(如图3(a)，由于其输入真实海底地形,泥沙异重流从黄河口流出后，便向周围扩散，这是由于实际上口门处地形最高，无论向北向南，还是向东均有一定的坡度。流动厚度变化图5同样显示，坡角越大，其流动厚度越大。试验13，1，2，3的0.8m等厚度线所括范围占整个计算区域的百分比分别为25.0, 66.6, 75.0, 83.3。悬沙含量变化图6显示，坡角越大，愈易保持高的悬沙含量，在河口以东2km处，试验13，1，2，3的悬浮含量分别为18.5，22.3，31.6，37.9g/l。2.2 总阻力系数cd的影响总阻力系数cd的大小主要与边界条件有关，对于不同边界不同性质的

7、泥沙异重流，其总阻力系数cd的值很不相同。cd值大小对泥沙异重流影响很大，cd越大，泥沙异重流在运动过程中动量损失越大，越易发生沿程沉降而导致快速消亡。本文设计了数值试验4、5、13、7来定量分析cd对泥沙异重流流动的影响。cd值分别是0.001、0.005、0.01和0.05，依次增大。图4显示随着cd值的增大，流速逐渐减小，在河口以东2km处，试验4的流速为1.36m/s，试验5，13的流速分别为1.05m/s，0.80m/s；试验7的cd值高达0.05，几乎是所量测到的cd值中的最大值，在如此大的阻力下，流速快速衰减，在2km处流速已降至0.32m/s。流动厚度变化图5 

8、60;  2007-04-23        表明，流动厚度的衰减速度随c3。2.3 水力挟带系数Ew的影响当泥沙异重流流动处于超临界流状态，水力挟带系数Ew大于零，这时有一定水量的卷入；当泥沙异重流处于亚临界流状态，水力挟带系数Ew等于零，甚至小于零出现负挟带。黄河口的泥沙异重流一般情况下，其流速小于1.5m/s,悬沙浓度小于60g/l,属低浓度泥沙异重流，通常处于亚临界状态，因此在本文的数值试验中，水力挟带系数Ew大部分取零值，只有试验6和试验8中的Ew分别取0.0015Ri和0.00034

9、Ri。三个试验的流速和悬沙含量基本类似，说明有无水力挟带项以及水力挟带系数的大小对流速和悬沙含量的分布影响不大；然而三个试验的流动厚度大相径庭(图7)，试验13的流动厚度沿程减小，而试验6，试验8的流动厚度沿程增大，Ew值越大，卷入的水量便越多，试验8的Ew略大于试验6的Ew，因此，试验8流动厚度的增加略快于试验6。图4 泥沙异重流流速(m/s)变化图(河口以东2km处的流速)Variation of velocity of sediment density current图5 泥沙异重流流动厚度变化图(0.8m等厚度线所括范围与整个计算区域之比)Variation of flow depth

10、 of sediment density current图6 泥沙异重流悬沙浓度（g/l）变化图(河口以东2km处的浓度)Variation of sediment concentration of sediment density current图7 Ew值大小对泥沙异重流流动厚度(m)的影响注：图中阿拉伯数字代表数值试验序号Variation of density current depth influenced by Ew2.4 泥沙饱和浓度Es的影响    2007-04-23      

11、;  泥沙饱和浓度Es是影响泥沙异重流数值模拟的关键因素，Es越大，泥沙异重流越易保持而不消亡。本文设计了试验9与试验13来定量分析Es不同对泥沙异重流所造成的影响。试验9的Es用水科院关于黄河口输沙力的公式来表达，试验13的Es用张青玉公式表达，除Es不同外，两个试验的其余参数均相同。流速变化图4清楚地显示，当Es值增大时，流速增大，口门以东2km处的流速从0.56m/s增到0.80m/s。在流动厚度变化图5中，试验9的0.8m等值线范围仅占整个区域的1/10，而试验13为1/4左右，即流动厚度增加了2.5倍左右；如悬沙含量图6所示，试验9的悬沙含量在河口以东不到2km处

12、已降至15kg/m33，说明泥沙饱和浓度越小，泥沙异重流越易发生沿程沉降而导致其快速消亡。2.5 泥沙沉速的影响本文进行了数值试验10、11、13和12来定量分析泥沙沉速的影响，四次试验的沉速取值分别是0.475cm/s、0.086cm/s、0.036cm/s和0.020 cm/s。泥沙颗粒沉速的大小对泥沙异重流流速分布几乎无影响，图4显示三个试验11，12，13的流速相差不多，只是试验10的流速与其它3个相比偏小，是由于试验10的沉速值是其它3个的1020倍相差很大造成的。而流动厚度和悬沙浓度对沉速值大小却相当敏感，沉速越大，流动厚度越小，相应空间点的悬沙浓度亦越小。尤其对于试验10，由于沉

13、速过大(只是一假设值)，悬沙含量在河口以东2km处已降至零，说明不到2km悬浮泥沙已全部沉降。3 结论通过数值试验及分析得出水下斜坡坡角对泥沙异重流影响最大，即泥沙异重流强度对斜坡坡角改变相当敏感，且呈非线性关系，类似于抛物线型。总阻力系数cd值大小对泥沙异重流影响次之。泥沙颗粒沉速大婿对泥沙异重流悬沙浓度有较大影响，对泥沙异重流的流速几乎无影响。当泥沙异重流流动处于超临界流状态，水力挟带系数Ew大于零，这时由于有一定水量的卷入，所以流动厚度沿程增加；而Ew对流速及悬沙浓度几乎无影响。泥沙饱和浓度Es是影响泥沙异重流数值模拟的关键因素，Es越大，泥沙异重流越易保持而不消亡。    参