一种海水溶解氧评估方法与流程

1.本发明涉及建筑能耗监测领域，尤其涉及一种海水溶解氧评估方法。


背景技术：

2.现有海水溶解氧评估方法：倒数方法、折线法；这两个方法是现有的海洋学中评价海水溶解氧的方法，但是仅适用于标准测量方式下的海水溶解氧数据，难以应用于其他方式测得含氧量。现有海水溶解氧评估方法存在结果与其他污染物指数不匹配的问题，难以考虑其他污染物指数对溶解氧的影响。
3.现有的可视图方法包括：基本可视图方法、水平可视图方法、有限穿越可视图方法、有向有限穿越可视图方法、过滤水平可视图方法、多尺度有限穿越水平可视图方法、时间相关有限穿越可视图、加权可视图。现有可视图系列方法均未同时考虑时间演变及变量间相互联系。
4.现有图相似性匹配方法：基于矩阵分解的图相似性匹配方法、基于深度学习的具有随机游走的图相似性匹配方法、基于深度学习的没有随机游走的图相似性匹配方法、基于生成模型的图相似性匹配方法、基于图核的图相似性匹配方法。
5.基于矩阵分解的图相似性匹配方法时空开销较大；基于深度学习的具有随机游走的图相似性匹配方法很难找到最优采样策略,没有随机游走的图相似性匹配方法计算成本较高；基于生成模型的图相似性匹配方法需大量训练数据；基于图核的图相似性匹配方法,主要利用局部子结构信息得到一维图特征向量,相比于前三类方法计算更快捷。其中子树模式图核的实例weisfeiler
‑
lehman子树核在解决图相似性匹配问题时更高效，但仅包含子图信息，考虑信息不周全。


技术实现要素：

