一种基于S-K特征的钻孔应变数据异常提取方法技术

本发明专利技术属于地震前兆观测数据异常检测领域，具体地来讲为一种基于S‑K特征的钻孔应变数据异常提取方法，首先是将同一台站的钻孔应变数据序列进行面应变换算，对换算后的数据进行调和分析处理；将处理后的面应变数据差分计算制成样本数据；并计算每一天的面应变样本数据的偏度和峰度；将得到的偏度和峰度画在一个平面内，并发现其呈现出的抛物线关系；基于这个关系，发现大多数天的平稳数据的偏度峰度都在零附近，故定义了一个地壳平稳时的背景并计算每一天与背景的偏移程度。通过本发明专利技术能够有效的对钻孔应变数据进行分析，对可能的地震前兆异常进行提取。

【技术实现步骤摘要】
一种基于S-K特征的钻孔应变数据异常提取方法
本专利技术属于地震前兆观测数据异常检测领域，具体地来讲为一种基于S-K特征的钻孔应变数据异常提取方法。
技术介绍
国的钻孔应变观测起步较早，特别是首创的水平四分量钻孔应变仪，设计理念具有国际领先水平。国家“十五”期间，这种四分量钻孔应变仪观测点发展较快，已经从原来屈指可数的几个观测点增加到50多个。其中一些观测点已经开始产出高质量的资料。中国地震局地震监测网络的首要目的是地震预测预报。地震预报研究的一个重要发展方向，是从“经验预报”向物理预报转变。连续地应变观测资料正是物理预报赖以进行的依据。地震预报的另一个发展方向，是解决短l临前兆问题。在数种主流大地测量观测方法(测震、GPS和钻孔应变观测)中，钻孔应变仪也正是在数月至数小时时段上观测精度最高。除此之外，钻孔应变观测在避免地面干扰、少占用地等方面也具有优势。CN106918836A公开了一种基于主成分分析的钻孔应变数据异常提取方法，首先是将同一台站的钻孔应变数据序列进行应变换算，对换算后的数据进行预处理；将预处理后的钻孔应变数据构造成一个矩阵；并对每一天的矩阵进行主成分分析，以获取每个矩阵的特征值和特征向量；将得到的特征值与计算出来的特征向量角度与地震事件相对应，以得到异常的检测结果。通过本专利技术能够有效的利用主成分分析的方法对钻孔应变数据进行分析，根据钻孔应变各测项的相关性，对可能的地震前兆异常进行提取。本专利技术钻孔应变数据异常提取方法，利用主成分分析中的特征值和特征向量角度分别将地壳的微弱变化表征出来；实现了在有较强背景干扰的情况下对钻孔应变数据异常的精确提取。根据长期的观测资料和经验，一般认为短周期的高频信息在无其他干扰因素的条件下，多表现为一种服从正态分布的随机信号。从统计的角度分析，在去除趋势变化后，且无其他影响因素存在的条件下，观测资料变化值服从正态分布。偏度和峰度就是统计学中描述数据分布偏离正态分布的两个参数，基于这一特性，本专利技术提出了一种提取钻孔应变数据中异常的方法。到目前为止，尚未见有关于偏度和峰度的方法对钻孔应变数据进行异常提取的报道。
技术实现思路
本专利技术所要解决的技术问题在于提供一种基于S-K特征的钻孔应变数据异常提取方法，能够有效的对钻孔应变数据进行分析，对可能的地震前兆异常进行提取。本专利技术是这样实现的，一种基于S-K特征的钻孔应变数据异常提取方法，包括如下步骤：a、对钻孔应变数据进行数据有效性的验证，若为有效，则进行下一步；b、对钻孔应变数据进行面应变换算与调和分析去低频信息；c、对处理后数据差分计算并制作样本数据；d、计算每天样本数据的偏度和峰度；e、将偏度和峰度画在S-K平面内；f、输入S-K平面散点图；g、定义一个地壳平静背景圆；h、计算每天样本数据偏离背景的程度，定义为S-K偏差；i、输出随时间变化的S-K偏差曲线。进一步地，步骤a所述的录入钻孔应变数据是选取一个台站的四分量钻孔应变数据，制作成按照分钟值采样的时间序列，数据有效性验证是按照不同分量，表示为S1,S2,S3,S4；根据由带孔平板弹性理论模型导出的圆孔径向变形与区域应力的关系式，对钻孔应变数据进行自洽分析。进一步地，步骤b所述对钻孔应变数据进行应变换算式根据下式(2)将四分量钻孔应变观测数据换算成面应变Sa：Sa(t)＝(S1(t)+S2(t)+S3(t)+S4(t))/2(2)。进一步地，步骤b中所述调和分析用来去除数据中固体潮低频信息的周期项,其调和分析函数S(t)表达式为：其中A0为时间序列的直流分量，m为谐波的次数，系数Am,Bm是权重因子，表示各次谐波对总序列的贡献。进一步地，差分表达式为：DSa(t)＝(Sa(t+1)-S(t+1))-(Sa(t)-S(t))，(4)，按天分段制成N天的样本数据DS(N)，其中DS(N)表示差分后的面应变。进一步地，计算每天样本数据的偏度和峰度，利用统计学上高阶累积量的算法计算，定义现有随机变量x，其均值是μ，方差是σ，则偏度和峰度表达式如下：其中，随机变量x取自样本数据，sk为偏度，ku为峰度，正态分布的偏度为0，若数据分布是对称的，偏度为0，若偏度>0,则分布为右偏；若偏度<0，则分布为左偏，同时偏度的绝对值越大，说明分布的偏移程度越严重。进一步地，定义一个地壳平静背景圆包括：S-K平面内所有的点中去掉一部分的极端点，剩余点求坐标参数的平均值作为圆点，坐标参数的标准差的平均值作为圆的半径。本专利技术与现有技术相比，有益效果在于：本专利技术公开的基于偏度和峰度的钻孔应变数据异常提取方法，第一次从数据分布上分析钻孔应变数据，并第一次把偏度峰度和数据的物理意义对应，深入的了解了钻孔应变数据；第一次发现钻孔应变数据的偏度和峰度成抛物线关系，并定义了合适背景，计算了每天面应变与背景的随时间的偏差，提出了有效的提取分布异常的方法。采用该方法能够对可能的地震前兆异常进行提取，根据以往发生地震前的数据分析，能否正确的提取地震前数据的异常。附图说明图1是基于S-K特征的钻孔应变数据异常提取方法流程图。图2是姑咱地震前兆监测台站与汶川地震震中位置示意图。图3是2007年1月1日到2009年7月1日的钻孔应变面应变的差分数据。图4是2007年1月1日到2009年7月1日的钻孔应变面应变的差分数据的每天的偏度(a)峰度(b)。图5是所有时间内的S-K平面图；图中灰色的点代表所有点，橙色的点是2007年1月1日到2009年7月1日的(S，K)，放大图中的灰色阴影圆是地壳平静背景；图6是2007年1月1日到2009年7月1日的S-K偏差曲线。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。实施例针对汶川地震，以四川地区地震前兆监测台站中姑咱台的钻孔应变数据为例。姑咱台站与汶川地震震中位置如图2所示。该数据由YRY四分量钻孔应变仪测得，采样一分钟一次，研究的时间段为2007年1月1日至2009年7月1日。参见图1所示：步骤一、录入钻孔应变数据,并进行数据有效性的验证，是，进行下一步，录入钻孔应变数据是选取一个台站的四分量钻孔应变数据，制作成按照分钟值采样的时间序列。数据有效性验证是按照不同分量，具体为：录入姑咱台站2007年1月1日至2013年7月1日钻孔应变分钟值时间序列，按照北南分量、东西分量、北东分量、北西分量的顺序将数据分别记为S1,S2,S3,S4；根据由带孔平板弹性理论模型导出的圆孔径向变形与区域应力的关系式，对钻孔应变数据进行自洽分析，其关系式为：S1+S3＝k(S2+S4)，(1)计算出的k值，选取自洽系数k在0.9以上的数据为有效数据，由k≥0.9可知姑咱台站2007年1月1日至2013年12月31日钻孔应变分钟值数据有效。步骤二、对钻孔应变数据进行面应变换算与调和分析处理：对钻孔应变数据进行应变换算是根据公式(2)将四分量钻孔应变观测数据换算成面应变Sa。Sa(t)＝(S1(t)+S2(t)+S3(t)+S4(t))/2，(2)然后，应用调和分析是用来去除数据中固体潮等低频信息的周期项。其调和分析函数S(t)表达式为：本文档来自技高网...

【技术保护点】
1.一种基于S‑K特征的钻孔应变数据异常提取方法，其特征在于，包括如下步骤：a、对钻孔应变数据进行数据有效性的验证，若为有效，则进行下一步；b、对钻孔应变数据进行面应变换算与调和分析去低频信息；c、对处理后数据差分计算并制作样本数据；d、计算每天样本数据的偏度和峰度；e、将偏度和峰度画在S‑K平面内；f、输入S‑K平面散点图；g、定义一个地壳平静背景圆；h、计算每天样本数据偏离背景的程度，定义为S‑K偏差；i、输出随时间变化的S‑K偏差曲线。

【技术特征摘要】
1.一种基于S-K特征的钻孔应变数据异常提取方法，其特征在于，包括如下步骤：a、对钻孔应变数据进行数据有效性的验证，若为有效，则进行下一步；b、对钻孔应变数据进行面应变换算与调和分析去低频信息；c、对处理后数据差分计算并制作样本数据；d、计算每天样本数据的偏度和峰度；e、将偏度和峰度画在S-K平面内；f、输入S-K平面散点图；g、定义一个地壳平静背景圆；h、计算每天样本数据偏离背景的程度，定义为S-K偏差；i、输出随时间变化的S-K偏差曲线。2.按照权利要求1所述的方法，其特征在于，步骤a所述的录入钻孔应变数据是选取一个台站的四分量钻孔应变数据，制作成按照分钟值采样的时间序列，数据有效性验证是按照不同分量，表示为S1,S2,S3,S4；根据由带孔平板弹性理论模型导出的圆孔径向变形与区域应力的关系式，对钻孔应变数据进行自洽分析。3.按照权利要求2所述的方法，其特征在于，步骤b所述对钻孔应变数据进行应变换算式根据下式(2)将四分量钻孔应变观测数据换算成面应变Sa：Sa(t)＝(S1(t)+S2(t)+S3(t)+S4(t))/2(2)。4.按照权利要求2所述的方法，其特征在于，...

【专利技术属性】
技术研发人员：朱凯光，于紫凝，池成全，樊蒙璇，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林,22