2. 中国地震局兰州地震研究所, 兰州 730000;
3. 中国地震台网中心, 北京 100045
2. Lanzhou Institute of Seismology, China Earthquake Administration, Lanzhou 730000, China;
3. China Earthquake Networks Center, Beijing, 100045, China
地球电场(简称“地电场”主要由固体地球外部的电流体系在地球内部感应产生的分布于地表的快变化电场和地下介质中的物理、化学作用引起的慢变化电场两部分组成.地电场观测在资源勘探、地质调查等方面应用广泛,利用固定台站观测地电场的主要目的之一是进行地震等灾害事件的监测预报和地球电磁环境研究.20世纪初以来,国内外学者开展了固定台站的地电场观测和研究,对地电场时空变化的分布规律取得了重要认识[1-8],吉什和鲁尼在总结长期观测资料的基础上绘制了第一幅全球涡旋电流场分布图[1] (1936年世界时18h).地电场观测在地震与火山等灾害事件监测预报实验中也得到广泛应用[9-18].我国从“九五”后期开始了大空间范围、定点、连续、规范化的地电场观测[19-20],截止2011年我国已建成由100多个台站组成的地电场观测台网,其中部分台站持续观测时间已超过10年以上.国内进行的地电场观测是目前国内外地球电磁学观测的一大特色,目前关于地电场在地震监测预报实验中的研究工作较多,但关于地电场时空变化规律的研究甚至于国外都很薄弱.文献[21-23]应用我国部分地电场台站的观测数据研究了地电场日变化(以下简称“地电日变化”的时频分布特征,但数据应用不具中国大陆大区域的特点,研究成果有限、且认识有所不同,对于认识地电日变化的时空分布规律来讲是不够的.鉴于我国地电场台网在某些区域密度高、数据积累丰富,本文在文献[21-23]的研究基础上,应用中国大陆东部、西部、中低纬度、中纬度各两条经度链、纬度链的地电场台站观测数据,研究地电日变化的时空分布规律,期望能获得对中国大陆地电日变化的广域性、普遍性特点的进一步认识.
2 数据选取我国地电场台站观测装置是水平布设NS、EW两个正交观测方向再加一个斜交方向,长、短极距共6个测道[20],每测道采样率为1次/min,仪器分辨率为10 μV,24 h连续观测.为了研究地电场在中国大陆大空间区域的时、空变化规律,选取我国地电场台网中两条经度链、两条纬度链台站的观测数据.选取原则为:①台站分布范围尽量大,空间距离尽可能均匀;②尽可能选取经度链、纬度链上及近距离台站; ③台站观测数据质量合格.如图 1所示,沿118°E经线选14个台站组成经度链a,由南向北台站依次为:福建泉州、邵武,江苏海安、高邮,安徽嘉山,山东郯城、无棣,河北大柏舍、昌黎,北京通州,内蒙古宝昌,吉林榆树,黑龙江绥化、满洲里;沿99°E经线选12个台站组成经度链b,由南向北依次为:云南腾冲、洱源,四川泸沽湖、甘孜,青海大武、都兰,甘肃寺滩、古浪、山丹、高台、嘉峪关,宁夏石嘴山;沿40°N纬线选7个台站组成纬度链c,由西向东依次为:新疆乌什,甘肃瓜州、古浪、嘉峪关,山西大同,北京延庆、通州;沿34°N选12个台站组成链d,由西向东依次为:新疆和田,青海都兰、大武,甘肃陇南、静宁,陕西周至、乾陵,河南洛阳、周口,山东郯城,江苏高邮、海安.由于台站观测环境影响,某些近距离台站未人选.
近年来,随着我国地电场台网的建设和台网规范化运行,执行了严格的观测质量监控程序,整体上数据质量有所提高,所以选取资料的时段主要为2010、2011年.个别台站受环境干扰,使用2008年以来其他年份的观测数据.
3 地电场日变化波形 3.1 波形变化根据亥姆霍兹波茨坦中心(GFZ)公布的磁情指数,尽可能选取每月磁静日、且日期接近天数的观测数据.数据处理采用时序迭加方法,突出有规律性的日变化周期、抑制不规则噪声.时序迭加公式如下:
(1) |
式(1)中Xn(t)表示所选时间范围内第n天的地电场数据,每天的数据长度为1440个(分钟值采样),k指第k天.
以经度链a上台站观测数据为主要研究对象. 图 2是链a中8个地电场台站每个台2010年2月3个磁静日(或一般磁静日)的观测数据时序迭加曲线.按纬度由低到高的顺序从上到下排列.由图明显可见:地电日变化主要表现为两峰一两谷变化波形,其中低纬度的泉州、海安、高邮台以及中低纬度的郯城、无棣和中高纬度的绥化台日变化波形表现为两峰一两谷波形;满洲里台表面上是中午前后1峰一1谷变化波形.
链a台站分布在25. 02°N-49. 57°N的纬度范围,由图 2似乎可见,总体上台站纬度高,日变化峰一谷幅度
由图 2可见,链a各台站在午时(即12h左右)日变化出现波峰和波谷.部分台站之间存在少许相位差,是由于链a台站并不全在118°E经线上,位于此经线附近,经度偏东的台站相位早于经度偏西台站.例如,绥化台经度为126. 9°E、满洲里台117.43°E,绥化台波峰、波谷极值时间早于满洲里台,而满洲里、郯城、泉州台(自东向西顺序)经度接近,基本无相位差.
选取链d台站观测数据,并按劳埃德(Loyd)天文季节的D季节(11、12、1、2月)、季节(5、6、7、8月)和E季节(3、4、9、10月)进一步研究日变化随着地方时的变化(即“经度效应”).在2011、2008、2010年中不同的天文季节中各选一个磁静日,比较链d台站在同一天的地电日变化波形(图 3).图中各台依从东到西的顺序上下排列,这些台站分布在79. 9°E-120.5°E之间,最大台间距超过3700 km.通过图中的辅助线可清晰地看到日变化相位有经度效应,随着经度偏西,日变化相位依次滞后.例如,2011年1月5日和2008年3月7日的日变化相位的经度效应最明显,2010年8月22日较明显,但不如前两个图(因数据质量影响).从1月5日地电场日变化可明显看出:工苏海安、新疆和田台第二个波谷到达的时差约为2.7h,与两台站经度差(约40. 6°)对应的地方时差(约2. 7h)相符(见图中辅助线);高邮一和田、郯城一和田波形相位约2. 6h,吻合台站当地的时差;其次,海安一周至、周至一都兰、都兰一和田台时差依次分别为0. 82、. 9、1h,日变化相位差依次分别为0. 82、. 08、1h.显然,地电场日变化两次起伏的相位存在明显的经度效应,相位差与台站之间当地的时差很吻合.
同时还注意到,地电场快变化成分是相当同步的,无时差(见2011年1月5日图 0~240 min和1320 min前后、2010年8月22日图 0~120min).实际上,在数据处理中注意到此类快变化现象是普遍存在的,特别是在地电暴期间.说明地电日变化相位存在经度效应,而快变化同步,似乎反应了二者的“源”因素有区别.
分析了经度链b、纬度链c台站地电日变化波形、经度效应,其结果与经度链a、纬度链d的结果类似.另需说明,链a_d个别台站地电日变化波形、经度效应不很符合前面的结果,这种情况可能与台站观测环境、技术系统运行等因素有关; 个别台站日变化波形还表现为1峰一1谷的形态.其次,极个别台站观测数据不可用,文中未应用.
4 地电场日变化周期成分 4.1 数据处理在地电场频率域变化的分析中,用快速傅里叶变换方法(FFT)处理了4条经、纬链上共14个台站2010-2011两年内116个月的地电场分钟值数据.由于日变化是频率极低的成分,为了客观获得日变化的周期成分,每个台站选取3~5个月的地电场北分量Ex和东分量Ey观测数据,且在选取数据时段的每个月各分量分钟值数据组成1个长度约为43200的时间序列,使时间序列长度内能足够包涵多个日变化周期成分.
4.2 周期特性图 4a是嘉山台2010年3月地电场Ex和Ey观测数据的FFT振幅谱.从图中可以清楚地看到, 优势周期的幅值从大到小排列依次为:12.4、8、24h、12、6、4.8h, 其中12.4h为最优势周期.在对其他台站数据的处理过程中, 较明显的优势周期中还出现了25.6、4h周期, 但其幅度明显小于前面几个周期成分的幅度.处理其他月份的观测数据, 其结果类似图 4a.图 4b为2010年全年嘉山台地电场Ex和Ey观测数据FFT振幅谱时频图.由图可见:① 6h、8、12.4、24h周期成分的谱幅度相对其他周期成分占优势, 其中12.4h周期成分的谱幅在全年显著, 24h周期成分在夏季比其他季节显著; ②夏季幅值大, 冬季幅值小, 特别是Ey分量更明显.
表 1是2011年1、2、3、7、10共5个月, 链a的泉州、海安、高邮、嘉山4个台站两个测道的地电场数据FFT振幅谱统计的优势周期.除表 1中的4个台外, 还对链a的郯城、无棣、绥化、满洲里4个台观测数据做了谱分析(篇幅所限未列出).8个台每台每月各分量谱分析结果记为一组, 共计80组, 在振幅谱中可明显识别的优势周期共有333个.表 1中每组数据依谱幅从大到小排列.由表可见, 12.4h、12h以及24h周期为优势周期, 而12.4h周期为优势周期中出现频次最高的周期成分.其次8、6h等周期成分也多次出现, 还间或出现约25.6h的周期成分.同时注意到, 除泉州台1、2、3月外, 其余台站Ex、Ey两个分量优势周期相同或接近, 其次不同月份优势周期也相同.
表 2是对链a、链b台站地电场数据FFT振幅谱的统计(链b取腾冲、泸沽湖、古浪、石嘴山、山丹、高台6个台2011年2、8、10月的数据).可以看出链a台站每月每台站出现8 h和12 h周期成分的频次最高,按频次统计组发生率分别为87. 5%和80%,6、24、12. 4h三种周期的组发生率相近,也都在60%以上,其次还有4.8、25. 6h和4h三种周期成分,但明显不及前5个周期成分出现的频次高.其他周期的组发生率的总和仅为13. 8%.因为地电场观测对环境干扰、电极稳定性、线路绝缘以及地下介质电性变化等因素敏感,易发生数据异常,少数不常见的周期成分可能是由于以上因素的影响.需指出,由于数据质量和谱估计误差,统计中提到的某一周期是指此周期±0.1 h范围(如12.4 h指12.4 h±0.1 h).与链a相比,链b台站每组FFT谱中可分辨的优势周期的个数平均为3. 11个,低于链a台站每组平均4. 16个,可能因为链b位于我国西部,台站环境艰苦、观测数据质量整体上不如链a台站.如云南腾冲、四川泸沽湖台存在观测场地环境干扰,可分辨的周期较少或无优势周期.链b台站主要周期按频次依次为12、8、24h,组发生率在60%以上,其中12h组发生率最高,为80. 6%;其余组发生率较高的周期成分为6、4 8h与链a不同的是,在链a台站出现12. 4h周期成分的频次高,组发生率高达62. 5%,而链b台站仅为2. 78%,25. 7h周期的频次也仅为2. 78%.
按最大FFT谱幅出现的频次统计(表 3),在链a的80组数据(即80个月数据)中,12. 4h周期成分的频次最高,为42. 5%,其余依次为12、8、24h.说明在链a的8个台站中,地电日变化主要是由12. 4、12、8h和24h周期成分组成,其中12. 4h周期幅值最显著,虽然4. 8、25. 7、4h周期波也较常见,但非普遍存在、振幅不大;连b台站36组数据中最大谱幅的周期成分出现频次最高为12h,其次是8h和24h周期成分.
总之,在两个经度链上地电日变化可分辨的优势周期按周期成分的振幅大小主要有12. 4、12h、8、24、6、4. 8、4、25. 7h,但两条经度链上各自的优势周期成分又有区别.链a台站12. 4h周期最显著,链b台站12h周期成分最优显著.
对郯城2010年1月6-8日三个磁静日的地电场观测数据进行时序迭加处理,提取其中振幅最大的前5阶谐波进行反变换拟合.图 5中灰色曲线为原始曲线,黑色为拟合曲线,拟合曲线很好地近似了原始曲线.前5阶谐波振幅依次为1. 67、1. 01、0.88、0.33和0. 26 mV,其中12h和8h周期成分的振幅最大.
文献[24]给出了用FFT处理的美国图森台1个月地磁场Sq日变化的主要周期,按幅值大小其主要周期成分依次为24 h、12 h和8 h周期等,说明地磁场Sq日变化主要周期成分为每日一次起伏.本文用每台、每测道1个月的数据组成长时间序列,其周期成分估算应准确,按幅值大小分别为12.4、12h、8、24 h,与文献[21-22]用最大熵方法估算的昌黎、宝坻、嘉峪关台地电日变化主要周期成分为12、24 h和8 h左右周期的结果吻合,说明地电日变化每日两次起伏为主,符合图 2-3中地电日变化波形.需说明:及个别台站日变化周期不是12. 4h/12h,例如满洲里台主要为8 h周期.
5 地电日变化幅度
由于地电日变化主要为每日两次起伏,故对日变化幅度的分析针对12 h/12.4 h半日波.选取链a台站2010年每个天文季节中1个月磁静日3天数据统计日变化峰一谷幅度
同时,由图 6(a-b)得到,J季节的红色曲线高于E、D季节的绿色和黑色曲线,而D季节黑色曲线最低.分别统计图 6(a-b)展示的三个季节各台站的
选取链b的12个台站观测数据,按链a方法统计
应用纬度链d自西到东11个台站的地电场Ex、Ey,分量观测数据,按链a方法分析了各台站的
吉什和鲁尼据全球地电场数据绘制了1936年世界时18时的全球涡旋地电流分布图的概貌[1],在北半球的中低纬度区、高纬度区各分布4个涡旋场,中低纬度涡旋中心在30°N附近,在高纬度区接近70°N附近,各涡旋分布相对于太阳的位置固定.其中,中低纬度涡旋中心大致在25°N、120°E和35°N、70°E的2个涡旋电流场覆盖或邻近中国大陆(在世界时18时),东部链a、西部链b两个子午链台站分布在25°N -49. 57°N的中低一中纬度范围,两个涡旋电流场覆盖链a、链b台站.由于各涡旋分布相对太阳的位置固定不变,在一天当中世界时的另外时刻,链a、链b台站还要经过另外两个中低纬度区的涡旋电流场•对于链a,泉州(25. 02°N、118. 51°E) -满洲里(49.57°N、117.43°E)共14个台站日变化幅度
地电场Sq和低磁情日变化主要为12. 4h (链a)、12h (链b)、8h和24h周期成分以及6、4.8h、4h等成分,其中12. 4h周期是主太阴半日分潮(M2波),12h周期是主太阳半日分潮(S2波),24h是太阳赤纬全日分潮(P1波),8、6、4. 8、4h也是潮汐调和分量.在这些周期成分中,24、12、8、6、4. 8 h成分与变化地磁场Sq谐波成分一致,文献[24]关于变化磁场中的上述谐波分量已做了详细研究,认为这些成分与中高层大气及电离层中的潮汐发电机理论息息相关.月球、太阳依次是地球潮汐现象的主要贡献者,链a台站地电日变化最主要是M2波周期,链b台站最主要是S2波周期,这些周期成分是主要的潮汐调和分量,因此认为地电日变化的主要周期与月日引潮力直接有关.半日潮每天两次高潮和两次低潮,全日潮每天一次循环,在引潮力作用下固体地球、大气层、电离层、海洋、湖泊等产生的潮汐现象引起了地面电流涡旋场强度发生周而复始的变化,引起了上述地电日变化的主要周期成分,即以12.4h、12、24h周期成分的幅度为最显著的日变化.
文献[23]研究了潮汐对岩石中裂隙走向的优势趋向作用和裂隙水渗流对地电场潮汐响应的影响.地壳介质中普遍存在裂隙诱发的各向异性[25],实际上该文讨论了介质中微裂隙走向优势趋向(或定向排列)构成的各向异性介质中传导电流引起的地电场变化的情况,导电机制主要是介质中的动电效应,导电溶液运移和微裂隙、孔隙是重要因素.地电场观测地表两个水平正交分量,本文应用上述微裂隙构成的方位各向异性介质进一步讨论月日引潮力对地电日变化的影响.在这种介质中欧姆定律是张量表示式,在微裂隙优势趋向方向及其垂直方向(即各向异性介质的两个水平电性主轴方向)的电场分量强度
另一方面,空间电流系也是影响地电日变化的重要因素.地电日变化的午前午后表现类似于变化磁场Sq,据文献[24, 27]变化磁场Sq午前午后的不对称现象起因于低纬电流系和高纬电流系.其次,比较分析地电日变化波形与钻孔应变、静水位、重力等观测的半日潮波形可知,地电日变化与应变等半日波存在明显的相位差和视周期差异.这种差异表明电离层、大气层等的电动力学过程也是影响地电日变化的重要原因.因此,本文认为月日引潮力引起的地面涡旋电流场电流强度的变化和太阳风等引起的空间电磁活动共同产生了地电日变化的半日波成分.
链a台站一般位于我国东部海拔低于500 m的区域,且邻海/湖等水系居多,日变化半日波周期为12. 4h,而链b台站一般位于海拔约1500 m及以上的高原区,半日波周期为12h.其周期差异的原因可能是链a台站地下介质内部水网络发育,月球引潮力的作用更明显,对于链b台站空间电磁活动的作用更明显.
7 结论
(1) 沿两个经度链、纬度链台站的地电日变化绝大多数表现为两峰一两谷的变化波形,波峰、波谷出现在当地时午前午后时段,一般午前日变化幅度
(2) 地电日变化周期成分按最大幅度出现的频次排序依次为12. 4、12、8、24、6、4. 8、4、25. 7h周期等,这些周期是潮汐调和分量.其中在118°E链台站以主太阴半日波(M2波)12. 4 h为最主要周期,99°E链以主太阳半日分波(S2波)12h为最主要周期.
(3) 纬度链地电日变化半日波相位差与台站之间的地方时的时差一致,表现为经度偏东台站早于偏西台站的经度效应.
(4) 地电日变化半日波周期成分的幅度有纬度效应,总体上,118°E链从福建泉州一内蒙古满洲里台日变化幅度逐渐减小,99°E链从约35°N开始大致随纬度分别增高和降低而日变化幅度减小;Loyd天文季节的夏季日变化幅度最大,春秋季节其次,冬季最小;绝大多数台站地电场北分量艮日变化幅度大于东分量Ex;地电日变化幅度与台址浅层电阻率存在正比关系.
(5) 月日引潮力引起的地面涡旋电流场的电流强度变化和太阳风引起的大气层、电离层电磁活动共同产生了地电日变化的半日波周期成分.
致谢十分感谢两位审稿专家和编辑部对本文付出的辛勤工作,审稿专家对文中逐句修改,并提出多条具体的修改意见和建议,体现了科学、严谨的学风.
[1] | 克拉耶夫.地电原理.张可迁, 陈培光, 张志诚等译.北京:地质出版社, 1954: 100-236. PAEB. Geoelectrics Principle (in Chinese). Zhang K Q, Chen P G, Zhang Z C, et al. Trans. Beijing: Geological Publishing House, 1954: 100-236. |
[2] | Lanzerotti L J, Sayres D S, Medford L V, et al. Response of large-scale geoelectric fields to identified interplanetary disturbances and the equatorial ring current. Adv. Space Res. , 2000, 26(1): 21-26. DOI:10.1016/S0273-1177(99)01021-2 |
[3] | Pirjola R. Effects of space weather on high-latitude ground systems. Adv. Space Res. , 2005, 36(12): 2231-2240. DOI:10.1016/j.asr.2003.04.074 |
[4] | Burns G B, Hesse M H, Parcell S K, et al. The geoelectric field at Davis station, Antarctica. Journal of Atmospheric and Terrestrial Physics , 1995, 57(14): 1783-1797. DOI:10.1016/0021-9169(95)00098-M |
[5] | Corney R C, Burns G B, Michael K, et al. The influence of polar-cap convection on the geoelectric field at Vostok, Antarctica. Journal of Atmospheric and Solar-Terrestrial Physics , 2003, 65(3): 345-354. DOI:10.1016/S1364-6826(02)00225-0 |
[6] | 傅承义, 陈运泰, 祁贵仲. 地球物理学基础. 北京: 科学出版社, 1985 : 203 -254. Fu C Y, Chen Y T, Qi G Z. Geophysical Fundamentals (in Chinese). Beijing: Science Press, 1985 : 203 -254. |
[7] | 孙正江, 王华俊. 地电概论. 北京: 地震出版社, 1990 : 3 -95. Sun Z J, Wang H J. Introduction of Geoelectricity (in Chinese). Beijing: Seismological Press, 1990 : 3 -95. |
[8] | 李金铭. 地电场与电法勘探. 北京: 地质出版社, 2005 : 36 -59. Li J M. Geoelectric Field and Electrical Exploration (in Chinese). Beijing: Geological Publishing House, 2005 : 36 -59. |
[9] | Uyeda S. Introduction to the VAN Method of Earthquake Prediction. Singapore: World Scientific Publishing Co., 1995 : 3 -28. |
[10] | Varotsos P, Alexopoulos K. Physical properties of the variation of the electric field of the earth preceding earthquake. Tectonophysics , 1984, 110(1-2): 73-98. DOI:10.1016/0040-1951(84)90059-3 |
[11] | Zlotnicki J, Kossobokov V, Le Mou l J L. Frequency spectral properties of an ULF electromagnetic signal around the 21 July 1995, M=5.7, Yong Deng (China) earthquake. Tectonophysics , 2001, 334(3-4): 259-270. DOI:10.1016/S0040-1951(00)00222-5 |
[12] | 李宁, 杜学彬, 谭大诚, 等. 松山观测台地震短临电磁现象. 地震 , 2007, 27(Supple): 103–111. Li N, Du X B, Tan D C, et al. Imminent electro-magnetic phenomenon related to earthquakes recorded at the Songshan Station. Earthquake (in Chinese) , 2007, 27(Supple): 103-111. |
[13] | 范莹莹, 杜学彬, ZlotnickiJ, 等. 汶川Ms8.0大震前的电磁现象. 地球物理学报 , 2010, 53(6): 2887–2898. Fan Y Y, Du X B, Zlotnicki J, et al. The electromagnetic phenomena before the Ms8.0 Wenchuan Earthquake. Chinese J. Geophys. (in Chinese) , 2010, 53(6): 2887-2898. |
[14] | 刘君, 杜学彬, ZlotnickiJ, 等. 几次大震前的地面和空间电磁场变化. 地球物理学报 , 2010, 54(11): 2885–2897. Liu J, Du X B, Zlotnicki J, et al. The changes of the ground and ionosphere electric/magnetic fields before several great earthquakes. Chinese J. Geophys. (in Chinese) , 2010, 54(11): 2885-2897. |
[15] | 汤吉, 詹艳, 王立凤, 等. 5月12日汶川8.0级地震强余震观测的电磁同震效应. 地震地质 , 2008, 30(3): 739–745. Tang J, Zhan Y, Wang L F, et al. Coseismic signal associated with aftershock of the Ms8.0 Wenchuan Earthquake. Seismology and Geololgy (in Chinese) , 2008, 30(3): 739-745. |
[16] | 黄清华, 刘涛. 新岛台地电场的潮汐响应与地震. 地球物理学报 , 2006, 49(6): 1745–1754. Huang Q H, Liu T. Earthquakes and tide response of geoelectric potential field at the Niijima station. Chinese J. Geophys. (in Chinese) , 2006, 49(6): 1745-1754. |
[17] | 郝锦绮, 钱书清, 高金田, 等. 岩石破裂过程中的超低频电磁异常. 地震学报 , 2003, 25(1): 102–111. Hao J Q, Qian S Q, Gao J T, et al. ULF electric and magnetic anomalies accompanying the cracking of rock sample. Acta Seismologica Sinica (in Chinese) , 2003, 25(1): 102-111. |
[18] | 钱书清, 吕智, 任克新. 地震电磁辐射前兆不同步现象物理机制的实验研究. 地震学报 , 1998, 20(5): 535–540. Qian S Q, Lü Z, Ren K X. Experimental study on the mechanism of non-synchronism of seismo-electromagnetic radiation precursors. Acta Seismologica Sinica (in Chinese) , 1998, 20(5): 535-540. |
[19] | 钱家栋, 顾左文, 赵家骝等.地震台站观测环境技术要求.第2部分:电磁观测(GB/T195312-2004).北京:地震出版社, 2004. Qian J D, Guo Z W, Zhao J L, et al. Technical Requirement for Observational Environment of Seismic. Stations-Part 2: Electromagnetic Observation (GB/T 19531.2-2004) (in Chinese). Beijing: Seismological Press, 2004. |
[20] | 杜学彬, 席继楼, 谭大诚等.地震台站建设规范地电台站第2部分:地电场台站(DB/T 18.2-2006).北京:地震出版社, 2004. Du X B, Xi J L, Tan D C, et al. Specification for the Construction of Seismic Station Geoelectrical Station Part 2: Geoelectrical Field Observatory (DB/T 18.2-2006) (in Chinese). Beijing: Seismological Press, 2006. |
[21] | 杜学彬, 叶青, 赵杰, 等. 地电场日变化研究. 地震 , 2007, 27(Supple): 121–130. Du X B, Ye Q, Zhao J, et al. Study on geoelectric field daily variation. Earthquake (in Chinese) , 2007, 27(Supple): 121-130. |
[22] | 叶青, 杜学彬, 周克昌, 等. 大地电场变化的频谱特征. 地震学报 , 2007, 29(4): 382–390. Ye Q, Du X B, Zhou K C, et al. Spectrum characteristics of geoelectric field variation. Acta Seismologica Sinica (in Chinese) , 2007, 29(4): 382-390. |
[23] | 谭大诚, 王兰炜, 赵家骝, 等. 潮汐地电场谐波和各向波形的影响要素. 地球物理学报 , 2011, 54(7): 1842–1853. Tan D C, Wang L W, Zhao J L, et al. Influence factors of harmonic waves and directional waveforms of tidal geoelectrical field. Chinese J. Geophys. (in Chinese) , 2011, 54(7): 1842-1853. |
[24] | 徐文耀. 地磁学. 北京: 地震出版社, 2003 : 221 -349. Xu W Y. Geomagnetism (in Chinese). Beijing: Seismological Press, 2003 : 221 -349. |
[25] | Negi J G, Saraf P D. Anisotropy in Geoelectromagnetism. Elsevier Science Publishing Company , 1989(1-88): 175-232. |
[26] | 杜学彬, 李宁, 叶青, 等. 强地震附近视电阻率各向异性变化的原因. 地球物理学报 , 2007, 50(6): 1802–1810. Du X B, Li N, Ye Q, et al. A possible reason for the anisotropic changes in apparent resistivity near the focal region of strong earthquake. Chinese J. Geophys. (in Chinese) , 2007, 50(6): 1802-1810. |
[27] | 赵旭东, 杜爱民, 徐文耀, 等. Sq电流系午前午后不对称性现象的来源. 地球物理学报 , 2008, 51(3): 643–649. Zhao X D, Du A M, Xu W Y, et al. The origin of the prenoon-postnoon asymmetry for Sq current system. Chinese J. Geophys. (in Chinese) , 2008, 51(3): 643-649. |