CN105404709B - 基于复杂网络的堤坝健康监测敏感测点分析方法 - Google Patents
基于复杂网络的堤坝健康监测敏感测点分析方法 Download PDFInfo
- Publication number
- CN105404709B CN105404709B CN201510691360.5A CN201510691360A CN105404709B CN 105404709 B CN105404709 B CN 105404709B CN 201510691360 A CN201510691360 A CN 201510691360A CN 105404709 B CN105404709 B CN 105404709B
- Authority
- CN
- China
- Prior art keywords
- measuring point
- dykes
- dams
- complex network
- measured value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Architecture (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了基于复杂网络的堤坝健康监测敏感测点分析方法,包括步骤有,根据堤坝健康监测项目,通过对应每一监测项目实测数据,得到反应工程健康状况的测值时间序列;运用等概率粗粒化的方法,将测值时间序列转化成符号序列;把符号序列编码在复杂网络的拓扑结构中,构造复杂网路;通过对复杂网络的特征参数进行分析,获取该堤坝健康监测项目的敏感测点及堤坝的性态演化规律。所提出的新方法具有有效识别重要测点和判别堤坝性态演化的能力,克服常规分析方法对数据误差的敏感性,可多尺度地发现各个测点的重要程度,实现各个测点之间相关关系的分析,并及时准确地掌握堤坝工程健康性态,实现堤坝性态的更有效预警。
Description
技术领域
本发明涉及一种堤坝健康监测分析方法,特别是涉及一种基于复杂网络的堤坝健康监测敏感测点分析确定方法。
背景技术
堤坝健康监测不仅有利于避免大的险情和灾害发生,对堤坝性态及微小变化的准确掌握也是反馈设计施工质量和检验科研效果的前提条件。及时、准确地掌握堤坝的性态,是保证堤坝安全运行的基础。
堤坝健康监测是河湖和水库安全运行的耳目,但目前,绝大部分实测资料分析方法对数据要求比较高,如要求数据满足平稳性、正态性和独立同分布等要求,而实际实测资料都具有非平稳性、相关性和异方差性,从而使得现有分析方法难以准确掌握测值规律,有时甚至得到错误结论。同时现有分析方法对敏感或重要测点分析确定往往无能为力,只能依靠专家经验。
实际上,堤坝健康监测的实际测值序列大多不平稳,尤其是堤坝危险期或非常工况期。而敏感或重要测点的确定往往同结构性态分析是相互关联的,是一个问题的两个方面。目前常规分析方法许多都只能适应平稳时间序列或对数据噪声敏感,从而既不能识别敏感测点,也不能对结构性态进行有效分析。
发明内容
本发明的主要目的在于,克服现有技术中的不足,充分利用复杂网络分析方法的多尺度分析功能和对误差不敏感的优势,提供一种基于复杂网络的堤坝健康监测敏感测点分析方法,克服常规分析方法对数据误差的敏感性,可发现各个测点的重要程度,实现各个测点之间相关关系的分析,并准确地掌握堤坝的工程健康性态,实现堤坝性态的有效预警。
为了达到上述目的,本发明所采用的技术方案是:
基于复杂网络的堤坝健康监测敏感测点分析方法,包括以下步骤:
1)根据堤坝健康监测项目和测点,应用布设于测点的传感器获取堤坝对应测点的实际测值时间序列;
2)运用等概率粗粒化的方法,将测值时间序列转化成符号序列;
3)把符号序列编码在复杂网络的拓扑结构中,构造复杂网路;
4)对复杂网络的特征参数进行分析,获取该堤坝健康监测项目的敏感测点及其性态演化规律,进行堤坝的工程健康性态判别。
本发明进一步设置为:所述步骤1)中堤坝健康监测项目包括环境量、渗流场、温度场和变形场。
本发明进一步设置为:所述步骤2)中等概率粗粒化的方法,具体包括步骤:2-1)设xmax和xmin分别是测值时间序列的最大值和最小值,转化后的符号共有N种,记为s1,s2,…,sN,定义
Si=Sj,xmin+(j-1)d≤xi<xmin+jd,j=1,…,N (1)
其中,N>2,d=(xmax-xmin)/N;由式(1)得到N种不同字符构成的符号序列;
2-2)记堤坝测值时间序列为{x(l),l=1,…,N},通过式(2)计算测值时间序列波动k(l),
其中,Δl为时间间隔;
2-3)通过式(3)计算不同波动值可能出现的概率p(k),
其中,Num(x)对应测值序列的波动模态x发生的次数;
2-4)把测值时间序列的波动k(l)分为5个区间,定义5个特征字符为式(4),
其中,t代表测值数值快速增加,r代表测值数值缓慢的增加,e代表测值数值不增不减,d代表测值数值缓慢的降低,f代表测值数值快速的降低;
从而将测值时间序列转换成符号序列为式(5),
S=(S1S2S3…),Si∈(t,r,e,d,f) (5)。
本发明进一步设置为:所述步骤3)中构造复杂网路,是在把符号序列编码在有向加权的复杂网络的拓扑结构之后,通过excle2pajek软件转换成pajek软件可处理的格式。
本发明进一步设置为:所述步骤4)中对复杂网络的特征参数进行分析,通过pajek软件进行。
本发明进一步设置为:所述步骤4)中复杂网络的特征参数包括平均最短路径长度、平均集聚系数和模块度系数。
本发明进一步设置为:所述环境量包括溃坝前上游水位测点、溃坝前下游水位测点和大气压力测点;所述渗流场包括渗压水位测点、渗流量测点和坝体水分测点;所述温度场采用分布式光纤测量,包括坝体表面气温测点、库水温测点、坝体内部温度测点、上游水温测点和下游水温测点;所述变形场包括坝体表面变形测点、坝体内部变形测点和倾度测点。
与现有技术相比,本发明具有的有益效果是:
根据溃坝试验实测数据,从复杂网络的角度揭示测值时间序列变化的动力学特征,首先运用等概率粗粒化的方法,将测值时间序列转化成符号序列,把符号序映射成网络参数构造复杂网路,然后通过对复杂网络的特征参数的分析,获取重要测点及其性态演化规律,从而实现堤坝健康性态分析,具有有效识别重要测点和判别堤坝性态演化的能力。其中,采用粗粒化方法,有利于克服数据误差对分析结果的影响,有利于实现不同尺度的分析;同时采用复杂网络分析方法,可以发现各个测点的重要程度,实现各个测点之间相关关系的分析,从而提取敏感测点,及时准确地掌握堤坝的工程健康性态,实现堤坝性态的更有效预警。
上述内容仅是本发明技术方案的概述,为了更清楚的了解本发明的技术手段,下面结合附图对本发明作进一步的描述。
附图说明
图1为本发明基于复杂网络的堤坝健康监测敏感测点分析方法的流程图;
图2为本发明实施例的堤坝试验过程中蓄水期的环境温度过程线;
图3为本发明实施例的堤坝试验过程中溃决过程的水位过程线;
图4为本发明实施例的溃坝前上、下游水位测值的复杂网络图;
图5为本发明实施例的整个过程渗流压力测点M1、M2对应的复杂网络图;
图6为本发明实施例的溃坝期渗流压力测点M1、M2对应的复杂网络图;
图7为本发明实施例的4个不同测点的光纤温度测值的复杂网络图;
图8为本发明实施例的溃坝期测点S1、S2测值对应的复杂网络图;
图9为本发明实施例的整个过程倾度测点K1、K2测值对应的复杂网络图;
图10为本发明实施例的溃坝期倾度测点K1、K2测值对应的复杂网络图;
图11为本发明实施例的整个过程的三组温度测值对应的复杂网络图。
具体实施方式
下面结合说明书附图,对本发明作进一步的说明。
如图1所示,本发明提供一种基于复杂网络的堤坝健康监测敏感测点分析方法,包括以下步骤:
1)根据堤坝健康监测所设置的项目应用监测传感器,获取该堤坝健康监测项目的测点的测值时间序列;其中,坝健康监测项目根据需要包括环境量、渗流场、温度场和变形场等。
2)运用等概率粗粒化的方法,将测值时间序列转化成符号序列;其中,等概率粗粒化的方法,具体包括步骤:
2-1)设xmax和xmin分别是测值时间序列的最大值和最小值,转化后的符号共有N种,记为s1,s2,…,sN,定义
Si=Sj,xmin+(j-1)d≤xi<xmin+jd,j=1,…,N (1)
其中,N>2,d=(xmax-xmin)/N;由式(1)得到N种不同字符构成的符号序列;
2-2)记堤坝测值时间序列为{x(l),l=1,…,N},通过式(2)计算测值时间序列波动k(l),
其中,Δl为时间间隔;
2-3)通过式(3)计算不同波动值可能出现的概率p(k),
其中,Num(x)对应测值序列的波动模态x发生的次数;
2-4)把测值时间序列的波动k(l)分为5个区间,定义5个特征字符为式(4),
其中,t代表测值数值快速增加,r代表测值数值缓慢的增加,e代表测值数值不增不减,d代表测值数值缓慢的降低,f代表测值数值快速的降低;
从而将测值时间序列转换成符号序列为式(5),
S=(S1S2S3…),Si∈(t,r,e,d,f) (5)。
3)把符号序列编码在有向加权的复杂网络的拓扑结构中,之后通过excle2pajek软件转换成pajek软件可处理的格式,构造复杂网路。
4)通过pajek软件对复杂网络的特征参数进行分析,获取该堤坝健康监测项目的敏感测点及其性态演化规律,进行堤坝的工程健康性态判别。
其中,复杂网络的特征参数包括平均最短路径长度、平均集聚系数和模块度系数。平均最短路径长度代表了一种模态向另一种模态转换所需要的时间,因此更短的平均最短路径长度可以表现出短程的相关性;而平均集聚系数则表示了波动模式之间转换的速率,更大的平均集聚系数表示转换的更频繁;同时模块度的数值表征了网络内部群簇相互之间联系的稀疏程度,更大的模块度表示具有更强的社团关系。
将每一个N字串视为网络的一个顶点,符号序列中相邻的两个N字串之间存在着一条边,由前一个N字串指向后一个N字串。如果两个顶点i和j之间存在着多条边,那么边数就是这两个顶点之间的边的权重Wij。假设网络分成了m个社团,构造一个m×m阶的矩阵c,cij表示连接两个社团i和社团j的边在复杂网络中所占有的比例,则模块度可以表示为:
其中,0≤M<1,M的值越接近1,意味着分割算法越好。
本发明实施例:
试验坝最大库容达10×104m3,拦河坝为粘性均质土坝,最大坝高达9.7m、坝顶长120m、宽3m,溢洪道位于大坝右侧、宽5m,左侧有孔径0.3m的放水涵洞。
堤坝健康监测项目包括环境量、渗流场、温度场和变形场,环境量包括溃坝前上游水位测点、溃坝前下游水位测点和大气压力测点,渗流场包括渗压水位测点、渗流量测点和坝体水分测点,温度场包括坝体表面气温测点、库水温测点、坝体内部温度测点、上游水温测点和下游水温测点,变形场包括坝体表面变形测点、坝体内部变形测点和倾度测点。
试验堤坝的试验过程包括大坝填筑期、首次蓄水期、水位稳定期、第二次蓄水期和溃决期等过程,工程历时46天,试验过程历时统计表如表1所示,试验过程中蓄水期的环境温度过程线如图2所示,试验过程中溃决过程的水位过程线如图3所示。
表1
其中,在大坝填筑期进行仪器埋设,共布设溃坝前上游水位测点1个、溃坝前下游水位测点1个、坝体表面气温测点1个、库水温测点1个、坝体内部温度测点1个、上游水温测点1个、下游水温测点1个、大气压力测点1个、渗压水位测点13个、渗流量测点1个、坝体水分测点2个、坝体表面变形测点6个、坝体内部变形测点9个;以及温度场采用分布式光纤测量,空间分辨率1m,坝体表面气温测点的温度测线20m,库水温测点的温度测线20m,坝体内部温度测点的温度测线343.1m。
根据堤坝健康监测项目,在试验堤坝的试验过程中分别获取其测值时间序列,将测值时间序列转换成符号序列,如2012/10/6 13:00:08到2012/10/620:00:07的坝内渗压M1测点的实测数据为43.25、43.25、43.25、43.25、43.25、43.26、43.26、43.26、43.26、43.26、43.28、43.28、43.28、43.31、43.3、43.29、43.29、43.29、43.22、43.22、43.22、43.22、43.22、43.13、43.13、43.13、43.13、43.13、43.14,则转换成的符号序列为:feeeereeeeteetddeefeeeefeeeer。
之后进行复杂网路构建和分析,获取不同测点中的敏感测点及其性态演化规律,进行堤坝的工程健康性态判别;以下根据不同监测项目的测点进行分别详述。
1、溃坝前上游水位测点和溃坝前下游水位测点的分析
将溃坝前后数据分开构建复杂网路,得出溃坝前上、下游水位监测数据的复杂网络图,如图4所示,图4中(a)为溃坝前上游水位测值的复杂网络图,图4中(b)为溃坝前下游水位测值的复杂网络图。
如图4-(a)所示的上游水位测值对应的复杂网络中,比较粗的线条有:ee→ee,ee→re,er→ee,re→ee,这些变化过程占了总数的69.3%,时间上这些变化出现在10月16日到11月16日,与蓄水期水位上升的实际情况相符。如图4-(b)所示的下游水位测值的复杂网络中,比较粗的线对应的边为:ee→ee,re→ee,ee→er,ee→re,这些变化占了变化总数的48.6%,时间上这些变化出现在10月14日到11月14日的蓄水期。
通过字频的统计发现r出现的频率都比较高。从以上统计可以看出,上游水位数据比下游水位数据的包含更多的信息。
对两组数据节点度进行统计,见表2所示的溃坝前上、下游水位测值节点度统计表。
表2
对表2中这些节点进行统计发现,上游水位主要的6个节点度占总数的98.85%,而下游水位主要的6个节点度占节点总数的92.58%,可见上游水位测值在某几个区间变化更为频繁。
溃坝期上游水位是不断降低的,而下游水位则在溃坝期瞬间水位有所上升。由于溃坝期水位数据相对较少,很难绘制成复杂网络图。对上、下游水位溃坝前数据的拓扑性质进行分析,如下表3所示的上、下游水位溃坝前数据的复杂网络指标。
表3
通过表3可以看出,上游水位测值对应的复杂网络在三个方面均大于下游测点,可见上游水位测点包含更多的信息,所以溃坝前上游水位测点为典型测点,应对其进行重点监测。
2、渗流压力的分析
渗压水位测点M1、M2在整个溃坝试验过程测值对应的复杂网络见图5,图5中(a)为整个过程渗压水位测点M1测值的复杂网络图,图5中(b)为整个过程渗压水位测点M2测值的复杂网络图。
首先对节点度进行分析,一个节点的度越大意味着这个节点在某种意义上就更为重要。表4、表5分别列举出了M1、M2的节点度大小的排序。
表4
表5
由表4和表5可以看出,节点er,re,ed,rd,de,dr的节点度数比较大,这说明这些节点所代表的波动模态在堤坝性态变化过程中起到了重要的关联作用,各种波动模态向这几种波动模态转变或者说被这几种波动模态转换的频率高。另外分别对发生这几种模态的时间段进行整理发现,M1的变化模态多集中在11月13日到11月17日,也就是第二次蓄水的末期和溃决期期间,明显的几种变化有:re→ee,ee→er,ee→rd,这些变化占到了变化总数的21.7%,而M2的变化模态多集中在11月15日到11月17日,比前者相应的要靠后,明显的变化有:ee→er,re→ee,de→ee,这些变化占变化总数的20.3%。通过节点的变化模态可以发现,M1的变化具有更早更集中的性质,但是M2节点的统计与M1亦相差不大,说明M1、M2在反应堤坝性态变化方面具有相同的特征,但前期M1作为坝体结构异常预警更加敏感。
渗压水位测点M1、M2在溃坝期测值对应的复杂网络见图6,图6中(a)为溃坝期渗压水位测点M1测值的复杂网络图,图6中(b)为溃坝期渗压水位测点M2测值的复杂网络图。
溃坝期与整个过程对应的复杂网络指标比较见表6,表6为渗流压力数据拓扑性质复杂网络特征值对比表。
表6
由表6可以看出,M2的复杂网络具有较大的平均集聚系数和较小的平均最短路径长度。平均最短路径长度表示了由一种模态向另一种模态转换所需要的时间,这个值接近于1,说明该模式的转换时间大约为2天左右,表现为了短程的相关性,这对堤坝健康测值的短期预测具有一定的参考价值。而较大的平均集聚系数则表明各种波动模态之间存在的小规模的群簇,而且群簇内部的关键较好,表现为波动的模式转换较为频繁。根据式(6)求得该网络模块度为0.10476,接近于0,因此该网络内部的群簇关系表现得较为稀疏,认为具有较弱的社团结构。并且M2测点溃坝期的数据具有更短的平均最短路径长度和较大的平均集聚系数,相比较M1测点的数据则不具有这种代表性。根据对这几种拓扑性质的分析,反映了M2测点测值序列波动的频繁,这些数据的变化反应出的就是各种波动模态之间的转换,可以更好的把握测值序列的变动,因而溃坝期应该对M2测点进行重点监测。
3、坝内温度的分析
如图7所示为4个不同测点的光纤温度测值的复杂网络图,图7中(a)为TF-1测点对应的复杂网络图,图7中(b)为TF-2测点对应的复杂网络图,图7中(c)为TF-3测点对应的复杂网络图,图7中(d)为TF-4测点对应的复杂网络图。
通过4个光纤温度测值对应的复杂网路的比较发现测点TF1在这些典型的模态变化中占有更好的比重,为:41.5%。另外根据统计发现rd→dr,dr→rd,rd→rr,rr→rd,4种变化出现的频率最高,体现在复杂网络图中就是所对应的线条也就是复杂网络的边更粗一些,并且这4种变化形式中测点TF-1分别占有40.22%,40.48%,42.59%,42.59%。对4个测点分别构建出图7所示的复杂网络图,然后进行拓扑性质的分析后与总体数据进行对比,见表7所示的光纤温度各测点测值复杂网络指标比较。
TF-1 | TF-2 | TF-3 | TF-4 | |
平均最短路径长度 | 1.00708 | 1.04545 | 1.04444 | 1.09091 |
平均集聚系数 | 0.98291 | 0.955222 | 0.96111 | 0.92015 |
模块度系数 | 0.00193 | 0.00165 | 0.00265 | 0.00638 |
表7
通过图7和表7可以发现,测点TF-1是变化最为集中的点,因此测点TF-1可以作为光纤测温的典型测点。下面对光纤温度数据的拓扑性质进行简要分析,测点TF-1数据的平均最短路径长度为1.00708,平均集聚系数为0.98291。由此可见,测点TF-1的复杂网络具有较大的平均集聚系数和较小的平均最短路径长度。平均最短路径长度表示了由一种模态向另一种模态转换所需要的时间,这个值接近于1,说明该模式的转换时间大约为2天左右,表现为了短程的相关性,而较大的平均集聚系数则表现了数据波动的频繁。根据式(6)求得汇总数据中该网络模块度为0.00168,具有相对很小的模块度,因此认为4个测点之间不存在明显的社团结构,群簇性表现得比较稀疏。从而根据拓扑性质的分析,得出TF-1为典型测点,建议在堤坝监测期间进行重点的监测。
4、含水率的分析
如图8所示为溃坝期测点S1、S2对应的复杂网络图,图8中(a)为溃坝期测点S1对应的复杂网络图,图8中(b)为溃坝期测点S2对应的复杂网络图。对图8所示的节点进行节点度的排序如表8和表9所示。
表8
表9
由表8和表9可以看出,节点rr,er,re,de,rd等节点的节点度比较大,说明这些节点代表的波动模态在模拟坝的溃坝的一系列数据变动之中起到了重要的关联作用,各种模态向这几种模态转换的频率高。通过图像可以直观的看出主要的变化形式有:ee→ee,rr→rr,ee→er,re→ee,这些变化占了变化总数的44.5%,另外对这些变化模态发生的时间段进行分析后发现,其中代表无变化的ee→ee中S1的变化模态多集中在10月15日到11月16日,对应的就是蓄水期和溃决期,相比较S2多集中在10月14日到11月16日,与S1差距不大。而代表连续缓慢上升的rr→rr变化形式,S1则多集中在10月30日到11月10日,而S2则集中在11月5日到11月13日,在宏观上差别亦不是很明显。将S1与S2两组数据分别绘制成复杂网络图8进行分析,并对比如表10所示的测点S1、S2溃坝期数据的拓扑性质与总体数据的比较。
复杂网络指标 | 整体 | S1 | S2 |
平均最短路径长度 | 1.16667 | 1.03636 | 1.10256 |
平均集聚系数 | 0.95411 | 0.96522 | 0.90857 |
模块度系数 | 0.16779 | 0.00167 | 0.00631 |
表10
通过比较明显可以看出,S1复杂网络具有较大的平均集聚系数和较小的平均最短路径长度。由此可见该网络模态之间转换更频繁,表现为短程的相关性,同时网络中各种波动模态之间存在小规模的群簇,并且群簇的内部关联较好,表现为波动的模式转换更为频繁。同时由式(6)求的网络的模块度为0.00167,通过比较发现,该网络具有很小的模块度,因此内部群簇关系比较稀疏,具有较弱的社团结构。根据上诉几种拓扑性质的分析,得出测点S1可以更好的反应测值序列的变化以及波动模态之间的转换,因此测点S1为典型测点,应该在日常监测中进行重点监测。
5、倾度的分析
倾度测点K1、K2在整个过程测值对应的复杂网络图见图9所示,图9中(a)为整个过程倾度测点K1对应的复杂网络图,图9中(b)为整个过程倾度测点K2对应的复杂网络图。对图9所示的节点度进行统计,如表11和表12所示。
表11
表12
由表11和表12可以看到,倾度测点K1、K2之间节点度差别比较大,K1的节点度统计相对更加集中一些,而K2的节点度的统计想对比较分散一些。由复杂网络图明显可以得到倾度数据K1中明显的变化为:ee→ee,re→ee,er→ee,ee→re,这些变化占到了变化总数的38.0%,对这些变化发生的时间进行统计发现多集中在11月13日到11月17日,也就是第二次蓄水的末期和溃决期。而K2的变化模态中明显的几种变化为:ft→rd,rf→tr,rf→tf,这些变化仅占变化总数的1.2%。因此初步判断K1的变化模态更加集中,是重要测点。
倾度测点K1、K2在溃坝期的测值数据得到的复杂网络图如图10所示,图10中(a)为溃坝期倾度测点K1测值的复杂网络图,图10中(b)为溃坝期倾度测点K2测值的复杂网络图。溃坝期与整个过程对应的复杂网络指标比较见表13,表13为倾度数据拓扑性质复杂网络特征值对比表。
拓扑性质/测点及时期 | K1整个过程 | K2整个过程 | K1溃坝期 | K2溃坝期 |
平均最短路径长度 | 1.04070 | 1.12944 | 1.28788 | 1.35498 |
平均集聚系数 | 0.97138 | 0.92783 | 0.83836 | 0.76318 |
模块度系数 | 0.00150 | 0.00773 | 0.12245 | 0.03991 |
表13
由表13可以清楚地看出K1测点数据所形成的复杂网络具有较大的平均集聚系数和较小的平均最短路径长度。平均最短路径长度值接近于1,说明该模式的转换时间大约为2天左右,表现为了短程的相关性,而较大的平均集聚系数则说明各种波动模态之间存在着关系较为紧密的一些群簇,它们之间的模式转换更为频繁。根据式(6)求得该网络模块度为0.00150,接近于0,因此认为该网络内部各群簇相互之间的关系表现的较为稀疏,不具有社团结构。根据上述复杂网络指标可知,K1测点为典型测点,在日常的监测中应当进行重点的监测。
6、上游水温、下游水温与气温的分析
整个试验过程的三组温度测值时间序列对应的复杂网络如图11所示,图11中(a)为上游水温测值的复杂网络图,图11中(b)为下游水温测值的复杂网络图,图11中(c)为气温测值的复杂网络图。
上游水温测值对应的复杂网络中比较粗的线条对应的变化有:ee→ee,de→ee,ee→ed,ee→de,这些变化占了总数的72.0%,对字频进行统计发现e,d出现的频率比较高,可以认为温度整体上呈下降的趋势,通过对时间的统计发现这些变化多出现在:10月18日到11月14日期间,亦符合自然规律的变化。
下游水位的温度数据中,复杂网络中比较粗的线对应的边所代表的有向线段为:ee→ee,de→ee,ee→ed,ee→er,er→de,这些变化占了变化总数的71.8%,对字频进行统计发现e,d出现的频率比较高,可以认为温度整体上呈下降的趋势,通过对时间的统计发现这些变化多出现在:10月18日到11月14日期间,亦符合自然规律的变化。
而坝体的温度计数据中节点度比较高的节点有dd(308),ff(213),tt(207),ee(185),可以看出节点随机性比较强,没有相对统一的变化趋势,对节点的变化形式进行统计发现,ff→ff,dd→dd,dd→de等几种变化发生次数最多,但是这些变化仅占到总变化数的7.8%,因而可以认为该组数据均不具有代表性,该网络的随机性比较强。
下面对它们的拓扑性质分别进行分析,如表14所示为上、下游水温与气温复杂网络指标。
表14
通过表14中各组数据进行比较分析,很显然,上游温度数据更具有代表性,同时上游的温度数据具有更小的模块度系数,证明其内部各群簇之间的关联比较稀疏,因此可以认为该组数据不具备明显的社团结构。综合以上分析,上游温度的数据包含更多信息,在日常的监测中应该进行重点的监测。
本发明的创新点在于,时间序列特征的复杂网络分析方法采用粗粒化方法,有利于克服数据误差对分析结果的影响,有利于实现不同尺度的分析;采用复杂网络分析方法,可以发现各个测点的重要程度,有利于分析各个测点之间的相关关系,而且可以提取敏感测点,从而及时准确地掌握堤坝的工程健康性态,有利于实现堤坝性态更有效预警。
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (6)
1.基于复杂网络的堤坝健康监测敏感测点分析方法,其特征在于,包括以下步骤:
1)根据堤坝健康监测项目和测点,应用布设于测点的传感器获取堤坝对应测点的实际测值时间序列;
2)运用等概率粗粒化的方法,将测值时间序列转化成符号序列;
3)把符号序列编码在复杂网络的拓扑结构中,构造复杂网路;
4)对复杂网络的特征参数进行分析,获取该堤坝健康监测项目的敏感测点及其性态演化规律,进行堤坝的工程健康性态判别;
所述步骤3)中构造复杂网路,是在把符号序列编码在有向加权的复杂网络的拓扑结构之后,通过excle2pajek软件转换成pajek软件可处理的格式。
2.根据权利要求1所述的基于复杂网络的堤坝健康监测敏感测点分析方法,其特征在于:所述步骤1)中堤坝健康监测项目包括环境量、渗流场、温度场和变形场。
3.根据权利要求1所述的基于复杂网络的堤坝健康监测敏感测点分析方法,其特征在于:所述步骤2)中等概率粗粒化的方法,具体包括步骤:
2-1)设xmax和xmin分别是测值时间序列的最大值和最小值,转化后的符号共有N种,记为s1,s2,…,sN,定义
Si=Sj,xmin+(j-1)d≤xi<xmin+jd,j=1,…,N (1)
其中,N>2,d=(xmax-xmin)/N;由式(1)得到N种不同字符构成的符号序列;
2-2)记堤坝测值时间序列为{x(l),l=1,…,N},通过式(2)计算测值时间序列波动k(l),
其中,Δl为时间间隔;
2-3)通过式(3)计算不同波动值可能出现的概率p(k),
其中,Num(x)对应测值序列的波动模态x发生的次数;
2-4)把测值时间序列的波动k(l)分为5个区间,定义5个特征字符为式(4),
其中,t代表测值数值快速增加,r代表测值数值缓慢的增加,e代表测值数值不增不减,d代表测值数值缓慢的降低,f代表测值数值快速的降低;
从而将测值时间序列转换成符号序列为式(5),
S=(S1S2S3…),Si∈(t,r,e,d,f) (5)。
4.根据权利要求1所述的基于复杂网络的堤坝健康监测敏感测点分析方法,其特征在于:所述步骤4)中对复杂网络的特征参数进行分析,通过pajek软件进行。
5.根据权利要求1所述的基于复杂网络的堤坝健康监测敏感测点分析方法,其特征在于:所述步骤4)中复杂网络的特征参数包括平均最短路径长度、平均集聚系数和模块度系数。
6.根据权利要求2所述的基于复杂网络的堤坝健康监测敏感测点分析方法,其特征在于:所述环境量包括溃坝前上游水位测点、溃坝前下游水位测点和大气压力测点;
所述渗流场包括渗压水位测点、渗流量测点和坝体水分测点;
所述温度场采用分布式光纤测量,包括坝体表面气温测点、库水温测点、坝体内部温度测点、上游水温测点和下游水温测点;
所述变形场包括坝体表面变形测点、坝体内部变形测点和倾度测点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510691360.5A CN105404709B (zh) | 2015-10-22 | 2015-10-22 | 基于复杂网络的堤坝健康监测敏感测点分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510691360.5A CN105404709B (zh) | 2015-10-22 | 2015-10-22 | 基于复杂网络的堤坝健康监测敏感测点分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105404709A CN105404709A (zh) | 2016-03-16 |
CN105404709B true CN105404709B (zh) | 2018-10-02 |
Family
ID=55470196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510691360.5A Active CN105404709B (zh) | 2015-10-22 | 2015-10-22 | 基于复杂网络的堤坝健康监测敏感测点分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105404709B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107944613A (zh) * | 2017-11-17 | 2018-04-20 | 湖南大学 | 一种基于符号网络的灾变预测方法 |
CN110069832B (zh) * | 2019-04-01 | 2021-07-27 | 南京航空航天大学 | 基于几何坐标转换的非连续板结构变形反演与拼接方法 |
CN111047079B (zh) * | 2019-11-25 | 2022-05-13 | 山东师范大学 | 一种风电场风速时间序列预测方法及*** |
CN113837674A (zh) * | 2021-11-29 | 2021-12-24 | 江西省水利科学院 | 一种基于权重模型的堤防险情敏感性分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104915555A (zh) * | 2015-05-28 | 2015-09-16 | 水利部南京水利水文自动化研究所 | 均质土坝溃坝分级预警指标的提取方法 |
-
2015
- 2015-10-22 CN CN201510691360.5A patent/CN105404709B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104915555A (zh) * | 2015-05-28 | 2015-09-16 | 水利部南京水利水文自动化研究所 | 均质土坝溃坝分级预警指标的提取方法 |
Non-Patent Citations (5)
Title |
---|
《利用复杂网络研究中国温度序列的拓扑性质》;周磊等;《物理学报》;20081115;第57卷(第11期);第7381页左栏第1段,第7381-7386页第2-3节 * |
《基于Pajek的航空网络特征评价与可生存性分析》;刘典勇等;《电子科技》;20151015;第28卷(第10期);第90-92页 * |
《基于等概率粗粒化的复杂度算法及其应用》;张佃中等;《***仿真学报》;20080805;第20卷(第15期);第4096-4098页 * |
《如何将excel形式数据转化为pajek的数据形式》;suigewu;《https://www.douban.com/note/200076634/》;20120212;第1页 * |
《径流时间序列复杂网络的动力学拓扑性质》;唐强等;《武汉纺织大学学报》;20140615;第27卷(第3期);第87-88页第1节 * |
Also Published As
Publication number | Publication date |
---|---|
CN105404709A (zh) | 2016-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105404709B (zh) | 基于复杂网络的堤坝健康监测敏感测点分析方法 | |
Howard et al. | Multivariate characterization of meandering | |
Sheu et al. | Extended maximum generally weighted moving average control chart for monitoring process mean and variability | |
CN104915555B (zh) | 均质土坝溃坝分级预警指标的提取方法 | |
CN110488218A (zh) | 一种电能表运行状态评估方法和评估装置 | |
CN112036042A (zh) | 一种基于变分模态分解的动力设备异常检测方法及*** | |
US20170249559A1 (en) | Apparatus and method for ensembles of kernel regression models | |
CN104035431A (zh) | 用于非线性过程监控的核函数参数的获取方法和*** | |
CN106951680A (zh) | 一种水文时间序列异常模式检测方法 | |
CN108469805A (zh) | 一种基于动态性最优选择的分散式动态过程监测方法 | |
Sahu et al. | Regional frequency analysis using L-moment methodology—a review | |
CN112213687B (zh) | 基于伪异常点辨识的关口电能表数据异常检测方法及*** | |
Qi et al. | Resilience changes in watershed systems: A new perspective to quantify long-term hydrological shifts under perturbations | |
Roushangar et al. | Exploring the multiscale changeability of precipitation using the entropy concept and self-organizing maps | |
CN112507591B (zh) | 一种基于分离式混合建模技术的重力坝坝顶位移监控方法 | |
CN114398786A (zh) | 压裂施工压力预测模型建立方法及装置、预测方法及装置 | |
CN112434887B (zh) | 一种结合网络核密度估计和svm的供水管网风险预测方法 | |
CN103389360A (zh) | 基于概率主元回归模型的脱丁烷塔丁烷含量软测量方法 | |
CN112801426A (zh) | 一种基于关联参数挖掘的工业过程故障融合预测方法 | |
Zidaoui et al. | Validation of wastewater data using artificial intelligence tools and the evaluation of their performance regarding annotator agreement | |
Mazaev et al. | Probabilistic leak localization in water distribution networks using a hybrid data-driven and model-based approach | |
Bai et al. | Evolving an information diffusion model using a genetic algorithm for monthly river discharge time series interpolation and forecasting | |
CN114818927B (zh) | 数据驱动的设备腐蚀预测方法 | |
Durante et al. | Cluster analysis of time series via Kendall distribution | |
Teoh et al. | Optimal design of the double sampling X̅ chart based on median run length |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |