[发明专利]一种基于ICA-RA的钻孔应变数据的分频段级联校正方法

摘要

本发明涉及钻孔应变领域，具体地而言为一种基于ICA‑RA的钻孔应变数据的分频段级联校正方法，目的是去除外界环境干扰引起的应变变化，首先是将同一观测站点的钻孔应变数据、气温数据、气压数据、降水数据和计算的理论固体潮数据进行4层小波分解；在分解后的5个频段中，逐频段进行多因素与应变的相关性分析；利用ICA‑RA级联校正方法，在各个频段内找到与钻孔应变数据相关性最大的影响因素。引入独立成分分析(ICA)方法，计算由于水位变化引起的应变变化并消除；对于非是气温、气压、和固体潮这类数据，利用回归分析(RA)方法，计算对应频段内主要因素产生的影响并校正。通过本发明能够有效的对钻孔应变数据各频段进行分析，对真实的应变情况有着重要的校正作用。

主权项

1.一种基于ICA-RA的钻孔应变数据的分频段级联校正方法，其特征在于，该方法包括：/na、录入钻孔应变数据，录入同站点的气温、气压和水位数据共4组；/nb、应用滑动平均法计算步骤a录入的4组数据的小时值数据；/nc、计算同一站点位置的等时间跨度的理论固体潮数据；/nd、对步骤b的4组数据的小时值数据和步骤c的理论固体潮数据，共5组数据均进行4层小波分解并重构，将数据分解到5个不同频段内，其中5个不同频段包括低频段作为第一频段，以及细节频段作为第二频段；/ne、选取N＝1频段作为第一频段时的各因素作为初始状态；/nf、计算各因素与对应频段应变的相关系数，并按降序排列；/ng、判断是否存在一个相关系数大于0.8，是，进行步骤h；否，进行步骤p；/nh、判断应变是否受水位数据影响最大，是，下一步；否，跳至步骤j，判断气温；/ni、利用独立成分分析方法计算水位变化引起的应变变化并校正；/nj、判断应变是否受气温数据影响最大，是，下一步；否，跳至步骤l，判断气压；/nk、利用回归分析计算气温变化引起的应变变化并校正；/nl、判断应变是否受气压数据影响最大，是，下一步；否，跳至步骤n，判断固体潮；/nm、利用回归分析计算气压变化引起的应变变化并校正；/nn、判断应变是否受固体潮数据影响最大，是，下一步；否，跳至步骤p；/no、利用回归分析计算固体潮变化引起的应变变化并校正；/np、判断是否最后一层已经校正结束，是，下一步；否，频段数N加1进入细节频段并重复步骤f～o；/nq、叠加各频段内剩余的应变。/n