[发明专利]一种基于逐次逼近法的河道断面自动插补方法

权利要求书

1.一种基于逐次逼近法的河道断面自动插补方法，其特征在于，包括以下步骤：

步骤1：以观测断面水位与流量水文要素的纳什系数之和最大为目标函数，以水位、流量误差和断面间距为约束条件，建立基于一维非恒定流数学模型的河道断面自动插补优化模型；

步骤2：根据河道的长度、宽度和比降，确定河道断面自动插补优化模型的参数取值范围；所述参数包括河道糙率、插补断面个数、插补断面位置、插补断面形状参数，其中插补断面形状参数包括断面底宽、河底高程和边坡系数；

步骤3：根据基于逐次逼近法的并行求解方法，采用以形定位计算模式或以位定形计算模式对河道断面自动插补优化模型进行求解；

所述步骤3中的以形定位计算模式是先优化插补断面形状，再确定插补断面位置：首先假定插补断面等距分布，然后按照目标函数对各插补断面的形状进行优化计算，得到断面底宽、河底高程和边坡系数，最后对各插补断面的断面间距进行寻优计算，得到各插补断面最优的间距；具体包括以下步骤：

步骤A1：对参数进行初始化，确定各参数的取值范围和各参数的个数k_i，i＝1,2,…,6；

步骤A2：确定河道糙率取值；

步骤A3：确定插补断面个数k₁，等间距布置插补断面，并对各插补断面按照从下游到上游的顺序进行编号；

步骤A4：令断面形状寻优迭代次数L＝1，插补断面j＝1；

步骤A5：将插补断面j的个形状参数组合方案代入模型，结合河道断面自动插补优化模型进行并行寻优计算，再转到步骤A6；

步骤A6：固定已优化的插补断面形状参数，再转到步骤A7；

步骤A7：判断插补断面j的寻优编号是否等于插补断面个数k₁，若否，令j＝j+1，再转到步骤A5；若是，则得到各插补断面形状优化结果和本次形状迭代目标函数值F_L，再转到步骤A8；

步骤A8：判断断面形状寻优迭代次数L是否大于1，若否，令L＝L+1，再转到步骤A5；若是，则转到步骤A9；

步骤A9：判断F_L和上次形状迭代目标函数值F_L-1两次迭代结果是否满足收敛条件，若否，令L＝L+1，再转到步骤A5；若是，则得到当前河道糙率和插补断面个数下的各插补断面形状的最优方案；

步骤A10：令断面位置寻优迭代次数M＝1，插补断面j＝1；

步骤A11：固定各插补断面形状，对插补断面j的k₂个位置参数方案代入模型，结合河道断面自动插补优化模型进行并行寻优计算；再转到步骤A12；

步骤A12：固定已优化的插补断面位置参数；再转到步骤A13；

步骤A13：判断插补断面j是否等于k₂，若否，令j＝j+1，再转到步骤A11；若是，则得到各插补断面位置优化结果和本次位置迭代目标函数值F_M，再转到步骤A14；

步骤A14：判断断面位置寻优迭代次数M是否大于1，若否，令M＝M+1，再转到步骤A11；若是，则转到步骤A15；

步骤A15：判断F_M和上次位置迭代目标函数值F_M-1两次迭代结果是否满足收敛条件，若否，令M＝M+1，再转到步骤A11；若是，则得到当前河道糙率和插补断面个数下的各插补断面形状和位置最优方案，再转到步骤A16；

步骤A16：依次改变插补断面数量和河道糙率取值，重复上述步骤A1～步骤A15，得到不同插补断面数量和河道糙率取值方案下的各插补断面形状和位置最优方案；

所述步骤3中的以位定形计算模式是先确定插补断面位置，再优化插补断面形状：首先设定数量足够多且形状相同的插补断面，并将断面形状设定为矩形，然后仅改变各插补断面的河底高程进行寻优计算，对河底高程相近的断面进行合并，以此选出具有代表性的位置，最后再对各插补断面形状进行寻优，得到各插补断面的最优形状参数；具体包括以下步骤：

步骤B1：对参数进行初始化，确定各参数的取值范围和各参数的个数k_i，i＝1,2,…,6，的取值；

步骤B2：确定河道糙率取值；

步骤B3：确定矩形插补断面个数k_矩，等间距布置插补断面，并对各插补断面按照从下游到上游的顺序进行编号；

步骤B4：令断面位置寻优迭代次数M＝1，插补断面j＝1；

步骤B5：将插补断面j的河底高程计算方案代入模型，结合河道断面自动插补优化模型进行并行寻优计算，再转到步骤B6；

步骤B6：固定已优化的插补断面河底高程，再转到步骤B7；

步骤B7：判断插补断面j的寻优编号是否等于矩形插补断面个数k_矩，若否，令j＝j+1，再转到步骤B5；若是，则得到各插补断面的河底高程优化结果和本次位置迭代目标函数值F_M，再转到步骤B8；

步骤B8：判断断面位置寻优迭代次数M是否大于1，若否，令M＝M+1，再转到步骤B5；若是，则转到步骤B9；

步骤B9：判断F_M和上次位置迭代目标函数值F_M-1两次迭代结果是否满足收敛条件，若否，令M＝M+1，再转到步骤A5；若是，则得到当前河道糙率和插补断面个数下的各插补断面的河底高程变化情况，合并河底高程相近的断面，选出具有代表性的位置；

步骤B10：令断面形状寻优迭代次数L＝1，插补断面j＝1；

步骤B11：确定插补断面个数k₁，将插补断面j的形状参数方案代入模型，结合河道断面自动插补优化模型进行并行寻优计算，再转到步骤B12；

步骤B12：固定已优化的插补断面形状参数，再转到步骤B13；

步骤B13：判断插补断面j是否等于k₁，若否，令j＝j+1，再转到步骤B11；若是，则得到各插补断面位置优化结果和本次形状迭代目标函数值F_L，再转到步骤B14；

步骤B14：判断断面形状寻优迭代次数L是否大于1，若否，令L＝L+1，再转到步骤B11；若是，则转到步骤B15；

步骤B15：判断F_L和上次形状迭代目标函数值F_L-1两次迭代结果是否满足收敛条件，若否，令L＝L+1，再转到步骤B11；若是，则得到当前河道糙率和插补断面个数下的各插补断面形状和位置最优方案；

步骤B16：依次改变插补断面数量和河道糙率取值，重复上述步骤B1～步骤B15，得到不同插补断面数量和河道糙率取值方案下的各插补断面形状和位置最优方案。