兴济地震台大地电场方位变化特征.pdf

李瑞卿，郭学增，李凤，等.兴济地震台大地电场方位变化特征 J.华北地震科学，2023,41（3）：91-96.doi:10.3969/j.issn.10031375.2023.03.013.LI Ruiqing，GUO Xuezeng，LI Feng，et al.Azimuth Variation Characteristics of Geoelectric Field at Xingji Seismic StationJ.North ChinaEarthquake Sciences，2023,41（3）：91-96.doi:10.3969/j.issn.10031375.2023.03.013.兴济地震台大地电场方位变化特征李瑞卿，郭学增，李 凤，张 蕾，信世民，尹宏伟（石家庄地震监测中心站,石家庄 050021）摘要：根据地电场成分分类取空间电磁环境较为平静且地电场观测数据日变形态较为稳定的时段对兴济台大地电场 ET进行提取，然后利用正交向两测道大地电场合成总的大地电场 ET 总。对总大地电场方位变化进行对比发现，兴济台午前午后大地电场方位角图像较平凉、大柏舍数据离散度高，无午前午后 3 个阶段的稳定形态变化特征；通过谐波分析及 F 检验发现，兴济台大地电场前 10 阶谐波拟合程度相对较低、高阶谐波成分显著，说明兴济台地电场成分复杂，可能是造成一天中大地电场方位角数据离散的原因。关键词：大地电场方位角；谐波分析；拟合程度；F 检验中图分类号：P315.722 文献标志码：A 文章编号：10031375（2023）03009106doi:10.3969/j.issn.10031375.2023.03.013 0 引言地电场是重要的地球物理场之一，其不仅能够反映空间电磁场的变化，还能反映与之有关的地壳和地表的电性结构及其变化。20 世纪中叶 Chapman提出，地电场由大地电场和自然电场组成1。地电场观测数据是观测方法、技术、场地环境、电磁环境等的综合体现2。随着对地电场的深入研究，有学者将自然电场 Esp与大地电场 ET进行提取及分析研究3-4，研究表明地电场波形日变化可能源于大地电场，其起源主要有空间电流系以及潮汐作用5-8；有学者对地电场稳定的峰-谷日变波形的研究探索出大地电场日变化的潮汐机理说5-6，认为地电场日变波形的物理过程是固体潮或空间电流导致裂隙水周期性渗流产生的过滤电场进而影响了地电场观测波形变化。裂隙水的渗透流动与地下裂隙结构密切相关，通过大地电场计算的优势方位角在一定程度上反映了地下裂隙的优势方位5-7。目前，许多震后分析研究表明优势方位角在震前存在不同形态的异常变化9-12。另外，基于谐波拟合计算的优势方位角对地电暴、雷电、高压直流输电等地电场干扰信号具有一定的抗干扰能力13-14。因此，大地电场优势方位角的研究对地震预测有一定的意义。地电场是矢量场，对于地电场的矢量变化研究发现，地电场在很多情况下是线性极化的，在地震发生前后地电场的极化方位会发生变化15-17。根据地电场极化的线性特征探究出极化斜率法18、垂直极化投影法19，并对地电场数据极化现象进行分析。大地电场是地电场主要构成之一，谭大诚等7研究表明，源于潮汐作用和空间电流系的大地电场其强度和方向受岩石裂隙结构影响，对大地电场矢量特性进行分析发现，平凉台地电场属于午前午后有稳定正弦波形的 TGF-B 型，其大地电场方位角午前午后存在 3 个阶段的稳定变化7。在对特殊布极方位的兴济台地电场研究过程中发现，同属于 TGF-B 型的兴济台大地电场方位角在午前午后数据离散，方位角不稳定。因此，需对这种不同的大地电场方位角变化形态研究分析，进一步了解兴济台地电场观测数据变化特征，为其他台站地电场观测数据分析及应用提供参考。1 大地电场 ET的提取及大地电场方位角计算地电场是较为稳定的自然电场 Esp和变化电场的叠加，其中变化的电场包括大地电场 ET和干扰收稿日期：2022-12-26基金项目：中国地震局监测、预报、科研三结合课题（3JH-202201056）第一作者简介：李瑞卿(1990)，男，工程师，主要从事地电观测与研究工作E-mail：第 41 卷第 3 期华 北 地 震 科 学Vol.41，No.32023 年 7 月NORTH CHINA EARTHQUAKE SCIENCESJul.，2023成分 Er（式 1）。当空间电磁环境较为平静、且地电场观测数据日变形态较为稳定时，视为干扰成分Er较少，那么地电场近似于自然电场 Esp和大地电场 ET的叠加。对地电场 E 进行 FFT 变换，前 10 阶谐波周期顺序是 2324 h、12 h、7.9 h、6 h、4.8 h、4 h、3.4 h、3 h、2.7 h、2.4 h3，将一整天 1 440 个分钟值数据求平均值（即日均值）可以消除大部分周期变化的大地电场 ET成分，进而得到稳定的自然电场 Esp（式 2）。使用地电场分钟值数据减去自然电场 Esp即可得到大地电场 ET分钟值数据（式 3）2-3。E=Esp+ET+Er(1)Esp=1439i=0Ei1440(2)ETi=EiEsp(3)提取各测向大地电场成分后可利用三角函数求解分钟值的总大地电场方位角。根据地电场强度分量值的方向约定2（以电流向北、向东为正，即地电场矢量“南北分量”的北指向南为正，“东西分量”的东指向西为正，“北东/北西分量”的偏北指向偏南为正），可以构建坐标系：以电场中心点为坐标系原点，指向正东方向为 x 轴正向，指向正北方向为 y 轴正向。根据大地电场值的正负判断大地电场方位，进而使用正交向大地电场及三角函数求解总大地电场方位角，此方位角 为 y 轴顺时针偏转的角度（图 1）。则一天中总大地电场 ET 方位角 计算公式为：i=180arcatan(|ETEWi|ETNSi|)(ETEWi0,ETNSi0,坐标系第一象限)360180arcatan(|ETEWi|ETNSi|)(ETEWi0,ETNSi0,坐标系第二象限)180+180arcatan(|ETEWi|ETNSi|)(ETEWi0,ETNSi0,坐标系第三象限)180180arcatan(|ETEWi|ETNSi|)(ETEWi0,ETNSi0,坐标系第四象限)(4)由于兴济台地电场布极方式较为特殊，布极方式为沿平行于沧东断裂带走向的 N30E 布极和垂直沧东断裂带走向的 N60W 布极（图 2）。对于布极方式比较特殊的兴济台，可以将N60W 及 N30E 两测向大地电场在东西、南北方向进行分解，求得东西、南北方向的大地电场，即：ET南北=ETN30Ecos30+ETN60Wcos60ET东西=ETN30Esin30+ETN60Wsin60(5)然后，再根据式 4 进行求解兴济台总大地电场方位角（式 5 与直接使用正交向 N30E 和 N60W 三角函数合成得到的总大地电场位置一致，此处为统一 y 轴指向正北的坐标系进行了转换）。2 大地电场 ET方位角对比谭大诚等7研究表明，平凉台大地电场方位角存在午前午后 3 个稳定的变化阶段，其相对稳定时刻的 3 个阶段为：（）6:0010:00 大地电场指向第三象限；（）其后大地电场方向发生反转，在 10:0016:00 大地电场指向第一象限；（）在 16:3019:30大地电场再次指向第三象限（图 3）。变化的大地电场潮汐机理认为：在固体潮汐作用或 Sq电流电磁感应下，裂隙水中电荷以周期性沿裂隙往返移动，进而造成大地电场 ET在不同阶段存在不同方位，此方位沿裂隙优势方向呈近 180反转。对兴济、平凉、大柏舍、昌黎 2021 年 10 月 28 日数据进行总大地电场方位角求解。兴济台总大地电场方位角散点图离散度高，而且没有午前午后3 个稳定阶段（图 4a）。平凉台总大地电场午前午后方位角数据稳定，其相对稳定时刻的 3 个阶段分别为：（）7:3011:00 大地电场指向第三象限，方位角为 247左右；（）11:3015:00 大地电场指向第 A1B1A2B2ONyET总x 图 1 大地电场方位角示意图Fig.1 Azimuth schematic diagram of the geoelectric field92华 北 地 震 科 学41 卷 一象限，方位角为 67左右；（）16:0018:30大地电场再次指向第三象限，方位角为 247左右（图 4b）。大柏舍台大地电场方位角与平凉台类似经历的 3 个阶段角度分别集中在 28、197、14左右（图 4c），昌黎台与兴济台类似一天中的大地电场方位角离散度高，无法辨别其稳定的形态（图 4d）。3 谐波分析拟合程度及 F 检验为了分析提取的大地电场成分，引入谐波分析。yt=y+n/2k=1Ckcos(2ktnk)=y+n/2k=1Akcos(2ktn)+Bksin(2ktn)(6)k=arcatan(BkAk),Ak0arcatan(BkAk),Ak02;Ak=0(7)AkBk其中，、由（8）式计算（k 表示谐波的阶数，同下）即Ak=2nnt=1ytcos(2ktn),Bk=2nnt=1ytsin(2ktn)(8)则各阶谐波周期及振幅为Tk=nk,Ck=A2k+B2k(9)yn/2k=1Ckcos(2ktnk)式 6 中：项为稳定的自然电场 Esp，即地电场日均值；项表示变化的电场：大地电场 ET及干扰电场 Er。116.0E116.5E117.0E117.5E118.0E39.0N39.0N39.0N39.0N图例兴济台直辖市断裂NB1B2A2A1O1、O2A1O1=300 mB1O1=300 mA1B1=424 mA2O2=150 mB2O2=150 mA2B2=212 mNa 构造地质 b 测区布极 图 2 兴济台构造地质和测区布极Fig.2 Geological Structure of Xingji Seismic Station and distribution of surveying area A1B1ONxyET(+)ET()第一象限 图 3 平凉台大地电场方位变化示意图Fig.3 Schematic diagram of azimuth variation of the geoelectric fieldon Pingliang Station3 期李瑞卿，等：兴济地震台大地电场方位变化特征93 理论上，使用谐波对地电场数据进行拟合，拟合使用的谐波阶数越多，拟合结果约接近真实的地电场观测变化。根据潮汐电场理论，地电场日变波形主体基本由前 10 阶谐波成分构成，使用谐波分析拟合程度及 F 检验分析兴济、平凉、大柏舍、昌黎地电场，可进一步分析大地电场成分。3.1谐波拟合程度Sk=12C2kni=1Sk=ni=112C2k选取电场空间电磁环境较为平静、且地电场观测数据日变形态较为稳定的 2021 年 10 月 16 日数据为例进行分析。各测向谐波拟合的解释率可以用各谐波的方差表示，即振幅平方的一半。用前 n 阶谐波来拟合观测数据则总的方差解释率为，则拟合程度 为解释率与观测k=ni=1Sk/序列（1 440 个分钟值）方差的比值：。兴济台 3 个测向前 10 阶谐波拟合程度均较低，大柏舍、平凉前 5 阶拟合程度就已经达到 90%以上。昌黎台绘制的一天中总大地电场方位角变化与兴济台相似，数据离散度高，其拟合程度同样较低（表 1）。3.2F 检验谐波函数可以检测数据样本的周期，第 k 阶对应周期 T=1 440/k，前 10 阶谐波拟合分别对应周期为 2324 h、12 h、7.9 h、6 h、4.8 h、4 h、3.4 h、3 h、2.7 h、和 2.4 h 的潮汐谐波。可构造 F 检验分析谐波拟合周期的显著程度。构造统计函数 F 检验：Fk=12C2k/2(s2C2k/2)/(n21)(10)表 1 各测向前 5、10、15、20、100 阶谐波拟合程度Table 1 The fitting degree of forward 5,10,15,20 and 100 order harmonics of each measurement 台站测向5/%10/%15/%20/%100/%方差兴济台N75E58.563.866.469.485.60.466N30E74.779.180.582.791.60.268N60W34.140.343.847.171.40.294大柏舍台NS91.795.896.296.598.81.886平凉台NS95.896.296.697.499.55.964昌黎NS56.963.567.769.589.71.983 0240480720960120014400100200300400b 平凉台优势方位角/()0100200300400优势方位角/()0100200300400优势方位角/()0100200300400优势方位角/()t/min024048072096012001440t/min024048072096012001440t/min024048072096012001440t/mina 兴济台c 大柏舍台d 昌黎台 图 4 2021 年 10 月 28 日大地电场方位角散点图Fig.4 Scatter diagram of geoelectric field azimuth on October 28,202194华 北 地 震 科 学41 卷 则 F 检 验 显 著 水 平 阈 值 为：fa=finv(0.95，2，1337)=3.0025当 Fk超过 F 阈值时说明第 k 谐波是显著的，即拟合 k 阶谐波周期成分有显著贡献度。对总大地电场方位角变化形式不同的 4 个台站数据进行 F 检验：一天中总大地电场方位角午前、午后数据稳定集中大柏舍、平凉台水平显著的k 阶谐波仅限制于前 10 阶中；而兴济台各测向地电场显著周期不仅存在于前 20 阶，更高阶谐波依然显著（如 N75E 测向在 79 阶依然显著，其对应周期为1 440/79 约为 18 分钟），说明未知的干扰成分叠加到大地电场之中。方位角数据离散的昌黎台地电场同样存在高阶的显著周期成分，表 2 中仅展示前30 阶谐波显著水平。表 2 各测向各阶谐波显著水平Table 2 Significant levels of each harmonics wave in each direction 台站测向12345678910兴济台N75E29.17434.29120.5413.495.4314.1410.553.452.767.24N30E10.01767.90146.5031.593.8015.014.465.335.381.45N60W46.20149.2976.773.415.2712.9314.323.190.5914.03大柏舍台NS67.50251.67646.3569.018.2915.456.243.583.810.27平凉台NS38.91253.38682.78116.994.462.991.306.270.940.36昌黎台NS69.5749.88235.43130.3811.540.9917.118.8314.876.15台站测向11121314151617181920兴济台N75E1.360.7811.323.911.439.105.713.132.930.90N30E0.770.773.981.463.769.071.481.731.271.91N60W1.571.1817.065.760.266.759.693.604.020.29大柏舍台NS1.021.170.180.430.051.270.190.390.490.01平凉台NS1.370.550.260.400.562.780.431.370.600.50昌黎台NS17.905.633.733.120.592.961.521.841.984.47台站测向21222324252627282930兴济台N75E1.592.196.230.715.120.240.273.481.110.04N30E0.741.854.180.812.960.181.250.923.040.59N60W2.892.136.322.275.660.600.277.651.280.16大柏舍台NS1.800.070.680.230.050.050.540.200.350.10平凉台NS1.121.232.090.150.280.560.870.810.010.12昌黎台NS2.944.660.811.170.322.2111.770.982.328.95 3.3讨论基于大地电场潮汐理论，这种午前、午后的数据稳定变化可能更多地受前 10 阶潮汐电场成分影响。平凉、大柏舍地电场前 10 阶就有较高的拟合程度，且 F 检验结果显著周期限制在前 10 阶谐波中，说明了平凉、大柏舍地电场数据波形变化基本由前 10 阶潮汐电场成分组成。而兴济、昌黎地电场数据则需要更多高阶谐波参与拟合来增加其拟合程度，说明这两个台站地电场成分中代表潮汐作用的前 10 阶周期成分外，高阶周期成分存在较多。通过 F 检验同样证实谐波拟合过程中高阶谐波成分显著。对大地电场提取过程进一步分析后认为：在式（1）（3）大地电场计算中提取过程中，对自然电场 Esp稳定过于理想化，实际中自然电场Esp同样在一天中不断变化的；提取过程中通过求平均值的方式来消除周期成分并不理想，周期信号可能存在残留成分（如 23 h 的日潮周期在 1 440个分钟值相加后存在残留，并不会完全消除）；噪声 Er 同样理想化，场地环境中的各种噪声（如轻轨干扰、施工干扰等）同样夹杂在求解得到的大地电场中。地电场本身成分复杂，在实际不同的场地环境下提取的大地电场效果不同。计算过程中，高阶周期成分、自然电场变化成份、周期残留成分、白噪成分都可能造成最终提取的结果并非完全是大地电场ET成分，这可能是求解得到的数据离散度高的原因。3 期李瑞卿，等：兴济地震台大地电场方位变化特征95 4 结论兴济台总大地电场方位角在一天中的数据变化离散，平凉、大柏舍总大地电场方位角存在午前午后数据稳定的 3 个阶段。通过谐波拟合程度及F 检验分析得到：总大地电场方位角数据离散的兴济、昌黎台地电场相对与平凉、大柏舍在谐波拟合过程中需高阶成分显著，推断在不同场地环境下，高阶周期成分、自然电场变化成份、周期残留成分、白噪成分都可能造成最终提取的大地电场ET成分复杂，进而造成兴济台总大地电场方位角数据离散。致谢感谢中国地震局兰州岩土地震研究所安张辉、谭大诚及中国地震台网中心叶青等老师的帮助和支持，感谢本文编辑和评审专家的辛苦劳动及有益帮助！参考文献：孙正江,王华俊.地电概论 M.北京:地震出版社,1984.1DB/T 18.2-2006,地震台站建设规范 地电观测台站 第 2 部分:地电场台站 S.北京:地震出版社,2006.2谭大诚,赵家骝,席继楼,等.青藏高原中强地震前的地电场变异及构成解析 J.地球物理学报,2012,55(3):875-885.3谭大诚,赵家骝,刘小凤,等.自然电场的区域性变化特征 J.地球物理学报,2014,57(5):1588-1598.4谭大诚,赵家骝,席继楼,等.潮汐地电场特征及机理研究 J.地球物理学报,2010,53(3):544-555.5谭大诚,王兰炜,赵家骝,等.潮汐地电场谐波和各向波形的影响要素 J.地球物理学报,2011,54(7):1842-1853.6谭大诚,辛建村,王建军,等.大地电场岩体裂隙模型的应用基础与震例解析 J.地球物理学报,2019,62(2):558-571.7黄清华,刘涛.新岛台地电场的潮汐响应与地震 J.地球物理学报,2006,49(6):1745-1754.8侯泽宇,安张辉,范莹莹,等.云南漾濞县 6.4 级地震前地电场异常回溯性分析 J.地震工程学报,2021,43(4):807-817.9辛建村,谭大诚.地电场多测向日变波形相位关联特征 J.地震学报,2017,39(4):604-614.10辛建村,张晨蕾.漾濞县 MS6.4、玛多县 MS7.4 地震前地电场变化分析 J.地震工程学报,2021,43(4):818-825.11艾萨伊斯马伊力,陈界宏,次旦,等.和田台地电场岩体裂隙优势方位角变化特征分析 J.内陆地震,2021,35(3):245-252.12张国苓,任印国,乔子云,等.华北、东北地区地电场优势方位角分析 J.地震地磁观测与研究,2021,42(6):89-97.13张波,谭大诚,罗娜,等.高压直流输电对大柏舍台地电场观测干扰的分析 J.地震,2020,40(1):140-150.14安张辉,杜学彬,谭大诚,等.四川芦山 Ms7.0 和汶川 Ms8.0 地震前地电场变化研究 J.地球物理学报,2013,56(11):3868-3876.15毛桐恩,席继楼,王燕琼,等.地震过程中的大地电场变化特征 J.地球物理学报,1999,42(4):520-528.16李旭升,张远富.平凉地震台地电场不极化现象探讨 J.高原地震,2018,30(1):40-45,51.17沈红会,冯志生,燕明芝,等.地电场震前变化的探讨 J.西北地震学报,2006,28(1):74-77.18阮爱国,赵和云.提取地震地电场异常的垂直极化方向投影法 J.地震学报,2000,22(2):171-175.19AzimuthVariationCharacteristicsofGeoelectricFieldatXingjiSeismicStationLI Ruiqing,GUO Xuezeng,LI Feng,ZHANG Lei,XIN Shimin,YIN Hongwei（Shijiazhuang Seismic center station,Shijiazhuang 050021,China）Abstract:Based on the classification of the geoelectric field,the ET of geoelectric field of Xingji seismic station isextracted at the time when the space electromagnetic environment is calm and the daily deformation of the geoelectricfield is stable.Then,the total ET of the geoelectric field is synthesized by using the orthogonal two-channelgeoelectrical field.The azimuth variation of the total geoelectrical field is analyzed:the azimuth angle is calculated andthe scatter plot is drawn.The azimuth of the geoelectric field in Xingji Seismic Station is more discrete than that inPingliang and Dabaisha.Tha data does not have the characteristics of three-stage steady-state variation.Harmonicfitting was used to analyze the harmonics of the electric field data of Xingji,Pingliang,Dabaishe and Changli stations.The differences of harmonic components of these stations were analyzed by harmonic fitting degree and f-test,furthermore,the reason of the instability of the azimuth data before and after noon in geoelectric field of Xingji sataionis analyzed.The low fitting degree of the first 10 order harmonics and the obvious high order harmonics in XingjiSeismic Station geoelectric field indicate that the complex composition of the field in Xingji Seismic Station may bethe cause of the dispersion of the azimuth data of the geoelectric field in a day.Finally,it provides a reference forunderstanding the composition and extraction of the geoelectric field.Keywords:geoelectric field azimuth；harmonic analysis；fitting degree；F test96华 北 地 震 科 学41 卷