Interference characteristics of new energy power generation on electromagnetic observation environment：Taking the electromagnetic observation of Xuzhuangzi station of Tianjin as an example
-
摘要: 利用天津市徐庄子地震台的电磁观测数据，针对台站周边观测区域内近年来出现的风电和光伏发电设施，通过野外试验测量，从多个角度深入剖析了不同观测环境中电磁观测数据在时间和空间上的变化形态，分别总结了风力发电、光伏发电的建设和运行状态对电磁观测环境干扰的影响特征。结果显示：风力发电机组对电磁观测环境的影响明显小于光伏发电设备，可视为一个干扰点，需规范其与观测设施的距离；若要确保光伏发电设备不影响电磁观测，需要将每个系统视为一个整体，并使二者保持足够大的距离。Abstract: New energy power generation is an important measure taken by the countries to replace non-renewable energy sources so as to achieve sustainable utilization of power resources under the sustainable development strategy. Wind power and photovoltaic power generation are clearly defined as the key development fields of renewable energy. These big power generation facilities will have a great impact on the original seismic observation environment. Based on the electromagnetic observation project of the Xuzhuangzi station of Tianjin, aiming at the wind power and photovoltaic power generation facilities in the observation area around the station in recent years we analyzed the temporal and spatial changes of electromagnetic observation data in different observation environments by the field experimental measurement, and summarized the influence characteristics of the construction and operation of wind power generation and photovoltaic power generation on electromagnetic observation environment interference. The results show that the impact of wind turbine is significantly less than that of photovoltaic power system on the electromagnetic observation environment, which can be regarded as an interference point to regulate the distance between wind turbine and observation facilities. In order to ensure that the photovoltaic power system does not affect the electromagnetic observation, it is necessary to treat each system as a whole and maintain a sufficiently large separation distance between them. This paper can provide a reference for the protection and evaluation of electromagnetic observation environment by building similar new energy power generation equipment in the seismic field, which is of important practical significance for protecting the observation environment of seismic stations.
-
引言
地震波的震相在地震学发展及地球结构研究中起着重要的作用。随着数值计算方法的发展和观测技术水平的提高，逐渐可以利用地震波形中尽可能多的信号来研究地球结构（崇加军，2013）。地震波是源自震源和震中区的弹性波，震相是地震波显示在地震记录图上的信号（赵荣国，1999）。震相分析对深入研究地球的内部构造和探查地震活动规律具有极为重要的意义，地震震相识别也是省级地震分析和编目工作的一项常规工作。地震震相识别中，人工识别的成分占较大比例，而且由于震源机制和传播介质的复杂性，同一种震相在不同地震中的形态均不相同（杨配新等，2004），因此常见的震相Pn，Pg，Sn，Sg较易识别，但非常见的震相识别难度较大。而康拉德面是不连续面，地震波速度在深度变化的过程中均匀变化，因此Pb震相不易被发现，难以得到地震分析人员的关注。
1923年，奥地利研究人员康拉德发现，大陆地壳内部还存在这样一个界面（Conrad，1925）：该界面以上，岩石密度较低，为硅铝层，由中酸性岩石组成；该界面以下，岩石密度较高，为硅镁层（铁镁层），由基性岩石组成。该界面被称为康拉德（Conrad）不连续面，简称康氏面或康拉德面（刘瑞丰等，2014）。此后许多地震学研究人员在大陆地壳内不同区域和不同尺度上也探测到了该界面（Berry，Fuchs，1973；胡德昭等，1989；郭杰等，2013），并对其进行了相关研究。胡德昭等（1989）在中国东南部地壳内探查到康拉德面；1970—1994年期间，苏联的超深钻井项目研究显示康拉德界面上的地震波速度跃变非常大（Richard，1989；Pavlenkova，1993）。但康拉德界面是否普遍存在仍有争议。相关研究认为，康拉德界面与研究区的地壳结构及其演化历史密切相关，在海洋地壳内或者有时在构造活动比较强烈的大陆地壳内，康拉德界面可能是缺失的（Litak，Brown，1989；Richard，1989；Pavlenkova，1993）。近期研究已在有些地质构造活动比较强烈的区域记录到的地震波形中发现有来自康拉德界面的强振幅震相（刘赛君等，2011；郭杰等，2013），且焦煜媛等（2017）已在青藏高原东北缘找到康拉德面存在的直接地震学证据。郭杰等（2013）在位于豫北东的濮阳市地震台的观测中发现了康拉德界面的反射波震相。
关于海南地区Pb震相的研究，也有一些相关的研究成果。范玉兰等（1990）在研究华南地区近震走时表时，描述到海口台、文昌台、定安台、那大台、琼中台等海南省台站在震中距大于151 km时，可记录到清晰的Pb震相；刘赛君等（2011）在海南岛西南海域地壳剖面海陆联合探测研究中提到Pb震相在50—130 km偏移距范围内可被追踪到。
这样看来，有些地区不仅有康拉德面存在的证据，而且前人已记录到海南地区存在Pb震相，因此有必要进一步研究海南地区的Pb震相，故本文一方面以云南省地震台网和广东省地震台网记录到的Pb震相为基础，采用折合走时分析，将其结果作为参考数据；另一方面重点分析海南省地震台网记录到的地震事件，并使用PTD方法（朱元清等，1990）重新定位地震事件的震源深度，以Pb震相理论到时结合实际波形，标注Pb震相，进而拟合标注的震相速度（林建民等，2008）。 此外，综合前人的研究成果（李志雄等，2008；刘赛君等，2011；黄海波等，2012；Kumar et al，2016 ）来分析Pb震相的物理特性，最后确定Pb震相。通过识别Pb震相，探查海南地区康拉德面的深度变化，对该地区一维速度模型的建立、地震定位精度的提高以及地震活动、地震定位深度等研究均具有重要意义。
1. 资料和方法
1.1 资料选取
本文研究区域为海南地区（106°E—114°E，17°N—23°N），选取2009年1月1日至2016年8月25日期间发生在研究区域内的地震，并按照以下标准严格筛选地震事件：① 至少被10个台站记录到；② 地震震级ML≥2.0；③ 地震事件含有Pn震相。最终统计筛选得到Pg震相1 549个，Pn震相588个，无Pb震相。由于研究区域的数据中不包含Pb震相，故以2009年至2016年云南省、广东省地震台网的数据为参考资料，按照上述标准，最终云南地区的震相筛选结果为：Pg震相49 521个，Pb震相3 406个，Pn震相10 781个；广东地区的震相筛选结果为：Pg震相5 038个，Pb震相23个，Pn震相2 240个。
1.2 方法
1.2.1 理论走时
由于近震波的传播路程短，受地球曲率影响小，因此在研究近震问题时，通常将地球表面及各层界面看作水平界面（刘瑞丰等，2014）。Pb震相是在康拉德面上的纵波性首波，其传播路径如图1所示。由该图所示的Pb震相走时路径，依据斯奈尔（Snell）定理和地震波走时方程，可求得Pb的理论走时为
${t_{{\rm{Pb}}}} {\text{＝}} \frac{\varDelta }{{{v_2}}} {\text{＋}} \left( {2{H_1} - h} \right)\frac{{\cos i}}{{{v_1}}},$
(1) 式中，tPb为Pb震相走时，Δ为震中距，H1为康拉德面深度，h为震源深度，v1为一维速度模型中第一层的Pg波平均速度，v2为第二层的Pb波速度，i为Pb波入射角。
1.2.2 折合走时
折合走时与介质厚度（康拉德界面或者莫霍面深度）、震源深度和波速相关（王莉婵等，2016），在同一震源深度的情况下，理论折合走时应为常数值b，震相走时可理解为横向走时和折合走时，因此将式（1）变换为
${t_{\rm redu}} {\text{＝}} t - \frac{\varDelta }{v},$
(2) 式中tredu为折合走时，t为理论震相或观测震相走时，v为波速。
1.2.3 PTD定位方法
PTD定位方法是采用初至波为Pg直达波的台站到时和初至波为Pn首波的台站到时，经转换后的到时差来确定地震的震源深度（朱元清等，1990）。测定震源深度的分辨率为深度每改变5 km，到时差改变0.7 s，换句话说，初动到时测量误差每增加0.1 s，深度误差则增加0.7 km （朱元清等，1997）。
本文可辨别出Pb震相的地震共计52次，由于原震源深度定位偏浅，8 km左右的深度占大部分，因此重新测定震源深度，便于下一步震相理论到时的计算。图2给出了使用PTD定位前后的震源深度比较，可见使用PTD方法测定深度可减小初动到时测量误差引起的深度误差，在深度测定方面有较大的优势。
![]()
图 2 使用PTD方法定位前、后震源深度沿纬度（上）和经度（下）方向的变化比较图Figure 2. Comparison of focal depths of the seismic events along latitude （upper） and longitude （lower） directions before and after location by PTD method1.2.4 模型选取
数据处理涉及理论值计算，速度模型是必不可少的部分，虽然目前大多采用IASP91模型，但由于本文涉及云南省（2015模型）、广东省（华南模型）和海南省（2015模型）的数据（朱元清等，2017），考虑到一维速度模型与到时数据之间的一致性（Wang，2014），本文采用各省的一维速度模型作为各自省份参与计算的理论到时，各省的速度模型列于表1。
表 1 云南、广东和海南各省份的地壳速度模型表Table 1. The crustal velocity model of Yunnan，Guangdong and Hainan Provinces省份 v1/（km·s−1） v2/（km·s−1） v3/（km·s−1） H1/km H/km 云南 6.01 6.60 7.89 20 41 广东 6.00 6.87 7.96 22 33 海南 6.00 6.84 7.97 21 30 注：v1和v2分别表示第一、二层的平均速度，v3表示莫霍面的平均速度，H1表示第一层的平均厚度，H表示莫霍面的平均厚度。 2. Pb震相识别
2.1 云南、广东Pb震相分析
采用折合走时分析云南、广东两省份的数据，结果显示：云南省的Pb震相数据基本处于理论震相范围内，标注的Pb震相接近该地区的康拉德面；广东省的Pb震相数据几乎全在理论震相范围内，标注的Pb震相在康拉德面上。这表明结果良好，与预期结果相符，故可用此方法分析海南的Pb震相数据。
2.2 海南Pb震相识别
震相识别的步骤为：① 对资料进行预处理，筛选出含有Pn震相的地震事件，共计281次；② 对281次地震事件使用PTD法重新定位其震源深度；③ 计算震相理论到时，结合人工分析识别出Pb震相；④ 对识别出的Pb震相进行速度拟合；⑤ 分析和佐证识别的Pb震相结果。
以近几年海南地区发生的影响较大的两次地震为例：事件1为2012年11月5日19时51分海南万宁ML4.1地震，那大台记录的波形见图3左；事件2为2014年7月28日10时10分海南儋州ML3.3地震，万宁台记录的波形见图3右。
那大台的震中距为169 km，PTD重新定位深度为10 km，计算其Pg，Pb，Pn的理论走时分别为28.21，27.26，26.06 s，在实际波形识别的相应结果分别为28.54，27.33，26.33 s。依据波形特性、振幅、周期等因素，Pb震相处于Pn震相与Pg震相之间，其振幅稍大于Pn，但小于Pg振幅，符合绕射波的动力学特征。实际人工识别Pb震相走时t2与理论走时t1相差无几，人工识别的Pb震相标注见图3。
以同样的方法，对万宁台记录到的事件2进行PTD定位计算，结果显示重定位深度为13 km，Pg，Pb，Pn震相的理论走时分别为21.60，21.18，20.71 s，实际波形识别的相应结果分别为21.61，21.22，21.04 s。
依据上述震相识别步骤，在281次地震事件中，从52次地震事件能识别出Pb震相57个，含有Pb震相的地震分布如图4所示。通过分析识别结果可知，能记录到Pb震相的台站共计16个，占全部台站个数（24）的67%，且分布均匀，统计台站记录到Pb震相的震中距处于66.75—234 km范围内，这与海南岛西南海域地壳剖面海陆联合探测研究中追踪到的Pb震相的震中距结果（刘赛君等，2011）基本相近。
采用速度线性拟合（图5）和折合走时（图6）两种方式分析Pb震相的速度，以佐证识别的震相是否与其物理意义相符，图6中折合走时理论值是基于海南模型计算而得。从图5可知，Pb速度为6.47 km/s，在对数据进行折合走时分析的过程中，扰动第二层的速度为6.65 km/s时，能达到图6的效果，这说明实际震相处于理论值范围内。海南地区在区域构造的伸展作用下，其地壳厚度相对正常大陆型地壳较薄，具有西南厚、东北薄的特点（黄海波等，2012），康拉德面作为上下地壳的分界面，穿过界面的波速由6.0 km/s激增至6.4 km/s （刘赛君等，2011）。
![]()
图 4 本文所用含有Pb震相的地震分布图Figure 4. Seismic events distribution with Pb phase used in this paper使用折合走时分析全部Pb震相，结果表明，实际震相基本处于理论震相范围内，尽管有些震相的实际走时与理论走时有偏差，这可能与一维模型单一、不能有效地反映局部地壳内部结构的复杂性和不均匀性有关。综合上述速度拟合分析，识别的Pb震相符合海南地区地震波在康拉德面滑行的物理意义。
3. 海南地区康拉德面
康拉德面是地球内部的次级不连续面，其深度介于10—40 km之间，在陆壳内的平均深度约为20 km （殷伟伟等，2017）。考虑到海南地区的地壳较薄，选取10—30 km作为反演海南地区的康拉德面深度Hi。当Hi以一定步长从10 km依次增加至30 km时，将记录到Pb震相的台站震中距、地震震源深度代入式（1），可得到理论走时t1，t1与实测走时t2之差的绝对值的最小值，即Pb震相的走时残差，为Y=min|t1−t2|，则57次地震的平均走时残差Si为
${S_i} = \frac{1}{{57}}\sum\limits_{i = 1}^{57} {{Y_i}} .$
(3) 为了较好地反映康拉德面的深度和速度情况，通过计算康拉德面的速度和深度变化，可得到平均走时残差Si随深度和速度变化的分布。取Pb波速度为6.40—6.84 km/s，以0.1 km/s的步长计算所有走时残差；深度范围取10—30 km，以1 km的步长计算。计算结果表明速度为6.60—6.80 km/s、深度为19—22 km时残差最小。基于该结果，按上述方法重复，进一步细算速度和深度变化，速度以0.02 km/s为步长，深度以0.5 km为步长，计算结果如图7所示。结果表明，海南地区的Pb波速度为6.60—6.72 km/s、深度为19—21.5 km较为合理。
4. 讨论与结论
本文利用海南省记录到的地震事件，识别了Pb震相并反演得到海南地区的康拉德面深度。结果显示，海南地区的Pb波速度介于6.60—6.72 km/s之间、康拉德面在19—21.5 km左右较合理。海南地区Pb波速度与前人的研究结果（范玉兰等，1990；刘赛君等，2011）基本吻合，首次在海南地区得出康拉德面的合理范围。由于初至Pb震相难识别，缺乏实际震相支持，今后会持续关注这方面的震相数据。
本文结果能对地震分析识别Pb震相起到辅助作用，可为建立海南地区的地壳速度模型提供参考资料。但是，由于本文研究的大部分震源位于上地壳，研究结果具有一定的局限性，震源位于下地壳时震相的识别有待进一步研究。
上海市地震局的朱元清研究员和海南省地震局的李志雄研究员、张慧高工对本研究给予了指导，江苏省地震局廖发军高工分别为本文提供了PTD软件，审稿专家提出了修改意见，作者在此一并表示衷心的感谢！
-
![]()
图 1 天津徐庄子台的基础观测环境
Figure 1. Basic observation environment of Xuzhuangzi station，Tianjin
![]()
图 2 风力发电工作原理示意图
Figure 2. Schematic diagram of working principle of the wind power generation
![]()
图 3 风力发电设备周围磁场变化曲线
（a） 风机静止时不同方向的磁场变化；（b） 风机静止和旋转发电状态下南侧的磁场变化
Figure 3. Magnetic field variation curve around the wind power equipment
（a） The magnetic field variation in different directions when the fan is stationary；（b） The magnetic field variation on the south side when the wind turbine is stationary and rotating，respectively
![]()
图 4 光伏发电工作原理示意图
Figure 4. Schematic diagram of the photovoltaic power generation working principle
![]()
图 5 光伏发电建设施工对徐庄子台地电场观测数据的干扰曲线
Figure 5. Interference curves of photovoltaic power generation construction on the geoelectric field observation data of Xuzhuangzi station
![]()
图 6 徐庄子台2017—2020年地电阻率观测日均值变化曲线
Figure 6. Variation curve of daily mean georesistivity observation at the station Xuzhuangzi from 2017 to 2020
![]()
图 7 光伏发电组件阵列周边不同距离的磁场变化曲线
Figure 7. Magnetic field variation curves at different distances around photovoltaic power generation module array
-
鲍海英,蒋延林,樊晓春,毕雪梅,卜玉菲,张秀霞. 2020. 高压直流输电对地电场观测的影响探讨[J]. 中国地震,36(3):607–619. doi: 10.3969/j.issn.1001-4683.2020.03.022 Bao H Y,Jiang Y L,Fan X C,Bi X M,Bu Y F,Zhang X X. 2020. Study on the impact of high voltage direct current (HVDC) transmission interference on geotechnical stations in Jiangsu Province[J]. Earthquake Research in China,36(3):607–619 (in Chinese).
毕金孟,马永,马朝晖. 2020. 徐庄子地震台地磁观测干扰因素剖析[J]. 地震地磁观测与研究,41(2):130–137. doi: 10.3969/j.issn.1003-3246.2020.02.015 Bi J M,Ma Y,Ma Z H. 2020. An analysis of the disturbance factors of geomagnetic observation at Xuzhuangzi seismic station[J]. Seismological and Geomagnetic Observation and Research,41(2):130–137 (in Chinese).
陈志刚,徐学恭,马永,康健,张明东. 2019. 京沪高铁干扰地电场观测的实验研究[J]. 国际地震动态,(12):36–42. doi: 10.3969/j.issn.0253-4975.2019.12.006 Chen Z G,Xu X G,Ma Y,Kang J,Zhang M D. 2019. Experimental study on interference of Beijing-Shanghai high-speed railway in geoelectric field observation[J]. Recent Developments in World Seismology,(12):36–42 (in Chinese).
冯志生,王建宇,蒋延林,吴晶,梅卫萍. 2001. 地磁垂直分量日变幅逐日比及其与地震关系的探讨[J]. 华南地震,21(2):20–27. doi: 10.3969/j.issn.1001-8662.2001.02.004 Feng Z S,Wang J Y,Jiang Y L,Wu J,Mei W P. 2001. Ratio of daily variation amplitude of geomagnetic vertical component and its relation with earthquake[J]. South China Journal of Seismology,21(2):20–27 (in Chinese).
冯志生,梅卫萍,张秀霞,庄明龙,张晓勇. 2004. 中强震前地磁谐波振幅比的趋势性变化特征初步研究[J]. 西北地震学报,26(1):50–56. Feng Z S,Mei W P,Zhang X X,Zhuang M L,Zhang X Y. 2004. Preliminary study on the characteristics of long and medium term variation for the amplitude ratio of geomagnetic humorous wave before moderate-strong earthquakes[J]. Northwestern Seismological Journal,26(1):50–56 (in Chinese).
蒋延林,张秀霞,杨冬梅,陈俊,赵卫红,张骞,王福才. 2014. 高压直流输电对地磁观测影响的特征分析[J]. 地震,34(3):132–139. doi: 10.3969/j.issn.1000-3274.2014.03.013 Jiang Y L,Zhang X X,Yang D M,Chen J,Zhao W H,Zhang Q,Wang F C. 2014. Influence characteristics of high voltage direct current transmission on geomagnetic observation[J]. Earthquake,34(3):132–139 (in Chinese).
雷晴,颜龙,史勇军. 2020. 基于数据跟踪分析的新疆区域电磁仪器环境干扰特征分析[J]. 地震工程学报,42(6):1507–1516. doi: 10.3969/j.issn.1000-0844.2020.06.1507 Lei Q,Yan L,Shi Y J. 2020. Environmental interference characteristics of electromagnetic instruments in Xinjiang area based on data tracking analysis[J]. China Earthquake Engineering Journal,42(6):1507–1516 (in Chinese).
马钦忠. 2018. 与2008 年汶川MS8.0 地震有关的成都台地电场异常信号[J]. 地震学报,40(3):351–363. Ma Q Z. 2018. Abnormal signals of geoelectric field related to Wenchuan MS8.0 earthquake recorded at the Chengdu station[J]. Acta Seismologica Sinica,40(3):351–363 (in Chinese).
马永,陈志刚. 2017. 地电场观测中常见干扰的对比分析[J]. 大科技,(20):299–300. Ma Y,Chen Z G. 2017. Comparative analysis of common interferences in geoelectric field observation[J]. Super Science,(20):299–300 (in Chinese).
牛微,李珊珊. 2005. 我国新能源发电技术应用现状及发展[J]. 沈阳工程学院学报(自然科学版),1(4):28–30. doi: 10.3969/j.issn.1673-1603.2005.04.008 Niu W,Li S S. 2005. Application situation and prospect of generating electricity technology by new energy in our country[J]. Journal of Shenyang Institute of Engineering (Natural Science)
,1(4):28–30 (in Chinese). 彭玉柱,张世民,李建章. 2015. 中强地震前周至台地电阻率的异常变化[J]. 震灾防御技术,10(2):455–463. doi: 10.11899/zzfy20150227 Peng Y Z,Zhang S M,Li J Z. 2015. Precursory anomalies of geo-resistance before moderate-strong earthquake based on observed data at Zhouzhi seismic station[J]. Technology for Earthquake Disaster Prevention,10(2):455–463 (in Chinese).
王立凤,赵国泽,陈小斌,汤吉,蔡军涛,詹艳,韩冰,程远志,赵凌强,肖骑彬. 2017. 芦山MS7.0地震余震期间大地电磁视电阻率变化[J]. 地震学报,39(1):64–77. Wang L F,Zhao G Z,Chen X B,Tang J,Cai J T,Zhan Y,Han B,Cheng Y Z,Zhao L Q,Xiao Q B. 2017. Variation of magnetotelluric apparent resistivity during the major aftershock period of Lushan earthquake in epicentral area[J]. Acta Seismologica Sinica,39(1):64–77 (in Chinese).
中国地震局. 2019. 地震标准汇编[M]. 北京: 中国标准出版社: 958–975. China Earthquake Administration. 2019. Seismic Standards Collection[M]. Beijing: China Standard Press: 958–975 （in Chinese）.
中国地震局监测预报司. 2002. 地震电磁数字观测技术[M]. 北京: 地震出版社: 120–126. Department of Monitoring and Forecasting, China Earthquake Administration. 2002. Earthquake Observation of Electromagnetic Technology[M]. Beijing: Seismological Press: 120–126 （in Chinese）.
钟李彬,林建,李雪浩,胡俊明,李国超,颜晓晔. 2017. 成灌高铁运行对成都地磁台观测数据干扰测试[J]. 地震地磁观测与研究,38(5):81–85. doi: 10.3969/j.issn.1003-3246.2017.05.015 Zhong L B,Lin J,Li X H,Hu J M,Li G C,Yan X Y. 2017. Testing and analysis of observation data to Chengdu geomagnetic observation on Chengguan high-speed rail[J]. Seismological and Geomagnetic Observation and Research,38(5):81–85 (in Chinese).
下载:
计量
- 文章访问数: 441
- HTML全文浏览量: 394
- PDF下载量: 56
