CN115270608A - 基于arima与lstm的海岸带地面沉降预测方法 - Google Patents
基于arima与lstm的海岸带地面沉降预测方法 Download PDFInfo
- Publication number
- CN115270608A CN115270608A CN202210802633.9A CN202210802633A CN115270608A CN 115270608 A CN115270608 A CN 115270608A CN 202210802633 A CN202210802633 A CN 202210802633A CN 115270608 A CN115270608 A CN 115270608A
- Authority
- CN
- China
- Prior art keywords
- time sequence
- arima
- lstm
- time
- sequence
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 title claims abstract description 49
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000012549 training Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000007476 Maximum Likelihood Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 3
- 230000004913 activation Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 2
- 238000013135 deep learning Methods 0.000 abstract description 6
- 210000004027 cell Anatomy 0.000 description 14
- 238000012544 monitoring process Methods 0.000 description 6
- 230000007704 transition Effects 0.000 description 6
- 238000013459 approach Methods 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 210000002569 neuron Anatomy 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000010801 machine learning Methods 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000015654 memory Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000306 recurrent effect Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- 239000004106 carminic acid Substances 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001373 regressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000006403 short-term memory Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Geometry (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于ARIMA与LSTM的海岸带地面沉降预测方法,涉及深度学习领域。本发明包括以下步骤:InSAR的SLC影像数据经过预处理得到干涉图;根据干涉图利用StaMPS得到第一时间序列;基于第一时间序列,利用ARIMA模型进行初步的预测,得到预测时间序列;所述预测时间序列与所述第一时间序列做差,得到第二时间序列;利用LSTM模型对第二时间序列进行学习,得到第三时间序列;将第三时间序列与预测时间序列相加,得到原始InSAR形变量时间序列的预测结果。本发明相对单一的预测方法得到的精度指标有所提升。
Description
技术领域
本发明涉及深度学习领域,更具体的说是涉及一种基于ARIMA与LSTM 的海岸带地面沉降预测方法。
背景技术
海陆过渡带是陆地与海洋之间的交界地带,其中包含了复杂的物理过程、化学过程、生物过程以及地质过程等,因而它成为了地球表面动态环境最复杂多变的区域之一。目前海陆过渡带区域有数十亿人生活着,是世界上最繁荣发达的区域之一,其演变过程与人类的可持续发展直接相关联,因而该区域也受到全球各个沿海国家的高度关注。
地面沉降是指地下支撑物的移动导致地面标高损失的一种地质灾害,它是自然因素与人为因素共同影响的地质过程。其中自然因素包括构造下沉、地震及火山活动、气候变化等,而人为因素包括开采地下水、油气资源、矿产资源等。海陆过渡带区域的长期地面沉降的累积,会导致洪涝灾害的加剧、建筑物寿命的缩短以及交通网络的损失等,这些都会严重威胁到人们的生命财产安全。因而能否有效监测地面沉降并给出相应的预警信息,是十分必要的。
地面沉降监测方法主要有水准测量、GPS(Global Positioning System,全球定位***)测量以及InSAR(Interferometric Synthetic Aperture Radar,合成孔径干涉雷达)测量等。其中InSAR技术具有全天时、全天候、高精度、实时性、连续性等优点,广泛应用于陆地地面沉降的监测。
而对于地面沉降预测方法,主要有基于物理机制的方法、基于数理统计的方法以及基于机器学习的方法。其中基于物理机制的方法需要大量的实地监测数据,提出基于岩石、水文特性演化的物理模型,并对大量复杂参数进行求解,并且受到时效性的严重影响;基于数理统计的方法则是根据大量历史测量数据得到统计规律,并进一步进行分析;基于机器学习的方法则是从更高的数据维度提取现有历史数据的特征,进而得到对现有数据的反馈,不会受限于物理参数的影响。
上述预测方法,其本质是对于时间序列的预测问题。而在时序预测方面,主要的模型包括:传统的ARIMA模型(Auto Regressive Integrated Moving Average model,差分整合移动平均自回归模型),机器学习的支持向量机、贝叶斯网络、矩阵分解、高斯过程、深度学习等。
目前,在对InSAR数据得到的地面沉降时间序列经过单一预测模型处理后,得到的预测效果并不理想,因此,如何解决这一问题是本领域技术人员亟需研究的。
发明内容
有鉴于此,本发明提供了一种混合了ARIMA和深度学习中的LSTM的预测方法,即对InSAR得到的形变量时间序列与ARIMA方法得到的预测时间序列作差,然后利用LSTM对得到的时间序列差值进行学习,并把LSTM 模型预测的结果加在ARIMA预测数据上,就得到了原始InSAR形变量时间序列的预测结果。
为了实现上述目的,本发明采用如下技术方案:
一种基于ARIMA与LSTM的海岸带地面沉降预测方法,包括以下步骤:
InSAR的SLC影像数据经过预处理得到干涉图;
根据干涉图利用StaMPS得到第一时间序列;
基于第一时间序列,利用ARIMA模型进行初步的预测,得到预测时间序列;
所述预测时间序列与所述第一时间序列做差,得到第二时间序列;
利用LSTM模型对第二时间序列进行学习,得到第三时间序列;
将第三时间序列与预测时间序列相加,得到原始InSAR形变量时间序列的预测结果。
可选的,干涉图的相位组成如式(1)所示:
可选的,根据干涉图利用StaMPS得到第一时间序列的具体步骤如下:
StaMPS利用三维解缠算法,在时间上计算每个PS点的相位差异;
设置参考点用最小二乘法在空间上进行解缠,得到第一时间序列。
可选的,ARIMA模型具体为:
其中,Xt是t时刻的随机变量,αi以及βi是权重系数,p以及q表示时间窗长度,εt则代表了t时刻的白噪声。
可选的,所述LSTM模型包括多个细胞单元,
其中,xt表示第t时刻该细胞单元的输入,ct表示t时刻该细胞单元的状态,ht表示该细胞单元的隐状态输出;
每个细胞单元中含有三个控制门,分别是遗忘门f、输入门i以及输出门 o;遗忘门决定上一时刻的细胞单元状态ct-1有多少成分保存到ct中;输入门决定当前时刻输入xt有多少成分保存到ct之中;输出门决定ct中有多少成分被传递到ht之中;故其前向计算方法如下所示:
ft=σ(Wf·[ht-1,xt]+bf)#(5)
it=σ(Wi·[ht-1,xt]+bi)#(6)
ct=ftct-1+ittanh(Wc·[ht-1,xt]+bc)#(7)
ot=σ(Wo·[ht-1,xt]+bo)#(8)
ht=ottanh(ct)#(9)
式中,W和b分别为对应的权重系数矩阵和偏置项;σ和tanh分别为 sigmoid函数和双曲正切激活函数。
可选的,所述第二时间序列求解过程如下:
利用StaMPS方法得到若干个PS点,每个PS点对应的形变量时间序列记作{Yt},其长度为T。
可选的,所述LSTM模型的损失函数求解过程如下:
采用ADAM优化算法,设定损失函数为绝对误差:
经由上述的技术方案可知,与现有技术相比,本发明公开提供了一种基于ARIMA与LSTM的海岸带地面沉降预测方法,相对单一的预测方法得到的精度指标有所提升。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明的整体流程示意图;
图2为本发明的LSTM隐藏层细胞结构图;
图3为本发明的主副影像时间基线与空间垂直基线关系图;
图4(a)-图4(b)为本发明年平均形变速率;
图5(a)为总形变量图;图5(b)年平均形变量分布直方图;
图6(a)-图6(b)选取的PS点;
图7(a)-图7(b)为时序预测结果对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种基于ARIMA与LSTM的海岸带地面沉降预测方法,是一种混合了ARIMA和深度学习中的LSTM(Long Short-Term Memory, 长短期记忆单元)的预测方法,即对InSAR得到的形变量时间序列与ARIMA 方法得到的预测时间序列作差,然后利用LSTM对得到的时间序列差值进行学习,并把LSTM模型预测的结果加在ARIMA预测数据上,就得到了原始 InSAR形变量时间序列的预测结果。具体如图1所示,包括以下步骤:
InSAR的SLC影像数据经过预处理得到干涉图;
根据干涉图利用StaMPS得到第一时间序列;
基于第一时间序列,利用ARIMA模型进行初步的预测,得到预测时间序列
预测时间序列与第一时间序列做差,得到第二时间序列;
利用LSTM模型对第二时间序列进行学习,得到第三时间序列;
将第三时间序列与预测时间序列相加,得到原始InSAR形变量时间序列的预测结果。
其中,InSAR的SLC(Single Look Complex,单视复数)影像数据经过欧洲航天局的SNAP(SentiNel Application Platform)软件进行预处理得到干涉图,从而进一步利用StaMPS(Stanford Method for Persistent Scatterers)得到特定点的时间序列。后续则利用ARIMA方法进行初步的预测,与原始序列作差,并利用深度学习对该差值序列进行学习,改进了ARIMA方法的误差项,在此基础上进行预测。
StaMPS方法是由Hooper发展而来的一种从合成孔径雷达采集的时间序列中提取地面形变量的算法软件包。相比于传统的PS-InSAR(Persistent ScattererInterferometric Synthetic Aperture Radar,永久散射体合成孔径干涉雷达)选取永久散射体点的方法,StaMPS方法的选取规则可以识别散射强度较低但稳定性较高的目标,因而可以改善在非城镇地区传统PS-InSAR方法选取 PS点不足的情况,更加适用于海陆过渡带区域。
干涉相位图的相位组成如式(1)所示:
StaMPS利用了三维解缠算法,在时间上计算每个PS点的相位差异,然后设置参考点用最小二乘法在空间上进行解缠[26],相比于传统的方法精度更高。
解缠后的相位通过时间上的高通滤波,空间上的低通滤波,得到大气延迟相位并减去,从而矫正了大气误差。最后经过相位反演得到形变信息。
经典的时间序列预测模型包括自回归模型(Auto Regressive,AR),移动平均模型(Moving Average,MA),以及两者相结合的自回归移动平均模型(Auto Regressive MovingAverage,ARMA)[8],进一步利用差分运算,就得到了 ARIMA模型。AR(p)模型如式(2)所示,MA(q)模型如式(3)所示,ARMA(p,q) 模型如式(4)所示:
其中Xt是t时刻的随机变量,αi以及βi是权重系数,p以及q表示时间窗长度,εt则代表了t时刻的白噪声。
上述三种模型是弱平稳的,适用于平稳时间序列预测,但涉及到非平稳时间序列,则需要结合差分运算得到平稳序列后再进行预测,即ARIMA模型。经过实验比对,实际上ARIMA模型的预测误差稳定性较好,因而可以用于后面与InSAR数据作差,进而可以使利用LSTM进行残差序列学习结果的收敛性得到保障。
循环神经网络(Recurrent Neural Network,RNN)是一种能够将当前的数据与之前的数据建立联系的一种深度神经网络。但由于非线性函数的存在,其中离现在数据较远的数据会受到梯度衰减的影响,进而造成梯度消失问题;但另一方面,如果离当前数据较近的数据记录过多的情况下,也会造成梯度***问题。
为解决RNN的潜在问题,LSTM模型被提出。它继承了RNN的记忆性,并增加了门控制单元对信息进行取舍,然后经过反向传播的梯度调整,从而完成模型的建立。其中以隐藏层最为重要,它由多个细胞单元组成,其细胞单元结构如图2所示。
其中xt表示第t时刻该细胞单元的输入,ct表示t时刻该细胞单元的状态, ht表示该细胞单元的隐状态输出。
每个细胞单元中含有三个控制门,分别是遗忘门f、输入门i以及输出门 o。其中,遗忘门决定了上一时刻的细胞单元状态ct-1有多少成分保存到ct中;输入门决定了当前时刻输入xt有多少成分保存到ct之中;输出门决定了ct中有多少成分被传递到ht之中。故其前向计算方法如下所示:
ft=σ(Wf·[ht-1,xt]+bf)#(5)
it=σ(Wi·[ht-1,xt]+bi)#(6)
ct=ftct-1+ittanh(Wc·[ht-1,xt]+bc)#(7)
ot=σ(Wo·[ht-1,xt]+bo)#(8)
ht=ottanh(ct)#(9)
式中,W和b分别为对应的权重系数矩阵和偏置项;σ和tanh分别为 sigmoid函数和双曲正切激活函数。
LSTM的模型训练过程与循环神经网络一致,误差项δ的反向传播也是包括两个方向:一个是沿时间的反向传播,计算每个时刻的误差项;一个是将误差项向上一层神经元传播。然后根据误差项计算每个权重的梯度,进而利用梯度优化算法更新权重完成训练。
利用LSTM对ARIMA的残差序列进行学习的预测方法:假设利用StaMPS 方法得到若干个PS点,每个PS点对应的形变量时间序列记作{Yt},其长度为T。
由最大似然估计得到ARIMA(p,d,q)的最优参数,进而得到了预测序列
然后利用LSTM模型对上述标准化残差序列进行训练。设测试集的长度为Ttest,训练用时间序列长度L。实践表明,单步预测的误差是最小的,因而对时间序列进行标注,使其成为监督数据时,只需要最后一个作为标签即可。故训练用的样本输入为(L-1)个。
LSTM的主要参数包括训练用时间序列长度L,神经网络层数K,隐藏层神经元个数N,训练次数epochs以及学习率η。可以采用网格搜索法得到最优参数,同时采用随训练次数增加逐步衰减的学习率。
确定LSTM参数后,就可以开始对其进行训练了。采用ADAM优化算法[29],设定损失函数为绝对误差:
式中,RMSE(Root Mean Squard Error)表示均方根误差,MAE(Mean AbsoluteError)表示平均绝对误差,MAP(Mean Absolute Practice)表示平均实际绝对量,μ表示平均预测精度。
本实施例选取杭州湾部分地区(30°N-30°30′N,120°45′E-120°55′E)为研究区域。由于水准数据采用的是南边的绍兴市上虞区数据,故通过相对稳定的嘉绍大桥沟通南北两岸,以验证北边的形变数据的准确度。
以2017年1月至2019年12月获取的Sentinel-1A影像的上升轨道,VV 极化方式的SLC影像为遥感数据源,利用PS-InSAR技术选取2018年7月16 日的影像作为主影像,共得到配准的主副影像干涉图83幅,其主副影像的时间基线以及空间垂直基线的关系如图3所示。
然后利用StaMPS方法筛选一共得到274926个高相干点(永久散射体),并进一步得到沿雷达视线方向的地形变化时间序列,其中每个点拥有83个形变记录,其年平均形变速率及其叠加在卫星图上的情况如图4(a)-图4(b)所示。从图中可以看到,沿杭州湾两岸处地形整体呈现相对下降的趋势,而处在更靠近陆地的地方则在整体上呈现相对上升的趋势。
对所有PS点的形变量数据、年平均形变量数据做统计分析,3年时间内,该地区的形变量在-71.98mm~90.25mm范围内,年平均形变量为 -23.65mm/year~10.93mm/year。图5(a)代表了所有PS点的形变量数据分布情况,可以看到其中的绝大多数形变量数据集中在-20mm~20mm之间。图5(b) 为年平均形变量直方图分布,可以看到平均形变速率大多集中在 -10mm/year~5mm/year的范围内。
针对地面沉降漏斗区域的监测情况,选取两处地面沉降较为明显的PS点,分别利用ARIMA、LSTM以及ARIMA-LSTM三种方法进行进一步预测以及对比。用深蓝色框标注了选取的两个点,并分别命名为A和B如图6(a)-图6(b) 所示。图7则分别表明了两个点各自对应的预测方法对比结果。其RMSE、 MAE以及μ指标结果如表1所示。
这里采用的是最大似然估计得到ARIMA的最优参数,网格搜索法得到 LSTM的最优参数。ARIMA方法的(p,d,q)参数最终选择为(5,1,0);LSTM方法则按照时间先后顺序划分训练集以及测试集,2017年和2018年的数据作为训练集,2019年作为测试集,epochs次数为100,隐藏神经元个数为10个; ARIMA-LSTM方法的以2019年的前2/3年数据作为训练集,后1/3年为测试集,epochs次数为50,隐藏层神经元个数为12。
从图7(a)-图7(b)中可以看到ARIMA和LSTM方法均会出现预测的结果偏后出现的情况,这源自于两种方法对过去信息的一种记忆方式。假设某时刻的实际值为yt,前一个时刻的实际值为yt-1,ARIMA方法在对t时刻的值进行预测的时候,往往与t时刻的前几个时刻的值有一定的关系,当信息数量不是很大的时候,前一个时刻的实际值yt-1所对应的权重较大,因而其预测结果会更偏向于yt-1,反映在图中的结果就会显示为滞后。同样的道理也适用于 LSTM方法,因为本质上两者都是利用函数模型去逼近。而ARIMA和LSTM 相结合的方法则会改善这种情况,并且在某些情况下会得到比较接近InSAR 原始形变量时间序列的预测点,例如B点在2019年的最后预测值。
表1 三种方法的精度指标对比
从上表中可以看到,ARIMA-LSTM方法相比于单一的ARIMA方法或者 LSTM方法有着明显的提升。RMSE至少减小了48.75%,MAE至少减小了 24.84%,μ至少提升了15.19%。借助于预测结果的可接受性,可以将LSTM 训练的模型合理外推至没有相应InSAR数据的未来年份,并及时对地面沉降结果进行适当预警。
本发明通过建立LSTM模型,对InSAR形变量时间序列与ARIMA预测形变量时间序列之间的差值序列进行预测,得到的均方根误差、平均绝对误差均比单一的LSTM、ARIMA模型预测的小,且平均预测精度更高。在数值上具体体现为RMSE至少减小了48.75%,MAE至少减小了24.84%,μ至少提升了15.19%。
而对于海陆过渡带地面沉降监测而言,及时发现出现的严重沉降区域是必要的。并且如果能够在未来一段时间内没有或者难以获取InSAR数据监测的情况下,合理的预测算法将具有重大的意义,它为潜在的地质灾害的预警提供相应的帮助。就目前的情况来看,通过ARIMA预测只能针对某个PS点。对大量的PS点而言,通过LSTM学习ARIMA与原始时间序列差值的方法高度依赖于ARIMA模型的建立,因而两者结合的方法有一定局限性。后续的工作可以寻找大范围预测中能够生成较为稳定的差值序列的方法,进而用LSTM 模型进行拟合,得到更加精确的预测结果。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (7)
1.一种基于ARIMA与LSTM的海岸带地面沉降预测方法,其特征在于,包括以下步骤:
InSAR的SLC影像数据经过预处理得到干涉图;
根据干涉图利用StaMPS得到第一时间序列;
基于第一时间序列,利用ARIMA模型进行初步的预测,得到预测时间序列;
所述预测时间序列与所述第一时间序列做差,得到第二时间序列;
利用LSTM模型对第二时间序列进行学习,得到第三时间序列;
将第三时间序列与预测时间序列相加,得到原始InSAR形变量时间序列的预测结果。
3.根据权利要求1所述的一种基于ARIMA与LSTM的海岸带地面沉降预测方法,其特征在于,根据干涉图利用StaMPS得到第一时间序列的具体步骤如下:
StaMPS利用三维解缠算法,在时间上计算每个PS点的相位差异;
设置参考点用最小二乘法在空间上进行解缠,得到第一时间序列。
5.根据权利要求1所述的一种基于ARIMA与LSTM的海岸带地面沉降预测方法,其特征在于,所述LSTM模型包括多个细胞单元,
其中,xt表示第t时刻该细胞单元的输入,ct表示t时刻该细胞单元的状态,ht表示该细胞单元的隐状态输出;
每个细胞单元中含有三个控制门,分别是遗忘门f、输入门i以及输出门o;遗忘门决定上一时刻的细胞单元状态ct-1有多少成分保存到ct中;输入门决定当前时刻输入xt有多少成分保存到ct之中;输出门决定ct中有多少成分被传递到ht之中;故其前向计算方法如下所示:
ft=σ(Wf·[ht-1,xt]+bf)#(5)
it=σ(Wi·[ht-1,xt]+bi)#(6)
ct=ftct-1+ittanh(Wc·[ht-1,xt]+bc)#(7)
ot=σ(Wo·[ht-1,xt]+bo)#(8)
ht=ottanh(ct)#(9)
式中,W和b分别为对应的权重系数矩阵和偏置项;σ和tanh分别为sigmoid函数和双曲正切激活函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210802633.9A CN115270608A (zh) | 2022-07-09 | 2022-07-09 | 基于arima与lstm的海岸带地面沉降预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210802633.9A CN115270608A (zh) | 2022-07-09 | 2022-07-09 | 基于arima与lstm的海岸带地面沉降预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115270608A true CN115270608A (zh) | 2022-11-01 |
Family
ID=83766022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210802633.9A Pending CN115270608A (zh) | 2022-07-09 | 2022-07-09 | 基于arima与lstm的海岸带地面沉降预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115270608A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115758089A (zh) * | 2022-11-08 | 2023-03-07 | 海南师范大学 | 一种钻孔应变数据预测方法 |
CN116068511A (zh) * | 2023-03-09 | 2023-05-05 | 成都理工大学 | 一种基于深度学习的InSAR大尺度***误差改正方法 |
-
2022
- 2022-07-09 CN CN202210802633.9A patent/CN115270608A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115758089A (zh) * | 2022-11-08 | 2023-03-07 | 海南师范大学 | 一种钻孔应变数据预测方法 |
CN116068511A (zh) * | 2023-03-09 | 2023-05-05 | 成都理工大学 | 一种基于深度学习的InSAR大尺度***误差改正方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11333796B2 (en) | Spatial autocorrelation machine learning-based downscaling method and system of satellite precipitation data | |
Bui et al. | Landslide susceptibility mapping at Hoa Binh province (Vietnam) using an adaptive neuro-fuzzy inference system and GIS | |
Mathew et al. | Weights of evidence modelling for landslide hazard zonation mapping in part of Bhagirathi valley, Uttarakhand | |
Li et al. | GIS‐Based Landslide Susceptibility Mapping Using Information, Frequency Ratio, and Artificial Neural Network Methods in Qinghai Province, Northwestern China | |
CN105243435B (zh) | 一种基于深度学习元胞自动机模型的土壤含水量预测方法 | |
Alkhasawneh et al. | Determination of important topographic factors for landslide mapping analysis using MLP network | |
Hu et al. | Satellite-based high-resolution mapping of ground-level PM2. 5 concentrations over East China using a spatiotemporal regression kriging model | |
CN115270608A (zh) | 基于arima与lstm的海岸带地面沉降预测方法 | |
CN110986747A (zh) | 一种滑坡位移组合预测方法及*** | |
CN102768361A (zh) | 基于遗传粒子滤波与模糊神经网络的gps/ins组合定位方法 | |
Muneer et al. | Modeling of runoff in the arid regions using remote sensing and geographic information system (GIS) | |
CN112446559A (zh) | 一种基于深度学习的大范围地面沉降时空预测方法和*** | |
Zandi et al. | Mapping of regional potential groundwater springs using logistic regression statistical method | |
Talaei | Landslide susceptibility zonation mapping using logistic regression and its validation in Hashtchin Region, northwest of Iran | |
CN114861840B (zh) | 一种多源降水数据融合方法 | |
CN105678417A (zh) | 一种施工隧道撑子面涌水量预测方法及装置 | |
Zhu et al. | Internal and external coupling of Gaussian mixture model and deep recurrent network for probabilistic drought forecasting | |
Yao et al. | Performance comparison of landslide susceptibility mapping under multiple machine-learning based models considering InSAR deformation: a case study of the upper Jinsha River | |
Wang et al. | Analysis and prediction of regional land subsidence with InSAR technology and machine learning algorithm | |
CN116910534A (zh) | 一种不同海域海洋环境要素时空智能预测方法和装置 | |
Chen et al. | Assimilation of SBAS-InSAR based vertical deformation into land surface model to improve the estimation of terrestrial water storage | |
Li et al. | An improved cyclic multi model-eXtreme gradient boosting (CMM-XGBoost) forecasting algorithm on the GNSS vertical time series | |
Shi et al. | Integrating SBAS-InSAR and LSTM for subsidence monitoring and prediction at Hong Kong international airport | |
CN109828270B (zh) | 一种表征地面沉降时序演变的方法 | |
Yue et al. | The spatio-temporal reconstruction of lake water levels using deep learning models: A case study on Altai mountains |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |