-
公开(公告)号:CN109242746B
公开(公告)日:2021-08-03
申请号:CN201811027603.5
申请日:2018-09-04
Applicant: 武汉大学
IPC: G06Q50/26
Abstract: 本发明提供基于应急监测数据的一维瞬时点源污染源溯源方法,其特征在于,包括以下步骤:步骤1.监测得到两个不同时刻t1和t2各观测点污染物浓度观测值并对浓度值C进行对数处理得到步骤2.将基本解析解进行对数变形,得到由组合参数表达的解析解的形式;步骤3.利用最小二乘法拟合数据,根据已经过对数处理后的数据求得法方程中的相应矩阵;步骤4.根据法方程组,求解出相应于第一个时刻t1浓度应急监测数据的组合参数α(α0,α1,α2)和相应于第二个时刻t2浓度应急监测数据的组合参数β(β0,β1,β2);步骤5.基于组合系数α(α0,α1,α2)和β(β0,β1,β2)以及步骤2中的组合参数的表达式,反解出各个基本的源项参数。
-
公开(公告)号:CN107818240A
公开(公告)日:2018-03-20
申请号:CN201711124429.1
申请日:2017-11-14
Applicant: 武汉大学
IPC: G06F19/00
CPC classification number: G06F19/00
Abstract: 本发明提供一种基于伴随同化的纵向离散系数反演方法,其特征在于,包括以下步骤:步骤1.给定上下游污染物浓度观测值C1obj(t)、C2obj(t)、河流水质模型相关参数流速u和时间t,假定纵向离散系数的初始值E0以及设定目标函数精度要求tol;步骤2.将E0代入河流水质模型中,得到计算值与观测值的差距J作为目标函数;步骤3.反向积分伴随方程,计算出目标函数关于纵向离散系数的梯度,根据该梯度采用最速下降法得到控制变量的搜索方向▽J与步长α,优化纵向离散系数变量;步骤4.将优化后的纵向离散系数代入河流水质模型方程进行计算得到下游观测断面x2处污染物浓度值分布C2(t),再次计算目标函数J,检查是否满足精度要求,否则返回步骤2,继续循环直到满足精度要求。
-
公开(公告)号:CN107895081A
公开(公告)日:2018-04-10
申请号:CN201711124437.6
申请日:2017-11-14
Applicant: 武汉大学
Abstract: 本发明提供一种多个瞬时污染源的溯源方法,其特征在于,包括以下步骤:步骤1.给定观测点污染物浓度观测值Cjobj、河流水质模型相关参数,假定i个瞬时污染源排放强度的初始值Mi(0)(i=1,2,3,…m1)、设定目标函数精度要求tol;步骤2.将Mi(0)代入河流水质模型中,得到计算值与观测值的差距J作为目标函数;步骤3.反向积分伴随方程,计算出目标函数关于瞬时污染源的排放强度Mi的梯度,采用最速下降法得到瞬时污染源的排放强度的搜索方向▽J与步长α,优化瞬时污染源的排放强度;步骤4.将优化后的各个瞬时污染源排放强度代入河流水质模型方程进行计算得到下游观测断面处污染物浓度值分布Cj,再次计算目标函数J,检查是否满足精度要求,否则返回步骤2,继续循环直到满足精度要求。
-
公开(公告)号:CN109242746A
公开(公告)日:2019-01-18
申请号:CN201811027603.5
申请日:2018-09-04
Applicant: 武汉大学
IPC: G06Q50/26
Abstract: 本发明提供基于应急监测数据的一维瞬时点源污染源溯源方法,其特征在于,包括以下步骤:步骤1.监测得到两个不同时刻t1和t2各观测点污染物浓度观测值并对浓度值C进行对数处理得到步骤2.将基本解析解进行对数变形,得到由组合参数表达的解析解的形式;步骤3.利用最小二乘法拟合数据,根据已经过对数处理后的数据求得法方程中的相应矩阵;步骤4.根据法方程组,求解出相应于第一个时刻t1浓度应急监测数据的组合参数α(α0,α1,α2)和相应于第二个时刻t2浓度应急监测数据的组合参数β(β0,β1,β2);步骤5.基于组合系数α(α0,α1,α2)和β(β0,β1,β2)以及步骤2中的组合参数的表达式,反解出各个基本的源项参数。
-
-
-