6.有鉴于现有技术的上述缺陷，本发明的目的是提供一种海水溶解氧评估方法，将时间序列转化为图，从而准确地评估海水溶解氧。
7.本发明的上述技术目的通过以下技术方案得以实现：
8.为实现上述目的本发明提供了一种海水溶解氧评估方法，其包括：
9.(s1)获取两个需要对比海域的海水溶解氧及其相关变量，采用tp有向可视图方法对海水溶解氧及其相关变量的时间序列进行转换，得到各海域的图；海水溶解氧的相关变量包括海温以及盐度；
10.(s2)对两个海域的图进行相似性匹配，得到匹配结果。
11.本发明的进一步改进在于，步骤(s1)中，采用tp有向可视图方法对海水溶解氧及其相关变量的时间序列进行转换得到各海域的图，其包括以下步骤：
12.(s11)对海水溶解氧及其相关变量的时间序列分别进行归一化操作，得到数值均处于[0,1]的多变量时间序列；
[0013]
(s12)利用皮尔逊相关系数的方法，进行海水溶解氧与各相关变量的相关性计算；
[0014]
(s13)对海水溶解氧以及每一相关变量对应时间序列应用有向可视图方法，得到图集；图集中海水溶解氧为主变量；
[0015]
(s14)结合有向概率转移矩阵、相关性计算结果与主变量图度分布的大小关系，判断各相关变量对海水溶解氧的影响，进而确定是否根据相关变量的图在主变量的图上增加边。
[0016]
本发明的进一步改进在于，步骤(s2)中采用gsn图相似性匹配方法进行相似性匹配，其具体包括以下步骤：
[0017]
(s31)利用weisfeiler
‑
lehman子树核计算子图层图间相似度；
[0018]
(s32)利用pagerank算法得到节点层图间相似度；
[0019]
(s33)利用图密度、平均度、全局聚集系数、平均聚集系数、平均最短路径长度得到全局层次图间相似度；
[0020]
(s34)利用规约余弦相似度归一化上述相似度结果，并等比例加权得到最后相似度结果。
[0021]
在本发明的较佳实施方式中，本发明提供的方法具有以下技术效果：tp有向可视图可以保留更多时间序列信息的转化成图，可更好的分析时间序列。gsn图相似性匹配从全局层、子图层、节点层综合计算图间相似度，使得相似度评估结果更加准确。
附图说明
[0022]
图1是本发明的海水溶解氧评估方法的流程图；
[0023]
图2是有向可视图生成过程的示意图；
[0024]
图3是weisfeiler
‑
lehman子树核计算子图层图间相似度的示意图；
[0025]
图4是数据对比图；从图4可知赤道太平洋较东北太平洋缺氧更严重；
[0026]
图5是tp有向可视图方法的实例示意图。左侧为海水溶解氧、温度、盐度三变量时间序列，右侧为应用tp有向可视图方法后得到的溶解氧图。
具体实施方式
[0027]
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需说明的是，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。
[0028]
如图1所示，本发明的实施例包括一种海水溶解氧评估方法，其包括：
[0029]
(s1)获取两个需要对比海域的海水溶解氧及其的相关变量，采用tp有向可视图方法对海水溶解氧及其相关变量的时间序列进行转换，得到各海域的图；海水溶解氧的相关变量包括海温以及盐度。
[0030]
在步骤(s1)中，分别对两个待比较的海域的海水溶解氧及其相关变量进行处理。本步骤中采用tp有向可视图方法对海水溶解氧及其相关变量的时间序列进行转换得到各海域的图，具体包括以下步骤：
[0031]
(s11)对海水溶解氧及其相关变量的时间序列分别进行归一化操作，得到数值均
处于[0,1]的多变量时间序列；
[0032]
(s12)利用皮尔逊相关系数的方法，进行海水溶解氧与各相关变量的相关性计算；
[0033]
(s13)对海水溶解氧以及每一相关变量对应时间序列应用有向可视图方法，得到图集；图集中海水溶解氧为主变量；
[0034]
对每个关变量对应的时间序列进行处理过程中，将某个相关变量以时间序列的方式进行表示，对于序列中的每个采样点，将其作为图的节点。对于每个采样点，在时间轴为x轴、采样值为y轴的直角坐标系中，作采样点与各后续采样点(后续采样点指的是采样时间在后的采样点)之间的有向线段，若该线段不经过其他采样点与x轴之间的垂线，则该有向线段作为有向可视图的边。得到每个相关变量的有向可视图之后，将各有向可视图的集合作为相应海域的图集。
[0035]
(s14)结合有向概率转移矩阵、相关性计算结果与主变量图度分布的大小关系，判断各相关变量对海水溶解氧的影响，进而确定是否根据相关变量的图在主变量的图上增加边。本步骤中，将各相关变量进行综合，得到的图可使得评估结果更加准确。
[0036]
步骤(s2)中，采用gsn图相似性匹配方法进行相似性匹配，其具体包括以下步骤：
[0037]
(s31)利用weisfeiler
‑
lehman子树核计算子图层图间相似度，如图3所示。weisfeiler
‑
lehman子树核方法是现有的，用于比较图间相似度，但相似度计算结果会存在可能无上限的问题，故本发明去除最后一部分的内积计算，改采用规约余弦相似度；同时weisfeiler
‑
lehman子树核方法仅利用了图的子图层信息，相似性计算结果不够全面，故增加(s32)(s33)两步。
[0038]
(s32)利用pagerank算法得到节点层图间相似度；
[0039]
(s33)利用图密度、平均度、全局聚集系数、平均聚集系数、平均最短路径长度得到全局层次图间相似度；
[0040]
(s34)利用规约余弦相似度归一化上述相似度结果(s31、s32、s33中得到的相似度结果)，并等比例加权得到最后相似度结果。
[0041]
如图4所示，实验数据中布朗序列hurst指数为0.60，经tp有向可视图方法及有向有限穿越可视图建图，对双对数坐标下的度分布图进行线性拟合，得到h＝0.70(tp有向可视图)、h＝0.81(有向有限穿越可视图)，可知tp有向可视图保留了更多原时间序列的信息。
[0042]
由2019年12月7日国际自然保护联盟发布的报告可知，赤道太平洋的缺氧程度比东北太平洋严重，验证方法的正确性，结合现有海水溶解氧评估方法可得到更好的评估结果。
[0043]
以死亡地带之墨西哥湾处海水溶解氧及相关参数的时间序列为例，构图，如图5所示。
[0044]
东北太平洋(enp,eastern north pacific)、东南太平洋(esp,eastern south pacific)、孟加拉湾(bb,bay of bengal)及赤道太平洋海域缺氧严重。在上述海域各选取一个浮标站点，利用gsn图相似性匹配方法与墨西哥湾对比得到相似度表格。
[0045][0046]
经表格数据可验证gsn图相似性匹配方法的正确性。
[0047]
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。