宽级配卵石推移质输移随机过程-分维研究

1、宽级配卵石推移质输移随机过程-分维研究            摘要：关键词：分维 R/S分析 推移质 泥沙输移随机过程   前言图1为1958年以来，长江上游寸滩沙量发展过程时间序列。由图可以看出，年径流量、沙量、平均含量逐年均有变化，杂乱无章，毫无规律可循。对于这种时间序列，分析其时间序列的特点、规律是现代数学方法分形几何学、R/S分析的研究内容。用新的手段方法尝试分析长江上游泥沙的来水来沙条件，无非是企图从不同角度促进长江三峡工程泥沙研究工作的深入。无疑引入现代

2、数学或其它现代科学的方法去探讨长江三峡工程泥沙问题，对于宽级配卵石推移质输移特性的研究具有理论意义和实用价值。 文中应用现代数学的新成果分形几何学、R/S分析方法进行分析，从理论上探讨长江上游来水来沙发展过程和宽级配卵石推移质输沙移机过程的特点，为推移质输移特性深入研究提供基础。分形理论还可用于长江上游河床地貌、河型分类和河型转化领域的研究中。Rn,S0,对于>0，定义HS(A)=infUiS    2007-04-23        (AUUi,Ui)，其中Ui表示U

3、I的直径，定义HS(A)=limHSR(A)，HS(A)称为集合A的豪斯道夫S维测度。可以证明，对于集合A，存在唯一的非负实数DH(A)，它满足：若0<S<DH(A)，则HS(A)=;若DH(A)<S<,则HS(A)=0。DH(A)叫做A的豪斯道夫维数。对于任意的集合A，恒有DH(A)>DT(A),DT(A)为A的拓扑维数。R/S分析是由赫斯特于1965年提出的一种时间序列统计法2。它在分形理论中有着重要的应用。下面就R/S分析的基本原理和方法、分维数与赫斯特数的关系及其物理意义作简单的介绍。考虑一个时间序列(t),t=1,2，对于任意整数>1，定义非负时均

4、序列<>=1/(t),=1,2；用X(t)表示累计离差：X(t,)= (u)-<>),1<t<；极差定义为：R()=maxX(t,)-minX(t,),(1<t<),=1,2)；标准差定义为：S()=(1/(t)-<>)2)1/2,=1,2,。现考虑比值R()/S()=R/S，赫斯特发现有如下经验标度关系：R/S(/2)H（或H），H称为赫斯特数。如果(t)相互独立，方差有限的随机过程序列，则由赫斯特和费勒证明了如下结果：R/S=(/2)1/2,即H=1/2。由分维定义和分式布朗运动的原理，可得出分维数DH和赫斯特数的关系为：DH=2

5、-H或H=2-DH。当H=0.5时，DH=1.5是通常的布朗运动的分维数，H<1/2时，呈现负相关，有两种情况，一则由“好”到“坏”，二则由“坏”到“好”。当其H>1/2时，呈现正相关，可能出现“良性循环”和“恶性循环”。决策者采取适当的措施，避“恶”扬“良”，达到较好的目的。显然，这在经济学中有极为重要的应用前景对国民经济宏观决策作出良好的综合评价。同时，若(t)为一学生成绩，则R/S分析结果可以评价学生成绩发展的综合情况。同样，通过赫斯特数的大小，可以判断江河水量、沙量和淤积量发展趋势。根据H>1/2或H<1/2，提出流域水土保持治理决策。  &

6、#160; 3 长江上游沙量发展的R/S分析长江上游50年代以来沙量过程如图2所示。用上述1958年以来的水沙量，含沙量资料进行R/S分析计算，其R/S和的关系如图2所示。长江上游寸滩推移质沙量过程及R/S分析结果如图3所示。 长江上游悬移质沙量发展过程的：R/S=0.87890.679，回归相关系数为0.9808。由此可得赫斯特数    2007-04-23        H=0.679，分维数DH=1.321。长江上游推移质沙量发展过程的R/S=0.93160.

7、622，回归相关系数为0.9566。由此可得赫斯特数H=0.622，分维数DH=1.378。图3 长江寸滩推移质沙量过程及R/S分析结果The process and R/S relation of bed-load transport in the upper Yangtze River由长江上游悬沙、推移质沙量过程的H>0.5可见，长江上游沙量发展呈正相关，过去的增长可能意味着将来的增长，过程具有持久记忆性，沙量发展可能出现“良性循环或恶性循环”。曼德尔布罗特曾对太阳黑子活动进行R/S分析，得到H=0.93，远大于1/2，这与太阳黑子活动周期性（11年一个周期）相一致。沙量发展过程的

8、H值大于0.5，有着记忆性和周期性。不能以某一时段的水沙规律说明其沙量是增加和减少，来制定长期治理的策略。长江上游沙量发展可能出现“良性循环或恶性循环”的可能，因此应依靠长久的水土保持和水库拦截等措施，达到良性循环，避免出现“恶性循环”的可能。 宽级配推移质输沙率随机过程的分维研究宽级配推移质输移的一个特点是输沙率具有强烈的脉动性（阵发性）。图、图、图和图均是岷江宝瓶口河段实测推移质输沙过程。同一流量所对应的输沙率值相差可达数十倍至上百倍，这种大幅度变化的周期，可长达数十分钟4,5,8。图4、图5为间隔89分钟采样的输沙率脉动过程。而图6、图7为间隔3分钟采样的输沙率脉动过程。下面分别对图4、

9、图5、图6和图7所示的输沙率脉动过程作分维计算，所得RS与的关系曲线如相应的图4、图5、图6和图7所示。    2007-04-23        图4所示的输沙率随机过程的：R/S=0.9490.5315。回归相关系数为0.9744。由此可得赫斯特数H=0.5315，分维数DH=1.4685。图7 卵石推移质输沙率随机过程及R/S分析结果The process and R/S relation of bed-load rate measured with 3 minutes

10、 intervals图5所示的输沙率随机过程的：R/S=1.0240.5355。回归相关系数为0.9766。由此可得赫斯特数H=0.5355，分维数DH=1.4645。由此可得间隔89分钟采样的输沙率脉动过程的分维数DH=1.47，接近马氏过程DH=1.5。这和以往假设泥沙跃移、休止时间等马氏过程是一致的。然而，间隔3分钟采样的输沙率脉动过程的分维数与上述结果就存在差别。图6所示的输沙率随机过程的：R/S=0.8060.6737。回归相关系数为0.9806。由此可得赫斯特数H=0.7637，分维数DH=1.3267。图7所示的输沙率随机过程的：R/S=0.81250.6605。回归相关系数为0

11、.978。由此可得赫斯特数H=0.6605，分维数DH=1.3395。由此可得间隔3分钟采样的输沙率脉动过程的分维数DH=1.331.34，H=0.660.67（大于0.5），表现出持久记忆性。两组过程的差别在于采样间隔不同，导致输沙率随机过程特性也不相同。这种由时间间隔不同，其测量的随机过程特性不同，与海洋波高研究成果也是一致的。在海洋波高预测研究中，若以=10天左右为间隔，则H=0.92，表现出持久记忆性；若以=20天左右为间隔，则H=0.52，表现出统计独立性。实际工程应用中，根据沙量发展具有持久记忆性特点制定长期治理措施。而推移质输沙率随机过程的独立性与持久记忆性取决于采样时间间隔，这对进一步探讨推移质输移脉动和采样具有指导