[发明专利]一种资料贫乏地区长系列泥沙数据综合重构方法

权利要求书

1.一种资料贫乏地区长系列泥沙数据综合重构方法，其特征在于：包括以下具体步骤，

1)采集目标河段水文、泥沙数据，具体包括：水文站长系列逐日流量Q_i,j,k，泥沙测验含沙量TS_i,n,k与测时流量q_i,n,k，水质测验中表层泥沙测验含沙量C_i,m,k和测时流量cq_i,m,k，i＝1～Nyear，Nyear为长系列年数；j＝1～Nday(i)，Nday(i)为i年总天数；k＝1～Nsta，Nsta为水文站数；n＝1～Nsam(i,k)，Nsam(i,k)为第i年k号水文站泥沙测验样本数；m＝1～NCsam(i,k)，NCsam(i,k)为第i年k号水文站水质测验表层泥沙测验样本数；

2)采用流量-含沙量年关系曲线法计算资料优良年份泥沙数据；

3)采用水质泥沙数据与水文站泥沙数据相关关系法计算仅水质泥沙测验年份泥沙数据；

4)采用临近站点同年流量-含沙量关系曲线法计算资料缺乏年份泥沙数据；

5)采用多年平均流量-含沙量关系曲线法计算其余年份泥沙数据；

所述步骤2)采用流量-含沙量年关系曲线法计算资料优良年份泥沙数据的具体步骤如下：

依据以下原则识别泥沙数据资料优良年份：①汛期各个月均有泥沙测验样本；②年内样本数量不少于某一特定数量Nsmin；

采用泥沙数据资料优良年份泥沙测验资料TS_i,n,k与q_i,n,k，建立年内流量-含沙量关系曲线：

TS_i,k＝f₁(q_i,k) (1)

采用显著性水平α＝0.05且p＜0.01检验所建立年内流量-含沙量关系曲线能否通过t检验，如不能，则将该年重新划出泥沙数据资料优良年份；反之，即可根据式(1)计算该年逐日含沙量数据：

S_i,k＝f₁(Q_i,k) (2)

并统计年输沙量特征值：

所述步骤3)采用水质泥沙数据与水文站泥沙数据相关关系法计算仅水质泥沙测验年份泥沙数据的具体步骤如下：

在泥沙数据资料优良年份中选出同时拥有水质测验表层泥沙测验样本的年份，采用这些年份表层泥沙测验含沙量C_i,m,k和测时流量cq_i,m,k，建立表层泥沙含沙量与流量的关系曲线：

C_i,k＝f₂(cq_i,k) (4)

根据式(4)计算该年表层泥沙年输沙量特征值

根据式(3)和(5)即可建立k水文站，泥沙测验资料优良且有水质表层泥沙测验资料年份，其年输沙量与年表层输沙量的关系曲线：

TQS_k＝f₃(TQC_k) (6)

根据式(6)和式(4)即可计算有水质表层泥沙测验资料但无优良泥沙测验资料年份的逐日泥沙数据：

S_i,j,k＝f₃(f_2,(Q_i,j,k)×Q_i,j,k)/Q_i,j,k (7)；

所述步骤4)采用临近站点同年流量-含沙量关系曲线法计算资料缺乏年份泥沙数据的具体步骤如下：

在相邻k和k+1两水文站各自泥沙数据资料优良年份中，选出两站同年均为资料优良的年份，根据式(3)建立相邻两站年输沙量的关系曲线：

TQS_k＝f₄(TQS_k+1) (8)

采用显著性水平α＝0.05且p＜0.01检验所建立关系曲线(8)能否通过t检验，如未通过检验，则相邻两站间并不存在明显的相关关系，无法应用该方法；如通过检验，则可由式(8)和k+1水文站资料优良年份水沙关系曲线式(2)，计算k站同年且无优良资料年份的泥沙数据：

S_i,j,k＝f₄(f_1,k+1(Q_i,j,k)×Q_i,j,k)/Q_i,j,k (9)；

所述步骤5)采用多年平均流量-含沙量关系曲线法计算其余年份泥沙数据的具体步骤如下：

采用所有资料优良年份泥沙测验资料TS_i,n,k与q_i,n,k，建立多年平均流量-含沙量关系曲线：

TS_k＝f₅(q_k) (10)

对于步骤2)～步骤4)三种方法均无法适用的年份，采用式(10)计算逐日含沙量数据：

S_k＝f₅(Q_k) (11)。