[发明专利]具有沉积信息释放率先验的源项反演方法及系统

权利要求书

1.一种具有沉积信息释放率先验的源项反演方法，其特征在于包括：

获取数据参数，并初始化校正系数矩阵W和源项释放率σ；

获取选定观测站点处的沉积数据；

基于观测站点处的沉积数据，制作考虑湿沉降的释放率先验；

基于带先验的联合偏差校正方法求解校正系数矩阵W和源项释放率σ直至收敛，获得输运矩阵H_new；

利用输运矩阵H_new求解源项释放率σ直至收敛，得到具有沉积信息释放率先验的源项反演结果。

2.根据权利要求1所述的具有沉积信息释放率先验的源项反演方法，其特征在于，所述数据参数包括沉积观测站点坐标向量组(x^o，y^o)、初始输运矩阵、空气浓度观测数据、正则化参数λ及收敛条件，o表示选取的站点个数，其中，收敛条件包括最大迭代次数和允许的最小误差。

3.根据权利要求2所述的具有沉积信息释放率先验的源项反演方法，其特征在于，所述获取选定观测站点处的沉积数据，包括：选定观测站点对应的降雨量、恒定释放σ_cⁿ下站点的沉积时序变化率dDep^o，t、仅空气浓度反演得到的站点沉积时序变化率aDep^o，t和地表累积沉积观测值mDep^o，其中，n代表模拟的时间步长，t为时间。

4.根据权利要求3所述的具有沉积信息释放率先验的源项反演方法，其特征在于，所述基于观测站点处的沉积数据，制作考虑湿沉降的释放率先验，包括：

获取湿沉降主导时段tR，并获得对应主导时段下恒定释放站点的沉积时序变化率dDep^o，tR和主导时段下仅空气浓度反演得到的站点沉积时序变化率aDep^o，tR；

折算各站点主导时段tR的考虑湿沉降的释放率σ_s^o，tR，并计算释放点处主导时段tR的考虑湿沉降的释放率先验σ_i^o，tR，制作考虑湿沉降的释放率先验σ_depⁿ。

5.根据权利要求4所述的具有沉积信息释放率先验的源项反演方法，其特征在于，所述折算各站点主导时段tR的考虑湿沉降的释放率σ_s^o，tR，并计算释放点处主导时段tR的考虑湿沉降的释放率先验σ_i^o，tR，制作考虑湿沉降的释放率先验σ_depⁿ，包括：

各站点主导时段tR的考虑湿沉降的释放率σ_s^o，tR的计算过程为：

式中，σ_c^tR为湿沉降主导时段tR下的恒定释放率；

释放点(x_r，y_r)处主导时段的考虑湿沉降的释放率先验σ_i^o，tR：

σ_i^tR＝σ_s^o，tR·weight^o  (5)

考虑湿沉降的释放率先验σ_depⁿ的计算公式为：

式中，σ_i^t为t＝i时的释放率，σ_c^t为恒定释放率。