第４３卷第３期２０２３年３月大地测量与地球动力学ＪｏｕｒｎａｌｏｆＧｅｏｄｅｓｙａｎｄＧｅｏｄｙｎａｍｉｃｓＶｏｌ．４３Ｎｏ．３Ｍａｒ．，２０２３收稿日期：２０２２－０４－２６项目来源：国家重点研发计划（２０１８ＹＦＣ１５０３８０３）；国家自然科学基金（４１７７２２５６）。第一作者简介：何案华，副研究员，主要从事地震前兆观测技术与理论研究，Ｅ－ｍａｉｌ：ｄｑｓ＿ｈａｈ＠１６３．ｃｏｍ。通讯作者：王言章，教授，博士生导师，主要从事弱磁信号检测、电磁探测传感器与数字信号处理研究，Ｅ－ｍａｉｌ：ｙａｎｚｈａｎｇ＠ｊｌｕ．ｅｄｕ．ｃｎ。ＤＯＩ：１０．１４０７５／ｊ．ｊｇｇ．２０２３．０３．０１７文章编号：１６７１－５９４２（２０２３）０３－０３１３－０５ＭＡＴＬＡＢ在地下流体数据分析中的基础应用（Ⅱ）何案华１，２，３王言章１，２１吉林大学地球信息探测仪器教育部重点实验室，长春市民主大街９３８号，１３００２６２吉林大学仪器科学与电气工程学院，长春市民主大街９３８号，１３００２６３应急管理部国家自然灾害防治研究院，北京市安宁庄路１号，１０００８５摘要：ＭＡＴＬＡＢ可用于地震地下流体数据求均值、去趋势性变化、频谱分析、自定义函数拟合、井孔水温梯度曲线绘制等操作。基于以上应用，本文进行井水位潮汐因子与气压因子计算，并讨论井水位中潮汐与气压成分的剔除问题。关键词：井水位；潮汐效应；气压效应；重力潮汐理论值；ＭＡＴＬＡＢ中图分类号：Ｐ３１５文献标识码：Ａ井水位观测值动态表达式一般为：yｎ＝ｘｎ＋Ｐｎ＋Ｔｎ＋Ｒｎ＋εｎ（１）式中，yｎ为水位观测值；ｘｎ为水位残差；Ｐｎ为气压效应；Ｔｎ为潮汐效应（体应变潮汐理论值）；Ｒｎ为降雨效应；εｎ为测量噪声，假定其为均值为０、方差为σ２的高斯白噪声［１］。为简化计算，暂不考虑井水位降雨效应。从已有研究来看，井水位中潮汐效应与气压效应并非是简单的一元线性关系，而是存在延时、滞后且持续一段时间。考虑实际过程，本文侧重讨论一元线性关系的计算，并对一元线性拟合后的数学问题进行简单讨论。１井水位潮汐因子与气压因子计算井水位潮汐与气压效应计算过程如图１所示，主要步骤如下：１）源数据预处理。①剔除明显错误数据；②缺数处理，由于ＭＡＴＬＡＢ许多计算过程要求数据完整，故采用一定数学方法对缺数进行插值；③降采样率，由于观测数据采样率多为１次／ｍｉｎ，甚至１次／ｓ，而地...
