MATLAB在地下流体数据分析中的基础应用(Ⅱ)_何案华.pdf

第 卷 第期 年月大 地 测 量 与 地 球 动 力 学 ，收稿日期：项目来源：国家重点研发计划（）；国家自然科学基金（）。第一作者简介：何案华，副研究员，主要从事地震前兆观测技术与理论研究，：。通讯作者：王言章，教授，博士生导师，主要从事弱磁信号检测、电磁探测传感器与数字信号处理研究，：。：文章编号：（）在地下流体数据分析中的基础应用（）何案华，王言章，吉林大学地球信息探测仪器教育部重点实验室，长春市民主大街 号，吉林大学仪器科学与电气工程学院，长春市民主大街 号，应急管理部国家自然灾害防治研究院，北京市安宁庄路号，摘要：可用于地震地下流体数据求均值、去趋势性变化、频谱分析、自定义函数拟合、井孔水温梯度曲线绘制等操作。基于以上应用，本文进行井水位潮汐因子与气压因子计算，并讨论井水位中潮汐与气压成分的剔除问题。关键词：井水位；潮汐效应；气压效应；重力潮汐理论值；中图分类号：文献标识码：井水位观测值动态表达式一般为：y（）式中，y为水位观测值；为水位残差；为气压效应；为潮汐效应（体应变潮汐理论值）；为降雨效应；为测量噪声，假定其为均值为、方差为的高斯白噪声。为简化计算，暂不考虑井水位降雨效应。从已有研究来看，井水位中潮汐效应与气压效应并非是简单的一元线性关系，而是存在延时、滞后且持续一段时间。考虑实际过程，本文侧重讨论一元线性关系的计算，并对一元线性拟合后的数学问题进行简单讨论。井水位潮汐因子与气压因子计算井水位潮汐与气压效应计算过程如图所示，主要步骤如下：）源数据预处理。剔除明显错误数据；缺数处理，由于 许多计算过程要求数据完整，故采用一定数学方法对缺数进行插值；降采样率，由于观测数据采样率多为次 ，甚至次，而地球固体潮汐、气压等周期都是以小时为单位，为减少计算的复杂性，进行适当降采样率处理，本文采用次 。）潮汐理论值计算。常规计算方法为采用专用软件，如 按指定要求进行计算，得到相应数据文件后再参与计算，但这会增加处理的复杂性。通过梳理重力固体潮汐理论，采用 自定义函数进行实时计算，即根据台站基础信息（经纬度）以及数据时间进行计算。）潮汐因子计算。采用低通滤波提取趋势变化，然后用源数据减去低通滤波的趋势性变化即得到仅保留潮汐成分的水位，利用一元线性回归计算潮汐因子，进而剔除水位中的潮汐成分。）气压因子计算。由于气压效应存在滞后效应，故先求出最大相关系数，并取得滞后时间；再根据滞后时间进行一元线性回归计算，求取气压因子并剔除气压干扰。图井水位潮汐与气压效应计算过程 大 地 测 量 与 地 球 动 力 学 年月根据上述过程，代码如下：（）；截止频率 （）；采样频率 ，（，）；（）；（，）；（，）；（日期（年），）；（水位（），）；（水位原始观测数据，水位低通滤波后数据，）；（）；（，）；（，）；（日期（年），）；（去掉低频波成分后水位数据，）；（水位（），）；将水位原始数据通过取均值方法降采样率为次 ，低通滤波取截止频率（），图为低通滤波后产品。图井水位低通滤波及残差 ；按台站信息求出重力潮汐理论值 ：（）（），）；（）；（）；（，水位（），重力潮汐（），日期（年），（，），）；（）；（，）；（，（）；，（，）；：；（）（）；（，）；（，：）；（重力潮汐（），）；（水位（），）；（，（），（）；（，）；采用自定义的 计算重力潮汐理论值。通过 自带的 完成一元线性回归，可通过 查询回归结果，在本实例中得到回归结果 ，说明回归效果非常理想。求得潮汐因子为 ，负号是由于水位为静水位数据，其数值大小与重力潮汐呈负相关（图）。图井水位潮汐因子计算 （）；（）（）；剔除潮汐干扰 （：）；计算水位与气压最大相关系数，确定气压效应滞后时间 ；：（：）；（，）；（）；（）；（：），（：），（：第 卷第期何案华等：在地下流体数据分析中的基础应用（），水位（），气压（），日期（年），（，），）；（：）；（：）；（，（）；，（，）；（）；（，）；：；（）（）；（，）；（气压（），）；（水位（），）；（，（），（）；（，）；（）（）；（）；（：），）；，（：），（，），）；（，）；（，），汤加火山，）；采用查找最大相关系数方法，先求取气压效应的滞后时间，本例中气压效应滞后时间 ，采样率为次 ，即气压效应滞后时间为 。利 用 相 同 线 性 回 归 求 得 气 压 因 子 为 。通过去除气压效应，可清晰观察到原来隐藏在观测数据中的汤加火山爆发事件（图）。图井水位气压因子计算及最终井水位动态 重力潮汐理论值计算重力潮汐理论值的计算原理与计算过程可参考文献，但上述文献未给出详细的 实现，本文在此提供。（，）；（，）（）；，（）；（，）；（）；（，）（，）；（，）；（，）（，）；（，）（，）；（，）（，）；（，）（，）；大 地 测 量 与 地 球 动 力 学 年月 ；（）（）（）（）（）（）；（）（）（）（）（）（）（）（）（）（）（）（）（）；（）（）（）（）（）（）（）；（）（）；（）（）；（）；（）（）（）（）（）；（）（）（）（）（）（）（）（）（）；（）（）；（）（）（）（）（）；（）（）；（）（）；（）；（，）（，）（，）（，）（，）（，）；结语通过上述过程，可自动计算井水位潮汐因子、气压因子。通过剔除井水位潮汐干扰与气压干扰，可清晰地观测到隐藏在井水位原始观测数据中由汤加火山爆发引发的气压扰动事件。从图最终产品可知，井水位动态中仍然可清晰地看到潮汐成分，这与井水位潮汐和气压响应原理有关，两者之间并非简单的线性关系，存在时间滞后，且响应不是一个时间点，而是一个作用过程。另外，气压动态自身包含潮汐作用力，因此将潮汐与气压分两步进行剔除会引入气压分量中的潮汐作用力，这在上述计算过程中无法避免。要解决该问题，需采用多元线性、偏最小二乘法、卷积等其他回归算法，在此不过多讨论，详细过程可参考文献。参考文献刘阳，何案华，赵刚，等地下水位中潮汐与气压效应分析大地测量与地球动力学，（）：（，（）：），（）：何案华，王言章 在地下流体数据分析中的基础应用大地测量与地球动力学，（）：（，（）：）董良，彭芳苹，杨涛，等利用新参数和软件改进重力固体潮计算程序地球物理学进展，（）：（，（）：）闫如玉，万永革，解朝娣，等基于 的全球重力固体潮可视化实现地震地磁观测与研究，（）：（，（）：）第 卷第期何案华等：在地下流体数据分析中的基础应用（）（），：，：；：，；，：，：，：（上接第 页）y ，：，：），；），：；：，；，；，：，：，：